ECOLE NATIONALE POLYTECHNIQUE

DEPARTEMENT: GENIE CIVIL

CHAPITRE:III

CALCUL DES STRUCTURES PAR LA METHODE DE

FLEXIBILITE BASEE SUR LES DEFORMATIONS
CONSISTANTES

M. DEMIDEM
[email protected]
3.1 Introduction:

Dans cette méthode, un certain nombre de définitions

basées sur le calcul matriciel sont développées.
f ∆
 Le sens positif est: ,

 La relation déplacement-force est donnée par: ∆= f.P

f: étant la flexibilité , ∆: déplacement , P: force

 Sous forme matricielle on écrit: ∆ =[f] 𝑃

3.2 Définition d’un coefficient de rigidité:

Un coefficient de rigidité 𝐾𝑖𝑗 se définit comme étant la

force au point (𝑖) résultant d’un déplacement unitaire
imposé au point 𝑗 , sans déplacer les autres points.

Poutre (a)

Poutre (b)

Poutre (c)
 𝑃1 = 𝐾11 ∆1 + 𝐾12 ∆2 + 𝐾13 ∆3
𝑃2 = 𝐾21 ∆1 + 𝐾22 ∆2 + 𝐾23 ∆3
𝑃3 = 𝐾31 ∆1 + 𝐾32 ∆2 + 𝐾33 ∆3

On rappelle que :
𝜕𝑃1 𝜕𝑃1
𝐾11 = , 𝐾12 = … 𝑒𝑡𝑐 …
𝜕∆1 𝜕∆2
𝜕𝑃𝑖
𝐾𝑖𝑗 =
𝜕∆𝑗
𝝏𝑼
or, le 1er th de Castigliano s’écrit: 𝑷𝒊 =
𝝏∆𝒊
𝝏 𝝏𝑼 𝝏𝟐 𝑼
d’où 𝑲𝒊𝒋 = =
𝝏∆𝒋 𝝏∆𝒊 𝝏∆𝒊 𝝏∆𝒋
Enfin, on écrit les relations ci-dessus sous la forme
matricielle suivante :

𝐾11 𝐾12 ⋯ 𝐾1𝑛 ∆1 𝑃1

⋮ 𝐾22 … 𝐾2𝑛 ⋮ ⋮
=
⋮ … … … ⋮ ⋮
𝐾𝑛1 ⋯ ⋯ 𝐾𝑛𝑛 ∆𝑛 𝑃𝑛

𝐾 ∆ = 𝑃
3.3 Définition d’un coefficient de flexibilité:

Par analogie avec la définition du coefficient de rigidité 𝐾𝑖𝑗 ,

un coefficient de flexibilité 𝑓𝑖𝑗 est défini comme étant la
déflection au point (i) résultant d’une charge unitaire
appliquée au point (j). A partir des 3 poutres (a),(b),(c) vues
précédemment, on peut écrire les relations suivantes:

∆1 = 𝑓11 + 𝑓12 + 𝑓13

∆2 = 𝑓21 + 𝑓22 + 𝑓23
∆3 = 𝑓31 + 𝑓32 + 𝑓33

Comme les déplacements sont proportionnels aux forces:

∆1 = 𝑓11 𝑃1 + 𝑓12 𝑃2 + 𝑓13 𝑃3
∆2 = 𝑓21 𝑃1 + 𝑓22 𝑃2 + 𝑓23 𝑃3
∆3 = 𝑓31 𝑃1 + 𝑓32 𝑃2 + 𝑓33 𝑃3

1 =1

(1) (2) (3)

11 21 31

2 =1

12 22 32

3 =M=1

13 23 33
On note que :
𝜕∆1 𝜕∆1
𝑓11 = , 𝑓12 = … 𝑒𝑡𝑐 …
𝜕𝑃1 𝜕𝑃2
𝜕∆𝑖
𝑓𝑖𝑗 =
𝜕𝑃𝑗
𝝏𝑼
or le 2eme th de Castigliano s’écrit: ∆𝒊 =
𝝏𝑷𝒊
𝝏 𝝏𝑼 𝝏𝟐 𝑼
d’où 𝒇𝒊𝒋 = =
𝝏𝑷𝒋 𝝏𝑷𝒊 𝝏𝑷𝒊 𝝏𝑷𝒋
Enfin, on écrit les relations ci-dessus sous la
forme matricielle suivante :

𝑓11 𝑓12 ⋯ 𝑓1𝑛 𝑃1 ∆1

⋮ 𝑓22 … 𝑓2𝑛 ⋮ ⋮
=
⋮ … … … ⋮ ⋮
𝑓𝑛1 ⋯ ⋯ 𝑓𝑛𝑛 𝑃𝑛 ∆𝑛

𝑓 𝑃 = ∆
3.4 Relation entre l’énergie U , [f] et [K]:

Soit, à titre d’exemple, un ressort de raideur K:

K
∆

∆= 𝑓. 𝑃 , sous forme matricielle: ∆ = 𝑓 𝑃 = 𝑓 𝐾 ∆

1 1
U= 𝐾∆ , sous forme matricielle U= {∆}𝑇 𝐾 ∆
2
2 2
1 1 1
U= {𝑃} [𝑓] 𝐾 ∆ = {𝑃} 𝑓 𝐾 ∆ = {𝑃}𝑇 [𝑓]
𝑇 𝑇 𝑇 𝑃
2 2 2
[𝑓]𝑇 = [𝑓] : matrice est Auto-adjointe
1 1
U= {𝑃}𝑇 𝑓 𝑃 = {∆}𝑇 𝐾 ∆ , 𝑓 𝐾=𝐼
2 2
Quelques exemples:

 Calculer le coefficient 𝑓22 de la console suivante:

E A L N

(1) (2)

1 𝐿 𝑁2 𝑀2 𝜆𝑇 2 𝑚𝑡2
𝑈= න + + + 𝑑𝑥
2 0 𝐸𝐴 𝐸𝐼 𝐺𝐴 𝐺𝒥
𝟏 𝑳 𝑵𝟐 𝑵𝟐 𝑳
𝑼= ‫׬‬ 𝒅𝒙 =
𝟐 𝟎 𝑬𝑨 𝟐𝑬𝑨
𝝏𝟐 𝑼 𝝏𝟐 𝑼
𝒇𝟐𝟐 = = = (calcul à faire…)
𝝏𝑷𝟐 . 𝝏𝑷𝟐 𝝏𝑵. 𝝏𝑵
 Calculer le coefficient de flexibilité à mi-travée de la
poutre simple suivante:
P

0.5 L 0.5 L

(1) EI , L (2)

1 𝐿 𝑀2
𝑈= ‫׬‬0
𝑑𝑥 , M=M(z) = 0.5Px
2 𝐸𝐼
1 0.5𝐿 𝑀2 𝑃2 𝐿3
𝑈=2 ‫׬‬ 𝑑𝑥 =
2 0 𝐸𝐼 96𝐸𝐼
𝜕2 𝑈
𝑓= = (calcul à faire …)
𝜕𝑃2
 Pour la poutre simple suivante, calculer la matrice de
flexibilité reliant les ddl (𝜃𝑖 ) aux efforts appliqués.

𝜃1 𝑓11 𝑓12 𝑀1
=
𝜃2 𝑓21 𝑓22 𝑀2

1 2

(1) EI , L (2)

1 𝐿 𝑀2
𝑈= ‫׬‬ 𝑑𝑥
2 0 𝐸𝐼
 𝑀1 +𝑀2 𝑀1 +𝑀2
M = -𝑀1 + x , 𝑅1 = et 𝑅2 = − 𝑅1
𝐿 𝐿

2
1 𝐿 𝑀1 + 𝑀2
𝑈= න −𝑀1 + x 𝑑𝑥
2𝐸𝐼 0 𝐿
𝜕2 𝑈
𝑓11 = =…
𝜕𝑀1 𝜕𝑀1
𝜕2 𝑈
𝑓12 = =…
𝜕𝑀1 𝜕𝑀2
𝜕2 𝑈
𝑓22 = =…
𝜕𝑀2 𝜕𝑀2

𝑎 𝑏 𝑀1 𝜃1
[𝑓𝑖𝑗 ] , =
𝑐 𝑑 𝑀2 𝜃2
 Pour la poutre ci-dessous, calculer la matrice de
flexibilité reliant les ddl(𝑣, 𝜃2 ) aux efforts appliqués.

𝑓11 𝑓12 𝑃 𝑣
= 𝜃
𝑓21 𝑓22 𝑀0 2

(1) EI , L (2)

1 𝐿 𝑀2
𝑈= න 𝑑𝑥
2 0 𝐸𝐼
En appliquant le principe de superposition des charges,
les deux réactions de la poutre sont données par:
𝑃 𝑀0 𝑃 𝑀0
𝑅1 = + , 𝑅2 = -
2 𝐿 2 𝐿
Il est plus approprié d’évaluer le moment global M à
partir des deux extrémités de la poutre considérée.
1 𝐿/2 𝑃 𝑀0 2
𝑈= ‫׬‬ + x 𝑑𝑥 +
2𝐸𝐼 0 2 𝐿
2
1 𝐿/2 𝑃 𝑀0
න − x + 𝑀0 𝑑𝑥
2𝐸𝐼 0 2 𝐿
1 𝑃2 𝐿3 𝑀2 𝐿 𝑃𝑀𝐿2
𝑈= ( + + )
𝐸𝐼 96 6 16
 𝜕2 𝑈
𝑓11 = =…
𝜕𝑃2
𝜕2 𝑈
𝑓12 = =…
𝜕𝑃 𝜕𝑀0
𝜕2 𝑈
𝑓22 = =…
𝜕𝑀02

𝑎 𝑏 𝑃 𝑣
[𝑓𝑖𝑗 ] , = 𝜃
𝑐 𝑑 𝑀0 2
 Soit la poutre hyperstatique suivante. On demande
d’utiliser la méthode de flexibilité basée sur les
déformations consistantes pour calculer les réactions
inconnues. Les équations à utiliser sont:
∆ 𝑟 = ∆ 𝑟𝑙 + [𝑓]𝑟𝑟 𝑅 = 0
𝑃 𝑠 = 𝑃 𝑠𝑙 + 𝑃 𝑠𝑟 𝑅

EI , L , q

(1) (2)

∆ 𝑟 : déplts dans les appuis surabondants

∆ 𝑟𝑙 : déplts dûs aux charges 𝐿 : load
[𝑓]𝑟𝑟 : matrice de flexibilité
𝑅 : forces redondantes inconnues
𝑃 𝑠 : réactions aux supports
𝑃 𝑠𝑙 : réactions dues aux charges 𝐿
𝑃 𝑠𝑟 : réactions dues aux forces unitaires
redondantes

(1) E,I,L,q (2)

Str isostatique de base {∆ }

 𝑃𝑠𝑙1 = 𝑞𝐿 , 𝑃𝑠𝑙2 = 𝑞𝐿2 /2
𝑃𝑠𝑟1 = −1 , 𝑃𝑠𝑟2 = −𝐿
𝑓2 : coefficient d’influence = 𝐿3 /3𝐸𝐼
1 1 𝑞𝐿4
∆2 = ∆ 𝑟𝑙 = ‫׬‬ − 𝑞𝑥 2 −𝑥 𝑑𝑥 =
𝐸𝐼 2 8𝐸𝐼

2 1

2 f22

1 =1
 −𝑞𝐿4
∆ 𝑟𝑙 = : car ce déplacement est vers le bas
8𝐸𝐼
𝑓22 =𝐿3 /3𝐸𝐼
∆ 𝑟 = ∆ 𝑟𝑙 + [𝑓]𝑟𝑟 𝑅
−𝑞𝐿4
0= + (𝐿3 /3𝐸𝐼) 𝑅
8𝐸𝐼
⇒ R= +3qL/8

𝑃 𝑠= 𝑃 𝑠𝑙 + 𝑃 𝑠𝑟 𝑅
+𝑞𝐿 −1 +5𝑞𝐿/8
𝑃 𝑠= 2 + (3qL/8)=
+𝑞𝐿 /2 −𝐿 +𝑞𝐿2 /8
C’est les réactions au support (1)
 Soit la poutre hyperstatique suivante. On demande
d’utiliser la méthode de flexibilité basée sur les
déformations consistantes pour calculer les réactions
inconnues. Les équations à utiliser sont:

∆ 𝑟 = ∆ 𝑟𝑙 + [𝑓]𝑟𝑟 𝑅 = 0
𝑃 𝑠 = 𝑃 𝑠𝑙 + 𝑃 𝑠𝑟 𝑅
2P P

(1) EI,L (2) EI,L (3)

Structure isostatique de base (SIB)

H=2, le choix des inconnues hyperstatiques: 𝑅2 et 𝑅3 .
2P P

(1) EI,L (2) EI,L (3)

2 =2.5PL 2P P

1 =3P (SIB)

Le calcul des déplacements ∆ 𝑟𝑙 càd (∆2 et ∆3 ) dûs aux

charges effectives appliquées à la (SIB) donne:
19𝑃𝐿3
∆2 = − 3𝑃
24𝐸𝐼
∆ 𝑟𝑙 = 103𝑃𝐿3
: vers le bas (-) ; 𝑃 𝑠𝑙 =
∆3 = − 2.5𝑃𝐿
48𝐸𝐼
Les coefficients de flexibilité sont: 𝑓22 =𝐿3 /3EI , 𝑓32 =5𝐿3 /6EI

(1) (2) (3)

2 =-L =1

1 =-1 (SIB)

Les coefficients de flexibilité sont: 𝑓33 =8𝐿3 /3EI , 𝑓23 =5𝐿3 /6EI

(1) (2) (3)

2 =-2L =1

1 =-1 (SIB)
 𝐿3 5𝐿3
3𝐸𝐼 6𝐸𝐼 −1 −1
𝑓 𝑟𝑟 = , 𝑃 𝑠𝑟 =
5𝐿3 8𝐿3 −𝐿 − 2𝐿
6𝐸𝐼 3𝐸𝐼

∆ 𝑟= ∆ 𝑟𝑙 + [𝑓]𝑟𝑟 𝑅 = 0

19𝑃𝐿3 𝐿3 5𝐿3
0 − 𝑅2 0
24𝐸𝐼 3𝐸𝐼 6𝐸𝐼
∆ 𝑟= = + =
0 103𝑃𝐿3 5𝐿3 8𝐿3 𝑅3 0
−
48𝐸𝐼 6𝐸𝐼 3𝐸𝐼

La solution donne: 𝑅2 =93P/56 et 𝑅3 =16P/56

 𝑅1
Calcul des réactions restantes:
𝑀1 𝑠
𝑃 𝑠= 𝑃 𝑠𝑙 + 𝑃 𝑠𝑟 𝑅

𝑅1 3𝑃 −1 − 1 𝑅2 59𝑃/56
= + =
𝑀1 2.5𝑃𝐿 −𝐿 − 2𝐿 𝑅3 15𝑃𝐿/56
2P P

(1) EI,L (2) EI,L (3)

1 =15PL/56 2P P

1 =59P/56 2 =93P/56 3 =16P/56