background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
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.171



SSNV171 - Comparaison POLY_CFC et
MONOCRISTAL




Résumé :

On présente ici un test ayant pour référence un autre modèle de comportement déjà validé dans Code_Aster,
qui est POLY_CFC. La géométrie traitée est un cube. Le but est de tester le concept MONOCRISTAL en utilisant
une loi de comportement équivalente à POLY_CFC.
background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
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
Module d'Young :
MPa
E 145200
=
Coefficient de Poisson :
3
.
0
=
Modélisation A :
Un monocristal de type POLY_CFC, orienté suivant
)
0
,
0
,
30
(
avec :
0
=
L
D
0
=
A
D
10
=
N
40
=
K
5
.
75
0
=
34
.
19
1
=
B
0
=
L
H
77
.
9
1
=
Q
27
.
33
2
-
=
Q
345
.
5
2
=
B
00001
.
0
1
=
C
68
.
36
1
=
D
0
2
=
C
X
Y
Z
P1
P1
P2
P3
P4
P5
P6
P7
P8
background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
Page :
3/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
Modélisation B :
ECOU_VISC2 avec :
TYPE_PARA='ECOU_VISC'
10
=
N
40
=
K
0
=
C
68
,
36
=
D
10
=
A
ECRO_ISOT2 avec :
TYPE_PARA='ECRO_ISOT'
5
.
75
0
=
R
34
.
19
1
=
B
345
.
5
2
=
B
77
.
9
1
=
Q
27
.
33
2
-
=
Q
0
=
H
ECRO_CINE1 avec :
TYPE_PARA='ECRO_CINE'
68
.
36
=
D
L'orientation choisie est la même que celle de la modélisation A,
)
0
,
0
,
30
(
.

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 :
35
=
FX
Le chargement est croissant de
0
=
FX
à
35
=
FX
, divisé en
100
incréments.


2
Solution de référence
2.1
Méthode de calcul
On compare la fonctionnalité de nouveaux concepts avec un concept déjà implémenté dans
Code_Aster, qui est POLY_CFC.

background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
Page :
4/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
3 Modélisation
A
3.1
Caractéristiques du maillage
Nombre de noeuds : 8.
Modélisation 3D : 1 élément de volume quadratique : HEXA8.
Texture propre a POLY_CFC, avec une famille de systèmes de glissement correspondant à la structure
cubique à faces centrées, en l'occurrence le système OCTAEDRIQUE. La structure contient un seul
grain, donc une fraction volumique égale à 1 et l'orientation choisie est
)
0
,
0
,
30
(
.


4
Résultats de la modélisation A
4.1 Fonctionnalités
testées

Commandes Options
DEFI_TEXTURE SYST_GLISSEMENT
PLAN
N,L
ANGL_NAUT,
PROPORTION
DEFI_MATERIAU POLY_CFC
TEXTURE,
DL,DA,N,K,TAU_0,B1,HL,Q1,Q2,B2,C1,D1,C2
STAT_NON_LINE COMP_INCR
RELATION='POLY_CFC'

4.2 Valeurs
testées

Identification Référence
Aster %
différence
xx
de SIEF_ELGA
-210 -210
-7.38E-05
xx
de EPSI_ELGA_DEPL
-0.0018955 -0.0018955
-4.30E-04
yy
de EPSI_ELGA_DEPL
0.000503377
0.000503377 9.11E-05
xx
de EPSP_ELGA
-0.000449212
-0.000449212 5.43E-05
yy
de EPSP_ELGA
0.0000694943
0.00006949337 -3.46E-05

background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
Page :
5/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
5 Modélisation
B
5.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 (fraction volumique égale à 1) et l'orientation choisie est
)
0
,
0
,
30
(
.


6
Résultats de la modélisation B
6.1 Fonctionnalités
testées

Commandes Options
DEFI_MATERIAU ELAS
ECOU_VISC2
ECRO_ISOT2
ECRO_CINE1
E,NU
TYPE_PARA='ECOU_VISC',
N,K,C,D,A
TYPE_PARA='ECRO_ISOT',
R_0,B1,B2,Q1,Q2,H
TYPE_PARA='ECRO_CINE',
D
DEFI_COMPOR MONOCRISTAL
MATER,
ECOULEMENT='ECOU_VISC2',
ECRO_ISOT='ECRO_ISOT2',
ECRO_CINE='ECRO_CINE1',
ELAS='ELAS',
FAMI_SYST_GLIS='OCTAEDRIQUE'
AFFE_CARA_ELEM MODELE
MASSIF
GROUP_MA='TOUT',
ANGL_REP
STAT_NON_LINE CARA_ELEM
COMP_INCR
RELATION='MONOCRISTAL',
COMPOR

6.2 Valeurs
testées

Identification Référence
Aster %
différence
xx
de SIEF_ELGA
-210 -210
-7.38E-05
xx
de EPSI_ELGA_DEPL
-0.00189 -0.00189
-4.30E-04
yy
de EPSI_ELGA_DEPL
0.0005033 0.0005033
9.15E-05
xx
de EPSP_ELGA
-0.000449 -0.000449
5.66E-05
yy
de EPSP_ELGA
0.0000695 0.0000695
-3.23E-05
background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
Page :
6/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
7 Modélisation
C
7.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 (fraction volumique égale à 1) et l'orientation choisie est
)
0
,
0
,
30
(
.
Le but de cette modélisation est la comparaison entre les valeurs des champs de contraintes et de
déformations obtenues avec le monocristal en implicite et en explicite ainsi que la comparaison pour
les deux modélisations des valeurs de la variable interne V43 qui est la déformation plastique cumulée
globale du monocristal.


8
Résultats de la modélisation C
8.1 Fonctionnalités
testées

Commandes Options
DEFI_MATERIAU ELAS
ECOU_VISC1
ECRO_ISOT1
ECRO_CINE1
E,NU
TYPE_PARA='ECOU_VISC',
N,K,C
TYPE_PARA='ECRO_ISOT',
R_0,B1,B2,Q1,Q2,H
TYPE_PARA='ECRO_CINE',
D
DEFI_COMPOR MONOCRISTAL
MATER,
ECOULEMENT='ECOU_VISC1',
ECRO_ISOT='ECRO_ISOT1',
ECRO_CINE='ECRO_CINE1',
ELAS='ELAS',
FAMI_SYST_GLIS='OCTAEDRIQUE'
AFFE_CARA_ELEM MODELE
MASSIF
GROUP_MA='TOUT',
ANGL_REP
STAT_NON_LINE CARA_ELEM
COMP_INCR
CONVERGENCE
RELATION='MONOCRISTAL',
COMPOR
RESO_INTE='IMPLICITE'
STAT_NON_LINE CARA_ELEM
COMP_INCR
CONVERGENCE
RELATION='MONOCRISTAL',
COMPOR
RESO_INTE='RUNGE_KUTTA_2'
background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
Page :
7/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
8.2 Valeurs
testées
En implicite
Identification Référence Aster %
différence
xx
de SIEF_ELGA
-210 -210
-5.97E-05
xx
de EPSI_ELGA_DEPL
-0.00189
-0.00189
1E-2
yy
de EPSI_ELGA_DEPL
0.0005056
0.0005059 6E-3
yy
de EPSP_ELGA
0.0000686
0.0000687 4.7E-3
43
V
de VARI_ELGA
0.0004846
0.0004846 5.8E-8
En explicite
Identification Référence Aster %
différence
xx
de SIEF_ELGA
-210 -210
-8.24E-05
xx
de EPSI_ELGA_DEPL
-0.00188
-0.00188 2.25E-4
yy
de EPSI_ELGA_DEPL
0.0005021
0.0005021 1.83E-6
yy
de EPSP_ELGA
0.0000682
0.0000682 2.22E-4
43
V
de VARI_ELGA
0.0004813
0.0004813 5.21E-4
background image
Code_Aster
®
Version
8.2
Titre :
SSNV171 - Comparaison POLY_CFC et MONOCRISTAL
Date :
02/11/05
Auteur(s) :
J. EL GHARIB, T. KANIT, J.M. PROIX
Clé
:
V6.04.171-B
Page :
8/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/05/005/A
9
Synthèse des résultats
Les résultats obtenus sont en parfait accord avec la solution de référence.