Modélisation H
==============

Caractéristiques de la modélisation
-----------------------------------

Simulation au point matériel. L'objectif est de valider la compatibilité de CALC_ESSAI_GEOMECA avec le modèle de comportement MCC. 
On réalise, pour ce faire, un essai triaxial drainé monotone à déplacement imposé ESSAI_TRIA_DR_M_D. 
Le nombre de pas de temps est de :math:`200`; la pression de confinement est de :math:`1` MPa; 
la déformation axiale maximale est égale à 5% en compression. Le critère de convergence est :math:`10^{-10}`.

Les paramètres du modèle MCC sont donnés dans le :numref:`v6.07.112-table_parametres_MCC`.

.. _v6.07.112-table_parametres_MCC:

.. list-table:: Paramètres du modèle MCC utilisés dans la modélisation H.
  
 * - **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  


Grandeurs testées et résultats
------------------------------

Deux tests de non-regression sont réalisés en fin d’essai sur les composantes SIG_LAT (contrainte latérale, SIXX ou SIYY) et SIG_AXI (contrainte axiale, SIZZ). 
À noter que ces contraintes sont positives (convention suivie ici de la mécanique des sols). Les réponses correspondent à une évolution avec écrouissage positif et contractance.

.. list-table::
  
 * - **Grandeur**
   - **Code Aster**
 * - SIG_LAT
   - :math:`1000000` Pa
 * - SIG_AXI 
   - :math:`2392637.9437635597` Pa
   
  
.. list-table:: 
   :width: 100%
   :class: borderless


   * - .. figure:: images/modelisationH_EPS_AXI_QDEV.png
          :align: center
          :width: 960
          :height: 360

          Réponses de la contrainte équivalente QDEV (Pa) avec la déformation axiale EPS_AXI.
         
     - .. figure:: images/modelisationH_EPS_AXI_EPS_VOL.png
          :align: center
          :width: 960
          :height: 360

          Réponses de la déformation volumique EPS_VOL avec la déformation axiale EPS_AXI.