Etude des modes propres d’un bâtiment en béton armé

Présenter par : Sous la direction de :

LIPOBO MINGONDZA Sarah Dr. PAMBOU NZIENGUI
Claude. F
NZENGUE Claille Darlin
Table des matières
Introduction............................................................................................................................................3
I. Objectif de l’étude dynamique........................................................................................................3
II. Résolution de la première question................................................................................................4
1. Données :....................................................................................................................................4
2. Propositions d’un système dynamique équivalent......................................................................4
3. Equation dynamique du système sous l’effet de l’excitation sismique d’accélération ẍ g............4
4. Calcul des masses M1 et M2.........................................................................................................5
5. Déterminons les rigidités des niveaux 1 et 2...............................................................................6
5.1 Calcul des inerties dans les deux directions..........................................................................6
5.2 Rigidités des niveaux 1 et 2...................................................................................................7
III. Calcul des modes et des fréquences normalisées dans les deux directions................................8
1. Fréquences de vibrations/Valeurs propres.................................................................................8
2. Modes de vibrations/vecteurs propres dans les deux directions..............................................11
VI. Détermination des participations modales (suivant OX et OY)........................................................14
V. Représentation de la décomposition dans l’espace modale :.......................................................17
VI. Conclusion.................................................................................................................................18
Annexe mathématique.........................................................................................................................19
Algèbre matricielle............................................................................................................................19
Bibliographie.........................................................................................................................................24

1
 Darlin - Sarah

Ecole Polytechnique de Masuku (EPM) Génie Civil 5A

2023 -2024

PROJET 2
Etude des modes propres d’un bâtiment en béton armé
On se propose d’étudier les vibrations suivant OX et OY d’un bâ timent R+1 réalisé en béton armé
(Fig. 1), de module d’Young E=30 GPa et de masse volumique 2500 =* ‫ م‬kg/m3. Le RDC est
composé de 8 poteaux de section (0,2×0,3 m2) et d’un voile central de dimension (2×0,2 m2). En
R+1, tous les poteaux sont continus, mais le voile est transformé en un poteau de dimension
(0,3×0,2 m2). L’ensemble poteaux et voile, supportent, au RDC et en R+1, des dalles pleines et
rigides d’épaisseurs e=20 cm.

(a) (b)

Fig. 1 : Bâ timent en Béton R+1 étudié : (a) vue en plan R+1 ; (b) système dynamique équivalent

1. Proposez un système dynamique équivalent ; donner l’équation dynamique du système

sous l’effet d’une excitation sismique ; calculez les masses permanentes M1 et M2, les
rigidités des niveaux 1 et 2. Donnez les matrices globales : de masses et de rigidités dans
les deux directions (OX et OY).
2. Calculez les modes et les fréquences propres normalisées dans les deux directions
(Indications : calculer le déterminant de la matrice (K-𝜆M))
3. Calculez les participations modales (suivant OX et OY)
4. Donnez une représentation de décomposition dans l’espace modale
5. Donnez vos remarques et conclusions sur le comportement dynamique du bâ timent
suivant les deux directions.
« Rien n’est impossible à celui qui croit »

2
CFPN, Février 2024, Dynamique des Structures…
 Introduction
Parmi les catastrophes naturelles qui affectent la surface de la terre, les secousses sismiques sont
sans doute celles qui ont le plus d’effets destructeurs dans les zones urbanisées. Face à ce risque
et à l’impossibilité de le prévoir, il est nécessaire de construire des structures pouvant résister à
de tels phénomènes, afin d’assurer au moins une protection acceptable des vies humaines, d’où
l’apparition de la construction parasismique. Cette dernière se base généralement sur une étude
dynamique des constructions agitées.

I. Objectif de l’étude dynamique

L’objectif initial de l’étude dynamique d’une structure est la détermination des caractéristiques
dynamiques propres de la structure lors de ses vibrations. Une telle étude pour notre structure
telle qu’elle se présente est souvent très complexe c’est pourquoi on fait souvent appelle à des
modélisations qui permettent de simplifier suffisamment les problèmes pour permettre
l’analyse.

3
 II. Résolution de la première question
1. Données :
E = 30 GPa

2500 = ‫ م‬kg/m3

Voile S1 = 2*0.2 m2 au RDC

Poteau S2 = 0.3*0.2 m2 au R+1

Epaisseur des dalles e = 0.2 m

2. Propositions d’un système dynamique équivalent

u 1(t) u2(t)

m1 m2

3. Equation dynamique du système sous l’effet de l’excitation sismique

d’accélération ẍg
Pour ce faire utilisons le principe fondamental de la dynamique (PFD)

 1 ère équation :

M1ẍ 1 = -k1u1-k2 (u1-u2)-c1u̇1-c2(u˙1-u˙2)

M1ẍ 1 + (c1+c2)u˙1 – c2u˙2 + (k1 + k2) u1 – k2u2 = 0

Or, u1 = x1 – xg

⇨ x1 = u1 + xg
Alors, M1u¨1 + (c1+c2)u˙1 – c2u˙2 + (k1 + k2) u1 – k2u2 = - M1ẍg (1)

 2 ème équation :

M2ẍ 2 = -c2(u˙1-u˙2) – k2 (u2 – u1)

M2ẍ 2 –c1u˙1 + c1u˙2 – k2u1 + k2u2 = 0

Or, u2 = x2 – xg

⇨ x2 = u2 + xg
Alors, M2u¨2 – c1u˙1 + c2u˙2 – k2u1 + k2u2 = - M2 ẍg (2)

4
  Sous forme matricielle nous avons :

M1 0 u¨1 c 1+ c 2 −c 2 u˙1 k 1+k 2 −k 2 u1 −M 1 u¨1

+ + =
0 M2 u¨2 −c 1 c 2 u˙2 −k 2 k2 u2 −M 2 u¨2

4. Calcul des masses M1 et M2

Pour calculer les masses permanentes M1 et M2 nous négligerons par hypothèse le poids propre des
poteaux. Ainsi :

 Masse du voile au RDC = M1

- Surface d’influence : SI1

SI1 = Lx*Ly = 4.5 * 5 = 22.5 m2 = V1 = SI1 * e = 4.5 m3

D’où : M1 = SI1 * 11250 = ‫ م‬kg

 Masse du voile au R+1 = M2

- Surface d’influence : SI2

SI2 = Lx*Ly = 4.5 * 5 = 22.5 m2, il vient que M1 = M2 = 11250 kg

5
 5. Déterminons les rigidités des niveaux 1 et 2

Dans le cadre de notre projet nous supposerons des poteaux bis – encastré au niveau des dalles.
D’après la méthode de VERECHAGUINE

12 EI
K= 3
h
Source : Polycopié de dynamique des structures cours et exercices, page 30

Avec : E = 30 GPa = 30*109 N/m2

5.1 Calcul des inerties dans les deux directions

a. Suivant OX

3
2∗0.2 1 4
 Au RDC : I 1 x = = m
12 750

 Au R + 1 :
3
0.3∗0.2 −4 4
I 2 x= =2∗10 m
12

6
 b. Suivant OY

 Au RDC :
3
0.2∗2 2 4
I 1 y= = m
12 15
 Au R + 1 :
3
0.2∗0.3 −4 4
I 2 y= =4.5∗10 m
12

5.2 Rigidités des niveaux 1 et 2

a. Suivant OX

 Au RDC :

12 E∗I 1 x
K 1 x= 3
h1

9
12∗30∗10 ∗1 6
K 1 x= 3
=3.84∗10 N /m
5 ∗750

 AR+1:

12 E∗I 2 x
K 2 x= 3
h2
9 −4
12∗30∗10 ∗2∗10 6
K 2 x= 3
=2.67∗10 N /m
3

b. Suivant OY

 Au RDC :

12 E∗I 1 y
K 1 y= 3
h1

9
12∗30∗10 ∗2 8
K 1 y= 3
=3.84∗10 N /m
5 ∗15

7
  Au R + 1 :

12 E∗I 2 y
K 2 y= 3
h2

9 −5
12∗30∗10 ∗45∗10 6
K 2 y= 3
=6∗10 N /m
3

D’où la matrice de [M] et rigidité [K] :

Matrice de masse :

1 1250 0
[M] =
0 11 250

Matrice de rigidité :

K 1 x + K 2 x −K 2 x 6.51 −2.67
[K]OX = = * 106 [N/M]
−K 2 x K 2x −2.67 2.67

K 1 y+ K 2 y −K 2 y 390 −6
[K]OY = = * 106 [N/M]
−K 2 y K2y −6 6

III. Calcul des modes et des fréquences normalisées dans les deux
directions

1. Fréquences de vibrations/Valeurs propres

a. Suivant OX
det | [K]OX - λ[M] |= 0

[M]-1 det | [K]OX - λ[M] |= [M]-1 * 0 en prémultipliant par [M]-1

Ainsi, det |[M]-1 [K]OX - λ[I]|= 0

Avec :

8
 −1 1 d −b
[ M] = *
dét [M ] −c a

−1
[ M] =
1
11 250
2 *
11 250
0
0
[
112 50 ]
−1 1 1 0
[ M] = * matrice de masse
11 250 0 1
inversée
Par ailleurs,
6
−1 10 1 0 6.51 −2.67
[ M ] ∗[K ] OX = * *
1 1250 0 1 −2.67 2.67

6
10 6.52 −2.67
¿ ∗
11 250 −2.67 2.67

Donc, det | [ M ] * [ K ]OX – λ[I] | = 0

−1

6
10 6.51−λ −2.67
⇨ det | ∗ ∨¿ 0
11 250 −2.67 2.67− λ

6
10
⇨ 2
∗( 6.51−λ ) ( 2.67−λ )−2.67 = 0
11 250

λ – 9.18λ + 10.2528 = 0
2

△ = 9.18 2 – 4 * 10.2528

△ = 43.2612 > 0
Nous avons alors une solution double ;
−b ± √ ∆
λ 1, 2
2a

9.18−√ 43.2612
 λ 1= =1.30
2
9.18+√ 43.2612
 λ 2= =7.88
2

9
 2
D’où : ω 1 x =λ 1=1.30
2
ω 2 x =λ 2=7.88
Nous avons ainsi, deux types/fréquences de vibrations.

b. Suivant OY

det | [K]Oy - λ[M] |= 0

⇨ det |[M]-1 [K]Oy - λ[I] |= 0 en prémultipliant par [M]-1

Avec,
6
10 1 0 390 6 λ 0
det | ∗ ∗ − ∨¿ 0
11 250 0 1 −6 6 0 λ

6
10 390−λ −6
det | ∗ ∨¿ 0
11 250 −6 6−λ

Cela entraine que :

6
10 2
∗( 390−λ ) ( 6−λ ) −6 = 0
11 250
(390 – λ)(6 – λ) - 6 2 = 0

λ – 396λ + 2340 = 0
2

△ = 147956 > 0
Nous avons alors une solution double ;
−b ± √ ∆
λ 1, 2
2a
Ainsi,
396−√ 147456
λ 1= =6
2
396 +√ 147456
λ 2= =390
2

10
 2
D’où , ω 1 y =λ1=6
2
ω 2 y =λ2 =390

2. Modes de vibrations/vecteurs propres dans les deux directions

Les modes de vibrations sont donnés par les valeurs propres.

a. Suivant OX

X = 0⃗
| [K]OX - λi[M] | ⃗ i

X = 0⃗
| [M]-1 [K]OX - λi [I] | ⃗ i

Pour : λi = λ1x

6
10 6.51 −2.67 1.30 0 ⃗ ( 1) ⃗
| ∗ − ∨X =0
11 250 −2.67 2.67 0 1.30

5 77.37 −2 37.33 x 11
⇔ ∗ =0
−2 37.33 2 36.03 x 12

5 77.37 x 11−2 37.33 x12=0(1)

⇨
−2 37.33 x 11 +2 36.03=0 (2)

(1) Donne : x 11 = 2.43 x 12

(2) Donne : x 11 = 0.99 x 12 ⩯ 1 x 12

1 1
Donc, ⃗
X = x 11 (
( 1)
) ou ⃗
X = x 12 (1)
( 1)
2.43

Pour : λi = λ2x = 7.88

X = 0⃗
| [M]-1 [K]OX – λ2 [X] | ⃗ i

6
10 6.51 −2.67 7.88 0 ⃗ ( 2)
⇔| ∗ − ∨ X =0
11 250 −2.67 2.67 0 7.88

11
 5 70.77 x 21−237.33 x 22=0(1)
⇨
−2 37.33 x 21+2 29.45 x 22=0(2)

(1) Donne : x 22 = 2.4 x 21

(2) Donne : x 21 = 0.99 x 22 ⩯ 1 x 22

1 1
Donc, ⃗
X = x 21 (
( 2)
) ou ⃗
X = x 22 (1)
( 2)
2.4
Nous pouvons représenter ses modes comme suit :

 ⃗
X (1) amplitudes = fonct (ω)

ω ω ou ω
suivant x 11 suivant x 12

 ⃗
X (2) amplitudes = fonct (ω)

ω ou ω
suivant x 21 suivant x 22

Nous reteindrons le cas de vibration en phase (2) soit :

1 1
⃗
X (1) = x 12 (1) et ⃗
X (2) = x 22 (1)

Deux modes de vibrations suivant OX.

b. Suivant OY

| [K]Oy – λi [M] | ⃗y i = 0⃗

| [M]-1 [K]Oy – λi [I] | ⃗y i = 0⃗

Pour : λ1y= w21x = 6

12
 6
10 390 −6 6 0 ⃗ (1) ⃗
| ∗ − ∨Y =0
11 250 −6 6 0 6

3 4660.67 y 11−533.33 y 12=0(1)

⇨
−533.33 y 11 +527.33 y 12=0(2)

(1) Donne : y 12 = 65 y 11
(2) Donne : y 11 = 0.99 y 12 ⩯ 1 y 12

1 1
⃗ (1 ) = y 11 (
Donc Y ⃗ (1 )= y 12 ( )
) ou Y
65 1
Pour : λ2y= w22y = 390

6
10 390 −6 390 0 ⃗ (2 ) ⃗
| ∗ − ∨Y =0
11 250 −6 6 0 390

3 4276.67 y 21−533.33 y22=0(1)

⇨
−533.33 y 21+143.33 y 22=0(2)

(1) Donne : y 22 = 64.27 y 21

(2) Donne : y 21 = 0.3 y 22

1 0. 3
Donc, ⃗
Y = y 21 (
(2 )
) ou ⃗
Y
(2 )
= y 22 ( )
64.2 7 1
Nous pouvons représenter ses modes comme suit :

 ⃗ (1 )
Y amplitudes = fonct (ω)

ω ou ω

suivant y 11 suivant y 12

 ⃗ ( 2)
Y amplitudes = fonct (ω)

ω ou ω

13
 suivant y 21 suivant y 22

Nous reteindrons le cas de vibration en phase (2) soit :

1 0. 3
⃗ ⃗ (2 ) = y 22 ( )
Y = y 12 (1) et Y
(1 )
1

Deux modes de vibrations suivant OY.

VI. Détermination des participations modales (suivant OX et OY)

Remarque :

Source : Dynamique des structures, 5.1.1 Oscillations non-amorties, page 49

En effet, d’après le processus de décomposition modale, un système à n = 2 degrés de liberté

peut être traité comme un oscillateur simple.

= ou z(t) est la coordonnée modale.

Connaissant la valeur de l’accélération du sol, l’instant à partir duquel on désir connaitre la

réponse et en nous imposant des conditions initiales, nous serions à même de déterminer à
partir de l’intégrale de Duhamel (car chargement quelconque), la réponse de chaque mode de
la façon suivante :
t
ẋ (0) 1
x (t )= sin (ω¿¿ n t)+ x ( 0 ) cos ( ωn t ) − ∫ ẍ g ( t ) sin ( t−τ ) dτ ¿ avec F(t) = -m ẍ g (t )
ωn ωn 0

La réponse finale étant la somme des réponses des différents modes.

Cependant, dans le cadre de notre étude, la condition d’orthogonalité des modes n’est pas
respectée. Et cela est certainement dû à la discontinuité du voile/poteau (inertie variable).

14
En effet :

[ 11 11][ m0 m0 ] [11 11]=[ 11 11] [mm mm]=m [ 22 22] n est pas diagonale
'

Soit A la matrice des vecteurs modaux, Zn les coordonnées modales s’exprimant à partir du
changement de variable.

2
x a a Z
X = 1= 11 12= 1= AZ=∑ an Z
x 2 a21 a22 Z 2 n=1
n

On a : F(t) = -m ẍ g (t)

En tenant compte de ce changement de variable et en prémultipliant par AT, on obtient l’équation

T T T T
généralisée suivante : A MA Z̈ + A CA Ż + A KAZ =− A M ẍ g (t ) e x

¿ ¿ ¿ T
Soit , M Z̈+C Ż+ K Z=− A M e x ẍ g ( t )=−r ẍ g (t)

T
Ou r =A M e x est le facteur de participation

a11 a12 m1 0 1 r 1
r= ∗ ∗ =
a21 a 22 0 m2 1 r 2
Une grandeur permet de quantifier à priori la contribution de chaque mode à la réponse totale. Il
s’agit de la masse modale qui est définie comme suit :
rn
mmod , n=( ¿ )2 m¿n
m n

Si la condition d’orthogonalité des modes était respectée nous aurions :

mmod , n=¿ ¿
¿ ¿ ¿
(Car la matrice M , C , K est diagonale).

a. Suivant OX

1 1
Nous avons retenu : ⃗
X = x 12 ( ) o ⃗
( 1)
X = x 22 ( )
( 2)
1 1

1 1
A=
1 1

¿ T 2 2
M ox = A MA=m
2 2

Avec : m1 = m2 = m = 11250 kg

15
 r ox = [ 11 11][ m0 m0 ][ 11]
1

2
, avec m 1=m2=m

r ox = [ 11 11][ mm ]=m [22]=[rr ] , 1

2
avec m = 11250 kg

On a :

2 2
rn rn ¿
mmod , n= ¿ 2 ∗m =
mn ¿∨m¿k ∨¿
n

||m1||=√ ( 2 m)2 +(2 m)²=2 √ 2∗m

||m1||=2 √2∗1 1250
2 2
r1 2 2500
 m mod ,1 = = =15909 . 90
∥ m¿1 ∥ 2 √ 2∗1 1 250
2
r2
 m mod , 2= =15909.90
∥ m¿2 ∥

| ¿| | ¿|
Car r 1=r 2 et |m1| = |m 2|

Ainsi suivant OX la contribution de chaque mode est la même

b. Suivant OY

1 0.2
Nous avons retenu : ⃗
Y (1 ) = Y 12 ( ) ou ⃗
Y ( 2) = Y 22 ( )
1 1

1 0.2 m1 0 1
r oy =
1 1 0 m2 1

1 0.2 m 0 1
r oy =
1 1 0 m1

1 0.2 m
r oy =
1 1 m

1.2 1.2 r
r oy =m =m = 1
2 2 r2

On a :

16
  Pour le voile :

2
r1 (1.2∗1 1250) ²
mmod ,1 = = =5727.56
¿∨m1∨¿ 2 √ 2∗1 1 250
¿

 Pour le poteau :

2
r2 (2∗11 250) ²
mmod , 2= = =15909.90
¿∨m2∨¿ 2 √2∗1 1250
¿

Suivant OY le mode 1 a une participation plus petite que le mode 2.

V. Représentation de la décomposition dans l’espace modale :

suivant OX suivant OY

17
 VI. Conclusion
Nous constatons que la déformation du bâ timent autour de l’axe OX est la même sous l’effet de
cette excitation sismique, ce qui se justifie d’une part par la représentation des amplitudes en
fonction des pulsations A = f (ω ) et d’autre part par la valeur des participations modales suivant
cette direction.
Par contre, autour de l’axe OY la déformation du voile est très petite devant celle du poteau. Ce
qui est clairement justifié par la représentation de A = f (ω ) mais plus précisément par les
participations modales suivant OY. En effet, m mod ,1 < mmod , 2.

Pour un changement d’inertie importante, la structure change difficilement de mode quelque soit
les directions.

18
 Annexe mathématique

Algèbre matricielle

 Notation matricielle

Pour un ensemble d’équations linéaires :

5x1 – 4x2 + x3 = 0
-4x1 + 6x2 -4x3 +x4 =1
X1 – 4x2 + 6x3 -4x4 =0
X2 – 4x3 + 5x4 =0

La notation matricielle de l’ensemble des équations, soit Ax = b, est :

 Matrice identité - (diagonale de 1)

19
 1 0 0
I =diag ( 1 ) , par exemple I = 0 1 0
0 0 1

 Opérations matricielles

o Addition, A+ B = C

A
[0.52 1 1
3 0
;B
] [
3 1 2
2 4 1
; C= A + B=
]
5 2 3
2.5 7 1 [ ]
o Multiplication par un scalaire ; kA = C

A
[0.52 1 1
3 0 ]
; k =2; C=kA=
4 2 2
1 6 0 [ ]
o Multiplication, AB = C
(A est une matrice p x m, B une matrice m x q et C une matrice p x q)
m
C ij =∑ air brj
r=1

[ ][ ]
5 3 1 1 5
Par exemple, si A= 4 6 2 ;B 2 4
10 3 4 3 2

Alors :

C 11=( 5 )( 1 ) + ( 3 ) ( 2 ) + ( 1 ) ( 3 )=14

C 21=( 4 )( 1 ) + ( 6 )( 2 )+ (2 )( 3 ) =22

C 31=( 10 ) (1 )+ ( 3 )( 2 ) + ( 4 ) ( 3 )=28 etc .

Donc :

[ ]
14 39
C= 22 48
28 70

o La multiplication matricielle n’est pas commutative

o A B≠ BA

[ 52 64 ] ; B=[ 18 93]
A=

AB=
[52 64 ][ 18 93]=[ 5334 6330]
20
 BA=
[ ][ ] [ ]
1 7 5 6
8 3 2 4
=
19 34
46 60

[ 5334 6330] ≠ [ 1946 3460 ]

 Transposition d’une matrice, AT ou échanger les lignes et les colonnes de la matrice initiale.
Par exemple :

[ ] [ ]
3 1 4 3 2 7
T
A= 2 5 3 ; A = 1 5 9
7 9 8 4 3 8

o Evaluer le produit V T AV

[ ] []
3 2 1 1
A= 2 4 2 ; V = 2
1 2 6 −1

[ ][ ] [ ]
3 2 1 1 6
X =A V = 2 4 2 2 = 8
1 2 6 −1 −1

[]
6
V A V =V X=[ 1 2 −1 ] 8 =23
T T

−1

 Transposition d’une matrice, AT ou échanger les lignes et les colonnes de la matrice

initiale. Par exemple :
3 1 4 3 2 7
T
2 5 3 ; A = 1 5 9
7 9 8 4 3 8
A =

o Evaluer le produit V T A V

21
 o
T

o La transposée d’un produit matriciel est égale au produit des matrices

transposées dans l’ordre inverse, (ABC)T = CTBTAT

 Matrice inverse, A-1

o Pour des produits matriciels

GAB = I où G = (AB)-1

Donc: (AB)-1 AB = I

 Determinant d’une matrice

n
Det A= ∑ (−1)1+ j a1 j det A 1 j
j=1

A=
[ a 11 a12
a 21 a22 ]
2 3
Det A=(−1) a11 det A 11 +(−1) a12 det A12 ou

det A 11=a11 et det A 12=a12

22
Alors: det A = a 11 a22 - a 12 a21

o Exemple de déterminant d’une matrice 2*2

A=
[−12 32]
det A = a 11 a22−a12 a21 =(−1 ) ( 3 )− ( 2 )( 2 ) =−7

o Exemple de déterminant d’une matrice 3*3

[ ]
2 1 0
A= 1 3 1
0 1 2

Det

2 3 4
A=(−1) a11 det A 11 +(−1) a12 det A12 +(−1) a 13 det A 13
Det

2
A=(−1 ) ( 2 ) det
[ 31 12] +(−1) ( 1) det [ 10 12]+(−1) (0)det [ 10 31]
3 4

det A = (2)(3)(2)-(1)(1)-(1)(2)-(0)(1) + 0 = 8

 Problèmes aux valeurs propres

Av = v
Av =i i vi

Il y a n valeurs propres correspondant aux n vecteurs propres

A I)v = 0
Cette équation a une solution non triviale seulement si

det(A- I) = 0
Les valeurs propres sont les racines de l’équation p( )
det(A- I)

o Exemple de valeurs propres et vecteurs propres

23
Il y a aussi deux vecteurs propres correspondants

24
 Bibliographie

[1]. Mohamed CHENNIT, « Dynamique des Structures cours et exercices », mémoire de Master,
Université des Sciences et de la technologie Houari Boumediene, Algérie, 2022

[2]. I. Smith & P. Lestuzzi, « Dynamique des Structures », mémoire de Master, Ecole Polytechnique
Fédérale de Lausanne, Russie, 2013/2014

25