Elementsfini 150517232857 Lva1 App6891

SEMESTRE 4
2006-2007
MODULE F412 :
MÉTHODES ÉNERGÉTIQUES
PROBLÈMES et CORRIGES
G.LHERMET-G.VESSIERE
TD-2A-S4-F412
SOMMAIRE
1. MÉTHODE DE CLAPEYRON ............................................................................................................................... 3

2. MÉTHODE DE MAXWELL-MOHR ET DE CASTIGLIANO.............................................................................. 17
3. MÉTHODE DES ÉLÉMENTS FINIS ................................................................................................................... 59
4. PROBLEMES DE SYNTHESE ............................................................................................................................. 81
ANNEXES................................................................................................................................................................ 90
ANNEXE A1. INTEGRALES DE MOHR............................................................................................................. 91
ANNEXE A2. MOMENTS FLÉCHISSANTS DANS QUELQUES CAS SIMPLES : ............................................ 98
IUTB-LYON1-GMP-DDS Page 2
TD-2A-S4-F412
1. MÉTHODE DE CLAPEYRON
PROBLÈME N°1
C Les barres AB et BC (de même section droite A et de même

matériau) sont articulées en A, B et C. Une force verticale
P est appliquée en B.
1°) Déterminez l'expression littérale de l'énergie de
déformation élastique emmagasinée par les deux barres en
fonction de P, L, E et A.
2°) En déduire l'expression littérale du déplacement
vertical du point B.
45°
B 3°) Peut-on avec cette méthode, calculer la déplacement
A
horizontal de B ?
P Application numérique : P=10kN, L=2m, E=200GPa et
A=1cm2. Calculez la valeur (en J) de l'énergie de
L
déformation élastique et la valeur du déplacement vertical
du point B en mm).
Réponses : (1 + 2 2 )P 2 L (1 + 2 2 )PL
Wdef = = 19,14J •VB = = 3.83 mm
2EA EA
RÈPONSES N°1
1 )
Nous devons pour appliquer la méthode de Clapeyron :
C
Wext = Wdéf calculer l’énergie de déformation emmagasinée
dans les deux barres, et l’identifier avec le travail
de la force extérieure
1 1
Wext = PD Wext = P δ vB
2 2
N X2 TY2 TZ2 M X2 M Y2 M Z2
∑∫
L
Wdéf = ( + + + + + )dx
barres
0 2EA 2GA 2GA 2GIG 2EIY 2EI Z 45° B
Compte tenu que les deux barres sont articulées en A, B, et C A
elles ne sont soumises qu’à des efforts normaux. P
N 2
N ( N )
AB 2
( BC
) 2
∑∫ )dx = ∫ ( )dx + ∫ (
L B C L
Wdéf = ( X X
)dx X
barres
0 2EA A 2EA B 2EA
TD-2A-S4-F412
Pour trouver les efforts normaux s’exercant dans les barres AB et CD, mettons en équilibre
le nœud B :
F2
Y
 2 45°
B
− F1 − F2 =0 F2 = P 2 F1
 2 ⇒ 
− P + F 2
=0  F1 = − P X P
 2
2
Après application du principe d’action mutuelle

N XAB = −P N XBC = P 2 P 2
C
d’où :
Wdéf = ∫
L (− P)
2
(P 2 ) )dx
2
∫
L 2
( )dx + (
0 2EA 0 2EA L 2
P2 L P2 L
2EA ∫0 EA ∫0
2
Wdéf = dx + dx
P2 P2 L
[x]0 + [x]0
P 2
= B
L 2
Wdéf A
2EA EA P P
L
Wdéf =
P2
2EA
L+
P2
EA
L 2 Wdéf =
P2 L
2EA
1+ 2 2 ( )
2 )
Appliquons la méthode de Clapeyron : Wext = Wdéf
Wext =
1
2
P δ vB W déf =
P2L
2 EA
1+ 2 2 ( )
En identifiant : δ vB =
PL
EA
(
1+ 2 2 )
3 )
Non, car l’identification implique qu’une seule charge.
Application numérique : P=30kN, L=2m, E=200GPa et A=1cm2.
( ) ( ) ( )
2
P2 L 10.103 2.103
Wdéf = 1+ 2 2 = 1 + 2 2 = 19142Nmm Wdéf=19,142 J
2EA 2 × 200.103 ×100
δ vB =
PL
EA
(
1+ 2 2 = )
10.103 × 2.103
200.103 × 100
(
1 + 2 2 = 3,828mm ) •vB=3,83 mm
TD-2A-S4-F412
ROBLÈME N°2
n2A 1°) Déterminez l'expression de l'énergie de déformation élastique en

A fonction de P, L, A, E et n.
P 2°) En déduire l'expression du raccourcissement de la barre. Déterminer
B
C l'expression du raccourcissement de la barre pour n=1 et n=2.
A
L/2 L/2
Réponses:
P 2L  1  PL  1  PL 5PL
Wdef = 1 + 2  •L = 1 + 2  •L = (n = 1) •L = (n = 2 )
4EA  n  2EA  n  EA 8EA
RÈPONSES N°2
Z
n2A
Y
A
G x
A P
B
C
L/2
L/2
1)
La poutre ABC est soumise à de la compression pure : N xAB = − P N xBC = − P
Wdéf = ∫
(N ) )dx + ( (N )
AB 2 BC 2
∫
B C
x x
( 2
)dx
A 2En A B 2EA
Wdéf = ∫
L
(N ) )dx + ( (N )
AB 2 BC 2 L
(− P)
2
(− P) 2
∫ )dx = ∫ ∫L (
L L
2
( x
2 L
x 2
( )dx + )dx
0 2En A 2 2EA 0 2En2 A 2 2EA
L L
P2 P2 L P2 P2 P2 L P2 L
[x]0 + [x]L = 2 +
2En2 A ∫0 2EA ∫2
= dx + L dx =
2 2 L
Wdéf
2En2 A 2EA 2 4En A 4EA
TD-2A-S4-F412
P2L  1
Wdéf = 1 + 2 
4 EA  n 
2)
Le déplacement associé à la charge P en C (•hC), est le raccourcissement de la barre •L :
1 1
Wext = PD Wext = P δ hC
2 2
1 P2 L  1 PL  1 
P δ hC = 1 + 2  δ hC = 1 + 2 
2 4 EA  n  2 EA  n 
Ce déplacement étant positif, il s’effectue dans le même sens que la charge . C’est donc bien
un raccourcissement de la barre.
PL  1
∆L = 1 + 2 
2 EA  n 
PL 5PL
Si n=1 ∆L = Si n=2 ∆L =
EA 8EA
PROBLÈME N°3
1°) Déterminez l’expression de l’énergie de

A
d 2 déformation élastique en fonction de C, L, G et d. En
déduire l’expression de la rotation de la section
L
d
droite C.
B 2°) L’énergie de déformation élastique admissible est
C C L/2 de 12J, L=1m, d=40mm, G=80 GPa. Calculez la valeur
du couple maximum admissible C (en Nm).et la valeur
de la rotation de la section C (en °).
12C 2 L 24CL
Réponses : W def = • = C = 802.1 Nm θC = 1.7°
G πd 4 Gπd 4
TD-2A-S4-F412
RÈPONSES N°3
d•2
A
Z
d x
G B
L
C Y
C L/2
1°)
Nous devons pour appliquer la méthode de Clapeyron : Wext = Wdéf ,calculer l’énergie de
déformation emmagasinée dans la poutre ABC, et l’identifier avec le travail exterieur du au
couple C.
Le déplacement associé au couple C en C (•xC), est la rotation suivant l’axe du couple de la
section C :
1 1
Wext = PD Wext = Cθ xC
2 2
∑∫
L
Wdéf = ( + + + + + )dx
barres
0 2EA 2GA 2GA 2GIG 2EIY 2EI Z
La poutre ABC est soumise à de la torsion pure : M xCB = −C M xBA = −C
Wdéf =∫ (
M x2
A B M
)dx = ∫ (
CB
x
A M
)dx + ∫ ( x(
BA 2
)2
)dx
( )
C 2GI CA C 2GI CB B 2GI BA
G G G
Wdéf = ∫
L
(− C)
2 3L
(− C ) 2
C 2 L2 C 2 32L
0
2
(
2GI GCB
)dx + ∫
L
2
2
(
2GI GBA
)dx =
2GI CB ∫0 dx + 2GI BA ∫L2 dx
G G
TD-2A-S4-F412
C2L C 2 L C 2 L  1 1  ( )4
Wdéf = + = + =
πd 4
et I GBA =
πd 2
2GI GBA 2G  2I GCB I GBA 
CB
avec I G
4GI CB
G
32 32
12C 2 L
Wdéf =
Gπd 4
12C 2 L
Wext
1
= Cθ xC et Wdéf = , identifions les 2 expressions :
2 Gπd 4
24CL
θ xC =
Gπd 4
2°) Application numérique : Wdéf=12J, L=1m, d=40mm, G=80 GPa
1
 Gπd Wdéf 
 = 80.10 × π × 40 ×12.10 = 802121Nmm
2 4 2 3 4 3
12C L 
Wdéf = ⇒ C =
Gπd 4  12L  12 ×1000
 
C=802,1 Nm
24CL 24 × 802121× 1000

θ xC = = = 0,0299rd •xC= 1°71
Gπd 4 80.103 × π × 404
PROBLÈME N°4
D
900
Φd = 26 P
A C0 En ne tenant compte que de l’énergie de déformation
P élastique due au moment de torsion et sachant que C0 =
63 Nm, G = 80 GPa, calculez :
Φd0 = 18 1°) L'énergie de déformation élastique (en J).
C
2°) La rotation de A (en d°).
600
Réponses : Wdef = 5.35 J ϕA = 9.72 °
B
P
56
20
TD-2A-S4-F412
RÈPONSES N°4
D
Ød=26mm
900mm A
C0=63Nm
56mm
600mm 20mm
56
C = C0 ×
C Ød0=18mm 20
B
C C0
C
56mm B
1°) 20mm
56
Les arbres AB et CD sont soumis à de la torsion pure : M xAB = C0 M xCD = C = C0 ×
20
Wdéf == ∫
B (M ) )dx + ( (M )
AB 2 CD 2
∫
D
x x
( AB CD
)dx
A 2GI G C 2GI G
C0
2
C2 C 02 600 C 2 900
Wdéf = ∫
600
∫0 2GI CD
900
AB ∫0 CD ∫0
( ) dx + ( ) dx = dx + dx
0 2GI GAB G
2GI G
2GI G
300C 02 450C 2 300C 0 450× 56 × C 0 C 0  300 450× 562 

2 2 2 2
Wdéf = + = + = +
GI GAB GI CD GI AB
202
GI CD
G  I AB 202 I CD 
G G G  G G 
Wdéf = + 
(
+
)
C 02  300 450× 562  63.103 2  300× 32 450× 562 × 32 
 =  = 5345Nmm

G  IG AB
20 I G  π × 80.103  184
2 CD 
202 × 264 
Wdéf=5,35 J
A
Le déplacement associé au couple C0 en A (• x ), est la rotation suivant l’axe du couple de la
section A :
1 1
Wext = P D dans notre cas : Wext = C0θ xA
2 2
Méthode de Clapeyron : Wext = Wdéf d’où : 1 C0θ xA = 5345Nmm ⇒ θ xA = 2 × 5345 = 0,170rd
2 63.103
•xA= 9°72
TD-2A-S4-F412
PROBLÈME N°5
L Déterminez l'expression littérale du

Z x Z P déplacement vertical de A en prenant en
G Y B x A compte le moment de flexion et l'effort
h
G tranchant.
b Montrez que si ν = 0.3 et L/h = 10
l'erreur commise en négligeant l'effort
tranchant est inférieure à 1%.
Réponses : PL3 PL PL3 1 + • h2 PL3 0.65 erreur 0.65%

•VA = + = (1 + 2
) = (1 + )
3EIY GA 3EIY 2 L 3EIY 100
RÈPONSES N°5
x Z
P
Y b
B
G x h
A
Détermination du travail extérieur :
Le déplacement associé à la charge P en A (•vA), est le déplacement vertical de la section A :
1 1
Wext = PD Wext = Pδ vA
2 2
Détermination de l’énergie de déformation :
∑∫
L
Wdéf = ( + + + + + )dx
barres
0 2EA 2GA 2GA 2GIG 2EIY 2EI Z
TD-2A-S4-F412
Nous devons chercher le Torseur de Section en G de la barre AB :

 0 0 
 
ℑ sec tionG =  0 P(L − X )
− P 
G 
0
La barre est donc soumise à un effort tranchant constantde B à A : TZ = − P ,
et à un moment de flexion linéaire : M Y = P (L − X )
d’où :
TZ2 M Y2 TZ2 M Y2
Wdéf = ∫ )dx = ∫ ( )dx + ∫ (
L L L
( + )dx = WdéfTZ + WdéfM Y
0 2GA 2EIY 0 2GA 0 2 EI
Y
Calculons l’énergie de déformation due à l’effort tranchant :
TZ2 2
P2 L P2 P2 L
[x]0 =
L P
Wdéf = ∫ )dx = ∫
2GA ∫0
L
dx = dx =
TZ L
(
0 2GA 0 2GA 2GA 2GA
Puis l’énergie de déformation due au moment fléchissant:
L (P(L − x )) P  (L − x ) 
2 3 L
MY2 P2 P2 L3
WdéfMY = ∫ )dx = ∫ ∫ (L − x) dx =
L L
dx =  =
2
( 
0 2EIY 0 2EIY 2EIY 0 2EIY  3  0 6EIY
D’où :
P 2 L P2 L3
Wdéf = Wdéf + Wdéf
TZ MY
= +
2GA 6EIY
Application numérique : si ν = 0.3 et L/h = 10
P2 L P2 L3 P 2 L  IY L2 
Wdéf = + =  +  G=
E
2GA 6EIY 2IY  GA 3E 
or
2(1 + υ )
 (1 + υ )  h 2  P 2 L3  (1 + 0,3)
(0,1)2 
P 2 L3 1 +
Wdéf =
    Wdéf = 1 +
6EIY  2 L  6EIY  2 
Wdéf =
P2 L3
6EIY
( −3
1 + 6,5.10 =
P2 L3
6EIY
) −2
1 + 0,65.10 = (
P 2 L3  0,65 
1 +
6EIY  100 
 )
Conclusion :
P 2 L3 P2 L3 0,65
WdéfM Y = et Wdéf =
TZ
×
6EIY 6EIY 100
TD-2A-S4-F412
L’énergie de déformation due à l’effort tranchant, ne représente que moins de 1 % de

l’énergie de déformation totale.
Nous en concluons que l’énergie de déformation due à l’effort tranchant est en

règle générale négligeable par rapport à l’énergie de déformation due au moment
de flexion.
Expression littérale du déplacement vertical de A
P L  0,65 2 3
PL3  0,65 
= Pδ vA et Wdéf = 1 +
1
Wext  En identifiant : δ =
A
1 + 
2 6EIY  100  v
3EI Y  100 
La flèche due à l’effort tranchant ne représente que 0,65% du déplacement

vertical de la section A.
PROBLÈME N°6
L
Déterminez en ne tenant compte que du moment
a b
C Z fléchissant l'expression littérale de la rotation de la
A x B
section droite C en fonction de C, L, a, b, E et IY.
C G
Réponses : C(a 3 + b 3 )
•C =
3EIY L2
TD-2A-S4-F412
RÈPONSES N°6
a b
C Z Y
G x
C
A B
L
Détermination du Torseur Externe: Voir module F112

a b
C C/L
C/L
A C B
L
Détermination du Torseur de Section: Voir module F112 (Annexes)

L
a b
C/L C/L
C Z Y
G x
A C B
TZ
C
TZ =
C/L L
+
A C B
MY
C
MY = X Ca/L
L
+ B
A C
-
Cb/L MY = −
C
(L − X )
L
Le déplacement associé au couple C en C (•YC), est la rotation suivant l’axe du couple de la

section C :
1 1
Wext = PD Wext = CθYC
2 2
TD-2A-S4-F412
2
N X2 TY2 TZ2 M X2 M Y2 M Z2 B T M Y2
∑∫ )dx Wdéf = ∫ (
L
Wdéf = ( + + + + + Z
+ )dx
barres
0 2EA 2GA 2GA 2GIG 2EIY 2EI Z A 2GA 2EIY
Nous devons négliger l’énergie de déformation due à l’effort tranchant par rapport celle
MY2
=∫
B
emmagasinée par la flexion : Wdéf dx
A 2 EI
Y
2 2
C   C 
2 2  X  − ( L − X )
dx = ∫   dx + B  L  dx
C M B M C L
Wdéf =∫ Y
dx + ∫ Y
∫
A 2 EI C 2 EI A 2 EI C 2EIY
Y Y Y
C2  a 2
X dx + ∫ (L − X ) dx  =
C 2   X 3  a  (L − X )3  L  C 2 a3 + b3
  ( )
2  ∫0
L
= +  =
2
2  
Wdéf
2EIY L  a  2EIY L  3 0  3  6 EI L2
 a Y
C a +b 2
( 3 3
) (
C a 3 + b3 )
= CθYC et Wdéf =
1
Wext En identifiant : θYC =
2 6EIY L2 3EIY L2
PROBLÈME N°7
Déterminez en ne tenant compte que du moment fléchissant

L l'expression littérale de la flèche verticale de la section droite C en
P fonction de P, L, a, b, E et IY.
a b
Z
A C x B
G
Réponses : Pa 2 b 2
•VC =
3EIY L
TD-2A-S4-F412
RÈPONSES N°7
a P b
Z Y
G x
C
A B
L
Détermination du Torseur Externe: Voir module F112
a b
P Pa/L
Pb/L
A C B
L
Détermination du Torseur de Section: Voir module F112 (Annexes)
a b
Pb/L Pa/L
P Z Y
G x
A C B
TZ Pb
TZ =
L
Pb/L
+
B
A C - TZ = −
Pa Pa/L
MY L
Pab/L
MY =
Pb
X MY =
Pa
(L − X )
L + L
B
A C
Le déplacement associé à la charge P en C (•vC), est le déplacement vertical de la section A :
1 1
Wext = PD Wext = Pδ vC
2 2
TD-2A-S4-F412
2
N X2 TY2 TZ2 M X2 M Y2 M Z2 B T M Y2
∑∫ )dx Wdéf = ∫ (
L
Wdéf = ( + + + + + Z
+ )dx
barres
0 2EA 2GA 2GA 2GIG 2EIY 2EI Z A 2GA 2EIY
Nous devons négliger l’énergie de déformation due à l’effort tranchant par rapport celle
MY2
=∫
B
emmagasinée par la flexion : Wdéf dx
A 2 EI
Y
2 2
 Pb   Pa 
2 2  X  ( L − X )
dx = ∫ 
L 
dx + ∫   dx
C M B M L
=∫ dx + ∫
C B
Y Y
Wdéf
A 2 EI C 2 EI A 2EIY C 2EIY
Y Y
P2  2  X 3 
3 L
2 (L − X )
  
a
P2  2 a 2 
( )
2EIY L2  ∫0 ∫a
L
=  + −  = +
2 2
 
 2EIY L2   3  0
Wdéf b X dx a L X dx b a
 3 a 
 
Wdéf =
( =
)
P2 b2a3 + a2b3 P2b2a 2 (a + b) P2b2a 2 L
= Wdéf =
P 2b2a2
6EIY L2 6EIY L2 6EIY L2 6EIY L
P 2b 2 a 2
Wext
1
= Pδ vC et W déf = En identifiant : δ =
C Pb 2 a 2
2 6 EI Y L v
3EI Y L
TD-2A-S4-F412
2. MÉTHODE DE MAXWELL-MOHR ET DE
CASTIGLIANO
PROBLÈME N°8
D=30
d=20 En appliquant la méthode de MAXWELL-MOHR
calculer les déplacements de A et de B.
C Q B P On donne : P = 30 kN, Q= 150 kN, E= 70 GPa.
A
Réponses : • A = −1,12 mm • B = 1,94 mm

800 600
RÈPONSES N°8
D=30
d=20
C Q
B P
A
800
600
Considérons l’état initial E0, et utilisons la notation charges généralisées Pi et déplacements

généralisés associés Di :
Etat initial
E0
D=30
d=20
C P2=Q
B P1=P
A
800
600
TD-2A-S4-F412
Le déplacement généralisé D1 associé à la charge généralisée P1 est le déplacement linéique

dans la direction de P1 de la section A.
Le déplacement généralisé D2 associé à la charge généralisée P2 est le déplacement linéique

dans la direction de P2 de la section B.
[D] = [f][P]  D1 
ou   = 
D
 f11 f12   P1 
f22  P2 
 2   f21
D1 = f11P1 + f12 P2
D2 = f21P1 + f22 P2
La matrice carrée [f] est la matrice de flexibilité de la structure.
Chaque coefficient d’influence fij est calculé par le théorème de Maxwell-Mohr
(nX )i (nX )j (tY )i (tY )j (tZ )i (tZ )j (mX )i (mX )j (mY )i (mY )j (mZ )i (mZ )j
∑∫
L
fij = ( + + + + )dx
barres
0 EA GA GA GIG EIY EIZ
Pour les déterminer nous devons décomposer l’état initial E0, en deux états unitaires E1 et E2
Etat unitaire
E1
Etat initial
E0
C
B P1=1
A
C P2 P1
B Etat unitaire
= A
E2
C
P2=1
E0 • P1×E1+P2×E2
B
A
La poutre ABC, dans les 3 états E0, E1 et E2 n’est soumise qu’à des efforts normaux :
(nX )i (nX ) j
fij = ∫
A
dx
C EA
A (nX )i (nX ) j B (nX )i (nX ) j A (nX )i (nX ) j
fij = ∫ dx = ∫ dx + ∫ dx
C EA C EACB B EABA
TD-2A-S4-F412
Pour calculer les fij, nous devons chercher les diagrammes des efforts normaux dans les deux
états unitaires :
Etat unitaire
f11 = ∫
800 (nX )1(nX )1 dx + 1400 (nX )1(nX )1 dx
∫ E1 D=30
0 EACB EABA
800 Z
d=20
800 1×1 1400 1×1
Y
f11 = ∫ dx + ∫ dx C
0 EA
CB
800 EA
BA
G x
B P1=1
800 A
1 800 1 1400
EACB ∫0 EABA ∫800
f11 = dx + dx 600
(nX)1
f11 =
800 600
+ =
800× 4
+
600× 4 + 1
EACB EABA 70.10 ×π × 30 70.103 ×π × 202
3 2
f11 = 4,345.10−5
f22 = ∫
800 (nX )2 (nX )2 dx + 1400 (0)2 (0)2 dx
Etat unitaire
E2 D=30 0 EACB 800 ∫
EABA
800 (−1) (−1)
Z
Y
P2=1 d=20 f22 = ∫ 2 2
dx
C 0 EACB
G x
B 800 800× 4
800 A f22 = =
600 EACB 70.103 ×π × 302
(nX)2
- 1
f22 = 1,616.10−5
f12 = ∫
800 (nX )1(nX )2 dx + 1400 (nX )1(0)2 dx
0 EACB 800 ∫
EABA
800 (+1) (−1)
f12 = ∫ 2 2
dx Matrice de flexibilité :
0 EACB
800× 4 4,345 − 1,616 
[ f ] =10−5 
800
f12 = − =−  mmN −1
EACB 70.103 ×π × 302  − 1,616 1,616 
f12 = −1,616.10−5
Nous en déduisons la matrice des déplacements généralisés :
[D] = [ f ][P ]  D1   f11

D  =  f
f12   P1 
f 22   P2 
 D1  −5  4,345 − 1,616   30  3
 D  =10 − 1,616 1,616  15010
 2   21  2   
D1 = 4,345 .10−5 × 30.103 − 1,616 .10 −5 × 150.103 = −1,12mm

D2 = −1,616 .10 −5 × 30.103 + 1,616 .10−5 × 150.103 = 1,94mm
D1 est négatif : le déplacement de la section A s’effectue en sens inverse de celui de
l’application de la charge P1. •A=-1,12 mm
D2 est positif : le déplacement de la section B s’effectue dans le même sens que celui
de l’application de la charge P2. •B=-1,94 mm
TD-2A-S4-F412
PROBLÈME N°9
Ce problème est le même que le problème 1 qui a été résolu partiellement par la
méthode de CLAPEYRON.
Les barres AB et BC (de même section droite A et de

C même matériau) sont articulées en A, B et C. Une
force verticale P1 est appliquée en B. On veut calculer
le déplacement vertical et le déplacement horizontal
du point B par la méthode de MAXWELL-MOHR.
Suggestion : Rajoutez en B une force horizontale

fictive P2 pour calculer le déplacement horizontal.
45°
B
P2 1°) Déterminer l’expression littérale de la matrice de
A
flexibilité de la structure (mettre L/EA en facteur).
P1 2°) En déduire l'expression littérale du déplacement
vertical et du déplacement horizontal du point B.
L
Application numérique : P1=10kN, L=2m, E=200GPa et A=1cm2. Calculez la valeur du

déplacement vertical et du déplacement horizontal du point B (en mm).
Réponses :
L ( 1 + 2 2 ) −1
[f ] =  
EA  −1 1 
2°) • = D = (1 + 2 2 )PL = 3.83 mm ↓ •HB = D2 = −
PL
= −1mm ←
VB 1
EA EA
RÈPONSES N°9
La méthode de Maxwell Mohr ou méthode de la charge

C
unité ne permet de calculer que les déplacements
des points d’applications des charges généralisées
Sur le problème initial qui nous est proposé, nous ne pouvons

donc déterminer que le déplacement vertical du point B.
Nous allons donc créer un problème plus général, dont l’initial

ne sera qu’un cas particulier.
45°
Nous désirons connaître le déplacement horizontal de la B
A
section B. Il est donc nécessaire qu’il existe en B
une charge horizontale.
P
Appliquons alors en B une charge « fictive » P2 : L
TD-2A-S4-F412
Etat initial Etat

C C intermédiaire
E0
E0’
E0 • E0’{ P1=P ; P2=0 }
45° 45°
A
B A
B P2
P P1
L L
Nous allons donc travailler sur l’état E0’. Décomposons ce dernier en deux états unitaires E1
et E2.
Etat unitaire
C E1
Etat initial
Etat intermédiaire
C E0 E0’
C
45°
B
A
L
P1=1
45° 45°
A
B B P2
A
L
P P1
L Etat unitaire
C E2
E0 • E0’ { P1=P ; P2=0 }
E0’ • P1×E1+P2×E2 45°

P2=1
A
B
L
TD-2A-S4-F412
1°) Expression littérale de la matrice de flexibilité de la structure :

Cherchons les efforts normaux s’exercant sur les deux barres, et pour les deux états
unitaires. Etat unitaire
C 2 E1
(n ) = −1 AB +
(nXBC)1
Etat 1 a
(n ) = + 2
x 1
BC
x 1 45°
A B
- P1=1
(nXAB)1 L 1
Etat unitaire
C E2
(n ) = +1
AB
Etat 2 a
(n ) = 0
x 2
BC
x 2
45°
P2=1
A
+ B
(nXAB)2 L 1

(nX )i (nX ) j (tY )i (tY ) j (tZ )i (tZ ) j (mX )i (mX ) j (mY )i (mY ) j (mZ )i (mZ ) j
∑∫ (
L
fij = + + + + )dx
barres
Dans notre cas :
(nX )i (nX ) j (n ) (n )
AB AB
(n ) (n )
BC BC
(n ) (n )
AB AB
(n ) (n )
BC BC
∑∫ (
L B C L L 2
fij = )dx = ∫ dx + ∫ dx = ∫ dx + ∫
X i X j X i X j X i X j X i X j
dx
barres
0 EA A EA B EA 0 EA 0 EA
(n ) (n )
AB AB
(n ) (n )
BC BC
(−1)1(−1)1 dx + ( 2) ( 2 ) dx = L (1+ 2 2)
f11 = ∫ dx + ∫ dx = ∫ ∫
L L 2 L L 2
X 1 X 1 X 1 X 1 1 1
0 EA 0 EA 0 EA 0 EA EA
(n ) (n )
AB AB
(n ) (n )
BC BC
(+1)2 (+ 1)2 dx = L
f 22 = ∫ dx + ∫ dx = ∫
L L 2 L
X 2 X 2 X 2 X 2
0 EA 0 EA 0 EA EA
(n ) (n )
AB AB
(n ) (n )
BC BC
(−1)1(+1)2 dx + ( 2 ) (0) dx = − L
f12 = ∫ dx + ∫ dx = ∫ ∫
L L 2 L L 2
X 1 X 2 X 1 X 2 1 2
0 EA 0 EA 0 EA 0 EA EA
L (1 + 2 2 ) − 1
[f ]=  
EA  − 1 1
TD-2A-S4-F412
2°) Expression littérale du déplacement vertical et du déplacement horizontal du

point B :
Dans l’état E0’ [D ]E = [ f ]E0 ' [P ]E0 '  D1  δ vB 

 D  =  B  = EA 
L (1 + 2 2 ) − 1  P1 
 
 2  E0 ' δ h  E0 '  −1
0'
1   P2  E0 '
δ vB  L (1 + 2 2 ) − 1  P PL (1 + 2 2 )
Dans l’état E0  B =    =  
δ h  E0 EA  − 1 1   0  E0 EA  − 1 
Application numérique : P1=10kN, L=2m, E=200GPa et A=1cm2.
PL 10.103 × 2.103
δ vB = (1 + 2 2 ) = (1 + 2 2 ) = 3,828mm •vB=3,83 mm
EA 200.103 ×100
PL 10.103 × 2.103
δ hB = − =− = −1mm •hB=-1 mm
EA 200.103 ×100
Le déplacement horizontal de B s’effectue en sens inverse de la charge « fictive » P2
PROBLÈME N°10
D=d 2 P En appliquant la méthode de MAXWELL-MOHR et en ne tenant

Z d compte que de l’énergie de déformation élastique due au moment
A G x fléchissant, déterminez en fonction de P, L, E et d l’expression
B C
littérale du déplacement vertical et de la rotation en C. En
déduire l’expression littérale de l’énergie de déformation.
Application numérique : P=1425N, E=210 GPa,
L L/2
L=1m, et d=40mm. Calculez Wdéf en Joules.
Réponses : 20L3 P 24L2 P 10L3 P 2

•VC = •C = Wdef = Wdef = 12J
•Ed 4 •Ed 4 •Ed 4
RÈPONSES N°10
Sur le problème qui nous est proposé, nous ne pouvons déterminer que le déplacement
vertical de la section C.
Nous allons donc créer un problème plus général, dont l’initial ne sera qu’un cas particulier.
Nous désirons connaître la rotation de la section C. Il est donc nécessaire qu’il existe en C
un moment « fictif » porté par l’axe de la rotation cherchée .
Appliquons alors en C un moment « fictif » P2 suivant l’axe Y:
TD-2A-S4-F412
Etat initial
Z E0
d•2
Y
A
G d
x
P
B
L
C
L/2
E0 • E0’{ P1=P ; P2=0 }

Etat
Z
intermédiaire
d•2
A E0’
G d
x
P1
B
L
C
P2
L/2
Constituons deux états unitaires de l’état E0’ et cherchons en les diagrammes des moments:
Z
d•2 Etat unitaire
Z Etat unitaire
Y E1 Y d•2 E2
A A
G x d G x d
P1=1
B B
3L/2 + C C P2=1
L (mY)1 L/2 (mY)2
L
+
L/2 L/2 1
(m ) y 1 =
3L
−X (m )
y 2 = +1
2 E0’ • P1×E1+P2×E2
TD-2A-S4-F412

(nX )i (nX ) j (tY )i (tY ) j (tZ )i (tZ ) j (mX )i (mX ) j (mY )i (mY ) j (mZ )i (mZ ) j
∑∫ (
L
fij = + + + + )dx
barres
Dans notre cas nous ne tenons compte que de l’énergie de déformation due au moment mY :
(mY )i (mY ) j
fij = ∫
C
dx
A EIY
B (mY )i (mY ) j C (mY )i (mY ) j
Compte tenu de la discontinuité de section en B : fij = ∫ A AB
EIY
dx + ∫
B EIYBC
dx
Calcul de f11 :
3L 3L 3L 3L
L
( − X )1 ( − X )1 3L ( − X )1 ( − X )1
f11 = ∫ 2 2 dx + ∫ 2 2 2 dx
0 EIYAB L EIYBC
Utilisons le tableau des Intégrales de Mohr pour calculer ces intégrales :
2 2
3L/2
L/2
L/2
En notation « symbolique » : f11 = 1/EIyAB + + 1/EIyBC
L
+ L/2
πd 4 πd 4
Or IYAB = et IYBC =
16 64
2 2
3L/2 L/2
L/2
f11 = 16/E•d4 + + 64/E •d4
L
+ L/2
16  1   3L  3L L  L    1  L 2 L  
2 2
f11 =    + × +  L + 4     
Eπd 4  3   2  2 2 2 

3 2  2 
  
20 L 3
f 11 =
Eπd 4
Calcul de f12 :
3L 3L
L
( − X )1 (+1)2 3L ( − X )1 (+1)2
f12 = ∫ 2 dx + ∫L2 2 dx
0 EIYAB EIYBC
3L/2
L/2 1 L/2 1
4
f12 = 16/E•d + x + + 64/E •d4 + L/2 +
L L
x
L/2
   3L  L  1 L L 
  (1) 
16 1
f12 =  × 1×    +  L + 4 
Eπd 4  2  2  2   2 2  2  
24 L2
f12 =
Eπd 4
TD-2A-S4-F412
Expression littérale du déplacement vertical C.
Dans l’état E0’ [D ]E = [ f ]E0 ' [P ]E0 '  D1 

D  =
δ vC 
 C = 
 f11 f12   P1 
0'
 2  E0 ' θY  E0 '  f 21 f 22   P2  E
0'
δ vC   f11 f12   P 
Dans l’état E0  C = 
θY  E0  f 21 f 22   0  E
0
20 PL 3
δ C
= f 11 P = •vC=20PL3/E•d4
E πd 4
v
Expression littérale de la rotation en C.

D  =
δ vC 
 C = 
 f11 f12   P1 
0'
 2  E0 ' θY  E0 '  f 21 f 22   P2  E
0'
δ vC   f11 f12   P 
Dans l’état E0  C = 
θY  E0  f 21 f 22   0  E
0
24 PL2
θYC = f 21P = •YC=24PL2/E•d4
Eπd 4
Expression littérale de l’énergie de déformation.
f f12   P1 
Wdéf =
1t
[P ][ f ][P ] = 1 [P1 P2 ]E0 '  11
f 22   P2  E
Dans l’état E0’
2 2  f 21 0'
f f12   P  1 10 P 2 L3
= [P 0] 11
1
= f11 P =
2
f 22   0  2
Dans l’état E0 Wdéf
2  f 21 Eπd 4
Application numérique : P=1425N, E=210 GPa, L=1m, et d=40mm.
10 P 2 L3 10 × 14252 ×10003
Wdéf = = Wdéf=12J
Eπd 4 210.10 3 × π × 40 4
TD-2A-S4-F412
PROBLÈME N°11
L/2 L/2 En appliquant la méthode de CASTIGLIANO et

en ne tenant compte que de l’énergie de
Z Z P déformation élastique due au moment
A x C B fléchissant :
G Y G a) Déterminez l'expression littérale de la
flèche maximum (en C) et des rotations aux
extrémités A et B en fonction de P, L, E et II.
b) Calculer numériquement la flèche (en mm) et les rotations (en°) pour P = 12kN, L=1m; la
section est un rectangle 80×40, E = 200 GPa.
Réponses : PL3 PL2 fmax = 0.73 mm •A = −•B = 2.210 −3 rd = 0,126 °

fmax = •A = −•B =
48EIY 16EIY
RÈPONSES N°11
Nous ne pouvons déterminer que le déplacement vertical de la section C.
Ce déplacement correspond à la flèche maximale, car le problème est symétrique en

l’abscence de charge horizontale.
L/2 L/2
P
Z Y
A G x C B
Nous désirons connaître la rotation de la section A. Il est donc nécessaire qu’il existe en A
Appliquons alors en A un moment « fictif » P2 suivant l’axe Y:
Nous désirons connaître la rotation de la section B. Il est donc nécessaire qu’il existe en B
TD-2A-S4-F412
Appliquons alors en B un moment « fictif » P3 suivant l’axe Y , mais compte tenu de la

symétrie :•A= -•B. L’approche de P3 n’est donc pas necessaire.
L/2 L/2
Etat initial
P E0
Z Y
A G x C B
L/2 L/2
P2 P1
Z Y Etat
intermédiaire
A G x C B E0’
E0 • E0’{ P1=P ; P2=0 }
Pour appliquer le théorème de Castigliano, nous devons au préalable calculer l’énergie de

déformation élastique dans l’état E0’ .
∑∫
L
Wdéf = ( + + + + + )dx
barres
0 2EA 2GA 2GA 2GIG 2EIY 2EIZ
(M ) 2
=∫
B Y E0 '
En ne tenant compte que de l’énergie due au moment fléchissant : Wdéf dx
A 2EIY
Nous devons donc chercher le moment My. Pour cela, décomposons l’état E0’ en deux états :
TD-2A-S4-F412
Etat Etat
L/2 L/2
intermédiaire E1
Z P1
E0’ Y
L/2 L/2 A G C B
x
P2 Z P1
Y
A G x
C B (MY)E1 -
P1L/4
L/2 L/2
P2 - P2 Z Y
(MY)E0’ P2/2+P1L/4 A G x
C B
E0’ • E1+ E2
Etat
P2
-
E2
(MY)E0’=(MY)E1+(MY)E2 (MY)E2
Wdéf = ∫
B (M ) 2
dx = ∫
B ((M ) + (M ) ) dx = (M ) + (M )
2
B
2 2
+ 2 MY ( ) (M )
∫
Y E0 ' Y E1 Y E2 Y E1 Y E2 E1 Y E2
dx
A 2 EIY A 2 EIY A 2EIY
(M ) 2
∫ (M ) ∫ (M ) ∫ (M ) (M )
B 1 B 1 B 1 B
Wdéf = ∫
2 2
dx = dx + dx +
Y E0 '
Y Y Y Y dx
A 2 EIY 2EIY A E1
2 EIY A E2
EIY A E1 E2
Pour calculer Wdéf nous pouvons soit faire le calcul avec l’état E0’, soit avec les états E1 et E2.
Première approche :
Utilisons le tableau des intégrales de Mohr avec l’état E0’.

B (M )
2
Wdéf = ∫
Y E0 '
dx
A 2EIY
2 2
P2/2+P1L/4 P2/2+P1L/4
P2
Wdéf = 1/2EIY - + 1/2E IY -
L/2 L/2
 1   P2 P1L   P2 P1 L   L  1  P2 P1L  L 
2 2

Wdéf =
1
 (P2 ) + P2 ×  +  +  +   +   +  
2

2 EIY  3  2 4  2 4   2  3 2 4  2 
P12 L3 P 2 L P P L2
Wdéf = + 2 + 1 2
96 EIY 6 EIY 16 EIY
TD-2A-S4-F412
Autre approche : à ‘aide des états E1 et E2
∫ (M ) ∫ (M ) ∫ (M ) (M )
1 B 2 1 B 2 1 B
Wdéf = Y dx + Y dx + Y Y dx
2EIY A E1
2EIY A E2
EIY A E1 E2
P1L/4 2 P2 2 P1L/4 P2
- - -
Wdéf = 1/2EIY + 1/2E IY + 1/E IY x -L
L L L
1  1  P1L  
2
1   1  P2  2   1  1  P1 L P2  
Wdéf =    L +     L  +   × L
2 EIY  3  4   2 EIY   3  2   EIY 4  4 2 
 
P12 L3 P22 L P1P2 L2

Wdéf = + +
96 EIY 6 EIY 16 EIY
a) Expression littérale de la flèche maximum (en C)
∂Wdéf
Appliquons le théorème de Castigliano : Di =
∂Pi
La dérivée partielle de l’énergie par rapport à P1 nous donne le déplacement associé à P1.
( E0 ' ) C ( E0 ' )
∂Wdéf P1 L3 P1 L2
D1 =δ = = +
∂P1
v
48EIY 16EIY
Revenons maintenant au problème initial, sachant que : E0 • E0’{ P1=P ; P2=0 }
C ( E0 )
P L3
δ v =
48EIY
a) Expression littérale des rotations aux extrémités A et B
La dérivée partielle de l’énergie par rapport à P2 nous donne le déplacement associé à P2.
( E0 ' ) A ( E0 ' )
∂WdéfP2 L P1L2
D2 =θ = = +
∂P2
y
3EIY 16EIY
TD-2A-S4-F412
Revenons maintenant au problème initial, sachant que : E0 • E0’{ P1=P ; P2=0 }
A ( E0 )
P L2 ( E0 )
θ y = = −θ yB
16EIY
b) Numériquement : P = 12kN, L=1m; la section est un rectangle 80×40, E = 200 GPa.
C ( E0 )
P L3 12.103 × 10003 ×12
δ = = = 0,73mm
48EIY 48 × 200.103 × 40 × 803
v
•vC=fmax=0,73 mm
A ( E0 ) 12.103 ×10002 ×12

P L2
θ = = = 2,197.10−3 rd ⇒ 0,126o
16EIY 16 × 200.10 × 40 × 80
y 3 3
•YA= -•YB =0,126°
PROBLÈME N°12
L L Considérons le système hyperstatique ci-

contre. On ne tient compte que du moment
Z Z P fléchissant pour le résoudre.
Y A x C B 1°) Déterminez les expressions littérales des
G G réactions. Tracer le diagramme des moments
fléchissants.
2°) Déterminez les expressions littérales de
la flèche du point C et des rotations en C et B.
Application numérique : P=12kN, L=1m, E=200GPa, la section est un rectangle de
80mmx37.5mm. Calculez la valeur du déplacement vertical de C (en mm) et les rotations de B
et C (en rd).
Réponses :
Les réponses sont données dans le repère XYZ.
1°) AZ = 11 P = 8.25kN 5P 3 PL
= 3.75kN MAY =
BZ = = 4500 Nm
16 16 8
3 PL
Moments fléchissants MY MYA = = 4500 Nm MYC = − 5 PL = -3750Nm MYB = 0
8 16
2
3 PL PL2
2°) δVC = − 7 PL = -2.7344mm θ YC = = 1.171910 −3 rd θYB = − = −4.687510 −3
rd
96EIY 32 EI Y 8 EIY
TD-2A-S4-F412
RÈPONSES N°12
L=1m
Z P=12kN
VA L=1m
MA Y
HA A b =37,5mm
G x VB
C
h =80mm
B
Problème plan (La géométrie est plane ainsi que le chargement) : H=i-3n .
3 degrés de liberté de supprimés en A (Encastrement) . Inconnues :HA ;VA ;MA
1 degré de liberté de supprimé en B (Appui glissant) . Inconnue : VB
Un solide pour appliquer le premier principe de Newton.
H= i-3n = 4-3 = 1 Le système est hyperstatique externe de degré 1 : Hext1°
Le « Mathématicien » ne sait pas résoudre un système ou il existe plus d’inconnues que

d’équations, le « Physicien » ne faisant pas mieux , il faut donc chercher des équations de
déformations, qui rajoutées à celle de statique permettront de trouver les 4 inconnues.
Notre problème est hyperstatique de degré 1. Nous allons donc choisir une inconnue
hyperstatique parmi les 4 (HA ;VA ; VB MA ) afin de traiter un problème isostatique, dont
l’hyperstatique ne constituera qu’un cas particulier. Plusieurs choix sont possibles pour
l’inconnue hyperstatique, donc plusieurs résolutions possibles. Le chargement de l’isostatique
sera celui de l’hyperstatique plus l’inconnue choisie :
TD-2A-S4-F412
Choix de l’inconnue hyperstatique : MA
Remplaçons l’encastrement en A par un appui simple.
Z
L
L Z
L
P P2 L
Y P1
Y
A G x A
C G x
C
B
B
Etat hyperstatique Etat isostatique

E0 E0 • E0’{ P1=P ; D2=0 }
E0’
Les 2 états seront équivalents lorsque la rotation de la section A de l’isostatique sera nulle.
Choix de l’inconnue hyperstatique : VB
Supprimons la liaison en B
Z
L L
Z
L L
P P1
Y
Y P2
A G x A
C G x
C
B B

E0 • E0’{ P1=P ; D2=0 }
E0 E0’
Les 2 états seront équivalents lorsque le déplacement vertical de la section B de létat

isostatique sera nul.
Il nous faut donc calculer des déplacements dans les états isostatiques.
Pour la suite du corrigé , conservons ce dernier choix de l’inconnue hyperstatique (VB), et

décomposons l’état isostatique E0’ en 2 états unitaires, afin de calculer D2, à l’aide de
Maxwell-Mohr.
TD-2A-S4-F412
Etat unitaire Z
L P1=1
E1 L
Y
A
G x
C
Z
L
L
+ B
P1 L (mY)1
Y P2
A
G x
C
B
Z
L
L P2=1
Etat isostatique Y
E0’ A
G x
C
B
E0’ • P1×E1+P2×E2 2L
+
Etat unitaire
(mY)2 L
E2
Dans l’état E0’ [D ]E = [ f ]E0 ' [P ]E0 '  D1  δ vC 

D  =  B  =  f
 f11 f12   P1 
0'
 2  E0 ' δ v  E0 '  21 f 22   P2  E
0'
Or dans l’état hyperstatique, le déplacement vertical de B est nul, P2 s’identifie alors avec la
valeur de la réaction VB de l’état E0.
f 21P1
D2 = δ vB = f 21P1 + f 22 P2 = 0 ⇒ P2 = −
f 22
(mY )i (mY ) j
fij = ∫
C
Calcul des coefficients d’influence dx à l’aide des intégrales de Mohr :
A EIY
Calcul de f11 :
2
L
f11 = ∫
C ((mY )1 )2 dx = 1 1 2
× ×L ×L f11 =
L3
f11 = 1/EIy +L A EIY EIY 3 3EIY
Calcul de f22 :
2
2L
f22 = ∫
B ((mY )2 )2 dx = 1 1
× × 4L2 × 2L f22 =
8L3
f22 = 1/EIy + 2L A EIY EIY 3 3EIY
TD-2A-S4-F412
(mY )1 (mY )2 B (m ) (m )
f12 = ∫ dx + ∫ Y 1 Y 2 dx
C
Calcul de f12 :
A EIY C EIY
2L
L L
L
f12 = 1/EIY + x +L + 1/EIY 0
x +
L L
L
1  5L3
 × L × (2(2 L ) + L )L 
1
f12 = f12 =
EIY 6  6 EIY
Calcul de la réaction VB :
5 L3
P1 5
P2 = −
f 21 P1
=−
6 EIY P2 = − P1
f 22 8 L3 16
3EIY Z
L P
L 5P/16
Compte tenu du choix arbitraire du sens
Y
de P2 que nous avons fait, et trouvant P2
négatif, nous en concluons que la réaction A G x
C
verticale en B dans l’état hyperstatique est B
verticale ascendante.
Sa norme étant de 5P/16.
Le problème initial est maintenant isostatique. Pour trouver les 3 inconnues externes
restantes HA ;VA ;MA , nous pouvons utiliser les 3 équations de statique, ou calculer sur
l’état E0’ en faisant P2=5P/16
E0 • E0’{ P1=P ; D2=0 } • P1×E1+P2×E2 • P×E1-(5P/16)×E2
HA0 = P×HA1-(5P/16)×HA2 HA0 = P×0-(5P/16)×0 = 0

VA0 = P×VA1-(5P/16)×VA2 VA0 = P×1-(5P/16)×1 = 11P/16
MA0 = P×MA1-(5P/16)×MA2 MA0 = P×L-(5P/16)×2L = 3PL/8
Diagramme des moments fléchissants :
Pour trouver le diagramme des moments fléchissants, nous utilisons les résultats obtenus sur
l’état isostatique E0’ au lieu de tout refaire le problème initial.
MY (E0) = P×mY (E1)-(5P/16)×mY (E2)
MYA(E0) = P×mYA(E1)-(5P/16)×mYA(E2) = P×L-(5P/16)×2L = 3PL/8

MYC(E0) = P×mYC(E1)-(5P/16)×mYC(E2) = P×0-(5P/16)×L = -5PL/16
MYB(E0) = P×mYB(E1)-(5P/16)×mYB(E2) = P×0-(5P/16)×0 = 0
TD-2A-S4-F412
1 Etat unitaire
E1
Z
L P1=1
L
Etat initial
P x E1-5P/16 x E2 Y
11P/16 A
G x
L C
Z
L
B
P
L 5P/16 +
Y L (mY)1
A G x
3PL/8 C Etat unitaire
B 1 E2
5PL/16 L
+ - Z
L P2=1
3PL/8 (MY)0 Y
A
G x
2L C
B
2L
+
(mY)2 L
Etat0 • P1×Etat1+P2×Etat2 • P×Etat1-(5P/16)×Etat2
Expressions littérales de la flèche du point C et des rotations en C et B
flèche du point C :
 5
δ vC   D1  L3 1 2   P1 
Dans l’état E0’ :  B =   = 5  
δ v  E0 '  0  E0 ' 3EIY  8   P2  E0 '
2 
 5
δ C  1  P 
Dans l’état E0 :  vB  =
L3 2   5P  δ vC
( E0 )
=
L3  5 5P 
P− ×  •vC=7PL3/96EIY
5  −
δ v  E0 3EI Y  8   16  3EIY  2 16 
2 
Le déplacement s’effectue dans le même sens que celui de P1.
Rotations en C et B :
Pour rechercher les rotations des sections C et B il faut qu’il existe en C et B des moments
portés par l’axe des rotations cherchées.
Compte tenu des résultats dont nous disposons, il sera plus simple de raisonner à partir de
l’état isostatique équivalent E0’ . Cherchons donc, dans un premier temps, les rotations des
sections C et B dans l’état isostatique E0’. Rajoutons sur ce dernier 2 couples P3, et P4 en C et
D.
TD-2A-S4-F412
Etat unitaire
Z
L P1=1 E1
L
Y
A
G x
C
(mY)1
+ B
L
Etat unitaire
Etat isostatique Z
L E2
E0’ L P2=1
Y
A
G x
L C
Z
L B
Y
P1
P2 2L
+
A (mY)2 L
G x
C Etat unitaire
Z
L
B L E3
P3
P4 Y
A
G x
C P3=1
B
+ 1
E0’ • P1×E1+P2×E2+ P3×E3+P4×E4
(mY)3
Z
L
L Etat unitaire
E4
Y
A
G x
C
E0 • E0’{ P1=P ; P2=-5P/16 ; P3=0 ; P4=0 } B
+ P4=1
(mY)4
1
Décomposons l’état E0’ en 4 états unitaires pour appliquer la méthode de Maxwell-Mohr
 L3 5 L3 
δ C
  D1   3EI 6 EIY
f13 f14  P
 1
 3Y  
v
   
δ
B
=  f 24   2 
v  =0 5L 8L3 P
Dans l’état E0’ : f 23
θ C   D3  6 EIY 3EIY P 

y
     3 
θ  f 31
B f 32 f33 f 34  P
y  E0 '  D4  E0 '  4  E0 '
 f 41 f 42 f 43 f 44 
TD-2A-S4-F412
 L3 5 L3 
δ C
  3EI 6 EIY
f13 f14   P 
 3Y   − 5P 
v
 
δ
B
=  f 24   16 
v  5L 8L3
Dans l’état E0 : f 23
θ C  6 EIY 3EIY  0 
   
f34  0 
y
θ  f 31
B f 32 f 33
y  E0   E0 '
 f 41 f 42 f 43 f 44 
Nous devons donc calculer seulement les coefficients f31, f32, f41, et f42 .
(mY )1 (mY )3 B (m ) (m )
f13 = ∫ dx + ∫ Y 1 Y 3 dx
C
Calcul de f13 :
A EIY C EIY
L 1
f13 = 1/EIY +L + + 1/EIY
x L
0
x 0
L L
1 1  L2
f13 =  × 1× L × L  f13 =
EIY 2  2 EIY
(mY )2 (mY )3 B (m ) (m )
f23 = ∫ dx + ∫ Y 2 Y 3 dx
C
Calcul de f23 :
A EIY C EIY
2L
L 1 L
0
F23 = 1/EIY + + + 1/EIY +
L
x L
x
L L
1 1  3L2
f 23 =  ×1× (2 L + L )L  f 23 =
EIY  2  2 EIY
(mY )1 (mY )4 B (m ) (m )
f14 = ∫ dx + ∫ Y 1 Y 4 dx
C
Calcul de f14 :
A EIY C EIY
L 1 1
f14 = 1/EIY +L + + 1/EIY +
x L
0
x
L L
1 1  L2
f14 =  × L × 1× L f14 =
EIY  2  2 EIY
TD-2A-S4-F412
(mY )2 (mY )4 B (m ) (m )
f24 = ∫ dx + ∫ Y 2 Y 4 dx
C
Calcul de f24 :
A EIY C EIY
2L
L 1 L 1
F24 = 1/EIY + + + 1/EIY + +
L
x L
x L
L
1  2 L2
 × 1× (2 L + L )L + × 1× L × L 
1 1
f 24 = f 24 =
EIY 2 2  EIY
 L3 5 L3 L2 L2 
 3EI 6 EIY 2 EIY 2 EIY 
 3Y
δ vC   5L 8 L3 3L2 2 L2   P 
 B  6 EI  − 5P 
δ v  EIY  
 16 
3EIY 2 EIY
Dans l’état E0 : =  2Y
θ yC   L 3L2  
 B f 33 f 34   0 
 2 EIY  0 
θ y  E0
2 EIY
 L2 2 L2   E0 '
 f 43 f 44 
 2 EIY EIY 
PL2 3L2 5 P PL2

θ yC = − = •yC=PL2/32EIY
2 EIY 2 EIY 16 32 EIY
PL2 2 L2 5 P PL2
θ yB = − =− •yB=-PL2/8EIY
2 EIY EIY 16 8EIY
Application numérique : P=12kN, L=1m, E=200GPa, la section est un rectangle de

80mmx37.5mm.
Dans le repère XYZ :
11P 5P 3PL
AZ = = 8.25kN BZ = = 3.75kN M AY = = 4500 Nm
16 16 8
3PL 5 PL
Moments fléchissants MY M YA = = 4500 Nm M YC = − = -3750Nm M YB = 0
8 16
7 PL3 PL2 PL2

δ VC = − = -2.7344mm θ YC = = 1.171910 −3 rd θYB = − = −4.687510 −3 rd
96 EI Y 32 EIY 8 EIY
TD-2A-S4-F412
L=1m P=12kN
Z
L=1m
4500 Nm Y
A
2,7344 mm
G x
b =37,5mm
C B
8,25 kN h =80mm
4,6875.10-3 rd
-3
1,1719.10 rd
3,75 kN
3750 Nm
A
-
+ C
B
(MY)
4500 Nm
PROBLÈME N°13
La poutre encastrée AB ci-contre, de section carrée 60×60, est haubanée par un

câble AC de diamètre 5 mm. L = 1 m, E = 200 GPa, P = 3 kN. Le système est
hyperstatique. Calculez :
1°) la tension dans le câble.
P 2°) le déplacement vertical du point A.
45° Z
x
A G B Réponses : N xAC = 2893 N δVA = 1 ,47 mm
TD-2A-S4-F412
RÈPONSES N°13
Deux solides (La barre AB et le cable AC),

C
3 inconnues en B (encastrement),
2 inconnues en A et en C (articulations) :
H= i-3n = 7-3x2 = 1
P
Le système est hyperstatique 45° Z

Y
externe de degré 1 : Hext1°
G x
Le cable AC étant articulé à ses deux extrémités A L B
est donc soumis qu’à un effort normal.
Choisissons donc comme inconnue hyperstatique

la réaction en C.

P2
E0 E0’
C C
E0 • E0’{ P1=P ; D2=0 }
P P1
45° Z 45° Z
Y Y
G x x
G
A B A B
L L
Décomposons l’état isostatique E0’ en deux états unitaires E1 et E2.
TD-2A-S4-F412
Etat unitaire Etat unitaire P2=1

E1 E2
C C
(nX)2
+
E0’ • P1×E1+P2×E2
P1=1
45° Z 45° Z
Y 1 Y
G x G x
A B A B
L L
+ (mY)2
-
L
(mY)1 L 2/2

D  =
 δ vA 
 C  =
 f11 f12   P1 
0'
 2  E0 ' δ AC  E0 '  f 21 f 22   P2  E
0'
Or dans l’état hyperstatique, le déplacement de C dans la direction de P2 est nul, P2

s’identifie alors avec la valeur de la réaction en C de l’état E0 .
f 21 P1
D2 = δ AC
C
= f 21 P1 + f 22 P2 = 0 ⇒ P2 = −
f 22
(nX )i (nX ) j (mY )i (mY ) j
∑∫ (
L
Calcul des coefficients d’influence fij = + )dx à l’aide des
barres
0 EA EIY
intégrales de Mohr :
Calcul de f11 :
2
L
f11 = ∫
B ((mY )1 )2 dx = 1 1 2
× ×L ×L f11 =
L3
f11 = 1/EIyAB +L A EIYAB EIYAB 3 3EI YAB
f22 = ∫
((mY )2 )2 dx + C((nX )2 )2 dx
∫
B
Calcul de f22 : A EIYAB A EAAC 2
2
L 2/2 1
f22 = 1/EIyAB - +1/EAAC +
L L 2
 1 1  2 
2
  1  L3 L 2
f 22 =  AB × ×  L × L +  AC ×12 × 2L f22 = AB
+ AC
 EIY 3  2    EA  6EIY EA
TD-2A-S4-F412
(mY )1 (mY )2
f12 = ∫
B
Calcul de f12 : dx
A EIYAB
L L 2/2
f12 = 1/EIy AB
+ x -L
L
1 1  L 2   2 L3
f12 = −  × L ×   L
  f12 = −
EI YAB  3  2   6 EI YAB
Calcul de la réaction en C :
2 L3
− P1
f 21 P1 6 EI YAB 2
P2 = − =− 3
P2 = P1
f 22 L L 2 6 2 I YAB
AB
+ 1+
6 EI Y EA AC L2 A AC
1°) Tension dans le câble.
Nous avons choisi arbitrairement l’inconnue de l’action en C de telle sorte que le cable AC
soit tendue. Nous trouvons P2 positif : Le cable est donc en traction.
2
L’effort normal dans le cable est : N XAC = P
6 2 I YAB
1+
L2 A AC
Application numérique : Section AB carrée 60×60, haubanée par un câble AC de diamètre 5 mm. L = 1 m,
E = 200 GPa, P = 3 kN.
2
N XAC = 3.103 = 2893 N NXAC=2893 N
6 2 × 60 4 × 4
1+
12 ×1000 2 × π × 52
2°) Déplacement vertical du point A.

D  =
 δ vA 
 C  =
 f11 f12   P1 
0'
 2  E0 ' δ AC  E0 '  f 21 f 22   P2  E
0'
δ   f11
A
f12  3000
 =
v
Dans l’état E0
 0   f 21 f 22   2893
L3 2 L3 L3  2 
δ vA = f11 3000 + f12 2893 = AB
3000 − AB
2893 = AB 
3000 − 2893 
3EI Y 6 EI Y 3EI Y  2 
1000 ×12 
3
2 
δ vA = 4 
3000 − 2893  = 1,473mm •vA=1,47 mm
3 × 200.10 × 60 
3
2 
TD-2A-S4-F412
PROBLÈME N°14
L
La poutre encastrée AB ci-contre a son extrémité
poids P libre qui repose en A sur un appui élastique
Z
x modélisé par un ressort dont la raideur est k. Le
A système set hyperstatique.
G B
En l'absence de la charge répartie, la poutre est
rectiligne et le ressort n'est pas sollicité.
On ne tient compte que du moment fléchissant pour résoudre le système.
1°) Déterminez l'expression littérale de la force qui comprime le ressort et en déduire

l'expression de la flèche du ressort.
2°) Calculez numériquement la force qui comprime le ressort et la flèche du ressort pour k =
0.1, 0.5, 1, 10 N/mm et k → ∞ La barre est à section carrée 10×10, ρ = 7.85 10-6 kg/mm3
(prendre g=10m/s2), L = 1 m, E = 210 GPa.
k (N/mm) FA (N) f (mm)
Réponses : 0.1 0.47 4.71
3P 1 FA 0.5 1.44 2.88

FA = ( ) f =
8 3 EIY k
1 + 1 1.93 1.93
kL3
10 2.80 0.28
∞ 2.94 0
RÈPONSES N°14
q Poids P
Z
Y
10
A G x B
10
Problème plan (La géométrie est plane ainsi que le chargement) : H=i-3n .
3 degrés de liberté de supprimés en A (Encastrement) . Inconnues :HB ;VB ;MB
1 degré de liberté imposé en B (Appui élastique) . Inconnue : VA
Un solide pour appliquer le premier principe de Newton.
H= i-3n = 4-3 = 1 Le système est hyperstatique externe de degré 1 : Hext1°
TD-2A-S4-F412
Choix de l’inconnue hyperstatique : VA
Supprimons la liaison en A
E0
E0 • E0’{ P1=P ; D2=•x } E0’
P2
L L
q Z
Poids P=qL Z P1=q
Y VB Y
10 VA 10
B A B
10 A G x HB 10
G x
MB
Les 2 états seront équivalents lorsque le déplacement vertical de la section A de létat

isostatique sera égal à la flèche •x du ressort.
Décomposons l’état isostatique E0’ en 2 états unitaires, afin de calculer D2=•x (flèche du
ressort), à l’aide de Maxwell-Mohr.
Etat unitaire Etat unitaire

E1 E2
E0’ • P1×E1+P2×E2
L P2=1
L
Z P1=1 Z
Y Y
10 10
A G x
B A B
10 G x
10
+ -
(mY)1 (mY)2
L2/2
L
Dans l’état E0’ D2 = δ vA = f 21 P1 + f 22 P2
L (mY )i (mY ) j
Calcul des coefficients d’influence fij = ∫(
0 EIY
)dx à l’aide des intégrales de Mohr :
TD-2A-S4-F412
f22 = ∫
B ((mY )2 )2 dx
Calcul de f22 :
A EIY
2
L
f22 = 1/EIy -
3 L
 1 1 2  L
f 22 =  × × L × L f22 =
 EIY 3  3EIY
(mY )1 (mY )2
f12 = ∫
B
Calcul de f12 : dx
A EIY L2/2
L
f12 = 1/EIy + x -
L L
1 1  L2   L4
f12 = −  × L ×   × L  f12 = −
EIY 4  2  8EIY
L4 P1 L3 P2
Dans l’état E0’ D2 = δ vA = f 21 P1 + f 22 P2 D2 = δ vA = − +
8 EIY 3EIY
L4 q L3 FA
Dans l’état E0 δ vA = − + or FA = − kf et P = qL
8 EIY 3EIY
L4 q L3 FA F 3P 1
− + =− A d’où : FA = ( )
8EIY 3EIY k 8 1 + 3EIY
kL3
3P 1
FA = ( )
8 1 + 3EIY
kL3
Application numérique : La barre est à section carrée 10×10, ρ = 7.85 10-6 kg/mm3 (prendre
g=10m/s2), L = 1 m, E = 210 GPa
3 × 10 × 10 × 1000 × 7,85 .10 −6 × 10 1

FA = ( )
8 3 × 210 .10 3 × 10 4
1+
12 × 1000 3 × k
TD-2A-S4-F412
 
 1 
FA = 2,94375 
 1 + 0,525 
 k 
Pour k=0 :
:
L4 q L3 P 1000 3 × 10 × 10 × 1000 × 7,85.10 −6 × 10 × 12
δ =−
A
=− =− = 5,60714mm
8 × 210.103 × 10 4
v
8 EIY 8 EIY
Pour k= ∞ :
3P
La force FA tend vers FA = = 2,94N , et la flèche f tend vers 0 :
8
Variation de la force et de la flèche du ressort en fonction de la

raideur k
6
5,6 mm
k FA (N) f (mm)
(N/mm) FA
5 FA (N) 0.1 0.47 4.71
f (mm) 0.5 1.44 2.88
1 1.93 1.93 f
Force(N) flèche(mm)
4
10 2.80 0.28
∞ 2.94 0
2,94375 N
3
0
k(N/mm)
0 5 10 15 20
TD-2A-S4-F412
PROBLÈME N°15
5m 5m
Avant d'appliquer la charge répartie, un espace de 12 mm existe entre la
20kN/m
Z poutre et l'appui C. Sachant que l'inertie de flexion vaut 217 106 mm4,
calculez (en kN) les réactions sur chaque appui (prendre E = 200 GPa). On
A x B ne tient compte que du moment fléchissant.
12mm C G
Réponses : A = B = 50 kN C = 100 kN
RÈPONSES N°15
5m 5m
Etat initial
q=20kN/m
E0
Z Y
G x
A C B
12mm
Le problème deviendra hyperstatique externe de degré 1 lorsque la flèche verticale de la

section C sera égale à 12 mm.
Calculons donc la part de la charge répartie q qui donne une flèche égale à 12mm.
Nous devons donc constituer un état E0’ pour lequel nous rajoutons une force verticale en C
Etat isostatique
E0’ E0 • E0’{ P1=q ; P2=0 }
5m P2 5m
P1=q
Z Y
G x
A C B
TD-2A-S4-F412
Décomposons en deux états unitaires :
Etat unitaire Etat unitaire

E1 E0’ • P1×E1+P2×E2 P2=1 E2
L/2 L/2 L/2 1/2 L/2 L/2 1/2

L/2
P1=1 Z
q
Z Y Y
A B A C B
G x C G x
- -
(mY)2
(mY)1 L/4
L2/8
Dans l’état E0’ D2 = δ vC = f 21 P1 + f 22 P2
L (mY )i (mY ) j
Calcul des coefficients d’influence fij = ∫(
0 EIY
)dx à l’aide des intégrales de Mohr :
f22 = ∫
B ((mY )2 )2 dx
Calcul de f22 :
A EIY
2
L
f22 = 1/EIy -
L/4
 1 1  L  2
L 3
f 22 =  × ×   × L f22 =
 EIY 3  4   48EIY
(mY )1 (mY )2
f12 = ∫
B
Calcul de f12 : dx
A EIY
L L
f12 = 1/EIy - x -
L2/8 L/4
1  5 L  L2   5L4
f12 = +  × ×   × L  f12 = +
EIY 12 4  8   384 EIY
5 L4 P1 L3 P2
Dans l’état E0’ D2 = δ vA = f 21 P1 + f 22 P2 D2 = δ = v
A
+
384 EIY 48 EIY
5L4 q 384 EIY δ vA 384 × 200.103 × 217.10 6 × 12

Dans l’état E0 δ vA = ⇒ q= = = 4 Nmm −1
384 EIY 5L4
5 × 10 4 4
( )
TD-2A-S4-F412
Un chargement réparti de 4N par mm, c'est-à-dire 4kN/m procure une flèche de 12mm/
Le chargement total est donc pour 10 m de 40kN. Chaque réaction vaut donc 20kN.
5m 5m
20kN 20kN
4kN/m
Z Y
A G x C B
12mm
Il subsiste donc une charge répartie de 16 kN supportée par les 3 appuis A, C, et B :
5m 5m
Etat hyperstatique
q’=16kN/m
Z Y
G x
A C B
HA
VA VB
VC
Le système est hyperstatique externe de degré 1. Choisissons comme inconnue

hyperstatique, la réaction verticale en C et constituons un système isostatique équivalent :
5m P2 5m
Etat isostatique
Equivalent si D2=0 P1=q’
Z Y
G x
A C B
Compte tenu des résultats obtenus précédemment :

5 L4 P1 L3 P2
D2 = δ = f 21 P1 + f 22 P2
v
A
D2 = δ = v
A
+ =0
384 EIY 48 EIY
5 L4 P1 48 EIY 5P L
P2 = − 3
=− 1
384 EIY L 8
TD-2A-S4-F412
La réaction verticale en C dans l’hyperstatique vaut donc :

5q ' L 5 × 16 × 10
VC = − =− = −100kN
8 8
Les deux autres réaction verticales en A et B sont égales (symétrie), et valent en écrivant
l’équilibre : VA=VB=-
5m
100kN 5m
30kN q’=16kN/m 30kN

Z Y
G x
A C B
Nous pouvons maintenant revenir au problème initial en superposant les résultats obtenus
respectivement sur l’état isostatique et sur l’hyperstatique.
Etat initial
E0
100kN (100kN +0kN)
50kN (30kN +20kN) 50kN (30kN +20kN)
q=20kN/m (16kN/m+4kN/m)
A B
12mm C
5m 5m
TD-2A-S4-F412
PROBLÈME N°16
Y
A x B Considérons le cadre fermé ci-contre. Le
G système est hyperstatique. On ne tient compte
que du moment fléchissant pour résoudre le
système.
P C P
Déterminez le moment de flexion maximum (en
Nmm), la contrainte de flexion maximum (en
L/2= 300 mm
MPa) et l'allongement du cadre (en mm).
La section droite du cadre est rectangulaire
8×25. Le matériau est un acier dont le module
E D
de YOUNG vaut 210 GPa. La force P vaut 50 N.
L= 600 mm
Réponses:
3 PL
MZC = = 5625 Nmm σ xx
C
= ±21.9 MPa
16
5 PL3
δ = = 1.26mm
192 EIII
RÈPONSES N°16
Le cadre est statiquement déterminé, et en équilibre sous l’action des deux charges P.
Par contre il est impossible de déterminer les actions internesdu Torseur de section (ou de
Y
Cohésion) compte tenu qu’il est fermé.
G x A B
En effet nous ne pouvons pas
mettre en équilibre un tronçon
Z
isolé, puisque nous ne pouvons
pas définir un Amont ou un Aval L/2=300mm
Etat
par rapport à une coupure hyperstatique
(Le cadre est fermé). E0
Y
P P
MY F C
MZ
TY Hint3°
TZ
Mx
G x L=600mm
Nx
Z
Le problème est plan : D
TZ=MX=MY=0 E
Les inconnues sont :
NX ; TY ; MZ.
Le système est Hyperstatique Interne de degré 3, car les inconnues que nous ne pouvons
calculer sont les sollicitations du Torseur de Section.
TD-2A-S4-F412
Rendons le système isostatique en y pratiquant un nombre de coupures simples, égal à son

degré d’hyperstaticité.
Le cadre étant Hint3° faisons une coupure totale TY Y
MZ
dans une section quelconque H
A G x B
(triple : NX ;TY ;MZ) H’’ H’ Nx
H devient H’ et H’’
Appliquons des sollicitations NX ;TY ;MZ Z
aux deux lèvres de la coupure, égales
et opposées et d’intensité telles que Etat isostatique
nous refermions exactement la coupure. Equivalent
δ hH '/ H '' = 0
P  H '/ H '' P
δ hH '/ H '' = 0 F δ v =0 C
θ H ' / H ''
=0
δ vH '/ H '' = 0  Z
θ ZH '/ H '' = 0 Hint3°

Ce sont des déplaçements relatifs
D
E
Nous pouvons simplifier la résolution en se souvenant qu’un axe de symétrie nous permet
d’abaisser le degré d’hyperstaticité du système.
Faisons intervenir l’axe de symétrie vertical AE en considérant la coupure en A :
TY=0 TY=0 Y
A devient A’ et A’’ MZ
Nx
G xB
Appliquons des sollicitations NX ;TY ;MZ A’’ A’
TY=0
aux deux lèvres de la coupure, égales
Z
et opposées.
L’effort tranchant TY dans la section A’’
est vertical descendant pour respecter Etat isostatique
le principe d’action mutuelle. Equivalent

Il doit être vertical ascendant pour P δ hA'/ A'' = 0
 A'/ A'' P
respecter la condition de symétrie. F θ Z =0 C
La seule valeur possible est donc TY=0 Hint2°

Pour refermer exactement la coupure :
δ hH '/ H '' = 0
θ ZH '/ H '' = 0 D
E
TD-2A-S4-F412
Faisons intervenir les 2 axes de symétrie (vertical et horizontal) en considérant la

coupure en A et en E:
Y
A devient A’ et A’’ MZ
Nx G xB
E devient E’ et E’’ A’’ A’
L’équilibre de l’une ou l’autre des Z

2 parties, nous permet de trouver
la valeur de l’effort normal dans Etat isostatique
les 2 sectionss : NX=P/2 Equivalent
θ ZA'/ A'' = 0 P
P
F C
Hint1°
Pour refermer exactement
la coupure :
θ ZA'/ A'' = θ ZE '/ E '' = 0

MZ
D
E’’ E’ Nx
Nous devons donc calculer la rotation relative deA’ par rapport à A’’ suivant l’axe Z.
Lorsque ce déplacement sera nul le moment s’identifiera avec le moment fléchissant dans les
sections A et E.
Y
MZ MZ
G x B
A’’ P/2 P/2A’
Etat isostatique
Equivalent
θ ZA'/ A'' = 0 P
P
F C
Hint1°
P/2 P/2
D
E’’ MZ MZ E’
TD-2A-S4-F412
Compte tenu des axes de symétrie, pour calculer ces rotations, nous allons raisonner sur le
quart du cadre.
P2 Y Etat
P2 isostatique
P1/2 E0‘
G x B
A’’ A’
P1/2
P1/2
P1/2
F C
E0 • E0’{ P1=P ; D2=0 }
P1/2 P1/2
L/2=300mm
P2 P2
D
E
P1/2 P1/2
L/2=300mm
Décomposons l’état isostatique E0’ en 2 états unitaires, afin de calculer D2, à l’aide de
Maxwell-Mohr.

D  =
 δ hC / F 
  =
 f11

f12   P1 
0'
 2  E0 ' θ Z  E0 '  f 21
A '/ A''
f 22   P2  E
0'
f 21 P1
D2 = θ ZA'/ A'' = f 21 P1 + f 22 P2 = 0 ⇒ P2 = −
f 22
(mZ )i (mZ ) j
∑∫ (
L
Calcul des coefficients d’influence fij = )dx à l’aide des intégrales de
barres
0 EIZ
Mohr :
TD-2A-S4-F412
G x B
1/2 A
Z
(mZ)1
Y
Etat E1 -
P2
1/2
C
G x B
A L/4
P1/2
Z
1
¼ Etat E0’ +
P1/2 1 Y (mZ)2
C
1
G x B
A
Z
+
Etat E2
 B((mZ )2 )2 C ((m ) )
2

Calcul de f22 : f22 = 4×  ∫
 dx + ∫ Z 2
dx
 EIZ EIZ 
A B
2 2
1 1
f22 = 4/EIZ + + 4/EIZ +
L/2 L/2
 1 L 4L
f 22 = 8 ×12 ×  f22 =
 EIZ 2 EIZ
 B(m ) (m ) C (m ) (m ) 
Calcul de f12 : f12 = 4 ×  ∫ Z 1 Z 2 dx + ∫ Z 1 Z 2 dx
A EIZ B EIZ 
L/4 1
f12 = 4/EIZ - x +
L/2 L/2
4 1 L L L2
f12 = −  × × 1×  f 12 = −
EI Z 2 4 2 4 EI Z
TD-2A-S4-F412
L2 P1 4 LP2
Dans l’état E0’ D2 = θ ZA'/ A'' = f 21 P1 + f 22 P2 D2 = θ ZA'/ A'' = − +
4 EI Z EI Z
PL PL
Dans l’état E0 θ ZA'/ A'' = 0 ⇒ P2 = ⇒ M ZA =
16 16
Pour trouver le diagramme des moments fléchissants, nous utilisons les résultats obtenus sur
l’état isostatique E0’ au lieu de tout refaire le problème initial.
MZ (E0) = P×mZ (E1)+(PL/16)×mZ (E2)
MZB(E0) = P×mZB(E1)+(PL/16)×mZB(E2) = P×0+(PL/16)×1 = PL/16
MZC(E0) = P×mZC(E1)+(PL/16)×mZC(E2) = P×(-L/4)+(PL/16)×1 = -3PL/16
Y
PL/16
PL/16 +
G x
+ A B +
Z
3PL/16
P F - 3PL/16
- C P
+ E D +
PL/16 + PL/16
Diagramme des moments fléchissants MZ

Le moment fléchissant maximum est de 3PL/16 dans les sections C et F.
3PL 3 × 50 × 600
M ZC = = = 5625mmN MZmax= 5,625 mN
16 16
La contrainte de flexion maximum vaut :
3PL
−
MZ P  h 9 PL 9 × 50 × 600
Y = − 163  ±  = ±
max
σ P
=− =± = ±21,09 MPa •xx = ±21,1 MPa
8 × 25 × 8
xx
IZ b.h  2  8bh 2 2
12
TD-2A-S4-F412
Allongement du cadre :
Il nous faut donc déterminer le déplacement relatif des sections C par rapport à F
horizontalement. Nous allons encore travailler avec l’état E0’.

D  =
 δ hC / F 
  =
 f11

f12   P1 
0'
 2  E0 ' θ Z  E0 '  f 21
A '/ A''
f 22   P2  E
0'
D1 = δ hC / F = f11 P1 + f12 P2
 PL 
Dans l’état E0 δ hC / F = f11P + f12  
 16 
Nous devons donc calculer f11.

 B((mZ )1 )2 C((m ) )
2

Calcul de f11 : f11 = 4 ×  ∫
 dx + ∫ Z 1
dx
 EIZ EIZ 
A B
L/4 2
f11 = 4/EIZ -
L/2
 1 1  L 2 L  L3
f11 = 4 × ×  ×  f11 =
 EIZ 3  4  2  24EIZ
L3 L2  PL  5 PL3 5 × 50 × 6003 × 12
δ hC / F = P− ×  δ hC / F = = = 1,26mm
24 EI Z 4 EI Z  16  192 EI Z 192 × 210.103 × 25 × 83
Allongement du cadre : •hC/F = 1,26 mm
TD-2A-S4-F412
3. MÉTHODE DES ÉLÉMENTS FINIS
PROBLÈME N°17
2A Considérons la poutre à section droite

A variable ci-contre. Le module de YOUNG
1 2 P 3 4 x E vaut 2.105 MPa, L=1 m. La poutre est
encastrée à ses deux extrémités et est
1 2 3 soumise à une charge axiale P=10kN.
L/2 L/2 L/2 La poutre est discrétisée en trois
éléments poutre et quatre nœuds. L’aire
A vaut 10 mm2.
1°) Calculez les matrices de rigidité élémentaire des trois éléments sous la forme :
4.103 [ ] .
2°) Calculez la matrice de rigidité globale de la structure (sous la forme : 4.103 [ ] ).
3°) Calculez les déplacements axiaux des nœuds 2 et 3 (en mm), puis les efforts axiaux
sur les nœuds 1 et 4 (en kN).
4°) Calculez le déplacement axial, la déformation axiale et la contrainte (en MPa) au
milieu de l’élément 2.
Réponses:
1°) [k 1 ] = [k 2 ] = 4.10 3  2 −2 1 −1

 2 
[k ] = 4.10
3 3
− 1 1 
2°)
− 2 
2 −2 0 0 
− 2 4 −2 0 
[K ] = 4.10 3 
0 −2 3 −1
 
0 0 −1 1 
3°) u2 = 0.9375 mm u3 = 0.625 mm Fx1 = −7.5 kN Fx4 = −2.5 kN
4°) u = 0.78125 mm ε = -6.25.10 -4 σ = -125MPa
RÈPONSES N°17
La poutre droite à section variable est encasrée à ses deux extrémités. La charge P étant
axiale, seules les réactions horizontales existent, et le problème est donc hyperstatique
externe de degré 1.
TD-2A-S4-F412
2
1
3
2 P
L/2 3
2A
L/2 4
A
L/2 x
1)
Discrétisons la poutre en 3 éléments et 4 nœuds . Dans notre cas les axes locaux et globaux
coïncident. Chaque élément, compte tenu du chargement, n’est soumis qu’à des efforts
normaux. Chaque nœud admet un déplacement u suivant x.
EA  1 − 1
La matrice de rigidité élémentaire en traction-compression k e = [ ]
L − 1 1 
.relie les
 Fxe1  EA  1 − 1  u1 
efforts d’extrémité de l’élément et les déplacements des nœuds  e  =
 Fx 2  L  − 1 1  u2 
1
2
2
2
L/2 2A P 3
3
3
2A
L/2 4
A x
L/2
TD-2A-S4-F412
Calculons les matrices de rigidité élémentaires.
Elément 1 :
− 1
[k ] = E2A
1 1
L − 1 1 

2 1
1
× 2 × 10  1 − 1
[ ] 2.10 5 F x1
k1 = − 1 1 
1000   1
2
− 2 F1x2
[ ] 
k 1 = 4.103 
2
 2
− 2 2 
L/2 2A x
 Fx11  3 2 − 2  u1 
 1  = 4.10   
 Fx 2   − 2 2  u2 
Elément 2 :
− 1
[k ] = E2A
2 1
L − 1 1 
 2
2 F2x2
2.105 × 2 × 10  1 − 1
[ ]
k2 =
1000 − 1 1  2
 
2
3
F2x3
 2 − 2
[ ]
k 2 = 4.103   2A
x
− 2 2  L/2
 Fx22  3 2 − 2 u2 
 2  = 4.10   
 Fx 3   − 2 2  u3 
Elément 3 :
− 1
[k ] = EA
3 1
L − 1 1 
 F3x3 3
2
3
2.105 × 10  1 − 1
[ ]
k3 =
1000  − 1 1 
4 F3x4
2
A
 1 − 1
[ ]
k 3 = 4.103   L/2
x
− 1 1 
 Fx33  3 1 − 1 u3 
 3  = 4.10   
 Fx 4  − 1 1  u4 
TD-2A-S4-F412
2)
Détermination de la matrice de rigidité K de la structure :
Procédons à l’assemblage des 3 matrices de rigidité élémentaires :
− 2
[k ] = 4.10 −22
1 3
2 
[k ] = 4.10 −22
2 3 − 2
[k ] = 4.10 −11
3 3 − 1
  2   1 
 2 −2 0 0
− 2 4 − 2 0 
[K ] = 4.103  
 0 − 2 3 − 1
 
0 0 −1 1 
3)
 FX 1   2 −2 0 0   u1 
  − 2 4 − 2 0  u 
D’où : [P ] = [K ][D ]  FX 2  = 4.103   2 
 FX 3   0 − 2 3 − 1 u3 
    
 FX 4  0 0 − 1 1  u4 
Appliquons les conditions limites :
 FX 1   2 −2 0 0  0 
  − 2 4 − 2 0  u 
 + P  = 4.103   2 
 0   0 − 2 3 − 1 u3 
    
 FX 4  0 0 −1 1   0 
La résolution des équations 2 et 3 du système nous permet de trouver
 P = 4.103 (4u 2 − 2u3 )


0 = 4.10 (− 2u2 + 3u3 )
3
u 2 = 0.9375mm u3 = 0.625mm
Ces valeurs reportées dans le système nous permettent de résoudre les équations 1 et 4 :
7,5 kN
FX 1 = 4.103 (− 2u2 ) = 4.103 (− 2 × 0,9375)

 FX 4 = 4.10 (− u3 ) = 4.10 (− 0,625)
3 3
P=10 kN
L/2 2,5 kN Fx1 = −7.5kN Fx4 = −2.5kN
L/2
x
L/2
TD-2A-S4-F412
4)
Déterminons les efforts d’extrémités de l’élément 2 :

2
2
 Fx22  − 2 u2  F
3 2
x2
 2  = 4.10   
 − 2 2  u3 
2
 Fx 3 
3
F2x3
F  2
3 2 − 2 0,9375
  = 4.10  x
x2
 
 − 2 2   0,625  2A
2
F  x3
L/2
 Fx22  3 2 − 2 0,9375  2 × 0,9375 − 2 × 0,625   2500 

 2  = 4.10     = 4.103  = 
 Fx 3   − 2 2   0,625  − 2 × 0,9375 + 2 × 0,625 − 2500
L’élément est donc soumis à un effort de compression de 2,5 kN
La contrainte vaut pour tout point: 2500 N 2
N x − 2500 2
σ xx = = = −125MPa
A 2 × 10
3
2500 N
La déformation associée vaut :
2A x
L/2
σ xx − 125
ε xx = = 5
= 625.10 −6
E 2.10
La matrice d’interpolation des déplacements [A] en traction-compression est :
[A] = 1 − X X
 L L 
[ ]
u ( x) = [ A] D 2
 X
u ( x) = 1 −
X  u2 
L  u3 
 L
L 1 1  u2 
u( ) =  = 0,78125mm
2 2 2  u3 
TD-2A-S4-F412
PROBLÈME N°18
1 Avant d'appliquer la charge horizontale P un jeu "j" existe entre
2
l'extrémité droite de la poutre et l'appui de droite. La poutre est
1 2
P 3 x modélisée par deux éléments "poutre" et trois nœuds.
A 1°) Déterminer sous forme littérale les matrices de rigidité

4A
j élémentaires et la matrice de rigidité totale en fonction de L, E et
L L/2
A (mettre 2EA/L en facteur).
2°) E=200GPa, A=5cm2, L=1m et P=100kN. Calculer:
a) Le jeu "j" pour que la poutre arrive juste en contact avec la liaison de droite.
b) Le jeu "j" pour qu'un l'effort de compression de 10kN s'exerce dans l'élément 2.
Réponses:
 2 −2 0 
[k ] = 2EA −22 −22[k 2 ] = 2EA −11 −11
1
[K ] =
2EA 
− 2 3 − 1 j = 0,25mm j = 0,175mm
L   L   L 
 0 − 1 1 
RÈPONSES N°18
2
P
L 4A
3
A
L/2
j x
1°)
Discrétisons la poutre en 2 éléments et 3 nœuds . Dans notre cas les axes locaux et globaux
coïncident. Chaque élément, compte tenu du chargement, n’est soumis qu’à des efforts
normaux. Chaque nœud admet un déplacement u suivant x.
EA  1 − 1
La matrice de rigidité élémentaire en traction-compression k e = [ ]
L − 1 1 
.relie les
 F e  EA  1 − 1  u1 
efforts d’extrémité de l’élément et les déplacements des nœuds  xe1  =
 Fx 2  L  − 1 1  u2 
TD-2A-S4-F412
Elément 1 :
− 1 1
[k ] = E4A
1 1
L − 1 1  F 1
x1

1
2EA  2 − 2 
[ ]
k1 =
L  − 2 2 
F1x2
 Fx11  2EA  2 − 2  u1  2
 1=
L − 2 2  u2  4
L x
 Fx 2 
Elément 2 :
F2x2 2
EA  1 − 1
[ ]2
k =
L  − 1 1 
2 2
2 EA  1 − 1
[ ]
k2 =
L  − 1 1 
3
F2x3
A x
L/2
 Fx22  2 EA  1 − 1 u 2 
 2=
 Fx 3  L − 1 1  u3 
− 2
[k ] = 2EA
1 2
L − 2 2 
[k ] = 2ELA −11
2 − 1
  1 
 2 −2 0 
2EA 
[K ] =  − 2 3 − 1
L
 0 − 1 1 
1°)  FX 1   2 − 2 0   u1 
  2EA   
[P ] = [K ][D ]  FX 2  = L − 2 3 − 1 u2 
D’où :  FX 3 
   0 − 1 1  u3 
TD-2A-S4-F412
a) Jeu "j" pour que la poutre arrive juste en contact avec la liaison de droite.
 FX 1   2 − 2 0  0 
  2 EA   
+ P = L − 2 3 − 1 u2 
 0   0 − 1 1  u3 
 
La résolution des équations 2 et 3 du système nous permet de trouver u2 et u3.

 P = L (3u2 − u3 )
2 EA
 2 × 200.103 × 5.102

100.103 =

(3u2 − u3 ) u2 = u3 = 0.25mm
1000
0 =
2 EA
(− u2 + u3 )  − u2 + u3 = 0
 L
J=0,25 mm
b) Le jeu "j" pour qu'un effort de compression de 10kN s'exerce dans l'élément 2.
 
 FX 1  2EA  2 −2 0  0 
+P= − 2 3 − 1  u 2 
 P  L   
−   0 −1 1   u 3 
 10 
La résolution des équations 2 et 3 du système nous permet de trouver u2 et u3.
  PL 
 P = L (3u2 − u3 )
2 EA 9 PL
 2 EA = 3u2 − u3  u2 =
40 EA
 P 2 EA  PL 
− = (− u 2 + u3 ) − = −u2 + u3 u3 =
7 PL
= 0,175mm
 10 L  20 EA  40 EA
J=0,175 mm
TD-2A-S4-F412
PROBLÈME N°19
Ce problème est le même que le problème 9 qui a été résolu par la méthode de
MAXWELL-MOHR.
Y Les barres (1) et (2) de même section droite A et de même

matériau) sont articulées en 1, 2 et 3. Une force verticale
3 P est appliquée en 2.
1°) Déterminez les matrices de rigidité élémentaire des

deux éléments en fonction de E, A et L et dans le repère
global X,Y. Écrire ces matrices sous la forme :
2
1
135° [k ] = EA
1
L
[] [k ] =
2 2 EA
4 L
[ ]
1 2
X
2°) Déterminez sous forme littérale la matrice de rigidité

P
L de la structure dans le repère global X, Y. Écrire cette
matrice sous la forme :
[K ] = EA [ ]
L
3°) Calculez (en mm) le déplacement horizontal et le déplacement vertical du noeud 2 pour
A= 1cm², E= 200 GPa, L= 2m, P= 10kN et la valeur (en kN) des réactions sur les nœuds 1 et 3.
4°) Calculez les efforts axiaux dans les éléments 1 et 2.
Réponses :
1°)
1 0 −1 0 1 −1 −1 1 
 0 0   − 1 
2 EA  − 1
[k ]1
=
EA  0
L − 1
0
0 1 0
[k ]
2
=
4 L − 1 1
1
1
1
−1 
   
0 0 0 0 1 −1 −1 1 
2°)
 1 0 −1 0 0 0
 0 0 0 0 0 0
 
− 1 2 2 2 2
0 (1 + ) − −
 4 4 4 4 
EA  2 2 2 2
[K ] =  0 0 − − 
L  4 4 4 4 
 2 2 2 2
 0 0 − −
 4 4 4 4 
 0 2 2 2 2
0 − −
 4 4 4 4 
TD-2A-S4-F412
3°) uX 2 = −1 mm uY 2 = −3 ,83 mm FX 1 = −10 kN FX 3 = −10 kN

FY 1 = 0 FY 3 = 10 kN
3°) Nx = −10 kN
AB
NxBC = 14.142 kN
RÈPONSES N°19
1
135°
1 2
X
P
L
1°) Matrices de rigidité élémentaire des deux éléments en fonction de E, A et L et

dans le repère global X,Y.
Le problème est à 2 degrés de liberté, uX suivant X et vY suivant Y ( X,Y étant le repère

global )
Elément 1 :
Dans les axes locaux xy
EA  1 − 1
[k ]
1
/xy =
L  − 1 1 
Y
y
 Fx11  EA  1 − 1  u x1 
 1 =
 Fx 2  /xy L − 1 1  u x 2 
1
Dans les axes globaux XY 1 2
1 0 − 1 0 x
 0 0 0  F 1
F 1
X
[k ] EA  0 x1 x2
1
= L
L − 1 0 1 0
/XY
 
0 0 0 0
TD-2A-S4-F412
Elément 2 :
Dans les axes locaux xy x Y
 1 − 1
[k ] 2
/xy =
EA
− 1 1 
L 2  F2x3
 Fx22  EA  1 − 1 u x 2 
 2 =    3
 Fx 3  /xy L 2 − 1 1  u x 3 
2
Dans les axes globaux XY

Nous devons calculer [λ] la matrice de
transformation,car pour l’élément 2, L 2
les axes locaux et globaux ne coïncident pas :
cosθ sin θ
[λ ] = 
0 0 
 0 0 cosθ sin θ 
cos 135° sin 135°

[λ ] = 
0 0  135°
 0 0 cos135° sin 135° 2 2
X
F x2
 − 1 1 0 0
[λ ] = 2
 0 0 − 1 1
y
2  
− 1 0 
 
[ ] [ ]
k 2 X ,Y = [λ ] k 2 x , y [λ ]
t
[ ]k 2 X ,Y =
2  1 0  EA  1 − 1 2 − 1 1
2  0 − 1 L 2 − 1 1  2  0 0
0 0
− 1 1
 
0 1
− 1 0 
 
EA 2  1 0   1 − 1 − 1 1 0 0
[ ]
k 2 X ,Y =
4L  0 − 1 − 1 1   0 0 − 1 1
 
0 1
− 1 1   1 −1 −1 1
 1 − 1 − 1 1 0 0  1 − 1
2 EA − 1 1
[ ]
k 2 X ,Y =
EA 2  
4L  1 − 1  0 0 − 1 1

[ ] k 2
X ,Y =
4 L − 1 1 1 − 1
   
− 1 1   1 −1 −1 1 
 2 2 2 2 
 − − 
 4 4 4 4 
− 2 2 2
−
2
[ ]
k 2 X ,Y =
EA
L
 4
 2
4
2
4
2
4 
2
− − 
 4 4 4 4 
 2 2 2 2 
 4 − −
4 4 4 
TD-2A-S4-F412
2°) Matrice de rigidité de la structure dans le repère global X, Y.

 2 2 2 2 
 − − 
 4 4 4 4 
1 0 −1 0  2 2 2 2
EA − 4 −
[k ]1
=

EA  0 0 0 0 [k ]
2
X ,Y =
L  2
4
2
4
2
4 
2
L − 1 0 − 
/XY
0 1 −
   4 4 4 4 
0 0 0 0  2 2 2 2 
 − − 
 4 4 4 4 
 1 0 −1 0 0 0
 0 0 0 0 0 0
 
− 1 0 (1 +
2
) −
2
−
2 2
 4 4 4 4 
EA  2
[K ]/ XY =  0 0 −
2 2 2
− 
L  4 4 4 4 
 2 2 2 2
 0 0 − −
4 4 4 4 
 
 0 2 2 2 2
0 − −
 4 4 4 4 
3°) Déplacement horizontal et déplacement vertical du noeud 2 :

 1 0 −1 0 0 0
 0 0 0 0 0 0  u 
 FX 1  
F  − 1 2 2 2 2  X 1 
0 (1 + ) − − v
 Y1   4 4 4 4  Y1 
 FX 2  EA  2  u X 2 
[P ] = [K ][D ]   =  0 0 −
2 2 2
−  
 FY 2  L  4 4 4 4   vY 2 
 FX 3   2 2 2 2  u X 3 
   0 0 − −
4 4 4 4   v 
 FY 3  / XY    Y 3  / XY
 0 2 2 2 2
0 − −
 4 4 4 4 
 1 0 −1 0 0 0
 0 0 0 0 0 0 0
 FX 1    
 F  − 1 0 (1 +
2
) −
2
−
2 2  
0
 Y1   4 4 4 4  
 0  EA  2 2 2 2  u X 2 
 4
=  0 0 − −  
 − 10  L  4 4 4 4   vY 2 
 FX 3   2 2 2 2  0 
   0 0 − −
4 4 4 4   0 
 FY 3  / XY     / XY
 0 2 2 2 2
0 − −
 4 4 4 4 
TD-2A-S4-F412
La résolution des équations 3 et 4 du système nous permet de trouver uX2 et vY2.
 EA  2 2    
0=  (1 + )u X 2 − vY 2   0 =  (1 +
2
)u X 2 −
2
vY 2 
 L  4 4    4 4 
 
− 10 4 = EA  − 2 u + 2   × ×
3
 L  4
X2 vY 2 
− 10 =
4 200 . 10 100 2
(− u X 2 + vY 2 )
 4  2.103 × 4
 2 2
 0 = (1 + )u X 2 − vY 2  u X 2 = −1mm
4 4
 
 − 10 × 2.10 × 4 = −u X 2 + vY 2 vY 2 = −3,828mm
4 3
 200.103 × 100 × 2
•H2= -1 mm
•V2= -3,83 mm
Réactions sur les nœuds 1 et 3 :
 1 0 −1 0 0 0
 0 0 0 0 0 0 0
 FX 1    
 F  − 1 0 (1 +
2
) −
2
−
2 2 
0 
 Y1   4 4 4 4  
 0   2 2 2 2  −1 
 4
= 10 4  0 0 − −  
 − 10   4 4 4 4   − 3,83
 FX 3   2 2 2 2  0 
   0 0 − −
 FY 3  / XY 4 4 4 4   0 
   / XY
 0 2 2 2 2
0 − −
 4 4 4 4 
La résolution des équations 1, 2, 5 et 6 du système nous permet de trouver les réactions.
 1 0 −1 0 0 0
 0 0 0 0 0 0 0
 FX 1      FX 1 = 104
 F  − 1 2  
0  
2 2 2
0 (1 + ) − − FY 1 = 0
 Y1   4 4 4 4 
 0   2 2 2 
2   − 1   F = 10 4  2 − 3,83 2 
 4
= 10 4  0 − −    X3  4 4 
0
 − 10   4 4 4 4  − 3,83  
 FX 3   2 2 2 2  0    2 2 
− − 
 0 4   0  FY 3 = 10  − 4 + 3,83 4 
  0 4
4 4 4
 FY 3  / XY    / XY   
 0 2 2 2 2
0 − −
 4 4 4 4 
TD-2A-S4-F412
10 kN
3
10 kN FX1= 10 kN
2
FY1= 0 kN
FX3= -10 kN
FY3= 10 kN
1
135°
10 kN 1 2
X
10 kN
L
4°) Efforts axiaux dans les éléments 1 et 2 :
Elément 1 :
Y
EA  1 − 1
[k ]
1
/xy =
L  − 1 1  y
 Fx11  EA  1 − 1  u x1  1
 1 =
 Fx 2  /xy L − 1 1  u x 2  1 2
x
1 1
Les axes locaux et globaux coïncidant : F x1 F x2 X
L
uX1=ux1 et uX2=ux2
 Fx11  EA  1 − 1  0  4 1 − 1  0   10 4 
 1 = =  − 1 1  − 1 =  4
L − 1 1   − 1
10
 Fx 2  /xy      − 10 
L’effort normal dans l’élément 1 est un effort axial de compression :
y
1
N 1x = −10kN 1 2
x
10 kN 10 kN
L
TD-2A-S4-F412
x Y
Elément 2 :
F2x3
−
[ ]
k 2 /xy =
EA 

1 1

L 2 − 1 1  3
 Fx 2 
2
EA  1 − 1 u x 2  2
 2 =   
 Fx 3  /xy L 2 − 1 1  u x 3 
Les axes locaux et globaux ne coïncident pas

L 2
il nous faut calculer :
ux2 et ux3
Utilisons la relation : D[ ]e
x, y [ ]
= [λ ] D e X ,Y
Avec [λ] trouvée précédement :

135°
− 1 1 0 0
[λ ] = 2 
X
2
2  0 0 − 1 1
2
F x2
u X 2  −1 
y
   
u x 2  2 − 1 1 0 0  vY 2  2 − 1 1 0 0  − 3,828 2 − 2,828
u  = 2  0 0 − 1 1 u  = 2  0 0 − 1 1  0  = 2  0 
 x3    X3    
   
 vY 3   0 
 Fx22  EA  1 − 1 u x 2  10 4  1 − 1 2 − 2,828  1 − 1 − 2,828
D’où :  2  =    =     = 5.103   
 Fx 3  /xy L 2  − 1 1  u x 3  2 − 1 1 2  0  − 1 1   0 
x
14,14 kN
 Fx22  − 14140 3
 2 = 
 Fx 3  /xy + 14140 2
L 2
L’effort normal dans l’élément 1 est un effort axial de traction :
N xBC = 14,14kN
2
y 14,14 kN
TD-2A-S4-F412
PROBLÈME N°20
Ce problème est le même que le problème 12 qui a été résolu par la méthode de
MAXWELL-MOHR.
P=12kN, L=1m, E=200GPa, la section est

y
un rectangle de 80mmx37.5mm.
L L
y
P 1°) Déterminez sous forme numérique (en
z 1 2 3 x N/m) la matrice de rigidité de la poutre
G
1 2 ci-contre (discrétiser en deux éléments
de longueur L et trois nœuds).
2°) En déduire les valeurs numériques de v2 (en mm), de θz2, θz3 (en rd) et des réactions sur
les nœuds 1 et 3 (en kN pour les forces, Nm pour les moments). Tracez le diagramme des
moments fléchissants sur les éléments 1 et 2.
3°) Calculez la flèche (en mm), la déformation axiale maximum et la contrainte de flexion
maximum (en MPa) au milieu de l’élément 2
RÈPONSES N°20
L=1m
Y P=12kN
VA L=1m
MA
HA A b =37,5mm
G
x VB
C
Z h =80mm
B
Le problème est hyperstatique externe de degré 1.
TD-2A-S4-F412
1°)
Discrétisons la poutre en 2 éléments et 3 nœuds :
FY1 L=1m
P=12kN L=1m
MZ1
Elément 1 FY3
Z
Nœud 1 Elément 2
X
Nœud 2
Nœud 3
Dans notre cas les axes locaux et globaux coïncident.
Chaque élément, compte tenu du chargement, n’est soumis qu’à un cisaillement de flexion.
Chaque nœud admet un déplacement v suivant Y, et une rotation • suivant Z .
La matrice de rigidité élémentaire en flexion dans le plan XY :
 12 6 L − 12 6 L 
 6 L 4 L2 − 6 L 2 L2 
[ ] EI
k e = 3z 
L − 12 − 6 L 12 − 6 L 
.
 2 
 6L 2L − 6L 4L 
2
Elle relie les efforts d’extrémité de l’élément et les déplacements des nœuds
 Fye1   12 6 L − 12 6 L   v1 
 e  6L 4L2 − 6 L 2 L2  θ 
 M z1  = EI z    z1 
 Fy 2  L − 12 − 6L 12 − 6 L  v2 
e 3
 e  2  
 M z 2   6L 2L − 6 L 4 L  θ z 2 
2
TD-2A-S4-F412

Y
Elément 1 :
6 L − 12 6 L  F1Y1
 12 L=1m
 6 L 4 L2 − 6L 2 L2 
[ ]
1
EI M F1Y2
k 1 = 3z   Z1
L − 12 − 6L 12 − 6 L
 2  1 M1Z2
 6 L 2 L − 6L 4 L 
2
Elément 1
2
Z X
 12 6 − 12 6   12 6 − 12 6 
.10 × (80.10 )
 6 4 −6 2    4 − 6 2 
[k ] = 200.10 × 37,512 [ ]
9 −3 −3 3
 5 6
1
; k = 3,2.10
1
− 12 − 6 12 − 6 − 12 − 6 12 − 6
   
 6 2 −6 4   6 2 −6 4 
 Fy11   12 6 − 12 6   v1 
 1  6
 M z1  = 3,2.105  4 − 6 2  θ z1 
 Fy12  − 12 − 6 12 − 6  v2 
 1    
 M z 2   6 2 − 6 4  θ z 2 
Elément 2 :
 12 6 L − 12 6 L  Y
 2 
EI 6 L 4 L − 6 L 2 L 
[ ]
2
k 2 = 3z 
L − 12 − 6 L 12 − 6 L F2Y2
 2  L=1m
 6L 2L − 6L 4L 
2
2
M Z2 F2Y3
2 M2Z3
 12 6 − 12 6  Elément 2
 6 4 − 6 2  3
[ ]
k = 3,2.10
2 5
− 12 − 6 12 − 6
Z X
 
 6 2 −6 4 
 Fy22   12 6 − 12 6   v2 
 2  6
 M z 2  = 3,2.105  4 − 6 2  θ z 2 
 Fy23  − 12 − 6 12 − 6  v3 
 2   
 M z 3   6 2 − 6 4  θ z 3 
TD-2A-S4-F412
 12 6 − 12 6   12 6 − 12 6 
 6  6 4 − 6 2 
4 − 6 2 
[ ]
k = 3,2.10
1 5
[k ]
2
= 3,2.105 
− 12 − 6 12 − 6
− 12 − 6 12 − 6
   
 6 2 −6 4   6 2 −6 4 
 12 6 − 12 6 0 0
 6 4 −6 2 0 0 

− 12 − 6 24 0 − 12 6 
[K ] = 3.2 105  
 6 2 0 8 −6 2 
 0 0 − 12 − 6 12 − 6
 
 0 0 6 2 − 6 4 
D’où : [P ] = [K ][D ]  FY 1   12 6 − 12 6 0 0   v1 
M   6 4 −6 2 0 0  θ Z 1 
 Z1  
 FY 2  
5 − 12 − 6 24 0 − 12 6   v2 
  = 3.2 10   
M Z 2   6 2 0 8 − 6 2  θ Z 2 
 FY 3   0 0 − 12 − 6 12 − 6  v3 
    
 M Z 3   0 0 6 2 − 6 4  θ Z 3 
2°)
 FY 1   12 6 − 12 6 0 0  0  (1)
M   6 4 −6 2 0 0   0  (2)
 Z1  
−P 
5 − 12 − 6 24 0 − 12 6   v2  (3)
  = 3.2 10   
 0   6 2 0 8 − 6 2  θ Z 2  (4)
 FY 3   0 0 − 12 − 6 12 − 6  0  (5)
    
 0   0 0 6 2 − 6 4  θ Z 3  (6)
La résolution des équations (3), (4),et (6) permet de trouver les déplacements :
v2 = −2.7344mm θ z 2 = −1.1719 10−3 θ z 3 = +4.6875 10−3 rd
TD-2A-S4-F412
 FY 1   12 6 − 12 6 0 0  0 
 M   6 4 −6 2 0 
0  0 
 Z1   
− 12000  − 12 − 6 24 0 − 12 6  − 2,7344.10 −3

 = 3.2 10 
5
  −3 
 0   6 2 0 8 − 6 2   − 1,1719.10 
 FY 3   0 0 − 12 − 6 12 − 6  0 
    
 0   0 0 6 2 − 6 4   4,6875.10−3 
La résolution des équations (1), (2),et (5) permet de trouver :
Fy1 = 8.25kN Fy 3 = 3.75kN M z1 = 4500Nm
Y
P=12kN
4500 Nm
1 2,7344 mm
Z 2 3
-3 x
4,6875.10 rd
8,25 kN
1,1719.10-3 rd
3,75 kN

Y
 Fy11   12 6 − 12 6   v1 
 1 F1Y1
 6 4 − 6 2  θ z1 
L=1m
 M z1  = 3,2.105  1
 Fy12  − 12 − 6 12 − 6  v2  M Z1 F1Y2
 1    
 M z 2   6 2 − 6 4  θ z 2  1 M1Z2
Elément 1
2
F 1
y1
  12 6 − 12 6   0 Z X
   6
M
1
z1  = 3,2.105  4 − 6 2   0  −3
10
F 1  − 12 − 6 12 − 6 − 2,7344
    
y2
 M
1
z2   6 2 − 6 4   − 1,1719 
TD-2A-S4-F412
 Fy11 = 3,2.105 (− 12 × (− 2,7344) + 6 × (− 1,1719))10−3  Fy11 = 8250 N

 1  1
 M z1 = 3,2.10 (− 6 × (− 2,7344) + 2 × (− 1,1719))10
−3
 M z1 = 4500mN
5
 1  1
 Fy 2 = 3,2.10 (12 × (− 2,7344) − 6 × (− 1,1719))10
−3
 Fy 2 = −8250N
5
M 1z 2 = 3,2.105 (− 6 × (− 2,7344) + 4 × (− 1,1719))10−3 M 1z 2 = 3750mN

Y
F1Y1=8250 N
F1Y2=-8250 N
1
M1Z1=4500 mN M1Z2=3750 mN
Elément 1
2
Z X
- + 3750 Nm
4500 Nm
(MZ)

Y
F 2
y2
  12 6 − 12 6   v2 
 
M
2
z2
 6
 = 3,2.105  4 − 6 2  θ z 2  F2Y2 L=1m
F 2  − 12 − 6 12 − 6  v3 

y3
    M2Z2 F2Y3
 M
2
z3   6 2 − 6 4  θ z 3 
2 M2Z3
 Fy22   12 6 − 12 6  − 2,7344Z Elément 2
3
 2  4 − 6 2   − 1,1719  −3
X
 M z 2  = 3,2.105  6 10
 Fy23  − 12 − 6 12 − 6  0 
 2   
 M z 3   6 2 − 6 4   4,6875 
 Fy22 = 3,2.105 (12 × (− 2,7344) + 6 × (− 1,1719) + 6 × (4,6875))10−3  Fy22 = −3750 N

 2  2
 M z 2 = 3,2.10 (6 × (− 2,7344) + 4 × (− 1,1719) + 2 × (4,6875))10
−3
M z 2 = −3750 mN
5
 2  2
 Fy 3 = 3,2.10 (− 12 × (− 2,7344) − 6 × (− 1,1719) − 6 × (4,6875))10
−3
 Fy 3 = 3750 N
5
 M z23 = 3,2.105 (6 × (− 2,7344) + 2 × (− 1,1719) + 4 × (4,6875))10−3  M z23 = 0

Y
F2Y2=-3750 N
F2Y3=3750 N
M2Z2=-3750 mN 2
Elément 2
Z 3
X
3750 Nm +
(MZ)
TD-2A-S4-F412
3°)
Flèche :
La matrice d’interpolation des déplacements vaut :
[A] = 1 − 3x2 x2 x3 

[ ]
2
2 x3 2 x 2 x3 3x 2 2 x3
+ x− + 2 − 3 − + v( x) = [A] D e
 L L3 L L L2 L L L2 
 v2 
 
 3x 2 2 x3 2x2 x3 3x 2 2 x3 x 2 x 3  θ z 2 
v( x) = 1 − 2 + 3 x− + 2 − 3 − + 2
 L L L L L2 L L L   v3 
 
θ z 3 
Pour x =
L
et L = 1m [A] =  1 1 1 1
− 
2 2 8 2 8
 − 2.7344 
 
1 1 1 1 1   − 1.1719  −3
v( ) =  −  10 = 2,1 .10 −3 m = 2.1mm
2 2 8 2 8 0 
 
 + 4.6875
Déformation:
La matrice d’interpolation des déformations vaut :
[B] = −6 12x

+ 3 −
4 6x
+
6 12x
− −
2 6x 
+
 L
2
L L L2 L2 L3 L L2 
D’où :
 v2 
 
 6 12 x 4 6x 6 12 x 2 6 x  θ z 2 
ε xx = − y − 2 + 3 − + − − + 2
 L L L L2 L2 L3 L L   v3 
 
θ z 3 
L
Pour : x= L = 1m y = ±40mm
2
− 2.7344
 − 1.1719 
ε xx = ±0.04[0 -1 0 1]  10−3 ε xx = ±2,3466 .10 −4
 0 
 
+ 4.6875
Contrainte :
σ xx = Eε xx = ±46.9 MPa
TD-2A-S4-F412
4. PROBLEMES DE SYNTHESE
PROBLÈME N°21
1ère partie L/2 L
P1
A B x
Considérons la poutre droite ABC hyperstatique ci-dessus. Toutes les liaisons sont des
articulations. La section droite de la poutre est constante et son moment d’inertie de
flexion est noté IY . Le module de YOUNG du matériau est E. La poutre est soumise à
un couple de flexion P1 appliqué en B.
Prenez comme inconnue hyperstatique la réaction en C (une force P2 dirigée vers le

haut ↑ ).
La méthode de résolution est au choix celle de CASTIGLIANO, ou celle de MAXWELL-
MOHR. Dans les deux méthodes l’énergie de déformation due à l’effort tranchant sera
négligée devant l’énergie de déformation due au moment fléchissant. Si vous choisissez
la méthode de MAXWELL-MOHR écrivez l’expression littérale de la matrice de
flexibilité sous la forme :
[f] = 1 [ ]
EIY
Quelle que soit la méthode choisie, CASTIGLIANO ou MAXWWELL-MOHR :
1°) Indiquez quelles sont, selon-vous, les différences essentielles entre les méthodes
de CASTIGLIANO et de MAXWELL-MOHR
2°) Déterminez l’expression littérale de l’inconnue hyperstatique P2 en C en fonction de
P1 et de L en précisant bien son sens ( ↑ ou ↓ ). Calculez sa valeur numérique en Newton si
P1= 144 Nm et L=1 m.
3°) Déterminez l’expression littérale de la rotation en B en fonction de P1, L, E et IY.
Calculer sa valeur numérique en radian si la section droite a un moment d’inertie IY=
0.25 cm4 et si E = 200 GPa. Y
L/2 L
ème
2 partie
Elément 1 Elément 2
X
C
Nœud 1 Nœud 2 Nœud 3
TD-2A-S4-F412
La poutre de la 1ère partie est discrétisée en deux éléments "poutre" et trois nœuds. La
longueur L vaut 1m, la section droite de la poutre est constante et son moment d’inertie
de flexion Iz vaut 0.25 cm4. Le module de YOUNG E du matériau vaut 200 GPa. La
poutre est soumise au couple de flexion C= 144 Nm appliqué sur le nœud 2.
Toutes les liaisons sont des articulations.
On donne la matrice de rigidité d’un élément de poutre de longueur L en flexion

dans le plan XY :
 12 6L − 12 6L 
 6L 4L2 − 6L 2L2 
[k e ] = EIz 
L − 12 − 6L 12 − 6L 
3


 
 6L 2L
2
− 6L 4L2 
Attention ! La longueur de l’élément 1 est L/2 !
1°) Calculez numériquement (utilisez les unités SI : le mètre et le Newton) les

matrices de rigidité élémentaire des deux éléments dans le repère global X,Y sous la
forme :
[k ] = EI [ ]
e
3
z
L
2°) Calculez numériquement (utilisez les unités SI : le mètre et le Newton) la matrice
de rigidité de la structure dans le repère global X,Y sous la forme :
[K ] = EI3z [ ]
L
3°) Calculez, dans le repère global X,Y les rotations en radian des nœuds 1, 2 et 3 et
les réactions en Newton s’exerçant sur les nœuds 1, 2 et 3.
4°) Calculer en mm la flèche de la section droite située au milieu de l’élément 2.

Calculez d’abord numériquement la matrice d’interpolation des déplacements [A],
puis la flèche.
On donne la matrice d’interpolation des déplacements d’un élément de poutre de

longueur L en flexion dans le plan XY :
 3x 2 2x 3 2x 2 x3 3x 2 2x 3 x2 x3 
[A] = 1 − 2 + 3 x − + 2 − − + 2
 L L L L L2 L3 L L 
Réponses:
1ère Partie
1°) La méthode de CASTIGLIANO est analytique, celle de MAXWELL-MOHR matricielle.
P1
2°) P2 = − = −48 N ↓
3L
TD-2A-S4-F412
P1 L
3°) D1 = = 0.032 rd
9 EI I
2ème partie
1°)
 96 24 − 96 24   12 6 − 12 6 
 24 − 24   6 −6 2 
[k ]
1
= 500 
8
− 96 − 24 96
4  2
− 24 
[ ]
k = 500 
− 12
4
−6 12 − 6
   
 24 4 − 24 8   6 2 −6 4 
2°)
 FY 1   96 24 − 96 240  0 
0
0   24 8 − 24 4 0 θZ 1 
0
    
FY 2  − 96 − 24 108 − 18 − 12 6  0 
  = 500   
144  24 4 − 18 12 − 6 2 θZ 2 
 FY 3   0 0 − 12 − 6 12 − 6  0 
    
 0   0 0 6 2 − 6 4 θ Z 3 
3°)
 0  8 4 0  θ Z 1 
    
144  = 500 4 12 2  θ Z 2 
 0  0 2 4  θ Z 3 
θ Z 1 = θ Z 3 = −0.016 rd
θ Z 2 = 0.032 rd
 FY 1 = 192   96 24 − 96 24 0 0  0 
 0   24 8 − 24 4 0 0  − 0.016
  
FY 2 = −144 − 96 − 24 108 − 18 − 12 6   0 
  = 500   
 144   24 4 − 18 12 − 6 2   0.032 
 FY 3 = −48   0 0 − 12 − 6 12 − 6  0 
    
 0   0 0 6 2 − 6 4  − 0.016
4°)
v ( x ) = [A ][D e ]
 2
2x3 2x2 x3 3x2 2x3 x2 x3 
[A ] = 1 − 3 x2 + x − + 2 − 3 − + 2 
 L L3 L L L2 L L L 
L
x = L = 1m
2
 0 
[A ] = 
1 1 1 1  0.032 
−  1 1 1 1 1
2 8 2 8 v( ) = −   = 0.006 m = 6 mm
2 2 8 2 8  0 
 
 − 0. 016 
TD-2A-S4-F412
RÈPONSES N°21
1ère partie
L/2 L
P1 Z Y
Etat E0
G x
A B C
1°) La méthode de CASTIGLIANO est analytique, celle de MAXWELL-MOHR matricielle.
L/2 L
P1 Z Y P2
Etat E0’
G x
A B C
L/2 L
P1=1 Z Y
Etat E1
G x
B
A C
(mY)E1
- 1
L/2 L
Z Y P2=1
Etat E2
G x
A B C
-
(mY)E2
L
TD-2A-S4-F412
L L2 
[D ]E = [ f ]E0 ' [P ]E0 '  D1 
0 =
1 6
 6   P1 
0'
  E0 ' EI Y  L2 L3   P2  E
 6 2 
0'
ou
2
P L P P L2 P 2 L3
Wdéf = + 1 2 + 2
1
12 EIY 6 EIY 2 EIY
2°) P2 = − P1 = −48 N ↓
3L
3°) D1 = P1 L = 0.032rd
9 EI I
ème
2 partie
Y
L/2 L
Nœud 1 Nœud 2 Nœud 3

C
X
A B C
Nœud 1 Y Nœud 2 Nœud 3
C=144 mN
•Z1 L/2 •Z2 L=1m •Z3

IZ=0,25cm4 IZ=0,25cm4 X
FY1 FY2 FY3
1°)
 96 24 − 96 24   12 6 − 12 6 
 24 − 24 4   6 −6 2 
[ ]
k = 500
1  8
− 96 − 24 96 − 24
[k ]2
= 500  4
− 12 − 6 12 − 6
   
 24 4 − 24 8   6 2 −6 4
TD-2A-S4-F412
 96 24 − 96 24 
 24 − 24 4 
[ ]
k 1 = 500
8
− 96 − 24 12
  96 6− 24−
 12 6
 24 4  −624 4 − 6
8 2 

 − 12 −6 12 −6 
[ ]
500 = k 2
 
 6 2 −6 4
 96 24 − 96 24 0 0
 24 8 − 24 4 0 0 

− 96 − 24 108 −18 −12 6 
[K ] = 500 
 24 4 − 18 12 − 6 2 
 0 0 −12 − 6 12 − 6
 
 0 0 6 2 − 6 4 
2°)
 FY 1   96 24 − 96 24 0  0 
0
0   24 8 − 24 4 0 θZ 1 
0
    
FY 2  − 96 − 24 108 − 18 − 12 6  0 
  = 500   
144  24 4 − 18 12 − 6 2 θZ 2 
 FY 3   0 0 − 12 − 6 12 − 6  0 
    
 0   0 0 6 2 − 6 4 θ Z 3 
3°)
 0  8 4 0  θ Z 1 
     θ Z 1 = θ Z 3 = −0.016 rd
144  = 500 4 12 2  θ Z 2 
θ Z 2 = 0.032 rd
 0  0 2 4  θ Z 3 
 FY 1 = 192   96 24 − 96 24 0 0  0 
 0   24 8 − 24 4 0 0  − 0.016
  
FY 2 = −144 − 96 − 24 108 − 18 − 12 6   0 
  = 500   
 144   24 4 − 18 12 − 6 2   0.032 
 FY 3 = −48   0 0 − 12 − 6 12 − 6  0 
    
 0   0 0 6 2 − 6 4  − 0.016
4°)
v ( x ) = [A ][D e ] avec [A ] = 1 − 3 x + 2 x x2 x3 
2 3
2x2 x3 3x2 2x3
 2 3
x − + 2 − 3 − + 2 
 L L L L L2 L L L 
[A ] = 
L 1 1 1 1
x = L = 1m −
2  2 8 2 8 
 0 
 0.032 
1 1 1 1 1
v( ) = −   = 0.006 m = 6 mm
2 2 8 2 8  0 
 
− 0.016 
TD-2A-S4-F412
PROBLÈME N°22
1ère partie
La poutre AB de longueur L et de section rectangulaire (b x h) repose sur deux appuis,

comme indiqué sur la figure ci-dessous. Une force P1 est appliquée au milieu de la
poutre.
On se propose de déterminer l’équation de la déformée de cette poutre en flexion par
la méthode de Maxwell-Mohr. (On négligera l’effet de l’effort tranchant devant celui
du moment fléchissant).
Pour cela, on applique une force verticale fictive P2 en un point D. La distance AD est
variable. Le paramètre correspondant est AD=x (0<x<L/2).
P1
P2
b
I A D B I
h
x 3
II L/2 II
L
1. Justifier le choix d’introduire cette force fictive afin de déterminer l’équation de la

déformée de la poutre.
2. En associant la force P1 à l’état #1 et la force P2 à l’état #2, écrire la relation entre
le vecteur des charges, le vecteur des déplacements et la matrice de flexibilité.
3. Calculer les termes de flexibilité nécessaires à la résolution du problème.
4. Déterminer l’équation de la déformée de la poutre.
5. En déduire la position et la valeur de la flèche maximale.
2ème partie
La poutre de la 1ère partie est désormais encastrée en A. Une force verticale P1 est
appliquée en son centre (point C). Un couple P2 porté par l’axe I est appliquée à son
extrémité (point B).
Données numériques : L = 500 mm E = 210 000 MPa

6
b= 20 mm h = 30 mm K = 10 N/m |P2| = 1000 Nm
TD-2A-S4-F412
P1
b P2
I A I
h
C B 3
II L/2 L/2 II
On utilisera la méthode de Castigliano. (On négligera l’effet de l’effort tranchant

devant celui du moment fléchissant).
1. Déterminer l’expression littérale de l’énergie de déformation stockée dans la poutre

ABC. Puis l’écrire sous la forme
wdef = α.P12 + β.P2 2 + γ.P1.P2
Dans la suite du problème, on utilisera les variables α, β, γ dans les
développements.
2. Calculer (forme littérale et application numérique) la flèche en B et C pour la même

poutre dans la configuration suivante
b P1 = 0 P2
I A I
h
B 3
II II
3. Calculer (forme littérale et application numérique) la flèche en B et la réaction en C

pour la même poutre dans la configuration suivante
P2
b
I A C
I
h
B 3
II L /2 L /2 II
4. Calculer (forme littérale et application numérique) la flèche en B et en C, ainsi que la

réaction en C pour la même poutre dans la configuration suivante (le ressort a une
raideur K).En dehors de tout chargement, le ressort n’exerce aucune force sur la
TD-2A-S4-F412
poutre)
P2
b
I A C
I
h
B 3
II L /2 L /2 II
Réponses: 1ère partie
2éme partie
RÈPONSES N°22
TD-2A-S4-F412
ANNEXES
TD-2A-S4-F412
ANNEXE A1. INTEGRALES DE MOHR
A1.1 DEFINITION:
Les intégrales de MOHR sont des intégrales du type:
L
I1 = ∫ y 1 ( x)y 2 ( x)dx
0
ou encore
L
I2 = ∫ 0
y 2 ( x)dx
Les résultats de ces intégrales sont donnés dans une table (voir pages suivantes) pour des
combinaisons de diverse fonctions y1(x) et y2(x) linéaires ou du second degré.
A1.2 UTILISATION DE LA TABLE:
1°) Intégrales de type MiMk.

Soit à calculer l'intégrale du type I1 du produit des deux fonctions y1 et y2 linéaires données
par leur représentation graphique ci-dessous. Les quantités Mi, Mk sont positives.
y
y 1(x)
Mi (y1) Mi
......... x
0 L
Mk (y2)
...........
y
y2(x) L 1
I1 = ∫ y1 y 2 dx = MiMkL
Mk
x 0 6
0 L
Vérification:
Mi
y1 ( x ) = − x + Mi
L
Mk
y2 ( x ) = x
L
TD-2A-S4-F412
L L Mi M L MM MM
I1 = ∫ y1 ( x ) y2 ( x )dx = ∫ ( − x + Mi ) k xdx = ∫ ( − i 2 k x 2 + i k x )dx
0 0 L L 0 L L
L
 MM MM  MM L MM L
I 1 = − i 2 k x 3 + − i k x 2  = − i k + i k
 3L 2L 0 3 2
MM L
I1 = i k
6
2°) Intégrales de type Mi2.

Soit à calculer l'intégrale du type I2 de la fonction linéaire y au carré donnée par sa
représentation graphique ci-dessous.
y
y(x) L 1
∫y
Mi
I2 = 2
dx = Mi2L
x 0 3
0 L
Vérification:
Mi
y1 ( x ) = x
L
L L Mi M2L
I2 = ∫0
y 2 dx = ∫ (
0 L
x ) 2 dx = i
3
3°) Cas hors table.
Exemple 1 :
y
0 L/2 x
-
-Mi Soit à calculer, à l'aide de la table, l'intégrale du type
y1(x) I1 du produit des deux fonctions y1 et y2 données par leur
représentation graphique ci-contre.
y
Mk y2(x)
+ x
0 L
TD-2A-S4-F412
L 3
Réponses: I1 = ∫ y1 y2 dx = − Mi Mk L
0 8
Exemple 2 :
y1
Soit à calculer, à l'aide de la table, l'intégrale du type
+ Mi I1 du produit des deux fonctions y1 et y2 données par leur
représentation graphique ci-contre.
y2
+ Mk
Mk/2
TD-2A-S4-F412
L 5
Réponses: I1 = ∫
0
y1 y2 dx =
12
Mi Mk L
TD-2A-S4-F412
A1.3 TABLE D’INTÉGRATION :
TD-2A-S4-F412
TD-2A-S4-F412
TD-2A-S4-F412
ANNEXE A2. MOMENTS FLÉCHISSANTS DANS QUELQUES

CAS SIMPLES :
MY
MY L 1N
1
+ 1Nm +
L x
L x
Z
+
G x MY
MY 1N
a L=a+b
1Nm L
x
- x -
-ab/L
-1
MY
MY 1 1N
L/2 L x
+ 1Nm
-
L x
-L/4
MY a/L
MY b
+ 1Nm L 1N
a - + x
x
a a+b=L
-b/L
MY 1/2
mY
+ 1Nm L 2
L /2 1N/m
L/2 - x + x
L
-1/2
MY MY
1 1N/m
+ 1Nm L
L/2 x
x -
a a+b=L
-L²/8
TD-2A-S4-F412

Elementsfini 150517232857 Lva1 App6891

Transféré par

Droits d'auteur :

Formats disponibles

Elementsfini 150517232857 Lva1 App6891

Transféré par

Informations du document

Titre original

Copyright

Formats disponibles

Partager ce document

Partager ou intégrer le document

Options de partage

Avez-vous trouvé ce document utile ?

Ce contenu est-il inapproprié ?

Droits d'auteur :

Formats disponibles

Elementsfini 150517232857 Lva1 App6891

Transféré par

Droits d'auteur :

Formats disponibles

SEMESTRE 4

1. MÉTHODE DE CLAPEYRON ............................................................................................................................... 3

C Les barres AB et BC (de même section droite A et de même

Après application du principe d’action mutuelle

Appliquons la méthode de Clapeyron : Wext = Wdéf

n2A 1°) Déterminez l'expression de l'énergie de déformation élastique en

Le déplacement associé à la charge P en C (•hC), est le raccourcissement de la barre •L :

1°) Déterminez l’expression de l’énergie de

La poutre ABC est soumise à de la torsion pure : M xCB = −C M xBA = −C

2°) Application numérique : Wdéf=12J, L=1m, d=40mm, G=80 GPa

24CL 24 × 802121× 1000

300C 02 450C 2 300C 0 450× 56 × C 0 C 0  300 450× 562 

L Déterminez l'expression littérale du

Réponses : PL3 PL PL3 1 + • h2 PL3 0.65 erreur 0.65%

Détermination du travail extérieur :

Le déplacement associé à la charge P en A (•vA), est le déplacement vertical de la section A :

Nous devons chercher le Torseur de Section en G de la barre AB :

Puis l’énergie de déformation due au moment fléchissant:

Application numérique : si ν = 0.3 et L/h = 10

L’énergie de déformation due à l’effort tranchant, ne représente que moins de 1 % de

Nous en concluons que l’énergie de déformation due à l’effort tranchant est en

Expression littérale du déplacement vertical de A

Appliquons la méthode de Clapeyron : Wext = Wdéf

La flèche due à l’effort tranchant ne représente que 0,65% du déplacement

Détermination du Torseur Externe: Voir module F112

Détermination du Torseur de Section: Voir module F112 (Annexes)

Détermination du travail extérieur :

Le déplacement associé au couple C en C (•YC), est la rotation suivant l’axe du couple de la

Détermination de l’énergie de déformation :

Appliquons la méthode de Clapeyron : Wext = Wdéf

Déterminez en ne tenant compte que du moment fléchissant

Détermination du Torseur Externe: Voir module F112

Détermination du Torseur de Section: Voir module F112 (Annexes)

Détermination du travail extérieur :

Le déplacement associé à la charge P en C (•vC), est le déplacement vertical de la section A :

Détermination de l’énergie de déformation :

Appliquons la méthode de Clapeyron : Wext = Wdéf

Réponses : • A = −1,12 mm • B = 1,94 mm

Considérons l’état initial E0, et utilisons la notation charges généralisées Pi et déplacements

Le déplacement généralisé D1 associé à la charge généralisée P1 est le déplacement linéique

Le déplacement généralisé D2 associé à la charge généralisée P2 est le déplacement linéique

La matrice carrée [f] est la matrice de flexibilité de la structure.

Chaque coefficient d’influence fij est calculé par le théorème de Maxwell-Mohr

[D] = [ f ][P ]  D1   f11

D1 = 4,345 .10−5 × 30.103 − 1,616 .10 −5 × 150.103 = −1,12mm

Les barres AB et BC (de même section droite A et de

Suggestion : Rajoutez en B une force horizontale

Application numérique : P1=10kN, L=2m, E=200GPa et A=1cm2. Calculez la valeur du

La méthode de Maxwell Mohr ou méthode de la charge

Sur le problème initial qui nous est proposé, nous ne pouvons

Nous allons donc créer un problème plus général, dont l’initial

Etat initial Etat

E0 • E0’{ P1=P ; P2=0 }

E0 • E0’ { P1=P ; P2=0 }

E0’ • P1×E1+P2×E2 45°

1°) Expression littérale de la matrice de flexibilité de la structure :

Chaque coefficient d’influence fij est calculé par le théorème de Maxwell-Mohr

2°) Expression littérale du déplacement vertical et du déplacement horizontal du

Dans l’état E0’ [D ]E = [ f ]E0 ' [P ]E0 '  D1  δ vB 