Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
1/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
Organisme(s) :
EDF/IMA/MMN
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire des structures volumiques
Document : V7.22.101
HSNV101 - Thermo-plasticité et métallurgie
découplées en traction simple
Résumé :
On traite la détermination de l'évolution mécanique d'un barreau cylindrique soumis à des évolutions thermiques
T t
( )
et métallurgiques
Z t
( )
connues et uniformes (la transformation métallurgique est de type bainitique).
Les éléments utilisés sont des éléments axisymétriques et la relation de comportement est la plasticité de
von Mises avec écrouissage isotrope linéaire (pour la modélisation B, on tient également compte de la plasticité
de transformation).
La limite élastique et la pente de la courbe de traction dépendent de la température et de la composition
métallurgique.
Le coefficient de dilatation
dépend de la composition métallurgique.
Les transformations métallurgiques ont lieu à
!
p
=
0
(c'est en ce sens que le test découple la plasticité de
transformation de la plasticité classique).
Les résultats fournis par le Code_Aster sont très satisfaisants avec des erreurs inférieures à 2 %.
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
2/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
1
Problème de référence
1.1 Géométrie
C
D
h
A
B
a
z
r
P
Rayon : a = 0.05 m.
Hauteur : h = 0.2 m.
1.2
Propriétés des matériaux
(
)
( )
( ) ( )
( )
( )
(
)
(
)
E
s
T T
H t
t E t
E t
t
H
H
T T
s
H
s
T T
H
H
yaust
oaust
aust
o
oaust
aust
oaust
aust
o
fbm
aust
oaust
aust
yfbm
ofbm
fbm
o
aust
ref
ofbm
fbm
o
fbm
=
=
+
-
=
-
=
=
=
+
-
=
=
=
=
=
+
-
=
=
=
=
-
-
-
-
-
-
-
°
°
°
°
200000.10 Pa
notons
0.3
400. 10 Pa
15. 10
C
0.5 10 Pa. C
1250 10 Pa
23.5 10
C
5 10 Pa. C
2.52 10
530. 10 Pa
6
6
6
1
6
1
6
6
1
6
1
3
6
(
)
fbm
fbm
o
ref
fbm
ofbm
fbm
m
T T
T
C
s
H
cp
k
+
-
=
°
=
=
=
=
=
=
°
-
°
°
-
°
-
-
-
-
-
-
-
900
1
0.5 10 Pa. C
50 10 Pa
2 000 000 J.m
C
9999.9 W.m
C
5 10 Pa C
1. 10
Pa
6
1
6
3
1
-1
1
6
1
10
*
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
fbm
=
déformation des phases ferritique, bainitique et martensitique à la température de 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).
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
3/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
TRC pour modéliser une évolution métallurgique de type bainitique, sur toute la structure, de la forme :
Z
fbm
=
0.
t
-
1
2
-
1
1.
si
t
1
si
1
t
<
2
si
t
2
1
=
60 s
2
=
112 s
Loi de plasticité de transformation :
( )
!
!
pt
fbm
fbm
fbm
K
F Z
Z
=
avec
F Z
fbm
( )
=
Z
fbm
2
-
Z
fbm
(
)
1.3
Conditions aux limites et chargements
·
u
Z
= 0 sur le côté AB (condition de symétrie).
·
traction imposée sur le côté CD,
p t
( )
=
p
o
t
360 10
6
Pa
pour
t
1
pour
t
1
p
o
=
6 10
6
Pa
1
=
60 s
.
·
T
T
t
o
=
+ µ
,
µ
= -5°C.s
-1
sur toute la structure.
1.4 Conditions
initiales
T
o
=
900
°
C
=
T
ref
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
4/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
Avant transformation, solution élastique pour
t
<
1
.
t
( )
=
p
o
t
zz
t
( )
=
zz
e
t
( )
+
zz
th
t
( )
=
t
( )
E
+
aust
T
-
T
o
(
)
La limite élastique est atteinte pour
1
'
=
oaust
p
o
-
s
aust
×
µ =
47.06 s .
Avant transformation, solution thermo-élasto-plastique,
1
'
t
1
,
1
= 60 s.
t
( )
=
p
o
t
zz
t
( )
=
zz
e
t
( )
+
zz
th
t
( )
+
zz
p
t
( )
zz
e
t
( )
=
t
( )
E
zz
th
t
( )
=
Z
aust
×
aust
T
-
T
o
(
)
zz
p
t
( )
=
t
( )
-
yaust
+
S
aust
µ
t
(
)
H
oaust
+
aust
µ
t
Pendant la transformation, solution thermo-élasto-métallurgique,
1
<
t
<
2
,
2
= 112 s.
t
( )
=
360 10
6
Pa
zz
t
( )
=
zz
e
t
( )
+
zz
th
t
( )
+
zz
pt
t
( )
+
zz
p
60
( )
zz
th
t
( )
=
Z
aust
×
aust
T
-
T
o
(
)
+
Z
fbm
×
fbm
T
-
T
o
(
)
+
Z
fbm
×
ref
fbm
zz
p
t
( )
=
k
fbm
F Z
fbm
( )
P
o
1
Après la transformation, solution thermo-élasto-plastique,
2
<
t
<
3
,
3
= 176 s.
t
( )
=
360 10
6
Pa
zz
t
( )
=
zz
e
t
( )
+
zz
th
t
( )
+
zz
p
t
( )
+
zz
pt
112
( )
zz
p
t
( )
=
t
( )
-
ofbm
+
s
fbm
µ
t
(
)
H
ofbm
+
fbm
µ
t
2.2
Résultats de référence
zz
p
,
,
et
zz
pour t = 47, 48, 64 et 114 secondes.
zz
p
pour t = 60 et 176 secondes.
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.
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
5/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation
N1
N2
N3
N4
N5
N6
N7
N8
N9
N10
N11
N12
N13
A
B
C
D
y
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]
STAT_NON_LINE
COMP_INCR
RELATION
META_EP
[U4.32.01]
CALC_ELEM
OPTION
EPSI_ELNO_DEPL
[U4.61.01]
RECU_CHAMP
NOM_CHAM
VARI_ELNO_ELGA
[U4.62.01]
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
6/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
4
Résultats de la modélisation A
4.1 Valeurs
testées
Identification
Référence
Aster
% différence
zz
p
t = 47 s
0
0
0
t = 47 s
0
0
0
t = 47 s
282. 10
6
282. 10
6
0
zz
t = 47 s
4.1125 10
3
4.1125 10
3
0
zz
p
t = 48 s
3.2653 10
3
3.26537 10
3
0.011
t = 48 s
1
1
0
t = 48 s
288. 10
6
288. 10
6
0
zz
t = 48 s
9.3469 10
4
9.34645 10
4
0.005
zz
p
t = 60 s
0.04
0.04
0
zz
p
t = 64 s
0.040
4.0 10
-2
0
t = 64 s
0
0
0
t = 64 s
360. 10
6
360. 10
6
0
zz
t = 64 s
3.4683 10
2
3.46908 10
2
0.023
zz
p
t = 114 s
0.04107
4.10688 10
2
+0.004
t = 114 s
1
1
0
t = 114 s
360. 10
6
360. 10
6
0
zz
t = 114 s
0.03684
3.68407 10
2
0
zz
p
t = 176 s
0.06206
6.20680 10
2
0.000
4.2 Remarques
Dans cette modélisation :
zz
pt
T , Z
(
)
=
0
4.3 Paramètres
d'exécution
Version : 4.02.14
Machine : CRAY C90
Encombrement mémoire :
8 mégamots
Temps CPU User :
109.3 secondes
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
7/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
5 Modélisation
B
5.1
Caractéristiques de la modélisation
N1
N2
N3
N4
N5
N6
N7
N8
N9
N10
N11
N12
N13
A
B
C
D
y
A = N4, B = N5, C = N13, D = N12.
5.2
Caractéristiques du maillage
Nombre de noeuds : 13
Nombre de mailles et types : 2 mailles QUAD8, 6 mailles SEG3
5.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_PT
STAT_NON_LINE
COMP_INCR
RELATION
META_EP_PT
[U4.32.01]
CALC_ELEM
OPTION
EPSI_ELNO_DEPL
[U4.61.01]
RECU_CHAMP
NOM_CHAM
VARI_ELNO_ELGA
[U4.62.01]
Code_Aster
®
Version
4.0
Titre :
HSNV101 - Thermo-plasticité et métallurgie
Date :
01/12/98
Auteur(s) :
A. RAZAKANAIVO, F. WAECKEL
Clé :
V7.22.101-B
Page :
8/8
Manuel de Validation
Fascicule V7.22 : Thermo-mécanique statique non linéaire
HI-75/98/040/A
6
Résultats de la modélisation B
6.1 Valeurs
testées
Identification
Référence
Aster
% différence
zz
p
t = 47 s
0
0
0
t = 47 s
0
0
0
t = 47 s
282. 10
6
282. 10
6
0
zz
t = 47 s
4.1125 10
3
4.1125 10
3
0
zz
p
t = 48 s
3.2653 10
3
3.26535 10
3
0.011
t = 48 s
1
1
0
t = 48 s
288. 10
6
288. 10
6
0
zz
t = 48 s
9.3469 10
4
9.34644 10
4
0.005
zz
p
t = 60 s
0.04
0.04
0
zz
p
t = 64 s
0.04
4.0 10
2
0
t = 64 s
0
0
0
t = 64 s
360. 10
6
359.99 10
6
0.004
zz
t = 64 s
4.00085 10
2
4.000268 10
2
0.015
zz
p
t = 114 s
0.041071
4.10751 10
2
+0.004
t = 114 s
1
1
0
t = 114 s
360. 10
6
360.01 10
6
0.000
zz
t = 114 s
0.072841
7.144112 10
2
1.915
zz
p
t = 176 s
0.06206
6.2066 10
2
0.000
6.2 Remarques
Dans cette modélisation, on prend en compte le terme dû à la plasticité de transformation :
( )
!
,
!
pt
T Z
Z
0
0
lorsque
6.3 Paramètres
d'exécution
Version : 4.02.14
Machine : CRAY C90
Encombrement mémoire :
8 mégamots
Temps CPU User :
114.11 secondes
7
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.025% sauf pour la déformation à l'instant 114 s où l'erreur atteint 2% pour la modélisation
B.