Spectrum Analysis of a Coaxial Dual-Rotor System with Coupling Misalignment

/e finite element model of a dual-rotor system was established by Timoshenko beam element. /e dual-rotor system is a coaxial rotor whose supporting structure is similar to that of an aero-engine rotor system. /e inner rotor is supported by three bearings, which makes it a redundantly supported rotor. /e outer rotor connects the inner rotor by an intershaft bearing. /e spectrum characteristics of the dual-rotor system under unbalanced excitation and misalignment excitation were analysed in order to study the influence of coupling misalignment of the inner rotor on the spectral characteristics of the rotor system. /e results indicate that the vibration caused by the misaligned coupling of the inner rotor will be transmitted to the outer rotor through the intershaft bearing. Multiple harmonic frequency components, mainly 1x and 2x, will be excited by the coupling misalignment. /e amplitudes of all harmonic frequencies increase with the misalignment in both the inner and outer rotors. /e vibration level of the outer rotor affected by the misalignment is lower than that of the inner rotor because it is far from the misaligned coupling. Harmonic resonance occurs when any harmonic frequencies of the misalignment response coincide with a natural frequency of the system. In order to verify the theoretical model, experiments are performed on a test rig. Both the experimental and simulation results are in good accordance with each other.


Introduction
Rotor misalignment is one of the most common difficulties in the operation of rotating machinery. A misalignment rotor negatively influences the rolling, sealing, and coupling parts and can also produce eccentricity in the air gap. erefore, a misaligned machine is more prone to failure due to the increased loads on bearings and couplings. In general, the misalignment in a rotor system is unavoidable and cannot be completely eliminated. Furthermore, it can disguise itself very well on industrial rotating machinery, which makes it not easy to be detected. Even if the alignment looks good when you do an offline check, a running misalignment may occur. What we witness is the secondary effects of misalignment as it slowly damages the machinery over long periods of time. Hence, the ability to clearly diagnose the presence of misalignment using vibration analysis can be vital in reducing costly machine unscheduled downtime. Vibration based identification of faults, such as rolling bearing fault [1], gear fault [2], crack [3], and rub [4], is welldeveloped and widely used in practice. However, because more attention is paid to the study of the misalignment of the series connection rotor (such as steam turbine rotor) than that of a coaxial rotor, the misalignment of a coaxial dual-rotor system remains as an outstanding area, where the basic understanding is somewhat lacking.
In the fault diagnosis of misalignment, Lees et al. [5] proposed a method that can reliably estimate both the rotor unbalance and misalignment from measured vibration during a single machine run-down. Lal and Tiwari [6] developed a model based approach which can estimate the bearing and coupling dynamic parameters of a laboratory test rig as a measure of misalignment. Rybczynski [7] indicated the potential of using trajectory patterns for diagnosing misalignment defects in rotating machines and suggested including sets of trajectory patterns to the knowledge base of a machine diagnostic system. Verucchi et al. [8]studied the techniques that can detect misalignment from the analysis of the motor electrical variable, and the results show that the conclusions drawn for a particular coupling are not necessarily applicable to others. Reddy and Sekhar [9] studied a technique that uses torque measurements to detect and monitor coupling misalignment. Qi et al. [10] proposed a novel method to analyse in time frequency domain under different misalignment conditions, and the experimental results show that the proposed method is able to detect misalignment fault characteristic frequency effectively. Patel et al. [11][12][13] investigated the influence of misalignment and its type on the forcing characteristics of flexible coupling by the full spectrum.
In terms of the dynamic characteristics of rotor misalignment, Hili et al. [14] and Fakhfakh et al. [15] developed a model with the shaft angular misalignment to identify misalignment from the vibration response. Li et al. [16] and Li et al. [17] discussed the nonlinear dynamic behaviours of a rotor system on misaligned journal bearings. Askarian and Hashemi [18] obtained the effects of unbalance and misalignment on the vibration of a multirotor system. Li et al. [19] presented a theoretical research on the periodic response of an offset disc rotor system with coupling misalignment. Ma et al. [20] discussed the stability of a rotorbearing-coupling system under misalignment and found that the instability rotating speed range increases with misalignment. AL-Hussain [21] indicated that an increase in angular misalignment or mechanical coupling stiffness terms leads to an increase of the model stability region. Yi et al. [22] present a simplified model of misaligned installing matched bearings and found that the amount of misalignment plays an important factor on both the global stability and vibration amplitude for a misaligned matched bearing rotor system. Zhang et al. [23] presented that the support force of a rotor system under misalignment faults appears double vibration and high harmonics frequencies. Nikolajsen [24] studied the relationship between radial bearing misalignment and stability of a multibearing rotor and found that neither perfect alignment nor equal load sharing between bearings can be relied upon to give optimum stability. Gibbons [25] depicted general formulas to calculate the misalignment reaction forces for both gear couplings and flexure couplings in 1976. Sekhar and Prabhu [26] derived the reaction force formulas for the angular misalignment of flexure couplings based on Gibbons formulas. Lee and Lee [27] developed a dynamic model for a misaligned rotor-ball bearing system by the Gibbons-Sekhar formula. Both the experimental and simulation results showed that, as the angular misalignment increases, the whirling orbits tend to collapse toward a straight line and the natural frequency of the misaligned rotor system associated with the misalignment direction increases largely. e reaction forces calculated from the Gibbons-Sekhar formula are the static load. For a rotating shaft, the static load can be decomposed as series dynamic misalignment forces. Prabhakar et al. [28], Jalan et al. [29], Wan et al. [30], and Ma [31] et al. built several rotor bearing system models using the dynamic misalignment forces. All these studies showed that 2x frequency components occur no matter for angular misalignment or parallel misalignment. Dewelland Mitchell [32] also gave a discussion of the harmonics arising inflexible couplings. Lees [33] illustrated the source of 2x harmonic excitation in misaligned rotors with a rigid coupling uses a purely linear model. Xuand Marangoni [34,35] developed a model of a motor-coupling-rotor system. e torque due to the misalignment of a universal joint was derived in their pair of papers, and their theoretical and experimental studies indicate that the forcing frequencies caused by shaft misalignment are even multiple frequencies of the motor rotation speed. e research object of the previous studies is the multispan rotor system. However, little information has been done on the misalignment of the coaxial rotor. A coaxial rotor of an aero-engine is different from a series connection rotor of steam turbine in twofold. For one thing, the coaxial shafts corotate or count-rotate at different speeds, for another, the inner rotor which is supported by more than two bearings is a redundant structure. ese special structures deserve to be paid more attention. Lu et al. [36] presented a theoretical research on the nonlinear response characteristics of an aero-engine dual-rotor-bearing system with flexible coupling misalignment faults in the inner rotor. Li et al. [37] described the mechanism and influencing factors of nonlinear properties of a misaligned inner rotor in an aero-engine through Lagrange equations. Wang et al. [38] and Wang et al. [39] et al. developed a multibearing rotor model with misalignment to investigate the dynamic response when the misalignment parameters are uncertain. Li et al. [40] also investigated the quantification of uncertainty effects on the dynamic responses and vibration characteristics of a multirotor bearing system with the fault of angular misalignment. Ren et al. [41] presented a dynamic model of a flexible rotor system with multisupport in order to solve the bearing misalignment problem of the rotor system in an aero-engine. Previous studies have given detailed discussions to understand the effects of misalignment on a dualrotor system of an aero-engine. However, most of them either took the inner rotor system instead of the whole dual rotors system or showed no experimental evidence to support the theoretical research. us, more researches are still required to illustrate the effects of misalignment on a coaxial dual-rotor system.
Spectrum analysis is one of the basic techniques used to conduct misalignment diagnosis. erefore, the spectrum analysis of an aero-engine multibearing dual-rotor system with coupling misalignment and disk unbalance has been studied by the finite element method (FEM) in the present study. A test rig has a similar bearing supporting structure like a real aero-engine was built to verify the theoretical studies. Experiments were performed on the test rig under different intentional misalignment conditions. e study could be useful in the detection of misalignment in the coaxial dual-rotor system diagnosis and in the design of an aero-engine rotor system.

Modelling of a Dual-Rotor-Bearing System with Coupling Misalignment
A structural diagram of a typical twin-spool jet engine is shown in Figure 1 [42]. LP Comp, HP Comp, HP Turb, and LP Turb denote low pressure compressor, high pressure compressor, high pressure turbine, and low pressure turbine, respectively. e jet engine is a coaxial dual-rotor system. e shaft of the inner rotor is solid and the shaft of the outer rotor is hollow. e two rotors are connected by an intershaft bearing. e intershaft bearing inner race is mounted on the inner rotor while the outer race is mounted on the outer rotor. e outer and inner rotors rotate independently, which can rotate at different rotating speeds and different directions. e rotor constitutive elements are shafts (circular section) and disks. And the inner rotor is supported by three bearings, which makes it a redundantly supported rotor. A coupling is designed to accommodate the misalignment of the inner rotor. e coupling is located at the end of the first span of the inner rotor near the second bearing, which is shown in Figure 2.

Finite Element Model of the Dual-Rotor System.
e dualrotor system shown in Figure 2 is modelled by the FEM [43]. To develop the equations of motion for the rotor bearing system, a stationary coordinate system is defined, which is shown in Figure 3. e stationary coordinate system is consisting of three mutually perpendicular axes-ox, oy, and oz-intersecting at the point o and fixed in space. e axis of the dual-shaft rotor, in equilibrium, is coincident with the axis oz. In the equilibrium position, the rotor lies on a straight line passing through the bearing canters. e axes ox, oy, and oz, in that order, form a right-handed set, which is defined as follows.
A rotation of a right-handed screw from ox to oy advances it along oz which is assumed to be horizontal. A rotation of a right-handed screw from oy to oz advances it along ox, which is perpendicular to the oyz-plane. A rotation of a right-handed screw from oz to ox advances it along oy which is assumed to be vertical.
Using the FEM, the rotor is divided into a number of shaft elements with nodes at both ends of each element. e finite shaft element of the dual-rotor system is accomplished by means of the Timoshenko beam. Each shaft element has two nodes, and the node location is shown in Table 1. Disks and bearings are assumed to be attached to the shaft at these nodes where required. e inner rotor is modelled with 15 elements and 16 nodes. e outer rotor is modelled with 8 elements and 9 nodes. us, the model has 23 elements and 25 nodes in total. e nodes of the inner rotor are denoted by blue dots and those of the outer rotor are denoted by red squares, which are shown in Figure 2. Only lateral vibrations are considered, so each node has four generalized coordinates (four degrees of freedom per node): transverse displacements in the x-and y-directions and rotations about the x-and y-axes. e translations of the shaft from the equilibrium position are u and v in the x-and y-directions, respectively. e rotations are defined such that θ is a positive rotation about the x-axis and φ about the y-axis.
e rotors and disks are made from steel. e rotor structure is simplified, and the change of some shaft sections is ignored. e diameters of the shafts are shown in Table 2.
All the five disks are identical. Disk 1 and disk 2 represent a low pressure compressor. Disk 3, disk 4, and disk 5 represent a low pressure turbine, a high pressure compressor, and a high pressure turbine, respectively. Disk 1, disk 2, disk 3, disk 4, and disk 5 are located at nodes 3, 5, 13, 20, and 22, respectively. e disk properties are listed in Table 3.
ere are five bearings in the system, namely, bearing 1, bearing 2, bearing 3, bearing 4, and bearing 5. Bearing 5 is an intershaft bearing. All bearings are assumed as to be linear and isotropic. Bearing 1, bearing 2, bearing 3, and bearing 4 are located at nodes 1, 7, 15, and 18, respectively. e intershaft bearing, bearing 5, is located at node 12 and node 24. e bearings properties are listed in Table 4.
According to the Lagrange equation, the dynamic equations of disk element and shaft element can be derived as follows: where M e d and M e s are the mass matrixes of the disk element and shaft element, respectively. q e d and q e s are the displacement vector of the disk element and shaft element, respectively. Since a disk element has one node, q e d is a 4dimensional vector. A shaft element has two nodes; q e s is an 8-dimensional vector. Ω d and Ω s are the rotating speeds of the disk element and shaft element, respectively. e rotating speed of the disk element and shaft element on the outer rotor is 1.6 times of that of the disk element and shaft element on the inner rotor. G e d and G e s are the gyroscopic matrixes of the disk element and shaft element, respectively. K e s is the stiffness matrix of a shaft element. F e d and F e s are the external forces on the disk element and shaft element, respectively. e dynamic equations of intershaft bearing and an ordinary bearing are the same, which can be expressed as where C e b is the damping matrix of bearings. K e b is the stiffness matrix of bearings. q e b is the displacement vector of bearings. F e b is the external force of bearings. e difference between the ordinary bearing and intershaft bearing is that the external force on the intershaft bearing should also include the interaction forces between the two nodes connected by the intershaft bearing. e interaction forces, F intershaft and -F intershaft , are shown in Figure 4. e intershaft bearing connects node j of the outer rotor and node i of the inner rotor. e inner and outer rotors interact with each other through the intershaft bearing. e interaction forces Shock and Vibration 3 are the action force and reaction force. Although these two forces are equal in magnitude and opposite in direction, they cannot be cancelled out. ese two forces should be calculated as the external forces of the corresponding node, because node i and node j where the intershaft bearing is located are not adjacent. Considering all the damping, stiffness, and mass matrices of shaft, disk, and bearing, the general dynamic equations of the dual-rotor system can be written as where M, C, and K are the mass, damping, and stiffness matrixes of the dual-rotor system, respectively. q is the displacement vector. Ω 1 and Ω 2 are the rotation speed of the inner rotor and outer rotor, respectively. G 1 and G 2 are the gyroscopic matrixes of the inner and outer rotor, respectively. F g is the gravity vector. F c is the reacting force and moment vector caused by the misaligned coupling, which   Figure 2: Schematic of the dual-rotor system, nodes, and elements. will be derived in Section 2.2. F u is the vectors of the unbalance excitation of the rotor system. F u can be written as where F i u are the out-of-balance forces at node i. m i is the mass of the disk and the shaft attached to node i. ε i is the eccentricity distance from the centre of mass to the shaft centre line at node i.Ω i is the rotation speed of the unbalance disk.

Reacting Forces and Moments Caused by Coupling
Misalignment.
e effects of misalignment on a rotor system depend on how the reaction forces accommodate the Outer rotor Inner rotor   Shock and Vibration 5 misalignment. ere are many factors that affect the reacting forces and moments of a misaligned coupling, and the structure of the coupling is one of the important factors.
[8] Figure 5 shows one kind of coupling used in an aero-engine rotor system [37]. e coupling is used to assemble a low pressure fan rotor and a low pressure turbine rotor in an aeroengine dual-rotor system. e coupling is a spline coupling having some radial locating surfaces and a hexagon nut. e special structure stiffens the spline coupling stronger than a common gear coupling and softer than the same diameter shaft. Moreover, a splined joint instead of a fixed joint makes the spline coupling accommodate misalignment. erefore, the misalignment of assembly errors and changing of operating temperature occurs indeed in the spline coupling.
According to the structural characteristics of the spline coupling, it can be regarded as a flexure coupling. e reaction forces and moments caused by the spline coupling misalignment can be calculated by the Gibbons-Sekhar formula [25,26], which are shown in Figure 6. e formula can be given as follows.
For parallel misalignment, where the misalignment angles θ 1 , φ 1 , θ 2 , and φ 2 can be obtained by the following equations: For angular misalignment, It can be observed that the reaction forces F x2 and F y2 and the reaction moments M x2 and M y2 are the functions of the angular displacements and lateral translational displacements of misalignment, which are displayed in Figure 6. e bending spring rate (k b ) and the centre of articulation (Z 3 ) of the coupling are also displayed in the figure. e parallel and angular misalignment parameters (ΔX 1 , ΔX 2 , ΔY 1 , ΔY 2 , and θ 3 ) can be measured using the reverse indicator method. Assuming that z 1 is the driving side, (+)T q is applied as shown in Figure 6, and the rotation is in the same direction as the applied torque. [44].
As already mentioned in the introduction section, the reaction forces and moments calculated from the Gibbons-Sekhar formula are static loads. For a rotating shaft, the static loads can be decomposed as series dynamic periodic forces which cause vibrations. e dynamic misalignment forces and moments are treated as excitations at the coupling node of the finite element model, and only the 1x, 2x, 3x, and 4x harmonic components are considered in this study. Assuming that all reacting forces and moments caused by the coupling misalignment are located at node j of the inner rotor, the forces and moments are given as follows [30,31]: F j c is the dynamic reaction forces and moments of the misaligned coupling at node j. Ω j is the rotation speed of the coupling.
For the dual-rotor system shown in Figure 2, the misalignment occurs in the coupling of the inner rotor. Because the inner rotor is supported by three bearings, the supporting structure can only be an angular misalignment. Suppose that bearing 1 is offset by △ millimetre corresponding to a misalignment angle of θ 3 degrees which is shown in Figure 7. erefore, only the effects of the angular misalignment on the spectral characteristics of the dualrotor system will be studied in this paper.

Simulation and Discussions
In this section, two case studies are performed to investigate the effects of coupling misalignment on the spectral characteristics of the dual-rotor system: (1) the dual-rotor system with only unbalanced excitations; (2) the system with both unbalanced excitations and misaligned excitations. e frequency responses of the two conditions are compared with each other to understand the mechanism of misalignment in the dual-rotor system. Assume that the inner and outer rotors rotate in the reverse direction; the rotation speed of the outer rotor is 1.6 times that of the inner rotor.
Regardless of the unbalance moment, assume that only the out-of-balance forces exist at both node 3 (disk 1) of the inner rotor and node 20 (disk 4) of the outer rotor; all reacting forces and moments caused by the coupling misalignment are located at node 6. And the time domain dynamic response of the system can be obtained by solving the dynamic equation of the dual-rotor system with the Newmark-β method. e frequency response of the system can be obtained by the FFT algorithm.

Frequency Response of Out-of-Balance.
A Campbell diagram is a map of natural frequency against rotor rotation speed. It is useful in identifying the forward and backward precession natural frequency components of the rotor system during spectral analysis. e Campbell diagram of the dual-rotor system, plotted against the inner rotor speed, is shown in Figure 8. Shown on the Campbell diagram are the lines giving the correspondence between the natural frequencies plotted on y-axis and the rotation speed of the inner rotor plotted on x-axis. As the rotor speed increases, each natural frequency pair diverges due to gyroscopic effects. FW in red and BW in black denote the forward and backward natural frequencies, respectively. erefore, the first-order forward and backward natural frequencies are expressed as 1st BFW and 1st BW, and so on. e first-, second-, third-, and fourth-order vibration of the inner rotor are represented as 1x, 2x, 3x, and 4x, respectively. e firstand second-order vibration of the outer rotor are denoted as 1.6x and 3.2x, because the rotating speed of the outer rotor is 1.6 times of that of the inner rotor. e spectrum diagrams of the dual-rotor system with only unbalance forces are shown in Figures 9 and 10 when the rotation speed of the inner rotor is 1260 rpm (21 Hz) and that of the outer rotor is 2016 rpm (33.6 Hz). e figures show the following.
(1) e dominant frequencies of the dual-rotor system are 1x and 1.6x when only unbalanced excitations exist in the system. e amplitude of 1.6x is always greater than that of 1x either in the inner rotor or the outer rotor although the two rotors have the same unbalanced mass. Due to the higher speed of the outer rotor, the unbalanced force of the outer rotor is greater than that of the inner rotor. (2) e vibration amplitudes in horizontal direction and vertical direction at the same bearing (Figures 9(a) M z1 (T q )       Shock and Vibration and 9(b) and Figures 10(a) and 10(b)) are at the same level, because the stiffness is isotropic. However, it will change obviously when the misalignment reacting force is applied, which will be discussed in Section 3.2. e amplitude of 1.6x in the outer rotor is greater than that in the inner rotor. e vibration attenuated when the unbalance force of the outer rotor passes through the intershaft bearing to the inner rotor. (3) Resonance occurs at bearing 1 when 1.6x (33.6 Hz) is close to the 2nd BW natural frequency (33.18 Hz) of the rotor system, which is shown in Figure 9, so does bearing 4 of the outer rotor, when 1x (21 Hz) is close to 1st FW (24.47 Hz), which is shown in Figure 10. is is different from the Jeffcott rotor whose unbalance force can only arouse the forward precession natural frequency. is is caused by the different rotation directions of the inner and outer rotors. e unbalanced forces on the outer rotor play a dominant role, which caused the outer rotor whirl in the forward precession while the inner rotor is in the backward precession.

Frequency Response of Misalignment.
In order to analyse the influence of misalignment on the spectrum characteristics of the dual-rotor system, it is assumed that the working conditions of the rotor system are exactly the same as the unbalanced conditions (as described in Section 3.1) except for the additional misalignment forces and moments. Bearing 1 is assumed to be up 0.3 mm (Δ � 0.3 mm) in the vertical direction, shown in Figure 7, and the inner rotor produces a misalignment angle of 0.21°. e reacting forces and moments are calculated and substituted into the dynamic equations of the system. e misalignment responses of the rotor system are obtained by solving the dynamic equation. e spectra of bearing 1 and bearing 4 are shown in Figures 11 and 12. e spectra of bearing 2 and bearing 3 are not separately listed here because they are similar to that of bearing 1 and bearing 4. e spectrum of bearing 5 is not shown separately because the vibration signal of the intershaft bearing cannot be verified by a direct test measurement.
Comparing Figures 11 and 12 with Figures 9 and 10, the following is shown: (1) e vibration spectrum of misalignment is more complicated than that of unbalance. Both the inner and outer rotors appear frequency components dominated by 1x, 1.6x, and 2x. In addition, 3x, 4x, and multiple forward and backward precession natural frequencies are also excited by the misalignment. In particular, both the 3rd FW and 3rd BW natural frequencies of the rotor are excited in the vicinity of 2x frequency. (2) Misalignment intensifies the vibration in the misalignment direction. Comparing Figures 11(a) and 11(b) indicates that the vibration amplitude of bearing 1 in the vertical direction is significantly higher than that in the horizontal direction, although the stiffness of bearing 1 is isotropic. It is due to that the intentional misalignment was set in the vertical direction. e force of misalignment is directional, which does not rotate with the rotor. However, the unbalance causes a rotating force, which is a onceper-revolution fault. e direction of the unbalanced force will rotate synchronously with the rotor. us, the huge gap between the two directions will not appear in the unbalanced response. is phenomenon is not obvious on bearing 4 of the outer rotor by comparing Figures 12(a) and 12(b), because bearing 4 is far from the misaligned coupling. (3) Note that the further away from the misalignment coupling, the less affected by misalignment. Misalignment increases the vibration level of both the inner and outer rotor, but the most obvious change is in the direction of misalignment near the misaligned coupling. e maximum vibration amplitude of bearing 1 is higher than that of bearing 4 under the condition of misalignment (Figures 11 and 12), because bearing 4 is further away from the misaligned coupling. However, as shown in Figures 9  and 10 in Section 3.1, the maximum vibration amplitude of bearing 1 is lower than that of bearing 4 when only an unbalanced force excites the system. e vibration of bearing 3 also proves this phenomenon. Due to the similar situation, the spectral of bearing 3 is not listed repetitively. e radar chart shown in Figure 13 is plotted in order to further illustrate the effects of misalignment on different bearings and different harmonic frequency components. As shown in Figure 13, each spoke represents the vibration direction of a bearing. e length of a spoke represents the amplitude of a frequency component. B1H represents the horizontal direction of bearing 1 and B1V stands for the vertical direction of bearing 1, and so on. e figures show the following: (1) Although bearing 2 is closest to the misaligned coupling, the maximum amplitude of the four bearings does not occur at the misalignment direction of bearing 2. Instead, the maximum amplitude occurs in the vertical direction of bearing 1 (B1V) because the inner rotor is a redundantly support rotor. e vibration of bearing 2 is suppressed because bearing 2 and bearing 3 make the inner rotor a stable structure.
(2) Although the stiffness of the four bearings is isotropic. e vibration amplitudes in the vertical direction are greater than that of the horizontal direction in both bearing 1 and bearing 2, which are near the misaligned coupling. However, bearing 3 and bearing 4, away from the coupling, show the opposite vibration level. significantly. e frequencies of 3x and 4x change a little. But that's not always the case when a resonance occurs, which will be explained in Section 3.3. e 1.6x radar diagram is nearly circular, indicating that the first-order component of the vibration of the outer rotor (1.6x) is almost unaffected by the misalignment of the inner rotor.
e effects of misalignment on rotor frequency characteristics have been qualitatively analysed. e quantitative influence of different degrees of misalignment on the harmonic frequencies of the system will be analysed too. Assuming that bearing 1 is raised by 0.1 mm, 0.2 mm, 0.3 mm, 0.4 mm, and 0.5 mm, respectively, the variation trend of the vertical amplitude of bearing 1 with the misalignment value is obtained, which is shown in Figure 14. It can be seen that the amplitude of all harmonic components increases with the increasing of the misalignment angle, but the frequency of 1x and 2x increases more obvious than 3x and 4x.

10
Shock and Vibration 3.3. Harmonic Resonance. As described in Section 3.1, resonance occurs when the frequency of the unbalance force is close to a natural frequency. Resonance occurs in the misalignment response too. Any harmonic frequency of the misalignment response may be coincided with a natural frequency of the system, which causes resonance. e resonance of bearing 1 in the vertical direction at different rotation speeds is shown in Figures 15-17. e spectrum of bearing 1 is shown in Figure 15 when the inner rotor speeds are 11.5 Hz and 13 Hz. It can be found that the amplitude of 3x at 11.5 Hz is almost six times that of 3x at 13 Hz although the spectrum components are the same at the two speeds. At the rotation speed of 11.5 Hz, the resonance occurs because the 2nd FW (33.52 Hz) is close to 3x (34.5 Hz). However, the amplitude decreases significantly when the rotation speed is increased to 13 Hz. Similarly, the frequency 2x coincides with the 2nd FW natural frequency when the rotation speed is 17 Hz. e resonance occurs which is shown in Figure 16(a). However, there is no resonance when the rotation speed is increased to 19 Hz, which is shown in Figure 16(b). e frequency 1x coincides with the 2ndBW natural frequency when the rotation speed is 32 Hz. Resonance occurs which is shown in Figure 17(a). No resonance occurs when the speed is reduced to 30 Hz, which is shown in Figure 17(b).
From the above analysis, it can be known that all harmonic frequencies excited by misaligned coupling, 1x, 2x, 3x and 4x, may cause resonance. Even the weak frequencies, like 3x and 4x, can cause severe vibration when resonance occurs. erefore, it is necessary to consider the harmonic resonance of misaligned coupling in the design of a dualrotor system.

Orbit
Analysis. An orbit is a powerful diagnostic tool which can provide important and useful information about the dynamic motion of the rotor. e orbits of the dual-rotor system with different degree misalignment are shown in Figures 18-20. According to the previous analysis, the misalignment has the greatest impact on bearing 1 and the weakest impact on bearing 4. erefore, bearing 1 and bearing 4 are taken as the research object to analyse the influence of misalignment on orbits. e influence of misalignment on the orbit of bearing 2 and bearing 3 is between bearing 4 and bearing 1 and is no longer listed separately. e orbit plots of rotor without and with misalignment are shown in Figures 18 and 19 when the inner rotor speed is 19 Hz and the outer rotor speed is 30.4 Hz. e orbit plots of rotor with different degree misalignment are shown in Figure 20 when the inner rotor speed is 11 Hz and the outer rotor speed is 17.6 Hz. A comparative analysis of the above figures shows the following: (1) e orbits of bearing 1 and bearing 4 are both in the shape of rings when only unbalanced excitation exists, as shown in Figure 18. e orbits are multilooped rings rather than perfect circle because the two rotors of the dual-rotor system running at different speeds and the vibration frequency components are more complex than the single 1x frequency. (2) e orbits of bearing 1 tend to collapse toward a straight line for the angular misalignment along the y-direction. Nikolajsens [24] also reported the same phenomenon of misalignment. Figures 19(a) and 20(a) show that the shape of orbits will change with the rotation speed when the misalignment is at the same degree. However, Figures 18(b) and 19(b) show that the orbit of bearing 4 has little change due to the weak influence of misalignment.
(3) Figure 20 shows that the shape of orbits change with the degree of misalignment. And the orbit diverges into a reticulated diamond shape from a straight line shape as the increase of misalignment (from 0.1 mm to 0.3 mm).

Experimental Study of the Rotor System with Coupling Misalignment
A dual-rotor test rig is built in order to verify the simulation results. Experiments were performed on the test rig under different intentional misalignment conditions.

Test Rig Configuration.
e configuration of the test rig is shown in Figure 21. e test rig mainly consists of the following four parts: (1) a pair of rotors, an inner rotor and an outer rotor, and they rotate in different speed and reverse direction; (2) a pair of frequency-variable motors: the inner motor and outer motor independently drive the inner rotor and outer rotor; (3) five identical disks: disk 1, disk 2, and disk 3, are mounted on the inner rotor and disk 4 and disk 5 are mounted on the outer rotor; (4) five bearings: the inner rotor is supported by three bearings, namely, bearing 1, bearing 2, and bearing 3. e outer rotor is supported by bearing 4 and bearing 5 (intershaft bearing). e bearing types are presented in Table 5. e high supporting stiffness of the five bearings makes the critical speed of the dual-rotor system exceed the rated speeds of the two motors. erefore, five squirrel cages are used to support the bearings to reduce the supporting stiffness. Bearing 1 is padded up by inserting gaskets (shown in Figure 22) under the bearing pedestal to produce intentional misalignment.
Vibration data are acquired from the four bearing pedestals by a B & K data collector. Both horizontal and vertical vibrations of the bearing pedestals are measured by accelerometers. e acceleration signal collected by the sensor is firstly converted into a displacement signal by double integration, and then the FFT algorithm is applied on   the displacement signal to obtain the spectrum. Because the low-frequency noise will be generated in this process, a highpass filter is adopted to filter low-frequency noise signals below 10 Hz.

Experimental Results and Discussions.
ree test plans are performed to verify the simulation results. e first experiment aims to confirm the spectral components of the misaligned dual-rotor system. When the spacer thickness is 0.3 mm and the rotation speed of the inner rotor is 21 Hz, the spectra of bearing 1 and bearing 4 are shown in Figures 23 and 24. e figures indicate the following.
(1) e misalignment excited the frequencies of 1x, 2x, 3x, 4x, and 1.6x and several natural frequencies, which are consistent with the theoretical analysis as shown in Figures 11 and 12. It should be pointed out that the second-order component of vibration of the outer rotor (denoted by 3.2x) appeared in the experimental spectrum. It does not appear in simulation. It may be caused by the misalignment between the inner and outer rotor. Further study about 3.2x will be conducted in the future. ) of the same bearing, it can be seen that the vibration amplitude of the vertical direction is larger than that of the horizontal direction, although the actual stiffness of bearing pedestal in the vertical direction is greater than that in the horizontal direction. is phenomenon confirms the simulation results that misalignment will intensify the vibration in the misalignment direction. e second experiment is conducted to check the effects of different degrees of misalignment on the amplitude of different frequency components. Gaskets with thickness of 0.2 mm, 0.3 mm, 0.4 mm, 0.5 mm, and 0.6 mm, shown in Figure 22, are customized to produce different misalignment angles. Because the 0.6 mm gasket is too thick, it is easy to cause bearing 1 and bearing 2 to heat and smoke. erefore, the 0.6 mm gasket was not used for the experiment. e 0.1 mm thick gasket is too thin to be processed. erefore, only 0.2 mm, 0.3 mm, 0.4 mm, and 0.5 mm gaskets are used for the experiments.
e spectra corresponding to different spacer thicknesses are tested at the same speed (Ω 1 � 21 Hz, Ω 2 � 33.6 Hz). Two sets of tests are performed on each gasket, recording the amplitude of each frequency component. en the amplitude of each frequency is averaged, and the experimental results are finally obtained as shown in     Figure 14 that the amplitude of each harmonic component increases with the degree of misalignment. is result is not completely consistent with Dewell's [32] research results, which showed that 1x does not change with the degree of misalignment. is may be due to the different structure of the coupling. Instead of a spline coupling, a disk coupling is the research object in his study. e third experiment is designed to verify the harmonic resonance caused by misalignment. e vibration of bearing 1 and bearing 4 at two rotation speeds is measured, as shown in Figures 26 and 27. It can be clearly seen from the two figures that the harmonic resonance occurs when the rotation speed of the inner rotor is 18 Hz, because 3x is close to 2nd FW. is is not the case at 11 Hz. e test is consistent with the theoretical analysis too.
Due to the limitation of experimental conditions, the orbits of the rotor cannot be measured directly. e orbits of the rotor are replaced by the orbits of bearings to analyse the influence of misalignment. e orbits of bearings are plotted by the vibration displacement of bearing pedestals in horizontal and vertical directions. Figure 28 shows the orbits without misalignment when the speeds of the inner rotor and outer rotor are 1300 rpm (21.6 Hz) and 2080 rpm (34.6 Hz). e orbits with misalignment (Δ � 0.2 mm) are shown in Figure 29. By comparing the two figures, it can be seen that the orbit of bearing 1 tend to collapse toward a straight line because of the reacting forces and moments caused by the misalignment. Figure 30 shows that the orbits diverge toward a shape of reticulated diamond with the severity of misalignment (from 0.2 mm to 0.4 mm). e experimental results confirm the simulation results very well.

Conclusions
In this paper, a finite element model for a dual-rotor system with a misaligned coupling is developed using the Timoshenko beam elements. e response of the system is generated using the Newmark-β method. e spectral response of the system with both unbalances and misalignment has been studied. e experiments are performed to illustrate the theoretically developed characteristic vibration frequencies for a misaligned coupling in the dual-rotor system. Experimental and theoretical studies show a good agreement with each other, which verifies the analytical procedure adopted in this study. Some conclusions drawn from the study can be summarized as follows: (1) e inner rotor with coupling misalignment makes both the inner and outer rotor appear 1x, 2x, 3x, and 4x harmonic components, among which 1x and 2x are the main frequency components.
(2) Although the vibration caused by misalignment of the coupling in the inner rotor can be transmitted to the outer rotor, the closer the coupling is to the bearing, the more severe the vibration is. And the vibration amplitude of the misalignment direction is larger than that of the direction without misalignment. (3) e amplitudes of each harmonic frequency increase with the increase of the misalignment. (4) Harmonic resonance occurs when any harmonic frequencies of the misalignment response coincide with a natural frequency of the system. And the vibration will be intensified. (5) e orbits of bearings tend to collapse toward a straight line for the angular misalignment. And the orbit diverges into a reticulated diamond shape from a straight line shape as the increase of misalignment.
Data Availability e data supporting this research article are available from the corresponding author or first author on reasonable request.

Conflicts of Interest
e authors declare that there are no conflicts of interest regarding the publication of this paper.