Abstract

This article deals with the analysis of the effects of passive control on the complex dynamics of a nonlinear damping gyroscope. After modeling the gyroscope dynamics under the influence of the control force, using the harmonic balance method, the amplitudes of the harmonic oscillations are determined. Subsequently, the Routh–Hurwitz criterion is used to analyze and determine the stability domains of the oscillations. The influence of the control force parameters on the amplitude of the oscillations is studied. The control of chaotic dynamics and the coexistence of gyroscope attractors are performed through bifurcation diagrams, Lyapunov exponents, phase portraits, and time series. Numerical simulations are used to confirm the effectiveness of the control force. This article revealed that the amplitude of the harmonic oscillations, the chaotic dynamics, and the coexistence of the attractors of the rotating gyroscope are better controlled when the latter vibrates in the opposite direction to the passive control force.

1. Introduction

The gyroscope is a device based on the principle of conservation of angular motion on one, two, or three axes with respect to an inert frame of reference. Gyroscopes are instruments used in a large number of civil and military applications, such as inertial navigation, guidance, and stabilization of platforms [1, 2]. They provide a measurement of the rotation, of the reference of the device with respect to an inertia reference frame. The first gyroscopes designed were mechanical.

Due to the moving mechanical parts that compose them, these gyroscopes are bulky, expensive, and require significant maintenance. Among the solutions developed to replace them, optical gyroscopes that do not contain moving parts have been a very advantageous option. Indeed, they offer a longer life, require less maintenance, are more compact and lighter, and better withstand shocks and strong accelerations. This system reached its full maturity at the end of the seventies, and the first demonstration of its technical competitiveness was its integration in 1978, in the navigation systems of Boeing. The number of gyroscopes used in the world is now exploding as they equip a growing number of smartphones.

They are used to precisely identify the position and orientation of the device in space in aviation for the artificial horizon, the heading maintainer, the coordinator or turn indicator, camera stabilization during a capture disturbed by the movement of the waves, and the pitch of an airplane [3]. One of the goals of technology, in general, and engineering, in particular, is therefore the miniaturization of devices and components. For industrial machines, whose complexity has constantly increased, understanding and diagnosing the vibratory phenomena involved require increasingly detailed simulations of their behavior [4]. In this context, the prediction of damping characteristics and their effects is fundamental for the design of rotating machines in order to provide an accurate idea of safe rotation ranges [5, 6]. For this reason, many studies in the field of dynamics have focused on the modeling of dissipative effects, the prediction of critical velocities, responses to unbalance, and finally on the prediction of instability thresholds [7, 8]. In general, the damping translates the energy dissipation of the vibratory system [9, 10]. It is therefore in the sense of understanding the influence of damping on the dynamics of gyroscopes that Chen modeled the dynamics of a gyroscope subjected to a cubic dissipation force [2]. The study of the rotating gyroscope under the influence of the nonlinear damping force has made it possible to understand that the gyroscope under given conditions has a complex dynamic [13, 913]. Of these behaviors, the most complex are chaotic movements and the phenomenon of multistability [1]. Several works have attempted to propose modes of control or synchronization of the chaotic gyroscopes according to the field of application of the gyroscope and its intended interest [26, 1117]. Multistability is one of the complex phenomena encountered in nonlinear systems. Thus, the dynamics of systems exhibiting this multistability phenomenon are difficult to predict with precision. Indeed, for the same value of a parameter of the system for which multistability appears, the system is in several states or at several vibration amplitudes, thus, making the system difficult to control. On the contrary, bistability reflects the coexistence of two attractors, while megastability designates the coexistence of an infinite number of attractors for the same system [1825].

Many researchers have endeavored to predict and control multistability because of its impact on the performance of systems that exhibit this phenomenon. Among the methods used for the control of complex phenomena in nonlinear dynamics, there are methods which include active, semiactive, passive, and feedback control. For example, recently, the control of multistability phenomenon of the famous hyperchaotic oscillator TNC has been achieved using a linear augmentation control scheme [26]. Similarly, in [27], the control of up to nine coexisting hidden chaotic attractors in a specific domain of the initial condition is successfully realized. This control is achieved by choosing initially an attractor as a designated survivor using the feedback term method. It is then found that the employed method can also be used to control multistability in a system without an equilibrium point.

From the literature, it should be noted, at least to our knowledge, that the control of vibration amplitude, chaotic dynamics, and coexistence of attractors for rotating gyroscopes has not yet been done using passive control force. However, during its use, the rotating gyroscope may be under the influence of a passive force when it is influenced by vibrating base velocity, for example. Thus, it would be necessary to study complex gyroscope dynamics and control them with passive forces [28, 29]. It is therefore in this perspective that this work proposes to model the dynamics of the rotating gyroscope by taking into account the passive control forces. It is also a question in this work of researching the influence of the control force on the amplitude of the harmonic vibrations of the gyroscope, on the one hand, and of analyzing the effects of each of the parameters of the passive control force on the chaotic behavior and on the coexistence of attractors phenomenon when the gyroscope is working. This paper, therefore, makes it possible to identify in which operating condition of the gyroscope the passive control would be effective.

The paper is structured as follows: In Section 2, we will give the mathematical modeling of the oscillations of the complex dynamics of nonlinear gyroscope and formulate the process of passive control. In Section 3, we will determine the amplitude of the harmonic oscillations. The stability limits are analyzed, and the effect of the different passive control force parameters is also studied on the amplitude of the harmonic resonance. Section 4 discusses the route to chaos and the effects of controlling force on chaotic dynamic states and the coexistence of attractors. The conclusion is presented in the last section.

2. Model and Equation of Oscillations

In this work, we consider the symmetrical gyroscope mounted on a vibrating base [1, 2].

The geometry of the considered problem is illustrated in Figure 1. The motion of a symmetrical gyroscope mounted on a vibrating base can be described by Euler angles (nutation), Փ (precession), and Ψ (rotation).

The Lagrangian of the system is given in the following equation [1, 2]:where and are the polar and equatorial moments of inertia of the symmetrical gyroscope, respectively, is the force of gravity, is the magnitude of the external excitation disturbance, and is the frequency of the external excitation disturbance.

Using Routh’s procedure and using the Routhian associated to hidden coordinates Փ et Ψ, the equation governing the dynamics of the gyroscope as a function of the angle θ is given by [1, 2]:where and . is parametric excitation, and are, respectively, linear and nonlinear damping, and is a nonlinear restoring force. In general, and are fixed parameters, while is the normalized amplitude, is the frequency of the external excitation, and is proportional to the speed of rotation.

In this work, we suppose that the support on which the gyroscope is fixed vibrates with a speed. The additional driving force that applies to the gyroscope can be defined by [29, 30] , where represents the control gain coefficient and indicates the direction of the passive control force . In the particular case where parameter , the control system is reduced to the original equation governing the movement of the gyroscope. For , the drag force follows the positive direction of vibrating base velocity, while directions of vibrating base velocity and forces are opposite when .

The control system can then be written as follows:

3. Harmonic Oscillations’ States and Stability

3.1. Harmonic Vibrations

Due to the highly nonlinear problem in gyroscopes, the nonlinear terms and are extended in power of [2, 4]. So, we getwith

Still for the same reason of the nonlinearity of the considered system, the differential equation of the control system is analytically solved using the harmonic balance method [30]. The general solution of equation (4) iswhere and are amplitudes of oscillations.

Inserting equation (6) into (4) and assuming that and are slowly varying functions with time, then ignoring the second derivatives and the terms of the second order, we obtain the following set of first-order coupled differential equations:withare expressions of amplitude and phase, respectively.

The steady-state vibrations are obtained for and . The characteristic geometric equation is thus given by

From equation (9) after some algebraic manipulations, we have

By combining equations (10) and (11), we get

Figures 24 show the effects of passive control on the amplitudes of the oscillations. These figures are obtained by numerically solving equation (11) with fixed parameters [13], ; ; ;

The effect of controller on system behavior is shown in Figure 2. This figure shows that the amplitude of the harmonic resonance is reduced for both directions of the passive control force. We investigated the effects of control force on the resonant state.

Figures 3 and 4 show the effects of , and on the amplitude of the harmonic oscillations, respectively, for and . It can be observed that the amplitude of the harmonic oscillations decreases when one of parameters and of the passive control force increases with . We can conclude that our passive control is practical because the amplitude of the oscillations of the complex dynamics of the nonlinear gyroscope is considerably reduced when the appropriate control parameters are varied.

3.2. Stability Analysis

An important question is the stability of the stationary solution in the parameter space. To determine which of the stationary solutions are physically acceptable, their stability properties must be studied [1]. For this purpose, we considerwhere and are the values of and at equilibrium states obtained by solving equation (11), and are the disturbances. Inserting (13) in the system of equation (9) and retaining only the linear terms, we obtain the system of linearized equations for the disturbance in the formwithwith of obtained at steady state, respectively, whose expressions are as follows:

Thereby,

Now, searching the solutions of system (17) areand substituting these solutions (18) into equation (17), we have the following characteristic equation:where and are functions of system parameters and are given by

According to the Routh–Hurwitz criterion, equilibrium states will be unstable if the characteristic equation (19) has, at least, one positive real root. Otherwise, they are stable. Thus, the stability domain is given by the equations and .

4. Control of Chaotic Dynamic States and Coexistence of Attractors

4.1. Control of Chaotic Dynamic States

The complex dynamics of the nonlinear gyroscope considered in this work can be chaotic. This chaotic dynamic does not guarantee stability in the movement of the gyroscope [1]. This can reduce the performance of the gyroscope. Our goal in this subsection is to reduce or remove this chaotic dynamics by using a passive control force [31]. To achieve this goal, we used the fourth-order Runge–Kutta algorithm [32] to numerically solve equation (3), and the resulting bifurcation diagrams and its corresponding Lyapunov exponents are plotted using the amplitude of the external excitation forces as the bifurcation parameter. For the simulations, the parameters used are ; the step for the bifurcation parameter is hf = 0.001; the step for time is h = 0.01; the number of iterations is 1000000. The bifurcation diagram and its corresponding Lyapunov exponent are obtained in the absence of the control force Figures 5(a) and 5(b) for . In the presence of the control force, we plotted for the same parameters as Figures 5(a), 5(c), and 5(d) for and Figures 5(e) and 5(f) for . It noticed through these figures that the chaos is accentuated for but totally reduced for [32].

Thus, chaos is well controlled when the gyroscope’s velocity is in the opposite direction than the passive control force. For example, to verify the predictions of Figure 5, we plotted for F = 45 the corresponding phase portraits and time series of the system for , , and (see Figure 6). It was noticed that the phase diagrams and the corresponding times stories confirm the dynamics of Figure 5 and the effectiveness of the control force in the direction opposite to that of the fluid, i.e., for . Figure 7 shows the effect of the control gain parameter and the fluid velocity on the dynamics of the gyroscope for . It emerges from the analysis of this figure, that in this direction, the parameter accentuates the chaotic behavior, while the large values of the fluid velocity make it possible to eliminate the chaotic dynamics of the gyroscope (see Figure 7(f), for example). This important result is easily seen in Figure 8, where it is clearly seen that the parameters and have inverse roles on the dynamics of the system when . Always on this figure, one can observe that the chaos is totally eliminated in the direction of the fluid velocity when the latter takes the value . The study of the effect of and on the dynamics of the gyroscope for revealed that and are very good control parameters because when is fixed, the parameter makes it possible to completely eliminate the chaotic behavior, and for fixed, the speed makes it possible to completely eliminate this undesirable behavior for the stability of the gyroscope. These different results are illustrated in Figure 9 where for , it is noted that chaos is totally eliminated for when is fixed at , while in the same direction, chaos is eliminated for when the gain parameter is set to . These results are confirmed by the phase portraits and their time series presented in Figure 10. It follows from these various analyses that the control of the chaotic dynamics of the gyroscope is very effective when it operates in the opposite direction of the fluid velocity under the influence of the passive control force.

4.2. Control of the Coexistence of Attractors

To further support the results of the theoretical investigations in this work, the research and control of the coexistence of attractors in the dynamic behavior of the gyroscope are introduced. Indeed, the phenomenon of coexistence of attractors is very complex and can affect the performance of the gyroscope. It is therefore, for this reason that we seek and propose the control of this complex phenomenon. For this, we study the influence of the passive control force on this behavior using bifurcation diagrams, Lyapunov exponents, and phase portraits. The bifurcation diagram and the Lyapunov exponents are obtained by choosing as the bifurcation parameter.

On these diagrams, the curve in blue represents the diagram obtained by varying the amplitude of the external excitation in the increasing direction, while the results obtained by varying in the decreasing direction while remaining in the same interval is the curve in red. When the two curves have the same nature and are exactly superimposed (merged) throughout the domain, then the gyroscope takes the same path on the outward as on the return side when parameter is varied in both directions. If the two curves do not overlap exactly over a given interval of , then the path followed by the gyroscope in its outward movement is not the same as that which it takes in the return for increasing and decreasing evolutions from . Attractors are said to coexist. In this condition, if the dynamics of the gyroscope are the same in this domain, we say that attractors of the same nature coexist; otherwise, we speak of the coexistence of attractors of different natures. Thus, for the same value of the bifurcation parameter, the vibrations of the gyroscope will have different amplitudes and of different natures depending on the case. It is therefore a very complex phenomenon because, in this state, it would be difficult to identify exactly and uniquely the amplitude and the nature of the vibrations of this gyroscope. This could have a negative impact on the performance of the rotating gyroscope because it would complicate the precise acceleration of Coriolis when the gyroscope will enter this phase.

Figure 11 represents the dynamics of the rotating gyroscope when for between 30 and 50. Thus, we note through this figure the coexistence of periodic attractors of period , of period , of chaotic attractors, and of attractors of period with chaotic attractors. In the presence of the passive control force (Figure 11), we note when for between 30 and 50, that the coexistence of attractors of period has disappeared giving way to the coexistence of attractors of period with those of period . Similarly, the coexistence of attractors of period with chaotic attractors is eliminated, leaving room for the coexistence of attractors of period with chaotic attractors. For , we note that the coexistence of attractors of different natures is totally eliminated and the domain of the coexistence of chaotic attractors is considerably reduced. We can therefore say that the control of the phenomenon of coexistence of attractors is more effective when the control force is applied in the opposite direction of the fluid velocity than in the same direction . Figure 12 shows the effects of the parameters and on the phenomenon of the coexistence of attractors for . It is observed that the parameter does not facilitate the elimination of this phenomenon, whereas for a high value of the velocity of the fluid, the coexistence of attractors of different natures as well as chaotic attractors is eliminated. We finally obtain, in this case, the coexistence of attractors of period . Finally, Figure 13 represents the effects of the parameters and on the phenomenon of the coexistence of attractors for . In this case, we note that the parameter of control and the velocity of the fluid eliminate the coexistence of the attractors but preserve the coexistence of the attractors of period . Figures 1416 represent the phase portraits of the system, respectively, for . These figures confirm the existence and control of the coexistence of the attractors predicted by Figures 1113.

In sum, from the analysis of the results obtained, it is clear that the control force makes it possible to eliminate the coexistence of attractors of different, chaotic natures but preserves the coexistence of attractors of period . We also note that this control is effective for the direction of than in the direction . Moreover, the speed of the fluid makes it possible more to control this complex phenomenon contrary to K.

5. Conclusions

In this work, we investigated the passive control of complex nonlinear gyroscope dynamics by passive control force. The model has been described, and the equation of the control system governing the corresponding movement is obtained. Using the method of harmonic balance, we found the amplitude of oscillatory states, the Routh–Hurwitz criterion is used to analyze and determine the stability limits of oscillations. Numerical simulations are used to validate and complete the results obtained by analytical methods. We noted a more robust stability of the complex dynamics of the gyroscope subjected to the passive control force in the negative direction of the fluid velocity than in the positive direction. In both directions, we observe the existence of instability domains. We also noticed the increase or decrease of these instability domains when the control parameters of the control force of the system increase or decrease. The effect of the control process on chaotic dynamic and on coexistence of attractors was effective with [32]. The high amplitude of harmonic oscillations, chaotic states, and coexistence of attractors was successfully controlled by the passive control investigated in this work. Advances in the design and manufacture of the rotating machines widely used today in many sectors of industry make it possible to increase both the performance and the efficiency of these machines by operating them in speed ranges increasingly high turnover. However, the forces generated are increasingly important and strongly solicit the overall dynamic behavior of the machines. Thus, the results obtained in this work are useful for the study of the control of the dynamics of rotating machines which is more relevant than ever. Precisely, this work has a capital importance in engineering because the results obtained allow to have a stable dynamics for a better performance of the gyroscope.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors would like to thank The World Academy of Sciences (TWAS) for financial support through Research Grant Agreement No. 20-307RG/PHYS/AF/ACG-FR3240314170.