v6.04.223 SSNV223 - Validation élémentaire de la loi ENDO_SCALAIRE et du pilotage PRED_ELAS pour la modélisation GRAD_VARI#
Résumé:
Ce test a pour but de valider l’algorithme d’intégration de la loi de comportement ENDO_SCALAIRE à gradient de variables internes ainsi que le pilotage PRED_ELAS disponible pour cette loi. Le problème étudié correspond à une sollicitation à déformation homogène imposée pour laquelle on peut obtenir une solution analytique.
Les différentes modélisations traitées sont suivantes:
Modélisation A (\(\mathrm{2D}\) ): On emploie la modélisation D_PLAN_GRAD_VARI.
Modélisation B (\(\mathrm{3D}\) ): On emploie la modélisation 3D_GRAD_VARI.
Solution de référence#
Méthode de calcul utilisée pour la solution de référence#
Ce problème admet une solution analytique. On détermine la relation qui à la déformation imposée \(\varepsilon\) associe le niveau d’endommagement homogène \(a\) (problème réciproque où plus généralement non-local n’admet plus de solution analytique simple). Le problème étant homogène, l’endommagement (en charge) et la déformation sont liés par la relation de cohérence (fonction seuil) :
\({A}^{'}(a)\Gamma (\varepsilon )+{\omega}^{'}(a)=0\)
avec
\(\omega (a)=\mathrm{ka},A(a)=\frac{{(1-a)}^{2}}{{(1-a)}^{2}+\mathrm{ma}(1+\mathrm{pa})}\text{et}\Gamma (\varepsilon )=\left[{c}_{T}\mathrm{tr}\varepsilon +\sqrt{{c}_{H}{\mathrm{tr}}^{2}\varepsilon +{c}_{S}{\varepsilon}_{\mathrm{eq}}^{2}}\right]\)
où les paramètres \({c}_{T},{c}_{H}\text{et}{c}_{S}\) se déduisent des données du problème par :
\({c}_{S}=\frac{E}{2}\left[(1-2\nu ){c}_{\mathrm{comp}}+(1+\nu )\sqrt{\frac{1-2\nu }{2(1+\nu )}{c}_{\mathrm{volu}}+1}\right];{c}_{T}={c}_{\mathrm{comp}}\sqrt{{c}_{S}};{c}_{H}=\frac{1+\nu }{2(1-2\nu )}{c}_{\mathrm{volu}}{c}_{S}\)
On adopte une déformation uniaxiale de la forme:
\(\varepsilon =\varepsilon n\otimes n\text{où}∥n∥=1\text{et}\varepsilon >0\)
Dans ce cas la fonction seuil en déformation s’écrit simplement : \(\Gamma (\varepsilon )=\left[{c}_{T}+\sqrt{{c}_{H}+{c}_{S}}\right]\varepsilon\)
Pour atteindre l’endommagement \(a\) donné, en sollicitant la déformation dans la direction \(n\) , il faut donc imposer une intensité de déformation :
\(\varepsilon =\frac{-k}{{A}^{'}(a)\left[{c}_{T}+\sqrt{{c}_{H}+{c}_{S}}\right]}\)
Pour la solution de référence nous adoptons donc une stratégie suivante : on se fixe le niveau d’endommagement est on vérifie par les calculs EF, que pour une déformation théorique estimée on atteint ce même niveau d’endommagement.
Résultats de référence#
En déformations planes, on adopte une direction de sollicitation \(n=(1/\sqrt{5,}2/\sqrt{5})\) . En \(\mathrm{3D}\) , il vaut \(n=(1/\sqrt{14,}2/\sqrt{14,}3/\sqrt{14})\) . On se fixe comme cible un endommagement \(a=0.6\) ; cela correspond à une intensité de sollicitation \(\varepsilon =9.574237\times {10}^{-4}\) d’après la solution de référence ci-dessus.
Le chargement est appliqué moyennant la technique de pilotage PRED_ELAS dans laquelle on fixe la borne maximale de sorte à atteindre le niveau de déformation \(\varepsilon\) ci-dessus. On vérifiera que l’endommagement correspondant atteint bien \(0.6\) .
Incertitudes sur la solution#
Néant
Références bibliographiques#
Sans objet
Modélisation A#
Caractéristiques de la modélisation#
Une modélisation D_PLAN_GRAD_VARI avec une maille unique, élément QUAD8.
Chargement dans la direction \(n=(1/\sqrt{5,}2/\sqrt{5})\) .
Caractéristiques du maillage#
Nombre de nœuds : 8
Nombre et types de mailles : 1 QUAD8, 4 SEG3
Grandeurs testées et résultats de la modélisation A#
On teste l’endommagement en trois nœuds de la maille, la valeur aux nœuds étant obtenue par extrapolation (CHAM_NO ‘VARI_NOEU’, composante V1).
Identification |
Référence |
Type |
Tolérance |
V1(\(X=2\) , \(Y=0\) ) |
0.6 |
ANALYTIQUE |
RELATIF– 0,1% |
V1(\(X=2\) , \(Y=1\) ) |
0.6 |
ANALYTIQUE |
RELATIF– 0,1% |
V1(\(X=2\) , \(Y=2\) ) |
0.6 |
ANALYTIQUE |
RELATIF– 0,1% |
Modélisation B#
Caractéristiques de la modélisation#
Une modélisation 3D_GRAD_VARI avec une maille unique, élément HEXA20.
Chargement dans la direction \(n=(1/\sqrt{14,}2/\sqrt{14,}3/\sqrt{14})\) .
Caractéristiques du maillage#
Nombre de nœuds : 20
Nombre et types de mailles : 1 HEXA20, 6 QUAD8, 8 SEG3
Grandeurs testées et résultats de la modélisation B#
On teste l’endommagement en trois nœuds de la maille, la valeur aux nœuds étant obtenue par extrapolation (CHAM_NO ‘VARI_NOEU’, composante V1).
Identification |
Référence |
Type |
Tolérance |
V1 (\(X=2\) , \(Y=0\) ) |
0.6 |
ANALYTIQUE |
RELATIF– 0,1% |
V1 (\(X=2\) , \(Y=1\) ) |
0.6 |
ANALYTIQUE |
RELATIF– 0,1% |
V1 (\(X=2\) , \(Y=2\) ) |
0.6 |
ANALYTIQUE |
RELATIF– 0,1% |
Synthèse des résultats#
Ce cas-teste est réalisé sur une seule maille, en conséquence c’est la réponse d’endommagement homogène qui est retrouvé numériquement. La solution de référence est obtenue en se mettant sur le seuil d’endommagement. On note un très bon accord entre la modélisation et la solution de référence. Partie non-locale de la loi néanmoins n’est pas testée.