v6.04.255 SSNV25 5 – Validation de la loi d’endommagement GTN à gradient avec des simulations d’unbarreau en déformations planes#
Résumé:
Ce cas-test permet de valider l’algorithme d’intégration de la loi d’endommagement de Gurson – Tvergaard – Needleman (\(\mathit{GTN}\) ) à gradient avec les éléments finis standard sous-intégrés ou mixtes en grandes déformations. On y modélise un barreau contenant 10 éléments finis.
Les différentes modélisations traitées sont:
Modélisation A (\(3D\) ): 3D_GRAD_VARI
Modélisation B (\(3D\) ): 3D_GRAD_INCO
Modélisation C (\(2D\) ): D_PLAN_GRAD_VARI
Modélisation D (\(2D\) ): D_PLAN_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 3D_GRAD_VARI.
Caractéristiques du maillage#
Le maillage est obtenu par SALOME.
Nombre de nœuds : 128.
Nombre et types de mailles : 10 HEXA20, 51QUAD8, 84 SEG3.
Grandeurs testées et résultats#
On récupère des valeurs suivantes à l’instant 0.5 et 0.8: la contrainte \({\sigma}_{xx}\) (“SIEF_ELGA”, ‘SIXX’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”) de deux points d’intégration. Le premier point d’intégration est celui le plus proche du nœud \(A\) (on le note \(G1\) ), le deuxième est celui le plus éloigné du nœud \(A\) (on le note \(G2\) ).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Point d’intégration |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
G1 |
SIXX |
731.0701 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
SIXX |
889.0125 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
SIXX |
135.7807 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
SIXX |
860.9395 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V1 |
1.407778 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V1 |
1.437888 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V1 |
1.643740 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V1 |
1.843401 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V2 |
0.06601665 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V2 |
0.007703500 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V2 |
0.2150436 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V2 |
0.02104307 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation B#
Caractéristiques de la modélisation#
Modélisation 3D_GRAD_INCO.
Caractéristiques du maillage#
Nombre de nœuds : 128.
Nombre et types de mailles : 10 HEXA20, 51QUAD8, 84 SEG3.
Grandeurs testées et résultats#
On récupère des valeurs suivantes à l’instant 0.5 et 0.8: la contrainte \({\sigma}_{xx}\) (“SIEF_ELGA”, ‘SIXX’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”) de deux points d’intégration. Le premier point d’intégration est celui le plus proche du nœud \(A\) (on le note \(G1\) ), le deuxième est celui le plus éloigné du nœud \(A\) (on le note \(G2\) ).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Point d’intégration |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
G1 |
SIXX |
726.2872 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
SIXX |
891.3942 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
SIXX |
128.4041 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
SIXX |
867.2785 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V1 |
1.407362 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V1 |
1.438084 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V1 |
1.641179 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V1 |
1.843949 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V2 |
0.06660734 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V2 |
0.007312952 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V2 |
0.2166102 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V2 |
0.01996904 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation C#
Caractéristiques de la modélisation#
Modélisation D_PLAN_GRAD_VARI.
Caractéristiques du maillage#
Nombre de nœuds : 53.
Nombre et types de mailles : 10QUAD8, 31SEG3.
Grandeurs testées et résultats#
On récupère des valeurs suivantes à l’instant 0.5 et 0.8: la contrainte \({\sigma}_{xx}\) (“SIEF_ELGA”, ‘SIXX’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”) de deux points d’intégration. Le premier point d’intégration est celui le plus proche du nœud \(A\) (on le note \(G1\) ), le deuxième est celui le plus éloigné du nœud \(A\) (on le note \(G2\) ).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Point d’intégration |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
G1 |
SIXX |
731.0701 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
SIXX |
889.0125 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
SIXX |
135.7807 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
SIXX |
860.9395 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V1 |
1.407778 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V1 |
1.437888 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V1 |
1.643740 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V1 |
1.843401 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V2 |
0.06601665 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V2 |
0.007703500 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V2 |
0.2150436 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V2 |
0.02104307 |
SOURCE EXTERNE |
RELATIF – 1% |
Modélisation D#
Caractéristiques de la modélisation#
Modélisation D_PLAN_GRAD_INCO.
Caractéristiques du maillage#
Nombre de nœuds : 53.
Nombre et types de mailles : 10QUAD8, 31SEG3.
Grandeurs testées et résultats#
On récupère des valeurs suivantes à l’instant 0.5 et 0.8: la contrainte \({\sigma}_{xx}\) (“SIEF_ELGA”, ‘SIXX’), la variable d’écrouissage \(\kappa ` ('VARI_ELGA', 'V1') et la porosité :math:`f\) (“VARI_ELGA”, “V2”) de deux points d’intégration. Le premier point d’intégration est celui le plus proche du nœud \(A\) (on le note \(G1\) ), le deuxième est celui le plus éloigné du nœud \(A\) (on le note \(G2\) ).
Le tableau suivant rassemble les valeurs obtenues par le logiciel de calcul par éléments finis Z-set.
INST |
Point d’intégration |
Identification |
Valeur de r éférence |
Type |
Tolérance |
0.5 |
G1 |
SIXX |
726.2872 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
SIXX |
891.3942 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
SIXX |
128.4041 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
SIXX |
867.2785 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V1 |
1.407362 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V1 |
1.438084 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V1 |
1.641179 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V1 |
1.843949 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G1 |
V2 |
0.06660734 |
SOURCE EXTERNE |
RELATIF – 1% |
0.5 |
G2 |
V2 |
0.007312952 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G1 |
V2 |
0.2166102 |
SOURCE EXTERNE |
RELATIF – 1% |
0.8 |
G2 |
V2 |
0.01996904 |
SOURCE EXTERNE |
RELATIF – 1% |
Synthèse des résultats#
Ce cas-test est réalisé sur un barreauen \(2D\) déformation plane ou en \(3D\) . 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.