Abstract

The dynamic investigation of helical planetary gears plays an important role in structure design as the vibration and noise are perceived negatively to the transmission quality. With consideration of the axial deformations of members, the gyroscopic effects, the time-variant meshing stiffness, and the coupling amongst stages, a three-dimensional dynamic model of the two-stage helical planetary gears is established by using of the lumped-parameter method in this paper. The model is applicable to variant number of planets in two stages, different planet phasing, and spacing configurations. Numerical simulation is conducted to detect the structured vibration modes of the equally spaced systems. Furthermore, the unique properties of these vibration modes are mathematically proved. Results show that the vibration modes of the two-stage helical planetary gears can be categorized as five classes: the rigid body mode, the axial translational-rotational mode, the radical translational mode, and the 1st-stage and the 2nd-stage planet mode.

1. Introduction

Planetary gears are widely used in aerospace and wind turbine industries because of their high reliability and compactness. As the higher load carrying capacity and higher reduction ratio compared with spur planetary gear sets, the helical planetary gear sets are the norm for all automotive applications [1]. The high-frequency dynamic meshing forces generated by the operated planetary gears are typically larger than the quasi-static forces. Besides, vibration and noise induced by the dynamic forces are potentially hazard for the reliability and durability of the gears and bearings. Hence, the dynamic modeling and analysis of planetary system have great significance for designers.

Many researches are conducted on the dynamics of planetary gear to reduce the vibration and noise of the system. Kahraman [2] extended a three-dimensional dynamic model of a single-stage planetary gear train to obtain the natural modes. Lin and Parker [35] established a translational-rotational dynamic model of a single-stage spur planetary gear and investigated the free vibration characteristics of systems with equally spaced and diametrically opposed planets, respectively. Eritenel and Parker [6] formulated the three-dimensional model of the single-stage helical planetary gears and mathematically categorized the structured modal properties. Sheng et al. [7] investigated the vibration modal properties of the double-helical planetary gear system with three-dimensional motion using numerical and analytical approach. Qian et al. [8] identified the vibration structures of a particular planetary gear used in a coal shearer and validated the model by comparing with the finite element model. Abousleiman and Velex [9] proposed an original hybrid FE/lumped-parameter model of planetary gear sets to simulate the three-dimensional dynamic behavior of planetary gears with tooth modifications and errors. Zhang et al. [10] investigated the vibration characteristics of the closed-form spur planetary gears. Most above dynamic models belong to lumped-parameter models [17, 1019] and finite element models [8, 9, 17, 2028]. The long computation times of the finite element model make lumped-parameter model more favorable. The precise modeling of the planetary gears gains more insights in the dynamics of the system. With these models, Al-Shyyab and Kahraman [2931] solved the equations of motion semianalytically using multiterm harmonic balance method (HBM) in conjunction with inverse discrete Fourier transform and Newton–Rapson method. Chen et al. [32] investigated the dynamic response of the planetary gears with consideration of the flexibility of the internal ring gears. Chapron et al. [33] analyzed optimum profile modifications in helical and double-helical planetary gears (PGTs) with regard to dynamic mesh forces. Ericson and Parker [34] proved the accuracy of lumped-parameter and finite element models by experimentation and highlighted several design and modeling characteristics of planetary gears.

Compared with the spur planetary gears, the helical planetary gears (HPGs) are more competent and smooth-operated. Nevertheless, for a HPG system, the axial dynamic forces generated by gear meshing exist ubiquitously. Thus, the axial deformations of gears should be considered in the dynamic behavior investigations. Also, numerous analytical models on single-stage helical gear dynamics can be found in the literature while research on multistage helical gears is quiet limited. The coupling stiffness amongst the adjacent stages should be considered.

In the present paper, the main contributions of this paper contain two aspects. Firstly, the equations of motion of two-stage HPG used in cranes with consideration of axial deformations are derived. And secondly, the vibration characteristics of two-stage systems are investigated. The lumped-parameter model with three-dimensional motions is proposed to investigate the free vibration characteristics of HPGs. Modal properties are observed and categorized through the numerical simulation. Characteristics of each type of vibration modes are summarized. The unique properties of these vibration modes are then mathematically proved.

2. Three-Dimensional Dynamic Model of Planetary Gear

A two-stage HPG shown in Figure 1 is considered in this study. It is composed of a 2K-H differential planetary gear train and a quasi-planetary gear set. The input power is split into two paths when it is transmitted to the differential gear set. One is the ring of the differential gear set and the other is the sun gear of the quasi-planetary gear set. Finally, the power is output by rings of two stages. The closed-form structure effectively extends the range of transmission ratio and efficiency.

The three-dimensional discrete model of the example system is proposed. Component bearings, gear meshes, and interactions between two stages are all modeled by linear springs. Three-dimensional translation and the axial rotation are considered for any member as shown in Figure 2. The damping is not shown in the model. By choosing reasonable stiffness, any stage of the system in Figure 1 can be represented by Figure 2. The stiffness which equals a bearing or spline support stiffness is dependent on the specific kinematic configuration of the stage () [35].

To obtain the relative displacements of members accurately, two kinds of three-dimensional coordinate systems shown in Figure 3 are created. One is the three-dimensional absolute coordinate system , and the other is the dynamic coordinate system that rotate about the absolute coordinate. and are the theoretical installation center of sun gear and carrier, respectively. These two points superpose each other when the installation error is ignored. The unit vectors of the absolute coordinate system and the dynamic coordinate system are signed as and , respectively. According to the coordinate transformation theory, three-dimensional translation relationships of members can be drawn from Figure 3. where is the three-dimensional component of the displacement vector in the absolute coordinate system and is the component of the displacement vector in the dynamic coordinate system. is the angular speed.

The corresponding acceleration of members in those two coordinate systems is

The component of absolute acceleration of planet in the dynamic coordinate system is

2.1. Formulation of the Equivalent Displacements

The structural member is denoted by the subscript of physical symbols and the superscript denotes the structural member that is in stage . The following analyses coincide with the rules that the torsional displacement is positive when it is counterclockwise and the relative displacements for gear meshes are defined positive when compressed.

The relative displacement of the sun-planet mesh normal to tooth surface is

Similarly, the relative displacement of the ring-planet mesh normal to tooth surface is

The relative compression between planet and carrier along the line of , , , and is defined as

Besides, the relative compression between planet and carrier along the line of and is given as

The ring gear of the 1st-stage PGT and the ring gear of 2nd-stage PGT are both a part of the shell, and they are coupled by the shell. The meshing position of each stage is different. Because the rotation speeds of planets of each stage differ, there should be a relative torsional displacement between the contact positions of rings in two stages. Similarly, the carrier of the 1st-stage PGT and the sun gear of the 2nd PGT are linked by the spline. The 3D model of the two-stage planetary gear set is shown in Figure 4. Relative displacements of the adjacent members should be analyzed. Thus, the equivalent displacements of members in adjacent stages are written as

2.2. Formulation of Meshing Forces

The normal meshing forces of sun-planet and ring-planet in stage i are where and are the piecewise linear backlash function of sun-planet and ring-planet, respectively, and they are defined as

2.3. Dynamic Equations of Motion

The displacement vector is written aswhere , ; .

The equations of motion of the two-stage helical planetary gear set under consideration of the rotational, translational, and axial displacements (a total of degrees of freedom) are derived by Newton’s second law and the Theorem of Moment of Momentum. The dynamic equations of motion of the example system are listed as follows.

(1) Equations of the carrier in the 1st stage are

(2) Equations of the ring gear in the 1st stage are

(3) Equations of the sun gear in the 1st stage are

(4) Equations of planets in the 1st stage are

(5) Equations of the carrier in the 2nd stage are

(6) Equations of the ring gear in the 2nd stage are

(7) Equations of the sun gear in the 2nd stage are

(8) Equations of the planet in the 2nd stage arewhere the summation index ranges from 1 to in equations of stage .

Assembling the dynamic equations in matrix form yieldswhere M is the positive definite mass matrix and is the velocity matrix, is the bearing stiffness matrix, and is the meshing stiffness matrix. and are the corresponding damping terms. G is the gyroscopic matrix and is the stiffness matrix induced by gyroscopic effect. The matrix components are given in Appendix.

3. Vibration Modes of the Example System

The damping and external excitation are ignored to investigate the free vibration characteristics of the system. Besides, Coriolis acceleration of gears is neglected and all gyroscopic effects are excluded when the carrier speed is low [3]. These assumptions are valid for the example system, since the system is used in cranes with low velocity of carriers. The simplified equations of motion and associated eigenvalue equations are where is the natural frequency and the corresponding modal vectors are expressed as . The mesh stiffness of the system is calculated by Ishikawa tooth deformation calculation formula and stiffness calculation formula. Because rings of two stages are coupled by the shell and the carrier of the 1st-stage PGT and the sun gear of the 2nd PGT are linked by the spline. The coupling stiffness amongst rings of two stages equals the stiffness of the shell between two rings and stiffness between the carrier of the 1st-stage PGT and the sun gear of the 2nd-stage PGT is dependent on the spline. The system parameters are listed in Table 1.

It is assumed that planets in both stages are equally spaced. The systems with and = 3~9 are simulated here to illustrate the structured vibration modes of two-stage HPGs. According to the modal theory, natural modes are grouped according to the multiplicity of the natural frequencies. The natural frequencies and mode types of the system are shown in Table 2.

From Table 2, the vibration modes can be recognized as 3 groups: the axial translational-rotational mode, the radical translational mode, and the second-stage planet mode (). Following results can be concluded.(a)The first-order natural frequency is and the associated vibration mode is the rigid body mode. The multiplicity of this mode is 1. As shown in Table 2, this mode can be concluded in the axial translational-rotational mode.(b)Twenty natural frequencies always have multiplicity for different and . The corresponding vibration mode is the axial translational-rotational mode. This mode type is mainly characterized by the motions of central members which only rotate and translate axially. Furthermore, deflections of planets in the same stage are identical. The typical axial translational-rotational mode with , , and  Hz is shown in Figure 5.(c)There are always fourteen pairs of natural frequencies with multiplicity for various and , and the corresponding vibration mode is the radical translational mode. In these modes, all central members have pure translation radically. Figure 6 shows the mode shape of radical translational mode with , , and  Hz.(d)As shown in Table 2, four degenerate natural frequencies with multiplicity are always calculated when . Associated modes are the 2nd-stage planet modes. In these modes, all central members are stationary and only planets in the 2nd stage have movements. The planets move in all four degrees of freedom (DOFs), and their modal deflection is a scalar multiple to that of the arbitrarily selected first planet. Figure 7 illustrates the typical second-stage planet mode.(e)It is concluded that the natural frequency varies monotonically as additional planets in stage two introduced.

In most planetary gears, the number of planets in any stage is more than 3. This structure brings the systems distinguishing advantage of power-split. In the following discussion, the eigensolution properties summarized in the example are analytically shown to be true for general HPGs.

4. Analytical Characterization of Vibration Modes

The eigenvalue matrix equations (21) are expanded to () groups of equations associated with the individual memberswhere 0 has dimension 4 × 1.

4.1. Axial Translational-Rotational Mode

According to the above vibration mode analysis, the axial translational-rotational mode has the following characteristics:

(a) The corresponding multiplicity of natural frequency is 1.

(b) The radical deflections of carrier, ring gear, and sun gear in two stages are zero

(c) The modal deflections of planets in the same stage are identical

The eigenvectors of the 1st stage and 2nd stage can be expressed as

Inserting (25)–((27a) and (27b)) into (23a) yields

Similarly, inserting (25)–(27a) into (23b)-(23c) yieldswhere , , , , , and . Equation (23d) can be simplified as

Inserting (25)–(27b) into (24a), we get

In the same way, (24b)-(24c) reduce to

Applying the mode characteristics (b) and (c), (24d) becomes

The above equations imply that (23a)–(24d) reduce to a 20-degree-of-freedom eigenvalue problem consisting of (28a)–(35). The modal vector q is determined by the reduced eigenvalue Twenty eigensolutions are derived from this reduced eigenvalue problem. Therefore, there are 20 axial translational-rotational modes and the corresponding natural frequency is distinct.

4.2. Radical Translational Mode

To sum up, the radical translational mode has the following characteristics:(a)The natural frequency multiplicity of the radical translational mode is 2. There is a pair of orthonormal vibration modes associated with each natural frequency(b)The axial translation and rotation of central members in both stages are zero. Besides, the relation of the modal deflections of central members in a pair of orthonormal vibration modes is(c)The modal deflections of planets of same stage in a pair of orthonormal vibration modes can be expressed as where and are the modal deflections of planets in a pair of orthonormal vibration modes.

Thus, a pair of candidate orthonormal radical translational modes of the 1st stage is given by

Analogous to the analyses of the axial translational-rotational mode, substituting (38)–(40b) into (23a) yields

Similarly, the following (42a)–(43b) are derived by inserting (38)–(40b) into (23b)-(23c).

Inserting and into (24d), we getwhere and are the introduced notations. Applying the mode characteristics of the radical translational modes (b) and (c), the above equations can be expressed as

The eigenvalue problems of the 1st stage can be reduced to fourteen independent equations: six equations from (41a)–(43b) and eight equations from (45a) and (45b). Thus, seven natural frequencies with multiplicity two are obtained. The reduced eigenvector of the 1st stage can be expressed as , . With similar operations on the eigenvalue matrix equations (24a), (24b), (24c), and (24d), seven natural frequencies with multiplicity two also can be derived from the reduced-order eigenvalue problems of the 2nd stage and the reduced eigenvector of the 2nd stage can be simplified as , . To sum up, there are always fourteen natural frequencies with multiplicity 2 for the two-stage HPGs, and the corresponding vibration modes are the radical translational modes.

4.3. Planet Mode

The planet mode is composed of the 1st-stage planet mode and the 2nd-stage planet mode. Taking the 1st-stage planet mode as an example, we analytically investigate the unique characteristics of the planet mode.

The typical characteristics of the 1st-stage planet mode are(a)The natural mode multiplicity of the 1st-stage planet mode is .(b)The modal deflections of central members in both stages and planets in the 2nd stage are zero. Thus, the displacement subvectors can be written as(c)Only planets in the 1st stage have motions. The deflections of diametrically opposed planets are identical. Besides, the motions of each planet are a scalar multiple of the motions of an arbitrarily selected first planet. The displacement subvector of planet in the 1st stage is given bywhere are constants and .

The eigenvectors of the 1st stage can be expressed aswhere the zero vectors are 4 × 1.

Substituting (46)-(47) into (23a), (23b), (23c), and (23d) gives

To satisfy (49)–(51), three independent constraint equations can be obtained.

It is obvious that there is no solution if . Four natural frequencies are calculated from (52). According to (48), independent solutions can be deduced with determined from (53) only if .

Similarly, in the 2nd-stage planet modes, only planets in stage two have modal deflections. Without loss of generality, we also analytically investigate the reduced-order eigenvalue problems of the 2nd-stage planet modes. With similar operations on the eigenvalue matrix equations of the 2nd-stage members, it is concluded that four natural frequencies with a multiplicity of are always obtained and the associated vibration modes are the 2nd-stage planet mode.

5. Conclusions

In this study, a three-dimensional dynamic model of the two-stage HPGs with consideration of axial and radical deformations is established. Two three-dimensional coordinate systems are created to analyze the relative displacements of members. Dynamic equations of motion which can be used to analyze HPGs with variant number of planets in two stages, different planet phasing, and spacing configurations are derived by Newton’s second law and the Theorem of Moment of Momentum. Subsequently, the natural frequencies and vibration modes of the two-stage HPGs with equally spaced planets are simulated numerically. Finally, the unique properties of these vibration modes are mathematically proved. The main properties of the two-stage HPGs are summarized below:(1)Vibration modes of the two-stage HPGs can be categorized as five classes: the rigid body mode, the axial translational-rotational mode, the radical translational mode, and the 1st-stage and the 2nd-stage planet mode.(2)Axial translational-rotational modes: twenty natural frequencies always have multiplicity for various and . The associated mode is the axial translational-rotational mode. All central members only rotate and translate axially in this vibration mode. Deflections of planets in the same stage are identical.(3)Radical translational modes: there are always fourteen pairs of natural frequencies with multiplicity for different and , and the corresponding vibration mode is the radical translational mode. In this mode, all central members have pure radical translations.(4)Planet modes: the planet modes of the two-stage HPGs are composed of two types, the 1st-stage planet mode and the 2nd-stage planet mode. Four natural frequencies with multiplicity are calculated only if . The associated vibration modes are the th stage planet mode in which only planets in the th stage have motions while other members remain stationary.(5)The natural frequency with multiplicity varies monotonically as additional planets in both stages are introduced.

Appendix

The related damping matrix C has similar form with stiffness matrix K by changing the stiffness parameters into corresponding damping parameters.

Nomenclature

:Translational displacement of component in the th stage in direction
:Translational velocity of component in the th stage in direction
:Translational acceleration of component in the th stage in direction
:Rotational displacement of component in the th stage in direction
:Rotational velocity of component in the th stage in direction
:Rotational acceleration of component in the th stage in direction
:Translational displacement of the planet in the th stage in direction
:Translational velocity of the planet in the th stage in direction
:Translational acceleration of the planet in the th stage in direction
:Rotational displacement of the planet in the th stage in direction
:Rotational velocity of the planet in the th stage in direction
:Rotational acceleration of the planet in the th stage in direction
:Number of planets in stage
:Moment of inertia of the component in stage
:Mass of the component in stage
:Circumferential angle of the planet in stage
:Backlash of stage
:Torsional angular displacement of component in the th stage
:Base circle radius
:Mesh stiffness (damping) of the th sun-planet or ring-planet
:Torsional stiffness (damping) of component in the th stage
:Bearing stiffness (damping) of component in the th stage in direction
(:Interstage coupling stiffness (damping) in direction ()
, :Pressure angle of sun-planet and ring-planet of stage
:Helix angle of stage
:Import (output) torque.
Subscripts
:Carrier gear
:Ring gear
:Sun gear
:th planet gear.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 51175299).