Euler's equations (rigid body dynamics)

This article is about rigid body dynamics. For other uses, see Euler equations (disambiguation).

Classical mechanics
${\vec {F}}=m{\vec {a}}$ Second law of motion
History Timeline
Branches Applied Celestial Continuum Dynamics Kinematics Kinetics Statics Statistical
Fundamentals Acceleration Angular momentum Couple D'Alembert's principle Energy kinetic potential Force Frame of reference Impulse Inertia / Moment of inertia Mass Mechanical power Mechanical work Moment Momentum Space Speed Time Torque Velocity Virtual work
Formulations Newton's laws of motion Analytical mechanics Lagrangian mechanics Hamiltonian mechanics Routhian mechanics Hamilton–Jacobi equation Appell's equation of motion Udwadia–Kalaba equation Koopman–von Neumann mechanics
Core topics Damping (ratio) Displacement Equations of motion Euler's laws of motion Fictitious force Friction Harmonic oscillator Inertial / Non-inertial reference frame Mechanics of planar particle motion Motion (linear) Newton's law of universal gravitation Newton's laws of motion Relative velocity Rigid body dynamics Euler's equations Simple harmonic motion Vibration
Rotation Circular motion Rotating reference frame Centripetal force Centrifugal force reactive Coriolis force Pendulum Tangential speed Rotational speed Angular acceleration / displacement / frequency / velocity
Scientists Galileo Newton Kepler Horrocks Halley Euler d'Alembert Clairaut Lagrange Laplace Hamilton Poisson Daniel Bernoulli Johann Bernoulli Cauchy

In classical mechanics, Euler's rotation equations are a vectorial quasilinear first-order ordinary differential equation describing the rotation of a rigid body, using a rotating reference frame with its axes fixed to the body and parallel to the body's principal axes of inertia. Their general form is:

{\mathbf {I}}\cdot {\dot {{\boldsymbol \omega }}}+{\boldsymbol \omega }\times \left({\mathbf {I}}\cdot {\boldsymbol \omega }\right)={\mathbf {M}}.

where M is the applied torques, I is the inertia matrix, and ω is the angular velocity about the principal axes.

In 3D principal orthogonal coordinates, they become:

{\begin{aligned}I_{1}{\dot {\omega }}_{{1}}+(I_{3}-I_{2})\omega _{2}\omega _{3}&=M_{{1}}\\I_{2}{\dot {\omega }}_{{2}}+(I_{1}-I_{3})\omega _{3}\omega _{1}&=M_{{2}}\\I_{3}{\dot {\omega }}_{{3}}+(I_{2}-I_{1})\omega _{1}\omega _{2}&=M_{{3}}\end{aligned}}

where M_k are the components of the applied torques, I_k are the principal moments of inertia and ω_k are the components of the angular velocity about the principal axes.

Motivation and derivation

Starting from Euler's second law, in an inertial frame of reference (subscripted "in"), the time derivative of the angular momentum L equals the applied torque

{\frac {d{\mathbf {L}}_{{{\text{in}}}}}{dt}}\ {\stackrel {{\mathrm {def}}}{=}}\ {\frac {d}{dt}}\left({\mathbf {I}}_{{{\text{in}}}}\cdot {\boldsymbol \omega }\right)={\mathbf {M}}_{{{\text{in}}}}

where I_in is the moment of inertia tensor calculated in the inertial frame. Although this law is universally true, it is not always helpful in solving for the motion of a general rotating rigid body, since both I_in and ω can change during the motion.

Therefore, we change to a coordinate frame fixed in the rotating body, and chosen so that its axes are aligned with the principal axes of the moment of inertia tensor. In this frame, at least the moment of inertia tensor is constant (and diagonal), which simplifies calculations. As described in the moment of inertia, the angular momentum L can be written

{\mathbf {L}}\ {\stackrel {{\mathrm {def}}}{=}}\ L_{{1}}{\mathbf {e}}_{{1}}+L_{{2}}{\mathbf {e}}_{{2}}+L_{{3}}{\mathbf {e}}_{{3}}=I_{{1}}\omega _{{1}}{\mathbf {e}}_{{1}}+I_{{2}}\omega _{{2}}{\mathbf {e}}_{{2}}+I_{{3}}\omega _{{3}}{\mathbf {e}}_{{3}}

where M_k, I_k and ω_k are as above.

In a rotating reference frame, the time derivative must be replaced with (see time derivative in rotating reference frame)

\left({\frac {d{\mathbf {L}}}{dt}}\right)_{{\mathrm {rot}}}+{\boldsymbol \omega }\times {\mathbf {L}}={\mathbf {M}}

where the subscript "rot" indicates that it is taken in the rotating reference frame. The expressions for the torque in the rotating and inertial frames are related by

{\mathbf {M}}_{{{\text{in}}}}={\mathbf {Q}}{\mathbf {M}},

where Q is the rotation tensor(not rotation matrix), an orthogonal tensor related to the angular velocity vector by

{\boldsymbol \omega }\times {\boldsymbol {v}}={\dot {{\mathbf {Q}}}}{\mathbf {Q}}^{{-1}}{\boldsymbol {v}}

for any vector v.

In general, L = I·ω is substituted and the time derivatives are taken realizing that the inertia tensor, and so also the principal moments, do not depend on time. This leads to the general vector form of Euler's equations

{\mathbf {I}}\cdot {\dot {{\boldsymbol \omega }}}+{\boldsymbol \omega }\times \left({\mathbf {I}}\cdot {\boldsymbol \omega }\right)={\mathbf {M}}.

If principal axis rotation

L_{{k}}\ {\stackrel {{\mathrm {def}}}{=}}\ I_{{k}}\omega _{{k}}

is substituted, and then taking the cross product and using the fact that the principal moments do not change with time, we arrive at the Euler equations in components at the beginning of the article.

Torque-free solutions

For the RHSs equal to zero there are non-trivial solutions: torque-free precession. Notice that if I is constant (because the inertia tensor is the 3×3 identity matrix, because we work in the intrinsic frame, or because the torque is driving the rotation around the same axis $\mathbf {\hat {n}}$ so that I is not changing) then we may write

{\mathbf {M}}\ {\stackrel {{\mathrm {def}}}{=}}\ I{\frac {d\omega }{dt}}{\mathbf {{\hat {n}}}}=I\alpha {\mathbf {{\hat {n}}}}

where

α is called the angular acceleration (or rotational acceleration) about the rotation axis

\mathbf {\hat {n}}

However, if I is not constant in the external reference frame (i.e. the body is moving and its inertia tensor is not the identity) then we cannot take the I outside the derivative. In this case we will have torque-free precession, in such a way that I(t) and ω(t) change together so that their derivative is zero. This motion can be visualized by Poinsot's construction.

Generalizations

It is also possible to use these equations if the axes in which

\left({\frac {d{\mathbf {L}}}{dt}}\right)_{{\mathrm {relative}}}

is described are not connected to the body. Then ω should be replaced with the rotation of the axes instead of the rotation of the body. It is, however, still required that the chosen axes are still principal axes of inertia. This form of the Euler equations is useful for rotation-symmetric objects that allow some of the principal axes of rotation to be chosen freely.

References

C. A. Truesdell, III (1991) A First Course in Rational Continuum Mechanics. Vol. 1: General Concepts, 2nd ed., Academic Press. ISBN 0-12-701300-8. Sects. I.8-10.
C. A. Truesdell, III and R. A. Toupin (1960) The Classical Field Theories, in S. Flügge (ed.) Encyclopedia of Physics. Vol. III/1: Principles of Classical Mechanics and Field Theory, Springer-Verlag. Sects. 166-168, 196-197, and 294.
Landau L.D. and Lifshitz E.M. (1976) Mechanics, 3rd. ed., Pergamon Press. ISBN 0-08-021022-8 (hardcover) and ISBN 0-08-029141-4 (softcover).
Goldstein H. (1980) Classical Mechanics, 2nd ed., Addison-Wesley. ISBN 0-201-02918-9
Symon KR. (1971) Mechanics, 3rd. ed., Addison-Wesley. ISBN 0-201-07392-7

This article is issued from Wikipedia - version of the 11/23/2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.

Euler's equations (rigid body dynamics)

Motivation and derivation

Torque-free solutions

Generalizations

See also

References