A Mathematical Model for Top Nutation Based on Inertial Forces of Distributed Masses

The main property of gyroscopic devices is maintaining the axis of a spinning rotor, a mathematical model formulated on the principle of the change in the angularmomentum.This principle is used formathematicalmodeling of themotions of a top at known publications. Nevertheless, practical tests of gyroscopic devices do not correspond to this analytical approach. Recent investigations have demonstrated that the origin of gyroscope properties is more complex than that represented in known publications. The applied torque on a gyroscope produces internal torques of the spinning rotor based on the action of the several inertial forces. These forces are the centrifugal, Coriolis, and common inertial forces as well as the change in the angular momentum generated by themass elements and center-mass of the spinning rotor.The action of these torquesmanifests itself in the resistance and precession torques of the gyroscopic devices. These inertial torques act simultaneously and interdependently around two axes and represent the fundamental principles of the gyroscope theory.The new inertial torques enable deriving mathematical models for the motions of well-known top that is the simplest form of gyroscopic devices. The novelty of the work is mathematical models for the motions of the top based on action of eight inertial forces acting around its two axes.The obtainedmathematical models for the top nutation and self-stabilization are represented in terms ofmachine dynamics and vibration analysis.The new analytical approach formotions of the well-balanced top and top with eccentricity of the center-mass definitely responds to the practical results.


Introduction
L. Euler first laid out the mathematical foundations for the gyroscope theory in his work on the dynamics of rigid bodies back in 1765.Since Industrial Revolution time, I. Newton, J-L.Lagrange, L. Poinsot, J. L. R. D' Alembert, P-S.Laplace, L. Foucault, and other brilliant scientists have investigated, developed, and added new interpretations of the gyroscope effects, which are in full display in the rotor's persistence in maintaining its plane of rotation.The applied theory of gyroscopes emerged mainly in the twentieth century in numerous publications that described the gyroscope effects [1].The gyroscope properties feature many engineering mechanisms and devices with rotating parts, which thus need to be computed for their proper functioning [2,3].Gyroscope effects are used in numerous gyroscopic devices in aerospace engineering, as well as on ships and other industries [4,5].The most basic textbooks of classical mechanics have chapters on the gyroscope theory [6,7] and consider motions with vibrations analysis [8].Numerous and valuable publications have dedicated themselves to gyroscope effects and their applications in engineering [9][10][11][12].Intuitively and without mathematical models, some researchers have noted that, on a gyroscope, the inertial forces that manifest the gyroscope effects are acting.However, in known publications, mathematical models for the gyroscope effects do not seem to match their practical applications in gyroscopic devices [13][14][15][16][17]. Therefore, researches have spawned artificial terms such as gyroscope resistance and gyroscope couple, as well as fantastical properties that contradict rules of classical mechanics.It is for this reason that the gyroscope theory still attracts many researchers who seek to discover true gyroscope theory [16,18].
Mechanically, a gyroscope is a spinning disc in which the axle is free to assume any orientation.The simplest gyroscope is represented by a top toy that is one of the most remarkable and widely recognized toys in the world.They still attract attention through their astonishing behavior and unusual gyroscopic properties.The motions of top toy have been described analytically in numerous publications and with complex numerical modeling [19] based on Lagrangian dynamics on Euler coordinates that are solved with computer software [20].However, all publications contain approximations, assumptions, and simplifications and explain gyroscope effects by the physical principle of the spinning rotor's angular momentum that generates the precession torque and motions.Gyroscope effects still represent a problem that remains to be solved.The origin of gyroscope effects is more complex than those represented in the theories known to date.Recent investigations into the physical principles of gyroscope motions demonstrate four classical inertial forces acting upon a spinning rotor to generate the gyroscope effects.Research shows that centrifugal, common inertial, and Coriolis forces, in addition to change in the angular momentum of spinning rotors, are the basis for all gyroscope effects and properties [21][22][23][24].Applied to the gyroscope, an external torque generates several torques based on the action of the aforementioned forces.In turn, the centrifugal and Coriolis forces of the rotating mass elements generate a resistance torque that counteracts the inclination of the rotor's location.Common inertial forces, rotating mass elements, and change in the angular momentum of a spinning rotor generate a precession torque.All torques are interrelated and occur simultaneously.New mathematical models for acting torques accurately describe the physics of gyroscope effects and are validated by tests [21].Table 1 represents four different torques generated by the inertial forces of the mass elements and centermass of the spinning rotor under the action of the external torque that was applied to the gyroscope.The action of these torques should be considered for mathematical models of the gyroscopic devices and the top's motions and nutation.
In Table 1,  . is a resistance torque generated by the centrifugal ( ct. ) and Coriolis ( cr. ) forces acting around axis ;  = ( 2 /2) is the mass moment of inertia of the disc around axis  or  [6];  . is a precession torque generated by the common inertial forces ( in. ) and change in the angular momentum ( am. ) acting around axis ; and  is the angular velocity of a spinning rotor.
When applied to a gyroscope, an external torque generates angular velocities of precessions around two axes and activates simultaneously and interdependently eight internal torques.The equations of these torques (Table 1) demonstrate that the centrifugal, common inertial, and Coriolis forces of the spinning rotor are real active physical components, just like the change in the angular momentum.The new equations for the internal torques demonstrate their proportional dependence on a rotor's mass moment of inertia and angular velocity of the spinning rotor, as well as on the angular velocity of its precession.New mathematical models for these torques clearly describe the physics of gyroscope properties and explain the unusual motions of gyroscopic devices.The aim of this paper is to represent mathematical models for the motions of the well-balanced top and with its eccentricity of the center-mass.This analytical approach is based on the action of centrifugal, common inertial, and Coriolis forces as well as the change in the angular momentum generated by the mass elements and center-mass of a spinning rotor.New mathematical model for the top motions describes the physics of forces acting on the top and its motions that represents the novelty of the work.The obtained mathematical models for top nutation and selfstabilization are represented in terms of machine dynamics and vibration analysis.

Materials and Methods
The new mathematical models for gyroscope motions will consider the simultaneous and interdependent actions of internal torques (Table 1).The external load torque generates eight internal resistance and precession torques acting around two axes and reveals new gyroscope properties [21].The interdependent action of internal torques is used to model the behavior of different gyroscopic devices, including the simplest top.
Observation of a fast-spinning top shows that it preserves the vertical position of the axis, remains steady on the supporting pivot, and avoids tilting or falling to the ground.If the axis of a spinning top toy is inclined from the vertical axis under the action of an external torque, the axis starts to describe the vertical circular cone of the precessed motion.Then the spinning top stabilizes itself and its spin axis goes to a vertical position.
However, the axis of a spinning top has free oscillation of the periodic complex motion, which is called nutation, and reveals itself as fast shivering of the precessing axis.This phenomenon of the top's gyroscope nutation represents small but fast oscillations of the spinning rotor's axis at about its mean position.The nature of gyroscope nutation is represented by the following sources: the center-mass is displaced from the axis of rotation, the shocked displacement of the center-mass is generated by the action of an external force and the imperfect surface geometry of the tip of the top's leg, and so forth.
The action of the center-mass, having been displaced from the spinning top, represents the continuing existence of a disturbing force that generates a forced oscillation.Any free oscillation of the spinning top is connected with an angular acceleration that eventually ceases because of a loss of energy.The amplitudes of free oscillated motions are very small in the case of a rapidly rotating top.Because of the inevitable presence of resistances, these oscillations for a well-balanced top toy are subject to asymptotically rapid decay.In most cases, the nutation of a well-balanced top toy is quickly damped by the action of inertial resistance forces and friction forces in the bearing, leaving uniform precession.
Gyroscopes with a high angular velocity represent overdamped or critically damped systems according to the theory of vibration [8].This makes it possible to neglect nutation in solving most engineering problems and to formulate a gyroscope theory that takes into account only the precessions.However, for a gyroscope with low angular velocity, nutation is a regular process and represents a simple damped system.
The nutation process of the top can be considered in terms of machine dynamics, particularly in terms of vibration analysis [8].The acting forces on the top are its weight and the force generated by centrifugal forces of the top toy's centermass, torques generated by the inertial forces of the rotating mass elements, and change in the angular momentum.The simulation of a free nutation process is presented by modeling it with a damping process and energy losses in the oscillation of a not so well-balanced spinning top.
Vibration analysis often takes into account the energy losses by means of a single factor called the damping factor.The nutation model of the spinning top consists of a not so well-balanced mass and a damper that represents the action of the top's weight and several internal torques.In the case of a gyroscope's free nutation process, the action of the single shock type torque on a spinning top toy represents one period of an oscillation with a short time for a single cycle.The damping resistance torques decay the action of one shock torque, resulting in the inclination of only the spinning top's axis.For well-balanced gyroscopic devices and a top with high angular velocities, the process of the damped oscillation is not displayed.
For a top that possesses an eccentric rotating mass, that is, the center-mass being displaced on some small distance from the top's axis, nutation is a regular process.The centrifugal force of the offset mass is asymmetric, which causes an angular displacement of the top's shaft.The spinning top toy is constantly being displaced and moved by these asymmetric forces.The eccentric rotating mass of the top represents harmonic nutation, which means that the top is forced to nutate at the frequency of excitation [8].A plot using the nutation amplitude along one axis and the forcing frequency along the other axis is then described as a performance or response curve for the system.The eccentric rotation of the mass at about the top's axis is modeled as a sinusoidal wave, as shown in Figure 1.Here, the function  =  sin  denotes the displacement of the eccentric mass that represents the excitation input, where  is the eccentricity of mass location.The frequency of this sine wave is also the frequency at which nutation occurs.
The analysis of a gyroscope motions and nutation is conducted using the example of a top that is not well balanced and is thus tilted (Figure 2).The motions of the spinning top are considered at around the point of support .If the axis of a top is adjusted on the angle  to the horizontal position and released, then, under the influence of the force of gravity, its axis should tilt in the direction of the action of this force.As a result, the top's axis will begin to precess at about the vertical and horizontal positions.The precessions of the top are also accompanied by visually undetectable nutational oscillations, which decay rapidly according to the action of the internal torques.Because of the internal torques' resistance forces, the proper rotation of the top gradually slows down, while the precession velocities correspondingly increase.When the angular velocity of the top becomes smaller, the resistance forces of the internal torques become weaker, and the top becomes unstable and falls.In the case of a slowly rotating top, the nutational oscillations are noticeable and manifest a substantial change in the pattern of the top toy's axial movement.In this case, the end of the top's axis can describe a clearly visible wavy or looped curve alternately departing from the vertical position.In the case of a highly rotating top, the nutational oscillations can decrease and manifest the stabilization of the top's axial movement.The top's axis can approach the vertical position.
The motions of a not so well-balanced spinning top are complex.The equation of nutation is formulated by the action of several forces on the top.These forces represent the top's weight, the centrifugal forces of eccentric mass and centermass, the centrifugal and Coriolis forces generated by the rotating mass elements, the common inertial forces, and the rate of change in the angular momentum of the spinning toy.These forces produce the torques acting on the top.These are demonstrated in Figure 2 and formulated by the equations presented in Table 1.
The centrifugal force of the eccentric mass generates the variable torque represented by the following equation [7]: where  =  2 is the centrifugal force generated by the rotated eccentric mass ,  is the length of the top's leg,  is the tilt angle of the top's axis,  is the angle that should be used for calculating the maximum magnitude, and   and   are the torques acting around axes  and , respectively.The weight of a top generates torque, for which the equation is as follows:  =  *  =  cos . ct. is the torque generated by the centrifugal force of the rotating top's center-mass m around axis oy.Other acting torques are presented in Table 1.
The parameters defined above allow for the formulation of a mathematical model for a top's motion around axes   and  in Euler's form.The torques acting on the top are similar to torques acting on the gyroscope suspended from a flexible cord [22,23].Both models consider the gyroscope with one-side support to be equal except for the added torque generated by the centrifugal force of the eccentric rotating center-mass.As such, the mathematical model for a top's motion around axes  and  is represented by the following system of equations: where   = ( 2 /4) +  2 is the top's mass moment of inertia around axis  [6], while the sign (±) denotes (+) and (−) as the positive and negative directions of the action, respectively, for the torque generated by the eccentric mass.
The torques generated by the centrifugal forces of the rotating top's center-mass around axis  are defined by the following equation: where  ct. =  cos   2 is the centrifugal force of the top's center-mass rotating around axis  and   is the angular velocity of the top around axis ; the other components are as specified above.
The methodology towards a solution for (3) and ( 4) by adding the ratio of the angular velocities of a gyroscope around axes   = (  ) is presented in Usubamatov [22,23].Substituting the defined parameters ((1), (2), and ( 5) and Table 1) into ( 3) and ( 4), transformation yields the equations of the top's motions as follows: where all parameters are as specified above.Solving ( 6), (7), and (8) yields the maximum and minimum angular precessions   and   which allow for defining the nutation amplitudes around axes  and , respectively.Changes in the angular velocities of precessions around two axes make it possible to depict the trajectory of a top's centermass nutation.The mathematical model for the motions of a tilted top is derived from ( 6), (7), and (8).Substituting the expression   from ( 8) into (6), transformation yields the following equation: where all parameters are as specified above.Around axis , =  2 /2 0.625 × 10 −5 Around axes  and  of the centre mass,  =  2 /4 0.3125 × 10 −5 Around axes  and  at the point of support,  =   =  2 /4 +  2  1.1125 × 10 −5 Self-Stabilization.The practical observation of a tilted, wellbalanced spinning top with a high angular velocity proves its capacity for self-stabilization.The axis of the tilted spinning top goes to the vertical position by the action of the internal torques, whose magnitudes are bigger than torques generated by the top's weight.The necessary condition for a top's selfstabilization is formulated in the following approaches.The expression of the variable torque generated by the eccentric mass is omitted from (9).The torques generated by the top's weight and by the centrifugal forces of the rotating top's center-mass around axis  are acting in an anticlockwise direction.By contrast, the torques based on the action of the centrifugal and Coriolis forces as well as the factored change in the angular momentum are acting in a clockwise direction.The torques acting in a clockwise direction account for the self-stabilization of the top, as expressed in the following equation: Analysis of (10) shows that equilibrium of the acting torques depends on two main components, that is, the angular velocity of the top and the angle of its inclination.The angular velocity of precession depends on the first one.The stabilization process is intensive when the radius of gyration and angular velocity of spinning both have big value in addition to the length of the top's leg being short, that is, the center-mass of the top is located towards the lower part.These components explain the unstable spinning of a top and other spinning objects with a long leg and small radius.This property is confirmed by ordinary tests.A detailed analysis of (10) demonstrates that the magnitude of its right component is always bigger than its left one except for the very low value of the top's angular velocity that does not support its rotation.This top property is validated by practical observation.The low critical angular velocity of the top, which leads to a loss of self-stabilization, is defined when both the right and left sides of (10) are equal.

Working Example
The example considers the nutation process of a disc-type top whose data is presented in Table 2, whose formulas of the mass moment of inertia are represented in Hibbeler [6].
The center-mass of the top is located at the plane of the disc (Figure 2).The top initially possesses a displaced mass, while spinning around a vertical axis is conducted with nutation.
The torque generated by the top's weight is as follows: The torque generated by the centrifugal force of the top's center-mass rotation around axis oy is as follows: The maximal torque generated by the centrifugal forces of the eccentric rotating center-mass around axis ox is as follows: where all parameters are as specified above.Substituting the defined parameters (Table 2, ( 11), ( 12)) into ( 9), transformation yields the following differential equation: ( Simplification and transformation of ( 14) yield the following equation: Separation of variables and transformation for the differential equation ( 15) yield the following expression: 2  − 1.354853169  + 0.027488981 ± 0.047442787 = 3320.982248327. ( The denominator on the left side of ( 16) is represented by an equation with a known solution: where the solutions are   = /2 ± √(/2) 2 −  and yield the following results:  = 1.354853169,  = 0.047442787(0.579413283± 1).
Substituting the new expressions into ( 16) and then converting to integral forms with definite limits yield the following expression: Integral equation ( 18) is transformed and represented by the following equation: Left integral of ( giving rise to the following: The next transformation yields the following result: The right component of (22) has small value of a high order that can be neglected.Substituting expressions  and  into simplified equation (22), transformation yields the following result: Mathematical Problems in Engineering 7 Solving (23) yields the maximal and minimal value of the precession angular velocity for the gyroscope around axis :  .max = +0.0577694142rad/s,  .min = −0.0145709463rad/s, (24) where the signs (+) and (−) represent motion in an anticlockwise and clockwise direction around axis , respectively.Equation ( 24) enables the angular velocity of a wellbalanced top precession around axes ox and oy to be defined, as represented by the following magnitudes: The precession's linear velocities of a balanced top's centermass around axes  and  have the following magnitudes: The distance of precessed motion for the center-mass of a balanced top per nutation of half oscillation time is as follows: =    = 14.494917406 × 0.015 = 0.217 mm.
The maximal and minimal value of the linear velocity of the top's center-mass around axis ox have the following magnitudes: where the signs (−) and (+) represent linear velocity in negative and positive direction along axis , respectively.The time spent on a half oscillation is as follows: The maximal and minimal value of the amplitudes of centermass nutation around axis  as well as the location of a balanced top are as follows: The total amplitude of oscillation around axis  is as follows: =        max      +   min = 0.0044 + 0.0011 = 0.0055 mm. (33) The maximal and minimal value of the angular velocity for top's precession around axis oy are defined by (8).Substituting the parameters defined above and in Table 2 into (26), transformation yields the following results: rotation in anticlockwise and clockwise directions.Hence, the linear motions of the top's center-mass will be in negative and positive directions along the axes.This situation is expressed in the diagrams of nutation amplitudes.This diagram of the top's nutation loop is presented in Figure 3 and contains two coordinate systems.The system with asterisks is accepted for the motions of a balanced top.The figure demonstrates that the actual loop of amplitudes for top's center-mass has stretched down along the axis .This is a natural result, because any action of the centrifugal force in the direction of top toy's weight increases the amplitude.Any action in the opposite direction decreases the amplitude of nutation.
The magnitudes of precessed motion's distance (29) and the amplitudes of oscillations around axes  (37) and  (32) make it possible for diagrams of top's nutation to be depicted.(39) Deployed diagram of top's nutation (illustrated in Figure 4).The deployed diagram of top's nutation can have loops if the negative amplitude of oscillation along axis  is bigger than the distance of top's motion for the same precessing time.Equations of top's motions with the given and calculated data of nutation parameters demonstrate that the amplitude and frequency of nutation are cyclic.The cycles are different, because the centrifugal force of top's eccentric center-mass acts at the first semicircle of rotation in the direction of precessed motions and at the second semicircle in a direction opposite to that of the precessed motions.The result of this example demonstrates that the axis of a spinning top toy precesses with a rapid tremor of small amplitude.For top with a big, eccentric, rotating center-mass, the tremor can be visible; however, for small eccentrics, this can hardly be noticeable to the naked eye.
Self-Stabilization.The presented data of a tilted spinning top allow the condition for its self-stabilization to be checked.To make it simple, this is considered a well-balanced top that does not have the eccentric mass.Substituting the data defined above and in Table 2

Results
The simplest form of a gyroscope is a top, whose motions are described by known publications in terms of mathematical models based on the change in the angular momentum.The analytical study of forces acting on a well-balanced top and on a top toy with the eccentricity of center-mass has formulated new mathematical models for its motions.These models are based on action of several internal forces generated by the mass elements and center-mass of a top.The action of a top's weight produces its own internal torques that interrelate and act at one time and expresses precession motions.A top rotating eccentric center-mass manifests a nutation process.The obtained mathematical models for the top's motions enable describing physical principles of acting forces.These models of motions for a well-balanced top allowed its minimal angular velocity for self-stabilization to be defined.The diagram of a top's nutation represents the motions of the top with eccentricity of its center-mass.The new mathematical models for the motions of the wellbalanced top and with eccentricity of its center-mass definitely respond to the practical results.

Discussion
New studies of the gyroscope effects have shown that the origin of the acting forces and motions in a gyroscope is more complex.The gyroscope effects result from action of the several inertial torques generated by the centrifugal, Coriolis, and common inertial forces as well as the change in the angular momentum of the spinning rotor.The action of the new internal torques clearly describes the physics of gyroscope's motions and changes the traditional presentation of gyroscope effects.This new study enables describing gyroscope properties that were unexplainable at former time.Today new mathematical models for acting internal torques can be applied to any gyroscopic devices and manually solve all technical problems.These new analytical solutions were applied for describing the forces acting on a well-balanced top and on a top with the eccentricity of center-mass.The new physical principles enabled formulating mathematical models of top processed motions, nutation, and its ability to rotate vertically.

Conclusion
In gyroscope theory, the top's motions are one of the most complex and intricate in terms of analytical solutions.Known mathematical models for the top motions are accepted with simplifications and do not adequately express a real picture of its motions.The new mathematical models for gyroscope torques consider the simultaneous and interdependent action of several inertial forces generated by the rotating mass of the spinning rotor.As a practical application, these new physical principles for gyroscope motions were used for modeling a top's motions that include micro oscillation and self-stabilization.These mathematical models are thus distinguishable from those in well-known publications, which tend to have complex numerical modeling that does not interpret the origin of gyroscope effects.The application of new mathematical models for the top's motions effectively and clearly demonstrates physical principles of acting forces and motions.In that regard, this is also a good example of educational processes.

Nomenclature
: Linear motion along axis , that is,   or   along axis  or , respectively  .max ;  .min : Maximal and minimal linear motion along axis , respectively : Gravity acceleration : C e n t r i f u g a l f o r c e  ct. : C e n t r i f u g a l f o r c e o f a Linear velocity along axis   .max ;  .min : Maximal and minimal value of the linear velocity along axis , respectively : T i m e : Angle for calculating the maximal value of the torques' magnitude : T i l ta n g l eo fat o p ' sa x i s : Angular velocity of a top   : Angular velocity of precession around axis   .max ;   min : Maximal and minimal value of the angular velocity of precession around axis , respectively.

Figure 1 :
Figure 1: The position of a top's eccentric motion according to sinusoidal law.

Figure 2 :
Figure 2: Forces acting on a spinning top.
The distance of the center-mass's motion along axes  and  is bigger than the minimal value of the amplitude.It means that the deployed diagram of top's nutation does not have loops of oscillated motions.The following represent the halfcyclic nutation of top along the axis :  .max=   +  .max = 0.217 + 0.609 = 0.826 mm,  .min=   −  .min = 0.217 − 0.153 = 0.064 mm.

Table 2 :
Technical data of a top.
am. ,  ct ,  cr. ,  in. : Torque generated by the change in the angular momentum, centrifugal, Coriolis, and common inertial forces acting around axis , respectively  . ,   : Resistance and precession torque acting around axis , respectively   ,   : Torques generated by eccentricity of the center-mass and action around axes  and , respectively   :