Next Article in Journal
Elemental Profiles of Wild Thymus L. Plants Growing in Different Soil and Climate Conditions
Next Article in Special Issue
Dynamic Analysis of a Wiper Blade in Consideration of Attack Angle and Clarification of the Jumping Phenomenon
Previous Article in Journal
Experimental Study of Drag Torque between Wet Clutch Discs
Previous Article in Special Issue
The Stability Analysis of a Vibrating Auto-Parametric Dynamical System Near Resonance
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

ALIPPF-Controller to Stabilize the Unstable Motion and Eliminate the Non-Linear Oscillations of the Rotor Electro-Magnetic Suspension System

by
Nasser A. Saeed
1,*,
Jan Awrejcewicz
2,
Abd Allah A. Mousa
3 and
Mohamed S. Mohamed
3
1
Department of Physics and Engineering Mathematics, Faculty of Electronic Engineering, Menoufia University, Menouf 32952, Egypt
2
Department of Automation, Biomechanics, and Mechatronics, Faculty of Mechanical Engineering, Lodz University of Technology, 90924 Lodz, Poland
3
Department of Mathematics and Statistics, College of Science, Taif University, P.O. Box 11099, Taif 21944, Saudi Arabia
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(8), 3902; https://doi.org/10.3390/app12083902
Submission received: 7 March 2022 / Revised: 8 April 2022 / Accepted: 10 April 2022 / Published: 12 April 2022
(This article belongs to the Special Issue Nonlinear Vibrations)

Abstract

:
Within this work, an advanced control algorithm was proposed to eliminate the non-linear vibrations of the rotor electro-magnetic suspension system. The suggested control algorithm is known as the Adaptive Linear Integral Positive Position Feedback controller (ALIPPF-controller). The ALIPPF-controller is a combination of first-order and second-order filters that are coupled linearly to the targeted non-linear system in order to absorb the excessive vibratory energy. According to the introduced control strategy, the dynamical model of the controlled rotor system was established as six non-linear differential equations that are coupled linearly. The obtained dynamical model was investigated analytically applying the asymptotic analysis, where the slow-flow equations were extracted. Based on the derived slow-flow equations, the bifurcation behaviors of the controlled system were explored by plotting the different bifurcation diagrams. In addition, the performance of the ALIPPF-controller in eliminating the rotor lateral vibrations was compared with the conventional Positive Position Feedback (PPF) controller. The acquired results illustrated that the ALIPPF-controller is the best control technique that can eliminate the considered system’s lateral vibrations regardless of the angular speed and eccentricity of the rotating shaft. Finally, to demonstrate the accuracy of the obtained analytical results, numerical validation was performed for all obtained bifurcation diagrams that were in excellent agreement with the analytical solutions.

1. Introduction

The rotor electro-magnetic suspension system has been introduced as an optimal solution to get rid of the drawbacks of the conventional bearings system such as the limited operational speed, the existence of friction between the rotor and the bearings, and the need for continuous lubrication, etc. The rotor electro-magnetic suspension system is a type of bearing that supports the rotating shafts using an electro-magnetic levitation mechanism. The main working principle of the electro-magnetic levitation is to apply controllable magnetic forces that support the rotor system in its hovering position without physical contact with the stator. The operation of the rotor with the contactless feature has earned the electro-magnetic suspension system many advantages over the other conventional bearing mechanisms such as the lack of friction between the moving parts and the stator, high-speed operation, less maintenance, clean working environment with no need for lubrication, high reliability, high durability, etc. Accordingly, the rotor electro-magnetic suspension system was and still is one of the main research subjects for scientists and engineers worldwide, where many researchers have focused their work on investigating the dynamical characteristics of the different electro-magnetic suspension configurations. On the other hand, some researchers have introduced advanced control algorithms to enhance the dynamical behaviors of electro-magnetic suspension systems regardless of their configurations. Ji et al. [1] studied non-linear dynamics and the bifurcation behaviors of the four-pole electro-magnetic suspension system. They employed the asymptotic perturbation method along with the normal form technique to analyze the considered system. The authors reported that the four poles electro-magnetic suspension system may lose its stability via different types of bifurcations such as the Hopf and saddle node. Saeed et al. [2,3] discussed the vibration control and bifurcation elimination of the six-pole electro-magnetic suspension system utilizing the conventional PD-control technique. The authors studied two different configurations for the introduced controller. The first PD-control strategy has been designed based on the lateral displacements and lateral velocities of the rotor system, while the second control methodology is implemented based on the rotor radial oscillations. According to the studied control techniques, they concluded that the first control method is more efficient in mitigating the system lateral vibrations, while the second technique is more robust against the system instability. Ji et al. [4,5] investigated the non-linear oscillations of the eight-pole electro-magnetic suspension system controlled via the PD-control algorithm at the primary and super-harmonic resonances. The authors illustrated that the system has both complex bifurcation behavior and multiple fixed-point solutions. El-Shourbagy et al. [6] introduced a non-linear position-velocity control algorithm to enhance the dynamical characteristics and mitigate the vibratory motion of the eight-pole suspension system. Saeed et al. [7] explored the dynamics of the eight-poles electro-magnetic suspension system when the rub and impact occur between the moving parts and the rotor housing, where the obtained results showed that the rotor system may perform periodic, quasiperiodic, or chaotic motion depending on the impact stiffness magnitude and dynamic friction coefficient. Zhang et al. [8,9,10,11,12,13] extensively investigated the eight-pole suspension system with time-varying stiffness coefficients, where many non-linear phenomena have been reported such as the multi-pulse chaotic motion, saddle-node, and Hopf bifurcations. El-Shourbagy et al. [14] studied the 12-pole electro-magnetic suspension system utilizing the PD-controller. They investigated the effect of the different control parameters on the rotor stability charts. The acquired results showed that the system stability can be controlled using the proportional and derivative gain of the applied controller. Saeed and Kandil [15] introduced both the radial and cartesian proportional-derivative control methods to stabilize the 16-pole electro-magnetic suspension system. They derived the equations of motion according to the suggested control methodologies. In addition, perturbation analysis was employed to explore the system’s steady-state oscillation. Based on the obtained results, the authors confirmed that the radial control technique is more robust than the Cartesian one. In addition, extensive investigations of the 16-pole electro-magnetic suspension system with time-varying stiffness have been presented by Zhang and his coworkers [16,17,18,19,20]. Due to the flexibility of the electro-magnetic suspension system in adjusting its dynamical behavior according to the applied control algorithm, many researchers employed this suspension system as an active actuator to reshape the dynamics of some rotating machinery [21,22,23,24]. Ishida and Inoue [21] implemented a non-linear magnetic absorber based on the push–pull control mechanism to mitigate the lateral vibrations of a non-linear Jeffcott rotor system. Saeed et al. [22,23,24] applied a non-linear position-velocity control algorithm to enhance the dynamical characteristics of an asymmetric rotor system utilizing the four-pole electro-magnetic suspension system as an actuator.
One of the feasible control techniques that are applied extensively to suppress the non-linear oscillations of different dynamical systems is known as Positive Position Feedback (PPF) controller [25,26,27,28,29]. The PPF-controller is a second-order filter that couples to the targeted vibratory system linearly to suppress its undesired oscillation. Saeed et al. [29] applied the PPF-controller to mitigate the non-linear oscillation of the eight-pole electro-magnetic suspension system. The authors reported that the PPF-controller can eliminate the lateral oscillation of the rotor system at the perfect resonance conditions (i.e., when the rotor angular velocity is equal to the natural frequency of the suspension system); however, the PPF-controller can add more excessive vibratory energy to the rotor system rather than suppress it, if the perfect resonance conditions have been lost. Besides the PPF-controller, there is another efficient control algorithm known as the Linear Integral Resonant (LIR) controller that has been utilized to dampen the non-linear vibration of many dynamical systems [30,31,32,33,34,35,36,37,38]. The LIR-controller is a first-order filter that has been coupled to the targeted system to modify its linear damping coefficients, which ultimately mitigates the undesired oscillation. Saeed et al. [38] studied the vibration control of the electro-magnetic suspension system utilizing the LIR-controller. The authors reported that the integration of the LIR-controller to the rotor system can eliminate the non-linear bifurcations and mitigate the system vibration to respond as a linear system.
With this work, a modified version of the control algorithms that have the advantage of both the PPF-controller and LIR-controller has been investigated for the first time. The proposed control strategy is a combination of first-order and second-order filters that have been couped to the considered suspension system to reduce the undesired lateral oscillation. The main task of the second-order filter is to behave as a PPF-controller to eliminate the catastrophic oscillations of the considered system at the perfect resonance conditions, while the essential role of the first-order filter is to work as an LIR-controller that increases the controlled system linear damping, which ultimately removes the drawback of the PPF-controller when the resonance condition is lost. According to this philosophy, the whole system mathematical model was derived. Then, the established dynamical model was investigated analytically using asymptotic techniques [39,40,41,42,43,44,45]. The performance of the proposed control method in eliminating the non-linear oscillations of the considered suspension system was explored by plotting the different response curves. The obtained analytical and numerical results showed that the introduced control strategy can eliminate the rotor vibrations regardless of both the angular speed and eccentricity of the rotor system.

2. Equations of Motion

The considered rotor system is illustrated in Figure 1, where the rotating shaft is supported via eight electro-magnetic poles that serve as active magnetic bearings. The rotor system is assumed to be a rigid body with mass m and eccentricity e and rotates with angular speed υ . Accordingly, the dynamical equations that govern the system’s lateral vibrations in X 1 and X 2 directions can be written as follows [46,47]:
m x ¨ 1 R X 1 = m e   υ 2 cos ( υ t )
m x ¨ 2 R X 2 = m e   υ 2 sin ( υ t )
where R X 1 and R X 2 denote the net restoring forces that are exerted on the rotating shaft due to the attractive electro-magnetic forces h j ( j = 1 , 2 , , 8 ). It is considered that the eight electro-magnetic poles are symmetric. Accordingly, the attractive magnetic force that is exerted on the rotor system due to each pole can be expressed as follows [47]:
h j = 1 4 μ 0 N 2 A cos ( θ ) I j 2 g j 2 = Δ I j 2 g j 2 , j = 1 8
where Δ = 1 4 μ 0 N 2 A cos ( θ ) is constant (i.e., μ 0 is the air permeability, N is the winding number of the electrical coil of each pole, A cos ( θ ) is the effective cross-sectional area of each pole), I j is the electrical current in the j pole, and g j is the air-gap size between the rotating shaft and the j th pole. For the small displacements x 1 ( t ) and x 2 ( t ) of the rotating shaft within the pole housing, the instantaneous air-gap size g j between the rotating shaft and the j pole can be expressed as follows:
g 1 ( x 1 , x 2 ) = s 0 + x 1 sin ( α ) x 2 cos ( α ) , g 2 ( x 1 , x 2 ) = s 0 + x 1 cos ( α ) x 2 sin ( α ) g 3 ( x 1 , x 2 ) = s 0 + x 1 cos ( α ) + x 2 sin ( α ) , g 4 ( x 1 , x 2 ) = s 0 + x 1 sin ( α ) + x 2 cos ( α ) g 5 ( x 1 , x 2 ) = s 0 x 1 sin ( α ) + x 2 cos ( α ) , g 6 ( x 1 , x 2 ) = s 0 x 1 cos ( α ) + x 2 sin ( α ) g 7 ( x 1 , x 2 ) = s 0 x 1 cos ( α ) x 2 sin ( α ) , g 8 ( x 1 , x 2 ) = s 0 x 1 sin ( α ) x 2 cos ( α ) }
where s 0 is the radial clearance between the rotating shaft and each electro-magnetic pole when the geometric center of both the rotor and pole housing coincided as shown in Figure 1a, and 2 α = 45 ° is the angle between every two consecutive electro-magnetic poles. Based on the system geometry shown in Figure 1a,b, the electrical currents of the eight-pole are designed based on the push–pull control mechanism as follows:
I 1 = I 8 = I 0 i X 2 , I 2 = I 3 = I 0 + i X 1 , I 4 = I 5 = I 0 + i X 2 , I 6 = I 7 = I 0 i X 1
where I 0 is a per-magnetizing constant current, i X 1 , and i X 2 are the control currents in X 1 and X 2 directions, respectively.
Within this article, a new control algorithm was proposed to eliminate the non-linear oscillations of the considered system beside the conventional PD-controller. The suggested new control technique is defined as the Linear Integral Positive Position Feedback controller (LIPPF-controller). Therefore, according to the introduced control strategy, the control currents i X 1 and i X 2 are designed as follows:
i X 1 = k 1 x 1 + k 2 x ˙ 1 + k 3 x 3 + k 4 x 5 , i X 2 = k 1 x 2 + k 2 x ˙ 2 + k 5 x 4 + k 6 x 6
where k 1 and k 2 are the control gains of the conventional PD-controller. k 3 and k 4 are the control signal gains of the LIPPF-controller that is connected to the system oscillation mode in the X 1 direction. k 5 and k 6 are the control signal gains of the LIPPF-controller that is connected to the system oscillation mode in the X 2 direction.
The LIPPF-controller is a combination of both first-order and second-order filters that are excited linearly by the rotor system displacement, where the equations of motion of the proposed LIPPF-controller are given as follows:
x ¨ 3 + c 1 x ˙ 3 + λ 1 x 3 = δ 1 x 1
x ¨ 4 + c 2 x ˙ 4 + λ 2 x 4 = δ 2 x 2
x ˙ 5 + λ 3 x 5 = δ 3 x 1
x ˙ 6 + λ 4 x 6 = δ 4 x 2
where c 1 ,   c 2 represent the linear damping coefficients, λ 1 ,   λ 2 denote the natural frequency, and λ 3 ,   λ 4 are the internal loop feedback gain of the suggested controller (i.e., LIPPF-controller). In addition, δ 1 ,   δ 2 ,   δ 3 , and δ 4 represent the feedback signal gains of the LIPPF-controller. The engineering implementation of the suggested control method is illustrated in Figure 2, where the instantaneous Cartesian displacements x 1 ( t ) and x 2 ( t ) of the rotor system are measured using proximity sensors. The measured signals are fed into a digital computer to manipulate them according to the proposed control strategy to form the control currents i X 1 = k 1 x 1 + k 2 x ˙ 1 + k 3 x 3 + k 4 x 5 and i X 2 = k 1 x 2 + k 2 x ˙ 2 + k 5 x 4 + k 6 x 6 . Then, the electrical currents of the eight-pole (i.e., I 1 ,   I 2 , I 8 ) are applied via a power amplifier network as illustrated in Figure 2 in detail. Now to obtain the full mathematical model that governs the rotor and the controller dynamics, let us substitute Equations (4)–(6) into Equation (3), we get the restoring forces of the eight-pole as follows:
h 1 = Δ ( I 0 k 1 x 2 k 2 x ˙ 2 k 5 x 4 k 6 x 6 s 0 + x 1 sin ( α ) x 2 cos ( α ) ) 2
h 2 = Δ ( I 0 + k 1 x 1 + k 2 x ˙ 1 + k 3 x 3 + k 4 x 5 s 0 + x 1 cos ( α ) x 2 sin ( α ) ) 2
h 3 = Δ ( I 0 + k 1 x 1 + k 2 x ˙ 1 + k 3 x 3 + k 4 x 5 s 0 + x 1 cos ( α ) + x 2 sin ( α ) ) 2
h 4 = Δ ( I 0 + k 1 x 2 + k 2 x ˙ 2 + k 5 x 4 + k 6 x 6 s 0 + x 1 sin ( α ) + x 2 cos ( α ) ) 2
h 5 = Δ ( I 0 + k 1 x 2 + k 2 x ˙ 2 + k 5 x 4 + k 6 x 6 s 0 x 1 sin ( α ) + x 2 cos ( α ) ) 2
h 6 = Δ ( I 0 k 1 x 1 k 2 x ˙ 1 k 3 x 3 k 4 x 5 s 0 x 1 cos ( α ) + x 2 sin ( α ) ) 2
h 7 = Δ ( I 0 k 1 x 1 k 2 x ˙ 1 k 3 x 3 k 4 x 5 s 0 x 1 cos ( α ) x 2 sin ( α ) ) 2
h 8 = Δ ( I 0 k 1 x 2 k 2 x ˙ 2 k 5 x 4 k 6 x 6 s 0 x 1 sin ( α ) x 2 cos ( α ) ) 2
According to the orientations of the eight attractive forces h j ( j = 1 , 2 , , 8 ) shown in Figure 1, the net restoring forces R X 1 and R X 2 in X 1 and X 2 directions can be expressed as:
R X 1 = ( h 6 + h 7 h 2 h 3 ) cos ( α ) + ( h 5 + h 8 h 1 h 4 ) sin ( α )
R X 2 = ( h 1 + h 8 h 4 h 5 ) cos ( α ) + ( h 2 + h 7 h 3 h 6 ) sin ( α )
Substituting Equations (11)–(20) into Equations (1) and (2), with introducing the dimensionless variables and parameters: y 1 = x 1 s 0 ,   y 2 = x 2 s 0 ,   y 3 = x 3 s 0 ,   y 4 = x 4 s 0 ,   y 5 = x 5 s 0 ,   y 6 = x 6 s 0 ,   τ = ω n t ,   ω n = Δ / m s 0 3 ,   p = s 0 I 0 k 1 ,   d = c 0 ω n I 0 k 2 ,   μ 1 = c 1 2 ω n ,   μ 2 = c 2 2 ω n ,   ω 1 = λ 1 ω n 2 ,   ω 2 = λ 2 ω n 2 ,   ρ 1 = λ 3 ω n ,   ρ 2 = λ 4 ω n ,   η 1 = 8 s 0 I 0 cos ( α ) k 3 ,   η 2 = 8 s 0 I 0 cos ( α ) k 4 ,   η 3 = 8 s 0 I 0 cos ( α ) k 5 ,   η 4 = 8 s 0 I 0 cos ( α ) k 6 ,   η 5 = δ 1 ω n 2 ,   η 6 = δ 2 ω n 2 ,   η 7 = δ 3 ω n ,   η 8 = η 4 ω n ,   f = e s 0 ,   Ω = υ ω n , we have from Equations (1), (2) and (7)–(10) the following dimensionless equations of motion:
y ¨ 1 + 2 μ y ˙ 1 + ω 2 y 1 ( α 1 y 1 3 + α 2 y 1 y 2 2 + α 3 y 1 2 y ˙ 1 + α 4 y ˙ 1 y 2 2 + α 5 y 1 y ˙ 2 2 + α 6 y 1 y ˙ 1 2 + α 7 y 1 y 2 y ˙ 2 + β 1 y 1 2 y 3 + β 2 y 1 y ˙ 1 y 3 + β 3 y 1 y 3 2 + β 4 y 1 y 2 y 4 + β 5 y 1 y 4 2 + β 6 y 1 y ˙ 2 y 4 + β 7 y 2 2 y 3 + β 8 y 1 y 3 y 5 + β 9 y 1 y ˙ 1 y 5 + β 10 y 1 y 4 y 6 + β 11 y 1 2 y 5 + β 12 y 1 y 6 2 + β 13 y 1 y 2 y 6 + β 14 y 1 y ˙ 2 y 6 + β 15 y 1 y 5 2 + β 16 y 2 2 y 5 ) = f Ω 2 cos ( Ω τ ) + η 1 y 3 + η 2 y 5
y ¨ 2 + 2 μ y ˙ 2 + ω 2 y 2 ( α 1 y 2 3 + α 2 y 2 y 1 2 + α 3 y 2 2 y ˙ 2 + α 4 y ˙ 2 y 1 2 + α 5 y 2 y ˙ 1 2 + α 6 y 2 y ˙ 2 2 + α 7 y 2 y 1 y ˙ 1 + γ 1 y 2 2 y 4 + γ 2 y 2 y ˙ 2 y 4 + γ 3 y 2 y 4 2 + γ 4 y 2 y 1 y 3 + γ 5 y 2 y 3 2 + γ 6 y 2 y ˙ 1 y 3 + γ 7 y 4 y 1 2 + γ 8 y 2 y 4 y 6 + γ 9 y 2 y ˙ 2 y 6 + γ 10 y 2 y 3 y 5 + γ 11 y 2 2 y 6 + γ 12 y 2 y 5 2 + γ 13 y 2 y 1 y 5 + γ 14 y 2 y ˙ 1 y 5 + γ 15 y 2 y 6 2 + γ 16 y 1 2 y 6 ) = f Ω 2 sin ( Ω τ ) + η 3 y 4 + η 4 y 6
y ¨ 3 + 2 μ 1 y ˙ 3 + ω 1 2 y 3 = η 5 y 1
y ¨ 4 + 2 μ 2 y ˙ 4 + ω 2 2 y 4 = η 6 y 2
y ˙ 5 + ρ 1 y 5 = η 7 y 1
y ˙ 6 + ρ 2 y 6 = η 8 y 2
where Equations (21) and (22) represent the equations of motion of the controller rotor system, while Equations (23)–(26) are the equations of the coupled LIPPF-controller. The coefficients of Equations (21)–(26) are given in the Appendix A.

3. Analytical Investigations

To evaluate the performance of the suggested control method (i.e., LIPPF-controller) in eliminating the non-linear oscillations of the studied rotor-magnetic bearings system, an analytical solution for Equations (21)–(26) is sought by applying the perturbation analysis as follows [39,40]:
y 1 ( τ , ε ) = y 10 ( T 0 , T 1 ) + ε y 11 ( T 0 , T 1 ) + O ( ε 2 )
y 2 ( τ , ε ) = y 20 ( T 0 , T 1 ) + ε y 21 ( T 0 , T 1 ) + O ( ε 2 )
y 3 ( τ , ε ) = y 30 ( T 0 , T 1 ) + ε y 31 ( T 0 , T 1 ) + O ( ε 2 )
y 4 ( τ , ε ) = y 40 ( T 0 , T 1 ) + ε y 41 ( T 0 , T 1 ) + O ( ε 2 )
y 5 ( τ , ε ) = ε y 50 ( T 0 , T 1 ) + ε 2 y 51 ( T 0 , T 1 ) + O ( ε 2 )
y 6 ( τ , ε ) = ε y 60 ( T 0 , T 1 ) + ε 2 y 61 ( T 0 , T 1 ) + O ( ε 2 )
where T 0 = τ , T 1 = ε τ , and ε is called a perturbation parameter used as a book-keeping parameter only through this work [40]. Accordingly, the ordinary derivatives d d τ and d 2 d τ 2 can be expressed in terms of T 0 and T 1 as follows:
d d τ = D 0 + ε D 1 ,   d 2 d τ 2 = D 0 2 + 2 ε D 0 D 1 ,   D j = T j , j = 0 , 1
To perform the perturbation procedure, the parameters of Equations (21)–(26) should be scaled as follows:
f = ε f ˜ , μ = ε μ ˜ , μ 1 = ε μ ˜ 1 , μ 2 = ε μ ˜ 2 , α j = ε α ˜ j , β j = ε β ˜ j , γ j = ε γ ˜ j , η k = ε η ˜ k , j = 1 , 2 , , 7 ; k = 1 , 3 , 5 , 6 , 7 , 8 .
Substituting Equations (27)–(34) into Equations (21)–(26) with expanding the resulting equations and comparing the coefficients of the same power of ε , yields O ( ε 0 ) :
( D 0 2 + ω 2 ) y 10 = 0
( D 0 2 + ω 2 ) y 20 = 0
( D 0 2 + ω 1 2 ) y 30 = 0
( D 0 2 + ω 2 2 ) y 40 = 0
O ( ε 1 ) :
( D 0 2 + ω 2 ) y 11 = 2 D 0 D 1 y 10 2 μ ˜ D 0 y 10 + α ˜ 1 y 10 3 + α ˜ 2 y 10 y 20 2 + α ˜ 3 y 10 2 D 0 y 10 + α ˜ 4 y 20 2 D 0 y 10 + α ˜ 5 y 10 ( D 0 y 20 ) 2 + α ˜ 6 y 10 ( D 0 y 10 ) 2 + α ˜ 7 y 10 y 20 D 0 y 20 + β ˜ 1 y 10 2 y 30 + β ˜ 2 y 10 D 0 y 10 y 30 + β ˜ 3 y 10 y 30 2 + β ˜ 4 y 10 y 20 y 40 + β ˜ 5 y 10 y 40 2 + β ˜ 6 y 10 D 0 y 20 y 40 + β ˜ 7 y 20 2 y 30 + β 8 y 10 y 30 y 50 + β 9 y 10 D 0 y 10 y 50 + β 10 y 10 y 40 y 60 + β 11 y 10 2 y 50 + β 12 y 10 y 60 2 + β 13 y 10 y 20 y 60 + β 14 y 10 D 0 y 20 y 60 + β 15 y 10 y 50 2 + β 16 y 20 2 y 50 + Ω 2 f ˜ cos ( Ω T 0 ) + η ˜ 1 y 30 + η 2 y 50
( D 0 2 + ω 2 ) y 21 = 2 D 0 D 1 y 20 2 μ ˜ D 0 y 20 + α ˜ 1 y 20 3 + α ˜ 2 y 20 y 10 2 + α ˜ 3 y 20 2 ( D 0 y 20 ) + α ˜ 4 y 10 2 ( D 0 y 20 ) + α ˜ 5 y 20 ( D 0 y 10 ) 2 + α ˜ 6 y 20 ( D 0 y 20 ) 2 + α ˜ 7 y 20 y 10 ( D 0 y 10 ) + γ ˜ 1 y 20 2 y 40 + γ ˜ 2 y 20 ( D 0 y 20 ) y 40 + γ ˜ 3 y 20 y 40 2 + γ ˜ 4 y 20 y 10 y 30 + γ ˜ 5 y 20 y 30 2 + γ ˜ 6 y 20 ( D 0 y 10 ) y 30 + γ ˜ 7 y 10 2 y 40 + γ 8 y 20 y 40 y 60 + γ 9 y 20 ( D 0 y 20 ) y 60 + γ 10 y 20 y 30 y 50 + γ 11 y 20 2 y 60 + γ 12 y 20 y 50 2 + γ 13 y 20 y 10 y 50 + γ 14 y 20 ( D 0 y 10 ) y 50 + γ 15 y 20 y 60 2 + γ 16 y 10 2 y 60 + f ˜ Ω 2 sin ( Ω T 0 ) + η ˜ 3 y 40 + η 4 y 60
( D 0 2 + ω 1 2 ) y 31 = 2 D 0 D 1 y 30 2 ε μ ˜ 1 D 0 y 30 + η ˜ 5 y 10
( D 0 2 + ω 2 2 ) y 41 = 2 D 0 D 1 y 40 2 μ ˜ 2 D 0 y 40 + η ˜ 6 y 20
( D 0 + ρ 1 ) y 50 = η ˜ 7 y 10
( D 0 + ρ 2 ) y 60 = η ˜ 8 y 20
The steady-state solutions of Equations (35)–(37), (43) and (44) can be written as follows:
y 10 ( T 0 , T 1 ) = A 1 ( T 1 ) e i ω T 0 + A ¯ 1 ( T 1 ) e i ω T 0
y 20 ( T 0 , T 1 ) = A 2 ( T 1 ) e i ω T 0 + A ¯ 2 ( T 1 ) e i ω T 0
y 30 ( T 0 , T 1 ) = B 1 ( T 1 ) e i ω 1 T 0 + B ¯ 1 ( T 1 ) e i ω 1 T 0
y 40 ( T 0 , T 1 ) = B 2 ( T 1 ) e i ω 2 T 0 + B ¯ 2 ( T 1 ) e i ω 2 T 0
y 50 ( T 0 , T 1 ) = ψ 1 A 1 ( T 1 ) e i ω T 0 + ψ ¯ 1 A ¯ 1 ( T 1 ) e i ω T 0
y 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 ,   A 1 ( T 1 ) ,   A 2 ( T 1 ) ,   B 1 ( T 1 ) , and B 2 ( T 1 ) are unknown functions that will be determined next. A ¯ 1 ( T 1 ) ,   A ¯ 2 ( T 1 ) ,   B ¯ 1 ( T 1 ) , and B ¯ 2 ( T 1 ) are the complex conjugate functions of A 1 ( T 1 ) ,   A 2 ( T 1 ) ,   B 1 ( T 1 ) , and B 2 ( T 1 ) , respectively. ψ 1 = ρ 1 i ω ρ 1 2 + ω 2 η ^ 7 ,   ψ ¯ 1 = ρ 1 + i ω ρ 1 2 + ω 2 η ^ 7 ,   ψ 2 = ρ 2 i ω ρ 2 2 + ω 2 η ^ 8 , and ψ ¯ 2 = ρ 2 + i ω ρ 2 2 + ω 2 η ^ 8 . Substituting Equations (45)–(50) into Equations (39)–(42), yields
( D 0 2 + ω 2 ) y 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 ) y 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 ) y 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 ) y 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 denotes the complex conjugate terms. To get a bounded solution for Equations (51)–(54), the small divisor and secular terms in these equations should vanish. Accordingly, to get the small divisor at the primary resonance cases, let σ ,   σ 1 , and σ 2 to describe the closeness of Ω ,   ω 1 , and ω 2 to the rotor natural frequency ω as follows:
Ω = ω + σ , ω 1 = ω + σ 1 , ω 2 = ω + σ 2
Substituting Equation (55) into Equations (51)–(54), one can obtain the solvability conditions of Equations (51)–(54) as follows:
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 investigate Equations (56)–(59), the unknown functions can be expressed in the polar forms as follows [41,42]:
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 a 3 ( T 1 ) e i θ 3 ( T 1 ) , B 2 ( T 1 ) = 1 2 a 4 ( T 1 ) e i θ 4 ( T 1 ) }
Substituting Equation (60) into Equations (56)–(59), and then separating the real and imaginary parts of the resulting equations, we have
a ˙ 1 = f 1 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = 1 2 ( 2 μ + η 2 η 7 ρ 1 2 + ω 2 ) a 1 + 1 8 ( α 3 2 β 11 η 7 ρ 1 2 + ω 2 2 ρ 1 β 15 η 7 2 ( ρ 1 2 + ω 2 ) 2 ) a 1 3 + 1 8 ( 2 α 4 β 13 η 8 ρ 2 2 + ω 2 ρ 2 β 14 η 8 ρ 2 2 + ω 2 2 β 16 η 7 ρ 1 2 + ω 2 ) a 1 a 2 2 + 1 8 ( α 4 + α 7 + 2 ρ 2 β 12 η 8 2 ( ρ 2 2 + ω 2 ) 2 β 13 η 8 ρ 2 2 + ω 2 + ρ 2 β 14 η 8 ρ 2 2 + ω 2 ) a 1 a 2 2 cos ( 2 φ 1 2 φ 2 ) + 1 8 ( α 2 ω α 5 ω + ( ρ 2 2 ω 2 ) β 12 η 8 2 ω ( ρ 2 2 + ω 2 ) 2 + ρ 2 β 13 η 8 ω ( ρ 2 2 + ω 2 ) ω β 14 η 8 ρ 2 2 + ω 2 ) a 1 a 2 2 sin ( 2 φ 1 2 φ 2 ) ( 1 2 ω η 1 a 3 + 1 8 ω β 1 a 1 2 a 3 + 1 4 ω β 7 a 2 2 a 3 ) sin ( φ 3 ) 1 8 ω β 3 a 1 a 3 2 sin ( 2 φ 3 ) + 1 8 ( β 2 2 β 8 η 7 ρ 1 2 + ω 2 ) a 1 2 a 3 cos ( φ 3 ) 1 8 ( β 10 η 8 ρ 2 2 + ω 2 ) a 1 a 2 a 4 cos ( φ 4 ) + 1 8 ω ( ρ 2 β 10 η 8 ρ 2 2 + ω 2 ) a 1 a 2 a 4 sin ( φ 4 ) + 1 8 ω β 7 a 2 2 a 3 sin ( 2 φ 1 2 φ 2 + φ 3 ) + 1 8 ω ( β 6 ω ω β 10 η 8 ρ 2 2 + ω 2 ) a 1 a 2 a 4 cos ( 2 φ 1 2 φ 2 φ 4 ) + 1 2 ω ( ω + σ ) 2 f sin ( φ 1 ) + 1 8 ω ( β 4 + ρ 2 β 10 η 8 ρ 2 2 + ω 2 ) a 1 a 2 a 4 sin ( 2 φ 1 2 φ 2 φ 4 ) + 1 8 ω β 5 a 1 a 4 2 sin ( 2 φ 1 2 φ 2 2 φ 4 )
a ˙ 2 = f 2 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = 1 2 ( 2 μ + η 4 η 8 ρ 2 2 + ω 2 ) a 2 + 1 8 ( α 3 2 γ 11 η 8 ρ 2 2 + ω 2 2 ρ 2 γ 15 η 8 2 ( ρ 2 2 + ω 2 ) 2 ) a 2 3 + 1 8 ( 2 α 4 γ 13 η 7 ρ 1 2 + ω 2 ω 3 γ 14 η 7 ρ 1 2 + ω 2 2 γ 16 η 8 ρ 2 2 + ω 2 ) a 2 a 1 2 + 1 8 ( α 4 + α 7 2 ρ 1 γ 12 η 7 2 ( ρ 1 2 + ω 2 ) 2 γ 13 η 7 ρ 1 2 + ω 2 + ρ 1 γ 14 η 7 ρ 1 2 + ω 2 ) a 2 a 1 2 cos ( 2 φ 2 2 φ 1 ) + 1 8 ( α 2 ω α 5 ω + ( ρ 1 2 ω 2 ) γ 12 η 7 2 ω ( ρ 1 2 + ω 2 ) 2 + ρ 1 γ 13 η 7 ω ( ρ 1 2 + ω 2 ) + ω γ 14 η 7 ρ 1 2 + ω 2 ) a 2 a 1 2 sin ( 2 φ 2 2 φ 1 ) ( 1 2 ω η 3 a 4 + 1 8 ω γ 1 a 2 2 a 4 + 1 4 ω γ 7 a 1 2 a 4 ) sin ( φ 4 ) 1 8 ω γ 3 a 2 a 4 2 sin ( 2 φ 4 ) + 1 8 ω ( γ 2 ω 2 ω γ 8 η 8 ρ 2 2 + ω 2 ) a 2 2 a 4 cos ( φ 4 ) 1 8 ω ( ω γ 10 η 7 ρ 1 2 + ω 2 ) a 2 a 1 a 3 cos ( φ 3 ) + 1 8 ω ( ρ 1 γ 10 η 7 ρ 1 2 + ω 2 ) a 2 a 1 a 3 sin ( φ 3 ) + 1 8 ω γ 7 a 1 2 a 4 sin ( 2 φ 2 2 φ 1 + φ 4 ) + 1 8 ω ( γ 6 ω ω γ 10 η 7 ρ 1 2 + ω 2 ) a 2 a 1 a 3 cos ( 2 φ 2 2 φ 1 φ 3 ) 1 2 ω ( ω + σ ) 2 f cos ( φ 2 ) + 1 8 ω ( γ 4 + ρ 1 γ 10 η 7 ρ 1 2 + ω 2 ) a 2 a 1 a 3 sin ( 2 φ 2 2 φ 1 φ 3 ) + 1 8 ω γ 5 a 2 a 3 2 sin ( 2 φ 2 2 φ 1 2 φ 3 )
a ˙ 3 = f 3 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = μ 1 a 3 + η 5 2 ( ω + σ 1 ) a 1 sin ( φ 3 )
a ˙ 4 = f 4 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = μ 2 a 4 + η 6 2 ( ω + σ 2 ) a 2 sin ( φ 4 )
φ ˙ 1 = f 5 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = σ + ρ 1 η 2 η 7 2 ω ( ρ 1 2 + ω 2 ) + 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ρ 1 β 11 η 7 ρ 1 2 + ω 2 + β 15 η 7 2 ( ρ 1 2 ω 2 ) ( ρ 1 2 + ω 2 ) 2 ) a 1 2 + 1 8 ω ( ρ 2 β 13 η 8 ρ 2 2 + ω 2 ω 2 β 14 η 8 ρ 2 2 + ω 2 + 2 ρ 1 β 16 η 7 ρ 1 2 + ω 2 + 2 α 2 + 2 α 5 ω 2 ) a 2 2 + 1 4 ω β 3 a 3 2 + 1 4 ω β 5 a 4 2 + 1 8 ω ( α 2 α 5 ω 2 + ( ρ 2 2 ω 2 ) β 12 η 8 2 ( ρ 2 2 + ω 2 ) 2 + ρ 2 β 13 η 8 ρ 2 2 + ω 2 + ω 2 β 14 η 8 ρ 2 2 + ω 2 ) a 2 2 cos ( 2 φ 1 2 φ 2 ) + 1 8 ω ( α 4 ω α 7 ω + 2 ρ 2 ω β 12 η 8 2 ( ρ 2 2 + ω 2 ) 2 + ω β 13 η 8 ρ 2 2 + ω 2 ρ 2 ω β 14 η 8 ρ 2 2 + ω 2 ) a 2 2 sin ( 2 φ 1 2 φ 2 ) + 1 ω a 1 ( 1 2 η 1 a 3 + 3 8 β 1 a 1 2 a 3 + 1 4 β 7 a 2 2 a 3 + 1 4 ρ 1 β 8 η 7 ρ 1 2 + ω 2 a 1 2 a 3 ) cos ( φ 3 ) + 1 8 ω β 3 a 3 2 cos ( 2 φ 3 ) 1 8 β 8 η 7 ρ 1 2 + ω 2 a 1 a 3 sin ( φ 3 ) 1 8 β 2 a 1 a 3 sin ( φ 3 ) + 1 8 ω ( 2 β 4 + ρ 2 β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 cos ( φ 4 ) + 1 8 ( 2 β 6 + β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 sin ( φ 4 ) + 1 8 ω a 1 β 7 a 2 2 a 3 cos ( 2 φ 1 2 φ 2 + φ 3 ) + 1 8 ω ( β 4 + ρ 2 β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 cos ( 2 φ 1 2 φ 2 φ 4 ) + 1 2 ω a 1 ( ω + σ ) 2 f cos ( φ 1 ) + 1 8 ( β 6 + β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 sin ( 2 φ 1 2 φ 2 φ 4 ) + 1 8 ω β 5 a 4 2 cos ( 2 φ 1 2 φ 2 2 φ 4 )
φ ˙ 2 = f 6 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = σ + ρ 2 η 4 η 8 2 ω ( ω 4 2 + ω 2 ) + 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ρ 2 γ 11 η 8 ρ 2 2 + ω 2 + γ 15 η 8 2 ( ρ 2 2 ω 2 ) ( ρ 2 2 + ω 2 ) 2 ) a 2 2 + 1 8 ω ( ρ 1 γ 13 η 7 ρ 1 2 + ω 2 ω 2 γ 14 η 7 ρ 1 2 + ω 2 + 2 ρ 2 γ 16 η 8 ρ 2 2 + ω 2 + 2 α 2 + 2 α 5 ω 2 ) a 1 2 + 1 4 ω γ 3 a 4 2 + 1 4 ω γ 5 a 3 2 + 1 8 ω ( α 2 α 5 ω 2 + ( ρ 1 2 ω 2 ) γ 12 η 7 2 ( ρ 1 2 + ω 2 ) 2 + ρ 1 γ 13 η 7 ρ 1 2 + ω 2 + ω 2 γ 14 η 7 ρ 1 2 + ω 2 ) a 1 2 cos ( 2 φ 2 2 φ 1 ) + 1 8 ω ( α 4 ω α 7 ω + 2 ρ 1 ω γ 12 η 7 2 ( ρ 1 2 + ω 2 ) 2 ω γ 13 η 7 ρ 1 2 + ω 2 + ρ 1 ω γ 14 η 7 ρ 1 2 + ω 2 ) a 1 2 sin ( 2 φ 2 2 φ 1 ) + 1 ω a 2 ( 1 2 η 3 a 4 + 3 8 γ 1 a 2 2 a 4 + 1 4 γ 7 a 1 2 a 4 + 1 4 ρ 2 γ 8 η 8 ρ 2 2 + ω 2 a 2 2 a 4 ) cos ( φ 4 ) + 1 8 ω γ 3 a 4 2 cos ( 2 φ 4 ) 1 8 γ 8 η 8 ρ 2 2 + ω 2 a 2 a 4 sin ( φ 4 ) 1 8 γ 2 a 2 a 4 sin ( φ 4 ) + 1 8 ω ( 2 γ 4 + ρ 1 γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 cos ( φ 3 ) + 1 8 ( 2 γ 6 + γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 sin ( φ 3 ) + 1 8 ω a 2 γ 7 a 1 2 a 4 cos ( 2 φ 2 2 φ 1 + φ 4 ) + 1 8 ω ( γ 4 + ρ 1 γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 cos ( 2 φ 2 2 φ 1 φ 3 ) + 1 2 ω a 2 ( ω + σ ) 2 f sin ( φ 2 ) + 1 8 ( γ 6 + γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 sin ( 2 φ 2 2 φ 1 φ 3 ) + 1 8 ω γ 5 a 3 2 cos ( 2 φ 2 2 φ 1 2 φ 3 )
φ ˙ 3 = f 7 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = σ 1 + 1 2 ( ω + σ 1 ) a 3 η 5 a 1 cos ( φ 3 ) 1 2 ω ρ 1 η 2 η 7 ω 3 2 + ω 2 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ρ 1 β 11 η 7 ρ 1 2 + ω 2 + ( ρ 1 2 ω 2 ) β 15 η 7 2 ( ρ 1 2 + ω 2 ) 2 ) a 1 2 1 8 ω ( 2 α 2 + 2 α 5 ω 2 + ρ 2 β 13 η 8 ρ 2 2 + ω 2 ω 2 β 14 η 8 ρ 2 2 + ω 2 + 2 ρ 1 β 16 η 7 ρ 1 2 + ω 2 ) a 2 2 1 4 ω β 3 a 3 2 1 4 ω β 5 a 4 2 1 8 ω ( α 2 α 5 ω 2 + ( ρ 2 2 ω 2 ) β 12 η 8 2 ( ρ 2 2 + ω 2 ) 2 + ρ 2 β 13 η 8 ρ 2 2 + ω 2 + ω 2 β 14 η 8 ρ 2 2 + ω 2 ) a 2 2 cos ( 2 φ 1 2 φ 2 ) 1 8 ( 2 ρ 2 β 12 η 8 2 ( ρ 2 2 + ω 2 ) 2 + β 13 η 8 ρ 2 2 + ω 2 ρ 2 β 14 η 8 ρ 2 2 + ω 2 + α 4 α 7 ) a 2 2 sin ( 2 φ 1 2 φ 2 ) 1 ω a 1 ( 1 2 η 1 a 3 + 3 8 β 1 a 1 2 a 3 + 1 4 β 7 a 2 2 a 3 + 1 4 ρ 1 β 8 η 7 ρ 1 2 + ω 2 a 1 2 a 3 ) cos ( φ 3 ) 1 8 ω β 3 a 3 2 cos ( 2 φ 3 ) + 1 8 β 8 η 7 ρ 1 2 + ω 2 a 1 a 3 sin ( φ 3 ) 1 8 β 2 a 1 a 3 sin ( φ 3 ) 1 8 ω ( 2 β 4 + ρ 2 β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 cos ( φ 4 ) 1 8 ( 2 β 6 + β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 sin ( φ 4 ) 1 8 ω a 1 β 7 a 2 2 a 3 cos ( 2 φ 1 2 φ 2 + φ 3 ) 1 8 ω ( β 4 + ρ 2 β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 cos ( 2 φ 1 2 φ 2 φ 4 ) 1 2 ω a 1 ( ω + σ ) 2 f cos ( φ 1 ) 1 8 ( β 6 + β 10 η 8 ρ 2 2 + ω 2 ) a 2 a 4 sin ( 2 φ 1 2 φ 2 φ 4 ) 1 8 ω β 5 a 4 2 cos ( 2 φ 1 2 φ 2 2 φ 4 )
φ ˙ 4 = f 8 ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = σ 2 + 1 2 ( ω + σ 2 ) a 4 η 6 a 2 cos ( φ 4 ) 1 2 ω ρ 2 η 4 η 8 ρ 2 2 + ω 2 1 8 ω ( 3 α 1 + α 6 ω 2 + 2 ρ 2 γ 11 η 8 ρ 2 2 + ω 2 + ( ρ 2 2 ω 2 ) γ 15 η 8 2 ( ρ 2 2 + ω 2 ) 2 ) a 2 2 1 8 ω ( 2 α 2 + 2 α 5 ω 2 + ρ 1 γ 13 η 7 ρ 1 2 + ω 2 ω 2 γ 14 η 7 ρ 1 2 + ω 2 + 2 ρ 2 γ 16 η 8 ρ 2 2 + ω 2 ) a 1 2 1 4 ω γ 3 a 4 2 1 4 ω γ 5 a 3 2 1 8 ω ( α 2 α 5 ω 2 + ( ρ 1 2 ω 2 ) γ 12 η 7 2 ( ρ 1 2 + ω 2 ) 2 + ρ 1 γ 13 η 7 ω 3 2 + ω 2 + ω 2 γ 14 η 7 ρ 1 2 + ω 2 ) a 1 2 cos ( 2 φ 2 2 φ 1 ) 1 8 ω ( α 4 ω α 7 ω + 2 ρ 1 ω γ 12 η 7 2 ( ρ 1 2 + ω 2 ) 2 ω γ 13 η 7 ρ 1 2 + ω 2 + ρ 1 ω γ 14 η 7 ρ 1 2 + ω 2 ) a 1 2 sin ( 2 φ 2 2 φ 1 ) 1 ω a 2 ( 1 2 η 3 a 4 + 3 8 γ 1 a 2 2 a 4 + 1 4 γ 7 a 1 2 a 4 + 1 4 ρ 2 γ 8 η 8 ρ 2 2 + ω 2 a 2 2 a 4 ) cos ( φ 4 ) 1 8 ω γ 3 a 4 2 cos ( 2 φ 4 ) + 1 8 γ 8 η 8 ρ 2 2 + ω 2 a 2 a 4 sin ( φ 4 ) + 1 8 γ 2 a 2 a 4 sin ( φ 4 ) 1 8 ω ( 2 γ 4 + ρ 1 γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 cos ( φ 3 ) 1 8 ( 2 γ 6 + γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 sin ( φ 3 ) 1 8 ω a 2 γ 7 a 1 2 a 4 cos ( 2 φ 2 2 φ 1 + φ 4 ) 1 8 ω ( γ 4 + ρ 1 γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 cos ( 2 φ 2 2 φ 1 φ 3 ) 1 2 ω a 2 ( ω + σ ) 2 f sin ( φ 2 ) 1 8 ( γ 6 + γ 10 η 7 ρ 1 2 + ω 2 ) a 1 a 3 sin ( 2 φ 2 2 φ 1 φ 3 ) 1 8 ω γ 5 a 3 2 cos ( 2 φ 2 2 φ 1 2 φ 3 )
where φ 1 = σ τ θ 1 ,   φ 2 = σ τ θ 2 ,   φ 3 = θ 1 θ 3 σ 1 τ , and φ 4 = θ 2 θ 4 σ 2 τ . Substituting Equations (45)–(50) and (60) into Equations (27)–(32), we can obtain the first-order approximated solution of the dynamical system given by Equations (21)–(26) as follows:
y 1 ( τ ) = a 1 ( τ ) cos ( Ω τ φ 1 ( τ ) )
y 2 ( τ ) = a 2 ( τ ) cos ( Ω τ φ 2 ( τ ) )
y 3 ( τ ) = a 3 ( τ ) cos ( Ω τ φ 1 ( τ ) φ 3 ( τ ) )
y 4 ( τ ) = a 4 ( τ ) cos ( Ω τ φ 2 ( τ ) φ 4 ( τ ) )
y 5 ( τ ) = η 7 a 1 ρ 1 2 + ω 2 ( ρ 1 cos ( Ω τ φ 1 ( τ ) ) + ω sin ( Ω τ φ 1 ( τ ) ) )
y 6 ( τ ) = η 8 a 2 ρ 2 2 + ω 2 ( ρ 2 cos ( Ω τ φ 2 ( τ ) ) + ω sin ( Ω τ φ 2 ( τ ) ) )
Based on Equations (69)–(74), a 1 ( τ ) and a 2 ( τ ) denote the oscillation amplitudes of the rotor system in X 1 and X 2 directions, respectively, while a 3 ( τ ) and a 4 ( τ ) are the oscillation amplitudes of the proposed LIPPF-controller. In addition, φ 1 ( τ ) ,   φ 2 ( τ ) ,   φ 3 ( τ ) , and φ 4 ( τ ) represent the phases of the motions. Moreover, the non-linear autonomous differential equations (i.e., Equations (61)–(68)) govern a j ( τ ) and φ j ( τ ) ( j = 1 , 2 , 3 , 4 ) as a function of the different control parameters (i.e., η j ,   j = 1 , , 8 ). Accordingly, setting a 1 = a 2 = a 3 = a 4 = φ 1 = φ 2 = φ 3 = φ 4 = 0 into Equations (61)–(68), yields eight non-linear algebraic equations that govern the vibration amplitudes of the controlled system at the steady-state condition as follows:
f j ( a 1 , a 2 , a 3 , a 4 , φ 1 , φ 2 , φ 3 , φ 4 ) = 0 ; j = 1 , 2 , , 8
By solving the non-linear algebraic system given by Equation (75) utilizing σ or f as a bifurcation parameter at the different values of the control gains (i.e., η 1 ,   η 2 , η 3 ,   η 4 ) and feedback gains (i.e., η 5 ,   η 6 , η 7 ,   η 8 ), one can explore the performance of the introduced controller (i.e., LIPPF) in suppressing the system oscillation amplitudes ( a 1 and a 2 ) in X 1 and X 2 directions as given in Section 3. In addition, the solution stability of Equation (75) was investigated by checking the eigenvalues of the linearized dynamical system given by Equations (61)–(68). Therefore, let the fixed-point solution of Equation (75) be ( a 10 ,   a 20 ,   a 30 ,   a 40 , φ 10 ,   φ 20 ,   φ 30 ,   φ 40 ), and suppose that ( a 11 ,   a 21 ,   a 31 ,   a 41 , φ 11 ,   φ 21 ,   φ 31 ,   φ 41 ) is a small variation regarding that fixed point solution. Accordingly, we can write:
a j = a j 0 + a j 1 , φ j = φ j 0 + φ j 1 , a ˙ j = a ˙ j 1 , φ ˙ j = φ ˙ j 1 ; j = 1 , 2 , , 8
Substituting Equation (76) into Equations (61)–(68) with expanding for the small variations, keeping the linear terms only, one can obtain the following variational differential equation:
a ˙ 11 = f 1 a 11 a 11 + f 1 a 21 a 21 + f 1 a 31 a 31 + f 1 a 41 a 41 + f 1 φ 11 φ 11 + f 1 φ 21 φ 21 + f 1 φ 31 φ 31 + f 1 φ 41 φ 41
a ˙ 21 = f 2 a 11 a 11 + f 2 a 21 a 21 + f 2 a 31 a 31 + f 2 a 41 a 41 + f 2 φ 11 φ 11 + f 2 φ 21 φ 21 + f 2 φ 31 φ 31 + f 2 φ 41 φ 41
a ˙ 31 = f 3 a 11 a 11 + f 3 a 21 a 21 + f 3 a 31 a 31 + f 3 a 41 a 41 + f 3 φ 11 φ 11 + f 3 φ 21 φ 21 + f 3 φ 31 φ 31 + f 3 φ 41 φ 41
a ˙ 41 = f 4 a 11 a 11 + f 4 a 21 a 21 + f 4 a 31 a 31 + f 4 a 41 a 41 + f 4 φ 11 φ 11 + f 4 φ 21 φ 21 + f 4 φ 31 φ 31 + f 4 φ 41 φ 41
φ ˙ 11 = f 5 a 11 a 11 + f 5 a 21 a 21 + f 5 a 31 a 31 + f 5 a 41 a 41 + f 5 φ 11 φ 11 + f 5 φ 21 φ 21 + f 5 φ 31 φ 31 + f 5 φ 41 φ 41
φ ˙ 21 = f 6 a 11 a 11 + f 6 a 21 a 21 + f 6 a 31 a 31 + f 6 a 41 a 41 + f 6 φ 11 φ 11 + f 6 φ 21 φ 21 + f 6 φ 31 φ 31 + f 6 φ 41 φ 41
φ ˙ 31 = f 7 a 11 a 11 + f 7 a 21 a 21 + f 7 a 31 a 31 + f 7 a 41 a 41 + f 7 φ 11 φ 11 + f 7 φ 21 φ 21 + f 7 φ 31 φ 31 + f 7 φ 41 φ 41
φ ˙ 41 = f 8 a 11 a 11 + f 8 a 21 a 21 + f 8 a 31 a 31 + f 8 a 41 a 41 + f 8 φ 11 φ 11 + f 8 φ 21 φ 21 + f 8 φ 31 φ 31 + f 8 φ 41 φ 41
Accordingly, the solution stability of the non-linear dynamical system given by Equations (61)–(68) can be investigated by exploring the eigenvalues of the linear variational differential equations (i.e., Equations (77)–(84)) (see Ref. [48]).

4. Bifurcation Analysis and Control Performance

The performance of the introduced LIPPF-controller in eliminating the rotor vibration amplitudes ( a 1 ,   a 2 ) is discussed within this section via solving Equation (75) numerically at different values of the control gains ( η 1 , η 2 , η 3 , η 4 ) and feedback gains ( η 5 , η 6 , η 7 , η 8 ) utilizing σ or f as the bifurcation parameter [49]. In addition, the solution stability of Equation (75) was explored according to the eigenvalues of the derived variational equations (i.e., Equations (77)–(84)), where the stable solution is plotted as a solid line and the unstable solution is illustrated as a red dotted line. To simulate the steady-state and transient oscillations of the studied controlled system, the parameters’ values were adopted as follows unless otherwise mentioned [29,39]: p = 1.22 , d = 0.005 , α = 22.5 ° , μ 1 = μ 2 = 0.01 , η 1 = η 2 = η 3 = η 3 = η 4 = η 5 = η 6 = 0.2 ,   η 7 = η 8 = ρ 1 = ρ 2 = 1 , Ω = ω + σ , ω 1 = ω + σ 1 ,   ω 2 = ω + σ 2 ,   σ 1 = σ 2 = 0 , and f = 0.015 . Before proceeding further, it is important to remember that the parameter σ is employed to represent the closeness of the angular speed ( Ω ) of the rotor system to its natural frequency ω = 2 p c o s ( α ) + p 3 as given in Equation (55). Accordingly, in the whole article σ is used as a bifurcation parameter to explore the nature of the rotor vibrations at the primary resonance case (i.e., when Ω is approaching ω from the right or the left side).
The rotor oscillation amplitudes ( a 1 and a 2 ) were obtained as shown in Figure 3 via solving Equation (75) using σ as the main bifurcation parameter at various levels of the rotor eccentricity ( f ) before activating the LIPPF-controller (i.e., when η j = 0 ,   j = 1 , 2 , , 8 ). It is clear from the figure that the uncontrolled system can behave as a linear dynamical system at the primary resonance case (i.e., σ 0 ) when the rotor eccentricity is too small (i.e., f = 0.001 ). However, the figure demonstrates that the increase in the rotor eccentricity beyond 0.001 is resulting in high oscillation amplitudes and complex oscillatory motion of the considered system. In addition, the figure illustrates that the uncontrolled system may have two or three stable periodic solutions in the vicinity of σ = 0 depending on the eccentricity magnitude. Moreover, the system may lose its stability when the rotor eccentricity f becomes 0.013 .
Based on the introduced dimensionless variables (i.e., y 1 = x 1 s 0 , and y 2 = x 2 s 0 ) given before Equation (21), one can deduce that y 1 and y 2 represent the relative oscillations of the rotor system in X 1 and X 2 directions with respect to the air-gap size ( s 0 ). In addition, it was illustrated that y 1 ( τ ) and y 2 ( τ ) are periodic functions as given in Equations (69) and (70). Accordingly, the actual instantaneous vibrational motion of the rotor in X 1 and X 2 directions can be expressed such that x 1 ( t ) = s 0 a 1 ( t ) cos ( Ω t φ 1 ( t ) ) and x 2 ( t ) = s 0 a 2 ( t ) cos ( Ω t φ 1 ( t ) ) , where s 0 a 1 and s 0 a 2 are the steady-state vibration amplitudes in X 1 and X 2 directions, respectively. Moreover, these oscillation amplitudes should be less than the air-gap size (i.e., s 0 a 1 < s 0 and s 0 a 2 < s 0 ) to avoid the catastrophic rub-impact force between the rotor and the magnetic pole housing. This means that the system can oscillate without rub-impact occurrence as long as a 1 < 1 and a 2 < 1 . Based on the above discussion, it is clear from Figure 3 that the uncontrolled system may be subjected to rub-impact force between the rotor and pole housing if the eccentricity magnitude became greater than 0.013 . Therefore, the main goal of this work is to mitigate the system lateral vibration close to zero even if the considered system has a large rotor eccentricity.
Based on the introduced control strategy, Equations (21) and (22) represent the equations of motion of the controlled rotor system, while Equations (23)–(26) are the equations of motion of the proposed controller (i.e., LIPPF-controller), where the LIPPF-controller is a combination of both the conventional PPF-controller (i.e., Equations (23) and (24)) and the linear integral resonant controller (i.e., Equations (25) and (26)). Therefore, we can explore the performance of the PPF-controller only by mitigating the non-linear oscillation of the rotor system when setting η 2 = η 4 = η 7 = η 8 = 0 . Figure 4 shows the oscillation amplitudes of both the rotor system (i.e., a 1 and a 2 ) and the connected PPF-controller (i.e., a 3 and a 4 ) along the σ -axis at different levels of the rotor eccentricity when η 1 = η 3 = η 5 = η 6 = 0.2 and η 2 = η 4 = η 7 = η 8 = 0 . Comparing Figure 3 and Figure 4, it is clear that the coupling of the conventional PPF-controller to the considered system eliminated the rotor vibrations in the vicinity of σ = 0 . However, the vibration amplitudes of the rotor system were amplified on both the right and left sides of σ = 0 , which is one of the drawbacks of the conventional PPF-controller.
To eliminate this undesired phenomenon of the PPF-controller, the LIPPF-controller was introduced as a new control strategy that possesses the advantages of both PPF and linear integral resonant controllers. Figure 5 shows the oscillation amplitudes of both the rotor system (i.e., a 1 and a 2 ) and the connected LIPPF-controller (i.e., a 3 and a 4 ) when η 1 = η 3 = η 5 = η 6 = 0.2 and η 2 η 7 = η 4 η 8 = 0.2 . It is clear from Figure 5 that the proposed LIPPF-controller eliminated the rotor vibrations in the vicinity of σ = 0 and suppressed the vibration amplitudes to a very small value on both sides of σ = 0 . The vibration suppression mechanism of the introduced LIPPF-controller can be explained simply based on the working principle of both the conventional PPF-controller and the IRC-controller. It is well known that the PPF-controller absorbs the non-linear vibrations of the connected oscillatory system with high efficiency via establishing an energy bridge when the controller’s natural frequency has the same value as the excitation frequency of the targeted system [25,26,27,28,29]. On the other hand, the IRC-controller can mitigate the non-linear oscillation via increasing the linear damping of the targeted system [30,31,32,33,34,35,36,37,38].
Figure 6 shows the steady-state vibration amplitudes of both the rotor system ( a 1 and a 2 ) and the connected LIPPF-controller ( a 3 and a 4 ) along the σ -axis at different values of σ 1 = σ 2 when f = 0.015 . It is clear from the figure that the optimum working condition of the LIPPF-controller occurs along the dashed red line shown in Figure 6a,b (i.e., when σ = σ 1 = σ 2 ), where a 1 = a 2 0 . Therefore, by substituting the condition σ = σ 1 = σ 2 into Equation (55), one can report that the optimum vibration suppression efficiency of the LIPPF-controller occurs when setting its natural frequencies ( ω 1 and ω 2 ) to become the same value as the rotor angular speed ( Ω ).
Based on Figure 6, the detuning parameters σ ,   σ 1 , and σ 2 were utilized as one bifurcation parameter via letting σ = σ 1 = σ 2 as shown in Figure 7, where the oscillation amplitudes of both the rotor system and the LIPPF-controller were plotted against the bifurcation parameter σ = σ 1 = σ 2 at different values of the rotor eccentricity f . It is clear from the figure that the rotor oscillation amplitudes ( a 1 and a 2 ) were suppressed to very small values along the σ = σ 1 = σ 2 axis regardless of the magnitude of the rotor eccentricity f , where the excessive vibration energy of the rotor system was channeled to the LIPPF-controller via the energy bridge that was established between the rotor and the controller when tuning ω 1 and ω 2 to be the same values of Ω . Accordingly, it is possible to eliminate the lateral vibrations of the considered system regardless of the rotor angular speed Ω and its eccentricity magnitude via making the introduced LIPPF-controller an adaptive controller. The Adaptive Linear Integral Positive Position Feedback controller (ALIPPF-controller) can be implemented from the engineering point of view via measuring the rotor angular speed Ω to be fed into the digital computer to update the controller’s natural frequency ω 1 and ω 2 to the same value of the measured Ω , instantaneously.

5. Numerical Validations and Temporal Oscillations

The bifurcation diagrams given in Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 were obtained based on the approximate analytical solution (i.e., Equations (69)–(74)) of the controlled system equations of motion (i.e., Equations (21)–(26)), where the steady-state oscillation amplitudes (i.e., a 1 ,   a 2 ,   a 3 , and a 4 ) of Equations (69)–(74) are governed by the non-linear algebraic system given by Equation (75). Accordingly, this section is dedicated to demonstrating the accuracy of that approximate solution via solving the Equations of motion (i.e., Equations (21)–(26)) numerically using the ODE45 MATLAB solver. The numerical simulation of the steady-state vibration amplitudes (i.e., a 1 ,   a 2 ,   a 3 , and a 4 ) was illustrated as small circles (when sweeping the bifurcation parameter ( σ or f ) forward) and as big dots (when sweeping the bifurcation parameter ( σ or f ) backward) as shown in Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12. Figure 8 shows the steady-state lateral oscillations of the uncontrolled rotor system which plotted against the parameter σ when f = 0.013 . In addition, the numerical validation is illustrated as small circles and big dots, where the figure demonstrates the accurate correspondence between the numerical and approximate solutions.
The applied numerical algorithm used to validate the accuracy of the obtained approximate solution can be explained through Figure 8 as follows:
  • Forward-sweep algorithm
    • Set σ as the bifurcation control parameter with an initial value σ i n i t i a l   = 0.3 , step-size Δ σ = 0.02 , and final value σ f i n a l   = 0.3 .
    • Set y j ( 0 ) = 0 ( j = 1 , 2 , , 6 ) as zero initial conditions of the first simulation step.
    • Set σ = σ i n i t i a l .
    • Set   f = 0.013 ,   p = 1.22 ,   d = 0.005 ,   α = 0.3927 ,   μ = 1 2 ( 2 d c o s ( α ) + d ) ,   ω = 2 p c o s ( α ) + p 3 )   , Ω = ω + σ ,   η 1 = η 2 = η 3 = η 4 = η 5 = η 6 = η 7 = η 8 = 0 ,   ρ 1 = ρ 2 = 1 ,   μ 1 = μ 2 = 0.01 ,   σ 1 = σ 2 = 0 ,   ω 1 = ω + σ 1 ,   ω 2 = ω + σ 2 .
    • Compute the system parameters α j , β k ,   γ k ( j = 1 ,   2 , , 8 ,   k = 1 ,   2 , , 16 ) as given in the Appendix A according to the values of the parameters given in step 4.
    • Solve the system equations of motion (i.e., Equations (21)–(26)) using ODE45 MATLAB solver on the time interval 0 τ 2000 to capture the steady-state motion.
    • Find the maximum oscillation amplitudes of y j ( τ ) ( j = 1 , 2 , , 6 ) on the time interval 1900 τ 2000 .
    • Set a j ( σ ) = max { y j ( 1900 τ 200 ) } , ( j = 1 , 2 , , 6 ).
    • Set the initial conditions for the next simulation step such that y j ( 0 ) = a j ( σ ) , ( j = 1 , 2 , , 6 ).
    • Increase the bifurcation parameter = σ + Δ σ .
    • If σ σ f i n a l go to step (4), else go to step 12.
    • Plot σ versus a j ( σ ) , ( j = 1 , 2 , , 4 ) as small circles.
    • End of forward-sweep algorithm.
  • Backward-sweep algorithm
    • Set σ as the bifurcation control parameter with an initial value σ i n i t i a l   =   0.3 , step-size Δ σ = 0.02 , and final value σ f i n a l   =   0.3 .
    • Set y j ( 0 ) = 0 ( j = 1 , 2 , , 6 ) as zero initial conditions of the first simulation step.
    • Set σ = σ i n i t i a l .
    • Set   f = 0.013 ,   p = 1.22 ,   d = 0.005 ,   α = 0.3927 ,   μ = 1 2 ( 2 d c o s ( α ) + d ) ,   ω = 2 p c o s ( α ) + p 3 )   , Ω = ω + σ ,   η 1 = η 2 = η 3 = η 4 = η 5 = η 6 = η 7 = η 8 = 0 ,   ρ 1 = ρ 2 = 1 ,   μ 1 = μ 2 = 0.01 ,   σ 1 = σ 2 = 0 ,   ω 1 = ω + σ 1 ,   ω 2 = ω + σ 2 .
    • Compute the system parameters α j , β k ,   γ k ( j = 1 ,   2 , , 8 ,   k = 1 ,   2 , , 16 ) as given in the Appendix A according to the values of the parameters given in step 4.
    • Solve the system equations of motion (i.e., Equations (21)–(26)) using ODE45 MATLAB solver on the time interval 0 τ 2000 to capture the steady-state motion.
    • Find the maximum oscillation amplitudes of y j ( τ ) ( j = 1 , 2 , , 6 ) on the time interval 1900 τ 2000 .
    • Set a j ( σ ) = max { y j ( 1900 τ 200 ) } , ( j = 1 , 2 , , 6 ).
    • Set the initial conditions for the next simulation step such that y j ( 0 ) = a j ( σ ) , ( j = 1 , 2 , , 6 ).
    • Increase the bifurcation parameter = σ + Δ σ .
    • If σ σ f i n a l go to step (4), else go to step 12.
    • Plot σ versus a j ( σ ) , ( j = 1 , 2 , , 4 ) as big dots.
    • End of backward-sweep algorithm
According to the above algorithms, the accuracy of the obtained bifurcation diagrams (obtained in the three-dimensional space as in Section 4) was validated numerically within this section as shown in Figure 9, Figure 10, Figure 11 and Figure 12, where Figure 9 shows the steady-state oscillation amplitudes of both the rotor system and the conventional PPF-controller only (i.e., when η 1 = η 3 = η 5 = η 6 = 0.2 ,   η 2 = η 4 = η 7 = η 8 = 0 ) when f = 0.013 . The figure illustrates that the connection of the PPF-controller to the rotor system can eliminate the system vibrations in the vicinity of σ = 0 . However, the analytical and numerical solutions confirm that the considered system may exhibit high oscillation amplitudes on the left and right sides of σ = 0 . On the other hand, Figure 10 shows the oscillatory behavior of the rotor system along the σ -axis at f = 0.013 when the LIPPF-controller or the ALIPPF-controller are coupled to the rotor system, independently. It is clear from the figure that the coupling of the LIPPF-controller (i.e., when η 1 = η 3 = η 5 = η 6 = 0.2 ,   η 2 η 7 = η 4 η 8 = 0.2 ) to the rotor system eliminated the lateral vibrations in the vicinity of σ = 0.0 and suppressed the oscillation amplitudes on both sides of σ = 0 . In addition, the numerical and analytical investigations demonstrated that the coupling of the ALIPPF-controller (i.e., when η 1 = η 3 = η 5 = η 6 = 0.2 ,   η 2 η 7 = η 4 η 8 = 0.2 ,   σ 1 = σ 2 = σ ) to the rotor system eliminated the oscillation amplitude regardless of the angular velocity value.
The performance of the ALIPPF-controller (i.e., η 1 = η 3 = η 5 = η 6 = η 2 η 7 = η 4 η 8 = 0.2 ,   σ 1 = σ 2 = σ = 0.15 ), LIPPF-controller (i.e., η 1 = η 3 = η 5 = η 6 = η 2 η 7 = η 4 η 8 = 0.2 ), and the conventional PPF-controller (i.e., η 1 = η 3 = η 5 = η 6 = 0.2 ,   η 2 = η 4 = η 7 = η 8 = 0 ) in eliminating the non-linear oscillation of the considered system was explored in Figure 11 and Figure 12 utilizing f as the main bifurcation parameter at two different values of the rotor angular speed Ω (i.e., when Ω = ω + σ ,   σ = 0 ,   0.15 ). Figure 11 shows that the three controllers (i.e., PPF, LIPPF, and ALIPPF) can suppress the rotor oscillation amplitudes ( a 1 and a 2 ) close to zero at σ = 0.0 , which agrees with Figure 9 and Figure 10. It is clear from Figure 11 that the three control techniques have the same performance in eliminating the rotor oscillation when σ = 0.0 along the f axis. On the other hand, Figure 12 shows that the conventional PPF-controller failed to suppress the rotor lateral oscillations when σ = 0.15 , where the system may lose its stability at the large eccentricity magnitudes. In addition, the figure demonstrates that LIPPF-controller can mitigate the rotor vibrations to a very small level along the f axis. Moreover, the figure confirms that the ALIPPF can eliminate the rotor system regardless of the eccentricity magnitude.
To investigate the transient and steady-state oscillatory motions of both the rotor system and the connected controllers, the system equations of motion (i.e., Equations (21)–(26)) were solved numerically using the ODE45 MATLAB solver according to Figure 12 (i.e., when σ = 0.15 ) at three different values of the rotor eccentricity f (i.e., f = 0.02 ,   0.05 , and 0.08 ) as shown in Figure 13, Figure 14 and Figure 15. Figure 13 compares the performance of the PPF-controller (i.e., η 1 = η 3 = η 5 = η 6 = 0.2 ,   η 2 = η 4 = η 7 = η 8 = 0 ), LIPPF-controller (i.e., η 1 = η 3 = η 5 = η 6 = η 2 η 7 = η 4 η 8 = 0.2 ), and the ALIPPF-controller ( η 1 = η 3 = η 5 = η 6 = η 2 η 7 = η 4 η 8 = 0.2 ,   σ 1 = σ 2 = σ = 0.15 ) when f = 0.02 and σ = 0.15 . Figure 14 and Figure 15 are repetitions of Figure 13, but when f = 0.05 and 0 .08 , respectively. Comparing Figure 13, Figure 14 and Figure 15, one can conclude that the conventional PPF-controller has failed to suppress the rotor vibrations in the three studied cases, where the system performs high amplitude periodic motion (at f = 0.02 ), quasiperiodic motion (at f = 0.05 ), and chaotic oscillation (at f = 0.08 ). On the other hand, Figure 13, Figure 14 and Figure 15 demonstrate the high vibration suppression efficiency of the LIPPF-controller and ALIPPF-controller.

6. Conclusions

Within this article, an advanced control technique (i.e., ALIPPF-controller) was presented to eliminate the non-linear vibrations of the rotor electro-magnetic suspension system. The introduced controller is a combination of both first-order and second-order filters that were coupled linearly to the electro-magnetic suspension system in order to reduce the undesired oscillations. Based on the suggested control method, the equations of motion that govern the dynamics of both the electro-magnetic suspension system and the proposed controller were derived as six non-linear ordinary differential equations. Then, an analytical approximate solution for the established dynamical model was obtained using the multiple time-scales perturbation procedure. The performance of the ALIPPF-controller in eliminating the lateral vibrations of the considered system was explored via plotting the different bifurcation diagrams. According to the introduced investigations given above, we can conclude the following important remarks:
  • The uncontrolled eight-pole electro-magnetic suspension system may respond as a linear system at the small rotor eccentricity ( f ) (i.e., f < 0.01 ).
  • At the large rotor eccentricity, the electro-magnetic suspension system may lose its stability to perform quasi-periodic or chaotic oscillations thanks to the complex bifurcation behaviors.
  • Integrating the conventional PPF-controller into the electro-magnetic suspension system can eliminate the rotor vibration amplitudes at the perfect resonance conditions (i.e., σ = 0 ). However, the PPF-controller can add more excessive vibratory energy to the rotor system rather than suppress it if the perfect resonance conditions have been lost.
  • The coupling of the LIR-controller to the eight-pole electromagnetic suspension system can eliminate the non-linear bifurcation and force the rotor to respond as a linear dynamical system. However, the controlled system may exhibit high oscillation at the perfect resonance conditions (i.e., σ = 0 ).
  • The integration of the LIPPF-controller to the considered system can eliminate rotor oscillations at the perfect resonance conditions as well as suppress the non-linear vibrations to very small magnitudes if the resonance conditions have been lost.
  • The ALIPPF-controller (i.e., when setting ω 1 = ω 2 = Ω ) can eliminate the undesired vibrations of the electro-magnetic suspension system to zero regardless of the angular speed and eccentricity of the rotating shaft.

Author Contributions

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

Funding

This research was funded by Researchers Supporting Project, Taif University, Taif, Saudi Arabia, under grant number TURSP-2020/160. 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 would like to thank Taif University, where this work was supported by the Taif University Researchers Supporting Project number (TURSP-2020/160), Taif, Saudi Arabia.

Conflicts of Interest

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

Nomenclature

y 1 ,   y ˙ 1 ,   y ¨ 1 Instantaneous displacement, velocity, and acceleration of rotor system in X 1 direction, respectively.
y 2 ,   y ˙ 2 ,   y ¨ 2 Instantaneous displacement, velocity, and acceleration of rotor system in X 2 direction, respectively.
y 3 ,   y ˙ 3 ,   y ¨ 3 Instantaneous displacement, velocity, and acceleration of the second-order filter that coupled to the rotor system in X 1 direction.
y 4 ,   y ˙ 4 ,   y ¨ 4 Instantaneous displacement, velocity, and acceleration of the second-order filter that coupled to the rotor system in X 2 direction.
y 5 ,   y ˙ 5 Instantaneous displacement and velocity of the first-order filter that coupled to rotor system in X 1 direction.
y 6 ,   y ˙ 6 Instantaneous displacement and velocity of the first-order filter that coupled to rotor system in X 2 direction.
μ Linear damping parameter of the rotor system in X 1 and X 2 directions.
μ 1 ,   μ 2 Linear damping parameters of the second-order filters that coupled to the rotor system in X 1 and X 2 directions, respectively.
ω The natural frequency of the rotor system in X 1 and X 2 directions.
ω 1 ,   ω 2 Natural frequencies of the second-order filters that coupled to the rotor system in X 1 and X 2 directions, respectively.
ρ 1 ,   ρ 2 Internal-loop feedback gains of the first-order filters that coupled to the rotor system in X 1 and X 2 directions, respectively.
Ω The angular speed of the rotor system.
f The rotor system eccentricity.
η 1 , η 2 ,   η 3 , η 4 Control signal gains of the ALIPPF-controller.
η 5 , η 6 ,   η 7 , η 8 Feedback signal gains of the ALIPPF-controller.
α j , ,   j = 1 , , 7 Non-linear coupling coefficients of the rotor system.
β j ,   γ j ,   j = 1 , , 16 Non-linear coupling coefficients of the ALIPPF-controller.
a 1 ,   a 2 Steady-state oscillation amplitudes of the rotor system in X 1 and X 2 directions, respectively.
a 3 ,   a 4 Steady-state oscillation amplitudes of the ALIPPF-controller.
σ The difference between the rotor angular speed ( Ω ) and its natural frequency ( ω ) : σ = Ω ω .

Appendix A

μ = 1 2 ( 2 d cos ( α ) + d ) , ω 2 = 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 2 d 2 cos 2 ( α ) ,
α 5 = 2 d 2 2 d 2 cos 2 ( α ) , α 6 = 2 d 2 cos 2 ( α ) + d 2 ,
α 7 = 12 d cos 3 ( α ) 12 d cos ( α ) + 4 p d 4 p d cos 2 ( α ) ,
β 1 = 2 p η 1 cos ( α ) 3 η 1 cos 2 ( α ) , γ 1 = 2 p η 3 cos ( α ) 3 η 3 cos 2 ( α ) ,
β 2 = 2 d η 1 cos ( α ) , γ 2 = 2 d η 3 cos ( α ) ,
β 3 = 1 8 η 1 2 , γ 3 = 1 8 η 3 2 ,
β 4 = 2 p η 3 tan ( α ) sin ( α ) 6 η 3 sin 2 ( α ) , γ 4 = 2 p η 1 tan ( α ) sin ( α ) 6 η 1 sin 2 ( α ) ,
β 5 = 1 8 η 3 2 tan 2 ( α ) , γ 5 = 1 8 η 1 2 tan 2 ( α ) ,
β 6 = 2 d η 3 tan ( α ) sin ( α ) , γ 6 = 2 d η 1 tan ( α ) sin ( α ) ,
β 7 = 3 η 1 sin 2 ( α ) , γ 7 = 3 η 3 sin 2 ( α ) ,
β 8 = 1 4 η 1 η 2 , γ 8 = 1 4 η 3 η 4 ,
β 9 = 2 d η 2 cos ( α ) , γ 9 = 2 d η 4 cos ( α ) ,
β 10 = 1 4 η 3 η 4 tan 2 ( α ) , γ 10 = 1 4 η 1 η 2 tan 2 ( α ) ,
β 11 = 2 p η 2 cos ( α ) 3 η 2 cos 2 ( α ) , γ 11 = 2 p η 4 cos ( α ) 3 η 4 cos 2 ( α ) ,
β 12 = 1 8 η 4 2 tan 2 ( α ) , γ 12 = 1 8 η 2 2 tan 2 ( α ) ,
β 13 = 2 p η 4 tan ( α ) sin ( α ) 6 η 4 sin 2 ( α ) , γ 13 = 2 p η 2 tan ( α ) sin ( α ) 6 η 2 sin 2 ( α ) ,
β 14 = 2 d η 4 tan ( α ) sin ( α ) , γ 14 = 2 d η 2 tan ( α ) sin ( α ) ,
β 15 = 1 8 η 4 2 , γ 15 = 1 8 η 4 2 ,
β 16 = 3 η 4 sin 2 ( α ) , γ 16 = 3 η 4 sin 2 ( α ) ,

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.; Mahrous, E.; Awrejcewicz, J. Non-linear dynamics of the six-pole rotor-AMBs under two different control configurations. Non-Linear Dyn. 2020, 101, 2299–2323. [Google Scholar] [CrossRef]
  3. 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]
  4. 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]
  5. 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]
  6. 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]
  7. 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]
  8. Zhang, W.; Zhan, X.P. Periodic and chaotic motions of a rotor-active magnetic bearing with quadratic and cubic terms and time-varying stiffness. Non-Linear Dyn. 2005, 41, 331–359. [Google Scholar] [CrossRef]
  9. 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]
  10. 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]
  11. 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]
  12. 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]
  13. 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]
  14. 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]
  15. 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]
  16. Wu, R.; Zhang, W.; Yao, M.H. Non-linear vibration 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]
  17. 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]
  18. 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]
  19. 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] [CrossRef]
  20. 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]
  21. 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]
  22. 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]
  23. 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]
  24. 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]
  25. 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]
  26. 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] [Green Version]
  27. Warminski, J.; Bochenski, M.; Jarzyna, W.; Filipek, P.; Augustyinak, M. Active suppression of non-linear composite beam vibrations by selected control algorithms. Commun. Non-Linear Sci. Numer. Simul. 2011, 16, 2237–2248. [Google Scholar] [CrossRef]
  28. El-Ganaini, W.A.; Saeed, N.A.; Eissa, M. Positive position feedback (PPF) controller for suppression of nonlinear system vibration. Nonlinear Dyn. 2013, 72, 517–537. [Google Scholar] [CrossRef]
  29. 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]
  30. 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]
  31. 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]
  32. 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]
  33. 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]
  34. 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. Non-Linear Sci. Numer. Simul. 2015, 22, 149–166. [Google Scholar] [CrossRef]
  35. Omidi, E.; Mahmoodi, S.N. Non-linear vibration suppression of flexible structures using non-linear modified positive position feedback approach. Non-Linear Dyn. 2015, 79, 835–849. [Google Scholar] [CrossRef]
  36. Saeed, N.A.; Awrejcewicz, J.; Alkashif, M.A.; Mohamed, M.S. 2D and 3D Visualization for the Static Bifurcations and Nonlinear Oscillations of a Self-Excited System with Time-Delayed Controller. Symmetry 2022, 14, 621. [Google Scholar] [CrossRef]
  37. 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]
  38. 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]
  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. Non-Linear Dyn. 2005, 40, 61–102. [Google Scholar] [CrossRef]
  41. De la Luz Sosa, J.; Olvera-Trejo, D.; Urbikain, G.; Martinez-Romero, O.; Elías-Zúñiga, A.; Lacalle, L.N.L.d. Uncharted Stable Peninsula for Multivariable Milling Tools by High-Order Homotopy Perturbation Method. Appl. Sci. 2020, 10, 7869. [Google Scholar] [CrossRef]
  42. Puma-Araujo, S.D.; Olvera-Trejo, D.; Martínez-Romero, O.; Urbikain, G.; Elías-Zúñiga, A.; López de Lacalle, L.N. Semi-Active Magnetorheological Damper Device for Chatter Mitigation during Milling of Thin-Floor Components. Appl. Sci. 2020, 10, 5313. [Google Scholar] [CrossRef]
  43. Urbikain, G.; Olvera, D.; López de Lacalle, L.N.; Beranoagirre, A.; Elías-Zuñiga, A. Prediction Methods and Experimental Techniques for Chatter Avoidance in Turning Systems: A Review. Appl. Sci. 2019, 9, 4718. [Google Scholar] [CrossRef] [Green Version]
  44. Urbikain, G.; Olvera-Trejo, D.; López de Lacalle, L.N.; Elías-Zúñiga, A. Spindle speed variation technique in turning operations: Modeling and real implementation. J. Sound Vib. 2016, 383, 384–396. [Google Scholar] [CrossRef]
  45. Urbikain, G.; Olvera-Trejo, D.; López de Lacalle, L.N.; Elías-Zúñiga, A. Stability and vibrational behaviour in turning processes with low rotational speeds. Int. J. Adv. Manuf. Technol. 2015, 80, 871–885. [Google Scholar] [CrossRef]
  46. 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] [CrossRef]
  47. Schweitzer, G.; Maslen, E.H. Magnetic Bearings: Theory, Design, and Application to Rotating Machinery; Springer: Berlin/Heidelberg, Germany, 2009. [Google Scholar] [CrossRef]
  48. Slotine, J.-J.E.; Li, W. Applied Non-Linear Control; Prentice Hall: Englewood Cliffs, NJ, USA, 1991. [Google Scholar]
  49. 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) The 3-D visualization of the eight-pole electro-magnetic suspension system, (b) schematic diagram of the eight-pole electro-magnetic suspension system when the rotor lies within its nominal position, (c) schematic diagram of the eight-pole electro-magnetic suspension system when the rotor makes displacements x 1 ( t ) and x 2 ( t ) in X 1 and X 2 directions, respectively.
Figure 1. (a) The 3-D visualization of the eight-pole electro-magnetic suspension system, (b) schematic diagram of the eight-pole electro-magnetic suspension system when the rotor lies within its nominal position, (c) schematic diagram of the eight-pole electro-magnetic suspension system when the rotor makes displacements x 1 ( t ) and x 2 ( t ) in X 1 and X 2 directions, respectively.
Applsci 12 03902 g001
Figure 2. The schematic diagram shows the interconnection between the proposed LIPPF-controller and the electro-magnetic suspension system.
Figure 2. The schematic diagram shows the interconnection between the proposed LIPPF-controller and the electro-magnetic suspension system.
Applsci 12 03902 g002
Figure 3. Vibration amplitudes of the rotor system without the LIPPF-controller using σ as a bifurcation parameter when f = 0.001 ,   0.004 ,   0.007 ,   0.01 , and 0.013 : (a) vibration amplitude in X 1   direction and (b) vibration amplitude in X 2 direction.
Figure 3. Vibration amplitudes of the rotor system without the LIPPF-controller using σ as a bifurcation parameter when f = 0.001 ,   0.004 ,   0.007 ,   0.01 , and 0.013 : (a) vibration amplitude in X 1   direction and (b) vibration amplitude in X 2 direction.
Applsci 12 03902 g003
Figure 4. Vibration amplitudes of the rotor system and PPF-controller using σ as a bifurcation parameter when f = 0.005 ,   0.01 ,   0.015 ,   0.02 , and 0.025 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller.
Figure 4. Vibration amplitudes of the rotor system and PPF-controller using σ as a bifurcation parameter when f = 0.005 ,   0.01 ,   0.015 ,   0.02 , and 0.025 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller.
Applsci 12 03902 g004
Figure 5. Vibration amplitudes of the rotor system and LIPPF-controller using σ as a bifurcation parameter when f = 0.005 ,   0.01 ,   0.015 ,   0.02 , and 0.025 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller.
Figure 5. Vibration amplitudes of the rotor system and LIPPF-controller using σ as a bifurcation parameter when f = 0.005 ,   0.01 ,   0.015 ,   0.02 , and 0.025 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller.
Applsci 12 03902 g005
Figure 6. Vibration amplitudes of the rotor system and LIPPF-controller using σ as a bifurcation parameter at various values of σ 1 = σ 2 when f = 0.015 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller.
Figure 6. Vibration amplitudes of the rotor system and LIPPF-controller using σ as a bifurcation parameter at various values of σ 1 = σ 2 when f = 0.015 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller.
Applsci 12 03902 g006
Figure 7. Vibration amplitudes of the rotor system and ALIPPF-controller using σ = σ 1 = σ 2 as a bifurcation parameter when f = 0.005 ,   0.01 ,   0.015 ,   0.02 , and 0.025 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller.
Figure 7. Vibration amplitudes of the rotor system and ALIPPF-controller using σ = σ 1 = σ 2 as a bifurcation parameter when f = 0.005 ,   0.01 ,   0.015 ,   0.02 , and 0.025 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller.
Applsci 12 03902 g007
Figure 8. Vibration amplitudes of the rotor system without the LIPPF-controller using σ as a bifurcation parameter when f = 0.013 : (a) vibration amplitude in X 1 direction and (b) vibration amplitude in X 2 direction.
Figure 8. Vibration amplitudes of the rotor system without the LIPPF-controller using σ as a bifurcation parameter when f = 0.013 : (a) vibration amplitude in X 1 direction and (b) vibration amplitude in X 2 direction.
Applsci 12 03902 g008
Figure 9. Vibration amplitudes of the rotor system and PPF-controller using σ as a bifurcation parameter when f = 0.013 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller.
Figure 9. Vibration amplitudes of the rotor system and PPF-controller using σ as a bifurcation parameter when f = 0.013 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller.
Applsci 12 03902 g009
Figure 10. Vibration amplitudes of the rotor system, LIPPF-controller, and ALIPPF-controller when f = 0.013 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller and ALIPPF-controller.
Figure 10. Vibration amplitudes of the rotor system, LIPPF-controller, and ALIPPF-controller when f = 0.013 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the LIPPF-controller and ALIPPF-controller.
Applsci 12 03902 g010
Figure 11. Vibration amplitudes of the rotor system, PPF-controller, LIPPF-controller, and ALIPPF-controller using the rotor eccentricity f a bifurcation parameter when σ = 0 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller, LIPPF-controller, and ALIPPF-controller.
Figure 11. Vibration amplitudes of the rotor system, PPF-controller, LIPPF-controller, and ALIPPF-controller using the rotor eccentricity f a bifurcation parameter when σ = 0 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller, LIPPF-controller, and ALIPPF-controller.
Applsci 12 03902 g011
Figure 12. Vibration amplitudes of the rotor system, PPF-controller, LIPPF-controller, and ALIPPF-controller using the rotor eccentricity f a bifurcation parameter when σ = 0.15 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller, LIPPF-controller, and ALIPPF-controller.
Figure 12. Vibration amplitudes of the rotor system, PPF-controller, LIPPF-controller, and ALIPPF-controller using the rotor eccentricity f a bifurcation parameter when σ = 0.15 : (a) vibration amplitude in X 1 direction, (b) vibration amplitude in X 2 direction, and (c,d) vibration amplitudes of the PPF-controller, LIPPF-controller, and ALIPPF-controller.
Applsci 12 03902 g012
Figure 13. Rotor electro-magnetic suspension time response, Poincare map, and the frequency spectrum with the PPF-controller, LIPPF-controller, and ALIPPF-controller at σ = 0.15 and f = 0.02 : (a,d) the temporal oscillations of the rotor system in X 1 and X 2 directions, (b,e) the corresponding Poincare map of the rotor system, (c,f) the corresponding frequency spectrum of the rotor, and (gj) the temporal oscillations of the controllers.
Figure 13. Rotor electro-magnetic suspension time response, Poincare map, and the frequency spectrum with the PPF-controller, LIPPF-controller, and ALIPPF-controller at σ = 0.15 and f = 0.02 : (a,d) the temporal oscillations of the rotor system in X 1 and X 2 directions, (b,e) the corresponding Poincare map of the rotor system, (c,f) the corresponding frequency spectrum of the rotor, and (gj) the temporal oscillations of the controllers.
Applsci 12 03902 g013aApplsci 12 03902 g013b
Figure 14. Rotor electro-magnetic suspension time response, Poincare map, and the frequency spectrum with the PPF-controller, LIPPF-controller, and ALIPPF-controller at σ = 0.15 and f = 0.05 : (a,d) the temporal oscillations of the rotor system in X 1 and X 2 directions, (b,e) the corresponding Poincare map of the rotor system, (c,f) the corresponding frequency spectrum of the rotor, and (gj) the temporal oscillations of the controllers.
Figure 14. Rotor electro-magnetic suspension time response, Poincare map, and the frequency spectrum with the PPF-controller, LIPPF-controller, and ALIPPF-controller at σ = 0.15 and f = 0.05 : (a,d) the temporal oscillations of the rotor system in X 1 and X 2 directions, (b,e) the corresponding Poincare map of the rotor system, (c,f) the corresponding frequency spectrum of the rotor, and (gj) the temporal oscillations of the controllers.
Applsci 12 03902 g014
Figure 15. Rotor electro-magnetic suspension time response, Poincare map, and the frequency spectrum with the PPF-controller, LIPPF-controller, and ALIPPF-controller at σ = 0.15 and f = 0.08 : (a,d) the temporal oscillations of the rotor system in X 1 and X 2 directions, (b,e) the corresponding Poincare map of the rotor system, (c,f) the corresponding frequency spectrum of the rotor, and (gj) the temporal oscillations of the controllers.
Figure 15. Rotor electro-magnetic suspension time response, Poincare map, and the frequency spectrum with the PPF-controller, LIPPF-controller, and ALIPPF-controller at σ = 0.15 and f = 0.08 : (a,d) the temporal oscillations of the rotor system in X 1 and X 2 directions, (b,e) the corresponding Poincare map of the rotor system, (c,f) the corresponding frequency spectrum of the rotor, and (gj) the temporal oscillations of the controllers.
Applsci 12 03902 g015
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.; Awrejcewicz, J.; Mousa, A.A.A.; Mohamed, M.S. ALIPPF-Controller to Stabilize the Unstable Motion and Eliminate the Non-Linear Oscillations of the Rotor Electro-Magnetic Suspension System. Appl. Sci. 2022, 12, 3902. https://doi.org/10.3390/app12083902

AMA Style

Saeed NA, Awrejcewicz J, Mousa AAA, Mohamed MS. ALIPPF-Controller to Stabilize the Unstable Motion and Eliminate the Non-Linear Oscillations of the Rotor Electro-Magnetic Suspension System. Applied Sciences. 2022; 12(8):3902. https://doi.org/10.3390/app12083902

Chicago/Turabian Style

Saeed, Nasser A., Jan Awrejcewicz, Abd Allah A. Mousa, and Mohamed S. Mohamed. 2022. "ALIPPF-Controller to Stabilize the Unstable Motion and Eliminate the Non-Linear Oscillations of the Rotor Electro-Magnetic Suspension System" Applied Sciences 12, no. 8: 3902. https://doi.org/10.3390/app12083902

APA Style

Saeed, N. A., Awrejcewicz, J., Mousa, A. A. A., & Mohamed, M. S. (2022). ALIPPF-Controller to Stabilize the Unstable Motion and Eliminate the Non-Linear Oscillations of the Rotor Electro-Magnetic Suspension System. Applied Sciences, 12(8), 3902. https://doi.org/10.3390/app12083902

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