Download Curvilinear Systems, Central Forces

Physics 430: Lecture 9
Curvilinear Systems,
Central Forces
Dale E. Gary
NJIT Physics Department
4.7 Curvilinear One-Dimensional
Systems
In addition to an object moving in one dimension in space, such as the
horizontal or vertical motions we have been dealing with, we can also treat
constrained motion along a curved path in terms of 1 dimension. An
example is a bead threaded onto a curved, rigid wire on which it can slide
without friction.
 Such a curve can be specified in terms of its three-dimensional (x, y, z)
coordinates, which might require three separate equations, but we could
instead use a single, parametric equation, i.e. specified the bead’s
location in terms of a single parameter, say its distance s along the wire
from some origin O. With such a choice, problems can be solved as if they
were one-dimensional.
2
 The speed (not velocity) is written s so the kinetic energy is T  12 ms
2
compared with 12 mx for a straight track. However, the force is a little more
complicated, since the wire exerts a force of constraint.
But notice that the constraint force does no work, since
O
it is not in the direction of motion, but only perpendicular. We only care about the tangential force Ftang, which is Ftang  ms

September 29, 2009
Potential Energy of Tangential
Forces





If all of the tangential forces are due to conservative forces, then we can
define a potential energy U(s) such that Ftang = -dU/ds. In this case, the
total energy E = T + U is constant, and we can use all of the ideas we have
discussed for the one-dimensional case.
In particular, points along s where U(s) is a minimum are stable equilibrium
points, and points where U(s) is a maximum are unstable.
There are many problems that at first appear to be more complicated than
a bead on a wire, but they are nevertheless one-dimensional and so can be
solved using these ideas.
A system composed of many bodies all in motion may seem multidimensional, but if they are connected by struts or strings so that the
position of the system can be described by a single parameter, then they
can be considered one-dimensional.
Let’s take a look at a couple of examples.
September 29, 2009
Example 4.7 Stability of a Cube
Balanced on a Cylinder

Statement of Problem:


A hard rubber cylinder of radius r is held fixed with its axis horizontal, and a wooden
cube of mass m and side 2b is balanced on top of the cylinder, with its center
vertically above the cylinder’s axis and four of its sides parallel to the axis. The cube
cannot slip on the rubber of the cylinder, but can of course rock from side to side, as
shown in the figure. By examining the cube’s potential energy, find out if the
equilibrium with the cube centered above the cylinder is stable or unstable.
Solution:


Note that the system is one-dimensional, since its position can
be parameterized by a single parameter, q. The constraint
forces are the normal and frictional forces of the cylinder on
the cube (because the cube does not slide), which do no work
(because they are perpendicular to the motion). The only
other force is gravity, which is a conservative force, so we can
define the potential energy as U = mgh, where h is the height
of the CM of the cube above the CM of the cylinder.
If there is any trick, it is to recognize that BC = rq.
C
B
b
q r
O
September 29, 2009
h
Example 4.7 Stability of a Cube
Balanced on a Cylinder

Solution, cont’d:



The angle BC makes from the horizontal is also q.
The length of h is BC sin q + (r+b) cos q = rq sin q + (r+b) cos q.
This let’s us write the potential energy of the system as
B
U  mgh  mg[( r  b) cos q  rq sin q ]



q
C
b
Remember that the equilibrium position will be a location
where the slope of U is zero (dU/dq = 0), and it will be stable if
the curvature is upward (d2U/dq 2 > 0).
q r
O
The value of q for which the derivative is zero is found by
dU
 mg[-(r  b) sin q  r sin q  rq cos q ]
dq
 mg[rq cos q - b sin q ]  0
Cube must be smaller
than cylinder
This is zero when q = 0, which confirms the obvious—the cube can be balanced
directly on top of the cylinder. But is the balance precarious? To find out:
d 2U
 mg[r cos q - rq sin q - b cos q ]  mg[( r - b) cos q - rq sin q ]
2
dq

At the equilibrium position q = 0, this is just mg(r - b), which is positive for r > b.
September 29, 2009
h
Another Example—Atwood Machine





The Atwood Machine provides another example of a system
with several moving parts that can be cast into a oneFT x
dimensional framework.
Here the pulley is massless (not essential to the problem)
m1
FT
and the cable connecting the two masses cannot stretch.
The position of the whole system can be determined by the m
m1g
2
position x of the mass below the center of the pulley.
m2g
Consider the energies of the two masses after a
displacement:
DT1  DU1  W1ten and DT2  DU 2  W2ten
where DT1 is the change in kinetic energy after the displacement, DU1 is the
change in potential energy due to the displacement, and W1ten is the work done
during the displacement against the tension force FT.
The tension force is exactly the same on the two masses, so the work done by
tension is equal and opposite, i.e. W1ten  -W2ten .
September 29, 2009
Another Example—Atwood Machine


Adding the two energies, then, the work done by tension
drops out and the total change in energy is zero:
D(T1  U1  T2  U 2 )  0.
FT x
In other words, the total mechanical energy
E  T1  U1  T2  U 2
FT
is conserved.
 It happens that many systems that contain several particles
that are constrained in some way (by strings, struts, or a
track) can be treated in this way.
m1
m2
m2g
m1g
The constraint forces are crucially important for determining how the system
moves, but they do no work on the system, so when using energy to solve
problems we can simply ignore the constraint forces.
N
 For each particle a, with potential energy Ua, the total energy E   (Ta  U a )
a 1
is constant.

September 29, 2009
Problem 4.36

Statement of the problem, part (a):


A metal ball (mass m) with a hole through it is threaded on a frictionless vertical
rod. A massless string (length l) attached to the ball runs over a massless,
frictionless pulley and supports a block of mass M, as show in the figure. The
postions of the two masses can be specified by the one angle q. (a) Write down
the potential energy U(q). (The PE is given easily in terms of the heights shown as
h and H. Eliminate these two variables in favor of q and the constants b and l.
Assume that the pulley and ball have negligible size.)
Solution to part (a):
 The PE is just U  - mgh - MgH . We need to write h
b
and H in terms of l, b, and q, but clearly h = (l - H) cos q,
while b = (l - H) sin q. Thus,
h  b / tan q , and H  (l - b / sin q ).

Plugging into the expression for U, we have
U (q )  -mgb / tan q  Mgb / sin q - Mgl.

H
We are free to drop the uninteresting constant Mgl,
which amounts to choosing a different zero for U(q), so
U (q ) 
l-H
M
gb
( M - m cos q ).
sin q
September 29, 2009
h
q
m
Problem 4.36

Statement of the problem, part (b):


(b) By differentiating U(q) find whether the system has an equilibrium position, and
for what values of m and M equilibrium can occur. Discuss the stability of any
equilibrium positions.
Solution to part (b):




Taking the derivative and setting it to zero, we have
dU
gb
gb
 - cos q
( M - m cos q ) 
(m sin q )
2
dq
sin q
sin q
gb cos q

(m cos q - M )  mgb  0
sin 2 q
Multiplying through by sin2q (note q cannot be zero), and
dividing by gb, we have m - M cosq  0.
m
So equilibria occur at
cos q 
M
,
b
H
l-H
h
q
m
M
m
which can only happen when m < M. In that case, there is an equilibrium at q o  arccos
M
Normally, to see if this is a stable equilibrium, take the second derivative,
evaluation at q = qo, and check if it is positive (stable) or negative (unstable). In
this case, that leads to a very nasty expression. Instead, imagine pulling the mass
m downward (q < qo), letting go, and asking what direction the net force acts. It is
clearly upward, tending to restore the equilibrium. Likewise, moving the mass
upward leads to a net force downward. Hence the equilibrium is stable.
September 29, 2009
4.8 Central Forces
We can treat central force problems (such as gravitation and electric force
problems) using some of the simplicity provided by one-dimensional problems,
because the symmetry of the problems mean that the dependence of the
motion can be cast into one involving only the radial component.
 Any central force can be written generally as


F(r)  f (r)rˆ
where f(r) gives the magnitude of the force as a function of position r.
Examples are the Coulomb force and the gravitational force
kQq
GMm
F(r )  2 rˆ ,
F(r )  - 2 rˆ.
r
r
These forces have two properties not necessarily shared by all central forces:
1) they are conservative, and 2) they are spherically symmetric (or
rotationally invariant). That is, the magnitude of the force does not depend
on the direction of r, but depends only on the magnitude of r, i.e. f (r )  f (r ).
 What is remarkable is that, for central forces, these two conditions always go
conservati ve  rotational ly invariant
together. That is,

September 29, 2009
Spherical Polar Coordinates
One way to prove the foregoing statement conservati ve  rotational ly invariant ,
we can use spherical polar coordinates. This gives us a good opportunity to
introduce them, which we do now.
 We already introduced polar coordinates and noted how the r and f unit
vectors work. We also extended that to a z component in cylindrical
coordinates. In spherical polar coordinates, in contrast, we introduce another
angular variable q, called the polar angle (or zenith angle when talking about
sky coordinates from some local position on the Earth). The relationships are
noted in the figure below.
z
r cosq
 By inspection of this figure, you should be able to see
r̂
φ̂
that
x  r sin q cos f , y  r sin q sin f , z  r cos q
θ̂
q
r
r sin q cos f
 Note that the unit vectors all point in the directions
y
f
of increase in the vectors. The direction of the r̂ unit
r sin q
vector is obvious. The θ̂ and φ̂ unit vectors are tangent xr sin q sin f
to the lines of “latitude” and “longitude” through the point.

September 29, 2009
Gradient in Spherical Polar Coordinates



We would like to express the gradient of a function f in spherical polar
coordinates, but before we do, we note that the dot product in these
coordinates behaves just as you would expect. If
z
ˆ
ˆ
ˆ
a  ar , aq , af   ar r  aq θ  af φ,
b  br , bq , bf   br rˆ  bq θˆ  bf φˆ . Circumference 2pr sin q
Then the dot product is
q
a b  ar br  aq bq  af bf . Circumference 2pr
f
Recall that the gradient of a function f in cartesian
coordinates was
f  xˆ
f
f
f
 yˆ
 zˆ .
x
y
z
Circumference 0
r
r̂ φ̂
θ̂
x

Recall also that the change in f due to a small displacement dr is
df  f  dr.

All we need to do, then, is express the displacement dr in spherical polar
coordinates.
dr  dr rˆ  rdq θˆ  r sin q df φˆ
dV=dr rdq r sin q df
= r2 dr sin q dq df
note, not just dq
note, not just rdf
September 29, 2009
y
Gradient in Spherical Polar Coordinates

We can rewrite df  f  dr completely generally as
df  f r dr  f q r dq  f f r sin q df ,
but df is also
f
f
f
df  dr 
dq 
df .
r
q
f

Equating these, term by term, we must have
f r  f ,
f q  1 f ,
r
r q
f f 
1 f
.
r sin q f

Writing this more compactly, we have

f ˆ 1 f
1 f
θ
 φˆ
.
r
r q
r sin q f
You can see that this is somewhat more complicated than in cartesian
coordinates. Likewise the other operators of vector calculus are also more
complicated, like the divergence and curl. The expressions for these are shown
on the inside back cover of the text.
f  rˆ
September 29, 2009
Conservative and Spherically
Symmetric Central Forces
Now we are ready to show that central forces that are conservative are also
spherically symmetric, and those that are spherically symmetric are also
conservative.
 Let’s first assume that the force F(r) is conservative and prove that means it
must be spherically symmetric. Because it is conservative, we can express it as
the gradient of a scalar potential energy.
U ˆ 1 U
1 U
F(r )  -U  -rˆ
-θ
- φˆ
.
r
r q
r sin q f

Because the force is central, only its radial component can be nonzero, so the
last two terms must be zero, i.e. U / q  U / f  0. That means the potential
energy is spherically symmetric, hence the force is too:
U
F(r )  -rˆ
.
r
 We will show the converse in a moment. The point, though, is that problems
involving central forces that are spherically symmetric are nearly as simple as
one-dimensional problems, as we will see in Chapter 8.

September 29, 2009
Problem 4.43

Statement of the Problem:


In Section 4.8 I claimed that a force F(r) that is central and spherically symmetric is
automatically conservative. Here are two ways to prove it: (a) Since F(r) is central
and spherically symmetric, it must have the form F(r)  f (r )rˆ.
Using cartesian coordinates, show that this implies that   F  0.
Solution to part (a):


The first thing we need to do is express r̂ in terms of x,y,z coordinates. To do this,
r ( x, y , z ) x
y
z
recall that
rˆ  
 xˆ  yˆ  zˆ .
r
r
r
r
r
f (r )
Thus, F(r ) 
( x xˆ  y yˆ  z zˆ )  g (r )( x xˆ  y yˆ  z zˆ ). Taking the curl, we have
r
xˆ

F 
x
g (r ) x

yˆ

y
g (r ) y
zˆ
 g (r )

g (r )   g (r )
g (r )   g (r )
g (r ) 
xˆ   x
zˆ
  z
-y
-z
-x
yˆ   y
z
z   z
x  
x
y 
 y
g (r ) z
g (r ) g (r ) r and
To deal with a term like g ( r ) , we use the chain rule

,
y
write r 
y
r y
x2  y2  z 2 .
September 29, 2009
Problem 4.43-cont’d

Solution to part (a), cont’d:

r
y
 , and likewise for the other
You can easily show that the partial derivative

y
r
partial derivatives of r. Thus, the first term is
  F x   z g (r ) - y g (r ) xˆ   z y g (r ) - y z g (r ) xˆ  0




z 
 r
r
r r 
Obviously, the same will be true for the other components, so indeed,   F  0.
We saw last time that this is enough to ensure that the work done by F is path
independent, hence, the force is conservative, which is what we needed to show.
Statement of the problem, part (b):


y
Even quicker, using the expression given inside the back cover for   F in
spherical polar coordinates, show that   F  0.
Solution to part (b):

Because the force is central, it has only a radial component fr, hence the expression
 

sin q ff  -  fq   θˆ  1  f r - 1  rff 
f   r sin q f
r r
 q

1 
  ˆ 1

1 
 φˆ  r fq  f r   θ
f r - φˆ
fr  0
r  r
f 
r sin q f
r f
  F  rˆ
1
r sin q
September 29, 2009