Chapitre 4 Cour
Chapitre 4 Cour
Chapitre 4 Cour
I
I.1 Introduction
L’étude dynamique consiste à déterminer les réponses de la structure sous les charges
sismiques.
À l’aide de cette étude nous pouvons connaître le comportement de la structure vis-à-vis
l’excitation sismique.
Toutes les structures sont susceptibles d’être soumises pendant leur durée de vie à des
chargements variables dans le temps. Ces actions peuvent être dues au vent, séisme ou aux
vibrations des machines, ce qui signifie que le calcul sous l’effet des chargements statiques
parait insuffisant, d’où la nécessité d’une étude dynamique qui nous permet de déterminer les
caractéristiques dynamiques d’une structure afin de prévoir son comportement (déplacement,
effort normauxet période…) sous l’effet du séisme.
I.2 Etude dynamique
I.2.1 Objectif de l’étude dynamique :
a) Modélisation de la rigidité:
La modélisation des éléments constituants le contreventement (rigidité) est effectué
comme suit :
Chaque poutre et chaque poteau a été modélisé par un élément fini de type
" Frame "(Poteaux- Poutres) à deux nœuds (6 d.d.l par nœud)
Les voiles sont modélisés par des éléments " Shell "à quatre nœuds.
b) Modélisation de la masse :
La masse volumique attribuée au matériau constituant les poteaux et les poutres est prise
égale à celle du béton armé à savoir 2,5t/m3.
La masse de l’acrotère et des murs extérieurs (maçonnerie) a été repartie linéairement
aux niveaux des poutres correspondants à leur emplacement.
I
possibilités pour l’analyse statique et dynamique, il est basé sur la méthode des éléments
finis.
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,…….etc).
De plus, de part sa spécificité pour le calcul des bâtiments, ETABS offre un avantage
certain par rapport au codes de calcul à utilisation plus étendue. En effet, grâce à ces
diverses fonctions il permet une descente de charges 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). Il permet également le transfert de données avec d’autres
logiciels (AUTOCAD, SAP2000).
3éme étape :
Spécification des conditions aux limites (appuis, encastrement… etc.) pour la structure
à modéliser.
4éme étape :
5éme étape :
6éme étape :
Exécution et analyse des résultats en choisissent le nombre de modes propre à prendre en
considération.
7éme étape :
Création des fichiers de sorties et exploitation des résultats pour la vérification des
élément et calcul de ferraillage
100
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
I.2.5 Modélisation de la vue en 3D
101
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
L’excentricité
Théorique
Est calculée dans les deux sens par les formules suivantes :
102
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
e x =| X CM −X CR| ; e y =|Y CM −Y CR|
Accidentelle
La valeur exigée par le RPA 99 / version 2003 ± 0,05 L.
L : Etant la dimension du plancher perpendiculaire à la direction de l’action sismique, cette
valeur doit être appliquée au niveau du plancher considéré suivant chaque direction.
e acc =0 , 05(L ¿ ¿ x ; L y ) ¿.
e accX =0 . 05 ×37 .5=1. 875 ;e accY =0. 05 ×14 =0.7
Caractéristiques massiques
Les résultats des centres de masses, de torsion, et d’excentricités sont résumés dans le tableau
ci-dessous :
Story Masse (t) XCM(m) YCM(m) XCR(m) YCR(m) ex(m) ey(m)
RDC 866,05417 18,5992 6,7844 18,5987 6,6798 0,0005 0,1046
Etage1 803,05936 18,6003 6,7763 18,5994 6,7042 0,0009 0,0721
Etage2 795,56215 18,6003 6,7756 18,6003 6,7348 0 0,0408
Etage3 787,6053 18,6003 6,7749 18,6008 6,7641 0,0005 0,0108
Etage 4 781,04419 18,6003 6,7743 18,6012 6,7907 0,0009 0,0164
Etage 5 774,02345 18,6003 6,7736 18,6015 6,814 0,0012 0,0404
Etage6 765,08505 18,5974 6,7581 18,6016 6,8301 0,0042 0,072
Etage7 706,79308 18,6 6,8894 18,6014 6,8293 0,0014 0,0601
∑ Mi 6279,22675
Tableau IV-1: Centre de masse, de torsion, et l’excentricité des différents niveaux
I
1 0,711 75,97 0.00037 0,00391 75,97 0,00037 0,00391
2 0,627 0,0003492 79.45 0 75,98 79,45 0,00391
3 0,461 0,003506 0 79,98 75,98 79,45 79,98
4 0,192 16,24 0.0000877 0 92,22 79,46 79,98
5 0,178 0,00007078 13.2 0 92,22 92,62 79,98
6 0,131 0,00005026 0 13,01 92,22 92,62 92,99
7 0,085 4,6 0.0000594 0,0003085 96,82 92,62 92,99
8 0,083 0 4.2 0 96,82 96,82 92,99
9 0,062 0,0001941 0 4,12 96,82 96,82 97,11
104
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
105
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
106
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
Figure IV-5: Mode de rotation. Auteure de l’axe Z-Z
Commentaire
Ce modèle présente une période fondamentale T = 0,711 sec.
Les deux premiers modes sont des modes de translation.
Le 3ème mode est un mode de rotation.
On doit retenir les 9 premiers modes, pour que le coefficient de participation massique
atteigne les 90% (selon le RPA99).
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,3 TRPA.
Tdyn= 0,711 sec < 1,3 (0,59) = 0,767 sec…………..............................…Condition vérifiée.
Niveau 1 F1 F2 F1 F2 F1 F2 F1 F2
1926.239 85.38 18.2 2.26 99.05 97.34 0.95 2.66
Niveau 1 F3 F3 F3 F3
7273.55 1673.86 76.98 23.01
Avec :
F1 : La charge horizontale suivant (x).
F2 : La charge horizontale suivant (y).
F3 : La charge verticale reprise par l’élément.
Conclusion
Les voiles reprennent plus de 20% des sollicitations dues aux charges verticales. On
considère que la sollicitation horizontale est reprise uniquement par les voiles (d’après
RPA99/V2003). Ce qui implique que la structure a un coefficient de comportement R= 3.5.
I.3 Etude sismique
I.3.1 Méthode de calcul
Selon les règles parasismiques Algériennes (RPA99/version2003) le calcul des forces
sismiques peut être mené suivant trois méthodes :
La méthode statique équivalente.
107
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
La méthode d’analyse modale spectrale.
La méthode d’analyse dynamique par accélérogramme.
I.3.2 Méthode statique équivalente (RPA)
Dans cette méthode les forces réelles dynamiques qui se développent dans la
construction sont remplacées par un système de forces statiques fictives dont les effets sont
considérés équivalents à ceux de l’action sismique.
I.3.3 Méthode d’analyse modale spectrale (RPA)
La méthode d’analyse modale spectrale peut être utilisée dans tous les cas, et en
particulier, dans le cas où la méthode statique équivalente n’est pas permise.
Par cette méthode, il est recherché 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.
I.3.4 Méthode d’analyse dynamique par accélérogramme(RPA)
La méthode d’analyse dynamique par accélérogrammes peut être utilisée au cas par cas
par un personnel qualifié, ayant justifié auparavant les choix des séismes de calcul et des lois
de comportement utilisées ainsi que la méthode d’interprétation des résultats et les critères de
sécurité à satisfaire.
I.4 Choix de la méthode de calcul
Le choix de la méthode de calcul ce fait en vérifiant les conditions d’application de ces
méthodes de calcul.
I.4.1 Conditions d’application de la méthode statique équivalente
La méthode statique équivalente peut être utilisée dans les conditions suivantes :
Le bâtiment ou bloc étudié, satisfaisait aux conditions de régularité en plan et en
élévation avec une hauteur au plus égale à 65m en zones I et II et à 30m en zones III.
Le bâtiment ou bloc étudié présente une configuration irrégulière tout en respectant,
outres les conditions précédentes, les conditions complémentaires suivantes :
Zone I : Tous les groupes d’usages.
Zone II : Groupe d’usage 3.
Groupe d’usage 2…………….…..…… H T ≤ 7 Niveaux ou 23 m.
Groupe d’usage 1B…………….….….. H T ≤ 5 Niveaux ou 17 m.
Groupe d’usage 1A…………………... H T ≤ 3Niveaux ou 10 m.
Zone III : Groupe d’usage 2 et 3………………… H T ≤ 5 Niveaux ou 17 m.
Groupe d’usage 1B………………..….. H T ≤ 3 Niveaux ou 10 m.
Groupe d’usage 1A………………….... H T ≤ 2 niveaux ou 08 m.
Dans notre cas, notre structure est implantée à wilaya de Sétif zone II a de groupe d’usage 1B,
avec une hauteur de 26.69 m qui est inférieure à 65 m selon la 1 ère condition, et supérieure à
17 m
Selon la 2ème condition, ce qui signifie que les conditions d’application de cette méthode
ne sont pas vérifiées, donc il faut passer à la méthode d’analyse modale spectrale en utilisant
le logiciel de calcul ETABS 2016.
Classification de l’ouvrage
Selon le rapport géotechnique de notre ouvrage et d’après la classification du sol du
RPA, notre sol est un sol meuble (catégorie S3).
Zone sismique : Zone IIa «Setif»
108
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
Groupe d’usage : groupe 1B
Site : « meuble S3 » les valeur de T1 et T2 sont respectivement 0,15et 0,5
Système structurel : voiles-portiques.
{
suivant :
( (
1 , 25 A 1+
T
T1
Q
))
2 , 5 SYMBOL 104 ¿ Symbol ¿ 12 h −1 … … … … ….0 ≤ T ≤ T 1
R
( )
Q
S A 2 ,5 SYMBOL 104 ¿ Symbol ¿ 12h ( 1 , 25 A ) R … … … … … … … … … … .. T 1 ≤T ≤T 2
=
( )( )
g 2
Q T2 3
2 ,5 SYMBOL 104 ¿ Symbol ¿ 12h ( 1 , 25 A ) … … … … … … … … .T 1 ≤ T ≤ 3 s
R R
( )( ) ( )
2/ 3 5 /3
Q T2 3
2 , 5 SYMBOL 104 ¿ Symbol ¿ 12 h ( 1, 25 A ) … … … . … … … . T >3 s
R 3 T2
Avec :
A : Coefficient d’accélération de zone (RPA 99. Tableau 4.1).
hh : Facteur de correction d’amortissement (quand l’amortissement est différent de 5% égale
à :¿
√ 7
(2+ SYMBOL 120 ¿Symbol ¿ 12 x)
≥ 0 , 7).
Nous avons un contreventement mixte voiles-portiques donc on prend xx=10%.
xx: Pourcentage d’amortissement critique (RPA99 tableau (4-7)).
R : Coefficient de comportement de la structure.
T1, T2 : Périodes caractéristiques associées à la catégorie de site (RPA art 4.7).
Q : Facteur de qualité (RPA 99. Tableau44).
g : Accélération de la pesanteur.
I.4.3.1 Coefficient d’accélération de zone A
Il est donné par le tableau (4.1 R.P. A V2003) suivant la zone sismique et le groupe
d'usage du bâtiment. Pour la zone IIa et groupe d’usage 1B on aura : A = 0,20.
109
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
I.4.3.2 Calcul du facteur d’amplification dynamique moyen « D »
{
2 , 5 η … … … … … … … … … … . 0 ≤ T ≤ T2
( )
2
T
2 , 5 SYMBOL 104 ¿ Symbol ¿ 12 h 2 3 … . … … .. … … T 2 ≤T ≤3 s
D= T
Ly 14 T y 0.59 D y 1.7
I
6
Q=1+ ∑ P q
q=1
Avec :
Pq : La pénalité à retenir selon que le critère de qualité (q) est satisfait ou non sa valeur est
donnée dans tableau (4.4 R.P. A V2003).
Le tableau ci-dessous présente les valeurs de pénalités ainsi que le facteur de qualité :
Critère Pq
Observé Non
observé
Condition minimale sur les files de contreventement / 0.05
Redondance en plan / 0.05
Régularité en plan Obs /
Régularité en élévation Obs /
Contrôle de la qualité des matériaux / 0.05
Contrôle de la qualité d’exécution / 0.05
Q (totale) 1,2
SchémaParamètres
de spectre de calcul A Dx Dy Q R W ET(KN)
Valeurs numériques 0.2 1.9 1.7 1,20 3.5 71014.17
111
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
I.4.4 Vérification de la résultante des forces sismiques de calcul
La résultante des forces sismiques à la base V t obtenue par combinaison des valeurs
modales doit vérifier la condition :V t ≥ 0 , 80 V pour une période fondamentale donnée par la
formule empirique appropriée.
Si 0.8 V stat <V daynil faudra augmenter tous les paramètres de la réponse (forces,
0 , 8V
déplacements, moments,...) dans le rapport .
V dyn
Tableau IV- 8: Vérification de la force sismique statique et dynamique
I
considérés comme indépendant si le rapport :
Ti
r= (AvecT i ≤ T j) vérifie la relation :
Tj
10
r ≤ Avec : ξi= ξj = 10%......... r ≤ 0 ,5 .
¿¿
Dans le cas où toutes les réponses modales retenues sont indépendantes les unes des
autres, la réponse totale est donnée par :
√∑
K
E=± Ei²
i=1
Pour :
E : effet de l’action sismique considérée.
Ei : Valeur modale de E selon le mode « i ».
K : nombre de mode retenus.
Dans le cas ou deux réponses modales ne sont pas indépendant, E1, E2 par exemple, la
réponse totale est donnée par :
√(| | | |)
K
E1 + E 2 + ∑ E i ² .
2
E=
i=1
Remarque
Les modes ne sont pas indépendants deux à deux donc on applique la loi :
√(| | | |)
K
E1 + E 2 + ∑ E i ² . E x =¿322058.93 KN , E y =¿193060.71KN
2
E=
i=3
I
Pour chaque niveau « K » est donné par :
δ k =R δ ek (Art 4 – 20 RPA99/2003).
Avec :
δ ek : Déplacement dû aux forces sismiques Fi.
R : Coefficient de comportement.
Déplacement relatif : au niveau « K » par rapport au niveau « K-1 » est donnée par :
∆ K =δ k −δ k−1(Art 4 – 20 RPA99/2003).
Déplacement admissible
∆ K =0 , 01 h0
Avec :
h 0: Hauteur d’étage.
Tableau IV- 12: Vérification des déplacements inter-étage dans le sens Y-Y
I
5 5 5 5 6
Niveau2 0,00716 3. 0,025091 0,015963 3.06 0,030 CV
9 5 5 5 0,009128 6
Niveau1 0,00456 3. 0,015963 0,007591 3.06 0,030 CV
1 5 5 5 0,008372 6
RDC 0,00216 3. 0,007591 0,007591 5.27 0,052 CV
9 5 5 / 5 7
Conclusion
Nous constatons que les déplacements inter-étage ne dépassent pas le déplacement
admissible donnée par le RPA99/V2003, donc la stabilité de l’ouvrage sous charges
horizontales est assurée.
{
F t= 0 , 07 VT siT >0 , 7 sec
0 siT <0 , 7 sec
(V −F t )W i hi
F i= n
∑ W j hj
j=1
La force concentrée Ft au sommet de la structure permet de tenir compte de l’influence
des modes supérieurs de vibration.
Elle doit être déterminée par la formule : Ft = 0,07 TV où T est la période fondamentale
de la structure (en secondes).
La valeur de Ftne dépassera en aucun cas 0,25 V et sera prise égale à 0 quand T est plus
petit ou égale à 0,7 secondes.
Avec :
T : la période fondamentale de la structure.
F i : Effort horizontal revenant au niveau i.
hi : Niveau de plancher où s’exerce la force F i.
h j : Niveau de plancher quelconque.
W i , W j : Poids revenant au plancher i, j.
Notre cas : { {
T x =0,484 s< 0 ,7 s Ftx =0 KN .
T y =0,493 s <0 , 7 s . Fty =0 KN .
Les résultats des forces sismiques sont présentés dans le tableau suivant :
Tableau IV- 13: Distribution de la résultante de la force sismique selon la hauteur (sens X-X)
I
Niveau1 92338,7808 8,33 769182,044 9252,13 5841124,28 7116572265 1218,35659
RDC 105989,4171 5,27 558564,228 9252,13 5841124,28 5167908852 884,745574
SOMME 9252,12999
Tableau IV- 14: Distribution de la résultante de la force sismique selon la hauteur (sens Y-Y)
116
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
Sens-Y-Y « transversale »
Tableau IV- 16: Vérification de l’effet (P-∆) sens Y-Y
Conclusion
On remarque que θ Kxet θ Kysont inférieures à 0,1 ce qui signifie que l’effet P-∆ peut être
négligé pour le cas de notre ouvrage.
Vérification au renversement
Le moment de renversement qui peut être causé par l’action sismique doit être calculé
par rapport au niveau de contact du sol fondation.
Le moment de stabilisation sera calculé en prenant en compte le poids total équivalent au
poids de la construction ( M s > M r ) avec :
M s : Moment stabilisant.
M r: Moment de renversement.
M stabilisateur =M W /o =W i ×bi
F i : Efforts sismique appliqués à chaque étage.
d i : La hauteur de niveau i mesuré à partir de la base.
b i : Les points les plus éloigné par apport au centre de gravité pour les deux directions.
117
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
W i : Poids de niveau i.
Momentstabilisateur
Il faut vérifier que : ≥ 1 ,5
Mo m entderenversement
NIVEAU Xi Yi Mi M i× Xi M i × yi XG YG
Niveau 7 18,6 6,8894 7067,9308 131463,513 48693,8025 18,6 6,79
Niveau 6 18,5974 6,7581 7650,8505 142285,927 51705,2128
Niveau 5 18,6003 6,7736 7740,2345 143970,684 52429,2524
Niveau 4 18,6003 6,7743 7810,4419 145276,562 52910,2766
Niveau 3 18,6003 6,7749 7876,053 146496,949 53359,4715
118
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
I
Niveau 2 18,6003 6,7756 7955,6215 147976,947 53904,109
Niveau 1 18,6003 6,7763 8030,5936 149371,45 54417,7114
RDC 18,5992 6,7844 8660,5417 161079,147 58756,5791
Somme 62792,2675 1167921,18 426176,415
X G =¿11,5135464m ; Y G =¿10,0710051m
b x =L x – X G=37.5−18.6=¿18.9m
b y =L y – Y G =14−6.79=¿7.21 m
Tableau IV- 20: Vérification du renversement
Conclusion
Moment stabilisateur
La condition ≥ 1 ,5 est vérifiée dans les deux sens donc notre
Moment de renversement
structure est stable au renversement.
Avec :
N d : L’effort normal maximal.
Bc : Section du poteau.
f c28 : Résistance caractéristique du béton.
Niveau Bc (c m )
2
N (KN ) V Observation
RDC 60×60 2452,5135 0,2725015 ≤ 0.3 C.N.V
Niveau 1 55×55 2071,9762 0,27398032 ≤ 0.3 C.N.V
Niveau 2 50×50 1730,9873 0,27695797 ≤ 0.3 C.N.V
Niveau 3 45×45 1405,4979 0,27762921 ≤ 0.3 C.V
Niveau 4 40×40 1090,6199 0,27265498 ≤ 0.3 C.V
Niveau 5 40×40 788,302 0,1970755 ≤ 0.3 C.V
Niveau 6 35×35 492,6627 0,16086945 ≤ 0.3 C.V
Niveau 7 35×35 288,0402 0,09405394 ≤ 0.3 C.V
On fait les vérifications pour les poteaux les plus défavorables, et les résultats sont présentés
dans le tableau qui suit :
119
CHAPITRE IV ETUDE DYNAMIQUE ET SISMIQUE
120