Code_Aster ®
Version
4.0
Titre :
Indicateurs de décharge et de perte de proportionnalité du chargement
Date :
06/04/98
Auteur(s) :
G. DEBRUYNE
Clé :
R4.20.01-A
Page :
1/6
Organisme(s) : EDF/IMA/MMN
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
Document : 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 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.
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
HI-74/98/008/A
Code_Aster ®
Version
4.0
Titre :
Indicateurs de décharge et de perte de proportionnalité du chargement
Date :
06/04/98
Auteur(s) :
G. DEBRUYNE
Clé :
R4.20.01-A
Page :
2/6
1 Introduction
1.1
Définition d'un chargement proportionnel
Considérant une structure soumise à un chargement thermo-mécanique dans l'intervalle de temps [0,t]
on dira que ce chargement est proportionnel (ou encore radial) au point matériel P si le champ de
contrainte représenté en ce point par le tenseur est proportionnel à un tenseur indépendant de
l'instant considéré, le coefficient de proportionnalité étant une fonction monotone du temps.
Formellement, cela s'exprimera par :
[0,t],(P,) =
(
)
( P),
( ) > 0
0
fonction monotone dans [0,t].
Cette définition implique, en particulier, que les directions principales des contraintes restent constantes,
au point considéré, tout au long du trajet de chargement (ces directions peuvent être bien sûr variables
d'un point à un autre).
1.2
Importance du chargement proportionnel et utilité d'indicateurs
Pour les matériaux plastiques, les champs mécaniques dépendent de toute l'histoire écoulée pendant
le trajet de chargement. Les lois d'écoulement ont donc un caractère incrémental et leur intégration
dépend de chaque cas de chargement. Une exception notable concerne justement le chargement
proportionnel pour lequel la loi d'écoulement peut être intégrée une fois pour toutes. Par exemple, la loi
de plasticité de Prandtl-Reuss basée sur le critère de Von Mises peut être remplacée par une loi
élastique non linéaire (appelée loi de Hencky-Mises). Les cas de chargement strictement proportionnel
sont assez rares. En effet, il faut remplir un grand nombre de conditions pour réaliser un tel cas [bib1]
et ces dernières sont peu souvent vérifiées pour les structures industrielles. On peut même dire que
pour des structures présentant des défauts géométriques tels que des fissures, ces conditions ne sont
jamais strictement vérifiées.
Quand les chargements sont multi-axiaux, cycliques, ou thermo-mécaniques transitoires, certains
tronçons du trajet de chargement peuvent être fortement non proportionnels. Il est alors utile de repérer
ces tronçons et d'évaluer l'importance de la perte de proportionnalité, afin par exemple d'ajuster la
discrétisation en temps du problème élastoplastique pour le tronçon considéré, ou de jauger la validité
de certains post-traitements (en mécanique de la rupture par exemple).
1.3
Différents types d'indicateurs de perte de proportionnalité
Il semble difficile de définir une grandeur unique et simple qui pourrait détecter à la fois des zones
spatiales de perte de proportionnalité et des tronçons de chemin de chargement (zones temporelles)
en un point matériel. C'est pourquoi nous proposons quatre grandeurs scalaires ayant chacune leur
spécificité : deux, définies par des champs mesurant en chaque point la décharge et la déviation des
contraintes entre deux pas de temps (indicateurs locaux), deux autres de nature plus globale,
caractérisant dans une zone donnée de la structure une histoire de chargement non proportionnelle.
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
HI-74/98/008/A
Code_Aster ®
Version
4.0
Titre :
Indicateurs de décharge et de perte de proportionnalité du chargement
Date :
06/04/98
Auteur(s) :
G. DEBRUYNE
Clé :
R4.20.01-A
Page :
3/6
2 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.
2.1
Indicateur de décharge
Cet indicateur mesure au point M et entre l'instant t et t + t
, la variation relative de la norme des
( M,t + t
) - ( M,t)
contraintes au sens de Von Misès. Il s'écrit formellement : I1 =
. Cette
( M,t + t
)
quantité est négative en cas de décharge locale au point M. La norme ( M ,t) peut s'écrire de
quatre façons différentes suivant le choix du modélisateur :
3
1)
( M,t)
D. D
=
, où D est la partie déviatorique du tenseur des contraintes (cette
2
norme est utile en plasticité avec écrouissage isotrope).
3
2)
( M,t) =
. , où l'on considère la totalité du tenseur des contraintes afin de détecter par
2
exemple les diminutions de pression hydrostatique.
3
3)
( M,t)
( D X ).( D
=
-
- X) , avec X le tenseur de la contrainte de rappel dans le
2
cas d'une loi élastoplastique avec un écrouissage cinématique.
3
4)
( M,t) =
( - X ).( - X )
2
2.2
Indicateur de perte de radialité
Cet indicateur mesure au point M et entre l'instant t et t + t
, la variation de la direction des
contraintes. Il s'écrit :
( M,t).
I2 = 1-
, où le produit scalaire « . » est associé à l'une des quatre normes
( M,t)
précédentes. Cette quantité est nulle lorsque la radialité est conservée durant l'incrément de temps.
Remarque :
Ces indicateurs sont étroitement liés à la discrétisation en temps du problème. En particulier, si
cette discrétisation est trop grossière, on peut très bien ne pas détecter la décharge ou la perte de
radialité se produisant pendant l'incrément de temps.
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
HI-74/98/008/A
Code_Aster ®
Version
4.0
Titre :
Indicateurs de décharge et de perte de proportionnalité du chargement
Date :
06/04/98
Auteur(s) :
G. DEBRUYNE
Clé :
R4.20.01-A
Page :
4/6
3 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).
3.1
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 S du domaine ) 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 :
p
1
.
I =
1
( -
) d
3
S
(
,
+ R( p)) p
S
Y
où 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 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]).
(pt )= (pt
1
2 )
(pt1)
(pt2)
Y
Y + R(p)
I3 = 0
I3 0
Figure 3.1-Erreur! Argument de commutateur inconnu. : Chemin de chargement entre t1 et t2
dans le plan déviatorique des contraintes
Remarque 2 :
Dans la formulation de cet indicateur interviennent trois ingrédients liés à la plasticité :
· l'écart entre la direction des contraintes et des déformations plastiques actuelles
( . p
p
),
· la position des contraintes par rapport au seuil actuel (
(
+ ( ))
Y
R p ),
· l'écart entre la norme de la déformation plastique actuelle et la déformation plastique cumulée
( p p ).
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
HI-74/98/008/A
Code_Aster ®
Version
4.0
Titre :
Indicateurs de décharge et de perte de proportionnalité du chargement
Date :
06/04/98
Auteur(s) :
G. DEBRUYNE
Clé :
R4.20.01-A
Page :
5/6
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 p 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
S entourant le défaut avec un faible voisinage afin d'obtenir une valeur significative de
l'indicateur.
3.2 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 :
1
I =
1
( -
) d
4
S S
,
t
où est la densité d'énergie de déformation définie par : (t) =
.
! d , et est la densité
0
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 :
1
2
2
µ
((t)) = Ktr () +
2 , si < ( + R
Y
(p) ,
2
3
1
2
2
R (P)
P
((t)) = Ktr () +
+
R q
2
6
( )dq, si = ( + R
Y
(p) ,
µ
0
avec K le module de compressibilité, µ le coefficient de cisaillement de Lamé, R le seuil de la
courbe de traction associé à la norme de la déformation plastique P
p
= (celui-ci peut être différent
du vrai seuil plastique, car P 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 ( = ) .
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
HI-74/98/008/A
Code_Aster ®
Version
4.0
Titre :
Indicateurs de décharge et de perte de proportionnalité du chargement
Date :
06/04/98
Auteur(s) :
G. DEBRUYNE
Clé :
R4.20.01-A
Page :
6/6
4
Environnement d'utilisation des différents indicateurs
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.
4.1 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_ELEM présente les options `DCHA_ELGA_SIGM' et
`DCHA_ELNO_SIGM' pour l'indicateur de décharge I1 évalué aux noeuds ou aux points de Gauss de
l'élément, et les options `RADI_ELGA_SIGM' et `RADI_ELNO_SIGM' pour l'indicateur de perte de
radialité I2 évalué aux noeuds ou aux points de Gauss de l'élément [bib2]. Le mot-clé NORME permet
d'utiliser l'une des quatre normes décrites au paragraphe 2.1, 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.
4.2 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 I3 et I4 . 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).
5 Bibliographie
[1]
J. LEMAITRE, J.-L.CHABOCHE, Mécanique des matériaux solides, Dunod 1985.
[2]
J. PELLET, Manuel d'Utilisation du Code_Aster. Document [U4.61.02], 24/07/96.
Manuel de Référence
Fascicule R4.20 : Méthodes d'analyse
HI-74/98/008/A
Document Outline