Chapitre 2

Partie des exercices avec solution

Exercice 1

a) Calcul la matrice de rigidité globale [K] de l'assemblage représenté sur la figure P3-1. Notez
que [K] devrait être en termes A1, A2; A3; E1; E2; E3; L1; L2 et L3.
b) A1 = A2 = A3 = A; E1 = E2 = E3 = E et L1 = L2 = L3 = L. Si les nœuds 1 et 4 sont fixés et
une force P agit au nœud 3 dans la direction positive x, trouve les déplacements des nœuds 2
et 3 en termes de A; E; L et P.

Figure. P3-1
a)

𝐴1 𝐸1 1 −1
[𝑘 (1) ]= [ ]
𝐿1 −1 1

𝐸2𝐴2 1 −1
[𝑘 (2) ]= [ ]
𝐿2 −1 1

𝐴3 𝐸3 1 −1
[𝑘 (3) ]= [ ]
𝐿3 −1 1

 A1 E1  A1 E 
 0 0 
 L1 L1 
  A1 E1 A1 E1 A2 E 2

 A2 E 2
0 

[k]=  L1 L1 L2 L2
 0 A2 E 2 A2 E 2 A2 E 2 AE 
  2 2
 L2 L2 L2 L2 
  A3 E3  A3 E3 
 0 0 
 L3 L3 

b)
𝐴1 𝐸1 𝐸2𝐴2 𝐴3 𝐸3 𝐴𝐸
= = =
𝐿1 𝐿2 𝐿3 𝐿

1 1 0 0
 1 2 1 0 

[𝑘]=  0  1
𝐴𝐸
1 2
𝐿
 
0 0 1 1

On sait que {F}=[k]{d}

 f 1x  ?  1 1 0 0  u1  0 
 f  0  0  u 2  ?
 2x  𝐴𝐸  1 2 1
→  = 𝐿  0  
 f 3 x p  1 2  1 u3  ? 
 
 f 4  ?  0 0 1 1  u 4  0

−𝐴𝐸 𝐴𝐸
→0= 𝑢2 - 𝑢3 → u3 = 2u2
𝐿 𝐿
−𝐴𝐸 2𝐴𝐸
P= 𝑢2 + 𝑢3
𝐿 𝐿
−𝐴𝐸 2𝐴𝐸
P= 𝑢2 + 𝑢3 (2u2 )
𝐿 𝐿

1 PL
u2 =
3 AE
1 PL 2 PL
u3 = 2 → u3 =
3 AE 3 AE
Exercice 2

Pour l’assemblage de barres représenté sur la figure P3-2, déterminer les déplacements
nodaux, les forces dans chaque élément et les réactions (utilisez la méthode de rigidité directe
pour ces problèmes).
 Figure P3-2.

Solution

Élément 1-2

[k1−2 ]=84*106 [ 1 −1]

−1 1
Élément 2-3

[k 2−3 ]=84*106 [ 1 −1]

−1 1

{f}=[k]{d} et u1 = 0

f1x =? 1 −1 0 u1 = 0
6
{ f2x = 0 }=84*10 [−1 2 −1] { u2 =? }
f3x = −500 0 −1 1 u3 =?

→ 2𝑢2 − u3 = 0 → u3 = 2u2

→ -5000= 84*106 [u2 + u3 ]

En remplaçant (1) en (2), nous avons

−5000
=-u2 + 2u2 → u2 = −0.595 ∗ 10−4 m
84∗106

u3 = −1.19 ∗ 10−4 m
Element 1-3
(1)
f 1 −1 0 𝑓1𝑥 = 5000𝑁
{ 1x }=84*106 [ ]{ −4 }→ (1)
f2x −1 1 −0.595 ∗ 10 𝑓2𝑥 = −5000𝑁
Element 2-3
(2)
f 1 −1 −0.595 ∗ 10−4 𝑓2𝑥 = 5000𝑁
{ 2x }=84*106 [ ]{ −4 }→ (2)
f3x −1 1 −1.19 ∗ 10 𝑓3𝑥 = −5000𝑁

𝐹1x = 84 ∗ 106 [1 −1 −0.595 ∗ 10−4 }→ 𝐹 = 5000N

0] { 1x
−1.19 ∗ 10−4

Exercice 3

Pour l’assemblage de barres représenté sur la figure P3-3, déterminer les déplacements
nodaux, les forces dans chaque élément et les réactions (utilisez la méthode de rigidité directe
pour ces problèmes).

Figure P3-3.
Solution
Elément 1-2

[k1−2 ]=3*106 [ 1 −1]

−1 1
Elément 2-3

[k 2−3 ]=1.2*106 [ 1 −1]

−1 1
Matrice de rigidité globale
3 −3 0
[k]=106 [−3 3 + 1.2 −1.2]
0 −1.2 1.2
Déplacement
f1x 3 −3 0 u1 = 0
6
= 8000}=10 [−3 3 + 1.2 −1.2] { u2 }
{f2x
f3x 0 −1.2 1.2 u3 = 0
u2 =1.905*10−3 mm
 f1x 3 −3 0 0
6
{ 2x }=10 [−3 3 + 1.2 −1.2] {1.905 ∗ 10−3 }
f
f3x 0 −1.2 1.2 0
f1x=−5715 N
→ f2x=−8000 N
f3x=−2286 N
(1)
f 1 −1 0 𝑓1𝑥 = −5715𝑁
{ 2x }=3*106 [ ]{ −3 }→ (1)
f3x −1 1 −1.905 ∗ 10 𝑓2𝑥 = −5715𝑁
(2)
f 1 −1 1.905 ∗ 10−3 𝑓2𝑥 = 2286 𝑁
{ 2x }=1.2*106 [ ]{ }→ (2)
f3x −1 1 0 𝑓3𝑥 = −2286 𝑁
Exercice 4

Pour l’assemblage de barres représenté sur la figure P3-4, déterminer les déplacements nodaux,
les forces dans chaque élément et les réactions (utilisez la méthode de rigidité directe pour ces
problèmes).

Figure P3-4
Solution

[k1−2 ]=[k 3−4 ]=1.5 ∗ 106 [ 1 −1

]
−1 1

[k 2−3 ]=5000[ 1 −1]

−1 1

f 1x  0  1500  1500 0 0  u1  0 

f2  0   1500 1505 5 0  u 2  ? 
=103   
f 3 x  10000  0 5 1505  1500   u 3  ? 
 
{ f 4 x  0 }  0 0  1500 1500  u 4  0
 0 1505 −5 u2 u2=2..21∗10−5 mm
{ }=103 [ ] {u }→
10000 −5 1505 3 u3 = 6.65 ∗ 10−3 𝑚𝑚
Les réactions

F1x =(1500*103 )(u3 ) → F1x = −33.15 N

F4x =(1500*103 )(u3 ) → F4x = −9975 N

(1)
f 1 −1 0 f1x = −33.15 N
{ 1x } = 1.5 x 106 [ ]{ } ⇒
f2x −1 1 2.21 x 10−5 (1)
f2x = 33.15 N

(2)
f 1 −1 2.21 x 10−5 f = −33.15 N
{ 2x } = 5000 [ ]{ } ⇒ 2x(2)
f3x −1 1 6.65 x 10−3 f3x = 33.15 N

(3)
f 1 −1 6.65 x 10−3 f3x = 9975 N
{ 3x } = 1.5 x 106 [ ]{ } ⇒ (3)
f4x −1 1 0 f4x = −9975 N

Exercice 5

Pour l’assemblage de barres représenté sur la figure P3-5, déterminer les déplacements nodaux,
les forces dans chaque élément et les réactions (utilisez la méthode de rigidité directe pour ces
problèmes).

Figure P3-5.
Solution
1 −1
[k (1) ] = [k (2) ] = 7000 [ ]
−1 1
1 −1
[k (2) ] = 2000 [ ]
−1 1
{F} = [k]{d}
𝐹1𝑋 =? 7 −1 0 0 𝑢1 = 0
𝐹 = 8𝑘𝑁
{ 2𝑥 } = 103 [−7 14 −7 0 ] {𝑢2 = ? }
𝐹3𝑥 = 0 0 −7 9 −2 𝑢3 =?
𝐹4𝑥 = ? 0 0 −2 2 𝑢4 = 0
⇒ 8 = 103 [14u2 − 7u3 ]

0 = 103 [−7u2 + 9u3 ]

7
⇒ u3 = u2
9

En remplaçant (2) en (1)

8 7
⇒ = 14 u2 − 7 x u2
103 9

⇒ u2 = 0.9351 x 10−3 m

⇒ u3 = 0.7273 x 10−3 m

Element (1)

f 1 −1 0 f1x −6.546
{ 1x } = 7 x 103 [ ]{ } ⇒ { }= { } kN
f2x −1 1 0.9351 x 10 −3 f2x 6.546
Element (2)

f 1 −1 0.9351 x 10−3 f2x 1.455

{ 2x } = 7 x 103 [ ]{ } ⇒ { } = { } kN
f3x −1 1 0.7273 x 10−3 f3x −1.455
Element (3)

f 1 −1 0.7273 x 10−3 f 1.455

{ 3x } = 2 x 103 [ ]{ } ⇒ { 3x } = { } kN
f4x −1 1 0 f4x −1.455
0
F1x = 103 [7 −7] { } = F1x = −6.546 kN
0.9351 x 10−3
−3
F4x = 103 [−2 2] {0.7273 x 10 } = F4x = −1.455 kN
0

Exercice 6

Pour le Treille supportée par le ressort au nœud 1 de la figure P3-33 (a), déterminer les
déplacements nodaux et les contraintes dans chaque élément. Soit E= 210 GPa et A = 5 104 m2
pour les deux éléments de treillis.
Solution

Elément 1-2 ; θ=135° C2=0.5, CS=-0.5, S2=0.5

0.5 −0.5 −0.5 0.5
(210 𝑋109 )(5𝑋10−4 ) −0.5 0.5 0.5 −0.5]
[𝑘1−2 ] = [
5 −0.5 0.5 0.5 −0.5
0.5 −0.5 −0.5 0.5

1 −1 −1 1
[𝑘1−2 ] = 105𝑋105 [−1 1 1 −1]
−1 1 1 −1
1 −1 −1 1

Elément 1-3 ; θ=180° C2=1, CS=0, S2=0

1 0 −1 0
(210 𝑋109 )(5𝑋10−4 )
[𝑘1−3 ] = [0 0 0 0]
10 −1 0 1 0
0 0 0 0
 1 0 −1 0
[𝑘1−3 ] = 105𝑋10 [ 05 0 0 0]
−1 0 1 0
0 0 0 0

Elément 1-2 ; θ=270° C2=0, CS=0, S2=1

0 0 0 0
[𝑘1−4 ] = 20 𝑋105 [ 0 1 0 −1]
0 0 0 0
0 −1 0 1

{𝐹} = [𝐾]{𝑑}
Les conditions aux limites
u2= v2= u3= v3= u4= v4 = 0
La matrice finale
𝐹1𝑥 = 0 210 −105 𝑢1
{𝐹 = −50} = 105 [ ]{ }
1𝑦 −105 125 𝑣1

⇒0 = 210 u1 – 105 v1 v1 = 2 u1

– 50000 = 105 [– 105 u1 + 125 (2 u1)]

⇒u1 = – 3.448 10–3 m

⇒v1 = – 6.896 10–3 m

−3.448𝑋10−3
210𝑋109 −3
𝜎1−2 = [0.707 −0.707 −0.707 0.707] {−6.896𝑋10 }
5𝑚 0
0

𝜎1−2 = 102.4𝑀𝑃𝑎 (𝑇)

9
−3.448𝑋10−3
210𝑋10 −3
𝜎1−3 = [1 0 −1 0] {−6.896𝑋10 }
10𝑚 0
0
 𝜎1−2 = −72.4 𝑀𝑃𝑎 (𝐶)

𝐾𝑁
𝐹𝑠 = (200 ) (−6.896𝑋10−3 𝑚)
𝑚

𝐹𝑠 = 13.792 𝐾𝑁

𝑓1−3 = 35.6 𝐾𝑁

𝑓1−2 = 51.2 𝐾𝑁

∑ 𝐹𝑦 = 0 = −50 + 13.79 + 36.198 = 0

Chapitre 3

Partie des exercices avec solution

Exercice 1
Pour la poutre représentée sur la figure P4-3, déterminer la rotation au support A et le
déplacement sous la charge P. Déterminer les réactions. Dessinez les schémas de force de
cisaillement et de flexion. EI est constant tout la long de la poutre.

Figure P4-3.

Solution

𝐿
1. On pose =𝑙
2

Elément 1-2

12 6𝐿 −12 6𝑙
𝐸𝑇 6𝐿 2 2
[𝑘1−2 ] = 3 [−12 4𝑙 −6𝑙 −2𝑙 ] ;
𝑙 −6𝐿 12 −6𝑙
2
6𝐿 2𝑙 −6𝑙 4𝑙 2
Elément 2-3
 12 6𝐿 −12 6𝑙
𝐸𝑇 6𝐿 2 2
[𝑘2−3 ] = 3 [−12 4𝑙 −6𝑙 −2𝑙 ] ;
𝑙 −6𝐿 12 −6𝑙
2 −6𝑙
6𝐿 2𝑙 4𝑙 2

F1y =? 0 0 v1=0
12 6l −12 6l
M1 = 0 6l 4l2 −6l 0 0 0 θ1 =?
F2y = −P EI −12 −6l 24 0 −12 6l v2=?
M2 = 0
=3
l 6l 2l2 0 8l2 −6l 2l2 θ2=? ?
F2y =? 0 0 −12 −6l 12 −6l v3 = 0
{ M3 =? } [ 0 0 6l 2l2 −6l 4l2 ] {θ3 = 0}

0
EI
4l2 −6l 2l2 θ1
{−P}=L3 [ 6l 24 0 ] { v2 }
0 2l2 0 8l2 θ2

0
ET
4l2 2l2 −6l θ1
{0}= L3 [ 2l2 8l2 0 ] { θ2 }
P −6l 0 24 v2

N = K ββ − K βα K −1
αα K αβ

2
= 3 [24 − [−6l 0] [4l2
EI 2l2 ] {−6l}] = 13.7148 EI
l 2l 8l2 0 l3

l3 −7l3 −7Pl3
dβ = N −1 F ⟹ v2 = (−P) = ⇒ v2 =
13.7148EI 96EI 768EI

θ1 2
2l2 ] {−6l} [−7Pl3] ⟹ {θ1 } =
{dα } = [K −1
αα ][K αβ ]{dβ } ⟹ {dα } = { } = [4l2
θ2 2l 8l2 0 96EI θ2
−pl2 −pl2

[ 8EI
pl2 ] = [ 32EI
pl2
]
32EI 8EI

En remplaçant l'équation dans la matrice globale, nous avons

 0
F1y 12 6l −12 6l 0 0 −Pl2
M1 6l 4l 2
−6l 2l2 −12 0 8EI
−7Pl3
F2y EI −12 −6l 24 0 −6l 6l
=3 96EI
M2 l 6l 2l2 0 8l2 12− 2l Pl2
F3y 0 0 −12 −6l 12 −6l 32EI
{ M3 } [ 0 0 6l 2l2 −6l 4l2 ] 0
{ 0 }
EI −6Pl3 7Pl3 6Pl3 EI 10Pl3 5P
F1y = 3
[ + + ]= ⇒ F1y =
l 8EI 8EI 32EI l3 32EI 16

De même
11P
M1 = 0 F3y =
16
−3PL
F2y = −P M3 =
16

M2 = 0

Exercice 2
Pour la poutre représentée sur la figure P4-7, déterminer la rotation au support A et la rotation
et le déplacement sous la charge P. Déterminer les réactions. Dessinez les schémas de force de
cisaillement et de flexion. EI est constant tout la long de la poutre.
 Figure P4-7

Solution

Elément 1-2

12 6l −12 6l
EI 4l2 −6l 2l2
[K1−2 ]= l3 [ ]
12 −6l
symmetry 4l2

Elément 2-3

12 6l −12 6l
EI 4l2 −6l 2l2
[K 2−3 ]= l3 [ ]
12 −6l
symmetry 4l2

Condition aux limite

v1 = θ1 = v3 = 0

EI 0 0 v2
[k] = [ ]{ }
l3 0 8l2 θ2

-0.0032142=24v2 ⟹ v2 = −1.34 ∗ 10−4 m

0.0064285=72θ2 ⟹ θ2 = 8.93 ∗ 10−5 rad

F1y 24 12l 0
−24 6l 0 0
M1 12l 8l2 2
−12l 4l 0 0 0
F2y (210 ∗ 109 )(4 ∗ 10−4 ) −12 −12l 24 −6l −12 6l −1.34 ∗ 10−4
=
M2 3 24 4l2 −6l 12l2 −6l 2l2 8.93 ∗ 10−5
F3y 0 0 −12 −6l 12 −6l 0
{ M3 } [ 0 0 6l 2l −6l 4l2 ] { 0 }

⇒ ℱ2𝑦 = 3.1 ∗ 106 (−12(1.34 ∗ 10−4 ) + 6(3)(8.93 ∗ 10−5 )) ⇒ ℱ1𝑦 = 10000𝑁

𝑀1 = 3.1 ∗ 106 (−6(3)(1.34 ∗ 10−4 ) + 2(3)2 (8.93 ∗ 10−5 )) ⇒ 𝑀1 = 12500𝑁. 𝑚

De même

𝐹2𝑦 = −10000𝑁

𝑀2 = 20000𝑁. 𝑚

𝐹3𝑦 = 1.87𝑁

𝑀3 = −2500𝑁. 𝑚

Element 1-2

F1y
12 6l −12 6l 0
m1 (210∗109) )(4∗10−4 ) 2 2 0
= [ 6l 4l −6l 2l ] { }
F2y (3)3 −12 −6l 12 −6l 1.34 ∗ 10−4
6l 2l2 −6l 4l2 8.93 ∗ 10−5
{m2 }

⇒ F1y = 10000N

M1 = 12500n. m

F2y = −10000N

M2 = 17500N. m

Element 2-3

ℱ2y
12 6l −12 6l 1.34 ∗ 10−4
m2 (210∗109) )(4∗10−4 ) 2
2l2 ] {8.93 ∗ 10−5 }
= [ 6l 4l −6l
ℱ3y (3)3 −12 −6l 12 −6l 0
6l 2l2 −6l 4l2 0
{m3 }

⇒ F2y = 1.87N

M2 = 2500n. m

F3y = 1.87N

M3 = −2500N. m
 Diagramme de cisaillement

Diagramme de moment

Exercice 3
Pour la poutre représentée sur la figure P4-10, déterminer la rotation au support A et la rotation
et le déplacement sous la charge P. Déterminer les réactions. Dessinez les schémas de force de
cisaillement et de flexion. EI est constant tout la long de la poutre.
 Figure P4-10.
Solution
Application de la symétrie

12 6l −12 6l
ET
4l2 6l 2l2
[K] = 12+KL2
L3 −6l
EI
[symmetry 4l2 ]

Application les conditions aux limites 𝑣1 = 0, 𝜃2 = 0

M1 = 0 (70∗109 )(2∗10−4 ) 4l2 6l θ1

{F = −6000N} = )[ 12+KL2 ] { }
2y 43 −6l v2
EI

6 6 6
⇒ 0 = 4l2 θ1 − 6lv2 ⇒ θ1 = 4l v2 ⇒ θ1 = 16 v2 − 6000 = 218750[−24 (16) v2 +
12.457v2 ] ⇒ v2=−7.9338∗10−3 m

6
θ1 = 16 (−7.9338 ∗ 10−3 ) ⇒ θ1 = −2.9752 ∗ 10−3 rad
 ℱ1y 12 6l −12 6l 2
M1 9 −4
(70∗10 )(2∗10 ) 6l 4l2 6l 2l2 −2.9752 ∗ 10−3
= ) [ ]{ }
ℱ2y 4 3
−12 −6l 12.457 −6l −79336 ∗ 10−3
{ M2 } 6l 2l2 −6l 4l2 0

F1y = 5.208KN ↑, M2 = 20.83KN. m ↺

F2y = 0KN ↓

Fspring = 1.587KN

Suivent la symétrie F3y = 5.208KN ↑

Exercice 4
Déterminez la flèche au point milieu de la poutre, les réactions et tracez les diagrammes de
force de cisaillement et de flexion pour la poutre fixé aux extrémités et soumis à une charge
uniformément répartie w représentée sur la figure P4-13. EI est constante. Comparez vos
réponses avec la solution classique.

Figure P4-13

Solution

Application de la condition aux limites

𝑣1 = ∅1 = 𝑣3 = ∅3 = 0

−wl ET 24 0 v2 ET 2
{ } = l3 [ 2 ] {∅ } ⇒ 0 = l3 (8l ∅2 ) ⇒ ∅2 = 0
0 0 8l 2
 ET −wl4
−wl = (24v2 ) ⇒ v2 =
l3 24ET

(6) wl (6) wl2 (6) (6)

F1y = , m1 = , F2y = −wl, m2 = 0
2 4

(6) wl −wl2
F2y = , m3 =
2 4

Ceux-ci sont obtenus à partir de l'équation matricielle suivante

(6)
ℱ1y
(6) 12 6l −12 6l 0 0 0
M1 0
4l2 −6l 2l2 0 0
(6) −wl4
ℱ2 ET 24 0 −12 6l 24ET
= )
M2
(6) l3 8l2 −6l 2l2 0
(6) 12 6l 0
ℱ3y [ 4l2 ] { 0 }
(6)
{M3 }
{F} = [K]{d} − {F0 }

wl −wl wl
F1y = −( ) = wl =
2 2 2

wl2 −wl2 wl2 wl2

M1 = −( )= =
4 12 3 12

F2y = −wl − (−wl) = 0, M2 = 0 − 0 = 0

wl −wl wl wl2 −wl2 wl2
F3y = ( )= , M3 = = =
2 2 2 12 12 12
1
wl4 w( )4 −wl4
2
v2 = 24ET = ⇒ v2 = 384ET
24ET
Chapitre 5

Partie des exercices avec solution

Exercice 1
Évaluez la matrice de rigidité pour les éléments représentés sur la figure P6-6. Les coordonnées
sont données en unités de millimètres. Supposons les conditions de contraintes planes. Soit E =
210 GPa, ν = 0. 25 et t = 10 mm.

Figure P6-6

Solution

a)
 βi = yj − ym = 30 − 120 = −90 γi = xm − xi = 50 − 80 = −30
βj = ym − yj = 120 − 30 = 90 γi = xi − xm = 20 − 50 = −30

βi = yi − yj = 120 − 30 = 90 γi = xj − xi = 80 − 20 = 60

2A = xi (yj − ym ) + xj (ym − yi ) + xm (yi − yj )

= 20(−90) + 80(90) + 50(0) = 4500 mm2

−90 0 90 0 0 0
1
[B] = [ 0 −30 0 −30 0 60]
5400
−30 −90 −30 90 60 0

1 0.25 0 1 0.25 0
210×109 11
[D] = [0.25 1 0 ] = 2.24 × 10 [0.25 1 0 ]
1−(0.25)2
0 0 0.375 0 0 0.375

[k] = tA[B]T [D][B]

−90 0 −30
0 −30 −90
5.4 × 10−3 1 90 0 −30 (2.24
[k] = (0.01) ( )( )
2 5.4 × 10−3 0 −30 90
0 0 60
[ 0 60 0 ]
1 025 0
× 1011 ) [0.25 1 0 ] [ B]
0 0 0.375
 −90 −22.5 −11.25
−7.5 −30 −33.75
−90 0 −30 0 −30 90
[k] = 1.12 × 109 90 22.5 −11.25 1
[ 0 −30 −90 0 0 60 ]
−7.5 −30 33.75 5.4×10−3
90 0 −30 0 60 0
0 0 22.5
[ 15 60 0 ]

8437.5 1687.5 −7762.5 −337.5 675 −1350

1687.5 3937.5 337.5 −2137.5 −2025 −1800
[k] = 2.074 × 105 −7762.5 337.5 8437.5 −1687.5 −675 1350
−337.5 −2137.5 −1687.5 3937.5 2025 −1800
−675 −2025 −675 2025 1350 0
[ −1350 −1800 1350 −1800 0 3600 ]

(b)

βi = −5 γi = 0

βj = 2.5 γj = −5

βm = 2.5 γm = 5

25.0 0 − 12.5 6.25 − 12.5 − 6.25

9.375 9.375 − 4.6875 − 9.375 − 4.6875
[k] = 4.48 × 107 15.625 − 7.8125 − 3.125 − 1.5625
27.343 1.5625 − 3.125
15.625 7.8125
[ 27.343 ]

t= 0.01

tA
A= A = 5 × 10−5
2

[k] = tA[B]T [D][B]

 [k]
1.225 × 109 3.5 × 108 −1.015 × 109 −7 × 107 −2.1 × 108 −2.8 × 108
3.5 × 108 7 × 108 7 × 107 −1.4 × 108 −4.2 × 108 −5.6 × 108
9 7
= −1.015 × 107 7 × 10 1.225 × 109 −3.5 × 108 −2.1 × 108 2.8 × 108
−7 × 10 −1.4 × 108 −3.5 × 108 7 × 108 4.2 × 108 −5.6 × 108
−2.1 × 108 −4.2 × 108 −2.1 × 108 4.2 × 108
4.2 × 108 0
[ −2.8 × 108 −5.6 × 108 2.8 × 108 −5.6 × 108
0 1.12 × 109 ]

Exercice 2
Pour les éléments de déformation plans représentés sur la figure P6-10, les déplacements
nodaux sont donnés comme suit : u1 = 0. 005 mm, v1 = 0. 002 mm, u2 = 0. 0 mm, v2 = 0. 0
mm, u3 = 0. 005 mm, v3 = 0. 0 mm.
Déterminez les contraintes σx , σy , τxy ,σ1 et σ2 et l'angle principal θp. E = 70 GPa et ν=0. 3, et
utilisez l'épaisseur de l'unité pour la déformation plane. Toutes les coordonnées sont en
millimètres.

Figure P6-10.
Solution
(a)
βi = −15 γi = −10

βj = 15 γj = 0

βm = 0 γm = 10

2A = 150 mm2 = 150 × 10−6 m2

σx 0.7 03 0
70 × 106
σ
{ y} = [0.3 0.7 0 ]
τxy (1 + 0.3)[1 − 2(0.3)]
0 0 0.2
5.2
2.0
1 −15 0 15 0 0 0
× [ 0 −10 0 0 0 10] 0 × 10−3
150 0
−10 −15 0 15 10 0
5.0
{0}

σx −52.5
σ
{ y } = [−32.8] MPa
τxy −5.38

−52.5 + (−32.8) −52.5 + 32.8 2

σ1.2 = √
± ( ) +(−5.38)2 = −42.65 ± 11.22
2 2

σ1 = −31.4 MPa σ2 = −53.9 MPa

2(−5.38)
2θp = tan−1 = 28.64°
−52.5— (−32.8)

θp = 14.32°

(b)
 0.005
σx 49.4 40.2 0 0 0 0 0.002
1 −15 0 15
σ
{ y } = 109 [ 94.0 0 ] × [ 0 15] 0
0 0 0 −15
τxy 225 × 10−6 0
26.9 0 −15 −15 15 15 0
0.005
{ 0 }
× 10−3

σx −31.4
σ
{ y } = {−13.5} MPa
τxy 5.38

−31.4 − 13.5 −31.4 + 13.5 2

σ1.2 = √
± ( ) +(−5.38)2 = −22.44 ± 10.46
2 2

σ1 = −11.98 MPa σ2 = −32.9 MPa

2(−5.38)
tan 2θp = = −15.5°
−31.4 + 13.46

(c)

σx 0.7 03 0 1
70
{ σy } = ( ) [0.3 0.7 0 ]
τxy 1.3 × 0.4 200
0 0 0.2

5
−10 0 10 0 0 0 2
×[ 0 −10 0 −10 0 20] 0 × 10−3
0
−10 −10 −10 10 20 0
5
{0}

σx −27.6
{ σy } = {−19.5} MPa
τxy 4.04
 −27.6 − 19.5 −27.6 + 19.5 2
σ1.2 = ± √( ) +4.042
2 2

σ1 = −17.9 MPa σ2 = −29.3 MPa

1 2 × 4.04
θp = tan−1 ( )
2 −27.6 + 19.5

θp = −22.5°

(d)

E = 70 × 109 v=3 t=1

γ1 = 0.015 γ2 = 0.005 γ3 = 0.015

x1 = 0.005 x2 = 0.015 x3 = 0.025

β1 = y2 − y3 β2 = y3 − y1 β3 = y1 − y2

γ1 = x3 − x2 γ2 = x1 − x3 γ3 = x2 − x1

1
A= [x (y − y3 ) + x2 (y3 − y1 ) + x3 (y1 − y2 )]
2 1 2

βi 0 βj 0 βm 0 1−v v 0
1 E v 1−v 0 ]
[B] = [ 0 γi 0 γj 0 γm ] [D] = [
2A (1+v)(1−2v) 1−2v
γi βi γj β j γm βm 0 0 2

[k] = tA[B]T [D][B]

u1 = 0.000005 u2 = 0 u3 = 0.000005

v1 = 0.000002 v2 = 0 v3 = 0
 u1
v1
u
{d} = v2
2
u3
{v3 }

{σ} = [D][B]{d}

σx 0.7 03 0 1
70
{ σy } = ( ) [0.3 0.7 0 ]
τxy 1.3 × 0.4 200
0 0 0.2

5
0 0 10 0 −10 0 2
×[ 0 −10 0 10 0 10 ] 0 × 10−3
0
−10 0 10 10 10 −10
5
{0}

σx −1.05
σ
{ } = { 0.70 } MPa
y
τxy 3.5

σx +σy σx −σy 2 σx +σy σx −σy 2

σ1 = + √( ) + τ2xy σ2 = + √( ) + τ2xy
2 2 2 2

−1.05 − 0.70 −1.05 + 0.07 2

σ1.2 = √
± ( ) +3.52
2 2

σ1 = −3.43 MPa σ2 = −3.78 MPa

1 2(−3.5)
θp = tan−1 ( )
2 −1.05 + 0.70

θp = −38°
Exercice 3

Déterminer les déplacements nodaux et les contraintes de l'élément, y compris les contraintes
principales, pour la plaque mince de la section 6.5 avec une charge de cisaillement uniforme
(au lieu d'une charge de traction) agissant sur le bord droit, comme le montre la figure P6-13.
Utilisez E = 30 106 N /mm2, ν = 0:30 et t = 1 mm.

Figure P6-14.

Solution

Reportez-vous à l’exercice 6.2 pour [K]

Avec u1=v1=0, u2=v2=0

0 48 0 − 28 14 u3
−5,000 75000 × 5 87 12 − 80 v3
{ }= [ 48 − 26 ] {u4 }
0 0.91
−5,000 Symmetry 87 v4
Résolution

u3 = 0.50 × 10−3 in. v3 = −0.275 × 10−2 in.

u4 = −0.609 × 10−3 in. v4 = −0.293 × 10−2 in.

Élément (1)

{σ} = [D][B]{d}

30 × 106 1 0.3 0 0 0 10 0 −10 0

{σ} = [0.3 1 0 ] [ 0 −20 0 0 0 20 ]
(0.91)(200)
0 0 0.35 −20 0 0 10 20 −10
0
0
0.5 × 10−3
×
−0.275 × 10−2
0
{ 0 }

σx 824
σ
{ y } = { 247 } psi
τxy −1586

824 + 247 824 + 247 2

σ1.2 = √
± ( ) +(−1586)2
2 2

σ1 = 2149 psi σ2 = −1077 psi

1 −2 × 1586
θp = tan−1 ( )
2 824 − 247

θp = −40°

Élément (2)

30 × 106 1 0.3 0 −10 0 10 0 0 0

{σ} = [0.3 1 0 ][ 0 0 0 −20 0 20]
(0.91)(200)
0 0 0.35 0 −10 −20 10 20 0
0
0
−0.609 × 10−3
×
−0.293 × 10−2
0.5 × 10−3
{−0.275 × 10−2 }
 σx −825
{ σy } = { 292 } psi
τxy −411

−825 + 292 −825 + 292 2

σ1.2 = √
± ( ) +(−411)2
2 2

σ1 = 426 psi σ2 = −960 psi

1 −2 × 411
θp = tan−1 ( )
2 −825 − 292

θp = −18.15°