EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 1.

Analyse sismique
Méthodes à disposition:

action statique dynamique

structure

élastique Forces de Spectre de

Remplacement Réponse

non-linéaire “Push-Over” Non-Linéaire

Dynamique

Equivalence pour la détermination des déplacements relatifs:

xa(t)
x(t) x(t)
..
m·xg(t)
m m

k k
xg(t)

..
xg(t)
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 2.

Systèmes à plusieurs degrés de liberté

Mouvement de la fondation (séismes)

..
x5(t) m5 xg(t) x5(t)
m5 m5
..
x4(t) m4 xg(t) x4(t)
m4 m4
..
x3(t) m3 xg(t) x3(t)
m3 m3
..
x2(t) m2 xg(t) x2(t)
m2 m2
..
x1(t) m1 xg(t) x1(t)
m1 m1

xg(t)

..
xg(t)

Equation du mouvement pour chaque degré de liberté:

m 1 ⋅ x˙˙1 + c ⋅ x˙1 + k 11 ⋅ x 1 + k 12 ⋅ x 2 + k 13 ⋅ x 3 + k 14 ⋅ x 4 + k 15 ⋅ x 5 = – m 1 ⋅ x˙˙g ( t )

m 2 ⋅ x˙˙2 + c ⋅ x˙2 + k 21 ⋅ x 1 + k 22 ⋅ x 2 + k 23 ⋅ x 3 + k 24 ⋅ x 4 + k 25 ⋅ x 5 = – m 2 ⋅ x˙˙g ( t )

m 3 ⋅ x˙˙3 + c ⋅ x˙3 + k 31 ⋅ x 1 + k 32 ⋅ x 2 + k 33 ⋅ x 3 + k 34 ⋅ x 4 + k 35 ⋅ x 5 = – m 3 ⋅ x˙˙g ( t )

m 4 ⋅ x˙˙4 + c ⋅ x˙4 + k 41 ⋅ x 1 + k 42 ⋅ x 2 + k 43 ⋅ x 3 + k 44 ⋅ x 4 + k 45 ⋅ x 5 = – m 4 ⋅ x˙˙g ( t )

m 5 ⋅ x˙˙5 + c ⋅ x˙5 + k 51 ⋅ x 1 + k 52 ⋅ x 2 + k 53 ⋅ x 3 + k 54 ⋅ x 4 + k 55 ⋅ x 5 = – m 5 ⋅ x˙˙g ( t )

Sous forme matricielle:

M ⋅ x˙˙ + C ⋅ x˙ + K ⋅ x = – M ⋅ e x ⋅ x˙˙g ( t ) (5.7)

(la numérotation des équations correspond à celle du polycopié)

EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 3.

où la matrice des masses doit être multipliée par le vecteur direction ex:

m1 0 0 0 0 m1
1
0 m2 0 0 0 1 m2
–M ⋅ ex = – 0 0 m 0 0 ⋅ 1 = – m3
3
0 0 0 m4 0 1 m4
1
0 0 0 0 m5 m5

et C est la matrice d’amortissement:

c 0 0 0 0
0 c 0 0 0
C = 0 0 c 0 0
0 0 0 c 0
0 0 0 0 c

Le système d’équations couplées (5.7) peut à nouveau être découplé et transformé en un

système d’oscillateurs simples (si M et K sont constantes dans le temps).
Les modes propres s’obtiennent en annulant le déterminant du système (solutions non
triviales):
2
K – ωn ⋅ M = 0 (5.8)

Les déplacements relatifs (x j ) s’expriment en coordonnées modales (z n ) par le

changement de variables suivant:

x1 A 11 A 12 A 13 A 14 A 15 z1
x2 A 21 A 22 A 23 A 24 A 25 z2 5
x = x3 = A 31 A 32 A 33 A 34 A 35 ⋅ z3 = A⋅z = ∑ An ⋅ zn (5.9)
x4 A 41 A 42 A 43 A 44 A 45 z4 n=1

x5 A 51 A 52 A 53 A 54 A 55 z5

Où Αn sont les vecteurs propres, z les coordonnées modales et Α la matrice des vecteurs
modaux.
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 4.

En tenant compte du changement de variables (5.9) et en prémultipliant par ΑT, on

obtient:
T T T T
A ⋅ M ⋅ A ⋅ z˙˙+ A ⋅ C ⋅ A ⋅ z˙ + A ⋅ K ⋅ A ⋅ z = – A ⋅ M ⋅ e x ⋅ x˙˙g ( t ) (5.10)
T
avec – A ⋅ M ⋅ e x ⋅ x˙˙g ( t ) = – r ⋅ x˙˙g ( t )
où r est le vecteur des facteurs de participation:

A 11 A 21 A 31 A 41 A 51 m1 r1
A 12 A 22 A 32 A 42 A 52 m2 r2
T
r = A ⋅ M ⋅ ex = A 13 A 23 A 33 A 43 A 53 ⋅ m3 = r3
A 14 A 24 A 34 A 44 A 54 m4 r4
A 15 A 25 A 35 A 45 A 55 m5 r5

r1 A 11 ⋅ m 1 + A 21 ⋅ m 2 + A 31 ⋅ m 3 + A 41 ⋅ m 4 + A 51 ⋅ m 5 ∑ Aj1 ⋅ mj
r2 A 12 ⋅ m 1 + A 22 ⋅ m 2 + A 32 ⋅ m 3 + A 42 ⋅ m 4 + A 52 ⋅ m 5 ∑ Aj2 ⋅ mj
r 3 = A 13 ⋅ m 1 + A 23 ⋅ m 2 + A 33 ⋅ m 3 + A 43 ⋅ m 4 + A 53 ⋅ m 5 = ∑ Aj3 ⋅ mj
r4 A 14 ⋅ m 1 + A 24 ⋅ m 2 + A 34 ⋅ m 3 + A 44 ⋅ m 4 + A 54 ⋅ m 5
∑ Aj4 ⋅ mj
A 15 ⋅ m 1 + A 25 ⋅ m 2 + A 35 ⋅ m 3 + A 45 ⋅ m 4 + A 55 ⋅ m 5
r5
∑ Aj5 ⋅ mj
En utilisant des grandeurs généralisées, comme la masse généralisée:

2 2 2 2 2
∑ Aj1 ⋅ mj
2
A 11 ⋅ m 1 + A 21 ⋅ m 2 + A 31 ⋅ m 3 + A 41 ⋅ m 4 + A 51 ⋅ m 5
m 1∗
2 2 2 2 2
∑ j2 ⋅ mj
2
m 2∗ A 12 ⋅ m 1 + A 22 ⋅ m 2 + A 32 ⋅ m 3 + A 42 ⋅ m 4 + A 52 ⋅ m 5 A

m 3∗ = A 2 ⋅ m + A 2 ⋅ m + A 2 ⋅ m + A 2 ⋅ m + A 2 ⋅ m =
∑ j3 ⋅ mj
2
13 1 23 2 33 3 43 4 53 5 A
m 4∗ 2 2 2 2 2
∑ Aj4 ⋅ mj
2
A 14 ⋅ m 1 + A 24 ⋅ m 2 + A 34 ⋅ m 3 + A 44 ⋅ m 4 + A 54 ⋅ m 5
m 5∗ 2 2 2 2 2
∑ Aj5 ⋅ mj
2
A 15 ⋅ m 1 + A 25 ⋅ m 2 + A 35 ⋅ m 3 + A 45 ⋅ m 4 + A 55 ⋅ m 5

Le système se simplifie:
∗ T
M∗ ⋅ z˙˙+ C∗ ⋅ z˙ + K ⋅ z = – A ⋅ M ⋅ e x ⋅ x˙˙g ( t ) = – r ⋅ x˙˙g ( t ) (5.11)
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 5.

Il faut noter que les grandeurs généralisées et les facteurs de participation dépendent de
la normalisation des vecteurs propres Αn.
La nième ligne de l’équation d’équilibre en coordonnées modales s’écrit:

m n∗ ⋅ z˙˙n + c n∗ ⋅ z˙n + k n∗ ⋅ z n = – r n ⋅ x˙˙g ( t ) (5.12)

Après division par mn*, cette nième ligne de l’équation d’équilibre en coordonnées
modales devient:

2 r n ˙˙
z˙˙n + 2 ⋅ ζ n ⋅ ω n ⋅ z˙n + ω n ⋅ z n = – ---------
- ⋅ xg ( t ) (5.13)
m n∗

A l’exception de la modification de l’amplitude de l’accélération du sol, c’est à nouveau

l’équation bien connue d’un oscillateur simple, dans laquelle:

2 k n∗
ω n = ---------
-
m ∗ n

Le processus de décomposition modale peut être illustré de la manière suivante: (Il faut
bien noter l’influence des facteurs de participation qui modifient l’amplitude de
l’accélération du sol de chaque mode)

x5(t)
m5
z1(t)
x4(t) z2(t)
m4
z3(t)
x3(t)
m3 m1* z4(t) z5(t)
m2*
x2(t)
m2 m3*
k1* m4*
x1(t)
k2* m5*
m1
k3*
k4* k5*

.. r1 .. r2 .. r3 .. r4 .. r5 ..
xg(t) xg(t) xg(t) xg(t) xg(t) xg(t)
m1* m2* m3* m4* m5*

La réponse modale s’obtient en résolvant (5.13) numériquement par une des méthodes
habituelles (Interpolation de l’excitation, Différence Centrale, Newmark, etc...).
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 6.

Les déplacements relatifs sont ensuite déterminés par la superposition des réponses
modales selon (5.9):

x1 A 11 A 12 A 13 A 14 A 15
x2 A 21 A 22 A 23 A 24 A 25
x3 = A ⋅ z1 ( t ) + A ⋅ z2 ( t ) + A ⋅ z3 ( t ) + A ⋅ z4 ( t ) + A ⋅ z5 ( t )
31 32 33 34 35
x4 A 41 A 42 A 43 A 44 A 45
x5 A 51 A 52 A 53 A 54 A 55

A51 A52 A53 A54 A55

A41 A42 A43 A44 A45

A31 A32 A33 A34 A35

A21 A22 A23 A24 A25

A11 A12 A13 A14 A15

z1(t) z2(t) z3(t) z4(t) z5(t)

25 25 25 25 25

20 20 20 20 20

15 15 15 15 15
temps [s]

10 10 10 10 10

5 5 5 5 5

0 0 0 0 0
-400 -200 0 200 400 -100 0 100 -100 0 100 -6 0 6 -6 0 6
déplacement [mm] déplacement [mm] déplacement [mm] déplacement [mm] déplacement [mm]
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 7.

Une grandeur permet de quantifier à priori la contribution de chaque mode à la réponse

totale, il s’agit de la masse modale qui est définie comme suit:
N 2
⎛ ⎞
⎛ rn ⎞ 2
⎜
⎝
∑
A jn ⋅ m j⎟
⎠
m mod, n = ⎜ ----------⎟ ⋅ m n∗ = --------------------------------------
j=1 (5.14)
⎝ m n∗⎠ N

∑
2
A jn ⋅ m j
j=1

Contrairement au facteur de participation et à la masse généralisée, la valeur de la masse

modale ne dépend pas de la normalisation des modes. Elle correspond à l’effort tranchant
à la base de l’oscillateur multiple dans chaque mode. De plus, la somme des masses
modales est égale à la somme des masses:
N N

∑ mmod, n = ∑ mn (5.15)
n=1 n=1
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 8.

Spectre de réponse
Outil de l’ingénieur:

élastique (coefficient d'amortissement fixé: ζ)

accélération [m/s ]
2

10
2
5 S =2.31 m/s
e
0
-5
-10
0 5 10 15 20
accélération [m/s ]
2

10
5
0
-5 2
S =7.00 m/s
e
-10 spectre de réponse de l'accélération
0 5 10 15 20

10
ζ=5%
S e [m/s2 ]

0
accélération [m/s ]

0.01 0.1 1 10
2

10 période [s]
5
0
-5 2
S =11.33 m/s
-10 e

0 5 10 15 20

séisme: accélérogramme
4
accélération [m/s ]
2

3 2
2 a =3.66 m/s
max
1
0
-1 2
-2 a = -2.62 m/s
min
-3
0 5 10 15 20
temps [s]
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 9.

Spectres de réponse
Définitions:
t
1
x··g(τ) ·e – ζω(t-τ) · sin ω(t-τ) dτ
Su(ω,ζ) = max |x(t)| ≈ ω ∫ max
0
t
· x··g(τ) ·e – ζω(t-τ) · cos ω(t-τ) dτ
Sv(ω,ζ) = max |x(t)| ≈ – ∫ max
0
t
·· ·· | ≈ ω x··g(τ) ·e – ζω(t-τ) · sin ω(t-τ) dτ ≈ ω Su
Se(ω,ζ) = max |x+x g ∫ max
2

0
t
1
x··g(τ) · e – ζω(t-τ) · sin ω(t-τ) dτ ≈ ω Se ≈ ω Su
Spv(ω,ζ) = – ∫ max
0
Cette dernière relation permet de réunir trois spectres (Su, Spv, Se) en une seule courbe
dans une double échelle logarithmique (pas pratique).

Ne pas confondre avec le spectre de Fourier:

4
3 2
accélération [m/s ]

a =3.66 m/s
2

2 max

1
0
-1 2
-2 a = -2.62 m/s
min
-3
0 5 10 15 temps [s] 20

Spectre de Fourier Spectre de Réponse

ζ=5%
amplitude [m/s /Hz]

10
2

S e [m/s2 ]

0 0
0.1 1 10 0.1 1 10
fréquence [Hz] fréquence [Hz]
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 10.

Spectres SIA 261

Formulation:

Paramètres en fonction de la classe de sol de fondation:

EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 11.

Spectres SIA 261

classes de sols de fondation
A B C D E
0m
terrain meuble

couche alluviale
h<5 m de classe de sol
C ou D
10 m

sables et sables et sables fins

graviers graviers non consolidés,
cimentés, norm. limons et argiles
sol préconsolidé consolidés,
matériau
morainique
20 m

couche de
classe de sol
A ou B

30 m

h>30 m h>30 m h>30 m 5 m<h<30 m

Graphiquement en fonction de la classe de sol de fondation:

4
sol A
3 sol B
Se /agd [-]

sol C
2
sol D
1
sol E

0
0.01 0.1 1 10
période [s]
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 12.

Méthode du spectre de réponse

Pour le dimensionnement par exemple, seule la réponse maximum est intéressante. Elle
peut être estimée rapidement à l’aide du spectre de réponse.

S e,2
3
Se [m/s2]

2 S e,3
S e,4 S e,1
S e,5

0
0.01 T5 T4 T3 0.1 T2 1 T1 10
période [s]
La réponse modale maximale s’obtient en utilisant le spectre de réponse:
rn rn
z n, max = ---------
- ⋅ S u, n ( ω n, ζ ) = --------------------- ⋅ S e, n ( ω n, ζ ) (5.16)
m n∗ 2
ωn ⋅ mn ∗

Les déplacements relatifs maximaux de chaque modes sont donnés par:

A 1n x 1n, max
A 2n x 2n, max
x n, max = A n ⋅ z n, max = A 3n ⋅ z n, max = x 3n, max (5.17)
A 4n x 4n, max
A 5n x 5n, max

De la même manière, la réponse totale s’obtient en superposant les réponses des

différents modes. Toutefois, la simple addition des réponses maximum conduit à des
valeurs trop défavorables car les maximum ne coïncident généralement pas dans le
temps. On utilise plutôt la règle de la racine de la somme des carrés (SRSS: square root
of the sum of the squares):
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 13.

N
2
x j, max, tot = ∑ x jn, max (5.18)
n=1

Les forces maximum ne peuvent pas être directement calculées avec les déplacements
relatifs maximum (Fmax,tot ≠ K·xmax,tot!!!). Il faut déterminer les forces maximales dans
chaque modes à partir de la matrice de rigidité,

F 1n, max x 1n, max k 11 k 12 k 13 k 14 k 15 A 1n

F 2n, max x 2n, max k 21 k 22 k 23 k 24 k 25 A 2n
F 3n, max = K ⋅ x 3n, max = k 31 k 32 k 33 k 34 k 35 ⋅ A 3n ⋅ z n, max
F 4n, max x 4n, max k 41 k 42 k 43 k 44 k 45 A 4n
F 5n, max x 5n, max k 51 k 52 k 53 k 54 k 55 A 5n

rn
- ⋅ S e, n ( ω n, ζ ) (5.19)
F n, max = K ⋅ x n, max = K ⋅ A n ⋅ z n, max = K ⋅ A n ⋅ -----------------------
2
ω ⋅m ∗ n n

pour les superposer ensuite selon la règle SRSS de l’équation (5.18).

N
2
F j, max, tot = ∑ F jn, max
n=1

F51,max 2 F52,max 2 F53,max 2 F5,max,tot

F41,max F42,max F43,max F4,max,tot

F31,max F32,max F33,max F3,max,tot

F21,max F22,max F23,max ... F2,max,tot

F11,max F12,max F13,max F1,max,tot

2
∑ Fj1,max ∑ Fj2,max ∑ Fj3,max ∑ (∑ Fjn,max) ∑ Fj,max,tot
EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 14.

Systèmes à plusieurs degrés de liberté

Matrice de flexibilité d’une console:
Fa2(3L-a)
Calcul des déformations avec la matrice de d=
6EI
flexibilité:

d=^
f·F F
L
a
si les hauteurs d’étage sont constantes:
^f h3 2
i,j = i (3j-i), j ≥ i
6EI

Matrice de flexibilité: Matrice de rigidité:

2 5 8 11 14

5 16 28 40 52
h3
f^ = 8 28 54 81 108 K = f^ -1
6EI
11 40 81 128 176

14 52 108 176 250

dj
Utilité pour le quotient de Rayleigh

F3=3f
n
h
1 ∑ Fj dj F2=2f
j=1
f =
2

√ n
∑ mj dj2
j=1
F1=f
h H

h
mj: masse de l’étage j
dj: déplacement de l’étage j
Fj: force de remplacement agissant à l’étage j
EPFL-ENAC-SGC Génie Parasismique semestre d’été 2007 P. Lestuzzi 15.

Méthode du spectre de réponse

Exemple numérique, bâtiment en béton armé de cinq étages:
6m

masse par étage: 300 t

4
CO I
b=0.3 m

4
spectre: SIA zone 3b C25/30

20 m
4
classe de sol C

4
4
4

8m
4
Matrice de flexibilité: Matrice des masses:

2 5 8 11 14 300 0 0 0 0

5 16 28 40 52 0 300 0 0 0
3
h
f^ = 8 28 54 81 108 M= 0 0 300 0 0
6EI
11 40 81 128 176 0 0 0 300 0

14 52 108 176 250 0 0 0 0 300

Matrice de rigidité:

3.1381 -1.9834 0.7956 -0.1989 0.0331

-1.9834 2.4420 -1.7845 0.6961 -0.1160

6EI
K = f^ -1 = 3 0.7956 -1.7845 2.3425 -1.5856 0.4309
h
-0.1989 0.6961 -1.5856 1.6464 -0.6077

0.0331 -0.1160 0.4309 -0.6077 0.2680

EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 16.

Systèmes à plusieurs degrés de liberté

Caractéristiques dynamiques du bâtiment de cinq étages:

m = 300 t, CO I ( γf = 1.0), classe de sols C

EI = 30% EI0 m
4m
1 2 3 4 5 m
f [Hz] 0.87 5.57 15.8 30.5 45.5 4m

m* [kNs2/m] 531 891 1011 714 868 m

4m
L = r/m* [–] 1.38 0.59 0.32 0.26 -0.14
m
Mmod [kNs2/m] 1018 310 105 49 18
4m
Se [m/s 2] 2.41 4.32 2.71 2.29 2.14
m
4m
m*: masse généralisée
r: facteur de participation
Mmod: masse modale

Modes propres normalisés:

1.00 -0.90 0.63 -0.31 -0.14

0.72 0.29 -0.97 0.85 0.52

0.45 1.00 -0.47 -0.72 -0.88

0.22 0.96 1.00 -0.21 1.00

0.06 0.39 0.90 1.00 -0.91

1er 2ème 3ème 4ème 5ème

EPFL-ENAC-SGC Génie parasismique semestre d’été 2007 P. Lestuzzi 17.

Bâtiment régulier de cinq étages

Forces maximales dans chaque mode [kN]:

1000.5 -683.9 166.2 -55.2 13.1

718.1 220.5 -255.7 154.6 -47.8

451.0 762.9 -122.8 -129.8 81.2

222.3 734.9 262.5 -37.2 -92.3

61.2 301.2 235.1 180.9 84.2

∑=2453 ∑=1336 ∑=285 ∑=113 ∑=38

Moments maximaux dans chaque modes [MNm]:

4.00 -2.74 0.66 -0.22 0.05

10.9 -4.59 0.31 0.18 -0.09

19.6 -3.39 -0.54 0.05 0.10

29.1 0.75 -0.34 -0.22 -0.08

38.9 6.09 0.80 0.24 0.07

1er 2ème 3ème 4ème 5ème