Vibration Control Techniques Overview
Vibration Control Techniques Overview
Yang, NTNU MT
Chapter 9
Vibration Control
9
台灣師範大學機電科技學系 -1-
C. R. Yang, NTNU MT
Chapter Outline
9.1 Introduction
9.2 Vibration Nomograph and Vibration Criteria
9.3 Reduction of Vibration at the source
9.4 Balancing of Rotating Machines
9.5 Whirling of Rotating Shafts
9.6 Balancing of Reciprocating Engines
9.7 Control of Vibration
9.8 Control of Natural Frequencies
9.9 Introduction of Damping
9.10 Vibration Isolation
9.11 Vibration Absorbers
台灣師範大學機電科技學系 -2-
C. R. Yang, NTNU MT
9.1
Introduction
9.1
台灣師範大學機電科技學系 -3-
C. R. Yang, NTNU MT
9.1 Introduction
台灣師範大學機電科技學系 -4-
C. R. Yang, NTNU MT
9.2
Vibration Nomograph and Vibration Criteria
台灣師範大學機電科技學系 -6-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -7-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -8-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -9-
C. R. Yang, NTNU MT
Example 9.1
Helicopter Seat Vibration Reduction
The seat of a helicopter, with the pilot, weights 1000N and is found to
have a static deflection of 10 mm under self-weight. The vibration of
the rotor is transmitted to the base of the seat as harmonic motion
with frequency 4 Hz and amplitude 0.2 mm.
a) What is the level of vibration felt by the pilot?
b) How can the seat be redesigned to reduce the effect of vibration?
台灣師範大學機電科技學系 -10-
C. R. Yang, NTNU MT
Example 9.1
Helicopter Seat Vibration Reduction
Solution
k 105
Natural frequency = ωn = 31.3209 rad/s 4.9849 Hz
m 101.9368
4.9849
Frequency ratio = r = 1.2462
n 4.0
Y
Amplitude of vibration felt by pilot: X
1 r 2
where Y is the amplitude of base displacement
台灣師範大學機電科技學系 -11-
C. R. Yang, NTNU MT
Example 9.1
Helicopter Seat Vibration Reduction
Solution
0. 2
X 2
0.3616 mm
1 1.2462
vmax 2fX 2 50.3616 9.0887 mm/s
amax 2f X 228.4074 mm/s 2 0.2284 m/s 2
2
台灣師範大學機電科技學系 -12-
C. R. Yang, NTNU MT
Example 9.1
Helicopter Seat Vibration Reduction
Solution
amax 10 mm/s 2 2f X 8
2 2
X 0.01583 mm
X 0.01583 1
2
or r 3.6923
Y 0.2 1 r
8 k
n 6.8068 rad/s
3.6923 3.6923 m
m 101.9368 kg k 4722.9837 N/m
台灣師範大學機電科技學系 -13-
C. R. Yang, NTNU MT
9.3
Reduction of Vibration at the Source
•
Some sources cannot be eliminated e.g. turbulence, engine
combustion instability, road roughness X
台灣師範大學機電科技學系 -15-
C. R. Yang, NTNU MT
9.4
Balancing of Rotating Machines
台灣師範大學機電科技學系 -17-
C. R. Yang, NTNU MT
• Single-Plane Balancing
Rotate the disc and let it come to rest. Mark the lowest point.
Repeat a few times.
台灣師範大學機電科技學系 -18-
C. R. Yang, NTNU MT
• Single-Plane Balancing
台灣師範大學機電科技學系 -19-
C. R. Yang, NTNU MT
• Single-Plane Balancing
a2 a
Measured bearing reactions: F1 mr 2 , F2 1 mr 2
l l
台灣師範大學機電科技學系 -20-
C. R. Yang, NTNU MT
• Single-Plane Balancing
台灣師範大學機電科技學系 -21-
C. R. Yang, NTNU MT
• Single-Plane Balancing
Turn the rotor and fire a stroboscopic light at the same frequency ω.
台灣師範大學機電科技學系 -22-
C. R. Yang, NTNU MT
• Single-Plane Balancing
Add a known trial weight W to the rotor and repeat the procedure.
The new angle of the marking, φ and the new amplitude Au+w are
noted.
台灣師範大學機電科技學系 -23-
C. R. Yang, NTNU MT
• Single-Plane Balancing
The difference vector AW AU W AU is the unbalance vector
due to trial weight W.
AW AU2 AU2 W 2 AU AU W cos
1 A 2
A 2
A 2
cos U U W U W
2 AU AW
Magnitude of original unbalance WO=(AU/AW)•W
台灣師範大學機電科技學系 -24-
C. R. Yang, NTNU MT
• Two-Plane Balancing
台灣師範大學機電科技學系 -25-
C. R. Yang, NTNU MT
• Two-Plane Balancing
台灣師範大學機電科技學系 -26-
C. R. Yang, NTNU MT
• Two-Plane Balancing
3
Thus m1 = m/3, m2 = 2m/3
Thus any unbalanced mass can be replaced by 2 unbalanced mass
at the end planes.
台灣師範大學機電科技學系 -27-
C. R. Yang, NTNU MT
• Two-Plane Balancing
Vibration analyzer
Replace unbalance weight by UL and UR as shown:
Measure vibration amplitude and phase VA and VB of original
unbalance at A and B
VA AALU L AARU R (9.13)
VB ABLU L ABRU R (9.14)
台灣師範大學機電科技學系 -28-
C. R. Yang, NTNU MT
• Two-Plane Balancing
Vibration analyzer
Add known trial weight WL in left plane at known position.
VA AAL U L WL AARU R
Eq 9.15
VB ABL U L WL ABRU R Eq 9.16
Subtract Eq. 9.13 and 9.14 from Eq. 9.15 and 9.16:
VA VA VB VB
AAL , ABL
WL WL
台灣師範大學機電科技學系 -29-
C. R. Yang, NTNU MT
• Two-Plane Balancing
Vibration analyzer
Remove WL and add known trial weight WR in right plane at known
position.
VA AAR U R WR AALU L
Eq 9.19
VB ABR U R WR ABLU L Eq 9.20
Subtract Eq. 9.13 and 9.14 from Eq. 9.19 and 9.20:
VA VA VB VB
AAR , ABR
WR WR
台灣師範大學機電科技學系 -30-
C. R. Yang, NTNU MT
• Two-Plane Balancing
analyzer
Vibration
Once Aij are known, Eq 9.13 and Eq 9.14 can be solved to find the
unbalance vectors.
ABRVA AARVB ABLVA AALVB
UL , UR
ABR AAL AAR ABL ABL AAR AAL ABR
Rotor can now be balanced by adding equal and opposite weights in
each plane.
BL U L , BR U R
台灣師範大學機電科技學系 -31-
C. R. Yang, NTNU MT
Example 9.2
Two-Plane Balancing of Turbine Rotor
台灣師範大學機電科技學系 -32-
C. R. Yang, NTNU MT
Example 9.2
Two-Plane Balancing of Turbine Rotor
Solution VA 0.260 0.1 i 0.1732
VB 0.16205 0.145 i 0.0676
VA 0.15125 0.086 i 0.1229
VB 0.1230 0.0643 i 0.0766
VA 0.1535 0.1229 i 0.086
VB 0.25160 0.2349 i 0.0855
WL 0.3270 0 i 0.3
WR 0.34180 0.34 i 0
台灣師範大學機電科技學系 -33-
C. R. Yang, NTNU MT
Example 9.2
Two-Plane Balancing of Turbine Rotor
Solution
VA VA 0.186 i 0.0503
AAL 0.1677 i 0.62
WL 0 i 0.3
VB VB 0.0807 i 0.009
ABL 0.03 i 0.269
WL 0 i 0.3
VA VA 0.0229 i 0.0872
AAR 0.0674 i 0.2565
WR 0.34 i 0
VB VB 0.0899 i 0.1531
ABR 0.2644 i 0.4503
WR 0.34 i 0
台灣師範大學機電科技學系 -34-
C. R. Yang, NTNU MT
Example 9.2
Two-Plane Balancing of Turbine Rotor
Solution
UL
0.1045 i0.0008 0.0271 i0.0327
0.2347 i0.2397 0.712 i0.0104
0.2565 i 0.1546
UR
0.0436 i 0.0321 0.663 i 0.0785
0.0712 i0.0104 0.2347 i0.2397
0.0875 i 0.1612
台灣師範大學機電科技學系 -35-
C. R. Yang, NTNU MT
Example 9.2
Two-Plane Balancing of Turbine Rotor
Solution
Thus the required balance weights are
BL U L 0.2565 i 0.1546 0.2994 31.0772
BR U R 0.0875 i 0.1612 0.1834 61.5242
台灣師範大學機電科技學系 -36-
C. R. Yang, NTNU MT
9.5
Whirling of Rotating Shafts
台灣師範大學機電科技學系 -38-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -39-
C. R. Yang, NTNU MT
• Equations of Motion
台灣師範大學機電科技學系 -40-
C. R. Yang, NTNU MT
• Equations of Motion
台灣師範大學機電科技學系 -41-
C. R. Yang, NTNU MT
• Equations of Motion
Substitute Fi , Fe , Fdi , Finto
de Eq. 9.26:
mx ci c x kx ciy m 2 a cos t Eq 9.32
my ci c y ky cix m 2 a sin t Eq 9.33
These equations are coupled.
Define w as w=x+iy
Add Eq 9.32 to Eq 9.33 and multiply by i:
台灣師範大學機電科技學系 -42-
C. R. Yang, NTNU MT
• Critical Speeds
k
n
m
When ω = ωn, rotor undergoes large deflections.
台灣師範大學機電科技學系 -43-
C. R. Yang, NTNU MT
Phase angle
c 1 2 r k c
tan 1 2
tan 2
r , n ,
k m 1 r n m 2 km
台灣師範大學機電科技學系 -44-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -45-
C. R. Yang, NTNU MT
• Stability Analysis
p2 = -m, p1 = q2 = 0, q1 = ci + c, p0 = k, q0 = -ωci
• Stability Analysis
k c
1 0
m ci
Internal and external friction can cause instability at rotating speeds
above the 1st critical speed.
台灣師範大學機電科技學系 -47-
C. R. Yang, NTNU MT
Example 9.3
Whirl Amplitude of a Shaft Carrying an Unbalanced Rotor
台灣師範大學機電科技學系 -48-
C. R. Yang, NTNU MT
Example 9.3
Whirl Amplitude of a Shaft Carrying an Unbalanced Rotor
Solution
12000 2
Forcing frequency of rotor: 400 1256.6371 rad/s
60
k 40 10 6
Natural frequency: n 894.4272 rad/s
m 50
1256.6371
Frequency ratio: r 1.4049
n 894.4272
台灣師範大學機電科技學系 -49-
C. R. Yang, NTNU MT
Example 9.3
Whirl Amplitude of a Shaft Carrying an Unbalanced Rotor
Solution
a) Steady-state amplitude: A
ar 2
0.0021.40492
1 r 2r
2 2 2
1 1.40492 2 2 0.11.40492
3.8950 10 3 m 3.8950mm
台灣師範大學機電科技學系 -50-
C. R. Yang, NTNU MT
9.6
Balancing of Reciprocating Engines
• Vibrations due to
Periodic variations of gas pressure in cylinder
Inertia forces associated with moving parts
台灣師範大學機電科技學系 -52-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -53-
C. R. Yang, NTNU MT
r cos
Resultant torque MQ=FhtanΦ where h
sin
Fr cos
MQ
cos
Torque induced at crank shaft is felt at engine support.
台灣師範大學機電科技學系 -54-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -55-
C. R. Yang, NTNU MT
Acceleration of Piston
xP r l r cos l cos
r l r cos t l 1 sin 2
But l sin r sin r sin t
r2
Hence cos 1 2 sin 2 t
l
Substituting into xP ,
r2
xP r l r cos t l 1 2 sin 2 t
l
台灣師範大學機電科技學系 -56-
C. R. Yang, NTNU MT
Acceleration of Piston
r 1
In general, and 1 1
l 4 2
r2
Hence xP r 1 cos t sin 2 t
2l
r r
or xP r 1 r cos t cos 2t
2l 4l
r
x P r sin t sin 2t
2l
r
xP r 2 cos t cos 2t
l
台灣師範大學機電科技學系 -57-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -58-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -60-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -61-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -62-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -63-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -64-
C. R. Yang, NTNU MT
9.7
Control of Vibration
9.7
台灣師範大學機電科技學系 -65-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -66-
C. R. Yang, NTNU MT
9.8
Control of Natural Frequencies
台灣師範大學機電科技學系 -68-
C. R. Yang, NTNU MT
9.9
Introduction of Damping
9.9
台灣師範大學機電科技學系 -69-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -70-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -71-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -72-
C. R. Yang, NTNU MT
9.10
Vibration Isolation
9.10
台灣師範大學機電科技學系 -73-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -74-
C. R. Yang, NTNU MT
impulsive forces
Protect system against motion of its foundation
台灣師範大學機電科技學系 -75-
C. R. Yang, NTNU MT
Fi t kxt cx t
Fi t mxt k xt y t c x t y t
台灣師範大學機電科技學系 -76-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -77-
C. R. Yang, NTNU MT
k m c
2 2 2 2
台灣師範大學機電科技學系 -78-
C. R. Yang, NTNU MT
FT k 2 2c 2
Transmissibility Tr
F0 k m c
2 2 2 2
1 2r
2
where r
1 r 2r
2 2 2
n
台灣師範大學機電科技學系 -79-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -80-
C. R. Yang, NTNU MT
Displacement transmissibility
1 2r
2
X
Td
Y 2 2
1 r 2r
2
台灣師範大學機電科技學系 -81-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -82-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -83-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -84-
C. R. Yang, NTNU MT
X 1 k m1 2 X 2 k F0
X 1k X 2 k m2 2 0
台灣師範大學機電科技學系 -85-
C. R. Yang, NTNU MT
0,
2 m1 m2 k
2
1 2
m1m2
Amplitude of m1 at steady-state:
X1
k m F 2
2
0
k m k m k
1
2
2
2 2
台灣師範大學機電科技學系 -86-
C. R. Yang, NTNU MT
Amplitude of m2 at steady-state:
kF0
X2
k m1 2 k m2 2 k 2
Force transmitted to structure:
m k 2
F0
Ft m 2 X 2
2 2
k m1 2 k m2 2 k 2
台灣師範大學機電科技學系 -87-
C. R. Yang, NTNU MT
Transmissibility of isolator:
Ft m2 k 2
Tr
F0 k m1 2 k m2 2 k 2
1 m2 1
m1 m2 m1 2 m1 m2 2
1 2
m2 k 2
where ω2 is the natural frequency of the system
Ft decrease as ω2 decrease
台灣師範大學機電科技學系 -88-
C. R. Yang, NTNU MT
Example 9.4
Spring Support for Exhaust Fan
台灣師範大學機電科技學系 -89-
C. R. Yang, NTNU MT
Example 9.4
Spring Support for Exhaust Fan
Solution
Transmissibility = 0.1
2
1 2
n
0.1 2
2
2
1 2
n n
1000 2
Forcing frequency 104.72 rad/s
60
台灣師範大學機電科技學系 -90-
C. R. Yang, NTNU MT
Example 9.4
Spring Support for Exhaust Fan
Solution
k 4K K
Natural frequency: n
m 40 3.1623
台灣師範大學機電科技學系 -91-
C. R. Yang, NTNU MT
Example 9.6
Isolator for Stereo Turntable
台灣師範大學機電科技學系 -92-
C. R. Yang, NTNU MT
Example 9.6
Isolator for Stereo Turntable
Solution
台灣師範大學機電科技學系 -93-
C. R. Yang, NTNU MT
Example 9.8
Isolation from Vibrating Base
台灣師範大學機電科技學系 -94-
C. R. Yang, NTNU MT
Example 9.8
Isolation from Vibrating Base
Solution
1 2
2
1 1
T
Setting ω=ωn, r or 0.1291
2 2 Tr 1 2 15
2
台灣師範大學機電科技學系 -95-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -96-
C. R. Yang, NTNU MT
X1
k Z X 2
k Z F0
k Z k m km
1
2
1
2
kF0
X2
Z k m1 2 km1 2
台灣師範大學機電科技學系 -97-
C. R. Yang, NTNU MT
kZ F0
Amplitude of transmitted force: Ft X 2 Z
Z k m1 2 km1 2
Ft kZ
Transmissibility of isolator Tr
F0
Z k m1 2 km1 2
台灣師範大學機電科技學系 -98-
C. R. Yang, NTNU MT
• Shock Isolation
Shock load is a force load applied for less than one natural time
period of the system
T
Impulse: F F t dt
0
F
Velocity imparted to the mass, v
m
i.e. application of shock load is equivalent to giving an initial
velocity to the system
台灣師範大學機電科技學系 -99-
C. R. Yang, NTNU MT
• Shock Isolation
ve nt
Free vibration solution: xt sin d t where d 1 2 n
d
d
cd
1
where tan
k c n
台灣師範大學機電科技學系 -100-
C. R. Yang, NTNU MT
Example 9.11
Isolation Under Shock
台灣師範大學機電科技學系 -101-
C. R. Yang, NTNU MT
Example 9.11
Isolation Under Shock
Solution
台灣師範大學機電科技學系 -102-
C. R. Yang, NTNU MT
Example 9.11
Isolation Under Shock
Solution
台灣師範大學機電科技學系 -103-
C. R. Yang, NTNU MT
Example 9.12
Isolation Under Step Load
台灣師範大學機電科技學系 -104-
C. R. Yang, NTNU MT
Example 9.12
Isolation Under Step Load
台灣師範大學機電科技學系 -105-
C. R. Yang, NTNU MT
Example 9.12
Isolation Under Step Load
Solution
台灣師範大學機電科技學系 -106-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -107-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -108-
C. R. Yang, NTNU MT
9.11
Vibration Absorbers
9.11
台灣師範大學機電科技學系 -109-
C. R. Yang, NTNU MT
• When the excitation freq coincides with the ωn, the system may
experience excessive vibration.
台灣師範大學機電科技學系 -110-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -111-
C. R. Yang, NTNU MT
Amplitude of masses:
X1
k2 m 2 2
F0
k1 k2 m1 2 k2 m2 2 k22
k 2 F0
X2
k1 k2 m1 2 k2 m2 2 k22
We want to reduce X1. Thus set numerator of X1 to zero.
k2 k
2 , 2 12 1
m2 m1
台灣師範大學機電科技學系 -112-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -113-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -114-
C. R. Yang, NTNU MT
At X1=0, ω= ω1, k1 F0
X 2 st
k2 k2
台灣師範大學機電科技學系 -115-
C. R. Yang, NTNU MT
4 2
2
2
m 2
1 1 2 2 1 0
2 1 2 m1 1
2 roots of the equation:
2
1 1 1 m2 2 m
2 2 2 2
2
1 1 2 2 4
2 m1 1 m1 1 1
2
2
2
2 2
2 1
台灣師範大學機電科技學系 -116-
C. R. Yang, NTNU MT
Example 9.15
Vibration Absorber for Diesel Engine
台灣師範大學機電科技學系 -117-
C. R. Yang, NTNU MT
Example 9.15
Vibration Absorber for Diesel Engine
Solution
6000
We have f 100Hz or 628.32 rad/s
60
m2 0.31665 kg
k 2 2 m2 628.32 0.31665 125009 N/m
2
台灣師範大學機電科技學系 -118-
C. R. Yang, NTNU MT
Example 9.16
Absorber for Motor-Generator Set
台灣師範大學機電科技學系 -119-
C. R. Yang, NTNU MT
Example 9.16
Absorber for Motor-Generator Set
台灣師範大學機電科技學系 -120-
C. R. Yang, NTNU MT
Example 9.16
Absorber for Motor-Generator Set
Solution k1 k2 m2
1 , 2 ,
m1 m2 m1
1 261.80 rad/s, 2 366.52 rad/s
1 261.80
r1 0.8333
2 314.16
2 366.52
r2 1.1667
2 314.16
2
r12 , r22 1 1 1
2 2
台灣師範大學機電科技學系 -121-
C. R. Yang, NTNU MT
Example 9.16
Absorber for Motor-Generator Set
Solution
r14 1
2 2 0.1345
r1
m2
m1 14.8699 kg
0.1345
1 209.44 rad/s
1 209.44
r1 0.6667
2 314.16
台灣師範大學機電科技學系 -122-
C. R. Yang, NTNU MT
Example 9.16
Absorber for Motor-Generator Set
Solution
台灣師範大學機電科技學系 -123-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -124-
C. R. Yang, NTNU MT
Equations of motion
x
Assume solution: j t X j e i t
, j 1,2
Steady-state solutions:
X1
F0 k 2 m2 2 ic2
k m
1 1
2
k 2 m2 2 m2 k 2 2 ic2 k1 m1 2 m2 2
X 1 k 2 ic2
X2
k 2 m2 2 ic2
台灣師範大學機電科技學系 -125-
C. R. Yang, NTNU MT
台灣師範大學機電科技學系 -126-
C. R. Yang, NTNU MT
X1
2g 2
2
g f
2 2
X2
2g 2 f 4
st 2g 2 g 2 1 g 2 2 f 2 g 2 g 2 1g 2 f2 2
台灣師範大學機電科技學系 -127-
C. R. Yang, NTNU MT
1 f 2 f 2 2 f 2
g 2 g
4 2
0
2 2
台灣師範大學機電科技學系 -128-
C. R. Yang, NTNU MT
3
2
2
for point A
81
3
台灣師範大學機電科技學系 -129-
C. R. Yang, NTNU MT
3
2
81
optimal 3
X1 X1 2
1
st optimal st max
台灣師範大學機電科技學系 -130-