Modeling of a Chaotic Gyrostat System and Mechanism Analysis of Dynamics Using Force and Energy
Incorporating both viscous friction torque and external torque, a mathematical model of a gyrostat system consisting of three rotors and a fixed outer frame is configured. This model is transformed into a Kolmogorov-type system for force analysis. The force field in the gyrostat system includes four different torques—inertial, internal, dissipative, and external. Correspondingly, four different energies are identified and the interconversion of energies is analyzed. The Casimir power being equivalent to the error power between the supplied power and the dissipative power is found and used to analyze the mechanism underlying the different dynamical behaviors. A four-wing chaotic attractor is found when the gyrostat is subject to a combination of inertial, internal, and dissipative torques as well as the addition of an external torque. The bifurcation of the Casimir power and leaps in Casimir energy level are used as indicators of the different definitive dynamics, which is demonstrated to be identical to that appearing in the state bifurcation and the Lyapunov exponent spectrum.
The dynamics of the gyrostat model have been extensively studied over the past two centuries, and they are important for applications in the gyrostat dynamics of satellites, spacecraft, aircraft, and robots. The earlier achievements concerning the equations of motion for a gyrostat and the analytical solutions of the equations governing the gyrostat dynamics were proposed by Wittenberg . Holmes and Marsden [2, 3] studied Smale’s horseshoe maps of Hamiltonian systems in two degrees of freedom with positive and negative damping and established the existence of such maps. From the research results of Holmes and Marsden , Krishnaprasad  established the Lie–Poisson structures for a dual-spin spacecraft with driven and free-spinning rotors and analyzed the asymptotic stability in the presence of damping. Liu, Tong, and coworkers [5–7] studied the global motion of an asymmetrical gyrostat with internal and external energy dissipations and the global motion of rigid bodies with small dissipative disturbance torques. By applying Melnikov’s method, the intersection of the stable and unstable manifolds of the system’s hyperbolic point was determined to highlight the ranges of the physical parameters corresponding to the occurrences of chaos . In studying the chaotic characteristics of the gyrostat system using Melnikov's method [8–14], the gyrostat system was analyzed by deriving its Melnikov equation and obtaining the transversal intersections of the stable and unstable manifolds. Usually, the Melnikov equation provides the necessary conditions for the existence of chaotic motion [15, 16].
However, these features of the model cannot identify the mechanism or explain the cause of the dynamic behavior generated and the energy exchange corresponding to the different force fields in the gyrostat system. The gyrostat motion is forced by comprehensive forces, and these forces are coupled. To find the mechanism underlying its definitive behavior, force and energy analysis is necessary. Arnold  presented a Kolmogorov model describing a dissipative-forced dynamic system or hydrodynamic instability with a Hamiltonian function. The importance is that the model decomposes the total forces (torques) into four classes: inertial, internal, dissipative, and external. A simple force (torque) corresponds to a simple behavior for the given system; a comprehensive force (torque) elicits complex behavior. If a system is transformed into the Kolmogorov model, its underlying mechanism is able to be uncovered by establishing the link between force (torque) and definitive behavior. Pasini and Pelino  transformed the Lorenz system into a Kolmogorov model, thereby providing the force analysis of the Lorenz system. Casimir energy plays a vital role in the energy analysis of rigid body systems. By the Casimir energy method of geometric mechanics, Pelino and coworkers  presented the energy cycle of the Lorenz system employing the Kolmogorov model. Qi and Liang [20, 21] studied the mechanics of the Qi four-wing chaotic system and the Qi chaotic system to perform a forcing analysis, where the variables of the systems are interpreted as angular momenta of a forced-dissipative rigid body. Furthermore, Qi and Zhang [22, 23] gave an analysis of energy cycling and bounds for the Qi and the Qi four-wing chaotic systems, the correspondence between torque and energy was found, and the links between the energetics and chaotic properties and the system definitive behaviors have been found. Qi and coworkers [24, 25] used the Kolmogorov model to perform a force and energy analysis on actual physical models such as the brushless DC motor and permanent magnet synchronous motor. In the dynamic analysis of various physical models [26–30], the bifurcation dynamics play an important role in the study of the motion mechanism of the dynamics. In the study of the plasma system model , Yang and Qi used the bifurcation of the Casimir power to evaluate the evolution of the system’s orbit.
The artificial satellite is one of the actual mechanical models of the gyrostats. The dynamics of different structural satellites, including the dual-spin satellite and the three-axis stabilized satellite, have received great attention in many studies. Based on the research of three-axis stabilized satellites [10, 32], we present a gyrostat system with a three-axis rotor. Thus, a mathematical model of the gyrostat system with the three-axis stabilized satellites as background is developed and is transformed into a Kolmogorov model. The object studied is a gyrostat system under external excitations and disturbances with the internal structural terms of the system being highly coupled. The error torque and power between the dissipative and external torques are found and used to reveal the cause of the different definitive behaviors of the system. A correlation is established between the error power and Casimir power.
The remainder of the paper is organized as follows. Sections 2–4 contain a description of the dynamic equations of the gyrostat and models for the external and disturbance torques, the Kolmogorov transformation of the gyrostat system model, and the four energy forms and their cycles. Section 5 analyzes the energy cycle between the different dynamic modes and reveals the causes of chaos and other definitive behaviors. In Section 6, the dynamic characteristics of the orbit of the system are analyzed using the bifurcation diagram, the Lyapunov exponent, and the bifurcation diagram of the Casimir power. Finally, some conclusions are presented in Section 7.
2. Description of the System Equations
The general dynamic equations of the gyrostat system required in subsequent sections are provided. Moreover, the relations concerning the disturbance torques are briefly reviewed.
2.1. Physical Model and Dynamic Equations
A rigid body consisting of rotors fixed to a main platform is often referred to as a gyrostat. We shall further specialize this configuration slightly and assume that three of the rotors are inertial and axisymmetric, each spinning about its axis with respect to the main body (see Figure 1). Euler’s dynamical equations of rotational motion for a disturbed gyrostat are where is the angular velocity vector; , and are the principal moments of inertia of the gyrostat (including rotors) in the body axis frame—they are constants depending on the shape and mass distribution of the body; , and are the constant components of the total angular momentum vector of the reaction rotors with respect to the body-fixed coordinate system; , and are components of the vector of disturbance torques; and , , and are external torques.
2.2. Disturbance Torques
The expressions specifying disturbance torques that exist in the gyrostat are ()where are the damping constants and is a numerical parameter. Assuming that any rotor damping relates only to its own axis (not to any other one) and remains constant, then, ; moreover, we set , and . is the damping coefficient with meaning that the energy is being dissipated from the system and meaning that the energy is supplied to the system . Here, means that the damping is negative, which is equivalent to a supplied force of the angular velocity along the -axis direction. In practice, to implement the negative damping, we can install a motor for rotor 2, and it rotates the axis with supplied torque .
The disturbance torque is expressed as follows: The focus here is to study the attitude motion of the gyrostat system for the particular case as studied in . Straightforwardly, as given in , we obtain from (1)-(3)
3. Force Field of the Gyrostat
Each variable or parameter in the gyrostat system has its own physical meaning. To study the mechanical properties of the gyrostat system, we introduce the Kolmogorov model and the Euler equation. Arnold  presented a Kolmogorov model describing dissipative-forced dynamical systems given in the formwhere , represents the Lie–Poisson structure of the kinetic energy part of the Hamiltonian function of a system, is a positive definite diagonal matrix, represents the dissipation (or disturbance torques), and represents the external torque. For a three-dimensional system, the Lie–Poisson bracket on functions of takes the form where , , is a Lie algebra and is the inner product. Hence, from (6) we haveIn the absence of external forces, the Euler equations for the rotational dynamics of a rigid body about its center of mass are where , is the principle moment of inertia for the group , and is the angular momentum satisfying , with being the angular velocity corresponding .
To satisfy the skew-symmetric property of for system (4), we set, , and . We introduce the reparameterizationinto the original gyrostat system (4) to obtainSettingwe then findUsing (9), the kinetic energy isSetting the internal forcing asthen the potential energy isThe Hamiltonian function isAs a Kolmogorov model, (10) becomeswhereTo study system (4) further, we rearrange (17). Substituting (9) into (17), we haveEquation (19) is equivalent to (4).
Remark 1. Through the transformation and derivation of the Kolmogorov model (17), the total coupled torque of the original system (4) is decomposed into four torques: inertial, internal, dissipative, and external. The decomposition is convenient to uncover certain simple dynamic behavior caused by the simple torque; the comprehensive behavior arising from the combined torques is to be described in Section 5.
4. Energy of the Gyrostat System
4.1. Energy Representation of the Gyrostat System
The energy decomposition is more important than the force decomposition in analyzing the dynamical behavior, because the force is a vector, whereas energy is a scalar. The former is difficult to analyze and to identify a relationship between force and definitive dynamic behavior of system, but the latter is relatively easier to quantify and therefore helpful in finding the mechanism underlying the dynamic behavior.
Normally, force and energy are in some form of coexistence. The four forms of energy (kinetic, potential, dissipated, and supplied) correspond, respectively, to the four forms of torques (inertial, internal, dissipative, and external) in the gyrostat system. There are internal energy exchanges between the kinetic and potential terms, and an internal–external energy exchange between the Hamiltonian and dissipated and supplied torques. Energy cycles take place among these four forms. However, the input energy is the external energy and the output energy is the dissipated energy for the system.
From (13) and (15), kinetic energy and potential energy have been given. Therefore, we need to derive a formula for the dissipated energy and supplied energy. Following , we set as the dissipative power, which is the changing rate of dissipated energy. From (9), we haveTerms , and are the damping power of the gyrostat system. Therefore, is really the total dissipative power. Likewise, we defineas the external power. Then the Kolmogorov model (17) is rewritten as
Remark 2. Differing from the original Kolmogorov model (5), in which only inertial torque and internal torque link with the kinetic and potential energies, the dissipative torque and external torque are expressed using the corresponding energies in (23).
We already have the external torque and dissipative torque . Settingas the error torque, as given in , we define and expandwhich gives the error power as the difference between the external power and total dissipative power.
4.2. Energy Transformation
From (6), we haveNext, we haveFrom (13), (23), (24), and (27), we then haveTherefore, the rate of change of kinetic energy is the sum of the exchange between the inside conservative part (kinetic and potential) and the exchange between kinetic and error energy. Likewise, we haveHence, the rate of change of potential energy is the negative exchange power between kinetic energy and potential energy and the exchange power between potential energy and error energy. For the Hamiltonian energy, we haveTherefore, the rate of change in the Hamiltonian energy is determined by the exchange power between Hamiltonian energy and error energy.
4.3. Casimir Function and Error Power
The Casimir function , like enstrophy or the potential vorticity in a fluid dynamical context, is very useful in analyzing stability conditions and the global description of a dynamical system . Given that the Casimir energy function is , then and, for any , we haveTherefore, is a constant of the motion of the Hamiltonian system. Analogously, define the Casimir energy function for the gyrostat chaotic system as
Remark 3. For a mechanical gyrostat system, the physical meanings of terms , and in (32) are stored mechanical rotational energies. Therefore, the Casimir energy is the total stored energy for the gyrostat system.
Calling the rate of change of the Casimir energy the Casimir power, then from (6), (24), and (31), the Casimir power isFrom (25), we obtainFrom (32) and (34), we have
Remark 4. The Casimir energy is the stored mechanical energy, but the Casimir power is the error power between the external and the dissipative.
Remark 5. The Casimir power sign, sign switch, and the extension of the power switch are important indicators of the definitive dynamical behavior. The leap in Casimir energy level is also an important indicator. Both the switch in Casimir power sign and the Casimir energy leap influence the dynamical bifurcation. We conclude the following:
) If holds for all times, the supplied power is greater than the dissipative power, and the Casimir energy keeps increasing, and the system diverges as a source.
) If holds for all times, then the orbit of the system converges to a sink.
) A nonzero constant Casimir function, i.e., with , signifies a conservative system that produces a periodic orbit.
) If converges and finally leaps periodically along the zero line, then the orbit is a limit cycle.
) If is bounded and irregularly vibrating around the zero line, then the orbit of the system is either chaotic or aperiodic.
) If converges into a zero asymptotically, the orbit converges to a sink.
) The bifurcation of Casimir power is identical with the bifurcation of the system’s states or the spectrum of the Lypapunov exponent. However, the former is analytic, which is helpful in revealing the mechanism of production of different dynamical behaviors.
The extremal points of the Casimir energy for the gyrostat chaotic system (4) are obtained by letting . They fall on a triaxle hyperboloid,
5. Dynamical Behaviors Influenced by Force Field and Energy
We now investigate the dynamic behavior corresponding to the different forms of torques in the gyrostat system (4) to find out the effects of each type of torque and energy on the system. We take , , , , , , , and and set the external torques , , and , initial conditions , and sampling time .
5.1. Gyrostat System Governed by Inertia and Internal Force Field
System (17) contains the inertial torque and internal torque, i.e., it has kinetic energy and potential energy. Correspondingly, the system is written aswhich means system (4) becomesThe rotation of the rotor relative to the carrier frame produces angular momenta and corresponding potential energy, both of which have an effect on the dynamics of the system. The principal moments of inertia of the gyrostat system, , determine the kinetic energy. A periodic orbit is generated [Figure 2(a)]. Why is the periodic orbit produced, and how do the two forms of energies influence the behavior? There is a transfer between kinetic energy and potential energy. From (28) and (29), we obtainTherefore, the amount of kinetic energy loss equals the potential energy gained, meaning the system is also conservative withThe exchange of kinetic energy and potential energy and the conservation of the total ensures the system behaves periodically. The periodic orbit is symmetrical along the red dotted line [Figure 2(a)]. The color coding indicates the kinetic energy value; a change from dark blue to yellow means the kinetic energy increases, . When the orbit reaches the red circle, the kinetic energy reaches its maximum. With , the kinetic energy decreases until reaching a minimum (black stars) upon which the cycle repeats. Figure 2(b) reflects the two forms of periodic energy cycling; also, the amount of kinetic energy loss equals the potential energy gained. The Hamiltonian is constant from (39) with . The Casimir energy is also constant with , which indicates that the system does not exchange Casimir energy with the outside.
Similarly, when the system contains only kinetic energy or potential energy, the system still undergoes a periodic orbit with energy conserved, but the two forms of orbits are simpler than the orbit produced by the sum of the two energies.
5.2. Gyrostat System Governed by Inertia, Internal, and Dissipative Torques
This governing equation of this gyrostat chaotic system is written asThe system contains kinetic, potential, and dissipative energies. System (4) becomesThe divergence of the system is where is the volume in phase space of the system, with solution . With the given parameter values, we obtain , indicating the system is dissipative, so the volume index of the phase space decreases to . In addition, the system is fractal for Lyapunov dimension .
There are two forms of dissipation: Rayleigh and non-Rayleigh . In physics, it frequently happens that the frictional (damping) force is proportional to the velocity, so that its component has the form . When all the friction coefficients, are positive, the dissipation is Rayleigh, otherwise, non-Rayleigh. Therefore, system (42) is a non-Rayleigh dissipative system with ; likewise, the original gyrostat system (4) is non-Rayleigh. A Rayleigh dissipative system without an external force never produces chaos for it evolves towards a sink . System (42) is non-Rayleigh dissipative, with term playing the role of a negative damping or external torque, so we cannot determine whether it produces a sink even despite having a negative divergence. However, we find the mechanism from the Casimir power, which isFrom (44), the Casimir power has an uncertainty in terms of sign; probably, the Casimir power varies and oscillates with the system orbiting. This is confirmed by Figure 3(a). The oscillations of Casimir power reflect that the system has a sudden energy exchange between input and output. Although system (42) does not contain the external torque, it contains a negative damping term , which plays the role of an external torque. Actually, the system generates a four-wing chaotic attractor [Figure 3(b)] with Lyapunov exponents , and fractal Lyapunov dimension . Although the phase-space volume of the system converges to zero, the energy does not.
From (4), the force field is the vector field, which is also the reaction torque , and because , , and , the force field vector is . We can plot the force field to show how the force impacts the orbit of system. Figure 3(c) presents the reaction torque in plane, which is the force field flow of the system. Figure 3(d) shows a trajectory phase diagram along with flowing Casimir energy in which the stored energy changes from a minimum to a maximum and continually cycles. Compared with Figure 3(d), the magnitude of the torque in the torque field determines the amount of energy stored in the orbit, and the direction of the orbit is also controlled by the direction of the torque [Figure 3(c)].
5.3. Gyrostat System Governed by All Four Force Fields
When the gyrostat system contains all four torques—inertial, internal, dissipative (disturbance torques), and external—the governing equation iswhich is fully equivalent to system (4). The Lyapunov exponents are found to be , and the Lyapunov dimension , when the external torque is set to . Comparing system (42) with system (4), we find that the positive Lyapunov exponent is larger than for a system without external torque. Additionally, the Lyapunov dimension is slightly bigger compared with the value of the system (42). Therefore, the existence of external torque complicates the dynamics and chaos of the gyrostat system [Figure 4(a)]. Note that external torque does not always make the system more chaotic or disordered, but it does add to the complexity of the system.
Figure 4(a) shows that the value of the Casimir energy (stored energy) varies with the orbit of the system; that is, the stored energy changes from a minimum to a maximum value and circulates constantly. From (33) and (34), the Casimir power (error power) isIn the middle panel of Figure 4(b), we obtain an error power between the dissipation and external energy that shifts in the range . The energy stored in the system is released and absorbed rapidly, causing the orbit of the system to become chaotic and folded [see top and bottom panels in Figure 4(b)].
The Casimir power is given by (33), and the intersected hyperboloid with chaotic attractor is given by in (36). The intersected hyperboloid is rewritten asIt is designed to analyze the changes in Casimir energy in the system trajectory [Figure 4(c)], from which we may find the maxima and minima of the Casimir energy (32). The orbit of the chaotic attractor of system (4) intersects the hyperboloid with a change in Casimir energy [Figure 4(c)]. When , i.e., the orbits are located inside the hyperboloid, the Casimir energy keeps increasing until it reaches a maximum (red circle). Conversely, when , i.e., the orbits are located outside the hyperboloid, the Casimir energy starts to decrease, eventually to a minimum (yellow stars), where the cycle then repeats. Figure 4(d) presents a projection of the chaotic attractor onto the plane with the Casimir energy (color coded) and signs: maxima (red circles) and minima (yellow stars). The analysis above demonstrates that the Casimir energy and its associated power are good indicators in understanding the folding and expanding of the orbit of system.
The equilibrium point of system (4) and its corresponding eigenvalues are presented in Table 1. Point is the saddle node, and the remaining equilibrium points are saddle focus nodes. The five equilibrium points are on the surface of [Figure 4(d)], with point marked by a green triangle, and , and marked by the four green pentagonal stars on the subcentral position of the four-wing chaotic attractor. We find that the four saddle focus nodes [Figure 4(d)] are also connection points of the maxima (red circles) and minima (yellow stars) of stored energy.
6. Bifurcation and Dynamics Analysis of a Gyrostat System
In the dynamic analysis of the gyrostat system, changes in system parameters affect the dynamics of the system. Because some parameters are determined by their structure, these system parameters are fixed at .
6.1. Dynamics Analysis of Gyrostat System by Varying Parameter
Our purpose is to study the effects of the dissipative damping parameters on the dynamic behavior of the system in the gyrostat system with full torques, i.e., system (4). In this section, we demonstrate the similarity between the bifurcation diagram of the state and the bifurcation of the Casimir power highlighting that the Casimir power is a good indicator in analyzing the dynamics of the system.
Figure 5(a) demonstrates the state bifurcation diagram of the gyrostat system with respect to damping parameter . The system undergoes rich bifurcation dynamics with sinks, periodic orbits, multiperiodic orbits, and chaos, and then the reverse process. When the gyrostat system operates in the chaotic state, the angular velocity oscillates tremendously and, suddenly, there occurs a motion which would probably disrupt the system in real situations. Therefore, the dissipative term or disturbance affects the system hugely. Other damping factors and other parameters most assuredly influence system dynamics but are not studied in the paper. Figure 5(a1) amplifies the state bifurcation in the range , in which the bifurcation process is clearer. The Lyapunov exponent spectrum with respect to the same parameter reflecting the bifurcation process is exactly the same as the state bifurcation diagram [Figure 5(a)]; Figure 5(b1) is contrasted with Figure 5(a1).
What makes the bifurcation take place? We are unable to find a clue from the state bifurcation diagram and not from the Lyapunov exponent spectrum. However, we may find a clue from the Casimir power, (33) and (34), which are rewritten asThe Casimir power reflects the error power between the supplied and the dissipated power. When the damping parameter is varied, the sign of the Casimir power changes; therefore, the pattern by which the system absorbs and expends energy always switches, resulting thereby in the bifurcation behavior. Similarly, varying other parameters impacts the bifurcation of the gyrostat system.
The bifurcations of the Casimir power with respect to parameter are shown in Figures 5(c) and 5(c1). As parameter varies, the sign of Casimir power switches. The positive (negative) implies the supplied power is greater (less) than the dissipative power, so the Casimir energy increases (decreases) [Figures 5(d) and 5(d1)]. When the Casimir energy level leaps from one level to another, a bifurcation of the definitive behavior of the system takes place. When the Casimir power converges to zero in the range , the system orbit converges to a sink [Figure 5(c1)], which is identical to the range of the sink in Figures 5(a1), 5(b1), and 5(d1) in accord with Remark 5(. When the Casimir power bifurcates from zero to switching periodically along the zero line in the range , the Casimir energy leaps when at which the system’s definitive behavior bifurcates to a periodic orbit in accord with Remark 5(. As the switching amplitude of the Casimir power reaches some threshold, i.e., the Casimir energy leaps to a new level at , the system behavior loses its stable periodicity and transits to the next definitive behavior, i.e., periodic-doubling. At point , the Casimir energy leaps to a new level, the system orbit loses stable periodic multidoubling behavior and transits to an irregular oscillation, i.e., chaos or a chaotic attractor is constructed in accordance with Remark 5(. The violent absorption and release of energy stored in the system causes the generation of chaos.
A chaotic gyrostat model was established that provides a reference model for researchers engaged in dynamical and gyrostat studies. The gyrostat was transformed into a Kolmogorov model to separate the torques. The conversion of energy and the cycles performed were revealed through identifying four forms of energy—kinetic, potential, dissipative, and supplied. The Casimir power is equivalent to the error power between the supplied power and the dissipative power. The Casimir power and energy leap were successfully used as indicators that assisted in evaluating the different dynamical features, such as sinks, periodic orbits, and chaos, and were identical to the indicators of the Lyapunov exponent spectrum and state bifurcation. However, the Casimir power is more convenient for the analysis of the mechanism in understanding different types of dynamics generation because of the analytical form. The idea and method may be applicable to the analyses of mechanisms for other dynamical systems. Therefore, how to find some direction information from the form of Casimir power prior to the numerical simulation is a significant work. Surely, there are many further investigations, such as Hopf bifurcation, pitchfork bifurcation, and visual demonstration, and hidden attractor can be conducted, which has provided some open topics for the other academic researchers.
The data used to support the findings of this study can be computed by the relative equations and parameters given in the article using ODE5 solver in Matlab and are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This work is supported by the National Natural Science Foundation of China (61873186) and Natural Science Foundation of Tianjin City (17JCZDJC38300). We thank Richard Haase, Ph.D., from Liwen Bianji, Edanz Group China (www.liwenbianji.cn/ac), for editing the English text of a draft of this manuscript.
J. Wittenburg, Dynamics of Systems of Rigid Bodies, Teubner, 1977.View at: MathSciNet
Y. Z. Liu and F. P. J. Rimrott, “Global motion of a dissipative asymmetric gyrostat,” Archive of Applied Mechanics, vol. 62, no. 5, pp. 329–337, 1992.View at: Google Scholar
W. Zhang, J. Zhang, and M. Yao, “The extended Melnikov method for non-autonomous nonlinear dynamical systems and application to multi-pulse chaotic dynamics of a buckled thin plate,” Nonlinear Analysis: Real World Applications, vol. 11, no. 3, pp. 1442–1457, 2010.View at: Publisher Site | Google Scholar
J. Zhang and W. Zhang, “An extended high-dimensional Melnikov analysis for global and chaotic dynamics of a non-autonomous rectangular buckled thin plate,” Science China Physics, Mechanics and Astronomy, vol. 55, no. 9, pp. 1679–1690, 2012.View at: Google Scholar