v1.01.391 PLEXU09 – Validation de VMIS_JOHN_COOK dans CALC_EUROPLEXUS#
Résumé:
Ce test valide l’utilisation de la loi de comportement VMIS_JOHN_COOK dans CALC_EUROPLEXUS. Il valide également la fonctionnalité VARI_INT =”OUI” de ETAT_INIT et l’utilisation en dehors de CALC_EUROPLEXUS de la macro-commande LIRE_EUROPLEXUS.
Modélisation A: reproduction du test EPX bm_ini_med_ecro_vmjc qui valide la prise en compte de variables internes en état initial (issu de Code_Aster) pour la loi VMIS_JOHN_COOK (VMJC de EPX).
Modélisation B: test de non régression pour valider la bonne traduction des paramètres matériaux.
Modélisation A#
But#
Le but de ce test est de valider l’envoie et la transformation (pour la loi VMIS_JOHN_COOK) d’un champ de variables internes en état initial d’un calcul EPX via CALC_EUROPLEXUS. Il valide également la transformation des variables internes dans le sens EPX vers Code_Aster pour la loi VMIS_JOHN_COOK et l’utilisation de LIRE_EUROPLEXUS en dehors de CALC_EUROPLEXUS.
Description#
Ce test est l’équivalent du «bench» EPX bm_ini_med_ecro_vmjc.
Il s’agit d’un cube formé par un seul élément. La face inférieure est encastrée et on applique un chargement (déplacement imposé selon X) sur la face supérieure. Ce calcul est fait avec l’opérateur STAT_NON_LINE afin de produire un état initial pour le calcul EPX.
On lance ensuite CALC_EUROPLEXUS avec cet état initial (déplacements + contraintes + variables internes) et sans chargement supplémentaire que ceux ayant permis d’obtenir l’état initial.
On vérifie alors que les variables internes initiales et en état final sont bien les mêmes que dans bm_ini_med_ecro_vmjc grâce au mot-clé COURBE (et à TEST_TABLE) et que les variables internes récupérées par Code_Aster en fin de calcul sont bien celles attendues.
Les variables internes ne doivent normalement pas évoluer puisque qu’aucun chargement supplémentaire n’a été ajouté. Cependant, certaines composantes qui ne dépendent pas des valeurs calculées sont recalculées totalement par EPX à chaque itération. C’est le cas de ECR5.
Principe de validation#
Comparaison avec le cas test EPX ou références analytiques.
Valeurs testées#
Dans ce premier tableau des numéros des variables internes sont ceux d’EPX.
Maille |
Instant |
Composante |
Point |
Valeur de référence |
Tolérance |
M1 |
\(\mathit{V1}=\mathit{ECR1}\) |
1 |
4.078819E+2 |
1E-4 |
|
M1 |
\(\mathit{V2}=\mathit{ECR2}\) |
2 |
1.420337E+2 |
1E-4 |
|
M1 |
\(\mathit{V3}=\mathit{ECR3}\) |
3 |
5.773488E-3 |
1E-4 |
|
M1 |
\(\mathit{V4}=\mathit{ECR4}\) |
4 |
1.561430E+2 |
1E-4 |
|
M1 |
0.004 |
\(\mathit{V1}=\mathit{ECR1}\) |
1 |
4.078819E+2 |
1E-4 |
M1 |
0.004 |
\(\mathit{V2}=\mathit{ECR2}\) |
2 |
1.420337E+2 |
1E-4 |
M1 |
0.004 |
\(\mathit{V3}=\mathit{ECR3}\) |
3 |
5.773488E-3 |
1E-4 |
M1 |
0.004 |
\(\mathit{V4}=\mathit{ECR4}\) |
4 |
1.561430E+2 |
1E-4 |
M1 |
0.004 |
\(\mathit{V5}=\mathit{ECR5}\) |
5 |
4019.1847623425019 |
1E-4 |
On teste ensuite la valeur des variables internes dans le résultat Aster en sortie de CALC_EUROPLEXUS.
Maille |
Instant |
Composante |
Point |
Valeur de référence |
Tolérance |
M1 |
\(\mathit{V1}\) |
3 |
5.773488E-3 |
1E-4 |
|
M1 |
0.004 |
\(\mathit{V1}\) |
3 |
5.773488E-3 |
1E-4 |
Enfin on teste la valeur des variables internes dans le résultat Aster en sortie de LIRE_EUROPLEXUS.
Maille |
Instant |
Composante |
Point |
Valeur de référence |
Tolérance |
M1 |
\(\mathit{V1}\) |
3 |
5.773488E-3 |
1E-4 |
|
M1 |
0.004 |
\(\mathit{V1}\) |
3 |
5.773488E-3 |
1E-4 |
Modélisation B#
But#
La modélisation A ne permet pas de s’assurer que les paramètres matériaux ont été traduits correctement. Le but de ce test est compléter cela.
Description#
Sur le même modèle que la modélisation A, on applique une pression sur la face supérieure du cube suffisamment forte pour sortir de la partie élastique de la loi. Le calcul est lancé sans état initial.
Principe de validation#
Après avoir vérifié que la traduction des paramètres matériaux était bien faite dans le fichier de commande EPX créé, on applique des tests de non-régression à ce calcul.
Valeurs testées#
Maille |
Composante |
Point |
Valeur de référence |
Tolérance |
M1 |
\(\mathit{SIXX}\) |
1 |
-150717574.84 |
1E-6 |
M1 |
\(\mathit{SIYY}\) |
3 |
-150717574.84 |
1E-6 |
M1 |
\(\mathit{SIZZ}\) |
5 |
-150729740.85 |
1E-6 |
Conclusion#
Les tests ont permis de montrer que:
la loi VMIS_JOHN_COOK était bien prise en compte,
un champ de variables internes pouvait être envoyé en plus des champs de déplacements et de contraintes dans l’état initial,
la transformation des variables internes de Code_Aster vers EPX était conforme à la demande,
la transformation des variables internes d’EPX vers Code_Aster aussi,
l’utilisation de LIRE_EUROPLEXUS en dehors de CALC_EUROPLEXUS est opérationnelle.