- •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
146 |
6 |
Motion in Two and Three Dimensions |
|
r(t ) and F(t ) . |
(6.25) |
Numerically, we will have a corresponding sequence of positions or forces at discrete times, ti :
ri = r(ti ) and Fi = F(ti ) . |
(6.26) |
Fortunately, it is simple to both represent and apply mathematical operations to an element in a sequence.
We generate a sequence of n vectors ri with x , y, and z coordinates by:
n = 10
r = zeros((n,3),float)
We can use mathematical vector operations directly on element in the sequence, as illustrated in the following example:
v = array([1.0,-2.0,2.0]) n = 10
r = zeros((n,3),float) r[0] = array([0,0,0]); dt = 0.1
r[1] = r[0] + v*dt
6.2 Description of Motion
The cheetah is the world fastest land animal. How can we characterize its motion as it chases its prey? How fast does it run and how fast does it turn?
Motion Diagram and Position Vector
Figure 6.2 shows a few frames from a movie of a cheetah chasing a Thompson gazelle. To quantify the motion we generate a motion diagram: We mark the position of the cheetah at regular time intervals and record the positions r(ti ) of the cheetah relative to the origin at time ti .
We are free to choose the origin and the axes of the coordinate system. The origin determines where we measure the positions from. In Fig. 6.2 we have chosen a stationary point on the ground as the origin. In addition to the origin, a coordinate system consists of a set of axes that we use to decompose the position vector. The directions of the axes indicate the positive direction of the corresponding unit vector. The position can be decomposed along the x , y, and z-axes respectively (Table 6.1):
r(t ) = x (t ) i + y(t ) j + z(t ) k , |
(6.27) |
where x (t ), y(t ), and z(t ) are lengths along the axes and hence have units of length. For example, the position at t = 0.5 s is:
6.2 Description of Motion
z
y
|
40 |
|
|
|
30 |
|
|
[ m ] |
20 |
|
|
y |
|
|
|
|
0.0s |
|
0.5s |
|
10 |
|
|
|
0 |
|
|
|
0 |
10 |
20 |
147
x
2.5s 3.0s
2.0s |
3.5s |
|
1.5s
4.0s
1.0s
30 |
40 |
50 |
60 |
70 |
80 |
90 |
100 |
|
|
x [ m ] |
|
|
|
|
|
Fig. 6.2 Illustration of a cheetah chasing a Thompson gazelle, and an illustration of the twodimensional motion and the two-dimensional motion diagram. (Illustration by S.B. Skattum)
Table 6.1 Positions of the running cheetah
ti |
(s) |
0.0 |
0.5 |
1.0 |
1.5 |
2.0 |
2.5 |
3.0 |
3.5 |
xi |
(m) |
0.0 |
16.2 |
29.9 |
38.9 |
46.6 |
57.2 |
71.7 |
84.8 |
yi |
(m) |
15.0 |
14.95 |
19.0 |
26.1 |
31.3 |
35.4 |
35.7 |
30.0 |
r(1.0 s) = x (1.0 s) i + y(1.0 s) j = 29.9 m i + 19.0 m j , |
(6.28) |
where we have skipped the z-coordinate since all the motion is in the x y-plane. You find the complete dataset in the file cheetah.d.2 Each line gives ti , xi , yi , where ti is measured in seconds, and xi and yi are measured in meters. We have tabulated the positions at t = 0.5 s intervals in the following table:
2http://folk.uio.no/malthe/mechbook/cheetah.d.
148 |
6 Motion in Two and Three Dimensions |
Velocity Vector |
|
Figure 6.2 shows how the position changes over a time interval t . The change in position is also a vector and is called the displacement. The displacement from t = 1.0 s to t = 2.0 s is denoted r(1.0 s):
r(1.0 s) = r(2.0 s) − r(1.0 s) = (46.6 m i + 31.3 m j) − (29.9 m i + 19.0 m j) = 16.7 m i + 12.3 m j .
(6.29) We can read the displacement directly from the motion diagram in Fig. 6.2 as the vector from r(1.0 s) to r(2.0 s). Because the displacement depends on a difference between two positions, it does not depend on the choice of origin.
We see from Fig. 6.2 that both the length and the direction of the displacement vectors are changing throughout the motion. The rate of change of the displacement, the velocity, must therefore also be a vector:
The average velocity from a time t = t0 to a time t = t0 + |
t is defined as: |
|||||||
v(t |
|
) |
= |
r(t0 + t ) − r(t0) |
= |
r(t0) |
. |
(6.30) |
|
|
|
||||||
¯ |
0 |
|
t |
t |
|
It is measured in meters per second, m/s.
We find the average velocity at t = 1.0 s using the data in the table above:
v¯ (1.0 s) = |
r(1.0 s) |
= |
16.7 m |
i + |
12.3 m |
j = 16.7 m/s i + 12.3 m/s j . (6.31) |
|
|
|
|
|
||||
1.0 s |
1.0 s |
1.0 s |
If we instead use a time interval t = 0.5 s to find the average velocity at t = 1.0 s we find:
r(1.0 s) = r(1.5 s) − r(1.0 s) = (38.9 m i + 26.1 m j) − (29.9 m i + 19.0 m j) = 9.0 m i + 7.1 m j ,
|
|
|
|
|
(6.32) |
v¯ (1.0 s) = |
r(1.0 s) |
= |
1 |
(9.0 m i + 7.1 m j) = 18.0 m/s i |
+ 14.2 m/s j . |
|
|
||||
t |
|
||||
|
|
0.5 s |
(6.33) |
||
|
|
|
|
|
|
We see that the average velocity depends on the time interval |
t , just as we saw |
in the one-dimensional case. Again, we can understand this better by studying the displacement at t = 1.0 s for smaller and smaller time intervals t , as shown in Fig. 6.3. We see that as t becomes smaller, the displacement also becomes smaller, but its direction approaches the tangent to the curve describing the motion around t = 1.0 s.
6.2 Description of Motion
Fig. 6.3 Illustration of the average velocity, r¯(2.0s) for decreasing time intervals t for the motion of the cheetah
|
|
|
|
|
|
|
|
|
149 |
|
23 |
|
|
|
|
|
|
|
|
|
22 |
|
|
|
|
|
|
|
|
|
21 |
|
|
|
|
|
|
|
|
] |
20 |
|
|
|
|
|
|
|
|
[ m |
19 |
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
18 |
|
|
|
|
|
|
|
|
|
17 |
|
|
|
|
|
|
|
|
|
16 |
|
|
|
|
|
|
|
|
|
15 |
|
|
|
|
|
|
|
|
|
18 |
20 |
22 |
24 |
26 |
28 |
30 |
32 |
34 |
|
|
|
|
|
x [ m ] |
|
|
|
|
The instantaneous velocity at the time t is defined as the limit of the average velocity when the time interval t goes to zero, that is, the time derivative of the position vector, r(t ).
v(t ) |
= |
lim |
0 |
r(t ) |
= |
lim |
r(t + |
t ) − r(t ) |
dr . |
(6.34) |
||
|
t |
→ |
t |
t |
→ |
0 |
t |
= dt |
|
|||
|
|
|
|
|
|
|
|
|
|
|
Whenever we use the term velocity, we will mean the instantaneous velocity. The velocity vector is tangential to the trajectory.
Speed
The magnitude of the velocity vector is called the speed, v, defined as:
v(t ) = |v(t )| . |
(6.35) |
We use the word velocity for the velocity vector, and the word speed for the magnitude of the vector velocity.
Time Derivatives of Vector Functions
How do we find the derivative of a vector function such as r(t )? The simplest approach is to write the vector in terms of the unit vectors for the coordinate system, and then take the derivative of each component:
|
d |
d |
d x d y |
d z |
||||||
v(t ) = |
|
r(t ) = |
|
(x (t ) i + y(t ) j + z(t ) k) = |
|
i + |
|
j + |
|
k , (6.36) |
|
|
|
|
|
dt |
dt |
dt |
dt |
dt |
150 |
|
|
|
6 Motion in Two and Three Dimensions |
|||
where3 we define the component-wise velocities as |
|
||||||
d x |
|
d y |
d z |
|
|||
vx (t ) = |
|
, vy (t ) = |
|
, vz (t ) = |
|
. |
(6.37) |
dt |
dt |
dt |
The velocity vector can therefore also be written:
v(t ) = vx (t ) i + vy (t ) j + vz (t ) k = vx (t ), vy (t ), vz (t ) . |
(6.38) |
Acceleration
It is clear from Fig. 6.2 that the average velocity of the cheetah is not constant throughout the motion—it is varying both in direction and magnitude. Just as we introduced velocity to characterize the change in the position vector, we introduce the acceleration vector to characterize the change in the velocity vector:
The average acceleration vector over a time interval |
t from t to t + t is |
||
defined as: |
|
|
|
a |
v(t + t ) − v(t ) |
. |
(6.39) |
¯ = |
t |
|
We define the instantaneous acceleration vector, or simply the instantaneous acceleration, as the limit of the average acceleration vector when the time interval approaches zero:
a(t ) |
|
lim |
|
v(t + |
t ) − v(t ) |
|
dv |
v . |
(6.40) |
|
= |
0 |
|
t |
= dt |
||||||
|
t |
→ |
|
= ˙ |
|
|||||
|
|
|
|
|
|
|
|
|
|
The acceleration vector is the time derivative of the velocity vector.
We find the acceleration in the vector component representation by componentwise derivation:
|
dvx |
d |
|
|
|
|
|
z |
|
|
|
|
|
|
||
a(t ) = |
d |
|
|
|
+ vy j + vz k |
|
|
|
|
|||||||
dt |
v = |
dt |
|
vx i |
|
|
|
(6.41) |
||||||||
|
|
|
dvy |
|
|
dv |
|
|
|
|
|
|
|
|||
= dt i + |
|
|
j + |
|
k = ax (t ) i + ay (t ) j + az |
(t ) k , |
|
|||||||||
dt |
dt |
|
||||||||||||||
where we see that: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
dvx |
, ay (t ) = |
dvy |
, az (t ) = |
dvz |
|
|
||||
|
|
ax (t ) = |
|
|
|
. |
(6.42) |
|||||||||
|
|
dt |
dt |
dt |
3Here we have implicitly assumed that the time derivatives of the unit vectors are zero. This is not necessarily the case: The unit vectors vary with time for rotating reference systems.
6.2 Description of Motion |
151 |
||
Since the velocity vector is the time derivative of the position vector |
|
||
v(t ) = |
d |
|
|
|
r(t ) , |
(6.43) |
|
dt |
we can write the acceleration vector as the second time derivative of the position vector:
a(t ) = |
d |
v = |
d d |
r = |
d2r |
(6.44) |
|||
|
|
|
|
|
, |
||||
dt |
dt |
dt |
dt 2 |
which can be written on component form:
|
dvx |
|
d2 x |
dvy |
d2 y |
dvz |
d2z |
|
|||||
ax (t ) = |
|
= |
|
, ay (t ) = |
|
= |
|
, az (t ) = |
|
= |
|
. |
(6.45) |
dt |
dt 2 |
dt |
dt 2 |
dt |
dt 2 |
Notation for Time Derivatives
In physics, we use both the differential form, d/dt and the dot notation for time derivatives:
v = |
dr |
= r˙ , |
(6.46) |
||
|
|||||
dt |
|||||
and |
|
|
|
||
a = |
dv |
= v˙ = r¨ , |
(6.47) |
||
|
|||||
dt |
but we do not use the “marked” notation, vx = x ′(t ), often used in mathematics because we often use x ′(t ) to mean the position x measured in another coordinate system. We therefore strongly recommend you to use the notations introduced here, either the d/dt notation or the dot-notation.
Interpretation of Motion Diagrams
It is often takes time to gain a good intuition for acceleration, in particular for twoand three-dimensional motions. Motion diagrams can help in developing that intuition by visualizing velocities and accelerations.
If the motion diagram is drawn using a constant time interval t , we can use the displacement vector as a visualization of the velocity, since the velocity is proportional to the displacement:
152 |
6 Motion in Two and Three Dimensions |
40 |
|
|
|
|
|
|
|
|
|
|
30 |
|
|
|
|
|
|
|
|
|
|
[m ] |
|
|
|
|
|
|
|
|
|
|
20 |
|
|
|
|
|
|
|
|
|
|
y |
|
|
|
|
|
|
|
|
|
|
10 |
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
0 |
10 |
20 |
30 |
40 |
50 |
60 |
70 |
80 |
90 |
100 |
|
|
|
|
|
x [ m ] |
|
|
|
|
|
Fig. 6.4 Motion diagrams for the cheetah with t = 0.5 s illustrating both the displacements, interpreted as velocities, and the change in displacements, interpreted as accelerations. The constructions of the accelerations are illustrated
v¯(ti ) = |
1 |
r(ti ) . |
(6.48) |
t |
The displacement vectors are illustrated by red vectors at intervals of 0.5 and 0.25 s in Fig. 6.4.
Notice that if we want to look at the change in velocity at the time t = 1.0 s, we would like to compare the velocity before the time t = 1.0 s with the velocity after the time t = 1.0 s. Now, the average velocity at the time t = 1.0 s is really the average velocity over the time interval from 1.0 to 1.5 s. And the average velocity at the time t = 0.5 s is the average velocity over the time interval from 0.5 to 1.0 s. Therefore, a reasonable way to characterize the change in velocity at t = 1.0 s is to characterize it as the change in velocity over the time interval from 0.5 t to 1.5 s:
v¯(1.0 s) = v¯ (1.0 s) − v¯ (0.5 s) . |
(6.49) |
We can interpret this as the average acceleration of the cheetah at t = 1.0 s, since it is (approximately) proportional to the average acceleration:
1
a¯ (1.0 s) |
|
(v¯ (1.0 s) − v¯(0.5 s)) . |
(6.50) |
|
t
This method therefore provides a way to use the motion diagram to find approximations for the acceleration vectors in each point on the motion diagram. Incidentally, this method is the same as the simplest numerical method to find the second order time derivative of the position vector.