More Problems

Problem 1:

A triple pendulum consists of masses am, m and m attached to a single light string at distances a, 2a, and 3a respectively from its points of suspension.  Consider only motion in a plane.
(a)  Determine the value of a such that one of the normal frequencies of the system will equal the frequency of a simple pendulum of length a/2 and mass m.  You may assume the displacements of the masses from equilibrium are small.
(b)  Find the mode corresponding to this frequency and sketch it.

Solution:

Concepts, principles, relations that apply to the problem:
Small oscillations, coupled oscillations, normal modes

Solutions of the form qj = Re(Ajeiwt) can be found.  We can find the w2 from det(kij - w2Tij) = 0.
Why do they apply?
We are asked to find the normal mode frequencies for small displacements.
How do they apply?

Let f1, f2, and f3 be the generalized coordinates.  Then
x1 = a sin
f1, x2 = a(sinf1 + sinf2), x3 = a(sinf1 + sinf2 + sinf3),
y1 = a
cosf1, y2 = a(cosf1 + cosf2), x3 = a(cosf1 + cosf2 + cosf3).
For small displacements we have sin
f = f, cosf = 1 - f2/2.
The kinetic energy of the system is
.
Let us only keep terms up to second order in
f and df/dt.  Then
.
.
The potential energy of the system is
U = -
amgy1 - mgy2 - mgy3.
.
We have:
.
T11 = ma2(
a + 2), T22 = 2ma2, T33 = ma2, T12 = T21 = 2ma2, T13 = T31 = T23 = T32 = ma2.
k11 = mga(
a + 2), k22 = 2mga, k33 = mga, kij(i ¹ j) = 0.

.

.
Details of the calculations:
For a simple pendulum with length a/2 we have w = (2g/a)1/2.
  We are looking for a mode with w2 = 2g/a, or l = 2.
.
The mode corresponding to this frequency is found from

.

A3 = -2A1, A2 = 0.

Problem 2:

A uniform thin rod of length 3l/2 and mass m is supported at one end A by a weightless string of length l under the gravitational force near the earth surface.

(a)  Calculate the normal modes and normal frequencies of small oscillations for motion in a vertical plane.
(b)  Now the end point A is slowly displaces by a small amount d, (without gaining any kinetic energy.)  The system is then released from rest and allowed to move freely.  What is the subsequent motion?

Solution:

Concepts, principles, relations that apply to the problem:
Small oscillations:  

Solutions of the form qj = Re(Ajeiwt) can be found.
For a particular frequency wa we solve



to find the Aja.  The most general solution for each coordinate qj is a sum of simple harmonic oscillations in all of the frequencies wa.

.
Why do they apply?
We are asked to find the normal modes for small displacements.

How do they apply?

(a)  The position of the CM of the rod is X = x1 + x2/2, Y = y1 + y2/2.
For the rod we have:
.
T is the kinetic energy of the motion of the CM and the motion about the CM.
U = -mg(y1 +y2/2) + constant.
For small oscillations x1 = lsinq1 @ lq1, x2 = (3/2)lsinq @ (3/2)lq.
y1 = lcosq1 @ l(1-q12/2), y2 = (3/2)lcosq @ (3/2)l(1-q2/2).
To second order in the small quantities we have:
T = (1/2)ml2((dq1/dt)2+(3/4)2(dq/dt)2+(3/2)(dq1/dt)(dq/dt)+(1/2)(3/16)ml2(dq/dt)2)
= (1/2)ml2((dq1/dt)2+(3/4)(dq/dt)2+(3/2)(dq1/dt)(dq/dt))
= (1/2)STijqiqj.
U = mgl[(q12/2) +(3/4)l(q2/2)] = (1/2)Skijqiqj.
q1 = q1, q2 = q, T11 = ml2, T22 = (3/4)ml2, T12 = T21 = (3/4)ml2.
k11 = mgl, k22 = (3/4)mgl, k12 = k21 = 0.
Details of the calculation:
.
We have a solution if
.
(3/4)(g-w2l)2 -(9/16)w4l2 = 0,  w4 - 8(g/l)w2 + 4g2/l = 0.
w2 = [4 ± (12)1/2]g/l.
mode 1: 
w12 = [4 + (12)1/2]g/l.
(g - [4
+ (12)1/2]g)A1 - (3/4)[4 + (12)1/2]gA2 = 0.  A1 = -0.866 A2.
mode 2:  w22 = [4 - (12)1/2]g/l.
(g - [4 - (12)1/2]g)A1 -
(3/4)[4 - (12)1/2]gA2 = 0.  A1 = 0.866 A2.

(b)  The most general solution is
.
q1 = Re(c1exp(iw1t) + c2exp(iw2t)),  q = (1/0.866)Re(-c1exp(iw1t) + c2exp(iw2t)).
Initial conditions:
q1(0) = d/l,  q(0) = 0,  dq1/dt = dq/dt = 0 at t = 0.
Re(c1 + c2) = d/l,  Re(c1 - c2) = 0,  Re(c1) = Re(c2) = d/(2l).
Re(iw1c1 + iwc2) = 0,  Re(-iw1c1 + iwc2) = 0.
Im(w1c1 + wc2) = 0,  Im(w1c1 - wc2) = 0,  Im(c1) = Im(c2) = 0.
q1(t) = (d/(2l))[cosw1t + cosw2t],  q(t) = (d/(2*0.866*l))[cosw1t - cosw2t].

Problem 3:

Two masses a and b are on a horizontal surface.  Mass b has a spring connected to it and is at rest.  Mass a has an initial velocity v0 along the x-axis and strikes the spring of constant k, compressing it and thus starting mass b in motion along the x-axis.

(a)  Find the maximum force exerted.
(b)  Find the resulting motion of mass a while in contact with the spring.

Solution:

Concepts, principles, relations that apply to the problem:
Small oscillations, coupled oscillations, normal modes

Solutions of the form qj = Re(Ajeiwt) can be found.  We can find the w2 from det(kij - w2Tij) = 0.
Why do they apply?
During the time interval that both masses are in contact with the spring, they can be treated as coupled harmonic oscillators.
How do they apply?
(a)  Assume that at t = 0 mass a makes contact with the spring.  At t = 0 both masses are at their equilibrium position.  Let xa(t) and xb(t) denote the their displacements from their equilibrium positions, respectively.
In the collision energy and momentum are conserved.
Ti = Tf + U, Pi = Pf.  The maximum force is exerted when the spring has maximum compression.  At that instant both masses move with the same velocity.
(1/2)mav02 = (1/2)(ma + mb)v2 + (1/2)k(xb - xa)min2.
mav0 = (ma + mb)v.
Combining the two equations we have (xb - xa)min = (mamb/k)1/2 v0/(ma + mb)1/2.
The maximum force is k(xb - xa)min = (kmamb)1/2 v0/(ma + mb)1/2.

(b)  T = (1/2)ma(dxa/dt)2 + (1/2)mb(dxb/dt)2
U = (1/2)k(xb - xa)2 = (1/2)k(xb2 + xa2 - xaxb - xbxa).
T11 = ma, T22 = mb, k11 = k22 = k, k12 = k21 = -k.
det(kij - w2Tij) = 0  --> w1 = 0, w22 = k(ma + mb)/(mamb).
For w1 = 0 we have A1 = A2.  For w22 = k(ma + mb)/(mamb) we have A2 = -(ma/mb)A1.
The most general solutions are:
xa = A1t + B1cos(wt + f),  xb = A1t - (ma/mb)B1cos(wt + f),
dxa/dt = A1 - wB1sin(wt + f),  dxb/dt = A1 + w(ma/ma)B1sin(wt + f).
Details of the calculations:
The initial conditions at t = 0 are xa = xb = 0, dxa/dt = v0, dxb/dt = 0.
This implies
B1cos(f) = 0, f = p/2.
v0
= A1 - wB1,  A1 + w(ma/mb)B1 = 0.
A1 = v0ma/(ma + mb),  B1 = -(1/w)v0mb/(ma + mb).
Therefore
xa = (v0/w)mb/(ma + mb)sinwt + mav0/(ma + mb)t,
dxa/dt =
v0mb/(ma + mb)coswt + mav0/(ma + mb),
as long as mass a is in contact with the spring.
[At t = 0 we have dxa(lab) /dt =
v0mb/(ma + mb) + mav0/(ma + mb) = v0.]

A different approach:

Solution:

Concepts, principles, relations that apply to the problem:
Simple harmonic motion, frame transformations, two interacting particles
Why do they apply?
In the center of mass frame we have an elastic collision between two particles with an interaction energy U(x) = (1/2)k(x-x0)2, where x-x0 is the distance between the particles, as long as particle a is in contact with the spring . 
How do they apply?
(a)  In the center of mass frame we have an elastic collision.  Assume mass a approaches the CM with with speed va from the left, and mass b approaches the CM with with speed vb from the right.  The CM is at rest and mass a recedes from the CM with with speed va to the right, while mass b recedes from the CM with with speed vb to the left after the collision.
We have mava = mbvb in the CM frame before the collision.  To find the relationship between the speeds in the CM and the lab frame we use
ma(v0 - vb) = mbvb.  vb = mav0/(ma + mb)),  va = mbv0/(ma + mb).
The maximum force is exerted when all the kinetic energy is converted into potential energy.
Fmax = k(xmin-x0).
k(xmin-x0)2 = mava2 + mbvb2 =  ma(mbv0/(ma + mb))2 + mb(mav0/(ma + mb))2 = mambv02/(ma + mb).
k(xmin-x0) = (kmamb)1/2 v0/(ma + mb)1/2.
Details of the calculations:
(b)  Assume that mass a makes contact with the spring at t = 0, at the origin of the CM and the lab frame.  (We assume that at t = 0 the origin of the two frames coincide.)  Assume that at t = 0 the distance between mass a and mass b is x0.
In the CM frame when mass a moves a distance xa to the right, mass b moves a distance xb = xa(ma/mb)  to the left.  The distance between the masses becomes x = x0 - (xa + xb) = x0 - xa(ma + mb)/mb = x0 - x'.
For t > 0 we have U = (1/2)kx'2 = (1/2)[k(ma + mb)2/mb2]xa2 until the time when contact is lost again.
The kinetic energy of the particles is
T = (1/2)ma(dxa/dt)2 + (1/2)mb(dxb/dt)2 = (1/2)(ma + ma2/mb)(dxa/dt)2.
L = T - U = (1/2)(ma(mb + ma)/mb)(dxa/dt)2 - 1/2)[k(ma + mb)2/mb2]xa2.
This is a Lagrangian for simple harmonic motion.
d2xa/dt2 = -k
(mb + ma)/(mamb)xa = -(k/m)xa, with m = mamb/(mb + ma).
We therefore have
xa = Asinwt, w = (k/m)1/2,
until t = p/w (1/2 period of the motion).
dxa/dt = wAcoswt
.
At t = 0 we have wA =
va = mbv0/(ma + mb),  A = mbv0/[w(ma + mb)]
In the lab frame we have xa(lab) = xa + vbt.

The motion of mass a is described by
xa(lab) = (v0/w)mb/(ma + mb)sinwt + mav0/(ma + mb)t,
dxa(lab) /dt =
v0mb/(ma + mb)coswt + mav0/(ma + mb),
as long as mass a is in contact with the spring.
[At t = 0 we have dxa(lab) /dt =
v0mb/(ma + mb) + mav0/(ma + mb) = v0.]