Lecture 3-ME2617 Forced Vibration of SDOF Systems
Lecture 3-ME2617 Forced Vibration of SDOF Systems
Lecture 3-ME2617 Forced Vibration of SDOF Systems
Dr Cristinel Mares
[email protected]
OBJECTIVES
• Determine the response of a vibrating SDOF system to an
external force
k c 𝜔 = ;𝑐 = 2 𝑘𝑚; 𝜁 =
𝐹(𝑡)
m 𝑥̈ (𝑡) + 2𝜁𝜔 𝑥̇ (𝑡) + 𝜔 𝑥(𝑡) =
𝑚
If we consider the underdamped free oscillation <1, then the response is given by
x ( t ) exp( n t )( A sin d t B cos d t ) exp( n t ) sin( d t )
To obtain the particular solution one must know the excitation F(t).
Let’s consider one type of oscillation only
𝑋 k − mω 𝑐𝑜𝑠𝜑 + 𝑐𝜔𝑠𝑖𝑛𝜑 = 𝐹 ⟹ 𝑋 =𝐹
𝐹
⟹𝑋 =
k − mω + 𝑐𝜔
One obtains
𝐹 𝑐𝜔
𝑋 = 𝑡𝑎𝑛𝜑 =
𝑘 − 𝑚𝜔 + 𝑐𝜔 𝑘 − 𝑚𝜔
𝜔 𝑐 𝜔 𝑘
𝑘 − 𝑚𝜔 + 𝑐𝜔 = 𝑘 1− + 2 𝑘𝑚 =𝑘 1−𝑟 + 2𝜁𝑟
𝜔 𝑐 𝜔 𝑚
𝜔 𝑐 𝐹
Let 𝑟= ; 𝜁= ; 𝛿 =
𝜔 2 𝑘𝑚 𝑘
𝑋 1 2𝜁𝑟 1
= 𝑡𝑎𝑛𝜑 = 𝑀=
𝛿 1−𝑟 + 2𝜁𝑟 1−𝑟 1−𝑟 + 2𝜁𝑟
M magnification factor
1
𝑀=
1−𝑟 + 2𝜁𝑟 M magnification factor
At = 1 the response
increases and the lower level
of damping the bigger the
response level. This
phenomenon is called
resonance.
2𝜁𝑟
𝑡𝑎𝑛𝜑 =
1−𝑟
𝜑 = 90° at resonance =1
𝑘
𝜔 = = 100𝑟𝑎𝑑/𝑠
𝑚
𝑐 = 2 𝑘𝑚 = 200𝑁𝑠/𝑚
𝜁= = 0.1
Brunel University London
Example
04/02/2021
1 2
𝑥 𝑡 = 𝑋 cos(𝜔𝑡 − 𝜑 = cos(50𝑡 − 𝐴𝑟𝑐𝑡𝑎𝑛 )
5 229 15
𝑥 𝑡 = 𝑥 (𝑡) + 𝑥 (𝑡
𝑥 𝑡 = 𝑥 (𝑡) + 𝑥 (𝑡
𝑥 (𝑡)
𝑥 (𝑡)
It is of the form
𝜔
Beyond = 2 the non-dimensional response is smaller than 1, the amplitude is smaller than
𝜔
the static deflection.
𝑥 𝑡 = 𝑥 𝑠𝑖𝑛 𝜔𝑡 − 𝜑
𝜋
𝑥̇ 𝑡 = 𝑥 𝜔𝑐𝑜𝑠 𝜔𝑡 − 𝜑 = 𝑥 𝜔𝑠𝑖𝑛 𝜔𝑡 − 𝜑 +
2
F0 c
x0 arctan 2
(k m2 )2 (c)2 k m
Let and
x0 1 2 r
arctan 2
F0 k (1 r 2 )2 (2 r)2 1 r
Stiffness control
Damping control
Inertia control
𝐹 𝛿 𝐹
𝑥 𝑡 = 𝑋 𝑐𝑜𝑠𝜔𝑡 = 𝑐𝑜𝑠𝜔𝑡 = ;𝛿 =
𝑘 − 𝑚𝜔 𝜔 𝑘
1− 𝜔
k
𝐹 𝑥̇ 𝑘
𝑥 𝑡 = 𝑥 − 𝑐𝑜𝑠𝜔 𝑡 + 𝑠𝑖𝑛𝜔 𝑡; 𝜔 =
𝑘 − 𝑚𝜔 𝜔 𝑚
m
x
𝐹 𝑥̇ 𝐹
𝑥 𝑡 = 𝑥 − 𝑐𝑜𝑠𝜔 𝑡 + 𝑠𝑖𝑛𝜔 𝑡 + 𝑐𝑜𝑠𝜔𝑡
𝑘 − 𝑚𝜔 𝜔 𝑘 − 𝑚𝜔
𝐹 𝑡 = 𝐹 𝑐𝑜𝑠𝜔𝑡
𝑋 1
=
𝛿 𝜔
1− 𝜔
𝑋 1
=
𝛿 𝜔
1−
𝜔
𝑋 1
=−
𝛿 𝜔
−1
𝜔
When 𝜔 = 𝜔 𝑟=1
the denominator 1 −
is zero and amplitude becomes
infinite, the system is at resonance.
𝑥̇ 𝛿 𝜔 𝑡
𝑥 𝑡 = 𝑥 𝑐𝑜𝑠𝜔 𝑡 + 𝑠𝑖𝑛𝜔 𝑡 + 𝑠𝑖𝑛𝜔 𝑡
𝜔 2
When 𝜔 ≅ 𝜔 𝑟 ≅ 1
the beating phenomenon occurs.
The equation of motion, beating
circular frequency and period are:
𝐹 1
𝑥 𝑡 =( 𝑠𝑖𝑛𝜀𝑡)𝑠𝑖𝑛𝜔𝑡
𝑚 2𝜀𝜔
2𝜋 2𝜋
𝜔 = 2𝜀 = 𝜔 − 𝜔; 𝑇 = =
2𝜀 𝜔 − 𝜔
𝑐 = 0; 𝜔 = 5𝑟𝑎𝑑/𝑠; 𝐹 = 100𝑁
The characteristics of the vibrating systems are: 𝐹 𝑡 = 𝐹 𝑐𝑜𝑠𝜔𝑡
𝑘
𝜔 = = 100𝑟𝑎𝑑/𝑠
𝑚
̇
𝑥 𝑡 = 𝐴 𝑐𝑜𝑠𝜔 𝑡 + 𝐵 𝑠𝑖𝑛𝜔 𝑡= 𝑥 − 𝑐𝑜𝑠𝜔 𝑡 + 𝑠𝑖𝑛𝜔 𝑡
=−0.000025 𝑐𝑜𝑠100𝑡 − 𝑠𝑖𝑛100𝑡
The particular solution
𝐹 4
𝑥 𝑡 = cos 𝜔𝑡 = cos 5𝑡
𝑘 − 𝑚𝜔 399
𝑥 𝑡 = 𝑥 (𝑡) + 𝑥 (𝑡
𝑥 𝑡 = 𝑥 (𝑡) + 𝑥 (𝑡
𝑥 (𝑡)
𝑥 (𝑡)
The system under base excitation 𝑦 𝑡 = 𝑌𝑠𝑖𝑛𝜔𝑡 The system is loaded by the harmonic force
What is the ratio between displacements X 𝐹 𝑡 = 𝐹 𝑠𝑖𝑛𝜔𝑡
of mass m and Y at the source (base) ? What is the ratio between the force 𝐹
transmitted to the base and 𝐹 at the source of
excitation ?
𝐹
k c
𝐹 𝑠𝑖𝑛𝜔𝑡
𝑚𝑥̈ + 𝑐 𝑥̇ − 𝑦̇ + 𝑘 𝑥 − 𝑦 = 0
𝑦 = 𝑌𝑠𝑖𝑛𝜔𝑡; 𝑦̇ = 𝑌𝜔𝑐𝑜𝑠𝜔𝑡;
𝑚𝑥̈ + 𝑐𝑥̇ + 𝑘𝑥 = 𝑐𝑦̇ + 𝑘𝑦
𝑚𝑥̈ + 𝑐𝑥̇ + 𝑘𝑥 = 𝑐𝑜𝑠𝜔𝑡 + 𝒌𝒀𝑠𝑖𝑛𝜔𝑡
Observe that sin(𝜔𝑡 + 𝜙 ) = sin𝜔𝑡𝑐𝑜𝑠𝜙 + sin𝜙 cos𝜔𝑡 and the phase 𝜙 must be determined
𝐴 𝐵
𝐴𝑠𝑖𝑛𝜔𝑡 + 𝐵𝑐𝑜𝑠𝜔𝑡 = 𝐴 +𝐵 𝑠𝑖𝑛𝜔𝑡 + 𝑐𝑜𝑠𝜔𝑡 = 𝐶𝑠𝑖𝑛 𝜔𝑡 + 𝜙
𝐴 +𝐵 𝐴 +𝐵
35 © 2011 Mechanical Vibrations Fifth Edition in SI Units
Brunel University London
RESPONSE UNDER HARMONIC MOTION OF THE BASE
𝑐𝜔
𝐹 = 𝑌 𝑘 + 𝑐𝜔 and 𝑡𝑎𝑛𝜙 = −
𝑘
it will respond harmonically at the same frequency and there will be a phase lag between
the force and the response:
𝑥 𝑡 = 𝑋 cos(𝜔𝑡 − 𝜙 )
36
Brunel University London
RESPONSE UNDER HARMONIC MOTION OF THE BASE
𝐹 𝑐𝜔
𝑋 = 𝑡𝑎𝑛𝜙 =
𝑘 − 𝑚𝜔 + 𝑐𝜔 𝑘 − 𝑚𝜔
37
Brunel University London
RESPONSE UNDER HARMONIC MOTION OF THE BASE
𝑥 𝑡 = 𝑋 sin[(𝜔𝑡 − 𝜙 ) − 𝜙 ]
𝑋 = =
𝑥 𝑡 = 𝑠𝑖𝑛 [𝜔𝑡 − 𝜙 + 𝜙 ]
𝑡𝑎𝑛𝜙 = − ; 𝜙 = 𝑡𝑎𝑛
with
( )
40
Brunel University London
RESPONSE UNDER HARMONIC MOTION OF THE BASE
𝐹 𝑚𝜔 𝑌 𝑘 + 𝑐𝜔 1 + 2𝜁𝑟
= =𝑟
𝑘𝑌 𝑘𝑌 𝑘 − 𝑚𝜔 + 𝑐𝜔 1−𝑟 + 2𝜁𝑟
𝑚𝑥̈ + 𝑐 𝑥̇ − 𝑦̇ + 𝑘 𝑥 − 𝑦 = 0
𝑦 𝑡 = 𝑌𝑠𝑖𝑛𝜔𝑡 ⇒ 𝑚𝑥̈ + 𝑐𝑥̇ + 𝑘𝑥 = 𝑘𝑌𝑠𝑖𝑛𝜔𝑡 + 𝑐𝜔𝑌𝑐𝑜𝑠𝜔𝑡
m𝜔 𝑋
𝑐𝜔𝑋
𝑐𝜔𝑌
𝑘𝑌
𝜑
𝑐𝜔𝑌 + 𝑘𝑌 𝑘𝑋
𝜔𝑡 − 𝜑
𝐹 = 𝑐𝑥̇ + 𝑘𝑥 𝐹 𝐹
k c
k c
m
𝑥
𝐹𝑠𝑖𝑛𝜔𝑡 𝐹
𝐹
𝐹
mx cx kx Fsint
𝑥 𝑡 = 𝐴𝑠𝑖𝑛 𝜔𝑡 − 𝜙 m
𝐹 𝑠𝑖𝑛𝜔𝑡
43 © 2011 Mechanical Vibrations Fifth Edition in SI Units
Brunel University London
TRANSMISSIBILITY
The forces acting on the mass are the spring force proportional to displacement, the damping force
proportional to velocity, the inertia force proportional to acceleration.
The system is connected to the foundation through the spring and damper only. Therefore, the force
transmitted is the sum of the reactions at the fixed end of the spring and dashpot and this can be seen better
on a vector (phasor) diagram:
m𝜔 𝑋
𝑐𝜔𝑋
𝐹
𝐹 = 𝑐𝑥̇ + 𝑘𝑥
𝜑
𝐹 𝑘𝑋
𝜔𝑡 − 𝜑
𝐹 = 𝑐𝑥̇ + 𝑘𝑥
𝐹
𝑥 𝑡 = 𝐴𝑠𝑖𝑛 𝜔𝑡 − 𝜙
Force transmissibility:
k c
𝐹 𝑘 + 𝑐𝜔 1 + 2𝜁𝑟
= =
𝐹 𝑘 − 𝑚𝜔 + 𝑐𝜔 1−𝑟 + 2𝜁𝑟
m
𝐹 𝑠𝑖𝑛𝜔𝑡
𝑋𝑀 𝑟
= The curves begin at zero amplitude the damping level
𝑚𝑒 1−𝑟 + 2𝜁𝑟
affecting the response amplitude.
for
Forced Vibrations
Frequency domain solution:
𝐹 𝛿
𝑋= = ;
𝑘 − 𝑚𝜔 + 𝑐𝜔 1−𝑟 + 2𝜁𝑟
𝜔 𝐹
𝑟= ; 𝛿 =
𝜔 𝑘
𝑐𝜔
𝜙 = 𝑡𝑎𝑛
𝑘 − 𝑚𝜔
Frequency Response Function:
1
𝐻 𝜔 =
𝑘 − 𝑚𝜔 + 𝑐𝜔