Code - Aster: Notice D'utilisation Des Éléments Plaques, Coques, Coques Volumiques SHB, Grilles Et Membranes

Code_Aster
Version
default
Titre : Notice dutilisation des lments plaques, coques,[...] Date : 15/05/2013 Page : 1/55
Responsable : Thomas DE SOZA Cl : U2.02.01 Rvision : 11028
Notice dutilisation des lments plaques, coques,
coques volumiques SHB, grilles et membranes
Rsum :
Ce document est une notice dutilisation pour les modlisations plaques, coques, coques volumiques SHB,
grilles et membranes.
Les lments de coques et de plaques interviennent dans la modlisation numrique des structures minces
surface moyenne, plane (modlisation plaque) ou courbe (modlisation coques). Les modlisations de grilles
interviennent pour la modlisation numrique des armatures et des cbles de prcontrainte dans les structures
en bton arm. La modlisation de membrane permet de modliser des structures minces dont la rigidit de
flexion est ngligeable.
Ils sont utilisables en mcanique linaire ou non linaire et en thermique.
Manuel d'utilisation Fascicule u2.02 : lments de structure
Document diffus sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)
Code_Aster
Version
default
Table des matires
1 Introduction .......................................................................................................................................... 4
2 Mcanique ........................................................................................................................................... 6
2.1 Capacits de modlisation ............................................................................................................ 6
2.1.1 Rappel des formulations ...................................................................................................... 6
2.1.2 Formulation des lments plaque, coques et coques volumiques ..................................... 9
2.1.3 Formulation des lments coques volumiques SHB ........................................................ 11
2.1.4 Comparaison entre les lments ....................................................................................... 12
2.2 Commandes utiliser ................................................................................................................. 15
2.2.1 Discrtisation spatiale et affectation dune modlisation : oprateur AFFE_MODELE .... 15
2.2.2 Caractristiques lmentaires : AFFE_CARA_ELEM ...................................................... 17
2.2.3 Matriaux : DEFI_MATERIAU .......................................................................................... 20
2.2.4 Chargements et conditions limites : AFFE_CHAR_MECA et AFFE_CHAR_MECA_F .... 21
2.3 Rsolution ................................................................................................................................... 26
2.3.1 Calculs linaires : MECA_STATIQUE et autres oprateurs linaires ............................... 26
2.3.2 Calculs non linaires : STAT_NON_LINE et DYNA_NON_LINE ..................................... 27
2.4 Calculs supplmentaires et post-traitements .............................................................................. 30
2.4.1 Calculs de matrices lmentaires : oprateur CALC_MATR_ELEM ................................ 30
2.4.2 Calculs par lments : oprateurs CALC_CHAMP et POST_CHAMP ............................ 30
2.4.3 Calculs aux nuds : oprateur CALC_CHAMP ............................................................... 32
2.4.4 Calculs de quantits sur tout ou partie de la structure : oprateur POST_ELEM ............. 32
2.4.5 Valeurs de composantes de champs de grandeurs : oprateur POST_RELEVE_T ........ 33
2.4.6 Impression des rsultats : oprateur IMPR_RESU ........................................................... 34
2.5 Exemples .................................................................................................................................... 34
2.5.1 Analyse statique linaire ................................................................................................... 34
2.5.2 Analyse modale en dynamique ......................................................................................... 37
2.5.3 Analyse statique non linaire matriau .............................................................................. 38
2.5.4 Analyse statique non linaire gomtrique ........................................................................ 39
2.5.5 Analyse en flambement dEuler ......................................................................................... 42
2.5.6 Raccords coques et autres lments mcaniques ............................................................ 43
3 Thermique ......................................................................................................................................... 45
3.1 Dfinition du problme ................................................................................................................ 45
3.1.1 Discrtisation spatiale et affectation dune modlisation : oprateur AFFE_MODELE .... 45
3.1.2 Caractristiques lmentaires : AFFE_CARA_ELEM ...................................................... 46
3.1.3 Matriaux : DEFI_MATERIAU .......................................................................................... 46
3.1.4 Chargements et conditions limites : AFFE_CHAR_THER et AFFE_CHAR_THER_F ..... 47
3.2 Rsolution ................................................................................................................................... 48
3.2.1 Calculs transitoires : oprateur THER_LINEAIRE ............................................................ 48
3.3 Calculs supplmentaires et post-traitements .............................................................................. 48
Code_Aster
Version
default
3.3.1 Calculs en post-traitement ................................................................................................. 48
3.4 Exemples ................................................................................................................................... 48
4 Chanage thermomcanique ............................................................................................................. 49
4.1 Formalisme ................................................................................................................................. 49
5 Conclusion et conseils dutilisation .................................................................................................... 51
6 Bibliographie ...................................................................................................................................... 53
Code_Aster
Version
default
1 Introduction
Les lments de coques et de plaques sont particulirement utiliss pour modliser des structures
minces o les rapports entre les dimensions (paisseur/longueur caractristique) sont trs infrieurs
1/ 10 (coques minces) ou de lordre de 1/ 10 (coques paisses).
Ces modlisations sont utilisables en mcanique linaire et non linaire, sous des hypothses de
petites dformations et petits dplacements ou bien des hypothses de grands dplacements et de
grandes rotations, suivant les modlisations. Une modlisation de coque mince est galement
disponible en thermique linaire transitoire.
Trois catgories dlments de structures minces sont dcrites dans ce document :
Les lments de plaques, qui sont plans, donc la courbure de la structure reprsenter nest
pas idalement prise en compte et il est ncessaire dutiliser un grand nombre dlments de
faon approcher correctement la gomtrie de la structure (aspect facettes).
Les lments de coques, qui sont courbes, donc la gomtrie de la structure est mieux
approche.
Les lments de coques volumiques ou tridimensionnelles SHB , qui sont en fait des
lments 3D enrichis, apprhendant au mieux les phnomnes pour de faibles paisseurs.
Concernant la nomenclature et les documentations de rfrence associes chacune des
modlisations :
Les lments de plaques plans triangle et quadrangle sont regroups sous les modlisations,
(documentation de rfrence [R3.07.03]) :
DKT : maille TRIA3 lment DKT, maille QUAD4 lments DKQ (linaire gomtrique) ;
DKTG : maille TRIA3 lment DKT, maille QUAD4 lments DKQ (linaire ou non linaire
gomtrique) ;
DST : maille TRIA3 lment DST, maille QUAD4 lment DSQ (linaire gomtrique) ;
Q4G : maille QUAD4 lment Q4G (linaire gomtrique) .
Q4GG : maille TRIA3 lment T3G, maille QUAD4 lments Q4G (linaire gomtrique) ;
Les lments de coques courbes issus de modles 3D avec une cinmatique de coque sont
regroups sous les modlisations :
COQUE_3D : maille TRIA7 et QUAD9, structure 3D gomtrie quelconque ([R3.07.04] en
linaire gomtrique, [R3.07.05] en non linaire gomtrique et [R3.03.07] avec des
pressions suiveuses) ;
COQUE_AXIS : maille SEG3, coques symtrie de rvolution autour de laxe 0Y
([R3.07.02] en linaire gomtrique) ;
COQUE_C_PLAN ou COQUE_D_PLAN : maille SEG3, coques gomtrie invariante le long de
laxe OZ ([R3.07.02] en linaire gomtrique) ;
Les lments de coques volumiques avec une cinmatique isoparamtrique comme les
lments 3D standards ( en linaire et non linaire gomtrique avec des pressions
suiveuses) sont regroups sous la modlisation (documentation de rfrence [R3.07.07a]) :
SHB8 : maille PENTA6 lment SHB6 (formulation linaire) ;
maille HEXA8 lment SHB8 (formulation linaire) ;
maille PENTA15 lment SHB15 ( formulation quadratique );
maille HEXA20 lment SHB20 (formulation quadratique).
lments de membrane classiques :
MEMBRANE : maille TRIA3, QUAD4, TRIA6 et QUAD8 (linaire gomtrique)
Ce sont des lments de membrane avec une simple rigidit membranaire (pas de degr de
libert de rotation). Il n'y a pas d'excentrement.
lments de plaques spcifiques pour reprsenter les nappes de ferraillage :
GRILLE_EXCENTRE : maille TRIA3 ou QUAD4 (linaire gomtrique)
La modlisation GRILLE_EXCENTRE correspond des lments de plaques DKT
orthotropes 3 ou 4 nuds excentrs par rapport au feuillet moyen en bton, (une seule
direction d'armatures). Ce dernier est modlis par des lments de plaques DKT ou DST
Code_Aster
Version
default
3 nuds. La structure en bton arm est alors reprsente par la superposition des
modlisations GRILLE_EXCENTRE et de celle utilise pour le bton (DKT ou DST).
GRILLE_MEMBRANE : mailles TRIA3, QUAD4, TRIA6, QUAD8 (doc [R3.08.07]): ce sont des
lments d'armature (une seule direction d'armatures) ne travaillant qu'en membrane
(pas de degr de libert de rotation). Il n'y a pas d'excentrement.
Code_Aster
Version
default
2 Mcanique

2.1 Capacits de modlisation
2.1.1 Rappel des formulations
2.1.1.1 Gomtrie des lments plaque et coques
Pour les lments plaques et coques on dfinit une surface de rfrence, ou surface moyenne, plane
(plan xy par exemple) ou courbe ( x et y dfinissent un ensemble de coordonnes curvilignes) et
une paisseur h( x , y) . Cette paisseur doit tre petite par rapport aux autres dimensions de la
structure modliser. La figure ci-dessous illustre ces diffrentes configurations.

Figure 2.1.1.1-a : Hypothse en thorie des plaques et des coques
On attache la surface moyenne
o
un repre local Oxyz diffrent du repre global OXYZ . La
position des points de la plaque ou de la coque est donne par les coordonnes curvilignes ((
1,
(
2
)
de la surface moyenne et llvation (
3
par rapport cette surface. Pour les plaques le systme de
coordonnes curvilignes est un systme de coordonnes cartsiennes locales.
Code_Aster
Version
default

Figure 2.1.1.1-b : Dfinition dune surface moyenne
Pour reprsenter des coques symtrie de rvolution autour dun axe (COQUE_AXIS) ou des coques
gomtrie invariante par translation (COQUE_C_PLAN en contraintes planes ou COQUE_D_PLAN en
dformations planes), la connaissance de la section de rvolution ou de la trace de la surface
moyenne est suffisante, comme la figure [Figure 2.1.2.1-a] nous le montre. Ces coques sappuient sur
un maillage linique et en un point m de la surface moyenne on dfinit un repre local (n , t , e
z
)
par :
t =
Om
, s
Om
, s
, nt=e
z
Lorsque lon souhaite modliser un solide de forme quelconque (non plan), on peut utiliser des
lments de coques pour rendre compte de la courbure, ou bien des lments de plaques. Dans ce
dernier cas, la gomtrie est approxime par un rseau de facettes.

Figure 2.1.1.1-c : Modlisation dun solide 3D quelconque par des lments de plaques ou coques
Pour la dfinition du repre intrinsque, on renvoie [U4.74.01] paragraphe 3.4.
2.1.1.2 Gomtrie des lments coques volumiques SHB
Code_Aster
Version
default
Les lments coques volumiques SHB sont des lments de gomtrie tridimensionnelle continue
dans lesquels une direction privilgie, appele paisseur, a t choisie. Cette paisseur est, en
gnral, petite par rapport aux autres dimensions de la structure modliser. Cette direction de
lpaisseur est dfinie par la faon de mailler. En gnral, on a besoin de deux surfaces pour gnrer
un volume tridimensionnel. La direction qui est perpendiculaire ces deux surfaces est la direction de
lpaisseur. Ces lments peuvent ainsi tre utiliss pour modliser les structures minces et pour
prendre en compte les phnomnes qui se dveloppent dans lpaisseur dans le cadre de la
mcanique des milieux continus tridimensionnelle. En approche linaire, on a llment prismatique
six nuds SHB6 et llment hexadrique huit nuds SHB8. En approche quadratique, on a
llment prismatique quinze nuds SHB15 et llment hexadrique vingt nuds SHB20. La
figure ci-dessous illustre ces diffrentes modlisations :

Solide 3D
X
Y
Z
h
L
b

Solide 3D
X
Y
Z
h
L
b

Figure 2.1.1.2-a : Gomtries des lments de rfrence et points dintgration

SHB8

9
10
11
12
13
14
15
16
17
18
19
20
1
2
3
4
5
11
12
13
14
15
6
7
8
9
10
16
17
18
19
20
SHB20

SHB6

11
15
12
9
10
14
13
1
2
3
4
5
11
12
13
14
15
6
7
8
9
10
SHB15
Code_Aster
Version
default
2.1.2 Formulation des lments plaque, coques et coques volumiques
2.1.2.1 Formulation en linaire gomtrique

Dans cette formulation, on suppose que les dplacements sont petits, on peut donc superposer la
gomtrie initiale et la gomtrie dforme. Ces lments (hors SHB) reposent sur la thorie des
coques selon laquelle :
les sections droites qui sont les sections perpendiculaires la surface de rfrence restent
droites; les points matriels situs sur une normale la surface moyenne non dforme
restent sur une droite dans la configuration dforme. Il rsulte de cette approche que les
champs de dplacement varient linairement dans lpaisseur de la plaque ou de la
coque. Si lon dsigne par u , v et w les dplacements dun point q( x , y , z) suivant x ,
y et z , on a ainsi :

u( x , y , z)
v( x , y , z )
w( x , y , z )
u( x , y)
v( x , y)
w( x , y )
+z
x
( x , y )
y
( x , y )
0

Le tenseur de dformation associ scrit alors :
e ( x , y , z )=e( x , y )+( x , y )+z. X( x , y )
.
Le premier terme e comprend les dformations de membrane (pour un lment de plaque
ce sont les dformations dans le plan de llment), le second
celles de cisaillement
transversal, et le troisime
z .X
les dformations de flexion, associes au tenseur de
courbure
X
. Pour les plaques ou les coques paisses les cisaillements transverses
sont
pris en compte suivant la formulation propose par Reissner, Hencky, Boll, Mindlin. Cette
formulation englobe lapproche sans cisaillement transverse, o le tenseur
est nul,
dveloppe par Kirchhoff pour les plaques ou les coques minces selon laquelle les points
matriels situs sur une normale n la surface moyenne non dforme restent sur la
normale la surface dforme.
La contrainte transversale c
zz
est nulle car considre comme ngligeable par rapport
aux autres composantes du tenseur des contraintes (hypothse des contraintes planes).
On ne dcrit pas la variation de lpaisseur ni celle de la dformation transversale e
zz
que
lon peut calculer en utilisant lhypothse prcdente de contraintes planes.
La prise en compte du cisaillement transverse dpend de facteurs de correction dtermins a
priori par des quivalences nergtiques avec des modles 3D, de faon ce que la rigidit
en cisaillement transverse du modle de plaque soit la plus proche possible de celle dfinie
par la thorie de llasticit tridimensionnelle. Pour les plaques homognes, le facteur de
correction de cisaillement transverse bas sur cette mthode est k=5/ 6 .
Remarque :
La dtermination des facteurs de correction repose pour Mindlin sur des quivalences de
frquence propre associe au mode de vibration par cisaillement transverse. On obtient
alors
k=n
2
/12
, valeur trs proche de 5/ 6 .
Ces lments font intervenir localement :
Cinq variables cinmatiques pour les lments plaques et coques quelconques ; les
dplacements de membrane u et v dans le plan de rfrence z=0 , le dplacement
transversal w et les rotations
x
et
y
de la normale la surface moyenne dans les plans
yz et xz respectivement.
Code_Aster
Version
default

Trois variables cinmatiques pour les lments liniques ; les dplacements u et v dans le
plan de rfrence z=0 et la rotation
n
de la normale la surface moyenne dans le plan
xy .

Figure 2.1.2.1-a : Variables cinmatiques pour les diffrents lments de plaques et de coques
Trois efforts rsultants de membrane nots Nxx , Nyy , Nxy et trois moments de flexion
nots Mxx , Myy , Mxy quel que soit llment de plaque ou de coque ; deux efforts
tranchants nots Vx et Vy dans le cas des lments de plaques et de coques
quelconques.

Figure 2.1.2.1-b : Efforts rsultants pour un lment de plaque ou de coque
Code_Aster
Version
default
2.1.2.2 Formulation en non-linaire gomtrique, Flambement dEuler
Dans la formulation en non-linaire gomtrique, on est en prsence de grands dplacements et de
grandes rotations, on ne peut pas superposer la gomtrie initiale et la gomtrie dforme.
La formulation, dcrite dans la documentation de rfrence [R3.07.05], est base sur une approche de
milieu continu 3D, dgnre par l'introduction de la cinmatique de coque de type
Hencky-Mindlin-Naghdi en contraintes planes dans la formulation faible de l'quilibre. La mesure des
dformations retenues est celle de Green-Lagrange, nergiquement conjugues aux contraintes de
Piola-Kirchhoff de deuxime espce. La formulation de l'quilibre est donc une formulation
Lagrangienne totale. Le cisaillement transverse est trait de la mme manire que dans le cas linaire
[R3.07.04].
L'lment retenu en non-linaire est un lment de coque volumique (COQUE_3D) de surface
moyenne courbe tel que prsent au paragraphe prcdent, dont les mailles supports sont des QUAD9
et des TRIA7.
Il est possible d'appliquer sur ces lments des pressions suiveuses, dont la formulation est dcrite
dans le document de rfrence [R3.03.07]. Ce chargement la particularit de suivre la gomtrie de
la structure au cours de sa dformation (par exemple, la pression hydrostatique reste toujours
perpendiculaire la gomtrie dforme).
Le flambement linaire appel aussi flambement d'Euler, dcrit dans la documentation de rfrence
[R3.07.05], se prsente comme un cas particulier du problme non-linaire gomtrique. Il est bas
sur une dpendance linaire des champs de dplacements, de dformations et de contraintes par
rapport au niveau de charge.
L'lment retenu en flambement linaire est l'lment de coque volumique (COQUE_3D) de surface
moyenne courbe tel que prsent au paragraphe prcdent, dont les mailles supports sont des QUAD9
et des TRIA7.
2.1.3 Formulation des lments coques volumiques SHB
Les lments coques volumiques SHB (SHB6, SHB8, SHB15 et SHB20) sont iso-paramtriques.
Leurs formulations sont tablies en sappuyant sur les mailles 3D volumiques standards
(successivement penta6, hexa8, penta15 et hexa20) avec seulement 3 degrs de libert de
dplacement chaque nud de llment.
Les formulations des ces lments sont dcrites dans les documentation de rfrence [R3.07.07] et
[R3.07.08]. Les quatre lments SHB reposent sur des formulations sous-intgres. Lide de base
consiste tout dabord sassurer quil y a suffisamment de points de Gauss dans lpaisseur pour
reprsenter correctement le phnomne de flexion, puis calculer des rigidits de stabilisation de
manire adaptative suivant ltat plastique de llment. Cela reprsente une amlioration certaine par
rapport aux formulations classiques pour les forces de stabilisation, car ces dernires reposent sur une
stabilisation lastique qui devient trop rigide lorsque les effets de la plasticit dominent la rponse de
la structure.
Le principe des formulations des lments SHB est assez identique, mais il existe quelques
diffrences entre les lments :
llment SHB8 est sous-intgr par cinq points de Gauss disposs dans la direction de
lpaisseur et stabilis par lintroduction dun champ de dformation postule assumed
strain ;
llment SHB6 est aussi sous-intgr par cinq points de Gauss disposs dans la direction de
lpaisseur. Cet lment ne prsentant pas de mode de hourglass, la stabilisation nest
donc pas ncessaire. Pour amliorer sa rapidit de convergence, on introduit des
modifications de type assumed strain sur loprateur gradient discrtis de llment ;
les lments quadratiques SHB15 et SHB20 sont sous-intgrs respectivement par quinze et
vingt points de Gauss dans la direction de lpaisseur. Ils nont pas besoin de stabilisation
ou modification de type assumed strain .
Code_Aster
Version
default

En non-linaire gomtrique, les lments SHB sont traits dans le cas de grands dplacements,
rotations faibles et petites dformations. On adopte pour cela une formulation lagrangienne mise
jour. Dans le cas des petits dplacements, on confond gomtrie en dbut et fin de pas, contrainte de
Cauchy et de Piola Kirchoff 2, de plus on utilise lexpression linaire des dformations.
Une particularit des lments SHB est que la matrice de comportement lastique
C
est choisie sous
la forme suivante (crite dans les repres locaux):
C=
|
\+2j

\ 0 0 0 0
\ \+2j 0 0 0 0
0 0 E 0 0 0
0 0 0 j 0 0
0 0 0 0 j 0
0 0 0 0 0 j

o E est le module dYoung,
+
le coefficient de Poisson, j=
E
2(1++)
le module de cisaillement et
\=
E +
1+
le coefficient de Lam modifi. Cette loi est spcifique aux lments SHB . Elle ressemble
celle que lon aurait dans le cas de lhypothse des contraintes planes, mis part le terme (3,3) . On
peut noter que ce choix entrane un comportement anisotropique artificiel. Ce choix permet de
satisfaire tous les tests sans introduire de blocage.
En non-linaire matriau, une mthode de construction particulire de la matrice tangente
C
T
est
utilise. Elle consiste supposer dabord que llment est en tat de contrainte plane dans le repre
local de chaque point de Gauss et les dformations hors plan sont lastiques. Cela entrane alors
immdiatement que les dformations totales hors plan sont gales aux dformations lastiques.
Appelons
C
CPT
la matrice tangente en contraintes planes. La matrice de comportement tangente
pour le comportement choisi est vidente et scrit :
C
T
=
|
C
xxxx
CPT
C
xxyy
CPT
0 C
xxxy
CPT
0 0
C
xyyx
CPT
C
yyyy
CPT
0 C
yyxy
CPT
0 0
0 0 E 0 0 0
C
xyxx
CPT
C
xyyy
CPT
0 C
xyxy
CPT
0 0
0 0 0 0 j 0
0 0 0 0 0 j

Ensuite, les contraintes hors plan sont calcules de faon lastique. Cette mthode permet ainsi de
connecter les lments SHB toutes les lois de comportement disponibles dans Code Aster .
2.1.4 Comparaison entre les lments
2.1.4.1 Les diffrences entre les lments plaques, coques et coques volumiques SHB
Les lments de coque sont des lments courbes alors que les lments de plaques sont plans. La
variation de mtrique de la gomtrie (cest--dire son rayon de courbure) en fonction de son
paisseur est prise en compte pour les lments de coques mais pas pour les lments de plaques.
Cette variation de mtrique implique un couplage entre les effets de membrane et de flexion pour des
structures non planes qui ne peut pas tre observ avec des lments de plaque plan pour un
matriau homogne (voir [bib1]).
Code_Aster
Version
default
Le choix des fonctions de forme pour la discrtisation de ces lments est diffrent car les lments
de coques courbes ont un nombre de degrs de libert plus important. Ainsi, les lments de plaques
sont des lments linaires en membrane alors que les lments de coques sont quadratiques.
Les lments coques volumiques SHB sappuient sur des mailles 3D. Il ny a pas ainsi de problmes
de raccord avec les lments solides 3D. Donc, les lments SHB sont utilisables pour modliser des
structures minces qui coexistent avec les lments solides 3D. En plus, il est plus pratique dutiliser
les lments SHB pour modliser des structures dont lpaisseur varie continument. La discrtisation
des lments SHB est soit linaire pour SHB6 et SHB8, soit quadratique pour SHB15 et SHB20 .
2.1.4.2 Les diffrences entre les lments plaques
On distingue les lments avec cisaillement transversal (DST, DSQ et Q4G) des lments sans
cisaillement transversal (DKT et DKQ). Les lments DST et DKT ont des mailles support triangulaires
3 nuds ( 35=15 degrs de libert) et les lments DKQ, DSQ et Q4G des mailles supports
quadrangulaires 4 nuds ( 45=20 degrs de libert).

Remarque importante :
Pour les lments de plaque 4 nuds ( DSQ , DKQ et Q4G ), les 4 nuds doivent tre
coplanaires pour que la thorie des plaques puisse tre valide. Cette vrification est
effectue systmatiquement par Code_Aster, et lutilisateur est alarm dans le cas o lun
des lments du maillage ne respecte pas cette condition.
Dans le cas des lments avec cisaillement transversal, pour viter le blocage des lments en
cisaillement transverse (surestimation de la rigidit pour des paisseurs trs faibles), une mthode
consiste construire des champs de cisaillement de substitution constants sur les bords de llment,
dont la valeur est lintgrale du cisaillement sur le bord en question. Dans le Code_Aster, les lments
de plaque et coque avec du cisaillement transverse utilisent cette mthode de faon ne pas bloquer
en cisaillement transverse. Ce blocage en cisaillement vient du fait que lnergie lastique de
cisaillement est un terme proportionnel h ( h tant lpaisseur de la plaque ou de la coque) bien
plus grand que le terme dnergie lastique de flexion qui est proportionnel en
h
3
. Quand lpaisseur
devient faible devant la longueur caractristique (le rapport h/ L est infrieur 1/ 20 ), pour
certaines fonctions de forme, la minimisation du terme prpondrant en h conduit une mauvaise
reprsentation des modes de flexion pure, pour lesquels la flche nest plus calcule correctement
(voir [bib1] page 295 avec h/ L=0.01 ).
Llment Q4G est un lment quadrilatral quatre nuds sans blocage en cisaillement transverse,
avec des fonctions de forme bilinaires en x et y pour reprsenter w ,
x
et
y
.
La principale diffrence entre la modlisation Q4G et DST vient du fait que lon utilise pour cette
dernire des fonctions de forme incluant un terme quadratique pour discrtiser les rotations
s
o
s=( x , y) .
La consquence pour l'lment Q4G est une approximation constante par morceaux des courbures qui
implique de mailler suffisamment fin dans les directions sollicites en flexion.
Remarque sur l'lment DST :
Des problmes ont t mis en vidence sur l'lment triangulaire DST dans le cas de
chargement en flexion. Les rsultats maillage quivalent sont nettement moins bons que ceux
du DSQ et dpendent fortement de la nature du maillage (voir cas-test ssls141 modlisations E
et F V3.03.141). On conseille aux utilisateurs :
de prfrer l'lment DSQ l'lment DST
sur les DST de prfrer les maillages libres aux maillages rgls (rsultats meilleurs)
d'viter les maillages disymtriques par rapport au chargement (s'il existe des symtries)
2.1.4.3 Les diffrences entre les lments coques
Code_Aster
Version
default
On distingue les lments de coques liniques COQUE_C_PLAN, COQUE_D_PLAN et axisymtrique
COQUE_AXIS des lments de COQUE_3D.
Les premiers sont utiliss pour modliser des structures invariantes selon laxe Oz ou de rvolution
daxe Oy et les seconds dans tous les autres cas. Dans le cas des coques invariantes suivant la
direction z , on distingue les coques libres en z (contraintes planes COQUE_C_PLAN) des coques
bloques en z (dformations planes COQUE_D_PLAN). Pour ces lments de coques, les mailles
supports sont liniques 3 nuds. Le nombre de degrs de libert de ces lments est de 9.

Les lments de coques quelconques COQUE_3D ont des mailles support triangulaires 7 nuds ou
quadrangulaires 9 nuds :
Dans le cas de mailles triangulaires, le nombre de degrs de libert pour les translations est
6 (les inconnues sont les dplacements aux nuds sommets et sur les milieux des cts du
triangle) et celui des rotations est 7 (les inconnues sont les 3 rotations aux points prcdents
et au centre du triangle). Le nombre de degrs de libert total de llment est donc de
Nddle=36+37=39 .
Dans le cas de mailles quadrangulaires 9 nuds, le nombre de degrs de libert pour les
translations est 8 (les inconnues sont les dplacements aux nuds sommets et sur les
milieux des cts du quadrangle) et celui des rotations est 9 (les inconnues sont les 3
rotations aux points prcdents et au centre du quadrangle). Le nombre de degrs de libert
total de llment est donc de Nddle=38+39=51 . Ces lments ont donc peu prs
deux fois plus de degrs de libert que les lments de plaque de la famille DKT
correspondants. Leur cot en temps, nombre gal, dans un calcul sera donc plus important.

Les lments de COQUE_3D prennent en compte automatiquement la correction de mtrique entre la
surface moyenne et les surfaces suprieure et infrieure. Pour les lments liniques, cette correction
doit tre active par lutilisateur (voir le paragraphe 17). La correction de mtrique apporte une
contribution en h/ L la contrainte et en (h/ L)
2
en dplacement (voir [V7.90.03]). Pour les plaques
cette correction est sans objet.

Pour les lments de coques le coefficient de correction de cisaillement k en comportement isotrope
peut tre modifi par lutilisateur. Ce coefficient de correction de cisaillement est donn dans
AFFE_CARA_ELEM sous le mot-cl A_CIS. Par dfaut, si lutilisateur ne prcise rien dans
AFFE_CARA_ELEM cela revient utiliser la thorie avec cisaillement de REISSNER ; le coefficient de
cisaillement est alors mis k=5/ 6 . Si le coefficient de cisaillement k vaut 1 on se place dans le
cadre de la thorie de HENCKY-MINDLIN-NAGHDI et sil devient trs grand (
k>10
6.
h/ L
) on se
rapproche de la thorie de LOVE-KIRCHHOFF.
En pratique il est conseill de ne pas changer ce coefficient. En effet, ces lments fournissent une
solution physiquement correcte, que la coque soit paisse ou mince, avec le coefficient k=5/ 6 .
Tous les lments de coque prsents dans Code_Aster sont protgs contre le verrouillage en
cisaillement comme les lments de plaques (voir 2.1.4.2). On utilise pour cela une intgration
slective [R3.07.03].
Cependant cette proprit est perdue si l'on projette un champ de dplacement sur un lment de
coque partir d'un calcul sur un lment de plaque (avec la commande PROJ_CHAMP). Il ne faut pas
s'tonner des contraintes de cisaillement obtenues ensuite par un calcul linaire partir des
dplacements (pour de faibles paisseurs).
2.1.4.4 La diffrence entre les coques volumiques SHB
Les lments SHB ont des mailles support et des formulations diffrentes :
llment SHB6 a pour maille support le pentadre six nuds. Sa formulation est linaire.
Son oprateur gradient discrtis est modifi pour amliorer sa vitesse de convergence ;
Code_Aster
Version
default
llment SHB8 a pour maille support lhexadre huit nuds. Sa formulation est linaire. Il
est stabilis par mthode assumed strain pour avoir de bonnes performances ;
llment SHB15 a pour maille support le pentadre quinze nuds. Sa formulation est
quadratique. Il na pas de stabilisation;
llment SHB20 a pour maille support lhexadre vingt nuds. Sa formulation est
quadratique. Il na pas de stabilisation.
2.2 Commandes utiliser
2.2.1 Discrtisation spatiale et affectation dune modlisation : oprateur
AFFE_MODELE
Dans cette partie, on dcrit le choix et laffectation dune des modlisations plaque ou coque ainsi que
les degrs de libert et les mailles associes. La plupart des informations dcrites sont extraites des
documentations dutilisation des modlisations ([U3.12.01] : Modlisation DKT - DST - Q4G,
[U3.12.02] : Modlisations COQUE_C_PLAN, COQUE_D_PLAN, COQUE_AXIS, [U3.12.05] :
Modlisation SHB8).

2.2.1.1 Degrs de liberts
Les degrs de libert de discrtisation sont en chaque nud de la maille support les composantes de
dplacement aux nuds de la maille support, sauf indication.

Modlisation Degrs de libert ( chaque
nud)
Remarques
COQUE_3D DX DY DZ DRX DRY DRZ
DRX DRY DRZ au nud central
Les nuds appartiennent au feuillet
moyen de la coque
DKT, DKTG DX DY DZ DRX DRY DRZ Les nuds appartiennent la facette
tangente au feuillet moyen de la coque
DST DX DY DZ DRX DRY DRZ Les nuds appartiennent la facette
Q4G DX DY DZ DRX DRY DRZ Les nuds appartiennent la facette
Q4GG DX DY DZ DRX DRY DRZ Les nuds appartiennent la facette
COQUE_C_PLAN DX DY DRZ Les nuds appartiennent la surface
moyenne de la coque
COQUE_D_PLAN DX DY DRZ Les nuds appartiennent la surface
moyenne de la coque
COQUE_AXIS DX DY DRZ Les nuds appartiennent la surface
moyenne de la coque
GRILLE_EXCENTRE DX DY DZ DRX DRY DRZ Les nuds appartiennent la facette
tangente au feuillet moyen de la
coque.
SHB8 DX DY DZ Degr de libert tous les nuds.
GRILLE_MEMBRANE DX DY DZ Degr de libert tous les nuds.
MEMBRANE DX DY DZ Degr de libert tous les nuds.
2.2.1.2 Mailles support des matrices de rigidit
Modlisation Maille lment fini Remarques
COQUE_3D TRIA7
QUAD9
MEC3TR7H
MEC3QU9H
Mailles non supposes planes
DKT, DKTG TRIA3
QUAD4
MEDKTR3
MEDKQU4
Mailles planes
DST TRIA3
QUAD4
MEDSTR3
MEDSQU4
Mailles planes
Q4G, QUAD4 MEQ4QU4 Mailles planes
Q4GG QUAD4 MEQ4GG4 Mailles planes
Code_Aster
Version
default
TRIA3 MET3GG3
COQUE_C_PLAN SEG3 METCSE3 Mailles non supposes planes
COQUE_D_PLAN SEG3 METDSE3 Mailles non supposes planes
COQUE_AXIS SEG3 MECXSE3 Mailles non supposes planes
GRILLE_EXCENTRE TRIA3,
QUAD4
MEGCTR3
MEGCQU4
Mailles planes
GRILLE_MEMBRANE TRIA3,
QUAD4,
TRIA6,
QUAD8
MEGMTR3
MEGMQU4
MEGMTR6
MEGMQU8
Mailles surfaciques quelconques
MEMBRANE TRIA3,
QUAD4,
TRIA6,
QUAD8
MEMBTR3
MEMBQU4
MEMBTR6
MEMBQU8
Mailles surfaciques quelconques
SHB8 PENTA6
HEXA8
PENTA15
HEXA20
MECA_SHB6
MECA_SHB8
MECA_SHB15
MECA_SHB20
Mailles 3D
Mailles 3D
Mailles 3D
Mailles 3D

La modlisation GRILLE_EXCENTRE utilise pour modliser les structures en bton arm a les
mmes caractristiques de maillage que la modlisation DKT

Remarque :
Dans un maillage, pour transformer des mailles TRIA6 en mailles TRIA7, ou QUAD8 en
mailles QUAD9 , on peut utiliser loprateur MODI_MAILLAGE .

2.2.1.3 Mailles support des chargements
Tous les chargements applicables aux facettes des lments utiliss ici sont traits par discrtisation
directe sur la maille support de llment en formulation dplacement. La pression et les autres forces
surfaciques ainsi que la pesanteur sont des exemples de chargements sappliquant directement sur
les facettes. Aucune maille spciale de chargement nest donc ncessaire pour les faces des
lments de plaques, de coques. Les lments SHB ont les mmes lments de bord que les
lments 3D standards.

Pour les chargements applicables sur les bords des lments, on a :

Modlisation Maille lment fini Remarques
COQUE_3D SEG3 MEBOCQ3
DKT, DKTG SEG2 MEBODKT
DST SEG2 MEBODST
Q4G SEG2 MEBOQ4G
Q4GG SEG2 MEBOQ4G
COQUE_C_PLAN POI1 Mailles supports rduites 1 point
COQUE_D_PLAN POI1 Mailles supports rduites 1 point
COQUE_AXIS POI1 Mailles supports rduites 1 point
GRILLE_EXCENTRE,
GRILLE_MEMBRANE,
MEMBRANE
Pas d'lment de bord affect par ces
modlisations.
SHB8 QUA4
TRI3
QUA8
TRI6
SEG2
SEG3
MECA_FACE4
MECA_FACE3
MECA_FACE8
MECA_FACE6
MECA_ARETE2
MECA_ARETE3
Les forces rparties, liniques, de traction, de cisaillement, les moments flchissant appliqus sur les
bords de structures coque entrent dans cette catgorie de chargements.
Code_Aster
Version
default
2.2.1.4 Modle : AFFE_MODELE
Laffectation de la modlisation passe par loprateur AFFE_MODELE [U4.41.01].
AFFE_MODELE Remarques
AFFE
PHENOMENE : 'MECANIQUE'
MODELISATION 'COQUE_3D'
'DKT'
'DST'
'DKTG'
'Q4G'
'Q4GG'
'COQUE_C_PLAN'
'COQUE_D_PLAN'
'COQUE_AXIS'
'GRILLE_MEMBRANE'
'GRILLE_EXCENTRE'
'MEMBRANE'
'SHB8'
Remarque :
Il est prudent de vrifier le nombre dlments affects.

2.2.2 Caractristiques lmentaires : AFFE_CARA_ELEM

Dans cette partie, les oprandes caractristiques des lments de plaques et de coques sont dcrites.
La documentation dutilisation de loprateur AFFE_CARA_ELEM est [U4.42.01].
AFFE_CARA_ELEM COQUE_3D DKT DKTG DST Q4G Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
COQUE
EPAIS
/ ANGL_REP
/ VECTEUR

A_CIS
COEF_RIGI_DRZ
EXCENTREMENT
INER_ROTA
MODI_METRIQUE
COQUE_NCOU
AFFE_CARA_ELEM
GRILLE_EXCENTRE GRILLE_MEMBRANE
GRILLE
SECTION
/ ANGL_REP
/ AXE,
ORIG_AXE

EXCENTREMENT
COEF_RIGI_DR
Z
Code_Aster
Version
default
AFFE_CARA_ELEM
MEMBRANE
MEMBRANE
/ ANGL_REP
/ AXE, ORIG_AXE
Les caractristiques pouvant tre affectes sur les lments de plaque ou de coque sont :
Lpaisseur EPAIS constante sur chaque maille, puisque le maillage ne reprsente que le
feuillet moyen.
Le coefficient de correction du cisaillement transverse A_CIS pour les coques courbes
isotropes.
La prise en compte de la correction de mtrique MODI_METRIQUE entre la surface moyenne
et les surfaces suprieure et infrieure (effective uniquement pour les COQUE_C_PLAN,
COQUE_D_PLAN, COQUE_AXIS).
Une direction de rfrence permettant de dfinir un repre local dans le plan tangent en tout
point d'une coque. La construction du repre local se fait soit l'aide des deux angles
nautiques
o
et
(fournis en degrs) qui dfinissent un vecteur v dont la projection

sur le plan tangent la coque fixe la direction x
l
. soit , si le mot clef VECTEUR est prsent,
par les 3 composantes du vecteur v On peut dfinir un unique vecteur V pour toute la
structure, ou bien un par zone (mots cls GROUP_MA / MAILLE). La construction du repre
local est dfinie dans AFFE_CARA_ELEM [U4.42.01]. On appellera par la suite ce repre local
repre utilisateur.
La dfinition de cet axe de rfrence est utile galement pour dfinir l'orientation des fibres
d'une coque multicouche ou orthotrope (Cf. oprateur DEFI_COMPOSITE [U4.42.03]).
Le nombre de couches COQUE_NCOU utilises pour lintgration dans l'paisseur de la coque,
dans les oprateurs STAT_NON_LINE et DYNA_NON_LINE (modlisations DKT, COQUE_3D,
COQUE_AXIS, COQUE_C_PLAN, COQUE_D_PLAN).

Une fonctionnalit de DEFI_GROUP permet de crer automatiquement un groupe de mailles
dont la normale est comprise dans un angle solide donn, daxe la direction de rfrence.

Cette commande peut tre utilise en pr-traitement pour affecter des donnes matriau non
isotropes ou en post-traitement aprs un calcul de coque.
Lexcentrement (constant pour tous les nuds de la maille) EXCENTREMENT de chacune
delles par rapport la maille support. Cette distance est mesure sur la normale de la maille
support. Dans le cas excentr les inerties de rotation sont obligatoirement prises en compte
et INER_ROTA est mis OUI.
COEF_RIGI_DRZ dfinit un coefficient de rigidit fictive (ncessairement petit) sur le degr
de libert de rotation autour de la normale la coque. Il est ncessaire pour empcher que la
matrice de rigidit soit singulire, mais doit tre choisi le plus petit possible. La valeur par
dfaut ( 1.E5 ) convient pour la plupart des situations (cest une valeur relative : la rigidit
autour de la normale est gale KRZ fois le plus petit terme diagonal de la matrice de
rigidit de llment).

Code_Aster
Version
default
Figure 2.2.2-a : Repre global et plan tangent

Pour les modlisations GRILLE_EXCENTRE et GRILLE_MEMBRANE,
Les donnes gomtriques suivantes sont ncessaires pour modliser la nappe darmatures :
SECTION = S
1
: section des armatures dans la direction 1. La section est donne par mtre
linaire. Elle correspond donc la section cumule sur une largeur de 1 mtre. S'il y a une
section s tous les 20.0cm , la section cumule est 5.s .
L'orientation des armatures est dfinie soit par
ANGL_REP, pour dfinir un vecteur projet sur l'lment
soit dans le cas dune coque cylindrique, par ORIG_AXE, AXE pour dfinir langle des
armatures, constant dans un repre cylindrique.
Lexcentrement (constant pour tous les nuds de la maille) de la nappe darmatures par rapport
la maille support (distance mesure sur la normale de la maille support), (modlisation
GRILLE_EXCENTRE uniquement).
COEF_RIGI_DRZ dfinit un coefficient de rigidit fictive (ncessairement petit) sur le degr de
libert de rotation autour de la normale la coque. Il est ncessaire pour empcher que la
matrice de rigidit soit singulire, mais doit tre choisi le plus petit possible. La valeur par dfaut
(1.E-5) convient pour la plupart des situations (cest une valeur relative : la rigidit autour de la
normale est gale KRZ fois le plus petit terme diagonal de la matrice de rigidit de
llment).
Pour dfinir une grille ou si la section des armatures dans le sens longitudinal et dans le transversal sont
diffrents, il faut crer deux couches d'lments (commande CREA_MAILLAGE, mot cl
CREA_GROUP_MA), une couche d'lment pour la direction longitudinale et une deuxime couche
d'lments pour la direction transversale :
Pour la modlisation MEMBRANE, seule l'orientation du comportement de la membrane est ncessaire.
Elle est dfinie soit par :
ANGL_REP,pour dfinir un vecteur projet sur l'lment
soit dans le cas dune coque cylindrique, par ORIG_AXE , AXE pour dfinir langle
des armatures, constant dans un repre cylindrique.
Le sens des normales chaque lment est un problme rcurrent concernant lutilisation de ce type
dlment, par exemple lorsque lon applique des chargements de type pression, ou bien pour dfinir
un excentrement ou un repre local.
Par dfaut pour les lments surfaciques lorientation est donne par le produit vectoriel 12^13 pour
un triangle numrot 123 (DKT,...) ou 1234567 (COQUE_3D) et 12^14 pour un quadrangle numrot
1234 (DKQ,...) ou 123456789 (COQUE_3D). Pour les coques liniques n est donne par la formule
du paragraphe 2.1.1.1 avec t orient dans le sens de parcours de la maille au niveau du maillage.
Gnralement, ces donnes sont accessibles en regardant dans le fichier de maillage, ce qui nest pas
trs pratique pour lutilisateur. En outre, il faut quil vrifie la cohrence de son maillage et sassurer
que toutes les mailles ont bien la mme orientation.

Lutilisateur peut modifier automatiquement lorientation des lments du maillage en imposant une
direction de normale, pour un maillage ou une partie de maillage utilisant des modlisations de coque
et quelle que soit le type de modlisation. La rorientation des lments se fait par le biais de
loprateur ORIE_NORM_COQU de la commande MODI_MAILLAGE [U4.12.05]. Le principe est le
suivant : on dfinit sous ORIE_NORM_COQU une direction par le biais dun vecteur et un nud
appartenant au groupe de mailles rorienter. Si le vecteur introduit nest pas dans le plan de la
maille slectionne par MODI_MAILLAGE, on en dduit automatiquement une direction de normale
obtenue comme le vecteur donn moins sa projection sur le plan de la maille. Toutes les mailles du
groupe connexes celles initialement slectionne auront alors la mme orientation de normale
automatiquement. Par ailleurs une vrification automatique de la mme orientation des mailles
Code_Aster
Version
default
connexes est effectue par le biais de loprateur VERI_NORM de la commande AFFE_CHAR_MECA
[U4.25.01].
Comme les coques, les lments coques volumiques SHB ont une direction privilgie, appele
paisseur. Cette direction de lpaisseur est dfinie par la faon de mailler. En gnral, on a besoin de
deux surfaces pour gnrer un volume tridimensionnel. La direction qui est perpendiculaire ces deux
surfaces est la direction de lpaisseur. Pour vrifier la bonne orientation des faces des lments
dsigns (compatibilit avec la direction privilgie), on utilise ORIE_SHB de loprateur
MODI_MAILLAGE . Si lorientation de llment nest pas bonne, loprateur MODI_MAILLAGE
renumrote les nuds de llment.
2.2.3 Matriaux : DEFI_MATERIAU
La dfinition du comportement d'un matriau seffectue laide de loprateur DEFI_MATERIAU
[U4.43.01].

DEFI_MATERIAU COQUE_3D
DKT,
DST,
Q4G
DKTG,
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_EXCENTRE,
GRILLE_MEMBRANE
SHB8 MEMBRANE
ELASTICITE .
LINEAIRE
ELAS
ELAS_ORTH
ELAS_ISTR

ELAS_COQUE
ELAS_MEMBRANE

DEFI_MATERIAU
COQUE_3
D
DKT
DST
Q4G
DKT
G
Q4GG
COQUE_
C_PLAN
COQUE_
D_PLAN
COQUE_
AXIS
GRILLE_EXCENTRE
GRILLE_MEMBRANE
SHB8 MEMBRANE
PAR COUCHE
comportements
disponibles
en C_PLAN

comportements
disponibles
en 1D
GLRC_DAMAGE
GLRC_DM,
KIT_DDI
Les matriaux utiliss avec lensemble des lments plaques ou coques peuvent avoir des
comportements lastiques en contraintes planes dont les caractristiques linaires sont constantes ou
fonctions de la temprature.

Tous les comportements non linaires en contraintes planes (soit directement, soit par la mthode de
De Borst [R5.03.03]) sont disponibles pour les modlisations DKT, SHB8 et coques. Pour plus
dinformations sur ces non linarits on peut se reporter au paragraphe [27].

Tous les comportements non linaires en 1D (soit directement, soit par la mthode de De Borst) sont
disponibles pour les modlisations GRILLE_EXCENTRE et GRILLE_MEMBRANE.
Les structures minces en matriaux composites ne peuvent tre traites actuellement que par les
modlisations plaques, en utilisant DEFI_COMPOSITE avec des caractristiques de matriaux
homogniss. On peut aussi introduire directement les coefficients de rigidit des matrices de
membrane, flexion et cisaillement avec ELAS_COQUE. Ces coefficients sont donns dans le repre
utilisateur de llment dfini par ANGL_REP. Il est noter que les termes de cisaillement ne sont pris
en compte avec le comportement ELAS_COQUE que pour les lments DST et Q4G. Ils ne sont pas pris
en compte avec les lments DKT.
Code_Aster
Version
default
Afin de faciliter la comprhension, nous avons reprsent sur la figure ci-dessous les diffrents
repres utiliss.

Figure 2.2.3-a : Repres utiliss pour la dfinition du matriau
Lexemple suivant est extrait du cas-test SSLS117B et illustre la syntaxe de DEFI_COMPOSITE :
MU2=DEFI_COMPOSITE(COUCHE=_F(EPAIS=0.2,
MATER=MAT1B,
ORIENTATION=0.0,),);
Dans cet exemple, on dfinit un composite multicouche dpaisseur 0.2 , le matriau tant dfinit par
MAT1B, et langle de la 1re direction d'orthotropie (sens longitudinal ou sens des fibres) tant nul. On
se reportera la documentation [U4.42.03] pour plus de dtails concernant lutilisation de
DEFI_COMPOSITE. (voir aussi [R4.01.01].
2.2.4 Chargements et conditions limites : AFFE_CHAR_MECA et AFFE_CHAR_MECA_F
Laffectation des chargements et des conditions aux limites sur un modle mcanique seffectue
laide de loprateur AFFE_CHAR_MECA, si les chargements et les conditions aux limites mcaniques
sur un systme sont des valeurs relles ne dpendant daucun paramtre, ou bien
AFFE_CHAR_MECA_F, si ces valeurs sont fonctions de la position ou de lincrment de chargement.
La documentation dutilisation dAFFE_CHAR_MECA et AFFE_CHAR_MECA_F est [U4.44.01].

2.2.4.1 Liste des mot-cls facteur d AFFE_CHAR_MECA

Code_Aster
Version
default
AFFE_CHAR_MECA COQUE_3D DKT,
DKTG
DST Q4G,
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_* SHB8
DDL_IMPO
FACE_IMPO
LIAISON_DDL
LIAISON_OBLIQUE
LIAISON_GROUP
CONTACT
LIAISON_UNIF
LIAISON_SOLIDE
LIAISON_ELEM
LIAISON_COQUE
FORCE_NODALE
DDL_IMPO Mot cl facteur utilisable pour imposer, des nuds ou des groupes de
nuds, une ou plusieurs valeurs de dplacement.
FACE_IMPO Mot cl facteur utilisable pour imposer, tous les nuds dune face dfinie
par une maille ou un groupe de mailles, une ou plusieurs valeurs de
dplacements (ou de certaines grandeurs associes).
LIAISON_DDL Mot cl facteur utilisable pour dfinir une relation linaire entre des degrs
de libert de deux ou plusieurs nuds.
LIAISON_OBLIQUE Mot cl facteur utilisable pour appliquer, des nuds ou des groupes de
nuds, la mme valeur de dplacement dfinie composante par
composante dans un repre oblique quelconque.
LIAISON_GROUP Mot cl facteur utilisable pour dfinir des relations linaires entre certains
degrs de libert de couples de nuds, ces couples de nuds tant obtenus
en mettant en vis--vis deux listes de mailles ou de nuds.
CONTACT Mot cl facteur utilisable pour notifier des conditions de contact et de
frottement entre deux ensembles de mailles.
LIAISON_UNIF Mot cl facteur permettant dimposer une mme valeur (inconnue) des
degrs de libert dun ensemble de nuds.
LIAISON_SOLIDE Mot cl facteur permettant de modliser une partie indformable dune
structure.
LIAISON_ELEM Mot cl facteur qui permet de modliser les raccords dune partie coque avec
une partie poutre ou d'une partie coque avec une partie tuyau (voir
paragraphe 2.2.4.5).
LIAISON_COQUE Mot cl facteur permettant de reprsenter le raccord entre des coques au
moyen de relations linaires.
FORCE_NODALE Mot cl facteur utilisable pour appliquer, des nuds ou des groupes de
nuds, des forces nodales, dfinies composante par composante dans le
repre GLOBAL ou dans un repre oblique dfini par 3 angles nautiques.

Code_Aster
Version
default
AFFE_CHAR_MECA
particuliers
COQUE_3D
DKT,
DKTG
DST,
Q4G
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_* SHB8
FORCE_ARETE
FORCE_COQUE
global
pres
local
tangent
PESANTEUR
PRES_REP
ROTATION
PRE_EPSI
FORCE_ARETE Mot cl facteur utilisable pour appliquer des forces liniques sur une arte
dun lment de coque. Pour les lments liniques lquivalent revient
appliquer une force nodale aux nuds supports de llment. Il ny a donc
pas de terme ddi particulier. En revanche, il ncessite des lments de
bords.
FORCE_COQUE Mot cl facteur utilisable pour appliquer des efforts surfaciques (pression par
exemple) sur des lments dfinis sur tout le maillage ou sur une ou
plusieurs mailles ou des groupes de mailles. Ces efforts peuvent tre donns
dans le repre global ou dans un repre de rfrence dfini sur chaque
maille ou groupe de mailles; ce repre est construit autour de la normale
llment de coque et dune direction fixe (voir paragraphe 2.2.2).
PESANTEUR Mot cl facteur utilisable pour un chargement de type pesanteur.
PRES_REP Mot cl facteur utilisable pour appliquer une pression sur une ou plusieurs
mailles, ou des groupes de mailles.
ROTATION Mot cl facteur utilisable pour calculer le chargement d la rotation de la
structure.
PRE_EPSI Mot cl facteur utilisable pour appliquer un chargement de dformation
initiale.
Remarque :
Les efforts de pression sexerant sur les lments de plaques peuvent sappliquer soit par
FORCE_COQUE ( pres ) soit par PRES_REP . Lutilisateur devra donc faire attention ne pas
appliquer deux fois le chargement de pression pour les lments concerns, surtout dans
les cas o les modlisations de plaques seraient mlanges avec dautres modlisations
utilisant PRES_REP .
Par ailleurs il faut noter que les efforts de pression, que ce soit avec FORCE_COQUE (pres) ou
PRES_REP sont tels quune pression positive agit dans le sens contraire celui de la normale
llment. Par dfaut, cette normale est dpendante du sens de parcours des nuds dun lment, ce
qui nest pas toujours trs facile pour lutilisateur.
En outre il faut que celui-ci sassure que tous ces lments sont orients de la mme manire. On
conseille donc dimposer lorientation de ces lments par le biais de loprateur ORIE_NORM_COQU
de la commande MODI_MAILLAGE (voir paragraphe [2.2.2]).
2.2.4.2 Liste des mot-cls facteur d AFFE_CHAR_MECA_F
Les mot-cls facteur gnraux de loprateur AFFE_CHAR_MECA_F sont identiques ceux de
loprateur AFFE_CHAR_MECA prsents ci-dessus.
Code_Aster
Version
default
AFFE_CHAR_MECA_F
particuliers
COQUE_3D DKT,
DKTG
DST Q4G,
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
FORCE_ARETE
FORCE_COQUE global
pres
local tangent
Les chargements de pression fonctions de la gomtrie peuvent tre renseigns par le biais de
FORCE_COQUE (pres).
2.2.4.3 Application dune pression : mot-cl FORCE_COQUE
Le mot-cl facteur FORCE_COQUE permet dappliquer des efforts surfaciques sur des lments de type
coque (DKT, DST, Q4G, ) dfinis sur tout le maillage ou sur une ou plusieurs mailles ou des groupes
de mailles. Suivant le nom de l'oprateur appel, les valeurs sont fournies directement
(AFFE_CHAR_MECA) ou par l'intermdiaire d'un concept fonction (AFFE_CHAR_MECA_F).
AFFE_CHAR_MECA
AFFE_CHAR_MECA_F
Remarques
FORCE_COQUE :
Repre
global
TOUT: 'OUI'
MAILLE
GROUP_MA
Lieu dapplication du chargement
FX
FY
FZ
MX
MY
MZ
Fourni directement pour
AFFE_CHAR_MECA, sous forme de
fonction pour AFFE_CHAR_MECA_F
PLAN 'MOY'
'INF'
'SUP'
'MAIL'
Permet de dfinir un torseur defforts
sur le plan moyen, infrieur, suprieur
ou du maillage (lments DKT et DST)
Repre
local
PRES F1
F2
F3
MF1
MF2
Fourni directement pour
AFFE_CHAR_MECA, sous forme de
fonction pour AFFE_CHAR_MECA_F
Nous renvoyons au paragraphe correspondant au mot-cl FORCE_COQUE du document dutilisation
des oprateurs AFFE_CHAR_MECA et AFFE_CHAR_MECA_F.
2.2.4.4 Conditions limites : mots-cls DDL_IMPO et LIAISON_*
Le mot-cl facteur DDL_IMPO permet dimposer, des nuds introduits par un (au moins) des mots
cls : TOUT, NOEUD, GROUP_NO, MAILLE, GROUP_MA, une ou plusieurs valeurs de dplacement (ou de
certaines grandeurs associes). Suivant le nom de l'oprateur appel, les valeurs sont fournies
directement (AFFE_CHAR_MECA) ou par l'intermdiaire d'un concept fonction (AFFE_CHAR_MECA_F).
Les oprandes disponibles pour DDL_IMPO, sont lists ci-dessous :
DX DY DZ Blocage sur la composante de dplacement en translation
DRX DRY DRZ Blocage sur la composante de dplacement en rotation
2.2.4.5 Raccords coques avec autres lments mcaniques
Code_Aster
Version
default
Ces raccords doivent satisfaire les exigences tablies dans [bib4] et que lon retrouve en particulier
dans le raccord 3D-POUTRE en [R3.03.03].

Les raccords disponibles avec les lments de plaques et de coques sont les suivants :

Raccord Poutre-Coque : il sagit d'tablir la liaison entre un nud extrmit dun lment
de poutre et un groupe de mailles de bord dlments de coques. Les thories de poutre et
de plaque ne connaissent que des coupures normales la fibre ou la surface moyenne. Les
raccords ne peuvent avoir lieu que suivant ces fibres ou surfaces moyennes. Le raccord
poutre-coque est ralisable pour des poutres dont la fibre neutre est orthogonale aux
normales aux facettes des plaques ou des coques. Ltendre dautres configurations
(une poutre arrivant perpendiculairement au plan dune plaque par exemple) demande une
tude de faisabilit car les lments de plaque ou de coque nont pas de rigidit associe
une rotation dans le plan perpendiculaire la normale la surface moyenne. Le raccord est
utilisable en utilisant le mot-cl LIAISON_ELEM : (OPTION : COQ_POU) de
AFFE_CHAR_MECA.
Raccord Coque-Tuyau : il sagit d'tablir la liaison entre un nud extrmit dun lment de
tuyau et un groupe de maille de bord dlments de coques. La formulation du raccord
coques-tuyaux est prsent dans le document de rfrence [R3.08.06]. Les thories de tuyau
et de plaque, ne connaissent que des coupures normales la fibre ou la surface moyenne.
Les raccords ne peuvent avoir lieu que suivant ces fibres ou surfaces moyennes. Le raccord
coque-tuyau est ralisable pour des tuyaux dont la fibre neutre est orthogonale aux normales
aux facettes des plaques ou des coques. Le raccord est utilisable en utilisant le mot-cl
LIAISON_ELEM : (OPTION : COQ_TUYAU) de AFFE_CHAR_MECA.

Figure 2.2.4.5-a : Raccords coques avec autres lments mcaniques
Raccord Coque 3D massif : le raccord coque-3D massif est ltude mais il sera limit
dans un premier temps aux cas o la normale au solide est orthogonale la normale lune
des facettes de llment de plaque ou de coque (voir [bib4]).
Raccord entre lments de Coques : pour raccorder deux lments de coques entre eux,
on utilise le mot-cl LIAISON_COQUE dAFFE_CHAR_MECA (_F) (documentation
[U4.44.01]). Ce raccord est effectu au moyen de relations linaires. Lapproche classique
admet que 2 plans maills en coques se coupent selon une droite qui appartient au maillage
de la structure. Afin dviter que le volume qui est lintersection des deux coques soit compt
deux fois, on arrte le maillage dune coque perpendiculaire une coque donne au niveau
de la peau suprieure ou infrieure de cette dernire. Sur la [Figure 2.2.4.5-b], la liaison
entre les 2 coques se fait par des liaisons de corps solide entre les nuds en vis--vis des
segments A
1
A
2
et B
1
B
2
.
Raccord Coques volumiques SHB 3D massif : ce raccord est naturel car les lments
coques volumiques SHB ont des mailles support 3D aussi. Il faut seulement assurer que tous
Code_Aster
Version
default
les nuds de la surface de raccord doivent appartenir llment coque volumique SHB et
3D massif.
Par consquent, les lments SHB de formulation linaire (SHB6, SHB8) ne peuvent pas se
raccorder aux lments 3D massifs de formulation quadratique (PENTA15, HEXA20 par
exemple) et inversement, les SHB15, SHB20 ne peuvent pas se raccorder aux lments 3D
massifs de formulation linaire.

Figure 2.2.4.5-b : Raccord entre lments de coques
Des cas-tests permettant de valider ces raccords sont disponibles dans la section exemples.

2.2.4.6 Variables de commandes

Les variables de commandes prises en compte par les diffrentes modlisations sont listes ici :
Variables
de
commandes
COQUE_3
D
DKT,
DST,
Q4G
DKTG Q4GG SHB8 COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_MEMBRANE
GRILLE_EXCENTRE
MEMBRANE
TEMP
autres :
SECH,
HYDR,
etc..
2.3 Rsolution
2.3.1 Calculs linaires : MECA_STATIQUE et autres oprateurs linaires
Les calculs linaires seffectuent en petites dformations. Plusieurs oprateurs de rsolution linaires
sont disponibles :
MECA_STATIQUE : rsolution dun problme de mcanique statique
linaire ([U4.51.01]) ;
MACRO_ELAS_MULT : calcule des rponses statiques linaires pour
diffrents cas de charges ou modes de Fourier.
([U4.51.02]).
MODE_ITER_SIMULT : calcul des valeurs et vecteurs propres par des
mthodes de sous espaces. ([U4.52.03]).
MODE_ITER_INV : calcul des valeurs et vecteurs propres par la
mthode ditrations inverses ([U4.52.04]).
MODE_ITER_CYCL : calcul des modes propres dune structure
symtrie cyclique ([U4.52.05]) ;
Code_Aster
Version
default
DYNA_LINE_TRAN : calcul de la rponse dynamique transitoire une
excitation temporelle quelconque ([U4.53.02]) ;
DYNA_TRAN_MODAL : calcul est effectu par superposition modale ou par
sous-structuration ([U4.53.21]) ;

2.3.2 Calculs non linaires : STAT_NON_LINE et DYNA_NON_LINE
2.3.2.1 Comportements et hypothses de dformations disponibles

Les informations suivantes sont extraites de la documentation dutilisation de loprateur
STAT_NON_LINE : [U4.51.03].
SHB8 COQUE
_3D
DKT DKTG DST,
Q4G
DKTG Q4GG GRILL
E_*
COQUE_
C_PLAN
COQUE_
D_PLAN
COQUE_
AXIS
COMP_INCR
(petites
dformations)
RELATION Toutes les relations
disponibles en
contraintes planes

Les relations 3D en
utilisant la mthode
de De Borst

Toutes les relations
disponibles en 1D
DEFORMATION:
'GROT_GDEP'
Coque_3D en grands
dplacements et
grandes rotations
disponible avec
comportements
incrmentaux non
linaires, mais en
petites dformations
DEFORMATION:
'PETIT' (ou
GROT_GDEP)
En petits ou grands
dplacements
disponible avec
comportements
incrmentaux non
linaires, mais en
rotations faibles et en
petites dformations

DEFORMATION
: PETIT ou
GROT_GDEP
En petits ou grands
dplacements
RELATION GLRC_DAMAGE
RELATION GLRC_DM, KIT_DDI
COMP_ELAS
(grands
dplacements,
grandes
rotations)
RELATION ELAS
DEFORMATION:
'GROT_GDEP'
TYPE_CHARGE:
'SUIV'
Pression suiveuse

Tous les comportements non linaires mcaniques de contraintes planes du code sont accessibles.
On distingue les relations de comportement incrmentales (mot cl facteur COMP_INCR) des relations
Code_Aster
Version
default
de comportements lastiques non linaires (mot cl facteur COMP_ELAS). La relation de
comportement relie les taux de dformation aux taux de contraintes.

Pour les modlisations GRILLE_EXCENTRE et GRILLE_MEMBRANE, pour des structures en bton
arm, les comportements non linaires 1D correspondent des comportements incrmentaux
particuliers dans STAT_NON_LINE (COMP_INCR) :
GRILLE_ISOT_LINE pour la plasticit crouissage isotrope,
GRILLE_ISOT_CINE pour la plasticit crouissage cinmatique bi linaire,
GRILLE_PINTO_MEN pour le comportement de Pinto Menegotto.
Les comportements 3D peuvent galement tre utiliss l'aide de la mthode de De Borst [R5.03.09].
La modlisation MEMBRANE n'est implmente que pour un comportement lastique en petites
dformations et petits dplacements.

Le concept RESULTAT de STAT_NON_LINE/DYNA_NON_LINE contient des champs de dplacements,
de contraintes et de variables internes aux points dintgration toujours calcules aux points de
gauss :
DEPL : champs de dplacements.
SIEF_ELGA : Tenseur des contraintes par lment aux points dintgration (COQUE_3D et
DKT) dans le repre utilisateur. Pour chaque couche, on stocke dans lpaisseur et pour
chaque paisseur sur les points dintgration surfacique. Ainsi si on veut des informations sur
une contrainte pour la couche NC, au niveau NCN (NCN = -1 si infrieur, NCN = 0 si milieu,
NCN = +1 si suprieur) pour le point dintgration surfacique NG, il faudra regarder la valeur
donne par le point dfini dans loption POINT tel que : NP = 3*(NC-1)*NPG+
(NCN+1)*NPG+NG o NPG est le nombre total de points dintgration surfacique de llment
de COQUE_3D (7 pour le triangle et 9 pour le quadrangle) et de llment DKT. Pour les
modlisations GRILLE_EXCENTRE,GRILLE_MEMBRANE, on stocke simplement une valeur
par point d'intgration : la composante SIXX dans la direction des armatures. Pour les
modlisations DKTG et Q4GG, SIEF_ELGA contient les 6 efforts gnraliss (efforts
membranaires, moments flchissants, efforts tranchants) par point de Gauss. Pour la
modlisation SHB8, SIEF_ELGA contient les 6 composantes des contraintes par point
d'intgration, en repre local, plus 12 composantes relatives aux termes de stabilisation.

VARI_ELGA : Champ de variables internes (DKT et COQUE_3D) par lment aux points
dintgration surfacique. Pour chaque point dintgration surfacique, on stocke les
informations sur les couches en commenant par la premire, niveau INF. Le nombre de
variables reprsentes vaut donc 2*NCOU*NBVARI o NBVARI reprsente le nombre de
variables internes.
Il peut tre enrichit des champs suivants, calculs en post-traitement par loprateur CALC_CHAMP :
EFGE_ELNO : active le calcul du tenseur des efforts gnraliss par lment aux nuds
(efforts membranaires, moments flchissants, efforts tranchants), dans le repre utilisateur
(dfini au paragraphe [2.2.2]).
VARI_ELNO : active le calcul du champ de variables internes par lment aux nuds dans
lpaisseur (par couche SUP/MOY/INF dans lpaisseur sauf indication).
2.3.2.2 Dtail sur les points dintgration
La numrotation des points de Gauss surfaciques pour les lments COQUE_3D est donne en
R3.07.04 4.7.1. Attention l'ordre des points de Gauss pour la formule 9 points, qui n'est pas le
mme que celui adopt pour les lments isoparamtriques.
Code_Aster
Version
default
Pour les modlisations DKT, COQUE_3D, COQUE_D_PLAN, COQUE_C_PLAN, COQUE_AXIS, dans le cas
des calculs non-linaires, la mthode dintgration pour les lments de plaque et de coques est une
mthode dintgration par couches, dont le nombre est dfini par lutilisateur. Pour chaque couche,
sauf modlisation GRILLE, on utilise une mthode de Simpson trois points dintgration, au milieu
de la couche et en peaux suprieure et infrieure de couche. Pour N couches le nombre de points
dintgration dans lpaisseur est de 2N+1 .

Pour traiter les non-linarits matriau, on conseille dutiliser de 3 5 couches dans lpaisseur pour
un nombre de points dintgration valant 7, 9 et 11 respectivement. Pour la rigidit tangente, on
calcule pour chaque couche, en contraintes planes, la contribution aux matrices de rigidit de
membrane, de flexion et de couplage membrane-flexion. Ces contributions sont ajoutes et
assembles pour obtenir la matrice de rigidit tangente totale. Pour chaque couche, on calcule ltat
des contraintes et lensemble des variables internes, au milieu de la couche et en peaux suprieure et
infrieure de couche. Ces informations sont disponibles dans VARI_ELGA et SIEF_ELGA. Le
comportement plastique ne comprend pas les termes de cisaillement transverses qui sont traits de
faon lastique, car le cisaillement transverse est dcoupl du comportement plastique.
Pour les modlisations GRILLE_EXCENTRE et GRILLE_MEMBRANE des structures en bton arm, il
ny a quun point dintgration par couche.

Les lments coques SHB n'ont pas de notion de couche comme les lments de plaque ou de
COQUE_3D. Ils sont intgrs sur tous leurs points de Gauss. Les contraintes sont calcules aux points
dintgration dans le repre local.
2.3.2.3 Comportement non-linaire gomtrique
Les calculs en non-linaire gomtrique (grands dplacements et grandes rotations), disponibles avec
la modlisation COQUE_3D, sont ralises l'aide de l'oprateur STAT_NON_LINE, en utilisant, sous le
mot-cl COMP_ELAS / COMP_INCR, DEFORMATION = 'GROT_GDEP'
Les calculs en non-linaire gomtrique (grands dplacements et petites rotations), disponibles avec
la modlisation SHB8 , sont ralises l'aide de l'oprateur STAT_NON_LINE , en utilisant sous le
mot-cl COMP_INCR , DEFORMATION = 'GROT_GDEP'.
Les calculs en non-linaire gomtrique (grands dplacements et petites dformations), disponibles
avec la modlisation DKTG , sont ralises l'aide de l'oprateur STAT_NON_LINE , en utilisant, sous
le mot-cl COMP_INCR , DEFORMATION = 'GROT_GDEP'.

Il est possible d'appliquer sur les lments de COQUE_3D et de SHB, des pressions suiveuses. Ce
chargement a la particularit de suivre la gomtrie de la structure au cours de sa dformation (par
exemple : la pression hydrostatique reste toujours perpendiculaire la gomtrie dforme). Pour
prendre en compte ce type de chargement, il faut prciser dans l'oprateur STAT_NON_LINE les
informations suivantes :
STAT_NON_LINE (
EXCIT =_F( CHARGE = pres
TYPE_CHARGE = 'SUIV' )
)

Le comportement non-linaire gomtrique des structures peut prsenter des instabilits (flambement,
snap-through/snap-back...). La dtermination et le passage de ces points limites, ne peuvent pas tre
obtenus en imposant le chargement, cependant les options de pilotage du chargement DDL_IMPO
ou LONG_ARC de l'oprateur STAT_NON_LINE permettent de franchir ces points critiques.
2.3.2.4 Flambement linaire
Les calculs en flambement linaire sont similaires la recherche de frquences propres et de modes
de vibration. Le problme a rsoudre sexprime sous la forme :
Code_Aster
Version
default
Trouver
(\, X)(,
N
)
tels que
AX=\ BX
o
A
est la matrice de rigidit
B
est la matrice de rigidit gomtrique (calcule avec loption RIGI_GEOM de
CALC_MATR_ELEM), disponible pour les modlisations COQUE_3D et SHB8
\
est la charge critique
X
est le mode de flambement associ a la charge critique

Les oprateurs MODE_ITER_INV [U4.52.04] et MODE_ITER_SIMULT [U4.52.03] sont utiliss
pour dterminer la charge critique et le mode de flambement associ.

2.4 Calculs supplmentaires et post-traitements
2.4.1 Calculs de matrices lmentaires : oprateur CALC_MATR_ELEM

Loprateur CALC_MATR_ELEM (documentation [U4.61.01]) permet de calculer des matrices
lmentaires, qui sont ensuite assemblables par la commande ASSE_MATRICE (documentation
[U4.61.22]).

Les options lmentaires de loprateur CALC_MATR_ELEM sont dcrites ci-dessous :

CALC_MATR_ELEM COQUE_3
D
DKT,
DKTG
DST Q4G,
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_
*
SHB8
AMOR_MECA
MASS_MECA
'MASS_INER'
RIGI_GEOM
RIGI_MECA
RIGI_MECA_HYST
AMOR_MECA : Matrice damortissement des lments calcule par combinaison linaire de la
rigidit et de la masse.
MASS_MECA : Matrice de masse.
MASS_INER: calcul des caractristiques inertielles (masse, centre de gravit)
RIGI_GEOM : Matrice de rigidit gomtrique ( pour les grands dplacements).
RIGI_MECA : Matrice de rigidit des lments.
RIGI_MECA_HYST : Rigidit hystrtique (complexe) calcule par le produit par un
coefficient complexe damortissement structural de la rigidit simple.
2.4.2 Calculs par lments : oprateurs CALC_CHAMP et POST_CHAMP
On prsente ci-aprs les options de post-traitement pour les lments de plaques et de coques. Elles
correspondent aux rsultats que peut obtenir un utilisateur aprs un calcul thermomcanique
(contraintes, dplacements, dformations, variables internes, etc...). Pour les structures modlises
par des lments de coques ou de poutres, il est particulirement important de savoir comment sont
prsents les rsultats de contraintes de faon pouvoir les interprter correctement. Lapproche
adopte dans Code_Aster consiste calculer les contraintes dans le repre utilisateur dfini dans
l'oprateur AFFE_CARA_ELEM.
Si l'on souhaite dpouiller ses rsultats dans un autre repre que le repre utilisateur, il faut utiliser la
commande MODI_REPERE .
Lorsqu'un post-traitement ne concerne qu'un sous-point , l'utilisateur dispose des mots-cls
NUME_COUCHE et NIVE_COUCHE du mot cl facteur EXTR_COQUE de la commande POST_CHAMP.
Les mots-cls sous EXTR_COQUE sont dcrits dans le tableau suivant :
Code_Aster
Version
default
OPTIONS COQUE_3D DKT DST,
Q4G
DKTG,
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_* SHB8
NUME_COUCHE
NIVE_COUCHE
Plus prcisment, dans le cas d'un matriau multicouche (coque multicouche dfinie par
DEFI_COMPOSITE), ou dun lment de structure avec comportement non linaire local,
intgr par couches, NUME_COUCHE est la valeur entire comprise entre 1 et le nombre de
couches, ncessaire pour prciser la couche o l'on dsire effectuer le calcul lmentaire.
Par convention, la couche 1 est la couche infrieure (dans le sens de la normale) dans le cas
des lments de coque.
Pour la couche nume dfinie par NUME_COUCHE, permet de prciser l'ordonne o l'on dsire
effectuer le calcul lmentaire : INF / MOY / SUP correspondent aux point d'intgration
situs en peau interne / moyenne / externe de la couche.
Le mot cl facteur COQU_EXCENT permet de modifier le plan de calcul des efforts gnraliss
( options EFGE_ELNO et EFGE_ELGA) pour un modle avec des lments de plaques (DKT, DST,
Q4G, DKTG) en tenant compte de lexcentrement (MODI_PLAN='OUI').
Les options de post-traitement disponibles sont :
OPTIONS COQUE_3D DKT DST,
Q4G
DKTG,
Q4GG
COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE_* SHB8
'ENEL_ELGA'
'ENEL_ELNO'
'ENEL_ELEM'

'ENER_ELAS'
'EPSI_ELGA'
'SIEQ_ELGA'
DEGE_ELGA
DEGE_ELNO
ECIN_ELEM
EFGE_ELGA
EFGE_ELNO
EPOT_ELEM
EPSI_ELNO
SIEQ_ELNO
SIEF_ELGA
SIEF_ELNO
SIGM_ELNO
VARI_ELNO
SIEF_ELGA : Calcul des efforts gnraliss par lment aux points dintgration de llment
partir des dplacements (utilisation uniquement en lasticit). Repre utilisateur.
SIGM_ELNO : Calcul des contraintes par lment aux nuds. Repre utilisateur.
SIEQ_ELNO : Contraintes quivalentes aux nuds, calcules en un point de lpaisseur
partir de SIGM_ELNO :
VMIS : Contraintes de Von Mises.
VMIS_SG : Contraintes de Von Mises signes par la trace des contraintes.
PRIN_1, PRIN_2, PRIN_3 : Contraintes principales.
Ces contraintes sont indpendantes du repre.
EFGE_ELGA : Calcul des efforts gnraliss par lment aux points dintgration de
llment partir des dplacements (utilisation uniquement en lasticit). Repre utilisateur.
EFGE_ELNO : Calcul des efforts gnraliss par lment aux nuds partir des
dplacements (utilisation uniquement en lasticit). Repre utilisateur.
EPSI_ELNO : Calcul des dformations par lment aux nuds partir des dplacements, en
un point de lpaisseur (utilisation uniquement en lasticit). Repre utilisateur.
Code_Aster
Version
default
EPSI_ELGA : Calcul des dformations par lment aux points d'intgration partir des
dplacements, en un point de lpaisseur (utilisation uniquement en lasticit). Repre
intrinsque.
DEGE_ELGA : Calcul des dformations gnralises par lment aux points dintgration de
llment partir des dplacements. Repre utilisateur.
DEGE_ELNO : Calcul des dformations gnralises par lment aux nuds partir des
dplacements. Repre utilisateur.
EPOT_ELEM : Calcul de lnergie lastique linaire de dformation par lment partir des
dplacements.
ENER_TOTALE : calcul de l'nergie de dformation totale intgre sur l'lment
ENER_ELAS : calcul de l'nergie de dformation lastique intgre sur l'lment
ENEL_ELGA / ENEL_ELNO : nergie lastique aux points d'intgration ou aux nuds
ENEL_ELEM : nergie lastique sur l'lment
ECIN_ELEM : Calcul de lnergie cintique par lment.
EFGE_ELNO : Option dactivation du calcul du tenseur des efforts gnraliss (voir
paragraphe []) par lment aux nuds, dans le repre utilisateur, par intgration des
contraintes SIEF_ELGA (en non-linaire).
EFGE_ELGA : Option dactivation du calcul du tenseur des efforts gnraliss (voir
paragraphe [ ]) par lment aux points de gauss de l'lment, dans le repre utilisateur, par
intgration des contraintes SIEF_ELGA (en non-linaire).
VARI_ELNO : Option dactivation du calcul du champ de variables internes par lment et par
couche aux nuds. Pour chaque point dintgration surfacique, on stocke les informations
sur les couches en commenant par la premire, niveau INF. Le nombre de variables
reprsentes vaut donc 3*NCOU*NBVARI o NBVARI reprsente le nombre de variables
internes.
NUME_COUCHE : Dans le cas dun matriau multicouche (composite ou coque en plasticit),
valeur entire comprise entre 1 et le nombre de couches, ncessaire pour prciser la couche
o lon veut effectuer le calcul lmentaire.
NIVE_COUCHE : Pour la couche n , on peut prciser lordonne o lon dsire effectuer le
calcul lmentaire. Un calcul en peau interne est indiqu par INF, en peau externe par
SUP et sur le feuillet moyen par MOY (suivant le sens de la normale).
2.4.3 Calculs aux nuds : oprateur CALC_CHAMP
OPTIONS COQUE_3D DKT DST Q4G COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE SHB8
FORC_NODA
REAC_NODA
_NOEU
Pour les lments de plaques et coques, loprateur CALC_CHAMP (documentation [U4.81.04]) permet
uniquement le calcul des forces et ractions (calcul des champs aux nuds par moyennation, option
_NOEU) :
partir des contraintes, l'quilibre : FORC_NODA (calcul des forces nodales partir des
contraintes aux points dintgration, lment par lment),
puis en enlevant le chargement appliqu : REAC_NODA (calcul des forces nodales de raction
aux nuds, partir des contraintes aux points dintgration, lment par lment) :
REAC_NODA = FORC_NODA - chargements appliqus ,
utiles pour vrification du chargement et pour calculs de rsultantes, moments, etc.
2.4.4 Calculs de quantits sur tout ou partie de la structure : oprateur POST_ELEM
Loprateur POST_ELEM (documentation [U4.81.22]) permet de calculer des quantits sur tout ou
partie de la structure. Les quantits calcules correspondent des options de calcul particulires de la
modlisation affecte.
Code_Aster
Version
default
OPTIONS Oprateur COQUE_3D DKT DST Q4G COQUE_C_PLAN
COQUE_D_PLAN
COQUE_AXIS
GRILLE SHB8
MASS_INER POST_ELEM
ENER_POT POST_ELEM
ENER_CIN POST_ELEM
ENER_ELAS POST_ELEM

MASS_INER : calcul des caractristiques gomtriques (volume, centre de gravit, matrice
d'inertie) pour les lments plaques et courbes.
ENER_POT : calcul de lnergie potentielle de dformation l'quilibre partir des
dplacements en mcanique linaire des milieux continus (2D et 3D) et en mcanique
linaire pour les lments de structures, ou bien lnergie dissipe thermiquement
l'quilibre en thermique linaire partir des tempratures (cham_no_TEMP_R).
ENER_CIN : calcul de lnergie cintique partir dun champ de vitesse ou partir dun
champ de dplacement et dune frquence (uniquement pour les lments de structure et les
lments 3D).
ENER_ELAS : calcul de lnergie de dformation lastique.

2.4.5 Valeurs de composantes de champs de grandeurs : oprateur POST_RELEVE_T

Loprateur POST_RELEVE_T (documentation [U4.81.21]) permet, sur un groupe de nuds, dextraire
des valeurs ou effectuer des calculs :
dextraire des valeurs de composantes de champs de grandeurs ;
deffectuer des calculs de moyennes et dinvariants :
Moyennes,
Rsultantes et moments de champs vectoriels,
Invariants de champs tensoriels,
Trace directionnelle de champs,
Dexpression dans les repres GLOBAL, LOCAL, POLAIRE, UTILISATEUR ou CYLINDRIQUE

Le concept produit est de type table.
Pour utiliser POST_RELEVE_T, il est ncessaire de dfinir trois concepts :
un lieu : loption NUD (exemple : N01 N045) ou bien loption GROUP_NO (exemple : APPUI) ;
un objet : au choix, loption RESULTAT (SD rsultat : EVOL_ELAS,...) ou bien loption
CHAM_GD (CHAM_NO : DEPL, ou CHAM_ELEM : SIGM_ELNO, ) ;
une nature : au choix, loption EXTRACTION (valeur, ) ou bien loption MOYENNE
(moyenne, maxi, mini, ).
Si on vient dune interface avec un mailleur), les nuds sont rangs par ordre numrique. Il
faut rordonner les nuds le long de la ligne de dpouillement. La solution est dutiliser
loprateur DEFI_GROUP avec loption NOEU_ORDO . Cette option permet de crer un
GROUP_NO ordonn contenant les nuds dun ensemble de mailles forms de segments (
SEG2 ou SEG3 ).
Un exemple dextraction de composante est donn dans le cas-test SSNL503 (voir la description au
paragraphe [2.5.3] page 39) :
TAB_DRZ=POST_RELEVE_T( ACTION=_F(
GROUP_NO = 'D',
INTITULE = 'TB_DRZ',
RESULTAT = RESUL,
NOM_CHAM = 'DEPL',
NOM_CMP = 'DRZ',
TOUT_ORDRE = 'OUI',
OPERATION = 'EXTRACTION'
)
Code_Aster
Version
default
)
Cette syntaxe a pour but :
dextraire : OPERATION = 'EXTRACTION'
sur le groupe de nuds D :
GROUP_NO = 'D'
la composante DRZ du dplacement :
NOM_CHAM = 'DEPL', NOM_CMP = 'DRZ',
pour tous les instants de calcul : TOUT_ORDRE = 'OUI'

2.4.6 Impression des rsultats : oprateur IMPR_RESU

Loprateur IMPR_RESU permet dcrire le maillage et/ou les rsultats d'un calcul sur listing au format
RESULTAT ou sur un fichier dans un format visualisable par des outils de post-traitement externes
Aster : format RESULTAT et ASTER (documentation [U4.91.01]), format CASTEM (documentation
[U7.05.11]), format IDEAS (documentation [U7.05.01]), format MED (documentation [U7.05.21]) ou
format GMSH (documentation [U7.05.32]).
Cette procdure permet d'crire au choix :
un maillage,
des champs aux nuds (de dplacements, de tempratures, de modes propres, de modes
statiques, ),
des champs par lments aux nuds ou aux points de GAUSS (de contraintes, defforts
gnraliss, de variables internes).
Les lments de plaque et coque tant traits de la mme manire que les autres lments finis, nous
renvoyons le lecteur aux notes dutilisation correspondant au format de sortie quil souhaite utiliser.
2.5 Exemples
Les cas-tests retenus ici sont des cas-tests classiques issus de la littrature et qui sont couramment
utiliss pour valider ce type dlments.
On rappelle que les modlisations DKT correspondent la thorie de Love-Kirchhoff et les
modlisations DST, Q4G la thorie avec nergie de cisaillement transverse (Reissner). Les rsultats
pour la modlisation COQUE_3D ne sont prsents que pour une thorie avec nergie de cisaillement
transverse.

2.5.1 Analyse statique linaire
Code_Aster
Version
default
SSLS20

Titre :: Coque cylindrique pince bords libres
Documentation V : [V3.03.020]
Modlisations :
SSLS20A DKT
SSLS20B COQUE_3D MEC3QU9H
SSLS20C COQUE_3D MEC3TR7H
SSLS100

Titre : Plaque circulaire encastre soumise une pression
uniforme.
Modlisations :
SSLS100K COQUE_3D MEC3QU9H
SSLS100L COQUE_3D MEC3TR7H
SSLS100B DKT
SSLS100E DKQ
SSLS100F DST
SSLS100G DSQ
SSLS100H Q4G
SSLS100I,J COQU_AXIS
SSLS101

Titre : Plaque circulaire pose soumise une pression
uniforme.
Modlisations :
SSLS101J COQUE_3D MEC3QU9H
SSLS101I COQUE_3D MEC3TR7H
SSLS101B DKT
SSLS101E DKQ
SSLS101F DST
SSLS101G DSQ
SSLS101H Q4G
SSLS101C SHB8 SHB8
SSLS101D SHB8 SHB20
SSLS101K SHB8 SHB6
SSLS101L SHB8 SHB15
Code_Aster
Version
default
SSLS104

Titre : Coque cylindrique pince avec diaphragme.
Modlisations :
SSLS104C COQUE_3D MEC3TR7H
SSLS104A DKT
SSLS105

Titre : Hmisphre doublement pinc.
Modlisations :
SSLS105A DKT
SSLS105C SHB8 SHB8
SSLS107

Titre : Panneau cylindrique soumis son propre poids.
Modlisations :
SSLS107A COQUE_3D MEC3QU9H
SSLS107B COQUE_3D MEC3TR7H

SSLS108

Titre : Coque hlicodale sous charges concentres.
Modlisations :
SSLS108A COQUE_3D MEC3QU9H
SSLS108B COQUE_3D MEC3TR7H
SSLS108C SHB8 SHB8
SSLS108D SHB8 SHB8
SSLS108E SHB8 SHB6
SSLS108F SHB8 SHB6
SSLS108G SHB8 SHB20
SSLS108H SHB8 SHB15
Remarques :
Utilisation dconseille avec des DKT/DKQ, sans
cisaillements transverses.
Dautres cas-tests sont dcrits plus brivement dans le tableau suivant :
Nom Modlisation Remarques
Code_Aster
Version
default
hpla100a

hpla100b
hpla100c
hpla100d

hpla100e
hpla100f

2D_AXIS

COQUE_AXIS
COQUE_3D
COQUE_3D

COQUE
COQUE
Titre : Cylindre creux thermolastique pesant en rotation uniforme.

Ce test a pour but dprouver les seconds membres correspondants
aux effets de pesanteur et une acclration due une rotation
uniforme.
Les solutions analytiques pour les COQUE_3D incluent la variation de
mtrique dans lpaisseur de la coque. Les solutions analytiques pour
les plaques sont sans correction de mtrique

hsls01a

hsls01b
DKT/DST/Q4G

COQUE_3D
Titre : Plaque mince encastre soumise un gradient thermique
dans lpaisseur.

hsns100a
hsns100b
COQUE_3D/DKT
COQUE_3D/DKT
Titre : Plaque soumise un gradient de temprature dans
lpaisseur.
Ce cas-test permet de tester deux faons dimposer le champ
thermique. Les rsultats obtenus en a et b doivent tre identiques,
mais les solutions de rfrence obtenues sont numriques.

ssll102a COQUE_C_PLAN Titre : Poutre encastre soumise des efforts unitaires
ssls501a

COQUE_D_PLAN Titre : Cylindre infiniment long soumis deux lignes de charges.

ssls114a
ssls114b
ssls114c

ssls114d
ssls114i
ssls114j
ssls114k
COQUE_3D
COQUE_3D
DKT/DST

DKQ/DSQ
COQUE_AXI
COQUE_C_PLAN
COQUE_D_PLAN
Titre : Mise sous pression dun quart de virole cylindrique.
Solution de rfrence analytique. Permet de tester le terme de
pression et lorientation des normales. On teste les rsultats en
dplacement radial et en contraintes radiales.
ssls124a
Ssls124b
SHB8
SHB8
Titre : Poutre en flexion avec divers lancements.
Solution de rfrence analytique. Ce test permet de monter les limites
des lments en terme dlancement, dune part, et de montrer leur
bonne convergence pour un maillage trs irrgulier, dautre part.
ssls123a SHB8 Titre : Sphre sous pression externe uniforme.
Code_Aster
Version
default
Solution de rfrence analytique. Ce test permet d'valuer la qualit
de la modlisation des forces de pression..

2.5.2 Analyse modale en dynamique
sdls01a
sdls01b

sdls01c
sdls01d
sdls01e
sdls01f

sdls01g
sdls01h

DKT
DKT

DKT
DKT
COQUE_3D
COQUE_3D

COQUE_3D
COQUE_3D
Titre : Plaque carre mince libre ou encastre sur un bord
Il sagit dun calcul modal et dun calcul de rponse harmonique.
Pour le calcul modal, il sagit de calculer les modes propres de
flexion dune plaque carre mince libre ou encastre sur un bord.

a - Artes de la plaque orientes selon les axes du repre.
b - Orientation quelconque de la plaque et rponse harmonique
pour la plaque encastre.
c - Calcul modal par sous-structuration dynamique classique et
cyclique.
d - Calcul modal suite une sous-structuration de Guyan.
e - Artes de la plaque orientes selon les axes du repre.
f - Artes de la plaque orientes selon les axes du repre.
g - Orientation quelconque de la plaque et rponse harmonique
h - Orientation quelconque de la plaque et rponse harmonique
Pour a et b la prcision sur les frquences propres est
infrieure 1% jusquau sixime mode de flexion
Pour c en sous-structuration, la qualit des rsultats peut
tre amliore par lutilisation dun maillage de sous-
structure plus fin.
Pour d, il est ncessaire afin dobtenir une prcision de 1%
sur les frquences propres de condenser aussi sur les
nuds milieu des bords.
Pour e, f, g et h, la prcision sur les frquences propres est
infrieure 1% jusquau sixime mode de flexion pour
les lments quadrangle et infrieure 2% pour les
lments triangles.

Llment de coque MEC3QU9H est performant compar llment
DKT qui est lui-mme plus performant que llment MEC3TR7H.
Code_Aster
Version
default
sdls109a
sdls109b,c
sdls109d,e
sdls109f
sdls109h
DKQ ( MEDKQU4 )
et
DSQ ( MEDSQU4 )
DKT ( MEDKTR3 ) et
DST ( MEDSTR3 )
COQUE_3D
( MEC3QU9H et
MEC3TR7H )
COQUE_C_PLAN
( METCSE3 )
SHB8

Titre : Frquences propres dun anneau cylindrique pais.
Ce test est inspir dune tude vibratoire ralise sur le collecteur
VVP des tranches N4. Ce collecteur est pais et prsente un rapport
maximum paisseur sur rayon moyen de 0,13. Cette valeur,
pouvant tre typique dune structure industrielle, est lgrement
suprieure la valeur limite de validit habituellement reconnue
pour les plaques et coques. Dans cette tude, la modlisation du
collecteur en coques est alors value par comparaison avec un
modle volumique sur un anneau.
Ce test permet dvaluer lalgorithme de recherche de valeurs
propres MODE_ITER_SIMULT [U4.52.03] avec les oprateurs de
rigidit et de masse.

2.5.3 Analyse statique non linaire matriau

SSNL501

Titre : Poutre encastre soumise une pression
uniforme.
Modlisations :
SSNL501E COQUE_3D MEC3QU9H
SSNL501D COQUE_3D MEC3TR7H
SSNL501B DKT
SSNL501C DKQ
SSNL501A COQUE_C_PLAN
ssnp15a
ssnp15b
ssnp15c

ssnp15d
3D
C_PLAN
DKT

COQUE_3D
Titre : Plaque carre en traction-cisaillement - Von Miss
(crouissage isotrope).
Une plaque, constitue dun matriau plastique crouissage
isotrope linaire, est soumise un effort de traction et un effort de
cisaillement. Mme si le test valide plutt la loi de comportement que
les lments sur lesquels elle sapplique, il permet de tester les
valeurs des contraintes, des efforts et des dformations dans le
repre dfini par lutilisateur (ANGL_REP).
ssnv115a

ssnv115b

D_PLAN

DKT

Titre : Tle ondule en comportement non linaire.
Ce test valide les comportements non linaires dans les
Code_Aster
Version
default
ssnv115c

ssnv115d

ssnv115e
ssnv115f
DKT

COQUE_3D

COQUE_3D
COQUE_D_PLAN
modlisations de plaques ou de coques minces. La modlisation A
(2D D_PLAN) sert de rfrence. Les valeurs des dplacements sont
testes.

La modlisation COQUE_D_PLAN fait apparatre des carts sur les
dplacements latraux de la tle de lordre de 13% par rapport aux
autres modlisations. Ceci est d la mthode dintgration dans
lpaisseur qui ne fait intervenir que 5 points de Gauss pour cette
modlisation, compars 19 points pour les lments DKT et DKQ et
8 points pour la modlisation D_PLAN.
2.5.4 Analyse statique non linaire gomtrique
SSNV138

Titre : Plaque cantilever en grandes rotations soumise
un moment.
Modlisations :
SSNV138 COQUE_3D MEC3QU9H
SSNV138 COQUE_3D MEC3TR7H
Remarque :
La plus grande rotation atteinte est lgrement infrieure
. Les rsultats obtenus sont trs satisfaisants, l'cart
maximum est infrieur 0.01%. Il est ncessaire
d'augmenter la valeur de COEF_RIGI_DRZ (10E-5 par
dfaut) 0.001 de faon pouvoir augmenter la valeur de
l'angle de rotation que l'on peut atteindre.
SSNV139

Titre : Plaque biaise.
Modlisations :
SSNV139 COQUE_3D MEC3QU9H
SSNV139 COQUE_3D MEC3TR7H
SSNL502

Titre : Poutre en flambement.
Modlisations :
SSNL502 COQUE_3D MEC3QU9H
SSNL502 COQUE_3D MEC3TR7H

SSNS501

Titre : Grands dplacements dun panneau cylindrique.
Modlisations :
SSNS501 COQUE_3D MEC3QU9H
SSNS501 COQUE_3D MEC3TR7H
2L
F
2

Code_Aster
Version
default

Titre : Claquage dun panneau cylindrique sous force
ponctuelle.
Modlisations :
SSNS101a,b,c,d SHB8 SHB8
SSNS101e SHB8 SHB20
SSNS101f SHB8 SHB6
SSNS101g SHB8 SHB15
Ce test de mcanique quasi - statique non linaire permet
de valider les lments SHB en non linaire gomtrique
et matriau.

Titre : Flambage dune coque cylindrique avec raidisseur.
Modlisations :
SSNS102a SHB8 SHB8
SSNS102b SHB8 SHB20
Ce test de mcanique quasi - statique non linaire permet
de valider les lments SHB8 en non linaire
gomtrique, avec ou sans prise en compte des pressions
suiveuses et en flambage dEuler. Il montre les capacits
de cet lment traiter des problmes de coques minces
avec raidisseur.
ssnv140a
ssnv140b
COQUE_3D
COQUE_3D
Titre : Panneau cylindrique encastr soumis une force surfacique.
Cette force est constante pour la modlisation a et suiveuse dans la
modlisation b. Le but de ce cas-test est de vrifier la modlisation
COQUE_3D non-linaire gomtrique en utilisant l'algorithme de mise
jour des grande rotations 3D GROT_GDEP de STAT_NON_LINE et de
vrifier le traitement des pressions suiveuses. Les donnes de ce
problme correspondent une coque mince h/ L=0.625% ce qui est
svre pour l'lment fini triangle MECQTR7H (cas de blocage au
cisaillement transverse).
ssnv141a COQUE_3D Titre : Calotte sphrique pince.
Les donnes de ce problme correspondent une coque mince
h/ L=0.4% ce qui est svre pour l'lment fini triangle MECQTR7H
SU
SUINF
SU90
z
y
x
Code_Aster
Version
default
(cas de blocage au cisaillement transverse). Il est ncessaire
d'augmenter la valeur du COEF_RIGI_DRZ qui attribue une rigidit
autour de la normale des lments de coque qui vaut par dfaut
10
5

la plus petite rigidit de flexion autour des directions dans le plan de la
coque de faon pouvoir augmenter la valeur de l'angle de rotation que
l'on peut atteindre. Des valeurs de ce coefficient jusqu'
10
3
restent
licites.
ssnv144a COQUE_3D Titre : Coude en flexion plane, lastique, encastr sur un cot et soumis
une force linique quivalente un moment de flexion.
Le but de ce cas-test est de vrifier que, pour les lments COQUE_3D,
les solutions quasi-statiques en linaire gomtrique (VMIS_ISOT_LINE
dans STAT_NON_LINE) et en non linaire gomtrique ( GROT_GDEP
dans STAT_NON_LINE) sont proches de la solution statique linaire
(MECA_STATIQUE) dans le domaine des petites perturbations.
ssnv145a
ssnv145b
COQUE_3D
COQUE_3D
Titre : Plaque cantilever en grandes rotations soumise a une pression
suiveuse.
Le but de ce cas-test est de vrifier la modlisation COQUE_3D (maille
TRIA7, QUAD9) en prsence d'une pression de type suiveuse.
2.5.5 Analyse en flambement dEuler
SSLS110

Titre : Stabilit dune plaque carre comprime.
Modlisations :
SSLS110 COQUE_3D MEC3QU9H
SSLS110 COQUE_3D MEC3TR7H
= z, w
y, v
x, u
L/2
Cylindre libre aux extrmits
B O
A
z
x
y
O C
D

Titre : Flambement dun cylindre libre sous pression
externe.
Modlisations :
SSLS125a SHB8 SHB8
SSLS125b SHB8 SHB6
SSLS125c SHB8 SHB20
SSLS125d SHB8 SHB15
Code_Aster
Version
default
SDLS504

Titre : Flambement latral dune poutre (dversement).
Modlisations :
SDLS504 COQUE_3D MEC3QU9H
SDLS504 COQUE_3D MEC3TR7H
SDLS505

Titre : Flambement dune enveloppe cylindrique sous
pression externe.
Modlisations :
SDLS505 COQUE_3D MEC3QU9H
SDLS505 COQUE_3D MEC3TR7H
2.5.6 Raccords coques et autres lments mcaniques
SSLX100

Titre : Mlange 3D-Coque-Poutre en flexion.
Modlisations :
SSLX100A 3D 1 MECA_HEXA20
DKT 4 MEDKTR3
POU_D_E 2 POU_D_E
SSLX100B 3D 1 MECA_HEXA20
DKT 4 MEDKTR3
POU_D_E 2 POU_D_E
SSLX100C 3D 1 MECA_HEXA20
DKT 4 MEDKTR3
COQUE_C_PLAN 1 MECPSE3
On teste les flches, contraintes, dformations axiales
et moments flchissants en 4 points de laxe de la
poutre.
SSLX102

Titre : Tuyauterie coude en flexion.
Modlisations :
SSLX102A DKT et TUYAU
(liaison COQUE_TUYAU)
SSLX102D COQUE et POUTRE
Code_Aster
Version
default
SSLX101A

Titre : Tuyau droit modlis en coques et en poutres
[V3.05.101].
Modlisations :
SSLX101A DIS_TR POI1
DKT MEDKQU4
POU_D_E 2 SEG2
Encastrement de la coque sur le bord P1 . Flexion et
traction en x1 . cart de 3 5% sur les
dplacements et rotations en P2 avec la solution
analytique, d au maillage coque avec des lments
plans.
SSLX101B DKT MEDKQU4, METUSEG3
TUYAU METUSEG3, MEDKQU4
DIS_TR POI1
Cette modlisation a pour objectif de tester la liaison
coque tuyau en prsence de chargements unitaires :
traction, flexion et de torsion. La solution de rfrence
est analytique (RDM). L'cart avec la solution
numrique s'explique par le fait que le maillage en
coques est en ralit constitu d'lments plans
(facettes). La gomtrie du tuyau est donc elle-mme
approche.
SSLX102A DKT MEDKQU4, METUSEG3
TUYAU MEDKQU4, METUSEG3
La modlisation A fait intervenir la liaison
coque_tuyau, la solution obtenue ( 2.7% d'cart en
flexion plane, et 0.4% en flexion hors plan, par
rapport la rfrence: maillage tout coques
(modlisation D) permet de tester le bon
fonctionnement du raccordement entre coque et
tuyau.
Code_Aster
Version
default
3 Thermique
Pour la rsolution de problmes thermomcaniques chans, on doit utiliser pour le calcul thermique
des lments finis de coque thermique [R3.11.01]. Ces lments sont des lments plaques, ou
liniques dans le cas de structures de rvolution ou de structures invariantes suivant un axe. La
courbure de la structure nest pas prise en compte en elle-mme. La mtrique du plan tangent de
chaque lment est calcule en supposant que tous les sommets sont coplanaires. Ces lments
supposent une rpartition a priori parabolique de temprature dans lpaisseur, qui rsulte dun
dveloppement asymptotique en thermique linaire pour une paisseur de coque faible, lorsque les
variations de tempratures ne sont pas trop importantes. Il est noter quun modle bas sur un
dveloppement du champ de temprature plus riche dans lpaisseur voit ses termes dordre
suprieur deux converger vers zro lorsque la coque est mince. On ne peut donc pas traiter les
problmes de chocs thermiques avec forte variation du profil de temprature dans lpaisseur avec
ces coques. Les modalits dutilisation de ces lments sont prsentes en [U1.22.01].
3.1 Dfinition du problme
3.1.1 Discrtisation spatiale et affectation dune modlisation : oprateur
AFFE_MODELE
3.1.1.1 Degrs de liberts
Les degrs de libert sont les tempratures TEMP_MIL (temprature sur la surface moyenne de la
coque), TEMP_INF (temprature sur la surface infrieure de la coque), et TEMP_SUP (temprature sur
la surface suprieure de la coque).
3.1.1.2 Mailles support des matrices de rigidit
Modlisation Maille Nature de la maille lment fini Remarques
COQUE QUAD9
QUAD8
QUAD4
TRIA7
TRIA6
TRIA3
plane
plane
plane
plane
plane
plane
THCOQU9
THCOQU8
THCOQU4
THCOTR7
THCOTR6
THCOTR3
nuds 3 coordonnes x , y , z
COQUE_PLAN SEG3 non suppose plane THCPSE3
nuds 2 coordonnes x , y
COQUE_AXIS SEG3 non suppose plane THCASE3
nuds 2 coordonnes x , y
Pour les THCOTRi, seuls les trois sommets sont exploits pour dfinir la gomtrie locale (plan
tangent, normale). Pour les THCOQUi, on considre que llment est plan et son plan tangent est
dfini par dfaut par 3 des 4 sommets de llment.
Code_Aster
Version
default
3.1.1.3 Mailles support des chargements
Modlisation Maille Elment fini Remarques
COQUE SEG2 THCOSE2 avec TRIA3 et QUAD4
COQUE SEG3 THCOSE3 avec TRIA6,7 et QUAD8 ,9
Tous les chargements applicables aux facettes des lments de coque sont traits par discrtisation
directe sur la maille support de llment en formulation temprature. Aucune maille de chargement
nest donc ncessaire pour les faces des lments de coques.
Pour les chargements applicables sur les bords des lments de coque, une maille support de type
SEG2 (lment THCOSE2) ou SEG3 (lment THCOSE3) doit tre utilise.
Pour les tempratures imposes les mailles support sont des mailles rduites un point.
3.1.1.4 Modle : AFFE_MODELE
Laffectation de la modlisation passe par loprateur AFFE_MODELE [U4.41.01].
AFFE_MODELE Remarques
AFFE
PHENOMENE : 'THERMIQUE'
MODELISATION 'COQUE'
'COQUE_PLAN'
'COQUE_AXIS'
3.1.2 Caractristiques lmentaires : AFFE_CARA_ELEM
Dans cette partie, les oprandes caractristiques des lments de plaques et de coques en thermique
sont dcrites. La documentation dutilisation de loprateur AFFE_CARA_ELEM est [U4.42.01].
AFFE_CARA_ELEM COQUE COQUE_PLAN COQUE_AXIS Remarques
COQUE
EPAIS
Les caractristiques affectes aux matriaux sont les mmes que pour un calcul mcanique. Il est
noter quil nest pas utile de dfinir un repre particulier pour lexploitation des rsultats du calcul
thermique car ceux-ci se limitent aux champs de temprature, grandeur scalaire, indpendante du
rfrentiel utilis.
3.1.3 Matriaux : DEFI_MATERIAU
DEFI_MATERIAU COQUE COQUE_PLAN COQUE_AXIS Remarques
THER
THER_FO
Les matriaux utiliss avec des lments plaques ou coques en thermique peuvent avoir des
caractristiques thermiques linaires constantes ou dpendantes de lincrment de chargement.
Code_Aster
Version
default
3.1.4 Chargements et conditions limites : AFFE_CHAR_THER et AFFE_CHAR_THER_F
Laffectation des chargements et des conditions aux limites sur un modle thermique seffectue
laide de loprateur AFFE_CHAR_THER, si les chargements et les conditions aux limites mcaniques
sur un systme sont des valeurs relles ne dpendant daucun paramtre, ou bien
AFFE_CHAR_THER_F, si ces valeurs sont fonctions de la position ou de lincrment de chargement.
La documentation dutilisation dAFFE_CHAR_THER et AFFE_CHAR_THER_F est [U4.44.02].
3.1.4.1 Liste des mot-cls facteur d AFFE_CHAR_THER
Les valeurs des chargements affectes sont relles et ne dpendent daucun paramtre.
AFFE_CHAR_THER
gnraux
COQUE COQUE_PLAN COQUE_AXIS Remarques
TEMP_IMPO
AFFE_CHAR_THER
particuliers
FLUX_REP sur les faces et les bords
des lments surfaciques
ECHANGE sur les faces et les bords
des lments surfaciques
TEMP_IMPO : Mot cl facteur utilisable pour imposer, sur des nuds ou des groupes de
nuds, une temprature.
FLUX_REP : Mot cl facteur utilisable pour appliquer des flux normaux une face de coque
thermique dfinie par une ou plusieurs mailles ou des groupes de mailles de type triangle ou
quadrangle.
ECHANGE : Mot cl facteur utilisable pour appliquer des conditions dchange avec une
temprature extrieure une face de coque, dfinie par une ou plusieurs mailles ou groupes
de mailles de type triangle ou quadrangle.
3.1.4.2 Liste des mot-cls facteur d AFFE_CHAR_THER_F

Les valeurs des chargements affectes peuvent tre fonction des coordonnes globales et du temps,
ou de la temprature en thermique non linaire (sauf en coques).
AFFE_CHAR_THER_F
gnraux
TEMP_IMPO
AFFE_CHAR_THER_F
particuliers
FLUX_REP sur les faces et les bords des
lments surfaciques
ECHANGE sur les faces et les bords des
lments surfaciques
Code_Aster
Version
default
3.2 Rsolution
3.2.1 Calculs transitoires : oprateur THER_LINEAIRE
Option de calcul
transitoire
CHAR_THER_EVOL
Il sagit ici du traitement des problmes de thermique volutive.
3.3 Calculs supplmentaires et post-traitements
3.3.1 Calculs en post-traitement
On prsente ci-aprs les options de post-traitement pour les lments de plaques et de coques
OPTIONS
lmentaires
FLUX_ELNO
FLUX_ELGA
FLUX_ELNO : Cette option effectue le calcul de flux de chaleur aux nuds partir de la
temprature.
FLUX_ELGA : Cette option effectue le calcul de flux de chaleur aux points dintgration
partir de la temprature.
3.4 Exemples
On donne ici la liste des cas-tests disponibles pour les coques thermiques. Ce sont des cas-tests de
thermique stationnaire. Les rsultats sont corrects pour lensemble de ces cas-tests, quel que soit
llment utilis.
Nom Modlisation lment Remarques
tplp301a COQUE THCOTR3 Titre : Plaque avec temprature impose
rpartie sinusodalement sur un ct.
Documentation : [V4.05.301]
tplp302a COQUE THCOTR6 Titre : Plaque rectangulaire avec temprature
impose sur les cts.

tpls100a
tpls100b
COQUE
COQUE_PLAN
THCOTR6/THCOTR3
THCPSE3
Titre : Plaque infinie soumise un couple de flux
de chaleur antisymtriques stationnaires sur ses
deux demi-faces.

La conduction est linaire, homogne et isotrope.
tpls101a
tpls101b
tpls101c
tpls101d
tpls101e
COQUE THCOTR6/THCOSE3
THCOQU4/THCOSE2
THCOQU8/THCOSE3
THCOQU9/THCOSE3
THCOTR7/THCOSE3
Titre : Plaque infinie soumise un couple de
conditions thermiques avec lextrieur,
symtriques par rapport au feuillet moyen.
Code_Aster
Version
default
La conduction est linaire, homogne et isotrope.
tpls302a
tpls302b
tpls302c
tpls302d
COQUE THCOQU8/THCOSE3
THCOQU4/THCOSE2
THCOQU9/THCOSE3
THCOTR7/THCOSE3
Titre : Plaque rectangulaire avec convection et
temprature impose
4 Chanage thermomcanique
4.1 Formalisme
Pour la rsolution de problmes thermomcaniques chans, on doit utiliser pour le calcul thermique
des lments finis de coque thermique [R3.11.01] dont le champ de temprature est rcupr comme
donne dentre du Code _Aster pour le calcul mcanique. Il faut donc quil y ait compatibilit entre le
champ thermique donn par les coques thermiques et celui rcupr par les coques mcaniques. Ce
dernier est dfini par la connaissance des 3 champs TEMP_INF, TEMP_MIL et TEMP_SUP donns en
peaux infrieure, milieu et suprieure de coque. Le tableau ci-dessous indique ces compatibilits :
Modlisation
THERMIQUE
Maille lment Maille lment Modlisation
MECANIQUE
COQUE QUAD9 THCOQU9 QUAD9 MEC3QU9H COQUE_3D
COQUE QUAD8 THCOQU8
COQUE QUAD4 THCOQU4 QUAD4 MEDKQU4
MEDSQU4
MEQ4QU4
DKT
DST
Q4G
COQUE TRIA7 THCOTR7 TRIA7 MEC3TR7H COQUE_3D
COQUE TRIA6 THCOTR6
COQUE TRIA3 THCOTR3 TRIA3 MEDKTR3
MEDSTR3
DKT
DST
COQUE_PLAN SEG3 THCPSE3 SEG3
SEG3
METDSE3
METCSE3
COQUE_D_PLAN
COQUE_C_PLAN
COQUE_AXIS SEG3 THCASE3 SEG3 MECXSE3 COQUE_AXIS
Remarques :
Les nuds des lments de coques thermiques et de plaques ou de coques
mcaniques doivent se correspondre. Les maillages pour la thermique et la mcanique
auront donc le mme nombre et le mme type de mailles.
Les lments de coques thermiques surfaciques sont traits comme des lments
plans par projection de la gomtrie initiale sur le plan dfini par les 3 premiers
sommets. Pour le chanage de calculs avec des lments courbes mcaniques il est
donc ncessaire que la gomtrie de la plaque ne soit pas trop loigne de celle de la
coque. Lorsque la structure est courbe, cela ncessite donc pour le calcul thermique de
la mailler de faon suffisamment fine de faon avoir des rsultats corrects en
prvision de la partie mcanique. Seuls les lments liniques de thermique sont
parfaitement associs aux lments liniques correspondants en mcanique car
prenant en compte la courbure de la structure maille.
Le chanage avec des matriaux multicouches nest pas disponible pour linstant.
Le chanage thermomcanique est aussi possible si lon connat, analytiquement ou par
des mesures exprimentales, la variation du champ de temprature dans lpaisseur de
la structure ou de certaines parties de la structure. Dans ce cas on travaille avec une
carte de temprature dfinie a priori ; le champ de temprature nest plus donn par les
trois valeurs TEMP_INF , TEMP_MIL et TEMP_SUP du calcul thermique obtenues par
EVOL_THER . Loprateur DEFI_NAPPE permet de crer de tels profils de tempratures
partir des donnes fournies par lutilisateur. Ces profils sont affects par la commande
CREA_CHAMP et CREA_RESU (cf. le cas-test hsns100b ). On notera quil nest pas
ncessaire pour le calcul mcanique que le nombre de points dintgration dans
lpaisseur soit gal au nombre de points de discrtisation du champ de temprature
Code_Aster
Version
default
dans lpaisseur. Le champ de temprature est automatiquement interpol aux points
dintgration dans lpaisseur des lments de plaques ou de coques.
L'volution thermique que l'on peut associer au champ de matriau par
AFFE_MATERIAU/AFFE_VARC doit tre prte tre utilise par les lments finis du
modle mcanique. Un problme se pose pour les lments de type coque ou tuyau qui
utilisent une temprature variant dans l'paisseur sur les diffrentes couches. Pour ces
lments, il est ncessaire de prparer le calcul de la temprature sur les couches en
amont de la commande AFFE_MATERIAU. Pour cela, l'utilisateur doit utiliser la
commande CREA_RESU avec l'une des oprations PREP_VRC1 ou PREP_VRC2 des
VaRiables de Commande") :
- OPERATION='PREP_VRC1' : calcul de la temprature dans les couches d'une coque
en partant d'une temprature TEMP= f(EPAIS,INST)
- OPERATION='PREP_VRC2' : calcul de la temprature dans les couches d'une coque
en partant d'une temprature calcule par Code_Aster avec un modle de coques
(TEMP_MIL/TEMP_INF/TEMP_SUP).
Code_Aster
Version
default
HPLA100

Il sagit dtudier un phnomne de dilatation thermique
o les champs de temprature sont calculs avec
THER_LINEAIRE par un calcul stationnaire :
- dilatation thermique :

T (j)T
ref
(j)=0.5(T
s
+T
i
)+2.(T
s
+Ti )(rR)/ h
avec :

T
s
=0.5 C, T
i
=0.5C , T
ref
=0C

T
s
=0.1C ,T
i
=0.1C , T
ref
=0 C

On teste les contraintes, les efforts et moments
flchissants en L et M . Les rsultats de rfrence
sont analytiques. Pour les modlisations COQUE_3D on
prend en compte la variation de mtrique avec
lpaisseur de la coque. Trs bons rsultats quel que soit
le type dlment considr.
Titre : Cylindre creux thermolastique pesant en
rotation uniforme
Modlisations :
HPLA100A
Thermique PLAN 32 THPLQU8
Mcanique AXIS 32 MEAXQU8
HPLA100B
Thermique COQUE_PLAN 10 THCPSE3
Mcanique COQUE_AXIS 10 MECXSE3
HPLA100C
Thermique COQUE 32 THCOQU9
Mcanique COQUE_3D 32 MEC3QU9H
HPLA100D
Thermique COQUE 64 THCOTR7
Mcanique COQUE_3D 64 MEC3TR7H
HPLA100E
Thermique COQUE 200 THCOQU4
Mcanique COQUE 200 MEDKQU4
HPLA100F
Thermique COQUE 400 THCOTR3
Mcanique COQUE 400 MEDKTR3

5 Conclusion et conseils dutilisation
Dans le tableau suivant, un rcapitulatif des possibilits offertes par les modlisations plaques et
coques est dcrit.
Modlisation DKT DST,
Q4G
DKTG,
Q4GG
COQUE_3
D
COQUE_AXIS,
COQUE_D_PLAN
COQUE_C_PLAN
SHB8 GRILLE_*
Domaine d'application
Statique linaire :
Matriau isotrope
X X X X X X X
Matriau orthotrope,
composite
X X
Statique non-linaire X X X X X X
Code_Aster
Version
default
matriau
Statique non-linaire
gomtrique
X X
Analyse dynamique X X X X X X X
Flambement d'Euler X X X
Sur la figure ci-dessous le domaine d'utilisation des plaques et des coques est schmatis.

h/L
0.1 (1/10) 0.05 (1/20)
DKT, DKQ, SHB8 DST, DSQ , Q4G, SHB8
COQUE_3D, COQUE_AXIS, SHB8
COQUE_C_PLAN, COQUE_D_PLAN
lments plans

lments courbes

Coques minces

Coques paisses

Figure 5-a : Domaines dutilisation des plaques, des coques et des coques SHB
Quelques recommandations concernant le domaine dutilisation de ces lments :
Structures minces : pour ces structures, dont le rapport h/L est infrieur 1/20, les effets de
cisaillement transverse peuvent tre ngligs et la thorie de Kirchhoff sapplique. On conseille
dutiliser pour ce type de structure des lments plaques DKT-DKQ ou des lments de coque
courbe (COQUE_3D, _AXIS, _C_PLAN, _D_PLAN) ou des lments coques volumiques SHB8,
SHB15, SHB20. Il est conseill dutiliser de prfrence les lments DKT et DKQ qui
donnent de trs bons rsultats sur les dplacements et plus approximatifs sur les contraintes (
recommander pour les analyses vibratoires). Mme si on doit utiliser un grand nombre de ces
lments, les temps dexcution restent raisonnables compars ceux des lments courbes.
Llment SHB6 est moins performant que les autres lments SHB pour les structures minces.
Nanmoins, llment SHB6 est indispensable pour mailler des gomtries quelconques. Une
modlisation mixte SHB6 et SHB8 donne des rsultats raisonnables en prcision et en temps de
calcul. Il est donc conseill de mailler la structure au maximum avec des lments SHB8, et
de mailler la partie restante avec des lments SHB6. Il est dconseill dutiliser llment
SHB6 pour mailler la totalit de la structure.

Structures paisses : pour ces structures, on utilisera des lments de plaque DST, DSQ et Q4G
qui tiennent compte du cisaillement transverse avec un facteur de correction de cisaillement
k=5/ 6 (thorie de Reissner) ou de prfrence des lments de coque courbe. Le facteur de
correction du cisaillement permet de passer dune thorie de Hencky-Mindlin-Naghdi pour k=1 ,
une thorie de Reissner pour k=5/ 6 . Le coefficient de cisaillement n'est ajustable que pour
les lments de coque courbe mais il est conseill de ne pas modifier sa valeur par dfaut.
Lorsque la modlisation Q4G a t privilgie, il faut raliser une petite tude de sensibilit
au maillage. Les tests montrent en effet que cette modlisation ncessite un maillage
suffisamment fin dans les directions sollicites en flexion pour obtenir des erreurs faibles.
Comme pour les structures minces, les recommandations dutilisation pour les lments SHB pour
les structures paisses sont les mmes. Les lments SHB8, SHB15 et SHB20 sont performants.
Code_Aster
Version
default
Il est donc conseill dutiliser llment SHB6 pour mailler les gomtries complexes quon ne
peut pas modliser par le seul lment SHB8.

Les lments DKT, DKQ, DST, DSQ et Q4G sont des lments plans, ils ne prennent pas en compte la
courbure des structures, il est donc ncessaire de raffiner le maillage dans le cas o la courbure est
importante si on veut viter les flexions parasites.

La variation de mtrique de la gomtrie (cest--dire son rayon de courbure) en fonction de son
paisseur est prise en compte :

automatiquement pour la modlisation COQUE_3D
dfini par l'utilisateur pour les modlisations COQUE_AXIS, COQUE_D_PLAN et
COQUE_C_PLAN.
Llment mcanique optimal en statique daprs lensemble des cas-tests du paragraphe
[2.5] est llment de coque 9 nuds MEC3QU9H, qui permet dobtenir de bons dplacements et
de bonnes contraintes grce son interpolation P2 en membrane. Cest un lment polyvalent qui
peut tre la fois utilis pour reprsenter des structures trs minces ( h/ L1/100 ) ou plus
paisses. Comme, par ailleurs, llment de coque 7 nuds MEC3TR7H est moins performant, il
est conseill lutilisateur de mailler sa structure en coques avec le plus grand nombre
possible de quadrangles.
Non-linarit matriau : les comportements non linaires (plasticit, etc.) en contraintes planes
sont disponibles pour les lments de coque courbe (COQUE_3D, COQUE_AXIS, COQUE_C_PLAN,
COQUE_D_PLAN) et les lments plaques DKT-DKQ uniquement. Le comportement plastique ne
prend pas les termes de cisaillement transverse qui sont trait de faon lastique, car le
cisaillement transverse est dcoupl du comportement plastique. Pour une bonne reprsentation
de la progression de la plasticit travers l'paisseur, on conseille d'utiliser pour l'intgration
numrique 3 5 couches dans l'paisseur pour un nombre de points de gauss valant
respectivement 3, 5 et 11. Tous les comportements non linaires sont disponibles pour les
lments de coques volumiques SHB (SHB6, SHB8, SHB15 et SHB20).

Non-linarit gomtrique : les non linairits gomtriques (grands dplacements , grandes
rotations) en contraintes planes sont disponibles pour les lments de coque courbe COQUE_3D
uniquement. Les comportements non linaires (grands dplacements, petites rotations) sont
disponibles pour les lments de coques volumiques SHB.
Flambement d'Euler : ce type d'analyse est disponible avec les lments de coque courbe
COQUE_3D et coques volumiques SHB.
Des lments correspondant aux lments mcaniques existent en thermique; les couplages
thermomcaniques sont donc disponibles la fois pour les lments de plaques et de coques. Pour
linstant ces couplages ne sont pas possibles pour des matriaux multicouches.
6 Bibliographie
1 J.L. Batoz, G. Dhatt "Modlisation des structures par lments finis: poutres et plaques"
Herms,Paris(1990)
2 J.L. Batoz, G. Dhatt "Modlisation des structures par lments finis: coques"
Herms,Paris(1992)
3 D. Bui "Evolution de AFFE_CARA_ELEM " CR MMN/97/004.
4 S. Andrieux "Raccords 3D/poutre, 3D/coques et autres fantaisies". Note paratre.
Code_Aster
Version
default
5 E. Lorentz "Grandes dformations plastiques. Modlisation dans Aster par
PETIT_REAC".EDF/DER CRMMN 1536/07.
6 P. Jetteur "Cinmatique non-linaire des coques". Rapport SAMTECH issu du contrat
PP/GC-134/96.
7 J. Argyris, P. Dunne, C. Malejannakis, E. Schelkie "A simple triangular facet shell element
with application to linear and non linear equilibrium and elastic stability problem ".
Comp. Meth. Appl. Mech. Eng.,vol 11, 1977.
8 F. Frey "Lanalyse statique non-linaire des structures par la mthode des lments finis et
son application la construction mtallique ". Thse de Doctorat, Lige, 1978.
9 A. B. Sabir and A. C. Lock "The application of finite elements to the large deflection
geometrically and non linear behavior of cylindrical shells " Variational methods in
Engineering, edited by Brebbia and Tottenham, Southampton, 1972.
10 G.S. Dhatt "Instability of thin shells by the finite elements method ". Proc. IASS Symp. ,Vol1,
Vienna 1970, pp1-36.
11 Raccord 3D-Poutre [R3.03.03].
12 Pression suiveuse pour les lments de coques volumiques [R3.03.07].
13 Coques thermolastiques axisymtriques et 1D [R3.07.02].
14 Elments de plaque DKT, DST, DKQ, DSQ et Q4 [R3.07.03].
15 Elments finis de coques volumiques [R3.07.04].
16 Elments de coques volumiques en non linaire gomtrique[R3.07.05].
17 Modle de thermique pour les coques minces [R3.07.11].
18 Elments finis de tuyau droit et courbe avec ovalisation, gonflement et gauchissement en
lasto-plasticit [R3.08.06].
19 Modle de thermique pour les coques minces [R3.11.01].
20 Intgration des relations lasto-plastiques [R5.03.02].
21 Relation de comportement lastique non linaire [R5.03.20].
22 Modlisation statique et dynamique des poutres en grandes rotations [R5.03.40].
23 Oprateur DEFI_MATERIAU [U4.23.01].
24 Oprateur DEFI_COMPOSITE [U4.23.03].
25 Oprateur AFFE_CARA_ELEM [U4.24.01].
26 Oprateur AFFE_CHAR_MECA et AFFE_CHAR_MECA_F [U4.25.01].
27 Oprateur AFFE_CHAR_THER et AFFE_CHAR_THER_F [U4.25.02].
28 Oprateur STAT_NON_LINE [U4.32.01].
29 Oprateur CALC_MATR_ELEM [U4.41.01].
30 Oprateur CALC_CHAMP [U4.81.04].
Code_Aster
Version
default
31 LEGAY A. and COMBESCURE A. : Elasto-plastic stability analysis of shells using the
physically stabilized finite element SHB8PS, International Journal for Numerical
Methods and Engineering, Vol. 57, 1299-1322, 2003.
32 ABED-MERAIM F. and COMBESCURE A. : SHB8PS a new adaptative assumed strain
continuum mechanics shell element for impact analysis , Computers and
Structures, Vol. 80, 791-803, 2002.
33 BELYTSCHKO T. and BINDEMAN L.P. : Assumed strain stabilization of the eight node
hexahedral elements, Computer Methods in Applied Mechanics and Engineering,
Vol. 105, 225-260, 1993.
34 FLANAGAN D.P. and BELYTSCHKO T. : A uniform strain hexahedron and equilateral with
orthogonal hourglass control, International Journal for Numerical Methods and
Engineering, Vol. 17, 679-706, 1981.
35 RIKS E. : An incremental approach to the solution of snapping and buckling problems,
International Journal of Solids and Structures, Vol. 15, 524-551, 1979.

Code - Aster: Notice D'utilisation Des Éléments Plaques, Coques, Coques Volumiques SHB, Grilles Et Membranes

Transféré par

Droits d'auteur :

Formats disponibles

Code - Aster: Notice D'utilisation Des Éléments Plaques, Coques, Coques Volumiques SHB, Grilles Et Membranes

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Code - Aster: Notice D'utilisation Des Éléments Plaques, Coques, Coques Volumiques SHB, Grilles Et Membranes

Transféré par

Droits d'auteur :

Formats disponibles

Code_Aster

(fournis en degrs) qui dfinissent un vecteur v dont la projection

Vous aimerez peut-être aussi