Modélisation D
==============


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

Modélisation DKT. On maille seulement la moitié du cylindre (symétrie par rapport au plan :math:`y=0` ) 30 mailles QUAD4 dans la hauteur et 60 sur la demi-circonférence.


.. image:: images/100000000000021A0000038E7D3C7A498A4E3653.png
    :width: 2.6902in
    :height: 4.55in


.. _RefImage_100000000000021A0000038E7D3C7A498A4E3653.png:


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

Nombre de nœuds : 1894

Nombre de mailles et types : 1800 QUAD4

Remarque:

Pour obtenir une solution précise de ce problème, il est nécessaire d'utiliser un maillage assez raffiné (ici 1800 QUAD4).

On observe les erreurs suivantes en fonction de la discrétisation:


.. csv-table::

    "Nombre d'éléments", "Erreur maximale sur le déplacement"
    "450 QUAD4", "0.5%"
    "1800 QUAD4", "0.2%"
    "900 TRIA3", "17%"
    "3600 TRIA3", "1.4%"


On voit que pour ce problème, le maillage en quadrangles est préférable.


Valeurs testées
---------------

**Matériau isotrope**


.. csv-table::

    "**Valeur**", "**Identification**", "**Référence**"
    ":math:`\mathrm{Ur}(z=0)`", ":math:`\mathrm{DY}(\mathrm{PM})`", "5.8018E–05"
    ":math:`\mathrm{Ur}(z=0)`", ":math:`\mathrm{DX}(\mathrm{A1})`", "5.8018E–05"
    ":math:`\mathrm{Ur}(z=0)`", ":math:`\mathrm{DX}(\mathrm{A2})`", "–5.8018E–05"
    ":math:`\mathrm{Uz}(z=L)`", ":math:`\mathrm{DZ}(\mathrm{A3})`", "–2.4429E–05"
    ":math:`\mathrm{Uz}(z=L)`", ":math:`\mathrm{DZ}(\mathrm{A4})`", "–2.4429E–05"
    ":math:`\mathrm{SigmaTT}(z=0)`", ":math:`\mathrm{SIZZ}(\mathrm{PM})`", "2.1375E+06"


**Matériau orthotrope**


.. csv-table::

    "**Valeur**", "**Identification**", "**Référence**"
    ":math:`\mathrm{Ur}(z=0)`", ":math:`\mathrm{DY}(\mathrm{PM})`", "5.8018E–05"
    ":math:`\mathrm{Ur}(z=0)`", ":math:`\mathrm{DX}(\mathrm{A1})`", "5.8018E–05"
    ":math:`\mathrm{Ur}(z=0)`", ":math:`\mathrm{DX}(\mathrm{A2})`", "–5.8018E–05"
    ":math:`\mathrm{Uz}(z=L)`", ":math:`\mathrm{DZ}(\mathrm{A3})`", "–6.10714E–06"
    ":math:`\mathrm{Uz}(z=L)`", ":math:`\mathrm{DZ}(\mathrm{A4})`", "–6.10714E–06"
    ":math:`\mathrm{SigmaTT}(z=0)`", ":math:`\mathrm{SIZZ}(\mathrm{PM})`", "2.1375E+06"