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 ex xg ( t ) (5.7)
(la numrotation des quations correspond celle du polycopi)
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
2.
o la matrice des masses doit tre multiplie par le vecteur direction ex: m1 0 0 m1 1 0 0 m2 1 0 0 1 = m3 1 m4 0 m4 1 m5 0 m5 0 0
0 m2 0 M ex = 0 0 m 3 0 0 0 0 0 0
et C est la matrice damortissement: c 0 C = 0 0 0 0 c 0 0 0 0 0 c 0 0 0 0 0 c 0 0 0 0 0 c
Le systme dquations couples (5.7) peut nouveau tre dcoupl et transform en un systme doscillateurs simples (si M et K sont constantes dans le temps). Les modes propres sobtiennent en annulant le dterminant du systme (solutions non triviales): K n M = 0 2
(5.8)
Les dplacements relatifs (x j ) sexpriment en coordonnes modales (z n ) par le changement de variables suivant: x1 x2 x = x3 x4 x5 = A 11 A 12 A 13 A 14 A 15 A 21 A 22 A 23 A 24 A 25 A 31 A 32 A 33 A 34 A 35 A 41 A 42 A 43 A 44 A 45 A 51 A 52 A 53 A 54 A 55 z1 z2 z3 z4 z5 = Az =
n=1
An zn
(5.9)
O n sont les vecteurs propres, z les coordonnes modales et la matrice des vecteurs modaux.
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
3.
En tenant compte du changement de variables (5.9) et en prmultipliant par T, on obtient:
T T T T A M A z + A C A z + A K A z = A M ex xg ( t )
(5.10)
avec o r est le vecteur des facteurs de participation: A 11 A 21 A 31 A 41 A 51 A 12 A 22 A 32 A 42 A 52 A 13 A 23 A 33 A 43 A 53 A 14 A 24 A 34 A 44 A 54 A 15 A 25 A 35 A 45 A 55
T A M ex xg ( t ) = r xg ( t )
m1 m2 m3 m4 m5
r1 r2 = r3 r4 r5
r = A M ex =
r1 r2
A 11 m 1 + A 21 m 2 + A 31 m 3 + A 41 m 4 + A 51 m 5 A 12 m 1 + A 22 m 2 + A 32 m 3 + A 42 m 4 + A 52 m 5
r 3 = A 13 m 1 + A 23 m 2 + A 33 m 3 + A 43 m 4 + A 53 m 5 = A 14 m 1 + A 24 m 2 + A 34 m 3 + A 44 m 4 + A 54 m 5 r4 r5 A 15 m 1 + A 25 m 2 + A 35 m 3 + A 45 m 4 + A 55 m 5
Aj1 mj Aj2 mj Aj3 mj Aj4 mj Aj5 mj Aj1 mj 2 A j2 m j 2 A j3 m j 2 Aj4 mj 2 Aj5 mj
2
En utilisant des grandeurs gnralises, comme la masse gnralise: m 1 m 2 A 11 m 1 + A 21 m 2 + A 31 m 3 + A 41 m 4 + A 51 m 5 A 12 m 1 + A 22 m 2 + A 32 m 3 + A 42 m 4 + A 52 m 5
2 2 2 2 2 2 2 2 2 2
m 3 = A 2 m + A 2 m + A 2 m + A 2 m + A 2 m = 13 1 23 2 33 3 43 4 53 5 2 2 2 2 2 m 4 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 A 15 m 1 + A 25 m 2 + A 35 m 3 + A 45 m 4 + A 55 m 5 Le systme 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
Gnie parasismique
semestre dt 2005
P. Lestuzzi
4.
Il faut noter que les grandeurs gnralises et les facteurs de participation dpendent de la normalisation des vecteurs propres n. La nime ligne de lquation dquilibre en coordonnes modales scrit: m n z n + c n z n + k n z n = r n x g ( t ) (5.12)
Aprs division par mn*, cette nime ligne de lquation dquilibre en coordonnes modales devient: r n 2 z n + 2 n n z n + n z n = --------- x g ( t ) m n (5.13)
A lexception de la modification de lamplitude de lacclration du sol, cest nouveau lquation bien connue dun oscillateur simple, dans laquelle: k n 2 n = --------m n
Le processus de dcomposition modale peut tre illustr de la manire suivante: (Il faut bien noter linfluence des facteurs de participation qui modifient lamplitude de lacclration du sol de chaque mode)
La rponse modale sobtient en rsolvant (5.13) numriquement par une des mthodes habituelles (Interpolation de lexcitation, Diffrence Centrale, Newmark, etc...).
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
5.
Les dplacements relatifs sont ensuite dtermins par la superposition des rponses modales selon (5.9): x1 x2 x3 x4 x5 A 11 A 21 A 12 A 22 A 13 A 23 A 14 A 24 A 15 A 25
= A z1 ( t ) + A z2 ( t ) + A z3 ( t ) + A z4 ( t ) + A z5 ( t ) 31 32 33 34 35 A 41 A 42 A 43 A 44 A 45 A 51 A 52 A 53 A 54 A 55
A51 A41 A31 A21 A11 A42 A32 A22 A12
A52
A53 A44
A54 A43 A33 A45 A34 A24 A14 A25
A55
A35
A23 A13
A15
z1(t)
25
z2(t)
25
z3(t)
25
z4(t)
25
z5(t)
25
20
20
20
20
20
10
10
10
10
10
0 -400 -200 0 200 400
-100
0 100
-100
0 100
0 -6 0 6 -6 0 6
dplacement [mm]
dplacement [mm]
dplacement [mm]
dplacement [mm]
dplacement [mm]
temps [s]
15
15
15
15
15
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
6.
Une grandeur permet de quantifier priori la contribution de chaque mode la rponse totale, il sagit de la masse modale qui est dfinie comme suit: A jn m j rn 2 j=1 m mod, n = --------- m n = -------------------------------------N m n 2 A jn m j
(5.14)
j=1
Contrairement au facteur de participation et la masse gnralise, la valeur de la masse modale ne dpend pas de la normalisation des modes. Elle correspond leffort tranchant la base de loscillateur multiple dans chaque mode. De plus, la somme des masses modales est gale la somme des masses:
Cette dernire relation permet de runir 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 2
acclration [m/s ]
3 2 1 0 -1 -2 -3 0 5
max
=3.66 m/s
min
= -2.62 m/s 10
15
temps [s]
20
Spectre de Fourier amplitude [m/s /Hz] 10 1 2
Spectre de Rponse =5% S e [m/s2 ]
0 0.1
10
0 0.1
10
frquence [Hz]
frquence [Hz]
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
9.
Spectre de dimensionnement Par exemple pour la Suisse, SIA 160, dition 1989.
Zones dala sismique:
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
10.
Spectres de dimensionnement En Suisse nouvelle norme SIA 261 partir de 2003 SIA 160/89 4 3 4 3
SIA 261
ah/g/as [-]
2 1 0 0.1 4 3
Se /agd [-] 1 10 100
2 1 0 0.1 4 3
frquence [Hz]
frquence [Hz]
10
100
ah/g/as [-]
2 1 0 0.01
Se /agd [-] 0.1 1 10
2 1 0 0.01
priode [s]
0.1
priode [s]
10
sols compacts sols semi-compacts
classe de sol A classe de sol B
classe de sol C classe de sol D classe de sol E
Zones de risque sismique:
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
11.
Spectres SIA 261
Formulation:
Paramtres en fonction de la classe de sol de fondation:
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
12.
Spectres SIA 261 A
terrain meuble
classes de sols de fondation B C D
0m
h<5 m
couche alluviale de classe de sol C ou D
10 m sables et graviers ciments, sol prconsolid sables et graviers norm. consolids, matriau morainique sables fins non consolids, limons et argiles
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 sol C sol D 1 0 0.01 sol E
Se /agd [-]
0.1
priode [s]
10
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
13.
Mthode du spectre de rponse
Pour le dimensionnement par exemple, seule la rponse maximum est intressante. Elle peut tre estime rapidement laide du spectre de rponse. 4
S e,2
Se [m/s2]
S e,5 1
S e,4
S e,3
S e,1
0 0.01
T5
T4
T3
0.1
T2
priode [s]
T1
10
La rponse modale maximale sobtient en utilisant le spectre de rponse: rn rn z n, max = --------- S u, n ( n, ) = --------------------- S e, n ( n, ) 2 m n n m n (5.16)
Les dplacements relatifs maximaux de chaque modes sont donns par: A 1n x n, max = x 1n, max x 2n, max (5.17)
n=1
An zn, max =
A 2n
A 3n z n, max = x 3n, max n=1 A x 4n, max 4n A 5n x 5n, max
De la mme manire, la rponse totale sobtient en superposant les rponses des diffrents modes. Toutefois, la simple addition des rponses maximum conduit des valeurs trop dfavorables car les maximum ne concident gnralement pas dans le temps. On utilise plutt la rgle de la racine de la somme des carrs (SRSS: square root of the sum of the squares):
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
14.
x j, max, tot =
n=1
x jn, max
(5.18)
Les forces maximum ne peuvent pas tre directement calcules avec les dplacements relatifs maximum (Fmax,tot Kxmax,tot!!!). Il faut dterminer les forces maximales dans chaque modes partir de la matrice de rigidit, F 1n, max F 2n, max x 1n, max x 2n, max k 11 k 12 k 21 k 22 k 13 k 14 k 23 k 24 k 33 k 34 k 43 k 44 k 15 k 25 k 35 k 45 A 1n A 3n z n, max n=1 A 4n A 5n
F 3n, max = K x 3n, max = k 31 k 32 F 4n, max x 4n, max k 41 k 42 F 5n, max x 5n, max k 51 k 52
A 2n
k 53 k 54 k 55
rn F n, max = K x n, max = K A n z n, max = K A n ----------------------- S e, n ( n, ) (5.19) 2 m
n n
pour les superposer ensuite selon la rgle SRSS de lquation (5.18). F j, max, tot =
n=1
F jn, max
F51,max F41,max F31,max F21,max F11,max
F52,max F42,max F32,max F22,max F12,max
F53,max F43,max F33,max F23,max F13,max
F5,max,tot F4,max,tot F3,max,tot
...
F2,max,tot F1,max,tot
Fj1,max
Fj2,max
Fj3,max
( Fjn,max)
Fj,max,tot
EPFL-ENAC-SGC
Gnie parasismique
semestre dt 2005
P. Lestuzzi
15.
Systmes plusieurs degrs de libert
Matrice de flexibilit dune console: Calcul des dformations avec la matrice de flexibilit: d= Fa2(3L-a) 6EI
d = ^ F f si les hauteurs dtage sont constantes: ^ h3 2 fi,j = i (3j-i), j i
