Modélisation D : FEM, plaque entière, maillage rayonnant
========================================================

On utilise un maillage rayonnant autour de la pointe du défaut. La partie rayonnante correspond à un disque de rayon :math:`a/10` . Ce disque est constitué de 8 couches d’éléments d’épaisseur constante. L’épaisseur d’une couche est donc :math:`h=a/80` . Des vues du maillage sont données sur les figures ci-dessous. Le maillage utilisé pour le calcul est quadratique, avec des éléments de Barsoum.


.. image:: images/maillage_rayonnant.png
    :width: 3.0437in
    :height: 5.4965in

Figure 6-1: Vue complète du maillage rayonnant de la plaque


.. image:: images/maillage_rayonnant_zoom.png
    :width: 6.3374in
    :height: 2.9701in

Figure 6-2: Zoom sur la fissure

Les conditions aux limites et le chargement sont identiques à ceux définis dans la modélisation A.

Pour le post-traitement en mécanique de la rupture avec la méthode G-thêta, on utilise une couronne d’intégration telle que :math:`\text{R\_INF}=2h` et :math:`\text{R\_SUP}=5h` .


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


Tests sur G
~~~~~~~~~~~


.. csv-table::

    "**Identification**", "**Référence**", "**% tolérance**"
    ":math:`G` de CALC_G, option G", "191.4", "0.6%"
    ":math:`{G}_{\mathit{Irwin}}` de CALC_G, option K", "191.4", "0.6%"
    ":math:`G` de POST_K1_K2_K", "191.4", "1.0%"



Tests sur KI
~~~~~~~~~~~~


.. csv-table::

    "**Identification**", "**Référence**", "**% tolérance**"
    ":math:`{K}_{I}` de CALC_G, option K", "6,646106", "0.3%"
    ":math:`{K}_{I}` de POST_K1_K2_K", "6,646106", "0.5%"