Nonlinear Responses of a Rotor-Bearing-Seal System with Pedestal Looseness

In this paper, a pedestal looseness fault model of a rotor-bearing-seal system is established. Under two working conditions of the same direction eccentricity (Working Condition 1) and reverse eccentricity (Working Condition 2), rotor orbits, vibration waveforms, spectrum cascade, and Poincaré maps are used to study the dynamic characteristics of the system when the sliding bearing support is loosened. (e influence of speed, the unbalance of two discs, the looseness clearance, and the mass of bearing support on system characteristics are analyzed. (e study found that Working Condition 2 is more likely to cause looseness of the bearing support. Moreover, for the rotor system in this paper, the pedestal looseness fault is more likely to occur near the secondorder critical speed. (rough analysis of the spectrum, it is found that the spectrum of the looseness fault will show multiple frequencies or continuous spectra, and the rotor orbits will appear “cylindrical.”


Introduction
In the actual working process of rotating machinery, looseness is a common fault, which usually includes the looseness of rotating parts and the looseness of the foundation. e looseness between bearing support and foundation is the most common fault in rotating machinery. Looseness fault may be caused by long-term severe vibration of the system due to low installation quality and faults such as oil-film instability. Once a looseness fault occurs, it is likely to cause a serious rub impact fault between the rotor and stator. As the rub impact continues to become severe, the loose end may be completely separated from the foundation, resulting in machine crashes and deaths [1][2][3][4]. erefore, the research on the looseness fault of rotating machinery is of great practical significance.
In the research of rotor system looseness fault, Chu and Tang [5] proposed a nonlinear mathematical model that considers the looseness of the rotor support. e periodic solution of the model was obtained by the shooting method, and its stability was analyzed by Floquet theory. Qin et al. [6] used the nonlinear finite element (FE) simulations method to calculate the time-varying stiffness caused by bolt looseness and introduced the calculated time-varying stiffness into the rotor's differential equation of motion, combined with the harmonic balance method to analyze the steady-state response of the system. Research on the nonlinear dynamics of the rotor system with coupled fault involving looseness has also attracted more and more attention [7,8]. Muszynska and Goldman [9] first analyzed the nonlinear dynamics of the rotor system considering the coupled faults of rub impact and looseness. Zhang et al. [10] proposed the wavelet energy entropy method to identify the connection tightness of bolts in the rotor system, and the results showed that the index monotonously increases when the bolt assembly tightness decreases. Jiang et al. [11] studied rotor loosening models considering piecewise linear stiffness, damping force, and nonlinear elastic force; established corresponding mathematical model; and evaluated different degrees of looseness using the nonlinearity-measure-based assessment method. Wei et al. [12] proposed a novel contact model of rotating part looseness that considers bilinear stiffness based on the new characteristics of the rotor rotating part looseness found in experimental tests. Yang et al. [13] used the Lagrangian method to numerically study the coupling faults in the rotor system caused by imbalance, loosening support, and rub impact of rotating and static parts. Yang et al. [14] also studied the coupling faults caused by the loosening of the support and the rub impact based on a dual-rotor test rig, analyzed the vibration of the system at different speeds, established the corresponding mathematical model, and carried out numerical analysis. Yang et al. [15] further studied the rubbing-looseness coupling fault in a rotor system considering geometrical nonlinearity. In addition, many scholars studied the rotor-bearing system with looseness fault and achieved rich research results [16][17][18][19][20].
In recent years, there are an increasing number of scholars who pay attention to the nonlinear vibration of the seal-bearingrotor system because of its strong nonlinear characteristics, including oil-film instability and other coupling faults, such as the airflow-induced vibration [21][22][23][24][25][26][27][28][29][30][31][32]. For example, Wang et al. [21] studied the effect of the real high-pressure methane gas properties on the nonlinear behavior of the seal-rotor system. It is found that when the real gas properties are considered, the axis trajectory of the rotor will increase. Based on the Muszynska nonlinear seal force model, Li et al. [22] studied the stability of the seal-rotor system with rigid support and sliding bearing support, respectively. Ma et al. [23] analyzed the influence of some sealing parameters such as length and diameter on system dynamics. Zhang et al. [24] researched the dynamic characteristics of the rotor-journal bearing system and considered the flow field features of the seal. Ma et al. [26] analyzed the effect of seal force on oil-film instability. Li et al. [27] proposed a new nonlinear model of rotor-bearing-seal system based on the Hamilton principle. It is found that, with the increase of the rotating speed, the system exhibits periodic, double-periodic, multiperiodic, quasiperiodic, and chaotic motion. Under the combined action of air force and oil-film force, the system vibration will become more serious. Luo et al. [29] established the finite element model of seal-ball bearing-rotor test bench and discussed the influence of the rubbing fault inside the seal on system dynamic behavior. In addition, based on the lumped mass method and the principle of displacement equilibrium of the rotor system, Zhou et al. [30] proposed a coupled lateraltorsional dynamic model describing the gear-rotor-seal-bearing (GRSB) system of high-speed centrifugal pumps, and then they analyzed the effects of transmission ratio on the system vibration response.
According to the existing literature, many scholars have studied the nonlinear vibration problem of rotor-bearing pedestal looseness and achieved fruitful results. However, the seal-bearing-rotor system is a strong nonlinear system because of the airflow force of the seal and the oil-film force of the bearing. In particular, both the oil-film instability and airflow-induced vibration are self-excited vibrations caused by fluids. Once they occur, the stable operation of the system will be damaged seriously, the vibration of the rotor system will be aggravated, and secondary faults of the rotor system will be further caused, such as loosening, cracking, and rubbing. However, there are few studies on the dynamics of rotor system with fault under twofluid excitation. Based on this practical background, this paper aims to theoretically reveal the vibration characteristics and influence rules of rotor system with loosening fault under the action of two fluids and to provide a certain basis for fault diagnosis and identification of this kind of rotating machinery.
In this paper, a nonlinear dynamic model of the rotorbearing-seal system considering the sliding bearing, labyrinth seal, and support looseness is first established. e lumped mass method is used to study the system under two working conditions: the in-phase unbalance and out-ofphase unbalance of two discs. e influence of the looseness fault of the right bearing support on the nonlinear dynamic characteristics of the system and the influence of unbalanced mass, loose clearance, and mass of the right bearing support on the system response are determined.

Mathematical Model of a Rotor-Bearing-Seal-Foundation System with Pedestal Looseness
A rotor-bearing-seal-foundation system model with oneend (right end) looseness is established in this paper, as shown in Figure 1. It can be shown in Figure 1 that the left bearing is a selflubricating graphite bearing, which is simulated by springdamping, and the right bearing is a sliding bearing, and the oil-film force adopts the Capone short bearing oil-film force model [33]. e sealing force adopts the Muszynska model [34,35]. According to [26,36], the governing equation of a rotor-bearing-seal system is where M is the mass matrix, G is the gyroscopic matrix, C is the damping matrix (including the left bearing damping), K is the stiffness matrix (including the left bearing stiffness), q is the displacement vector, and F e , F s , F b , and F g are unbalanced force vector, nonlinear sealing force vector, oil-film force vector, and external excitation vector of gravity, respectively.
in which ω is the rotating speed of the rotor. In equations (2) and (3), m i , J pi , and J di (i � 1, 2, 3, 4, 5) are lumped mass, polar moment of inertia of lumped mass points, and diametral moment of inertia of lumped mass points, respectively. i is the number of lumped mass points.
where k blx and k bly represent the stiffnesses of the left bearing in x and y directions, respectively. e matrix elements of K x and K y are as shown as follows: The first loading condition The second loading condition

Shock and Vibration
in which a 1i � 12EI/l 3 i , a 2i � (1/2)l i a 1i , and a 3i � (1/6)l 2 i a 1i . e subscript i represents 1, 2, 3, and 4, respectively. l i is the distance between every two adjacent lumped mass points; E is Young's modulus of elasticity; E � 2.1 × 10 11 Pa. I is the area moment of inertia.
where x i , y i , θ xi , θ yi (i � 1, 2, 3, 4, 5) are displacements in x and y directions and angles of orientation associated with the x and y axes, respectively.
in which C 1 is the rotor damping matrix and C 2 is the bearing damping matrix. For the rotor studied in this paper, the manufacturing material is single, and the properties of the material can be regarded as isotropic. en the damping characteristics of the system are distributed evenly. erefore, the proportional damping is used to simulate the rotor damping C 1 . e Rayleigh damping model, as a kind of proportional damping model, is widely used in engineering practice. e mathematical expression is written as follows: where α and ß are the Rayleigh damping coefficients, and their values can be determined experimentally. According to the basic theory of vibration, equation (8) is transformed into modal coordinates, and the damping matrix C N in canonical coordinate is a diagonal matrix, which can be expressed as [37] In (9), ω j represents the j-th natural frequency. If the first and second natural frequencies ω 1 and ω 2 (rev/min) of the system are known, then α + βω 1 2 � 2ξ 1 ω 1 ，α + βω 2 2 � 2ξ 2 ω 2 . erefore, the first and second modal damping ratios ξ 1 and ξ 2 can be expressed as follows, respectively: e damping coefficients a and ß can be obtained by solving the binary linear equations in (10). e formula is expressed as follows: where c blx and c bly are the dampings of the left bearing in x and y directions, respectively. For the convenience of calculation, the dimensionless transformations are defined as ∆ � ωt and q � q/c.
Equation (1) can be rewritten as a dimensionless form where c represents bearing clearance, and where m 3 e 1 and m 4 e 2 denote eccentricity of unbalance mass of the left and right discs. φ 1 and φ 2 represent initial phase angles of eccentricity in left and right discs.
where F bx and F by represent nonlinear oil-film forces of the right bearing in x and y directions, respectively. e Capone oil-film force model is used in this study. e equations are as follows [33]: where Shock and Vibration 5 where σ � ηω(D/2)L(D/2c) 2 (L/D) 2 is Sommerfeld coefficient. η, L, D, and c denote lubricant viscosity, bearing length, bearing diameter, and bearing radial clearance, respectively.
In (12), the expressions of gravity vector and seal force vector are as follows: where g is the gravity acceleration, which is equal to 9.8 m/s 2 . F f1x , F f1y , F f2x and F f1y are the components of sealing force in x and y directions on nodes 3 and 4, respectively. e Muszynska seal force model [35] is adopted in this study because it is based on a large number of experiments and has a clear physical meaning, which can better reflect the nonlinear characteristics of the seal force. Furthermore, its accuracy has been generally recognized and validated in practical applications. e circumferential angular velocity of fluid near the rotor is ω, and that near the stator is reduced to zero. Further, t is introduced to express the circumferential average velocity ratio of fluid, and τω represents the average flow velocity in the seal cavity. e specific expression of the seal force is as follows: In (18), K f , D f , and m f represent the equivalent stiffness, equivalent damping, and equivalent mass, respectively.
Here, K f , D f , and t are the nonlinear functions related to the disturbance displacements x and y, which can be expressed as follows: In (19), the parameters n f , b 0 , and τ 0 describe the specific seal structure, and e and c represent the relative eccentricity of the rotor and the seal clearance, respectively. Furthermore, K 0 , D 0 , and m f can be expressed as follows [38]: e expressions of parameters presented in (20) are listed as follows: Next, consider the looseness fault in the right bearing. e schematic of the connecting bolt looseness fault is shown in Figure 2. 6 Shock and Vibration At this time, the governing (1) can be expressed as y 6 is the displacement of the bearing support in vertical direction. k f and c f are the equivalent stiffness and equivalent damping between the bearing support and the foundation, respectively, which behave as piecewise linear under different bearing displacement y 6 . e expressions are as follows: where k f1 is the tensile stiffness of the bolt, k f2 is the stiffness when the bolt is loosened, k f3 is the stiffness when the bolt is not loosened, and δ is the clearance between the bearing support and the foundation. e simulation model parameters of the rotor and bearing are shown in Table 1.
e seal parameters related to the rotor-bearing-sealing system are shown in Table 2.
e Newmark-β integral method is adopted to solve (23). e loosening fault occurs at the bearing support, the journal position changes more obviously, so the right end journal and bearing support are taken as research objects. Under two working conditions (Working Conditions 1 and 2 are the same direction eccentricity and reverse eccentricity, respectively), the influences of speed, unbalance mass, looseness clearance, and bearing support mass on the dynamic characteristics of the system are studied. en the dynamic response is observed using the bifurcation diagrams, spectrum cascades, Poincaré maps, FFT spectra, vibration waveform, and rotor orbits. e Poincaré map is adopted to indicate the nature of the system motion. It is a stroboscopic diagram of motion in a phase plane, which is composed of the time series at a constant interval of T r . e rotation period of rotor T r � 2π/ω, where ω is the rotor speed (rad/s).

Dynamic Response under Looseness
According to the dynamic model established in this paper, the dynamic response of the looseness fault rotor system at different speeds under two working conditions is studied. e rotating frequency caused by the rotor unbalance is defined as f r , and its value is one-sixtieth of the speed (r/min). e first-order oil-film oscillation frequency of the sliding bearing is defined as f n1 , and its value is about the first-order natural frequency of the system. e other parameters of the system are the unbalance of the two discs m 3 e 1 � m 4 e 2 � 1.1838 × 10 -4 kg·m, the looseness clearance δ � 1 mm, the tensile stiffness of bolt k f1 � 9 × 10 7 N/m, the stiffness of bolt when the bolt is loosened k f2 � 0 N/m, the   Shock and Vibration 7 stiffness of the bolt when the bolt is not loosened k f3 � 2 × 10 9 N/ m, and the mass of the bearing support m 6 � 1.1 kg.

In-Phase Unbalances of Two Discs.
e spectrum cascade of the right journal (mass 5) in the vertical direction (ydirection) under Working Condition 1 is shown in Figure 3. When the rotating speed is 600 ≤ ω < 6900 r/min, the rotor response is synchronous positive precession. In addition, the frequency f r , 2f r , and 3f r also appear. When the speed reaches 6900 r/min, the first-order oil-film oscillation begins to be violent and reaches the instability speed. At the same time, the lock frequency of the first-order oil-film oscillation (f n1 ) is about 32.2 Hz. Figures 4 and 5 are, respectively, the orbit, the Poincaré maps, and the vibration waveforms of the right journal at various speeds under Working Condition 1.
When ω � 2800 r/min, according to Figure 4, it can be seen that the orbit is an elliptic, and the Poincaré maps are an isolated point. According to Figure 5(a), it can be seen that the vibration waveforms show a sine curve under the unbalanced excitation. ese characteristics all indicate that the system is in period 1 motion. In addition, according to the spectrum cascade in Figure 3, the first-order oil-film oscillation frequency f n1 does not appear. When ω � 5000 r/ min, the orbit, Poincaré maps, and the vibration waveforms all show that the system is still operating stably. When ω � 7500 and 12000 r/min, due to the occurrence of oil-film oscillation, there are many inner "8"s in the orbits. e Poincaré maps are closed loops, indicating a quasiperiodic motion.

Out-of-Phase Unbalances of Two Discs.
In Working Condition 2, the spectrum cascade of the right journal (mass 5) in the vertical direction (y-direction) is shown in Figure 6. When the speed is 6000 ≤ ω ≤ 7200 r/min, it can be clearly seen that, in addition to f r , 2f r , 3f r , and other multipliers, f r /3, f r /2, 2f r /3, 4f r /3, 5f r /3, etc., also appear. is is because in this speed range, it is close to the second-order critical speed, and the vibration is more severe, which is more likely to cause the bearing support to loosen. In addition, when ω ≥ 13800 r/ min, there exists a first-order oil-film oscillation in the system. At this time, the system still has a frequency lock phenomenon at around 32.2 Hz, which is the same as Working Condition 1. e rotor orbits and Poincaré maps of the right journal at various speeds under Working Condition 2 are shown in Figure 7. When the speed ω � 5700 r/min, 6600 r/min, and 6900 r/min, which are close to the second-order critical speed of the system, and the looseness is obvious, so the orbits are "cylindrical," which is a typical orbit characteristic of looseness. At this time, it can be seen from Figure 7(b) that the Poincaré maps show that the system motions are chaotic. Figures 8 and 9 show the vibration waveforms of the right journal and bearing support in the vertical direction (ydirection) under Working Condition 2. It can be seen from Figure 8 that when the speed ω � 5700 r/min and 6900 r/min, the amplitude of the vibration waveforms increases significantly due to the occurrence of looseness. It can be seen from Figure 9 that since the looseness clearance used in the simulation is large (δ �1 mm), the support will not be restricted by the bolt due to the loosening gap, so the vibration waveforms only have single-sided clipping.

The Influence of Loose Rotor System Parameters on the Dynamic
For the rotor system studied in this paper, looseness fault is more likely to be caused under Working Condition 2. For this reason, this section only studies the influence of system parameters on the dynamic for Working Condition 2.

Influence of Unbalance.
Taking the unbalance (disc eccentricity e) as the control parameter, other parameters of the system are as follows: the rotating speed ω � 7100 r/min, the looseness clearance δ � 1 mm, the tensile stiffness of bolt k f1 � 9 × 10 6 N/m, the stiffness of bolt when the bolt is loosened k f2 � 0 N/m, the stiffness of the unloosened bolt k f3 � 2 × 10 9 N/m, and the mass of the bearing support m 6 � 1.1 kg. e spectrum cascade of the right journal (mass 5) and the right bearing support (mass 6) in the vertical direction (y-direction) is as shown in Figure 10. When the eccentricity is 0.1 ≤ e < 0.22 mm, since the support does not loosen, there is no fraction frequency in the vibration response of the journal and the right bearing support at this time. At the same time, it can be seen from Figure 10(b) that with the increase of the eccentricity e, the amplitude of the rotating frequency f r of the right bearing support is increased continuously. is is because the vertical displacement of the support at this time is much smaller than the looseness clearance. When the eccentricity is 0.22 ≤ e ≤ 0.43 mm, a slight looseness has just occurred. In the spectrum cascade, a small fraction frequency appears, as shown in Figure 10. When the eccentricity e > 0.43 mm, because the looseness of the support is obvious, there is an obvious f r /2 fraction frequency, and its amplitude is greater than the amplitude of the frequency f r (see Figure 10(b)). e orbits and Poincaré maps of the right journal with different disc eccentricities under Working Condition 2 are shown in Figure 11. It can be seen that when the eccentricity e is small (e � 0.1 mm), no looseness occurs, and the orbit of the right journal is an ellipse. e Poincaré maps show that the motion is period 1. When the eccentricity e � 0.3 mm, the orbit of the right journal presents an obvious "cylindrical," and the Poincaré maps show chaos. When the eccentricity e � 0.9 mm, the "cylindrical" orbit of the right journal is more obvious, indicating that the increase of the eccentricity e has aggravated the looseness. e vibration waveforms of the right journal and the right bearing support in the vertical direction (y-direction) with different eccentricities under Working Condition 2 are shown in Figures 12 and 13, respectively. It can be seen that, with the increase of the eccentricity e, the amplitude of the vibration waveforms of the journal and the support in the ydirection becomes larger, which also shows that the looseness becomes severe. 8 Shock and Vibration

Influence of Loose Clearance.
Taking the looseness clearance as the control parameter, other parameters of the system are as follows: the rotating speed ω � 6900 r/min, the unbalance m 3 e 1 � m 4 e 2 � 1.1838×10 -4 kg·m, the tensile stiffness of bolt k f1 � 9 × 10 7 N/m, the stiffness of bolt when the bolt is loosened k f2 � 2 × 10 3 N/m, the stiffness of unloosened bolt k f3 � 2 × 10 10 N/m, and the mass of the bearing support m 6 � 1.1 kg. e spectrum cascade of the right journal (mass 5) and the right bearing support (mass 6) in the vertical direction (y-direction) is shown in Figure 14. When the looseness clearance is 1×10 -6 ≤ δ < 0.2674×10 -3 m, because the looseness clearance in this section is small, it is easy to cause the bearing support to contact the top of the bolt, and the displacement of bearing support will be limited by the bolt. erefore, in this looseness clearance interval, the stiffness value k f will change three times, including the loose stiffness k f2 of the bolt, the stiffness k f3 when the bearing support is pressed against the foundation, and the tensile stiffness k f1 of the bolt. As shown in Figures 14(a) and 14(b), there is no fraction frequency, but multiple frequency is the main component. Because the looseness clearance is small at this

Shock and Vibration
time, the upward movement of the support will be restricted by the bolt, so the amplitude of the support response will improve as the looseness clearance increases, as shown in Figure 14(b).
With the increase of the looseness clearance, when δ � 0.2674×10 -3 m, the upward movement of the support sometimes exceeds the loosening clearance and is restricted by the bolt and sometimes does not exceed the looseness clearance. At this time, the number of changes in the stiffness value k f converts from three to two, so in the three-dimensional spectrum, it is continuous, as shown in Figure 14.
When the looseness clearance is 0.3007×10 -3 ≤ δ ≤ 1×10 -3 m, the looseness clearance is large at this time, so the bearing support will not contact the top of the bolt. erefore, the stiffness value k f only changes twice, namely the bolt looseness stiffness k f2 and the stiffness k f3 when the bearing support is compressed with the foundation. It can be seen from Figure 14 that the fraction frequency such as f r /2, 3f r /2, and 5f r /2 appear, and the spectrum of the bearing support shows that the amplitude of the fraction frequency f r /2 is very large, almost equivalent to the amplitude of the rotating frequency f r . Because the motion of the support will not be restricted by the bolt at this time, so the vibration amplitude of the bearing support will not improve when the looseness clearance increases, as shown in Figure 14 are shown in Figure 15. It can be seen that with the increase of the looseness clearance, the orbit changes from the insignificant "cylindrical" when δ � 0.1 × 10 -3 m to the obvious "cylindrical" when δ � 1 × 10 -3 m. e Poincaré maps also show that the motion state of the system gradually changes from quasiperiodic to chaotic, as the looseness clearance increases. e vibration waveforms of the right journal and the right bearing support in the vertical direction (y-direction) with different looseness clearance under Working Condition  It can be seen from Figure 17(b) that when the looseness clearance δ � 0.27×10 -3 m, the upward movement of the support sometimes exceeds the looseness clearance and is restricted by the bolt but sometimes does not exceed the looseness clearance. e number of changes in the stiffness k f is between two and three, so the upper part of the vibration waveforms of the support is not always clipping. It can be seen from Figure 17(c) that when the looseness clearance δ � 1×10 -3 m, due to the large looseness clearance, the displacement of the support will not exceed the looseness clearance and will not be restricted by the tensile force of the bolt, so the vibration waveforms occur clipping only once and only in the lower part.

Influence of the Mass of the Right Bearing Support.
Taking the bearing support m 6 as the control parameter, other parameters of the system are as follows: the rotating speed ω � 6900 r/min, the unbalance m 3 e 1 � m 4 e 2 � 1.1838×10 -4 kg·m, the looseness clearance δ � 1 mm, the tensile stiffness of bolt k f1 � 9 × 10 7 N/m, the stiffness of bolt when the bolt is loosened k f2 � 2 × 10 3 N/m, and the unloosened bolt stiffness k f3 � 2 × 10 10 N/m. e spectrum cascade of the right journal (mass 5) and the right bearing support (mass 6) in the vertical direction (y-direction) is shown in Figure 18. When the mass of the bearing support is 0.5 ≤ m 6 < 0.6 kg, the looseness is more severe due to the small mass of the support, and a large number of continuous spectra appear in the spectrum. As can be seen from Figure 18(b), the amplitude of rotating frequency f r is the largest, which also shows that the looseness is obvious. When the bearing support mass m 6 � 0.6 kg, the amplitude of the fraction frequency f r /2 is the largest. With the increase of the support mass, the amplitude of f r /2 gradually decreases. After m 6 ≥ 1.25 kg, f r /2, 3f r /2, and 5f r /2 components disappeared.
is also shows that the larger bearing pedestal mass can restrain the occurrence of looseness. e orbits and Poincaré maps of the right journal with different right bearing support mass m6 under Working Condition 2 are shown in Figure 19. It can be seen that when the bearing mass m 6 is small (m 6 � 0.55, 1 kg), the "cylindrical" orbits of the right journal are very obvious, and the Poincaré maps show that the system is chaotic. With the increase of the support mass m 6 (m 6 � 1.5 kg), the "cylindrical" orbits of the right journal gradually change from obvious to inconspicuous. When m 6 � 2 kg, the orbits are no longer "cylindrical," and the Poincaré maps show that the system motion is almost in period 1.
e vibration waveforms of the right journal and bearing support in the vertical direction (y-direction) with different right bearing support mass m 6 under Working Condition 2 are shown in Figures 20 and 21, respectively. It can be seen     20 Shock and Vibration that as the right bearing support mass m 6 increases, the amplitude of the vibration waveforms of the journal and the support in the y-direction becomes smaller. is also shows that the effect of looseness is getting less and less obvious, and the larger bearing support mass can restrain the occurrence of looseness fault.

Conclusions
Based on the lumped mass model of the single-span doubledisc rotor-bearing-seal system, the dynamic characteristics of the system when the oil-film bearing support is loosened under the two working conditions of the same direction and reverse eccentricity of the two discs are studied. It mainly analyzes the influence of rotating speed, unbalance, looseness clearance, and bearing support mass on nonlinear vibration characteristics. e main conclusions obtained are as follows: (1) e first-order oil-film oscillation locking frequency of the rotor system under different working conditions is the same, which is the first-order natural frequency of the system. Compared with Working Condition 1, Working Condition 2 is more likely to cause the bearing support to loosen. And for the rotor system studied in this paper, the looseness is more likely to occur near the second critical rotating speed. (2) Looseness faults usually behave as fraction frequencies or continuous spectra in the spectrum. e rotor orbit will appear "cylindrical". Due to the stiffness of the unloosened end being larger, the displacement in the negative vertical direction is limited. From the vibration waveforms, the performance of the limited displacement can be clearly seen. When the displacement of the support exceeds the looseness clearance, the vibration waveforms will be limited by the foundation as well as the bolt. (3) For the rotor system in this paper, when the looseness clearance is small, the frequency components of the system are mainly multiple frequency. When the looseness clearance is large, there are mainly fraction frequencies. e larger mass of the bearing support can restrain the looseness fault to a certain extent, but it will increase the cost. erefore, careful consideration is needed in the design.
Data Availability e datasets supporting the conclusions of this work are included within the article.