Skip to main content

Planar Rigid-Body Motion

Planar rigid-body motion combines translation of a body with rotation about an axis perpendicular to the plane. A wheel rolling, a link rotating, a sliding block connected to a crank, and a plate moving in a mechanism all require this topic. The kinematics describe velocity and acceleration relationships across the body; the kinetics relate forces and moments to the motion.

The central benefit of the rigid-body model is that distances between points remain fixed. Once the motion of one point and the angular motion are known, the velocity and acceleration of every other point can be related by vector equations. Once the acceleration of the mass center and angular acceleration are known, planar force and moment balance determine or check the required forces.

Definitions

A rigid body in planar motion has angular velocity

ω=ωk\boldsymbol{\omega}=\omega\mathbf{k}

and angular acceleration

α=αk.\boldsymbol{\alpha}=\alpha\mathbf{k}.

For two points AA and BB fixed in the same rigid body, the relative velocity equation is

vB=vA+ω×rB/A.\mathbf{v}_B=\mathbf{v}_A+\boldsymbol{\omega}\times\mathbf{r}_{B/A}.

The relative acceleration equation is

aB=aA+α×rB/A+ω×(ω×rB/A).\mathbf{a}_B=\mathbf{a}_A+\boldsymbol{\alpha}\times\mathbf{r}_{B/A} +\boldsymbol{\omega}\times(\boldsymbol{\omega}\times\mathbf{r}_{B/A}).

The second added term is tangential relative acceleration; the third is normal relative acceleration toward point AA.

For a rigid body with mass center GG, planar kinetics are

Fx=maGx,\sum F_x=ma_{Gx}, Fy=maGy,\sum F_y=ma_{Gy}, MG=IGα.\sum M_G=I_G\alpha.

Moments may also be taken about an accelerating point, but extra terms are needed unless the point is fixed or the equation is written carefully. Taking moments about the mass center is usually safest.

For rolling without slipping between a wheel of radius RR and a fixed surface,

vG=Rω,v_G=R\omega, aG=Rαa_G=R\alpha

along the rolling direction when the radius is constant and the surface is fixed.

Key results

Rigid-body velocity fields are determined by one point velocity and angular velocity. In 2D, the cross product can be evaluated as

ω×(xi+yj)=ωyi+ωxj.\boldsymbol{\omega}\times(x\mathbf{i}+y\mathbf{j})= -\omega y\mathbf{i}+\omega x\mathbf{j}.

This makes it easy to write components:

vBx=vAxωyB/A,v_{Bx}=v_{Ax}-\omega y_{B/A}, vBy=vAy+ωxB/A.v_{By}=v_{Ay}+\omega x_{B/A}.

The instantaneous center of zero velocity is a point, possibly outside the body, whose velocity is zero at a given instant. If located, velocities of all points are perpendicular to lines from that center and have magnitudes v=ωrv=\omega r. This is a velocity tool only; the instantaneous center generally does not have zero acceleration.

Acceleration analysis cannot be done by using the instantaneous center as if it were a fixed pivot unless the center is actually fixed. The acceleration equation must include both tangential and normal relative terms:

aB/A=α×rB/Aω2rB/A.\mathbf{a}_{B/A}=\boldsymbol{\alpha}\times\mathbf{r}_{B/A} -\omega^2\mathbf{r}_{B/A}.

For planar kinetics, the mass center translates according to the resultant force, while the body rotates according to the resultant moment about the mass center. The equations are coupled by constraints. For example, in rolling without slipping, friction may be unknown but the constraint aG=Rαa_G=R\alpha links translation and rotation.

The planar rigid-body kinetic energy formula is

T=12mvG2+12IGω2.T=\frac{1}{2}mv_G^2+\frac{1}{2}I_G\omega^2.

This works for any planar rigid-body motion. If the body rotates about a fixed point OO, it may also be written as

T=12IOω2T=\frac{1}{2}I_O\omega^2

using the parallel-axis theorem, because vG=ωrG/Ov_G=\omega r_{G/O}.

Visual

SituationKinematic relationKinetic relation
Pure translationω=0\omega=0F=maG\sum\mathbf{F}=m\mathbf{a}_G, no rotational acceleration
Fixed-axis rotationv=rωv=r\omega, at=rαa_t=r\alpha, an=rω2a_n=r\omega^2MO=IOα\sum M_O=I_O\alpha if OO fixed
General planar motionvB=vA+ω×r\mathbf{v}_B=\mathbf{v}_A+\omega\times rMG=IGα\sum M_G=I_G\alpha
Rolling without slippingvG=Rωv_G=R\omega, aG=Rαa_G=R\alphaFriction supplies torque as needed

Worked example 1: Velocity and acceleration of a rotating bar

Problem. A rigid bar ABAB is 0.90.9 m long and rotates counterclockwise in the plane. At an instant, point AA has velocity vA=1.2i\mathbf{v}_A=1.2\mathbf{i} m/s and acceleration aA=0.5j\mathbf{a}_A=0.5\mathbf{j} m/s2^2. The vector from AA to BB is rB/A=0.9i\mathbf{r}_{B/A}=0.9\mathbf{i} m. The bar has ω=4\omega=4 rad/s counterclockwise and α=2\alpha=2 rad/s2^2 counterclockwise. Find vB\mathbf{v}_B and aB\mathbf{a}_B.

Method. Use relative velocity and acceleration equations.

  1. Angular velocity vector:
ω=4k.\boldsymbol{\omega}=4\mathbf{k}.
  1. Relative velocity:
ω×rB/A=4k×0.9i=3.6j.\boldsymbol{\omega}\times\mathbf{r}_{B/A} =4\mathbf{k}\times0.9\mathbf{i} =3.6\mathbf{j}.
  1. Velocity of BB:
vB=1.2i+3.6j m/s.\mathbf{v}_B=1.2\mathbf{i}+3.6\mathbf{j}\ \text{m/s}.
  1. Tangential relative acceleration:
α×rB/A=2k×0.9i=1.8j.\boldsymbol{\alpha}\times\mathbf{r}_{B/A} =2\mathbf{k}\times0.9\mathbf{i} =1.8\mathbf{j}.
  1. Normal relative acceleration:
ω×(ω×rB/A)=ω2rB/A=(42)(0.9i)=14.4i.\boldsymbol{\omega}\times(\boldsymbol{\omega}\times\mathbf{r}_{B/A}) =-\omega^2\mathbf{r}_{B/A} =-(4^2)(0.9\mathbf{i}) =-14.4\mathbf{i}.
  1. Acceleration of BB:
aB=aA+1.8j14.4i.\mathbf{a}_B=\mathbf{a}_A+1.8\mathbf{j}-14.4\mathbf{i}. aB=14.4i+(0.5+1.8)j.\mathbf{a}_B=-14.4\mathbf{i}+(0.5+1.8)\mathbf{j}. aB=14.4i+2.3j m/s2.\mathbf{a}_B=-14.4\mathbf{i}+2.3\mathbf{j}\ \text{m/s}^2.

The checked answer is

vB=1.2i+3.6j m/s,aB=14.4i+2.3j m/s2.\boxed{\mathbf{v}_B=1.2\mathbf{i}+3.6\mathbf{j}\ \text{m/s},\quad \mathbf{a}_B=-14.4\mathbf{i}+2.3\mathbf{j}\ \text{m/s}^2.}

The large negative xx acceleration is the centripetal part toward point AA.

Worked example 2: Rolling cylinder pulled by a horizontal force

Problem. A solid cylinder of mass 1212 kg and radius 0.200.20 m rolls without slipping on a horizontal surface. A horizontal force P=30P=30 N is applied at its center to the right. Find the acceleration of the center and the friction force. For a solid cylinder, IG=12mR2I_G=\frac{1}{2}mR^2.

Method. Draw forces: PP right at the center, friction FF at the ground, normal NN, and weight mgmg. Use translation, rotation about GG, and rolling constraint.

  1. Horizontal translation:

Assume friction positive to the right:

Fx=maG:P+F=ma.\sum F_x=ma_G:\quad P+F=ma.
  1. Rotation about GG. Force PP passes through GG, so it creates no moment. Friction at the bottom creates the torque. With rightward rolling requiring clockwise angular acceleration, take clockwise positive for rotation. If FF acts left, it creates clockwise torque. The algebra can instead assume friction left with magnitude ff:
Pf=ma.P-f=ma.

Moment about GG:

fR=IGα.fR=I_G\alpha.
  1. Rolling constraint:
a=Rα.a=R\alpha.
  1. Substitute α=a/R\alpha=a/R into moment equation:
fR=IGaR.fR=I_G\frac{a}{R}. f=IGR2a.f=\frac{I_G}{R^2}a.

For a solid cylinder,

IG=12mR2,I_G=\frac{1}{2}mR^2,

so

f=12ma.f=\frac{1}{2}ma.
  1. Substitute in translation:
P12ma=ma.P-\frac{1}{2}ma=ma. P=32ma.P=\frac{3}{2}ma. a=2P3m=2(30)3(12)=1.667 m/s2.a=\frac{2P}{3m}=\frac{2(30)}{3(12)}=1.667\ \text{m/s}^2.
  1. Friction magnitude:
f=12ma=12(12)(1.667)=10.0 N.f=\frac{1}{2}ma=\frac{1}{2}(12)(1.667)=10.0\ \text{N}.

Because we assumed ff left, the friction force is 1010 N left. The checked answer is

aG=1.67 m/s2 right,f=10.0 N left.\boxed{a_G=1.67\ \text{m/s}^2\ \text{right},\qquad f=10.0\ \text{N left}.}

Friction acts left even though the cylinder accelerates right because friction supplies the clockwise torque needed for rolling.

Code

import numpy as np

def planar_cross_omega_r(omega, r):
x, y = r
return np.array([-omega * y, omega * x])

vA = np.array([1.2, 0.0])
aA = np.array([0.0, 0.5])
rBA = np.array([0.9, 0.0])
omega = 4.0
alpha = 2.0

vB = vA + planar_cross_omega_r(omega, rBA)
aB = aA + planar_cross_omega_r(alpha, rBA) - omega**2 * rBA
print("vB =", vB)
print("aB =", aB)

m = 12.0
P = 30.0
a = 2.0 * P / (3.0 * m)
f = 0.5 * m * a
print(f"rolling acceleration = {a:.3f} m/s^2")
print(f"friction magnitude = {f:.2f} N left")

Common pitfalls

  • Using the instantaneous center to compute acceleration as if it were a fixed point.
  • Dropping the normal relative acceleration term ω2r-\omega^2\mathbf{r}.
  • Confusing angular velocity sign with angular acceleration sign.
  • Taking moments about a moving point without accounting for extra terms.
  • Assuming friction always opposes the center's motion in rolling problems.
  • Forgetting the rolling constraint changes sign depending on the chosen positive angular direction.
  • Using IGI_G when the body rotates about a fixed point and the correct equation needs IOI_O or a moment about GG.

Connections