v6.04.250 SSNV250 – Validation élémentaire de la loi d’endommagement GTN dans le cas a xisymétrique#
Résumé:
Ce cas-test permet de valider l’algorithme d’intégration de la loi d’endommagement de Gurson – Tvergaard – Needleman (\(\mathit{GTN}\) ) local et à gradient avec les éléments finis standard sous-intégrés ou mixtes en grandes déformations. On y modélise un élément volume en traction simple.
Les différentes modélisations traitées sont:
Modélisation A (\(2D\) ): AXIS_SI.
Modélisation B (\(\mathrm{2D}\) ): AXIS_GRAD_VARI.
Modélisation C (\(\mathrm{2D}\) ): AXIS_GRAD_INCO.
Modélisation D (\(3D\) ): 3D_SI.
Modélisation E (\(\mathrm{3D}\) ): 3D_GRAD_VARI.
Modélisation F (\(\mathrm{3D}\) ): 3D_GRAD_INCO.
Solution de référence#
Résultats de référence#
Les solutions de référence sont obtenues en réalisant les mêmes calculs dans le code de calcul par éléments finis Z-set (développé par Mines ParisTech et l’ONERA).
Modélisation A#
Caractéristiques de la modélisation#
Modélisation AXIS_SI.
Caractéristiques du maillage#
Le maillage est obtenu par SALOME.
Nombre de nœuds : 8.
Nombre et types de mailles : 1 QUAD8, 4 SEG3.
Grandeurs testées et résultats#
Le problème considéré ici est homogène, les contraintes (ou les variables internes) sont donc identiques aux différents points de Gauss.
On récupère les valeurs suivantes à l’instant 0.5 et 0.7: la composante \({\sigma}_{yy}\) du tenseur de contrainte (“SIEF_ELGA”, ‘SIYY’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
SIYY |
551.8862 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIYY |
267.7706 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V1 |
1.542931 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V1 |
1.714632 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V2 |
0.09350500 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V2 |
0.1735777 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation B#
Caractéristiques de la modélisation#
Modélisation AXIS_GRAD_VARI.
Caractéristiques du maillage#
Nombre de nœuds : 8.
Nombre et types de mailles : 1 QUAD8, 4 SEG3.
Grandeurs testées et résultats#
Le problème considéré ici est homogène, les contraintes (ou les variables internes) sont donc identiques aux différents points de gauss.
On récupère les valeurs suivantes à l’instant 0.5 et 0.7: la composante \({\sigma}_{yy}\) du tenseur de contrainte (“SIEF_ELGA”, ‘SIYY’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
SIYY |
552.0310 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIYY |
268.0311 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V1 |
1.542972 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V1 |
1.714746 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V2 |
0.0946505 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V2 |
0.1735008 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation C#
Caractéristiques de la modélisation#
Modélisation AXIS_GRAD_INCO.
Caractéristiques du maillage#
Nombre de nœuds : 8.
Nombre et types de mailles : 1 QUAD8, 4 SEG3.
Grandeurs testées et résultats#
Le problème considéré ici est homogène, les contraintes (ou les variables internes) sont donc identiques aux différents points de gauss.
On récupère les valeurs suivantes à l’instant 0.5 et 0.7: la composante \({\sigma}_{yy}\) du tenseur de contrainte (“SIEF_ELGA”, ‘SIYY’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
SIYY |
551.9953 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIYY |
267.9747 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V1 |
1.542963 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V1 |
1.714722 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V2 |
0.09347332 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V2 |
0.1735167 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation D#
Caractéristiques de la modélisation#
Modélisation 3D_SI.
Caractéristiques du maillage#
Nombre de nœuds : 20.
Nombre et types de mailles : 1 HEXA20, 6QUAD8, 12SEG3.
Grandeurs testées et résultats#
Le problème considéré ici est homogène et identique au problème \(A\) (modélisation \(A\) ). Les contraintes (ou les variables internes) sont donc identiques aux différents points de gauss.
On récupère les valeurs suivantes à l’instant 0.5 et 0.7: la composante \({\sigma}_{yy}\) du tenseur de contrainte (“SIEF_ELGA”, ‘SIYY’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
SIYY |
551.8862 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIYY |
267.7706 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V1 |
1.542931 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V1 |
1.714632 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V2 |
0.09350500 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V2 |
0.1735777 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation E#
Caractéristiques de la modélisation#
Modélisation 3D_GRAD_VARI.
Caractéristiques du maillage#
Nombre de nœuds : 20.
Nombre et types de mailles : 1 HEXA20, 6QUAD8, 12SEG3.
Grandeurs testées et résultats#
Le problème considéré ici est homogène et identique au problème \(B\) (modélisation \(B\) ). Les contraintes (ou les variables internes) sont donc identiques aux différents points de gauss.
On récupère les valeurs suivantes à l’instant 0.5 et 0.7: la composante \({\sigma}_{yy}\) du tenseur de contrainte (“SIEF_ELGA”, ‘SIYY’),la composante \({\sigma}_{yy}\) du tenseur de contrainte moyenné par élément (“SIMY_ELGA”, ‘SIYY’), la variable d’écrouissage \(\kappa\) (“VARI_ELGA”, “V1”) et la porosité \(f\) (“VARI_ELGA”, “V2”).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
SIEF_ELGA, SIYY |
552.0310 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIEF_ELGA, SIYY |
268.0311 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
SIMY_ELGA, SIYY |
552.0310 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIMY_ELGA, SIYY |
268.0311 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V1 |
1.542972 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V1 |
1.714746 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V2 |
0.09346505 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V2 |
0.1735008 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation F#
Caractéristiques de la modélisation#
Modélisation 3D_GRAD_INCO.
Caractéristiques du maillage#
Nombre de nœuds : 20.
Nombre et types de mailles : 1 HEXA20, 6QUAD8, 12SEG3.
Grandeurs testées et résultats#
Le problème considéré ici est homogène et identique au problème \(C\) (modélisation \(C\) ). Les contraintes (ou les variables internes) sont donc identiques aux différents points de gauss.
On récupère les valeurs suivantes à l’instant 0.5 et 0.7: la composante \({\sigma}_{yy}\) du tenseur de contrainte (“SIEF_ELGA”, ‘SIYY’),la composante \({\sigma}_{yy}\) du tenseur de contrainte moyenné par élément (“SIMY_ELGA”, ‘SIYY’), la variable d’écrouissage \(\kappa\) (“VARI_ELGA”, “V1”) et la porosité \(f\) (“VARI_ELGA”, “V2”).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
SIEF_ELGA, SIYY |
551.9953 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIEF_ELGA, SIYY |
267.9747 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
SIMY_ELGA, SIYY |
551.9953 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
SIMY_ELGA, SIYY |
267.9747 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V1 |
1.542963 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V1 |
1.714722 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
V2 |
0.09347332 |
SOURCE EXTERNE |
RELATIF – 1% |
0.7 |
V2 |
0.1735167 |
SOURCE EXTERNE |
RELATIF – 1% |
Synthèse des résultats#
Ce cas-test est réalisé sur une seule maille en \(2D\) axisymétrique ou en \(3D\) . Avec les conditions aux limites et les chargements donnés, la réponse devient purement homogène. Les solutions de référence sont obtenues en réalisant les mêmes calculs dans le code de calcul par éléments finis Z-set (développé par Mines ParisTech et l’ONERA). On a unebonne concordance entre les résultats calculés et les solutions de référence.