ENSA Tétouan

Filière Génie Civil

Dynamique des structures

Cours 8:
Réponse sismique d'un bâtiment

Abdellatif Khamlichi

2020/2021 1
 Plan

• Notions de spectre de déplacement et de pseudo-spectre d'accélération

• Modélisation d'un bâtiment
• Exemple

Ce cours est adapté du cours d'Alain Pecker Ecole des Ponts Paris Tech
2
Notions de spectre de déplacement et de pseudo-spectre d'accélération

u(t)

F = ku + cu

a(t)
4
Accélération (m/s2)

a(t) 0

-2

-4
0 5 10 15 20 25 30
Temps (s)

Oscillateur à 1 ddl soumis à une excitation par le support 3

Notions de spectre de déplacement et de pseudo-spectre d'accélération

v(t) = u(t ) + a( t )
y(t) u(t)
f i = mv

+




= m(u a)

Fs = ku
Fd = cu
fi

O x

+


m(u a ) + cu
+ ku = 0

Oscillateur à 1 ddl soumis à une excitation par le support

4
Notions de spectre de déplacement et de pseudo-spectre d'accélération

m
u
+ cu
+ ku = −m

a

u
+ 2ξωu
+ ω2 u = −

a

Pulsation propre: k
ω=
m
Fréquence propre: 1 k
f=
2π m

Période propre: 1 m
T= = 2π
f k
c c
Coefficient d'amortissement: ξ= =
2 km 2mω

Chargement effectif: f eff (t) = − m

a( t)
5
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Comme les conditions initiales correspondent au repos, la réponse est donnée

par l'intégrale de Duhamel:

1 t
u(t) = −
ωd ∫0 a(τ)e −ξω(t −τ) sin [ ωd (t − τ)] dτ

2
où ωd = ω 1 − ξ
60

40
Force effective (kN)

20

−ma(t) 0

-20

-40

-60
0 5 10 15 20 25 30
Temps (s)
Chargement effectif pour m=15 t
6
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Séisme El Centro composante EW

4

Accélération (m/s2)
2

a(t) 0

-2

-4
0 5 10 15 20 25 30
Temps (s)

4πξ 4π 2 u(t) = u(ξ, T, t)

u+ u
+ 2 u = −


a
T T
Dans la suite, on prend:
592180
m = 15000kg , k =
T2 7
Notions de spectre de déplacement et de pseudo-spectre d'accélération

0.5

u(m)
T = 1s 0

-0.5
0 5 10 15 20 25 30
t (s)
0.5
u(m)

T = 1.5s 0

-0.5
0 5 10 15 20 25 30
t (s)
0.5
u(m)

T = 3s 0

-0.5
0 5 10 15 20 25 30
t (s)

Influence de la période de l'oscillateur sur la réponse; amortissement fixe ξ = 2%

sous le séisme El Centro composante EW
8
Notions de spectre de déplacement et de pseudo-spectre d'accélération

0.2

u (m)
ξ = 2% 0

-0.2
0 5 10 15 20 25 30
t (s)
0.2
u (m)

ξ = 5% 0

-0.2
0 5 10 15 20 25 30
t (s)
0.2
u (m)

0
ξ = 8%
-0.2
0 5 10 15 20 25 30
t (s)

Influence de la période de l'oscillateur sur la réponse; période fixe T = 1.5s

sous le séisme El Centro composante EW
9
Notions de spectre de déplacement et de pseudo-spectre d'accélération

0.2

0.15
u max =10.56cm
0.1

0.05
u (m)

-0.05

-0.1

-0.15

-0.2
0 5 10 15 20 25 30
t (s)

Réponse maximale en déplacement pour T=1.5s et ξ=5%

sous le séisme El Centro composante EW 10
Notions de spectre de déplacement et de pseudo-spectre d'accélération

0.5
5%
2%
0.4 8%
SD (T)

0.3
umax (m)

0.2

0.1

0
0 5 10 15 20 25 30
T (s)

Spectre de réponse déplacement sous le séisme

El Centro composante EW 11
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Effort élastique: fs (t) = ku(t)

k = mω ⇒ fs (t ) = m ω u ( t )
2 2

Force statique équivalente: fsmax = mSA

2
SA =ω SD
Spectre de pseudo-accélération

12
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Déplacement relatif: SD (T,ξ)

2π
Pseudo-vitesse: SV (T,ξ) =ωSD (T,ξ) = SD (T,ξ)
T
2 4π2
Pseudo-accélération: SA (T,ξ) =ω SD (T,ξ) = 2 SD (T,ξ)
T

Remarque:

Spectre de pseudo-vitesse est différent de spectre de vitesse

1 1 1

=−
u(t) ωu(t ) −


u(t) −


a(t)
2ξ 2ξω 2ξω
Spectre de pseudo-accélération est différent de spectre d'accélération
u(t) = −





−ω2 u(t)
a(t) − 2ξωu(t)
13
Notions de spectre de déplacement et de pseudo-spectre d'accélération

10
2%
5%
8
8%
amax (m/s2)
6

0
0 3 6 9 12 15
T (s)

Pseudo-spectre d'accélération sous le séisme El Centro composante EW

14
Notions de spectre de déplacement et de pseudo-spectre d'accélération

0.6

0.4

0.2
Accélération (m/s2)

-0.2

-0.4

-0.6

-0.8
0 5 10 15 20 25 30
Temps (s)

Un accélérogramme du séisme Kobe

15
Notions de spectre de déplacement et de pseudo-spectre d'accélération

40
2%
5%
amax (m/s2)
30 8%

20

10

0
0 2 4 6 8
T (s)

Pseudo-spectre d'accélération sous le séisme Kobe

16
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Spectre normalisé d'accélération; Eurocode 8

17
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Spectre normalisé de déplacement; Eurocode 8

18
Notions de spectre de déplacement et de pseudo-spectre d'accélération

Spectre normalisé RPS 2011 ; amortissement à 5%

Za numéro de la zone d'accélération entre 1 et 5

Zv numéro de la zone d'accélération entre 0 et 4 19
Modélisation d'un bâtiment

Table vibrante E-Defense Japon

20
Modélisation d'un bâtiment

Déplacements

Relatifs Absolus

 u1   v1  i
    mi
U =  ui  V =  vi 
   
   
u N  v N 

∆a s (t) u i (t)
V = U +∆a s (t) vi (t)

∆ Vecteur de composantes 0 ou 1.
Il donne la direction de l'accélération sismique

21
Modélisation d'un bâtiment

+ CU
MU
+ KU = −M∆


a s (t)

+ CU
MU
+ KU = P ( t )
eff

Chargement effectif: Peff (t) = −M∆

a s (t)

22
Modélisation d'un bâtiment

Méthode des éléments finis:

* Matrice de rigidité à structure bande; couplage entre nœuds voisins

uniquement.

* Matrice de masse:
- Matrice diagonale; masses concentrées aux nœuds;
- Matrice consistante (gain en précision pas nécessairement
significatif).

23
Modélisation d'un bâtiment

K −ωi2 M Φ i = 0
  N modes propres

Les sont définis à une constante multiplicative près (normalisation).

Φ Tj K −ωi2 M Φ i = 0

Orthogonalité des vecteurs modaux par rapport aux matrice de masse et de rigidité.

Φ Tj KΦ i
ωi2 =
Φ Tj MΦ i

24
Modélisation d'un bâtiment
Vibrations libres: solution générale

N
U = ∑ Ai Φ i sin (ωi t +ψ i )
i=1

Ai et ψ i sont déterminés par les conditions initiales

N
Déplacement initial: U(0) = ∑ Ai sin(ψi )Φ i
i=1

N
Vitesse initiale: V(0) = ∑ Ai cos(ψ i )ωi Φ i
i=1

Deux systèmes linéaires de N équations à N inconnues

25
Modélisation d'un bâtiment

Vibrations forcées amorties

+ CU
MU
+ KU = P ( t )
eff

N
Décomposition sur la base des modes propres: U = ∑qi (t) Φi
i=1

N N N
M ∑


q i (t) Φ i + C∑q
i (t) Φ i + K ∑q i (t) Φ i = P(t )
i=1 i=1 i=1

N
Φ Tj MΦ j


q j (t) + ∑ j ii
Φ T
C Φ

q (t) +Φ T
j K Φ q
j j (t) =Φ T
j P(t )
i=1

N
µ j
q
j (t ) + ∑Φ Tj C Φ i q
i (t ) + κ jq j (t ) = p j ( t )
i=1

Couplage des équations

26
Modélisation d'un bâtiment
La matrice d'amortissement

• L'amortissement est complexe

• Il ne peut être calculée

• L'énergie dissipée dépend:

- du frottement aux liaisons (acier boulonné);
- de l'ouverture/fermeture microfissures;
- de la déformation des éléments secondaires;
- du frottement structure/éléments secondaires;
- du comportement anélastique des matériaux;
-…

• Tous ces mécanismes sont fonction du niveau de déformation de la structure.

27
Modélisation d'un bâtiment
Définition de la matrice C

La matrice C peut être mesurée ou identifiée (mais essais complexes).

Souvent, on l'estime a priori via l'hypothèse de l'amortissement modal:

C est diagonale par rapport aux modes propres

Φ Tj C Φ i = c jδij

Cette hypothèse permet le découplage des équations

28
Modélisation d'un bâtiment
Mesure expérimentale de l'amortissement; bibliothèque Millikan (Caltech)

• Mesures des vibrations avec excitateur balayage en sinus

• Excitation par enregistrements sismiques:
- Lytle Creek (M=5.4, pga=0.05 g) 29
- San Fernando (M=6.4; pga=0.20 g)
Modélisation d'un bâtiment

Instrumentation de la bibliothèque Millikan; 36 canaux distribués

sur 9 étages 30
Modélisation d'un bâtiment

Mesure expérimentale de l'amortissement; bibliothèque Millikan

(Caltech) 31
Modélisation d'un bâtiment

32
Modélisation d'un bâtiment

Résultats de mesure expérimentale de l'amortissement; bibliothèque

Millikan (Caltech)

33
Modélisation d'un bâtiment
Niveau de déformation

Matériau <0.5 limite élastique Près limite élastique

3-5 7-10
Béton armé
2-3 7-10
Béton précontraint

Acier soudé 2-3 5-7

Acier boulonné 5-7 10-15

Bois cloué 5-7 15-20

Valeurs typiques d'amortissement en %

34
Modélisation d'un bâtiment

f d = cu

Force fd
u = u sin(ωt)
ceq  max
Ed
E c = E d
u
f d,max  max,c = u max,d = u max
Déplacement u

u max

T Ed
Ed = ∫ 2


= πceq ωu max
(ceq u)udt ⇒ ceq = 2
0 πωu max

Amortissement linéaire équivalent

35
Modélisation d'un bâtiment

Ed Ed k
fs,stat 2E s,stat
ceq = = k=
πωu 2max 2πE s,stat ω u 2max

Es,stat
ceq Ed
ξ= =
2mω 4πE s,stat

u max u

Taux d'amortissement linéaire équivalent

36
Modélisation d'un bâtiment

Système avec amortissement différent dans différentes zones

37
Modélisation d'un bâtiment

E si énergie élastique de déformation du sous système i

E di énergie dissipée par le sous système i

E di
ξi amortissement du sous système i ξi =
4πE si,stat

E d = ∑ E di = 4π∑ξi E si,stat ⇒ ξeq =

∑ ξi E si,stat
∑Esi,stat

Amortissement modal pondéré

38
Modélisation d'un bâtiment

Dans le cas d'un amortissement quelconque (non classique), en particulier si

le système présente un amortissement localisé, une méthode d'intégration pas
à pas est nécessaire.

En revanche, si l'amortissement est de type Rayleigh ou plus généralement

Caughey, on pourra pratiquer l'assemblage de matrices d'amortissement modal
(amortissement orthogonal).

Φ t CΦ = ci δij

Amortissement orthogonal: décomposition sur base des modes propres non

amortis.

Amortissement non orthogonal (quelconque): décomposition sur base des

modes propres amortis (modes complexes) ou intégration pas à pas.

39
Modélisation d'un bâtiment

Amortissement de Rayleigh

C = aM + bK

µ j


q j (t) + (aµ j + bκ j )q
j (t) + κ jq j (t) = p j (t )

aµ j + bκ j 1 a 
ξ j =  + bω j 
2ω jξ j = 2  ω j 

mj

Amortissement proportionnel

40
Modélisation d'un bâtiment

Généralisation de l'amortissement; amortissement de Caughey

N−1 
C = M ∑a  M K
−1

=0

N nombre de degrés de liberté du système

a M M −1K  0 = a M
 0   0 1 L−1
a M M −1K1 = a K ξ j = ∑a  ω2j −1
2 =0
 1   1
avec L ≤ N
a M M −1K  2 = a KM −1K
 2   2



41
Modélisation d'un bâtiment

Calcul de la matrice d'amortissement

Φ = Φ1 Φ 2  Φ j  Φ N 

Φ T CΦ = c jδij ⇔ Φ T CΦ =Ξ

c1 
 c2 
  0 
Ξ= 
 cj  c j = 2m jω jξ j
 0  
 c 
 N

42
Modélisation d'un bâtiment
On peut calculer a posteriori la matrice d'amortissement sous l'hypothèse
d'un amortissement modal (Caughey en général).

Φ T CΦ =Ξ ⇒ C = Φ − T ΞΦ −1
Φ −1 = µ  −1 Φ T M
  j
Φ T MΦ = µ j  ⇒  −1
−T
Φ = MΦ µ j  
−1 −1
C = MΦ µ j  Ξ µ j  Φ T M

 N 2ξ jω j 
C= M∑ T
Φ j ⊗Φ j  M
 j=1 m j 
 

D = A ⊗ B ⇔ Dij = Ai B j
43
Modélisation d'un bâtiment

Solution générale en vibrations forcées avec amortissement modal

µ j


q j (t) + 2ξ jω jµ jq
j (t) + κ jq j (t ) = p j (t)

avec

µ j =Φ Tj MΦ j

κ j =Φ Tj KΦ j

p j =Φ Tj P

1
q j (t) + 2ξ jω jq
j (t) +ω2j q j (t ) =


p j (t)
µj
44
Modélisation d'un bâtiment
Solution générale en vibrations forcées avec amortissement modal

q j (t) ⇒ Intégrale de Duhamel

1 t −ξ jω j (t −τ )
q j (t) =
µ jωdj ∫0
p j ( τ) e sin ωdj (t −τ) dτ

ωdj =ω j 1−ξ 2j

D'où N
U = ∑Φ jq j (t)
j=1

45
Modélisation d'un bâtiment

Solution générale en vibrations forcées en cas d'excitation du support

Relatifs Absolus

 u1   v1  i
    mi
U =  ui  V =  vi 
   
   
u N  v N 

∆a s (t) u i (t)
V = U +∆a s (t) vi (t)

∆ Vecteur de composantes 0 ou 1.
Il donne la direction de l'accélération sismique
46
Modélisation d'un bâtiment

Solution générale en vibrations forcées en cas d'excitation du support

1
q j (t) + 2ξ jω jq
j (t) +ω2j q j (t ) =


p j (t) p j (t) =Φ Tj P(t) = −Φ Tj M∆


a s (t)
µj

1 T Φ Tj M∆
Posons: a j = Φ j M∆ = T facteur de participation du mode j.
µj Φ j MΦ j Il dépend du mode de normalisation de Φj
alors
q j (t) + 2ξ jω jq
j (t) +ω2j q j (t) = −a j




a s (t )

D'où q j (t) = a jq j (t) et U j (t) = a jΦ jq j (t)

1 t −ξ ω (t −τ )
avec q j (t) = − ∫

a s (τ)e j j sin ωdj (t −τ) dτ
ωdj 0

47
Modélisation d'un bâtiment

Propriétés des facteurs de participation

N
∑ a jΦ j = ?
j=1

N
Φ iT M∆
ΦiT M ∑ a jΦ j =ΦiT MΦi a i =Φ iT MΦ i T
Φi MΦi
= Φ T
i M∆
j=1

N 
ΦiT M  ∑a jΦ j − ∆ = 0 pour tout i∈[1, N]
 j=1 

D'où
N
∑ a jΦ j = ∆
j=1

48
Modélisation d'un bâtiment
Calcul des efforts et des déplacements

Effort élastique:

Fj = KU j = KΦ jq j (t)

mais KΦ j =ω2j MΦ j Fj =ω2j MΦ jq j (t)

Fj = a jω2j MΦ jq j (t)

Solution globale (cumul des modes):

N N
F = ∑ Fj U = ∑U j
j=1 j=1
49
Modélisation d'un bâtiment
Cumul des modes

Fj (t), U j (t) Connus par leurs représentations temporelles

(intégrale de Duhamel)

Pas de difficulté: sommation à chaque t

Fj (t), U j (t) Connus par leurs valeurs maximales (spectre de réponse)

Règle de cumul

50
Modélisation d'un bâtiment

Valeur maximale par spectre de réponse

q j (t) + 2ξ jω jq
j (t) +ω2j q j (t) = −


a s (t)

 1 t −ξ jω j (t −τ ) 
SD (ω j ,ξ j ) = max q j (t) = max  − ∫


a (τ)e sin ωdj (t −τ) dτ 
t t  ω 0 s 
 dj 

51
Modélisation d'un bâtiment

Cumul des participations modales

Cumul en valeur absolue

N
R =∑ R j
j=1

Cumul SRSS (square root of the sums of the squares )

N
R= ∑R 2j
j=1

Formule non applicable lorsque les modes propres ont des fréquences voisines:

ωi
∈[0.9,1.1]
ωj

52
Modélisation d'un bâtiment

Cumul quadratique complet CQC

N N
R= ∑∑ρijR i R j
i=1 j=1

ρij Coefficient de corrélation entre les modes i et j

8 ξi ξ jωi ω j (ωi ξi +ω jξ j )ωi ω j

ρij =
(ωi2 −ω2j ) 2 + 4ξi ξ j (ωi2 +ω2j )ωi ω j + 4(ξi2 +ξ2j )ωi2 ω2j

53
Modélisation d'un bâtiment

Coefficient de corrélation
54
Modélisation d'un bâtiment

Choix du nombre de modes

• Système à N degrés de liberté possède N modes

• Les modes d'ordre supérieur contribuent "généralement" faiblement

à la réponse

Combien de modes retenir pour la réponse?

55
Modélisation d'un bâtiment
Masse modale
2

Posons: m*j =
( Φ Tj M∆ ) Masse modale dans la direction ∆
Φ Tj MΦ j

Elle ne dépend pas de la normalisation de Φj

2
2
 Φ j M∆ 
T
1 (Φ Tj M∆ ) 2 m*j
a j = T  = T =
 Φ j MΦ j  T
Φ j MΦ j Φ j MΦ j µ j
 
D'où

m*j
aj= et m*j =µ ja 2j
µj

56
Modélisation d'un bâtiment

Masse modale dans la direction ∆

M T = ∆ T M∆

Décomposition sur la base des modes: ∆ =ΦZ

M T = ZT Φ T MΦZ = ZT µ j  Z = ZT µ j Z j 

N
soit M T = ∑µ j Z2j
j=1

or Φ Tj M∆Φ Tj MΦZ Φ Tj MΦ j Z j
aj= T = T = T = Zj
Φ j MΦ j Φ j MΦ j Φ j MΦ j
D'où
N N
MT = ∑ µ ja 2j = ∑ m*j
j=1 j=1
57
Modélisation d'un bâtiment

Erreur commise

L'erreur commise sur les efforts élastiques lorsqu'une troncature modale

à n modes est opérée:
N
εF = ∑ Fj
j=n +1

avec

m*j
Fj = a jω2j MΦ jq j (t) = ω2j q j (t)MΦ j
µj
mj

On voit que l'erreur est fonction:

- des caractéristiques dynamiques de la structure {ω j ,Φ j}

- des caractéristiques du chargement q j (t)

- des masse modales m*j

58
Modélisation d'un bâtiment

L'ordre de troncature convenable n dépend de la variable d'intérêt:

χ = U,U,U,M,V,N, σ,...
En général:

χ = U ⇒ n << N petit

χ = V ⇒ n << N élevé
Le critère de convergence modale impose en général (c'est le cas du RPS2011
et de l'Eurocode 8) de satisfaire suivant toutes les directions sismiques
(deux au moins) que:
n
∑ j ≥ 0.9M T
m*

j=1

59
Modélisation d'un bâtiment
Prise en compte du déficit de masse modale (modes résiduels)
n N
déficit = M T − ∑ m*j = ∑ j
m*

j=1 j=n+1

Pour les fréquences > à la fréquence de coupure du spectre (30 à 50 Hz

en général), on admet que l'amplification dynamique est nulle. Cela veut dire
que les contributions modales de ces modes sont prises statiques:

q j (t) + 2ξ jω jq
j (t) +ω2j q j (t) = −


a s (t)

1
q j ( t ) = − 2


a s (t ) pour j>n
ωj

60
Modélisation d'un bâtiment

Spectre de réponse; fréquence de coupure

(dans le RPS 2011 et l'Eurocode 8 c'est 33Hz)
61
Modélisation d'un bâtiment

Prise en compte des modes résiduels

n N n N aj
U = ∑a jΦ jq j (t) + ∑ a jΦ jq j (t) U = ∑a jΦ jq j (t) − ∑ Φ j


a s (t)
2
j=1 j=n +1 j=1 j=n+1 ω j

n aj N aj
Posons: W = ∑ ω2 Φ j + ∑ ω2 Φ j
j=1 j j=n +1 j

Cette quantité s'interprète comme étant la déformé statique dans un champ

d'accélération unitaire


a s (t) = −1
n aj N aj
W −∑ 2
Φj = ∑ ω2 Φ j
j=1 ω j j=n +1 j

n  n a 
U = ∑a jΦ jq j (t) −


a s (t)  W − ∑ 2 Φ j 
j
 ω 
j=1  j=1 j 
Pseudo-mode 62
Exemple

1 u1
600
1.5 u2
1200
2 u3
1800

 1 −1 0  1 0 0
K = 600 −1 3 −2 M = 0 1.5 0
 0 −2 5  0 0 2
   

Modèle discret à 3 ddls

63
Exemple

Modèle discret: Modèle EF poutre plane:

v1
θ1
u1 1
u1 0
v2 0
θ2
u2  
u2 1
v3 ∆ = 0
θ3 0
u3 1
u3
 
0
0

1 avec les matrices de masse

∆ = 1 et de rigidité assemblées
1 à partir des matrices élémentaires
 du modèle Euler-Bernoulli
64
Exemple

600 −ω2 −600 0 

K −ω2 M =  −600 1800 −1.5ω2 −1200 

  
 0 −1200 3000 − 2ω2 

----> Polynôme du troisième degré en ω2

ω1 =14.52rd/s ω2 = 31.05rd/s ω3 = 46.1rd/s

600 −ωi2 −600 0   ϕ  0

K −ωi2 M Φ i =  −600 1800 −1.5ωi2   i1   
−1200  ϕi2 = 0
      
 0 −1200 3000 − 2ωi2  ϕi3  0

On normalise les modes en imposant: ϕi1 =1

65
Exemple

ω1 =14.52rd/s ω2 = 31.05rd/s ω3 = 46.1rd/s

T1 = 0.43s T2 = 0.2s T3 = 0.14s
µ1 =1.8 µ 2 = 2.46 µ3 = 23.11
κ1 = 380.05 κ 2 = 2367.2 κ3 = 49110

Modes propres 66
Exemple

Caractéristiques modales

1 u1 Mode 1 Mode 2 Mode 3

1.5 u2 ωj 14.5 31.1 46.1

2 u3 aj 1.424 -0.51 0.09

3.669 0.647 0.186

m*j 81.5% 14.4% 4.1%

67
Exemple

Vibrations libres: mouvement du nœud 2 pour D T0 = [0.5 0.5 0.5] 68

Exemple

Détermination de l'amortissement proportionnel 69

Exemple

Détermination de l'amortissement proportionnel Caughey

70
Exemple

Enregistrement lake Hughes EW

71
Exemple

1 u1

1.5 u2

2 u3

Déplacement relatif Nœud 2

72
Exemple

Force au nœud 2
73
Exemple

Forces 74
Modélisation d'un bâtiment

Valeur maximale par spectre de réponse

75
Exemple

T1

Déplacement 76
Exemple

T3 T2 T1

Pseudo-accélération 77
Exemple

Spectres de réponse

Déplacement relatif: SD (ω,ξ)

Pseudo-accélération: SA (ω,ξ) =ω2SD (ω,ξ)

U jmax = a jΦ jSD (ω j ,ξ j )

Fjmax = a jMΦ jSA (ω j ,ξ j )

78
Exemple

F1 = a1MΦ1SA (ω1 ,ξ1 ) F2 = a 2 MΦ 2SA (ω2 ,ξ 2 ) F3 = a 3MΦ 3SA (ω3 ,ξ3 )

Forces modales 79
Exemple

∑ Fi SRSS CQC

1 mode 7.43 7.43 7.43

2 modes 11.5 8.477 7.99

12.59 8.54 8.63

3 modes

Cumul de l'effort tranchant en pied; la valeur exacte est 7.99

80