Skip to article frontmatterSkip to article content
Site not loading correctly?

This may be due to an incorrect BASE_URL configuration. See the MyST Documentation for reference.

Newton's laws

Newton’s second law for 1D motion

Let’s first derive Newton’s second law for a set of particles, and then we’ll specialize it to 1D motion. Each particle comprising the rigid body satisfies the standard version of Newton’s second law for a particle:

mir¨i=Fiext+jFij,m_i \ddot \r_i = \F^\textsf{ext}_i + \sum_j \F_{ij},

where Fiext\F^\textsf{ext}_i is the total external force applied to particle ii, and Fij\F_{ij} is the internal force applied to particle ii from particle jj.

Define the center of mass as rcm:=1mimiri\r_\textsf{cm} := \frac{1}{m}\sum_i m_i \r_i, where m:=imim := \sum_i m_i is the total mass. differentiating twice and plugging in (1), we obtain

mr¨cm=imir¨i=iFiext+ijFij0m\,\ddot \r_\textsf{cm} = \sum_i m_i \ddot\r_i = \sum_i \F^\textsf{ext}_i + \underbrace{\sum_i\sum_j \F_{ij}}_0

The internal forces cancel because of Newton’s third law: Fij=Fji\F_{ij} = -\F_{ji}. We therefore obtain Newton’s second law for a system of particles:

mr¨cm=iFiext\boxed{m\,\ddot \r_\textsf{cm} = \sum_i \F^\textsf{ext}_i}

In other words, no matter what the particles are doing (translating, rotating, moving apart), the center of mass behaves like a particle of mass mm with all external forces acting at the center of mass. This is particularly useful when we have gravity acting on a rigid body, for example; we can just assume the gravity is acting at the center of mass with a force equal to the total gravity force mgmg.

Special case: 1D rigid body motion with no rotation. In this case, all particles of the body have the same velocity and acceleration, and it is a scalar quantity due to the 1D nature of the motion. So if x(t)x(t) is the position of any point on the rigid body, we have:

mx¨=iFiext.\boxed{ m\,\ddot x = \sum_i F_i^\textsf{ext}}.

Newton’s second law for planar rotation

Let’s first derive a more general result for rotational motion of rigid bodies, and then we’ll specialize it to rotation about a fixed axis. As with translational motion, we can think of the rigid body as a collection of particles with positions ri\r_i and masses mim_i. Each particle satisfies Newton’s second law:

mir¨i=Fiext+jFij,m_i \ddot \r_i = \F^\textsf{ext}_i + \sum_j \F_{ij},

where Fiext\F^\textsf{ext}_i is the external force on particle ii and Fij\F_{ij} is the internal force on particle ii from particle jj. To extract rotational dynamics, we take the moment (torque) of each particle’s equation of motion about a reference point rref\r_\textsf{ref}. We will assume an arbitrary (possibly moving) reference point for now. Taking the cross product with (rirref)(\r_i - \r_\textsf{ref}), position with respect to the reference point, we obtain

(rirref)×mir¨i=(rirref)×Fiext+j(rirref)×Fij.(\r_i - \r_\textsf{ref}) \times m_i \ddot{\r}_i = (\r_i - \r_\textsf{ref}) \times \F_i^{\textsf{ext}} + \sum_j (\r_i - \r_\textsf{ref}) \times \F_{ij}.

If we sum over all particles ii, the internal forces cancel as in the linear case. However, the reason is more subtle this time. Newton’s third law says that Fji=Fij\F_{ji} = -\F_{ij}, so we have:

(rirref)×Fij+(rjrref)×Fji=(rirref)×Fij(rjrref)×Fij=(rirj)×Fij=0\begin{aligned} &\hspace{-1cm} (\r_i - \r_\textsf{ref}) \times \F_{ij} + (\r_j - \r_\textsf{ref}) \times \F_{ji} \\ &= (\r_i - \r_\textsf{ref}) \times \F_{ij} - (\r_j - \r_\textsf{ref}) \times \F_{ij} \\ &= (\r_i-\r_j) \times \F_{ij} \\ &= 0 \end{aligned}

The last step is true because we assume the force between particles ii and jj acts along the vector joining particles ii and jj (also known as a central force). Therefore the vector rirj\r_i-\r_j is parallel to the vector Fij\F_{ij} and so their cross product is zero. This assumption is valid for ordinary rigid bodies and classical force interactions. Systems that transmit internal moments (e.g., via friction or actuation) require additional torque terms.

Summing over all particles ii in Eq. (6) therefore gives

i(rirref)×mir¨i=i(rirref)×FiextTiext,\sum_i (\r_i - \r_\textsf{ref}) \times m_i \ddot{\r}_i = \sum_i \underbrace{(\r_i - \r_\textsf{ref}) \times \F_i^{\textsf{ext}}}_{\T_i^\textsf{ext}},

where we defined Tiext\T_i^\textsf{ext} as the torque about the reference point caused by the external force Fiext\F_i^\textsf{ext}. Let us also define the angular momentum of a particle measured with respect to the reference point:

Li:=(rirref)×mir˙i\L_i := (\r_i - \r_\textsf{ref}) \times m_i \dot\r_i

Taking the time derivative and applying the product rule, we have:

L˙i=(r˙ir˙ref)×mir˙i+(rirref)×mir¨i\dot\L_i = (\dot\r_i -\dot\r_\textsf{ref}) \times m_i \dot\r_i + (\r_i-\r_\textsf{ref}) \times m_i \ddot\r_i

Rearrange and use the fact that the cross-product of parallel vectors is zero:

(rirref)×mir¨i=L˙i+r˙ref×mir˙i(\r_i-\r_\textsf{ref}) \times m_i \ddot\r_i = \dot\L_i \,+\, \dot\r_\textsf{ref} \times m_i \dot\r_i

Substitute (11) into (8) and obtain

L˙+r˙ref×imir˙i=iTiext,\dot\L \,+\, \dot\r_\textsf{ref} \times \sum_i m_i \dot\r_i = \sum_i \T_i^\textsf{ext},

where we defined the total angular momentum L:=iLi\L := \sum_i\L_i. Recall our definition of the center of mass rcm:=1mimiri\r_\textsf{cm} := \frac{1}{m} \sum_i m_i \r_i, so we can rewrite (12) as

L˙+r˙ref×mr˙cm=iTiext\dot\L \,+\, \dot\r_\textsf{ref} \times m \dot\r_\textsf{cm} = \sum_i \T_i^\textsf{ext}

The extra term r˙ref×mr˙cm \dot\r_\textsf{ref} \times m \dot\r_\textsf{cm} will vanish if we choose our reference rref\r_\textsf{ref} to be either a fixed point (zero velocity), or the center of mass (so the cross product is zero). In either of these cases, we obtain Newton’s second law for rotational motion of a rigid body:

L˙=iTiext\boxed{\dot\L = \sum_i \T_i^\textsf{ext}}

So far, we assumed rref\r_\textsf{ref} is either fixed in space or it’s the center of mass of the body. We want the moment of inertia of the body with respect to our reference to be constant, so let’s further assume that our reference is fixed to the body. So we are now assuming that

  1. rref\r_\textsf{ref} is fixed to the body and fixed in space (for example, if the object is pinned at a point and rotates about that point), or

  2. rref\r_\textsf{ref} is the center of mass of the object.

Since the body is rigid, the distance between any pair of points fixed to the body is constant. This means that any point ri\r_i can only rotate with respect to rref\r_\textsf{ref}, so there must exist some vector ω\bomega such that for any particle ii,

r˙i=r˙ref+ω×(rirref)\dot\r_i = \dot\r_\textsf{ref} + \bomega \times (\r_i - \r_\textsf{ref})

substituting this into (14) with the definition (9), we obtain

L=i(rirref)×mir˙i=i(rirref)×mi(r˙ref+ω×(rirref))=i(rirref)×mi(ω×(rirref))+m(rcmrref)×r˙ref0\begin{aligned} L &= \sum_i (\r_i-\r_\textsf{ref}) \times m_i \dot\r_i \\ &= \sum_i (\r_i-\r_\textsf{ref}) \times m_i \bigl( \dot\r_\textsf{ref} + \bomega \times (\r_i - \r_\textsf{ref}) \bigr) \\ &= \sum_i (\r_i-\r_\textsf{ref}) \times m_i \bigl( \bomega \times (\r_i - \r_\textsf{ref}) \bigr) + \underbrace{m(\r_\textsf{cm}-\r_\textsf{ref})\times \dot\r_\textsf{ref}}_0 \end{aligned}

The rightmost term is zero because either r˙ref=0\dot\r_\textsf{ref}=0 (fixed reference), or rref=rcm\r_\textsf{ref}=\r_\textsf{cm}. The remainder of the expression is a linear function of ω\bomega, so we can write it as L=Iω\L = \I \bomega, where I\I is called the inertia tensor.

In general, I\I can be a function of time, so when we plug in L=Iω\L = \I\bomega into (14), we need to use the product rule. This doesn’t happen in our particular setting, because we have a rigid body rotating around a fixed axis.

Suppose the axis of rotation is k^\khat and the body rotates with angular speed θ˙(t)\dot\theta(t). We therefore have ω=θ˙k^\bomega = \dot\theta\,\khat. We only care about the component of L\L in the k^\khat direction, so we have:

k^L=k^imiθ˙(rirref)×(k^×(rirref))=imik^(rirref)2J  θ˙\begin{aligned} \khat \cdot \L &= \khat\cdot \sum_i m_i \dot\theta\, (\r_i-\r_\textsf{ref}) \times \bigl( \khat \times (\r_i - \r_\textsf{ref}) \bigr) \\ &= \underbrace{\sum_i m_i \norm{ \khat \cdot (\r_i-\r_\textsf{ref})}^2}_J \; \dot\theta \end{aligned}

where we used the cyclic property of the triple product to simplify the expression:

a(b×c)=b(c×a)=c(a×b)\a \cdot (\b \times \c) = \b \cdot (\c \times \a) = \c \cdot (\a \times \b)

So we have k^L=Jθ˙\khat \cdot \L = J\dot\theta, where JJ is the moment of inertia of our body about an axis parallel to our axis of rotation and passing through rref\r_\textsf{ref}. Substituting into (14), we obtain Newton’s second law for a rigid body rotating about a fixed axis, where the reference point is either the axis or the center of mass.

Jθ¨=iTiext\boxed{ J\ddot \theta = \sum_i T_i^\textsf{ext}}