Université Nationale des Sciences,

Technologies, Ingénierie et Mathématiques

(UNSTIM)

Ecole Nationale Supérieure des travaux Publics

(ENSTP)

COURS D’ANALYSE
DYNAMIQUE DES STRUCTURES

Cycle d’Ingénieurs

Enseignant : O. Farouk DJIBRIL

Docteur-Ingénieur de Génie Civil
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Introduction
La dynamique est une partie de la mécanique des constructions qui traite
des méthodes de calcul des constructions soumises à des forces dynamiques
(variables). C’est-à-dire des forces dont l’intensité, la direction ou la position
varie avec le temps. L’analyse dynamique des structures est donc une des
parties les plus difficiles ou complexes de la mécanique de construction.
Sous l’effet des sollicitations dynamiques, l’importance des masses de
construction et celle de la force d’inertie n’est plus à démontrer. En effet, on
connaît des cas où une construction calculée avec un grand coefficient de
sécurité pour une charge statique s’est rompue sous l’effet de forces
périodiques relativement petites. Souvent, une construction rigide et très
résistante est inadéquate en présence de forces variables, alors qu’une
construction identique plus légère et à première vue moins résistante supporte
ces efforts facilement. C’est pourquoi, les questions d’oscillation et en général,
le comportement des systèmes élastiques sous l’action de charges variables
exigent du constructeur une attention particulière.

GENERALITES SUR LA DYNAMIQUE DES STRUCTURES

1. Charges ou forces dynamiques
1.1. Définitions
Les charges dynamiques sont les charges pour lesquelles la
vitesse de croissance est si rapide qu’on ne peut négliger les
forces d’inerties qui apparaissent pendant la déformation du
système. Les chocs sont aussi classés parmi les charges
dynamiques mêmes si quelques fois, on peut négliger pour eux les
forces d’inertie.
Les contraintes et les déformations dues aux charges dynamiques
sont appelées contraintes et déformations dynamiques.

2. Différentes sortes de forces dynamiques

2.1. Force immobile périodique

ENSTP 2021-2022 2
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Il s’agit de forces fixes évoluant selon une période bien déterminée. Elle
peut être continue ou discontinue. Une force périodique qui varie selon une loi
sinusoïdale est appelée harmonique ou vibratoire.

Exemple d’un groupe, d’un moteur.

3. Force impulsive (non périodique)

Qui surgit et s’éteint très vite. Il s’agit là des forces explosives (bombe,
grenade).

4. Force de battage (pieux, enclume, chute d’un corps, etc)

Celle-ci se produit en un lieu précis et est caractérisée par le changement
brusque de vitesse (force de frappe instantanée). De telle force peut être
périodique.

ENSTP 2021-2022 3
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

5. Charge roulante de valeur constante ou variable

Elle change de position en fonction du temps.
Exemple : train, automobile sur un pont, pont roulant, grue, etc.

6. Force sismique
Caractérisée par un grand désordre. Il est difficile de maîtriser l’évolution de
celle-ci.
Tremblement de terre – effets sismiques
On appelle effets sismiques, des sollicitations dues aux oscillations au
déplacement de la terre pendant le tremblement. Ces actions sont transmises
au bâtiment par le biais des fondations secouées. Il apparaît alors une
nouvelle force qu’il convient de prendre en compte dans le calcul de ces
bâtiments. Il s’agit là de la force d’inertie qui émane de la masse des
constructions.

Le vent
Le vent constitue une force dynamique. En effet, il est indispensable de
connaître la forme des oscillations horizontales libres des bâtiments à étages,
car celles-ci sont nécessaires pour le calcul d’un bâtiment soumis à l’action du
vent (coups du vent).

ENSTP 2021-2022 4
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Les caractéristiques dynamiques sont fonctions du schéma constructif du

bâtiment et du schéma théorique correspondant. L’inertie des masses d’étage
joue dans ce cadre un grand rôle.
Pour les bâtiments à étages, une vérification de l’accélération à laquelle se
produisent les oscillations du sommet du bâtiment sous l’effet des coups du
vent devient nécessaire.
Exemple : Pour assurer l’état normal des personnes se trouvant dans le
bâtiment, cette accélération ne doit pas dépasser 15 cm/s2.

1. 2. 3.
Formes d’oscillations

7. Types d’oscillations
Les oscillations sont les manifestations des charges dynamiques les plus
répandues. Elles sont classées en plusieurs types selon un certain nombre de
critères :
 Selon la loi de variation dans le temps, on a les types suivants
d’oscillations (voir fig 1.1)

Fig 1.1
ENSTP 2021-2022 5
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

- Les oscillations sinusoïdales ou périodiques simples (fig. 1,a) ;

- Les oscillations amorties (fig1.1,b) ;
- Les oscillations croissantes (fig.1.1,c) ;
- Les oscillations périodiques complexes fig(1.1,d) ;
- Les oscillations désordonnées (fig.1.1,e)
 Selon la liberté du mouvement, on a :
- Des oscillations libres ;
- Des oscillations forcées ou non libres.
 Selon la résistance (l’opposition au mouvement) du milieu, on a :
- Des oscillations amorties ;
- Des oscillations non amorties ;
 Selon le degré de liberté (DDL), on a :
- Des oscillations à un degré de liberté (1 DDL) ;
- Des oscillations à n degrés de liberté (nDDL) ;
- Des oscillations à une infinité de degré de liberté (∞DDL).
 Selon la direction de la charge dynamique, on a (voir fig.1.2) ;
- Des oscillations transversales (fig.1.2,a) ;
- Des oscillations longitudinales (fig.1.2,b).

Fig 1.2
Objectifs d’une étude dynamique
Les forces dynamiques entraînent l’oscillation des constructions sur lesquelles
elles agissent. Alors l’étude dynamique a pour but :
- de vérifier la résistance des constructions,
- de déterminer les déformations,
- de connaître la vitesse et l’accélération issues des sollicitations
dynamiques afin de les comparer aux valeurs admissibles.
Il existe deux méthodes fondamentales pour la résolution des problèmes de la
dynamique :
- Méthode statique basée sur l’équilibre du système. Celle-ci diffère des
équations de la statique de part un facteur complémentaire issu du principe
de Dalembert : la force d’inertie ;

ENSTP 2021-2022 6
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

d 2x
X  m 2  0
dt X  mx  0
d2y 
Y  m 2  0 ou
Y  my  0  l ' équation d ' équilibre
dt Z  mz  0 
d 2z
Z  m 2  0
dt
d 2x d2y d 2z
mX  m ; mY  m ; mZ  m
dt 2 dt 2 dt 2
Sont les forces d’inerties x, y, z déplacement du corps dans différentes
directions.
En résumé nous dirons que la méthode statique ou cinéto- statique consiste
à considérer le corps en équilibre sous l’action des charges dynamiques,
des forces d’inerties et des masses déposées sur elle.
- méthode énergétique est fondée sur la recherche des conditions d’équilibre
de l’énergie totale du système en tenant compte des forces d’inerties, soit
sur la recherche de l’énergie cinétique et élastique pendant le movement du
corps. En particulier, quand les forces de résistance au mouvement ne sont
pas tenues en compte, elle est fondée sur la loi de la conservation de
l’énergie.
C  U  const (1.2)
avec mv
2
v 2 ( x) (1.3)
C  i 1    m ( x) dx
2
   2 
masses ponctuelle s uniforméme nt réparties

1  M dx 2
N 2 dx T 2 dx 
2   EI  EA  GA 
U         (1.4)

Degré de liberté
On entend par degré de liberté, le nombre de coordonnées indépendantes
déterminant la position du système.
Le DDL d’un système élastique est constitué des DDL des différentes masses
qui sont disposées sur lui (voir fig.1.3)
On néglige le poids propre de la poutre.

ENSTP 2021-2022 7
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Fig 1.3.
a. Une masse rigide liée à un ressort a un degré de liberté puisque sa
position est définie par une seule coordonnée comptée à partir d’un
certain point (si la masse du ressort est négligeable).

b. Si on néglige la masse de l’arbre, on peut dire que le système a deux

degrés de liberté : il faut deux coordonnées angulaires déterminant la
rotation des disques rigides.

c. La position de la charge oscillante dans le plan de la figure est

déterminée par trois variables indépendantes :
- deux coordonnées du centre de gravité ;
- angle de rotation de la masse par rapport au centre de gravité.
Par conséquent, le degré de liberté est trois. Cependant, en négligeant le
moment d’inertie de la masse, on peut dire que le système a deux degrés de
liberté.
Trois degrés de liberté Quatre degrés de liberté

ENSTP 2021-2022 8
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

OSCILLATIONS DE SYSTEMES A UN
Chapitre II : DEGRE DE LIBERTE (étude de mouvement)

Tout système élastique bougé de son état d’équilibre par une sollicitation
extérieure effectue des mouvements oscillatoires. Il s’agira pour nous de faire
une étude dynamique de ces mouvements et de déboucher sur certaines
caractéristiques indispensables lors de l’analyse de ces systèmes.
Soit le système élastique suivant. Faisons le bilan de forces qui s’exercent sur
lui :

- Selon la seconde loi fondamentale (Newton) des équations de la

dynamique, la force Q s’exprime comme le produit de la masse m  par g
où
Q – force de pesanteur ; Q=mg
m – la masse kg.s 2 / cm en système technique 

g  981 cm / s 2 accélération de la pesanteur .

Sous l’action de Q, la console subit à elle seule un déplacement y ST par

rapport à la position d’équilibre. En déplaçant le système de sa position

d’équilibre statique, y ST tout en ramenant à Y  et en le relâchant, tout le

système subit des oscillations.

- La force de rétablissement de mouvement ou force réactive élastique du

système : R
La force R tend à ramener le système à la position d’équilibre statique. Il
s’agit donc d’une force négative. Son expression est la suivante :

ENSTP 2021-2022 9
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

R =-ry où
r est un coefficient de proportionnalité caractérisant la réaction de la poutre
au point fixé (en étude) pour un déplacement égal à l’unité.
1
y  r  1  r 

 - déplacement du point sous l’action d’une force unitaire.
Exemple du cas en examen.   B11 (méthode des forces) ;    1P

.Q 2 1 3
  ..   B11
2 3 EI 3 EI
1 3 EI
r  1  r   3
 

- Force de résistance : F
Elle est née de la résistance du milieu environnant c’est-à-dire l’eau, l’air…etc.
Elle dépend aussi du frottement interne des matériaux. Pour contourner les
difficultés mathématiques dues à la prise en compte de cette force, on suppose
qu’elle est proportionnelle à la vitesse des oscillations v  y ' , et qu’elle est
appliquée à la masse m son signe est négatif.
dy
F  k   ky ' où
dt
k le coefficient de proportionnalité au sens physique des oscillations.

- Force d’inertie : X
Selon le principe de Dalembert la force d’inertie X qui est égale au produit de
la masse par l’accélération des oscillations doit être prise en compte dans le
bilan. Cette force est dirigée dans le sens contraire de l’accélération du
mouvement donc est négatif.
d2y
X  m 2  my"
dt
L’équation de l’équilibre dynamique du système s’écrit.
Y  X  F  R  0

NB : La force Q ne figure pas dans cette équation comme elle est déjà prise en
compte dans l’équilibre statique.
ENSTP 2021-2022 10
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

En remplaçant les forces par leur valeur, on obtient une équation

différentielle linéaire et homogène.

k r k r
y" y ' y  0 ou ÿ  y  y  0 (1)
m m m m

I. Oscillations propres non amorties : F=0

r  1
L’équation (1) devient ÿ  y  0 r 
m  

1 1 1
ÿ . y  0 posons  2   
m m m
Cette équation différentielle admet pour solution de mouvement oscillatoire

y  A sin t  B cos t (2)

Alors la vitesse y  A cos t  B sin t

Utilisons les conditions à l’instant t  o à travers les paramètres initiaux pour
déterminer A et B.
y o  y (o )
to
y o  y (o)  v (o)

y o  y (o )  y (o )  B
y (o) v (o)
y o  y (o)  y (o)  A  A  
 
L’équation du mouvement devient alors

v (o )
y  y (0) cos t  sin t (4)


En posant A  C cos  et B  C sin  dans (2) il vient

y  C sin t    (5) où

C est l’amplitude des oscillations
 l’angle de phase
De même à travers les paramètres initiaux t  o

v 2 (0) B y (0)
C A2  B 2  y 2 (0)  ; tg  
 2
A v (o )

ENSTP 2021-2022 11
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

à travers (2) ; (4) et (5) il est évident que le mouvement de la masse est
harmonique. La période des oscillations T se détermine facilement de la
condition qu’en ajoutant à temps t le temps T, l’expression sous le signe sinus
augmente de 2  .
 (t  T )     t     T

T  2  T  2  période des oscillations



 c’est la fréquence circulaire c’est la pulsation, c’est –à-dire le nombre

d’oscillations en 2π seconde (donnée technique)

2
  2  1
  f  La fréquence
T T

On peut donc conclure que dans un système de degré de liberté 1 et non

amorti les oscillations sont toujours propres et ont une période et fréquence ne
dépendant pas des conditions initiales de mouvement. Cependant, l’amplitude
et la phase sont intimement liées aux conditions initiales.

 Q 
2 y ST m  
- La période T   2  m   2 en s ;  g  Q force de pesanteur
 g  Q  y 
 ST 

1  1 1 g 1
- La fréquence  f    en où s 1
T 2 2 m 2 y ST s

 La fréquence c’est le nombre d’oscillations en une seconde

ENSTP 2021-2022 12
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

 La fréquence circulaire, c’est le nombre d’oscillations en 2 

seconde
NB :
- Plus la rigidité d’un système est grande, plus  est petit ce qui augmente,
la fréquence des oscillations propres ;
- Le nombre d’oscillations en une minute devient
60 g
n 300
2 y ST  n y ST en cm
y ST
g  981  10 
2

Cette caractéristique est très utilisée en pratique.

NB : Dans le cas d’une masse oscillante, la position de la masse est

déterminée par un seul paramètre : l’angle de rotation ø.
Selon le principe de Dalembert le moment devient J m . ; J m étant le moment

d’inertie de la masse par rapport à l’axe de l’articulation. L’équation du

1
mouvement étant    2  0 avec    '.I m . et  2 
' I m
 : angle de torsion
' étant l’angle de rotation imprimé à la masse par l’application d’un moment
unitaire.

II. Oscillations propres amorties F≠0

L’équilibre dynamique est régi par l’équation différentielle suivante :
k r
ÿ y  y  0
m m
k 1 r
Posons  2 ;  2   ; l’équation différentielle transformée prend la
m m m

forme suivante :

ÿ  2  y   2 . y  0

Equation caractéristique r12  2  r1   2  0 ;

r1     2   2 ; r1'     2 ω 2

Etudions les différents cas de  2 2

ENSTP 2021-2022 13
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

1.     r1    i 1 ; r1'    i 1 avec 1   2   2 (fréquence

circulaire des oscillations amorties)
La solution complexe de l’équation différentielle dans ce cas serait :

y  e  t  A sin 1 t  B cos 1 t 

1 k k  2
où   ;   ; 1  2− 
2

m 2m 2
2 2
la période des oscillations T1  
1 2  2
2
si   0  T   T1 (diminution de la force de résistance)

Déterminons les coefficients A et B en se servant des conditions initiales au
temps t  o .
a.) y o  y (0)  y o  B

   
b.) y o  v (o)   e  t  A sin 1 t  B cos 1 t   e  t  1 A cos 1t  1 B sin 1t 
     
 B   1 A 

a. )  y o  B

t  o v (o )   B v (o )   y o
b. )   B  1 A  y o  v (o)  A  1

1


yo  B
v (o )   y o
A
1

L’équation prend alors la forme

 v (o )   y o 
y  e  t  . sin 1t  y o cos 1t 
 1 
  est le coefficient caractérisant la capacité d ' amortissement du système oscillatoire
Cette équation peut être améliorée en posant :
A  C cos  ; B  C sin 

y  e  t C cos  . sin 1t  C sin  . cos 1 t 

y  C e  t cos  . sin 1t  C sin  . cos 1t 
L’analyse de l’équation montre bien que l’amplitude des oscillations décroît
avec le temps : les oscillations s’amortissent.

ENSTP 2021-2022 14
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

y  C e  t sin 1t   

A 2  C 2 cos 2  sin  B
 C A 2  B 2 tg  
B  C . sin 
2 2 2
cos  A
2
 k 
2  vo  y o . 
 v (o )   y o   2m 
C y (0)  
2
  yo 
2

 1   1 
 
 
2
1  k  1
1      2 2
  avec r
m  2 m  

Il vient :
2
r  k 
1   
m  2 m 
B yo y o .1
tg   
A v (o )   y o k
v (o )  . yo
1 2m
y o .1
  arc tg
k
v (o )  . yo
2m

La dimension de l’amplitude suit une progression géométrique. En effet, pour

t  0 ; t  T ; t  2 T etc…. les amplitudes ont respectivement les valeurs

t  o  Ce
t  T  C.e  T
t  2 T  C.e  2  T

ENSTP 2021-2022 15
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Cherchons le logarithme néprien du rapport des amplitudes successives

(entre une période T)
kt
 k
a e 2m t n  1  t n 
  n n  n ktn 1  n.e 2 m
a n 1 
e 2m
kT
kT
 n e 2m
 T
2m

   T1

Cette valeur qui caractérise la vitesse d’amortissement du mouvement

oscillatoire s’appelle décrément logarithmique des oscillations, c’est la
caractéristique principale des oscillations.
 c’est le coefficient d’amortissement

 2  m.  2   2
k  2 m  2 m
T1 2

Cette formule montre bien que l’amortissement est bien lié à la fréquence des
oscillations.
De tout ce qui précède, on peut conclure que :
1. la force résistante fait diminuer la fréquence des oscillations propres ;
2. la fréquence des oscillations ne dépend pas de l’amplitude ;
3. les oscillations varient progressivement jusqu’à s’éteindre (s’amortir)

2ème cas   

L’équation caractéristique s’écrit :

r1     2  0 ; r1'     2  0
y  C e  2
 t

 D e   2  e C e  D e
 t t 2  t 2

ou
y  e  t  A sh 2 t  B ch  2 t 
Dans ce cas, le mouvement ne sera pas oscillatoire c’est-à-dire que lorsqu’on
bouge le système de sa position d’équilibre le système tendra à reprendre
progressivement sa position d’équilibre.
Au plus le système peut dépasser une fois sa position d’équilibre et y revenir.
Ce genre de mouvement est non périodique : aucune vibration n’a eu lieu
après la perturbation.

ENSTP 2021-2022 16
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

3ème cas :   w

r1   ; r1'  

y  e  t  A  Bt 
Ce genre de mouvement est aussi non périodique.

ENSTP 2021-2022 17
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

OSCILLATIONS FORCEES D’UN

Chapitre III :
SYSTEME A UN DEGRE DE LIBERTE
(AVEC AMORTISSEMENT) F ≠ 0

Soit le système suivant soumis à une force dynamique harmonique

P (t )  P sin  t

Le bilan des forces tient compte :

- R force réactive élastique ;
- F ≠ 0 oscillation amortie ;
- X force d’inertie ;
- P(t) la force due à la vibration où P est sa valeur maximum et  sa
fréquence circulaire.
La prise en compte de la force perturbatrice P(t) donne l’équation d’équilibre
dynamique suivante :
P (t ) P sin t
ÿ  2 y   2 y  
m m
qui présente un second membre. La solution complète de cette équation est la
somme de la solution de l’équation homogène sans le second membre et d’une
solution particulière de l’équation avec le second membre.
Pour ce qui est de l’équation homogène différentielle (1er membre), sa solution
a été examinée au chapitre précédent. Il reste donc à déterminer la solution
particulière de l’équation d’équilibre dynamique.

Posons y1  C1 sin  t  C 2 cos  t ; les constantes C1 et C 2 sont à déterminer.

P sin  t
Substituons cette expression dans ÿ  2 y   2 y  (1)
m
y  C1  cos  t  C 2  sin  t
ÿ  C1  2 sin  t  C 2  2 cos  t
P sin  t
 C1  2 sin  t  C 2  2 cos t  2 C1 cos  t  C 2  sin  t    2 C1 sin  t  C 2 cos  t  
m

ENSTP 2021-2022 18
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

P sin  t
  
C1   2 sin  t  2  cos  t   2 sin  t  C 2   2 cos  t  2  sin  t  w 2 cos  t   m
( 2)

C1 et C 2 seront choisies par identification des équations (1) et (2)

   
C1 sin  t  2   2  2  cos  t  C 2 cos  t  2   2  2   sin  t      P sin  t
m

  
sin  t C1  2   2  2  C 2   P sin  t
m
(1) '

cos  t C 2    C 
1 2
2
 2  0 (2)'

(2)'  C1 
 C2   2
2

2 
(1)'  
C2  2   2 
.  2   2  2  C 2 
P  
2 m
 2  2
  C2 
2
  P
 2  

 2   m
 

  2   2 2  4 2 2 
 C2  P 
 2   m
 
P 2 
C2 
m
 2

 2   2  4 2  2  
C1  
P 2  2  2
 
P 2  2    
.

m  2   2 2  4 2  2
.
2  
.
m  2   2 2  4 2  2   
En substituant C1 et C 2 dans

y1  C1 sin  t  C 2 cos  t , on obtient

P
y1  m

 
 2   2 sin  t

 2  cos  t


 2
 2   
 4 2 2   2   2 2  4 2 2   2
 2 
2 2 
 4  
2

2  2
Posons cos  ' 
 2
 2 2
 4 2 2

puisque la somme des carrées des coefficients sin  t et cos  t est égale à 1 on
peut écrire que :
2 
sin  ' 
 2
 2  2
 4 2 2

alors
P
y1  m .cos  '.sin  t  cos  t. sin  '
 2
 2 
2
 4 2 2

ENSTP 2021-2022 19
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

P
y1  m . sin  t   '
 2
 
2 2
 4 
2 2

1 r 1
2   
m m m 2
en définitive

P
y1  . sin  t   '
2
  2
 4  2 2
1  2  
   4

P
Par conséquent, la solution générale de l’équation ÿ  2 y   2 y  sin  t
m
s’écrira

y  C e  t sin 1 t     y1

On voit que dans le cas envisagé le système participe simultanément à deux

mouvements oscillatoires. Le premier représente le mouvement oscillatoire
propre, et son amplitude et sa phase sont déterminées par les conditions
initiales. Ces oscillations sont amorties et s’évanouissent pratiquement au
bout d’un certain temps. Le second mouvement oscillatoire s’opère avec la
fréquence de la force perturbatrice et son déphasage est  ' . Or, il ne s’amortit
pas, mais dure aussi longtemps que la force perturbatrice. On dit qu’on a des
oscillations forcées. L’amplitude des oscillations forcées A s’écrit :
P
A forcées    .P.
2
  2
 4  2 2
k
1  2   . 2  
   4 m

Le produit P. représente le déplacement qui serait communiqué à la masse

m si on lui appliquait une force statique P . Par conséquent, le coefficient
1   est aussi appelé coefficient 
 dynamique du système 
 2 
2
4 2 2  
1  2  
   4

ENSTP 2021-2022 20
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Commentaire du graphique

Le coefficient  dépend de deux grandeurs : du rapport des fréquences et

de la constante d’amortissement  .
Lorsque   0 , c’est-à-dire en l’absence d’amortissement, la quantité 
devient infinie. Quand les fréquences circulaires des oscillations propres et des
 
oscillations forcées coïncident   1 . Cela signifie qu’en de telles conditions,
 
l’amplitude des oscillations forcées croît indéfiniment.

k
 
m .

Lorsqu’il y a amortissement (pour   0 ),  reste borné, mais sa valeur passe

par un maximum dans la zone où les fréquences coïncident.
Le phénomène d’accroissement de l’amplitude lorsque les fréquences des
oscillations propres et de la force perturbatrice coïncident est dit RESONANCE.
Le phénomène de résonance est d’une importance primordiale dans les calculs
de résistance dynamique. Le fait est que dans la plupart des cas, les lois de
variation des forces perturbatrices ont un caractère périodique. Ainsi, les
parties mobiles mal équilibrées d’un moteur créent des forces variables
périodiques. Un train se mouvant à vitesse constante reçoit des chocs
périodiques à la jonction des rails. Les pièces des appareils installés sur des
fondations vibrantes (dans un avion, une voiture) reçoivent des chocs à la
fréquence de la fondation vibrante. On pourrait multiplier ces exemples. Dans
tous ces cas, il se pose la question du degré de danger des forces
perturbatrices pour le fonctionnement du système élastique et il importe de
ENSTP 2021-2022 21
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

savoir si les vibrations ne sont pas trop excessives, car elles provoquent la
destruction prématurée des appareils.
Un tel problème se résout en premier lieu en comparant les fréquences des
oscillations propres et de la force perturbatrice. Si ces fréquences sont
nettement différentes, on peut être sûr qu’il n’y aura pas résonance et que les
conditions de fonctionnement des éléments élastiques sont favorables. Il est
alors possible de calculer sans difficulté l’amplitude des oscillations forcées
sans connaître a priori la grandeur de la constante d’amortissement  . Les
courbes  diffèrent notament l’une de l’autre seulement dans la zone de
résonance. Dès que la fréquence  est une fois et demie ou deux fois plus
grande ou plus petite que w, on peut poser que les courbes tracées coïncident
pratiquement et la constante d’amortissement  n’intervient pas. On peut la
poser égale à 0, ce qui renforce le coefficient de sécurité. Alors, l’expression
donne :
1
' 
2
1 2

Ainsi pour un système de DDL=1 sous l’action d’une charge dynamique,
l’amplitude de vibration augmente de  ’ fois; de même que l’effort interne, les
contraintes et les déformations augmentent de  ’ fois.
Si le coefficient  ' a été calculé, les contraintes dans les éléments élastiques du
système oscillant se déterminent facilement :
   st  '

 st étant la contrainte qui aurait lieu dans le système lors de l’application

statique de la valeur maximum de la force perturbatrice.

s ' Il résulte de la comparaison des fréquences  et w qu’il y a danger de
résonance, on s’efforce en modifiant la construction, selon les circonstances,
de changer telle ou telle fréquence. Alors, le plus avantageux est de changer

les fréquences de manière à augmenter le rapport pour réduire au

maximum le coefficient  . Le plus simple, à cet effet, est d’assouplir la
suspension, c’est-à-dire de réduire la rigidité des éléments élastiques du
système oscillant. Si, en vertu des circonstances, le constructeur n’a pas la
possibilité de varier les fréquences, lorsqu’il y a danger de résonance, le plus

ENSTP 2021-2022 22
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

simple est d’amortir le système, c’est-à-dire d’installer des dispositifs

élevant la dissipation de l’énergie pendant les oscillations. La constante
d’amortissement  croît alors et l’amplitude dans la zone de résonance est
réduite, le rapport des fréquences étant le même.
1

2
  2
 4  2 2
1  2  
   4

1

 2    
2

1  2    2 2 
     

1

2
   k  
2 2

1  2    . 2 
   m  
1 1
 
2 2
 2   k   2   
2 2

1  2    .  1  2     
    m..       

k

m.

ENSTP 2021-2022 23
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

MESURES POUR CONTRER LE

Chapitre IV :
PHENOMENE DE RESONANCE

I- Caractéristiques physiologiques et autres influences de la vibration

Les oscillations peuvent avoir des répercussions nuisibles sur l’organisme de

l’homme et sur le fonctionnement des installations (équipements).
En effet, l’oscillation des constructions selon des fréquences et amplitudes
données peut entraîner diverses influences sur l’homme.
Sur la base des expériences réalisées en laboratoire, on peut conclure que le
critère de sensibilité de l’homme aux oscillations de petites fréquences (entre 1
et 10 Hz, c’est-à-dire le nombre d’oscillations par seconde) est lié à
l’accélération des oscillations ; et pour les oscillations de grandes fréquences
(nombre d’oscillations / seconde >10) est lié à la vitesse des oscillations.
Le degré d’influence sur l’homme des oscillations harmoniques d’amplitude
allant jusqu’à 1mm est consigné à titre d’exemple dans le tableau ci-dessous.
DEGRE PHYSIOLOGIQUE DE L’INFLUENCE DES OSCILLATIONS
HARMONIQUES
Influence des Fréquences comprises Fréquences comprises
oscillations sur l’homme entre 1 et 10 Hz entre 10 et 100 Hz
Accélération limite Accélération limite
mm / s 2 mm / s 2
- non ressentie 10 0,16
- peu ressentie 40 0,24
- bien ressentie 125 2
- très bien ressentie 400 6,4
(gène)
- nuisible 1.000 16
- très nuisible > 1.000 > 16

Lors des calculs dynamiques, au cas où les différents critères du constructeur

ne sont pas vérifiés à savoir : la résistance, la rigidité ou l’influence des
oscillations sur l’homme, et les équipements, il convient de prendre des

ENSTP 2021-2022 24
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

dispositions ou mesures pour lutter contre la vibration. Ces mesures sont

diverses et effectives pourvu que leur choix soit judicieusement fait. Les plus
importantes sont classées en deux groupes :
- la réduction du phénomène de résonance ;
- isolation contre les vibrations.

II- Réduction du phénomène de résonance

II-1- Installation judicieuse des machines

Exemple : - en installant la machine à côté des appuis de la poutre porteuse on
peut augmenter la fréquence des oscillations propres et cependant réduire le
risque de résonance.
- Lorsqu’une machine crée une force d’inertie horizontale, il faut la disposer
de façon que cette force agisse dans la direction de grande rigidité du
Bâtiment.

II-2- Changement de la rigidité de la construction

Il s’agira d’élever la rigidité de la construction soit en augmentant les
dimensions transversales, soit en réduisant les travées en créant des appuis
intermédiaires. Tout ceci aidera à faire élever la fréquence des oscillations
propres et accentuera le décalage entre la fréquence des oscillations propres et
celle de la force perturbatrice.
Les dispositions contraires à celles indiquées sont possibles dans de rare cas
raison pour laquelle elles ne sont pas recommandées.

II-3- Les amortisseurs

Les amortisseurs dynamiques sont des dispositions complémentaires qu’on
applique aux constructions qui subissent des oscillations forcées dans
l’objectif de neutraliser les réactions émanant de la force perturbatrice.

III- Isolation contre les vibrations

L’isolation contre les vibrations est la manière la plus simple de lutter contre
les oscillations et donc de réduire considérablement les influences nuisibles.
On distingue notamment :

ENSTP 2021-2022 25
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

- l’isolation active contre les vibrations utilisées pour diminuer la force

dynamique transmise par la machine aux structures portantes ;
- l’isolation passive contre les vibrations destinées à isoler les machines et la
construction.
Le calcul de l’isolation active contre les vibrations se conduit sur la base du
schéma suivant.

La masse m sur les amortisseurs possède un degré

de liberté (déplacement vertical)

L’efficacité de l’isolation contre la vibration est appréciée à travers un

coefficient de transmission de force dont l’expression sous sa forme simplifiée
1
s’écrit :  
 
2

  1
 
Pour une bonne efficacité de l’isolation contre la vibration, le rapport entre la
fréquence  de la force perturbatrice et la fréquence des oscillations propres 
ne devrait pas être inférieur à 4. Dans ce cas, la grandeur  ne devrait pas
1
être supérieure à =0,67 et l’efficacité de l’isolation contre la vibration
15
avoisinerait 93,5 %, c’est bien le cas de sous amortissement.
- Pour   (très petite) ;   1 et l’isolation contre la vibration est inutile (pas
nécessaire).

- Pour  1,41 ;   1 et l’isolation contre la vibration peut avoir des effets

nuisibles (par exemple : l’instabilité de la construction). L’amortissement
dans ce cas est au-delà du raisonnable.

- Pour  1 intervient la résonance qui peut être dangereuse pour les

ressorts d’amortissement.
Pour baisser les oscillations pour des fréquences moyennes 6  25 Hz  on
utilise le plus souvent des ressorts isolateurs en acier qui supportent de grands
déplacements statiques ce qui permet de diminuer la fréquence des
oscillations propres de la masse.

ENSTP 2021-2022 26
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Chapitre V : OSCILLATIONS LIBRES DE SYSTEMES A

PLUSIEURS DEGRES DE LIBERTE

Les chapitres précédents parlent des systèmes à un degré de liberté. Leurs

études ont porté essentiellement sur la caractéristique fondamentale des
oscillations à savoir : les fréquences des oscillations propres et forcées.
Dans ce chapitre, nous nous intéresserons à la même étude pour des
systèmes à plusieurs degrés de liberté.

Soit le système suivant :

si l’on néglige la masse de la poutre
et si l’on considère les charges
comme étant concentrées, le
système aura deux degrés de
liberté  deux caractéristiques 1 et

 2 à déterminer. Désignons par y1

le déplacement de la 1ère masse et
par y 2 celui de la 2ème masse. On
déduit les équations suivantes :

  
 y1  B11  m1 y1"  B12  m2 y 2" 

 (I )
y  B
 2   
21  m1 y1  B 22  m 2 y 2
" "


En se servant de l’expression des déplacements sous la forme canonique

y1  B11 X 1  B12 X 2  ..............B1n X n

NB : Au lieu d’une «équation différentielle, nous avons à présent un système de

deux équations comportant y 1 et y 2 . L’application d’une force à la première
masse ayant pour résultat de déplacer la seconde, ces deux équations sont
liées B12  0 .

ENSTP 2021-2022 27
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Les coefficients canoniques B11 ; B12 ; B21 ; B22 peuvent être déterminés soit
par la méthode de MOHR ou de VERECHAGUINE (produit des épures
unitaires).
Choisissons de chercher la solution des équations formées sous la forme
y1  A1 sin t    ; y 2  A2 sin t   
Alors
y1'   A1  cos t    ; y 2'   A2  cos t   
y1"   A1  2 sin t    ; y 2"   A2  2 sin t   

y1"    2 . y1 ; y 2"   2 y 2
Substituons ces fonctions dans les équations suivantes :
  
y1   B11 m1 y1"  B12 m2 y 2"   B11 m1   2 . y1  B12 m2   2 y 2 
  
y 2   B21 m1 y1"  B22 m2 y 2"   B21 m1   2 y1  B22 m2   2 y 2 

   
sin t    B11 m1  2  1 A1  sin t    B12 m2  2 A2  0



21 1 1 22 
sin t    B m  2 A  sin t    B m  2  1 A  0
 2 2 
Simplifions le système d’équation par  2 . sin t   

 1 
 B11 m1  2  A1  B12 m 2 A 2  0
 
 (II)

 B21 m1  A1   B22 m 2  2  A 2  0
1
   
1
Les quantités correspondant à chaque yi s’appellent valeurs propres ou
2
caractéristiques propres. On les désigne dans certaines littératures sous la
1
forme i = 2
i
On obtient après substitutions un système de deux équations homogènes par
rapport à A1 et A2 . (Amplitudes des oscillations en fonction des masses
m1 ; m2 ... ).
Ce système aura une solution non nulle si son déterminant est nul :

ENSTP 2021-2022 28
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

 1  
  B21 m1 B12 m2   0
1
 B11 m1  2   B22 m2  2
    

B12  B21

B11  
m1  2  1 B22 m2  2  1  − B12 m1 B12 m2  2  0
2

B11 B22 m1 m2 4  B11 m1  2  B22 m2  2  1  B12 .B21 m1 m2  4  0

 
 4 m1 m2 B11 .B22  B122   2 B11 m1  B22 m2   1  0

Cette équation détermine deux valeurs pour la fréquence des oscillations

propres :

  2
 B11 m1  B22 m2    B11 m1  B22 m2 2 - 
4 m1 m2 B11 B22  B122 

2 m1 m2 B11 .B22  B 2
12 
4 3 7 3
B11  B22  ; B12  B21  Véréchaguine 
9 EI 18 EI

Posons pour la simplicité de l’étude

6 EI 90 E I
m1  m2  m . On a alors 12  ;  22 
5m 3
5 m 3

1- Analyse et interprétation des résultats

Si le déterminant du système d’équations est nul, l’une des équations se
déduit de l’autre par combinaison linéaire.
En donnant par exemple à la fréquence w une des valeurs indiquées, on
annule le déterminant et l’une des équations (II) peut être considérée comme
conséquente de l’autre.
Il faudrait souligner ici ; qu’en fait nous n’avons pas deux équations, mais une
seule qui permet de déterminer le lien entre A1 et A2 .

Si nous posons la première équation de (II) en y substituant 12 à  2 , on a

pour m1  m2  m A2   A1 et si l’on substitue  22 on aurait A2=-A1.

Par conséquent, les fonctions (II) ont pour expression :
- pour   1

y1  A sin 1 t    ; y 2  A sin  2 t   

ENSTP 2021-2022 29
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

- pour    2

y1  A sin 1 t    ; y 2   A sin  2 t   

Par conséquent, il existe pour un système à deux degrés de liberté, deux
formes d’oscillations.
- Pour la petite fréquence 1 les masses m1 et m2 oscillent en phase étant
donné que les amplitudes ont le même signe.
- Pour la deuxième forme d’oscillations (fréquence  2 ), les amplitudes sont
de signes différents.
Les oscillations s’effectuent en opposition de phases.

2- Représentation des formes d’oscillations

- Oscillation en phase

- Oscillation en opposition de phase

L’expression générale du déterminant pour un système à n degrés de liberté.

1
B11m1  B12 m 2 .......... B1n mn
2
1
D  B21m1 B22 m 2  .......... B2 n m n  0
2
1
Bn1m1 Bn 2 m 2 .......... Bn n m n  2


ENSTP 2021-2022 30
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

OSCILLATIONS DE SYSTEMES A PLUSIEURS

Chapitre VI :
DEGRES DE LIBERTE SOUS L’ACTION DE
SOLLICITATIONS

Sous l’action de sollicitations extérieures, les systèmes à plusieurs degrés de

liberté ont un comportement différent de celui des sollicitations statiques. En
effet, les moments, les efforts tranchants et normaux varient en fonction du
temps. Ainsi la recherche du cas défavorable, demandent un peu plus
d’efforts.
- Si  est la fréquence des forces vibratoires agissant sur un système, alors
la force d’inertie ainsi que   , Q  et  N  atteignent leur plus grande valeur
en un même temps.
- Pour les systèmes à plusieurs degrés de liberté, on peut noter plusieurs
occasions de résonance qui se déclenchent dès que coïncident les
fréquences avec telle ou telle fréquence des oscillations libres desdits
systèmes.
- En gros, l’étude dynamique de ces systèmes consiste à déterminer
l’amplitude des forces internes et des contraintes et à vérifier les systèmes
à la résonance.
 Pour cette étude, on peut se servir soit de la méthode des forces
ou des déplacements.
Soit le système ci-dessous muni de masses ponctuelles sur lesquelles
agissent les forces vibratoires Pk t   Pk sin  t

ENSTP 2021-2022 31
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Le déplacement (déformation) de chaque masse mi au moment t

quelconque est déterminé par y i  Bi1  1  Bi2  2  .....+Bi i i  .....+Bn n  n   i p (I) où

-  1 ,  2 ....., i ...... n désignent les forces d’inertie liées à chaque masse

- Bi1 ; Bi2 ; Bi i désignent les déplacements selon la direction de la force

d’inertie  i émanant des forces unitaires 1 ;  2 , ……..,  i ……  n

-  i p déplacement du point i dû à l’action des forces vibratoires.

Pour des oscillations harmoniques forcées de fréquence  , le déplacement

de la masse mi s’écrit : y i  Ai sin  t ; y i"   Ai  2 sin  t   y i  2

La force d’inertie de masse mi s’exprime par  i   mi y i"  mi y i  2 d’où

i
yi  introduisant yi dans l’équation canonique (I) on obtient
mi  2

Bi 1  1  Bi 2  2  ......  Bi*i  i  ......  Bi n  n   i p  0

1
où Bi*i  Bi i 
mi  2
Ainsi de cette manière on peut composer un système d’équation en fonction
des inconnues qui sont ici des forces d’inertie.
La détermination de ces forces favorise la connaissance des amplitudes.
Pour un système à n équations on aura

 B11*  1  B12  2  ......  B1 n  n   1 p  0


 B21  1  B22  2  ......  B2 n  n   2 p  0
*


..................................................................
 B   B   .......  B *     0
 n1 1 n2 2 nn n np

Le système obtenu est analogue à celui de la méthode des forces (équation

canonique). Cependant, les inconnues ici ne sont pas des réactions superflues
mais des forces d’inertie de masse qu’on peut appréhender aussi bien dans
des systèmes isostatiques que hyperstatiques.
Après la détermination des forces d’inertie  1 ,  2 ..... n , on peut procéder à la

correction des épures unitaires  1 ,  2 ,  n .

L’épure dynamique des moments fléchissant est donnée par l’expression
   1  1   1  2  ......   n  n   p

ENSTP 2021-2022 32
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Selon la procédure que nous connaissons (multiplication des épures et leur

addition).

NB : L’exploitation de la symétrie en dynamique des structures est possible si

la disposition des masses est symétrique.

 Si la disposition des forces vibratoires est antisymétrique et celle des

forces d’inertie est symétrique et antisymétrique, le système d’équation
canonique éclate en deux groupes d’équations non liées. Dans ce cas,
les coefficients canoniques principaux sont déterminés par
1 2
Bi*i  Bi i   Bi i 
1
2
mi  2  mi  2

 Une disposition symétrique des forces vibratoires peut entraîner la

résonance seulement s’il y a coïncidence de sa fréquence avec celle des
oscillations libres symétriques.
 De même, une disposition antisymétrique des forces vibratoires peut
entraîner la résonance seulement s’il y a coïncidence de sa fréquence
avec celle des oscillations libres antisymétriques.

ENSTP 2021-2022 33
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

CALCUL DES SYSTEMES

Chapitre VII : HYPERSTATIQUES

- Méthode de déplacement
La méthode fondamentale de calcul de systèmes hyperstatiques soumis à des
sollicitations dynamiques est celle de la méthode de déplacement.
Dans cette méthode, les inconnues dynamiques de déplacement sont de la
forme Z i t   Z i sin  t

Ex : Pour un système de degré 3, on aura

Z 1 t   Z 1 sin  t ; Z 2 t   Z 2 t   Z 2 sin  t ; Z 3 t   Z 3 sin  t

Les conditions d’équilibre d’un tel système se traduisent par les équations
canoniques suivantes :
r11 z1  r12 z 2  r13 z 3  R1 p  0, 

r21 z1  r22 z 2  r23 z 3  R2 p  0,

r31 z1  r32 z 2  r33 z 3  R3 p  0 

Les inconnues dans ce système d’équations sont les amplitudes des

déplacements vibratoires z1 , z 2 , z 3

ri j - valeur des amplitudes issues de l’action des déplacements unitaires

vibratoires
Ri p - valeur des amplitudes dues à l’effort des forces dynamiques extérieures

La différence fondamentale entre le calcul statique et dynamique par la

méthode des déplacements est la prise en compte des forces d’inerties
ponctuelles ou réparties des différentes masses des éléments portiques.

ENSTP 2021-2022 34
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Ceci est possible grâce aux formules issues des calculs statiques de la
méthode de déplacement plus l’introduction de fonction.
Si la masse des éléments de portiques est uniformément répartie sur toute sa
longueur, alors les amplitudes des réactions dynamiques sont données par les
formules du tableau (schéma+réactions) ; voir cours antérieur.
L’influence de la force d’inertie de la masse répartie est prise en compte grâce
à la fonction  i (u ) ajoutée aux réactions statiques.

m. 3 . 2
L’argument u de la fonction est donné par u  4 ( I ) où
Ei

 - la longueur des éléments

m - masse unitaire d’élément

i - rigidité élémentaire d’élément de portiques
 - fréquence des oscillations forcées ou fréquence des forces vibratoires
(extérieures)
Pour calculer la fréquence des oscillations propres, on se sert toujours de la
formule (I) en remplaçant  par  et en annulant les facteurs Ri p Ri p  0  dans

les équations canoniques.

Ensuite, on annule le déterminant du système.

NB :
 Si les éléments de portiques n’ont pas la même rigidité, alors l’argument
u doit nécessairement changer. L’essentiel c’est de pouvoir uniformiser
l’argument u aux différents éléments.
 La résultante de l’épure dynamique fléchissante est toujours donnée par
M res  M i Z i  M p

 La simplification des calculs passe souvent par la création de système

symétrique et antisymétrique.

ENSTP 2021-2022 35
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

Exercices d’application
N°1
Déterminer les fréquences propres du portique indiqué ci-dessous en
supposant la masse uniformément répartie.

Résolution

La création de système symétrique et antisymétrique réduit de moitié les

coefficients, ainsi, nous obtenons deux équations indépendantes
r11 z1  0 et r22 z 2  0
L’objectif étant de déterminer  i ; Ri p  0
r11 z1  0  r11  0  z1  0 
Cette équation transcendante se
r22 z 2  0  r22  0  z 2  0  résout progressivement jusqu’à
la vérification de l’égalité. On
trouve alors
r11  2 3 i 1 (u )  4 i  2 (u )  4 i  2 (u )  2 i  3 (u )   0
u1  3,34 ; u 2  4,25 ; u 3  4,73
en reportant les ui dans la
31 (u )  8  2 (u )  2  3 (u )  0
formule (I) en changeant  en 

ENSTP 2021-2022 36
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

2
i u EI
on obtient :   u 2
 
m 3
 m

Pour   6 m
2
 3,34  EI EI
1     0,318
 6  m m

2
 4,25  EI EI
2     0,501
 6  m m

2
 4,73  EI EI
3     0,622
 6  m m
r 22  2 3 i  1 ( u )  4 i  2 ( u )  4 i  2 ( u )  2 i  3 ( u )   0

3 1 (u )  8  2 (u )  2  3 (u )  0

 '  3 , 60 ; u 2'  4 , 53

2
 3 , 60  E I E I
 1'     0 , 36
 6  m m

2
 4 , 53  E I E I
 '
2     0 , 57
 6  m m

Suite de l’étude dynamique

Maintenant que la plus petite fréquence propre est connue, on peut évoluer
vers la construction de l’épure dynamique.
Soit   0,6 1 la relation entre  et 1

;  2  0,19 .
EI EI 2 EI
  0,61 1  0,6.0,318  0,19
m m m
alors

m 0,19 
2 EI
m 2 m  6 4 0,19 2  2,61
u  4 6
4

EI EI

ENSTP 2021-2022 37
 COURS D'ANALYSE DYNAMIQUE DES STRUCTURES

à l’aide du tableau on peut déterminer enfin 1 (u )  0,6484 ;  2 (u )  0,8819 ;

3 (u )  1,1813 .
On peut alors résoudre l’équation  z1 
r11 z1  R1 p  0 ; connaissant z1 on construit l’épure rectifiée et l’épure
dynamique des moments fléchissant est donnée par :

   1 Z1   p

ENSTP 2021-2022 38