Published: 17 October 2022

Rubbing vibration characteristics of double-rotor system under wave load

Qiang Wang1
Zhikai Xing2
Yongbao Liu3
Mo Li4
1, 2, 3, 4College of Power Engineering, Naval University of Engineering, Wuhan, 430000, China
Corresponding Author:
Zhikai Xing
Views 96
Reads 34
Downloads 357

Abstract

Structural characteristics of a certain type of gas turbine double rotor and its hull are analyzed. The wave of a regular form is considered, and the wave load is simplified to the form of sinusoidal force. The rubbing force and nonlinear supporting force of the bearing are deduced as per the Hertz contact theory and Coulomb friction law. The mechanical model of the double rotor system under the combined action of the wave load and rotor friction impact is established, which is solved by the fourth-order Runge-Kutta method. The bifurcation diagram, phase diagram, Poincare section diagram, axis trajectory diagrams, effects of rotor eccentricity, and friction coefficient on the nonlinear dynamic characteristics of the system under rub impact fault are studied. The obtained results show that with an increase in eccentricity, the system enters into a chaotic state through the period doubling and paroxysmal bifurcations, and successively experiences nonlinear behaviors such as periodic 2 bifurcation, periodic 4 bifurcation, and quasi periodic and chaos. In the double rotor system, the nonlinear dynamic behavior of the single period and semi Nermark-Sacker bifurcation occurs with increasing friction coefficient. The present findings can provide a theoretical basis for condition monitoring, fault diagnosis, and design optimization of this type of gas turbine double rotor system.

1. Introduction

With the development of the economy, gas turbine, whose operational structure is mainly a double rotor system, is more and more widely used in the marine field. Considering the performance and economy of gas turbine, the radial clearance between the rotor blade and casing is designed to be smaller and smaller. Research shows that reducing the tip clearance among the compressor blade, turbine blade and casing can effectively reduce the air intake loss and fuel consumption, and improve the work efficiency of the compressor and turbine. The structure of a certain type of ship has the characteristics of small tonnage, high speed, small structural stiffness, and flat ship shape. When sailing at a high speed, the pitch and heave cause periodic slamming between the waves and bow, resulting in hull flutter. At the same time, the gas turbine rotor rotates at a high speed for a long time in the environment of high temperature and high pressure, which may easily cause the imbalance fault of the rotor, aggravate the vibration of the gas turbine rotor, and cause the rub impact fault between the rotor and the gearbox. Therefore, this paper takes the gas turbine double rotor system as the research object. The nonlinear dynamic behavior of the double rotor system under wave load is studied, which provides a reference for the design of gas turbine.

Scholars have conducted a lot of research on the rotor rub impact. In reference [1], some models of the rub impact force were established based on the Coulomb friction law and generalized Coulomb law by considering cubic dependence of the relative slip speed and the effects of two different models on the bifurcation and chaos characteristics of the rotor. The phenomena of period doubling bifurcation, Hopf bifurcation and saddle point bifurcation of the Jeffcott rotor during rub impact were studied in reference [2]. Literature [3] studied the bifurcation and chaos characteristics of the Jeffcott rotor when the rub impact fault occurs, and three ways to cause chaos were found. In reference [4], the nonlinear behavior of grazing bifurcation into chaos was studied when the friction between the disk and shell was caused by increasing rotor speed. In reference [5], considering the rotor static rub impact caused by sudden impact excitation and sudden unbalanced excitation, the nonlinear response of the rotor system under different rub impact forms were studied by using finite element method. Literature [6] studied the friction impact between the rotor and stator caused by blade falling off, and analyzed the important role of the friction coefficient and diaphragm model. Reference [7] proposed a new Poincare section method that can identify the nonlinear behavior in the process of rub impact more accurately. The influence of bearing waviness on the nonlinear dynamics of a low pressure gas turbine rotor system was studied in reference [8]. Reference [9] studied the effects of speed and bearing clearance on the nonlinear response of three supported double span rotor bearing system. According to the linear contact force and Coulomb friction force, a model of friction force was established in reference [10], and a new continuous shooting method was proposed to get the way for the system to enter and leave chaos. In reference [11], a sudden unbalance was added to the double rotor system, resulting in rub impact-θ, where a new way of entering into chaos through torus intersection was found. In reference [12], a finite element dynamic model of steam turbine rotor was established by considering the effect of air flow excitation force to study the changes of the nonlinear response of bearing when local rub impact occurs at different positions. In reference [13], the effects of oil film force and rub impact force on the bifurcation and chaos characteristics of single span double disk rub impact rotor system were studied by a numerical method, and experiments were designed to verify the correctness of the model. Reference [14] established a dynamic model of the double rotor system by considering the rub impact fault, studied the nonlinear dynamic response characteristics of the system under three different conditions: slight rub impact, heavy rub impact, and serious rub impact. Some achievements have been made by M. Sayed, et al in reference [15, 16].

In the dynamic modeling of the rotor system in the above research, the influence of the actual working environment of the rotor was ignored, and the bearing was considered as a linear support with stiffness and damping so the model was not practical enough. Therefore, in this paper, the wave load is simplified as a simple harmonic force acting on the double rotor system, and the nonlinear bearing force of the bearing and the rub impact force of the rotor are deduced according to the elastic Hertz theory. The dynamic model of the double rotor system with wave load rotor rub impact coupling is also established, and the bifurcation and chaos characteristics of the system when the rotor eccentricity changes are studied. The present findings provide a theoretical basis for condition monitoring and fault diagnosis of the actual gas turbine rotor rub impact fault.

2. Dynamic model of double rotor system

2.1. Model diagram

In this paper, the dual rotor system of a marine gas turbine is taken as the research object. The dual rotor system model is obtained by the centralized mass method as shown in Fig. 1, where Ob1, Ob2, Ob3 and Ob4 are the centroids of bearings 1, 2, 3 and 4, respectively. Op1, Op2 and Op3 are the centroids of discs 1, 2 and 3, respectively. Oc1, Oc2 and Oc3 are the centroids of discs 1, 2 and 3, respectively. e1, e2 and e3 are the eccentricities of discs 1, 2 and 3, respectively. Assuming that the structural materials of shafts 1 and 2 are the same, the stiffness relationship between each shaft section is obtained according to the plane beam structure relationship, where k1=k/125, k2=4k, k5=k/2 and k3=k4=k6=k. The inner and outer rotors are connected through intermediate bearing 2. Bearings 1, 3 and 4 connect the rotor with the gearbox. See Table 1 for structural parameters of double rotors and Table 2 for bearing parameters.

Fig. 1Schematic diagram of rubbing fault of double rotors

Schematic diagram of rubbing fault of double rotors

Table 1Structural parameters of double rotors

Parameter name
Parameter symbol
Dimension
Bearing quality
mb1, mb2, mb3, mb4
Kg
Quality of disc
mp1, mp2, mp3
Kg
Bearing damping
cb1, cb2, cb3, cb4
N∙s/m

Table 2Parameters of bearings

Parameter
Dimension
bearing 1
bearing 2
bearing 3
bearing 4
Internal diameter Ri
mm
22
25
40
30
External diameter Ro
mm
50
47
64
62
Number of rolling elements Z
number
7
15
8
12
Radial clearance Gr
m
12×10-6
12×10-6
12×10-6
12×10-6
Contact stiffness kb
N/m3/2
13.34×109
3×108
7.055×109
7.055×109

2.2. Model of external load, bearing force and rubbing force

Waves can be divided into regular waves and irregular waves. The waveform contours of regular waves can be expressed by a cosine function as shown in Fig. 2.

Fig. 2Schematic diagram of regular waves

Schematic diagram of regular waves

The wave surface of cosine wave propagating along O0x0 direction can be expressed by Eq. (1):

1
ξ=ξAcosk0x0-ωt,

where, ξ is the rise of wave surface; ξA is the wave amplitude that refers to the distance between the wave crest and trough from the horizontal plane with the wave height ξW=2ξA; k0 is the wave number equal to 2π times the number of waves appearing on the unit distance; ω is the wave frequency; α is the wave inclination that indicates the inclination of the wave surface; α0 is the maximum wave inclination which is conventionally called the wave inclination; and λ represents the wavelength equal to the distance between two adjacent peaks and troughs. Yao [15] studied the calculation method for the load of the fully padded hovercraft under regular waves, and deduced the dynamic pressure of the hovercraft as in Eq. (2):

2
Pt=Asinωet+Bsinωet,

where A and B are the matrices related to the geometric characteristics of the ship and wave parameters, respectively, and ωe is the wave frequency. In this paper, the wave surface is considered as the ideal state of a regular wave, and the wave load is simplified as a sinusoidal force. If the period is 1 s and the maximum amplitude of load is 100 N, and the action of the wave load only along the Y direction is considered, then Fvy1=Fvy2=Fvy3=100sin2πt. It acts on the three support bearing positions connected with the gearbox of the double rotor system as shown in Fig. 1.

In this paper, the influence of the bearing oil film force on the system is ignored. According to the Hertz contact theory, the elastic contact force between the bearing rolling element and raceway is considered to reflect the influence of the bearing nonlinear force on the system. Fig. 3 shows the rolling bearing model, where Ri is the inner diameter of the bearing, Ro is the outer diameter of the bearing, Gr represents the radial clearance of the bearing, ωi is the inner ring speed, ωo is the outer ring speed, j represents the jth rolling element, θj is the position angle of the jth rolling body, and the coordinates after the motion of the axis diameter centroid are x and y.

Fig. 3Model diagram of rolling bearing

Model diagram of rolling bearing

The common angular speed of the rolling element:ωc=ωiRi+ωoRo/Ri+Ro, angular position of the jth rolling element: θj=ωct+2πZj-1, j=1,2,3Z, and contact deformation of the jth rolling element: χj=xcosθj+ysinθj-Gr. A contact can produce a positive contact force only. When the contact deformation is χj0, there will be no pressure. Therefore, the Heaviside function is introduced as follows: Hχj=1, χj>0,0, χj0.

The nonlinear bearing force of a rolling bearing is deduced according to the Hertz contact theory, and it is expressed by Eq. (3), where BK is the nonlinear Hertz contact stiffness, n= 3/2 for ball bearings 1 and 4, and N= 10/9 for roller bearings 2 and 3:

3
fbx=j=1Zfj,x=j=1Zfjcosθj=kbj=1ZHχjχjncosθj,fby=j=1Zfj,y=j=1Zfjsinθj=kbj=1ZHχjχjnsinθj.

Fig. 4The friction model

The friction model

The rub impact model is shown in Fig. 4. The deformation relationship caused by rub impact is treated according to the contact deformation of the disc and circular hole with certain width and thickness, which simplifies the problem to the contact relationship of two inscribed circles. The collision force acts along the common normal line of the inscribed circles, and the friction force is perpendicular to it. According to the Hertz contact theory, there is a 3/2 nonlinear relationship between the elastic force generated by collision deformation and the amount of deformation. Local elastic deformation takes place during the collision and extrusion of two inscribed circles. The collision force FN is obtained as FN=0, δ0,δa32, δ<0,. The friction between the rotating stators meets the Coulomb’s law of friction. The friction coefficient is a fixed parameter and it is independent of relative speed. The relation Fr=fFN of friction Fr is obtained. According to the assumptions, the friction force can be expressed as in Eq. (4):

4
PxPy=FNr-1-ff-1xy,
a=0.89161-μs2Es+1-μr2Er31Rs+1Rr13, δ=δ0-r, r=x2+y2,

where μs, Es, Rs, μr, Er and Rr are the Poisson’s ratio, elastic modulus, and radius of the stator and rotor, respectively. f is the friction coefficient between the rotor and stator. Px and Py are the friction forces along the X direction and Y direction, respectively. δ0 is the initial clearance between the rotor and stator, r is the radial displacement of the rotor, and δ is the elastic deformation of the rotating static contact. The parameters are as follows: the compressor blade material is titanium alloy TC4, whose elastic modulus Er= 1.05×1011N/m2 and Poisson’s ratio μs= 0.34. The casing material is cast steel, whose elastic modulus is Es= 2×1011 N/m2, Poisson’s ratio μs= 0.3, and Rr= 0.1 m.

2.3. Differential equation of motion of system

The dynamic model of marine gas turbine double rotor system established in this paper has 7 lumped masses and 14 degrees of freedom. Assuming that the rotating shaft is an elastic shaft without mass, only the transverse vibration of the rotor is considered by ignoring both torsional vibration and axial vibration. Accordingly, the differential equation of motion of the system is derived from the Lagrange equation of motion as expressed by Eq. (5):

5
MX¨+CX˙+KX=Fbx+Fex+Px,MY¨+CY˙+KY=Fby+Fey+Fvy+Py-Mg,

where:

K=k1-k100000-k1k1+k2+k6-k2000-k60-k2k2+k3-k300000-k3k30000000k4-k400000-k4k4+k5-k50-k6000-k5k5+k6,
M=diagmb1,mb2,mp1,mb3,mb4,mp2,mp3,C=diagcb1,cb2,cp1,cb3,cb4,cp2,cp3,X=xb1,xb2,xp1,xb3,xb4,xp2,xp3T,Y=yb1,yb2,yp1,yb3,yb4,yp2,yp3T,Fbx=Fb1x,Fb2x,0,Fb3x,Fb4x,0,0T,Fby=Fb1y,Fb2y,0,Fb3y,Fb4y,0,0T,Fvy=Fvy1,0,0,Fvy3,Fvy2,0,0T,Px=0,0,0,0,0,px,0T,Py=0,0,0,0,0,py,0T,Fex=0,0,mp1e1ω12cosω1t,0,0,mp2e2ω22cosω2t,mp3e3ω22cosω2tT,Fey=0,0,mp1e1ω12sinω1t,0,0,mp2e2ω22sinω2t,mp3e3ω22sinω2tT.

In Eq. (5), M is the system mass matrix, C is the system damping matrix, K is the system stiffness matrix, Fbx and Fby are the bearing force matrices, Fex and Fey are the rotor unbalance force matrices, Fvy is the wave load simple harmonic force matrix, and Px and Py are the rub impact force matrices. Further, ω1 is the speed of rotor 1, ω2is the speed and α=ω2/ω1is the peed ratio of the rotor. Dimensionless parameter ξ is the reference value of the dimensionless length. Other introduced dimensionless parameters are as follows:

t=w1t, Xi=xix, Yi=yix, X˙i=x˙ixwi, Y˙i=y˙ixwi, X¨i=x¨ixwi2, Y¨i=y¨ixwi2, Ki=kimiwi2,Kb=kbxhmiwi2, Ci=cimiw, G=gxw2, Ei=eix, P~x=Pxξmp2ω2, P~y=Pyξmp2ω2.

For ball bearings, H= 1/2; and for cylindrical roller bearings, H= 1/9. According to the dimensionless parameters and Eq. (5), the dimensionless equations are obtained as expressed by Eq. (6):

6
X~¨+CX~˙Mω+KXMω2=F~bx+F~ex+P~xMγω2,Y~¨+CY~˙Mω+KYMω2=F~by+F~ey+F~vy+P~y-MgMγω2.

According to the high-speed operational characteristics of the rotor of a gas turbine generator, reasonable simplification is carried out to facilitate the modeling of the rub impact fault of the rotor system. The following assumptions are made: (1) nonlinear elastic collision occurs at the contact position when the rotating stator collides with friction; (2) the friction satisfies the Coulomb's law of friction; (3) the thermal effect caused by friction is ignored; and (4) the damping effect in the process of collision deformation is ignored. Based on the above assumptions, the rub impact force model of gas turbine and gas generator rotor is established as per the Hertz contact theory and Coulomb friction law.

3. Bifurcation and chaos characteristics of two rotor system under rub impact fault

The double rotor structure characteristics of a marine gas turbine are analyzed. The rotor speed of the gas generator is high, and the blade length of the compressor is greater than that of the turbine. When the rotor system vibration is too large, the compressor is more likely to collide with the casing, making the nonlinear characteristics of the system more complex. Therefore, the compressor rotor and disk 2 are also taken as the research object, and the influence of rotor unbalance and friction coefficient on the nonlinear characteristics of the system is analyzed.

3.1. Influence of eccentricity on Bifurcation and chaos characteristics of double rotor system

Rotor unbalance is one of the major reasons for abnormal vibration and excessive response amplitude of the rotor, which leads to the rub impact fault between the rotor and stator. It has a significant impact on the nonlinear response characteristics of the double rotor system. Various parameter values are set here as follows: speed of rotor 2, ω2= 1680 rad/s; initial clearance, δ0= 20 μm; and friction coefficient f= 0.1. Taking the eccentricity ep2 of disk 2 as the bifurcation parameter, ep2 is studied in the range of 0 to 30 μm. When m is changed, the bifurcation and chaos characteristics of the system also change.

A bifurcation diagram of the variation of the horizontal displacement Xp2 of disc 2 with eccentricity ep2 is shown in Fig. 5, where it can be seen that disc 2 exhibits rich nonlinear dynamic behavior with increasing eccentricity. When the eccentricity is small, disk 2 is in a chaotic state. When the eccentricity is increased to 7 μm, disk 2 enters into cycle 2 through inverted bifurcation and moves at 10 μm to enter into the chaotic state through the period doubling bifurcation road. When the eccentricity continues to increase to 13 μm, disk 2 enters into the state of periodic four motion through the inverted bifurcation path, and then enters into the chaos through the period doubling bifurcation path. When the eccentricity is increased to 20 μm, disk 2 enters into the quasi periodic motion through the inverted bifurcation path, and then enters into the chaos through the paroxysmal bifurcation path. The bifurcation and chaos characteristics of disk 2 are analyzed by selecting different values of eccentricity.

When eccentricity ep2 = 10 μm as shown in Fig. 6, the phase trajectory of disk 2 forms two closed circular orbits, which are mapped into two isolated point sets through Poincare. The axis trajectory is in the shape of "8". In the spectrum diagram, it is mainly frequency conversion and 1/2 frequency doubling. At this moment, disk 2 moves in cycle two. When eccentricity ep2 = 13 μM for disc 2 as shown in Fig. 7, its phase trajectory forms four closed orbits. Due to the action of rubbing, the amplitude of disc 2 becomes 20 μm. The “cutting” phenomenon occurs at M. Poincare is mapped into four isolated points. The conversion frequency is the main excitation frequency of disk 2. There are also 0.25 octave, 0.5 octave, 0.75 octave, 1.25 octave and 1.5 octave. At this moment, disk 2 makes periodic four motion. When eccentricity ep2 = 20 μM for disc 2 as shown in Fig. 8, its phase trajectory produces a “cutting” phenomenon due to the friction impact. Poincare maps to a certain section to form an unsealed semicircular curve, and the axis trajectory presents a regular state of periodic motion. In the spectrum diagram, frequency conversion is the main, with 2-fold, 3-fold and 5-fold frequency. At this moment, disk 2 makes quasi periodic motion.

Fig. 5Bifurcation diagram of disk 2 with eccentricity

Bifurcation diagram of disk 2 with eccentricity

Fig. 6Nonlinear dynamic characteristic diagram of disk 2 at ep2= 10 μm

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 10 μm

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 10 μm

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 10 μm

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 10 μm

d) Spectrum diagram

Fig. 7Nonlinear dynamic characteristic diagram of disk 2 at ep2= 13 μm

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 13 μm

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 13 μm

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 13 μm

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 13 μm

d) Spectrum diagram

Fig. 8Nonlinear dynamic characteristic diagram of disk 2 at ep2= 20 μm

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 20 μm

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 20 μm

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 20 μm

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at ep2= 20 μm

d) Spectrum diagram

3.2. Influence of friction coefficient on Bifurcation and chaos characteristics of double rotor system

By analyzing the rub impact force model, the friction generated by rub impact is determined by the contact displacement and friction coefficient between the rotor and stator. When rub impact occurs between the rotor and stator, the material characteristics of the rotor-stator contact surface change due to the thermal effect and other factors, and the friction coefficient also changes. Various parameters are set here as follows: speed of rotor 2, ω2= 1000 rad/s; initial clearance δ0= 30 μm; and ep2= 10 μm. Taking the friction coefficient f between the rotor and stator as the bifurcation parameter, the influence of the change in f on the bifurcation and chaos characteristics of the system is studied.

As shown in Fig. 9, the horizontal displacement of disk 2 changes with the friction coefficient. When the amplitude of disk 2 is too large, it will rub with the stator. At this moment, the “cutting” phenomenon occurs in the bifurcation diagram. It is seen in the figure that the displacement response Xp2 of disk 2 is only half, and the vibration amplitude of the other half cannot continue to increase due to the rubbing. When the friction coefficient is small, the rubbing force has little effect on the nonlinear characteristics of disc 2. With increasing friction coefficient, the bifurcation and chaos characteristics of disk 2 change, and the nonlinear behavior becomes more complex.

Fig. 9Bifurcation diagram of disk 2 with friction coefficient

Bifurcation diagram of disk 2 with friction coefficient

When the friction coefficient is small, the friction effect is also small, and it will not have a significant impact on the nonlinear response of the system. As shown in Fig. 10, when the friction coefficient is 0.01, the nonlinear dynamic characteristics of disc 2 are the same as those at the same speed as stated in section 3.4.1, and they are in the state of a single cycle motion. At this moment, the impact of friction is small. When the friction coefficient is 0.1 as shown in Fig. 11, disc 2 is in the state of a quasi-periodic motion. When the friction coefficient is 0.2 as shown in Fig. 12, the friction is large, the impact of the rub impact fault is enhanced, the phase trajectory in the phase diagram of disk 2 forms a more complex track pattern, the scattered points on the Poincare section tend to be irregularly distributed, the axis trajectory gets the characteristics of periodic motion, the transfer frequency in the spectrum diagram is the main frequency, and the 3.5 times frequency is the secondary frequency. At this moment, disk 2 tends to transform from quasi periodic motion to chaos. When the friction coefficient is 0.25 as shown in Fig. 13, the phase trajectories of disc 2 are staggered, the points in Poincare section are scattered, the axis trajectory is in the state of an irregular motion, the transfer frequency in the spectrum is the main frequency, and a continuous frequency peak appears near the third harmonic frequency, indicating that the degree of the rub impact fault is serious at this moment, causing disc 2 to enter into a chaotic state.

Fig. 10Nonlinear dynamic characteristic diagram of disk 2 at f= 0.01

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.01

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.01

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.01

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.01

d) Spectrum diagram

Fig. 11Nonlinear dynamic characteristic diagram of disk 2 at f= 0.1

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.1

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.1

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.1

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.1

d) Spectrum diagram

Fig. 12Nonlinear dynamic characteristic diagram of disk 2 at f= 0.2

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.2

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.2

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.2

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.2

d) Spectrum diagram

Fig. 13Nonlinear dynamic characteristic diagram of disk 2 at f= 0.25

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.25

a) Phase diagram

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.25

b) Poincare section

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.25

c) Axis trajectory

Nonlinear dynamic characteristic diagram of disk 2 at f= 0.25

d) Spectrum diagram

4. Conclusions

In this paper, the rub impact problem of rotor-stator is transformed into the contact problem of two inscribed circles. The rub impact force model is deduced according to the Hertz contact theory and Coulomb friction law. The mechanical model of a double rotor system is established by coupling wave load with the rub impact fault of the rotor. It is found that when the eccentricity of disk is small, disk 2 is in a chaotic state. When the eccentricity is increased to 10 μm, disk 2 enters into cycle 2 through inverted bifurcation, and moves at 10 μm enters the chaotic state through the period doubling bifurcation road. If the eccentricity continues to increase to 13 μm, disk 2 enters into the state of periodic four motion through the inverted bifurcation path, and then enters into chaos through the period doubling bifurcation path. Upon further increasing the eccentricity to 20 μm, disk 2 enters into the quasi periodic motion through the inverted bifurcation path, and then enters into chaos through the paroxysmal bifurcation path. These two paths to chaos can provide a theoretical basis for condition monitoring and fault diagnosis of a double rotor system.

References

  • N. Vlajic, A. R. Champneys, and B. Balachandran, “Nonlinear dynamics of a Jeffcott rotor with torsional deformations and rotor-stator contact,” International Journal of Non-Linear Mechanics, Vol. 92, pp. 102–110, Jun. 2017, https://doi.org/10.1016/j.ijnonlinmec.2017.02.002
  • S.-K. Choi and S. T. Noah, “Mode-locking and chaos in a Jeffcott rotor with bearing clearances,” Journal of Applied Mechanics, Vol. 61, No. 1, pp. 131–138, Mar. 1994, https://doi.org/10.1115/1.2901387
  • F. Chu and Z. Zhang, “Bifurcation and chaos in a rub-impact Jeffcott rotor system,” Journal of Sound and Vibration, Vol. 210, No. 1, pp. 1–18, Feb. 1998, https://doi.org/10.1006/jsvi.1997.1283
  • W. Qin, H. Su, and Y. Yang, “Grazing bifurcation and chaos in response of rubbing rotor,” Chaos, Solitons and Fractals, Vol. 37, No. 1, pp. 166–174, Jul. 2008, https://doi.org/10.1016/j.chaos.2006.08.018
  • H. Ma, Q. Zhao, X. Zhao, Q. Han, and B. Wen, “Dynamic characteristics analysis of a rotor-stator system under different rubbing forms,” Applied Mathematical Modelling, Vol. 39, No. 8, pp. 2392–2408, Apr. 2015, https://doi.org/10.1016/j.apm.2014.11.009
  • S. Roques, M. Legrand, P. Cartraud, C. Stoisser, and C. Pierre, “Modeling of a rotor speed transient response with radial rubbing,” Journal of Sound and Vibration, Vol. 329, No. 5, pp. 527–546, Mar. 2010, https://doi.org/10.1016/j.jsv.2009.09.016
  • J.-D. Jeng, Y. Kang, and Y.-P. Chang, “An alternative Poincaré section for high-order harmonic and chaotic responses of a rubbing rotor,” Journal of Sound and Vibration, Vol. 328, No. 1-2, pp. 191–202, Nov. 2009, https://doi.org/10.1016/j.jsv.2009.07.018
  • Y.-B. Liu et al., “Study on nonlinear Dynamics of low pressure rotor system under waviness fault,” (in Chinese), Journal of Naval Engineering University, Vol. 32, No. 1, pp. 50–56, 2020, https://doi.org/10.7495/j.issn.1009-3486.2020.01.009
  • M. Li, Y.-B. Liu, and Q. Wang, “Nonlinear dynamic Analysis of low pressure rotor-bearing system of Marine Gas Turbine,” (in Chinese), Journal of Bearing, Vol. 11, pp. 38–44, 2019, https://doi.org/10.19533/j.issn1000-3762.2019.11.010
  • H.-B. Zhang et al., “Bifurcation and Stability of Rub-impact response of Impeller rotor,” (in Chinese), Journal of Vibration Engineering, Vol. 32, No. 4, pp. 635–643, 2019, https://doi.org/10.16385/j.cnki.issn.1004-4523.2019.04.010
  • T. Sun, W.-Y. Qin, and H. Xiang, “Nonlinear response analysis of rub-impact of engine counter-rotating double-rotor system,” (in Chinese), Journal of Mechanical Strength, Vol. 32, No. 4, pp. 635–643, 2019, https://doi.org/10.16579/j.issn.1001.9669.2019.02.001
  • Zhang Ya et al., “Study on rub-impact dynamic characteristics of steam turbine rotor system considering flow-induced vibration,” (in Chinese), Journal of Beijing University of Chemical Technology (Natural Science Edition), Vol. 46, No. 1, pp. 76–83, 2019, https://doi.org/10.13543/j.bhxbzr.2019.01.012
  • Z.-Q. Lu, Y.-H. Wei, and J.-J. Jiao, “Dynamic simulation and experimental verification of rotor rub-impact fault,” Journal of Shenyang University of Technology, Vol. 37, No. 4, pp. 62–67, 2018.
  • Y.-Z. Jin, D.-Y. Wang, and B.-C. Wen, “Simulation study on rub-impact fault of aeroengine double-rotor system,” Journal of Harbin Engineering University, Vol. 38, No. 12, pp. 1872–1876, 2017.
  • M. Sayed, A. A. Mousa, and I. Mustafa, “Stability and bifurcation analysis of a buckled beam via active control,” Applied Mathematical Modelling, Vol. 82, pp. 649–665, Jun. 2020, https://doi.org/10.1016/j.apm.2020.01.074
  • M. Sayed and Y. S. Hamed, “Stability analysis and response of nonlinear rotor-seal system,” Journal of Vibroengineering, Vol. 16, No. 8, pp. 4152–4170, Dec. 2014.

Cited by

Nonlinear Dynamic Analysis on Dual‐Rotor‐Bearing‐Casing System for Marine Gas Turbine
Zhikai Xing | Qiang Wang | Yongbao Liu | Mo Li | Xin Zhang | Giacomo Innocenti
(2024)

About this article

Received
08 April 2022
Accepted
05 September 2022
Published
17 October 2022
SUBJECTS
Mechanical vibrations and applications
Keywords
double rotor
non-linear characteristic
rubbing vibration
chaos
Acknowledgements

This work is supported by National Major Science and Technology Projects of China No. J2019-Ⅳ-0021-0089 and No. J2019-Ⅰ-0012-0012, and Hubei Provincial Central Guidance Local Science and Technology, No. 2020CFB536.

Data Availability

The datasets generated during and/or analyzed during the current study are available from the corresponding author on reasonable request.

Conflict of interest

The authors declare that they have no conflict of interest.