Université Privée de Marrakech

CHAPITRE 4 : APPLICATION ET
ETUDE DE CAS

Z. KAHOT Année Universitaire: 2023/2024

 CHAPITRE 4 : APPLICATION ET ETUDE DE CAS

I. Description du projet :

Construction simple avec un RDC et un étage

RDC H=3m Q = 250 kg/m² G = 800 kg/m²

PREMIER ETAGE H=3m Q = 250 kg/m² G = 800 kg/m²

• Type d'acier utilisé : Barres haute adhérence Fe500

• La résistance à la compression à 28j égale à 25 MPa.

• Poteau carré avec a = 30cm.

Les données sismiques :

Coefficient d'accélération de la zone 3 A = 0.16

Coefficient de site 2 S = 1.2
Facteur d'amplification dynamique D = 2.5
Coefficient de priorité I=1
Facteur de comportement K=2
Coefficient 𝚿 Ψ = 0.2
II. Détermination des pulsations et des modes propres :

Pour ce qui concerne l’analyse modale, la détermination des modes propres doit être faite dans le cas
d’un système libre (F=0) non-amorti (C=0). Cela revient à résoudre le système suivant :

𝑴𝑿̈ + 𝑲𝑿 = 𝟎 (1)

On cherche une solution sous la forme :

𝑿(𝒕) = 𝝋𝒆𝒊𝝎𝒕

Remplaçons cette solution dans (1) et on arrive au système suivant :

(𝑲 − 𝝎𝟐𝑴)𝝋 = 𝟎

𝝎² = 𝝀 = 𝒗𝒂𝒍𝒆𝒖𝒓 𝒑𝒓𝒐𝒑𝒓𝒆

L’équation caractéristique de ce système :

|𝑲 − 𝝀𝑴| = 𝟎

𝒌𝟏 + 𝒌𝟐 −𝒌𝟐 𝒎𝟏 𝟎
|( ) − 𝝀( 𝟎
−𝒌𝟐 𝒌𝟐 𝒎𝟐)| = 𝟎

𝟏𝟐𝐄𝐈
𝑲𝟏 = 𝑲2 = 𝑲 = 𝟒 𝐱 = 𝟑8. 𝟔 𝑴𝑵/𝒎
𝐇𝟑

Avec :
𝟑
𝑬 = 𝟏𝟏𝟎𝟎0 √𝐟𝐜 𝟐𝟖 = 𝟑𝟐𝟏𝟔𝟒 𝑴𝑷𝒂

𝐚𝟒
𝑰= = 6.75 x 𝟏𝟎−𝟒 𝐦𝟒
𝟏𝟐

Et 𝒎𝟏 et 𝒎𝟐 les masses des deux planchers 1 et 2 (on néglige le poids des poutres devant la dalle) :

𝒎𝟏 = 𝒎2 = 𝒎 = 𝑮 ∗ 𝑺 + 𝚿 ∗ 𝑸 ∗ 𝑺 = 𝟖𝟎𝟎 ∗ 𝟐𝟎 + 𝟎. 𝟐 ∗ 𝟐𝟓𝟎 ∗ 𝟐𝟎 = 𝟏𝟕𝟎𝟎𝟎𝑲𝒈

|(𝟐𝒌 −𝒌) − 𝝀 (𝒎 𝟎 )| = 𝟎
−𝒌 𝒌 𝟎 𝒎

𝒎𝟐𝝀𝟐 − 𝟑𝒌𝒎𝝀 + 𝒌𝟐 = 𝟎

D’où :

{ 𝝀𝟏 = 𝟖𝟔𝟕. 𝟐𝟖 ⇔ {𝝎𝟏 = 𝟐𝟗. 𝟒𝟒𝒓𝒂𝒅/𝒔 ⇔ { 𝒇𝟏 = 𝟒. 𝟔𝟖𝑯𝒛 𝑻𝟏 = 𝟎. 𝟐𝟏𝒔

𝝀𝟐 = 𝟓𝟗𝟒𝟒. 𝟒𝟕 , 𝝎𝟐 = 𝟕𝟕. 𝟏𝒓𝒂𝒅/𝒔 , 𝒇𝟐 = 𝟏𝟐. 𝟐𝟕𝑯𝒛 ⇔ {𝑻𝟐 = 𝟎. 𝟎𝟖𝒔

D’après ce système on détermine les modes propres :

(𝑲 − 𝝀𝒊𝑴)𝝓𝒊 = 𝟎

 Mode 1 :
6
−38.6 106 ) − (λ1 . 17 0 𝜑11 0
⌊( 77.2 10 6 ) 103 ⌋ (𝜑 ) = ( )
−38.6 10 38.6 106 0 λ1 . 17 21 0

𝜑
Φ1 = (𝜑11 ) = (0.618)
21 1

 Mode 2 :

6
−38.6 106 ) − (λ2 . 17 0 𝜑21 0
⌊( 77.2 10 6 ) 103 ⌋ (𝜑 ) = ( )
−38.6 10 38.6 106 0 λ2 . 17 22 0

𝜑
Φ2 = (𝜑21 ) = (−1.618)
22 1

Représentation graphique des modes :

Mode 1 Mode 2

4
III. Normalisation des masses et des modes propres :

1- Normalisation de la matrice de masse M :

2- Normalisation des modes propres :

5
IV. Réponse de vibration libre non-amorti :

Considérons l’équation de mouvement de notre système à deux degré de liberté avec les conditions
initiales suivantes :

(2)

On fait un changement de variable de coordonnée :

Soit {𝝓i} une base orthogonale

X (𝒕) = 𝝓 𝒒 (𝒕) = ∑𝟐𝒊=𝟏 𝝋𝒊𝒒𝒊(𝒕) (3)

𝒒 (𝒕) : Vecteur des coordonnées modales ou cordonnées généralisées.

𝝓 : La matrice modale.

Remplaçons (3) dans (2) et on obtient :

𝝓𝑻 𝑴𝝓 𝒒̈ + 𝝓𝑻 𝑲𝝓 𝒒 = 𝟎

𝑴 ̃𝒒̈ +𝑲 ̃ 𝒒 = 𝟎

La solution de cette équation différentielle est :

𝒒𝒊(𝒕) = 𝑨𝒊 𝐜𝐨𝐬(𝝎𝒊𝒕) + 𝑩𝒊𝐬𝐢𝐧 (𝝎𝒊𝒕)

Avec :

6
Après les calculs on arrive aux résultats suivants :

Par la superposition modale on trouve le vecteur déplacement suivant en fonction du temps :

V. Réponse sismique :

Le calcul de la réponse du bâtiment sous l’effet de séisme peut se faire par deux méthodes :

 L’analyse modale temporelle peut être effectuée à partir des accélérogrammes adaptés au site de
la construction, le calcul est en fonction du temps.

 L’analyse modale spectrale basé sur des courbes de réponse des constructions dans chaque site,
il nous donne les réponses maximums de la construction (vitesse, accélération ou déplacement).

1- Analyse modale temporelle :

L’équation du mouvement de notre bâtiment dans le cas d’une excitation extérieur F (force de séisme)
est :

(4)
Avec : F(𝒕) = −M𝑿̈𝒈 (𝒕)

7
𝑿̈ 𝒈 (𝒕) : est l’accélération du sol.

La solution de (4) est donnée par l’intégral de Duhamel :

Avec :
𝝎𝑎 = 𝝎𝟎√1 − ξ : Est la pulsation propre du système amorti (pseudo pulsation).

𝝃 = 5% : Est l’amortissement du bâtiment.

Par le changement de coordonnée :

𝑿(𝒕) = 𝝓𝒒(𝒕) = ∑𝟐𝒊=𝟏 𝝋𝒊𝒒𝒊(𝒕)

On arrive à l’expression suivante :

Avec 𝑳 = 𝝓𝑻 𝑴 𝒓 est le vecteur de participaient modale.

Sur certaines conditions sur l’expression de la matrice d’amortissement C, les propriétés d’orthogonalités
sont applicables à C (𝜙𝑖 𝑇 𝐶 𝜙𝑗 = 0 𝑝𝑜𝑢𝑟 𝑖 ≠ j). On va supposer la proportionnalité de C par rapport à M
et K.

𝑪 = 𝜶𝑴 + 𝜷𝑲

Et on arrive à :

𝜶 et 𝜷 : sont deux constante à déterminer lors de l’expérience.

1-1 Calcul des déplacements :

Pour le mode i :

(5)

8
La solution de (5) est donnée aussi par l’intégral de Duhamel :

Par la superposition modale :

 Mode 1 :

𝒎̃𝟏 = 𝟏 on travaille dans la base orthogonale

 Mode 2 :

9
1-2 Calcul des forces dynamiques:

Ce sont des forces élastiques correspond au déplacement relatif :

𝒇 = 𝑲𝑿 = 𝝓𝑲𝒒(𝒕)

Pour le mode i :

𝒇𝒊 = 𝑴𝝓𝒊 𝝎𝒊²𝒒𝒊 (𝒕) = 𝑴𝝓𝒊 𝝎𝒊² 𝐦̃𝐋𝐢𝐢 𝑫𝒊 (𝒕) = 𝑴 𝝓𝒊 𝐦̃𝐋𝐢𝐢 𝑨𝒊 (𝒕) Avec : 𝑨𝒊(𝒕) = 𝝎𝒊² D𝒊(𝒕) c’est la pseudo-
accélération

 Mode 1 :

 Mode 2 :

D’où l’effort tranchant V à la base est :

10
 2- Analyse spectrale :

Avec :
• A: coefficient d'accélération de la zone.
• 𝜼: facteur de correction d'amortissement 𝜼 = √3.5 + ξ ≥ 𝟎.𝟕
• 𝝃: Coefficient d'amortissement critique.
• R: coefficient de comportement de la structure.
• T1et T2 : périodes caractéristiques à la catégorie du site.
• Q: facteur de qualité.

Les valeurs de T1 et T2 :

Site S1 S2 S3 S5
T1(s) 0.15 0.15 0.15 0.15
T2(s) 0.30 0.40 0.50 0.70

Et tel que notre structure situe au site 2 donc : T1 = 0.15s et T2 = 0.40s

Facteur de qualité Q :

La valeur de Q est déterminée par la formule suivante :

𝟓
𝑸 = 𝟏 + ∑ 𝒑𝒒
𝟏
Pq est la pénalité à retenir selon que le critère de qualité q est satisfait ou non sa valeur est donnée par le
tableau suivant :

11
 Pq
Critère q Observé N/observé
Condition minimales sur les files de 0 0.05
contreventement
Redondance en plan 0 0.05
Régularité en plan 0 0.05
Régularité en élévation 0 0.05
Contrôle de qualité des matériaux 0 0.05
Contrôle 0 0.10

Et pour notre cas tous ces critères sont observés donc Pq = 0 donc le facteur de qualité égale à 1.
Q=1

Coefficient de comportement de la structure R :

Sa valeur est unique en fonction du système de contreventement. Notre exemple est structure en béton
armé avec un système portique sans remplissage en maçonnerie la valeur de R est :

R=5

• Spectre d'accélération du mode 1 :

0.15≤ 𝑇 = 0.21𝑠 ≤ 0.4

Avec :

S𝒂𝟏 = 𝟏. 𝟖𝟖 𝒎/𝒔²

• Spectre d'accélération du mode 2 :

0 ≤ 𝑇 = 0.08𝑠 ≤ 0.15

12
Avec :

S𝒂𝟐 = 𝟏. 𝟗𝟑 𝒎/𝒔²

2-1 Calcul des déplacements :

 Mode 1 :

 Mode 2 :

Et par l'application du CQC :

13
2-2 Calcul des forces dynamiques :

Pour le mode i :

 Mode 1 :

𝒎̃𝟏 = 𝟏 on travaille dans la base orthogonale

 Mode 2 :

En appliquant la combinaison quadratique complète pour trouver la force de cisaillement fdyn à la base :

VI. Méthode statique équivalente :

Cette méthode est basée sur le premier mode seulement, elle est valable pour les bâtiments dont la
réponse est entièrement déterminer de mode fondamentale.

14
 Les forces sismiques horizontales agissant sur les masses de la structure sont représentées par la force
équivalente de cisaillement à la base agissant dans la direction du calcul qui est donnée par la formule
suivante :

La charge W de la structure correspond à la totalité des charges permanentes G et une fraction

q des charges d’exploitation Q en fonction de la nature des charges et leur durée. On prend :

𝑊 = 𝐺 + Ψ𝑄 = 340𝐾𝑁

𝑽 = 𝟖𝟏. 𝟔𝑲𝑵

VII. Interprétation des résultats :

On constate que le mode 1 qui nous a donné les valeurs totales des déplacements, ça c’est grâce à ce que
la construction est régulière, donc on peut appliquer la méthode statique équivalente pour la
détermination des déplacements qui est basée sur le premier mode seulement.

Pour les valeurs des forces dynamiques fdyn finales n’égalent pas aux valeurs des forces trouvées par le
mode 1, à cause des modes supérieurs qui contribuent à la réponse de la force donc on a besoin d’aller
plus que de deux modes pour trouver les forces.

La valeur de l’effort latéral sismique servant au calcul ne doit pas être inférieure à 0.90 fois la valeur
obtenue par l’approche statique équivalente :

𝒇𝒅𝒚𝒏 ≤ 𝟎. 𝟗𝑽 ⇔ 𝟔𝟐. 𝟓𝟗 ≤ 𝟕𝟑. 𝟒𝟒  Vérifiée

L’origine de ce pourcentage de 90% provient essentiellement du pourcentage de l’erreur e des masses

immobiles, c.à.d. que lors de mouvement du bâtiment le pourcentage des masses immobiles doit
être ≤ 10%.

Ce qui applique, que les masses qui vont être en mouvement (les masses dynamiques 𝑀𝑘∗ ) par rapport à
la masse totale Mtot de la construction, doivent présenter un pourcentage ≥ 90%.

15
Revenons maintenant à notre cas, la masse modale qui représente la masse dynamique est donnée par
cette formule :

 Mode 1 :

 Mode 2 :

Donc le pourcentage cumulé des deux modes est :

16