Mechanical Vibration Solved Examples

Mechanical Vibrations: 4600-431 Example Problems
December 20, 2006
Contents
1 Free Vibration of Single Degree-of-freedom Systems 1
2 Frictionally Damped Systems 33
3 Forced Single Degree-of-freedom Systems 42
4 Multi Degree-of-freedom Systems 69
1 Free Vibration of Single Degree-of-freedom Systems
Problem 1:
In the figure, the disk and the block have
mass m and the radius of the disk is r.
a) Find the equations of motion for this

system.
c (m, r)
b) What are the natural frequency and
k m
damping ratio of the system in terms
of m, c, and k?
c) If the block is displaced 18 cm to the

right and released from rest find the re-
sulting angular displacement of the disk
with
k 4k
m = 3 kg, r = 9 cm,
k = 21 N/m, c = 63 N · s/m,
1
Problem 2: x
For the system shown to the right, the disk of
mass m rolls without slip and x measures the
displacement of the disk from the unstretched
k
position of the spring.
a) Find the equations of motion.

c
(m, r)
b) With
c = 16 N/(m/s), m = 2 kg,
r = 0.10 m
for what value of the spring stiffness k

is the damping ratio of the system one-
half of the critically damped value, so
that ζ = ζcr /2?
c) With these parameter values, find the
displacement of the disk if it is rolled
20 cm to the right (from static equilib-
rium) and released from rest.
Problem 3:
From the figure shown to the right
k2
a) find the equations of motion in terms of
the angular rotation of the disk;
r k1
b) what are the damping ratio and natural
frequency of the system in terms of the
parameters m, b, k1 , and k2 ; r
m 2
c) can you draw an equivalent spring-
mass-damper system?
2
Problem 4:
For the mechanical system shown to the
right, the uniform rigid bar has mass m and ̂
pinned at point O. For this system:
k
a) find the equations of motion;
ı̂
b) Identify the damping ratio and natural
frequency in terms of the parameters m
m, c, k, and ℓ.
ℓ
2 y
c) For: ℓ
θ 2
m = 1.50 kg, ℓ = 45 cm,

O
c = 0.125 N/(m/s), k = 250 N/m,
2m
find the angular displacement of the bar x
θ(t) for the following initial conditions:
k c
θ(0) = 0, θ̇(0) = 10 rad/s.
Assume that in the horizontal position the

system is in static equilibrium and that all
angles remain small.
Solution:
a) In addition to the coordinate θ identified in
the original figure, we also define x and y as −k y ̂
the displacment of the block and end of the
bar respecively. The directions ı̂ and ̂ are
defined as shown in the figure.
A free body diagram for this system is
shown to the right. Note that the tension in −T ̂
the cable between the bar and the block is
unknown and represented with T while the FR T ̂
reaction force FR is included, although both
its magnitude and direction are unspecified.
In terms of the identifed coordinates, the
angular acceleration of the bar αβ/F and the
F
linear acceleration of the block aG are
−k x ̂ −c ẋ ̂
F
αβ/F = θ̈ k̂, aG = ẍ ̂.
We can also relate the identified coordinates as

ℓ
x = θ, y = ℓ θ.
2
The equations of motion for this system can be obtained with linear momentum balance
applied to the block and angular momentum balance aout O on the bar. These can be
written as
X ³ ´
F
F = m aG −→ T − k x − c ẋ ̂ = 2 m ẍ ̂,
m ℓ2
µ ¶
X ℓ
M O = I O αβ/F −→ − T − k y ℓ k̂ = θ̈ k̂.
2 3
3
Solving the first equation for T and substituting into the second equation yields
ℓ m ℓ2
− (2 m ẍ + k x + c ẋ) −kyℓ = θ̈.
2 3
Using the coordinate relations we can obtain the equation of motion as
5 m ℓ2 c ℓ2 5 k ℓ2
θ̈ + θ̇ + θ = 0.
6 4 4
b) In the above equation the equivalent mass, damping, and stiffness are
5 m ℓ2 c ℓ2 5 k ℓ2
meq = , beq = , keq = .
6 4 4
From these the damping ratio and natural frequency are

c ℓ2
√
beq 3c
ζ= p = q 4 = √ ,
2 keq meq 2 5 k ℓ 2 5 m ℓ2 2 50 km
4 6
s s 2
5kℓ
r
keq 4 3k
ωn = = 5 m ℓ2
=
meq 6
2 m
c) Evaluating the damping ratio and natural frequency we find that for the given values
of the parameters
ωn = 15.8 rad/s, ζ = 7.91 × 10−4 .
Therefore the system is underdamped and the general solution can be written as
³ ³ p ´ ³ p ´´
θ(t) = e−ζ ωn t a sin ωn 1 − ζ 2 t + b cos ωn 1 − ζ 2 t ,
where a and b are arbitrary constants used to fit the initial conditions. Evaluating θ(t)
and θ̇(t) at t = 0 yields
p
θ(0) = b = 0, θ̇(0) = −ζ ωn b + ωn 1 − ζ 2 a = 10 rad/s,
so that the general solution becomes

³ ´
θ(t) = 0.632 e−0.0125 t sin 15.8 t .
4
Problem 5:
In the figure, the disk has mass m, radius r,
2
and moment of inertia IG = mr 2 about the
mass center G and is assumed to roll without k g
slip. The identified coordiante θ measures the
rotation of the disk with respect to the equi-
librium position. c
G
system. k (m, r)
b) If the disk is released from rest with
θ(0) = − π2 rad, find the resulting an-
gular displacement θ(t) for
m = 3 kg, r = 15 cm,
k = 36 N/m, c = 3 N · s/m,
c) What is the force in the upper spring

during this motion?
Problem 6:
right, the uniform rigid bar has mass m and
pinned at point O. For this system: ℓ
k 2
m θ
b) Identify the damping ratio and natural z
frequency in terms of the parameters
ℓ
m, c, k, and ℓ. 2 m
c) For: ̂
x
k c
m = 2 kg, ℓ = 25 cm,
c = 0.25 N/(m/s), k = 50 N/m,
ı̂
find the angular displacement of the bar

θ(t) for the following initial conditions:
θ(0) = 0, θ̇(0) = 10 rad/s.
d) for this motion, find the tension in the

cable connecting the rod and the block
as a function of time.
Assume that the system is in static equilib-

rium at θ = 0, and that all angles remain
small.
Solution:
a) We identify the coordinates x and z as shown above, which are related to the angular
5
displacement θ as:
ℓ ℓ
x= θ, z= θ.
2 2
An appropriate free-body diagram is

k z ̂
shown to the right. Applying linear mo-
mentum balance on the block yields
−T ̂
X F
F = m aG ,
T ̂
(T − k x − c ẋ) ̂ = m ẍ ̂. FR
Likewise, angular momentum balance on
the bar provides
X
MO = I O αβ/F ,
−k x ̂ −c ẋ ̂
m ℓ2
µ ¶
ℓ ℓ
−T − k z k̂ = θ̈ k̂.
2 2 12
Combining these equations and eliminating the tension, the equation of motion can be
written as
7m
θ̈ + c θ̇ + 2 k θ = 0.
6
b) For the above equation the equivalent mass, damping, and stiffness are
7m
meq = , beq = b, keq = 2 k,
6
and the natural frequency and damping ratio are

s r √
keq 12 k beq 3b
ωn = = , ζ= p =√ .
meq 7m 2 keq meq 28k m
Problem 7: ̂
The block shown to the right rests on a fric-
tionless surface. Find the response of the sys-
tem if the block is displaced from its static x
equilibrium position 15 cm to the right and k1 ı̂
released from rest. k2
m
m = 4.0 kg, b = 0.25 N/(m/s),
b
k1 = 1.5 N/m, k2 = 0.50 N/(m/s).
Solution:
An appropriate free-body diagram is shown to the
right. Notice that the two springs are effectively in −k1 x ı̂
parallel, as the displacement across each spring is
−k2 x ı̂
identical. Linear momentum balance on this block
provides
X F −b ẋ ı̂
F = m aG ,
(−k1 x − k2 x − b ẋ) ı̂ = m ẍ ı̂,
6
or, writing this in standard form
m ẍ + b ẋ + (k1 + k2 ) x = 0.
Further, the system is released from rest so that the initial conditions are
x(0) = x0 = 15 cm, ẋ(0) = 0 cm/s.
Problem 8:
displacement of the disk from the unstretched k x
position of the spring. The surface is inclined
z
at an angle of φ with respect to vertical.
a) find the equations of motion. Do not

neglect gravity; c
b) if the system is underdamped, what is φ (m, r) θ

the frequency of the free vibrations of
this system in terms of the parameters
k, c, and m; ̂ ê2 C
c) for what value of the damping constant ı̂

c is the system critically damped;
d) what is the static equilibrium displace- ê1

ment of the disk?
Solution:
a) In addition to x, the displacement of the center of the disk, we identify the coordinates
z and θ, the displacement across the spring and the rotation of the disk respectively.
These additional coordinates are related to x as

−k z ê1
z = 2 x, x = −r θ.
−c ẋ ê1 −m g ̂
An appropriate free-body diagram is shown to the
right. We note that (ı̂, ̂) are related to the direc-
tions (ê1 , ê2 ) as
ı̂ = cos φ ê1 + sin φ ê2 ,

̂ = − sin φ ê1 + cos φ ê2 . fr ê1
N ê2
The moment produced by gravity about point C is

Mgravity = rGC × (−m g ̂),
= (r ê2 ) × (−m g ̂) = −m g r sin φ k̂.
Angular momentum balance about the contact point C yields

X
MC = I C αD/F ,
3 m r2
³ ´ µ ¶
(2 r) k z + r c ẋ − m g r sin φ k̂ = θ̈ k̂.
2
7
Eliminating the coordinates z and θ, we can write the equation of motion in terms of x
as
3m
ẍ + c ẋ + 4 k x = m g sin φ.
2
Since the gravitational force has been included in the development of this equation of
motion, the coordinates are measured with respect to the unstretched position of the
spring.
b) Assuming
p the system is underdamped, the frequency of the free vibrations is ωd =
ωn 1 − ζ 2 , where
s r
keq 8k beq c
ωn = = , ζ= p = √ ,
meq 3m 2 keq meq 2 6km
so that r r
8k c2
ωd = 1− ,
3m 24 k m
c) The system is critically damped when ζ = 1, which corresponds to a damping coefficient

of √
ccr = 2 6 k m.
d) The system is stationary in static equilibrium, so that x ≡ x0 = constant—both ẋ and

ẍ vanish, and the equation of motion reduces to
4 k x0 = m g sin φ.
Solving for x0 , the equilibrium displacement is

m g sin φ
x0 = .
4k
8
Problem 9:
In the figure shown to the right, in the ab- z2
sence of gravity the springs are unstretched k2
in the equilibrium position.
a) Determine the deflection of each spring r2

from its unstretched length when the
system shown is in equilibrium. θ
̂ r1
b) If the system is released from the un- I
stretched position of the springs, what
is the maximum angular velocity of the
disk during the resulting motion? ı̂
m z1
k1
Solution:
a) We define the coordinates x, θ, z1 , and z2
−k2 z2 ı̂
as shown in the figure, which are related
as
x = −r1 θ, z1 = r1 θ, z2 = −r2 θ.
Notice that because of these coordinate

definitions, a rotation with positive θ
gives rise to a negative value in both x
FR
and z2 . Likewise, we see that x = −z1 .
−T ̂ k1 z1 ̂
Using the free-body diagram shown to
the right, linear momentum balance on
T ̂
the block provides
X F
F = m aG ,
(T − m g) ̂ = m ẍ ̂,
−m g ̂
while angular momentum balance on the
disk yields
X
MO = I O αD/F ,
(T r1 − k1 z1 r1 + k2 z2 r2 ) k̂ = I θ̈ k̂.
Eliminating the unknown tension T from these equations and using the coordinate
relations, the equation of motion becomes
³ ´ ³ ´
I + m r12 θ̈ + k1 r12 + k2 r22 θ = m g r1 .
9
The equilibrium rotation of the disk thus is found to be
m g r1
θeq = .
k1 r12 + k2 r22
With this, the equilibrium deflection of each spring is found to be
m g r12
z1,eq = r1 θeq = ,
k1 r12 + k2 r22
m g r1 r2
z2,eq = −r2 θeq =− .
k1 r12 + k2 r22
b) The general free response of the disk can be expressed as
θ(t) = θeq + A sin(ωn t) + B cos(ωn t),
where θeq is given above, A and B are arbitrary constants, and

s
k1 r12 + k2 r22
ωn = .
I + m r12
The system is released with the initial conditions:
θ(0) = 0, θ̇(0) = 0,
so that solving for the arbitrary constants
A = 0, B = −θeq .
Therefore the solution is

Ã Ãs !!
m g r1 k1 r12 + k2 r22
θ(t) = θeq (1 − cos(ωn t)) = 1 − cos t .
k1 r12 + k2 r22 I + m r12
The angular velocity of the disk becomes

³ ´
θ̇(t) = θeq ωn sin(ωn t),
which has amplitude

m g r1
Ω = θeq ωn = p
(k1 r12 + k2 r22 ) (I + m r12 )
10
Problem 10: (Spring 2003)
For the spring-mass-damper system shown to
the right, x is measured from the static equi-
librium position. If the surface is assumed to x
be frictionless:
a) determine the governing equations of m

motion; 6k b
b) what is the frequency of oscillation of
the system;
c) what value of the damping coefficient b

corresponds to critical damping?
d) if k = 2 N/m, b = 4 N/(m/s), and m =

3 kg, find the displacement of the mass
x(t) when the system is started with the
initial conditions:
x(0) = 0.10 m, ẋ(0) = 0 m/s.

right, the uniform rigid bar has mass m and k

z2 m
b) what value of the damping constant c
gives rise to a critically damped sys- ℓ
tem? k 2
m θ
z1
ℓ
2 ̂
c
ı̂
Solution:
a) In addition to θ, we define two additional coordinates, z1 and z2 , to measure the deflec-
tion at the left and right ends of the bar.
11
These coordinates are related as:
−k z2 ̂
ℓ ℓ
z1 = θ, z 2 = θ = z1 .
2 2
A free-body diagram for this system is
shown to the right. Applying angular mo- −T ̂
mentum balance on the bar eliminates the
appearance of the reaction force and leads T ̂
k z1 ̂
to: G
X
MG = I G θ̈ k̂,
FR
³ ´ℓ m ℓ2
T − k z1 − cż1 k̂ = θ̈ k̂. c ż1 ̂
2 12
Likewise, application of linear momentum balance on the block yields:

X F
F = m aG ,
³ ´
− k z2 − T ̂ = mz̈2 ̂
Eliminating the unknown tension T and solving for z1 in terms of z2 , the equation of
motion becomes:
m ℓ2 c ℓ2 k ℓ2
θ̈ + θ̇ + θ = 0.
3 4 2
b) A critically damped system occurs when ζ = 1. For this system:

√
3c
ζ=√ .
32 k m
Solving for ccr yields: r

32 k m
ccr = .
3

Find the response of the system shown to the
right if the block is pulled down by 15 cm and ̂
k1
released form rest.
b ı̂
m = 2.0 kg, b = 0.5 N/(m/s),
k1 = 0.5 N/m, k2 = 0.25 N/(m/s). k2
x
Solution:
For this system, the two springs in series may be replaced by an equivalent spring, with
constant:
1 k1 k2
keq = 1 1 = k +k .
k1 + k2 1 2
12
Therefore, the free-body diagram is shown
to the right. Applying linear momentum keq x ̂ b ẋ ̂
balance to the block yields:
X F
F = m aG ,
³ ´
keq x + b ẋ ̂ = −m ẍ ̂,
which can finally be written as:
m ẍ + b ẋ + keq x = 0.
With the numerical values given above, this becomes:

µ ¶ µ ¶
1 N 1 N 3
(2 kg) ẍ + ẋ + x = 0, x(0) = m, ẋ(0) = 0 m/s.
2 m/s 6 m 20
With this, the damping ratio and natural frequency are:

r √
1 −1 3
ωn = s , ζ= .
12 4
Therefore, the system is underdamped and the general response can be written as:
³ ´
x(t) = e−ζ ωn t A cos(ωd t) + B sin(ωd t) .
Using the initial conditions to solve for A and B, we find:

Ã Ã√ ! r Ã√ !!
3 −t/8 13 3 13
x(t) = e cos √ t + cos √ t .
20 8 3 13 8 3

In the figure shown to the right, in the ab- x2
sence of gravity the springs are unstretched
in the equilibrium position. k2
a) Determine the deflection of each spring

from its unstretched length when the r2
system shown is in equilibrium.
θ
b) If the system is released from the un- r1
stretched position of the springs, what
is the maximum angular velocity of the
disk during the resulting motion?
̂
ı̂ m
x1
k1
13
Solution:
a) We define θ, x1 , and x2 as indicated in the above figure. In particular, x1 and x2
represent the displacement in their respective springs as measured from their unstretched
position. These coordinates are related through the following transformations:
x1 = −r1 θ, x2 = −r2 θ.
An appropriate free-body diagram for this −k2 x2 ı̂

system is shown to the right. Notice that
the gravitational force must be included to
determine the equilibrium deflection in the
system. To eliminate the reaction force on G
the disk, angular momentum balance is ap-
plied about the center, yielding:
FR
X
MG = I G θ̈ k̂, −T ̂
³ ´
T r1 + k2 r2 x2 k̂ = I G θ̈ k̂.
T ̂
Also, applying linear momentum balance to
the block yields: −m g ̂
X F −k1 x1 ̂
F = m aG ,
³ ´
T − k1 x1 − m g ̂ = m ẍ1 ̂.
Finally, eliminating the unknown tension from these equations and using the above
coordinate transformations, this single-degree-of-freedom system can be modeled with
the equation: ³ ´ ³ ´
I G + m r12 θ̈ + k1 r12 + k2 r22 θ = (m g r1 ).
This equation of motion determines the equilibrium position θeq (with θ̈eq = 0) to be:
m g r1
θeq = .
k1 r12 + k2 r22
Therefore, the equilibrium displacements in each spring are:
m g r12 m g r1 r2
x1,eq = −r1 θeq = , x2,eq = −r2 θeq = .
k1 r12 + k2 r22 k1 r12 + k2 r22
b) Define new coordinates z1 and z2 , which measure the displacement in springs 1 and 2
with respect to the static equilibrium position, that is:
z1 = x1 − x1,eq , z2 = x2 − x2,eq .
Likewise, let φ represent the angular displacement of the disk from the static equilibrium
position:
φ = θ − θeq .
14
Therefore, the potential energy of this system can be written as:
1 1
V = k1 z12 + k2 z22 ,
2 2
1³ ´
= k1 r12 + k2 r22 φ2 .
2
Also, the kinetic energy becomes:

1 G 2 1
T = I φ̇ + m ż12 ,
2 2
1³ G ´
= I + m r12 φ̇2 .
2
If the system is released from rest at the unstretched position of the springs, then:
m g r1
φ(0) = −θeq = − , φ̇(0) = 0.
k1 r12 + k2 r22
At this initial state, the potential and kinetic energies become:
(m g r1 )2
T0 = 0, V0 = .
2(k1 r12 + k2 r22 )
Because this system is conservative, the total energy, E = T + V remains constant.

Therefore, when the kinetic energy is maximal, the potential energy is minimal, that is:
1³ G ´
V1 = 0, T1 = I + m r12 φ̇2max .
2
Finally, conservation of energy implies that V0 = T1 , and solving for φ̇max we find that:
s
(m g r1 )2
φ̇max = .
(I + m r12 )(k1 r12 + k2 r22 )
G

For the system shown to the right, the disk of x ̂
3k z
displacement of the disk from the unstretched
ı̂
a) find the equations of motion; θ

b) if the system is underdamped, what is k m c
the frequency of the free vibrations of
this system in terms of the parameters
k, c, and m;
Solution:
a) We define the three coordinates as shown as the figure, related as:
x = −r θ, z = −2 r θ, z = 2 x.
15
A free-body diagram for this system is −3k z ı̂
shown to the right. Notice that the force in
the upper spring depends on z, rather than G
x, while the friction force has an unknown
magnitude f . Because the disk is assumed −k x ı̂ C −c ẋ ı̂
to roll without slip, we are unable to specify
the value of f , but instead can relate the dis- f ı̂
placement and rotation of the disk through
the coordinate relations above.
The equations of motions can be developed directly with angular momentum balance
about the contact point, so that:
X
MC = I C θ̈ k̂,
³ ´ 3 m r2
(3k z) 2r + (k x) r + (c ẋ) r k̂ = θ̈ k̂.
2
Finally, writing this equation in terms of a single coordinate, we obtain:
3 m r2
µ ¶
θ̈ + (c r2 ) θ̇ + 13 k r2 θ = 0.
¡ ¢
2
p
b) For an underdamped response, the frequency of oscillation is ωd = ωn 1 − ζ 2 . With
this system, we find that:
r
26 k c
ωn = , ζ=√ ,
3m 78 k m
so that: r
26 k 2 c2
ωd = − .
3m 9 m2
Problem 15:
In the system shown to the right, the pulley r
has mass m and radius r, so that the moment
2
of inertia about the mass center is IG = mr
2 .
a) Find the governing equations of mo-

r
tion; 2
b) Find the frequency of oscillation for free

vibrations of the system;
c) For what value of the damping constant m

is the system critically damped?
m
k
k c
Solution:
16
a) We choose coordinates (x1 , x2 , θ), where x1 measures the displacement of the first mass
in the −̂ direction, x2 measures the displacement of the second mass in the ̂ direction,
and θ measures the angular rotation of the wheel in the k̂ direction. The kinetic and
potential energies for this system are:
1 1 1 mr2 2
T = m ẋ21 + m ẋ22 + θ̇ ,
2 2 2 2
1 1
V = k x21 + k x22 .
2 2
However, these three coordinates are dependent through the transformations:

r
x1 = θ, x2 = r θ.
2
Therefore, the Lagrangian reduces to:

1 7mr2 1 5kr2
µ ¶ µ ¶
L=T −V = θ̇2 − θ2 .
2 4 2 4
Further, the generalized force resulting from the viscous damper becomes:
cr2
Qθ = − θ̇,
4
so that the equation of motion for this system can be reduced to:
(7m) θ̈ + c θ̇ + (5k) θ = 0.
b) For this system the natural frequency and damping ratio are:
r
5k c
ωn = , ζ= √ .
7m 2 35 km
Therefore, the damped natural frequency becomes:

r r
p
2
5k c2
ωd = ωn 1 − ζ = 1− .
7m 140 km
c) For a critically damped system, ζ = 1, so that we may solve for c = ccr to yield:
√
ccr = 2 35 km.
Problem 16:
2
mass center G, and is attached to a block of b
mass m which rolls across the surface. If the G
disk rolls without slip (µ is sufficiently large), k m
while the block moves without friction:
a) find the equations of motion for this

system;
b) for m = 2 kg, b = 0.5 (N · s)/m, and

k = 8 N/m, find the frequency of oscil-
lation for the system.
17
Solution:
a) Let x denote the translational displacement of G in the ı̂ direction, while θ denotes the
angular displacement of the disk in the k̂ direction. If T is the tension between the disk
and the mass and f is the frictional force acting on the disk, the equations of motion
on the disk and the mass are:
(−kx − bẋ + T + f ) ı̂ = mẍ ı̂,

mr2
(rf ) k̂ = θ̈ k̂,
2
(−T ) ı̂ = mẍ ı̂.
In addition, x and θ are related through the kinematic constraint:
x = −rθ.
Eliminating (θ, T, f ) we obtain a single degree-of-freedom system on x of the form:

5m
ẍ + bẋ + kx = 0.
2
b) From the above equation we identify the undamped natural frequency and damping
ratio as: r
2k b
ωn = , ξ=√ .
5m 10km
Thus, the damped natural frequency is:

r
p 10km − b2
ωd = ωn 1− ξ2 = .
25m2
For the given values of the parameters, this reduces to ωd = 1.264 rad/s.
Problem 17:
We model a nonuniform beam as a single-
degree-of-freedom system in the form:
mẍ + bẋ + kx = 0,
x1 = 1.00 m
bc
and experimentally measure the mass as m =
kg. In free vibration we experimentally de- x2 = 0.75 m
termine the equivalent spring constant to be bc
k = 4 N/m, and we measure the response as

shown.
a) Determine the equivalent damping con-
stant.
b) What is the exponential decay rate of

the transient solution?
Solution:
18
a) We use the logarithmic decrement so that:
µ ¶ µ ¶
δ1 x1 0.75 m
ζ= p , δ1 = − ln = − ln = 0.288,
4π 2 + δ12 x0 1.00 m
and we find that ζ = 0.0457. With m = 1 kg, the damping constant b is given as:
√
b = 2ζ km = 0.183 N/(m/s).
We note that with m given, we can also find the period of the oscillations to be:
2π 2π
T = = p ,
ωd ωn 1 − ζ 2
√ p 2
m 4π + δ12
= √ = 3.145 s.
k
b) The exponential decay rate is σ = ζωn , which is found to be σ = 0.0914 s−1 .
Problem 18:
The rigid beam (mass m = 2 kg, length l =
1.5 m) is supported by an elastic spring (k =
4 N/m) and damper (b = 2 N/(m/s)), and is
pinned to the ground.
k
a) Find the linearized equations of motion
in terms of z, the relative displacement (m, l)
between the end of the beam and the z
ground;
b) what is the frequency of the resulting b

motion;
c) if the mass of the spring is taken to be

mspring = 1 kg, find the new frequency
of the oscillations.
Solution:
a) We will define the inclination of the bar from the horizontal position as θ, so that the
angular acceleration is θ̈ k̂. Therefore, using angular momentum balance about the
point of rotation, we find:
mbar l2
θ̈ + bl ż + kl z = 0.
3
We assume that the coordinates θ and z are related by z = l sin θ, which for small
rotations reduces to z = lθ. Using this constraint to eliminate θ, the equation of motion
reduces to:
mbar
z̈ + b ż + k z = 0.
3
For the parameter values given above, this becomes:

µ ¶
2 ³ ´ ³ ´
kg z̈ + 2 N/(m/s) ż + 4 N/m z = 0.
3
19
For this system, the damping ratio and natural frequency can be expressed as:
r √
3k 3b
ωn = , ζ = √
mbar 2 kmbar
√
r
3
= 6 rad/s = .
8
p
b) The frequency of oscillation, ωd = ωn 1 − ζ 2 , reduces to:
s µ ¶2 √
3k 3b 15
ωd = − = rad/s = 1.94 rad/s.
mbar 2mbar 2
c) If the mass of the spring is considered, it is treated as an additional equivalent mass

meq = mspring /3 located at the end of the bar. Therefore the new moment of inertia of
the bar about the point of rotation is:
mbar l2
IO = + meq l2 ,
3
mbar l2 ³ mspring ´ 2
= + l ,
3 3
(mbar + mspring ) l2
= .
3
Therefore, the new frequency of oscillation is:

s ¶2
√
µ
3k 3b
ωd = − = 3 rad/s = 1.73 rad/s.
(mbar + mspring ) 2(mbar + mspring )
Problem 19:
The non-uniform beam supports an end-mass
of m = 10 kg. If the response of the system
is such that: m
x1 = 0.25 m x2 = 0.20 m,
t1 = 1.00 s t2 = 4.00 s,
as shown in the figure, find the equivalent (t1 , x1 )

stiffness and equivalent damping of the beam
bc
(assume that the beam is massless).
(t2 , x2 )
bc
Solution:
20
a) The logarithmic decrement is defined as δ = ln |x2 /x1 |. In terms of this quantity and
the period of oscillation T , the damping ratio and natural frequency are defined as:
p
(2π)2 + δ 2 δ
ωn = , ζ=p .
T (2π)2 + δ 2
Therefore, the stiffness and damping constants reduce to:
(2π)2 + δ 2 2δ
k = m ωn2 = m , b = m 2ζωn = m .
T2 T
For this system, we find that T = 4 s and δ = 0.223, and the stiffness and damping
constant reduce to:
k = 24.7 N/m, b = 1.12 N/(m/s).
Problem 20:
For the single-degree-of-freedom mechanical 2r k
system shown in the figure:
a) determine the governing equations of

motion;
r
b) what are the damping ratio and un-
damped natural frequency of this sys-
tem;
c) find the response of the system x(t) m

subject to the initial conditions x(0) =
x0 , ẋ(0) = 0,
c
when m = 1 kg, b = 12 (N · s)/m, and k =
9 N/m. Assume the pulley is massless and
neglect the effects of gravity.
Solution:
a) Let x represent the displacement of mass m in the vertical direction and θ measure the
angular displacement of the pulley, both measured from the static equilibrium position.
If T represents the tension in the cable supporting mass m, then T = kr 2 θ, where θ is
the angular displacement of the massless pulley from static equilibrium. In addition,
x
we find that θ is related to the linear displacement of mass m as θ = 2r . As a result, in
the ̂ direction, the equation governing the motion of the mass is:
k
mẍ + bẋ + x = 0.
4
As a result, we find that the the equivalent spring constant is keq = k4 .
b) The damping ratio and natural frequency are simply:

r √
b b keq k
ζ= p =√ , ωn = = √ .
2 keq m km m 2 m
21
In terms of the given parameters, we find ζ = 4.0 and ωn = 3 rad/s. We note that the
damping ratio has no units.
c) For an overdamped system, the general solution is:

³ ³ p ´ ³ p ´´
x(t) = exp(−ζωn t) c1 exp (ωn ζ 2 − 1)t + c2 exp −(ωn ζ 2 − 1)t ,
and with these initial conditions this reduces to:

Ã p
ζ + ζ2 − 1 ³ p ´
x(t) = x0 exp (−ζωn t) p exp (ωn ζ 2 − 1)t +
2 ζ2 − 1
p !
ζ − ζ2 − 1 ³ p
2
´
p exp −(ωn ζ − 1)t ,
2 ζ2 − 1
Ã √ √ !
−6t 4 + 15 3 √15t 4 − 15 − 3 √15t
= x0 e √ e2 + √ e 2 ,
2 15 2 15
= x0 e−6.0t 1.016e5.81t + 0.01640e−5.81t .
¡ ¢
Problem 21:
librium position, and the surface is friction- x
less.
a) Determine the governing equations of m

motion. 6k 2c
b) What is the period of each oscilla-
tion in terms of the system parameters
(m, k, c)?
c) For what value of c is the system criti-

cally damped?
d) If the system is released with the initial

conditions:
x(0) = 0.0 m, ẋ(0) = 5.0 m/s,
find the resulting solution x(t) if m =

2 kg, k = 48 N/m and c = 4 N/(m/s).
Solution:
a) In terms of x, the spring and damping forces can be written as:
Fspring = −6k x ı̂, Fdamper = −2c ẋ ı̂.
Using linear momentum balance on the block, we find that:

X F
F = (−6k x − 2c ẋ) ı̂ = mẍ ı̂ = m aG ,
22
and the equation of motion can be written in standard form as:
r
6k
ẍ + (2ζωn ) ẋ + (ωn2 ) x = 0, with ωn = ,
m
c
ζ=√ .
6km
b) In p
terms of ζ and ωn , the period of oscillation is simply T = (2π)/ωd , where ωd =
ωn 1 − ζ 2 . Therefore, for this system, the period reduces to:
2π m
T =r s µ ¶2 = 2π √6km − c2 ,
6k c
1− √
m 6km
provided c2 < 6km.
c) For a critically damped system, ζ = 1, and solving for c, yields:

√
ccritical = 6km.
d) With these parameter values, the equation of motion reduces to:
ẍ + 4 ẋ + 144 x = 0,
so that ωn = 12 rad/s, and ζ = 1/6. For these values, the general solution can be
written as: ³ √ ´
x(t) = A e−2 t sin 2 35 t + φ ,
where A and φ are arbitrary constants used to fit the initial conditions. Solving for A
and φ, we find: r
5
A= m, φ = 0 rad,
28
so that the general solution can be written as:

r
5 −2 t ³ √ ´
x(t) = e sin 2 35 t m.
28
23
Problem 22:
For the system shown to the right, x is mea-
sured from the unstretched position of the x
spring. Each block has mass m and the disk
has moment of inertia I and radius r. If the k (I, r)
gravitational constant is g: m
a) find the equations of motion which de-

termine x(t);
b) what is the period of the free oscilla-

tions?
Problem 23:
For the single-degree-of-freedom mechanical

motion;
b) what are the damping ratio and un- k 2k b

damped natural frequency of this sys- g
tem;
c) what is the stretch in the spring when

the system is in equilibrium?
m
Solution:
We assume that the ı̂ and ̂ directions are standard orthonormal basis in the horizontal
and vertical directions respectively.
a) The forces due to the springs in parallel and damping are:
Felastic = 3k · x̂, Fdamping = b · ẋ̂.
With the inclusion of the gravitational force, the equation of motion for this system can
be written:
X
F = m (−ẍ̂) ,
3k · x̂ + b · ẋ̂ − mg̂ = −mẍ̂.
Taking components in the ̂ direction, we obtain the governing equation of motion:
mẍ + bẋ + 3kx = mg.
24
b) Dividing through by the mass m, we find:
b 3k
ẍ + ẋ + x = g,
m m
so that:
b 3k
2ζωn = , ωn2 = ,
m m
which can be solved to yield:

r
3k b
ωn = , ζ= √ .
m 2 3mk
c) In equilibrium, the system is stationary, so that ẋeq = 0 and ẍeq = 0. Substitution into
the governing equations yields:
mg
3kxeq = mg, → xeq = .
3k
Problem 24:
a) find the linearized governing equations k

of motion for small θ;
b) find the frequency of oscillation for the (m, l)

G
response; θ
c) if the initial velocity is zero, and θ(0) =

θ0 , determine the time response of the
system. 4k
Solution:
a) The displacement of each end of the bar in the ̂ direction is:
l sin θ
x± = ± .
2
As a result, the total moment produced by the springs about the center of mass G is:
X l2
MG = − (4k sin θ + k sin θ) k̂,
4
5kl2
= − sin θk̂.
4
Thus angular momentum balance about G provides:

X
MG = I G θ̈k̂,
5kl2 ml2
− sin θk̂ = θ̈k̂.
4 12
25
For small angular displacements sin θ ∼ θ, and the governing equations of motion are
therefore:
15k
θ̈ + θ = 0.
m
b) This system is undamped. Therefore the frequency of the oscillation is equal to the
undamped natural frequency: r
15k
ω = ωn = .
m
c) This system possesses the general solution:
θ(t) = c1 sin ωt + c2 cos ωt,
which, for the initial conditions given above, yields the solution:
r
15k
θ(t) = θ0 cos t.
m
Problem 25:
We obtain the differential equation:
mẍ + bẋ + kx = 0,
as a model for a spring-mass-damper system with:
m = 2, k = 18.
a) Identify the damping constant b that gives rise to critical damping;
b) If, instead, b = 24, approximately how long will it take for the the amplitude of free vibration
to be reduced to within 2% of zero?
Solution:
a) Written in nondimensional form, the equations of motion are:
b k
ẍ + ẋ + x = 0,
m m
and so we identify the damping ratio and natural frequency as:

r
b k
ζ= √ , ωn = .
2 km m
ζ = 1 corresponds to critical damping, so that:

√
bcritical = 2 km = 12.
b) For b = 24 we find:
r
b k
ζ= √ = 2, ωn = = 3.
2 km m
26
Therefore the eigenvalues of this system are:
³ p ´ √
λ1,2 = ωn ζ ± ζ 2 − 1 = 3(−2 ± 3).
√
The dominant eigenvalue is λ = −6 + 3 3, and so the time t = τ required for the
amplitude of the free vibration to be reduced to within 2% of zero is:
−4
τ∼ √ ∼ 5.0
−6 + 3 3
Problem 26:
system shown in the figure, the bar has mass
2
m and length l (so that IO = ml3 ). If the k
spring is unstretched when θ = 0:
a) find the linearized governing equations ℓ

of motion for small θ; 2
ℓ
θ
2
b) find the frequency of oscillation for the (m, l)
free response;
c
Neglect gravity.
Solution:
a) Using angular momentum balance about the fixed point O, we find:
X
MO = IO αβ/F F,
l2 ml2
µ ¶
−kl2 sin θ − c θ̇ cos θ k̂ = θ̈k̂,
4 3
so that the equation of motion can be written as:
ml2 l2
θ̈ + c θ̇ cos θ + k l2 sin θ = 0.
3 4
Linearizing this equation about θ = 0, we obtain:

3c 3k
θ̈ + θ̇ + θ = 0.
4m m
b) The frequency of oscillation, that is, the damped natural frequency, is given as:
p
ωd = ωn 1 − ζ 2 ,
where ωn is the undamped natural frequency and ζ is the damping ratio. For this system
we find: r √
3k 3c
ωn = , ζ= √ ,
m 8 km
27
so that the damped natural frequency is:
s µ ¶2
3k 3c
ωd = − .
m 8m
Problem 27:
2
mass center G. g

x
system assuming that the disk rolls
without slip. c
b) If the disk is released from rest with G
initial displacement x(0) = x0 , find the
minimum value of the coefficient of fric- k (m, r)
tion for which the disk does not slip.
Problem 28:
librium position. If the surface is assumed to x
be frictionless:

motion; 6k b
b) what is the period of each oscillation;
c) what value of the damping coefficient b

corresponds to critical damping?
d) if k = 1 N/m and m = 4 kg, find the

displacement of the mass x(t) if the sys-
tem is critically damped and started
with the initial conditions x(0) = 0,
ẋ(0) = ẋ0 ;
Problem 29:
has moment of inertia I. If the gravitational k (I, r)
constant is g: m
a) what is the displacement of the spring

at static equilibrium;
b) find the kinetic energy of the system in

terms of the coordinate x (and/or its
velocity).
28
Problem 30:
2
2 .

r
tion; 2
b) What is the equivalent mass of the sys-

tem;
c) Find the frequency of oscillation for free m

vibrations of the system?
m
k
k c
Solution:
a) Define θ as the angular displacement of the disk in the −k̂ direction (clockwise), and
x1 and x2 as the displacements of the two blocks so that:
r
x1 = θ, x2 = rθ.
2
We define the tension in the left and right cable as T1 and T2 respectively. Thus, linear
momentum balance on the two blocks, and angular momentum balance on the disk
yield:
r mr2
T1
− T2 r = − θ̈,
2 2
T1 − kx1 − cẋ1 = mẍ1 ,
T2 + kx2 = −mẍ2 .
Eliminating the two unknown tensions from these three equations, we find that the
equation of motion (on θ) can be reduced to:
7mr2
µ 2¶
5kr2
µ ¶ µ ¶
cr
θ̈ + θ̇ + θ = 0.
4 4 4
2
b) Examination of the above equation shows that meq = 7mr 4 . Note that this answer is
not unique, but depends on the equation of motion. For example, if we had written the
equation of motion as 7m θ̈ + c θ̇ + k θ = 0, the equivalent mass would be meq = 7m.
c) The frequency of oscillation is given by ωd , which reduces to:

r r
p 5k c2
ωd = ωn 1 − ζ 2 = 1− .
7m 140km
29
Problem 31:
In the spring-mass-damper system shown,
the block slides with no friction. With m =
2 kg, k = 18 N/m, b = 13 (N · s)/m: x
a) Find the resulting solution if the sys-
b
tem is released from the unstretched
position with initial velocity ẋ(0) = m
−10.0 m/s k
b) Identify the damping constant b that
gives rise to critical damping;
Solution:
For this system the differential equation of motion can be written as:
mẍ + bẋ + kx = 0,
b k
ẍ + ẋ + x = 0,
m m
subject to specified initial conditions.
a) With the given values for the mass, stiffness, and damping coefficient, the above equation
becomes:
ẍ + (6.5 kg/s)ẋ + (9 kg/s2 )x = 0,
whose characteristic equation reduces to:

13
λ2 + λ + 9 = 0,
2
which has two real solutions of the form:

−13 ± 5
λ= .
4
Thus the system has two purely real eigenvalues and the resulting system is overdamped
and decays exponentially with no sustained oscillations. The general solution is given
as:
9
x(t) = c1 e− 2 t + c2 e−2t .
For the initial conditions x(0) = 0 m, ẋ(0) = 10.0 m/s, we find that:
c1 = 4 m, c2 = −4 m,
and the general solution to this equation, subject to these initial conditions, becomes:
³ ´ ³ −1 −1
´
x(t) = 4 m e−(9/2 s )t − e−(2 s )t .
b) For this system, the damping ratio is:

b
ζ= √ .
2 km
30
Thus, for a critically damped system ζ = 1, and solving for b with m = 2 kg and
k = 18 N/m, we find:
bcr = 12 (N · s)/m.
Problem 32:
2
mass center G. g

x
system assuming that the disk rolls
without slip. b
b) What value of b correspond to critical G
damping?
k (m, r)
c) Find the displacement of the center of
the disk when the system is critically
damped and released from rest with
x(0) = x0 .
Solution:
We define θ as the angular displacement of the disk in the k̂ direction as measured
from the equilibrium position of the disk (when the spring is unstretched). Assuming
that the disk rolls without slip, the rotation and translation of the disk can be related
through the constraint equation:
x = −rθ.
a) The frictional force, which is unknown, is defined as f = fı̂, while the forces due to the
spring and damper are:
Fspring = −kx ı̂, Fdamper = −bẋ ı̂.
Using linear and angular momentum balance on the disk, we find that:
X F
F = (f − kx − bẋ) ı̂ = mẍ ı̂ = m aG ,
X mr2
MG = (f r) k̂ = θ̈ k̂ = IG αβ/F F.
2
Eliminating the unknown frictional force, and using the kinematic constraint, we find
the equation of motion is: µ ¶
3m
ẍ + bẋ + kx = 0.
2
b) We can write this differential equation in standard form, that is:

r
2k
ẍ + (2ζωn )ẋ + (ωn2 )x = 0, with ωn = ,
3m
b
ζ=√ .
6km
31
If the system is critically damped, then this implies that the damping ratio is unity.
Therefore, solving for b when ζ = 1, we find:
√
bcritical = 6km.
c) When the system is critically damped, the general solution takes the form:
x(t) = c1 + c2 t e−ωn t ,
¡ ¢
while this can be differentiated with respect to time to obtain the velocity:
ẋ(t) = (c2 − ωn c1 ) − (ωn c2 ) t e−ωn t .

¡ ¢
If the system is released from rest when a known initial displacement, then the initial
conditions are (x(0), ẋ) = (x0 , 0). Thus, returning these to the general solution, we find:
x0 = x(0) = c1 ,
0 = ẋ(0) = c2 − ωn c1 .
Solving for c1 and c2 , the solution to these initial conditions becomes:
x(t) = x0 1 + ωn t e−ωn t ,
¡ ¢
p
where recall that ωn = (2k)/(3m).
Problem 33:
and moment of inertia IG about the mass cen-
ter, and the applied moment has a constant g
magnitude M k̂. If the disk rolls without slip
(µ is sufficiently large): x
M k̂
a) determine the governing equations of c
motion;
G
b) what are the equivalent mass, stiffness,
and damping of the system; k (m, r)
c) what is the stretch in the spring when

the system is in equilibrium?
Solution:
a) The governing equations of motion are:
µ ¶
I M
m + 2 ẍ + cẋ + kx = − .
r r
b) The equivalent mass, damping, and stiffness are:
I
meq = m + , ceq = c, keq = k.
r2
32
c) When the system is in equilibrium, the displacement of the disk is:
M
xeq = − .
kr
2 Frictionally Damped Systems
Problem 34:
The block shown to the right rests on a rough ̂
surface with coefficient of friction µ and
m = 6 kg, k = 128 N/m. ı̂
a) If the block is displaced 3 cm to the g

right and released, for what values of µ
will the block remain in that position? x z
A B
b) With µ = 0.50, if the block is displaced k
30 cm to the right and released from k
m
rest, how long will it take the block to
come to rest?
µ
Solution:
In addition to the variable x identified in the problem statement, we also define z to
be the stretch in the spring parallel with the cable system. As a one degree-of-freedom
system, the variables x and z are directly related. The relative velocity across the spring
can be identified as
F F
vB − vA = ż ı̂,
= (−ẋ ı̂) − (ẋ ı̂) ,
so that ż = −2 ẋ. Therefore the kinematic −m g ̂

relationship becomes
z = −2 x. T ı̂
−k x ı̂
An appropriate free-body diagram for this
T ı̂
system is shown to the right. Note that the
unknown friction force is denoted as fr ı̂ and fr ı̂
the tension in the cable is T . Finally, exam-
N ̂
ining the spring in the cable, the tension T
and the displacement z are related as
−T ı̂ T ı̂
T = k z = −2 k x.
Applying linear momentum balance to the block yields

X F
F = (2 T − k x + fr ) ı̂ + (N − m g) ̂ = m ẍ ı̂ = m aG ,
and in terms of x the equation of motion becomes

m ẍ + 5 k x = fr .
33
If the block slips then fr = −µ m g sgn(ẋ) while is sticking occurs |fr | ≤ µ m g.
a) If the block is in static equilibrium at a displacement x = xeq , then ẍeq ≡ 0 and the
equation of motion reduces to
5 k xeq = fr ,
so that equilibrium is maintained provided
|fr | = |5 k xeq | ≤ µ m g.
This inequality is satisfied provided

µmg
|xeq | ≤ .
5k
Problem 35:
For the system shown to the right, the block
slides on a rough surface (coefficient of fric-
tion µ) inclined at an angle of φ with respect
to vertical. If the block is subject to a peri-
odic force of the form g
k
F (t) = F0 sin(ω t),

neglect gravity; φ m
F (t)
b) find the amplitude of the steady-state
response using Mc when
m = 1.25 kg, k = 20 N/m,

φ = 30◦ , µ = 0.125,
F0 = 4 N, ω = 2.00 rad/s,
34
Problem 36:
[(Spring 2003)] The system shown in the fig-
ure has mass m and rests on a plane inclined
at an angle φ. The coefficient of friction for
the rough surface is µ and the system is re-
leased from rest at the unstretched position k x ̂2
of the spring (with stiffness k).
a) If µ = 0, what is the equilibrium dis-
µ ı̂2
placement of the mass (as measured
̂ m
from the unstretched position)?
θ
b) For µ > 0, at what angle θ does the ı̂
block begin to slip?
c) Find the value of θ so that the system

comes to rest after one full cycle exactly
at the equilibrium position of the sys-
tem found in part a (so that the friction
force vanishes when the system comes
to rest), with:
m = 2 kg, k = 32 N/m,
µ = 0.35,
Solution:
The unit directions ı̂2 and ̂2 are defined to be coincident with the inclined plane and
the coordinate x represents the displacement of the mass from the unstretched position
of the spring, as shown in the figure.

shown to the right. Notice that the force in −k x ı̂2 −m g ̂
the upper spring depends on z, rather than
x, while the friction force has an unknown
magnitude f . Because the disk is assumed
to roll without slip, we are unable to specify
the value of f , but instead can relate the dis-
placement and rotation of the disk through
f ı̂2
the coordinate relations above. N ̂2
Therefore, linear momentum balance yields the following equations:

X F
F = m aG ,
³ ´ ³ ´ ³ ´
f − k x ı̂2 + N ̂2 − m g ̂ = m ẍ ı̂2 ,
³ ´ ³ ´
f − k x + m g sin θ ı̂2 + N − m g cos θ ̂2 = m ẍ ı̂2 .
Therefore, this leads to the following scalar equations in the ı̂2 and ̂2 directions:
m ẍ + k x = f + m g sin θ,
N = m g cos θ.
a) If µ = 0, then the friction forces vanishes and the first of the above equations reduces
35
to:
m ẍ + k x = m g sin θ.
The equilibrium displacement of the mass, xeq , then can be found to be:
mg
xeq = sin θ.
k
b) With µ 6= 0, an equilibrium state is maintained provided:

¯ ¯
¯f ¯ ≤ µ N,
¯ ¯
¯ ¯
¯k x − m g sin θ¯ ≤ µ m g cos θ.
¯ ¯
Therefore, if the system is released from x = 0, the block begins to slide when:
tan θ = µ.
c) Define z to be a new coordinate measuring the displacement of the system from static
equilibrium:
mg
z =x− sin θ,
k
so that the equations of motion become:

m z̈ + k z = f, N = m g cos θ.
with initial displacement z(0) = −(m g sin θ)/k. Over one complete cycle of motion,
for a frictionally damped system the amplitude decreases by a value:
4µN
∆A=− .
k
Therefore, if the system comes to rest at exactly the equilibrium position, then this
decrease in amplitude must exactly match the initial displacement. That is:
¯ ¯
¯ 4 µ N ¯ ¯¯ m g ¯
|∆ A| = ¯¯− ¯ = ¯− sin θ ¯ = |z(0)|.
¯
k ¯ k
Solving for θ:
tan θ = 4 µ.
Problem 37:
The spring mass system rests on a surface
with coefficient of friction µ and x is mea-
g
spring.. If the initial conditions of the sys-
tem are chosen to be x(0) = 0 and ẋ(0) = ẋ0 ,
find the range of ẋ0 so that the system comes
m
to rest after exactly one cycle of motion. k
36
Problem 38:
For the spring-mass system with Coulomb
damping:
x
motion;
m
b) what is the period of each oscillation.
k k
Solution:
a) We measure the displacement of the mass from the static equilibrium of the frictionless
F
system, i.e., µ = 0, so that the acceleration of the block is aG = ẍı̂. Thus linear
momentum balance yields:
mẍı̂ = Fspring ı̂ + fµ ı̂ + (N − mg)̂.
The spring force is Fspring = −2kx, while the force due to sliding friction opposes the
velocity and is simply:
ẋ
fµ = −µmg ,
|ẋ|
since the normal force balances the gravitational force, i.e., N = mg. If the block
is stationary the magnitude of the frictional force is less than µmg. Therefore, the
governing equations of motion are:
fµ = −µmg |ẋẋ| |ẋ| =

6 0,
mẍ + 2kx = fµ , with
|fµ | ≤ µmg ẋ = 0.
b) Coulombic damping does not effect the frequency of oscillation, which is simply:
r
2k
ω= .
m
Therefore the period of the oscillation is:

√
2π 2π m
T = = √ .
ω 2k
37
Problem 39:
librium position. If the coefficient of friction x
is µ:

motion; 6k 2k
b) what is the period of each oscillation; µ
c) if the system is released from rest with
x(0) = x0 > 0, what is the minimum
value of x0 so that the block slips;
d) find the range of initial displacements

so that the system comes to rest after
one complete cycle.
Problem 40:
For the spring-mass system with Coulombic
damping, x is measured from the unstretched x
position of the spring. If the coefficient of
friction is µ and the gravitational constant is k I=0
g: m

motion; µ
b) if the system is released from rest in the

unstretched position (x(0) = 0, ẋ(0) =
0), for what values of µ will the system
move;
m
c) what is the displacement (from the un-
stretched position) of the upper block
when it first comes to rest?
Solution:
Notice that x describes the displacement of both masses and, since the pulley is massless,
the tension in the string connecting the masses is constant, say T . Also, notice that in
part c we ask for the displacement from the unstretched position of the spring, rather
than from equilibrium. Therefore we include the gravitational force which will influence
this result.
a) With the frictional force defined as F = f ı̂, linear momentum balance on the upper
and lower block yields:
mẍ + kx = f + T,
mẍ = mg − T,
where the frictional force is defined as:

ẋ
f = −µmg |ẋ| , ẋ 6= 0,
|f | ≤ µmg, ẋ = 0.
38
Eliminating the unknown tension T , the equation of motion is given as:
2m ẍ + kx = f + mg,
where f is defined as above and depends on the motion of the system, that is, the value
of ẋ.
b) If the system is released from rest in the unstretched position, it will remain there
provided the magnitude of the frictional force is less than µmg—the transition to move-
ment occurs when |f | = µmg. Thus the system does not move is |f | ≤ µmg and ẍ = 0.
Substituting these conditions into the equations of motion we find:
|f | = |kx − mg| ≤ µmg,
which, solving for µ with x(0) = 0, implies that the system does not move if µ ≥ 1.
Therefore, the system does move when:
µ < 1.
c) The displacement of the upper block when it first comes to rest is:
2(1 − µ)mg
x1 = .
k
This can be found by either solving the equations of motion explicitly, or through a
work-energy analysis. Since the initial and final kinetic energy is zero, the work done
by the frictional force balances out the change in potential energy from the spring and
gravity.
Problem 41:
For the spring-mass system with Coulombic
damping, x is measured from the unstretched
position of the spring. If the coefficient of x
g
friction is µ and the gravitational constant is
g:

motion; k
b) if the system is released from rest, so

that ẋ(0) = 0, for what range of initial µ
displacements (from the unstretched
position) will the block come to rest
when the block first comes again to rest
(ẋ(t1 ) = 0 for t1 > 0)?
Solution:
a) The equations of motion can be written as:
mẍ + kx = f,
39
where f is the force due to friction, modeled by Coulomb’s law of friction as:
ẋ
f = −µmg , ẋ 6= 0,
|ẋ|
|f | ≤ µmg, ẋ = 0.
b) If the system is released from rest, the initial displacement must be sufficiently large so
that the block slides, rather than remaining at rest. Sliding does not occur if the force
due to friction is sufficient to balance the elastic force, that is, µmg ≥ f = kx(0). Thus,
solving for x(0) we find, that for sliding to occur:
µmg
|x(0)| > .
k
However, if |x(0)| is too large, the system will undergo multiple reversals as the ampli-
tude of the motion decays. Consider the block sliding to the left (ẋ < 0), released from
rest with initial displacement x(0) > µmgk . Thus the equation of motion becomes:
mẍ + kx = µmg,
which has the general solution:

 
µ ¶
µmg k µmg
x(t) = x(0) −
  cos t + .
k m k
Therefore, when the block comes again to rest at time t1 (unknown), the mass is at the
position:
µmg
x(t1 ) = 2 − x(0).
k
µmg
At this point, the block sticks if and only if |x(t1 )| ≤ k . Therefore, solving for x(0),
we find that:
µmg
x(0) ≤ 3 .
k
So for a block with x(0) > 0, the allowable range for x(0) is µmg k < x(0) ≤ 3 µmg
k .
Together with an identical argument for x(0) < 0 yields the total allowable range as:
µmg µmg
< |x(0)| ≤ 3 .
k k
40
Problem 42:
has moment of inertia I and radius r. The k (I, r)
coefficient of friction between the upper block m
and the table is µ. If the gravitational con-
stant is g:
µ
a) find the equations of motion which de-
termine x(t);
b) what is the minimum value of µ so that

the system slips when release from rest
with x(0) = 0; m
c) what is the period of the free oscilla-

tions?
d) if the system is released from rest, what

is the range of initial displacements
x(0) so that the systems comes to rest
after exactly one complete cycle?
Solution:
We begin by defining two additional coordinates, θ, which describes the rotation of
the disk in the −k̂ direction (clockwise), and y which measures the displacement of
the hanging mass in the −̂ direction. These additional coordinates are related to the
displacement of the upper mass by the constraint equations:
x
y = x, θ= .
r
a) On each mass the equations of motion can be written as:

X ³ ´ ³ ´
F
F = − k x + T1 + f ı̂ + N − mg ̂ = mẍ ı̂ = m aG1 ,
X ³ ´
MO = T1 r − T2 r k̂ = −I O θ̈k̂ = I O αD/F F,
X ³ ´
F
F = T2 − mg ̂ = −mÿ ı̂ = m aG2 ,
Notice that I O = I 6= 0, so that provided θ̈ 6= 0 the tensions T1 and T2 are not equal.
Taking the components of these equations and eliminating the unknowns (T1 , T2 ), while
using the constraint equations, we find that the equation of motion for this system
reduces to: µ ¶
I
2m + 2 ẍ + k x = f + mg,
r
where f the value of the frictional force in the ı̂ direction, can be written as:
½ ẋ
−µmg |ẋ| , ẋ 6= 0
f=
f0 , |f0 | ≤ µmg, ẋ = 0.
41
b) The minimum value for slip is simply µmin = 0. If we would like to find the range of
µ for which slip occurs, we resort to the value of f at static equilibrium. Assuming
(ẋ, ẍ) = (0, 0), the equations of motion reduce to:
kx = fstatic + mg,
where fstatic represents the force required to maintain static equilibrium. Solving for
this quantity and using the frictional inequality, we find:
|fstatic | = |kx0 − mg| ≤ µmg.
Therefore, solving for µ yields:

¯ ¯ ¯ ¯
¯ kx0 ¯ ¯ kx0 ¯¯
µ ≥ ¯¯ − 1¯¯ = ¯¯1 − ,
mg mg ¯
which provides a necessary condition for sticking at x = x0 . So for sliding to occur for
x0 = 0, this implies that µ < 1.
c) The period of oscillation for a frictionally damped system is identical to that of an

undamped system. Therefore:
s
2π 2m + rI2
T = = 2π .
ωn k
d) Let δ describe the displacement of the system from equilibrium. The amplitude of
oscillation will decay by a value of ∆ = 4µmg/k over one cycle of motion. Therefore,
Since the system will come to rest within the range:
µmg µmg
− < |δfinal | ≤ ,
k k
the initial displacement δ0 from the equilibrium in the absence of friction must be in
the range:
µmg 3µmg 5µmg µmg
− +∆= < |δ0 | ≤ = + ∆.
k k k k
However, the equilibrium position corresponds to xeq = mg

k , and so the allowable range
of x0 is:
3µmg ¯¯ mg ¯¯ 5µmg
< ¯x0 − ¯≤ .
k k k
3 Forced Single Degree-of-freedom Systems
42
Problem 43:
For the system shown to the right the bar
of length ℓ has mass m and is subject to a
time-dependent moment of the form
M (t) = M0 sin(ω t) k̂. M (t)

k
a) Find the equations of motion for this G
rotation of the bar.
b) What is the steady-state amplitude of (m, ℓ)

k b
the forced response, with
m = 3 kg, b = 16 N/(m/s),
1
k = 8 N/m, ℓ = m,
4
1
ω = 4 rad/s, M0 = N·m
16
Problem 44:
The block of mass m = 20 kg shown to the
right rests on a rigid foundation and is sub- f (t) ̂
ject to a time-dependent load
F (t) = f0 sin(ω t) ̂.
a) Design an undamped foundation to m

achieve isolation ≥ 33% for all forcing
frequencies ω > 3 π rad/s;
b) If, for the isolator that you designed,

the damping ratio was measured to be
ζ = 0.125 (rather than ζ = 0 as as-
sumed above), what is the minimum
isolation achieved over this frequency
range?
43
Problem 45:
In the figure shown to the right, the disk is
subject to a time dependent moment of the k
form M (t) k̂
M (t) = M0 sin(ω t).
r
a) Find the equations of motion for the
angular displacement of the disk.
b) With I r/2
k = 280 N/m, b = 12 N/(m/s),

m = 4 kg,
I = 0.40 kg · m2 , r = 0.10 m
M0 = 3 N · m, ω = 5 rad/s,
m
Determine the steady-state response of
disk as a function of time. k b
Problem 46:
The unbalanced rotor shown in the figure is
pinned to a frame and supported by a spring
and damper. The total mass of the structure
is m, with a small rotating component (10% k b
of the total mass) offset by a distance r from
the center of rotation C.
a) Find the distance between the center r

of mass of the system and the center of
rotation; C
b) what is the damped natural frequency;
m
c) determine the steady-state amplitude
of vibration when the rotor spins at an
angular speed of ω = 5 rad/s with:
k = 256 N/m, b = 12 N/(m/s),

m = 5 kg, r = 10 cm
44
Problem 47:
2
and moment of inertia IG = mr 2 about the g
mass center G and is assumed to roll without
slip. The attached plate undergoes harmonic
motion of the form u(t)
c
u(t) = u0 sin(ω t).
G
a) Find the equation of motion in terms of
the displacement between the moving (m, r) k
plate and the center of the disk;
b) What are the damping ratio and natu-

ral frequency for this system?
c) If the system is critically damped, find

the amplitude of the relative displace-
ment of the disk for
m = 3 kg, r = 0.10 m
k = 36 N/m, c = 3 N/(m/s)
u0 = 0.05 m, ω = 5 rad/s
Problem 48:
right, the uniform rigid bar is massless and
pinned at point O while a force is applied at ℓ
A of the form k 3 B
O
f (t) = t e−σ t .
θ
For this system: A
2ℓ
3 m

frequency in terms of the parameters f (t) k c
m, c, k, and ℓ.
c) With
m = 2 kg, ℓ = 30 cm,
c = 0.25 N/(m/s), k = 50 N/m,
−1
σ = 2.00 s ,
find the convolution integral for the re-

sponse of the system. You need not
evaluate the integral.

small.
45
Problem 49:
The block shown to the right rests on a rough
surface with coefficient of friction µ (assume
that any cables can support compression and x
tension). Find the amplitude of the vibra-
tions of the block if f (t) = f0 cos(ω t), with
k
m
m = 4.0 kg, k = 64 N/m,
µ = 0.05,
f0 = 40 N, ω = 4 rad/s.
−f (t) ̂
Problem 50:
and damper. The total mass of the structure
is m, with a small rotating component (10% k b
of the total mass) offset by a distance r from
the center of rotation C.
a) find the damped natural frequency; r

b) what is the steady-state amplitude of
vibration when the rotor spins at an an- C
gular speed of ω = 5 rad/s with:
m
k = 108 N/m, b = 9 N/(m/s),
m = 3 kg, r = 7.5 cm
46
Problem 51:
length ℓ, and is pinned at point O. A har-
monic force is applied at A. For this system:
k O
B θ
ℓ
frequency in terms of the parameters A 2
m, c, k, and ℓ.
c) For:
4m
m = 6 kg, ℓ = 25 cm,
f0 sin(ω t)
c = 0.50 N/(m/s), k = 80 N/m,
c
f0 = 2.00 N, ω = 10 rad/s,
find the steady-state amplitude of the

displacement of the block.

small.
Problem 52:
2
and moment of inertia IG = mr 2 about the g
mass center G and is assumed to roll without
slip. The attached plate undergoes harmonic
motion of the form u(t)
c
u(t) = u0 sin(ω t).
G
a) Find the equation of motion in terms of
the angular rotation of the disk; k (m, r)
b) What are the damping ratio and natu-

ral frequency for this system?
c) If the system is critically damped, find

the amplitude of the rotation of the disk
for
m = 3 kg, k = 36 N/m,
u0 = 10 cm, ω = 3 rad/s
47
Problem 53:
and damper. The total mass is measured as
m = 200 kg. When the system is operated at k b
ω = 25 rad/s the phase φ of the response with
respect to the rotation of the unbalanced disk
is measured to be π/2 rad and the steady-
state vibration amplitude is X = 2.00 cm. ε
When the rotation rate of the disk is much G
larger than this value the amplitude reduces C
to X = 0.50 cm.
Find the stiffness and damping constant for m
the foundation and the distance between the
center of rotation C and the mass center G.
Problem 54:
right, the uniform rigid bar is massless and
pinned at point O while a harmonic force is ℓ
applied at A. For this system: k 3 B
O
θ
b) Identify the damping ratio and natural A
frequency in terms of the parameters 2ℓ
3 m
m, c, k, and ℓ.
c) For:
f0 sin(ω t) k c
m = 2 kg, ℓ = 30 cm,
c = 0.25 N/(m/s), k = 50 N/m,
f0 = 2.00 N, ω = 10 rad/s,
find the steady-state displacement of

the block.
d) What is the magnitude of the force
transmitted to the ground through the
sping and damper attached to the
block? (Do not include the spring at-
tached at A.)

small.
48
Problem 55:
The block shown to the right rests on a rough
surface with coefficient of friction µ and the
block is subject to a compressive force of N = x
20 N (do not include gravity, just this normal
load and assume that any cables can support
compression and tension). k
m (m, r)
a) If f (t) = f0 = constant, find the range
of initial displacements for which the
block will remain stationary if released
−f (t) ı̂
from rest (it will stick).
b) Find the amplitude of the vibrations of

the block if f (t) = f0 cos(ω t), with
m = 4.0 kg, k = 64 N/m,

r = 12.5 cm, µ = 0.50,
f0 = 40 N, ω = 4 rad/s.
Problem 56:
displacement of the disk from the unstretched k
position of the spring. The surface is inclined x
at an angle of φ with respect to vertical.

neglect gravity;
b) what is the static equilibrium displace- φ (m, r) M (t)

ment of the disk?
c) if the disk is subject to a periodic mo-

ment
M (t) = M0 sin(ω t)k̂,
find the displacement of the center of

the disk if it is released from rest at the
unstretched position of the spring.
49
For the system shown to the right, the disk of x
mass m rolls without slip and the inner hub ̂
has radius ρ/2. z
2k ı̂
a) Find the equations of motion (in terms
of the given parameters—do not substi- θ
tute in numerical values yet);
b G
b) If the applied moment takes the form: k
M (t) (m, ρ)
M (t) = (2 N · m) sin(4 t),
C
find the steady-state amplitude of the
translation of the center of the disk
when:
k = 16 N/m, b = 2 N/(m/s),
m = 2 kg, ρ = 0.125 m
c) Determine the steady state amplitude

of the friction force.
Solution:
a) We identify the three coordinates x, z, and −m g ̂
θ as shown in the figure above. These are
related as: M (t) k̂
3
x = −ρ θ, z = x. −2 k z ı̂
2
An appropriate free-body diagram for this −b ẋ ı̂ −k x ı̂
system is shown to the right. Since the disk
is assumed to roll without slip, the equation
of motion can be directly obtained with an-
gular momentum balance about the contact fr ı̂
point C
X N ̂
MC = I C αD/F ,
which yields
3 m ρ2
µ ¶
3ρ
ρ (k x + b ẋ) + (2 k z) + M (t) k̂ = θ̈ k̂.
2 2
Using the above coordinate transformations this equation can be written as
3m 11 k M (t)
ẍ + b ẋ + x=− .
2 2 ρ
b) For the numerical values given above (with consistent units), this equation reduces to
3 ẍ + 2 ẋ + 88 x = −8 sin(4 t),
50
from which we can identify the appropriate parameters as:
r r
88 1 6
ωn = , ζ=√ , r= .
3 264 11
Therefore, the amplitude of the translational oscillations becomes

F 8 1 1
X= M(r, ζ) = r =√ .
k 88 ¡ 6 2
¢ ³
1
q
6
´ 2 26
1 − 11 + 2 √264 11
Likewise, the phase shift of the response is:

q
1 6
2ζ r 2 √264 11 1
tan φ = = 6 = ,
1 − r2 1− 11
5
so that φ = 0.20 rad = 11.3◦ .
c) In the development of the equation of motion, the friction force was eliminated by
summing moments about C. Using linear momentum balance we can reintroduce the
friction force as
X ³ ´ ³ ´
F
F = fr − k x − 2 k z − b ẋ ı̂ + N − m g = m ẍ ı̂ = m aG .
Therefore, solving for fr we find that
fr = m ẍ + b ẋ + 4 k x.
With x(t) represented as x(t) = X sin(ω t − φ), where X and φ are found above, the
friction force becomes
³ ´ ³ ´
fr (t) = 4 k − m ω 2 X sin(ω t − φ) + b ω X cos(ω t − φ).
The magnitude of the friction force is then found to be

° ° r³ ´2 ³ ´2
°f ° = X 4 k − m ω2 + b ω .
° °
For these parameter values kf k = 6.47 N.
51
and damper. If the total mass is m while the
mass center G is located at an eccentricity of k b
ε from the the center of rotation O,
a) find the damped natural frequency;

ε
vibration when the rotor spins at an an- G
gular speed of ω = 8 π rad with: O
k = 32 N/m, b = 16 N/(m/s), m
m = 4 kg, ε = 2.5 cm
c) If the system is undamped (i.e., b =

0 N/(m/s)) for what range of operat-
ing speeds (ω) will the amplitude of the
force transmitted to the ground FT be
less than 1 N.

a) find the equations of motion (in terms f (t) = sin(ω t)

of the given parameters—do not substi-
tute in numerical values yet);
ℓ
2
b) if c = 0.25 N/(m/s), k = 32 N/m,
k θ
m = 2 kg, and ℓ = 0.25 m, find the m
amplitude of the force transmitted to z
the ground through combination of the
spring and damper when ω = 4 rad/s. ℓ
2 ̂
c) if c = 0, m = 2 kg, and ℓ = 0.25 m, c
find the value of the stiffness k so that ı̂
the bar’s amplitude of oscillation is less
than π/6 rad for all forcing frequencies
greater than 20 rad/s.
Solution:
a) In addition to θ, we define the additional coordinate z, which measures the deflection
at the left end of the bar, with θ and z related as:
52
ℓ
z= θ
2
A free-body diagram for this system is f (t) ̂
shown to the right. Applying angular mo- k z1 ̂
mentum balance on the bar eliminates the G
appearance of the reaction force and leads
to:
FR
c ż1 ̂
X
MG = I G θ̈ k̂,
³ ´ℓ m ℓ2
f (t) − k z − cż k̂ = θ̈ k̂.
2 12
Solving for z in terms of θ, the equation of motion becomes:
m ℓ2 c ℓ2 k ℓ2 ℓ
θ̈ + θ̇ + θ = f (t).
12 4 4 2
b) In standard form, this equation of motion can be written as:

µ ¶ µ ¶ µ ¶
3c 3k 6
θ̈ + θ̇ + θ= sin(ω t),
m m mℓ
so that: r √
3k 3c 6
ωn = , ζ= √ , M0 = .
m 2 km mℓ
The amplitude of the moment transmitted to the ground can be written as:
p
1 + (2 ζ r)2
MT = (meq M0 ) p ,
(1 − r2 )2 + (2 ζ r)2
q ¡ ¢2
1 + c kω
p
ℓ ℓ k 2 + (c ω)2
= q¡ = q¡ .
2 2 2 2 2 2 ¢2 2
1 − m3 ωk + c kω k − m 3ω
¢ ¡ ¢
+ (c ω)
The amplitude of the force transmitted to the ground is then FT = MT /(ℓ/2), or:
p
k 2 + (c ω)2
F T = q¡ .
2 ¢2 2
k − m 3ω + (c ω)
Substituting in the numerical values given in the problem statement, we find that:
FT = 1.50 N
c) The amplitude of the steady-state vibrations can be written as:
M0 1
Θ = ,
ωn2
p
(1 − r ) + (2 ζ r)2
2 2
2 1 2 1
= q¡ = q¡ .
kℓ m ω 2 ¢ 2 ¡ c ω
¢ 2 ℓ m ω2 2
¢ 2
1− 3 k + k k− 3 + (c ω)
53
Substituting in the numerical values given in the problem statement, we find that:
8
Θ = ¯¡ 2 ω2
¢¯
¯ k− ¯
3
Therefore, if the amplitude of vibration is less than π/6:

8 π
¯ 2 ω2 ¯
¯ ≤ ,
¯k −
3
6
2 ω 2 ¯¯
¯ ¯
48 ¯
≤ ¯¯k − .
π 3 ¯
This inequality has two solutions:
48 2 ω 2 2 ω2 48
k≥ + , k≤ − .
π 3 3 π
Since this condition must be satisfied for all ω ≥ 20 rad/s, we take the second inequality
and find that:
2 (20)2 48
k≤ − = 251.
3 π

For the system shown to the right, the disk of x ̂
displacement of the disk from the unstretched 3k z
ı̂
a) find the equations of motion; θ f (t)
b) if the forcing takes the form: k m

½
f0 , 0 ≤ t < t0 ,
f (t) =
f0 /2, t0 ≤ t,
find the response of the system with

zero initial conditions.
Solution:
a) We define the three coordinates as shown as the figure, related as:
x = −r θ, z = −2 r θ, z = 2 x.
A free-body diagram for this system is −3k z ı̂

shown to the right. Notice that the force
in the upper spring depends on z, rather G
than x, while the friction force has an un-
known magnitude fr . Because the disk is −k x ı̂ C f (t) ı̂
assumed to roll without slip, we are unable
to specify the value of fr , but instead can fr ı̂
relate the displacement and rotation of the
disk through the coordinate relations above.
54
The equations of motions can be developed directly with angular momentum balance
about the contact point, so that:
X
MC = I C θ̈ k̂,
³ ´ 3 m r2
(3k z) 2r + (k x) r − f (t) r k̂ = θ̈ k̂.
2
Finally, writing this equation in terms of a single coordinate, we obtain:
3 m r2
µ ¶
θ̈ + 13 k r2 θ = r f (t).
¡ ¢
2
In standard form: µ ¶
26 k 2 f (t)
θ̈ + θ= .
3m 3mr
b) We use the convolution integral to determine the response, so that:

Z t
θ(t) = F (τ ) h(t − τ ) dτ,
0
and for this system:

r Ãr !
f (t) 2 f (t) 1 3m 26 k
F (t) = = , h(t) = sin(ωn t) = sin t .
meq 3mr ωn 26 k 3m
Because the forcing function changes abruptly at t = t0 , the solution must be written
separately for 0 < t ≤ t0 , and t > t0 :
2 f0 t sin(ωn (t − τ ))
Z
x(t) = dτ, 0 < t ≤ t0 ,
3m 0 ωn
2 f0 t0 sin(ωn (t − τ ))
Z t
f0 sin(ωn (t − τ ))
Z
x(t) = dτ + dτ, t > t0 .
3m 0 ωn 3 m t0 ωn
Evaluating these integrals, the solution becomes:

2 f0 ³ ´
x(t) = 1 − cos(ωn t) , 0 < t ≤ t0 ,
26 k
f0 ³ ´
x(t) = 1 + cos(ωn (t − t0 )) − 2 cos(ωn t) , t > t0 .
26 k
55
̂
The system shown to the right is subject to
base excitation. Find the steady-state re- m
sponse of the system in terms of z, with: ı̂
m = 2.0 kg, b = 4.0 N/(m/s),

k1 k1
k1 = 3.00 N/m, k2 = 12.00 N/m.
b z
x
k2
u(t) = 0.50 sin(2 t) m
Solution:
a) We define the addition coordinate x which measures the absolute displacement of the
mass with respect to the ground, so that:
x = z + u(t).
Notice that the collection of springs can be

replaced by a single equivalent spring, with: m
1 2 k1 k2
keq = 1 1 = = 4 N/m.
2 k1 + k2
2 k1 + k2 keq b z
The new equivalent system is shown to the

right. u(t)
An appropriate free-body diagram is shown

to the right. In terms of the identified coor-
dinates, the acceleration of the mass center
is:
F
aG = ẍ ̂ = (z̈ + ü) ̂, −keq z ̂ −b ż ̂
with ü(t) = −(u0 ω 2 ) sin(ω t).
Therefore, linear momentum balance on the mass yields:

X F
F = m aG ,
³ ´
− keq z − b ż ̂ = m ẍ ̂
Writing this in terms of z, the equation of motion is:
m z̈ + b ż + keq z = −m ü(t),
and in standard form:
z̈ + 2 ζ ωn ż + ωn2 z = u0 ω 2 sin(ω t),
56
with: r
keq b
ωn = , ζ= p .
m 2 keq m
Therefore the steady state response of this system becomes:
z(t) = Z sin(ω t − ψ),
with Z = u0 Λ(r, ζ), and:
r2 2ζ r ω
Λ= p , tan ψ = , and r =
(1 − r2 )2 + (2 ζ r)2 1 − r2 ωn
For the numerical values of this problem:

√ 1 √
ωn = 2, ζ=√ , r= 2,
2
so that:
2 2
Λ= √ , tan ψ =
5 −1
Recall that the phase shift ψ must be positive, so that is tan ψ is negative, then ψ is in
the second quadrant, so that ψ = 3.03 rad. Finally:
1
z(t) = √ sin(2 t − 2.03).
5
Problem 62:
2
mass center G, subject to an applied moment g
of the form:
x
M (t) = M0 sin(ωt)k̂. M (t)
c
a) If the coefficient of friction is µ, find the
condition which determines if the disk G
rolls with or without slip, and find the (m, r)
governing equations of motion when the k
disk rolls with and without slip;
b) If µ is sufficiently large so that the

disk rolls without slipping, what are the
equivalent mass, stiffness, and damping
of the system;
c) If again the disk rolls without slipping,

find the amplitude of the steady-state
motion of G, that is, x(t), when the sys-
tem is critically damped, with M0 =
4 N · m, m = 1 kg, r = 0.1 m, and
k = 2 N/m.
57
Problem 63:
and moment of inertia IG = α m about the
mass center G. The disk is subject to a time- g
varying force f (t) = 4 sin(ω t).
a) Find the equations of motion for this x

system assuming that the disk rolls c
without slip.
f (t)
b) After the transient solution decays, find G
the amplitude of the force transmit- k (m, r)
ted to the ground through the spring-
damper element.
c) For what value of the damping ratio is

this transmitted force less than twice
the applied load for all values of the
forcing frequency?
Problem 64:
The disk shown in the figure rolls without
slip and is subject to a time-varying moment
M (t) = sin(t). g

tion; x
M (t) k̂
b) Find the frequency of oscillation for the b
unforced response, i.e. M (t) = 0; G
c) What is the steady-state amplitude of (m, r)
k
the forced response?
d) Determine the amplitude of the fric-

tional force during the steady-state mo-
tion.
Solution:
a) We assume that (ı̂, ̂, k̂) represent the standard orthonormal basis, while the transla-
tional displacement of G is xı̂ and the angular displacement is θk̂. Linear and angular
momentum balance on the disk yield:
mẍ = −kx − bẋ − fµ ,
mr2
θ̈ = −fµ r − M (t),
2
where fµ is the unknown frictional force. If the disk rolls without slipping, we find the
kinematical relation x = −rθ, and, eliminating fµ from the above balance laws, we find
the governing equation of motion can be expressed as:
3m M (t)
ẍ + bẋ + kx = .
2 r
b) For M (t) = 0, the frequency of oscillation is the damped natural frequency ωd =
58
p
ωn 1 − ζ 2 , where:
r r
2k 2 b 1
ωn = = , ζ=√ =√ ,
3m 3 6km 2
so that the damped natural frequency is ωd = √1 = 0.577.

3
c) In nondimensional form, the system is represented as:
ẍ + 2ζωn ẋ + ωn2 x = F0 sin(t),
where ζ and ωn are given above, and:

2 2
F0 = = .
3m 3
Recall that the magnification factor for a harmonically driven system is:
1
M= p ,
(1 − r2 )2 + (2ζr)2
q
where r = ωωn is the frequency ratio. For this system we find r = 3
2, so that M = √2 .
13
As a result, the steady-state amplitude is:
F0 2
A= M = √ = 0.555.
ωn2 13
d) The response of the disk is x(t) = A sin(t + φ), where A is given above and:
√
−2ζr − 3
tan φ = = .
1 − r2 − 21
From the expression for angular momentum balance, we find:
m M (t)
fµ = ẍ − ,
2 r
1
= − √ sin(t + φ) − sin(t),
13
·µ ¶ µ ¶ ¸
1 1
= − √ sin φ cos(t) + √ cos φ + 1 sin(t) ,
13 13
√
√1 2 3
where cos φ = 13
and sin φ = √
13
, so that the amplitude of the frictional force is:
4
|fµ | = √ .
13
59
Problem 65:
2
mass center G, subject to an applied moment g
of the form:
x
M (t) = M0 sin(ωt)k̂. M (t)
b
If the disk rolls without slip (µ is sufficiently
large): G
a) find the equations of motion for this k (m, r)

system;
b) for m = 2 kg, b = 0.5 (N · s)/m, and

k = 8 N/m, find the steady-state am-
plitude of the rotation of the disk?
Solution:
In addition to x, defined in the figure, we define θ as the angular displacement of the
disk from the unstretched position in the k̂ direction. If the disk rolls without slip, x
and θ are related as:
x = −rθ.
a) The frictional force, which is unknown, is defined as f = fı̂, while the forces due to the
spring and damper are:
Fspring = −kx ı̂, Fdamper = −bẋ ı̂.
Using linear and angular momentum balance on the disk, we find that:
X F
F = (f − kx − bẋ) ı̂ = mẍ ı̂ = m aG ,
X mr2
MG = (M0 sin(ωt) + f r) k̂ = θ̈ k̂ = IG αβ/F F.
2
Eliminating the unknown frictional force, and using the kinematic constraint, we find
the equation of motion is:
µ ¶
3m M0
ẍ + bẋ + kx = − sin(ωt).
2 r
b) For the given values of the parameters, we find that:

8 1 M0
ωn2 = , ζ= √ , F0 = − .
3 8 6 3r
Thus, the steady-state amplitude may be easily found as:

F0
X = M,
ωn2
−M 0
3r
= q¡ ¢2 ¡ ω ¢ 2 .
8
3 − ω2 + 6
60
Problem 66:
The mass m = 2 kg is supported by an elas- ̂
tic cantilever beam attached to a founda-
tion which undergoes harmonic motion of the m ı̂
form:
u(t) = 4 sin(ωt) m,
If the beam has length l = 20 cm, while x z
AE = 16 N:
a) find the equations of motion in terms

of z, the relative displacement between u(t) = 4 sin(ω t) m
the mass and the foundation (assume
the beam has zero mass);
b) what is the amplitude of the result-

ing motion in terms of the forcing fre-
quency ω?
Solution:
a) The equivalent spring for this cantilever beam is:
AE 16 N
keq = = = 80 N/m.
l 0.2 m
F
The acceleration of the block with respect to the ground is aG = (ü + z̈)̂, so that, in
terms of z, the equations of motion become:
mz̈ + keq z = −mü,

z̈ + 40z = 4ω 2 sin(ωt).
b) For this undamped system, the amplitude of the resulting steady-state motion is:
4ω 2 1
X = q¡ ¢2 ,
40 ω2
1− 40
4ω 2
= .
|40 − ω 2 |
61
Problem 67:
̂
ı̂

ε
vibration when the rotor spins at this G
angular speed. O x
m
Solution:
We define x(t) as the vertical displacement of the geometric center of the rotor as
measured from static equilibrium. As a result, the mass center G is described by the
position z(t) = x(t) + ε sin(ωt). Note that ε measures the eccentricity of the mass
center, not the location of the mass imbalance. Consequently, the governing equations
of motion can be written:
mz̈ = −kx − bẋ,

2
mẍ − εω sin(ωt) = −kx − bẋ,
or in more standard form:

b k
ẍ + ẋ + x = εω 2 sin(ωt).
m m
q
k √b ,
a) In the above system we find ωn = m and ζ = 2 km
so that the damped natural
frequency can be written:
r
p k b2
ωd = ωn 1− ζ2 = − .
m 4m2
b) For an arbitrary forcing frequency the amplitude of oscillation is A = εΛ, where:
ω2
Λ= p ,
(ωn2 − ω 2 )2 + (2ζωωn )2
which, with ω = ωd , reduces the amplitude to:
1 − ζ2
A=ε p ,
ζ 4 − 3ζ 2
with ζ defined above.
62
Problem 68:
The system shown in the figure is supported
by a foundation that undergoes harmonic mo- ̂
tion of the form: g
m
u(t) = 4 sin(ωt) m, ı̂
If the mass and stiffness are m = 2 kg, and

k = 8 N/m: k z
x
of z, the relative displacement between
the mass and the base (assume the u(t) = 4 sin(ω t) m
spring has zero unstretched length);
b) what is the amplitude of the resulting

motion in terms of the frequency ratio
r;
c) for what forcing frequencies is the re-

sulting amplitude of the steady-state
motion Z ≤ 8?
Solution:
a) We construct a free-body diagram as shown. The only forces acting on the mass arise
from the gravitational force and the spring force.
However, the acceleration of the mass with respect to

the ground is written as: −m g ̂
³ ´
F
aG = ü + z̈ ̂.
Therefore, linear momentum balance takes the form:

X F
F = m aG , −k z ̂
³ ´ ³ ´
− kz − mg ̂ = m z̈ + ü ̂.
Substituting in the expression for u(t), and taking the component in the ̂ direction, we
obtain the governing equation of motion:
mz̈ + kz = −mg − mü,

z̈ + (4 s−2 )z = −9.81 m/s2 + 4ω 2 sin(ωt) m,
Notice that each term has units of acceleration, that is, m/s2 . In what follows we will
dispense with the explicit inclusion of the units. For this system ωn = 2, ζ = 0, and the
amplitude of the forcing is F = 4ω 2 . Thus we have frequency-squared excitation.
b) Using the above values for the natural frequency and the damping ratio, we find that
63
the forced response can be written as z(t) = Z sin(ωt + φ), where the amplitude Z is:
r2
Z = U · Λ(r, ζ) = Up ,
(1 − r2 )2 + (2ζr)2
4r2
= .
|1 − r2 |
c) If Z < 8 then this implies that:

4r2
Z<p .
(1 − r2 )2
Defining Zcr as:

4r2
Zcr = p ,
(1 − r2 )2
We may solve for r to yield:

s
2 ± 4Z
Zcr cr
r = 2 − 16
,
Zcr
(r r )
Zcr Zcr
= , .
Zcr − 4 Zcr + 4
Thus for Zcr = 8, we find that Z < 8 in the range:

√
r
2
0<r< , or r > 2.
3
Problem 69:
The single-degree-of-freedom system shown is
subject to a harmonic force. If the natural
frequency is ωn = 4 rad/s and m = 1 kg: x F (t) = sin(ω t) N
a) determine the spring and damping con-
stants when the system is critically m
damped; k b
b) determine the amplitude of the total
force transmitted to the ground un-
der steady-state oscillations when ω =
1 rad/s.
Solution:
The equation of motion for this system takes the form:
b k F (t)
ẍ + ẋ + x = .
m m m
a) In terms of the system parameters, the damping ratio and natural frequency can be
written as: r
b k
ζ= √ , ωn = .
2 km m
64
So, for a critically damped system ζ = 1, and solving for k and b we find:
k = 16 N/m, b = 8 N · s/m.
b) The amplitude
p of the total force transmitted to the ground, which is defined as Ft , is
Ft = Xωn2 1 + (2ζr)2 , where X is the amplitude of the response and r = ωωn = 14 is
the frequency ratio. Thus, for this system:
Ã !
F0 1 p
Ft = 2
p · ωn2 1 + (2ζr)2 ,
ωn (1 − r2 )2 + (2ζr)2
p √
F0 1 + (2ζr)2 8 5
= p = = 1.05.
(1 − r2 )2 + (2ζr)2 17
Problem 70:
A constant force is applied to the undamped
single degree-of-freedom system for a dura-
tion of t1 , at which point it is removed, that x
is:
½ F (t)
F0 ı̂, 0 ≤ t < t1 , m
F (t) =
0, t ≥ t1 . 4k
If the system starts with zero initial condi-
tions, determine the resulting displacement
of the mass x(t). You may either use the
convolution integral or you may try to solve
this explicitly.
Solution:
The equation of motion for this system is:
4k F (t)
ẍ + x= ,
m m
which has an impulse response of the form:

r
sin ωn t k
h(t) = , ωn = 2 .
mωn m
Therefore, using the convolution integral, the response of the system is:
Z t
x(t) = F (ξ)h(t − τ ) dτ,
0
 Z t
F0 sin(ωn (t − τ ))


 0 m dτ, 0 ≤ t < t1 ,
ωn

=
 Z τ F sin(ω (t − τ ))
 0 n

 dτ, t ≥ t1 .
0 m ωn
F0 ³
 ´
 1 − cos(ω n t) , 0 ≤ t < t1 ,
mωn2



=
 F0 ³ ´
cos(ω (t − t )) − cos(ω t) , t ≥ t1 .


n 1 n
mωn2

65
Problem 71:
A mass m = 2 kg is rigidly connected to a
rigid massless bar of length ℓ = 40 cm, which
is pinned to a wall. The system is supported
by a combination of springs and damper as ̂ b
shown in the figure, and subjected to a time- M (t)
m
dependent moment M (t) = 12 sin(t). If k = ı̂
12 N/m and b = 6 N/(m/s):
ℓ
a) find the undamped natural frequency 2 z2
and the damping ratio of the system; θ z1
b) what are the steady-state amplitude
and phase of the forced response? 6k
k 2k
Problem 72:
The system shown in the figure is supported
by a foundation that undergoes an exponen-
tially decaying motion of the form: g
m
u(t) = 16e−t/4 .
If the mass and stiffness are m = 5 kg, and

k = 45 N/m: k z

of z, the relative displacement between
the mass and the base (assume the u(t)
spring has zero unstretched length);
b) find the resulting solution z(t) if the

mass is started from rest and the spring
is initially unstretched (assume g =
10 m/s2 and use the convolution inte-
gral);
c) what is the amplitude of the oscillations
as t increases (i.e. u(t) → 0)?
66
Problem 73:
In the system shown at right, the disk is as- 2r k
sumed to be massless while:
m = 2 kg, k = 4 N/m,
b = 14 N/(m/s), r = 1 cm. r
a) Determine the governing equations of
motion;
b) What is the period of oscillation;
c) If f (t) = sin t, find the amplitude and m

phase of the resulting motion as t → ∞.
f (t)
b
Problem 74:
In the mechanical system shown each spring
is identical, with spring constant k. If m = ̂
f (t)
1 kg:
ı̂
a) determine the spring constant k and
damping constant b so that the un-
forced system is critically damped and m
the exponential rate of decay is τ =
2 s−1 ;
b x
b) with the initial conditions:
x(0) = 0, ẋ(0) = 2 m/s,

k
find the resulting solution of x(t) for the

unforced problem, i.e., f (t) = 0;
c) with f (t) = sin(2t), determine the am-

plitude of the force transmitted to the
supporting structure.
67
Problem 75:
̂
ı̂
b) for:
ε
G
m = 4 kg, ε = 0.1 m,
O x
b = 0.5 (N · s)/m, k = 8 N/m
m
what is the steady-state amplitude of
vibration when the rotor spins at this
angular speed? 2k
Problem 76:
The disk shown in the figure has mass m =
4 kg, and is subject to a time-dependent mo- 4k
ment M (t) = M0 sin(t) k̂. If k = 2 N/m and
r = 0.2 m:
IG
a) Find the steady-state amplitude of the r
response when: b
b = 1 N/(m/s), M0 = 12 N · m.
b) What is the amplitude of the oscil- r

lations when the system is critically 2
damped?
2k
68
Problem 77:
For the system shown to the right, the block
slides on a rough surface (coefficient of fric-
tion µ) inclined at an angle of φ with respect
to vertical. If the block is subject to a peri-
odic force of the form g
k
F (t) = F0 sin(ω t),

neglect gravity; φ m
F (t)
b) find the amplitude of the steady-state
response using Mc when
m = 1.25 kg, k = 20 N/m,

φ = 30◦ , µ = 0.125,
F0 = 4 N, ω = 2.00 rad/s,
4 Multi Degree-of-freedom Systems
Problem 78:
For the system shown to the right
a) use Lagrange’s equations to determine

the equations of motion;
b) find the mode shapes and natural fre-

m 4m
quencies of the system;
2k k F (t) ı̂
c) determine the modal equations for this
system.
69
)
,d
Problem 79:
m
(4
b) are your equations statically coupled,

dynamically coupled, both or neither?
2 d)
k
k (m,
f (t)
m
2k b
Problem 80:
a) use Lagrange’s equations to determine

the equations of motion;
4k k
b) find the mode shapes and natural fre-
3m
quencies of the system. Normalize the
mode shapes so that u T M u = 1. k k
m
Problem 81:
For the system shown to the right the bar
of length ℓ is massless and the block on the
right is subject to a time-dependent force of
the form k
k
F (t) = F0 sin(ω t).
G
The intermediate springs are located a dis-
x2 k
tance 2ℓ from the center of the bar.
x1
a) Find the equations of motion in terms
of x1 and x2 . m
m
b) Determine the natural frequencies and
mode shapes of this system.
c) Determine the modal equations. F (t)
d) What is the steady-state amplitude of

the in-phase motion?
70
right, the uniform rigid bar is supported by
identical springs. For this system:
a) find the equations of motion in terms k

of x and θ (and in terms of the given k
parameters—do not substitute in nu- x
merical values yet); θ ℓ
G 3
b) if k = 64 N/m, m = 2 kg, and ℓ = m, ℓ
0.25 m, find the natural frequencies and
mode shapes for this system;

For the system shown to the right: F0 sin(ω t)
a) find the equations of motion (and in m βm

terms of the given parameters—do not k αk
substitute in numerical values yet);
b) for α = β = 2, find the natural frequen-

cies and mode shapes of the system and
normalize the mode shapes by M ;
c) determine the equations of motion that

describe the response of each mode.

In the figure shown to the right, in the ab-
sence of gravity the springs are unstretched x2
in the equilibrium position. Determine the k2
equations of motion for this system.
m IG
r2
b θ
r1
x3
z
k1
̂
ı̂ m
x1
k1
71
Solution:
Because we can, we define five different coordinates to describe the dynamical behavior
of this two degree-of-freedom system, leading to the following transformations:
x2 = −r2 θ, x3 = −r1 θ, z = x3 − x1 .

shown to the right. We develop three equa-
−(k2 x2 + b ẋ2 ) ı̂
tions of motion based on linear momentum
−T ı̂
balance on both blocks and angular momen-
tum balance on the disk: T ı̂
Block 1:
X F
F = m a G1 ,
FR
³ ´
k1 z − k1 x1 ̂ = m ẍ1 ̂
Block 2: −k1 z ̂
X F
F = m a G2 ,
³ ´ k1 z ̂
T − b ẋ2 − k2 x2 ı̂ = m ẍ2 ı̂
Disk:
X
MG = I G αD/F ,
³ ´ −k1 x1 ̂
T r2 + k1 r1 z k̂ = I G θ̈ k̂
From the equations on block 2 and the disk, we eliminate the unknown tension T from
the system to obtain:
I G θ̈ − m r2 ẍ2 − b r2 ẋ2 − k2 r2 x2 − k1 r1 z = 0.
From this equation, and the equation on block 1, we eliminate the coordinates z and
x2 , and obtain the equations of motion to be:
m ẍ1 + 2 k1 x1 + k1 r1 θ = 0,
³ ´ ³ ´
I G + m r22 θ̈ + b r22 θ̇ + k1 r1 x1 + k1 r12 + k2 r22 θ = 0.
Problem 85:
For the system shown in the figure: x
a) what is the degree-of-freedom for this k k (I, r)
system?
2m
b) using Lagrange’s equations, determine
the differential equations that govern
the motion.
72
Solution:
a) This system contains three masses which are each allowed to move in only one direc-
tion. The upper mass slides horizontally with displacement x1 , while the disk rotates
through an angle θ. Finally, the suspended mass moves vertically and its position can
be described by the coordinate x2 . However, because the disk and the suspended mass
are connected by an inextensible string, their motion can be related by:
x2 = rθ.
So this system has only two independent coordinates and therefore it is a two-degree-
of-freedom system.
b) We utilize the coordinates x1 , x2 , and θ, as measured from the unstretched position of

the two springs. Therefore, the kinetic and potential energies are written as:
1 1 1
T = (2m) ẋ21 + m ẋ22 + I θ̇2 ,
2 2 2
1 1
V = k x21 + k (rθ − x1 )2 − m g x2 .
2 2
However, x2 and θ are related by the above relationship. Thus eliminating θ, the energies
become:
1 1 1 I 2
T = (2 m) ẋ21 + m ẋ22 + ẋ ,
2 2 2 r2 2
1 1
V = k x21 + k (x2 − x1 )2 − m g x2 ,
2 2
which, using Lagrange’s equations, yields the equations of motion:

2m ẍ1 + 2 k x1 − k x2 = 0,
µ ¶
I
m + 2 ẍ2 − k x1 + k x2 = m g.
r
Problem 86:
For the system shown in the figure:
k2
a) find the mass and the stiffness matrix;
k1 m2
b) is your system of equations dynamically
coupled, statically coupled, or both? m1 k3 k3
m3
Solution:
We choose the coordinates (x1 , x2 , x3 ), which represent the positions of the three masses.
a) To determine the stiffness matrix, we use the stiffness influence coefficients. Maintaining
a unit displacement of each mass in turn requires forces of the form:
(x1 , x2 , x3 ) = (1, 0, 0) → f = (k1 + k2 , −k2 , 0)T ,
(x1 , x2 , x3 ) = (0, 1, 0) → f = (−k2 , k2 + 2k3 , −2k3 )T ,
(x1 , x2 , x3 ) = (0, 0, 1) → f = (0, −2k3 , 2k3 )T .
73
Therefore, with these coordinates the stiffness matrix is:
 
k1 + k2 −k2 0
K =  −k2 k2 + 2k3 −2k3 
0 −2k3 2k3
Alternatively, we can define the potential energy of the system as:

1 1 1
V = (k1 )(x1 )2 + k2 (x2 − x1 )2 + (2k3 )(x3 − x2 )2 ,
2 2 2
1 1 1
= (k1 + k2 )(x1 )2 − (2k2 )(x1 x2 ) + (k2 + 2k3 )(x2 )2
2 2 2
1 1
− (2k3 )(x2 x3 ) + (2k3 )(x3 )2 ,
2 2
which leads to the same stiffness matrix.
To determine the mass matrix, we could use the inertia influence coefficients, but, for
variety, we determine the kinetic energy as:
1 1 1
T = m1 ẋ21 + m2 ẋ22 + m3 ẋ23 .
2 2 2
Therefore, the mass matrix is:

 
m1 0 0
M = 0 m2 0 
0 0 m3
b) With this choice of coordinates, the mass matrix is diagonal and the stiffness matrix
contains nonzero off-diagonal terms. Thus, the system is statically coupled but dynam-
ically uncoupled.
Problem 87:
For the system shown at right:
a) determine the governing equations of k2

motion; k1 k1
m1 k2 m1
b) is the system statically or dynamically m2
coupled, or both;
c) find the matrix M −1 K .
74
Problem 88:
In the system shown to the right, the mass
and stiffness matrix are:
· ¸
m1 0
M = ,
0 m2
· ¸
k1 + k2 −k2
K =
−k2 k2
while the damping is proportional, i.e., C =

βM +αK . If m1 is subject to harmonic forc-
ing F (t) = sin t, with:
m1 = 2.0 kg, m2 = 1.0 kg, β = 0.1,

k1 = 6.0 N/m k2 = 2.0 N/m, α = 0.1
a) determine the eigenvectors of this sys-

tem, normalized by the kinetic energy
inner product;
b) what are the damped natural frequen-

cies of this system;
c) find the steady-state amplitude of vi-

bration of the mode of vibration with
the lowest natural frequency.
Problem 89:
The two-degree-of-freedom system shown is
subject to a harmonic force applied to the
f (t)
block of mass 2m, of the form:
2m 2m
¡ ¢
f (t) = (2 sin(t)) N ı̂.
k k k
If the system is subject to proportional damp-
ing with α = 0.25, m = 2 kg, and k =
4 N/m, find:
a) the mass, damping, and stiffness matri-

ces;
b) the forced, damped equation (single-

degree-of-freedom) that describes the
motion of each mode;
c) the steady-state response of the system.

Solution:
a) Using influence coefficients, we find that the mass and stiffness matrices are:
· ¸ · ¸
2m 0 2k −k
M = , K = .
0 2m −k 2k
Therefore, with proportional damping, the damping matrix becomes C = αK , so that
75
we find: · ¸ · ¸
1 0 2 −1
M = (4 kg) , K = (4 N/m) ,
0 1 −1 2
· ¸
2 −1
C = (1 N/(m/s)) .
−1 2
Problem 90:
In the multi-degree-of-freedom system shown
in the figure, the block with mass 4m slides x1
on a smooth, frictionless surface. If the pulley z
is massless:
2k k r
a) using Lagrange’s equations, determine
the differential equations governing the 4m
motion, as measured from static equi-
librium;
b) with m = 1 kg and k = 16 N/m,

find the natural frequencies and mode
shapes for the free vibration of this sys-
tem. Normalize the mode shapes so x2
that with respect to the mass matrix
2m
the amplitude of each mode is one;
c) find the general solution to these equa-

tions for the above values of m and k.
Solution:
a) We identify the three coordiantes x1 , x2 , and z, with
z = x2 − x1 .
Measuring the response from static equilibria and neglecting the gravitational potential
energy, the kinetic and potential energies for this system can be written as
1¡ 1¡
4 m ẋ21 + 2 m ẋ22 ,
¢ ¢
T =
2 2
1¡ ¢ 2 1¡ ¢ 2
V = 2 k x1 + k z ,
2 2
In terms of x1 and x2 , the potential energy becomes

1 ¡ ¢ 2 1 ¡ ¢¡ ¢2 1¡ ¢ 2 1¡ ¢ 1¡ ¢ 2
V= 2 k x1 + k x2 − x1 = 3 k x1 + − 2 k x1 x2 + k x2 .
2 2 2 2 2
Therefore the equations of motion become
4 m ẍ1 + 3 k x1 − k x2 = 0,
2 m ẍ2 − k x1 + k x2 = 0,
or in matrix form
· ¸ · ¸ · ¸ · ¸
4 0 ẍ1 3 −1 x1
m +k = 0.
0 2 ẍ2 −1 1 x2
76
b) The corresponding eigenvalue problem for the above system is
· 3
− 14
¸
¡ −1 ¢ k 4
k ¡ ¢
M K u =λu −→ = βu
m − 21 1
2 m
The characteristic equation is

µ ¶µ ¶
3 1 1 5 1
−β − β − = β 2 − β + = 0,
4 2 8 4 4
with the solution

½ ¾
5±3 k k
β= −→ λ = ω2 = , .
8 4m m
Returning this to the eigenvalue problem, the mode shapes are defined by the equation
3 1
ui1 + ui2 = βi ui1 ,
4 4
so that
· ¸ · ¸ · ¸
1 1 1
ui = ci u1 = c1 , u2 = c2 .
4 βi − 3 −2 1
Normalizing ui by the mass matrix implies that

· ¸ · ¸
4m 0 1
1 = uiT M ui = c2i
£ ¤
1 (4 βi − 3) .
0 2m (4 βi − 3)
Solving for ci
s r r
1 1 1
ci = −→ c1 = , c2 =
2 m (2 + (4 βi − 3)2 ) 12 m 6m
Finally, the normalized eigenpairs are

" # " #
√ 1 √1
r r
k 12 m k 6m
ω1 = , u1 = 2 , ω2 = , u2 = √1
.
4m − √12 m
m 6m
c) With the above mode shapes and natural frequencies the general solution becomes
· ¸ 2
x1 (t) X
q (t) = = (Ai sin(ωi t) + Bi cos(ωi t)) ui ,
x2 (t)
i=1
Ã Ãr ! Ãr !! " #
k k √ 1
= A1 sin t + B1 cos t 12 m
2
4m 4m − √12 m
Ã Ãr ! Ãr !! " #
√ 1
k k 6m
+ A2 sin t + B2 cos t √1
.
m m 6m
77
Problem 91:
For the system shown in the figure:
a) find the mass and the stiffness matrices; k3

k1 m3
b) is your system of equations dynamically k2 k1
coupled, statically coupled, or both? m1
m2
Solution:
a) For coordinates we choose (x1 , x2 , x3 ) as the displacements of each mass with respect
to inertial space. Using influence coefficients, we find that:
   
m1 0 0 k1 + k2 −k2 0
M =  0 m2 0  , K =  −k2 k1 + 2k2 −k2 
0 0 m3 0 −k2 k2
Alternatively, if we choose coordinates (x1 , x2 , z), where z represents the stretch in the
spring connecting m2 and m3 , we can determine the mass and stiffness matrices from
the Lagrangian. The kinetic and potential energies are:
1 1 1
T = m1 ẋ21 + m2 ẋ22 + m3 (ẋ2 − ẋ3 )2 ,
2 2 2
1 1 1 1
V = k1 x1 + k2 (x2 − x1 )2 + k2 z 2 + k1 x22 .
2
2 2 2 2
Thus, with these coordinates the mass and stiffness matrices become:
   
m1 0 0 k1 + k2 −k2 0
M =  0 m2 + m3 −m3  , K =  −k2 k1 + k2 0 
0 −m3 m3 0 0 k2
b) With the first choice of coordinates, the mass matrix is diagonal while the stiffness
matrix is not, the system is statically coupled but dynamically uncoupled. With the
latter coordinates neither matrix is diagonal so that the system is both statically and
dynamically coupled.
78
Problem 92:
2
2 .
a) What is the degree-of-freedom for this

r
system? 2
b) Find the governing equations of mo-

tion;
k
c) If m = 1 kg and k = 4 N/m, what
are the frequencies of oscillation for
the motion and the corresponding mode m
shapes, normalized by the kinetic en- m
ergy inner product?
k
Solution:
a) Let the displacement of the left block, disk, and right block be described as (−x1 ̂),
(θ k̂), and (x2 ̂) respectively. Although x1 and θ are related by the following constraint:
r
x1 = θ,
2
x2 is independent from the above two coordinates. Therefore, the system has two
degrees-of-freedom.
b) With the above coordinates, the kinetic and potential energies can be written as:
1 1 1 mr2 2
T = mẋ21 + mẋ22 + θ̇ ,
2 2 2 2
1 2 1
V = kx + k(x2 − rθ)2 ,
2 1 2
Thus, using the above kinematic constraint to eliminate θ, the Lagrangian becomes:
L = T − V,
· ¸
1 £ 2 2
¤ 1 2 2
= m 3ẋ1 + ẋ2 − k 5x1 − 4x1 x2 + x2 .
2 2
Using Lagrange’s equations of motion, the governing equations are:
3m ẍ1 + 5k x1 − 2k x2 = 0,
m ẍ2 − 2k x1 + k x2 = 0.
c) From the above equations, the mass and stiffness matrices can be written as:
· ¸ · ¸
3 0 5 −2
M =m , K =k
0 1 −2 1
79
The characteristic matrix, A = M −1 K becomes:
· 5
− 32
¸
k 3
A= ,
m −2 1
and the characteristic equation can be written as:

8 1
β2 − β + = 0,
3 3
k
where, if β is a solution to this equation, λ = m β is an eigenvalue of the characteristic
matrix A. This quadratic equations has solutions of the from:
√
4 ± 13
β= .
3
To determine the eigenvectors, we return to the characteristic matrix A, so that Au =

λu. The elements of u then satisfy the equation:
5 2
u 1 − u 2 = β u1 .
3 3
Thus, if u1 = 1, this yields:

√ √
1∓ 13 4± 13
u2 = for β= .
2 3
Normalizing by the kinetic energy inner product, we find that:

u u
û = p =p
(u, u)M 3u21 + u22
Problem 93:
For the system shown in the figure, the sur-
face is assumed to be frictionless.
a) Using Lagrange’s equations, determine k

the differential equations governing the 4m
motion (measured from static equilib- (m, r)
rium);
b) Find the mode shapes and natural fre-

quencies of the resulting motion.
k
80
Problem 94:
The multi-degree-of-freedom system shown in
the figure is subject to harmonic forcing of
the form f (t) = sin(4 t).
a) Find the equations of motion in terms k

of the coordinates x1 , x2 , and x3 , and
identify the mass and stiffness matrices.
m
b) If m = 2 kg and k = 4 N/m, the natu-
ral frequencies of the system are found
to be:
2k 2k
ω1 = 0.38 rad/s, ω2 = 2.62 rad/s,
ω3 = 1.41 rad/s.
m m
Find the corresponding mode shapes
and normalize them so that (ui , ui )M =
1. f (t)
c) Determine the forced, uncoupled equa-

tions of motion for the modal ampli-
tudes Qi (t).
Problem 95:
For the system shown in the figure, the sur-
face is assumed to be frictionless. If each
block is displaced by a distance d (down and x1
to the right), find the resulting motion of the
system.
k (m, r)
m
θ
k z
m
x2
Solution:
We define the coordinates x1 , x2 , θ, and z as shown in the figure, so that
x1 = −r θ, z = r θ − x2 .
81
With these coordinates, the kinetic and potential energies can be written as
1 ¡ ¢ 2 1 m r2
µ ¶
1¡ ¢ 2
T = m ẋ1 + θ̇2 + m ẋ2 ,
2 2 2 2
1¡ ¢ 2 1¡ ¢ 2
V= k x1 + k z .
2 2
Expressing these only in terms of the coordinates x1 and x2 , we obtain

µ ¶
1 3m 1¡ ¢ 2
T = ẋ21 + m ẋ2 ,
2 2 2
1¡ ¢ 2 1¡ ¢ 2 1¡ ¢ 2 1¡ ¢ 1¡ ¢ 2
V= k x1 + k (−x1 − x2 ) = 2 k x1 + 2 k x1 x2 + k x2 .
2 2 2 2 2
Therefore, the mass and stiffness matrix can be identifed as

· 3 ¸ · ¸ · ¸
0 2 1 x1 (t)
M =m 2 , K =k with q (t) =
0 1 1 1 x2 (t)
and the equations of motion are:

· 3 ¸ · ¸ · ¸ · ¸
2 0 ẍ1 2 1 x1
m +k = 0.
0 1 ẍ2 1 1 x2
The solution to this equation requires the solution of an eigenvalue problem of the
form · 4 2 ¸
k
(M −1 K ) u = 3 3 u = λ u,
m 1 1
which is determined from the characteristic equation

·µ ¶ ¸
³
−1
´ k 4 2
det M K − λ I = − β (1 − β) − = 0,
m 3 3
k
with λ = m β. This quadratic equation has the solution
½ ¾ (r r )
7±5 1 k 2k
β= = , 2 −→ ω = ,
6 3 3m m
With this, the eigenvectors are determined by returning to the original eigenvalue equa-
tion · 4 2 ¸ · ¸ · ¸
3 3 ui1 ui1 3 βi − 4
= βi , −→ ui2 = ui1
1 1 ui2 ui2 2
In addition, normalizing the eigenvectors by the mass matrix

¤ 3m 0 6 + (3 βi − 4)2
· ¸ · ¸ µ ¶
ui1
1 = uiT M ui = ui1 ui2 = m u2i1
£
2 .
0 m ui2 4
Solving for ui1 the normalized eigenvectors are

" # " #
√ 2 √ 2
u1 = 15 m , u2 = 10 m .
3 √ 2
− √15 m 10 m
82
One can easily verify that both u1T M u2 = 0 and u2T M u1 = 0.
The general solution is written as

2 ³
X ´
q (t) = Ai sin(ωi t) + Bi cos(ωi t) ui
i=1
subject to the initial conditions

· ¸ · ¸
d 0
q (0) = , q̇ (0) = .
−d 0
Premultiplying by uiT M yields
uiT M q (0) = uiT M ui Bi , uiT M q̇ (0) = uiT M ui (Ai ωi ).

¡ ¢ ¡ ¢
Since uiT M ui = 1 from our normalization, the constants are directly solved to be
r r
m m
A1 = 0, B1 = 6 d , A2 = 0, B2 = d
15 10
Finally, the solution to the specific initial conditions becomes

· ¸ ( Ãr !· ¸ Ãr !· ¸)
x1 (t) d k 4 2k 1
q (t) = = cos t + cos t .
x2 (t) 5 3m −6 m 1
Problem 96:
The two-degree-of-freedom system shown is
subject to a harmonic force applied to the
f (t)
block of mass 2m, of the form:
f (t) = (2 sin(t) N) ı̂. 2m m

k k 2k
If the mass and stiffness of the system are
assumed to be m = 2 kg, and k = 4 N/m, x1 z
find:
a) the equations of motion; x2
b) the forced, damped equation (single-

degree-of-freedom) that describes the
motion of each mode;
Solution:
a) We define the coordinates x1 , x2 , and z, which are related as
z = x2 − x1 .
The equations of motion become

· ¸ · ¸ · ¸
2 0 2 −1 f (t)
m +k = .
0 1 −1 3 0
83
b) The corresponding eigenvalue problem can be written as
1 − 12
· ¸ µ ¶
¡ −1 ¢ k k
M K = u= β u
m −1 3 m
and the characteristic equation becomes

1 5
(1 − β) (3 − β) − = β 2 − 4 β + = 0.
2 2
This quadratic equation has the solutions

r ( √ √ √ √ )
3 2 ( 8 − 3) k ( 8 + 3) k
β =2± , −→ ω = √ , √ .
2 2m 2m
Returning to the eigenvalue equation, the eigenvectors satisfy the equation

1
ui1 − ui2 = βi ui1 −→ ui2 = 2 (1 − βi ) ui1
2
so that
s √ √ · ¸
( 8 − 3) k √ 1
ω1 = √ , u1 = ,
2m 6−2
s √ √ · ¸
( 8 + 3) k √1
ω2 = √ , u2 = .
2m −( 6 + 2)
For each eigenvector the kinetic energy inner product is

√ √
u1T M u1 = 14 − 4 6, u2T M u2 = 14 + 4 6,
Finally, the modal equation for the first mode can be written as
¡ T
u1 M u1 Q̈1 + u1T K u1 Q1 = u1T f (t),
¢ ¡ ¢
u1T f (t)
Q̈1 + ω12 Q1 = u1T M u1
,
³ √ √ ´
( 8− 3) k f (t)
Q̈1 + √
2m
Q1 = 14−4 6
√ ,
while the response of the second mode is governed by

¡ T
u2 M u2 Q̈2 + u2T K u2 Q2 = u2T f (t),
¢ ¡ ¢
u2T f (t)
Q̈2 + ω22 Q2 = u2T M u2
,
√ √
³ ´
Q̈2 + ( 8+
√
2m
3) k
Q2 = f (t)
14+4 6
√ ,
84

Mechanical Vibration Solved Examples

Uploaded by

Copyright:

Available Formats

Mechanical Vibration Solved Examples

Uploaded by

Document Information

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Mechanical Vibration Solved Examples

Uploaded by

Copyright:

Available Formats

What are the main types of vibration problems discussed in the document?

What are the main types of vibration problems discussed in the document?

What are the steps to derive the equations of motion for a multi degree-of-freedom system?

What are the steps to derive the equations of motion for a multi degree-of-freedom system?

Mechanical Vibrations: 4600-431 Example Problems

December 20, 2006

2 Frictionally Damped Systems 33

3 Forced Single Degree-of-freedom Systems 42

4 Multi Degree-of-freedom Systems 69

1 Free Vibration of Single Degree-of-freedom Systems

a) Find the equations of motion for this

c) If the block is displaced 18 cm to the

a) Find the equations of motion.

for what value of the spring stiffness k

m = 1.50 kg, ℓ = 45 cm,

Assume that in the horizontal position the

We can also relate the identified coordinates as

Using the coordinate relations we can obtain the equation of motion as

From these the damping ratio and natural frequency are

so that the general solution becomes

c) What is the force in the upper spring

find the angular displacement of the bar

θ(0) = 0, θ̇(0) = 10 rad/s.

d) for this motion, find the tension in the

Assume that the system is in static equilib-

An appropriate free-body diagram is

and the natural frequency and damping ratio are

a) find the equations of motion. Do not

b) if the system is underdamped, what is φ (m, r) θ

c) for what value of the damping constant ı̂

d) what is the static equilibrium displace- ê1

These additional coordinates are related to x as

ı̂ = cos φ ê1 + sin φ ê2 ,

The moment produced by gravity about point C is

Angular momentum balance about the contact point C yields

c) The system is critically damped when ζ = 1, which corresponds to a damping coefficient

d) The system is stationary in static equilibrium, so that x ≡ x0 = constant—both ẋ and

Solving for x0 , the equilibrium displacement is

a) Determine the deflection of each spring r2

Notice that because of these coordinate

With this, the equilibrium deflection of each spring is found to be

b) The general free response of the disk can be expressed as

θ(t) = θeq + A sin(ωn t) + B cos(ωn t),

where θeq is given above, A and B are arbitrary constants, and

The system is released with the initial conditions:

so that solving for the arbitrary constants

Therefore the solution is

The angular velocity of the disk becomes

which has amplitude

a) determine the governing equations of m

c) what value of the damping coefficient b

d) if k = 2 N/m, b = 4 N/(m/s), and m =

x(0) = 0.10 m, ẋ(0) = 0 m/s.

Problem 11: (Spring 2003)

a) find the equations of motion;

Likewise, application of linear momentum balance on the block yields:

b) A critically damped system occurs when ζ = 1. For this system:

Solving for ccr yields: r

Problem 12: (Spring 2003)