CHAPITRE 4

REPONSE SISMIQUE DE L'OSCILLATEUR A UN DEGRE DE LIBERTE

4.1 INTRODUCTION

Dans les chapitres précédents, on s'est intéressé à la réponse d'un oscillateur à un degré de
liberté soumis à des sollicitations appliquées directement à la masse M. Une autre forme de
sollicitation est celle constituée par un mouvement du support. C'est en particulier le cas de
sollicitations provenant par le sol comme les vibrations ambiantes ou les sollicitations
sismiques; dans ce cas, le support de l'oscillateur simple représente le sol. Ce chapitre est
plus particulièrement consacré à l'étude de la réponse sismique puisque la réponse aux
vibrations ambiantes relève du domaine de la dynamique stochastique.

On rappelle brièvement qu'une sollicitation sismique provient d'une rupture s'initiant dans le
manteau terrestre; cette rupture donne naissance à des ondes qui se propagent dans le milieu
et atteignent éventuellement la surface du sol où elles se manifestent par une vibration
ressentie dans les trois directions de l'espace. Les appareils enregistreurs, les sismographes,
recueillent ces vibrations sous la forme d'accélérations du sol en fonction du temps. C'est
sous cette forme que la sollicitation du support est définie dans ce chapitre.

4.2 MISE EN EQUATION

La figure 4.1 schématise le système à étudier. La masse M repose sur un support soumis à
une accélération y  (t) fonction du temps. On supposera, dans un premier temps, que

l'accélération y( t ) est connue par sa variation en fonction du temps.

La masse est reliée à son support par un élément développant une force F(u, u ) fonction du
déplacement et de la vitesse relatifs de la masse par rapport à celui ci. On considère que cette
liaison est de type viscoélastique linéaire et s'écrit :

(4.1) F = k u + c u

où u est le déplacement de la masse dans un référentiel lié au support.

55
 u(t)

M
•
F(u, u)

ÿ(t)

ÿ(t)

Figure 4.1 : Excitation sismique de l'oscillateur

L'équation d'équilibre est obtenue par exemple en identifiant les forces s'exerçant sur la
masse (figure 4.2). Ces forces sont constituées par :

• la force de rappel élastique fS = ku

• la force de rappel de l'amortisseur fD = cu
• la force d'inertie fI

y(t) u(t)

fD
fI
fS

X
O

Figure 4.2 : Forces appliquées

La force d'inertie s'exerçant sur la masse est égale au produit de la masse par l'accélération
absolue de celle-ci.

Dénotant v le déplacement de la masse dans un référentiel fixe, la règle de composition des

mouvements s'écrit :

(4.2) v =u+y

où y est le déplacement du support.

56
L'équilibre du système s'obtient en écrivant la nullité de la résultante des efforts, soit :

(4.3) fS + fD + fI = 0

qui en remplaçant les efforts par leurs expressions respectives devient :

(4.4)  + y
M (u  ) + c u + k u = 0

Cette équation peut s'écrire sous la forme réduite introduite au chapitre 2 :

(4.5)  + 2ξω u + ω2 u = − y
u  ( t)

où ω et ξ représentent les caractéristiques de l'oscillateur : pulsation propre et pourcentage

d'amortissement critique.

En comparant l'équation (4.5) à l'équation (2.27) du chapitre 2, on constate que la réponse

sismique de l'oscillateur est analogue à sa réponse lorsqu'il est sollicité directement en
introduisant un chargement équivalent :

(4.6) p eff ( t) = −M y
 (t)

4.3 REPONSE TEMPORELLE DE L'OSCILLATEUR

L'équivalence donnée par l'équation (4.6) indique que les solutions obtenues au chapitre 2
restent applicables en remplaçant p(t) par − My( t ) .

Ainsi les méthodes d'analyse fréquentielle ou temporelle donnent :

• analyse fréquentielle

2i π k
∞
2πk t
u(t) = ∑ ck H(
Tp
(4.7) )e
k =−∞ Tp

où Tp est la période fondamentale de la sollicitation y( t ) rendue périodique et Ck les

coefficients de Fourier de cette fonction périodique; H(•) est la fonction de transfert,
complexe, de l'oscillateur donnée par l'équation (2.82).

• analyse temporelle

L'intégrale de Duhamel donnant la réponse devient :

57
 1 t
(4.8) u( t) = −
ωD ∫ 0
 ( τ) e −ξ ω(t−τ) sin [ ωD(t − τ)] dτ
y

où ωD est la pulsation propre amortie.

Exemple :

Considérons l'accélérogramme enregistré à la station de Lake Hughes (composante E-W)

pendant le séisme de Northridge en 1994 (figure 4.3) et examinons la réponse d'un
oscillateur de pourcentage d'amortissement critique fixe, égal à 2%, et de période propre, T,
variable. Cette réponse, calculée par l'intégrale de Duhamel (4.8) est donnée sur la figure 4.4.

3
Accélération (m/s2)

-1

-2
Temps (s)
-3
0 5 10 15 20

Figure 4.3 :Accélérogramme de Lake Hughes (Northridge 1994)

L'examen de la figure 4.4 montre que l'amplitude, la durée et la forme du signal dépendent
fortement de la période propre de l'oscillateur. Celui-ci répond de façon quasi harmonique
avec une période égale à celle de l'oscillateur. Ce résultat est très général pour une
sollicitation possédant un large contenu fréquentiel et peut être démontré par la théorie des
vibrations aléatoires.

On observe également que parmi les trois oscillateurs, celui avec la période la plus longue
possède la réponse en déplacement la plus élevée. Il ne faut cependant pas faire une règle
générale de cette constatation qui n'est valable que pour une certaine gamme de périodes,
comme on le verra au paragraphe suivant. La comparaison des réponses pour T = 0.5 s et
T = 1 s dément d'ailleurs cette proposition.

58
 0.4

T=2s
0.2 ξ=2%

Déplacement
0

-0.2

temps
-0.4
0 10 20 30 40

0.4 0.4
T=1s T = 0.5 s
0.2 ξ=2% 0.2 ξ=2%
Déplacement

0 0

-0.2 -0.2

temps temps
-0.4 -0.4
0 10 20 30 40 0 10 20 30 40

Figure 4.4 : Influence de la période propre sur la réponse de l'oscillateur

La figure 4.5 présente la réponse de l'oscillateur de période fixe égale à 0.5 s et de

pourcentage d'amortissement critique variable entre 0% et 5%. Pour le système non amorti,
au bout d'une durée correspondant à la phase transitoire de la réponse, la réponse devient
stationnaire avec une amplitude constante, la plus élevée des trois systèmes étudiés. Pour un
amortissement non nul l'amplitude maximale de la réponse diminue et le mouvement
s'atténue d'autant plus rapidement après la phase d'excitation forcée (t>10-15s) que le
pourcentage d'amortissement critique est élevé. On notera également, qu'hormis le cas ξ = 0,
les réponses présentent une certaine similarité.

0.4
T = 0.5 s
0.2 ξ=0%
Déplacement

-0.2

temps
-0.4
0 10 20 30 40

0.4 0.4
T = 0.5 s T = 0.5 s
0.2 ξ=2% 0.2 ξ=5%
Déplacement

0 0

-0.2 -0.2

temps temps
-0.4 -0.4
0 10 20 30 40 0 10 20 30 40

Figure 4.5 : Influence de l'amortissement sur la réponse de l'oscillateur

59
4.4 CALCUL DES EFFORTS

Une fois la réponse de l’oscillateur déterminée, par exemple par l’équation (4.8), l’effort
élastique dans le système s’obtient simplement par :

(4.9) fS (t) = k u(t)

L'expression de la raideur k en fonction de la masse de l’oscillateur et de sa pulsation propre

permet d’écrire (4.9) sous la forme :

(4.10) fS = M ω2 u( t) = M A( t)

où A(t) représente une accélération, appelée pseudo-accélération de l’oscillateur, à ne pas

confondre avec son accélération. Sous la forme (4.10), on reconnaît que l’effort est donné
par le produit de la masse par une accélération.

La pseudo-accélération A(t) de l’oscillateur est peu différente en valeur absolue de

l’accélération absolue de la masse ; la différence provient des forces d’amortissement.
Réécrivant l’équation (4.5) sous la forme :

(4.11)  + y
(u  ) = −ω2 u − 2ξ ω u

on reconnaît l’accélération absolue ( x  ) de la masse et sa pseudo-accélération ω 2 u . Ces

 + y
deux grandeurs ne sont égales, et de signe opposé, que lorsque le pourcentage
d’amortissement critique est nul.

De façon équivalente, on peut définir une pseudo-vitesse relative de la masse V(t), donnée
par :

(4.12) V (t) = ω u(t)

Cette vitesse diffère également de la vitesse relative u de la masse, comme le montre la

relation (4.5) réécrite sous la forme :

1 1
(4.13) u = − ωu −  + y
(u  )
2ξ 2ξ ω

Les pseudo-vitesse et pseudo-accélération sont reliées au déplacement relatif de la masse

par :

(4.14) A = ω V = ω2 u

60
4.5 RÉPONSE MAXIMALE DE L’OSCILLATEUR

L’équation (4.8) fournit la réponse de l’oscillateur à tout instant, à partir de laquelle les
efforts se développant dans le système peuvent être calculés. Pour le dimensionnement d’une
structure, la connaissance de la variation temporelle de l’effort n’est pas obligatoirement
nécessaire et seule la valeur maximale de cette grandeur est requise.

 ( t) permet d’accéder
L’introduction de la notion de spectre de réponse de la sollicitation y
directement à cette grandeur.

Reprenons l’exemple de la figure 4.4 et examinons la réponse de l’oscillateur de période

propre 2s et de pourcentage d’amortissement critique 2% ; cette réponse est donnée sur la
figure 4.6.

Umax = 0.25 dm

0.4
T=2s
ξ=2%
dm))
Déplacement ((dm

0.2

-0.2

-0.4
0 10 20 30 40
Temps

Figure 4.6 : Réponse de l’oscillateur T = 2s, ξ = 2%

La valeur maximale du déplacement de la masse M intervient au temps t = 8.66 s et vaut

2.5 cm ; rappelons que ce déplacement représente le déplacement relatif de la masse par
rapport à son support et non son déplacement absolu.

Il est possible de répéter l’opération précédente pour une collection d’oscillateurs de

périodes propres T et de pourcentages d’amortissement critique variables. Pour chaque
réponse, on ne retient que la valeur maximale de déplacement que l’on notera SD.

(4.15) SD (T, ξ) = max u(T, ξ, t)

t

En reportant les valeurs SD en fonction de la période T pour des valeurs fixées de ξ, on

obtient les courbes de la figure 4.7, appelées spectres de réponse en déplacement, sous-

61
entendu d’un oscillateur simple soumis à la sollicitation de la figure 4.3. Ces courbes,
établies une fois pour toutes pour une sollicitation ÿ(t) donnée, permettent d’obtenir
directement, sans autre calcul, la réponse maximale de tout oscillateur simple de
caractéristiques (T, ξ ) ; ceci est illustré sur la figure 4.7 pour l’oscillateur T = 2s, ξ = 2%
de la figure 4.6.

On a vu au paragraphe 4 que les efforts élastiques dans le système pouvaient être obtenus à
partir du produit de la masse par la pseudo-accélération. Il est donc intéressant d’introduire
un spectre de réponse en pseudo-accélération SA (T, ξ ), défini par :

(4.16) SA (T, ξ) = max ⎡⎣ω2 u(T, ξ, t) ⎤⎦ = ω2 SD (T, ξ)

t

0.12
Déplacement relatif (m)

ξ = 0%
0.1
ξ = 2%
0.08 ξ = 5%
0.06

0.04

0.02 SD = 0.25 dm

0
0 1 2 3 4 5
Période (s)

Figure 4.7 : Spectres de réponse en déplacement

De même, on peut définir un spectre de pseudo-vitesse relative SV (T, ξ ). Ces grandeurs

sont reliées par :

(4.17) SA (T, ξ) = ω SV (T, ξ) = ω2 SD (T, ξ)

Les spectres de la figure 4.7 sont présentés sur la figure 4.8 en termes de pseudo-accélération
calculée par la relation (4.17).

Il convient de réaliser que les figures 4.7 et 4.8 contiennent la même information qui
caractérise la sollicitation ÿ(t) de la figure 4.3.

Une fois connu le spectre de réponse en pseudo-accélération, l’effort élastique maximal

développé dans l’oscillateur est donné par :

(4.18) fS max = MSA (T, ξ)

62
 25

Pseudo--accélération (m/s2) 20 ξ = 0%

15

10
Pseudo

5
ξ = 5%
0
0 0.2 0.4 0.6 0.8 1
Période (s)

Figure 4.8 : Spectres de réponse en pseudo-accélération

Le produit MSA est la force statique équivalente qu’il convient d’appliquer à l’oscillateur
pour l’obtention de l’effort élastique maximal qui s’y développera sous l’action de
l’accélérogramme de la figure 4.3.

4.6 SPECTRES DE RÉPONSE NORMALISÉS

Le paragraphe précédent a montré que si l’on ne s’intéresse qu’à la réponse maximale de la

structure afin de déterminer les efforts maximaux qui s’y développent, la connaissance de la
variation temporelle y  ( t) de l’accélération du sol n’est pas nécessaire. Seul le spectre de
réponse en pseudo-accélération (ou déplacement relatif) est requis.

L’introduction de cette grandeur permet de définir, dans une zone sismique, une sollicitation
enveloppe pour le système, ce que ne permet pas de faire la donnée y  ( t) . En effet, tant du
point de vue du contenu fréquentiel que du point de vue des amplitudes, la variabilité du
mouvement sismique est grande et parmi un ensemble de sollicitations y  i ( t) , une
 j ( t) peut être la plus défavorable pour une structure, alors que pour une autre
sollicitation y
 k ( t) . Par contre, en
structure la sollicitation la plus défavorable sera représentée par y
calculant les spectres de réponse SAi (T, ξ ) de ces sollicitations, il est possible de définir un
spectre enveloppe qui, en vertu de l’équation (4.18), conduira aux efforts les plus élevés. Il
est cependant important de réaliser que ce spectre enveloppe n’est représentatif d’aucun
mouvement sismique y  i ( t) et il est vain de rechercher dans une banque de données de
séismes enregistrés, un enregistrement possédant un spectre de réponse identique, ou même
voisin, d’un spectre enveloppe.

63
Par contre, la définition d’un spectre enveloppe est très utile pour les règlements de
construction parasismique. Ces spectres y sont définis à partir de traitements statistiques des
spectres d’enregistrement réels et possèdent généralement tous des formes identiques
(figure 4.9).

10

ξ =2%
Pseudo--accélération (m/s2)

8
ξ =5%

6 ξ =10%

ξ =20%
4

2
Pseudo

TB TC TD
0
0 1 2 3 4 5
Période (s)

Figure 4.9 : Spectres de réponse normalisés

• Pour les faibles périodes, la pseudo-accélération croît linéairement jusqu’à une

valeur d’amplification maximale atteinte pour une période TB ;

• La pseudo-accélération est constante et maximale pour des périodes comprises entre

TB et au plus TC ;

• Entre les périodes TC et TD, la pseudo-vitesse relative est constante, ce qui signifie
que la pseudo-accélération décroît en 1/T (équation 4.17) ;

• Au delà de TD, le déplacement relatif est constant, ce qui signifie que la pseudo-
accélération décroît en 1/T2 (équation 4.17).

Il est intéressant d’examiner les deux extrémités du spectre, à période nulle et à période
infinie :

• A période nulle, la pseudo-accélération est égale à l’accélération maximale du

support, comme le montre l’équation (4.5) lorsque ω = ∞ ; le déplacement
relatif u devient nul. Quelle que soit la valeur de l’amortissement, tous les
 ( t) ⎤⎦ ;
spectres se raccordent à période nulle à l’accélération maximale max ⎡⎣ y

• A période infinie, le déplacement relatif SD (∞, ξ) tend vers le déplacement

maximal du support max ⎡⎣ y ( t) ⎤⎦ . Ce résultat s’explique aisément : un
oscillateur de période infinie correspond à un système avec une raideur nulle :

64
 tout se passe comme si la masse n’était plus connectée au support; par suite, lors
de l’excitation de celui-ci, elle reste immobile. Son mouvement relatif par
rapport au support est donc égal au mouvement du support. En pratique,
l’examen des spectres de réponse d’enregistrements réels montre que la réponse
devient égale au déplacement du sol pour des périodes TF de l’ordre de 10s.

0.40
Déplacement relatif (m)

ξ =2%

0.30 ξ =5%
ξ =10%
0.20 ξ =20%

0.10

TC TD TE TF
0.00
0 2 4 6 8 10 12
Période (s)

Figure 4.10 : Spectres de réponse en déplacement

Les spectres en pseudo-accélérations de la figure 4.9 peuvent être tracés en termes

de déplacement relatif par la transformation (4.17) (figure 4.10). Sous cette forme,
les considérations générales développées au paragraphe 3.0, relatives à la réponse de
divers oscillateurs, s’explicitent ; on constate en particulier que les oscillateurs de
périodes élevées, comprises entre TD et TE présentent les valeurs maximales.

65
66