Inno2018 Emt2502 Vibrations Undamped Mdof
Inno2018 Emt2502 Vibrations Undamped Mdof
Inno2018 Emt2502 Vibrations Undamped Mdof
By
Inno Odira
c Copyright by Inno Odira, 2018
Table of Contents
Bibliography 6
REFERENCES 6
iii
1
0.1 MDOF Undamped Forced Response
Example:
Consider a simple 2DOF mass- spring system Fig. 1
For m1
For m2
In Matrix form:
" #" # " #" # " #
m1 0 ẍ1 k1 + k2 −k2 x1 f1
+ = (1)
0 m2 ẍ2 −k2 k2 + k3 X2 f2
| {z } | {z } | {z } | {z } | {z }
M Ẍ K X F
In a familiar form
M Ẍ + KX = F (2)
Assume The forcing function is harmonic:
Equivalent to
[K − ω 2 M ]X̄ = F̄ (5)
To have a solution [K − ω 2 M ] must be invertible.
Thus
[K − ω 2 M ]X̄ = F̄ → X̄ = [K − ω 2 M ]−1 F (6)
For our example problem.
X̄ = [K − ω 2 M ]−1 F
" #−1
ki + k2 − ω 2 m1 −k2
= F̄
−k2 k2 + k3 − ω 2 m1
" #" #
1 k2 + k3 − ω 2 m1 k2 f¯1
=
det[K − ω 2 M ] k2 ki + k2 − ω 2 m1 f¯2
Let
∆(ω) = det [K − ω 2 M ] = (k1 + k2 − ω 2 m1 )(k2 + k3 − ω 2 m2 ) − k22
| {z }
characteristic polynomial
2. If numerator num 6= 0 at natural frequencies, then the response goes ∞ at each natural
frequency as shown in Fig . 2
3
For m2
m2 ẍ2 = k1 (x1 − x2 ) ⇒ m2 ẍ2 + k1 x2 − k1 x1 = 0
In Matrix form:
" #" # " #" # " #
m1 0 ẍ1 k1 + k2 −k2 x1 f
+ = (7)
0 m2 ẍ2 −k2 k1 X2 0
" #
f¯
f (t) = f¯eiωt so F = eiωt F = F̄ eiωt
0
We can assume a solution x(t) to match the form of f (t):
X̄ = [K − ω 2 M ]−1 F
" #−1 " #
ki + k2 − ω 2 m1 −k2 f¯
X̄ =
−k2 k2 + −ω 2 m2 0
" #" #
1 k2 − ω 2 m2 k2 f¯
=
det[K − ω 2 M ] k2 ki + k2 − ω 2 m1 0
Let
∆(ω) = det [K − ω 2 M ] = (k1 + k2 − ω 2 m1 )(k2 − ω 2 m2 ) − k22
| {z }
characteristic polynomial
Look at the FBD 3, If the Spring force Fsp2 is exactly opposite of f (t) for all time, then
there are no external forces on m1 (Assume x1 = 0)
Let us look at the equation of motion
q q
Assume f (t) = f cos(ωn t) where ωn = m1 = mk22
k1
With
k1 k2
∆(ω) = (k1 +k2 −ω 2 m1 )(k2 −ω 2 m2 )−k22 = (k1 +k2 − m1 )(k2 − m2 )−k22 = (k2 )(0)−k22
m1 m2
Thus
(k2 )f¯
X̄2 =
−k22
Then
−f¯
X̄2 (t) = cos ωn t
k2
And the resulting spring force
−f¯
Fsp2 = k2 (X2 ) = k2 cos ωn t = −f¯ cos(ωn t)
k2 | {z }
Exactly cancels the excitation force
REFERENCES