- •Preface
- •Contents
- •1 Introduction
- •1.1 Physics
- •1.2 Mechanics
- •1.3 Integrating Numerical Methods
- •1.4 Problems and Exercises
- •1.5 How to Learn Physics
- •1.5.1 Advice for How to Succeed
- •1.6 How to Use This Book
- •2 Getting Started with Programming
- •2.1 A Python Calculator
- •2.2 Scripts and Functions
- •2.3 Plotting Data-Sets
- •2.4 Plotting a Function
- •2.5 Random Numbers
- •2.6 Conditions
- •2.7 Reading Real Data
- •2.7.1 Example: Plot of Function and Derivative
- •3 Units and Measurement
- •3.1 Standardized Units
- •3.2 Changing Units
- •3.4 Numerical Representation
- •4 Motion in One Dimension
- •4.1 Description of Motion
- •4.1.1 Example: Motion of a Falling Tennis Ball
- •4.2 Calculation of Motion
- •4.2.1 Example: Modeling the Motion of a Falling Tennis Ball
- •5 Forces in One Dimension
- •5.1 What Is a Force?
- •5.2 Identifying Forces
- •5.3.1 Example: Acceleration and Forces on a Lunar Lander
- •5.4 Force Models
- •5.5 Force Model: Gravitational Force
- •5.6 Force Model: Viscous Force
- •5.6.1 Example: Falling Raindrops
- •5.7 Force Model: Spring Force
- •5.7.1 Example: Motion of a Hanging Block
- •5.9.1 Example: Weight in an Elevator
- •6 Motion in Two and Three Dimensions
- •6.1 Vectors
- •6.2 Description of Motion
- •6.2.1 Example: Mars Express
- •6.3 Calculation of Motion
- •6.3.1 Example: Feather in the Wind
- •6.4 Frames of Reference
- •6.4.1 Example: Motion of a Boat on a Flowing River
- •7 Forces in Two and Three Dimensions
- •7.1 Identifying Forces
- •7.3.1 Example: Motion of a Ball with Gravity
- •7.4.1 Example: Path Through a Tornado
- •7.5.1 Example: Motion of a Bouncing Ball with Air Resistance
- •7.6.1 Example: Comet Trajectory
- •8 Constrained Motion
- •8.1 Linear Motion
- •8.2 Curved Motion
- •8.2.1 Example: Acceleration of a Matchbox Car
- •8.2.2 Example: Acceleration of a Rotating Rod
- •8.2.3 Example: Normal Acceleration in Circular Motion
- •9 Forces and Constrained Motion
- •9.1 Linear Constraints
- •9.1.1 Example: A Bead in the Wind
- •9.2.1 Example: Static Friction Forces
- •9.2.2 Example: Dynamic Friction of a Block Sliding up a Hill
- •9.2.3 Example: Oscillations During an Earthquake
- •9.3 Circular Motion
- •9.3.1 Example: A Car Driving Through a Curve
- •9.3.2 Example: Pendulum with Air Resistance
- •10 Work
- •10.1 Integration Methods
- •10.2 Work-Energy Theorem
- •10.3 Work Done by One-Dimensional Force Models
- •10.3.1 Example: Jumping from the Roof
- •10.3.2 Example: Stopping in a Cushion
- •10.4.1 Example: Work of Gravity
- •10.4.2 Example: Roller-Coaster Motion
- •10.4.3 Example: Work on a Block Sliding Down a Plane
- •10.5 Power
- •10.5.1 Example: Power Exerted When Climbing the Stairs
- •10.5.2 Example: Power of Small Bacterium
- •11 Energy
- •11.1 Motivating Examples
- •11.2 Potential Energy in One Dimension
- •11.2.1 Example: Falling Faster
- •11.2.2 Example: Roller-Coaster Motion
- •11.2.3 Example: Pendulum
- •11.2.4 Example: Spring Cannon
- •11.3 Energy Diagrams
- •11.3.1 Example: Energy Diagram for the Vertical Bow-Shot
- •11.3.2 Example: Atomic Motion Along a Surface
- •11.4 The Energy Principle
- •11.4.1 Example: Lift and Release
- •11.4.2 Example: Sliding Block
- •11.5 Potential Energy in Three Dimensions
- •11.5.1 Example: Constant Gravity in Three Dimensions
- •11.5.2 Example: Gravity in Three Dimensions
- •11.5.3 Example: Non-conservative Force Field
- •11.6 Energy Conservation as a Test of Numerical Solutions
- •12 Momentum, Impulse, and Collisions
- •12.2 Translational Momentum
- •12.3 Impulse and Change in Momentum
- •12.3.1 Example: Ball Colliding with Wall
- •12.3.2 Example: Hitting a Tennis Ball
- •12.4 Isolated Systems and Conservation of Momentum
- •12.5 Collisions
- •12.5.1 Example: Ballistic Pendulum
- •12.5.2 Example: Super-Ball
- •12.6 Modeling and Visualization of Collisions
- •12.7 Rocket Equation
- •12.7.1 Example: Adding Mass to a Railway Car
- •12.7.2 Example: Rocket with Diminishing Mass
- •13 Multiparticle Systems
- •13.1 Motion of a Multiparticle System
- •13.2 The Center of Mass
- •13.2.1 Example: Points on a Line
- •13.2.2 Example: Center of Mass of Object with Hole
- •13.2.3 Example: Center of Mass by Integration
- •13.2.4 Example: Center of Mass from Image Analysis
- •13.3.1 Example: Ballistic Motion with an Explosion
- •13.4 Motion in the Center of Mass System
- •13.5 Energy Partitioning
- •13.5.1 Example: Bouncing Dumbbell
- •13.6 Energy Principle for Multi-particle Systems
- •14 Rotational Motion
- •14.2 Angular Velocity
- •14.3 Angular Acceleration
- •14.3.1 Example: Oscillating Antenna
- •14.4 Comparing Linear and Rotational Motion
- •14.5 Solving for the Rotational Motion
- •14.5.1 Example: Revolutions of an Accelerating Disc
- •14.5.2 Example: Angular Velocities of Two Objects in Contact
- •14.6 Rotational Motion in Three Dimensions
- •14.6.1 Example: Velocity and Acceleration of a Conical Pendulum
- •15 Rotation of Rigid Bodies
- •15.1 Rigid Bodies
- •15.2 Kinetic Energy of a Rotating Rigid Body
- •15.3 Calculating the Moment of Inertia
- •15.3.1 Example: Moment of Inertia of Two-Particle System
- •15.3.2 Example: Moment of Inertia of a Plate
- •15.4 Conservation of Energy for Rigid Bodies
- •15.4.1 Example: Rotating Rod
- •15.5 Relating Rotational and Translational Motion
- •15.5.1 Example: Weight and Spinning Wheel
- •15.5.2 Example: Rolling Down a Hill
- •16 Dynamics of Rigid Bodies
- •16.2.1 Example: Torque and Vector Decomposition
- •16.2.2 Example: Pulling at a Wheel
- •16.2.3 Example: Blowing at a Pendulum
- •16.3 Rotational Motion Around a Moving Center of Mass
- •16.3.1 Example: Kicking a Ball
- •16.3.2 Example: Rolling down an Inclined Plane
- •16.3.3 Example: Bouncing Rod
- •16.4 Collisions and Conservation Laws
- •16.4.1 Example: Block on a Frictionless Table
- •16.4.2 Example: Changing Your Angular Velocity
- •16.4.3 Example: Conservation of Rotational Momentum
- •16.4.4 Example: Ballistic Pendulum
- •16.4.5 Example: Rotating Rod
- •16.5 General Rotational Motion
- •Index
10.3 Work Done by One-Dimensional Force Models |
285 |
|
|
|
|
|
|
|
Fig. 10.7 The measured force F (y) as a function of the position y of the top of the cushion. Notice that y is negative when the cushion is compressed
10.3.2 Example: Stopping in a Cushion
You have installed a new soft cushion from SoftCush technology to prevent falling damage. According to the producers, the cushion has been developed to produce a normal force that depends on the compression of the cushion only (and not on the speed of deformation). You do not know the functional form of the force response, F (y), of the cushion, but you have measured the response in a controlled experiment where you pushed the cushion down to a position y and measured the corresponding force F (y). The results are provided in the file cushionforce.d,3 and shown in Fig. 10.7.
Problem: If you jump from a height h and onto the cushion, how far will it compress before you stop? What is the maximum of the force from the cushion on you?
Approach: The cushion reaches its maximum compression when you stop, that is, when your kinetic energy is zero. We will first use the work-energy theorem to find the kinetic energy you have after falling a distance h, and then we will use the workenergy theorem during the contact with the cushion, from you touch the cushion and until you stop.
Sketch and Identify: Figure 10.8 shows an illustration of the forces acting on the person during the jump and a free-body diagram of the person.
Model: We divide the motion into two phases, phase 1 from x0 = h to x1 = 0, where the person is only affected by gravity, G = mg, and phase 2 from x1 = 0 to x2 = −d, where the person is affected by gravity, G, and the force from the cushion, F .
Phase 1: Free fall: In phase 1 the person starts from rest, v0 = 0 and K0 = 0, at a height x0 = h, and reaches the height, x1 = 0 with a velocity, v1, and a kinetic energy, K1.
3http://folk.uio.no/malthe/mechbook/cushionforce.d.
286 |
10 Work |
10000 |
|
|
|
|
|
|
5000 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
-5000 |
|
|
|
|
|
|
0 |
0.2 |
0.4 |
0.6 |
0.8 |
1 |
1.2 |
10 |
|
|
|
|
|
|
5 |
|
|
|
|
|
|
0 |
|
|
|
|
|
|
-5 |
|
|
|
|
|
|
0 |
0.2 |
0.4 |
0.6 |
0.8 |
1 |
1.2 |
Fig. 10.8 The measured force F (y) as a function of the position y of the top of the cushion. Notice that y is negative when the cushion is compressed
We find K1 by applying the work-energy theorem from 0 to 1:
W0,1 |
= |
x0 |
1 |
G d x = h |
0 |
−mg d x = mgh = K1 − K0 = K1 , |
(10.46) |
|
|
x |
|
|
|
|
|
where we have used that K0 = 0. We see that K1 = (1/2)mv12 = mgh, so that we can determine v1 if needed.
Phase 2: Contact with cushion: In phase 2 the person starts with velocity v1 and kinetic energy K1 at the height y1 = 0 and stops with velocity v2 = 0 and kinetic energy K2 = 0 at height y2 = −d. You are in contact with the SoftCusion, and the net force affecting you is therefore:
F net = F (x ) − mg , |
(10.47) |
We apply the work-energy theorem from 1 to 2 to determine d, the stopping distance:
W1,2 |
= |
x1 |
2 |
F net d x = 0 |
(F (x ) − mg) d x = K2 − K1 = −mgh , (10.48) |
|
|
x |
|
|
−d |
where we have used that K2 = 0 and K1 = mgh from (10.46).
The idea is that the value for d that satisfies (10.48) gives us the maximum compression of the cushion. Unfortunately, (10.48) is not an explicit equation in d that
10.3 Work Done by One-Dimensional Force Models |
287 |
we can solve since the unknown d is the upper limit in the integral. If we knew how to solve the integral analytically, we might have been able to solve the equation, but in this case we only know how to solve the integral numerically. How can we solve such an equation and find the d that satisfies (10.48)?
What we need to do, is to calculate the work:
x
W1,2 = F (x )d x − mgx , (10.49)
0
as a function of the position x of the jumper. Our plan is to vary x until we find a value for x that satisfies the equation, which corresponds to −d. We need to vary x systematically. We start from x = 0, where the jumper comes in contact with the cushion, and then gradually decrease x (remember that the jumper is moving down during the contact with the cushion) until the equation is satisfied, that is, until we find the value for x for which W1,2 is closest to −mgh. More precisely, we make
a sequence of x values starting from 0, decreasing in small steps, |
x : x = 0, |
|
|
|
0 |
x = x + x , and we calculate the work for each of these numbers until the work |
||
1 |
0 |
|
is equal to −mgh. |
|
|
|
Can we use any step size x in this sequence of x values? If we knew the force |
|
|
i |
|
F (x ) for any position x , we could calculate the work at any resolution |
x . But in this |
case we only know the forces F (xi ) for the values xi where it has been measured. Therefore, the best possible resolution we can get is to use the measured xi values as our sequence of xi -values. This means that we start by calculating the work for x1 = x1. Then we calculate the work for x2 = x2, then for x3 = y3 and so on, until, for some value i , the work is approximately equal to −mgh.
How do we calculate the work for x = xk , where k is a number in the sequence of xi values in the datafile? The integral in (10.49) can be calculated using the trapezoidal rule even if we do not know the underlying function—it works well also for a measured dataset:
0 |
xk |
k |
(xi +1 − xi ) |
2 |
F (xi ) + F (xi +1) |
. |
(10.50) |
F (x ) d x i =1 |
|||||||
|
|
|
|
1 |
|
|
|
Numerically, this is done using the function trapz. Here, we only want to do the sum over the first k values of xi , that is, we only want to do the sum from x1 to xk in the list of n such x -values. This corresponds to doing the sum over the first [0:k-1] elements in the position array x:
xval, F = loadtxt(’cushionforce.d’,usecols=[0,1],unpack=True) k = 10
I = trapz(F[0:k-1],x=xval[0:k-1])
We calculate the work for all possible values of k, and store the results in an array work, which gives us the work as a function of k:
from pylab import * m = 80.0 # kg
g = 9.8 # m/sˆ2
288 |
10 Work |
h = 6.0 |
# m |
x, F = loadtxt(’cushionforce.d’,usecols=[0,1],unpack=True) n = length(x)
work = zeros(n,float) for k in range(1,n-1):
I = trapz(F[0:k],x=-x[0:k]) work[k] = I - m*g*x[k]
plot(y,work,’-b’) xlabel(’x (m)’) ylabel(’W(x) (J)’)
Figure 10.9 shows the work W1,2 as a function of the position x calculated using this program. In the same plot we have also shown the value for −mgh. We see that the work starts at zero and decreases gradually as the deformation increasesa and x becomes a large negative number because the work done is negative: The net work acts to reduce the kinetic energy. When the work reaches the value −mgh the jumper stops: the kinetic energy is now zero. How can we find the x -value this occurs at from our numerical data? Based on the plot, we realize that the x value
when W1,2 |
= −mgh can be found as the first value in the sequence of xi ’s for |
which W1,2 |
< −mgh. This value of x will correspond to a work that is slightly |
smaller than −mgh, but this is a good first approximation. How can we find this value numerically? We use the function find:
i2 = min(find(work<mgh))
x2 = x(i2)
The function find(work<-mgh) returns a list of all the indexes of work that satisfies the condition work<-mgh. We need to find the first of the indexes in this sequence, therefore we take the minimum index in this list, and find the corresponding y-value. We have found the maximum compression of the cushion!
What is the cushion force at this value? This is the corresponding cushion force: F (y2) = 7050 N found by:
>>F2 = F[i2] >>print F2
F2 = |
7050 |
What would happen if the jumper started from h = 9 m instead? We plot the corresponding value for −mgh in Fig. 10.9. You can read the maximum compression directly from this graph, or determine it numerically as we did above.
Fig. 10.9 Plot of the work W1,2 as a function of the position x of the jumper
10.4 Work Done in Twoand Three-Dimensional Motions |
289 |
10.4 Work Done in Twoand Three-Dimensional Motions
The work-energy theorem was demonstrated for a three-dimensional motion. We have so far studied one-dimensional motions only. How does the application of the theorem change for three-dimensional motion?
Work of Tangential and Normal Forces
Figure 10.10 illustrates a general three-dimensional motion from 0 to 1. For this motion, the net work is
t1
W = Fnet · v dt = K1 − K0 . (10.51)
t0
Because v points in the tangential direction, we notice that it is only the tangential component of the force that does any work! We can demonstrate this by introducing a local coordinate system along the path of motion, with unit vectors uˆ T in the tangential direction and uˆ N in the normal direction. The velocity vector points along the tangential unit vector:
v = v uˆ T (t ) . |
(10.52) |
If we decompose a force, F j , in the tangential and normal directions along the path of motion:
F j = Fj,T (t ) uˆ T (t ) + Fj,N (t ) uˆ N (t ) ,
we see that the work done by the force F j is:
W j = |
t0 |
F j · dt dt |
|
|
t1 |
|
dr |
t1
= F
t0
Fig. 10.10 Illustration of the path followed by an object moving from position 0 to position 1
j,T (t )uˆ T (t ) + Fj,N (t )uˆ N (t ) · dr dt dt
F
FN
uT FT
uN
(10.53)
(10.54)
(10.55)
r(t1)
r(t0) r(t)
y
x
290 10 Work
= |
t0 |
Fj,T (t )uˆ T (t ) · dt |
dt + |
t0 |
Fj,N (t ) uˆ N (t ) · |
dt |
dt (10.56) |
||||||
|
t1 |
|
dr |
|
t1 |
|
|
|
|
dr |
|
||
|
t1 |
|
dr |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
=0 |
|
|
|
|||
= t0 |
Fj,T (t )uˆ T (t ) · |
|
dt . |
|
|
|
|
|
|
|
(10.57) |
||
dt |
|
|
|
|
|
|
|
where we have used that the normal unit vector, uˆ N , is normal to the velocity vector, v. Consequently, it is only the tangential component of the force that contributes to the work.
This is also intuitive, since we know that normal forces only contribute to a change in the direction of the velocity vector and not to a change in the speed—it is only the tangential force that can cause a tangential acceleration which causes a change in speed.
Work of a Constant Force in Two and Three Dimensions
The work done by a constant force has a particularly simple solution. Let us address the work done on an object by a constant force F as the object is moved from r(t0) to r(t1) as illustrated in Fig. 10.11. The work done on the object is
WF = t0t1 |
F · vdt = 0 |
1 |
F · dr , |
(10.58) |
Since the force is constant, we move it outside the integration:
WF = F · |
0 |
1 dr |
|
|
= F · (r(t1) − r(t0)) . |
(10.59) |
||||
|
(r |
|
|
|
|
|
|
) |
|
|
= |
1 |
0 |
) |
|
||||||
(t |
)−r(t |
|
|
We call s = r = r(t1) − r(t0) the displacement. The work of a constant force is therefore:
W = F · r . |
(10.60) |
|
|
Fig. 10.11 Work done by a constant force F as the object is moved from point 0 at r(t0) to point 1 at r(t1)
F r(t) F
F F
1
0
y |
r ’(t) |
x