Université Moulay Ismail

École Nationale Supérieure d’Arts et Métiers

Travaux Pratiques sur le Calcul

Parasismique des structures
Sous le thème :

CALCUL SISMIQUE D’UN BÂTIMENT

Réalisé par
HAFIDI Oussama
AIT SI YOUSSEF Fatimezzahra
RABENANDRASANA Geraldo

Encadré par
M. BAHI
M. MHARZI

Établissement : ENSAM Meknès, Marjane 2, B.P. 15290 El Mansour Meknès 50500

 Table des matières
1 Présentation générale 3

2 Calcul des efforts sismiques selon le RPS2011 4

2.1 Type de bâtiment et classification . . . . . . . . . . . . . . . . . . . . . . 4
2.2 Localisation sismique . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.3 Capacité de dissipation d’énergie (facteur k) . . . . . . . . . . . . . . . . 4
2.4 Conditions du sol . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
2.5 Calcul de la période propre de vibration . . . . . . . . . . . . . . . . . . 4
2.6 Détermination du facteur dynamique D . . . . . . . . . . . . . . . . . . . 5
2.7 Charge totale appliquée à la structure . . . . . . . . . . . . . . . . . . . 5
2.8 Calcul de la force sismique horizontale équivalente . . . . . . . . . . . . 5
2.9 Distribution des efforts sismiques selon les niveaux . . . . . . . . . . . . 6

3 Évaluation des déplacements inter-étages 6

3.1 Détermination des moments d’inertie et des rigidités latérales . . . . . . 6
3.2 Rigidité globale par niveau . . . . . . . . . . . . . . . . . . . . . . . . . . 8
3.3 Calcul des déplacements entre étages . . . . . . . . . . . . . . . . . . . 8

4 Position des centres de masse et de rigidité 9

4.1 Localisation du centre de masse . . . . . . . . . . . . . . . . . . . . . . 9
4.2 Centre de rigidité de l’étage . . . . . . . . . . . . . . . . . . . . . . . . . 10
4.3 Position du centre de rigidité – Étages bas (1er et 2e ) . . . . . . . . . . . 10
4.4 Position du centre de rigidité – Étages hauts (3e et 4e ) . . . . . . . . . . 11

5 Moments de torsion induits par les excentricités 12

6 Vérification de la stabilité globale de la structure 13

7 Simulation numérique avec Robot Structural Analysis 14

8 Analyse comparative des résultats 14

1
 Table des figures
1 Impact des séismes sur les bâtiments . . . . . . . . . . . . . . . . . . . 3
2 Surfaces du batiment . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
3 Illustration de l’excentricité horizontale . . . . . . . . . . . . . . . . . . . 12

Liste des tableaux

1 Données sismiques pour Al Hoceima . . . . . . . . . . . . . . . . . . . . 4
2 Charges permanentes et variables par étage . . . . . . . . . . . . . . . 5
3 Répartition des charges et efforts sismiques par niveau . . . . . . . . . 6
4 Rigidité latérale des poteaux par étage . . . . . . . . . . . . . . . . . . . 7
5 Rigidités latérales des voiles selon x . . . . . . . . . . . . . . . . . . . . 7
6 Rigidités latérales des voiles selon y . . . . . . . . . . . . . . . . . . . . 7
7 Rigidité totale en x et y par niveau . . . . . . . . . . . . . . . . . . . . . 8
8 Déplacements inter-étages par direction et par niveau . . . . . . . . . . 8
9 Calcul du centre de masse à partir des surfaces pondérées . . . . . . . 9
10 Calcul du centre de rigidité Xr . . . . . . . . . . . . . . . . . . . . . . . 10
11 Calcul du centre de rigidité Yr . . . . . . . . . . . . . . . . . . . . . . . . 11
12 Calcul du centre de rigidité Xr basé sur les raideurs selon y . . . . . . . 11
13 Calcul du centre de rigidité Yr à partir des raideurs selon x . . . . . . . . 12
14 Moments de torsion pour un séisme dans la direction Y . . . . . . . . . 13
15 Moments de torsion pour un séisme dans la direction X . . . . . . . . . 13
16 Calcul des indices de stabilité θsx et θsy par étage . . . . . . . . . . . . . 13

2
Travaux pratiques sur le calcul parasismique des structures

1 Présentation générale
Dans le cadre de ce travail pratique, nous analysons le comportement d’un bâ-
timent situé à Al Hoceima, une région connue pour son activité sismique. L’édifice
étudié comprend quatre niveaux, avec une utilisation différenciée : les deux premiers
étages sont destinés à la restauration, tandis que les deux derniers sont réservés à
l’habitation.
L’objectif principal est de vérifier la stabilité de la structure face aux actions sis-
miques, conformément aux exigences du Règlement de construction parasismique
RPS 2000. Pour ce faire, nous appliquons la méthode statique équivalente, en pre-
nant comme base un sol rocheux de bonne qualité. L’étude prend en considération les
paramètres géométriques du bâtiment ainsi que les charges permanentes et d’exploi-
tation appliquées aux différents niveaux.

F IGURE 1 – Impact des séismes sur les bâtiments

Rapport de projet 3 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

2 Calcul des efforts sismiques selon le RPS2011

2.1 Type de bâtiment et classification

Le bâtiment étudié sert à plusieurs usages (habitation et commerce), ce qui le
classe dans la **catégorie 3** du RPS2011.
— Coefficient d’importance : I = 1

2.2 Localisation sismique

Le projet se situe dans la ville de Al Hoceima, zone sismique importante. Les
données suivantes sont extraites du zonage sismique pour une probabilité de dépas-
sement de 10% en 50 ans :

Paramètre Valeur
Vitesse sismique v 17 cm/s
Zone en vitesse Zv 4
Accélération sismique g 18%
Zone en accélération Za 4

TABLE 1 – Données sismiques pour Al Hoceima

2.3 Capacité de dissipation d’énergie (facteur k)

Ce facteur tient compte de la capacité de la structure à absorber l’énergie d’un
séisme grâce à sa déformation. Pour notre cas :
— Classe 3 avec 0,10 < v ≤ 0,2 : système mixte voile-portique
— Coefficient de ductilité : k = 2

2.4 Conditions du sol

Le sol est un terrain rocheux dur d’épaisseur inférieure à 30 mètres. Cela permet
d’utiliser :
— Coefficient de site S = 1

2.5 Calcul de la période propre de vibration

La période fondamentale T est estimée par la formule suivante, adaptée aux struc-
tures non contreventées :
H
T = 0,09 √
L
Avec :
— Hauteur totale H = 14 m

Rapport de projet 4 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

— Longueurs dans les deux directions : Lx = 15 m, Ly = 15 m

Tx = 0,325 s Ty = 0,349 s

2.6 Détermination du facteur dynamique D

Le facteur dynamique est lié à la période propre. Comme Tx et Ty sont compris
entre 0,25 et 0,5 :
Za
=1
Zv

Dx = −2,4Tx + 3,1 = 2,32 Dy = −2,4Ty + 3,1 = 2,26

2.7 Charge totale appliquée à la structure

La charge sismique W comprend les charges permanentes G et une part des
charges d’exploitation Q pondérées selon l’usage. Le bâtiment a 4 niveaux :

Étage Usage G (kN/m²) ψ

1 Restauration 5,27 0,4
2 Restauration 5,27 0,4
3 Habitation 5,27 0,2
4 Habitation 6,29 0,2

TABLE 2 – Charges permanentes et variables par étage

La charge totale est donc :

Wtot = G1 + ψQ1 + G2 + ψQ2 + G3 + ψQ3 + G4 + ψQ4 = 24,57 kN/m2

2.8 Calcul de la force sismique horizontale équivalente

La force sismique latérale est donnée par la formule :
v · S · D · I ·Wtot
F=
K
On l’applique dans les deux directions :

Fx = 17·1·2,32·1·24,57

2 = 901,2 kN
17·1·2,26·1·24,57
Fy = 2 = 877,9 kN

Rapport de projet 5 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

2.9 Distribution des efforts sismiques selon les niveaux

Pour répartir la force sismique totale entre les différents étages du bâtiment, on
utilise la formule suivante :
 
Wn × hn
Fn = (F − Ft ) ×
∑(Wi × hi )
Avec :
— Fn : force horizontale appliquée à l’étage n
— Wn : charge gravitaire à l’étage n
— hn : hauteur de l’étage n par rapport au sol
— Ft : effort appliqué en tête de structure. Ici, comme Tx et Ty sont inférieurs à 0,7
s, on a Ft = 0
Ainsi, la formule devient :
  
 Fnx = Fx · Wn ·hn
 ∑(Wi ·hi ) 
 Fny = Fy · Wn ·hn
∑(Wi ·hi )

Le tableau ci-dessous résume les résultats obtenus pour chaque étage, dans les
deux directions.

Niveau G (kN) Q (kN) totale wi (kN) Hauteur hi (m) wi · hi (kN) Effort Fx (kN) Effort Fy (kN)
0 – – – – – 901,21 877,90
1 5,27 2,50 1166,22 4 4664,88 99,59 97,02
2 5,27 2,50 1166,22 8 9329,76 199,18 194,03
3 5,27 1,50 1036,02 11 11396,22 243,30 237,01
4 6,26 1,00 1201,56 14 16821,84 359,13 349,85

TABLE 3 – Répartition des charges et efforts sismiques par niveau

3 Évaluation des déplacements inter-étages

3.1 Détermination des moments d’inertie et des rigidités latérales

1. Étude des poteaux

Les poteaux ont une section rectangulaire de 30 × 60 cm2 , soit :

— Largeur b = 30 cm, hauteur h = 60 cm
Le calcul des moments d’inertie dépend de la disposition des poteaux :
— Pour les 4 poteaux du rez-de-chaussée :

h · b3 b · h3
Ix = Iy =
12 12
Rapport de projet 6 ENSAM-MEKNES
Travaux pratiques sur le calcul parasismique des structures

— Pour les autres poteaux, la section est retournée :

b · h3 h · b3
Ix = Iy =
12 12

Les rigidités horizontales Kx et Ky sont données par :

12 · E · Iy 12 · E · Ix
Kx = Ky =
Li3 Li3
Le tableau ci-dessous résume les résultats pour les 19 poteaux, étage par étage :

Étage Ei (GPa) Li (m) Ix,i (mm4 ) Iy,i (mm4 ) Kx,i (N/mm) Ky,i (N/mm)
1 45 4 8,64e10 4,185e10 353109,38 729000
2 45 4 8,64e10 4,185e10 353109,38 729000
3 45 3 8,64e10 4,185e10 837000 1728000
4 45 3 8,64e10 4,185e10 837000 1728000

TABLE 4 – Rigidité latérale des poteaux par étage

2. Étude des voiles

Deux voiles ont été modélisés avec une épaisseur commune de 30 cm. Les hau-
teurs considérées sont :
— Voile 1 : hv1 = 270 cm
— Voile 2 : hv2 = 440 cm
Les formules utilisées sont les mêmes que pour les poteaux :

e · h3 e · b3
Ix = Iy =
12 12
Les coefficients de rigidité pour chaque voile sont ensuite calculés.

Étage Ei (GPa) Li (m) Iv1,x (mm4 ) Iv2,x (mm4 ) Kx,v1 (N/mm) Kx,v2 (N/mm)
1 45 4 4,92075 · 1011 2,1296 · 1012 45562,5 74250
2 45 4 4,92075 · 1011 2,1296 · 1012 45562,5 74250
3 45 3 4,92075 · 1011 2,1296 · 1012 108000 176000
4 45 3 4,92075 · 1011 2,1296 · 1012 108000 176000

TABLE 5 – Rigidités latérales des voiles selon x

Étage Ei (GPa) Li (m) Iv1,y (mm4 ) Iv2,y (mm4 ) Ky,v1 (N/mm) Ky,v2 (N/mm)
1 45 4 6,075 · 109 9,9 · 109 3690562,5 15972000
2 45 4 6,075 · 109 9,9 · 109 3690562,5 15972000
3 45 3 6,075 · 109 9,9 · 109 8748000 37859555,56
4 45 3 6,075 · 109 9,9 · 109 8748000 37859555,56

TABLE 6 – Rigidités latérales des voiles selon y

Rapport de projet 7 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

3.2 Rigidité globale par niveau

Pour chaque niveau, on additionne les contributions des poteaux et des voiles afin
d’obtenir les rigidités totales dans les directions x et y :
(
Kx_tot = Kx,poteaux + Kx,v1 + Kx,v2
Ky_tot = Ky,poteaux + Ky,v1 + Ky,v2
Le résultat est résumé dans le tableau suivant :

Étage Kx,total (N/m) Ky,total (N/m)

1 487898,44 22849383
2 487898,44 22849383
3 1156500 54161500
4 1156500 54161500

TABLE 7 – Rigidité totale en x et y par niveau

3.3 Calcul des déplacements entre étages

Pour estimer le déplacement relatif entre deux niveaux successifs, on se base sur
le lien classique entre force et déplacement :

F = K ·∆
Dans cette formule :
— F est la somme des efforts appliqués à l’étage considéré,
— K représente la rigidité globale du niveau,
— ∆ correspond au déplacement inter-étage, noté Vi,i+1 .
Les formules deviennent alors :
F

Vi,i+1_x = x cumulée

 Kx total
Fy cumulée
Vi,i+1_y =


Ky total

Les résultats calculés pour chaque niveau sont présentés ci-dessous :

Déplacement inter-étage selon x (mm) Déplacement inter-étage selon y (mm)

1,8471 0,0384
1,6430 0,0342
0,5209 0,0108
0,3105 0,0065

TABLE 8 – Déplacements inter-étages par direction et par niveau

Rapport de projet 8 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

4 Position des centres de masse et de rigidité

4.1 Localisation du centre de masse

Le centre de masse d’un étage correspond au barycentre pondéré des surfaces
qui le composent. Lorsque la masse est répartie uniformément, on peut remplacer la
masse par les surfaces dans le calcul :

∑Sj · Xj

XCM =


∑Sj
∑ S j ·Y j
 YCM =


∑Sj
Ici, S j représente la surface de la zone j, et (X j ,Y j ) les coordonnées de son centre.
L’image suivante montre les surfaces prises en compte dans le calcul :

F IGURE 2 – Surfaces du batiment

Surface S (m2 ) X j (m) Y j (m) S · Xj S ·Y j XCMi YCMi

A 12 2 1,5 24 18
B 12 6 1,5 72 18
C 12 10 1,5 120 18
D 20 2 5,5 40 110
E 20 6 5,5 120 110
F 20 10 5,5 200 110
G 15 13,5 5,5 202,5 82,5
H 20 2 10,5 40 210
I 20 6 10,5 120 210
J 20 10 10,5 200 210
K 15 13,5 10,5 202,5 157,5 7,21 6,74
Total 186 — — 1341 1254

TABLE 9 – Calcul du centre de masse à partir des surfaces pondérées

Rapport de projet 9 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

4.2 Centre de rigidité de l’étage

Le centre de rigidité représente le point d’un étage où l’application d’une force ho-
rizontale engendre une translation sans torsion. Il est obtenu à partir de la distribution
des rigidités des éléments structuraux (poteaux, voiles, etc.).
Les coordonnées de ce point sont calculées avec :

∑ Ky, j · X j

Xcr =


∑ Ky, j
∑ Kx, j ·Y j
 Ycr =


∑ Kx, j
À noter :
— Les calculs sont faits séparément pour les étages 1 et 2, et pour les étages 3 et
4, car leurs hauteurs diffèrent, influençant les rigidités.

4.3 Position du centre de rigidité – Étages bas (1er et 2e )

Pour déterminer le centre de rigidité de chaque étage, on applique la méthode des
moments pondérés selon les rigidités des éléments porteurs (poteaux et voiles). Le
principe est de calculer une moyenne pondérée des coordonnées de chaque élément.
Le tableau ci-dessous regroupe les données nécessaires au calcul du centre de
rigidité pour les deux premiers niveaux :

Poteaux et voiles X j (m) Ky j (N/mm) Ky j · X j

P1 0 11390,625 0
P2 4 11390,625 45562,5
P3 8 11390,625 91125
P4 12 11390,625 136687,5
P5 0 45600 0
P6 4 45600 182250
P7 8 45600 364500
P8 12 45600 546750
P9 15 45600 683437,5
P10 0 45600 0
P11 4 45600 182250
P12 8 45600 364500
P13 12 45600 546750
P14 15 45600 683437,5
P15 0 45600 0
P16 4 45600 182250
P17 8 45600 364500
P18 12 45600 546750
P19 15 45600 683437,5
V1 12 4151882,81 49822593,75
V2 8 17968500 143748000
Total — — 199174781,3
Centre de rigidité Xr 8,72 m

TABLE 10 – Calcul du centre de rigidité Xr

Rapport de projet 10 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

Poteaux et voiles Y j (m) Kx j (N/mm) Kx j ·Y j

P1 0 45600 0
P2 0 45600 0
P3 0 45600 0
P4 0 45600 0
P5 3 11390,625 34171,875
P6 3 11390,625 34171,875
P7 3 11390,625 34171,875
P8 3 11390,625 34171,875
P9 3 11390,625 34171,875
P10 8 11390,625 91125
P11 8 11390,625 91125
P12 8 11390,625 91125
P13 8 11390,625 91125
P14 8 11390,625 91125
P15 13 11390,625 148078,125
P16 13 11390,625 148078,125
P17 13 11390,625 148078,125
P18 13 11390,625 148078,125
P19 13 11390,625 148078,125
V1 1.425 51257,8125 73042,38281
V2 5.5 883531,25 4594218,75
Total — — 18993439,258
Centre de rigidité Yr 3,89 m

TABLE 11 – Calcul du centre de rigidité Yr

4.4 Position du centre de rigidité – Étages hauts (3e et 4e )

Le même raisonnement est appliqué pour les deux étages supérieurs, où la rigidité
des poteaux est augmentée :

Poteaux et voiles X j (m) Ky j (N/mm) Ky j · X j

P1 0 27000 0
P2 4 27000 108000
P3 8 27000 216000
P4 12 27000 324000
P5 0 108000 0
P6 4 108000 432000
P7 8 108000 864000
P8 12 108000 1296000
P9 15 108000 1620000
P10 0 108000 0
P11 4 108000 432000
P12 8 108000 864000
P13 12 108000 1296000
P14 15 108000 1620000
P15 0 108000 0
P16 4 108000 432000
P17 8 108000 864000
P18 12 108000 1296000
P19 15 108000 1620000
V1 12 9841500 118098000
V2 8 42592000 340736000
Total — — 472118000
Centre de rigidité Xr 8,72 m

TABLE 12 – Calcul du centre de rigidité Xr basé sur les raideurs selon y

Rapport de projet 11 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

Poteaux et voiles Y j (m) Kx j (N/mm) Kx j ·Y j

P1 0 108000 0
P2 0 108000 0
P3 0 108000 0
P4 0 108000 0
P5 3 27000 81000
P6 3 27000 81000
P7 3 27000 81000
P8 3 27000 81000
P9 3 27000 81000
P10 8 27000 216000
P11 8 27000 216000
P12 8 27000 216000
P13 8 27000 216000
P14 8 27000 216000
P15 13 27000 351000
P16 13 27000 351000
P17 13 27000 351000
P18 13 27000 351000
P19 13 27000 351000
V1 1.425 121500 173137.5
V2 5.5 198000 1089000
Total — — 4502137.5
Centre de rigidité Yr 3,89 m

TABLE 13 – Calcul du centre de rigidité Yr à partir des raideurs selon x

5 Moments de torsion induits par les excentricités

Lorsque le centre de masse et le centre de rigidité ne coïncident pas, l’application
d’une force horizontale génère un moment de torsion. Ce moment est donné par :

Mt = F × e
— F : effort horizontal appliqué
— e : distance entre le centre de masse et le centre de rigidité (excentricité)
Afin de tenir compte de l’effet de torsion accidentelle, deux valeurs d’excentricité
sont calculées :

e1 = 0,5 · e + 0,05 · L e2 = 0,05 · L

F IGURE 3 – Illustration de l’excentricité horizontale

Rapport de projet 12 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

Cas d’un séisme selon l’axe Y

Dans ce cas, l’excentricité est mesurée selon l’axe X :

Mt1 = Fy · (ex + e1 ) Mt2 = Fy · (ex − e2 )

Étage Fy (kN) ex (m) e1 (m) e2 (m) Mt1 (kN·m) Mt2 (kN·m)

4 349,845 1,507 1,503589 0,75 1053,302512 264,8952629
3 237,008 1,507 1,503589 0,75 713,5763481 179,4574608
2 194,032 1,507 1,503589 0,75 584,1845866 146,9167004
1 97,016 1,507 1,503589 0,75 292,0922933 73,45835021

TABLE 14 – Moments de torsion pour un séisme dans la direction Y

Cas d’un séisme selon l’axe X

L’excentricité est cette fois-ci mesurée selon l’axe Y :

Mt1 = Fx · (ey + e1 ) Mt2 = Fx · (ey − e2 )

Étage Fx (kN) ey (m) e1 (m) e2 (m) Mt1 (kN·m) Mt2 (kN·m)

4 359,133 2,849 2,074518 0,65 1768,211069 789,7466297
3 240,300 2,849 2,074518 0,65 1183,130261 528,4285073
2 199,183 2,849 2,074518 0,65 980,6884509 438,0107173
1 99,592 2,849 2,074518 0,65 490,3466872 219,0064582

TABLE 15 – Moments de torsion pour un séisme dans la direction X

6 Vérification de la stabilité globale de la structure

Pour s’assurer que le bâtiment reste stable sous l’effet des forces horizontales, on
calcule un indice noté θs , qui compare les déplacements relatifs à la force sismique
appliquée.
K ·Ws · ∆x

θsx =

 Vx · h
K ·W s · ∆y
θsy =


Vy · h
— K : coefficient de comportement sismique
— Ws : poids cumulé des étages supérieurs
— ∆ : déplacement relatif mesuré
— Vx , Vy : efforts horizontaux
— h : hauteur de l’étage

Étage K Ws (kN) U (mm) V (mm) Fx (kN) Fy (kN) hs (m) θsx θsy

1 2 4570,02 1,847122013 0,038421206 901,207944 877,900842 4 0,004683372 0,000100003
2 2 3403,8 1,642998549 0,03417532 801,616425 780,884966 8 0,001744113 3,72417E-05
3 2 2237,58 0,52091084 0,010835247 602,433387 586,853213 11 0,000351779 7,51147E-06
4 2 1201,56 0,310534417 0,006459295 359,1330533 349,84513 14 0,000148423 3,16925E-06

TABLE 16 – Calcul des indices de stabilité θsx et θsy par étage

Rapport de projet 13 ENSAM-MEKNES

Travaux pratiques sur le calcul parasismique des structures

Interprétation des résultats

L’indice θs permet de juger si une structure reste stable sous l’effet d’un séisme.
Voici les seuils d’interprétation :
• Si θs ≤ 0,10, alors la structure est considérée comme **stable**. Aucune mesure
corrective n’est nécessaire.
• Si 0,10 < θs ≤ 0,20, les **effets du second ordre** doivent être pris en compte
dans l’analyse.
• Si θs > 0,20, la structure présente une **instabilité importante**. Un renforcement
est alors indispensable.
Dans notre cas θs est toujours inférieure à 0,10, Alors la structure est stable.

7 Simulation numérique avec Robot Structural Analy-

sis

8 Analyse comparative des résultats

Rapport de projet 14 ENSAM-MEKNES