Exercice RDM

R
esistance des mat

eriaux :
elasticit
e,
m
ethodes
energ
etiques,
m
ethode des
el
ements nis
Rappels de cours
et exercices avec solutions
Yves Debard
Institut Universitaire de Technologie du Mans
D
epartement G
enie M
ecanique et Productique
http://iut.univ-lemans.fr/ydlogi/index.html
20 juin 2011
Table des mati`

eres
1 Elasticit
e
1.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . .
1.1.1 Deplacements et deformations . . . . . . . .
1.1.2 Contraintes . . . . . . . . . . . . . . . . . .
1.1.3 Loi de comportement ou loi constitutive . .
1.1.4 Cas particulier : etat de contraintes planes .
1.1.5 Formules mathematiques . . . . . . . . . .
1.2 Exercices . . . . . . . . . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
2 M
ethode des
el
ements nis : approche r
esistance des mat
eriaux
2.1 Rappels : resolution dun problème stationnaire . . . . . . . . . . . .
2.1.1 Partition des degres de liberte . . . . . . . . . . . . . . . . . .
2.1.2 Calcul des deplacements inconnus . . . . . . . . . . . . . . .
2.1.3 Calcul des reactions dappui . . . . . . . . . . . . . . . . . . .
2.2 Poutre soumise à un eort normal . . . . . . . . . . . . . . . . . . .
2.2.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.2.2 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.3 Treillis plans à nuds articules . . . . . . . . . . . . . . . . . . . . .
2.3.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.3.2 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.4 Poutre soumise à un moment de torsion . . . . . . . . . . . . . . . .
2.4.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.4.2 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
2.5 Flexion des poutres à plan moyen : modèle de Bernoulli . . . . . . .
2.5.1 Rappels : exion dans le plan {xy} . . . . . . . . . . . . . . .
2.5.2 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
3 M
ethodes
energ
etiques : poutres
3.1 Rappels . . . . . . . . . . . . . . . . . . . . . . .
3.1.1 Expression de lenergie de deformation en
mule de Clapeyron . . . . . . . . . . . . .
3.1.2 Theorème de reciprocite de Maxwell-Betti
3.1.3 Theorème de Castigliano . . . . . . . . .
3.1.4 Theorème de Menabrea . . . . . . . . . .
3.1.5 Energie
de deformation dune poutre . . .
3.1.6 Formules mathematiques utiles . . . . . .
3.2 Exercices . . . . . . . . . . . . . . . . . . . . . .
. . . . .
fonction
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
. . . . .
. . . . . .
des forces
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
. . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
. . . . . . .
appliquees :
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
. . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
1
1
1
2
2
3
4
6
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
25
25
25
26
26
27
27
28
44
44
45
54
54
55
58
58
60
. . .
for. . .
. . .
. . .
. . .
. . .
. . .
. . .
83
83
83
83
83
84
84
84
85
II
Exercices de resistance des materiaux
4 M
ethode des
el
ements nis
4.1 Rappels . . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.1.1 Energie
de deformation . . . . . . . . . . . . . . . .
4.1.2 Energie
cinetique . . . . . . . . . . . . . . . . . . . .
4.1.3 Energie potentielle et elements nis . . . . . . . . .

4.1.4 Modes propres . . . . . . . . . . . . . . . . . . . . .
4.2 Exercices . . . . . . . . . . . . . . . . . . . . . . . . . . . .
4.2.1 Assemblage . . . . . . . . . . . . . . . . . . . . . . .
ement de poutre droite soumis à un eort normal
4.2.2 El
4.2.3 Exercice : mise en equation . . . . . . . . . . . . . .
4.2.4 Exercice : mise en equation . . . . . . . . . . . . . .
4.2.5 Exercice : contraintes et energie de deformation . . .
4.2.6 Exercice : contraintes et energie de deformation . . .
ement de poutre droite soumis à un eort normal
4.2.7 El
4.2.8 Exercice : modes propres . . . . . . . . . . . . . . .
ement ni de torsion . . . . . . . . . . . . . . . . .
4.2.9 El
ement ni de exion : modèle de Bernoulli . . . .
4.2.10 El
4.2.11 Exercice : elasticite plane . . . . . . . . . . . . . . .
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
.
121
121
121
122
123
123
124
124
126
129
131
132
134
137
140
143
144
148
Chapitre 1
Elasticit
e
1.1
Rappels
Les deplacements et les deformations sont petits.
1.1.1
D
eplacements et d
eformations
Vecteur deplacement :
u(x, y, z)
{u} = v(x, y, z)
w(x, y, z)
u = M0 M ,
(1.1.1)
Tenseur des deformations :
1
xy
2
xx
1
[] = xy
2
1
xz
2
xy =
xx =
u
x
u v
+
y
x
yy
1
yz
2
,
1
xz
2
yz
zz
yy =
xz =
v
y
[]T = []
(1.1.2)
w
z
(1.1.3a)
zz =
u w
+
z
x
yz =
w v
+
y
z
(1.1.3b)

nx
Allongement unitaire en M dans la direction {n} = ny :

nz
(M, n) = {n}T [(M )] {n}
= n2x xx + n2y yy + n2z zz + nx ny xy + nx nz xz + ny nz yz
(1.1.4)
Glissement en M dans les directions orthogonales na et nb :

(M, na , nb ) = 2 {nb }T [(M )] {na }
{nb }T {na } = 0
(1.1.5)
Variation relative de volume :

V (M ) = tr [] = xx + yy + zz
(1.1.6)
1.1.2
Contraintes
Vecteur contrainte sur la facette n en M :

T (M, n) = n n + n
n = n T (M, n)
(1.1.7a)
T (M, n)2 = n2 + n 2
(1.1.7b)

nx
Soit {n} = ny un vecteur unitaire en M . Le vecteur contrainte sur la facette n en M est

nz
donne par la formule de Cauchy :

xx yx zx nx
Tx
Ty = xy yy zy ny

nz
xz yz zz
Tz
{T } = [(M )] {n}
(1.1.8)
o`
u [(M )] est le tenseur des contraintes en M .
Le tenseur des contraintes est symetrique :
[] = []T
soit xy = yx , xz = zx , yz = zy
(1.1.9)
La contrainte normale sur la facette n est :

n = {n}T [] {n}
= n2x xx + n2y yy + n2z zz + 2 nx ny xy + 2 nx nz xz + 2 ny nz yz
(1.1.10)
Soient 1 , 2 et 3 les trois contraintes principales en un point M dun solide. Les critères de
Rankine, Von Mises et de Tresca secrivent :
R = max(|1 |, |2 |, |3 |) E
VM =
(1.1.11a)
1
((1 2 )2 + (1 3 )2 + (3 2 )2 ) E
2
T = 2 max = max ( 1 , 2 , 3 ) min ( 1 , 2 , 3 ) E
1.1.3
(1.1.11b)
(1.1.11c)
Loi de comportement ou loi constitutive
Si le materiau est isotrope, la loi de comportement secrit :

1
(xx (yy + zz ))
E
1
= (yy (xx + zz ))
E
1
= (zz (xx + yy ))
E
xx =
yy
zz
xy =
xy
G
xz =
xz
G
yz =
yz
G
(1.1.12a)
G=
E
2 (1 + )
o`
u E et sont respectivement le module de Young et le coecient de Poisson du materiau.
(1.1.12b)

Elasticit
e
1.1.4
Cas particulier :
etat de contraintes planes
Le tenseur des contraintes se reduit à :
xx xy 0
[] = xy yy 0
0
0 0
(1.1.13)
do`
u lexpression du tenseur des deformations :
xx
1
[] =
xy
2
0
1
xy
2
yy
0
(1.1.14)
zz
et de la loi de comportement :
E
E
(xx + yy ) , yy =
(yy + xx )
2
1
1 2
E
= (xx + yy ) , xy = G xy , G =
E
2 (1 + )
xx =
zz
(1.1.15)
Les contraintes et les deformations principales sont :

}
xx + yy
1
1
2
=
(xx yy )2 + 4 xy
, 3 = 0
2
2
2
}
xx + yy
1
1
2
=
(xx yy )2 + xy
, 3 = zz
2
2
2
Les directions principales sont :
cos 1
sin 1
cos 1
{n1 } = sin 1
, {n2 } =
0
0

0
{n3 } = 0

1
avec
tan 1 =
(1.1.16)
(1.1.17)
1 xx
xy
(1.1.18)
Les critères de Rankine, Von Mises et de Tresca se reduisent à :

R = max( |1 | , |2 | ) E
VM = 12 + 22 1 2 E
T = 2 max = max( |1 2 | , |1 | , |2 | ) E

nx
Lallongement unitaire en M dans la direction {n} = ny se reduit à :

0
(M, n) = {n}T [(M )] {n} = n2x xx + n2y yy + nx ny xy
(1.1.19a)
(1.1.19b)
(1.1.19c)
(1.1.20)
1.1.5
Formules math
ematiques
Valeurs et vecteurs propres dune matrice symetrique de dimension deux à coecients reels :
Considerons la matrice symetrique [S] :
[
]
Sxx Sxy
[S] =
Sxy Syy
( [S]T = [S] )
(1.1.21)
Les valeurs propres Sn=1,2 et les vecteurs propres {n} sont les solutions de lequation :
[S]{n} = Sn {n}
{ }
]{ }
[
nx
nx
Sxx Sxy
= Sn
ny
ny
Sxy Syy
soit :
Sxx Sn
Sxy
Sxy
Syy Sn
avec n2x + n2y = 1
]{ } { }
0
nx
=
ny
0
(1.1.22)
(1.1.23)
Cette equation na de solution autre que la solution triviale nx = ny = 0 que si et seulement si :

[
]
Sxx Sn
Sxy
det
=0
Sxy
Syy Sn
(1.1.24)
do`
u lequation caracteristique :
2
=0
Sn2 (Sxx + Syy ) Sn + Sxx Syy Sxy
{z
}
|
|
{z
}
tr [S]=S1 +S2
(1.1.25)
det[S]=S1 S2
et les valeurs propres :

S1
S2
}
=
Sxx + Syy
1
2
(Sxx Syy )2 + 4 Sxy

2
2
Les vecteurs propres associes sont :

{
}
cos 1
{n1 } =
sin 1
avec :
tan 1 =
{
} {
}
cos 2
sin 1
{n2 } =
=
sin 2
cos 1
S1 Sxx
Sxy
tan 2 =
S2 Sxx
Sxy
(1.1.26)
(1.1.27)
(1.1.28)
Remarque : les deux directions principales sont orthogonales :

| 1 2 |=
tan 2 1 = tan 2 2 =
2 Sxy
Sxx Syy
tan 1 tan 2 = 1
Determinant dune matrice carree symetrique de dimension 3 :
S11 S12 S13

det S21 S22 S23
S31 S32 S33
[
]
[
]
[
]
S22 S23
S12 S13
S12 S13
= S11 det
S21 det
+ S31 det
S32 S33
S32 S33
S22 S23
2
2
2
= S11 S22 S33 S11 S23
S33 S12
S22 S13
+ 2 S12 S13 S23
(1.1.29)
(1.1.30)

Elasticit
e
Formules trigonometriques :
tan =
sin
cos
sin() = sin
cos() = cos ,
cos(1 + 2 ) = cos(1 ) cos(2 ) sin(1 ) sin(2 )
(1.1.32)
sin(1 + 2 ) = sin 1 cos 2 + cos 1 sin 2
(1.1.33)
1 cos 2
sin 2
sin2 =
, sin cos =
2
2
(
)
cos
= sin
2
1
2
1
3
cos 45 = sin 45 = =
, cos 60 =
, sin 60 =
2
2
2
2
1
3
cos 120 =
, sin 120 =
2
2
Si x et y sont petits devant lunite :
cos2 =
(1.1.31)
1 + cos 2
2
| x | 1 ,
(1.1.34)
(1.1.35)
(1.1.36)
(1.1.37)
| y | 1
(1.1.38a)
(1 + x)(1 + y) 1 + x + y
(1.1.38b)
on a les relations :
x
1+x1+
2
1
1x
1+x
sin x x ,
cos x 1
x2
2
tan x x
(1.1.38c)
1.2
Exercices
ELA 1 : vecteur contrainte sur une facette

En un point M dun solide, dans le repère orthonorme {, , k}, le tenseur des contraintes a pour
valeur :
100 40 20
[(M )] = 40 60 50 MPa
20
50 40
1. Faire un dessin qui montre la signication physique des composantes du tenseur des contraintes.
2. Soit le vecteur unitaire n de composantes :

1
1
2
{n} =
3
2
Sur la facette n :
(a) Calculer les composantes du vecteur contrainte T (M, n).
(b) Calculer la contrainte normale n .
(c) Calculer les composantes du vecteur cisaillement n , puis le module n du cisaillement.
Solution
Repr
esentation graphique des composantes du tenseur des contraintes
Les composantes en MPa du tenseur des contraintes dans le repère {, , k} :
composantes sur

k
T (M, ) T (M, ) T (M, k)
100
40
20
40
60
50
20
50
40
sont representees sur la gure ci-dessous.
Facette n
Les composantes du vecteur contrainte

100 40
T x
Ty = 40 60

20
50
Tz
sont (formule de Cauchy : {T (M, n)} = [(M )] {n}) :

20
1
60 20
1
1
50
2 =
60 = 20
MPa

3 3
66.67
40
2
200

Elasticit
e
On en deduit la contrainte normale (n = n T (M, n)) :

1
340
n = (60 120 + 400) =
= 37.78 MPa
9
9
les composantes du vecteur cisaillement (n = T (M, n) n n) :

nx 1 60 340 1 1 7.41
ny =
60
2 = 45.19 MPa
3
9 3
nz
200
2
41.48
et le module du cisaillement :
n = n =
7.412 + (45.19)2 + 41.482 = 61.78 MPa
Programme Scilab 5.3.1

// v e c t e u r c o n t r a i n t e s u r une f a c e t t e
// t e n s e u r d e s c o n t r a i n t e s dans l e r e pè r e { x y z }
sigma = [ 1 0 0 , 4 0 , 2 0 ; 4 0 , 6 0 , 5 0 ; 2 0 , 5 0 , 4 0 ]
// f a c e t t e {n}
n=[1/3;2/3;2/3]
T=sigma n
sigma n=n T
ta u n=Tsigma n n
norm( t a u n )
//
//
//
//
v e c t e u r c o n t r a i n t e : f o r m u l e de Cauchy
c o n t r a i n t e normale
vecteur cisaillement
module du c i s a i l l e m e n t
ELA 2 :
etat de contrainte en un point
En un point M dun solide, dans le repère orthonorme
valeur :
100 40
[(M )] = 40 80
0
0
{, , k}, le tenseur des contraintes a pour
0
0 MPa
0
2. Soit le vecteur unitaire n de composantes :

1
1
2
{n} =
5 0
Sur la facette n :
(a) Calculer les composantes du vecteur contrainte T (M, n).
(b) Calculer la contrainte normale n .
(c) Calculer les composantes du vecteur cisaillement n , puis le module n du cisaillement.
(d) Faire un dessin qui montre la facette, le vecteur contrainte, la contrainte normale et le
vecteur cisaillement.

3. Calculer les contraintes et les directions principales.
4. Faire un dessin qui montre la signication physique des contraintes et des directions principales.
5. Calculer les contraintes equivalentes de Von Mises et Tresca.
Solution
Les composantes en MPa du tenseur des contraintes dans le repère {, , k} :
composantes sur

k
T (M, ) T (M, ) T (M, k)
100
40
0
40
80
0
0
0
0
sont representees sur la gure ci-dessous.
Facette n
Les composantes du vecteur contrainte

100 40
T x
Ty = 40 80

Tz
0
0
sont (formule de Cauchy : {T (M, n)} = [(M )] {n}) :

0
1
20 8.94
1
1
0
2 =
120 = 53.67 MPa
5 0
5 0 0
0
On en deduit la contrainte normale (n = n T (M, n)) :

} 1 20
1 {
1 2 0
120 = 52.00 MPa
n =
5
5 0
le vecteur cisaillement (n = T (M, n) n n) :

20
1
32 14.31
nx
1
1
1
ny =
120 52
2 =
16
7.16
=
MPa

5
5 0
5 0 0
nz
0
et le module du cisaillement :
n = ||n || =
1
(322 + 162 ) = 16.00 MPa
5
Le vecteur contrainte T (M, n), la contrainte normale n et le cisaillement n sont reprsentees sur la
gure ci-dessous.

Elasticit
e
Contraintes et directions principales

Laxe k est direction principale :
n3 = k
3 = 0
Les contraintes principales 1 et 2 sont les solutions de lequation (det([ ] n [ I ]) = 0) :

]
[
100 n
40
det
= (100 n ) (80 n ) (40)2 = n2 180 n + 6400 = 0
40
80 n
do`
u:
1
2
}
=
180 1 2
180 4 6400
2
2
1 = 131.23 MPa ,
2 = 48.77 MPa
La position angulaire des directions principales n1 et n2 est denie par :

tan 1 =
131.23 100
40
tan 2 =
48.77 100
40
do`
u:
1 = 37.98 ,
2 = 52.02
Les contraintes et directions principales sont representees sur la gure ci-dessous.
Dans le repère principal, le tenseur des contraintes secrit :
131.23
0
0
48.77 0 MPa
[(M )]{n1 ,n2 ,n3 } = 0
0
0
0
Les contraintes de Tresca et de Von Mises sont :
T = max (|131.23 48.77|, |131.23 0|, |48.77 0|) = 131.23 MPa
1
VM =
((131.23 48.77)2 + (131.23 0)2 + (48.77 0)2 ) = 114.89 MPa
2
10
Remarque : calcul des contraintes et des directions principales à laide du cercle de Mohr.

sxx =100
syy=80
sxy=40
// t e n s e u r d e s c o n t r a i n t e s dans l e r e pè r e { x y z }
sigma =[ sxx , sxy , 0 ; sxy , syy , 0 ; 0 , 0 , 0 ]
// f a c e t t e {n}
n=[1/ sqrt ( 5 ) ; 2 / sqrt ( 5 ) ; 0 ]
T=sigma n
// v e c t e u r c o n t r a i n t e : f o r m u l e de Cauchy
sigma n=n T
// c o n t r a i n t e normale
ta u n=Tsigma n n // v e c t e u r c i s a i l l e m e n t
norm( ta u n )
// module du c i s a i l l e m e n t
// c o n t r a i n t e s p r i n c i p a l e s
d =0.5( sxx+syy ) ;
r =0.5 sqrt ( ( sxxsyy )2+4 sxy 2 ) ;
s i g 1=d+r
s i g 2=dr
s i g 3 =0
// a n g l e s en d e g re s d e s d i r e c t i o n s p r i n c i p a l e s a v e c l axe x
t 1=atan ( ( s i g 1 sxx ) / sxy ) 1 8 0 / %pi
// c o n t r a i n t e s de Von Mises e t Tresca
Von Mises=sqrt ( s i g 1 2+ s i g 2 2 s i g 1 s i g 2 )
Tr e sc a=max( s i g 1 , s i g 2 ,0) min( s i g 1 , s i g 2 , 0 )
ELA 3 : contraintes et directions principales

valeur :
80 40
[(M )] = 40 120
0
0
0
0 MPa
0
2. Calculer les contraintes principales et les directions principales.

Elasticit
e
11
3. Faire un dessin qui montre la signication physique des contraintes et des directions principales.
4. Calculer les contraintes equivalentes de Von Mises et Tresca.
Solution
Les composantes du tenseur des contraintes sont representees sur la gure ci-dessous.

n3 = k
3 = 0
Les contraintes principales 1 et 2 sont les solutions de lequation (det([ ] n [ I ]) = 0) :

[
]
80 n
40
det
= (80 n ) (120 n ) (40)2 = n2 200 n + 8000 = 0
40
120 n
do`
u:
1
2
}
=
200 1 2
200 4 8000
2
2
1 = 144.72 MPa ,
2 = 55.28 MPa

tan 1 =
144.72 80
40
tan 2 =
55.28 80
40
do`
u:
1 = 58.28 ,
2 = 31.72
Les contraintes et directions principales sont representees sur la gure ci-dessous.
Dans le repère principal, le tenseur des contraintes secrit :
144.72
0
0
55.28 0 MPa
[(M )]{n1 ,n2 ,n3 } = 0
0
0
0
12
Les contraintes equivalentes de Tresca et de Von Mises sont :

T = max (|144.72 55.28|, |144.72 0|, |55.28 0|) = 144.72 MPa
1
VM =
((144.72 55.28)2 + (144.72 0)2 + (55.28 0)2 ) = 126.49 MPa
2
Remarque : calcul des contraintes et des directions principales à laide du cercle de Mohr.

// c o n t r a i n t e s p l a n e s : c o n t r a i n t e s e t d i r e c t i o n s p r i n c i p a l e s
sxx=80
syy =120
sxy=40
// t e n s e u r d e s c o n t r a i n t e s
d =0.5( sxx+syy ) ;
r =0.5 sqrt ( ( sxxsyy )2+4 sxy 2 ) ;
s i g 1=d+r
s i g 2=dr
s i g 3 =0
// a n g l e s en d e g re s d e s d i r e c t i o n s p r i n c i p a l e s a v e c l axe x
Tr e sc a=max( s i g 1 , s i g 2 ,0) min( s i g 1 , s i g 2 , 0 )
ELA 4 : contraintes et directions principales

expression :
[(M )] = 0

Calculer les contraintes et les directions principales.

Elasticit
e
13
Solution
Les contraintes principales n sont les solutions de lequation :
n
= n (n2 2 ) ( n 2 ) + ( 2 + n )
det
n
= n3 + 3 n 2 + 2 2 = 0
do`
u:
1 = 2
2 = 3 =
La direction principale n1 est solution de lequation :

1 = 2
n1x 0
n1y = 0
1 = 2

n1z
0
1 = 2
do`
u:

1
1
1
{n1 } =
3 1
Les directions principales n2 et n3 sont les solutions de lequation :

2 = 3 =
nx 0
ny = 0
2 = 3 =
2 = 3 =
nz
0
do`
u:

1
1
1
{n2 } =
2 0

1
1
1
{n3 } =
6 2
Remarque : tout vecteur n = a n2 +b n3 avec a2 +b2 = 1 est direction principale associee à la contrainte
principale .
La contrainte de Von Mises et la contrainte de Tresca sont egales à :
VM = T = 3 | |
ELA 5 : d
eformations autour dun point
En un point M dun solide, dans le repère orthonorme {, , k}, le tenseur des deformations est egal
a:
`
150 120 0
[(M )] = 120 200 0 106
0
0
0
Un tel etat de deformation est dit etat de deformation plane par rapport à laxe k.
1. Faire un dessin qui montre la signication physique des composantes du tenseur des deformations.
14

2. Soient les deux directions orthogonales na et nb :

1
2
1
1
2
1
{na } =
, {nb } =
5 0
5 0
na nb = 0
(a) Calculer lallongement unitaire dans les directions na et nb et le glissement dans ces deux
directions.
(b) Faire un dessin qui montre la signication physique de ces deformations.
(c) Expression du tenseur des deformations dans le repère {na , nb , k}.
3. Calculer les deformations et les directions principales. Faire un dessin qui montre la signication
physique des deformations principales.
4. Calculer la variation relative de volume V .
Solution
Les composantes du tenseur des deformations dans le repère {, , k} sont representees sur la gure
ci-dessous.
D
eformations dans les directions na et nb
Lallongement unitaire dans la direction na est :
(M, na ) = {na }T [(M )] {na }

150 120 0
1
{
}
1
6 1
2 = 286 106
= 1 2 0 120 200 0 10

5
5
0
0
0
0
Lallongement unitaire dans la direction nb est :
(M, nb ) = {nb }T [(M )] {nb }

2
} 150 120 0
1
1 {
1 = 64 106
= 2 1 0 120 200 0 106
5
5 0
0
0
0
Le glissement dans les deux directions orthogonales na et nb est :
(M, na , nb ) = 2 {nb }T [(M )] {na }

150
120
0
1
}
1
1 {
2 = 104 106
= 2 2 1 0 120 200 0 106
5
5 0
0
0
0
Dans le repère {na , nb , k}, le tenseur des deformations secrit :
286 52 0
[(M )]{na ,n ,k} = 52 64 0 106
b
0
0 0

Elasticit
e
15
Les composantes du tenseur des deformations dans le repère {na , nb , k} sont representees sur la gure
ci-dessous.
D
eformations et directions principales
n3 = k
3 = 0
Les deformations principales 1 et 2 sont les valeurs propres de la matrice :

[
]
150 120
106
120 200
do`
u:
1
2
(
=
150 + 200 1
(150 200)2 + 4 (120)2

2
2
1 = 298 106
106
2 = 52 106

tan 1 =
298 150
120
tan 2 =
52 150
120
do`
u:
1 = 50.88 ,
2 = 39.12
Dans le repère principal, le tenseur des deformations secrit
298 0
[(M )]{n1 ,n2 ,n3 } = 0 52

0
0
0
0 106
0
Les composantes du tenseur des deformations dans le repère principal sont representees sur la gure
ci-dessous.
La variation relative de volume est :

V (M ) = tr [(M )] = (150 + 200 + 0) 106 = (286 + 64 + 0) 106 = (298 + 52 + 0) 106 = 350 106
16

// de f o r m a t i o n s a u t o u r d un p o i n t
// de f o r m a t i o n s E6 dans l e r e pè r e { x y z }
exx=150
eyy=200
gxy=240
// t e n s e u r d e s de f o r m a t i o n s dans l e r e pè r e { x y z }
e p s =[ exx , gxy / 2 , 0 ; gxy / 2 , eyy , 0 ; 0 , 0 , 0 ]
// d i r e c t i o n s o r t h o g o n a l e s {n} e t { t }
n=[1/ sqrt (5); 2/ sqrt ( 5 ) ; 0 ]
t=[n ( 2 ) ; n ( 1 ) ; 0 ]
// a l l o n g e m e n t s u n i t a i r e s dans l e s d i r e c t i o n s {n} e t { t }
enn=n e p s n
e t t=t e p s t
// g l i s s e m e n t d e s d i r e c t i o n s o r t h o g o n a l e s {n} e t { t }
gnt=2t e p s n
// t e n s e u r d e s de f o r m a t i o n s dans l e r e pè r e { n t z }
e p s n t z =[enn , gnt / 2 , 0 ; gnt / 2 , e t t , 0 ; 0 , 0 , 0 ]
// de f o r m a t i o n s p r i n c i p a l e s
d=(exx+eyy ) / 2 ;
r=sqrt ( ( exxeyy )2+ gxy 2 ) / 2 ;
e p s 1=d+r
e p s 2=dr
// d i r e c t i o n s p r i n c i p a l e s ( a n g l e s en d e g re s )
an1=atan ( 2 ( eps1exx ) / gxy ) 1 8 0 / %pi
an2=atan ( 2 ( eps2exx ) / gxy ) 1 8 0 / %pi
// de f o r m a t i o n v o l u m i q u e
eps V=trace ( e p s )
ELA 6 :
etat de contrainte uniaxial
Une poutre de section droite constante est soumise à une contrainte de traction .
1. Donner lexpression du tenseur des contraintes dans le repère {, , k}.

2. Considerons la facette n en M .
(a) Calculer les composantes du vecteur contrainte sur cette facette.
(b) Calculer la valeur du cisaillement maximal.

Elasticit
e
17
(c) Calculer la valeur de pour laquelle lallongement unitaire est nul.

Application numerique : = 0.3.
dV
3. Calculer la variation relative de volume de la barre
en fonction de E, et .
V
Solution
Dans le repère {, , k}, les composantes du tenseur des contraintes et du tenseur des deformations
sont :
0 0
1 0
0
0 0
[] = 0 0 0 , [] =
E
0 0 0
0 0
Les composantes des vecteurs n et t sont :
cos
{n} = sin
sin
cos
{t} =
Le vecteur contrainte sur la facette n est egal à (formule de Cauchy) :
cos
0
{T } = [] {n} =
0
do`
u lexpression du cisaillement :
1
n = {t}T {T } = cos sin = sin(2 )
2
Le cisaillement est maximal pour tel que :
dn
= 0 soit
d
cos(2 ) = 0
do`
u:
,
2
Lallongement unitaire dans la direction n est :
=
nn = {n}T [] {n} =
max =
||
2
(cos2 sin2 )
E
et sannule pour tel que :

cos2 sin2 = 0
soit
tan2 =
Pour = 0.3, on obtient = 61.29.

La deformation volumique est :
V = tr [] =
1 2
18
ELA 7 : rosette `
a 45 degr
es
La rosette à 45 degres representee sur la gure ci-contre est collee sur un
reservoir au point critique.
Les caracteristiques du materiau sont :
module de Young : E = 210000 MPa
coecient de Poisson : = 0.27
limite elastique : E = 300 MPa
On mesure :
a = 100 106
= 450 106
b
c = 400 106
1. Calculer les composantes du tenseur des deformations.

2. En deduire les composantes du tenseur des contraintes.
3. Calculer zz .
4. Calculer les contraintes principales 1 , 2 et 3 .
5. Calculer les deformations principales 1 , 2 et 3 .
6. Calculer la position angulaire 1 de la direction principale n1 avec la jauge a.
7. La pression à linterieur du reservoir est egale à p = 50 MPa. Calculer la pression maximale
admissible pmax (on utilisera les critères de Tresca et Von Mises).
Solution
En M , letat de contrainte est plan par rapport à laxe k. Les allongements unitaires a , b et c
secrivent en fonction des composantes du tenseur des deformations :
jauge a :
jauge b :
jauge c :
On en deduit :

1
{na } = 0

0
2/2
{nb } =
2/2

0
{nc } = 1

0
a = (M, na ) = xx
b = (M, nb ) =
1
1
1
xx + yy + xy
2
2
2
c = (M, nc ) = yy

Elasticit
e
19
les composantes non nulles du tenseur des deformations :

xx = a = 100 106
yy = c = 400 106
xy = 2 b a c = 2 (450 106 ) (100 106 ) 400 106 = 1200 106
zz =
0.27
(xx + yy ) =
(100 106 + 400 106 ) = 111 106
1
1 0.27
les composantes non nulles du tenseur des contraintes :

xx =
yy =
E
210000
(xx + yy ) =
(100 + 0.27 400) 106 = 1.81 MPa
2
1
1 0.272
210000
E
(yy + xx ) =
(400 + 0.27 (100)) 106 = 84.49 MPa
2
1
1 0.272
E
210000
xy =
xy =
(1200 106 ) = 99.21 MPa
2 (1 + )
2 (1 + 0.27)
Les composantes du tenseur des deformations et du tenseur des contraintes dans

sont :
100 600
0
1.81
99.21
6
0
[(M )] = 600 400
10
, [(M )] = 99.21 84.49
0
0
111
0
0
le repère {, , k}
0
0 MPa
0

n3 = k
3 = 0
Les contraintes principales 1 et 2 sont :

}
xx + yy
1
=
2
2
1.81 + 84.49
=
2
3 = 111 106
1
2
(xx yy )2 + 4 xy
2
1
(1.81 84.49)2 + 4 (99.21)2

2
soit :
1 = 150.63 MPa ,
2 = 64.33 MPa
La position angulaire de la direction principale n1 est donnee par lequation :

tan 1 =
150.63 1.81
1 xx
=
xy
99.21
do`
u:
1 = 56.31
Les deformations principales 1 et 2 sont :
1 =
1
1
(1 2 ) =
(150.63 0.27 (64.33)) = 800 106
E
210000
20

2 =
1
1
(2 1 ) =
(64.33 0.27 150.63) = 500 106
E
210000

T = max(|1 2 |, |1 |, |2 |) = 150.63 (64.33) = 214.96 MPa
VM =
12 + 22 1 2 =
150.632 + (64.33)2 150.63 (64.33) = 191.10 MPa
La pression maximale admissible est donc :

critère de Tresca : pmax = p
critère de Von Mises : pmax
E
300
= 50
= 69.78 MPa
T
214.96
E
300
=p
= 50
= 78.49 MPa
VM
191.10

Le programme se trouve dans le cours delasticite.
ELA 8 : rosette `
a 120 degr
es
La rosette à 120 degres representee sur la gure ci-contre est collee en un
point M dune structure.
Les caracteristiques du materiau sont :
module de Young : E = 210000 MPa
coecient de Poisson : = 0.3
On mesure :
a = 900 106
= 400 106
b
c = 50 106
1. Calculer les composantes du tenseur des deformations.

2. En deduire les composantes du tenseur des contraintes.
3. Calculer les contraintes principales 1 , 2 et 3 .
4. Calculer les deformations principales 1 , 2 et 3 .
5. Calculer la position angulaire 1 de la direction principale n1 avec la jauge a.
6. Calculer la contrainte equivalente de Von Mises et la contrainte equivalente de Tresca.
Solution
En M , letat de contrainte est plan par rapport à laxe k. Les allongements unitaires a , b et c
secrivent en fonction des composantes du tenseur des deformations :

Elasticit
e
21
jauge a :

1
{na } = 0

0
jauge b :
1/2
{nb } =
3/2
jauge c :
a = (M, na ) = xx
1
3
3
b = (M, nb ) = xx + yy
xy
4
4
4
1/2
{nc } = 3/2
1
3
3
c = (M, nc ) = xx + yy +
xy
4
4
4
On en deduit :
les composantes non nulles du tenseur des deformations :
xx = a = 900 106
yy =
1
(2 b + 2 c a ) = 533 106
3
2
xy = (c b ) = 520 106
3
(xx + yy ) = 157 106

zz =
1
les composantes non nulles du tenseur des contraintes :

xx =
E
210000
(xx + yy ) =
(900 + 0.3 533) 106 = 170.77 MPa
2
1
1 0.32
yy =
E
210000
(yy + xx ) =
(533 + 0.3 (900)) 106 = 60.77 MPa
2
1
1 0.32
xy =
E
210000
xy =
(520 106 ) = 41.97 MPa
2 (1 + )
2 (1 + 0.3)
Les composantes du tenseur des deformations et du tenseur des contraintes dans

sont :
170.77 41.97
900 260 0
60.77
0 106
, [(M )] = 41.97
[(M )] = 260 533
0
0
0
0
157

n3 = k
3 = 0
3 = 157 106
le repère {, , k}
0
0 MPa
0
22
Les contraintes principales 1 et 2 sont :

}
xx + yy
1
1
2
=
(xx yy )2 + 4 xy
2
2
2
170.77 + 60.77 1
=
(170.77 60.77)2 + 4 (41.97)2
2
2
soit :
1 = 68.14 MPa ,
2 = 178.14 MPa
La position angulaire de la direction principale n1 est donnee par lequation :

tan 1 =
68.14 (170.77)
1 xx
=
xy
41.97
do`
u:
1 = 80.04
Les deformations principales 1 et 2 sont :
1 =
1
1
(1 2 ) =
(68.14 0.27 (178.14)) = 579 106
E
210000
2 =
1
1
(2 1 ) =
(178.14 0.27 68.14) = 946 106
E
210000

T = max(|1 2 |, |1 |, |2 |) = 68.14 (178.77) = 246.28 MPa
VM =
12 + 22 1 2 =
68.142 + (178.14)2 68.14 (178.14) = 220.27 MPa

// De p o u i l l e m e n t d une r o s e t t e d e x t e n s o m e t r i e
// a n g l e de l a r o s e t t e : 45 , 60 ou 120 d e g re s
p h i =120
i f phi =[45 ,60 ,120]
then halt ( p h i d o i t e t r e e g a l à 4 5 , 60 ou 120 d e g re s )
end
// m a te r i a u
E=210000
nu=0.3
// module de Young
// c o e f f i c i e n t de P o i s s o n
// a l l o n g e m e n t s u n i t a i r e s mesure s
e p s a=900E6
epsb =400E6
e p s c=50E6
// c a l c u l d e s de f o r m a t i o n s e x x , eyy , e z z e t gxy
// l axe x e s t l axe de l a j a u g e a
r p h i=p h i %pi / 1 8 0 ;
// a n g l e en r a d i a n s
c=cos ( r p h i ) ; s=sin ( r p h i ) ;
Methode des elements finis
23
exx=e p s a ;
eyy=( e p sa ( 2 c 21)+2epsb (1 2 c 2 )+e p s c ) / ( 2 s 2 ) ;
gxy=( e p s a (1 4 c 2 ) +4c 2 epsbe p s c ) / ( 2 s c ) ;
e z z=nu ( exx+eyy ) /(1nu ) ;
// c a l c u l d e s c o n t r a i n t e s Sxx , Syy e t Sxy
c=E/(1nu 2 ) ;
sxx=c ( exx+nu eyy ) ;
syy=c ( eyy+nu exx ) ;
G=E/2/(1+nu ) ; // module d e l a s t i c i t e t r a n s v e r s a l
sxy=G gxy ;
// t e n s e u r d e s c o n t r a i n t e s
d =0.5( sxx+syy ) ;
r =0.5 sqrt ( ( sxxsyy ) 2+4 sxy 2 ) ;
s i g 1=d+r
s i g 2=dr
s i g 3 =0
// a n g l e en d e g r e s de l a d i r e c t i o n p r i n c i p a l e n1 a v e c l a j a u g e a
t 1=atan ( ( s i g 1 sxx ) / sxy ) 180/ %pi
// t e n s e u r d e s de f o r m a t i o n s
e p s i l o n =[ exx , 1 / 2 gxy , 0 ; 1 / 2 gxy , eyy , 0 ; 0 , 0 , e z z ]
// de f o r m a t i o n s p r i n c i p a l e s
ep s 1 =( s i g 1 nu s i g 2 ) /E
ep s 2 =( s i g 2 nu s i g 1 ) /E
ep s 3=e z z
// de f o r m a t i o n v o l u m i q u e dV/V
epsV=trace ( e p s i l o n )
Tr e s c a=max( s i g 1 , s i g 2 , 0 )min( s i g 1 , s i g 2 , 0 )
//
e d i t i o n s
=====================================================
p r i n t f ( \n\ n R o s e t t e d e x t e n s o me t r i e à %3d d e g re s \n\n , p h i ) ;

p r i n t f ( Module de Young : %8. 0 f MPa\n\n ,E)
p r i n t f ( C o e f f i c i e n t de P o i s s o n : %5. 2 f \n\n , nu )
p r i n t f ( Mesures : ep s a = %5. 0 f E6 , epsb = %5. 0 f E6 , e p s c = %5. 0 f E6\n\n ,
e p s a 1E6 , epsb 1E6 , e p s c 1E6 )
p r i n t f ( De f o r m a t i o n s p r i n c i p a l e s : ep s 1 = %5. 0 f E6 , e p s 2 = %5. 0 f E6 , e p s 3
= %5. 0 f E6\n\n , e p s 1 1E6 , e ps 2 1E6 , e p s 3 1E6 )
p r i n t f ( De formation volumique : epsV = %5. 0 f E6\n\n , epsV 1E6 )
p r i n t f ( C o n t r a i n t e s p r i n c i p a l e s : s i g 1 = %5. 2 f MPa , s i g 2 = %5. 2 f MPa , s i g 3 =
0\n\n , s i g 1 , s i g 2 )
p r i n t f ( Angle de l a d i r e c t i o n p r i n c i p a l e n1 avec l a j a u g e a : %5. 2 f d e g re s \n\n ,
t1 )
p r i n t f ( C o n t r a i n t e s de Von Mises = %5. 2 f MPa\n\n , Von Mises )
p r i n t f ( C o n t r a i n t e s de Tresc a = %5. 2 f MPa\n\n , T r e s c a )
Chapitre 2
M
ethode des
el
ements nis : approche
r
esistance des mat
eriaux
2.1
Rappels : r
esolution dun probl`
eme stationnaire
Les equations dequilibre de la structure secrivent :

[K] {U } = {Fnod } + {F }
(2.1.1)
o`
u:
[K] est la matrice de rigidite.
{U } est le vecteur deplacement.
{Fnod } est le vecteur des forces nodales.
{F } est le vecteur force equivalent aux charges reparties et aux gradients thermiques.
2.1.1
Partition des degr

es de libert
e
Ce système dequations lineaires na pas de solution : la matrice de rigidite est une matrice singulière
(det[K] = 0). Imposons un certain nombre de deplacements et numerotons les degres de liberte de
facon à obtenir la partition suivante :
{
}
{UL } = ?
{U } =
(2.1.2)
{UP }
o`
u:
{UL } est le vecteur des deplacements inconnus.
{UP } celui des deplacements connus (nuls ou non nuls).
Cette partition du vecteur deplacement induit une partition des vecteurs force et de la matrice de
rigidite :
[
]{
} {
} {
}
[KLL ] [KLP ] {UL } = ?
{Fnod,L }
{FL }
=
+
(2.1.3)
[KP L ] [KP P ]
{UP }
{Fnod,P } = ?
{FP }
o`
u:
26

{FL } et {Fnod,L } sont les vecteurs des forces connues.
{Fnod,P } est le vecteur des reactions dappui.
2.1.2
Calcul des d
eplacements inconnus
Du système dequations precedent, on deduit :

[KLL ] {UL } + [KLP ] {UP } = {Fnod,L } + {FL }
(2.1.4)
[KLL ] {UL } = {Fnod,L } + {FL } [KLP ] {UP }
(2.1.5)
do`
u:
Le deuxième membre de lequation ci-dessus est connu. Si le nombre de liaisons est susant, la
matrice [KLL ] nest pas singulière (det[KLL ] = 0) et les deplacements inconnus sont :
{UL } = [KLL ]1 ( {Fnod,L } + {FL } [KLP ] {UP } )
2.1.3
(2.1.6)
Calcul des r
eactions dappui
Tous les deplacements etant connus, les reactions dappui sont calculees à laide de lequation :
{Fnod,P } = [KP L ] {UL } + [KP P ] {UP } {FP }
(2.1.7)
2.2
2.2.1
27
Poutre soumise `
a un eort normal
Rappels
Soit un element de poutre droite i j, daxe x, de longueur L, de section droite constante, soumis à
un eort normal N (x).
A est laire de la section droite.

E et sont respectivement le module de Young et le coecient de dilatation du materiau.
Lelement porte sur toute sa longueur une force lineairement repartie dintensite lineique (px , 0, 0)
avec :
(
x)
x
x
px (x) = pxi 1
(2.2.1)
+ pxj = pxi + ( pxj pxi )
L
L
L
et subit une variation de temperature T constante dans lelement.
ui et uj sont les deplacements nodaux.
Ni et Nj sont les eorts normaux dans les sections i et j.
Lintegration de lequation dequilibre :
dN
+ px = 0
dx
avec N (0) = Ni
(2.2.2)
conduit à lexpression de leort normal :

N (x) = Ni pxi x ( pxj pxi )
x2
2L
(2.2.3)
La contrainte normale dans la poutre est egale à :

xx (x) =
N (x)
A
(2.2.4)
Lintegration de la relation de comportement (equation de la deformee) :

xx =
du
N
=
+ T
dx
EA
avec u(0) = ui
conduit à lexpression de la deformee :

(
)
1
x2
x3
u(x) = ui +
Ni x pxi
( pxj pxi )
+ T x
EA
2
6L
(2.2.5)
(2.2.6)
Des conditions aux limites N (L) = Nj et u(L) = uj , on deduit :

{fnod } = [ k ] {u} {f } {fth }
(2.2.7a)
28
avec
{
} {
}
N (0)
Ni
=
N (L)
Nj
{
}
L 2 pxi + pxj
{f } =
6 pxi + 2 pxj
{fnod } =
[k] =
[
]
EA 1 1
L 1 1
{ }
ui
uj
{u} =
{fth } = EA T
(2.2.7b)
{ }
1
1
(2.2.7c)
Remarque : lequilibre de lelement secrit :

Ni + Nj +
2.2.2
L
(pxi + pxj ) = 0
2
Exercices
MEF NOR 1
La poutre representee sur la gure est discretisee en trois elements de meme longueur L.
Soient E et A respectivement le module dYoung du materiau et laire de la section droite.

Les sections 1 et 4 sont encastrees.
La section 2 porte une force dintensite 2 F et la section 3 une force dintensite F .
1. Eectuer letude elementaire.
2. Eectuer lassemblage.
3. Calculer les deplacements nodaux.
4. Calculer les actions de liaison et verier lequilibre de la structure.
Solution
Lassemblage conduit à la relation [K] {U } = {Fnod } :
F1x =?
u1 = 0
1
1
u2 2 F
EA
1
1
+
1
1
=
u3
1 1 + 1 1
F
L
F4x =?
u4 = 0
1
1
Les deplacements inconnus sont les solutions de lequation [KLL ] {UL } = {Fnod,L } :
[
]{ } { }
EA 2 1 u2
2F
=
1
2
u
F
L
3
(2.2.8)
29
do`
u:
u2 =
5 FL
3 EA
u3 =
4 FL
3 EA
Les actions de liaisons sont les solutions de lequation [KP L ] {UL } = {Fnod,P } :
[
]{ } {
}
EA 1
u2
F1x
=
1 u3
F4x
L
do`
u:
5
F1x = F
3
4
F4x = F
3
Lequilibre de la poutre est verie :

5
4
F1x + F2x + F3x + F4x = F + 2 F + F F = 0
3
3
Leort normal N (x) et le champ de deplacements u(x) sont representes sur la gure ci-dessous.
MEF NOR 2
Considerons la poutre daxe x representee sur la gure. Soit E le module de Young du materiau. La
rigidite lineique est egale à EA entre les nuds 1 et 2, 2 EA entre les nuds 2 et 4.
La poutre est encastree en 1 et 4. Elle est sollicitee au nud 2 par une force (F, 0, 0) et au nud 3
par une force (3 F, 0, 0) avec F > 0.
1. Calculer les deplacements nodaux et les actions de liaison.
2. Representations graphiques : tracer u(x), N (x) et xx (x)
Solution
Les matrices de rigidite elementaires sont :
[
]
[
]
EA 1 1
2 EA 1 1
[k12 ] =
, [k23 ] =
1 1
L 1 1
L
[
]
2 EA 1 1
[k34 ] =
3 L 1 1
30

1
1
u1 = 0
F1x =?
u2 =? F
EA
1
1
+
2
2
=
2 2 + 2/3 2/3
u3 =?
3F
L
2/3
2/3
u4 = 0
F4x =?
[
]{ } { }
EA 3 2
u2
F
=
3F
L 2 8/3 u3
do`
u:
u2 =
13 F L
6 EA
u3 =
11 F L
4 EA
Les actions de liaisons sont les solutions de lequation [KP L ] {UL } = {Fnod,P } :
[
]{ } {
}
EA 1
u2
F1x
=
2/3 u3
F4x
L
do`
u:
F1x =
13
F
6
F4x =
11
F
6

F1x + F2x + F3x + F4x =
13
11
F + F + 3F
F =0
6
6
Leort normal N (x) et la contrainte normale xx (x) sont representes sur la gure ci-dessous.
MEF NOR 3
La poutre representee sur la gure est constituee de deux troncons de meme longueur L.
31
Soit E le module de Young du materiau.

Laire de la section droite est egale à 2A entre les nuds 1 et 2 et A entre les nuds 2 et 3.
La section 1 est encastree.
La section 2 porte une force dintensite 2 P et la section 3 une force dintensite P .
2. Pour chaque element, calculer les eorts nodaux.
3. Representations graphiques : tracer N (x), xx (x) et u(x).
Solution
[
]
2 EA 1 1
[k12 ] =
1 1
L
[
]
EA 1 1
[k23 ] =
L 1 1

2
2
u1 = 0 F1x =?
EA
2 2 + 1 1
u2
2P
=
L
1
1
u3
P
[
]{ } { }
EA 3 1 u2
2P
=
u3
P
L 1 1
do`
u:
3 PL
5 PL
, u3 =
2 EA
2 EA
La reaction dappui est solution de [KP L ] {UL } = {Fnod,P } :
{ }
] u2
{
}
EA [
1
= F1x
u
L
3
u2 =
do`
u:
F1x = 3 P
F1x + F2x + F3x = 3 P + 2 P + P = 0
Leort normal, la contrainte normale et le champ de deplacements sont representes sur la gure
ci-dessous.
32
MEF NOR 4
Considerons la poutre daxe x representee sur la gure. Soit E le module de Young du materiau. La
rigidite lineique est egale à 2 EA entre les nuds 1 et 2, EA entre les nuds 2 et 3 et à 3 EA entre
les nuds 3 et 4.
La poutre est encastree en 1 et 4. Elle est sollicitee au nud 2 par une force (F, 0, 0) et au nud 3
par une force (2 F, 0, 0) avec F > 0.
2. Pour chaque element, calculer leort normal aux deux extremites.
Solution
[
]
[
]
2 EA 1 1
EA 1 1
[k12 ] =
, [k23 ] =
1 1
L
2 L 1 1
[
]
3 EA 1 1
[k34 ] =
1 1
L
F
=?
2
2
u
=
0
1x
1

u2 F
EA
2
2
+
1/2
1/2
=
2F
u
1/2 1/2 + 3 3
L

3
u4 = 0
3
F4x =?
3
[
]{ } { }
EA 5/2 1/2 u2
F
=
u3
2F
L 1/2 7/2
do`
u:
u2 =
9 FL
17 EA
u3 =
11 F L
17 EA
33
Les actions de liaisons sont les solutions de [KP L ] {UL } = {Fnod,P } :

[
]{ } {
}
EA 2
u2
F1x
=
3 u3
F4x
L
do`
u:
F1x =
18
F
17
F4x =
33
F
17

F1x + F2x + F3x + F4x =
18
33
F + F + 2F
F =0
17
17
MEF NOR 5
Considerons la poutre daxe x representee sur la gure. Soit E le module de Young du materiau.
Laire de la section droite est egale à A entre les nuds 1 et 2, 2 A entre les nuds 2 et 3 et à 3 A
entre les nuds 3 et 4. La poutre est encastree en 1 et 4.
La poutre porte entre les nuds 1 et 2 une force uniformement repartie dintensite lineique (p, 0, 0)
avec p > 0 et au nud 3 une force (F = 2 pL, 0, 0).
Solution
Assemblage
34
Lassemblage conduit à la relation [K] {U } = {Fnod } + {F } :

F1x =?
u1 = 0
1/2
1/2
pL
EA
1/2
1/2
+
2
2
u
0
pL
2
+
=
2 pL
2
2 + 3 3
u
0
L

3
u4 = 0
F4x =?
0
3
3
Calcul des d
eplacements inconnus et des actions de liaison
Les deplacements inconnus sont les solutions de lequation [KLL ] {UL } = {Fnod,L } + {FL } :
[
]{ } {
} { }
EA 5/2 2 u2
0
pL
=
+
u3
2 pL
0
L 2 5
do`
u:
18 pL2
14 pL2
, u3 =
17 EA
17 EA
Le champ de deplacements est represente sur la gure ci-dessous.
u2 =
Les actions de liaisons sont les solutions de [KP L ] {UL } = {Fnod,P } + {FP } :
[
]{ } {
} { }
EA 1/2
u2
F1x
pL
=
+
3
u
F
0
L
3
4x
do`
u:
F1x =
26
pL ,
17
F4x =
42
pL
17

F1x + F2x + F3x + F4x + 2 pL =
26
42
pL + 0 + 2 pL
pL + 2 pL = 0
17
17
Calcul des eorts

el
ementaires et repr
esentations graphiques
Les eorts elementaires sont :
element 1-2 :
{
}
[
]{ }
{ }
{
}
EA 1 1 u1
pL 26
N1
1
=
pL
=
N2
u2
1
8
2 L 1 1
17
Le deplacement est maximal pour x = xm tel que N (xm ) = N1 p xm = 0, do`

u:
xm =
26
L = 1.52 L ,
17
u(xm ) =
676 pL2
pL2
= 1.17
578 EA
EA
35
element 2-3 :
{
}
[
]{ }
{ }
2 EA 1 1 u2
pL
N2
8
=
=
N3
1 1
u3
L
17 8
element 3-4 :
N3
N4
[
]{ }
{
}
3 EA 1 1 u3
pL
8
=
=
1 1
u4
L
17 42
MEF NOR 6
Considerons la poutre daxe x representee sur la gure. Soit E le module de Young du materiau.
Laire de la section droite est egale à A.
La poutre est encastree en 1.

La poutre est soumise à une force uniformement repartie sur toute sa longueur dintensite lineique p.
1. Calculer les deplacements nodaux et laction de liaison.
3. Representations graphiques : tracer N (x) et u(x).
Solution
Assemblage
u1 = 0 F1x =? pL/4
1
1
2 EA
0
u2
1 1 + 1 1
+ pL/2
=
L
pL/4
0
u3
1
1
Calcul des d
eplacements inconnus et de laction de liaison
36
Les deplacements inconnus sont les solutions de lequation [KLL ] {UL } = {FL } :
]{ } {
[
}
2 EA 2 1 u2
pL/2
=
1 1
u3
pL/4
L
do`
u:
3 pL2
1 pL2
, u3 =
8 EA
2 EA
La reaction dappui est solution de lequation [KP L ] {UL } = {Fnod,P } + {FP } :
u2 =
2 EA
pL
(u2 ) = F1x +
L
4
do`
u:
F1x = pL
F1x + F2x + F3x + pL = pL + 0 + 0 + pL = 0
Calcul des eorts

el
ementaires et repr
element 1-2 :
{
}
[
]{ }
{ } {
}
2 EA 1 1 u1
pL 1
N1
pL
=
=
N2
1 1
u2
pL/2
L
4 1
element 2-3 :
{
}
[
]{ }
{ } {
}
2 EA 1 1 u2
pL 1
N2
pL/2
=
=
N3
1 1
u3
0
L
4 1
MEF NOR 7
La poutre representee sur la gure est discretisee en deux elements de meme longueur L.
Soient E et A respectivement le module de Young du materiau et laire de la section droite.
37
Les sections 1 et 3 sont encastrees.

La poutre porte :
entre les nuds 1 et 2 une force uniformement repartie dintensite lineique 2 p.
entre les nuds 2 et 3 une force uniformement repartie dintensite lineique p.
1. Calculer le deplacement du nud 2 et les actions de liaison.
Solution
Assemblage

1
1
u1 = 0 F1x =? pL
EA
1 1 + 1 1
u2
0
=
+ 3 pL/2
L
1
1
u3 = 0
F3x =?
pL/2
Calcul du d
eplacement inconnu et des actions de liaison
Le deplacement inconnu est solution de lequation [KLL ] {UL } = {FL } :
EA
(2 u2 ) = 3 pL/2
L
do`
u:
u2 =
3 pL2
4 EA
Les reactions dappui sont les solutions de lequation [KP L ] {UL } = {Fnod,P } + {FP } :
[ ]
{
} {
}
EA 1 { }
F1x
pL
u2 =
+
F3x
pL/2
L 1
do`
u:
7
F1x = pL
4
5
F3x = pL
4

7
5
F1x + F2x + F3x + 3 pL = pL + 0 pL + 3 pL = 0
4
4
Calcul des eorts

el
ementaires et repr
38

element 1-2 :
{
[
]{ }
{ }
}
{ }
pL 7
EA 1 1 u1
N1
1
=
pL
=
N2
u2
1
L 1 1
4 1
Le deplacement est maximal pour x = xm tel que N (xm ) = N1 2 p xm = 0, do`

u:
xm =
element 2-3 :
7
L ,
8
u(xm ) =
49 pL2
pL2
= 0.766
64 EA
EA
}
{
[
{ }
{ }
]{ }
EA 1 1 u2
pL 1
pL 1
N2
=
=
N3
u3
L 1 1
2 1
4 5
MEF NOR 8
La poutre representee sur la gure est discretisee en deux elements. Soient E le module de Young du
materiau et A laire de la section droite. La poutre est encastree en 1 et 3.
La poutre est soumise entre les nuds 2 et 3 à une force repartie dont lintensite lineique varie entre
les valeurs 0 et p.
Solution
Assemblage
u1 = 0 F1x =?
0
1
1
pL
EA
0
u2
1
1 1 + 1 1
+
=
L
6
F3x =?
u3 = 0
2
1
1
39
Calcul du d
eplacement inconnu et des actions de liaison
Le deplacement inconnu est solution de lequation [KLL ] {UL } = {FL } :
EA
pL
(2 u2 ) =
L
6
do`
u:
1 pL2
12 EA
Les reactions dappui sont les solutions de lequation [KP L ] {UL } = {Fnod,P } + {FP } :
}
[ ]
{
{ }
pL 0
EA 1 { }
F1x
u2 =
+
F3x
2
L 1
3
u2 =
do`
u:
F1x =
1
pL
12
F3x =
5
pL
12

F1x + F2x + F3x +
1
5
1
1
pL = pL + 0
pL + pL = 0
2
12
12
2
Calcul des eorts

el
ementaires et repr
element 1-2 :
{
}
[
]{ }
{ }
EA 1 1 u1
pL 1
N1
=
=
N2
u2
1
L 1 1
12
element 2-3 :
{
}
[
]{ }
{ }
{ }
EA 1 1 u2
pL 1
pL 1
N2
=
=
N3
u3
L 1 1
6 2
12 5
Le deplacement est maximal pour x = xm tel que N (xm ) = N2 p x2m /(2 L) = 0, do`
u:
2 + 3 6 pL2
pL2
= 0.106
xm = L/ 6 , u(xm ) =
EA
36 6 EA
40
MEF NOR 9
La poutre representee sur la gure est encastree à ses deux extremites. E et A sont respectivement le
module de Young du materiau et laire de la section droite.
La poutre porte une force uniformement repartie dintensite lineique p entre les nuds 1 et 2, 2 p
entre les nuds 2 et 3 et 3 p entre les nuds 3 et 4 avec p > 0.
Solution
Assemblage
Lassemblage conduit à
1
EA
1
la relation [K] {U } = {Fnod } + {F } :
pL/2
F1x =?
u1 = 0
1
u2 =? F2x = 0 pL/2 + pL
1 + 1 1
+
=
u =?
1 1 + 1 1
pL + 3 pL/2
F3x = 0

3
3 pL/2
F4x =?
u4 = 0
1
1
Calcul des d
[
]{ } {
}
EA 2 1 u2
3 pL/2
=
u3
5 pL/2
L 1 2
do`
u:
11 pL2
13 pL2
, u3 =
6 EA
6 EA
[
]{ } {
} {
}
EA 1
u2
F1x
pL/2
=
+
1 u3
F4x
3 pL/2
L
u2 =
do`
u:
7
F1x = pL
3
F4x =
11
pL
3
11
7
pL + 6 pL = 0
F1x + F2x + F3x + F4x + 6 pL = pL + 0
3
3
41
Calcul des eorts

el
ementaires et repr
element 1-2 :
{
}
[
]{ }
{ }
{ }
EA 1 1 u1
pL 1
pL 7
N1
=
=
N2
u2
4
L 1 1
2 1
3
element 2-3 :
}
{
{ }
[
]{ }
{ }
EA 1 1 u2
pL 4
N2
1
=
pL
=
N3
u3
1
L 1 1
3 2
Le deplacement est maximal pour x = xm tel que N (xm ) = N2 2 p xm = 0, do`

u:
xm = 2 L/3
element 3-4 :
u(xm ) =
41 pL2
18 EA
}
{
[
]{ }
{ }
{
}
EA 1 1 u3
3 pL 1
pL
N3
2
=
=
N4
u4
1
L 1 1
2
3 11
MEF NOR 10 : poutre soumise `

a une variation de temp
erature
La poutre representee sur la gure est encastree à ses deux extremites et a une section droite constante
dont laire est A et le moment quadratique par rapport à z, Iz .
Soient E et repectivement le module de Young et le coecient de dilatation du materiau.

La poutre est soumise à une variation de temperature T > 0.
1. Calculer leort normal et la contrainte normale xx dans la poutre.
42

2. Application numerique. On donne :
L = 1 m, E = 200000 MPa, = 105 K1 , T = 100 K.
La poutre est un rond plein de diamètre D.
Calculer la contrainte normale xx et leort normal N .
Sachant que la force critique dEuler pour une poutre bi-encastree est :
FEuler =
4 2 EIz
L2
Calculer D.
Solution
Leort normal dans la poutre est :
N = EA T
do`
u la containte normale
xx = E T = 200000 105 100 = 200 MPa
Lequilibre de la poutre est stable si :
N = EA T <
avec :
A=
soit :
D2 >
D2
4
4 2 EIz
L2
Iz =
D4
64
4 L2 T
4 105 10002 100
=
2
2
do`
u:
D > 20.13 mm
Remarque : lallongement unitaire dans les directions et k est egal à :
yy = zz =
xx + T = (1 + ) T
E
Si = 0.3, on obtient :
yy = zz = (1 + 0.3) 105 100 = 1300 106
43
MEF NOR 11 : poutre soumise `

a une variation de temp
erature
La poutre representee sur la gure est constituee dun element de longueur L.
A est laire de la section droite.

E et sont respectivement le module de Young et le coecient de dilatation du materiau.
La section 1 est encastree. La section 2 repose sur un appui elastique de raideur k. Lorsque la structure
nest pas chargee, le ressort est au repos.
La poutre est soumise à une variation de temperature T > 0.
1. Eectuer letude elementaire.
3. Calculer le deplacement du nud 2 et leort normal N dans la poutre.
4. Representer graphiquement lallure de u2 et N en fonction de la raideur k du ressort.
Solution
Lassemblage conduit à la relation [K] {U } = {Fnod } + {Fth } :
[
]{
} {
}
{ }
EA 1 1 u1 = 0
F1x =?
1
=
+ EA T
1
1
u
=?
F
=
k
u
1
L
2
2x
2
u2 est solution de lequation :
)
EA
+ k u2 = EA T
L
do`
u:
u2 = L T
1
kL
1 + EA
et :
N = F2x = k u2 = EA T
k
k + EA
L
44
2.3
2.3.1
Treillis plans `
a nuds articul
es
Rappels
Soit (i j) un element de treillis plan de section droite

constante. i est le noeud origine et j est le noeud extremite.
Soient L la longueur de lelement et A laire de la section
droite.
Soient E et respectivement le module de Young et le
coecient de dilatation du materiau.
Lelement est soumis à une variation de temperature egale à T constante dans lelement.
(ui , vi ) et (uj , vj ) sont les deplacements nodaux.
Soit n de composantes (nx , ny ) le vecteur unitaire (n2x + n2y = 1) porte par laxe de la poutre :
nx =
xj xi
= cos
L
ny =
yj yi
= sin
L
L2 = (xj xi )2 + (yj yi )2
(2.3.1)
Lallongement unitaire suivant n est :

1
(nx (uj ui ) + ny (vj vi ))
L
(2.3.2)
N
= E (n T )
A
(2.3.3)
EA
( nx (uj ui ) + ny (vj vi ) ) EA T
L
(2.3.4)
n =
De la relation de comportement :
n =
on deduit lexpression de leort normal :
N=
puis la relation :
{fnod } = [ k ] {u} {fth }

o`
u
nx
ny
{fnod } = N
n
x
ny
nx ny
n2x nx ny
n2x
EA
nx ny
n2y
n2y
nx ny
[k] =
2
2
nx ny
nx
L nx nx ny
nx ny
n2y
nx ny
n2y

ui

vi
{u} =
u
j
vj
ny
{fth } = E A T
n
x
ny
(2.3.5a)
(2.3.5b)
(2.3.5c)
La force critique dEuler est egale à :

2 E Iz
(2.3.6)
L2
o`
u Iz est le moment quadratique de la section droite par rapport à laxe Gz (G est le centre de gravite
de la section).
FEuler =
2.3.2
45
Exercices
MEF TRE 1
Le treillis plan à nuds articules represente sur la gure est compose de trois poutres de meme nature
(module de Young E, limite elastique E ) et de meme section droite (aire A).
Les nuds 1 , 3 et 4 sont articules. Le nud 2 porte une force de composantes (P, 3P, 0) avec P > 0.
1.
2.
3.
4.
Calculer les deplacements des nuds.

Calculer les actions de liaison et verier lequilibre de la structure.
Calculer les eorts normaux dans les poutres.
On donne : L = 0.8 m , P = 25000 N , E =200000 MPa , E = 300 MPa.
Les poutres sont des carres pleins de cote c. Determiner c pour que les critères de traction,
compression et ambement (critère dEuler) soient satisfaits.
Solution
Etude
el
ementaire et assemblage
Les caracteristiques elementaires sont :
element 1 2 :
L ,
element 2 3 :
2L ,
nx = 1
ny = 0
element 4 2 :
3L ,
nx = 0
ny = 1
do`
u les matrices de rigidite elementaires :
1 0 1
EA
0 0 0
[k12 ] =
L 1 0 1
0 0 0
nx = 1
0
0
0
0
ny = 0
EA 0
[k42 ] =
3 L 0
0
1
EA
0
[k23 ] =
2 L 1
0
0 0 0
1 0 1
0 0 0
1 0 1
0 1 0
0 0 0
0 1 0
0 0 0
46
1
0
EA
0
L
0
0
0
0
F1x =?
u1 = 0
1
0
F
=?
v
=
0
0
0
1y
1
u
=?
F
=
P
3/2
0
1/2 0 0
0
2x
2
F
=
3
P
0
1/3
0
0 0 1/3
v
=?
2y
2
=
u3 = 0 F3x =?
1/2
0
1/2 0
v3 = 0
F3y =?
0
0
0
0
F
=?
0
0
0
0
u
=
0
4x
4
F4y =?
0
1/3
0 1/3
v4 = 0
Calcul des d
[
]{ } {
}
EA 3/2 0
u2
P
=
3 P
0 1/3 v2
L
do`
u:
2 PL
PL
, v2 = 9
3 EA
EA
Les actions de liaison sont les solutions de [KP L ] {UL } = {Fnod,P } :
u2 =
2 P/3
F1x
1
0
0
F
0
1y
{ }
EA 1/2
P/3
F3x
0 u2
=
=
0
F3y
0
0
v2
0
F4x
0
3P
F4y
0
1/3
Lequilibre de la structure est verie :

F1x + F2x + F3x + F4x = 2 P/3 + P P/3 + 0 = 0
F1y + F2y + F3y + F4y = 0 3 P + 0 + 3 P = 0
Calcul des eorts normaux

Les eorts normaux sont :
N12 =
EA
(1 (u3 u2 ) + 0 (v3 v2 )) = P/3 < 0 (compression)
2L
EA
=
(0 (u2 u4 ) + 1 (v2 v4 )) = 3 P < 0 (compression)
3L
N23 =
N42
Dimensionnement
EA
(1 (u2 u1 ) + 0 (v2 v1 )) = 2 P/3 > 0 (traction)
L
47
Le dimensionnement en contrainte secrit :

max =
soit :
A = c2 >
3P
< E
A
3P
3 250000
=
E
300
do`
u:
c > 50 mm
La condition de stabilite des elements secrit :
|N23 | =
P
2 EIz
<
3
4 L2
soit :
Iz =
On en deduit :
c4 >
|N42 | = 3 P <
2 EIz
9 L2
c4
27 P L2
>
12
2E
12 27 P L2
12 27 250000 8002
=
2E
2 200000
do`
u:
c > 71.59 mm
MEF TRE 2
Le treillis plan à nuds articules represente sur la gure est compose de deux poutres de meme section
droite (aire A) et de meme materiau (module de Young E, limite elastique E ).
Les nuds 1 et 3 sont articules. Le nud 2 porte une force de composantes (0, P, 0) avec P > 0.
1. Calculer les deplacements des nuds.
3. Calculer les eorts normaux dans les poutres.
4. On donne : L = 0.2 m , P = 10 kN , E = 200000 MPa , E = 300 MPa.
Les poutres sont des carres pleins de cote c. Determiner c pour que les critères de traction,
compression et ambement (critère dEuler) soient satisfaits.
Solution
Etude
el
48

element 1-2 :
L 2 ,
element 2-3 :
E
,
,
A
A
,
nx = ny = 1/ 2
nx = 1 ,
do`
1
1 1 1
EA 1
1 1 1
[k12 ] =
1
2 2L 1 1 1
1 1 1
1
ny = 0
1
1
1
1
1
1
1
1
EA 1 1 2 2 + 1 1 2 2
1
1
0
2 2L
1 1
2 2
0
2 2
0
0
0
1
EA
0
[k23 ] =
L 1
0
0 1 0
0 0 0
0 1 0
0 0 0
u1 = 0
F1x =?
v1 = 0
F1y =?
0 u2 =?
F2x = 0
=
0
v2 =?
F2y = P
u
=
0
F3x =?
0 3
v3 = 0
F3y =?
0
Calcul des d
[
]{ } {
}
EA 2 2 + 1 1 u2
0
=
P
1
1 v2
2 2L
do`
u:
PL
PL
, v2 = (1 + 2 2)
EA
EA
1
1 { }
F1x
EA
1
1
u2
F1y
F3x
2 2 L 2 2 0 v2
F3y
0
0
u2 =
do`
u:
F1x = F1y = P
F3x = P
F3y = 0

F1x + F2x + F3x = P + 0 P = 0 ,
F1y + F2y + F3y = P P + 0 = 0

N12
EA
=
L 2
N23 =
1
1
(u2 u1 ) + (v2 v1 ) = 2 P < 0
2
2
(compression)
EA
(1 (u3 u2 ) 0 (v3 v2 )) = P < 0 (compression)
L
49
Dimensionnement
max
soit :
2P
=
< E
A
2P
2 10000
A=c >
=
E
300
2
do`
u:
c > 6.87 mm
La condition de stabilite des elements secrit :
2 EIz
2 EIz
, N23 = P <
|N12 | = 2P <
2
2L
L2
On en deduit :
c4
2 2 P L2
2 2 10000 2002
Iz =
>
=
12
2E
2 200000
do`
u:
c > 9.11 mm
MEF TRE 3
Le treillis plan represente sur la gure est compose de deux poutres articulees entre elles. Les nuds 1
et 3 sont articules. Soient E et E respectivement le module de Young et la limite elastique du
materiau. Laire des sections droites est egale à A.
Le nud 2 porte une force de composantes (P, 2P, 0) avec P > 0.

3. Calculer les eorts normaux dans les poutres.
Solution
Etude
el
50

element 1-3 :
element 2-3 :
L 2 ,
nx = 1 ,
do`
u les matrices de rigidite elementaires
1 0 1
EA
0 0 0
[k13 ] =
L 1 0 1
0 0 0
ny = 0
nx = 1/ 2
ny = 1/ 2
0
0
0
0
1 1 1 1
EA 1 1
1 1
[k23 ] =
1 1
2 2L 1 1
1 1 1 1

2 2 0
2 2
0
u1 = 0
F1x =?
0
v1 = 0
0
0
0
F1y =?
EA
1 1
1
1 u2 = 0
F2x =?
=
1 1
1
F2y =?
2 2L

v2 = 0
2 2 0 1 1 2 2 + 1 1
u
=?
F
3
3x = P
v3 =?
F3y = 3 P
0
0 1 1
1
1
Calcul des d
[
]{ } {
}
EA 2 2 + 1 1 u3
P
=
v3
3 P
1
1
2 2L
do`
u:
PL
PL
, v3 = (4 + 6 2)
EA
EA
F1x
2 2 0 { }
EA
u3
F1y
0
0
=
F
1 v3
2 2 L 1
2x
F2y
1
1
u3 = 4
do`
u:
F1x = 4 P
F1y = 0
F2x = 3 P
F2y = 3 P

F1x + F2x + F3x = 4 P P 3 P = 0
F1y + F2y + F3y = 0 + 3 P 3 P = 0

EA
(1 (u3 u1 ) + 0 (v3 v1 )) = 4 P < 0 (compression)
L
(
)
1
1
EA
(u3 u2 ) (v3 v2 ) = 3 2 P > 0 (traction)
=
L 2
2
2
N13 =
N23
51
MEF TRE 4
Le treillis plan represente sur la gure est compose de trois poutres. La structure est en acier de
module de Young E.
Laire des sections droites est A (poutres 3-1 et 1-2) et 2 A (poutre 3-2).
Le nud 1 est articule et le nud 3 repose sur un appui simple dont la normale est horizontale. Le
nud 2 porte une force de composantes (P, 2 P, 0) avec P > 0.
3. Calculer les tensions dans les poutres.
4. On donne : P = 25000 N , L = 0.6 m , E = 200000 MPa , E = 300 MPa
Les poutres
3-1 et 1-2 sont des ronds pleins de diamètre D ; la poutre 3-2 est un rond plein de
diamètre 2 D. Calculer D.
Solution
Etude
el
element 1-2 :
L ,
nx = 1
ny = 0
element 3-1 :
L ,
nx = 0
ny = 1
element 3-2 :
L 2 ,
2A
do`
1 0 1
EA
0 0 0
[k12 ] =
L 1 0 1
0 0 0
nx = ny = 1/ 2
0
0
0
0
0 0
EA
0 1
[k31 ] =
L 0 0
0 1
1
1 1 1
EA
1 1 1
1
[k32 ] =
1
2 L 1 1 1
1 1 1
1
0 0
0 1
0 0
0 1
52
2
2
0
0
0
F
=?
u
=
0
1x
1
F
=?
v
=
0
2
0
0
0
1y
1
EA 2
F2x = P
0
2 + 1 1 1
1 u2 =?
=
F2y = 2 P
v2 =?
0
1
1 1
1
2L
F
=?
u
=
0
0
1
1
1
1
3x
3
F3y = 0
v3 =?
0
2
1
1 1
2+1
Calcul des d
2+1 1
1
u2 P
EA
1
1 1 v2 = 2 P

2L
v3
0
2+1
1
1
do`
u:
PL
PL
PL
, v2 = (5 + 2 2)
, v3 = 2
EA
EA
EA
u2 = 3

0
2 0
u2 F1x
EA
0
0 2 v2 = F1y

2 L 1 1
v3
F3x
1
do`
u:
F1x = 3 P
F1y = 2 P
F3x = 2 P

F1x + F2x + F3x = 3 P + P + 2 P = 0
F1y + F2y + F3y = 2 P 2 P + 0 = 0

N12 =
N32
EA
(1 (u2 u1 ) + 0 (v2 v1 )) = 3 P > 0 (traction)
L
EA
(0 (u3 u1 ) + 1 (v3 v1 )) = 2 P > 0 (traction)
N31 =
L
(
)
2 EA
1
1
(u2 u3 ) + (v2 v3 ) = 2 2 P < 0 (compression)
=
L 2
2
2
Les contraintes normales sont :

12
3P
=
A
31
2P
=
A
32
2P
=
A
53
Dimensionnement
max =
soit :
A=
do`
u:
D2 >
3P
< E
A
D2
3P
>
4
E
12 P
12 25000
=
E
300
soit :
D > 17.85 mm
La condition de stabilite de lelement 3 2 secrit :
|N32 | = 2
soit :
Iz32
On en deduit :
4
D >
64
2 EIz32
2P <
2 L2
( 2D)4
4 2 P L2
=
>
64
2E
2 P L2
64 2 25000 6002
=
3E
3 200000
do`
u:
D > 19.03 mm
54
2.4
2.4.1
Poutre soumise `
a un moment de torsion
Rappels
Considerons lelement de poutre droite (i j), de longueur L et de section droite constante, soumis
à un moment de torsion Mt (x).
J est la constante de torsion de Saint-Venant et G =
E
est le module delasticite transversal.
2 (1 + )
Lelement porte sur toute sa longueur un couple lineairement reparti dintensite lineique (mx , 0, 0)
avec :
(
x
x)
x
mx x) = mxi 1
+ mxj = mxi + (mxj mxi )
(2.4.1)
L
L
L
xi et xi sont les deplacements nodaux.
Mti et Mtj sont les moments de torsion dans les sections i et j.
Lintegration de lequation dequilibre :
dMt
+ mx = 0 avec Mt (0) = M ti
dx
conduit à lexpression du moment de torsion :
Mt (x) = Mti mxi x ( mxj mxi )
x2
2L
(2.4.2)
(2.4.3)
Lintegration de la relation de comportement :

Mt = GJ
donne la deformee :
1
x (x) = xi +
GJ
dx
dx
avec x (0) = xi
(
)
x2
x3
Mti x mxi
( mxj mxi )
2
6L
(2.4.4)
(2.4.5)
Des conditions aux limites :

Mt (L) = Mtj
x (L) = xj
on deduit :
{fnod } = [ k ] {u} {f }
avec :
} {
}
{ }
Mt (0)
Mti
xi
{fnod } =
=
, {u} =
Mt (L)
Mtj
xj
[
]
{
}
GJ 1 1
L 2 mxi + mxj
[k] =
, {f } =
L 1 1
6 mxi + 2 mxj
(2.4.6a)
(2.4.6b)
(2.4.6c)

Mti + Mtj +
L
(mxi + mxj ) = 0
2
(2.4.7)
2.4.2
55
Exercices
MEF TOR 1
Considerons larbre (1, 2, 3, 4) represente sur la gure ci-contre. Soit G le module delasticite de glissement du materiau. La rigidite de torsion est egale à GJ entre les nuds 1 et 2, 3 GJ entre les nuds 2
et 3 et à 2 GJ entre les nuds 3 et 4.
La poutre est encastree en 1 et 4. Elle est sollicitee au nud 2 par un couple (C, 0, 0) avec C > 0.
3. Representations graphiques : tracer Mt (x) et x (x).
Solution.
[
]
[
]
GJ 1 1
3 GJ 1 1
[k12 ] =
, [k23 ] =
L 1 1
2 L 1 1
[k34 ] =
[
]
2 GJ 1 1
1 1
L
M1x =?
x1 = 0
1
1
x2 =? C
GJ
1
1
+
3/2
3/2
=
x3 =?
3/2 3/2 + 2 2
0
L

M4x =?
2
2
x4 = 0
[
]{ } { }
GJ 5/2 3/2 x2
C
=
x3
0
L 3/2 7/2
do`
u:
7 CL
3 CL
, x3 =
13 GJ
13 GJ
Les actions de liaisons sont les solutions de [KP L ] {UL } = {Fnod,P } :
[
]{ } {
}
GJ 1
x2
M1x
=
2 x3
M4x
L
x2 =
do`
u:
M1x =
7
C
13
M4x =
6
C
13

M1x + M2x + M3x + M4x =
6
7
C +C +0
C=0
13
13
56
Le moment de torsion Mt (x) et la rotation des sections droites x (x) sont representes sur la gure
ci-dessous.
MEF TOR 2
La poutre de longueur 3 L et de section constante representee sur la gure ci-contre est encastree
à ses deux extremites. Soient J la constante de torsion de Saint-Venant et G le module delasticite
transversal.
Elle porte entre les nuds 1 et 2 un couple de torsion uniformement reparti dintensite lineique m.
Le nud 2 est soumis à un couple de torsion dintensite mL.
3. Representations graphiques : tracer Mt (x) et x (x).
Solution.
Les matrices elementaires sont :
[
]
GJ 1 1
[k12 ] =
2 L 1 1
[
]
GJ 1 1
[k23 ] =
L 1 1
{ }
1
{f12 } = mL
1

1/2
1/2
x1 = 0 M1x =? mL
GJ
mL
1/2 1/2 + 1 1 x2 =? =
+ mL
L
x3 = 0
M3x =?
0
1
1
Le deplacement inconnu est solution de lequation lequation [KLL ] {UL } = {Fnod,L } + {FL } :
GJ 3
x2 = 2 mL
L 2
do`
u:
x2 =
4 mL2
3 GJ
57
[
]
{
} {
}
GJ 1/2 { }
M1x
mL
x2 =
+
1
M3x
0
L
do`
u:
5
M1x = mL
3
4
M3x = mL
3
5
4
M1x + M2x + M3x + 2 mL = mL + mL mL + 2 mL = 0
3
3
element 1-2 :
}
{
[
]{ }
{ }
{ }
GJ 1 1 x1
mL 5
Mt1
1
=
mL
=
Mt2
x2
1
1
2 L 1 1
3
Le deplacement est maximal pour x = xm tel que Mt (xm ) = Mt1 m xm = 0, do`

u:
xm =
element 2-3 :
5
L ,
3
x (xm ) =
25 mL2
mL2
= 1.39
18 GJ
GJ
{
}
[
]{ }
{ }
GJ 1 1 x2
mL 4
Mt2
=
=
Mt3
x3
4
L 1 1
3
Le moment de torsion Mt (x) et la rotation des sections droites x (x) sont representes sur la gure
ci-dessous.
58
2.5
2.5.1
Flexion des poutres `

a plan moyen : mod`
ele de Bernoulli
Rappels : exion dans le plan {xy}
Pour chaque poutre etudiee :

laxe x est la bre moyenne de la poutre.
le plan {x, y} est un plan de symetrie des sections de la poutre.
laxe z forme avec les axes x et y un trièdre direct.
Lelement de poutre (i j) represente ci-contre, de longueur L, de moment quadratique Iz constant et

de module de Young E est soumis sur toute sa longueur à une force uniformement repartie dintensite
lineique py .
EIz est la rigidite de exion lineique de la poutre.
Ty et Mfz sont leort tranchant et le moment echissant.
v et z sont respectivement le deplacement suivant y (èche) et la rotation autour de laxe z dune
section droite.
(Tyi , Mfzi ) et (Tyj , Mfzj ) sont les eorts nodaux.
(vi , zi ) et (vj , zj ) sont les deplacements nodaux.
Lintegration des deux equations dequilibre :
dTy
+ py = 0
dx
avec
Ty (0) = Tyi
(2.5.1)
et
dMfz
+ Ty = 0 avec Mfz (0) = Mfzi
dx
conduit à lexpression de leort tranchant et du moment echissant :
Ty (x) = Tyi py x
x2
Mfz (x) = Mfzi Tyi x + py
2
Remarque : si py <> 0, leort tranchant sannule pour xm deni par :
Ty (xm ) = 0
soit xm =
Tyi
py
Si 0 < xm < L, le moment echissant passe par une valeur extremale :

Mfz (xm ) = Mfzi
2
Tyi
2 py
(2.5.2)
(2.5.3)
59
La contrainte normale dans la poutre est egale à :

xx (x, y) =
La quantite Wel.z =
y
Mfz (x)
Iz
(2.5.4)
Iz
est le module de exion elastique.
|y|max
La relation de comportement secrit :

EIz
dz
= Mfz
dx
avec
z (0) = zi
(2.5.5)
Au cours de la mise en charge, une section droite reste plane et perpendiculaire à la deformee de la
bre moyenne (hypoth`
ese de Bernoulli) do`
u la relation cinematique :
z =
dv
dx
avec
v(0) = vi
(2.5.6)
La deformee et la rotation des sections droites sont donnees par les formules :
(
)
1
x2
x3
z (x) = zi + EI Mfzi x Tyi 2 + py 6

z
(
)
x2
x3
x4
v(x) = vi + zi x + 1
Mfzi
Tyi
+ py
EIz
2
6
24
(2.5.7)

Ty (L) = Tyj
Mfz (L) = Mfzj
v(L) = vj
z (L) = zj
(2.5.8)
on deduit lexpression des eorts nodaux en fonction des deplacements nodaux :

{fnod } = [ k ] {u} {f }
(2.5.9a)

vi

zi
, {u} =
v
zj
p L/2
y
2
py L /12
py L
L
=
{f } =
py L/2
6
12
L
2
p L /12
Tyi
Ty (0)
Mfzi
Mfz (0)
=
{fnod } =
T (L)
Tyj
y

Mfzj
Mfz (L)
avec
12
6 L 12 6 L
EIz 6 L 4 L2 6 L 2 L2
[k] = 3
L 12 6 L 12 6 L
6 L 2 L2 6 L 4 L2
(2.5.9b)
(2.5.9c)
Tyi + Tyj + py L = 0
Mfzi + Mfzj + L Tyj + 1 py L2 = 0
2
(2.5.10)
60
2.5.2
Exercices
MEF FLE 1 : exion dans le plan {x, z}

Soit une poutre droite sollicitee en exion simple dans le plan {x, z}. Cette poutre porte une force repartie dintensite lineique (0, 0, pz ) et un couple reparti dintensite lineique (0, my , 0). Les eorts resultants dans la section droite dabscisse x sont leort tranchant Tz (x) et le moment echissant Mfy (x).
1. Ecrire
lequilibre du troncon de poutre compris entre les abscisses x et x + dx.
2. En deduire les equations dierentielles que doivent satisfaire Tz et Mfy .
Solution
Lequilibre du morceau de poutre compris entre les sections droites dabscisses x et x + dx secrit au
premier ordre près :
dTz
dx + pz dx = 0
Tz + Tz +
dx
dMfy
Mfy + Mfy +
dx dx Tz + my dx = 0
dx
Après simplication, on obtient les deux equations dequilibre :
dTz
+ pz = 0
dx
dMfy
Tz + my = 0
dx
MEF FLE 2
La poutre droite representee sur la gure est en acier de module de Young E et de limite elastique E .
Elle a une section constante de moment quadratique Iz .
La poutre est encastree en 1 et repose sur un appui simple en 2 et 3.

La section 3 porte un couple dintensite (0, 0, C) avec C > 0.
61

3. Tracer le diagramme des eorts tranchants et le diagramme des moments echissants. Determiner la position et la valeur du moment echissant maximal.
4. Application numerique : E = 300 MPa , C = 1000 daN.m. La poutre est un rond plein de
diamètre D. Calculer D.
Solution
Assemblage
Lassemblage conduit à la relation [K] {U } = {Fnod }
12
6 L 12 6 L
0
6 L 4 L2 6 L 2 L2
0
EIz
12
6
L
24
0
12
2
0
8 L2 6 L
L3
6L 2L
0
0
12 6 L 12
0
0
6 L 2 L2 6 L
F
=?
0
v
=
0
1y
1
M
=?
0
=
0
1z
z1
F
=?
6L
v
=
0
2y
2
=
M2z = 0
2 L2
=?
z2
F
=?
6 L
v
=
0
3y
3

2
M3z = C
4L
z3 =?
Calcul des d
eplacements inconnus et des actions de liaisons
[
]{ } { }
EIz 8 L2 2 L2 z2
0
=
2 4 L2
3
2
L
C
L
z3
do`
u:
2 CL
1 CL
, z3 =
14 EIz
7 EIz
F
6L
0
1y
{ }
2
z2
EIz
M
2
L
0
1z
=
F
6 L z3
L3 0
2y
6 L 6 L
F3y
z2 =
do`
u:
3C
1
12 C
9C
, M1z = C , F2y =
, F3y =
7 L
7
7 L
7 L
(
)
C 3 12 9
F
+
F
+
F
=
+
=0
1y
2y
3y
L
7
7
7
(
)
1
12 18
=0
M1z + M2z + M3z + LF2y + 2 LF3y = C + 0 + 1 +
7
7
7
F1y =
Calcul des eorts

el
ementaires et repr
62
T
6
L
3
C/L
= F1y
y1
EI
2
1
Mfz1
C
= M1z
z 2L
{z2 } =
= 3
T
6
L
3
C/L
L
7
y2
Mfz2
4 L2
2 C
element 1-2 :
element 2-3 :
6L
6L { }
Ty2
9 C/L
EI
2
1 2C
2 L2
Mfz2
z2
z 4L
= 3
=
6 L 6 L z3
= F3y
Ty3
9 C/L
L
7
2
2
Mfz3
2L
4L
7C
Leort tranchant Ty (x) et le moment echissant Mfz (x) sont representes sur la gure ci-dessous.
Dimensionnement
Le moment echissant est maximal dans la section 3 :
max(|Mfz (x)|) = C
Le dimensionnement en contrainte secrit donc :
Wel.z =
do`
u:
D3 >
C
D3
>
32
E
32 107
32 C
=
E
300
D > 69.76 mm
MEF FLE 3
La poutre droite representee sur la gure est encastree en 1 et repose sur un appui simple en 2 et 3.
Soient E le module de Young du materiau et Iz le moment quadratique de la section droite.
La poutre porte en 2 et 3 un couple ponctuel de composantes (0, 0, C) avec C > 0.

63
3. Tracer le diagramme des eorts tranchants et le diagramme des moments echissants.

4. Application numerique : E = 300 MPa , C = 10000 N.m. Calculer Wel.z .
Solution
Assemblage
Lassemblage conduit à la relation [K] {U } = {Fnod }
12
6 L 12 6 L
0
6 L 4 L2 6 L 2 L2
0
EIz
0
12
12 6 L 24
2
2
6
L
2
L
0
8
L
6
L
L3
0
0
12 6 L 12
0
0
6 L 2 L2 6 L
0
v1 = 0
F1y =?
z1 = 0
0
M1z =?
6 L v2 = 0
F2y =?
=
2 L2
=?
M2z = C
z2
6 L
v
=
0
F3y =?
3
2
4L
z3 =?
M3z = C
Calcul des d
[
]{ } { }
EIz 8 L2 2 L2 z2
C
=
2
2
3
2L 4L
z3
C
L
do`
u:
1 CL
3 CL
, z3 =
14 EIz
14 EIz
F1y
6L
0
{
}
2
z2
EIz
M
2
L
0
1z
=
F
6 L z3
L3 0
2y
6 L 6 L
F3y
z2 =
do`
u:
3C
1
9C
12 C
, M1z = C , F2y =
, F3y =
7 L
7
7 L
7 L
(
)
C 3 9 12
+
=0
F1y + F2y + F3y =
L 7 7
7
(
)
1
9 24
+1+1+
=0
M1z + M2z + M3z + LF2y + 2 LF3y = C
7
7
7
F1y =
Calcul des eorts

el
ementaires et repr
element 1-2 :
= F1y
3 C/7 L
Ty1
6L
EI
2
= M1z
C/7
Mfz1
2
L
z
{z2 } =
= 3
3 C/7 L
6 L
Ty2
2 C/7
Mfz2
4 L2
64

element 2-3 :
12 C/7 L
Ty2
6L
6L { }
EI
2
5 C/7
z2
Mfz2
2 L2
z 4L
=
= 3
12 C/7 L
= F3y
6 L 6 L z3
Ty3
L
C
Mfz3
2 L2 4 L2
Dimensionnement
Le moment echissant est maximal dans la section 3 :
max(|Mfz (x)|) = C
Wel.z >
C
10000 103
=
= 33333 mm3 34 cm3
E
300
MEF FLE 4 : d
eplacement impos
e
La poutre droite representee sur la gure est encastree en 1. Soient E le module de Young du materiau
et Iz le moment quadratique de la section droite.
Le nud 2 subit un deplacement vertical egal à v2 = d avec d > 0.

Calculer la pente en 2 et les actions de liaison.
Solution
F1y =?
v1 = 0
12
6 L 12 6 L

2
2
EIz
M
=?
=
0
6
L
4
L
6
L
2
L
1z
z1
=
v = d
F2y =?
L3 12 6 L 12 6 L

2
M2z = 0
z2 =?
6 L 2 L2 6 L 4 L2
65
La rotation z2 est solution de lequation [KLL ] {UL } = [KLP ] {UP } :

{ }
] d
EIz [
2
6 L 4 L
= {0}
z2
L3
do`
u:
z2 =
3 d
2L
Les actions de liaison sont donnees par [KP L ] {UL } + [KP P ] {UP } = {Fnod,P } :
12 6 L { } F1y
EIz
d
6 L 2 L2
= M1z
z2
L3
F2y
12 6 L
do`
u:
F1y = F2y = 3
EIz d
L3
M1z = 3
EIz d
L2
MEF FLE 5 : d
eplacement impos
e
La poutre droite representee sur la gure ci-dessous est en acier de module de Young E et de limite
elastique E . Soit Iz le moment quadratique de la section droite.
La poutre est encastree en 1 et repose sur un appui simple en 2. Le nud 3 subit un deplacement
vertical : v3 = d avec d > 0.
3. Tracer les diagrammes de leort tranchant et du moment echissant. Determiner la position et
la valeur du moment echissant maximal.
4. Quelle inegalite doit satisfaire la grandeur v = |y|max (intervenant dans le module de exion
Wel.z ).
5. Application numerique :
E = 200000 MPa , E = 300 MPa , L = 90 cm , d = 6 mm
Calculer v.
Solution
Assemblage
66
F1y =?
12
6 L 12 6 L
0
0
v1 = 0
6 L 4 L2 6 L 2 L2
0
0
=
0
M
=?
1z
z1
EIz 12 6 L 24
F2y =?
0
12 6 L v2 = 0
=
2
M2z = 0
0
8 L2 6 L 2 L2
z2 =?
L3
6L 2L
F
=?
0
12 6 L 12 6 L
v
=
d
3y
3
2
2
0
0
6L 2L
6 L 4 L
z3 =?
M3z = 0
Calcul des d
Les deplacements inconnus sont les solutions de lequation [KLL ] {UL } = [KLP ] {UP } :

[ 2
] z2 { }
2
EIz 8 L 6 L 2 L
0
d =
0
L3 2 L2 6 L 4 L2
z3
do`
u:
3 d
9 d
, z3 =
7L
7L
Les actions de liaison sont les solutions de [KP L ] {UL } + [KP P ] {UP } = {Fnod,P } :
F1y
6L
0
0
z2
2
EIz
0
0
2L
d = M1z
F
12 6 L
L3 0
z3
2y
F3y
6 L 12 6 L
z2 =
do`
u:
18 EIz d
6 EIz d
30 EIz d
12 EIz d
, M1z =
, F2y =
, F3y =
7 L3
7 L2
7 L3
7 L3
(
)
EIz d
18 30 12
=0
F1y + F2y + F3y =
L3
7
7
7 (
)
6
30 24
EIz d
+0+0+
=0
M1z + M2z + M3z + LF2y + 2 LF3y =
L2
7
7
7
F1y =
Calcul des eorts

el
ementaires et repr
element 1-2 :
Ty1
6L
EI
Mfz1
2 L2
z
{z2 } = EIz d
= 3
T
6 L
L
7 L3
y2
2
Mfz2
4L
18
= F1y
= M1z
6 L
18
12 L
element 2-3 :
Ty2
6 L 12 6 L
EI
z2
Mfz2
4 L2 6 L 2 L2
z
d = EIz d
= 3
L 6 L 12 6 L
7 L3
Ty3
z3
2
2
Mfz3
2L
6 L 4 L
12
12 L
= F3y
12
67
Dimensionnement
Le moment echissant maximal est :
max(|Mfz (x)|) =
12 EIz d
7 L2

Wel.z =
do`
u:
v<
Iz
12 EIz d
>
v
7 L2 E
7 9002 300
7 L2 E
=
12 Ed
12 200000 6
v < 118.1 mm
MEF FLE 6
Soit Iz le moment quadratique de la section droite.
La poutre est encastree en 1 et repose sur un appui simple en 2 et 3. Elle porte entre les nuds 1 et 2
une force uniformement repartie. Soit (0, p, 0) avec p > 0 la force par unite de longueur.
4. Application numerique : E = 300 MPa , L = 1.7 m , p = 2000 daN/m. La poutre est un rond
plein de diamètre D. Valeur de D ?
Solution
Assemblage
68
pL/2
F
=?
v
=
0
12
6 L 12 6 L
0
0
1y
1
2 /12
6 L 4 L2 6 L 2 L2
M
=?
pL
0
0
=
0
1z
z1
v2 = 0
EIz
pL/2
F
=?
12
6
L
24
0
12
6
L
2y
+
=
2
pL2 /12
M2z = 0
z2 =?
0
8 L2 6 L 2 L2
L3
6L 2L
0
F3y =?
v3 = 0
0
12 6 L 12 6 L

2
2
M3z = 0
0
0
0
6L 2L
6 L 4 L
z3 =?
Calcul des d
[
]{ }
{ }
pL2 1
EIz 8 L2 2 L2 z2
=
L3 2 L2 4 L2 z3
12 0
do`
u:
z2 =
1 pL3
84 EIz
z3 =
1 pL3
168 EIz
Les actions de liaison sont les solutions de [KP L ] {UL } = {Fnod,P } + {FP } :
pL/2
F
6L
0
1y
{ }
2 /12
2
z2
EIz
pl
M
2
L
0
1z
+
=
pL/2
F
6 L z3
L3 0

2y
0
F3y
6 L 6 L
do`
u:
F1y =
4
pL
7
M1z =
3
pL2
28
F2y =
13
pL
28
F3y =
1
pL
28
(
)
4 13
1
F
+
F
+
F
pL
=
pL
+
1
=0
1y
2y
3y
7 28 28
(
)
pL2
3
13
2
1
= pL
+0+0+
=0
M1z + M2z + M3z + LF2y + 2 LF3y
2
28
28 28 2
Calcul des eorts

el
ementaires et repr
element 1-2 :
6
16 pL
Ty1
6L
= F1y
EI

2 = M
2
1
pL
L
3
pL
Mfz1
z 2L
1z
{ }
=
= 3
6 L z2
6
12 pL
Ty2
L
12
28
L
pL2
Mfz2
4 L2
Leort tranchant sannule pour xm = 4 L/7 do`
u Mfz (xm ) = 11 pL2 /196 0.056 pL2 .
69
element 2-3 :
T
6
L
6
L
pL
y2
{ }
EI
2
2
1
Mfz2
2 L2
pL
z 4L
z2
= 3
=
T
6 L 6 L z3
pL
= F3y
L
28
y3
2
2
Mfz3
2L
4L
0
Dimensionnement
max(|Mfz (x)|) = 3 pL2 /28
Wel.z =
do`
u:
D3 >
3 pL2
D3
>
32
28 E
24 pL2
24 20 17002
=
7 E
7 300
D > 59.47 mm
MEF FLE 7
La poutre droite representee sur la gure ci-contre est en acier de module de Young E et de limite
elastique E . Soit Iz le moment quadratique de la section droite.
La poutre est encastree en 1 et repose en 2 sur un appui simple. Elle porte sur toute sa longueur une
force uniformement repartie. Soit (0, p, 0) avec p > 0 la force par unite de longueur.
70

4. Application numerique : E = 300 MPa , L = 0.8 m , p = 2000 daN/m . La poutre est un rond
plein de diamètre D. Valeur de D ?
Solution
Assemblage
12
6L
6 L 4 L2
EIz
12 6 L
2
L3
6L 2L
0
0
0
0

12 6 L
0
0
v1 = 0
F1y =?
pL/2
z1 = 0
6 L 2 L2
0
0
M1z =?
pL2 /12
24
0
12 6 L v2 = 0
F2y =?
pL
=
+
0
8 L2 6 L 2 L2
z2 =?
0

M2z = 0

v3 =?

F3y = 0

12 6 L 12 6 L
pL/2
2
2
2
6L 2L
6 L 4 L
z3 =?
M3z = 0
pL /12
Calcul des d
Les deplacements inconnus sont les solutions
8 L2 6 L
EIz
6 L 12
L3
2 L2 6 L
de lequation [KLL ] {UL } = {FL } :

2 L2 z2 0
6 L v3 = pL/2
2
4 L2
z3
pL /12
do`
u:
11 pL4
13 pL3
5 pL3
, v3 =
, z3 =
48 EIz
48 EIz
48 EIz

6L
0
0 z2 F1y pL/2
EIz 2
2L
0
0 v3 = M1z + pL2 /12

L3
pL
F2y
0
12 6 L
z3
z2 =
do`
u:
1
17
1
F1y = pL , M1z = pL2 , F2y =
pL
8
8
8
(
)
1 17
+02 =0
F1y + F2y + F3y 2 pL = pL +
8
8
(
)
1
17
2
2
M1z + M2z + M3z + LF2y + 2 LF3y 2 pL = pL + 0 + 0 +
2 =0
8
8
Calcul des eorts

el
ementaires et repr
element 1-2 :
pL/8
6
Ty1
6L
= F1y
EI

2 /8 = M
2
pL
pL
L
Mfz1
2
L
z
1z
= 3
6 L {z2 } 12 6 = 9 pL/8
T
L
y2
pL2 /2
L
Mfz2
4 L2
71
element 2-3 :
Ty2
6 L 12 6 L
6
pL
z2 pL 2
EI
2
Mfz2
6 L 2 L2
L
pL /2
z 4L
v3
= 3
=
T
6 L 12 6 L 12
6
0
z3
y3
Mfz3
2 L2 6 L 4 L2
L
0
Dimensionnement
max(|Mfz (x)|) = pL2 /2
Wel.z =
do`
u:
D3 >
pL2
D3
>
32
2 E
32 pL2
32 20 8002
=
2 E
2 300
D > 60.12 mm
MEF FLE 8
Elle a une section constante de moment quadratique Iz . La poutre est encastree en 1 et repose sur un
appui simple en 2 et 3. Elle porte sur toute sa longueur une force uniformement repartie dintensite
lineique (0, p, 0) avec p > 0.

4. Application numerique : E = 300 MPa , L = 0.6 m , p = 1500 daN/m. Calculer le module de
exion elastique Wel.z .
72
Solution
Assemblage
12
6L
6 L 4 L2
EIz
12 6 L
2
L3
6L 2L
0
0
0
0

12 6 L
0
0
v1 = 0
F1y =?
pL/2
2 /12
6 L 2 L2
0
0
=
0
M
=?
pL
z1
1z
24
0
12 6 L v2 = 0
F2y =?
pL
=
+
0
8 L2 6 L 2 L2
z2 =?
M2z = 0
0
v3 = 0

F3y =?

12 6 L 12 6 L
pL/2
6 L 2 L2 6 L 4 L2
z3 =?
M3z = 0
pL2 /12
Calcul des d
{ }
[
]{ }
EIz 8 L2 2 L2 z2
pL2 0
=
L3 2 L2 4 L2 z3
12 1
do`
u:
1 pL3
1 pL3
, z3 =
168 EIz
42 EIz
pL/2
F
6L
0
1y
{ }

2 /12
2
z2
EIz
pl
M
2
L
0
1z
+
=
F2y
6 L z3
pL
L3 0
pL/2
F3y
6 L 6 L
z2 =
do`
u:
1
8
11
13
pL , M1z =
pL2 , F2y = pL , F3y =
pL
28
14
7
28
(
)
13 8 11
+ +
2 =0
F1y + F2y + F3y 2 pL = pL
28 7 28
(
)
1
8 22
2
2
+0+0+ +
2 =0
M1z + M2z + M3z + LF2y + 2 LF3y 2 pL = pL
14
7 28
F1y =
Calcul des eorts

el
ementaires et repr
element 1-2 :
6
13 pL
Ty1
6L
= F1y
EI

2 =M
2
1
pL
L
2
pL
Mfz1
z 2L
1z
{ }
=
= 3
6 L z2
6
15 pL
Ty2
L
12
28
L
3 pL2
Mfz2
4 L2
u Mfz (xm ) = 57 pL2 /1568 0.036 pL2 .
73
element 1-2 :
Ty2
6L
6L { }
6
17 pL
EI
2
2
pL
1
2 L2
Mfz2
L
3
pL
z 4L
z2
= 3
=
Ty3
6 L 6 L z3
6
11 pL
= F3y
L
12
28
Mfz3
2 L2 4 L2
L
0
u Mfz (xm ) = 121 pL2 /1568 0.077 pL2 .
Dimensionnement
max(|Mfz (x)|) = 3 pL2 /28
Wel.z >
3 15 6002
3 pL2
=
= 1929 mm3 = 1.93 cm3
28 E
28 300
MEF FLE 9 : appui

elastique
La poutre droite representee sur la gure est encastree en 1 et repose sur un appui elastique de
raideur k en 2. Soient E le module de Young du materiau et Iz le moment quadratique de la section
droite.
La poutre porte sur toute sa longueur une force repartie dintensite lineique p avec p > 0.
1. Calculer les deplacements inconnus et les actions de liaison.
74

2. La poutre est un carre plein de cote c. On donne :
L = 800 mm , c = 50 mm , E = 210000 MPa , p = 10 N/mm , k = 20000 N/mm
Solution
12
EIz
6L
L3 12
6L

6 L 12 6 L
v1 = 0
F1y =?
6

pL
2
2
4L
6 L 2 L z1 = 0
L
M1z =?
=
+
6 L 12 6 L
6
v =?
F2y = kv2
12
2 L2 6 L 4 L2
L
z2 =?
M2z = 0
Les deplacements inconnus sont les solutions de lequation :

(
[
] [
]) { } {
}
EIz 12 6 L
k 0
v2
pL/2
=
+
z2
pL2 /12
0 0
L3 6 L 4 L2
do`
u:
v2 =
1
pL4
8 EIz 1 + C/3
z2 =
pL3 1 C/24
6 EIz 1 + C/3
C=
kL3
EIz
Remarque : la rotation de la section 2 est nulle pour C = 24 soit :

k = 24
EIz
L3
Les actions de liaisons sont :

F1y = pL
1 + 5 C/24
1 + C/3
M1z =
pL2 1 + C/12
2 1 + C/3
F2y = kv2 =
pL
C
8 1 + C/3
La gure suivante represente la variation de v2 , z2 et F2y en fonction de k.
Application numerique :
v2 = 1.453 101 mm ,
z2 = 7.027 104 rad
MEF FLE 10
La poutre droite representee sur la gure est encastree en 1 et 3. Elle est en acier de module de
Young E et de limite elastique E . Soit Iz le moment quadratique de la section droite.
75
La poutre porte entre les nuds 2 et 3 une force uniformement repartie. Soit (0, p, 0) avec p > 0 la
force par unite de longueur.
E = 300 MPa , L = 30 cm , p = 2500 daN/m
La poutre est un rond plein de diamètre D. Calculer D.
Solution
Assemblage
12
6L
6 L 4 L2
EIz
12 6 L
2
L3
6L 2L
0
0
0
0
0
F1y =?
v1 = 0
12 6 L
0
0
0
M
=?
=
0
6 L 2 L2
0
0
1z
z1
pL/2
F2y = 0
24
0
12 6 L v2 =?
+
=
pL2 /12
M2z = 0
z2 =?
0
8 L2 6 L 2 L2
pL/2
F3y =?
v3 = 0
12 6 L 12 6 L
pL2 /12
M3z =?
z3 = 0
6 L 2 L2 6 L 4 L2
Calcul des d
[
]{ } {
}
EIz 24
0
v2
pL/2
=
pL2 /12
L3 0 8 L2 z2
do`
u:
1 pL4
1 pL3
, z2 =
48 EIz
96 EIz
0
F1y
12 6 L { }
2
EIz
0
6 L 2 L v2 = M1z +
pL/2
F
L3 12 6 L z2

2
3y
pL /12
M3z
6 L 2 L2
v2 =
do`
u:
3
5
13
11
pL , M1z =
pL2 , F3y =
pL , M3z = pL2
16
48
16
48
(
)
3
13
+0+
1 =0
F1y + F2y + F3y pL = pL
16
16
(
)
3 2
5
11
26 3
+0
+0+
=0
M1z + M2z + M3z + LF2y + 2 LF3y pL = pL
2
48
48
16 2
F1y =
76
Calcul des eorts

el
ementaires et repr
element 1-2 :
12 6 L { }
3 pL/16
Ty1
= F1y
EI
2
2 /48 = M
v
5
pL
Mfz1
z 6 L 2 L
2
1z
=
= 3
12 6 L z2
3 pL/16
Ty2
L
Mfz2
6 L 4 L2
pL2 /12
element 2-3 :
Ty2
12
6L { }
6
3 pL/16
EI

2 /12
pL
4 L2
v
Mfz2
L
pL
z 6L
2
= 3
=
12 6 L z2
Ty3
6
= F3y
13 pL/16
L
12
2
2
Mfz3
6L 2L
L
11 pL /48
= M3z
u Mfz (xm ) = 155 pL2 /1536 0.101 pL2 .
Dimensionnement
max(|Mfz (x)|) = 11 pL2 /48
Wel.z =
do`
u:
D3 >
D3
11 pL2
>
32
48 E
22 pL2
22 25 3002
=
3 E
3 300
D > 25.97 mm
MEF FLE 11
La poutre droite representee sur la gure est encastree en 1 et 3. Elle est en acier de module de
Young E et de limite elastique E . Soit Iz le moment quadratique de la section droite.
77
La poutre porte entre les nuds 1 et 2 une force uniformement repartie dintensite (0, 3 p, 0) et entre
les nuds 2 et 3 une force uniformement repartie dintensite (0, p, 0) avec p > 0.
4. Application numerique : E = 230 MPa , L = 1.2 m , p = 3000 daN/m. Calculer le module de
resistance à la exion Wel.z .
Solution
Assemblage
3 pL/2
F1y =?
v1 = 0
12
6 L 12 6 L
0
0
z1 = 0
6 L 4 L2 6 L 2 L2
pL2 /4
M1z =?
0
0
EIz
2
pL
F
=
0
v
=?
12
6
L
24
0
12
6
L
2y
2
+
=
2
pL2 /6
M2z = 0
z2 =?
0
8 L2 6 L 2 L2
L3
6L 2L
pL/2
F3y =?
v3 = 0
0
12 6 L 12 6 L
2
2
2
pL /12
M3z =?
z3 = 0
0
0
6L 2L
6 L 4 L
Calcul des d
[
]{ } {
}
EIz 24
0
v2
2 pL
=
pL2 /6
L3 0 8 L2 z2
do`
u:
v2 =
1 pL4
12 EIz
z2 =
1 pL3
48 EIz
3 pL/2
F1y
12 6 L { }
2
2
EIz
6 L 2 L v2 = M1z + pL /4
L3 12 6 L z2

pL/2
F3y
pL2 /12
M3z
6 L 2 L2
do`
u:
F1y =
21
pL ,
8
M1z =
19 2
pL
24
F3y =
11
pL ,
8
M3z =
13 2
pL
24
78
(
)
21
11
F
+
F
+
F
4
pL
=
pL
+
0
+
4
=0
1y
2y
3y
8
8
(
)
3 2 3 2
19
13
22
+0
+0+
3 =0
M1z + M2z + M3z + LF2y + 2 LF3y pL pL = pL
2
2
24
24
8
Calcul des eorts

el
ementaires et repr
element 1-2 :
12 6 L { }
Ty1
6
21 pL/8
= F1y
EI

2
2 /24 = M
pL
Mfz1
6
L
2
L
v
L
19
pL
z
2
1z
= 3
12 6 L z2 4 6 = 3 pL/8
T
L
y2
Mfz2
6 L 4 L2
L
pL2 /3
Leort tranchant sannule pour
element 2-3 :
Ty2
12
EI
Mfz2
6L
z
= 3
12
T
L
y3
Mfz3
6L
xm = 7 L/8 do`
u Mfz (xm ) = 137 pL2 /384 0.357 pL2 .
3 pL/8
6
6L { }

2 /3
pL
pL
L
4 L2
v
2
= F3y
6
6 L z2
11 pL/8
12

= M3z
13 pL2 /24
L
2 L2
Dimensionnement
max(|Mfz (x)|) = 19 pL2 /24
Wel.z >
19 pL2
19 30 12002
=
24 E
24 230
Wel.z > 148696 mm3 = 148.7 cm3
MEF FLE 12
Soit Iz le moment quadratique de la section droite. La poutre est encastree en 1 et repose sur un
appui simple en 2 et 3.
79
La poutre porte sur le troncon 2 3 une force uniformement repartie. Soit (0, p, 0) avec p > 0 la
force par unite de longueur.
E = 300 MPa , L = 100 cm , p = 1500 daN/m
La poutre a une section carree de cote c. Calculer c.
Solution
Assemblage
0
F1y =?
v1 = 0
12
6 L 12 6 L
0
0
6 L 4 L2 6 L 2 L2
0
M
=?
=
0
0
0
1z
z1
EIz 12 6 L 24
pL/2
F2y =?
0
12 6 L v2 = 0
=
+
2
2
2
2
pL /12
M2z = 0
z2 =?
0
8L
6 L 2 L
L3
6L 2L
pL/2
F
=?
0
12 6 L 12 6 L
v
=
0
3y
3
2
2
2
pL /12
M3z = 0
z3 =?
0
0
6L 2L
6 L 4 L
Calcul des d
[
]{ }
{ }
EIz 8 L2 2 L2 z2
pL L
=
L
L3 2 L2 4 L2 z3
12
do`
u:
z2 =
1 pL3
56 EIz
z3 =
5 pL3
168 EIz
0
F
6L
0
1y
{ }

2
z2
EIz
0
M
2
L
0
1z
+
=
pL/2
F
6 L z3
L3 0

2y
pL/2
F3y
6 L 6 L
80
do`
u:
3
1
19
3
pL , M1z = pL2 , F2y =
pL , F3y = pL
28
28
28
7
(
)
3
19 3
+ 1 =0
F1y + F2y + F3y pL = pL +
28 28 7
(
)
3 2
1
19 6 3
+
=0
M1z + M2z + M3z + LF2y + 2 LF3y pL = pL + 0 + 0 +
2
28
28 7 2
F1y =
Calcul des eorts

el
ementaires et repr
element 1-2 :
T
3
pL/28
6
L
y1
= F1y
EI
2
2 /28 = M
Mfz1
pL
z 2L
1z
{z2 } =
= 3
T
6
L
3
pL/28
y2
Mfz2
4 L2
pL2 /14
element 2-3 :
4 pL/7
Ty2
6
6L
6L { }
EI
2
pL L pL2 /14
Mfz2
2 L2
z2
z 4L
= 3
=
= F3y
Ty3
6
6 L 6 L z3
3 pL/7
L
12

0
Mfz3
L
2 L2 4 L2
u Mfz (xm ) = 9 pL2 /98.
Dimensionnement
max(|Mfz (x)|) = 9 pL2 /98
Wel.z =
do`
u:
c3 >
9 pL2
c3
>
6
98 E
27 pL2
27 15 10002
=
49 E
49 300
c > 30.20 mm
Methodes energetiques
81
MEF FLE 13 : appui

elastique
La poutre droite representee sur la gure, de module de Young E et de moment quadratique Iz , repose
sur un appui elastique de raideur k en 1 et 2.
Elle porte sur toute sa longueur une force repartie dintensite lineique p avec p > 0.
1. Calculer les deplacements nodaux.
2. Application numerique.
L = 1 m , E = 200000 MPa , k = 2 106 N/m , p = 15000 N/m
La poutre est un carre plein de cote c = 50 mm.
Calculer v1 , v2 , 1z et 2z .
Solution

F1y = kv1
v1
12
6 L 12 6 L
pL

2
2
EIz
M
=
0
L
6
L
4
L
6
L
2
L
1z
z1
+
=
F = kv2
v
6
L3 12 6 L 12 6 L
12
2y

2

M2z = 0
6 L 2 L2 6 L 4 L2
z2
L
La symetrie du système implique :
v1 = v2
z1 = z2
do`
u les deplacements :
v1 = v2 =
z1 = z2 =
pL
15000 1
=
= 3.75 103 m
2k
2 2 106
1 pL3
15 10003
=
= 6 103 rad
24 EIz
14 200000 (504 /12)
Chapitre 3
M
ethodes
energ
etiques : poutres
3.1
3.1.1
Rappels
Expression de l
energie de d
eformation en fonction des forces appliqu
ees :
formule de Clapeyron
Soit un solide soumis au système de charges {F }. {U } est le vecteur deplacement associe à {F }. Si le

comportement du materiau est elastique et lineaire, lenergie de deformation du solide est :
1
{F }T {U }
2
1
1
= {F }T [C] {F } = {U }T [K] {U }
2
2
Edef =
(3.1.1)
o`
u:
[C] est la matrice de souplesse ou matrice de exibilite.
[K] est la matrice de rigidite.
et :
{F } = [K] {U } ,
{U } = [C] {F }
[K] = [C]1
(3.1.2)
Remarque : lenergie de deformation sexprime en joule (1 J = 1 N.m = 1 kg.m2 .s2 ).
3.1.2
Th
eor`
eme de r
eciprocit
e de Maxwell-Betti
Si un solide est soumis à deux systèmes de forces {F } et {F }, le travail de {F } quand on applique {F }

est egal au travail de {F } quand on applique {F }.
On en deduit :
[C] = [C]T
[K] = [K]T
(3.1.3)
La matrice de exibilite et la matrice de rigidite sont symetriques.
3.1.3
Th
eor`
eme de Castigliano
Lenergie de deformation Edef etant exprimee en fonction des forces appliquees {F }, le deplacement Ui
de la composante Fi est egal à la derivee de Edef par rapport à Fi :
Ui =
Edef
Fi
(3.1.4)
84
3.1.4
Th
eor`
eme de M
enabr
ea
Soit un syst`
eme hyperstatique. Si r est le nombre de reactions dappui et p le nombre dequations
fournies par la statique, le degre dhyperstaticite est egal au nombre n = r p de reactions surabondantes (inconnues hyperstatiques).
Les valeurs que prennent les inconnues hyperstatiques Ri=1,...,n rendent stationnaire lenergie de deformation :
Edef
= 0 i = 1, . . . , n
(3.1.5)
Ri
3.1.5
Energie
de d
eformation dune poutre
Lenergie de deformation emmagasinee dans un troncon de poutre de longueur dx est egale à :

dEdef =
N2
dx
(eort normal)
2 EA
Ty2
Mfz2
dx +
dx
(exion dans le plan {x, y})
+
2 EIz
2 GAky
Mfy2
Tz2
+
dx +
dx
(exion dans le plan {x, z})
2 EIy
2 GAkz
Mt2
+
dx
(torsion)
2 GJ
avec :
G=
E
2 (1 + )
(materiau isotrope)
(3.1.6)
(3.1.7)
o`
u:
E et sont respectivement le module de Young et le coecient de Poisson du materiau.
A est laire de la section.
J est la constante de torsion de Saint-Venant.
Iy et Iz sont les moments quadratiques centraux principaux de la section.
ky et kz sont les coecients de cisaillement.
N est leort normal.
Mfy et Mfz sont les moments echissants.
Ty et Tz sont les eorts tranchants.
M t est le moment de torsion.
3.1.6
Formules math
ematiques utiles
Identite remarquable :
(a + b)2 = a2 + b2 + 2 ab
(3.1.8)
(un ) = n un1 u
(3.1.9)
Si u est une fonction de x :

sin (u) = cos(u) u
On en deduit :
cos (u) = sin(u) u
(3.1.10)
]L
[
Ln+1
(L x)n+1
=
(L x)n dx =
n+1
n+1
0
(3.1.11)
85
Relations trigonometriques :
cos2 =
1 + cos 2
2
sin2 =
1 cos 2
2
cos sin =
sin 2
2
Linverse dune matrice carree [M ] de dimension 2 :

[
]
M11 M12
[M ] =
M21 M22
est :
1
[M ]
3.2
]
[
1
M22 M12
=
det [M ] M21 M11
det [M ] = M11 M22 M12 M21
(3.1.12)
(3.1.13a)
(3.1.13b)
Exercices
ENER 1
La poutre representee sur la gure est encastree en 1. Soit E le module de Young du materiau. Laire
de la section droite est 3 A entre les sections 1 et 2 et A entre les sections 2 et 3.
Elle porte en 2 une force de composantes (F2 , 0, 0) et en 3 une force de composantes (F3 , 0, 0).
1. Expression de leort normal N (x).
2. Calculer lenergie de deformation.
3. Calculer les deplacements u2 et u3 des sections 2 et 3.
4. En deduire les matrices de exibilite et de rigidite de la poutre.
Solution
Leort normal est egal à :
N12 = F2 + F3
N23 = F3
do`
u lenergie de deformation :
L
2L
)
1
1
L ( 2
2
Edef =
(F2 + F3 ) dx +
F32 dx =
F2 + 4 F32 + 2 F2 F3
6 EA 0
2 EA L
6 EA
et le deplacement horizontal des sections 2 et 3 :
u2 =
Edef
L
=
(F2 + F3 ) ,
F2
3 EA
u3 =
Edef
L
=
(F2 + 4 F3 )
F3
3 EA
On en deduit la matrice de exibilite :

[
]{ }
{ }
L
1 1 F2
u2
=
u3
3 EA 1 4 F3
|
{z
}
[C]
86
et la matrice de rigidite :
{ }
[
]{ }
EA 4 1 u2
F2
=
F3
u3
L 1 1
|
{z
}
[K]=[C]1
ENER 2
La poutre representee sur la gure est encastree en 1. Soit E le module de Young du materiau. Laire
de la section droite est 2 A entre les sections 1 et 2 et A entre les sections 2 et 3.
Elle porte en 3 une force de composantes (F, 0, 0).

3. Calculer le deplacement u3 de la section 3.
Solution
D
eplacement de la section 3
N (x) = F
do`
Edef
1
=
4 EA
1
F dx +
2 EA
F 2 dx =
0
3 F 2L
4 EA
et le deplacement horizontal de la section 3 :

u3 =
Edef
3FL
=
F
2 EA
D
Introduisons une force auxiliaire Q en 2. Leort normal est alors egal à :
N12 = F + Q ,
N23 = F
87
do`
u:
N12
=1 ,
Q
N23
=0
Q
Lenergie de deformation est egale à :

Edef
1
=
4 EA
L
2
N12
1
dx +
2 EA
L
2
N23
dx
do`
u le deplacement horizontal de la section 2 :

L
L

Edef
1
N12
FL
1

u2 =
=
N12
F dx =
dx
=

Q Q=0
2 EA 0
Q
2 EA 0
2 EA
Q=0
ENER 3
La poutre representee sur la gure est encastree en 1. Soient A laire de la section droite et E le
module de Young du materiau.
La poutre porte en 2 une force de composantes (F, 0, 0) et sur toute sa longueur, une force uniformement repartie dintensite lineique p.
4. Calculer le deplacement u(x) de la section dabscisse x.
Solution
D
N (x) = F + p (L x)
do`
u:
N
=1
F
Lenergie de deformation de la poutre est :

Edef =
1
2 EA
N 2 dx
0
do`
L
L
L
1
N
F
p
FL
pL2
Edef
=
N
dx =
dx +
(L x) dx =
+
u2 =
F
EA 0
F
EA 0
EA 0
EA 2 EA
88
D
eplacement de la section dabscisse x
Introduisons une force auxiliaire Q dans la section dabscisse x.
Leort normal est alors egal à :

{
F + Q + p (L s)
N=
F + p (L s)
do`
u:
N
=
Q
1
0
0<s<x
x<s<L
0<s<x
x<s<L

Edef
1
=
2 EA
1
N ds +
2 EA
N 2 ds
x
do`
u le deplacement horizontal de la section dabscisse x :

x
x
Edef
1
N
1
u(x) =
ds
=
N
=
(F + p (L s)) ds
Q Q=0
EA 0
Q Q=0 EA 0
soit :
u(x) =
Fx
px
+
(2 L x)
EA 2 EA
ENER 4
La representee sur la gure est encastree en 1 et 2. Soient A laire de la section droite et E le module
de Young du materiau.
Elle porte en 2 une force de composantes (F, 0, 0).
1. Ecrire
lequation dequilibre.
2. Calculer les actions de liaisons.
Solution
Lequilibre de la poutre secrit :
F1x + F3x + F = 0
89
o`
u F1x et F3x sont les reactions dappui. Le système est donc hyperstatique dordre 1. Le système
isostatique associe est represente sur la gure ci-dessous.
F3x est linconnue hyperstatique.

N12 = F + F3x
do`
u:
N12
=1 ,
F3x
N23 = F3x
N23
=1
F3x

Edef
1
=
2 EA
L
2
N12
1
dx +
2 EA
2L
2
N23
dx
Linconnue hyperstatique est solution de lequation (theorème de Menabrea) :

Edef
=0
F3x
do`
u:
1
EA
N12
1
dx +
F3x
EA
N12
0
2L
N23
0
N23
dx = 0
F3x
soit :
(F + F3x ) L + F3x 2 L = 0
On en deduit linconnue hyperstatique F3x :
1
F3x = F
3
et laction de liaison F1x :
2
F1x = F F3x = F
3
do`
u:
N12 =
2
F
3
1
N23 = F
3

Edef
1
=
2 EA
4 2
1
F dx +
9
2 EA
2L
do`
u2 =
2 FL
Edef
=
F
3 EA
1 2
F 2L
F dx =
9
3 EA
90
ENER 5 :
el
ement de poutre de section variable
La poutre representee sur la gure est un rond plein dont le diamètre varie lineairement entre D1
et D2 . La section 1 est encastree. Soit E le module de Young du materiau.
La section 2 porte une force de composantes (F, 0, 0).

3. En deduire la matrice de rigidite de lelement de poutre 1 2.
Solution
D
N (x) = F
Lenergie de deformation est :
Edef =
0
F2
dx
2 EA
A=
)2
((
x)
x
1
D1 + D2
4
L
L
do`
u le deplacement de la section 2 :
Edef
F
u2 =
=
F
E
L
4 F L
1
1
((
)
dx =
x)
x
A
E D2 D1
1
D1 + D2
L
L
0
soit :
u2 =
4FL
ED1 D2
On en deduit :
F = k u2
ED1 D2
k =
4L
Matrice de rigidit
e de l
el
ement 1 2
Considerons lelement 1 2 (gure ci-dessous) dont les deplacements nodaux sont u1 et u2 .
91
Decomposons le mouvement de cet element en un mouvement de corps rigide (R) et un mouvement

avec deformation (D) :
Lenergie de deformation de lelement est :

Edef = Edef (R) +Edef (D)
| {z }
=0
do`
u:
Edef =
La quantite u2 u1 secrit :
1
1
{u}T [ k ] {u} = k (u2 u1 )2
2
2
{ }
[
] u1
[
]
u2 u1 = 1 1
= 1 1 {u}
u2
do`
Edef
[ ]
[
]
1
T 1
1 1 {u}
= {u} k
1
2
|
{z
}
[k]
ED1 D2
[k] =
4L
]
1 1
1 1
ENER 6
On considères la poutre daxe x representee sur la gure. Soit G le module delasticite transversal du
materiau. La constante de torsion est J entre les sections 1 et 2 et 3 J entre les sections 2 et 3.
La section 1 est encastree. Les sections 2 et 3 portent un couple de torsion dintensites respectives C2
et C3 .
92

1. Expression du moment de torsion Mt (x).
3. Calculer les rotations x2 et x3 des sections 2 et 3.
4. En deduire les matrices de exibilite et de rigidite de la poutre.
Solution
Le moment de torsion dans la poutre est :
Mt12 = C2 + C3
Mt23 = C3
do`
2L
L
1
L
1
2
2
Mt12
dx +
Mt23
dx =
(6 C22 + 7 C32 + 12 C2 C3 )
Edef =
2 GJ 0
6 GJ 0
6 GJ
et la rotation des sections 2 et 3 :
x2 =
2L
Edef
=
(C2 + C3 ) ,
C2
GJ
x3 =
Edef
L
=
(6 C2 + 7 C3 )
C3
3 GJ

[
]{ }
{ }
L
6 6 C2
x2
=
x3
3 GJ 6 7 C3
|
{z
}
[C]
C2
C3
[
]{ }
GJ 7 6 x2
=
x3
2 L 6 6
|
{z
}
[K]=[C]1
ENER 7 :
el
ement de poutre de section variable
La poutre representee sur la gure est un rond plein dont le diamètre varie lineairement entre D1
et D2 . La section 1 est encastree. Soit G le module delasticite transversal du materiau.
La section 2 porte un couple de composantes (C, 0, 0).

2. Calculer la rotation x2 de la section 2.
3. En deduire la matrice de rigidite de lelement de poutre 1 2.
93
Solution
Rotation de la section 2
Le moment de torsion dans la poutre est :
Mt (x) = C
L
C2
Edef =
dx
0 2 GJ
J=
)4
((
x)
x
1
D1 + D2
32
L
L
do`
u la rotation de la section 2 :
x2 =
Edef
C
=
C
G
1
32 C
L
dx =
((
J
G
3 (D2 D1 )
soit :
x2 =
x)
x )3
1
D1 + D2
L
L
0
32 CL D12 D1 D2 + D2
3 G
D13 D23
On en deduit :
C = k x2
3 G
D13 D23
k =
32 L D12 + D1 D2 + D22
Matrice de rigidit
e de l
el
ement 1 2
Considerons lelement 1 2 (gure ci-dessous) dont les deplacements nodaux sont x1 et x2 .
Decomposons le mouvement de cet element en un mouvement de corps rigide (R) et un mouvement

avec deformation (D) :
Lenergie de deformation de lelement est :

Edef = Edef (R) +Edef (D)
| {z }
=0
do`
u:
Edef =
1
1
{u}T [ k ] {u} = k (x2 x1 )2
2
2
94
La quantite x2 x1 secrit :
x2 x1
{ }
[
] x1
[
]
= 1 1
= 1 1 {u}
x2
do`
Edef
[ ]
[
]
1
T 1
1 1 {u}
= {u} k
1
2
|
{z
}
[k]
D13 D23
3 G
[k] =
32 L D12 + D1 D2 + D22
]
1 1
1 1
ENER 8
La poutre 1 2 de longueur L et de section constante (moment quadratique : Iz ) est encastree en 1.
Elle est en acier de module de Young E. Elle porte en 2 un couple ponctuel de composantes (0, 0, C).
On neglige linuence de leort tranchant : modèle de Bernoulli.

1. Expression du moment echissant Mfz (x).
3. Calculer la rotation z2 de la section 2.
4. Calculer le deplacement v2 de la section 2.
Solution
Le moment echissant dans la poutre est C, do`
L
1
C 2L
Edef =
C 2 dx =
2 EIz 0
2 EIz
et la pente en 2 :
z2 =
D
eplacement vertical de la section 2
Introduisons en 2 une force auxiliaire (0, F, 0).
Edef
CL
=
C
EIz
95
Le moment echissant est alors egal à :

Mfz (x) = C + F (L x)
do`
u:
Mfz
=Lx
F
La èche en 2 est egale à :

L
L

Edef
1
Mfz
C
CL2

v2 =
=
Mf
dx
=
(L
x)
dx
=
z

F F =0
EIz 0
F
EIz 0
2 EIz
F =0
ENER 9
La poutre 1 2 3 de longueur L et de section constante (moment quadratique : Iz ) repose en 1 et 3
sur un appui simple. Elle est en acier de module de Young E. Elle porte en son milieu 2 une force de
composantes (0, F, 0).

3. Calculer la èche en 2.
Solution
Le moment echissant dans la poutre est :
Mfz (1 2) =
F
x ,
2
Mfz (2 3) =
F
(L x)
2

Edef = 2 Edef (1 2) =
do`
u la èche en 2 :
v2 =
1
2 EIz
L/2
F2 2
F 2 L3
x dx =
4
96 EIz
Edef
F L3
=
F
48 EIz
96
ENER 10
La poutre 1 2 de longueur L et de section constante (moment quadratique : Iz ) repose en 1 et 2
sur un appui simple. Elle est en acier de module de Young E. Elle porte en 1 un couple de composantes (0, 0, C).

Solution
Le moment echissant est egal à :
Mfz (x) = C
do`
u:
x)
L
(
x)
Mfz
= 1
C
L

Edef =
do`
z1
Edef
1
=
=
C
EIz
1
2 EIz
Mfz2 dx
0
Mfz
C
Mfz
dx =
C
EIz
L(
Introduisons un couple auxiliaire (0, 0, M ) en 2.
Le moment echissant devient :

Mfz (x) = C
(
x
x)
+M
1
L
L
x )2
CL
dx =
L
3 EIz
97
do`
u:
Mfz
x
=
M
L
et la rotation de la section 2 :
z2

L

Edef
Mfz
1
=
Mfz
=
dx

M M =0
EIz 0
M
M =0
soit :
z2
C
=
EIz
L(
x) x
CL
dx =
L L
6 EIz
ENER 11
La poutre 1 2 de longueur L et de section constante (moment quadratique : Iz ) repose en 1 et 2
sur un appui simple. Elle est en acier de module de Young E. Elle porte en 2 un couple de composantes (0, 0, C).

Solution
x
L
Mfz (x) = C
do`
u:
x
Mfz
=
C
L

Edef
1
=
2 EIz
Mfz2 dx
0
do`
z2
1
Edef
=
=
C
EIz
Mfz
C
Mfz
dx =
C
EIz
L(
x )2
CL
dx =
L
3 EIz
98
Introduisons un couple auxiliaire (0, 0, M ) en 1.
Le moment echissant devient :

Mfz (x) = M
do`
u:
x)
x
+C
L
L
(
Mfz
x)
= 1
M
L
z1

L

Mfz
Edef
1
Mfz
dx
=
=

M M =0
EIz 0
M
M =0
soit :
z1
C
=
EIz
x(
x)
CL
1
dx =
L
L
6 EIz
ENER 12
Soit la poutre (12) de longueur L et de section droite constante. Soient A, ky et Iz les caracteristiques
de la section droite. E et sont respectivement le module de Young et le coecient de Poisson du
materiau.
La poutre est encastree en 1. Elle porte en 2 une force de composantes (0, F, 0) et un couple de composantes (0, 0, C).
1. Donner lexpression de leort tranchant Ty (x) et du moment echissant Mfz (x).
2. Calculer les deplacements v2 et z2 de la section 2.
3. En deduire lexpression des matrices de souplesse et de rigidite de la structure.
Solution
Leort tranchant et le moment echissant dans la poutre sont :
Ty (x) = F
Mfz (x) = C + F (L x)
do`
u:
99
Ty
=1
F
Ty
=0
C
Mfz
=Lx ,
F
Mfz
=1
C
Mod`
ele de Bernoulli
1
=
2 EIz
Edef
Mfz2 dx
0
do`
u les deplacements de la section 2 :
1
Edef
=
F
EIz
v2 =
z2 =
Edef
1
=
C
EIz
Mfz
Mfz
dx
F
Mfz
Mfz
dx
C
Ces equations secrivent :

L
C
(L x) dx +
(L x) dx
EIz 0
0
L
L
F
C
=
(L x) dx +
dx
EIz 0
EIz 0
F
v2 =
EIz
z2
do`
u:
v2 =
F L3
CL2
+
3 EIz
2 EIz
z2 =
CL
F L2
+
2 EIz
EIz

{ }
[ 2
]{ }
L
v2
2L 3L F
=
z2
6
C
6 EIz 3 L
|
{z
}
[C]
]{ }
{ }
[
EIz 12 6 L
v2
F
= 3
2
z2
C
6 L 4 L
L
|
{z
}
[K]=[C]1
Mod`
ele de Timoshenko
Edef
1
=
2 EIz
Mfz2
0
1
dx +
2 GAky
Ty2 dx
0
do`
u les deplacements de la section 2 :
Edef
1
v2 =
=
F
EIz
Mfz
1
Mfz
dx +
F
GAky
Ty
0
Ty
dx
F
100
z2
Edef
1
=
=
C
EIz
L
0
Mfz
1
dx +
Mfz
C
GAky
Ty
0
Ty
dx
C

F
v2 =
EIz
C
(L x) dx +
EIz
z2 =
F
EIz
F
(L x) dx +
GAky
(L x) dx +
C
EIz
dx
0
dx
0
do`
u:
v2 =
F L3
CL2
FL
+
+
3 EIz
2 EIz
GAky
z2 =
F L2
CL
+
2 EIz
EIz

(
)
{ }
{ }
L 2 L2 1 + y
v2
3 L F
=
4
C
z2
6 EIz
3L
6
{z
}
|
y =
12 EIz
L2 GAky
[C]
{ }
[
]{ }
EIz
F
12
6 L
v2
= 3
C
L (1 + y ) 6 L (4 + y ) L2 z2
|
{z
}
[K]=[C]1
Remarque : de lexpression de [K], on deduit la matrice de rigidite dun element de poutre pour la
exion dans le plan {x, y} (http://iut.univ-lemans.fr/ydlogi/cours/exion xy.pdf).
ENER 13
La poutre 1 2 de longueur L et de section constante repose en 1 et 2 sur un appui simple. Soient A,
ky et Iz les caracteristiques de la section droite. E et sont respectivement le module de Young et
le coecient de Poisson du materiau. Elle porte en 1 un couple de composantes (0, 0, C1 ) et en 2 un
couple de composantes (0, 0, C2 ).
1. Expression de leort tranchant Ty (x) et du moment echissant Mfz (x).

2. Calculer les rotations z1 et z2 des sections 1 et 2.
3. En deduire lexpression de la matrice de souplesse et de la matrice de rigidite de la structure.
101
Solution
Leort tranchant et le moment echissant sont :
Ty (x) =
do`
u:
C1 + C2
L
Ty
Ty
1
=
=
C1
C2
L
(
x)
x
Mfz (x) = C1 1
+ C2
L
L
(
x)
Mfz
= 1
C1
L
x
Mfz
=
C2
L
Mod`
ele de Bernoulli
1
=
2 EIz
Edef
Mfz2 dx
0
do`
u la rotation des sections 1 et 2 :
z1
1
Edef
=
=
C1
EIz
z2 =
Edef
1
=
C2
EIz
Mfz
Mfz
dx
C1
Mfz
Mfz
dx
C2

z1
z2
L(
x )2
C2
1
dx
L
EIz
x(
x)
1
dx
L
0
0 L
L(
L( )
C2
C1
x) x
x 2
dx +
dx
=
1
EIz 0
L L
EIz 0
L
C1
=
EIz
do`
u:
C1 L
C2 L
3 EIz
6 EIz
C1 L
C2 L
=
+
6 EIz
3 EIz
z1 =
z2
{ }
[
]{ }
L
z1
2 1 C1
=
C2
z2
6 EIz 1 2
|
{z
}
On en deduit :
[C]
et :
C1
C2
[
]{ }
2 EIz 2 1 z1
=
1 2 z2
L
|
{z
}
[K]=[C]1
Mod`
ele de Timoshenko
Edef
1
=
2 EIz
Mfz2
0
1
dx +
2 GAky
Ty2 dx
0
102
do`
z1
Edef
1
=
=
C1
EIz
1
Edef
=
C2
EIz
z2 =
Mfz
1
Mfz
dx +
C1
GAky
Mfz
0
1
Mfz
dx +
C2
GAky
Ty
Ty
dx
C1
Ty
Ty
dx
C2

L(
L (
L
x
C1 + C2
C1
x )2
C2
x)
1
dx
z1 =
1
dx
1
dx +
EIz 0
L
EIz 0 L
L
GAky 0
L2
z2
C1
=
EIz
L(
x) x
C2
1
dx +
L L
EIz
do`
u:
z1 =
L(
x )2
1
dx +
L
GAky
C1 + C2
dx
L2
C2 L
C1 + C2
C1 L
+
3 EIz
6 EIz
LGAky
z2 =
C2 L
C1 + C2
C1 L
+
+
6 EIz
3 EIz
LGAky
y
{ }
L 2+ 2
z1
=
y
z2
6 EIz
1 +
{z2
|
y { }
1 +
2 C1
y C 2
2+
2 }
y =
12 EIz
L2 GAky
[C]
[
]{ }
{ }
EIz
4 + y 2 y
C1
z1
=
C2
z2
L (1 + y ) 2 y 4 + y
|
{z
}
[K]=[C]1
Remarque : de lexpression de [K], on deduit la matrice de rigidite dun element de poutre pour la
exion dans le plan {x, y} (http://iut.univ-lemans.fr/ydlogi/cours/exion xy.pdf).
ENER 14
La poutre 1 2 de longueur L et de section constante (moment quadratique : Iz ) repose en 1 et 2 sur
un appui simple. Elle est en acier de module de Young E. Elle porte sur toute sa longueur une force
uniformement repartie dintensite lineique (0, p, 0) avec p > 0.
103
Calculer les rotations z1 et z2 des sections 1 et 2.

Solution
Introduisons un couple auxiliaire (0, 0, C1 ) en 1 et un couple auxiliaire (0, 0, C2 ) en 2.

Mfz (x) =
(
pL (
x)
x
x)
x 1
C1 1
+ C2
2
L
L
L
do`
u:
(
Mfz
x)
x
Mfz
= 1
=
,
C1
L
C2
L
L
1
Mfz2 dx
Edef =
2 EIz 0
do`

L
Edef
1
Mfz
z1 =
=
C1 C1 =C2 =0
EIz 0

L
Edef
1
Mfz
z2 =
=
C2 C1 =C2 =0
EIz 0

Mfz
dx
C1
C1 =C2 =0

Mfz
dx
C2
C1 =C2 =0

z1 =
z2
(
x)(
x)
pL3
x 1
1
dx =
L
L
24 EIz
0
L (
pL
x) x
pL3
=
x 1
dx =
2 EIz 0
L L
24 EIz
pL
2 EIz
ENER 15
La poutre (1 2) de longueur L et de section droite constante est encastree en 1.
Soient A , Iz et ky les caracteristiques de la section droite. E et sont respectivement le module de

Young et le coecient de Poisson du materiau.
104
Elle porte en 2 une force de composantes (0, F, 0).

On prend en compte linuence de leort tranchant : modèle de Timoshenko.
5. La section droite est rectangulaire : A = b h
Iz =
b h3
12
ky =
5
6
= 0.3
Etudier
linuence de leort tranchant en fonction de L/h = 20, 10, 5, 4, 3.
Solution
D
eplacement vertical de la section 2
Ty (x) = F
Mfz (x) = F (L x)
On en deduit :
Ty
Mfz
=1 ,
=Lx
F
F
L
L
1
1
2
Edef =
Mfz dx +
T 2 dx
2 EIz 0
2 GAky 0 y
On en deduit le deplacement vertical de la section 2 :
L
L
Ty
1
Mfz
1
Edef
=
Mfz
dx +
Ty
dx
v2 =
F
EIz 0
F
GAky 0
F
soit :
v2 =
v2 secrit :
v2 =
F
EIz
F L3
3 EIz
(L x)2 dx +
1
GAky
F dx =
0
(
)
(
)
y
3 EIz
F L3
1+ 2
=
1+
L GAky
3 EI
4
| {z z}
Bernoulli
Introduisons un couple auxiliaire C en 2.
F L3
FL
+
3 EIz
GAky
,
y =
12 EIz
L2 GAky
105

Ty (x) = F
On en deduit :
Mfz (x) = F (L x) + C
Ty
=0 ,
C
Mfz
=1
C
La rotation de la section 2 est :

z2

L

Edef
1
Mfz
F L2

=
=
Mf
dx
=
z

C C=0
EIz 0
C
2 EIz
C=0
Application num
erique
Pour une section droite rectangulaire, il vient :
12
y =
(1 + )
5
( )2
h
L
do`
u la variation de y en fonction de L/h pour = 0.3 :
L/h
y
20
0.0078
10
0.0312
5
0.125
4
0.195
3
0.347
ENER 16
materiau.
La poutre est encastree en 1. Elle porte sur toute sa longueur une force uniformement repartie dont
lintensite lineique est (0, p, 0).
Calculer les deplacements v2 et z2 de la section 2.
Solution
Introduisons dans la section 2 une force auxilaire (0, F, 0) et un couple auxiliaire (0, 0, C).
106

Ty (x) = F + p (L x) ,
do`
u:
Ty
=1
F
Mfz (x) = C + F (L x) +
Ty
=0
C
Mfz
=Lx ,
F
p (L x)2
2
Mfz
=1
C

Edef
1
=
2 EIz
Mfz2
0
1
dx +
2 GAky
Ty2 dx
0
On en deduit :
la èche en 2 :

L
L

Ty
Edef
1
Mfz
1

v2 =
=
Mfz
dx
+
Ty
dx

F F =C=0
EIz 0
F
GAky 0
F
F =C=0
F =C=0
L
L
p
p
(L x)3 dx +
(L x) dx
v2 =
2 EIz 0
GAky 0
v2 =
pL4
pL2
+
8 EIz
2 GAky
Remarque : la èche en 2 secrit :

(
)
(
)
y
4 EIz
pL4
pL4
1+ 2
=
1+
v2 =
8 EIz
L GAky
8 EIz
3
y =
12 EIz
L2 GAky
la rotation de la section 2 :

L
L

Ty
Edef
Mfz
1
1

z2 =
Mf
dx
T
dx
=
+
z
y

C F =C=0
EIz 0
C
GAky 0
C
F =C=0
F =C=0
L
p
(L x)2 dx
z2 =
2 EIz 0
z2 =
pL3
6 EIz
ENER 17
materiau.
La poutre est encastree en 1. Elle porte :
107
en 2 une force de composantes (0, F, 0) et un couple de composantes (0, 0, C).

sur toute sa longueur une force uniformement repartie dintensite lineique (0, p, 0).
2. Calculer les deplacements v2 et z2 de la section 2.
Solution
Ty (x) = F + p (L x) ,
do`
u:
Ty
=1
F
Mfz (x) = C + F (L x) +
Ty
=0
C
Mfz
=Lx ,
F
p (L x)2
2
Mfz
=1
C
La èche en 2 est egale à :

L
L
Ty
Edef
1
Mfz
1
v2 =
=
Mfz
dx +
Ty
dx
F
EIz 0
F
GAky 0
F
L
L
L
C
F
p
2
=
(L x) dx +
(L x) dx +
(L x)3 dx
EIz 0
EIz 0
2 EIz 0
L
L
F
p
+
dx +
(L x) dx
GAky 0
GAky 0
soit :
v2 =
CL2
F L3
pL4
FL
pL2
+
+
+
+
2 EIz
3 EIz
8 EIz
GAky
2 GAky
Remarque : la èche en 2 secrit :

F L3
CL2
v2 =
+
2 EIz
3 EIz
(
)
(
)
y
y
pL4
1+
+
1+
4
8 EIz
3
y =
12 EIz
L2 GAky
La rotation de la section 2 est :

z2
L
Edef
1
Mfz
=
=
Mfz
dx
C
EIz 0
C
L
L
L
C
F
p
=
dx +
(L x) dx +
(L x)2 dx
EIz 0
EIz 0
2 EIz 0
soit :
z2 =
F L2
CL
pL3
+
+
2 EIz
EIz
6 EIz
108
ENER 18
Soit la poutre (1 2) de longueur L et de section droite constante. Iz est le moment quadratique de
la section droite. E est le module de Young du materiau.
La poutre est encastree en 1 et repose sur un appui simple en 2. La section 2 porte un couple de
composantes (0, 0, C).
1. Ecrire
les equations dequilibre.
Solution
Les actions de liaison sont : F1y , F2y et M1z do`
u les equations dequilibre de la poutre :
F1y + F2y = 0
M1z + LF2y + C = 0
Le système est hyperstatique dordre 1. Le système isostatique associe est represente sur la gure
ci-dessous.
F2y est linconnue hyperstatique.

Mfz (x) = C + F2y (L x)
do`
u:
Mfz
=Lx
F2y
Edef
1
=
F2y
EIz
Mfz
0
Mfz
dx = 0
F2y
soit :
109
(L x) dx + F2y
(L x)2 dx = C
L2
L3
+ F2y
=0
2
3
On en deduit F2y :
F2y =
3C
2 L
puis F1y et M1z à laide des equations dequilibre :

F1y = F2y
1
M1z = C
2

Mfz (x) = C + F2y (L x) =
C
(3 x L)
2L
do`
Edef
1
=
2 EIz
2z =
C2
C 2L
2
(3
x
L)
dx
=
4 L2
8 EIz
Edef
CL
=
C
4 EIz
ENER 19
La poutre est encastree en 1 et repose sur un appui simple en 2. Elle porte sur toute sa longueur une
force uniformement repartie dont lintensite lineique est (0, p, 0).
1. Ecrire
Solution
Les actions de liaison sont : F1y , F2y et M1z do`
u les equations dequilibre de la poutre :
F1y + F2y + pL = 0
M1z + LF2y +
1 2
pL = 0
2
Le système isostatique associe est represente sur la gure ci-dessous.
110
F2y est linconnue hyperstatique.

Mfz (x) = F2y (L x) + p
do`
u:
(L x)2
2
Mfz
=Lx
F2y

1
Edef
=
F2y
EIz
soit :
p
(L x) dx +
2
F2y
0
Mfz
0
L
Mfz
dx = 0
F2y
(L x)3 dx = F2y
L3 pL4
+
=0
3
8
On en deduit F2y :
puis F1y et M1z
3
F2y = pL
8
à laide des equations dequilibre :
5
F1y = pL ,
8
1
M1z = pL2
8
ENER 20
Les sections 1 et 2 sont encastrees. Elle porte sur toute sa longueur une force repartie dont lintensite
lineique suivant y varie lineairement entre p et 0.
1. Ecrire
111

Solution
Les actions de liaison sont : F1y , F2y , M1z et M2z do`
u lequilibre de la poutre :
1
1
pL = 0 , M1z + M2z + LF2y + pL2 = 0
2
6
Le système isostatique associe est represente sur la gure ci-dessous.
F1y + F2y +
F2y et M2z sont les inconnues hyperstatiques.

Mfz (x) = M2z + F2y (L x) +
do`
u:
Mfz
=Lx ,
F2y
p
(L x)3
6L
Mfz
=1
M2z
Les inconnues hyperstatiques sont les solutions des equations (theorème de Menabrea) :
L
L
Edef
1
Edef
1
Mfz
Mfz
=
dx = 0 ,
=
dx = 0
Mfz
Mfz
F2y
EIz 0
F2y
M2z
EIz 0
M2z
qui secrivent :
L
M2z
dx + F2y
0
et :
M2z
p
(L x) dx +
6L
(L x) dx + F2y
(L x)2 dx +
(L x)3 dx = M2z L + F2y
p
6L
(L x)4 dx = M2z
L2 pL2
+
=0
2
24
L2
L3 pL5
+ F2y
+
=0
2
3
30
On en deduit F2y et M2z :

3
pL ,
20
à laide des equations dequilibre :
F2y =
puis F1y et M1z
M2z =
1
pL2
30
7
1
pL , M1z = pL2
20
20
Application `
a la methode des elements finis : pour un element de poutre droite de section droite
constante qui porte sur toute sa longueur L une force repartie dont lintensite lineique varie lineairement entre p et 0, le vecteur force est :
F1y
21
pL
M1z
3L
=
{f } =
F2y
9
60
M2z
2 L
F1y =
112
ENER 21
Larc (1 2) de centre O et de rayon moyen R est encastre en 1. Il a une section droite constante.
Soit Iz le moment quadratique de la section droite. E est le module de Young du materiau.
Larc est soumis en 2 à un couple de composantes (0, 0, C).

On neglige linuence de leort normal et de leort tranchant.
1. Expression du moment echissant Mfz ().
4. Calculer les deplacements u2 et v2 de la section 2.
Solution
Le moment echissant dans la section droite est C, do`
Edef
1
=
2 EIz
2z =
/2
C 2R
C 2 |{z}
ds =
4 EIz
R d
Edef
CR
=
C
2 EIz
Translation de la section 2
Introduisons dans la section 2 une force auxiliaire de composantes (Fx , Fy , 0).
Le moment echissant dans la section droite est :

Mfz () = C RFx (1 cos ) + RFy sin
113
do`
u:
Mfz
= R (1 cos ) ,
Fx
Mfz
= R sin
Fy
On en deduit les deplacements de la section 2 :

/2

Edef
1
Mfz

u2 =
=
Mf
R
d

z

Fx Fx =Fy =0
EIz 0
Fx
Fx =Fy =0
/2
(
2
2
CR
CR
) CR2
/2
=
(1 cos ) d =
[ sin ]0 = 1
EIz 0
EIz
2 EIz

/2

Edef
1
Mfz

v2 =
=
Mfz
R d

Fy Fx =Fy =0
EIz 0
Fy
Fx =Fy =0
/2
2
2
CR
CR
CR2
/2
=
sin d =
[ cos ]0 =
EIz 0
EIz
EIz
ENER 22
Soit Iz le moment quadratique de la section droite. E est le module de Young du materiau.
Larc est soumis en 2 à une force de composantes (F, 0, 0).

On neglige linuence de leort normal et de leort tranchant.
1. Expression du moment echissant Mfz ().
Solution
D
eplacement horizontal de la section 2
Mfz () = RF (1 cos )
114
do`
/2
1
F 2 R3 /2
2
Edef =
Mfz |{z}
ds =
(1 cos )2 d
2 EIz 0
2 EIz 0
F 2 R3
=
2 EIz
(
=
R d
/2
[
]/2
F 2 R3 3 sin(2 )
(1 + cos 2 cos ) d =
+
2 sin
| {z }
2 EIz 2
4
0
2
3
2
4
1+cos(2 )
2
F 2 R3
2 EIz

Edef
u2 =
=
F
)
3
F R3
2
4
EIz
D
eplacement vertical et rotation de la section 2
Introduisons dans la section 2 une force auxiliaire de composantes (0, Fy , 0) et un couple auxiliaire de
composantes (0, 0, C).

Mfz () = C RF (1 cos ) + RFy sin
do`
u:
Mfz
= R sin ,
Fy
Mfz
=1
C
On en deduit le deplacement vertical de la section 2 :

/2

1
Mfz
Edef

=
Mf
R
d
v2 =

z

Fy Fy =C=0
EIz 0
Fy
Fy =C=0
F R3
EIz
/2
(sin cos sin ) d

| {z }
sin(2 )
2
[
]
F R3
cos(2 ) /2
F R3
=
cos +
=
EIz
4
2 EIz
0

/2

1
Edef
Mfz

=
2z =
Mf
R
d

z

C Fy =C=0
EIz 0
C
(
F R2
) F R2
/2
=
[ sin ]0 = 1
EIz
2 EIz
Fy =C=0
F R2
=
EIz
/2
(1 cos ) d
115
ENER 23
Larc (1 2) de centre O et de rayon moyen R et de section droite constante. Il est articule en 1 et
repose sur un appui simple de normale verticale en 2. Soit Iz le moment quadratique de la section
droite. E est le module de Young du materiau.
Larc est soumis en 2 à une force de composantes (F, 0, 0).

On neglige linuence de leort tranchant (modèle de Bernoulli).
1. Expression du moment echissant Mfz () et de leort normal N ().
4. Donner la position et la valeur du moment echissant maximal.
5. Calculer la rotation 2z de la section 2.
6. La section droite est un rectangle plein de base b et de hauteur h. On donne :
E = 210000 MPa , R = 100 mm , b = 21 mm , h = 4 mm , F = 200 N
Calculer u2 et 2z .
Solution
D
eplacement horizontal de la section 2
Les reactions dappui sont F1x , F1y et F2y .
Les equations dequilibre secrivent :

F1x + F = 0
F1y + F2y = 0 ,
RF + RF2y = 0
116
do`
u:
F1x = F
F1y = F
F2y = F
Le moment echissant et leort normal dans la section droite dabscisse sont :

Mfz () = RF sin + RF2y (1 cos ) = F R (sin + cos 1)
N () = F sin F2y cos = F (sin + cos )
do`
/2
/2
1
1
2
Mfz |{z}
ds +
N 2 ds
Edef =
2 EIz 0
2 EA 0
=
F 2 R3
2 EIz
R d
/2
F 2R
(sin + cos 1) d +
2 EA
(sin + cos )2 d
0
F 2 R /2
(2 + sin(2 ) 2 sin 2 cos ) d +
(1 + sin(2 )) d
2 EA 0
0
[
]/2
[
]
F 2 R3
cos(2 )
F 2R
cos(2 ) /2
=
2
+ 2 cos 2 sin
+
2 EIz
2
2 EA
2
0
0
(
)
2
3
2
F R
F R
= ( 3)
+
+1
2 EIz
2
2 EA
F 2 R3
=
2 EIz
/2
/2

u2 =
) FR
F R3 (
Edef
= ( 3)
+
+1
F
EIz
2
EA
Remarque : le moment echissant est maximal pour = /4 et :
Mfmax = ( 2 1) F R
Introduisons en 2 un couple auxiliaire C.
Laction de liaison F2y est alors egale à :

F2y = F
C
R
do`
u lexpression du moment echissant :
Mfz () = RF sin + RF2y (1 cos ) + C = F R (sin + cos 1) + C cos
117
et de leort normal :
N () = F sin F2y cos = F (sin + cos ) +
On en deduit :
Mfz
= cos
C
La rotation de la section 2 est egale à :
2z
C
cos
R
Mfz
cos
=
C
R

/2
/2

Edef
1
Mf
1
N

=
=
R
d
R
d
Mf
+
N

C C=0
EIz 0
C
EA 0
C
C=0
C=0
/2
F
F R2 /2
=
(sin cos + cos2 cos ) d +
(sin cos + cos2 ) d
| {z }
| {z }
EIz 0
EA 0
sin(2 )
2
1+cos(2 )
2
[
[
]
]/2
F R2
cos(2 ) sin(2 )
cos(2 ) sin(2 ) /2
F
=
+ +
sin
+
+ +
EIz
4
2
4
EA
4
2
4
0
0
(
)
(
)
2
1 FR
1
F
=
+
+
4 2 EIz
4 2 EA
Application num
erique
Avec RDM-Ossatures , larc etant discretise en elements de poutre droite de meme longueur, on
obtient :
nombre delements
3
5
10
20
40
solution analytique
u2 en mm
0.986
-18.31 %
1.124
-6.87 %
1.186
-1.73 %
1.202
-0.41 %
1.206
-0.08 %
1.207
2z en degres
1.206
-13.32 %
1.323
-4.91 %
1.374
-1.24 %
1.387
-0.31 %
1.390
-0.09 %
1.391
118
ENER 24
Soient A , Iz et k les caracteristiques de la section droite. E et sont respectivement le module de
Young et le coecient de Poisson du materiau.
Larc est soumis en 2 à une force horizontale de composantes (F, 0, 0).

On prend en compte linuence de leort normal et de leort tranchant.
1. Donner lexpression de leort normal N (), de leort tranchant T () et du moment echissant Mfz ().
4. La section droite est rectangulaire :
A = bh
Iz =
b h3
12
k=
5
6
On donne :
R = 60 mm
b = 10 mm
E = 200000 MPa ,
= 0.25
h = 6 mm
,
F = 30 N
Calculer u2 .
Solution
Les eorts dans larc sont :
N () = F sin ,
T () = F cos ,
Mfz () = RF sin
do`
3/2
3/2
3/2
1
1
1
2
2
N |{z}
ds +
T |{z}
ds +
Mfz2 |{z}
ds
Edef =
2 EA 0
2 GAk 0
2 EIz 0
=
F 2R
2 EA
R d
3/2
F 2R
sin2 d +
| {z }
2 GAk
1cos(2 )
2
soit :
Edef =
R d
3/2
0
F 2 R3
cos2 d +
| {z }
2 EIz
1+cos(2 )
2
3 F 2R 3 F 2R
3 F 2 R3
+
+
8 EA
8 GAk
8 EIz
R d
3/2
sin2 d
| {z }
1cos(2 )
2
119
On en deduit le deplacement horizontal de la section 2 :

u2 =
Edef
3 FR
3 FR
=
+
+
F
4
EA
GAk}
| {z }
| 4 {z
eort normal
eort tranchant
3 F R3
4 EI
| {z z}
moment
echissant
Application numerique :
u2 = 0.00035
| {z } + 0.00106
| {z } + 0.42412
| {z } = 0.42553 mm
0.08 %
0.25 %
99.67 %
Chapitre 4
M
ethode des
el
ements nis
4.1
4.1.1
Rappels
Energie
de d
eformation
Expression de l
energie de d
eformation en fonction des contraintes et des d
eformations
Si le comportement du materiau est elastique et lineaire, lenergie de deformation contenue dans le
volume dV est (formule de Clapeyron) :
1
(xx xx + yy yy + zz zz + xy xy + xz xz + yz yz ) dV
2
1
1
= {}T {} dV = {}T [D] {} dV
, dV = dx dy dz
2
2
dEdef =
(4.1.1)
avec les relations cinematiques :

u
v
w
, yy =
, zz =
x
y
z
u v
u w
=
+
, xz =
+
,
y
x
z
x
xx =
xy
yz
v w
=
+
z
y
(4.1.2)
et la relation de comportement :
{} = [D] {}
(4.1.3)
o`
u:
u(x, y, z; t), v(x, y, z; t) et w(x, y, z; t) sont les composantes du champ de deplacements.
}T
{
est le vecteur contrainte.
{}T = xx yy zz xy xz yz
{
}T
{}T = xx yy zz xy xz yz
est le vecteur deformation.
[D] est la matrice des coecients elastiques. Si le materiau est isotrope, elle ne depend que du
module de Young E et du coecient de Poisson (cours delasticite).
Remarque : lenergie de deformation sexprime en joule (1 J = 1 N.m = 1 kg.m2 .s2 ).
122
Energie
de d
eformation dun
el
ement
Lenergie de deformation dun element secrit :
Edef =
1
{u}T [k] {u}
2
(4.1.4)
o`
u {u} et [k] sont respectivement le vecteur des deplacements nodaux et la matrice de rigidite.
Application aux poutres droites
Lenergie de deformation accumulee dans un troncon de poutre, de longueur dx, soumis à un eort
normal N est egale à :
( )2
1 u
1
u
dEdef = N
dx = EA
dx
(4.1.5)
2 x
2
x
o`
u u(x) est le deplacement suivant x.
Lenergie de deformation accumulee dans un troncon de poutre, de longueur dx, soumis à un moment
de exion Mfz est egale à :
dEdef
1
z
1
= Mfz
dx = EIz
2
x
2
z
x
)2
dx
(4.1.6)
o`
u z (x) et v(x) sont respectivement la rotation suivant z et le deplacement suivant y. Si on adopte
lhypothèse de Bernoulli :
v
z =
(4.1.7)
x
lenergie de deformation secrit :
dEdef
1
= EIz
2
2v
x2
)2
dx
(4.1.8)
Lenergie de deformation accumulee dans un troncon de poutre, de longueur dx, soumis à un moment
de torsion Mt est egale à :
dEdef
1
x
1
= Mt
dx = GJ
2
x
2
x
x
)2
dx
(4.1.9)
o`
u x (x) est la rotation suivant x.
4.1.2
Energie
cin
etique
Expression g
en
erale
Lenergie cinetique contenue dans le volume dV est :
(( )
( )2 (
) )
)
1
u 2
v
w 2
1 (
dEcin =
+
+
dV = u 2 + v 2 + w 2 dV
2
t
t
t
2
(4.1.10)
o`
u u(x, y, z; t), v(x, y, z; t) et w(x, y, z; t) sont les composantes du deplacement du point de coordonnees (x, y, z) et la masse volumique du materiau.
123
Energie
cin
etique dun
el
ement
Lenergie cinetique dun element secrit :
Ecin =
1
{u}
T [m] {u}
(4.1.11)
o`
u {u}
et [m] sont respectivement le vecteur des vitesses nodales et la matrice de masse.
Application aux poutres droites
Lenergie cinetique dun troncon de poutre, de longueur dx, soumis à un eort normal N est egale à :
( )2
u
1
1
dx = A u 2 dx
dEcin = A
(4.1.12)
2
t
2
o`
u u(x; t) est le deplacement suivant x.
Lenergie cinetique de translation dun troncon de poutre, de longueur dx, soumis à un moment de
exion Mfz est egale à :
( )2
1
v
1
dEcin = A
dx = A v 2 dx
(4.1.13)
2
t
2
o`
u v(x; t) est le deplacement suivant y (èche).
4.1.3
Energie
potentielle et
el
ements nis
Lenergie potentielle est :

1
Epot = Edef Wext = {U }T [K]{U } {U }T {F } {U }T {Fnod }
{z
}
|2
{z
} |
Edef
(4.1.14)
Wext
o`
u Wext est le travail des forces exterieures pour le deplacement {U }.
Compte-tenu de la relation :
[K]{U } = {F } + {Fnod }
(4.1.15)
1
Epot = {U }T [K]{U } = Edef
2
(4.1.16)
Epot (elements nis) Epot (exacte)
(4.1.17)
lenergie potentielle secrit :
On montre que :
4.1.4
Modes propres
En labsence damortissement, les equations du mouvement libre dune structure sont :

L } + [KLL ]{UL } = 0
[MLL ]{U
(4.1.18)
o`
u {UL } est le vecteur des deplacements inconnus. En posant :
L } sin t
{UL (t)} = {U
(4.1.19)
L } est independant du temps, il vient :

o`
u {U
L } = 2 [MLL ]{U
L }
[KLL ]{U
(4.1.20)
124
L } le vecteur propre associe.

o`
u est une pulsation propre de la structure et {U
Les pulsations propres sont les solutions de lequation :
(
)
det [KLL ] 2 [MLL ] = 0
4.2
(4.1.21)
Exercices
4.2.1
Assemblage
Le vecteur deplacement {u(e) } dun element secrit en fonction du vecteur deplacement de la structure {U } :
{u(e) } = [L(e) ] {U }
Considerons, par exemple, la poutre droite soumise à un eort normal representee sur la gure cidessous.
Le vecteur deplacement {u(2) } de lelement (2) est :
u 1
{ } [
]
u 2
u
0
0
0
1
0
4
u3 = [L(2) ] {U }
{u(2) } =
=
u2
0 1 0 0 0
u4

u5

1
{U }T [K] {U }
2
(e)
=
Edef =
Edef =
el
ements
1
{U }T
2
el
ements
(
|
1 (e) T (e)
{u } [k ] {u(e) } =
2
el
ements
1
{U }T [L(e) ]T [k (e) ] [L(e) ] {U }
|
{z
}
2
[K (e) ]
[K (e) ] {U }
el
ements
{z
[K]
do`
u lexpression de la matrice de rigidite :
[K] =
[K (e) ]
el
ements
Dans lexemple ci-dessus, on obtient pour la matrice [K (2) ] :
0
0
[K (2) ] =
0
1
0
0 0
0
(2)
[
]
[
]
1 k (2) k (2)
0 k22
0
0
0
1
0
11
12
0
= 0 0
(2)
(2)
0
1
0
0
0
(2)
0 k21 k22
0 k12
0
0 0
0 0
(2)
0 k21
0 0
(2)
0 k11
0 0
0
0
125
et la matrice de rigidite [K] de la poutre est :
(1)
(1)
k11
0
0
k12
0
(2)
(3)
0 k (2) + k (3)
0
k
k
22
11
21
12
(4)
(4)
[K] =
0
k22
0
k21
0
(1)
(2)
(1)
(2)
k
k12
0 k22 + k11
0
21
(3)
(4)
(3)
(4)
0
k21
k12
0
k22 + k11
Lenergie cinetique est :
1
Ecin = {U }T [M ] {U }
2
(e)
=
Ecin =
el
ements
1 T
{U }
2
el
ements
(
|
1 (e) T (e)
{u } [m ] {u (e) } =
2
el
ements
1 T (e) T (e) (e)

{U } [L ] [m ] [L ] {U }
|
{z
}
2
[M (e) ]
[M (e) ] {U }
el
ements
{z
[M ]
do`
u lexpression de la matrice de masse :
[M ] =
[M (e) ]
el
ements
Le travail des forces exterieures appliquees aux elements pour le deplacement u(x) est :
Wext = {U }T {F }
(e)
Wext =
=
el
ements
el
ements
(
= {U }T
|
{u(e) }T {f (e) } =
{U }T [L(e) ]T {f (e) }
|
{z
}
el
ements
{F (e) }
el
ements
{z
{F }
do`
u lexpression du vecteur force :
{F } =
{F (e) }
el
ements
Dans lexemple ci-dessus, on obtient pour le vecteur {F (2) } :
0 0
(2)
0 1 { (2) }
f2
f
1
(2)
0
=
{F } = 0 0
(2)
(2)
1 0 f2
0 0
0
et le vecteur force {F } est :
(1)
(2)
(3)
f
+
f
2
(4)
{F } =
f2
(1)
(2)
f2 + f1
f (3) + f (4)
2
1
{F e }
126
4.2.2
ement de poutre droite soumis `

El
a un eort normal
On considère un element à deux nuds dont les coordonnees nodales sont x1 = 0 et x2 = L. A est
laire de la section droite. E et sont respectivement le module de Young et la masse volumique du
materiau.
Les deplacements nodaux sont u1 (t) et u2 (t).

Champ de d
eplacements et fonctions dinterpolation
Le deplacement suivant x des sections droites est de la forme :
u(x) = a + b x
o`
u a et b sont des constantes.
u(0) = u1 = a
u(L) = u2 = a + bL
on deduit :
u2 u1
L
do`
u lexpression du champ de deplacement en fonction des deplacements nodaux :
(x)
(
x)
u1 +
u2
u(x) = 1
L}
| {z L }
| {z
a = u1
b=
N1 (x)
soit :
N2 (x)
{ }
u1
u(x) = N1 (x) u1 + N2 (x) u2 = N1 (x) N2 (x)
= [N ] {u}
u2
[
o`
u N1 (x) et N2 (x) sont les fonctions dinterpolation ou fonctions de forme. Ces fonctions verient les
relations :
{
1 si i = j
Ni (xj ) = ij =
i, j = 1, 2
,
N1 (x) + N2 (x) = 1
0 si i = j
o`
u les xj sont les coordonnees nodales. Le champ de deplacements et les fonctions dinterpolation sont
representes sur la gure ci-dessous.
D
eformations et contraintes
127
On en deduit lallongement unitaire :

[
]
]
u
1 [
u2 u1
N1 N2
1 1 {u} = [B] {u} =
xx =
=
{u} =
x
L
L
x
x
et la contrainte normale :
u2 u1
L
dans lelement. Ces quantites sont constantes dans lelement (gure ci-dessous).
xx = E [B] {u} = E
Energie
de d
eformation et matrice de rigidit
e
Edef
1
=
2
EA
0
u
x
)2
1
dx =
2
EA 2xx dx
0
avec :
xx = [B] {u} = {u}T [B]T
do`
u:
Edef
1
= {u}T
2
)
EA [B] [B] dx {u}
{z
}
T
[k]
On en deduit lexpression de la matrice de rigidite :

[ ]
[
]
L
L
L
]
1 B1 [
1
1 1
T
B1 B2 dx =
EA 2
EA [B] [B] dx =
[k] =
EA 2
dx
L B2
L 1 1
0
0
0
[
soit :
k k
[k] =
k k
E
avec k = 2
L
A dx
0
Cas particuliers :
Si la section droite est constante, la matrice de rigidite est egale à :
[
]
EA 1 1
[k] =
L 1 1
Si la section droite est un rond plein dont le diamètre varie lineairement entre D1 et D2 , on a :
D2 ((
x)
x )2
A=
1
D1 + D2
4
L
L
do`
u:
L(
L( )
L(
x) x
D2 E
x 2
x )2
2
2
+D
+2
D
D
1
k =
D
dx
dx
dx
1
1 2
2
4 L2 1 0
L
L
L
L
{z
}
{z
}
{z
}
|
|0
|0
L/3
L/3
L/6
128

et :
[
]
) 1 1
E ( 2
2
[k] =
D1 + D1 D2 + D2
1 1
12 L
La valeur exacte de k obtenue à laide du theorème de Castigliano ( Methodes energetiques :

page 90) est :
E
k =
D1 D2
4L
La table ci-dessous donne lerreur (kmef kexact )/kexact en fonction de D2 /D1 .
D2 /D1
erreur
1
0
1.125
0.46 %
1.25
1.67 %
1.5
5.56 %
1.75
10.71 %
2
16.67 %
Remarque : en general, k est evalue numeriquement.
Energie
cin
etique et matrice de masse
Ecin =
1
2
A u 2 dx
0
avec :
u = [N ] {u}
= [N ]T {u}
T
do`
u:
Ecin
1
T
= {u}
2
)
A [N ] [N ] dx {u}
{z
}
T
[m]
On en deduit lexpression de la matrice de de masse :
[m] =
0
[ ]
]
N1 [
N1 N2 dx
A
N2
soit :
m11 m12
[m] =
m21 m22
avec :
A N12 dx =
m11 =
0
A N22
m22 =
m12 = m21 =
dx =
(
x )2
A 1
dx
L
A N1 N2 dx =
0
( x )2
L
dx
(
x) x
A 1
dx
L L
Remarques :
la somme des termes de la matrice de masse est egale à la masse de lelement :
me =
A dx
0
129
si la section droite est constante, la matrice de masse

[
me 2
[m] =
6 1
est :
]
1
2
Vecteur force d
u`
a une force r
epartie
Le travail de la force repartie dintensite lineique p pour le deplacement u(x) = [N ] {u} = {u}T [N ]T
est :
( L
)
L
T
T
Wext =
p u dx = {u}
p [N ] dx
0
0
|
{z
}
{f }
do`
{f } =
0
N
p 1
N2
{ }
f1
dx =
f2
Si la force est repartie lineairement entre p1 et p2 :

(
x)
x
p(x) = p1 1
+ p2
L
L
il vient :
L(
p N1 dx = p1
x )2
1
dx + p2
L
x(
L
x)
L
1
dx = p1 + p2
L
3
6
0
0
0 L
L
L(
L ( )2
x) x
x
L
L
f2 =
p N2 dx = p1
1
dx + p2
dx = p1 + p2
L L
L
6
3
0
0
0
f1 =
do`
u:
L
{f } =
6
2 p1 + p2
p1 + 2 p2
La discretisation de la force repartie est representee sur la gure ci-dessous.
4.2.3
Exercice : mise en
equation
La poutre daxe x represente sur la gure est discretisee en trois elements à deux nuds. Soient E
et le module de Young et la masse volumique du materiau. Laire de la section droite est egale à A
entre les nuds 1 et 2, 2 A entre les nuds 2 et 3 et à 3 A entre les nuds 3 et 4. La section 1 est
encastree.
130
La poutre porte entre 1 et 2 une force uniformement repartie dintensite lineique p.

La section 4 est soumise à une force F4x = F sin t.
Les conditions initiales sont u(x; t = 0) = u(x;
t = 0) = 0.
1.
2.
3.
4.
Eectuer lassemblage.
Construire lequation dont les deplacements inconnus sont les solutions.
Donner lexpression de la reaction dappui en fonction des deplacements inconnus.
Calculer les pulsations propres.
Solution
[
]
EA 1 1
[k12 ] =
,
2 L 1 1
[
]
[
]
2 EA 1 1
3 EA 1 1
[k23 ] =
, [k34 ] =
1 1
1 1
L
L
[
]
[
]
{ }
2 AL 2 1
3 AL 2 1
pL
[m12 ] = [m23 ] =
, [m34 ] =
, {f12 } =
1 2
1 2
pL
6
6
} + [K] {U } = {Fnod } + {F } :
Lassemblage conduit à la relation [M ] {U
=
0
u
=
0
4
2
1/2
1/2
1
1
EA
u
u2 =?
AL
2
4
+
4
2
1/2
1/2
+
2
2
2
+
u

u =?
2
4 + 6 3
2
2 + 3 3
6
L
3
3
u
4
u4 =?
3
6
3
3

F1x =?
pL
pL
F2x = 0
+
=
F3x = 0

F4x = F sin t
Les deplacements inconnus sont les solutions de lequation :

8 2 0 u
2
5/2 2 0 u2 pL
AL
EA
2 10 3 u
3 +
2 5 3 u3 =
0

6
L
0 3 6
u
4
3 3
u4
F sin t
avec les conditions initiales : u2 (0) = u3 (0) = u4 (0) = u 2 (0) = u 3 (0) = u 4 (0) = 0.
La reaction dappui est donnee par la relation :
AL
EA
u
2
u2 pL
3
2L
Les pulsations propres sont les solutions de lequation caracteristique :
5/2 2 0
8 2 0
EA
AL
2 5 3 2
2 10 3 = 0
det
L
6
3 3
0 3 6
F1x =
do`
u:
= 0.2780 , 1.4736 , 3.1711
E
L2
4.2.4
131
Exercice : mise en
equation
La poutre daxe x represente sur la gure est discretisee en deux elements à deux nuds. Soient E
et le module de Young et la masse volumique du materiau. Laire de la section droite est egale à A.
La poutre porte sur toute sa longueur une force uniformement repartie dintensite lineique p.
La section 1 est soumise à un deplacement impose u1 = a sin t.
Les conditions initiales sont u(x; t = 0) = u(x;
t = 0) = 0.
2. Construire lequation dont les deplacements inconnus sont les solutions.
3. Donner lexpression de la reaction dappui en fonction des deplacements inconnus.
Solution
[k12 ] = [k23 ] =
[
]
EA 1 1
L 1 1
[m12 ] = [m23 ] =
[
]
AL 2 1
1 2
6
{ }
pL 1
{f12 } = {f23 } =
2 1
} + [K] {U } = {Fnod } + {F } :
Lassemblage conduit à la relation [M ] {U
2
1
u
1 = a 2 sin t
1
1
u1 = a sin t
EA
AL
1 2 + 2 1
u
2
1 1 + 1 1
u2 =?
+
6
L
1
2
u
3
1
1
u3 =?
F1x =? pL 1
1+1
= F2x = 0 +
2
F3x = 0
1
Les deplacements inconnus u2 (t) et u3 (t) sont les solutions de lequation :
[
]{ }
[
]{ }
EA 2 1 u2
AL 4 1 u
2
+
1 2 u
3
u3
6
L 1 1
{ }
{ 2
}
{
}
pL 2
AL a sin t
EA a sin t
=
+
+
0
0
2 1
6
L
avec les conditions initiales : u2 (0) = u3 (0) = u 2 (0) = u 3 (0) = 0.
La reaction dappui est donnee par la relation :
F1x =
AL
EA
1
(2 a 2 sin t + u
2 ) +
(a sin t u2 ) pL
6
L
2
132
4.2.5
Exercice : contraintes et
energie de d
eformation
Considerons la poutre daxe x representee sur la gure. Soient E le module de Young du materiau
et A laire de la section droite supposee constante.
La section 1 est encastree. Elle est soumise à une force uniformement repartie sur toute sa longueur
dintensite lineique p.
1. Pour plusieurs discretisations, calculer les contraintes dans les elements et lenergie de deformation de la poutre.
2. Representer graphiquement la variation de lenergie de deformation en fonction de la discretisation.
Solution exacte
N = p (L x)
do`
L
L
1
p2
p2 L3
p2 L3
Edef =
N 2 dx =
(L x)2 dx =
= 0.167
2 EA 0
2 EA 0
6 EA
EA
Le champ de deplacements est solution de lequation :
du
N
= xx =
dx
EA
do`
u:
u(x) =
avec u(0) = 0
px
(2 L x)
2 EA
Solution
el
ements nis
Remarque : lenergie de deformation est egale à :
1
1 L
Edef =
xx xx dV =
A xx xx dx
2 V
2 0
1. La poutre est discretisee en un element à deux nuds. Les deplacements nodaux sont :
u1 = 0 ,
u2 =
1 pL2
2 EA
La deformation et la contrainte elementaires sont :

xx (1 2) =
1 pL
u2 u1
=
L
2 EA
xx (1 2) =
1 pL
2 A
133

Edef
p2 L2
=
2 EA
( )2
5
1 p2 L 3
p2 L 3
L=
= 0.125
6
8 EA
EA
(erreur : 25.00 %)
2. La poutre est discretisee en deux elements à deux nuds de meme longueur L/2. Les deplacements nodaux sont :
3 pL2
1 pL2
u1 = 0 , u2 =
, u3 =
8 EA
2 EA
Les deformations et les contraintes elementaires sont :
xx (1 2) =
u2 u1
3 pL
=
L/2
4 EA
xx (1 2) =
3 pL
4 A
u3 u2
1 pL
=
L/2
4 EA
xx (2 3) =
xx (2 3) =
1 pL
4 A

Edef
p2 L 2
=
2 EA
(( )
( )2 )
3 2
1
L
5 p2 L3
p2 L3
+
=
= 0.156
4
4
2
32 EA
EA
(erreur : 6.25 %)
3. La poutre est discretisee en trois elements à deux nuds de meme longueur L/3. Les deplacements nodaux sont :
u1 = 0
u2 =
5 pL2
18 EA
u3 =
4 pL2
9 EA
u4 =
1 pL2
2 EA

xx (12) =
u2 u1
5 pL
=
L/3
6 EA
xx (1 2) =
5 pL
6 A
xx (23) =
u3 u2
3 pL
=
L/3
6 EA
xx (2 3) =
3 pL
6 A
xx (34) =
xx (3 4) =
u4 u3
1 pL
=
L/3
6 EA
1 pL
6 A
134

(( )
( )2 ( )2 )
p2 L 2
5 2
3
1
L
35 p2 L3
p2 L3
Edef =
+
+
=
= 0.162
2 EA
6
6
6
3
216 EA
EA
(erreur : 2.78 %)
Conclusion
Le graphe ci-dessous represente la variation de lenergie de deformation en fonction de la discretisation.
On a :
Edef (1 element) < Edef (2 elements) < Edef (3 elements) < Edef (exacte)
4.2.6
Exercice : contraintes et
energie de d
eformation
Considerons la poutre daxe x representee sur la gure. Soient E le module de Young du materiau
et A laire de la section droite supposee constante.
Les deux extremites de la poutre sont encastrees. Elle est soumise à une force uniformement repartie
sur toute sa longueur dintensite lineique p.
1. Pour plusieurs discretisations, calculer les contraintes dans les elements et lenergie de deformation de la poutre.
2. Representer graphiquement la variation de lenergie de deformation en fonction de la discretisation.
135
Solution exacte.
N=
pL
px
2
do`
Edef
1
=
2 EA
N 2 dx =
0
p2 L3
p2 L 3
= 0.0.0417
24 EA
EA
Le champ de deplacements est solution de lequation :

du
N
= xx =
dx
EA
do`
u:
u(x) =
avec u(0) = 0
px
(L x)
2 EA
Solution
el
ements nis
Remarque : lenergie de deformation est egale à :
Edef =
1
2
xx xx dV =
V
1
2
A xx xx dx
0
1. La poutre est discretisee en deux elements à deux nuds de meme longueur L/2. Les deplacements nodaux sont :
1 pL2
u1 = 0 , u2 =
, u3 = 0
8 EA
xx (1 2) =
1 pL
u2 u1
=
L/2
4 EA
xx (1 2) =
1 pL
4 A
xx (2 3) = xx (1 2)
xx (2 3) = xx (1 2)

Edef =
p2 L 2
2
2 EA
( )2
1
L
1 p2 L3
p2 L 3
=
= 0.0313
4
2
32 EA
EA
(erreur : 25.00 %)
136
2. La poutre est discretisee en trois elements à deux nuds de meme longueur L/3. Les deplacements nodaux sont :
1 pL2
1 pL2
, u3 =
,
9 EA
9 EA
1 pL
u3 u2
u2 u1
=
xx (2 3) =
=0
xx (1 2) =
L/3
3 EA
L/3
u1 = 0
xx (1 2) =
1 pL
3 A
u2 =
xx (2 3) = 0 ,
u4 = 0
xx (3 4) = xx (1 2)
xx (3 4) = xx (1 2)

( )2
1
1 p2 L 3
p2 L3
p2 L2
L
2
=
= 0.0370
Edef =
2 EA
3
3
27 EA
EA
(erreur : 11.11 %)
3. La poutre est discretisee en quatre elements à deux nuds de meme longueur L/4. Les deplacements nodaux sont :
1 pL2
3 pL2
3 pL2
u3 =
u4 =
u6 = 0
32 EA
8 EA
32 EA
u2 u1
3 pL
u3 u2
1 pL
xx (1 2) =
=
xx (2 3) =
=
L/4
8 EA
L/4
8 EA
u1 = 0
xx (12) =
3 pL
8 A
u2 =
xx (3 4) = xx (2 3)
xx (4 5) = xx (1 2)
1 pL
xx (23) =
xx (34) = xx (23)
xx (45) = xx (12)
8 A

( ( )
( )2 )
3 2
1
L
5 p2 L 3
p2 L3
p2 L2
2
+2
=
= 0.0390
Edef =
2 EA
8
8
4
128 EA
EA
(erreur : 6.25 %)
137
Conclusion
Le graphe ci-dessous represente la variation de lenergie de deformation en fonction de la discretisation.
On a :
Edef (2 element) < Edef (3 elements) < Edef (4 elements) < Edef (exacte)
4.2.7
ement de poutre droite soumis `

El
a un eort normal
On considère lelement à trois nuds equidistants represente sur la gure. A est laire de la section
droite. E et sont respectivement le module de Young et la masse volumique du materiau.
Les deplacements nodaux sont : u1 , u2 et u3 .

Champ de d
Le deplacement suivant x des sections droites est de la forme :
u(x) = a + b x + c x2
o`
u a, b et c sont des constantes.
u(L/2) = u1 = a
1
1
bL + cL2
2
4
u(0) = u2 = a
u(L/2) = u3 = a +
on deduit :
1
1
bL + cL2
2
4
1
2
(u3 u1 ) , c = 2 (u1 + u3 2 u2 )
L
L
do`

[
] u1
u2 = [N ] {u}
u(x) = N1 (x) u1 + N2 (x) u2 + N3 (x) u3 = N1 (x) N2 (x) N3 (x)

u3
a = u2
b=
o`
u N1 (x), N2 (x) et N3 (x) sont les fonctions dinterpolation ou fonctions de forme :
[N ] =
1 [
x (2 x L) L2 4 x2
L2
x (2 x + L)
138
Ces fonctions verient les relations :

{
1 si i = j
Ni (xj ) = ij =
0 si i = j
i, j = 1, 2, 3
Ni (x) = 1
i=1
o`
u les xj sont les coordonnees nodales. Les fonctions dinterpolation sont representes sur la gure
ci-dessous.
D
On en deduit lallongement unitaire :
u
=
=
x
xx
avec :
N1
x
N2
x
N3
x
]
{u} = [B] {u}
[
]
]
1 [
[B] = B1 B2 B3 = 2 4 x L 8 x 4 x + L
L
do`
u:
xx =
u3 u1 4 x
+ 2 (u1 + u3 2 u2 )
L
L
Les deformations nodales sont :

2 = xx (0) =
u3 u1
L
avec :
1 = xx (L/2) = 2 d
3 = xx (L/2) = 2 + d
2
(u1 + u3 2 u2 )
L
= E [B] {u} est lineaire dans lelement (gure ci-dessous).
d=
La contrainte normale xx = E xx
Energie
de d
e
Edef
1
=
2
L/2
L/2
EA 2xx
1
dx = {u}T
2
L/2
)
EA [B] [B] dx {u}
T
L/2
{z
[k]
139
do`

B1 [
]
T
[k] =
EA [B] [B] dx =
EA B2 B1 B2 B3 dx
L/2
L/2
B3
L/2
L/2
Si la section droite est constante, on obtient :
7 8 1
EA
8 16 8
[k] =
3L
1 8 7
Energie
cin
etique et matrice de masse
Ecin
1
=
2
L/2
1
A u dx = {u}
T
2
L/2
L/2
)
A [N ] [N ] dx {u}
L/2
{z
[m]
do`
u lexpression de la matrice de de masse :

N1 [
]
[m] =
A N2 N1 N2 N3 dx
L/2
N3
L/2
4
2 1
me
2 16 2
[m] =
30
1 2
4
o`
u me = AL est la masse de lelement.
Remarque : la somme des termes de la matrice de masse est egale à me .
Vecteur force d
u`
a une force r
epartie
Le travail de de la force repartie dintensite lineique p pour le deplacement u(x) = [N ] {u} = {u}T [N ]T
est :
(
)
L/2
Wext =
L/2
L/2
p u dx = {u}T
p [N ]T dx
|
L/2
{z
{f }
do`

N1
{f } =
p N2 dx
L/2
N3
L/2
Si la force est uniformement repartie sur toute la longueur de lelement, on obtient :

1
pL
4
{f } =
6
1
140
La discretisation de cette force uniformement repartie est representee sur la gure ci-dessous.
4.2.8
Exercice : modes propres
La poutre representee sur la gure est encastree à ses deux extremites et sollicitee en tractioncompression.
A est laire de la section droite. E et sont respectivement le module de Young et la masse volumique
du materiau.
Etudier
les modes propres de la poutre avec les modelisations suivantes :
1. La poutre est discretisee en deux elements à deux nuds de longueur L/2.
2. La poutre est discretisee en trois elements à deux nuds de longueur L/3.
3. La poutre est discretisee en quatre elements à deux nuds de longueur L/4.
4. La poutre est discretisee en un element à trois nuds equidistants.
Solution analytique
Lequation dequilibre secrit :
A
u
x
(
)
u
EA
=0
x
avec les conditions aux limites : u(x = 0; t) = u(x = L; t) = 0.

La recherche de la solution harmonique :
u(x; t) = u
(x) sin(t)
conduit à lequation :
E
2u
+ 2 u
=0
x2
avec u
(x = 0) = u
(x = L) = 0
qui admet comme solutions :

(
x)
u
n (x) = a sin n
L
n = n
do`
u:
1 = 3.14
E
L2
n = 1, 2, . . .
E
L2
2 = 6.28
E
L2
3 = 9.42
E
L2
141
Solution
1. La poutre est discretisee en deux elements à deux nuds de longueur L/2.
[
]
2 EA 1 1
[k12 ] = [k23 ] =
1 1
L
[
]
AL 2 1
[m12 ] = [m23 ] =
12 1 2
La pulsation propre 1 est solution de lequation :

4
EA
1
12 AL = 0
L
3
do`
u:
1 = 3.46
E
L2
(erreur = 10.27 %)
La deformee modale est representee sur la gure ci-dessous.
2. La poutre est discretisee en trois elements à deux nuds de longueur L/3.

Les modes propres sont les solutions de lequation :
(
[
]
[
])
3 EA 2 1
2 AL 4 1
L } = {0}
{U
1 2
L
18 1 4
do`
u:
1 = 3.29
E
L2
(erreur = 4.61 %)
L1 } = 1
{U
2
(erreur = 16.95 %)
L2 } = 1
{U
2
2 = 7.35
E
L2
{ }
1
1
{
}
1
1
Les deformees modales sont representees sur la gure ci-dessous.
3. La poutre est discretisee en quatre elements à deux nuds de longueur L/4.
142

Les modes propres sont les solutions de lequation :
2 1
4 1
4
EA
AL
L } = {0}
1 2 1 2
1 4 1 {U
L
24
1 2
1 4
do`
u:
E
L2
1 = 3.22
(erreur = 2.59 %)
(erreur = 10.27 %)
(erreur = 19.46 %)
2 = 6.93
E
L2
3 = 11.26
E
L2

1
L1 } = 1
{U
2
2
1

1
L2 } = 1
0
{U
2 1

L3 } = 1 2
{U
2
1
Les deformees modales sont representees sur la gure ci-dessous.
4. La poutre est discretisee en un element à trois nuds equidistants.

7 8 1
EA
8 16 8
[k123 ] =
3L
1 8 7
4
2 1
AL
2 16 2
[m123 ] =
30
1 2
4
La pulsation propre 1 est solution de lequation :
do`
u:
16 EA
8
12
AL = 0
3 L
15
E
1 = 3.16
(erreur = 0.66 %)
L2
La deformee modale est (gure ci-dessous) :

u
1 (x) =
4x
(L x)
L2
4.2.9
143
ement ni de torsion
El
On considère lelement à deux nuds represente sur la gure. J est la constante de torsion et G est
le module delasticite transversal du materiau.
Les deplacements nodaux sont x1 et x2 .

Champ de d
La rotation des sections droites suivant x est de la forme :
x (x) = a + b x
o`
u a et b sont des constantes.
x (0) = x1 = a ,
x (L) = x2 = a + bL
on deduit :
x2 x1
L
do`
(
(x)
x)
x (x) = 1
x1 +
x2
L}
| {z L }
| {z
a = x1
b=
N1 (x)
soit :
x (x) = N1 (x) x1 + N2 (x) x2
N2 (x)
[
]
= N1 (x) N2 (x)
{ }
x1
= [N ] {u} = {u}T [N ]T
x2
Energie
de d
e
Edef =
avec :
[
x
N1
=
x
x
N2
x
1
2
GJ
0
do`
u:
Edef
1
= {u}T
2
{u} =
(
x
x
)2
dx
]
1 [
1 1 {u} = [B] {u}
L
)
GJ [B] [B] dx {u}
{z
}
T
[k]

[ ]
[
]
L
L
L
]
1 B1 [
1
1 1
T
B1 B2 dx =
[k] =
GJ [B] [B] dx =
GJ 2
GJ 2
dx
L B2
L 1 1
0
0
0
144

[
soit :
k k
[k] =
k k
G
avec k = 2
L
J dx
0
Cas particuliers :
Si la section droite est constante, la matrice de rigidite est egale à :
[
]
GJ 1 1
[k] =
L 1 1
Si la section droite est un rond plein dont le diamètre varie lineairement entre D1 et D2 , on a :
)4
((
x)
x
J=
1
D1 + D2
32
L
L
do`
u:
G
k =
32 L2
J dx =
0
)
G ( 4
D1 + D13 D2 + D12 D22 + D1 D23 + D24
160 L
[
]
1
1
[ k ] = k
1 1
et :
La valeur exacte de k obtenue à laide du theorème de Castigliano ( Methodes energetiques :

page 92) est :
3 G
D13 D23
k =
32 L D12 + D1 D2 + D22
La table ci-dessous donne lerreur (kmef kexact )/kexact en fonction de D2 /D1 .
D2 /D1
erreur
4.2.10
1
0
1.125
1.86 %
1.25
6.80 %
1.5
23.73 %
1.75
48.59 %
2
80.83 %
ement ni de exion : mod`

El
ele de Bernoulli
On considère un element à deux nuds de longueur L. E et sont respectivement le module de Young

et la masse volumique du materiau. Iz est le moment quadratique de la section droite.
Les deplacements nodaux sont v1 , z1 , v2 et z2 .

Champ de d
Le deplacement suivant y (èche) est de la forme :
v(x) = a + b x + c x2 + c x2 + d x3
145
o`
u a, b, c et d sont des constantes.
La rotation des sections droites est (hypothèse de Bernoulli) :
z (x) =
v
= b + 2 cx + 3 dx2
x

v1 = v(0) = a ,
z1 = z (0) = b
v2 = v(L) = a + bL + cL2 + dL3
z2 = z (L) = b + 2 cL + 3 dL2
on deduit :
a = v1
3
1
(v2 v1 ) (2 z1 + z2 )
2
L
L
Le champ de deplacement secrit donc :
c=
b = z1
d=
2
1
(v1 v2 ) + 2 (z1 + z2 )
3
L
L
[
v(x) = N1 (x) v1 + N2 (x) z1 + N3 (x) v2 + N4 (x) z2 = N1 N2 N3
v1
] z1
N4
= [N ] {u}
v
2
z2
o`
u les fonctions dinterpolations sont :
(
x2
x3
N1 = 1 3 2 + 2 3
L
L
N3 = 3
x2
x3
2
L2
L3
N2 = L
x2
x
x3
2 2 + 3
L
L
L
(
,
N4 = L
x3
x2
L3 L2
Remarque : le champ de deplacements :

v(x) = [N ] u
est le champ de deplacements exact pour un element de poutre à section constante et non charge ; en
eet, v(x) est solution de lequation :
d4 v
1
=
py = 0
4
dx
EIz
La rotation des sections droites est :

[
u
N1
=
z (x) =
x
x
N2
x
N3
x
]
N4
{u} = [N ] {u}
x
146
avec :
N1 =
N3
(
)
x
x2
x2
x
+ 2
, N2 = 1 4 + 3 2
L L
L
L
(
)
6 x
x2
x2
x
=
2
, N4 = 3 2 2
L L L
L
L
6
L
Energie
de d
e
Edef
avec :
[
z
N1
=
x
x
1
=
2
EIz
0
z
x
)2
dx
]
[
]
N3 N4
{u} = B1 B2 B3 B4 {u} = [B] {u}
x
x
)
)
(
2( x
6
x
, B2 =
3 2
B1 = 2 2 1
L
L
L
L
)
(
)
(
6
x
2
x
B3 = 2 1 2
, B4 =
3 1
L
L
L
L
N2
x
do`
u:
1
Edef = {u}T
2
(
|
)
EIz B T B dx {u}
{z
}
[k]
L
T
EIz [B] [B] dx =
[k] =
B1
B2 [
]
B1 B2 B3 B4 dx
EIz
B3
B4
Le coecient kij de cette matrice est egal à :
kij =
EIz Bi Bj dx
0
Pour une poutre à section constante :

le coecient k11 est egal à :
k11 = EIz
0
36 EIz
B1 B1 dx =
L4
L(
2
0
)2
12 EIz
x
1 dx =
L
L3
la matrice de rigidite est :
12
6 L 12 6 L
EIz 6 L 4 L2 6 L 2 L2
[k] = 3
L 12 6 L 12 6 L
6 L 2 L2 6 L 4 L2
147
Energie
cin
etique de translation et matrice de masse
La vitesse dune section droite suivant y est :
v = [N ] {u}
= {u}
T [N ]T
do`
u lexpression de lenergie cinetique de translation :
Ecin
1
=
2
1
A v dx = {u}
T
2
A [N ] [N ] dx {u}
{z
}
T
[m]
et celle de la matrice de masse :
N1
L
L
N2 [
]
[m] =
A [N ]T [N ] dx =
A
N3 N1 N2 N3 N4 dx
0
0
N4
156
22 L
54
13 L
AL
4 L2
13 L 3 L2
22 L
[m] =
54
13 L
156 22 L
420
13 L 3 L2 22 L 4 L2
Vecteur force d
u`
a une force transversale r
epartie sur l
el
ement
Le travail de la force transversale dintensite lineique p pour le deplacement v(x) = [N ] {u} =
{u}T [N ]T est :
L
L
T
Wext =
p v dx = {u}
p [N ]T dx
0
|0
{z
}
{f }
On en deduit lexpression du vecteur force :

N1
L
N2
{f } =
p
N3 dx
0
N4
Si la force est repartie uniformement sur toute la longueur de lelement, on obtient :

6
pL
L
{f } =
6
12

L
La discretisation de la force repartie est representee sur la gure ci-dessous.
148
Vecteur force d
u`
a un couple r
eparti sur l
el
ement
Le travail du couple reparti dintensite lineique m pour le deplacement z (x) = [N ] {u} = {u}T [N ]T
est :
Wext =
m z dx = {u}T
|
m [N ]T dx
{z
}
{f }
On en deduit lexpression du vecteur force :
N1
L
N2
{f } =
m
N3 dx
0
N4
Si le couple est reparti uniformement sur toute la longueur de lelement, on obtient :
0
{f } =
m
4.2.11
Exercice :
elasticit
e plane
On considère lelement triangle à trois nuds depaisseur t represente sur la gure ci-dessous.
Les deplacements nodaux sont : (u1 , v1 ) ,
(u2 , v2 ) ,
(u3 , v3 ).
1. Donner lexpression des fonctions dinterpolation.

2. Calculer les deformations et les contraintes.
3. Donner sans la calculer lexpression de la matrice de rigidite.
Solution
Fonctions dinterpolation
Le deplacement suivant x est de la forme :
u(x, y) = a + b x + c y
o`
u a, b et c sont des constantes.
149

u(0, 0) = u1 = a ,
u(L, 0) = u2 = a + bL ,
on deduit :
a = u1
b=
1
(u2 u1 ) ,
L
c=
u(0, H) = a + cH
1
(u3 u1 )
H
Le deplacement u(x, y) secrit donc :

1
1
u(x, y) = u1 + (u2 u1 ) x +
(u3 u1 ) y
L
H
)
(
y
x
x
y
= 1
u1 +
u2 +
u3
L
H
L
H
|
{z
}
|{z}
|{z}
N1 (x,y)
N2 (x,y)
N3 (x,y)
Remarque : les fonctions dinterpolation verient les relations :

{
3
1 si i = j
Ni (x, y) = 1
,
Ni (xj , xj )) = ij =
0 si i = j
i=1
i, j = 1, 2, 3
o`
u les (xj , yj ) sont les coordonnees nodales.
De meme, on a :
(
y)
x
y
x
v1 + v2 +
v3
v(x, y) = 1
L H
L
H
D
Les deformations sont :
u
u2 u1
xx =
=
x
L
soit :
yy =
L
xx
yy = 0

xy
1
0
1
H
1
v
v3 v1
=
y
H
1
L
0
0
1
L
1
H
xy =
u v
u3 u1 v2 v1
+
=
+
y
x
H
L

u 1

0
v1

1
u
2
v
H
u
0
3
v3
{} = [B] {u}
Les contraintes {} = [D] {} sont constantes dans lelement.
Energie
de d
e
(
)
1
1
T
T
T
{} [D] {} dV = {u}
[B] [D] [B] dV {u}
Edef =
|{z}
2 V | {z }
2
V
{u}T [B]T
[B] {u}
do`
[k] =
[B]T [D] [B] dV =
t [B]T [D] [B] dA
V
o`
u t et A sont respectivement lepaisseur et laire de lelement. Si lepaiseur de lelement est constante,
il vient :
[ k ] = [B]T [D] [B] A t

Exercice RDM

Transféré par

Droits d'auteur :

Formats disponibles

Exercice RDM

Transféré par

Informations du document

Description originale:

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Exercice RDM

Transféré par

Droits d'auteur :

Formats disponibles

R

esistance des mat

Table des mati`

Exercices de resistance des materiaux

4.1.3 Energie potentielle et elements nis . . . . . . . . .

Les deplacements et les deformations sont petits.

Tenseur des deformations :

Glissement en M dans les directions orthogonales na et nb :

Variation relative de volume :

Exercices de resistance des materiaux

Vecteur contrainte sur la facette n en M :

La contrainte normale sur la facette n est :

T = 2 max = max ( 1 , 2 , 3 ) min ( 1 , 2 , 3 ) E

Loi de comportement ou loi constitutive

Si le materiau est isotrope, la loi de comportement secrit :

Le tenseur des contraintes se reduit `a :

Les contraintes et les deformations principales sont :

Les crit`eres de Rankine, Von Mises et de Tresca se reduisent `a :

Exercices de resistance des materiaux

avec n2x + n2y = 1

Cette equation na de solution autre que la solution triviale nx = ny = 0 que si et seulement si :

et les valeurs propres :

(Sxx Syy )2 + 4 Sxy

Les vecteurs propres associes sont :

Remarque : les deux directions principales sont orthogonales :

Determinant dune matrice carree symetrique de dimension 3 :

S11 S12 S13

cos(1 + 2 ) = cos(1 ) cos(2 ) sin(1 ) sin(2 )

sin(1 + 2 ) = sin 1 cos 2 + cos 1 sin 2

Exercices de resistance des materiaux

ELA 1 : vecteur contrainte sur une facette

T (M, ) T (M, ) T (M, k)

sont representees sur la gure ci-dessous.

sont (formule de Cauchy : {T (M, n)} = [(M )] {n}) :

On en deduit la contrainte normale (n = n T (M, n)) :

7.412 + (45.19)2 + 41.482 = 61.78 MPa

Programme Scilab 5.3.1

{, , k}, le tenseur des contraintes a pour

Exercices de resistance des materiaux

T (M, ) T (M, ) T (M, k)

sont representees sur la gure ci-dessous.

sont (formule de Cauchy : {T (M, n)} = [(M )] {n}) :

On en deduit la contrainte normale (n = n T (M, n)) :

Contraintes et directions principales

Les contraintes principales 1 et 2 sont les solutions de lequation (det([ ] n [ I ]) = 0) :

La position angulaire des directions principales n1 et n2 est denie par :

Les contraintes et directions principales sont representees sur la gure ci-dessous.

Dans le rep`ere principal, le tenseur des contraintes secrit :

Exercices de resistance des materiaux

Programme Scilab 5.3.1

ELA 3 : contraintes et directions principales

{, , k}, le tenseur des contraintes a pour

Laxe k est direction principale :

Les contraintes principales 1 et 2 sont les solutions de lequation (det([ ] n [ I ]) = 0) :

La position angulaire des directions principales n1 et n2 est denie par :

Les contraintes et directions principales sont representees sur la gure ci-dessous.

Dans le rep`ere principal, le tenseur des contraintes secrit :

Exercices de resistance des materiaux