4.forced Vibrations
4.forced Vibrations
.. .
mx cx kx
k C
m
F Sin wt
m
x
Free Body diagrams
F Sinwt
Forced Vibration
b cw
F
Velocity x w A cos(wt ) w A sin wt and
kA 2
wt-w
wt x aw 2 A sin(
Acceleration t ) w 2 A sin wt
Substituting in (Reference
i), the equilibrium
axis equation becomes,
mw A sin wt + cw A sin wt + kA sin(wt ) = Fo sin wt
2
2
The above forces may be represented by vectors drawn in order such that
they form a closed polygon. Damping force is perpendicular to spring force
& the inertia force is perpendicular to damping force.
Vector representation From the triangle oab in the fig,
of forces F 2 (kA mw 2 A) 2 (cw A) 2
A2 (k mw 2 ) 2 (cw ) 2
F
A
A
mw (k mw 2 ) 2 (cw ) 2
b cw
F Taking k outside the sq root,
kA F
A k
wt a wt-
k
2
1 mw cw
2 2
k
Reference axis
Also k mwn2 , c 2 mwn
F
A k
2
w 2w 2
2
1 w w
n n
Steady state amplitude & phase angle
of forced vibrations :
A
mw
F
b cw
A
F k
kA
wt a wt- 1 r 2 r
2 2 2
2ζr
From the fig , phase angle between force & displacement is = tan 2
-1
1 - r
(i) At resonance, w = wn r 1 A
F
k , = tan -1 900
2
1 w wn
n
Ast A 1
A =
A
2 2 2
2 2 2
1 r 2 r st 1 - r + 2ζr
A
The dimensionless ratio is called the Magnification factor (MF)
Ast
& it is a function of frequency ratio & damping factor.
Plot of MF vs. r for different values of
Important observations of MF vs. Frequency ratio
• Any amount of damping ( >0) reduces the magnification factor (MF) for all
values of forcing frequency.
• The amplitude of vibration is infinite at resonant frequency & zero damping
factor.
• The peak amplitude for any amount of damping occurs slightly before r=1
• For >1/√2, the graphs of MF decreases with increasing values of r.
Frequency corresponding to peak amplitude
It is observed that the peak amplitude for any amount of damping occurs
slightly before r 1 or resonant frequency.
Ast
A . Max value of A occurs when the denominator
1 r 2 r
2 2 2
A
mw
b cw
F
kA
wt a wt-
Reference axis
1) For ζ >0 and 0<r<1 the phase angle is given by 00<<900, implying that the
response lags excitation.
2) For ζ >0 and r>1, the phase angle is given by 900<< 1800, implying that the
response leads excitation.
3) For ζ >0 and r=1, the phase angle is =900 implying that the phase difference
between the excitation and response is 900.
4) For ζ >0 and large values of r, the phase angle approaches 1800 implying that
the response and excitation are out of phase.
Numerical problem 1
A
F
k =
30 2039 0.0154 m
(1 r ) (2 r )
2 2 2
(1 0.88 ) (2 0.526 0.88)
2 2 2
A'
F
k =
30
2039 0.0652 m
(1 r 2 ) (1 0.882 )
Increase in amplitude upon removal of damper
( A ' A) (0.0652 0.0154) 0.05 m = 50 mm
Numerical Problem 2
A
F
k =
10012269.4 0.022 m
(1 r ) (2 r )
2 2 2
(1 1.318 ) (2 11.318)
2 2 2
Numerical Problem 3
2
1- r 1-1.58
-ve sign indicates that the force lags the displacement.
(iii) Amplitude at resonance : w = wn r 1
A
F 150
k 10000 0.0822 m = 82.2 mm
2 2 0.0912
Rotating Unbalance
• Unbalance in rotating machines is a common source of
vibration excitation.
• We consider here a spring (k) & mass (M) constrained to
move in the vertical direction and excited by a rotating
machine that is unbalanced, as shown in Fig.
Rotating Unbalance
• The unbalance is represented by an eccentric mass m
with eccentricity e that is rotating with angular velocity
w.
• By letting x be the displacement of the non rotating
mass (M - m) from the static equilibrium position, the
displacement of m is (x+esinwt)
Rotating Unbalance
l
l where m mass of piston,
stroke length
e wt e crank radius
2
l length of connecting rod ,
w angular velocity of crank .
M x
e
As e l , may be neglected,
l
harmonic force is (mw 2e) sin wt F = (mω2 e)
k/2 c k/2
cw
2
A k
(1 r 2 ) 2 2 r (1 r 2 ) 2 2 r (1 r 2 ) 2 2 r
2 2 2
1 + 2ζr
2
Amplitude of body A
Motion transmissibility = =
Amplitude of support B (1 - r ) + 2ζr
2 2 2
2 r
Phase angle between A & B = ( - ) = tan 2
tan 1
2 r
1
1 r
Relative Motion
mx k ( x y ) c( x y )
Relative motion of a mass is its motion w.r.t support. Let the absolute
motion of mass be x = A sinωt & the support motion be y = B sinωt.
Hence the relative motion is z ( x - y ) x ( z y )
The equation of motion can be written as; mx k ( x y ) c( x y ) 0
m( z y ) cz kz 0 Or mz + cz + kz = -my (i)
Substituting y B sin wt & y w 2 B sin wt in Eqn (i),
mz cz kz mw 2 B sin wt Comparing mz cz kz mw 2 B sin wt with
standard equation of forced vibration mz cz kz F sin wt ,
k 2 k w
we see that; F mw B 2 w B Bkr wn
2 2 2
and r
wn m wn
Also we know that the the steady state amplitude of forced vibration
F Bkr 2
k k Br 2
Z
(1 r 2 ) 2 2 r (1 r 2 ) 2 2 r (1 r 2 ) 2 2 r
2 2 2
2 r
Phase angle between Z & B = = tan 2
1
1 r
Vibration Isolation & Force Transmissibility
• Vibrations are produced in machines having unbalanced masses
& these are transmitted to the foundation upon which machines
are installed.
• To diminish these undesirable transmission of forces, machines
are usually mounted on vibration isolating material such as
rubber, felt cork, metallic springs. etc. which provide stiffness &
damping.
• Transmissibility (denoted by TR or e) is the ratio of the force
transmitted to the foundation to the force applied.
• It is a measure of effectiveness of the vibration isolating material.
Expression for Transmissibility Ratio (TR or e )
F
FT
cw
kA
Reference axis
As the transmitted force ( FT ) is the vector sum of the spring force (kA)
& the damper force (cw A) which are perpendicular to each other,
cw
2
F
F 1
2 r
2
FT k k 1 2 r
2
(1 r 2 ) 2 2 r 2
(1 r ) 2 r
2 2 2
Impressed Force
Force transmissibility ratio TR (or e ) =
Transmitted force
1 + 2ζr
2
FT
e =
F (1 - r 2 )2 + 2ζr
2
2 r
Phase angle between F & FT = ( - ) = tan 1 2
tan 1
2 r
1 r
Note : The expressions for motion & force transmissibility are same.
Transmissibility vs. frequency ratio for different values of damping ratio
Transmissibility vs. frequency ratio for different values of damping ratio
(i ) When r 2, e 1. i.e. the transmitted force
A 1 1
When 0, , As r 1, A 2
B (1 r )
2
(r 1)
A 1
Amplitude of main mass A 0.884 mm
5 (2.58 1)
2
(ii) When damping ratio ( ) 0.2 :
1 2 0.2 2.58
2
A
When 0.2,
B (1- 2.58 ) 2 0.2 2.58
2 2 2
A 1 1.0322
0.25
B 5.6562 1.0322
A
0.25 A 1.25 mm
5
Amplitude of main mass A 1.25 mm
Numerical Problem 2
(Based on support excitation)
Determine the critical speed when an automobile trailer is travelling
over a road with sinusoidal profile of wavelength 15 meters and an
amplitude of 75 mm. The springs of the automobile are compressed
0.125 m under its own weight. Also determine the amplitude of
vibration at 50 kmph.
x
m
k
y
V Kmph
15 m
Data : 0.125 m, B 75 mm , v 50 kmph, Wave length 15 m
Solution :
(i) Critical speed of trailer :
g 9.81
Natural frequency wn 8.86 rad / sec
0.125
2
For critical speed, w wn 8.86 f n 1.41 Hz
wn
We know that Linear velocity = frequency × wave length
i.e. v = f n λ = 1.41 15 = 21.15 m / sec
3600
To convert m / sec into kmph, multiply by
1000
3600
Critical Velocity of the automobile = 21.15 76.14 kmph
1000
(ii) Amplitude of vibration at 50 kmph :
1000
Velocity in m / sec v 50 13.89 m / sec
3600
We know that v = f λ 13.89 = f 15 f = 0.926 Hz
w 2 f 5.82 rad / sec. Also wn 46.67 rad / sec
w 5.82
Frequency ratio r 0.657
wn 8.86
A 1 A 1
When 0, , As r 1,
B (1 r )2
B (1 r 2 )
A 1
Amplitude of trailer A 132 mm
75 (1 0.657 )
2
Numerical Problem 3
(Based on support excitation)
A trailer has 1000 Kg mass when fully loaded and 250 Kg
when empty. The suspension has a stiffness of 350 kN/m.
The damping factor is 0.5. The speed of the trailer is 100
Km/hr. The road varies sinusoidally with a wave length of 5
m. Determine the amplitude ratio of the trailer:
1. When fully loaded.
2. When empty.
Data :
m 1000 kg ( full load ) & m 250kg ( Empty ), v 100 kmph, 0.5
Wavelength of road suface 5 m, k 300 KN / m 300 103 N / m
Solution :
1000
Velocity of the vehicle v = 100 kmph = 100 27.78 m / sec
3600
We know that Linear velocity (v) = frequency (f) × wave length ( )
i.e. 27.78 = f 5 f = 5.556 Hz
w 2 f 2 5.556 34.91 rad / sec
(i) Amplitude ratio when vehicle is fully loaded : (m = 1000 kg)
1 2 0.5164 1.352
2
A
When 0.5164,
B (1-1.3522 ) 2 2 0.5164 1.352
2
A 1.7175
1.058 Amplitude of main mass A 0.106 m m
0.1 1.623
(ii) Dynamic load on isolators : W .K .T Relative amplitude ratio
Z r2 1.3522
B (1- r 2 ) 2 2 r (1-1.352 2 ) 2 2 0.5164 1.352
2 2
Z
1.126 Z 1.126 0.1 0.1126 mm 0.1126 10 3 mm
B
Dynamic load on each isolator FD = Z k 2 (cw ) 2
FD = 0.1126 10 3 (30000) 2 (400 104.72) 2 5.8N
Numerical Problem 1
Force Transmissibility problems
A reciprocating machine of mass 75 Kg is mounted on springs of
stiffness 11.76x105 N/m and a damper of damping factor 0.2.
The slider of mass 2 Kg within the machine has a reciprocating
motion with a stroke of 0.08 m. The speed is 3000 rpm.
Assuming the motion of the piston to be harmonic.
1. Amplitude of vibration of the machine.
2. Transmissibility ratio.
3. Force transmitted to the foundation.
4. Is vibration isolation achieved? If so how?
Data : M 75 kg , m 2 kg , N 3000 rpm, 0.2, k 11.76 105 N / m
Stroke of piston 0.08 m crank radius e 0.04 m
Solution :
k 11.76 105
Natural frequency wn 125.22 rad / sec
M 75
2 N 2 3000
Forcing frequency w 314.16 rad / sec
60 60
w 314.16
Frequency ratio r 2.51
wn 125.22
Unbalanced force F = mw 2e 2 (314.16) 2 0.04 7896 N
F
(i) Amplitude of vibration : A k
(1- r 2 ) 2 2 r
2
7896
5
6.714 103
11.76 10
A 1.245 103 m
(1- 2.51 ) 2 0.2 2.51
2 2 2
28.09 1.008
Amplitude of vibration of machine A 1.245 m m
(ii) Transmissibility ratio :
1 2 r 1 2 0.2 2.51
2 2
TR
(1- r 2 ) 2 2 r (1- 2.512 ) 2 2 0.2 2.51
2 2
1.417
TR (or e ) 0.2627
5.394
(iii) Force transmitted to the foundation :
Force transmitted FT FT
As TR 0.2627
Force impressed F 7896
Force transmitted to the foundation FT 2074.3 N
TR
(1- r 2 ) 2 2 r (1- 4.582 ) 2 2 0.057 4.58
2 2
FT
TR (or e ) 0.0564 =
F
Force transmitted to foundation at 800 rpm
FT 0.0564 F 0.0564 772 43.55 N
(iii) Force transmitted at resonance with damping ( = 0.057) :
resonance, w wn r 1,
1 2 1 2 0.057
2 2
TR 8.83
2 2 0.057
Also, at resonance, Fres mwn2e 2.2 (18.28) 2 0.05 36.76 N
i.e. Force transmitted t FT 8.83 Fres 8.83 36.76 325 N
(iv) Amplitude of vibration at resonance : (r 1)
F Fres 36.76
A k Ares k 26740 0.012 m
(1- r ) 2 r
2 2 2 2 2 0.057