Rigid Body Dynamics Article Index for
Rigid 		Website Links For
Body 		 

Information About

Rigid Body Dynamics
  CATEGORIES ABOUT RIGID BODY DYNAMICS
   
rigid bodies
   
rotational symmetry
   
APPAREL
BABY
BEAUTY
BOOKS
CAR TOYS
CELL PHONES
DVD'S
ELECTRONICS
GOURMET FOOD
GROCERIES
HEALTH & PERSONAL
HOME & GARDEN
JEWELRY
MUSIC
MUSIC INSTRUMENTS
OFFICE PRODUCTS
SOFTWARE
SPORTING GOODS
TOOLS & HARDWARE
TOYS
VIDEO GAMES
SHOPPING HOME
MORE SHOPPING...



In Physics , Rigid Body Dynamics differs from Particle Dynamics in that the body takes up Space and can Rotate , which introduces other considerations.

Note: This article has much overlap with the Rigid Rotor and Rigid Body articles. Articles should eventually be merged.


Equations from particle dynamics can be generalized to rigid body dynamics as follows:


RIGID BODY LINEAR MOMENTUM


The equation for particle Linear Momentum is

: rac{\mathrm{d}(m v)}{\mathrm{d}t}=\sum_{i=1}^N f_i

where:
  • ''m'' is the particle's mass.

  • ''v'' is the particle's velocity.

  • ''f''''i'' is one of the ''N'' forces acting on the particle.


Assuming constant mass, this reduces to

:m rac{\mathrm{d}v}{\mathrm{d}t}=\sum_{i=1}^N f_i.

To generalize assume a body of finite mass and size is composed of such particles. There exist internal forces, acting between any two particles, and external forces, acting only on the outside of the mass. Each particle has:
  • a mass \mathrm{d}m.

  • a position vector ''r''.


Thus, the linear momentum equation of any given particle would look like this:

:\mathrm{d}m rac{\mathrm{d}^2r}{\mathrm{d}t^2}= \sum_{i=1}^M f_{i, ext{internal}} + \sum_{j=1}^N f_{j,\mathrm{external}}.

If the equation for each particle were added together, the internal forces would cancel out, since by Newton's third law, any such force would have opposite magnitudes on the two particles. Also, the left side would become an integral over the entire body, and the second derivative operator could come out of the integral, leaving

: rac{\mathrm{d}^2}{\mathrm{d}t^2} \int r\, \mathrm{d}m = \sum_{j=1}^N f_{j,\mathrm{external}}.

Letting ''M'' be the total mass, the left side can be multiplied and divided by M without changing the validity:

: M rac{\mathrm{d}^2 rac{\int r\, \mathrm{d}m}{M}}{\mathrm{d}t^2} = \sum_{j=1}^N f_{j,\mathrm{external}}

However, rac{\int r\, \mathrm{d}m}{M} is the formula for the position of Center Of Mass . Denoting this by ''r''cm, the equation reduces to

: M rac{\mathrm{d}^2 r_{cm}}{\mathrm{d}t^2} = \sum_{j=1}^N f_{j,\mathrm{external}}.

Thus, linear momentum equations can be extended to rigid bodies by denoting that they describe the motion of the center of mass of the body.


RIGID BODY ANGULAR MOMENTUM


The most general equation for rotation of a rigid body in three dimensions about an arbitrary origin ''O'' with axes ''x'', ''y'', ''z'' is

:M b_{G/O} imes rac{\mathrm{d}^2 R_O}{\mathrm{d}t^2} + rac{\mathrm{d}(\mathbf{I}\boldsymbol{\omega})}{\mathrm{d}t} = \sum_{j=1}^N au_{O,j}

where:
:
\mathbf{I} =
\begin{pmatrix}
\int (y^2+z^2)\, \mathrm{d}m & -\int xy\, \mathrm{d}m & -\int xz\, \mathrm{d}m\ -\int xy\, \mathrm{d}m & \int (x^2+z^2)\, \mathrm{d}m & -\int yz\, \mathrm{d}m \ -\int xz\, \mathrm{d}m & -\int yz\, \mathrm{d}m & \int (x^2+y^2)\, \mathrm{d}m
\end{pmatrix}
\quad \hbox{and} \quad \boldsymbol{\omega} =
\begin{pmatrix} \omega_x \ \omega_y \ \omega_z
\end{pmatrix}.

Here \mathbf{I} is the Moment Of Inertia Tensor and \boldsymbol{\omega} is the Angular Velocity (a vector). Based on this, a theorem states that any rigid body is equivalent when moving to a Poinsot's Ellipsoid .

Further
  • ωq is the angular velocity about axis ''q''.

  • ''M'' is the total mass.

  • ''b''''G''/''O'' is the vector from ''O'' to the body's center of mass.

  • ''R''''O'' is the position of ''O''.

  • ''t'' is time.

  • \int \, \mathrm{d}m is an integral over the mass of the body.


  • au_{O,j} is one of the ''N'' moments about ''O''.


This equation follows from equation for linear momentum of a particle and kinematics; no additional observations of nature are necessary to arrive at it.

There are many special cases that simplify this equation. The first term goes to zero if any of three conditions are met:

  • O is a fixed point (since its second derivative would be zero).

  • A set of axes is chosen with its origin attached to the body's center of mass (since this would reduce the vector b to zero).

  • The vector b always points in the direction of the acceleration of O (since the cross product of parallel vectors is zero).


Also, if the axes are chosen are the principal axes (i.e., the moments about the ''xy'', ''xz'', and ''yz'' planes is zero), the off-diagonal terms of the matrix are zero. This case is further discussed by Euler's Equations .

When learning about angular motion, students are generally first exposed to the case of rotation only in the x-y plane and a fixed axis or axis at the center of mass with constant rotational inertia. That equation is

:\int \left(x^2+y^2 ight)\, \mathrm{d}m rac{\mathrm{d}\omega_z}{\mathrm{d}t} = \sum_{j=1}^N au_{O,j,z}.


ANGULAR MOMENTUM AND TORQUE


Similarly, the Angular Momentum \mathbf{L} for a system of particles with linear momenta p_{i} and distances r_{i} from the rotation axis is defined

:
\mathbf{L} = \sum_{i=1}^{N} \mathbf{r}_{i} imes \mathbf{p}_{i} =
\sum_{i=1}^{N} m_{i} \mathbf{r}_{i} imes \mathbf{v}_{i}


For a rigid body rotating with angular velocity \omega about the rotation axis \mathbf{\hat{n}} (a Unit Vector ), the velocity vector \mathbf{v}_{i} may be written as a Vector Cross Product

:
\mathbf{v}_{i} = \omega \mathbf{\hat{n}} imes \mathbf{r}_{i} \ \stackrel{\mathrm{def}}{=}\ \boldsymbol\omega imes \mathbf{r}_{i}


where
:angular velocity vector \boldsymbol\omega \ \stackrel{\mathrm{def}}{=}\ \omega \mathbf{\hat{n}}
:\mathbf{r}_{i} is the shortest vector from the rotation axis to the point mass.

Substituting the formula for \mathbf{v}_{i} into the definition of \mathbf{L} yields

:
\mathbf{L} =
\sum_{i=1}^{N} m_{i} \mathbf{r}_{i} imes (\boldsymbol\omega imes \mathbf{r}_{i}) =
\boldsymbol\omega \sum_{i=1}^{N} m_{i} r_{i}^{2} =
I \omega \mathbf{\hat{n}}


where we have introduced the special case that the position vectors of all particles are perpendicular to the rotation axis (''e.g.,'' a Flywheel ): \boldsymbol\omega \cdot \mathbf{r}_{i} = 0.

The Torque \mathbf{N} is defined as the rate of change of the angular momentum \mathbf{L}

:
\mathbf{N} \ \stackrel{\mathrm{def}}{=}\ rac{d\mathbf{L}}{dt}


If I is constant (because the inertia tensor is the identity, because we work in the intrinsecal frame, or because the torque is driving the rotation around the same axis \mathbf{\hat{n}} so that \mathrm{I} is not changing) then we may write

:
\mathbf{N} \ \stackrel{\mathrm{def}}{=}\ I rac{d\omega}{dt}\mathbf{\hat{n}} =
I \alpha \mathbf{\hat{n}}


where
:\alpha is called the ''angular acceleration'' (or ''rotational acceleration'') about the rotation axis \mathbf{\hat{n}}.

Notice that if I is not constant in the external reference frame (ie. the three main axes of the body are different) then we cannot take the I outside the derivate. In this cases we can have Torque-free Precession .


APPLICATIONS


Computer Physics Engine s use rigid body dynamics to increase interactivity and realism in Video Games .


SEE ALSO



EXTERNAL LINKS