CALCUL DYNAMIQUE DES

STRUCTURES

⚫ CALCUL DYNAMIQUE
APPLIQUE AUX STRUCTURES
DE BETON ARME SOUMISE A
L’ACTION SISMIQUE

1
Azzouz Bouguerba-:[Link]- 22/02/2023
 INTRODUCTION

⚫ Calcul dynamique des structures

⚫ Connaissance du calcul isostatique et
hyperstatique notamment calcul des
déformations
⚫ S’adresse aux ingénieurs pour les bases de
calcul dynamique des structures

2
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
PROGRAMME DE FORMATION

⚫ Rappels fondamentaux de la RDM

⚫ Eléments de la rigidité
⚫ Construction de la Matrice de Rigidité
⚫ Notion de dynamique
⚫ Equation de la Dynamique
⚫ Exemple de calcul

3
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE RAIDEUR(4)
⚫ ) Barre biarticulée :
48EI
k= 3 L(EI)
L
⚫ ) Barre encastrée aux 2 extrémités
⚫ Déplacement unitaire à l’origine

FB
d=1
Ma Mb
A B
EI
FA h
4
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE RAIDEUR (6)
⚫ Barre encastrée - articulée soumise à un déplacement
unitaire
Ma=0 F = k *d
FB
d=1
Mb
d =1
EI
A B 3 EI
h k =
FA h3

5
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE RAIDEUR (7)
⚫ Les portiques : Cas articulation
d=1 d=1

I
I1 = I 2 =
2
d =1
(h,E,I/2) (h,E,I/2)
3 E ( I1 + I 2 )
k =
h3
3 EI
k =
h3

6
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE RAIDEUR (7)
⚫ Les portiques : Cas encastré

d=1 d=1
I
I1 = I 2 =
2
d =1
(h,E,I/2) (h,E,I/2) 12 E ( I1 + I 2 )
k =
h3
12 EI
k =
h3

7
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE RAIDEUR (8)
⚫ LA RAIDEUR EST :

K=
3
h
⚫α :Coefficient dépendant du
mode de liaison de la barre
8
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE RAIDEUR (9)
LA RAIDEUR (k) DEPEND :
1. Des liaisons

2. De l’Inertie des sections

3. Du Module E → Matériau

L
4. De la Longueur
9
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 10
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE DYNAMIQUE(1)
⚫ Les éléments de structure,pour son analyse
dynamique,sont :
⚫ Sa masse
⚫ Son élasticité ou sa raideur
⚫ Le matériau constituant cette structure

⚫ Sous l’action sismique , il y a apparition :

⚫ Des forces d’inertie des masses
⚫ Des forces de rappel par les raideurs
⚫ Des forces d’amortissement par la capacité,du
matériau utilisé, à stocker puis à dissiper
l’énergie.
11
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE DYNAMIQUE(2)
⚫ Système à un seul degré de liberté :1ddl
⚫ 1)-Formulation des forces
x
fr
Ressort/Raideur :k

k fr
x

F(t)−mx=0
2
F(t)−m d 2x =0
⚫ 2) Force d’inertie dt

fi =mx
accélération

fi m
Solide F(t)
d 2x
indéformable
x x= =accélération
dt 2 12
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
ELEMENTS DE DYNAMIQUE(3)
⚫ 3) Force d’amortissement
c fa =cx

fa
m  = dx =vitesse
x
Solide F(t) dt
indéformable
x

c
m fi
F(t)
m
k Solide fa
indéformable F(t)
x fr
x

13
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE DYNAMIQUE(4)
⚫ Formulation de l’équation
fi + fa + fr = F(t) mx+cx + Kx= F(t)
⚫ 1) Oscillations libres : F(t) = 0
= K
⚫ m
mx+cx+Kx=0 x+ c x

m
 + K x=0
m
Posons

2 = c
m

x+ 2x + x=0 

14
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 ELEMENTS DE DYNAMIQUE(5)
⚫ 2) Oscillations libres non amorties (C=0)
x+  x=0 x(t)= [Link] t + [Link] t
x=− ([Link] t +[Link] t) x= −  x
⚫ Les conditions aux limites sont :Pour t=0 => x=0
0
B= x0 x(0)= x0 A= x

⚫ Soit l’équation du mouvement :
𝑥ሶ
𝑥(𝑡) = sin 𝜔 𝑡 + 𝑥0. cos 𝜔 𝑡
𝜔

15
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 RIGIDITE-MATRICE DE RIGIDITE(1)
⚫ SYSTEME A 1DDL DIAPHRAGME RIGIDE

Δ
Δ x

L U
N U
4 EI 6 EI
E MA =  6 EI
 = M B
N MA = 
h
h2 h2
6 EI R D
 O 6 EI E
h2 B  = M B
2 EI T 2 EI h2 P 6 EI
B  = M B A MB =  L MB = 
h T h A h2
I h C
h O 6 EI E
N TA = A M TA =
12EI

6 EI h2 A E h3
  E
N
h2 N 12 EI
A
− 6 EI  = M A T
4 EI TB = A
12EI
TA = 3   h2 − 12EI
 = M A TB = 
h A h2 h h3
A

16
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 RIGIDITE-MATRICE DE RIGIDITE(2)
⚫ Des rotationsθA en (A) , θB en (B) et une translation ΔA en (A) induisent en (B) un moment de :
4 EI 2 EI 6 EI
MB =  +  + 2 
h h h
⚫ Ainsi une force horizontale appliquée au diaphramme,qui se déplace horizontalement d’une valeur Δ par
rapport à sa position d’équilibre, induit aux poteaux les sollicitations suivantes :

FH
Δ
B C L’équilibre est obtenu par :
12EI 12EI
1 FH − 3
 − 3
=
6 EI 6 EI h h
MB =  MC = 
h2 h2
4 EI 2 EI F x = 0
MB =  + C
h h M = 0
4 EI 2 EI
MA =
6 EI
 MC = C + 
12 EI

h2 h h
h3
Par exemple , l’équilibre des moments en B est :
12 EI
 6 EI 6 EI
A h3 D MD =
6 EI
h2
 MB + MC − 2 − 2  =
Rigidité du portique : h h
24EI La rigidité K du portique est la force obtenue
FH = 3 
h pour un déplacement unitaire (Δ=1) , soit :
24 EI
K ( = FH ) =
* (  = ) 17
Azzouz Bouguerba-:[Link]-Universiapolis
h 3 Agadir 22/02/2023
 RIGIDITE-MATRICE DE RIGIDITE(3)
Structure à deux degrés de liberté (2DDL)
K12 K 22 Δ
=1
Δ=1
I2 I2 h2

Δ=1
K11 K 21

Δ=1
I1 I1 h1

L
Déplacement unitaire au 1er niveau Déplacement unitaire au 2ème niveau

K 11 K 12
La matrice de rigidité est :
K=
K 21 K 22
Azzouz Bouguerba-:[Link]-Universiapolis Agadir
18
22/02/2023
 RIGIDITE-MATRICE DE RIGIDITE(4)
⚫ Pour des forces appliquées H1 et H2 aux 1er et 2ème niveaux donnant des
déplacements horizontaux x1 et x2 aux 1er et 2ème niveaux sont liés selon la forme
matricielle suivante :

H⚫1 K11 K12 x1 H 1 = K 11 .x1 + K 12 .x 2

= F = K.d
H 2 K 21 K 22 x2
H 2 = K 21 .x1 + K 22 .x 2
K12

I 21 I 22 I 21 I 22
h2 x1=1.00 h2 x1=1.00
K11

I11 I12 I11 I12

h1 h1

19
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 RIGIDITE-MATRICE DE RIGIDITE(5)
I 11 = I 12 = I 21 = I 22 = I 48EI
h1 = h2 = h 12EI 22 K 11 = 3
K12 h23 h
Pour le deuxième niveau , l’équilibre est :
12EI 21
h23 12EI 21 12EI 22
K 12 + + =0
I 21 I 22
3 3
h 2 h 2
12 EI 21 12EI 22
h23
24EI
3
h2

K11
Δ=1
K 12 = − 3
12EI 11 12EI 12 h
h 3 h13
1
I11 I12

20
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
RIGIDITE-MATRICE DE RIGIDITE(6)
( Rappel : I21=I22=I et h2=h)
12EI 21 12EI 22
Δ=1 12EI 22 K 22 − − =0
h23 h23
K22 h23
12EI 21 24EI
h23 12EI 22 K 22 = 3
K21
h 3
2 h
12 EI 21
24EI
K 21 = − 3
3
h2 12EI 21 12EI 22
K 21 + 3
+ 3
=0
h h
2 2
h

48EI 24EI
K=
K 11
3
− 3 K 12
K 21 K= h h K 22
24EI 24EI
D’où la matrice de rigidité : − 3 3
h h 21
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 DETERMINATION DU MODE FONDAMENTAL(1)
⚫ L’équation dynamique pour la vibration libre non amortie est :

Mx + Kx = 0 M : matrice des masses

K : matrice de rigidité de la structure
:
x accélération
x : déplacement dans le sens du degré de liberté

L’équation des déplacements est de la forme : x = C sin(t ) C=

a1
a2
x = −C  sin(t ) − C  sin(t ) + CK sin(t ) = 0
Soit : Relation vérifiée si
x = −C 2 x
(−  M + K ) = 0 (D28a) Déterminant de la matrice est nul :

− M + K C = 0 Par exemple pour le système à 2DDl , ci-dessus , on a :

m 0 −   m1 0
M= 1 K=
K 11 K 12 −  M =
0 m2 K 21 K 22 0 −   m2

22
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 DETERMINATION DU MODE FONDAMENTAL(2)
K 11 −   m1 K 12
Le déterminant est : =0
K 21 K 22 −  m 2


Le développement du déterminant donne :

( K11 −   m1 )( K 22 −   m2 ) − K12 * K 21 = 0
La résolution de l’équation permet d’obtenir 2 pulsations (1 et   )
pour les oscillations libres de la structures . Soit

1er mode de pulsation K 11 − 12 m1

::
K 12 a1
=0
K 21 K 22 − 1 m 2 a 2
2

K 11 −  22 m1 K 12 a1
2ème mode de pulsation =0
K 21 K 22 −  2 m2 a 2
2

23
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 DETERMINATION DU MODE FONDAMENTAL(3)
⚫ Exemple <Portique à 2 niveaux> :
m2 I 21 = I 22 = I 2
I 11 = I 12 = I 1
30*45 0.30 * 0.45 3 27337.5 *10 − 6 4
I2 = = m (Poteau de 30 * 45)
I2 3.00 12 12
I2 m1

0.30 * 0.60 3 64800*10 − 6 4

I1 = = m (Poteau de 30 * 60)
12 12
4.00
30*60
I1 Ei = 3.10 *10 6 (module du béton)
I1

24
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 DETERMINATION DU MODE FONDAMENTAL(4)
Détermination de la matrice de rigidité :

K 11 K 12
K=
K 21 K 22
12EI 21 12EI 22 12EI 11 12EI 12 12EiI 21 12EiI 22 − 24EiI 2
K 11 = 3 + 3 + 3 + 3 K 12 = − − K 12 =
h2 h2 h1 h1 h23 h23 h23
24EiI 2 24EiI 1
K 11 = + 3 − 24 * 3.10 *10 6 * 27337.50 *10 −6
3
h2 h1 K 12 =
I2 I1
12 * 27
K 11 = 24Ei ( +
= −6277.50t / m
)
h 3
2 h13
−6 −6
K 12
27337.5 *10 64800*10
K 11 = 24 * 3.10 *10 6 ( +
K 21 = −6277.50t / m
3
)
12 * 3 12 * 4 3
24 * 3.10 27337.5 64800.00
K 11 = ( + )

K 11
12
= 12555t / m
27 64
K 22 = 6277.50t / ml

12555 − 6277.50
K=
− 6277.50 6277.50
Azzouz Bouguerba-:[Link]-Universiapolis Agadir
25
22/02/2023
DETERMINATION DU MODE FONDAMENTAL(5)
20 0 
La matrice des masses : M M = 
 0 15
PREMIERE METHODE :
Résolution par développement du déterminant :
 12555 − 6277.50 20 0 
K − M = 
− 6277.50 6277.50 
 −    (2 − 4 ) − 1.00 
 0 15
K − M =  
 − 1.00 (1 − 3 ) 
 2.00 − 1.00  4 0
K −   M = 6277.50  − 5    Le déterminant de cette matrice est :
− 1.00 1.00   0 3
( −  ) − 1.00
 2.00 − 1.00  4 0
=
5
− 1.00 (1 − 3 )
K − M =   − 
− 1.00 1.000  6277.50 0 3

 =   − 10 + 
 = ( −  )( −  ) − 1.00
 2.00 − 1.00  4 0
K − M =   −   Les pulsations des vibrations libres sont pour :
− 1.00 1.00   0 3
K − M = 0
  = 0.128
5 =
= 
6277.50  
− 10 +  = 0   = 0.705
5 
=

Azzouz Bouguerba-:[Link]-Universiapolis Agadir
26
22/02/2023
DETERMINATION DU MODE FONDAMENTAL(6)
Détermination des pulsations :

5 52
 = 2 =
 
5 52
0.128 = 0.705 =
 
 = 160.693 2 = 885.1275
 = 12.68rd / s 2 = 29.75rd / s

2 2 * 3.14
Détermination des périodes
T1 = =  T1 = 0.495s
 12.68
2
Ti =
i 2 2 * 3.14
T2 = =  T1 = 0.211s
2 29.75
27
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (1)

m3=25.0T

30*45
I3
I3
3.0
m2=28.0T

30*50 3.0
I2
I2

m3=30.0T

30*60 3.0
I1
I1

28
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (2)
h1 = h2 = h3 = 3.00m ELEMENTS MATRICE DE RIGIDITE
12EI 21 12EI 22 12EI 11 12EI 12 12EI 21 12EI 22
K 11 = 3
+ 3
+ 3
+ K 21 = − −
h h h h13
h23 h23
INERTIES 2 2 1

K21
12E
64800*10 −6 12E
= 3 (2 I 2 + 2 I 1 ) K11 K 12 = − 3 (2 I 2 )
I 11 = I 12 = I 1 = m4 K 11
h
12 h
12 * 3.10 *10 6 * 2 * 37500*10 − 6
12 * 3.10 *10 6 * 2(64800 + 37500) *10 − 6
K 12 =−
37500*10 − 6 K 11 = 27 *12
I 21 = I 22 = I 2 = m4 27 *12
12 2 * 3.10
K 12 =− * 37500 t / m
2 * 3.10 27 K12
27337.50 *10 −6 K 11 = *102300 t / m
I 31 = I 32 = I 3 = m4 27
12
2 * 3.10 * 27337.50
K 23 = − t / m K23
12EI 31 12EI 32 12EI 21 12EI 22 27
K 22 = + + +
h33 h33 h23 h23 Matrice des masses
12E
K 22 = (2 I 3 + 2 I 2 )
h3 1.20 0 0
12 * 3.10 *10 6 * 2(27337.50 + 37500) *10 − 6
K 22 =
27 *12
M = 25 0 1.12 0
K 22 =
2 * 3.10 * 64837.50
t/m K22 0 0 1.00
27
29
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (3)
3.742 − 1.372 0
MATRICE DE RIGIDITE K =   − 1.372 2.372 − 1.00
0 − 1.00 1.00
Résolution du système libre non amorti :
avec :
Le déterminant de la relation doit être nul : soit 6.20 * 27337.50
0 =
12
3.742 − 1.372 0 1.20 0 0
 =   − 1.372 2.372 − 1.00 − 25  0 1.12 0
0 − 1.00 1.00 0 0 1.00

− M + K C = 0 3.742 − 1.372 0 1.20 0 0

 = − 1.372 2.372 − 1.00 −   0 1.12 0
 = K − M 0 − 1.00 1.00 0 0 1.00
et
= Avec :
27
  = 25 
6.20 * 27377.50

1 =
 30
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (4)
ON OBTIENT :
(3.742 − 1.20  ) −  
= −  ( −   ) − 
 −  ( −   )
DEVELOPPEMENT DU DETERMINANT DONNE

(3.742 − 1.20  )(2.372 − 1.12  )(1 −   ) − 1.00 + 1.372− 1.372(1.00 −   ) = 0

  − 6.4275  + 9.84  − 2.494 = 0
Les solutions sont : 
= 0.315
 () = 0.315 
 = 8.894 rd/s
 () = 1.862
 
 () = 4.250 = 1.862

  = 21.622 rd/s
LES PULSATIONS SONT  
= 4.250

  = 32.668 rd/s
31
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (5)
DETERMINATION DES VECTEURS PROPRES PREMIER MODE → MODE 1

⚫ La matrice correspondant à (ω1=8.894 rd/s ) est :

(3.742 − 1.20  ) − 1.372 0 
 
K − M =  − 1.372 ( 2.372 − 1.12  ) −  

  −  ( −   )  VECTEUR
(3.742 − 1.20 * 0.315) − 1.372 0 
PROPRE
K − M =  − 1.372 ( 2.372 − 1.12 * 0.315) −  
  −  ( − 0.315) 

 3.364 − 1.372 0   c1 
K − M = − 1.372

Avec :
 
2.0192
− 
− 
0.685  
C1 =  c2  

 = 

c3 = 1


2
8.894
  = = 0.315
251.10

32
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (6)
( K − M )C1 = 0  3.364 − 1.372 0   c1 
− 1.372 2.0192 −   c  = 0 − 1.00c 2 + 1.00 * 0.685 = 0
  2  c 2 = 0.685
  −  0.685 1.00

Développement de la 1ère ligne donne : 3.364c1 − 1.372c 2 = 0

c1 = 0.279
soit

0.279

C1 = 0.685 
 1.00 

33
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (7)
DEUXIEME MODE : MODE 2   = 1.862

(3.742 − 1.20 12 ) − 1.372 0  1.5076 − 1.372 0   c1 

 
K − M 2 = 

− 1.372 (2.372 − 1.12 1 ) −  
2
− 1.372 0.287 −    c  = 0
  −  ( −  2 )   2 
   −  − 0.862 1.00
(3.742 − 1.20 *1.862) − 1.372 0 

K − M =  − 1.372 (2.372 − 1.12 *1.862) −  
2

  −  ( − 1.862)

1.5076 − 1.372 0 

K − M 2 = − 1.372 0.287 −  
  −  − 0.862
Avec :

 2 =
 2

c 2 = −0.862
 =
2

21.622 2
251.10
= 1.862
 − 0.784 c1 = −0.874

C 2 = − 0.862  
 1.00  34
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (8)
Troisième mode : mode 3  3 = 4.25

(3.742 − 1.20 13 ) − 1.372 0  − 1.358 − 1.372 0   c1 


K − M 3 =  − 1.372 ( 2.372 − 1.12 13 ) −  
 − 1.372 − 2.388 −   c  = 0

  2 
 −  ( −  3 ) 
   −  − 3.25 1.00
(3.742 − 1.20 * 4.25) − 1.372 0 
K − M =  
3 − 1.372 (2.372 − 1.12 * 4.25) −  
  −  ( − 4.25)

− 1.358 − 1.372 0 
K − M = − 1.372

3 − 2.388 − 
  −  − 3.25
Avec :

 =
  3 c 2 = −3.25
 3.284 


2 c1 = 3.284
 
32.688
 3
 = = 4.25

C3 = − 3.25
251.10

 1.00 
35
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE PORTIQUE A 3 NIVEAUX : 3 DDL (9)

Mode : 1 Mode : 2 Mode : 3

1.00 1.00 1.00

-0.862 -3.25
0.685

0.279 -0.784 3.284

36
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
m M = 100 t
E = 3.45 *10 KN/m 7 2
EI

H
I = 0.0185 m 4

EI
H

37
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
DETERMINATION DES ELEMENTS DE LA MATRICE DE RIGIDITE

k11 k12
K=
k21 k22
Les termes de la matrice de rigidité sont
déterminés en appliquant un déplacement
unitaire à chaque degré de liberté et en
calculant les forces élastiques induites.
38
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
k11 V1=1

Déplacements : v1 = 1 et v2 = 0
k11 k 21
k 21
V2=0

H H

k21.H
3
k11.H
2
1
39
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
v1 = 1.00 6 EI
k11 = 2
k11H 3 3k11H 3 5k11H 3 5k 21 H3
1= + + +
3EI 2 EI 6 EI 6 EI 6 EI
6 EI
k0 =
= 16k11 + 5k 21 (1) H3
3
H
v2 = 0.00 k11 = 2k 0
5
k21 = − k11
0=
5k11H 3 k 21H 3
+ 2 k 21 = −5k 0
6 EI 3EI
0 = 5k11 + 2k 21 (2)

40
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
k12 V1=0
Déplacements : v1 = 0 et v2 = 1
6 EI
3
= 16k11 + 5k21 (1)
H k 22 V2=1
L’équation (1) devient en remplaçant :
k11 par k12
Et sachant que : v1 = 0
k 21 par k 22

0 = 16k21 + 5k22 (1' )

Par symétrie on a :
k21 = k12

k22 = 16k0
41
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
Matrice de Masse Matrice de Rigidité
2 −5
m 0 K = k0
M = − 5 16
6 EI
0 m k0 =
7H 3

42
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
PULSATIONS PROPRES DE VIBRATIONS
Pour les oscillations libres non amorties en ne tenant compte que des forces
d’inerties et des forces élastiques : Les pulsations propres :
−5
 + Kv = 0
Mv k0
2
−
m 0
=0
−5 16 0 m
 = −  v
v ( 2k 0 − m  ) − 5k 0

=0
(16k 0 − m )
Kv − M  v = 0 − 5k 0

K − M  = 0
Les solutions sont :

k0
 = 0.63
6 EI 6 * 3.45 *107 * 0.0185 m
k0 = = = 4377
7H 3 7 * 53 k0
4377
2 = 4.195
 = 0.63   = 4.168rd / s m
100 2 2
T =  T1 = = 1.50s
4377  
 = 4.195   = 27.75rd / s
100 2
T2 = = 0.226s
 43
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
Vecteurs propres Mode 1
Le vecteur propre : Φ1
 = aD1
K − M D 2
1 1 =0
1.00 0.32
01 m
( 2k 0 − m ) − 5k 0 d11 D1 M
t
0 m1
=0 a= t =
− 5k 0 (16k 0 − m ) d12 D1 MD1 m 0 1.00
Soit : 1.00 0.32
0 m 0.32
1.00
D1 = 1.32
0.32 a= = 1.197
1.1024
 
1 =   1 =
 

44
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
Vecteurs propres Mode 2
K − M D 2
2 2 =0  2 = aD2
( 2k 0 − m 2 ) − 5k 0 d 21 m 01

=0 − 0.32 1.00
− 5k 0 (16k 0 − m 2 ) 22
d
D M
t
0 m1
Soit : a = t2 =
D2 MD2 m 0 − 0.32
− 0.32 − 1.32 1.00
D2 = 0 m 1.00
1.00
a = 0.617
− 0.32 − 0.
 2 = 0.617  2 =
1.00 0.617

45
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
MODES PROPRES

1.197 -0.197

0.383 0.617

Mode :1 Mode :2
46
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023
 EXEMPLE NUMERIQUE
Pour 7% d’ amortissement ,le spectre
REPONSES MODALES correspondant donne les pseudo accélérations
Accélérations 1er Mode  T1 = 1.50s    = 0.047 g
des masses 2eme Mode  T2 = 0.23s    = 0.11g

 =   
Charges dynamiques 6.02t
Mode 1 :
er
1 Mode : F11 = 100 * 0.056g = 5.60t
1.197 0.056g F21 = 100 * 0.018g = 1.80t
0.047 g =
0.383 0.018g Mode 2 :
2eme Mode : F12 = 100 * ( −0.022g ) = −2.20t 7.03t
− 0.197 − 0.022g F22 = 100 * 0.068g = 6.80t
0.11g =
0.617 0.068g
F1 = 5.602 + ( −2.20) 2 → F1 = 6.02t
F2 = 1.802 + 6.802 → F2 = 7.03t

47
Azzouz Bouguerba-:[Link]-Universiapolis Agadir 22/02/2023