The equations of motion
I have started the Dzhanibekov effect series with quaternions, but that is because of my personal attraction to quaternions. Some years ago I coauthored a paper on Quaternions and Monopoles. The monopoles there are the exotic “magnetic monopoles” that are lake fairies. Some physicist say they have seen them in the woods, but the majority of physicists take it as a joke. Quaternions are also somewhat exotic. I do not know if magnetic monopoles have any relation to Dzhanibekov effect, though, in fact they may have one – but that is for the physics of the future. For the physics of today, for the physics of a cosmic spinning and flipping nut, we need the physics of a rigid body. That is part of classical mechanics.
According to Wikipedia
In physics, a rigid body is an idealization of a solid body in which deformation is neglected. In other words, the distance between any two given points of a rigid body remains constant in time regardless of external forces exerted on it.
I should say that the concept of a rigid body becomes somewhat iffy when we want to take into account what Albert Einstein taught us about first special, and then general relativity. I did not see yet a satisfactory relativistic theory of Dzhanibekov effect. What I know is that some physicists trying to sell their fantastic ideas to the investors (example G. Shipov) try to relate Dzhanibekov effect to “torsion fields” and “4D gyroscopes“. I have published some comments on this subject in International Journal of Unconventional Science. “Comments on Chapter 5 of G. I. Shipov’s “A Theory of Physical Vacuum”. Part I” is available in English. The second part, dealing with Shipov’s “4D gyroscopes” is for now available only in Russian, but it will be translated into English soon.
But that is all in the exotic sphere. Interesting, for sure, but my aim is to describe the physics and mathematics of Dzhanibekov effect using old and good classical mechanics of Newton and Euler.
We consider a free rigid body, observed form an inertial reference system, rotating with a fixed point at the center of its mass. Let be an orthonormal frame corotating with the body, and aligned with its principal axes, and let be an inertial laboratory frame, both centered at the center of mass of the body. The two frames are related by time-dependent orthogonal matrix
Then
is often called the attitude matrix. For a rotating body, if are coordinates of a fixed point in the body, then its coordinates in the laboratory system change in time:
Differentiating we get
The matrix is orthogonal, therefore, by differentiating, the matric is antisymmetric. Every antisymmetric matrix can be written as
Writing
(*)
we find that
The vector is the angular velocity vector in the laboratory frame.
We denote by the body representative of the angular velocity vector where we skip the time dependence, with components
Thus
and, from Eq. (*)
where In the body frame the inertia tensor is diagonal The angular momentum vector in the body frame is then given by the formula
and the Euler equations for the free rotation read:
When written componentwise in terms of the angular velocity components they take the well known form
(1)
(2)
Now we need to solve explicitly Euler’s equations, which is just the first part of the task. Then we have to find the attitude matrix. Then we have to use these explicit solutions in order to produce realistic animations of Dzhanibekov’s effect. It is a lot of work, especially if one is going to reveal all the details. So, let us see where it will lead us to.
Just testing if I can post an image here
My cosmic 4D gyroscopes!
Test of Registeration. My question was Qusetion:- what is the difference between “inertial reference system” and “inertial laboratory frame”? I guess one tumbles with the object.
I now understand after re-reading. The “inertial reference system” is rotating with the object and the “inertial laboratory frame” is fixed.
I am using here “inertial reference frame” and “inertial laboratory frame” to mean the same – a frame attached to a body that is not rotating nad not accelerating with respect to the fixed stars, or, if you wish, to “Aether”. Well, approximately.
From Wikipedia:
…
Newton viewed the first law as valid in any reference frame that is in uniform motion relative to the fixed stars;[32] that is, neither rotating nor accelerating relative to the stars.[33] Today the notion of “absolute space” is abandoned, and an inertial frame in the field of classical mechanics is defined as:[34][35]
An inertial frame of reference is one in which the motion of a particle not subject to forces is in a straight line at constant speed.
….
On the other hand there is another frame, tumbling with the object. That frame is not an inertial frame (there is a centrifugal force acting within this frame). I will often call it the rotating frame. I am assuming that this frame is oriented so that its axes point in the directions of “principal axes” see Moment of inertia in Wikipedia
If you would like to have further clarification – do not hesitate to ask.
Thank you for the above. Oddly I had been meaning to study up on Inertia matrices for a good while now. My nest question is on the form of inputs to Q or A. I am treating it as a rotoation matrix made up from RotZ . RotY .RotX where the individual matrices are of the form RotZ=[[cos(omega1 t),-sin(omega1 t),0],[sin(omega1 t),cos(omega1 t),0],[0,0,1]] ,omega2 for Y and omega3 for X. Difficult to copy and paste anything nice looking. Following the process it produced the antisymmetric matrix fine. Have produced the 3 Euler eqns. Can you explain where A1 to m come from.
I find it takes a lot of thinking about to understand each step