Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
1/12
Organisme(s) : EDF/IMA/MMN
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
Document R7.02.03
Taux de restitution de l'énergie en thermo-élasticité
non-linéaire
Résumé :
On présente le calcul du taux de restitution de l'énergie par la méthode théta en 2D ou en 3D pour un problème
thermo-élastique non-linéaire. La relation de comportement élastique non linéaire est décrite en [R5.03.20].
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
2/12
Table des matières
1 Calcul du taux de restitution de l'énergie par la méthode théta en thermo-élasticité non linéaire ......... 3
1.1 Relation de comportement .............................................................................................................. 3
1.2 Energie potentielle et relations d'équilibre ....................................................................................... 4
1.3 Expression Lagrangienne du taux de restitution de l'énergie.......................................................... 4
1.4 Implantation de G en thermo-élasticité non linéaire dans le Code_Aster ....................................... 5
1.5 Mise en garde.................................................................................................................................. 7
2 Calcul du taux de restitution de l'énergie par la méthode théta en grandes transformations ................ 7
2.1 Relation de comportement .............................................................................................................. 7
2.2 Energie potentielle et relations d'équilibre ....................................................................................... 8
2.3 Expression Lagrangienne du taux de restitution de l'énergie en thermo-élasticité non-linéaire et
en grandes transformations ............................................................................................................ 9
2.4 Implantation dans le Code_Aster .................................................................................................. 11
2.5 Restriction...................................................................................................................................... 11
3 Bibliographie ........................................................................................................................................ 12
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
3/12
1
Calcul du taux de restitution de l'énergie par la méthode
théta en thermo-élasticité non linéaire
1.1
Relation de comportement
On considère un solide fissuré occupant le domaine de l'espace R2 ou R3. Soit :
·
u le champ de déplacement,
·
T le champ de température,
·
f le champ de forces volumiques appliquées sur ,
·
g le champ de forces surfaciques appliquées sur une partie S de ,
·
U le champ de déplacements imposés sur une partie Sd de .
S
f
g
Sd
Le comportement du solide est supposé élastique non-linéaire tel que la relation de comportement
coïncide avec la loi élastoplastique de Hencky-Von Mises (écrouissage isotrope) dans le cas d'un
chargement qui induit une évolution radiale et monotone en tout point. Ce modèle est choisi dans les
commandes CALC_G_THETA [U4.63.03] et CALC_G_LOCAL [U4.64.04] par l'intermédiaire du mot-clé
RELATION:'ELAS_VMIS_LINE' ou 'ELAS_VMIS_TRAC' sous le mot-clé facteur COMP_ELAS
[R5.03.20].
On désigne par :
·
le tenseur de déformations,
·
° le tenseur des déformations initiales,
·
le tenseur des contraintes,
·
° le tenseur des contraintes initiales,
·
(,, °,°°,,T) la densité d'énergie libre.
est relié au champ de déplacement u par :
1
u
( ) =
u
(
)
2 i, j + uj,i
La densité d'énergie libre (,
, °, °°,,T) est une fonction convexe et différentiable, connue pour un
état donné [R5.03.20 éq 3]. La relation de comportement du matériau s'écrit sous la forme :
ij =
(, °, °,T)
ij
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
4/12
Elle dérive du potentiel énergie libre. Pour cette relation de comportement hyperélastique, on sait
donner un sens au taux de restitution d'énergie dans le cadre de l'approche globale en mécanique de la
rupture. Ce n'est pas le cas pour une relation de comportement plastique.
1.2
Energie potentielle et relations d'équilibre
On définit les espaces des champs cinématiquement admissibles V et Vo .
V =
v
{ admissibles, v = U sur Sd}
Vo =
v
{ admissibles, v = O sur Sd}
Avec les hypothèses du paragraphe [§1.1] (avec ° = ° = 0 ), les relations d'équilibre en formulation
faible sont :
u
V
ij v
v
v
i, j d
=
fi
i d +
gi
S
i d , v Vo
Elles sont obtenues en minimisant l'énergie potentielle globale du système :
W v
( ) =
( (v),T)d- f
v
v
i
i d -
gi
S
i d
La démonstration est identique à celle en élasticité linéaire [R7.02.01 §1.2].
1.3
Expression Lagrangienne du taux de restitution de l'énergie
Soit m la normale unitaire à o située dans le plan tangent à dans .
o
m
Soit le champ tel que :
=
µ
{ tels que µn = 0 sur }
en notant n la normale à .
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
5/12
Le taux de restitution d'énergie G est solution de l'équation variationnelle :
G m dS
G () ,
=
o
où G () est défini par :
G () =
ij iu,p p,j - k,k -
T,k d
k
T
1
1
+
th
ij -
°
°
°
,
,
ij
2
ij kk - ij - ij -
ij
2
ij k kd
+
fi i
u k,k +
fi,k k i
u d
+
gi,k u + g u k,k -
nk
d
S
k
i
i
i
nk
- ij n U d
S
j
i,k
k
d
La démonstration est identique à celle du calcul de G en élasticité linéaire [R7.02.01]. L'expression est
la même, le post-traitement est donc identique.
1.4 Implantation de G en thermo-élasticité non linéaire dans le
Code_Aster
Les types d'éléments et de chargements, l'environnement nécessaire sont les mêmes que pour
l'implantation de G en thermo-élasticité linéaire [R7.02.01 §2.4].
Pour le calcul des différents termes de G
( ), à un état donné, on récupère la densité d'énergie libre
,
( , T), les déformations et les contraintes , calculées pour la relation de comportement
non-linéaire (routine NMELNL).
On suppose que ° = ° = 0 (terme identique en thermo-élasticité linéaire ou non-linéaire). La densité
d'énergie libre s'écrit alors [R5.03.20 §1.5] :
· en thermo-élasticité linéaire :
1
2
2µ
(,T) = K
3
2
2 (
- (T - T
kk
réf ) +
eq
3
avec
3
2
D D
3
1
1
eq =
ij ij =
ij - kkij ij
kkij
2
2
3
-
3
3
2
1
2
eq =
ijij - kk
2
3
· en thermo-élasticité non-linéaire (2µ eq y ) :
( ) 1
2
1
2
3
2 (
(
)
p
,T
K
T T
R p
eq
=
-
-
+
+
R s ds
kk
réf
6µ
( ( eq)
( )
( )
0
avec R( (
p eq ) : fonction d'écrouissage.
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
6/12
Pour un écrouissage isotrope linéaire (RELATION : `ELAS_VMIS_LINE') on a :
E E
R(p) = + p
T
= + a p
y
E - E
y
T
eq - y
E E
p =
avec a
T
=
µ
3 + a
E - ET
(
p
1
1
p
A p) =
R(s)ds = + a p2 = p
yp
y
+
0
( +a p
y
)
2
2
2
(
p
A p) = ( + R
y
(p)
2
Le post-traitement est ensuite identique au problème en élasticité linéaire sauf pour le terme
thermique :
THER = -
,
T k
T
k
Si les coefficients de LAME (T) et µ T
( ) sont indépendants de la température, ce terme est nul.
Dans le cas contraire, il faut calculer
( , T) à un instant donné.
T
Pour un écrouissage isotrope linéaire, on a :
(
1 dK T
d
T
,T)
( )
=
kk -
3 T - T
- 3K
réf
+
T - Tréf kk -
3 T - T
T
2 dT (
(
)
( )(
) (
réf )
dT
R(p)
dR(P) dµ(T)
dA p
+
2µ
-
R
2
(p)
( )
6µ
dT
dT
+ dT
dR(p) d y (T) d a(T)
d (
p T)
avec
=
+
p + a
dT
dT
dT
dT
d a(T)
1
dET (T)
d E
2
(T) 2
=
E -
E
dT
(
2
T
E - E
T )
dT
dT
d (
p T)
1
dµ T
dA T
d T
=
( -
y
eq )
( )
( )
3
+
( µ
3
a)
( )
dT
(
µ
2
y
3 + a)
dT
dT -
+
dT
dA(p) 1 (
dp T)
=
(
1
d T
dRp T
y
y + Rp )
( )
( )
+ p
+
dT
2 dT
2 dT
dT
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
7/12
1.5
Mise en garde
Attention ! Par définition, dans le cas général :
( , T) :
Bien qu'il soit possible d'effectuer un calcul élasto-plastique suivi du calcul de G en élasticité non
linéaire, il faut bien savoir que cela n'a aucun sens thermo-dynamique et qu'il est normal que le résultat
dépende du champ .
2
Calcul du taux de restitution de l'énergie par la méthode
théta en grandes transformations
On étend la relation de comportement du [§1] à de grands déplacements et de grandes rotations, dans
la mesure où elle dérive d'un potentiel (loi hyperélastique). Cette fonctionnalité est déclenchée par le
mot-clé DEFORMATION: 'GREEN' dans les commandes CALC_G_THETA [U4.63.03] et
CALC_G_LOCAL [U4.64.04].
2.1
Relation de comportement
On désigne par :
·
E le tenseur de déformations de Green-Lagrange,
·
S le tenseur des contraintes de Piola-Lagrange appelé encore deuxième tenseur de
Piola-Kirchoff,
·
E
( ) la densité d'énergie interne.
Le comportement du solide est supposé hyperélastique, à savoir que :
·
E est relié au champ de déplacement u mesuré par rapport à la configuration de référence
o par :
1
E
(
)
ij u
( ) =
u
2 i, j + uj,i + uk,i uk, j
·
S est relié au tenseur des contraintes de cauchy T par :
S
-1
-1
ij
= det F
( ) Fik Tkl Fjl
F étant le gradient de la transformation qui fait passer de la configuration de référence o à
la configuration actuelle , relié au déplacement par :
F
(
)
ij
= ij +ui,j
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
8/12
La relation de comportement d'un matériau hyperélastique s'écrit sous la forme :
Sij =
=
= S
E
ji
ij
Ejj
Cette relation décrit un comportement élastique non-linéaire, analogue à celle du [§1.1]. Elle offre la
possibilité de traiter les problèmes de mécanique de la rupture sans y intégrer la plasticité. Et dans le
cas d'un chargement radial monotone, elle permet d'obtenir des déformations et des contraintes de la
structure semblables à celles que l'on obtiendrait si le matériau présentait un écrouissage isotrope. Le
matériau hyperélastique a un comportement mécanique réversible, c'est-à-dire que tout cycle de
chargement n'engendre aucune dissipation.
Ce modèle et choisi dans les commandes CALC_G_THETA [U4.63.03] et CALC_G_LOCAL [U4.64.04]
par l'intermédiaire du mot-clé :
RELATION : 'ELAS'
pour une relation élastique "linéaire", c'est-à-dire que la relation entre les déformations et les
contraintes considérées est linéaire,
RELATION : 'ELAS_VMIS_LINE' ou 'ELAS_VMIS_TRAC'
pour une relation de comportement élastique "non linéaire" (loi de HENCKY-VON MISES à
écrouissage isotrope linéaire).
Une telle relation de comportement permet en toute rigueur de prendre en compte de grandes
déformations et de grandes rotations. Toutefois, on se cantonne à de petites déformations pour assurer
l'existence d'une solution et pour être identique à un comportement élastoplastique sous un
chargement radial monotone [R5.03.20 §2.1].
2.2
Energie potentielle et relations d'équilibre
Le chargement considéré se réduit à une densité surfacique non suiveuse R appliquée sur une partie
o du bord de o (hypothèse des charges mortes [R5.03.20 §2.2]).
On définit un espace des champs cinématiquement admissible V :
V
=
v
{ admissibles, v = 0 sur o}
Les relations d'équilibre en formulation faible sont :
F
ik Skj vi,j d =
R
i vi d , vV
o
Elles peuvent être obtenues en minimisant l'énergie potentielle globale du système :
W v
( ) =
( v())d - R
v
i i d
o
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
9/12
En effet, si cette fonctionnelle est minimale pour le champ de déplacement u , alors :
W =
E
v
ij d -
Ri i d = 0 , v V
o E
ij
=
S
(
)
ij
1 v
d - R
v
i,j + vj,i + vp,i up, j + u p,i vp, j
i
i d
o
2
=
S
(
)
ij
v
v
ip + up,i
p,i d -
Ri i d
o
=
F
pi S
v
ij vp,j d -
Ri i d
o
=
F
ik S
v
kj vi, j d -
Ri i d = 0 , v V
o
Nous retrouvons donc les équations d'équilibre et la relation de comportement en ayant posé :
Sij =
E
ij
2.3 Expression Lagrangienne du taux de restitution de l'énergie en
thermo-élasticité non-linéaire et en grandes transformations
Par définition, le taux de restitution d'énergie G est défini par l'opposé de la dérivée de l'énergie
potentielle à l'équilibre par rapport au domaine [bib1]. Il est calculé par la méthode théta, qui est une
méthode lagrangienne de dérivation de l'énergie potentielle [bib4] et [bib2]. On considère des
transformations F : M M + M
( ) du domaine o en un domaine qui correspondent à
des propagations de la fissure. A ces familles de configuration de référence ainsi définies
correspondent des familles de configurations déformées où la fissure s'est propagée. Le taux de
restitution de l'énergie G est alors l'opposé de la dérivée de l'énergie potentielle W u
(
( )) à l'équilibre
par rapport à l'évolution initiale du fond de fissure :
W u
(
( ))
G = -
d
=0
On note comme dans [bib 4] par . la dérivation lagrangienne dans une propagation virtuelle de fissure
de vitesse . Soit (, M) un champ quelconque ( réel positif et M appartenant au domaine
o), nous noterons :
(, M) = (
, F (M)) et
! = =0
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
10/12
L'énergie potentielle définie sur est ramenée sur o , R est supposé indépendant de , la
dérivation par rapport au paramètre de propagation est alors aisée et le taux de restitution de
l'énergie dans cette propagation est solution de l'équation variationnelle :
G m dS
G () ,
=
o
avec :
" #
$ %
$
- G () = ((E,T) + (E,T)k, d
k
- Ri !ui + Ri,k u
k
i + R u
i
i k ,k -
nk
d
o
nk
Or :
(" #
$ %
$
(E,T) =
!E +
!T
E
ij
T
ij
Par la suite, nous ne considérerons que le terme
!
E , le terme thermique étant traité de la même
E
ij
ij
façon qu'en petits déplacement [R7.02.01].
Et d'après la proposition 2 de [bib4] :
1
!Eij = 2(!ui,j + !uj,i + !uk, u
i
k , + u
j
k ,i !
uk, j )
1
-
2 (ui,p p, + u
j
j,p
p, + u
i
k , p
p, u
i
k , + u
j
k , u
i
k , p
p, j )
On peut éliminer !u de l'expression de G comme en petites déformations en remarquant que !u est
cinématiquement admissible (cf. [bib3] pour les problèmes de régularité) et en utilisant l'équation
d'équilibre :
"#%
((E))
d -
R u!
d
i
i
=
o
1 u! , +u! , +u! u
,
, + u
u
, !
d
,
o
E 2 ( i j
j i
k i
k j
k i
k j )
ij
1
- R u
i !i
d -
u ,
, + u ,
, + u
u
,
,
, + u
u
d
,
,
,
o
E 2 ( i p p j
j p p i
k p p i
k j
k i
k p p j )
ij
= -
S u ,
, + u
u
d
ij( i p p j k i, k,p p,j)
o
= -
S
(
k i, )
ij
+ u
u
d
ki
k ,p
p, j
o
= -
S
ij kFi uk,p p,j d
o
= - ik
F Skj i
u ,p p, j d
o
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
11/12
Finalement, on obtient :
G () =
F S
, ,
E ,
,
,
(u
)- ( ) d + R u + R u k k -
d
ik
kj
i p
p j
k k
i k
k
i
i
i
o
n
k
L'expression complète pour les chargements suivants :
· densité surfacique non suiveuse R appliquée sur une partie du bord de o,
· densité volumique non suiveuse f appliquée sur le domaine ,
et en tenant compte de la thermique :
G () = F S
ik
kj
E
(ui,p p,j)- ( ) k,k - ,T d
o
T
k
k
+ f u
i
i k ,k + fi, u d
k
k
i
o
+ Ri, u
k
k
i + R u
i
i k ,k -
nk d
nk
2.4
Implantation dans le Code_Aster
La comparaison des formules de G
( ) du [§1.3] et du [§1.4] montre que les termes de G
( ) sont très
proches. L'introduction des grandes transformations nécessite peu de modification en post-traitement.
La présence du mot-clé DEFORMATION:'GREEN' sous le mot-clé facteur COMP_ELAS des
commandes CALC_G_THETA et CALC_G_LOCAL indique qu'il est nécessaire de récupérer le tenseur
des contraintes de Piola-Lagrange S et le gradient de la transformation F (routines NMGEOM et
NMELNL).
Les types d'éléments finis sont les mêmes qu'en élasticité linéaire [R7.02.01 §2.4]. Ce sont les
éléments isoparamétriques 2D et 3D.
Les chargements supportés sont ceux supportés en élasticité linéaire à condition que ce soient des
charges mortes : typiquement une force imposée est une charge morte tandis que la pression est un
chargement suiveur puisqu'il dépend de l'orientation de la surface, donc de la transformation.
2.5 Restriction
Avec la relation de comportement précisée au 2, on a une formulation de G valable pour de grandes
déformations pour des matériaux hyper-élastiques, mais... si l'on souhaite une cohérence avec le
matériau réel qui, rappelons-le, est élasto-plastique, il est impératif de se cantonner à des déformations
petites, les déplacements et les rotations pouvant être grandes.
Les conditions de chargements proportionnels et monotones, indispensables pour assurer la
cohérence du modèle avec le matériau réel, conduisent à des restrictions importantes du champ des
problèmes à même d'être traités par cette méthode (le thermique en particulier peut conduire à des
décharges locales). Il ne peut donc s'agir que d'une solution palliative avant d'être en mesure de
donner un sens au taux de restitution d'énergie dans le cadre de comportements plastiques.
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Code_Aster ®
Version
4.0
Titre :
Taux de restitution de l'énergie en thermo-élasticité non-linéaire
Date :
08/10/97
Auteur(s) :
E. VISSE
Clé :
R7.02.03-B
Page :
12/12
3 Bibliographie
[1]
BUI H.D., Mécanique de la rupture fragile, Masson, 1977.
[2]
DESTUYNDER Ph., DJAOUA M., Sur une interprétation de l'intégrale de Rice en théorie de la
rupture fragile, Mathematics Methods in the Applied Sciences, Vol. 3, pp. 70-87, 1981.
[3]
GRISVARD P., "Problèmes aux limites dans les polygones", Mode d'emploi - EDF - Bulletin de
la Direction des Etudes et Recherches, Série C, 1, 1986 pp. 21-59.
[4]
MIALON P. , "Calcul de la dérivée d'une grandeur par rapport à un fond de fissure par la
méthode théta", EDF - Bulletin de la Direction des Etudes et Recherches, Série C, n3, 1988,
pp. 1-28.
[5]
MIALON P., Etude du taux de restitution de l'énergie dans une direction marquant un angle
avec une fissure, note interne EDF, HI/4740-07, 1984.
[6]
SIDOROFF F., Cours sur les grandes déformations, Ecole d'été, Sophia-Antipolis, 8 au
10 septembre 1982.
[7]
LORENTZ E., Une relation de comportement hyperélastique non linéaire, Note interne EDF
DER HI-74/95/011/0, 1995.
Manuel de Référence
Fascicule R7.02 : Mécanique de la rupture
HI-75/97/025/A
Document Outline