# Modeling of a Chaotic Gyrostat System and Mechanism Analysis of Dynamics Using Force and Energy

Authors: Guoyuan Qi; Xiaogang Yang
Journal: Complexity
 (2019)

---

## Abstract

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.


---

## Body

## 1. Introduction

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 [1]. 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 [3], Krishnaprasad [4] 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 [6]. 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 [17] 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 [18] 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 [19] 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 [31], 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. Sections2–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 Figure1). Euler’s dynamical equations of rotational motion for a disturbed gyrostat are [1](1)Ixdωxdt=Iy-Izωyωz-ωyhz+ωzhy-Tx+Lx,Iydωydt=Iz-Ixωzωx-ωzhx+ωxhz-Ty+Lx,Izdωzdt=Ix-Iyωxωy-ωxhy+ωyhx-Tz+Lx,where ω=ωxωyωzT is the angular velocity vector; Ix,Iy, and Iz 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; hx,hy, and hz are the constant components of the total angular momentum vector of the reaction rotors with respect to the body-fixed coordinate system; Tx,Ty, and Tz are components of the vector of disturbance torques; and Lx, Ly, and Lz are external torques.Figure 1
Gyrostat system model.


### 2.2. Disturbance Torques

The expressions specifying disturbance torques that exist in the gyrostat are ([9])(2)Tx=εμxxωx+μxyωy+μxzωz,Ty=εμyxωx+μyyωy+μyzωz,Tz=εμzxωx+μzyωy+μzzωz,where μij(i,j=x,y,z) 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, μxy=μxz=μyx=μyz=μzy=μzx=0; moreover, we set μxx=μx,μyy=-μy,μzz=μz, and ε=1. μij(i,j=x,y,z) is the damping coefficient with μij>0 meaning that the energy is being dissipated from the system and μij<0 meaning that the energy is supplied to the system [33]. Here, μyy=-μy means that the damping is negative, which is equivalent to a supplied force of the angular velocity ωy along the y-axis direction. In practice, to implement the negative damping, we can install a motor for rotor 2, and it rotates the axis y with supplied torque -μyωy.The disturbance torque is expressed as follows:(3)TxTyTz=μx000-μy000μzωxωyωz.The focus here is to study the attitude motion of the gyrostat system for the particular case hx=0,hy≠0,hz≠0,Ix≠Iy≠Iy as studied in [12]. Straightforwardly, as given in [12], we obtain from (1)-(3) (4)Ixdωxdt=Iy-Izωyωz-ωyhz+ωzhy-μxωx+Lx,Iydωydt=Iz-Ixωzωx+ωxhz+μyωy+Ly,Izdωzdt=Ix-Iyωxωy-ωxhy-μzωz+Lz.


## 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 Figure1). Euler’s dynamical equations of rotational motion for a disturbed gyrostat are [1](1)Ixdωxdt=Iy-Izωyωz-ωyhz+ωzhy-Tx+Lx,Iydωydt=Iz-Ixωzωx-ωzhx+ωxhz-Ty+Lx,Izdωzdt=Ix-Iyωxωy-ωxhy+ωyhx-Tz+Lx,where ω=ωxωyωzT is the angular velocity vector; Ix,Iy, and Iz 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; hx,hy, and hz are the constant components of the total angular momentum vector of the reaction rotors with respect to the body-fixed coordinate system; Tx,Ty, and Tz are components of the vector of disturbance torques; and Lx, Ly, and Lz are external torques.Figure 1
Gyrostat system model.


## 2.2. Disturbance Torques

The expressions specifying disturbance torques that exist in the gyrostat are ([9])(2)Tx=εμxxωx+μxyωy+μxzωz,Ty=εμyxωx+μyyωy+μyzωz,Tz=εμzxωx+μzyωy+μzzωz,where μij(i,j=x,y,z) 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, μxy=μxz=μyx=μyz=μzy=μzx=0; moreover, we set μxx=μx,μyy=-μy,μzz=μz, and ε=1. μij(i,j=x,y,z) is the damping coefficient with μij>0 meaning that the energy is being dissipated from the system and μij<0 meaning that the energy is supplied to the system [33]. Here, μyy=-μy means that the damping is negative, which is equivalent to a supplied force of the angular velocity ωy along the y-axis direction. In practice, to implement the negative damping, we can install a motor for rotor 2, and it rotates the axis y with supplied torque -μyωy.The disturbance torque is expressed as follows:(3)TxTyTz=μx000-μy000μzωxωyωz.The focus here is to study the attitude motion of the gyrostat system for the particular case hx=0,hy≠0,hz≠0,Ix≠Iy≠Iy as studied in [12]. Straightforwardly, as given in [12], we obtain from (1)-(3) (4)Ixdωxdt=Iy-Izωyωz-ωyhz+ωzhy-μxωx+Lx,Iydωydt=Iz-Ixωzωx+ωxhz+μyωy+Ly,Izdωzdt=Ix-Iyωxωy-ωxhy-μzωz+Lz.


## 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 [17] presented a Kolmogorov model describing dissipative-forced dynamical systems given in the form(5)x˙=x,H-Λx+f,where x=x1x2x3T, ·,· represents the Lie–Poisson structure of the kinetic energy part of the Hamiltonian function H of a system, Λ is a positive definite diagonal matrix, -Λx represents the dissipation (or disturbance torques), and f represents the external torque. For a three-dimensional system, the Lie–Poisson bracket on functions of x takes the form [34](6)F,Gx=-x,∇F×∇G,where F, G∈C∞(g∗), g is a Lie algebra and , is the inner product. Hence, from (6) we have(7)x,H=-x·∇x×∇H=x×∇H.In the absence of external forces, the Euler equations for the rotational dynamics of a rigid body about its center of mass are [34](8)L˙1=I2-I3ω2ω3=Π3-Π2L2L3,L˙2=I3-I1ω1ω3=Π1-Π3L1L3,L˙3=I1-I2ω1ω2=Π2-Π1L1L2,where Πi=Ii-1i=1,2,3, Ii is the principle moment of inertia for the group SO(3), and Li is the angular momentum satisfying Li=Iiωi, with ωi being the angular velocity corresponding Ii.To satisfy the skew-symmetric property ofx,H for system (4), we setω-=ω-xω-yω-zT, ω=ωxωyωzT, and M=diagIxIyIz. We introduce the reparameterization(9)ω-=Mωinto the original gyrostat system (4) to obtain(10)dω-xdt=Iy-IzIxIyIzω-yω-z-hzIxIyω-y+hyIxIzω-z-μxIxω-x+LxIx,dω-ydt=Iz-IxIxIyIzω-xω-z+hzIxIyω-x+μyIyω-y+LyIy,dω-zdt=Ix-IyIxIyIzω-xω-y-hyIxIzω-x-μzIzω-z+LzIz.Setting(11)Π3-Π2=Iy-IzIxIzIy,Π1-Π3=Iz-IxIxIzIy,Π2-Π1=Ix-IyIxIzIy,we then find(12)Π1Π2Π3T=11+Ix-IyIxIzIy1+Ix-IzIxIzIyT.Using (9), the kinetic energy is(13)K=12Π1ω-x2+Π2ω-y2+Π3ω-z2=12Π1Ixωx2+Π2Iyωy2+Π3Izωz2.Setting the internal forcing as(14)β=0-hyIxIz-hzIxIyT,then the potential energy is(15)U=ω-,β=-hzIxIyω-z-hyIxIzω-y=-hzIzIxIyωz-hyIyIxIzωy.The Hamiltonian function is(16)H=K+U.As a Kolmogorov model, (10) becomes(17)
ω
-
˙
=
ω
-
,
K
+
ω
-
,
U
-
Λ
ω
-
+
f
=
ω
-
×
∇
K
+
ω
-
×
∇
U
-
Λ
ω
-
+
f
=
ω
-
×
Π
ω
-
+
ω
-
×
β
-
Λ
ω
-
+
f
=
I
y
-
I
z
I
x
I
y
I
z
ω
-
y
ω
-
z
I
z
-
I
x
I
x
I
y
I
z
ω
-
x
ω
-
z
I
x
-
I
y
I
x
I
y
I
z
ω
-
x
ω
-
y
+
h
y
I
x
I
z
ω
-
z
-
h
z
I
x
I
y
ω
-
y
h
z
I
x
I
y
ω
-
x
-
h
y
I
x
I
z
ω
-
x
-
μ
x
I
x
ω
-
x
-
μ
y
I
y
y
ω
-
y
μ
z
I
z
ω
-
z
+
L
x
I
x
L
y
I
y
L
z
I
z
,where(18)Λ=diagμxIx-μyIyμzIz,f=LxIxLyIyLzIzT.To study system (4) further, we rearrange (17). Substituting (9) into (17), we have(19)
I
x
ω
˙
x
I
y
ω
˙
y
I
z
ω
˙
z
=
I
y
-
I
z
I
x
ω
y
ω
z
I
z
-
I
x
I
y
ω
x
ω
z
I
x
-
I
y
I
z
ω
x
ω
y
+
h
y
I
x
ω
z
-
h
z
I
x
ω
y
h
z
I
y
ω
x
-
h
y
I
z
ω
x
-
μ
x
I
x
ω
x
-
μ
y
I
y
ω
y
μ
z
I
z
ω
z
+
L
x
I
x
L
y
I
y
L
z
I
z
.Equation (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 [24], we set (20)D=12ω-,Λω-,as the dissipative power, which is the changing rate of dissipated energy. From (9), we have(21)D=12ω-,Λω-=12μxωx2-μyωy2+μzωz2.Terms μxωx2,-μyωy2, and μzωz2 are the damping power of the gyrostat system. Therefore, 2D=μxωx2-μyωy2+μzωz2 is really the total dissipative power. Likewise, we define(22)G=ω-,f=Lxωx+Lyωy+Lzωz,as the external power. Then the Kolmogorov model (17) is rewritten as(23)ω-˙=ω-,H-∇D+∇G.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 torquef and dissipative torque Λω-. Setting(24)e=f-Λω-,as the error torque, as given in [25], we define and expand(25)EP=ω-,e=ω-,f-ω-,Λω-=G-2D,which gives the error power as the difference between the external power and total dissipative power.


### 4.2. Energy Transformation

From (6), we have(26)∇K,ω-,H=∇K,ω-×∇H=-ω-,∇K×∇H=K,H.Next, we have(27)∂K∂ω-ω-,H=K,H,∂U∂ω-ω-,H=U,H.From (13), (23), (24), and (27), we then have(28)K˙=∂K∂ω-ω-˙=∂K∂ω-ω-,H+∇K,e=K,K+K,U+Πω-,e=hyωzωxΠ1-Π3+hzωyωxΠ2-Π1-Π1μxωx2-Π2μyωy2+Π3μzωz2+Π1Lxωx+Π2Lyωy+Π3Lzωz.Therefore, 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 have(29)U˙=∂U∂ω-ω-˙=∂U∂ω-ω-,H+∇U,e=-K,U+β,e=hyωzωxΠ3-Π1+hzωyωxΠ1-Π2-hyμyIxIyIzωy+hzμzIxIyIzωz-hyLyIxIyIz-hzLzIxIyIz.Hence, 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 have(30)H˙=∇H,ω-˙=H,H+∇H,e=Πω-,e+β,e.Therefore, 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 functionC, 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 [34]. Given that the Casimir energy function is C(x)=1/2x2, then ∇C(x)=x and, for any F∈f(R3), we have(31)C,Fx=-x·∇C×∇F=-x·x×∇F=-∇F·x×x=0.Therefore, C is a constant of the motion of the Hamiltonian system. Analogously, define the Casimir energy function for the gyrostat chaotic system as(32)C=12ω-,ω-=12Ixωx2+Iyωy2+Izωz2.Remark 3.
For a mechanical gyrostat system, the physical meanings of termsIxωx2,Iyωy2, and Izωz2 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 is(33)C˙=∂C∂ω-ω-,H+∇C,e=C,H+∇C,e=∇C,e=ω-,e=Lxωx+Lyωy+Lzωz-μxωx2-μyωy2+μzωz2.From (25), we obtain(34)C˙=EP.From (32) and (34), we have(35)Ct=∫0tEPτdτ+C0.Remark 4.
The Casimir energy is the stored mechanical energy, but the Casimir power is the error power between the external and the dissipative.Based on the conclusion concerning Casimir power [24, 25, 31], we summarize as follows.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:
( 1) If C˙(x)>0 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.
( 2) If C˙(x)<0 holds for all times, then the orbit of the system converges to a sink.
( 3) A nonzero constant Casimir function, i.e., C˙≡0 with C≠0, signifies a conservative system that produces a periodic orbit.
( 4) If C˙(x) converges and finally leaps periodically along the zero line, then the orbit is a limit cycle.
( 5) If C˙(x) is bounded and irregularly vibrating around the zero line, then the orbit of the system is either chaotic or aperiodic.
( 6) If C˙(x) converges into a zero asymptotically, the orbit converges to a sink.
( 7) 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 C˙=0. They fall on a triaxle hyperboloid,(36)Ξ0:μxωx-Lx2μx2-μyωy+Ly2μy2+μzωz-Lz2μz2=Lx24μx-Ly24μy+Lz24μz.


## 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 [24], we set (20)D=12ω-,Λω-,as the dissipative power, which is the changing rate of dissipated energy. From (9), we have(21)D=12ω-,Λω-=12μxωx2-μyωy2+μzωz2.Terms μxωx2,-μyωy2, and μzωz2 are the damping power of the gyrostat system. Therefore, 2D=μxωx2-μyωy2+μzωz2 is really the total dissipative power. Likewise, we define(22)G=ω-,f=Lxωx+Lyωy+Lzωz,as the external power. Then the Kolmogorov model (17) is rewritten as(23)ω-˙=ω-,H-∇D+∇G.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 torquef and dissipative torque Λω-. Setting(24)e=f-Λω-,as the error torque, as given in [25], we define and expand(25)EP=ω-,e=ω-,f-ω-,Λω-=G-2D,which gives the error power as the difference between the external power and total dissipative power.


## 4.2. Energy Transformation

From (6), we have(26)∇K,ω-,H=∇K,ω-×∇H=-ω-,∇K×∇H=K,H.Next, we have(27)∂K∂ω-ω-,H=K,H,∂U∂ω-ω-,H=U,H.From (13), (23), (24), and (27), we then have(28)K˙=∂K∂ω-ω-˙=∂K∂ω-ω-,H+∇K,e=K,K+K,U+Πω-,e=hyωzωxΠ1-Π3+hzωyωxΠ2-Π1-Π1μxωx2-Π2μyωy2+Π3μzωz2+Π1Lxωx+Π2Lyωy+Π3Lzωz.Therefore, 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 have(29)U˙=∂U∂ω-ω-˙=∂U∂ω-ω-,H+∇U,e=-K,U+β,e=hyωzωxΠ3-Π1+hzωyωxΠ1-Π2-hyμyIxIyIzωy+hzμzIxIyIzωz-hyLyIxIyIz-hzLzIxIyIz.Hence, 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 have(30)H˙=∇H,ω-˙=H,H+∇H,e=Πω-,e+β,e.Therefore, 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 functionC, 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 [34]. Given that the Casimir energy function is C(x)=1/2x2, then ∇C(x)=x and, for any F∈f(R3), we have(31)C,Fx=-x·∇C×∇F=-x·x×∇F=-∇F·x×x=0.Therefore, C is a constant of the motion of the Hamiltonian system. Analogously, define the Casimir energy function for the gyrostat chaotic system as(32)C=12ω-,ω-=12Ixωx2+Iyωy2+Izωz2.Remark 3.
For a mechanical gyrostat system, the physical meanings of termsIxωx2,Iyωy2, and Izωz2 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 is(33)C˙=∂C∂ω-ω-,H+∇C,e=C,H+∇C,e=∇C,e=ω-,e=Lxωx+Lyωy+Lzωz-μxωx2-μyωy2+μzωz2.From (25), we obtain(34)C˙=EP.From (32) and (34), we have(35)Ct=∫0tEPτdτ+C0.Remark 4.
The Casimir energy is the stored mechanical energy, but the Casimir power is the error power between the external and the dissipative.Based on the conclusion concerning Casimir power [24, 25, 31], we summarize as follows.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:
( 1) If C˙(x)>0 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.
( 2) If C˙(x)<0 holds for all times, then the orbit of the system converges to a sink.
( 3) A nonzero constant Casimir function, i.e., C˙≡0 with C≠0, signifies a conservative system that produces a periodic orbit.
( 4) If C˙(x) converges and finally leaps periodically along the zero line, then the orbit is a limit cycle.
( 5) If C˙(x) is bounded and irregularly vibrating around the zero line, then the orbit of the system is either chaotic or aperiodic.
( 6) If C˙(x) converges into a zero asymptotically, the orbit converges to a sink.
( 7) 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 C˙=0. They fall on a triaxle hyperboloid,(36)Ξ0:μxωx-Lx2μx2-μyωy+Ly2μy2+μzωz-Lz2μz2=Lx24μx-Ly24μy+Lz24μz.


## 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 Ix=0.85Kg.m2, Iy=0.45Kg.m2, Iz=0.2Kg.m2, hy=0.57416Kg.m2/s, hz=2.38Kg.m2/s, μx=7.6N.m.s, μy=6.42N.m.s, and μz=5.8N.m.s and set the external torques Lx=0N.m, Ly=0N.m, and Lz=22.8N.m, initial conditions ωx0=ωz0=ωy0=2.5rad/s, and sampling time τ=0.0001s.


### 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 as(37)ω-˙=ω-,K+ω-,U,which means system (4) becomes(38)Ixω˙xIyω˙yIzω˙z=Iy-IzωyωzIz-IxωxωzIx-Iyωxωy+hyωz-hzωyhzωx-hyωx.The rotation of the rotor relative to the carrier frame produces angular momenta hx,hy,hz 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, Ix,Iy,Iz, 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 obtain(39)K˙=K,U=hyωzωxΠ1-Π3+hzωyωxΠ2-Π1,U˙=-K,U.Therefore, the amount of kinetic energy loss equals the potential energy gained, meaning the system is also conservative with(40)H˙=H,H=0.The 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, K˙>0. When the orbit reaches the red circle, the kinetic energy reaches its maximum. With K˙<0, 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 H=1.5522(J). The Casimir energy is also constant with C=C0=4.6875(J), which indicates that the system does not exchange Casimir energy with the outside.Figure 2
System governed by inertial and internal torques: (a) a periodic orbit with a kinetic energy cycle, (b) kinetic, potential, and Hamiltonian energy.
(a) (b)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 as(41)ω-˙=ω-,K+ω-,U-∇D.The system contains kinetic, potential, and dissipative energies. System (4) becomes(42)Ixω˙xIyω˙yIzω˙z=Iy-IzωyωzIz-IxωxωzIx-Iyωxωy+hyωz-hzωyhzωx-hyωx-μxωx-μyωyμzωz.The divergence of the system is (43)divV=∂ω˙x∂ωx+∂ω˙y∂ωy+∂ω˙z∂ωz=-μxIx+μyIy-μzIz,where V is the volume in phase space of the system, with solution V=V0exp-μx/Ix-μy/Iy+μz/Iz. With the given parameter values, we obtain -μx/Ix+μy/Iy-μz/Iz=-23.6745, indicating the system is dissipative, so the volume index of the phase space decreases to 0. In addition, the system is fractal for Lyapunov dimension Ld<3.There are two forms of dissipation: Rayleigh and non-Rayleigh [23]. In physics, it frequently happens that the frictional (damping) force is proportional to the velocity, so that its x component has the form Fxi=-kxiνxi. When all the friction coefficients, kxi are positive, the dissipation is Rayleigh, otherwise, non-Rayleigh. Therefore, system (42) is a non-Rayleigh dissipative system with ∂ω˙y/∂ωy=μx/Ix>0; 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 [20]. System (42) is non-Rayleigh dissipative, with term μyωy 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 is(44)C˙=-μxωx2+μyωy2-μzωz2.From (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 μyωy, which plays the role of an external torque. Actually, the system generates a four-wing chaotic attractor [Figure 3(b)] with Lyapunov exponents L1=3.8719,L2=0,L3=-27.5464, and fractal Lyapunov dimension Ld=2.1406<3. Although the phase-space volume of the system converges to zero, the energy does not.Figure 3
System governed by inertia, internal, and dissipative torques: (a) Casimir energy and Casimir power, (b) four-wing chaotic orbit, (c) phase plane of force field, and (d) phase plane orbit along with flowing Casimir energy.
(a) (b) (c) (d)From (4), the force field is the vector field, which is also the reaction torque Ixω˙x,Iyω˙y,Izω˙zT, and because Lx=Ixωx, Ly=Iyωy, and Lz=Izωz, the force field vector is L˙x,L˙y,L˙zT=Ixω˙x,Iyω˙y,Izω˙zT. We can plot the force field to show how the force impacts the orbit of system. Figure 3(c) presents the reaction torque in L˙x-L˙y 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 is(45)ω-˙=ω-,H-∇D+∇G,which is fully equivalent to system (4). The Lyapunov exponents are found to be L1=3.8975,L2=0,L3=-27.572, and the Lyapunov dimension Ld=2.1411<3, when the external torque is set to Lx=0,Ly=0,Lz=22.18. Comparing system (42) with system (4), we find that the positive Lyapunov exponent L1=3.8975 is larger than L1=3.8719 for a system without external torque. Additionally, the Lyapunov dimension Ld=2.1411 is slightly bigger compared with the value Ld=2.1406 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
System under full torques. (a) 3D orbits with flowing stored energy, (b) time series of the stored energy, error power, and angular velocityωx, (c) intersection of the chaotic system trajectory and hyperboloid of the Casimir function, the maxima and minima of which are marked by red circles and yellow stars, respectively, and (d) signs of extrema and equilibrium points on the four-wing chaotic attractor.
(a) (b) (c) (d)Figure4(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) is(46)C˙=EP=Lxωx+Lyωy+Lzωz-μxωx2-μyωy2+μzωz2.In the middle panel of Figure 4(b), we obtain an error power between the dissipation and external energy that shifts in the range -1.974×10-5,1.185×105. 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 Ξ0 in (36). The intersected hyperboloid Ξ0 is rewritten as(47)Ξ0:μxωx-Lx2μx2-μyωy-Ly2μy2+μzωz-Lz2μz2=Lx24μx-Ly24μy+Lz24μz.It 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 C˙>0, i.e., the orbits are located inside the hyperboloid, the Casimir energy keeps increasing until it reaches a maximum (red circle). Conversely, when C˙<0, 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 ωx-ωy 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 S3 is the saddle node, and the remaining equilibrium points are saddle focus nodes. The five equilibrium points are on the surface of Ξ0 [Figure 4(d)], with point S3 marked by a green triangle, and S1,S2,S4, and S5 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.Table 1
Equilibrium points and corresponding eigenvalues of system (4).
Equilibrium point Eigenvalue(λ1,λ2,λ3) S 1   ( 11.6173 , - 17.9707 , - 11.6170 ) - 33.2194 4.7725 + j 18.9858 4.7725 - j 18.9858 S 2   ( - 12.5942,17.9807 , - 10.4397 ) - 32.9280 4.6267 + j 19.5380 4.6267 - j 19.5380 S 3   ( 0.2935 , 0.0071,3.9021 ) - 9.0583 - 28.8997 14.2835 S 4   ( 12.2571,23.8915,22.9136 ) - 37.6678 6.9966 + j 19.8438 6.9966 - j 19.8438 S 5   ( - 11.5737 ,  - 25.6311,25.5351 ) - 38.9770 7.6512 + j 20.4982 7.6512 - j 20.4982


## 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 as(37)ω-˙=ω-,K+ω-,U,which means system (4) becomes(38)Ixω˙xIyω˙yIzω˙z=Iy-IzωyωzIz-IxωxωzIx-Iyωxωy+hyωz-hzωyhzωx-hyωx.The rotation of the rotor relative to the carrier frame produces angular momenta hx,hy,hz 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, Ix,Iy,Iz, 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 obtain(39)K˙=K,U=hyωzωxΠ1-Π3+hzωyωxΠ2-Π1,U˙=-K,U.Therefore, the amount of kinetic energy loss equals the potential energy gained, meaning the system is also conservative with(40)H˙=H,H=0.The 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, K˙>0. When the orbit reaches the red circle, the kinetic energy reaches its maximum. With K˙<0, 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 H=1.5522(J). The Casimir energy is also constant with C=C0=4.6875(J), which indicates that the system does not exchange Casimir energy with the outside.Figure 2
System governed by inertial and internal torques: (a) a periodic orbit with a kinetic energy cycle, (b) kinetic, potential, and Hamiltonian energy.
(a) (b)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 as(41)ω-˙=ω-,K+ω-,U-∇D.The system contains kinetic, potential, and dissipative energies. System (4) becomes(42)Ixω˙xIyω˙yIzω˙z=Iy-IzωyωzIz-IxωxωzIx-Iyωxωy+hyωz-hzωyhzωx-hyωx-μxωx-μyωyμzωz.The divergence of the system is (43)divV=∂ω˙x∂ωx+∂ω˙y∂ωy+∂ω˙z∂ωz=-μxIx+μyIy-μzIz,where V is the volume in phase space of the system, with solution V=V0exp-μx/Ix-μy/Iy+μz/Iz. With the given parameter values, we obtain -μx/Ix+μy/Iy-μz/Iz=-23.6745, indicating the system is dissipative, so the volume index of the phase space decreases to 0. In addition, the system is fractal for Lyapunov dimension Ld<3.There are two forms of dissipation: Rayleigh and non-Rayleigh [23]. In physics, it frequently happens that the frictional (damping) force is proportional to the velocity, so that its x component has the form Fxi=-kxiνxi. When all the friction coefficients, kxi are positive, the dissipation is Rayleigh, otherwise, non-Rayleigh. Therefore, system (42) is a non-Rayleigh dissipative system with ∂ω˙y/∂ωy=μx/Ix>0; 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 [20]. System (42) is non-Rayleigh dissipative, with term μyωy 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 is(44)C˙=-μxωx2+μyωy2-μzωz2.From (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 μyωy, which plays the role of an external torque. Actually, the system generates a four-wing chaotic attractor [Figure 3(b)] with Lyapunov exponents L1=3.8719,L2=0,L3=-27.5464, and fractal Lyapunov dimension Ld=2.1406<3. Although the phase-space volume of the system converges to zero, the energy does not.Figure 3
System governed by inertia, internal, and dissipative torques: (a) Casimir energy and Casimir power, (b) four-wing chaotic orbit, (c) phase plane of force field, and (d) phase plane orbit along with flowing Casimir energy.
(a) (b) (c) (d)From (4), the force field is the vector field, which is also the reaction torque Ixω˙x,Iyω˙y,Izω˙zT, and because Lx=Ixωx, Ly=Iyωy, and Lz=Izωz, the force field vector is L˙x,L˙y,L˙zT=Ixω˙x,Iyω˙y,Izω˙zT. We can plot the force field to show how the force impacts the orbit of system. Figure 3(c) presents the reaction torque in L˙x-L˙y 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 is(45)ω-˙=ω-,H-∇D+∇G,which is fully equivalent to system (4). The Lyapunov exponents are found to be L1=3.8975,L2=0,L3=-27.572, and the Lyapunov dimension Ld=2.1411<3, when the external torque is set to Lx=0,Ly=0,Lz=22.18. Comparing system (42) with system (4), we find that the positive Lyapunov exponent L1=3.8975 is larger than L1=3.8719 for a system without external torque. Additionally, the Lyapunov dimension Ld=2.1411 is slightly bigger compared with the value Ld=2.1406 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
System under full torques. (a) 3D orbits with flowing stored energy, (b) time series of the stored energy, error power, and angular velocityωx, (c) intersection of the chaotic system trajectory and hyperboloid of the Casimir function, the maxima and minima of which are marked by red circles and yellow stars, respectively, and (d) signs of extrema and equilibrium points on the four-wing chaotic attractor.
(a) (b) (c) (d)Figure4(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) is(46)C˙=EP=Lxωx+Lyωy+Lzωz-μxωx2-μyωy2+μzωz2.In the middle panel of Figure 4(b), we obtain an error power between the dissipation and external energy that shifts in the range -1.974×10-5,1.185×105. 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 Ξ0 in (36). The intersected hyperboloid Ξ0 is rewritten as(47)Ξ0:μxωx-Lx2μx2-μyωy-Ly2μy2+μzωz-Lz2μz2=Lx24μx-Ly24μy+Lz24μz.It 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 C˙>0, i.e., the orbits are located inside the hyperboloid, the Casimir energy keeps increasing until it reaches a maximum (red circle). Conversely, when C˙<0, 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 ωx-ωy 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 S3 is the saddle node, and the remaining equilibrium points are saddle focus nodes. The five equilibrium points are on the surface of Ξ0 [Figure 4(d)], with point S3 marked by a green triangle, and S1,S2,S4, and S5 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.Table 1
Equilibrium points and corresponding eigenvalues of system (4).
Equilibrium point Eigenvalue(λ1,λ2,λ3) S 1   ( 11.6173 , - 17.9707 , - 11.6170 ) - 33.2194 4.7725 + j 18.9858 4.7725 - j 18.9858 S 2   ( - 12.5942,17.9807 , - 10.4397 ) - 32.9280 4.6267 + j 19.5380 4.6267 - j 19.5380 S 3   ( 0.2935 , 0.0071,3.9021 ) - 9.0583 - 28.8997 14.2835 S 4   ( 12.2571,23.8915,22.9136 ) - 37.6678 6.9966 + j 19.8438 6.9966 - j 19.8438 S 5   ( - 11.5737 ,  - 25.6311,25.5351 ) - 38.9770 7.6512 + j 20.4982 7.6512 - j 20.4982


## 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 atIx=0.85Kg.m2,Iy=0.45Kg.m2,Iz=0.2Kg.m2,hy=0.57416Kg.m2/s,hz=2.38Kg.m2/s.


### 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.Figure5(a) demonstrates the state bifurcation diagram of the gyrostat system with respect to damping parameter μx∈[1.2,13.4]. 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 μx∈[1.2,2.5], in which the bifurcation process is clearer. The Lyapunov exponent spectrum with respect to the same parameter μx∈[1.2,13.4] 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).Figure 5
Bifurcation of the gyrostat system (4) with respect to μx with μy=6.42N.m.s,μz=5.8N.m.s,Lx=Ly=0,Lz=22.8N.m, (a) bifurcation diagram of state ωx, (b) Lyapunov exponent spectrum, (c) bifurcation of Casimir power, and (d) bifurcation of Casimir energy.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 as(48)C˙=Ep=Lxωx+Lyωy+Lzωz-μxωx2-μyωy2+μzωz2.The Casimir power reflects the error power between the supplied and the dissipated power. When the damping parameter μx 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μx are shown in Figures 5(c) and 5(c1). As parameter μx 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 μx∈[1.2,1.39], 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(6). When the Casimir power bifurcates from zero to switching periodically along the zero line in the range μx∈(1.39,2.04], the Casimir energy leaps when μx=1.39 at which the system’s definitive behavior bifurcates to a periodic orbit in accord with Remark 5(4). As the switching amplitude of the Casimir power reaches some threshold, i.e., the Casimir energy leaps to a new level at μx=2.04, the system behavior loses its stable periodicity and transits to the next definitive behavior, i.e., periodic-doubling. At point μx=2.205, 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(6). The violent absorption and release of energy stored in the system causes the generation of chaos.


## 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.Figure5(a) demonstrates the state bifurcation diagram of the gyrostat system with respect to damping parameter μx∈[1.2,13.4]. 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 μx∈[1.2,2.5], in which the bifurcation process is clearer. The Lyapunov exponent spectrum with respect to the same parameter μx∈[1.2,13.4] 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).Figure 5
Bifurcation of the gyrostat system (4) with respect to μx with μy=6.42N.m.s,μz=5.8N.m.s,Lx=Ly=0,Lz=22.8N.m, (a) bifurcation diagram of state ωx, (b) Lyapunov exponent spectrum, (c) bifurcation of Casimir power, and (d) bifurcation of Casimir energy.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 as(48)C˙=Ep=Lxωx+Lyωy+Lzωz-μxωx2-μyωy2+μzωz2.The Casimir power reflects the error power between the supplied and the dissipated power. When the damping parameter μx 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μx are shown in Figures 5(c) and 5(c1). As parameter μx 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 μx∈[1.2,1.39], 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(6). When the Casimir power bifurcates from zero to switching periodically along the zero line in the range μx∈(1.39,2.04], the Casimir energy leaps when μx=1.39 at which the system’s definitive behavior bifurcates to a periodic orbit in accord with Remark 5(4). As the switching amplitude of the Casimir power reaches some threshold, i.e., the Casimir energy leaps to a new level at μx=2.04, the system behavior loses its stable periodicity and transits to the next definitive behavior, i.e., periodic-doubling. At point μx=2.205, 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(6). The violent absorption and release of energy stored in the system causes the generation of chaos.


## 7. Conclusion

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.



---

*Source: 5439596-2019-07-02.xml*