HOW TO TWIRL A HULA-HOOP

A.P. Seyranian*, A.O. Belyakov Institute of Mechanics Moscow State Lomonosov University Michurynski pr. 1, Moscow 119192, Russia *[email protected] [email protected] ABSTRACT. We consider twirling of a hula-hoop when the waist of a sportsman moves along an elliptic trajectory close to a circle. For the case of the circular trajectory, two families of exact solutions are obtained. Both of them correspond to twirling of the hula-hoop with a constant angular speed equal to the speed of the excitation. We show that one family of solutions is stable, while the other one is unstable. These exact solutions allow us to obtain approximate solutions for the case of an elliptic trajectory of the waist. We demonstrate that in order to twirl a hula-hoop one needs to rotate the waist with a phase difference lying between / 2 and . An interesting effect of inverse twirling is described when the waist moves in opposite direction to the hula-hoop rotation. The approximate analytical solutions are compared with the results of numerical simulation. I. INTRODUCTION. A hula-hoop is a popular toy a thin hoop that is twirled around the waist, limbs or neck. In the last decades it is widely used as an implement for fitness and gymnastic performances, Fig. 1. To twirl a hula-hoop the waist of a sportsman carries out a periodic motion in the horizontal plane. For the sake of simplicity we assume that the waist is a circle and its center moves along an elliptic trajectory close to a circle. Earlier the simple case when a hula-hoop is treated as a pendulum with the pivot oscillating along a line has been considered. 1,2 The stationary rotations of a hula-hoop excited in two directions have been studied by an approximate method of separate motions. 3 Similar problem on the spinner mounted loosely on a pivot with a prescribed bi-directional motion has been treated numerically and experimentally. 4 Here we derive the exact solutions in the case of a circular trajectory of the waist center and approximate solutions in the case of an elliptic trajectory. We also check the condition of keeping contact with the waist during twirling. The paper differs from our previous one 5 by new

approximate solutions found analytically and compared with the numerical simulation results.

Fig. 1. Twirling of hula-hoops in the circus performance, Moscow 2008

(photo by K. Stukalov) .

II. MAIN RELATIONS AND EXACT SOLUTIONS. We assume that the center O ' of a sportsmans waist moves in time t according to the elliptic law x = a sin t , y = b cos t with the amplitudes a , b and the excitation frequency > 0 , Fig 2. The equations of motion in the waist-fixed coordinate system with allowance for linear viscous damping take the form

&& & I c + k = FT R && m( R r ) = m ( &&sin + && cos ) + FT x y

& m( R r ) 2 = N + m ( && cos && sin ) x y

(1) (2) (3)

Fig.2. A hula-hoop with the radius R twirling with the angle around a circular waist (shaded) with the radius r . The center O ' of the waist moves along the elliptic curve x = a sin t , y = b cos t with the fixed center O . where is the rotation angle around center of mass C , I c = mR is the central moment of inertia of the hula-hoop, is the angle between axis x and radius CO ' , FT is the friction force, m and R are the mass and radius of the hula-hoop, k is the coefficient of viscous friction, r is the radius of the waist, and N is the normal force. Equation (1) describes rotation of the hula-hoop around C , while (2) and (3) are equations of motion of the hulahoop in the longitudinal and transverse directions to the radius CO ' . Assuming that slipping at the point of contact is absent we obtain the kinematic relation
2

& & ( R r ) = R

(4)

We exclude from equations (1) (2) the force FT and with relation (4) obtain the equation

2 ( a sin t sin + b cos t cos ) k && & + + =0 2 mR 2 2( R r )

(5)

From equation (3) we find the normal force and imply the condition N > 0 as (6) which means that the hula-hoop during its motion keeps contact with the waist of the sportsman. We introduce new time = t and non-dimensional parameters

& m ( R r ) 2 + m 2 (a sin t cos b cos t sin ) > 0

a b k ,= , = 2( R r ) 2( R r ) 2mR2

(7) (8)

Then equation (5) and inequality (6) take the form

(9) where the dot means differentiation with respect to the time . These equations contain three non-dimensional parameters: the damping coefficient , and the amplitudes and along the axes x and y. Relation between the amplitudes and determines the form of ellipse the trajectory of the waist center. For = 0 the trajectory is a line, and for = it is a circle. For the equal excitation amplitudes = and arbitrary damping coefficient equation (8) takes the form && & + + cos ( ) = 0 (10)

&& & + + sin sin + cos cos = 0 & 2 + 2( sin cos cos sin ) > 0

= + (11) with the constant initial phase given by the equation cos = (12) Therefore, solution (11) exists only under the condition (13) From equation (12) we find = arccos ( / ) + 2 k , k = 0,1, 2,... (14) Solutions (11), (14) correspond to the rotation of the hula-hoop with the constant angular velocity equal to the excitation frequency . Let us investigate the stability of the obtained solutions. For this purpose we take the angle in the form = + + ( ) where ( ) is a small

which has the exact solution 5

quantity, and substitute it into equation (10). Taking linearization with respect to and with the use of (12) we obtain a linear equation && & + sin = 0 (15) Without loss of generality we assume > 0 since the case < 0 is reduced to the previous one by the time transformation ' = + in equation (8). According to Lyapunovs theorem on the stability based on a linear approximation 6 solution (11), (14) is asymptotically stable if all the eigenvalues of linearized equation (15) have negative real parts. From Routh-Hurwitz criterion 6 due to the assumption > 0 we obtain the stability conditions as > 0 , sin < 0 (16) From conditions (16) and relation (14) we find that for 0 < < solution (11) with = arccos ( / ) + 2 k , k = 0,1, 2,... (17) is asymptotically stable, and solution (11) with = arccos ( / ) + 2 k , k = 0,1, 2,... (18) is unstable. For negative damping < 0 both solutions (11), (17) and (11), (18) are unstable. It is necessary to verify for the exact solutions (11), (14) the condition of twirling without losing contact (9) which takes the form 1 2 sin > 0 (19) Due to (16) and the assumption > 0 condition (19) is satisfied for the stable solution (11), (17). Thus, solution (11), (17) under the condition 0 < < provides asymptotically stable twirling of the hula-hoop with the constant angular velocity without losing contact with the waist of the sportsman. The phase of this solution belongs to the interval / 2 plus 2 k to both sides of this inequality, and for vanishing damping +0 the phase tends to / 2 . Below we will show that this phase inequality also holds for the approximate solutions. This is the way how to twirl a hula-hoop! III. APPROXIMATE SOLUTIONS. For close but not equal amplitudes we introduce notation = ( ) / 2 , = ( + ) / 2 . For the sake of simplicity we assume which means, 0 0 . Then equations (8), (9) take the form

(21) Taking as a small parameter we apply perturbation method assuming that solution of (20) can be expressed in a series s ( ) = * ( ) + o( ) , * ( ) = + 0 + 1( ) , (22) where * ( ) is the first approximation function. After substitution of series (22) in (20) and grouping the terms by equal powers of we derive the following chain of equations

& 2 2 sin ( ) + 2 sin ( + ) > 0

&& & + + cos ( ) = cos ( + )

(20)

0 : + cos 0 = 0 && & 1 : 1 + 1 sin(0 )1 = cos(0 + 2 )

Taking solution of the first equation (23) 0 = arccos ( / ) + 2 k

(23) (24)

corresponding to the stable zero order approximation ( ) = + 0 , we write the second equation (23) as

&& & 1 + 1 + 2 2 1 = cos(0 + 2 )

It has a unique periodic solution 1 ( ) = C sin ( 2 + 0 ) + D cos ( 2 + 0 ) where

(25) (26)

C=

+ 3 8 + 16
2 2 2 2

, D=

4 + 2 2

2 + 3 2 8 2 2 + 16

(27)

We see that the approximate solutions * ( ) with (24), (26), (27) differ from the exact solutions (11), (17) by small vibrating terms and phase shift. Note that the approximate solutions were obtained with the assumption that the excitation amplitudes and damping are not small. To find the stability conditions for solution (22), (24), (26) we take a small variation to the solution of (22) = s + u and substitute this expression into (20). After linearization with respect to u and keeping only terms of first order we obtain a linear equation && & u + u + [ (sin 0 + 1 cos 0 ) + sin(2 + 0 )]u = 0 (28)

where 0 is given by expression (24). Equation (28) can be written in the form of damped Mathieu-Hill equation as && & u + u + [ p + (2 )] u = 0 (29) where (30) Then the stability condition (absence of parametric resonance at all frequencies p ) is given by the inequalities 6

p = 2 2 , = ( C + 1)sin(2 + 0 ) + D cos(2 + 0 )
2 ( C + 1) 2 + 2 D 2

<

(31)

with C and D defined in (27). This is the inequality to the problem parameters , and . The condition of twirling without losing contact (21) takes the following form

1+ 2 2 2 < 2

2 + 3 2 8 2 2 + 16 2 + 8 2 12 2 2 + 36

(32)

Conditions (31), (32) imply restrictions to , i.e. how much the elliptic trajectory of the waist center differs from the circle.
=1.4, =1, =1.2, =0.2, =1, (0) =2.5389, d/d(0) =0.90802 1.1 numerical asymptotic

angular velocity d/d

0.9 0

time

2*pi

Fig. 3. Comparison between the approximate analytical and numerical results.

In Fig. 3 the approximate analytical solution is presented and compared with the results of numerical simulation for the case when the excitation amplitudes , and damping coefficient are not small. IV. SMALL EXCITATION AMPLITUDES AND DAMPING. It is interesting to consider the case when the excitation amplitudes and damping coefficient are small having the same order as . Then we introduce new % parameters = / and % = / and assume that the solution has the form ( ) = + 0 ( ) + 1( ) + o ( ) , (33)
where is the angular velocity of rotation, and the functions 0 ( ) , 1 ( ) are supposed to be limited in time. We substitute expression (33) into (20) and equating terms of the same powers of obtain the following equations

&& 0 : 0 = 0 (34) && %& % % 1 : 1 = cos(0 + + ) cos(0 + ) 0 From the first equation (34) we get 0 = const . Then the second equation (34) can have non-growing solutions only when takes the values -1, 0, 1.
Thus, besides clockwise rotation = 1 we have also counterclockwise rotation = 1 , and no rotational solution = 0 . The letter case is not interesting, so we omit it. For clockwise rotation = 1 from equations (33) and (34) we obtain in the first approximation the solution ( ) = + 0 + 1( ) , (35)

1 cos 0 = (36) 4 First expression (36) is the special case of (26), (27) for = = 0 . To verify the stability conditions for solution (35) we use damped Mathieu-Hill equation (29) and for the case of small damping and excitation amplitudes get > 0 , sin 0 < 0 .6 These conditions are similar with inequalities (16) derived for the undisturbed exact solution. Thus, the stable solution (35) with 0 = arccos ( / ) + 2 k exists for 0 < < . For solution (35) condition (21) of keeping contact in the first approximation reads

1 ( ) = cos (0 + 2 ) ,

1+ 2 2 2 < (37) 3 For counterclockwise rotation = 1 we obtain in the first

approximation the solution

( ) = + 0 +

cos (0 2 ) ,

cos 0 =

with the stability conditions > 0 , sin 0 > 0 . Thus, the stable counterclockwise rotation (38) with 0 = arccos ( / ) + 2 k exists for

(38)

0 < < . For this case condition (21) takes the form similar with (37)

Inequalities (37), (39) are satisfied for sufficiently small , and . Stable clockwise and counterclockwise rotations (35), (36) and (38) coexist if 0 < < min{ , } . Such possibility is demonstrated in the circus performance, see Fig. 1. Coexisting clockwise and counterclockwise rotations are illustrated in Figs. 4 and 5.

1+ 2 2 2 < 3

(39)

Fig. 4. Stable twirling of the hula-hoop for the cases: a) direct twirling b) inverse
twirling.

Fig. 4 demonstrates phase difference for stable rotations in both directions.

=0.1875, =0.125, =0.15625, =0.03125, =0.025, =1.7301, d/d =0.98403 numerical asymptotic 1.01 0.95 =0.1875, =0.125, =0.15625, =0.03125, =0.025, =2.4696, d/d =0.95212 numerical asymptotic

angular velocity d/d

angular velocity d/d

2*pi

0.99

1.05

time

time

2*pi

Fig. 5. Comparison between the approximate analytical and numerical results for small excitation amplitudes and damping coefficient for clockwise and counterclockwise rotations. In Fig. 5 the approximate analytical solutions for rotations in both directions are presented and compared with the results of numerical simulation for the case of small excitation amplitudes , and the damping coefficient . The values of , correspond to the dimensional parameters a = 15 cm, b = 10 cm, r = 10 cm, R = 50 cm. V. CONCLUSION. We have derived simple explicit relations showing how to twirl a hula-hoop stably in both directions without losing contact with the waist.

REFERENCES
1. T.K. Caughey, Hula-hoop: an example of heteroparametric excitation, Am. J. Phys. 28 (2), 104-109 (1960). 2. T. Horikawa and Y. Tsujioka, Motion of hula-hoop and its stability, Keio Science and Technology Reports, 40 (3), 27-39 (1987). 3. I.I. Blekhman, Vibrational Mechanics, (Fizmatlit, Moscow, 1994). 4. J.F. Wilson, Parametric spin resonance for a spinner with an orbiting pivot. Int. J.Non-linear Mech, 33 (2), 189-200 (1998). 5. A.O. Belyakov and A.P Seyranian, The hula-hoop problem, Doklady Physics 55 (3), 99-104 (2010). 6. A.P.Seyranian and A.A.Mailybaev, Multiparameter Stability Theory with Mechanical Applications, (World Scientific, New Jersey, 2003).