Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
1/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
Organisme(s) :
EDF/AMA
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
Document : V6.04.156
SSNV156 - Colonne sous chargement volumique.
Loi élastoplastique a gradient
Résumé :
Outre la loi élastoplastique à gradient qu'il ne teste que dans sa version à écrouissage linéaire, ce test a surtout
pour objet de valider l'algorithme d'intégration des lois de comportement à gradient de variables internes. En
effet, le problème proposé, la mise en traction d'une colonne sous forces volumiques, conduit à une solution
non homogène qui active les différentes composantes de l'algorithme (Newton, recherche linéaire, BFGS) et
pour laquelle on peut obtenir une expression analytique. Les résultats obtenus sont en accord avec celle-ci, et
ce avec une grande précision.
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
2/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
1
Problème de référence
1.1 Géométrie
La structure est un cylindre, au sens le plus général, de hauteur
mm
2
=
L
et dont la forme de la
section n'influe pas sur la solution. De ce fait, on adoptera une section carrée pour les modélisations
en déformations planes et en 3D (côté
mm
1
,
0
=
a
) ainsi qu'une section circulaire pour la
modélisation axisymétrique (rayon
mm
1
,
0
=
R
).
1.2
Propriétés du matériau
Le matériau suit une loi de comportement élastoplastique à gradient dont l'écrouissage est isotrope et
linéaire. Les caractéristiques du matériau, respectivement le module de Young
E
, le coefficient de
Poisson
NU
, la limite d'élasticité
SY
, la pente d'écrouissage
D_SIGM_EPSI
et la longueur
caractéristique
LONG_CARA
, sont égales à :
MPa
000
100
=
E
3
,
0
=
MPa
100
=
y
MPa
000
10
=
T
E
mm
707
982
,
0
=
b
L
1.3
Conditions aux limites et chargement
Les déplacements verticaux sont bloqués sur sa face supérieure, les déplacements horizontaux sont
bloqués sur les faces latérales et la face inférieure est libre de toute condition cinématique. Par
ailleurs, la structure est soumise à une force volumique verticale et dirigée vers le bas d'intensité
croissante
t
f
t
f
0
)
(
=
où
3
0
N/mm
1
=
f
et
t
est un paramètre de chargement (sans unité).
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
3/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
2
Solution de référence
2.1
Méthode de calcul de la solution de référence
Ce problème admet une solution analytique. Les grandeurs recherchées, à savoir les contraintes
,
les déformations
, les déformations plastiques
p
et la déformation plastique cumulée
p
, ne
dépendent que du niveau de chargement
t
et de la cote
z
de la section considérée, où
0
=
z
désigne la face inférieure (libre) du cylindre et
L
z
=
sa face supérieure (bloquée).
Du fait du blocage des faces latérales, les sections ne se déforment pas horizontalement, si bien que
les champs de déplacements et de déformations s'écrivent :
0
)
,
(
)
,
(
)
,
(
)
,
(
)
,
(
=
=
=
t
L
u
t
z
t
z
t
z
u
t
z
avec
z
z
z
e
e
e
u
éq 2.1-1
Quant au tenseur de contraintes, il est diagonal. Sa composante verticale
est fixée par l'équation
d'équilibre tandis que ses composantes horizontales, identiques dans les deux directions, dépendent
de la loi de comportement (effet de bridage) :
(
)
)
(
)
,
(
)
,
(
)
,
(
)
,
(
t
f
z
t
z
t
z
t
z
s
t
z
=
+
+
=
avec
z
z
y
y
x
x
e
e
e
e
e
e
éq
2.1-2
On peut remarquer que l'évolution du déviateur des contraintes est radiale :
(
)
s
eq
eq
-
=
+
+
-
=
et
z
z
y
y
x
x
e
e
e
e
e
e
3
2
3
1
D
éq 2.1-3
Ayant proposé un champ de déplacements cinématiquement admissible et un champ de contraintes
statiquement admissible, il ne reste plus qu'à montrer qu'ils sont liés par la loi de comportement. Bien
qu'il s'agisse d'un modèle non local, la loi d'écoulement conserve sa forme usuelle :
(
)
+
+
-
=
=
z
z
y
y
x
x
p
p
e
e
e
e
e
e
2
1
2
3
D
p
p
eq
&
&
éq
2.1-4
Il en va de même de la relation contrainte - déformation, où
et
µ
sont les coefficients de Lamé qui
se déduisent du module de Young et du coefficient de Poisson :
( )
(
)
p
Id
-
µ
+
=
2
tr
éq
2.1-5
En reportant [éq 2.1-1], [éq 2.1-2] et [éq 2.1-4] dans [éq 2.1-5], on en déduit l'expression des
contraintes et des déformations en fonction de la déformation plastique cumulée uniquement :
p
E
f
z
f
z
p
E
s
-
-
+
-
-
=
+
-
=
1
1
1
2
1
2
1
1
2
éq 2.1-6
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
4/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
La condition de cohérence permet alors de déterminer
p
. Dans la mesure où l'évolution du déviateur
des contraintes est radiale et monotone, on peut directement déterminer l'état actuel sans avoir à
intégrer le taux de déformation plastique. En effet, on peut distinguer deux zones dans la structure :
l'une,
b
z
0
, dans laquelle la déformation plastique est nulle et où le seuil de plasticité n'est pas
atteint, et l'autre,
L
z
b
, où la déformation plastique est non nulle et le seuil atteint. La frontière
b
entre ces deux zones est une nouvelle inconnue du problème. Ainsi :
=
=
0
F
0
0
L
z
b
p
b
z
éq
2.1-7
D'après [bib1], le seuil de plasticité a pour expression :
( )
13
4
,
F
2
b
T
T
y
eq
L
h
c
E
E
E
E
h
p
c
p
h
p
=
-
=
+
-
-
=
et
où
éq
2.1-8
En outre, toujours d'après [bib1], le champ de déformation plastique cumulée
p
est
1
C
et de dérivée
nulle au bord de la structure. Cela implique en particulier :
( )
( )
( )
0
0
0
=
=
=
L
p
b
p
b
p
éq
2.1-9
Compte-tenu de [éq 2.1-7] et de [éq 2.1-8],
p
vérifie une équation différentielle linéaire du second
ordre sur le domaine
L
z
b
:
(
)
y
f
z
z
p
h
E
z
p
c
-
-
-
=
+
-
+
1
2
1
)
(
1
2
)
(
éq
2.1-10
La donnée des 3 conditions aux limites [éq 2.1-9] permet alors de déterminer totalement
p
ainsi que
la position
b
de la frontière libre.
Pour alléger les expressions, on introduit les notations suivantes :
(
)
(
)
-
=
-
=
-
=
-
-
=
=
=
+
-
=
c
c
y
b
c
L
L
z
z
L
L
z
z
H
z
z
f
H
H
h
L
H
c
L
h
E
H
sinh
)
S(
cosh
)
C(
)
(
1
2
1
13
4
1
2
éq 2.1-11
Alors, après quelques calculs, on obtient l'expression suivante pour
p
:
-
=
-
=
+
+
=
c
c
L
B
b
b
b
L
A
z
B
z
A
z
z
p
)
C(
)
(
)
S(
où
)
S(
)
C(
)
(
)
(
éq 2.1-12
Quant à la frontière libre, elle est déterminée en fonction du niveau de chargement par l'équation
suivante (ou inversement, le niveau de chargement correspondant à une certaine position de la
frontière libre) :
(
)
)
C(
1
)
S(
)
S(
b
L
b
b
b
H
c
y
-
+
=
éq
2.1-13
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
5/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
2.2
Résultats de référence
On examine la déformation plastique cumulée
p
, la déformation totale
, la contrainte équivalente
eq
et la contrainte horizontale
s
au point
L
z
=
pour différents niveaux de chargement qui
correspondent à différentes positions de la frontière libre.
L
b
(
)
3
N/mm
f
p
(
)
MPa
eq
(
)
MPa
s
0, 75
104, 811 963
1, 165 975 E-4 1, 623 833 E-3 111, 456 702
98, 167 224
0, 50
146, 159 407
6, 125 415 E-4 2, 521 534 E-3 123, 286 355
169,032 459
0, 25
250, 078 993
1, 905 213 E-3 4, 804 152 E-3 149, 717 896
350, 440 090
0, 00
875, 079 453
9, 693 407 E-3 1, 854 027 E-2 307, 704 531 1 442, 454 356
2.3
Incertitudes sur la solution
Il s'agit d'une solution analytique.
2.4 Références
bibliographiques
[1]
Lorentz E., Andrieux S. : A variational formulation for nonlocal damage models. Int. J. Plas.,
15, pp. 119-138 (1999)
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
6/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
3 Modélisation
A
3.1
Caractéristiques de la modélisation
Il s'agit d'une modélisation 3D. La section du cylindre est carrée. L'écrouissage linéaire est caractérisé
dans un premier temps par une courbe d'écrouissage (
VMIS_ISOT_TRAC
) puis par sa limite
d'élasticité et son module (
VMIS_ISOT_LINE
), ce qui permet de tester les deux implantations.
3.2
Caractéristiques du maillage
Le maillage est réalisé par GMSH. Les mailles sont de plus petite taille dans la zone où l'on examine
les résultats (face supérieure). Au total, on dénombre 1749 TETRA10.
3.3 Fonctionnalités
testées
Commandes
DEFI_MATERIAU
NON_LOCAL
LONG_CARA
AFFE_MODELE
AFFE
MODELISATION = '3D_GRADIENT'
STAT_NON_LINE
LAGR_NON_LOCAL
4
Résultats de la modélisation A
4.1 Valeurs
testées
Les différentes valeurs sont testées à la cote
L
z
=
, comme présenté au [§2.2]. Il s'agit de la
déformation plastique cumulée (composante `
V1'
du
CHAM_NO
`
VARI_NOEU_ELGA
' et composante
`
VANL'
du
CHAM_NO
`
DEPL
'), de la déformation verticale (composante `
EPZZ
' du
CHAM_NO
`
EPSI_NOEU_DEPL
'), de la contrainte de von Mises (composante `
VMIS
' du
CHAM_NO
`
EQUI_NOEU_SIGM
') et enfin de la contrainte horizontale (composante `
SIXX
' du
CHAM_NO
`
SIEF_NOEU_DEPL
'). Pour rappel, les valeurs analytiques sont les suivante.
(
)
3
N/mm
f
p
(
)
MPa
eq
(
)
MPa
s
104, 811 963
1, 165 975 E-4
1, 623 833 E-3
111, 456 702
98, 167 224
146, 159 407
6, 125 415 E-4
2, 521 534 E-3
123, 286 355
169,032 459
250, 078 993
1, 905 213 E-3
4, 804 152 E-3
149, 717 896
350, 440 090
875, 079 453
9, 693 407 E-3
1, 854 027 E-2
307, 704 531
1 442, 454 356
Les résultats obtenus diffèrent extrêmement peu de la solution analytique (différence relative inférieure
à la précision par défaut de 0, 1 %).
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
7/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
5 Modélisation
B
5.1
Caractéristiques de la modélisation
Il s'agit d'une modélisation 3D. La section du cylindre est carrée. L'écrouissage linéaire est caractérisé
dans un premier temps par une courbe d'écrouissage (
VMIS_ISOT_TRAC
) puis par sa limite
d'élasticité et son module (
VMIS_ISOT_LINE
), ce qui permet de tester les deux implantations.
5.2
Caractéristiques du maillage
Le maillage est réalisé par GMSH. Les mailles sont de plus petite taille dans la zone où l'on examine
les résultats (face supérieure). Au total, on dénombre 169 TRIA6.
5.3 Fonctionnalités
testées
Commandes
DEFI_MATERIAU
NON_LOCAL
LONG_CARA
AFFE_MODELE
AFFE
MODELISATION = 'D_PLAN_GRADIENT'
STAT_NON_LINE
LAGR_NON_LOCAL
6
Résultats de la modélisation B
6.1 Valeurs
testées
Les différentes valeurs sont testées à la cote
L
z
=
, comme présenté au [§2.2]. Il s'agit de la
déformation plastique cumulée (composante `
V1'
du
CHAM_NO
`
VARI_NOEU_ELGA
' et composante
`
VANL
' du
CHAM_NO
`
DEPL
'), de la déformation verticale (composante `
EPYY
' du
CHAM_NO
`
EPSI_NOEU_DEPL
'), de la contrainte de von Mises (composante `
VMIS
' du
CHAM_NO
`
EQUI_NOEU_SIGM
') et enfin de la contrainte horizontale (composante `
SIXX
' du
CHAM_NO
`
SIEF_NOEU_DEPL
'). Pour rappel, les valeurs analytiques sont les suivante.
(
)
3
N/mm
f
p
(
)
MPa
eq
(
)
MPa
s
104, 811 963
1, 165 975 E-4
1, 623 833 E-3
111, 456 702
98, 167 224
146, 159 407
6, 125 415 E-4
2, 521 534 E-3
123, 286 355
169,032 459
250, 078 993
1, 905 213 E-3
4, 804 152 E-3
149, 717 896
350, 440 090
875, 079 453
9, 693 407 E-3
1, 854 027 E-2
307, 704 531
1 442, 454 356
Les résultats obtenus diffèrent extrêmement peu de la solution analytique (différence relative inférieure
à la précision par défaut de 0, 1 %).
Code_Aster
®
Version
6.0
Titre :
SSNV156 - Colonne sous chargement volumique
Date
:
14/10/02
Auteur(s) :
E. LORENTZ
Clé
:
V6.04.156-A
Page :
8/8
Manuel de Validation
Fascicule V6.04 : Statique non linéaire des structures volumiques
HT-66/02/001/A
7 Modélisation
C
7.1
Caractéristiques de la modélisation
Il s'agit d'une modélisation 2D axisymétrique. La section du cylindre est nécessairement circulaire.
L'écrouissage linéaire est caractérisé dans un premier temps par une courbe d'écrouissage
(
VMIS_ISOT_TRAC
) puis par sa limite d'élasticité et son module (
VMIS_ISOT_LINE
), ce qui permet
de tester les deux implantations.
7.2
Caractéristiques du maillage
Le maillage est réalisé par GMSH. A la différence des deux modélisations précédentes, il s'agit d'un
maillage réglé. Les mailles sont de plus petites tailles dans la zone où l'on examine les résultats (face
supérieure). Au total, on dénombre 50 QUAD8.
7.3 Fonctionnalités
testées
Commandes
DEFI_MATERIAU
NON_LOCAL
LONG_CARA
AFFE_MODELE
AFFE
MODELISATION = 'AXIS_GRADIENT'
STAT_NON_LINE
LAGR_NON_LOCAL
8
Résultats de la modélisation C
8.1 Valeurs
testées
Les différentes valeurs sont testées à la cote
L
z
=
, comme présenté au [§2.2]. Il s'agit de la
déformation plastique cumulée (composante `
V1'
du
CHAM_NO
`
VARI_NOEU_ELGA
' et composante
`
VANL
' du
CHAM_NO
`
DEPL
'), de la déformation verticale (composante `
EPYY
' du
CHAM_NO
`
EPSI_NOEU_DEPL
'), de la contrainte de von Mises (composante `
VMIS
' du
CHAM_NO
`
EQUI_NOEU_SIGM
') et enfin de la contrainte horizontale (composante `
SIXX
' du
CHAM_NO
`
SIEF_NOEU_DEPL
'). Pour rappel, les valeurs analytiques sont les suivante.
(
)
3
N/mm
f
p
(
)
MPa
eq
(
)
MPa
s
104, 811 963
1, 165 975 E-4
1, 623 833 E-3
111, 456 702
98, 167 224
146, 159 407
6, 125 415 E-4
2, 521 534 E-3
123, 286 355
169,032 459
250, 078 993
1, 905 213 E-3
4, 804 152 E-3
149, 717 896
350, 440 090
875, 079 453
9, 693 407 E-3
1, 854 027 E-2
307, 704 531
1 442, 454 356
Les résultats obtenus diffèrent extrêmement peu de la solution analytique (différence relative inférieure
à la précision par défaut de 0, 1 %).
9
Synthèse des résultats
Du fait du caractère positif de l'écrouissage, le problème n'exhibe pas d'instabilités, ce que confirme la
bonne convergence des calculs. La validation porte ainsi sur la loi de comportement proprement dite et
sur l'algorithme d'intégration des lois non locales, dont les différentes composantes sont activées
(itérations primales et duales observées). La remarquable concordance entre les valeurs de référence
et celles calculées démontre les capacités de l'algorithme lorsqu'une précision fine est requise (critères
de convergence sévères) tandis que la taille des problèmes traités, particulièrement en 3D, semble
prouver sa robustesse.