PC235W13 Assignment5 Solutions

PC235 Winter 2013
Classical Mechanics
Assignment #5 Solutions
#1 (5 points) JRT Prob. 5.26
An undamped oscillator has period 0 = 1.000s, but I now add a little
damping so that its period changes to 1 = 1.001s. What is the damping
factor ? By what factor will the amplitude of oscillation decrease after 10
cycles? Which effect of damping would be more noticeable, the change of
period or the decrease of the amplitude?
Solution
p
The damping changes the frequency to 1 = 02 2 , which we can solve
to give
s
s
12
2
(1)
= 0 1 2 = 0 1 02
0
1
= 0 1 0.998 = 0.04470 = 0.281s1 .

After a time t = 101 100 , the amplitude will have changed by a factor of
et e100 = e20/0 = e20(0.0447) = 0.060.
(2)
In other words, the amplitude will have diminished by a factor of 1/0.060

17. Clearly, the change of amplitude by a factor of 17 is more noticeable
than the change of period by a factor of 0.1%.

As the damping on an oscillator is increased, there comes a point when
the name oscillator seems barely appropriate.
(a) To illustrate this, prove that a critically damped oscillator can never
pass through the origin x = 0 more than once.
(b) Prove the same for an overdamped oscillator.
1
Solution
(a) If the oscillator is critically damped, = 0 , and x(t) = et (C1 +C2 t).
This vanishes if and only if t = C1 /C2 . Since C1 and C2 are constants,
there is at most one such value for t (it may never vanish; for example
if the motion starts at t = 0 and C1 /C2 < 0).
(b) If the oscillator is overdamped, > 0 and
x(t) = et (C1 et + C2 et ) = e()t (C1 + C2 e2t ),
(3)
p
where = 2 02 . This vanishes if and only if e2t = C1 /C2 , and
since e2t is a monotonic function (always decreasing), this happens
at most once.

Consider an underdamped oscillator (such as a mass on the end of a
spring) that is released from rest at position x0 at time t = 0.
(a) Find the position x(t) at later times in the form
x(t) = et [B1 cos(1 t) + B2 sin(1 t)] .
(4)
That is, find B1 and B2 in terms of x0 .

(b) Now show that if you let approach the critical value 0 , your solution
automatically yields the critical solution.
(c) Using appropriate graphing software, plot the solution for 0 t 20,
with x0 = 1, 0 = 1, and = 0, 0.02, 0.1, 0.3, and 1.
Solution
(a) If we write x(t) in the suggested form, it is easy to see that x(0) = B1
and x(0)
= 1 B2 B1 . Therefore, B1 = x0 and B2 = x0 /1 , so that

t
sin 1 t .
(5)
cos 1 t +
x(t) = x0 e
1
p
(b) When 0 , 1 = 02 2 approaches zero. The cos 1 t term thus
approaches one. To deal with the (/1 ) sin 1 t term, we then make
use of the limit
2
sin 1 t
t sin 1 t
= tsinc(0) = t.
(6)
= lim
0
1
1 t
1
Therefore, the solution of part (a) becomes
x(t) = x0 et (1 + t),
(7)
which is precisely the critically-damped solution satisfying the initial
conditions (see eq. 6.34 in the class notes.)
lim
1 0
1
0
10
t
=0.1
15
1
0
20
1
0
(c) The requested plots are shown below.
10
t
=1
15
20
10
t
15
20
1
0
10
t
=0.3
15
20
10
t
15
20
1
0
1
0
Fig. 1: Plots for question #3

Another interpretation of the Q of a resonance comes from the following: Consider the motion of a driven damped oscillator after any transients
have died out, and suppose that it is being driven close to resonance, so you
can set = 0 .
(a) Show that the oscillators total energy (kinetic plus potential) is E =
1
m 2 A2 .
2
(b) Show that the energy Edis dissipated during one cycle by the damping
force Fdmp is 2mA2 .
3
(c) Hence, show that Q is 2 times the ratio of Edis .

Solution
(a) Since x = A cos(t ), the total energy is
1
1
1
1
E = mx 2 + kx2 = m 2 A2 cos2 (t ) + kA2 sin2 (t ). (8)
2
2
2
2
Because 0 , we can replace k = m02 with m 2 , and then, with
cos2 + sin2 = 1, we get the required E = 21 m 2 A2 .
(b) The rate at which the damping force dissipates energy is Fdmp v = bv 2 =
2mv 2 . Therefore, the energy dissipated in one period is
Z
Z
2
2
2 2
bv = 2mv dt = 2m A
Edis =
sin2 (t ) dt.
(9)
0
The integral evaluates to /, and we are left with Edis = 2mA2 .

(c) Combining the results of parts (a) and (b), we find that
1
m 2 A2
0
Q
E
2
=
=
,
=
2
Edis
2mA
4
2
(10)
where we have again used = 0 . Therefore, the ratio of total energy

to energy lost per cycle is Q/2.

Repeat problem 6.1 but find the length L of a path on a cylinder of
radius R, using cylindrical polar coordinates (, , z). Assume that the path
is specified in the form = (z).
Solution
Consider an infinitesimal section of path on the cylinder, in which increases
by d and z increases by dz. This carries us a distance Rd around the cylinder and a distance dz up it. The distance ds along the path is therefore
p
p
(11)
ds = (Rd)2 + (dz)2 = R2 (z)2 + 1dz.
4
Therefore, the total path length is

Z
Z 2
ds =
L=
1
2
1
R2 (z)2 + 1dz.
(12)

A ray of light travels from point P1 in a medium of refractive index n1
to a point P2 in a medium of index n2 , by way of the point Q on the plane
interface between the two media. Show that Fermats principle implies that,
on the actual path followed, Q lies in the same vertical plane as P1 and P2
and obeys Snells law, that n1 sin 1 = n2 sin 2 .
Solution
The lengths of the paths P1 Q and QP2 are
q
q
2
2
2
P1 Q = x + h1 + z and QP2 = (x2 x)2 + h22 + z 2 .
(13)
The time for light to traverse each path is the path length divided by the
local speed of light, v = c/n. Thus, the total time is

q
q
1
2
2
2
2
2
2
(14)
n1 x + h1 + z + n2 (x2 x) + h2 + z .
t=
c
To find where this is minimum, we must set
first of these conditions gives
t
1
=
z
c
n1 z
t
z
and
n2 z
t
x
equal to zero. The
p
+p
(x2 x)2 + h22 + z 2
x2 + h21 + z 2
= 0,
(15)
which is satisfied only when z = 0. That is, Fermats principle requires that
Q lie in the plane containing P1 and P2 and containing the interface. The
second condition gives
!
t
n1 x
1
1
n2 (x2 x)
p
p
=
(n1 sin 1 n2 sin 2 ) = 0,
=
+
x
c
c
(x2 x)2 + h22 + z 2
x2 + h21 + z 2
(16)
5
which is satisfied only if n1 sin 1 = n2 sin 2 ; this is Snells law.

Rx p
Show that the path y = y(x) for which the integral x12 x 1 y 2 dx
is stationary is an arcsinh function.
Solution
p
The integrand is f (y, y , x) = x 1 y 2 . Since this is independent of y,
f /y = 0, p
and the E-L equation implies simply that f /y is a
constant;
that is, xy / 1 y 2 = k. This can be rearranged to give y = k/ x2 + k 2 ,
which integrates to give y = k sinh1 (x/k) + C, where C is a constant of
integration (to perform the integral, make the substitution x/k = sinh u).
A more pleasing representation of this path is found by rearranging to find
x = k sinh
yC
.
k
(17)

An aircraft whose airspeed is v0 has to fly from town O (at the origin)
to town P , which is a distance D due east. There is a steady gentle wind
shear, such that vwind = V yx, where x and y are measured east and north,
respectively. Find the path y = y(x), which the plane should follow to
minimize its flight time, as follows:
(a) Find the planes ground speed in terms of v0 , V, (the angle by which
the plane heads to the north of east), and the planes position.
RD
(b) Write down the time of flight as an integral of the form 0 f dx. Show
that if we assume that y and both remain small (as is certainly the
case if the wind speed is a small fraction of the planes airspeed), then
the integrand f takes the approximate form f = (1 + y 2 /2)/(1 + ky)
(times an uninteresting constant), where k = V /v0 .
(c) Write down the Euler-Lagrange equation that determines the best path.
To solve it, make the intelligent guess that y(x) = x(D x), which
passes through the two towns. Show that it satisfies the E-L equation,
6

provided that = ( 4 + 2k 2 D2 2)/(kD2 ). How far north does this
path take the plane, if D = 2000 miles, v0 = 500 mph, and the wind
shear is V = 0.5 mph/mile? How much time does the plane save by
following this path?
Solution
(a) If the plane is at position (x, y) and is heading at an angle north
of east, its velocity relative to the air is (v0 cos , v0 sin ), while the
winds velocity is (V y, 0). The planes velocity relative to the ground
is the sum of these two vectors, so its ground speed is
p
v =
(v0 cos + V y)2 + (v0 sin )2
(18)
q
v02 + 2v0 V y cos + V 2 y 2
(19)
=
v0 + V y = v0 (1 + ky),
(20)
where we have made the approximation cos 1 (a zeroth-order

approximation), and where k = V /v0 .
(b) The time of flight is
p
Z D
Z
1 + 12 y 2
1 + y 2 dx
1 D
t=
dx =
f dx,
v
v0 0
0
0
0 v0 (1 + ky)
(21)
where we used the binomial approximation in the numerator (assuming
y is small) and where f = (1 + y 2 /2)/(1 + ky).
Z
ds
=
v
(c) A bit of algebra shows that the E-L equation is

1
y (1 + ky) ky 2 + k = 0.
2
(22)
We dont possess the skills necessary to solve this nonlinear 2nd-order

differential equation, so we will substitute the proposed guess, y =
x(D x), which results in
kD2 2 + 4 2k = 0,
(23)
which has the solution = ( 4 + 2k 2 D2 2)/(kD2 ) (the other solution

is negative, which makes no physical sense since the wind direction is
7
always opposing the desired direction of travel for y < 0).

The maximum displacement to the north is ymax = D2 /4 = 366 miles.
If we solve the integral in part (b), we find for the time of flight t = 3.556
hours. The flight time along the direct path is 4 hours (theres no wind
on the x-axis), so the time saved is 0.444 hours, or 27 minutes.
#9 (10 Bonus points) JRT Prob. 5.53

An oscillator is driven by the periodic force shown in the figure below,
which has period = 2.
(a) Find the long-term motion x(t), assuming the following parameters:
natural period 0 = 2 (that is, 0 = ), damping parameter = 0.1,
and maximum drive strength fmax = 1. Find the coefficients in the
Fourier series for x(t) and plot the sum of the first four terms in the
series for 0 t 6.
(b) Repeat, except with natural period equal to 3.
Fig. 2: Periodic driving force of question #9
Solution
Because the function f (t) is even, all of the Fourier coefficients bm are zero.
We need only to find the coefficients am . From eq. (6.68) from the class
notes,
Z
1 /2
fmax
a0 =
(24)
f (t)dt =
/2
2
while for m 1 (and noticing that the portion of the integrand to the left
and to the right of t = 0 contribute equally to the integral),
Z
Z 1
4 /2
cos(mt)f (t) dt = 2fmax
cos(mt)(1 t) dt.
(25)
am =
0
0
8
This integral is evaluated using integration by parts, yielding

0
m even
2fmax
1
.
am =
[cos(mt]0 =
4fmax
2
m odd
(m)
(m)2
(26)
From eqns. (6.71-6.72) of the class notes, we know that the motion of the
oscillator is described by
x(t) =
X
n=0
An cos(nt n )
(27)
with
An = p
fn
(02 n2 2 )2 + 4 2 n2 2
, n = tan
2n
02 n2 2
(28)
(for the case n = 0 this just becomes A0 = 21 2 and 0 = 0.) From the figure
0
illustrating the driving force, we see that the driving frequency = 2/ = ,
where is the period of the driving force.
(a) For 0 = 2, we have 0 = . For = 0.1 and fmax = 1, this means
that the first four coefficients An (n = 0, 1, 2, 3) are 0.0507, 0.6450, 0,
and 0.0006. Note that the term A1 dominates the sum, because this
term is on resonance with the natural frequency of the oscillator. The
oscillations are plotted below.
(b) For 0 = 3, we have 0 = 2/3, and the first four coefficients An
(n = 0, 1, 2, 3) are 0.1140, 0.0734, 0, and 0.0005. In this case, note
that the constant term A0 is the largest, meaning that the oscillations
are significantly biased toward positive x. The oscillations are plotted
below.
0.8
0.6
0.4
x(t)
0.2
0
0.2
0.4
0.6
0
3
time
Fig. 3: Question #11 - motion of the oscillator. Blue: 0 = 2. Green: 0 = 3.
10

PC235W13 Assignment5 Solutions

Uploaded by

Copyright:

Available Formats

PC235W13 Assignment5 Solutions

Uploaded by

Document Information

Original Description:

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

PC235W13 Assignment5 Solutions

Uploaded by

Copyright:

Available Formats

PC235 Winter 2013

= 0 1 0.998 = 0.04470 = 0.281s1 .

In other words, the amplitude will have diminished by a factor of 1/0.060

#2 (5 points) JRT Prob. 5.27

#3 (5 points) JRT Prob. 5.32

That is, find B1 and B2 in terms of x0 .

= 1 B2 B1 . Therefore, B1 = x0 and B2 = x0 /1 , so that

(c) The requested plots are shown below.

Fig. 1: Plots for question #3

#4 (5 points) JRT Prob. 5.44

(c) Hence, show that Q is 2 times the ratio of Edis .

The integral evaluates to /, and we are left with Edis = 2mA2 .

where we have again used = 0 . Therefore, the ratio of total energy

#5 (5 points) JRT Prob. 6.2

Therefore, the total path length is

#6 (10 points) JRT Prob. 6.4

equal to zero. The

which is satisfied only if n1 sin 1 = n2 sin 2 ; this is Snells law.

#7 (5 points) JRT Prob. 6.12

#8 (10 points) JRT Prob. 6.23

where we have made the approximation cos 1 (a zeroth-order

(c) A bit of algebra shows that the E-L equation is

We dont possess the skills necessary to solve this nonlinear 2nd-order

which has the solution = ( 4 + 2k 2 D2 2)/(kD2 ) (the other solution

always opposing the desired direction of travel for y < 0).

#9 (10 Bonus points) JRT Prob. 5.53

Fig. 2: Periodic driving force of question #9

This integral is evaluated using integration by parts, yielding

Fig. 3: Question #11 - motion of the oscillator. Blue: 0 = 2. Green: 0 = 3.

You might also like