Anomalie #371
ouvertgrandeur EPAISSEUR_FINALE : interrogation sur sa valeur et son accessibilité en fonction nD
Description
Gérard,
voici le contexte :
je vais (enfin) commencer à voir comment compiler Herezh pour ensuite rendre accessible la grandeur EPAISSEUR_FINALE pour les fonctions nD. J'aurais besoin de cette valeur pour connaitre la déformation III dans le cas de loi 2D_C et l'utiliser pour piloter des lois (Hypoélastique notamment).
J'ai donc commencer à regarder la valeur de EPAISSEUR_FINALE. Dans mon esprit, il y a deux grandeurs différentes dans Herezh :
- EPAISSEUR_MOY_FINALE : l'épaisseur moyenne sur l'élément
- EPAISSEUR_FINALE : l'épaisseur à un point d'intégration
Déjà, deux questions importantes pour moi :
1) est-ce que c'est bien cela qui est sensé être implanté dans Herezh ?
2) et si non, quelle serait la(les) grandeur(s) à exploiter dans une fonction nD pour connaitre la déformation III dans une loi 2D_C ?
La réponse à ces questions va me permettre de commencer dès maintenant à travailler sur le code Herezh à titre de premier exercice => rendre accessible des nouvelles grandeurs dans les fonctions nD.
Et ensuite, de quoi débugger si besoin :
Par le test, je me rends compte qu'en fait ces deux grandeurs EPAISSEUR_FINALE et EPAISSEUR_MOY_FINALE donnent exactement la même chose. Quelque soit le cas (chargement, loi de comportement), j'obtiens toujours la même valeur EPAISSEUR_FINALE sur les 4 points d'intégration. Même en imposant des déplacements conduisant à un champ de déformation hétérogène => on a une égalité stricte au digit près pour tous les points d'intégration.
Dans l'exemple joint, il s'agit d'un quadrangle 1x1 (épaisseur 0.038) avec des conditions de blocage de type traction uniaxiale (bord X=0 bloqué selon UX, bord Y=0 bloqué selon Y) et j'impose un déplacement X différent sur les 2 noeuds en X=1. En conséquence, on obtient un champ de déformation hétérogène, qui normalement devrait donner une épaisseur différente aux divers points d'intégration.
Dans le .info, j'ai proposé plusieurs lois. Mais je me concentre sur la loi HYPOELAS 3D (MAT_hypo), car c'est la seule pour laquelle je suis capable de calculer l'épaisseur actuelle en fonction soit de la déformation III (DEF_EPAISSEUR), soit de la déformation principale I et II connaissant le coef de Poisson 0.4.
Le fichier gnuplot gnu affiche pour le premier point d'intégration :
- EPAISSEUR_FINALE (points rouge)
- EPAISSEUR_MOY_FINALE (points verts)
- épaisseur finale calculée en fonction de def principale I et II + coef Poisson 0.4 (trait bleu)
- épaisseur finale calculée en fonction de DEF_EPAISSEUR (trait rose)
Voilà ce qui est obtenu pour le point d'intégration 1 :
Et comme autre résultat, on peut également constater que EPAISSEUR_FINALE est la même partout dans l'élément malgré un état de déformation variable. Par exemple au dernier incrément, on a dans l'ordre (EPAISSEUR_FINALE, def principale I, def_principale II) :
point 1 : (3.690643220818e-02, 7.881446996732e-02, -3.545869816873e-02)
point 2 : (3.690643220818e-02, 8.101757499525e-02, -2.690764388003e-02)
point 3 : (3.690643220818e-02, 5.556599910222e-02, -3.533466140407e-02)
point 4 : (3.690643220818e-02, 5.734173928020e-02, -2.582961933585e-02)
Et pour finir, si on calcule la vraie épaisseur aux 4 points d'intégration et qu'on en fait la moyenne, on retombe bien sur EPAISSEUR_MOY_FINALE et EPAISSEUR_FINALE :
merci
Julien
Fichiers