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................................................................................ 5
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 .............................................. 10
VI. Détermination des participations modales (suivant OX et OY) ........................................................ 13
V. Représentation de la décomposition dans l’espace modale :........................................................ 16
VI. Conclusion .................................................................................................................................. 17
Annexe mathématique ........................................................................................................................... 18
B.Algèbre matricielle .......................................................................................................................... 19
Bibliographie........................................................................................................................................... 23

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 GP

‫ = م‬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

u1(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)-c1𝑢̇ 1-c2(𝑢1̇-𝑢2̇)

M1ẍ1 + (c1+c2) 𝑢1̇ – c2𝑢2̇ + (k1 + k2) u1 – k2u2 = 0

Or, u1 = x1 – xg

⇨ x1 = u1 + xg
Alors, M1 𝒖𝟏̈ + (c1+c2) 𝒖𝟏̇ – c2𝒖𝟐̇ + (k1 + k2) u1 – k2u2 = - M1ẍg (1)

• 2 ème équation :

M2ẍ2 = -c2(𝑢1̇-𝑢2̇) – k2 (u2 – u1)

M2ẍ2 –c1𝑢1̇ + c1𝑢2̇ – k2u1 + k2u2 = 0

Or, u2 = x2 – xg

⇨ x2 = u2 + xg
Alors, M2 𝒖𝟐̈ – c1 𝒖𝟏̇ + c2𝒖𝟐̇ – k2u1 + k2u2 = - M2 ẍg (2)

• Sous forme matricielle nous avons :

4
𝑀1 0 𝑢1̈ + 𝑐1 + 𝑐2 −𝑐2 𝑢1̇ + 𝑘1 + 𝑘2 −𝑘2 𝑢1
= −M1 𝑢1
̈
0 𝑀2 𝑢2̈ −𝑐1 𝑐2 𝑢2̇ −𝑘2 𝑘2 𝑢2 −M2 𝑢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 * ‫ = م‬12250 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 = 12250 kg

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

5
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 𝐸𝐼
𝐾=
ℎ3
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

• Au RDC :
2 ∗ 0.23 1
𝐼1𝑥 = = 𝑚4
12 750

• Au R + 1 :

0.3 ∗ 0.23
𝐼2𝑥 = = 2 ∗ 10−4 𝑚4
12
b. Suivant OY

6
 • Au RDC :

0.2 ∗ 23 2
𝐼1𝑦 = = 𝑚4
12 15
• Au R + 1 :

0.2 ∗ 0.33
𝐼2𝑦 = = 4.5 ∗ 10−4 𝑚4
12

5.2 Rigidités des niveaux 1 et 2

a. Suivant OX

• Au RDC :

12𝐸 ∗ 𝐼1𝑥
𝐾1𝑥 =
ℎ1 3

12 ∗ 30 ∗ 109 ∗ 1
𝐾1𝑥 = = 3.84 ∗ 106 𝑁/𝑚
53 ∗ 750

• AR+1:

12𝐸 ∗ 𝐼2𝑥
𝐾2𝑥 =
ℎ2 3
12 ∗ 30 ∗ 109 ∗ 2 ∗ 10−4
𝐾2𝑥 = = 2.67 ∗ 106 𝑁/𝑚
33

b. Suivant OY

• Au RDC :

12𝐸 ∗ 𝐼1𝑦
𝐾1𝑦 =
ℎ1 3

12 ∗ 30 ∗ 109 ∗ 2
𝐾1𝑦 = = 3.84 ∗ 108 𝑁/𝑚
53 ∗ 15

• Au R + 1 :

7
 12𝐸 ∗ 𝐼2𝑦
𝐾2𝑦 =
ℎ2 3

12 ∗ 30 ∗ 109 ∗ 45 ∗ 10−5
𝐾2𝑦 = = 6 ∗ 106 𝑁/𝑚
33

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

Matrice de masse :

12250 0
[M] =
0 12250

Matrice de rigidité :

𝐾1𝑥 + 𝐾2𝑥 −𝐾2𝑥 6.51 −2.67

[K]OX = = * 106 [N/M]
−𝐾2𝑥 𝐾2𝑥 −2.67 2.67

𝐾1𝑦 + 𝐾2𝑦 −𝐾2𝑦 390 −6

[K]OY = = * 106 [N/M]
−𝐾2𝑦 𝐾2𝑦 −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 :

1 𝑑 −𝑏
[M]−1 = 𝑑é𝑡[M] *
−𝑐 𝑎

1 12250 0
[M]−1 = * [ ]
122502 0 1250

8
 1 1 0
[M]−1 = 12250 * matrice de masse inversée
0 1
Par ailleurs,
106 1 0 6.51 −2.67
[M]−1 ∗ [K]𝑂𝑋 = 12250 * *
0 1 −2.67 2.67

106 6.52 −2.67

= ∗
12250 −2.67 2.67

Donc, det | [M]−1 * [K]𝑂𝑋 – λ[I] | = 0

6
6.51 − 𝜆 −2.67
⇨ det |12250
10
∗ | =0
−2.67 2.67 − 𝜆

6
⇨ 12250
10
∗ (6.51 − 𝜆)(2.67 − 𝜆) − 2.672 = 0

𝜆2 – 9.18λ + 10.2528 = 0

△ = 9.182 – 4 * 10.2528
△ = 43.2612 > 0
Nous avons alors une solution double ;

−𝑏 ± √∆
𝜆1,2
2𝑎

9.18− √43.2612
✓ 𝜆1 = 2
= 1.30
9.18+ √43.2612
✓ 𝜆2 = 2
= 7.88

D’où :
2
𝜔1𝑥 = 𝜆1 = 1.30
2
𝜔2𝑥 = 𝜆2 = 7.88

Nous avons ainsi, deux types/fréquences de vibrations.

b. Suivant OY

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

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

Avec,

106 1 0 390 6 𝜆 0
det | ∗ ∗ − |=0
12250 0 1 −6 6 0 𝜆

106 390 − 𝜆 −6
det | ∗ |=0
12250 −6 6−𝜆

Cela entraine que :

106
∗ (390 − 𝜆)(6 − 𝜆) − 62 = 0
12250

(390 – λ)(6 – λ) - 62 = 0

𝜆2 – 396λ + 2340 = 0
△ = 147956 > 0
Nous avons alors une solution double ;

−𝑏 ± √∆
𝜆1,2
2𝑎
Ainsi,

396 − √147456
𝜆1 = =6
2
396 + √147456
𝜆2 = = 390
2

D’où,
2
𝜔1𝑦 = 𝜆1 = 6
2
𝜔2𝑦 = 𝜆2 = 390

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

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

10
 a. Suivant OX

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

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

Pour : λi = λ1x

106 6.51 −2.67 1.30 0

| ∗ − |𝑋⃗ (1) = ⃗0⃗
12250 −2.67 2.67 0 1.30

530.13 −217.96 𝑥11

⇔ ∗ 𝑥 =0
−217.96 216.66 12

530.13𝑥11 − 217.96𝑥12 = 0 (1)

⇨
−217.96𝑥11 + 216.66𝑥12 = 0 (2)
(1) Donne : 𝑥11 = 2.43𝑥12
(2) Donne : 𝑥11 = 0.99𝑥12 = 1𝑥12

1 1
Donc, 𝑋⃗ (1) = 𝑥11 ( ) ou 𝑋⃗ (1) = 𝑥12 ( )
2.43 1

Pour : λi = λ2x = 7.88

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

106 6.51 −2.67 7.88 0

⇔ | 12250 ∗ − |𝑋⃗ (2) = 0
−2.67 2.67 0 7.88

523.55𝑥21 − 217.96𝑥22 = 0 (1)

⇨
−217.96𝑥21 + 210.08𝑥22 = 0 (2)

(1) Donne : 𝑥22 = 2.4𝑥21

(2) Donne : 𝑥21 = 0.99𝑥22 = 1𝑥22

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

• 𝑋⃗ (1) amplitudes = fonct (ω)

11
 ω ω ou ω
suivant 𝑥11 suivant 𝑥12

• 𝑋⃗ (2) amplitudes = fonct (ω)

ω ou ω
suivant 𝑥21 suivant 𝑥22

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

1 1
𝑋⃗ (1) = 𝑥12 ( ) et 𝑋⃗ (2) = 𝑥22 ( )
1 1

Deux modes de vibrations suivant OX.

b. Suivant OY

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

| [M]-1 [K]Oy – λi [I] | 𝑦⃗ 𝑖 ⃗⃗

=0

Pour : λ1y= w21x = 6

106 390 −6 6 0 ⃗⃗ (1)

| ∗ − |𝑌 = ⃗0⃗
12250 −6 6 0 6

31860.73𝑦11 − 489.79𝑦12 = 0 (1)

⇨ −489.79𝑦11 + 483.79𝑦12 = 0 (2)

(1) Donne : 𝑦12 = 65𝑦11

(2) Donne : 𝑦11 = 0.99𝑦12 = 1𝑦12

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

12
 106 390 −6 390 0 ⃗⃗ (2) ⃗⃗
| ∗ − |𝑌 = 0
12250 −6 6 0 390

31446.73𝑦21 − 489.79𝑦22 = 0 (1)

⇨ −489.79𝑦21 + 99.79𝑦22 = 0 (2)

(1) Donne : 𝑦22 = 64.2𝑦21

(2) Donne : 𝑦21 = 0.2𝑦22 = 1𝑦22

⃗⃗ (2) = 𝑦21 ( 1 ) ou 𝑌
Donc, 𝑌 ⃗⃗ (2) = 𝑦22 (0.2)
64.2 1
Nous pouvons représenter ses modes comme suit :

• ⃗⃗ (1) amplitudes = fonct (ω)

𝑌

ω ou ω

suivant 𝑦11 suivant 𝑦12

• ⃗⃗ (2) amplitudes = fonct (ω)

𝑌

ω ou ω

suivant 𝑦21 suivant 𝑦22

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

⃗⃗ (1) = 𝑦12 (1) et 𝑌

𝑌 ⃗⃗ (2) = 𝑦22 (0.2)
1 1

Deux modes de vibrations suivant OY.

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

13
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 :
𝑥̇ (0) 1 𝑡
𝑥(𝑡) = 𝑠𝑖𝑛(𝜔𝑛 𝑡) + 𝑥(0) cos(𝜔𝑛 𝑡) − ∫ 𝑥̈ (𝑡) sin(𝑡 − 𝜏) 𝑑𝜏 avec F(t) = -m𝑥̈𝑔 (𝑡)
𝜔𝑛 𝜔𝑛 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).
En effet :
1 1 𝑚 0 1 1 1 1 𝑚 𝑚 2 2 ′
[ ][ ][ ]= [ ][ ] = 𝑚[ ] n est pas diagonale
1 1 0 𝑚 1 1 1 1 𝑚 𝑚 2 2

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

𝑥1 𝑎11 𝑎12 𝑍1 2
𝑋= 𝑥 = 𝑎
2 21 𝑎22 = 𝑍2 = 𝐴𝑍 = ∑𝑛=1 𝑎𝑛 𝑍𝑛

On a : F(t) = - 𝑚 𝑥̈𝑔 (t)

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

généralisée suivante :

14
 𝐴𝑇 𝑀𝐴𝑍̈ + 𝐴𝑇 𝐶𝐴𝑍̇ + 𝐴𝑇 𝐾𝐴𝑍 = −𝐴𝑇 𝑀 𝑥̈𝑔(𝑡) 𝑒𝑥

𝑺𝒐𝒊𝒕, 𝑀∗ 𝑍̈ + 𝐶 ∗ 𝑍̇ + 𝐾 ∗ 𝑍 = −𝐴𝑇 𝑀𝑒𝑥 𝑥̈𝑔 (𝑡) = −𝑟𝑥̈𝑔 (𝑡)

Ou 𝑟 = 𝐴𝑇 𝑀𝑒𝑥 est le facteur de participation

𝑎11 𝑎12 𝑚1 0 1 𝑟1
𝑟 = 𝑎 𝑎22 ∗ 0 ∗ = 𝑟
21 𝑚2 1 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 :
𝑟
𝑚𝑚𝑜𝑑,𝑛 = ( 𝑚𝑛∗ )2 𝑚𝑛∗
𝑛

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

(∑𝑛𝑚=1 𝑎𝑚𝑛 𝑚𝑚 )²
𝑚𝑚𝑜𝑑,𝑛 =
∑𝑛𝑛=1 𝑎𝑚𝑛
2 𝑚
𝑛

(Car la matrice 𝑴∗ , 𝑪∗ , 𝑲∗ est diagonale).

a. Suivant OX

1 1
Nous avons retenu : 𝑋⃗ (1) = 𝑥12 ( ) o 𝑋⃗ (2) = 𝑥22 ( )
1 1

1 1
𝐴=
1 1

∗ 2 2
𝑀𝑜𝑥 = 𝐴𝑇 𝑀𝐴 = 𝑚
2 2

Avec : m1 = m2 = m = 12250 kg

1 1 𝑚1 0 1
𝑟𝑜𝑥 = [ ][ ] [ ] , avec 𝑚1 = 𝑚2 = 𝑚
1 1 0 𝑚2 1

1 1 𝑚 2 𝑟1
𝑟𝑜𝑥 = [ ] [ ] = 𝑚 [ ] = [𝑟 ], avec m = 12250 kg
1 1 𝑚 2 2

On a :

𝑟𝑛2 ∗
𝑟𝑛2
𝑚𝑚𝑜𝑑,𝑛 = ∗ 𝑚 𝑛 =
𝑚𝑛∗2 ||𝑚𝑘∗ ||

||𝑚1 || = √(2𝑚)2 + (2𝑚)² = 2√2 ∗ 𝑚

||𝑚1 || = 2√2 ∗ 12250

𝑟2 245002
❖ 𝑚𝑚𝑜𝑑,1 = ∥𝑚1∗ ∥ = 2 = 17324.12
1 √2∗12250

15
 𝑟2
❖ 𝑚𝑚𝑜𝑑,2 = ∥𝑚2∗ ∥ = 17324.12
2

Car 𝑟1 = 𝑟2 𝑒𝑡 ||𝑚1∗ || = ||𝑚2∗ ||

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

b. Suivant OY

⃗⃗ (1) = 𝑌12 (1) ou 𝑌

Nous avons retenu : 𝑌 ⃗⃗ (2) = 𝑌22 (0.2)
1 1

1 0.2 𝑚1 0 1
𝑟𝑜𝑦 =
1 1 0 𝑚2 1

1 0.2 𝑚 0 1
𝑟𝑜𝑦 =
1 1 0 𝑚 1

1 0.2 𝑚
𝑟𝑜𝑦 =
1 1 𝑚

1.2 1.2 𝑟1
𝑟𝑜𝑦 = 𝑚 =𝑚 = 𝑟
2 2 2

On a :

➢ Pour le voile :

𝑟12 (1.2 ∗ 12250)²

𝑚𝑚𝑜𝑑,1 = ∗ = = 6236.68
||𝑚1 || 2√2 ∗ 12250

➢ Pour le poteau :

𝑟22 (2 ∗ 12250)²
𝑚𝑚𝑜𝑑,2 = ∗ = = 17324.12
||𝑚2 || 2√2 ∗ 12250

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 :

16
 suivant OX suivant OY

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

17
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, 𝒎𝒎𝒐𝒅,𝟏 < 𝒎𝒎𝒐𝒅,𝟐 .

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

Annexe mathématique

18
 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)

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

▪ Opérations matricielles

o Addition, A+ B = C

2 1 1 3 1 2 5 2 3
𝐴[ ]; 𝐵[ ];𝐶 = 𝐴 + 𝐵 = [ ]
0.5 3 0 2 4 1 2.5 7 1

o Multiplication par un scalaire ; kA = C

2 1 1 4 2 2
𝐴[ ] ; 𝑘 = 2; 𝐶 = 𝑘𝐴 = [ ]
0.5 3 0 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)
𝐶𝑖𝑗 = ∑𝑚
𝑟=1 𝑎𝑖𝑟 𝑏𝑟𝑗

5 3 1 1 5
Par exemple, si 𝐴 = [ 4 6 2 ] ; 𝐵 [2 4]
10 3 4 3 2

19
Alors :

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

𝐶21 = (4)(1) + (6)(2) + (2)(3) = 22
𝐶31 = (10)(1) + (3)(2) + (4)(3) = 28 𝑒𝑡𝑐.
Donc :
14 39
𝐶 = [22 48]
28 70

o La multiplication matricielle n’est pas commutative

o 𝐴𝐵 ≠ BA
5 6 1 9
𝐴=[ ]; 𝐵 = [ ]
2 4 8 3
5 6 1 9 53 63
𝐴𝐵 = [ ][ ]=[ ]
2 4 8 3 34 30
1 7 5 6 19 34
𝐵𝐴 = [ ][ ]=[ ]
8 3 2 4 46 60
53 63 19 34
[ ]≠[ ]
34 30 46 60

▪ Transposition d’une matrice, 𝐴𝑇 ou échanger les lignes et les colonnes de la matrice initiale.
Par exemple :
3 1 4 3 2 7
𝐴 = [2 5 3] ; 𝐴𝑇 = [1 5 9]
7 9 8 4 3 8

o Evaluer le produit 𝑉 𝑇 𝐴𝑉

3 2 1 1
𝐴 = [ 2 4 2] ; 𝑉 = [ 2 ]
1 2 6 −1
3 2 1 1 6
𝑋 = 𝐴𝑉 = [2 4 2] [ 2 ] = [ 8 ]
1 2 6 −1 −1
6
𝑉 𝑇 𝐴𝑉 = 𝑉 𝑇 𝑋 = [1 2 −1] [ 8 ] = 23
−1

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

initiale. Par exemple :

20
 A =

o Evaluer le produit 𝑉 𝑇 𝐴𝑉

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

Det 𝐴 = ∑𝑛𝑗=1(−1)1+𝑗 𝑎1𝑗 𝑑𝑒𝑡𝐴1𝑗

𝑎11 𝑎12
𝐴 = [𝑎 𝑎22 ]
21

Det𝐴 = (−1)2 𝑎11 𝑑𝑒𝑡𝐴11 + (−1)3 𝑎12 𝑑𝑒𝑡𝐴12 ou

𝑑𝑒𝑡𝐴11 = 𝑎11 et 𝑑𝑒𝑡𝐴12 = 𝑎12

Alors: det𝐴 = 𝑎11 𝑎22 - 𝑎12 𝑎21

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

21
 −1 2
𝐴=[ ]
2 3
det𝐴 = 𝑎11 𝑎22 − 𝑎12 𝑎21 = (−1)(3) − (2)(2) = −7

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

2 1 0
𝐴 = [1 3 1]
0 1 2
Det𝐴 = (−1)2 𝑎11 𝑑𝑒𝑡𝐴11 + (−1)3 𝑎12 𝑑𝑒𝑡𝐴12 +

(−1)4 𝑎13 𝑑𝑒𝑡𝐴13

3 1
Det𝐴 = (−1)2 (2)𝑑𝑒𝑡 [ ]+
1 2
1 1 1 3
(−1)3 (1)𝑑𝑒𝑡 [ ] + (−1)4 (0)𝑑𝑒𝑡 [ ]
0 2 0 1
det 𝐴 = (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

Il y a aussi deux
vecteurs
propres correspondants

22
Bibliographie

23
[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

24