00 MDS 1920 PDF

Centre des Métiers d'Education et de Formation
Settat
Cycle Préparatoire d’Agrégation

Sciences Industrielles de l’Ingénieur
Mécanique des structures
M'hammed EL ALLAMI
[email protected]
Année universitaire 2019 – 2020

Table des matières
Table des matières ................................................................................................ 1

Chapitre 1 Généralités et notions de base. ........................................................... 6
1.1 Buts de la mécanique des structures ................................................................................. 6
1.2 Hypothèses ....................................................................................................................... 7
1.2.1 Le matériau ................................................................................................................... 7
1.2.2 La disposition de la matière ......................................................................................... 7
1.2.3 Les actions mécaniques extérieures ............................................................................. 8
1.2.4 Les déformations .......................................................................................................... 8
1.3 Principe fondamental de la statique ................................................................................. 9
1.4 Torseur des efforts de cohésion ...................................................................................... 10
1.4.1 Définition ................................................................................................................... 10
1.4.2 Composantes du torseur de cohésion ......................................................................... 11
1.5 Tenseur des contraintes .................................................................................................. 12
1.5.1 Vecteur contrainte en un point M ............................................................................... 12
1.5.2 Contrainte normale/ contrainte tangentielle : ............................................................. 13
1.5.3 Tenseur des contraintes .............................................................................................. 13
1.5.4 Torseur de cohésion : écriture locale.......................................................................... 15
Chapitre 2 Traction / compression des poutres. ................................................. 16
2.1 Définition ....................................................................................................................... 16
2.2 Essai d'traction ............................................................................................................... 16
2.3 Caractéristiques mécaniques d'un matériau.................................................................... 17
2.4 Contraintes et déformations ........................................................................................... 19
2.5 Déformation d'un tronçon infiniment petit ..................................................................... 19
2.6 Torseur des efforts intérieurs.......................................................................................... 20
2.7 Condition de résistance .................................................................................................. 20
2.8 Concentration de contraintes .......................................................................................... 21
2.9 Triés ................................................................................................................................ 22
2.9.1 Efforts sur la structure ................................................................................................ 22
M. El Allami Mécanique des structures
1
2.9.2 Efforts internes: effort normal dans les barres ........................................................... 23
2.9.3 Barre à effort normal nul N=0 .................................................................................... 23
2.9.4 Déformations des nœuds ............................................................................................ 24
2.9.5 Nœud articulé et nœud rigide ..................................................................................... 24
2.9.6 Hyperstaticité ............................................................................................................. 24
2.9.7 Résolution des triés isostatiques ................................................................................. 26
2.9.8 Résolution des triés hyperstatiques ............................................................................ 28
2.10 Mesures par jauges de contraintes .................................................................................. 28
Chapitre 3 Cisaillement des poutres. .................................................................. 30
3.1 Définition ....................................................................................................................... 30
3.2 Essai de cisaillement ...................................................................................................... 31
3.3 Caractéristiques mécaniques d'un matériau.................................................................... 32
3.3.1 Contrainte tangentielle limite élastique 𝛕𝒆 ou Rge .................................................... 32
3.3.2 Contrainte tangentielle de rupture .............................................................................. 32
3.4 Contrainte tangentielle de cisaillement .......................................................................... 32
3.5 Relation entre contrainte et déformation ........................................................................ 33
3.6 Condition de résistance .................................................................................................. 33
Chapitre 4 Torsion des poutres cylindriques...................................................... 34
4.1 Hypothèse supplémentaire ............................................................................................. 34
4.2 Définition ....................................................................................................................... 34
4.3 Essai de torsion............................................................................................................... 35
4.4 Déformations élastiques : distorsion .............................................................................. 36
4.5 Contraintes et déformations de torsion........................................................................... 36
4.6 Conditions de résistance ................................................................................................. 38
4.7 Condition de rigidité ou de déformation ........................................................................ 38
Chapitre 5 Flexion plane simple des poutres. .................................................... 41
5.1 Définition ....................................................................................................................... 41
5.2 Relation effort tranchant/moment fléchissant ................................................................ 41
5.3 Essai de flexion (domaine élastique) .............................................................................. 42
5.4 Déformation longitudinale et Contrainte normale ......................................................... 44
5.5 Déformations transversales : Équation de la déformée .................................................. 46
5.6 Contraintes tangentielles ................................................................................................ 46

2
5.7 Conditions de résistance ................................................................................................. 49
5.8 Condition de rigidité....................................................................................................... 50
5.10 Structure hyperstatique ................................................................................................... 50
Chapitre 6 Sollicitations combinées des structures. ........................................... 53
6.1 Repère de Frenet............................................................................................................. 53
6.2 Equilibre local ................................................................................................................ 54
6.3 Technique des fonctions de singularité .......................................................................... 55
6.4 Formules de Bresse ........................................................................................................ 56
6.4.1 Utilité.......................................................................................................................... 56
6.4.2 Torseur de petits déplacements .................................................................................. 56
6.4.3 Torseur de petites déformations ................................................................................. 57
6.4.4 Lois de comportement ................................................................................................ 57
6.4.5 Formules de Bresse .................................................................................................... 57
6.5 Théorèmes énergétiques ................................................................................................. 58
6.5.1 Travail des efforts extérieurs : première formule de Clapeyron ................................ 58
6.5.2 Énergie de déformation : forme intégrale .................................................................. 60
6.5.3 Egalité de Clapeyron .................................................................................................. 61
6.5.4 Théorèmes de réciprocité de Betti- Rayleigh et de Maxwell- Betti ........................... 61
6.5.5 Travail des efforts extérieurs : seconde formule de Clapeyron .................................. 63
6.5.6 Théorème de Castigliano ............................................................................................ 64
6.5.7 Théorème de Ménabréa : Pour les structures hyperstatiques ..................................... 66
Chapitre 7 Flambage .......................................................................................... 68
7.1 Essai de flambage ........................................................................................................... 68
7.2 Charge critique d’Euler .................................................................................................. 69
7.3 Cas possibles pour la poutre ........................................................................................... 71
7.4 Influence des liaisons aux appuis ................................................................................... 71
7.5 Calcul de résistance (Méthode d’Euler) ......................................................................... 72
7.5.1 Élancement : ............................................................................................................... 72
7.5.2 Coefficient de sécurité :.............................................................................................. 72
7.5.3 Élancement critique 𝛌𝐜............................................................................................... 72
7.5.4 Condition de résistance au flambage :........................................................................ 73
7.5.5 Critères de résistance .................................................................................................. 73

3
Chapitre 8 Vibration de systèmes continus. ....................................................... 74
8.1 Vibrations longitudinales ............................................................................................... 75
8.1.1 Vibrations libres ......................................................................................................... 75
8.1.2 Fréquences propres et modes propres ........................................................................ 76
8.1.3 Exemples .................................................................................................................... 77
8.1.4 Orthogonalité des modes ............................................................................................ 78
8.1.5 Normalisation ............................................................................................................. 79
8.1.6 Exemple...................................................................................................................... 79
8.1.7 Vibrations forcées ...................................................................................................... 80
8.2 Vibrations de torsion ...................................................................................................... 81
8.3 Vibrations transversales de flexion ................................................................................ 81
8.3.1 Fréquences et modes propres ..................................................................................... 82
8.3.2 Exemples .................................................................................................................... 83
8.3.3 Orthogonalité des modes ............................................................................................ 84
8.3.4 Normalisation ............................................................................................................. 85
8.3.5 Exemples .................................................................................................................... 85
8.3.6 Vibrations forcées ...................................................................................................... 86
8.4 Méthodes approchées ..................................................................................................... 87
8.4.1 Energie cinétique et énergie potentielle ..................................................................... 87
8.4.2 Méthode de Rayleigh : Approximation de la fréquence fondamentale ...................... 87
8.4.3 Méthode de Rayleigh- Ritz : Approximation des premières fréquences ................... 90
Chapitre 9 Calcul de structures par la méthode des éléments finis ................... 92
9.1 Les principes .................................................................................................................. 92
9.1.1 Définitions .................................................................................................................. 92
9.1.2 Famille d’éléments ..................................................................................................... 93
9.1.3 Exemples .................................................................................................................... 93
9.1.4 Approximation des déplacements .............................................................................. 94
9.1.5 Déformations .............................................................................................................. 94
9.1.6 Contraintes ................................................................................................................. 94
9.2 Résolution d’un problème de statique 1D ...................................................................... 95
9.2.1 Ecriture matricielle en utilisant l’énergie de déformation .......................................... 95
9.2.2 Analyse globale et assemblage ................................................................................... 95
Résolution............................................................................................................................... 96

4
9.2.3 Exemple1 : Analyse statique d'un système constitué de ressorts ............................... 97
9.2.4 Elément de barre en traction- compression ................................................................ 98
9.2.5 Elément de poutre en flexion ..................................................................................... 99
9.2.6 Exemples .................................................................................................................. 101
9.3 Résolution d’un problème 2D ...................................................................................... 101
9.3.1 Ecriture matricielle par utilisation du principe des travaux virtuels PTV ................ 101
9.3.2 Différents types d’éléments en 2D ........................................................................... 103
9.3.3 Résolution d’un problème plan ................................................................................ 103
9.3.4 Résolution................................................................................................................. 106
9.3.5 Exemple 2 ................................................................................................................. 106
Annexe A Jauges de déformations ................................................................... 107

5
Chapitre 1
Généralités et notions de base.
L’objectif de ce premier chapitre est de mettre en place et définir les notions de base permettant
d’aborder les chapitres suivants traitant de la mécanique des structures, plus communément
appelée Résistance des Matériaux.
1.1 Buts de la mécanique des structures

Établir, en vue de construire dans les meilleures conditions, les relations qui existent entre les
dimensions, les contraintes et les déformations d’un solide sous l’action de :
 forces extérieures qui lui sont appliquées,
 et compte tenu des propriétés des matériaux dont il est constitué,
La résistance des matériaux est utilisée pour 4 raisons :
Calcul de résistance
Il faut déterminer les dimensions d’une pièce pour :
 qu’elle résiste aux diverses forces externes
 et que les sollicitations internes ne dépassent pas certaines limites (en tout point).
Calcul de rigidité
Les déformations doivent être inférieures à une déformation limite que l’on s’est fixée.
Calcul de vérification
Certaines dimensions peuvent être imposées par le montage, il faudra donc vérifier si ces
dimensions conviennent :
 soit en fonctions des sollicitations
 Soit en fonction des déformations.
Choix des matériaux

Dans certains problèmes ou les dimensions ou la masse ou les déformations sont limitées
(industrie aéronautique par exemples), c’est la nature du matériau qui sera l’inconnue. Il faut
choisir le matériau ou le vérifier d’après les charges appliquées.

6
1.2 Hypothèses
Sont de 4 types :
1.2.1 Le matériau
Le matériau est :
 Homogène : il a les mêmes propriétés en tout point
 isotrope : en un point donné, il a les mêmes propriétés dans toutes les directions.
Le comportement du matériau est :
 linéaire : les relations entre les contraintes et les déformations sont linéaires
 Élastique : le solide reprend sa forme initiale si les forces appliquées sont supprimées.
Figure 1 –Différents comportement du matériau
1.2.2 La disposition de la matière

La RDM étudie des pièces dont les formes sont relativement simples. Ces pièces sont désignées
sous le terme de « poutres ».
(S ) G
FIG .1 (C ) L i g n e m o y e n n e
Figure 2 –Poutre

7
Poutre
On appelle poutre (voir figure) un solide engendré par une surface plane (S) dont le centre de
surface G décrit une courbe plane (C) appelée ligne moyenne.
Les caractéristiques de la poutre sont :
 ligne moyenne : droite ou à grand rayon de courbure.
 section droite (S) : constante ou variant progressivement.
 grande longueur par rapport aux dimensions transversales.
 existence d'un plan de symétrie.
1.2.3 Les actions mécaniques extérieures

Plan de symétrie : les forces extérieures seront situées dans le plan de symétrie de la poutre ou
alors disposées symétriquement par rapport à ce plan.
Types d'actions mécaniques extérieures : deux types d'actions mécaniques peuvent s'exercer
sur la poutre : charges concentrées ( F1 ou moment MC ) et les charges réparties p sur DE.
(exprimées en N/m).
Mc
F1
A D E
C B
p
Figure 3 –Types d'actions mécaniques extérieures
Les déformations étant petites devant les dimensions de la poutre, les actions s'exerçant sur celle-
ci seront calculées à partir du principe fondamental de la statique.
Les supports des forces : seront considérés comme constants après déformation.
O A
A'
fi g .3
Figure 4 –Les supports des forces sont considérés constants après déformation.
1.2.4 Les déformations

Hypothèse des petites déformations HPD
Les déformations sont petites par rapport aux dimensions de la poutre.

8
Navier & Bernoulli
Les sections planes normales aux fibres avant déformation demeurent planes et normales aux
fibres après déformation.
Barré de St Venant
Les résultats obtenus par la RDM ne s'appliquent valablement qu'à une distance suffisamment
éloignée de la région d'application des efforts concentrés.
O A
fig.4
Figure 5 –Poutre après déformation
1.3 Principe fondamental de la statique

Un solide indéformable E, sollicité par un système de forces extérieures, est en équilibre si et
seulement si :

E  E  
 F ext  0 1.1
 M F ext / M  0
NB: le choix du point M est uniquement guidé par des raisons de commodité.
Equations d’équilibre
Système 3D (Oxyz) Système 2D (Oxy)
 proj. sur Ox des Fext=0  proj. sur Ox des Fext=0

R=0  proj. sur Oy des Fext=0  proj. sur Oy des Fext=0
 proj. sur Oz des Fext=0
 moments/ Ox des Fext=0

MR/M=0  moments/ Oy des Fext=0
 moments/ Oz des Fext=0  moments/ Oz des Fext=0
Tableau 1. Equations d’équilibre

9
Les actions de liaison
Efforts Efforts
Nom Ddl Ddl interdits transmis transmis Schema
3D 2D
Cas de flexion
𝑋𝐴 𝐿𝐴 𝑋𝐴 0
3 translations
Encastrement 0 { 𝑌𝐴 ; 𝑀𝐴 } { 𝑌𝐴 ; 0 }
3 rotations
𝑍𝐴 𝑁𝐴 𝐴 0 𝑚𝐴 𝐴
𝑋𝐴 𝐿𝐴 𝑋𝐴 0
3 translations
Articulation 1 { 𝑌𝐴 ; 𝑀𝐴 } { 𝑌𝐴 ; 0}
2 rotations
𝑍𝐴 0 𝐴 0 0𝐴
0 0 0 0
Appui simple 5 1 translation {𝑌𝐴 ; 0} {𝑌𝐴 ; 0}
0 0𝐴 0 0𝐴
Tableau 2. Actions de liaison
1.4 Torseur des efforts de cohésion

1.4.1 Définition
Soit une poutre (E) en équilibre sous l'action de n actions extérieures. On associe à cette poutre
un repère R (x,y,z) dont l'axe x coïncide avec la ligne moyenne de la poutre.
Coupons la poutre (E) par un plan (P) orthogonal à sa ligne moyenne, situé à l'abscisse x. On
définit ainsi deux portions de poutre (E1) et (E2).
Y
x
G X
Z E1 E2
Figure 6 –Poutre (E) coupée en deux portions (E1) et (E2).

(E) étant en équilibre, on peut écrire : E  E  0
(E1) étant en équilibre, on peut écrire : E  E1   E 2  E1  0

10
(E2) étant en équilibre, on peut écrire : E  E 2   E1  E 2  0
On en déduit :
 E 2  E1  E  E1  E  E 2

Le torseur  E 2  E1 qui traduit l'action de contact de (E2) sur (E1) est le torseur de cohésion
ou torseur des efforts intérieurs.
Ces efforts intérieurs permettent à la poutre de ne pas se "disloquer" sous l'effet d'actions
extérieures.
La RDM vise en particulier à vérifier qu'en aucun point de la poutre, les efforts de cohésion ne
soient supérieurs aux capacités du matériau.
On note :
Cohésion  E  E 2  E  E1 1.2
1.4.2 Composantes du torseur de cohésion

Dans le torseur de cohésion, on peut faire apparaître la résultante et le moment qui dépendent de
la position de la section (x).
Ty
R My
MG
G G
N MT
Tz Mz
Figure 7 – Résultante et moment du torseur de cohésion dans une section droite.

 N Mt
 R( x )   
 Cohésion     Avec : R Ty et MG Mfy
G  M G ( x ) 
R
Tz R
Mfz
 N : effort normal dans la section S,

 Ty et Tz : efforts tranchants dans la section S,
 MT : moment de torsion dans la section S,
 Mfy et Mfz : moments de flexion dans la section S,
Suivant les éléments de réduction non-nuls du torseur de cohésion : N, Ty, Tz, MT, Mfy et Mfz,
on peut alors identifier le type de sollicitation que subit la poutre, à savoir :

11
Composantes Sollicitation
>0 Extension (traction)
N
<0 Compression
Ty
Cisaillement
Tz
Mt Torsion
Mfy
Flexion
Mfz
Tableau 3. Type de sollicitation dans une poutre
Lorsque l’on a une seule de ces sollicitations on parle de sollicitation simple, sinon on a un
problème de sollicitations composées.
1.5 Tenseur des contraintes

1.5.1 Vecteur contrainte en un point M
Soit la poutre suivante :
(a) (b) (c)
E2
E1 E1 E1
Figure 8 –coupure, facette et vecteur contrainte dans un solide.

Effectuons une coupure du solide, par une surface S, en deux parties E1 et E2 (Figure 8. a)
La partie E1, par exemple, est en équilibre sous l’action de (Figure 8. b):
 des forces extérieures qui lui sont directement appliquées
 et des forces intérieures reparties sur la coupure.
Considérons un point M de S et soient dS un élément de la surface infiniment petit entourant M
et 𝑛⃗ le vecteur unitaire, normal en M à S et dirigé vers l’extérieur de la partie E1 (Figure 8. c).
⃗⃗⃗⃗ la force qui s’exerce sur la facette 𝑛⃗ en M.
Soit dF
On appelle vecteur contrainte sur la facette 𝑛⃗ en M, la quantité :

12
𝑑𝐹
⃗ (𝑀, 𝑛⃗) =
𝑇 1.3
𝑑𝑆
Une contrainte s’exprime en pascal (1 Pa = 1 N/m2). Dans la pratique, on utilise le Mégapascal
(1 MPa = 106 Pa = 1 N/mm2)
1.5.2 Contrainte normale/ contrainte tangentielle :

⃗ (𝑀, 𝑛⃗) qui s'exerce en M sur la surface ds de normale 𝑛⃗.
Soit la contrainte 𝑇
Figure 9 –Contrainte normale et contrainte tangentielle

⃗ (𝑀, 𝑛⃗) se décompose en :
Le vecteur contrainte 𝑇
 ⃗ (𝑀, 𝑛⃗) sur la normale 𝑛⃗ de la facette,
Contrainte normale 𝜎 : c'est la projection de 𝑇
 Contrainte tangentielle 𝜏 , appelée aussi contrainte de cession ou contrainte de cisaillement :
⃗ (𝑀, 𝑛⃗) sur le plan de la facette.
c'est la projection de 𝑇
On a alors :
𝜎𝑛 = 𝑇⃗ (𝑀, 𝑛⃗). 𝑛⃗
⃗ (𝑀, 𝑛⃗) = 𝜎 𝑛⃗ + 𝜏
𝑇 Avec : 1.4
⃗ (𝑀, 𝑛⃗)‖² − 𝜎𝑛 ²
‖𝜏𝑛 ‖² = ‖𝑇
 Si 𝜎𝑛 est positive, on a localement un état de traction,
 Si 𝜎𝑛 est négative, on a localement un état de compression.
1.5.3 Tenseur des contraintes
̿(M) sur les trois facettes du point M

Figure 10 –les composantes du tenseur de contraintes σ

13
Considérons le cube infiniment petit d’un solide, centré sur le point M, construit sur les axes x, y
et z.
Le tenseur des contraintes de Cauchy en M, dans le repère { ı, ȷ, k}:
𝜎𝑥𝑥 𝜎𝑥𝑦 𝜎𝑥𝑧
𝜎̿(𝑀) = (𝜎𝑦𝑥 𝜎𝑦𝑦 𝜎𝑦𝑧 ) 1.5
𝜎𝑧𝑥 𝜎𝑧𝑦 𝜎𝑧𝑧 (𝑥,𝑦,𝑧)
Contraintes en un point sur les facettes 𝐢, 𝐣, 𝐤

Les colonnes du tenseur des contraintes 𝜎̿(𝑀) sont les vecteurs contrainte sur les facettes
⃗ telle que (Figure 10) :
𝑖, 𝑗 , 𝑘
sur la facette 𝑖 : 𝑇 ⃗ = 𝜎̿(𝑀). 𝑖

⃗ (𝑀, 𝑖) = 𝜎𝑥𝑥 𝑖 + 𝜎𝑦𝑥 𝑗 + 𝜎𝑧𝑥 𝑘
sur la facette 𝑗 : 𝑇 ⃗ = 𝜎̿(𝑀). 𝑗

⃗ (𝑀, 𝑗) = 𝜎𝑥𝑦 𝑖 + 𝜎𝑦𝑦 𝑗 + 𝜎𝑧𝑦 𝑘 1.6
⃗ :
sur la facette 𝑘 ⃗ ) = 𝜎𝑥𝑧 𝑖 + 𝜎𝑦𝑧 𝑗 + 𝜎𝑧𝑧 𝑘
⃗ (𝑀, 𝑘
𝑇 ⃗ = 𝜎̿(𝑀). 𝑘
⃗
Soit en notation indicielle, le vecteur contrainte sur la facette 𝑖 :
𝑇𝑖 = 𝜎𝑗𝑖 𝑛𝑗 1.7
⃗⃗
Contraintes en un point sur une facette quelconque 𝒏
Figure 11 –Contraintes sur une facette 𝑛⃗ quelconque.

En général, sur une facette 𝑛⃗, le vecteur contrainte est (Figure 11) :
⃗ (𝑀, 𝑛⃗) = [𝜎(𝑀)]𝑛⃗
𝑇 1.8
La contrainte normale 𝜎𝑛 sur la facette 𝑛⃗ en M est égale à :

14
⃗ (𝑀, 𝑛⃗) = {𝑛}𝑇 [𝜎(𝑀)]{𝑛}
𝜎𝑛 = 𝑛⃗. 𝑇 1.9
La contrainte tangentielle ⃗⃗⃗⃗

𝜏𝑛 sur la facette 𝑛⃗ en M est telle que :
⃗ (𝑀, 𝑛⃗)‖ = 𝜎𝑛 ² + 𝜏𝑛 ²
‖𝑇 1.10
1.5.4 Torseur de cohésion : écriture locale

Le torseur des actions mécaniques s’exerçant sur dS au point M et au point G :
𝑑𝐹 (𝐸2 → 𝐸1) 𝑑𝐹 (𝐸2 → 𝐸1) ⃗ (𝑀, 𝑛⃗)𝑑𝑆

𝑇
{ } ={ } ={ } 1.11
⃗0 𝑀
⃗⃗⃗⃗⃗⃗
𝐺𝑀 ∧ 𝑑𝐹 (𝐸2 → 𝐸1) 𝐺 ⃗⃗⃗⃗⃗⃗
𝐺𝑀 ∧ 𝑇 ⃗ (𝑀, 𝑛⃗)𝑑𝑆
𝐺
Le torseur des efforts intérieurs sur la surface totale S :
⃗ (𝑀, 𝑛⃗)𝑑𝑆
∬ 𝑇
𝑅⃗ (𝐸2 → 𝐸1) 𝑆
{𝑇𝑖𝑛𝑡 } = {𝑇(𝐸2→𝐸1) } = { } = 1.12
⃗⃗⃗⃗⃗
𝑀𝐺 (𝐸2 → 𝐸1) 𝐺
∬ ⃗⃗⃗⃗⃗⃗
𝐺𝑀 ∧ 𝑇 ⃗ (𝑀, 𝑛⃗)𝑑𝑆
{ 𝑆 }𝐺
Les actions mécaniques de cohésion sont modélisées par le torseur des efforts intérieurs de
cohésion caractérisé au point G, centre de la section droite.
Le torseur des efforts intérieurs de cohésion ne représente qu’une vision globale sur la section
droite de toutes les actions mécaniques qui s’appliquent localement en chaque point de la surface.

15
Chapitre 2
Traction / compression des poutres.
2.1 Définition
Une poutre est sollicitée à la traction (compression) simple lorsqu'elle est soumise à deux forces
directement opposées et qui tendent à l'allonger (comprimer).
B B
A A
R y
(S)
A R
A x
Figure 12 – Poutre sollicitée à la traction

Les éléments de réduction en G du torseur des efforts de cohésion s'expriment par :
N 0
Cohésion  0 
0 Avec N > 0 pour la traction et N < 0 pour la compression
0 
G
0( x , y , z )
2.2 Essai d'traction
Figure 13 –Eprouvette d'essai de traction

16
Une éprouvette en acier est sollicitée à l'traction par une machine d'essai, qui permet de
déterminer l'allongement de l'éprouvette en fonction de l'effort qui lui est appliqué.
Courbe d’essai de traction
F(N)
A
B
D
O
l (m m )
Figure 14 –Courbe d’essai de traction
Zone OA :
C'est la zone des déformations élastiques. Si l'on réduit la valeur de F jusqu'à une valeur nulle,
l'éprouvette retrouve sa longueur initiale.
Dans cette zone, l'allongement est proportionnel à l'effort d'traction.
Des essais effectués avec des éprouvettes de dimensions différentes permettent de constater que
pour un même matériau, l'allongement relatif ( ∆l / l0) est proportionnel à l'effort surfacique (F /
S0).
Les sections droites et planes de l'éprouvette restent droites et planes pendant l'essai.
Zone ABCD :
C'est la zone des déformations permanentes. Si l'on réduit la valeur de F jusqu'à une valeur nulle,
l'éprouvette ne retrouve pas sa longueur initiale.
On ne s'intéressera (pour l’instant) qu'à la zone des déformations élastiques.
2.3 Caractéristiques mécaniques d'un matériau

Contrainte limite élastique en traction 𝝈𝒆 ou 𝑹𝒆
C'est la valeur limite de la contrainte dans le domaine élastique, appelée aussi limite d'élasticité.
𝐹𝑒
𝑅𝑒 = 2.1
𝑆
Pour l'acier, 𝑅𝑒 est voisine de 300 MPa.

17
Contrainte limite de rupture en traction 𝝈𝐫 ou Rr
C'est la valeur limite de la contrainte avant rupture de l'éprouvette, appelée aussi résistance à la
traction.
𝐹𝑟
𝑅𝑟 = 2.2
𝑆
Pour l'acier, cette valeur est voisine de 480 MPa.
Allongement A%
L' L
A%  *100 2.3
L
Avec
L : longueur initiale de l'éprouvette.
L' : longueur de l'éprouvette à sa rupture.
Pour l'acier, on constate des valeurs de A% voisines de 20%.
Loi de HOOKE : Module d'élasticité longitudinal ou module de Young E

La propriété constatée lors de l'essai de traction pour différents matériaux a permis d'établir la la
loi de Hooke :
𝑁 𝐿′ − 𝐿 ∆𝐿
𝜎 = 𝐸. 𝜀 Avec : 𝜎= 𝜀= = 2.4
𝑆 𝐿 𝐿
N : effort normal de traction en Newton. S : section en mm2,
L=L'- L : Allongement de la poutre en en mm,  : Allongement relatif,
𝜎: contrainte normale (traction 𝜎 > 0, compression 𝜎 < 0) en MPa.
2
E en MPa (N/mm ) : module d'élasticité longitudinal ou module de Young.
Matériau Fontes Aciers Cuivre Aluminium Tungstène
E (MPa) 60000à160000 200000 120000 70000 400000
Tableau 4. Exemples de valeurs du module de Young.
Coefficient de Poisson 
Lors de cet essai, on met aussi en évidence une autre caractéristique de l’élasticité ; il existe un
rapport constant entre la contraction relative transversale (∆𝑐 / c) et l'allongement relatif
longitudinal (∆𝐿 / L) qui s’écrit :
∆𝑐 ∆𝐿
= −𝜈 Ou bien encore 𝜀𝑇 = −𝜈. 𝜀𝐿 2.5
𝑐 𝐿
𝜈 : coefficient de Poisson (sans unité), il est de l'ordre de 0,3 pour les métaux.

18
Figure 15 –Mise en évidence du rapport entre la contraction relative transversale et l'allongement
relatif longitudinal.
2.4 Contraintes et déformations

E1
y
(S )
R = N .x
F R
x
G
fi g .8
z
Figure 16 –Petit tronçon de poutre de longueur dx.
Lorsqu’une poutre est sollicitée à la traction, la contrainte tangentielle est nulle. La contrainte de
traction devient (équation 1.4 ).
⃗ (𝑀, 𝑥) = 𝜎𝑥
𝑇 2.6
Le tenseur des contraintes et de déformation en M, dans le repère (𝑥, 𝑦, 𝑧) :
𝜎 0 0 𝜀 0 0
𝜎̿(𝑀) = (0 0 0) 𝜀̿(𝑀) = (0 −𝜈𝜀. 0 ) 2.7
0 0 0 (𝑥,𝑦,𝑧) 0 0 −𝜈𝜀. (𝑥,𝑦,𝑧)
2.5 Déformation d'un tronçon infiniment petit

Soit une poutre en traction de longueur L.
Prenons un petit tronçon de poutre infiniment petit de longueur dx et situé à une distance x de O :
L'allongement relatif (x) du petit tronçon de la poutre de longueur dx est :
𝑢(𝑥 + 𝑑𝑥) − 𝑢(𝑥) 𝑑𝑢(𝑥)
𝜀(𝑥) = = donc 𝑑𝑢(𝑥) = 𝜀(𝑥)𝑑𝑥 2.8
𝑑𝑥 𝑑𝑥
L’allongement total de la poutre est :
𝐿 𝐿 𝐿
∆𝐿
∆𝐿 = ∫ 𝑑𝑢 = ∫ 𝜀(𝑥)𝑑𝑥 = 𝜀 ∫ 𝑑𝑥 = 𝜀𝐿 D’ou 𝜀=
0 0 0 𝐿

19
Figure 17 –Petit tronçon de poutre de longueur dx.
2.6 Torseur des efforts intérieurs

Le torseur des efforts intérieurs se calcule : (équation 1.12)
⃗ (𝑀, 𝑛⃗)𝑑𝑆
∬ 𝑇 ∬ 𝜎𝑥𝑑𝑆
𝑁𝑥 𝑆 𝑆
{𝑇𝑖𝑛𝑡 } = { } = =
⃗0 𝐺
∬ ⃗⃗⃗⃗⃗⃗
𝐺𝑀 ∧ 𝑇 ⃗ (𝑀, 𝑛⃗)𝑑𝑆 ∬ ⃗⃗⃗⃗⃗⃗
𝐺𝑀 ∧ 𝜎𝑥𝑑𝑆
{ 𝑆 }𝐺 { 𝑆 }𝐺
Si  est contant sur S, on a alors :
𝑁
𝑁𝑥 = ∬ 𝜎𝑥 𝑑𝑆 = 𝜎𝑥 ∬ 𝑑𝑆 = 𝜎𝑆𝑥 soit 𝜎=
𝑆 𝑆 𝑆
Encore, on vérifie que :
⃗⃗⃗⃗⃗⃗ ∧ 𝜎𝑥𝑑𝑆 = 0
∬ 𝐺𝑀
𝑆
Au bilan, La seule contrainte non nulle dans la section droite de normale x , de surface S, est la
contrainte normale.
2.7 Condition de résistance

Pendant toute la durée de son service, une pièce doit conserver un comportement élastique. Cette
condition s'exprime par :
𝜎𝑚𝑎𝑥𝑖 ≤ 𝜎𝑒 = 𝑅𝑒 2.9
L'incertitude sur la valeur de la contrainte Re nous amène à réécrire la condition de résistance :
𝜎𝑒 𝑅𝑒
𝜎𝑚𝑎𝑥𝑖 ≤ 𝜎𝑝𝑒 = = 𝑅𝑝𝑒 = 2.10
𝑠 𝑠
 𝑅𝑒 ou 𝜎𝑒 : résistance limite élastique en MPa
 𝑅𝑝𝑒 ou 𝜎𝑝𝑒 : résistance pratique de limite élastique en MPa
 s : coefficient de sécurité (s>1)
20
2.8 Concentration de contraintes
Des pièces industrielles peuvent posséder des singularités de formes (perçages, gorges, rainures,
filetages…). Ces singularités font augmenter la contrainte.
Figure 18 –Concentration de contraintes au niveau du perçage sur la poutre

On définit un coefficient de concentration de contrainte Kt tel que la contrainte maximale a pour
valeur :
𝜎𝑚𝑎𝑥𝑖 ≤ 𝑘𝑡 . 𝜎𝑛𝑜𝑚 2.11
 𝜎𝑚𝑎𝑥𝑖 = contrainte atteinte au voisinage de la singularité

 𝜎𝑛𝑜𝑚 = contrainte moyenne nominale calculée
Exemple 1
Pour un filetage triangulaire ISO : Kt = 2.5
Exemple 2
Arbre de section circulaire épaulé.

21
2.9 Triés
Les treillis sont des structures formés de barres assemblées de façon à former des triangles
indéformables.
Figure 19 –Tout polygone autre qu'un triangle est déformable.

 Permettent d'alléger l'ensemble d'une construction tout en assurant une plus grande stabilité,
 Largement utilisés en construction : fermes de toiture, grues, ponts roulants, pylônes, etc.
 Chaque articulation ou nœuds de la structure est une rotule parfaite : on parle de nœud
articulé ou rotulé, par opposition au nœud rigide;
 Les charges extérieures sont supposées appliquées aux nœuds de la structure : Les charges
sur les éléments du treillis sont nulles ou faibles devant les charges nodales.
 Les éléments du treillis ne travaillent donc qu'en traction- compression.
 Dans le cas contraire il faudra prendre en compte la flexion des éléments, nous parlerons
dans ce cas de portiques.
2.9.1 Efforts sur la structure
Figure 20 –Efforts externes (___) et Efforts internes (- - -) dans la structure

 Efforts externes sont les charges à supporter : le poids propre de la structure, le poids de
la neige, le trafic, les réactions d'appuis, ...

22
 Efforts internes : efforts dans les barres,
Ces efforts doivent être déterminés pour pouvoir choisir les matériaux requis dans la réalisation
des constructions.
Les pièces trop grosses ne sont pas économiques et les pièces trop petites ne sont pas sécuritaires.
2.9.2 Efforts internes: effort normal dans les barres

L’équilibre d’une barre du treillis (AB) montre :
 Que les barres sont soumises seulement à un effort normal (Traction ou compression)
 l'absence de moments fléchissant et d'efforts tranchants dans les barres..
Figure 21 –Mise en évidence de l'effort normal dans les barres.
2.9.3 Barre à effort normal nul N=0

(a) (b)
Figure 22 –Barres à effort normal.

Deux cas de barres ou l'effort normal est nul :
Cas 1 (Figure 22 .a) : Si deux barres concourent en un nœud non chargé, l'effort normal est nul
dans ces barres.
Considérons par exemple la barre BC et raisonnons par l'absurde en supposant que N 0.
 N (BC) possèdera une composante perpendiculaire à la barre AB selon la ligne en pointillés.
 Or, comme le nœud n'est pas chargé, cette composante n'est équilibrée par aucune force.

23
 L'effort correspondant est donc forcément nul.
Cas 2 (Figure 22 .b) : Une barre (3) joignant, en un nœud i non chargé, deux autres barres
alignées (1) et (2), est le siège d'un effort normal N=0. De plus Na=Nb.
Cette propriété s'explique de la même façon que dans le cas précédent
2.9.4 Déformations des nœuds

La déformation au nœud k du treillis peut être calculée en utilisant le théorème énergétique de
Castigliano (voir 6.5) :
𝑏
∂Wd 𝑁 ∂𝑁 Ni L i 𝜕Ni
δk = [ ] =∫ ds = ∑ n Avec ni = 2.12
∂Pk Pk =0
𝑡𝑟𝑒𝑖𝑙𝑙𝑖𝑠 𝐸𝑆 ∂Pk E i Si i 𝜕Pk
𝑆i Pk fictive 𝑖=1
N.B : Prendre Pk = 0 seulement si Pk est une charge fictive.
Figure 23 –Déformation des nœuds.
2.9.5 Nœud articulé et nœud rigide

La plupart des treillis sont construits avec des nœuds rigides plutôt qu'articulés : Nœuds rigides
obtenus par soudage ou par assemblages boulonnés. Dans ce cas :
 des contraintes de flexion apparaissent car les barres doivent fléchir pour suivre le
déplacement des nœuds au lieu de tourner librement autour de leurs extrémités.
 les efforts normaux sont quasi identiques;
 les déplacements des nœuds sont quasi identiques; mais l'allure de la déformée est différente,
 des contraintes parasitaires de flexion se produisent et peuvent être importantes dans certains
cas.
2.9.6 Hyperstaticité
Un treillis est une structure discrète constituée de b barres reliées entre elles aux n noeuds de la
structure.
24
Le degré d'hyperstaticité de la structure h est :
h = Ni - Ne
 Ne : nombre d'équations d'équilibre = 2n (si 2D) ou 3n (si 3D), n: nbre de nœuds,
 Ni : nombre d'inconnues = b +r +m
o b : les efforts normaux dans les b barre du treillis,
o r : Les réactions aux appuis,
o m : le nombre de mouvements d'ensemble (s'il y en a).
Le degré d'hyperstaticité de la structure h devient :
h = b +r +m - 2n 2.13
Exemples
(a) (b)
Figure 24 –Exemples de triés
Cas de la Figure 24.a

h = Ni – Ne=0
Cette structure possède 4 nœuds, le problème est plan, soit : Ne= 8
Il y a 3 mouvements rigides, et 5 barres donc 5 efforts intérieurs inconnus, soit Ni =8
Cette structure est isostatique, les efforts dans les barres ne dépendent que de la géométrie, ils
sont indépendants des caractéristiques mécaniques des barres (section, matériau)
Cas de la Figure 24.b

Si on ajoute une barre sur l'autre diagonale, la structure sera hyperstatique de degré 1.
h = Ni – Ne=1
Sa rigidité sera plus élevée, et la distribution des efforts interne dépendra des caractéristiques
mécaniques des barres.

25
2.9.7 Résolution des triés isostatiques
Il existe plusieurs méthodes de résolution des treillis :
Méthode des sections ou de Ritter

La méthode des sections ou méthode de Ritter (Allemagne 1779-1859) : Consiste à isoler des
morceaux de structure bien choisis et à écrire leurs équations d'équilibre;
Pour les morceaux isolés : Il faut au maximum trois efforts inconnus, soit Trois barres au
maximum sont coupées.
Exemple :
Figure 25 –Principe de la méthode de Ritter
3 efforts inconnus 3 équations d’équilibre disponibles : On obtient :
N1+ N2 cos  + N3 + Rx = 0
N2 =R/sin 
-N2 sin  + Ry = 0
N1, N2, N3 N3 = -R d/h
N3 h+ R d = 0
N1 =R d/h -R/tg(a)
Avec : Rx = 0, Ry = R

26
Méthode des nœuds ou de Cremona
La méthode des nœuds ou de Cremona (Italie1830- 1903) : Consiste à écrire les équations
d’équilibre pour chaque nœud; l'un après l'autre.
Pour chaque nœud : Il faut au maximum deux efforts inconnus,
Exemple : Nœud O
Figure 26 –Principe de la méthode de Cremona

2 efforts inconnus Equations d’équilibre du Nœud O On obtient :
N1 cos(a)+N2=0 N1 =-R/sin(a)
N1, N2
N1 sin(a)+R=0 N2= -N1 cos(a)= R/tg(a)
Méthode générale
Figure 27 – Principe de la méthode Générale

Exprime l'équilibre de tous les nœuds. Nécessite la résolution d'un système de 2n équations (2
équations par nœud si 2D). Cette méthode est plus intéressante lors d'un calcul par ordinateur.
Prenons le noeud i (i=1 à n) reliant b barres (j=1 à b) de la Figure 27.
L'équilibre du nœud i s’écrit si pas d’efforts extérieurs ni réaction d’appui :

27
𝑥𝑗 − 𝑥𝑖
∑ 𝑁𝑗 𝑐𝑜𝑠 𝛼𝑗 = 0 ∑ 𝑁𝑗 =0
𝐿𝑗
𝑗 𝑗
soit
𝑦𝑗 − 𝑦𝑖
∑ 𝑁𝑗 sin 𝛼𝑗 = 0 ∑ 𝑁𝑗 =0
𝐿𝑗
{ 𝑗 { 𝑗
Si présence d’efforts extérieurs et réaction d’appui, on obtient les 2n équations suivantes :
𝑥𝑗 − 𝑥𝑖
𝑅𝑥𝑖 + 𝐹𝑥𝑖 + ∑ 𝑁𝑗 =0
𝐿𝑗
𝑗
𝑦𝑗 − 𝑦𝑖
𝑅𝑦𝑖 + 𝐹𝑦𝑖 + ∑ 𝑁𝑗 =0
𝐿𝑗
{ 𝑗
2.9.8 Résolution des triés hyperstatiques

Méthode des déplacements
Pour les treillis hyperstatiques de degré h, on applique le théorème de Ménabréa (voir 6.5) aux
nœuds fixes pour ajouter les h équations manquantes :
𝑏
𝜕Wd 𝑁 𝜕𝑁 Ni Li
δk = =∫ ds = ∑ n =0 Avec k varie de 1 à h 2.14
𝜕Xk 𝑡𝑟𝑒𝑖𝑙𝑙𝑖𝑠 𝐸𝑆 𝜕Pk Ei Si i
𝑖=1
Méthode des éléments finis

Voir cours de Calcul de structures par la méthode
des éléments finis (Chapitre 9).
2.10 Mesures par jauges de contraintes
Figure 28 – Mesure de déformations par jauges de contraintes.

En points A et A', on colle deux jauges (j) et (j’) identiques, dans la direction axiale X’X, et on en
fait les résistances R1 et R4 du demi pont de Wheatstone.

28
Les résistances R2 et R3 sont constituées par deux jauges de compensation thermique (identiques
à (j) et (j’)), le pont est insensible aux variations de température. Vi désignant la tension
d’alimentation et VO la tension de déséquilibre du pont qui apparaît lorsque les jauges se dilatent
sous l’action du chargement. Les déformations se calculent (voir annexe A) :
𝑉0 2 𝑉0 𝑁
𝜀1 + 𝜀4 = 4 ou puisque : 𝜀1 = 𝜀4 𝜀1 = = 2.15
𝐾𝑉𝑖 𝐾 𝑉𝑖 𝐸𝑆
Connaissant le facteur de jauge et la mesure de la tension 𝑉0, on peut déduire la déformation
axiale 𝜀1 et l'effort normal 𝑁.

29
Chapitre 3
Cisaillement des poutres.
3.1 Définition
Une poutre subit une sollicitation de cisaillement simple lorsqu'elle est soumise à deux systèmes
d'action de liaison qui se réduisent dans un plan (P) perpendiculaire à la ligne moyenne à deux
forces directement opposées.
(E) (P)
F
F'
B
Figure 29 –Poutre sollicitée au cisaillement simple.

Sous l'action de ces deux forces la poutre tend à se séparer en deux tronçons E1 et E2 glissant l'un
par rapport à l'autre dans le plan de section droite (P).
y (S)
(E1 )
F
E1
T x
E2
G
z
F'
(P)
Figure 30 –Glissement du tronçon E2 dans le plan de la section droite.

 0 0
 
Cohésion Ty 0
Tz 0   
G ( x , y , z )

30
3.2 Essai de cisaillement
Il est physiquement impossible de réaliser du cisaillement pur au sens de la définition précédente.
Les essais et résultats qui suivent permettent toutefois de rendre compte des actions tangentielles
dans une section droite et serviront ainsi dans le calcul de pièces soumises au cisaillement.
On se gardera cependant le droit d'adopter des coefficients de sécurités majorés pour tenir compte
de l'imperfection de la modélisation.

Considérons une poutre (E) parfaitement encastrée et appliquons-lui un effort de cisaillement F
uniformément réparti dans le 0)
d'encastrement (voir fig.).
y
x
B
(E 1 ) G (E 2 )
A x
(S ) (S 0 )
(P )
Figure 31 –Essai de cisaillement

Si l'on isole (E1), on trouve alors le torseur de cohésion suivant :
 0 0 
 
Cohésion  F 0 
 0 F .x 
G ( x , y , z )
Analyse de la courbe obtenue

F(N)
B
x
y
C
(S 0 )
A
(S )
F
O
y (mm)
Figure 32 –Courbe d'essai de cisaillement

31
Zone OA
C'est la zone des déformations élastiques. Si l'on réduit la valeur de F jusqu'à une valeur nulle,
l'éprouvette retrouve sa forme initiale
Zone ABC
C'est la zone des déformations permanentes. Si l'on réduit la valeur de F jusqu'à une valeur nulle,
l'éprouvette ne retrouve pas sa forme initiale (déformations plastiques).
3.3 Caractéristiques mécaniques d'un matériau

3.3.1 Contrainte tangentielle limite élastique 𝛕𝒆 ou Rge
C'est la valeur limite de la contrainte dans le domaine élastique.
𝐹𝑒
𝜏𝑒 = 3.1
𝑆
Pour l'acier, cette valeur est comprise entre 250 MPa et 600 MPa.
3.3.2 Contrainte tangentielle de rupture

𝐹𝑟
𝜏𝑟 = 3.2
𝑆
C'est la valeur limite de la contrainte avant rupture de l'éprouvette.
3.4 Contrainte tangentielle de cisaillement

Lorsqu’une poutre est sollicitée au cisaillement, la contrainte normale est nulle. La contrainte de
cisaillement devient (équation 1.4 ).
⃗ (𝑀, 𝑛⃗) = 𝜏 = 𝜏𝑦
𝑇 3.3
Le tenseur des contraintes de Cauchy en M, dans le repère (𝑥, 𝑦, 𝑧) :
0 𝜏𝑥𝑦 0
𝜎̿(𝑀) = (𝜏𝑥𝑦 0 0) 3.4
0 0 0 (𝑥,𝑦,𝑧)
La contrainte tangentielle moyenne de cisaillement vaut :
𝑇
𝜏= 3.5
𝑆
 : Valeur moyenne de la contrainte tangentielle de cisaillement en MPa.
T : effort tranchant en Newton.
S : aire de la section droite (S) en mm2.

32
3.5 Relation entre contrainte et déformation
L'essai précédent a permis pour différents matériaux d'établir la relation :
𝑇 ∆𝑦
=𝐺 3.6
𝑆 ∆𝑥
F en Newton, S en mm , G en MPa, ∆𝑦 et ∆𝑥 en mm.
2
∆𝑦
𝛾 = ∆𝑥 : est appelé glissement relatif.
G est une caractéristique appelée module d'élasticité transversal ou module de Coulomb.
D'après cette relation de l'essai, on peut écrire la loi de comportement entre la contrainte et la
déformation :
𝜏=𝐺𝛾 3.7
Valeurs du module de Coulomb pour quelques matériaux.
Matériau Fontes Aciers Laiton Duralumin Plexiglas
G (MPa) 40000 80000 34000 32000 11000
Tableau 5. Exemples de valeurs du module de Coulomb.
3.6 Condition de résistance

à une valeur limite
p.
𝑇𝑀𝐴𝑋 𝜏𝑒
𝜏𝑀𝐴𝑋 = < 𝜏𝑝 = 3.8
𝑆 𝑠
s est un coefficient de sécurité qui varie de 1,1 à 10 selon le domaine d'application.

33
Chapitre 4
Torsion des poutres cylindriques.
4.1 Hypothèse supplémentaire

Les poutres étudiées sont des cylindres de révolution à base circulaire.
Pour toute section qui n’est pas circulaire :
 les sections ne restent pas planes et se gauchissent,
 la contrainte de cisaillement ne peut pas être tangente au contour non circulaire de la section.
Figure 33 –Section circulaire et section non circulaire.
4.2 Définition
Une poutre est sollicitée en torsion simple lorsqu'elle est soumise à ses deux extrémités à des
liaisons dont les torseurs associés se réduisent à deux torseurs couples opposés dont les moments
sont parallèles à l'axe du cylindre.
MG1 G2 MG2
G1
R y
(S)
MG1
MG
x
G1

34
Figure 34 –Cylindre soumis à la torsion.
0 Mt 
Cohésion 0 0  4.1
0 0    
G ( x , y , z )
4.3 Essai de torsion

Un dispositif permet d'effectuer un essai de torsion sur une poutre encastrée à son extrémité G 1 et
soumise à un torseur couple à son extrémité G2.
Cette machine permet de tracer le graphe du moment appliqué en G2 en fonction de l'angle de
rotation d'une section droite.
Figure 35 –Essai de torsion.
avec x, l'abscisse de la section droite étudiée

Ou bien encore :
𝛼𝐿 𝛼 𝛼𝑖
= = = ⋯ = 𝑐𝑡𝑒
𝐿 𝑥 𝐿𝑖
On pose :
𝛼 𝛼 + 𝑑𝛼 𝑑𝛼
𝜃= = = = 𝑐𝑡𝑒 4.2
𝑥 𝑥 + 𝑑𝑥 𝑑𝑥
 𝜃 : angle de torsion unitaire (rad/mm)

 𝛼 : angle total de torsion de (S)/(S ) (rad)
0
 x : distance entre (S) et (S ) (mm)

0

35
Graphe des moments en fonction de l'angle de rotation
M G 2 (m N ))
O 
Figure 36 –Graphe de l'essai de torsion.
Zone OA :
C'est la zone des déformations élastiques. Si l'on réduit la valeur du moment jusqu'à une valeur
nulle, l'éprouvette retrouve sa forme initiale.
Dans cette zone, l'angle de torsion α est proportionnel au couple appliqué.
Les sections droites et planes de l'éprouvette restent droites et planes pendant l'essai.
Zone AB :
C'est la zone des déformations permanentes.
L'éprouvette ne retrouve pas sa forme initiale après déformation.
4.4 Déformations élastiques : distorsion

On pose :
𝑀𝐿 𝑀′𝐿 𝑀𝑀′ 𝑀𝑀′
= = = 𝑡𝑔𝛾 ≈ 𝛾
𝑀𝐿 𝑀0 𝑀𝑀0 𝑥
Encore :
𝑀𝑀′ 𝜌𝛼1 𝜌𝛼
𝛾= = = = 𝜌𝜃 4.3
𝑥 𝐿1 𝑥
 : représente l’angle de glissement de (S)/(S0) ou distorsion.
4.5 Contraintes et déformations de torsion

On coupe le cylindre en une section (S) et on exprime que la partie isolée est en équilibre sous
l’action du moment de torsion MG ou Mt et des forces de cohésion dans la section (S).
dS : élément de surface situé à une distance r de l’axe Gx, soumis à une contrainte de torsion 𝜏𝑀 .

36
(E1)
y
R (S) M a x
MG1 dF   .dS
M
MG M
G1 x 
G
G v
(S )
z
Figure 37 –Contrainte de cisaillement dans la section (S).
La contrainte de cisaillement devient (équation 1.4 ).
⃗ (𝑀, 𝑛⃗) = 𝜏 = 𝜏𝑒⃗⃗⃗⃗𝜃
𝑇 4.4
Le tenseur des contraintes et de déformation en M, dans le repère (𝑟, 𝜃, 𝑥) :
0 0 0
0 0 0 𝛾𝑥𝜃
𝜎̿(𝑀) = (0 0 𝜏𝑥𝜃 ) 0 0 𝜀𝑥𝜃 =
𝜀̿(𝑀) = 2 4.5
0 𝜏𝑥𝜃 0 (𝑟,𝜃,𝑥) 𝛾𝑥𝜃
0 𝜀𝑥𝜃 = 0
( 2 )(𝑟,𝜃,𝑥)
Avec : 𝜏𝑥𝜃 = 𝜏.
La loi de comportement permet d'écrire :
𝜏𝑥𝜃 = 2𝜇𝜀𝑥𝜃 = 2𝐺𝜀𝑥𝜃 = 𝐺𝛾𝑥𝜃
Dans la zone élastique, l'angle de torsion 𝛼 est proportionnel au couple appliqué. Egalement, la
distorsion 𝛾 est proportionnelle à la contrainte de cisaillement 𝜏. Cela se traduit par la loi de
comportement :
𝜏 = 𝐺. 𝛾 = 𝐺. 𝜃. 𝜌 4.6
G: module d’élasticité transversal ou module de Coulomb

L’effort élémentaire de cisaillement dF vaut donc : 𝑑𝐹 = 𝜏. 𝑑𝑠.
L’équilibre de l’élément isolé s’écrit donc:
𝑀𝑡 = ∬ 𝜌2 Gθds
𝑆
Comme 𝐺. 𝜃 est identique pour chaque section, on obtient finalement :
𝑀𝑡 = ∬ 𝜌2 Gθds = Gθ ∬ 𝜌2 ds = Gθ𝐼0
𝑆 𝑆
𝐼0 est le moment d’inertie polaire de (S)/ à G en mm 4
On peut écrire alors :

37
𝑀𝑡 = Gθ𝐼0 4.7
ou encore
𝑀𝑡 𝑀𝑡
𝜏= 𝜌=
𝐼0 𝐼0 4.8
𝜌
𝐼0
est module de torsion en mm3
𝜌
La contrainte de cisaillement est donc proportionnelle au rayon et elle est maxi pour 𝜌 = 𝑅
𝑀𝑡
𝜏 𝑀𝑎𝑥𝑖 =
𝐼0 4.9
𝑅
Figure 38 –Contrainte de cisaillement dans la section (S) en fonction du rayon.
4.6 Conditions de résistance

ure à une valeur limite
appelée contrainte pratique p :
𝜏𝑒
𝜏𝑀𝐴𝑋 𝑟é𝑒𝑙𝑙𝑒 = 𝑘𝑡 𝜏𝑀𝐴𝑋 < 𝜏𝑝 = 4.10
𝑠
𝑘𝑡 est le coefficient de concentration de contraintes,

𝜏𝑒 est la contrainte de limite élastique
Et s est le coefficient de sécurité
4.7 Condition de rigidité ou de déformation

La condition de résistance donne des arbres de petits diamètres très déformables et qui sont
source de vibrations et ne réalisent pas des transmissions rigides. La plus part des arbres se
calculent par une condition de déformation :
Mt M t .
  lim ou bien :    lim 4.11
G.I0 G.I0

38

Considérons l’arbre cylindrique de rayon Re et P un point de la surface latérale.
Les fibres axiales et circonférentielles ne subissent aucune déformation (ε = 0).
Les plus grandes déformations linéaires sont les déformations principales 𝜀𝐼 et 𝜀𝐼𝐼 des fibres
hélicoïdales inclinées à ± 45° par rapport aux génératrices.
En effet, la direction radiale (ou normale) est direction principale associée à la valeur propre zéro
et les directions 𝑋𝐼 et 𝑋𝐼𝐼 des fibres hélicoïdales inclinées à ± 45° par rapport aux génératrices
sont directions principales.
Les tenseurs de contraintes et de déformation s'écrivent au point P de la base principale
(𝑟, 𝑋𝐼 , 𝑋𝐼𝐼 ) :
0 0 0
0 0 0 𝛾𝑥𝜃
𝜎̿(𝑃) = (0 𝜏𝐼 0) 0 𝜀𝐼 = 0
𝜀̿(𝑃) = 2 4.12
0 0 𝜏𝐼𝐼 (𝑟,𝑋𝐼 ,𝑋𝐼𝐼) 𝛾𝑥𝜃
0 0 𝜀𝐼𝐼 = −
( 2 )(𝑟,𝑋 ,𝑋
𝐼 𝐼𝐼 )
Les déformations principales 𝜀𝐼 et 𝜀𝐼𝐼 des fibres hélicoïdales s'écrivent en utilisant la loi de
comportement (équation 4.6) :
𝛾𝑥𝜃 𝜏 𝑀𝑡
𝜀𝐼 = = = .𝑅
2 2𝐺 2𝐺𝐼0 𝑒
𝛾𝑥𝜃 𝜏 𝑀𝑡 4.13
𝜀𝐼𝐼 = − =− =− .𝑅
2 2𝐺 2𝐺𝐼0 𝑒
En pratique, pour une mesure précise de la distorsion 𝛾, on utilise un montage de quatre jauges en
pont complet qui élimine les dilatations pouvant provenir éventuellement de la température, de
l’effort normal N, du moment de flexion Mf et de l’effort tranchant T.
Un tel pont est constitué de quatre jauges identiques, notées j1, j2, j3, j4, de résistance nominale
R et de facteur de jauge K.
Les jauges j1 et j2 sont collées en P suivant les directions perpendiculaires 𝑋𝐼 et 𝑋𝐼𝐼 . Les jauges j4
et j3 sont collées en P’, diamétralement opposé à P, suivant les directions principales 𝑋𝐼 et 𝑋𝐼𝐼 de
ce point.
39
Elles sont câblées de façon à constituer les quatre résistances du pont complet de Wheatstone,
Au repos la tension VO est nulle. Sous l’action du chargement, l’arbre se déforme et les jauges
subissent des allongements relatifs ε1, ε2, ε3, ε4. Les déformations se calculent (voir annexe A) :
𝑉0 ou puisque : 𝛾𝑥𝜃 1 𝑉0 𝑀𝑡
𝜀1 − 𝜀2 − 𝜀3 + 𝜀4 = 4 𝜀1 = = = .𝑅 4.14
𝐾𝑉𝑖 𝜀1 = 𝜀4 = −𝜀2 − 𝜀3 2 𝐾 𝑉𝑖 2𝐺𝐼0 𝑒
Connaissant le facteur de jauge et la mesure de la tension 𝑉0, on peut déduire la distorsion 𝛾𝑥𝜃 et
le moment de torsion 𝑀𝑡 .

40
Chapitre 5
Flexion plane simple des poutres.
5.1 Définition
Une poutre est sollicitée à la flexion simple si les éléments de réduction au centre de gravité de
chaque section des forces de cohésion sont un effort tranchant et un moment de flexion.
y
ligne moyenne
z Mf
G
(S)
T
Plan de symétrie longitudinal et plan des charges
Figure 40 –Poutre soumise à la flexion plane simple.

0 0 
Cohésion Ty 0 

5.1
 0 Mfz    
G ( x , y , z )
5.2 Relation effort tranchant/moment fléchissant
Figure 41 –Tronçon de poutre de longueur dx.

41
Isolons le tronçon de poutre de longueur dx. Le bilan des actions mécaniques extérieures sur le
tronçon donne :
En G En G’
−𝑇𝑦 𝑦 (𝑇𝑦 + 𝑑𝑇𝑦 )𝑦
{𝑇𝑖𝑛𝑡 } = { } {𝑇𝑖𝑛𝑡 } = { }
−𝑀𝑓𝑧 𝑧 (𝑀𝑓𝑧 + 𝑑𝑀𝑓𝑧 )𝑧
𝐺 𝐺′
Le Principe Fondamental de la Statique au tronçon de poutre :
1
−𝑇𝑦 + 𝑇𝑦 + 𝑑𝑇𝑦 − 𝑝𝑑𝑥 = 0 −𝑀𝑓𝑧 + 𝑀𝑓𝑧 + 𝑑𝑀𝑓𝑧 + 𝑇𝑦 𝑑𝑥 + 𝑝𝑑𝑥² = 0
2
Qui donne :
𝑑𝑇𝑦 𝑑𝑀𝑓𝑧
=𝑝 et = −𝑇𝑦 5.2
𝑑𝑥 𝑑𝑥
5.3 Essai de flexion (domaine élastique)

Un dispositif représenté ci-dessous permet d'effectuer un essai de flexion
 plane simple sur une
F
poutre reposant sur deux appuis A et B et soumise en C à une force . Un comparateur placé en
D permet de mesurer la flèche lorsque F varie.
y
y
C y
x z
G
D
A B
(S)
F
l
Figure 42 –Essai de flexion plane simple.
Constatations
 La flèche est proportionnelle à l'effort F appliqué et ceci quel que soit le point D choisi.
 Pour une même valeur de F, la flèche est maximum lorsque D est au milieu de la poutre.
 On observe, en effectuant l'essai avec différentes poutres, que la flèche en D est inversement
proportionnelle au moment quadratique IGz de la section.
 Les fibres longitudinales situées au-dessus de la ligne moyenne se raccourcissent et celles
situées en dessous de la ligne moyenne s'allongent.
 Les fibres appartenant au plan (G,x,z) ne changent pas de longueur.

42
 Les allongements et raccourcissement relatifs ( l/l ) sont proportionnels à la distance y de la
fibre considérée au plan (G,x,z).
 Les sections planes normales aux fibres restent planes et normales aux fibres après
déformation.
y
(S')
raccourcissement relatif
y
x
G
plan des fibres neutres allongement relatif
(S)

Figure 43 – Constatations après l'essai de flexion plane simple.

Les constats précédents amènent à deux conséquences.
 – Les fibres s’allongent ou se raccourcissent et sont donc soumises à des contraintes
normales.
Figure 44 – Présence des contraintes normales dues à la flexion plane simple.

 – Entre chaque fibre, on a des variations de longueur qui induisent des contraintes
tangentielles
On a à la fois des :
 contraintes tangentielles longitudinales dans le plan (z;x)
 et des contraintes tangentielles transversales dans le plan (z;y)
Figure 45 – Présence des contraintes tangentielles longitudinales et transversales dues à la

flexion plane simple

43
5.4 Déformation longitudinale et Contrainte normale
On considère un tronçon de poutre infiniment petit, de longueur dx soumis à la flexion plane
simple.
Figure 46 – Tronçon de poutre soumis à la flexion plane simple

Après déformation, le tronçon se transforme en un arc de cercle de longueur dx' et de de rayon de
courbure R. La déformation longitudinale du tronçon dx s'écrit d'après l'équation Erreur !
Source du renvoi introuvable. :
𝑑𝑥 ′ − 𝑑𝑥
𝜀(𝑥) =
𝑑𝑥
Encore l'angle de courbure s'exprime par :
𝑑𝑥 𝑑𝑥′
𝑑𝛼 = =
𝑅 𝑅−𝑦
Donc la courbure s'écrit :
1 𝑑𝛼
= 5.3
𝑅 𝑑𝑥
Ce qui permet d'écrire pour la déformation longitudinale du tronçon :
𝑅−𝑦 1 𝑑𝛼
𝜀(𝑥) = − 1 = −𝑦 = −𝑦 5.4
𝑅 𝑅 𝑑𝑥
Contrainte normale
Figure 47 – Répartition linéaire des contraintes normales dans la section droite

 Répartition linéaire des contraintes normales dans la section droite
44
 Traction / compression de part et d’autre de la ligne moyenne
La contrainte normale devient (équation 1.4 ).
⃗ (𝑀, 𝑥 ) = 𝜎𝑥
𝑇 5.5
Le tenseur des contraintes de Cauchy en M, dans le repère (𝑥, 𝑦, 𝑧) :
𝜎 0 0
𝜎̿(𝑀) = ( 0 0 0) 5.6
0 0 0 (𝑥,𝑦,𝑧)
Loi de Hooke
D'après l'équation de la déformation longitudinale (5.4) :
1
𝜎 = 𝐸. 𝜀 = −𝐸. 𝑦 5.7
𝑅
Relation contrainte normale/moment fléchissant
Le vecteur contrainte dans la section droite s'exprime d'après l'équation (5.7) :
1
𝜎𝑥 = −𝐸. 𝑦 𝑥
𝑅
Le moment résultant du torseur de cohésion :
𝑦 𝐸 𝐸 𝜎
⃗⃗⃗⃗⃗⃗ = ∫ 𝐺𝑀
𝑚𝑓 ⃗⃗⃗⃗⃗⃗ ∧ −𝐸 𝑥 . 𝑑𝑠 = 𝑧 ∫ 𝑦 2 𝑑𝑠 = 𝑧.
⃗⃗ 𝐼𝐺𝑧 = − 𝐼𝐺𝑧 . 𝑧
𝑆 𝑅 𝑅 𝑆 𝑅 𝑦
𝐼𝐺𝑧 = ∫𝑆 𝑦 2 𝑑𝑠 : est le moment quadratique de la section S par rapport à l’axe (G;z).

Finalement
𝑚𝑓 𝑚𝑓𝑀𝐴𝑋
𝜎=− 𝑦 et 𝜎𝑀𝐴𝑋 = 𝑦𝑀𝐴𝑋 5.8
𝐼𝐺𝑧 𝐼𝐺𝑧
Remarques
 La contrainte σ est petite pour les sections ou le moment quadratique IGz est grand.
 La contrainte σ grande pour les points de la section les plus éloignés de la ligne moyenne
Déformation longitudinale
D'après l'expression de la contrainte normale (équation 5.8) et la loi de Hooke (équation Erreur !
Source du renvoi introuvable.)la déformation longitudinale peut être écrite :
𝑚𝑓
𝜀=− 5.9
𝐸𝐼𝐺𝑧
 Les fibres supérieures (y > 0) se raccourcissent, elles sont comprimées ε < 0.

45
 Les fibres inférieures (y < 0) s’allongent, elles sont tendues ε > 0.
 La fibre neutre ne subit pas de déformation y=0.
5.5 Déformations transversales : Équation de la déformée

Reprenons le tronçon de poutre infiniment petit, de longueur dx soumis à la flexion plane simple.
Déformée 𝑦(𝑥)
𝑑𝛼
Figure 48 – Tronçon de poutre soumis à la flexion plane simple

La ligne moyenne se déforme et devient la déformée.
On remarque d'après la Figure 48:
𝑦̇ (𝑥) = 𝑡𝑔(𝛼) ≈ 𝛼(𝑥)
Par suite, il vient d'après les équations (5.3), (5.7) et (5.8) :
𝑑𝛼 𝑚𝑓(𝑥)
𝑦̈ (𝑥) = = Ou encore : 𝑚𝑓 = 𝐸𝐼𝐺𝑧 𝑦̈ 5.10
𝑑𝑥 𝐸𝐼𝐺𝑧
C’est l’équation de la déformée. Il faut alors procéder à deux intégrations successives.

Les constantes d'intégration s'obtiennent grâce aux conditions aux limites aux liaisons.
Appui simple y=0
Encastrement y=0 y' = 0
Tableau 6.Exemples de conditions aux limites :
5.6 Contraintes tangentielles

Mise en évidence expérimentale
On considère deux poutres de sections globales identiques, faites d’un même matériau, soumises
au même chargement. Une des deux poutres est constituée d’un empilement de barres.

46
Figure 49 – Mise en évidence expérimentale des contraintes tangentielles.
Après le chargement, on constate :
 Glissement longitudinal des éléments constituant la poutre composée.
 Poutre monobloc est moins déformée car pas de glissement longitudinal
On déduit :
 existence de forces internes longitudinales,
 donc existence de contraintes tangentielles longitudinales
Théorème de Cauchy
Considérons un petit élément de poutre carré, soumise à la flexion simple, de côté da, situé dans
le plan xOy de la poutre.
Figure 50 Equilibre d'un carré élément de poutre, soumise à la flexion simple..

Ce carré est en équilibre sous l’action :
 des forces volumiques négligées ;
 des résultantes des contraintes agissant sur ses côtés.
La somme des moments /O est nul :
𝑑𝑎 𝜕𝜎𝑥 𝑑𝑎 𝑑𝑎 𝜕𝜎𝑦 𝑑𝑎 𝜕𝜏𝑥𝑦 𝑑𝑎
𝜎𝑥 − (𝜎𝑥 + 𝑑𝑎) − 𝜎𝑦 + (𝜎𝑦 + 𝑑𝑎) + (𝜏𝑥𝑦 + 𝑑𝑎)
2 𝜕𝑥 2 2 𝜕𝑦 2 𝜕𝑥 2
𝜕𝜏𝑦𝑥 𝑑𝑎
− (𝜏𝑦𝑥 + 𝑑𝑎) =0
𝜕𝑥 2
Si l’on néglige les termes en da², cette équation s’écrit :
(𝜏𝑥𝑦 − 𝜏𝑦𝑥 )𝑑𝑎 = 0
Soit encore :
𝜏𝑥𝑦 = 𝜏𝑦𝑥 5.11
C'est le théorème de Cauchy : Les contraintes tangentielles s’exerçant sur deux faces
perpendiculaires sont égales et dirigées vers l’arête commune.

47
Expression de contrainte tangentielle
Considérons un tronçon élémentaire de poutre SS’, de longueur dx et soit le prisme élémentaire
ABCA’B’C’ situé au-dessus du plan ABA’B’.
Le tronçon est en équilibre, le prisme supérieur l’est donc également.
Ecrivons que la résultante des contraintes suivant Gx est nulle :
Figure 51 Equilibre d’un prisme élémentaire de poutre.

τ(y) est la contrainte de cisaillement agissant le long de la corde AB dans la section S, d’ordonnée
y et de longueur b(y).
Bilan des contraintes qui s’exercent sur ABCA’B’C’ suivant Gx :
𝑀 𝑀
Sur S : 𝑅=∬ 𝜎𝑑𝑠 = ∬ 𝑌𝑑𝑠 = 𝑚(𝑦)
𝐴𝐵𝐶 𝐼𝐺𝑧 𝐴𝐵𝐶 𝐼𝐺𝑧
𝑀′ ′ 𝑀′
Sur S : 𝑅′ = ′
𝑚 (𝑦) = 𝑚(𝑦)
𝐼 𝐺𝑧 𝐼𝐺𝑧
Sur ABA’B’ : 𝑅𝑆 = 𝜏(𝑦)𝐴𝑖𝑟𝑒(𝐴𝐵𝐴′ 𝐵 ′ ) = 𝜏(𝑦)𝑏(𝑦)𝑑𝑥
m’(y) et m(y) sont les moments statiques des sections S(A’B’C’) et S(ABC)
Les moments statiques m’(y) = m(y) et les moments quadratiques 𝐼′𝐺𝑧 = 𝐼𝐺𝑧 car les sections
S(A’B’C’)=S(ABC)
La résultante des contraintes suivant Gx est nulle : R-R’+RS = 0 donne
𝑀 𝑀′
𝑚(𝑦) − 𝑚(𝑦) + 𝜏(𝑦)𝑏(𝑦)𝑑𝑥 = 0
𝐼𝐺𝑧 𝐼𝐺𝑧
En fin :
48
𝑇 𝑚(𝑦)
𝜏(𝑦) = 5.12
𝐼𝐺𝑧 𝑏(𝑦)
τ(y) est la contrainte de cisaillement agissant le long de la corde AB de la section S, d’ordonnée y

et de longueur b(y).
T : Effort tranchant dans la section S(x),
m(y) : moment statique de la partie S(y) de la section totale S(x) située au dessus de l’ordonnée y.
𝐼𝐺𝑧 : Moment quadratique de la section totale S(x)
Facteur de forme en cisaillement
Tableau 7. Quelque valeurs de .

On définit le facteur de forme en cisaillement  tel que :
𝑇
𝜏𝑀𝑎𝑥 = 𝜆. 𝜏𝑀𝑜𝑦 = 𝜆. 5.13
𝑆
5.7 Conditions de résistance

𝜎𝑒 𝑚𝑓
Contraintes
𝜎𝑀𝐴𝑋 < 𝜎𝑝 = 𝜎=
𝐼𝐺𝑧 5.14
normales 𝑠 𝑦
Contraintes 𝜏𝑒 𝑇 𝑇. 𝑚(𝑦)
𝜏𝑀𝐴𝑋 < 𝜏𝑝 = 𝜏 = 𝜏𝑀𝑜𝑦 = ou 𝜏= 5.15
tangentielles 𝑠 𝑆 𝐼. 𝑏(𝑦)
Effet prépondérant
En général, une poutre sollicitée à la flexion simple est plus sensible aux contraintes normales
qu’aux contraintes tangentielles. Donc le calcul se fera selon la contrainte normale.

49
5.8 Condition de rigidité
Flèche ou Rotation
5.16
𝑓 ≤ 𝑓𝑙𝑖𝑚𝑖𝑡𝑒 𝛼 ≤ 𝛼𝑙𝑖𝑚𝑖𝑡𝑒

On suppose que la section droite S admette GZ comme axe de symétrie.
En deux points symétriques sur la poutre A et A’ d’ordonnées a et – a, on colle deux jauges
identiques j et j’ dans la direction axiale X’X et on en fait les résistances R1 et R3 d’un demi pont
de Wheatstone. On donne de plus aux résistances étalonnées R2 et R4 des valeurs égales pour
assurer l’équilibre du pont.
Sous l’action du chargement, on a sur S un moment de flexion Mfz et les jauges subissent les
déformations linéaires relatives opposées (équation 5.9) :
𝑚𝑓 𝑚𝑓
𝜀=− 𝑎 et 𝜀′ = 𝑎
𝐸𝐼𝐺𝑧 𝐸𝐼𝐺𝑧
Les déformations se calculent (voir annexe A) :

𝑉0 ou puisque : 2 𝑉0 𝑚𝑓
𝜀1 − 𝜀3 = 4 𝜀1 = = 𝑎 5.17
𝐾𝑉𝑖 𝜀1 = −𝜀3 𝐾 𝑉𝑖 𝐸𝐼𝐺𝑧
5.10 Structure hyperstatique

Exemple de structure hyperstatique
Déterminer les actions en A et B pour le cas suivant (Figure 53):
Le PFS (équation 1.1) appliqué à cette structure permet d'écrire les trois équations suivantes :

50
𝐵𝑥 = 0
𝐴𝑦 + 𝐵𝑦 = 𝑝𝐿
𝐿2
{ 𝑀 𝐵 = 𝑝
2
Le nombre d'équation (p=3) est insuffisant pour déterminer les inconnues de liaisons (n=4) en A
et B : 𝐴𝑦 , 𝐵𝑥 , 𝐵𝑦 et 𝑀𝐵 .
h=n-p=1>0
Le système est hyperstatique d'ordre 1. Il faudra introduire de nouvelles équations issues de la
RDM pour résoudre le problème.
Figure 53 Exemple d'une structure hyperstatique.
h est le degré d'hyperstatisme.
Méthodes de résolution de structures hyperstatiques

 Équations de déformation supplémentaires,
 Principe de superposition des effets de forces
 Formules de Bresse (voir 6.4)
 Méthodes énergétiques (voir 6.5)
Equations de déformation supplémentaires
Etapes à suivre pour résoudre le problème (Exemple Figure 53)

1. Déterminer les actions en A et B,
2. Trouver l’équation du moment de flexion,
3. Etablir l’équation de la déformée y(x)
4. Trouver l’équation manquante en écrivant : y(x=L)=0
5. Déduire les actions en A et B
Principe de superposition des effets de forces
Enoncé
L’étude d’un système hyperstatique se ramène à l’addition, ou la superposition, de deux (ou
plusieurs) systèmes isostatiques 1 et 2.

51
Figure 54 Principe de superposition des effets de forces
Conséquences
F = F1+F2  = 1+2
Charges extérieures Contraintes
q = q1+q2  = 1+2
A = A1+A2
Actions des appuis B = B1+B2 Déformée y(x) = y1(x)+y2(x)
m = m1+m2
Effort tranchant T = T1+T2 Flèche f = f1+f2
Moment de flexion mf = mf1+mf2 Rotation  = 1+2
Méthode de résolution
Utiliser l’équation f = f1+f2 pour trouver les actions des appuis
Etapes à suivre pour résoudre le problème (Exemple Figure 53)

52
Chapitre 6
Sollicitations combinées des structures.
6.1 Repère de Frenet
Figure 55 –Repère de Frenet lié à la ligne moyenne G0G1

Soit une poutre curviligne, définie par sa ligne moyenne G0G1 et soit un repère orthonormé
d’origine G, lié à la ligne moyenne G0G1.
⃗ ,𝑏⃗} est le repère de Frenet de la courbe𝑡 G0G1 au point G :

Le repère orthonormé direct {G, 𝑡, 𝑛
 𝑡 premier vecteur de la base : unitaire, tangent à la courbe G0G1 en G et orienté dans le
⃗⃗⃗⃗⃗⃗⃗⃗⃗
𝑑𝑂𝐺
sens des abscisses curvilignes croissantes : 𝑡 =
𝑑𝑠
 𝑛⃗ second vecteur de la base : unitaire, normal à 𝑡, dirigé vers le centre de courbure, nommé
vecteur normal principal,.
 ⃗ = 𝑡 ∧ 𝑛⃗, se nomme binormal.

𝑏⃗ troisième vecteur de la base défini par : 𝑏

53
Torseur des efforts intérieurs exprimé dans le repère de Frenet
Figure 56 –Torseur des efforts intérieurs exprimé dans le repère de Frenet

𝑅⃗ = 𝑁𝑡 + 𝑇𝑛 𝑛⃗ + 𝑇𝑏 𝑏⃗
Cohésion  E  E 2 E  E1= { 6.1
⃗⃗ 𝐺 = 𝑚𝑡 𝑡 + 𝑚𝑓𝑛 𝑛⃗ + 𝑚𝑓𝑏 𝑏⃗
𝑚
6.2 Equilibre local

Considérons une tranche élémentaire de poutre de longueur ds, comprise entre les sections droites
voisines S’(d’abscisses curvilignes s et de centres G’) et S (d’abscisses curvilignes s + ds et de
centres G).
Supposons que les forces et couples extérieurs appliqués au niveau de cette tranche soient tous
répartis. Notons p
⃗ (𝑠). ds et ⃗⃗⃗⃗⃗
mr (𝑠). ds leur résultantes en G.
Figure 57 –Tranche élémentaire de poutre de longueur ds.
Etude de l’équilibre du tronçon

L'équilibre du tronçon permet d'écrire les deux équations de l'équilibre local :
𝑑𝑅⃗(𝑠)
Équation de la résultante + 𝑝(𝑠) = ⃗0
𝑑𝑠 6.2
𝑑𝑚
⃗⃗ (𝑠)
Équation du moment en G + 𝑡 ∧ 𝑅⃗(𝑠) + 𝑚 ⃗
⃗⃗ 𝑟 (𝑠) = 0
𝑑𝑠
54
6.3 Technique des fonctions de singularité
C’est une technique basée sur les fonctions de singularités ou à domaines limité, utilisée pour
déterminer la déformée des poutres en flexion.
L'utilisation de cette technique est très pratique lorsqu’on a plus d’une zone d’étude dans la
même poutre et ceci pour le cas de charges concentrées et réparties.
Macaulay a appliqué cette technique en 1919 [1] pour effectuer la méthode de la double
intégration et donc portera après son nom. Mais en réalité, cette méthode a été développé
auparavant par Clebsch en 1862 [2]
[1] W. H. Macaulay, "A note on the deflection of beams", Messenger of Mathematics, 48 (1919),
129.
[2] J. T. Weissenburger, ‘Integration of discontinuous expressions arising in beam theory’, AIAA
Journal, 2(1) (1964), 106–108.
Principales fonctions de singularité
Figure 58 –Principales fonctions de singularité :
Propriétés
 pour n0
0 𝑠𝑖 𝑥 < 𝑎
⟨𝑥 − 𝑎⟩𝑛 = { ∫ 𝑃⟨𝑥 − 𝑎⟩𝑛 𝑑𝑥 = 𝑃
⟨𝑥 − 𝑎⟩𝑛+1
(𝑥 − 𝑎)𝑛 𝑠𝑖 𝑥 > 𝑎 et + 𝐶𝑚 6.3
𝑛+1
 pour n=-1
⟨𝑥 − 𝑎⟩−1 = 0 0 𝑠𝑖 𝑥 < 𝑎
et ∫ 𝑃⟨𝑥 − 𝑎⟩−1 𝑑𝑥 = 𝑃⟨𝑥 − 𝑎⟩0 + 𝐶𝑚 = 𝐶𝑚 + { 6.4
1 𝑠𝑖 𝑥 > 𝑎
 pour n=-2
⟨𝑥 − 𝑎⟩−2 = 0 et ∫ 𝑃⟨𝑥 − 𝑎⟩−2 𝑑𝑥 = 𝑃⟨𝑥 − 𝑎⟩−1 + 𝐶𝑚 = 𝐶𝑚 + 0 6.5

55
6.4 Formules de Bresse
6.4.1 Utilité
Les Formules de Bresse permettent de résoudre les déplacements et les actions de liaisons pour
des poutres dont la ligne moyenne n'est pas rectiligne ou dont les sollicitations sont simples ou
composées et pour des problèmes 2D ou 3D.
6.4.2 Torseur de petits déplacements

Hypothèses cinématiques
 Les dimensions transversales d’une poutre sont faibles devant les dimensions longitudinales,
 Chaque section droite S(s) d’abscisse s se comporte comme un solide,
Figure 59 –Tranche élémentaire de poutre de longueur ds.

Pour chaque section, on a donc :
𝜔⃗ (𝑆(𝑠))
{𝑑𝑒𝑝𝑆/𝑅 } = { } et 𝑢
⃗ (𝑀 ∈ 𝑆(𝑠)) = 𝑢
⃗ (𝐺 ∈ 𝑆(𝑠)) + 𝜔 ⃗⃗⃗⃗⃗⃗
⃗ (𝑆(𝑠))Λ𝐺𝑀 6.6
𝐺 𝑢
⃗ (𝐺 ∈ 𝑆(𝑠)) 𝐺
Avec G : centre de la surface S(s) et M : point de la surface S(s).
Exemples
E1
y
(S )
R = N .x
⃗0
{𝑑𝑒𝑝𝑆/𝑅 } = { } F R
x
⃗ (𝐺) = 𝑢. 𝑥 𝐺
𝑢
G
Traction 𝐺
⃗ (𝑀) = 𝑢
𝑢 ⃗ (𝐺) = 𝑢. 𝑥 M G1
G1 z G2
fi g .8
M G2
𝜔
⃗ (x) = 𝜔𝑥 R y
(S)
{𝑑𝑒𝑝𝑆/𝑅 } = { } avec ∶ 𝜔 =∝= 𝜃. 𝑥 M G1
𝐺 ⃗0
MG
x
G1
𝐺 G
Torsion 𝑢 ⃗ (𝑥))⃗⃗⃗⃗⃗⃗
⃗ (𝑀) = 𝜔 𝐺𝑀 = 𝜔𝑥 (𝑦𝑦 + 𝑧𝑧) z
= 𝜔(𝑦𝑧 − 𝑧𝑦) = 𝜃. 𝑥. (𝑦𝑧 − 𝑧𝑦)

56
⃗⃗⃗⃗⃗⃗⃗⃗⃗⃗ = 𝜔𝑧
𝜔(𝑥)
{𝑑𝑒𝑝𝑆/𝑅 } = { }
𝐺 ⃗ (𝐺) = 𝑢. 𝑦 𝐺
𝑢
Flexion ⃗ (𝑀) = 𝑢
𝑢 ⃗ (𝐺) + 𝜔 ⃗ (𝑥))⃗⃗⃗⃗⃗⃗
𝐺𝑀
= 𝑢. 𝑦 + 𝜔𝑧(𝑦𝑦 + 𝑧𝑧)
= 𝑢. 𝑦 + 𝜔. 𝑦. 𝑥
6.4.3 Torseur de petites déformations

Prenons les deux sections voisines S et S(s+ds) de la poutre Figure 59 :
Pour chaque section, on peut écrire :
𝜔⃗ (𝑠) 𝜔⃗ (𝑠 + 𝑑𝑠) = 𝜔⃗ (𝑠) + 𝑑𝜔⃗ (𝑠)

{𝑑𝑒𝑝𝑆/𝑅 }𝐺 = { } et {𝑑𝑒𝑝𝑆′/𝑅 } = { }
𝑢
⃗ (𝑠) 𝐺 𝐺′ ⃗ (𝑠 + 𝑑𝑠) = 𝑢
𝑢 ⃗ (𝑠) + 𝑑𝑢
⃗ (𝑠) 𝐺′
Le torseur de taux de déformation :
1 𝑑 𝛾 (𝑠)
{𝐷𝑒𝑓𝑆/𝑅 }𝐺 = ({𝑑𝑒𝑝𝑆′/𝑅 }𝐺 − {𝑑𝑒𝑝𝑆/𝑅 }𝐺 ) = {𝑑𝑒𝑝𝑆/𝑅 }𝐺 = { }
𝑑𝑠 𝑑𝑠 𝜀 (𝑠) 𝐺
𝑑𝜔
⃗ (𝑠)
𝛾 (𝑠) =
{𝐷𝑒𝑓𝑆/𝑅 } = 𝑑𝑠 6.7
𝐺 𝑑𝑢
⃗ (𝑠)
{𝜀 (𝑠) = 𝑑𝑠 + 𝑡 𝜔 ⃗}
𝐺
6.4.4 Lois de comportement

Ce sont les relations qui lient les contraintes aux déformations :
1 1
𝑁 𝑇𝑛 𝑇𝑏 0 0 0 0
𝐺𝐼0
𝜀 (𝑠) = [𝐴]. 𝑅⃗ (𝑠) = 𝑡+ 𝑛⃗ + 𝑏⃗ 𝐸𝑆
1
𝐸𝑆 𝐺𝑆 𝐺𝑆 [𝐴] = 0
1
𝑚𝑡 𝑚𝑓𝑛 𝑚𝑓𝑏 0 , [𝐵] = 0 0 6.8
𝐺𝑆 𝐸𝐼𝑛
𝛾 (𝑠) = [𝐵]. 𝑚
⃗⃗ (𝑠) = 𝑡+ 𝑛⃗ + 𝑏⃗ 1 1
𝐺𝐼0 𝐸𝐼𝑛 𝐸𝐼𝑏
[0 0
𝐺𝑆] [
0 0
𝐸𝐼𝑏 ]
6.4.5 Formules de Bresse

Soit une poutre curviligne, définie par sa ligne moyenne G1G2 et soit le repère orthonormé
d’origine G, lié à la ligne moyenne G1G2.
Si on intégrons les équations 6.7 entre les points G1 et G2, on retrouve les deux formules de
Bresse :
𝐺2
𝜔 ⃗ (𝐺1) + ∫ 𝛾 (𝑠) 𝑑𝑠
⃗ (𝐺2) = 𝜔 1FB
𝐺1
6.9
𝐺2 𝐺2
𝑢
⃗ (𝐺2) = 𝑢
⃗ (𝐺1) + 𝜔 ⃗⃗⃗⃗⃗⃗⃗⃗⃗
⃗ (𝐺1) 𝐺 ⃗⃗⃗⃗⃗⃗⃗
1 𝐺2 + ∫ 𝜀 (𝑠) 𝑑𝑆 + ∫ 𝛾 (𝑠) 𝐺𝐺2 𝑑𝑆 2FB
𝐺1 𝐺1

57
Figure 60 – Poutre définie par sa ligne moyenne G1G2 et le repère de Frenet défini en G.
6.5 Théorèmes énergétiques
Figure 61 –Déformation élastique d’une poutre
Par définition :
Énergie de déformation :
Wd = - Wint = Wext 6.10
Wd = Énergie de déformation.
Wint = Travail des efforts intérieurs de cohésion.
Wext : Travail des efforts extérieurs.
6.5.1 Travail des efforts extérieurs : première formule de Clapeyron

Considérons un corps élastique soumis à un système de n forces Pi agissant aux points Ai leur
imposant des déplacements i

58
Figure 62 –Corps élastique soumis à un système de n forces Pi
Le travail élémentaire fournit par chaque force Pi : 𝑑𝑊𝑖 = 𝑃()𝑑
Le travail total fournit par chaque force Pi est l'air total du triangle :
1
𝑊𝑖 = 𝑃𝛿 6.11
2 𝑖 𝑖
Le travail total fournit par l’ensemble des forces Pi :
1 𝑛
𝑊𝑒𝑥𝑡 = ∑ 𝑃𝑖 𝛿𝑖 6.12
2 𝑖=1
C’est la première formule de Clapeyron
Exemple1 : Sollicitation de traction
Hypothèses :
 Effort de traction F variable de 0 à F
 Proportionnalité entre l’effort et

l’allongement
Le travail de l’effort de traction est l'aire du

1
triangle OAB : 𝑊𝑒𝑥𝑡 = 2 𝐹. 𝑥
Figure 63 –Poutre sollicitée à la traction.
Exemple2 : Sollicitation quelconque

1 1
𝑊𝑒𝑥𝑡 = 𝐹𝑖 𝑈𝑖 . + 𝐶𝑖 𝑖
2 2
(𝑈𝑖 , 𝜃𝑖 ) sont les déplacements et rotation dus à
(𝐹𝑖 , 𝐶𝑖 ).
Figure 64 –Poutre avec une sollicitation quelconque.

59
6.5.2 Énergie de déformation : forme intégrale
Cas de la traction
Etudions l'équilibre d’un tronçon de poutre de longueur dx sollicité à la traction :
F=N
Figure 65 –Équilibre d’un tronçon de longueur dx sollicité à la traction
Soit ∆(dx) l’allongement du tronçon dx.

D'après la loi de Hooke (2.4), on peut écrire :
𝑁. 𝑑𝑥
𝑑𝑢 =
𝐸𝑆
L'énergie de déformation élémentaire du tronçon dx est donc d'après l'équation (6.11) :
1 1 𝑁2
𝑑𝑊𝑑 = 𝑁. 𝑑𝑢 = . 𝑑𝑥
2 2 𝐸𝑆
L'énergie de déformation totale s'obtient en intégrant sur toute la poutre :
1 𝐿 𝑁2 1 𝐿 1 𝐿
𝑊𝑑 = ∫ . 𝑑𝑥 = ∫ 𝜎𝑖𝑗 . 𝜀𝑖𝑗 . 𝑆𝑑𝑥 = ∫ {𝜏𝑐 }𝐺 . {𝐷𝑒𝑓}𝐺 𝑑𝑥
2 0 𝐸𝑆 2 0 2 0
Cas général :
1 𝐿 1 𝐿
𝑊𝑑 = ∫ 𝜎𝑖𝑗 . 𝜀𝑖𝑗 . 𝑆𝑑𝑥 = ∫ {𝜏𝑐 }𝐺 . {𝐷𝑒𝑓}𝐺 𝑑𝑥 6.13
2 0 2 0
Le torseur de cohésion (équation 6.1Erreur ! Source du renvoi introuvable.) et le torseur des
petites déformations (6.8) permettent de réécrire l'équation (6.13) :
1 𝑁2 1 𝑇2 1 𝑀𝑡 2 1 𝑀𝑓𝑧 2
𝑊𝑑 = ∫ . 𝑑𝑥 + ∫ . 𝑑𝑥 + ∫ . 𝑑𝑥 + ∫ . 𝑑𝑥 6.14
2 𝑠𝑡𝑟 𝐸𝑆 2 𝑠𝑡𝑟 𝐺𝑆 2 𝑠𝑡𝑟 𝐺𝐼0 2 𝑠𝑡𝑟 𝐸𝐼𝐺𝑧
Avec :
Sollicitation simple Energie de déformation

1 𝑁2
N : Traction/ compression 𝑊𝑑 = ∫ . 𝑑𝑥
2 𝑠𝑡𝑟 𝐸𝑆
60
1 𝑇2
Ty ou Tz : Cisaillement 𝑊𝑑 = ∫ . 𝑑𝑥
2 𝑠𝑡𝑟 𝐺𝑆
1 𝑀𝑡 2
Mt : Torsion 𝑊𝑑 = ∫ . 𝑑𝑥
2 𝑠𝑡𝑟 𝐺𝐼0
1 𝑀𝑓𝑧 2
Mfy ou Mfz : Flexion 𝑊𝑑 = ∫ . 𝑑𝑥
2 𝑠𝑡𝑟 𝐸𝐼𝐺𝑧
Tableau 8. Energie de déformation pour les sollicitations simples.
6.5.3 Egalité de Clapeyron

L’égalité (6.10) entre énergie de déformation et travaux des efforts extérieurs permet d’écrire
enfin : l'égalité de Clapeyron :
1 𝑁2 1 𝑇2 1 𝑀𝑡 2 1 𝑀𝑓𝑧 2 1 𝑛
∫ . 𝑑𝑥 + ∫ . 𝑑𝑥 + ∫ . 𝑑𝑥 + ∫ . 𝑑𝑥 = ∑ 𝑃𝑖 𝛿𝑖 6.15
2 𝑠𝑡𝑟 𝐸𝑆 2 𝑠𝑡𝑟 𝐺𝑆 2 𝑠𝑡𝑟 𝐺𝐼0 2 𝑠𝑡𝑟 𝐸𝐼𝐺𝑧 2 𝑖=1
Limite de l'égalité de Clapeyron

Cette égalité ne peut s’appliquer que si les conditions suivantes sont toutes satisfaites :
 Système isostatique,
 Appuis fixes,
 Une seule force extérieure localisée est appliquée,
 Le déplacement cherché est celui du point d’application de la seule force appliquée.
Cherchons une autre forme pour la formule de Clapeyron qui pourra être appliquée pour toutes
les conditions.
6.5.4 Théorèmes de réciprocité de Betti- Rayleigh et de Maxwell- Betti

Réciprocité des travaux
Considérons un corps élastique soumis à un système de n forces Pi (Figure 66- a) :
Figure 66 –Corps élastique soumis à un système de n forces Pi puis un système de m forces Qj

61
1
L’énergie élastique correspondante au système Pi: 𝑊(𝑃𝑖 ) = ∑𝑛𝑖=1 𝛿𝑖 𝑃𝑖
2
On applique ensuite un système de m forces Qj. (Figure 66- b)

1
L’énergie élastique correspondante au système Qj : 𝑊(𝑄𝑗 ) = ∑𝑚
𝑗=1 𝑗 𝑄𝑗
2
Le système Qj provoque aussi des nouveaux déplacements aux points Ai. L’énergie élastique
correspondante : 𝑊′(𝑃𝑖 ) = ∑𝑛𝑖=1 𝛿′𝑖 𝑃𝑖
1 1
L’énergie de déformation totale : 𝑊(𝑃𝑖 , 𝑄𝑗 ) = ∑𝑛𝑖=1 𝛿𝑖 𝑃𝑖 + ∑𝑚 𝑛
𝑗=1 𝑗 𝑄𝑗 + ∑𝑖=1 𝛿′𝑖 𝑃𝑖
2 2
Procédons cette fois de manière inverse en appliquant le système Qj et ensuite le système Pi
(Figure 67) :
1
L’énergie élastique correspondante au système Qj: 𝑊(𝑄𝑗 ) = ∑𝑚
𝑗=1 𝑗 𝑄𝑗
2
On applique ensuite le système Pi.
1
L’énergie élastique correspondante : 𝑊(𝑃𝑖 ) = ∑𝑛𝑖=1 𝛿𝑖 𝑃𝑖
2
Le système Qj provoque aussi des nouveaux déplacements aux points Ai. L’énergie élastique
correspondante : 𝑊′(𝑄𝑗 ) = ∑𝑚
𝑗=1 ′𝑗 𝑄𝑗
1 1
L’énergie de déformation totale: 𝑊(𝑄𝑗 , 𝑃𝑖 ) = ∑𝑚
𝑗=1 𝑗 𝑄𝑗 + ∑𝑛𝑖=1 𝛿𝑖 𝑃𝑖 + ∑𝑚
𝑗=1 ′𝑗 𝑄𝑗
2 2
Figure 67 –Corps élastique soumis à un système de m forces Qj puis un système de n forces Pi
Théorèmes de réciprocité
1. Les deux expressions de l’énergie de déformation totale sont égales puisque l’énergie ne
dépend que de l’état final du système. On déduit donc :
′
∑𝑛𝑖=1 𝛿 ′ 𝑖 𝑃𝑖 = ∑𝑚
𝑗=1  𝑗 𝑄𝑗
C’est le théorème de réciprocité de Betti- Rayleigh

Le travail du système de forces Pi lors de la déformation imposée par le système Qj est égale au
travail du système de forces Qj lors de la déformation imposée par le système Pi
2. Le théorème de réciprocité de Betti- Rayleigh peut s’exprimer aussi par l’égalité des
coefficients d’influence réciproques : 𝑎𝑖𝑗 = 𝑎𝑗𝑖

62
3. Du théorème de Betti-Rayleigh on peut déduire un cas particulier :
Le déplacement du point A provoqué par une force appliquée au point B est égal au déplacement
du point B provoqué par la même force appliquée au point A.
C’est le théorème de Maxwell ou théorème de réciprocité de Maxwell- Betti
Exemple :
𝑃
𝑦𝐴 = (2𝑙 + 𝑎)(𝑙 − 𝑎)²
6𝐸𝐼
𝑃
𝑦𝐵 = (2𝑙 + 𝑎)(𝑙 − 𝑎)²
6𝐸𝐼
Le déplacement du point A provoqué par P au point B de la poutre 1, est égal au déplacement du

point B provoqué par P (la même force ) appliquée au point A de la poutre 2.
C’est le théorème de Maxwell ou théorème de réciprocité de Maxwell- Betti
6.5.5 Travail des efforts extérieurs : seconde formule de Clapeyron

Considérons un corps élastique subissant un système de n forces 𝑃𝑖 {𝑃1 , 𝑃2 , … , 𝑃𝑛 } appliqués aux
points 𝐴𝑖 {𝐴1 , 𝐴2 , … , 𝐴𝑛 } causant les déplacements 𝛿𝑖 {𝛿1 , 𝛿2 , … , 𝛿𝑛 }.
Figure 68 –Corps élastique soumis à un système de n forces Pi
 Lorsque la force 𝑃1 est appliquée en 𝐴1 , le déplacement en 𝐴1 est

𝛿1 = 𝛿11
 Si en plus de la force 𝑃1 en 𝐴1 , on applique la force 𝑃2 en 𝐴2 ,
Le déplacement en 𝐴1 est :
𝛿1 = 𝛿11 + 𝛿12
Et le déplacement en 𝐴2 est :
𝛿2 = 𝛿21 + 𝛿22
 Si l’ensemble de toutes les forces 𝑃𝑖 est appliquée, le déplacement de tout point Ai s’écrit :
63
𝑛
𝛿𝑖 = 𝛿𝑖1 + 𝛿𝑖2 + ⋯ + 𝛿𝑖𝑗 + ⋯ + 𝛿𝑖𝑛 = ∑ 𝛿𝑖𝑗
𝑗=1
𝛿𝑖𝑗 est la contribution de la force 𝑃𝑗 au déplacement du point 𝐴𝑖 dans la direction de 𝑃𝑖
Figure 69 –Contribution de la force Pj au déplacement du point Ai dans la direction de Pi
La proportionnalité entre forces et déplacements permet d’écrire :

𝑛
𝛿𝑖 = 𝑎𝑖1 𝑃1 + 𝑎𝑖2 𝑃2 + ⋯ + 𝑎𝑖𝑗 𝑃𝑗 + ⋯ + 𝑎𝑖𝑛 𝑃𝑛 = ∑ 𝑎𝑖𝑗 𝑃𝑗
𝑗=1
𝑎𝑖𝑗 est le coefficient d’influence : il représente la projection sur la direction 𝑃𝑖 appliquée au point
𝐴𝑖 du déplacement provoqué en ce point par une force unité appliquée au point 𝐴𝑗 dans la
direction de 𝑃𝑗
D’après la première formule de Clapeyron (6.12), on peut écrire la seconde forme de la formule
de Clapeyron :
1 𝑛 𝑛
𝑊𝑒𝑥𝑡 = ∑ ∑ 𝑎𝑖𝑗 𝑃𝑖 𝑃𝑗 6.16
2 𝑖=1 𝑗=1
6.5.6 Théorème de Castigliano

𝜕𝑊𝑑
𝛿𝑘 = 6.17
𝜕𝑃𝑘
Dans un corps élastique, le déplacement k d’une force Pk est égal à la dérivée partielle de
l’énergie de déformation par rapport à cette force.
Pour démontrer le théorème, reprenons la deuxième formule de Clapeyron de l’énergie :
1 𝑛 𝑛
𝑊𝑒𝑥𝑡 = ∑ ∑ 𝑎𝑖𝑗 𝑃𝑖 𝑃𝑗
2 𝑖=1 𝑗=1
Puis isolons les termes contenant Pk :

1 𝑛 1 𝑛 1
𝑊𝑒𝑥𝑡 = ∑𝑖=1 𝑎𝑖𝑘 𝑃𝑖 𝑃𝑘 + ∑𝑗=1 𝑎𝑘𝑗 𝑃𝑘 𝑃𝑗 + 𝑎𝑘𝑘 𝑃²𝑘 + 𝑊′
2 2 2
𝑖≠𝑘 𝑗≠𝑘
W’ : Terme Indépendant de 𝑃𝑘 .
Dérivons l'énergie de déformation par rapport à 𝑃𝑘 :

64
𝜕𝑊𝑑 𝜕𝑊𝑒𝑥𝑡 1 𝑛 1 𝑛
= = ∑ 𝑎𝑖𝑘 𝑃𝑖 + ∑𝑗=1 𝑎𝑘𝑗 𝑃𝑗 + 𝑎𝑘𝑘 𝑃𝑘
𝜕𝑃𝑘 𝜕𝑃𝑘 2 𝑖=1 2
𝑖≠𝑘 𝑗≠𝑘
1 𝑛 1 𝑛
= ∑𝑖=1 𝑎𝑖𝑘 𝑃𝑖 + ∑𝑗=1 𝑎𝑘𝑗 𝑃𝑗
2 2
1 𝑛 1 𝑛
= ∑𝑖=1 𝑎𝑘𝑖 𝑃𝑖 + ∑𝑗=1 𝑎𝑘𝑗 𝑃𝑗
2 2
Car 𝑎𝑖𝑘 = 𝑎𝑘𝑖 . D’où :

𝜕𝑊𝑑 1 1
= 𝛿𝑘 + 𝛿𝑘 = 𝛿𝑘
𝜕𝑃𝑘 2 2
Théorème de Castigliano : Intégrales de Mohr au cas où seul le moment de flexion est

considéré
∂Wd 1 s2 ∂mf 1 s2
δk = = ∫ mf ds = ∫ mf m ̅ f ds 6.18
∂Pk EI str ∂Pk EI str
mf : Moment de flexion réelle de la structure
m
̅ f : Moment de flexion de la structure ou PK =1 est la seule force exercée.
Souvent les répartitions de mf sont de forme constante, linéaire (triangle ou trapèze) ou
parabolique : on utilise le tableau des intégrales de Mohr pour ce genre de calcul (Tableau 9).
Théorème de Castigliano pour un point non chargé : Application d’une charge fictive
Exemple : Poutre sur 2 appuis
Flèche en G ?
Figure 70 –Poutre sur 2 appuis et charge fictive au point G.
On applique au point G, une charge fictive Q travaillant dans le déplacement yG . On applique

ensuite le théorème de Castigliano (6.17) tout en rendant Q=0 au résultat définitif.
𝜕𝑊𝑑
𝐺 = [ ] 6.19
𝜕𝑄 𝑄=0
Théorème de Castigliano pour un point non chargé : Intégrales de Mohr au cas où seul le
moment de flexion est considéré
1 𝑠2
𝛿𝐺 = ∫ 𝑚𝑓 𝑚 ̅ 𝑓 𝑑𝑠 6.20
𝐸𝐼 𝑠1
mf : Moment de flexion réelle de la structure (sans Q)
𝑚
̅ 𝑓 : Moment de flexion de la structure ou Q=1 est la seule force exercée

65
Le tableau des intégrales de Mohr peut être utilisé aussi pour ce genre de calcul (Tableau 9).
6.5.7 Théorème de Ménabréa : Pour les structures hyperstatiques

Théorème :
La dérivée partielle de l’énergie de déformation par rapport à chacune des inconnues
surabondantes Ii est nulle à condition que les points d’application des forces ne bougent pas (Ui =
0) ou que les sections ne tournent pas (qi = 0)
Pour une structure hyperstatique, suivre les 4 étapes ci-dessous :
Etape1 : Fixer les n inconnues de liaisons et les p équations d'équilibre,
Etape2 : n - p = k est le degré d'hyperstaticité. Il faut k équations à ajouter.
Etape3 : On fixe k inconnues surabondantes (𝐈𝟏 , 𝐈𝟐 , … , 𝐈𝐤 ) puis on écrit les inconnues restantes
(𝐈𝐤+𝟏 , 𝐈𝐤+𝟐 , … , 𝐈𝐧 ) et l'énergie de déformation (𝐖𝐝 ) en fonction des k inconnues surabondantes :
Ik+1 , Ik+2 , … , In = 𝑓(I1 , I2 , … , Ik )
{
Wd = Wd (I1 , I2 , … , Ik )
Etape4 : Les k équations à ajouter sont :
∂Wd
[ ] =0 6.21
∂Ii 𝑖=1 ..𝑘

66
Tableau 9. Intégrales de Mohr.

67
Chapitre 7
Flambage
7.1 Essai de flambage
Figure 71 – Poutre bi-articulée en essai de flambage.
Hypothèses
 Poutre parfaitement rectiligne, homogène et isotrope.
 Les actions aux points A et B sont parfaites.
 La charge P est parfaitement axiale et croit progressivement.
Observations
 On observe que la poutre reste rectiligne dans un premier temps : elle est simplement
comprimée.
 Lorsque la charge P croit et dépasse une certaine valeur Pc appelée charge critique, la poutre
se déforme en fléchissant dans une direction perpendiculaire à l'axe de compression.
Ce phénomène s’appelle le flambement ou flambage.

68
Causes du flambage
Dans la pratique, les hypothèses énoncées sont irréalisables et la charge P provoque vite un
moment de flexion.
7.2 Charge critique d’Euler
Soit la poutre AB sous l’influence de l’effort Normal P en A

et B. La ligne moyenne se déforme et prend une très légère
courbure.
y(x) est la déformée au point G.
Cette déformation influe sur l’équilibre statique et provoque

un moment de flexion 𝑀𝑧 .
Avec :
𝑀𝑧 (𝑥) = −𝑃 𝑦(𝑥)
Or d'après les résultats de la RDM avec l'hypothèse des

petites déformations :
𝑀𝑧 (𝑥) = 𝐸. 𝐼. 𝑦′′(x) Figure 72 – poutre en flambage.
Ce qui donne : 𝑃
𝑦 ′′ + 𝑦=0 7.1
𝐸𝐼
La solution générale de l'équation 7.1 est de la forme :
𝑦 = 𝐴 cos(𝜔𝑥) + B sin(𝜔𝑥)
Avec : P
ω2 = 7.2
E𝐼
A et B sont des constantes déterminées par les conditions aux limites :
 Pour x = 0, y(0) = 0 ; donc A = 0
 Pour x = L, y(L) = 0 ; donc B sin ωL = 0 : B  0, sin ω L = 0
Ce qui donne : ωL = n π Avec n entier 7.3
1er mode de flambage (n = 1 )

D'après l'équation 7.3 pour n = 1 :
π
ω=
L
Et la déformée : π
𝑦 = B sin 𝑥 7.4
L
69
N.B : La constante B reste indéfini. Elle représente la flèche maximale lorsque x=L/2.
D'après les équations 7.2 et 7.3 :
P π²
ω2 = =
E I L²
La charge critique d’Euler (n = 1 )
π²E I
Pc =
L²
Autres modes de flambement (n = 2, 3, …,n )
π²E. I π²E. I
Pc(2) = 4 = 4Pc1 Pc(n) = n² = n²Pc1
L² L²
N.B : Seule la charge critique Pc(1) est considérée, en effet la poutre sera déjà ruinée avant que la
charge n'atteigne Pc(2) ou Pc(n) .
Remarques
 I est le moment quadratique le plus faible. Pour les deux sections ci- contre, le flambement se
produit autour de l’axe Gz (I𝐆𝐳 I𝐆𝐲).
z z
y y
G G
 La section du triangle équilatéral résiste le mieux au flambement en comparaison avec

d’autre forme de section de même aire.
 L'équation (5) suppose que la poutre est dans le domaine élastique : il faut donc vérifier que :
Pc π²E I
𝜎c = ≤ 𝜎ec ou encore : ≤ 𝜎ec
𝑆 S L²
𝐿 𝐸 𝐸
Soit : ≥ 𝜋√ ou encore : 𝜆 ≥ 𝜋√
𝜎ec
√𝐼 𝜎ec
𝑆
Avec :
𝐿 : Elancement (sans unité)
𝜆=
𝜌
𝐼 : Rayon de giration
𝜌=√
𝑆

70
𝜎ec : Limite élastique à la compression
7.3 Cas possibles pour la poutre

Plusieurs cas sont possibles pour la poutre :
 P < Pc : compression simple, la poutre reste droite, elle est en équilibre stable.
 P = Pc : la poutre peut rester droite ou fléchir (flamber) avec une flèche égale à B, elle est en
équilibre neutre.
A noter que B = ymaxi est en général petit.
 P > Pc : forte tendance au flambement de la poutre, elle est en équilibre instable.
B augmentera très rapidement avec un léger accroissement de N.
7.4 Influence des liaisons aux appuis

La formule donnant Pc n’est valable que si la poutre est bi– articulée et sans frottement. Les
résultats établis pour la poutre bi-articulée peuvent être généralisés pour des poutres dont les
conditions d’appuis sont différentes (Tableau 10). La charge critique d’Euler sera :
𝜋²𝐸. 𝐼
𝑃𝑐 = 7.5
Lf ²
Lf : longueur libre de flambage.
longueur libre et
Cas de poutre Poutre équivalente
charge critique
𝐿𝑓 = 2𝐿
Encastrée-
libre 1 𝜋²𝐸. 𝐼
𝑃𝑐 =
4 𝐿²
𝐿
𝐿𝑓 =
Encastrée- 2
encastrée 𝜋²𝐸. 𝐼
𝑃𝑐 = 4
𝐿²
√2
𝐿𝑓 = 𝐿 = 0.7 𝐿
Encastrée- 2
articulée 𝜋²𝐸. 𝐼
𝑃𝑐 = 2
𝐿²
Tableau 10. Longueur libre et charge critique pour différents cas d’appuis.
71
7.5 Calcul de résistance (Méthode d’Euler)
7.5.1 Élancement :
Lf
𝜆= 7.6
𝜌
I : moment quadratique le plus faible,
S : aire de la section droite.
7.5.2 Coefficient de sécurité :

Pc 𝜎e
𝑘= = 2𝑠 = 2
𝑃 𝜎pc
𝜎e : contrainte de limite élastique,
𝜎pc : contrainte pratique de compression
s : coefficient de sécurité de contraintes.
Donc :
Pc Pc 𝜎pc
𝑃= = = P 7.7
𝑘 2𝑠 2𝜎e c
7.5.3 Élancement critique 𝛌𝐜

 Lorsque Mf = 0, y’’ = 0 : seule la contrainte de compression est à considérer :
𝑃
𝜎comp = 7.8
𝑆
 Lorsque Mf est Maxi : on vérifie expérimentalement que la contrainte normale est le double
de la contrainte de compression et doit alors être égale à la contrainte pratique de
compression 𝜎pc :
𝑃
2= 𝜎pc 7.9
𝑆
Or d'après les équations 7.7 et 7.9, on peut écrire :
Pc 𝜎pc 𝜎pc
𝑃= = Pc = 𝑆
𝑘 2𝜎e 2
Donc :
Pc = 𝜎e . 𝑆 7.10
Or d'après les équations 7.5 et 7.10, on peut écrire :
𝜋 2 E. 𝐼 D’où : Lf ². 𝑆 Lf ² 𝜋²E
𝑃𝑐 = = 𝜎e . 𝑆 = = 𝜆² =
L2f 𝐼 𝜌² 𝜎e
À partir de cette valeur de l’élancement appelée élancement critique, la poutre doit être calculée
au flambage.

72
𝜋²E
𝜆c ² = 7.11
𝜎e
Remarques
L'élancement critique ne dépend que des caractéristiques du matériau.
 Pour acier A33 : E=2. 105 Mpa, 𝜎e =155 Mpa, 𝜆c =113
 Pour acier 45SCD6 : E=2. 105 Mpa, 𝜎e =1 370 Mpa, 𝜆c =38
Valeurs moyennes adoptées

Configuration Profilés en acier Poteaux en fonte Poteaux en bois
𝝀𝐜 100 60 70
7.5.4 Condition de résistance au flambage :

Or d'après les équations 7.5 et 7.7 on peut écrire :
𝜎pc π²E I
𝑃=
2𝜎e 𝐿f ²
Et sachant les équations 7.11 et 7.6 , la force admissible d'Euler peut s'écrire :
𝑆. 𝜎pc
𝑃=
𝜆 2 7.12
2( )
𝜆c
𝑆. 𝜎pc : Force axiale appliquée en compression simple.
𝜆 2
(𝜆 ) : Terme de correction qui tient compte de l’élancement 𝜆.
c
7.5.5 Critères de résistance

Pour les poutres en acier (selon l'élancement), des corrections sont à apporter :
Poutres courtes Poutres moyennes Poutres élancées

𝜆 ≤ 20 20 ≤ 𝜆 ≤ 100 𝜆 ≥ 100
Compression simple Rankine Euler
7.13
𝜎𝑝𝑐 . 𝑆 𝜎𝑝𝑐 . 𝑆
𝑃= 𝑃=
𝑃 = 𝜎𝑝𝑐 . 𝑆 𝜆 2 𝜆 2
1+( ) 2( )
𝜆𝑐 𝜆𝑐

73
Chapitre 8
Vibration de systèmes continus.
Introduction
En 1850 à Angers (France) : Un pont suspendu s'est effondré par le seul passage au pas cadencé
de soldats d'un régiment peu nombreux.
En 1940 au Tacoma État de Washington (USA) : Un pont suspendu s'est effondré par le seul
effet de rafales de vent régulières sans être particulièrement violentes (60 km/h).
Question légitime : Comment la simple marche d’un régiment peu nombreux ou des rafales de
vent peu violentes peuvent-elles détruire un pont ?
Réponses
 Il y a résonance lorsque la fréquence imposée par un agent extérieur (rafales de vent ou le
pas cadencé) devient égale à la fréquence propre du système mécanique (oscillation libre).
 Ce phénomène se manifeste par des oscillations très fortes pouvant entraîner la destruction
du système.
 Les oscillations induites sont bien plus fortes que les oscillations du phénomène qui impose
sa fréquence.
Depuis lors :
 Les codes militaires (monde entier) interdisent aux troupes de marcher au pas sur un pont,
 Dans la construction d’un pont (gratte-ciel aussi), les constructeurs tiennent compte de la
fréquence propre de l’ouvrage en lui donnant une valeur différente de celles de rafales de
vent.
Ce chapitre a pour objet l’étude des :

 Vibrations libres de poutres afin de trouver les Fréquences propres et les Modes propres.
 Vibrations forcées de poutres pour éviter le phénomène de résonnance.
Equations aux dérivées partielles : Rappels mathématiques

Equations aux dérivées partielles Solution générale
𝑑2 𝑇
+ 𝜔2 𝑇 = 0 𝑇(𝑡) = 𝐴 sin 𝜔𝑡 + 𝐵 cos 𝜔𝑡
2
𝑑𝑡
74
𝑑2 𝑇
− 𝜔2 𝑇 = 0 𝑇(𝑡) = 𝐴𝑒 𝜔𝑡 + 𝐵𝑒 −𝜔𝑡
𝑑𝑡2
𝑑4 𝑇
4
− 𝛽4 𝑇 = 0 𝑇(𝑡) = 𝐴 𝑐ℎ𝛽𝑡 + 𝐵 𝑠ℎ𝛽𝑡 + 𝐶 cos 𝛽𝑡 + 𝐷 sin 𝛽𝑡
𝑑𝑡
Avec :
A, B, C et D sont des constantes calculées à partir des conditions initiales et des conditions aux
limites.
𝑒 𝑡 + 𝑒 −𝑡 𝑒 𝑡 − 𝑒 −𝑡
𝑐ℎ𝛽𝑡 = et 𝑠ℎ𝛽𝑡 =
2 2
8.1 Vibrations longitudinales

Considérons une barre d’axe Ox en vibrations longitudinales dont la section S peut être fonction
de l’abscisse x. p=p (x, t) est une densité de force suivant l’axe Ox appliquée à la barre.
Figure 73 –Barre en vibrations longitudinales.
Nous proposons d’établir l’équation aux dérivées partielles du déplacement longitudinal u (x, t)
de la section d’abscisse x pendant la vibration longitudinale de la poutre.
Pour le tronçon élémentaire de longueur dx, les variables considérées sont :
𝜕𝑢 𝜕𝑢
𝜀= 𝑁 = −𝐸𝑆𝜀 = −𝐸𝑆 8.1
𝜕𝑥 𝜕𝑥
L'équilibre du tronçon élémentaire de longueur dx s'écrit :
𝑑𝑁 𝜕²𝑢
𝑝− = 𝜌𝑆
𝑑𝑥 𝜕𝑡²
Soit en introduisant l'équation (8.1) :
𝑝 𝜕 𝜕𝑢 𝑆 𝜕²𝑢
+ (𝑆 ) = 𝜌 8.2
𝐸 𝜕𝑥 𝜕𝑥 𝐸 𝜕𝑡²
8.1.1 Vibrations libres

Les vibrations libres sont les mouvements en l’absence de forces extérieures, donc p = 0. La
fonction u (x, t) est donc une intégrale de l’équation aux dérivées partielles :

75
𝜕 𝜕𝑢 𝜌𝑆 𝜕²𝑢
(𝑆 ) = 8.3
𝜕𝑥 𝜕𝑥 𝐸 𝜕𝑡²
Si la barre possède une section constante : S=cte, alors l'équation précédente s'écrira :
𝜕²𝑢 1 𝜕²𝑢 𝐸
= 𝑐² = 8.4
𝜕𝑥² 𝑐² 𝜕𝑡² 𝜌
Conditions aux limites possibles

 déplacement imposé nul aux extrémités :
𝑢(0, 𝑡) = 0 et / ou 𝑢(𝐿, 𝑡) = 0
 – effort imposé nul aux extrémités :
𝜕𝑢 𝜕𝑢
(0, 𝑡) = 0 et / ou (𝐿, 𝑡) = 0
𝜕𝑥 𝜕𝑥
Conditions initiales possibles
 déplacement imposé quel que soit x :
𝑁. 𝑥
𝑢(𝑥, 0) = 𝜑(𝑥) =
𝐸. 𝑆
 vitesse imposée nulle quel que soit x :
∂u
[ ] =0
∂t t=0
8.1.2 Fréquences propres et modes propres

On effectue une séparation des variables, le déplacement u (x, t) s'écrit :
𝑢(𝑥, 𝑡) = 𝑈(𝑥). 𝑇(𝑡) 8.5
L’équation d’équilibre devient :
1 𝑑²𝑈 1 1 𝑑²𝑇
= = 𝑐𝑡𝑒 8.6
𝑈 𝑑𝑥² 𝑐² 𝑇 𝑑𝑡²
On pose :
1 𝑑²𝑈 1 1 𝑑²𝑇
= 𝑓1 et / ou = 𝑓2
𝑈 𝑑𝑥² 𝑐² 𝑇 𝑑𝑡²
𝑓1 et 𝑓2 sont deux fonctions à variables indépendantes.
La constante cte est choisie négative pour assurer la stabilité de la solution en temps :
𝑑2 𝑈 𝜔 2
+ ( ) 𝑈=0
1 𝑑²𝑈 1 𝑑²𝑇 𝜔² 𝑑𝑥 2 𝑐
= =− Ou encore : 8.7
𝑈 𝑑𝑥² 𝑐² 𝑑𝑡² 𝑐² 𝑑²𝑇
{ 𝑑𝑡² + 𝜔²𝑇 = 0
Ces équations donnent les solutions suivantes :

76
𝜔 𝜔
𝑈(𝑥) = 𝐴 𝑠𝑖𝑛 𝑥 + 𝐵 𝑐𝑜𝑠 𝑥
{ 𝑐 𝑐 8.8
𝑇(𝑡) = 𝐶 𝑠𝑖𝑛 𝜔𝑡 + 𝐷 𝑐𝑜𝑠 𝜔𝑡
 Les constantes A, B, C et D sont calculées à partir des conditions initiales et des conditions
aux limites.
 Les i solutions des équations (8.8) sont les pulsations propres
 Les Ui sont les modes propres associés
 L'équation (8.3) est l’équation caractéristique des modes propres.
8.1.3 Exemples
Pour chaque cas de barre, appliquer les conditions aux limites en x=0 et x=L puis:
 Retrouver les pulsations propres et les modes propres puis dessiner leur allure,
 Retrouver la solution générale des déplacements.
Barre libre aux

deux extrémités
𝜋𝑐 𝜋 𝐸
Pulsations propres 𝜔𝑖 = 𝑖 =𝑖 √
𝐿 𝐿 𝜌
𝑖𝜋
Modes propres 𝑈𝑖 (𝑥) = 𝑐𝑜𝑠 𝑥
𝐿
∞
𝑖𝜋𝑥
Solution générale 𝑢(𝑥, 𝑡) = ∑ 𝑐𝑜𝑠 (𝐶𝑖 𝑠𝑖𝑛 𝜔𝑖 𝑡 + 𝐷𝑖 𝑐𝑜𝑠 𝜔𝑖 𝑡)
𝐿
𝑖=1
Barre encastrée-
libre
𝜋𝑐 𝜋 𝐸
Pulsations propres 𝜔𝑖 = (2𝑖 − 1) = (2𝑖 − 1) √
2𝐿 2𝐿 𝜌
𝜋
Modes propres 𝑈𝑖 (𝑥) = 𝑠𝑖𝑛(2𝑖 − 1) 𝑥
2𝐿
∞
𝜋𝑥
Solution générale 𝑢(𝑥, 𝑡) = ∑ 𝑠𝑖𝑛(2𝑖 − 1) (𝐶𝑖 𝑠𝑖𝑛 𝜔𝑖 𝑡 + 𝐷𝑖 𝑐𝑜𝑠 𝜔𝑖 𝑡)
2𝐿
𝑖=1
Barre encastrée-
encastrée

77
𝜋𝑐 𝜋 𝐸
Pulsations propres 𝜔𝑖 = 𝑖 =𝑖 √
𝐿 𝐿 𝜌
𝑖𝜋
Modes propres 𝑈𝑖 (𝑥) = 𝑠𝑖𝑛 𝑥
𝐿
∞
𝑖𝜋𝑥
Solution générale 𝑢(𝑥, 𝑡) = ∑ 𝑠𝑖𝑛 (𝐶𝑖 𝑠𝑖𝑛 𝜔𝑖 𝑡 + 𝐷𝑖 𝑐𝑜𝑠 𝜔𝑖 𝑡)
𝐿
𝑖=1
8.1.4 Orthogonalité des modes

L’équation caractéristique (8.3) pour les modes propres 𝑈𝑖 peut s'écrire en considérant les
solutions des exemples précédents :
(𝐸𝑆𝑈𝑖 ′)′ = −𝜔𝑖 ²𝜌𝑆𝑈𝑖
En multipliant chaque membre par le mode Uj et en intégrant sur la barre, on obtient :
𝐿 𝐿
∫ (𝐸𝑆𝑈𝑖 ′)′𝑈𝑗 𝑑𝑥 = −𝜔𝑖 ² ∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥 8.9
0 0
On procède par une intégrale par parties du premier terme de la dernière équation, on obtient :
𝐿 𝐿
𝐿
[𝐸𝑆𝑈𝑖 ′𝑈𝑗 ]0 − ∫ 𝐸𝑆𝑈𝑖 ′𝑈𝑗 ′𝑑𝑥 = −𝜔𝑖 ² ∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥
0 0
Le terme entre crochets est nul par ce qu'aux extrémités (x=0 et x=L), la barre est soit encastrée
(𝑈𝑗 = 0) ou elle est libre (𝑈𝑖′ = 0). Il reste :
𝐿 𝐿
∫ 𝐸𝑆𝑈𝑖 ′𝑈𝑗 ′𝑑𝑥 = 𝜔𝑖 ² ∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥 8.10
0 0
On prend cette fois, l’équation caractéristique (8.3) pour les modes propres 𝑈𝑗 :
(𝐸𝑆𝑈𝑗 ′)′ = −𝜔𝑗 ²𝜌𝑆𝑈𝑗
De la même manière que pour les modes Ui en multipliant chaque membre par le mode Ui et en
intégrant sur la barre, on obtient :
𝐿 𝐿
∫ 𝐸𝑆𝑈𝑖 ′𝑈𝑗 ′𝑑𝑥 = 𝜔𝑗 ² ∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥 8.11
0 0
Les deux équations (8.10) et (8.11) permettent d'écrire :
𝐿
(𝜔𝑖2 − 2
𝜔𝑗 ) ∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥 =0 8.12
0
Pour 𝑖 ≠ 𝑗, 𝜔𝑖2 ≠ 𝜔𝑗2 , alors :

𝐿
∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥 = 0 ∀𝑖 ≠ 𝑗 8.13
0
Deux modes différents sont orthogonaux par rapport à l’opérateur de masse 𝜌𝑆.

78
L'équation (8.13) associée aux équations (8.10) ou (8.11) et (8.9) permettent d'écrire :
𝐿 𝐿
∫ 𝐸𝑆𝑈𝑖 ′𝑈𝑗 ′𝑑𝑥 = − ∫ (𝐸𝑆𝑈𝑖 ′)′𝑈𝑗 𝑑𝑥 = 0 𝑖≠𝑗 8.14
0 0
Deux modes différents sont orthogonaux par rapport à l’opérateur de raideur ES.
Pour 𝜌𝑆 =cte et 𝐸𝑆 =cte et pour deux modes différents on a :
𝐿 𝐿 𝐿
∫ 𝑈𝑖 𝑈𝑗 𝑑𝑥 = ∫ 𝑈𝑖 ′𝑈𝑗 ′𝑑𝑥 = ∫ 𝑈𝑖 ′′𝑈𝑗 𝑑𝑥 = 0 𝑃𝑜𝑢𝑟 𝜌𝑆 = 𝑐𝑡𝑒, 𝐸𝑆 = 𝑐𝑡𝑒 𝑒𝑡 𝑖 ≠ 𝑗 8.15
0 0 0
8.1.5 Normalisation
Lorsque i =j, on normalise le mode Ui de manière à ce que :
𝐿
∫ 𝜌𝑆𝑈𝑖2 𝑑𝑥 = 1 8.16
0
La condition d’orthonormalité des modes peut donc s’écrire :

𝐿
∫ 𝜌𝑆𝑈𝑖 𝑈𝑗 𝑑𝑥 = 𝛿𝑖𝑗 8.17
0
Avec cette normalisation équation (8.17), les équations (8.9) et (8.10) s'écrivent pour i=j :
𝐿 𝐿
∫ 𝐸𝑆𝑈𝑖′2 𝑑𝑥 = − ∫ (𝐸𝑆𝑈𝑖 ′)′𝑈𝑖 𝑑𝑥 = 𝜔𝑖2 8.18
0 0
8.1.6 Exemple
Reprenons l'exemple de la Barre encastrée-libre :
On applique les conditions initiales puis en utilisant les propriétés de l'orthogonalité des modes,
on retrouve la solution générale des déplacements :
∝
𝟖𝑭𝑳 (−𝟏)𝒊 𝝅𝒙
𝒖(𝒙, 𝒕) = − 𝟐 ∑ 𝟐
𝒔𝒊𝒏(𝟐𝒊 − 𝟏) 𝒄𝒐𝒔(𝝎𝒊 𝒕)
𝝅 𝑬𝑺 (𝟐𝒊 − 𝟏) 𝟐𝑳
𝒊=𝟏
Vérification :
En statique, le déplacement de la barre 𝑢(𝑥 = 𝐿, 𝑡 = 0) est :
𝐹𝐿
𝑢(𝑥 = 𝐿, 0) =
𝐸𝑆
En dynamique, calculer le déplacement de la barre 𝑢(𝑥 = 𝐿, 0) pour les 4 premiers modes :
Le premier mode 8𝐹𝐿 𝐹𝐿

𝑢(𝐿, 𝑡) = − 2
= 0.81
𝜋 𝐸𝑆 𝐸𝑆
Les deux premiers modes 𝐹𝐿 1 𝐹𝐿
𝑢(𝐿, 𝑡) = 0.810 (1 + 2 ) = 0.900
𝐸𝑆 3 𝐸𝑆
79
Les trois premiers modes 𝐹𝐿 1 1 𝐹𝐿
𝑢(𝐿, 𝑡) = 0.81 (1 + 2 + 2 ) = 0.933
𝐸𝑆 3 5 𝐸𝑆
Les quatre premiers modes 𝐹𝐿 1 1 1 𝐹𝐿
𝑢(𝐿, 𝑡) = 0.81 (1 + 2 + 2 + 2 ) = 0.95
𝐸𝑆 3 5 7 𝐸𝑆
Les n premiers modes 𝐹𝐿 1 1 1 𝐹𝐿
𝑢(𝐿, 𝑡) = 0.81 (1 + 2 + 2 + 2 + ⋯ ) =
𝐸𝑆 3 5 7 𝐸𝑆
8.1.7 Vibrations forcées

Lorsque qu’on force la vibration par un effort f(x, t), l’équation d’équilibre (8.3) devient :
𝜕 𝜕𝑢 𝜕 2𝑢
(𝐸𝑆 ) − 𝜌𝑆 2 + 𝑓(𝑥, 𝑡) = 0 8.19
𝜕𝑥 𝜕𝑥 𝜕𝑡
On cherche une solution décomposée dans la base modale :
𝒖(𝒙, 𝒕) = ∑ 𝒒𝒋 (𝒕)𝑼𝒋 (𝒙) 8.20

𝒋=𝟎

′
∞ ∞
𝜌𝑆 (∑ 𝑞̈ 𝑗 (𝑡)𝑈𝑗 (𝑥)) − (∑ 𝑞𝑗 (𝑡)𝐸𝑆𝑈𝑗′ (𝑥)) = 𝑓(𝑥, 𝑡) 8.21

𝑗=0 𝑗=0
En multipliant chaque membre de l'équation (8.21) par Ui et en intégrant le long de la barre on

obtient :
′
𝐿 ∞ 𝐿 ∞ 𝐿
∫ 𝜌𝑆𝑈𝑖 (∑ 𝑞̈ 𝑗 (𝑡)𝑈𝑗 (𝑥)) 𝑑𝑥 − ∫ 𝑈𝑖 (∑ 𝑞𝑗 (𝑡)𝐸𝑆𝑈𝑗′ (𝑥)) 𝑑𝑥 = ∫ 𝑈𝑖 𝑓(𝑥, 𝑡)𝑑𝑥 8.22
0 𝑗=0 0 𝑗=0 0
On pose :
𝐿
𝑄𝑖 (𝑡) = ∫ 𝑈𝑖 𝑓(𝑥, 𝑡)𝑑𝑥
0
𝑄𝑖 (𝑡) est la projection de la force imposée sur le mode i.

En utilisant les propriétés d’orthonormalité des modes, il reste :
𝑞̈ 𝑖 (𝑡) + 𝜔𝑖2 𝒒𝒊 (𝒕) = 𝑄𝑖 (𝑡) 𝑖 = 0… ∞ 8.23
La résolution du problème de vibrations forcées se ramène à la résolution d’un ensemble de
systèmes scalaires à un degré de libertés indépendant.

80
8.2 Vibrations de torsion
Considérons un cylindre d’axe Ox en vibrations de torsion dont la section S peut être fonction de
l’abscisse x. M=M (x, t) est une densité de couple suivant l’axe Ox appliquée au cylindre.
Figure 74 –Cylindre en vibrations de torsion.
Nous proposons d’établir l’équation aux dérivées partielles de l'angle de rotation  (x, t) autour
de l'axe Ox de la section d’abscisse x pendant la vibration de torsion du cylindre.
𝜕𝛼
𝜃= 𝐶 = 𝑀𝑡 = 𝐺𝐼0 𝜃 8.24
𝜕𝑥
L'équilibre du tronçon élémentaire de longueur dx s'écrit :

𝑑𝑀𝑡
= 𝜌𝐼0 𝛼̈ 8.25
𝑑𝑥
Soit :
𝜕 𝜕𝛼 𝜌𝐼0 𝜕²𝛼
(𝐼0 ) = 8.26
𝜕𝑥 𝜕𝑥 𝐺 𝜕𝑡²
Cette forme est identique à celle obtenue pour un problème de vibration longitudinale de barre.
Les solutions sont du même type.
8.3 Vibrations transversales de flexion

Considérons une poutre d’axe Ox en vibrations de flexion dont la section S peut être fonction de
l’abscisse x.
Figure 75 –Poutre en vibrations de flexion.

81
Nous proposons d’établir l’équation aux dérivées partielles du déplacement transversal v (x, t) de
la section d’abscisse x pendant la vibration de flexion de la poutre.
Le déplacement transversal : 𝑣(𝑥, 𝑡) 8.27
𝜕𝑣
La rotation de la section (hypothèse de Bernoulli) : 𝜃(𝑥, 𝑡) = 8.28
𝜕𝑥
𝜕𝜃
La courbure : 𝑣̈ (𝑥, 𝑡) = 8.29
𝜕𝑥
Le moment de flexion : 𝑀𝑓 = 𝐸𝐼𝑣̈ 8.30
𝜕𝑀𝑓
Effort tranchant : 𝑇=− 8.31
𝜕𝑥
Théorème de la résultante dynamique appliqué sur le tronçon élémentaire de longueur dx :
𝜕𝑇 𝜕²𝑣
− 𝜌𝑆 8.32
𝜕𝑥 𝜕𝑡²
L'équation (8.31) dans l'équation (8.32) donne :
𝜕²𝑀𝑓 𝜕²𝑣
+ 𝜌𝑆 =0 8.33
𝜕𝑥² 𝜕𝑡²
En introduisant l'équation (8.30), on obtient finalement :
𝜕 2𝑣 𝜕 2𝑣 𝜕²𝑣
2
(𝐸𝐼 2
) + 𝜌𝑆 =0 8.34
𝜕𝑥 𝜕𝑥 𝜕𝑡²
Conditions aux limites possibles

 déplacement imposé nul aux extrémités :
𝑣(0, 𝑡) = 0 et / ou 𝑣(𝐿, 𝑡) = 0
 – Rotation imposée nulle aux extrémités :
𝜕𝑣 𝜕𝑣
(0, 𝑡) = 0 et / ou (𝐿, 𝑡) = 0
𝜕𝑥 𝜕𝑥
 – Moment imposé nul aux extrémités :
𝜕²𝑣 𝜕²𝑣
(0, 𝑡) = 0 et / ou (𝐿, 𝑡) = 0
𝜕𝑥² 𝜕𝑥²
 – effort imposé nul aux extrémités :
𝜕 3𝑣 𝜕 3𝑣
(0, 𝑡) = 0 et / ou (𝐿, 𝑡) = 0
𝜕𝑥 3 𝜕𝑥 3
8.3.1 Fréquences et modes propres

On effectue une séparation des variables :
𝑣(𝑥, 𝑡) = 𝑉(𝑥). 𝑇(𝑡) 8.35
L’équation d'équilibre (8.34) devient :
82
𝐸𝐼 1 𝑑4 𝑉 1 𝑑²𝑇
4
=− = 𝑐𝑡𝑒 8.36
𝜌𝑆 𝑉 𝑑𝑥 𝑇 𝑑𝑡²
On a égalité de deux fonctions de variables indépendantes. Les deux fonctions sont donc égales à
une constante choisie positive pour assurer la stabilité de la solution en temps :
𝐸𝐼 1 𝑑4 𝑉 1 𝑑2 𝑇
4
=− 2
= +𝜔2 8.37
𝜌𝑆 𝑉 𝑑𝑥 𝑇 𝑑𝑡
Ou encore :
𝑑4𝑉
− 𝛽4𝑉 = 0
𝑑𝑥 4 4
𝜔2 𝜌𝑆 2
Avec 𝛽 = 2= 𝜔 8.38
𝑑²𝑇 𝛼 𝐸𝐼
{ 𝑑𝑡² + 𝜔²𝑇 = 0
Ces équations donnet les solutions suivantes :
𝑉(𝑥) = 𝐴 𝑐ℎ𝛽𝑥 + 𝐵 𝑠ℎ𝐵𝑥 + 𝐶 𝑐𝑜𝑠𝛽𝑥 + 𝐷 𝑠𝑖𝑛𝛽𝑥 𝜔2 𝜌𝑆 2
{ Avec 𝛽4 = = 𝜔 8.39
𝑇(𝑡) = 𝐸 𝑐𝑜𝑠 𝜔𝑡 + 𝐹 𝑠𝑖𝑛 𝜔𝑡 𝛼 2 𝐸𝐼
 Les constantes A, B, C et D sont calculées à partir des conditions initiales et des conditions
aux limites.
 Les i solutions des équations (8.39) sont les pulsations propres
 Les Vi sont les modes propres associés
 L'équation (8.34) est l’équation caractéristique des modes propres.
8.3.2 Exemples
Pour chaque cas de barre, appliquer les conditions aux limites en x=0 et x=L puis:
 Retrouver les pulsations propres et les modes propres puis dessiner leur allure,
 Retrouver la solution générale des déplacements.
Poutre en appuis
simples
𝑖𝜋 𝐸𝐼
Pulsations propres 𝛽𝑖 = 𝜔𝑖 = 𝑖 2 𝜋 2 √
𝐿 𝜌𝑆𝐿4
𝑖𝜋
Modes propres 𝑉𝑖 (𝑥) = 𝑠𝑖𝑛 𝑥
𝐿
∞
𝑖𝜋𝑥
Solution générale 𝑣(𝑥, 𝑡) = ∑ 𝑠𝑖𝑛 (𝐸𝑖 𝑐𝑜𝑠 𝜔𝑖 𝑡 + 𝐹𝑖 𝑠𝑖𝑛 𝜔𝑖 𝑡)
𝐿
𝑖=1

83
2
1.5 m5 m4
0.5
Barre encastrée- 0
libre -0.5
-1
-1.5 m1
m2
m3
-2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
Pulsations propres
𝑉𝑖 (𝑥) = 𝐴𝑖 ch𝛽𝑖 x+𝐵𝑖 sh𝛽𝑖 x-𝐴𝑖 cos𝛽𝑖 x-𝐵𝑖 .sin𝛽𝑖 x

Modes propres = −𝐴𝑖 (cos𝛽𝑖 x+𝑎𝑖 sin𝛽𝑖 x-ch𝛽𝑖 x-𝑎𝑖 sh𝛽𝑖 x)
𝛼𝑖 -0.7340 -1.0184 -0.9992 -1.0000 -0.9999
∞
𝑢(𝑥, 𝑡) = ∑ = −𝐴𝑖 (𝑐𝑜𝑠𝛽𝑖 𝑥 + 𝑎𝑖 𝑠𝑖𝑛𝛽𝑖 𝑥 − 𝑐ℎ𝛽𝑖 𝑥

Solution générale
𝑖=1
− 𝑎𝑖 𝑠ℎ𝛽𝑖 𝑥)(𝐸𝑖 𝑐𝑜𝑠 𝜔𝑖 𝑡 + 𝐹𝑖 𝑠𝑖𝑛 𝜔𝑖 𝑡)
8.3.3 Orthogonalité des modes

L’équation caractéristique (8.34) pour les modes propres 𝑉𝑖 peut s'écrire :
(𝐸𝐼𝑉𝑖 ′′ )′′ = 𝜔𝑖 ²𝜌𝑆𝑉𝑖
En multipliant chaque membre par le mode 𝑉𝑗 et en intégrant sur la barre, on obtient :
𝐿 𝐿
∫ (𝐸𝐼𝑉𝑖 ′′ )′′𝑉𝑗 𝑑𝑥 = 𝜔𝑖 ² ∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥 8.40
0 0
On procède par une double intégrale par parties du premier terme de la dernière équation, on
obtient :
𝐿 𝐿
𝐿 𝐿
[(𝐸𝐼𝑉𝑖 ′′ )′𝑉𝑗 ]0 − [𝐸𝐼𝑉𝑖 ′′ 𝑉𝑗 ′]0 + ∫ 𝐸𝐼𝑉𝑖 ′′ 𝑉𝑗 ′′𝑑𝑥 = 𝜔𝑖 ² ∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥
0 0
Les termes dans le premier crochet sont nuls par ce qu'aux extrémités (x=0 et x=L), la barre est
soit bloquée (𝑉𝑗 = 0) soit libre d’effort (𝑉′′′𝑖 = 0).
Les termes dans le deuxième crochet sont aussi nuls par ce qu'aux extrémités (x=0 et x=L), la
barre est soit à rotation bloquée (𝑉′𝑗 = 0) soit libre de moment (𝑉′′𝑖 = 0).
Il reste :
𝐿 𝐿
′′
∫ 𝐸𝐼𝑉𝑖 𝑉𝑗 ′′𝑑𝑥 = 𝜔𝑖 ² ∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥 8.41
0 0
On prend cette fois, l’équation caractéristique (8.34) pour les modes propres 𝑉𝑗 :

84
(𝐸𝐼𝑉𝑗 ′′ )′′ = 𝜔𝑖 ²𝜌𝑆𝑉𝑗
De la même manière que pour les modes Vi en multipliant chaque membre par le mode Vi et en
intégrant sur la barre, on obtient :
𝐿 𝐿
′′
∫ 𝐸𝐼𝑉𝑖 𝑉𝑗 ′′𝑑𝑥 = 𝜔𝑗 ² ∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥 8.42
0 0
Les deux équations (8.41) et (8.42) permettent d'écrire :
𝐿
(𝜔𝑖2 − 𝜔𝑗2 ) ∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥 = 0 8.43
0
Pour 𝑖 ≠ 𝑗, 𝜔𝑖2 ≠ 𝜔𝑗2 , alors :

𝐿
∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥 = 0 ∀𝑖 ≠ 𝑗 8.44
0
Deux modes différents sont orthogonaux par rapport à l’opérateur de masse 𝜌𝑆.
L'équation (8.44) associée aux équations (8.40) et (16.42) ou (8.42) permettent d'écrire :
𝐿 𝐿
∫ (𝐸𝐼𝑉𝑖 ′′ )′′𝑉𝑗 𝑑𝑥 = ∫ 𝐸𝐼𝑉𝑖 ′′ 𝑉𝑗 ′′𝑑𝑥 = 0 𝑖≠𝑗 8.45
0 0
Deux modes différents sont orthogonaux par rapport à l’opérateur de raideur EI.
Pour 𝜌𝑆 =cte et 𝐸𝐼 =cte et pour deux modes différents on a :
𝐿 𝐿 𝐿
′′′
∫ 𝑉𝑖 𝑉𝑗 𝑑𝑥 = ∫ 𝑉𝑖 ′ 𝑉𝑗 𝑑𝑥 = ∫ 𝑉𝑖 ′′ 𝑉𝑗 ′′𝑑𝑥 = 0 𝑃𝑜𝑢𝑟 𝜌𝑆 = 𝑐𝑡𝑒, 𝐸𝐼 = 𝑐𝑡𝑒 𝑒𝑡 𝑖 ≠ 𝑗 8.46
0 0 0
8.3.4 Normalisation
Lorsque i =j, on normalise le mode Ui de manière à ce que :
𝐿
∫ 𝜌𝑆𝑉𝑖2 𝑑𝑥 = 1 8.47
0
La condition d’orthonormalité des modes peut donc s’écrire :
𝐿
∫ 𝜌𝑆𝑉𝑖 𝑉𝑗 𝑑𝑥 = 𝛿𝑖𝑗 8.48
0
Avec cette normalisation équation (8.48) les équations (8.40) et (8.41) s'écrivent pour i=j :
𝐿 𝐿
′′2
∫ 𝐸𝐼𝑉𝑖 𝑑𝑥 = ∫ (𝐸𝐼𝑉𝑖 ′′ )′′𝑉𝑖 𝑑𝑥 = 𝜔𝑖2 8.49
0 0
8.3.5 Exemples
Reprenons l'exemple de la Poutre en appuis simples :

85
On applique les conditions initiales puis en utilisant les propriétés de l'orthogonalité des modes,
on retrouve la solution générale des déplacements :
∝
2𝐹𝐿3 1 𝑖𝜋 𝑖𝜋𝑥
𝑣(𝑥, 𝑡) = 4 ∑ 4 sin . sin cos 𝜔𝑖 𝑡
𝜋 𝐸𝐼 𝑖 2 𝐿
𝑖=1
Vérification :
En statique, la déformée en milieu de la poutre est :
𝐿 𝐹𝐿3
𝑣 (𝑥 = , 0) =
2 48𝐸𝐼
En dynamique, calculer le déplacement de la barre 𝑢(𝑥 = 𝐿, 0) pour les 3 premiers modes :
Le premier mode 𝐿 𝐹𝐿3

𝑢 ( , 𝑡) =
2 48.7𝐸𝐼
Les deux premiers modes 𝐿 𝐹𝐿3
𝑢 ( , 𝑡) =
2 48.11𝐸𝐼
Les trois premiers modes 𝐿 𝐹𝐿3
𝑢 ( , 𝑡) =
2 48.03𝐸𝐼
Les n premiers modes 𝐿 2𝐹𝐿3 1 1
𝑣 (𝑥 = , 𝑡) = 4
(1 + 4 + 4 + … )
2 𝐸𝐼𝜋 3 5
8.3.6 Vibrations forcées

Lorsque qu’on force la vibration par un effort f(x, t), l’équation d’équilibre (8.34) devient :
𝜕 2𝑣 𝜕 2𝑣 𝜕²𝑣
2
(𝐸𝐼 2
) + 𝜌𝑆 = 𝑓(𝑥, 𝑡) 8.50
𝜕𝑥 𝜕𝑥 𝜕𝑡²
On cherche une solution décomposée dans la base modale :
𝒗(𝒙, 𝒕) = ∑ 𝒒𝒋 (𝒕)𝑽𝒋 (𝒙) 8.51

𝒋=𝟎

′′
∞ ∞
𝜌𝑆 (∑ 𝑞̈ 𝑗 (𝑡)𝑉𝑗 (𝑥)) + (∑ 𝑞𝑗 (𝑡)𝐸𝐼𝑉𝑗′′ (𝑥)) = 𝑓(𝑥, 𝑡) 8.52

𝑗=0 𝑗=0
En multipliant chaque membre de l'équation (8.21) par Vi et en intégrant le long de la barre on

obtient :

86
′′
𝐿 ∞ 𝐿 ∞ 𝐿
∫ 𝜌𝑆𝑉𝑖 (∑ 𝑞̈ 𝑗 (𝑡)𝑉𝑗 (𝑥)) 𝑑𝑥 + ∫ 𝑉𝑖 (∑ 𝑞𝑗 (𝑡)𝐸𝐼𝑉𝑗′′ (𝑥)) 𝑑𝑥 = ∫ 𝑉𝑖 𝑓(𝑥, 𝑡)𝑑𝑥 8.53
0 𝑗=0 0 𝑗=0 0
On pose :
𝐿
𝑄𝑖 (𝑡) = ∫ 𝑉𝑖 𝑓(𝑥, 𝑡)𝑑𝑥
0
𝑄𝑖 (𝑡) est la projection de la force imposée sur le mode i.

En utilisant les propriétés d’orthonormalité des modes, il reste :
𝑞̈ 𝑖 (𝑡) + 𝜔𝑖2 𝒒𝒊 (𝒕) = 𝑄𝑖 (𝑡) 𝑖 = 0… ∞ 8.54
La résolution du problème de vibrations forcées se ramène à la résolution d’un ensemble de
systèmes scalaires à un degré de libertés indépendant.
8.4 Méthodes approchées

Souvent un problème de vibrations est difficile à résoudre. Dans de tels cas, on fait recours aux
méthodes approchées telles que la méthode de Rayleigh, la méthode de Rayleigh- Ritz ou la
méthode des éléments finis.
8.4.1 Energie cinétique et énergie potentielle

Énergie cinétique :
1 𝐿 𝜕𝑢 2 𝐿
1
𝐸c = ∫ 𝜌𝑆 ( ) 𝑑𝑥 = q̇ 2 (t) ∫ 𝜌𝑆 𝑈² 𝑑𝑥 8.55
2 0 𝜕𝑡 0 2
Énergie potentielle ou énergie de déformation :

L'énergie potentielle ou l'énergie de déformation s'exprime par Equation (6.14) :
1 𝐿 𝑁2
𝐸p = ∫ . 𝑑𝑥 Cas de vibrations longitudinales
2 0 𝐸𝑆
8.56
1 𝐿 𝑀𝑓𝑧 2
𝐸p = ∫ . 𝑑𝑥 Cas de vibrations longitudinales de flexion
2 0 𝐸𝐼𝐺𝑧
8.4.2 Méthode de Rayleigh : Approximation de la fréquence fondamentale

Un mode de vibration s’exprime sous la forme :
u(x; t) =q(t).U (x) Avec q(t)=A sint + B cost 8.57
 est la pulsation propre et U(x) est le mode propre.
Dans ce cas, l'énergie cinétique s'exprimera par :
1 𝐿 𝜕𝑢 2 𝐿
1
𝐸c = ∫ 𝜌𝑆 ( ) 𝑑𝑥 = q̇ 2 (t) ∫ 𝜌𝑆 𝑈² 𝑑𝑥 8.58
2 0 𝜕𝑡 0 2
Et l'énergie potentielle ou l'énergie de déformation s'exprimera par :
87
1 𝐿 𝑁2 2 (t)
𝐿
1 𝜕𝑈 2
𝐸p = ∫ . 𝑑𝑥 = q ∫ 𝐸𝑆 ( ) 𝑑𝑥 Cas de vibrations longitudinales
2 0 𝐸𝑆 0 2 𝜕𝑥
𝐿𝑀 2 𝐿 2 8.59
1 𝑓𝑧 2 (t)
1 𝜕²𝑉 Cas de vibrations longitudinales de
𝐸p = ∫ . 𝑑𝑥 = q ∫ 𝐸𝐼 ( ) 𝑑𝑥
2 0 𝐸𝐼𝐺𝑧 0 2 𝜕𝑥² flexion
Dans un système conservatif, l’énergie mécanique est constante :
𝜕
𝐸M = 𝐸p + 𝐸c = 𝑐𝑡𝑒 Ou encore (𝐸 + 𝐸c ) = 0 8.60
𝜕𝑥 p
Ce qui permet d’écrire :
𝐿1
∫0 2 𝐸𝑆(𝑈 ′ )2 𝑑𝑥
𝜔2 = = 𝑅(𝑈) Cas de vibrations longitudinales
𝐿1
∫0 2 𝜌𝑆 𝑈² 𝑑𝑥
𝐿1
8.61
∫0 2 𝐸𝐼(𝑉 ′′ )2 𝑑𝑥
𝜔2 = = 𝑅(𝑈) Cas de vibrations longitudinales de flexion
𝐿1
∫0 2 𝜌𝑆 𝑉² 𝑑𝑥
R(U) est le quotient de Rayleigh de U. Il ne dépend que de la forme de U(x)
Approximation de Rayleigh
On considère Ua une approximation du mode exacte U(x).
L’approximation Ua de U(x) doit être cinématiquement admissible : satisfait aux conditions aux
limites géométriques du problème.
Dans ce cas la pulsation approchée 𝜔0′ est :
𝜔0′2 = 𝑅(𝑈𝑎 ) 8.62
Remarques
 On montre que la pulsation fondamentale (premier mode) du système 𝜔0 est inférieure ou
égale à la pulsation approchée 𝜔0′ :
𝜔02 = 𝑅(𝑈0 ) ≤ 𝜔0′2 = 𝑅(𝑈𝑎 ) 8.63
U0 est le mode propre fondamentale du système (premier mode).
 La fréquence fondamentale du système f0 est inférieure ou égale à la fréquence approchée 𝑓0′
√𝑅(𝑈𝑎 )
𝑓0 ≤ 𝑓0′ = 8.64
2𝜋
 Le choix de la forme de 𝑈𝑎 joue un rôle important dans l’approximation.
 𝑈𝑎 doit avoir un sens physique
 Des choix possibles de forme approchée Ua sont :
o Solution exacte d’un problème voisin,
o Déformée statique de la structure sous son poids propre

88
Exemples
Vibrations longitudinales Barre + ressort
Figure 76 –Vibrations longitudinales de Barre + ressort
𝐸 𝐸
Pulsations exactes (numériques) 𝜔0 = 2.029√ 𝜔1 = 4.913√
𝜌𝐿2 𝜌𝐿2
Rayleigh : premier mode de la même
poutre libre à droite (∆𝑟 = 4%) 𝐸
𝜋𝑥 𝜔0 ′ = 2.11√
𝑈𝑎 (𝑥) = 𝑠𝑖𝑛 𝜌𝐿2
2𝐿
Vibrations transversale étude d’un pont
Figure 77 –Vibrations transversale d’un pont

Un pont est composé d’une simple dalle de béton de largeur b = 2m, longueur L = 20m épaisseur
h = 50cm. Elle est encastrée à gauche dans une zone rocheuse et appuyée à droite sur un sol mou.
Ce sol est considéré comme élastique de raideur k = 106N/m, E = 35GPa; = 2500kg/m3;
Modèle1 : poutre encastré- libre
𝐸𝐼
𝜔0 ² = 1.8754
Solution exacte : 𝜌𝑆𝐿4
𝑓0 = 0.755 𝐻𝑧
Approximation très grossière (∆𝑟 = 27%) 𝜔0′ = 6.04 𝑟𝑑/𝑠
𝑣(𝑥) = 𝑥 2 𝑓0 ′ = 0.961 𝐻𝑧
Approximation RDM (∆𝑟 = 1.45%)
𝐹 𝐿𝑥 2 𝑥 3 𝜔0′′ = 4.816 𝑟𝑑/𝑠
𝑣(𝑥) = ( − ) 𝑓0 ′′ = 0.766 𝐻𝑧
⏟ 2
𝐸𝐼 6
=1
Modèle2 : poutre encastrée-libre avec poids uniformément réparti
Approximation RDM poids réparti (∆r = 0.4%)

𝜌𝑆 𝐿²𝑥 2 𝑥 4 𝐿𝑥 3 𝜔0′′′ = 4.766 𝑟𝑑/𝑠
𝑣(𝑥) = ( + − ) 𝑓0 ′′′ = 0.758 𝐻𝑧
⏟
2𝐸𝐼 2 12 3
=1

89
Modèle3 : poutre encastré- libre + appui élastique à droite
Approximation (RDM) (∆𝑟 = 2.70%)

𝜌𝑆 𝐿²𝑥 2 𝑥 4 𝐿𝑥 3 𝑓0 (𝐸𝐹) = 1.554 𝐻𝑧
𝑣(𝑥) = ( + − ) 𝑓0 = 1.596 𝐻𝑧
⏟
2𝐸𝐼 2 12 3
=1
Approximation (RDM). Appui trop raide k =
107N/m (∆𝑟 = 65%)
𝑓0 (𝐸𝐹) = 2.890 𝐻𝑧
𝜌𝑆 𝐿²𝑥 2 𝑥 4 𝐿𝑥 3
𝑣(𝑥) = ( + − ) 𝑓0 = 4.506 𝐻𝑧 ≫ 𝑓0 (𝐸𝐹)
⏟
2𝐸𝐼 2 12 3
=1
8.4.3 Méthode de Rayleigh- Ritz : Approximation des premières fréquences

La méthode de Rayleigh-Ritz consiste à chercher une approximation des modes de vibration
dans un espace de dimension N engendré par N fonctions i choisies.
On chercher des solutions de la forme :
𝒒𝟏 𝝓𝟏
𝐍
. .
𝑾(𝒙, 𝒚, 𝒛) = ∑ 𝒒𝒊 𝝓𝒊 (𝒙, 𝒚, 𝒛) = {𝒒}𝑻 {𝝓} avec {𝑞}={ . } et {𝝓}={ } 8.65
.
𝒊=𝟏 𝒒𝑵 𝝓𝑵
Le quotient de Rayleigh dépend uniquement des qi et peut être mis sous la forme d’un rapport de
deux formes quadratiques :
{𝑞}𝑇 𝐾{𝑞}
𝑅(𝑊) = = 𝜔2 8.66
{𝑞}𝑇 𝑀{𝑞}
Où K est appelée matrice de raideur et M matrice de masse.
On minimise le quotient de Rayleigh par rapport aux paramètres qi :
𝜕𝑅 𝐾{𝑞} − 𝜔2 𝑀{𝑞}
=0 Ce qui donne =0 8.67
𝜕{𝑞} {𝑞}𝑇 𝑀{𝑞}
Les quantités i qui assurent le minimum du quotient de Rayleigh sont donc les termes qui
annulent le déterminant :
[𝐾 − 𝜔2 𝑀] = 0 8.68
On peut montrer que i est une approximation par l’excès de la ième fréquence propre du
système.
Les vecteurs {qi}, associés aux i, qui sont solutions du problème aux valeurs propres se
déduisent par l'équation :
[𝐾 − 𝜔2 𝑀]{𝑞} = 0 8.69
Les vecteurs {qi} trouvés donnent les approximations des N premiers modes propres du système
(équation 8.65).

90
Approximation de Rayleigh- Ritz
 La qualité de l’approximation des modes décroit avec i,
 La qualité de l’approximation des pulsations est meilleure que celle de l’approximation des
modes,
 La qualité de l’approximation augmente avec la dimension de la base N.
Exemples
Vibrations longitudinales Barre + ressort (Figure 76)

𝐸 𝐸
Pulsations exactes (numériques) 𝜔0 = 2.029√ 𝜔1 = 4.913√
𝜌𝐿2 𝜌𝐿2
Rayleigh : premier mode de la même 𝐸
poutre libre à droite (∆𝑟 = 4%) 𝜔0 ′ = 2.11√
𝜋𝑥 𝜌𝐿2
𝑈𝑎 (𝑥) = 𝑠𝑖𝑛
2𝐿 (∆𝑟 = 4%)
Rayleigh-Ritz 2D par les 2 premiers 𝐸 𝐸
modes de la poutre libre à droite : 𝜔1 " = 2.066√ 2 𝜔2 " = 4.940√
𝜋𝑥 3𝜋𝑥 𝜌𝐿 𝜌𝐿2
𝝓𝟏 (𝑥) = 𝑠𝑖𝑛 𝝓 (𝑥) = 𝑠𝑖𝑛
2𝐿 𝟐
2𝐿 (∆𝑟 = 1.82%) (∆𝑟 = 0.55%)
Remarque
La méthode du quotient de Rayleigh en 2D permet une meilleur approximation par rapport à la
méthode en 1D

91
Chapitre 9
Calcul de structures par la méthode
des éléments finis
Ce cours utilise des calculs éléments finis sous le logiciel RDM6 Gratuit à télécharger depuis :
http://iut.univ-lemans.fr/ydlogi/index.html
9.1 Les principes

La Méthode des Eléments Finis consiste à :
 découper la structure en des subdivisions de formes simples
 et à choisir une approximation du déplacement sur chaque subdivision.
C’est une méthode de Ritz ’par morceaux’ qui s’adapte aux géométries les plus complexes.
Figure 78 –Découpage de la structure en des subdivisions de formes simples
9.1.1 Définitions
Figure 79 – Eléments finis et nœuds d'une structure.
 Les subdivisions sont appelées : Eléments finis,

92
 et les connexions entre éléments sont appelées : Nœuds.
9.1.2 Famille d’éléments

Il existe 3 grandes familles d’éléments. Les éléments uni-dimensionnels, bi-dimensionnels, tri-
dimensionnels.
 Les éléments uni-dimensionnels En 1D,
Le modèle uni-dimensionnel s’appuie sur les hypothèses restrictives de la RDM.
(Voir chapitre 2)
 Les éléments bi-dimensionnels En 2D et tri-dimensionnels En 3D,
Le modèle bi-dimensionnels ou tri-dimensionnels utilise, selon les cas, les théories des
membranes, des plaques et des coques ou exploite les relations de l’élasticité.
(Voir chapitre 3)
9.1.3 Exemples
Figure 80 – Elément triangles à 3 nœuds
 (ui ,vi ,wi) sont les composantes des déplacements pour le nœud i.
 Si l’élément introduit 3 nœuds (i, j, k) alors le nombre d’inconnus en 2D est 6. On range ces
inconnus dans une colonne appelée vecteur des déplacements nodaux de l’élément e.
𝑢𝑖
𝑣𝑖
𝑢𝑗
⟨𝑈𝑒 ⟩ =
𝑣𝑗
𝑢𝑘
[𝑣𝑘 ]
 Si le maillage introduit n nœuds alors le nombre d’inconnus est 2n (en 2D) ou 3n(en 3D). On
range ces inconnus dans une colonne appelée vecteur des déplacements globaux de la
structure.
𝑢1
𝑣1
𝑢2
⟨𝑈⟩ = 𝑣2
..
𝑢𝑛
[𝑣𝑛 ]

93
9.1.4 Approximation des déplacements
Le déplacement d'un point M autre que les nœuds :
𝑢𝑖
𝑣𝑖
𝑢 𝑁11 𝑁12 𝑁13 … 𝑁1 2𝑛 𝑢𝑗
⃗ ∗ (𝑀) = [ ] = [
𝑢 ] = [𝑁𝑒 (M)]⟨𝑈𝑒 ⟩ 9.1
𝑣 𝑁21 𝑁22 𝑁13 … 𝑁2 2𝑛 𝑣𝑗
𝑢𝑘
[𝑣𝑘 ]
⃗ ∗ (𝑀) est l’approximation des déplacements réels 𝑢
𝑢 ⃗ (𝑀)
[𝑁𝑒 (𝑀)] : est la matrice d’interpolation élémentaire.
Nij (M) : sont les fonctions d’interpolation élémentaires.
⟨𝑈𝑒 ⟩ : regroupe les composantes des déplacements des nœuds de l’élément (e).
9.1.5 Déformations
𝑢𝑖
𝜀11 𝜕𝑢/𝜕𝑥 𝑣𝑖
𝐵11 𝐵12 𝐵13 … 𝐵1 2𝑛
𝑢𝑗
𝜀 (𝑀) = [ 𝜀22 ] = [ 𝜕𝑣/𝜕𝑦 ] = [𝐵21 𝐵22 𝐵13 … 𝐵2 2𝑛 ] = [𝐵𝑒 (M)]⟨𝑈𝑒 ⟩ 9.2
2𝜀12 𝑣𝑗
𝜕𝑢/𝜕𝑦 + 𝜕𝑣/𝜕𝑥 𝐵31 𝐵32 𝐵33 … 𝐵3 2𝑛
𝑢𝑘
[𝑣𝑘 ]
9.1.6 Contraintes
Lois de comportement :
𝜎11 𝑐11 𝑐12 𝑐13 𝜀11
𝜀
𝜎(𝑀) = [𝜎22 ] = [𝑐21 𝑐22 𝑐23 ] [ 22 ] = [C][𝐵𝑒 (M)]⟨𝑈𝑒 ⟩ 9.3
𝜎12 𝑐31 𝑐32 𝑐33 2𝜀12
Pour trouver la solution recherchée en déplacements, efforts, contraintes et déformations, on
utilise l’énergie de déformation ou la formulation intégrale du principe des travaux virtuels. Ces
deux méthodes utilisent toutes les deux un calcul en fonction des déplacements nodaux.
Si le comportement est linéaire, on obtient un système linéaire du type (voir démonstration aux
paragraphes 9.2 et 9.3):
[𝑀]⟨𝑈̈⟩ + [𝐾]⟨𝑈⟩ = ⟨𝐹⟩ 9.4
Avec :
[𝑀] : Matrice de masse globale,
[𝐾] : Matrice de raideur globale,
⟨𝑈⟩ : Vecteur des déplacements nodaux.
⟨𝐹⟩ : Vecteur de forces nodales.
Ou encore si le problème est statique, le système s’écrit :

94
[𝐾]⟨𝑈⟩ = ⟨𝐹⟩ 9.5
La solution d’un problème de structure mécanique ⟨𝑈⟩ s’obtient-en résolvant le système linéaire.
9.2 Résolution d’un problème de statique 1D

9.2.1 Ecriture matricielle en utilisant l’énergie de déformation
Soit un corps élastique soumis à n forces Pi aux points Ai leur imposant des déplacements Ui
Figure 81 –Corps élastique soumis à n forces Pi aux points Ai
Le travail total fournit ou l’énergie élastique 1 𝑛 1

𝑊= ∑ 𝑃𝑖 𝑈𝑖 = ⟨𝑃⟩𝑇 ⟨𝑈⟩
correspondante (1F de Clapeyron) 2 𝑖=1 2
Le déplacement d’un point Ai s’écrit : 𝑛

𝑎𝑖𝑗 sont les coefficients d’influence ou de 𝑈𝑖 = ∑ 𝑎𝑖𝑗 𝑃𝑗 = [𝑎]⟨𝑃⟩
𝑗=1
souplesse, et [a] matrice de souplesse
1
2F de Clapeyron : W= ⟨𝑃⟩𝑇 [𝑎]⟨𝑃⟩
2
Introduisons la matrice de raideur [k]
⟨𝑃⟩ = [𝐾]⟨𝑈⟩
[K] est carrée et symétrique et : [𝐾] = [𝑎]−1
1 1
Enfin l’expression de l’énergie : 𝑊 = ⟨𝑈⟩𝑇 ⟨𝐾⟩𝑇 [𝑎][𝑘]⟨𝑈⟩ = ⟨𝑈⟩𝑇 [𝐾]⟨𝑈⟩
2 2
9.2.2 Analyse globale et assemblage

L'énergie de déformation de la structure est la somme des énergies de déformations élémentaires :
1 𝑇 [𝐾]⟨𝑈⟩
𝑛
𝑊= ⟨𝑈⟩ = ∑ 𝑊𝑒 9.6
2 𝑒=1
⟨U⟩ est le vecteur déplacement global (ensemble des déplacements nodaux) et [K] est la matrice
de rigidité globale. La construction de [K] à partir des matrices de rigidité élémentaires s'appelle
assemblage. Pour réaliser cette opération, on exprime les énergies élémentaires en fonction des
variables globales :

95
1 n 1 n 1
W= ⟨U⟩T [𝐾]⟨U⟩ = ∑ ⟨U e ⟩T [k e ]⟨U e ⟩ = ∑ ⟨U⟩T [K e ]⟨U⟩
2 e=1 2 e=1 2
1 T
n
= ⟨U⟩ (∑ [K e ]) ⟨U⟩
2 e=1
D’où l’expression de la matrice de rigidité globale :

n
[𝐾] = ∑ [K e ] 9.7
e=1
Dans la matrice [K e ] obtenue par expansion de [k e ], les seuls termes non nuls sont les
termes associés aux degrés de liberté de l'élément (e).
De même, le travail des forces extérieures s’écrit :
1 n
e
n 1 n 1
Wext = ⟨U⟩T ⟨F⟩ = ∑ Wext =∑ ⟨U e ⟩T ⟨𝑓 e ⟩ = ∑ ⟨U⟩T ⟨F e ⟩
2 e=1 e=1 2 e=1 2
1 n
= ⟨U⟩T (∑ ⟨F e ⟩)
2 e=1
D’où l’expression du vecteur force global :

n
⟨F⟩ = ∑ ⟨F e ⟩ 9.8
e=1
L’égalité W = Wext permet d’écrire enfin l’équation du système :

[𝑘]⟨𝑈⟩ = ⟨𝐹⟩ 9.9
C'est l'équation du système linéaire cherchée (équation 9.4 ou 9.5)
Résolution
Partition des degrés de liberté
Méthode directe ou résolution par blocs (ou de partition), nécessite de modifier (réarranger) les
termes des matrices (élimination de ligne - colonne). Elle est utilisée dans le cas où les ddl sont
limités.
[𝑘11 ][𝑘12 ] 𝑈𝑖 𝐹
[ ] { } = { 𝑐} 9.10
[𝑘21 ][𝑘22 ] 𝑈𝑐 𝐹𝑖
𝑈𝑖 : sont les déplacements libres (inconnus)
𝑈𝑐 : sont les déplacements imposés (connus)
𝐹𝑐 : Forces correspondant aux déplacements libres (connues)
𝐹𝑖 : Forces correspondantes aux déplacements imposés (inconnues) réactions d’appuis
Ce qui donne :
{𝑈𝑖 } = [𝑘11 ]−1 ({𝐹𝑐 } − [𝑘12 ]{𝑈𝑐 })
9.11
{𝐹𝑖 } = [𝑘21 ]{𝑈𝑖 } + [𝑘22 ]{𝑈𝑐 }
Cas particulier fréquent :
Si {𝑈𝑐 } = {0} :

96
{𝑈𝑖 } = [𝑘11 ]−1 {𝐹𝑐 }
9.12
{𝐹𝑖 } = [𝑘21 ]{𝑈𝑖 }
Seules les termes [𝑘11 ] et [𝑘21 ] de la matrices de raideur globale sont à calculer.
9.2.3 Exemple1 : Analyse statique d'un système constitué de ressorts
Figure 82 –Système constitué de ressorts
Bilan :
Nœuds Ddl Eléments
u1
3 3 : 1 ddl/nœud (u1,u2,u3) : 〈𝑈〉 = {u2} 5
u3
Système de 3 équations à 3 inconnues
Figure 83 –Système constitué de ressorts idéalisé
Equations d'équilibre pour chaque élément

Écriture des matrices élémentaires avec l'ensemble des ddl.
Pour chaque élément, on écrit l'équilibre sous la forme de l'équation (9.5) :
 Elément 1
97
matrice de rigidité de
Equation d'équilibre Sous forme matricielle
l'élément
(1)
𝑘1 0 0 u1 𝐹1 𝑘1 0 0
(1)
k1u1 =F1 (1)
[0 0 0] {u2} = { 0 } [𝑘](1) = [ 0 0 0]
0 0 0 u3 0 0 0 0
 De même, on obtient pour les autres éléments :
Elément1 Elément2 Elément3 Elément4 Elément5
𝑘1 0 0 u1 𝑘2 −𝑘2 0 u1 𝑘3 −𝑘3 0 u1 𝑘4 0 −𝑘4 u1 0 0 0 u1
[0 0 0] u2 [−𝑘2 𝑘2 0] u2 [−𝑘3 𝑘3 0] u2 [ 0 0 0 ] u2 [0 𝑘5 −𝑘5 ] u2
0 0 0 u3 0 0 0 u3 0 0 0 u3 −𝑘 4 0 𝑘4 u3 0 −𝑘5 𝑘5 u3
Tableau 11. Matrice de rigidité des cinq éléments du système.
Assemblage des caractéristiques élémentaires

La matrice de rigidité globale est obtenue à partir de l'équation (9.7) :
𝑘1 + 𝑘2 + 𝑘3 + 𝑘4 −𝑘2 − 𝑘3 −𝑘4
[ −𝑘2 − 𝑘3 𝑘2 + 𝑘3 + 𝑘5 −𝑘5 ]
−𝑘4 −𝑘5 𝑘4 + 𝑘5
Résolution du problème
Les rigidités et les forces externes étant connues, il suffit de résoudre le système linéaire obtenu.
9.2.4 Elément de barre en traction- compression
Figure 84 –Elément barre en traction- compression
E:module d'Young (N/m2) – S:section (m2) – L:longueur(m)
Bilan :
 2 nœuds
 2 ddl : 1 ddl par nœud (ui,uj)
 1 élément
𝑢𝑖
Les déplacements nodaux de l’élément barre : ⟨𝑈𝑒 ⟩ = ( )
𝑢𝑗
𝑎0
Soit l'approximation nodale : 𝑈 ∗ (𝑥) = 𝑎0 + 𝑎1 𝑥 = 〈1, 𝑥〉 (
)
𝑎1
Pour construire l'approximation nodale, nous devons identifier aux nœuds la valeur de
l'approximation et les déplacements nodaux :
98
𝑢𝑖
⟨𝑈𝑒 ⟩ = ( ) = [1 0] (𝑎0) = [𝑃𝑒 (M)⟨𝐴⟩
𝑢𝑗 1 𝐿 𝑎1
1 0 𝑎0
Avec : [𝑃𝑒 (M) = [ ] et ⟨𝐴⟩ = ( )
1 𝐿 𝑎1
Encore, on peut écrire :
⟨𝐴⟩ = [𝑃𝑒 (M)]−1 ⟨𝑈𝑒 ⟩
D’où l'expression du déplacement :
𝑈 ∗ (𝑥) = 〈1, 𝑥〉⟨𝐴⟩ = 〈1, 𝑥〉[𝑃𝑒 (M)]−1 ⟨𝑈𝑒 ⟩ = [𝑁𝑒 (M)]⟨𝑈𝑒 ⟩
Avec :
[𝑁𝑒 (M)] = 〈1, 𝑥〉[𝑃𝑒 (M)]−1
Enfin l'expression du déplacement :
𝑥 𝑥 𝑢𝑖
𝑈 ∗ (𝑥) = 〈1 − , 〉 ( ) 9.13
𝐿 𝐿 𝑢𝑗
Le calcul de l’énergie de déformation permet d’avoir :
𝐸𝑆 1 −1
𝐾= [ ] 9.14
𝐿 −1 1
Transformation des matrices
Lorsque l’orientation d'un élément est différente de celle définie dans le repère de référence.
Figure 85 –Elément barre d'orientation quelconque.
𝑢𝑖
𝑣𝑖 𝐸𝑆 𝐴 −𝐴 2
⟨𝑈⟩ = ( ) 𝐾𝑅2 = [ ] 𝐴 = [ 𝑐𝑜𝑠 𝜃 𝑐𝑜𝑠𝜃𝑠𝑖𝑛𝜃]
𝑢𝑗 𝐿 −𝐴 𝐴 𝑐𝑜𝑠𝜃𝑠𝑖𝑛𝜃 𝑠𝑖𝑛2 𝜃
𝑣𝑗
9.2.5 Elément de poutre en flexion

Matrice de rigidité élémentaire (poutre type Bernoulli : pas de cisaillement transverse)

99
Figure 86 –Elément de poutre en flexion.
E:module d'Young (N/m2) – I:inertie de flexion (m4) – L:longueur(m)

Bilan :
 2 nœuds
 4 ddl : 2 ddl/nœud (vi, i, vj, j)
 1 élément
𝑣𝑖
𝜃𝑖
Les déplacements nodaux de l’élément poutre : ⟨𝑈⟩ = (𝑣𝑗 )
𝜃𝑗
Soit l'approximation nodale :
𝑣∗ (𝑥) 𝑣∗ (𝑥) = a0 + a1 x + a2 x² + a3 𝑥3
( ) = ( )
𝜃∗ (𝑥) 𝜃∗ (𝑥) = a1 + 2a2 x + 3a3 x²
Ou encore sous forme matricielle :
𝑎0
𝑣∗ (𝑥) 1, 𝑥, 𝑥2 , 𝑥3 𝑎1
( ∗ )=( 2 ) (𝑎2)
𝜃 (𝑥) 0, 1, 2𝑥, 3𝑥
𝑎3
Les conditions aux limites permettent d’écrire :
𝑣𝑖 1 0 0 0 𝑎0
𝜃𝑖 0 1 0 0 𝑎1
⟨𝑈𝑒 ⟩ = (𝑣𝑗 ) = [ 2 3 ]( ) = [𝑃𝑒 (M)] ⟨𝐴⟩
1 𝐿 𝐿 𝐿 𝑎2
𝜃𝑗 0 1 2𝐿 3𝐿2 𝑎3
Avec :
⟨𝐴⟩ = [𝑃𝑒 (M)]−1 ⟨𝑈𝑒 ⟩
D’où
𝑣∗ (𝑥) 1, 𝑥, 𝑥2 , 𝑥3 −1
( ∗ )=( 2 ) [𝑃𝑒 (M)] ⟨𝑈𝑒 ⟩
𝜃 (𝑥) 0, 1, 2𝑥, 3𝑥
Et :

100
avec 𝑥2 𝑥3
𝑁1 = 1 − 3 2 + 2 3
𝐿 𝐿
𝑣𝑖 𝑥 𝑥2 𝑥3
𝑁2 = 𝐿 ( − 2 2 + 3 )
𝜃𝑖 𝐿 𝐿 𝐿
𝑣 ∗ (𝑥) = 〈𝑁1 𝑁2 𝑁3 𝑁4 〉 (𝑣𝑗 ) = 〈𝑁𝑒 〉〈𝑈𝑒 〉 9.15
𝑥2 𝑥3
𝜃𝑗 𝑁3 = 3 2 − 2 3
𝐿 𝐿
𝑥2 𝑥3
𝑁4 = 𝐿 (− 2 + 3 )
{ 𝐿 𝐿
Le calcul de l’énergie de déformation permet d’avoir :
12 6𝐿 −12 6𝐿
𝐸𝐼 6𝐿 4𝐿2 −6𝐿 2𝐿2
𝐾 = 3[ ] 9.16
𝐿 −12 −6𝐿 12 −6𝐿
6𝐿 2𝐿2 −6𝐿 4𝐿2
9.2.6 Exemples
Analyser les structures suivantes :
 Structure de 1 élément de type "poutre en flexion"
Calculer les déplacements et des réactions aux liaisons.

 Structure de 2 éléments de type "poutre en flexion"
Calculer les déplacements et des réactions aux liaisons.
9.3 Résolution d’un problème 2D

9.3.1 Ecriture matricielle par utilisation du principe des travaux virtuels PTV
Le principe des travaux virtuels (PTV) pour un problème de déformation de structures continue D
⃗ ∗ , sous forme de l''équation :
s’écrit quel que soit le champ de déplacement virtuel 𝑢
⃗ ∗ 𝑑𝑣 = ∫ ⃗⃗⃗
∫ 𝜌𝛾 𝑢 𝑓𝑣 𝑢⃗ ∗ 𝑑𝑣 + ∫ ⃗⃗𝑓𝑠 𝑢
⃗ ∗ 𝑑𝑠 − ∫ 𝜎̿ : 𝜀̿∗ 𝑑𝑣 9.17
𝐷 𝐷 𝜕𝐷 𝐷
Les déplacements réels et virtuels ont la même forme (équation 9.1) :
⃗ (𝑀) = [𝑁𝑒 (𝑀)]⟨𝑢𝑒 ⟩
𝑢 et ⃗ ∗ (𝑀) = [𝑁𝑒 (𝑀)]⟨𝑢∗ ⟩
𝑢

101
Le premier terme de l'équation du PTV :
∫ 𝜌⟨𝑢̈ 𝑒 ⟩𝑇 [𝑁𝑒 (M)]𝑇 [𝑁𝑒 (M)]⟨𝑢∗ ⟩𝑑𝑣 = ⟨𝑢̈ 𝑒 ⟩𝑇 (∫ 𝜌[𝑁𝑒 (M)]𝑇 [𝑁𝑒 (M)]𝑑𝑣) ⟨𝑢∗ ⟩
𝐷 𝐷
∗
⟩𝑇
= ⟨𝑢̈ 𝑒 [𝑀𝑒 ]⟨𝑢 ⟩
Le deuxième et troisième terme de l''équation du PTV :
∫ ⃗⃗⃗
𝑓𝑣 𝑢⃗ ∗ 𝑑𝑣 + ∫ ⃗⃗𝑓𝑠 𝑢
⃗ ∗ 𝑑𝑠 = (∫ [𝑓𝑣 ]𝑇 [𝑁𝑒 (M)] 𝑑𝑣 + ∫ [𝑓𝑠 ]𝑇 [𝑁𝑒 (M)]𝑑𝑠) ⟨𝑢∗ ⟩
𝐷 𝜕𝐷 𝐷 𝜕𝐷
= [𝐹𝑒 ]𝑇 ⟨𝑢∗ ⟩
Le quatrième terme de l''équation du PTV :
𝑇
∫ 𝜎̿ : 𝜀̿∗ 𝑑𝑣 = ∫ [[C][𝐵𝑒 (M)]⟨𝑈𝑒 ⟩ ] [𝐵𝑒 (M)]⟨𝑢∗ ⟩𝑑𝑣 = ⟨𝑈𝑒 ⟩𝑇 (∫ 𝐵𝑒 (M)𝑇 [C]𝑇 [𝐵𝑒 (M)]𝑑𝑣) ⟨𝑢∗ ⟩
𝐷 𝐷 𝐷
= ⟨𝑈𝑒 ⟩𝑇 [𝐾𝑒 ]⟨𝑢∗ ⟩

En remplaçant les quatre termes de (9.17), le principe des travaux virtuels s'écrira donc quel que
⃗ ∗:
soit le champ de déplacement virtuelle 𝑢
⟨𝑢̈ 𝑒 ⟩𝑇 [𝑀]⟨𝑢∗ ⟩ = [𝐹]𝑇 ⟨𝑢∗ ⟩ − ⟨𝑈𝑒 ⟩𝑇 [𝐾]⟨𝑢∗ ⟩
Ou encore :
⟨𝑢̈ 𝑒 ⟩𝑇 [𝑀𝑒 ] = [𝐹𝑒 ]𝑇 − ⟨𝑈𝑒 ⟩𝑇 [𝐾𝑒 ]
Par construction [𝑀𝑒 ] et [𝐾𝑒 ] sont symétriques. Finalement on peut écrire :
[𝑀𝑒 ]⟨𝑢̈ 𝑒 ⟩ + [𝐾𝑒 ]⟨𝑈𝑒 ⟩ = [𝐹𝑒 ] 9.18
C'est l'équation du système linéaire cherché pour l'élément e (équation 9.4 ou 9.5). Avec :
[𝑀𝑒 ] = ∫ 𝜌[𝑁𝑒 (M)]𝑇 [𝑁𝑒 (M)]𝑑𝑣 Matrice masse de l’élément e 9.19

𝐷
[𝐾𝑒 ] = ∫ [𝐵𝑒 (M)]𝑇 [C]𝑇 [𝐵𝑒 (M)]𝑑𝑣 Matrice raideur de l’élément e 9.20
𝐷
[𝐹𝑒 ]𝑇 = ∫ [𝑓𝑣 ]𝑇 [𝑁𝑒 (M)] 𝑑𝑣 + ∫ [𝑓𝑠 ]𝑇 [𝑁𝑒 (M)]𝑑𝑠 Vecteur force de l’élément e 9.21
𝐷 𝜕𝐷
Ce calcul permet après identification de trouver la matrice de masse [𝑀𝑒 ] et de raideur [𝐾𝑒 ] de
l’élément e.
Les matrices globales de masse [𝑀] et de raideur [𝐾] de la structure s’obtiennent par assemblage
des matrices de masse [𝑀𝑒 ] et de raideur [𝐾𝑒 ] de tous les éléments e.

102
9.3.2 Différents types d’éléments en 2D
Triangles
Figure 87 –Elément Triangle à 3 et 6 nœuds.
Triangles de degré 1 : triangles à 3 nœuds, fonctions linéaires

𝑢(𝑥, 𝑦) = 𝑎0 + 𝑎1𝑥 + 𝑎2𝑦
𝑢
⃗ (𝑀) = [ ]
𝑣(𝑥, 𝑦) = 𝑏0 + 𝑏1𝑥 + 𝑏2𝑦
Triangles de degré 2 : triangles à 6 nœuds, polynômes de degré 2
𝑢(𝑥, 𝑦) = 𝑎0 + 𝑎1𝑥 + 𝑎2𝑦 + 𝑎3𝑥𝑦 + 𝑎4𝑥² + 𝑎5𝑦²
𝑢
⃗ (𝑀) = [ ]
𝑣(𝑥, 𝑦) = 𝑏0 + 𝑏1𝑥 + 𝑏2𝑦 + 𝑏3𝑥𝑦 + 𝑏4𝑥² + 𝑏5𝑦²
Quadrilatères
Figure 88 –Elément Quadrilatère ou carré à 4, 8 et 9 nœuds.
Quadrilatères de degré 1 : carrés à quatre nœuds, fonctions linéaires)

𝑢(𝑥, 𝑦) = 𝑎0 + 𝑎1𝑥 + 𝑎2𝑦 + 𝑎3𝑥𝑦
𝑢
⃗ (𝑀) = [ ]
𝑣(𝑥, 𝑦) = 𝑏0 + 𝑏1𝑥 + 𝑏2𝑦 + 𝑏3𝑥𝑦
Quadrilatères de degré 2 : carrés à 8 ou 9 nœuds, polynômes de degré 2
𝑢(𝑥, 𝑦) = 𝑎0 + 𝑎1𝑥 + 𝑎2𝑦 + 𝑎3𝑥𝑦 + 𝑎4𝑥² + 𝑎5𝑦² + 𝑎6𝑥²𝑦 + 𝑎7𝑥𝑦² + 𝑎8𝑥²𝑦²
⃗ (𝑀) = [
𝑢 ]
𝑣(𝑥, 𝑦) = 𝑏0 + 𝑏1𝑥 + 𝑏2𝑦 + 𝑏3𝑥𝑦 + 𝑏4𝑥² + 𝑏5𝑦² + 𝑏6𝑥²𝑦 + 𝑏7𝑥𝑦² + 𝑏8𝑥²𝑦²
9.3.3 Résolution d’un problème plan

Par raison de simplification des calculs, on utilisera dans cette partie des éléments triangles à 3
nœuds.

103
Figure 89 –Utilisation des éléments Triangle à 3 nœuds.
Pour l'élément Triangle à 3 nœuds, le déplacement d'un point M s'écrit :

𝑢(𝑥, 𝑦) = 𝑎0 + 𝑎1𝑥 + 𝑎2𝑦
⃗ (𝑀) = [
𝑢 ]
𝑣(𝑥, 𝑦) = 𝑏0 + 𝑏1𝑥 + 𝑏2𝑦
Ou encore sous forme matricielle :
𝑎0
𝑎1
𝑢(𝑥, 𝑦) 1 𝑥 𝑦 0 0 0 𝑎2
⃗ (𝑀) = [
𝑢 ]=[ ]
𝑣(𝑥, 𝑦) 0 0 0 1 𝑥 𝑦 𝑏0
𝑏1
[𝑏2]
Déplacement et Matrice d’interpolation élémentaire

Les conditions aux limites permettent d’écrire aux points i, j et k :
𝑢𝑖 1 𝑥𝑖 𝑦𝑖 0 0 0 𝑎0
𝑣𝑖 0 0 0 1 𝑥𝑖 𝑦𝑖 𝑎1
𝑢𝑗 1 𝑥𝑗 𝑦𝑗 0 0 0 𝑎2
⟨𝑈𝑒 ⟩ = = = [𝑃𝑒 (M)]⟨𝐴⟩
𝑣𝑗 0 0 0 1 𝑥𝑗 𝑦𝑗 𝑏0
𝑢𝑘 1 𝑥𝑘 𝑦𝑘 0 0 0 𝑏1
[𝑣𝑘 ] [0 0 0 1 𝑥𝑘 𝑦𝑘] [𝑏2]
Avec :
⟨𝐴⟩ = [𝑃𝑒 (M)]−1 ⟨𝑈𝑒 ⟩
𝑢(𝑥, 𝑦) 1 𝑥 𝑦 0 0 0
⃗ (𝑀) = [
𝑢 ]=[ ] [𝑃 (𝑀)]−1 ⟨𝑈𝑒 ⟩ = [𝑁𝑒 (𝑀)]⟨𝑈𝑒 ⟩
𝑣(𝑥, 𝑦) 0 0 0 1 𝑥 𝑦 𝑒
Avec :
1 𝑥𝑖 𝑦𝑖 0 0 0 𝑎𝑖 0 𝑎𝑗 0 𝑎𝑘 0
0 0 0 1 𝑥𝑖 𝑦𝑖 𝑏𝑖 0 𝑏𝑗 0 𝑏𝑘 0
1 𝑥𝑗 𝑦𝑗 0 0 0 1 𝑐𝑖 0 𝑐𝑗 0 𝑐𝑘 0
[𝑃𝑒 (M) = et [𝑃𝑒 (M)]−1 =
0 0 0 1 𝑥𝑗 𝑦𝑗 ∆𝑒 0 𝑎𝑖 0 𝑎𝑗 0 𝑎𝑘
1 𝑥𝑘 𝑦𝑘 0 0 0 0 𝑏𝑖 0 𝑏𝑗 0 𝑏𝑘
[0 0 0 1 𝑥𝑘 𝑦𝑘] [ 0 𝑐𝑖 0 𝑐𝑗 0 𝑐𝑘 ]
Avec :

104
1 𝑥𝑖 𝑦𝑖
∆𝑒 = |1 𝑥𝑗 ⃗⃗ ∧ ⃗⃗⃗
𝑦𝑗 | = 2 ∗ 𝑆𝑖𝑗𝑘 = 2 ∗ 𝑎𝑖𝑟𝑒 𝑑𝑢 𝑡𝑟𝑖𝑎𝑔𝑙𝑒(𝑖, 𝑗, 𝑘) = ‖𝑖𝑗 𝑖𝑘 ‖
1 𝑥𝑘 𝑦𝑘
Et :
𝑎𝑖 = 𝑥𝑗𝑦𝑘 − 𝑥𝑘𝑦𝑗
{ 𝑏𝑖 = 𝑦𝑗 − 𝑦𝑘
𝑐𝑖 = 𝑥𝑘 − 𝑥𝑗
Les autres termes (indice j et k) sont obtenus par permutation circulaire des indices (i,j,k).
Et :
𝑎𝑖 0 𝑎𝑗 0 𝑎𝑘 0
𝑏𝑖 0 𝑏𝑗 0 𝑏𝑘 0
1 𝑥 𝑦 0 0 0 1 1 𝑥 𝑦 0 0 0 𝑐𝑖 0 𝑐𝑗 0 𝑐𝑘 0
[𝑁𝑒 (M)] = [ ] [𝑃 (M)]−1 = [ ]
0 0 0 1 𝑥 𝑦 𝑒 ∆𝑒 0 0 0 1 𝑥 𝑦 0 𝑎𝑖 0 𝑎𝑗 0 𝑎𝑘
0 𝑏𝑖 0 𝑏𝑗 0 𝑏𝑘
[ 0 𝑐𝑖 0 𝑐𝑗 0 𝑐𝑘 ]
1 𝑁1 0 𝑁2 0 𝑁3 0
= [ ]
𝑆𝑖𝑗𝑘 0 𝑁1 0 𝑁2 0 𝑁3
Avec :
𝑁1 = 𝑎𝑖 + 𝑏𝑖𝑥 + 𝑐𝑖𝑦
{ 𝑁2 = 𝑎𝑗 + 𝑏𝑗𝑥 + 𝑐𝑗𝑦
𝑁3 = 𝑎𝑘 + 𝑏𝑘𝑥 + 𝑐𝑘𝑦
Finalement :
𝑢(𝑥, 𝑦) 1 𝑁1 0 𝑁2 0 𝑁3 0
⃗ (𝑀) = [
𝑢 ]= [ ] ⟨𝑈𝑒 ⟩ 9.22
𝑣(𝑥, 𝑦) 𝑆𝑖𝑗𝑘 0 𝑁1 0 𝑁2 0 𝑁3
Déformations
𝑢𝑖
𝜀11 𝜕𝑢/𝜕𝑥 𝑏𝑖 0 𝑏𝑗 0 𝑏𝑘 0 𝑣𝑖
𝑢𝑗
𝜀 (𝑀) = [𝜀22 ] = [ 𝜕𝑣/𝜕𝑦 ] = [ 0 𝑐𝑖 0 𝑐𝑗 0 𝑐𝑘 ] = [𝐵𝑒 (M)]⟨𝑈𝑒 ⟩ 9.23
𝜀12 𝑣𝑗
𝜕𝑢/𝜕𝑦 + 𝜕𝑣/𝜕𝑥 𝑐𝑖 𝑏𝑖 𝑐𝑗 𝑏𝑗 𝑐𝑘 𝑏𝑘
𝑢𝑘
[𝑣𝑘 ]
Matrice raideur [Ke]

D'après l'équation (9.20) :
[𝐾𝑒 ] = ∫ [𝐵𝑒 (M)]𝑇 [C]𝑇 [𝐵𝑒 (M)]𝑑𝑣 = [𝐵𝑒 (M)]𝑇 [C]𝑇 [𝐵𝑒 (M)] ∫ 𝑑𝑣
𝐷 𝐷
𝑑𝑣 = 𝑒 𝑆𝑖𝑗𝑘 et e : épaisseur de la plaque

D’où :

105
[𝐾𝑒 ] = [𝐵𝑒 (M)]𝑇 [C]𝑇 [𝐵𝑒 (M)]. 𝑒. 𝑆𝑖𝑗𝑘
Ou encore :
𝑏𝑖 0 𝑐𝑖
0 𝑐𝑖 𝑏𝑖 𝑐 𝑐 𝑐 𝑏𝑖 0 𝑏𝑗 0 𝑏𝑘 0
11 11 11
𝑏𝑗 0 𝑐𝑗 𝑐
[𝐾𝑒 ] = [ 21 𝑐21 𝑐21 ] [ 0 𝑐𝑖 0 𝑐𝑗 0 𝑐𝑘 ] 𝑒 𝑆𝑖𝑗𝑘 9.24
0 𝑐𝑗 𝑏𝑗 𝑐 𝑐 𝑐
11 11 11 𝑐𝑖 𝑏𝑖 𝑐𝑗 𝑏𝑗 𝑐𝑘 𝑏𝑘
𝑏𝑘 0 𝑐𝑘
[0 𝑐𝑘 𝑏𝑘 ]
[𝐾𝑒 ] est une matrice (6,6)=(6,3)*(3,3)*(3,6)
9.3.4 Résolution
On utilise la méthode directe ou résolution par blocs (ou de partition) vue au chap2 dans le cas où
les ddl sont limités.
Au cas de plusieurs ddl, on utilise des méthodes numériques de résolution. Ces méthodes sont
vues en analyse Numérique. Pour les problèmes de statique on trouvera dans les codes EF deux
types de méthodes :
Les méthodes directes : Élimination de Gauss, décomposition LDU, ou Cholesky …
Les méthodes itératives de type Gauss-Seidel.
9.3.5 Exemple 2
Analyse d'une structure de type "poutre 2D en flexion"
F=10.000 N, Longueur L=1 m

E=210.000 MPa Hauteur h=0,2m
=0,3 Epaisseur e= 1m
Calculer les déplacements et des réactions aux liaisons en utilisant :

 Un maillage triangulaire de 3 nœuds,
 Un maillage rectangulaire de 8 nœuds,
Comparer puis conclure.

106
Annexe A
Jauges de déformations
 Les jauges de déformations sont des capteurs électriques collées sur des pièces élastiques
pour mesurer les déformations.
 La résistance électrique d'une jauge de déformation change légèrement lorsqu’une force
extérieure l'étire ou la comprime.
 Ce changement de résistance est directement proportionnel au déplacement
(déformation).
Les jauges de déformations standards

 Les jauges de déformations sont constituées par une grille conductrice en métal laminé très
fin en forme de zigzag.
 Leur épaisseur est seulement de quelques microns et elles sont solidaires d’un support isolant
assurant une stabilité mécanique et une isolation électrique.
 Les utilisateurs collent les jauges à la surface d’une pièce mécanique où ils veulent
déterminer les contraintes lorsque la pièce est sous charge.
 La jauge de déformation va s’allonger ou se compresser avec la surface de la pièce sur
laquelle elle est collée.
 Une valeur négative de la déformation correspond à une déformation en compression alors
qu’une valeur positive correspond à une déformation en traction.
Représentation réaliste d’une jauge Diagramme simplifié et agrandi d’une jauge.

107
Les rosettes de jauges de déformation
 Une rosette est un ensemble de deux jauges ou plus déjà fixées sur le même support.
 Utilisées quand vous devez mesurer les déformations dans deux ou plus directions
différentes
 L’utilisation de rosettes vous permet ainsi de gagner du temps et d’obtenir des résultats plus
justes et précis.
Jauges de déformation pour mesure de couple et cisaillement

Comme indiqué, les jauges de déformation fonctionnent mieux quand elles sont alignées avec la
déformation qu’elles mesurent.
Pour le montage en torsion :
 la contrainte et la déformation sont à 45° par rapport à l’axe de la barre (voir cours
d’élasticité).
 Des rosettes spéciales à deux jauges placées à 45° permettent de mesurer cette déformation.

108
Pont de Wheatstone : circuit de base
C’est un circuit le mieux adapté pour la mesure de petites variations électriques.
Il est constitué de quatre résistances identiques connectées entre elles selon une forme en losange.
Une tension Vi est appliquée en entrée entre deux connexions opposée. La tension Vo en sortie
est mesurée entre les deux autres connexions.
La tension Vo en sortie est exprimée en fonction de la tension Vi appliquée en entrée et les quatre
résistances selon la relation suivante :
𝑅3 𝑅4
𝑉𝑂 = 𝑉𝑖 ( − )
𝑅1 + 𝑅3 𝑅2 + 𝑅4
 Si les résistances R1, R2, R3 et R4 sont toutes égales, alors la tension en sortie Vo est nulle,
 Si une ou plusieurs des résistances est différente des autres, alors la tension Vo changera en
fonction de cette résistance (R1).
Influence Du nombre de jauges utilisé

La tension de sortie V0 varie en fonction de la variation des résistances :
Si une seule résistance varie Si deux résistances varient Si quatre résistances varient
V0 n’est pas linéaire, mais V0 est deux fois plus V0 est quatre fois plus
elle donnera de bons importante que dans le cas importante que dans le cas
résultats pour ∆R faibles de 1seule jauge, mais V0 est de 1seule jauge, et V0 est
toujours non linéaire linéaire

109
Facteur de jauge
∆𝑅/𝑅 ∆𝑅/𝑅
𝑘𝐴 = =
∆𝐿/𝐿 𝜀𝐴
𝑘𝐴 : est le facteur de jauge ou encore sensibilité du matériau, relie la déformation de la jauge 𝜀𝐴
au changement de sa résistance ∆𝑅/𝑅.
Équation du pont de jauges

La relation qui relie la déformation 𝜺 à la tension mesurée 𝑽𝟎 est de la forme :
𝑉0
𝜀=4
𝐾𝑉𝑖 𝑁
Avec :
𝜺 : la déformation
𝑽𝟎 : la tension OUT mesurée à travers le pont (V)
𝑽𝒊 : la tension IN imposé au pont (V)
K (GF) : le facteur de jauge
N : le nombre de bras actifs (nombre de jauges connectées)
La tension est alors multipliée par 106 pour donner un résultat en (micro déformation)
(déformation x 10-6)
Influence de la température
 Une jauge de déformation est une résistance. Elle s’échauffe légèrement quand vous la
raccorder dans un circuit électrique.
 Sa résistance change légèrement si sa température augmente. Ceci signifie que la mesure ne
sera pas correcte tant que sa température ne sera pas stabilisée.
 Il faut donc toujours laisser au minimum une minute aux jauges de l’appareil pour
s’échauffer avant de faire le zéro.

110

00 MDS 1920 PDF

Transféré par

Droits d'auteur :

Formats disponibles

00 MDS 1920 PDF

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

00 MDS 1920 PDF

Transféré par

Droits d'auteur :

Formats disponibles

Centre des Métiers d'Education et de Formation

Cycle Préparatoire d’Agrégation

Mécanique des structures

Année universitaire 2019 – 2020

Table des matières ................................................................................................ 1

M. El Allami Mécanique des structures

M. El Allami Mécanique des structures

M. El Allami Mécanique des structures

M. El Allami Mécanique des structures

Généralités et notions de base.

1.1 Buts de la mécanique des structures

Choix des matériaux

M. El Allami Mécanique des structures

Figure 1 –Différents comportement du matériau

1.2.2 La disposition de la matière

M. El Allami Mécanique des structures

1.2.3 Les actions mécaniques extérieures

1.2.4 Les déformations

M. El Allami Mécanique des structures

Figure 5 –Poutre après déformation

1.3 Principe fondamental de la statique

Système 3D (Oxyz) Système 2D (Oxy)

 proj. sur Ox des Fext=0  proj. sur Ox des Fext=0

 moments/ Ox des Fext=0

Tableau 1. Equations d’équilibre

M. El Allami Mécanique des structures

Tableau 2. Actions de liaison

1.4 Torseur des efforts de cohésion

Figure 6 –Poutre (E) coupée en deux portions (E1) et (E2).

(E1) étant en équilibre, on peut écrire : E  E1   E 2  E1  0

 E 2  E1  E  E1  E  E 2

1.4.2 Composantes du torseur de cohésion

Figure 7 – Résultante et moment du torseur de cohésion dans une section droite.

 N : effort normal dans la section S,

M. El Allami Mécanique des structures

1.5 Tenseur des contraintes

(a) (b) (c)

Figure 8 –coupure, facette et vecteur contrainte dans un solide.

M. El Allami Mécanique des structures

1.5.2 Contrainte normale/ contrainte tangentielle :

Figure 9 –Contrainte normale et contrainte tangentielle

1.5.3 Tenseur des contraintes

̿(M) sur les trois facettes du point M

M. El Allami Mécanique des structures

Contraintes en un point sur les facettes 𝐢, 𝐣, 𝐤

sur la facette 𝑖 : 𝑇 ⃗ = 𝜎̿(𝑀). 𝑖

sur la facette 𝑗 : 𝑇 ⃗ = 𝜎̿(𝑀). 𝑗

Soit en notation indicielle, le vecteur contrainte sur la facette 𝑖 :

Figure 11 –Contraintes sur une facette 𝑛⃗ quelconque.

M. El Allami Mécanique des structures

La contrainte tangentielle ⃗⃗⃗⃗

1.5.4 Torseur de cohésion : écriture locale

𝑑𝐹 (𝐸2 → 𝐸1) 𝑑𝐹 (𝐸2 → 𝐸1) ⃗ (𝑀, 𝑛⃗)𝑑𝑆

M. El Allami Mécanique des structures

Traction / compression des poutres.

Figure 12 – Poutre sollicitée à la traction

2.2 Essai d'traction

Figure 13 –Eprouvette d'essai de traction

M. El Allami Mécanique des structures