Modélisation L
==============


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

On utilise une modélisation 3D avec le type d’élément PENTA15 et la méthode LAC de traitement du contact.


Caractéristiques du maillage
----------------------------

On réalise une étude de convergence avec la finesse du maillage de la solution calculée vers la solution analytique. Une suite de maillages obtenus par raffinement uniforme à l'aide de la commande MACR_ADAP_MAIL est utilisée pour le type d’élément PENTA15.

Cependant, en utilisant l’option DECOUPE_HEXA= ‘PYRA’ dans la commande CREA_MAILLAGE, les mailles esclaves (PENTA15) sont découpés en TETRA10 et PYRAM13 pour le traitement du contactpar la méthode LAC :


*   maillage 0 : 8 éléments de type TRIA6, 4 éléments de type QUAD8, 2 éléments de type TETRA10, 6 éléments de type PYRAM13 et 0 élément de type PENTA15


*   maillage 1: 32éléments de type TRIA6, 16éléments de type QUAD8, 8 éléments de type TETRA10, 24éléments de type PYRAM13et 8élément de type PENTA15


*   maillage 2 : 128éléments de type TRIA6, 64éléments de type QUAD8, 32 éléments de type TETRA10, 96 éléments de type PYRAM13et 96 élément de type PENTA15


*   maillage 3 : 512éléments de type TRIA6, 256éléments de type QUAD8, 128 éléments de type TETRA10, 384 éléments de type PYRAM13et 896 élément de type PENTA15



On note que, par rapport aux modélisations E et F, la surface :math:`\text{MAITRE}` courbée est maillée avec 2 TRIA6au lieu d’un unique QUAD8.


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

On teste la vitesse de convergence avec la finesse du maillage de la solution calculée vers la solution analytique en norme :math:`{L}_{2}` :


*   le plus grand réel :math:`{\alpha}_{U}>0` tel que :math:`{\mathrm{\parallel }{\underline{U}}^{\text{calc}}-{\underline{U}}^{\text{ref}}\mathrm{\parallel }}_{0,\Omega }<{C}_{U}\times {h}^{{\alpha}_{U}}` où :math:`{C}_{U}` est indépendant de :math:`h` pour le déplacement;


*   le plus grand réel :math:`{\alpha}_{p}>0` tel que :math:`{\Vert {p}^{\text{calc}}-{p}^{\text{ref}}\Vert }_{0,{\Gamma}_{C}}<{C}_{p}\times {h}^{{\alpha}_{p}}` où :math:`{C}_{p}` est indépendant de :math:`h` pour la pression de contact.


.. csv-table::

    "**Identification**", "**Type de référence**", "**Valeur de référence**"
    ":math:`{\alpha}_{U}`", "'ANALYTIQUE'", "2.5"


.. csv-table::

    ":math:`{\alpha}_{p}`", "'ANALYTIQUE'", "1.0"