Study on Clearance-Rubbing Dynamic Behavior of 2-DOF Supporting System of Magnetic-Liquid Double Suspension Bearing

Jianhua Zhao 1,2,3,*, Weidong Yan 1,*, Ziqi Wang 1,*, Dianrong Gao 1,* and Guojun Du 2,* 1 Fluid Power Transmission and Control Laboratory, Yanshan University, Qinhuangdao 066004, China 2 College of Civil Engineering and Mechanics, Yanshan University, Qinhuangdao 066004, China 3 Jiangsu Provincial Key Laboratory of Advanced Manufacture and Process for Marine Mechanical Equipment, Zhenjiang 212003, China * Correspondence: zhaojianhua@ysu.edu.cn (J.Z.); yanweidongysu@163.com (W.Y.); wangziqi174518@163.com (Z.W.); gaodr@ysu.edu.cn (D.G.); dugj2002@ysu.edu.cn (G.D.)


Introduction
The active electromagnetic bearing (AMB) has many defects, such as the insufficient electromagnetic attraction caused by the magnetic pole magnetic saturation, the higher temperature rise in the magnetic pole/coil caused by the copper loss and the eddy current loss [1,2] (the faults of AMB are shown as Figure 1), the bearing characteristics of operation stability of AMB can be limited, and then it becomes "technical bottleneck", which restricts the further development and application promotion of AMB [3,4]. The hydrostatic bearing concept is introduced into AMB to form the novel suspension bearing-Magnetic-Liquid Double Suspension Bearing (MLDSB) [5]. It is supported by the electromagnetic suspension and supplemented by the hydrostatic supporting, and it has the advantages of an electromagnetic system and hydrostatic system [6]. The bearing capacity and stiffness can be improved drastically [7].
MLDSB is composed of a bracket, motor, coupling, multi-diameter shaft, journal bearing unit, axial bearing unit, journal loading motor, axial loading motor, etc., as in Figure 2  The journal bearing unit is composed of a rotor, magnetic sleeve, bearing pocket, magnetic pole, inlet hole, shell, coil, etc., as in Figures 3 and 4 [9]. The materials of magnetic pole and magnetic sleeve are cold-rolled non-oriented silicon-steel, and chromium layer is coated in order to solve the problem of matrix corrosion in the lubricants for a long time [10]   The journal bearing unit is composed of a rotor, magnetic sleeve, bearing pocket, magnetic pole, inlet hole, shell, coil, etc., as in Figures 3 and 4 [9]. The materials of magnetic pole and magnetic sleeve are cold-rolled non-oriented silicon-steel, and chromium layer is coated in order to solve the problem of matrix corrosion in the lubricants for a long time [10].  The hydrostatic bearing concept is introduced into AMB to form the novel suspension bearing-Magnetic-Liquid Double Suspension Bearing (MLDSB) [5]. It is supported by the electromagnetic suspension and supplemented by the hydrostatic supporting, and it has the advantages of an electromagnetic system and hydrostatic system [6]. The bearing capacity and stiffness can be improved drastically [7].
MLDSB is composed of a bracket, motor, coupling, multi-diameter shaft, journal bearing unit, axial bearing unit, journal loading motor, axial loading motor, etc., as in Figure 2  The journal bearing unit is composed of a rotor, magnetic sleeve, bearing pocket, magnetic pole, inlet hole, shell, coil, etc., as in Figures 3 and 4 [9]. The materials of magnetic pole and magnetic sleeve are cold-rolled non-oriented silicon-steel, and chromium layer is coated in order to solve the problem of matrix corrosion in the lubricants for a long time [10]   The adapting principle of MLDSB is shown as Figure 5. The magnetic suspension system adopts a PD control mode, and the hydrostatic system adopts a constant-pressure supply system [11].  With the development of a bearing-rotor system towards being large-scale and high-power, the design gap between the rotor and stator becomes smaller and smaller, which makes the liquid film thickness of MLDSB reduce by 10 times compared with the traditional electromagnetic air gap [12]. (Air gap = Liquid film + Chromium coating, Liquid film is (30, 50) μm, Air gap is (100, 300) μm).
Moreover, the clearance-rubbing fault between magnetic pole and magnetic sleeve can be caused easily by the eccentricity, crack, bending, improper assembly of rotor; the coating cracks and falls off, and then the operation stability and reliability of MLDSB can be reduced significantly [13].
At present, many scholars have researched the clearance-rubbing fault of a Bearing-Rotor system in detail, and achieved fruitful results. Braut et al. [14] proposed the method and experiment of using Variational Mode Decomposition (VMD) to diagnose partial rotor-stator rubbing. The results showed that the VMD method can effectively analyze the partial rubbing phenomenon in a constant and slightly variable speed operation. Bularevich et al. [15] analyzed regularities of the possible separation of additives and products of tribochemical reactions under sliding friction. The results showed that nonregular oscillations of the force of friction in the tribosystems under study are interpreted as the result of the rupture-healing of adsorbed antifriction layers. Kostyuk et al. [16] studied the rubbing between the absolute rigid rotor and the absolute rigid stator. The results showed that the system model has a simple analytic solution, which deepened the understanding of rubbing phenomenon of readers. Zhang et al. [17] researched nonlinear dynamics of a 5-DOF active The adapting principle of MLDSB is shown as Figure 5. The magnetic suspension system adopts a PD control mode, and the hydrostatic system adopts a constant-pressure supply system [11]. The adapting principle of MLDSB is shown as Figure 5. The magnetic suspension system adopts a PD control mode, and the hydrostatic system adopts a constant-pressure supply system [11].  With the development of a bearing-rotor system towards being large-scale and high-power, the design gap between the rotor and stator becomes smaller and smaller, which makes the liquid film thickness of MLDSB reduce by 10 times compared with the traditional electromagnetic air gap [12]. (Air gap = Liquid film + Chromium coating, Liquid film is (30, 50) μm, Air gap is (100, 300) μm).
Moreover, the clearance-rubbing fault between magnetic pole and magnetic sleeve can be caused easily by the eccentricity, crack, bending, improper assembly of rotor; the coating cracks and falls off, and then the operation stability and reliability of MLDSB can be reduced significantly [13].
At present, many scholars have researched the clearance-rubbing fault of a Bearing-Rotor system in detail, and achieved fruitful results. Braut et al. [14] proposed the method and experiment of using Variational Mode Decomposition (VMD) to diagnose partial rotor-stator rubbing. The results showed that the VMD method can effectively analyze the partial rubbing phenomenon in a constant and slightly variable speed operation. Bularevich et al. [15] analyzed regularities of the possible separation of additives and products of tribochemical reactions under sliding friction. The results showed that nonregular oscillations of the force of friction in the tribosystems under study are interpreted as the result of the rupture-healing of adsorbed antifriction layers. Kostyuk et al. [16] studied the rubbing between the absolute rigid rotor and the absolute rigid stator. The results showed that the system model has a simple analytic solution, which deepened the understanding of rubbing phenomenon of readers. Zhang et al. [17] researched nonlinear dynamics of a 5-DOF active With the development of a bearing-rotor system towards being large-scale and high-power, the design gap between the rotor and stator becomes smaller and smaller, which makes the liquid film thickness of MLDSB reduce by 10 times compared with the traditional electromagnetic air gap [12]. (Air gap = Liquid film + Chromium coating, Liquid film is (30, 50) µm, Air gap is (100, 300) µm).
Moreover, the clearance-rubbing fault between magnetic pole and magnetic sleeve can be caused easily by the eccentricity, crack, bending, improper assembly of rotor; the coating cracks and falls off, and then the operation stability and reliability of MLDSB can be reduced significantly [13].
At present, many scholars have researched the clearance-rubbing fault of a Bearing-Rotor system in detail, and achieved fruitful results. Braut et al. [14] proposed the method and experiment of using Variational Mode Decomposition (VMD) to diagnose partial rotor-stator rubbing. The results showed that the VMD method can effectively analyze the partial rubbing phenomenon in a constant and slightly variable speed operation. Bularevich et al. [15] analyzed regularities of the possible separation of additives and products of tribochemical reactions under sliding friction. The results showed that nonregular oscillations of the force of friction in the tribosystems under study are interpreted as the result of the rupture-healing of adsorbed antifriction layers. Kostyuk et al. [16] studied the rubbing between the absolute rigid rotor and the absolute rigid stator. The results showed that the system model has a simple analytic solution, which deepened the understanding of rubbing phenomenon of readers. Zhang et al. [17] researched nonlinear dynamics of a 5-DOF active magnetic bearings-rotor system and analyzed the rotor motion stability under different conditions. The results showed that there is the phenomenon of stable periodic motion, critical quasi-periodic motion and instable chaotic motion when the system is under different parameter conditions. Zhu et al. [18] analyzed the support damping and radial clearance of the flexible backup bearing on the transient behavior of the rotor. The results showed that the rolling element backup bearing with flexiblity and more damping can greatly improve the dynamic behavior of a rotor drop on rolling element backup bearing.
In sum, the traditional rubbing dynamics research method takes Bearing Body as the rigid supporting unit to analyze dynamic characteristics of a Bearing-Rotor system. MLDSB bears the severe Clearance-Rubbing compound impact load and becomes the weakest link for Bearing-Rotor system [19].
The Clearance-Rubbing fault is the external representation, and dynamic behavior and law is the internal mechanism of MLDSB. Therefore, the 2-DOF clearance-rubbing dynamics equation was established and transformed into Taylor series form with the nonlinear component. Dimensionless processing was adopted through the use of a Dimensional Unitary Method [20]. The displacement response of a rotor is numerical, simulated under different eccentricity ratios and speeds [21].

Rubbing Force Model
Due to the small Clearance-Rubbing time, the assumptions can be shown as follows [22]: (1) There is local Elastic Collision and Coulomb Friction between magnetic pole and magnetic sleeve; (2) The small gap between two magnetic poles can be ignored; (3) The frictional thermal effect between magnetic pole and magnetic sleeve can be ignored; (4) The winding leakage, marginal magnetic flux, vortex loss, magnetic material saturation, coupling effect between magnetic poles can be ignored; (5) The inertial force and viscous pressure characteristics of the liquid are ignored.
Clearance-Rubbing phenomenon of MLDSB can be shown as Figure 6. o, o 1 , P T , P N , α, w s , e are, respectively, center of stator, center of rotor, tangential rubbing force, radial rubbing force, angle between rubbing point and x axis, rotation angle, radial displacement [23]. magnetic bearings-rotor system and analyzed the rotor motion stability under different conditions. The results showed that there is the phenomenon of stable periodic motion, critical quasi-periodic motion and instable chaotic motion when the system is under different parameter conditions. Zhu et al. [18] analyzed the support damping and radial clearance of the flexible backup bearing on the transient behavior of the rotor. The results showed that the rolling element backup bearing with flexiblity and more damping can greatly improve the dynamic behavior of a rotor drop on rolling element backup bearing. In sum, the traditional rubbing dynamics research method takes Bearing Body as the rigid supporting unit to analyze dynamic characteristics of a Bearing-Rotor system. MLDSB bears the severe Clearance-Rubbing compound impact load and becomes the weakest link for Bearing-Rotor system [19].
The Clearance-Rubbing fault is the external representation, and dynamic behavior and law is the internal mechanism of MLDSB. Therefore, the 2-DOF clearance-rubbing dynamics equation was established and transformed into Taylor series form with the nonlinear component. Dimensionless processing was adopted through the use of a Dimensional Unitary Method [20]. The displacement response of a rotor is numerical, simulated under different eccentricity ratios and speeds [21].

Rubbing Force Model
Due to the small Clearance-Rubbing time, the assumptions can be shown as follows [22]: (1) There is local Elastic Collision and Coulomb Friction between magnetic pole and magnetic sleeve; (2) The small gap between two magnetic poles can be ignored; (3) The frictional thermal effect between magnetic pole and magnetic sleeve can be ignored; (4) The winding leakage, marginal magnetic flux, vortex loss, magnetic material saturation, coupling effect between magnetic poles can be ignored; (5) The inertial force and viscous pressure characteristics of the liquid are ignored.
Clearance-Rubbing phenomenon of MLDSB can be shown as Figure 6. o, o1, PT, PN, α, ws, e are, respectively, center of stator, center of rotor, tangential rubbing force, radial rubbing force, angle between rubbing point and x axis, rotation angle, radial displacement [23].  The mathematical expression of rubbing force P N and P T can be shown as follows [24] where h 0 is liquid film thickness, k c is radial stiffness of stator, f is friction coefficient. Decompose it in the x and y axis to where f x is rubbing force of x axis, f y is rubbing force of y axis. Equation (2) can be transformed into dimensionless form as follows [25] where f X is dimensionless rubbing force of x axis, f Y is dimensionless rubbing force of y axis, X is dimensionless displacement of x axis, Y is dimensionless displacement of y axis.

Dynamic Equation of 2-DOF Supporting System
Due to the same supporting structure of x axis and y axis, the force diagram of a y-direction single DOF supporting system can be shown as Figure 7 [26].
where h0 is liquid film thickness, kc is radial stiffness of stator, f is friction coefficient. Decompose it in the x and y axis to where fx is rubbing force of x axis, fy is rubbing force of y axis. Equation (2) can be transformed into dimensionless form as follows [25]  where fX is dimensionless rubbing force of x axis, fY is dimensionless rubbing force of y axis, X is dimensionless displacement of x axis, Y is dimensionless displacement of y axis.

Dynamic Equation of 2-DOF Supporting System
Due to the same supporting structure of x axis and y axis, the force diagram of a y-direction single DOF supporting system can be shown as Figure 7  The dynamic model of supporting system can be shown as follows [27] where m is mass of rotor, ΔFx and ΔFy are magnetic force of x and y axis, ΔPx and ΔPy are hydrostatic force of x and y axis, ux and uy are the disturbing force of x and y axis.
(1) ΔFx, ΔFy can be transformed into an increment form as follows [28]  The dynamic model of supporting system can be shown as follows [27] m ..
where m is mass of rotor, ∆F x and ∆F y are magnetic force of x and y axis, ∆P x and ∆P y are hydrostatic force of x and y axis, u x and u y are the disturbing force of x and y axis. (1) ∆F x , ∆F y can be transformed into an increment form as follows [28] where i 0 is bias current, k is electromagnetic coefficient, k = µ 0 SN 2 /4, θ is angle between electromagnetic force and y axis, δ 0 is air gap, µ 0 is air permeability, S is area of magnetic pole, N is number of coil, F x0 , F y0 are external load of x, y axis, P x , D x are proportionality and differential coefficient of x axis, P y , D y are proportionality and differential coefficients of y axis; (2) ∆P x , ∆P y can be transformed into increment form as follows [29], where A is hydrostatic supporting coefficient, A = 2p s A e cosθ, p s is system pressure, A e , A b are supporting area and extrusion area, R h0 is initial hydraulic resistor, β is throttle ratio, x = x cos θ/h 0 , y = y cos θ/h 0 ; (3) Disturbing force can be shown as follows [30] u x = me x w 2 s cos(w s t) + me y w 2 s sin(w s t) u y = −me x w 2 s sin(w s t) + me y w 2 s cos(w s t) where, e x , e y are eccentricity of x and y axis, w s is rotation speed of rotor. ∆F x , ∆F y , ∆P x , ∆P y can be transformed to Taylor Series with nonlinear part and substituted into Equation (4), and the dimensionless equation can be shown as follows [31,32] where w n is nature frequency, U x , U y are eccentricity ratio of x and y axis.

Simulation of Clearance-Rubbing of 2-DOF Supporting System
The initial design parameters of MLDSB can be shown as Table 1. U y = 0.1, the influence of design parameters (eccentricity ratio U x and rotation speed n) on rotor displacement and rubbing force can be shown as follows. In Section 3, all the parameters in the figure are dimensionless except for the rotation speed.

Influence of U x on Rotor Displacement
Rotor displacement under different eccentricity ratios, U x , can be simulated with the fourth-order Runge-Kutta method, as in Figure 8 [33].
According to Figure 8a,b, the rotor keeps stable and its displacement is single trajectory when U x < 0.2. The displacement has a positive value under the slow rotation speed. As rotation speed increases, the displacement of rotor is translated from a positive value into a negative value, and it reaches the maximum value when n = 9.5 krpm. The rotor displacement reduces to equilibrium state when rotation speed increases again.
According to Figure 8c,d, the rotor changes in the order of single-period, quasi-period, double-period, chaos and single-period when U x ≥ 0.3 and n ∈ (4.5, 18). As U x increases, the unstable bifurcation interval and double-period interval increase, as in Table 2.

Influence of Rotation Speed on Rotor Displacement
Take U x = 0.3, for example, where the bifurcation rule of rotor can be analyzed as in Figure 8c. The displacement of rotor is the stable single curve when n ∈ (4.5, 9.5). Rotor changes are in the order of quasi-period, double-period, chaos when n ∈ (9.5, 12.5). The displacement of rotor becomes the stable single curve again when n ∈ (12.5, 18), as in Table 3. U x ∈ (0.2, 0.3) is the transitional interval from a stable single-period motion to an unstable motion, and then it can be analyzed in detail in Section 3.3 in order to determine the critical U x . According to Figure 8a,b, the rotor keeps stable and its displacement is single trajectory when Ux < 0.2. The displacement has a positive value under the slow rotation speed. As rotation speed increases, the displacement of rotor is translated from a positive value into a negative value, and it reaches the maximum value when n = 9.5 krpm. The rotor displacement reduces to equilibrium state when rotation speed increases again.   Table 3, the displacement of the rotor shows the rich dynamic characteristics. Therefore, the rule of the rotor can be analyzed at the different rotational speeds, as shown in the Figures below.
According to Figure 9, when n = 5.5 krpm, the displacement of rotor is the regular sinusoidal motion, as shown in Figure 9a, the single discrete peak occurs in the range of one frequency, as shown in Figure 9b, the mass center of rotor is single closed elliptical trajectory, as shown in Figure 9c, the mapping of Poincare is single point, as shown in Figure 9d, and the dimensionless rubbing force is zero without clearance-rubbing phenomenon, as shown in Figure 9e.

Influence of Rotation Speed on Rotor Displacement
Take Ux = 0.3, for example, where the bifurcation rule of rotor can be analyzed as in Figure 8c. The displacement of rotor is the stable single curve when n ∈ (4.5, 9.5). Rotor changes are in the order of quasi-period, double-period, chaos when n ∈ (9.5, 12.5). The displacement of rotor becomes the stable single curve again when n ∈ (12.5, 18), as in Table 3.  Table 3, the displacement of the rotor shows the rich dynamic characteristics. Therefore, the rule of the rotor can be analyzed at the different rotational speeds, as shown in the Figures below. According to Figure 9, when n = 5.5 krpm, the displacement of rotor is the regular sinusoidal motion, as shown in Figure 9a, the single discrete peak occurs in the range of one frequency, as shown in Figure 9b, the mass center of rotor is single closed elliptical trajectory, as shown in Figure  9c, the mapping of Poincare is single point, as shown in Figure 9d, and the dimensionless rubbing force is zero without clearance-rubbing phenomenon, as shown in Figure 9e.  According to Figure 10, when n = 10 krpm, the displacement of rotor is the quasi-period motion (the irregular small-amplitude sinusoidal motion), as shown in Figure 10a, the no common-divisor spectrum occurs in the range of double frequency, as shown in Figure 10b, the mass center of rotor is multi-closed elliptical trajectory, as shown in Figure 10c, the mapping of Poincare is closed loop, as shown in Figure 10d, and the dimensionless rubbing force is 0.8 with clearance-rubbing phenomenon, as shown in Figure 10e.
According to Figure 10, when n = 10 krpm, the displacement of rotor is the quasi-period motion (the irregular small-amplitude sinusoidal motion), as shown in Figure 10a, the no common-divisor spectrum occurs in the range of double frequency, as shown in Figure 10b, the mass center of rotor is multi-closed elliptical trajectory, as shown in Figure 10c, the mapping of Poincare is closed loop, as shown in Figure 10d, and the dimensionless rubbing force is 0.8 with clearance-rubbing phenomenon, as shown in Figure 10e. According to Figure 11, when n = 10.7 krpm, the displacement of rotor is the double-period motion (the irregular sinusoidal motion and the transitional state from quasi-period to chaos state), as shown in Figure 11a, the finite discrete spectrum occurs in the range of double frequency, as shown in Figure 11b, the mass center of rotor is 3-closed elliptical trajectory, as shown in Figure 11c, the mapping of Poincare is 3 points, as shown in Figure 11d, the dimensionless rubbing force is 1.5 with clearance-rubbing phenomenon, and the peak number of rubbing-force decreases, as shown in Figure  11e. According to Figure 11, when n = 10.7 krpm, the displacement of rotor is the double-period motion (the irregular sinusoidal motion and the transitional state from quasi-period to chaos state), as shown in Figure 11a, the finite discrete spectrum occurs in the range of double frequency, as shown in Figure 11b, the mass center of rotor is 3-closed elliptical trajectory, as shown in Figure 11c, the mapping of Poincare is 3 points, as shown in Figure 11d, the dimensionless rubbing force is 1.5 with clearance-rubbing phenomenon, and the peak number of rubbing-force decreases, as shown in Figure 11e.
According to Figure 12, when n = 11.5 krpm, the displacement of rotor is the chaos motion, as shown in Figure 12a, the continuous spectrum occurs in the range of double frequency, as shown in Figure 12b, the mass center of rotor is multi irregular overlapped elliptical trajectory, as shown in Figure 12c, the mapping of Poincare is cloud-shaped closed graph, as shown in Figure 12d, the dimensionless rubbing force is 3.0 with clearance-rubbing phenomenon, and the peak number of rubbing-force decreases to the lowest level, as shown in Figure 12e.
According to Figure 13, when n = 13.5 krpm, the displacement of rotor is the single-period sinusoidal motion, as shown in Figure 13a, the single discrete spectrum occurs in the range of one frequency, as shown in Figure 13b, the mass center of rotor is the regular elliptical trajectory, as shown in Figure 13c, the mapping of Poincare is one point, as shown in Figure 13d, and the dimensionless rubbing force is zero without clearance-rubbing phenomenon, as shown in Figure 13e.

Influence of Rotor Displacement on Eccentricity Ratio and Rotation Speed
Maximum displacement, number of point and maximum rubbing-force are selected as the indexes of the x-t curve, Poincare mapping and clearance-rubbing phenomenon, and then the motion rule of rotor can be analyzed in the range of U x ∈ (0.2, 0.3), as shown in the Figures below.
According to Figure 14, rotor presentsa single-period trajectory and its amplitude is in the range of (−0.6, 0.3) without clearance-rubbing fault within the largest range, as in the white zone. A series of complicated motions of rotor, such as single-period, double-period and chaos, occurs in turn, and its displacement scope is (−1, 1.5) when n ∈ (9, 13) krpm and U x ∈ (0.27, 0.4), as in the color zone. The color depth presents the severity of bifurcating phenomenon as color zone of Figure 14b. According to Figure 11, when n = 10.7 krpm, the displacement of rotor is the double-period motion (the irregular sinusoidal motion and the transitional state from quasi-period to chaos state), as shown in Figure 11a, the finite discrete spectrum occurs in the range of double frequency, as shown in Figure 11b, the mass center of rotor is 3-closed elliptical trajectory, as shown in Figure 11c, the mapping of Poincare is 3 points, as shown in Figure 11d, the dimensionless rubbing force is 1.5 with clearance-rubbing phenomenon, and the peak number of rubbing-force decreases, as shown in Figure  11e According to Figure 12, when n = 11.5 krpm, the displacement of rotor is the chaos motion, as shown in Figure 12a, the continuous spectrum occurs in the range of double frequency, as shown in Figure 12b, the mass center of rotor is multi irregular overlapped elliptical trajectory, as shown in Figure 12c, the mapping of Poincare is cloud-shaped closed graph, as shown in Figure 12d, the dimensionless rubbing force is 3.0 with clearance-rubbing phenomenon, and the peak number of rubbing-force decreases to the lowest level, as shown in Figure 12e.  According to Figure 12, when n = 11.5 krpm, the displacement of rotor is the chaos motion, as shown in Figure 12a, the continuous spectrum occurs in the range of double frequency, as shown in Figure 12b, the mass center of rotor is multi irregular overlapped elliptical trajectory, as shown in Figure 12c, the mapping of Poincare is cloud-shaped closed graph, as shown in Figure 12d, the dimensionless rubbing force is 3.0 with clearance-rubbing phenomenon, and the peak number of rubbing-force decreases to the lowest level, as shown in Figure 12e. According to Figure 13, when n = 13.5 krpm, the displacement of rotor is the single-period sinusoidal motion, as shown in Figure 13a, the single discrete spectrum occurs in the range of one frequency, as shown in Figure 13b, the mass center of rotor is the regular elliptical trajectory, as shown in Figure 13c, the mapping of Poincare is one point, as shown in Figure 13d, and the dimensionless rubbing force is zero without clearance-rubbing phenomenon, as shown in Figure   Figure 12. Rotor displacement when n = 11.5 krpm.

Influence of Rotor Displacement on Eccentricity Ratio and Rotation Speed
Maximum displacement, number of point and maximum rubbing-force are selected as the indexes of the x-t curve, Poincare mapping and clearance-rubbing phenomenon, and then the motion rule of rotor can be analyzed in the range of Ux ∈ (0.2, 0.3), as shown in the Figures below.
According to Figure 14, rotor presentsa single-period trajectory and its amplitude is in the range of (−0.6, 0.3) without clearance-rubbing fault within the largest range, as in the white zone. A series of complicated motions of rotor, such as single-period, double-period and chaos, occurs in turn, and its displacement scope is (−1, 1.5) when n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4), as in the color zone. The color depth presents the severity of bifurcating phenomenon as color zone of Figure  14b.

Influence of Rotor Displacement on Eccentricity Ratio and Rotation Speed
Maximum displacement, number of point and maximum rubbing-force are selected as the indexes of the x-t curve, Poincare mapping and clearance-rubbing phenomenon, and then the motion rule of rotor can be analyzed in the range of Ux ∈ (0.2, 0.3), as shown in the Figures below.
According to Figure 14, rotor presentsa single-period trajectory and its amplitude is in the range of (−0.6, 0.3) without clearance-rubbing fault within the largest range, as in the white zone. A series of complicated motions of rotor, such as single-period, double-period and chaos, occurs in turn, and its displacement scope is (−1, 1.5) when n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4), as in the color zone. The color depth presents the severity of bifurcating phenomenon as color zone of Figure  14b.  According to Figure 15, rotor presents a single-period state without bifurcating phenomenon, and the point number of Poincare mapping is one within the largest range, as in the white zone. A series of points of Poincare mapping occurs and its maximum is up to 15 when n ∈ (9, 13) krpm and U x ∈ (0.27, 0.4), as in the color zone. The color depth presents the point number of Poincare mapping as the color zone of Figure 15b.
According to Figure 15, rotor presents a single-period state without bifurcating phenomenon, and the point number of Poincare mapping is one within the largest range, as in the white zone. A series of points of Poincare mapping occurs and its maximum is up to 15 when n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4), as in the color zone. The color depth presents the point number of Poincare mapping as the color zone of Figure 15b According to Figure 16, there isn't clearance-rubbing phenomenon and its rubbing force is zero within the most range as white zone. The clearance-rubbing phenomenon occurs and its rubbing-force is −5.7 when n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4) as color zone. The color depth presents the magnitude of rubbing force as color zone of Figure 16b. In sum, in order to ensure the stable operation of MLDSB, the rotor should be avoided in the range of n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4) for a long time.

Conclusions
The traditional AMB air gap is 0.1 mm. The bearing will fail when clearance-rubbing fault occurs. Therefore, protective bearings are often added. The hydrostatic bearing concept is introduced into AMB to form the novel suspension bearing MLDSB without protective bearings. In this paper, its operation characteristics and clearance-rubbing fault rules are studied.
When the system rotation speed is close to its natural frequency, it will cause strong resonance, and the eccentricity, crack, bending of the rotor and assembly error will lead to machine faults. Firstly, the clearance-rubbing dimensionless dynamic equation of 2-DOF system of MLDSB is established. Then, rotor displacement under different eccentricity ratios Ux and rotating speeds can be simulated with the fourth-order Runge-Kutta method. Finally, the results show that when the According to Figure 16, there isn't clearance-rubbing phenomenon and its rubbing force is zero within the most range as white zone. The clearance-rubbing phenomenon occurs and its rubbing-force is −5.7 when n ∈ (9, 13) krpm and U x ∈ (0.27, 0.4) as color zone. The color depth presents the magnitude of rubbing force as color zone of Figure 16b.
According to Figure 15, rotor presents a single-period state without bifurcating phenomenon, and the point number of Poincare mapping is one within the largest range, as in the white zone. A series of points of Poincare mapping occurs and its maximum is up to 15 when n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4), as in the color zone. The color depth presents the point number of Poincare mapping as the color zone of Figure 15b According to Figure 16, there isn't clearance-rubbing phenomenon and its rubbing force is zero within the most range as white zone. The clearance-rubbing phenomenon occurs and its rubbing-force is −5.7 when n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4) as color zone. The color depth presents the magnitude of rubbing force as color zone of Figure 16b. In sum, in order to ensure the stable operation of MLDSB, the rotor should be avoided in the range of n ∈ (9, 13) krpm and Ux ∈ (0.27, 0.4) for a long time.

Conclusions
The traditional AMB air gap is 0.1 mm. The bearing will fail when clearance-rubbing fault occurs. Therefore, protective bearings are often added. The hydrostatic bearing concept is introduced into AMB to form the novel suspension bearing MLDSB without protective bearings. In this paper, its operation characteristics and clearance-rubbing fault rules are studied.
When the system rotation speed is close to its natural frequency, it will cause strong resonance, and the eccentricity, crack, bending of the rotor and assembly error will lead to machine faults. Firstly, the clearance-rubbing dimensionless dynamic equation of 2-DOF system of MLDSB is established. Then, rotor displacement under different eccentricity ratios Ux and rotating speeds can be simulated with the fourth-order Runge-Kutta method. Finally, the results show that when the In sum, in order to ensure the stable operation of MLDSB, the rotor should be avoided in the range of n ∈ (9, 13) krpm and U x ∈ (0.27, 0.4) for a long time.

Conclusions
The traditional AMB air gap is 0.1 mm. The bearing will fail when clearance-rubbing fault occurs. Therefore, protective bearings are often added. The hydrostatic bearing concept is introduced into AMB to form the novel suspension bearing MLDSB without protective bearings. In this paper, its operation characteristics and clearance-rubbing fault rules are studied.
When the system rotation speed is close to its natural frequency, it will cause strong resonance, and the eccentricity, crack, bending of the rotor and assembly error will lead to machine faults. Firstly, the clearance-rubbing dimensionless dynamic equation of 2-DOF system of MLDSB is established. Then, rotor displacement under different eccentricity ratios U x and rotating speeds can be simulated with the fourth-order Runge-Kutta method. Finally, the results show that when the eccentricity ratio is greater than 0.27, the rotor displacement exceeds its liquid film limit during the whole acceleration process and clearance-rubbing fault occurs, so the rotation speed range of resonance is (9, 13) krpm. Simultaneously, the value of rotor eccentricity ratio is not required with a working speed lower than