Methode Des Elements Finis Approche Pratique
Methode Des Elements Finis Approche Pratique
Methode Des Elements Finis Approche Pratique
Mthode
deslments finis
Approche pratique
en mcanique des structures
9782100544639.indb 3
09/02/10 14:53
9782100544639.indb 4
09/02/10 14:53
1 Introduction
15
19
22
27
31
31
4.2 Rsolution
33
43
15
45
45
45
46
47
54
61
61
63
72
78
89
9782100544639.indb 5
09/02/10 14:53
7 lments isoparamtriques
7.1 Problmatique du maillage
101
101
104
120
123
8 lments de membrane
127
128
140
155
176
9 lments de plaque
185
185
192
208
10 lments de coque
215
215
218
221
227
227
240
254
265
275
Bibliographie
281
Complments en ligne
283
1. Notice dutilisation
283
285
Index alphabtique
287
VI
9782100544639.indb 6
09/02/10 14:53
1 Introduction
9782100544639.indb 1
09/02/10 14:53
1 Introduction
ainsi proposs. Quand cela savre ncessaire, ceux-ci sont accompagns de rappels sur les thories des poutres, des plaques ou des coques. Leurs rsultats seront
dailleurs, comme indiqu prcdemment, trs souvent utiliss comme rfrentiels.
Ceci tant, cet ouvrage ne prtend pas couvrir la mthode de manire exhaustive
car les techniques numriques abordes font partie des plus courantes dans les
codes de calcul commerciaux.
Partant de pr requis en mathmatiques et mcanique du solide, le principe dapproximation lmentaire utilis par cette mthode est tout dabord appliqu en
statique aux structures filaires en barres et poutres. La problmatique du maillage
et de la validation des modles de calcul est ensuite aborde lors de ltude des
modlisations surfaciques avec des lments membranes, plaques ou coques. Enfin,
ces lments sont ensuite utiliss dans le cadre de ltude de la non linarit gomtrique et des mthodes de rsolution associes telles que celles de Newton-Raphson
ou de longueur darc. Ces diffrents aspects sont ensuite appliqus aux instabilits
des structures comme le flambement, le dversement ou le voilement des mes.
Afin de bien dcrire la mthodologie utilise, la grande majorit de ces exemples
est traite pas pas par des calculs manuels ou semi automatiques avec le logiciel
Mathcad dvelopp par la socit Parametric Technology Corporation et dont
les rsultats sont recoups avec les codes de calcul Advance Structure/Effel ou
A
baqus dits respectivement par les socits Graitec SA et Dassault Systmes
Simulia Corporation. Nous remercions dailleurs trs sincrement ces trois diteurs
pour nous avoir permis dutiliser leurs logiciels pour illustrer nos exemples.
9782100544639.indb 2
09/02/10 14:53
A
Le calcul par lments finis ncessitant le maniement de nombreuses valeurs numriques, il est plus ais dexprimer celles-ci sous forme matricielle.
En regroupant des termes de mme nature au sein dune seule et mme variable,
cette criture plus synthtique permet en effet une meilleure comprhension des
diffrentes phases de construction de la mthode.
Ceci ncessite nanmoins la matrise des oprations de base associes ce type de
calcul : laddition ou le produit de plusieurs matrices, la rsolution de systmes
linaires, etc.
v1 = v ( 0 ) = b0
1 = v ( 0 ) = b1
2 + b L3
v
v
L
b
b
L
b
L
=
=
+
+
(
)
2
0
1
2
3
2 = v ( L ) = b1 + 2b2 L + 3b3 L2
(2.3)
9782100544639.indb 3
09/02/10 14:53
Anticipant sur les rgles relatives au produit des matrices (cf. 2.2.2), On a donc:
v1 1
0
1
v
=
{ } =
v 2 1
2 0
0 0
0 b0
1 0
0 b1
= [ R ] {b}
L L2 L3 b2
1 2 L 3L2 b3
a11
a
21
a31
[ A] = .
ai1
.
a
n1
a12
a22
a32
.
ai 2
.
an 2
a13
a23
a33
.
ai 3
.
an 3
.
.
.
.
.
.
.
a1 j
a2 j
a3 j
.
aij
.
anj
.
.
.
.
.
.
.
a1m
a2m
a3m
.
aim
.
anm
(2.4)
9782100544639.indb 4
09/02/10 14:53
c13
c 23
c33
.
ci 3
.
cn 3
.
.
.
.
.
.
.
c1 j
c2 j
c3 j
.
cij
.
cnj
.
.
.
.
.
.
.
c1m
c 2m
c 3m
. =
cim
.
cnm
c11 c12
c
c
21 22
c31 c32
[C ] = . .
c i1 c i 2
.
.
c
n1 cn 2
a11 + b11
a + b
21 21
a31 + b31
ai1 + bi1
a + b
n1 n1
a12 + b12
a22 + b22
a32 + b32
.
ai 2 + bi 2
.
an 2 + bn 2
a13 + b13
a23 + b23
a33 + b33
.
ai 3 + bi 3
.
an 3 + bn 3
. a1 j
. a2 j
. a3 j
.
. aij
.
. anj
+ b1 j
+ b2 j
+ b3 j
.
+ bij
.
+ bnj
. a1m
. a2m
. a3m
.
. aim
.
. anm
+ b1m
(2.5)
+ b2m
+ b3m
.
= [ A] + [B ]
+ bim
+ bnm
(2.6)
1 0 2 3
3 4 7 2
Exemples: Soit les matrices [ A ] =
et [ B ] =
3 1 0 5
1 6 2 3
4 4 9 5
2 4 5 1
[C ] = [ A ] [ B ] =
.
7 2 8
2 5 2 2
[C ] = [ A ] + [ B ] = 4
2.2.2 Produit
Produit dune matrice par un scalaire
9782100544639.indb 5
1 0 2 3 3 0 6 9
=
1 0 5 9 3 0 15
[C ] = 3 [ A ] = 3 3
09/02/10 14:53
Produit de 2 matrices
cij =
k =m
aik bkj
(2.8)
k =1
7 5
1 4 + 0 3 + 2 1 + 3 7 1 9 + 0 1 + 2 6 + 3 5 27 36
[C ] = [ A ] [ B ] = 3 4 + 1 3 + 0 1 + 5 7 3 9 + 1 1 + 0 6 + 5 5 = 50 53
Produit de 3 matrices
4
3
1 0 2 3
B
=
,
Exemple: Soit les matrices [ A ] =
[
]
1
3
1
0
5
9
27 36
1
et [C ] =
6
50 53
4 9
1 0 2 3 3 1 27 36 27 36 27
([ A ] [ B ]) [C ] = 3 1 0 5 1 6 50 53 = 50 53 50
7
5
4 9
1 0 2 3 3 1 27 36 27 36 27 36 2529
([ A ] [ B ]) [C ] = 3 1 0 5 1 6 50 53 = 50 53 50 53 = 4000
7
5
1. Le produit
[A][B] [B][A].
Gnralement,
36 2529 28880
=
53 4000 4609
28880
4609
9782100544639.indb 6
09/02/10 14:53
558 621
=
=
6 50 53 3 1 0 5 327 354 4000 4609
5
439 3117
4 9
558 621
1 0 2 3 3 1 27 36
1 0 2 3 131 161
2529 2880
[ A ] ([ B ] [C ]) = 3 1 0 5 1 6 50 53 = 3 1 0 5 327 354 = 4000 4609
7 5
439
31
1
7
1 0 2 3 3
[ A ] ([ B ] [C ]) = 3 1 0 5 1
1
0
=
2
3
1
.
0
([ A ] [ B ])
= [B ] [ A]
(2.9)
La matrice identit, note [ I ] , est une matrice carre dont les termes diagonaux
sont gaux 1, tous les autres tant nuls. De ce fait, le produit de la matrice identit par une matrice [ A ] quelconque (ou inversement) est gal la matrice [ A ]
elle-mme.
9782100544639.indb 7
1
0
[ I ] = 0
.
0
0
1
0
.
0
0
0
1
.
0
.
.
.
.
.
0
0
0 et [ I ] [ A ] = [ A ] [ I ] = [ A ]
.
1
(2.10)
09/02/10 14:53
[ A] [ A] = [ A] [ A]
Calcul du dterminant
2 dimensions:
a11
det
a21
a12 a11
=
a22 a21
a12
= a11 a22 a12 a21
a22
(2.12)
3 dimensions:
a11
det a21
a31
a12
a22
a32
a13 a11
a23 = a21
a33 a31
a12
a22
a32
a13
a23
a33
(2.13)
a12
a22
a32
a13
a23
a33
a11
a21
a12
a22
a13
a33
9782100544639.indb 8
09/02/10 14:53
2 dimensions:
21
[ A]
a12
a22
Com [ A ] =
a22
a12
a22
1
=
det [ A ] a21
a12
a11
a21
a11
(2.15)
3 dimensions:
a21 a23
a21 a22
a22 a23
+
+
a31 a33
a31 a32
a322 a33
a11 a12 a13
a12 a13
a11 a12
a11 a13
[ A ] = a21 a22 a23 Com [ A ] = a a + a a a a
32 33
31
33
31
32
a31 a32 a33
a21 a23
a21 a22
a22 a23
a22 a33 a23 a32 a23 a31 a21a33 a21a32 a22 a31
1
1
a13 a32 a12 a33 a11a33 a13 a31 a12 a31 a11a32
[ A] =
det [ A ]
a12 a23 a13 a22 a13 a21 a11a23 a11a22 a12 a21
(2.16)
a11
[ A] = a
(2.17)
[I ]
9782100544639.indb 9
09/02/10 14:53
lastiques comportant plusieurs cas de charges, leur utilisation apparat moins pertinente dans la mesure o celles-ci ncessitent une rsolution complte chaque
changement dtat de charges.
En effet et dans le cas des mthodes directes appliques aux problmes linaires
lastiques, seul le premier cas de charges ncessite une inversion de la matrice de
rigidit, les rsultats des cas suivants tant obtenus par linarit aprs stockage dela
matrice inverse en mmoire (uniquement si les conditions dappui ou de temprature ne varient pas).
Rsolution par la mthode par limination de Gauss
.
.
.
.
.
.
.
.
.
.
(2.18)
q1 =
1
( F1 k12 q2 k13q3 ................ k1n qn )
k11
(2.19)
avec
1 q + k1 q + ................. + k1 q = F 1
k22
2
23 3
2n n
2
.
.
.
.
.
.
.
.
.
.
.
.
1
1
1
kn 2 q2 + kn 3 q3 + ................. + knn qn = Fn1
(2.20)
k s 1ksjs1
s
s 1 is
k
=
k
ij
ij
ksss1
s 1 F s 1
k
s
is
F s = F s 1
i
i
ksss1
9782100544639.indb 10
09/02/10 14:53
1 q + k1 q + ................. + k1 q = F 1
k22
2
23 3
2n n
2
2
2
2
.
.
.
.
n
1q = F n 1 q
knn
n
n
n
(2.21)
3 0 1 1 0 U 2 0
0 2 1 0 0 V 0
2
ES
1 1 2 1 1 V3 = P
2L
1 0 1 3 0 U 4 0
0 0 1 0 2 V4 0
(2.22)
U2 =
V3 + U 4
3
(2.23)
1
1
ES 0 1 2
1
3
3
2L
1
1
0 0 1
3
3
3
0 0
1
0
9782100544639.indb 11
0
3 0 1 1 0
0 2 1 0 0 U 2 0
U
2
0
0
V 2
V2
5 2
1 V P
1 V ES 0 1
3=
3=
3 3
2L
U
0 0 2 8 0 U 4 0
0 4
V
V 0
3 3
4
0 0 1 0 2 4
2
(2.24)
11
09/02/10 14:53
2e tape: limination de V2
Pour liminer V2 , on rpte lopration en considrant cette fois la deuxime quation de (2.24) do:
V
V2 = 3
(2.25)
2
Le systme devient alors:
3
0
ES 0
2L
0
1 0
3 0 1 1 0
0 2 1 0 0 U 2 0
0 0 U 2
V2
V 2
7 2
2
1 V P
1 V ES 0 0
3=
3=
6 3
3
2L
U
8
0 0 2 8 0 U 4 0
0 4
V
V 0
3 3
3
4
0 0 1 0 2 4
0 2
0
2
1
1
5 1
0
3 2
2
0
3
0
1
(2.26)
Et ainsi de suite
3e tape: limination de V3
6 2 PL 2
V3 =
U 4 + V4
7 ES 3
3
0
0
ES
2L
0
3
0
ES 0
2L
0
0 1
2 1
7
0
6
0
0
0
2
0
0
0
0
1
0
0
0 U
2
0
2
1 V2 P
3
V3 =
8 2 6 2 2 6 2 6 P
0
+
U4
3 3 7 3 3 7 3 7
V 6
6 2
6 4 P
0
2
7
7 3
7
1 1 0
0
1 0 0 U 2
7 2
1 V2 P
6 3
V3 = 4
16 4 P
U4
0
7
7 7
V4 6
4 8
7 P
0
7 7
(2.27)
(2.28)
12
9782100544639.indb 12
09/02/10 14:53
4e tape: limination de U4
U4 =
0 1 1
2 1 0
7 2
0
6 3
16
0 0
7
0
(2.29)
0
0
V
2
1
P
V
=
3
4
4
P
U
7
4
7
V4 6
1
8 4 1
7 P 7 P
+
7 7 4
(2.30)
3
0
0
ES
2L
0
PL V4
2 ES 4
Soit finalement:
3
0
ES 0
2L
0
0 1 1
2 1 0
7 2
0
6 3
16
0 0
7
0
0
0
0 U 2
V 2
P
1
V
=
3
4
4
P
U4
7
7
1
V4 6
8 1
7 P 7 P
7 7
(2.31)
On dduit de la 5e quation:
ES
2 PL
V4 = P V4 =
2L
ES
(2.32)
Do partir de (2.23), (2.25), (2.27) et (2.29), les valeurs des autres inconnues:
PL V4 PL
=
2 ES 4 ES
6 2 PL 2
4 PL
V3 =
U 4 + V4 =
ES
7 ES 3
V3
2 PL
V2 =
=
ES
2
V3 + U 4
PL
U2 =
=
ES
3
U4 =
9782100544639.indb 13
(2.33)
13
09/02/10 14:53
Bien videmment, la mthode par limination de Gauss savre dans ce cas bien
complique compare une approche plus classique telle que celle dcrite ci-dessous.
Reprenant le systme (2.22), les 2e et 5e quations permettent de dduire respectivement que 2V2 = V3 et 2V4 = V3 do V2 = V4 .
De la somme des 1re et 4e quations rsulte que 2U 2 + 2U 4 = 0 U 2 = U 4 .
Do partir de la 1re: 3U 2 V3 U 4 = 0 V3 = 4U 2 .
En remplaant ces diffrents rsultats dans la 3e quation, on obtient finalement:
2 PL
2 PL
4 PL
2U 2 + V3 =
V3 =
ES
ES
ES
2 PL
PL
.
U 2 = U 4 =
et V2 = V4 =
ES
ES
U 2 V2 + 2V3 + U 4 V4 =
Soit
14
9782100544639.indb 14
09/02/10 14:53
A
3.1 Les contraintes
3.1.1 Notion de contrainte
Son tat dquilibre permet dtablir que la somme vectorielle des forces issues
de
1 et agissant
sectionS.
dF
On dfinit donc le vecteur contrainte f1 comme tant la limite de
lorsque la
dS
surface dS tend vers zro. Quand dS est considre dans le plan yz, ce vecteur f1
peut tre dcompos en trois composantes: une normale cette surface et deux
dF
dans son plan. La premire, gale xx = x = lim x , est appele contrainte
dS 0 dS
dFy
normale alors que les deux autres, notes respectivement xy = lim
et
dS 0 dS
dF
xz = lim z , sont dites de cisaillement.
dS 0 dS
9782100544639.indb 15
15
09/02/10 14:53
xy
xz
dF
j
z
dS
i
k
( )
( )
xx yx zx
f1 = xy ; f 2 = yy ; f 3 = zy
(3.1)
xz
zz
yz
La convention de signes la plus souvent retenue dans les logiciels lments finis est
dassocier une contrainte normale positive une traction.
3.1.2 quations internes
Laxe z
M /Oz = 0
xy
yx
dy
dx
dx
dy
dx dydz + xy dydz yx +
dy dxdz yx dxdz = 0
xy +
x
2
y
2
M /Oz = 0
2
2
xy
yx
dy
dx
dx
dy
+
dx dydz + xy dydz yx +
dy dxdz yx dxdz = 0
2
x
y
2
2
2
xy dx yx dy
xy yx +
=0
y 2
x
2
16
9782100544639.indb 16
09/02/10 14:53
yy = yy +
xz
zy zy
fz v
zz
yx
dy
yx
yx = yx +
xy = xy +
xy
x
zz
zx
dy
xy
fyv
yz
xx
yy
fxv
O zy
zx
yz
yx
y
dy
dx
xx = xx +
xx
dx
x
xz
dz
y
x
dx
yy
zx = zx +
zx
dz
z
do:
xy = yx
(3.2)
yz = zy et zx = xz
(3.3)
{ }T = { xx yy zz xy yz xz }
(3.4)
F/ x = 0 xx +
yx
xx
dx dydz + yx +
dy dxdz
x
y
xx yx zx
+
+ f xv = 0
+
y
z
x
9782100544639.indb 17
17
09/02/10 14:53
(3.5)
zz
xz
zx
xx
xy
zy
yz
yx
yy
C
z
Figure 3.4 quilibre de surface.
Soit un ttradre OABC infiniment petit dont laire ABC, note ds , constitue la
surface du milieu solide.
f xs
Sachant quune charge rpartie f s de composantes f ys est applique sur cette
f s
z
18
9782100544639.indb 18
09/02/10 14:53
surface et que les cosinus directeurs de la normale au plan ABC sont nots l, m etn,
lquilibre du ttradre se rsume un systme de trois quations:
f xs ds = l ds xx + m ds yx + n ds zx
f ys ds = l ds xy + m ds yy + n ds zy
f zs ds = l ds xz + m ds yz + n ds zz
f xs xx
f xs = l xx + m yx + n zx
f ys = l xy + m yy + n zy f ys = xy
f s
f zs = l xz + m yz + n zz
z xz
yx
yy
yz
l
zx l
zy m = [ ] m (3.6)
n
zz n
Laction de forces extrieures sur un solide dformable entrane pour chacun des
points de sa gomtrie un mouvement pouvant tre dcompos en trois phases
distinctes: une translation densemble (mouvement de corps rigide en translation),
une rotation densemble (mouvement de corps rigide en rotation) et enfin une
dformation pure.
La dformation qui correspond un mouvement relatif par rapport un rfrentiel donn, peut suivre plusieurs hypothses comme celles des petites ou grandes
dformations mais galement tre associe des petits ou grands dplacements.
Considrant le cas le plus courant, seules les hypothses de petites dformations
et petits dplacements seront dveloppes dans ce qui suit, les deux restantes tant
abordes au chapitre 11.
3.2.2 Dplacements
q = u i + v j + w k
(3.7)
On dfinit par ailleurs la norme du dplacement comme tant:
q = u2 + v 2 + w 2
(3.8)
9782100544639.indb 19
19
09/02/10 14:53
q
w
z
Figure 3.5 Dplacement dun point.
v u
+
xy = lim B AD =
dx 0
2
x y
dy 0
(3.10)
La gnralisation de cette approche bidimensionnelle aux deux autres plans permettra dobtenir:
u
x
v
=
y
w
=
z
v u
+
x y
v w
=
+
z y
u w
=
+
z x
xx =
xy =
yy
yz
zz
zx
(3.11)
20
9782100544639.indb 20
09/02/10 14:53
(3.12)
{ } = xx yy zz xy yz xz
Lexpression (3.11) est en fait une version linarise de celle tablie lors de ltude
de la non linarit gomtrique au chapitre 11 (cf. quations (11.3) et (11.4)).
Par ailleurs, (3.11) pourra galement tre exprime sous la forme matricielle suivante:
u+
z
0
(3.13)
C
D
C
B
v
x
B
A
u
dy
y
D
dy
u
z
v
0 w
xy
2
u
y
v
v + dy
y
xx 0
yy 0
zz
{ } = =
xy
yz y
xz 0
u
dx
u+
v+
v
dx
x
u
dx
x
x
O
Figure 3.6 Dformations planes.
9782100544639.indb 21
21
09/02/10 14:53
Secon S
x
y
F
Figure 3.7 Allongement dun barreau.
u
e
xx =
F
S
xx =
E=
xx
xx
Zone linaire
lasque
Zone
plasque
Raffermissement
L
L
Rupture
22
9782100544639.indb 22
09/02/10 14:53
avec:
e et u limites dlasticit et de rupture du matriau,
E: module dlasticit longitudinal ou module de Young (pour lacier E =
2.11011N/m2).
La zone lastique est caractrise par la loi de Hooke:
xx = E xx
(3.14)
liant linairement contraintes et dformations. En plus de la dformation axiale
et toujours dans le domaine lastique, une dformation associe la contraction
transversale du barreauest observe la fois suivant y:
h
= xx = xx
h
E
(3.15)
zz =
= xx = xx
d
E
(3.16)
yy =
et suivant z:
n qui correspond au coefficient de Poisson doit tre, tout en tant positif, stricte-
9782100544639.indb 23
23
09/02/10 14:53
1
xx yy + zz
E
1
yy = yy ( xx + zz )
(3.17)
E
1
zz = zz xx + yy
E
De plus, il est galement possible via un essai de cisaillement pur de montrer quil
existe des relations linaires entre contraintes et dformations de cisaillement telles
que:
xy
yz
E
xy =
; yz =
; xz = xz avec G =
(3.18)
2 (1 + )
G
G
G
xx =
1
E
=
D
D
avec
{} [ ]{ }
[ ] E
E
1
E
E
1
E
0
0
1
G
1
G
(3.19)
(1 + ) (1 2 ) xx
E
(1 ) + ( xx + zz )
yy =
(1 + ) (1 2 ) yy
E
(1 ) + xx + yy
zz =
(1 + ) (1 2 ) zz
xy = G xy
xx =
(3.20)
yz = G yz
xz = G xz
24
9782100544639.indb 24
09/02/10 14:53
{ } = [ H ] { }
E
0
[H ] =
(1 + ) (1 2 )
0
0
ce qui permet de dmontrer que
1
0
0
0
[H ]
0
0
0
1 2
2
(3.21)
0
0
0
0
0
0
0
0
0 = [ D ]1
1 2
0
0
0
2
1 2
0
0
0
2
nest pas dfinie pour = 0.5 .
avec
Considrant un solide de section transversale constante, lhypothse de dformation plane consiste considrer un comportement plan identique quelle que soit
la profondeur. En dautres termes, le solide peut tre dcompos en tranches
dpaisseur unitaire, le calcul se limitant ltude dune seule de ces tranches.
y
yy
zz = 0
xx
z
Figure 3.10 tat plan de dformation.
9782100544639.indb 25
25
09/02/10 14:53
En consquence, toutes les dformations associes laxe transversal, z en loccurrence, seront prises gales zro.
On a donc:
zz = xz = yz = 0 et xz = yz = 0
(3.22)
1
xx
E
1
yy =
(1 + ) (1 2 )
xy
0
0
0 xx
0 yy
1 2 xy
(3.23)
zz = ( xx + yy )
(3.24)
Contrainte plane
zz = xz = yz = 0 et xz = yz = 0
(3.25)
zz = 0
xx
yy
y
Figure 3.11 tat plan de contrainte.
zz =
( xx + yy )
E
(3.26)
26
9782100544639.indb 26
09/02/10 14:53
1
0 xx
0 yy
1
1 xy
0 0
(3.27)
xx
E
yy =
2
1
xy
F
F
dW
dy
y
x
dz
dx (1+ x)
u
du
9782100544639.indb 27
27
09/02/10 14:53
1
1
xx xx dxdydz = xx xx dV
2V
2V
(3.30)
Lextension du cas uniaxial prcdent aux deux autres axes y et z permettra dobtenir par superposition:
W =
1
xx xx + yy yy + zz zz dV
2 V
(3.31)
xy
xy dx
dy
xy
xy
Figure 3.13 nergie de dformation de cisaillement.
dW =
1
1
xy dydz xy dx = xy xy dxdydz
2
2
(3.32)
1
xx xx + yy yy + zz zz + xy xy + yz yz + xz xz dV (3.33)
2 V
W =
T
1
{ } { } dV
2V
(3.34)
28
9782100544639.indb 28
09/02/10 14:53
1
2 + 2 + 2 dV
xx
yy
zz
2 E V
(3.35)
1
( xx yy + yy zz + xx zz ) dV +
( xy2 + 2yz + xz2 ) dV
EV
2 G V
9782100544639.indb 29
29
09/02/10 14:53
9782100544639.indb 30
09/02/10 14:53
4 Principes de la mthode
deslmentsfinisenstatique
lments
Nuds
Figure 4.1 (a) Solide (Poutre en I); (b) Modle lments finis.
9782100544639.indb 31
31
09/02/10 14:53
4 Principes de la mthode
deslmentsfinisenstatique
X ( x, y, z ) = N i ( x, y, z ) X i
i =1
= N1 ( x , y , z ) X 1 + N 2 ( x , y , z ) X 2 + ..... + N n ( x , y , z ) X n
(4.1)
o les:
X i dsignent les inconnues en dplacement associes aux degrs de libert
considrs,
N i ( x , y , z ) dsignent les fonctions de forme de llment permettant dobtenir
les dplacements de celui-ci en nimporte quel point de sa gomtrie et ce toujours partir des dplacements nodaux.
titre dexemple, prenons un lment barre travaillant uniquement en traction
ou en compression. Le matriau utilis tant suppos linaire, la dformation est
de fait constante sur la hauteur de sa section. La fonction de dplacement u( x ) est
alors forcment linaire. On a donc:
xx =
du
= Cte u( x ) = a0 + a1 x
dx
(4.2)
x
x
x
u( x ) = (1 ) u1 + u2 = 1
L
L
L
soit sous une autre forme(figure 4.2):
N1
N2
x
1
u1
u (x)
u2
L
Figure 4.2 Fonctions de forme dun lment barre.
32
9782100544639.indb 32
09/02/10 14:53
Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
2
x
x
u( x ) = (1 ) u1 + u2 = N1 u1 + N 2 u2 = N i ui = [ N ] {qe } (4.3)
L
L
i =1
{u} = [ N ] {qe }
4.2 Rsolution
E = W T =
T
1
{ } { } dV {u}T { f v } dV {u}T { f s } dS (4.5)
2
V
S
V
W
9782100544639.indb 33
33
09/02/10 14:54
4 Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
1
x
L
W
u
Y
X
F
Figure 4.3 Barreau en traction.
Posant que xx =
du u2 u1 u 0
F
=
=
et xx = , la relation contrainte-dfordx
L
L
S
mation scrit:
xx = E xx
F
u
ES
=E F =
u = k u
S
L
L
(4.6)
ES
correspondant dans le domaine lastique la rigidit du barreau et
L
la pente de la droite F = k u . La dformation xx tant indpendante de x, lnergie de dformation devient:
avec k =
2
1
1
2 dS dx = E xx dS dx
W = xx xx dV = E xx
2V
2V
2 S
0
S
(4.7)
E S L u
1
= k u2
2
2
L
E = W T =
1
k u2 F u
2
(4.8)
34
9782100544639.indb 34
09/02/10 14:54
Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
Notant la relation (3.13) sous une forme plus synthtique { } = [ ]{u} , le vecteur des dformations pourra tre exprim en fonction des dplacements nodaux
partir de:
][ N ] {qe } = [ B ] {qe }
{ } = [ ]{u} = [
(4.10)
[B ]
avec [ ] : oprateur drivation dont les dimensions dpendent de celles de llment tudi.
Soit pour llment barre: xx =
du 1 1 u1
=
. Le champ de contrainte
u
dx
L
L
2
[B ]
T
T
1
1
{ } { } dVe = {qe } [ B ]T [ H ] [ B ] {qe } dVe
2V
2V
e
1
T
= {qe } [ ke ] {qe }
2
(4.11)
[ke ] = [ B ]T [ H ] [ B ] dVe
avec
(4.12)
Ve
[ke ] = [ B ] [ H ] [ B ] dVe = 0
Ve
9782100544639.indb 35
1
L L2
[ke ] = ES 0 1
L2
1
L
1
dS
e 1 E L
s
L
S
1
2
L dx = ES 1 1
1
L 1 1
L2
1
dx
L
(4.13)
35
09/02/10 14:54
4 Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
En traduisant les actions des lments connects e sous la forme dun vecteur
forces { pe } , Te scrit:
Te =
}T {p
e}
{u} { f v } dVe + {u} { f s } dSe + {qe
Ve
Se
(4.14)
Te =
T
T
T
{qe } [ N ] { f v } dVe + {qe } [ N ] { f s } dSe + {qe } { pe } (4.15)
Ve
Se
Do finalement:
T
Te = {qe }
avec:
({ f ev } + { f es } + { pe })
(4.16)
(4.17)
{ f es } = [ N ]T { f s } dSe
(4.18)
Ve
Se
Les caractristiques globales sont obtenues par assemblage des diffrentes quantits
We et Te tablies en (4.11) et (4.16). Lnergie potentielle totale de la structure est
alors gale la somme des diffrences entre nergies de dformation et travaux des
forces de chacun des lments:
E = W T = We Te
1
T
T
= {qe } [ ke ] {qe } {qe } ({ f ev } + { f es } + { pe })
e 2
e
(4.19)
Changement de repres
Cependant, les lments nayant pas forcment le mme repre local, la sommation
des nergies de dformation et des travaux des forces extrieures ne peut tre directement applique. Il est donc ncessaire de les exprimer dans un rfrentiel unique
appel repre global. Pour ce faire, un changement de base doit tre effectu du
repre global vers le repre local de chacun des lments.
Quand on se limite des problmes plans, les formules de changement de base
permettant de passer du repre global au repre local scrivent:
x = X cos ( ) + Y sin ( )
y = X sin ( ) + Y cos ( )
(4.20)
36
9782100544639.indb 36
09/02/10 14:54
Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
Y
y
titre dexemple, le changement de base dun lment barre plan deux nuds
revient appliquer la relation (4.20) ses deux nuds i et j, soit:
0
0 U i
ui cos ( ) sin ( )
v sin cos
0
0 Vi
( )
( )
i
(4.21)
=
u
0
cos ( ) sin ( ) U j
j 0
v j 0
0
sin ( ) cos ( ) V j
Dune manire gnrale, ce changement de base permet de calculer les dplacements des nuds de llment e en repre local {qe } partir de ceux du repre
global {Qe } via une matrice de passage nomme [ Re ] .
{qe } = [ Re ] {Qe }
(4.22)
Vi
Dunod La photocopie non autorise est un dlit.
Vj
u (x)
ui
uj
Uj
Ui
9782100544639.indb 37
(4.23)
37
09/02/10 14:54
4 Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
1
{Qe }T [ K e ] {Qe }
2
Do lexpression de la rigidit lmentaire en repre global:
We =
Soit pour llment barre:
(4.24)
[ K e ] = [ Re ]T [ke ] [ Re ]
cos2
ES cos sin
[K e ] =
cos2
L
cos sin
cos sin
sin 2
cos sin
sin 2
(4.25)
cos2
cos sin
cos2
cos sin
ES 1 1 (ui ) ES 0
=
avec daprs (4.13) [ ke ] =
L 1 1 (u j ) L 1
cos sin
sin 2
(4.26)
cos sin
sin 2
0 1
0 0
0 1
0 0
0 (ui )
0 (vi )
0 (u j )
0 (v j )
0
0
cos ( ) sin ( )
0
0
sin ( ) cos ( )
et (4-21) [ Re ] =
0
0
cos ( ) sin ( )
0
sin ( ) cos ( )
0
Ce qui donne pour le barreau du chapitre 4.2.1 (nud 1 vers 2 =-90 ):
0 1
1 0
[ R1 ] = 0 0
0 0
0 0
0 0
et
0 1
1 0
0 0
ES 0 1
[ K 1 ] = 0 0
L
0 1
0 0 (U1 )
0 1 (V1 )
0 0 (U 2 )
0 1 (V2 )
En appliquant (4.22) aux vecteurs charges, le travail lmentaire des forces (4.16)
scrit:
T
Te = {qe }
= {Qe }
= {Qe }
({ f ev } + { f es } + { pe }) = ([ Re ] {Qe })T ({ f ev } + { f es } + { pe })
T
T
T
v
s
Re ] { f e } + [ Re ] { f e } + [ Re ] { pe }
[
{Pe }
{Fev }
{Fes }
(4.27)
38
9782100544639.indb 38
09/02/10 14:54
Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
1
T
T
= {Qe } [ K e ] {Qe } {Qe }
2
e
(4.28)
{Qe } = [ Ae ] {Q }
Exemple: soit deux lments barres de connectivit 1, 2 et 3
Les matrices de rigidit et vecteurs forces des diffrents lments tant exprims
dans le repre global, il reste nanmoins positionner chacun dentre eux dans une
matrice de rigidit [ K ] et un vecteur charges {F } caractrisant lensemble de la
structure.
Considrant une structure n nuds, le vecteur global des dplacements {Q }
aura un nombre de lignes gal n fois le nombre de degrs de libert par nud.
Il en est dailleurs de mme pour le vecteur {F } . Pour ce qui concerne la matrice
[ K ] , matrice carre de mmes dimensions, lassemblage seffectuera en positionnant chacune des matrices de rigidit lmentaires [ K e ] en fonction des degrs de
libert des nuds de llment considr.
Ce positionnement est mis en uvre au niveau informatique via une matrice boolenne dassemblage [ Ae ] qui permet de placer les degrs de libert associs un
lment donn dans le vecteur global des dplacements {Q } . On a donc:
(4.29)
Y
1
U1
U 2
V
V
1
2
{Q1} = , {Q2 } =
U 2
U 3
V2
V3
et du vecteur global
9782100544639.indb 39
U1
V
1
U 2
{Q } = ,
V2
U 3
V3
39
09/02/10 14:54
4 Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
on dduit pour [ A1 ] et [ A2 ] :
U1
U1 1 0 0 0 0 0 V1
V 0 1 0 0 0 0 U
1
2
{Q1} = =
U 2 0 0 1 0 0 0 V2
V2 0 0 0 1 0 0 U 3
[ A1 ]
V3
U1
U 2 0 0 1 0 0 0 V1
V 0 0 0 1 0 0 U
2
2
{Q2 } = =
U 3 0 0 0 0 1 0 V2
V3 0 0 0 0 0 1 U 3
[ A2 ]
V3
En introduisant (4.29) dans lexpression de lnergie potentielle (4.28), celle-ci
devient:
E = W T
1
T
T
T
T
= {Q } [ Ae ] [ K e ][ Ae ]{Q } {Q } [ Ae ]
e 2
E=
(4.30)
1
{Q }T [ Ae ]T [ K e ] [ Ae ] {Q }
2
e
T
T
{Q } [ Ae ]
e
T
Sachant que les termes [ Ae ] [ K e ] [ Ae ] et [ Ae ] ({Fev } + {Fes } + {Pe }) correspondent respectivement aux positionnements de [ K e ] dans [ K ] et des {Fe } dans
{F } , on obtient alors:
1
T
T
E = W T = We Te = {Q } [ K ] {Q } {Q } {F } (4.32)
2
e
e
avec:
T
[ Ae ] [ K e ] [ Ae ] = [ K ] : matrice de rigidit de la structure.
e
T
[ Ae ] {Pe } = {P } : vecteur des charges nodales concentres dfinies dans le
e
repre global.
40
9782100544639.indb 40
09/02/10 14:54
Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
Comme lquilibre est caractris par une variation nulle dnergie potentielle par
rapport au dplacement, le systme scrit finalement:
= {dQ }
(4.33)
([ K ] {Q } {F }) = 0 [ K ] {Q } = {F }
dE
T
T
= {dQ } [ K ] {Q } {dQ } {F }
dQ
Pour fixer la structure au sol, des conditions dappui sont introduites de telle
manire assurer sa stabilit densemble.
Suivant la nature des degrs de libert bloqus, ces supports peuvent reconstituer
des appuis ponctuels, des appuis rotuls ou des encastrements mais galement correspondre des dplacements nuls1 ou imposs2.
Regroupant les dplacements connus3 lis aux appuis au sein dun vecteur {Q B }
et ceux inconnus dans {Q L } , le systme [ K ] {Q } = {F } peut tre transform de
telle manire sparer les degrs de libert libres et bloqussoit:
[ K ] [ K LB ] {Q L } {Fext }
[ K ] {Q } = {F } LL
=
[ K BL ] [ K BB ] {Q B } {R }
(4.34)
(4.35)
Avec:
{R } : vecteur des ractions dappui,
{Fext } : vecteur des charges associes aux degrs de libert libres.
Les dplacements inconnus {Q L } peuvent donc tre dtermins en rsolvant le
systme linaire:
Pour tre inversible et donc ne pas tre singulire, [ K LL ] ne devra pas comporter
de termes diagonaux nuls (cf. 2.3.3). Ceux-ci ne pourront pas non plus tre
T
ngatifs du fait de son caractre dfinie positive (i.e. {Q L } [ K LL ] {Q L } > 0 ).
Dans le cas contraire, on parlera de pivot ngatif ou nul, terme qui dailleurs est
largement repris par les logiciels lorsque le calcul est interrompu pour ces raisons.
Une fois {Q L } dtermin, les ractions dappui peuvent tre calcules partir de:
{R } = [ K BL ] {QL } + [ K BB ] {QB }
(4.36)
9782100544639.indb 41
41
09/02/10 14:54
4 Principes de la mthode
deslmentsfinisenstatique
4.2 Rsolution
ES 0 1
[ K ] {Q } = [ K 1 ] {Q } = 0 0
L
0 1
0 0 U1 RH1
0 1 V1 RV1
=
0 0 U 2 0
0 1 V2 F
do
RH1 = 0
U 2 = 0
et
FL
ES
V2 = ES
RV1 = L V2 = F
On notera par ailleurs que la force F est une force nodale et que le dplacement
V2 est ngatif car exprim dans le repre global.
u1 0 1
v 1 0
1
De plus et comme {q1} = [ R1 ] {Q1 } =
u2 0 0
v2 0 0
dduit:
FL
u2 =
puisque u2 = V2
ES
0 0 U1
0 0 V1
, on
0 1 U 2
1 0 V2
Une fois les dplacements en repre global connus et pour dterminer les efforts
internes aux nuds, il suffit de calculer pour chacun des lments, les dplacements associs en repre local puis de les injecter dans la relation de rigidit lmentaire. On a donc:
(4.37)
=
L 1 1 u2 N1 F
42
9782100544639.indb 42
09/02/10 14:54
Principes de la mthode
deslmentsfinisenstatique
Construction de la matrice
de rigidit [k e ] (repre local)
m lments
e=1
Construction du vecteur
de charges [ f e ] (repre local)
e=1
Calcul de la matrice de passage
m lments
Assemblage de [K ]
Assemblage de {F } en prenant en compte
les ventuelles charges nodales
Dunod La photocopie non autorise est un dlit.
m lments
e=1
Calcul des efforts internes en repre
local par la relation [ke ] {qe }
FIN
9782100544639.indb 43
43
09/02/10 14:54
9782100544639.indb 44
09/02/10 14:54
A
5.1 lment de barre
Y
Vj
E, S
Vi
i
ui
uj
Uj
Ui
9782100544639.indb 45
45
09/02/10 14:54
On a alors:
1
ES 1 1 (ui ) ES 0
[ke ] = 1 1 (u ) = 1
L
L
j
0 1
0 0
0 1
0 0
0 (ui )
0 (vi )
0 (u j )
0 (v j )
(5.1)
[ K e ] = [ Re ]
cos
sin
[ ke ] [ Re ] avec [ Re ] =
0
sin
cos
0
0
0
0
cos
sin
0
0
sin
cos
(5.2)
cos2
ES cos sin
[K e ] =
cos2
L
cos sin
cos sin
sin 2
cos sin
sin 2
cos2
cos sin
cos2
cos sin
cos sin
sin 2
(5.3)
cos sin
sin 2
Vj
Vi u
i
i
uj
Uj
Ui
9782100544639.indb 46
09/02/10 14:54
1
0
1 1 (ui )
[ke ] = k 1 1 (u ) = k 1
0 1
0 0
0 1
0 0
0 (ui )
0 (vi )
0 (u j )
0 (v j )
(5.4)
y
x
E, S, L
k=
4
E,
45
ES
L
2S , 2 L
Y
x
1
Nud i
Nud j
1 (barre)
2L
2S
2 (barre)
3 (ressort)
Application numrique :
9782100544639.indb 47
= 10 m,
Longueur
= 0.0001 m2,
Section
= 10000 N
47
09/02/10 14:54
E 2S 0
[k1 ] =
2L 1
0
1
ES 0
k
=
[ 2]
L 1
0 1
0 0
0 1
0 0
0 1
0 0
0 1
0 0
0 (u1 )
0 (v1 )
0 (u3 )
0 (v3 )
(5.5)
0 (u2 )
0 (v2 )
0 (u3 )
0 (v3 )
Llment 2: nuds 2 3, q = 0
ES 0
[ K 2 ] = 1
L
0 1
0 0
0 1
0 0
0 (U 2 )
0 (V2 )
0 (U 3 )
0 (V3 )
(5.6)
Llment 1: nuds 1 3, q = 45
1
2
1
ES
[ K 1 ] = 21
L
2
1
2
1
2
1
2
1
2
1
1
2
1
2
1
2
1
2
1
2
1 (U1 )
(V )
2 1
1 (U 3 )
2 (V3 )
1
2
(5.7)
9782100544639.indb 48
09/02/10 14:54
1
2
1
2
ES 0
[ K ] {Q } = {F } 0
L
1
2
1
1
2
1
2
0
0
1
2
1
1
0
0
0
1 0
0
1
2
1
2
1
0
1
1+
2
1
2
1
2
1 U1 RH1
2 V1 RV1
0 U 2 RH 2
=
(5.8)
0 V2 RV2
1 U 3 0
2 V3 F
1
Dans la mesure o les appuis sont infiniment rigides et comme indiqu en 4.2.5,
lintroduction des conditions aux limites aux nuds 1 et 2 revient barrer la ligne
et la colonne du degr de libert considr.
9782100544639.indb 49
1
1+
ES 2
L 1
2
1
2 U 3 = ES 3 1 U 3 = 0
1 V3 2 L 1 1 V3 F
2
(5.9)
49
09/02/10 14:54
FL
U 3 = ES = 4.762 mm
V3 = 3 FL = 14.286 mm
ES
do
(5.10)
Il suffit alors dinjecter les rsultats (5.10) dans (5.8) pour obtenir les ractions aux
appuis:
RV1 = RH1 =
ES U 3 V3
ES
= F ; RH 2 =
( U 3 ) = F ; RV2 = 0
2
L
L 2
(5.11)
N2
E, S, L
N1
N2
2
4
E,
45
k=
ES
L
2S , 2 L
Y
x
N1
La dtermination des efforts dans les barres est base sur lapplication de la relation
(4.37). Il est nanmoins ncessaire de calculer pralablement les dplacements dans
le repre local de chacun des lments. Ainsi, nous avons pour:
Llment 1: nuds 1 3, q = 45
1
u1
v
2 1
1
Comme daprs (4.22) {q1} = [ R1 ] {Q1 } =
u3 2 0
v3
0
1 0
1 0
0 1
0 1
0 U1
0 V1
,
1 U 3
1 V3
50
9782100544639.indb 50
09/02/10 14:54
E 2S 0
=
2 L 1
0 1
0 0
0 1
0 0
0
1
0 2 1
0 2 0
0
0
1 0
1 0
0 1
0 1
0 U1
0 V1
=
1 U 3
1 V3
2P
2P
0
(5.12)
N1
0
= [ k1 ] {q1} = [ k1 ] [ R1 ] {Q1 }
N1
0
ES 2
(U 3 + V3 ) = 2 P = 14142 N
L 2
N1
xx1 =
= 100 MPa
2S
N1 =
(5.13)
Comme indiqu prcdemment, lhypothse retenue pour les efforts normaux est
dassocier une valeur positive une traction. Cela revient dans notre cas retenir
comme convention, laction du nud sur la barre do le signe appliqu
leffort normal du nud de dpart.
Llment 2: nuds 2 3, q = 0
Les repres local et global tant confondus, il est possible de dterminer directement leffort dans la barre n2 car [ R2 ] = [ I ] (i.e. u3 = U 3 ). On a donc:
ES 1 1 u2 N 2 F
=
=
L 1 1 u3 N 2 F
N2
N 2 = F = 10000 N xx2 =
= 100 MPa
S
(5.14)
Llment ressort tant orient de4 vers 3, on peut crire daprs (5.4) que:
9782100544639.indb 51
ES 0
[k3 ] = 1
L
k
0
0 1
0 0
0 1
0 0
0 (u4 )
0 (v4 )
0 (u3 )
0 (v3 )
(5.15)
51
09/02/10 14:54
ES 0 1
[ K 3 ] = 0 0
L
0 1
0 0 (U 4 )
0 1 (V4 )
0 0 (U 3 )
0 1 (V3 )
(5.16)
1
1
1 + 2
2
1
1
1+
2
2
0
0
0
ES 0
1
L 1
2
2
1
1
2
2
0
0
0
0
1
0
0
0
1 0
0
0
0
0
0
1
1
0 0
2
2
U1 RH 1
1
1
0 0 V R
1 V1
2
2
1
0
0 0 U 2 RH 2
0
0
0 0 V2 Rv 2
U = 0 (5.17)
1
1
1+
0 0 3
2
2
V3 F
1
1
1+
0 1 U 4 RH 4
2
2
V R
0
0
1 0 4 V 4
1 0 1
0
1
1
1+
ES
2
2 U 3 = ES 3 1 U 3 = 0
1 V
2 L 1 3 V3 F
L 1
1+ 3
2
2
(5.18)
FL
U 3 = 4 ES = 1.190 mm
(5.19)
do
3
FL
V3 =
= 3.571mm
4 ES
En effet, la seule diffrence se situe au niveau de la rigidit associe V3 qui comES
supplmentaire d la prsence du ressort.
porte un terme
L
52
9782100544639.indb 52
09/02/10 14:54
Les ractions peuvent alors tre calcules en injectant (5.19) dans (5.17)do:
ES U 3 V3 F
RV1 = RH1 = L 2 2 = 4
R = ES U = F
( 3)
H
2
4
L
RV = RH = 0
4
2
S
3F
E
R =
V3 ) =
(
V
4
L
4
(5.20)
Llment 1: nuds 1 3, q = 45
1
u1
v
2 1
1
Comme {q1} = [ R1 ] {Q1 } =
u3 2 0
v3
0
9782100544639.indb 53
1 0
1 0
0 1
0 1
0 U1
0 V1
1 U 3
1 V3
53
09/02/10 14:54
On a:
N1
0
= [ k1 ] {q1} = [ k1 ] [ R1 ] {Q1 }
N1
0
E 2S 0
=
2 L 1
0 1
0 0
0 1
0 0
0
1
0 2 1
0 2 0
0
0
1 0
1 0
0 1
0 1
0 U1
0 V1
=
1 U 3
1 V3
(5.21)
2F
4
0
2F
4
0
do:
N1 =
2F
ES 2
N
= 3536 N xx1 = 1 = 25 MPa (5.22)
(U 3 + V3 ) =
4
L 2
2S
Llment 2: nuds 2 3, q = 0
1
1
u
N
2 2 4
ES
F
1 1 u = N = F N 2 = = 2500 N
4
L
3 2
4
N
xx2 = 2 = 25 MPa
S
(5.23)
Llment 3: nuds 3 4, q = 0
N 3 = k V3 =
3F
4
On notera quil est possible de trouver directement leffort dans le ressort grce la
relation liant rigidit et dplacement du nud 3.
9782100544639.indb 54
09/02/10 14:54
E, S, H
E, S, H
E, S, H
E,
S
, 2H
E,
2
S
, 2H
2
PX
4
Y
PY
Nud i
Nud j
Longueur
Section
1 (barre)
2 (barre)
3 (barre)
2H
4 (barre)
2H
5 (barre)
En adaptant la relation (5.1) aux spcificits de chacune des barres, les matrices en
repre local scriventpour:
Les lments 1, 2 et 5
9782100544639.indb 55
ES 1 1
H 1 1
(5.24)
55
09/02/10 14:54
Les lments 3 et 4
[k3 ] = [k4 ] =
S
2
2H
1 1 ES 1 1
1 1 =
2 H 1 1
(5.25)
PX
4
X
PY
Toujours en utilisant la forme gnrique (5.3), les rigidits en repre global scrivent pour:
Llment 1: nuds 1 2, q = 0
ES 0
[ K 1 ] = 1
H
0 1
0 0
0 1
0 0
0 (U1 )
0 (V1 )
0 (U 2 )
0 (V2 )
(5.26)
Llment 2: nuds 2 3, q = 0
[K 2 ] = [K1 ]
Llment 3: nuds 1 4, q = 45
1 1 1 1 (U1 )
ES 1 1 1 1 (V )
[ K 3 ] = 1 1 1 1 (U1 )
4H
4
1
1
1
1
(V
V
4)
(5.27)
56
9782100544639.indb 56
09/02/10 14:54
ES 1 1 1 1 (V3 )
[K 4 ] =
4 H 1 1 1 1 (U 4 )
V4 )
1 1 1 1 (V
(5.28)
0 0
ES 0 1
K
=
[ 5]
H 0 0
0 1
0 0 (U 2 )
0 1 (V2 )
0 0 (U 4 )
0 1 (V4 )
lment 5: nuds 2 4, q = 90
(5.29)
4
4
0 4+4 0
0
0
0 U 2 0
4 V2 0
0
0
4
0
0
0
ES 0
=
(5.30)
1 U 3 RH 3
4
1
1
0
0 4 +1 1
4H 0
1
1 V3 RV 3
0
0
0
1
1
0
1 1
1 1 1 + 1
1 + 1 U 4 PX
0
0
1
4 1 1 1 + 1 4 + 1 + 1 V4 PY
0
1
9782100544639.indb 57
5 1 4 0 0 0 1
1 1 0 0 0 0 1
4 0 8 0 4 0 0
4 0 0 0
ES 0 0 0
4 H 0 0 4 0 5 1 1
0 0 0 0 1 1 1
1 1 0 0 1 1 2
1 1 0 4 1 1 0
1 U1 RH 1
1 V1 RV 1
0 U 2 0
4 V2 0
=
1 U 3 RH 3
1 V3 RV 3
0 U 4 PX
6 V4 PY
(5.31)
57
09/02/10 14:54
8 0
ES 0 4
4 H 0 0
0 4
0 0 U 2 0
0 4 V2 0
=
2 0 U 4 PX
0 6 V4 PY
(5.32)
U 2 = 0
V = V
4
2
2P H
U 4 = X = 0.952 103 m
ES
2 PY H
= 1.905 103 m
V4 =
ES
(5.33)
58
9782100544639.indb 58
09/02/10 14:54
(5.34)
ES 2 PX H 2 PY H
PX + PY
= 150 kN
RH 1 = 4 H ES ES =
2
ES 2 PX H 2 PY H PX + PY
= 150 kN
RV 1 = 4 H ES + ES =
2
RH 3 = ES 2 PX H + 2 PY H = PX + PY = 50 kN
2
4H
ES
ES
RV 3 = ES 2 PX H + 2 PY H = PX + PY = 50 kN
4H
2
ES
ES
1 1
u1
v
2 1 1
1
=
u4 2 0 0
v4
0 0
1 1
u1
v
2 1 1
1
=
u4 2 0 0
v4
0 0
0 0 U1
0 0 V1
2
2 2 ( PX + PY ) H
U 4 V4 ) =
do d o u4 =
(
1 1 U 4
2
2
ES
1ce qui
1 permet
V
4 de trouver grce (4-37):
1 1 u1 N 3
1 1 u = N
4 3
(5.35)
2 ( PX + PY )
2 ( PX + PY ) H
= 212132 N
=
2
ES
ES
2H
N3 =
Dunod La photocopie non autorise est un dlit.
0 0 U1
0 0 V1
d o u = 2 (U V ) = 2 2 ( PX + PY
4
4
4
1 1 U 4
2
2
ES
1 1 V4
ES
2H
v
2 1 1 0 0 V3
2
2 2 ( PX
3
d o u4 =
U 4 V4 ) =
(
2
2
u4 2 0 0 1 1 U 4
v4
0 0 1 1 V4
u3
1 1 0 0 U 3
2 1 1 0 0 V3
2
2 2 ( PX + PY ) H
3
=
d o u4 =
U 4 V4 ) =
(
do
u4 2 0 0 1 1 U 4
2
2
ES
V
0
0
1
1
4
4
9782100544639.indb 59
59
09/02/10 14:54
1 1 u3 N 4
1 1 u = N
4 4
(5.36)
2 ( PX + PY ) H
2 ( PX + PY )
= 70711 N
=
2
ES
ES
2H
N4 =
ES
2H
60
9782100544639.indb 60
09/02/10 14:54
A
6.1 quation gnrale des poutres planes
y
G
Axe neutre
Les sections droites dune poutre en flexion restant droites aprs dformations
(principe de Navier-Bernoulli), ltude des dformations longitudinales sur un
tronon ABGH permet dtablir une relation linaire entre la dformation longitudinale et le rayon de courbure . En effet et de par sa dfinition (3.9), la dformation lordonne y est gale :
E F EF
xx =
(6.1)
EF
9782100544639.indb 61
61
09/02/10 14:54
(6.2)
xx y dS M z = 0
S
xx = E xx
1
M
E
M z = y 2 dS z =
E
S
EI z
= y
Iz
(6.3)
De plus et sous leffet des charges extrieures, le point situ sur laxe neutre labscisse x subira un dplacement vertical v( x ) et une rotation ( x ) . Par ailleurs et
d 1
daprs la figure 6.2, on sait que dx = d
mais galement que
=
dx
dv
=
ce qui permet de dduire:
dx
1 d d 2v
=
=
(6.4)
dx dx 2
En galisant (6.3) et (6.4), les relations moment-courbure et dformation-courbure deviennent:
1
M
d 2v
= z = 2
(6.5)
EI z dx
xx =
y
d 2v
= y 2
dx
(6.6)
dx
dT y
= q y ( x ) , on obtient finalement:
dx
d 2v
EI z 2 = M z ( x )
dx
d 3v
EI z 3 = T y ( x )
dx
d 4v
EI z 4 = q y ( x )
dx
(6.7)
(6.8)
(6.9)
62
9782100544639.indb 62
09/02/10 14:54
Les relations (6.7) et (6.9) sont les formes les plus connues de lquation gnrale
des poutres.
O
A
y
d
q y (x)
v(x)
d
ds
dx ds
9782100544639.indb 63
63
09/02/10 14:55
v(x)
y
i =
vi
dvi
dx
j =
vj
dv j
dx
u (x)
ui
i
E, S , I
uj
X
Figure 6.3 lment de poutre plane.
u( x )
En labsence de flexion, llment du fait de son raccourcissement ou de son allongement, a un comportement identique celui dune barre. u( x ) est donc encore
une fois une fonction linaire de la forme u( x ) = a0 + a1 x . La dformation associe
du
= a1 .
est bien videmment homogne sur la section et vaut comme en (4.2) xx =
dx
De plus et comme u ( 0 ) = ui en i et u ( L ) = u j en j, on obtient comme pour llment barre, la mme fonction de dplacement que celle tablie en (4.3):
x
x
u( x ) = (1 ) ui + u j
L
L
(6.10)
64
9782100544639.indb 64
09/02/10 14:55
v( x )
En flexion, chacun des nuds de cet lment possde deux degrs de libert, un
en translation dans la direction transversale v et un associ la rotation autour
de laxe z perpendiculaire au plan xy. Quatre conditions aux limites vi , i , v j , j
peuvent donc tre utilises pour dfinir sa fonction dapproximation ce qui explique
lexpression de v( x ) sous la forme dun polynme de degr 3.
v( x ) = b0 + b1 x + b2 x 2 + b3 x 3 = {1 x
x2
b0
b
1
x 3}
b2
b3
(6.11)
Cette approche revient en fait intgrer la forme (6.9) de lquation gnrale des
poutres avec une charge rpartie nulle. Lquation du moment flchissant sera alors
daprs (6.7) une fonction linaire.
partir des conditions aux limites de v( x ) (i .e . v(0) = vi , v (0) = i , v( L ) = v j , v ( L )
= j ) , le vecteur des dplacements nodaux scrit:
vi 1
0
i
v =
j 1
j 0
0 0
0 b0
b0
1 0
0 b1
= [ R ] 1
L L2 L3 b2
b2
2
b3
1 2 L 3L b3
(6.12)
v( x ) = {1 x
x2
Dunod La photocopie non autorise est un dlit.
v( x ) = {1 x
x2
vi
1 i
3
x } [R ]
v j
j
0
0
1
0
1
0
x 3}
2
3 / L 2 / L 3 / L2
3
1 / L2 2 / L3
2/L
0 vi
0 i
1 / L v j
1 / L2 j
(6.13)
Do finalement,
3x 2 2 x 3
v( x ) = 1 2 + 3
L
L
9782100544639.indb 65
2x 2 x 3
+ 2
L
L
3x 2 2 x 3
3
L2
L
vi
x 2 x 3 i
+ 2 (6.14)
L L v j
j
65
09/02/10 14:55
Cest dailleurs grce cette relation que certains logiciels arrivent tracer la dforme entre les nuds. Ils calculent les valeurs de v( x ) en certains points intermdiaires1 et reconstituent ensuite la dforme. Ceci permet de se limiter la
gomtrie de la structure et donc de rduire le nombre de nuds.
Champ de dformation
=
v( x )
x
0
1 L
3x 2 2 x 3
0
1 2 + 3
L
L
x
L
0
x
2x 2 x 3
+ 2
L
L
0
3x 2 2 x 3
3
L
L2
ui
v
i
0
i
x 2 x 3 u (6.15)
+ 2 j
L L v j
j
u( x )
= [ N ] {qe }
v( x )
xx =
du
d 2v
y 2
dx
dx
(6.16)
xx =
6 12 x
4 6x
2 3 y 2 y
L
L
L L
xx = [ B ] {qe }
1
L
ui
v
i
2 6 x i
6 12 x
2 + 3 y 2 y u
L
L L j
L
v j
j
(6.17)
1. Au quart, la moiti et aux trois quarts de la longueur de llment par exemple.
66
9782100544639.indb 66
09/02/10 14:55
2V
2V
e
1
{qe }T [ke ] {qe }
2
(6.18)
T
1
1
xx xx dVe = ([ B ] {qe }) E [ B ] {qe } dVe
2V
2V
e
[ke ]
We =
We =
1
{qe }T
2
T
E [ B ] [ B ] dVe {qe }
(6.19)
T
T
k
=
E
B
dV
=
E
[ e] [ ] [ ] e
dSe [ B ] [ B ] dx
0 Se
Ve
6 12 x y
L2 L3
4 6x
L
L L2 y
= E dSe
1
0 Se
L
6 12 x
2 + 3 y
L
L
2 6x
2 y
L L
9782100544639.indb 67
6 12 x
4 6x
2 3 y 2 y
L
L
L L
1
L
(6.20)
6 12 x
2 6x
2 + 3 y 2 y dx
L
L
L L
67
09/02/10 14:55
Rigidit
du nud i
ES
L
0
[ke ] = ES
Facteur de transmission
du nud i vers j
0
12 EI
L3
6 EI
L2
6 EI
L2
4 EI
L
12 EI
L3
6 EI
L2
6 EI
L2
2 EI
L
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
Facteur de transmission
du nud j vers i
avec:
dSe = S (section constante).
se
y dSe = 0
6 EI
L2
2 EI
6 EI
2
L
4 EI
L
0
(6.21)
Rigidit
du nud j
y 2 dSe = I
(inertie constante).
se
se
On notera que les termes en EI sont bien lis aux degrs de libert associs la
flexion de llment poutre.
6.2.2 Matrices de rigidit lmentaire en repre global
Llment de poutre plane possdant deux nuds trois degrs de libert, la matrice
de changement de base scriten vertu de (4.22):
n ( )
ui cos ( ) sin
v sin cos
( )
( )
i
i 0
0
{qe } = [ Re ] {Qe } u =
0
j 0
v j 0
0
0
j 0
0
0
0
0
0
0
1
0
0
0 cos ( ) sin ( )
0 sin ( ) cos ( )
0
0 U i
0 Vi
0 i
(6.22)
0 U j
0 V j
1 j
On notera cependant lajout du 1 associ aux rotations autour des axes z et Z qui
sont colinaires dans le cas plan. La matrice de rigidit lmentaire en repre global
est ensuite obtenue en appliquant (4.25), soit:
68
9782100544639.indb 68
09/02/10 14:55
12 EI 2
ES 2
L cos + L3 sin
ES 12 EI
L L3 cos sin
6 EI
2 sin
L
ES
cos2 12 EI sin 2
L
L3
ES 12 EI
+ 3 cos sin
L
L
6 EI
2 sin
ES 12 EI
3 cos sin
L
L
ES 2
12 EI
sin + 3 cos2
L
L
6 EI
cos
L2
ES 12 EI
+ 3 cos sin
L
L
12 EI
ES 2
sin 3 cos2
L
L
6 EI
cos
L2
9782100544639.indb 69
ES 2
12 EI
cos 3 sin 2
L
L
ES 12 EI
+ 3 cos sin
L
L
6 EI
sin
L2
ES 2
12 EI
cos + 3 sin 2
L
L
ES
EI
12
3 cos sin
L
L
6 EI
sin
L2
(6.23)
6 EI
sin
L2
6 EI
cos
L2
4 EI
L
6 EI
sin
L2
6 EI
2 cos
L
2 EI
L
[ K e ] = [ Re ] [ ke ][ Re ] =
ES 12 EI
+ 3 cos sin
L
L
ES
12 EI
sin 2 3 cos2
L
L
6 EI
2 cos
L
ES 12 EI
3 cos sin
L
L
ES 2
12 EI
sin + 3 cos2
L
L
EI
6
2 cos
L
6 EI
sin
L2
6 EI
cos
2
L
2 EI
6 EI
sin
L2
6 EI
2 cos
L
EI
4
69
09/02/10 14:55
qy
i
{ fe }
se
[ N ] { f s } dSe .
Se
T
{ f e } = { f es } = [ N ] { f s } dSe = 0 [ N ]T q dx
y
Se
3x 2 2 x
1 2 + 3
L
L
Fiy
2x 2 x 3
M
x
+ 2
L
iz
L
L
=
=
0
2
3
F
x
x
3
2
jy
L2
M jz
L3
x2 x3
+
L L3
qyL
q y L2
0
12
dx
q
qyL
y
2
q L2
12
(6.24)
Dune manire gnrale et dans le cas de forces orientes dans le sens y ngatif,
le vecteur charges dun lment de poutre plane soumis un systme de charges
transversales est gal linverse des ractions et moments dencastrement de la
poutre bi-encastre subissant le mme chargement. En dautres termes, ce vecteur
traduit les actions nodales quivalentes au chargement appliqu sur la poutre (cf.
tableau 6.1 des charges nodales quivalentes).
70
9782100544639.indb 70
09/02/10 14:55
9782100544639.indb 71
L/2
c
L
Type de chargement
P
2
qL
2
3qL
20
Pb2
( 2a + L )
L3
Fiy
( a + c )3 a3
q
L2
( a + c )4 a4
+
2L3
PL
8
(
(
(
qL2
30
)
)
)
6L2 ( a + c )2 a2
q
3
8L ( a + c ) a3
12L2
+3 ( a + c )4 a4
qL2
12
Pab2
L2
Miz
P
2
qL
2
7qL
20
Pa2
( L + 2b )
L3
Fjy
( a + c )3 a3
L2
q
4
( a + c ) a4
2L3
q
12L2
(
(
4L ( a + c )3 a3
4
3 ( a + c ) a4
qL2
20
qL2
12
Pa2b
L2
PL
8
Mjz
71
09/02/10 14:55
Une fois les dplacements en repre global connus, il suffit de calculer pour chacun
des lments, les dplacements associs en repre local puis de les injecter dans la
relation de rigidit lmentaire pour obtenir les efforts aux nuds correspondants.
On a donc:
Ni
T
i
Mi
N = [ ke ] {qe } { f e }
j
Tj
M j
ES
L
N
i 0
T
i
Mi 0
N =
j ES
Tj L
M j 0
12 EI
L3
6 EI
L2
6 EI
L2
4 EI
L
12 EI
L3
6 EI
L2
6 EI
L2
2 EI
L
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
(6.25)
6 EI
ui
L2
2 EI vi
L i f
{ e}
u j
0
vj
6 EI
2 j
L
4 EI
L
0
avec:
{qe } = [ Re ] {Qe } .
72
9782100544639.indb 72
09/02/10 14:55
y
i
z
Figure 6.5 lment de poutre 3D.
Lapproche sera alors identique celle utilise pour llment barre, do:
x ( x ) = a0 + a1 x
Dunod La photocopie non autorise est un dlit.
avecen x = 0 , x ( 0 ) = ix et x = L , x ( L ) = jx .
ce qui donne:
x
x
x ( x ) = (1 ) ix + jx
L
L
(6.26)
(6.27)
Les sections droites restant planes (torsion sans gauchissement), on peut tablir
pour la dformation de torsion que:
9782100544639.indb 73
=r
d x
1
= r
dx
L
1 ix
L jx
(6.28)
73
09/02/10 14:55
T
1
{ } { } dVe
2V
e
WT =
1
r ix
2 V
WT =
1
ix
2
1
L
1
jx G r
1
L
L
ix
[ kT ]
jx
1 ix
dVe
L jx
(6.29)
jx }
dx
r
dx
0 (centre de torsion)
Figure 6.6 Relation entre angle et dformation de torsion.
[kT ] =
GJ
L
1
L2
1
dx = GJ
1
L
0
L2
L
1
L2 dx
1
L2
(6.30)
1 1
1 1
9782100544639.indb 74
09/02/10 14:55
Bien videmment, cette matrice de rigidit est trs similaire celle de llment
barre (5.1).
6.3.2 Prise en compte de la dformation deffort tranchant
Considrant une poutre console deux nuds, la relation de rigidit dans le plan
xy scrit:
0
12 EI z
L3
6 EI z
L2
6 EI z
L2
4 EI z
L
12 EI z
L3
6 EI z
L2
6 EI z
L2
2 EI z
L
ES
L
0
12 EI z
L3
6 EI
2z
L
0
ES
L
0
12 EI z
L3
6 EI
2z
L
0
0
6 EI z
ui N i
L2
2 EI z vi Ti
L i = Mi
N
u
0 j j
v j T j
6 EI
2 z j M j
L
4 EI z
ES
L
ES
L
(6.31)
xy
1
L3
EI z 12 6 L v j T j v j 3 EI z
= =
L3 6 L 4 L2 j M j j L2
2 EI
z
L2
2 EI z T j
L M j
EI z
(6.32)
xy = G xy
(6.33)
T y = xy S1, y
(6.34)
9782100544639.indb 75
75
09/02/10 14:55
xy =
Ty
GS1, y
(6.35)
Ty L
GS1, y
(6.36)
v j 3 EI z GS1, y
=
L2
j
2 EI z
12 EI z
3
Tj L 1+ y
=
M j 6 EI z
L2 1 + y
avec y =
)
)
L2
2 EI z T j
L M j
EI z
6 EI z
2
L 1+ y v
j
EI z 4 + y j
L 1 + y
(
(
(
(6.37)
12 EI z
.
GS1, y L2
76
9782100544639.indb 76
09/02/10 14:55
ES
L
ES
L
9782100544639.indb 77
(1 + y )
6 EI z
(1 + y )
12 EI z
(1 + y )
6 EI z
L3
L2
L2
L3 1 + y
12 EI z
(1 + z )
6 EI y
L2
(1 + z )
6 EI y
L3 (1 + z )
12 EI y
L2
L3 (1 + z )
GJ
L
GJ
L
12 EI y
(6.38)
( 2 y ) EI z
L (1 + y )
(1 + y )
6 EI z
( 2 z ) EI y
L (1 + z )
L2
ES
L
L2 (1 + z )
6 EI y
( 4 + y ) EI z
L (1 + y )
ES
L
( 4 + z ) EI y
L (1 + z )
L2 (1 + z )
L2 1 + y
6 EI z
6 EI y
12 EI z
(1 + y )
6 EI z
(1 + y )
L2
(1 + y )
6 EI z
12 EI z
L2
L3 1 + y
L3
L2
(1 + z )
6 EI y
L3 (1 + z )
12 EI y
(1 + z )
6 EI y
L3 (1 + z )
L2
12 EI y
GJ
L
GJ
L
L2 (1 + z )
( 4 + z ) EI y
L (1 + z )
L2 (1 + z )
6 EI y
( 2 z ) EI y
L (1 + z )
6 EI y
6 EI z
L2 1 + y
2 y EI z
L 1+y
6 EI z
L2 1 + y
4 + y EI z
L 1 + y
0
77
09/02/10 14:55
z =
12 EI y
GS1,z L2
12 EI y
12 EI z
et
=
z
GS1,z L2
GS1, y L2
deffort tranchant. Ces deux paramtres sont lis aux sections deffort tranchant
S1, y et S1,z respectivement associes aux axes y et z.
titre dexemple, la section deffort tranchant dune section rectangulaire est gale
5/6 de sa section axiale, celle dune section circulaire 9/10.
Si ces sections sont nulles, la formulation de llment revient celle de llment
de poutre classiquedit lanc, y et z tant pris gaux zro dans ce cas. On sera
dailleurs dans la mme configuration lorsque L est grand.
Y
1
ky
L/2
Nud i
Nud j
Section
Inertie
1 (poutre)
2 (poutre)
3 (ressort)
Application numrique: E = 2.1 1011 N/m2, L = 10 m, S = 0.00285 m2, I = 0.00001943 m4 (IPE200), P = 1000 N.
9782100544639.indb 78
09/02/10 14:55
PL
8
PL
8
i
L
P
2
P
2
6.4.1 ky = 0
Dans ce cas, il sagit dune poutre sur deux appuis simples ( v1 = v2 = 0 ) subissant
une charge centre P. Le porte faux ntant pas charg, il est possible de considrer uniquement llment 1 en posant que V3 = L 2 .
Le systme rsoudre se rduit donc :
6 EI
12 EI
6 EI
P
3
2
2
2
L
L
L
v1
PL RV1
4 EI
6 EI
2 EI
2
L
L
L 1 = 8 + 0
6 EI v2 P RV2
6 EI 12 EI
(6.39)
2
L3
L2 2 2 0
L
PL
2 EI
6 EI
4 EI
2
L
L
L
8
PL2
1 =
= 2 = 0.00153rad
16 EI
12 EI
L3
6 EI
L2
12 EI
3
L
6 EI
L2
9782100544639.indb 79
P 6 EI
RV1 = 2 + L2 ( 1 + 2 ) =
RV = P 6 EI ( 1 + 2 ) =
2 2 L2
P
2
P
2
(6.40)
79
09/02/10 14:55
N1
T
1
M1
et les efforts dans llment 1 = [ k1 ] {q1} { f1} qui est quivalent :
N2
T2
M 2
N1
T
1
M1
=
N2
T2
M 2
ES
0
L
12 EI
0
L3
6 EI
0
L2
ES
0
L
12 EI
0
3
6 EI
0
L2
0
6 EI
L2
4 EI
L
0
6 EI
L2
2 EI
L
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
P
6 EI
u
1
L2 2
2 EI v1 PL
L 1 8
u2 0
0
P
v2
6 EI 2
2 2
L
PL
4 EI
8
L
0
(6.41)
Ce qui permet de montrer que les moments sur appuis sont bien nuls.
2 EI
PL 4 EI
PL 2 EI
PL PL
M1 = 8 + L 1 + L 2 = 8 + L 1 = 8 8 = 0
(6.42)
M 2 = PL + 2 EI 1 + 4 EI 2 = PL + 2 EI 2 = PL + PL = 0
8
8
8
L
L
8
L
Le moment flchissant au milieu de la premire trave peut tre dtermin partir
de lquation gnrale des poutres (6.7) et de la double drivation de (6.14) ce qui
donne:
d 2v
PL
L
M = EI 2 =
dx
8
2
(6.43)
80
9782100544639.indb 80
09/02/10 14:55
PL PL
PL
L PL
MA = M
=
=
= 2500 Nm
8
8
4
2 8
(6.44)
A
Figure 6.10 Exemple 3: diagramme (Effel) du moment flchissant ky = 0.
L
partir de la
2
relation (6.14), do :
v1
PL3
L 1 L 1
L 1 L
v( x = ) =
= ( 1 2 ) =
(6.45)
2 2 8 2
8 v2 8
64 EI
2
auquel il faut ajouter le dplacement au mme endroit de la poutre bi-encastre
PL3
correspondante
soit finalement au point A:
192 EI
vA =
PL3
PL3
PL3
=
= 5.106 mm
64 EI 192 EI
48EI
De plus et comme v3 = L 2 , on a v3 =
(6.46)
PL3
= 15.318 mm .
16 EI
6.4.2 ky
La rigidit tant infinie lappui 3, la poutre est dsormais sur trois appuis simples
(v1 = v2 = v3 = 0 ) . Il sagit donc dune poutre continue deux traves. Lassemblage des lments 1 et 2 au nud 2 revient additionner les rigidits des degrs
12 EI 12 EI
de libert correspondants ce qui permet de trouver des termes en 3 + 3 et
L
L
9782100544639.indb 81
81
09/02/10 14:55
4 EI 4 EI
6 EI 6 EI
+
respectivement associs v2 et 2 alors que ceux en 2 2
L
L
L
L
sannulent. Une fois assembl, le systme a pour expression:
6 EI
12 EI
12 EI
3
L3
2
L
L
6
EI
4
EI
6
EI
L2
L
L2
12 EI
6 EI 12 EI 12 EI
+ 3
2
3
L3
L
L
L
2 EI
6 EI 6 EI
6 EI
2
L2
L
L2
L
12
EI
0
0
3
6 EI
0
0
L2
P
2
RV
PL 1
8 0
P RV
= + 2
2
0
PL
RV
8 3
0 0
6 EI
L2
2 EI
L
6 EI 6 EI
2
L2
L
4 EI 4 EI
+
L
L
6 EI
2
L
2 EI
L
0
0
12 EI
L3
6 EI
2
L
12 EI
L3
6 EI
2
L
0 v1
6 EI 1
L2 v2
2 EI 2
L v3
6 EI
2
3
L
4 EI
(6.47)
L
0
Assemblage au nud 2
2 EI
L
2 EI
L
8EI
L
2 EI
L
d
0
8 1
64 EI
2 EI PL
PL2
= 0.00077 rad (6.48)
2 =
2 =
32 EI
L 8
4 EI 3 0
PL2
= 0.00038 rad
3
L
64 EI
82
9782100544639.indb 82
09/02/10 14:55
RV1 =
RV2
RV3
(6.49)
N1
T
1
M1
Les efforts dans llment 1 sont obtenus en posant = [ k1 ] {q1} { f1} soit,
N2
T2
M 2
ES
0
0
L
12 EI
6 EI
N
1 0
3
L
L2
T
1
6 EI
4 EI
M1 0
2
L
L
=
N
ES
2
0
0
T2 L
12 EI
6 EI
3
2
M 2 0
L
L
6 EI
2 EI
0
L2
L
N1 13 P 0
T 32 406.25
1 0
M1
0
= 0 =
0
N2
T2 19 P 593.75
32
M 2 3 PL 937.50
32
9782100544639.indb 83
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
P
6 EI
u1
L2 2
2 EI v1 PL
L 1 8
0
u
0 2
P
v2
6 EI 2
2
2
L
PL
4 EI
8
(6.50)
83
09/02/10 14:55
N2
T
2
M 2
et ceux de llment 2 partir de = [ k2 ] {q2 } { f 2 } soit,
N3
T3
M 3
ES
0
0
L
12 EI
6 EI
N2 0
3
L
L2
T
6 EI
4 EI
2
M 2 0
2
L
L
=
N
ES
3
0
0
T3 L
12 EI
6 EI
3
2
M 3 0
L
L
6 EI
2 EI
0
L2
L
3P
0
32 93.75
3 PL
937.50
= 32 =
0 0
93.75
3P
32 0
0
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
6 EI
u2
L2
2 EI v2
L 2
u
0 3
v3
6 EI
2 3
L
4 EI
(6.51)
L
0
84
9782100544639.indb 84
09/02/10 14:55
Pour la trave 1
v( x =
9782100544639.indb 85
L 1
)=
2 2
L
8
1
2
v1
23 PL3
L 1 PL3
=
= 3.67 mm
8 v2 192 EI
1536 EI
2
(6.52)
85
09/02/10 14:55
Pour la trave 2
v( x =
L 1
)=
2 2
1
2
L
8
v2
L 2 3 PL3
= 1.436 mm
=
8 v3 512 EI
3
(6.53)
6 EI
12 EI
L3
L2
4 EI
6 EI
2
L
L
12 EI
6 EI
2
3
L
L
2 EI
6 EI
L2
L
0
0
0
0
P
2
RV
PL 1
8 0
P RV
= + 2
2
0
PL
RV
8 3
0 0
12 EI
L3
6 EI
2
L
24 EI
L3
12 EI
L3
6 EI
L2
6 EI
L2
2 EI
L
0
8EI
L
6 EI
2
L
2 EI
L
0
0
12 EI
L3
6 EI
2
L
12 EI
+ ky
L3
6 EI
2
L
0 v1
6 EI 1
L2 v2
2 EI 2
L v3
6 EI
2
3
L
4 EI
(6.54)
L
0
86
9782100544639.indb 86
09/02/10 14:55
2 EI
L
8EI
L
6 EI
2
L
2 EI
L
6 EI
L2
12 EI
+ ky
L3
6 EI
2
L
PL
1 8
2 EI
L 2 = PL
8
6 EI v
2 3 0
L 3
0
4 EI
(
(
)
)
3 PL2 2 EI + k y L3
1 =
= 0.00129 rad
32 EI 3 EI + 2k y L3
PL2 3 EI + k y L3
= 0.00106 rad
2 =
16 EI 3 EI + 2k y L3
3 PL3
v =
3 16 3 EI + 2k L3 = 0.005816 m
y
PL2 6 EI + k y L3
3 =
= 0.00034 rad
3
EI
EI
+
k
L
32
3
2
y
(
(
4 EI
L
2 EI
L
(6.55)
(
(
)
)
v( x =
L 1
)=
2 2
L
8
1
2
v1
L 1 PL3
= 4.215mm
8 v2 192 EI
2
(6.56)
9782100544639.indb 87
87
09/02/10 14:55
RV2
3
6 EI
P P 24 EI + 13k y L
+
=
= 441.84 N
(
)
1
2
2
L2
16 3 EI + 2k y L3
(
)
3
P P (12 EI + 11k y L )
6 EI
= 2 ( 1 + 3 ) + =
= 616.31 N
L
2
8 ( 3 EI + 2k y L3 )
RV3 = k y v3 =
3 PL3 k y
16 3EI + 2k y L3
(6.57)
= 58.16 N
Pour llment 1
0
3
P 24 EI + 13k y L
N1 16 3EI + 2k L3
y
T
1
0
M1
0
=
N
2
P 24 EI + 19k L3
y
T2
16 3 EI + 2k y L3
M 2
3 PL4 k y
3
16 3 EI + 2k y L
(6.58)
441.84
0
=
0
558.16
581.56
3 PL3 k y
N 2 16 3EI + 2k L3
y
T
2
4
3PL k y
M 2
= 16 3 EI + 2k y L3
N3
0
T3
3 PL3 k y
M 3
16 3 EI + 2k y L3
0
58.16
581.56
=
0
58.16
(6.59)
Pour llment 2
9782100544639.indb 88
09/02/10 14:55
y
q
x
2
E , S , I ( IPE 270)
E , S , I ( IPE 270)
Y
X
E , S , I ( IPE 270)
x
y
9782100544639.indb 89
89
09/02/10 14:55
Nud i
Nud j
Longueur (m)
Section (m2)
Inertie (m4)
1 (poutre)
2 (poutre)
3 (poutre)
Application numrique: L = 10 m, q = 10 kN/m, E = 2.1 1011 N/m2, S = 0.00459 m2, I = 0.0000579 m4.
0
0
0
0
L
L
EI
EI
12
6
12 EI
6 EI
0
0
3
L3
L2
L
L2
6 EI z
2 EI
6 EI
4 EI
0
L
L2
L (6.60)
L2
k
=
k
=
k
=
[ 1] [ 2 ] [ 3 ]
ES
ES
0
0
0
0
L
L
EI
EI
12
6
12 EI
6 EI
0
3
2
0
L
L
L3
L2
2 EI
6 EI
4 EI
6 EI
0
L2
L
L2
L
6.5.2 Matrices de changement de repres
[ R1 ] = [ R2 ] = [ I ]
(6.61)
0
[ R3 ] = 0
sin
cos
0
0
0
0
0
0
0
0
1
0
0 cos
0 sin
0
0
0
0
0
sin
cos
0
0 0
0 1
0 0
=
0 0
0 0
1 0
1
0
0
0
0
0
0
0
1
0
0
0
0
0
0
0
0
1
0
0
0
(6.62)
0
0 1 0 0
0 0 0 1
90
9782100544639.indb 90
09/02/10 14:55
12 EI
L3
6 EI
L2
6 EI
L2
4 EI
L
12 EI
L3
6 EI
L2
6 EI
L2
2 EI
L
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
6 EI
L2
2 EI
6 EI
2
L
4 EI
L
0
A
(6.63)
ES
L
0
K
=
k
=
[ 1 ] [ 1 ] ES
Rigidit du nud 2
ES
L
0
[ K 2 ] = [k2 ] = ES
9782100544639.indb 91
12 EI
L3
6 EI
L2
6 EI
L2
4 EI
L
12 EI
L3
6 EI
L2
6 EI
L2
2 EI
L
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
6 EI
L2
2 EI
L (6.64)
6 EI
2
L
4 EI
L
0
91
09/02/10 14:55
Pour llment 3
12 EI
L3
6 EI
2
[ K 3 ] = [ R3 ]T [k3 ][ R3 ] = 12LEI
L3
6 EI
2
L
0
ES
L
0
0
ES
L
0
6 EI
L2
12 EI
L3
4 EI
L
6 EI
L2
6 EI
L2
12 EI
L3
2 EI
L
6 EI
L2
ES
L
0
0
ES
L
0
6 EI
L2
2 EI
L (6.65)
6 EI
L2
4 EI
Rigidit du nud 2
6.5.4 Vecteurs charges
La charge rpartie q applique chacun des deux lments est quivalente aux
charges nodales suivantes (cf. tableau 6.1):
0
qL
2
qL2
12
{ f1} = { f 2 } = 0
qL
2
qL2
12
(6.66)
Les repres locaux des lments 1 et 2 tant confondus avec le repre global, on a
bien sr:
92
9782100544639.indb 92
09/02/10 14:55
0
qL
2
qL2
12
{F1} = { f1} =
0
qL
2
qL2
12
0
qL
2
qL2
12
{F2 } = { f 2 } =
(6.67)
0
qL
2
qL2
Charges nodales
12
du nud 2
De par les conditions dappui, seul le nud 2 est libre. Il suffit donc dadditionner
les parties de rigidit et de vecteurs charges encadres aux chapitres 6.5.3 et 6.5.4:
ES ES 12 EI
L + L + L3
6 EI
L2
qL qL
=
2
2
2
2
qL
qL
12 12
0
12 EI 12 EI ES
+ 3 +
L3
L
L
6 EI 6 EI
2
L2
L
6 EI
2
L
U 2
6 EI 6 EI
V2
2
2
L
L
4 EI 4 EI 4 EI B2
+
+
L
L
L
(6.68)
2 ES 12 EI
L + L3
6 EI
L2
0
24 EI ES
+
L3
L
0
6 EI
L2 U 2 0
0 V2 = qL
12 EI B2 0
(6.69)
9782100544639.indb 93
93
09/02/10 14:55
qL
= 1.034 mm
(6.70)
24 EI ES
+
L3
L
La structure tant symtrique par rapport un axe passant par X = L , il aurait t
possible de rduire le problme aux seuls lments 1 et 3. Dans cette hypothse,
il aurait fallu en plus des conditions dappui, bloquer le nud 2 horizontalement
mais galement en rotation (U 2 = B2 = 0 ) .
Cependant et pour tre totalement quivalente lapproche globale, ces calculs
devraient prendre en compte la moiti de linertie I pour llment 3, celui-ci
tant dans laxe de symtrie. Dune manire gnrale, les caractristiques des entits
(charges, appuis lastiques, etc.) situes sur laxe de symtrie doivent tre divises
par deux.
Enfin, un moyen mnmotechnique simple permet de se souvenir des degrs de
libert bloquer: on bloque la translation suivant laxe de symtrie passant par
une valeur de X,Y ou Z puis les rotations suivant les deux autres axes ce qui se
rsume par le tableau suivant:
V2 =
Translation
bloque
Rotation bloque
autour de
Rotation bloque
autour de
94
9782100544639.indb 94
09/02/10 14:55
Les dplacements aux milieux des deux traverses peuvent galement tre calculs
partir de (6.45) et (6.46) ce qui donne:
v( x =
L 1
)=
2 2
1
2
L
8
v1
qL4
L 1
= 21.935mm
8 v2 384 EI
2
(6.71)
N1
T
1
M1
= [ k1 ] {q1} { f1}
N2
T2
M 2
ES
L
N1 0
T
1
M1 0
=
N 2 ES
T2 L
M 2 0
12 EI
L3
6 EI
L2
6 EI
L2
4 EI
L
12 EI
L3
6 EI
L2
6 EI
L2
2 EI
L
ES
L
0
0
ES
L
0
0
0
12 EI
L3
6 EI
2
L
0
12 EI
L3
6 EI
2
L
qL
6 EI
u1
L2 2
2 EI v1 qL2
L 1 12 (6.72)
0
u
0 2
v2 qL
6 EI
2 2 2
L
qL2
4 EI
12
L
0
N1 0
u1
U1
T 50151
v
V
1
1
1
M1 84088
1
1
avec = [ I ] do =
.
N2 0
u2
U 2
T2 49849
v2
V2
M 2 82579
2
2
9782100544639.indb 95
95
09/02/10 14:55
96
9782100544639.indb 96
09/02/10 14:55
M 4
k
q
=
[ 3 ] { 3} { f3}
N2
T2
M 2
ES
0
0
L
12 EI
6 EI
N4 0
3
L
L2
T
6 EI
4 EI
4
M 4 0
2
L
L
=
N
ES
2
0
0
T2 L
12 EI
6 EI
3
2
M 2 0
L
L
6 EI
2 EI
0
L2
L
99698
N
N
=
T
0
4
M 4 0
N 2 = + N 99698
T2
0
M 2 0
u4
U 4 0
v
V 1
4
4
4
B4 0
avec = [ R3 ] =
u2
U 2 0
v2
V 2 0
2
B2 0
9782100544639.indb 97
1
0
0
0
0
0
ES
L
0
12 EI
L3
6 EI
2
L
0
ES
L
0
12 EI
L3
6 EI
2
L
0
0
0 0
0 0
1 0
0 0
0 1
0 0
0
0
0
1
0
0
6 EI
u4
L2
2 EI v4
L 4
u
0 2
v2
6 EI
2 2
L
4 EI
L
(6.73)
0
0 U 4
0 V4
u = V
2
0 B4 2
v
U
=
2
0 U 2 2
= B2
0 V 2 2
1 B2
97
09/02/10 14:56
Les contraintes peuvent ensuite tre dduites des efforts grce aux relations issues
de la rsistance des matriaux, savoir pour:
Les contraintes normales
xx =
N My
M
z z y
S
Iy
Iz
(6.74)
y =
Ty
Sy
et z =
Tz
Sz
(6.75)
98
9782100544639.indb 98
09/02/10 14:56
9782100544639.indb 99
99
09/02/10 14:56
9782100544639.indb 100
09/02/10 14:56
7 lments isoparamtriques
A
7.1 Problmatique du maillage
Les modlisations filaires abordes aux chapitres 5 et 6 ont montr que le strict respect des singularits gomtriques (positions des appuis, des charges, des attaches,
etc.) suffisait obtenir en statique des rsultats identiques ceux de la thorie des
poutres. Ceci sexplique par la conformit des fonctions de forme choisies avec la dite
thorie. Hormis quelques cas comme les charges (forces ou moments) ponctuelles,
la discrtisation intermdiaire des lments savre rarement ncessaire en statique.
Il nen est pas de mme pour les modlisations surfaciques ou volumiques qui
ncessitent une discrtisation beaucoup plus fine. En effet et afin de coller au mieux
la gomtrie, cette reconstitution requiert surtout au niveau des singularits, des
lments de petites tailles.
Cette opration est prise en charge par un outil appel mailleur dont le rle va tre
dtablir automatiquement, dans le respect de la gomtrie tudie, la forme des
lments, leur connectivit et les coordonnes des nuds. Ceci tant et comme on
peut le voir la figure 7.1, la taille et la forme des lments peuvent varier significativement. Enfin, deux grands types de maillage sont possibles partir de bases
quadrangulaires (carrs, hexadres, etc.) ou triangulaires (triangles, ttradres, etc.).
9782100544639.indb 101
101
09/02/10 14:56
7 lments isoparamtriques
Les lments pouvant tre de tailles et de formes quelconques, il nest pas envisageable, pour tablir leurs caractristiques lmentaires, de reprendre stricto sensu
la mthodologie utilise pour les poutres et barres. Les matrices de rigidit et vecteurs charges variant systmatiquement, une mthodologie permettant de calculer
ces caractristiques quelles que soient les gomtrie et configuration des lments,
savrera beaucoup plus rentable au niveau calcul.
Celle-ci consistera, pour chaque type dlment, dfinir un lment de rfrence de gomtrie conventionnelle (appel galement lment parent) de telle
manire obtenir la gomtrie de nimporte quel lment rel de forme semblable
partir dune transformation gomtrique biunivoque (i.e. bijective).
En dautres termes, pour chacun des points P ( , ) de llment parent dfini
dans un repre unitaire correspondra, via la transformation e , un point et un seul
de llment rel P ( x , y ) .
Y
2
(1,1)
P' ( , )
(-1,-1)
(-1,1)
k
P ( x, y )
j
l
(1,-1)
X
4
lment rel
lment parent
x ( , )
e : ,
y ( , )
(7.1)
(
(
)
)
x , , xi , x j , x k , xl
e : ,
y , , yi , y j , yk , yl
(7.2)
102
9782100544639.indb 102
09/02/10 14:56
7 lments isoparamtriques
Y
3
2
(-1,1)
(1,1)
(-1,-1)
(1,-1)
4
lments rels
lment parent
Du fait du caractre biunivoque de la transformation, chacun des points de llment de rfrence devra concider avec un point de llment rel et un seul. Il en
rsulte une correspondance entre nuds gomtriques des lments rel et de rfrence. En dautres termes, il nest pas par exemple possible dassocier un lment
parent triangulaire un lment rel carr.
Chaque point de coordonnes x, y, z de llment rel peut donc tre repr partir
des coordonnes de ses nuds et de variables N i ( x , y , z ) , similaires aux fonctions
de forme vues prcdemment, soit:
n1
x = N i ( x , y , z ) xi
i =1
n1
y = N i ( x , y , z ) yi
(7.3)
i =1
n1
z = N i ( x , y , z ) zi
i =1
avec:
ij = 1 si i = j
N i ( x , y , z ) = ij
.
ij = 0 si i j
n1: nombre de nuds gomtriques de llment rel.
Par ailleurs, la notion de fonctions de forme aborde au chapitre 4.1 nous a permis
dtablir pour les dplacements que:
9782100544639.indb 103
103
09/02/10 14:56
7 lments isoparamtriques
n2
u ( x , y , z ) = N i ( x , y , z ) ui
i =1
n2
v ( x , y , z ) = N i ( x , y , z ) vi
(7.4)
i =1
n2
w ( x , y , z ) = N i ( x , y , z ) wi
i =1
du
= Cte u( ) = a0 + a1
d
104
9782100544639.indb 104
09/02/10 14:56
7 lments isoparamtriques
x = 1 u ( 1) = u1
x = 1 u (1) = u2
do:
i =2
1
1
u( ) = (1 ) u1 + (1 + ) u2 = N i ui
2
2
i =1
N1
1
-1
N2
u1
(7.5)
0
u ( )
u2
Ces nouvelles fonctions de forme sont tout fait quivalentes celles trouves en
x
(4.3), le changement de variable = 2 1 permettant de passer des unes aux
L
autres. Ceci nous permet galement dvoquer les proprits trs particulires des
fonctions de forme. Celles-ci sont en effet gales lunit aux nuds auxquels elles
sont rattaches, zro ailleurs1. De plus, la somme des fonctions de forme est dans
ce cas gale 12 ce qui permet de dmontrer la condition de mouvement de corps
rigide (MCR), soit:
Translation densemble
u0 u1 = u2 = u0
do
u( ) = N1 u0 + N 2 u0 = u0
(7.6)
1. En i Ni = 1, en i j Ni = 0.
2.
Ni = 1 .
i =1
9782100544639.indb 105
105
09/02/10 14:56
7 lments isoparamtriques
7.2.1 Famille C0
u( , ) = ai j k
(7.7)
i =1
2
2
3
2
22
32
Cubique
Parabolique
3
3
23
33
Linaire
2
2
3
Cubique
Parabolique
3
3
106
9782100544639.indb 106
09/02/10 14:56
7 lments isoparamtriques
Linaire
2
3
Parabolique
Cubique
u( , ) = a0 + a1 + a2 + a3 = N i ui
(7.8)
i =1
(-1,1)
(-1,-1)
(1,1)
(1,-1)
4
Figure 7.8 lment Q4.
Il sagit en fait de lquation dun plan. Le champ de dplacement vertical peut tre
4
Pour chacun des segments de llment, il est possible de dfinir lquation dune
droite fonction des valeurs nodales do son caractre bilinaire. Ses fonctions de
forme ont donc pour expressions:
N1 =
1
1
(1 + ) (1 + ) N 2 = (1 ) (1 + )
4
4
N3 =
1
1
(1 ) (1 ) N 4 = (1 + ) (1 )
4
4
(7.9)
9782100544639.indb 107
107
09/02/10 14:56
7 lments isoparamtriques
On notera que ces expressions peuvent tre galement obtenues en effectuant les
produits croiss des fonctions de forme de llment barre (7.5). De manire plus
synthtique, celles-ci peuvent sexprimer en fonction de coordonnes i ,i des
nuds de llmenten posant:
1
N i = (1 + i ) (1 + i )
(7.10)
4
Quadrangle quadratique complet Q9
Cet lment est quadratique, chacun de ces segments comportant trois nuds.
Ceux-ci permettent une approximation parabolique sur chacun de ses cts.
Nanmoins, la fonction dapproximation naurait pu tre totalementparabolique
sans le nud milieu, le terme 2 2 tant obtenu grce au nud n9.
Cest pourquoi, il est dit complet ce qui signifie quil couvre tous les termes du
degr de la fonction dapproximation vise (parabolique).
v
4
On a donc:
9
u( , ) = N i ui
i =1
= a0 + a1 + a2 + a3 + a4 2 + a5 2
+ a6 2 + a7 2 + a8 2 2
(7.11)
= {a} { X }
avec:
{a}T = {a0 a1 a2 a3 a4 a5 a6 a7 a8 }
{ X }T = {1 2 2 2 2 2 2 }
108
9782100544639.indb 108
09/02/10 14:56
7 lments isoparamtriques
Pour tablir les fonctions de forme de cet lment, il suffit de reprendre la mthodologie dcrite au chapitre 6.2.1 en posant que:
u (1,1)
u ( 0,1)
u ( 1,1)
u ( 1, 0 )
u ( 1, 1) = [ R ] {a}
u ( 0, 1)
u (1, 1)
u 1, 0
( )
u ( 0, 0 )
(7.12)
avec
1
1
1
1
[ R ] = 1
1
1
1
0
1
1
1
0
1
1
0
1
1
1
0
1
1
1
0
0
1
0
1
0
1
0
1
0
0
1
0
1
1
1
0
1
1
0
1
1
1
0
1
1
1
0
0
T
1
0
1
0
1
0
1
0
0
1
0
1
0
1
0
1
0
0
0
1
0
1
0
1
N1 =
N3
N5
N7
N9
9782100544639.indb 109
1
(1 2 ) (1 + )
2
1
N 4 = (1 ) (1 2 )
2
1
N 6 = (1 2 ) (1 )
2
1
N 8 = (1 + ) (1 2 )
2
N2 =
(7.13)
109
09/02/10 14:56
7 lments isoparamtriques
v
4
u( , ) = N i ui
i =1
u( , ) = a0 + a1 + a2 + a3 + a4 2
+ a5 2 + a6 2 + a7 2
u( , ) =
(7.14)
{a}T { X }
avec:
{a}T = {a0 a1 a2 a3 a4 a5 a6 a7 }
{ X }T = {1 2 2 2 2 }
Comme,
u (1,1)
u ( 0,1)
u ( 1,1)
u ( 1, 0 )
= [ R ] {a}
,
1
1
u
(
)
u ( 0, 1)
u (1, 1)
u (1, 0 )
(7.15)
110
9782100544639.indb 110
09/02/10 14:56
7 lments isoparamtriques
avec
1
0
1
1
1
0
1
1
1
1
1
0
1
1
1
0
1
0
1
0
1
0
1
0
1
0
1
1
1
0
1
1
1
1
1
0
1
1
1
0
1
0
1
0
1
0
1
0
1
0
0
1
0
1
1
1
1
[ R ] = 1
1
1
1
(1 2 ) (1 + )
2
1
N 4 = (1 ) (1 2 )
2
(7.16)
1
2
N 6 = (1 ) (1 )
2
1
N 8 = (1 + ) (1 2 )
2
N1 =
N2 =
Hexadre trilinaire H8
8
u( , , ) = N i ui
i =1
w
3
u
4
1
9782100544639.indb 111
111
09/02/10 14:56
7 lments isoparamtriques
1
(1 + i ) (1 + i ) (1 + i )
8
(7.17)
u( , ) = a0 + a1 + a2 = N i ui
i =1
N1( , ) = 1
N 2 ( , ) =
N 3 ( , ) =
(7.18)
(0,1)
(1,0)
u( , ) = a0 + a1 + a2 + a3 2 + a4 2 + a5
6
u( , ) = N i ui
i =1
N1( , ) = (1 ) (1 2 (1 ) )
N 2 ( , ) = 4 (1 )
N 3 ( , ) = (1 2 )
N 4 ( , ) = 4
(7.19)
N 5 ( , ) = (1 2 )
N 6 ( , ) = 4 (1 )
112
9782100544639.indb 112
09/02/10 14:56
7 lments isoparamtriques
5
4
u( , , ) = N i ui
i =1
(7.20)
w
u
7.2.2 Famille C1
9782100544639.indb 113
113
09/02/10 14:56
7 lments isoparamtriques
dui dui
dw
,
,..., i . De ce fait, les fonctions de forme doivent tre au minimum de
dx dy
dz
degr 2 pour satisfaire au critre de la dformation constante ; la condition de
mouvement de corps rigide tant quant elle assure par les deux premiers termes
de cette mme fonction.
Ces lments ont une continuit dite dordre 1, la continuit des dplacements
tant assure la fois sur llment et ses frontires si celui-ci est conforme (cf.
8.3.3). Llment poutre tel que celui dcrit au chapitre 6.2.1 fait partie de cette
famille.
Plaque Q4 (semi C1)
wi
wi
Sachant que cet lment trois degrs de libert par nud {wi } =
, la
wi
fonction dapproximation scrit:
w( , ) = a0 + a1 + a2 + a3 + a4 2 + a5 2 + a6 2
+ a7 2 + a8 3 + a9 3 + a10 3 + a11 3
T
= {a} { X }
(7.21)
avec:
9782100544639.indb 114
09/02/10 14:56
7 lments isoparamtriques
w
=
(7.22)
a1 + a3 + 2 a4 + 2 a6 + a7 2 + 3a8 2 + 3a10 2 + a11 3
w ( , ) =
do:
w(1,1)
w (1,1)
w (1,1)
w( 1,1)
w ( 1,1)
w ( 1,1)
= [ R ] {a}
w( 1, 1)
w ( 1, 11)
w ( 1, 1)
w(1, 1)
w (1, 1)
w (1, 1)
w
=
(7.23)
a2 + a3 + 2a5 + a6 2 + 2a7 + 3a9 2 + a10 3 + 3a11 2
w ( , ) =
(7.24)
avec:
1
0
1
0
0
[ R ] = 1
0
0
0
0
9782100544639.indb 115
1
1
0
1
1
1
0
1
1
0
1
1
1
1
1
1
2
0
1
2
1
0
2
1
0
1
2
1
1
2
1
1
2
1
1
1
3
0
1
3
1
0
3
1
0
1
3
1
1
3
0
1
1
1
1
1
0
1
2
1
1
1
2
1
0
1
3
1
1
1
1
0
1
0
1
1
1
1
1
2
0
1
0
2
1
2
1
1
1
2
1
3
0
1
0
3
1
3
1
1
1
0
0
1
1
1
2
0
0
2
2
1
1
2
3
0
0
3
3
1
1
1
1
1
3
1
1
3
1
1
3
115
09/02/10 14:56
7 lments isoparamtriques
1 ( + 1) (1 + ) ( 2 + 2 2)
8
2
1
( 1) ( + 1) (1 + )
2
1
( 1) (1 + ) ( + 1)
8
2
2
8 ( 1) (1 + ) ( + + 2)
2
1
( + 1) ( 1) (1 + )
2
1
( 1) (1 + ) ( 1)
8
T
{N } =
1
2
2
(
)
(
)
(
)
+
+
+
( + 1) ( 1) ( 1)
8
2
1
(
)
(
)
(
)
1
+
1 ( + 1) ( 1) ( 2 + + 2 2)
8
2
1
( 1) ( + 1) ( 1)
8
2
1
(1 + ) ( 1) ( + 1)
(7.25)
116
9782100544639.indb 116
09/02/10 14:56
7 lments isoparamtriques
Plaque Q4 (C1)
(7.26)
T
Avec:
{ X }T = {1 2 2 2 2 3 3 3 3 3 2 2 3 2 2 3 3 }
wi
w
i
Cet lment ncessite donc quatre degrs de libert par nud {wi } = wi
2
wi
cequi explique lintroduction de la drive seconde par rapport et .
= a1 + a3 + 2 a4 + 2 a6 + a7 2 + 3a8 2
w ( , ) =
(7.27)
= a2 + a3 + 2 a5 + a6 2 + 2 a7 + 3a9 2
w ( , ) =
Dunod La photocopie non autorise est un dlit.
(7.28)
2w
w ( , ) =
(7.29)
9782100544639.indb 117
117
09/02/10 14:56
7 lments isoparamtriques
w(1,1)
w (1,1)
w (1,1)
w (1,1)
w( 1,1)
= [ R ] {a}
.
w(1, 1)
w (1, 1)
w (1, 1)
w (1, 1)
(7.30)
Avec:
1
0
0
1
0
0
0
[ R ] = 1
0
0
1
0
0
0
1
1
1
0
1
1
1
2
1
0
1
2
1
1
1
3
1
0
1
3
1
1
1
3
1
2
1
2
0
0
1
1
0
0
1
1
0
1
0
1
0
1
0
1
0
1
1
1
1
1
1
1
1
1
1
0
0
1
2
0
0
1
2
0
2
0
1
0
2
0
1
0
2
1
2
1
2
1
2
1
2
1
2
2
1
1
2
2
1
1
2
0
0
1
3
0
0
1
3
0
3
0
1
0
3
0
1
0
3
1
3
1
3
1
3
1
3
1
3
3
1
1
3
3
1
1
3
2
6
1
3
2
6
1
3
2
3
6
1
2
3
6
1
2
3
2
4
1
2
2
4
1
2
2
0
1
1
0
0
0
1
0
1
0
1
1
1
1
1
0
1
2
0
0
0
1
0
2
0
2
1
2
1
2
2
1
1
2
2
0
1
3
0
0
0
1
0
3
0
3
1
3
1
3
3
1
1
3
3
6
1
3
2
6
6
1
2
3
6
4
1
2
2
4
1
3
9
1
3
3
9
1
3
3
9
1
3
3
9
118
9782100544639.indb 118
09/02/10 14:56
7 lments isoparamtriques
9782100544639.indb 119
1 ( 2) ( + 1) 2 ( 2) (1 + )2
16
1 ( 1) ( + 1 )2 ( 2) (1 + )2
16
1
2
2
( 2) ( + 1 ) ( 1) (1 + )
16
1 ( 1) ( + 1) 2 ( 1) (1 + )2
16
2
2
1 ( + 2 ) ( 1 ) ( 2) (1 + )
16
2
2
1
16 ( + 1 ) ( 1 ) ( 2) (1 + )
1 ( + 2) ( 1) 2 ( 1) (1 + )2
16
1
2
2
( + 1) ( 1) ( 1) (1 + )
16
{N }T =
2
2
1
16 ( + 2) ( 1) (2 + ) ( 1)
1 ( + 1) ( 1) 2 (2 + ) ( 1)2
16
1
2
2
( + 2) ( 1) (1 + ) ( 1)
16
1 ( + 1) ( 1) 2 (1 + ) ( 1)2
16
2
2
1 ( 2) ( + 1 ) (2 + ) ( 1)
16
2
2
1
16 ( 1) ( + 1) (2 + ) ( 1)
1 ( 2) ( + 1) 2 (1 + ) ( 1)2
16
1
2
2
( 1) ( + 1) (1 + ) ( 1)
16
(7.31)
119
09/02/10 14:56
7 lments isoparamtriques
(7.32)
Ve
La relation (7.32) fait intervenir une matrice [ B ] dduite des relations (3.13) et
(4.4). Par dclinaison de (4.10), celle-ci peut galement tre dveloppe en une
srie de sous matrices [ Bi ] associes chacun des nuds de llment, soit:
{ } = [ B ] {qe } = B1 B2 . Bi . Bn {qe }
Ni
x
avec [ Bi ] = N
i
Ni
(7.33)
Ni
0
y
N i
0
z
. La dtermination de [ B ] ncessitera de faire
Ni
0
x
Ni
0
x
Ni Ni
z
y
Ni Ni Ni
le calcul des diffrentes drives partielles
,
,
associes chacun
x
y
z
0
des nuds. Ceci tant et comme vu en 7.2, les fonctions de forme N i sont dfinies
dans un repre unitaire ( , , ) et non dans le repre local ( x , y , z ) . La transformation gomtrique e issue de (7-1) permet nanmoins dcrire que :
Ni Ni x Ni y Ni z
=
x
y
z
Ni Ni x Ni y Ni z
x y z
(7.34)
Ni Ni x Ni y Ni z
=
x
y
z
120
9782100544639.indb 120
09/02/10 14:56
7 lments isoparamtriques
z Ni
Ni
z Ni
Ni
= [ J ]
y
y
Ni
z Ni
z
z
(7.35)
x
Ni
1 N i
(7.36)
= [J ]
y
Ni
Ni
Ni x
Ni x
=
Ni x
x
=
N i xi
i =1
z = N i zi
i =1
9782100544639.indb 121
n Ni
xi
i =1
n N
[ J ] = i xi
i =1
n N
i
xi
i =1
N1
N
[ J ] = 1
N1
Ni
yi
i =1
Ni
yi
i =1
Ni
yi
i =1
Ni
Ni
.
Ni
.
Ni
zi
i =1
Ni
zi
i =1
n
N
i zi
i =1
Nn x
1
.
Nn
.
xi
.
Nn
.
x
n
y1
.
yi
.
yn
(7.37)
z1
.
zi
.
zn
121
09/02/10 14:56
7 lments isoparamtriques
Pour calculer lintgrale sur le volume Ve et sachant que [ B ] est dsormais fonction des coordonnes unitaires, il reste effectuer un changement de variables entre
les repres ( x , y , z ) et ( , , ) . Ce changement de variables reprsente en fait les
relations gomtriques liant, via la transformation e , lments rel et de rfrence.
Ces relations sont obtenues partir du jacobien de la transformation. La matrice
rigidit lmentaire exprime dans le repre ( , , ) a donc pour expression:
[ke ] = [ B ]T [ H ] [ B ] dxdydz
1 1 1
[ B ] [ H ] [ B ] det [ J ] d d d
(7.38)
1 1 1
Avec:
dx
d
dy = [ J ] d et dxdydz = det [ J ] d d d
dz
d
En appliquant ces rsultats la matrice de rigidit de llment barre (4.13), on
saperoit que le jacobien ou son inverse reprsentent le rapport des dimensions des
lments rel et de rfrence. On a donc daprs le changement de variable tabli
en 7.2:
x
2
= 2 1 d = dx
(7.39)
L
L
ce qui donne comme prvu pour [ ke ] :
T
[ke ] = [ B ] [ H ] [ B ] dVe
Ve
= ES [ B ] [ B ] dx
0
1
1 1 2
= ES 2
d
1 2 2 L
1
2
1
1
1
2 ES 4
4 d = ES 1 1
=
L 1 1 1
L 1 1
4 4
1
(7.40)
Pour les forces de volume (4.17), la dmarche est similaire celle suivie pour la
rigidit lmentaire, le changement de base seffectuant toujours via le dterminant
du jacobien.
1 1 1
1 1 1
122
9782100544639.indb 122
09/02/10 14:56
7 lments isoparamtriques
Pour ce qui concerne les forces de surface (4.18), le calcul consistera exprimer la
quantit dSe en fonction des coordonnes normes des artes ou facettes des lments surfaciques ou volumiques.
f ( ) d = wi f (i )
Le calcul des matrices de rigidit des lments barre (7.40) et poutre (6.21) a permis dtablir celles-ci en fonction de paramtres tels que longueur, section, inerties,
etc. Ces matrices sont gnralement codes de manire explicite dans les logiciels
de calcul par lments finis.
Cette approche nest bien sr pas envisageable pour les lments surfaciques et
volumiques dont les formes et configurations varient systmatiquement. La solution numrique simpose alors. Ceci tant, il est indispensable de disposer dune
mthode permettant une intgration exacte des quantits recherches.
La mthode de Gauss qui est la plus couramment utilise, rpond cette condition. Bien que sa matrise ne soit pas fondamentale pour les utilisateurs courants,
il est nanmoins important de bien assimiler la notion de points dintgration ou
points de Gauss pour comprendre les diffrentes approches utilises pour le
postraitement des rsultats. Cest pourquoi, lintgration numrique 1 dimension
est dveloppe ci-aprs.
Soit une fonction f dune seule variable , lintgrale de f sur lintervalle [ a,b ]
pourra tre value numriquement par la mthode de Gauss en posant que:
(7.42)
i =1
a f ( ) d . On a donc:
f ( ) ( ) = a1 + a2 + ....... + a2r 2r 1
a ( ) d
(7.43)
a f ( ) d = a1 a d + a2 a d + ..... + a2r a 2r 1 d
9782100544639.indb 123
123
09/02/10 14:56
7 lments isoparamtriques
2
lorsque est pair,
1+
lquation (7-44) peut alors scrire sous la forme dun systme linaire de 2r quations:
2 = w1 + w2 + ........ + wr
0 = w11 + w2 2 + ........ + wr r
2
= w112 + w2 22 + ........ + wr r2
(7.45)
3
.
.
Sachant que
1 d
wi f (i ) = w1 f (1 ) + w2 f ( 2 )
1 f ( ) d =
i =1
Cette intgrale sera value de faon exacte si le polynme est de degr infrieur ou gal 2r 1 = 3 .
Les coefficients sont obtenus en rsolvant le systme linaire suivant:
2 = w1 + w2
0 = w1 1 + w2 2 w1 = w2 = 1
1
2
= w1 12 + w2 22 1 = 2 =
3
3
0 = w1 13 + w2 23
Tableau 7.1 Tables de Gauss-Legendre (1 dimension).
r
xi
wi
degr maximum
dupolynme intgrer
1
3
0
3
5
5/9
32 6 /5
7
1
1
+
2 6 6/5
3+2 6 /5
7
1
1
2 6 6/5
8/9
5
124
9782100544639.indb 124
09/02/10 14:56
7 lments isoparamtriques
1 1
f ( , ) d d = i j wi w j f (i , j )
(7.46)
1 1
9782100544639.indb 125
125
09/02/10 14:56
9782100544639.indb 126
09/02/10 14:56
8 lments de membrane
A
Les lments membranes sont destins modliser des structures planes travaillant
uniquement dans leur plan. Leurs nuds possdent uniquement deux degrs de
libert u et v . Ils ont une paisseur constante note e et suivent lhypothse
decontrainte plane dcrite au chapitre 3.3.4. La relation contrainte-dformation
issue de (3.27) scrit alors:
xx
E
yy =
2
1
xy
1
0 xx
0 yy = [ H ] { }
1
1 xy
0 0
(8.1)
1
0
E
1
Avec [ H ] =
0 .
2
1
1
0 0
Le vecteur dformation scrit quant lui en adaptant (3.13) au cas plan, soit:
xx x
{ } = yy = 0
xy
y
u
u
=
[
]
y v
v
x
(8.2)
9782100544639.indb 127
127
09/02/10 14:56
8 lments de membrane
1
[H ] =
(1 + ) (1 2 )
0
0
0
1 2
(8.3)
2
y
q = 1 107 N/m
2h
x
Y
X
6
h
9782100544639.indb 128
09/02/10 14:56
8 lments de membrane
Llment utilis est un quadrangle de type Q4. Ses fonctions de forme sont gales
celles tablies en (7.9), soit:
1
1
(1 ) (1 + ) N1 = (1 + ) (1 + )
4
4
1
1
N 3 = (1 ) (1 ) N 4 = (1 + ) (1 )
4
4
N2 =
lment
(8.4)
9782100544639.indb 129
N1
[ J ] = N
N2
N2
N3
N3
N 4 x1
x2
N 4 x3
x 4
y1
y2
y3
y4
(8.5)
129
09/02/10 14:56
8 lments de membrane
[J]=
[J]=
[J]=
x1
(1 ) x2
(1 + ) x3
x4
1 (1 + ) (1 + ) (1 )
4 (1 + ) (1 ) (1 )
1 x1 x2 x3 + x 4 + ( x1 x2 + x3 x 4 )
4 x1 + x2 x3 x 4 + ( x1 x 2 + x3 x 4 )
1 x1 x2 x3 + x 4 + ( x1 x2 + x3 x 4 )
4 x1 + x2 x3 x 4 + ( x1 x 2 + x3 x 4 )
y1
y2
y3
y4
(8.6)
y1 y2 y3 + y 4 + ( y1 y2 + y3 y 4 )
y1 + y2 y3 y 4 + ( y1 y2 + y3 y 4 )
y1 y2 y3 + y 4 + ( y1 y2 + y3 y 4 )
(8.7)
y1 + y2 y3 y 4 + ( y1 y2 + y3 y 4 )
Bien que la connectivit des lments rel et de rfrence soit identique, les coordonnes xi , yi correspondent bien celles de llment rel.
h
0 + (0)
h + h + ( 0 )
1
= 2
[ J ] =
h h h h
4 0 + (0)
+ + + + ( 0 ) 0
2 2 2 2
0
h
2
(8.8)
h2
Le dterminant et linverse de [ J ] sont donc respectivement gaux det [ J ] =
4
et:
h
1
2 1 0
(8.9)
[ J ]1 = h 2 2 h = 0 1
h
4
2
Ce qui permet de dduire les relations entre les
Ni Ni
Ni Ni
,
et
,
:
x y
Ni 2 Ni
=
x
h
Ni 2 Ni
=
h
y
(8.10)
9782100544639.indb 130
09/02/10 14:57
8 lments de membrane
1 1
[k1 ] = e [ B ]T [ H ] [ B ] det [ J ] d d
(8.11)
1 1
[k1 ] =
1 1
Ee
T
[B ]
2
1 1 1
(8.12)
1
0
h2
0 [ B ] d d
1
4
1
0
0
u
N1
= [ N ] {q1} =
v
0
0
N1
N2
0
0
N2
N3
0
0
N3
N4
0
N4
x
u1
v
1
u2
0 v2
(8.13)
N 4 u3
v3
u4
u4
[ B ] = [ ][ N ]
N
1
= 0
N
1
N1
y
N1
x
N2
x
N2
y
N2
y
N2
x
N3
x
0
N3
y
N3
y
N3
x
N4
y
N4
y
N4
(8.14)
[B ] =
1 2
4 h
(1 + ) 0 (1 + )
(1 )
0
0
0
(1 )
0
0
(1 )
0
(1 + ) (8.15)
(1 )
0 (1 + )
(1 + ) (1 + ) (1 ) (1 + ) (1 ) (1 ) (1 + ) (1 )
9782100544639.indb 131
131
09/02/10 14:57
8 lments de membrane
1 +
8
3
12
1 + 3
8
[k1 ] = Ee 2
1 3 +
12
1
8
1 3
8
1 +
8
3
6
1 3
8
6
1
8
3 +
12
1 + 3
8
3
12
3
12
1 3
8
3
6
1
8
6
1 + 3
8
3 +
12
1 +
8
1 + 3
8
6
1
8
3
6
1 3
8
3
12
1 +
8
3 +
12
3 +
12
1
8
6
1 3
8
3
6
1 +
8
3
12
1 + 3
8
1
8
3 +
12
1 + 3
8
3
12
1 +
8
3
6
1 3
8
6
1 + 3
8
3 +
12
1 +
8
3
12
1 3
8
3
6
1
8
1 3
8
3
12
1 +
8
3 +
12
(8.16)
1 + 3
8
1
8
3
6
1
1 1
,
3 3
1 1
,
3 3
u
1 1
,
3 3
1 1
,
3 3
4
Figure 8.3 Points de Gauss.
132
9782100544639.indb 132
09/02/10 14:57
8 lments de membrane
Pour = 1 3 , les deux termes k11 et k12 de la matrice de rigidit sont donc gaux
:
k11 =
Ee
1 2
1 1
1 + 2 1 1 + 2
Ee 3
4 + 2 4 d d = 1 2 6 = 1.05 1010 N/m
1 1
2
2
2
2
1
1
1
1
1 + 3 1
1 3 1
1+ 3
1+ 3
1 1
1 1
1 1
+
+
+ 1 1
2
2
4
4
4
4
Ee
=
2
2
1 2
1
1
1
1
1
+
1
1
3 1
3 1
3
3
+1 1
1 1
1 1
+
+
+ 1 1
2
2
4
4
4
4
2
2
2
2
1
1
1
1
+
1
1
1
1
+
+
3 1
3
3 1
3
1 1
1 1
1 1
+
+
+ 1 1
2
2
4
4
4
4
Ee
=
(8.17)
2
2
2
1
1
1
1
1
1 3 1
1 + 3 1
1 3
1 3
+1 1
1 1
1 1
+
+
+ 1 1
4
4
2
4
4
2
k12 =
Ee
1 2
1 1
Ee 1 +
1 + 1 + 1 +
= 3.9375 109 N/m
d d =
2 8
2
4
4
1
1 1
1
1
1
1
1+
1+
1
1+
1 +
1 +
3
3
3
3
1 1
+1 1
2 4 4
2 4 4
Ee
(8.18)
=
1
1
1
1
1 2
1
1
1
1+
1 +
1 +
3
3
3
3
+1 1
+1 1
2 4 4
2 4 4
9
= 3.9375 10 N/m
9782100544639.indb 133
133
09/02/10 14:57
8 lments de membrane
} dS1 .
Considrant uniquement la face charge ( = 1) de llment 1, seules les fonctions de forme N1 et N 4 interviennent dans le calcul de { f1s } , les charges nodales
associes aux nuds 2 et 3 tant forcment nulles. valant tout le temps 1 le long
h
de la face 1-4 et comme dy = d , on a:
2
F1x
F
1y
F2 x
F
s = 2y
f
{ 1 } F
3x
F3 y
F4 x
F4 y
F
1x = 1 N 1
F1 y 1 0
F
4x = 1 N 4
F4 y 1 0
qh
2
0
0
0
=
0
0
qh
2
0
134
9782100544639.indb 134
qh
qh
1
4 (1 + 1) (1 + ) 2 d = 2
0
1
F2 x 0
=
F2 y 0
F3 x 0
F3 y 0
qh
1
0 q h
qh
1
d = (1 + 1) (1 ) d = 2
N 4 0 2
4
2
0
1
0 q h
d =
N1 0 2
(8.19)
09/02/10 14:57
8 lments de membrane
[ K ] {Q } = [k1 ] {Q } = { f1s } = {F }
Avec:
(8.20)
[K ] =
3
6
1 +
8
12
1 + 3
Ee 8
1 2 3 +
12
1
1 3
8
1 +
8
3 1 + 3 3 + 1
12
8
12
8
3
6
1 3
8
1 3
8
3
6
1
8
1
8
3
6
1 3 3
8
12
1 +
8
1
8
1 3
8
3
6
1 +
8
3
12
3 + 1 + 3 3
12
8
12
1 +
8
3
6
1 3
8
1 + 3 3 +
8
12
1 +
8
3 1 3
12
8
3
6
3
12
3 + 1 + 3
12
8
1 +
8
1 3 + 1 + 3
8
12
8
1 + 3 33 +
12
8
1
8
1 3
8
3
12
1 +
8
3 +
12
(8.21)
1 + 3
8
6
1
8
3
6
9782100544639.indb 135
3
6
Ee 1 +
1 2 8
1 +
8
3
6
1 + 3
8
qh
6 U
1
2
1 + 3
V1 = 0
8
3 U 4 qh
2
6
(8.22)
135
09/02/10 14:57
8 lments de membrane
8.1.5 v = 0
Dplacements
1
2
1
Ee
8
1
8
1
2
1
8
qh
0
2
U
1
1
V1 = 0
8
1 U 4 qh
2
2
(8.23)
qh
U1 = Ee = 0.476 mm
V1 = 0
qh
U 4 =
= 0.476 mm = U1
Ee
(8.24)
En injectant les dplacements (8.24) dans (8.20), on obtient pour les ractions:
136
9782100544639.indb 136
09/02/10 14:57
8 lments de membrane
E e 3
1 3
3 +
U1 +
V1 +
U 4
2
8
12
(1 ) 12
qh
EeU1
U V U
R X 2 = Ee 1 + 1 4 =
=
8
4
2
2
4
E e 1 + 3
1 +
RY 2 =
U
+
V
+
U
1
1
4
6
8
(1 2 ) 8
RX 2 =
U
U
RY 2 = Ee 1 + 0 + 4 = 0
8
8
(8.25)
E e 3 +
3
1
RX 3 =
V +
U
U +
(1 2 ) 12 1 8 1 12 4
qh
EeU1
U V U
R X 3 = Ee 1 1 4 =
=
8
4
2
2
4
E e 1
3 +
1 3
RY 3 =
U
+
V
+
1
1
U 4
1
8
12
8
(
)
U
U
RY 3 = Ee 1 + 0 + 4 = 0
8
8
Les dformations pouvant tre calcules grce aux expressions (4.10)et (8.15),
u1
v
1
u2
v2
(8.26)
{ } = [ B ] {q1} = [ B ]
u3
v3
u4
u4
On dduit les dformations et contraintes dans llment 1 en appliquant (8.2)et
(8.1) :
9782100544639.indb 137
u1
xx
h
(1 + ) u1 + (1 ) u4
2 1
{ } = yy =
(1 + ) v1
= 0 (8.27)
h 4 (1 + ) u + (1 + ) v (1 + ) u 0
1
1
4
xy
137
09/02/10 14:57
8 lments de membrane
8.1.6 v = 1/3
q
xx 100 MPa
{ } = yy = 0 = 0
xy 0
(8.28)
Dplacements
3 1
Ee
16 2
1
0
16
Ce qui donne pour les dplacements:
1
qh
16 U
1
2
0 V1 = 0
1 U 4 qh
2
2
28 qh
U1 = 27 Ee = 0.494 mm
7 qh
= 0.185 mm
V1 =
18 Ee
47 qh
U 4 = 54 Ee = 0.414 mm
(8.29)
(8.30)
138
9782100544639.indb 138
09/02/10 14:57
8 lments de membrane
Ractions
1 3
3 +
U1 +
V1 +
U 4
8
12
13
Ee
R X 2 = [ 10U1 + 0 8U 4 ] = qh = 5.417 MN
32
24
E e 1 + 3
1 +
RY 2 =
U1 + V1 +
U 4
2
6
8
(1 ) 8
E e 3
(1 2 ) 12
RX 2 =
9 Ee V1 U 4
= 1.389 MN
0+ +
8 18 6
(8.31)
E e 3 +
1
3
V +
U
=
U +
(1 2 ) 12 1 8 1 12 4
RY 2 =
RX 3
9 Ee 2U1 V1 5U 4
= 4.583 MN
8 9
6
18
E e 1
3 +
1 3
=
U 1 +
V1 +
2
12
8
(1 ) 8
RX 3 =
RY 3
RY 3 =
U 4
9 Ee U1 2V1
+ 0 0
8 6
9
La raction R X 3 est encore une fois multiplier par deux en raison de la symtrie.
Celle au nud 3 dans la direction Y est non nulle. Cependant, cette valeur sannule
avec celle de llment limitrophe qui grce la symtrie permet de bien trouver
RY 3 gale zro.
Contraintes
9782100544639.indb 139
xx
(1 + ) u1 + (1 ) u4
2 1
{ } = yy =
(1 + ) v1
h 4 (1 + ) u + (1 + ) v (1 + ) u
1
1
4
xy
103 + 9
108
7 (1 + )
q
=
Eh
36
( 4 + 7 3 )
36
(8.32)
139
09/02/10 14:57
8 lments de membrane
( 96 + 9 7 )
96
xx
q ( 40 + 9 63 )
{ } = yy =
288
e
xy
( 4 + 7 3 )
96
(8.33)
= 1; = 1
xx 116.7 MPa
{ } = yy = 38.9 MPa
14.6 MPa
xy
= 1; = 1
xx 102 MPa
{ } = yy = 4.86 MPa
8.33 MPa
xy
= 0; = 0
1
xx
100 MPa
q 40
{ } = yy =
= 13.89 MPa
e 288 4.17 MPa
xy
4
96
= 1; = 1
xx 97.9 MPa
{ } = yy = 32.6 MPa
xy
(8.34)
= 1; = 1
xx 83.3 MPa
{ } = yy = 11.1 MPa
6.25 MPa
xy
9782100544639.indb 140
09/02/10 14:57
8 lments de membrane
q = 1 107 N/m
Y
X
4
x
y
y
1
x
2
h
Les fonctions de forme de llment de rfrence associ sont daprs (7.18) gales :
N1( , ) = 1
N 2 ( , ) =
N ( , ) =
3
(8.35)
3 (0,1)
v
u
1
(1,0)
2
9782100544639.indb 141
141
09/02/10 14:57
8 lments de membrane
N1
[ J1 ] = N
N2
N2
N3
x1
x
N3 2
x
4
x1
y1
1 1 0
x
y2 =
1 0 1 2
x4
y 4
y1
y2 (8.36)
y 4
y2 y1 h 0
=
et det [ J1 ] = h 2
y 4 y1 0 h
(8.37)
ce qui donne:
x 2 x1
4 x1
[ J1 ] = x
x
2
y3
y4
y2
(8.38)
soit:
x 4 x3
2 x3
[ J2 ] = x
y 4 y3 h 0
=
et det [ J2 ] = h 2
y2 y3 0 h
(8.39)
1 1 0
h 0 1
(8.40)
do finalement:
[ J1 ]1 = [ J 2 ]1 =
Ni Ni
Ni Ni
,
et
,
sont gales :
x y
Ni 1 Ni
x = h
N
1 Ni
i
=
h
y
(8.41)
On notera par ailleurs quen (8.37) et (8.39) le dterminant du jacobien est gal
deux fois laire du triangle ce qui constitue une des particularits de llment
triangulaire trois nuds.
142
9782100544639.indb 142
09/02/10 14:57
8 lments de membrane
N2
x
N1
y
N1
x
N2
y
N2
x
N2
y
N3
x
0
N3
y
N3
y
N3
x
(8.42)
N
1
[B ] = 0
N1
y
1 0 1 0 0 0
1
= 0 1 0 0 0 1
h
1 1 0 1 1 0
[ B1 ]
[k1 ] =
Les
1
0
0 [ B1 ] det [ J1 ] d d
1
0 0
2
T
eEh 2 1 1 1 1
B
[
]
1
(1 2 ) h h 0 0
(8.43)
Ni
Ni
et
tant indpendants de , , il est possible de considrer chacun
des termes de la matrice [ k1 ] comme une constante et dintgrer une seule fois le
produit d d , soit:
1 1
0 0
1
0
d d = [ ]
0
1
2
d = (1 ) d = = (8.44)
0
2 0 2
do finalement:
3
1 +
2
Ee
[k1 ] =
2
4 (1 ) 1 +
1 +
1 +
3
2
1 +
1 +
2
2
2
2
0
0
2
1 +
1 +
0
1
1
0
1 +
1 +
0
1
1
0
2
2
2
(8.45)
0
0
qui est aussi gale [ k2 ] puisque les deux lments sont semblables.
9782100544639.indb 143
143
09/02/10 14:57
8 lments de membrane
cos
sin
0
[ Re ] = 0
sin
cos
0
0
0
0
0
0
cos
sin
0
0
0
0
sin
cos
0
0
0
0
0
0
cos
sin
0
0
0
0
sin
cos
(8.46)
lment 1: q = 0
[ R1 ] = [ I ] [ K 1 ] = [ R1 ]T [k1 ] [ R1 ] = [k1 ]
(8.47)
[ R2 ] = [ I ] [ K 2 ] = [ R2 ]T [k2 ] [ R2 ] = [k2 ]
(8.48)
lment 2: q = p
} dS2 .
La face 3-4 de llment 2 tant la seule charge, seules les fonctions de forme
N1 et N 2 seront concernes pour le calcul. En posant le changement de base
dx = h d , on obtient:
0 0
F3 x
qh
N1 0
qh
F
N
0
0
1
1
2 2
3y =
F4 x 0 N 2 0 q h d = 0 0
=
{ f 2s } = F4 y 0 N 2
qh qh
2 2
0
F2 x 0
=
F2 y 0
(8.49)
144
9782100544639.indb 144
09/02/10 14:57
8 lments de membrane
0
qh
2
0
{F2s } = qh
2
0
(8.50)
K 11
[ K 1 ] = K 21
K 41
K 12
K 22
K 42
K 33
K 14
K 24 et [ K 2 ] = K 43
K 23
K 44
K 34
K 44
K 24
K 32
K 42
K 22
(8.51)
( K 11 )1
(K )
K
=
[ ] 021 1
( K 41 )1
( K 12 )1
( K 14 )1
( K 22 )1 + ( K 22 )2 ( K 23 )2 ( K 24 )1 + ( K 24 )2
(8.52)
( K 32 )2
( K 33 )2
( K 34 )2
( K 42 )1 + ( K 42 )2 ( K 43 )2 ( K 44 )1 + ( K 44 )2
do:
[K ] =
Dunod La photocopie non autorise est un dlit.
3
1 +
Ee
1 +
4 (1 2 ) 0
0
1 +
1 +
2
1 +
3
2
1 +
2 2 + 1
0
1 +
0
1 + 2
0
1 +
2
0
1 +
2
1 +
0
1 + 2
2 2 + 1
0
0
0
1 +
2
3
1 +
2
1 +
0
1 +
2
0
1 +
2
1 +
0
2 + 1
2 1 + 2
0
1 +
2
1 +
3
2
1 +
2 1 + 2
0
1 +
0
2 + 1
(8.53)
9782100544639.indb 145
145
09/02/10 14:57
8 lments de membrane
2 2
Ee
1 + 1 +
[K ] =
2
0
4 (1 ) 0
0
0
1 + 1 +
2
2
2
2
3
0
1 +
1 +
0
1 +
1 +
1 +
0
3
2
2
1 +
0
0
0
1 +
2
3
1 +
2
1 +
0
0
1 +
2
1 +
3
2
1 +
1 +
1 +
0
1 +
2
2
3
0
2
2
1 +
0
(8.54)
11 +
1 +
0
Ceci tant, il aurait t tout fait possible dutiliser les matrices boolennes dassemblagevues en 4.2.4. Il suffirait alors de poser que [ A1 ] et [ A2 ] sont respectivement gales :
1 0 0 0 0 0 0 0
0 0 0 0 1 0 0 0
0 1 0 0 0 0 0 0
0 0 0 0 0 1 0 0
0 0 1 0 0 0 0 0
0 0 0 0 0 0 1 0
[ A1 ] = 0 0 0 1 0 0 0 0 [ A2 ] = 0 0 0 0 0 0 0 1 (8.55)
0 0 0 0 0 0 1 0
0 0 1 0 0 0 0 0
0 0 0 0 0 0 0 1
0 0 0 1 0 0 0 0
et de calculerensuite:
[ K ] = [ A1 ]T [ K 1 ] [ A1 ] + [ A2 ]T [ K 2 ] [ A2 ]
(8.56)
Les nuds 1,2 et 4 tant bloqus dans les directions horizontale et verticale, le
vecteur charges {F } en repre global scrit:
RX 1
R
Y1
RX 2
RY 2
(8.57)
{F } = 0
qh
RX 4
qh
RY 4
2
146
9782100544639.indb 146
09/02/10 14:57
8 lments de membrane
3 1 + U 3
E
(8.58)
1 + 3 V = qh
2
4 (1 )
3
2
ce qui donne pour les dplacements au nud 3:
A
(8.59)
qh (1 + )
= 0.201mm
U 3 =
4 Ee
V = qh (1 + ) ( 3 ) = 0.418 mm
3
4 Ee
Ractions
Les ractions aux nuds bloqus peuvent tre obtenues directement partir de la
matrice de rigidit du systme (8.54) et de la relation:
0
0
0
0
RX 1
0 0 qh (1 )
RY 1
0
0
RX 2
0 0 qh ( 3 + ) 875000 N
RY 2
0 0
4125000 N
8
0 = [ K ] U
=
(8.60)
=
0
0
qh
3
0
0
0
V3 2
R X 4
0 0 qh (1 ) 875000 N
R
0 qh
5875000 N
8
Y4
2 qh ( 5 )
8
147
9782100544639.indb 147
09/02/10 14:57
8 lments de membrane
Contraintes
Les nuds 1, 2 et 4 de llment 1 tant tous bloqus, les contraintes dans cet
lment sont bien videmment toutes nulles. Concernant llment 2 et afin dvaluer les dplacements dans son repre local, il est tout dabord ncessaire de faire le
changement de base inverse suivant:
(8.61)
{q2 } = [ R2 ] {Q2 } = [ I ] {Q2 } = {Q2 }
Avant de calculer les contraintes proprement dites, il faut valuer les dformations
au moyen de la relation (4.10) en posantque { } = [ B ] {q2 } , [ B ] correspondant
dans ce cas lexpression trouve en (8.42). On a:
u3
v
3
1 0 1 0 0 0
1
u4
(8.62)
{ } = 0 1 0 0 0 1
v4
h
1 1 0 1 1 0
u
2
v
2
do les contraintesdans llment:
q (1 )
4e 17.5 MPa
1
xx
0
q (3 + )
1
0
=
{ 2 } = yy =
{
}
1
e
4
17.5 MPa
xy
q (1 )
0 0
4e
Y1
RX 2
RY 2
(8.64)
{F } = R X 3
qh
RY 3 2
qh
148
9782100544639.indb 148
09/02/10 14:57
8 lments de membrane
3
E
4 (1 2 ) 0
0
0 U 4
= qh
3 V4
2
(8.66)
2qh (1 2 )
= 0.321mm
V
=
Ee ( 3 )
(8.65)
9782100544639.indb 149
qh
0
qh
1111111 N
RX 1
0 0
3703703 N
Y
1
1
qh
+
(
)
2407407 N
RX 2
0 0
2
3
(
)
0
RY 2
0 0
0
(8.67)
=
= [K ]
=
R X 4
0 qh qh (1 ) 1296296 N
RY 4
0 2
6296296 N
2 ( 3 )
0
0
U 4 0 qh ( 2 )
0
V4 qh
0
149
09/02/10 14:57
8 lments de membrane
Contraintes
Le repre local de llment 1 tant confondu avec le repre global, on sait que
{q1} = [ R1 ] {Q1} = [ I ] {Q1} = {Q1} do pour les contraintes dans cet lment:
u1
v
1
xx
1 0 1 0 0 0
u
1
2
{1} = yy = [ H ] { } = [ H ] 0 1 0 0 0 1 (8.68)
v
h
1 1 0 1 1 0 2
u
xy
4
v4
2q
e ( 3 )
xx
22 MPa
2q
do { 1 } = yy =
= 74 MPa .
e
3
(
)
xy
0
La dmarche de calcul tant identique celle tablie en 8.2.6, les contraintes dans
llment 2 sont gales:
u3
v
3
xx
1 0 1 0 0 0
1
u4
{ 2 } = yy = [ H ] { } = [ H ] 0 1 0 0 0 1 (8.69)
v
h
1 1 0 1 1 0 4
u
xy
2
v
2
xx 0 0
soit { 2 } = yy = 0 = 0 .
q 1 26 MPa
)
xy (
e ( 3 )
8.2.8 Proprits de llment triangle
Les dplacements obtenus aux chapitres 8.2.6 et 8.2.7 mettent en vidence une
dissymtrie de comportement trange, les deux modlisations devant amener normalement aux mmes rsultats. Ceci est d au fait que llment triangle ne rpartit
pas uniformment les rigidits. Pour viter cet inconvnient, il est donc conseill
150
9782100544639.indb 150
09/02/10 14:57
8 lments de membrane
de raliser des maillages le plus symtriques possibles ou dfaut dutiliser des lments de petites tailles.
Considrant maintenant le modle prcdent appuy uniquement aux nuds 1
et 2, on saperoit mme que llment triangle introduit des dissymtries sur des
modles symtriques. Pour calculer ces nouveaux dplacements, il suffit de reprendre
la matrice de rigidit [ K ] (8.54) et de rsoudre le systme [ K LL1 ] {Q L1 } = {Fext }
avec:
3
1 +
Ee
1 +
4 (1 2 ) 0
0
1 +
1 +
3
2
1 +
0
0
1 +
2
2
2
3
0
1 +
1 +
0
1 +
1 +
1 +
0
3
2
2
1 +
0
0
0
1 +
2
3
1 +
2
1 +
0
0
1 +
2
1 +
3
2
1 +
1 + 2
1 + 2
0
1 +
1 +
0
2 11 + (8.70)
2 1 +
3
0
0
3
[K ] =
[ K LL1 ]
0
qh
4qh (1 + )
= 0.0989 mm
U 3 =
Ee ( 2 2 7 )
qh (1 + ) ( 2 4 + 7 )
V3 =
= 0.486 mm
Ee ( 2 2 7 )
[ K LL1 ] {QL1} = {Fext }
2
U = 2qh (1 ) = 0.0346 mm
4 Ee ( 2 2 7 )
qh ( + 7 ) (1 2 )
= 0.421mm
V4 =
Ee ( 2 2 7 )
9782100544639.indb 151
(8.71)
151
09/02/10 14:57
8 lments de membrane
La mme structure modlise avec un lment quadrangle Q4 donnerait, a contrario de llment triangle, des rsultats symtriques. Pour ce faire, on reprend la
matrice de rigidit trouve en (8.21):
[K ] =
3
6
1 +
8
3
12
1 + 3
Ee 8
1 2 3 +
12
1
6
1 3
8
1 +
8
3
6
1 3
8
6
1
8
3 +
12
1 + 3
8
3
12
3
12
1 3
8
3
6
1
8
6
1 + 3
8
3 +
12
1 +
8
1 + 3
8
6
1
8
3
6
1 3
8
3
12
1 +
8
3 +
12
3 +
12
1
8
6
1 3
8
3
6
1 +
8
3
12
1 + 3
8
1
8
3 +
12
1 + 3
8
3
12
1 +
8
3
6
1 3
8
6
1 + 3
8
33 +
12
1 +
8
3
12
1 3
8
3
6
1
8
1 3
8
3
12
1 +
8
3 +
12 (8.72)
1 + 3
8
6
1
8
3
6
[ K LL2 ]
et on rsout [ K LL 2 ] {Q L 2 } = {Fext } do:
152
9782100544639.indb 152
09/02/10 14:57
6qh (1 2 )
U
= 0.0956 mm
=
3
Ee ( 6 2 + 9 )
qh ( 9 ) (1 2 )
V3 =
= 0.462 mm
Ee ( 6 2 + 9 )
2
U = 6qh (1 ) = 0.0956 mm
4 Ee ( 6 2 + 9 )
qh ( 9 ) (1 2 )
= 0.462 mm
V4 =
Ee ( 6 2 + 9 )
(8.73)
8 lments de membrane
9782100544639.indb 153
153
09/02/10 14:57
8 lments de membrane
Bien videmment, la diminution de la taille des lments triangulaires fera disparatre peu peu la dissymtrie, lidal tant bien sr de raliser un maillage le plus
homogne possible (cf. figure 8.13).
Dans le cas contraire et comme il est possible de le constater la figure 8.12, les
rsultats pourront tre parfois assez curieux, ceci se traduisant dans notre cas par
des diffrences entre les dplacements latraux et verticaux. Bien videmment, tout
rentrera dans lordre avec un maillage plus homogne.
En consquence, lutilisation de triangles htrognes peut perturber assez sensiblement les rsultats. En effet et sachant quen thorie, le dplacement vertical de cette
membrane doit tre gal :
qh
v qh he
=
= =
v =
= 0.476 mm 1
h
E
E
Ee
1. On nglige dans ce calcul leffet d au coefficient de Poisson. De ce fait, on retrouve le rsultat (8.24).
154
9782100544639.indb 154
09/02/10 14:57
8 lments de membrane
0,1 m
0,2 m
0,05 m
A
q
Y
X
9782100544639.indb 155
155
09/02/10 14:57
8 lments de membrane
Dans notre cas dtude, cette modlisation se rduira la cration dune surface
rectangulaire de 0.6 par 0.2 m et la soustraction dun disque de 0.05 m de rayon.
Cependant, cette squence pourra varier en fonction du logiciel utilis, certains
prenant en charge automatiquement la cration de la plaque troue ds gnration
du rectangle et du cercle. Finalement et comme indiqu la figure 8.16, le modle
comportera une surface troue et cinq lignes, quatre pour le rectangle, une pour
le cercle.
8.3.2 La problmatique du maillage
Pour piloter le mailleur plus prcisment, une deuxime technique consiste dfinir un nombre de divisions ou des tailles minimales dlments sur les lignes gomtriques. Le mailleur suit alors ces informations tout en essayant dassurer une
certaine progressivit dans les tailles dlments. Ainsi et pour obtenir le maillage
de la figure 8.18, les lignes verticales, horizontales et le cercle ont t scindes respectivement en 12, 40 et 26 divisions.
Il est galement possible dimposer des tailles dlments sur certaines zones tout
en laissant le mailleur grer le complmentaire avec une taille minimale dlment.
Enfin, cette approche automatique a pour grand intrt dautoriser la modification, laugmentation ou la diminution des tailles de mailles volont.
156
9782100544639.indb 156
09/02/10 14:58
9782100544639.indb 157
8 lments de membrane
157
09/02/10 14:58
8 lments de membrane
Lapproximation polynomiale dcrite au chapitre 4.1, base de la mthode des lments finis fait que la solution obtenue nest pas forcment exacte. En dautres
termes, les fonctions de forme des lments utiliss reconstituent la forme des
dplacements de la structure avec plus ou moins de prcision. La solution exacte
est donc sense tre approche lorsque la taille des lments diminue.
Nanmoins, cette convergence monotone du champ de dplacement ne peut tre
effective qu condition de satisfaire la condition de mouvement de corps rigide
(MCR), au critre de la dformation constante, et de garantir la continuit des
dplacements sur llment et aux frontires de celui-ci (cf. 7.2). Dans ce cas,
llment est qualifi de conforme.
A contrario, les lments non conformes qui nassurent pas la continuit aux frontires ont une convergence lente et non assure. Cependant, certains dentre eux
ont des performances parfois suprieures leurs homologues conformes. Ceci doit
tre vrifi par la technique du patch test. De plus, [3] prcise que llment doit
possder une base polynomiale complte jusqu lordre n pour garantir la convergence des drives de u, v ou w dordren.
8.3.4 Discontinuit du champ de contraintes-dformations
La continuit du champ de dplacement tant garantie dans la majorit des cas (les
lments utiliss sont gnralement conformes et convergents), il nen est pas de
mme pour les champs de dformations, de contraintes et defforts. Ceux-ci sont
en effet continus sur llment mais en aucune faon sur ses frontires. Dans ce cas
de figure, les quations dquilibre de volume (3.5) et de surface (3.6) ne sont pas
forcment satisfaites.
Deux modes de visualisation des contraintes, des dformations ou des efforts sont
gnralement disponibles dans les codes de calcul pour aider lutilisateur vrifier
lquilibre de son modle: les modes liss ou non liss1.
1. On utilise galement parfois le terme de moyenn.
158
9782100544639.indb 158
09/02/10 14:58
8 lments de membrane
( xx )i = e =1
(8.74)
Gnralement, ce calcul est effectu partir des valeurs calcules aux points de
Gauss. Bien videmment, ce mode de calcul ne permet plus destimer les discontinuits des contraintes, des dformations ou des efforts. Il a de plus comme principal inconvnient dcrter ces valeurs en cas de fortes discontinuits.
Le deuxime mode de visualisation en valeurs non lisses ou non moyennes
revient afficher lvolution des contraintes, des dformations ou des efforts sur
les lments pris un un. On peut se servir dans ce cas de la mthodologie de
calcul utilise pour les exemples 5 et 6 aux chapitres 8.1.6 ou 8.2.6. Certains
logiciels se limitent cependant une seule valeur calcule au centre de gravit de
llment.
2
( xx )i e
4
( xx )1 =
[( ) ]
e =1
xx 1 e
Points de Gauss
xx
au nud 1.
9782100544639.indb 159
159
09/02/10 14:58
8 lments de membrane
En labsence de solution connue, la mise au point dun modle de calcul par lments finis passe en premier lieu par une vrification de la convergence des dplacements. Dans un deuxime temps, on compare les solutions lisses et non lisses
pour valider les contraintes et estimer leurs discontinuits.
Pour valider notre modle, trois densits de maillage imposant respectivement 24,
32 et 48 divisions au niveau du trou circulaire ont t retenues. Choisissant dans
un premier temps un lment linaire Q4, les dplacements obtenus avec le logiciel
Abaqus se rsument aux graphiques pages suivantes.
160
9782100544639.indb 160
09/02/10 14:58
8 lments de membrane
9782100544639.indb 161
161
09/02/10 14:58
8 lments de membrane
162
9782100544639.indb 162
09/02/10 14:58
8 lments de membrane
9782100544639.indb 163
163
09/02/10 14:58
8 lments de membrane
On retiendra donc une valeur de 0.36 mm pour le dplacement maximal. Laffichage en mode non liss met en vidence lvolution des discontinuits de
contraintes. Llment utilis tant un Q4, lvolution des contraintes sur celui-ci
est linaire (cf. 8.1.6). On reconstitue en fait le champ de contraintes partir
de tronons de droites do une convergence plus lente que pour son homologue
quadratique (Q8). En effet et en utilisant ce type dlment, on obtient pour les
trois types de discrtisation:
164
9782100544639.indb 164
09/02/10 14:58
8 lments de membrane
9782100544639.indb 165
165
09/02/10 14:58
8 lments de membrane
166
9782100544639.indb 166
09/02/10 14:58
8 lments de membrane
On saperoit dans ce cas que les contraintes lisses et non lisses sont tout fait
quivalentes, les variations entre les deux tant ngligeables. On notera galement
que la continuit des contraintes non lisses est quasi parfaite quelle que soit la
densit de maillage.
La contrainte maximale retenue sera donc comprise entre 433 (Q4) et 438 MPa
(Q8), valeur tout fait comparable celle calcule partir du coefficient de concentration de contrainte ( K = 2.16 ) issu des courbes de Peterson (cf. [4]).
0,6 m
B
qh
1 10 6 0,2
= 200MPa
A
A = A' = K AB =
= 2,16 AB = 432 MPa
0,05 m
h = 0,2 m
d = 0,1 m
AB =
Finalement, les deux types de maillage Q4 et Q8 donnent des rsultats tout fait
similaires, le deuxime convergeant plus rapidement du fait de sa formulation mais
avec bien sr un nombre de nuds plus important nombre de divisions gal.
Tout est en fait une affaire de compromis entre rapidit de convergence et temps
de calcul.
Maillage en triangles
9782100544639.indb 167
167
09/02/10 14:58
8 lments de membrane
168
9782100544639.indb 168
09/02/10 14:58
8 lments de membrane
9782100544639.indb 169
169
09/02/10 14:58
8 lments de membrane
170
9782100544639.indb 170
09/02/10 14:58
8 lments de membrane
9782100544639.indb 171
171
09/02/10 14:58
8 lments de membrane
172
9782100544639.indb 172
09/02/10 14:58
8 lments de membrane
9782100544639.indb 173
173
09/02/10 14:58
8 lments de membrane
Lutilisateur ayant gnralement peu de rfrences pour valider son modle, celui-ci
sera immanquablement amen sinterroger sur la qualit des rsultats obtenus.
Pour ce faire, il devra dvelopper une mthodologie lui permettant de sassurer de
leur stabilit.
En pralable, deux options soffrent lui, soit un maillage en quadrangles, soit celui
en triangles. Gnralement et densit de maillage constante, ce dernier converge
plus lentement et gnre plus dlments ce qui est un lger handicap. Comme
indiqu prcdemment, il a cependant pour principal avantage de librer lutilisateur des oprations gomtriques1 ncessites par les mailleurs hexadriques. Il
sagit en fait dun compromis entre le temps pass pour la modlisation et celui
ncessit par le calcul.
Pour valider son modle, lutilisateur devra dans un premier temps contrler la
stabilit des dplacements. Les modles lments finis convergeant assez vite dans
ce cas2, une variation infrieure 5% entre deux maillages successifs constitue un
critre dacceptation.
0,362
0,361
0,360
Q4
0,359
Q8
T3
0,358
T6
0,357
0,356
0,355
24
32
40
Maillage (nombre de divisions)
48
1. Ces oprations ont pour but dobtenir un maximum de surfaces 4 cts ou de cubes 6 faces.
2. Les lments utiliss sont trs souvent conformes.
174
9782100544639.indb 174
09/02/10 14:58
Gnralement, les dplacements obtenus avec des lments quadrangles sont assez
rapidement reprsentatifs, ceux associs des triangles voluant plus lentement.
Ceci vrifi, il devra examiner les tracs des contraintes non lisses et statuer sur
leur continuit. Bien videmment, plus llment utilis est riche, plus cette continuit est tablie rapidement. Pour les lments linaires, celle-ci peut voluer assez
lentement et contraindre1 lutilisateur se contenter de maillages de densit insuffisante. Dans ce cas, il sera prfrable de retenir les contraintes non lisses celles
lisses, ces dernires pouvant tre altres ventuellement par de fortes variations
entre lments.
Une autre approche consiste comparer les maximums des contraintes lisses et
non lisses puis vrifier leur convergence. Ceci tant, on se mfiera dans ce cas
de figure des discontinuits gomtriques ou des changements de matriaux qui
peuvent provoquer des concentrations de contraintes non convergentes2.
8 lments de membrane
440
430
420
410
Q4
400
Q8
390
T3
380
T6
370
Objecf
360
350
340
330
320
24
32
40
48
9782100544639.indb 175
175
09/02/10 14:58
8 lments de membrane
450
440
430
Q4
Q8
420
T3
T6
410
Objecf
400
390
380
24
32
40
48
h = 0,5 m
B
b = 0,2 m
L/2
L = 10 m
Figure 8.38 Exemple 8: poutre console.
176
9782100544639.indb 176
09/02/10 14:58
8 lments de membrane
x
2
ES
ES
0
0
0
0
L
0
qL
EI
EI
EI
EI
12
6
12
6
0
u
H
0 3
L3
L2
L
L2 1 2 1
V
v
6 EI
4 EI
6 EI 2 EI 1 qL2 1
0
0
M
L2
L
L2
L 1 = 12 + 1 (8.75)
ES
u2
0
ES
0
0
0
0
0
L
L
v2 qL 0
12 EI
6 EI
12 EI
6 EI 2 2 0
0
0
3
2
3
2
L
L
L
L
qL
12
EI
EI
EI
EI
6
2
6
4
0
0 2
2
L
L
L
L
Comme u2 est forcment nul du fait de lhypothse de petits dplacements, on
obtient aprs introduction des conditions dappui:
qL
12
6
L
v
2 2
EI
=
L3 6 L 4 L2 2 qL2
12
do les dplacements du nud 2 (point C):
avec I =
9782100544639.indb 177
qL4
= 28.57 mm
8EI
qL3
2 =
= 0.003809 rad
6 EI
v2 =
(8.76)
(8.77)
bh 3
= 0.002083 m 4 .
12
177
09/02/10 14:58
8 lments de membrane
qL
qL2
5
; M2 =
; S1, y = bh = 0.08333 m 2 .
2
12
6
avec G =
L
L3
+
v2 3 EI z GS1, y
=
L2
2
2 EI z
L2
2 EI z T2 28.65 mm
=
L M 2 0.003809 rad
EI z
(8.78)
E
.
2 (1 + )
Leffet de leffort tranchant est donc quasi ngligeable ce qui tait prvisible compte
L
Compte tenu de la gomtrie de cette poutre console, il est bien videmment prfrable de sorienter vers une modlisation en quadrangles. De plus et comme pour
lexemple prcdent, celle-ci sera effectue dans un premier temps avec des lments membranes Q41 ce qui reste tout fait conforme la thorie des poutres
puisque la variation des contraintes sur cet lment est linaire (cf. 8.1.6).
Enfin et afin de mettre en vidence linfluence du nombre dlments sur la prcision des rsultats, deux discrtisations diffrentes sont envisages verticalement et
horizontalement: 2, 4, 8 ou 16 lments suivant y; 20, 40, 80, 120 ou 160 lments par rapport x. Le nombre dlments fix horizontalement peut cependant
paratre lev. Pour lexpliquer, il est ncessaire de rappeler que llment membrane utilis na que deux degrs de libert de translation par nud ce qui est insuffisant pour modliser un problme de flexion tel que le ntre. En effet, les rotations
nodales qui en thorie des poutres sont issues de la drive de la ligne lastique sont
reconstitues dans notre cas par des dplacements diffrentiels.
1. Dpaisseur b.
178
9782100544639.indb 178
09/02/10 14:58
8 lments de membrane
30
20
2
4
15
10
16
Objecf
Dplacement VC (mm)
25
5
0
20
40
60
80
100
120
140
160
lanalyse des dplacements, on notera tout dabord une relative faible influence
de la discrtisation verticale, 4 ou 8 lments tant largement suffisants pour approcher les rsultats du modle filaire.
Figure 8.41 Dforme (Effel) avec lments Q4 8 divisions verticales et 160 horizontales.
Ceci tant et comme pour lexemple prcdent, on constate une convergence lente
de llment Q4 vers la valeur objective, celle-ci ntant atteinte moins de 1%
quau bout de 160 divisions horizontales et 8 verticales. Bien videmment, un
lment quadratique aurait permis davoir une bien meilleure prcision plus rapidement. Par exemple, un lment Q9 avec 40 divisions horizontales et 4 verticales
permet dobtenir28.6 mm en extrmit de console soit une erreur de 0.07 % :
9782100544639.indb 179
179
09/02/10 14:58
8 lments de membrane
14
12
2
10
16
Objecf
4
2
0
20
40
60
80
100
120
140
160
xx
Par contre celles associes au point A divergent pour atteindre une valeur nettement suprieure 60 MPa2. On est en fait en prsence dun effet local li lencastrement et d au coefficient de Poisson. En effet, le raccourcissement au niveau
de lencastrement est rendu impossible par le blocage des nuds ce niveau. Il
en rsulte la cration dune contrainte verticale qui via le coefficient de Poisson se
projette horizontalement et explique lexcs de contrainte constate.
Pour sen convaincre, il suffit de refaire le calcul avec un coefficient de Poisson nul.
1. Cette comparaison a t faite uniquement avec des contraintes non lisses.
2. 68.85 MPa pour 16 divisions verticales et 160 horizontales.
180
9782100544639.indb 180
09/02/10 14:58
8 lments de membrane
50,00
40,00
60,00
16
30,00
Objecf
20,00
10,00
0,00
20
40
60
80
100
120
140
Maillage horizontal (nombre d'lments)
xx
160
80
68,85
60,31
60
52,01 54,63
44,20 45,09
40
36,63 36,81
29,19 29,07
21,83 21,63
20
14,52 14,35
7,25 7,15
0
0,01 0,00
7,16
0
0,1
0,2 7,26 0,3
0,4
0,5
14,53 14,35
20
21,84 21,63
29,19 29,08
36,63 36,80
40
44,20 45,08
51,99
54,61
60
60,26
68,79
80
Hauteur h
0,3
0
9782100544639.indb 181
181
09/02/10 14:58
8 lments de membrane
182
9782100544639.indb 182
09/02/10 14:58
8 lments de membrane
{ } = [ H ] ({ } { 0 })
(8.80)
avec
{ 0 }T = { x T y T z T 0 0 0}
T
T
1
1
{ } { } dVe = {qe } [ B ]T [ H ] ([ B ] {qe } { 0 }) dVe (8.81)
2V
2V
e
do,
We =
T
1
1
{qe }T [ke ] {qe } {qe } [ B ]T [ H ] { 0 } dVe
2
2V
e
1
1
T
T
We = {qe } [ ke ] {qe } {qe } { f e0 }
2
2
(8.82)
Une fois calcules, ces forces viennent sajouter purement et simplement celles
associes aux charges appliques sur llment (4.16). Le schma de rsolution est
alors identique celui vu au chapitre 4.2.
Ces forces nodales sauto quilibrant avec les efforts calculs lors de la rsolution, il
est assez dlicat de vrifier la bonne entre des donnes1. En dautres termes, quelles
que soient les tempratures et coefficients de dilatation introduits, les ractions
sont systmatiquement nulles. Une solution pour rsoudre ce problme consiste
isoler une partie du modle et intgrer quand cest possible, les efforts afin de
vrifier que ceux-ci correspondent bien aux donnes introduites.
Ceci tant, deux cas extrmes peuvent se prsenter aprs calcul. Soit le modle est
compltement libre de se dformer et dans ce cas, la dformation thermique se
traduira uniquement par des dplacements sans apparition de contraintes.
1. Gnralement, on vrifie la correspondance entre charges appliques et ractions.
9782100544639.indb 183
183
09/02/10 14:58
8 lments de membrane
On notera dans ce cas que leffort normal est nul dans llment.
Inversement et si on applique la mme charge thermique une poutre bi-encastre
de mmes caractristiques, on nobtiendra aucun dplacement mais un effort normal de compression de:
0 = T N = xx S = E 0 S = E T S
N = 2.1 1011 1.2 105 20 0.5 0.2 = 5040 kN
soit une contrainte normale xx de compression dans llment gale 50.4 MPa.
T = 20 .
Il est donc trs important de ne pas perdre de vue ces aspects lorsque lon effectue
un calcul thermo-lastique car suivant la nature des appuis ou des rigidifications
lies la gomtrie de la structure, on peut voir apparatre des contraintes parfois
trs importantes dans certaines zones brides ou au contraire observer une dformation libre dans certaines autres.
184
9782100544639.indb 184
09/02/10 14:58
9 lments de plaque
A
9.1 Rappels sur les thories des plaques
9.1.1 Thorie des plaques minces
La thorie de Kirchhoff relative aux plaques minces qui revient ne pas prendre en
compte le cisaillement transverse, est applicable lorsque le rapport de la plus petite
de leurs dimensions sur lpaisseur est suprieur ou gal 201 (cf. [5]).
Relations moments-courbures
z, w( x, y)
y, v (x, y)
x, u (x, y)
1.
a
b
ou 20 .
e
e
9782100544639.indb 185
185
09/02/10 14:58
9 lments de plaque
w
<0
x
w
<0
y
(9.2)
La relation (9.2) est analogue lexpression (6.6) issue de la thorie des poutres
(cf. 6.1). La thorie de Kirchhoff sapparente en fait celle des poutres lances.
Le comportement des structures planes en flexion relevant de lhypothse de
contrainte plane (cf. 3.3.4), la relation entre contraintes et dformations normales
scrit daprs (3.27):
xx =
yy
E
E z 2w
2w
x + y =
2 + 2
2
2
1
1 x
y
E z 2w
E
2w
y
x
1 2 y 2
1 2
x 2
(9.3)
e
2
e xx z dzdy
M yy dx =
2
e
2
(9.4)
e yy z dzdx
2w
2w
M xx = D 2 + 2
y
x
2w
2w
M yy = D 2 + 2
x
y
(9.5)
186
9782100544639.indb 186
09/02/10 14:58
9 lments de plaque
E
1 2
z 2 dz =
Ee 3
.
12 (1 2 )
avec D =
e
2
v u
2w
xy = G xy = G + = 2G z
(9.6)
x y
x y
On sait par ailleurs que ces contraintes de cisaillement, intgres sur lpaisseur,
squilibrent avec un moment de torsion not M xy :
M xy =
e
2
e xy z dz
(9.7)
do
M xy =
Ge 3 2w
6 x y
(9.8)
soit en fonction de D :
M xy = D (1 )
2w
x y
(9.9)
Relations moments-contraintes
Les relations (9.3) et (9.5) ayant des termes communs, les contraintes normales
peuvent tre directement relies aux moments de flexion M xx et M yy par:
12 z M xx
xx =
e3
(9.10)
12 z M yy
yy =
e3
e
xx =
yy =
9782100544639.indb 187
6 M xx
e2
6 M yy
(9.11)
e2
187
09/02/10 14:58
9 lments de plaque
Lexpression (9.11) est de plus tout fait quivalente celle utilise pour la vrification des poutres en flexion simple. En effet, (9.10) peut tre aisment retrouve
partir de (8.79) en posant que linertie I est gale e 3 12 .
De mme et pour la contrainte de cisaillement, (9.6) et (9.8) permettent de
dduire:
12 z M xy
xy =
(9.12)
e3
e
soit pour z = :
2
6 M xy
xy = 2
(9.13)
e
Flexion des plaques sous charges latrales rparties
q
y
Ty +
T
Tx + x dx
x
T y
y
dy
Tx
Ty
dx
dy
Figure 9.3 Sommation des forces par rapport z.
T y
Tx
dx dy Tx dy + T y +
dy dx T y dx + q dx dy = 0
Tx +
x
y
do aprs simplification:
Tx T y
+
= q
x
y
(9.14)
188
9782100544639.indb 188
09/02/10 14:58
9 lments de plaque
q
y
M xy
M xx +
M xx
M xx
dx
x
M xy +
M xy
y
Tx
dy
dx
dy
Figure 9.4 Sommation des moments par rapport laxe y.
M xy
M xx
Tx dy dx = 0
dx dy M xx dy + M xy +
dy dx M xy dx T
M xx +
x
y
do aprs simplification:
M xx M xy
(9.15)
+
= Tx
x
y
M xy
x
M yy
y
= Ty
(9.16)
M xx M xy
+
x x
y
M xy M yy
+
+
y
y x
+q =0
(9.17)
Soit
2 M xy 2 M yy
2 M xx
2
+
+
= q
x 2
x y
y 2
(9.18)
9782100544639.indb 189
189
09/02/10 14:58
9 lments de plaque
4w
4w
4w q
2
+
+
= = ( w )
x 4
x 2 y 2 y 4 D
(9.19)
Lexpression (9.19) est en fait une extension de la relation (6.9) relative la thorie
des poutres lances.
9.1.2 Thorie des plaques paisses
y , v ( x, y )
a
x, u ( x, y )
{ f }
u
x
xx
y
v
= yy =
= z
y
y
xy v u
y x
+
y
x
x y
(9.21)
a b
1. 4 < ; < 20 .
e e
190
9782100544639.indb 190
09/02/10 14:58
9 lments de plaque
Supposant les dformations lies au cisaillement transverse constantes sur lpaisseur, il est alors possible dcrire que:
w
+y
yz y
{ c } = =
xz w +
x
(9.22)
de (9.20).
Lhypothse de contrainte plane nous permet donc de dduire:
Pour la flexion
{ f }
1
xx
0
y
Ez
= yy =
1
0
= H f f
1
xy
0 0
y x
2
+
y
x
{ }
(9.23)
Pour le cisaillement
+y
yz
1 0 y
E
{ c } = =
= [ H c ] { c }
xz 2 (1 + ) 0 1 w +
x
x
(9.24)
Lnergie de dformation lastique totale sera alors gale la somme des nergies
de flexion et de cisaillement:
Dunod La photocopie non autorise est un dlit.
W =
T
1
1
{ } { } dV = f
2V
2V
{ } { f } dV + k2 { c }
T
{ c } dV
(9.25)
9782100544639.indb 191
1
f
2 V
{ }
T
k
H f f dS dz + { c } [ H c ] { c } dS dz
2
{ }
(9.26)
191
09/02/10 14:58
9 lments de plaque
1
f
2 S
{ }
T
k
H 1 f f dS + { c } [ H 1c ] { c } dS (9.27)
2
{ }
1
0
Ee 3
H 1 f =
0
1
2
12 (1 )
1
0 0
avec
[ H 1c ] =
Eek 1 0
2 (1 + ) 0 1
y
Appuis simples
2a
Appuis simples
Appuis simples
Appuis simples
b
x
Figure 9.6 Exemple 9: plaque simplement appuye en flexion.
9782100544639.indb 192
09/02/10 14:58
9 lments de plaque
x (y = 0)
( w ) x =0 = 0
2w
2w
2w
2w
=
0
+
2
2
y 2 x =0
y 2 x =0
x
x
( M xx )x =0 = 0 = D
(9.28)
y (x = 0)
2
( M yy ) y =0 = 0 = D yw2 + xw2
y =0
2w
2w
=0
2 + 2
x y =0
y
(9.29)
( w ) y =0 = 0
Pour rsoudre lquation gnrale des plaques (9.19) avec ce type de conditions
dappui, une mthode consiste dcomposer les fonctions de dplacement et de
charge en sries de fonctions sinus. On pose alors que:
m x
n y
w ( x , y ) = amn sin
sin
2b
2a
m =1 n =1
(9.30)
m x
n y
q ( x , y ) = qmn sin
sin
2b
2a
m =1 n =1
(9.31)
n y
m x
sin
sin
16 q
2b
2a
w ( x, y ) = 6
2
D m =1 n =1
m2
n2
m n
+
( 2 a )2 ( 2b )2
(9.32)
M xx ( x , y ) = 6
m =1 n =1 ( 2a )2
( 2b )2
9782100544639.indb 193
n y
m x
sin 2 a sin 2b (9.33)
m2
n2
m n
+
( 2 a )2 ( 2b )2
193
09/02/10 14:58
9 lments de plaque
2w
2w
M yy ( x , y ) = D 2 + 2
x
y
n 2 2
16 q
m 2 2
+
M yy ( x , y ) = 6
m =1 n =1 ( 2b )2
( 2 a )2
n y
m x
sin 2 a sin 2b (9.34)
m2
n2
m n
+
( 2 a )2 ( 2b )2
194
9782100544639.indb 194
09/02/10 14:58
9 lments de plaque
Il est cependant important de prciser que les moments et autres efforts sont indiqus par unit de longueur le long de laxe considr et non autour de laxe comme
cela a t le cas lors de ltude des systmes de poutres.
9.2.2 Plaque rectangulaire de 2 m 4 m
w ( a,b ) = 0.01053 m
M xx ( x , y ) .
9782100544639.indb 195
195
09/02/10 14:58
9 lments de plaque
M yy ( x , y ) .
9.2.3 Maillage 2 2
196
9782100544639.indb 196
09/02/10 14:58
9 lments de plaque
w( , ) = N i wi
i =1
4
x ( , ) = N i iy
(9.35)
i =1
4
y ( , ) = N i ix
i =1
1
(1 + i ) (1 + i ) .
4
avec N i =
x1
(1 ) x2
(1 + ) x3
x4
y1
b
y2 2
=
y3
0
y 4
0
(9.36)
a
2
a 0
ab
1
det [ J ] =
[ J ] =
4
0 2
b
do:
Ni 2 Ni
=
x
a
Ni 2 Ni
=
b
y
(9.37)
(9.38)
1
f
2 S
{ }
T
k
H 1 f f dS + { c } [ H 1c ] { c } dS (9.39)
2
{ }
[k1 ] = k1 f + [k1c ] =
1 1
B f H 1 f B f det [ J ] d d +
1 1
9782100544639.indb 197
1 1
T
[ Bc ] [ H 1c ] [ Bc ] det [ J ] d d
(9.40)
1 1
197
09/02/10 14:58
9 lments de plaque
avec:
1
0
Ee 3
0
1
H 1 f =
2
12 (1 )
1
0 0
[ H 1c ] =
Eek 1 0
2 (1 + ) 0 1
{ }
f = B f {q1} = B f 1 B f 2 B f 3
{ c } = [ Bc ] {q1} = Bc1
Bc 2
Bc 3
w1
1x
0
0
1 y
2 Ni
.
B f 4
; B fi = 0
a
.
w4
0 2 N i
b
4 x
4 y
Bc 4 {q1}
2 Ni
a
; [ Bci ] =
2 Ni
b
N i
0
2 Ni
b
2 Ni
a
Ni
Llment dcrit prcdemment qui est bas sur la thorie de Henky-Mindlin, est
de par sa conception, adapt au calcul des plaques paisses. Il devrait donc se rvler peu efficace dans la mesure o notre plaque a un lancement gomtrique gal
100. Lexprience montre dailleurs que les rsultats donns par ce type dlment
se dtriorent trs rapidement lorsque lallongement augmente.
Il savre en fait que la part de rigidit relative au cisaillement [ k1c ] est trs survalue dans ce cas do la notion de blocage en cisaillement1. Pour pallier ce problme et afin de permettre une utilisation sans restriction de ce type dlment,
une technique consiste sous intgrer [ k1c ] . On parlera alors dintgration rduite.
De plus et dans la mesure o les termes relatifs la flexion k1 f sont intgrs de
manire exacte (2 2 points de Gauss pour notre Q4), elle est dite slective car
rserve uniquement au cisaillement.
1. Shear locking en anglais. Ce verrouillage en cisaillement doit tre analys avec prcaution dans la mesure
o celui-ci peut amener des erreurs grossires.
198
9782100544639.indb 198
09/02/10 14:58
9 lments de plaque
Dans notre cas, les termes associs x et y seront intgrs avec deux points de
Gauss dans la direction considre et un seul perpendiculairement. Par exemple,
ceci donnera pour le terme 2,2 de la matrice [ k1c ] :
1 1
Eek ab
N12 d d
2 (1 + ) 4 1 1
=1
3
2
2
2
Eek ab
1
1
2
2
1
2 1 1 +
1 + 0 ) + 1
1 + 0)
=
(
(
2 (1 + ) 4
3
3
4
(9.41)
k1c 22 =
{ f zs } = s [ N ] { f s } dS
2
1 1
N1
N2
N3
1 1
N 4 ( q ) det [ J ] d d
(9.42)
qab
.
4
[ K ]{q} = {F }
Le repre local de llment tant confondu avec le repre global, il na donc pas
lieu deffectuer de changement de base. Le systme peut donc tre rsolu directementen posant[K] = [k1, rduite] matrice de rigidit de llment 1 obtenue par
intgration rduite et {F } = { f zs } .
Plaque carre
Dans le cas de la plaque carre, lcriture de cette matrice [K] = [k1, rduite] est considrablement simplifie par le fait de la prise en compte de = 0 . On obtient aprs
intgration:
9782100544639.indb 199
199
09/02/10 14:58
40a
40a
160
24e 2 + 20a 2
6e 2
24e 2 + 20a 2
Ee
576
SYM
9782100544639.indb 200
6e 2
40a
6e 2
10a 2
6e 2
40a
24e 2 + 20a 2
20a
40a
160
20a
80
40a
6e 2
12e 2 + 5a 2
20a
40a
40
6e 2
24e 2 + 20a 2
20a
40
20a
20a
80
160
40a
40a
160
10a 2
6e 2
6e 2
20a
12e 2 + 5a 2 20a
6e 2
20a
40
40a
20a
40
20a
12e 2 + 10a 2
6e 2
6e 2
10a 2
20a
20a
12e 2 + 5a 2
6e 2
6e 2
12e 2 + 5a 2
(9.43)
20a
40a
2
2
10a
6e
6e 2
12e 2 + 10a 2
40a
40a
24e 2 + 20a 2
6e 2
2
2
24e + 20a
40a
9 lments de plaque
9.2 Exemple 9: plaque
simplement appuye sur 4 cts
200
09/02/10 14:58
9 lments de plaque
20a 2 x
20a 2 y
40a 3 y
160 w4
0
0
0
=
0
qa 2
A
(9.44)
10a 2
6e 2
6e 2
24e 2 + 20a 2
10a 2
24e 2 + 20a 2
6e 2
6e 2
Ee
6e 2
6e 2
24e 2 + 20a 2
10a 2
576
6e 2
10a 2
24e 2 + 20a 2
6e 2
40a
20a
20a
40a
6 a 3 (11e 2 + 5a 2 )
e2
1x
42 a 3
3
2 x
42 a
q
=
2y
2 + 25a 2
3
2
2
Ee
48
e
(
)
6 a (11e + 5a )
3y
2
e
w4
2
4
2 2
4
3a ( 72e + 110a e + 25a )
5e 2
(9.45)
soit
w4 = 0.357 102 m ; 1x = 3 y = 0.714 102 rad ; 2 x = 2 y = 0.399 105 rad .
Dunod La photocopie non autorise est un dlit.
Plaque rectangulaire
9782100544639.indb 201
201
09/02/10 14:59
9 lments de plaque
{ f }
1
xx
0
y
Ez
= yy =
1
0
2
y
1
1
xy
0 0
y x
2
+
y
x
(9.46)
{ f } = H f { f } = H f B f {q1}
e
avec z correspondant aux fibres infrieure et suprieure de la plaque, soit z = .
2
xx
M xx
12 z
M xx
M yy = H 1 f B f {q1}
M
xy
(9.47)
Pour lexemple de la plaque carre, les expressions obtenues pour les moments
sont:
M xx
e 2 ( 4 18 ) + 5a 2 ( 1)
2
qa
2
2
e ( 4 + 18 ) 5a ( + 1)
M yy =
2 + 25a 2
4
48
e
(
) 2e 2 ( 4 4 8) + 10a 2 ( 2 )
M
xy
soiten fonction des coordonnes x, h:
= 1; = 1
M xx 0.56 N.m
{M } = M yy = 0.56 N.m
M 999 N.m
xy
= 1; = 1
M xx 0.56 N.m
{M } = M yy = 1000 N.m
M 499 N.m
xy
= 0 ; = 0
M xx 500 N.m
{M } = M yy = 500 N.m
M 499 N.m
xy
202
9782100544639.indb 202
= 1; = 1
M xx 1000 N.m
{M } = M yy = 0.56 N.m
M 499 N.m
xy
= 1; = 1
M xx 1000 N.m
{M } = M yy = 1000 N.m
M
0 N.m 09/02/10
xy
14:59
= 0 ; = 0
M xx 500 N.m
{M } = M yy = 500 N.m
9 lments de plaque
9.2 Exemple 9: plaque
M 499 N.m
simplement appuye sur 4 cts
xy
= 1; = 1
= 1; = 1
M xx 1000 N.m
{M } = M yy = 0.56 N.m
M 499 N.m
xy
M xx 1000 N.m
{M } = M yy = 1000 N.m
M
0 N.m
xy
e
do les valeurs des contraintes pour z = + :
2
= 1; = 1
xx 0
{ } = yy = 15 MPa
7.5 MPa
xy
xx 0
{ } = yy = 0
15 MPa
xy
= 1; = 1
= 0 ; = 0
xx 7.5 MPa
{ } = yy = 7.5 MPa
7.5 MPa
xy
= 1; = 1
xx 15 MPa
{ } = yy = 0
7.5 MPa
xy
= 1; = 1
xx 15 MPa
{ } = yy = 15 MPa
0
xy
+y
yz
1 0 y
E
{ c } = =
= [ H c ] { c } = [ H c ] [ Bc ] {q1} (9.48)
0 1
2
1
w
(
) +
xz
x
x
On pose que leffort tranchant dans la direction x et y est gal lintgrale des
contraintes de cisaillement sur lpaisseur e . On a donc aprs correction des
contraintes de cisaillement par le facteur k:
e
Txz
2
{T } = = k e { c } dz = k e { c } = [ H 1c ] { c } = [ H 1c ] [ Bc ] {q1} (9.49)
T yz
2
9782100544639.indb 203
203
09/02/10 14:59
9 lments de plaque
avec [ H 1c ] =
Eek 1 0
.
2 (1 + ) 0 1
Les rsultats obtenus au chapitre prcdent ntant pas conformes avec la thorie des plaques, il est alors ncessaire daugmenter le maillage afin dapprocher au
mieux la solution objective. Pour ce faire, chacune des lignes dappui a t divise
en un nombre de segments pair.
Tableau 9.1 wA en fonction du maillage.
Dplacement au centre de la plaque (mm)
Maillage
wA
Objectif
Erreur
22
3.57
4.64
23.06%
44
4.52
4.64
2.59%
66
4.60
4.64
0.86%
88
4.62
4.64
0.43%
10 10
4.63
4.64
0.22%
12 12
4.64
4.64
0.00%
14 14
4.64
4.64
0.00%
Il apparat aprs calculs que les dplacements sont trs rapidement proches de la
valeur escompte et ce ds un maillage 4 4. Ceci tait nanmoins prvisible
puisque lquation gnrale des plaques (9.19) suppose une fonction de dplacement w du quatrime degr. Notre lment de plaque ayant une fonction dapproximation linaire en w , cinq points sont ncessaires pour la reconstituer soit
quatre lments.
Il est donc recommand, avec ce type dlment, davoir au minimum quatre divisions entre lignes dappui. Bien videmment, cette rgle, qui na rien dabsolu,
varie en fonction du degr de sophistication de llment choisi. En effet, nous
verrons au chapitre 10 quil sera possible dobtenir des rsultats similaires ceux du
maillage 4 4 avec un lment dit intgration rduite Q4R. Cet lment dont
les caractristiques sont similaires celles du Q4 de base se diffrencie par une intgration rduite un point des termes de rigidit de flexion, ceci ayant pour effet de
rduire le caractre trop rigide du Q4.
Pour ce qui concerne les moments M xx et M yy qui sont gaux dans ce cas, le
mme constat est valable pour le maillage. Ceci tant, il faut distinguer les rsultats
204
9782100544639.indb 204
09/02/10 14:59
9 lments de plaque
lisss et non lisss qui nvoluent pas de la mme faon. En effet et comme constat
au chapitre 8.3.5, on saperoit quil est prfrable de retenir les efforts non lisss
avec des maillages peu denses.
22
0
44
66
Maillage
8 8 10 10 12 12 14 14
400
Moments lisss
au centre de la plaque (N.m)
200
600
Mxx
800
Objecf Mxx
1 000
1 200
1 400
1 600
1 800
Maillage
22
0
44
66
8 8 10 10
12 12 14 14
200
400
600
Mxx
800
Objecf Mxx
1 000
1 200
1 400
1 600
9782100544639.indb 205
205
09/02/10 14:59
9 lments de plaque
Objectif
Erreur
22
8.23
10.53
21.84%
44
10.39
10.53
1.33%
66
10.44
10.53
0.85%
88
10.48
10.53
0.47%
10 10
10.50
10.53
0.28%
12 12
10.51
10.53
0.19%
14 14
10.52
10.53
0.09%
Les conclusions pour la plaque rectangulaire sont tout fait similaires celles tablies partir des rsultats du tableau 9.1, les dplacements se rapprochant de la
thorie ds un maillage quatre divisions. Concernant les moments, la progression est toujours lgrement plus lente puisque les rsultats sont conformes ds six
mailles pour les valeurs lisses et quatre pour celles non lisses. On notera encore
une fois une meilleure convergence des rsultats non lisss.
206
9782100544639.indb 206
09/02/10 14:59
9 lments de plaque
22
0
44
66
Maillage
8 8 10 10
12 12 14 14
500
Moments lisss
au centre de la plaque (N.m)
1 000
1 500
Myy
Objecf Myy
Mxx
Objecf Mxx
2 000
2 500
3 000
3 500
4 000
4 500
5 000
9782100544639.indb 207
207
09/02/10 14:59
9 lments de plaque
Maillage
22
0
44
66
8 8 10 10
12 12 14 14
500
1 000
Myy
Objecf Myy
Mxx
Objecf Mxx
1 500
2 000
2 500
3 000
3 500
4 000
4 500
a
a/2
a
A
Appuis simples
Appuis simples
2a
Appuis simples
Figure 9.19 Exemple 10: plancher dalle.
Application numrique:
a = 2 m, b = 4 m.
Module longitudinal E = 1.962 1010 N/m2.
Coefficient de Poisson n = 0.1.
paisseur de la dalle e = 0.2 m.
Charge rpartie transversale q = 25000 N/m2.
Conditions dappui: appuis simples aux extrmits + 4 appuis simples (poteaux).
208
9782100544639.indb 208
09/02/10 14:59
9 lments de plaque
Sachant que la taille retenue pour chacun des lments est de 0.5 m, les rsultats
obtenus sont:
9782100544639.indb 209
209
09/02/10 14:59
9 lments de plaque
Avec llment Q4 dcrit au chapitre 9.2.3 et des mailles de 0.5 mtre de ct, les
rsultats obtenus montrent une lgre divergence entre les modles surfaciques et
poutres, le dplacement vertical en B tant de 2.21 mm pour 2.11 attendu. Ceci
tant, il faut relativiser ce constat dans la mesure o la flexion transversale nest pas
prise en compte par la modlisation filaire. Enfin, il y a galement un effet minime
d linfluence du coefficient de Poisson.
En intgrant le moment sur une coupe correspondant la demi-largeur de la dalle,
on trouve respectivement 43685 N.m et 44023 N.m pour les moments rsultants
sur appui et en milieu de trave. Compare au modle poutres, la rsultante en
trave reste conforme aux attentes. Par contre, celle sur appui diverge sensiblement
de la valeur objective. Ceci peut sexpliquer par le caractre linaire de llment
Q4 utilis qui engendre une progression galement linaire du moment sur son
domaine. De ce fait et en raison de la forte variation sur appui, celui-ci est incapable de reconstituer correctement cette courbe. Deux options soffrent alors
nous, soit augmenter la densit de maillage dans la direction longitudinale, soit
utiliser un lment plus performant.
Tableau 9.3 Moments rsultants obtenus avec diffrentes tailles dlments.
Dplacement vertical
en B (mm)
Moment en trave
(N.m)
Q4 0.50
2.21
44023
43685
Q4 0.25
2.28
43975
49789
Q4 0.17
2.29
43996
51846
Q9 0.50
2.35
45293
55433
210
9782100544639.indb 210
09/02/10 14:59
9 lments de plaque
9782100544639.indb 211
211
09/02/10 14:59
9 lments de plaque
Largeur a
0
0,25
0,5
0,75
1,25
1,5
1,75
2 500
5 000
7 500
Q9 - 0,5
10 000
Q4 - 0,50
12 500
Q4 - 0,25
15 000
Q4 - 0,17
17 500
20 000
22 500
25 000
45 000
40 000
Q9 - 0,50
35 000
Q4 - 0,50
30 000
Q4 - 0,25
25 000
Q4 - 0,17
20 000
15 000
10 000
5 000
0
0
0,25
0,5
0,75
1,25
1,5
1,75
Largeur a
Figure 9.26 Coupe sur appui.
212
9782100544639.indb 212
09/02/10 14:59
9 lments de plaque
9782100544639.indb 213
213
09/02/10 14:59
9782100544639.indb 214
09/02/10 14:59
10 lments de coque
A
10.1 Aspects thoriques
Une coque est un lment surfacique tridimensionnel prenant en compte simultanment les effets de membrane et de flexion. Les surfaces modliser peuvent
tre planes ou ventuellement gauches. On distinguera ainsi deux grandes familles
dlments relatives aux coques; les coques planes ou courbes qui sont elles-mmes
divises en deux sous catgories: les coques peu profondes (shallow shell) ou profondes (deep shell). Du fait de son caractre tridimensionnel, llment coque possde six degrs de libert par nud.
Pour les structures formes de surfaces planes telles que les btiments constitus de
voiles en bton arm ou les ponts caissons, le recours un lment de coque plane
se justifie amplement. Ceci tant, la mme dmarche est gnralement adopte
pour les surfaces gauches car les lments de coques courbes qui sont difficiles
formuler, sont dune utilisation dlicate.
La surface gauche est donc, dans ce cas, reconstitue par un assemblage de surfaces
planes. Cette approche prsente nanmoins plusieurs inconvnients. Le premier
concerne la prcision de cette reconstitution qui suppose des lments de petites
tailles et donc un nombre de facettes important. Les charges appliques perpendiculairement sur ces facettes ne sont pas non plus parfaitement reprsentatives de la
ralit du fait de la non prise en compte de la courbure dans le calcul des vecteurs
charges.
Comme pour les arcs maills avec des lments de poutres droites, on constate des
pousses au vide lies aux efforts de membrane. Bien videmment, ces pousses
nexistent pas dans la ralit puisque ces efforts sont toujours tangents la courbure
de la coque.
Le principe gnralement le plus rpandu pour la formulation dun lment de
coque plane consiste considrer la superposition dun lment de membrane
celui dune plaque de mme gomtrie quadrangulaire ou triangulaire. Toutefois,
cette approche suppose que ces deux effets sont dcoupls.
Llment membrane amenant deux degrs de libert en translation u et v , la
plaque un seul en translation w et deux en rotation x et y , le sixime associ
la normale au plan de la coque z est donc absent ce qui pose un rel problme
9782100544639.indb 215
215
09/02/10 14:59
10 lments de coque
pour les lments de coques coplanaires puisque les termes de rigidit associs sont
nuls1. Il en rsulte une singularit de la matrice de rigidit.
w x
0 0
0 0
0
0
0 0
0 0
0 0
0
0
u v
Membrane
Flexion
u
v
w
x
y
z
Pour rsoudre ce problme, une solution propose par [8] consiste associer des
termes de rigidit fictifs faibles la rotation z assurant ainsi la stabilit numrique du systme sans en perturber la solution. On pose alors que lnergie de
dformation associe llment est gale :
We = We + We z
(10.1)
z = N 1 z1 + N 2 z 2 + N 3 z 3 + N 4 z 4
(10.3)
1
(1 + i ) (1 + i ) .
4
Les dformations associes tant directement dpendantes des variations de z
par rapport x et y, on peut crire pour un lment n nuds daprs (7.36) que:
avec N i =
N i
N n
z
N1
z
.
.
1
1
{ ze } (10.4)
=[J]
= [ J ]
N
N
N
i
n
z
1
.
.
y
B z
k z = B z H z B z dSe
(10.5)
Se
1. Certains logiciels utilisent cependant des lments de coque 5 DDL par nud.
216
9782100544639.indb 216
09/02/10 14:59
10 lments de coque
Ee 3 1 0
avec H z =
.
12 0 1
Ee 3
provient du fait que cette dmarche sinspire de celle rete12
nue pour le calcul de la part de rigidit associe la flexion. Le coefficient qui
varie en fonction de la prcision de lordinateur utilis, doit tre choisi suffisamment petit1 pour ne pas affecter les rsultats.
Pour tablir la matrice de rigidit dun lment de coque plane quatre nuds, il
suffira donc de reprendre lexpression (8.12) obtenue pour la membrane quadrangulaire au chapitre 8.1 et de ladditionner celle de la plaque (9.40) du chapitre
9.2.3.2.
On notera cependant que dans ce cas la formulation retenue pour la partie flexion,
est celle de la plaque paisse. Cet lment de coque paisse reprend donc intgralement les dveloppements relatifs au blocage en cisaillement (cf. 9.2.3) et
la prise en compte du cisaillement transverse. Il serait galement tout fait possible
de dvelopper un lment de coque mince se basant la thorie de Kirchhoff.
Le choix du terme
u1 v1 w1 x1y1 z1 u2 v2 w2 x2 y2 z2 u3 v3 w3 x3 y3 z3 u4 v4 w4 x4 y4 z4
u1
v1
w1
x 1
y 1
z 1
u2
v2
w2
x 2
y 2
z 2
u3
v3
w3
x 3
y 3
z 3
u4
v4
w4
x 4
y 4
z 4
1. De lordre de 103-104.
9782100544639.indb 217
217
09/02/10 14:59
10 lments de coque
Appuis simples
2a
qh
qh
Appuis simples
Appuis simples
b
x
Figure 10.3 Exemple 11: coque comprime flchie.
1 y = 3 x = 4 x = 4 y = w1 = w2 = w3 = 0
De plus et pour simplifier le problme, on considrera uniquement le dplacement
u4 non nul. Il reste donc six inconnues.
y
1
218
9782100544639.indb 218
09/02/10 14:59
10 lments de coque
Considrant les mmes hypothses de calcul que celles utilises au chapitre 9.2.1
(E = 2.1 1011 N/m2, n = 0.0, q = 1 104 N/m2, a = b = 1 m, e = 0.02 m), on
ajoute celles-ci une charge uniformment rpartie qh = 1 107 N / m applique
dans lpaisseur de la coque.
tant donn quil sagit dune simple superposition des effets de membrane et de
flexion, on trouve bien videmment pour la partie flexion les mmes rsultats que
ceux rsultants de (9.45) soit:
Pour ce qui concerne les effets de membrane, le dplacement u3 peut tre directement dduit des rsultats de lexemple 5 pour = 0 (cf. 8.1.5) en posant que les
dplacements du nud 1 dans (8.23) sont nuls. On a donc:
q a
u3 = h = 0.00238 m .
E e
Ceci tant, lutilisation dun lment de coque avec intgration rduite (Q4R) pour
les parties membrane et flexion amnera une amlioration sensible des rsultats
pour le dplacement w et les rotations:
w4 = 0.447 102 m ; 1x = 3 y = 0.893 102 rad ; 2 x = 2 y = 0.286 105 rad
mais galement une augmentation de ceux associs aux effets de membrane
( u3 = 0.00317 m ) .
Bien videmment, laugmentation du maillage amnera encore une fois la superposition des rsultats obtenus en 9.2.4 ceux des effets de membrane proprement
dits.
9782100544639.indb 219
219
09/02/10 14:59
10 lments de coque
220
9782100544639.indb 220
09/02/10 14:59
10 lments de coque
L/2
L
Figure 10.8 Exemple 12: tude dune poutre en I.
Sachant quil est admis de prendre la section deffort tranchant gale la section de
lme (30 1.15 cm), la dforme au centre de la poutre sera donne par:
L
v = vM
2
flexion
tranchant
3
PL
PL
+ vT =
+
48EI z 4GS1, y
(10.6)
9782100544639.indb 221
221
09/02/10 14:59
10 lments de coque
vT = xy
P
Ty L
L
L
PL
=
= 2
=
2 GS1, y 2 GS1, y 2 4GS1, y
(10.7)
PL3
PL
L
+
= 1.838 + 0.538 = 2.376 mm .
do v =
2 48EI z 4GS1, y
Leffet de leffort tranchant nest donc pas ngligeable. Par ailleurs et daprs la
thorie des poutres, les contraintes normales dans les semelles devraient tre gales
en milieu de trave :
PL
M
(10.8)
max/min = z = 4 = 82.35MPa
Wz
Wz
10.3.2 Modlisation surfacique
h-t
Plans moyens
des lments
coques
b
Figure 10.9 Plans moyens.
Llment retenu pour le calcul est un lment de coque paisse carr quatre
nuds de ct 5 cm et de type Q4 tel que celui dcrit au chapitre 10.1. Le modle
comporte un peu plus dun millier de nuds et autant dlments.
Pour se conformer au mieux avec la gomtrie filaire, les conditions dappui sont
introduites au niveau de la fibre moyenne de la poutre. Le nud correspondant
lappui gauche est donc bloqu dans les trois directions de lespace alors que celui
de lappui droit est laiss libre horizontalement suivant x.
1. h-tf au lieu de h-2tf.
222
9782100544639.indb 222
09/02/10 14:59
10 lments de coque
Une premire analyse de la dforme de la figure 10.10 met en vidence un dplacement lgrement suprieur (2.5 mm) celui obtenu avec le modle filaire (2.38
mm). Il convient nanmoins de relativiser ce rsultat et donc de ne pas perdre de
vue leffet local d la charge ponctuelle mais galement linfluence du coefficient
de Poisson. En effet et en recalculant le modle avec n gal zro, on obtient un
dplacement en milieu de poutre de 2.37 mm (cf. figure 10.12).
Une fois les dplacements valids, il est alors possible de passer la comparaison des
contraintes normales. Il faudra nanmoins prendre garde bien exprimer celles-ci
au niveau des fibres extrmes des semelles. Les fibres infrieure et suprieure des
lments coques et plaques tant repres par rapport leurs normales sortantes, il
faudra bien vrifier que le postraitement seffectue au niveau de la bonne fibre et ce
tout particulirement pour les coques paisses. Pour ce faire, on pourra afficher les
repres locaux des lments coques ou plaques et ainsi vrifier leurs orientations.
Dans notre cas dtude, les contraintes normales au niveau des semelles sont gales
89 MPa (cf. figure 10.11), valeurs lgrement suprieures celles escomptes.
Il faut encore une fois relativiser ce rsultat dans la mesure o leffet du coefficient
de Poisson nest pas pris en compte dans la rfrence. Pour sen assurer, il suffira
de recalculer le modle avec une valeur de n nulle. Dans notre exemple, cet effet
reprsente une majoration des contraintes de lordre de 4 MPa1.
1. La contrainte obtenue dans ce cas est de 85.35 MPa (cf. figure 10.13).
9782100544639.indb 223
223
09/02/10 14:59
10 lments de coque
224
9782100544639.indb 224
09/02/10 14:59
10 lments de coque
9782100544639.indb 225
225
09/02/10 14:59
9782100544639.indb 226
09/02/10 14:59
A
11.1 Aspects thoriques
11.1.1 Notion de grandes dformations
u
dx
x
dx
v+
v
x
v
dx
x
dx
v
dx
x
v
x
u v
dx
x x
x
Figure 11.1 Calcul de exx.
9782100544639.indb 227
227
09/02/10 14:59
u v
u v
dx = dx + dx + dx + dx
x x
x x
soit dx = dx 1 + 2
(11.1)
2
u u v
v u u v
+ + + 2 + .
x
x
x
x
x x x
0
xx =
2
2
u 1 u v
+ + = xx 0 + xxL
x 2 x x
(11.2)
(11.3)
u
dformation correspondant lhypothse des petites dformations.
x
On pourrait galement dmontrer (cf. [9]) que:
avec xx 0 =
yy
xy
2
2
v 1 u v
=
+ +
y 2 y y
v u u u v v
+
=
+
+
x y x y x y
(11.4)
9782100544639.indb 228
0
N1
N2
0
u1 x
0
0 v1 1 L
=
N 2 u2
x
0
1
v2
L
x
L
0
u
1
0 v
1
x u2
L v2
(11.5)
09/02/10 14:59
v2
u2
u1
v1
Les variations de u et v par rapport x pourront alors tre relies aux dplacements nodaux par la matrice [G ] en posant que:
du dN1
dx dx
=
dv 0
dx
1
L
=
0
dN 2
dx
0
dN1
dx
0
1
L
0
1
L
0
u
1
0 v
1
dN 2 u2
dx v2
u
1
0 v
1
= [G ] {qe }
1 u2
L v2
(11.6)
La matrice [G ] tant associe aux drives de 1er ordre des fonctions de forme,
celle-ci aura pour expression dans le cas plan:
u
x x
u
y y
=
v 0
x
v
0
y
0
x
0
u
y
=
v
0
0
y
[ N ] {qe } = [G ] {qe }
(11.7)
Le terme xxL issu de (11.3) qui est li aux grandes dformations, peut alors
scrire:
9782100544639.indb 229
229
09/02/10 14:59
xxL
u
v x
x v
x
2
2
1 u v 1 u
= + =
2 x x 2 x
xxL =
(11.8)
1
1
{qe }T [G ]T [G ] {qe } = BL ( q ) {qe }
2
2
xx = xx 0 + xxL = [ B0 ] + BL ( q ) {qe }
2
(11.9)
1
1
0
0 {qe } correspondant aux petites dforavec xx 0 = [ B0 ] {qe } =
L
L
{ } = { 0 } + { L } = [ B0 ] + BL ( q ) {qe }
(11.10)
avec [ B0 ] et BL ( q ) correspondant respectivement aux petites et grandes dformations/grands dplacements. On notera cependant que la matrice BL ( q ) est
elle-mme fonction des dplacements.
Le problme tant non linaire, il nest pas possible dutiliser directement lexpression de lnergie potentielle (4.5). Il est donc ncessaire, pour caractriser lquilibre, dutiliser sa forme dite variationnelle soit:
dE = dW dT
=
T
T
{d } { } dV {du} { f v } dV {du} { f s } dS = 0
Ve
Ve
(11.11)
Se
d
Figure 11.3
fonction de
230
9782100544639.indb 230
09/02/10 14:59
{d } = [ B0 ] + BL ( q ) {dqe } = [ B ] {dqe }
do lexpression de lquilibre lmentaire:
(11.12)
T
T
{dqe } [ B ] { } dVe {dqe } ({ f ev } + { f es }) = 0
(11.13)
Ve
{dqe }T
T
[ B ] { } dVe { f ext } = 0
Ve
{ f ext } , on
(11.14)
Soit:
T
[ B ] { } dVe { f ext } = 0
(11.15)
Ve
Le terme
[ B ] { } dVe
Ve
la diffrence
[ B ] { } dVe { f ext }
{ fi } ,
Ve
{re ( q )} dans llment. Bien videmment et en vertu de (11.15), ces forces rsi-
(11.16)
Ve
Lors dun faible accroissement des dplacements, on devra malgr tout encore vrifier {re ( q + dq )} 0 . Comme {re ( q + dq )} peut tre exprim sous la forme:
On a de mme:
(11.17)
(11.18)
est gal au
(11.19)
Ve
Ve
9782100544639.indb 231
Ve
09/02/10 14:59
(11.21)
T
avec:
(11.22)
T
T
[ B ] [ H ] [ B ] dVe = [ B0 ] [ H ] [ B0 ] dVe +
Ve
Ve
(11.23)
[ke 0 ]
T
T
T
[ B0 ] [ H ] [ BL ] dVe + [ BL ] [ H ] [ B0 ] dVe + [ BL ] [ H ] [ BL ] dVe
Ve
Ve
Ve
[keL ]
T
T
[dBL ] { } dVe = [G ] [G ] {dqe } { } dVe = [ke ] {dqe }
Ve
(11.24)
Ve
[ke ] = [G ] [ ] [G ] dVe
(11.25)
Ve
De part sa formulation, la matrice de rigidit lmentaire [ ke0 ] est associe lhypothse des petites dformations (tablie au 4.2.2) alors que [ keL ] qui est fonction des dplacements nodaux, correspond celle des grandes dformations/grands
dplacements. [ ke ] est appele matrice des contraintes initiales et correspond
ltat de rigidification dune structure due son tat de contrainte.
Cet effet de rigidification couple la gomtrie dune structure ltat de ses dplacements transversaux. Elle ncessite dtre prise en compte tout particulirement
dans le cas des structures fines et lances. [ ke ] vaut par exemple pour un lment
barre:
[ke ] = [G ] [ ] [G ] dVe = [G ]T
Ve
N
S
[G ] dS dx (11.26)
=S
1 0 1 0
N 0 1 0 1
(11.27)
[ke ] = 1 0 1 0
L
0 1 0 1
Les rgles dassemblage tant identiques celles tablies au chapitre 4, il est possible de gnraliser les dveloppements prcdents lensemble de la structure en
posant:
[ K T ] {dq} = {F } ([ K 0 ] + [ K L ] + [ K ]) {dq} = {F }
(11.28)
232
9782100544639.indb 232
09/02/10 14:59
avec:
grands dplacements.
[ K ] : Matrice de rigidit globale associe aux contraintes initiales.
[ K T ] {dq} = {F } ([ K 0 ] + [ K ]) {dq} = {F }
11.1.3 Instabilit
(11.29)
{Fcritique } = {F }
(11.30)
9782100544639.indb 233
233
09/02/10 14:59
Nous verrons au chapitre 11.3 quil est galement possible dvaluer le comportement des structures dans le domaine post critique soit au-del de la charge critique.
La grande difficult sera alors de pouvoir saffranchir des singularits de la matrice
tangente [ K T ] .
11.1.4 Les mthodes de calcul non linaire
Incrment
n itraons
5 itraons
3 itraons
F1
234
9782100544639.indb 234
09/02/10 14:59
Le principe consistera mettre en uvre un processus itratif permettant dquilibrer les charges statiques et ainsi de trouver une suite dincrments de dplacements satisfaisant la relation dquilibre:
[ K T i ] {qi } = {F } {Fi }
[ K T ] {dq} = {R }
(11.32)
avec {qi } = {qi +1} {qi } et {Fi } : efforts internes calculs litration i.
Dans labsolu, larrt du processus itratif seffectue lorsque le reste litration i
{Ri } cest--dire la diffrence entre {F } et {Fi } tend vers 0. On parlera alors de
convergence en forces. Il est possible galement denvisager une convergence en
dplacements en posant que la diffrence {qi } = {qi +1} {qi } tend vers 0 ou
celle en nergie qui correspond en fait un mixe des deux prcdentes. Ces trois
approches pourront dailleurs tre utilises de manire simultane.
F
Fi+2
Fi+1
Fi
q
qi
qi+1
qi+2
9782100544639.indb 235
235
09/02/10 14:59
La rigidit tangente [ K T ] tant dfinie partir des rigidits induites par les petits
dplacements [ K 0 ] , les grandes dformations ou grands dplacements [ K L ] et les
contraintes initiales [ K ] , les matrices [ K L ] et [ K ] restent nanmoins dpendantes de ltat de dformation de la structure. Elles doivent donc tre mises jour
suivant une frquence dpendant de la mthode de Newton Raphson utilise.
Dans le cas de lutilisation dune mthode de Newton Raphson complte, ces deux
matrices seront actualises chaque itration de convergence. Si cette remise
jour induit des temps de calcul trop importants par rapport au nombre ditrations ncessites par la convergence, il est possible dutiliser la mthode de Newton
Raphson modifie qui consiste ne redfinir ces matrices quau dbut de chaque
incrment. Ceci revient utiliser la mme rigidit tangente sur lincrment de
charges considr.
Fi+2
Fi+1
Fi
q
qi qi+1 qi+2
Figure 11.6 Mthode Newton Raphson modifie.
9782100544639.indb 236
09/02/10 15:00
peuvent ventuellement tre utilises pour optimiser les temps de calcul. De plus
et en prsence de courbes non linaires prsentant des inflexions, asymptotes ou
extremums, il est parfois ncessaire lorsque lon ne dispose pas de la mthode de
longueur darc de piloter la monte en charge en dplacement impos. Cependant, ce type de solution nest applicable que dans des cas trs spcifiques.
Mthode de longueur darc
{qi }T {qi } + i 2 = l 2
{ di } = [ K T i ]1 {Ri }
(11.34)
{ Lqi } = [ K T i ]1 {F }
(11.35)
{ qi } = { di } + i { Lqi }
(11.36)
Enfin, la suite du processus consiste couper la courbe avec une surface de forme
variable. Riks [14] et Wempner [15] proposent ds les annes 70 une projection
perpendiculaire la matrice tangente. Cest celle que nous retiendrons dans les
exemples qui suivent.
Plus rcemment, Crisfield [11] [12] a dvelopp une mthode base sur un arc
sphrique qui prend en compte uniquement le terme D qi dans le calcul de la longueur darc soit:
{qi }T {qi } = l 2
(11.37)
Au dbut de litration i, lquilibre des forces est caractris par la diffrence entre
charges extrieures pondres par le facteur i et les efforts internes {Fi } . Ce
calcul revient valuer le reste:
{Ri } = i {F } {Fi }
(11.38)
qui doit dans labsolu tendre vers 0. Si tel nest pas le cas, il est alors ncessaire de
dterminer un nouvel incrment de charge D i . Reprenant la mthode du plan
9782100544639.indb 237
237
09/02/10 15:00
normal (cf. [16]), on pose tout dabord que les vecteurs t0 et ni sont respective qi
q0
ment gaux
. Ces deux vecteurs ont de fait un produit scalaire
et
i
0
nul soit:
{q0 }T { qi } + 0 i
=0
(11.39)
d1
1
Plan de coupure
t0
0
n1
n0
q
q1
q2
q1
q0
i =
{q0 } { d i }
{q0 }T { Lqi } + 0
(11.40)
Les diffrentes variables sont ensuite mises jour pour le calcul de litration suivante:
i +1 = i + i
(11.41)
Cependant et pour dmarrer le processus, lutilisateur devra estimer la valeur initiale D 0 partir dun lincrment de longueur darc. Sachant que la variation de
dplacement litration 0 peut tre estime par:
238
9782100544639.indb 238
09/02/10 15:00
{q0 } = 0 { Lq0 } = 0 [ K T 0 ]1 {F }
(11.42)
l
T
1 + { Lq0 } { Lq0 }
(11.43)
(11.44)
De plus et en vertu de (11.43), il sera ncessaire de prciser le signe de D 0 . Gnralement, trois critres [17] sont utiliss pour dterminer la direction du processus:
Le signe du dterminant de la matrice tangente [ K T ] .
T
{Ri } {F }
{qi } {q0 }
(11.45)
Ces critres ntant pas exhaustifs, ceux-ci peuvent varier suivant la nature du code
de calcul utilis.
9782100544639.indb 239
239
09/02/10 15:00
x
3
L0
L0
X
Figure 11.8 Exemple 13: shed symtrique.
P = F/2
V2
h
Lj
L0
Figure 11.9 Exemple 13: modlisation barre.
240
9782100544639.indb 240
09/02/10 15:00
nanmoins de bloquer horizontalement le nud 2 et de ne considrer que la moiti de la charge applique. Deux cas de figure sont ensuite envisageren fonction
de la nature de sa hauteur h . Si h est significatif par rapport la base L0 du
shed, la longueur L du barreau sera directement fonction du dplacement V2 . Si
au contraire h est suppos petit par rapport L0 alors L pourra tre considr
comme invariant. On aura alors quel que soit lincrment de charge j: Lj L.
Considrant lhypothse h petit dans un souci de simplification, les dimensions et caractristiques du shed prises en compte dans les calculs seront donc: L0
= 2000 mm, h = 20 mm, S = 200 mm2, E = 210000 N/mm2, P 8.5 N. La
longueur initiale L du barreau est ainsi gale 2000.1 mm. Les dplacements du
nud 1 tant nuls, la dformation dans le barreau est gale daprs (11.3):
xx =
2
2
2
2
u 1 u v u2 1 u2 v2
+ +
+ + =
x 2 x x L 2 L L
(11.46)
v2
V2
u2
2
Figure 11.10 Dcomposition de V2
9782100544639.indb 241
xx =
u2 1 v2 V22 + 2h V2
+ =
L 2 L
2 L2
(11.48)
241
09/02/10 15:00
N
= E xx do
S
(11.49)
( )
N sin j P = 0
soit:
(11.50)
h + V2
PL
P =0 N =
L
h + V2
(11.51)
ES
(V23 + 3h V22 + 2h 2 V2 )
2 L3
(11.52)
9
8
7
6
5
4
3
2
1
0
1
2
3
4
5
6
7
8
9
P(N)
D
E
45
40
35
30
25
20
V2 (mm)
15
10
La courbe non linaire P = f (V2 ) de la figure11.11 met en vidence deux extrmums en A et C ainsi quun point de changement de courbure en B qui correspond
la position horizontale du barreau. Ensuite, le barreau continue sa course pour
atteindre le point D symtrique de celui dorigine par rapport Y. Enfin, on notera
les valeurs de charge nulle associes B et D.
242
9782100544639.indb 242
09/02/10 15:00
[k0 ] = [ B0 ]
0 1
0 0
0 1
0 0
0
0
0
(11.53)
1 0 1 0
N 0 1 0 1
[ ] [G ] dVe =
L 1 0 1 0
0 1 0 1
(11.54)
Ve
ES 0
[ H ] [ B0 ] dVe =
L 1
[k ] = [G ]
Ve
[k0 ]
v2
L2
u2
L2
v2
L2
(11.55)
{ } = { 0 } + { L } = [ B0 ] + BL ( q ) {q} .
[kL ] = [ B0 ]T [ H ] [ BL ] dVe
Ve
T
T
[ BL ] [ H ] [ B0 ] dVe + [ BL ] [ H ] [ BL ] dVe
Ve
9782100544639.indb 243
(11.56)
Ve
243
09/02/10 15:00
Soit
u2 ( 2L + u2 ) v2 ( L + u2 ) u2 ( 2 L + u2 ) v2 ( L + u2 )
v2 ( L + u2 )
v22
v22
ES v2 ( L + u2 )
(11.57)
[kL ] = 3
L u2 ( 2 L + u2 ) v2 ( L + u2 ) u2 ( 2 L + u2 ) v2 ( L + u2 )
vv22
v2 ( L + u2 )
v22
v2 ( L + u2 )
La rigidit tangente exprime dans le repre local tant gale la somme des trois
matrices (11.53), (11.54) et (11.57), on obtient pour finir:
NL2
2
2
+
+ ( L + u2 ) v2 ( L + u2 )
+
+
L
u
v
L
u
(
)
(
)
2
2
2
ES
ES
NL
NL2
2
2
v2 ( L + u2 )
+ v2
v2 ( L + u2 )
+ v2
ES
ES
ES
3
2
2
L
NL
NL
2
2
+ ( L + u2 ) v2 ( L + u2 )
+ ( L + u2 )
v2 ( L + u2 )
ES
ES
2
2
NL
NL
v2 ( L + u2 )
+ v22
+ v22
v2 ( L + u2 )
ES
ES
(11.58)
La dernire tape concerne lexpression de cette matrice en repre global qui
consiste, comme en 5.2, effectuer le changement de base suivant:
T
[ K T ] = R ( ) [ kT ] R ( )
cos
sin
avec R ( ) =
0
sin
cos
0
0
0
0
cos
sin
(11.59)
0
0
.
sin
cos
K T ( N ,V2 ) =
N ES
2
+ 3 ( L sin + V2 )
L L
(11.60)
h
, on obtient pour finir:
L
K T ( N ,V2 ) =
N ES
2
+ 3 ( h + V2 )
L L
(11.61)
244
9782100544639.indb 244
09/02/10 15:00
h + V2
correspond au sinus de langle j (cf. figure11.9) obtenu
L
lors du dplacement lincrment de charge j. Ceci permet de retrouver lexpression tablie par Bathe (cf. [18]) :
En fait, le terme
K T ( N ,V2 ) =
N ES 2
+
sin j
L
L
( )
(11.62)
qui dcoule dun changement de base effectu uniquement sur les matrices [ k0 ] et
[k ] . Cette approche revient en fait calculer la structure dans son tat dform
de manire successive en posant que:
[ K T ] = R ( j ) ([k0 ] + [k ]) R ( j )
(11.63)
Le chapitre 11.4.1 a permis de prsenter deux types de mthodes de NewtonRaphson. Celle consistant rvaluer la rigidit tangente chaque itration de
convergence, dite complte, est trs souvent la plus efficace mais galement la plus
coteuse en termes de temps de calcul. Les fortes variations de la dforme de
notre cas dtude justifient nanmoins lutilisation de cette mthode qui assure une
convergence rapide du processus. En application de (11.61), la rigidit tangente
fonction de leffort normal N i et du dplacement V2i , associs litration i, sera
donc remise jour systmatiquement et vaudra cette mme itration:
K Ti ( N i ,V2i ) =
N i ES
2
+ 3 ( h + V2 i )
L
L
V 2 + 2 h V2 i
avec daprs (11.49) N i (V2i ) = ES xx = ES 2i
2 L2
(11.64)
P valant toujours 8.5 N, on dcide deffectuer une monte en charge en 100 incrments de longueur gale et de tolrer un maximum de 50 itrations par incrment.
P
. La convergence
Au jme incrment, la force applique vaudra donc: F j = j
100
retenue est celle en forces. Le reste litration i aura pour expression daprs
(11.50):
h + V2 i
Ri = F j N i
(11.65)
L
Le programme Mathcad de calcul et lorganigramme correspondant sont indiqus pages suivantes. On notera nanmoins quil serait assez simple de transformer
celui-ci en mthode de Newton-Raphson modifie en transfrant le calcul de la
rigidit tangente dans la boucle j et en utilisant celle-ci chacune des itrations de
convergence.
9782100544639.indb 245
245
09/02/10 15:00
Variables:
Nsubstp : nombre dincrments (100)
Niter : nombre ditrations maxi par incrment (50)
e
: tolrance (0.001 par exemple)
: dplacement initial (0 en gnral)
v0
: effort normal initial (0 en gnral)
N0
Conv : convergence (0 en forces, 1 en dplacements)
Type : Newton Raphson 0 Complet, 1 Modifi
Newt onRaphson (NSubst p , NIt er, , v0 , N0, Conv , Ty pe) :=
V v0
0
N N0
0
for j 1 .. NSubst p
CV "Diverge"
j
j P
F
j
NSubst p
i0
q V
i
j 1
( i)
N1 EFFN q
i
K KTR N1 , q
j
if Ty pe
while i NIt er
K K if Ty pe
1
h + qi
R F EFFN q
i
j
i
L
if Conv 0 R < F j
i
( )
V q
j
N EFFN V
( j)
CV "Converge"
j
break
K KTR EFFN q , q
i
i+ 1
q i q
if Conv
( i) i)
+q
Ki
i+ 1
if Ty pe
1 q i < q 0
Vj q
i+ 1
( )
N j EFFN Vj
CVj "Converge"
break
return ( CV V F N ) if i
ii+ 1
NIt er
( CV V F N )
246
9782100544639.indb 246
09/02/10 15:00
Mthode de Newton
Raphson Complte
j =1
Fj = j
P
100
i = 0 , V2 i = V2 j 1
V 2 + 2hV
N i (V2i ) = ES xx = ES 2i 2 2i
2L
Ri = F j N i
h + V2i
L
Ri F j
j = j +1
V2j = V2i
oui
non
N i ES
2
+ 3 (h + V2i )
L
L
KTi ( N i ,V2i ) =
V2i +1 =
Ri
+ V2i
K Ti
i = i +1
oui
oui
i 50
non
FIN
j 100
non
FIN
Organigramme 11-1
9782100544639.indb 247
247
09/02/10 15:00
9
8
7
6
5
4
3
2
1
0
1
2
3
4
5
6
7
8
9
P(N)
Rfrence
Newton -Raphson (0.001)
45
40
35
30
25
20
15
10
V2 (mm)
Figure 11.12 Exemple 13: mthode de Newton Raphson avec e = 0.001.
Avec une tolrance de 0.001, le processus sarrte partir du 96e incrment soit
approximativement peu avant le point A (figure 11.12). ce niveau, la rigidit
tangente devient peu peu horizontale, rendant impossible une intersection avec
la suite de la courbe. Au mieux et comme indiqu la figure11.13, le processus
continue aprs une projection au niveau du point E.
P(N)
9
8
7
6
5
4
3
2
1
0
1
2
3
4
5
6
7
8
9
Rfrence
Newton-Raphson (0,01)
45
40
35
30
25
20
V2 (mm)
15
10
248
9782100544639.indb 248
09/02/10 15:00
1re trajectoire
possible
2e trajectoire
possible
q
q
Figure 11.14 Snap back / Bifurcation.
P(N)
9
8
7
6
5
4
3
2
1
0
1
2
3
4
5
6
7
8
9
Rfrence
Longueur d'arc
45
40
35
30
25
20
15
10
V2 (mm)
Figure 11.15 Exemple 13: mthode de longueur darc.
9782100544639.indb 249
249
09/02/10 15:00
Variables:
Nsubstp : nombre dincrments
Niter : nombre ditrations maxi par incrment
e
: tolrance (0.001 par exemple)
: dplacement initial (0 en gnral)
v0
: effort normal initial (0 en gnral)
N0
Conv : convergence (0 en forces, 1 en dplacements)
lmin : incrment de longueur darc mini (= incrment de dpart)
lmax : incrment de longueur darc maxi
lmax : longueur darc maxi
250
9782100544639.indb 250
09/02/10 15:00
9782100544639.indb 251
251
09/02/10 15:00
Fj = j P
oui
V2 i = V2 j 1
K Ti (N i , V2i ) =
non
i=0
Ri = ( j + ) P N i
N i ES
2
+ 3 (h + V2i )
L
L
h + V2i
L
oui
Ri P
j +1 = j +
V2 j = V2 i
non
P
LV2i =
K Ti
0 =
V 2 + 2hV
N i (V2i ) = ES xx = ES 2i 2 2i
2L
l
1 + LV2i
K Ti (N i , V2i ) =
LV2i =
V2i = 0 LV2i
V2i +1 = V2i + V2i
i =
= 0
dV2i =
Ri
K Ti
dV20 dV2i
dV20 LV2i + 0
i = i +1
oui
P
K Ti
N i ES
2
+ 3 (h + V2i )
L
L
i3
oui
= + i
FIN
j 220
non
FIN
Organigramme 11-2
252
9782100544639.indb 252
09/02/10 15:00
Lutilisation dune mthode de longueur darc ncessite tout dabord de dfinir lincrment darc Dl du premier incrment de charge (j = 1). Nous avons choisi pour
cet exemple 0.01. Cet incrment est nanmoins mis lchelle du problme la 1re
itration de convergence (i = 0) lors du calcul de D 0 . Ensuite, la procdure suit
celle dcrite au chapitre 11.1.4 qui consiste en une augmentation de lincrment de
charge . La monte en charge est donc proportionnelle son volution.
La convergence choisie tant toujours celle en forces ( = 0.001) , le reste a pour
expression:
h + V2 i
L
(11.66)
Ri = j + P N i
(11.67)
l j = {qi +1} {qi +1} soit dans notre cas l j = V2i +1 V2i +1 = V2i +1 .
Il est galement possible dincrmenter automatiquement cette quantit suivant le
rapport entre nombre ditrations dsir j+1 et celui obtenu au prcdent, soit:
n
l j +1
Niterd j +1
= l j
Niter j
(11.68)
9782100544639.indb 253
253
09/02/10 15:00
9
7
B
C
3
1
1
A
E
P (N)
KT (N/mm)
5
7
45
40
35
30
25
20
15
10
V2 (mm)
Figure 11.16 volution de KT et P en fonction de V2.
Nous avons vu au chapitre 6.1 que le moment flchissant pouvait tre exprim
en fonction de la courbure de la dforme partir de la relation (6.7) :
d 2v
EI z 2 = M z ( x ) . Par ailleurs, le moment flchissant vaut galement en tout
dx
254
9782100544639.indb 254
09/02/10 15:00
y
v(x)
EI z
d 2 v( x )
d 2 v( x )
=
M
(
x
)
=
v
(
x
)
EI
+ P v( x ) = 0 (11.69)
z
z
dx 2
dx 2
En posant =
P
, (11.69) devient:
EI z
d 2 v( x )
+ 2 v( x ) = 0
dx 2
(11.70)
v( x ) = A sin ( x ) + B cos ( x )
(11.71)
avec A et B constantes, et que les conditions aux limites permettent dcrire en:
x = 0, v(0) = 0 B = 0 ,
x = L, v(l ) = 0 A sin ( L ) = 0 .
Lquilibre se caractrise alors par:
A sin ( L ) = 0
ce qui permet de dduire:
L = k =
k
L
(11.72)
(11.73)
do
P
k
k 2 2 EI z
=
P=
EI z
L
L2
(11.74)
(11.75)
L2
9782100544639.indb 255
255
09/02/10 15:00
v ( x ) = {0 N 1 N 2 0 N 3 N 4 } {qe }
(11.77)
[G ] = 0
dN 1 dN 2
dN 3 dN 4
0
dx
dx
dx
dx
(11.78)
[k ] = [G ]
Ve
0 0
6
0
5
L
0
N
10
[ ] [G ] dVe =
0
0
L
0 6
L
0
10
0
L
10
2 L2
15
0
L
10
L2
30
0
0
0
0
0
0
0
6
5
L
10
0
6
5
L
10
0
L
10
L2
30
0
L
10
2 L2
15
(11.79)
[ K T ] = R ( ) ([k0 ] + [k ]) R ( )
(11.80)
256
9782100544639.indb 256
09/02/10 15:00
cos ( ) sin ( )
sin ( ) cos ( )
0
0
avec R ( ) =
0
0
0
0
0
0
0
0
0
0
0
0
1
0
0
0 cos ( ) sin ( )
0 sin ( ) cos ( )
0
0
0
0
0
.
0
0
Nous avons vu au chapitre 11.1.3 quil est possible de calculer la charge critique
dune structure au moyen dune analyse modale de flambement. Cette analyse
base sur une recherche de valeurs et vecteurs propres, permet de dterminer le
plus petit coefficient critique pondrant le chargement appliqu. Si ce coefficient
est infrieur 1, la charge critique est dj atteinte. La structure est donc considre
dans ce cas comme instable.
Bien que la prcision de ces coefficients soit directement lie la qualit du maillage et
donc celle des dformes modales, il est quand mme possible dvaluer correctement
la charge critique dune telle poutre avec un seul lment. Une premire solution
consisterait effectuer lanalyse modale en ne considrant que les degrsde libert
de rotation des extrmits libres. Cette approche donnant gnralement depitres
rsultats, il est prfrable dassocier les modes des degrs de libert de translation.
Aussi et en raison de la symtrie de la dforme, il est possible de rduire le modle
un seul lment de longueur gale la moiti de celle de la figure11.17.
V2
Y
x
2
U1
h
1
Dunod La photocopie non autorise est un dlit.
L1
Figure 11.18 Exemple 14-1: poutre simplement appuye (L1 = L/2).
9782100544639.indb 257
257
09/02/10 15:00
Reprenant les expressions de [k0] et [ks], respectivement (6.21) et (11.79), le systme1 aprs introduction des conditions aux limites se rduit :
4 EI
det L
6 EI
L2
2 L2
6 EI
L2 P 15
12 EI
L L
L
10
10 = 0
6
5
(11.81)
PCritique =
2 EI z 2 EI z
=
= 7772 N
L2
( 2L1 )2
(11.82)
Il serait nanmoins tout fait possible de retrouver prcisment ce rsultat en augmentant le maillage de la poutre.
Une fois le chargement critique dtermin, il est possible alors de complter lanalyse par un calcul non linaire gomtrique bas sur la mthode de longueur darc,
le dfaut h tant fix suivant le cas 1 mm ou 0.1 mm. Ce dfaut est en effet
ncessaire pour initier le processus non linaire et correspond en fait une pr
dformation de la poutre suivant la dforme du 1er mode. Lalgorithme utilis est
identique celui tabli pour lexemple prcdent. La tolrance e est cependant fixe
0.01 dans ce cas spcifique. Le calcul des efforts internes suit galement la mme
dmarche que celle du chapitre 6.2.4, la seule diffrence se situant dans la prise en
compte des grands dplacements dans le calcul de leffort normal, soit:
N = E xx
u2 u1 1 v2 v1 2
= E
+
2 L
L
(11.83)
258
9782100544639.indb 258
09/02/10 15:00
8 000
7 000
6 000
P (N)
5 000
4 000
Dfaut 1 mm
3 000
Dfaut 0,1 mm
2 000
1 000
0
0,00
0,20
0,40
0,60
0,80
1,00
1,20
U1 (mm)
Figure 11.19 Exemple 14-1: dplacement horizontal.
Ceci tant, il est parfois impossible datteindre cette valeur de charge critique.
Deux raisons peuvent expliquer ce constat. La premire concerne les effets dus aux
grands dplacements/grandes dformations. En effet et comme indiqu au chapitre
11.1.3, ceux-ci sont ngligs dans lanalyse modale de flambement entranant parfois une survaluation du coefficient critique. Lanalyse non linaire gomtrique
permettra alors de mesurer ces effets et donc de rvaluer le coefficient critique
gnralement la baisse. La deuxime est lie au phnomne de bifurcation qui
correspond une multiplicit des solutions dquilibre. Ce phnomne sera abord
dans lexemple 14-2.
8 000
7 000
6 000
P (N)
5 000
4 000
Dfaut 1 mm
3 000
Dfaut 0,1 mm
2 000
1 000
0
0,00
5,00
10,00
15,00
20,00 25,00
V2 (mm)
30,00
35,00
40,00
45,00
9782100544639.indb 259
259
09/02/10 15:01
V2
V1
2
y
U1
k2
h
X
k1
L
Lanalyse du comportement dune telle poutre sur appuis lastiques amne bien
videmment sinterroger sur linfluence des ressorts. Une premire approche,
certes un peu triviale, consiste prendre en compte des ressorts de rigidit infinie.
Dans ce cas, le calcul revient tudier la poutre sur appuis simples de lexemple
14-1. A contrario, sa dforme correspondra celle dune dnivellation dappui.
Les rotations nodales, supposes petites, seront alors gales langle correspondant
au dplacement diffrentiel des nuds 1 et 2. De ce fait, on pourra considrer
la tangente la dforme nulle aux extrmits de la poutre. Deux autres conditions aux limites peuvent galement tre tablies en 1 et 2 en posant que leffort
tranchant est gal, au signe prs, la force dans le ressort. Celui-ci est galement
quivalent la somme de la drive du moment flchissant et de la projection de la
charge P (cf. [9]). On a donc:
Ty =
dM z
dv
P
dx
dx
(11.84)
d 2v
Do daprs (6.7) EI z 2 = M z ( x ) :
dx
T y = EI z
d 3v
dv
P
3
dx
dx
(11.85)
260
9782100544639.indb 260
09/02/10 15:01
P
, (11.85) scrit:
EI z
v ( x ) + 2 v ( x ) =
v IV
(x ) + 2
Ty (x )
dT y
(11.86)
= q y ( x ) , on obtient:
dx
v ( x ) =
EI z
q y (x )
EI z
= 0
A
(11.87)
En posant =
qui est gale zro puisque la poutre nest pas charge transversalement. Finalement et en intgrant deux fois (11.87), cette quation diffrentielle devient:
v ( x ) + 2 v( x ) = R1 x + M1
(11.88)
Avec R1 et M1 constantes.
R x+M
La solution de (11.88) tant de la forme v( x ) = A sin ( x ) + B cos ( x ) + 1 2 1
R x+M
( x ) + B cos ( x ) + 1 2 1 avec A et B constantes, on dduit pour ses drives:
R
v ( x ) = A cos ( x ) B sin ( x ) + 12
(11.89)
v ( x ) = A 2 sin ( x ) B 2 cos ( x )
(11.90)
v ( x ) = A 3 cos ( x ) + B 3 sin ( x )
(11.91)
v ( x ) + 2 v( x ) = R1 =
T y ( x )
EI z
(11.92)
Comme Ty(x) vaut aux appuis k1v(0) et k1v(L), les deux premires conditions
aux limites scrivent:
Au nud 1
(11.93)
Au nud 2
R L + M
k
avec ki* = i .
EI z
9782100544639.indb 261
261
09/02/10 15:01
Les tangentes horizontales aux nuds tant considres comme nulles aux appuis,
on a galement:
En x = 0
v (0) = A +
R1
= 0
2
(11.95)
En x = L
v ( L ) = A cos ( L ) B sin ( L )
R1
= 0
2
(11.96)
1
cos ( L ) sin ( L )
2
0
1
k1*
*
k2* L
k2 sin ( L ) k2* cos ( L ) 1 2
A 0
B 0
= [C ] =
R1 0
M1 0
0
A
0
B
k1* R1
2 M1
k*
22
(11.97)
Le dterminant de [C] doit donc tre diffrent de 0 pour que (11.97) admette des
solutions non nulles. Dans le cas contraire, les racines de lquation det[C] = 0
correspondent aux valeurs de charge critique et donc aux diffrentes instabilits de
la structure. Dans notre cas dtude, cette rsolution revient poser que:
(11.98)
2
PCritique
=
k1 k2 L
k1 + k2
(11.99)
262
9782100544639.indb 262
09/02/10 15:01
.
L3
Point de bifurcation
Par dfinition, une bifurcation correspond la possibilit de suivre plusieurs trajectoires quilibres. Les solutions calcules aprs ce point de bifurcation peuvent
nanmoins traduire des tats stables ou instables. Nous avons vu au chapitre 11.2.4
que ces diffrents tats taient lis aux signes de la matrice de rigidit tangente.
Crisfield dveloppe dans (cf. [19]) les mthodologies permettant de calculer avec
prcision les positions de ces diffrents points singuliers: extrmums ou points de
bifurcation.
Sachant que lincrment de dplacement {Dq} peut tre exprim sous la forme
n
il tablit que:
A1 0 et Ai = 0 "i 1 pour les extrmums. Dans ce cas, il arrive montrer que
la longueur darc vaut {Dq}T{Dq} = Dl2 = A12. Ce rsultat a dailleurs t utilis
dans nos diffrents exemples pour prdire la variation de longueur darc
ncessaire lincrment de charge suivant (cf. programme Mathcad 11.2).
A1 = 0 et Ai 0 "i 1 pour les points de bifurcation.
Par ailleurs et dans ces deux cas, le dterminant de [KT] est daprs les conclusions
du chapitre 11.2.4 nul. Celle-ci admet donc ces endroits une valeur propre w1 de
mme valeur. Toujours daprs [19], la diffrentiation entre extrmums et points
de bifurcation seffectue partir du produit du vecteur des charges extrieures {F}
et du vecteur propre associ {f1}. {F}T {f1} est ainsi nul pour les points de bifurcation et non nul pour les extrmums.
Reprenant notre cas dtude, le point de bifurcation peut tre valu partir du
rsultat tabli au chapitre prcdent, soit:
k k L 2 EI
PCritique = min 1 2 ; 2
L
k1 + k2
9782100544639.indb 263
263
09/02/10 15:01
l 1T { d i }
l 1T { Lqi } + 0
(11.102)
P(N)
1 200
1 000
Dfaut 5 mm
800
Dfaut 1 mm
600
Dfaut 0.1 mm
400
200
0
0,00
0,10
0,20
0,30
U1 (mm)
0,40
0,50
2 000
1 800
1 600
1 200
1 000
Dfaut 5 mm
P(N)
1 400
800
Dfaut 1 mm
600
Dfaut 0.1 mm
400
200
0
7,50
6,50
5,50
4,50
3,50
V (mm)
2,50
1,50
0,50
264
9782100544639.indb 264
09/02/10 15:01
2 000
1 800
1 600
1 400
1 000
800
Dfaut 5 mm
600
Dfaut 1 mm
400
Dfaut 0,1 mm
P(N)
1 200
200
0
0,00
5,00
10,00
15,00
V2 (mm)
20,00
25,00
30,00
y
a
yy
xy
xx
xx
xy
yy
Figure 11.25 Plaque soumise un tat de contrainte plane.
9782100544639.indb 265
265
09/02/10 15:01
1
2 + 2 dV dV + 1
xx
xy2 dV
yy
xx
yy
E V
2E V
2 G V
(11.103)
yy
xx
dV
+
3
3
e
e
1 V
W =
2
2E
12 z M xy
12 z M xx 12 z M yy
2
1
dV
+
+
dV
(
)
2
e3
e3
e3
V
V
(11.104)
soit aprs transformation en une intgrale de surface:
W =
2
2
12 M xx + M yy dxdy 2 M xx M yy dxdy
(11.105)
2 E e 3 +2 (1 + ) M xy2 dxdy
12
avec dV = dxdydz et 3
e
e
2 2
z2
12
dz = 3
e
h 3 12
= 3 .
12 e
En remplaant maintenant les moments Mxx, Myy et Mxy par leurs expressions (9.5)
et (9.9), on obtient:
2 w 2 2 w 2
2w 2w
2
+
+
x 2 y 2
x 2 y 2
1
W = D
dxdy
2
2
2
+2 (1 ) x y
(11.106)
Reprenant la mme dmarche que celle de lexemple 9, nous savons que la surface
lastique dune plaque rectangulaire simplement appuye peut tre reprsente par
la double srie:
266
9782100544639.indb 266
09/02/10 15:01
m x
n y 1
w ( x , y ) = amn sin
sin
a
b
m =1 n =1
(11.108)
Aprs calcul des diffrentes drives secondes de (11.108) ncessaires la dtermination de (11.107), on obtient pour lnergie de dformation dune plaque de ce
type:
2
2
2
2
2
ab
2 m + n
W = D amn
8
b 2
a2
m =1 n =1
(11.109)
car:
2w 2w 2w 2
2 2
dx dy = 0 .
x
y
x
y
0 0
ab
ab
m x
2 m y dx dy = ab
.
sin
a
4
b
sin2
00
(11.110)
W =
a b
m2 2 2
D A2
+ 2
8
b
a2
(11.111)
Timoshenko tablit par ailleurs dans [9] que le travail de la force de compression
N xx = xx e est gal:
T =
ab
2
xx e w
dxdy
2 0 0 x
e m
= A 2 xx
2 a
2 ab
m x
m y
cos2 a sin2 b dxdy
00
(11.112)
9782100544639.indb 267
267
09/02/10 15:01
ab
Comme
m x
2 m y dxdy quivaut galement ab , on
sin
a
4
b
cos2
00
xx e m 2 2 b
8a
(11.113)
2
2
a b
m2 2 2
2 xx e m b (11.114)
=
D A2
+
A
8
8a
b 2
a2
Lquilibre critique tant caractris par la drive seconde du potentiel total par
rapport A, on obtient pour la contrainte critique de flambement:
2
2
2 E
crit = k E e
0
xx
12 (1 2 ) b
A 2
(11.115)
m
a
avec k = + ; = 1.
b
m
La variable m reprsentant le nombre dondes suivant x de la dforme de flambement, on recherchera la valeur minimale de k correspondant la plus petite des
contraintes critiques. Pour ce faire et en considrant m constant et f variable, on
drive la fonction k par rapport f soit:
m 1 m
dk
= 2 + 2 = 0 m2 = 2
d
m m
(11.116)
268
9782100544639.indb 268
09/02/10 15:01
10
9
8
7
6
5
4
3
m=1
m=2
m=3
m=4
m=5
m=6
2
1
0
0
9782100544639.indb 269
269
09/02/10 15:01
a
=2.
b
Le calcul revient appliquer les rsultats du chapitre 11.1.3 et donc rechercher les
valeurs et vecteurs propres du systme: det ([ K 0 ] + .[ K ]) = 0 . Llment utilis
est un Q4 sans intgration rduite.
Le but tant de faire une comparaison prcise des rsultats avec la solution thorique, la taille retenue pour les lments est de 10 cm de ct. Il nest bien videmment pas ncessaire davoir systmatiquement une telle densit de maillage.
Cependant, il est important de noter que cette prcision dpendra de la qualit de
reconstitution de la dforme modale et ce dautant plus que llment choisi est
pauvre.
270
9782100544639.indb 270
09/02/10 15:01
x
Figure 11.30 Exemple 15-1: plaque en compression simple.
La plaque tant simplement appuye, le contour de celle-ci est bloqu verticalement suivant z. Deux lignes dappuis suivant x et y sont galement ajoutes afin
dassurer la stabilit du modle. Enfin et pour terminer, une pression unitaire q de
100 MPa est applique au bord libre suivant y. Il est nanmoins important de noter
que celle-ci est introduite sous la forme dune charge linaire de compression gale
qe soit 2106 N/m.
La premire valeur propre1 trouve (0.759) par lanalyse modale faite avec Abaqus
correspond bien non seulement une contrainte de 75.9 MPa mais galement un
mode deux ondes suivant x.
9782100544639.indb 271
271
09/02/10 15:01
272
9782100544639.indb 272
09/02/10 15:01
xy
xy
Figure 11.34 Exemple 15-2: plaque en cisaillement.
k = 5.35 +
k=
4
2
5.35
+4
2
(11.117)
si 1
si 1
(11.118)
D 2
= 120 MPa
e b2
a
q
q
9782100544639.indb 273
273
09/02/10 15:01
Les conditions aux limites restant inchanges, on applique sur chacun des cts du
contour (cf. figure11.35) une charge linaire transverse de 2 .106 N/m quivalente
une pression q de 100 MPa.
274
9782100544639.indb 274
09/02/10 15:01
Le dversement appel galement flambement latral est un phnomne dinstabilit qui se caractrise, sous leffet de charges verticales, par une flexion latrale
accompagne dune rotation de torsion1. Les forces verticales produisant des dformations latrales2, il est possible de calculer les valeurs critiques associes cette
instabilit.
P1
P2
w
z
L
x
w
z
9782100544639.indb 275
275
09/02/10 15:01
y
y
P
h
b
L/2
L
Figure 11.39 Exemple 16: dversement dune poutre charge en son milieu.
PCR
3 2 4 2 EI y I L2 GJ 4 3
4 3
=
+
+
zg
zg
2
L3 I y 2 EI y 6 + 2
6 +2
6
3 2 4 2 EI y I L2 GJ 4 3
4 3
=
+
+
zg
zg
2
L3 I y 2 EI y 6 + 2
6 +2
6
(11.119)
h
position du point dapplication de la charge P.
2
Le phnomne de dversement tant associ de la flexion-torsion, on prfrera
utiliser la notion de moment critique de dversement qui dans notre cas est gal :
avec z g =
MCR =
PCR L
4
(11.120)
1. Appuis simples pour la flexion mais rotation de torsion et dplacements latraux empchs.
276
9782100544639.indb 276
09/02/10 15:01
3m
320 mm
300 mm
20 mm
11.5 mm
2.1 1011 N/m2
0.30
8.0769 1010 N/m2
152.2 cm2
29144 cm4
9003 cm4
171 cm4
2025000 cm6
L=
Hauteur h =
Largeur b =
Epaisseur semelles tf =
Epaisseur me tw =
Module dlasticit longitudinal E =
Coefficient de Poisson n =
Module dlasticit transversal G =
Section S =
Inertie Iz =
Inertie Iy =
Inertie de torsion J =
Inertie de gauchissement Iw =
+h/2 m
MCR =
PCR L
= 2.94 103 kN.m et PCR = 3919 kN
4
(11.121)
(11.122)
1. La charge ultime peut en effet tre atteinte bien avant la charge critique. On a tendance parfois loublier.
9782100544639.indb 277
277
09/02/10 15:01
Figure 11.40 Exemple 16: dforme (Abaqus) de 1er mode de flambement {f1}.
278
9782100544639.indb 278
09/02/10 15:01
Les coordonnes des nuds du modle sont ensuite modifies de telle manire
obtenir:
(11.123)
avec m choisi suffisamment petit pour ne pas perturber les rsultats du calcul1.
Le calcul qui tient compte de la non linarit gomtrique, est bas sur une mthode
de longueur darc.
Les tracs (figures 11.42 et 11.43) de lvolution de la charge P en fonction des
dplacements des jonctions me-semelle de la section mdiane2 mettent en vidence une croissance jusqu un maximum de 3858 kN, valeur qui reste proche
de la charge critique prvue en (11.121) et (11.122). Au-del, celle-ci diminue
progressivement avec laccroissement du dplacement latral. Le passage de cet
extrmum traduit le changement prononc de comportement de la poutre qui
passe de la flexion simple une flexion-torsion. On le vrifie bien en visualisant la
dforme aprs la charge critique qui traduit bien le dplacement transversal li au
flambement de la semelle suprieure.
Lanalyse des courbes des figures 11.42 et 11.43 permet galement de mettre
envidence les diffrences de comportement des semelles suprieure et infrieure.
P (kN)
4 000
3 500
3 000
2 500
2 000
Semelle sup
1 500
Semelle inf
1 000
500
0
0,01
0,01
0,03
0,05
0,07
0,09
9782100544639.indb 279
279
09/02/10 15:01
4 000
3 500
2 500
2 000
P (kN)
3 000
1 500
Semelle
sup
1 000
500
0
0,07
0,06
0,05
0,04
0,03
0,02
0,01
Comme prdit par la dforme du premier mode de flambement, la semelle infrieure reste presque immobile transversalement compar au dplacement latral
de la semelle suprieure. Cependant, on remarquera que le phnomne de dversement sinitie en fait bien avant le moment critique et ce bien videmment en raison
de la pr-dformation de la poutre.
Le choix du coefficient m contenu dans (11.123) est dailleurs la principale difficult de ce type danalyse puisquil doit initier la dforme sans trop la perturber.
280
9782100544639.indb 280
09/02/10 15:01
Bibliographie
[2] Bazergui, A., Bui-Quoc, T., Biron, A., McIntyre, G., Laberge, C. Rsistance des Matriaux 2e dition. s.l.: cole Polytechnique de Montral, 1993.
[3] Dhatt, G. et Touzot, G. Une prsentation de la mthode des lment finis. s.l.:
Maloine, 1981.
[4] Pilkey, Walter D. et Pilkey, Deborah F. Petersons Stress Concentration Factors 3rd Edition. s.l.: Wiley, 2008.
[6] Timoshenko, S. Theory of Plates and Shells 2nd Edition. s.l.: McGraw-Hill,
1959.
[7] Imbert, J.F. Analyse des structures par lments finis. s.l.: Cepadues, 1995.
[8] Batoz J.L., Dhatt G. Modlisation des structures par lment finis Vol.3
Coques. s.l.: Herms, 1992.
[9] Timoshenko, S.P. et Gere, J.M. Theory of Elastic Stability. s.l.: McGraw-Hill,
1963.
[10] Zienkiewicz, O.C. et Taylor, R.L. The Finite Element Method For Solid and
Structural Mechanics 6th Edition. s.l.: Elsevier Butterworth Heinemann,
2005.
[11] Crisfield, M.A. Non-Linear Finite Element Analysis of Solids and Structures
Essentials Vol.1. s.l.: John Wiley & Sons, 1991.
281
9782100544639.indb 281
09/02/10 15:01
Bibliographie
[15] Wempner, G.A. Discrete approximation related to non linear theory of solids.
International Journal of Solids and Structures. 1971, 7, pp. 1581-1599.
[17] De Souza Neto, E.A. et Feng, Y.T. On the determination of the path direction for arc-length methods in the presence of bifurcations and snap-backs.
Computer Methods in Applied Mechanics and Engineering. 1999, 179, pp.
81-89.
[19] Crisfield, M.A. Non-Linear Finite Element Analysis of Solids and Structures
Advanced Topics Vol. 2. s.l.: John Wiley & Sons, 1997.
[20] Vlassov, B.Z. Pices longues en voiles minces. s.l.: Eyrolles, 1962.
[21] Prat, M. Calcul des ouvrages gnraux de construction. s.l.: Herms, 1997.
282
9782100544639.indb 282
09/02/10 15:01
Complments en ligne
A
Toute utilisation des ressources mises disposition sur le site de Dunod
implique la prise de connaissance des conditions gnrales dutilisation.
1. Notice dutilisation
Vous venez de tlcharger les fichiers de donnes relatifs aux 16 exemples de louvrage Mthode des lments Finis, paru chez Dunod. Ces exemples ont t tablis
partir de trois logiciels: Abaqus1 6.8.2, Advance Structure2 3.1, et Mathcad3
versions 13 et 14.
Les rsultats de louvrage peuvent ainsi tre reproduits par les dtenteurs des
logiciels mentionns ci-dessus. Ceci tant, les utilisateurs dEffel, dit galement par Graitec, pourront utiliser directement ces diffrents modles, puisque
ceux-ci ont t construits partir de ce logiciel. Ceux dAdvance devront les
importer partir dEffel.
Les exemples se divisent en trois catgories:
exemples 1 4: lments barres et poutres;
exemples 5 12: lments membranes, plaques et coques;
exemples 13 16: analyse non linaire gomtrique.
Pour diffrentes raisons, chacun des logiciels ne couvre pas lintgralit des exemples.
Notamment, les n10, 12, 15 et 16 ne peuvent tre traits quavec Abaqus et/ou
Effel en raison de leur taille. Pour la mme raison, Mathcad ne permet de rsoudre
que les n1, 2, 3, 4, 5, 6, 9, 10, 13 et 144.
Les fichiers de donnes ont des extensions:
.cae pour Abaqus (Abaqus CAE) rpertoire ./FichiersAbaqus/,
.eff pour Effel/Advance rpertoire ./FichiersEffel/,
.xmcd pour Mathcad rpertoire ./FichiersMathcad/.
1 Abaqus est un logiciel de simulation numrique dit par la socit Dassault Systmes Simulia Corporation lien internet: www.simulia.com/products/abaqus_fea.html.
2 Advance Structure est un logiciel de calcul dit par la socit Graitec et orient vers les mtiers du Gnie Civil principalement lien internet: http://www.graitec.com/fr/ad.asp.
3 Mathcad est un logiciel de simulation mathmatique dit par la socit Parametric Technology Corporation lien
internet: http://www.ptc.com/appserver/mkt/products/home.jsp?k=3901.
4 Certains exemples Mathcad comme les n13 et 14 ncessitent une mmoire RAM importante ( 4 Go).
9782100544639.indb 283
283
09/02/10 15:01
Complments en ligne
1. Notice dutilisation
Concernant ces derniers, un gestionnaire (cf. figure 1) permettant de les excuter directement peut tre ouvert dans Mathcad partir du fichier MEF.hbk.
Chacun dentre eux comporte gnralement une version analytique (= rsultats
analytiques) et une version numrique intgrant les hypothses de calcul retenues
dans le livre.
284
9782100544639.indb 284
09/02/10 15:01
Complments en ligne
Article 1: responsabilits
La mise en ligne des diffrents fichiers de donnes mentionns dans la notice dutilisation ntant quune facilit offerte au lecteur pour suivre pas pas les dveloppements faits dans louvrage Mthode des lments Finis, ni lauteur, ni lditeur ne
pourront tre rendus responsables de dysfonctionnements lis lutilisation de ces
fichiers.
Les contenus mis disposition ont t mis en ligne aprs avoir t jugs raisonnablement conformes ltat de lart. Les rsultats correspondants restent nanmoins
soumis au jugement critique de lutilisateur. Lauteur et lditeur dclinent toute
responsabilit quant aux dommages ou pertes pouvant survenir suite lutilisation
dun de ces fichiers.
Article 2: licence
Lutilisation de ces fichiers suppose que le lecteur possde une licence dun, de deux
Article 3: portabilit
9782100544639.indb 285
285
09/02/10 15:01
9782100544639.indb 286
09/02/10 15:01
Index alphabtique
A
Analyse de flambement
257, 277
Approximation
polynomiale 33, 65,
158, 204
Arc (longueur d) 2, 234,
237, 249, 258, 263, 279
B
Batoz 237
Bifurcation 249, 254,
259, 263
Blocage en cisaillement.
VoirIntgration rduite
C
CAO 155
Charge critique 233, 255,
258
Charges nodales
quivalentes 70
Cisaillement transverse
185, 190
Coefficient
critique 233, 257, 259,
271, 277
de concentration 167
de Poisson 23, 140,
180, 210, 223
Continuit 158, 167
Contrainte
critique 273
dans les plaques 187
9782100544639.indb 287
D
Dformation plane 25,
127
Dversement 2, 275, 280
E
lancement 178, 198, 273
lment
barre 35, 38, 45, 47, 55
conforme 158
coque 222
isoparamtrique 123,
128, 140
plaque 196
poutre 63, 68, 76, 78,
89, 256
Q4 129, 152, 160, 167,
178, 196, 204, 210,
222, 270
Q8 164
Q9 179
ressort 46, 51
T3 141, 153, 167
T6 167, 171, 172
nergie de dformation
28, 33, 35, 67, 74, 191,
266
quation gnrale
des plaques 190, 193,
204
des poutres 63, 65, 80
quations dquilibre
de surface 19, 31
de volume 18, 31, 158
Extrmums 242, 254,
263
F
Flambement 2, 233, 254,
259, 268, 275, 279
Fonctions de forme 32,
228
G
Gauchissement 73, 275
Gauss (limination de) 10
Grandes dformations
228, 236
Grands dplacements 228
I
Intgration
numrique 123
rduite 198, 204, 219,
270
287
09/02/10 15:01
Index alphabtique
Isoparamtrique.
Voirlment
J
Jacobien 122, 129, 142,
197
K
Kirchhoff 185
M
Maillage 33, 151, 154,
160, 167, 174, 192,
204, 257
Matrice de rigidit
contraintes initiales
232, 256
lmentaire 35, 130
grandes dformations
232, 236, 243
membrane 135, 145
petites dformations
232
plaque 197
poutre 3D 76
poutre plane 67
ressort 46
N
Newton Raphson
(Mthode) 2, 234, 239,
245
Non linarit gomtrique
2, 279
P
Peterson (Courbes de)
167
Pivot 11, 41
Points de bifurcation.
VoirBifurcation
Points de Gauss 123, 199
Post critique 234
Post-processeurs 1
Pr-processeur 1
Processus de convergence
234
R
Rigidit tangente 233,
244, 253, 256
S
Section deffort tranchant
75, 78, 191, 209, 221
Snap back 249
Snap through 249
Solveurs 1
Striction 140, 180
T
Timoshenko 191, 267,
276
Torsion 73, 187, 275
Transformation
gomtrique 120
Travail des forces 33, 38,
239, 267
V
Voilement 2
288
9782100544639.indb 288
09/02/10 15:01