Published: 31 October 2023

Rigid-flexible coupling modeling of the dual-rotor system for aero-engine

Pingping Ma1
Yuehui Dong2
Hangqi Zhao3
Yang Li4
Muge Liu5
1School of Machinery and Automation, Weifang University, Weifang, 26100, China
2Dongfeng Machine Tool Plant Co., Ltd., Shiyan, 442022, China
3State units, Shenyang, 110000, China
4Air Force Engineering University, Xian, 261000, China
5School of Mechanical Engineering, Liaoning University of Technology, Jinzhou, 121000, China
Corresponding Author:
Pingping Ma
Views 66
Reads 26
Downloads 142

Abstract

The dual-rotor system is the core component of advanced aero-engine. Establishing a reasonable, accurate, and efficient dynamics model is the key to studying the dynamics and vibration of the rotor system for aero-engine. This manuscript takes a representative aircraft engine dual-rotor system as a prototype, considers the rig-id-flexible coupling characteristics of different stiffness elastic supports and rotor structures, and establishes an analytical dynamic model of the dual-rotor system. Based on the established dynamic model, the natural char-acteristics of the dual-rotor system are analyzed. The model was validated using two different research methods: the rigid-flexible coupling multi-body system dynamics simulation platform ADAMS, and finite element analy-sis. The dynamic model of the dual-rotor system established in this paper can meet the requirements of hierar-chical rigid-flexible coupling of system and structure, overall mass distribution, and stiffness distribution. In particular, it can also effectively realize the simulation of multi-facet and multi-phase unbalanced vibration of the rotor system. The research methods of this paper can further enrich the basic theory of dynamics and vibration of the aero-engine rotor system.

Rigid-flexible coupling modeling of the dual-rotor system for aero-engine

Highlights

  • The dual-rotor system is the core component of advanced aero-engine
  • Accurate dynamic model is the key to correctly predict rotor system vibration
  • Elastic strain energy will be excessively introduced for parts with high stiffness under elastic assumption
  • Rigid-flexible coupling model can more accurately express the dynamic characteristics of aero-engine dual-rotor system

1. Introduction

The dual-rotor system of aero-engine has multiple elastic supports, complex and special structures, and many factors are affecting its dynamic characteristics. Exploring a simplified model of the dual-rotor system that can reflect the structural characteristics and dynamic behavior of the aero-engine is the key to establishing an accurate and efficient dynamic model of rotor system for aero-engine. The model of the double rotors and multi-supports system used in the current research mainly includes the following types: the dual-rotor system with one disc for the high-pressure rotor and one disc for the low-pressure rotor [1-4], the dual-rotor system with two discs for the high-pressure rotor and two discs for the low-pressure rotor [5-7], and the dual-rotor system with multiple-discs [4, 8-10]. The rigid-flexible coupling of rotor system structure is not considered in dynamic modeling. When analyzing parts with large stiffness, the elastic strain energy will be too much included in the analysis based on the elastic hypothesis. The calculation error of the dynamics analysis and vibration prediction for the dual-rotor system is too large, so it is difficult to realize the accurate analysis and vibration prediction of the dynamic characteristics of an aero-engine rotor.

Aiming at the theoretical analysis requirements of dynamics and vibration characteristics of the dual-rotor system for aero-engine, considering multiple supports and the coupling of internal and external rotor structures, the energy method is used to establish an analytical dynamic model. The model has four lumped mass discs (low-pressure compressor, low-pressure turbine, high-pressure compressor, and high-pressure turbine) and five supports with different stiffness characteristics, and the low-pressure turbine shaft is considered as a flexible shaft. Based on the analytical model, the inherent characteristics of the dual-rotor system are analyzed and compared with the finite element results and multi-body dynamics simulation results, which verifies the effectiveness and accuracy of the analytical model. The accurate modeling of the dual-rotor system for aero-engine is realized, which lays a foundation for the subsequent vibration characteristics analysis of the dual-rotor system.

2. Dynamic modeling of the dual-rotor system

2.1. Mechanical model of the dual-rotor system

This section takes a typical aviation engine dual-rotor system as a prototype and establishes a simplified mechanical model of the dual-rotor system based on the principle of dynamic similarity. Through the qualitative study of the simplified model, this paper explores the important parameters that affect the regularity of the motion behavior of the system. According to the structure and working characteristics of the dual rotor system, a simplified model that can reflect the main dynamic behavior of the dual-rotor system for aero-engine is established, as shown in Fig. 1. The modeling process follows the dynamic similarity principle, and the obtained model has the same mass ratio as the prototype rotor system.

Fig. 1Simplified mechanical model of the dual-rotor system

Simplified mechanical model of the dual-rotor system

In the simplified model of the dual-rotor system, the compressor and turbine disc are simplified to a disc with equivalent mass and moment of inertia. The actual multi-disc structure is simplified into a few disc structures. The specific simplification method is as follows: the four-stage fan discs in the low-pressure rotor are simplified as disc D1 (LPC), the low-pressure turbine disc is equivalent to disc D2 (LPT), the nine-stage compressor discs of the high-pressure rotor is simplified as disc D3 (HPC), and the high-pressure turbine disc is equivalent to disc 4 (HPT).

The supports of the dual-rotor are simplified into five springs with different stiffness. The support stiffness Kb1, Kb3 and Kb5 are the combined equivalent stiffness of the bearing, squirrel cage and support at the front fulcrum of the fan rotor, the front fulcrum of the high-pressure rotor and the rear fulcrum of the low-pressure turbine respectively. Kb2 is the combined stiffness of support 2 and support 3. The intermediate fulcrum adopts the same roller bearing as the prototype, and the stiffness is Kb4.

2.2. Rigid-flexible coupling dynamic model of the dual-rotor system

Taking the simplified model of the dual-rotor system shown in Fig. 1 as the research object, the motion analytical equation of the dual-rotor system is deduced with the help of the Lagrange method. The support is considered to be elastic and simulated with equivalent stiffness. The rotating speed of the rotating shaft is considered to be constant, and the disc, fan shaft, and high-pressure shaft are considered to be rigid. In particular, the elastic characteristics of the slender shaft of the low-pressure turbine are considered.

2.2.1. Kinetic energy of the dual-rotor system

The kinetic energy of the dual-rotor system includes the translational kinetic energy and rotational kinetic energy of the four discs, which can be recorded as follows:

1
T=Tt1+Tr1+Tt2+Tr2+Tt3+Tr3+Tt4+Tr4,

where, Ttiand Tri (i=1~4) are the translational kinetic energy and rotational kinetic energy of disc respectively.

Based on the method of establishing a coordinate system and determining the generalized displacement and velocity of the rotor system, the kinetic energy of the dual-rotor system is obtained. The mass, polar moment of inertia, and diameter moment of inertia of the rotor system are represented by mi, Jpi, and Jdi (i=1:4) respectively. The translational and rotational kinetic energy of turntables 1-4 are:

2
Tti=12mi(x˙pi2+y˙pi2+z˙pi2) =12mix˙i2+y˙i2+z˙i2+ei2Ωi2-2y˙ieiΩisin(Ψi+Ψi0)+2z˙ieiΩicos(Ψi+Ψi0),
3
Tri=12ΩiTJiΩi =12Ω1+θ˙ziθyi-θ˙zisinΨi+θ˙yicosΨiθ˙zicosΨi+θ˙yisinΨiTJpi000Jdi000JdiΩi+θ˙ziθyi-θ˙zisinΨi+θ˙yicosΨiθ˙zicosΨi+θ˙yisinΨi =12JpiΩi2+2JpiΩiθ˙yiθzi+Jdiθ˙yi2+Jdiθ˙zi2.

2.2.2. Potential energy of the dual-rotor system

The potential energy of the dual-rotor system is shown in Eq. (4), including the elastic potential energy of five fulcrums, the potential energy of the coupling and the elastic potential energy of the low-pressure turbine shaft:

4
U=U1+U2+U3+U4+Uc+U5+Us.

It includes the elastic potential energy of five supports and the potential energy of coupling:

5
Ui=12qBiTKBiqB i, i=1~4,

where, qBi (i=1~4) is the displacement vector of the rotor pivot Bi (i=1~4) in the generalized coordinate system. KBi is the stiffness matrix of support Bi (i=1~4)

The elastic potential energy of the intermediate bearing is as follows:

6
U5=12qCH-qCLTKB5qCH-qCL.

Among them, qCH and qCL are the generalized displacement vectors of the intermediate support points on the high-pressure and low-pressure rotors, respectively. KB5 is the stiffness matrix of the intermediate support.

The elastic potential energy at the coupling is as follows:

7
UC=12qC2-qC1TKCqC2-qC1,

where, qC1 and qC2 are the generalized displacement vectors at coupling connection points C1 and C2, respectively. KC is the equivalent stiffness of the coupling.

Potential energy of the elastic shaft of a low-pressure turbine:

8
Us=12FyFzMyMzTδ1100δ140δ22δ2300δ32δ330δ4100δ44FyFzMyMz.

Among them, Fy and Fz are the forces acting on the low-pressure turbine shaft in the y and z directions. My and Mz are the torque of the low-pressure turbine shaft in the y and z directions. δij (i=1~4, j=1~4) is the flexibility coefficient.

2.2.3. Dynamic equation of the dual-rotor system

The dynamic model of the dual-rotor system is established through the Lagrange equation:

9
m1y¨1-kb1yθz1a1+kb1yy1+12EcIcLc3θz12+θz22Lc+θz1c1+θz2c2+y1-y2=0,
m1z¨1+kb1zθy1a1+kb1zz1+-12EcIcLc3θy12+θy22Lc+θy1c1+θy2c2-z1+z2=0,
Jd1θ¨y1+Jp1Ω1θ˙z1+4EcIcθy1+θy2/2Lc+kb1zθy1a12+kb1zz1a1 +12EcIycθy1c1+θy22c1+12θy2c2-z12+z22Lc+c1θy1c1+θy2c2-z1+z2Lc3=0,
Jd1θ¨z1-Jp1Ω1θ˙y1+4θz1+θz2/2EcIcLc+kb1yθz1a12-kb1yy1a1 +12EcIycθz1c1+θz2c12+12θz2c2+y12-y22Lc+c1θz1c1+θz2c2+y1-y2Lc3=0,
m2y¨2+-kb2ya2-kb4ya41+kb5ya5θz2+y2-a34θz4-y4kb4y -6EcIcθz1+θz2Lc+12EcIcθz1c1+θz2c2+y1-y2Lc3+kb5y+kb2yy2 +3IEsLla23a2-Ll3-y2(3a22-3a2Ll+Ll2)+θz2a2(2a22-3a2Ll+Ll2)=0,
m2z¨2+kb2za2+kb4za41-kb5za5θy2+z2+a34θy4-z4kb4z +6EcIcθy1+θy2Lc+12EcIcθy1c1+θy2c2-z1+z2Lc3+kb5z+kb2zz2 +3IEsLla23a2-Ll3-z2(3a22-3a2Ll+Ll2)+θy2a2(2a22-3a2Ll+Ll2)=0,
Jd2θ¨y2+Jp2Ω1θ˙z2+kb2za22+kb4za412+kb5za52+kb4θyθy2
+z2+a34θy4-z4kb4za41+kb2za2-kb5za5z2-kb4θyθy4
+12Ecc2Iycθy1c1+θy2c2-z1+z2+6EcIcθy1c1+θy1c2+2θy2c2-z1+z2Lc+2EcIcθy1+2θy2Lc2Lc3
+3IEsLlz2(2a2-Ll)a22a2-Ll2-3IEsLlθy2a2a2-Ll=0,
Jd2θ¨z2-Jp2Ω1θ˙y2+a22kb2y+a52kb5y+a412kb4y+kb4θzθz2
+a41-y2+a34θz4+y4kb4y+-a2kb2y+a5kb5yy2-kb4θzθz4
+12 EcIcc2c1θz1+c2θz2+y1-y2+2 EcIcθz1+2 θz2Lc2+3 c1θz1+3 c2θz1+6 c2θz2+3 y1-3 y2LcLc3
+3IEsLly22a22a2-Ll22a2-Ll+3IEsLly22a2-Lla22a2-Ll2-3IEsLlθz2a2a2-Ll=0,
m3y¨3 -a3θz3kb3y+y3kb3y+6 EgIgθz3+θz4Lh+2 y3-2 y4Lh3=0,
m3z¨3 +a3θy3kb3z+z3kb3z-6EgIg θy3+θy4Lh-2 z3+2 z4Lh3=0,
Jd3θ¨y3+Jp3Ω3θ˙z3 +a32kb3zθy3+a3kb3zz3+4EgIgθy3+12θy4Lh-32(z3-z4)Lh2=0,
Jd3θ¨z3 - Jp3Ω3θ˙ y3 +a32kb3yθz3-a3kb3yy3+4EgIg θz3+12 θz4Lh+32 y3-y4Lh2=0,
m4y¨4+(y4+a34θz4+θz2a41-y2)kb4y-6 EgIgθz3+θz4Lh+2 y3-2 y4Lh3=0,
m4z¨4-a34θy4+θy2a41+z2-z4kb4z+6 EgIgθy3+θy4Lh-2 z3+2 z4Lh3=0,
Jd4θ¨y4+Jp4Ω3θ˙z4+a34a34θy4+θy2a41+z2-z4kb4z -kb4θyθy2+kb4θyθy4+2 EgIgθy3+2 θy4Lh-3 z3+3 z4Lh2=0,
Jd4θ¨z4-Jp4Ω3θ˙y4+a34(y4+a34θz4+θz2a41-y2)kb4y -kb4θzθz2+kb4θzθz4+2 EgIgθz3+2 θz4Lh+3 y3-3 y4Lh2=0,

where, ai (i=1~5) and ci (i=12) represents the distance between two points in Fig. 1, in units of m. a1=B1D1, a2=B2D2, a3=B3D3, a4=B4D4, a5=B5D2, c1=C1D1, c2=B5C2. Additionally, a41 is the distance between two points B4andD2. ah is the distance between two points B4 and D4, expressed as ah=a4-a34. Lc is the length of the coupling, expressed as Lc=C1C2. ll is the total length of the low-pressure turbine shaft, expressed as ll=B2B5. Esand Ec is the elastic modulus, in units of kg/(m·s2). Is and Ic is the moment of inertia, in units of kg·m·s2.

3. The inherent characteristics of the dual-rotor system

Mastering the natural frequency and vibration mode is the basis of studying the dynamic response of the dual-rotor system. The inherent characteristics of the dual-rotor system are obtained based on the above dynamic Eq. (9).

The first three natural frequencies and vibration modes of the simplified model are shown in Fig. 2. Campbell diagram is shown in Fig. 3, and the first two critical speeds are shown in Table 1.

Fig. 2Analytical results for the first three orders natural frequencies and vibration modes of the dual-rotor system

Analytical results for the first three orders natural frequencies  and vibration modes of the dual-rotor system

a) The first order (58.8 Hz)

Analytical results for the first three orders natural frequencies  and vibration modes of the dual-rotor system

b) The second order (75.4 Hz)

Analytical results for the first three orders natural frequencies  and vibration modes of the dual-rotor system

c) The third order (127 Hz)

Table 1The first two critical speeds of the dual-rotor system

Order
The first order
The second order
Value of the critical speeds
3528 r/min
4524 r/min

It can be seen from Fig. 2 that the first vibration mode is represented by the bending of the low-pressure turbine rotor and the translation of the high-pressure rotor. The second vibration mode is characterized by the low-pressure turbine rotor bending and the high-pressure rotor pitching. The third-order modes are low-pressure fan rotor pitch, low-pressure turbine rotor bending, and high-pressure rotor pitch.

Fig. 3Campbell diagram for the dual-rotor system

Campbell diagram for the dual-rotor system

4. Model validation

4.1. Model verification of the dual-rotor system based on multi-body dynamics

In this section, the analytical dynamic model of the dual-rotor system obtained in Section 2 is verified by the multi-body dynamics simulation analysis method. The rigid-flexible coupling simulation model of the dual-rotor system is carried out based on the software ADAMS, and the natural characteristics are analyzed. The results obtained by simulation are compared with the analytical results.

Based on the theory of multi-body dynamics, ADAMS can combine the mechanical boundary, load boundary, and geometric conditions of complex rotor systems for aero-engine in a concise simulation analysis model. ADAMS can realize the rigid-flexible coupling modeling of the aero-engine rotor system and can simulate the system motion and dynamic characteristics under actual working conditions. The rigid-flexible coupling simulation model of the aero-engine dual-rotor system is shown in Fig. 4.

Fig. 4Rigid-flexible coupling dynamic model of the dual-rotor system based on ADAMS

Rigid-flexible coupling dynamic model of the dual-rotor system based on ADAMS

Where the low-pressure turbine shaft is set as a flexible body, and other components adopt the rigid assumption. The supports are defined by stiffness and damping, simulated by the "bushing" element. The coupling is simplified to a short shaft structure. The rotation of the dual-rotor system is realized by adding a drive to the center of the rotating shaft of the low-pressure rotor and the high-pressure rotor respectively.

The inherent characteristics of the dual-rotor system are calculated by the simulation method. The first three natural frequencies and vibration modes of the dual-rotor system are shown in Fig. 5.

Fig. 5The first three orders natural frequencies and vibration modes of the dual-rotor system based on ADAMS

The first three orders natural frequencies and vibration modes  of the dual-rotor system based on ADAMS

a) The first order (56.4 Hz)

The first three orders natural frequencies and vibration modes  of the dual-rotor system based on ADAMS

b) The second order (79 Hz)

The first three orders natural frequencies and vibration modes  of the dual-rotor system based on ADAMS

c) The third order (117.1 Hz)

The first vibration mode is represented by the bending of the low-pressure turbine rotor and the translation of the high-pressure rotor; The vibration mode is characterized by the low-pressure turbine rotor bending and the high-pressure rotor pitching. The third-order modes are low-pressure fan rotor pitch, low-pressure turbine rotor bending, and high-pressure rotor pitch. The first three modes are consistent with the results obtained in the previous Section 3 based on the analytical dynamic model.

4.2. Model verification of the dual-rotor system based on finite element method

In this section, the analytical dynamic model of the dual-rotor system obtained in section 2 is verified by the finite element method. The finite element model of the dual-rotor system is carried out, and the natural characteristics are analyzed. The results obtained by simulation are compared with the analytical results. In addition, the results of the dual-rotor system with four discs and five supports are compared with the prototype model.

Fig. 6Finite element model of the dual-rotor system with four discs and five supports

Finite element model of the dual-rotor system with four discs and five supports

4.2.1. Finite element model of the dual-rotor system with four discs and five supports

During the establishment of the finite element model of the dual-rotor system, the following simplifications are made: (1) The shaft section is considered to be elastic. The mass and moment of the inertia of the shaft section are placed on the nodes at both ends. (2) The supports are considered as spring elements without cross stiffness and damping. (3) The fan disc, low-pressure turbine disc, high-pressure compressor disc, and high-pressure turbine disc are simplified as hollow discs, which are regarded as rigid bodies with rotation effects. It is concentrated at the center of the mass of each disc in the form of concentrated mass and moment of inertia.

The finite element model of the dual-rotor system is established, as shown in Fig. 6. The low-pressure rotor of the dual-rotor system with four discs and five supports is discretized with 30 nodes and 29 shaft segments. The high-pressure rotor system is discretized with 12 nodes and 11 shaft segments.

The first three natural frequencies and vibration modes of the low-pressure rotor system are calculated as shown in Fig. 7. The first vibration mode is the bending of the low-pressure turbine rotor. The second vibration mode is the bending of the fan rotor and the low-pressure turbine rotor. The third vibration mode is the bending of the low-pressure turbine section.

Fig. 7The first three orders natural frequencies and modes of the low-pressure rotor obtained by finite element method

The first three orders natural frequencies and modes  of the low-pressure rotor obtained by finite element method

a) The first order (76.9 Hz)

The first three orders natural frequencies and modes  of the low-pressure rotor obtained by finite element method

b) The second order (122.6 Hz)

The first three orders natural frequencies and modes  of the low-pressure rotor obtained by finite element method

c) The third order (150.6 Hz)

Fig. 8The first three orders natural frequencies and modes of the high-pressure rotor obtained by finite element method

The first three orders natural frequencies and modes  of the high-pressure rotor obtained by finite element method

a) The first order (60.4 Hz)

The first three orders natural frequencies and modes  of the high-pressure rotor obtained by finite element method

b) The second order (96.4 Hz)

The first three orders natural frequencies and modes  of the high-pressure rotor obtained by finite element method

c) The third order (388.2 Hz)

The first three natural frequencies and vibration modes of the high-pressure rotor are calculated as shown in Fig. 8. The first vibration mode is translational of the high-pressure rotor, the second vibration mode is the pitching of the high-pressure rotor, and the third vibration mode is the bending of the high-pressure rotor.

The first three natural frequencies and vibration modes of the dual-rotor system with four discs and five supports are obtained, as shown in Fig. 9. The first vibration mode is the bending of the low-pressure turbine rotor and the translational of the high-pressure rotor. The second vibration mode is the bending of the low-pressure turbine rotor and the pitching of the high-pressure rotor. The third vibration mode is the bending of the fan rotor and the low-pressure turbine rotor, the pitching of the high-pressure rotor.

Fig. 9The first three orders natural frequencies and modes of the dual-rotor system obtained by finite element method

The first three orders natural frequencies and modes  of the dual-rotor system obtained by finite element method

a) The first order (55.9 Hz)

The first three orders natural frequencies and modes  of the dual-rotor system obtained by finite element method

b) The second order (75.5 Hz)

The first three orders natural frequencies and modes  of the dual-rotor system obtained by finite element method

c) The third order (122.0 Hz)

4.2.2. Finite element model of the dual-rotor system for the prototype

The finite element model of the dual-rotor system of the prototype is established, as shown in Fig. 10. The low-pressure rotor system of the prototype is discretized with 54 nodes and 53 shaft segments. The high-pressure rotor is discretized with 47 nodes and 46 shaft segments.

Fig. 10Finite element model for the dual-rotor system of the prototype

Finite element model for the dual-rotor system of the prototype

Fig. 11The first three orders natural frequencies and modes for the low-pressure rotor of the prototype

The first three orders natural frequencies and modes for the low-pressure rotor of the prototype

a) The first order (89.6 Hz)

The first three orders natural frequencies and modes for the low-pressure rotor of the prototype

b) The second order (122.3 Hz)

The first three orders natural frequencies and modes for the low-pressure rotor of the prototype

c) The third order (156.8 Hz)

The first three natural frequencies and vibration modes of the low-pressure rotor system of the prototype are obtained as shown in Fig. 11. The first vibration mode of the low-pressure rotor of the prototype is the bending of the low-pressure turbine rotor. The second vibration mode is the pitching of the fan rotor. The third vibration mode is the pitching of the fan rotor and the bending of the low-pressure turbine rotor.

The first three natural frequencies and vibration modes of the high-pressure rotor of the prototype are shown in Fig. 12. The first vibration mode is translational of the high-pressure rotor, the second vibration mode is the pitching of the high-pressure rotor, and the third vibration mode is the bending of the high-pressure rotor.

Fig. 12The first three orders natural frequencies and modes for the high-pressure rotor of the prototype

The first three orders natural frequencies and modes for the high-pressure rotor of the prototype

a) The first order (59.8 Hz)

The first three orders natural frequencies and modes for the high-pressure rotor of the prototype

b) The second order (106.4 Hz)

The first three orders natural frequencies and modes for the high-pressure rotor of the prototype

c) The third order (587.7 Hz)

The first three natural frequencies and vibration modes of the dual-rotor system for the prototype are obtained as shown in Fig. 13. The first vibration mode is the bending of the low-pressure turbine rotor and the translational of the high-pressure rotor. The second vibration mode is the bending of the low-pressure turbine rotor and the pitching of the high-pressure rotor; The third vibration mode is the pitching of the fan rotor, bending of the low-pressure turbine rotor, and translation of the high-pressure rotor.

4.2.3. Comparative analysis of results

The first three natural frequencies and vibration modes obtained from the two finite element models of the dual-rotor system are compared, and the following conclusions are obtained: 1) Among the first three vibration modes of the low-pressure rotor system, the low-pressure turbine rotor are the same, and the fan rotor shows bending and pitching in the simplified and prototype model respectively. The relative errors of the first three natural frequencies are 14 %, 0.2 %, and 3.9 % respectively. 2) The first three modes of the high-pressure rotor system are consistent, and the relative errors of the first three natural frequencies are 1 %, 9 %, and 33 % respectively. 3) Among the first three vibration modes of the dual-rotor system, the vibration modes of the low-pressure turbine rotor and the high-pressure rotor are the same. The fan rotor shows the bending vibration mode in the simplified model and the pitching vibration mode in the prototype model. The relative errors of the first three natural frequencies are 0.3 %, 5.2 %, and 8.6 % respectively.

The first three vibration modes of the dual-rotor system obtained in Section 2 based on the analytical dynamic model consist with the results of the prototype model. It further shows that it is reasonable to consider the fan rotor as rigid and the low-pressure turbine shaft as flexible in the rigid-flexible dynamic model of the dual-rotor system.

Fig. 13The first three orders natural frequencies and modes for the dual-rotor system of the prototype

The first three orders natural frequencies and modes for the dual-rotor system of the prototype

a) The first order (55.7 Hz)

The first three orders natural frequencies and modes for the dual-rotor system of the prototype

b) The second order (79.7 Hz)

The first three orders natural frequencies and modes for the dual-rotor system of the prototype

c) The third order (133.5 Hz)

5. Conclusions

In response to the urgent demand for a rotor system model that can reasonably and accurately express the dynamics and vibration characteristics of the rotor system in the study of vibration in aircraft engine rotor systems, this article takes a representative aero-engine dual-rotor system as the research object, established a simplified mechanical model which can reflect the structural characteristics and motion behavior of the prototype model is introduced. Considering the rigid-flexible coupling characteristics of an elastic support with different stiffness and rotor structure, and analytical dynamic model of the dual-rotor system is established. The model contains the kinetic energy of four equivalent discs and the elastic potential energy of five supports and a low-pressure turbine shaft, with 16 degrees of freedom.

The dynamic model of the rigid-flexible coupled dual-rotor system is verified and confirmed by the rigid-flexible coupling multi-body dynamics simulation method and finite element method. The feasibility of using the analytical dynamic model of the dual-rotor system to study the vibration is clarified, which lays a model foundation for the subsequent research on the vibration of the dual-rotor system. It shows that it is reasonable to consider the fan rotor as rigid and the low-pressure turbine shaft as flexible in the rigid-flexible dynamic model of the dual-rotor system.

References

  • Xun Hu, Guihuo Luo, and Deping Gao, “Numerical analysis and experiment of counter-rotating dual-rotor’s steady-state response,” (in Chinese), Journal of Aeronautical Power, Vol. 22, No. 7, pp. 1044–1049, 2007, https://doi.org/10.3969/j.issn.1000-8055.2007.07.004
  • J. Liu, C. Wang, and Z. Luo, “Research nonlinear vibrations of a dual-rotor system with nonlinear restoring forces,” Journal of the Brazilian Society of Mechanical Sciences and Engineering, Vol. 42, No. 9, pp. 1–20, Sep. 2020, https://doi.org/10.1007/s40430-020-02541-w
  • Peng Gao, Lei Hou, and Yushu Chen, “Nonlinear vibration characteristics of a dual-rotor system with inter-shaft bearing,” Shock and Vibration, Vol. 38, No. 15, pp. 1–10, 2019.
  • Guihuo Luo, Hailun Zhou, Fei Wang, and Xiguan Yang, “Dynamic response of co-rotating and counter-rotating dual-rotor system supported on ball bearing,” Journal of Aeronautical Power, Vol. 27, No. 8, pp. 1887–1894, 2012.
  • Guo Chen, “Vibration modeling and analysis for dual-rotor aero-engine,” (in Chinese), Journal of Vibration Engineering, Vol. 24, No. 6, pp. 619–632, 2011.
  • Z. Lu, X. Wang, L. Hou, Y. Chen, and X. Liu, “Nonlinear response analysis for an aero engine dual-rotor system coupled by the inter-shaft bearing,” Archive of Applied Mechanics, Vol. 89, No. 7, pp. 1275–1288, 2019.
  • Xiguan Yang, Guihuo Luo, Zhenyuan Tang, and Fei Wang, “Modeling method and dynamic characteristics of high-dimensional counter-rotating dual rotor system,” Journal of Aeronautical Power, Vol. 29, No. 3, pp. 585–595, 2014.
  • Z.-J. Feng and G.-H. Luo, “Analysis on the dynamic characteristics of the dual-rotor structures of a certain aero-engine,” Vibroengineering Procedia, Vol. 5, pp. 469–474, Sep. 2015.
  • Dayi Zhang, Yehui Liu, Zhichao Liang, Yanhong Ma, and Jie Hong, “Prediction for critical speed of double spools system in aero engines,” Journal of Propulsion Technology, Vol. 36, No. 2, pp. 292–298, 2015.
  • N. Wang, H. Xu, and D. Jiang, “Dynamic model and fault feature research of dual-rotor system with bearing pedestal looseness,” Mathematical Problems in Engineering, Vol. 2016, pp. 1–18, 2016, https://doi.org/10.1155/2016/3817405

About this article

Received
13 July 2023
Accepted
15 September 2023
Published
31 October 2023
Keywords
dual-rotor system
dynamic model of the dual-rotor system
rigid-flexible coupling dynamic simulation model
Acknowledgements

Study on Dynamic Optimization Method of Large Squeeze Film Damper Driven by Data (12072069).

Data Availability

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

Author Contributions

Pingping Ma is responsible for providing main ideas and writing articles. Yuehui Dong provides information collection and data analysis. Muge Liu is responsible for revising the format and grammar of the paper. Qingkai Han is responsible for reviewing the content and results.

Conflict of interest

The authors declare that they have no conflict of interest.