Problème de référence ===================== Géométrie --------- .. figure:: images/Object_1.svg :align: center :width: 480 :height: 360 Hauteur: :math:`h=1` m, largeur : :math:`l=1` m, épaisseur: :math:`e=1` m. .. _RefImage_Object_1.svg: Propriétés du matériau ---------------------- Modèle de comportement de CAM_CLAY ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Les paramètres spécifiques à la loi élastoplastique CAM_CLAY [:ref:`r7.01.14 `], en ce qui concerne les modélisations A, B et C, sont : * :math:`\mu =6` MPa, :math:`PORO=0.66`, :math:`\lambda =0.25`, :math:`\kappa =0.05`, :math:`M=0.9`, :math:`P^0_{cr}=300` kPa, :math:`K_{cam}=0` Pa, :math:`P_{trac}=0` Pa. Le fluide saturant est suposé avoir un module de compressibilité :math:`K_w` infini. La masse spécifique initiale est :math:`\rho_{w}^0=1000\text{kg/m}^3`, et la porosité initiale :math:`\varphi_0=PORO=0.66`. Le couplage poromécanique utilise un coefficient de Biot :math:`b=1`. Modèle de comportement de MCC ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^ Les paramètres spécifiques à la loi élastoplastique MCC [:ref:`r7.01.10 `] sont donnés dans le :numref:`v7.31.122-table_parametres_CSSM` suivant concernant la modélisation D. .. _v7.31.122-table_parametres_MCC: .. list-table:: Paramères tu modèle MCC dans la modélisation D. * - **Effet** - **Appellation** - **Définition** - **Symbole** - **Valeur** * - Élasticité non-linéaire isotrope - | *BulkModulus* | *ShearModulus* | *SwellingIndex* - | Module de compressibilité | Module de cisaillement | Indice de non-linéarité élastique - | :math:`K` | :math:`\mu` | :math:`\kappa` - | 160 MPa | 100 MPa | 50 * - Domaine d'élasticité initial - | *InitCritPress* | *CritStateSlope* | *TensileYieldStress* - | Pression critique initiale | Pente d'état critique | Limite d'élasticité en traction isotrope - | :math:`p_{c0}` | :math:`M` | :math:`\sigma_0` - | 1 MPa | 1 | 10 kPa * - Écrouissage cinématique-isotrope - | *IncoPlastIndex* - | Indice d'incompressibilité plastique - | :math:`\beta` - | 50 Le module de compressibilité du fluide saturant est :math:`K_w=10^3` GPa, sa masse spécifique initiale :math:`\rho_{w}^0=1000\text{ kg/m}^3`, et la porosité initiale :math:`\varphi_0=0.3`. Le couplage poromécanique utilise un coefficient de Biot :math:`b=1`. Modèle de comportement CSSM ^^^^^^^^^^^^^^^^^^^^^^^^^^^ Les paramètres spécifiques au modèle CSSM [:ref:`r7.01.44 `] sont regroupés dans le :numref:`v7.31.122-table_parametres_CSSM` concernant la modélisation E. .. _v7.31.122-table_parametres_CSSM: .. list-table:: Paramètres du modèle CSSM utilisés dans la modélisation E. * - **Intervention** - **Appellation** - **Définition** - **Symbole** - **Valeur** * - Élasticité - | *BulkModulus* | *ShearModulus* | *ShearModulusRatio* - | Module de compressibilité total | Module de cisaillement total | Rapport du module de cisaillement du composant 1 sur le module de cisaillement total - | :math:`K` | :math:`\mu` | :math:`\rho` - | 516 MPa | 238 MPa | 0.1 * - Composant 1 - | *CritStateSlope* | *InitCritPress* | *IncoPlastIndex* | *IsoHardRatio* | *IsoHardIndex* - | Pente d'état critique | Pression critique initiale | Indice d'incompressibilité plastique | Rapport de réduction homothétique du domaine d'élasticité initial | Indice d'écrouissage par agrandissement homothéTique du domaine d'élasticité initial - | :math:`M` | :math:`p_{c0}` | :math:`\beta` | :math:`\eta` | :math:`\omega` - | :math:`1.38` | :math:`100` kPa | :math:`30` | :math:`0.99` | :math:`32` * - Composant 2 - | *HypDistortion* | *HypExponent* | *MinCritPress* - | Distorsion de référence de la relation "hyperbolique modifiée" | Paramètre de courbure de la relation "hyperbolique modifiée" | Pression minimale pour laquelle l'état critique est atteignable - | :math:`\gamma_{\mathrm{hyp}}` | :math:`n_{\mathrm{hyp}}` | :math:`C` - | :math:`2.10^{-4}` | :math:`0.78` | :math:`448` kPa Le module de compressibilité du fluide saturant est :math:`K_w=10^3` GPa, sa masse spécifique initiale :math:`\rho_{w}^0=1000\text{ kg/m}^3`, et la porosité initiale :math:`\varphi_0=0.3`. Le couplage poromécanique utilise un coefficient de Biot :math:`b=1`. **Remarque :** Le coefficient de Biot `b`, le module de compressibilité de la phase solide :math:`K_s`, et le module de compressibilité drainé du squelette :math:`K_0`, celui calculé par le mot-clé ELAS, sont supposés être liés par la relation suivante : .. math:: b=1-\frac{K_0}{K_s} Donc :math:`b=1\rightarrow K_0/K_s=0`, et ce quelque soit la valeur renseignée de :math:`K_0` qui, *a priori*, est un paramètre indépendant de l'élasticité non-linéaire à renseigner dans les mots-clés CAM_CLAY ou MCC. Plus généralement, on pourra se réferer à l'influence des paramètres :math:`b,K_s,K_0` sur l'évolution de la porosité [:ref:`r7.01.10 `] dans le cas plus général où :math:`b<1`. Conditions aux limites et chargements ------------------------------------- Pour les modélisations A, B et C, concernant le modèle CAM_CLAY, le premier trajet de chargement est une compression isotrope, en condition drainéee (:math:`PRE1=0`), jusqu'à la pression :math:`{P}_{sup}`. Pour le second trajet, la pression :math:`{P}_{sup}` est maintenue constante sur les faces latérales et, simultanément, un déplacement imposé vertical de compression est imposé. Ce second chargement est non drainé, ce qui correspond à un flux hydrostatique nul sur toutes les faces. * Pour la modélisation A : :math:`{P}_{sup}={P}_{conso}=2{P}^0_{cr}-P_{trac}=600` kPa (état final contractant) ; * Pour la modélisation B : :math:`{P}_{sup}={P}^0_{cr}=300` kPa (état final critique à variation volumique nulle) ; * Pour la modélisation C : :math:`{P}_{sup}=220` kPa :math:``]). Pour cela, un calcul purement élastique en faisant évoluer la pression de :math:`0` à :math:`1` Pa est effectué. De ce calcul est extrait le champ de contraintes aux points de gauss. Ce champ de contraintes issu du calcul élastique est considéré comme l’état initial de la contrainte hydrostatique nécessaire à la loi CAM_CLAY du calcul suivant (compression isotrope en conditions drainées jusque :math:`{P}_{sup}`).