Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
1/20
Organisme(s) : EDF/EP/AMV
Manuel de Référence
Fascicule R4.05 : Analyse sismique
Document : R4.05.01

Réponse sismique par analyse transitoire
Résumé
Les méthodes les plus fréquemment utilisées pour l'analyse sismique des structures sont les méthodes
spectrales et les méthodes transitoires.
Les méthodes transitoires (directe linéaire ou non, par synthèse modale) permettent de calculer la réponse de
structures sous l'effet de séismes imposés : excitation unique (identique en chacun des points d'ancrage de la
structure) ou multiple et de prendre en compte leur éventuel comportement non linéaire.
En ce qui concerne les méthodes spectrales, on calcule la réponse maximale, pour chaque mode de vibration,
de chaque point d'ancrage. La réponse maximale de l'ensemble de la structure est alors déterminée par
combinaison des réponses maximales des modes. Ce type d'analyse est explicité dans la documentation de
référence [R4.05.03].
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
2/20
Table des matières
1 Comportement sismique d'une structure ............................................................................................... 3
1.1 Définitions........................................................................................................................................ 3
2 Réponse sismique d'un système à un degré de liberté ......................................................................... 3
3 Réponse sismique d'un système à plusieurs degrés de liberté ............................................................. 5
3.1 Equations du mouvement dans le repère absolu ............................................................................ 5
3.2 Equations du mouvement dans le repère relatif.............................................................................. 6
3.2.1 Décomposition du mouvement absolu ................................................................................... 6
3.2.2 Excitation simple ou multiple.................................................................................................. 7
3.2.3 Modélisation de l'amortissement ............................................................................................ 9
3.2.4 Equation fondamentale de la dynamique ............................................................................... 9
3.3 Calcul du chargement sismique ...................................................................................................... 9
3.4 Chargement de type onde incidente.............................................................................................. 10
4 Réponse sismique transitoire par synthèse modale ............................................................................ 11
4.1 Description de la méthode............................................................................................................. 11
4.2 Choix de la base modale............................................................................................................... 12
4.3 Calcul de la réponse dynamique de la structure étudiée par synthèse modale ............................ 12
4.4 Prise en compte des modes négligés par correction statique....................................................... 13
4.5 Prise en compte du caractère multi - supporté d'une structure..................................................... 14
4.6 Post traitements ............................................................................................................................ 14
5 Réponse sismique transitoire directe................................................................................................... 15
5.1 Prise en compte d'un amortissement équivalent à l'amortissement modal .................................. 15
5.2 Prise en compte d'une sollicitation multi appuis avec restitutions des champs relatifs et absolus15
6 Interaction sol-structure ....................................................................................................................... 16
6.1 Impédance d'une fondation ........................................................................................................... 16
6.2 Prise en compte d'un amortissement modal calculé selon la règle du RCC-G............................. 17
6.3 Répartition des raideurs et amortissement de sol......................................................................... 18
6.4 Prise en compte d'une frontière absorbante ................................................................................. 19
7 Bibliographie ........................................................................................................................................ 20
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
3/20
1
Comportement sismique d'une structure
1.1 Définitions
L'analyse du comportement sismique d'une structure consiste à étudier sa réponse à un mouvement
imposé : une accélération, en ses différents appuis. L'accélération imposée est un signal temporel
(t) appelé accélérogramme (cf. [Figure 1.1-a]).
Figure 1.1-a : Accélérogramme LBNS
Le mouvement sismique considéré dans le calcul est soit un accélérogramme réel connu et lu par
l'opérateur LIRE_FONCTION [U4.21.08] soit un accélérogramme synthétique calculé directement dans
le code, par exemple avec la procédure FORMULE [U4.21.11].
2
Réponse sismique d'un système à un degré de liberté
Soit un oscillateur simple constitué d'une masse m reliée à un point fixe par un ressort k et un
amortisseur c pouvant se déplacer dans une seule direction x (cf. [Figure 2-a]). Cet oscillateur à un
degré de liberté est soumis à un accélérogramme (t) horizontal en son support (point A).
repère galiléen
support
k
m
xe
A
c
xr
xa
Figure 2-a : oscillateur simple soumis à une sollicitation sismique
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
4/20
Les déplacements de l'oscillateur sont mesurés ou calculés, soit dans un repère relatif lié au point A :
déplacement relatif xr , soit dans un repère absolu (Ra) : déplacement absolu xa . Le déplacement
absolu xa se décompose en un déplacement d'entraînement uniforme en translation xe et en un
déplacement relatif xr :
x = x + x
a
r
e
éq 2-1
On en déduit par dérivation la relation entre les accélérations :
!
x =
+
avec (t) = x
a
!
x
(t)
r
! e!
éq 2-2
La masse est soumise à une force horizontale de rappel qui est proportionnelle au déplacement relatif :
F = -k. x
r
r et à une force horizontale d'amortissement supposée proportionnelle à la vitesse relative :
F = -c. x
v
!r .
L'équation du mouvement de la masse s'écrit alors : - k. x - c. x! = m.x
r
r
! a
! .
Soit, compte tenu des équations [éq 2-1] et [éq 2-2] :
m. x! + c.x! + k.x = -m. (t) = p(t
r
r
r
)
éq 2-3
Remarque :
L'étude de la réponse sismique d'un oscillateur à un degré de liberté dans le repère relatif consiste
donc en l'étude de la réponse d'un oscillateur à une force
(
p t) de forme quelconque. La solution
de l'équation de mouvement [éq 2-3] est alors fournie par l'intégrale de Duhamel :
1
t
x =
p( ).e-..(t-) .sin
.
.

[ D(t - )]d
r
m D 0
avec :
p(t) = -m. (t)
k
c
=
, =
= . 1- 2
et
m
.
2 m
D
.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
5/20
3
Réponse sismique d'un système à plusieurs degrés de
liberté

3.1
Equations du mouvement dans le repère absolu
L'équilibre d'un système mécanique consiste à écrire, quelque soit l'instant de calcul considéré, que la
somme des forces internes, d'inerties et d'amortissement est égale aux forces extérieures imposées
sur ce dit système : F
+ F
+ F
= F
iner
amo
int
ext .
Dans le cas d'un comportement linéaire, su le système est représenté par un modèle d'éléments finis
ou d'éléments discrets, on a (après discrétisation) :
F
iner = M X!a
F
int = K Xa
·
Xa est le vecteur des déplacements nodaux de la structure discrétisée, dans le repère
absolu ;
·
M est la matrice masse de la structure ;
·
K est la matrice raideur de la structure ;
·
F
= F - F
ext
e
c est le vecteur des forces imposées sur la structure étudiée, Fc celui des
éventuelles forces de choc (cf. [R5.06.03]).
Pour simplifier la présentation, on considère que la structure est uniquement sollicitée par des
déplacements imposés au niveau de ses différents supports. Ainsi, Fe = 0 .
Dans le but de simplifier la présentation, on sépare généralement les degrés de liberté en deux, en
fonction de leur type :
· les degrés de liberté de structure non soumis à un mouvement imposé - également appelés
degrés de liberté actifs - ce sont les inconnues du problème ;
· les degrés de liberté de structure soumis à un mouvement imposé - également appelés
ddl_impo - ce sont les conditions aux limites en déplacement du problème (conditions limites
de Dirichlet).
Sur les bords de la structure où les déplacements Xs sont imposés, on a : B X = X
a
s . B est la
matrice de passage de tous les degrés de liberté de la structure aux degrés de liberté de structure
soumis à un mouvement imposé.
L'équilibre du système s'écrit alors, quelque soit v appartenant à l'espace des déplacements
cinématiquement admissible c'est à dire, quelque soit v tel que B v = 0 :
M X!a + Famo + K Xa - Fext,v = 0
B X

a = Xs
Soit :
M
X!
T
a + Famo + K Xa - Fext = -B .

éq 3.1-1
B X

a = Xs
F
B
a = - T . est le vecteur des forces de réactions exercées par les appuis sur la structure.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
6/20
En tenant compte de la partition des degrés de liberté, le vecteur des déplacements dans le repère
x
a
m
m xs
absolu s'écrit : Xa = . Les opérateurs décrivant la structure deviennent : M =
,
x


s
m
m
sx
ss
k
k xs
K =
avec m
mT
=
et k
kT
=
et le vecteur des forces extérieures appliquées à
k
k

sx
xs
sx
xs
sx
ss
- fc
la structure s'écrit : Fext =
.
0
L'équation fondamentale de la dynamique dans le référentiel absolu s'écrit alors, en tenant compte de
la partition des degrés de liberté :
m

m x
!
k
k x
- f
xs
a

. + F
xs
a
c
+
. =

m

m
x
amo

!
k
k
x
f
sx
ss s
sx
ss s
a
Soit, en ne considérant que les degrés de liberté actifs :
m x! + F
+ k x = -f - m x
a
amo
a
c
xs !
- k x
s
xs
s
Cette approche nécessite la connaissance des déplacements et vitesses absolus associés à
l'accélérogramme (t) or les enregistreurs mesurent soit des accélérations soit des vitesses. On peut
remonter aux déplacements par simple ou double intégration avec la commande CALC_FONCTION
[U4.62.04]. Cependant, les incertitudes de mesure donnent des dérives qu'il convient de corriger : les
déplacements sont donc moins bien connus que les vitesses et accélérations. On gardera en mémoire
les ordres de grandeur des amplitudes maximales suivants :
· quelques dixièmes de « g » pour les accélérations ;
· quelques dizaines de cm/s pour les vitesses ;
· quelques dizaines de cm pour les déplacements.
On s'assurera également qu'en fin de séisme les vitesse et déplacement sont réalistes c'est-à-dire au
plus de quelques dizaines de cm pour le déplacement, nulle pour la vitesse.
3.2
Equations du mouvement dans le repère relatif
3.2.1 Décomposition du mouvement absolu
Les sollicitations subies par une structure lors d'un séisme sont classées en deux types dans les règles
de construction (ASME, RCC-M) :
· les contraintes induites par le mouvement relatif de la structure par rapport à sa déformée
statique ou contraintes primaires. Ces sollicitations sont dues aux effets inertiel du séisme ;
· les contraintes induites par les déplacements différentiels des ancrages ou contraintes
secondaires.
Le plus souvent, on décompose donc l'étude des structures en l'étude de la déformée statique due aux
mouvements des supports (c'est le mouvement d'entraînement) et en l'étude des vibrations induites par
les accélérations des supports autour de cette déformée (c'est le mouvement relatif).
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
7/20
X

Xa
s déplacement imposé des supports
Xe
mouvement absolu Xa
Xr
mouvement d'entrainement Xe
mouvement relatif Xr
(Ra)
Le déplacement absolu de tout point M de la structure, non soumis à un déplacement imposé, est égal
à la somme du déplacement relatif et du déplacement d'entraînement de ce point :
X (M) = X (M) + X (M
a
r
e
)
éq 3.2.1-1
Soit :
·
Xa, le vecteur des déplacements dans le référentiel absolu ;
·
Xr , le vecteur des déplacements relatifs définis comme le vecteur des déplacements de la
structure par rapport à la déformée qu'elle aurait sous l'action statique des déplacements
imposés au niveau des supports. Xr est donc nul aux points d'ancrage : B Xr = 0 ;
·
Xe , le vecteur des déplacements d'entraînement définis comme les déplacements de la
structure sollicitée statiquement par les déplacements imposés des supports
B
X = X
e
s

= r +
X
.X .
K X = -B

T
e
e =
s
e
. avec

e
est la matrice des modes statiques. Les modes statiques représentent, en l'absence de
forces extérieures, la réponse de la structure à un déplacement unitaire imposé à chaque
degré de liberté de liaison (les autres étant bloqués).
3.2.2 Excitation simple ou multiple
Pour expliciter plus en détail l'approche en mouvement relatif, et plus particulièrement le calcul des
composantes d'entraînement, il est nécessaire d'introduire la notion d'excitation simple ou multiple.
3.2.2.1 Excitation
simple
On considère que le mouvement sismique imposé est un mouvement de corps solide. On dit
généralement que la structure est mono - supportée.
Xa(t)
Xr(t)
Xs(t)
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
8/20
Le déplacement absolu de tout point M de la structure, non soumis à un déplacement imposé se
décompose donc en un déplacement relatif par rapport à un repère mobile lié au support où est
imposé le mouvement sismique et en un déplacement d'entraînement rigide.
Dans ce cas, les modes statiques correspondent aux six modes de corps rigide. Comme la structure
est linéaire élastique, on étudie séparément les effets des six composantes du mouvement sismique.
Pour chaque direction sismique, on écrit les forces d'inerties induites par le séisme simplement sous la
forme suivante :
P(t) = -M . X
! s = - (t).M
·
(t) est l'accélérogramme du mouvement sismique dans une direction ;
·
est le mode de corps solide et unitaire dans cette direction ;
· Les sismographes ne mesurent que des signaux de translation. Considérer que la structure
étudiée est mono supportée revient à supposer que tous ses points d'appuis subissent la
même translation. Dans ce cas, les composantes de [] valent 1 pour les degrés de liberté
qui correspondent à des déplacement dans la direction sismique considérée et 0 pour les
degrés de liberté qui correspondent à des déplacement dans des directions sismiques
perpendiculaires à celle considérée ou à des rotations.

Cependant, vu la taille des modèles, l'analyse sismique complète d'un équipement s'effectue
généralement en plusieurs étapes. L'analyse sismique détaillée de l'équipement considéré
utilise alors comme excitations, les accélérations calculées lors de la première étape. Elles se
composent des six accélérogrammes de translation et de rotation. On calcule donc les trois
modes correspondant aux déplacements imposés de translation et les trois modes
correspondant aux déplacements imposés de rotation. Si le mouvement sismique est une
"
"
rotation imposée, en un point M, M = OM pour les degrés de liberté qui
correspondent à des déplacement de translation et pour les degrés de liberté qui
correspondent à rotations.
3.2.2.2 Excitation
multiple
On ne peut pas toujours considérer que :
· les accélérations subies par l'ensemble des points d'ancrage de la structure étudiée sont
identiques et en phase ;
· les supports sont indéformables et animés d'un même mouvement de corps rigide.

Dans ce cas, on dit que la structure est multi - supportée. Les modes statiques =
Id
correspondent alors aux 6.nb_supports modes statiques (ou 3.nb_supports modes) où nb_supports est
le nombre d'accélérogrammes différents subis simultanément par la structure. Ils sont calculés par
l'opérateur MODE_STATIQUE [U4.52.04] avec l'option DDL_IMPO. Ils sont solution de l'équation
suivante :
X = X
e
s
k

k 0

xs
soit
. =

éq 3.2.2.2-1
K X = -B



T
e
. e
k
k
Id
f
sx
ss

a
e
Soit, en ne considérant que les degrés de liberté actifs : k. + k .Id = 0 .
xs
Les forces d'inerties induites par le séisme s'écrivent alors simplement :
nb_supports
P(t)= -
M. .X
m ! s (t )
m
m=1
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
9/20
3.2.3 Modélisation de l'amortissement
On considère que l'amortissement dissipé par la structure est de type visqueux c'est-à-dire que la force
d'amortissement est proportionnelle à la vitesse relative de la structure :
F
= C X
amo
! r où C est la matrice d'amortissement de la structure.
Cela revient à négliger l'effet de la vitesse imposée. En effet, on peut écrire plus généralement :
F
= C X! = C X! + C .X
amo
a
r
! s .
Dans le cas d'une excitation uniforme à la base (cas du mono-appui), l'amortissement n'intervient que
sur les déplacements relatifs (les forces d'amortissement sont nulles pour un déplacement rigide).
Dans le cas d'une excitation multiple où la solution statique n'est plus un déplacement rigide,
considérer que la force d'amortissement est proportionnelle à la vitesse relative de la structure est une
hypothèse simplificatrice.
3.2.4 Equation fondamentale de la dynamique
L'équation fondamentale de la dynamique [éq 3.1-1], dans le repère relatif, s'écrit alors, compte tenu
des équations [éq 3.2.1-1] et [éq 3.2.2.2-1] :
M X
! + C X! + K X = -M..X! + F - BT .
r
r
r
s
ext
r
éq 3.2.4-1
Soit, en partitionnant les degrés de liberté :
m

m x
! c
c x
! k

k x
- f
.
.

xs
r
xs
r
xs
r
c
(m + m I
xs
)dx!s
. +
. +
. =
-

m

m

0
c
c
0
k

k
0
f
sx
ss
sx
ss
sx
ss

a
r
(m . + m .I

sx
ss
)dx!s
avec c
= cT
sx
xs
Soit, en ne considérant que les degrés de liberté actifs :
m.x! + c.x
r
! + k.x = -f
r
r
c - (m. + m .I
xs
)dx!s
Les principaux avantages de l'approche en déplacement relatif par rapport à celle en déplacement
absolu sont les suivants :
· il n'est pas nécessaire d'intégrer l'accélérogramme (t) ;
· les déplacements relatifs obtenus permettent de déterminer directement les contraintes primaires
induites par le séisme.
3.3
Calcul du chargement sismique
Le chargement sismique (cf. [§3.2]) - M. soit - (m. + m .I
xs
)dx!s sur les degrés de liberté
actifs est construit par l'opérateur CALC_CHAR_SEISME [U4.43.01]. Il est utilisable directement lors
d'une analyse transitoire directe avec DYNA_LINE_TRAN [U4.54.01] ou d'une analyse transitoire par
synthèse modale avec DYNA_TRAN_MODAL [U4.54.03]. En revanche, lors d'une analyse transitoire
directe non linéaire avec DYNA_NON_LINE [U4.32.02], il faut le transformer en un concept de type
charge. Ceci est réalisé à partir de l'opérateur AFFE_CHAR_MECA [U4.25.01] de la façon suivante :
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
10/20
char_sei = CALC_CHAR_SEISME (...) ;
charge = AFFE_CHAR_MECA (MODELE :..., VECT_ASSE : char_sei) ;
dyna_nlin = DYNA_NON_LINE (
excit(CHARGE: con_lim, CHARGE: cham_no, FONC_MULT: acceler)
...) ;
Dans le cas d'une structure mono supportée, il suffit d'indiquer la direction du séisme :
mono_x = CALC_CHAR_SEISME (MATR_ASSE : masse,
DIRECTION(...), MONO_APPUI :'OUI') ;
Dans le cas d'une structure multi supportée, il faut au préalable avoir calculé les modes statiques. On
calcule autant de chargements sismiques que de supports qui subissent une accélération différente.
multi_xi = CALC_CHAR_SEISME (MATR_ASSE : masse, DIRECTION(...),
NOEUD : NOI, MODE_STAT : mode_stat,) ;
3.4
Chargement de type onde incidente
Il est également possible d'imposer un chargement sismique par onde plane par l'intermédiaire de la
commande AFFE_CHAR_MECA et le mot clé facteur ONDE_PLANE. Cela correspond aux chargements
classiquement rencontrés lors des calculs d'interaction sol-structure par les équations intégrales.
En harmonique, une onde plane élastique est caractérisée par sa direction, sa pulsation et son type
(onde P pour les ondes de compression, ondes SV ou SH pour les ondes de cisaillement). En
transitoire, la donnée de la pulsation, correspondant à une onde stationnaire en temps, doit être
remplacée par la donnée d'un profil de déplacement dont on va prendre en compte la propagation au
cours du temps dans la direction de l'onde.
Plus précisément, on caractérise :
· une onde P par la fonction u(x,t) = f (k.x - C t)k
p
· une onde S par la fonction u(x,t) = f (k.x - C t) k
s
Avec :
·
k , vecteur unitaire de direction
·
f représente alors le profil de l'onde donné selon la direction k .
O
Front d'onde « principal »
k
correspondant à l'origine
du profil
H
Fonction f
H0 est la distance du front d'onde principal à l'origine O, portée par le vecteur directeur de l'onde à
l'instant initial du calcul, H la distance du front d'onde principal à l'origine O, à un instant quelconque.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
11/20
Remarque :
Ce type de charge est disponible dans un calcul transitoire direct linéaire DYNA_LINE_TRAN ou
non DYNA_NON_LINE.
L'utilisation de ce type de chargement sera détaillé dans une note spécifique.

4
Réponse sismique transitoire par synthèse modale
4.1
Description de la méthode
La méthode de recombinaison modale consiste à décomposer le mouvement relatif de la structure sur
la base des modes propres. Comme celui-ci est nul au niveau des supports, on projette l'équation de la
dynamique sur la base des modes propres bloqués (modes propres obtenus en bloquant tous les
degrés de liberté de liaison).
X
Q
r = .
·
est la matrice des modes propres bloqués ;
·
Q le vecteur des inconnues généralisées sur la base des modes propres bloqués.
Les modes propres bloqués sont solution de :
(
0
K -- 2i.M) i = où F
F

i sont les réactions modales aux points d'appuis.
i
L'équation du mouvement projetée sur la base des modes dynamiques s'écrit alors :
M Q
! (t) + C ! (
Q t) + K Q(t)
T
= - .M..X
T
! s + .F
T
ext - .B T
G
G
G
.r
M G , CG et K G sont les matrices de masse, d'amortissement et raideur généralisées. Pour
simplifier, on considère qu'elles sont diagonales. La matrice d'amortissement généralisée CG aussi
car on suppose que l'hypothèse de Basile est vérifiée (la matrice d'amortissement est une combinaison
linéaire des matrices masse et raideur).
Soit, en ne considérant que les degrés de liberté actifs :
m .! (
q t) + c . !(
q t) + k . (
q t)
T
= - .f
T
c - .(m. + m . I
xs
)d x
G
G
G
! s
En l'absence de choc, on est donc conduit à résoudre un ensemble d'équations découplées (il y en a
autant que de modes propres).
Remarque :
Il est possible de calculer une base modale avec des matrices non diagonales. Il suffit de le
préciser lors de la construction de la numérotation généralisée par le mot clé
STOCKAGE : `PLEIN' de la commande NUME_DDL_GENE [U4.55.07].

Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
12/20
4.2
Choix de la base modale
Pour l'analyse sismique d'une structure linéaire, il faudrait en principe retenir tous les modes dont les
fréquences propres sont inférieures à la fréquence de coupure (généralement de l'ordre de 33 Hz).
Dans la pratique, on se contente souvent de ne conserver dans la base modale que les modes qui
contribuent de façon significative à la réponse. On conserve alors uniquement les modes dont la masse
effective unitaire dans une direction est supérieure à 1 et on s'assure également que, pour
l'ensemble de ces modes retenus, la masse effective unitaire cumulée dans chaque direction est peu
différente de la masse totale de la structure (supérieure à 90%). Le critère de cumul des masses
modales effectives est atteint en enchaînant les opérateurs suivants :
· Calcul de la masse totale de la structure : POST_ELEM [U4.61.04]
masse_in = POST_ELEM (MASS_INER : (TOUT : `OUI'))
· Calcul des modes propres dynamiques bloqués
: ils sont calculés dans l'opérateur
MODE_ITER_SIMULT [U4.52.02] ou dans MODE_ITER_INV [U4.52.01] selon la méthode choisie.
mode = MODE_ITER_SIMULT(); ou mode = MODE_ITER_INV();
· Normalisation des modes par rapport à la masse généralisée : NORM_MODE [U4.64.02]
NORM_MODE(MODE : mode, NORME :'MASSE_GENE', MASSE_INER : masse_in) ;
· Extraction de la base modale des modes dont la masse effective unitaire dépasse un certain seuil
(1 par exemple) et vérification que les modes extraits représentent au moins 90% de la masse
totale : EXTR_MODE [U4.64.03]
EXTR_MODE(
FILTRE_MODE(MODE: mode, CRIT_EXTRE: `MASSE_EFFE_UN', SEUIL :1.e-3)
IMPRESSION(CUMUL :'OUI') ;
Remarque :
La macro commande MACRO_MODE_MECA [U4.52.05] permet d'enchaîner directement l'ensemble
des trois dernières commandes précédentes.

Attention, certaines réponses locales (dans le cas particulier de non linéarités localisées) peuvent être
fortement influencées par des modes d'ordre supérieur dont la fréquence est au delà de la fréquence
de coupure et dont la masse modale effective est faible (inférieure à 1). Le mot clé VERI_CHOC de la
commande DYNA_TRAN_MODAL [U4.54.03] permet de vérifier a posteriori que la base modale choisie
est suffisante . Si ce n'est pas le cas, on conseille vivement de la compléter.
4.3 Calcul de la réponse dynamique de la structure étudiée par
synthèse modale
Après avoir calculé la base des modes propres dynamiques et construit une numérotation généralisée
par NUME_DDL_GENE [U4.55.07], on projette ensuite les matrices de masse, d'amortissement et de
raideur, sur cette même base avec l'opérateur PROJ_MATR_BASE [U4.55.01], les vecteurs second
membre avec PROJ_VECT_BASE [U4.55.02].
Remarque :
La macro commande MACRO_PROJ_BASE [U4.55.11] permet d'enchaîner directement l'ensemble
des trois opérations.

Les matrices et vecteurs ainsi projetés, on calcule la réponse généralisée du système mono ou
multi-excité à l'aide de l'opérateur DYNA_TRAN_MODAL [U4.54.03].
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
13/20
4.4
Prise en compte des modes négligés par correction statique
Lors du calcul de la réponse généralisée d'une structure mono excitée, il est possible de prendre en
compte, a posteriori, l'effet statique des modes négligés. Dans ce cas, une fois revenu sur la base
physique on corrige la valeur du déplacement relatif calculé (respectivement la vitesse relative et
l'accélération relative) par la contribution d'un pseudo-mode. Le pseudo-mode est défini par la
différence entre le mode statique associé au chargement unitaire de type accélération uniforme
imposée et la projection sur les modes dynamiques calculés du déplacement (respectivement la
vitesse relative et l'accélération relative).
On a alors :


p

X
= X +
f (t).
i
i
-
.
r_corrigé
r


j
j
i


j=1




p


!X
f t

r_corrigé = !
X + ! ( ).
r
i
i - ! .
j
j

i


j=1




p

!X
f t

r_corrigé = !X + ! ( ).
r
i
i - ! .

j
j
i


j=1



Les fonctions multiplicatives du temps f (t
i
) correspondent à l'accélérogramme imposé i(t) dans
chaque direction i considérée.
La démarche à suivre est la suivante :
· Calcul du chargement unitaire de type force imposée (accélération uniforme) dans la direction du
séisme : AFFE_CHAR_MECA [U4.25.01]. On fera attention à permuter le signe de la direction
puisque la force d'inertie sismique est de la forme (
P t) = -M .
. X!s
cham_no = AFFE_CHAR_MECA(MODELE:modèle, PESANTEUR:(VALE,DIRECTION)) ;
· Calcul de la réponse statique linéaire de la structure au précédent cas de charge :
MACRO_ELAS_MULT [U4.31.03].
mode_cor = MACRO_ELAS_MULT(CHAR_MECA_GLOBAL : con_lim, ...
CAS_CHARGE : (NOM_CAS :'xx', CHAR_MECA : cham_no)) ;
On notera qu'il y a autant de cas de charge que de direction de séisme
· Calcul des dérivées premières et seconde de l'accélérogramme : CALC_FONCTION [U4.62.04].
deri_pre et deri_sec = CALC_FONCTION ( OPTION : DERIVE) ;
· Calcul de la réponse généralisée en prenant en compte les modes négligés par correction
statique :
dyna_mod = DYNA_TRAN_MODAL (MASS_GENE : , RIGI_GENE :
MODE_CORR : mode_cor
EXCIT (CORR_STAT : `OUI'
D_FONC_DT : deri_pre, D_FONC_DT2 : deri_sec, ..)
...) ;
Remarque :
Dans le cas d'une structure multi excitée, la prise en compte des modes négligés par correction
statique n'est pas développée. On postraite le déplacement absolu dans ce cas.

Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
14/20
4.5
Prise en compte du caractère multi - supporté d'une structure
On a vu précédemment (cf. [§3.3]) que pour calculer le chargement sismique dans le cas d'une
structure multi supportée, il faut au préalable avoir calculé les modes statiques. Si on veut pouvoir
restituer les grandeurs calculées dans le repère absolu ou si on veut pouvoir prendre en compte des
non linéarités localisées, il faut également préciser dans DYNA_TRAN_MODAL que la structure étudiée
est multi excitée. En effet, dans ce dernier cas, on compare à chaque instant, le vecteur des
déplacements absolus de chacun des points de choc considérés, afin de déterminer si il y a choc et de
calculer les forces de choc correspondantes.
La démarche à suivre est la suivante :
· Calcul des modes statiques : MODE_STATIQUE [U4.52.04].
mode_stat = MODE_STATIQUE(DDL_IMPO : (...)) ;
· Calcul de la réponse généralisée en prenant en compte la composante d'entraînement :
dyna_mod = DYNA_TRAN_MODAL (MASS_GENE : , RIGI_GENE :
MODE_STAT : mode_stat
EXCIT (MULT_APPUI : `OUI'
ACCE : accelero, VITE : vitesse, DEPL : deplace
DIRECTION : (...), NOEUD :NO1
...)
...) ;
4.6 Post
traitements
Les opérateurs REST_BASE_PHYS [U4.64.01] ou RECU_FONCTION [U4.62.03] peuvent alors restituer
dans l'espace physique les évolutions calculées :
· l'opérateur
REST_BASE_PHYS restitue globalement (le champ complet) les déplacements,
vitesses et accélérations ;
· l'opérateur
RECU_FONCTION restitue localement (évolution temporelle d'un degré de liberté)
les déplacements, vitesses et accélérations.
On peut restituer les grandeurs relatives en précisant (MULT_APPUI : `NON') ou les grandeurs
absolues par (MULT_APPUI : `OUI').
On obtient alors les déplacements d'entraînement nécessaires au calcul des grandeurs secondaires en
soustrayant aux déplacements absolus les déplacements relatifs. Ceci est effectué par la commande
CALC_FONCTION [U4.62.04] option COMB.
A partir des évolutions précédentes, on peut également extraire les valeurs max. et RMS et calculer le
spectre de réponse d'oscillateur associé. Ceci est effectué par la commande CALC_FONCTION options
MAX, RMS et SRO.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
15/20
5
Réponse sismique transitoire directe
L'intégration directe est réalisable soit avec des hypothèses de comportement linéaire : opérateur
DYNA_LINE_TRAN [U4.54.01] soit avec des hypothèses de comportement non linéaire : opérateur
DYNA_NON_LINE [U4.32.02]. Mise à part la façon de prendre en compte le chargement sismique
(cf. [§3.3]), les syntaxes de DYNA_NON_LINE et DYNA_LINE_TRAN sont identiques.
5.1 Prise en compte d'un amortissement équivalent à l'amortissement
modal
Généralement, l'information la plus précise que l'on ait sur l'amortissement provient des essais de
vibration qui permettent de déterminer, pour une fréquence de résonance donnée fi , la largeur de
résonance correspondante et donc l'amortissement réduit i à cette résonance. Il est donc
nécessaire de pouvoir prendre en compte, dans un calcul transitoire direct, un amortissement
équivalent à l'amortissement modal
.
A partir du développement spectral de la matrice identité :
n_mod es X XT K n_modes X XT K
Id =
i
i
=
i
i
T

2
i 1
X K X
i 1
M
=
i
i
=
G_i .i
on montre :
· qu'on peut développer la matrice d'amortissement de la structure C en série de modes
propres :
n_mod es
T
C = ai .(K. i
)(K. i
)
i=1
· et que, compte tenue de la définition du pourcentage d'amortissement critique :

T
i
i . .
C i = .
2 MG_i .i . i .ai = .
2 KG_i.i
Il est donc conseillé à l'utilisateur de spécifier (les syntaxes de DYNA_NON_LINE et DYNA_LINE_TRAN
sont identiques), les valeurs des amortissements modaux pour chaque fréquence propre par
l'intermédiaire du mot clé facteur AMOR_MODAL.
Cela revient à imposer une force d'amortissement proportionnelle à la vitesse relative de la structure :
n_mod es

F
= C X
i
T
amo
! r avec C = 2.
.(K. i
)(K. i
)
K
i=
.
G i
1
_
i
5.2 Prise en compte d'une sollicitation multi appuis avec restitutions
des champs relatifs et absolus
Par défaut, les grandeurs sont calculées dans le repère relatif. Dans DYNA_NON_LINE et
DYNA_LINE_TRAN, on utilise une syntaxe identique à celle de DYNA_TRAN_MODAL (présence des mots
clés MODE_STAT et MULT_APPUI : `OUI') pour les calculer dans le repère absolu.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
16/20
6 Interaction
sol-structure
Le comportement sismique d'un bâtiment dépend des caractéristiques du sol sur lequel il est posé
puisqu'il dépend du mouvement sismique imposé au sol et du comportement dynamique du bâtiment et
de ses fondations. L'interaction sol-structure contribue le plus fréquemment à diminuer la réponse de la
structure étudiée.
6.1
Impédance d'une fondation
Soit une fondation rigide superficielle sans masse, soumise à une force harmonique de pulsation :
P t
it
( ) = P .e
0
Elle est donc animée d'un mouvement X (t) de même fréquence. On appelle
impédance de la fondation, le nombre complexe (
K ) , fonction de la fréquence tel que :
(
P(t)
K ) =
.
X (t)
Plusieurs méthodes analytiques ou numériques permettent de calculer l'impédance d'une fondation
selon la complexité de la fondation et du sol sur lequel elle est posée ou partiellement enfouie. Parmi
les plus fréquemment utilisées, on cite :
· les méthodes analytiques des ressorts de WOLF ou de DELEUZE où l'on suppose que le
radier est circulaire, rigide et posé sur un sol homogène. La fondation doit être superficielle ;
· la méthode numérique du code CLASSI où l'on suppose que le radier est de forme
quelconque, rigide et posé sur un sol éventuellement stratifié. La fondation doit être
superficielle ;
· la méthode numérique du code MISS3D où le radier peut être de forme quelconque,
éventuellement déformable et posé sur un sol éventuellement stratifié.
Il est possible de traiter l'interaction sol-fondation par la méthode fréquentielle de couplage (prise en
compte de la dépendance en fréquence de la matrice d'impédance) en réalisant un calcul couplé
MISS3D / Code_Aster. Ce type de calculs n'est pas détaillé dans cette documentation de référence. On
présente ici uniquement le cas le plus courant où l'interaction sol-fondation est traitée par la méthode
des ressorts de sol
(on considère que les termes de la matrice d'impédance sont indépendants de la
fréquence).
Dans le cas d'une fondation rigide superficielle, l'impédance est calculée au centre de gravité de la
surface en contact dans un repère lié aux axes principaux d'inertie de cette surface. Pour chaque
fréquence, elle s'exprime sous la forme d'une matrice de dimension (6, 6). On ajuste ensuite la valeur
de chaque terme en fonction d'un mode propre particulier du bâtiment étudié en base bloquée :
· fréquence du premier mode de balancement 0 pour les raideurs horizontales
Kx(0), Ky(0) et de rotation Krx(0), Kry(0) ;
· fréquence du premier mode de pompage

1 pour la raideur verticale Kz( 1) et de torsion
Krz(1).
Comme les fréquences propres du bâtiment dépendent des raideurs de sol, le calcul des valeurs
globales des six ressorts de sol résulte d'un processus itératif illustré figure [Figure 6.1-a]. Les
premières raideurs de sol Kx(0 ), Ky(0), Kz(1), Krx(0), Kry(0 ) et Krz(1) sont choisies
en fonction des premières fréquences propres de balancement (0) et de pompage (1) de la
structure en base bloquée. Les raideurs de sols sont ensuite ajustées aux premières fréquences
propres significatives de la structure sur ressort jusqu'à correspondance des fréquences auxquelles les
fonctions d'impédance sont calculées avec les valeurs des fréquences propres du système couplé
sol-bâtiment.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
17/20
Sélection des raideurs de sol
kv(j) et k(i) dans la matrice
d'impédance
Code_Aster
Calcul modal avec structure sur
ressort de sol.
Non
Les fréquences propres
du modèles sur ressort
correspondent aux
fréquences i et j
Oui
Fin
Figure 6.1-a : Processus d'ajustement des raideurs de sol
6.2 Prise en compte d'un amortissement modal calculé selon la règle
du RCC-G
On décompose l'amortissement dû au sol en une partie d'origine matérielle et en une partie
géométrique : amortissement dû à la réflexion des ondes élastiques dans le sol.
La règle du RCC_G consiste à sommer, pour chaque mode, les amortissements de chaque sous
structure constitutive du bâtiment considéré et les amortissements structurel et géométrique du sol
pondérés par leur taux respectifs d'énergie potentielle par rapport à l'énergie potentielle totale :
Eki.k + Esi.si

k
s
i =
Eki + Esi
k
s
avec :
·
i , l'amortissement réduit moyen du mode i ;
·
k , l'amortissement réduit du kème élément de la structure ;
·
si , l'amortissement réduit du ressort de sol s pour le mode i ;
·
Eki , l'énergie potentielle du kème élément de la structure pour le mode i ;
· et
Esi , l'énergie potentielle du ressort de sol s pour le mode i.
Dans le règlement, l'amortissement modal est limité à une valeur maximale de 0,3.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
18/20
La partie d'origine matérielle de l'amortissement du sol est calculée en pondérant l'amortissement de
chaque sous structure par le rapport : taux d'énergie potentielle sur énergie potentielle totale. Quant à
la partie géométrique de l'amortissement, elle est calculée en répartissant les valeurs d'amortissement
pour chaque direction (trois translations et trois rotations) pondérées par le taux d'énergie potentielle
dans le sol de la direction. Les valeurs d'amortissement directionnelles sont obtenues en interpolant,
pour chaque fréquence propre calculée, les fonctions d'amortissement directionnelles issue d'un code
d'interaction sol-structure (PARASOL, CLASSI ou MISS3D). Le rapport de la partie imaginaire sur deux
Im(K())
fois la partie réelle de la matrice d'impédance :
, fournit les valeurs de cet amortissement
.R
2 e(K())
radiatif.
La démarche à suivre est la suivante :
· Calcul de l'énergie potentielle dissipée dans la structure étudiée : POST_ELEM [U4.61.04]
Ek = POST_ELEM (ENER_POT: (TOUT : `OUI')) ;
· Calcul de l'amortissement modal par la règle du RCC_G : CALC_AMOR_MODAL [U4.64.04]
l_amor = CALC_AMOR_MODAL (
ENER_SOL : (MODE_MECA : base_modale, GROUP_NO_RADIER : ,
Kx(
Ky(
Kz(1)
0 )
0 )
KX :
, KY :
, KZ :
,
Krx(
Kry(
Krz(1)
0 )
0 )
KRX :
, KRY :
, KRZ :
)) ;
AMOR_INTERNE : (GROUP_MA :, ENER_POT : E k , AMOR_REDUIT : k )
Im(K())
AMOR_SOL : (FONC_AMOR_GEO :
.R
2 e(K()) )
) ;
Le calcul de la contribution du sol à l'énergie potentielle Es (mot clé facteur ENER_SOL) est calculée à
partir des valeurs d'impédance de sol déterminées précédemment (cf. [§6.1]). Elle peut être calculée
selon deux méthodes différentes selon que l'on moyenne les efforts modaux (mot clé RIGI_PARASOL)
ou les déplacements modaux aux noeud du radier.
L'amortissement réduit du ressort de sol s s (mot clé facteur AMOR_SOL) est calculé à partir des
valeurs de l'amortissement radiatif.
6.3
Répartition des raideurs et amortissement de sol
Si l'on veut étudier l'effet d'un séisme sur le décollement éventuel du radier par exemple, on peut être
amené à modéliser le sol non plus par un ressort unique au centre de gravité de l'interface sol-bâtiment
mais par un tapis de ressorts. Ceci est possible grâce à la commande AFFE_CARA_ELEM [U4.24.01]
option RIGI_PARASOL.
La démarche consiste à calculer en chaque noeud du maillage du radier les raideurs élémentaires
(k ,k ,k ,kr ,kr ,kr
x
y
z
x
y
z ) à appliquer à partir des valeurs globales des trois ressorts de translations :
kx, ky, kz et des trois ressorts de rotations : krx, kry, krz issues d'un code d'interaction sol-structure
(ou calculées analytiquement).
On suppose que les raideurs élémentaires de translation sont proportionnelles à la surface S(P)
représentée par le noeud P et à une fonction de répartition f (r) dépendant de la distance r du noeud
P au centre de gravité du radier O :
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
19/20

K = k
.
x
x ( P) = k
S( P). f ( P
x
O )

P
P

K = k
y
y ( P) = k . S( P). f (OP
y
)

P
P
K = k
z
z ( P) = k . S( P). f (OP
z
)


P
P
On en déduit alors k x puis k ( P
x
) à partir du calcul :
S( P). f (OP)
kx (P) = k .S(P). f ( P
O ) = K
x
x .
.
S( P). f (OP)
P
On en déduit de même k ( P
y
) et k (P
z
).
Pour les raideurs élémentaires de rotation, on répartit ce qui reste après avoir enlevé les contributions
dues aux translations de la même façon que les translations :

K = k
2
2
2
2
rx
rx (P) + [k y (P).zO + k
P
z (P). yOP ] = k . S( P). f (OP
rx
) + [ky(P).zO + k
P
z (P). y P
O ]

P
P
P
P

K = k
2
2
2
2
ry
ry (P) + [kx (P).zO + k
P
z (P). xOP ] = k . S( P). f (OP
ry
) + [kx(P).zO + k
P
z (P). x P
O ]

P
P
P
P
K = k
2
2

=
. ( ). (O ) +
2
2

k
S P f
P
[kx(P).y PO +ky(P).x
rz
OP ]
rz
rz (P) +
[kx(P).yO +k
P
y (P). x P
O ]

P
P
P
P
On en déduit alors krx puis k ( P
rx
) à partir du calcul :
krx (P) = k .S(P). f (OP
rx
)

S(P). f (OP)
= K
. 2
. 2
rx - [k y (P) zO + k
P
z (P) y P
O ] .

S(P). f (OP
P
)
P
On en déduit de même k ( P
ry
) et k (P
rz
) .
Remarque :
Par défaut, on considère que la fonction de répartition est constante et unitaire c'est-à-dire que
chaque surface est affectée du même poids.

On peut répartir de même six valeurs globales d'amortissement, analytiques ou calculées par un
code d'interaction sol-structure.

6.4
Prise en compte d'une frontière absorbante
Si l'on veut calculer la réponse sismique d'un barrage, il faut, entre autre, pouvoir prendre en compte la
non réflexion des ondes dans la vallée. Ceci est possible grâce à des éléments à frontière absorbante :
option IMPE_ABSO dans DYNA_NON_LINE et DYNA_LINE_TRAN. Cette fonctionnalité n'est pas
détaillée dans le présent document. Elle fera l'objet d'une note spécifique.
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Code_Aster ®
Version
5.0
Titre :
Réponse sismique par analyse transitoire
Date :
28/09/99
Auteur(s) :
Fe WAECKEL
Clé :
R4.05.01-A
Page :
20/20
7 Bibliographie
[1]
R.W. CLOUGH, J. PENZIEN - Mc GRAW-HILL : « Dynamics of structures » - (1975).
[2]
« Génie parasismique » - ouvrage collectif sous la direction de V. DAVIDOVICI - Presses de
l'E.N.P.C. (1995).
[3]
P. LABBE, A. PECKER, J.P.
TOURET
: «
Comportement sismique des structures
industrielles » - cours IPSI du 20 au 22 septembre 1994.
[4]
Fe WAECKEL : « Méthode pour le calcul par superposition modale de la réponse sismique
d'une structure multi-supportée » - Note HP62/95.017B (09/95).
[5]
V. GUYONVARH, G. DEVESA : « Méthode pour le calcul des excitations sismiques aux
ancrages du CPP N4 » - Note HP52/99.006 (09/99).
[6]
« Calcul sismique des bâtiments » - Règles RCC-G, annexe A
Manuel de Référence
Fascicule R4.05 : Analyse sismique
HP-52/99/017 - Ind A

Document Outline