EPFL Cours de Dynamique des structures Série 09

IMAC-IS-ENAC Prof. I. Smith / Dr. P. Lestuzzi -Corrigé-

Corrigé de la série d'exercices N° 9

Exercice 1

m, E , I x

A B
x
L=14m
Choix de la forme de la déformée approximant le premier mode et qui satisfait les conditions de bord :

𝑥
𝜓(𝑥) = −sin (𝜋 )
𝐿

𝜋 𝑥
Avec : 𝜓′(𝑥) = − cos (𝜋 )
𝐿 𝐿

𝜋 2 𝑥
𝜓′′(𝑥) = ( ) sin (𝜋 )
𝐿 𝐿
𝐿 𝑥 𝐿
Rappel : ∫0 sin2 (𝜋 𝐿 ) 𝑑𝑥 = 2

Données : m =6 t/m ; L=14 m ; Ix=40.5*104 cm4 ; M=20 t ;

6 2
h=1 m ; E=21*10 kN/m

1. Fréquence propre de l’ouvrage

𝐿
𝐸𝐼𝜋 4 𝐿 2 𝑥 𝐸𝐼𝜋 4 𝑘𝑁
𝑘 ∗ = ∫ 𝐸𝐼(𝑥)[𝜓′′(𝑥)]2 𝑑𝑥 = 4
∫ sin (𝜋 ) 𝑑𝑥 = 3
≅ 1.5 ∙ 103
0 𝐿 0 𝐿 2𝐿 𝑚

𝐿 𝐿
𝑥 𝑚
̅𝐿
𝑚∗ = ∫ 𝑚(𝑥)[𝜓(𝑥)]2 𝑑𝑥 = 𝑚
̅ ∫ sin2 (𝜋 ) 𝑑𝑥 = ≅ 42𝑡
0 0 𝐿 2

𝑘∗
𝜔𝑛 = √ ≅ 6𝑟𝑎𝑑/𝑠
𝑚∗

2. Fréquence propre de l’ouvrage modifié

𝑘𝑁
𝑘 ∗ = 1.5 ∙ 103
𝑚
𝐿
𝑚
̅𝐿
𝑚∗ = ∫ 𝑚(𝑥)[𝜓(𝑥)]2 𝑑𝑥 + 𝑀[𝜓(𝑥 = 𝐿/2)]2 = + 𝑀 ≅ 62𝑡
0 2

𝑘∗
𝜔𝑛 = √ ≅ 4.9𝑟𝑎𝑑/𝑠
𝑚∗

3. Force statique équivalente

Choc mou : 𝐹𝑐ℎ𝑜𝑐 = 𝜔𝑛 𝑚1 𝑣1 = 4.9 ∙ 3 ∙ √2 ∙ 9.81 ∙ 1 ≅ 65𝑘𝑁

Force due au camion : 𝐹𝑐 = (𝑀 − 𝑚1 )𝑔 = (20 − 3) ∙ 9.81 ≅ 167𝑘𝑁

Force statique équivalente : 𝐹𝑚𝑎𝑥 = 𝐹𝑐ℎ𝑜𝑐 + 𝐹𝑐 = 65 + 167 = 232𝑘𝑁

1/4
This study source was downloaded by 100000842064575 from [Link] on 02-19-2022 08:11:02 GMT -06:00

[Link]
 EPFL Cours de Dynamique des structures Série 09
IMAC-IS-ENAC Prof. I. Smith / Dr. P. Lestuzzi -Corrigé-

4. Ressort sur appui gauche

𝑥
On suppose que la forme de la déformée reste la même 𝜓(𝑥) = −sin (𝜋 ). Pour le calcul de 𝑘 ∗ , il faut en plus
𝐿
considérer le travail virtuel du ressort :
𝜋2
𝛿𝑤𝑟𝑒𝑠𝑠𝑜𝑟𝑡 = 𝐾[𝜓′(𝑥 = 0)]2 = 2 𝐾
𝐿

Où 𝜓′(𝑥 = 0) est la rotation interne du ressort. Ainsi :

𝐿
𝐸𝐼𝜋 4 𝜋 2 𝑘𝑁
𝑘 ∗ = ∫ 𝐸𝐼(𝑥)[𝜓′′(𝑥)]2 𝑑𝑥 + 𝛿𝑤𝑟𝑒𝑠𝑠𝑜𝑟𝑡 = 3
+ 2 𝐾 ≅ 2.0 ∙ 103
0 2𝐿 𝐿 𝑚

𝑚∗ = 42𝑡

𝑘∗
𝜔𝑛 = √ ≅ 6.9𝑟𝑎𝑑/𝑠
𝑚∗

En réalité, l’approximation de la déformée par la fonction sinusoïdale n’est pas valable pour des rigidités plus
grandes. Deux méthodes sont possibles pour obtenir une meilleure approximation.

Méthode 1 :
1) Déterminer la déformée sous une charge répartie 𝑞 (par exemple : poids propre) :
𝑞(𝐿 − 𝑥)𝑥[𝐾𝐿(3𝐿 − 2𝑥)𝑥 + 6𝐸𝐼(𝐿2 + 𝐿𝑥 − 𝑥 2 )]
𝜓(𝑥) = −
48𝐸𝐼(3𝐸𝐼 + 𝐾𝐿)
2) Déterminer 𝛿𝑤𝑓𝑙𝑒𝑥𝑖𝑜𝑛 , 𝛿𝑤𝑟𝑒𝑠𝑠𝑜𝑟𝑡 et 𝛿𝑤𝑖𝑛𝑒𝑟𝑡𝑖𝑒 et ainsi 𝑘 ∗ et 𝑚∗ .
3) Finalement,
𝑘∗ 𝐸𝐼(24𝐸 2 𝐼 2 + 11𝐸𝐼𝐾𝐿 + 𝐾 2 𝐿2 )
𝜔𝑛 = √ = 18√14√
𝑚∗ 𝑚̅ 𝐿4 (1116𝐸 2 𝐼 2 + 285𝐸𝐼𝐾𝐿 + 19𝐾 2 𝐿2 )
Méthode 2 :
Si l’équation de la déformée est trop compliquée à obtenir et que l’on dispose d’un programme d’éléments finis :
1) Discrétiser le système en 𝑁 intervalles. Appliquer 𝑁 + 1 forces 𝐹𝑖 uniformément répartie sur le
système et déterminer la masse par intervalle 𝑚𝑖 .
2) Calculer les déplacements 𝑑𝑖 au droit des forces 𝐹𝑖 avec la méthode des éléments finis.
3) Finalement,
∑𝑁+1
𝑖=1 𝐹𝑖 𝑑𝑖
𝜔𝑛 = √ 𝑁+1
∑𝑖=1 𝑚𝑖 𝑑𝑖2
Plus le nombre d’intervalles augmentent, plus la solution se rapproche la valeur exacte de la pulsation propre
(par le haut). Si 𝑁 → ∞, la solution est égale à celle de la méthode 1.

Voici la solution pour différentes valeurs de 𝐾 :

𝐾 [kNm/rad] 0 102 103 104 105 106 109
ère
1 approximation [rad/s] 5.99 6.00 6.09 6.92 12.48 35.14 1095
Méthode 1 ou 2 [rad/s] 5.99 6.00 6.10 6.77 8.53 9.28 9.39

Avec la 1ère approximation (déformée sinusoïdale), la pulsation diverge fortement de l’approximation selon la
méthode 1 ou 2 pour 𝐾 ≥105. Il faut donc être prudent sur le choix de la forme de la déformée. La meilleure
approximation en cas de doute est la déformée sous charge répartie.

Exercice 2
La réponse u(x, t) de la structure peut être écrite comme :
u  x, t     x  z  t 

2/4
This study source was downloaded by 100000842064575 from [Link] on 02-19-2022 08:11:02 GMT -06:00

[Link]
 EPFL Cours de Dynamique des structures Série 09
IMAC-IS-ENAC Prof. I. Smith / Dr. P. Lestuzzi -Corrigé-

Avec   x  une fonction de forme approximant la déformée du premier mode et qui satisfait les conditions de
bord.
On peut choisir une fonction parabolique du type :   x   a x
c 2

Avec les deux conditions de bord :   x  0   0 et   x  H   1

x2
On obtient :   x   2
H

u(x, t)

F*(t) z(t)
m*

q(t) H

k*

La structure étudiée peut alors être assimilé à un oscillateur simple équivalent dont les caractéristiques sont les
suivantes :
𝐻
4𝐸𝐼
𝑘 ∗ = ∫ 𝐸𝐼(𝑥)[𝜓′′(𝑥)]2 𝑑𝑥 = 3 = 1097.4 𝑘𝑁/𝑚
0 𝐻

𝐻
𝑚
̅𝐻
𝑚∗ = ∫ 𝑚(𝑥)[𝜓(𝑥)]2 𝑑𝑥 = = 2520 𝑘𝑔
0 5

𝑘∗
𝜔𝑛 = √ = 20.86 𝑟𝑎𝑑/𝑠
𝑚∗
𝐻
𝑞0 𝐻
𝐹 ∗ (𝑡) = ∫ 𝑞(𝑥, 𝑡)𝜓(𝑥) dx = sin 𝜔𝑡
0 3

Le système équivalent est un oscillateur simple non amorti soumis à une force d’excitation sinusoïdale. Sa
réponse s’écrit donc de la manière suivante :

q0 H
z (t )  3k * sin t   0.566sin 15t 
1   2 n2

La correction par rapport à la solution particulière est relativement limitée si les conditions initiales sont nulles
même sans amortissement.

Et la réponse de la structure s’écrit donc :

q0 H q0
2
x
u  x, t     x  z  t  = 2 . 
3k *
sin  t   3Hk * x 2 sin t 
H 1   2 n2 1   2 n2

Application numérique :

3/4
This study source was downloaded by 100000842064575 from [Link] on 02-19-2022 08:11:02 GMT -06:00

[Link]
 EPFL Cours de Dynamique des structures Série 09
IMAC-IS-ENAC Prof. I. Smith / Dr. P. Lestuzzi -Corrigé-

u  x, t   1.753 x 2 sin 15t 

Exercice 3
Détermination de l’équation du mouvement

+
MI2
Convention: +
+

MI1

FI1

Sous l'effet d'un déplacement x(t), les éléments apportent les contributions suivantes :

Incrément du Contribution
Type Désignati Valeur de l’effort associé x (t ) déplacement en au travail
d’action on
fonction de x virtuel
Force due
FS 1 Kx x Kx x
au réssort
Force due à
l’amortisseu FD1 Cx x Cx x
r
x 3 x 3 x
FI 1 M 0 3L M Lx M 0Lx
2 2 0 2 2 2
(3L)2 x 9L2 3 x 3 x
Effets de MI1 I0 M 0 3L M0 x M 0L2x M 0L2x
l’inertie 12 3L 12 4 3L 4 3L
2
L x L x L 2 x
MI 2 I0 M M x M x
8 L 4 L 4 L
2 2
Force x x
F 2L P(t) P (t )2L
externe 3 3

Travail effectué lors d'un déplacement virtuel compatible avec les liaisons :
wi we
3 x 3 x L 2 x x
M 0Lx M 0L2x M x Kx x Cx x 0 P (t )2L
2 2 4 3L 4 L 3
Equation du mouvement :

4/4
This study source was downloaded by 100000842064575 from [Link] on 02-19-2022 08:11:02 GMT -06:00

[Link]
 EPFL Cours de Dynamique des structures Série 09
IMAC-IS-ENAC Prof. I. Smith / Dr. P. Lestuzzi -Corrigé-

1 2 1 2
M 0Lx Mx Kx Cx P(t ) L 0 (M 0L M )x Cx Kx P(t ) L
2 3 2 3

Rappel : Les travaux des forces d’inertie sont toujours positifs. Les travaux des ressorts et des amortisseurs sont
toujours négatifs.

5/4
This study source was downloaded by 100000842064575 from [Link] on 02-19-2022 08:11:02 GMT -06:00

[Link]
Powered by TCPDF ([Link])