u4.84.44 Opérateur POST_LIQUEFACTION#

Syntaxe#

Détail de la syntaxe
evol_noli = POST_LIQUEFACTION(
    ◆ AXE = / "X",
            / "Y",
            / "Z",
    ◆ RESULTAT = evol_noli,
    ◆ CRITERE = / "DP_SIGV_REF",
                / "DP_SIGM_REF",
                / "DP",
                / "P_SIGM",
    # Si: equal_to('CRITERE', 'DP_SIGV_REF')
        ◆ RESU_REF = evol_noli,
        ◆ INST_REF = float,
    # Si: equal_to('CRITERE', 'DP_SIGM_REF')
        ◆ RESU_REF = evol_noli,
        ◆ INST_REF = float,
    # Si: equal_to('CRITERE', 'DP')
        ◆ RESU_REF = evol_noli,
        ◆ INST_REF = float,
)


◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles

Description#

La commande POST_LIQUEFACTION permet de calculer la valeur d’un critère évaluant les zones de liquéfaction dans un ouvrage géotechnique, modélisé par un milieu poreux saturé. Quatre formules sont disponibles pour calculer ce critère:

  • DP_SIGV_REF: le critère est basé sur les valeurs en tout point matériel de la pression interstitielle \({p}_{h}\) à l’instant \(t\) , de la pression interstitielle à l’instant de référence \({t}_{\mathrm{réf}}\) , de la composante du champ de contrainte effective \({\sigma}_{v}^{\mathrm{eff}}\) à l’instant de référence, selon la direction de la gravité, associée au poids des terres. Il est admis que lorsque la valeur du critère dépasse 0,8, on en déduit un risque de liquéfaction locale: consulter par exemple le document DGPR, «Risque sismique et sécurité des ouvrages hydrauliques», version définitive d’octobre 2014. Son expression est la suivante, à l’instant \(t\) et au point \(x\) :

(202)#\[{r}_{u}(t,\mathrm{x})=\frac{{p}_{h}(t,\mathrm{x})–{p}_{h}({t}_{\mathit{réf}},\mathrm{x})}{{\sigma}_{v}^{\mathit{eff}}({t}_{\mathit{réf}},\mathrm{x})}\]
  • DP_SIGM_REF: le critère est calculé à partir des valeurs en tout point matériel de la pression interstitielle \({p}_{h}\) à l’instant \(t\) , de la pression interstitielle à l’instant de référence \({t}_{\mathit{réf}}\) et de la contrainte moyenne effective \({\sigma}_{m}^{\mathit{eff}}\) à l’instant de référence. Son expression est la suivante, à l’instant \(t\) et au point \(x\) :

(202)#\[{r}_{u}(t,\mathrm{x})=\frac{{p}_{h}(t,\mathrm{x})–{p}_{h}({t}_{\mathit{réf}},\mathrm{x})}{{\sigma}_{m}^{\mathit{eff}}({t}_{\mathit{réf}},\mathrm{x})}\]

avec la contrainte moyenne définie par

(202)#\[{\sigma}_{m}^{\mathit{eff}}({t}_{\mathit{réf}},\mathrm{x})=\frac{({\sigma}_{xx}^{\mathit{eff}}({t}_{\mathit{réf}},\mathrm{x})+{\sigma}_{yy}^{\mathit{eff}}({t}_{\mathit{réf}},\mathrm{x})+{\sigma}_{zz}^{\mathit{eff}}({t}_{\mathit{réf}},\mathrm{x}))}{3}\]
  • DP: c’est une différence de pression interstitielle \({p}_{h}\) entre un instant \(t\) et un instant de référence \({t}_{\mathit{réf}}\) . Son expression est la suivante, à l’instant \(t\) et au point \(\mathrm{x}\) :

(202)#\[{r}_{u}(t,\mathrm{x})={p}_{h}(t,\mathrm{x})–{p}_{h}({t}_{\mathit{réf}},\mathrm{x})\]
  • P_SIGM: le critère se calcule à partir de la pression interstitielle \({p}_{h}\) et la contrainte moyenne effective à l’instant \(t\) . Son expression est la suivante, à l’instant \(t\) et au point \(x\) :

(202)#\[ \begin{align}\begin{aligned}{r}_{u}(t,\mathrm{x})=\frac{{p}_{h}(t,\mathrm{x})}{{\sigma}_{m}^{\mathit{eff}}(t,\mathrm{x})}\\ avec la contrainte moyenne définie par\end{aligned}\end{align} \]
(202)#\[{\sigma}_{m}^{\mathit{eff}}(t,\mathrm{x})=\frac{({\sigma}_{xx}^{\mathit{eff}}(t,\mathrm{x})+{\sigma}_{yy}^{\mathit{eff}}(t,\mathrm{x})+{\sigma}_{zz}^{\mathit{eff}}(t,\mathrm{x}))}{3}\]

Ce calcul de post-traitement est indépendant du modèle de comportement mécanique du milieu poreux considéré. Précisons qu’une seule formule peut être associée à la commande. Si l’utilisateur veut calculer plusieurs critères, il doit utiliser plusieurs fois la commande. Le cas-test zzzz403a [V1.01.403] illustre une mise en œuvre de la commande sur un cas fictif et la documentation [U2.04.09] présente les calculs de liquéfaction.

Mots-clefs#

Concept créé par la commande#

Cette commande crée en sortie un concept de type resultat EVOL_NOLI, dans lequel on extraira le champ “SIEF_ELGA” , puis la composante SIPXX, ou SIPYY ou SIPZZ selon la direction de la gravité dans le maillage (voir §4.4).

Mot-clefCRITERE#

Ce mot-clef permet de choisir la formule de calcul du critère.

♦    CRITERE    =    'DP_SIGV_REF'ou'DP_SIGM_REF'ou'DP'ou'P_SIGM'

Mots-clefsRESU_REF et RESULTAT#

♦    RESULTAT    =    resu

Ce mot-clef obligatoire permet de renseigner le nom du concept de typeresultat, pour lequel le calcul de l’ouvrage a été mené, dans lequel on récupère la pression interstitielle à l’instant considéré, sous forme d’un champ “SIEF_ELGA”.

Pour les critères 'DP_SIGV_REF'/'DP_SIGM_REF'/'DP':
♦    RESU_REF    =    resu_ref
Cet opérande obligatoire pour les critères 'DP_SIGV_REF'/'DP_SIGM_REF'/'DP'permet de renseigner le nom du concept de type resultat, pour lequel le calcul de l'ouvrage a été mené, dans lequel on récupère la pression interstitielle à l'instant de référence, sous forme d'un champ 'SIEF_ELGA'.

Remarque:

Seule une structure résultat de type evol_nolipeut être renseignée dans ces mots-clefs.

Mot-clefINST_REF#

Pour les critères 'DP_SIGV_REF'/'DP_SIGM_REF'/'DP'.
♦    INST_REF = inst_ref

Ce mot-clef est obligatoire pour les critères “DP_SIGV_REF”/”DP_SIGM_REF”/”DP”. Il permet de renseigner l’instant de référence où est prise la valeur de la composante du champ de contrainte effective \({\sigma}_{v}^{\mathrm{eff}}\) et de la pression interstitielle \({p}_{h}\) .

Mot-clefAXE#

♦    AXE = /'X'

/”Y” /”Z”

Ce mot-clef obligatoire permet de renseigner la direction :math:`x`, :math:`y`ou :math:`z`du maillage selon laquelle s'exerce la gravité.