SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

UNIVERSITE BADJI-MOKHTAR ANNABA

FACULTE DES SCIENCES DE L’INGENIEUR
DEPARTEMENT DE GENIE-CIVIL
3èmeLicence Génie CivilS6
SERIE N°1ENERGIE POTENTIELLE DE DEFORMATION
EXERCICE N° 1
Soit la poutre suivante (Figure 1)

1. Calculer l’énergie potentielle de déformation totale de la poutre

2. Comparer les énergies potentielles dues au moment fléchissant et aux efforts
tranchant.Quelle conclusion en tirez-vous ?
𝟔
AN :E= 2,1 𝟏𝟎𝟖 𝒌𝑵/𝒎𝟐 G= 0,8 𝟏𝟎𝟕 𝒌𝑵/𝒎𝟐 k = 𝟓 (coef de section réduite pour
une section rectangulaire)
EXERCICE N°2:
Soit la structure de la figure 2.
1 En appliquant le théorème de Castigliano, déterminer le déplacement horizontal en E δEH
2 Déterminer la rotation E en E en appliquant les intégrales de Mohr
3 Calculer le déplacement vertical δEVen E,par la méthode graphique de Véréchaguine
4 Calculer le déplacement horizontal HC en C en utilisant la multiplication des diagrammes
par les tableaux de Mohr

EXERCICE N° 3 :Soit la structure ci-dessous figure 3

1. Déterminer le déplacement horizontal en B HB

Figure 3

Pr Hacene Chaouche Abdelmadjid Page 1

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

Solution série N°1 énergie potentielle de déformation

EXERCICE N° 1

1. Calculer l’énergie potentielle de déformation totale de la poutre

𝒏 𝑳𝒊 𝑵𝟐𝒊 𝒏 𝑳𝒊 𝑴𝟐𝒊 𝒏 𝑳𝒊 𝒌𝑻𝟐𝒊 𝒏 𝑳𝒊 𝑴𝟐𝒕𝒊
𝐔𝐭𝐨𝐭 = 𝒊=𝟏 𝒐 (𝟐𝑬𝑨) dx + 𝒊=𝟏 𝒐 (𝟐𝑬𝑰) dx + 𝒊=𝟏 𝒐 (𝟐𝑮𝑨) dx + 𝒊=𝟏 𝒐 (𝟐𝑮𝑰) dx
𝒊 𝒊 𝒊 𝑷𝒊

 Calcul des réactions d’appui

𝑭𝒙 = 𝟎 𝑯𝑨 = 𝑷 =𝟓 𝒌𝑵
𝑴
𝑴𝑩 = 𝟎𝑹𝑨 𝟔 + 𝑴 = 𝟎𝑹𝑨 = − = −𝟏, 𝟔𝟕 𝒌𝑵
𝟔
𝑭𝒚 = 𝟎 𝑹𝑨 + 𝑹𝑩 =0𝑹𝑩 = −𝑹𝑨 = 𝟏, 𝟔𝟕 𝒌𝑵
Tronçon AB : 0≤ 𝐱 ≤ 𝟏 𝐦
𝑵𝑨𝑪 = - 𝑯𝑨 = −𝟓 𝒌𝑵𝑵𝟐𝑨𝑪 = 𝟐𝟓 𝒌𝑵𝟐
𝑴𝑨𝑪 = 𝑹𝑨 𝒙 = −𝟏, 𝟔𝟕𝒙  𝑴𝑨𝑪 𝟐 = 2,78 𝒙𝟐
𝑻𝑨𝑪 = 𝑹𝑨 = −𝟏, 𝟔𝟕 𝒌𝑵  𝑻𝑨𝑪 𝟐 = 𝟐, 𝟕𝟖 𝒌𝑵𝟐
Tronçon CB : 0≤ 𝐱 ≤ 𝟓 𝐦
𝑵𝑩𝑪 = - P= −𝟓 𝒌𝑵𝑵𝟐𝑩𝑪 = 𝟐𝟓 𝒌𝑵𝟐
𝑴𝑪𝑩 = 𝑹𝑩 𝒙 = 𝟏, 𝟔𝟕𝒙𝑴𝑨𝑪 𝟐 = 2,78 𝒙𝟐
𝑻𝑪𝑩 = −𝑹𝑩 = −𝟏, 𝟔𝟕𝒌𝑵  𝑻𝑨𝑪 𝟐 = 𝟐, 𝟕𝟖 𝒌𝑵𝟐
 Calcul des caractéristiques géométriques
𝒃𝒉𝟑 𝟎,𝟏.𝟎,𝟐𝟑
A= 2 𝟏𝟎−𝟐 𝒎𝟐 ; Iz= 𝟏𝟐 = 𝟏𝟐 =6,667 𝟏𝟎−𝟓 𝒎𝟒
E= 2,1 𝟏𝟎𝟖 𝒌𝑵/𝒎𝟐 et G= 0,8 𝟏𝟎𝟕 𝒌𝑵/𝒎𝟐
𝒏 𝑳𝒊 𝑵𝟐𝒊 𝒏 𝑳𝒊 𝑴𝟐𝒊 𝒏 𝑳𝒊 𝒌𝑻𝟐𝒊
𝐔𝐭𝐨𝐭 = 𝒊=𝟏 𝒐 (𝟐𝑬𝑨) 𝐝𝐱 + 𝒊=𝟏 𝒐 (𝟐𝑬𝑰) dx + 𝒊=𝟏 𝒐 (𝟐𝑮𝑨) dx
𝒊 𝒊 𝒊
𝒏 𝑳𝒊
𝑵𝟐𝒊 𝟏
𝟐𝟓 𝟓
𝟐𝟓
𝐔𝐭𝐨𝐭𝐍 = 𝒅𝒙 = 𝒅𝒙 + 𝒅𝒙 = 𝟏, 𝟕𝟖𝟓 𝟏𝟎−𝟓 𝒌𝑵 . 𝒎
(𝟐𝑬𝑨)𝒊 𝟐𝑬𝑨 𝟐𝑬𝑨
𝒊=𝟏 𝒐 𝒐 𝒐

𝒏 𝑳𝒊
𝑴𝟐𝒊 𝟏
𝟐, 𝟕𝟖 𝒙𝟐 𝟓
𝟐, 𝟕𝟖 𝒙𝟐
𝐔𝐭𝐨𝐭𝐌 = 𝒅𝒙 = 𝒅𝒙 + 𝒅𝒙
(𝟐𝑬𝑰)𝒊 𝟐𝑬𝑰 𝟐𝑬𝑰
𝒊=𝟏 𝒐 𝒐 𝒐
𝟐,𝟕𝟖 𝟏 𝟏𝟐𝟓
𝐔𝐭𝐨𝐭𝐌 = + = 𝟒, 𝟏𝟕 𝟏𝟎−𝟑kNm
𝟐𝑬𝑰 𝟑 𝟑
𝒏 𝑳𝒊 𝒌𝑻𝟐𝒊
𝐔𝐭𝐨𝐭𝐓 = + 𝒊=𝟏 𝒐 (𝟐𝑮𝑨) dx
𝒊
𝟏 𝟔 𝟐,𝟕𝟖 𝟓 𝟔 𝟐,𝟕𝟖
𝐔𝐭𝐨𝐭𝐓 = 𝒐 𝟓 𝟐𝑮𝑨
𝒅𝒙 + 𝒐 𝟓 𝟐𝑮𝑨
𝒅𝒙 = +𝟎, 𝟎𝟎𝟎𝟎𝟏𝟎𝟒 + 𝟎, 𝟎𝟎𝟎𝟎𝟓𝟐𝟏𝟐𝟓 = 𝟎, 𝟎𝟎𝟎𝟎𝟔𝟐𝟓 kNm

𝑳 𝒊 𝑴𝟐 𝑳𝒊 𝒌𝑻𝟐 𝑳𝒊 𝑵𝟐𝒊
𝐔𝐭𝐨𝐭 = + 𝒏 𝒊
𝒊=𝟏 𝒐 (𝟐𝑬𝑰) dx +
𝒏 𝒊
𝒊=𝟏 𝒐 (𝟐𝑮𝑨) 𝐝𝐱 + 𝒏
𝒊=𝟏 𝒐 𝒅𝒙 =
𝒊 𝒊 (𝟐𝑬𝑨)𝒊

Pr Hacene Chaouche Abdelmadjid Page 2

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

𝐔𝐭𝐨𝐭 = 𝟏. 𝟕𝟖𝟓 𝟏𝟎−𝟓 + 𝟒. 𝟏𝟕 𝟏𝟎−𝟑 + 𝟔, 𝟐𝟓 𝟏𝟎−𝟓 =0,00425kNm

On peut conclure que les énergies potentielles de déformation dues à l’effort tranchant et
l’effort normal sont st négligeables
par rapport à celle du moment fléchissant.
EXERCICE N°2 :

1. Détermination du déplacement horizontal en E δEH en appliquant le théorème de

Castigliano.
Le théorème de Castigliano s’énonce comme suit : La dérivée partielle de l’énergie
potentielle de déformation d’une structure ou d’un corps solide par rapport à une
force est égale au déplacement du point d’application de cette force suivant sa ligne
d’action.
On néglige les énergies potentielles dues aux efforts tranchant et aux efforts normaux

𝑼 𝒏 𝑳𝒊 𝑴 𝒊  𝑴 𝒊
𝛅 𝐄𝐇 = = 𝒊=𝟏 𝟎 (𝑬𝑰) 𝑷 dx
𝑷 𝒊

Ce portique présente quatre tronçons, en l’occurrence, les tronçons AB, BC, CD et DE.
TronçonED :0≤ 𝐱 ≤ 𝐋
𝑴𝑬𝑫
MED= - Px =-x
𝑷
Tronçon DC : 0≤ 𝐱 ≤ 𝐋
𝑴𝑫𝑪
MDC= - PL  =- L
𝑷
Tronçon CB : 0≤ 𝐱 ≤ 𝟐𝐋
𝑴𝑫𝑪
MCB= - P(L – x)  =-
𝑷
(L – x)
Tronçon BA : 0≤ 𝐱 ≤ 𝟐𝐋
𝒒𝒙𝟐 𝑴𝑩𝑨
MBA= - PL −  = -L
𝟐 𝑷
𝟏 𝑳 𝟏 𝑳 𝟏 𝟐𝑳 𝟏 𝟐𝑳 𝒙𝟐
𝛅 𝐄𝐇 = 𝐏𝐱 𝟐 dx +𝑬𝑰 ( 𝐏𝐋𝟐 ) dx +𝑬𝑰 [ 𝐏(𝐱 − 𝐋)𝟐 dx+𝟐𝑬𝑰 𝐏𝐋𝟐 + 𝐏 dx
𝟐𝑬𝑰 𝟎 𝟎 𝟎 𝟎 𝟐

Pr Hacene Chaouche Abdelmadjid Page 3

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

𝑷𝑳𝟑 𝑷𝑳𝟑 𝟐𝑷𝑳𝟑 𝟑 𝟕𝑷𝑳𝟑

𝛅 𝐄𝐇 = + + + 𝟓𝑷𝑳 = (Même sens que la charge)
𝟔𝑬𝑰 𝑬𝑰 𝟑𝑬𝑰 𝟑𝑬𝑰 𝟐𝑬𝑰

2. Détermination de la rotation E en E en appliquant les intégrales de Mohr

On applique d’abord un moment fictif  en E puis on calcule la rotation E en dérivant
l’énergie potentielle totale de la structure par rapport à ce moment fictif, et faisant =0
dans les intégrales de Mohr, on obtient alors l’expression permettant de calculer la
rotation E.
L’énergie potentielle de déformation est donnée par :
𝒏 𝑳𝒊
𝑴𝟐𝒊
𝑼=
𝒐 (𝟐𝑬𝑰)𝒊 𝒊=𝟏
𝑼 𝒏 𝑳𝒊 𝑴𝒊 𝑴𝒊
E= = 𝒊=𝟏 𝟎 (𝑬𝑰) 
 𝒊

TronçonED :0≤ 𝐱 ≤ 𝐋
𝑴𝑬𝑫
MED= - Px- =-1

Tronçon DC : 0≤ 𝐱 ≤ 𝐋
𝑴𝑫𝑪
MDC= - PL -  =- 1

Tronçon CB : 0≤ 𝐱 ≤ 𝟐𝐋
𝑴𝑫𝑪
MCB= - P(L –x)- =- 1

Tronçon BA : 0≤ 𝐱 ≤ 𝟐𝐋
𝒒𝒙𝟐 𝑴𝑩𝑨
MBA=-PL - + =1
𝟐 
𝐔 𝟏 𝐋 𝟏 𝐋 𝟏 𝟐𝐋 𝟏 𝟐𝐋 𝐪𝐱 𝟐
𝑬 =  (=𝟎) = 𝟐𝐄𝐈 𝟎
(−𝐏𝐱)(−𝟏) dx + 𝟎 (− 𝐏𝐋)(−𝟏) dx + 𝟎 [ −𝐏 𝐋 − 𝐗 (−𝟏)
𝐄𝐈 𝐄𝐈
dx+
𝟐𝐄𝐈 𝟎
= (−𝐏𝐋 −
𝟐
) dx
𝟏 𝐋 𝟏 𝐋 𝟏 𝟐𝐋 𝟏 𝟐𝐋 𝐏𝐱 𝟐
𝑬 = 𝟐𝐄𝐈 𝟎
(𝐏𝐱)dx + 𝟎 𝐏𝐋 dx + 𝟎 [ 𝐏 𝐋 − 𝐱 dx+
𝐄𝐈 𝐄𝐈 𝟐𝐄𝐈 𝟎
(−𝐏𝐋 − ) dx
𝟐𝐋

𝟓
E= − 𝟏𝟐𝐄𝐈 𝐏𝐋𝟐 (rotation dans le sens inverse de )
3. Calcul du déplacement vertical δEV en E, par la méthode graphique de
Véréchaguine
Si la barre est à section constante (EI= constante)

Pr Hacene Chaouche Abdelmadjid Page 4

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

Supposons que l’on doive intégrer le produit deux fonctions f1(x) et f2(x),
respectivement non linéaire et linéaire.
I= 𝐟𝟏 (𝐱) 𝐟𝟐 (x) dx = 𝐟𝟏 (𝐱) (ax +b) dx = a 𝒙𝐟𝟏 𝐱 𝐝𝐱 +b 𝐟𝟏 (𝐱) dx = a xG1 A1 +b A1
Donc : I = A1(axG1+b)= A1 y2
A1 : Aire de la fonction non linéaire
y2= ordonnée mesuré sur le diagramme de la fonction linéaire au niveau du centre de gravité
du diagramme de la fonction non linéaire
Ce résultat peut s’énoncer de la manière suivante :
L’intégrale du produit de deux fonctions, dont l’une est linéaire est égal au produit de
l’aire de l’épure non linéaire par l’ordonnée de son centre de gravité mesurée sur
l’épure linéaire.

 Donc on a besoin d’un diagramme unitaire et d’un diagramme de charge, or dans

le but de simplifier les calculs et en vertu du principe de superposition, on trace
autant de diagramme de charge dus au moment fléchissant que de charges
appliquée sur la structure.
 Dans ce cas précis on a deux charges appliquées sur le portique, donc il est
impératif de tracer deux diagrammes du moment fléchissant dus au chargement
appliqués, en l’occurrence, les charges P et q

Pr Hacene Chaouche Abdelmadjid Page 5

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

Les diagrammes :

δEV=δEVP+δEvq(principe de superposition la somme des déplacements dus à la charge P et

la charge répartie q)
𝟏 𝟏 𝑷𝑳𝟑
δEVP=− 𝑬𝑰 (PL L) (𝟐L) + 0 +0= − 𝟐𝑬𝑰 ( le produit du diagramme de P Mp par le diagramme
unitaire m)
𝟏
δEvq=𝟐𝑬𝑰 𝑨𝒑𝒂𝒓𝒂𝒃𝒐𝒍𝒆 𝒚𝒍𝒊𝒏é𝒂𝒊𝒓𝒆 ( le produit du diagramme de q Mq par le diagramme unitaire m)

𝟏 𝟏 𝑷𝑳 𝟏 𝟏 𝑷𝑳 𝑳 𝟏 𝟑𝑷𝑳 𝟑𝑳 𝟕
𝛅𝐄𝐯𝐪 = 𝟐𝑬𝑰 𝑳 𝑳 -𝟐𝑬𝑰 ( 𝟐 𝑳) 𝟐 + (𝟑 𝑳) = − 𝟐𝟒 𝑬𝑰 𝑷𝑳𝟑
𝟑 𝟐 𝟒 𝟐 𝟒
𝑷𝑳𝟑 𝟕 𝟏𝟗𝑷𝑳𝟑
δEV=(-𝟐𝑬𝑰 - 𝟐𝟒 𝑬𝑰 𝑷𝑳𝟑 )= − (dirigé vers le bas) (sens inverse de la charge unitaire)
𝟐𝟒𝑬𝑰

4. Calculedu déplacement horizontal HC en C en utilisant la multiplication des

diagrammes par les tableaux de Mohr

Pr Hacene Chaouche Abdelmadjid Page 6

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

 Il s’agit en fait de trouver graphiquement le résultat des intégrales de Mohr.

Les diagrammes de charge étant déjà tracés pour le compte de la 3ème question
Il ne manque que le diagramme unitaire du à une charge appliquée horizontalement en C

δCH=δCHP+δCHq(principe de superposition la somme des déplacements dus à la charge P et

la charge répartie q)
𝑳 𝟏 𝑳 𝟏 𝟐𝑳 𝟖𝑷𝑳𝟑
δCHP(− 𝑬𝑰 𝟔 . 𝑷𝑳. 𝑳)+𝑬𝑰 𝟔 . 𝑷𝑳(𝑳 + 𝟐. 𝟐𝑳)+𝟐𝑬𝑰 . 𝟐𝑳. 𝑷𝑳= 𝟑𝑬𝑰
𝟐𝑳 𝟏 𝟒𝑷𝑳𝟑
δCHq=𝟐𝑬𝑰 𝟑.2PL. 2L= 𝟑𝑬𝑰
𝟒𝑷𝑳𝟑
δCH= (dirigé vers la droite)
𝑬𝑰

EXERCICE N° 3 :
.1. Déterminationdu déplacement horizontal en B HB
𝑵𝒊𝒋 𝒏𝑰𝑱 𝑳𝒊𝒋
HB=
𝑬𝑨𝒊𝒋

𝐍𝐢𝐣 ∶ 𝐞𝐟𝐟𝐨𝐫𝐭 𝐧𝐨𝐫𝐦𝐚𝐥 𝐝𝐚𝐧𝐬 𝐛𝐚𝐫𝐫𝐞 𝐢𝐣 𝐝𝐮 𝐚𝐮 𝐜𝐡𝐚𝐫𝐠𝐞𝐦𝐞𝐧𝐭 𝐚𝐩𝐩𝐥𝐢𝐪𝐮é

𝒏𝐢𝐣 : 𝐞𝐟𝐟𝐨𝐫𝐭 𝐧𝐨𝐫𝐦𝐚𝐥 𝐝𝐚𝐧𝐬 𝐛𝐚𝐫𝐫𝐞 𝐢𝐣 𝐝𝐮 à 𝐮𝐧𝐞 𝐜𝐡𝐚𝐫𝐠𝐞 𝐮𝐧𝐢𝐭𝐚𝐢𝐫𝐞 𝐚𝐩𝐩𝐢𝐪𝐮é 𝐚𝐮 𝐩𝐨𝐢𝐧𝐭 𝐞𝐭 𝐬𝐮𝐢𝐯𝐚𝐧𝐭 𝐥𝐚 𝐝𝐢𝐫𝐞𝐜𝐭𝐢𝐨
𝐋𝐢𝐣 : 𝐥𝐨𝐧𝐠𝐮𝐞𝐮𝐫 𝐝𝐞 𝐥𝐚 𝐛𝐚𝐫𝐫𝐞 𝐢𝐣
𝐄𝐀 𝐢𝐣 : 𝐫𝐢𝐠𝐢𝐝𝐢𝐭é à 𝐥′ 𝐞𝐟𝐟𝐨𝐫𝐭 𝐧𝐨𝐫𝐦𝐚𝐥𝐞 𝐝𝐞 𝐥𝐚 𝐛𝐚𝐫𝐫𝐞 𝐢𝐣

 Calcul des réactions

 Calcul des réactions du treillis sous chargement
𝐅𝐱 = −𝐇𝐀 + 𝟐𝐏 = 𝟎 𝐇𝐀 = 𝟐𝐏
𝐌𝐁 = 𝐑 𝐀 . 𝟑𝐋 – P. 2L+ 2P. L= 0 𝐑 𝐀 = 𝟎

Pr Hacene Chaouche Abdelmadjid Page 7

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

𝐅𝐲 = 𝐑 𝐀 + 𝐑 𝐁 = P  𝐑 𝐁 = P

 Calcul des réactions du treillis sous charge unitaire

𝐅𝐱 = 𝐇𝐀 − 𝟏 = 𝟎 𝐇𝐀 = 𝟏
𝟏
𝐌𝐁 = 𝐑 𝐀 . 𝟑𝐋 +1 . L= 0 𝐑 𝐀 = − 𝟑
𝟏
𝐅𝐲 = 𝐑 𝐀 + 𝐑 𝐁 = 0  𝐑 𝐁 =- 𝐑 𝐀 =𝟑

 Les efforts dans les barres sont calculés par la méthode graphique de Cremona
(Figures ci-dessus)

Barres 𝑵𝒊𝒋 𝒏𝑰𝑱 𝑳𝒊𝒋 𝑬𝑨𝒊𝒋 𝑵𝒊𝒋 𝒏𝒊𝒋 𝑳𝒊𝒋

𝑬𝑨𝒊𝒋
AB= C1 0 0,333 L 2EA 0
AC=1a 2P 1 L EA 𝟐𝑷𝑳
𝑬𝑨
BC=21 1,425 P 0 1,414 L EA 0
BD=d2 -3,175 P -1,08 1,054 L EA 𝟑, 𝟔𝟏𝑷𝑳
𝑬𝑨
DE =d3 -3,175 P -1,08 2,108 L EA 𝟕, 𝟐𝟐𝑷𝑳
𝑬𝑨
DC= 32 0 0 0,666 L 2EA 0
CE=e3 3P 1 2L EA 𝟔𝑷𝑳
𝑬𝑨
𝑵𝒊𝒋 𝒏𝒊𝒋 𝑳𝒊𝒋 𝟏𝟖, 𝟖𝟑𝑷𝑳
𝑯𝑩 =
𝑬𝑨
𝑬𝑨𝒊𝒋 Vers la droite
(même sens
que la charge
unitaire)

Pr Hacene Chaouche Abdelmadjid Page 8

 SERIE N° 1 ENERGIE POTENTIELLE DE DEFORMATION

Pr Hacene Chaouche Abdelmadjid Page 9