3A103 Vibrations 1617
3A103 Vibrations 1617
3A103 Vibrations 1617
Vibra&ons
26 h de cours
26 h de TD
8 h de TP
Evalua&on :
Examen (80 %) :
TP (20%)
Enseignants
Cours : C. Maurini (E&S) S. Rgnier (Vibra&ons)
corrado.maurini@upmc.fr & stephane.regnier@upmc.fr
Le cours sur powerpoint a t ini&alement produit par F. Ollivier qui est remerci ici
F. Ollivier S. Rgnier 16/17
page 1
Vibra&ons
page 2
Vibra&ons
1. ExcitaJon harmonique
(4 h)
2. ExcitaJon priodique
3. ExcitaJon quelconque Mthode de Laplace
page 3
Vibra&ons
FroZements
Liaisons
Sta&que
Cinma&que
Cin&que
Solides dformables
(ou las6ques)
Sta&que et dynamique
des Systmes mul&corps
Dynamique
page 4
Vibra&ons
Avant (Calcul)
Dimensionner
Modliser
Rles de lingnieur
Construire
Aprs (Mesure)
Surveiller
DiagnosJquer
Modier
page 5
Vibra&ons
Exemples de Modles
Corps humain
Modle discret
Rail-Ballast
Modle con2nu
Y
V0t
V0
Globes
oculaires
Bras
Bras
B
Thorax
Automobile
modle 2ddl
z
Colonne
vertbrale
F0(x,t)
G
2L
Turbine de gnrateur
modle 1ddl
Bassin
c
Jambes
Pied
Pied
r
y
B2ment
n tages
modle n ddl
r
Fx
r
x
page 6
Vibra&ons
page 7
Vibra&ons
page 8
Vibra&ons
page 9
Vibra&ons
Modles con2nus
nombre de ddl
(cordes - poutres - plaques)
Ecriture : analyJque
Solu6on : analyJque exacte
page 10
Vibra&ons
Problme
Quelles causes ?
Quel mouvement ?
Quels paramtres ?
ObservaJons
Mouvement :
OscillaJon autour
d une posiJon
d quilibre staJque
Paramtres inuents :
Iner&e
masse, moment
Elas&cit
Amor&ssement
page 11
Vibra&ons
Systme 1DDL
page 12
Vibra&ons
c
k, c
En transla&on
Paramtrage)
Dplacement+
Inertie+
Raideur+
Amortissement+
Pendulaire
En rota&on
Translation*
Rotation*
Pendulaire*
Longitudinal+:+x+
Angulaire+:++
Angulaire+:++
Moment+dinertie+:+
Masse+:+M+
Masse+:+M+
I+
Rsistance++
Rsistance++
Pesanteur+
+lallongement+:+k+
+la+torsion+:+K+
Frottements+visqueux+:+c+ou+C+
ObjecJf :
Dterminer la forme du mouvement de la structure en foncJon :
- du temps : x(t), (t)
- de la frquence de l excitaJon : x(), ()
page 13
Vibra&ons
Masse sismique
Base
!
x(t)
!
! (t)
!
x
Bielle
M
Pivot
c
Axe vilebrequin
Bloc moteur
page 14
Vibra&ons
Arbre
et hlice
I
K, C
!
y
BJment
!
Fx
m
k
!
x
r
x (t )
(t )
Passerelle
!
x
c
F. Ollivier S. Rgnier 16/17
k
page 15
Vibra&ons
La cinmaJque
La cinJque
page 16
Vibra&ons
page 17
Vibra&ons
page 18
Vibra&ons
Une Jge (T) numrote corps 1 de longueur L a une liaison pivot daxe z0
avec le bJ
Un disque (D) numrote 2 de rayon R a une liaison pivot daxe z0 avec la
Jge (T)
Calculer le torseur cinma&que du disque par rapport au b& au point P
F. Ollivier S. Rgnier 16/17
page 19
Vibra&ons
page 20
Vibra&ons
page 21
Vibra&ons
page 22
Vibra&ons
page 23
Vibra&ons
page 24
Vibra&ons
page 25
Vibra&ons
page 26
Vibra&ons
page 27
Vibra&ons
La matrice dinerJe principale scrit : pour une Jge de longueur L et daxe z
page 28
Vibra&ons
Torseur cinJque
Le torseur cin&que au point A scrit donc :
page 29
Vibra&ons
Energie cinJque
On appelle nergie cin&que dun systme matriel con&nu la quan&t
scalaire exprime en J
page 30
Vibra&ons
Torseur dynamique
LacclraJon dun point M est note :
page 31
Vibra&ons
page 32
Vibra&ons
Le systme est compos dun disque plein (D) plein et homogne et dune Jge
homogne (T).
La liaison entre le disque et la Jge en C est une liaison pivot parfaite sans
frotement.
Le repre R0 est le repre xe li lobservateur et considr comme galilen.
Le mouvement du systme seectue dans le plan (O, x0, y0)
F. Ollivier S. Rgnier 16/17
page 33
Vibra&ons
page 34
Vibra&ons
(Em ) = Pe
(Ep +Ec ) = Pe
page 35
Vibra&ons
Equations de Lagrange
Les quaJons de Lagrange sont un ouJl trs praJque pour obtenir les
quaJons de mouvement notamment dans le cas de systmes
Il sagit de lexpression du PFD sous une forme nergJque.
On peut dmontrer que :
d L L Ed
+
dt x x x
!
Le Lagrangien L scrit comme :
Pe
=
x
L = Ec Ep
!
F. Ollivier S. Rgnier 16/17
page 36
Vibra&ons
page 37
Vibra&ons
(en transla6on)
Conserva2on de l nergie (Lagrange) :
2 mthodes
Bilan des forces (PFD) :
Force d excitaJon :
F(t)
Force de rappel lasJque :
-k x(t)
Force de frotement visqueux : -c x! ( t )
Forces = m
Energie dissipe :
Ed =
Puissance extrieure :
F kx cx! = m!x!
c
!
F
m
k
nergie cinJque :
nergie potenJelle
lasJque :
E c = 12 mx2
!
Ep = 12 kx2
!
x
1 2
c x
2
Pe = F x
d L L Ed Pe
+
=
dt x x x x
!
mx +kx + cx = F
m
x
+
c
x
+kx
=F
!
page 38
Vibra&ons
I
K, C
m
k
!
x(t)
!
! (t)
M
F. Ollivier S. Rgnier 16/17
page 39
Vibra&ons
Rponse
VibraJon naturelle de la structure
sans amor2ssement :
avec amor2ssement :
Rponse
VibraJon force de la structure
Mouvement accord lexcita2on :
Amplitude et phase
selon le spectre de lexcita2on
+
page 40
Vibra&ons
quaJon du mouvement :
Solu6on :
3 critures
quivalentes
j t
j t
x(t) = X1e 0 + X 2e 0
x(t) = A sin(0t) + Bcos(0t)
x(t) = X cos( 0 t )
Condi6ons ini6ales :
x(0) = x0
x (0) = v0
k
x(t) = 0
m
k
0 =
m
0 : pulsa6on propre,
Rq : foncJon du rapport
de l lasJcit l inerJe.
Constantes
A, B ou X1, X2 ou X,
v0
x(t) =
sin( 0 t) + x 0 cos( 0 t) = X cos( 0 t )
0
X : amplitude,
: phase
F. Ollivier S. Rgnier 16/17
page 41
Vibra&ons
page 42
Vibra&ons
x(t) = X cos( 0t )
C est la mme criture
2
X = A +B
tg =
A
B
page 43
Vibra&ons
x(t)
x0
v0
X cos(0t)
t =
tg =
v0
0 x 0
= 0 t
page 44
Vibra&ons
x(t) = X cos(0 t )
w0 X
= 0 X cos( 0 t + )
2
w02 X
t
Acclration
x(t) = 0 2 X cos( 0 t )
= 0 2 X cos( 0 t + )
Dplacement et vitesse
en quadrature
Dplacement et acclraJon
en opposiJon de phase
F. Ollivier S. Rgnier 16/17
page 45
Vibra&ons
m
x(t) + cx (t) + kx(t) = 0
x(t) +
x(t) = Xest
Solutions de la forme :
quation caractristique :
Racines :
c
x (t) + 02 x(t) = 0
m
c
2
s + s + 0 = 0
m
2
c
c
2
2
s=
1
0
0
0
2m
2m
c
=
2m 0
: facteur damortissement
3 cas damortissement selon
F. Ollivier S. Rgnier 16/17
page 46
Vibra&ons
page 47
Vibra&ons
v0
x(t) = 0 2 1
x(t) = Aet + Be t e
0 t
t
= X1 ch(t) + X 2 sh(t) e 0
0t
0 x 0 + v 0
= x 0 ch(t) +
sh(t) e
v0>0
v0<0
x0 exp(-0t)
t
F. Ollivier S. Rgnier 16/17
page 48
Vibra&ons
A = x0
B = x 0 0 + v 0
F. Ollivier S. Rgnier 16/17
page 49
Vibra&ons
x(t) = [ x0 (1 + 0t) + v0 t ] e0 t
v0<0
x0 exp(-w0t)
t
F. Ollivier S. Rgnier 16/17
page 50
Vibra&ons
page 51
Vibra&ons
x(t) = X e0 t cos(dt )
Pseudo pulsation
x(0)
d = 0 1 2
Td = 2/d
X exp(-0t)
X = x0
v + 0 x 0
+ 0
v 0 + x 00
tg =
dx 0
F. Ollivier S. Rgnier 16/17
page 52
Vibra&ons
x(t)
Td
Xn
Xn+1
t
x(t 0 )
Xn
Xe 00
T
=
=
=e 0 d
Xn+1 x(t 0 + Td ) X e 0 (t0 +Td )
Dcrment logarithmique : = ln
0
= 0Td = 2
d
=
2
1
xn
xn+1
Rq : On mesure plus
1 Xn
= ln
prcisment avec des
p Xn+p
maxima spars de pTd :
4 +
page 53
Vibra&ons
Spara&on sta&que/dynamique
La soluJon peut galement tre dcompose en une parJe staJque et une
parJe dynamique :
page 54
Vibra&ons
m
k
!
! (t)
!
x(t)
!
x1(t)
!
x2(t)
La posiJon x(t) est repre par rapport la x2(t) est repre par rapport la base
base, elle mme mobile avec une
la base du ressort
acclraJon
La vitesse et lacclra&on sont absolues
La vitesse et le dplacement lis au ressort/amor&sseur sont rela&fs
F. Ollivier S. Rgnier 16/17
page 55
Vibra&ons
K, C
page 56
Vibra&ons
Acclromtre
c
m
k
!
! (t)
!
x(t)
page 57
Vibra&ons
M
er
page 58
Vibra&ons
page 59
Vibra&ons
t
t 0 t
x(t) = Ae + Be e
amorJssement criJque : = 1
x(t) = [ x0 (1 + 0t) + v0 t ] e
0 t
x(t) = X e0 t cos(dt )
Pseudo pulsaJon
d = 0 1 2
page 60
Vibra&ons
m
k
!
x
Excita2on harmonique
(nota6on complexe) :
m
x(t) + cx (t) + kx(t) = F(t)
F(t) = F0 cos(t) = F0 Re e jt
page 61
Vibra&ons
m
k
!
x
m
x(t) + kx(t) = F0 cos(t) (1)
ou
x t (t) = X t cos( 0 t t )
x p (t) = Xp cos(t)
En la subsJtuant dans (1) on a lamplitude : X =
p
Le mouvement forc permanent scrit :
F0
k m
xp (t) =
F0 m
02 2
F0 m
cos ( t )
2
2
0
F. Ollivier S. Rgnier 16/17
page 62
Vibra&ons
CondiJons iniJales
et la soluJon complte (trans. + perm.) :
A = x0
F0 k
1 ( 0 )2
v0
B=
0
En ralit lamorJssement ntant pas nul - on nest plus alors avec un systme
conservaJf - la rponse due aux condiJons iniJales disparat au bout dun
certain temps. Pour cete raison, on ne sintresse dans ce qui suit quau
rgime tabli, qui correspond la parJe de lexcitaJon uniquement due
lexcitaJon extrieure.
F. Ollivier S. Rgnier 16/17
page 63
Vibra&ons
Rgime tabli
F0 m
F0
F0
1
x p (t) = 2
cos t =
cos t = H() cos t
2
2
k 1 ( 0 )
k
0
( )
( )
( )
x stat
H()
Frquence adimensionnelle
Frquence adimensionnelle
F. Ollivier S. Rgnier 16/17
page 64
Vibra&ons
Rgime tabli
F0
F0
1
x p (t) =
cos t = H() cos t
2
k 1 ( 0 )
k
( )
( )
page 65
Vibra&ons
Rgime tabli
F0
F0
1
x p (t) =
cos t = H() cos t
2
k 1 ( 0 )
k
( )
( )
page 66
Vibra&ons
F(t)
m
k
!
x
m
x(t) + cx (t) + kx(t) = F0 cos(t)
x(t) + 2 0 x (t) + 02 x(t) =
F0
cos(t)
m
(2)
x t (t) = X t e0 t cos(dt t )
xp (t) = Xp cos( t p )
La soluJon complte :
x(t) = [ Xt cos(dt t )] e
Rponse transitoire
0 t
+ Xp cos ( t p )
Rponse permanente
F. Ollivier S. Rgnier 16/17
page 67
Vibra&ons
m
x(t) + cx (t) + kx(t) = F0 cos(t)
F(t)
m
Xp () =
=
0
F0
m
F0
cos(t)
m
2 0
p () = Arctg 2
2
x(t) + 2 0 x (t) + 02 x(t) =
!
x
( ) + 4 0
2
0
F0
Xp () =
k
La soluJon complte :
2 2
(1 ) + 4
2
2 2
2
p () = Arctg
1 2
0 t
x(t) = X t cos( dt t ) e
+ Xp cos t p
Rponse transitoire
Rponse permanente
F. Ollivier S. Rgnier 16/17
page 68
Vibra&ons
F0
cos(t)
m
)
F0
cos t
m
(1) 2 (X1 cos t + X2 sin t) + 20( X1 sin t + X2 cos t) + 0 (X1 cos(t) + X2 sin t) =
F0
2
2
X
+
1
X
=
2
1
soit :
k
1 2 X 2 X = 0
2
1
X2 =
F0
cos t
k
) = Fm cos t
0
F0
2
2
2
X
et
2
+
1
X
=
1
1
1 2
1 2
k
F0
1 2
X1 =
k 4 2 2 + 1 2
et
X2 =
F0
2
k 4 2 2 + 1 2
page 69
Vibra&ons
F0
cos(t)
m
F0
1! " 2
F0
2#"
X1 =
et
X
=
2
k 4# 2 " 2 + (1! " 2 )2
k 4# 2 " 2 + (1! " 2 )2
Pour avoir lamplitude et la phase de x(t), il faut lidentifier la forme
x( t ) = X cos(t + )
X = X1 + X 2
tan =
F0
F0 1 2 + 4 2 2
X
=
=
k
k 4 2 2 + 1 2 2
X2 2
=
2
X1 1
(1 )
2 2
+ 4 2 2
2
= Arctg
2
1
page 70
Vibra&ons
F0
cos(t)
m
x(t) = Xcos ( t + )
!!
x(t) = !"2 X cos "t + #
et
F0
cos t
m
F
1 2 cos ( t + ) 2 sin ( t + ) X = 0 cos t
k
F
1 2 (cos t cos sin t sin ) 2 ( sin t cos + sin cos t ) X = 0 cos t
k
F
cos t 1 2 cos 2 sin sin t 1 2 sin + 2 cos X = 0 cos t
k
F0
cos
sin
X
=
(
)
#2$%
k ! tg" =
2
1#
%
2
(1 ) sin + 2 cos = 0
((
((
page 71
Vibra&ons
F0
cos(t)
m
+
%
2$# ( F0
2
X
cos!
1"
#
+
2$#
=
'
2*
1" # ) k
&
,
2
%
( F0
1"
#
- X sin! " 1" # 2
" 2$# * =
'
2$#
&
) k
.
F0
1 2
Xcos =
2
2
k
1 + 422
2
2
do X = ( Xcos ) + ( X sin )
2
X sin = F0
k 1 2 2 + 422
F0
X=
k
1! " 2 + 4# 2" 2
et
2
= Arctg
2
1
page 72
Vibra&ons
F0
cos(t)
m
(1) + 2j 0 + Xe
2
0
1 + 2j Xe
2
j t
F0
=
Xe j
k
F
do X = 0
k
(1 )
2
F0 jt
F0
2
2
j
= e + 2j 0 + 0 Xe =
m
m
F0
F0 1 2 2j
1
j
=
Xe =
k 1 2 + 2j
k 1 2 2 + 4 22
et
+ 422
2
= Arctg
1 2
page 73
Vibra&ons
Exemples de rponses
des excita6ons
harmoniques
x(t) = Xt cos( dt t )e
cos t p
F0
+
m ( 2 2 )2 + 2 2
( 0 )
0
page 74
Vibra&ons
F0
Xp =
m
( )
( ) + 2 0
2
0
2 2
et
2 0
p = Arctg 2
2
0
( )
k X ( )
H ( ) = H =
= X ( ) =
0
F0
X st
H ( ) =
1
(1 2 )2 + 422
page 75
Vibra&ons
X/Xst
3.5
3
2.5
X ( )
1
H ( ) =
=
X st ( = F0 k )
(1 2 )2 + 422
0
0.1
0.15
0.25
0.5
0.707
1
2
180
160
140
120
100
80
2
tg ( ) =
2
1
60
1.5
40
20
0
0.5
1.75
1.5
2.5
0.5
0
0
0.25
0.5
0.75
1.25
1.5
2.25
2.5
2.75
3
page 76
Vibra&ons
H ( ) =
3.5
1
(1 ) + 4
2 2
2 2
Mthode de
lamplica2on maximum
0.15
0.25
3
2.5
2
Q = Hmax H (1) =
1.5
1
On en dduit
0.5
1
2
1
2Q
Noter que :
0
0
0.25
0.5
0.75
1.25
1.5
1.75
Q=
mk
1 m 0
=
=
2
c
c
page 77
Vibra&ons
H ( ) =
3.5
3
1
(1 ) + 4
2 2
3dB
Mthode de la
bande passante -3dB
2 2
0.15
0.25
Hmax
2
H2 ( 3dB ) =
2.5
1
(1 3dB2 )2 + 42 3dB2
3dB2 = 1 22 4 1 2
3dB 1 si << 1
1.5
1
2 1 2 1 3dB
=
=
2
20
20
0.5
b1
0
0
0.25
0.5
0.75
b2
1
1.25
Noter que :
1.5
1.75
Q=
0
1
=
2 3dB
page 78
Vibra&ons
Admittance
Impdance
Exemple pour = 0.05 k=10
Admitance et Impdance
X 1
1
1
= H =
F k
k (1 2 )2 + 4 22
()
AdmiZance ou mobilit :
suscepJbilit au mouvement sous
leet dune force :
X
M() = ()
F
Impdance :
rsistance au mouvement
sous leet dune force
Z() =
F
()
X
page 79
Vibra&ons
Admitance et Impdance
Aux basses frquences
X 1
<< 0 : M() =
F k
La mobilit du systme dpend de la raideur
Autour de la rsonance
X
1
= 0 : M() =
=
F 2k
La mobilit du systme dpend de lamor&ssement
En haute frquence
02
X
1
>> 0 : M =
2
2
F
k
m
La mobilit du systme dpend de la masse
F. Ollivier S. Rgnier 16/17
page 80
Vibra&ons
m
x(t) = c( x (t) x s (t)) k(x(t) x s (t))
( )
x s (t) = X s cos t
m
x(t) + c( x (t) x s (t)) + k(x(t) x s (t)) = 0
x(t) + 2 0 x (t) + 02 x(t) = 2 0 x s (t) + 02 x s (t)
A
tg =
B
C = A 2 + B2
page 81
Vibra&ons
Donc ici
(1 ) + 4
2
X=
2 2
F0
cos(t)
m
F0
= Xs 02 4 2 2 + 1
m
mX s 02 1+ 4 2 2
k
x(t) = X cos(t )
+ 4
2 2
= Xs
1+ 4 2 2
+ 4 2 2
page 82
Vibra&ons
Transmission
intgrale
Isolation
xs(t)
c
On value le transfert
vibratoire du support vers
lobjet par la
Transmission
amplifie
Transmissibilit :
X
T=
=
Xs
1 + 4 2 2
+ 4 2 2
page 83
Vibra&ons
IsolaJon vibratoire
M = 1 kg et k = 105 N/m
soit
o = 300 rad/s ou fo= 50 Hz
avec = 0.1
k/200
et
*5
page 84
Vibra&ons
F1(t)
x1(t)
F2(t)
SL
x2(t)
.. Fn(t)
SL
xn(t)
a1 F1(t)
a2 F2(t)
a3 F3(t)
.
SL
....
an Fn(t)
F. Ollivier S. Rgnier 16/17
page 85
Vibra&ons
F0
t
Une excitaJon de priode T (telle que F(t+nT) = F(t)) peut scrire comme une srie :
F(t) =
F (t) = A + F
n=0
n=1
cos nt n
page 86
Vibra&ons
F0
t
Une excitaJon de priode T (telle que F(t+nT) = F(t)) peut scrire comme une srie :
F(t) =
F (t) = A + F
n
n=0
n=1
n=1
Coecients de Fourier
2
An = F(t)cos(nt)dt
T0
et
Valeur moyenne
A0 =
1
F(t) dt
T0
2
Bn = F(t)sin(nt)dt
T0
Fn = An + Bn
2
Bn
tg n =
An
F. Ollivier S. Rgnier 16/17
page 87
Vibra&ons
F0
t
Une excitaJon de priode T (telle que F(t+nT) = F(t)) peut scrire comme une srie :
F(t) =
F (t) = A + F
n=0
n=1
cos nt n
2
n Bn = F(t)sin(nt)dt= 0
T0
A0 = 0
F. Ollivier S. Rgnier 16/17
page 88
Vibra&ons
A0 = 0
n An = 0
2
=
T
T/2
T
1
T/2
T
2
cos(nt) + cos(nt)
Bn= sin(nt)dt sin(nt)dt =
0
T/2
T 0
n
T/2
1
2
1 cos(n)
Bn =
cos(n) + 1+ cos(2n) cos(n) =
n
n
0 si n est pair
Bn = 4
si n est impair
n
F. Ollivier S. Rgnier 16/17
page 89
Vibra&ons
A0 = 0
B2p = 0 B2p+1=
n An = 0
2
=
T
4
(2p + 1)
4
sin(3t) sin(5t)
f(t)= sin(t) +
+
+ ... = B2p+1 sin((2p + 1)t)
3
5
p=0
page 90
Vibra&ons
page 91
Vibra&ons
x(t) = Xp cos t p
m
x(t) + cx (t) + kx(t) = F0 cos(t)
F0
Xp () =
m
F0
Xp () =
k
( ) + 4 0
2
0
2 2
+ 4 2 2
2 0
p () = Arctg 2
2
0
2
p () = Arctg
1 2
F. Ollivier S. Rgnier 16/17
page 92
Vibra&ons
Fn
Avec Xn (n) =
m
et
( (n) ) + 4 n 0
2
0
2 2
2n
n (n) = Arctg
1 (n)2
page 93
Vibra&ons
et
Fn
Xn (n) =
k
(1 (n) ) + 4 n
2
2 2
2 0n
n (n) = Arctg 2
2
(n)
0
page 94
Vibra&ons
n =0
n =0
page 95
Vibra&ons
4
B2p+1=
(2p + 1)
( )
f1(t) = B1 sin t
B1
X1() =
m
x1(t) = X1 sin t 1
( () ) + 4 0
2
0
2 2
2 0
1() = Arctg 2
2
page 96
Vibra&ons
4
B2p+1=
(2p + 1)
f3 (t) = B3 sin 3t
B3
X 3 (3) =
m
x 3 (t) = X 3 sin 3t 3
( (3) ) + 4 3 0
2
0
2 2
2 0 3
3 (3) = Arctg 2
2
(3)
0
page 97
Vibra&ons
4
B2p+1=
(2p + 1)
f5 (t) = B5 sin 5t
B5
X 3 (5) =
m
x 5 (t) = X 5 sin 5t 5
( (5) ) + 4 5 0
2
0
2 2
2 0 5
5 (5) = Arctg 2
2
0 (5)
F. Ollivier S. Rgnier 16/17
page 98
Vibra&ons
4
B2p+1=
(2p + 1)
x(t) =
p=0
X
p=0
X 2p+1((2p + 1)) =
2p+1
B2p+1
2 2
2 0 (2p + 1)
2p+1((2p + 1)) = Arctg 2
2
((2p
+
1))
0
page 99
Vibra&ons
4
B2p+1=
(2p + 1)
x(t) =
X
p=0
X 2p+1((2p + 1)) =
2p+1
B2p+1
2 2
2(2p + 1)
2p+1((2p + 1)) = Arctg
1 ((2p + 1))2
F. Ollivier S. Rgnier 16/17
page 100
Vibra&ons
x(t)
F(t)
m
!
x
t
W= 0.6 w0 et x = 0.05
X1 = 1.5 F1/k
X2 = 2.2 F1/k
X3 = 0.5 F1/k
F1
wo
2W
3W
X2
X1
X3
W
wo
2W
3W
t
F. Ollivier S. Rgnier 16/17
page 101
Vibra&ons
f(t)e st dt
page 102
Vibra&ons
f(t)e st dt
Proprits uJles :
Thorme de Linarit
f(t t 0 )
e
ts0
f(t)
F(s)e
st 0
F(s + s0 )
page 103
Vibra&ons
sin t
Echelon
s2 + 2
(t)
1
s
cos t
Impulsion
s2 + 2
1
t
2
sinh t
s
2
2
s
1
at
e
s
sa
cosh t
2
2
n!
s
tn
n+1
s
( )
( )
( )
( )
page 104
Vibra&ons
F(t)
m
!
x
m
x(t) + cx (t) + kx(t) = F(t)
F(t) F(s);
x(t) X(s); x (t) sX(s) x 0 ;
x(t) s2 X(s) sx 0 v 0
s2 + 20 s + 0 2 X(s) (s + 20 )x 0 v 0 = F(s) m
F. Ollivier S. Rgnier 16/17
page 105
Vibra&ons
(s + 2 0 )x 0 + v 0
F(s)
1
X(s) =
+ 2
2
2
m s + 2 0s + 0 s + 2 0s + 02
s2 + 20 s + 0 2
(s + 2 0 )x 0 + v 0
F(s)
1
X(s) =
+
m s s1 s s2
s s1 s s2
)(
Rponsepermanente Xp (s)
) (
)(
X(s) x(t)
F. Ollivier S. Rgnier 16/17
page 106
Vibra&ons
F(s)
1
X(s) =
m s2 + 2 0s + 02
s2 + 20 s + 0 2
F(s)
1
X(s) =
avec s1,2 = 0 0 2 1
m s s1 s s2
)(
X(s) x(t)
F. Ollivier S. Rgnier 16/17
page 107
Vibra&ons
F0
1
X(s) =
avec s1,2 = 0 0 2 1
m s s1 s s2
)(
s1 s2 = 2 0 2 1
(
(
si > 1
s2 + s1 = 2 0
F0 1/ s1 s2 1/ s2 s1
X(s) = (
+
)
m
s s1
s s2
1
e
sa
at
()
F0
1
st
s t
x(t) =
e 1 e 2
m s1 s2
page 108
Vibra&ons
()
F0
1
st
s t
x(t) =
e 1 e 2
m s1 s2
s1 s2 = 2 0 2 1
F0
1
t
(
x(t) =
e 0 e 0
m 2 2 1
0
F0
x(t) =
m
1
1
2
= o 1
0 t
2 1)t
( 0 2 1)t
sinh( 0 2 1)t
Fo 0t
x(t) =
e
sinh(t)
m
F. Ollivier S. Rgnier 16/17
page 109
Vibra&ons
ts0
sinh t 2
s 2
( )
f(t) F(s + s0 )
F0
x(t) =
m
1
1
2
0 t
sinh( 0 2 1)t
page 110
Vibra&ons
F0
1
X(s) =
avec s1,2 = 0 j 0 1 2
m s s1 s s2
)(
s1 s2 = 2j 0 2 1
(
(
si 1
s2 + s1 = 2 0
F0 1/ s1 s2 1/ s2 s1
X(s) = (
+
)
m
s s1
s s2
()
F0
1
st
s t
x(t) =
e 1 e 2
m s1 s2
)
F. Ollivier S. Rgnier 16/17
page 111
Vibra&ons
()
F0
1
st
s t
x(t) =
e 1 e 2
m s1 s2
s2 s1 = 2j 0 1 2
F0
1
t
j(
x(t) =
e 0 e 0
m 2j 1 2
0
F0
x(t) =
m
1
0
d = o 1 2
0 t
2 1)t
j( 0 2 1)t
sin( 0 1 2 )t
Fo 0t
x(t) =
e
sin( dt)
m d
F. Ollivier S. Rgnier 16/17
page 112
Vibra&ons
Fo
1
m ( s s1 )( s s2 )
avec
s1,2 = 0 0 2 1
d = o 1 2
Fo t
e 0 sin( dt) si < 1
m d
F
t
x(t) = o .t.e 0
si = 1
m
Fo 0t
sh(t)
si > 1
m e
2
= o 1
m = 10 kg
k = 105N/m
Fo = 1000 N
page 113
Vibra&ons
A
O
page 114
Vibra&ons
u(t) est un chelon et marque le fait que lorsque t<0 le signal est nul et
est sans ralit physique
F. Ollivier S. Rgnier 16/17
page 115
Vibra&ons
page 116
Vibra&ons
u(t)
1
s
( )
cos t
f(t t 0 )
F(s)e
st 0
s
s2 + 2
page 117
Vibra&ons
page 118
Vibra&ons
(Em ) = Pe
(Ep +Ec ) = Pe
page 119
Vibra&ons
page 120
Vibra&ons
Equa&ons de Lagrange
Les quaJons de Lagrange sont un ouJl trs praJque pour obtenir les
quaJons de mouvement notamment dans le cas de systmes plusieurs
degrs de libert
Il sagit de lexpression du PFD sous une forme nergJque.
page 121
Vibra&ons
Le premier terme Ep(0) de ce potenJel est une constante quon peut la choisir
arbitrairement. On le prend en gnral 0. Si la posiJon dquilibre considre est
stable, le second terme est nul.
page 122
Vibra&ons
page 123
Vibra&ons
P
Qi =
q! i
F. Ollivier S. Rgnier 16/17
page 124
Vibra&ons
Ed
Dans le cas du frotement visqueux, ce potenJel de dissipaJon scrit
sous une forme quadraJque des vitesses :
Ed
Dans ce cas, la matrice C est appele matrice damorJssement ; elle est
symtrique dnie posiJve.
page 125
Vibra&ons
(en transla6on)
2 mthodes
Conserva2on de l nergie (Lagrange) :
1 seul degr de libert q=x
Bilan des forces (PFD) :
2
1
nergie cinJque : Ec =
2 mx
!
Force d excitaJon :
F(t)
nergie potenJelle
Ep = 12 kx2
Force de rappel lasJque :
-k x(t)
lasJque :
1 2
Force de frotement visqueux : -c x! ( t )
E
=
c x!
Energie dissipe :
d
Forces = m
Puissance extrieure :
F kx cx! = m!x!
c
!
F
m
k
!
x
2
Pe = F x
d L L Ed Pe
+
=
dt x x x x
!
mx +kx + cx = F
m
x
+
c
x
+kx
=F
!
page 126
Vibra&ons
F2
M2
M1
k
1- Modlisa&on
k
k
= IdenJcaJon :
- des 2 ddl : x1, x2
= q1,q2 = Coordonnes gnralises
- des iner&es
= Mij = Masses, Mnts dinerJe
- des raideurs
= Kij
- des forces
= F1, F2 = Eorts gnralises
2- Ecriture des 2 Equa&ons du mouvement
- Par le PFD
souvent complexe
- Par la conservaJon de lnergie
1
d
dt
(Ec + Ep ) = Pe
Remarque importante :
M11!!
q1 + K11q1 + M12!!
q2 + K12q2 = F1
q2 + K22q2 + M21!!
q1 + K21q1 = F2
M22!!
M21 = M12
et
K21 = K12
F. Ollivier S. Rgnier 16/17
page 127
Vibra&ons
!
1 2 1 2
1 2 1
Ec = mx 1 + mx 2
Ep = kx1 + k(x2 x1 )2
2
2
2
2
!
!
F2
m
page 128
Vibra&ons
L = Ec Ep
!
1 2 1 2 1 2 1
2
mx1 + mx2 kx1 k(x2 x1 )
2
2
2
2
L L
=
=
q1 x1
!
x1
d L L Pe
=
dt q i q q i
!
d L
= mx1 !!
= mx1
dt x1
1 2 1 2 1 2 1
2
mx 1 + mx 2 kx1 k(x2 x1 )
2
2
2
2
Libre
x1
P
Q1 =
= F1(t)
x! 1
= kx1 +k(x2 x1 )
Force
page 129
Vibra&ons
L L
=
=
q x2
! 2
1 2 1 2 1 2 1
2
mx1 + mx2 kx1 k(x2 x1 )
2
2
2
2
L L
=
=
q2 x2
!
Libre
x2
dt q i
!
P
Q2 =
= F2 (t)
x! 2
q i
d L
= mx2 !!
= mx2
dt x2
1 2 1 2 1 2 1
2
mx 1 + mx 2 kx1 k(x2 x1 )
2
2
2
2
x2
= k(x2 x1 )
Force
page 130
Vibra&ons
F2
m
M11
q1 +K11q1 +M12
q2 +K12q2 = F1
M22
q2 +K22q2 +M21
q1 +K21q1 = F2
!
M21 = M12
et
K21 = K12
F. Ollivier S. Rgnier 16/17
page 131
Vibra&ons
M11
q1 +K11q1 +M12
q2 +K12q2 = 0
M22
q2 +K22q2 +M21
q1 +K21q1 = 0
qi (t) = Qiert
(Q1, Q2 et r ?)
donnent 1 systme
dquations
caractristiques :
ra = ja
et
rb = jb
page 132
Vibra&ons
(
( e
) (
) + Q ( e
ja t
+ ae jat
b
2
jb t
+ be jbt
(
) + Q1bBcos (bt )
q2 (t) = Qa2 A cos ( a t ) + Qb2Bcos ( b t )
q1(t) = Q1a A cos a t
page 133
Vibra&ons
(
)
(
)
q2 (t) = Qa2 A cos ( a t ) + Qb2Bcos ( b t )
q1(t) = Q1a A cos a t + Q1bBcos b t
Q
+
K
Q
=0
(
)
(
)
12
12 a,b
1
22
22 a,b
2
Comme Det(S) = 0, les deux quaJons sont quivalentes et
donnent
2
2
K
a,b
a,b
a,b
11 a,b
12
12 a,b
Q 2a,b = 11
Q
=
Q
=
Q
1
1
a,b 1
2
2
K12 M12a,b
K 22 M22a,b
F. Ollivier S. Rgnier 16/17
page 134
Vibra&ons
(
) Bcos (bt )
q2 (t) = a A cos ( a t ) + b Bcos ( b t )
q1(t) =
A cos a t +
page 135
Vibra&ons
F1
k
F2
m
) ( )
( ) (
)
mr2 + 2k X + k X = 0
1
2
rt (S) =
x i (t) = Xie
2
k
X
+
mr
+k X2 = 0
1
!
)(
) ( )
3 5 k
ra = j a = j
2 m
!
3+ 5 k
rb = j b = j
2 m
!
F. Ollivier S. Rgnier 16/17
page 136
Vibra&ons
F2
m
3 5 k
3+ 5 k
a =
b =
2 m !
2 m
!
(
) + X1bBcos (bt )
x 2 (t) = Xa2 A cos ( a t ) + Xb2Bcos ( b t )
x1(t) = X1a A cos a t
) ( )
( ) (
)
2k m 2 X a + k X a = 0
a
1
2
a
(S) =
a
2
a
k
X
+
k
m
X
=0
1
a
2
!
2
2
2k
m
k
a
a
a a
X2a =
X
=
X
=
X
1
a 1
k2
k m 2a 1
F. Ollivier S. Rgnier 16/17
page 137
Vibra&ons
F2
m
3 5 k
3+ 5 k
a =
b =
2 m !
2 m
!
) ( )
( ) (
)
2k m 2 Xb + k Xb = 0
b
1
2
(S)b =
b
2
b
k X1 + k m b X2 = 0
!
2
2
2k
m
k
b
b
b b
On choisit arbitrairement:
Xb2 =
X
=
X
=
X
1
b 1
k2
k m b2 1
Xb1 = X1a = 1
!
!
(
) Bcos (bt )
x 2 (t) = a A cos ( a t ) + b Bcos ( b t )
x1(t) =
A cos a t +
page 138
Vibra&ons
F1
k
F2
m
3 5 k
3+ 5 k
a =
b =
2 m !
2 m
!
(
) Bcos (bt )
x 2 (t) = a A cos ( a t ) + b Bcos ( b t )
x1(t) =
A cos a t +
page 139
Vibra&ons
M22
q2 +K22q2 +M21
q1 +K21q1 = F2
!
M
K
q F
M
q
K
12
12
11
1 + 11
1 = 1
M12 M22
q2 K12 K22 q2 F2
!
M
K
q
F
q + Kq = F
! M
= Matrice dinerJe ou de masse
= Matrice de raideur
= Vecteur des coordonnes gnralises
= Vecteur des eorts gnralises
page 140
Vibra&ons
M22
q2 +K22q2 +M21
q1 +K21q1 = F2
!
M
K
q F
M
q
K
12
12
11
1 + 11
1 = 1
M12 M22
q2 K12 K22 q2 F2
!
m 0 x1 2k k x1 F1
0 m x + k k x = F
2
2 2
!
F. Ollivier S. Rgnier 16/17
page 141
Vibra&ons
Ep = 12 Kijqiqj
i,j
e = Fi q i
i
page 142
Vibra&ons
1 2 1 2
1 2 1
Ec = mx 1 + mx 2 Ep = kx1 + k(x2 x1 )2
2
2
2
2
!
!
M
K
M
q
K12
11
12
1
11
+
M12 M22
q2 K12 K22
!
q F
1 = 1
q2 F2
m 0 x1 2k k x1 F1
0 m x + k k x = F
2
2 2
!
F1
k
F2
m
page 143
Vibra&ons
jt
+ Kq = 0
LquaJon matricielle du mouvement libre scrit : !Mq
(S)
Q1
(K M ) Q = 0
2
2
Det (K M2 ) =
K22 M22
=0
page 144
Vibra&ons
(
)
(
)
q2 (t) = Qa2 A cos ( a t ) + Qb2Bcos ( b t )
Q
2
1
(K Ma ) Q a = 0
2
On obtient
1
Q
a
a
et
Q
2
1
(K Mb ) Q b = 0
2
et
1
Q
b
b
page 145
Vibra&ons
F1
m 0 x1 2k k x1 F1
0 m x + k k x = F
2
2 2
!
q + Kq = 0
LquaJon matricielle du mouvement libre scrit : !M
X
(S)
K M 2 1 = 0
x i (t) = Xie jt
X2
!
2
2k
m
k
2
2
2
Det K M =
=
(2k
m
)(k
m
)k
=0
2
k
k m
!
On obJent une quaJon du 2nd degr en w4 : 4 racines wa et wb
2
3 5 k
3+ 5 k
a =
b =
2 m !
2 m
!
F. Ollivier S. Rgnier 16/17
page 146
Vibra&ons
F2
m
(
) + X1bBcos (bt )
x 2 (t) = Xa2 A cos ( a t ) + Xb2Bcos ( b t )
x1(t) = X1a A cos a t
(
)(
)
Pour dterminer les modes propres on rsout :
Xb
Xa
K M 2a 1a = 0
K M b2 1b = 0
et
X2
X2
!
On obtient
1
1
a
b
et
X
X
a
b
m2 2a 2 b2 2
page 147
Vibra&ons
M
! q + Kq = F
F1 jt
Avec par exemple : F = Re e W : pulsaJon dexcitaJon
0
Q1 jt
q ( t ) = Re e
Q2
Q
Q1 jt F1 jt
F1
2 1
(K -M ) Q e = 0 e (K -M ) Q = 0
2
2
2
page 148
Vibra&ons
F K M 2 K M 2 Q F
1
1
11
12
12
= 1 11
=
2
2
0 K12 M12 K22 M22 Q 2 0
F1 K22 M222
Q 1 () =
et
Q 2 () =
Det(S())
!
Les racines de Det(S) sont wa et wb, do :
Q 1 ()
!
F1 K22 M222
2
a
)(
2
b
!!Q 2 ()
F1 K12 M122
Det(S())
F1 K12 M122
)(
2a 2 b2 2
F. Ollivier S. Rgnier 16/17
)
page 149
Vibra&ons
F1
k
F2
X F
X F
2
k
1
1
K #M!2 1 = 1 2k m
=
2
X2 0
k
k m X2 0
!
!
k
1
Det
2
F1 (k m2 )
0 k m
X1 () =
=
2
2
2
2k m2
(2k
m
)(k
m
)k
k
Det
2
k
k
m
!
2k m2 F
1
Det
F1k
k
0
X2 () =
=
2
2
2
2k m2
(2k
m
)(k
m
)k
k
Det
2
k
k
m
F. Ollivier S. Rgnier 16/17
!
page 150
Vibra&ons
F1
On a donc :
)(
X1 () =
!
X2 () =
!
F1 (k m2 )
)(
m2 2a 2 b2 2
F1k
)(
m2 2a 2 b2 2
)
F. Ollivier S. Rgnier 16/17
page 151
Vibra&ons
1 2 1 2
Ec = mx 1 + mx 2
2
2
!
kc
M
k
P = F1(t)x1(t)
1 2 1
1
2
Ep = kx1 + k(x2 x1 ) + k(x2 )2 = kx12 +k(x2 x1 )+kx22
2
2
2
!
m 0 x1 2k k x1 F1
0 m x + k 2k x = 0
2
2
!
F. Ollivier S. Rgnier 16/17
page 152
Vibra&ons
X1 () =
X2 () =
!
m 0 x1 2k k x1 F1
0 m x + k 2k x = 0
2
2
kc
F1 K22 M222
Det(S())
)=
F1 K12 M122
Det(S())
F1 (2k m2 )
)=
2
a
)( )
2
b
F1k
)(
m2 2a 2 b2 2
page 153
Vibra&ons
X1 () =
X2 () =
!
kc
F1 (2k m2 )
m
2
2
a
)( )
2
b
F1k
)(
m2 2a 2 b2 2
)
F. Ollivier S. Rgnier 16/17
page 154
Vibra&ons
kc
page 155
Vibra&ons
k1
k2
m1
m2
x2(t)
x1(t)
F. Ollivier S. Rgnier 16/17
page 156
Vibra&ons
k1
m1
x1(t)
1
1
H1(f) =
m1 f 2 f12
m1
k2
k1
m2
r
x
H1 '(f)
(f
f2 2 )
(f 2 fa 2 )(f 2 fb 2 )
page 157
Vibra&ons
Cordes
F
T
Barres et poutres
Ondes Transverses Flexion
Membranes
Ondes
de torsion
Plaques
F. Ollivier S. Rgnier 16/17
page 158
Vibra&ons
Exemples de problmes
Rail-Ballast
Modle continu
lment
du viaduc
de Millau
Y
V0t
V0
F0(x,t)
Modlisation
de la guitare acoustique
page 159
Vibra&ons
Programme
Mthodes exactes
structures simples
Ondes dans les cordes
Mthodes approches
pour lingnieur
analyse rapide
des structures composes
Mthode de Rayleigh :
approximation 1DDL
Mthode de Ritz :
approximation 2 ou 3 DDL
lments finis
Simulation des structures complexes
Principes de mise en uvre
Mthodes exprimentales
Analyse du signal vibratoire
Analyse modale
Expriences
tude des cordes vibrantes
tude des ondes de flexion
dans les poutres
tude dune plaque mince
F. Ollivier S. Rgnier 16/17
page 160