# Response Characteristics of Looseness-Rubbing Coupling Fault in Rotor-Sliding Bearing System

Authors: Yang Liu; Zengyuan Xue; Lei Jia; Tuo Shi; Hui Ma
Journal: Mathematical Problems in Engineering
 (2017)

---

## Abstract

For the diagnosis of looseness-rubbing coupling fault of rotor-sliding bearing system caused by bolt looseness fault, the mechanical model and finite element model of dual-disc rotor system with looseness-rubbing coupling fault are established based on the nonlinear finite element method, nonlinear oil film force, looseness stiffness model, and Hertz contact theory. With the augmented Lagrange method, contact constraint conditions are dealt with to ensure that the rotary disk and casing contact with each other meeting boundary penetrating depth within the prescribed tolerance range. And then the dynamics characteristics of the health rotor system supported by sliding bearing are studied. Combined with experimental study and simulation analysis, it is found that the looseness-rubbing coupling fault is often characterized by rubbing fault, the lower part of the time-domain fluctuated shape is denser, while the upper part is relatively loose, and multiple nested half ellipse is shown in orbit diagram. Because of the loosing stiffness and rubbing force, the phenomenon of unstable oil film is depressed. The appearance of the first- and second-order oil film oscillation phenomenon is delayed. It could be used as a theoretical basis for diagnosing looseness-rubbing coupling fault of rotor-sliding bearing system.


---

## Body

## 1. Introduction

Oil film oscillation is a common phenomenon in the rotor-sliding bearing system. The oil film oscillation in the rotating machinery makes the shafting amplitude increase rapidly, which easily initiates the bolt looseness between the bearing chock and the foundation, and furthermore initiates collision and friction. It leads to a periodical rubbing phenomenon, which directly influences the safe operation of unit. Therefore, it has an important significance to research the complicated dynamics characteristic and its evolution regulation of the looseness-rubbing coupling fault of rotor-sliding bearing system to promote the stability and safety of the rotating machinery.In recent years, the research on the looseness-rubbing coupling fault and the system instability problem which is caused by oil film oscillation is always an important subject which the scholars pay attention to. In order to simulate the nonlinearity of the sliding bearing, some nonlinear oil film force models have been proposed, such as Capone G [1, 2]. The Muszynska model [3–5] was introduced as a simple model of nonlinear fluid dynamic forces generated in bearings as well as in seals based on the results of a series of experiments. A parameter called the fluid average circumferential velocity ratio is used to describe the characteristic of the fluid motion as a whole. The model has been used to analyze antiswirl to prevent rotor/seal instability correctly and effectively. de Castro et al. [6] researched the system instability threshold influenced by the unbalance, rotor arrangement form, and bearing parameters. Jing et al. [7] use the oil film force which is obtained from bearing theory to model a continuum rotor-bearing system by the finite element method. Based on Capone model, Adiletta et al. [8] analyzed the possible chaotic motions resulting from the nonlinear response of bearings. Ding and Leung [9] considered an isotropic flexible shaft which is activated by nonlinear fluid-induced forces generated from oil-lubricated journal bearings and hydrodynamic seal.Considering the coupled effect of the nonlinear oil film force and fault force such as looseness or rubbing fault, many researchers studied the dynamic behaviors of the rotor-bearing system. The dynamic behavior in the rubbing of the rotor and stator based on the oil film force of the short bearing is researched by Chu and Zhang [10]. After the analysis of the research, it is found to be an obvious bifurcation and chaotic phenomenon. Liu et al. [11] study the stability of the periodic motion on the looseness-rubbing coupling fault of rotor-sliding bearing system and analyze the change of the bifurcation in different parameter regions of the system. Despite the response feature research of the rotor system in the looseness-rubbing coupling fault, Chen [12] finds that looseness fault can lead the amplitude of the system to increase rapidly and chaotic phenomenon appears. Muszynska and Goldman [13] study the nonlinear dynamic behavior in the rotor system of the base-looseness and rotor-stator rubbing fault and combine the vibration characteristics of the rotor system to complete the relevant experiment. Su et al. [14] use the model of short bearing with unsteady oil film force to study the random response problem of the rotor system where rubbing and support looseness exist. Luo et al. [15] study the periodic motion of the dual-span rotor with the looseness-rubbing coupling fault and point out that the response of the rotor system mainly consists of chaotic motion. Xiang et al. [16] used dynamics theory to establish the mathematical model of an asymmetric double-disc rotor-bearing system which was considering nonlinear oil film force and rub-impact force and pointed out that the phenomenon of the simulation would have a considerable value as diagnostic tools in settling oil film instability. Hu et al. [17] designed a dedicated experimental setup to validate the interaction between rubs and oil film vibration and put forward a model considering nonlinear rub-impact forces and nonlinear oil film forces to analyze the complicated nonlinear vibrations of the rotor-bearing system. The research discloses the worthy energy exchange phenomenon of multifault system and presents the nonlinear dynamical characteristics of the interaction between rub-impact and oil film instability for such a rotor-bearing system.The dynamic finite element model of dual-disc and one-span rotor system with looseness-rubbing coupling fault is supported by oil film force, and based on the looseness stiffness model theory, Hertz contact theory, and the oil film force model of the short bearing, the model is established. And then the research has been conducted to deal with the contact constraint condition through the method of augmented Lagrange to make sure that the permeation depth of the boundary meets prescriptive tolerance range when the disc and the casing contact with each other. Furthermore, by the Newmark-β numerical iteration solution method which is combined with the Newton-Raphson numerical iteration solution method, the excitation response result on the corresponding nodes of the fault rotor system is got. In the meanwhile, the relevant characteristics in the periodic and aperiodic motion of the looseness-rubbing coupling fault of rotor-sliding bearing system are analyzed with experimental study and simulation analysis. The research result represents that the looseness-rubbing coupling fault is often characterized by rubbing fault. And its frequency components are relatively rich. The rotation frequency is as the main frequency for a long time. The lower part of the time-domain fluctuated shape is denser, while the upper part is relatively loose, and multiple nested half ellipse is shown in orbit diagram. This could be used as a theoretical basis for diagnosing looseness-rubbing coupling fault of rotor-sliding bearing system.


## 2. Dynamic Model of Rotor-Bearing System with the Fault

The looseness-rubbing process of the rotor system which is supported by the oil film force is very complex. Because it is the main problem, the following assumption should be defined: (1) the shaft segment of the rotor system is simulated with Timoshenko beam element and only used two nodes (A and B) of end face in the beam element. In the meanwhile, the rotor system neglects the deformation on the axis direction and torsion deformation. Figure1 is the schematic diagram of the shaft segment; (2) the oil film force from the sliding bearing uses the short bearing oil film force model which is in the literature [18]. And the effect of base loosening on oil film force is neglected; (3) only the basic elastic on the looseness direction (y direction) is considered.Figure 1
Timoshenko beam element.Nodes A and B in the Timoshenko beam element each have six DOF.xA,yA,zA,θxA,θyA,θzA are the displacement and angle of the x,y,z direction on the node A; xB,yB,zB,θxB,θyB,θzB are the displacement and angle of the x,y,z direction on the node B. The DOF of the Timoshenko beam element range as u=xA,yA,zA,θxA,θyA,θzA,xB,yB,zB,θxB,θyB,θzBT.As shown in Figure2, the rotor system with looseness-rubbing coupling fault supported by the sliding bearing has 26 nodes in all and divided into 24 shaft segments. To simulate the looseness fault of the support, one node which is named 26 is especially added. It is supported by the oil film force between node 26 and node 25. A segmentation linear spring-damper element is used to simulate between node 26 and the foundation and its equivalent loosening stiffness is shown in (2). The diameters of the two disks are 60 mm. Their location is between shafts 7 and 18. The left bearing and right bearing are, respectively, on node 1 and node 25. The main parameters of the shaft segment are shown in Table 1.Table 1
The main parameters of rotor system shaft segment.
Property
Data
Shaft segment
1
2
3
4
5
6
7
8
9
10
11
12
Length/mm
10
40
40
40
40
40
20
40
40
40
40
40
Diameter/mm
20
20
20
20
20
20
60
20
20
20
20
20
Shaft segment
13
14
15
16
17
18
19
20
21
22
23
24
Length/mm
40
40
40
40
40
20
40
40
40
40
40
10
Diameter/mm
20
20
20
20
20
60
20
20
20
20
20
20Figure 2
Dynamic model of fault rotor system.


### 2.1. The Model of Loosening Stiffness

As shown in Figure3, the tension stiffness of the loosening bolt is kf1. The tension stiffness of the nonloosening bolt is kf2. The stiffness of the foundation (ground) is kf3. When the rotor has a small vibration on the radial direction, the bearing pedestal and the foundation contact with each other. The foundation supports the bearing pedestal. Then the equivalent stiffness of the loosening side is the stiffness of the foundation, kf=kf3. Because of the phenomenon of the bolt looseness, when amplitude is large, the bearing pedestal could get rid of the constraint of the loosening bolt and be separated from the foundation temporarily. At this time, the nonloosening bolt is in the condition of elastic deformation. The equivalent stiffness of the loosening side is kf=kf2; if the displacement of the bearing pedestal is larger than the loosening bolt gap, then the equivalent stiffness kf≈kf1+kf2. So the equivalent stiffness of the loosening [19] side can be summarized as follows:(1)kf=kf1+kf2,y>δ,kf2,0≤y≤δ,kf3,y<0.Figure 3
The schematic diagram of pedestal looseness.An assumption is given that the foundation stiffnesskf3 is approximate with the sum of the tension stiffness of all the loosening bolt. kf3≈kf1+kf2. Equation (1) can be simplified to(2)kf=kf2,0≤y≤δ,kf3,therest.


### 2.2. The Rubbing Force Model Based on Contact Judgment

For the safety, the disc (or vane) of the rotor system should be protected by the casing. It is assumed that the rotary disk 1 is protected by the casing. When the system is at rest, the radius clearance between the rotor and the casing (stator) is shown in Figure4. It is assumed that the cross section of the rotary disk is always keeping in the plane xoy. O is the center of the rotor. o1 is still the center of the stator. ω is the rotational velocity. This paper completes the contact detection based on the Kuhn-Tucker condition and deals with the contact constraint condition through the method of augmented Lagrange [20] to make sure that the permeation depth of the boundary meets prescriptive tolerance range when the rotary disk and the casing contact with each other. So the rubbing force Frub can be divided into the rubbing force on the normal direction Fr and the tangential direction Ft. Its expression is(3)Fr=r-ckr+λri+1,Ft=fFr,r≥c,where kr is the rubbing stiffness; f is the friction coefficient at the friction point; c is the radius gap between the rotor and stator; r is the normal direction of the rotor; λri+1 is the Lagrange multiplier in the iteration process on the i+1 step, and its expression is(4)λri+1=λri+krδr,δr≥ς,λri,δr<ς,where δr=r-c; ς is the allowed permeation tolerance value; when the value is over ς, the rubbing force on the normal direction Fr expands through the Lagrange multiplier λri+1 and iterates continuously until the permeation value is smaller than the maximum allowed value ς.Figure 4
The model of the rubbing force.In thexoy reference frame, the rubbing force along the x direction and y direction can be expressed as(5)FrxFry=-Hδrδrδrkr+λri+1r1-ff1x1y1,where x1 is the displacement of the rotor axis on the x and y direction. y1 is the displacement of the rotor axis on the y direction. Hδr is the step function.(6)Hδr=1,δr≥0,0,δr<0.The meaning of the expression above is that whiler<c, the rubbing phenomenon does not happen and then the rubbing force is 0. However r≥c, the rubbing phenomenon happens and the rubbing force appears.


### 2.3. The Finite Element Model of the Fault Rotor System

As known above, under the support of the sliding bearing, the bolt looseness and the rotary disk rubbing with the blades are considered in the rotor-bearing system of which finite element model is shown below.(7)Mu¨+C+Gu˙+Ku=Fe+Fb+Frub,where M,G,C,K,Fe,Fb,Frub are, respectively, the mass matrix, gyroscope matrix, damping matrix, stiffness matrix, the vector of the unbalance force, the vector of the oil film force, and the vector of the rubbing force. Because the normal deformation on the node of the model is much larger than the deformation on the axial direction and the magnitude has a huge difference, neglect the displacement vector on the axial direction and after the torsion deformation. The displacement vector  u=xA,yA,θxA,θyA,xB,yB,θxB,θyBT. The damping in the model uses the expression of the Rayleigh damping.(8)C=αM+βK,where α is quality proportion coefficient and β is stiffness proportion coefficient(9)α=2ξ2/ωn2-ξ1/ωn11/ωn22-1/ωn12,β=2ξ2ωn2-ξ1ωn1ωn22-ωn12,ωn1 is the first-order critical speed, ωn2 is the second-order critical speed, ξ1 is the first-order modal damping ratio; ξ2 is the second-order modal damping ratio.The two support bearings use the model of the short bearing oil film force and then get a corresponding promotion and normalization processing.Fbx and Fby are, respectively, the nonlinear oil film force (Fb) of the bearing on the x,y direction vector, and its expression is(10)FbxFby=σfbxfby,fbxfby=x~-2y~˙2+y~+2x~˙21/21-x~2-y~2×3x~Vx~,y~,α-sin⁡αGx~,y~,α-2cos⁡αSx~,y~,α3y~Vx~,y~,α+cos⁡αGx~,y~,α-2sin⁡αSx~,y~,α,where(11)Vx~,y~,α=2+y~cos⁡α-x~sin⁡αGx~,y~,α1-x~2-y~2,Sx~,y~,α=x~cos⁡α+y~sin⁡α1-x~cos⁡α+y~sin⁡α2,Gx~,y~,α=21-x~2-y~21/2π2+arctan⁡y~·cos⁡α-x~·sin⁡α1-x~2-y~21/2,α=arctan⁡y~+2x~˙x~-2y~˙-π2sign⁡y~+2x~˙x~-2y~˙-π2sign⁡y~+2x~˙,where σ is the Sommerfeld coefficient; and its expression is(12)σ=ηωD2LD2τ2LD2,where L is bearing width; η is oil viscosity; D is bearing diameter; τ is oil film gap.The Newmark-β numerical algorithm combined with the Newton-Raphson iteration method has been used to solve the nonlinear differential equation (7) which is established formerly. And then the vibrating response data on the relevant nodes of the finite element model are got to analyze the fault of rotor system.


## 2.1. The Model of Loosening Stiffness

As shown in Figure3, the tension stiffness of the loosening bolt is kf1. The tension stiffness of the nonloosening bolt is kf2. The stiffness of the foundation (ground) is kf3. When the rotor has a small vibration on the radial direction, the bearing pedestal and the foundation contact with each other. The foundation supports the bearing pedestal. Then the equivalent stiffness of the loosening side is the stiffness of the foundation, kf=kf3. Because of the phenomenon of the bolt looseness, when amplitude is large, the bearing pedestal could get rid of the constraint of the loosening bolt and be separated from the foundation temporarily. At this time, the nonloosening bolt is in the condition of elastic deformation. The equivalent stiffness of the loosening side is kf=kf2; if the displacement of the bearing pedestal is larger than the loosening bolt gap, then the equivalent stiffness kf≈kf1+kf2. So the equivalent stiffness of the loosening [19] side can be summarized as follows:(1)kf=kf1+kf2,y>δ,kf2,0≤y≤δ,kf3,y<0.Figure 3
The schematic diagram of pedestal looseness.An assumption is given that the foundation stiffnesskf3 is approximate with the sum of the tension stiffness of all the loosening bolt. kf3≈kf1+kf2. Equation (1) can be simplified to(2)kf=kf2,0≤y≤δ,kf3,therest.


## 2.2. The Rubbing Force Model Based on Contact Judgment

For the safety, the disc (or vane) of the rotor system should be protected by the casing. It is assumed that the rotary disk 1 is protected by the casing. When the system is at rest, the radius clearance between the rotor and the casing (stator) is shown in Figure4. It is assumed that the cross section of the rotary disk is always keeping in the plane xoy. O is the center of the rotor. o1 is still the center of the stator. ω is the rotational velocity. This paper completes the contact detection based on the Kuhn-Tucker condition and deals with the contact constraint condition through the method of augmented Lagrange [20] to make sure that the permeation depth of the boundary meets prescriptive tolerance range when the rotary disk and the casing contact with each other. So the rubbing force Frub can be divided into the rubbing force on the normal direction Fr and the tangential direction Ft. Its expression is(3)Fr=r-ckr+λri+1,Ft=fFr,r≥c,where kr is the rubbing stiffness; f is the friction coefficient at the friction point; c is the radius gap between the rotor and stator; r is the normal direction of the rotor; λri+1 is the Lagrange multiplier in the iteration process on the i+1 step, and its expression is(4)λri+1=λri+krδr,δr≥ς,λri,δr<ς,where δr=r-c; ς is the allowed permeation tolerance value; when the value is over ς, the rubbing force on the normal direction Fr expands through the Lagrange multiplier λri+1 and iterates continuously until the permeation value is smaller than the maximum allowed value ς.Figure 4
The model of the rubbing force.In thexoy reference frame, the rubbing force along the x direction and y direction can be expressed as(5)FrxFry=-Hδrδrδrkr+λri+1r1-ff1x1y1,where x1 is the displacement of the rotor axis on the x and y direction. y1 is the displacement of the rotor axis on the y direction. Hδr is the step function.(6)Hδr=1,δr≥0,0,δr<0.The meaning of the expression above is that whiler<c, the rubbing phenomenon does not happen and then the rubbing force is 0. However r≥c, the rubbing phenomenon happens and the rubbing force appears.


## 2.3. The Finite Element Model of the Fault Rotor System

As known above, under the support of the sliding bearing, the bolt looseness and the rotary disk rubbing with the blades are considered in the rotor-bearing system of which finite element model is shown below.(7)Mu¨+C+Gu˙+Ku=Fe+Fb+Frub,where M,G,C,K,Fe,Fb,Frub are, respectively, the mass matrix, gyroscope matrix, damping matrix, stiffness matrix, the vector of the unbalance force, the vector of the oil film force, and the vector of the rubbing force. Because the normal deformation on the node of the model is much larger than the deformation on the axial direction and the magnitude has a huge difference, neglect the displacement vector on the axial direction and after the torsion deformation. The displacement vector  u=xA,yA,θxA,θyA,xB,yB,θxB,θyBT. The damping in the model uses the expression of the Rayleigh damping.(8)C=αM+βK,where α is quality proportion coefficient and β is stiffness proportion coefficient(9)α=2ξ2/ωn2-ξ1/ωn11/ωn22-1/ωn12,β=2ξ2ωn2-ξ1ωn1ωn22-ωn12,ωn1 is the first-order critical speed, ωn2 is the second-order critical speed, ξ1 is the first-order modal damping ratio; ξ2 is the second-order modal damping ratio.The two support bearings use the model of the short bearing oil film force and then get a corresponding promotion and normalization processing.Fbx and Fby are, respectively, the nonlinear oil film force (Fb) of the bearing on the x,y direction vector, and its expression is(10)FbxFby=σfbxfby,fbxfby=x~-2y~˙2+y~+2x~˙21/21-x~2-y~2×3x~Vx~,y~,α-sin⁡αGx~,y~,α-2cos⁡αSx~,y~,α3y~Vx~,y~,α+cos⁡αGx~,y~,α-2sin⁡αSx~,y~,α,where(11)Vx~,y~,α=2+y~cos⁡α-x~sin⁡αGx~,y~,α1-x~2-y~2,Sx~,y~,α=x~cos⁡α+y~sin⁡α1-x~cos⁡α+y~sin⁡α2,Gx~,y~,α=21-x~2-y~21/2π2+arctan⁡y~·cos⁡α-x~·sin⁡α1-x~2-y~21/2,α=arctan⁡y~+2x~˙x~-2y~˙-π2sign⁡y~+2x~˙x~-2y~˙-π2sign⁡y~+2x~˙,where σ is the Sommerfeld coefficient; and its expression is(12)σ=ηωD2LD2τ2LD2,where L is bearing width; η is oil viscosity; D is bearing diameter; τ is oil film gap.The Newmark-β numerical algorithm combined with the Newton-Raphson iteration method has been used to solve the nonlinear differential equation (7) which is established formerly. And then the vibrating response data on the relevant nodes of the finite element model are got to analyze the fault of rotor system.


## 3. Dynamics Characteristic Analysis of Coupling Fault of Rotor System

The main parameters of the looseness-rubbing coupling fault of rotor system are shown below. The shaft diameter is 20 mm, the shaft length is 860 mm, the disc diameter is 60 mm, the disc thickness is 20 mm, the unbalance quality is156e-6 kg·m, the material density ρ=7850 kg/m3, the elastic modulus E=2.06×1011 Pa, the Poisson ratio υ=0.3, the rubbing gap is u=200 μm, the rubbing stiffness kr=2×107 N/m, the friction coefficient f=0.3, and the left and right bearing use the same parameter of the oil film force. The bearing width L=10 mm, the viscosity η=0.04 Pa·s, the bearing diameter D=20 mm, the oil film gap c=0.3 mm, the bearing damping d=2×103 N/(m/s), and the eccentric magnitude is 0.1 mm.


### 3.1. The Influence of the Oil Film Force on the Dynamics Characteristic of Rotor System

In the role of oil film force, the system will show oil film instability phenomenon [21–25]. To clearly observe the dynamics characteristic of rotor system under the oil film force, then specially draw the three-dimensional spectrum on the right neck (node 25) of the rotor system without fault, shown in Figure 5. When the speed is 0≤ω<2830 r/min, the rotor system represents a synchronous positive precession. And then, there are twice frequency 2fr, three times frequency 3fr, and other higher-order times frequency with little peak besides the rotation frequency fr. While the speed is 2830≤ω≤3460 r/min, the eddy turbulence of the rotor system appears, and the frequency is nearly half the rotation frequency fr. When the rotational speed continues to increase to 3460<ω<4810 r/min, the oil film instability of the rotor is suppressed by the rapid increase in the frequency component, and the oil film vortex phenomenon disappears temporarily. There are fr,2fr, and 3fr on the spectrum at this time. When the speed is ω≥4810 r/min, the first-order oil film oscillation phenomenon appears, and the frequency is 27.5 Hz.Figure 5
The spectrum waterfall of the rotor system without fault.Through the evolution process analysis of the dynamics characteristic of rotor system supported by oil film force, and comparing the simulation result with the Muszynska experiment result in the literature [18], it is shown that the simulation code of the oil film force in this paper can fit the experiment result in the literature [18] well to some extent and verify the accuracy of the simulation code of the oil film force.


### 3.2. Dynamics Characteristic Analysis of Rotor System with Fault

The frequency spectrogram of rotor system with looseness-rubbing coupling fault under the nonlinear oil film force is shown in Figure6. From Figure 6(a), while the speed is a little low, the system under the oil film force shows some high frequency components which are not obvious. It is mainly the twice frequency 2fr and three times frequency 3fr. While the speed is above 4000 r/min, the situation of the rotor-stator rubbing is gradually serious and the friction radian increases. Under the increasing rubbing force, the phenomenon of oil film oscillation is depressed and the main character of the system is the rubbing fault. In the meanwhile, because of the existence of the base loosing fault, the half frequency 1/2fr, one-third frequency 1/3fr, and even a quarter of frequency 1/4fr appear in the vibration signal. The rotating frequency amplitude gradually decreases instead. The low-order harmonic becomes the dominant frequency instead of the rotating frequency.Figure 6
The frequency spectrogram of rotor system with coupling fault.
(a)
(b)The radial runout can be caused by the base loosing fault of the rotor system, and the low frequency vibration can be excited. In the meanwhile, the vibration of multitimes fundamental frequency exists under the rubbing force and the oil film force. However, the twice-fundamental frequency is dominant. The more serious the system fault is, the larger the proportion the frequency division and frequency doubling occupy. By studying the peak trend in the vibrating frequency of the looseness-rubbing coupling fault of rotor system which has the representative half frequency1/2fr, the rotating frequency fr, and the twice frequency 2fr, it is helpful to judge the severity of the fault. Figure 6(b) is the peak curves spectrogram with the changing speed of the coupling fault system which has the half frequency 1/2fr, the rotating frequency fr, and the twice frequency 2fr. Shown from the figure, while the speed is about 2500 r/min, the peak of the rotating frequency fr and the twice frequency 2fr reach the highest point. (fr is 8.233, 2fr is 3.176), and then there is no half frequency 1/2fr. The rotating frequency is the dominant frequency. While the speed increases to 3500 r/min, the half frequency 1/2fr appears. Its peak gradually surpasses the rotating frequency to be the dominant frequency and then reaches the top point (1/2fr is 6.953) near the speed of 4500 r/min. Subsequently it begins to decrease and the peak of the half frequency 1/2fr trends to 0 while the speed interval is between 6000 r/min and 8500 r/min. When the speed is above 8600 r/min, the half frequency 1/2fr increases fast. And then the rotating frequency fr and the twice frequency 2fr enter the growing step; however they do not surpass the frequency division to be the dominant frequency. From the total trend of the every vibrating frequency peak, the trend of the rotating frequency fr and the twice frequency 2fr are nearly the same and are both the decreasing trend with increasing the speed. However, from the trend of the half frequency 1/2fr, its speed quickly increases and surpasses the rotating frequency and the twice frequency on the speed interval 3000 r/min–6000 r/min and 8500 r/min–10000 r/min. Combined with Figure 6(a), on the step of the 1/2fr peak (above 6000 r/min) there are the 1/3 harmonic and 1/4 harmonic and their peaks are very large and surpass the rotating frequency. In total, the fault degree of rotor system is gradually serious with the speed increasing. The peaks of the rotating frequency and the twice frequency are gradually decreasing. A series of the frequency division harmonic gradually increases and its peak quickly increases. All the above illustrate that the constantly increasing speed leads to the rubbing fault, and then the system enters the step of coupling fault. The severity of the fault develops further.It is assumed that the eccentric mass only has an effect on the two rotary disks, and the rubbing fault happens on the right of the rotary disk which is on rotary disk 2. To distinguish the dynamics characteristic of looseness-rubbing coupling fault of rotor system under the support by oil film force, in the condition of foundation stiffnesskf3=1×109 N/m, the bolt tensile stiffness (loosing stiffness) kf2=8×105 N/m, and loosing gap δ=10 µm, Figures 7(a)–7(c) are the time-domain graph, orbit diagram, frequency spectrogram, and Poincare graph which are drawn by the steady response data on rotary disk 2 of the system with the speed ω being 3000 r/min, 4500 r/min, and 8000 r/min.Figure 7
Time-domain graph, orbit diagram, frequency spectrogram, and Poincare graph of rotor system with different rotational speed.
(a)
Rotational speedω=3000 r/min
(b)
Rotational speedω=4500 r/min
(c)
Rotational speedω=8000 r/minAs shown in Figure7, the time-domain waveforms of the rotor system produce distortion under the loosing stiffness and oil film force. And then the shape of time-domain gradually appears, and the lower part of the time-domain fluctuated shape is denser, while the upper part is relatively loose. The axis orbit is a semiellipse. With the speed increasing, the rotor-stator rubbing fault is gradually obvious. The axis orbit changes from a single semiellipse to many semiellipses which nest each other. Because of the influence on the eddy turbulence phenomenon, the twice frequency 2fr appears and then lots of frequency components appear. The half frequency 1/2fr, the 1/3 harmonic, 1/4 harmonic, and 1/5 harmonic appear in turn and the frequency division surpasses the rotation frequency to be the dominant frequency in the speed region. The Poincare graph gradually changes from one isolated point to two points with the speed increasing and then develops to five points. It unfolds the vibrating characteristic of frequency division of the system. The system is from periodic motion 1 to periodic motion 2 and then bifurcates to periodic motion 5.In total, because of the existence of the loosing fault in the rotor system under the support of the oil film force, the rubbing fault of the rotary disk and the casing is caused with the speed increasing, which leads to the fault degree of the rotor system gradually increased, a series of frequency division components are increasing, and the phenomenon of oil film oscillation is lagging.


### 3.3. Experiment Study of the Fault Rotor System

Build a fault rotor test-bed based on the simulation data above. Choose the portable analytic instrument BK3560D made in the B&K company to acquire and analyze the data. The sampling frequency is3.2×103 Hz. To prevent the eddy current sensor from disturbing each other, the test points on the x direction are 15 mm away from the test points on the y direction. The position of the test point is nearest to the number 25 node of the rotor system in Figure 2, so admit that the displacement sensor is installed at node 25 (right bearing of the model). Release the bolt on the bearing base a little, and then the loosing fault of the bolt appears. Pay attention to the loosing situation of the bolt to avoid the accident in the experiment. Choose the brass as the rubbing stator material and the C45 as the rotor material. The bearing type is made of oil-bearing graphite plain bearings. The rotor test-bed is shown in Figure 8.Figure 8
The schematic diagram of the rotor test-bed with fault.Through the experiment, the three-dimensional frequency spectrum near the test point of the right bearing on the vertical direction (y direction) can be obtained. Shown in Figure 9, there is only the rotation frequency in the low speed region. With the speed increasing, because of the existence of looseness fault, high power frequency division begins to appear. The first-order oil film oscillation appears in the following speed region, and then the oil film oscillation disappears. With the speed continuously increasing, the first-order oil film oscillation happens. While the speed continuously increases, the peak of the first-order oil film oscillation temporarily reduces and the second-order oil film oscillation happens. In the meanwhile, because of the existence of the looseness-rubbing fault, the proportion of the vibrating energy occupied by the frequency division and frequency doubling also increases.Figure 9
The frequency spectrum of the rotor test-bed with fault.Contrasting the combined frequency characteristics in every speed region acquired by the simulation with the result by the experiment, the comparison result is shown in Table2. From the table we can see that the test results and simulation results in the overall evolution of the trend are basically the same.Table 2
The combined frequency comparison of the simulation and the experiment.
Rotational speed/(r/min)
Frequency component/Hz
Simulation
100
≤
ω
<
990
f
r
1,000 ≤ω ≤ 2,750
f
r, 2fr, 3fr
2,750 <ω < 5,930
1/2fr, fr, 3/2fr, 2fr, 5/2fr, 3fr
5,930 ≤ω ≤ 7,980
f
n1, fr-fn1, fr, fr+fn1, 3/2fr, 2fr, 5/2fr, 3fr
7,980 <ω < 8,460
f
n1, 1/3fr, fr-fn1, 2/3fr, fr, fr+fn1, 4/3fr, 5/3fr, 2fr, 3fr, 1/3fr spectrum series, etc.
8,460 <ω < 9,520
1/3fr, 1/2fr, fr-fn1, 2/3fr, fr, fr+fn1, 4/3fr, 3/2fr, 5/3fr, 2fr, 3fr, 1/3fr spectrum series, 1/2fr spectrum series, etc.
9,520 <ω < 10,000
1/4fr, 1/3fr, 1/2fr, fr-fn1, 2/3fr, 3/4fr, fr, fr+fn1, 5/4fr, 4/3fr, 3/2fr, 5/3fr, 7/4fr, 2fr, 3fr, 1/4fr spectrum series, 1/3fr spectrum series, 1/2fr spectrum series, etc.
Experiment
100 ≤ω < 1,000
f
r
1,000 ≤ω ≤ 2,910
f
r, 2fr, 3fr
2,910 <ω < 6,070
1/2fr, fr, 3/2fr, 2fr, 5/2fr, 3fr
6,070 ≤ω ≤ 6,960
f
n1, fr, 3/2fr, 2fr, 5/2fr, 3fr
6,960 <ω < 7,980
1/2fr, fr, 3/2fr, 2fr, 5/2fr, 3fr
7,980 <ω < 8,580
f
n1, 1/3fr, 1/2fr, 2fn1, fn2, fr-fn1, 2/3fr, fr, fr+fn1, 4/3fr, 3/2fr, 5/3fr, 2fr, 3fr, 1/3fr spectrum series, 1/2fr spectrum series, etc.
8,580 <ω < 10,000
1/4fr, fn1, 1/3fr, 1/2fr, 2fn1, fn2, fr-fn1, 2/3fr, 3/4fr, fr, fr+fn1, 5/4fr, 4/3fr, 3/2fr, 5/3fr, 7/4fr, 2fr, 3fr, 1/4fr spectrum series, 1/3fr spectrum series, 1/2fr spectrum series, etc.The comparison of simulation and experiment shows that, in the simulation, the first- and second-order oil film oscillations occur after about 5000 r/min. The frequency of first-order oil film oscillation gradually replaces the rotation frequency as the main frequency. At the same time, the relevant combination frequency also appears. However, in the experiment, due to external interference and the presence of noise signals, the frequency division appears earlier, and the peak of frequency becomes higher. Meanwhile, the first- and second-order oil film oscillations are delayed, and the peak of the frequency is not obvious in many frequency components [26].


## 3.1. The Influence of the Oil Film Force on the Dynamics Characteristic of Rotor System

In the role of oil film force, the system will show oil film instability phenomenon [21–25]. To clearly observe the dynamics characteristic of rotor system under the oil film force, then specially draw the three-dimensional spectrum on the right neck (node 25) of the rotor system without fault, shown in Figure 5. When the speed is 0≤ω<2830 r/min, the rotor system represents a synchronous positive precession. And then, there are twice frequency 2fr, three times frequency 3fr, and other higher-order times frequency with little peak besides the rotation frequency fr. While the speed is 2830≤ω≤3460 r/min, the eddy turbulence of the rotor system appears, and the frequency is nearly half the rotation frequency fr. When the rotational speed continues to increase to 3460<ω<4810 r/min, the oil film instability of the rotor is suppressed by the rapid increase in the frequency component, and the oil film vortex phenomenon disappears temporarily. There are fr,2fr, and 3fr on the spectrum at this time. When the speed is ω≥4810 r/min, the first-order oil film oscillation phenomenon appears, and the frequency is 27.5 Hz.Figure 5
The spectrum waterfall of the rotor system without fault.Through the evolution process analysis of the dynamics characteristic of rotor system supported by oil film force, and comparing the simulation result with the Muszynska experiment result in the literature [18], it is shown that the simulation code of the oil film force in this paper can fit the experiment result in the literature [18] well to some extent and verify the accuracy of the simulation code of the oil film force.


## 3.2. Dynamics Characteristic Analysis of Rotor System with Fault

The frequency spectrogram of rotor system with looseness-rubbing coupling fault under the nonlinear oil film force is shown in Figure6. From Figure 6(a), while the speed is a little low, the system under the oil film force shows some high frequency components which are not obvious. It is mainly the twice frequency 2fr and three times frequency 3fr. While the speed is above 4000 r/min, the situation of the rotor-stator rubbing is gradually serious and the friction radian increases. Under the increasing rubbing force, the phenomenon of oil film oscillation is depressed and the main character of the system is the rubbing fault. In the meanwhile, because of the existence of the base loosing fault, the half frequency 1/2fr, one-third frequency 1/3fr, and even a quarter of frequency 1/4fr appear in the vibration signal. The rotating frequency amplitude gradually decreases instead. The low-order harmonic becomes the dominant frequency instead of the rotating frequency.Figure 6
The frequency spectrogram of rotor system with coupling fault.
(a)
(b)The radial runout can be caused by the base loosing fault of the rotor system, and the low frequency vibration can be excited. In the meanwhile, the vibration of multitimes fundamental frequency exists under the rubbing force and the oil film force. However, the twice-fundamental frequency is dominant. The more serious the system fault is, the larger the proportion the frequency division and frequency doubling occupy. By studying the peak trend in the vibrating frequency of the looseness-rubbing coupling fault of rotor system which has the representative half frequency1/2fr, the rotating frequency fr, and the twice frequency 2fr, it is helpful to judge the severity of the fault. Figure 6(b) is the peak curves spectrogram with the changing speed of the coupling fault system which has the half frequency 1/2fr, the rotating frequency fr, and the twice frequency 2fr. Shown from the figure, while the speed is about 2500 r/min, the peak of the rotating frequency fr and the twice frequency 2fr reach the highest point. (fr is 8.233, 2fr is 3.176), and then there is no half frequency 1/2fr. The rotating frequency is the dominant frequency. While the speed increases to 3500 r/min, the half frequency 1/2fr appears. Its peak gradually surpasses the rotating frequency to be the dominant frequency and then reaches the top point (1/2fr is 6.953) near the speed of 4500 r/min. Subsequently it begins to decrease and the peak of the half frequency 1/2fr trends to 0 while the speed interval is between 6000 r/min and 8500 r/min. When the speed is above 8600 r/min, the half frequency 1/2fr increases fast. And then the rotating frequency fr and the twice frequency 2fr enter the growing step; however they do not surpass the frequency division to be the dominant frequency. From the total trend of the every vibrating frequency peak, the trend of the rotating frequency fr and the twice frequency 2fr are nearly the same and are both the decreasing trend with increasing the speed. However, from the trend of the half frequency 1/2fr, its speed quickly increases and surpasses the rotating frequency and the twice frequency on the speed interval 3000 r/min–6000 r/min and 8500 r/min–10000 r/min. Combined with Figure 6(a), on the step of the 1/2fr peak (above 6000 r/min) there are the 1/3 harmonic and 1/4 harmonic and their peaks are very large and surpass the rotating frequency. In total, the fault degree of rotor system is gradually serious with the speed increasing. The peaks of the rotating frequency and the twice frequency are gradually decreasing. A series of the frequency division harmonic gradually increases and its peak quickly increases. All the above illustrate that the constantly increasing speed leads to the rubbing fault, and then the system enters the step of coupling fault. The severity of the fault develops further.It is assumed that the eccentric mass only has an effect on the two rotary disks, and the rubbing fault happens on the right of the rotary disk which is on rotary disk 2. To distinguish the dynamics characteristic of looseness-rubbing coupling fault of rotor system under the support by oil film force, in the condition of foundation stiffnesskf3=1×109 N/m, the bolt tensile stiffness (loosing stiffness) kf2=8×105 N/m, and loosing gap δ=10 µm, Figures 7(a)–7(c) are the time-domain graph, orbit diagram, frequency spectrogram, and Poincare graph which are drawn by the steady response data on rotary disk 2 of the system with the speed ω being 3000 r/min, 4500 r/min, and 8000 r/min.Figure 7
Time-domain graph, orbit diagram, frequency spectrogram, and Poincare graph of rotor system with different rotational speed.
(a)
Rotational speedω=3000 r/min
(b)
Rotational speedω=4500 r/min
(c)
Rotational speedω=8000 r/minAs shown in Figure7, the time-domain waveforms of the rotor system produce distortion under the loosing stiffness and oil film force. And then the shape of time-domain gradually appears, and the lower part of the time-domain fluctuated shape is denser, while the upper part is relatively loose. The axis orbit is a semiellipse. With the speed increasing, the rotor-stator rubbing fault is gradually obvious. The axis orbit changes from a single semiellipse to many semiellipses which nest each other. Because of the influence on the eddy turbulence phenomenon, the twice frequency 2fr appears and then lots of frequency components appear. The half frequency 1/2fr, the 1/3 harmonic, 1/4 harmonic, and 1/5 harmonic appear in turn and the frequency division surpasses the rotation frequency to be the dominant frequency in the speed region. The Poincare graph gradually changes from one isolated point to two points with the speed increasing and then develops to five points. It unfolds the vibrating characteristic of frequency division of the system. The system is from periodic motion 1 to periodic motion 2 and then bifurcates to periodic motion 5.In total, because of the existence of the loosing fault in the rotor system under the support of the oil film force, the rubbing fault of the rotary disk and the casing is caused with the speed increasing, which leads to the fault degree of the rotor system gradually increased, a series of frequency division components are increasing, and the phenomenon of oil film oscillation is lagging.


## 3.3. Experiment Study of the Fault Rotor System

Build a fault rotor test-bed based on the simulation data above. Choose the portable analytic instrument BK3560D made in the B&K company to acquire and analyze the data. The sampling frequency is3.2×103 Hz. To prevent the eddy current sensor from disturbing each other, the test points on the x direction are 15 mm away from the test points on the y direction. The position of the test point is nearest to the number 25 node of the rotor system in Figure 2, so admit that the displacement sensor is installed at node 25 (right bearing of the model). Release the bolt on the bearing base a little, and then the loosing fault of the bolt appears. Pay attention to the loosing situation of the bolt to avoid the accident in the experiment. Choose the brass as the rubbing stator material and the C45 as the rotor material. The bearing type is made of oil-bearing graphite plain bearings. The rotor test-bed is shown in Figure 8.Figure 8
The schematic diagram of the rotor test-bed with fault.Through the experiment, the three-dimensional frequency spectrum near the test point of the right bearing on the vertical direction (y direction) can be obtained. Shown in Figure 9, there is only the rotation frequency in the low speed region. With the speed increasing, because of the existence of looseness fault, high power frequency division begins to appear. The first-order oil film oscillation appears in the following speed region, and then the oil film oscillation disappears. With the speed continuously increasing, the first-order oil film oscillation happens. While the speed continuously increases, the peak of the first-order oil film oscillation temporarily reduces and the second-order oil film oscillation happens. In the meanwhile, because of the existence of the looseness-rubbing fault, the proportion of the vibrating energy occupied by the frequency division and frequency doubling also increases.Figure 9
The frequency spectrum of the rotor test-bed with fault.Contrasting the combined frequency characteristics in every speed region acquired by the simulation with the result by the experiment, the comparison result is shown in Table2. From the table we can see that the test results and simulation results in the overall evolution of the trend are basically the same.Table 2
The combined frequency comparison of the simulation and the experiment.
Rotational speed/(r/min)
Frequency component/Hz
Simulation
100
≤
ω
<
990
f
r
1,000 ≤ω ≤ 2,750
f
r, 2fr, 3fr
2,750 <ω < 5,930
1/2fr, fr, 3/2fr, 2fr, 5/2fr, 3fr
5,930 ≤ω ≤ 7,980
f
n1, fr-fn1, fr, fr+fn1, 3/2fr, 2fr, 5/2fr, 3fr
7,980 <ω < 8,460
f
n1, 1/3fr, fr-fn1, 2/3fr, fr, fr+fn1, 4/3fr, 5/3fr, 2fr, 3fr, 1/3fr spectrum series, etc.
8,460 <ω < 9,520
1/3fr, 1/2fr, fr-fn1, 2/3fr, fr, fr+fn1, 4/3fr, 3/2fr, 5/3fr, 2fr, 3fr, 1/3fr spectrum series, 1/2fr spectrum series, etc.
9,520 <ω < 10,000
1/4fr, 1/3fr, 1/2fr, fr-fn1, 2/3fr, 3/4fr, fr, fr+fn1, 5/4fr, 4/3fr, 3/2fr, 5/3fr, 7/4fr, 2fr, 3fr, 1/4fr spectrum series, 1/3fr spectrum series, 1/2fr spectrum series, etc.
Experiment
100 ≤ω < 1,000
f
r
1,000 ≤ω ≤ 2,910
f
r, 2fr, 3fr
2,910 <ω < 6,070
1/2fr, fr, 3/2fr, 2fr, 5/2fr, 3fr
6,070 ≤ω ≤ 6,960
f
n1, fr, 3/2fr, 2fr, 5/2fr, 3fr
6,960 <ω < 7,980
1/2fr, fr, 3/2fr, 2fr, 5/2fr, 3fr
7,980 <ω < 8,580
f
n1, 1/3fr, 1/2fr, 2fn1, fn2, fr-fn1, 2/3fr, fr, fr+fn1, 4/3fr, 3/2fr, 5/3fr, 2fr, 3fr, 1/3fr spectrum series, 1/2fr spectrum series, etc.
8,580 <ω < 10,000
1/4fr, fn1, 1/3fr, 1/2fr, 2fn1, fn2, fr-fn1, 2/3fr, 3/4fr, fr, fr+fn1, 5/4fr, 4/3fr, 3/2fr, 5/3fr, 7/4fr, 2fr, 3fr, 1/4fr spectrum series, 1/3fr spectrum series, 1/2fr spectrum series, etc.The comparison of simulation and experiment shows that, in the simulation, the first- and second-order oil film oscillations occur after about 5000 r/min. The frequency of first-order oil film oscillation gradually replaces the rotation frequency as the main frequency. At the same time, the relevant combination frequency also appears. However, in the experiment, due to external interference and the presence of noise signals, the frequency division appears earlier, and the peak of frequency becomes higher. Meanwhile, the first- and second-order oil film oscillations are delayed, and the peak of the frequency is not obvious in many frequency components [26].


## 4. Conclusions

(1) In the healthy rotor system under the support of the sliding bearing (oil film force), with the speed increasing, the phenomenon of oil film oscillation gradually appears. A series of harmonic components of the frequency division and frequency doubling are very rich. And the frequency division gradually becomes the dominant frequency instead of the rotation frequency. At the same time, the phenomenon of oil film oscillation intensifies.(2) While the looseness-rubbing coupling fault of rotor system happens, because of the loosing stiffness and rubbing force, the phenomenon of unstable oil film is depressed. The appearance of the first- and second-order oil film oscillation phenomenon is delayed. In the meanwhile, the rubbing phenomenon of the system intensifies and gradually bifurcates from cycle 1 motion to cycleN motion with the speed increasing.(3) Through studying the dynamics characteristic of the rotor system in different speed condition, according to the research on the dynamics characteristic of the rotor system with different rotational speed, it was found that the looseness-rubbing coupling fault is often characterized by rubbing fault, and the lower part of the time-domain fluctuated shape is denser, while the upper part is relatively loose, and multiple nested half ellipses were shown in orbit diagram. It could be used as a theoretical basis for diagnosing looseness-rubbing coupling fault of rotor-sliding bearing system.(4) In this paper, the nonlinear dynamic behavior of the coupled fault flexible rotor system supported by the sliding bearing is studied. It is the first time to reveal the relationship between the oil film instability and the rotating speed of the rotor system. It could be found that when the speed reached the certain value, the oil film oscillation occurred. And the frequency of the oil film instability is submerged by the frequency division component caused by the fault. It is not the same as healthy rotor system. And the frequency component of rotor system is dominated by the fault frequency.



---

*Source: 8742468-2017-12-05.xml*