Bulletin de lInstitut Scientifique, Rabat, section Sciences de la Terre, 2004, n26, 1-10.

Analyse spectrale non itrative des oscillateurs simples

sous laction des tremblements de terre
Youssef BELMOUDEN
Ecole Polytechnique Fdrale de Lausanne, Laboratoire dinformatique et de mcanique appliques la construction (IMAC), GC G1
494, ENAC-IS-IMAC, CH-1015 Lausanne, Suisse. e-mail : youssef.belmouden@epfl.ch
Rsum. Les tremblements de terre ou sismes constituent une activit gologique naturelle. Ils induisent des destructions importantes. Se
doter de moyens de prvention contre leurs effets tels que les codes parasismiques est une entreprise vitale. La nouvelle gnration de
codes parasismiques permet de dimensionner les structures et de vrifier la vulnrabilit des constructions existantes contre laction des
sismes. Lanalyse spectrale non itrative est une procdure de dimensionnement et de vrification des structures inlastiques modlises
par des oscillateurs simples sur la base de la mthode des dplacements. Cette mthode adopte lintersection entre le diagramme de
capacit obtenu par une analyse statique monotone et le spectre de rponse inlastique convertie en format acclration-dplacement. La
mthode est gnralise pour rendre compte de leffet de lcrouissage en considrant un modle damortissement bas sur la relation entre
les coefficients de rduction et la ductilit.
Mots cls: Sismes, mthode non itrative, spectre inlastique, coefficient de rduction, diagramme de capacit, oscillateur simple.
A non-iterative spectrum analysis of single degree of freedom systems under earthquake action.
Abstract. Earthquakes are a natural geological activity. The seismic activity generates massive destructions. The aseismic codes represent
an efficient tool against their consequences. The seismic design methodologies for the next generation of codes, safeguards against the
damage produced by earthquakes and permits an assessment of vulnerability of existing buildings. The capacity spectrum analysis is a
direct displacement-based seismic design procedure of inelastic structures modelled as single degree-of-freedom. This method uses the
intersection of the capacity diagram from a pushover analysis and an inelastic response spectrum converted to acceleration-displacement
format. The method is then extended to hardening behaviour with a damping model constructed based on relations between reduction
factors and ductility.
Key words: Earthquakes, non iterative method, inelastic spectrum, reduction factor, capacity diagram, single degree of freedom system.

INTRODUCTION
Dans les codes ou rglements parasismiques, ltude de la
rponse des structures des btiments aux effets des
tremblements de terre est conduite en utilisant des mthodes
simplifies bases sur les spectres de rponse. Dans le
dimensionnement ou la vrification parasismique des
structures soumises laction des sismes, une nouvelle
mthode a merg ces dernires annes et est de plus en
plus utilise. La mthode consiste utiliser une procdure
de dimensionnement directe par la mthode des
dplacements en dfinissant ds le dpart le dplacement
limite des btiments au lieu de les dimensionner sur la base
de la mthode des forces, habituellement adopte par les
codes parasismiques.
Lanalyse spectrale non itrative est une mthode de
dimensionnement des structures non linaires sous la forme
doscillateurs simples, en considrant directement un
diagramme de rponse inlastique.
Dans cet article, on propose une gnralisation de la
mthode spectrale non itrative aux systmes non linaires
avec crouissage et suramortis. Lamortissement est
considr variable, fonction de la ductilit et des
coefficients de rduction. Enfin, une procdure de
dimensionnement et de vrification simplifie des structures
de btiments est propose.

jeu pendant le sisme est value par la magnitude. Par

consquent, les mouvements de la crote terrestre gnrent
des ruptures ayant pour consquence la mise en vibration
du sol. Aprs avoir parcouru la distance ncessaire depuis le
foyer du sisme, au cours de laquelle il y a attnuation et
transformation des signaux transmis, les ondes de surface
mettent en vibration les constructions se trouvant dans
lenvironnement proche de lpicentre (quelques dizaines
de km).
Les consquences de ces effets seront prises en compte
dans les calculs (choix des spectres par exemple), puisquils
entranent des effets directs (effets de site) sous forme de
mouvements forts ( proximit de lpicentre) ou de
mouvements faibles (loin de lpicentre), et des effets
induits sous forme de grandes dformations conduisant
des liqufactions et des mouvements de terrain
(glissements, boulements, tassements et affaissements.).
Par consquent, le sisme lchelle humaine est une
source de destruction. En effet, les btiments se trouvant
dans les zones touches par un tremblement de terre
subissent, selon les amplitudes et la dure d'oscillation, le
type de frquence, mais aussi en fonction de leur mode de
construction, des dgts plus ou moins lourds, jusqu'
l'croulement.
Action sismique

Le sisme
Le sisme se traduit essentiellement par une rupture brutale
au sein de lcorce terrestre selon diffrents types de failles
et de mcanismes, et par une mission dondes lastiques
(longitudinales ou de cisaillement) dont lnergie mise en

Laction dun tremblement de terre sur une structure de

btiment (Fig. 1) se traduit par un dplacement du sol,
variable dans le temps dg(t), qui implique des composantes
de translation U(t) et de rotation (t) alternes (Fig. 1).
Laction sismique peut tre aussi dtermine laide de

Y. Belmouden Analyse spectrale des oscillateurs simples

Structure de
Btiment

Sol

Ux

Uy

PRINCIPE DE LA PARASISMIQUE

z
Uz
Figure 1. Action du tremblement de terre sur une structure de
btiment.

d g (ou a g )

t (sec)
Temps

Figure 2. Enregistrement du mouvement du sol.

lacclration du sol (ou acclrogramme) ag(t) en

translation et en rotation. Les composantes de rotation ont
gnralement des effets ngligeables. La composante de
translation verticale est plus faible que les composantes
horizontales. Le principal effet du tremblement de terre est
donc un mouvement horizontal du sol.
En sismologie, la caractrisation des tremblements de terre
seffectue de diverses manires. La Magnitude (chelle de
Richter) reprsente lnergie libre au foyer ; lIntensit
permet de caractriser les dgts observs suivant une
chelle qualitative (Mercalli). Pour les ingnieurs, les
paramtres les plus utiles sont soit les lois de dplacements
du sol dg(t) et des acclrogrammes ag(t) (Fig. 2), soit des
donnes plus globales, appeles spectres de rponse,
tablies partir de ces acclrogrammes.
ANALYSE DES STRUCTURES
SOUS LACTION SISMIQUE
Lanalyse des structures de btiment sous laction des
sismes consiste dterminer les forces maximales et
surtout les dformations maximales quelles peuvent subir.
Ces informations permettent de les dimensionner afin de
rsister aux tremblements de terre. Diffrentes approches
existent :

lapproche dynamique directe, base sur les lois

temporelles dg(t) et ag(t) ; trop coteuse, elle est rserve
aux cas de structures irrgulires et douvrages dart (art.
1.2, note (b), commentaire 1.2 du RPS2000) ;
lapproche statique (quivalente), base sur les spectres
de rponses ; moins coteuse, elle est adopte par les codes
et rglements de calcul parasismique dans le cas de
structures courantes (structures de btiments rgulires
faible et moyenne hauteur).

Concevoir une construction pour quelle rsiste aux grands

sismes sans dommages est possible, mais un tel cas ne sera
ni conomique, ni justifi cause de la faible probabilit
doccurrence de tels sismes. Il est donc plus sens de
calculer les structures en fonction des sismes les plus
importants de la rgion.
Les rglements parasismiques adoptent le principe qui fait
appel la rigidit et la rsistance des structures aux
sismes faibles, mais la ductilit (capacit subir des
dformations importantes sans risque de ruine prmature)
et les capacits dabsorption-dissipation des nergies
injectes dans la structure, sont les principaux moyens de
rsister aux sismes moyens et forts.
Le principe de la construction parasismique peut donc se
rsumer par les objectifs suivants :
confrer aux ouvrages un comportement global
satisfaisant en vue dassurer la scurit des personnes, en
ralisant une faible probabilit de ruine par effondrement
par le maintien de lintgrit structurale et de la capacit
portante rsiduelle ;
limiter les dommages conomiques en augmentant la
probabilit de rparabilit des ouvrages (augmenter la
chance de rutilisation) par la limitation des dgts
(fiabilit, non fragilit).
Lapplication des rgles parasismiques rduit donc la
vulnrabilit dune construction ou dune population de
constructions et par consquent le cot des dommages par
rapport au cot de la construction.
Le principe des codes parasismiques est dassurer une
protection nominale dont lefficacit se juge en fonction des
rsultats globaux observs sur un grand nombre de
construction. La protection nominale dun ouvrage donn
parmi un ensemble de constructions situes dans une zone
dtermine par un niveau dagression impos (carte
sismique ou microzonage), doit prsenter :
une grande probabilit pour quil napparaisse pas de
dsordre structural ncessitant rparation (sisme de faible
intensit, distance picentrale importante) (Fig. 3, d3) ;
une trs grande probabilit pour que le btiment sinistr
mais non effondr soit rendu rparable, donc rcuprable
(sisme de moyenne intensit, distance picentrale
moyenne) (Fig. 3, d2) ;
une probabilit encore plus grande de non effondrement
grave de telle faon que les chances de rparabilit soient
faibles, voire inexistantes, mais favorisant la sauvegarde

Y. Belmouden Analyse spectrale des oscillateurs simples

lments qui la constituent et de la distribution des masses

dans la structure.
La mise en uvre des spectres de rponse est base sur le
principe suivant :
Epicentre

d1

d2

d3

Foyer

Figure 3. Degr de protection des structures de btiments en

fonction de la distance picentrale.

des vies humaines (sisme de forts intensit, distance

picentrale faible) (Fig. 3, d1).
Ceci conduit adopter les principes suivants (art. 2.2 du
RPS2000) :
pour les sismes de faible intensit, on cherche limiter
les dplacements afin dviter les dommages architecturaux
dans les btiments par rigidit des structures ; il sagit du
comportement lastique linaire ;
pour les sismes moyens, les structures doivent avoir de
la rigidit ainsi que de la rsistance pour rester dans le
domaine lastique et limiter par consquent les dommages
importants ; les matriaux cassants ou fragiles sont donc
viter ;
pour les sismes de forte intensit (majeurs), les structures
doivent prsenter de la rigidit, de la rsistance, mais aussi
de la ductilit densemble, qui permettront de rsister aux
mouvements forts sans effondrement, en acceptant des
dformations irrversibles et des destructions dlments
non-structuraux dans certaines limites, afin de laisser des
chances de sauvetage et dvacuation des lieux ; il sagit du
comportement inlastique ou non linaire.
Lun des paramtres les plus influents dans les rgles
parasismiques est larbitrage conomique, qui conduit
fixer le niveau de protection en choisissant la scurit de
laction sismique contre laquelle on entend se prmunir.
Ce choix rsulte dun compromis entre des considrations
physiques (ala sismique sur lequel il ny a aucune
possibilit dagir), des considrations conomiques (rapport
du surcot rsultant de la protection de louvrage aux
bnfices que lon peut en tirer) et des considrations
politiques (intrt accord la protection de louvrage et
lintrt que lon attache sa conservation ; art. 3.2.2 du
RPS2000).
NOTION DE SPECTRE DE REPONSE
Elaboration des spectres de rponse
La mise en mouvement dg(t) du sol support de la structure
dun btiment entrane sa mise en mouvement. Pour une
structure infiniment flexible, le mouvement de chaque point
de la structure est celui du sol. En ralit, pour une structure
flexible, le mouvement de chaque point dpend de sa
rigidit globale elle-mme fonction des rigidits des

1- La structure du btiment est modlise par un oscillateur

simple quivalent. Un oscillateur simple est un systme
matriel un seul degr de libert de translation form
dune masse M et li au sol laide dune tige de rigidit K
(cf. Fig. 10). Ce systme est aussi caractris par son
amortissement C. En effet, selon les principes de la
mcanique, lorsquune structure est sollicite par une action
dynamique, telle leffet des ondes sismiques, sa raction fait
intervenir :
la distribution des masses au sein de la structure du
btiment ; de ce point de vue, leffet des masses en hauteur
est plus difficile matriser que celui des mmes masses
dans la partie basse ;
la rigidit locale et globale de la structure du btiment ;
la possibilit de dissiper de lnergie injecte par le
sisme, par frottement et interactions entre les diffrents
lments constitutifs de la structure ou entre le sol et la
structure.
Lorsquon applique ce systme un sisme sous forme
dun acclrogramme ag(t), lquation de mouvement
doscillations forces est la suivante (Principe fondamental
de la dynamique) :
M(t)+C(t)+KU(t)=Mag(t) (1)
o (t), (t) et U(t) sont respectivement lacclration, la
vitesse et le dplacement de loscillateur simple linstant t.
Loscillateur simple est caractris par sa priode propre
T=2(M/K)1/2 (sa frquence doscillation f=1/T ou sa
pulsation =2f) et par son amortissement relatif =C/C0,
avec C0 tant lamortissement critique tel que si C>C0
loscillateur retourne au repos sans oscillation
(amortissement important), si C<C0 loscillateur samortit
avec des oscillations priodiques sinusodales. Lquation
(1) peut donc tre rsolue (dtermination de (t), (t) et
U(t)) soit analytiquement (intgrale de Duhamel) soit
numriquement.
2- Ltape suivante consiste rechercher lamplitude
maximale de la rponse de loscillateur simple en
acclration, en vitesse ou en dplacement. Les rsultats
sont ensuite reports dans un graphique en fonction de la
priode T ou de la frquence f.
3- Ce calcul est ensuite appliqu de multiples oscillateurs
simples, en faisant varier leurs priodes propres T (Fig. 4).
La srie de valeurs discrtes S(Ti) obtenues reprsente le
spectre de rponse dune variable physique donne S
correspondant lacclrogramme ag(t). Lorsque la variable
S est une acclration, on note A(T) ou (A-T).
4- Dans les codes parasismiques, on tient compte des
incertitudes lies lacclrogramme rel dun site, en
considrant plusieurs acclrogrammes enregistrs dans un
site ou dans une rgion donne. Ensuite, on retient la courbe
enveloppe des spectres correspondants. Le spectre de
3

Y. Belmouden Analyse spectrale des oscillateurs simples

Rponse de loscillateur simple

Spectre de rponse

A Acclration spectrale

S (T )

Spectre de rponse lastique

T1

T2

Ti

Tn 1

Tn

T =0

Priodes

Spectre de rponse lastolasto-plastique

T =
Rduction de lordre de

Oscillateurs
simples

Sisme

2 1

a g (t ) a g (t )

a g (t )

a g (t )

HF

a g (t )

Hautes frquences

Figure 4. Construction dun spectre de rponse lastique sous

laction dun tremblement de terre.

Spectre de code:
spectre de calcul

S max

T
Tb

Tc

Priodes

Figure 5. Spectre de calcul (ou spectre de code). Tb et Tc sont les

priodes de transition.

rponse ainsi obtenu est parfois appel dans les codes

spectre de calcul (Fig. 5).
5- Par consquent, il est possible dobtenir directement,
pour une structure ayant une frquence ou une priode
propre connue, la rponse maximale de la structure sous
laction du sisme ou encore leffet maximum d un
tremblement de terre par simple lecture sur le spectre de
rponse ainsi obtenu. Ces spectres sont aussi appels
spectres de rponses lastiques. Le comportement des
oscillateurs est considr comme lastique linaire et est
dtermin pour un amortissement relatif constant et gal
5% (art. 3.4 et art. 5.2.3.3 note (a) du RPS2000).
Si loscillateur est infiniment rigide (K, T=0) :
loscillateur subit alors lacclration maximale du sol Amax
sans dformation. Si loscillateur est infiniment flexible
(K=0, T), loscillateur ne subit aucune acclration
mais subit exactement le dplacement maximal du sol Dmax
(Fig. 4) (mouvement dun corps rigide).
Si le spectre est dacclration, de vitesse ou de
dplacement, on note Smax respectivement par PSA, PSV ou
PSD. Dans lanalyse spectrale, dans le domaine lastique,
les relations entre les acclrations A, les vitesses V et les
dplacements D scrivent :
A=(2/T)V=(2/T)2D et V=(2/T)D.
4

0,5

Rduction de lordre de

BF

T (s )

Basses frquences

Figure 6. Spectre de rponse lastique et spectre de rponse

inlastique en acclration A(T).

Elaboration des spectres de rponse inlastique

Dans la nouvelle gnration de codes parasismiques, on
permet aux structures de grandes dformations et par
consquent un comportement non linaire au lieu de rester
lastique linaire. Par consquent, sous leffet des
tremblements de terre, les structures des btiments vont
subir des acclrations et des sollicitations rduites, mais en
contrepartie ils vont subir des dformations et des
dplacements plus importants. Ce concept permet une
grande dissipation de lnergie injecte par les sismes dans
les structures. Dans de telles situations, on utilise un spectre
de rponse inlastique. Ce dernier se construit partir des
spectres de rponses lastiques en appliquant un ou
plusieurs facteurs de rduction des acclrations.
Spectre de rponse inlastique de Newmark & Hall
(1982)
Lintrt du comportement non linaire dun systme un
seul degr de libert ou oscillateur simple a t montr par
Newmark & Hall (1982). Ces auteurs ont tudi un grand
nombre doscillateurs simples linaires et non linaires avec
lhypothse dun comportement lastique parfaitement
plastique. Suite une tude des spectres lastiques et
inlastiques (lastoplastiques) rsultant du sisme dEl
Centro en Californie, du 18 mai 1940 (Fig. 6), ces auteurs
ont observ que :
dans les basses frquences (BF), les systmes lastique et
inlastique ont approximativement le mme dplacement
maximal ;
dans les frquences moyennes, le principe de
conservation des nergies donne une bonne et suffisante
approximation ;
dans les hautes frquences (HF), les systmes lastique et
lastoplastique ont la mme force.
Ces rsultats peuvent tre rsums en termes des frquences
telles que :
principe dgalit des dplacements maximaux Dp=DE,
valide pour f<fc1 (oscillateurs souples) ; dans ce cas, le
facteur de rduction R (coefficient de comportement) est

Y. Belmouden Analyse spectrale des oscillateurs simples

Principe dEgalit des dplacements

Principe dEgalit des nergies

l
as
tiq

ue

A
AE

Co
m

po
rte
m

en
t

Comportement lastique

Ay = R AE

AE

Dp
Dy

Comportement
lastoplastique

Ay = R AE

D p = DE

Figure 7. Schmatisation du principe dgalit des dplacements

des oscillateurs simples ; AE acclration lastique maximale, Ay
acclration de plastification, DE dplacement lastique maximal,
Dy dplacement de plastification, Dp dplacement inlastique
maximal.

gal 1/ (Fig. 7), tant le facteur de ductilit (>1).

En effet,
Dp=DE Dy=DE Ay/k=AE/k Ay =AE
RAE=AE R=1/

principe
dgalit
des
nergies
dissipes,
Aire(OCG)=Aire(OBEF), valide pour f>fc2 (oscillateurs
assez rigides). Le facteur de rduction R est gal 1/(2
1)1/2 (Fig. 8) ; en effet, de lgalit des nergies dissipes on
a
(1/2) AE DE =(1/2)AyDy+Ay (DpDy),
do AEDE=AyDy(21) (AE)2=(Ay)2(21)
R=1/(21)1/2
principe dgalit des forces, valide pour f>33 Hz ; le
facteur de rduction R est gal lunit.
Le troisime principe nest pas toujours applicable aux
btiments. Les deux frquences de transition fc1 et fc2
dpendent de la ductilit du systme et des
caractristiques du spectre lastique. Les valeurs usuelles de
fc1 et fc2 sont respectivement 1,4 Hz et 2,0 Hz. Les
frquences intermdiaires sont obtenues par interpolation.
Il est important de diffrencier entre le facteur de
comportement (art. 3.3.4 du RPS2000) et le facteur de
ductilit. Le facteur de comportement reprsente le rapport
de leffort ou lacclration que peut supporter la structure
si celle-ci tait parfaitement lastique, leffort ou
lacclration effectivement supports par la structure relle
pour une mme dformation. Par contre, on dfinit la
ductilit (en dplacement) , appele aussi facteur de
ductilit disponible (art. 3.3.1 du RPS2000), comme tant la
capacit de la structure subir, avant la ruine, des
dformations importantes qui sont la source de
consommation dnergie leve. Par consquent, plus la
structure montre un comportement ductile, plus les
mcanismes de sa rsistance toute action sismique bass
sur la capacit dabsorber et de dissiper lnergie induite
par sa rponse globale est efficace.

Dy

DE

D
D p = Dy
F

Figure 8. Schmatisation du principe dgalit des nergies

dissipes par les oscillateurs simples.
3
Facteur d'amplification
dynamique D(T)

Dy

Dy

k
A

Dp

Comportement lastoplastique

k
o

Site 1
Site 2
Site 3

2.5
2
1.5
1
0.5
0
0

0.25

0.5

0.75

1
1.25
Priodes T (sec)

1.5

1.75

Figure 9. Facteur damplification dynamique D(T) : spectre du

code RPS2000.

Spectre de calcul du code parasismique RPS2000

Lvaluation de laction du tremblement de terre sur une
structure de btiment, selon le rglement parasismique
RPS2000, seffectue laide des paramtres suivants :
lacclration maximale du sol Amax obtenue partir du
zonage sismique, (art. 5.2.2 du RPS2000) : Amax= 0,01 g
(Zone 1), Amax= 0,08 g (Zone 2), Amax= 0,16 g (Zone 3),
avec une probabilit dapparition de 10% en 50 ans ;
un spectre de rponse en termes dacclration pour le
mouvement horizontal relatif un type de site (A-T) (art.
5.2.3 du RPS2000) ;
un spectre de rponse du mouvement vertical dduit du
spectre horizontal par un coefficient de 2/3.
Le spectre de rponse lastique (A-T) ou spectre de calcul
est dfini dans le RPS2000 pour un amortissement relatif
=5% par : A(T)=(Amax/g)D(T) o g est lacclration de la
pesanteur et D(T) le facteur damplification dynamique (art.
5.2.3.3 du RPS2000). Il reprsente lamplification des
acclrations dans la structure par rapport celle du sol.
Dautres paramtres interviennent dans le spectre de
rponse lastique pour tenir compte de leffet de site (S, art.
5.2.3.3 du RPS2000), du coefficient de priorit (I, art. 3.2.3
du RPS2000) et de lamortissement (=(5/)0,4, art. 5.2.3.3
note (d) du RPS2000). Le spectre de rponse inlastique est
obtenu par lapplication dun facteur de rduction des
5

Y. Belmouden Analyse spectrale des oscillateurs simples

Structure

Modlisation
Oscillateur simple quivalent
Acclration

K ,

Dplacement

Charge statique
monotone croissante
Dplacement
au sommet

UN

UN
Diagramme

Masse en
mouvement

Vbase

Vbase

Structure
de btiment

de capacit

rk

Ay

k
Vbase

Dy

Figure 10. Conversion du diagramme de capacit en format (A-D) pour un oscillateur simple. Dy, Ay dplacement et acclration ltat
limite lastique de la structure, r, taux dcrouissage ; k, la rigidit lastique initiale (k=Ay / Dy), rk est la rigidit plastique avec
crouissage.

acclrations, appel facteur de comportement (K, art.

3.3.4. du RPS2000) tel que A(T)=(Amax/g)D(T)/K (Fig. 9).
La priode de transition Tc (Fig. 5) est indique dans le
tableau 5.4 du RPS2000 : Tc= 0,4 s (Site 1) ; Tc= 0,6 s (Site
2) ; Tc= 1,0 s (Site 3). Cette priode dpend de la magnitude
du sisme, de son contenu frquentiel, de sa dure, et des
rapports entre les valeurs maximales des caractristiques du
mouvement (Fig. 5) : lacclration Amax, la vitesse Vmax et
le dplacement Dmax. Ces rapports varient avec la distance
entre le site et la source du sisme (commentaire 5.2.3 du
RPS2000).
METHODE DANALYSE SPECTRALE
NON ITERATIVE
La mthode danalyse spectrale non itrative est une
mthode simplifie pour lvaluation de la rponse
inlastique des structures, sans avoir recours lanalyse
dynamique directe non linaire (Chopra 1999, Fajfar 2000,
Xue 2001). Lanalyse des structures sur la base de la
mthode des forces (RPS2000), consiste tablir le spectre
de rponse inlastique ou spectre de calcul (art. 5.2.3 du
RPS2000), dterminer les forces agissant sur la structure
(art. 6.2.1.3 et 6.2.1.4 du RPS2000) et vrifier
ladmissibilit des dformations rsultantes (art. 3.1 et 8.4
du RPS2000).
Si ces dformations dpassent les limites fixes, il faut
refaire les calculs en agissant sur les paramtres considrs.
Par contre, dans lanalyse base sur la mthode des
dplacements, le dplacement limite ou cible est considr
ds le dpart. Ainsi, on travaille sur un diagramme
acclration-dplacement A(D) (ou (A-D)) au lieu de
considrer un spectre en acclrationpriode A(T). Cette
mthode est base sur lintersection du diagramme de
capacit et du diagramme de rponse inlastique sous
laction des sismes pour estimer le dplacement maximum
de la structure.
Diagramme de capacit
Le diagramme de capacit est obtenu partir de la courbe
dplacement latral en tte de la structure (UN) - effort
6

tranchant la base (Vbase) (Fig. 10) sous chargement

statique monotone croissant F. Cette courbe traduit les
performances non linaires de la structure. Elle peut tre
obtenue laide dune modlisation des structures par
lments finis en adoptant la technique dapproximation par
lments finis multicouches (Belmouden 2003, Belmouden
& Elharif 2003a, b, c). Ces lments finis sont capables de
tenir compte des caractristiques mcaniques des matriaux
constitutifs des structures (dgradations, fissurations,
crouissage) ainsi que les mcanismes susceptibles de se
dvelopper dans les structures (distorsions, glissement
dadhrence, dcollement).
Le diagramme rsultant est ensuite converti au format
acclrationdplacement (A-D) correspondant un
oscillateur simple quivalent de masse M, de rigidit K et
damortissement (Fig. 10). Enfin, ce diagramme est
idalis en forme bilinaire (Chopra 1995, Chopra & Goel
1999, 2001, Xue 2001, Chopra & Chintanapakdee 2003).
Diagramme de rponse inlastique
Le diagramme de rponse inlastique, pour une ductilit
donne , est une courbe en format (A-D) traduisant
lacclration spectrale A en fonction du dplacement
spectral D. Il est obtenu par conversion du spectre de
rponse lastique (A-T)el en format (A-D)el laide des
relations indiques dans le tableau I et dans la figure 11. Le
diagramme (A-D)el est transform ensuite en diagramme de
rponse inlastique (A-D)inel par des coefficients de
rduction appropris.
Formulation du diagramme de rponse inlastique
Considration de lcrouissage et de la variation de
lamortissement
Dans la dfinition du coefficient de comportement K
(facteur de rduction) dans le code RPS2000, on voque la
capacit de dissipation de lnergie vibratoire de la structure
transmise par les secousses sismiques en fonction du type
de la structure de btiment et du niveau de ductilit choisi
(art. 3.3.4 du RPS2000). Cette capacit est intimement lie
lcrouissage et lamortissement propre la structure.

Y. Belmouden Analyse spectrale des oscillateurs simples

Tableau I. Formules de conversion du spectre dacclration lastique en diagramme de rponse lastique

Points
A
B
C
D

Acclration
Aa=Amax
Ab=PSA
Ac=PSA
Ad=(2/Td)PSV

Priode
Ta=0
Tb
Tc=(2/PSA)PSV
Td=(2/PSV)PSD

Dplacement
Da=0
Db=(Tb/2)2PSA
Dc=(Tc/2)2PSA
Dd=PSD

Rigidit
kb=(42M/(Tb)2)
kc=(42M/(Tc)2)
kd=(42M/(Td)2)

Diagramme de rponse lastique

B

Ac

en acclration - dplacement

ED

Ap

Acclration

Ay

Aa
Ad

kb

kc

Dy

kd

Db

Dc

Es

Dd

Figure 11. Diagramme de rponse lastique en termes

dacclrationdplacement (A-D)el.

Figure 12. Amortissement quivalent dans les systmes non

linaires avec crouissage.

Ac=(2/Tc)Vc RAPSA=(2/Tc)RVPSV
Tc=(2/PSA)PSV(RV/RA)

Les coefficients sont donns pour les systmes non crouis

(r=0) (Newmark & Hall 1982) tels que :
systme avec un amortissement =5% on a
RA=1/(21)1/2 , RV=1/,
systmes suramortis avec 5% on a
RA=(3,210,68 ln())/2,12 ; RV=(2,31-0,41ln())/1,65
o les dnominateurs reprsentent les facteurs
damplification correspondant un amortissement de 5%.
Le modle damortissement pour les systmes avec
crouissage adopt pour lvaluation de est bas sur le fait
que lnergie dissipe par un systme inlastique est gale
celle dissipe par un systme lastique linaire quivalent
tel que : (,r)=0.05+eq (Chopra 1995, Chopra & Goel
2001) avec

eq=ED/(4 ES)=2(AyDpDyAp)/(ApDp)
=2(1)(1r)/[(1+rr)]
o Dp, Ap sont le dplacement et lacclration ltat
limite de performance de la structure, Dp=Dy,
Ap=Ay+r(DpDy)=Ay+r Dy(1)=Ay(r-r+1), ES est
lnergie de dformation lastique quivalent, ED est
lnergie dissipe dans le systme inlastique (Fig. 12). En
adoptant ce modle, lamortissement nest plus constant, il
est donc fonction de la ductilit (dformabilit) et de
lcrouissage (rsistance) du systme.
Elaboration du diagramme de rponse inlastique
Llaboration du diagramme de rponse inlastique (AD)inel partir du diagramme de rponse lastique (A-D)el est
base sur la formulation suivante (Xue 2001) (Fig. 13):
Point c: Ac=RAAc=RAPSA,

Dp

Dplacement

do Tc=(RV/RA)Tc,
Dc=(Tc/2)Vc=(Tc/2)RVPSV=((RV)2/RA)(Tc/2)PSV
do Dc=((RV)2/RA)Dc,
Point d:
Vd=(2/Td)Dd RVVd=(2/Td)RDDd
RVPSV=(2/Td)RDPSD
do Td=(RD/RV)Td ,
Ad=(2/Td)Vd Ad=(2/Td) (RV/RD)RVVd
Ad=((RV)2/RD)Ad,
RA, RV, RD sont les coefficients de rductions spectrales ; Ac,
Ac, Ad, Ad les acclrations spectrales aux points c, c, d,
d; Dc, Dc, Dd, Dd les dplacements aux mmes points, et
Tc, Tc, Td, Td les priodes correspondantes.
La prise en compte du comportement ductile seffectue
comme suit :
Point c* : Ac*= Ac=RAAc, Dc*=Dc=((RV)2/RA)Dc,
Tc*=(()0,5RV/RA)Tc,
Point d*: Ad*= Ad, Dd*=Dd=RDDd=RDPSD,
Td*=()0,5(RD/RV)Td,
o : Ac*, Ad* sont les acclrations spectrales aux points c*,
d*; Dc*, Dd* les dplacements aux mmes points, Tc*, Td* les
priodes correspondantes (Fig. 13).
Mthode danalyse non itrative des oscillateurs simples
sous laction des tremblements de terre
Pour un oscillateur simple modlisant une structure de
btiment sous laction sismique dfini par sa masse M et sa
7

Y. Belmouden Analyse spectrale des oscillateurs simples

Dimensionnement dune nouvelle structure contre les

sismes

Ac
Ac '

PSA
RA PSA c '

Diagramme de rponse lastique

c*

Diagramme de rponse inlastique

d
d*

d'

Dd '

Dc ' Dc*

Dd *

Figure 13. Elaboration du diagramme de rponse inlastique (AD)inel partir du diagramme de rponse lastique (A-D)el.

Diagramme de rponse lastique

Ac

Diagramme de rponse inlastique

Diagramme de capacit
bilinaire avec crouissage

c*

Ac*

AP
Ay

Vrification dune structure existante vis--vis des

secousses sismiques

rK
Branche descendante

D y DP

Figure 14. Point dintersection P dans la branche descendante du

spectre de rponse inlastique.

Diagramme de rponse lastique

Ac
P

AP
Ay

La procdure de dimensionnement des structures de

btiments contre les sismes consiste calculer leurs
caractristiques
lastiques
(initiales)
et
plus
particulirement leurs rigidits (art. 1.2 du RPS2000) :
tape 1 : on tablit un spectre de rponse lastique (PSA,
PSV et Tc) en fonction du site considr, puis conversion du
spectre en diagramme inlastique (A-D)inel ;
tape 2 : tant donn le taux dcrouissage r, on fixe le
dplacement limite admis pour la structure du btiment sous
laction des sismes Dp et sa ductilit ;
tape 3 : le dplacement lastique limite Dy est dduit
partir de Dp et ;
tape 4 : on superpose le diagramme inlastique (A-D)inel
et le diagramme de capacit bilinaire ; le point
dintersection peut se trouver soit dans la branche
descendante (Fig. 14) soit au niveau du plateau
dacclration constante du diagramme inlastique (Fig.
15) ; on dduit lacclration limite lastique Ay ;
tape 5 : on en dduit la rigidit initiale K=AyM/Dy et la
priode doscillation de la structure T=2(M/K)1/2 ; la
priode ltat limite de dformation est T=Tp[(r
r+1)/]1/2.

Diagramme de capacit
bilinaire avec crouissage

c*

Diagramme de rponse
inlastique

Cette procdure consiste vrifier les performances limites

de la structures de btiment existant sous laction des
tremblements de terre :
tape 1 : identique celle de la procdure de
dimensionnement ;
tape 2 : le diagramme de capacit tant connu (le
comportement lastique et le comportement inlastique de
la structure sont connus) : r , Ay, Dy et K ;
tape 3 : superposer le diagramme inlastique (A-D)inel et
le diagramme de capacit bilinaire (rsolution graphique) ;
en dduire le point dintersection dfinissant ltat limite de
dformation de la structure sous laction des sismes Ap, Dp
et ;
tape 4 : calculer la priode au point de performance Tp et
la frquence correspondante.

rK
APPLICATION

Dy

DP = P D y

Figure 15. Point dintersection P sur le plateau du spectre de

rponse inlastique.

rigidit K, on prsente une procdure de dimensionnement

et de vrification ltat limite de performance.
8

Pour mettre en vidence la mthode spectrale non itrative,

on tudie six structures de btiments rguliers en bton
arm modlises par des oscillateurs simples quivalents
sous laction du sisme dEl Centro de 1940 (PSA=2,71 g,
PSV=280,42 cm/s, Tc=0,66 s). On fixe pour chaque
structure le dplacement limite admis et sa ductilit (Tabl.
II).
Ltude consiste calculer les caractristiques lastiques
des oscillateurs simples quivalents aux structures. Les
caractristiques lastiques ncessaires au dimensionnement
sont donc obtenues par la mthode non itrative propose et
sont compares aux rsultats obtenus par Chopra & Goel

Y. Belmouden Analyse spectrale des oscillateurs simples

Tableau II. Caractristiques des oscillateurs simples soumis au sisme dEl Centro de 1940.
Systmes

DP [cm]
D y [cm]

22,30

22,30

19,46

44,61

44,61

44,61

3,72

5,58

9,70

7,44

11,16

22,32

0c

0,0

0,0

0,0

0,0

0,0

Tableau III. Rsultats compars de la mthode itrative et non itrative.

Mthode spectrale non itrative
Systme

RV

RA

Ay
(cm/s/s)

Mthode spectrale itrative

(s)

Tc*

Ay

(s)

(cm/s/s)

(s)

0,166

0,3015

586,71

0,5

0,893

588,1

0,49

0,25

0,3779

880,77

0,5

0,873

882,1

0,49

0,50

0,5773

1534,75

0,5

0,808

1532,3

0,49

0,166

0,3015

293,335

1,0

0,893

294,0

0,9

0,25

0,3779

440,38

1,0

0,873

441,0

0,9

0,50

0,5773

880,77

1,0

0,808

882,1

0,9

1999) en utilisant la mthode itrative (Tabl. III). Lexamen

des rsultats obtenus montre que lcart relatif (de lordre
de 2%) entre les deux rsultats est trs satisfaisant.
CONCLUSION
Dans ce travail, a t prsente une mthode danalyse
spectrale non itrative pour le dimensionnement et la
vrification directe des structures de btiments sous laction
des tremblements de terre. Elle est base sur la mthode
des dplacements, en adoptant un diagramme acclration
dplacement au lieu dun spectre acclrationpriode. La
formulation dlaboration dun tel diagramme est propose.
Dans le cadre de lamlioration des codes de construction
parasismique, on a aussi propos une amlioration pour
tenir compte de la variation de lamortissement dans une
structure sous laction dun sisme en fonction de la
ductilit et de lcrouissage. Leffet de lcrouissage a t
aussi introduit dans ltude pour pouvoir considrer
plusieurs types de structure de btiments. Une application
de la mthode compare des rsultats obtenus laide de
la mthode itrative montre la capacit de la mthode
obtenir des rsultats trs satisfaisants en vitant les calculs
itratifs. La nature volutive des codes et rglements de
calcul parasismique doit nous inciter constamment
amliorer ces rglements et suivre les progrs dans ce
domaine tout en conservant un certain degr de simplicit
dapplication.
Remerciements
Ce travail entre dans le cadre de la prparation dun projet de
recherche en collaboration avec l'Office Fdral des Eaux et de la
Gologie (OFEG) qui consistera tablir les courbes de
vulnrabilit sismique des btiments existants type de Suisse pour
le compte de la Fondation de Prvention des Etablissements

Cantonaux d'Assurance (anne 2004/2005). Mes remerciements

vont aux membres du laboratoire dinformatique et mcanique
appliques la construction (IMAC) de lEcole Polytechnique
Fdrale de Lausanne en la personne du Pr. Ian F. C. Smith, chef
du laboratoire IMAC, et du Dr. P. Lestuzzi chef du projet. Je tiens
aussi remercier vivement MM. M. Harnafi et J. Sebbani,
enseignants-chercheurs lInstitut Scientifique de Rabat pour
leurs prcieux conseils et leurs commentaires sur le prsent
manuscrit. Jexprime galement ma reconnaissance et mon amiti
Monsieur T.-E. Cherkaoui, sismologue lInstitut Scientifique
de Rabat, pour son soutien et sa disponibilit constante.

Rfrences
Belmouden Y. 2003. Analyse par lments finis multicouches des
performances non linaires des structures en bton arm vis-vis des secousses sismiques. Thse de Doctorat, Universit
Mohammed V-Agdal, Fac. Sci. Rabat, N2163, 101 p.
Belmouden Y. & Elharif A. 2003a. Modlisation semi-locale des
parois porteuses en bton arm par lments finis
multicouches. 6me Congrs de Mcanique, 15-18 avril,
Tanger, Socit Marocaine des Sciences Mcaniques, Tome 1,
pp. 164-165.
Belmouden Y. & Elharif A. 2003b. Modlisation numrique de la
tenue aux sismes des structures en bton arm. Bull. Lab.
Central Ponts & Chausses, 242, 49-65.
Belmouden Y. & Elharif A. 2003c. Modlisation des murs
porteurs en bton arm par lments finis multicouches. Rev.
Europ. Elments Finis, 12, 907-932.
Chopra A.K. & Chintanapakdee C. 2003. Inelastic deformation
ratios for design and evaluation of structures: single degree of
freedom bilinear systems. Earthq. Eng. Res. Center, Univ.
California, Berkeley, Report 2003-09, 81 p.
Chopra A.K. & Goel R.K. 1999. Capacity-demand diagram
methods for estimating seismic deformation of inelastic
structures: SDF systems. Pacific Earthq. Eng. Res. Center,
Univ. California, Berkeley, report 1999/02, 67 p.

Y. Belmouden Analyse spectrale des oscillateurs simples

Chopra A.K. & Goel R.K. 2001. A modal pushover analysis

procedure to estimate seismic demands for buildings: theory
and preliminary evaluation. Pacific Earthq. Eng. Res. Center,
College of Eng., Univ. California, Berkeley, report 2001/03,
87 p.
Chopra A.K. 1995 Dynamics of structures, theory and application
to earthquake engineering. Prentice Hall, Englewood Cliffs,
New Jersey.
Fajfar P. 2000 A non linear analysis method for performance
based seismic design. Earthq. Spectra, 16, 573-592.
Newmark N.M. & Hall W. J. 1982 Earthquake spectra and design.
Earthquake Engineering Institute.

10

RPS2000. 2001. Rglement de construction parasismique

RPS2000 (applicable aux btiments). Ministre de
lamnagement du territoire, de lurbanisme, de lhabitat et de
lenvironnement, Secrtariat dEtat lHabitat, 78 p.
Xue Q. 2001. A direct displacement-based seismic design
procedure of inelastic structures. Eng. Str. J.

Manuscrit reu le 2 janvier 2004

Version dfinitive accepte le 9 juin 2004