Take the 2-minute tour ×
Physics Stack Exchange is a question and answer site for active researchers, academics and students of physics. It's 100% free, no registration required.

I am in search of a simplified version of the derivation of Newton-Euler equations of motion (both translational and rotational) for a rigid body (3D block) that has a body fixed frame and where the center of mass of the body is not at the center of gravity. I can find elementary derivations for the same system when the center of mass is at the center of gravity, but not for my system in question.

I am using the derivation as background research for a rotordynamics project that I'm working on.

Any help and or references would be greatly appreciated! :)

share|improve this question
add comment

1 Answer

up vote 1 down vote accepted

It depends if you already have the mass moment of inertia tensor defined or not.

If you know the body inertia is $I_{body}$ and the 3x3 rotation matrix is $E$ then the angular momentum vector at the center of gravity C is

$$ \vec{H}_C = \left( E I_{body} E^\top \right) \vec{\omega} $$

and the linear momentum vector is $$\vec{L} = m \vec{v}_C$$

The mass moment of inertia tensor along the world coordinates on the center of gravity is $I_C = E I_{body} E^\top $ which transforms the rotational velocity $\vec \omega$ into local coordinates, multiplies by $I_{body}$ and transforms backs into world coordinates.

Now the equations of motion on the center of gravity are defined from the sum of forces and moments equals the rate of change of momentum

$$ \sum \vec{F} = \dot{\vec{L}} $$ $$ \sum \vec{M}_C = \dot{\vec{H}}_C $$

or

$$ \sum \vec{F} = m \vec{a}_C $$ $$ \sum \vec{M}_C = I_C \vec{\alpha} + \vec{\omega} \times I_C \vec{\omega} $$

since the time derivative of angular momentum on a rotating frame is $\dot{\vec{H}_C} = \frac{\partial \vec{H}_C}{\partial t} + \vec{\omega} \times \vec{H}_C $

Note that $\dot{\vec{v}} = \vec{a} $ and $\dot{\vec{\omega}} = \vec{\alpha} $.

Now to describe the equations on a frame A not on C use the following transformations (with relative position of the c.g. $ \vec{c} =\vec{r}_C - \vec{r}_A $.

$$ \vec{a}_C = \vec{a}_A + \vec{\alpha} \times \vec{c} + \vec{\omega} \times \vec{\omega} \times \vec{c} $$

$$ \sum \vec{M}_C = \sum \vec{M}_A + \vec{c} \times \sum \vec{F} $$

So finally the equations of motion of a rigid body, as described by a frame A not on the center of gravity C is (rather messy)

$$ \boxed{ \begin{aligned} \sum \vec{F} &= m \vec{a}_A - m \vec{c}\times \vec{\alpha} + m \vec{\omega}\times\vec{\omega}\times\vec{c} \\ \sum \vec{M}_A &= I_C \vec{\alpha} + m \vec{c} \times \vec{a}_A - m \vec{c} \times \vec{c} \times \vec{\alpha} +\vec{\omega} I_C \vec{\omega} + m \vec{c} \times \left( \vec{\omega} \times \vec{\omega} \times \vec{c} \right) \end{aligned} } $$

Which is why people use the spatial notation (look up screw theory) to compact the above into

$$ \sum \bf{f}_A = I_A \bf{a}_A + \bf{p} $$

$$ \begin{pmatrix} \sum \vec{F} \\ \sum \vec{M}_A \end{pmatrix} = \begin{bmatrix} m & -m \vec{c}\times \\ m \vec{c}\times & I_C - m \vec{c}\times \vec{c}\times \end{bmatrix} \begin{pmatrix} \vec{a}_A \\ \vec{\alpha} \end{pmatrix} + \begin{bmatrix} 1 & 0 \\ \vec{c}\times & 1 \end{bmatrix} \begin{pmatrix} m \vec{\omega}\times\vec{\omega}\times\vec{c} \\ \vec{\omega}\times I_C \vec{\omega} \end{pmatrix} $$

Notice the above the $0$ and $1$ are 3x3 matrices and $\vec{c}\times$ is the 3x3 cross product operator defined by

$$ \vec{ \begin{pmatrix} x\\y\\z \end{pmatrix} }\times = \begin{vmatrix} 0 & -z & y \\ z & 0 & -x \\ -y & x & 0 \end{vmatrix} $$

Now the big 6x6 matrix multiplying the acceleration term is the spatial inertia at A. More here and here.

share|improve this answer
    
~ that is exactly what I was looking for. I found the EOMs in Wikipedia, but no derivation. Thank you sooo much!! –  Lanae Oct 12 '13 at 1:59
    
Phew - my graduate degree was not in vain. –  ja72 Oct 12 '13 at 2:00
    
~I am working on an applied project for my Masters in Aerospace Engineering. The simplest things from my math days just elude me sometimes. The actual equations that I have to work with are daunting and I wanted to be able to see the basics first. :) –  Lanae Oct 12 '13 at 2:03
    
@Lanae - Yes the Wiki pages are "screwed" up because multiple people from various backgrounds have edited them (like molecular dynamics) and each uses different notations and conventions. You are welcome to use my thesis presentation to get some ideas. –  ja72 Oct 12 '13 at 2:07
    
~ you mention a 6X6 matrix multiplying the acceleration term being the spatial inertia at A. ?? My end game is to simulate the new I_g, randomize it, solve the EOMs and get the response (natural frequencies). Therefore, I need to find the new a_g and how it relates to the system since C does not equal A. –  Lanae Oct 12 '13 at 2:16
show 1 more comment

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.