Bâtiment En Cisaillement–SHEAR BUILDING
REPONSE SISMIQUE MODALE
ENCADRÉ PAR:
RÉALISÉ PAR:
DAKDAKI Soukaina
[Link] ESSALIH Adam
Oussama BHIJA Samira
EL KHADIR Morad
KONATÉ Mory
PLAN
Construction des matrices Réponse sismique
01 02 03 04 05
Modélisation physique Analyse modale Validation
But et idée générale:
On remplace le bâtiment par un système masse-
ressort-amortisseur à un degré de liberté par étage :
chaque étage est une masse concentrée mi, reliée
horizontalement aux autres par des ressorts de
raideur kiet amortisseur Ci . On cherche M, K, C et
l’équation qui gouverne les déplacements horizontaux
X(t).
L’objectif final est :
prédire comment un bâtiment
à plusieurs étages se déforme
et oscille lorsqu’il subit un
tremblement de terre.
Hypothèses
-Masses concentrées &
-Cisaillement pur : uniformes :
déplacements horizontaux Ressorts linéaires :
masse d’un étage concentrée en
uniquement (chaque étage se relations linéaires de Hooke
un nœud (pas de distribution
déplace comme un corps rigide F=k Δx
continue).
horizontalement).
m1=m2=m3=m4=m.
Base (sol) : Amortissement
Raideurs uniformes : le sol se déplace avec une proportionnel (Rayleigh) :
k1=k2=k3=k4=k accélération donnée u¨ (base C=αM+βK pratique et usuel pour
excitation). l’analyse modale.
Équations du mouvement:
1. Choix des coordonnées
déplacement horizontal relatif de l’étage i par
On définit :
rapport au sol.
Le vecteur des déplacements :
2. Énergie cinétique (inertie)
Chaque étage possède une masse mi, donc son énergie cinétique est :
Sous forme matricielle :
Cette expression donne après dérivation :
Ce terme correspond à l’inertie du système.
Équations du mouvement:
3. Énergie potentielle élastique (raideur)
Chaque liaison ressort entre les étages i et
i+1 stocke une énergie :
En sommant tous les ressorts :
La dérivée par rapport à x donne :
Ce terme représente les forces de raideur, qui
s’opposent au déplacement.
Équations du mouvement:
4. Damping (amortissement)
On modélise les pertes d’énergie par un
amortissement visqueux :
La fonction de dissipation de Rayleigh :
donne : terme d’amortissement proportionnel à la vitesse
5. Forces extérieures : mouvement du sol
Le bâtiment repose sur un support soumis à une
accélération du sol.
Chaque masse subit une force inertielle :
Équations du mouvement:
Sous forme matricielle :
6. Application de l’équation de Lagrange
Formulation :
Remplacement des termes :
PLAN
Construction des matrices Réponse sismique
01 02 03 04 05
Modélisation physique Analyse modale Validation
Construction de la matrice de masse M :
La matrice de masse est diagonale car chaque étage est
modélisé par une masse concentrée indépendante. Il n’existe
pas de couplage inertiel entre étages, donc seuls les termes
mi apparaissent sur la diagonale.
Matrice de raideur K :
La matrice K est obtenue en appliquant la loi
de Hooke sur chaque liaison : chaque ressort
produit une force proportionnelle au
déplacement relatif entre deux étages.
La diagonale contient la somme des raideurs
connectées à l’étage, et les termes −ki hors
diagonale représentent le couplage entre
étages.
Amortissement Rayleigh (choix de α et β)
L’amortissement du bâtiment est modélisé avec le modèle de Rayleigh :
C=αM+βK Ce modèle suppose que l’amortissement dépend :
de la masse du système (effets visqueux globaux) → terme αM
de la raideur du système (frottements internes dans la structure) →
terme βK
Les coefficients α et β permettent de
représenter l’amortissement du bâtiment en
combinant les effets liés à la masse et à la
raideur.
Comment on choisit α et β ?
On choisit α et β pour obtenir un taux
d’amortissement modal souhaité sur un ou
plusieurs modes :
équivalemment
Dans notre projet, on a choisi :
α= 3.57326e-08
β= 3.06147e-04
PLAN
Construction des matrices Réponse sismique
01 02 03 04 05
Modélisation physique Analyse modale Validation
Analyse Modale:
pour comprendre comment le bâtiment va vibrer naturellement On
calcule les fréquences propres et les modes propres.
Matrices données
On résout le problème généralisé
ou, équivalemment
Pour qu’il existe une solution non triviale
Analyse Modale:
Calculer le déterminant:
Comme la matrice est tridiagonale, le déterminant se
calcule par expansion successives
On note :
Analyse Modale:
Donc la matrice devient :
On utilise le développement récursif :
Ici :
C’est le principe général des matrices
tridiagonales.
À la fin on obtient un polynôme du 4ᵉ degré en λ
Analyse Modale:
Voici les résultats obtenus à partir de l’analyse modale de la structure:
[ ]
-0.0158 0.0207 0.0158 -0.0086
= -0.0158 0.0086 -0.0158 0.0207
-0.0158 -0.0086 -0.0158 -0.0207
-0.0158 -0.0207 0.0158 0.0086
1.0 0 0 Les résultats montrent un mode rigide à ω = 0 (translation
uniforme), ce qui signifie que la base n’est pas encastrée dans ce
modèle — il faut vérifier la condition aux limites. Le premier mode
0.2929 54.1196 8.6134 propre non nul a une fréquence de 8.61 Hz et représente la
déformation globale du bâtiment ; c’est généralement celui qui
domine la réponse si le séisme contient de l’énergie proche de cette
1.0000 100.0000 15.9155
fréquence. Les modes suivants (≈ 15.9 Hz et 20.8 Hz) sont des
modes d’ordre supérieur présentant plus de nœuds et se manifestent
pour excitations à plus haute fréquence.
1.7071 130.6563 20.7946
Analyse Modale:
Voici les résultats obtenus à partir de l’analyse modale de la structure:
[ ]
-0.0158 0.0207 0.0158 -0.0086
= -0.0158 0.0086 -0.0158 0.0207
-0.0158 -0.0086 -0.0158 -0.0207
-0.0158 -0.0207 0.0158 0.0086
1.0 0 0 Les résultats montrent un mode rigide à ω = 0 (translation
uniforme), ce qui signifie que la base n’est pas encastrée dans ce
modèle — il faut vérifier la condition aux limites. Le premier mode
0.2929 54.1196 8.6134 propre non nul a une fréquence de 8.61 Hz et représente la
déformation globale du bâtiment ; c’est généralement celui qui
domine la réponse si le séisme contient de l’énergie proche de cette
1.0000 100.0000 15.9155
fréquence. Les modes suivants (≈ 15.9 Hz et 20.8 Hz) sont des
modes d’ordre supérieur présentant plus de nœuds et se manifestent
pour excitations à plus haute fréquence.
1.7071 130.6563 20.7946
PLAN
Construction des matrices Réponse sismique
01 02 03 04 05
Modélisation physique Analyse modale Validation
Réponse sismique
La structure est soumise à une accélération du sol représentative
d’un séisme, et nous étudions la réponse dynamique obtenue afin
d’évaluer les déplacements et l’amplification au niveau des étages
Nous avons utilisé Python pour simuler la réponse du bâtiment
Le graphe représente l’évolution du déplacement
horizontal du dernier étage du bâtiment pendant 20
secondes d’excitation sismique.
Réponse sismique
Le spectre (Fast Fourier Transform - FFT) de l'accélération du sol montre
plusieurs pics d'amplitude importants à très basse fréquence (proche de 0 l'accélération du sol
Hz) et deux pics nets à des fréquences plus élevées, suivis d'un bruit de
fond.
Réponse sismique
Analyse et Interprétation des Modes Propres de la Structure Cette diapositive confirme que seul le
mode 1 contribue réellement à la réponse
dynamique. Son amplitude dépasse 150 m
alors que les autres modes sont
négligeables. La structure se comporte
donc comme un corps rigide qui glisse
Réponse sismique
Analyse fréquentielle du déplacement : dominance du
Contributions modales au déplacement du dernier étage corps rigide
Le spectre montre un pic massif à très basse fréquence (≈ 0 Hz).
Cela confirme que le déplacement est dominé par un mouvement
de glissement (mode rigide).
PLAN
Construction des matrices Réponse sismique
01 02 03 04 05
Modélisation physique Analyse modale Validation
Validation
1)Description du modèle
Type de structure :Bâtiment à ossature (4 étages)
Matériau : [béton/acier – préciser]
Type d’éléments : [coques, poutres, volumes]
Conditions limites : Encastrement au niveau du sol
2)Objectif de la simulation
Vérifier le comportement dynamique du bâtiment
Déterminer les fréquences propres et les modes de vibration
Identifier les zones sensibles (déformations maximales)
3) Paramètres d’analyse Index Description Valeur
Fréquence
(cycles/time)
Méthode utilisée : Analyse modale 0 Base State – –
Nombre de modes extraits: 10 modes 1 Mode 1 0.97986 0.15754
2 Mode 2 0.97986 0.15754
4)Résultats principaux 3 Mode 3 1.7583 0.21104
Fréquences propres : 4 Mode 4 8.1546 0.45449
Mode 1 : ~ [X] Hz (translation) 5 Mode 5 10.584 0.51778
Mode 2 : ~ [X] Hz (translation 6 Mode 6 10.584 0.51778
perpendiculaire) 7 Mode 7 17.742 0.67038
Mode 3 : ~ [X] Hz (torsion) 8 Mode 8 21.316 0.73481
Mode 4 : ~ [X] Hz (flexion 9 Mode 9 37.315 0.97222
10 Mode 10 37.315 0.97222
supérieur)
La visualisation présente l’animation des
modes propres obtenus, permettant
Les résultats numériques ont ensuite été
d’observer la déformation dynamique et
exportés au format .csv (fréquences modales,
la distribution du déplacement à travers
vecteurs propres, déplacements nodaux…) afin
la structure.
de les exploiter sous Python pour le tracé des
courbes et l’analyse comparative des modes.
les déplacement modale
Mode 1 – Translation globale
Déplacement uniforme de tous les étages.
Mouvement latéral global, sans changement
de phase.
Mode dominant dans la réponse sismique.
Capte la plus grande part de la masse
structurelle.
Mode 2 – Translation perpendiculair
Déformation plus courbée que le mode 1.
Apparition d’une première zone d’inflexion.
Contribution modérée à la dynamique globale.
Influence importante sur les efforts internes.
Mode 3 – Torsion
Forme comportant plusieurs courbures.
Étages pouvant entrer en opposition relative.
Contribution faible à la masse globale.
Indicateur des zones sensibles et effets locaux.
Mode 4 – Flexion supérieure (haut ordre)
Déplacement global sans inversion de phase.
Distribution quasi uniforme le long de la hauteur.
Représente une flexion d’ordre supérieur dans la
même direction que le mode 1.
Contribue faiblement à la masse mais important pour
les effets locaux
Comparaison des fréquences propres (critère principal)
Mode Fréquence Python (Hz) Fréquence Abaqus (Hz) Écart (%)
1 0 0.15754 non défini
2 6.0906 0.15754 97 ,41
3 11.254 0.21104 98 ,13
4 14.704 0.45449 96 ,91
Les fréquences propres d’Abaqus sont plus élevées que celles du modèle masse-ressort Python
Cette différence est normale et due à :
la simplification du modèle masse-ressort (masses ponctuelles, rigidité réduite),
la modélisation complète en 3D dans Abaqus (poutres, colonnes, inerties, torsion),
la rigidité des planchers et la masse distribuée,
les conditions d’encastrement plus réalistes.
Résultats obtenus
Même ordre des modes (translation, translation perpendiculaire, torsion, flexion).
Formes modales similaires entre Python et Abaqus.
PLAN
Construction des matrices Réponse sismique
01 02 03 04 05
Modélisation physique Analyse modale Validation
Reponse Sismique
La structure est soumise à une accélération du sol représentative
d’un séisme, et nous étudions la réponse dynamique obtenue afin
d’évaluer les déplacements et l’amplification au niveau des étages
Reponse Sismique
Contribution Modale Tendance décroissante continue (une
dérive négative progressive) sur
laquelle sont superposées de petites
oscillations périodiques.
Reponse Sismique
Déplacements Modaux
Cette figure montre la réponse sismique
d’un bâtiment à 4 étages, représentée
par les déplacements modaux
q1(t)q_1(t)q1( t), q2(t)q_2(t)q2( t),
q3(t)q_3(t)q3( t) et q4(t)q_4(t)q4( t)
sur 10 secondes :
Reponse Sismique
Déplacements Modaux
Mode 1 (0 Hz) : déplacement lent et
croissant, correspondant au
mouvement rigide global du
bâtiment.
Reponse Sismique
Déplacements Modaux
Mode 2 (≈ 2.18 Hz) : oscillations
régulières de grande amplitude,
dominant la réponse.
Reponse Sismique
Déplacements Modaux
Mode 3 (≈ 4.03 Hz) : oscillations
plus rapides, d’amplitude moyenne et
décroissante.
Reponse Sismique
Déplacements Modaux
Mode 4 (≈ 5.26 Hz) : oscillations
très rapides, de faible amplitude,
fortement amorties.
Reponse Sismique
Analyse modale
Ce graphe Représente les formes
modales mass-normalisées d’un
bâtiment de 4 étages lors d’une
analyse de réponse sismique.
Chaque courbe montre la manière dont
chaque étage se déforme dans un mode
de vibration donné
Reponse Sismique
Analyse fréquentielle Le bâtiment amplifie surtout les
basses fréquences, là où se
trouvent ses modes de vibration,
alors que les hautes fréquences
sont fortement filtrées.
Reponse Sismique
Acceleration sismique Le graphe représente
l’accélération sismique du sol
utilisée comme excitation
d’entrée. Le signal est
irrégulier et reflète les
variations rapides typiques d’un
séisme, qui induisent les
vibrations et réponses dynamiques
des différents étages.
Reponse Sismique
Animation 3D
Cette
accélération
sismique du sol
sert d’entrée
dynamique au
modèle du bâtiment
à 4 étages et
provoque la
vibration
successive des
étages selon leur
rigidité et leurs
modes propres.
Reponse Sismique
Résumé des résultats de la simulation
--- Parametres generaux --- --- Amortissement Rayleigh ---
Nombre d'etages : 4 alpha (M) : 8.46888e-09
Masse par étage (kg) : 25000.0 beta (K) : 1.21015e-03
Taux cible zeta : 0.020 (≈ 2.0%)
Raideur k (N/m) : 8.000e+06
Pas de temps (s) : 0.0010 --- Deplacements ---
Duree simulation (s) : 10.00 Max |x| dernier étage : 9.5838e-01 m
Dérive finale (x(T)) : -9.5838e-01 m
Rapport drift/max : 1.000
--- Donnees modales ---
Freq. propres (Hz) : 0.000 2.179 4.026 5.261
Omega (rad/s) : 0.000 13.691 25.298 33.054
--- Contribution modale au dernier étage ---
Mode dominant (forme) : Mode 2 (2.179 Hz)
Energie relative modes : [ 1.00 0.00 0.00 0.00 ]
--- Verification matrices ---
Norme(M) : 2.500e+04
Norme(K) : 2.731e+07
Norme(C) : 3.305e+04
Conclusion
Ce projet avait pour objectif d’étudier la réponse sismique d’un
bâtiment. Après la modélisation géométrique sur CATIA, AutoCAD
et 3D Experience, nous avons développé un modèle dynamique
masse-ressort et réalisé la simulation numérique sous Python et
MATLAB. Les résultats montrent une amplification des
déplacements vers les étages supérieurs et mettent en évidence le
rôle du premier mode de vibration.
THANK
YOU