Peeter Joot's (OLD) Blog.

Math, physics, perl, and programming obscurity.

Rigid body motion (notes from an audited phy354 classical mechanics course (Taught by Prof. Erich Poppitz))

Posted by peeterjoot on March 8, 2012

[Click here for a PDF of this post with nicer formatting and figures if the post had any (especially if my latex to wordpress script has left FORMULA DOES NOT PARSE errors.)]

Rigid body motion.

Setup

We’ll consider either rigid bodies as in the connected by sticks figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig1}) or a body consisting of a continuous mass as in figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig2})

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig1}
\caption{Rigid body of point masses.}
\end{figure}

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig2}
\caption{Rigid solid body of continuous mass.}
\end{figure}

In the first figure our mass is made of discrete particles

\begin{aligned}M = \sum m_i\end{aligned} \hspace{\stretch{1}}(1.1)

whereas in the second figure with mass density \rho(\mathbf{r}) and a volume element d^3\mathbf{r}, our total mass is

\begin{aligned}M = \int_V \rho(\mathbf{r}) d^3\mathbf{r}\end{aligned} \hspace{\stretch{1}}(1.2)

Degrees of freedom

How many numbers do we need to describe fixed body motion. Consider figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig3})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig3}
\caption{Body local coordinate system with vector to a fixed point in the body}
\end{figure}

We will need to use six different numbers to describe the motion of a rigid body. We need three for the position of the body \mathbf{R}_{CM} as a whole. We also need three degrees of freedom (in general) for the motion of the body at that point in space (how our local coordinate system at the body move at that point), describing the change of the orientation of the body as a function of time.

Note that the angle \phi hasn’t been included in any of the pictures because it is too messy with all the rest. Picture something like figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig5})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig5}
\caption{Rotation angle and normal in the body}
\end{figure}

Let’s express the position of the body in terms of that body’s center of mass

\begin{aligned}\mathbf{R}_{CM} = \frac{\sum_i m_i \mathbf{r}_i}{\sum_j m_j},\end{aligned} \hspace{\stretch{1}}(1.3)

or for continuous masses

\begin{aligned}\mathbf{R}_{CM} = \frac{\int_V d^3 \mathbf{r}' \mathbf{r}' \rho(\mathbf{r}')}{ \int d^3 \mathbf{r}'' \rho(\mathbf{r}'')}.\end{aligned} \hspace{\stretch{1}}(1.4)

We consider the motion of point \mathbf{P}, an arbitrary point in the body as in figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig4}), whos motion consists of

\begin{enumerate}
\item displacement of the CM \mathbf{R}_{CM}
\item rotation of \mathbf{r} around some axis \hat{\mathbf{n}} going through CM on some angle \phi. (here \hat{\mathbf{n}} is a unit vector).
\end{enumerate}

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig4}
\caption{A point in the body relative to the center of mass.}
\end{figure}

From the picture we have

\begin{aligned}\boldsymbol{\rho} = \mathbf{R}_{CM} + \mathbf{r}\end{aligned} \hspace{\stretch{1}}(1.5)

\begin{aligned}d\boldsymbol{\rho} = d\mathbf{R}_{CM} + d\boldsymbol{\phi} \times \mathbf{r}\end{aligned} \hspace{\stretch{1}}(1.6)

where

\begin{aligned}d\boldsymbol{\phi} = \hat{\mathbf{n}} d\phi.\end{aligned} \hspace{\stretch{1}}(1.7)

Dividing by dt we have

\begin{aligned}\frac{d{{\boldsymbol{\rho}}}}{dt} = \frac{d{{\mathbf{R}_{CM}}}}{dt} + \frac{d{{\boldsymbol{\phi}}}}{dt} \times \mathbf{r}.\end{aligned} \hspace{\stretch{1}}(1.8)

The total velocity of this point in the body is then

\begin{aligned}\mathbf{v} = \mathbf{V}_{CM} = \boldsymbol{\Omega}_{CM} \times \mathbf{r}\end{aligned} \hspace{\stretch{1}}(1.9)

where

\begin{aligned}\boldsymbol{\Omega}_{CM} = \frac{d{{\boldsymbol{\phi}}}}{dt} = \frac{d{{(\hat{\mathbf{n}} \phi)}}}{dt} = \text{angular velocity of the body}.\end{aligned} \hspace{\stretch{1}}(1.10)

This circular motion is illustrated in figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig6})

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig6}
\caption{circular motion.}
\end{figure}

Note that \mathbf{v} is the velocity of the particle with respect to the unprimed system.

We’ll spend a lot of time figuring out how to express \boldsymbol{\Omega}_{CM}.

Now let’s consider a second point as in figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig7})
\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig7}
\caption{Two points in a rigid body.}
\end{figure}

\begin{aligned}\boldsymbol{\rho} &= \mathbf{R} + \mathbf{r} \\ \boldsymbol{\rho} &= \tilde{\mathbf{R}} + \tilde{\mathbf{r}} \\ \tilde{\mathbf{r}} &= \mathbf{r} + \mathbf{a}\end{aligned} \hspace{\stretch{1}}(1.11)

we have

\begin{aligned}\mathbf{v}_p&=\frac{d{{\boldsymbol{\rho}}}}{dt} \\ &= \frac{d{\mathbf{r}}}{dt} + \frac{d{{\boldsymbol{\phi}}}}{dt} \times \mathbf{r} \\ &= \frac{d{\mathbf{r}}}{dt} + \frac{d{{\boldsymbol{\phi}}}}{dt} \times (\tilde{\mathbf{r}} - \mathbf{a}) \\ &= \frac{d{\mathbf{r}}}{dt} - \frac{d{{\boldsymbol{\phi}}}}{dt} \times \mathbf{a} + \frac{d{{\boldsymbol{\phi}}}}{dt} \times \tilde{\mathbf{r}} \end{aligned}

\begin{aligned}\boxed{\mathbf{v}_p = \mathbf{V}_{CM} - \boldsymbol{\Omega}_{CM} \times \mathbf{a}+\boldsymbol{\Omega}_{CM} \times \tilde{\mathbf{r}}}\end{aligned} \hspace{\stretch{1}}(1.14)

Have another way that we can use to express the position of the point

\begin{aligned}\frac{d{{\boldsymbol{\rho}}}}{dt} = \frac{d{{\tilde{\mathbf{R}}}}}{dt} + \frac{d{{\tilde{\boldsymbol{\rho}}}}}{dt} \times \tilde{\mathbf{r}}_p \end{aligned} \hspace{\stretch{1}}(1.15)

or

\begin{aligned}\mathbf{v}_p = \mathbf{V}_A + \boldsymbol{\Omega}_A \times \tilde{\mathbf{r}}_p.\end{aligned} \hspace{\stretch{1}}(1.16)

Equating with above, and noting that this holds for all \tilde{\mathbf{r}}_p, and noting that if \tilde{\mathbf{r}}_p = 0

\begin{aligned}\mathbf{V}_A = \mathbf{V}_{CM} - \boldsymbol{\Omega}_{CM} \times \mathbf{a}\end{aligned} \hspace{\stretch{1}}(1.17)

hence

\begin{aligned}\boldsymbol{\Omega}_{CM} \times \tilde{\mathbf{r}}_p = \boldsymbol{\Omega}_A \times \tilde{\mathbf{r}}_p\end{aligned} \hspace{\stretch{1}}(1.18)

or

\begin{aligned}\boldsymbol{\Omega}_{CM} = \boldsymbol{\Omega}_A.\end{aligned} \hspace{\stretch{1}}(1.19)

The moral of the story is that the angular velocity \boldsymbol{\Omega} is a characteristic of the system. It does not matter if it is calculated with respect to the center of mass or not.

See some examples in the notes.

Kinetic energy

For all P in the body we have

\begin{aligned}\mathbf{v}_p = \mathbf{V}_A + \boldsymbol{\Omega} \times \mathbf{r}_p\end{aligned} \hspace{\stretch{1}}(2.20)

here \mathbf{V}_A is an arbitrary fixed point in the body as in figure (\ref{fig:rigidBodyMotion:rigidBodyMotionFig8})

\begin{figure}[htp]
\centering
\includegraphics[totalheight=0.2\textheight]{rigidBodyMotionFig8}
\caption{Kinetic energy setup relative to point A in the body.}
\end{figure}

The kinetic energy is

\begin{aligned}T &= \sum_a \frac{1}{{2}} m_a \boldsymbol{\rho}_a \\ &= \sum_a \frac{1}{{2}} \mathbf{v}_a^2 \\ &= \sum_a \frac{1}{{2}} \left(\mathbf{V}_A + \boldsymbol{\Omega} \times \mathbf{r}_a\right)^2 \\ &=\sum_a \frac{1}{{2}} \left( \mathbf{V}_A^2 +2 \mathbf{V}_A \cdot (\boldsymbol{\Omega} \times \mathbf{r}_a)+ (\boldsymbol{\Omega} \times \mathbf{r}_a)^2\right)\end{aligned}

We see that if we take A to be the center of mass then our cross term

\begin{aligned}\sum_a m_a \mathbf{V}_A \cdot (\boldsymbol{\Omega} \times \mathbf{r}_a)&=\mathbf{V}_A \cdot (\boldsymbol{\Omega} \times \sum_a m_a \mathbf{r}_a) \\ &=\mathbf{V}_A \cdot (\boldsymbol{\Omega} \times \mathbf{R}_{CM} )\end{aligned}

which vanishes. With

\begin{aligned}\mu = \sum_a m_a\end{aligned} \hspace{\stretch{1}}(2.21)

we have

\begin{aligned}T = \frac{1}{{2}} \mu \mathbf{V}_{CM}^2+ \frac{1}{{2}} \sum_a (\boldsymbol{\Omega} \times \mathbf{r}_a) \cdot (\boldsymbol{\Omega} \times \mathbf{r}_a) \end{aligned} \hspace{\stretch{1}}(2.22)

with

\begin{aligned}(\mathbf{A} \times \mathbf{B}) \cdot (\mathbf{C} \times \mathbf{D}) = (\mathbf{A} \cdot \mathbf{C}) (\mathbf{B} \cdot \mathbf{D}) -(\mathbf{A} \cdot \mathbf{D}) (\mathbf{B} \cdot \mathbf{C}),\end{aligned} \hspace{\stretch{1}}(2.23)

or

\begin{aligned}(\mathbf{A} \times \mathbf{B}) \cdot (\mathbf{A} \times \mathbf{B}) = \mathbf{A}^2 \mathbf{B}^2 - (\mathbf{A} \cdot \mathbf{B})^2.\end{aligned} \hspace{\stretch{1}}(2.24)

Forgetting about the \mu dependent term for now we have

\begin{aligned}T = \frac{1}{{2}} \sum_a m_a \left( \boldsymbol{\Omega}^2 \mathbf{r}_a^2 - (\boldsymbol{\Omega} \cdot \mathbf{r}_a)^2\right)\end{aligned} \hspace{\stretch{1}}(2.25)

Expanding this out with

\begin{aligned}\mathbf{r}_a = ( r_{a_1} r_{a_2}, r_{a_3} ) = \{ r_{a_i} \}\end{aligned} \hspace{\stretch{1}}(2.26)

and

\begin{aligned}\boldsymbol{\Omega} = ( \Omega_{1} \Omega_{2}, \Omega_{3} ) = \{ \Omega_{i} \}\end{aligned} \hspace{\stretch{1}}(2.27)

we have

\begin{aligned}T = \frac{1}{{2}} \sum_a m_a \left( \Omega_k \Omega_k r_{a_j} r_{a_j} - (\Omega_k r_{a_k})^2\right)\end{aligned} \hspace{\stretch{1}}(2.28)

Advertisements

Leave a Reply

Fill in your details below or click an icon to log in:

WordPress.com Logo

You are commenting using your WordPress.com account. Log Out / Change )

Twitter picture

You are commenting using your Twitter account. Log Out / Change )

Facebook photo

You are commenting using your Facebook account. Log Out / Change )

Google+ photo

You are commenting using your Google+ account. Log Out / Change )

Connecting to %s

 
%d bloggers like this: