Etablissement de la courbe de fragilit par la mthode pushover

Hamou kada *, Ramdane kheir eddine *

*Dpartement de gnie civil U.S.T.O U.S.T.O BP1505EL MNAOUER ORAN ALGERIE. 1E-MAIL : [email protected]

RSUM.

Durant ces dernires dcennies, plusieurs sismes de diffrentes intensits ont frapp plusieurs pays et

occasionn des pertes considrables en vies humaines et dimportants dgts matriels. Ces taux de dgts ont encourag les chercheurs dans le domaine de gnie parasismique chercher des solutions afin de rduire ce risque sismique, cest pour cela plus de travaux on t ralis dans cet axe afin de comprendre rellement comment un sisme affecte une structure.La prvision ou bien lestimation du dommage est un axe du gnie parasismique qui nous donne des informations sur les endommagements qui peuvent tre subis par une structure, et les lments qui vont tre touchs lors dun futur sisme, et partir de ces prvisions.Un outil performant souvent utilis dans ce domaine est fourni par les courbes de fragilit qui dcrivent la probabilit datteindre ou de dpasser un tat de dommage provoqu par un sisme.
MOTS-CLS : analyse

pushover, courbe de capacit, spectre de rponse, comportement lastique, comportement inlastique,

vulnrabilit, sisme.

1. Introduction Une valuation fiable du comportement dun btiment face un tremblement de terre, en termes de dommages, requiert un outil permettant une analyse des structures au-del du domaine lastique. La mthode danalyse de type pushover, base sur les courbes de capacit (ou bien les courbes pushover), reprsente une nouvelle approche destimation de dommages (structurels et non structurels), incluant notamment les effets du comportement post-lastique.

7me Colloque National AFPS 2007 Ecole Centrale Paris

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

Les techniques danalyse de structures consistent essentiellement comparer un paramtre de demande de dplacement avec un paramtre de capacit. Leffort tranchant a t le paramtre utilis traditionnellement pour la conception parasismique des btiments. Lingnieur calcule la sollicitation (la force) provoque par un sisme donn (ou plusieurs) la base de la structure, et la compare avec la rsistance de ce btiment. Dans le calcul traditionnel, les forces sont rduites dune manire artificielle (c a d avec des coefficients de pondration), afin de garder la conception dans le domaine lastique. Alors que dans le domaine post - lastique, le dommage apparat progressivement, en diverses parties du btiment, provoquant la plastification de certains lments, aussi, il apparat par consquent une

redistribution des efforts, la demande dpendant ainsi du comportement de chaque composant. De cette manire, lendommagement est plus sensible, en gnral, au dplacement qu la force. Ces facteurs ont conduit au dveloppement des outils danalyse non linaire, tels que la mthode de la scante ou lanalyse time history non linaire [1].Toutefois, les deux mthodes sont relativement complexes pour un usage gnralis. Pour faciliter laccs ce type danalyse, plusieurs approches simplifies ont t dveloppes, telles que la mthode des coefficients de dplacement, du dplacement quivalent ou encore du spectre de capacit [1] [2] [3]. (Dans cette tude, nous allons dvelopper des courbes moyennes de capacit pour une construction dans le but dappliquer la mthode pushover afin destimer ltat dendommagement du btiment). 2. la Mthode pushover : 2.1 Dfinition de lanalyse pushover : Lanalyse 'pushover' est une procdure statique non-linaire dans laquelle la structure subite des charges latrales suivant un certain modle prdfini en augmentant lintensit des charges jusqu ce que les modes de ruine commencent apparatre dans la structure, [4] Les rsultats de cette analyse sont reprsents sous forme de courbe (voire figure1) qui relie leffort tranchant la base en fonction du dplacement du sommet de la structure.

Figure 1 : Signification physique de la courbe de capacit [2].

7me Colloque National AFPS 2007 Ecole Centrale Paris 3

Figure2: Niveaux dendommagement dcrits par une courbe de capacit [5]. Daprs la figure 2 on remarque que la courbe est compose de quatre segment, chaque segment correspond une tape dendommagement. a) Le premier niveau correspond au comportement lastique de la structure et reprsente le niveau de conception parasismique habituel. Il indique par consquent un tat dendommagement superficiel (ou bien de non endommagement). b) Le deuxime niveau dendommagement correspond un niveau de dommage contrl. La stabilit de la structure nest pas en danger, mais toutefois un endommagement mineur est susceptible de se dvelopper. c) Le troisime niveau reprsente un tat dendommagement avanc, sa stabilit tant en danger. Au del de ce niveau, la structure est susceptible leffondrement, ne prsentant aucune capacit de rsistance.

2.2 Origine de lanalyse pushover : Lanalyse statique pushover est base sur lhypothse que la rponse de la structure qui peut tre assimile la rponse dun systme un seul degr de libert quivalent, ce qui implique que la rponse est fondamentalement contrle par un seul mode de vibration et la forme de ce mode demeure constante durant la dure du sisme. Les chercheurs ont montr que ces hypothses donnent de bons rsultats concernant la rponse sismique (dplacement maximale) donne par le premier mode de vibration de la structure simul un systme linaire quivalent [6].

2.3 But de lanalyse pushover : Le but de lanalyse pushover est de dcrire le comportement rel de la structure et dvaluer les diffrents

paramtres en termes de sollicitations et dplacements dans les lments de la structure. Lanalyse pushover est suppose fournir des informations sur plusieurs caractristiques de la rponse qui ne peuvent tre obtenues par une simple analyse lastique [6], on cite : Lestimation des dformations dans le cas des lments qui doivent subir des dformations inlastiques afin de dissiper de lnergie communique la structure par le mouvement du sol.

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

La dtermination des sollicitations relles sur les lments fragiles, telles que les sollicitations sur les assemblages de contreventements, les sollicitations axiales sur les poteaux, les moments sur les jonctions poteau-poutre, les sollicitations de cisaillement. Les consquences de la dtrioration de la rsistance des lments sur le comportement global de la structure ce qui permet de dterminer les points forts et les points faibles de notre structure. Lidentification des zones critiques dans lesquelles les dformations est suppose tre grandes. Lidentification des discontinuits de rsistance en plan et en lvation qui entranent des variations dans les caractristiques dynamiques dans le domaine inlastique. Lestimation des dplacements inter-tage qui tiennent compte des discontinuits de la rigidit et de la rsistance qui peut tre utiliss dans le contrle de lendommagement.

2.4 Le Systme linaire quivalent : La rponse aux sismes des systmes non lastiques peut tre estime par mthode analytique approximative dans les quelle le systme non-linaire est remplac par un systme linaire dun seul degr de libert quivalent [5]. Une grande partie du travail fondamental a t accompli il y a deux dcennies par [Hudson, 1965 ; Jennings, 1968 ; Iwan et portes, 1979a)]. En gnral, les mthodes approximatives pour dterminer les paramtres du systme linaire quivalent sont classes en deux catgories : Mthodes bases sur la rponse harmonique et mthodes bases sur la rponse alatoire (Six mthodes sont disponibles dans la premire catgorie, et trois dans la deuxime catgorie). La priode normale de vibration et la constante d'amortissement sont disponibles pour chaque mthode (Iwan et portes, 1979a).D'une manire gnrale, les mthodes bases sur la surestimation harmonique de rponse prennent considrablement en compte le dcalage de priode, tandis que les mthodes considrant la rponse alatoire donnent beaucoup plus dimportance lvaluations de la priode (Iwan et portes, 1979b).

La mthode scante de rigidit (Jennings, 1968) est utilis dans la mthode de spectre de rponse pour vrifier l'adquation d'une conception structurale ( Freeman, 1978 ; Deierlein et Hsieh, 1990 ; Reinhorn et autres, 1995) et a t adapt pour dvelopper "le procd statique non-linaire" dans le rapport ATC-40 [1] et le rapport FEMA-274 [6]. Bas sur la rponse harmonique, ces deux mthodes sont connues pour ne pas tre assez prcise comme mthodes bases sur la rponse alatoire (Iwan et portes, 1979a). Plusieurs recherches rcentes ont tais faites par plusieurs chercheurs, (Saiidi et Sozen (1981), Fajfar et Fischinger (1988), Qi et Moehle (1991), Miranda (1991), Lawson et al. (1994)) ont montr que le

7me Colloque National AFPS 2007 Ecole Centrale Paris 5

dplacement dun systme plusieurs degrs de libert peut tre approximativement gal celui dun systme un seul degr de libert quand la rponse dominante est celle du premier mode [6]. 2.4.1 Hypothse : Lhypothse de base commune entre toutes les approches, est que la dforme du systme plusieurs degrs de libert peut tre reprsente par une forme de vecteur {! } qui reste constante durant le mouvement. 2.4.2 Mise en quation : Le vecteur de dplacement relatif x sous la forme : x= {! }.xt (xt : Dplacement au sommet de la structure). Lquation diffrentielle du mouvement du systme sera donc crite comme suit : (1)

[ M ]{ } &&t + [C ]{ } xt + {f }= ![ M ]{} &&g " x " & 1 x

[M] : matrice de masse. [C] : matrice damortissement. {f} : vecteur des forces dans les niveaux de la structure.

(2)

&& g : Lacclration du sol. x

x* sera not comme tant le dplacement du systme un seul degr de libert quivalent, lexpression suivante :
T { } [ M ]{ } ! ! = x T { } [ M ]{} t ! 1

donn par

"

(3)

En multipliant lq.(1) par {! }T, et en utilisant lq.(2), on obtient lquation diffrentielle du systme quivalent :

& M * &&* + C * x * + f * = ! M * &&g x x

(4)

M*, C*, p* dsignent les proprits du systme quivalent et elles sont donnes par les formules suivantes avec {! n = 1 , n : niveau du toit} :

M * = { } [ M ]{} ! 1

(5)

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

{ }T [ M ]{} ! 1 = { } [C ]{ } T ! ! { } [ M ]{ } ! !
T

(6)

f * = { } {f } !

(7)

La forme du vecteur {! } est connue et est suppose reprsenter la dform du systme plusieurs degrs de libert durant sa rponse au sisme dans le premier mode. Plusieurs formulations localisent le nud pour lequel le dplacement est calcul au sommet du btiment et la forme du vecteur {! } est normalise pour avoir une amplitude gale 1 au sommet [5]. On dduit alors leffort tranchant Vy et le dplacement correspondant xt,e , comme le montr par la figure (3) ci-dessous :

Figure 3 : courbe pushover. De la courbe bi-linaire reliant ces deux paramtres, qui peuvent tre calculs avec Lq. (2) et lq. (4) pour calculer les paramtres du systme un seul degr de libert et qui sont donnes par les expressions ci aprs :

" y

{ }T [M ]{ } ! ! = x { }T [M ]{} t ,e ! 1
T

(8)

f y* = { } { y } ! f
{fy} : vecteur des forces dans chaque tage la limite lastique (Vy = {1}T {fy}).

(9)

Ayant calcul ces caractristiques, on peut dterminer la priode lastique du systme quivalent comme suit :

Teq = 2!

x* M * y f Y*

(10)

7me Colloque National AFPS 2007 Ecole Centrale Paris 7

3. Application de la mthode pushover sur une structure en bton arm : 3.1-Introduction : La mthode pushover qui a t dcrite, qui se base sur laugmentation de la charge jusqu' ce que les modes de ruine apparaissent. Dans cette partie nous allons essayer de vrifier la validit de cette mthode qui sera illustre par un exemple dune structure en bton arm et de dfinir son point de performance ainsi le positionnement des rotules plastiques.

Ainsi nous avons pris trois types de sisme diffrent (sisme de faible, moyen et forte magnitude), pour essayer de comparer leurs manires dinfluencer sur le comportement de notre structure. Avant de commenc lapplication, il faut not certains hypothses qui ont tais pris en compte dans la modlisation :
-Le plancher est considr comme tant un diaphragme rigide. -Les poteaux sont considrs qui travaille a la compression avec une rigidit de flexion pris gal 0,7.E c .I g . -Les poutres sont considres qui travaille la traction avec une rigidit de initial 0,5.E c .I g . -L'interaction de sol structure n'a pas t considre dans le calcul.

3.2-Description de la structure : Il sagit dune structure R+5 en bton arm, dalles en corps creux de 20 cm dpaisseur,et les dimensions des poutres principales sont de (30x60) cm, et les poutres secondaires sont de dimensions de (30x40) cm, ainsi les poteaux sont de section (30x30) cm, la hauteur de chaque niveau est de 3m.

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

Figure 4 : caractristique verticaux de la structure.

Figure 5 : caractristique horizontaux de la structure.

7me Colloque National AFPS 2007 Ecole Centrale Paris 9

Daprs la figure (6), on peut visualiser notre structure en une vue en trois dimensions (3D).

Figure5 : vue en 3D de notre structure. Proprit des matriaux : a) le bton : - contrainte de compression.. ! b = 25 MPa. -Module dlasticit...E =32164 MPa. b) Lacier : - la limite lastique ..fy = 413 MPa. -Module dlasticit. Es=248211.28 MPa. 3.3-Lanalyse Pushover : 3.3-1.Distribution des forces latrales : La distribution des forces est prise celle qui a la forme triangulaire donne par lquation :

F px =

1 C1 .C 2 .C 3

!F
i=x

wx

!w
i=x

(11)
i

&63 0 0 0 0 # $ 0 63 0 0 0 ! $ ! M= $0 0 63 0 0 ! [t] $ ! $ 0 0 0 63 0 ! $ 0 0 0 0 60! % "

Figure 6: matrice de masse. 3.3-2.Dtermination de la courbe pushover : Nous avons utilis pour lobtention de la courbe pushover le logiciel SAP2000 non-linaire version8.08.Ce logiciel t dvelopper par le CSI (Computers and Structures, Incorporation) [9] [10]. En premier lieu les rsultats sont donns sous forme de courbe non linaire qui relie leffort tranchant la base et le dplacement au sommet de la structure voire la figure (6).

10

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

Nous allons dterminer cette courbe pour voire quelle est le sens le plus vulnrable. a) Lanalyse pushover Dans le sens x-x : En dtermine alors la courbe de fragilit, est elle est donne sous forme :

Figure6 : dplacement en fonction de leffort tranchant la base. Dans la figure (7), on la forme gnrale de notre courbe, on ajoutant le spectre de rponse.

Figure7 : courbe de fragilit selon le logiciel suivant le sens x-x.

7me Colloque National AFPS 2007 Ecole Centrale Paris 11

En plus ce logiciel nous permet de connatre la possibilit de la formation des rotules plastiques et leurs positionnements dans la structure, et leurs positionnements au niveau de la structure tudier, ce que lon nous montre au dessous la figure (8) :

Figure 8 : Formations des rotules plastiques. Commentaire : Daprs lanalyse pushover nous pouvons conclure que pour un effort tranchant dune valeur de 199,102 KN qui correspond un dplacement de 45cm, provoque lapparition des rotules plastique de ruine. Pour cette valeur de leffort tranchant qui peut caus la ruine de la structure on peut dire que notre structure est fragile.Pour la conserv ou amlior son niveau de performance il faut intervenir il faut donc utilis ou remdier par lutilisons des systmes de renforcement existantes.

b) Lanalyse pushover Dans le sens y-y : La figure (9), nous montre les rsultats trouvs suivant le sens y - y :

12

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

Figure9 : Courbe de fragilit suivant le sens y-y.

Figure 10 : Formations et des rotules plastiques.

Commentaire : Daprs lanalyse suivant le sens y-y, on remarque quil y a apparition de rotule du typa LS au niveau des poutre qui ne causeront pas un grand dommage. Conclusion : On comparons les deux courbe de vulnrabilit calculer dans les deux sens (x-x) Et le sens (y-y), on constat bien que la structure est plus vulnrable dans le sens (x-x). Donc pour augmenter la rsistance de la structure, il est prjudiciaux que limplantation des systmes de renforcement soit dune manire a augment la rsistance dans le sens vulnrable c d suivant le sens (x-x). 4. Etude de lInfluence des sismes sur les courbes de capacits : Dans cette partie nous avons pris trois type de sisme, pour essayer de comparer leurs influence sur le comportement de notre structure, ainsi de voire ou prdire le point de performance de la structure pour chaque sisme.

7me Colloque National AFPS 2007 Ecole Centrale Paris 13

Nous avons essay de prendre diffrents sismes : 1)-Sisme de faible magnitude (LOW earthquake). 2)- Sisme de moyenne magnitude (Moderate earthquake). 3)- Sisme de forte magnitude (HIGH earthquake). 4.1) Structure ltat initial : a) seime de faible magnitude (LOW): On incrmentant les donnes suivantes : CA=0,126 et CV=0,126. La figure (11) nous montre la courbe de Pushover obtenus pour ces donnes.

Figure 11 : dtermination du point de performance de la structure sous sisme de

faible magnitude.

b) sisme de moyenne magnitude (Moderate): On incrmentant les donnes suivantes : CA=0,2.et CV=0,2. La figure (12) nous montre la courbe de Pushover obtenus pour ces donnes.

14

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

Figure 12 : dtermination du point de performance de la structure sous sisme de magnitude modr. c) sisme de forte magnitude (High) : On incrmentant les donnes suivantes : CA=0,308 et CV=0,518. La figure (13) nous montre la courbe de Pushover obtenus pour ces donnes.

Figure 13 : dtermination du point de performance de la structure sous sisme de forte magnitude.

7. Conclusion : Les rsultats quon peut conclure peuvent tres reprsent par la figure (14) suivant :

7me Colloque National AFPS 2007 Ecole Centrale Paris 15

Figure1 : reprsentation des rsultats.

Le degr de fragilit augmente, lorsque la sismicit dun sisme augmente aussi.

8. Bibliographie [1] ATC, ATC 40, Seismic evaluation and retrofit of concrete buildings , Applied Technology Council, 1996.

16

7me Colloque National AFPS 2007 Ecole Centrale Paris - 013

[2] Chopra, A.K. and Goel, R.K. (2001). A Modal Pushover Analysis Procedure to Estimate Seismic Demands for
Buildings: Theory and Preliminary Evaluation.Tech. Rep.2001/3, Pacific Earthquake Engineering Research Center, University of California, Berkeley, CA.

[3] P.Fajfar. A nonlinear analysis method for performance based seismic design. Earthquake spectra, Vol.16, N.3, pp.573592, 2000.

[4] H.Krawinkler, Procedure and construction of pushover analysis of seismic performance evaluation engineering structure, vol: 20, edition: Elsever science, Department of civil engineering standford university U.S.A 1998. [5] M.Inel, M.Erich, B.Edgar, F.Black, A. Mark. USEE 2001.Utility Software for Earthquake Engineering. Report and Users Manual.Civil and Environmental Engineering.University of Illinois at Urbana.. [6] FEMA273 Federal Emergency Management Agency. NEHRP recommended Provisions for Seismic Regulations for New Buildings and Other Structures. Chapitre3, Washington, D.C. [7] P.Fajfar. Capacity spectrum method based on inelastic demand spectra. Earthquake engineering and structural dynamics
.1999.

[8] J.A.Castillo Hiyane, Study on Criteria for Seismic Evaluation of Existing Peruvian School RC building, Earthquake Engineering Course, faculty of Civil engineering, Lima-Peru, July 4th, 2003. [9] SAP2000, web tutorial1-quick pushover analysis tutorial (1999).computer and structures, Inc.Berkeley, California. [10] SAP2000, web tutorial2-quick pushover analysis tutorial (1999).computer and structures, Inc.Berkeley, California.