r5 03 21

Code_Aster
Version
default
Titre : Modlisation lasto(visco)plastique avec crouissa[...] Date : 19/01/2011 Page : 1/16
Responsable : Renaud BARGELLINI Cl : R5.03.21 Rvision : 5129
Modlisation lasto(visco)plastique avec
crouissage isotrope en grandes dformations
Rsum
On dcrit ici une relation de comportement thermolastoplastique avec crouissage isotrope crite en grandes
dformations et propose par Simo et Miehe. Ce modle est disponible dans la commande STAT_NON_LINE
par lintermdiaire du mot-cl RELATION : VMIS_ISOT_TRAC, VMIS_ISOT_PUIS ou
VMIS_ISOT_LINE sous le mot-cl facteur COMP_INCR et avec le mot-cl DEFORMATION :
SIMO_MIEHE. Une version visqueuse de ce modle est galement propose : VISC_ISOT_TRAC et
VISC_ISOT_LINE.
Ce modle est implant pour les modlisations tridimensionnelles (3D), axisymtrique (AXIS) et en
dformations planes (D_PLAN).
On prsente lcriture et le traitement numrique de cette loi, ainsi que la formulation variationnelle associe. Il
sagit dune formulation variationnelle eulrienne, avec ractualisation de la gomtrie et qui tient compte de la
rigidit de comportement et de la rigidit gomtrique.
Manuel de rfrence Fascicule r5.03 : Mcanique non linaire
Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Version
default
Table des matires
1 Introduction .......................................................................................................................................... 3
2 Notations ............................................................................................................................................. 4
3 Rappels sur les grandes dformations ................................................................................................ 5
3.1 Cinmatique .................................................................................................................................. 5
3.2 Contraintes .................................................................................................................................... 6
3.3 Objectivit ..................................................................................................................................... 6
4 Prsentation du modle de comportement ......................................................................................... 7
4.1 Aspect cinmatique ...................................................................................................................... 7
4.2 Relations de comportement .......................................................................................................... 8
4.3 Choix de la fonction dcrouissage ............................................................................................... 9
4.4 Contraintes et variables internes ................................................................................................ 10
4.5 Domaine dutilisation ................................................................................................................... 10
4.6 Intgration de la loi de comportement ........................................................................................ 10
5 Formulation variationnelle ................................................................................................................. 13
5.1 Cas du milieu continu ................................................................................................................. 13
5.2 Discrtisation par lments finis ................................................................................................. 13
5.3 Expression de la matrice tangente du comportement ................................................................ 14
6 Bibliographie ...................................................................................................................................... 16
7 Description des versions du document .............................................................................................. 16
Code_Aster
Version
default
1 Introduction
Nous prsentons ici une loi de comportement thermolastoplastique crite en grandes dformations
propose par SIMO J.C. et MIEHE C. [bib1] qui tend en petites dformations vers le modle de
comportement lastoplastique crouissage isotrope et critre de Von Mises, dcrit en [R5.03.02].
Les choix cinmatiques permettent, comme avec la simple ractualisation disponible via le mot-cl
PETIT_REAC, de traiter des grands dplacements et des grandes dformations mais galement des
grandes rotations de manire exacte.
Les spcificits de ce modle sont les suivantes :
tout comme en petites dformations, on suppose lexistence dune configuration relche, cest-
-dire localement libre de contrainte, qui permet de dcomposer la dformation totale en une
partie thermolastique et une partie plastique,
la dcomposition de cette dformation en des parties thermolastique et plastique nest plus
additive comme en petites dformations (ou pour les modles grandes dformations crits en
taux de dformation avec par exemple une drive de Jaumann) mais multiplicative,
les dformations lastiques sont mesures dans la configuration actuelle (dforme) tandis que
les dformations plastiques sont mesures dans la configuration initiale,
comme en petites dformations, les contraintes dpendent uniquement des dformations
thermo-lastiques,
les dformations plastiques se font volume constant. La variation de volume est alors
uniquement due aux dformations thermo-lastiques,
ce modle conduit lors de son intgration numrique un modle incrmentalement objectif (cf.
[3.3]) ce qui permet dobtenir la solution exacte en prsence de grandes rotations.
Une version visqueuse de ce modle est galement disponible (loi en sinus hyperbolique comme dans
le cas du modle de Rousselier ROUSS_VISC, cf. [R5.03.07]).
Par la suite, on rappelle brivement quelques notions de mcanique en grandes dformations, puis on
prsente les relations de comportement du modle et son intgration numrique pour traiter les
quations dquilibre.
On propose une formulation variationnelle eulrienne, avec ractualisation de la gomtrie. A ce titre,
on exprime le travail des efforts intrieurs et sa variation (dans le but dune rsolution par la mthode
de Newton) pour le problme continu, qui fournissent respectivement aprs discrtisation par
lments finis le vecteur des forces intrieures et la matrice tangente.
Nota Bene :
On trouvera dans [bib2] ou [bib3] une prsentation approfondie sur les grandes dformations.
Ce document est extrait de [bib4] o lon fait une prsentation plus dtaille du modle
lastoplastique, de son intgration numrique et o lon donne quelques exemples de
validation.
Code_Aster
Version
default
2 Notations
On notera par :
Id
matrice identit
tr A trace du tenseur A
A
T
transpos du tenseur A
det A dterminant de A
X partie positive de X
A

partie dviatorique du tenseur A dfinie par
A=A(
1
3
tr A) Id
:
produit doublement contract :
A: B=
i , j
A
ij
B
ij
=tr ( AB
T
)

produit tensoriel : ( AB)
ijkl
=A
ij
B
kl
A
eq

valeur quivalente de von Mises dfinie par
A
eq
=
.
3
2

A:

A
X
A
gradient :
X
A=
A
X
div
x
A
divergence : ( div
x
A)
i
=
j
A
ij
x
j
\, j, E ,+ , K
coefficients de llasticit isotrope
c
y

limite dlasticit
o
coefficent de dilatation thermique
T
temprature
T
ref
temprature de rfrence
Par ailleurs, dans le cadre dune discrtisation en temps, toutes les quantits values linstant
prcdent sont indices par
, les quantits values linst

t +At
ant ne sont pas indices et les
incrments sont dsigns par A . On a ainsi :
AQ=QQ

Code_Aster
Version
default

3 Rappels sur les grandes dformations
3.1 Cinmatique
Considrons un solide soumis des grandes dformations. Soit D
0
le domaine occup par le solide
avant dformation et
D( t )
le domaine occup linstant t par le solide dform.
Confi gurati on actuel l e dforme
Confi gurati on i ni ti al e
F
0
( ) t

Figure 3.1-a : Reprsentation de la configuration initiale et dforme
Dans la configuration initiale D
0
, la position de toute particule du solide est dsigne par X
(description lagrangienne). Aprs dformation, la position linstant t de la particule qui occupait
la position X avant dformation est donne par la variable x (description eulrienne).
Le mouvement global du solide est dfini, avec u le dplacement, par :
x= x( X,t )=X+u
Pour dfinir le changement de mtrique au voisinage dun point, on introduit le tenseur gradient de la
transformation F :
F=
x
X
=Id+
X
u
Les transformations de llment de volume et de la masse volumique valent :
d D=Jd D
o
avec J=det F=
j
o
j

o j
o
et
j
sont respectivement la masse volumique dans les configurations initiale et actuelle.
Diffrents tenseurs de dformations peuvent tre obtenus en liminant la rotation dans la
transformation locale. Par exemple, en calculant directement les variations de longueur et dangle
(variation du produit scalaire), on obtient :
E=
1
2
( CId) avec C=F
T
F
A=
1
2
( Idb
1
) avec b=FF
T
E et A sont respectivement les tenseurs de dformation de Green-Lagrange et dEuler-
Almansi et C et b , les tenseurs de Cauchy-Green droit et gauche respectivement.
Code_Aster
Version
default
En description lagrangienne, on dcrira la dformation par les tenseurs C ou E car ce sont
des quantits dfinies sur D
0
, et en description eulrienne par les tenseurs b ou A (dfinis
sur
D
).

Remarque :
Soit un solide subissant deux transformations successives, par exemple la premire
transformation fait passer le solide de la configuration initiale D
0
une configuration D
1

(tenseur gradient F
1/ 0
et vecteur dplacement u
1/ 0
), puis la seconde transformation de la
configuration D
1
D
2
(tenseur gradient F
2/ 1
et vecteur dplacement u
2/ 1
).
0

1

2
F
1/0
F
2/1
F
2/0

Le passage de la configuration D
0
D
2
est donn par le tenseur gradient F
2/ 0

(dplacement u
2/ 0
=u
2/ 1
+u
1/ 0
) tel que :
F
2/ 0
=F
2/ 1
F
1/ 0

On obtient alors, par exemple, pour le tenseur de dformation de Green-Lagrange E
E
2/ 0
=F
1/ 0
T
E
2/ 1
F
1/ 0
+E
1/ 0

o E
2/ 0
, E
1/ 0
et E
2/ 1
sont les dformations de Green-lagrange des configurations D
2

par rapport D
0
associes F
2/ 0
, D
1
par rapport D
0
associes F
1/ 0
et D
2
par
rapport D
1
associes F
2/ 1
, respectivement.
Ceci constitue une des difficults rencontres lors de lcriture dune loi de comportement en
grandes dformations car on ne peut plus crire une formule analogue celle crite en
petites dformations, savoir e
2/ 0
=e
2/ 1
+e
1/ 0
o
e
est le tenseur de dformation
totale linaris.
Pour retrouver e
2/ 0
=e
2/ 1
+e
1/ 0
en petites dformations partir de lexpression de E
2/ 0
, il
faut ngliger tous les termes dordre 2 de u
2/ 0
, u
1/ 0
et u
2/ 1
. Dans ce cas, on a
E
2/ 0
e
2/ 0
, E
1/ 0
e
1/ 0
et F
1/ 0
T
E
2/ 1
F
1/ 0
e
2/ 1
.
3.2 Contraintes
Pour le modle dcrit ici, le tenseur des contraintes utilis est le tenseur eulrien de Kirchhoff
t

dfini par :
J c=t
o
c
est le tenseur eulrien de Cauchy. Le tenseur
t
rsulte donc dune mise lchelle
par la variation de volume du tenseur de Cauchy
c
; ceci nest pas le cas dautres tenseurs de
contraintes utiliss (premier et second tenseur de Piola-Kirchhoff).
En description eulrienne, les quations dquilibre sont donnes par :
div
x
c+jf =0 sur D
c. n=t
d
sur D
f

Code_Aster
Version
default
o f est la force volumique applique sur le domaine
D
, n la normale extrieure la
frontire
D
f
et
D
f
la partie de la frontire du domaine
D
o sont appliques les forces
surfaciques
t
d
.
3.3 Objectivit
Lorsquon crit une loi de comportement en grandes dformations, on doit vrifier que cette loi est
objective, cest--dire invariante par tout changement de rfrentiel spatial de la forme :
x
*
=c( t )+Q( t ) x
o Q est un tenseur orthogonal qui traduit la rotation du rfrentiel et c un vecteur qui traduit
la translation.
Plus concrtement, si on ralise un essai de traction dans la direction e
1
, par exemple, suivi dune
rotation de 90 autour de e
3
, ce qui revient effectuer un essai de traction selon e
2
, alors le danger
avec une loi de comportement non objective est de ne pas retrouver un tenseur des contraintes
uniaxial dans la direction e
2
(ce qui est notamment le cas avec la cinmatique PETIT_REAC).
4 Prsentation du modle de comportement
4.1 Aspect cinmatique
Ce modle suppose, tout comme en petites dformations, lexistence dune configuration relche
D
r
, cest--dire localement libre de contrainte, qui permet alors de dcomposer la dformation totale
en des parties lastique et plastique, cette dcomposition tant multiplicative.
Par la suite, on notera par F le tenseur gradient qui fait passer de la configuration initiale D
0
la
configuration actuelle
D( t )
, par F
p
le tenseur gradient qui fait passer de la configuration D
0
la
configuration relche
D
r
, et F
e
de la configuration D
r

D( t )
. Lindice p se rfre la
partie plastique, lindice e la partie lastique.
Confi gurati on i ni ti al e
Confi gurati on actuel l e
Confi gurati on rel che
F
e F
p
F
0

( ) t
T T
ref =
=
0

Figure 4.1-a : Dcomposition du tenseur gradient F en une partie lastique F
e
et plastique F
p
Par composition des mouvements, on obtient la dcomposition multiplicative suivante :
F=F
e
F
p

Code_Aster
Version
default
Les dformations lastiques sont mesures dans la configuration actuelle avec le tenseur eulrien de
Cauchy-Green gauche b
e
et les dformations plastiques dans la configuration initiale par le tenseur
G
p
(description lagrangienne). Ces deux tenseurs sont dfinis par :
b
e
=F
e
F
eT
, G
p
=( F
pT
F
p
)
1
do b
e
=FG
p
F
T

Le modle prsent est crit de facon distinguer les termes isochores des termes de changement de
volume. On introduit pour cela les deux tenseurs suivants :
F=J
1/3
F
et
b
e
=J
2/ 3
b
e
avec J=det F
Par dfinition, on a : det

F=1 et
det
b
e
=1
.
4.2 Relations de comportement
La loi prsente est un modle thermolasto(visco)plastique avec crouissage isotrope qui tend sous
lhypothse des petites dformations vers le modle [R5.03.02] avec critre de Von Mises (il sagit du
modle plastique). Les dformations plastiques se font volume constant si bien que :
J
p
=det F
p
=1
do
J=J
e
=det F
e
Les relations de comportement sont donnes par :
relation contrainte - dformation thermolastique :
t=j
b
e

tr t=
3K
2
( J
2
1)
9K
2
o( TT
ref
)( J+
1
J
)
seuil de plasticit (on admet quil sexprime avec les contraintes de Kirchhoff) :
f =t
eq
R( p)c
y

o R est la variable dcrouissage isotrope, fonction de la dformation plastique cumule p.
lois dcoulement :
F

G
p
F
T
=
\
3
t
eq
t
b
e
=3

\
(
1
3
tr
b
e
+
t
j
)
t
t
eq

p=
\

Pour le modle de plasticit, le multiplicateur plastique est obtenu en crivant la condition de
cohrence
f =0
et on a :
p0, f 0 et p f =0
Dans le cas visqueux, on prend p gale :
p=
e
0
|
sh
(
f
0
)
m

Code_Aster
Version
default
o e
0
, c
0
et m sont les coefficients de viscosit. Signalons que cette loi se rduit une loi de
type Norton lorsque les 2 paramtres matriaux e
0
et c
0
sont trs grands.
On rappelle que :
b
e
=J
2/ 3
b
e

F=J
1/3
F

et que la partition des dformations scrit :
b
e
=
FG
p
F
T

Pour des matriaux mtalliques o le rapport t
eq
/ j est petit devant 1, lexpression de la loi
dcoulement peut tre approche par :
F

G
p
F
T
=
\tr
b
e
t
t
eq
+O
(
t
eq
j
)
q.4.2-1
o
O
(
t
eq
j
)
est ngligeable devant le premier terme.
Cest cette dernire expression qui est implante dans le Code_Aster.
Remarque :
Si les dformations sont petites, on a :
J 1+tr e
b
e
Id+2e
e
G
p
Id2e
p

o
e
est la dformation totale,
e
e
la dformation lastique et
e
p
la dformation
plastique en petites dformations.
En remplaant ces trois expressions dans les quations de la loi de comportement prsente ici, on
retrouve bien le modle classique thermo-lasto-plastique avec crouissage isotrope et critre de Von
Mises.
4.3 Choix de la fonction dcrouissage
Cette relation de comportement est disponible dans loprateur STAT_NON_LINE, sous le mot-cl
facteur COMP_INCR et largument SIMO_MIEHE du mot-cl facteur DEFORMATION. On peut choisir
pour la fonction dcrouissage, un crouissage linaire ou bien fournir une courbe de traction. Cinq
relations peuvent tre utilises.
RELATION = / VMIS_ISOT_TRAC
/ VMIS_ISOT_PUIS
/ VMIS_ISOT_LINE
/ VISC_ISOT_TRAC
/ VISC_ISOT_LINE
Pour un comportement purement thermolastique, lutilisateur choisit largument VMIS_ISOT_LINE
par exemple, avec SY trs grand (le comportement est alors hyperlastique) ; pour un crouissage
isotrope donn par une courbe de traction, lutilisateur choisit largument VMIS_ISOT_TRAC dans le
Code_Aster
Version
default
cas plastique ou VISC_ISOT_TRAC dans le cas visqueux et pour un crouissage isotrope linaire,
largument VMIS_ISOT_LINE dans le cas plastique ou VISC_ISOT_LINE dans le cas visqueux.
Pour un comportement lastoplastique dont la courbe dcrouissage (courbe de traction rationnelle)
est donne par une loi en puissance, de la forme
R( p)=c
y
+c
y
(
E
ac
y
p
)
1
n
, lutilisateur choisit largument VMIS_ISOT_PUIS.
Les diffrentes caractristiques du matriau sont renseignes dans loprateur DEFI_MATERIAU
([U4.23.01]) sous les mots-cls :
ELAS quelque soit la loi (on donne le module dYoung, le coefficient de Poisson et
ventuellement le coefficient de dilatation thermique),
TRACTION pour VMIS_ISOT_TRAC et VISC_ISOT_TRAC (on donne la courbe de traction
rationnelle),
ECRO_PUIS pour VMIS_ISOT_PUIS (on donne les paramtres de la loi puissance),
ECRO_LINE pour VMIS_ISOT_LINE et VISC_ISOT_LINE (on donne la limite dlasticit
et la pente dcrouissage),
VISC_SINH pour VISC_ISOT_TRAC et VISC_ISOT_LINE (on donne les trois coefficients
de viscosit).
Remarque :
Lutilisateur doit bien sassurer que la courbe de traction exprimentale utilise, soit
directement, soit pour en dduire la pente dcrouissage est bien donne dans le plan
contrainte rationnelle
c=F/ S
- dformation logarithmique ln( 1+Al / l
0
) o l
0
est la
longueur initiale de la partie utile de lprouvette,
Al
la variation de longueur aprs
dformation, F la force applique et S la surface actuelle. On remarquera que
c=F/ S=
F
S
0
l
0
l
1
J
do t=J c=
F
S
0
l
0
l
. En gnral, cest bien la quantit
F
S
0
l
0
l
qui
est mesure par les exprimentateurs et ceci donne directement la contrainte de Kirchhoff
utilise dans le modle de Simo et Miehe.
4.4 Contraintes et variables internes
Les contraintes sont les contraintes de Cauchy
c
, calcules donc sur la configuration actuelle (six
composantes en 3D, quatre en 2D).
Les variables internes produites dans le Code_Aster sont :
V1, la dformation plastique cumule p ,
V2 V7 : loppos de la dformation lastique
b
e
V8 lindicateur de plasticit (0 si le dernier incrment calcul est lastique, 1 sinon).

Remarque :
Si lutilisateur veut rcuprer ventuellement des dformations en post-traitement de son
calcul, il faut tracer les dformations de Green-Lagrange E , qui reprsente une mesure
des dformations en grandes dformations (options EPSG_ELGA ou EPSG_ELNO ). Les
dformations linarises
e
classiques mesurent des dformations sous lhypothse des
petites dformations et nont pas de sens en grandes dformations.
4.5 Domaine dutilisation
Le choix dune cinmatique DEFORMATION : PETIT_REAC permet galement de traiter une loi de
comportement thermo-lastoplastique avec crouissage isotrope et critre de Von Mises en grandes
dformations. La loi est crite en petites dformations et la prise en compte des grandes dformations
se fait en ractualisant la gomtrie.
Entre la loi prsente ici (SIMO_MIEHE ) et PETIT_REAC,
Code_Aster
Version
default
il ny a pas de diffrence si les dformations sont petites
il ny a pas de diffrence si les dformations sont grandes mais les rotations petites
il y a des diffrences si les rotations sont importantes.
En particulier, la solution obtenue avec la cinmatique PETIT_REAC peut scarter notablement de la
solution exacte en prsence de grandes rotations et ce quelle que soit la taille des pas de temps
choisie par lutilisateur, contrairement la cinmatique SIMO_MIEHE.
4.6 Intgration de la loi de comportement
Dans le cas dun comportement incrmental, mot-cl facteur COMP_INCR, connaissant la contrainte
c
, la dformation plastique cumule

p
, la trace divise par trois du tenseur de dformations

lastiques
1
3
tr
b
e
, les dplacements
u
et
Au
et les tempratures
T
et T , on cherche
dterminer ( c , p ,
1
3
tr
b
e
) .
Les dplacements tant connus, les gradients de la transformation de D
0

D
, not
F
, et de
D

D
, not
AF
, sont connus.
La discrtisation implicite de la loi donne :
F=AFF

J=det F
F=J
1/3
F

b
e
=
FG
p
F
T

J c=t

t=j
b
e

1
3
tr t=
1
2
K ( J
2
1)
3
2
Ko( TT
ref
)( J+
1
J
)
f =t
eq
R( p
+A p)c
y

F( G
p
G
p
)
F
T
=tr
b
e t
t
eq
A p do
b
e
=
FG
p

F
T
tr
b
e t
t
eq
A p
Dans le cas plastique :
Ap0, f 0 et f A p=0
Dans le cas visqueux : t
eq
R( p
+Ap)c
y
c
0
sh
1
|
(
A p
e
0
At
)
1
m
=0
Remarque :
Cette formulation est incrmentalement objective car la seule quantit tensorielle
incrmentale qui intervient dans la discrtisation est
G
p
. Comme
G
p
et
G
p
sont
mesures sur la mme configuration, cest--dire la configuration initiale, la discrtisation de
G
p
, soit
AG
p
=G
p
G
p
, est incrmentalement objective.
On introduit
t
Tr
, le tenseur de Kirchhoff qui rsulte dune prdiction lastique (Tr : trial, en anglais
essai) :
t
Tr
=j
b
eTr

o
b
eTr
=
FG
p

F
T
=A
b
e
A
F
T
,
A
F=( AJ )
1/3
AF
et
AJ =det ( AF)
Code_Aster
Version
default
On obtient
b
e
partir des contraintes
t
par la relation contrainte - dformation thermolastique

et de la trace du tenseur des dformations lastiques.
b
e
=
t
+
1
3
tr
b
e

Remarque :
Lintrt de cette formulation est quil nest pas ncessaire de calculer la dformation
plastique
G
p
qui nous obligerait inverser le gradient de la transformation
F . On a
besoin uniquement de connatre
FG
p
F
T
.
Si t
eq
Tr
R( p
)+c
y
, on reste lastique. Dans ce cas, on a :
p=p
, t= t
Tr
+
1
3
tr t
Tr
Id et
1
3
tr
b
e
=
1
3
tr
b
eTr
sinon, on obtient :
tr
b
e
=tr
b
eTr
, grce la simplification sur la loi dcoulement :
b
e
=
b
e
Tr
tr
b
e
t
t
eq
Ap
En prenant les parties dviatoriques de cette quation, et en les multipliant par
j
on aboutit :
t
Tr
= t( 1+
jtr
b
eTr
A p
t
eq
)
En calculant la contrainte quivalente, on se ramne une quation scalaire non linaire en
A p
:
t
eq
Tr
t
eq
jtr

b
eTr
A p=0
Dans le cas plastique : t
eq
=c
y
+R( p
+Ap) , ce qui conduit

Ap
solution de lquation :
t
eq
Tr
c
y
R( p
+Ap )jtr
b
eTr
Ap=0
Dans le cas viscoplastique : t
eq
=c
y
+R( p
+A p)+c
0
sh
1
|
(
A p
e
0
At
)
1
m
, ce qui conduit
Ap

solution de lquation :
t
eq
Tr
c
y
R( p
+Ap )c
0
sh
1
|
(
A p
e
0
At
)
1
m
jtr

b
eTr
A p=0
Dans le cas o lcrouissage est linaire, ou bien donn par une courbe de traction dfinie point par
point, donc affine par morceaux, lquation rsoudre est linaire. On obtient directement la solution .
A p
Dans les autres cas, la rsolution est effectue dans Code_Aster par une mthode des
scantes avec intervalle de recherche (cf. [R5.03.05]). Lintgration peut tre contrle par les
paramtres RESI_INTE_RELA et ITER_INTE_MAXI.
Une fois
Ap
connu, on peut alors en dduire le tenseur de Kirchhoff, soit :
t=
t
eq
t
eq
Tr
t
Tr
+
|
K
2
( J
2
1)
3K
2
( TT
ref
)( J +
1
J
)
Id
Code_Aster
Version
default
Une fois calcule la dformation plastique cumule, le tenseur des contraintes et la matrice tangente,
on effectue une correction sur la trace du tenseur des dformations lastiques
b
e
pour tenir compte
de lincompressibilit plastique, qui nest pas conserve avec la simplification faite sur la loi
dcoulement [q 4.2.1]. Cette correction seffectue en utilisant une relation entre les invariants de
b
e
et
b
e
et en exploitant la condition dincompressibilit plastique
J
p
=1
(ou de manire
quivalente
det
b
e
=1
). Cette relation scrit :
x
3
J
2
e
x( 1
J
3
e
)=0
avec
J
2
e
=
1
2
b
eq
e
2
=
t
eq
2
2j
2
,
J
3
e
=det

b
e
=det
t
j
et x=
1
3
tr

b
e
La rsolution de cette quation du troisime degr permet dobtenir
tr
b
e
et par consquent la
dformation thermolastique
b
e
au pas de temps suivant. Dans le cas o cette quation admet
plusieurs solutions, on prend la solution la plus proche de la solution du pas de temps prcdent. Cest
dailleurs pourquoi on stocke dans une variable interne
1
3
tr
b
e
.
5 Formulation variationnelle
Dans la mesure o les contraintes fournies par la loi de comportement sont eulriennes, on choisit
une formulation variationnelle crite sur la configuration actuelle (eulrienne) et non sur la
configuration initiale, soit :
D
c
x
6v d D
_
F
int
. 6v
=
D
jf 6v d D+
D
f
t
d
6v dS
_
F
ext
. 6 v

6v
Cinmatiquement admissible
Nous nous intressons uniquement ici au travail des forces intrieures et sa variation dans loptique
dune rsolution par la mthode de Newton. On trouvera dans [bib 4] la dmonstration des
expressions prsentes.
5.1 Cas du milieu continu
On rcrit ici le travail des efforts intrieurs sous forme indicielle, soit :
F
int
.6 v=
D
c
ij
6v
i
x
j
d D
Nous avons besoin galement dexprimer la variation des efforts intrieurs dans la configuration
actuelle
D
soit :
6F
int
. 6u. 6 v=
D
|
c
ij
6u
p
x
p
c
ik
6u
j
x
k

6v
i
x
j
d D rigidit gomtrique
+
D
|
c
ij
AF
pq
6u
p
x
q

6 v
i
x
j
d D rigidit de comportement
i

o
x
sont les coordonnes dun point sur la configuration

D
.
Code_Aster
Version
default
5.2 Discrtisation par lments finis
On discrtise les dplacements u et les dplacements virtuels v par lments finis. Les
notations sont les suivantes, en adoptant la convention de sommation des indices rpts :
u
i
( x )=N
n
( x) U
i
n
u
i
x
j
=D
j
n
( x ) U
i
n
u
i
x
j
=D
j
n
( x) U
i
n

o :
N
n
( x )
est la fonction de forme associe au nud n
U
i
n
, la composante i du dplacement nodal du nud n
D
j
n
( x ) , les composantes du gradient des fonctions de forme sur la configuration
D
D
j
n
( x ) , les composantes du gradient des fonctions de forme sur la configuration D
On obtient pour le vecteur des forces intrieures :

( F
int
)
i
n
=
D
c
ij
D
j
n
d D

et pour la matrice tangente, qui nest pas a priori symtrique :
K
i p
nm
=
D
|
D
p
m
c
ij
D
j
n
D
k
m
c
ik
D
p
n
d D
+
D
|
D
q
n
c
ij
AF
pq
D
j
n
d D

Dans le cas dune modlisation bidimensionnelle (dformation plane), les expressions du vecteur des
forces intrieures et de la matrice tangente sont identiques ceci prt que les indices correspondants
aux composantes varient de 1 2 uniquement.
Dans le cas dune modlisation axisymtrique, en numrotant les axes dans lordre
(r , z , 0 ) ,
le
vecteur des forces intrieures scrit :
( F
int
axi
)
o
n
=
D
|
c
o
D
n
+c
33
N
n
r
6
o1
d D , o1,2 , 1,2
et la matrice tangente :
| K
axi
=| K +| K
corr

avec :
|
K(1 )
corr
1
nm
=
D
N
n
r
c
m
d D+
D
N
n
r
AF
33
D
m
d D
|
K( 2)
corr
o1
nm
=
D
D
o
n
c
33
N
m
r
d D+
D
D
- n
c
33
AF
o
N
m
r
d D
|
K( 3)
corr
11
nm
=
D
N
n
r
c
33
AF
33
N
m
r

Dun point de vue algorithmique, la matrice lmentaire tangente K nest pas symtrique a priori.
La rsolution globale se fera donc par dfaut avec un solveur non symtrique. Il est toutefois possible
de symtriser la matrice tangente globale avant rsolution (mot-cl SOLVEUR), ce qui permet
dconomiser du temps calcul mais peut dgrader la convergence.
Code_Aster
Version
default
5.3 Expression de la matrice tangente du comportement
On donne ici lexpression de la matrice tangente (option FULL_MECA au cours des itrations de
Newton, option RIGI_MECA_TANG pour la premire itration). Celle-ci est obtenue en linarisant le
systme dquations qui rgit la loi de comportement. Nous donnons ici le rsultat final de cette
linarisation. On trouvera dans [bib4] le dtail de ce calcul.
On pose :
J=det F ,
J
=det F
et
AJ =det AF
Pour loption FULL_MECA , on a :
A=
c
AF
=
( AJ )
1/3
J
H
1
3 J AJ
( HA
F)B
J
J
2
tB
+
J
J
|
KJ
3
2
K o( TT
ref
)( 1J
2
)
IdB

o B vaut :
B
11
=AF
22
AF
33
AF
23
AF
32
B
22
=AF
11
AF
33
AF
13
AF
31
B
33
=AF
11
AF
22
AF
12
AF
21
B
12
=AF
31
AF
23
AF
33
AF
21
B
21
=AF
13
AF
32
AF
33
AF
12
B
13
=AF
21
AF
32
AF
22
AF
31
B
31
=AF
12
AF
23
AF
22
AF
13
B
23
=AF
31
AF
12
AF
11
AF
32
B
32
=AF
13
AF
21
AF
11
AF
23

H et
HA
F
sont donns par :
Dans le cas lastique (f < 0) :
H
ijkl
=j
(
6
ik
b
lp
e
A
F
jp
+A
F
ip
b
pl
e
6
jk
)
2j
3
6
ij
A
F
kp
b
lp
e

et
HA
F=2j
b
eTr

Dans le cas plastique ( f =0 ) ou viscoplastique :
H
ijkl
=
j
a
( 6
ik
b
lp
e
A

F
jp
+A
F
ip
b
pl
e
6
jk
)
2j
|
6
ij
3a
+
RAp t
ij
t
eq
(

R+jtr
b
eTr
)
F
kp
b
lp
e
+
3j
2
tr
b
eTr
(

RA pt
eq
)
at
eq
3
(
R+mu tr
b
eTr
)
t
ij
t
kq
A
F
qp
b
lp
e

et
Code_Aster
Version
default
HA
F=
2j
a

b
eTr
2jtr
b
eTr
|
Id
3a
+

RA p t
t
eq
(
R+jtr
b
eTr
)
+
3j
2
tr

b
eTr
(
RA pt
eq
)
at
eq
3
(
R+jtr

b
eTr
)
( t:
b
eTr
) t
o a=
t
eq
Tr
t
eq

et
R=R' ( p)+c
0
(
1+
(
Ap
e
0
At
)
2
)
1
2
1
m
( e
0
At
)
1
m
( Ap)
1
m
1
_
uniquement cas visqueux
,
R
'
( p)
tant la drive de lcrouissage isotrope par rapport la dformation plastique cumule
p .
Pour loption RIGI_MECA_TANG, il sagit des mmes expressions que celles donnes pour
FULL_MECA mais avec
Ap=0
et avec toutes les variables et coefficients du matriau pris
linstant
t
(en principe, il faudrait dans le cas visqueux, prendre les expressions de FULL_MECA
dans le cas lastique, toutes les variables tant prises linstant
t
). En particulier, on aura
A
F=Id
.

6 Bibliographie
1) SIMO J.C., MIEHE C. : "Associative coupled thermoplasticity at finite strains :
Formulation, numerical analysis and implementation", Comp. Meth. Appl.
Mech. Eng., 98, pp. 41-104, North Holland, 1992.
2) SIDOROFF F. : "Formulations lasto-plastiques en grandes dformations", Rapport
Greco n29, 1981.
3) SIDOROFF F. : "Cours sur les grandes dformations", Rapport Greco n51, 1982.
4) CANO V., LORENTZ E. : "Introduction dans le Code_Aster dun modle de
comportement en grandes dformations lastoplastiques avec crouissage
isotrope", Note interne EDF DER, HI-74/98/006/0, 1998
7 Description des versions du document
Version
Aster
Auteur(s) ou contributeur(s),
organisme
Description des modifications
4.4 V.Cano, E.Lorentz EDF/R&D/AMA Texte initial
6.3 V.Cano, E.Lorentz EDF/R&D/AMA Fiche 6396
7.4 S.Michel-Ponnelle EDF/R&D/AMA Fiche 8000 : ajout de la viscosit
9.4 J.M.Proix EDF/R&D/AMA Ajout de VMIS_ISOT_PUIS.

r5 03 21

Transféré par

Droits d'auteur :

Formats disponibles

r5 03 21

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

r5 03 21

Transféré par

Droits d'auteur :

Formats disponibles

Code_Aster

, les quantits values linst

, la dformation plastique cumule

, la trace divise par trois du tenseur de dformations

par la relation contrainte - dformation thermolastique

+Ap) , ce qui conduit

sont les coordonnes dun point sur la configuration

On obtient pour le vecteur des forces intrieures :

Vous aimerez peut-être aussi