Code_Aster
®
Version
6.1
Titre :
SSNV152 - Traction élastique. Calcul des contraintes de Cauchy
Date
:
19/08/02
Auteur(s) :
P. MASSIN, G. BERTRAND
Clé
:
V6.04.152-A
Page :
1/6
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
HT-66/02/001/A
Organisme(s) :
EDF/AMA, CS SI
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
Document : V6.04.152
SSNV152- Traction élastique. Calcul des contraintes
de Cauchy
Résumé
Le but de ce test est de valider le calcul des contraintes de Cauchy dans la commande
CALC_ELEM
par l'option
SIGM_ELNO_COQUE
.
Code_Aster
®
Version
6.1
Titre :
SSNV152 - Traction élastique. Calcul des contraintes de Cauchy
Date
:
19/08/02
Auteur(s) :
P. MASSIN, G. BERTRAND
Clé
:
V6.04.152-A
Page :
2/6
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
HT-66/02/001/A
1
Problème de référence
1.1 Géométrie
La géométrie de ce test est une plaque carrée dans le plan (x,y) tournée de 30° par rapport à x autour
de z.
local
x
global
x
On appelle
l
la longueur de la plaque déformée, on notera x, y, z, les coordonnées de la configuration
déformée et X,Y,Z, les coordonnées de la configuration initiale
1.2
Propriétés des matériaux
On prend
0
et
.
000
200
=
=
MPa
E
1.3
Conditions aux limites et chargements mécaniques
On bloque les noeuds No1, No5 et No2 de sorte que DX=DY=DZ=DRX=DRY=DRZ=0,
et on impose un déplacement local Dx=100. sur les noeuds No3, No4 et No7.
.
1000
0
=
l
No2
No1
No5
No4
No7
No3
M1
z
y
x
Code_Aster
®
Version
6.1
Titre :
SSNV152 - Traction élastique. Calcul des contraintes de Cauchy
Date
:
19/08/02
Auteur(s) :
P. MASSIN, G. BERTRAND
Clé
:
V6.04.152-A
Page :
3/6
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
HT-66/02/001/A
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
La solution de référence est analytique.
Passage de l'état initial à l'état déformé :
Z
b
b
z
Y
a
a
y
X
l
x
0
0
0
,
,
1
=
=
=
où
a est la longueur de la déformée de la plaque suivant Y,
a
0
est la longueur initiale de la plaque,
b est l'épaisseur de la plaque déformée,
b
0
est l'épaisseur initiale de la plaque.
Du fait que
0
=
et des hypothèses de coque, on a
0
0
,
b
b
a
a
=
=
Tenseur Green-Lagrange :
Par définition du tenseur de Green-Lagrange, on a
+
+
=
j
k
i
k
i
j
j
i
ij
X
u
X
u
X
u
X
u
E
2
1
Avec
X
l
l
l
X
x
u
0
0
-
=
-
=
, on a donc
(
)
-
=
-
+
-
+
-
=
2
0
2
0
2
2
0
2
0
0
0
11
2
1
2
1
l
l
l
l
l
l
l
l
l
l
l
l
E
En remplaçant , on a
105
.
0
1000
1000
1100
2
1
2
2
2
11
=
-
=
E
Gradient de déformation :
Par définition :
=
=
1
0
0
0
1
0
0
0
0
l
l
dZ
dz
dY
dz
dX
dz
dZ
dy
dY
dy
dX
dy
dZ
dx
dY
dx
dX
dx
F
Soit
0
det
l
l
F
J
=
=
Code_Aster
®
Version
6.1
Titre :
SSNV152 - Traction élastique. Calcul des contraintes de Cauchy
Date
:
19/08/02
Auteur(s) :
P. MASSIN, G. BERTRAND
Clé
:
V6.04.152-A
Page :
4/6
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
HT-66/02/001/A
Contraintes de Piola-Kirchhoff de seconde espèce :
Soit S la contrainte de PK2, dans notre cas,
21000
105
.
0
200000
.
11
11
=
×
=
=
E
E
S
Contrainte de Cauchy
Soit s le tenseur de contraintes de Cauchy, on a la relation
(
)
T
F
S
F
F
s
.
.
det
1
=
, on en déduit alors
que
23100
21000
.
1000
1100
.
.
.
1
11
0
0
11
0
0
=
=
=
=
S
l
l
l
l
S
l
l
l
l
s
xx
2.2
Résultats de référence
On calcule des déplacements DX et DY au noeud NO3, les contraintes de PK2 et les contraintes de
Cauchy sur la maille M1.
2.3
Incertitude sur la solution
Résultat analytique.
2.4 Références
bibliographiques
Néant.
Code_Aster
®
Version
6.1
Titre :
SSNV152 - Traction élastique. Calcul des contraintes de Cauchy
Date
:
19/08/02
Auteur(s) :
P. MASSIN, G. BERTRAND
Clé
:
V6.04.152-A
Page :
5/6
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
HT-66/02/001/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation
On utilise des éléments COQUE_3D
3.2
Caractéristiques du maillage
Les coordonnées des principaux noeuds :
Noeud Coor_x
Coor_y
Coor_z
N01 500
866.025
0.
N02
0 0 0.
N03 866.025
500
0.
N04 366.025
1366.025
0.
Les mailles utilisées sont :
1 maille QUAD9
2 mailles TRIA7
3.3
Fonctionnalités testées
Commandes
Option
CALC_ELEM
SIGM_ELNO_COQUE
No2
No1
No4
No3
Code_Aster
®
Version
6.1
Titre :
SSNV152 - Traction élastique. Calcul des contraintes de Cauchy
Date
:
19/08/02
Auteur(s) :
P. MASSIN, G. BERTRAND
Clé
:
V6.04.152-A
Page :
6/6
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des éléments volumiques
HT-66/02/001/A
4
Résultats de la modélisation A
4.1 Valeurs
testées
Identification Référence
Aster
Différence
DX(No4) 8.66025
E+01
8.66025 E+01
4.66 E-05%
DY(No4) 50.0
50.0
0%
PK2-SIXX(M1) 21000.
21000.
2.04
E-08%
Cauchy-SIXX(M1) 23100.
23100.
2.14
E-08%
5
Synthèse des résultats
Les résultats trouvés sont en accord avec la solution analytique.