Modélisation C
==============

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

Idem à la modélisation A, avec néanmoins sur la maille *HAUT*, :math:`-\sigma_{zz}` augmente depuis 0 jusqu'à 10 MPa à l'instant :math:`t=1` avant de diminuer jusque 10 kPa à l'instant :math:`t=2` en 100 pas de temps (uniformément distribués).

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

Les solutions sont calculées au point NO8 à l'instant :math:`t=1`. Elles sont données en termes déformation totale :math:`\epsilon_{zz}` et déformation plastique volumique :math:`\epsilon_v^p`. Elles sont comparées à des solutions issues d’un calcul *mtest* au point matériel. Celles-ci sont synthétisées par les deux tableaux ci-dessous.

.. list-table::

 * - **Grandeur**
   - **MTest**
   - **Code Aster**
 * - :math:`\epsilon_{zz}`
   - :math:`-0.054672`
   - :math:`-0.0546715`
 * - :math:`\epsilon_v^p`
   - :math:`-0.030513`
   - :math:`-0.0305132`

La correspondance est satisfaisante. La figure suivante montre l'évolution de la déformation totale :math:`\epsilon_{zz}` avec la contrainte :math:`\sigma_{zz}` appliquée. On y observe l'écrouissage positif, en cas d'écoulement plastique (lorsque :math:`Ip=1`), ainsi que la non-linéarité élastique (lorsque :math:`Ip=0`).

.. figure:: images/figure_modelisation_C.svg
 :align: center
 :width: 960
 :height: 360

 Réponse de la déformation totale :math:`\epsilon_{zz}` avec la contrainte :math:`\sigma_{zz}` appliquée. L'indicateur plastique :math:`Ip` est égal à 1 en cas de chargement plastique, et est égal à 0 autrement.