CorrigΓ© 10 - IMAC

EPFL
IMAC-IS-ENAC
Cours de Dynamique des structures
Prof. I. Smith
EXERCICE 10
- Corrigé -
CORRIGE EXERCICE No 10
A.
m , Ix
x
L = 14 m
Choix de la forme de la déformée approximant le premier mode et qui satisfait les conditions
de bord :
π‘₯
πœ“(π‘₯) = βˆ’sin (πœ‹ )
𝐿
πœ‹
π‘₯
𝐿
𝐿
πœ“β€²(π‘₯) = βˆ’ cos (πœ‹ )
Avec :
πœ‹ 2
π‘₯
πœ“β€²β€²(π‘₯) = ( ) sin (πœ‹ )
𝐿
𝐿
𝐿
Données :
a)
π‘₯
𝐿
∫0 sin2 (πœ‹ 𝐿 ) 𝑑π‘₯ = 2
Rappel :
m=6 t/m ;
h=1 m ;
L=14 m ;
Ix=40.5*104 cm4 ;
E=21*106 kN/m2
M=20 t ;
Fréquence propre de l’ouvrage
𝐿
π‘˜ βˆ— = ∫ 𝐸𝐼(π‘₯)[πœ“β€²β€²(π‘₯)]2 𝑑π‘₯ =
0
πΈπΌπœ‹ 4 𝐿 2 π‘₯
πΈπΌπœ‹ 4
π‘˜π‘
∫
sin
(πœ‹
)
𝑑π‘₯
=
β‰… 1.5 βˆ™ 103
4
3
𝐿
𝐿
2𝐿
π‘š
0
𝐿
𝐿
π‘₯
π‘š
̅𝐿
π‘šβˆ— = ∫ π‘š(π‘₯)[πœ“(π‘₯)]2 𝑑π‘₯ = π‘š
Μ… ∫ sin2 (πœ‹ ) 𝑑π‘₯ =
β‰… 42𝑑
𝐿
2
0
0
πœ”π‘› = √
b)
π‘˜βˆ—
β‰… 6π‘Ÿπ‘Žπ‘‘/𝑠
π‘šβˆ—
Fréquence propre de l’ouvrage modifié
π‘˜ βˆ— = 1.5 βˆ™ 103
π‘˜π‘
π‘š
𝐿
π‘šβˆ— = ∫ π‘š(π‘₯)[πœ“(π‘₯)]2 𝑑π‘₯ + 𝑀[πœ“(π‘₯ = 𝐿/2)]2 =
0
πœ”π‘› = √
c)
π‘š
̅𝐿
+ 𝑀 β‰… 62𝑑
2
π‘˜βˆ—
β‰… 4.9π‘Ÿπ‘Žπ‘‘/𝑠
π‘šβˆ—
Force statique équivalente
Choc mou :
Force due au camion :
πΉπ‘β„Žπ‘œπ‘ = πœ”π‘› π‘š1 𝑣1 = 30.32 βˆ™ 3 βˆ™ √2 βˆ™ 9.81 βˆ™ 1 β‰… 65π‘˜π‘
𝐹𝑐 = (𝑀 βˆ’ π‘š1 )𝑔 = (20 βˆ’ 3) βˆ™ 9.81 β‰… 167π‘˜π‘
Page 1/3
EPFL
IMAC-IS-ENAC
Cours de Dynamique des structures
Prof. I. Smith
πΉπ‘šπ‘Žπ‘₯ = πΉπ‘β„Žπ‘œπ‘ + 𝐹𝑐 = 65 + 167 = 232π‘˜π‘
Force statique équivalente :
d)
EXERCICE 10
- Corrigé -
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 :
π›Ώπ‘€π‘Ÿπ‘’π‘ π‘ π‘œπ‘Ÿπ‘‘ = 𝐾[πœ“β€²(π‘₯ = 0)]2 =
πœ‹2
𝐾
𝐿2
Où πœ“β€²(π‘₯ = 0) est la rotation interne du ressort. Ainsi :
𝐿
π‘˜ βˆ— = ∫ 𝐸𝐼(π‘₯)[πœ“β€²β€²(π‘₯)]2 𝑑π‘₯ + π›Ώπ‘€π‘Ÿπ‘’π‘ π‘ π‘œπ‘Ÿπ‘‘ =
0
πΈπΌπœ‹ 4 πœ‹ 2
π‘˜π‘
+ 2 𝐾 β‰… 2.0 βˆ™ 103
2𝐿3
𝐿
π‘š
π‘šβˆ— = 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 𝐹𝑖 𝑑𝑖
2
βˆ‘π‘+1
𝑖=1 π‘šπ‘– 𝑑𝑖
Plus le nombre d’intervalles augmentent, plus la solution se rapproche de la valeure 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 𝐾 :
0
102
103
𝐾 [kNm/rad]
1ère approximation [rad/s] 5.99
6.00
6.09
Méthode 1 ou 2 [rad/s]
5.99
6.00
6.10
104
6.92
6.77
105
12.48
8.53
106
35.14
9.28
109
1095
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.
Page 2/3
EPFL
IMAC-IS-ENAC
Cours de Dynamique des structures
Prof. I. Smith
EXERCICE 10
- Corrigé -
B. 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 :
Type d’action
Valeur de l’effort associé
Désignation
Incrément du
déplacement en
x (t )
fonction de
Force due au
réssort
Force due à
l’amortisseur
FS 1
Kx
x
Kx x
FD 1
Cx
x
Cx x
FI 1
M 0 3L
MI 1
Effets de
l’inertie
I0
M 0 3L
x
2
3
M 0Lx
2
(3L)2 x
12 3L
I0
MI 2
M
M0
2
L x
8 L
2
9L2
12
3
x
M
4
M 0L2x
x
2
3
x
M 0Lx
2
2
x
3L
3
x
M L2x
4 0
3L
x
L
x
4
L
M
2
x
3
2L P(t )
F
Force externe
x
Contribution au
travail virtuel
P (t )2L
Travail effectué lors d'un déplacement virtuel compatible avec les liaisons :
wi
we
3
x
M 0Lx
2
2
3
x
M 0L2x
4
3L
M
L 2 x
x
4 L
Kx x
Cx x 0
P (t )2L
x
3
Equation du mouvement :
M 0Lx
1
Mx
2
Kx
Cx
2
P(t ) L
3
0
(M 0L
1
M )x
2
Cx
Kx
L 2 x
x
4 L
2
P (t ) L
3
Rappel : Les travaux des forces d’inertie sont toujours positifs. Les travaux des
ressorts et des amortisseurs sont toujours négatifs.
Page 3/3
x
3