Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
1/10
Organisme(s) : EDF/IMA/MMN
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
Document V6.04.121
SSNV121 - Rotation et traction hyper-élastique
d'un barreau
Résumé :
Ce test de mécanique quasi-statique consiste à faire tourner de 90° un barreau parallélépipédique, à le
soumettre à une traction importante pour finalement le laisser revenir dans un état déchargé. On valide ainsi la
cinématique des grandes déformations hyper-élastiques (commande STAT_NON_LINE [U4.32.01], mot-clé
COMP_ELAS), et donc en particulier les grandes rotations, pour une relation de comportement élastique linéaire.
Le barreau est modélisé par un élément volumique (HEXA8, modélisation A) ou bien plan (QUAD4, hypothèse
de déformations planes, modélisation B).
Les résultats obtenus par le Code_Aster ne diffèrent pas de la solution théorique.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
2/10
1
Problème de référence
1.1 Géométrie
y
1 000 (mm)
1
2
3
4
x
1 000 (mm)
1.2
Propriétés de matériaux
Comportement hyper-élastique de St Venant - Kirchhoff :
E
E
E = 200 000.MPa
S
= (
tr
+
1 + )(1-
2 )
(E) 1
E
1 +
= 0 3
.
1.3
Conditions aux limites et chargements
Le chargement est appliqué en deux temps : tout d'abord, une rotation d'ensemble de la structure,
suivie par une traction dans la nouvelle configuration :
Rotation d'ensemble (0 < t < 1 s)
Traction (1 s < t < 2 s)
T
2'
4'
1
2
2'
4'
3
4
1'
1'
3
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
3/10
2
Solution de référence
2.1
Méthode de calcul utilisée pour la solution de référence
Il s'agit d'un problème plan. On peut chercher la solution sous la forme d'une rotation rigide et d'un
allongement d'un facteur dans la direction y .
- X - Y
U( X ,Y, Z) =
(1+ )X -Y
0
Le gradient de la transformation et la déformation de Green-Lagrange sont alors :
0
-1
0
e 0
0
( + 2)
F
=
1 +
0
0 E
= 0 0
0
avec e =
2
0
0
1
0 0 0
La relation de comportement conduit alors à un tenseur de contraintes lagrangiennes diagonal :
(1-) E
S
=
xx
(1+)(1- 2) e
E
S
= S
=
yy
zz
(1+)(1- 2) e
La condition aux limites de l'équation d'équilibre nous permet alors de déterminer la valeur de
l'allongement :
(1-)
+1 +
= (
2
FS)
= (1+ )
E
(
)(
)
T
S
(
=
1 + )(1- 2 )
T
yx
xx
2
La contrainte de Cauchy est donnée par :
S yy
1
xx
=
zz
=
=
F S F T
1 +
Det F
yy
= (1+ )
S xx
Enfin, la force exercée sur les faces :
· [2,4] : F
=
=
y
yy
[
S 2,4]
yy So[2,4]
· [4,3] : F
=
= 1+
x
xx
[
S 4, ]3
xx (
) So[4, ]3
· [1,2,3,4] : F
=
= 1+
z
zz
[
S 1,2,3,4]
zz (
) So[1,2,3,4]
où S [
o ]
représentent les surfaces initiales des faces.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
4/10
2.2
Résultats de référence
On adopte comme résultats de référence les déplacements, la contrainte de Cauchy et la force
exercée sur les faces [2,4] et [4,3].
Au temps t = 2 s :
On cherche T tel que l'allongement
= 01.
==> T = 31 096.154 MPa.
La contrainte de Cauchy est alors :
=
=
xx
zz
11 013.986 MPa
=
yy
31 096.154 MPa
Les forces exercées sont :
F
= 12 115 385
.
×
x
So[4, ] N
3
F
= 31 096.154 ×
y
So[2,4] N
F
= 12 115.385 ×
z
So[ ,12, ,34] N
Au temps t = 3 s :
Le barreau est revenu dans son état initial :
= 0
= 0
F =
0
2.3
Incertitude sur la solution
Solution analytique.
2.4 Références
bibliographiques
[1]
Eric LORENTZ "Une relation de comportement hyperélastique non linéaire" Note interne
EDF/DER HI-74/95/011/0
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
5/10
Page laissée intentionnellement blanche.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
6/10
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Modélisation volumique :
1 maille HEXA 8
1 maille QUAD4
z
5
6
y
7
8
1
2
1 000 (mm)
3
4
x
Conditions aux limites :
R t
( )
(3,7) : DX = 0
DY = 0
1
(1,5) : DX = - 1 000 R(t)
DY = - 1 000 R(t)
(2,6) : DX = - 2 000 R(t)
(4,8) : DX = - 1 000 R(t)
0
1
2
3
t
Chargement : Traction sur la face [2,4,8,6]
maille [2,4,8,6] (QUAD4) : FY = 31 096.154 f(t) MPa
f t
( )
1
0
1
2
3
t
3.2
Caractéristiques du maillage
Nombre de noeuds : 8
Nombre de mailles : 2
1 HEXA8
1 QUAD4
3.3 Fonctionnalités
testées
STAT_NON_LINE
COMP_ELAS
DEFORMATION : 'GREEN'
CALC_NO
OPTION : 'FORC_NODA'
GEOMETRIE : 'DEFORMEE'
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
7/10
4
Résultats de la modélisation A
4.1 Valeurs
testées
Identification
Référence
Aster
% différence
t = 2 Déplacement DX (NO2)
100
100
0
t = 2 Déplacement DY (NO4)
1100
1100
0
t = 2 Contraintes SIGXX (PG1)
11013.986
11013.986
0
t = 2 Contraintes SIGYY (PG1)
31096.154
31096.154
0
t = 2 Contraintes SIGZZ (PG1)
11013.986
11013.986
0
t = 2 Contraintes SIGXY (PG1)
0
109
/
t = 2 Contraintes SIGXZ (PG1)
0
1010
/
t = 2 Contraintes SIGYZ (PG1)
0
1010
/
t = 3 Déplacement DX10 (NO2)
0
1011
/
t = 3 Déplacement DY (NO4)
0
1012
/
t = 3 Contraintes SIGXX (PG1)
0
109
/
t = 3 Contraintes SIGYY (PG1)
0
1010
/
t = 3 Contraintes SIGZZ (PG1)
0
109
/
t = 3 Contraintes SIGXY (PG1)
0
1011
/
t = 3 Contraintes SIGXZ (PG1)
0
1010
/
t = 3 Contraintes SIGYZ (PG1)
0
1011
/
t = 2 Force nodale DX (NO8)
3.0289 109
3.0288 109
0.002%
t = 2 Force nodale DY (NO8)
7.774 109
7.774 109
0
t = 2 Force nodale DZ (NO8)
3.0289 109
3.0288 109
0.002%
4.2 Remarques
Calcul de la force nodale :
La force appliquée F sur une face décrite par une maille linéaire se répartit par :
1/4
1/4
F
= 1
noeud
F
4
1/4
1/4
4.3 Paramètres
d'exécution
Version : 3.03.30
Machine : CRAY C90
Encombrement mémoire :
8 MW
Temps CPU User :
62.3 secondes
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
8/10
5 Modélisation
B
5.1
Caractéristiques de la modélisation
Modélisation 2D déformations planes
y
1
2
3
4
x
Conditions aux limites :
R t
( )
3 :
DX = 0
DY = 0
1
1 :
DX = - 1 000 R(t)
DY = - 1 000 R(t)
2 :
DX = - 2 000 R(t)
4 :
DX = - 1 000 R(t)
0
1
2
3
t
Chargement :
Traction sur la face [2,4]
maille [2,4] : FY = 31 096.154 f(t) MPa
f t
( )
1
0
1
2
3
t
5.2
Caractéristiques du maillage
Nombre de noeuds : 4
Nombre de mailles : 2
1 QUAD4
1 SEG2
5.3 Fonctionnalités
testées
STAT_NON_LINE
COMP_ELAS
DEFORMATION : 'GREEN'
CALC_NO
OPTION : 'FORC_NODA'
GEOMETRIE : 'DEFORMEE'
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
9/10
6
Résultats de la modélisation B
6.1 Valeurs
testées
Identification
Référence
Aster
% différence
t = 2 Déplacement DX (NO2)
100
100
0
t = 2 Déplacement DY (NO4)
1100
1100
0
t = 2 Contraintes SIGXX (PG1)
11013.986
11013.986
0
t = 2 Contraintes SIGYY (PG1)
31096.154
31096.154
0
t = 2 Contraintes SIGZZ (PG1)
11013.986
11013.986
0
t = 2 Contraintes SIGXY (PG1)
0
10-10
/
t = 3 Déplacement DX (NO2)
0
1012
/
t = 3 Déplacement DY (NO4)
0
1012
/
t = 3 Contraintes SIGXX (PG1)
0
1010
/
t = 3 Contraintes SIGYY (PG1)
0
1010
/
t = 3 Contraintes SIGZZ (PG1)
0
1010
/
t = 3 Contraintes SIGXY (PG1)
0
1010
/
t = 2 Force nodale DX (NO4)
6.0577 106
6.0577 106
0
t = 2 Force nodale DY (NO4)
15.5481 106
15.5481 106
0
6.2 Remarques
Calcul de la force nodale :
La force appliquée F sur une face décrite par une maille linéaire se répartit par :
1/2
1/2
F
= 1
noeud
F
2
6.3 Paramètres
d'exécution
Version : NEW 3.03.30
Machine : CRAY C90
Encombrement mémoire :
8 MW
Temps CPU User :
48.3 secondes
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A
Code_Aster ®
Version
3
Titre :
SSNV121 Rotation et traction hyper-élastique d'un barreau
Date :
23/07/99
Auteur(s) :
E. LORENTZ
Clé :
V6.04.121-A Page :
10/10
7
Synthèse des résultats
Il apparaît à l'issue de ce test que la solution numérique coïncide remarquablement avec la solution
analytique. On remarquera cependant que la forte non linéarité due aux grandes rotations nécessite
une discrétisation en temps relativement fine, sans être pénalisante sur la précision puisque,
contrairement à une loi de comportement incrémentale, les erreurs ne se cumulent pas d'un pas de
temps sur l'autre.
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HI-75/96/044 - Ind A