CHAPITRE IV Etudes dynamique

OBJECTIF DE L'ETUDE DYNAMIQUE :

L’objectif initial de l’étude dynamique d’une structure est la détermination des caractéristiques
dynamiques propres de la structure lors de ces Vibrations Libres Non Amorties (VLNA).
L’étude dynamique d’une structure telle qu’elle se présente, est souvent très complexe. C’est pour cela
qu'on fait souvent appel à des modélisations qui permettent de simplifier suffisamment le problème pour
pouvoir l’analyser.

IV.1. PRESENTATION DU LOGICIEL DE CALCUL ETABS :

On a fait appel pour la modélisation de notre bâtiment à un logiciel appelé ETABS pour déterminer :
- Les périodes propres
- Les coefficients de participation.
- Les déplacements.
Les caractéristiques principales d’ETABS sont :
ETABS est un logiciel de calcul conçu exclusivement pour le calcul des bâtiments. Il permet de
modéliser facilement et rapidement tous types de bâtiments grâce à une interface graphique unique. Il offre
de nombreuses possibilités pour l’analyse statique et dynamique.
Ce logiciel permet la prise en compte des propriétés non-linéaires des matériaux, ainsi que le calcul et le
dimensionnement des éléments structuraux suivant différentes réglementations en vigueur à travers le monde
(Euro code, UBC, ACI.etc.). De plus de par sa spécificité pour le calcul des bâtiments, ETABS offre un
avantage certain par rapport aux codes de calcul à utilisation plus étendue. En effet, grâce à ces diverses
fonctions il permet une décente de charge automatique et rapide, un calcul automatique du centre de masse et
de rigidité, ainsi que la prise en compte implicite d’une éventuelle excentricité accidentelle. De plus, ce
logiciel utilise une terminologie propre au domaine du bâtiment (plancher, dalle, trumeau, linteau etc).
ETABS permet également le transfert de donnée avec d’autres logiciels (AUTOCAD, SAP2000 et SAFE).

IV.1.1 modélisation de la structure :

Une structure est composée de membrures flexibles reliées par des nœuds et supportée par une fondation
rigide. Cette définition implique que le modèle choisi n’est qu’une idéalisation mathématique de la structure
réelle. L’étude de la réponse dynamique d’une structure (période et mode propre) nécessite le choix d’un
modèle dynamique adéquat traduisant le plus concrètement possible la nature du système réel.
Le modèle dynamique que nous avons adopté pour le calcul de notre ouvrage est une console verticale
encastrée à sa base. Les masses seront considérées concentrées au niveau du centre de gravité de chaque
niveau.
IV.1.2 Etapes de modélisation :
Pour modéliser notre bâtiment nous avons considéré les étapes suivantes:

- Introduction de la géométrie du modèle (position des nœuds, connectivité des éléments).

- Spécification et assignement des caractéristiques mécaniques et géométriques aux différents éléments.
- Définition des différents cas de charges (statique et dynamique). Assignement de ces charges aux
nœuds et aux éléments.
- Spécification des conditions aux limites (appuis, diaphragmes,……).
- Lancement de l’analyse du problème, apporter des corrections au modèle s’il y a lieu.
- Visualisation et interprétation des résultats (à l’écran, sur fichier, etc.….).

U.A.M.O.BOUIRA (2016-2017) Page 69

CHAPITRE IV Etudes dynamique

IV .1.3 Proposition de la structure :

Dans cette partie, nous allons aborder l’analyse du comportement dynamique de deux structures cela afin
d’aboutir à une meilleur conception parasismique, dans ce volet, nous serons ramenés à déterminer les
périodes propres des différentes structures aux valeurs obtenues en appliquant les formules empirique du
RPA99V2003

a) structure mixte.

Fig (IV.1) : Modélisation de la structure mixte (Vue en 2D)

Fig (IV.2) : Modélisation de la structure mixte (Vue en 3D)

U.A.M.O.BOUIRA (2016-2017) Page 70

CHAPITRE IV Etudes dynamique

b) structure entierement en voiles.

Fig (IV.3) : Modélisation de la structure entièrement en voiles (Vue en 2D)

Fig (IV.4) : Modélisation de la structure entièrement voiles (Vue en 3D)

U.A.M.O.BOUIRA (2016-2017) Page 71

CHAPITRE IV Etudes dynamique

IV .1.4. Caractéristiques géométriques de la structure :

a) Centre de masse :

Le centre de masse est par définition le point d’application de la résultante de l’effort sismique. Les
coordonnées (XG, YG) du centre de masse sont données par les formules suivantes :
n n

 M i  YGi M i  X Gi
XG  i 1
n
YG  i 1
n

 Mi
i 1
M
i 1
i

Mi : Est la masse de l’élément "i".

Xi , Yi : Sont les cordonnées de l’élément considérée.

b) Centre de rigidité :
Le centre de torsion est le barycentre des rigidités des éléments dans le contreventement du bâtiment.
Autrement dit, c’est le point par lequel passe la résultante des réactions des voiles, des poteaux. En général
deux cas se présentes:

- Si le centre de masse et le centre de torsion sont confondus; les efforts horizontaux (séisme, vent…etc.) ne
provoquent sur la structure qu’une translation.
- Si par contre le centre de masse est excentré par rapport au centre de torsion, la structure subie donc une
translation et une rotation en même temps.

Dans le cas où il est procédé à une analyse tridimensionnelle, en plus de l'excentricité théorique calculée,
une excentricité accidentelle égale à ±0.05L, (L : la dimension du plancher perpendiculaire à la direction de
l'action sismique) doit être appliquée au niveau du plancher considéré et suivant chaque direction.

c) Excentricité accidentelle :

L’excentricité accidentelle est donnée par la formule suivante : eacc  0.05 L

L : La plus grande dimension de bâtiment

Sens X :eacc(x) = 0,05 . 28,45 =1,42m

Sens Y : eacc(y) = 0.05 . 18,4 =0.92m

IV .1.5. Calcul de la période fondamentale :

La valeur de la période fondamentale (T) de la structure peut être estimée à partir de formules empiriques
ou calculée par des méthodes analytiques ou numériques.
La formule empirique (4-6 donnée par le RPA99/version2003) est la suivante :

T= Ct hN 3/ 4
Avec:
- hN : hauteur mesurée en mètre a partir de la base de la structure jusqu’au dernier niveau (N) :
hN = 30.6 m
U.A.M.O.BOUIRA (2016-2017) Page 72
CHAPITRE IV Etudes dynamique

- CT : est un coefficient qui est fonction du système de contreventement et de type de remplissage ce

paramètre est donné par le tableau 4.6 du RPA99/version2003

Pour un contreventement assuré partiellement ou totalement par des voiles en béton armé: Tab4.6 (RPA99)
D’où : CT = 0.05
Donc : T = 0,05 × (30.6)3/4 = 0.65 sec
T =0,65 sec
N.B : on remarque que la même période fondamentale pour les deux structures.
IV .1.6. Nombre de modes à considérer :

 Le nombre de modes à prendre en compte est tel que la somme des coefficients de participation de
ces modes est au moins égale à 90% (  i  90%) , ou que tous les modes ayant une masse
modale effective supérieure à 5% de la masse totale de la structure .
 Le minimum de mode à retenir est de trois dans chaque direction considérée (Art 4.3.4
RPA 2003).

IV .2. ETUDE DE LA STRUCTURE MIXTE POTEAUX -VOILES :

La variante que nous avons analysée est une structure mixte (poteaux, voiles de contreventement 20cm),
disposées selon le schéma montré dans la figure ci-dessous.

Fig (IV.5) : schéma d’emplacement des voiles (Vue en 2D)

U.A.M.O.BOUIRA (2016-2017) Page 73

CHAPITRE IV Etudes dynamique

 Caractéristique massique :

Story Masse (ton) XCM YCM XCR YCR ex(M) ey(m)

STORY1 552.4512 14.035 9.134 14.225 8.787 0.19 -0.347
STORY2 547.3009 14.095 9.125 14.225 8.734 0.13 -0.391
STORY3 538.71 14.093 9.124 14.225 8.74 0.132 -0.384
STORY4 531.1975 14.091 9.123 14.225 8.746 0.134 -0.377
STORY5 523.4257 14.089 9.122 14.225 8.75 0.136 -0.372
STORY6 516.7323 14.087 9.121 14.225 8.754 0.138 -0.367
STORY7 509.7796 14.085 9.12 14.225 8.758 0.14 -0.362
STORY8 503.9053 14.084 9.119 14.225 8.762 0.141 -0.357
STORY9 498.9523 14.082 9.062 14.225 8.763 0.143 -0.299
STORY10 458.6016 14.225 9.059 14.225 8.757 0 -0.302
Poids (KN) 50826.1633

Tab (IV.1): Caractéristique massique de la structure mixte poteaux voiles

 Caractéristique dynamique :

Mode Période UX UY SumUX SumUY RZ

1 0.823682 72.9236 0 72.9236 0 0.0219
2 0.723525 0.0001 65.9988 72.9237 65.9988 0.0047
3 0.506606 0.0303 0.0042 72.954 66.003 66.053
4 0.262948 13.4139 0 86.3679 66.003 0.0053
5 0.169446 0 20.6928 86.3679 86.6957 0.0017
6 0.13466 5.3722 0 91.7401 86.6957 0.0001
7 0.120167 0.0006 0.0017 91.7407 86.6974 20.6966
8 0.083132 2.9974 0 94.7381 86.6974 0.0002
9 0.076762 0 7.062 94.7381 93.7595 0.0005
Tab (IV.2): Caractéristique dynamique de la structure mixte poteaux voiles
Commentaire :

D’après les résultats du tableau, on constate que :

Tdyn : la période obtenue par ETABS (Tdyn = 0.824 sec)
T : la période fondamentale calculée par les méthodes empiriques (TRPA = 0.65 sec)
- La période fondamentale Tdyn = 0.824 sec
- La période donnée par ETABS ne doit pas dépasser celle estimée à partir de la formule empirique de
plus de 30%, Tdyn < 1.3TRPA
Tdyn=0.824 sec < 1.3 (0.65) =0.845 sec  condition vérifiée.
- Le premier mode : C’est un mode de translation selon(x-x) avec un coefficient de participation
modale égale à 72.92 %

U.A.M.O.BOUIRA (2016-2017) Page 74

CHAPITRE IV Etudes dynamique

- le deuxième mode : C’est un mode de translation selon (y-y) avec un coefficient de participation
modale égale à 65.99 %.
Le troisième mode : C’est un mode de torsion de période avec un coefficient de participation modale égale à
0.03 % selon (x-x) et 0.004 % selon le sens (y-y)

 Schéma des modes de vibrations :

Fig (IV.6) : schéma de 1 er mode de vibration (Vue en 2D)

Fig (IV.7) : schéma de 2eme mode de vibration (Vue en 2D)

U.A.M.O.BOUIRA (2016-2017) Page 75

CHAPITRE IV Etudes dynamique

Fig (IV.8) : schéma de 3eme mode de vibration (Vue en 2D)

IV.2.1 étude sismique :

Généralité :

Les tremblements de terre présentent depuis toujours un des plus graves désastres pour l’humanité. Leur
apparition brusque et surtout imprévue, la violence des forces mises en jeu, l’énormité des pertes humaines et
matérielles ont marqué la mémoire des générations.
Parfois, les dommages sont comparables à un budget national annuel pour un pays, surtout si ses ressources
sont limitées.
Une catastrophe sismique peut engendrer une grave crise économique, ou au moins une brusque
diminution du niveau de vie dans toute une région pendant une longue période.
L’une des questions qui se posent est : Comment limiter les endommagements causés aux constructions
par le séisme ?
Il est donc nécessaire de construire des structures résistant aux séismes en zone sismique

IV.2.1.1 Introduction :
Selon le RPA 2003 notre structure est implantée en zone sismique (IIa) et vu les dégâts énormes que
peut engendrer un séisme, une étude sismique s’impose. Ainsi, le calcul d’un bâtiment vis à vis du séisme
repose sur l’évaluation des charges susceptibles d’être engendrées dans le système structural lors du séisme.
Dans le cadre de notre projet, la détermination de ces efforts est conduite par le logiciel ETABS en utilisant
une approche dynamique basée sur le principe de la superposition modale.

U.A.M.O.BOUIRA (2016-2017) Page 76

CHAPITRE IV Etudes dynamique

IV.2.1.2. Choix de la méthode de calcul :(R.P.A 2003)

Le calcul des forces sismiques peut être mené suivant trois méthodes :
 Par la méthode statique équivalente.
 Par la méthode d’analyse modale spectrale.
 Par la méthode d’analyse dynamique par accélérogrammes.
Dans le cas de notre bâtiment étudié on à l’irrégularité en plan, nous utiliserons la méthode d’analyse
modale spectrale pour l’analyse sismique. Cette dernière d’après le RPA99/version 2003 peut être utilisée
dans tous les cas et en particulier, dans le cas où la méthode statique équivalente n’est pas permise.

IV.2.1.3 Méthode d’analyse modale spectrale

a) Principe :

Rappelons que le principe de la méthode d’analyse modale spectrale est de rechercher pour chaque mode
de vibration, le maximum des effets engendrés dans la structure par les forces sismiques représentées par un
spectre de réponse de calcul. Ces effets sont par la suite combinés pour obtenir la réponse de la structure.
L'action sismique est simulée grâce à un spectre de réponse. Le comportement de la structure est supposé
élastique, permet le calcul des modes propres.

b) Spectre de réponse de calcul :

Le RPA99/2003 impose un spectre de réponse de calcul défini par la fonction suivante :

  T Q 
Cas1 : 1.25 A  1   2.5  1   Si 0  T  T1
  T1  R 
 Q
Cas2 : 2.5  1.25 A  Si T1  T  T2
Sa  R
 2/3
g  Q  T2 
Cas3 : 2.5   1.25 A    Si T2  T  3.0 s
 RT 


2/3 5/3
Q  T2   3 
Cas4 : 2.5  1.25 A  R  3   T  Si T  3.0 s
    

Avec :
- g : accélération de la pesanteur.
- A : coefficient d’accélération de zone, donné par le tableau (4.1) de RPA 99/version 2003 en fonction de la
zone sismique et du groupe d’usage du bâtiment
Dans notre cas, on est dans une Zone de type (IIa) et un Groupe d’usage 2.
Nous trouverons : A = 0.15
.
-  : Facteur de correction d’amortissement (quand l’amortissement est différent de 5%).
7
  0.7
(2   )
Où (%) est le pourcentage d’amortissement critique, fonction du matériau constitutif, du type de structure
et de l’importance des remplissages.
 est donnée par le tableau 4.2 du RPA99 :

U.A.M.O.BOUIRA (2016-2017) Page 77

CHAPITRE IV Etudes dynamique

Nous avons une structure mixte (portique et voiles) avec un remplissage dense, Donc  = 7 %
7
  0.882 0.7
9
- R : coefficient de comportement de la structure donnée par le tableau 4.3 du RPA99,
Dans notre cas on a un système de contreventement (mixte portiques/voiles avec interaction) en béton armé,
donc d’après le tableau R=5

- T1, T2 : périodes caractéristiques associées à la catégorie de site, elles sont données par le tableau 4.7 du
RPA99/ version 2003
Notre structure est construite dans un site meuble (S3) donc T1= 0.15sec, T2= 0.50sec
- Q : est le facteur de qualité il est fonction de :
 Conditions minimales sur les files de contreventement → non observée
 La redondance en plan. → non observée
 La régularité en plan. → non observée
 La régularité en élévation → observée
 Contrôle de la qualité des matériaux. → observée
 Contrôle de la qualité de l’exécution. → observée

 Sa valeur est déterminée par la formule : Q = 1 +  Pq

Avec Pq est la pénalité à retenir selon que le critère de qualité "q" est satisfait ou non.
À partir du tableau 4.4 de RPA 99 version 2003 on trouve : Q = 1.15

Pour le calcul du spectre de réponse on utilise le logiciel Spectre 99 V 2003

Fig (IV.9) : Spectre de réponse (structure mixte)

U.A.M.O.BOUIRA (2016-2017) Page 78

CHAPITRE IV Etudes dynamique

c) Calcul de l’effort tranchant à la base :

Sa
D’après le (RPA99/V 2003) : Vi    i  Wt
g

Avec : Vi : l’effort tranchant modal à la base.

 i : Coefficient de participation du mode i.

Wt : Poids total de la structure.

 Sens longitudinal (x-x) :

Mode Période Cas αi % Wt(KN) Sa/g Vi(KN)

1 0.823682 3 72.9236 50826.16 0.068 252037.01
2 0.723525 3 0.0001 50826.16 0.074 0.3761136
3 0.506606 3 0.0303 50826.16 0.094 144.76307
4 0.262948 2 13.4139 50826.16 0.095 64768.818
5 0.169446 2 0 50826.16 0.095 0
6 0.13466 1 5.3722 50826.16 0,104 28397.023
7 0.120167 1 0.0006 50826.16 0.113 3.4460136
8 0.083132 1 2.9974 50826.16 0.136 20719.101
9 0.076762 1 0 50826.16 0.14 0

Tab (IV.3): l’effort tranchant à la base sens (x-x)

 Sens transversal (y-y) :

Mode Période Cas αi % Wt(KN) Sa/g Vi(KN)

1 0.823682 3 0 50826.16 0.068 0
2 0.723525 3 65.9988 50826.16 0.074 248230.45
3 0.506606 3 0.0042 50826.16 0.094 20.066168
4 0.262948 2 0 50826.16 0.095 0
5 0.169446 2 20.6928 50826.16 0.095 99914.879
6 0.13466 1 0 50826.16 0,104 0
7 0.120167 1 0.0017 50826.16 0.113 9.7637053
8 0.083132 1 0 50826.16 0.136 0
9 0.076762 1 7.062 50826.16 0.14 50250.808

Tab (IV.4): l’effort tranchant à la base sens (y-y)

U.A.M.O.BOUIRA (2016-2017) Page 79

CHAPITRE IV Etudes dynamique

d) Combinaisons des réponses modales: (Art 4.3.5 ; RPA 2003)

1- Les réponses de deux modes de vibration i et j de périodes Ti, Tj d’amortissement i , j sont considérées
comme indépendantes si le rapport
Ti
r Avec Ti  T j
Tj
10
r Avec :  I   J  7%
10   I  J
2- Dans le cas où tous les réponses modales retenues sont indépendantes les unes des autres, la réponse totale
est donnée par
K
E =  E
i 1
i
2

Avec :
E : effet de l’action sismique considéré
Ei : valeur modale de E selon la mode « i »
K : nombre de modes retenues.

3- Dans le cas ou deux réponses modales ne sont pas indépendant, la réponse totale est donnée par :
K
E = ( E1  E2 )2   E i2
i 3

Les résultats sont représentés dans le tableau suivant :

Mode Période r

1 0.823682 0.878 0,993

2 0.723525 0.700 0,993
3 0.506606 0.519 0,993
4 0.262948 0.644 0,993
5 0.169446 0.795 0,993
6 0.13466 0.892 0,993
7 0.120167 0.692 0,993
8 0.083132 0.923 0,993
9 0.076762 0.878 0,993
Tab (IV.5): réponses modales
Commentaire

Toutes les réponses modales retenues sont indépendantes les unes des autres, donc la réponse totale est
K
donnée par : E E
i 1
i
2

Ex=2625.89 KN , EETABS =2627.59 KN

Ey=2722.62 KN , EETABS =2744.08 KN

U.A.M.O.BOUIRA (2016-2017) Page 80

CHAPITRE IV Etudes dynamique

IV.2.1.4. Méthode statique équivalente :

a) Principe :
Les forces réelles dynamiques qui se développent dans la construction sont remplacées par un système de
forces statique fictives dont les effets sont considérés équivalents à ceux de l’action sismique

b) Calcul de force sismique total par la méthode statique équivalant.

D’après l’art 4.2.3 de RPA99/version 2003, la force sismique totale V qui s’applique à la base de la
structure, doit être calculé successivement dans deux directions horizontales et orthogonales selon la
formule :
A.D.Q
V .W
R

- D : est un facteur d’amplification dynamique moyen qui est fonction de la catégorie de site, du facteur
d’amortissement () et de la période fondamental de la structure (T).

2.5 .............................0  T  T
 2

 2
 T 3
D  2.5  2  ...................T1  T  3.0 s
 T 
 2 5
  T2  3  3  3
2.5  3    T  ........T  3.00 s
  

Dans le cas où la structure est contreventée par des voiles en béton armé en plus de la formule (4.6 R.P.A
2003), on peut également utiliser aussi la formule 4.7 R.P.A 2003
0, 09hN
T 
D
- D : est la dimension du bâtiment mesurée à sa base dans la direction de calcul.
a) - sans longitudinale : Dx = 28.45 m TX = 0.51 sec
b) - sens transversale : Dy = 19.00 m TY = 0.65 sec

D’après RPA99/version 2003, il y a lieu de retenir dans chaque direction considérée la plus petite des
deux valeurs données respectivement par (4.6) et (4.7).
Donc: TX = min (T, Tx) = 0.51 sec
TY = min (T, Ty) = 0.65 sec

On a: T1 = 0.15 sec
T2 = 0.50 sec

2/3
 T2 
0.15 < TX < 3s donc D = 2.5  
T 
2/3
 T2 
0.15 < TY < 3s donc D = 2.5  
T 

U.A.M.O.BOUIRA (2016-2017) Page 81

CHAPITRE IV Etudes dynamique

D’où :
2
0,5 3
𝑫𝑿 = 2,5 . 0,882 ( ) = 2,18
0,51
0,5 2
𝑫𝒀 = 2,5 . 0.882( )3 = 1,85
0,65

W : poids total de la structure.

W est égal à la somme des poids Wi calculés à chaque niveau (i).

n
W  Wi Avec : Wi  WGi  0.2WQi (formule 4-5 R.P.A2003)
i 1

WGi : Poids du aux charges permanents.

WQi : Charge d’exploitation.

 : Coefficient de pondération, fonction de la nature et de la durée de la charge d’exploitation et donné

par le tableau (4.5 du RPA99 version 2003).
Dans notre cas, (le bâtiment à usage d’habitation)  =0,20.
Donc à chaque niveau : Wi  WGi  0.2WQi
La masse totale calculée par ETAPS :
WT=50826,1633 ton

Conclusion :

Paramètres A  Q Dx Dy R WT(KN) Vx(KN) Vy(KN)

Valeurs 0.15 0.882 1.15 2.18 1.85 5 50826,1633 3822,64 3244

c) Vérification de la résultante des forces sismiques de calcul :

D’après le RPA99/version2003 (art 4.3.6),
 la résultante Vt des forces obtenues à la base par combinaison des valeurs modales ne doit
pas être inférieure à 80% de la résultante des forces sismiques déterminés par la méthode
statique équivalente V.
 si Vt < 0.8V il faudra augmenter tous les paramètres de la réponse (force, déplacements,
moment …..)dans le rapport 0.8V/ Vt .
Ceci nous donne :
Vt (KN) V(KN) Vt/ V Observation
Sens longitudinal (x-x) 2625.89 3822,64 0.69 Non Vérifier

Sens transversal (y-y) 2722.62 3244 0.84 Vérifier

Tab (IV.6): Vérification de la résultante des efforts tranchants a la base
Commentaire :

Dans notre cas le sens (x-x) on doit calcule le rapport 0.8V/ Vt =1.16 → majoré (force, moment,
déplacement…) avec 16 %

U.A.M.O.BOUIRA (2016-2017) Page 82

CHAPITRE IV Etudes dynamique

IV.2.1.5. Vérifications :
a) Vérification des déplacements latéraux inter-étage :

Le RPA99/2003 exige de vérifier que les déplacements entres niveaux ne doivent pas dépasser 1% de la
hauteur d'étage.
Le déplacement horizontal à chaque niveau " k " de la structure est calculé comme suit.

 k  R .ek (Art 4 - 19 – R.P.A 2003).

 ek : Déplacement horizontal dû aux forces sismiques Fi (y compris l’effet de torsion).

R : Coefficient de comportement.
Le déplacement relatif au niveau " k " par rapport au niveau " k-1" est égal à :

k   k   k 1 (Art 4 – 20 RPA99/2003).
Les résultats des déplacements inter-étage sont présentés dans le tableau suivant :

 k (cm) R *  k (cm)  k (cm) 1% hk

NIV R hk (cm)
( 1,16)x  x y R *x R*  y x y (cm)

10 1.856 1.47 5 9.28 7.35 / / / /

9 1.72 1.29 5 8.6 6.45 0.68 0.9 306 3.06

8 1.55 1.1 5 7.75 5.5 0.85 0.95 306 3.06

7 1.357 0.91 5 6.78 4.55 0.97 0.95 306 3.06

6 1.148 0.72 5 5.74 3.6 1.04 0.95 306 3.06

5 0.916 0,54 5 4.58 2.7 1.16 0.9 306 3.06

4 0.68 0.38 5 3.4 1.9 1.18 0.8 306 3.06

3 0,46 0.24 5 2.3 1.2 1.1 0.7 306 3.06

2 0,24 0,12 5 1.2 0.6 1.1 0.6 306 3.06

1 0,079 0,04 5 0.39 0.2 0.81 0.4 306 3.06

Base 0 0 5 0 0 0.39 0.2 306 3.06

Tab (IV.7): vérification des déplacements inter-étage

Commentaire
Les résultats obtenus montrent que les déplacements relatifs inter-étages sont inférieurs à 1% de la
hauteur d’étage   k  1%hk

U.A.M.O.BOUIRA (2016-2017) Page 83

CHAPITRE IV Etudes dynamique

b) Vérification de l'effet P-  :

L'effet P-  représente l’augmentation des déplacements horizontaux dû à l’application simultanée de

charges verticales et d’une petite charge horizontale.
Le RPA99/2003 estime que les effets de seconde ordre peuvent être négligés si la condition suivantes
PK . K
est vérifié:   0,10 (article 5.6.RPA99/2003).
V K .hK
 Si 0,10 <  K < 0,20 : il faut augmenter les effets de l’action sismique calculés au moyen
1
d’une analyse élastique du 1° ordre par un facteur
(1   k )
 Si  K>0.20 : la structure est potentiellement instable et doit être redimensionnée.
PK : Poids total de la structure et des charges d'exploitation situés au-dessus du niveau K.
VK : effort tranchant d'étage au niveau " K ".
 K : déplacement relatif de niveau " K " par rapport au niveau " K-1 "
HK : hauteur d'étage.

Niv hk(cm) pk(KN) (1.16)×Vx(KN) Vy (KN) Δx(cm) Δy(cm) θx θy

10 306 4498.88 623.38 537.4 0.68 0.9 0.016 0.025

9 306 9393.6 1127.5 972 0.85 0.95 0.023 0.03
8 306 14291.9 1547.52 1334.07 0.97 0.95 0.029 0.033
7 306 19292.8 1911.76 1648 1.04 0.95 0.034 0.036
6 306 24362 2218.37 1912.4 1.16 0.9 0.042 0.037
5 306 29496.81 2485.08 2142.31 1.18 0.8 0.046 0.036
4 306 34707.86 2708.52 2334.93 1.1 0.7 0.046 0.034
3 306 39992.61 2878.23 2481.23 1.1 0.6 0.045 0.032
2 306 45361.63 2972.17 2584.5 0.81 0.4 0.04 0.023
1 306 50826.16 3048 2726.59 0.39 0.2 0.021 0.012

Tab (IV.8): vérification de l’effet P-Δ

Commentaire
La condition   0,10 est vérifiée.
Donc : l’effet P-Delta peut être négligé pour le cas de notre structure.

c) Justification vis à vis de l’équilibre d’ensemble:(RPA99/VERS03 Art.5.5.)

Cette condition d’équilibre se réfère à la stabilité d’ensemble du bâtiment ou de l’ouvrage, soumis à des
effets de renversement et/ou de glissement.

n
- M Re  M Fi / O  F d
i 1
i i

- M stab  MW / O  W  b

U.A.M.O.BOUIRA (2016-2017) Page 84

CHAPITRE IV Etudes dynamique

Fig (IV.10) : schéma isostatique des moments

Il faut vérifier que :

Moment stabilisateur
 1.5
Moment de renverssement

Wt(KN) b(m) Ms(KN.m) Mr(KN) Ms/Mr

Sens longitudinal 50826.16 14.22 722747.99 (1.16) x 57021.576 10.93

Sens transversal 50826.16 9.2 467600.67 54788.198 8.53

Tab (IV.9): vérification d’équilibre d’ensemble

Commentaire
L’ouvrage est donc stable au renversement (Ms > 1,5 Mr). De plus le poids des fondations et la butée par
les terres le rendent encore plus stable.

d) Vérification de l'effort normal réduit :

L'article (7-4-3-1) du R.P.A 2003 exige la vérification de l'effort normal réduit pour éviter la rupture
fragile de la section de béton.et comme cette vérification est vérifiée sous charges statiques donc on doit
refaire la vérification sous charge dynamique, La vérification s’effectue par la formule suivante :
N
Nrd =  0.3
B. f c 28

Où : N : l'effort normal maximal.

B : section du poteau.
Fc28 : résistance caractéristique du béton

U.A.M.O.BOUIRA (2016-2017) Page 85

CHAPITRE IV Etudes dynamique

Niveau Section (1.16)×N(KN) Nrd remarque

Story10 40x40 666.9768 0.167 Vérifier
Story7 45x45 1112.44 0.22 Vérifier
Story5 50x50 1601.032 0.256 Vérifier
Story3 55x55 2110.3532 0.28 Vérifier
Story1 60x60 2370.0424 0.26 Vérifier

Tab (IV.10): Vérification de l'effort normal réduit

IV .3. ETUDE DE LA STRUCTURE ENTIEREMENT VOILES :

La structure que nous avons analysée est entièrement en voiles, disposées selon le schéma montré dans la
figure ci-dessous

Fig (IV.11) : schéma d’emplacement des voiles (Vue en 2D)

U.A.M.O.BOUIRA (2016-2017) Page 86

CHAPITRE IV Etudes dynamique

 Caractéristique massique :

Masse XCM YCM XCR YCR ex ey

STORY1 560.2326 14.318 9.86 14.397 9.06 -0.079 0.8
STORY2 560.2326 14.318 9.86 14.395 9.011 -0.077 0.849
STORY3 560.2326 14.318 9.86 14.393 8.997 -0.075 0.863
STORY4 560.2326 14.318 9.86 14.392 9.006 -0.074 0.854
STORY5 560.2326 14.318 9.86 14.391 9.025 -0.073 0.835
STORY6 560.2326 14.318 9.86 14.39 9.048 -0.072 0.812
STORY7 560.2326 14.318 9.86 14.389 9.069 -0.071 0.791
STORY8 560.2326 14.318 9.86 14.389 9.086 -0.071 0.774
STORY9 556.9783 14.314 9.835 14.388 9.096 -0.074 0.739
STORY10 535.1342 14.345 9.911 14.388 9.096 -0.043 0.815
poids (KN) 54680.678

Tab (IV.11): Caractéristique massique de la structure entièrement voiles

 Caractéristique dynamique :

Mode Period UX UY SumUX SumUY RZ

1 0.337102 69.267 0.0005 69.267 0.0005 4.219
2 0.291319 4.6132 0.0245 73.8802 0.025 64.6225
3 0.239046 0.0004 69.6676 73.8806 69.6926 0.0245
4 0.095926 17.1712 0 91.0518 69.6926 0.2224
5 0.068762 0.2148 0.0117 91.2666 69.7043 19.6851
6 0.058584 0 19.4331 91.2666 89.1374 0.0122
7 0.049198 4.5628 0 95.8294 89.1374 0.0681
8 0.033748 2.0238 0 97.8533 89.1375 0.0198
9 0.031968 0.0001 0.0043 97.8534 89.1418 5.7532
10 0.027758 0 5.4937 97.8534 94.6355 0.0046
Tab (IV.12): Caractéristique dynamique de la structure entièrement voiles
Commentaire :

D’après les résultats du tableau, on constate que :

Tdyn : la période obtenue par ETABS (Tdyn = 0.317 sec)
T : la période fondamentale calculée par les méthodes empiriques (T RPA = 0.65 sec)
- La période fondamentale Tdyn = 0.337sec
- La période donnée par ETABS ne doit pas dépasser celle estimée à partir de la formule empirique de
plus de 30%, Tdyn < 1.3TRPA
Tdyn=0.337 sec < 1.3 (0.65) =0.845 sec  condition vérifiée.

U.A.M.O.BOUIRA (2016-2017) Page 87

CHAPITRE IV Etudes dynamique

IV.3.1 étude sismique :

IV.3.1.1 Méthode d’analyse modale spectrale

a) Spectre de réponse de calcul :

Le RPA99/2003 impose un spectre de réponse de calcul défini par la fonction suivante :

  T Q 
Cas1 : 1.25 A  1   2.5  1   Si 0  T  T1
  T1  R 
 Q
Cas2 : 2.5  1.25 A  Si T1  T  T2
Sa  R
 2/3
g  Q  T2 
Cas3 : 2.5   1.25 A    Si T2  T  3.0 s
 RT 


2/3 5/3
Q  T2   3 
Cas4 : 2.5  1.25 A  R  3   T  Si T  3.0 s
    
Avec :
 g : accélération de la pesanteur.
 A : coefficient d’accélération de zone, donné par le tableau (4.1) de RPA 99/version 2003 en
fonction de la zone sismique et du groupe d’usage du bâtiment
Dans notre cas, on est dans une Zone de type (IIa) et un Groupe d’usage 2.
Nous trouverons : A = 0.15
.
  : Facteur de correction d’amortissement (quand l’amortissement est différent de 5%).
7
  0.7
(2   )
Où (%) est le pourcentage d’amortissement critique, fonction du matériau constitutif, du type de
structure et de l’importance des remplissages.
 est donnée par le tableau 4.2 du RPA99 :
Nous avons une structure en voile (voiles en béton arme), Donc  = 10 %
η =0.764
 R : coefficient de comportement de la structure donnée par le tableau 4.3 du RPA99,
Dans notre cas on a un système de contreventement (voiles porteurs) en béton armé, donc d’après le tableau
R=3.5
 T1, T2 : périodes caractéristiques associées à la catégorie de site, elles sont données par le tableau 4.7
du RPA99/ version 2003
Notre structure est construite dans un site meuble (S3) donc T1= 0.15sec, T2= 0.50sec

 Q : est le facteur de qualité il est fonction de :

- Conditions minimales sur les files de contreventement → non observée
- La redondance en plan. → non observée
- La régularité en plan. → non observée
- La régularité en élévation → observée
- Contrôle de la qualité des matériaux. → observée
- Contrôle de la qualité de l’exécution. → observée

U.A.M.O.BOUIRA (2016-2017) Page 88

CHAPITRE IV Etudes dynamique

Sa valeur est déterminée par la formule : Q = 1 +  Pq

Avec Pq est la pénalité à retenir selon que le critère de qualité "q" est satisfait ou non.
À partir du tableau 4.4 de RPA 99 version 2003 on trouve : Q = 1.15.
Pour le calcul du spectre de réponse on utilise le logiciel Spectre 99 V 2003

Fig (IV.12) : Spectre de réponse (structure entièrement voiles)

b) Calcul de l’effort tranchant à la base :

 Sens longitudinal (x-x) :

Mode Période Cas αi % Wt(KN) Sa/g Vi(KN)

1 0.337102 2 0.69267 54680.678 0.115 4355.7015

2 0.291319 2 0.046132 54680.678 0.115 290.09084

3 0.239046 2 0 54680.678 0.115 0

4 0.095926 1 0.171712 54680.678 0.1428 1340.7961

5 0.068762 1 0.002148 54680.678 0.1555 18.264112

6 0.058584 1 0 54680.678 0.16 0

7 0.049198 1 0.045628 54680.678 0.165 411.67005

8 0.033748 1 0.020238 54680.678 0.172 190.33994

9 0.031968 1 0 54680.678 0.172 0

10 0.027758 1 0 54680.678 0.175 0

Tab (IV.13): l’effort tranchant à la base sens (x-x)

U.A.M.O.BOUIRA (2016-2017) Page 89

CHAPITRE IV Etudes dynamique

 Sens transversal (y-y) :

Mode Période Cas αi % Wt(KN) Sa/g Vi(KN)

1 0.337102 2 0.000005 54680.678 0.115 0.0314414

2 0.291319 2 0.000245 54680.678 0.115 1.5406281

3 0.239046 2 0.696676 54680.678 0.115 4380.8923

4 0.095926 1 0 54680.678 0.1428 0

5 0.068762 1 0.000117 54680.678 0.1555 0.9948329

6 0.058584 1 0.194331 54680.678 0.16 1700.1841

7 0.049198 1 0 54680.678 0.165 0

8 0.033748 1 0 54680.678 0.172 0

9 0.031968 1 0.000043 54680.678 0.172 0.4044183

10 0.027758 1 0.054937 54680.678 0.175 525.69867

Tab (IV.14): l’effort tranchant à la base sens (y-y)

c) Combinaisons des réponses modales: (Art 4.3.5 ; RPA 99)

1- Les réponses de deux modes de vibration i et j de périodes Ti, Tj d’amortissement i , j sont considérées
comme indépendantes si le rapport
Ti
r Avec Ti  T j
Tj
10
r Avec :
ξi = ξj=10%
10   I  J
2- Dans le cas où tous les réponses modales retenues sont indépendantes les unes des autres, la réponse
totale est donnée par
K
E =  E i 1
i
2

Avec :
E : effet de l’action sismique considéré
Ei : valeur modale de E selon la mode « i »
K : nombre de modes retenues.
Ex=4589.12kn Exetabs=4747.47KN
Ey=4728.55kn Eyetabs=4887.54kn

U.A.M.O.BOUIRA (2016-2017) Page 90

CHAPITRE IV Etudes dynamique

Les résultats sont représentés dans le tableau suivant :

mode Période r
1 0.337102 0.86 0,99
2 0.291319 0.82 0,99
3 0.239046 0.40 0,99
4 0.095926 0.72 0,99
5 0.068762 0.85 0,99
6 0.058584 0.84 0,99
7 0.049198 0.69 0,99
8 0.033748 0.95 0,99
9 0.031968 0.87 0,99
10 0.027758 0.86 0.99
Tab (IV.15): réponses modales

IV.3.1.2. Méthode statique équivalente :

a) Calcul de force sismique total par la méthode statique équivalant.
On a
a) - sans longitudinale : Dx = 28.8 m TX = 0.51 sec
b) - sens transversale : Dy = 19 m TY = 0.62 sec

D’après RPA99/version 2003, il y a lieu de retenir dans chaque direction considérée la plus petite des deux
valeurs données respectivement par (4.6) et (4.7).
Donc: TX = min (T, Tx) = 0.51 sec
TY = min (T, Ty) = 0.62 sec

On a: T1 = 0.15 sec
T2 = 0.50 sec
2/3
T 
0.15 < TX < 3s donc D = 2.5  2 
T 
2/3
T 
0.15 < TY < 3s donc D = 2.5  2 
T 
D’où :
2
0,5 3
𝑫𝑿 = 2,5 . 0,882 ( ) = 2,18
0,51
0,5 2
𝑫𝒀 = 2,5 . 0.882( )3 = 1,85
0,65

W : poids total de la structure.

W : est égal à la somme des poids Wi calculés à chaque niveau (i).

n
W  Wi Avec : Wi  WGi  0.2WQi (formule 4-5 R.P.A2003)
i 1

WGi : Poids du aux charges permanents.

U.A.M.O.BOUIRA (2016-2017) Page 91

CHAPITRE IV Etudes dynamique

WQi : Charge d’exploitation.

 : Coefficient de pondération, fonction de la nature et de la durée de la charge d’exploitation et donné

par le tableau (4.5 du RPA99 version 2003).
Dans notre cas, (le bâtiment à usage d’habitation)  =0,20.
Donc à chaque niveau : Wi  WGi  0.2WQi
La masse totale calculée par ETAPS :

WT=54680.678 KN

Conclusion :

Paramètres A  Q Dx Dy R WT(KN) Vx(KN) Vy(KN)

Valeurs 0.15 0.764 1.15 2.18 1.85 3.5 54680.678 5875.05 4985.7

IV.3.1.3. Vérifications :
a) Vérification de la résultante des forces sismiques de calcul :

D’après le RPA99/version2003 (art 4.3.6),

 la résultante Vt des forces obtenues à la base par combinaison des valeurs modales ne doit
pas être inférieure à 80% de la résultante des forces sismiques déterminés par la méthode
statique équivalente V.
 si Vt < 0.8V il faudra augmenter tous les paramètres de la réponse (force, déplacements,
moment …..)dans le rapport 0.8V/ Vt .
Ceci nous donne :

Vt (KN) V(KN) Vt/ V Observation

Sens longitudinal (x-x) 4747.47 5875.05 0.81 Vérifiée
Sens transversal (y-y) 4887.54 4987 0.94 Vérifiée

Tab (IV.16): Vérification de la résultante des efforts tranchants a la base

b) Vérification des déplacements latéraux inter-étage :

Le RPA99/2003 exige de vérifier que les déplacements entres niveaux ne doivent pas dépasser 1% de la
hauteur d'étage.
Le déplacement horizontal à chaque niveau " k " de la structure est calculé comme suit.

 k  R .ek (Art 4 - 19 – R.P.A 2003).

 ek : Déplacement horizontal dû aux forces sismiques Fi (y compris l’effet de torsion).

R : Coefficient de comportement.
Le déplacement relatif au niveau " k " par rapport au niveau " k-1" est égal à :

U.A.M.O.BOUIRA (2016-2017) Page 92

CHAPITRE IV Etudes dynamique

k   k   k 1 (Art 4 – 20 RPA99/2003).
Les résultats des déplacements inter-étage sont présentés dans le tableau suivant :

 k (cm) R *  k (cm)  k (cm)

hk
NIV R 1% hk (cm)
x y R *x R*  y x y (cm)

10 0.4502 0.2429 3.5 1.5757 0.85015 / / / /

9 0.4064 0.2155 3.5 1.4224 0.75425 0.153 0.096 306 3.06

8 0.3608 0.1868 3.5 1.2628 0.6538 0.159 0.100 306 3.06
7 0.3118 0.1574 3.5 1.0913 0.5509 0.171 0.103 306 3.06
6 0.2606 0.1278 3.5 0.9121 0.4473 0.179 0.1036 306 3.06
5 0.2085 0.0989 3.5 0.72975 0.34615 0.180 0.101 306 3.06
4 0.1571 0.0714 3.5 0.54985 0.2499 0.179 0.096 306 3.06
3 0.1082 0.0464 3.5 0.3787 0.1624 0.171 0.087 306 3.06
2 0.0638 0.025 3.5 0.2233 0.0875 0.155 0.075 306 3.06
1 0.0262 0.0087 3.5 0.0917 0.03045 0.131 0.057 306 3.06
Base 0.4502 0.2429 3.5 1.5757 0.85015 0.092 0.030 306 3.06

Tab (IV.17): vérification des déplacements inter-étage

Commentaire
Les résultats obtenus montrent que les déplacements relatifs inter-étages sont inférieurs à 1% de la
hauteur d’étage   k  1%hk

c) Vérification de l'effet P-  :

L'effet P-  représente l’augmentation des déplacements horizontaux dû à l’application simultanée de

charges verticales et d’une petite charge horizontale.
Le RPA99/2003 estime que les effets de seconde ordre peuvent être négligés si la condition suivantes
PK . K
est vérifié:   0,10 (article 5.6.RPA99/2003).
V K .hK
 Si 0,10 <  K < 0,20 : il faut augmenter les effets de l’action sismique calculés au moyen
1
d’une analyse élastique du 1° ordre par un facteur
(1   k )
 Si  K>0.20 : la structure est potentiellement instable et doit être redimensionnée.
PK : Poids total de la structure et des charges d'exploitation situés au-dessus du niveau K.
VK : effort tranchant d'étage au niveau " K ".
 K : déplacement relatif de niveau " K " par rapport au niveau " K-1 "
HK : hauteur d'étage.

U.A.M.O.BOUIRA (2016-2017) Page 93

CHAPITRE IV Etudes dynamique

Niv hk(cm) pk(KN) Vx(KN) Vy (KN) Δx(cm) Δy(cm) θx θy

10 306 5249.66 980.87 1084.78 0.153 0.096 0.0027 0.0015
9 306 10713.62 1803.53 1960.4 0.159 0.100 0.0031 0.0018
8 306 16209.50 2464.54 2624.65 0.171 0.103 0.0037 0.0021
7 306 21705.39 3013.62 3154.31 0.179 0.1036 0.0042 0.0023
6 306 27201.27 3473.31 3604.82 0.180 0.101 0.0046 0.0025
5 306 32697.15 3861.56 3990.96 0.179 0.096 0.0050 0.0026
4 306 38193.03 4188.7 4318.77 0.171 0.087 0.0051 0.0025
3 306 43688.91 4450.71 4592.98 0.155 0.075 0.0050 0.0023
2 306 49184.79 4643.93 4794.08 0.131 0.057 0.0045 0.0019
1 306 54680.68 4747.47 4887.54 0.092 0.030 0.0034 0.0011

Tab (IV.18): vérification de l’effet P-Δ

Commentaire
La condition   0,10 est vérifiée.
Donc : l’effet P-Delta peut être négligé pour le cas de notre structure.

d) Justification vis à vis de l’équilibre d’ensemble:(RPA99/VERS03 Art.5.5.)

Cette condition d’équilibre se réfère à la stabilité d’ensemble du bâtiment ou de l’ouvrage, soumis à des
effets de renversement et/ou de glissement.

n
- M Re  M Fi / O  F d
i 1
i i

- M stab  MW / O  W  b

Il faut vérifier que :

Moment stabilisateur
 1.5
Moment de renverssement

Wt(KN) b(m) Ms(KN.m) Mr(KN.m) Ms/Mr

Sens longitudinal 54680.678 14.22 777559.24 102432.77 7.59

Sens transversal 54680.678 9.2 503062.24 99367.754 5.06

Tab (IV.19): vérification d’équilibre d’ensemble

Commentaire
L’ouvrage est donc stable au renversement (Ms > 1,5 Mr). De plus le poids des fondations et la butée par
les terres le rendent encore plus stable.

U.A.M.O.BOUIRA (2016-2017) Page 94