r4.20.01 Indicateurs de décharge et de perte de proportionnalité du chargement en élastoplasticité#
Résumé
On présente un ensemble de paramètres scalaires appelés indicateurs, permettant d’apprécier une décharge ou une perte de proportionnalité d’un chargement au cours de son histoire. Deux types d’indicateurs sont proposés: des indicateurs se présentant sous la forme de champs scalaires permettant de détecter les zones de la structure subissant des décharges ou des chargements non radiaux, et des indicateurs globaux intégrés sur une zone de la structure choisie par l’utilisateur. Ces derniers sont plus spécialement destinés à l’évaluation de la validité du taux de restitution de l’énergie en mécanique de la rupture non linéaire.
Les indicateurs décrits dans ce document sont disponibles:
pour les indicateurs locaux sous la commande CALC_CHAMP, options DERA_ELGA et DERA_ELNO;
pour les indicateurs globaux, sous la commande POST_ELEM, options INDIC_SEUIL et INDIC_ENER.
Indicateurs locaux#
Le but de ces indicateurs est de cerner les zones de la structure où, à un instant particulier, se produit soit une décharge ou une perte de radialité du champ de contraintes. Ils sont produits en post‑traitement d’un calcul statique ou dynamique, 2D ou 3D, utilisant une loi de comportement élastique ou non. Ils se présentent sous la forme de champs de scalaires dont le dépouillement peut être effectué en traçant leurs isovaleurs par un post-processeur graphique.
Indicateurs de décharge#
Indicateur local de décharge totale#
Cet indicateur mesure au point \(M\) et entre l’instant \(t\) et \(t+\Delta t\) , la variation relative de la norme des contraintes au sens de Von Misès. Il s’écrit formellement: \({I}_{1}=\frac{\parallel \sigma (M,t+\Delta t)\parallel -\parallel \sigma (M,t)\parallel }{\parallel \sigma (M,t+\Delta t)\parallel }\) . Cette quantité est négative en cas de décharge locale au point M. La norme \(\parallel \sigma (M,t)\parallel\) peut s’écrire de quatre façons différentessuivant le choix du modélisateur :
\(\parallel \sigma (M,t)\parallel =\sqrt{\frac{3}{2}{\sigma}^{D}.{\sigma}^{D}}\) , où \({\sigma}^{D}\) est la partie déviatorique du tenseur des contraintes (cette norme est utile en plasticité avec écrouissage isotrope).
\(\parallel \sigma (M,t)\parallel =\sqrt{\frac{3}{2}\sigma .\sigma }\) , où l’on considère la totalité du tenseur des contraintes afin de détecter par exemple les diminutions de pression hydrostatique.
\(\parallel \sigma (M,t)\parallel =\sqrt{\frac{3}{2}({\sigma}^{D}-X).({\sigma}^{D}-X)}\) , avec \(X\) le tenseur de la contrainte de rappel dans le cas d’une loi élastoplastique avec un écrouissage cinématique.
\(\parallel \sigma (M,t)\parallel =\sqrt{\frac{3}{2}(\sigma -X).(\sigma -X)}\)
Indicateur local de décharge élastoplastique#
Cet indicateur permet de savoir si la décharge reste élastique ou s’il y aurait un risque de plastification si on utilisait un écrouissage cinématique pur. C’est un indicateur «extrême», sachant que, pour les métaux, les écrouissages isotropes et cinématiques sont tous deux présents.
L’option DERA_ELGA calcule donc (en plus des composantes DCHA_V et DCHA_T) les composantes :
I_decha=IND_DCHA :
IND_DCHA=0 valeur initiale sans contrainte.
IND_DCHA=1 si charge élastique
IND_DCHA=2 si charge plastique
IND_DCHA=-1 si décharge élastique licite (quelque soit le type d’écrouissage)
IND_DCHA=-2 si décharge abusive (on aurait plastifié avec un écrouissage cinématique).
VAL_DCHA : indique le proportion de sortie du critère (voir plus loin).
Le fonctionnement est le suivant : pour les lois VMIS_ISOT * seulement en chaque point d’intégration, à chaque instant \(t\) , à partir du tenseur des contraintes \(\sigma (t)\) , de la déformation plastique équivalente cumulée \(p(t)\) , et de la courbe d’écrouissage isotrope \(R(p(t))\) ,
Initialisation : IND_DCHA=0, VAL_DCHA=0.
tant que \(p(t)=0\) , IND_DCHA=1 (élasticité),
si IND_DCHA=-2 (critère de décharge abusive atteint), on ne fait plus rien
si \(\Delta p(t)>0\) :
si on est en charge, donc si l’angle entre le tenseur accroissement de contraintes et le tenseur contraintes totale est «petit» : , \(\frac{\tilde{\sigma}(t-\Delta t)\cdot \Delta \tilde{\sigma}}{∣∣\tilde{\sigma}(t-\Delta t)∣∣∣∣\Delta \tilde{\sigma}∣∣}>0\) IND_DCHA=2 ; (points entre A et B sur la figure ).
si on est en décharge, IND_DCHA=-2 (cas rare : cela signifie qu’en un pas de temps, on franchit la surface de charge dans une direction éloignée de celle du pas précédent)
si \(\Delta p(t)=0\) :
si IND_DCHA>-1, on calcule le tenseur \(X\) (centre du domaine d’élasticité si on était en cinématique pur, donc si la surface de charge initiale, représentée par un cercle dans le plan déviateur, s’était translatée jusqu’au point courant) par :
\(X=\sigma (t)\frac{R(p)-R(0)}{R(p)}\)
si IND_DCHA=-1, on utilise le tenseur \(X\) déjà calculé.
on calcule le critère « cinématique » \({(\sigma (t)-X)}_{\mathrm{eq}}\le R(0)\)
Si ce critère est vérifié, alors la décharge serait élastique aussi avec un écrouissage cinématique, elle est donc «licite» : IND_DCHA=-1 (tout point entre B et E sur la figure 2). Pour appliquer le critère au prochain instant, on stocke \(X\) (6 composantes)
sinon, IND_DCHA=-2 et garde cette valeur (car la suite du calcul serait modifiée si l’écrouissage était cinématique). (points entre E et F, ou entre C et C” sur la figure suivante). \(\text{VAL\_DCHA}=\frac{∥\sigma -X∥}{{R}_{0}}\)
Figure 2.1.2-a: indicateur de décharge élastoplastique
Commentaires :
L’inconvénient de cette méthode vient de ce qu’elle n’est pas systématique : il faut que l’utilisateur «pense» à appeler l’option DERA_ELGA.
Le critère obtenu est relativement sévère : il suppose que l’on remplace l’écrouissage isotrope pur par un écrouissage cinématique pur. Une possibilité d’évolution consisterait à définir un rayon de surface seuil plus grand que \({R}_{0}\) .
Indicateurs de perte de radialité#
Indicateur local de radialité du chargement#
Cet indicateur mesure au point \(M\) et entre l’instant \(t\) et \(t+\Delta t\) , la variation de la direction des contraintes. Il s’écrit:
\({I}_{2}=1-\frac{\mid \sigma (M,t).\Delta \sigma \mid }{\parallel \sigma (M,t)\parallel \parallel \Delta \sigma \parallel }\) , où le produit scalaire «.» est associé à l’une des quatre normes précédentes. Cette quantité est nulle lorsque la radialité est conservée durant l’incrément de temps. Ce critère peut aussi être interprété comme la quantité \(1-\cos(\theta )\) , où \(\theta\) est l’angle entre \(\sigma\) et \(\Delta \sigma\) . Cet indicateur à valeur réelle évolue entre 0 pour les chargements radiaux et 1. Il est utile en particulier pour déterminer la validité d’une solution élasto-plastique en mécanique de la rupture,
Indicateur d’erreur due à la discrétisation temporelle#
Il fournit une mesure de l’erreur \(\eta\) due à la discrétisation en temps, directement reliée à la rotation de la normale à la surface de charge. On calcule l’angle entre \({n}^{-}\) , la normale au critère de plasticité au début du pas de temps (instant t-), et \({n}^{+}\) , la normale au critère de plasticité calculée à la fin du pas de temps de la façon suivante : \({I}_{\eta}=\frac{1}{2}∣∣\Delta n∣∣=\frac{1}{2}∣∣{n}^{+}-{n}^{-}∣∣=∣\sin(\frac{\alpha}{2})∣\) . Cet indicateur est directement liée à la norme de la variation de la normale au convexe deplasticité (ceci peut se généraliser facilement à toute loi élastoplastique à écoulement normal), et elle peut également s’interpréter comme le sinus de la moitié de l’angle les deux normales.Cela fournit une mesure de l’erreur (voir [4]). Ce critère est opérationnel pour les comportements élastoplastiques de Von Mises à écrouissage isotrope, cinématique linéaire et mixte. Il peut être utilisé pour piloter la subdivision automatique du pas de temps.
Indicateurs globaux#
Ces indicateurs sont destinés à détecter si, au cours de l’histoire de la structure et jusqu’à l’instant actuel \(t\) , et pour une zone de la structure choisie par le modélisateur, il y a eu perte de proportionnalité du chargement (ces indicateurs laissent donc une trace de l’histoire contrairement aux indicateurs locaux qui sont instantanés). Ils sont uniquement utilisables dans le cadre d’un comportement élastoplastique avec écrouissage isotrope (en 2D ou 3D).
Indicateur sur les paramètres de plasticité#
Cette quantité permet, dans le cas de la plasticité de Von Mises avec écrouissage isotrope, d’une part de savoir (en moyenne sur une zone \({\Omega}_{S}\) du domaine \(\Omega\) ) si les contraintes et les déformations plastiques ont les mêmes directions et si le seuil plastique est atteint à l’instant actuel, et d’autre part si au cours de l’histoire la déformation plastique a changé de direction. Cette quantité s’écrit:
\({I}_{3}=\frac{1}{{\Omega}_{S}}{\int}_{{\Omega}_{S}}(1-\frac{\mid \sigma .{\varepsilon}^{p}\mid }{({\sigma}_{Y}+R(p))p})d\Omega\) ,
où \({\sigma}_{Y}\) est le seuil plastique initial, \(R\) l’extension de la surface de charge liée à l’écrouissage et \(p\) la déformation plastique cumulée. Le produit scalaire est associé à la norme au sens de Von Misès. Cet indicateur est normalisé et a une valeur comprise entre zéro et un. Il est nul si le chargement a conservé son caractère proportionnel en chaque point de \({\Omega}_{S}\) tout au long de l’histoire écoulée.
Remarque 1:
L’indicateur n’est pas affecté si au cours de l’histoire, il y a eu des décharges puis des recharges élastiques sans changement de direction des contraintes lorsque l’on revient sur le seuil (cf.[Figure 3.1-a]).
Figure 3.1-a : Chemin de chargement entre \({t}_{1}\) et \({t}_{2}\) dans le plan déviatorique des contraintes
Remarque 2:
Dans la formulation de cet indicateur interviennent trois ingrédientsliés à la plasticité :
l’écart entre la direction des contraintes et des déformations plastiques actuelles ( \(\sigma .\varepsilon p\ne \parallel \sigma \parallel \parallel {\varepsilon}^{p}\parallel\) ),
la position des contraintes par rapport au seuil actuel ( \(\parallel \sigma \parallel \le ({\sigma}_{Y}+R(p))\) ),
l’écart entre la norme de la déformation plastique actuelle et la déformation plastique cumulée ( \(\parallel {\varepsilon}^{p}\parallel \ne p\) ).
Une perte de proportionnalité a pu se produire au cours de l’histoire sans que l’indicateur ne le détecte par l’intermédiaire des deux premiers ingrédients (c’est-à-dire que l’on peut avoir en fin de chargement coïncidence des directions des contraintes et des déformations plastiques et être sur le seuil plastique). En revanche, on aura , \(\parallel {\varepsilon}^{p}\parallel \ne p\) et l’indicateur sera obligatoirement supérieur à zéro, par conséquent l’utilisateur sera averti de la perte de proportionnalité.
Remarque 3:
Si l’indicateur détecte obligatoirement une perte de proportionnalité dans une zone, dans la pratique il faut que cette dernière contienne suffisamment de points matériels à chargement non proportionnel. En effet, si l’on choisit une zone très vaste avec peu de points concernés, la normalisation de l’indicateur effectuée avec la division par le volume de la zone implique un certain «écrasement» vers zéro de la valeur de la quantité. Typiquement, pour une structure contenant un défaut source de non proportionnalité, on a intérêt à choisir une zone d’intégration \({\Omega}_{S}\) entourant le défaut avec un faible voisinage afin d’obtenir une valeur significative de l’indicateur.
Indicateur énergétique#
Cet indicateur a la même fonction que le précédent, mais est fondé sur la densité d’énergie. Il s’écrit:
\({I}_{4}=\frac{1}{{\Omega}_{S}}{\int}_{{\Omega}_{S}}(1-\frac{\psi}{\omega})d\Omega\) ,
où \(\omega\) est la densité d’énergie de déformation définie par: \(\omega (t)={\int}_{0}^{t}\sigma .\dot{\varepsilon}d\tau\) , et \(\psi\) est la densité d’énergie élastique associée à la courbe de traction si on considérait le matériau élastique non linéaire. Plus précisément, cette quantité s’écrit:
\(\begin{array}{}\psi (\varepsilon (t))=\frac{1}{2}{\text{Ktr}}^{2}(\varepsilon )+\frac{2\mu }{3}{\parallel \varepsilon \parallel }^{2},\text{si}\parallel \sigma \parallel <({\sigma}_{Y}+R(p)),\\ \psi (\varepsilon (t))=\frac{1}{2}{\text{Ktr}}^{2}(\varepsilon )+\frac{{R}^{2}(P)}{6\mu }+{\int}_{0}^{P}R(q)\text{dq},\text{si}\parallel \sigma \parallel =({\sigma}_{Y}+R(p)),\end{array}\)
avec \(K\) le module de compressibilité, \(\mu\) le coefficient de cisaillement de Lamé, \(R\) le seuil de la courbe de traction associé à la norme de la déformation plastique \(P=\parallel {\varepsilon}^{p}\parallel\) (celui-ci peut être différent du vrai seuil plastique, car \(P\ne p\) si le chargement est non proportionnel). Cet indicateur est également normalisé entre 0 et 1. Il est nul pour un chargement ayant toujours gardé son caractère proportionnel \((\psi \ne \omega )\) .
Fonctionnalités et vérification#
Les indicateurs présentés ici sont utilisables en post-traitement d’un calcul mécanique et sont disponibles pour les éléments finis des milieux continus en 2D (mode de déformations planes, contraintes planes ou axisymétrique, mailles triangulaires ou quadrangulaires) ou 3D (mailles hexaédriques, tétraédriques, pentaédriques ou pyramides). Les éléments filaires, poutres, plaques et coques sont exclus de cette application.
Indicateurs locaux#
Ces indicateurs sont accessibles après un calcul statique ou dynamique quelle que soit la loi de comportement du matériau. L’opérateur CALC_CHAMP présente les options “DERA_ELGA” et “DERA_ELNO” pour l’indicateur de décharge \({I}_{1}\) et l’indicateur de perte de radialité \({I}_{2}\) évalués aux nœuds ou aux points de Gauss de l’élément [2] . Dans le cas général, les deux premières normes décrites au paragraphe 2.1 sont calculées, les deux dernières n’étant utilisées que si l’on a réalisé au préalable un calcul élastoplastique avec écrouissage cinématique.
Les composantes du champ ainsi calculé sont:
DCHA_T: l’indicateur de décharge \({I}_{1}\) calculé à partir du tenseur des contraintes;
DCHA_V: l’indicateur de décharge \({I}_{1}\) calculé à partir du déviateur des contraintes;
IND_DCHA: l’indicateur de décharge \({I}_{\mathrm{decha}}\)
IND_DCHA=0 valeur initiale sans contrainte
IND_DCHA=1 si charge élastique
IND_DCHA=2 si charge plastique
IND_DCHA=-1 si décharge élastique licite (quelque soit le type d’écrouissage)
IND_DCHA=-2 si décharge abusive (on aurait plastifié avec un écrouissage cinématique).
VAL_DCHA: l’indicateur de la proportion de sortie du critère
RADI_V: l’indicateur de perte de radialité \({I}_{2}\) calculé à partir du déviateur des contraintes.
ERR_RADI: l’indicateur d’erreur \({I}_{\eta}\) due à la discrétisation temporelle
Ces indicateurs sont vérifiés par les cas tests suivants:
Indicateurs globaux#
Ces indicateurs ne sont accessibles qu’après un calcul élastoplastique avec écrouissage isotrope. L’opérateur POST_ELEM présente les options ‘INDIC_SEUIL’ et ‘INDIC_ENER’ correspondant respectivement aux indicateurs globaux \({I}_{3}\) et \({I}_{4}\) . Ceux-ci sont évalués sur un groupe de maille précédemment défini par l’utilisateur (par exemple par la commande DEFI_GROUP).
Ces indicateurs sont vérifiés par les cas tests suivants:
Bibliographie#
LEMAITRE, J.-L.CHABOCHE, Mécanique des matériaux solides, Dunod 1985.
Manuel d’Utilisation de Code_Aster . Document [U4.61.02].
CR-AMA-11.035 : «Projet LOCO - lot B2 - Indicateurs locaux de décharge et de radialité» A.Foucault , 24/01/2011
CR-AMA-11.30 «Etude d’un critère de radialité» J.M.Proix
Description des versions du document#
Version Aster |
Auteur(s) Organisme(s) |
Description des modifications |
4 |
G. DEBRUYNE (EDF/IMA/MMN) |
Texte initial |
\(11.1\) |
J.M.PROIX EDF/R&D/AMA |
Ajout des nouveaux indicateurs de décharge. |
\(11.2\) |
J.M.PROIX EDF/R&D/AMA |
Ajout d’un nouveau critère de radialité. |