Tennis racket theorem

Last updated November 29, 2024

The tennis racket theorem or intermediate axis theorem, is a kinetic phenomenon of classical mechanics which describes the movement of a rigid body with three distinct principal moments of inertia. It has also been dubbed the Dzhanibekov effect, after Soviet cosmonaut Vladimir Dzhanibekov, who noticed one of the theorem's logical consequences whilst in space in 1985.^[1] The effect was known for at least 150 years prior, having been described by Louis Poinsot in 1834^[2]^[3] and included in standard physics textbooks such as Classical Mechanics by Herbert Goldstein throughout the 20th century.

This can be demonstrated by the following experiment: Hold a tennis racket at its handle, with its face being horizontal, and throw it in the air such that it performs a full rotation around its horizontal axis perpendicular to the handle (ê₂ in the diagram), and then catch the handle. In almost all cases, during that rotation the face will also have completed a half rotation, so that the other face is now up. By contrast, it is easy to throw the racket so that it will rotate around the handle axis (ê₁) without accompanying half-rotation around another axis; it is also possible to make it rotate around the vertical axis perpendicular to the handle (ê₃) without any accompanying half-rotation.

The experiment can be performed with any object that has three different moments of inertia, for instance with a (rectangular) book, remote control, or smartphone. The effect occurs whenever the axis of rotation differs – even only slightly – from the object's second principal axis; air resistance or gravity are not necessary.^[4]

Theory

Dzhanibekov effect demonstration in microgravity, NASA.

The tennis racket theorem can be qualitatively analysed with the help of Euler's equations. Under torque–free conditions, they take the following form:

${\begin{aligned}I_{1}{\dot {\omega }}_{1}&=-(I_{3}-I_{2})\omega _{3}\omega _{2}~~~~~~~~~~~~~~~~~~~~{\text{(1)}}\\I_{2}{\dot {\omega }}_{2}&=-(I_{1}-I_{3})\omega _{1}\omega _{3}~~~~~~~~~~~~~~~~~~~~{\text{(2)}}\\I_{3}{\dot {\omega }}_{3}&=-(I_{2}-I_{1})\omega _{2}\omega _{1}~~~~~~~~~~~~~~~~~~~~{\text{(3)}}\end{aligned}}$

Here $I_{1},I_{2},I_{3}$ denote the object's principal moments of inertia, and we assume $I_{1}<I_{2}<I_{3}$ . The angular velocities around the object's three principal axes are $\omega _{1},\omega _{2},\omega _{3}$ and their time derivatives are denoted by ${\dot {\omega }}_{1},{\dot {\omega }}_{2},{\dot {\omega }}_{3}$ .

Stable rotation around the first and third principal axis

Consider the situation when the object is rotating around the axis with moment of inertia $I_{1}$ . To determine the nature of equilibrium, assume small initial angular velocities along the other two axes. As a result, according to equation (1), $~{\dot {\omega }}_{1}$ is very small. Therefore, the time dependence of $~\omega _{1}$ may be neglected.

Now, differentiating equation (2) and substituting ${\dot {\omega }}_{3}$ from equation (3),

${\begin{aligned}I_{2}{\ddot {\omega }}_{2}&=-(I_{1}-I_{3})\omega _{1}{\dot {\omega }}_{3}\\I_{3}I_{2}{\ddot {\omega }}_{2}&=(I_{1}-I_{3})(I_{2}-I_{1})(\omega _{1})^{2}\omega _{2}\\{\text{i.e. }}~~~~{\ddot {\omega }}_{2}&={\text{(negative quantity)}}\cdot \omega _{2}\end{aligned}}$

because $I_{2}-I_{1}>0$ and $I_{1}-I_{3}<0$ .

Note that $\omega _{2}$ is being opposed and so rotation around this axis is stable for the object.

Similar reasoning gives that rotation around the axis with moment of inertia $I_{3}$ is also stable.

Unstable rotation around the second principal axis

Now apply the same analysis to the axis with moment of inertia $I_{2}.$ This time ${\dot {\omega }}_{2}$ is very small. Therefore, the time dependence of $~\omega _{2}$ may be neglected.

Now, differentiating equation (1) and substituting ${\dot {\omega }}_{3}$ from equation (3),

${\begin{aligned}I_{1}I_{3}{\ddot {\omega }}_{1}&=(I_{3}-I_{2})(I_{2}-I_{1})(\omega _{2})^{2}\omega _{1}\\{\text{i.e.}}~~~~{\ddot {\omega }}_{1}&={\text{(positive quantity)}}\cdot \omega _{1}\end{aligned}}$

Note that $\omega _{1}$ is not opposed (and therefore will grow) and so rotation around the second axis is unstable. Therefore, even a small disturbance, in the form of a very small initial value of $\omega _{1}$ or $\omega _{3}$ , causes the object to 'flip'.

Matrix analysis

If the object is mostly rotating along its third axis, so $|\omega _{3}|\gg |\omega _{1}|,|\omega _{2}|$ , we can assume $\omega _{3}$ does not vary much, and write the equations of motion as a matrix equation: ${\frac {d}{dt}}{\begin{bmatrix}\omega _{1}\\\omega _{2}\end{bmatrix}}={\begin{bmatrix}0&-\omega _{3}(I_{3}-I_{2})/I_{1}\\-\omega _{3}(I_{1}-I_{3})/I_{2}&0\end{bmatrix}}{\begin{bmatrix}\omega _{1}\\\omega _{2}\end{bmatrix}}$ which has zero trace and positive determinant, implying the motion of $(\omega _{1},\omega _{2})$ is a stable rotation around the origin—a neutral equilibrium point. Similarly, the point $(\omega _{1},0,0)$ is a neutral equilibrium point, but $(0,\omega _{2},0)$ is a saddle point.

Geometric analysis

During motion, both the energy and angular momentum-squared are conserved, thus we have two conserved quantities: ${\begin{cases}2E=\sum _{i}I_{i}\omega _{i}^{2}\\L^{2}=\sum _{i}I_{i}^{2}\omega _{i}^{2}\end{cases}}$ and so for any initial condition $\omega (0)$ , the trajectory of $\omega (t)$ must stay on the intersection curve between two ellipsoids defined by ${\begin{cases}\sum _{i}I_{i}\omega _{i}^{2}=\sum _{i}I_{i}\omega _{i}(0)^{2}\\\sum _{i}I_{i}^{2}\omega _{i}^{2}=\sum _{i}I_{i}^{2}\omega _{i}(0)^{2}\end{cases}}$ This is shown on the animation to the left.

By inspecting Euler's equations, we see that ${\dot {\omega }}(t)=0$ implies that two components of $\omega (t)$ are zero—that is, the object is exactly spinning around one of the principal axes. In all other situations, $\omega (t)$ must remain in motion.

By Euler's equations, if $\omega (t)$ is a solution, then so is $c\omega (ct)$ for any constant $c>0$ . In particular, the motion of the body in free space (obtained by integrating $c\omega (ct)dt$ ) is exactly the same, just completed faster by a ratio of $c$ .

Consequently, we can analyze the geometry of motion with a fixed value of $L^{2}$ , and vary $\omega (0)$ on the fixed ellipsoid of constant squared angular momentum. As $\omega (0)$ varies, the value of $2E$ also varies—thus giving us a varying ellipsoid of constant energy. This is shown in the animation as a fixed orange ellipsoid and increasing blue ellipsoid.

For concreteness, consider $I_{1}=1,I_{2}=2,I_{3}=3$ , then the angular momentum ellipsoid's major axes are in ratios of $1:1/2:1/3$ , and the energy ellipsoid's major axes are in ratios of $1:1/{\sqrt {2}}:1/{\sqrt {3}}$ . Thus the angular momentum ellipsoid is both flatter and sharper, as visible in the animation. In general, the angular momentum ellipsoid is always more "exaggerated" than the energy ellipsoid.

Now inscribe on a fixed ellipsoid of $L^{2}$ its intersection curves with the ellipsoid of $2E$ , as $2E$ increases from zero to infinity. We can see that the curves evolve as follows:

For small energy, there is no intersection, since we need a minimum of energy to stay on the angular momentum ellipsoid.
The energy ellipsoid first intersects the momentum ellipsoid when $2E=L^{2}/I_{3}$ , at the points $(0,0,\pm L/I_{3})$ . This is when the body rotates around its axis with the largest moment of inertia.
They intersect at two cycles around the points $(0,0,\pm L/I_{3})$ . Since each cycle contains no point at which ${\dot {\omega }}=0$ , the motion of $\omega (t)$ must be a periodic motion around each cycle.
They intersect at two "diagonal" curves that intersects at the points $(0,\pm L/I_{2},0)$ , when $2E=L^{2}/I_{2}$ . If $\omega (t)$ starts anywhere on the diagonal curves, it would approach one of the points, distance exponentially decreasing, but never actually reach the point. In other words, we have 4 heteroclinic orbits between the two saddle points.
They intersect at two cycles around the points $(\pm L/I_{1},0,0)$ . Since each cycle contains no point at which ${\dot {\omega }}=0$ , the motion of $\omega (t)$ must be a periodic motion around each cycle.
The energy ellipsoid last intersects the momentum ellipsoid when $2E=L^{2}/I_{1}$ , at the points $(\pm L/I_{1},0,0)$ . This is when the body rotates around its axis with the smallest moment of inertia.

The tennis racket effect occurs when $\omega (0)$ is very close to a saddle point. The body would linger near the saddle point, then rapidly move to the other saddle point, near $\omega (T/2)$ , linger again for a long time, and so on. The motion repeats with period $T$ .

The above analysis is all done in the perspective of an observer which is rotating with the body. An observer watching the body's motion in free space would see its angular momentum vector ${\vec {L}}=I{\vec {\omega }}$ conserved, while both its angular velocity vector ${\vec {\omega }}(t)$ and its moment of inertia $I(t)$ undergo complicated motions in space. At the beginning, the observer would see both ${\vec {\omega }}(0),{\vec {L}}$ mostly aligned with the second major axis of $I(0)$ . After a while, the body performs a complicated motion and ends up with $I(T/2),{\vec {\omega }}(T/2)$ , and again both ${\vec {L}},{\vec {\omega }}(T/2)$ are mostly aligned with the second major axis of $I(T/2)$ .

Consequently, there are two possibilities: either the rigid body's second major axis is in the same direction, or it has reversed direction. If it is still in the same direction, then ${\vec {\omega }}(0),{\vec {\omega }}(T/2)$ viewed in the rigid body's reference frame are also mostly in the same direction. However, we have just seen that $\omega (0)$ and $\omega (T/2)$ are near opposite saddle points $(0,\pm L/I_{2},0)$ . Contradiction.

Qualitatively, then, this is what an observer watching in free space would observe:

The body rotates around its second major axis for a while.
The body rapidly undergoes a complicated motion, until its second major axis has reversed direction.
The body rotates around its second major axis again for a while. Repeat.

This can be easily seen in the video demonstration in microgravity.

With dissipation

When the body is not exactly rigid, but can flex and bend or contain liquid that sloshes around, it can dissipate energy through its internal degrees of freedom. In this case, the body still has constant angular momentum, but its energy would decrease, until it reaches the minimal point. As analyzed geometrically above, this happens when the body's angular velocity is exactly aligned with its axis of maximal moment of inertia.

This happened to Explorer 1, the first satellite launched by the United States in 1958. The elongated body of the spacecraft had been designed to spin about its long (least-inertia) axis but refused to do so, and instead started precessing due to energy dissipation from flexible structural elements.

In general, celestial bodies large or small would converge to a constant rotation around its axis of maximal moment of inertia. Whenever a celestial body is found in a complex rotational state, it is either due to a recent impact or tidal interaction, or is a fragment of a recently disrupted progenitor.^[5]

Related Research Articles

<span class="mw-page-title-main">Angular momentum</span> Conserved physical quantity; rotational analogue of linear momentum

Angular momentum is the rotational analog of linear momentum. It is an important physical quantity because it is a conserved quantity – the total angular momentum of a closed system remains constant. Angular momentum has both a direction and a magnitude, and both are conserved. Bicycles and motorcycles, flying discs, rifled bullets, and gyroscopes owe their useful properties to conservation of angular momentum. Conservation of angular momentum is also why hurricanes form spirals and neutron stars have high rotational rates. In general, conservation limits the possible motion of a system, but it does not uniquely determine it.

<span class="mw-page-title-main">Precession</span> Periodic change in the direction of a rotation axis

Precession is a change in the orientation of the rotational axis of a rotating body. In an appropriate reference frame it can be defined as a change in the first Euler angle, whereas the third Euler angle defines the rotation itself. In other words, if the axis of rotation of a body is itself rotating about a second axis, that body is said to be precessing about the second axis. A motion in which the second Euler angle changes is called nutation. In physics, there are two types of precession: torque-free and torque-induced.

A gyrocompass is a type of non-magnetic compass which is based on a fast-spinning disc and the rotation of the Earth to find geographical direction automatically. A gyrocompass makes use of one of the seven fundamental ways to determine the heading of a vehicle. A gyroscope is an essential component of a gyrocompass, but they are different devices; a gyrocompass is built to use the effect of gyroscopic precession, which is a distinctive aspect of the general gyroscopic effect. Gyrocompasses, such as the fibre optic gyrocompass are widely used to provide a heading for navigation on ships. This is because they have two significant advantages over magnetic compasses:

Kinematics is a subfield of physics and mathematics, developed in classical mechanics, that describes the motion of points, bodies (objects), and systems of bodies without considering the forces that cause them to move. Kinematics, as a field of study, is often referred to as the "geometry of motion" and is occasionally seen as a branch of both applied and pure mathematics since it can be studied without considering the mass of a body or the forces acting upon it. A kinematics problem begins by describing the geometry of the system and declaring the initial conditions of any known values of position, velocity and/or acceleration of points within the system. Then, using arguments from geometry, the position, velocity and acceleration of any unknown parts of the system can be determined. The study of how forces act on bodies falls within kinetics, not kinematics. For further details, see analytical dynamics.

<span class="mw-page-title-main">Angular velocity</span> Direction and rate of rotation

In physics, angular velocity, also known as angular frequency vector, is a pseudovector representation of how the angular position or orientation of an object changes with time, i.e. how quickly an object rotates around an axis of rotation and how fast the axis itself changes direction.

In fluid dynamics, a vortex is a region in a fluid in which the flow revolves around an axis line, which may be straight or curved. Vortices form in stirred fluids, and may be observed in smoke rings, whirlpools in the wake of a boat, and the winds surrounding a tropical cyclone, tornado or dust devil.

<span class="mw-page-title-main">Moment of inertia</span> Scalar measure of the rotational inertia with respect to a fixed axis of rotation

The moment of inertia, otherwise known as the mass moment of inertia, angular/rotational mass, second moment of mass, or most accurately, rotational inertia, of a rigid body is defined relative to a rotational axis. It is the ratio between the torque applied and the resulting angular acceleration about that axis. It plays the same role in rotational motion as mass does in linear motion. A body's moment of inertia about a particular axis depends both on the mass and its distribution relative to the axis, increasing with mass & distance from the axis.

Rotational energy or angular kinetic energy is kinetic energy due to the rotation of an object and is part of its total kinetic energy. Looking at rotational energy separately around an object's axis of rotation, the following dependence on the object's moment of inertia is observed: $where$

In the physical science of dynamics, rigid-body dynamics studies the movement of systems of interconnected bodies under the action of external forces. The assumption that the bodies are rigid simplifies analysis, by reducing the parameters that describe the configuration of the system to the translation and rotation of reference frames attached to each body. This excludes bodies that display fluid, highly elastic, and plastic behavior.

Dynamical simulation, in computational physics, is the simulation of systems of objects that are free to move, usually in three dimensions according to Newton's laws of dynamics, or approximations thereof. Dynamical simulation is used in computer animation to assist animators to produce realistic motion, in industrial design, and in video games. Body movement is calculated using time integration methods.

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 angular velocity ω whose axes are fixed to the body. They are named in honour of Leonhard Euler. Their general vector form is

<span class="mw-page-title-main">Differential wheeled robot</span> Robot with a particular driving system

A differential wheeled robot is a mobile robot whose movement is based on two separately driven wheels placed on either side of the robot body. It can thus change its direction by varying the relative rate of rotation of its wheels and hence does not require an additional steering motion. Robots with such a drive typically have one or more caster wheels to prevent the vehicle from tilting.

In rotordynamics, the rigid rotor is a mechanical model of rotating systems. An arbitrary rigid rotor is a 3-dimensional rigid object, such as a top. To orient such an object in space requires three angles, known as Euler angles. A special rigid rotor is the linear rotor requiring only two angles to describe, for example of a diatomic molecule. More general molecules are 3-dimensional, such as water, ammonia, or methane.

Screw theory is the algebraic calculation of pairs of vectors, also known as dual vectors – such as angular and linear velocity, or forces and moments – that arise in the kinematics and dynamics of rigid bodies.

In theoretical physics a Coriolis field is one of the apparent gravitational fields felt by a rotating or forcibly-accelerated body, together with the centrifugal field and the Euler field.

<span class="mw-page-title-main">Rotation around a fixed axis</span> Type of motion

Rotation around a fixed axis or axial rotation is a special case of rotational motion around an axis of rotation fixed, stationary, or static in three-dimensional space. This type of motion excludes the possibility of the instantaneous axis of rotation changing its orientation and cannot describe such phenomena as wobbling or precession. According to Euler's rotation theorem, simultaneous rotation along a number of stationary axes at the same time is impossible; if two rotations are forced at the same time, a new axis of rotation will result.

In classical mechanics, Poinsot's construction is a geometrical method for visualizing the torque-free motion of a rotating rigid body, that is, the motion of a rigid body on which no external forces are acting. This motion has four constants: the kinetic energy of the body and the three components of the angular momentum, expressed with respect to an inertial laboratory frame. The angular velocity vector $of the rigid rotor is not constant, but satisfies Euler's equations. The conservation of kinetic energy and angular momentum provide two constraints on the motion of .$

In classical mechanics, the Hannay angle is a mechanics analogue of the geometric phase. It was named after John Hannay of the University of Bristol, UK. Hannay first described the angle in 1985, extending the ideas of the recently formalized Berry phase to classical mechanics.

The angular velocity tensor is a skew-symmetric matrix defined by:

In general relativity, Lense–Thirring precession or the Lense–Thirring effect is a relativistic correction to the precession of a gyroscope near a large rotating mass such as the Earth. It is a gravitomagnetic frame-dragging effect. It is a prediction of general relativity consisting of secular precessions of the longitude of the ascending node and the argument of pericenter of a test particle freely orbiting a central spinning mass endowed with angular momentum $.$

References

↑ Эффект Джанибекова (гайка Джанибекова), 23 July 2009 (in Russian). The software can be downloaded from here
↑ Poinsot (1834) Theorie Nouvelle de la Rotation des Corps, Bachelier, Paris
↑ Derek Muller (September 19, 2019). The Bizarre Behavior of Rotating Bodies, Explained. Veritasium. Retrieved February 16, 2020.
↑ Levi, Mark (2014). Classical Mechanics with Calculus of Variations and Optimal Control: An Intuitive Introduction. American Mathematical Society. pp. 151–152. ISBN 9781470414443.
↑ Efroimsky, Michael (March 2002). "Euler, Jacobi, and Missions to Comets and Asteroids". Advances in Space Research. 29 (5): 725–734. arXiv: astro-ph/0112054 . Bibcode:2002AdSpR..29..725E. doi:10.1016/S0273-1177(02)00017-0. S2CID 1110286.

External links

Dan Russell (5 March 2010). "Slow motion Dzhanibekov effect demonstration with table tennis rackets" . Retrieved 2 February 2017– via YouTube.
zapadlovsky (16 June 2010). "Dzhanibekov effect demonstration" . Retrieved 2 February 2017– via YouTube. on Mir International Space Station
Viacheslav Mezentsev (7 September 2011). "Djanibekov effect modeled in Mathcad 14" . Retrieved 2 February 2017– via YouTube.
Louis Poinsot, Théorie nouvelle de la rotation des corps, Paris, Bachelier, 1834, 170 p. OCLC 457954839 : historically, the first mathematical description of this effect.
"Ellipsoids and The Bizarre Behaviour of Rotating Bodies". YouTube . 24 July 2020. - intuitive video explanation by Matt Parker
The "Dzhanibekov effect" - an exercise in mechanics or fiction? Explain mathematically a video from a space station,
The Bizarre Behavior of Rotating Bodies, Veritasium

This page is based on this Wikipedia article
Text is available under the CC BY-SA 4.0 license; additional terms may apply.
Images, videos and audio are available under their respective licenses.

[1] Эффект Джанибекова (гайка Джанибекова), 23 July 2009 (in Russian). The software can be downloaded from here

[2] Poinsot (1834) Theorie Nouvelle de la Rotation des Corps, Bachelier, Paris

[3] Derek Muller (September 19, 2019). The Bizarre Behavior of Rotating Bodies, Explained. Veritasium. Retrieved February 16, 2020.

[4] Levi, Mark (2014). Classical Mechanics with Calculus of Variations and Optimal Control: An Intuitive Introduction. American Mathematical Society. pp. 151–152. ISBN 9781470414443.

[5] Efroimsky, Michael (March 2002). "Euler, Jacobi, and Missions to Comets and Asteroids". Advances in Space Research. 29 (5): 725–734. arXiv: astro-ph/0112054 . Bibcode:2002AdSpR..29..725E. doi:10.1016/S0273-1177(02)00017-0. S2CID 1110286.

[1]

[2]

[3]

[4]

[5]