- •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
15.3 Calculating the Moment of Inertia
Fig. 15.9 A thin plate with sides a and b and thickness h. The coordinate system is placed at the center of mass, with the axis z in the direction of the thickness of the plate
h
469
z
y
ρ
O
x b
a
Since the place is homogeneous, the mass density is uniform— constant in space. The distance ρ from the axis to the point r is simply the length of the horizontal (x y-plane) component of r:
ρ2 = x 2 + y2 . |
(15.29) |
The moment of inertia is found from the integral:
Iz = ρM |
−a/2 |
−b/2 |
0 |
x 2 + y2 d z d y d x = ρM h |
−a/2 |
−b/2 x 2 |
+ y2 d y d x , |
|
a/2 |
b/2 |
h |
|
a/2 |
b/2 |
|
|
|
|
|
|
|
|
(15.30) |
where we integrate the two parts of the sum inidividually:
Iz = ρM h −a/2 −b/2 x 2 d y d x + ρM h |
−a/2 |
−b/2 y2 d y d x |
|
||||||||||||||||||||||
|
|
|
|
|
|
|
a/2 |
|
|
b/2 |
|
|
|
|
|
a/2 |
|
|
b/2 |
|
|
|
|
||
|
|
= ρM hb |
−a/2 x 2d x + ρM h a −b/2 y2 d y |
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
|
|
a/2 |
|
|
|
|
|
|
b/2 |
|
|
|
|
− |
|
|
|
|
||
|
|
= |
|
M |
3 2 − 2 |
+ 2 |
|
2 |
|
||||||||||||||||
|
|
|
|
|
1 |
|
|
a |
3 |
−a |
3 |
|
|
b |
|
3 |
|
−b |
|
3 |
|
||||
|
|
|
ρ |
hb |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
, |
(15.31) |
|
which we simplify to: |
+ |
|
|
= |
|
|
|
a2 + b2 |
|
= |
|
|
a2 + b2 |
|
|||||||||||
Iz = |
3 |
ρM hba |
4 |
4 |
12 |
ρM V |
|
12 M |
, (15.32) |
||||||||||||||||
|
1 |
|
|
|
|
a2 |
|
|
b2 |
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
where we used that V = abh and M = ρM V .
15.4 Conservation of Energy for Rigid Bodies
We have found the kinetic energy of a rigid body that is either rotating around a fixed axis or that is rotating around its moving center of mass. If the solid body is only affected by conservative forces, we may use energy conservation principles to
470 |
15 Rotation of Rigid Bodies |
determine the relation between angular velocity and position, just as we have done with translational velocity and position previously. But in order to employ energy conservation, we need to determine the potential energy of a rigid body.
Potential Energy for a Constant Gravitational Force
What is the potential energy for a rigid body due to the gravitational force? We found that for a point particle with mass mi at the position ri , the potential energy due to a constant gravitational force directed along the y-axis is:
Ui = mi gyi . |
(15.33) |
We often call a constant gravitational force a homogeneous gravity, since gravity is the same everywhere. The total potential energy for multiparticle system of N particles is therefore:
NN
|
|
|
|
||||
U = |
Ui = |
mi gyi = g mi yi = M g Y , |
(15.34) |
||||
i =1 |
i =1 |
i =1 |
|
||||
|
|
|
= |
|
|
|
|
|
|
|
|
M Y |
|
where Y is the position of the center of mass. This result is general, and is also valid for a rigid body:
Potential energy of a rigid body for a constant gravitational force
U = M g Y . |
(15.35) |
Potential Energy Due to a Spring Force
A rigid body may also be affected by a force modelled by a spring force. This force may have several origins. It may represent a contact force due to an elastic contact; a contact force due to a rope or string attachment; an approximation for a more complex force acting on a single part of the rigid body; or as a force acting on all parts of the rigid body.
First, we start from the simplest case, where a spring is acting on a particular point on the rigid body. We illustate this situation by the contact between a sphere and an elastic floor for example while the rod is bouncing off the floor (see Fig. 15.10). We can then model the force from the floor on the sphere using a spring model, where
15.4 Conservation of Energy for Rigid Bodies
y
F
F
yf |
|
|
P |
P |
|
|
x |
|
|
|
Q
θ
P
R
471
O
D
Fig. 15.10 Left Illustration of contact between a sphere and a surface. Right Illustration of a rod that can rotate around point O, but it is attached with a spring at point P. The other end of the spring is attached at the point Q
we assume that the force depends on the position of a particular point, the point P , on the sphere. If the surface is horizontal at y = y f , the spring models a vertical normal force from the floor:
F = −k( py − y f ) , |
(15.36) |
when the sphere is in contact with the floor, that is when py < y f , where py is the vertical position of the point P . The potential energy of the sphere due to this spring force is then:
U ( py ) = |
1 |
py − y f 2 . |
(15.37) |
2 k |
For the sphere, we see that the position of point P is related to the center of mass of the sphere, py = Y − R, where Y is the vertical position of the center of mass, and R is the radius of the sphere.
This result is not only valid for a sphere, but for another body affected by a spring force. The potential energy of a spring force is related to the elongation of the spring. We therefore need to relate the motion of the rigid body to the elongation of the spring by relating for example the position of the center of mass of the body or the rotation angle θ of the body to the elongation of the spring. This is illustration by the rod on the right in Fig. 15.10: The rod can rotate freely around the point O, but is affected by a spring attached to the rod at the point P and to the wall at the point Q, so that when the rod is parallel to the wall, corresponding to θ = 0, the spring is in its equilibrium position. The potential energy in the spring then depends on the distance between P and Q, which for small rotation angles θ can be approximated as Q P Rθ , where R is the distance from P to O. The potential energy due to
the spring force affecting the rotating rod is therefore: |
|
||||||
Us (θ ) = |
1 |
k ( |
Q P)2 = |
1 |
k R2θ 2 . |
(15.38) |
|
|
|
|
|||||
2 |
2 |
472 |
15 Rotation of Rigid Bodies |
If the rod also is subject to gravity, we must also add the potential energy due to gravity:
Ug (θ ) = −Mg D sin θ , |
(15.39) |
where D is the distance from O to the center of mass of the rod. The total potential energy of the rod is then U = Us + Ug .
If you want to use a spring force model for a rigid body, you must therefore consider how to describe it in detail. We will introduce several ways to model such interactions in the following, but generally you are left to your own ingenuity.
Conservation of Energy for a Rigid Body
How can we now put the pieces together and use conservation of energy of a rigid body relate its velocity (angular and translational) to its position (rotational and translational)? As long as the system is only subject to conservative forces, the total energy of a rigid body is:
E = K + U = |
1 |
M V 2 + |
1 |
Icm ω2 + U . |
(15.40) |
|
|
||||
2 |
2 |
|
|
For a rigid body the potential energy can only depend on the translational and rotational motion of the body and not on the internal deformation. We have therefore only included one term for the potential energy: The potential energy due to external forces.
If the object is rotating around a fixed axis—an axis that does not move—we can simplify the kinetic energy, getting
|
1 |
|
E = |
2 IO ω2 + U . |
(15.41) |
These two equations ((15.40) and (15.41)) forms the basis for using energy conservation to solve problems with rigid bodies.
15.4.1 Example: Rotating Rod
Problem: Find the angular velocity for a rod that rotates without friction about an attachment point at one of its ends. The rod starts in a horizontal position. You can neglect air resistance.
15.4 Conservation of Energy for Rigid Bodies |
473 |
Approach: We use an energy argument to relate the angular velocity of the rod to its position based on the kinetic energy of a rotating rigid body and the potential energy of a solid body affected by gravity.
Solution: How can we address this problem with what we have learned so far? Can we use Newton’s laws of motion? Not really, because we do not know the forces acting in the attachment point. Normal forces will be present in this point—but we do not have a model to determine their magnitude.
What about energy conservation? We know how the kinetic energy of the rod depends on the angular velocity, and we know the potential energy of the rod in the gravity field. But is the mechanical energy conserved in this case? There are no frictional forces and no air resistance. The only external forces acting are the normal force in the attachment point and gravity. The force acting on the attachment point does no work, because this point is not moving! And we know that the work done by gravity can be expressed as a potential energy. We can therefore use energy conservation to solve this problem!
Identify: We start from a sketch of the system in Fig. 15.11. The rod has a length L. Its position is given by the angle θ , and the rod starts at θ = 0. Our task is to find the angular velocity as a function of θ . We choose a positive rotational direction, which means that the angle θ in Fig. 15.11 is negative!
Kinetic energy: The kinetic energy of a rigid rod rotating around the z-axis through O is given as:
|
1 |
|
K = |
2 Iz ω2 , |
(15.42) |
where Iz is the moment of inertia for the rod around the axis z through O. For a rod rotating around a fixed axis, this is the only term for the kinetic energy.
Caution: Notice that for an object rotating around a fixed axis you use the term K = (1/2)I ω2 for the kinetic energy, where ω is the angular velocity around the fixed axis. If you instead want to seperate the kinetic energy into the motion of the
Fig. 15.11 A thin rod of length L and mass M attached at the point O. The orientation of the rod is given by the angle θ . The position of the center of mass of the rod is found from the geometry to be
Y = −(L/2) sin θ
|
y |
|
|
x |
|
O |
cm |
|
θ |
||
|
Y=(L/2)sinθ
cm
474 |
15 Rotation of Rigid Bodies |
center of mass and the rotation around the center of mass, you have to figure out what the angular velocity around the center of mass is, which may not be the same as the angular velocity around O. A frequent mistake is the confuse the two cases of rotation around a fixed axis and rotation around the center of mass.
Potential energy: From (15.35) we found that the potential energy of a solid body due to gravity is
U = M g Y , |
(15.43) |
where Y is the vertical position of the center of mass of the object. For a homogenous rod, the center of mass is located at its geometric center, that is at a distance L/2 from the end. What is the vertical coordinate of this position? From Fig. 15.11 we see that of the center of mass is:
Y = |
L |
sin θ . |
(15.44) |
2 |
This is a negative number as long as θ is negative (and larger than −π ). The potential energy of the rod is therefore:
U = −Mg |
L |
sin θ . |
(15.45) |
2 |
Solve: The total energy of the rod is conserved since all the forces are conservative (or not doing any work on the body), the total energy, E = K + U , is conserved. The total energy is therefore the same in the initial position, 0, when θ = θ0 = 0, and in the position 1, when the angle is θ :
|
|
|
|
|
E0 = K0 + U0 = E1 = K1 + U1 , |
|
|
(15.46) |
||||||||||||||||||
where |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
E0 = K0 + U0 = |
1 |
Iz ω02 + MgY0 = 0 , |
|
|
(15.47) |
||||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||||||||||
|
|
|
|
|
2 |
|
|
|||||||||||||||||||
and |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
1 |
Iz ω12 + MgY1 |
|
|
|
1 |
Iz ω2 − |
|
L |
|
|
|
|
|
||||||||
|
|
E1 = K1 + U1 = |
|
= |
|
|
|
|
Mg sin θ . |
|
|
(15.48) |
||||||||||||||
2 |
2 |
2 |
|
|
||||||||||||||||||||||
Energy conservation, E0 = E1, gives: |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
2 Iz ω2 |
− |
2 |
Mg sin θ = 0 |
2 Iz ω2 |
= |
2 |
|
Mg sin θ |
ω = ± |
|
|
|
|
|||||||||||||
|
|
Iz |
sin θ , |
|||||||||||||||||||||||
1 |
|
|
L |
|
|
|
1 |
|
|
|
|
|
L |
|
|
|
|
|
|
|
Mg L |
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
(15.49) |