background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
1/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
Organisme(s) :
EDF-R&D/AMA















Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
Document : V6.04.172



SSNV172 ­ Comportements viscoplastiques
monocristallins




Résumé :

Ce test permet de valider les comportements mono-cristallins dans une situation uni-axiale. La géométrie
traitée est un cube, l'état de contraintes et de déformations est homogène. Afin de tester au mieux les lois de
comportement, indépendamment de la définition des systèmes de glissement, on utilise ici un système de
glissement particulier, non physique, qui représente un glissement dans une seule direction. Ceci permet de
comparer deux des comportements viscoplastiques mono-cristallins, (avec écrouissage cinématique défini par
une variable de rappel) aux comportements viscoplastiques macroscopiques de Chaboche. Pour les autres
comportements on vérifie seulement la non régression des résultats.

Tous ces tests sont réalisés dans une seule modélisation.
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
2/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
1
Problème de référence
1.1 Géométrie
On définit un élément MA1, contenant les noeuds P1, P2, P3, P4, P5, P6, P7 et P8.
1.2
Propriétés de matériaux
Comportement élastique avec :
Module d'Young :
MPa
E 145200
=
Coefficient de Poisson :
3
.
0
=
Premier calcul (référence)
Comportement VISC_CIN1_CHAB
CIN1_CHAB=_F(R_0=75.5 MPa
R_I=85.27 MPa
B=19.34,
C_I=10.0 MPa
K=1.0,
W=0.0,
G_0=36.68,
A_I=1.0,),
LEMAITRE=_F(N=10.0,
UN_SUR_K=0.025 Mpa
-1
UN_SUR_M=0.0,),

Deuxième calcul : mono-cristallin, avec système de glissement UNIAXIAL
Le comportement du monocristal est défini par :
Type d'écoulement : ECOU_PLAS1 dont les paramètres sont :
MPa
K
n
MPa
c
40
,
10
,
10
=
=
=
X
Y
Z
P1
P1
P2
P3
P4
P5
P6
P7
P8
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
3/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
Type d'écrouissage isotrope : ECRO_ISOT1 dont les paramètres sont :
MPa
R
5
.
75
0
=
34
.
19
=
b
MPa
Q
77
.
9
=
0
=
h
Type d'écrouissage cinématique : ECRO_CINE1 dont les paramètres sont :
68
.
36
=
d
La famille des systèmes de glissement est : UNIAXIAL
Les paramètres utilisés ici, en uni-axial, correspondent à ceux utilisés pour VISC_CIN1_CHAB, en
remarquant simplement que
0
R
R
Q
I
-
=
.
Troisième calcul :
Le comportement du monocristal est défini par :
Type d'écoulement : ECOU_VISC2 dont les paramètres sont :
0
10
10
40
10
10
=
=
=
=
=
a
d
MPa
c
MPa
k
n
Type d'écrouissage isotrope : ECRO_ISOT2 dont les paramètres sont :
5
.
75
0
=
R
34
.
19
1
=
b
0
2
=
b
77
.
9
1
=
Q
0
2
=
Q
0
=
h
Type d'écrouissage cinématique : ECRO_CINE2 dont les paramètres sont :
68
.
36
=
d
10
10
=
M
0
=
m
0
=
c
La famille des systèmes de glissement est : UNIAXIAL
Là encore, les paramètres utilisés ici, en uni-axial, correspondent à ceux utilisés pour
VISC_CIN1_CHAB, en remarquant simplement que Q=R_I-R_0.
Quatrième calcul :
Le comportement du monocristal est identique à celui du deuxième calcul (la différence vient de la
résolution : implicite dans ce cas, et explicite dans le deuxième calcul).
Cinquième calcul :
Le comportement du monocristal est cette fois multiple. De plus son orientation par rapport au
repère global de définition des coordonnées est définie par un angle de 30 degrés autour de Z.
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
4/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
Pour la famille de systèmes de glissement BASAL, les lois sont :
Type d'écoulement : ECOU_PLAS1 dont les paramètres sont :
MPa
K
n
MPa
c
40
,
10
,
10
=
=
=
Type d'écrouissage isotrope : ECRO_ISOT1 dont les paramètres sont :
MPa
R
5
.
75
0
=
34
.
19
=
b
MPa
Q
77
.
9
=
0
=
h
Type d'écrouissage cinématique : ECRO_CINE1 dont les paramètres sont :
68
.
36
=
d
Pour la famille de systèmes de glissement PRISMATIQUE , le comportement du monocristal est
défini par :
Type d'écoulement : ECOU_VISC2 dont les paramètres sont :
0
10
10
40
10
10
=
=
=
=
=
a
d
MPa
c
MPa
k
n
Type d'écrouissage isotrope : ECRO_ISOT2 dont les paramètres sont :
5
.
75
0
=
R
34
.
19
1
=
b
0
2
=
b
77
.
9
1
=
Q
0
2
=
Q
0
=
h
Type d'écrouissage cinématique : ECRO_CINE2 dont les paramètres sont :
68
.
36
=
d
10
10
=
M
0
=
m
0
=
c
Pour la famille de systèmes de glissement CUBIQUE1 , le comportement du monocristal est défini
par :
Type d'écoulement : ECOU_VISC3 dont les paramètres sont :
1
10
10
40
0
=
=
=
=
V
MPa
k
µ
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
5/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
Type d'écrouissage isotrope : ECRO_ISOT1 dont les paramètres sont :
5
.
75
0
=
R
34
.
19
=
b
77
.
9
=
Q
0
=
h
Type d'écrouissage cinématique : ECRO_CINE2 dont les paramètres sont :
68
.
36
=
d
10
=
M
5
.
1
=
m
1
=
c
Pour la famille de systèmes de glissement MACLAGE , le comportement du monocristal est défini
par :
Type d'écoulement : ECOU_VISC3 dont les paramètres sont :
1
10
10
40
0
=
=
=
=
V
MPa
k
µ
Type d'écrouissage isotrope : ECRO_ISOT2 dont les paramètres sont :
5
.
75
0
=
R
34
.
19
1
=
b
0
2
=
b
77
.
9
1
=
Q
0
2
=
Q
0
=
h
Type d'écrouissage cinématique : ECRO_CINE1 dont les paramètres sont :
68
.
36
=
d

1.3
Conditions aux limites et chargements
·
Noeud P4
:
0
=
= DY
DX
·
Noeud P8
:
0
=
=
=
DZ
DY
DX
·
Noeuds P2 et P6
:
0
=
DX
·
Noeuds P1, P3, P5 et P7
:
25
=
FX
Le chargement est croissant de
0
=
FX
à
755
.
0
*
25
=
FX
N, en un incrément, ce qui conduit à un
état de contraintes uniaxial de 75.5 MPa (limite de linéarité)
Le chargement croit ensuite jusqu'à
955
.
0
*
25
=
FX
N en
100
incréments.
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
6/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
2
Solution de référence
2.1
Méthode de calcul
Le premier calcul (comportement VISC_CIN1_CHAB) sert de référence. Les valeurs obtenues au
dernier incrément sont (en n'importe quel point de la maille) :
Contrainte SIXX
-9.55E+01
Déformation totale
EPXX
-6.7594E-04
Déformation plastique EPXX
-1.823E-05
Les deuxième, troisième et quatrième calculs seront validés par comparaison à ce résultat. Par contre,
le cinquième calcul n'a pas de valeur de référence. Les tests sont donc de non-régression.


3 Modélisation
A
3.1
Caractéristiques du maillage
Nombre de noeuds : 8.
Modélisation 3D : 1 élément de volume quadratique :
HEXA8
.
La structure contient un seul grain, l'état de contraintes et de déformations est homogène.
Les deuxième, troisième et quatrième calculs seront validés par comparaison à ce résultat.
Le deuxième et troisième calculs utilisent une intégration explicite (RUNGE_KUTTA).
Le quatrième calcul utilise une intégration implicite.
Le cinquième calcul utilise une intégration implicite, l'orientation choisie est
)
0
,
0
,
30
(
.


4
Résultats de la modélisation A
4.1 Fonctionnalités
testées
Commandes Options
DEFI_MATERIAU ECOU_VISC1
ECOU_VISC2
ECOU_VISC3
ECRO_ISOT1
ECRO_ISOT2
ECRO_CINE1
ECRO_CINE2
DEFI_COMPOR MONOCRISTAL FAMI_SYST_GLIS
STAT_NON_LINE COMP_INCR
RELATION='MONOCRISTAL'
COMPOR
RESO_INTE
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
7/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
4.2 Valeurs
testées
Deuxième calcul (explicite, ECOU_VISC1, ECRO_ISOT1, ECRO_CINE1)

Identification Référence
Aster %
différence
xx
de SIEF_ELGA
-95.5 -95.5
0
xx
de
EPSI_ELGA_DEPL
-6.7594E-04
-6.7575E-04 0.03
xx
de EPSP_ELGA
-1.823E-05 -1.804E-05
1.04


Troisième calcul (explicite, ECOU_VISC2, ECRO_ISOT2, ECRO_CINE2)

Identification Référence
Aster %
différence
xx
de SIEF_ELGA
-95.5 -95.5
0
xx
de EPSI_ELGA_DEPL
-6.7594E-04 -6.7575E-04
0.03
xx
de EPSP_ELGA
-1.823E-05 -1.804E-05
1.04

Quatrième calcul (implicite, ECOU_VISC1, ECRO_ISOT1, ECRO_CINE21)

Identification Référence
Aster %
différence
xx
de SIEF_ELGA
-95.5 -95.5
0
xx
de EPSI_ELGA_DEPL
-6.7594E-04
-6.7624E-04 0.04
xx
de EPSP_ELGA
-1.823E-05 -1.8528E-05
1.6

Cinquième calcul (implicite, tous les comportements, 4 familles de systèmes de glissement)
Test de non régression

Identification Référence
Aster %
différence
xx
de SIEF_ELGA
-95.5 -95.5
0
xx
de EPSI_ELGA_DEPL
-6.7571E-04 -6.7571E-04
0
background image
Code_Aster
®
Version
7.4
Titre :
SSNV172 ­ Comportements viscoplastiques monocristallins
Date
:
02/11/05
Auteur(s) :
J.M. PROIX, T. KANIT
Clé
:
V6.04.172-A
Page :
8/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
5
Synthèse des résultats
Les résultats obtenus sont en bon accord avec la solution de référence. La différence constatée
(1.6% au maximum) provient des schémas numériques différents, sensibles à la discrétisation
temporelle.