The Twisting Tennis Racket
The Twisting Tennis Racket
The Twisting Tennis Racket
1, 1991
KEY WORDS: Euler equations; rigid body motion; Euler angles; Eulerian
wobble.
1980 AMS(MOS) MATHEMATICS SUBJECT CLASSIFICATIONS: 70E15,
58F05.
1. I N T R O D U C T I O N
The classical treatments of the dynamics of a tennis racket about its inter-
mediate axis fail to describe a remarkable aspect of its motion which is
revealed in the following experiment. Mark the faces of the racket so that
they can be distinguished. Call one rough and the other smooth. Hold the
racket horizontally by its handle with the smooth face up. Toss the racket
into the air attempting to make it rotate about the intermediate axis
(namely, the axis in the plane of the face which is perpendicular to the
handle). After one rotation, catch the racket by the handle. The rough face
will almost always be up! In other words, the racket typically makes a half-
twist about its handle.
The experiment above was shown to one of us (R.C.) by Professor
W. Burke of the University of California at Santa Cruz. The twisting
phenomenon seems to be new. It is not mentioned in a recent article on the
Eulerian wobble (Colley, 1987), in general texts on classical mechanics
(Arnol'd, 1978; Goldstein, 1950; Landau and Lifschitz, 1976), or in
~'2
Fig. 1. The principal axes of a tennis racket.
dM~=dt - (~2 - l l
dM2=(%_ 1)
dt \11 ~ MIM3 (1)
dM3=dt - (~ - ~) M1M2
Here I1, 12, 13 are the principal moments of inertia of the racket, which we
assume satisfy
0<11 <12<13 (2)
11 '~I2 (3)
and
Ilq-12~I3 4 (4)
4 For a standard tennis racket such as the Wilson T-2000 the values of I1,/2, and 13
are I1 =0.00121 kg-m 2, I2=0.01638 kg-m 2, and 13=0.01748 kg-m 2 as measured by Brody
(1985). These values are used in all of our numerical examples.
70 Ashbaugh, Chicone, and Cushman
/ ( 1 3 -- I2)(M 2 -- 211 E)
Zl=/V ~ +c=Bt+c
/I2(213E-M z) /I2(M2-211 E)
A2= V I3_12 ' A'2= V L ~ '
t M x/(I3 - I2)(I 2 - I1)
722= ~ 2 -- 1113 + c = Bot + c
:/I3(M2--211 E) M/I3(I2-I1)
A3 V ' A;=
^ A
). ^e 3 C C
^
e2
Fig. 2. Solutions of Euler's equations for a fixed value of energy and angular momentum.
and
m = k 2 = (12 - 11 )(213 E - M 2) 1
(I3-I2)(M2-211E) rn'
If 6 < O, then
1, if M1 > 0
e= -1, if M1 < 0
They are the stable and unstable manifolds of the hyperbolic equilibrium
points B: (0, 2x/~2E, 0) and B': (0, - 2x/~-E2E, 0) (see Fig. 2). The choice of
e2 = _1 determines whether we are on the stable or unstable manifold of
B, while the choice of el determines which of the two branches of these
manifolds we are on.
The Eulerian picture of the m o t i o n of the tennis racket has two
short-comings. The first one is that it is based on the noninertial frame
(el, e2, e3) which is corotating with the body. To remedy this, we recall
that M is a constant of the m o t i o n of the tennis racket when referred to an
inertial frame (X, Y, Z) fixed in space. Thus, a p p a r e n t motions of M along
72 Ashbaugh, Chicone, and Cushman
The full Euler equations for Euler angles I are obtained by first solving the
equations
M1 = ~ sin 0 cos ~, - 0 sin
11
M2 = ~ cos 0 + ~b
12
Z /x
e3 ^
\N
II
Fig. 3. Two choicesof Euler angles. N is the line of nodes: the intersection of the ~1-~3plane
in 1 (the E2 k3 in H) with the X-Y plane
The Twisting Tennis Racket 73
for the angular velocities 0, ~, ~ about N, 82, and the Z-axis, respectively,
and then using (9). This gives
- - = M - s i n 0 sin ~p cos ~p
dt
m 2
12
&b=MsinZO+)TC~ ~' (12)
dMOd-t( ~ Msin2
OI2 Mc~c ~ ~ )
[see Goldstein (1950) and Landau and Lifschitz (1976) for more details].
since ~(t)= 17z/2-O(t)l. Therefore we need only look for the maximum
value of IMl(t)l. The argument now breaks up into three cases depending
on whether 6 is positive, negative, or zero.
Suppose that 6 > 0. Then from the second column in Table I and the
fact that Icn(t, m)l ~< 1, we find that the maximum value of sin ~(t) is
A1 1 / I I ( 2 1 3 E - M 2)
(14)
M - M ~/ 13-11
which, after a bit of manipulation, gives (13).
Suppose that 6 < 0. Then from the fourth column in Table I and the
fact that 1 - m ~<dn(t, m) ~< 1, we find that the maximum of sin ~(t) is again
given by (14).
Suppose that 6 = 0. Then from the third column in Table I and the fact
that sech(t) ~< 1, we find that the maximum of sin e(t) is given by
A'~ /11(13-12)
M2 = 212 E and /3 ~ 11 + / 2
~~--1 ~ I1
9 2 2 ~ -
using (3). Therefore c~(t) is small for a tennis racket whose principal
moments of inertia satisfy (2), (3), and (4). 5
Let fl(t) be the angle that the vector ~2 makes with the Z-axis (in Euler
angles I this is 0). Then the amount of twist that the tennis racket makes
about its handle in time t is
5 For the standard tennis racket, substituting M 2 = 212E and the values of Ia,/2, 13 into (13),
we find that c~(t) is at most 4.05 ~
The Twisting Tennis Racket 75
The approximation is valid because the face of the racket is initially nearly
horizontal, fl(0)~ 0. The following theorem shows that the handle rotates
nearly uniformly. Therefore, the amount of twist r(T) after one revolution
is a quantity which can be reproducibly measured. More precisely, we show
that the projection of ~l(t) on the X - Y plane moves almost uniformly.
Using Euler angles II, this means that ~b(t) (which measures the rotation
of the handle) increases nearly linearly with time. We now prove the
following.
Theorem 2. Let O(t), (J(t), t~(t) be solutions of the full Euler equations
(12). Then for all t >>.0,
M M
s-~ t <. ~(t) - ~o <. --I2 t (15)
Moreover, there are constants Jo, J1, 61 >>-O, 62 >/0 such that
dq) M ( M 2 - 2 I I E ) + ( 2 1 3 E - M 2 ) s n Z ( B t + c ' m )
d---/= I3(M 2 - 2IIE) + I~(213E- M z) sn2(Bt + c, m) (18)
from (18) we see that d4/dt is periodic of period 2K/B. Let Jo be the
average of the right-hand side of (18) over a period, then a straightforward
change of the variable of integration gives
J _Mf~C ( M Z _ Z l a E ) + ( 2 1 3 E _ M Z ) s n 2 ( s , m )
o - K Jo I3(M 2 - 211E) + I1(213E- M 2) sn2(s, m) ds
76 Ashbaugh, Chicone, and Cushman
=~(t)+jo N(t+s)ds-Jo t+
2K 2K
= ~(t/+ Jo ~ - J o t - Jo ~
=r
Let J1 be the average of ~bl over a period and put ~b2(t)= ~bl(t)- Jl. Then
~b2 is a continuous periodic function with average 0. Hence there are
constants 6~, 62/> 0 such that - 6 1 ~<~b2(t) ~<62 for all t. This establishes
(16). When 6 < 0, the argument can be carried out in essentially the same
manner. When 6 = 0, we have
dqJ=M ( I z - I 1 ) + ( 1 3 - I 2 ) t a n h 2 ( B t + c)
dt I3(12- I~) + Ia(I3- I2) tanhZ(Bt + c)
Jo - 2K 3o dt
Although the theorems in this section do not prove that the tennis
racket twists, they do show that the handle's motion is close to being
The Twisting Tennis Racket 77
uniform periodic rotation in the X-Y plane. This allows us to make the
definition of half-twist presented in Section 1 meaningful and to formulate
a simple criterion for idealizing the experiment of tossing a tennis racket,
viz., the toss is completed when the handle returns close to its starting
point, as compared to an ideal toss of a racket rotating uniformly in a
plane. Once it is known that the handle returns close to its starting point,
the twist angle is measured by the Euler angle 0 (in Euler angles I).
These facts constitute our explanation for the twisL Here (a) states
that the twist exist, (b) that the handle will have enough time to twist, and
(c) that the racket is likely to be caught after a half-twist.
We begin by discussing (a). To establish that a half-twist occurs as M
traverses the unstable manifold of B to B', we use Euler angles/. From (8)
and (9) note that the stable and unstable manifolds of B are given by
+/13(12--11)
(20)
tan ~ = --X/I1(13 12)
M
m ~ m
(21)
o
dt
78 Ashbaagh, Chicone, and Cushman
2 212/" 1113
Bo M 4-(I3-I2)(I2-I,)
212 /12+11 (using the approximation 13 ,~ I1 + I2)
~---M ~/12--11
212
-- (using 11 ~ 12)
M
COSO
B.8
~.?
B.5
0,3-
0.2:
~.~
4.1
-B2
-0.3
-OA,
-0.5~
-0,7
-8.8
-0.9
--B'6
1 . BIj , 1 ' I ' I ' I ' 1 ~ 1 ' 1 ' I ' 1 ' I ~ t
-1.8 -8.8 -ft.6 -I],4 -92 fl,I] fl,2 8.4 0.6 0,8 1.0
Fig. 4. The graph of cos 0 = - t a n h Bot. Here the characteristic time for the twist is
0.35836 sec when E = 0.32333 J. The time for one revolution of the handle is 1 sec.
The Twisting Tennis Racket 79
m l = '~ aM3
(22)
M3 = ~ bM1
where
1 1
(22) becomes
(23)
__ AM 3 M3=~ ~,
\ T }_o.,,
# ~ / 1
-axis
! I I\ --~IM3=~-b M 1
at B at B
2 in xf2 d
TN - ~ ~)
5. N U M E R I C A L RESULTS
To show for a large percentage of suitable initial conditions that the
tennis racket does perform a near-half-twist, we have simulated the tennis
racket experiment described in Section 1 on a computer by integrating the
full Euler equations (11). This allows us to simulate performing a large
number of experiments with a variety of initial conditions. We use the
values of 11,/2, and 13 of the standard racket and the value 0.32333 for E.
We choose the initial ~bo from [0, ~) and take small initial values of 00. The
total angular momentum M is then determined by
M2=2E( c~176176 +Tc~ 00 -~ sin2 ~9~- sin2
1 10~3
using (5) and (9). Finally, we choose ~bo so as to make ~1 (that is, the
racket handle) as near to the direction of the positive Y-axis as possible for
the given values of 0o and ~9o. This is accomplished by taking ~bo to be the
angle between (cos ~bo cos 0o, - s i n ~bo) and (1, 0), that is,
( - sin ~bo
~o = arctan \cos ~o cos OoJ
The Twisting Tennis Racket 81
tan - 1 F /I1(I3-I2!]~007073 r a d i a n t 4 . 0 5 ~
LVI3(/2 - L)J "
( - sin 0 cos ~b- cos 0 cos 0 sin ~, - sin 0 sin ~b+ cos ~ cos 0 cos ~b)
Finally, we checked the value of O(t) at the stopping time. If ~2(t) was
within 27 ~ of the negative Z-axis, i.e., within 15% of completing a half-
twist, then we considered that a near-half-twist had occurred; otherwise it
had not. The results of the numerical experiments are presented in Fig. 6.
Figure 6 suggests that, to a good approximation, the region of no
near-half-twist is a strip of constant width at an angle 0 = 1.5001 radians
and centered on the stable manifold through B. Assuming this to be correct
we can derive the half-width 0.004426 for the strip from our knowledge
of how the twist occurs along the unstable manifold. This compares
reasonably well with the rigorous bound 0~< 0.00019 for the no near-half-
twist region found in the Appendix. In addition, one can use the initial con-
ditions (00, 0o) = (0.0044, 0) to show by numerical integration that the
near-half-twist occurs. This confirms Fig. 6.
Moreover, we have computed the expected percentage of the time that
the twist occurs, using the initial conditions in Fig. 6, our twist criterion,
and the formula
Iql3
M1
Fig. 6. The black dots are those points (0 cos ~,, 0 sin 0) corresponding to initial conditions
(0, 0) in [-0,0.025] x [-0, 27z) where the racket makes a near-half-twist according to our
criterion.
APPENDIX
[M1M31 <~'~ M 1+
where
/I1 (13 -- I s)
(r = X//--~2
" 11)
1 1
(26)
By (15) the largest time required for the projection of the racket handle to
have completed one revolution is
2rcI3
t 0 --
M
which is the intersection of the energy ellipsoid (5) and the plane { m 2 = 0 } ,
then the solution M(t) of Euler's equations on the energy ellipsoid does not
cross the {M2 = 0} plane. Thus no near-half-twist can occur. Since
13 - / 2 ~<I2 -- 11
I 3 J \ 11 + 13 /t
Therefore if
sn00)sin0o
The Twisting Tennis Racket 85
where (0o, ~bo) are the initial values of (0, ~b). Combining (29) and (30), we
see that no near-half-twist occurs if
Msinuo< /
/211E(I3-12)
~ exp[-2~/
13/(I3-I2)(I2-I,!]
/~ _~ (3i)
This bound is rigorous when M2<212E. This covers 95% of all initial
conditions due to the small size of the acute angle sectors at B in Fig. 2.
Even when M 2 > 212Ethe corrections to this estimate will be very small.
ACKNOWLEDGMENTS
We are grateful to Howard Brody of the University of Pennsylvania,
Enoch Durban of Princeton University, and Doug Wisbey of Wilson
Sporting Goods for sharing their knowledge of the moments of inertia of
various tennis rackets with us. We would also like to thank William Burke
of the University of California at Santa Cruz and Tom Kane of Stanford
University for their helpful criticism of the first draft of this paper.
Our computer experiments were carried out on an IBM 4381-R14
using the IMSL routine DGEAR.
This research was supported by The Air Force Office of Scientific
Research under Grant AF-AFOSR-89-0078 (C.C.C.).
REFERENCES
Abramowitz, M., and Stegun, I. A. (eds.) (1964). Handbook of Mathematical Functions,
Applied Mathematics Series, Vol. 55, National Bureau of Standards, Washington, D.C.
Arnol'd V. I. (1978). Mathematical Methods of Classical Mechanics, Springer-Verlag,
New York.
Brody, H. (1985). The moment of inertia of a tennis racket. Phys. Teach. April, 213-216.
Colley, S. J. (1987). The tumbling box. Am. Math. Month. 94, 62-68. [See also (1987). Letter
to the editor. Am. Math. Month. 94, 646.]
Goldstein, H. (1950). Classical Mechanics, Addison-Wesley, Reading, Mass.
Gradshteyn, I. S., and Ryzhik, I. M. (1980). Table of Integrals, Series and Products, corrected
and enlarged ed., Academic Press, New York.
Klein, F., and Sommerfeld, A. (.1897-1910). (Yber die Theorie des Kreisels (4 vols,), B.G.
Teubner, Leipzig.
Landau, L. D., and Lifshitz, E. M. (1976). Mechanics, 3rd ed., Pergamon Press, Oxford.
Rauch, H., and Lebowitz, A. (1973). Elliptic Functions, Theta Functions, and Riemann
Surfaces, Williams and Wilkins, Baltimore.
Tricomi, F. G. (1953). Differential Equations, Blackie and Son, London.
Webster, A. G. (1920). The Dynamics of Particles and of Rigid, Elastic, and Fluid Bodies,
Stechert-Hafner, New York.
865/3/1-6