A Hostetler Handbook
v 0.82
Contents
Preface v
1 Newton’s Laws1 1.1 Many Particles ...... 2 1.2 Cartesian Coordinates...... 3 1.3 Polar Coordinates ...... 4 1.4 Cylindrical Coordinates ...... 7 1.5 Air Resistance, Friction, and Buoyancy...... 8 1.6 Charged Particle in a Magnetic Field...... 17 1.7 Summary: Newton’s Laws...... 20
2 Momentum and Center of Mass 23 2.1 Rocket with no External Force ...... 24 2.2 Multistage Rockets...... 25 2.3 Rocket with External Force...... 25 2.4 Center of Mass...... 26 2.5 Angular Momentum of a Single Particle...... 28 2.6 Angular Momentum of a System of Particles ...... 29 2.7 Angular Momentum of a Continuous Mass Distribution ...... 30 2.8 Summary: Momentum and Center of Mass ...... 34
3 Energy 36 3.1 General One-dimensional Systems ...... 40 3.2 Atwood Machine...... 43 3.3 Spherically Symmetric Central Forces ...... 44 3.4 The Energy of a System of Particles ...... 45 3.5 Summary: Energy ...... 47
4 Oscillations 49 4.1 Simple Harmonic Oscillators...... 49 4.2 Damped Harmonic Oscillators...... 52 4.3 Driven Oscillators ...... 57 4.4 Summary: Oscillations...... 69
5 Calculus of Variations 72 5.1 The Euler-Lagrange Equation...... 72 5.2 The Brachistochrone Problem...... 77 5.3 General Parametrization...... 80 5.4 Summary: Calculus of Variations...... 83
6 Lagrangian Mechanics 84 6.1 2D Particle in Cartesian Coordinates...... 84 6.2 2D Particle in Polar Coordinates...... 86 iv Contents
6.3 Unconstrained Particles in 3D...... 87 6.4 Constrained Particles in 3D...... 87 6.5 Noether’s Theorem...... 95 6.6 The Hamiltonian...... 96 6.7 Lagrange Multipliers...... 98 6.8 Summary: Lagrangian Mechanics...... 100
7 Orbits 102 7.1 The Kepler Problem...... 105 7.2 Kepler’s Laws...... 109 7.3 Transfer Orbits...... 109 7.4 Summary: Orbits...... 113
8 Noninertial Frames 115 8.1 Frame with Linear Acceleration...... 115 8.2 Rotating Frames ...... 116 8.3 Centrifugal Force...... 118 8.4 Coriolis Force...... 119 8.5 Summary: Noninertial Frames...... 123
9 Rigid Rotations 125 9.1 Rotation About a Fixed Axis...... 126 9.2 The Inertia Tensor...... 128 9.3 Principal Axes ...... 133 9.4 Torque-free Motion...... 136 9.5 Summary: Rigid Rotations ...... 139
10 Coupled Oscillators 141 10.1 Two Masses and Three Springs...... 141 10.2 The General Problem ...... 145 10.3 Double Pendulum...... 147 10.4 Triatomic Molecule...... 149 10.5 Parallel and Series Springs...... 151 10.6 Summary: Coupled Oscillators ...... 152
11 Collision and Scattering 154 11.1 Hard Sphere Scattering ...... 155 11.2 The General Case...... 156 11.3 Rutherford Scattering ...... 158
12 Math Reference 160 12.1 Complex Numbers ...... 160 12.2 Vectors ...... 160 12.3 Curvilinear Coordinates...... 160 12.4 Differential Equations ...... 163 12.5 Taylor Series ...... 163 12.6 Approximations...... 164
Index 165 Preface
About These Notes
Cover image: NASA (The International Space Station) These are my class notes from two courses on classical mechanics (PHY 3221 and PHY 4222) that I took at Florida State University. Our textbook was Classical Mechanics by John R. Taylor. My class notes can be found at www.leonhostetler.com/classnotes Please bear in mind that these notes will contain errors. Any errors are certainly my own. If you find one, please email me at [email protected] with the name of the class notes, the page on which the error is found, and the nature of the error. This work is currently licensed under a Creative Commons Attribution-NonCommercial- NoDerivatives 4.0 International License. That means you are free to copy and distribute this document in whole for noncommercial use, but you are not allowed to distribute derivatives of this document or to copy and distribute it for commercial reasons.
Updates
Last Updated: July 10, 2019 Version 0.82: (Jul. 10, 2019) Updated layout, spell-checked Version 0.81: (Dec. 21, 2017) First upload.
Conventions
I will represent a vector by a bold letter topped with an arrow. For example: F~. I will represent a unit vector (i.e. a direction vector with magnitude 1) by a bold leter topped with a hat. For example: xˆ. In physics, a time derivative is often represented with a dot. For example, instead dx d2x of writing dt or dt2 , one might writex ˙ orx ¨. This convention will generally be used in these notes.
Chapter 1
Newton’s Laws
Mechanics is the study of the motion of material bodies in space as a function of time. We’ll be studying classical mechanics which assumes v << c (i.e. non-relativistic) and Planck’s constant is treated as h = 0 (i.e. non-quantum). First Law: In the absence of forces, a body moves with uniform velocity. It’s the ability to visualize the ideal case that enabled Newton to achieve this break- through since in the real world, moving bodies tend to slow down. Remarkably, a state- ment about an ideal situation can allow us to make arbitrarily precise predictions. Second Law: F~ = m~a This law does not hold relativistically. Here, m is the inertial mass, which is an attribute of the body. Inertial mass quantifies a body’s resistance to force. F~ is not an attribute of the body, but rather it quantifies the interaction of the body with its environment. Newton also introduced a “quantity of motion” that is known today as momentum, ~p = m~v. Using this, we can write Newton’s law in a more general form that also holds relativistically: d~p F~ = = ~p˙. dt Third law: For every action, there’s an equal but opposite reaction. The notation F~ ij denotes the force exerted on particle i by particle j. Read it as “force on i by j”. According to Newton’s third law, for two particles i and j,
F~ ij = −F~ ji.
ext If a system of two particles is isolated with no outside forces, that is, F~ = 0, then the only force on any of the two particles is the force exerted by the other particle. Since ~ ~ ~ ˙ ~ ˙ F ij = −F ji and F 12 = ~p1 and F 21 = −~p2, we have that d (~p + ~p ) = 0, dt 1 2 or conservation of momentum. This is what the third law is really telling us. The weight force of a body is ~w = m~g. The mass m in this case is gravitational mass. Although gravitational mass and inertial mass are defined differently, for any object ever measured, the two always had the same value. In fact, Einstein’s principle of equivalence states that
gravitational mass = inertial mass.
Some important concepts implied by Newton’s laws are 2 Newton’s Laws
Inertial frames: An inertial frame is any frame in which Newton’s first law holds. Fic- titious forces arise when not in an inertial frame. Galilean invariance: Any frame in uniform motion with respect to an inertial frame is also an inertial frame. Equation of motion: The equation of motion for bodies is Newton’s second law. The second law can be written as a second order differential equation with respect to position F~ = m~r¨. To get the position equation for a body, we integrate this twice ˙ with respect to time—getting two constants of integration ~r0 and ~r0 = ~v0. These constants are often fixed by initial conditions.
1.1 Many Particles
Suppose you have many particles spread throughout some space. Each particle, in this case, is experiencing an external force from outside the system as well as a force from every other particle. If you define some origin O, then the position of the ith particle can be represented by ~ri, and its momentum by ~pi = mi~ri. Newton’s second law implies d X ext ~p = F~ + F~ . dt i ij i j6=i What this states is that the time derivative of the momentum of the ith particle is equal to the sum of all the forces due to the other particles plus the external force on the particle. The index j 6= i simply means to sum over all j not equal to i. That is, sum over all other particles—ignoring the case where F~ ii since a particle doesn’t exert a force on itself. The total momentum of the system, which we denote with a capital P~ is the sum of the momenta of all the particles ~ X P = ~pi. i Since we can differentiate sums term by term, we know that ~˙ X ˙ P = ~pi. i ˙ Replacing ~pi with our result from above gives us ~˙ X X ~ ~ ext X X ~ X ~ ext P = F ij + F i = F ij + F i . i j6=i i j6=i i Looking at the first term on the right, we can expand this as X X X X X X F~ ij = F~ ij + F~ ij. i j6=i i ji We can write the right side more compactly by combining the two terms and switching the subscripts on one of the force vectors X X X X F~ ij = F~ ij + F~ ji . i j6=i i j
This simplifies to ˙ ext P~ = F~ .
In other words, the rate of change of the total momentum of the system depends only on the external forces on the system—not on the internal forces. This also implies that if the external forces on a system of particles is zero, then the total momentum is a constant (the indefinite integral of 0 is a constant).
1.2 Cartesian Coordinates
A law of motion is the second order differential vector equation for the position vector ~r as a function of time F~ = m~r¨.
In Cartesian coordinates, the position vector is
~r = xxˆ + yyˆ + zzˆ = hx, y, zi, where x, y, and z are generally functions of t. Differentiating the position function is clean and easy since the unit vectors do not depend on time—they are just constants when we differentiate. Since acceleration is the second derivative of the position vector with respect to time, we have that
~r¨ = ~a =x ¨xˆ +y ¨yˆ +z ¨zˆ.
Since the components in the three directions are separable in Cartesian coordinates, we can write a force vector in the form
F~ = Fxxˆ + Fyyˆ + Fzzˆ.
This means we can treat Newton’s second law in 3D as three separate versions of Newton’s second law in 1D F = mx¨ x ~ F = m~a =⇒ Fy = my¨ Fz = mz¨
These three ordinary differential equations (ODEs) are separable, so finding the solutions x(t), y(t), and z(t) is a simple matter of separating variables and integrating twice. For example, integrating Fx = mx¨, gives us the velocity function and a constant—the initial speed vx(t = 0) = vx0 . Integrating the result of that gives us the position function and a second constant—the initial position x(t = 0) = x0
F a =x ¨ = x = constant x m v (t) =x ˙ = x¨ dt = a dt = a t + v x ˆ x ˆ x x0 1 x(t) = (a t + v ) dt = a t2 + v t + x . ˆ x x0 2 x x0 0
We can use the same process in the y and z directions to find vy(t), vz(t), y(t), and z(t). 4 Newton’s Laws
1.3 Polar Coordinates
φˆ rˆ If one or more of the forces involved is a central force, it may be easier to use polar • coordinates. Recall the transformation equations when going from polar to Cartesian or vice versa ~r x = r cos φ, y = r sin φ φ p y r = x2 + y2, φ = tan−1 . x Figure 1.1: Polar coordinates where r is the distance to the particle, and φ is the angle between ~r and the positive x-axis. The position vector in 2D polar coordinates, i.e., z = 0 is
~r = hx, y, zi = hr cos φ, r sin φ, 0i.
To find the velocity and acceleration vectors, we differentiate the position vector. This time it’s more complicated because r and φ are both functions of time, so we have to use the product and chain rules to differentiate the position vector to get
~v = ~r˙ = hr˙ cos φ − rφ˙ sin φ, r˙ sin φ + rφ˙ cos φ, 0i =r ˙hcos φ, sin φ, 0i + rφ˙h− sin φ, cos φ, 0i.
Differentiating again to get the acceleration vector gives us the following after rearranging
~a = (¨r − rφ˙2)hcos φ, sin φ, 0i + (2r ˙φ˙ + rφ¨)h− sin φ, cos φ, 0i.
This is a lot messier than it was with Cartesian coordinates. We can write these in a more condensed manner by noting that hcos φ, sin φ, 0i is the unit vector rˆ in the direction of ~r. Furthermore, h− sin φ, cos φ, 0i, is a unit vector that is at right angles to rˆ and points in the direction of increasing φ. We can verify that they’re at right angles to each other by noting that their dot product is zero. We know it’s pointing in the direction of increasing φ because when φ = 0, it is h0, 1, 0i. Therefore, we call this unit vector φˆ. So our unit vectors in polar coordinates are
rˆ = hcos φ, sin φ, 0i φˆ = h− sin φ, cos φ, 0i.
This allows us to rewrite our position, velocity, and acceleration in terms of the unit vectors rˆ and φˆ
~r = r rˆ ~v = ~r˙ =r ˙ rˆ + rφ˙ φˆ ~a = ~r¨ = r¨ − rφ˙2 rˆ + 2r ˙φ˙ + rφ¨ φˆ.
Notice, that unlike the unit vectors in Cartesian coordinates, these are changing with time. It’s as if the particle is carrying them along. How are the unit vectors changing with time? Differentiating them gives us d rˆ = φ˙h− sin φ, cos φ, 0i = φ˙φˆ dt d φˆ = φ˙h− cos φ, − sin φ, 0i = −φ˙rˆ. dt 1.3. Polar Coordinates 5
So the time derivatives of the unit vectors are
rˆ˙ = φ˙φˆ φˆ˙ = −φ˙rˆ.
This gives us an alternate way of writing ~v by replacing φ˙φˆ with rˆ˙ d ~v = ~r =r ˙ rˆ + rφ˙ φˆ =r ˙ rˆ + rrˆ˙, dt and an alternate way of writing ~a by differentiating the new form of ~v d d ~a = ~v = (r ˙ rˆ + rrˆ˙) =r ¨rˆ + 2r ˙rˆ˙ + rrˆ¨. dt dt We can write Newton’s second law in component form as
F~ = Frrˆ + Fφφˆ, where Fr is the force in the direction of rˆ and Fφ is the force in the direction of φˆ. This is not usually convenient since the force may be changing with time, since it depends on φ and r. From the law F~ = m~a and the equation ~a = (¨r − rφ˙2)rˆ + (2r ˙φ˙ + rφ¨)φˆ, we can separate Newton’s law into two components ˙2 Fr = m r¨ − rφ F~ = m~a =⇒ ˙ ¨ Fφ = m 2r ˙φ + rφ
Example 1.3.1
Using polar coordinates, analyze the motion of a mass m being twirled around on a string of length R. From Newton’s second law in polar coordinates, we have that
˙2 ˙ ¨ F~ = Frrˆ + Fφφˆ = m r¨ − rφ rˆ + m 2r ˙φ + rφ φˆ
Since the length of the string R is fixed, we have the r = R is constant. This implies thatr ˙ =r ¨ = 0, giving us F~ = −m Rφ˙2 rˆ + m Rφ¨ φˆ.
We know that the centripetal radial force is being applied by the string tension. Since the tension force points in the direction opposite of rˆ, we have that Fr = −T . We also know that there is no force in the φˆ direction in this idealized situation. So we have that F~ = −T~ T~ = mRφ˙2rˆ. ¨ ¨ ˙ Since Fφ = mRφ = 0, we know that φ = 0, which implies that φ is a constant. We know this constant as the angular velocity ω, so we have that T = mRω2, where Rω2 is the centripetal acceleration associated with uniform circular motion. This tells us that the rφ˙2 from the equation for ~a is just our familiar Rω2. 6 Newton’s Laws
Example 1.3.2
Using polar coordinates, analyze the motion of a simple pendulum with a weight of mass m on a string of length R. From Newton’s second law in polar coordinates, we have that
˙2 ˙ ¨ F~ = Frrˆ + Fφφˆ = m r¨ − rφ rˆ + m 2r ˙φ + rφ φˆ
The forces on the mass are T , the tension force in the string and ~w = m~g the weight. Using a free body diagram, we find that
Fr = mg cos φ − T
Fφ = −mg sin φ.
Notice that T is negative since the tension force is pointing in the direction opposite of increasing r. Also, g is negative in Fφ since it points in the direction of decreasing φ. Equating these components with the ones we get from Newton’s law gives us
˙2 Fr = m r¨ − rφ = mg cos φ − T ˙ ¨ Fφ = m 2r ˙φ + rφ = −mg sin φ.
Since our string length is fixed, we have that r = R is a constant, sor ˙ =r ¨ = 0, giving us
˙2 Fr = −mRφ = mg cos φ − T ¨ Fφ = mRφ = −mg sin φ.
Solving Fr for T gives us T = m Rφ˙2 + g cos φ .
¨ Solving Fφ for φ gives us g φ¨ = − sin φ, R which is a second order nonlinear differential equation. For |φ| << 1, we can use the small angle approximation sin φ ≈ φ to write g φ¨ ≈ − φ. R or the more familiar φ¨ ≈ −ω2φ, where ω = pg/R is the angular frequency. We know the general solution to this second order differential equation is
φ(t) = A cos ωt + B sin ωt. 1.4. Cylindrical Coordinates 7
z 1.4 Cylindrical Coordinates ρ In cylindrical coordinates, P 0 is the projection of a point P onto the xy-plane. Then φ is the angle between this projection and the positive x-axis, and ρ is the length of this • projection. Notice that ρ is not the distance from a particle at P to the origin, rather, it P is the distance from the particle P to the z-axis. The transformations between Cartesian and cylindrical coordinates are y x = ρ cos φ, y = ρ sin φ, z = z ρ p y ρ = x2 + y2, φ = tan−1 . φ • 0 x x P Figure 1.2: Cylindrical coordi- We can define directions in cylindrical coordinates in terms of the unit vectors ρˆ, nates φˆ and zˆ. These vectors are orthonormal (i.e. mutually orthogonal and of unit length). The unit vector ρˆ points in the direction of increasing ρ. The unit vector φˆ is orthogonal to ρˆ and points in the direction of increasing φ. Notice that φˆ and ρˆ are parallel to the xy-plane. The unit vector zˆ is orthogonal to both of the others and points parallel to the z-axis in the direction of increasing z. Their components in terms of Cartesian coordinates are
ρˆ = hcos φ, sin φ, 0i φˆ = h− sin φ, cos φ, 0i zˆ = h0, 0, 1i.
The position vector for a particle in cylindrical coordinates can be written as
~r = ρρˆ + zzˆ.
Differentiating this twice gives us
~r˙ =ρ ˙ρˆ + ρρˆ˙ +z ˙zˆ ~r¨ =ρ ¨ρˆ + 2ρ ˙ρˆ˙ + ρρˆ¨ +z ¨zˆ.
From 2-dimensional polar coordinates, we know that ρˆ˙ = φ˙φˆ and differentiating this gives us ρˆ¨ = φ¨φˆ + φ˙φˆ˙, and we know that φˆ˙ = −φ˙ρˆ, so we can rewrite the above two equations as
~r˙ = ~v = ρρˆ + ρφ˙φˆ +z ˙zˆ ~r¨ = ~a = ρ¨ − ρφ˙2 ρˆ + ρφ¨ + 2ρ ˙φ˙ φˆ + (¨z) zˆ.
Example 1.4.1
A particle of mass m is constrained to move at a distance of ρ = R from the vertical z-axis. Using Newton’s second law, find the equations of motion for the particle. Newton’s second law in cylindrical components can be broken into components
F~ = Fρρˆ + Fφφˆ + Fzzˆ.
Using F~ = m~a and the equation for ~a in cylindrical coordinates, we can write the 8 Newton’s Laws
force components as
˙2 Fρ = maρ = m ρ¨ − ρφ ¨ ˙ Fφ = maφ = m ρφ + 2ρ ˙φ
Fz = maz = m (¨z) .
From the description of the problem, the only forces acting on the particle are the weight force ~w = m~g which is acting entirely in the negative z direction and the centripetal force, which is directed toward the z-axis.
v2 F = m ρ¨ − ρφ˙2 = −m T = −mρω2 ρ ρ ¨ ˙ Fφ = m ρφ + 2ρ ˙φ = 0
Fz = m (¨z) = −mg.
We know that ρ = R is a constant which implies thatρ ˙ =ρ ¨ = 0. Simplifying the equations gives us
φ˙ = ω φ¨ = 0 z¨ = −g.
Integrating the first equation gives us
φ(t) = ωt + φ0.
Integrating the third equation gives us
z˙ = vz(t) = −gt + vz0 .
The constant of integration is C = vz0 since vz(t = 0) = vz0 . Integrating again gives us 1 z(t) = − gt2 + v t + z . 2 z0 0
1.5 Air Resistance, Friction, and Buoyancy
Incorporating air resistance is a good way to introduce differential equations into our equations of motion. We know from experience that air resistance is a function of velocity, f(v). The faster you go through the air, the more air resistance you experience. Empirically, the magnitude of the force f~(v) of air resistance is
f(v) = bv + cv2,
2 where flin = bv is the linear air resistance and fquad = cv is the quadratic air resistance. The form of this function makes sense since any function f(v) will have a Taylor series expansion of the form f(v) = a + bv + cv2 + ··· , where a, b, and c are constants. In the case of air resistance, the leading constant a is zero since there is no air resistance when an object is not moving. 1.5. Air Resistance, Friction, and Buoyancy 9
We can also assume that the direction of the air resistance is opposite the direction of the object’s velocity. This is exactly true for non-rotating, spherical objects, but is ~v at best an approximation valid for non-rotating or non-spherical objects. Consider an airplane wing for example. A good portion of the air resistance force, the lift, is not in • the direction opposite of the plane’s motion. In our case, we’ll assume that
f~ = −f(v)vˆ. ~ f = −f(v)vˆ ~w = m~g Recall that vˆ = ~v/v. We’ll assume the object is a sphere. The constants b and c depend on the character- Figure 1.3: Air resistance istics of the object and the fluid. Empirically, we find that linear air resistance is due to the viscosity (i.e. viscous drag) of the fluid and the size of the object. ~v b ∝ viscosity · D, where D is the diameter of the sphere. It turns out from Stokes’ law that • b = 3πηD, flin = −b~v ~w = m~g where η is the viscosity of the fluid. We also find that quadratic air resistance is due to density of the fluid and the Figure 1.4: Linear air resistance cross-sectional area of the sphere c ∝ density · D2. We can understand this as meaning that quadratic air resistance is due to the fact that the object must accelerate the mass of the fluid that it collides with. Obviously, then, it will be related to both the density of the fluid and the cross-sectional area (D2) of the object. For a sphere in air at STP, we find that b = 1.6 × 10−4 (kg/ms) D c = 0.25 (kg/m3) D2.
Looking at the ratio of fquad and solving for v, we find that if flin 6 × 10−4 m2/s v << , then f dominates D lin 6 × 10−4 m2/s v >> , then f dominates. D quad
For typical everyday objects moving at everyday velocities, fquad dominates.
Linear Air Resistance For a spherical object, the linear air resistance is
~ f lin = −b~v. The sum of forces acting on the particle are
F~ = m~g − b~v.
Linear air resistance is simple because we can still separate the components when using Cartesian coordinates
Fx = −bvx
Fy = −mg − bvy. 10 Newton’s Laws
Horizontal Motion with Linear Air Resistance
~v We can use the same process to consider friction for an object sliding on a plane (as shown in Fig. (1.5)) when the friction depends on the velocity. We will use this idea to f~ = −b~v m investigate the horizontal component of linear air resistance. Applying Newton’s law, we get for the x component, the separable differential equa- tion
Figure 1.5: 1D horizontal mo- F = mv˙ = −bv tion with friction. x x x dv b x = −kv , k = dt x m 1 1 vx(t) − dvx = dt k ˆ vx ˆ 1 t = − ln v + C. k x
vx0 Then we fix the arbitrary constant of integration C using the initial conditions vx(t =
0) = vx0 . In other words, when t = 0, vx = vx0 . Plugging these values in gives us 1 t 0 = − ln vx0 + C 0 τ k Figure 1.6: The exponential de- 1 C = ln vx , cay of speed for an object sliding k 0 with linear friction. so our complete equation is 1 1 1 vx t = − ln vx + ln vx0 = − ln . k k k vx0 x(t) What we actually want is the velocity in terms of t instead of the other way around. To invert the equation, we solve the right side for the natural logarithm, then raise both sides as powers to e to get
−kt − t v (t) = v e = v e τ , x0 + vx0 τ x x0 x0 where 1 m τ = = , t k b is an example of a characteristic time. x0 τ Figure 1.7: The position func- Notice that our function vx(t), shown in Fig. (1.6), is an exponential function. A tion of an object sliding with related function is the radioactive decay function where the characteristic time is the linear friction. decay constant. Another related function appears in circuits, where for example, the time constant might give the time it takes for a capacitor in a circuit to discharge to 1/e = 36.8% of its initial charge. In our case, τ is the time it takes for the velocity of the object to decay to 36.8% of its initial value. Integrating vx(t) gives us
− t − t x(t) = v e τ dt = −v τe τ + C. x0 ˆ x0
Plugging in the initial condition x(t = 0) = x0, gives us x0 = −vx0 τ +C, so C = x0 +vx0 τ and our complete equation is
t − τ x(t) = x0 + vx0 τ 1 − e .
Notice that as t → ∞, x − x0 → vx0 τ, that is, vx0 τ is how far the object can travel. This is illustrated in Fig. (1.7). 1.5. Air Resistance, Friction, and Buoyancy 11
Check various limits of x(t) to make sure the solution behaves the way you would expect it to. If vx0 = 0 then x(t) = x0. This makes sense. If the object has no initial velocity, it stays at its starting place. What does the solution tell us if air resistance is negligible? When b → 0, τ = m/b → ∞ and t/τ → 0. Then t t2 x(t) = x + v τ 1 − 1 − + + ··· ' x + v t. 0 x0 τ 2τ 2 0 x0
Notice that we replaced e−t/τ with its Taylor expansion. When the air resistance is negligible, our result reduces to the familiar equation for 1D motion with constant velocity. −t/τ If we look at an infinitely viscous fluid, that is, τ → 0, we see that e → 0, so vx = 0 and x = x0. That is, the object goes from v = vx0 to v = 0 instantaneously. In other words, it comes to a stop instantaneously. No matter how complicated the air resistance function is, the object will eventually slow until flin dominates, and when it does, it has a finite stopping distance—vx0 τ. In general, for 1D motion with only a velocity-dependent force F (v), Newton’s second law gives us dv m v dv0 m = F (v) =⇒ dt = dv =⇒ t = m 0 . dt F (v) ˆv0 F (v ) Doing the integration gives you t as a function of v, t = t(v). You then need to invert the function to get v as a function of t, v = v(t). Once you have v(t), you can calculate the position function by integrating
t 0 0 x − x0 = v(t ) dt . ˆ0
Vertical Motion with Linear Air Resistance We turn now to vertical motion with linear air resistance—as with an object falling through the air. With vertical motion, the forces acting on the body are weight and air resistance. From Newton’s second law, for upward motion, we get
Fy = mv˙y = −mg − bvy. Notice that we chose the upward direction to be positive y. This becomes the first order separable differential equation dv b y = −g − kv , k = . dt y m Separating the differentials and integrating gives us dv dt = − y g + kvy t vy 1 dt0 = − dv0 ˆ ˆ g + kv0 y 0 vy0 y 1 vy 1 g + kv 0 y t = − ln(g + kvy) = − ln k vy0 k g + kvy0
We now have t in terms of v, but we want v in terms of t, so solving for vy gives us
g + kvy −kt − t = e = e τ , g + kvy0 which simplifies to g g − t v (t) = − + + v e τ . y k k y0 12 Newton’s Laws
−t/τ Notice that for t >> τ, e → 0, so vy(t) → −g/k = −vter, the terminal velocity of the −vy(t) body, which can be written mg g v = = = gτ, (linear air resistance). ter b k Another way to find the terminal velocity is to realize that it occurs when the acceleration v ter of a falling body is zero, that is, when the weight force is balanced by the air resistance, so bvter = mg =⇒ vter = mg/b = gτ. The expression for terminal velocity allows us to rewrite the equation for vy(t) as
t − τ t vy(t) = −vter + (vter + vy0 ) e . 0 τ If vy0 = 0, then the velocity of the body approaches the terminal velocity in the manner Figure 1.8: Velocity of an ob- shown in Fig. (1.8). ject falling against linear air re- To get the vertical position function of the object, we just integrate vy(t) sistance. t t 0 t 0 0 0 − t y − y = v (t ) dt = −v t − τ (v + v ) e τ . 0 y ter ter y0 ˆ0 0 0 This simplifies to t − τ y(t) = −vtert − τ (vter + vy0 ) e − 1 .
If we calculate the terminal velocity of an object assuming that flin is the dominant resistive force, we should check that our assumption is valid. Our assumption is valid if
fquad << flin, when calculated at the terminal velocity. For an object falling in a dense fluid, such as a liquid, another factor comes into play—that of buoyancy. The buoyant force on such an object is equal to the weight of the displaced fluid
Fb = ρV g, where ρ is the density of the fluid, V is the volume of the object, and g is the local gravitational acceleration. The buoyant force is directed opposite the gravitational force. The net force on an object falling vertically in a dense fluid and experiencing only linear drag, is dv F = m y = 3πηDv + ρV g − mg. y dt y We’ve done horizontal and vertical linear resistance, and since they decouple, we can put them together to get the entire equation of projectile motion with linear air resistance. For the horizontal and vertical motions with x0 = y0 = 0, we have that
t − τ x(t) = vx0 τ 1 − e
t − τ y(t) = −vtert + τ (vter + vy0 ) 1 − e .
Recall that the range for a projectile, when not including air resistance, is 2v v R = x0 y0 . vac g What is the range when air resistance is taken into account? To find the shape of the
trajectory, we solve x(t) for t to get t = −τ ln(1 − x/(vxo τ)), which we plug in to y(t) and simplify to get x vter + vy0 y(x) = vterτ ln 1 − + x. vxo τ vx0 1.5. Air Resistance, Friction, and Buoyancy 13
Fig. (1.9) shows the trajectory of a body with linear air resistance (solid line) and no air resistance (dashed line). Looking at the equation for x(t), we see that the limit as y x = vx0 τ t → ∞ is x(t) = vx0 τ, so the line x = vx0 τ is our vertical asymptote. Remember that x will never actually reach vx0 τ, but it asymptotically approaches it. The range R of the projectile under linear air resistance is the x-value when y(x) = 0, so R vter + vy0 0 = vterτ ln 1 − + R. x vxo τ vx0 y0 We can’t find a simple closed-form expression for R but we can approximate it by finding the first correction term for the limit of weak air resistance. In the limit of weak air resistance, τ → ∞, so Figure 1.9: The trajectory of a R projectile under linear air resis- << 1. tance. vx0 τ Using the Taylor expansion
1 1 ~v ln(1 − x) = − x + x2 + x3 + ··· , 2 3 • and replacing vter with gτ, we get ~ 2 ! f quad = −cv vˆ ~w = m~g R 1 R 2 1 R 3 gτ + v 0 = −gτ 2 + + + ··· + y0 R vxo τ 2 vxo τ 3 vxo τ vx0 Figure 1.10: Quadratic air resis- gτR gR2 gR3 gτR v R tance = − − − − · · · + + y0 v 2v2 3v3 τ v v xo xo xo x0 x0 gR2 gR3 v R gR2 gR3 v R = − − − · · · + y0 ≈ − − + y0 2v2 3v3 τ v 2v2 3v3 τ v xo xo x0 xo xo x0 2v v 2 R ≈ x0 y0 − R2. g 3vxo τ We now have a quadratic equation in R, and we could solve for R using the quadratic formula, but there’s no need to. We have only an approximation here, and there’s no need to solve an equation exactly if the equation is itself only approximate. We can approximate it by iterating on a calculator. We start by calculating the first order approximation of R, that is, we ignore the second term in the equation above and then we improve that approximation by plugging that result in for R
2v v 2 2v v 2 2v v 4v R ≈ x0 y0 − x0 y0 ≈ x0 y0 1 − y0 . g 3vxo τ g g 3vter Note: This kind of approximation has its limitations. A quadratic equation can have two solutions. This method only works when there are two solutions and there is a large difference between the two solutions. This approximation method will then converge to the solution that is closest to zero.
Quadratic Air Resistance Quadratic air resistance is given by
~ 2 f quad = −cv vˆ = −cv~v.
From Newton’s second law, we get
d~v F = m = m~g − cv~v. dt 14 Newton’s Laws
q 2 2 The complication is the fact that v = vx + vy depends on both the x and y components. Newton’s second law for each component gives us c q v˙ = − v2 + v2v x m x y x c q v˙ = − v2 + v2v − g. y m x y y Now they are nonlinear and coupled, so we can’t solve them separately as we did with linear air resistance. Fortunately, if we focus only on one-dimensional motion, we don’t have to worry about the orthogonal component.
Horizontal Motion with Quadratic Air Resistance From Newton’s second law, we get dv F = m = −cv2. dt We solve this in the usual way, integrating and then solving for v in terms of t to get
v0 m v(t) = t , τ = . 1 + τ cv0
Notice that the characteristic time τ for quadratic air resistance is different from the one for linear air resistance. This one depends on the initial speed. In both cases, though, the characteristic time gives us the time frame in which the object slows appreciably. If we plot v(t), we find that there’s a slower decay than for linear air resistance. That is, an object under quadratic air resistance, doesn’t slow as quickly as an object under linear air resistance. At time τ, the object has a speed exactly half of its initial speed. To find the position function, we integrate v(t) as before to get
t x(t) = x + v τ ln 1 + . 0 0 τ
Notice that this function grows logarithmically, surpassing any finite value of x. So under just quadratic air resistance, an object would never slow to a complete stop. In the real world, any object under quadratic air will eventually slow enough that linear air resistance takes over and brings it to a complete stop in a finite distance.
Vertical Motion with Quadratic Air Resistance For vertical motion, we have two forces—the weight force pointing downward and the quadratic air resistance pointing opposite the direction of motion. From Newton’s second law, for an object going up, we have dv F = m = −cv2 − mg, dt if we take the upward direction to be positive y. For an object going down, we have dv F = m = cv2 − mg, dt Because of the squared term v2, our equation won’t change sign, as it should, when the object goes from upward motion to downward motion. Because of that, we need to treat the two directions separately. In other words, for a single object thrown upwards, we 1.5. Air Resistance, Friction, and Buoyancy 15 need two equations—one for the upward portion of the trajectory and the other for the downward portion. For downward motion, the object reaches terminal speed when the air resistance equals the gravitational force. So we get that terminal velocity under quadratic air resistance is rmg v = (quadratic air resistance). ter c Looking at downward motion, we have that
2 2 dv c 2 cv v = v − g = −g 1 − = −g 1 − 2 dt m mg vter 1 dt = − dv v2 g 1 − 2 vter t 1 v 1 dt0 = − dv0. v2 ˆ0 g ˆv0 1 − 2 vter
0 0 On the right side, we can make the substitution u = v /vter then du = (1/vter)dv and 0 dv = vter du. v vter vter 1 t(v) = − 2 du. g ˆ v0 1 − u vter We could use partial fraction decomposition, 1/(1 − u2) = 1/(2(1 − u)) + 1/(2(1 + u)), to rewrite and solve it as logarithms, but it gets very messy. From a table of integrals, we find that the integral is the inverse of the hyperbolic tangent, so
v v v v v v t = − ter tanh−1 u ter = − ter tanh−1 − tanh−1 0 . g v0 g v v vter ter ter Solving for v gives us −1 v0 gt vd(t) = vter tanh tanh − . vter vter
For v0 = 0, that is, the object is dropped from rest, this simplifies to
gt vd(t) = −vter tanh . vter
What happens in the limit of weak air resistance? For weak air resistance vter >> gt, since we know that for no air resistance, the terminal velocity is infinite. From the Taylor expansion tanh x = x − x3/3 + ··· , we see that tanh x ≈ x when x is small. So
gt v(t) = −vter + ··· ≈ −gt, vter and we know that v(t) = −gt when there is no air resistance, so this limit checks out. Finally, we can find the position function for an object dropped from rest by inte- grating v(t). The result is
2 vter gt yd(t) = y0 − ln cosh . g vter
For the case of upward motion with quadratic air resistance, we’ll use a slightly different trick. For upward motion, the air resistance and the gravitational force are both 16 Newton’s Laws
opposite the direction of the velocity. We start again with Newton’s second law, giving us dv m = −mg − cv2 dt dv c v2 = −g 1 + v2 = −g 1 + . dt mg vter Notice that we are again incorporating the terminal velocity despite this being upward motion. Terminal velocity is the asymptotic velocity for downward motion, but it’s also a useful quantity when dealing with upward motion. It is still defined in terms of downward motion. We want to find the velocity v as a function of height y (instead of time t). We could solve for v(t) and y(t), eliminate t and plug one into the other to get v(y), but it’s easier to use the chain rule dv dv dy dv = · = v · . dt dy dt dy The key is in noticing that dy/dt is just the velocity v. Using this, we get dv v2 v · = −g 1 + 2 dy vter 1 v dy = − 2 2 dv g 1 + v /vter y v 0 1 v 0 dy = − 02 2 dv . ˆy0 g ˆv0 1 + v /vter 02 2 We can make the substitution u = 1 + v /vter, then
2 2 2 2 2 1+v /vter 2 2 2 2 v 1 v 1+v /vter v 1 + v /v ter ter ter 0 ter y − y0 = − du = − ln u = ln 2 2 . 2g ˆ 2 2 u 2g 1+v2/v2 2g 1 + v /v 1+v0 /vter 0 ter ter So 2 2 2 vter 1 + v0/vter y(v) = y0 + ln 2 2 . 2g 1 + v /vter To find v(y), we just have to solve for v. What is the maximum height of the object if it’s initial upward velocity is v0? Since the maximum height occurs when v = 0, we can just replace v with 0 in the equation above to get 2 2 vter v0 ymax = y(0) = y0 + ln 1 + 2 . 2g vter
For weak air resistance, vter >> v0. The Taylor expansion of ln(1 + x) is 1 1 1 ln(1 + x) = x − x2 + x3 − x4 + ··· 2 3 4 so 2 " 2 2 2 # 2 2 vter v0 1 v0 v0 1 v0 ymax = y0 + 2 − 2 + ··· = y0 + 1 − 2 + ··· . 2g vter 2 vter 2g 2 vter
So for an object starting at y0 = 0, we have that 2 2 v0 1 v0 ymax ≈ 1 − 2 . 2g 2 vter Notice that the first term is our familiar result for no air resistance and the second term (i.e. the first correction term) is negative, which is what we’d expect since with air resistance, the object will not travel quite as high as without air resistance. 1.6. Charged Particle in a Magnetic Field 17
General Motion with Quadratic Air Resistance In general, for motion that is not just horizontal or vertical, Newton’s second law gives
m~v˙ = m~g − cv2vˆ = m~g − cv~v.
The components are then q ˙ 2 2 m~vx = −c vx + vy vx q ˙ 2 2 m~vy = −mg − c vx + vy vy.
Notice that the x component of the velocity depends also on the y-component and vice versa, i.e., the equations are coupled. These equations can’t be solved analytically, but they can be solved numerically given specific initial conditions. We know, however, that after a long time, vx → 0 and vy → −vter. Plugging these values into the pair of equations gives us their limit after a long time ˙ m~vx = −cvtervx ˙ m~vy = −mg − cvtervy.
This tells us that if we wait a long time, the coupled equations eventually behave like the uncoupled case of linear air resistance. Motion in the x-direction asymptotically approaches a finite value.
1.6 Charged Particle in a Magnetic Field
Consider a particle with charge q moving in a magnetic field in the positive z-direction, so B~ = Bzˆ. The force on the charged particle is
F~ = q~v × B~ .
Since B~ has only a component in the z direction, the cross product is easily calculated as
xˆ yˆ zˆ ~ F = q vx vy vz = q (vyBxˆ − vxByˆ) .
0 0 B
Newton’s second law gives us
mv˙x = qBvy, mv˙y = −qBvx, mv˙z = 0.
We can immediately see that the acceleration in the z-direction is 0 and therefore, the velocity in that direction is the constant initial velocity in that direction. This gives us the equation of motion in the z direction
z(t) = z0 + vz0 t.
Notice that v˙x = ωvy, v˙y = −ωvx, where ω = qB/m is the cyclotron frequency. If we differentiate the first equation and plug it into the second equation, we get
2 v¨x = −ω vx, 18 Newton’s Laws
a second order differential equation, which we recognize as the equation for oscillatory motion, and it has the solution
vx = A cos ωt + B sin ωt.
A and B are the two constants of integration, which are determined by the initial condi- tions. If we let A = a cos δ and B = a sin δ, then we can write it in the form
vx = A cos ωt + B sin ωt = a cos δ cos ωt + a sin δ sin ωt = a cos(ωt − δ), √ where a = A2 + B2 and δ = tan−1(B/A). In this form, we call δ the phase angle. Differentiating this result for vx and plugging it into the equation we got forv ˙x and solving for vy, gives us vy = −a sin(ωt − δ). All together, the velocity vector for the charged particle in the magnetic field is
~v = ha cos(ωt − δ), −a sin(ωt − δ), vz0 i .
A Complex Approach We will now solve the same problem of the charged particle in a magnetic field in a more elegant manner involving a use of complex numbers that is a standard approach in physics. Recall that
v˙x = ωvy, v˙y = −ωvx. √ If we multiply the second equation by i = −1 and add it to the first equation, we get
d (v + iv ) = ωv − ωiv . dt x y y x If we factor out −i from the right side, we get
d (v + iv ) = −iω(v + iv ). dt x y x y
If we let η = vx + ivy, we can write our first order differential equation as
η˙ = −iωη,
which has the solution η(t) = Ae−iωt. Notice that this has only a single integration constant A. Where’s the second one? The number A is a complex number so there’s a real part and an imaginary part, and that’s where the integration constants lie. Using Euler’s formula eiθ = cos θ + i sin θ, we can write
η(t) = Ae−iωt = A(cos ωt − i sin ωt).
Since the real component of eiθ is cos θ and the imaginary component is sin θ, the complex number eiθ lives on the unit circle of the complex plane. The constant A is a complex number, but we can write it in polar form as
A = aeiδ, 1.6. Charged Particle in a Magnetic Field 19 where a is the magnitude of A (i.e. its distance from the origin in the complex plane) and δ is the angle it makes with the positive x-axis. Using this definition for A, we can rewrite the solution η(t) as =(η) η(t) = Ae−iωt = aeiδe−iωt = aei(δ−ωt). η(0) = A a This shows us that η(t) lies on the circle with the same radius as A and the angle between δ A and η(t) is ωt. So at t = 0, η(t) is at A, and at some later time t = t, η(t) makes an <(η) angle ωt with A. This is all illustrated in Fig. (1.11). Notice that η(t) moves clockwise ωt around the circle of radius a as time advances. In the complex form of the solution, the constants a and δ are the same as the • constants from our previous solution. Using Euler’s formula again, we can write η(t) = aei(δ−ωt) η(t) = aei(δ−ωt) = a (cos(ωt − δ) − i sin(ωt − δ)) . Figure 1.11: The complex form Recalling that η = v + iv , we can easily substitute between the complex solution and x y of the solution on the complex our earlier solution. plane. q 2 2 Notice that a = vx + vy = v is the magnitude of the transverse velocity of the charged particle. Recall that radial velocity is the velocity toward or away from the origin when dealing with circular motion in a polar coordinate system, and the transverse velocity is the velocity component perpendicular to the radial velocity. The transverse velocity is perpendicular to B, which points in the z direction since we are currently only considering the x and y components of the velocity. To find the position function of the particle, we start by introducing another complex number ξ = x + iy, such that ξ = η(t) dt, ˆ since position is the integral of the velocity function. So that we only have to do a single integral (albeit over the complex numbers) we integrate the exponential form of η(t) rather than the expanded form containing sine and cosine. A ξ = η(t) dt = Ae−iωt dt = − e−iωt + C. ˆ ˆ iω Keep in mind that since we integrated over the complex numbers, C is a complex number. If we let D = −A/(iω) and C = X + iY , and we let the point (X,Y ) be our new origin, then our position function becomes ξ(t) = De−iωt. Because of the function being in complex polar form, we know that ξ lives on a circle with radius D in the complex plane, and at t = 0, ξ(0) = D = x0 + iy0. At some later time, t, ξ(t) has moved an angle ω clockwise around the circle. So for the position function, we have that A ia a ξ(t) = De−iωt = − e−iωt = ei(δ−ωt) = (sin(ωt − δ) + i cos(ωt − δ)) . iω ω ω The cyclotron orbit of the particle has a radius
q 2 2 a vx + vy = , ω ω where ω is the familiar angular frequency. Our results tell us that the particle moves in a circle in the complex plane, since we have x+iy, but due to the way we introduced/defined the complex function, (x, y) is also the real position of the particle in a regular rectangular coordinate system. Since x and y describe a circle as t advances, the particle really is moving in a circle. 20 Newton’s Laws
1.7 Summary: Newton’s Laws
Skills to Master • Know Newton’s laws of motion • Be able to write down the position, velocity, and acceleration vectors for a particle in Cartesian, polar, cylindrical, and spherical coordinates • Be able to derive the unit vectors of polar, cylindrical, and spherical coordinates and write them in Cartesian form • Understand the time dependence of the unit vectors in different coordinate systems • Obtain the equations of motion for a system by applying Newton’s second law in the most conve- nient coordinate system • Solve for the equations of motion of a projectile under linear air resistance • Calculate the terminal speed of an object under linear or quadratic air resistance • Solve for the equation of motion of a projectile under quadratic air resistance for the special case of 1D motion (either vertical or horizontal)
Newton’s laws of motion are: Newton’s laws imply that the time rate of change of the total momentum of a system of particles is equal First law: In the absence of forces, a body moves with to the net external force acting on the system uniform velocity. Second law: The acceleration of a body is propor- ~˙ ~ ext tional to the applied force and the proportional- P = F . ity constant is the inertial mass of the body F~ = m~a. To identify the equation of motion for a particle, use this procedure: Third law: For every action, there’s an equal but op- 1. Start with the general form of Newton’s second posite reaction law F~ = m~a 2. Identify the forces acting on the particle ~ ~ F ij = −F ji. 3. Decide which coordinate system to use and write down the acceleration ~a for that coordinate sys- Conservation of momentum is implied by the tem third law. 4. Equate component-wise the net forces from step Newton’s second law is a second order differential 2 and m~a from step 3 in the form F~ = m~a equation in the position of the moving body 5. Simplify by applying any constraints such as r˙ =r ¨ = 0 ~ ¨ F = m~r. 6. Solve for the coordinate components, integrating This ODE is called the “equation of motion,” and its as needed solution is the position function ~r(t) of the moving body. The total momentum of a system of particles is simply the sum of the momenta of the individual par- Newton’s Laws in Cartesian Coordinates ticles ~ X In Cartesian coordinates, the unit vectors xˆ, yˆ, and zˆ P = ~pi. i do not depend on time, and the position, velocity, and 1.7. Summary: Newton’s Laws 21 acceleration vectors are Newton’s Laws in Cylindrical Coordinates The transformations between Cartesian and cylindrical ~r = xxˆ + yyˆ + zzˆ coordinates are ~r˙ =x ˙xˆ +y ˙yˆ +z ˙zˆ x = ρ cos φ, y = ρ sin φ, z = z ¨ ~r =x ¨xˆ +xy ¨ yˆ +xz ¨ zˆ. p y ρ = x2 + y2, φ = tan−1 . x Then Newton’s second law, aka the equation of motion, In terms of Cartesian coordinates, the unit vectors in separates as cylindrical coordinates are F = mx¨ ρˆ = hcos φ, sin φ, 0i x ~ ˆ F = m~a =⇒ Fy = my¨ φ = h− sin φ, cos φ, 0i Fz = mz¨ zˆ = h0, 0, 1i. The position, velocity, and acceleration of a particle in These can be integrated to get the velocity functions cylindrical coordinates is x˙(t),y ˙(t), andz ˙(t) and integrated a second time to get the position functions x(t), y(t), and z(t). ~r = ρρˆ + zzˆ ~r˙ = ~v = ρρˆ + ρφ˙φˆ +z ˙zˆ Newton’s Laws in Polar Coordinates ~r¨ = ~a = ρ¨ − ρφ˙2 ρˆ + ρφ¨ + 2ρ ˙φ˙ φˆ + (¨z) zˆ. The transformations to switch back and forth between Cartesian and polar coordinates are Air Resistance Empirically, the magnitude of the force f~(v) of air re- x = r cos φ, y = r sin φ sistance is p y 2 r = x2 + y2, φ = tan−1 . f(v) = bv + cv , x where flin = bv is the linear air resistance and fquad = In polar coordinates, the unit vectors are cv2 is the quadratic air resistance.
rˆ = hcos φ, sin φ, 0i Linear Air Resistance φˆ = h− sin φ, cos φ, 0i. For a spherical object, the linear air resistance is ~ The position, velocity, and acceleration can be written f lin = −b~v. as with b = 3πηD, ~r = r rˆ and D is the diameter of the sphere and η is the vis- ~v = ~r˙ =r ˙ rˆ + rφ˙ φˆ cosity of the fluid. For a projectile in the presence of linear air resis- ¨ ˙2 ˙ ¨ ˆ ~a = ~r = r¨ − rφ rˆ + 2r ˙φ + rφ φ. tance, the net force acting on it is F~ = m~g − b~v. In polar coordinates, Newton’s second law, separates as This can be separated into Cartesian components ˙2 Fr = m r¨ − rφ Fx = −bvx F~ = m~a =⇒ ˙ ¨ Fφ = m 2r ˙φ + rφ Fy = −mg − bvy. 22 Newton’s Laws
For the horizontal component, we can separate which gives us the pair of nonlinear, coupled ODES variables and integrate to get c q t v˙ = − v2 + v2v − τ x x y x vx(t) = vx0 e , m c q v˙ = − v2 + v2v − g. where τ = 1/k = m/b is a characteristic time. Inte- y m x y y grating again gives us the position function These can be solved in the 1D cases of horizontal mo- t − τ x(t) = x0 + vx0 τ 1 − e . tion only or vertical motion only. For horizontal motion, we can integrate to get For the vertical component, we can again separate variables and integrate to get v0 v(t) = t , 1 + τ g g − t v (t) = − + + v e τ . y k k y0 where τ = m/cv0. Integrating again gives us the posi- Integrating gives us the position function tion function
− t y(t) = −v t − τ (v + v ) e τ − 1 . t ter ter y0 x(t) = x + v τ ln 1 + . 0 0 τ Terminal velocity occurs when the downward weight force equals the upward drag force. For linear For vertical motion, we have to treat the upward air resistance, it is and downward cases separately. For an object travel- mg ing downward under quadratic air resistance, Newton’s v = . ter b law is dv F = m = cv2 − mg. For an object falling in a dense fluid, we have to dt include the buoyant force. The buoyant force on such an object is equal to the weight of the displaced fluid For an object dropped from rest, this has solution Fb = ρV g, gt vd(t) = −vter tanh . vter where ρ is the density of the fluid, V is the volume of the object, and g is the local gravitational acceleration. Integrating gives us the position function The buoyant force is directed opposite the gravitational force. The net force on an object falling vertically in a 2 vter gt dense fluid and experiencing only linear drag, is yd(t) = y0 − ln cosh . g vter
dvy Fy = m = 3πηDvy + ρV g − mg. Terminal velocity for an object falling under dt quadratic air resistance is Quadratic Air Resistance rmg v = . Quadratic air resistance is given by ter c f~ = −cv2vˆ = −cv~v. quad For an object traveling upward under quadratic Newton’s second law with quadratic air resistance is air resistance, Newton’s law is d~v dv F = m = m~g − cv~v, F = m = −cv2 − mg. dt dt Chapter 2
Momentum and Center of Mass
In an earlier section, we discussed the conservation of momentum. The momentum of the ith particle is ~pi = mi~vi, and the total momentum of a system of particles is
n ~ X P = ~pi. i=1 From Newton’s third law which states that every action has an equal but opposite reaction, we know that the time derivative of the total momentum is equal to the external force on the system (i.e. the internal forces cancel each other out). If the external force is zero, then the time rate of change of the total momentum is zero, so the total momentum is constant X P~ = mi~vi = constant. i This conservation of momentum greatly simplifies the solutions to many probes. Consider a perfectly inelastic collision between a particle with mass m1 moving at ~v1 and a particle with mass m2 moving at ~v2. Because the particles stick together after the collision, the total mass is m1 + m2, and they move with the same final velocity ~vf . Using conservation of momentum before and after the collision, we know that
m1~v1 + m2~v2 = (m1 + m2)~vf , which allows us to calculate the final velocity as
m1~v1 + m2~v2 ~vf = . m1 + m2 That we can easily calculate the velocities of the particles after the collision is remarkable given the complex chemistry of an inelastic collision. Another application is that of a person of mass mp standing on a cart with wheels and mass mc, which is itself on a frictionless surface. If the person on the cart begins to run in one direction, the cart beneath him will move in the opposite direction. If the person on the cart runs toward the opposite end of the cart with speed vp relative to the cart, then the cart will move in the opposite direction with some speed vc relative to the ground. Relative to a stationary observer, the person is moving with speed vp − vc. From conservation of momentum, we know that
mcvc = mp(vp − vc).
Solving for vc gives us the speed of the cart
mp vc = vp. mc + mp 24 Momentum and Center of Mass
As always, it is a good idea to check our solution by looking at the limiting cases. If Tip mc → 0, we see that the equation reduces to vc = vp. Since your speed relative to the The tricky part of rocket stationary observer is vp − vc, this shows that you will remain stationary with respect to problems is that the mass the stationary observer. If mc → ∞, then vc → 0, so the cart remains stationary and the depends on time. person runs off of it at a speed vp.
2.1 Rocket with no External Force
A rocket, which has nothing to push against and nothing to be pushed by, moves by expelling mass. The conservation of linear momentum causes the rocket to move forward with the same but opposite momentum as the mass that is pushed out the back. A similar situation would occur if you were standing on a frictionless surface such as ice. With nothing to push against, how do you move to shore? Your only option is to throw something, such as a shoe, in the direction opposite the shore as hard as you can. To analyze a rocket, we’ll consider the system at a time t and at a later time t + dt. We don’t care what occurred before t or how massive the rocket and its fuel were before t. At time t, the rocket (and its fuel) has mass m, and the rocket is moving forward with speed v. At time t + dt, some amount of mass has been expelled in the form of exhaust. If we say the mass at this time is m + dm, where dm is less than zero, then the mass that has been expelled is −dm. The speed of the rocket (relative to a stationary observer) is now v + dv. If the speed of the exhaust (i.e. the expelled mass) is vex relative to the rocket, then the speed of the exhaust relative to a stationary observer is vex − v − dv. At time t, the total momentum of the system is
P (t) = mv.
At time t + dt, the total momentum of the system is
P (t + dt) = mom. of rocket + mom. of exhaust
= (m + dm)(v + dv) − (−dm)(vex − v − dv) = mv + m dv + vex dm.
Notice that the double infinitesimals dm dv cancel out, but even if they hadn’t we could have neglected them. Because a single infinitesimal is arbitrarily small, two of them multiplied together is effectively zero. So the change in momentum is
dP = P (t + dt) − P (t) = mv + m dv + vex dm − mv = m dv + vex dm.
Given that there are no external forces, the total momentum is conserved, meaning the change in momentum is zero
dP = m dv + vex dm = 0.
Dividing through by dt gives us
P˙ = mv˙ + vexm˙ = 0
mv˙ = −mv˙ ex.
The force −mv˙ ex on the right side is called thrust. Since the time rate of change of the mass is negative, the thrust ends up being a positive quantity. Returning to
m dv + vex dm = 0, 2.2. Multistage Rockets 25 we can find the speed as a function of the mass by integrating.
m dv = −vex dm 1 dv = −v dm ex m v m 0 1 0 dv = −vex 0 dm ˆv0 ˆm0 m m v − v = v ln 0 . 0 ex m
Notice that vex is assumed to be constant.
2.2 Multistage Rockets
Notice from the velocity function for a rocket m v = v + v ln 0 , 0 ex m that to maximize v, we need to maximize vex and the ratio m0/m. To do that, we need a multistage rocket. To analyze a two stage rocket (where one stage is emptied of fuel and then jettisoned), we need to analyze it in parts. Let m0 be the initial mass of the rocket and v0 its initial speed. Then let m1 be the mass after the first stage’s fuel has been completely burned and v1 its speed. Then let m2 be the mass after the first stage is jettisoned. Since the first stage is basically disconnected rather than blown off, the rocket’s speed is still v1. Finally, let m3 be the mass after the second stage’s fuel has been completely burned and v2 its speed. After the first stage fuel is burned but before the first stage is jettisoned, we know from the rocket equation and the given quantities that the rocket’s speed is m0 v1 = v0 + vex ln . m1
For the second stage, we begin with an initial mass of m2, a final mass of m3, and an initial speed of v1. So the speed after the second stage fuel has been burned is m2 v2 = v1 + vex ln . m3
Plugging in our value for v1 gives us m0 m2 m0 m2 v2 = v0 + vex ln + vex ln = v0 + vex ln . m1 m3 m1 m3
Notice that if the first stage is massless, that is, m1 = m2, then we recover the form of the original equation for the single stage rocket.
2.3 Rocket with External Force
Consider a rocket near the surface of the earth where the external force on the rocket is no longer zero. We now have that
ext mv˙ = −mv˙ ex + F . In the case of vertical ascent in the presence of gravity, we have that
mv˙ = −mv˙ ex − mg, 26 Momentum and Center of Mass
where we takev ˙ to be in the positive direction. Assuming a constant burn rate, the rate of mass decrease is constant, som ˙ = −k. Then
mv˙ = kvex − mg dv k = v − g dt m ex v t 0 k 0 dv = vex − g dt ˆv0 ˆ0 m
However, we know that m depends on t according to m = m0 −kt, so we need to substitute this in to get
t k 0 m0 v − v0 = 0 vex − g dt = vex ln − gt. ˆ0 m0 − kt m0 − kt From this, we can find the height of the rocket as a function of time, y(t), by inte- grating v with respect to time, and we can find the acceleration a(t) by differentiating v with respect to time.
2.4 Center of Mass
For a group of particles mi with positions ~ri, the center of mass is
1 X R~ = m ~r , cm M i i i
P where M = i mi is the total mass of the particles. The center of mass can be decoupled into the components
1 X 1 X 1 X X = m x ,Y = m y ,Z = m z , cm M i i cm M i i cm M i i i i i
then R~ cm = hXcm,Ycm,Zcmi. ~ P ˙ Recall that total momentum is given by P = i mi~ri, which tells us that the total momentum of a system of particles is the total mass times the velocity of the center of mass ˙ MR~ cm = P~ . Also recall that the time derivative of the total momentum is the total external force ~˙ P ~ ext ~ ext P = i F i = F , which tells us that the total external force acting on a system of particles is equal to the total mass times the acceleration of the center of mass
¨ ext MR~ cm = F~ .
In other words, the center of mass of a system of particles obeys Newton’s second law F~ = m~a, as if it were a point particle located at the center of mass. Consider a projectile that is following a typical parabolic projectile path when it explodes into many pieces mi. The only external force being considered is the force of ~ ext P gravity, so F = i mi~g = M~g. From the fact that the total external force equals ¨ the total mass times the acceleration of the center of mass, we have that MR~ cm = M~g 2.4. Center of Mass 27
¨ and so R~ cm = ~g. In other words, the acceleration of the center of mass of the particles following the explosion is still just the acceleration due to gravity. This tells us that after the explosion, where the original projectile would have followed a perfect parabolic trajectory, now the center of mass of all the pieces follows that perfect parabolic trajectory. This allows us to predict the distribution of the pieces after they fall to the ground. Suppose you have two particles m1 and m2 with m1 at the origin of your coordinate system. If R~ cm is the vector from the origin to the center of mass of the two particles, and ~r is the vector from the origin to m2, then the center of mass is given by
m1~0 + m2~r m2 R~ cm = = ~r. m1 + m2 m1 + m2
So the distance between m1 and the center of mass is d1 = |R~ cm| = m2r/(m1 + m2), and the distance between m2 and the center of mass is d2 = |~r| − |R~ cm| = m1r/(m1 + m2), so d m 1 = 2 . d2 m1 We consider now a continuous mass distribution. We start by examining an elemental mass dm at ~r. To find the center of mass, we must sum over all elemental masses, each multiplied by their position vector. Our sum becomes the integral
1 R~ = ~r dm. cm M ˆ
The mass density of this mass distribution depends on the position, so we denote it ρ(~r). Recall that density is mass divided by volume, so we can rewrite the elemental mass as dm = ρ(~r) dV , where dV is the volume element. Then
1 R~ = ~rρ(~r) dV, cm M ˆ where the total mass is just
M = ρ(~r) dV. ˆ
Example 2.4.1
Find the center of mass of a hemisphere of radius a and constant density ρ lying on the xy-plane and centered at the origin. The center of mass occurs at some point R~ = hX,Y,Zi. Due to the symmetry we can separate it into components and calculate the x-component, y-component, and z-component of the center of mass separately. However, whenever you do an integral, pause and see if you really have to do the integral. In this case, we can simplify the problem to finding only the z-component of the center of mass. Because of the symmetry, we know that the x and y components are just zero. So all we have to find is 1 R = zρ dV. z M ˆ We’ll do the integral in spherical coordinates. Recall that the transformation relations are
x = r sin θ cos φ, y = r sin θ sin φ, z = r cos θ. 28 Momentum and Center of Mass
The volume element, which we should know or be able to derive quickly (via the Jacobian or a geometric picture) is
dV = r2 sin θ dθ dφ dr.
Our integral becomes
π 2π 2 a 1 3 Rz = ρr cos θ sin θ dr dθ dφ. M ˆ0 ˆ0 ˆ0 Since ρ doesn’t depend on any of the integrals, we can pull it outside. We can separate the integral into three, and since the φ integral doesn’t have anything in the integrand, it simplifies to 2π.
π 2 a 1 3 Rz = ρ2π cos θ sin θ dθ r dr. M ˆ0 ˆ0 Evaluating the integral with respect to θ can be done with a straightforward 1 substitution u = sin θ after which it evaluates to 2 . Evaluating the other integral and simplifying, gives us ρπa4 R = . z 4M To calculate the total mass M, we could set up another triple integral, to evaluate ρ dV , but since we’re dealing with a hemisphere with constant density, we know that´ its volume is V = 2πa3/3, and so its mass is M = ρV = 2ρπa3/3. So Rz = 3a/8. The center of mass is then at
3 R~ = 0, 0, a . cm 8
No matter how this oddly shaped object is tumbling through a vacuum, its center of mass will exactly follow the smooth parabolic path.
2.5 Angular Momentum of a Single Particle
Linear momentum doesn’t care where the origin is, but angular momentum does. When working with angular momentum, we have to specify an origin, since the angular momen- tum of a particle is about some point. For a particle with position vector ~r and linear momentum ~p, its angular momentum about the origin is
~l = ~r × ~p.
The time rate of change of the angular momentum of the particle is
d~l d = (~r × ~p) = ~r˙ × ~p + ~r × ~p˙. dt dt Notice that we can expand the derivative of the cross product using the regular product rule. The first term on the right can be expressed as ~v × m~v, and since the cross product of parallel vectors is zero, that term drops out. Since ~p˙ = F~ we can rewrite the rightmost term as ~r × ~p˙ = ~r × F~, which is equal to the torque, ~Γ, due to F~ about the origin, so
d~l = ~r × F~ = ~Γ. dt 2.6. Angular Momentum of a System of Particles 29
If there is no torque on the particle, that is, ~Γ = 0, then the angular momentum ~l is constant, which gives us conservation of angular momentum. If F~ = 0, then ~v ~ ~ ~ ~ necessarily Γ = 0 and l is conserved. If F is a central force then ~r × F = 0 since the m cross product of parallel vectors is zero. Again, ~l is conserved. A central force is just a force parallel to ~r, and one reason they are special is that they don’t affect angular momentum. F~ Consider a planet of mass m orbiting the sun with mass M. Since the center of mass is close to the center of the system, we assume that the sun is at the origin. Let ~r point ~r from the sun to the planet. The force keeping the planet in orbit is the gravitational force F~ = −GMmrˆ/r2. Since the only force on the planet is a central force, the torque is zero and angular momentum is conserved. From ~l = ~r × ~p = ~r × (m~v), we know that ~r and ~v are both perpendicular to ~l. If M we let ~l point in the z direction, then ~r and ~v, which define the planet’s orbit, both lie in the xy-plane. Figure 2.1: A planet orbiting In a short time dt, the planet moves a small distance ~v dt, and the position vector ~r the sun. sweeps out a small area dA. Recall that the area of a parallelogram defined by two vectors is the magnitude of their cross product. This means that the area of a triangle defined by two vectors is half ~v dt of the magnitude of the cross product of the two vectors. So 1 1 1 dA dA = |~r × ~v dt| = |~r × m~v| dt = |~r × ~p| dt 2 2m 2m dA 1 = ~l . ~r dt 2m Since we know that ~l is constant, we know that dA/dt is constant, which means that the O area swept out by ~r is constant with respect to time. In other words, the radius sweeps out equal areas in equal time. This is Kepler’s second law. Figure 2.2: Kepler’s second law. We could also have used polar coordinates to describe the orbit of the planet with φ being the angle between ~r and some x-axis. Then
~r = rrˆ ~r˙ =r ˙rˆ + rφ˙φˆ
So
~l = ~r × m~r˙ = (rrˆ) × m(r ˙rˆ + rφ˙φˆ) = mrr˙(rˆ × rˆ) + mr2φ˙(rˆ × φˆ) = mr2φ˙(rˆ × φˆ) = mr2φ˙ zˆ.
2.6 Angular Momentum of a System of Particles
For a system of particles, the angular momentum of the ith particle about the origin is ~ li = ~ri × ~pi. Then the total angular momentum of the system of particles about the origin is ~ X~ X L = li = ~ri × ~pi, i i and the time rate of change of the total angular momentum is
~˙ X ˙ X ˙ X X ˙ X ˙ L = ~ri × ~pi + ~ri × ~pi = ~vi × m~vi + ~ri × ~pi = ~ri × ~pi. i i i i i Notice that the first summation is zero since ~v k m~v, and the cross product of parallel vectors is zero. 30 Momentum and Center of Mass
Recall from Newton’s second law that the time rate of change of the momentum of the ith particle is the sum of the forces on i due to all the other particles plus the external force on i ˙ ~ ext X ~ ~pi = F i + F ij. j6=i Using this, we can rewrite the time rate of change of the total momentum as ~˙ X ~ ext X ~ X ~ ext X X ~ L = ~ri × F i + F ij = ~ri × F i + ~ri × F ij. i j6=i i i j6=i
~ ext Since ~ri × F i is the external torque on the ith particle, the first term is just the total external torque on the system about the origin
~ ext X ~ ext Γ = ~ri × F i . i The second term consists of pairs of terms of the form
~ri × F~ ij + ~rj × F~ ji,
where the first term is the torque on particle i due to particle j and the second term is the torque on particle j due to particle i. If the internal forces obey Newton’s third law, then F~ ji = −F~ ij and we can write the above as
~ri × F~ ij + ~rj × F~ ji = (~ri − ~rj) × F~ ij.
In a system of particles where the internal forces are all central forces and Newton’s third law holds, we have that (~ri − ~rj) × F~ ij = 0, that is, all the internal torques vanish since ~ri − ~rj is parallel to F~ ij. This would not be the case if the internal forces were not central forces. For example, we could have a system of particles in which the forces between pairs of particles are equal and opposite but not pointing toward each other. In such a case, there would be internal torque and the angular momentum would increase. In our case, the internal torques cancel, so we have that the time rate of change of the total momentum of the system is equal to the total external torque on the system.
˙ ext L~ = ~Γ .
If the external torque is zero, then the angular momentum of the system is conserved.
2.7 Angular Momentum of a Continuous Mass Distribution
Consider a rigid, extended body rotating about the z axis with angular velocity ω. We can think of the rigid body as a system of particles where each particle is frozen into a configuration such that their distances from each other and the angles between them are constant. If we consider the system of rigid particles in cylindrical coordinates, then the ith particle has coordinates (ρi, φi, zi). In rigid rotation about the z axis, ρi and zi are ˙ constant. Only the φi component changes with time. Thenρ ˙i = 0,z ˙i = 0, and φi = ω, and these are the same for all particles. So this is a very restricted motion compared to, say, a gas of particles in which the particles have the freedom to move in all directions. The position of the ith particle is given by the vector
~ri = ρiρˆ + zizˆ. 2.7. Angular Momentum of a Continuous Mass Distribution 31
Differentiating, gives us the velocity vector Tip Many problems can be ~r˙ =ρ ˙ ρˆ + ρ ρˆ˙ +z ˙ zˆ =ρ ˙ ρˆ + ρ φ˙ φˆ +z ˙ zˆ = 0 + ρ ωφˆ + 0 = ρ ωφˆ. i i i i i i i i i i i solved using conservation of angular momentum of Then the angular momentum of the ith particle about the origin ~l = ~r × (m ~r˙ ) can be i i i a system when the exter- written as nal torque on the system is h i zero. Furthermore, the an- ~l = (ρ ρˆ + z zˆ) × (m ρ ωφˆ) = m ρ ω ρ (ρˆ × φˆ) + z (zˆ × φˆ) = m ρ2ωzˆ − m ρ ωz ρˆ. i i i i i i i i i i i i i i gular momentum of a rigid and symmetric body can be The cross products of the unit vectors are obtained using the right-hand rule. We’re only found easily by first calcu- concerned with the z-component of the angular momentum so lating its moment of inertia, then its angular momentum 2 liz = miρi ω, about its axis of rotation is L = Iω. then the z-component of the total angular momentum of the system of rigid particles is
X X L = l = m ρ2ω = Iω, z iz i i Tip i i For the moment of inertia, where I = ρ2 dm you must keep X 2 I = miρi , in mind´ that ρ is the dis- i tance of an elemental mass from the axis of rotation— is called the moment of inertia of the rigid body. The moment of inertia doesn’t depend not from the origin. This on the angular speed ω at which the body is rotating, it is only a parameter associated makes calculating the mo- with the mass distribution of the rigid body. The moment of inertia makes it easy to ment of inertia in spherical calculate the angular momentum of a rigid and symmetric body rotating about a given coordinates, for example, a axis little more tricky. Lz = Iω.
Depending on the shape and mass distribution of a body, a body rotating about the z axis will have angular momentum components in the x and y directions, but we are only considering the component in the z direction. For the continuous case, we consider an elemental mass dm instead of a particle mi, and our summation turns into an integral
I = ρ2 dm = ρ2 ρ dV, ˆ ˆ dens where ρ is the distance from the elemental volume to the axis of rotation, and ρdens is the mass density of the rigid body. Example 2.7.1
Calculate the moment of inertia of a ring with mass M and radius a. This should be immediately obvious as
I = Ma2,
since the moment of inertia is just the distance squared of an elemental mass from the axis of rotation times its mass, summed over all elemental masses. Since all the elemental masses are the same distance, a, from the axis of rotation, I is just Ma2. 32 Momentum and Center of Mass
Example 2.7.2 mi 0 ~ri Calculate the moment of inertia of a uniform disk with mass M and radius a. C.O.M. We know that ~ri I = r2 dm, ˆ R~ where r is the distance from an elemental mass dm to the axis of rotation. We know that the surface density of the disk is σ = M/A where the total mass is M and the surface area is A, so dm = σ dA O
I = σr2 dA. ˆ To convert to 2-dimensional polar coordinates, we make the substitution dA = r dr dθ, then our integral becomes
2π a a 1 I = σr3 dr dθ = σ2π r3 dr = σπa4. ˆ0 ˆ0 ˆ0 2 Since σ = M/A = M/(πa2), we have that 1 I = Ma2. 2
What if you want to know the angular momentum of a system about its center of mass when the center of mass is not at the origin of your coordinate system? Recall that the total linear momentum of a collection of particles is the total mass times the velocity of the center of mass. Is there an analogous thing for angular momentum? If we have a system of particles mi and associated positions ~ri, then their center of ~ 1 P mass is given by Rcm = M i mi~ri. We’ll start by redefining our coordinate system so that the center of mass is at the new origin, then
0 ~ri = R~ cm + ~ri ,
0 0 where ~ri are the new position vectors for the particles. That is, ~ri points from the center of mass to the ith particle. Keep in mind that the center of mass won’t necessarily be in an inertial reference frame, since it may be accelerating relative to the original origin. Differentiating gives us ˙ ˙ ˙ 0 ~ri = R~ + ~ri . The angular momentum of the system about the original origin is
~ X X ˙ L = ~ri × ~pi = ~ri × mi~ri. i i
Substituting in the equations above gives us
X 0 ˙ ˙ 0 L~ = (R~ + ~ri ) × mi(R~ + ~ri ) i X h ˙ ˙ 0 0 ˙ 0 ˙ 0 i = mi (R~ × R~ ) + (R~ × ~ri ) + (~ri × R~ ) + (~ri × ~ri ) i X ˙ X ˙ 0 X 0 ˙ X 0 ˙ 0 = mi(R~ × R~ ) + mi(R~ × ~ri ) + mi(~ri × R~ ) + mi(~ri × ~ri ). i i i i 2.7. Angular Momentum of a Continuous Mass Distribution 33
P 0 P ˙ 0 The middle two sums drop out because i mi~ri = 0, and so i mi~ri is also zero. P 0 The sum i mi~ri gives the center of mass of the particles with respect to the center of 0 mass since it uses the ~ri vectors, and so the sum is zero. Our equation simplifies to
˙ X 0 ˙ 0 L~ = R~ cm × MR~ cm + ~ri × mi~ri . i Notice that the first term on the right is the angular momentum associated with the motion of the center of mass itself, and the second term is the total angular momentum of the system about the center of mass. Solving for the second term on the right, we have that L~ cm = L~ − R~ cm × P~ , so the angular momentum of a system about the center of mass is the angular momentum of the system about the origin minus the angular momentum of the center of mass about the origin. For example, if the sun is taken to be the origin of our coordinate system, but we want to know the angular momentum of Earth about its own center of mass, then we could calculate the angular momentum of the particles of Earth about the sun and subtract from that the angular momentum of Earth’s center of mass about the sun. ~˙ ~ ext P ~ ext ~˙ ~ ext P ~ ext Recall that P = F = i F i and that L = Γ = i Γi . Differentiating the equation above, gives us
d ˙ ˙ ˙ ext ext L~ = L~ − R~ × P~ − R~ × P~ = ~Γ − R~ × F~ dt cm cm cm cm X ~ ext X ~ ~ ext = ~ri × F − Rcm × F i i i X ~ ~ ext X 0 ~ ext ~ ext = (~ri − Rcm) × F = ~ri × F = Γcm . i i So the time rate of change of the angular momentum of a system about its center of mass is equal to the external torque on the system about the center of mass. This result holds even if the center of mass is accelerating. Recall that Newton’s laws are only valid in inertial frames. 34 Momentum and Center of Mass
2.8 Summary: Momentum and Center of Mass
Skills to Master • Use conservation of linear and angular momentum to solve problems • Calculate the center of mass of objects and systems of objects • Apply Newton’s laws to systems of objects • Calculate the torque acting on an object • Be able to derive Kepler’s second law from angular momentum conservation • Calculate the moment of inertia for various mass distributions
Linear Momentum The time rate of change of the angular momen- tum of a system about its center of mass is equal to The total momentum of a system of n particles is the external torque on the system about the center of n mass X P~ = ~p . d ~ ~ ext i Lcm = Γcm . i=1 dt The center of mass of a continuous mass distribu- If the external force is zero, then the time rate of change tion with density ρ(~r) is of the total momentum is zero, so the total momentum 1 1 is constant R~ = ~r dm = ~rρ(~r) dV, cm M ˆ M ˆ ~ X P = mi~vi = constant. where the total mass is i M = dm = ρ(~r) dV. This conservation of momentum greatly simplifies the ˆ ˆ solutions to many probes. When integrating over a mass to find the center of mass we can often use symmetry and integrate in cylindrical Center of Mass or spherical coordinates to greatly simplify the prob- lem. To do that, we need to remember the volume For a group of particles mi with positions ~ri, the center of mass is elements in those coordinate systems. ~ 1 X Rcm = mi~ri, • Cylindrical coordinates (x, y, z) → (ρ, φ, z), M i dV = ρ dρ dφ dz P • Spherical coordinates (x, y, z) → (r, θ, φ), dV = where M = i mi is the total mass of the particles. The center of mass of a system of particles obeys r2 sin θ dr dθ dφ Newton’s second law as if it were a point particle at ˙ R~ cm. That is, MR~ cm = P~ , where P~ is the total mo- Angular Momentum ~¨ ~ ext ~ ext mentum and MRcm = F where F is the net The angular momentum of a particle about the origin external force. is The angular momentum of a system about its cen- ~l = ~r × ~p. ter of mass is the angular momentum of the system The torque on a particle is the time rate of change of about the origin minus the angular momentum of the its angular momentum center of mass about the origin d~l ~ ~ ~ ~ ~Γ = = ~r × F~. Lcm = L − R × P . dt 2.8. Summary: Momentum and Center of Mass 35
If F~ is a central force, then it applies no torque since where ρ is the distance to the elemental mass from ~r × F~ = 0. If either F~ = 0 or ~Γ = 0, then angular the axis of rotation (not the origin) and ρdens is the momentum is conserved. density. Some common moments of inertia are The total angular momentum of a system of parti- cles is just the sum of the individual angular momenta. • Point mass at distance r: I = mr2 The time rate of change of the total angular momen- • Hoop or cylindrical shell of radius R: I = mR2 tum, provided that all internal forces are central forces, • Rod of length L: I = mR2/12 is just the net external torque on the system • Disk or solid cylinder of radius R: I = mR2/2 • Sphere of radius R: I = 2mR2/5 ˙ ext L~ = ~Γ . Recall that the moment of inertia is the rotational If there’s no net external torque on the system, then analogue of mass and that the rotational analogue of its angular momentum is conserved. Newton’s second law is For a rigid body rotating about the z axis, the z-component of the angular momentum is Γ = Iα,
Lz = Iω, where Γ is the torque and α is the angular accelera- tion. Likewise, L = Iω is the rotational analogue of ˙ z where ω = φ is the angular velocity of rotation and p = mv.
X 2 When solving problems, determine if angular mo- I = miρi . mentum is conserved. If it is, then you know that, i ~ ~ For a continuous mass distribution, Li = Lf
~ri × ~pi = ~rf × ~pf 2 2 I = ρ dm = ρ ρdens dV, ˆ ˆ ripi sin θi = rf pf sin θf . Chapter 3
Energy
For a particle of mass m traveling at speed v, the kinetic energy is
1 T = mv2. 2
Notice that v2 = ~v · ~v. Substituting this allows us to take the time derivative of the kinetic energy as dT 1 = m ~v˙ · ~v + ~v · ~v˙ = m~v˙ · ~v. dt 2 From Newton’s second law, we know that m~v˙ = F~, and ~v dt = d~r is just the infinitesimal displacement along the particle’s path, which gives us dT = F~ · d~r. This is the work done by the force F~ in the displacement d~r, and we write dW = F~ · d~r. So we can think of work as the differential of the kinetic energy with respect to displace- ment. For a particle moving from point 1 to point 2 along a path, the total work done by the force is calculated by summing over all the differential d~r’s
2 W1→2 = F~· d~r. ˆ1 Since we are summing over a path that is not necessarily a straight line, the integral above is a line integral. Earlier, we showed that dT = F~ · d~r. This implies that
2 ∆T = T2 − T1 = F~· d~r. ˆ1 So the work done in moving the particle from 1 to 2 is the same as the change in kinetic energy. This is called the work-energy theorem. A force is conservative if it satisfies the following two conditions • The force depends only on position and not on velocity or time or any other pa- rameter. This excludes air resistance, for example. • The work done in moving a particle from point 1 to point 2 is independent of the path taken. This excludes friction, for example, since a longer path would obviously mean more work is done. Kinetic friction, for example, fails both of the above since it depends on velocity and the path. Energy 37
Example 3.0.1
Consider a particle being moved from point 1 to point 2 in a uniform gravita- ~ 2 ~ 2 tional field F = −mgzˆ. The work done is W1→2 = 1 F · d~r = 1 −mgzˆ· d~r. Since ´ 2 ´ zˆ·~r is just the z component of d~r, this simplifies to 1 −mg d~z = −mg(z2 − z1) = −mg∆z, so the work done depends only on the change´ in the height of the particle and not on its path. Since gravitational force depends only on position, it is a conservative force.
If we have a conservative force, we can define a potential energy for that force. To do that, we first have to define some reference point ~r0 where we define the potential energy to be zero, then the potential energy of the particle at position ~r is the negative of the work done for the force to move the particle from the reference point to the new point
~r ~ ~0 U(~r) = −W~r0→~r = − F · dr . ˆ~r0
Potential energy and work are related in that if a conservative force does work on a particle, the result is an increase in potential energy, and if the force does negative work, the result is a decrease in potential energy. A simple example is that of dropping an object. The force of gravity does negative work when the object drops and there is a resulting decrease in the object’s potential energy. We know that for a particle being taken from point 1 to point 2, the change in kinetic energy is related to the work done via T2 − T1 = W1→2. But for a conservative force, this equals the work done in taking the particle from point 1 to a reference point 0 plus the work done taking the particle from point 0 to point 2, so T2 − T1 = W1→0 + W0→2. The work done in taking the particle from point 0 to point 1 is the negative of the work done taking the particle from point 1 to point 0, so we can write T2 − T1 = −W0→1 + W0→2. Now on the right side we have two potential energy functions—the potential energy at point 1 relative to the reference point and the potential energy of the particle at point 2 relative to the reference point, so we can write T2 − T1 = U1 − U2. This gives is the important result that for conservative forces, mechanical energy, the sum of the kinetic and potential energy, is conserved
E = T + U.
We have related the potential energy function (a scalar function) to the integral of a vector function F~. Now, we want a way to find the force F~ given the potential energy function U(~r). We know that dW = F~ · d~r = Fx dx + Fy dy + Fz dz. For a conservative force, we also know that dW = −dU. We can write the differential dU as dU = ∇~ U · d~r. This tells us that dW = F~ · d~r = −dU = −∇~ U · d~r, which implies that
F~ = −∇~ U.
That is, the force on a particle is the negative of the gradient of the particle’s potential energy function. Example 3.0.2
Calculate the force if the potential energy function is
U(~r) = mgz. 38 Energy
Since there is only a z and no x or y, the gradient is simply
∇~ U = mg zˆ,
so the force is F~ = −∇~ U = −mg zˆ.
Example 3.0.3
Using Cartesian coordinates, find the force on a particle if its potential energy varies as the inverse of r. We have that α α U(~r) = = , r px2 + y2 + z2 where α is an arbitrary constant. We know that ∂U ∂U ∂U F~ = −∇~ U = − xˆ − yˆ − zˆ. ∂x ∂y ∂z Taking the partial derivative with respect to x, we have that
∂U 1 3 αx 2 2 2 − 2 = − α(x + y + z ) (2x) = − 3 . ∂x 2 (x2 + y2 + z2) 2
Notice that we can rewrite this as ∂U αx = − . ∂x r3 The other two partials have the same form, so αx αy αz α α α F~ = xˆ + yˆ + zˆ = (xxˆ + yyˆ + zzˆ) = ~r = rˆ. r3 r3 r3 r3 r3 r2 Gravity is an example of this kind of force. The potential function varies with inverse distance and the force varies with inverse distance squared.
Example 3.0.4
Using spherical coordinates, find the force on a particle if its potential energy varies as the inverse of r. We have that α U(~r) = , r where α is an arbitrary constant. We know that ∂U 1 ∂U 1 ∂U F~ = −∇~ U = − rˆ − θˆ − φˆ. ∂r r ∂θ r sin θ ∂φ Since our potential function is only a function of r, we only have to worry about the first term in the gradient. The rest are zero. ∂U α ∇~ U = = − rˆ. ∂r r2 So α F~ = −∇~ U = rˆ. r2 Energy 39
Given a force, can we definitely determine if it is conservative or not? The answer is yes, and it derives from Stokes’ theorem. A force F~ is conservative if and only if its curl is zero ∇~ × F~ = ~0.
Example 3.0.5
Determine if F~ = hay, bx, 0i is a conservative force. Taking the curl of F~, we find that
∇~ × F~ = (b − a) zˆ,
so the force is conservative if a = b.
Once we know that a force is conservative, it’s straightforward to do the line integral
~r U(~r) = − F~· dr~0. ˆ~r0 Since we know the force is conservative, we can pick any convenient path to integrate over. We can also let the reference point be the origin ~r0 = h0, 0i, and let the final point be an arbitrary position ~r = hx, yi. Then the most convenient path would probably be the path from the origin along the x-axis to (x, 0) and from there, vertically to (x, y). Then we can easily split the line integral into a pair of regular integrals—one with respect to x and the other with respect to y. Since y = 0 along the entire path integrated with respect to x, this integral will be 0 if Fx contains a factor of y. Example 3.0.6
Find the potential function for
F~ = hay, ax, 0i.
Earlier, we confirmed that this force is conservative, so we know that we can find a unique potential function
~r U(~r) = − F~· dr~0. ˆ~r0 We let our reference point be the origin and the final point be (x, y). Then we integrate on a path from the origin to (x, 0), and from there to (x, y), so we can split the integral into
x y x y 0 0 0 0 U = − Fx dx − Fy dy = − ay dx − ax dy . ˆ0 ˆ0 ˆ0 ˆ0 Since y = 0 along the entire path integrated with respect to x, the first integral is zero, and we have y U = − ax dy0 = −axy. ˆ0 To confirm, we check that −∇~ U = F~, and it does.
Any force that depends only on x and not on the velocityx ˙ is automatically conser- vative because there’s only one path. If a force depends only on position, then the work done in moving the particle from point 1 to point 2 is the same even if the particle takes 40 Energy
x2 0 0 a curvy path. The work done in that case is W1→2 = F (x ) dx and the potential x1 Tip x 0 0 energy is U(x) = − F (x ) dx , where x0 is the reference´ point. x0 “One-dimensional” in me- ´ chanics generally means 3.1 General One-dimensional Systems that a system can be completely characterized A second case in which the force is guaranteed to be conservative is in a one-dimensional by a single parameter. system. In this case, 1-dimensional, doesn’t necessarily mean that the particle is moving It doesn’t mean the mo- in a single spatial dimension, but that the particle’s position can be described using a tion is in a single spatial single parameter. For example, a bead sliding on a curvy wire is a 1D system since its dimension. position along the wire can be specified by a single number—the length of wire between it and where it started. Of course there is a normal force exerted by the wire on the bead, but since this force is always perpendicular to the bead’s motion, it does no work. It is U(x) just a constraint force. E Another example of a 1-dimensional system is an Atwood machine with a frictionless and massless pulley. In this case, the tension force is a constraint force and on both strings it does work, but since they cancel, there is no net work done. This is true in general—constraint forces do no net work. x In general, for a one dimensional system (where the force depends only on the par- ticle’s position x) in which the forces are all conservative, the potential energy is x U(x) = − F (x0) dx0, Figure 3.1: A plot of the po- x ˆx0 tential and total energies for a generic 1D system. the force is dU F = − , x dx Tip and the total mechanical energy is conserved and 1 It is often helpful to con- E = mx˙ 2 + U(x) = constant. sider extreme examples 2 when trying to intuit the This is a separable differential. After solving forx ˙ 2 and taking the square root of both behavior of a specific case. sides, we can separate variables and solve: For example, what happens if a box in equilibrium r 2 x˙ = ± pE − U(x) on top of a cylinder is m nudged. We can consider rm 1 the extreme case where dt = ± dx 2 p the box is much smaller E − U(x) r x than the cylinder, and the m 1 t(x) = ± dx0. other extreme where the p 0 2 ˆx0 E − U(x ) box is much larger than the cylinder. Clearly, the The sign on the right side is fixed by the initial velocity of the particle. Is it moving left behavior depends on their or right? The initial position x0 is x(t = 0). To find x(t), the position as a function of relative sizes. t, we only have to invert the result found above. To find the velocity of the particle we would then differentiate x(t). It is often helpful to plot U(x) and the total energy E on the same graph as shown in Fig. (3.1). From a graph like this, there are several things we can say qualitatively about the motion of the object. From the conservation of mechanical energy, we know that the particle will have the largest velocity when U(x) is at its least. We know it will have the lowest velocity when U(x) is at its greatest. This graph shows the potential energy of an oscillator. If the particle is heading to the right of x = 0 its velocity will be decreasing since the potential energy is increasing. Eventually it will reach a stopping point and then accelerate back toward x = 0. Since it is impossible for the potential energy U(x) to be greater than the total energy E, we know that regions on either side of the intersections of U(x) and E are 3.1. General One-dimensional Systems 41 forbidden zones as shown in Fig. (3.2). These intersections are turning points where the U(x) velocity of the particle reaches exactly zero and changes sign. Clearly, this value of x is the amplitude of the oscillator. Consider the potential energy function shown in Fig. (3.3). In Fig. (3.4) is the graph with the turning points noted and the forbidden regions identified. Notice that with this potential energy, a particle with initial position within the first allowed region will oscillate with a large amplitude in that region, and a particle with initial position within the second region will oscillate with small amplitude in that region alone. The behavior −A A of the system is different based on the initial conditions. Figure 3.2: Forbidden regions. dU Recall that the force on the particle is the negative of dx , so we can infer just by looking at any point on the graph, which direction the force is pointing and what its magnitude is. At any point x on the graph, the force is the negative of the first derivative U(x) dU of U(x). This means for points where dx = 0, such as at local minimums and maximums, E the force is zero. A particle placed exactly at any of those points will just sit there at rest. We call these equilibrium points. Notice that there are two kinds of equilibria. x At a local minimum, if a particle is nudged in either direction, the force on it will point d2U back toward the equilibrium. So when dx2 > 0 at that point, it is a stable equilibrium point. At a local maximum, if a particle is nudged in either direction, the force on it will point away from the equilibrium and the particle will fly away from that point. So when d2U Figure 3.3: Potential and total dx2 < 0 at that point, it is an unstable equilibrium point. d2U energy. If we get dx2 = 0, this doesn’t tell us whether the equilibrium is stable or not. In such a case, we may be able to expand U(x) as a Taylor series, differentiate it, then examine the higher terms.
Example 3.1.1
Given the potential function U(x) = 2 cos x + x sin x, determine the stability at x = 0. The first derivative is U 0 = − sin x + x cos x. Evaluating the first derivative at x = 0 shows that U(x) has an equilibrium point there. To know if it’s a stable or unstable equilibrium point, we evaluate the second derivative at that point. The Figure 3.4: Turning points and second derivative is U 00 = −x sin x. This is 0 at x = 0, the second derivative test forbidden zones. doesn’t tell us if this is a stable or unstable equilibrium point. To answer this, we can expand U(x) about the point 0 to get U(x) x2 x4 x3 U(x) = 2 cos x + x sin x = 2 1 − + − · · · + x x − + ··· E1 2 24 6 x x4 = 2 − + O(x6). 12 E2 From this, we can see that the second derivative of U(x) near x = 0 is negative as zero is approached from either direction since U 00 = −x2 + O(x4). So the equilibrium point is unstable although barely since for x near x = 0, the quantity Figure 3.5: Electron potential x2 will be very small. well.
The plot in Fig. (3.5) represents a more interesting potential energy function. If the particle’s energy is E2 then it is constrained to oscillate in a small allowed region. If its energy is E1 it becomes unbounded and escapes to infinity. Think of the behavior of an electron near its host atom at x = 0. Notice that the electron can never actually reach the nucleus because that region is forbidden. If its energy is low, it is bound to its host atom and will oscillate close to it indefinitely. If its energy is large, however, it will escape this potential well and leave the atom. 42 Energy
Example 3.1.2
Consider an object in free fall, disregarding air resistance. We know that the gravitational force is conservative, and the object’s potential energy function is
U(y) = −mgy,
where +y is taken to be the downward direction. If the mass is dropped from rest, and we take the reference point for the potential energy function to be the y position where the object was dropped from, then the object’s total energy is 0. This is because its potential energy at y = 0 is taken to be zero and its kinetic energy is zero at the moment it is dropped from rest. Then using the formula we derived earlier
rm y 1 r 1 y 1 r 1 y r2y t(y) = √ dy0 = √ dy0 = 2py0 = . 0 0 2 ˆ0 mgy 2g ˆ0 y 2g 0 g Inverting this gives us the familiar result 1 y(t) = gt2. 2
Example 3.1.3
Analyze the energy characteristics of a harmonic oscillator, such as a spring, obeying Hooke’s law. Hooke’s law gives us the force on the object
F (x) = −kx.
Notice that this is a first-order (i.e. linear) restoring (notice the negative sign) force. From this, we should already be able to picture the potential energy func- tion. If the object is at its equilibrium position then a nudge in either direction will make it want to return to equilibrium, so the potential energy function should be an upward opening graph with a local minimum at the equilibrium position x = 0. Calculating the potential energy function, we have that
x x 0 0 0 1 2 U(x) = − Fx(x ) dx = k x dx = kx . ˆx0 ˆx0 2 We see that it is an upward opening quadratic centered at the equilibrium position as shown in Fig. (3.6). For a given energy E, we can determine the amplitude A of the object’s oscillations by seeing where E intersects U(x) as shown in Fig. (3.7). When the object reaches x = A, its velocity will be zero, so its kinetic energy will be 0, and so its total energy must equal its potential energy at x = A, that 1 2 is, E = 2 kA . We can now plug E and U(x) into the formula we derived earlier to get its equation of motion
rm x 1 rm x 1 t(x) = ± dx0 = ± √ dx0. p 0 2 02 2 ˆx0 E − U(x ) k ˆx0 A − x
Notice that the time it takes the object to go from x = 0 to x = A can be calculated by integrating this integral from x0 to A. The period (i.e. one cycle), τ, of the oscillator is then 4 times this quantity. In one cycle, the object goes from x = 0 3.2. Atwood Machine 43
to x = A to x = −A and back to x = 0. So U(x)
rm A 1 τ = 4 √ dx0. 2 02 k ˆ0 A − x Making the change of variables u = x0/A lets us express the integral in dimension- x less variables rm 1 1 τ = 4 √ du. 2 k ˆ0 1 − u Figure 3.6: Harmonic oscillator From this, we can see that the period of a harmonic oscillator does not depend potential well. on its amplitude. Noting that the integral above is the arcsin, we can go on to evaluate it to get rm U(x) τ = 2π . k E
If you have an integral, it is always a good idea to try to write it in terms of dimen- sionless variables. This allows you to see the relations between the different variables. In the example above, we made the substitution u = x/A and since both x and A have units of length, the result u is dimensionless. This allowed us to determine that the period of x a harmonic oscillator does not depend on the amplitude without us having to actually evaluate the integral. In general, try to extract dimensionful variables from integrals by making a natural change of variables. Usually then, we don’t even need to know the ac- −A A Figure 3.7: Harmonic oscillator tual value of the integral—it is just some number. What we care about is the physically potential well. interesting relations between the dimensionful variables. This cannot be overemphasized.
3.2 Atwood Machine
Here, we analyze an ideal (frictionless, massless pulley and massless string) Atwood ma- chine using an energy approach. Calculating the potential energy of an Atwood machine can seem tricky given that one mass will be moving down while the other mass will be moving up. However, recall that potential energy is just the negative of the work done, and the work being done is clear. We can ignore the work done by the tension forces since the net work done by the tensions is zero. We can ignore the tension force from the beginning because we’re dealing with a 1-dimensional system, and the tension is just a constraint force. All we have to calculate is the work done by the gravitational force. In the case of an Atwood machine, the work done can be calculated as just the gravitational force times the distance since the angle between the force and the direction of motion is zero. If block M moves downward a distance H, the work done against the gravitational force is −MgH since the object is falling. At the same time, block m moves up a distance H, so the work done on that block is mgH. The total work done is
W = −MgH + mgH = (m − M)gH, and in general, if M moves down a distance x, then
W = (m − M)g(H − x).
Since potential energy is the negative of the work, we have that
U(x) = g(M − m)(H − x).
Notice that if M = m, then U(x) = 0. This makes sense since if the two masses are the same, the system won’t move—it has zero potential energy. If M > m, then in the 44 Energy
Tip first Atwood machine in Fig. (3.8), we have that U(x) = g(M − m)H. This makes sense since the potential energy for an object in a gravitational field is just U = mgH, and in For 1 dimensional systems, this case, the small mass m effectively reduces the total mass. In other words, this is the you can extract an equa- potential energy of a single object at a height H with mass M − m. In the third Atwood tion of motion by differ- machine shown in Fig. (3.8), we have that U(x) = g(M − m)(H − x) = 0 since x = H. entiating the total energy Since we’ve defined M > m, we know that with M resting on the floor, the system is with respect to time. When at rest—its lowest potential energy. The same equations still work if m > M, the only doing this, be careful to difference is that the potential energy instead of going from some positive maximum to use the product rule where zero will go from a maximum at zero to a negative value. you need to... don’t make We can also define the potential energy of the Atwood machine in terms of the length the mistake of differentiat- of the string using the pulley as the reference point rather than the floor. In this case, ing with respect to position the potential energy can be calculated as −Mgx − mg(L − x) where L is the length of when you need to differen- the string. The signs are negative since we’re measuring downward rather than upward. tiate with respect to time. The total energy of the Atwood machine is just the sum of the kinetic and potential Previously, we only knew energies 1 1 how to generate an equa- E = Mx˙ 2 + mx˙ 2 − Mgx − mg(L − x). tion of motion by looking at 2 2 the net forces on an object We know that mechanical energy is conserved, that is, E = constant, so we can differen- and applying Newton’s sec- tiate the energy with respect to time to get us the equation of motion ond law to it. dE 1 1 = M(2x ˙x¨) + m(2x ˙x¨) − Mgx˙ + mgx˙ = 0 dt 2 2 0 = (M + m)x ˙x¨ − (M − m)gx˙ 0 = (M + m)¨x − (M − m)g M − m x¨ = g. M m M + m x Notice that when differentiatingx ˙ 2 with respect to time, we had to use the product rule. If the pulley is massive, the only difference is that you must include the pulley’s H M m rotational kinetic energy Iω2/2 in the total energy.
x m M 3.3 Spherically Symmetric Central Forces Another whole class of forces that are automatically conservative are spherically-symmetric Figure 3.8: Atwood machine. central forces. In spherical coordinates, a general potential energy function has the form U(~r) = U(r, θ, φ), that is, it is a function of r, θ, and φ. A spherically symmetric central force is the special case in which U depends only on the distance r from the origin and not on either of the angles. That is, U = U(r). In this case, the force can be calculated as F~ = −∇~ U(r). Applying the gradient operator for spherical coordinates, we get
∂U F~ = − rˆ, ∂r
∂U ∂U since ∂θ = ∂φ = 0. We can see that it is a central force since it has a component only in the r direction. Any spherically symmetric central force has the form
∂U F~ = f(r) rˆ, f(r) = − . ∂r A force is a central force if it can be written in the form f(~r). It is, further, spherically symmetric if it can be written in the form f(r). A central force is conservative if and only if it is spherically symmetric. 3.4. The Energy of a System of Particles 45
Given U(r), we can calculate F~ using the equation above. Given F~ = f(r) rˆ, we can calculate U(r) as r U(r) = − f(r0) dr0. ˆr0
3.4 The Energy of a System of Particles M m For a system of particles, we can calculate the total kinetic energy by summing the indi- vidual kinetic energies. If the external force acting on a system of particles is conservative, H − x x that is, it depends only on the position of a particle, then the external force on particle i is related to the potential energy via
~ ext ~ ext Figure 3.9: Atwood machine. F (~ri) = −∇iUi (~ri).
Consider the case of two interacting particles. If the position of particle 1 is ~r1 and the position of particle 2 is ~r2, then the vector pointing from particle 2 to particle 1 is just ~r = ~r1 − ~r2. The force on particle 1 due to particle 2 is then a function of ~r − ~r , that is F~ = 1 2 12 x L − x F~ 12(~r1 − ~r2). In other words, the force on particle 1 due to particle 2 could depend on the position of both particles, but only in the form ~r − ~r . Since the force between two 1 2 m particles is the same regardless of their position (provided that their relative positions stay M the same), we can conveniently translate the system so that ~r coincides with the origin 2 Figure 3.10: Atwood machine. of our coordinate system. This property is called translational invariance. We can now treat the system as a single particle under an external force, that is F~ 12 = F~ 12(~r1). Provided that the force between the particles is conservative, we can now define a Tip potential energy function for the pair of particles such that One way to get the equation F~ 12 = −∇~ 1U(~r1), of motion for a 1D system is to differentiate the total when particle 2 is at the origin, and energy with respect to time.
F~ 12 = −∇~ 1U(~r1 − ~r2), when not at the origin. Note that ∇~ 1 is the gradient with respect to particle 1 ∂ ∂ ∂ ∇~ 1 = xˆ + yˆ + zˆ. ∂x1 ∂y1 ∂z1 In words, the force on particle 1 due to particle 2 is the negative of the gradient (with respect to particle 1) of the potential energy function of the system. For the force on particle 2, we can say that F~ 21 = −∇~ 2U(~r1 − ~r2), that is, the force on particle 2 due to particle 1 is the gradient at particle 2’s position of the same potential energy function. From Newton’s third law, we know that F~ 12 = −F~ 21, so −∇~ 1U(~r1 − ~r2) = ∇~ 2U(~r1 − ~r2). We can now write the total potential energy for a system of two particles as
ext ext U = U1 (~r1) + U2 (~r2) + U12(~r1 − ~r2). It is the sum of the potential energy of particle 1 due to the external force, the potential energy of particle 2 due to the external force, and the internal potential energy of the system. Notice that we don’t have to include a second internal potential energy function 46 Energy
since both particles experience the same internal potential energy. To confirm, we can take the negative gradient of U at particle 1 to get the total force on particle 1
~ ~ ext ~ ~ ext ~ −∇1U = −∇1U1 (~r1) − ∇1U12(~r1 − ~r2) = F 1 + F 12,
which we know is the total force on particle 1 in the system of two particles. Notice that ~ ext the term −∇1U2 (~r2) drops out. The total potential energy for a system of n particles is
n n n X ext X X U = Ui (~ri) + Uij(~ri − ~rj), i i j>i
where the first term is the total potential energy due to external forces and the second term is the total potential energy due to the internal forces of the system. It may seem that we’re ignoring half of the potential energies in the second term, but keep in mind that Uij(~ri − ~rj) is the potential energy of the interaction forces between particles i and j and is experienced by both particle i and particle j. The work-energy theorem also holds for many particles, that is, the total kinetic energy and the total potential energy is constant. The total kinetic energy of a system of particles is just the sum of the individual kinetic energies X 1 T = m~r˙ 2. 2 i i Putting these two results together, we have that the total energy of a system of n particles is n n n X 1 X X X E = m~r˙ 2 + U ext(~r ) + U (~r − ~r ), 2 i i i ij i j i i i j>i where the first term is the total kinetic energy, the second term is the total potential energy due to external forces and the last term is the total potential energy due to the internal forces. An important special case of this is that of a rigid body. For a rigid body, the internal forces (i.e. the interatomic forces) are central forces and so the potential energy does not depend on the direction of (~ri − ~rj). Furthermore, by definition of a rigid body, the quantity |~ri − ~rj| does not change for any given particles i and j. Therefore, for rigid bodies, the potential energy due to the internal forces is just some constant, which we can neglect. For a rigid body, we can also split the total kinetic energy into the translational 2 2 kinetic energy Mvcm/2 and the rotational kinetic energy Icmω /2. So the total energy of a rigid body can be expressed as
1 1 E = Mv2 + I ω2 + U ext + constant. 2 cm 2 cm 3.5. Summary: Energy 47
3.5 Summary: Energy
Skills to Master • Calculate the kinetic energy of a moving object • Identify whether a force is conservative and describe what this implies • Calculate the work done in moving an object against a force • Use conservation of mechanical energy to solve problems • Calculate the force given the potential energy and vice versa • Determine the qualitative behavior of a particle in 1D from a graph of its potential energy
Kinetic energy is Typically, we choose the reference point to be the origin ~r0 = (0, 0), and the final point to be ~r = (x, y), then we 1 1 T = mv2 = m~v · ~v. can break this line integral into a pair of straight line 2 2 integrals, the first from (0, 0) to (x, 0) and the second The work done in moving a particle from 1 to 2 is from (x, 0) to (x, y), then x y 2 0 0 U(~r) = U(x, y) = − Fx dx − Fy dy . W1→2 = F~· d~r. ˆ0 ˆ0 ˆ1 If F~ is in 3-dimensions, we have to add a third in- The work-kinetic energy theorem states that tegral for the path in the z-direction. A force is a spherically-symmetric central force if U = U(r), that ∆T = T2 − T1 = W1→2. is, if the potential energy depends only on r. Then us- ing the gradient operator for spherical coordinates, we To calculate the work done, if given the initial and final get F~ = − ∂U rˆ. Given F~, we can calculate U(r) using speeds, we can just use the work energy theorem. ∂r U(r) = − r f(r0) dr0. A force F~ is conservative if and only if any of the r0 The force´ vector function can be calculated from following are true the scalar potential function by taking the gradient • It depends only on position F~ = −∇~ U. • The work done is independent of the path taken • It’s a one dimensional system, that is, the posi- Given a conservative force field F~, we can learn a tion is described by a single parameter. lot from it: • It is a spherically-symmetric central force 1. Verify that it is conservative. • ∇~ × F~ = ~0 2. Calculate the potential function using U(~r) = − ~r F~· d~r. ~r0 A conservative force implies that 3. For´ 1-dimensional systems, we can go further: • A potential energy function can be defined a) Plot U(x). • Mechanical energy is conserved, that is, E = b) Characterize the equilibrium points of U(x) T + U = constant. by differentiating once to find them and a second time to determine whether they’re Given a conservative force F~, we can define a po- stable or unstable. If d2U/dx2 < 0 it is un- tential energy function stable. If d2U/dx2 > 0 it is stable. If the second derivative test is inconclusive, use ~r Taylor series to expand U(x) and then take U(~r) = − F~· d~r. the second derivative of that. ˆ~r0 48 Energy
c) For small oscillations about a stable equilib- forces and the last term is the total potential energy rium point, d2U/dx2 evaluated at the equi- due to the internal forces. librium point gives the effective spring con- The total energy of a rigid body is stant k. Then the angular frequency of the p oscillations is ω = k/m and the period of 1 2 1 2 ext E = Mv + Icmω + U + constant. the oscillations is τ = 2π/ω. 2 cm 2 d) Write t(x) as an integral by solving E = mx˙ 2/2 + U(x) forx ˙ and integrating. If A is For a one-dimensional system, find the position given, you can calculate E since at x = A, function x(t) by differentiating the total energy and ˙ ˙ v = 0, so E = U(A). Then solve the inte- solving the resulting differential equation T + U = 0. gral for t(x). If given the amplitude A, find When differentiating, remember to differentiate with the period τ by integratingx ˙ from 0 to the respect to t rather than x, using the product rule where amplitude A and multiplying by 4. appropriate. e) Invert t(x) to get x(t). The total energy of an ideal Atwood machine with f) Calculate the velocity functionx ˙(t) by dif- masses M and m and string length L is just the sum ferentiating the position function x(t). of the kinetic and potential energies The total energy of a system of n particles is 1 2 1 2 n n n E = Mx˙ + mx˙ − Mgx − mg(L − x). X 1 X X X 2 2 E = m~r˙ 2 + U ext(~r ) + U (~r − ~r ), 2 i i i ij i j i i i j>i If the pulley is massive, the only difference is that where the first term is the total kinetic energy, the sec- you must include the pulley’s rotational kinetic energy ond term is the total potential energy due to external Iω2/2 in the total energy. Chapter 4
Oscillations
4.1 Simple Harmonic Oscillators
Consider a generic 1-dimensional potential energy due to a conservative force. We are typically interested in the behavior of the system in the small oscillations near a stable equilibrium. If the generic potential energy function U(x) has a stable equilibrium at Tip some point, then we redefine the coordinate system so that this equilibrium point is at x = 0. We can Taylor expand U about the equilibrium point x = 0 as Any system near a stable equilibrium point can be dU 1 d2U approximated as a simple U(x) = U(0) + x + x2 + ··· dx 0 2 dx2 0 harmonic oscillator. We know that dU/dx evaluated at x = 0 is 0 because of the equilibrium point there. The first term, U(0), is just some irrelevant constant that we can neglect by redefining the reference point of the system. If we let
d2U k = , dx2 0 (an effective spring constant), then
1 U(x) = kx2, 2 for small oscillations about the equilibrium point. This is the potential energy of a simple harmonic oscillator. The key point is that we can approximate any system near a stable equilibrium point as a simple harmonic oscillator. Since force is related to potential energy via F (x) = −dU/dx, we get a linear restoring force F (x) = −kx, called Hooke’s law. Applying Newton’s second law gives us the equation of motion for the oscillating object r k x¨ = −ω2x, ω = . 0 0 m
The quantity ω0 is the angular frequency of the oscillation. The period is τ0 = 2π/ω0. The solution of this second order differential equation is
0 0 x(t) = A cos(ω0t) + B sin(ω0t).
Taking the second derivative of x(t) and plugging it intox ¨ given above, verifies the solution. Since this form of the solution contains both a cosine and a sine term, the 50 Oscillations
overall motion of the oscillator isn’t easy to visualize. However, two particular cases are U(x) easy to visualize. If the oscillator is pulled to x = x0 and released at t = 0, then the sine term vanishes, and the motion is completely described by the cosine term. If the oscillator is at x = 0 and is given a sharp impulse at t = 0, then the cosine term vanishes, a x and the motion is completely described by the sine term. Using U0 1 1 cos x = (eix + e−ix), sin x = (eix + e−ix), 2 2i we can also write the solution in the form Figure 4.1: The Morse poten- iω0t −iω0t tial. x(t) = C1e + C2e ,
0 0 0 0 where C1 = A /2 + B /(2i) and C2 = A /2 − B /(2i). The benefit of this form of the solution is that the exponentials are easy to work with—easy to integrate and differentiate. However, it isn’t easy to visualize the motion of the oscillator in this form. We can write x(t) a third way if we define an angle δ. Consider the right triangle 0 0 with angle δ, hypotenuse A, adjacent√ side A = A cos δ, and opposite side B = A sin δ. Then δ = tan−1(B0/A0) and A = A02 + B02. Then using the trigonometry identity cos(α + β) = cos α cos β − sin α sin β, we can write
0 0 x(t) = A cos(ω0t) + B sin(ω0t) = A cos δ cos(ω0t) + A sin δ sin(ω0t)
= A (cos(ω0t) cos δ + sin(ω0t) sin δ) = A cos(ω0t − δ).
The nice thing about this form of the solution is that A is the real amplitude of the motion and the phase angle δ is the real phase shift of the motion. With this form of the solution, the motion of the oscillator is very easy to visualize—it is simply a cosine curve with amplitude A and shifted to the right by δ. From here, we can go to a third form of the solution by noting that
i(ω0t−δ) Ae = A cos(ω0t − δ) + Ai sin(ω0t − δ).
On the left side, we have Aei(ω0t−δ) = Aeiω0te−iδ. Taking the real part of both sides gives us iω0t Re Ce = A cos(ω0t − δ), where C = Ae−iδ is a complex constant. The significance of this solution is that the particle is undergoing uniform circular motion in the complex plane, but its projection on the real axis is undergoing simple harmonic motion. From here, we can go back to any of the other solutions again. For example, using Euler’s formula again and the exponential form of the cosine function, we get
A Re Ceiω0t = A cos(ω t − δ) = ei(ω0t−δ) + e−i(ω0t−δ) 0 2 A A = e−iδeiω0t + eiδe−iω0t = C eiω0t + C e−iω0t. 2 2 1 2
−iδ iδ where C1 = Ae /2 and C2 = Ae /2. So for a simple (i.e. undamped and undriven) harmonic oscillator with angular oscillation frequency ω, we can write the displacement function x(t) in four different ways:
x(t) = A0 cos(ωt) + B0 sin(ωt) x(t) = A cos(ωt − δ) iωt −iωt x(t) = C1e + C2e x(t) = Re Ceiωt 4.1. Simple Harmonic Oscillators 51
Example 4.1.1
The Morse potential
2 − x−a U(x) = U0 1 − e δ − U0,
as graphed in Fig. (4.1), is a simple model for 1D motion of a diatomic molecule. The quantity x is the distance between the atoms. The quantity U0 is a convenient constant, U(x = a) = −U0. To find the equilibrium point(s), we set the derivative of U equal to zero. Differentiating gives us
dU 2U0 − x−a − x−a = 1 − e δ e δ . dx δ Setting it equal to zero, we find that an equilibrium point occurs at x = a. This tells us that a is the average bond length. That is, if the diatomic molecule is disturbed, it will oscillate and eventually settle to a distance of x = a between the atoms. Recall that F (x) = −dU/dx. By examining, the equation for dU/dx, we can see that for x < a =⇒ x − a < 0 =⇒ −(x − a)/δ > 0 =⇒ a strong repulsive force between the atoms. On the other hand, x > a =⇒ x − a > 0 =⇒ −(x − a)/δ < 0 =⇒ an attractive force between the atoms. Near x = a, the diatomic molecule acts like a harmonic oscillator, and to find the proportionality constant k (from Hooke’s law), we start by taking the second derivative of U(x) since we know that
dU d2U F = − = −kx, implies = k. dx dx2 x=0 We can also think of it in terms of a Taylor expansion about x = a. Taylor expanding, we get that 2 − x−a x − a 1 x − a e δ = 1 − + + ··· . δ 2 δ So for |(x − a)/δ| << 1, we have that
x − a2 x − a2 U U(x) ≈ U 1 − 1 − − U ≈ U − U ≈ 0 (x − a)2 − U . 0 δ 0 0 δ 0 δ2 0
2 2 2 Notice that the term U0(x − a) /δ is analogous to kx /2, so we have that 2U k = 0 . δa Differentiating U(x) gives us the force dU 2U F (x) = − = − 0 (x − a), dx δ2 which is just a translated version of Hooke’s law. If we pretend that we’re dealing with a single mass now (e.g. one atom might be much larger than the other), then we can calculate an equation of motion by applying Newton’s second law. We have that 2U mx¨ = − 0 (x − a). δ2 52 Oscillations
If we make the substitution z = x − a, then we get the straightforward differential equation 2 z¨ = −ω0z, p where ω0 = k/m is the angular frequency of the oscillation—something which can be measured with a spectrometer, and k is the value given above. The solution to this ODE is z(t) = A cos(ω0t − δ), and when we back-substitute, we get
x(t) = A cos(ω0t − δ) + a.
Differentiating x(t) = A cos(ω0t − δ), gives us the velocity of the oscillator as a function of time v =x ˙ = −ω0A sin(ω0t − δ). Then its kinetic energy as a function of time is 1 1 T = mv2 = mω2A2 sin2(ω t − δ). 2 2 0 0 Its potential energy is 1 1 U = kx2 = kA2 cos2(ω t − δ). 2 2 0 Its total mechanical energy is 1 1 T + U = mω2A2 sin2(ω t − δ) + kA2 cos2(ω t − δ) 2 0 0 2 0 1 1 = kA2 sin2(ω t − δ) + kA2 cos2(ω t − δ). 2 0 2 0 This simplifies to 1 T + U = kA2. 2 So the total energy of a harmonic oscillator depends only on its amplitude. An interesting and perhaps counter-intuitive note on oscillators is that an impulse applied to an oscillator does not change the oscillator’s frequency. It only changes its amplitude. For example, if you kick an oscillating mass attached to a spring, the mass won’t oscillate faster or slower, it will just oscillate with a larger amplitude. This is because an impulse changes the kinetic energy of the oscillator, therefore, it changes the amplitude, but not the frequency. The frequency depends only on the mass and the spring constant, and these are system characteristics rather than initial conditions.
4.2 Damped Harmonic Oscillators
All simple harmonic oscillators (i.e. undamped) have a natural frequency ω0. No matter how far such an oscillator is displaced from the equilibrium position, it will oscillate p at its natural frequency. This is because ω0 = k/m depends only the physical properties of the oscillator, k and m, and not on the amplitude A or the initial position or velocity of the oscillator. We will now deal with damped harmonic oscillators, which will no longer be oscillating at their natural frequencies. In the real world, oscillators are subject to friction forces, which damp the oscillations. The oscillations of a block on a spring and a pendulum will both die out in the long-term due to friction forces. The sound of a ringing bell (produced by the oscillations of the metal) will also die out after a time. Damping removes energy from the system. 4.2. Damped Harmonic Oscillators 53
A harmonic oscillator subject to a velocity dependent friction force is a damped oscillator. The equation of motion becomes
mx¨ = −kx − bx.˙ This is usually written in the form
2 x¨ + 2βx˙ + ω0x = 0, (4.1) p where ω0 = k/m is the natural frequency or the regular “undamped” angular fre- quency, and β = b/2m is the damping parameter. We can also write this more compactly in operator notation as Dx = 0, where D is the linear differential operator d2 d + 2β + ω2. dt2 dt 0 The solution to Eq. (4.1) is rich with multiple interesting regimes. Since it is a second order differential equation, there will be two independent solutions. If we have two solutions, then by the principle of superposition, the linear combination of the two solutions will also be a solution. Since it is a second order ODE, we need two linearly independent solutions, which we combine to get the general solution. The following method will work for any second order linear homogeneous differential equation. We start by letting x(t) = ert, thenx ˙ = rert andx ¨ = r2ert. Substituting these into the equation of motion gives us
2 x¨ + 2βx˙ + ω0x = 0 2 rt rt 2 rt r e + 2βre + ω0e = 0 rt 2 2 e r + 2βr + ω0 = 0 2 2 r + 2βr + ω0 = 0. This last equation is called the auxiliary equation. In words, x(t) = ert is a solution if 2 rt rt 2 rt 2 2 r e + 2βre + ω0e = 0, and this is true if the quadratic equation r + 2βr + ω0 = 0. Using the quadratic formula, we see that this equation is zero in two cases: q 2 2 r1 = −β + β − ω0 q 2 2 r2 = −β − β − ω0. The general solution has the form
r1t r2t x(t) = A1e + A2e . In our case, the general solution is √ √ −βt β2−ω2t − β2−ω2t x(t) = e A1e 0 + A2e 0 .
There are three distinct cases: 2 2 1. β − ω < 0, that is, β < ω0. In this case, the exponents on e are imaginary. This is the regime of underdamping. 2 2 2. β − ω = 0, that is, β = ω0. In this case, the exponents on e are zero. This is the regime of critical damping. 2 2 3. β − ω > 0, that is, β > ω0. In this case, the exponents on e are real. This is the regime of overdamping. Notice that if β = 0, we recover the undamped solution, as expected. 54 Oscillations
Underdamped Oscillator x(t) A The general solution for damped oscillations is √ √ −βt β2−ω2t − β2−ω2t x(t) = e A1e 0 + A2e 0 . t 2 2 In the underdamped case, β − ω0 < 0, which means the exponents in the solution are p 2 2 p 2 2 imaginary. By factoring a negative out of the radical, we get β − ω0 = i ω0 − β , −A and now the radical is a real quantity. We define this quantity, which is a frequency lower p 2 2 than the natural frequency ω0, as ω1 = ω0 − β . Now our solution can be written as x(t) = e−βt A eiω1t + A e−iω1t . Figure 4.2: Position function of 1 2 an underdamped oscillator. This is a complex quantity, but in the end, this doesn’t matter because we can always just take the real part, which is also a solution. A complex function of time is real for all time if it equals its complex conjugate. The complex conjugate of x(t) is
∗ −βt ∗ −iω1t ∗ iω1t x(t) = e A1e + A2e . ∗ ∗ In order for x(t) = x(t) , it must be that A1 = A2. That is, one constant must be the ∗ ∗ complex conjugate of the other. So we can define the new constants C = A1 = A2 and ∗ C = A2 = A1, so that x(t) = e−βt C∗eiω1t + Ce−iω1t . The pair of terms in parentheses are now complex conjugates, and we know that when complex conjugates are added, the imaginary parts disappear. However, we can write this in a more convenient form. Any complex number can be written in exponential form containing a real amplitude and a real angle. If we write C = Aeiδ/2, then C∗ = Ae−iδ/2. Substituting these values gives us A x(t) = e−βt ei(ω1t−δ) + e−i(ω1t−δ) , 2 and so the general solution for underdamped oscillations can be written as q −βt 2 2 x(t) = Ae cos(ω1t − δ), ω1 = ω0 − β (underdamped). (4.2)
This tells us that in the underdamped case, the oscillator is oscillating with a magnitude −βt Ae and a frequency ω1, which is less than the natural frequency ω0. Notice that the amplitude Ae−βt is no longer constant in time. In fact, it is decaying exponentially, and we call Ae−βt the envelope function since the oscillations stay between Ae−βt and −Ae−βt. This can be seen in Fig. (4.2). Strictly speaking, this motion is non-periodic, but we can define the period as the time to go from peak to peak. This period is 2π τ1 = , ω1 and can be found either by differentiating Eq. (4.2) and looking at twice the time intervals between extrema or by looking at twice the time period between consecutive zeros of Eq. (4.2). The decrement of motion refers to the rate at which the amplitudes decrease and can be found as the ratio of the amplitude of a peak to the amplitude of the next peak Ae−βt . Ae−β(t+τ1)
Notice that the envelope function Ae−βt drops to 1/e of its initial value in time 1/β. This time is called the decay time for the oscillator. 1 decay time = . β 4.2. Damped Harmonic Oscillators 55
Critically Damped Oscillator
2 2 x(t) When β − ω0 = 0, that is, β = ω0, the system is said to be critically damped. In this A case, the general solution for damped oscillations √ √ −βt β2−ω2t − β2−ω2t t x(t) = e A1e 0 + A2e 0 , reduces to the single solution x(t) = e−βt. −A However, a general solution to a second order differential equation is the sum of two solutions, so we have to find a second solution. It turns out that Figure 4.3: Position function of critically damped oscillator. x(t) = te−βt, is also a solution in the special case that β = ω0. This can be verified by taking the first and second derivatives and plugging them intox ¨ + 2βx˙ + β2x = 0. Adding the two solutions together, we get the general solution for the critically damped case
x(t) = (A + Bt)e−βt (critically damped).
This is shown in Fig. (4.3). Critical damping brings the object to the equilibrium position as fast as possible without going past the equilibrium position. Critical damping is usually optimal if you’re designing something in which you want the oscillations damped as quickly as possible.
Overdamped Oscillator Perhaps unexpectedly, the oscillations of a system that is more than critically damped will decay more slowly than if it was critically damped. In a way, overdamping gets in the way of itself. 2 2 When β − ω0 > 0, that is, β > ω0, the system is said to be overdamped. In this case, the general solution for damped oscillations √ √ −βt β2−ω2t − β2−ω2t x(t) = e A1e 0 + A2e 0 , becomes q −βt ω2t −ω2t 2 2 x(t) = e A1e + A2e , ω2 = β − ω0 (overdamped).
Note that ω2 is a new frequency. After a long time,
−(β−ω2)t x(t) ≈ A1e .
An overdamped system returns to equilibrium more slowly than an underdamped or a critically damped system. Unlike an underdamped system, there are no oscillations in an overdamped system—the system is so strongly damped that if displaced from equilibrium it returns to equilibrium only after infinite time.
Energy Dissipation Recall that damping is removing energy from the oscillating system. Since damping isn’t conservative, the mechanical energy is not conserved, so dE/dt is no longer zero. With damping, the equation of motion for a harmonic oscillator is
mx¨ = −kx − bx,˙ 56 Oscillations
x(t) the energy is 1 1 A E = mx˙ 2 + kx2, 2 2 and the time derivative is t dE 1 2 3 =x ˙(mx¨ + kx). dt We applied the chain rule here. Plugging in mx¨ from the equation of motion gives us −A dE = −bx˙ 2, dt Figure 4.4: Undamped oscilla- tor from the example. which is the rate of energy dissipation. The negative sign indicates that damping is always taking energy out of the system. We can also look at the rate of work done by the resistive force F = −bx˙ x(t) A dW = −bx˙ dx. Dividing by dt gives us t dW dx = −bx˙ = −bx˙ 2. 1 2 3 dt dt
−A Example 4.2.1 Consider a typical harmonic oscillator in the form of a block connected to a Figure 4.5: Damped oscillator spring all resting on a frictionless surface. The other end of the spring is connected from the example. to the wall and the whole experiment is being done in a vacuum. The system is motionless when you pull the block (at the end of the spring) 1 meter past its equilibrium point (motionless resting position) and release it. It then begins −1 oscillating with an angular frequency ω0 = 2π s with a period of τ0 = 2π/ω0 = 1 second. That is, the block is going from equilibrium to 1 meter past equilibrium, back through equilibrium to 1 meter past in the other direction, and returning to equilibrium once every second. The position equation for an undamped harmonic oscillator is
x(t) = A cos(ω0t − δ).
−1 In our case, the amplitude is A = 1m, the natural frequency is ω0 = 2π s , and the phase shift is δ = 0. The position equation is simply
x(t) = 1 m cos(2π s−1 t).
A plot of the displacement is shown in Fig. (4.4). Let’s consider now a more realistic version of the experiment. We let air back into our experiment and account for friction between the tabletop and the oscillating block. The general solution for the underdamped harmonic oscillator is q −βt 2 2 x(t) = Ae cos(ω1t − δ), ω1 = ω0 − β .
−1 p 2 2 −1 If we have β = 0.5 s , then ω1 = 4π − (0.5) ≈ 6.26 s . Then our position function is −1 x(t) = (1 m) e−0.5 s t cos(6.26 s−1t). As we can see from the graph in Fig. (4.5), after only 3 seconds, the block is reaching only about a fourth as far from the equilibrium point of x = 0 as it did when it started. −1 What happens if the friction is adjusted perfectly such that β = ω0 = 2π s ? That is, when the system becomes critically damped? For a critically damped 4.3. Driven Oscillators 57
system, the position equation is x(t) A x(t) = (A + Bt)e−βt.
We know that at t = 0, x = 1m, andx ˙ = 0 m/s. This allows us to compute t A = 1m and B = 2π m/s. Our position function is 1 2 3
−1 x(t) = (1 m + 2π m/s t)e−2π s t. −A
Notice from the graph in Fig. (4.6) that the block is practically at rest at the Figure 4.6: Critically damped equilibrium position in a very short time—just a little over 1 second. At t = τ0 = 1 s, the block’s position is oscillator from the example.
−1 x(t = 1 s) = (1 m + 2π m/s 1 s)e−2π s 1 s = (1 + 2π)e−2πm ≈ 0.014 m. x(t) So in the time it takes the undamped oscillator to complete one period, the crit- A ically damped oscillator has decayed to 0.014 m. This is true in general. If you wait the period of the undamped oscillator, the critically damped one will have t decayed to 0.014x0, and you can’t get better performance than this. 1 2 3 Finally, what happens if we increase the friction beyond the point of critical damping? For overdamping, the position function is −A
−βt ω2t −ω2t x(t) = e A1e + A2e , Figure 4.7: Overdamped oscil- p 2 2 lator from the example. where ω2 = β − ω0. Setting x(t = 0) andx ˙(t = 0) to zero, we can get a system of equations in A1 and A2. Using the fact that x(t = 0) = 1m and x˙(t = 0) = 1 m/s, we find that A1 = 1/2 + β/(2ω2) and A2 = 1/2 − β/(2ω2). If −1 β = 20, then ω2 = 19.0 s , A1 = 1.026 and A2 = 0.0263, so our position function becomes x(t) = e−20t 1.026e19t + 0.0263e−19t . In the graph in Fig. (4.7), we can see that the damping (i.e. friction) is so high that it takes the block a lot longer to move back to the equilibrium position than when it was critically damped.
4.3 Driven Oscillators
The equation of motion for a damped harmonic oscillator with a sinusoidal driving force is
mx¨ = −kx − bx˙ + F0 cos ωt, where −bx˙ is the damping force and F0 cos ωt is the driving force. Frequently, we write the equation of motion in the form
2 x¨ + 2βx˙ + ω0x = f0 cos ωt,
p where β = b/2m is the damping parameter, ω0 = k/m is the natural frequency (i.e. the frequency in the absence of damping), ω is the frequency of the driving force, and f0 = F0/m. This is now an inhomogeneous second order differential equation. Inhomogeneous since it has a term f0 cos ωt that does not depend on x or any of its derivatives. In operator notation, we can write the equation of motion as
Dx = f0 cos ωt, 58 Oscillations
d2 d 2 where D = dt2 + 2β dt + ω0 as before. After some time, the motion associated with the natural frequency ω0 will die down x(t) because of the damping force, and all that will be left will be motion due to the driving frequency ω.
A t Driven Undamped Oscillators −0.5 0.5 1 1.5 2 −A In the undamped case, β = 0, and our equation becomes 2 x¨ + ω0x = f0 cos ωt. 2 In this case, if ω << ω0, then the ω0x term dominates. If ω >> ω0, then the driving Figure 4.8: Driven undamped term dominates the motion. oscillator. For the particular, or steady-state solution, we guess xs(t) = A1 cos ωt since both 2 x andx ¨ will contain cosine. Taking the second derivative, we getx ¨ = −A1ω cos ωt. Plugging this in gives us 2 x¨ + ω0x = f0 cos ωt 2 2 −A1ω cos ωt + ω0A1 cos ωt = f0 cos ωt f0 A1 = 2 2 . ω0 − ω 2 2 So xs(t) = A1 cos ωt is the unique steady-state solution and A1 = f0/(ω0 − ω ). Notice 2 2 that the sign of A depends on which of ω0 or ω is larger. We can write the steady-state solution as xs(t) = A1 cos(ωt − δ1),
if we redefine A1 as ( f0 0 if ω < ω0 A1 = 2 2 , δ1 = |ω0 − ω | π if ω > ω0.
The δ1 serves to change the sign of the function. The key thing to remember is that in the driven undamped case, the amplitude is 1 A ∝ 2 2 . |ω0 − ω |
We know that the homogeneous solution is xh(t) = A2 cos(ω0t − δ2) where A2 is the p amplitude in the undriven and undamped case, ω0 = k/m is the natural frequency of the oscillator, and δ2 is the phase shift in the undriven and undamped case. The general solution is the sum of the homogeneous and steady-state solutions
x(t) = xs(t) + xh(t) = A1 cos(ωt − δ1) + A2 cos(ω0t − δ2). Consider a simple harmonic oscillator, say a block at the end of a spring on a fric- tionless table. The oscillator is released/driven from rest, has a natural frequency of ω0 = 10π, and the driving force has frequency ω = 2π and amplitude f0 = 1000. Since ω < ω0, we have that δ1 = 0. From the amplitude of the driving force and the two frequencies, we calculate A1 = 1.055. We chose such a large value for f0 in order to get a convenient unit for A1. Plugging in what we have so far, gives us
x(t) = 1.055 cos(2πt) + A2 cos(10πt − δ2).
By solving x(t) andx ˙(t) for the initial conditions, we get δ2 = 0, and A2 = −1.055. So the position function is x(t) = 1.055 cos(2πt) − 1.055 cos(10πt). The dashed line in the graph shown in Fig. (4.8) shows the driving force cos(2πt), and the solid line shows the resulting motion x(t). 4.3. Driven Oscillators 59
Driven Damped Oscillators We now consider the damped case of the driven oscillator. The equation of motion is
2 x¨ + 2βx˙ + ω0x = f0 cos ωt, where the first term is the acceleration of the oscillator, the second term is the damping force, the third term is the restoring force, and the right side is the driving force. Using the operator notation defined earlier, we can also write this as
Dx = f0 cos ωt. This is a second order inhomogeneous linear differential equation, so we expect two constants. We only have to find one particular solution of Dx = f0 cos ωt, then we can add any solution of the homogeneous equation Dx = 0. We know that the solution for r1t r2t the homogeneous case is xh(t) = A1e + A2e because we solved it when studying the undriven damped oscillator. The general solution of Dx = f0 cos ωt is just the sum of the particular solution and the homogeneous solution
x(t) = xp(t) + xh(t).
r1t r2t Remember that in the long-term, the homogeneous solution xh(t) = A1e + A2e goes to zero because of the damping. So in the long-term the particular solution xp(t) is the only surviving part of x(t), so we call it the steady state solution. The steady state part of the solution survives in the long-term because of the driving force. The terms from the homogeneous part of x(t) go down to zero in the long-term, so we call them the transients. 2 To find the particular solution forx ¨ + 2βx˙ + ω0x = f0 cos ωt, it is convenient to consider 2 y¨ + 2βy˙ + ω0y = f0 sin ωt. We do this because we know that for any solution x(t), there is a corresponding solution y(t) with the cosine replaced by the sine term because they are only shifted in time. If we multiply the above equation by i and add it to the one containing x, we get
2 iωt x¨ + iy¨ + 2β (x ˙ + iy˙) + ω0 (x + iy) = f0 cos ωt + if0 sin ωt = f0e . If we let z = x + iy, we can write this as
2 iωt z¨ + 2βz˙ + ω0z = f0e . Now we solve for z(t) and then x(t) = Re (z(t)). If we try z(t) = Ceiωt where C is a complex constant, thenz ˙ = iωCeiωt andz ¨ = −ω2Ceiωt. Plugging them in gives us
2 iωt iωt 2 iωt iωt (−ω Ce ) + 2β(iωCe ) + ω0(Ce ) = f0e f0 C = 2 2 . ω0 − ω + i2ωβ Any complex number can be written in exponential form with a real-valued amplitude A and a real-valued phase angle δ. We write C in exponential form as
f0 −iδ C = 2 2 = Ae . ω0 − ω + i2ωβ Recall that multiplying a complex number by its complex conjugate gives the magnitude squared of the complex number, so 2 ∗ −iδ iδ f0 f0 A = CC = (Ae )(Ae ) = 2 2 2 2 ω0 − ω + i2ωβ ω0 − ω − i2ωβ 2 f0 = 2 2 2 2 2 . (ω0 − ω ) + 4ω β 60 Oscillations
Remember that taking the complex conjugate of even a complicated expression is very easy—just change the sign of every i that appears in the expression. From C = Ae−iδ, x(t) we see that A/C = eiδ, and we can rewrite eiδ as cos δ + i sin δ, so plugging C back in gives us A A 2 2 A cos δ + i sin δ = = (ω0 − ω + i2ωβ). t C f0 −0.5 0.5 1 1.5 2 −A Separating the right side into real and imaginary parts and equating the real and imagi- nary parts implies that A sin δ = 2ωβ f0 Figure 4.9: Driven damped har- A 2 2 monic oscillator. cos δ = (ω0 − ω ). f0 Dividing the first equation by the second gives us 2ωβ tan δ = 2 2 ω0 − ω −1 2ωβ δ = tan 2 2 . ω0 − ω Recall that our goal was to find the solution z(t) and then take the real part of it. Our solution is z(t) = Ceiωt. Taking the real part gives us our particular solution