Next Article in Journal
Hybrid VOF–Lagrangian CFD Modeling of Droplet Aerobreakup
Next Article in Special Issue
Assessing the Seismic Demands on Non-Structural Components Attached to Reinforced Concrete Frames
Previous Article in Journal
Method for the Analysis of Health Personnel Availability in a Pandemic Crisis Scenario through Monte Carlo Simulation
Previous Article in Special Issue
Measurement of Orthotropic Material Constants and Discussion on 3D Printing Parameters in Additive Manufacturing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Resonant Vibrations Control of the Nonlinear Rotor Active Magnetic Bearing Systems

1
Department of Physics and Engineering Mathematics, Faculty of Electronic Engineering, Menoufia University, Menouf 32952, Egypt
2
Department of Basic Science, Higher Technological Institute, Tenth of Ramadan City 44629, Egypt
3
Mathematics Department, Faculty of Science, Al-Azhar University, Nasr City 11651, Egypt
4
Department of Automation, Biomechanics and Mechatronics, Faculty of Mechanical Engineering, Lodz University of Technology, 90924 Lodz, Poland
5
Department of Mathematics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
*
Author to whom correspondence should be addressed.
Submission received: 4 July 2022 / Revised: 5 August 2022 / Accepted: 16 August 2022 / Published: 19 August 2022

Abstract

:
Nonlinear vibration control of the twelve-poles electro-magnetic suspension system was tackled in this study, using a novel control strategy. The introduced control algorithm was a combination of three controllers: the proportional-derivative ( P D ) controller, the integral resonant controller ( I R C ) , and the positive position feedback ( P P F ) controller. According to the presented control algorithm, the mathematical model of the controlled twelve-poles rotor was established as a nonlinear four-degree-of-freedom dynamical system coupled to two first-order filters. Then, the derived nonlinear dynamical system was analyzed using perturbation analysis to extract the averaging equations of motion. Based on the extracted averaging equations of motion, the efficiency of different control strategies (i.e., P D ,   P D + I R C ,   P D + P P F ,   and   P D + I R C + P P F ) for mitigating the rotor’s undesired vibrations and improving its catastrophic bifurcation was investigated. The acquired analytical results demonstrated that both the P D and P D + I R C controllers can force the rotor to respond as a linear system; however, the controlled system may exhibit the maximum oscillation amplitude at the perfect resonance condition. In addition, the obtained results demonstrated that the P D + P P F controller can eliminate the rotor nonlinear oscillation at the perfect resonance, but the system may suffer from high oscillation amplitudes when the resonance condition is lost. Moreover, we report that the combined control algorithm ( P D + I R C + P P F ) has all the advantages of the individual control algorithms (i.e., P D ,   P D + I R C ,   P D + P P F ), while avoiding their drawbacks. Finally, the numerical simulations showed that the P D + I R C + P P F controller can eliminate the twelve-poles system vibrations regardless of both the excitation force magnitude and the resonant conditions at a short transient time.

1. Introduction

Vibration analysis and control of the electro-magnetic suspension system are among the most important research topics for scientists and engineers worldwide. The importance of this suspension system is due to its many industrial applications, including its use in rotor dynamics and in the automobile industries. The rotor electro-magnetic suspension system is a special type of active bearing that is used to support the rotating shafts without any physical contact with the stator parts of the system. The working principle of rotor electro-magnetic suspension is the application of controllable electro-magnetic attractive forces to support the rotating shafts in their hovering positions via compensating for the external loads that are exerted on these shafts. The operation of the rotating shafts without physical contact with the stators gives this suspension system many preferable features when compared with conventional bearings systems, such as less maintenance, no need for lubrication between the rotors and stators, a clean working environment, high operational speed, high reliability, and high durability. Accordingly, many research articles have investigated the dynamical characteristics of different configurations of the rotor electro-magnetic suspension system.
Different control algorithms have been proposed to enhance the vibratory characteristics and eliminate the undesired nonlinear bifurcation behaviors of this suspension system. Ji et al. [1] studied the nonlinear dynamics and motion bifurcations of a rotor electro-magnetic suspension system consisting of a four-poles configuration. They established the mathematical model that governs the rotor lateral vibrations as a two-degree-of-freedom nonlinear dynamical system. Then, they investigated the derived equations of motion using the multiple-time scales perturbation method. Based on their analysis, they reported that the rotor system may lose its stability either via saddle-node or Hopf bifurcations.
Saeed et al [2] investigated the vibratory characteristics of a rotor supported by a six-poles electro-magnetic suspension system. They introduced two control strategies utilizing the P D -control algorithm. The first control technique was established based on the Cartesian displacements and velocities of the rotor in the horizontal and vertical directions, while the second control technique was designed according to the radial oscillations of the rotor in the direction of the six poles. Based on their analysis, they reported that the rotor system may lose its stability and exhibit unbounded oscillation in the case of the radial control technique at a specific value of the proportional gain. In addition, they showed that the system may perform either a quasi-periodic or chaotic response in the case of the Cartesian control strategy at a strong excitation force.
Ji and Hansen [3,4] studied the nonlinear dynamics of a rotor supported by an eight-poles electro-magnetic suspension system. They applied the Cartesian PD-control strategy to improve the vibratory characteristics of that system at both primary [3] and super-harmonic resonance conditions [4]. They reported that the eight-poles system has bi-stable and tri-stable solutions. In addition, they showed that the system may be exposed to a multi-jump when the rotor angular speed crosses its first critical speed.
El-Shourbagy et al. [5] introduced a nonlinear PD-control algorithm to enhance the nonlinear lateral vibrations of a rotor supported by an eight-poles electro-magnetic suspension system. Saeed et al. [6] explored numerically the motion bifurcations of a rotor system supported by the eight poles when the rub-impact force between the rotor and stator occurs. They illustrated that the rotor may execute either full annular rub mode or rub-impact motion, depending on both the impact stiffness and the dynamic friction coefficients. In addition, Zhang et al. [7,8,9,10,11,12] introduced detailed investigations of the eight-poles electro-magnetic suspension system with variable stiffness coefficients. The nonlinear dynamical behaviors of the twelve-poles electro-magnetic suspension system were investigated utilizing the P D -control algorithm for the first time by El-Shourbagy et al. [13]. They reported that proportional control gain can play an important role in reshaping the system dynamics. In addition, they demonstrated that the twelve-poles system may lose its stability at a strong excitation force. Saeed et al. [14] explored the dynamical characteristics of the sixteen-poles system with constant stiffness coefficients utilizing the conventional P D -control algorithm. Zhang et al. [15,16,17,18] introduced extensive investigations for the sixteen-poles rotor system with time-varying stiffness coefficients. Due to the controllability and flexibility of the rotor electro-magnetic suspension system, it was used as an active actuator to control the dynamical behaviors of some rotating machines [19,20,21,22,23].
The positive position feedback ( P P F ) control algorithm has been applied extensively to eliminate the resonant vibrations of many dynamical systems [24,25,26,27,28]. Saeed et al. [28] utilized the P P F -control strategy with a P D controller to mitigate the undesired vibrations of the eight-poles rotor system for the first time. They concluded that the P P F controller can eliminate the system’s lateral vibration at the perfect resonance condition. However, the main drawback of this control strategy was that the controller may add excessive vibratory motion to the rotor system if the tuning condition was lost. In addition, the integral resonance controller ( I R C ) was one of the feasible control methods that was applied to mitigate the undesired vibrations and eliminate the nonlinear bifurcations of different dynamical systems [29,30,31,32,33,34,35,36]. Recently, Saeed et al. [36] introduced the I R C -control algorithm for the first time to mitigate the unwanted vibrations of the eight-poles rotor system. They reported that the I R C controller can reduce the system’s vibrations and suppress the corresponding catastrophic bifurcations. However, the main drawback of this control method was that the I R C -controller could not eliminate the rotor vibrations at a resonance condition close to zero.
In the present work, a new control strategy is introduced to eliminate the nonlinear lateral vibrations of the twelve-poles rotor system. The proposed controller is a combination of the three control algorithms: P D ,   I R C , and P P F . Accordingly, the whole-system mathematical model is derived as a four-degree-of-freedom dynamical system that is coupled to two first-order differential equations. Then, the system dynamical model is analyzed, and the corresponding slow-flow modulation equations are extracted. Based on the obtained slow-flow modulation equation, the performance of the suggested control technique is explored. The obtained analytical results showed that the P D ,   I R C , and P D + I R C controllers can mitigate the nonlinear oscillation of the system and force the rotor to respond as a linear system. but the main drawback of these types of controllers (i.e., P D ,   I R C , and P D + I R C ) is that the controlled system may perform the maximum oscillation amplitude at the resonant condition. In addition, we found that the coupling of the P D + P P F controller to the system can eliminate the rotor’s undesired oscillation at the perfect resonance, but the system may suffer from high oscillation amplitudes if the resonance condition is lost. Moreover, the acquired analytical and numerical investigations demonstrated that the P D + I R C + P P F controller has all the advantages of the individual control algorithms (i.e., P D ,   P D + I R C , and I R C + P P F ), while avoiding their drawbacks.

2. Equations of Motion

The studied rotor system is assumed to be a rigid body with a two-degree-of-freedom system that has mass m and eccentricity e and rotates with angular velocity ψ , as shown in Figure 1. In addition, this rotor system is supported in its nominal position via the restoring forces f x and f y that are generated by twelve electro-magnetic poles. Therefore, the system equations of motion can be expressed as follows [37,38]:
m x ¨ f x = m e ψ 2 cos ( ψ t )
m y ¨ f y = m e ψ 2 sin ( ψ t )
where f x and f y represent the resultant restoring forces of the twelve poles in both the X and Y directions, respectively. In this study, the attractive forces f j ( j = 1 , 2 , ,   6 ) are designed so that each adjacent pair of the poles generates a push-pull attractive force. Therefore, f j ( j = 1 , 2 , ,   6 ) can be expressed according to the electro-magnetic theory, as follows [38]:
f j = Θ [ ( I 0 I j ) 2 ( c 0 δ j ) 2 ( I 0 + I j ) 2 ( c 0 + δ j ) 2 ] , j = 1 , 2 , , 6
where Θ = 1 4 μ 0 N 2 Acos ( φ ) is constant, I 0 is constant current defined as a bias current, I j ( j = 1 , 2 , , 6 ) is the control currents that will be defined later according to the purposed control algorithm, c 0 is the nominal air-gap size between the rotor and the twelve poles, and δ j is the radial deviation of the rotor away from the geometric center O in the direction of the j t h pole.
Based on the system’s geometry as shown in Figure 1, for the small temporal Cartesian displacements x ( t ) and y ( t ) of the rotor in both the X and Y directions, one can express the radial displacements δ j ,   ( j = 1 , 2 , , 6 ) of the rotor system as follows:
δ 1 ( x , y ) = x ( t ) cos ( α ) y ( t ) sin ( α ) , δ 2 ( x , y ) = x ( t ) , δ 3 ( x , y ) = x ( t ) cos ( α ) + y ( t ) sin ( α ) , δ 4 ( x , y ) = x ( t ) sin ( α ) + y ( t ) cos ( α ) , δ 5 ( x , y ) = y ( t ) , δ 6 ( x , y ) = x ( t ) sin ( α ) + y ( t ) cos ( α ) }
where α is the angle between every two consecutive poles (i.e., α = 360 / 12 = 30 ). In this work, the control currents were designed so that the control forces f 1 ,   f 2 , and f 3 depend on the horizontal displacement x ( t ) , while the forces f 4 ,   f 5 , and f 6 depend on the vertical displacement y ( t ) . Accordingly, the control currents I j ( j = 1 , 2 , , 6 ) are selected as follows:
I X = I 1 = I 2 = I 3 , I Y = I 4 = I 5 = I 6
where I X is the control current that is responsible for eliminating the nonlinear oscillations of the rotor system in the X direction, while I Y is the control current that is responsible for eliminating the nonlinear oscillations of the rotor system in the Y direction. Accordingly, to eliminate the undesired vibrations of the system, an advanced control strategy was introduced. The suggested control method is a combination of three control algorithms: the P D controller, the I R C controller, and the P P F controller. Therefore, the control laws (i.e., control currents I X and I Y ) are designed as follows:
I X = k 1 x + k 2 x ˙ k 3 u 1 + k 4 u 2 , I Y = k 1 y + k 2 y ˙ k 5 v 1 + k 6 v 2
where k 1 and k 2 are the control gains of the P D controller, k 3 and k 5 denote the control gains of the P P F controller, and k 4 and k 6 represent the control gains of the I R C controller. Accordingly, k 1 x + k 2 x ˙ and k 1 y + k 2 y ˙ are the components of the control currents ( I X and I Y ) due to the P D controller in the X and Y directions, respectively, k 3 u 1 and k 5 v 1 are the components of the control currents due to the P P F controller in the X and Y directions, respectively, while + k 4 u 2 and + k 6 v 2 denote the control current components due to the I R C controller in the X and Y directions, respectively. The equations of motion that describe the oscillatory behaviors of the P P F controllers are provided as follows [24,25,26,27,28]:
u ¨ 1 + c 1 u ˙ 1 + λ 1 u 1 = L 1 x
v ¨ 1 + c 2 v ˙ 1 + λ 2 v 1 = L 2 y
where c 1 and c 2 denote the damping coefficients of the P P F controllers, λ 1 and λ 2 represent the controller’s natural frequencies, and L 1 and L 2 are the feedback signals gains. In addition, the dynamical behaviors of the I R C controllers are governed by first-order differential equations that are provided as follows [29,30,31,32,33,34,35,36]:
u ˙ 2 + λ 3 u 2 = L 3 x
v ˙ 2 + λ 4 v 2 = L 4 y
where λ 3 and λ 4 denote the internal feedback gain of the I R C controller, and L 3 and L 4 represent the feedback signals gains. The interconnection between the twelve-poles system and the proposed control algorithm (i.e., the PD + IRC + PPF controllers) is illustrated schematically in Figure 2, where the temporal Cartesian oscillations (i.e., x ( t ) and y ( t ) ) of the rotor in both the X and Y directions can be measured using two position sensors that may be fixed on the poles-housing in the + X and + Y directions, as shown in Figure 1a. Then, the measured signals, x ( t ) and y ( t ) , are fed into a digital computer on which the control algorithm (i.e., the PD + IRC + PPF controller) is implemented. According to the programmed algorithm, the controller computes the control currents I X = k 1 x + k 2 x ˙ k 3 u 1 + k 4 u 2 and I Y = k 1 y + k 2 y ˙ k 5 v 1 + k 6 v 2 , as shown in Figure 2. Finally, the computed control currents are applied to a power amplifiers network to energize the twelve-poles electrical coils in order to generate the electro-magnetic forces ( f 1 ,   f 2 ,   ,   f 6 ), which in turn try to mitigate the lateral oscillations, x ( t ) and y ( t ) , of the rotor system.
Now, to investigate the performance of the proposed closed-loop system, the whole-system model should be obtained and then analyzed to report the optimum working conditions of this system. Therefore, by substituting Equations (4) to (6) into Equation (3), we have the following:
f 1 = Θ ( ( I 0 k 1 x k 2 x ˙ + k 3 u 1 k 4 u 2 ) 2 ( c 0 x cos ( α ) + y sin ( α ) ) 2 ( I 0 + k 1 x + k 2 x ˙ k 3 u 1 + k 4 u 2 ) 2 ( c 0 + x cos ( α ) y sin ( α ) ) 2 )
f 2 = Θ ( ( I 0 k 1 x k 2 x ˙ + k 3 u 1 k 4 u 2 ) 2 ( c 0 x ) 2 ( I 0 + k 1 x + k 2 x ˙ k 3 u 1 + k 4 u 2 ) 2 ( c 0 + x ) 2 )
f 3 = Θ ( ( I 0 k 1 x k 2 x ˙ + k 3 u 1 k 4 u 2 ) 2 ( c 0 x cos ( α ) y sin ( α ) ) 2 ( I 0 + k 1 x + k 2 x ˙ k 3 u 1 + k 4 u 2 ) 2 ( c 0 + x cos ( α ) + y sin ( α ) ) 2 )
f 4 = Θ ( ( I 0 k 1 y k 2 y ˙ + k 5 v 1 k 6 v 2 ) 2 ( c 0 x sin ( α ) y cos ( α ) ) 2 ( I 0 + k 1 y + k 2 y ˙ k 5 v 1 + k 6 v 2 ) 2 ( c 0 + x sin ( α ) + y cos ( α ) ) 2 )
f 5 = Θ ( ( I 0 k 1 y k 2 y ˙ + k 5 v 1 k 6 v 2 ) 2 ( c 0 y ) 2 ( I 0 + k 1 y + k 2 y ˙ k 5 v 1 + k 6 v 2 ) 2 ( c 0 + y ) 2 )
f 6 = Θ ( ( I 0 k 1 y k 2 y ˙ + k 5 v 1 k 6 v 2 ) 2 ( c 0 + x sin ( α ) y cos ( α ) ) 2 ( I 0 + k 1 y + k 2 y ˙ k 5 v 1 + k 6 v 2 ) 2 ( c 0 x sin ( α ) + y cos ( α ) ) 2 )
Based on the system geometry, as shown in Figure 1, the resultant attractive forces f x and f y in the X and Y directions due to the forces f 1 ,   f 2 ,   ,   f 6 can be expressed as follows:
f x = f 2 + ( f 1 + f 3 ) cos ( α ) + ( f 4 f 6 ) sin ( α )
f y = f 5 + ( f 4 + f 6 ) cos ( α ) + ( f 3 f 1 ) sin ( α )
To simplify the rational form of the attractive forces f 1 ,   f 2 ,   ,   f 6 , Equations (11) to (16) were expanded, using the Maclaurin series, up to the third order approximation, as provided in Appendix A. Now, substituting the expanded Equations (A1) to (A6) that are provided in Appendix A into Equations (17) and (18), then inserting the resulting equations into Equations (1) and (2) and introducing the dimensionless parameters: t * = ϑ t , z 1 = x c 0 , z 1 ˙ = x ϑ c 0 ˙ , z 1 ¨ = x ¨ ϑ 2 c 0 , z 2 = y c 0 , z 2 ˙ = y ϑ c 0 ˙ , z 2 ¨ = y ¨ ϑ 2 c 0 , z 3 = u 1 c 0 , z 3 ˙ = u 1 ϑ c 0 ˙ , z 3 ¨ = u 1 ¨ ϑ 2 c 0 , z 4 = v 1 c 0 , z 4 ˙ = v 1 ϑ c 0 ˙ , z 4 ¨ = v 1 ¨ ϑ 2 c 0 , z 5 = u 2 c 0 , z 5 ˙ = u 2 ϑ c 0 ˙ , z 6 = v 2 c 0 , z 6 ˙ = v 2 ϑ c 0 ˙ , ω 1 = λ 1 ϑ 2 , ω 2 = λ 2 ϑ 2 , ω 3 = λ 3 ϑ , ω 4 = λ 4 ϑ , Ω = ψ ϑ , f = e c 0 , p = c 0 I 0 k 1 , d = c 0 ϑ I 0 k 2 , μ 1 = c 1 2 ϑ , μ 2 = c 2 2 ϑ , η 1 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 3 , η 2 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 4 , η 3 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 5 , η 4 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 6 , η 5 = L 1 ϑ 2 , η 6 = L 2 ϑ 2 , η 7 = L 3 ϑ , η 8 = L 4 ϑ , ϑ = Θ / m c 0 3 , one can obtain the following dimensionless equations of motion that govern the nonlinear dynamics of the proposed closed-loop system:
z ¨ 1 + 2 μ z ˙ 1 + ω 2 z 1 ( α 1 z 1 3 + α 2 z 1 z 2 2 + α 3 z 1 2 z ˙ 1 + α 4 z ˙ 1 z 2 2 + α 5 z 1 z ˙ 2 2 + α 6 z 1 z ˙ 1 2 + α 7 z 1 z 2 z ˙ 2 + β 1 z 1 2 z 3 + β 2 z 1 z ˙ 1 z 3 + β 3 z 1 z 3 2 + β 4 z 1 z 2 z 4 + β 5 z 1 z 4 2 + β 6 z 1 z ˙ 2 z 4 + β 7 z 3 z 2 2 + β 8 z 1 z 3 z 5 + β 9 z 1 z ˙ 1 z 5 + β 10 z 1 z 4 z 6 + β 11 z 1 2 z 5 + β 12 z 1 z 6 2 + β 13 z 1 z 2 z 6 + β 14 z 1 z ˙ 2 z 6 + β 15 z 1 z 5 2 + β 16 z 2 2 z 5 ) = Ω 2 f cos ( Ω t ) + η 1 z 3 + η 2 z 5
z ¨ 2 + 2 μ z ˙ 2 + ω 2 z 2 - ( α 1 z 2 3 + α 2 z 2 z 1 2 + α 3 z 2 2 z ˙ 2 + α 4 z ˙ 2 z 1 2 + α 5 z 2 z ˙ 1 2 + α 6 z 2 z ˙ 2 2 + α 7 z 2 z 1 z ˙ 1 + γ 1 z 2 2 z 4 + γ 2 z 2 z ˙ 2 z 4 + γ 3 z 2 z 4 2 + γ 4 z 2 z 1 z 3 + γ 5 z 2 z 3 2 + γ 6 z 2 z ˙ 1 z 3 + γ 7 z 4 z 1 2 + γ 8 z 2 z 4 z 6 + γ 9 z 2 z ˙ 2 z 6 + γ 10 z 2 z 3 z 5 + γ 11 z 2 2 z 6 + γ 12 z 2 z 5 2 + γ 13 z 2 z 1 z 5 + γ 14 z 2 z ˙ 1 z 5 + γ 15 z 2 z 6 2 + γ 16 z 1 2 z 6 ) = Ω 2 f sin ( Ω t ) + η 3 z 4 + η 4 z 6
z ¨ 3 + 2 μ 1 z ˙ 3 + ω 1 2 z 3 = η 5 z 1
z ¨ 4 + 2 μ 2 z ˙ 4 + ω 2 2 z 4 = η 6 z 2
z ˙ 5 + ω 3 z 5 = η 7 z 1
z ˙ 6 + ω 4 z 6 = η 8 z 2
Equations (19) and (20) represent the dimensionless equations of motion of the controlled twelve-poles system, while Equations (21) and (22) are the dimensionless equations of motion of the P P F controller. In addition, Equations (23) and (24) are the dimensionless equations of motion of the I R C controller. Accordingly, the suggested closed-loop system is governed by six-coupled nonlinear ordinary differential equations, four of which are of the second order and the other two of which are of the first order, where the coefficients of the above six equations are provided in Appendix B.

3. Analytical Investigations

Many analytical methods have been introduced in the literature to investigate both the linear and nonlinear vibration problems [39,40,41]. Accordingly, to explore the efficiency of the introduced closed-loop system, we sought an approximate solution for the system equations of motions (i.e., Equations (19) to (24)) within this section, in the form of a first-order perturbation series as follows [39,40]:
z 1 ( t , ε ) = z 1 0 ( T 0 , T 1 ) + ε z 1 1 ( T 0 , T 1 )
z 2 ( t , ε ) = z 2 0 ( T 0 , T 1 ) + ε z 2 1 ( T 0 , T 1 )
z 3 ( t , ε ) = z 3 0 ( T 0 , T 1 ) + ε z 3 1 ( T 0 , T 1 )
z 4 ( t , ε ) = z 4 0 ( T 0 , T 1 ) + ε z 4 1 ( T 0 , T 1 )
z 5 ( t , ε ) = ε z 5 0 ( T 0 , T 1 ) + ε 2 z 5 1 ( T 0 , T 1 )
z 6 ( t , ε ) = ε z 6 0 ( T 0 , T 1 ) + ε 2 z 6 1 ( T 0 , T 1 )
where T 0 = t , T 1 = ε t , and ε is the perturbation parameter that was used as a bookkeeping coefficient during this analysis [40]. According to the introduced two-time scales (i.e., T 0 ,   T 1 ), the ordinary derivatives d d t and d 2 d t 2 should be re-written as follows:
d d t = D 0 + ε D 1 , d 2 d t 2 = D 0 2 + 2 ε D 0 D 1 , and   D j = T j , j = 0 , 1
In addition, to perform the perturbation analysis using ε as a bookkeeping coefficient, the parameters of Equations (19) to (24) should be re-scaled as follows:
f = ε f ˜ , μ = ε μ ˜ , μ 1 = ε μ ˜ 1 , μ 2 = ε μ ˜ 2 , α j = ε α ˜ j , β j = ε β ˜ j , γ j = ε γ ˜ j , η k = ε η ˜ k ; j = 1 , , 7 , k = 1 , 3 , 5 , 6 , 7 , 8
Then, by substituting Equations (25) to (32) into Equations (19) to (24), we have O ( ε 0 ) :
( D 0 2 + ω 2 ) z 1 0 = 0
( D 0 2 + ω 2 ) z 2 0 = 0
( D 0 2 + ω 1 2 ) z 3 0 = 0
( D 0 2 + ω 2 2 ) z 4 0 = 0
O ( ε 1 ) :
( D 0 2 + ω 2 ) z 1 1 = 2 D 0 D 1 z 1 0 2 μ ˜ D 0 z 1 0 + α ˜ 1 z 10 3 + α ˜ 2 z 1 0 z 20 2 + α ˜ 3 z 10 2 D 0 z 1 0 + α ˜ 4 z 20 2 D 0 z 1 0 + α ˜ 5 z 1 0 ( D 0 z 2 0 ) 2 + α ˜ 6 z 1 0 ( D 0 z 1 0 ) 2 + α ˜ 7 z 1 0 z 2 0 D 0 z 2 0 + β ˜ 1 z 10 2 z 30 + β ˜ 2 z 1 0 D 0 z 1 0 z 30 + β ˜ 3 z 1 0 z 30 2 + β ˜ 4 z 1 0 z 2 0 z 40 + β ˜ 5 z 1 0 z 40 2 + β ˜ 6 z 1 0 D 0 z 2 0 z 40 + β ˜ 7 z 20 2 z 30 + β 8 z 1 0 z 30 z 50 + β 9 z 1 0 D 0 z 1 0 z 50 + β 10 z 1 0 z 40 z 60 + β 11 z 10 2 z 50 + β 12 z 1 0 z 60 2 + β 13 z 1 0 z 2 0 z 60 + β 14 z 1 0 D 0 z 2 0 z 60 + β 15 z 1 0 z 50 2 + β 16 z 20 2 z 50 + Ω 2 f ˜ cos ( Ω t ) + η ˜ 1 z 30 + η 2 z 50
( D 0 2 + ω 2 ) z 2 1 = 2 D 0 D 1 z 20 2 μ ˜ D 0 z 20 + α ˜ 1 z 20 3 + α ˜ 2 z 20 z 10 2 + α ˜ 3 z 20 2 D 0 z 20 + α ˜ 4 z 10 2 D 0 z 20 + α ˜ 5 z 20 ( D 0 z 10 ) 2 + α ˜ 6 z 20 ( D 0 z 20 ) 2 + α ˜ 7 z 20 z 10 D 0 z 10 + γ ˜ 1 z 20 2 z 40 + γ ˜ 2 z 20 D 0 z 20 z 40 + γ ˜ 3 z 20 z 40 2 + γ ˜ 4 z 20 z 10 z 30 + γ ˜ 5 z 20 z 30 2 + γ ˜ 6 z 20 D 0 z 10 z 30 + γ ˜ 7 z 10 2 z 40 + γ 8 z 20 z 40 z 60 + γ 9 z 20 D 0 z 20 z 60 + γ 10 z 20 z 30 z 60 + γ 11 z 20 2 z 60 + γ 12 z 20 z 50 2 + γ 13 z 20 z 10 z 50 + γ 14 z 20 D 0 z 10 z 50 + γ 15 z 20 z 60 2 + γ 16 z 10 2 z 60 + Ω 2 f ˜ s i n ( Ω t ) + η ˜ 3 z 40 + η 4 z 60
( D 0 2 + ω 1 2 ) z 31 = 2 D 0 D 1 z 30 2 ε μ ˜ 1 D 0 z 30 + η ˜ 5 z 10
( D 0 2 + ω 2 2 ) z 41 = 2 D 0 D 1 z 40 2 μ ˜ 2 D 0 z 40 + η ˜ 6 z 20
( D 0 + ω 3 ) z 50 = η ˜ 7 z 10
( D 0 + ω 4 ) z 60 = η ˜ 8 z 20
The steady-state periodic solutions of Equations (33) to (36), (41), and (42) can be written as follows:
z 10 ( T 0 , T 1 ) = A 1 ( T 1 ) e i ω T 0 + A ¯ 1 ( T 1 ) e i ω T 0
z 20 ( T 0 , T 1 ) = A 2 ( T 1 ) e i ω T 0 + A 2 ¯ ( T 1 ) e i ω T 0
z 30 ( T 0 , T 1 ) = B 1 ( T 1 ) e i ω 1 T 0 + B 1 ¯ ( T 1 ) e i ω 1 T 0
z 40 ( T 0 , T 1 ) = B 2 ( T 1 ) e i ω 2 T 0 + B 2 ¯ ( T 1 ) e i ω 2 T 0
z 50 ( T 0 , T 1 ) = ρ 1 A 1 ( T 1 ) e i ω T 0 + ρ ¯ 1 A ¯ 1 ( T 1 ) e i ω T 0
z 60 ( T 0 , T 1 ) = ρ 2 A 2 ( T 1 ) e i ω T 0 + ρ ¯ 2 A ¯ 2 ( T 1 ) e i ω T 0
where i = 1 , ρ 1 = ω 3 i ω ω 3 2 + ω 2 η ˜ 7 , ρ ¯ 1 = ω 3 + i ω ω 3 2 + ω 2 η ˜ 7 , ρ 2 = ω 4 i ω ω 4 2 + ω 2 η ˜ 8 , ρ ¯ 2 = ω 4 + i ω ω 4 2 + ω 2 η ˜ 8 . A 1 ( T 1 ) ,   A 2 ( T 1 ) ,   B 1 ( T 1 ) and B 2 ( T 1 ) are unknown that will be defined later. A ¯ 1 ( T 1 ) ,   A ¯ 2 ( T 1 ) ,   B ¯ 1 ( T 1 ) , and B ¯ 2 ( T 1 ) are the complex conjugate forms of A 1 ( T 1 ) ,   A 2 ( T 1 ) ,   B 1 ( T 1 ) , and B 2 ( T 1 ) , respectively. Inserting Equations (43) to (48) into Equations (37) to (40), we have the following:
( D 0 2 + ω 2 ) z 11 = ( 2 i ω D 1 A 1 2 i μ ˜ ω A 1 + 3 α ˜ 1 A 1 2 A ¯ 1 + 2 α ˜ 2 A 1 A 2 A ¯ 2 + α ˜ 2 A ¯ 1 A 2 2 + i α ˜ 3 ω A 1 2 A ¯ 1 + 2 i α ˜ 4 ω A 1 A 2 A ¯ 2 i α ˜ 4 ω A ¯ 1 A 2 2 + 2 α ˜ 5 ω 2 A 1 A 2 A ¯ 2 α ˜ 5 ω 2 A ¯ 1 A 2 2 + α ˜ 6 ω 2 A 1 2 A ¯ 1 + i α ˜ 7 ω A ¯ 1 A 2 2 + 2 β ˜ 3 A 1 B 1 B ¯ 1 + 2 β ˜ 5 A 1 B 2 B ¯ 2 + 2 β 11 ρ 1 A 1 2 A ¯ 1 + β 12 ρ 2 2 A ¯ 1 A 2 2 + β 13 ρ 2 A 1 A 2 A ¯ 2 + β 13 ρ 2 A ¯ 1 A 2 2 i β 14 ω ρ 2 A 1 A 2 A ¯ 2 + i β 14 ω ρ 2 A ¯ 1 A 2 2 + β 15 ρ 1 2 A 1 2 A ¯ 1 + 2 β 16 ρ 1 A 1 A 2 A ¯ 2 + η ˜ 2 ρ 1 A 1 ) e i ω T 0 + η ˜ 1 B 1 e i ω 1 T 0 + ( α ˜ 1 A 1 3 + α ˜ 2 A 1 A 2 2 + i α ˜ 3 ω A 1 3 + i α ˜ 4 ω A 1 A 2 2 α ˜ 5 ω 2 A 1 A 2 2 α ˜ 6 ω 2 A 1 3 + i α ˜ 7 ω A 1 A 2 2 + i β 9 ω ρ 1 A 1 3 + β 11 ρ 1 A 1 3 + β 12 ρ 2 2 A 1 A 2 2 + β 13 ρ 2 A 1 A 2 2 + i β 14 ω ρ 2 A 1 A 2 2 + β 15 ρ 1 2 A 1 3 + β 16 ρ 1 A 1 A 2 2 ) e 3 i ω T 0 + ( β ˜ 1 A 1 2 B 1 + i β ˜ 2 ω A 1 2 B 1 + β ˜ 7 A 2 2 B 1 + β 8 ρ 1 A 1 2 B 1 ) e i ( 2 ω + ω 1 ) T 0 + ( β ˜ 1 A 1 2 B ¯ 1 + i β ˜ 2 ω A 1 2 B ¯ 1 + β ˜ 7 A 2 2 B ¯ 1 + β 8 ρ 1 A 1 2 B ¯ 1 ) e i ( 2 ω ω 1 ) T 0 + ( 2 β ˜ 1 A 1 A ¯ 1 B 1 + 2 β ˜ 7 A 2 A ¯ 2 B 1 + β 8 ρ 1 A 1 A ¯ 1 B 1 ) e i ω 1 T 0 + β ˜ 3 A 1 B 1 2 e i ( ω + 2 ω 1 ) T 0 + β ˜ 3 A ¯ 1 B 1 2 e i ( 2 ω 1 ω ) T 0 + ( β ˜ 4 A 1 A 2 B 2 + i β ˜ 6 ω A 1 A 2 B 2 + β 10 ρ 2 A 1 A 2 B 2 ) e i ( 2 ω + ω 2 ) T 0 + ( β ˜ 4 A 1 A 2 B ¯ 2 + i β ˜ 6 ω A 1 A 2 B ¯ 2 + β 10 ρ 2 A 1 A 2 B ¯ 2 ) e i ( 2 ω ω 2 ) T 0 + ( β ˜ 4 A 1 A ¯ 2 B 2 + β ˜ 4 A ¯ 1 A 2 B 2 i β ˜ 6 ω A 1 A ¯ 2 B 2 + i β ˜ 6 ω A ¯ 1 A 2 B 2 + β 10 ρ 2 A ¯ 1 A 2 B 2 ) e i ω 2 T 0 + β ˜ 5 A 1 B 2 2 e i ( ω + 2 ω 2 ) T 0 + β ˜ 5 A ¯ 1 B 2 2 e i ( 2 ω 2 ω ) T 0 + 1 2 Ω 2 f ˜ e i Ω T 0 + c c
( D 0 2 + ω 2 ) z 21 = ( 2 i ω D 1 A 2 2 i μ ˜ ω A 2 + 3 α ˜ 1 A 2 2 A ¯ 2 + 2 α ˜ 2 A 2 A 1 A ¯ 1 + α ˜ 2 A ¯ 2 A 1 2 + i α ˜ 3 ω A 2 2 A ¯ 2 + 2 i α ˜ 4 ω A 2 A 1 A ¯ 1 i α ˜ 4 ω A ¯ 2 A 1 2 + 2 α ˜ 5 ω 2 A 2 A 1 A ¯ 1 α ˜ 5 ω 2 A ¯ 2 A 1 2 + α ˜ 6 ω 2 A 2 2 A ¯ 2 + i α ˜ 7 ω A ¯ 2 A 1 2 + 2 γ ˜ 3 A 2 B 2 B ¯ 2 + 2 γ ˜ 5 A 2 B 1 B ¯ 1 + 2 γ 11 ρ 2 A 2 2 A ¯ 2 + γ 12 ρ 1 2 A ¯ 2 A 1 2 + γ 13 ρ 1 A 2 A 1 A ¯ 1 + γ 13 ρ 1 A ¯ 2 A 1 2 i γ 14 ω ρ 1 A 2 A 1 A ¯ 1 + i γ 14 ω ρ 1 A ¯ 2 A 1 2 + γ 15 ρ 2 2 A 2 2 A ¯ 2 + 2 γ 16 ρ 2 A 2 A 1 A ¯ 1 + η ˜ 4 ρ 2 A 2 ) e i ω T 0 + η ˜ 3 B 2 e i ω 2 T 0 + ( α ˜ 1 A 2 3 + α ˜ 2 A 2 A 1 2 + i α ˜ 3 ω A 2 3 + i α ˜ 4 ω A 2 A 1 2 α ˜ 5 ω 2 A 2 A 1 2 α ˜ 6 ω 2 A 2 3 + i α ˜ 7 ω A 2 A 1 2 + i γ 9 ω ρ 2 A 2 3 + γ 11 ρ 2 A 2 3 + γ 12 ρ 1 2 A 2 A 1 2 + γ 13 ρ 1 A 2 A 1 2 + i γ 14 ω ρ 1 A 2 A 1 2 + γ 15 ρ 2 2 A 2 3 + γ 16 ρ 2 A 2 A 1 2 ) e 3 i ω T 0 + ( γ ˜ 1 A 2 2 B 2 + i γ ˜ 2 ω A 2 2 B 2 + γ ˜ 7 A 1 2 B 2 + γ 8 ρ 2 A 2 2 B 2 ) e i ( 2 ω + ω 2 ) T 0 + ( γ ˜ 1 A 2 2 B ¯ 2 + i γ ˜ 2 ω A 2 2 B ¯ 2 + γ ˜ 7 A 1 2 B ¯ 2 + γ 8 ρ 2 A 2 2 B ¯ 2 ) e i ( 2 ω ω 2 ) T 0 + ( 2 γ ˜ 1 A 2 A ¯ 2 B 2 + 2 γ ˜ 7 A 1 A ¯ 1 B 2 + γ 8 ρ 2 A 2 A ¯ 2 B 2 ) e i ω 2 T 0 + γ ˜ 3 A 2 B 2 2 e i ( ω + 2 ω 2 ) T 0 + γ ˜ 3 A ¯ 2 B 2 2 e i ( 2 ω 2 ω ) T 0 + ( γ ˜ 4 A 2 A 1 B 1 + i γ ˜ 6 ω A 2 A 1 B 1 + γ 10 ρ 1 A 2 A 1 B 1 ) e i ( 2 ω + ω 1 ) T 0 + ( γ ˜ 4 A 2 A 1 B ¯ 1 + i γ ˜ 6 ω A 2 A 1 B ¯ 1 + γ 10 ρ 1 A ¯ 2 A 1 B ¯ 1 ) e i ( 2 ω ω 1 ) T 0 + ( γ ˜ 4 A 2 A ¯ 1 B 1 + γ ˜ 4 A ¯ 2 A 1 B 1 i γ ˜ 6 ω A 2 A ¯ 1 B 1 + i γ ˜ 6 ω A ¯ 2 A 1 B 1 + γ 10 ρ 1 A 2 A 1 B 1 ) e i ω 1 T 0 + γ ˜ 5 A 2 B 1 2 e i ( ω + 2 ω 1 ) T 0 + γ ˜ 5 A ¯ 2 B 1 2 e i ( 2 ω 1 ω ) T 0 1 2 i Ω 2 f ˜ e i Ω T 0 + c c
( D 0 2 + ω 1 2 ) z 31 = 2 i ω 1 D 1 B 1 e i ω 1 T 0 2 i μ ˜ 1 ω 1 B 1 e i ω 1 T 0 + η ˜ 5 A 1 e i ω T 0 + c c
( D 0 2 + ω 2 2 ) z 41 = 2 i ω 2 D 1 B 2 e i ω 2 T 0 2 i μ ˜ 2 ω 2 B 2 e i ω 2 T 0 + η ˜ 6 A 2 e i ω T 0 + c c
where c c in Equations (49) to (52) denote the complex conjugate term. To obtain the periodic solutions of Equations (49) to (52), the resonance conditions should be eliminated. Therefore, let σ ,   σ 1 , and σ 2 represent the closeness of the rotor angular speed ( Ω ) and the controller natural frequencies ( ω 1 and ω 2 ) to the rotor system natural frequency ( ω ), as follows:
Ω = ω + σ , ω 1 = ω + σ 1 , ω 2 = ω + σ 2
Inserting Equation (53) into Equations (49) to (52), one can extract the following solvability conditions:
2 i ω D 1 A 1 2 i μ ˜ ω A 1 + 3 α ˜ 1 A 1 2 A ¯ 1 + 2 α ˜ 2 A 1 A 2 A ¯ 2 + α ˜ 2 A ¯ 1 A 2 2 + i α ˜ 3 ω A 1 2 A ¯ 1 + 2 i α ˜ 4 ω A 1 A 2 A ¯ 2 i α ˜ 4 ω A ¯ 1 A 2 2 + 2 α ˜ 5 ω 2 A 1 A 2 A ¯ 2 α ˜ 5 ω 2 A ¯ 1 A 2 2 + α ˜ 6 ω 2 A 1 2 A ¯ 1 + i α ˜ 7 ω A ¯ 1 A 2 2 + 2 β ˜ 3 A 1 B 1 B ¯ 1 + 2 β ˜ 5 A 1 B 2 B ¯ 2 + 2 β 11 ρ 1 A 1 2 A ¯ 1 + β 12 ρ 2 2 A ¯ 1 A 2 2 + β 13 ρ 2 A 1 A 2 A ¯ 2 + β 13 ρ 2 A ¯ 1 A 2 2 i β 14 ω ρ 2 A 1 A 2 A ¯ 2 + i β 14 ω ρ 2 A ¯ 1 A 2 2 + β 15 ρ 1 2 A 1 2 A ¯ 1 + 2 β 16 ρ 1 A 1 A 2 A ¯ 2 + η ˜ 1 B 1 e i ε σ ˜ 1 T 0 + η ˜ 2 ρ 1 A 1 + ( β ˜ 1 A 1 2 B ¯ 1 + i β ˜ 2 ω A 1 2 B ¯ 1 + β ˜ 7 A 2 2 B ¯ 1 + β 8 ρ 1 A 1 2 B ¯ 1 ) e i ε σ ˜ 1 T 0 + ( 2 β ˜ 1 A 1 A ¯ 1 B 1 + 2 β ˜ 7 A 2 A ¯ 2 B 1 + β 8 ρ 1 A 1 A ¯ 1 B 1 ) e i ε σ ˜ 1 T 0 + ( β ˜ 4 A 1 A 2 B ¯ 2 + i β ˜ 6 ω A 1 A 2 B ¯ 2 + β 10 ρ 2 A 1 A 2 B ¯ 2 ) e i ε σ ˜ 2 T 0 + β ˜ 3 A ¯ 1 B 1 2 e 2 i ε σ ˜ 1 T 0 + ( β ˜ 4 A 1 A ¯ 2 B 2 + β ˜ 4 A ¯ 1 A 2 B 2 i β ˜ 6 ω A 1 A ¯ 2 B 2 + i β ˜ 6 ω A ¯ 1 A 2 B 2 + β 10 ρ 2 A ¯ 1 A 2 B 2 ) e i ε σ ˜ 2 T 0 + β ˜ 5 A ¯ 1 B 2 2 e 2 i ε σ ˜ 2 T 0 + 1 2 ( ω + ε σ ˜ ) 2 f ˜ e i ε σ ˜ T 0 = 0
2 i ω D 1 A 2 2 i μ ˜ ω A 2 + 3 α ˜ 1 A 2 2 A ¯ 2 + 2 α ˜ 2 A 2 A 1 A ¯ 1 + α ˜ 2 A ¯ 2 A 1 2 + i α ˜ 3 ω A 2 2 A ¯ 2 + 2 i α ˜ 4 ω A 2 A 1 A ¯ 1 i α ˜ 4 ω A ¯ 2 A 1 2 + 2 α ˜ 5 ω 2 A 2 A 1 A ¯ 1 α ˜ 5 ω 2 A ¯ 2 A 1 2 + α ˜ 6 ω 2 A 2 2 A ¯ 2 + i α ˜ 7 ω A ¯ 2 A 1 2 + 2 γ ˜ 3 A 2 B 2 B ¯ 2 + 2 γ ˜ 5 A 2 B 1 B ¯ 1 + 2 γ 11 ρ 2 A 2 2 A ¯ 2 + γ 12 ρ 1 2 A ¯ 2 A 1 2 + γ 13 ρ 1 A 2 A 1 A ¯ 1 + γ 13 ρ 1 A ¯ 2 A 1 2 i γ 14 ω ρ 1 A 2 A 1 A ¯ 1 + i γ 14 ω ρ 1 A ¯ 2 A 1 2 + γ 15 ρ 2 2 A 2 2 A ¯ 2 + 2 γ 16 ρ 2 A 2 A 1 A ¯ 1 + η ˜ 3 B 2 e i ε σ ˜ 2 T 0 + η ˜ 4 ρ 2 A 2 + ( γ ˜ 1 A 2 2 B ¯ 2 + i γ ˜ 2 ω A 2 2 B ¯ 2 + γ ˜ 7 A 1 2 B ¯ 2 + γ 8 ρ 2 A 2 2 B ¯ 2 ) e i ε σ ˜ 2 T 0 + ( 2 γ ˜ 1 A 2 A ¯ 2 B 2 + 2 γ ˜ 7 A 1 A ¯ 1 B 2 + γ 8 ρ 2 A 2 A ¯ 2 B 2 ) e i ε σ ˜ 2 T 0 + ( γ ˜ 4 A 2 A 1 B ¯ 1 + i γ ˜ 6 ω A 2 A 1 B ¯ 1 + γ 10 ρ 1 A 2 A 1 B ¯ 1 ) e i ε σ ˜ 1 T 0 + γ ˜ 3 A ¯ 2 B 2 2 e 2 i ε σ ˜ 2 T 0 + ( γ ˜ 4 A 2 A ¯ 1 B 1 + γ ˜ 4 A ¯ 2 A 1 B 1 i γ ˜ 6 ω A 2 A ¯ 1 B 1 + i γ ˜ 6 ω A ¯ 2 A 1 B 1 + γ 10 ρ 1 A ¯ 2 A 1 B 1 ) e i ε σ ˜ 1 T 0 + γ ˜ 5 A ¯ 2 B 1 2 e 2 i ε σ ˜ 1 T 0 1 2 i ( ω + ε σ ˜ ) 2 f ˜ e i ε σ ˜ T 0 = 0
2 i ( ω + ε σ ˜ 1 ) D 1 B 1 e i ε σ ˜ 1 T 0 2 i μ ˜ 1 ( ω + ε σ ˜ 1 ) B 1 e i ε σ ˜ 1 T 0 + η ˜ 5 A 1 = 0
2 i ( ω + ε σ ˜ 2 ) D 1 B 2 e i ε σ ˜ 2 T 0 2 i μ ˜ 2 ( ω + ε σ ˜ 2 ) B 2 e i ε σ ˜ 1 T 0 + η ˜ 6 A 2 = 0
To obtain the autonomous dynamical system that describes the oscillatory behaviors of the considered closed-loop system, let us express the unknown functions A 1 , A 2 ,   B 1 , and B 2 in the polar form as follows:
A 1 ( T 1 ) = 1 2 a 1 ( T 1 ) e i θ 1 ( T 1 ) , A 2 ( T 1 ) = 1 2 a 2 ( T 1 ) e i θ 2 ( T 1 ) B 1 ( T 1 ) = 1 2 b 1 ( T 1 ) e i θ 3 ( T 1 ) , B 2 ( T 1 ) = 1 2 b 2 ( T 1 ) e i θ 4 ( T 1 ) }
According to Equation (58), we have the following:
D 1 A 1 = d ε d t A 1 = 1 2 ε ( a ˙ 1 e i θ 1 + i a 1 θ ˙ 1 e i θ 1 ) , D 1 A 2 = d ε d t A 2 = 1 2 ε ( a ˙ 2 e i θ 2 + i a 2 θ ˙ 2 e i θ 2 ) D 1 B 1 = d ε d t B 1 = 1 2 ε ( b ˙ 1 e i θ 3 + i b 1 θ ˙ 3 e i θ 3 ) , D 1 B 2 = d ε d t B 2 = 1 2 ε ( b ˙ 2 e i θ 4 + i b 2 θ ˙ 4 e i θ 4 ) }
Substituting Equations (58) and (59) into Equations (54) to (57) and separating the real and imaginary part yields the following:
a ˙ 1 = F 1 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = 1 2 ( 2 μ + η 2 η 7 ω 3 2 + ω 2 ) a 1 + 1 8 ( α 3 2 β 11 η 7 ω 3 2 + ω 2 2 ω 3 β 15 η 7 2 ( ω 3 2 + ω 2 ) 2 ) a 1 3 + 1 8 ( 2 α 4 β 13 η 8 ω 4 2 + ω 2 ω 4 β 14 η 8 ω 4 2 + ω 2 2 β 16 η 7 ω 3 2 + ω 2 ) a 1 a 2 2 + 1 8 ( α 4 + α 7 + 2 ω 4 β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 β 13 η 8 ω 4 2 + ω 2 + ω 4 β 14 η 8 ω 4 2 + ω 2 ) a 1 a 2 2 cos ( 2 ϕ 1 2 ϕ 2 ) + 1 8 ( α 2 ω α 5 ω + ( ω 4 2 ω 2 ) β 12 η 8 2 ω ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω ( ω 4 2 + ω 2 ) ω β 14 η 8 ω 4 2 + ω 2 ) a 1 a 2 2 sin ( 2 ϕ 1 2 ϕ 2 ) + ( 1 2 ω η 1 b 1 1 8 ω β 1 a 1 2 b 1 1 4 ω β 7 a 2 2 b 1 ) sin ( ϕ 3 ) 1 8 ω β 3 a 1 b 1 2 sin ( 2 ϕ 3 ) + 1 8 ( β 2 2 β 8 η 7 ω 3 2 + ω 2 ) a 1 2 b 1 cos ( ϕ 3 ) 1 8 ( β 10 η 8 ω 4 2 + ω 2 ) a 1 a 2 b 2 cos ( ϕ 4 ) + 1 8 ω ( ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 1 a 2 b 2 sin ( ϕ 4 ) + 1 8 ω β 7 a 2 2 b 1 sin ( 2 ϕ 1 2 ϕ 2 + ϕ 3 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 1 a 2 b 2 cos ( 2 ϕ 1 2 ϕ 2 ϕ 4 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 1 a 2 b 2 sin ( 2 ϕ 1 2 ϕ 2 ϕ 4 ) + 1 8 ω β 5 a 1 b 2 2 sin ( 2 ϕ 1 2 ϕ 2 2 ϕ 4 ) + 1 2 ω ( ω + σ ) 2 f sin ( ϕ 1 )
a ˙ 2 = F 2 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = 1 2 ( 2 μ + η 4 η 8 ω 4 2 + ω 2 ) a 2 + 1 8 ( α 3 2 γ 11 η 8 ω 4 2 + ω 2 2 ω 4 γ 15 η 8 2 ( ω 4 2 + ω 2 ) 2 ) a 2 3 + 1 8 ( 2 α 4 γ 13 η 7 ω 3 2 + ω 2 ω 3 γ 14 η 7 ω 3 2 + ω 2 2 γ 16 η 8 ω 4 2 + ω 2 ) a 2 a 1 2 + 1 8 ( α 4 + α 7 2 ω 3 γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 γ 13 η 7 ω 3 2 + ω 2 + ω 3 γ 14 η 7 ω 3 2 + ω 2 ) a 2 a 1 2 cos ( 2 ϕ 2 2 ϕ 1 ) + 1 8 ( α 2 ω α 5 ω + ( ω 3 2 ω 2 ) γ 12 η 7 2 ω ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω ( ω 3 2 + ω 2 ) + ω γ 14 η 7 ω 3 2 + ω 2 ) a 2 a 1 2 sin ( 2 ϕ 2 2 ϕ 1 ) + ( 1 2 ω η 3 b 2 1 8 ω γ 1 a 2 2 b 2 1 4 ω γ 7 a 1 2 b 2 ) sin ( ϕ 4 ) 1 8 ω γ 3 a 2 b 2 2 sin ( 2 ϕ 4 ) + 1 8 ω ( γ 2 ω 2 ω γ 8 η 8 ω 4 2 + ω 2 ) a 2 2 b 2 cos ( ϕ 4 ) 1 8 ω ( ω γ 10 η 7 ω 3 2 + ω 2 ) a 2 a 1 b 1 cos ( ϕ 3 ) + 1 8 ω ( ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 2 a 1 b 1 sin ( ϕ 3 ) + 1 8 ω γ 7 a 1 2 b 2 sin ( 2 ϕ 2 2 ϕ 1 + ϕ 4 ) + 1 8 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 2 a 1 b 1 cos ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 2 a 1 b 1 sin ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) + 1 8 ω γ 5 a 2 b 1 2 sin ( 2 ϕ 2 2 ϕ 1 2 ϕ 3 ) 1 2 ω ( ω + σ ) 2 f cos ( ϕ 2 )
b ˙ 1 = F 3 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = μ 1 b 1 + 1 2 ( ω + σ 1 ) η 5 a 1 sin ( ϕ 3 )
b ˙ 2 = F 4 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = μ 2 b 2 + 1 2 ( ω + σ 2 ) η 6 a 2 sin ( ϕ 4 )
ϕ ˙ 1 = F 5 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = σ + 1 2 ω ( ω 3 η 2 η 7 ω 3 2 + ω 2 ) + 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ω 3 β 11 η 7 ω 3 2 + ω 2 + ( ω 3 2 ω 2 ) β 15 η 7 2 ( ω 3 2 + ω 2 ) 2 ) a 1 2 + 1 8 ω ( 2 α 2 + 2 α 5 ω 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 ω 2 β 14 η 8 ω 4 2 + ω 2 + 2 ω 3 β 16 η 7 ω 3 2 + ω 2 ) a 2 2 + 1 4 ω β 3 b 1 2 + 1 4 ω β 5 b 2 2 + 1 8 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 2 2 cos ( 2 ϕ 1 2 ϕ 2 ) + 1 8 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 2 2 sin ( 2 ϕ 1 2 ϕ 2 ) + ( 1 2 η 1 b 1 + 3 8 β 1 a 1 2 b 1 + 1 4 β 7 a 2 2 b 1 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 a 1 2 b 1 ) cos ( ϕ 3 ) ω a 1 + 1 8 ω β 3 b 1 2 cos ( 2 ϕ 3 ) 1 8 β 8 η 7 ω 3 2 + ω 2 a 1 b 1 sin ( ϕ 3 ) 1 8 β 2 a 1 b 1 sin ( ϕ 3 ) + 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 cos ( ϕ 4 ) + 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 sin ( ϕ 4 ) + 1 8 ω a 1 β 7 a 2 2 b 1 cos ( 2 ϕ 1 2 ϕ 2 + ϕ 3 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 cos ( 2 ϕ 1 2 ϕ 2 ϕ 4 ) + 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 sin ( 2 ϕ 1 2 ϕ 2 ϕ 4 ) + 1 8 ω β 5 b 2 2 cos ( 2 ϕ 1 2 ϕ 2 2 ϕ 4 ) + 1 2 ω a 1 ( ω + σ ) 2 f cos ( ϕ 1 )
ϕ ˙ 2 = F 6 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = σ + 1 2 ω ( ω 4 η 4 η 8 ω 4 2 + ω 2 ) + 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ω 4 γ 11 η 8 ω 4 2 + ω 2 + ( ω 4 2 ω 2 ) γ 15 η 8 2 ( ω 4 2 + ω 2 ) 2 ) a 2 2 + 1 8 ω ( 2 α 2 + 2 α 5 ω 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 ω 2 γ 14 η 7 ω 3 2 + ω 2 + 2 ω 4 γ 16 η 8 ω 4 2 + ω 2 ) a 1 2 + 1 4 ω γ 3 b 2 2 + 1 4 ω γ 5 b 1 2 + 1 8 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 1 2 cos ( 2 ϕ 2 2 ϕ 1 ) + 1 8 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 1 2 sin ( 2 ϕ 2 2 ϕ 1 ) + ( 1 2 η 3 b 2 + 3 8 γ 1 a 2 2 b 2 + 1 4 γ 7 a 1 2 b 2 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 a 2 2 b 2 ) cos ( ϕ 4 ) ω a 2 + 1 8 ω γ 3 b 2 2 cos ( 2 ϕ 4 ) 1 8 γ 8 η 8 ω 4 2 + ω 2 a 2 b 2 sin ( ϕ 4 ) 1 8 γ 2 a 2 b 2 sin ( ϕ 4 ) + 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 cos ( ϕ 3 ) + 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 sin ( ϕ 3 ) + 1 8 ω a 2 γ 7 a 1 2 b 2 cos ( 2 ϕ 2 2 ϕ 1 + ϕ 4 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 cos ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) + 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 sin ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) + 1 8 ω γ 5 b 1 2 cos ( 2 ϕ 2 2 ϕ 1 2 ϕ 3 ) + 1 2 ω a 2 ( ω + σ ) 2 f sin ( ϕ 2 )
ϕ ˙ 3 = F 7 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = σ 1 + 1 2 ( ω + σ 1 ) b 1 η 5 a 1 cos ( ϕ 3 ) 1 2 ω ( ω 3 η 2 η 7 ω 3 2 + ω 2 ) 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ω 3 β 11 η 7 ω 3 2 + ω 2 + ( ω 3 2 ω 2 ) β 15 η 7 2 ( ω 3 2 + ω 2 ) 2 ) a 1 2 1 8 ω ( 2 α 2 + 2 α 5 ω 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 ω 2 β 14 η 8 ω 4 2 + ω 2 + 2 ω 3 β 16 η 7 ω 3 2 + ω 2 ) a 2 2 1 4 ω β 3 b 1 2 1 4 ω β 5 b 2 2 1 8 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 2 2 cos ( 2 ϕ 1 2 ϕ 2 ) 1 8 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 2 2 sin ( 2 ϕ 1 2 ϕ 2 ) ( 1 2 η 1 b 1 + 3 8 β 1 a 1 2 b 1 + 1 4 β 7 a 2 2 b 1 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 a 1 2 b 1 ) cos ( ϕ 3 ) ω a 1 1 8 ω β 3 b 1 2 cos ( 2 ϕ 3 ) + 1 8 β 8 η 7 ω 3 2 + ω 2 a 1 b 1 sin ( ϕ 3 ) 1 8 β 2 a 1 b 1 sin ( ϕ 3 ) 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 cos ( ϕ 4 ) 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 sin ( ϕ 4 ) 1 8 ω a 1 β 7 a 2 2 b 1 cos ( 2 ϕ 1 2 ϕ 2 + ϕ 3 ) 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 cos ( 2 ϕ 1 2 ϕ 2 ϕ 4 ) 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 2 b 2 sin ( 2 ϕ 1 2 ϕ 2 ϕ 4 ) 1 8 ω β 5 b 2 2 cos ( 2 ϕ 1 2 ϕ 2 2 ϕ 4 ) 1 2 ω a 1 ( ω + σ ) 2 f cos ( ϕ 1 )
ϕ ˙ 4 = F 8 ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = σ 2 + 1 2 ( ω + σ 2 ) b 2 η 6 a 2 cos ( ϕ 4 ) 1 2 ω ( ω 4 η 4 η 8 ω 4 2 + ω 2 ) 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ω 4 γ 11 η 8 ω 4 2 + ω 2 + ( ω 4 2 ω 2 ) γ 15 η 8 2 ( ω 4 2 + ω 2 ) 2 ) a 2 2 1 8 ω ( 2 α 2 + 2 α 5 ω 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 ω 2 γ 14 η 7 ω 3 2 + ω 2 + 2 ω 4 γ 16 η 8 ω 4 2 + ω 2 ) a 1 2 1 4 ω γ 3 b 2 2 1 4 ω γ 5 b 1 2 1 8 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 1 2 cos ( 2 ϕ 2 2 ϕ 1 ) 1 8 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 1 2 sin ( 2 ϕ 2 2 ϕ 1 ) ( 1 2 η 3 b 2 + 3 8 γ 1 a 2 2 b 2 + 1 4 γ 7 a 1 2 b 2 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 a 2 2 b 2 ) cos ( ϕ 4 ) ω a 2 1 8 ω γ 3 b 2 2 cos ( 2 ϕ 4 ) + 1 8 γ 8 η 8 ω 4 2 + ω 2 a 2 b 2 sin ( ϕ 4 ) + 1 8 γ 2 a 2 b 2 sin ( ϕ 4 ) 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 cos ( ϕ 3 ) 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 sin ( ϕ 3 ) 1 8 ω a 2 γ 7 a 1 2 b 2 cos ( 2 ϕ 2 2 ϕ 1 + ϕ 4 ) 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 cos ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 1 b 1 sin ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) 1 8 ω γ 5 b 1 2 cos ( 2 ϕ 2 2 ϕ 1 2 ϕ 3 ) 1 2 ω a 2 ( ω + σ ) 2 f sin ( ϕ 2 )
where ϕ 1 = σ t θ 1 ,   ϕ 2 = σ t θ 2 ,   ϕ 3 = θ 1 θ 3 σ 1 t , and ϕ 4 = θ 2 θ 4 σ 2 t . By inserting Equations (43) to (48) and (58) into Equations (25) to (30), one can extract an approximate solution for the closed-loop system given by Equations (19) to (24), as follows:
z 1 ( t ) = a 1 ( t ) cos ( Ω t ϕ 1 ( t ) )
z 2 ( t ) = a 2 ( t ) cos ( Ω t ϕ 2 ( t ) )
z 3 ( t ) = b 1 ( t ) cos ( Ω t ( ϕ 1 ( t ) + ϕ 3 ( t ) ) )
z 4 ( t ) = b 2 ( t ) cos ( Ω t ( ϕ 2 ( t ) + ϕ 4 ( t ) ) )
z 5 ( t ) = η 7 a 1 ω 3 2 + ω 2 ( ω 3 cos ( Ω t ϕ 1 ( t ) ) + ω sin ( Ω t ϕ 1 ( t ) ) )
z 6 ( t ) = η 8 a 2 ω 4 2 + ω 2 ( ω 4 cos ( Ω t ϕ 2 ( t ) ) + ω sin ( Ω t ϕ 2 ( t ) ) )
It is clear from Equations (68) to (71) that a 1 ( t ) and a 2 ( t ) are the steady-state oscillation amplitudes of the twelve-poles rotor system, while ϕ 1 ( t ) and ϕ 2 ( t ) represent the phase angles of the controlled rotor. In addition, b 1 ( t ) and b 2 ( t ) represent the oscillation amplitudes of the PPF controllers and ϕ 1 ( t ) + ϕ 3 ( t ) ,   ϕ 2 ( t ) + ϕ 4 ( t ) are the corresponding phase angles. In addition, Equations (72) and (73) show that the dynamical characteristics of the IRC controller depend on the dynamics of the rotor system (i.e., z 5 ( t ) depends on a 1 ( t ) ,   ϕ 1 ( t ) and z 6 ( t ) depends on a 2 ( t ) ,   ϕ 2 ( t ) ) . Moreover, the derived nonlinear autonomous system that is provided by Equations (60) to (67) governs the evolution of the oscillation amplitudes ( a 1 ,   a 2 ,   b 1 ,   b 2 ) and the corresponding phase angles ( ϕ 1 ,   ϕ 2 ,   ϕ 3 , ϕ 4 ) of the closed-loop system as a function of the different system parameters (i.e., f ,   σ ,   σ 1 ,   σ 2 ,   p ,   d ,   η 1 ,   η 2 ,   η 3 ,   η 4 ,   η 5 ,   η 6 ,   η 7 ,   η 8 , , etc . ). Accordingly, the dynamical characteristics of the closed-loop system can be explored by investigating the nonlinear dynamical system provided by Equations (60) to (67). Therefore, one can explore the steady-state dynamics of the closed-loop system by inserting a ˙ 1 = a ˙ 2 = b ˙ 1 = b ˙ 4 = ϕ ˙ 1 = ϕ ˙ 2 = ϕ ˙ 3 = ϕ ˙ 4 = 0 into Equations (60) to (67), which results in the following nonlinear algebraic system:
F j ( a 1 , a 2 , b 1 , b 2 , ϕ 1 , ϕ 2 , ϕ 3 , ϕ 4 ) = 0 ; j = 1 , 2 , , 8
Solving Equation (74), utilizing σ as a bifurcation parameter at the different values of the system and control parameters ( f ,   σ 1 ,   σ 2 ,   p ,   d ,   η 1 ,   η 2 ,   η 3 ,   η 4 ,   η 5 ,   η 6 ,   η 7 ,   η 8 ,   ,   etc . ), we can explore the efficiency of the introduced control technique (i.e., P D + I R C + P P F controller). In addition, to investigate the stability of the solution of Equation (74), one can check the eigenvalues of the Jacobian matrix of the dynamical system given by Equations (60) to (67), which can be obtained via letting ( a 10 ,   a 20 ,   b 10 ,   b 20 ,   ϕ 10 ,   ϕ 20 ,   ϕ 30 ,   ϕ 40 ) be the solution of Equation (74) and ( a 11 ,   a 21 ,   b 10 ,   b 20 ,   ϕ 11 ,   ϕ 21 ,   ϕ 31 ,   ϕ 41 ) be a small deviation about this solution. Therefore, one can write
a j = a j 0 + a j 1 , b j = b j 0 + b j 1 , ϕ k = ϕ k 0 + ϕ k 1 , a ˙ j = a ˙ j 1 , b ˙ j = b ˙ j 1 , ϕ ˙ k = ϕ ˙ k 1 ; j = 1 , 2 ; k = 1 , 2 , , 4 .
Inserting Equation (75) into Equations (60) to (67) and expanding for the small deviations ( a 11 ,   a 21 ,   b 10 ,   b 20 ,   ϕ 11 ,   ϕ 21 ,   ϕ 31 ,   ϕ 41 ), retaining the linear terms only, one can derive the following linearized dynamical system:
a ˙ 11 a ˙ 21 b ˙ 11 b ˙ 21 ϕ ˙ 11 ϕ ˙ 21 ϕ ˙ 31 ϕ ˙ 41 = J 11 J 12 J 13 J 14 J 15 J 16 J 17 J 18 J 21 J 22 J 23 J 24 J 25 J 26 J 27 J 28 J 31 J 32 J 33 J 34 J 35 J 36 J 37 J 38 J 41 J 42 J 43 J 44 J 45 J 46 J 47 J 48 J 51 J 52 J 53 J 54 J 55 J 56 J 57 J 58 J 61 J 62 J 63 J 64 J 65 J 66 J 67 J 68 J 71 J 72 J 73 J 74 J 75 J 76 J 77 J 78 J 81 J 82 J 83 J 84 J 85 J 86 J 87 J 88 a 11 a 21 b 11 b 21 ϕ 11 ϕ 21 ϕ 31 ϕ 41
where J m n ( m = 1 , 2 , , 8 ,   n = 1 , 2 , , 8 ) are provided in Appendix C. Accordingly, the stability of the dynamical system provided by Equations (60) to (67) has been studied by examining the eigenvalues of the linearized system provided by Equation (76) (see [42]), where the stable solution was illustrated as a solid-line, while the unstable solution was plotted as a dotted-line, as shown in the different bifurcation diagrams in Section 4.

4. Steady-State Oscillation and Bifurcation Analysis

Based on both the derived mathematical model of the closed-loop system provided by Equations (19) to (24) and the nonlinear algebraic Equation (74), one can investigate the efficiency of the proposed control technique (i.e., P D + I R C + P P F controller) in improving the oscillatory characteristics and eliminating the catastrophic bifurcation behaviors of the studied twelve-pole system. As Equation (74) governs the steady-state vibration amplitudes ( a 1 , a 2 , b 1 , b 2 ) and the corresponding phase angles (i.e., ϕ 1 ,   ϕ 2 ,   ϕ 3 ,   ϕ 4 ), we can investigate the steady-state oscillatory behaviors of both the rotor system and the connected controller via solving Equation (74) numerically using the Newton–Raphson predictor–corrector algorithm (see [43]), utilizing σ as the bifurcation parameter, where the stable solution is plotted as a solid-line and the unstable solution is shown as a dotted-line. In addition, to validate the accuracy of the derived analytical solution (i.e., Equation (74)), as well as to investigate the full system response (i.e., steady-state and transient response) of the closed-loop system, one can simulate the system’s temporal equations of motion (i.e., Equations (19) to (24)) numerically using the Rung–Kutta method of fourth order. Accordingly, the following values of the parameters have been used to simulate the system dynamics [29,36]: f = 0.013 , p = 1.5 , d = 0.005 , η 1 = η 2 = η 3 = η 4 = η 5 = η 6 = 0.2 , η 7 = η 8 = 1 , μ 1 = μ 2 = 0.01 , σ = σ 1 = σ 2 = 0 , Ω = ω + σ , ω 1 = ω + σ 1 , ω 2 = ω + σ 2 , ω 3 = ω 4 = 1 , α = 30 ° , where the other system parameters μ , ω , α j , β k , γ k , ( j = 1 , 2 , , 7 ;   k = 1 , 2 , , 16 ) are defined below Equation (24). Before proceeding further, let us go back first to the normalized equations of motion (Equations (19) to (24)), where the rotor normalized temporal displacements in the X and Y directions are defined such that z 1 ( t ) = x ( t ) c 0 and z 2 ( t ) = y ( t ) c 0 , where c 0 is the nominal air-gap size between the rotor and the poles-housing and x ( t ) , y ( t ) are the actual temporal displacement of the rotor in the X and Y directions, respectively. Accordingly, for safe working conditions for the rotor system without the occurrence of rub and/or impact forces between the rotor and the pole housing, x ( t ) and y ( t ) should be smaller than the air-gap size c 0 (i.e., | x ( t ) c 0 | = | z 1 ( t ) | < 1 and | y ( t ) c 0 | = | z 2 ( t ) | < 1 ). Therefore, for the safe operation of the rotor system without rub and/or impact between the rotor and the stator, | z 1 ( t ) | = | a 1 ( t ) cos ( Ω t ϕ 1 ( t ) ) | and | z 2 ( t ) | = | a 2 ( t ) cos ( Ω t ϕ 2 ( t ) ) | should be smaller than unity, which implies that | a 1 | and | a 2 | must be lower than unity (i.e., | a 1 | < 1 &   | a 2 | < 1 ). In addition, the parameter σ is defined in Equation (53) such that Ω = ω + σ . Accordingly, σ is used in the whole article as a bifurcation control parameter to describe the rotor dynamics when the system angular speed ( Ω ) is close to or equal to the rotor’s natural frequency ( ω ).

4.1. System Dynamics in the Case of P D -Control Algorithm

The parameters P = c 0 I 0 k 1 and d = c 0 ϑ I 0 k 2 denote the normalized proportional gain and derivative gain of the P D -control algorithm, respectively. In addition, the parameters η 1 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 3 and η 3 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 5 are the normalized control gains the P P F -control algorithm that is connected to the rotor system, while η 2 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 4 and η 4 = ( 1 + 2 cos ( α ) ) c 0 I 0 k 6 represent the normalized control gains of the I R C -control algorithm (see Equations (6), (19) and (20)). Moreover, the parameters η 5 = L 1 ϑ 2 and η 6 = L 2 ϑ 2 are the normalized feedback gains of the P P F -control algorithm, while η 7 = L 3 ϑ , η 8 = L 4 ϑ denote the feedback gains of the I R C -control algorithm (see Equations (7)–(10) and (21)–(24)). Accordingly, one can investigate the influence of the P D controller only on the rotor dynamics via setting η k = 0 ( k = 1 , 2 , , 4 ).
This section is dedicated to investigating the rotor dynamics in the case of the P D -control algorithm only, at different levels of the excitation force f , as shown in Figure 3. The figure was obtained by solving the nonlinear system provided by Equation (74) when η k = 0 ,   ( k = 1 , 2 , , 8 ) at f = 0.004 ,   0.007 ,   0.01 , and 0.013 . Figure 3a,b shows the rotor steady-state vibration amplitudes in both the X and Y directions at four different values of the excitation force f , while Figure 3c shows the evolution of the phase angles ϕ 1 and ϕ 2 versus σ when f = 0.004 . In addition, Figure 3d illustrates the phase angles ϕ 1 and ϕ 2 at f = 0.013 . It is clear from Figure 3a,b that the vibration amplitudes ( a 1 and a 2 ) of the twelve-poles system is a monotonic increasing function of the excitation force, where the rotor system may be subjected to rub and/or impact force between the rotating disk and the pole-housing if f 0.013 (i.e., the rotor may exhibit vibration amplitudes a 1 > 1 and/or a 2 > 1 if f 0.013 ). Accordingly, one can conclude that the considered system can work properly without a catastrophic rub and/or impact between the rotor and stator, as long as the excitation force f is smaller than 0.013 when only the P D -control algorithm is applied. In addition, Figure 3c,d depicts that the phase angle ϕ 2 is always greater than ϕ 1 , which means that the rotor system performs a forward whirling motion only (according to Equations (68) and (69)) along the σ axis, regardless of the excitation force magnitude. By examining Figure 3, one can note that the rotor system has symmetric oscillation amplitudes in both the X and Y directions (i.e., a 1 = a 2 ) and the phase difference ϕ 2 ϕ 1 is always π 2 , which demonstrates that the rotor system performs a circular forward whirling motion along the σ axis, regardless of the excitation force magnitude.

4.2. System Dynamics in the Case of the P D + P P F -Control Algorithm

The rotor dynamics at four different magnitudes of the excitation force were investigated when both the P D - and P P F -control algorithms were applied simultaneously. Figure 4a–c shows the steady-state vibration amplitudes of both the rotor system ( a 1 and a 2 ) and the P P F -control algorithm ( b 1 and b 2 ) when p = 1.5 , d = 0.005 , η 1 = η 3 = 0.2 , and η 2 = η 4 = 0.0 at the four excitation force amplitudes f = 0.0025 ,   0.0075 ,   0.0125 , and 0.0175 . In addition, Figure 4e,f illustrates the evolution of the rotor phase angles ϕ 1 and ϕ 2 at f = 0.0025 and f = 0.0175 , respectively. By comparing Figure 3a,b with Figure 4a,b, one can deduce that the integration of the PPF-control algorithm to the twelve-poles rotor has suppressed the system’s vibrations at the perfect resonance condition (i.e., it has suppressed the system’s vibrations at σ = 0.0 ). However, two undesired resonant peaks appeared on both sides of σ = 0.0 . In addition, Figure 4a,b demonstrates that the rotor system may work safely without rub and/or impact between the rotor and the poles-housing, as long as the excitation force f < 0.0175 (i.e., a 1 < 1 and a 2 < 1 as long as f < 0.0175 ). Moreover, Figure 4e,f shows that the phase difference ( ϕ 2 ϕ 1 ) of the rotor lateral oscillations in the X and Y directions is always constant, so that ϕ 2 ϕ 1 = π 2 , which implies that the system exhibits only a forward circular whirling motion, regardless of both the angular speed and the excitation force magnitude. Generally, Figure 4 shows that the integration of the P P F -control algorithm to the system with the P -controller suppressed the rotor’s undesired vibrations at the perfect resonance condition (i.e., when Ω = ω + σ , σ = 0.0 ), regardless of the excitation force amplitude; however, the system may suffer from high oscillation, especially if Ω > ω . Accordingly, the P P F -control algorithm acts as a notch filter that eliminates the system’s vibrations at a specific frequency band.

4.3. System Dynamics in the Case of the P D + I R C -Control Algorithm

The oscillatory behaviors of the system were explored when the I R C -control algorithm was coupled to the rotor system with the P D controller, while the P P F controller was turned off. Accordingly, Figure 5 shows the motion bifurcation of the rotor system when p = 1.5 , d = 0.005 , η 1 = η 3 = 0 , and η 2 = η 4 = 0.2 at four different magnitudes of the excitation force (i.e., f = 0.02 ,   0.04 ,   0.06 and 0.08 ). It is clear from Figure 5a,b that the I R C -control algorithm forced the twelve-poles system to respond like the linear system, even at the strong excitation forces. Moreover, Figure 5c,d demonstrates that the system can perform only a circular forward whirling motion along the σ axis, regardless of the excitation force magnitude, where ϕ 2 ϕ 1 = π 2 and a 1 , a 2 are symmetric on the interval 0.3 < σ < 0.3 . By examining Figure 5, we can deduce that the system can rotate safely without rub and/or impact force between the rotor and the stator, even at the strong excitation forces (i.e., f = 0.08 ), compared with the case of the P D -control algorithm only, as shown in Figure 3. Therefore, coupling the I R C -control algorithm to the system increased the rotor linear damping coefficients, which ultimately decreased the lateral vibrations even at the large excitation forces. However, the IRC controller could not eliminate the system vibrations close to zero at the resonance condition (i.e., at σ = 0 ), as in the case of the PPF-control technique, but the maximum vibration occurred at σ = 0 . Therefore, utilizing the P P F - and I R C -control techniques as a one-control algorithm, along with the PD controller, may have the advantages of both the P P F and I R C controllers, as illustrated in the next subsection.

4.4. System Dynamics in the Case of P D + I R C + P P F -Control Algorithm

The dynamical behaviors of the considered twelve-poles rotor system were explored when the three control algorithms (i.e., P D + I R C + P P F -control algorithms) were activated simultaneously. Figure 6 shows the nonlinear dynamics of the controlled rotor system when P = 1.5 , d = 0.005 , η 1 = η 2 = η 3 = η 4 = 0.2 when the extinction force f = 0.025 ,   0.045 ,   0.065 , and 0.085 . Figure 6a,b illustrates the evolution of the system’s lateral vibrations ( a 1 , a 2 ) against σ , while Figure 6c,d shows the vibration amplitudes of the PPF controller against the detuning parameter σ . In addition, Figure 6e,f illustrates the phase difference of the rotor’s lateral vibrations in both the X and Y directions when f = 0.025 and f = 0.085 , respectively. It is clear from Figure 6a,b that the vibration amplitudes ( a 1 , a 2 ) of the twelve-poles system was close to zero, regardless of the excitation force magnitude, as long as σ 0 , due to the effect of the P P F -control algorithm. In addition, the resonant peaks that appeared on both sides of σ = 0 (as in Figure 4a,b) were mitigated, due to the effect of the I R C -control algorithm. Moreover, Figure 6a,b shows that the rotor system worked properly without impact occurrence between the rotor and stator, as long as f 0.085 . It was also clear from Figure 6e,f that the controlled rotor system performed a circular forward whirling motion, as long as 0.3 σ 0.3 , where the phase difference was ϕ 2 ϕ 1 = π / 2 . Based on Figure 3 to Figure 6, one can conclude that the integration of the P D , I R C , and P P F -control algorithms to act as a single controller can provide for the safe operation of the considered rotor system with small oscillation amplitudes at the resonant conditions (i.e., when σ = 0.0 ), even if the excitation force is strong.

4.5. Sensitivity Analysis of the P D + I R C + P P F -Control Algorithm

As the combined control algorithm (i.e., P D + I R C + P P F ) has many advantages over the individual three control techniques, this subsection explores the sensitivity of this control method to the variation of different control gains. The effect of increasing the proportional gain ( P ) on the vibration suppression efficiency of the control algorithm is illustrated in Figure 7. The figure shows that the increase in P increases the oscillation amplitudes ( a 1 , a 2 ) of the twelve-poles system and degrades the control algorithm’s efficiency. Therefore, the P gain should be kept at the small possible value to guarantee the high performance of the proposed control technique. Based on the system parameters provided below Equation (24), the natural frequency of the rotor system ω is defined as = 2 P cos ( α ) + P 3 . Therefore, the minimum value of P should be selected in a way that guarantees that ω > 0 . On the other hand, the effect of the P P F -control gains (i.e., η 1 and η 3 ) on the whole-system dynamics is depicted in Figure 8. The figure demonstrates that the increase of η 1 and η 3 (i.e., η 1 = η 3 = 0.5 ) enhanced the controller performance in eliminating the rotor oscillations at the perfect resonance condition (i.e., when σ = 0.0 ), as well as widening the frequency band at which the system could work properly with small vibration amplitudes. In addition, Figure 9 demonstrates that the increase in the I R C -control gains (i.e., η 2 = η 4 = 0.5 ) decreased the resonant peaks that appeared on both sides of σ = 0.0 , and improved the controller efficiency in suppressing the twelve-poles rotor vibrations along the σ axis (i.e., the controller was able to eliminate the rotor oscillations at any angular speed Ω = ω + σ , 0.3 σ 0.3 ). Finally, the best tuning conditions between natural frequencies of both the rotor system ( ω ) and the suggested control technique ( ω 1 and ω 2 ) are shown in Figure 10, where the rotor vibration amplitudes ( a 1 and a 2 ) are plotted in 3D space against the variables σ and σ 1 = σ 2 . By examining Figure 10a,b, we deduced that the smallest oscillation amplitudes of the rotor system (i.e., a 1 = a 2 0 ) occurred along the dashed line that had the equation σ = σ 1 = σ 2 . Therefore, the best working condition of the introduced control algorithm occurred if σ = σ 1 = σ 2 . Accordingly, one can conclude from Equation (53) that the optimum tuning conditions (i.e., σ = σ 1 = σ 2 ) occurred when adjusting the controller’s natural frequencies ( ω 1 and ω 2 ) had the same value of rotor angular speed ( Ω ). Accordingly, the combined control algorithm eliminated the rotor vibrations close to zero, regardless of the excitation force amplitude and its angular speed, if the control gains and the tuning condition were applied, as discussed above.

5. Numerical Simulations and Comparative Study

Numerical validations for all of the obtained results in Section 4 were validated numerically via solving the temporal equations of the closed-loop system (i.e., Equations (19) to (24)), using the Rung–Kutta method. In addition, the performances of the different control algorithms in eliminating the twelve-poles system vibrations were compared. It is worth mentioning that the small circles illustrated in Figure 11 represents the steady-state numerical solution of Equations (19) to (24). This numerical solution was obtained via solving Equations (19) to (24) numerically, using the ODE MATLAB solver for a long time-period until reaching the steady-state response at the different values of σ (notice that Ω = ω + σ ). Then, the maximum temporal vibration amplitudes at steady-state were captured as the steady-state vibration amplitudes (i.e., a 1 = max   ( z 1 ( t ) ) , a 2 = max ( z 2 ( t ) ) , a 3 = max ( z 3 ( t ) ) , a 4 = max ( z 4 ( t ) ) ).
Figure 11 compares the rotor dynamics in the case of both the P D - and the P D + P P F -control algorithms, when the excitation force f = 0.013 . The excellent correspondence between the numerical solutions (i.e., small circles) obtained by solving Equations (19)–(22) and the analytical solutions (i.e., solid and dotted lines) obtained by solving the algebraic system provided by Equation (74) is clear. In addition, the figure demonstrates that the coupling of the P P F -control algorithm with the PD controller eliminated the strong vibration amplitudes of the rotor at the resonance condition (i.e., when σ 0 ); however, two resonant peaks appeared on both sides of σ = 0 . Accordingly, we concluded that the P D + P P F -control algorithm had high efficiency in eliminating the rotor’s undesired vibrations at the perfect resonance case (i.e., when σ = 0 or, in other words, when the angular speed Ω was equal to the system’s natural frequency ω , Ω = ω + σ ). However, if the resonant condition was lost (i.e., Ω ω ), the controller may pump excessive vibratory energy to the rotor system, rather than suppress it (see, for example, Figure 11a,b at σ = 0.1 ).
The instantaneous oscillations of the controlled twelve-poles system in the case of both the P D - and the P D + P P F -control algorithms are simulated in Figure 12, according to Figure 11, at σ = 0.0 , f = 0.013 , and Ω = ω . The figure was obtained by solving Equations (19) to (22) numerically on the time interval 0 t < 500 and turning off the P P F -control algorithm (i.e., with setting η 1 = η 3 = η 5 = η 6 = 0 ); then, at the instant t = 500 , the P P F -control algorithm was turned on via setting η 1 = η 3 = η 5 = η 6 = 0.2 along the period 500 t 1000 . Figure 12a,b illustrates the instantaneous oscillations of the twelve-poles system in the case of both the P D -control algorithm on the time interval 0 t < 500 and the P D + P P F -control algorithm on the time interval 500 t 1000 , while Figure 12c shows the rotor whirling orbit before and after turning on the PPF-control algorithm. Figure 12d compares the vibration amplitude of the rotor system in the case of both the P D and P D + P P F -control algorithms. In addition, Figure 12e,f illustrates the temporal oscillations of the P P F controller. It is clear from the figure that the high instantaneous oscillations of the rotor system (i.e., z 1 ( t ) and z 2 ( t ) ) in the case of the P D controller only were suppressed close to zero when the P P F -control algorithm was activated at the time instant t = 500 , where the rotor vibration energy was channeled to P P F controller.
Figure 13 illustrates the instantaneous oscillatory behaviors of the twelve-poles system in the case of both the PD-control algorithm only and the P D + P P F -control algorithm, according to Figure 11, when σ = 0.1 (i.e., when the perfect resonance condition is lost, Ω = ω + 0.1 ). Therefore, Figure 13 is a repetition of Figure 12, but σ = 0.1 . It is clear from Figure 13a,b that the twelve-poles system exhibited small vibration amplitudes on the time interval 0 t < 500 , as long as the P D controller only was activated. However, the figures demonstrate that the activation of the P P F controller along with P D controller on the time interval 500 t 1000 increased the rotor lateral vibration rather than suppressing it, which agrees with Figure 11 at σ = 0.1 .
The steady-state oscillatory motion of the rotor system in the case of both the P D -control algorithm only and the P D + I R C -control algorithm is compared in Figure 14, when f = 0.013 . It is clear from the figure that the high oscillation amplitudes that occurred at the resonance case (i.e., when σ 0 ) in the case of the PD-control algorithm only was mitigated to small lateral oscillations when the combined P D + I R C -control algorithm was activated. However, while the P D + I R C -control algorithm can mitigate the rotor vibrations along the σ axis, it cannot eliminate the rotor vibration close to zero at the resonant condition, as in the case of the P D + P P F -control algorithm.
Numerical simulations for the instantaneous lateral vibrations of the rotor system (i.e., z 1 ( t ) and z 2 ( t ) ) and the IRC-control algorithm (i.e., z 5 ( t ) and z 6 ( t ) ) are illustrated in Figure 15 and Figure 16, according to Figure 14 when σ = 0.0 and σ = 0.1 , respectively. The two figures were obtained via solving Equations (19), (20), (23), and (24) using ODE45 MATLAB solver on the time interval 0 t < 700 and deactivating the I R C -control algorithm (i.e., when η 2 = η 4 = η 7 = η 8 = 0 ), while at t = 700 the I R C controller was turned on by setting η 2 = η 4 = 0.2 , η 7 = η 8 = 1 on the interval 700 t 1000 . One can note from Figure 15 that the strong instantaneous vibrations of the system (i.e., z 1 ( t ) and z 2 ( t ) ) in the case of the P D -control technique at σ = 0 was reduced to small values (but not close to zero) when the I R C -control algorithm was turned on at t = 700 and the rotor vibration energy was partially transferred to the I R C controller. On the other hand, Figure 16 shows that the I R C -control algorithm also reduced the rotor vibrations to a small value when σ = 0.1 rather than pumping more excess energy to the system, as in the case of the P D + P P F -control technique (see Figure 13).
Finally, the rotor dynamics in the case of both the P D - and the P D + I R C + P P F -control algorithms are compared in Figure 17, when f = 0.013 . It is clear from the figure that the high oscillation amplitudes of the rotor system in the vicinity of σ = 0 in the case of the P D -control algorithm only have been eliminated close to zero, when the P D + I R C + P P F -control algorithm is considered. In addition, the resonant peaks that appeared in Figure 11 (i.e., in the case of P D + P P F -control) were also suppressed, as shown in Figure 17. In other words, the P D + I R C + P P F -control algorithm had all the advantages of the individual control algorithms P D ,   I R C , and P P F , while avoiding their drawbacks.
Figure 18 and Figure 19 compare the instantaneous oscillations of the rotor system in the care of both the P D - and the P D + I R C + P P F -control algorithms, according to Figure 17, when f = 0.013 at σ = 0 and σ = 0.1 , respectively. Figure 18 was obtained by solving Equations (19) to (24) numerically, using the ODE45 solver along the time interval 0 t < 700 and activating the P D controller only (i.e., P = 1.5 , d = 0.005 , and η k = 0 ,   k = 1 , 2 , , 4 ); then, at the time instant t = 700 , the I R C + P P F -control algorithm was turned on, along with the P D controller, via setting η 1 = η 2 = η 3 = η 4 = 0.2 on the time interval 700 t 1000 . Figure 19 is a repetition of Figure 18, but when σ = 0.1 rather than σ = 0.0 . By examining Figure 18a–c, one can notice that the high oscillation amplitudes of the twelve-poles system were eliminated close to zero at a very small transient time as soon as the P D + I R C + P P F controller was turned on. In addition, Figure 19 demonstrates that the P D + I R C + P P F -control algorithm did not add excessive energy to the rotor system when the resonant condition was lost (i.e., when = 0.1 ).

6. Conclusions

In this article, three different control techniques were introduced to eliminate the undesired vibrations of the twelve-poles electro-magnetic suspension system. The introduced control algorithms were the P D , I R C , and P P F controllers and their different combinations (i.e., P D + I R C ,   P D + P P F ,   P D + I R C + P P F ). Relying on the classical mechanics’ principle, the dynamical model that governs the controlled twelve-poles rotor was established as a nonlinear four-degree-of-freedom system that is coupled to two first-order filters. Then, an approximate analytical solution for the controlled system mathematical model was obtained using the asymptotic analysis. Based on the derived analytical solution, the efficiency of the different control algorithms in suppressing the undesired vibrations and improving the bifurcation characteristics of the considered twelve-poles system was explored. In addition, numerical simulations were performed to confirm the accuracy of the obtained analytical investigations, as well as to explore the transient oscillatory behaviors of the rotor system with the different control strategies. Based on our analysis and the discussions above, we reached the following conclusions:
  • The rotor system responds as a linear dynamical system with small vibration amplitudes in the case of the PD-control algorithm, as long as the excitation force f < 0.004 .
  • When only the P D -control algorithm is activated, the twelve-poles rotor behaves like a hardening duffing oscillator, and the nonlinearities dominate its response when the rotor is exposed to a considerable excitation force amplitude (i.e., f > 0.004 ) at the resonance condition. In addition, the electro-magnetic suspension system may suffer from rub and/or impact force between the rotor and the stator if f > 0.013 in the case of P D -control algorithm.
  • Integrating the P P F -control algorithm with P -controller can eliminate the rotor’s undesired vibrations at the resonance condition (i.e., when Ω ω ,   σ 0 ) to negligible oscillation amplitudes, regardless of the excitation force magnitude, but two undesired resonant peaks appear on both sides of σ = 0.0 that may result in high vibrations for the rotor system if the resonant condition is lost (i.e., if Ω ω ).
  • The I R C + P D -control algorithm can mitigate the undesired vibrations and eliminate the nonlinear bifurcation behaviors of the twelve-poles system. However, the main drawback of this controller is that the rotor may perform high oscillation amplitude at the perfect resonance (i.e., when Ω ω ,   σ 0 ).
  • Utilizing the three control algorithms (i.e., P D + I R C + P P F ) as one control strategy eliminated the high oscillation amplitudes of the rotor system close to zero at the perfect resonance conditions. In addition, the resonant peaks that appeared in the case of P D + P P F controller were also suppressed close to zero.
  • The P D + I R C + P P F -control algorithm has all the advantages of the individual control algorithms, P D ,   P D + I R C and P D + P P F , while avoiding their drawbacks.
  • Although both the P D + P P F and P D + I R C + P P F -control algorithms can eliminate the nonlinear vibrations of the twelve-poles system at the perfect resonance condition, the P D + I R C + P P F has the advantage of having the short transient time in suppressing this undesired motion.
  • Tuning the natural frequencies ( ω 1 and ω 2 ) of the P D + I R C + P P F -control algorithm to be close to or equal to the rotor angular speed ( Ω ) guarantees the elimination of the system’s lateral vibrations, regardless of the excitation force magnitude.

Author Contributions

Conceptualization, N.A.S., J.A. and K.A.G.; data curation, N.A.S.; formal analysis, S.M.E.-S. and J.A.; funding acquisition, J.A. and K.A.G.; investigation, N.A.S., S.M.E.-S., J.A. and K.A.G.; methodology, N.A.S., S.M.E.-S. and J.A.; project administration, M.K., K.R.R., J.A. and K.A.G.; resources, N.A.S., M.K. and J.A.; software, N.A.S., S.M.E.-S. and J.A.; supervision, N.A.S., M.K., K.R.R., J.A. and K.A.G.; validation, N.A.S., S.M.E.-S. and J.A.; visualization, N.A.S., S.M.E.-S. and J.A.; writing—original draft, N.A.S., S.M.E.-S. and J.A.; writing—review and editing, N.A.S., J.A. and K.A.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Researchers Supporting Project, Taif University, Taif, Saudi Arabia, under grant number TURSP-2020/16. This work was also supported by the National Science Centre, Poland, under grant number OPUS 14 No. 2017/27/B/ST8/01330.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors thank to Taif University Researchers for Supporting Project number (TURSP-2020/16), Taif, Saudi Arabia.

Conflicts of Interest

The authors declared no potential conflict of interest with respect to the research, authorship, and/or publication of this article.

Abbreviations

z 1 ,   z ˙ 1 ,   z ¨ 1 Normalized displacement, velocity, and acceleration of the twelve-poles system in the X direction.
z 2 ,   z ˙ 2 ,   z ¨ 2 Normalized displacement, velocity, and acceleration of the twelve-poles system in the Y direction.
z 3 ,   z ˙ 3 ,   z ¨ 3 Normalized displacement, velocity, and acceleration of the P P F -control algorithm that connected to the twelve-poles system in the X direction.
z 4 ,   z ˙ 4 ,   z ¨ 4 Normalized displacement, velocity, and acceleration of the P P F -control algorithm that connected to the twelve-poles system in the Y direction.
z 5 ,   z ˙ 5 Normalized displacement, and velocity of the I R C -control algorithm that connected to the twelve-poles system in the X direction.
z 6 ,   z ˙ 6 Normalized displacement, and velocity of the I R C -control algorithm that connected to the twelve-poles system in the Y direction.
μ Normalized damping parameter of the twelve-poles rotor system.
μ 1 ,   μ 2 Normalized damping parameters of the P P F -control algorithms.
ω The normalized natural frequency of the twelve-poles rotor system.
ω 1 ,   ω 2 Normalized natural frequencies of the P P F -control algorithms.
ω 3 ,   ω 4 Normalized Internal-loop feedback gains of the I R C -control algorithms.
Ω The normalized angular speed of the twelve-poles rotor system.
f Normalized excitation force of the twelve-poles rotor system.
P ,   d Normalized proportional and derivative control gains of the P D -control algorithm, respectively.
η 1 ,   η 3 Normalized control gains of the P P F -control algorithms.
η 2 ,   η 4 Normalized control gains of the I R C -control algorithms.
η 5 , η 6 Normalized feedback gains of the P P F -control algorithms.
η 7 , η 8 Normalized feedback gains of the I R C -control algorithms.
α j ,   j = 1 , , 7 Normalized nonlinear coupling coefficients due to the P D -control algorithm.
β j ,   j = 1 , , 16 Normalized nonlinear coupling coefficients due to both the I R C and P P F control algorithms in the X direction.
γ j ,   j = 1 , , 16 Normalized nonlinear coupling coefficients due to both the I R C and P P F control algorithms in the Y direction.
a 1 ,   a 2 Normalized vibration amplitudes of the twelve-poles rotor system in the X and Y directions, respectively.
ϕ 1 ,   ϕ 2 Phase angles of the twelve-poles rotor system in the X and Y directions, respectively.
a 3 ,   a 4 Normalized vibration amplitudes of the P P F -control algorithms in the X and Y directions, respectively.
ϕ 3 ,   ϕ 4 Phase angles of the P P F -control algorithms in the X and Y directions, respectively.
σ Difference between the angular speed ( Ω ) and the normalized natural frequency ( ω ) : σ = Ω ω .

Appendix A

Expanding Equations (11) to (16), using the Maclaurin series up to the third-order approximation, yields the following:
f 1 4 c 0 5 μ 0 N 2 A cos ( φ ) [ ( I 0 2 c 0 2 sin ( α ) ) y + ( k 1 I 0 c 0 3 I 0 2 c 0 2 cos ( α ) ) x + ( 2 I 0 2 cos 3 ( α ) k 1 2 c 0 2 cos ( α ) + 3 k 1 I 0 c 0 cos 2 ( α ) ) x 3 + ( 2 I 0 2 sin ( α ) 2 I 0 2 sin ( α ) cos 2 ( α ) ) y 3 + ( 3 k 1 I 0 c 0 3 k 1 I 0 c 0 cos 2 ( α ) 6 I 0 2 cos ( α ) + 6 I 0 2 cos 3 ( α ) ) x y 2 + ( 6 I 0 2 cos 2 ( α ) sin ( α ) 6 k 1 I 0 c 0 cos ( α ) sin ( α ) + k 1 2 c 0 2 sin ( α ) ) x 2 y + ( 3 k 2 I 0 c 0 3 k 2 I 0 c 0 cos 2 ( α ) ) y 2 x ˙ + ( k 2 2 c 0 2 sin ( α ) ) y x ˙ 2 + ( k 2 2 c 0 2 cos ( α ) ) x x ˙ 2 + ( k 2 I 0 c 0 3 ) x ˙ + ( 3 k 2 I 0 c 0 cos 2 ( α ) 2 k 1 k 3 c 0 2 cos ( α ) ) x 2 x ˙ + ( 6 k 2 I 0 c 0 cos ( α ) sin ( α ) + 2 k 1 k 2 c 0 2 sin ( α ) ) x x ˙ y + ( 2 k 2 k 3 c 0 2 cos ( α ) + 3 k 3 I 0 c 0 cos 2 ( α ) ) x 2 u 1 + ( 2 k 2 k 4 c 0 2 cos ( α ) + 3 k 4 I 0 c 0 cos 2 ( α ) ) x 2 u 2 + ( 3 k 3 I 0 c 0 3 k 3 I 0 c 0 cos 2 ( α ) ) y 2 u 1 + ( 3 k 4 I 0 c 0 3 k 4 I 0 c 0 cos 2 ( α ) ) y 2 u 2 + ( 2 k 2 k 3 c 0 2 sin 2 ( α ) 6 k 3 I 0 c 0 sin ( α ) cos ( α ) ) x y u 1 + + ( 2 k 2 k 4 c 0 2 sin 2 ( α ) 6 k 4 I 0 c 0 sin ( α ) cos ( α ) ) x y u 2 + ( k 3 2 c 0 2 cos ( α ) ) x u 1 2 + ( k 4 2 c 0 2 cos ( α ) ) x u 2 2 + ( k 3 2 c 0 2 sin ( α ) ) y u 1 2 + ( k 4 2 c 0 2 sin ( α ) ) y u 2 2 + ( 2 k 3 k 4 c 0 2 cos ( α ) ) x u 1 u 2 + ( 2 k 3 k 4 c 0 2 sin ( α ) ) y u 1 u 2 + ( 2 k 1 k 3 c 0 2 cos ( α ) 2 k 1 k 4 c 0 2 cos ( α ) ) x x ˙ u 1 + ( 2 k 1 k 3 c 0 2 sin ( α ) + 2 k 1 k 4 c 0 2 sin ( α ) ) y x ˙ u 1 + ( k 3 I 0 c 0 3 ) u 1 + ( k 4 I 0 c 0 3 ) u 2 ]
f 2 4 c 0 5 μ 0 N 2 A cos ( φ ) [ ( k 1 I 0 c 0 3 I 0 2 c 0 2 ) x + ( 2 I 0 2 + k 1 2 c 0 2 3 k 1 I 0 c 0 ) x 3 + ( k 2 2 c 0 2 ) x x ˙ 2 + ( k 2 I 0 c 0 3 ) x ˙ + ( 3 k 2 I 0 c 0 + 2 k 1 k 3 c 0 2 ) x 2 x ˙ + ( 2 k 2 k 3 c 0 2 + 3 k 3 I 0 c 0 ) x 2 u 1 + ( 2 k 2 k 4 c 0 2 + 3 k 4 I 0 c 0 ) x 2 u 2 + ( k 3 2 c 0 2 ) x u 1 2 + ( k 4 2 c 0 2 ) x u 2 2 + ( 2 k 1 k 3 c 0 2 ) x x ˙ u 1 + ( 2 k 1 k 4 c 0 2 ) x x ˙ u 2 + ( 2 k 3 k 4 c 0 2 ) x u 1 u 2 + ( k 3 I 0 c 0 3 ) u 1 + ( k 3 I 0 c 0 3 ) u 2 ]
f 3 = 4 c 0 5 μ 0 N 2 A cos ( φ ) [ ( I 0 2 c 0 2 sin ( α ) ) y + ( k 1 I 0 c 0 3 + I 0 2 c 0 2 cos ( α ) ) x + ( 2 I 0 2 cos 3 ( α ) + k 1 2 c 0 2 cos ( α ) 3 k 1 I 0 c 0 cos 2 ( α ) ) x 3 + ( 2 I 0 2 sin ( α ) 2 I 0 2 sin ( α ) cos 2 ( α ) ) y 3 + ( 3 k 1 I 0 c 0 + 3 k 1 I 0 c 0 cos 2 ( α ) + 6 I 0 2 cos ( α ) 6 I 0 2 cos 3 ( α ) ) x y 2 + ( 6 I 0 2 cos 2 ( α ) sin ( α ) 6 k 1 I 0 c 0 cos ( α ) sin ( α ) + k 1 2 c 0 2 sin ( α ) ) x 2 y + ( 3 k 2 I 0 c 0 + 3 k 2 I 0 c 0 cos 2 ( α ) ) y 2 x ˙ + ( k 2 2 c 0 2 sin ( α ) ) y x ˙ 2 + ( k 2 2 c 0 2 cos ( α ) ) x x ˙ 2 + ( k 2 I 0 c 0 3 ) x ˙ + ( 3 k 2 I 0 c 0 cos 2 ( α ) + 2 k 1 k 3 c 0 2 cos ( α ) ) x 2 x ˙ + ( 6 k 2 I 0 c 0 cos ( α ) sin ( α ) + 2 k 1 k 2 c 0 2 sin ( α ) ) x x ˙ y + ( 2 k 2 k 3 c 0 2 cos ( α ) + 3 k 3 I 0 c 0 cos 2 ( α ) ) x 2 u 1 + ( 2 k 2 k 4 c 0 2 cos ( α ) + 3 k 4 I 0 c 0 cos 2 ( α ) ) x 2 u 2 + ( 3 k 3 I 0 c 0 3 k 3 I 0 c 0 cos 2 ( α ) ) y 2 u 1 + ( 3 k 4 I 0 c 0 3 k 4 I 0 c 0 cos 2 ( α ) ) y 2 u 2 + ( 2 k 2 k 3 c 0 2 sin 2 ( α ) + 6 k 3 I 0 c 0 sin ( α ) cos ( α ) ) x y u 1 + ( 2 k 2 k 4 c 0 2 sin 2 ( α ) + 6 k 4 I 0 c 0 sin ( α ) cos ( α ) ) x y u 2 + ( k 3 2 c 0 2 cos ( α ) ) x u 1 2 + ( k 4 2 c 0 2 cos ( α ) ) x u 2 2 + ( k 3 2 c 0 2 sin ( α ) ) y u 1 2 + ( k 4 2 c 0 2 sin ( α ) ) y u 2 2 + ( 2 k 3 k 4 c 0 2 cos ( α ) ) x u 1 u 2 + ( 2 k 3 k 4 c 0 2 sin ( α ) ) y u 1 u 2 + ( 2 k 1 k 3 c 0 2 cos ( α ) 2 k 1 k 4 c 0 2 cos ( α ) ) x x ˙ u 1 + ( 2 k 1 k 3 c 0 2 sin ( α ) 2 k 1 k 4 c 0 2 sin ( α ) ) y x ˙ u 1 + ( k 3 I 0 c 0 3 ) u 1 + ( k 4 I 0 c 0 3 ) u 2 ]
f 4 = 4 c 0 5 μ 0 N 2 A cos ( φ ) [ ( I 0 2 c 0 2 sin ( α ) ) x + ( k 1 I 0 c 0 3 I 0 2 c 0 2 cos ( α ) ) y + ( 2 I 0 2 cos 3 ( α ) k 1 2 c 0 2 cos ( α ) + 3 k 1 I 0 c 0 cos 2 ( α ) ) y 3 + ( 2 I 0 2 sin ( α ) + 2 I 0 2 sin ( α ) cos 2 ( α ) ) x 3 + ( 3 k 1 I 0 c 0 3 k 1 I 0 c 0 cos 2 ( α ) 6 I 0 2 cos ( α ) + 6 I 0 2 cos 3 ( α ) ) y x 2 + ( 6 I 0 2 cos 2 ( α ) sin ( α ) + 6 k 1 I 0 c 0 cos ( α ) sin ( α ) k 1 2 c 0 2 sin ( α ) ) y 2 x + ( 3 k 2 I 0 c 0 3 k 2 I 0 c 0 cos 2 ( α ) ) x 2 y ˙ + ( k 2 2 c 0 2 sin ( α ) ) x y ˙ 2 + ( k 2 2 c 0 2 cos ( α ) ) y y ˙ 2 + ( k 2 I 0 c 0 3 ) y ˙ + ( 3 k 2 I 0 c 0 cos 2 ( α ) 2 k 1 k 3 c 0 2 cos ( α ) ) y 2 y ˙ + ( 6 k 2 I 0 c 0 cos ( α ) sin ( α ) 2 k 1 k 2 c 0 2 sin ( α ) ) y y ˙ x + ( 2 k 2 k 5 c 0 2 cos ( α ) + 3 k 5 I 0 c 0 cos 2 ( α ) ) y 2 v 1 + ( 2 k 2 k 6 c 0 2 cos ( α ) + 3 k 6 I 0 c 0 cos 2 ( α ) ) y 2 v 2 + ( 3 k 5 I 0 c 0 3 k 5 I 0 c 0 cos 2 ( α ) ) x 2 v 1 + ( 3 k 6 I 0 c 0 3 k 6 I 0 c 0 cos 2 ( α ) ) x 2 v 2 + ( 2 k 2 k 5 c 0 2 sin 2 ( α ) + 6 k 5 I 0 c 0 sin ( α ) cos ( α ) ) x y v 1 + ( 2 k 2 k 6 c 0 2 sin 2 ( α ) + 6 k 6 I 0 c 0 sin ( α ) cos ( α ) ) x y v 2 + ( k 5 2 c 0 2 cos ( α ) ) y v 1 2 + ( k 6 2 c 0 2 cos ( α ) ) y v 2 2 + ( k 5 2 c 0 2 sin ( α ) ) x v 1 2 + ( k 6 2 c 0 2 sin ( α ) ) x v 2 2 + ( 2 k 5 k 6 c 0 2 cos ( α ) ) y v 1 v 2 + ( 2 k 5 k 6 c 0 2 sin ( α ) ) x v 1 v 2 + ( 2 k 1 k 6 c 0 2 cos ( α ) 2 k 1 k 5 c 0 2 cos ( α ) ) y y ˙ v 1 + ( 2 k 1 k 6 c 0 2 sin ( α ) + 2 k 1 k 5 c 0 2 sin ( α ) ) x y ˙ v 1 + ( k 5 I 0 c 0 3 ) v 1 + ( k 6 I 0 c 0 3 ) v 2 ]
f 5 = 4 c 0 5 μ 0 N 2 A cos ( φ ) [ ( k 1 I 0 c 0 3 I 0 2 c 0 2 ) y + ( 2 I 0 2 + k 1 2 c 0 2 3 k 1 I 0 c 0 ) y 3 + ( k 2 2 c 0 2 ) y y ˙ 2 + ( k 2 I 0 c 0 3 ) y ˙ + ( 3 k 2 I 0 c 0 + 2 k 1 k 3 c 0 2 ) y 2 y ˙ + ( 2 k 2 k 5 c 0 2 + 3 k 5 I 0 c 0 ) y 2 v 1 + ( 2 k 2 k 6 c 0 2 + 3 k 6 I 0 c 0 ) y 2 v 2 + ( k 5 2 c 0 2 ) y v 1 2 + ( k 6 2 c 0 2 ) y v 2 2 + ( 2 k 1 k 5 c 0 2 ) y y ˙ v 1 + ( 2 k 1 k 6 c 0 2 ) y y ˙ v 2 + ( 2 k 5 k 6 c 0 2 ) y v 1 v 2 + ( k 5 I 0 c 0 3 ) v 1 + ( k 6 I 0 c 0 3 ) v 2 + ( O ) 3 ]
f 6 = 4 c 0 5 μ 0 N 2 A cos ( φ ) [ ( I 0 2 c 0 2 sin ( α ) ) x + ( k 1 I 0 c 0 3 + I 0 2 c 0 2 cos ( α ) ) y + ( 2 I 0 2 cos 3 ( α ) + k 1 2 c 0 2 cos ( α ) 3 k 1 I 0 c 0 cos 2 ( α ) ) y 3 + ( 2 I 0 2 sin ( α ) + 2 I 0 2 sin ( α ) cos 2 ( α ) ) x 3 + ( 3 k 1 I 0 c 0 + 3 k 1 I 0 c 0 cos 2 ( α ) + 6 I 0 2 cos ( α ) 6 I 0 2 cos 3 ( α ) ) y x 2 + ( 6 I 0 2 cos 2 ( α ) sin ( α ) + 6 k 1 I 0 c 0 cos ( α ) sin ( α ) k 1 2 c 0 2 sin ( α ) ) y 2 x + ( 3 k 2 I 0 c 0 + 3 k 2 I 0 c 0 cos 2 ( α ) ) x 2 y ˙ + ( k 2 2 c 0 2 sin ( α ) ) x y ˙ 2 + ( k 2 2 c 0 2 cos ( α ) ) y y ˙ 2 + ( k 2 I 0 c 0 3 ) y ˙ + ( 3 k 2 I 0 c 0 cos 2 ( α ) + 2 k 1 k 3 c 0 2 cos ( α ) ) y 2 y ˙ + ( 6 k 2 I 0 c 0 cos ( α ) sin ( α ) 2 k 1 k 2 c 0 2 sin ( α ) ) y y ˙ x + ( 2 k 2 k 5 c 0 2 cos ( α ) + 3 k 5 I 0 c 0 cos 2 ( α ) ) y 2 v 1 + ( 2 k 2 k 6 c 0 2 cos ( α ) + 3 k 6 I 0 c 0 cos 2 ( α ) ) y 2 v 2 + ( 3 k 5 I 0 c 0 3 k 5 I 0 c 0 cos 2 ( α ) ) x 2 v 1 + ( 3 k 6 I 0 c 0 3 k 6 I 0 c 0 cos 2 ( α ) ) x 2 v 2 + ( 2 k 2 k 5 c 0 2 sin 2 ( α ) 6 k 5 I 0 c 0 sin ( α ) cos ( α ) ) x y v 1 + ( 2 k 2 k 6 c 0 2 sin 2 ( α ) 6 k 6 I 0 c 0 sin ( α ) cos ( α ) ) x y v 2 + ( k 5 2 c 0 2 cos ( α ) ) y v 1 2 + ( k 6 2 c 0 2 cos ( α ) ) y v 2 2 + ( k 5 2 c 0 2 sin ( α ) ) x v 1 2 + ( k 6 2 c 0 2 sin ( α ) ) x v 2 2 + ( 2 k 5 k 6 c 0 2 cos ( α ) ) y v 1 v 2 + ( 2 k 5 k 6 c 0 2 sin ( α ) ) x v 1 v 2 + ( 2 k 1 k 6 c 0 2 cos ( α ) 2 k 1 k 5 c 0 2 cos ( α ) ) y y ˙ v 1 + ( 2 k 1 k 6 c 0 2 sin ( α ) + 2 k 1 k 5 c 0 2 sin ( α ) ) x y ˙ v 1 + ( k 5 I 0 c 0 3 ) v 1 + ( k 6 I 0 c 0 3 ) v 2 ]

Appendix B

μ = 1 2 ( 2 d cos ( α ) + d ) , ω = 2 p cos ( α ) + p 3 ,
α 1 = 6 p cos 3 ( α ) 3 p + 2 p 2 cos 2 ( α ) + 6 + p 2 + 8 cos 4 ( α ) 8 cos 2 ( α ) ,
α 2 = 2 p 2 2 p 2 cos 2 ( α ) + 24 cos 2 ( α ) 24 cos 4 ( α ) 18 p cos ( α ) + 18 p cos 3 ( α ) ,
α 3 = 3 d + 2 p d + 4 p d cos 2 ( α ) 6 d cos 3 ( α ) , α 4 = 6 d cos 3 ( α ) 6 d cos ( α ) ,
α 5 = 2 d 2 sin 2 ( α ) , α 6 = d 2 ( 1 + 2 cos 2 ( α ) ) ,
α 7 = 4 d ( 3 cos 3 ( α ) 3 cos ( α ) + p sin 2 ( α ) ) , β 1 = η 1 1 + 2 cos ( α ) ( 3 2 p + 6 cos 3 ( α ) 4 p cos 2 ( α ) ) ,
γ 1 = η 3 1 + 2 cos ( α ) ( 3 2 p + 6 cos 3 ( α ) 4 p cos 2 ( α ) ) , β 2 = η 1 1 + 2 cos ( α ) ( 2 d 4 d cos 2 ( α ) ) ,
γ 2 = η 3 1 + 2 cos ( α ) ( 2 d 4 d cos 2 ( α ) ) , β 3 = η 1 2 ( 1 + 2 cos ( α ) ) 2 ( 1 + 2 cos 2 ( α ) ) ,
γ 3 = η 3 2 ( 1 + 2 cos ( α ) ) 2 ( 1 + 2 cos 2 ( α ) ) , β 4 = 4 η 3 1 + 2 cos ( α ) ( p sin 2 ( α ) + 3 cos ( α ) sin 2 ( α ) ) ,
γ 4 = 4 η 1 1 + 2 cos ( α ) ( sin 2 ( α ) + 3 cos ( α ) sin 2 ( α ) ) , β 5 = 2 η 3 2 sin 2 ( α ) ( 1 + 2 cos ( α ) ) 2 ,
γ 5 = 2 η 1 2 sin 2 ( α ) ( 1 + 2 cos ( α ) ) 2 , β 6 = 4 d η 3 sin 2 ( α ) 1 + 2 cos ( α ) ,
γ 6 = 4 d η 1 sin 2 ( α ) 1 + 2 cos ( α ) , β 7 = 6 η 1 cos ( α ) sin 2 ( α ) 1 + 2 cos ( α ) ,
γ 7 = 6 η 3 cos ( α ) sin 2 ( α ) 1 + 2 cos ( α ) , β 8 = η 1 η 2 ( 1 + 2 cos ( α ) ) 2 ( 2 4 cos 2 ( α ) ) ,
γ 8 = η 3 η 4 ( 1 + 2 cos ( α ) ) 2 ( 2 4 cos 2 ( α ) ) , β 9 = η 2 1 + 2 cos ( α ) ( 2 d + 4 d cos 2 ( α ) ) ,
γ 9 = η 3 1 + 2 cos ( α ) ( 2 d + 4 d cos 2 ( α ) ) , β 10 = 4 η 3 η 4 sin 2 ( α ) ( 1 + 2 cos ( α ) ) 2 ,
γ 10 = 4 η 1 η 2 sin 2 ( α ) ( 1 + 2 cos ( α ) ) 2 , β 11 = η 2 1 + 2 cos ( α ) ( 3 + 2 p + 4 p cos 2 ( α ) 6 cos 3 ( α ) ) ,
γ 11 = η 4 1 + 2 cos ( α ) ( 3 + 2 p + 4 p cos 2 ( α ) 6 cos 3 ( α ) ) , β 12 = 2 η 4 2 sin 2 ( α ) ( 1 + 2 cos ( α ) ) 2 ,
γ 12 = 2 η 2 2 sin 2 ( α ) ( 1 + 2 cos ( α ) ) 2 , β 13 = η 4 1 + 2 cos ( α ) ( 4 p sin 2 ( α ) 12 cos ( α ) sin 2 ( α ) ) ,
γ 13 = η 2 1 + 2 cos ( α ) ( 4 p sin 2 ( α ) 12 cos ( α ) sin 2 ( α ) ) , β 14 = 4 d η 4 sin 2 ( α ) 1 + 2 cos ( α ) ,
γ 14 = 4 d η 2 sin 2 ( α ) 1 + 2 cos ( α ) , β 15 = η 2 2 ( 1 + 2 cos ( α ) ) 2 ( 1 + 2 cos 2 ( α ) ) ,
γ 15 = η 4 2 ( 1 + 2 cos ( α ) ) 2 ( 1 + 2 cos 2 ( α ) ) , β 16 = 6 η 2 cos ( α ) s i n 2 ( α ) 1 + 2 cos ( α ) ,
γ 16 = 6 η 4 cos ( α ) s i n 2 ( α ) 1 + 2 cos ( α ) .

Appendix C

J 11 = F 1 a 11 = 1 2 ( 2 μ + η 2 η 7 ω 3 2 + ω 2 ) + 3 8 ( α 3 2 β 11 η 7 ω 3 2 + ω 2 2 ω 3 β 15 η 7 2 ( ω 3 2 + ω 2 ) 2 ) a 10 2 + 1 8 ( 2 α 4 β 13 η 8 ω 4 2 + ω 2 ω 4 β 14 η 8 ω 4 2 + ω 2 2 β 16 η 7 ω 3 2 + ω 2 ) a 20 2 + 1 8 ( α 4 + α 7 + 2 ω 4 β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 β 13 η 8 ω 4 2 + ω 2 + ω 4 β 14 η 8 ω 4 2 + ω 2 ) a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) + 1 8 ( α 2 ω α 5 ω + ( ω 4 2 ω 2 ) β 12 η 8 2 ω ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω ( ω 4 2 + ω 2 ) ω β 14 η 8 ω 4 2 + ω 2 ) a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω β 1 a 10 b 10 sin ( ϕ 30 ) 1 8 ω β 3 b 10 2 sin ( 2 ϕ 30 ) + 1 4 ( β 2 2 β 8 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( ϕ 30 ) 1 8 ( β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( ϕ 40 ) + 1 8 ω ( ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( ϕ 40 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 12 = F 1 a 21 = 1 4 ( 2 α 4 β 13 η 8 ω 4 2 + ω 2 ω 4 β 14 η 8 ω 4 2 + ω 2 2 β 16 η 7 ω 3 2 + ω 2 ) a 10 a 20 + 1 4 ( α 4 + α 7 + 2 ω 4 β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 β 13 η 8 ω 4 2 + ω 2 + ω 4 β 14 η 8 ω 4 2 + ω 2 ) a 10 a 20 cos ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ( α 2 ω α 5 ω + ( ω 4 2 ω 2 ) β 12 η 8 2 ω ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω ( ω 4 2 + ω 2 ) ω β 14 η 8 ω 4 2 + ω 2 ) a 10 a 20 sin ( 2 ϕ 10 2 ϕ 20 ) 1 2 ω β 7 a 20 b 10 sin ( ϕ 30 ) 1 8 ( β 10 η 8 ω 4 2 + ω 2 ) a 10 b 20 cos ( ϕ 40 ) + 1 8 ω ( ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 b 20 sin ( ϕ 40 ) + 1 4 ω β 7 a 20 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 10 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) ,
J 13 = F 1 b 11 = ( 1 2 ω η 1 1 8 ω β 1 a 10 2 1 4 ω β 7 a 20 2 ) sin ( ϕ 30 ) 1 4 ω β 3 a 10 b 10 sin ( 2 ϕ 30 ) + 1 8 ( β 2 2 β 8 η 7 ω 3 2 + ω 2 ) a 10 2 cos ( ϕ 30 ) + 1 8 ω β 7 a 20 2 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) ,
J 14 = F 1 b 21 = 1 8 ( β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 cos ( ϕ 40 ) + 1 8 ω ( ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 sin ( ϕ 40 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 a 10 b 20 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 15 = F 1 ϕ 11 = 1 4 ( α 4 + α 7 + 2 ω 4 β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 β 13 η 8 ω 4 2 + ω 2 + ω 4 β 14 η 8 ω 4 2 + ω 2 ) a 10 a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ( α 2 ω α 5 ω + ( ω 4 2 ω 2 ) β 12 η 8 2 ω ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω ( ω 4 2 + ω 2 ) ω β 14 η 8 ω 4 2 + ω 2 ) a 10 a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 a 10 b 20 2 cos ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) + 1 2 ω ( ω + σ ) 2 f cos ( ϕ 10 ) ,
J 16 = F 1 ϕ 21 = 1 4 ( α 4 + α 7 + 2 ω 4 β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 β 13 η 8 ω 4 2 + ω 2 + ω 4 β 14 η 8 ω 4 2 + ω 2 ) a 10 a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) 1 4 ( α 2 ω α 5 ω + ( ω 4 2 ω 2 ) β 12 η 8 2 ω ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω ( ω 4 2 + ω 2 ) ω β 14 η 8 ω 4 2 + ω 2 ) a 10 a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ω β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω β 5 a 10 b 20 2 cos ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 17 = F 1 ϕ 31 = ( 1 2 ω η 1 b 10 1 8 ω β 1 a 10 2 b 10 1 4 ω β 7 a 20 2 b 10 ) cos ( ϕ 30 ) 1 4 ω β 3 a 10 b 10 2 cos ( 2 ϕ 30 ) 1 8 ( β 2 2 β 8 η 7 ω 3 2 + ω 2 ) a 10 2 b 10 sin ( ϕ 30 ) + 1 8 ω β 7 a 20 2 b 10 cos ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) ,
J 18 = F 1 ϕ 41 = 1 8 ( β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 sin ( ϕ 40 ) + 1 8 ω ( ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 cos ( ϕ 40 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 10 a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω β 5 a 10 b 20 2 cos ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 21 = F 2 a 11 = 1 4 ( 2 α 4 γ 13 η 7 ω 3 2 + ω 2 ω 3 γ 14 η 7 ω 3 2 + ω 2 2 γ 16 η 8 ω 4 2 + ω 2 ) a 20 a 10 + 1 4 ( α 4 + α 7 2 ω 3 γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 γ 13 η 7 ω 3 2 + ω 2 + ω 3 γ 14 η 7 ω 3 2 + ω 2 ) a 20 a 10 cos ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ( α 2 ω α 5 ω + ( ω 3 2 ω 2 ) γ 12 η 7 2 ω ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω ( ω 3 2 + ω 2 ) + ω γ 14 η 7 ω 3 2 + ω 2 ) a 20 a 10 sin ( 2 ϕ 20 2 ϕ 10 ) 1 2 ω γ 7 a 10 b 20 sin ( ϕ 40 ) 1 8 ω ( ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 b 10 cos ( ϕ 30 ) + 1 8 ω ( ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 b 10 sin ( ϕ 30 ) + 1 4 ω γ 7 a 10 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 8 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) ,
J 22 = F 2 a 21 = 1 2 ( 2 μ + η 4 η 8 ω 4 2 + ω 2 ) + 3 8 ( α 3 2 γ 11 η 8 ω 4 2 + ω 2 2 ω 4 γ 15 η 8 2 ( ω 4 2 + ω 2 ) 2 ) a 20 2 + 1 8 ( 2 α 4 γ 13 η 7 ω 3 2 + ω 2 ω 3 γ 14 η 7 ω 3 2 + ω 2 2 γ 16 η 8 ω 4 2 + ω 2 ) a 10 2 + 1 8 ( α 4 + α 7 2 ω 3 γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 γ 13 η 7 ω 3 2 + ω 2 + ω 3 γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) + 1 8 ( α 2 ω α 5 ω + ( ω 3 2 ω 2 ) γ 12 η 7 2 ω ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω ( ω 3 2 + ω 2 ) + ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω γ 1 a 20 b 20 sin ( ϕ 40 ) 1 8 ω γ 3 b 20 2 sin ( 2 ϕ 40 ) + 1 4 ω ( γ 2 ω 2 ω γ 8 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( ϕ 40 ) 1 8 ω ( ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 b 10 cos ( ϕ 30 ) + 1 8 ω ( ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( ϕ 30 ) + 1 8 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 ( cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 23 = F 2 b 11 = 1 8 ω ( ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 cos ( ϕ 30 ) + 1 8 ω ( ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 sin ( ϕ 30 ) + 1 8 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 sin ( 2 ϕ 2 2 ϕ 1 ϕ 3 ) + 1 4 ω γ 5 a 20 b 10 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 24 = F 2 b 21 = ( 1 2 ω η 3 1 8 ω γ 1 a 20 2 1 4 ω γ 7 a 10 2 ) sin ( ϕ 40 ) 1 4 ω γ 3 a 20 b 20 sin ( 2 ϕ 40 ) + 1 8 ω ( γ 2 ω 2 ω γ 8 η 8 ω 4 2 + ω 2 ) a 20 2 cos ( ϕ 40 ) + 1 8 ω γ 7 a 10 2 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) ,
J 25 = F 2 ϕ 11 = 1 4 ( α 4 + α 7 2 ω 3 γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 γ 13 η 7 ω 3 2 + ω 2 + ω 3 γ 14 η 7 ω 3 2 + ω 2 ) a 20 a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) 1 4 ( α 2 ω α 5 ω + ( ω 3 2 ω 2 ) γ 12 η 7 2 ω ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω ( ω 3 2 + ω 2 ) + ω γ 14 η 7 ω 3 2 + ω 2 ) a 20 a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω γ 7 a 10 2 b 20 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 4 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ω γ 5 a 20 b 10 2 cos ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 26 = F 2 ϕ 21 = 1 4 ( α 4 + α 7 2 ω 3 γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 γ 13 η 7 ω 3 2 + ω 2 + ω 3 γ 14 η 7 ω 3 2 + ω 2 ) a 20 a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ( α 2 ω α 5 ω + ( ω 3 2 ω 2 ) γ 12 η 7 2 ω ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω ( ω 3 2 + ω 2 ) + ω γ 14 η 7 ω 3 2 + ω 2 ) a 20 a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω γ 7 a 10 2 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 4 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ω γ 5 a 20 b 10 2 cos ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) + 1 2 ω ( ω + σ ) 2 f sin ( ϕ 20 ) ,
J 27 = F 2 ϕ 31 = 1 8 ω ( ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 sin ( ϕ 3 ) + 1 8 ω ( ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 cos ( ϕ 30 ) + 1 8 ω ( γ 6 ω ω γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 20 a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 8 ω γ 5 a 20 b 10 2 cos ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 28 = F 2 ϕ 41 = ( 1 2 ω η 3 b 20 1 8 ω γ 1 a 20 2 b 20 1 4 ω γ 7 a 10 2 b 20 ) cos ( ϕ 40 ) 1 4 ω γ 3 a 20 b 20 2 cos ( 2 ϕ 40 ) 1 8 ω ( γ 2 ω 2 ω γ 8 η 8 ω 4 2 + ω 2 ) a 20 2 b 20 sin ( ϕ 40 ) + 1 8 ω γ 7 a 10 2 b 20 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) ,
J 31 = F 3 a 11 = 1 2 ( ω + σ 1 ) η 5 sin ( ϕ 30 ) ,
J 32 = F 3 a 21 = 0 ,
J 33 = F 3 b 11 = μ 1 ,
J 34 = F 3 b 21 = 0 ,
J 35 = F 3 ϕ 11 = 0 ,
J 37 = F 3 ϕ 31 = 1 2 ( ω + σ 1 ) η 5 a 10 cos ( ϕ 30 ) ,
J 38 = F 3 ϕ 41 = 0 ,
J 41 = F 4 a 11 = 0 ,
J 42 = F 4 a 21 = 1 2 ( ω + σ 2 ) η 6 sin ( ϕ 40 ) ,
J 43 = F 4 b 11 = 0 ,
J 44 = F 4 b 21 = μ 2 ,
J 45 = F 4 ϕ 11 = 0 ,
J 46 = F 4 ϕ 21 = 0 ,
J 47 = F 4 ϕ 31 = 0 ,
J 48 = F 4 ϕ 41 = 1 2 ( ω + σ 2 ) η 6 a 20 cos ( ϕ 40 ) ,
J 51 = F 5 a 11 = 1 4 ω ( 3 α 1 + α 6 ω 2 + 2 ω 3 β 11 η 7 ω 3 2 + ω 2 + ( ω 3 2 ω 2 ) β 15 η 7 2 ( ω 3 2 + ω 2 ) 2 ) a 10 + ( 1 2 a 10 2 η 1 b 10 + 3 8 β 1 b 10 1 4 a 10 2 β 7 a 2 2 b 10 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 b 10 ) cos ( ϕ 30 ) ω 1 8 β 8 η 7 ω 3 2 + ω 2 b 10 sin ( ϕ 30 ) 1 8 β 2 b 10 sin ( ϕ 30 ) 1 8 ω a 10 2 β 7 a 20 2 b 10 cos ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) 1 2 ω a 10 2 ( ω + σ ) 2 f cos ( ϕ 10 ) ,
J 52 = F 5 a 21 = 1 4 ω ( 2 α 2 + 2 α 5 ω 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 ω 2 β 14 η 8 ω 4 2 + ω 2 + 2 ω 3 β 16 η 7 ω 3 2 + ω 2 ) a 20 + 1 4 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 20 cos ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 20 sin ( 2 ϕ 10 2 ϕ 20 ) + 1 4 β 7 a 2 2 b 1 ω a 1 cos ( ϕ 30 ) + 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) b 20 cos ( ϕ 40 ) + 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) b 20 sin ( ϕ 40 ) + 1 4 ω a 10 β 7 a 20 b 10 cos ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) b 2 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) b 2 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) ,
J 53 = F 5 b 11 = 1 2 ω β 3 b 10 + ( 1 2 η 1 + 3 8 β 1 a 10 2 + 1 4 β 7 a 20 2 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 a 10 2 ) cos ( ϕ 30 ) ω a 10 + 1 4 ω β 3 b 10 cos ( 2 ϕ 30 ) 1 8 β 8 η 7 ω 3 2 + ω 2 a 10 sin ( ϕ 30 ) 1 8 β 2 a 10 sin ( ϕ 30 ) + 1 8 ω a 10 β 7 a 20 2 cos ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) ,
J 54 = F 5 b 21 = 1 2 ω β 5 b 20 + 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 cos ( ϕ 40 ) + 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) a 20 sin ( ϕ 40 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 b 20 cos ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 55 = F 5 ϕ 11 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω a 10 β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) 1 4 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) 1 2 ω a 10 ( ω + σ ) 2 f sin ( ϕ 10 ) ,
J 56 = F 5 ϕ 21 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ω a 10 β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) + 1 4 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 57 = F 5 ϕ 31 = 1 ω a 10 ( 1 2 η 1 b 10 + 3 8 β 1 a 10 2 b 10 + 1 4 β 7 a 20 2 b 10 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 a 10 2 b 10 ) sin ( ϕ 30 ) 1 4 ω β 3 b 10 2 sin ( 2 ϕ 30 ) 1 8 β 8 η 7 ω 3 2 + ω 2 a 10 b 10 cos ( ϕ 30 ) 1 8 β 2 a 10 b 10 cos ( ϕ 30 ) 1 8 ω a 10 β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) ,
J 58 = F 5 ϕ 41 = 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( ϕ 40 ) + 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( ϕ 40 ) + 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 61 = F 6 a 11 = 1 4 ω ( 2 α 2 + 2 α 5 ω 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 ω 2 γ 14 η 7 ω 3 2 + ω 2 + 2 ω 4 γ 16 η 8 ω 4 2 + ω 2 ) a 10 + 1 4 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 10 cos ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 sin ( 2 ϕ 20 2 ϕ 10 ) + 1 2 γ 7 a 10 b 20 cos ( ϕ 40 ) ω a 20 + 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) b 10 cos ( ϕ 30 ) + 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) b 10 sin ( ϕ 30 ) + 1 4 ω a 20 γ 7 a 10 b 20 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) ,
J 62 = F 6 a 21 = 1 4 ω ( 3 α 1 + α 6 ω 2 + 2 ω 4 γ 11 η 8 ω 4 2 + ω 2 + ( ω 4 2 ω 2 ) γ 15 η 8 2 ( ω 4 2 + ω 2 ) 2 ) a 20 + ( 1 2 a 20 2 η 3 b 20 + 3 8 γ 1 b 20 1 4 a 20 2 γ 7 a 10 2 b 20 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 b 20 ) cos ( ϕ 40 ) ω 1 8 γ 8 η 8 ω 4 2 + ω 2 b 20 sin ( ϕ 40 ) 1 8 γ 2 b 20 sin ( ϕ 40 ) 1 8 ω a 20 2 γ 7 a 10 2 b 20 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) 1 2 ω a 20 2 ( ω + σ ) 2 f sin ( ϕ 20 ) ,
J 63 = F 6 b 11 = 1 2 ω γ 5 b 10 + 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 cos ( ϕ 30 ) + 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 sin ( ϕ 30 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ω γ 5 b 10 cos ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 64 = F 6 b 21 = 1 2 ω γ 3 b 20 + ( 1 2 η 3 + 3 8 γ 1 a 20 2 + 1 4 γ 7 a 10 2 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 a 20 2 ) cos ( ϕ 40 ) ω a 20 + 1 4 ω γ 3 b 20 cos ( 2 ϕ 40 ) 1 8 γ 8 η 8 ω 4 2 + ω 2 a 20 sin ( ϕ 40 ) 1 8 γ 2 a 20 sin ( ϕ 40 ) + 1 8 ω a 20 γ 7 a 10 2 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) ,
J 65 = F 6 ϕ 11 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ω a 2 γ 7 a 10 2 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 4 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 66 = F 6 ϕ 21 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω a 2 γ 7 a 10 2 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) 1 4 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) + 1 2 ω a 20 ( ω + σ ) 2 f cos ( ϕ 20 ) ,
J 67 = F 6 ϕ 31 = 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( ϕ 30 ) + 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( ϕ 30 ) + 1 4 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) + 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) ,
J 68 = F 6 ϕ 41 = ( 1 2 η 3 b 20 + 3 8 γ 1 a 20 2 b 20 + 1 4 γ 7 a 10 2 b 20 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 a 20 2 b 20 ) sin ( ϕ 40 ) ω a 20 1 4 ω γ 30 b 20 2 sin ( 2 ϕ 40 ) 1 8 γ 8 η 8 ω 4 2 + ω 2 a 20 b 20 cos ( ϕ 40 ) 1 8 γ 2 a 20 b 20 cos ( ϕ 40 ) 1 8 ω a 2 γ 7 a 10 2 b 20 s i n ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) ,
J 71 = F 7 a 11 = 1 2 ( ω + σ 1 ) b 10 η 5 cos ( ϕ 30 ) 1 4 ω ( 3 α 1 + α 6 ω 2 + 2 ω 3 β 11 η 7 ω 3 2 + ω 2 + ( ω 3 2 ω 2 ) β 15 η 7 2 ( ω 3 2 + ω 2 ) 2 ) a 10 ( 1 2 a 10 2 η 1 b 10 + 3 8 β 1 b 10 + 1 4 β 7 b 10 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 b 10 ) cos ( ϕ 30 ) ω + 1 8 β 8 η 7 ω 3 2 + ω 2 b 10 sin ( ϕ 30 ) 1 8 β 2 b 10 sin ( ϕ 30 ) 1 2 ω a 10 2 ( ω + σ ) 2 f cos ( ϕ 10 ) ,
J 72 = F 7 a 21 = 1 4 ω ( 2 α 2 + 2 α 5 ω 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 ω 2 β 14 η 8 ω 4 2 + ω 2 + 2 ω 3 β 16 η 7 ω 3 2 + ω 2 ) a 20 1 4 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 20 cos ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 20 sin ( 2 ϕ 10 2 ϕ 20 ) 1 2 β 7 a 20 b 10 cos ( ϕ 30 ) ω a 10 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) b 20 cos ( ϕ 40 ) 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) b 20 sin ( ϕ 40 ) 1 4 ω a 10 β 7 a 20 b 10 cos ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) ,
J 73 = F 7 b 11 = 1 2 ( ω + σ 1 ) b 10 2 η 5 a 10 cos ( ϕ 30 ) 1 2 ω β 3 b 10 ( 1 2 η 1 + 3 8 β 1 a 10 2 + 1 4 β 7 a 20 2 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 a 10 2 ) cos ( ϕ 30 ) ω a 10 1 4 ω β 3 b 10 cos ( 2 ϕ 30 ) + 1 8 β 8 η 7 ω 3 2 + ω 2 a 10 sin ( ϕ 3 ) 1 8 β 2 a 10 sin ( ϕ 30 ) 1 8 ω a 10 β 7 a 20 2 cos ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) ,
J 74 = F 7 b 21 = 1 2 ω β 5 b 20 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 cos ( ϕ 40 ) 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) a 20 sin ( ϕ 40 ) 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω β 5 b 20 cos ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 75 = F 7 ϕ 11 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ω a 1 β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) + 1 4 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 4 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) + 1 2 ω a 10 ( ω + σ ) 2 f sin ( ϕ 10 ) ,
J 76 = F 7 ϕ 21 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 4 2 ω 2 ) β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω 4 β 13 η 8 ω 4 2 + ω 2 + ω 2 β 14 η 8 ω 4 2 + ω 2 ) a 20 2 sin ( 2 ϕ 10 2 ϕ 20 ) + 1 4 ω ( α 4 ω α 7 ω + 2 ω 4 ω β 12 η 8 2 ( ω 4 2 + ω 2 ) 2 + ω β 13 η 8 ω 4 2 + ω 2 ω 4 ω β 14 η 8 ω 4 2 + ω 2 ) a 20 2 cos ( 2 ϕ 10 2 ϕ 20 ) 1 4 ω a 1 β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) 1 4 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) 1 8 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 77 = F 7 ϕ 31 = ( 1 2 η 1 b 10 + 3 8 β 1 a 10 2 b 10 + 1 4 β 7 a 20 2 b 10 + 1 4 ω 3 β 8 η 7 ω 3 2 + ω 2 a 10 2 b 10 ) sin ( ϕ 30 ) ω a 10 + 1 4 ω β 3 b 10 2 sin ( 2 ϕ 30 ) + 1 8 β 8 η 7 ω 3 2 + ω 2 a 10 b 10 cos ( ϕ 30 ) 1 8 β 2 a 10 b 10 cos ( ϕ 30 ) + 1 8 ω a 10 β 7 a 20 2 b 10 sin ( 2 ϕ 10 2 ϕ 20 + ϕ 30 ) ,
J 78 = F 7 ϕ 41 = 1 8 ω ( 2 β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( ϕ 40 ) 1 8 ( 2 β 6 + β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( ϕ 40 ) 1 8 ω ( β 4 + ω 4 β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 sin ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 8 ω ( β 6 ω + ω β 10 η 8 ω 4 2 + ω 2 ) a 20 b 20 cos ( 2 ϕ 10 2 ϕ 20 ϕ 40 ) + 1 4 ω β 5 b 20 2 sin ( 2 ϕ 10 2 ϕ 20 2 ϕ 40 ) ,
J 81 = F 8 a 11 = 1 4 ω ( 2 α 2 + 2 α 5 ω 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 ω 2 γ 14 η 7 ω 3 2 + ω 2 + 2 ω 4 γ 16 η 8 ω 4 2 + ω 2 ) a 10 1 4 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 10 cos ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 sin ( 2 ϕ 20 2 ϕ 10 ) 1 2 γ 7 a 10 b 20 cos ( ϕ 40 ) ω a 20 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) b 10 cos ( ϕ 30 ) 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) b 1 sin ( ϕ 30 ) 1 4 ω a 2 γ 7 a 10 b 20 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) ,
J 82 = F 8 a 21 = 1 2 ( ω + σ 2 ) b 20 η 6 cos ( ϕ 40 ) 1 4 ω ( 3 α 1 + α ^ 6 ω 2 + 2 ω 4 γ 11 η 8 ω 4 2 + ω 2 + ( ω 4 2 ω 2 ) γ 15 η 8 2 ( ω 4 2 + ω 2 ) 2 ) a 20 ( 1 2 a 20 2 η 3 b 20 + 3 8 γ 1 b 20 1 4 a 20 2 γ 7 a 10 2 b 20 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 b 20 ) cos ( ϕ 40 ) ω + 1 8 γ 8 η 8 ω 4 2 + ω 2 b 20 sin ( ϕ 40 ) + 1 8 γ 2 b 20 sin ( ϕ 40 ) + 1 8 ω a 20 2 γ 7 a 10 2 b 20 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 2 ω a 20 2 ( ω + σ ) 2 f sin ( ϕ 20 ) ,
J 83 = F 8 b 11 = 1 2 ω γ 5 b 10 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 cos ( ϕ 30 ) 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 sin ( ϕ 30 ) 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ω γ 5 b 10 cos ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 84 = F 8 b 21 = 1 2 ( ω + σ 2 ) b 20 2 η 6 a 20 cos ( ϕ 40 ) 1 2 ω γ 3 b 20 ( 1 2 η 3 + 3 8 γ 1 a 20 2 + 1 4 γ 7 a 10 2 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 a 20 2 ) cos ( ϕ 40 ) ω a 20 1 4 ω γ 3 b 20 cos ( 2 ϕ 40 ) + 1 8 γ 8 η 8 ω 4 2 + ω 2 a 20 sin ( ϕ 40 ) + 1 8 γ 2 a 20 sin ( ϕ 40 ) 1 8 ω a 20 γ 7 a 10 2 cos ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) ,
J 85 = F 8 ϕ 11 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω a 20 γ 7 a 10 2 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) 1 4 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 86 = F 8 ϕ 21 = 1 4 ω ( α 2 α 5 ω 2 + ( ω 3 2 ω 2 ) γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 + ω 3 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 sin ( 2 ϕ 20 2 ϕ 10 ) 1 4 ω ( α 4 ω α 7 ω + 2 ω 3 ω γ 12 η 7 2 ( ω 3 2 + ω 2 ) 2 ω γ 13 η 7 ω 3 2 + ω 2 + ω 3 ω γ 14 η 7 ω 3 2 + ω 2 ) a 10 2 cos ( 2 ϕ 20 2 ϕ 10 ) + 1 4 ω a 2 γ 7 a 10 2 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) + 1 4 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 4 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) 1 2 ω a 20 ( ω + σ ) 2 f cos ( ϕ 20 ) ,
J 87 = F 8 ϕ 31 = 1 8 ω ( 2 γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( ϕ 30 ) 1 8 ( 2 γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( ϕ 30 ) 1 8 ω ( γ 4 + ω 3 γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 sin ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) + 1 8 ( γ 6 + γ 10 η 7 ω 3 2 + ω 2 ) a 10 b 10 cos ( 2 ϕ 20 2 ϕ 10 ϕ 30 ) 1 4 ω γ 5 b 10 2 sin ( 2 ϕ 20 2 ϕ 10 2 ϕ 30 ) ,
J 88 = F 8 ϕ 41 = ( 1 2 η 3 b 20 + 3 8 γ 1 a 20 2 b 20 + 1 4 γ 7 a 10 2 b 20 + 1 4 ω 4 γ 8 η 8 ω 4 2 + ω 2 a 20 2 b 20 ) sin ( ϕ 40 ) ω a 20 + 1 4 ω γ 3 b 20 2 sin ( 2 ϕ 40 ) + 1 8 γ 8 η 8 ω 4 2 + ω 2 a 20 b 20 cos ( ϕ 40 ) + 1 8 γ 2 a 20 b 20 cos ( ϕ 40 ) + 1 8 ω a 2 γ 7 a 10 2 b 20 sin ( 2 ϕ 20 2 ϕ 10 + ϕ 40 ) .

References

  1. Ji, J.C.; Yu, L.; Leung, A.Y.T. Bifurcation behavior of a rotor supported by active magnetic bearings. J. Sound Vib. 2000, 235, 133–151. [Google Scholar] [CrossRef]
  2. Saeed, N.A.; Awwad, E.M.; El-Meligy, M.A.; Nasr, E.S.A. Radial Versus Cartesian Control Strategies to Stabilize the Non-linear Whirling Motion of the Six-Pole Rotor-AMBs. IEEE Access 2020, 8, 138859–138883. [Google Scholar] [CrossRef]
  3. Ji, J.C.; Hansen, C.H. Non-linear oscillations of a rotor in active magnetic bearings. J. Sound Vib. 2001, 240, 599–612. [Google Scholar] [CrossRef]
  4. Ji, J.C.; Leung, A.Y.T. Non-linear oscillations of a rotor-magnetic bearing system under superharmonic resonance conditions. Int. J. Non-Linear Mech. 2003, 38, 829–835. [Google Scholar] [CrossRef]
  5. El-Shourbagy, S.M.; Saeed, N.A.; Kamel, M.; Raslan, K.R.; Abouel Nasr, E.; Awrejcewicz, J. On the Performance of a Non-linear Position-Velocity Controller to Stabilize Rotor-Active Magnetic-Bearings System. Symmetry 2021, 13, 2069. [Google Scholar] [CrossRef]
  6. Saeed, N.A.; Mahrous, E.; Abouel Nasr, E.; Awrejcewicz, J. Non-linear dynamics and motion bifurcations of the rotor active magnetic bearings system with a new control scheme and rub-impact force. Symmetry 2021, 13, 1502. [Google Scholar] [CrossRef]
  7. Zhang, W.; Zhan, X.P. Periodic and chaotic motions of a rotor-active magnetic bearing with quadratic and cubic terms and time-varying stiffness. Nonlinear Dyn. 2005, 41, 331–359. [Google Scholar] [CrossRef]
  8. Zhang, W.; Yao, M.H.; Zhan, X.P. Multi-pulse chaotic motions of a rotor-active magnetic bearing system with time-varying stiffness. Chaos Solitons Fractals 2006, 27, 175–186. [Google Scholar] [CrossRef]
  9. Zhang, W.; Zu, J.W.; Wang, F.X. Global bifurcations and chaos for a rotor-active magnetic bearing system with time-varying stiffness. Chaos Solitons Fractals 2008, 35, 586–608. [Google Scholar] [CrossRef]
  10. Zhang, W.; Zu, J.W. Transient and steady non-linear responses for a rotor-active magnetic bearings system with time-varying stiffness. Chaos Solitons Fractals 2008, 38, 1152–1167. [Google Scholar] [CrossRef]
  11. Li, J.; Tian, Y.; Zhang, W.; Miao, S.F. Bifurcation of multiple limit cycles for a rotor-active magnetic bearings system with time-varying stiffness. Int. J. Bifurc. Chaos 2008, 18, 755–778. [Google Scholar] [CrossRef]
  12. Li, J.; Tian, Y.; Zhang, W. Investigation of relation between singular points and number of limit cycles for a rotor–AMBs system. Chaos Solitons Fractals 2009, 39, 1627–1640. [Google Scholar] [CrossRef]
  13. El-Shourbagy, S.M.; Saeed, N.A.; Kamel, M.; Raslan, K.R.; Aboudaif, M.K.; Awrejcewicz, J. Control Performance, Stability Conditions, and Bifurcation Analysis of the Twelve-pole Active Magnetic Bearings System. Appl. Sci. 2021, 11, 10839. [Google Scholar] [CrossRef]
  14. Saeed, N.A.; Kandil, A. Two different control strategies for 16-pole rotor active magnetic bearings system with constant stiffness coefficients. Appl. Math. Model. 2021, 92, 1–22. [Google Scholar] [CrossRef]
  15. Wu, R.; Zhang, W.; Yao, M.H. Analysis of non-linear dynamics of a rotor-active magnetic bearing system with 16-pole legs. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Cleveland, OH, USA, 6–9 August 2017. [Google Scholar] [CrossRef]
  16. Wu, R.Q.; Zhang, W.; Yao, M.H. Non-linear dynamics near resonances of a rotor-active magnetic bearings system with 16-pole legs and time varying stiffness. Mech. Syst. Signal Process. 2018, 100, 113–134. [Google Scholar] [CrossRef]
  17. Zhang, W.; Wu, R.Q.; Siriguleng, B. Non-linear Vibrations of a Rotor-Active Magnetic Bearing System with 16-Pole Legs and Two Degrees of Freedom. Shock. Vib. 2020, 2020, 5282904. [Google Scholar]
  18. Ma, W.S.; Zhang, W.; Zhang, Y.F. Stability and multi-pulse jumping chaotic vibrations of a rotor-active magnetic bearing system with 16-pole legs under mechanical-electric-electro-magnetic excitations. Eur. J. Mech. A/Solids 2021, 85, 104120. [Google Scholar] [CrossRef]
  19. Ishida, Y.; Inoue, T. Vibration suppression of non-linear rotor systems using a dynamic damper. J. Vib. Control. 2007, 13, 1127–1143. [Google Scholar] [CrossRef]
  20. Saeed, N.A.; Awwad, E.M.; El-Meligy, M.A.; Nasr, E.S.A. Sensitivity analysis and vibration control of asymmetric non-linear rotating shaft system utilizing 4-pole AMBs as an actuator. Eur. J. Mech. A/Solids 2021, 86, 104145. [Google Scholar] [CrossRef]
  21. Saeed, N.A.; Eissa, M. Bifurcation analysis of a transversely cracked non-linear Jeffcott rotor system at different resonance cases. Int. J. Acoust. Vib. 2019, 24, 284–302. [Google Scholar] [CrossRef]
  22. Saeed, N.A.; Awwad, E.M.; El-Meligy, M.A.; Nasr, E.S.A. Analysis of the rub-impact forces between a controlled non-linear rotating shaft system and the electromagnet pole legs. Appl. Math. Model. 2021, 93, 792–810. [Google Scholar] [CrossRef]
  23. Srinivas, R.S.; Tiwari, R.; Kannababu, C. Application of active magnetic bearings in flexible rotordynamic systems—A state-of-the-art review. Mech. Syst. Signal Processing 2018, 106, 537–572. [Google Scholar] [CrossRef]
  24. Shan, J.; Liu, H.; Sun, D. Slewing and vibration control of a single-link flexible manipulator by positive position feedback (PPF). Mechatronics 2005, 15, 487–503. [Google Scholar] [CrossRef]
  25. Ahmed, B.; Pota, H.R. Dynamic compensation for control of a rotary wing UAV using positive position feedback. J. Intell. Robot. Syst. 2011, 61, 43–56. [Google Scholar] [CrossRef]
  26. Warminski, J.; Bochenski, M.; Jarzyna, W.; Filipek, P.; Augustyinak, M. Active suppression of non-linear composite beam vibrations by selected control algorithms. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 2237–2248. [Google Scholar] [CrossRef]
  27. Omidi, E.; Mahmoodi, S.N. Non-linear vibration suppression of flexible structures using non-linear modified positive position feedback approach. Nonlinear Dyn. 2015, 79, 835–849. [Google Scholar] [CrossRef]
  28. Saeed, N.A.; Kandil, A. Lateral vibration control and stabilization of the quasiperiodic oscillations for rotor-active magnetic bearings system. Nonlinear Dyn. 2019, 98, 1191–1218. [Google Scholar] [CrossRef]
  29. Diaz, I.M.; Pereira, E.; Reynolds, P. Integral resonant control scheme for cancelling human-induced vibrations in light-weight pedestrian structures. Struct. Control Health Monit 2012, 19, 55–69. [Google Scholar] [CrossRef]
  30. Al-Mamun, A.; Keikha, E.; Bhatia, C.S.; Lee, T.H. Integral resonant control for suppression of resonance in piezoelectric micro-actuator used in precision servomechanism. Mechatronics 2013, 23, 1–9. [Google Scholar] [CrossRef]
  31. Omidi, E.; Mahmoodi, S.N. Non-linear integral resonant controller for vibration reduction in non-linear systems. Acta Mech. Sin 2016, 32, 925–934. [Google Scholar] [CrossRef]
  32. MacLean, J.D.J.; Sumeet, S.A. A modified linear integral resonant controller for suppressing jump phenomenon and hysteresis in micro-cantilever beam structures. J. Sound Vib. 2020, 480, 115365. [Google Scholar] [CrossRef]
  33. Omidi, E.; Mahmoodi, S.N. Sensitivity analysis of the Non-linear Integral Positive Position Feedback and Integral Resonant controllers on vibration suppression of non-linear oscillatory systems. Commun. Nonlinear Sci. Numer. Simul. 2015, 22, 149–166. [Google Scholar] [CrossRef]
  34. Saeed, N.A.; El-Shourbagy, S.M.; Kamel, M.; Raslan, K.R.; Aboudaif, M.K. Nonlinear Dynamics and Static Bifurcations Control of the 12-Pole Magnetic Bearings System Utilizing the Integral Resonant Control strategy. J. Low Freq. Noise Vib. Act. Control 2022. [Google Scholar] [CrossRef]
  35. Saeed, N.A.; Moatimid, G.M.; Elsabaa, F.M.; Ellabban, Y.Y.; Elagan, S.K.; Mohamed, M.S. Time-Delayed Non-linear Integral Resonant Controller to Eliminate the Non-linear Oscillations of a Parametrically Excited System. IEEE Access 2021, 9, 74836–74854. [Google Scholar] [CrossRef]
  36. Saeed, N.A.; Mohamed, M.S.; Elagan, S.K.; Awrejcewicz, J. Integral Resonant Controller to Suppress the Non-linear Oscillations of a Two-Degree-of-Freedom Rotor Active Magnetic Bearing System. Processes 2022, 10, 271. [Google Scholar] [CrossRef]
  37. Ishida, Y.; Yamamoto, T. Linear and Non-Linear Rotordynamics: A Modern Treatment with Applications, 2nd ed.; Wiley-VCH Verlag GmbH & Co. KGaA: New York, NY, USA, 2012. [Google Scholar]
  38. Schweitzer, G.; Maslen, E.H. Magnetic Bearings: Theory, Design, and Application to Rotating Machinery; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  39. Nayfeh, A.H.; Mook, D.T. Non-Linear Oscillations; Wiley: New York, NY, USA, 1995. [Google Scholar]
  40. Nayfeh, A.H. Resolving Controversies in the Application of the Method of Multiple Scales and the Generalized Method of Averaging. Nonlinear Dyn. 2005, 40, 61–102. [Google Scholar] [CrossRef]
  41. Vlase, S.; Năstac, C.; Marin, M.; Mihălcică, M. A method for the study of the vibration of mechanical bars systems with symmetries. Acta Tech. Napoc. 2017, 60, 539–544. [Google Scholar]
  42. Slotine, J.-J.E.; Li, W. Applied Non-Linear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  43. Yang, W.Y.; Cao, W.; Chung, T.; Morris, J. Applied Numerical Methods Using Matlab; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2005. [Google Scholar]
Figure 1. (a) twelve-poles rotor system at its nominal position, (b) twelve-poles rotor system with small displacements x ( t ) and y ( t ) in the horizontal and vertical directions, respectively.
Figure 1. (a) twelve-poles rotor system at its nominal position, (b) twelve-poles rotor system with small displacements x ( t ) and y ( t ) in the horizontal and vertical directions, respectively.
Applsci 12 08300 g001
Figure 2. The engineering implementation of the combined control algorithm (i.e., P D + I R C + P P F controller).
Figure 2. The engineering implementation of the combined control algorithm (i.e., P D + I R C + P P F controller).
Applsci 12 08300 g002
Figure 3. Vibration amplitudes and phase angles of the twelve-poles rotor in the case of the PD-control algorithm only: (a,b) vibration amplitudes ( a 1 ,   a 2 ) when f = 0.004 ,   0.007 ,   0.01 , and 0.013 , (c) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.004 , (d) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.013 .
Figure 3. Vibration amplitudes and phase angles of the twelve-poles rotor in the case of the PD-control algorithm only: (a,b) vibration amplitudes ( a 1 ,   a 2 ) when f = 0.004 ,   0.007 ,   0.01 , and 0.013 , (c) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.004 , (d) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.013 .
Applsci 12 08300 g003
Figure 4. Vibration amplitudes of the twelve-poles rotor and the P P F controller in the case of the P D + P P F -control algorithm when f = 0.0025 ,   0.0075 ,   0.0125 , and 0.0175 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.0025 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.0175 .
Figure 4. Vibration amplitudes of the twelve-poles rotor and the P P F controller in the case of the P D + P P F -control algorithm when f = 0.0025 ,   0.0075 ,   0.0125 , and 0.0175 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.0025 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.0175 .
Applsci 12 08300 g004aApplsci 12 08300 g004b
Figure 5. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C -control algorithm when f = 0.02 ,   0.04 ,   0.06 , and 0.08 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.02 , and (d) phase angles when f = 0.08 .
Figure 5. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C -control algorithm when f = 0.02 ,   0.04 ,   0.06 , and 0.08 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.02 , and (d) phase angles when f = 0.08 .
Applsci 12 08300 g005
Figure 6. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm when f = 0.025 ,   0.045 ,   0.065 , and 0.085 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.025 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.085 .
Figure 6. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm when f = 0.025 ,   0.045 ,   0.065 , and 0.085 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.025 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when f = 0.085 .
Applsci 12 08300 g006aApplsci 12 08300 g006b
Figure 7. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 ,   η 1 = η 2 = η 3 = η 4 = 0.2 when p = 1.2 ,   1.6 , and 2.0 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Figure 7. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 ,   η 1 = η 2 = η 3 = η 4 = 0.2 when p = 1.2 ,   1.6 , and 2.0 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Applsci 12 08300 g007
Figure 8. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 ,   P = 1.5 ,   η 2 = η 4 = 0.2 when η 1 = η 3 = 0.1 ,   0.3 , and 0.5 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Figure 8. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 ,   P = 1.5 ,   η 2 = η 4 = 0.2 when η 1 = η 3 = 0.1 ,   0.3 , and 0.5 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Applsci 12 08300 g008
Figure 9. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 ,   η 1 = η 3 = 0.2 , when η 2 = η 4 = 0.1 ,   0.3 , and 0.5 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Figure 9. Vibration amplitudes of the twelve-poles rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 ,   η 1 = η 3 = 0.2 , when η 2 = η 4 = 0.1 ,   0.3 , and 0.5 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Applsci 12 08300 g009
Figure 10. Vibration amplitudes of the twelve-poles s rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 , and η 1 = η 2 = η 3 = η 4 = 0.2 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Figure 10. Vibration amplitudes of the twelve-poles s rotor in the case of the P D + I R C + P P F -control algorithm at f = 0.013 , and η 1 = η 2 = η 3 = η 4 = 0.2 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller.
Applsci 12 08300 g010aApplsci 12 08300 g010b
Figure 11. Vibration amplitudes and phase angles of the twelve-poles rotor in the case of both the P D -control only and the P D + P P F -control algorithms, when f = 0.013 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F -controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 2 = η 3 = η 4 = 0.0 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 3 = 0.2 ,   η 2 = η 4 = 0.0 .
Figure 11. Vibration amplitudes and phase angles of the twelve-poles rotor in the case of both the P D -control only and the P D + P P F -control algorithms, when f = 0.013 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F -controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 2 = η 3 = η 4 = 0.0 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 3 = 0.2 ,   η 2 = η 4 = 0.0 .
Applsci 12 08300 g011aApplsci 12 08300 g011b
Figure 12. Time response of the rotor system according to Figure 11 when σ = 0 (i.e., when Ω = ω ) in the case of the P D -control algorithm and the P D + P P F -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, and (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the PPF controller.
Figure 12. Time response of the rotor system according to Figure 11 when σ = 0 (i.e., when Ω = ω ) in the case of the P D -control algorithm and the P D + P P F -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, and (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the PPF controller.
Applsci 12 08300 g012aApplsci 12 08300 g012b
Figure 13. Time response of the rotor system according to Figure 11 when σ = 0.1 (i.e., when Ω = ω + 0.1 ) in the case of the P D -control algorithm and the P D + P P F -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, and (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the PPF controller.
Figure 13. Time response of the rotor system according to Figure 11 when σ = 0.1 (i.e., when Ω = ω + 0.1 ) in the case of the P D -control algorithm and the P D + P P F -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, and (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the PPF controller.
Applsci 12 08300 g013aApplsci 12 08300 g013b
Figure 14. Vibration amplitudes and phase-angles of the twelve-poles rotor in the case of both the P D -control only and the P D + I R C -control algorithm when f = 0.013 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c) phase angles ( a 1 ,   a 2 ) when η 1 = η 2 = η 3 = η 4 = 0.0 , and (d) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 3 = 0.0 ,   η 2 = η 4 = 0.2 .
Figure 14. Vibration amplitudes and phase-angles of the twelve-poles rotor in the case of both the P D -control only and the P D + I R C -control algorithm when f = 0.013 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, (c) phase angles ( a 1 ,   a 2 ) when η 1 = η 2 = η 3 = η 4 = 0.0 , and (d) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 3 = 0.0 ,   η 2 = η 4 = 0.2 .
Applsci 12 08300 g014aApplsci 12 08300 g014b
Figure 15. Time response of the rotor system according to Figure 14 when σ = 0.0 (i.e., when Ω = ω ) in the case of the P D -control algorithm and the P D + I R C -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, and (e,f) the temporal oscillations z 5 ( t ) and z 6 ( t ) of the I R C controller.
Figure 15. Time response of the rotor system according to Figure 14 when σ = 0.0 (i.e., when Ω = ω ) in the case of the P D -control algorithm and the P D + I R C -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, and (e,f) the temporal oscillations z 5 ( t ) and z 6 ( t ) of the I R C controller.
Applsci 12 08300 g015aApplsci 12 08300 g015b
Figure 16. Time response of the rotor system according to Figure 14 when σ = 0.1 (i.e., when Ω = ω + 0.1 ) in the case of the P D -control algorithm and the P D + I R C -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, (e,f) the temporal oscillations z 5 ( t ) , and z 6 ( t ) of the I R C controller.
Figure 16. Time response of the rotor system according to Figure 14 when σ = 0.1 (i.e., when Ω = ω + 0.1 ) in the case of the P D -control algorithm and the P D + I R C -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, (e,f) the temporal oscillations z 5 ( t ) , and z 6 ( t ) of the I R C controller.
Applsci 12 08300 g016aApplsci 12 08300 g016b
Figure 17. Vibration amplitudes and phase angles of the twelve-poles rotor in the case of both the P D -control algorithm only and the P D + I R C + P P F -control algorithm when f = 0.013 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 2 = η 3 = η 4 = 0.0 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 3 =   η 2 = η 4 = 0.2 .
Figure 17. Vibration amplitudes and phase angles of the twelve-poles rotor in the case of both the P D -control algorithm only and the P D + I R C + P P F -control algorithm when f = 0.013 : (a,b) vibration amplitudes ( a 1 ,   a 2 ) of the rotor, and (c,d) vibration amplitudes ( b 1 ,   b 2 ) of the P P F controller, (e) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 2 = η 3 = η 4 = 0.0 , and (f) phase angles ( ϕ 1 ,   ϕ 2 ) when η 1 = η 3 =   η 2 = η 4 = 0.2 .
Applsci 12 08300 g017
Figure 18. Time response of the rotor system according to Figure 17 when σ = 0.0 (i.e., when Ω = ω ) in the case of the P D -control algorithm and the P D + I R C + P P F -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the P P F controller, and (g,h) the temporal oscillations z 5 ( t ) and z 6 ( t ) of the I R C controller.
Figure 18. Time response of the rotor system according to Figure 17 when σ = 0.0 (i.e., when Ω = ω ) in the case of the P D -control algorithm and the P D + I R C + P P F -control algorithm: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the P P F controller, and (g,h) the temporal oscillations z 5 ( t ) and z 6 ( t ) of the I R C controller.
Applsci 12 08300 g018aApplsci 12 08300 g018b
Figure 19. Time response of the rotor system according to Figure 17 when σ = 0.1 (i.e., when Ω = ω + 0.1 ) when the P D -control algorithm and the P D + I R C + P P F -control algorithm are applied: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the P P F controller, and (g,h) the temporal oscillations z 5 ( t ) and z 6 ( t ) of the I R C controller.
Figure 19. Time response of the rotor system according to Figure 17 when σ = 0.1 (i.e., when Ω = ω + 0.1 ) when the P D -control algorithm and the P D + I R C + P P F -control algorithm are applied: (a,b) the temporal oscillations z 1 ( t ) and z 2 ( t ) of the rotor system, (c) the rotor whirling orbits, (d) the rotor frequency spectrum, (e,f) the temporal oscillations z 3 ( t ) and z 4 ( t ) of the P P F controller, and (g,h) the temporal oscillations z 5 ( t ) and z 6 ( t ) of the I R C controller.
Applsci 12 08300 g019aApplsci 12 08300 g019b
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Saeed, N.A.; El-Shourbagy, S.M.; Kamel, M.; Raslan, K.R.; Awrejcewicz, J.; Gepreel, K.A. On the Resonant Vibrations Control of the Nonlinear Rotor Active Magnetic Bearing Systems. Appl. Sci. 2022, 12, 8300. https://0-doi-org.brum.beds.ac.uk/10.3390/app12168300

AMA Style

Saeed NA, El-Shourbagy SM, Kamel M, Raslan KR, Awrejcewicz J, Gepreel KA. On the Resonant Vibrations Control of the Nonlinear Rotor Active Magnetic Bearing Systems. Applied Sciences. 2022; 12(16):8300. https://0-doi-org.brum.beds.ac.uk/10.3390/app12168300

Chicago/Turabian Style

Saeed, Nasser A., Sabry M. El-Shourbagy, Magdi Kamel, Kamal R. Raslan, Jan Awrejcewicz, and Khaled A. Gepreel. 2022. "On the Resonant Vibrations Control of the Nonlinear Rotor Active Magnetic Bearing Systems" Applied Sciences 12, no. 16: 8300. https://0-doi-org.brum.beds.ac.uk/10.3390/app12168300

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop