Chapitre VI Etude dynamique

1. Introduction :

La dynamique des structures reprsente une branche de la mcanique

applique ayant pour objet le calcul et ltude du comportement des structures soumises
des vibrations ou aux chocs.

Gnralement on comprend par action dynamique la sollicitation produite par des

charges qui varient rapidement dans le temps (sismes, vent, machines) et qui contribuent
lapparition de forces dinertie.

Le but du prsent chapitre est de dterminer les modes propres de vibration de la

structure, ceci permettra la dtermination des efforts appliqus la structure lorsque cette
dernire est sollicite par un chargement d au vent ou au sisme.

2. Modlisation de la structure :

La modlisation consiste rechercher un modle simplifi qui nous rapproche le plus

possible du comportement rel de la structure en tenant compte le plus correctement possible
de la masse et de la rigidit du systme dorigine, en se basant sur les hypothses suivantes :

- la structure est considre comme une console encastre la base et les masses sont
considres concentres chaque niveau de la structure.
- Chaque masse sera considre avec un seul degr de libert correspondant au
dplacement horizontal (xk).
- Les planchers sont supposs infiniment rigides dans leurs plans.

3. Mthode de calcul :

Gnralement, le calcul exact des caractristiques propres de vibration dune

structure par les mthodes mathmatiques traditionnelles savre extrmement laborieux.
Cest pour cela que des mthodes approches ont t labores (mthodes bases sur
quelques hypothses ou simplifications).

Les mthodes approximatives les plus utilises sont :

- Mthode de LORD RAYLEIGH .

- Mthode de VIANELLO-STODOLA .
- Mthode de HOLZER .
- Mthode de JAKOBI .

Le choix de la mthode sera fait en fonction de la technique de calcul dont on

dispose, de limportance et du type de structure analyse, de la prcision dsire et du
nombre de modes propres qui doivent tre dtermins.

83
Chapitre VI Etude dynamique

Nous utiliserons dans le cadre de notre tude la mthode de VIANELLO-

STODOLA pour ses avantages dapplicabilit, de simplicit et de prcision.

3.1. Principe de la mthode :

On utilise la mthode de STODOLA pour le calcul des formes et des pulsations

propres des systmes ayant nimporte quel nombre de degrs de libert dynamiques. On
suppose que le systme oscillant, ayant n degrs de libert excute un mouvement vibratoire
autour de sa position dquilibre statique.

La mthode de STODOLA est base sur le calcul itratif ; elle consiste faire une
hypothse initiale sur lallure du mode recherch et de lamliorer par un processus
ditrations successives jusqu' lobtention dune approximation satisfaisante du mode.

On considre lquation suivante dun systme n degrs de libert en vibration libre

non amortie :

. { } + []
[] . {} = {} .........(1)

Avec :

[]
, []
: Matrices rigidit et masse.

{ }, {} : Vecteurs acclration et dplacement

Lquation (1) a pour solution la fonction harmonique :

{()} = {} . (. )

Do :

{ ()} = . {} (. )

{ ()} = . {} (. )

84
Chapitre VI Etude dynamique

En substituant dans (1), on obtient ;

. ( . {}
[] (. )) + [] . (. )) = {}
. ({}

En sachant que : -1 (. ) 1

On obtient donc ;

{} []
. ([] ) = {}...............(2)

On dfinit la matrice [F] comme tant la matrice de flexibilit (souplesse) qui est linverse
de la matrice de rigidit [K]

[F] = [K]-1

En multipliant lquation (2) par [F], on obtient :

[] []
([] ) . {}
= {}

([] []
. [] []
) . {}
= {}

([] []
. [] []
) . {}
= {}

([] []
[]
) . {}
= {}..............(3)

On pose : =

En multipliant lquation (3) par , on obtient ;

[]
([] []) . {}
= {}...................(4)

85
Chapitre VI Etude dynamique

[]
On pose : [] = []

En remplaant dans lquation (4), on obtient ;

{}
[] = [] {}
................(5)

On pose : {} {}
= []

En remplaant dans lquation (5), on obtient ;

= . {}
{} ............(6)

Cest partir de lquation (6) que les prochains calculs seront mens.

3.2. tapes de calcul :

- Dtermination de la matrice masse [M].

- Dtermination de la matrice flexibilit [F].
- Dtermination de la matrice dynamique [D]= [M] * [F]
- Calcul des trois premiers modes propres.
- Dtermination du coefficient de participation ().

3.2.1. Dtermination de la matrice masse :

Cest une matrice diagonale qui est du mme ordre que le nombre de masses
supposes concentres de la structure, dans notre cas ceci correspond une matrice [8x8]

Les masses dtages ayant dj t calcul prcdemment, on obtient les rsultats suivants :

456.07 0 0 0 0 0 0 0
0 400 0 0 0 0 0 0
0 0 400 0 0 0 0 0
0 0 0 400 0 0 0 0
0 0 0 0 400 0 0 0
0 0 0 0 0 400 0 0
0 0 0 0 0 0 414.05 0
0 0 0 0 0 0 0 358.56

Les masses sont donnes en tonnes (t)

86
Chapitre VI Etude dynamique

3.2.2. Dtermination de la matrice flexibilit :

Cest une matrice symtrique dont les termes sont nots ij , tel que :

ij : dplacement en i d une force unitaire applique en j.

En considrant le thorme de MAXWELL suivant :

Le dplacement produit dans la section 1 par une force P agissant dans la section 2 est gal
au dplacement produit dans la section 2 si la mme force P agit dans la section 1

On conclue que la matrice flexibilit [F] est une matrice symtrique avec : ij = ji

En considrant leffet de tous les efforts, lexpression du dplacement dans le cas gnral est
donne comme suit :

. . . .
= = [ . + . + . + . ]
. . . .

Avec ;
, , , : Effort normal, moment de flexion, effort tranchant et moment de torsion
donns par des charges relles.

, , , : Effort normal, moment de flexion, effort tranchant et moment de torsion

donns par des charges de valeur unitaire agissant dans la section o le dplacement est
recherch.

Les dplacements dus au moment flchissant sont prpondrant par rapport aux
autres efforts, donc lexpression du dplacement devient comme suit :
.
= = [ . ]
.

Il existe plusieurs mthodes pour le calcul des dplacements, nous utiliserons la

mthode de MAXWELL-MOHR en considrant que la structure est une console encastre
sa base.

Selon la mthode, on utilise les deux systmes de force suivants :

- Le premier systme est constitu par les forces relles agissant sur la structure ce qui
correspond dans notre cas une force unitaire applique en j .

87
Chapitre VI Etude dynamique

- Le deuxime systme est form par une seule force unitaire qui agit selon la
direction du dplacement recherch ce qui correspond dans notre cas une force
unitaire applique en i .

Il faut par la suite tracer les diagrammes des efforts correspondants pour chaque
systme de force puis effectuer les intgrales donnes par la formule dans le cas gnral

On trouve alors les rsultats suivants :

1
Sens transversal : = 2,5.10-9 KN.m2
.

8.058 16.116 24.175 32.233 40.291 48.349 56.407 64.466

16.117 37.306 59.690 82.074 104.458 126.842 149.226 171.610
24.174 59.690 102.369 146.241 190.114 233.986 277.858 321.731
[S] (m/KN)= 32.234 82.074 146.241 217.571 290.095 362.619 435.142 507.666
40.291 104.458 190.114 290.095 397.239 505.577 613.915 722.253 .10-8
48.350 126.842 233.987 362.619 505.578 655.698 807.013 958.328
56.408 149.226 277.858 435.143 613.915 807.012 1007.274 1208.729
64.465 171.610 321.731 507.665 722.253 958.328 1208.729 1466.293

1
Sens longitudinal : = 1.8.10-9 KN.m2
.

5.80 11.59 17.39 23.19 28.98 34.78 40.58 46.37

11.59 26.84 42.94 59.04 75.14 91.25 107.35 123.45
17.39 42.94 73.64 105.20 136.76 168.32 199.88 231.44
[S] (m/KN) = 23.19 59.04 105.20 156.51 208.69 260.86 313.03 365.20 .10-8
28.98 75.14 136.76 208.69 285.76 363.70 441.63 519.57
34.78 91.25 168.32 260.86 363.70 471.69 580.54 689.39
40.58 107.35 199.88 313.03 441.63 580.54 724.60 869.52
46.37 123.45 231.44 365.20 519.57 689.39 869.52 1054.81
88
Chapitre VI Etude dynamique

3.2.3. Dtermination de la matrice dynamique :

La matrice dynamique [D] est obtenue en multipliant la matrice flexibilit [F] par la
matrice masse [M], et le rsultat est donn comme suit :

Sens transversal :

3.675 6.447 9.670 12.893 16.116 19.340 23.355 23.115

7.350 14.923 23.876 32.830 41.783 50.737 61.787 61.532
11.025 23.876 40.948 58.497 76.045 93.594 115.047 115.360
[D] = 14.701 32.830 58.496 87.028 116.038 145.047 180.170 182.029 . 10-5
18.375 41.783 76.046 116.038 158.896 202.231 254.190 258.971
22.051 50.737 93.595 145.048 202.231 262.279 334.142 343.618
25.726 59.690 111.143 174.057 245.566 322.805 417.060 433.401
29.400 68.644 128.692 203.066 288.901 383.331 500.472 525.753

Sens longitudinal :

2.644 4.637 6.956 9.275 11.594 13.912 16.801 16.628

5.288 10.735 17.176 23.617 30.058 36.498 44.447 44.264
7.931 17.176 29.456 42.081 54.705 67.329 82.761 82.986
[D] = 10.575 23.617 42.081 62.606 83.474 104.343 129.609 130.946 .10-5
13.219 30.058 54.705 83.474 114.305 145.479 182.856 186.295
15.863 36.498 67.329 104.343 145.479 188.676 240.371 247.188
18.506 42.939 79.953 125.211 176.653 232.216 300.020 311.775
21.150 49.380 92.577 146.079 207.827 275.756 360.024 378.211

3.2.4. Calcul des modes propres :

La mthode consiste faire des approximations successives jusqu convergence

vers le mode dominant (1er mode, 2me mode,.....etc.)

1er mode :
- Premire itration :
On note le vecteur propre comme suit : {i}j

Avec i : reprsente le mode

J : reprsente litration

89
Chapitre VI Etude dynamique

On commence par choisir arbitrairement le vecteur propre de la premire approximation

{1}1

On calcule ensuite le vecteur {y} comme suit :

{ } = []. { }

Puis, on calcule comme suit :

{ }
= { }
en choisissant les derniers termes des deux vecteurs

Et nous pouvons ainsi dterminer le vecteur propre de la deuxime approximation

{ }
{ } =

- Deuxime itration :

On prend pour la premire approximation le vecteur propre {1 }2 et on rpte le mme

processus itratif jusqu convergence, cest- dire jusqu atteindre :

{ } = { }

La pulsation propre et la priode sont dtermines comme suit :

1 2
12 = ; =
1

2me mode :

Pour que le processus itratif converge vers le 2me mode, on doit supprimer le 1er mode,
pour cela on utilise la matrice de balayage suivante :

0 a12 . a1i a1n

0 1 0
[S1] = 0 1 0

0 1 0
1
0

Avec,
.
=
.

m1, mi : masses des niveaux 1 et i respectivement

90
Chapitre VI Etude dynamique

Et on refait les calculs avec : [D1] = [D]* [S1] jusqu convergence du mode 2 et ainsi nous
aurons le vecteur propre , la pulsation 2 et la priode T2.

3me mode :

Pour que le processus itratif converge vers le 3me mode, on doit supprimer le 1er et le 2me
mode, pour cela on utilise la matrice de balayage suivante :

1 0 . 0 0

0 0 b23 b2i b2 n
[S2] = 0 1 0 0

0 1 0
1
0

Avec :

(2 . 11 1 . 12 )
2 =
(22 . 11 21 . 12 )2

Et on refait les calculs avec : [D2] = [D1]* [S2] jusqu convergence du mode 3 et ainsi nous
aurons le vecteur propre , la pulsation 3 et la priode T3.

3.2.5. Dtermination du coefficient de participation modale :

Le coefficient de participation modale correspond un seul mode de vibration, il dfinit le

pourcentage dnergie absorbe durant ce mode. Il est donn par la formule suivante :

(1 . )2 1
= 2
1 . 1

Daprs larticle 4.3.4 du RPA, le nombre de modes considrer est tel que :

> %

91
Chapitre VI Etude dynamique

Les rsultats des vecteurs propres, des pulsations propres, des priodes et des coefficients de
participation sont rsums dans les schmas suivants pour chaque mode et dans chaque
sens :

= 95.55%
1

92
Chapitre VI Etude dynamique

= 95.52%
1

93
Chapitre VI Etude dynamique

4. Vrification de lorthogonalit :

On doit vrifier que : { } . []. { } = 0

- Sens longitudinal :

{1 } . []. {2 } = 0.. 1 1 3.12*10-5 0

{2 } . []. {3 } = 0..2 Aprs calcul : 2 5.33*10-5 0
{1 } . []. {3 } = 0.. 3 3 2.45*10-5 0

Les conditions sont vrifies.

- Sens transversal :

{1 } . []. {2 } = 0.. 1 1 2.55*10-5 0

{2 } . []. {3 } = 0 ..2 Aprs calcul : 2 4.39*10-5 0
{1 } . []. {3 } = 0.. 3 3 2.02*10-5 0

Les conditions sont vrifies.

94