Version

Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 1/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

Éléments de plaque : Modélisation Q4GG

Résumé

Ce document présente la formulation des éléments de plaque épaisse, utilisés dans la modélisation Q4GG.

Pour cette modélisation, deux types éléments finis sont disponibles, selon les mailles :

• Élément fini quadrangulaire (Q4G)

• Élément fini triangulaire (T3G)

Cette modélisation d’élément de plaque épaisse, comme les modélisations DKT, DKQ, DST, DSQ, DKTG et Q4G
[7], est destinée aux calculs en petites déformations et grands déplacements de structures minces courbes ou
planes. Ce sont des éléments plans qui ne prennent pas en compte la courbure géométrique des structures,
contrairement aux éléments de coque (COQUE_3D) qui sont courbes: il en résulte des flexions parasites qui
peuvent être réduites en utilisant plus d’éléments de façon à pouvoir approcher correctement les géométries
courbes. La formulation en est donc simplifiée et le nombre de degrés de liberté réduit.

Note:
La formulation de l'élément quadrangulaire Q4G utilisée dans la modélisation Q4GG est décrite dans le
document [7]. Par conséquent, on se limite ici au cas de l'élément triangulaire T3G.

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 2/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

Table des Matières

1 Introduction........................................................................................................................................... 3
2 Formulation.......................................................................................................................................... 3
2.1 Géométrie des éléments plaques.................................................................................................. 3
2.2 Cinématique................................................................................................................................... 4
2.3 Support élément fini....................................................................................................................... 4
2.4 Membrane-Flexion......................................................................................................................... 5
2.5 Distorsion transverse..................................................................................................................... 6
2.5.1 Mise en évidence d'un verrouillage en cisaillement transverse...........................................6
2.5.2 Remède aux verrouillages par l'approche T3G....................................................................6
2.6 Calcul des contraintes et des forces internes..............................................................................10
2.6.1 Élasticité linéaire................................................................................................................ 10
2.6.2 Comportement global non-linéaire de type GLRC..............................................................11
2.7 Matrice de masse « lumpée »...................................................................................................... 13
2.8 Second membre correspondant au travail extérieur....................................................................15
3 Implantation des éléments de plaque dans Code_Aster....................................................................16
3.1 Description................................................................................................................................... 16
3.2 Utilisation et développements introduits......................................................................................16
3.3 Calcul en élasticité linéaire.......................................................................................................... 17
3.4 Calcul non-linéaire....................................................................................................................... 17
4 Validation de la modélisation Q4GG.................................................................................................. 18
4.1 SDLS123 – Poutre droite avec amortissement de Rayleigh (V2.03.123)....................................18
4.2 SDNS108 – Réponse dynamique d'une dalle en béton armé (modèle GLRC_DAMAGE)
appuyée sur 4 côtés soumise à une charge concentrée : régime de 4.2 plaque élasto-plastique
(V5.06.108).................................................................................................................................. 18
5 Bibliographie....................................................................................................................................... 19

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 3/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

1 Introduction
On présente ici les détails de la formulation de l'élément de plaque épaisse T3G, telle qu’elle a été
proposée dans la publication [6]. Cet élément est construit d’une manière analogue à l'élément Q4G
[7].

2 Formulation
2.1 Géométrie des éléments plaques
Pour les éléments de plaque on définit une surface de référence, ou surface moyenne, plane (plan
x y par exemple) et une épaisseur h  x , y  . Cette épaisseur doit être petite par rapport aux autres
dimensions (extensions, rayons de courbure) de la structure à modéliser. La figure ci-dessous illustre
notre propos.

On attache à la surface moyenne  un repère orthonormé local Oxyz associé au plan tangent de la
structure différent du repère global OXYZ . La position des points de la plaque est donnée par les
coordonnées cartésiennes  x , y  de la surface moyenne et l’élévation z par rapport à cette surface.

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 4/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

2.2 Cinématique
La cinématique de base de l'élément T3G est celle de Hencky-Mindlin (voir [7] pour les détails), où le
tenseur de déformations 3D,  , est défini comme suit :

 xx =e xx z  xx  yy=e yy z  yy  xy =e xyz  xy

 xz = x  yz = y (1)
où
∂u ∂v 1 ∂u ∂v
e xx= e yy= e xy=   
∂x ∂y 2 ∂y ∂x

∂ x ∂ y 1 ∂ ∂
 xx=  yy =  xy=  x  y 
∂x ∂y 2 ∂ y ∂x

 x = y  y =− x

∂w ∂w
 x = x   y = y  (2)
∂x ∂y

e   étant le tenseur d’extensions membranaires,    le tenseur de variation de courbure,  x ,  y

les variables auxiliaires définies ci-dessus et  x ,  y les distorsions transverses. Les variables
principales sont les déplacements u , v , w et les rotations  x ,  y et  z . Dans les formulations
considérées ici la rotation normale à la surface moyenne,  z , est introduite exclusivement pour des
raisons de compatibilité entre les degrés de liberté nodaux pour des géométries non-planes.

2.3 Support élément fini

L’élément fini a comme support une maille T3, avec 6 degrés de liberté par nœud, ce qui fait 18
degrés de liberté au total. Les caractéristiques principales sont données dans le tableau ci-dessous.
On définit les mêmes fonctions de forme pour les déplacements, u , v , w  , que pour les rotations,
 x ,  x  .

TRIA3 N i i=1, n

i=1 à 3 • Trois points de collocation

(4,5,6) pour déterminer  x et
N 1  ,  ==1−− y

N 2  ,  = • Un point d'intégration au

1 1
3
N  ,  = = ; = 
3 3

Tableau 2.1 : Élément fini de référence T3G utilisé dans la formulation Q4GG à base triangulaire

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 5/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

2.4 Membrane-Flexion
La présente formulation est basée sur 5 degrés de liberté en négligeant tout effet de la rotation
perpendiculaire au feuillet moyen. Cette troisième rotation, représentée par  z dans le repère local,
n'est pas présentée dans cette formulation, voir [7] pour plus de détails.

L’extension membranaire et la variation de courbure sont discrétisées d’une manière équivalente. On

utilise les mêmes fonctions de forme pour les déplacements et les rotations, ce qui conduit à :

U i= N r u ri ,  i =N r ir r =1 à 3
r
u i et  ri
étant les valeurs nodales des déplacements et des rotations (transformées),
respectivement. En suivant les équations (1 ) et (2) on a :
e ij =B rijk u rk ,  ij =B rijk  rk ,

r 1 ∂ N r −1 −1
Bijk =  J   J li  jk 
2 ∂ l lj ik

Où on introduit le jacobien de la transformation sur l'élément de référence :

∂ xi ∂ N r r
J ij = = xi
∂ j ∂  j

Ainsi, pour les fonctions définies dans le Tableau 2.1 et en choisissant le repère local d’un élément tel
que x 1= y 1= y 2 =0 , on obtient :

J=
 
x2
0
x3
y3
; J−1=

1 y3 −x3
x2 y3 0 x2 
∂ N r −1 −1
∂ k
= 1
0
0
1 
et au final le tenseur de déformation membranaire uniforme sur l'élément :


u1
v1

   
e 11 − y3 0 y3 0 0 0
1 u2
e 22 = 0 x 3−x 2 0 −x 3 0 x2 (3)
x 2 y3 v2
2 e 12 x 3−x 2 − y 3 −x3 y 3 x 2 0
u3
3
v

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 6/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2


 1x
1y

  
 11 0 − y3 0 y3 0 0
1  2x
 22 = − x 3− x 2  0 −x 3 0 −x 2 0 (4)
x2 y3 2
2  12 y3 x 3− x 2 − y 3 −x 3 0 x 2 3y
x
 3y

2.5 Distorsion transverse

2.5.1 Mise en évidence d'un verrouillage en cisaillement transverse
Pour éviter le verrouillage en cisaillement, on doit judicieusement choisir la discrétisation de la
distorsion transverse. Par exemple, en imposant,

∂w ∂w
 x = x  et  y = y  (5)
∂x ∂y
Soit en fonction des degrés de liberté :

  

∂w   ∂ N −1 
 y N  y J k1 w


x
y
=
− x 
∂ x
∂w
∂y
=
−N   x 
∂ k
∂ N  −1 
J k2 w
∂ k

La limite des plaques minces, où h  0 et par conséquent,  x  0 ,  y  0 mène aux équations

suivantes en tout point  ,  :

 
1
1y   w2 −w 1 2y − 1y  3y − 1y 
x
y
=
− 
1
1
x
x2
 x −x  w − x w  x w − 2x − 1x − 3x −1x 
x2 y3 3 2 1 3 2 2 3
0
= 0 , ∀, (6)

Les conditions dans l'équation 6 conduisent directement aux relations  3x = 2x = 1x et

 3y = 2y = 1y . Les rotations de flexion sont donc toutes égales et forcément nulles si les mouvements
de corps rigide sont correctement supprimés du système. On observe ainsi le verrouillage en
cisaillement : une plaque mince ( h≈0 ) avec la cinématique de la distorsion transverse décrite dans
l'équation 5 ne peut pas se déformer en flexion.

2.5.2 Remède aux verrouillages par l'approche T3G

Pour remédier aux verrouillages et de façon similaire à la formulation Q4G [7], la formulation T3G
utilise une interpolation de la distorsion transverse indépendante des déplacements et des rotations.
Les relations (5) sont relâchées, grâce à trois points de collocation. Primo, on suppose que la variation
spatiale de la distorsion transverse est linéaire avec une composante tangentielle constante sur
chaque arrête. Secundo, comme dans Q4G, la compatibilité est imposée aux milieux des trois arrêtes
et seulement dans la direction tangentielle.

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 7/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

Les valeurs de la distorsion transverse le long des arrêtes peuvent être calculées directement comme
suit :

w j −w i 1 ki kj
 ks =   s  s  (7)
Lk 2

où  ks est la composante tangentielle de la distorsion transverse sur l’arrête k , qui se trouve entre
ki kj
les nœuds i et j . L k est la longueur de l’arrête k ;  s et  s sont les valeurs projetées des
rotations (transformées) sur la tangente de l’arrête k , correspondant aux nœuds i et j ,
respectivement :

 kis = ix cos  k iy sin  k

où  k représente l’angle de l’arrête k par rapport à l’abscisse du repère local. L’équation (7) peut
être réécrite comme suit :


1
x
 1y

 
1 1 1 1 1 1 w1
C S1 − C S1 0 0 0
2 1 2 L1 2 1 2 L1


1
 s  2x
2 1 1 1 1 1 1
 = 0
s 0 0 C S2 − C S2  2y (8)
3
2 2 2 L2 2 2 2 L2
 w2
1
s 1 1 1 1 1
C S 0 0 0 C S3 −  3x
2 3 2 3 L3 2 3 2 L3
 3y
w3

ou C k =cos  k et S k =sin  k

D’autre part, on impose la forme suivante de la variation spatiale dans le repère de référence de la
distorsion transverse :

 x = a x b x c x 
(9)
 y = a y b y c y 

où a x , a y , b x , b y , c x , c y sont les paramètres à déterminer à partir des hypothèses selon

lesquelles les composantes tangentielles de  sont constantes sur chaque arrête.

On obtient donc :

 1s= x | =0=a x b x   a x = 1s , b x =0

 2s = x C 2 y S 2 | =1=[a x C 2a y S 2c x C 2b y S 2−c x C 2−c y S 2 ]

 a x C 2a y S 2c x C 2= 2s ; b y S 2−c x C 2−c y S 2=0

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 8/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

 3s=  x C 3 y S 3 =0=a x C 3a y S 3 c x C 3c y S 3 

 a x C 3a y S 3= 3s ; c x C 3c y S 3=0

Ce qui mène au résultat final

C3 1 1 3
a x = 1s , a y =−  
S 3 s S3 s

b x =0 , 1 1 1 2 1 3
b y =−  − 
QS 2 s S 2 s S 3 s
S C3 1 C3 2 C S
c x =− 1sQ  2s −Q 2  3s , c y =  s −Q  sQ 3 2 2  3s (10)
S3 S3 S3 S3

C3
où Q=1/C 2−S 2 
S3

A partir des équations 9 et 10 on peut écrire :

  
S2
1− Q − Q  1s

x
y
=
C
− 3−
1
S 3 QS 2
C
 3 
S3
1
S2
C
− 3 Q 
S3
1 1
S3
C S
s
2

−  3 2 2 Q   3s
S3 S3
(11)

S3

En combinant les équations 8 et 11 on obtient :

  
S2
1− Q −Q
x S3
= ×
y C3 1 C 1 C 1 1 C3 S 2
− −  3  − 3 Q  −  2 Q 
S 3 QS 2 S3 S2 S3 S3 S3 S3


 1x
 1y

 
1 1 1 1 w1
0 − 0 0 0 0
2 L1 2 L1
 2x
1 1 1 1 1 1
0 0 0 C S − C S2  2y (12)
2 2 2 2 L2 2 2 2 L2
w2
1 1 1 1 1 1
C S 0 0 0 C S3 −  3x
2 3 2 3 L3 2 3 2 L3
 3y
w3

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 9/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

Les expressions ci-dessus ont été légèrement simplifiées en utilisant le même repère local que dans
[§2.4], c’est-à-dire celui où  1=0 et donc C 1=1 et S 1=0 . La relation entre les degrés de liberté
de l’élément, U i et la distorsion transverse en tout point  ,  , de l'équation 12 s’écrit sous la
forme :

 i =Bikct   , U k
Où

 
S2
1− Q −Q
S3
B ct = ×
C3 1 C 1 C 1 1 C 3 S2
− −  3  − 3 Q  −  2 Q 
S 3 QS 2 S3 S2 S3 S3 S3 S3

 
1 1 1 1
0 − 0 0 0 0
2 L1 2 L1
1 1 1 1 1 1
0 0 0 C S − C S2
2 2 2 2 L2 2 2 2 L2
1 1 1 1 1 1
C S 0 0 0 C S3 −
2 3 2 3 L3 2 3 2 L3

et
T
U =  1x  1y w 1  2x  2y w2  3x  3y w3 

Contrairement à la discrétisation exposée dans [§2.5.1], et dans l'équation 5 , la limite h  0

appliquée à la cinématique indépendante conduit à trois équations qui couplent d’une manière plus
équilibrée les déplacements w , et les rotations  x ,  y , suivant l'équation 8 .


 1x
 1y

 
1 1 1 1 1 1 w1
C S1 − C C2 0 0 0
2 1 2 L1 2 1 2 L1
 2x


0
1 1 1 1 1 1
0= 0 0 0 C S2 − C S2  2y (13)
2 2 2 L2 2 2 2 L2
0 w2
1 1 1 1 1 1
C 0 0 0 C S −  3x
2 3 S3 L3 2 3 2 3 L3
 3y
w3

Les conditions dans l'équation 13 ne provoquent aucun blocage sur les degrés de liberté de
rotation. Par ailleurs, dans l'équation 13 on a utilisé l'équation 11 pour déduire que

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 10/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

  
x
y
=
0
0

en tout point, seulement si :

  
 1s 0
 2s = 0
 3s 0

2.6 Calcul des contraintes et des forces internes

On présente ici l’application de la cinématique définie au [§2.4] au calcul des forces internes. Deux
cas différents seront considérés :

• Élasticité linéaire,
• Modèle global de type GLRC.

Pour l'élasticité linéaire, un seul point d’intégration dans l’épaisseur est nécessaire, car les
déformations et les contraintes varient linéairement dans l’épaisseur.

L’objectif principal de la modélisation Q4GG est l’utilisation du modèle GLRC_DAMAGE [R7.01.31] dans
le cadre de l’élément T3G et Q4G. Un seul point d’intégration dans l'épaisseur est nécessaire, car la
non-linéarité est prise en compte directement en termes d'efforts et de déformations généralisés.

2.6.1 Élasticité linéaire

Pour un comportement élastique linéaire on définit l’énergie potentielle de la manière suivante :

h
1
 = ∫S ∫−h2 ij C elijkl  kl dz dS ext
L
2 2

1 M M CT ext
2 ∫S ij ijkl kl ij ijkl kl i ij j
= [e H e  H   H  ]dS  (14)

el
où  et C sont, respectivement, le tenseur 3D de la déformation et le tenseur élastique
correspondant. Pour un matériau linéaire, l’intégrale selon z s’effectue analytiquement. ext
représente la contribution due aux conditions limites. Les tenseurs e ,  , représentent l’extension
membranaire et la variation de courbure et le vecteur  la distorsion transverse, comme introduit
dans [§2.4]. Dans l'équation 14 , on fait aussi l’hypothèse que la coque est symétrique par rapport
au feuillet moyen, sinon un terme de couplage membrane-flexion supplémentaire est obtenu. Les
détails du calcul dans l'équation 14 sont fournis dans [7], où on trouve les expressions suivantes
pour les matrices H :

   
1  0 1  0

 
3
Eh  1 0 Eh  1 0 kEh 1 0
HM= H F= H CT =
1− 2 0 1− 12 1− 2 0 1− 21  0 1
0 0
2 2

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 11/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

où h est l’épaisseur, E le module de Young,  le coefficient de Poisson et k un paramètre qui

peut varier selon la théorie appliquée. Le plus souvent on utilise k =
5 (modèle Reissner).
6
Ensuite, on applique la cinématique introduite aux [§2.3] et [§2.4] :


1x

 
 1y
u1  1x
w1
v1  1y

 
e xx  xx  2x


2 2
M u F x x CT
e yy =B  yy =B =B  2y
v
2
 2y y
e xy  xy w2
u3  3x
v3 3x
 3y
 3y
3
w

pour obtenir à partir de l'équation 14 :

L 1 T M M M 1 T F F F
 = U  [∫S B ij  H ijkl B kl  dS ] U  U  [∫S B ij  H ijkl Bkl  dS ] U  
2  2 
K M K F 
1 T CT CT CT ext
U  [∫S Bij  H ijkl B kl  dS ] U  
2 
K CT


En ayant recours au principe d’énergie libre minimale, ∂ L=0 , ∀  U  , on obtient :

K M K F K CT
 
U  =F ext

int
F

où Fint et Fext sont, respectivement, les forces internes et externes.

L’intégration de toutes les contributions à la matrice de rigidité K M , K F et K CT s’effectue avec un

seul point d’intégration dans le plan du triangle (centre de gravité). Ceci donne une intégration exacte
pour les termes K M et K F car H M , H F , B M , B F sont constantes par élément (voir [§4.5]). En
revanche, B CT est une fonction linéaire (voir [§2.4]) et l’intégration exacte de K CT exigerait trois
points de Gauss.

Comme l’élément T3G, tel que proposé dans [3] et discuté dans [12], est sous-intégré sur la partie du
cisaillement transverse, il a un mode parasite, c’est-à-dire un mode à énergie zéro, qui n’est pas un
mode de corps rigide. À partir de deux éléments, le modèle T3G retrouve cependant un « rang »
correct.

2.6.2 Comportement global non-linéaire de type GLRC

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 12/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

Lors de l’utilisation d’un modèle global, la relation entre les contraintes et les déformations est définie
en terme des variables généralisées,  G ,  G :

 
e N
G G
 =   = M
 T
où les déformations généralisées sont l’extension membranaire e , la variation de courbure  , et la
distorsion transverse  (voir [§2.4] pour des définitions plus précises). Les contraintes généralisées
sont l’effort membranaire N , le moment fléchissant M , et le cisaillement transverse T :

    
N xx h  xx M xx h  xx
1 1
N= N yy = ∫ 2h  yy dz M= M yy = ∫ 2h z  yy dz
2 −2 2 −2
N xy  xy M xy  xy

   
h
Tx 
T= =∫ 2h xz dz
Ty −
2  yz

Le modèle GLRC_DAMAGE prend en compte l’élasto-plasticité en membrane-flexion, principalement

due aux aciers, et l’endommagement en flexion, principalement dû au béton. L’énergie libre
correspondant au modèle est définie comme :

G 1 P M p 1 P Fd p
 =
2 ∫S
e−e ij H ijkl e−e kl dS ∫S − ij H ijkl d 1, d 2− kl dS
2
(15)

1 CT

2 ∫S
 ij H ijkl  kl dS

où e p et  p sont, respectivement, l’extension membranaire plastique et la variation de courbure

plastique ; H M et H CT sont les matrices élastiques pour la membrane et le cisaillement transverse,
les mêmes que celles définies au [§2.6.1]. Pour la partie de flexion on définit la matrice H Fd , qui est
une fonction des variables d’endommagement, d 1 et d 2 [8]. On fait remarquer qu’aucune non-
linéarité n’est considérée pour la partie du cisaillement transverse. Cette hypothèse n’est pas vraiment
justifiée, car l’endommagement du béton induit forcément une diminution de la raideur en cisaillement
transverse. Puisque le modèle a été conçu dans le cadre de l’hypothèse des plaques minces, il est
déconseillé d’utiliser le modèle GLRC_DAMAGE dans des situations où l’énergie de cisaillement
transverse devient comparable aux contributions de la membrane ou de la flexion.

En imposant ∂ G=0 , on obtient:

  
N xx M xx
G  M T


N xy

  F T
  = U ∫S  B   N yy dS   ∫S  B   M yy  B
M xy





 CT T H CT  x dS
y (16)

M F
F int Fint

 w ∫s  B CT

T
w  H
CT  x

y
dS

F CT
int

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 13/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

où on utilise les relations suivantes :

h h h
B M =∫ 2
h
M
B dz B F =∫ 2
h
F
z B dz B CT =∫ 2
h BCT dz
− − −
2 2 2

La loi de comportement, dont pratiquement toutes les caractéristiques découlent de la définition de

l’énergie libre [éq .4.7.3-1], peut être résumée, suivant l'équation 16 , en terme de relations :

N=Ne ,   M=M e ,   T=H CT 

Les détails de l’application du modèle GLRC_DAMAGE dans un cadre « éléments finis » sont donnés
dans [4] et [8].

2.7 Matrice de masse « lumpée »

Dans ce paragraphe on présente l’approche utilisée pour construire la matrice de masse « lumpée ».
On se focalise sur les termes dus à la flexion, les termes correspondant à la membrane étant obtenus
d’une manière classique appliquée aussi aux éléments 2D . La matrice de masse M est définie à
partir de l’énergie cinétique E cin et de la cinématique proposée au [§2.2]:


h
̇ x
1
E cin= ∫ ∫  u̇2 dz dS = ̇x ̇y ẇ  M   ̇y  E cin
2 S −
2
h memb
2 
ẇ

où E cin
memb
représente la contribution de la membrane,  étant la masse volumique et la matrice de
masse M étant égale à :

 
h3
0 0
12
M   =m  h3
0 0 (17)
12
0 0 h

m  =∫s  N  N  dS (18)

Une manière de calculer l’intégrale 18 passe par la quadrature de Gauss, qui conduit à la matrice
cohérente couplant les contributions des différents degrés de liberté. Il est généralement reconnu que
le calcul d’une matrice de masse cohérente ne se justifie pas pour les calculs en dynamique rapide,
pour lesquels la matrice « lumpée » offre toujours un meilleur rapport précision-coût.

On utilise la même matrice de masse pour le T3G que pour le DKT. En plus de l’approximation du
découplage des différents degrés de liberté, conduisant à l’utilisation de la matrice lumpée, on fait des
approximations sur les termes d’inertie afin d’augmenter le pas de temps de stabilité.

Dans [5] on propose de construire la matrice lumpée à partir de l'équation 18 en utilisant la
quadrature de Lobatto, dont les variantes sont le schéma trapézoïdal et le schéma de Simpson, où les
points d’intégration coïncident avec les nœuds. La construction de la matrice de masse se fait à

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 14/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

travers un élément fini de poutre, linéaire et à deux nœuds, en utilisant le schéma trapézoïdal,
conduisant à :

 
I
0 0 0
A
1
M 0pout =  LA 0 1 0 0
2 I (19)
0 0 0
A
0 0 0 1

T
pour laquelle le vecteur de degrés de liberté s’écrit comme 1 w 1  2 w 2 . A , L et I sont l’aire de
la section, la longueur et le moment d’inertie de l’élément poutre, respectivement. L’utilisation de la
matrice eq. 19 semblant trop restrictive par rapport à la condition de stabilité du schéma
d'intégration explicite en temps, on propose dans [5] plutôt :

 
 0 0 0
pout 1 0 1 0 0
M 0 =  LA
2 0 0  0
0 0 0 1

où le paramètre  est introduit afin que son ajustement puisse maximiser le pas de temps de
1 2
stabilité. Selon [5] sa valeur optimale serait = L . En appliquant directement ces résultats par
8
analogie aux plaques on remplace la matrice de l'équation 19 par :

 
hAe
0 0
8
M   =m  hAe (20)
0 0
8
0 0 h

où Ae est l’aire de l’élément considéré. Dans le passage de la poutre à la plaque on a supposé une
certaine équivalence entre la longueur de l’élément poutre et l’aire de l’élément plaque, de sorte que
Ae ≈ L 2 . On rappelle que la démarche proposée dans [5] n’est pas rigoureuse du point de vue
géométrique et qu’elle se focalise sur la maximisation du pas de stabilité. Dans la version implantée
dans Code_Aster on s'assure de l’effet souhaité de la modification de 18 en 19 en utilisant
(comme cela a été fait pour les éléments de la famille DKT [7]) :

 
h3 hAe
max  ,  0 0
12 8
M   =m  h
3
hA
e
(21)
0 max  ,  0
12 8
0 0 h

car l'équation 20 n’est intéressante que pour les maillages grossiers, a priori plus courants,
tandis que l'équation 17 devient favorable pour des maillages très fins. D’ailleurs, dans
Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 15/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

Code_Aster on décide d'intégrer (18) par la quadrature de Gauss à trois points en ne gardant que
les termes diagonaux :

m  =  ∫s  N  N  dS (22)

Remarque : la matrice de masse cohérente classique n'a pas été développée pour cet élément fini
dans Code_Aster.

2.8 Second membre correspondant au travail extérieur

La formulation variationnelle du travail extérieur et la discrétisation du potentiel des efforts extérieurs
sont présentés dans §4.8 de [R3.07.03].
Remarque : Le chargement thermique n'est pas traité.

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 16/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

3 Implantation des éléments de plaque dans Code_Aster

3.1 Description
Ces éléments (de noms MEQ4GG4, MET3GG3) s'appuient sur des mailles QUAD4 et TRIA3 planes. Ces
éléments ne sont pas exacts aux nœuds et il faut mailler avec plusieurs éléments pour obtenir des
résultats corrects.

3.2 Utilisation et développements introduits

Ces éléments s'utilisent de la façon suivante :

• AFFE_MODELE (MODELISATION='Q4GG',...) pour le triangle et le quadrangle.

• AFFE_CARA_ELEM ( COQUE=_F(EPAISSEUR='EP'
ANGL_REP = ( '  ' '  ' )
COEF_RIGI_DRZ = 'CTOR')

Pour faire des post-traitements (efforts généralisés,...) dans un repère choisi par l’utilisateur qui n’est
pas le repère local de l’élément, on donne une direction de référence d définie par deux angles
nautiques dans le repère global. La projection de cette direction de référence sur le plan de la plaque
fixe une direction X1 de référence. La normale au plan en fixe une seconde et le produit vectoriel
des deux vecteurs précédemment définis permet de définir le trièdre local dans lequel seront
exprimés les efforts généralisés et les contraintes. L’utilisateur devra veiller à ce que l’axe de
référence choisi ne se retrouve pas parallèle à la normale de certains éléments de plaque du modèle.
Par défaut cette direction de référence est l’axe X du repère global de définition du maillage.

La valeur CTOR correspond au coefficient que l’utilisateur peut introduire pour le traitement des termes
de rigidité et de masse suivant la rotation normale au plan de la plaque. Ce coefficient doit être
suffisamment petit pour ne pas perturber le bilan énergétique de l’élément et pas trop petit pour que
les matrices de rigidité et de masse soient inversibles. Une valeur de 10−5 est mise par défaut.
Pour un comportement :

• Élastique isotrope homogène dans l’épaisseur on utilise le mot-clé ELAS dans

DEFI_MATERIAU où l’on définit les coefficients E module d’Young,  coefficient de
Poisson,  la masse volumique :
ELAS =_F (E = young NU = nu , RHO = rho ... )

• Élastoplastique endommageable de type GLRC_DAMAGE, ce modèle global de plaque en

béton armé est capable de représenter son comportement jusqu'à la ruine. Les
caractéristiques du béton et des armatures sont données dans DEFI_GLRC.

DEFI_GLRC ( RELATION = “GLRC_DAMAGE“,

BETON = _F(MATER = ... , EPAIS= ...) ,
NAPPE = _F(MATER = ... , OMX = ...))

AFFE_CHAR_MECA (

DDL_IMPO =_F ( DX =.. DY =.. DZ =.. DRX =.. DRY =.. DRZ =.. degré de
liberté de plaque dans le repère global.

FORCE_COQUE =_F (FX =.. FY =.. FZ =.. MX =.. MY =.. MZ =.. ) Il s’agit des
efforts surfaciques (membrane et flexion) sur des éléments de plaque. Ces efforts peuvent
être donnés dans le repère global ou dans le repère utilisateur défini par ANGL_REP.
Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 17/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

FORCE_NODALE =_F (FX =.. FY =.. FZ =.. MX =.. MY =.. MZ =.. ) Il s’agit
des efforts de coque dans le repère global.

3.3 Calcul en élasticité linéaire

La matrice de rigidité et la matrice de masse (respectivement les options RIGI_MECA et MASS_MECA)
sont intégrées numériquement. On ne vérifie pas si la maille est plane ou non. Le calcul tient compte
du fait que les termes correspondant aux degrés de liberté de plaque sont exprimés dans le repère
local de l’élément. Une matrice de passage permet de passer des degrés de liberté locaux aux degrés
de liberté globaux.

Les calculs élémentaires disponibles actuellement correspondent aux options :

DEGE_ELNO : qui donne les déformations généralisées par élément aux nœuds à partir des
déplacements dans le repère utilisateur : EXX, EYY, EXY, KXX, KYY, KXY, GAX,
GAY.
EFGE_ELNO : qui donne les efforts (contraintes généralisées) par élément aux nœuds à partir des
déplacements : NXX, NYY, NXY, MXX, MYY, MXY, QX, QY.
SIEF_ELGA : qui donne les efforts (contraintes généralisées) par élément aux points de Gauss à
partir des déplacements : NXX, NYY, NXY, MXX, MYY, MXY, QX, QY.
EPOT_ELEM : qui donne l’énergie élastique de déformation par élément à partir des
déplacements.
ECIN_ELEM : qui donne l’énergie cinétique par élément.

Enfin on calcule aussi l’option FORC_NODA de calcul des forces nodales pour l’opérateur
CALC_CHAMP.

3.4 Calcul non-linéaire

La matrice de rigidité est là aussi intégrée numériquement. On fait appel à l’option de calcul
STAT_NON_LINE dans laquelle on définit au niveau du comportement non linéaire le nombre de
couches à utiliser pour l’intégration numérique.

STAT_NON_LINE (....
COMPORTEMENT =_F (RELATION ='GLRC_DAMAGE'
....)

Pour la modélisation Q4GG, les seules lois de comportement utilisées sont des lois globales (puisqu’il
n’y a qu’un point d’intégration dans l’épaisseur), reliant les déformations généralisées aux contraintes
généralisées.

Les calculs élémentaires disponibles actuellement correspondent aux options:

• DEGE_ELNO qui fournit les déformations généralisées par élément aux nœuds dans le repère
utilisateur à partir des déplacements.
• SIGM_ELNO qui permet d’obtenir le champ de contraintes dans l’épaisseur par élément aux
nœuds pour tous les sous-points (toutes les couches et toutes les positions : en peau
inférieure, au milieu et en peau supérieure de couche).
• EFGE_ELNO qui permet d’obtenir les efforts (contraintes généralisées) par élément aux
nœuds dans le repère utilisateur.
• VARI_ELNO qui calcule le champ de variables internes et les contraintes par élément aux
nœuds pour toutes les couches, dans le repère local de l’élément.

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 18/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

4 Validation de la modélisation Q4GG

4.1 SDLS123 – Poutre droite avec amortissement de Rayleigh (V2.03.123)
Ce test représente le calcul dynamique élastique avec amortissement de Rayleigh d'une poutre droite,
encastrée aux deux extrémités et soumise à une pression constante. Les résultats obtenus sont
comparés aux résultats obtenus avec Europlexus.
Deux modélisations sont effectuées :
• Modélisation A permet de tester le modèle Q4GG avec des QUAD4,
• Modélisation B permet de tester le modèle Q4GG avec des TRIA3.

4.2 SDNS108 – Réponse dynamique d'une dalle en béton armé (modèle

GLRC_DAMAGE) appuyée sur 4 côtés soumise à une charge
concentrée : régime de 4.2 plaque élasto-plastique (V5.06.108)
Ce test représente le calcul d’une dalle en béton armé, en flexion, soumise à une charge concentrée.
Il permet de valider la modélisation Q4GG utilisant le modèle matériau global GLRC_DAMAGE en la
comparant à un calcul analogue avec Europlexus. La dalle est simplement appuyée sur les quatre
côtés.
Deux modélisations sont effectuées :
• Modélisation A permet de tester le modèle Q4GG avec des QUAD4,
• Modélisation B permet de tester le modèle Q4GG avec des TRIA3.
P la q u e
Z Y
B é to n a rm é
dy

dx
D C

ez
h
A B
a r m a t u r e d i a m è t r e  2
f e u ille t m o y e n
a r m a t u r e d i a m è t r e  1 N a p p e d ’a r m a tu r e s é q u iv a le n t e f e u ille t m o y e n

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)

 Version
Code_Aster default
Titre : Éléments de plaque : Modélisation Q4GG Date : 25/09/2013 Page : 19/19
Responsable : POTAPOV Serguei Clé : R3.07.09 Révision :
974be9fbc1a2

5 Bibliographie
1. J.L. Batoz, G.Dhatt, « Modélisation des structures par éléments finis : poutres et plaques »,
Hermès, Paris, 1992.

2. T.A. Rock, E. Hinton, « A finite element method for the free vibration of plates allowing for
transverse shear deformation », Computers and Structures, Vol.6, p.37-44,1976.

3. T.J.R. HUGHES, R.L. TAYLOR, « The linear triangle bending elements», dans The
Mathematics of Finite Element and Application IV , MAFELAP 1981 , Academic Press,
London, 1982, p. 127-142.

4. P. KOECHLIN, « Modèle de comportement membrane-flexion et critère de perforation pour

l’analyse de structures minces en béton armé sous choc mou », Thèse de doctorat Paris VI,
2007.

5. T.J.R. HUGHES, M. COHEN, M. HAROUN, « Reduced and selective integration techniques in

the finite element analysis of plates», Nuclear Engineering and Design, vol. 46 (1978), p. 203-
222.

6. D. MARKOVIC, «Implantation d'un nouvel élément fini de coque épaisse (T3GS) dans
Europlexus », H-T62-2008-00080-FR, février 2008.

7. R3.07.03 « Élément de plaque: modélisations DKT, DST, DKTG et Q4G ».

8. R7.01.31 « Loi de comportement de plaques en béton armé GLRC_DAMAGE ».

9. V2.03.123 « Poutre droite avec amortissement de Rayleigh (Comportement élastique) ».

10. V5.06.108 « Réponse dynamique d'une dalle en béton armé (Modèle GLRC_DAMAGE)
appuyée sur 4 cotés soumise à une charge concentrée: régime élastoplastique ».

11. T. BELYTSCHKO, I. LEVIATHAN « Physical stabilization of the 4-node shell element with one
point quadrature », Comp. Methods Appl. Engrg., vol. 113 (1994), p. 321-350.

12. R. AYAD, G. DHATT, J. L. BATOZ, « A new hybrid-mixed variational approach for Reissner-
Mindlin plates. The MiSP model », Int. J. Numer. Meth. Engng., vol. 42 (1998), p. 1149-1179.

Manuel de référence Fascicule r3.07: Eléments mécaniques à surface moyenne

Document diffusé sous licence GNU FDL (http://www.gnu.org/copyleft/fdl.html)