Research Article  Open Access
Zhipeng Wang, Yunbo Yuan, Zhiyong Wang, Wei Liu, Yibin Guo, Donghua Wang, "LateralTorsional Coupling Characteristics of a TwoStage Planetary Gear Rotor System", Shock and Vibration, vol. 2018, Article ID 4293475, 15 pages, 2018. https://doi.org/10.1155/2018/4293475
LateralTorsional Coupling Characteristics of a TwoStage Planetary Gear Rotor System
Abstract
Planetary gears are one part of the whole transmission chain, and the dynamics and vibration characteristics of them are strongly coupled with external rotors. In this paper, to demonstrate the interaction between multistage planetary gears and external rotors as well as investigate the lateraltorsional coupling characteristics of them, a coupling model of a twostage planetary gear rotor system is proposed. In such a model, the twostage planetary gear subsystem is established as a lumpedparameter model and the external rotor subsystem is established as a finite element model. The vibration mode distribution properties and lateraltorsional coupling characteristics are both analyzed by modal strain energy. Three different conditions are considered: uncoupled, partially coupled, and fully coupled. The results indicate that the coupling among multiple subsystems and the lateraltorsional coupling mainly exist in the lowmode region. Natural frequencies dominated by the twostage planetary gear subsystem are sensitive to coupled conditions, whereas natural frequencies dominated by the input rotor subsystem are remarkably insensitive to coupled conditions. Furthermore, the natural frequency of the first torsional mode can be obtained only in the fully coupled condition. Experiments are implemented to obtain natural frequencies, and the experiment results validate the numerical results.
1. Introduction
In the fields of helicopter, marine, automotive, and wind power, variablespeed devices with high power density, heavy carrying capacity, high torquetoweight ratio, and efficiency are in demand. Compared with parallelaxis gears and singlestage planetary gears, multistage planetary gears are capable of transmitting the same amount of power with smaller volume and lighter weight, which makes them an important way of compact transmission [1]. Despite their benefits, multistage planetary gears often have more complex noise and vibration problems than singlestage planetary gears [2].
The mathematical model and vibration mode properties of simple, singlestage planetary gears have been studied extensively. Kahraman [3] derived a 2D planar lumpedparameter dynamic model, and he studied the vibration mode properties. Lin and Parker [4] analyzed the free vibration properties under the coordinates relative to the rotating carrier basis. Three types of modes were defined: translational, rotational, and planet modes. Cooley and Parker [5] showed that the singlemode vibration of spur planetary gear has welldefined geometry and frequency content due to its vibration mode structure. Lin and Parker [6] investigated the vibration of planetary gears with unequally spaced planet configurations. The planet modes were proved to be independent of planet spacing form, while the translational and rotational modes were considerably coupled and lost their distinct structure.
For multistage planetary gears, Kahraman [2] studied the vibration characteristics with different configurations and speed ratios using a purely rotational model. Guo and Parker [7] derived a purely rotational model of compound planetary gears to investigate the rotational vibration mode properties. They showed that there are exactly two vibration mode types, called overall modes and planet modes. Zhang et al. [8] presented a translationalrotational dynamic model of a twostage closedform planetary gear to investigate the vibration characteristics. It was found that the couplingtwist stiffness only had a significant impact on translational modes while the couplingtranslational stiffness only affected rotational modes.
However, most of the relative published papers are about planetary gears themselves and few researches on planetary gears and their external rotors have been found. Vibration mode properties of singlestage planetary gears and their external rotors were investigated by August and Kasuba [9] and Saada and Velex [10]. But there are very few literatures handling vibration mode properties of multistage planetary gears and their external rotors. Additionally, it is essential to develop a coupling mathematical model when the rigidity of external rotors is so large that input and output units have significant influence on multistage planetary gears. Through this way, one can get a more correct and more reliable vibration characteristics of multistage planetary gears.
Accordingly, the first objective of the present study is to develop a coupling mathematical model of multistage planetary gears and their external rotors; such a model takes translational and rotational degrees of freedom into account simultaneously. Then the coupling interaction between multistage planetary gear subsystem and external rotor subsystem will be investigated. A convenient and easy method to demonstrate the lateraltorsional coupling characteristics of vibration modes will be provided, since it is not easy to make mode shapes of complex systems readable.
In this paper, a lateraltorsional coupling mathematical model is proposed to investigate the free vibration characteristics of a twostage planetary gear rotor system (the simplest form of the multistage planetary gear rotor system); three different coupling conditions between each two adjacent subsystems are considered in such a model. However, the same methodology can be applied to other forms of multistage planetary gear rotor system, such as threestage and fourstage planetary gear rotor system. Furthermore, the modal strain energy (MSE) method is used to investigate the lateraltorsional coupling characteristics of vibration modes, and the relative deviation of natural frequencies dominated by the twostage planetary gear subsystem. Finally, spinning test and modal test are implemented on a test rig, and Harmonic Response Analysis and Experimental Modal Analysis are performed in order to compare the experimental natural frequencies with the numerical results.
2. Coupled Modelling and Equations of Motion
In this paper, the research subject is shown in Figure 1, a twostage planetary gear rotor system. The gear rotor system consists of a drive motor, a planetary gearbox, a load motor, input and output shafts, three couplings, two tilting pad bearings, and a balancing disk. The planetary gearbox is a twostage reducer.
As shown in Figure 2, the twostage planetary gear is connected in series by two simple, singlestage planetary gears with standstill ring gears. Its input unit is the sun gear of the firststage, and its output unit is the carrier of the secondstage. The carrier of the firststage and the sun gear of the secondstage are connected together by an involute profile inner spline.
2.1. SingleStage Planetary Gear Dynamic Model
The 2D planar lumpedparameter dynamic model of simple, singlestage planetary gears is represented in Figure 3. Two translational and one rotational degrees are attributed to every central component. In total, the model has ( is the number of planets) degrees of freedom. The absolute coordinate is established in the center of the sun gear, and the rotating coordinates are fixed and rotating with the carrier. The axis is directed toward the equilibrium position of planet 1. Translational coordinates and fixed to planet are radial and tangential coordinates, respectively.
Differing from the literature [4], the rotational coordinates are defined to be rather than , for the sake of being more convenient to the following study. Hence, the governing equations for simple, singlestage planetary gears in matrix form arewhere , , , , are the mass, gyroscopic, elastic support stiffness, tooth mesh stiffness, and centrifugal stiffness matrices, respectively. is the vector of applied external forces and torques, and represents the static transmission error excitation vector.
For further simplicity, it defines that , . Therefore, (1a) can be rewritten as follows:
2.2. Coupled Modelling of TwoStage Planetary Gear
For the twostage planetary gear in this paper, the carrier of firststage is connected with the sun gear of secondstage by an involute profile inner spline. The boundary conditions of these two components differ from those of simple, singlestage planetary gears, which is shown in Figure 4.
Define the radial coupling stiffness as , and the torsional coupling stiffness as ; then equations of motion of the carrier of firststage and the sun gear of secondstage can be derived as (3) and (4), respectively.
Extracting the coupling elements from (3) and (4), the coupling stiffness submatrix for the firststage and the secondstage planetary gears is given as follows:
Hence, the equations of motion of the twostage planetary gear can be obtained aswhere is the vector of generalized displacements, and is the vector of applied forces and torques. The assembly of the mass, gyroscopic, and stiffness matrices is illustrated in Figure 5.
(a)
(b)
(c)
2.3. Finite Element Model of External Rotor
The external rotor subsystem is analyzed by employing Timoshenko beam theory. A generalized displacement vector of each node can be defined with respect to the global coordinate system as , where , represent translational degrees, , represent bending degrees, and represents rotational degree. For the sake of simplicity, all support bearings are supposed to be linear. Only the direct stiffness is considered and the damping effect is ignored.
Figure 6(a) shows a schematic of the input rotor subsystem and Figure 6(c) shows a schematic of the output rotor subsystem. The corresponding finite element models of these two subsystems are illustrated in Figures 6(b) and 6(d). For the input rotor subsystem, there are 39 nodes, 35 shaft elements, 6 disk elements, 5 bearing elements, and 3 coupling elements. Similarly, the output rotor subsystem is divided into 3 nodes, including 1 shaft element, 3 disk elements, 1 bearing element, and 2 coupling elements, which has a total of 15 degrees of freedom.
(a)
(b)
(c)
(d)
Equations of motion of the input rotor subsystem and the output rotor subsystem in matrix form are as follows:where and are the mass matrices, and are the gyroscopic matrices, and and are the stiffness matrices. The details of these matrices are given in [11].
2.4. Coupled Modelling of Planetary Gear and External Rotor
With the matrices of the twostage planetary gear subsystem, the input rotor subsystem, and the output rotor subsystem available, the governing equations of the twostage planetary gear rotor system are constructed systematically as
The generalized displacement vector, the mass matrix, the gyroscopic matrix, and the applied force vector are given in (8b)–(8e), respectively.where subscripts , , and denote the input rotor, twostage planetary gear, and output rotor subsystem, respectively.
The stiffness matrix in (8a), the overall stiffness matrix, is demonstrated in Figure 7. It is assembled by the input rotor block submatrix, the planetary gear block submatrix, the output rotor block submatrix, and three coupling submatrices.
3. Natural Frequencies and Vibration Modes
3.1. Eigenvalue Problem and System Parameters
To determine the natural frequencies, the linear, timeinvariant representation is considered. All mesh stiffness is considered to be constant and equal to the corresponding average stiffness, which can be calculated as in [12]. Assuming the system operates at lower speed, the gyroscopic and centrifugal terms can be neglected [4]. For such a special case, (8a) can be rewritten as
The associated eigenvalue problem derived from iswhere and are the natural frequencies and the corresponding vibration modes, respectively.
For the numerical study, essential parameters of the twostage planetary gear subsystem are given in Tables 1(a)1(b). The stiffness of tilting pad bearings and rigid couplings of the external rotor subsystem are listed in Table 2.
(a)  
 
(b)  


With the parameters given in Tables 1 and 2, the natural frequencies of this twostage planetary gear rotor system can be determined, and the first 72 natural frequencies are listed in Table 3.

3.2. Modal Strain Energy of Vibration Modes
In virtue of the complexity of the system studied in this paper, vibration modes cannot be expressed in figure intuitively. For the sake of convenience, MSE method is used to weigh which subsystem is the most dominant one in the entire system for each vibration mode. MSE of the entire system is given aswhere are MSE of the input rotor, firststage planetary gear, secondstage planetary gear, and output rotor subsystem, respectively.
Furthermore, in order to investigate the lateraltorsional coupling characteristics of each vibration mode clearly, the other form of is expressed aswhere , , and are MSE associated with the lateral, torsional, and lateraltorsional coupling vibration, respectively.
With definitions of MSE in (11) and (12), MSE of each subsystem are separately calculated; meanwhile, the MSE associated with the lateral, torsional, and lateraltorsional coupling modes are determined. MSE distribution of each mode associated with each subsystem is depicted in Figure 8. And MSE distribution of each mode associated with the lateral, torsional, and lateraltorsional coupling modes is illustrated in Figure 9.
(a)
(b)
(c)
(d)
(a)
(b)
(c)
(d)
In this paper, for the convenience of the following discussion, we arbitrarily define that the region containing the first onequarter modes of a system is the lowmode region and the region containing the last threequarters modes of a system is the highmode region. Accordingly, the first 61st modes are in the lowmode region and those modes higher than 61st mode are in the highmode region, since there are 246 modes in all of this twostage planetary gear rotor system.
As shown in Figure 8, for the modes such as 2nd, , and 27th, MSE is distributed among multiple subsystems. For the modes such as 1st, , and 25th, MSE is distributed mainly in two adjacent subsystems. However, for the modes such as 3rd, 4th, 11th, and , MSE is only distributed in the input rotor subsystem and in the output rotor subsystem. Moreover, one can particularly notice that MSE mainly concentrates in one subsystem for the 61st mode as well as for the modes higher than 61st mode. The above phenomena indicate that the coupling modes among subsystems only exist in the lowmode region, whereas, in the highmode region, the modes are uncoupling and independent modes of the external rotor subsystem.
As shown in Figure 9, for the modes such as 1st, 2nd, , and , MSE is stored in three components including the lateral vibration component, the torsional vibration component, and the lateraltorsional coupling component. For the modes such as 3rd, 4th, 11th, 12th, and 23rd, MSE is only stored in the lateral vibration component; it indicates that these modes are the lateral vibration modes. Conversely, for the modes such as 26th, 47th, and 60th, MSE is only stored in the torsional vibration component; it indicates that these modes are the torsional vibration modes. However, one can particularly notice that MSE is only stored in the lateral vibration component or in the torsional vibration component for the 58th mode as well as for the modes higher than 58th mode, with no coupling components existing. It can be concluded from Figure 9 that there is no coupling effect between the lateral and torsional vibration in the highmode region.
Meanwhile, Figures 8 and 9 clearly indicate that those modes dominated by the external rotor subsystem in the highmode region are particularly the lateral and torsional decoupling modes. A main reason for this phenomenon is that no interaction between lateral and torsional vibration is taken into consideration for the finite element model of the external rotor subsystem.
Figures 8 and 9 also show that the coupling modes among subsystems in the lowmode region are particularly the lateraltorsional coupling modes, such as the and the modes. The mesh interactions of the sunplanet and the ringplanet gear pair contribute to the lateral vibration coupling with the torsional vibration. The interactions here are the torsional torque term contained in the lateral vibration equations, and, on the other hand, the torsional vibration equations contain the lateral force term as well.
3.3. Relative Deviation of Natural Frequency between Coupled and Uncoupled Conditions
In this section, the relative deviation of the natural frequencies dominated by the planetary gear subsystem and the input rotor subsystem is analyzed, respectively. Three different conditions are considered: (i) no coupling of each two adjacent subsystems is considered (uncoupled), (ii) only the coupling between the firststage and the secondstage planetary gear is considered, without interaction between the twostage planetary gear subsystem and the external rotor subsystem (partially coupled), and (iii) the coupling of each two adjacent subsystems is taken into account simultaneously, namely, the interaction between the input rotor subsystem and the firststage planetary gear and the firststage and the secondstage planetary gear, as well as the secondstage planetary gear and the output rotor subsystem being all considered simultaneously (fully coupled).
3.3.1. Natural Frequency Dominated by TwoStage Planetary Gear
Herein, to select the natural frequencies dominated by the twostage planetary gear subsystem under the above three conditions, MSE defined in (11) and (12) is separately calculated for each natural frequency. Then, the relative deviation of such natural frequencies is analyzed under partially coupled condition and under fully coupled condition relative to uncoupled condition, respectively.
The relative deviation of the natural frequencies dominated by the twostage planetary gear subsystem is shown in Figure 10. However, two rigidbody modes for the firststage and the secondstage planetary gear under uncoupled condition are not included in Figure 10. The relative deviation of such two modes will be discussed later.
(a)
(b)
As indicated in Figure 10(a), for the 9th, 10th, and modes, the natural frequencies obtained under fully coupled condition have a significant increase compared with those obtained under uncoupled condition. However, these natural frequencies almost keep constant under partially coupled condition compared to those under uncoupled condition. Compared with uncoupled condition, the 16th, 17th, 20th, 21st, 27th, and 28th natural frequencies have an increase of about 5% in partially coupled condition, and a more remarkable increase is found in fully coupled condition. However, for the 34th, 35th, and modes, the relative deviation is almost the same either under fully coupled condition or under partially coupled condition compared to uncoupled condition.
Combining the analysis in Section 3.2, one can observe that these modes, such as the 34th, 35th, and modes, are particularly dominated by secondstage planetary gear. It clearly indicates that the coupling (interaction) between firststage planetary gear and secondstage planetary gear has more significant influence on the natural frequencies dominated by secondstage planetary gear compared with the coupling (interaction) between secondstage planetary gear and the output rotor subsystem. Moreover, it can be concluded from Figures 8 and 10(a) that the natural frequencies dominated by firststage planetary gear are mainly affected by the input rotor subsystem rather than secondstage planetary gear.
As for the torsional vibration modes, Figure 10(b) shows that coupling has more impact on the relative deviation of natural frequency with the order of mode increases. Moreover, Figure 10(b) also illustrates that the natural frequencies dominated by firststage planetary gear are mainly affected by the input rotor subsystem (see the 36th mode) while the natural frequencies dominated by secondstage planetary gear are mainly affected by firststage planetary gear (see the 22nd and 33rd modes).
For the two torsional modes of the rigidbody at 0 Hz in uncoupled condition, one of them jumps to a torsional mode of the twostage planetary gear subsystem at 263.50 Hz under partially coupled condition. Furthermore, this mode jumps to a torsional mode for the entire gear rotor system at 56.60 Hz in fully coupled condition, which is tabulated in Table 4.

Under partially coupled condition, since the interaction between firststage planetary gear and secondstage planetary gear is considered, the twostage planetary gear set is considered to be a whole system. For such a system, there is only one rigidbody torsional mode, while the other one turns into a torsional vibration mode of the twostage planetary gear set at 263.50 Hz. Moreover, the interaction between each two adjacent subsystems is considered in fully coupled condition, which makes the first torsional frequency of the entire gear rotor system further decrease (see the last column of Table 4). One potential cause for this phenomenon is that the whole gear rotor system is a slender flexible rotor system.
3.3.2. Natural Frequency Dominated by Input Rotor Subsystem
Similar to the previous case, the natural frequencies dominated by the input rotor subsystem under uncoupled, partially coupled, and fully coupled condition are selected, respectively. Moreover, the relative deviation of such natural frequencies is determined in a similar way as that used in the previous case.
Excluding the torsional rigidbody mode at 0 Hz, the relative deviation of the natural frequency of the first 72 modes is illustrated in Figure 11. Figure 11(a) clearly indicates that the formulations of partially coupled analysis and fully coupled analysis contribute minimal impact on the lateral frequencies dominated by the input rotor subsystem, with the biggest relative variation ratio about 1%. Figure 11(b) illustrates that, for the torsional frequencies dominated by the input rotor subsystem, the fully coupled analysis has a more significant impact on the relative deviation compared to the partially coupled analysis, with a 4.2% reduction of the 6th torsional frequency and a 1.2% increase of the 40th torsional frequency.
(a)
(b)
What is more, as shown in Figure 11(b), the relative deviation of the natural frequencies of the 5th, 6th, and 40th modes is more significant than that of the 26th, 47th, 61st, and 63rd modes. Combining the analysis in Section 3.2, one can observe that though the former modes are dominated by the input rotor subsystem, they are also coupling modes of the input rotor subsystem and firststage planetary gear. However, the 26th, 47th, 61st, and 63rd modes are all particularly local modes of the input rotor subsystem, which indicates that the restriction effect is faint between the input rotor subsystem and firststage planetary gear for these modes.
For the torsional rigidbody modal frequency of the input rotor subsystem, it remains constant under partially coupled condition. But it jumps to 56.60 Hz under fully coupled condition, turning into the first torsional natural frequency.
4. Experimental Validation
The laboratory test rig is essentially a twostage planetary gear rotor system shown in Figure 12. Its main components have been described in Section 2. To verify the lateraltorsional coupling model proposed in this paper which considers the coupling interaction between each two adjacent subsystems, both the torsional and the lateral natural frequencies are determined on the test rig.
Complex systems, like the one considered in this study, are generally characterized by a high density of resonances and mode shapes in a relatively narrow frequency range [13]. In such an instance, it is not easy to correctly associate the numerically computed natural frequencies and mode shapes to the measured ones. To this purpose, the use of Experimental Modal Analysis or other proper strategies may be required [14]. In this present study, Harmonic Response Analysis is performed in order to compare the experimental torsional natural frequencies to the numerical ones, and Experimental Modal Analysis is performed in order to compare the experimental lateral natural frequencies to the numerical ones.
4.1. Torsional Natural Frequency Test and Validation
To gain torsional natural frequencies, the torsional vibration response is measured. The schematic of the entire experimental system is shown in Figure 13, and the experimental photo is illustrated in Figure 14. A magnetoelectric sensor is mounted in the vicinity of toothdisk coupling to collect pulse signals, and its probe is along the axis in a radial direction. Spinning Test is used in this experiment. The vibration response is measured every 50 rpm from 400 rpm to 3400 rpm under idling condition. Vibration signals are processed by a LMS Frontend (LMS SCADAS SCM05) controlled by software LMS Test.Lab to obtain the torsional frequencies.
Harmonic Response Analysis is performed in order to obtain the experimental torsional natural frequencies. The curve of the first twelve harmonics of the torsional vibration response can be obtained which is shown in Figure 15. Torsional vibration amplitude varying with the frequency is plotted, and resonant peaks appear at the natural frequencies regardless of the harmonic number. Figure 15 shows four natural frequencies of 41.7 Hz, 55.8 Hz, 145.8 Hz, and 259.6 Hz. However, according to the author’s further study, the natural frequency of 41.7 Hz is associated with a mode of the base of the test rig.
(a)
(b)
(c)
Table 5 compares the torsional natural frequencies identified in spinning test with the numerical prediction. The results show that the torsional natural frequencies determined by numerical prediction under fully coupled condition correlate with the experimental ones well, and all errors are less than 5%. Although the errors of the natural frequency of the 2nd and 3rd torsional modes are less than 5% under uncoupled and partially coupled conditions, the natural frequency of the first torsional mode is missing under these two conditions.

4.2. Lateral Natural Frequency Test and Validation
To acquire lateral natural frequencies, modal test is carried out. The schematic of the modal test is shown in Figure 16, and the picture of the test rig is illustrated in Figure 17. There are totally 18 measurement points along the input shaft. Such a modal test is implemented by the “roving hammer” approach with two accelerometers (type: Brüel & Kjær 4526; sensitivity: 9.888 mV/ms^{−2} and 10.09 mV/ms^{−2}) located at point 7. The impact force is imparted through an impact hammer (type: SINOCERA CLYD303; sensitivity: 4.06 pC/N) with both the nylon and the rubber hammer tip. Signals are sampled by a LMS Frontend (LMS SCADAS SCM05); such a Frontend is controlled by software LMS Test.Lab in a laptop. The LMS Test.Lab is adopted for signal processing, data management, and parameter estimation purposes.
Experimental Modal Analysis is performed in order to acquire the experimental lateral natural frequencies. Figure 18 shows the direction sum FRFs and the direction sum FRFs of the input rotor system. As shown in Figure 18(a), only five natural frequencies are found below 1024 Hz of direction; however, the similarity of the two corresponding mode shapes of and is higher than 70%, which indicates that these two modes are the same mode. As shown in Figure 18(b), six potential natural frequencies are found below 1024 Hz of direction; however, the similarities of mode shapes of and as well as and are both higher than 70%. Thus, the real lateral natural frequencies are 82.9 Hz, 352.5 Hz, 454.2 Hz, and 802.7 Hz of direction and 83.5 Hz 368.7 Hz, 435.5 Hz, and 794.6 Hz of direction.
(a)
(b)
Table 6 compares the lateral natural frequencies identified in impact test with the numerical prediction. The results indicate that the numerical lateral natural frequencies under three conditions agree with the experimental ones well, and most errors are less than 5%, except for the 18th and 19th mode. Furthermore, there is a marginal difference of these eight lateral natural frequencies among these three conditions; in particular this phenomenon confirms the correctness of the conclusions drawn in Section 3.3.2.

5. Conclusion
In this paper, the influence of external rotors has been investigated on the vibration characteristics of planetary gears and the lateraltorsional coupling characteristics of the entire system, and a coupled modelling methodology of a twostage planetary gear rotor system is presented. The proposed fully coupled model offers significant advantages in identifying natural frequencies and vibration mode properties. The conclusions can be summarized as follows:(1)A fully coupled model for the twostage planetary gear rotor system is presented, where the twostage planetary gear subsystem is modelled by lumpedparameter method and the external rotor subsystem by finite element method.(2)The natural frequencies dominated by the twostage planetary gear subsystem are sensitive to coupling effect. However, the coupling does not obviously affect the natural frequencies dominated by the input rotor subsystem. The first order of the torsional natural frequency is absent unless the fully coupled condition (the interaction between each two adjacent subsystems) is considered.(3)The coupled modes among multiple subsystems or each two adjacent subsystems mainly appear in the lowmode region; however, in the highmode region, there are particularly uncoupled and independent modes of the external rotor subsystem. Furthermore, in the highmode region, there is no coupling effect between the lateral and torsional vibration.(4)Spinning test and modal test have been done to obtain lateral and torsional natural frequencies, and the experiment data validate the numerical results.
Data Availability
All the data supporting the conclusions of the study are available on readers’ request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The research work is supported by the Heilongjiang Province Funds for Distinguished Young Scientists (Grant no. JC 201405), China Postdoctoral Science Foundation (Grant no. 2015M581433), and Postdoctoral Science Foundation of Heilongjiang Province (Grant no. LBHZ15038).
References
 P. Lynwander, Gear Drive Systems: Design and Application, CRC Press, 1983.
 A. Kahraman, “Free torsional vibration characteristics of compound planetary gear sets,” Mechanism and Machine Theory, vol. 36, no. 8, pp. 953–971, 2001. View at: Publisher Site  Google Scholar
 A. Kahraman, “Load sharing characteristics of planetary transmissions,” Mechanism and Machine Theory, vol. 29, no. 8, pp. 1151–1165, 1994. View at: Publisher Site  Google Scholar
 J. Lin and R. G. Parker, “Analytical characterization of the unique properties of planetary gear free vibration,” Journal of Vibration and Acoustics, vol. 121, no. 3, pp. 316–321, 1999. View at: Publisher Site  Google Scholar
 C. G. Cooley and R. G. Parker, “The geometry and frequency content of planetary gear singlemode vibration,” Mechanical Systems and Signal Processing, vol. 40, no. 1, pp. 91–104, 2013. View at: Publisher Site  Google Scholar
 J. Lin and R. G. Parker, “Structured vibration characteristics of planetary gears with unequally spaced planets,” Journal of Sound and Vibration, vol. 233, no. 5, pp. 921–928, 2000. View at: Publisher Site  Google Scholar
 Y. Guo and R. G. Parker, “Purely rotational model and vibration modes of compound planetary gears,” Mechanism and Machine Theory, vol. 45, no. 3, pp. 365–377, 2010. View at: Publisher Site  Google Scholar
 L. Zhang, Y. Wang, K. Wu, R. Sheng, and Q. Huang, “Dynamic modeling and vibration characteristics of a twostage closedform planetary gear train,” Mechanism and Machine Theory, vol. 97, pp. 12–28, 2016. View at: Publisher Site  Google Scholar
 R. August and R. Kasuba, “Torsional vibrations and dynamic loads in a basic planetary gear system,” Journal of Vibration, Acoustics, Stress, and Reliability in Design, vol. 108, no. 3, pp. 348–353, 1986. View at: Publisher Site  Google Scholar
 A. Saada and P. Velex, “An extended model for the analysis of the dynamic behavior of planetary trains,” Journal of Mechanical Design, vol. 117, pp. 241–247, 1995. View at: Publisher Site  Google Scholar
 J. Rao S, T. Shiau N, and J. Ghang R, “Theoretical analysis of lateral response due to torsional of geared rotors,” Mechanism and Machine Theory, vol. 33, no. 6, pp. 761–783, 1998. View at: Google Scholar
 X. H. Liang, M. J. Zuo, and T. H. Patel, “Evaluating timevarying mesh stiffness of a planetary gear set using potential energy method,” Journal of Mechanical Engineering Science, vol. 228, no. 3, pp. 535–547, 2014. View at: Publisher Site  Google Scholar
 T. M. Ericson and R. G. Parker, “Natural frequency clusters in planetary gear vibration,” Journal of Vibration and Acoustics, vol. 135, no. 6, Article ID 061002, 2013. View at: Publisher Site  Google Scholar
 T. M. Ericson and R. G. Parker, “Planetary gear modal vibration experiments and correlation against lumpedparameter and finite element models,” Journal of Sound and Vibration, vol. 332, no. 9, pp. 2350–2375, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Zhipeng Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.