Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
1/8
Organisme(s) : EDF/IMA/MMN
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
Document : V7.22.104

HSNV104 - Thermo-plasticité et métallurgie en
déformations planes avec restauration
d'écrouissage

Résumé :
On traite la détermination de l'évolution mécanique d'un parallélipipède rectangle en déformations planes
soumis à des évolutions thermique T t
( ) et métallurgique Z t() connues et
uniformes (la transformation métallurgique est de type bainitique).
Les éléments utilisés sont des éléments bidimensionnels en déformations planes et la relation de comportement
est la plasticité de von Mises avec écrouissage isotrope linéaire. On tient compte de la restauration
d'écrouissage, mais non de la plasticité de transformation.
Le coefficient de dilatation dépend de la composition métallurgique.
La solution de référence est obtenue par la résolution analytique du problème.
Les résultats fournis par le Code_Aster sont très satisfaisants avec des erreurs inférieures à 0,8%.
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
2/8
1
Problème de référence
1.1 Géométrie
y
Largeur : a = 0.05 m.
Hauteur : h = 0.2 m.
C
D
h
r
x
A
B
a
1.2
Propriétés des matériaux
T
6
-6
-1
E T
E
fbm
= 200000.10 Pa
= 20.10 ° C
notons H(t)
( ) ( )
= E(T)-(T)
aust
-6
-1
o = 20.10 °C
= 0.3
Haust =
6
-
2000.10 Pa
3
ref
= 2.52 10
fbm
T ref =
°
900 C
H fbm = 2000.106 Pa
cp =
-
2000000 J.m 3
-
° C 1
= 9999.9 W.m
-1
-
° C 1
*aust
=
caractéristiques relatives à la phase austénitique,
*fbm
=
caractéristiques relatives aux phases ferritique, bainitique et martensitique,
fbm = coefficient de dilatation thermique des phases ferritique, bainitique et martensitique,
aust = coefficient de dilatation de la phase austénitique
ref
=
déformation des phases ferritique, bainitique et martensitique à la température de
fbm
référence, l'austénite étant considérée comme non déformée à cette température :
traduit la différence de compacité entre les structures cristallographiques cubiques à
faces centrées (austénite) et cubiques centrées (ferrite).
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
3/8
TRC pour modéliser une évolution métallurgique de type bainitique, sur toute la structure, de la forme :
0.

si t 1
1 = 60 s
t -
Z
1
fbm
=
si


1 t < 2
2 = 112s
2 - 1

si t
1.
2

Loi de plasticité de transformation : ! pt
fbm
= K
F(Z fbm) !Z fbm
avec F Z
( )
(
)
fbm = Z fbm Z - Z fbm
on ne tient pas compte de la plasticité de transformation donc on prend K fbm = 0
Notations : T
( 1)= T1
T
( 2)= T2
1.3
Conditions aux limites et chargements
·
uY = 0 sur le côté AB ; uX = 0 en A.
·
T
T o
=
+ µt , µ = ­5°C.s­1 sur toute la structure.
· Le chargement sur la structure est dû aux phénomènes de dilatation thermique et
métallurgique contraints dans la direction z par la condition de déformations planes.
1.4 Conditions
initiales
T o = 900°C = Tref
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
4/8
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
Avant transformation, solution thermo-élastique jusqu'en t1 tel que en t1 :
- aust

y
zz = - E th = y
T - To =
= -100°C
E + a
t1 = 20s
donc pour t t1
zz = -E T - To
(
)
Avant transformation, et pour t t1, solution thermo-élasto-plastique telle que :



th
p
zz = zz + zz + = 0
E

p

zz = R0zz + y
aust
1
1
y
0
p
-
d' où (
+ ) =
- (T - T ) et zz = p =
R aust
E
aust
aust
0
R0
R0
Pendant la transformation, on est en régime élastique, on a donc une solution thermo-élastique avec
changement de phase.
= -E (T - T0) + Z
p
réf
+ (
fbm
zz
1)
[
]
Après la transformation, on a toujours une solution thermo-élastique jusqu'en t2 .
En t2 :
(
)
zz = R T ,Z , eff
(
)+ y T,Z
A cause de la restauration d'écrouissage et du fait qu'on était en régime élastique pendant toute la
transformation : R = 0 avant replastification.
On a donc en t2 :
fbm + E(
+ p (

p
fbm
y
réf fbm
zz
1)
[
]
zz = - E (T - T 0 ) + réf
+ (
(T - T0) = -
fbm
zz
1)
[
]=y
E
(T - T0) = -624°C
t2 125s
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
5/8
Pour T < 276°C on a une solution thermo-élasto-plastique telle que :


p
zz = th + + (t)
E
zz

p
p
zz = R0 zz (t) - zz (1)


[
]+y
fbm
1
1

d' où (
+ ) = y
-(T - T0) -
- p (
R fbm
E
fbm
réf fbm
zz
1)
0
R0
2.2
Résultats de référence

eff
eff
zz ,
et à t = 60s

eff
eff
zz ,
et à t = 89s

eff
eff
zz ,
et à t = 112s

eff
eff
zz ,
et à t =176s
2.3 Bibliographie
[1]
DONORE A.M. - WAECKEL F. - Influence des transformations structurales dans les lois de
comportement élasto-plastiques Note HI-74/93/024.
[2]
DONORE.A.M. - WAECKEL.F. - RAZAKANAIVO.A. - Doc. Aster [R4.04.02].
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
6/8
3 Modélisation
3.1
Caractéristiques de la modélisation
y
C
D
N13
N11
N12
N9
N10
N7
N6
N8
N1
N3
N4
N2
N5
A
B
A = N4, B = N5, C = N13, D = N12.
3.2
Caractéristiques du maillage
Nombre de noeuds : 13.
Nombre de mailles et types : 2 mailles QUAD8, 6 mailles SEG3.
3.3 Fonctionnalités
testées
Commandes
Clés
DEFI_MATERIAU
META_THER
[U4.23.01]
THER_LINEAIRE
OPTION
META_ELGA_TEMP
[U4.23.05]
DEFI_MATERIAU
META_MECA_FO
[U4.23.01]
META_RE
STAT_NON_LINE
COMP_INCR
RELATION
META_EP_RE
[U4.32.01]
CALC_ELEM
OPTION
EPSI_ELNO_DEPL
[U4.61.01]
RECU_CHAMP
NOM_CHAM
VARI_ELNO_ELGA
[U4.62.01]
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
7/8
4 Résultats
4.1 Valeurs
testées
Identification
Référence
Aster
% différence
zz t = 60s
4.0792E8
4.080E8
0.02%
eff
t = 60s
3.9604E­3
3.9599E­3
­0.02%
eff
t = 60s
0.
0.
0.%
zz t = 89s
7.0684E8
7.0152E8
­0.752%
eff
t = 89s
3.9604E­3
3.9599E­3
­0.01%
eff
t = 89s
0.
0.
0.%
zz t = 112s
9.4392E8
9.44E8
0.008%
eff
t = 112s
0.
0.
0.
eff
t = 112s
0.
0.
0.
zz t = 176s
12.101E8
12.102E8
0.012%
eff
t = 176s
0.
0.
0.
eff
t = 176s
5.068921E­3
5.0688E­3
­0.002%
4.2 Remarques
Dans cette modélisation :
ptzz T
( , Z) = 0
L'erreur sur zz à 89 secondes provient en fait de l'erreur commise sur la description numérique de la
transformation métallurgique qui est, à cet instant, d'environ 56%.
4.3 Paramètres
d'exécution
Version : 4.02.14
Machine : CRAY C90
Encombrement mémoire :
8 mégamots
Temps CPU User :
188.88 secondes
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Code_Aster ®
Version
4.0
Titre :
HSNV104 - Thermo-plasticité et métallurgie en déformations planes
Date :
01/12/98
Auteur(s) :
F. WAECKEL, A. RAZAKANAIVO
Clé :
V7.22.104-A Page :
8/8
5
Synthèse des résultats
Les résultats trouvés avec le Code_Aster sont très satisfaisants, avec des pourcentages d'erreur
inférieurs à 0.8%.
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
HI-75/98/040/A

Document Outline