## The Solution of Euler's Equations## Dynamic EquationsThe Euler equations are
The total angular momentum and initial conditions are
The Euler angles , , and using the transformation matrices give the equations
For display purposes, the initial Euler angles will be chosen so that output plots of the motion have the net angular momentum directed along the
Solve two equations for two initial angles, the angle being indeterminate.
The first solution in the list will be used, using a small initial with positive value, so that there is correspondence with the analytic solution found elsewhere for the case of the symmetric top when .
## Initial Conditions
The task is to solve the six parameters , , , and as a function of time, starting with initial values for each of the parameters, using the three principle moments of inertia for the Mir station. The principal axes of inertia of Mir roughly coincide with the module axes originating from the center section, known as the node. Although the Mir's moments of inertia cannot be known at any one time precisely because of the uncertain arrangement of internal stores throughout the Mir, the basic fixed hardware results in moments of inertia whose ratios can be characterized by setting . The initial values for are denoted to be respectively. The rate about the axis is denoted as
NDSolve. List the Mir physical properties as a single list, called dynamics.
Check the consistency of our parameters, by finding the minimum value of for
This should always be less than the value specified in
## Numerical Solutions
Using the parameters given, the Euler equations are solved simultaneously with the three equations for the Euler angles, applying our geometric constraint
To understand our solution from the point of view of the crew members on board Mir, plot the angular rates as observed about the three principle axes of Mir. These can be measured directly by the crew, utilizing the stars as an inertial reference. It is apparent that the motion is that of an irregular body, with most angular momentum about its middle axis of inertia, causing an unstable exchange of angular momentum between the principle axes. Although the Mir starts spinning initially about the axis, after the direction of rotation is reversed and oscillates with a period of roughly . The spin about the axis, initially small, grows periodically, but never reverses.
Solid curve is , short dashed curve is , long dashed curve is
No physical intuition is readily apparent when the Euler angles are plotted, but they determine the behavior of the physical transformation matrix used to do the animation of the Mir model. It should be noted that the roughly linear growth of with time corresponds to the continuous spin of the Mir complex about the inertial Z axis, along which the total angular momentum was constrained as an initial condition using anglerule
Solid curve is , short dashed curve is , long dashed curve is Converted by Mathematica
October 6, 1999
[Prev Page][Next Page] |