Abstract

Dynamic excitation caused by time-varying meshing stiffness is one of the most important excitation forms in gear meshing process. The mesh phase relations between each gear pair are an important factor affecting the meshing stiffness. In this paper, the mesh phase relations between gear pairs in an encased differential gear train widely used in coaxial twin-rotor helicopters are discussed. Taking the meshing starting point where the gear tooth enters contact as the reference point, the mesh phase difference between adjacent gear pairs is analyzed and calculated, the system reference gear pair is selected, and the mesh phase difference of each gear pair relative to the system reference gear pair is obtained. The derivation process takes into account the modification of the teeth, the processing, and assembly of the duplicate gears, which makes the calculation method and conclusion more versatile. This work lays a foundation for considering the time-varying meshing stiffness in the study of system dynamics, load distribution, and fault diagnosis of compound planetary gears.

1. Introduction

Planetary gear transmission is widely used in aviation, automobile and other industrial fields because of its strong load-carrying capacity and small size. However, the vibration and noise caused by the dynamic meshing force and dynamic supporting force of the system affect its reliability and service life. Equally distributing the load to each planet and avoiding excessive vibration are the goals of design and manufacture of planetary gear train. The mesh phase relations between gear pairs are an important factor affecting the meshing stiffness of the system. It has a great influence on the load distribution and vibration of the system and has always been the focus of researchers.

Based on an improved potential energy method, a time-varying mesh stiffness calculation model considering crack growth was proposed. The meshing stiffness of normal gear pair and gear pair with different levels of crack is calculated by traditional method, proposed method, and ANSYS method, respectively. The influence of the crack level of the gear on the dynamic response was studied by using a one-stage gear system dynamic model [1]. Cui calculated the gear tooth thickness decrements by straight-affecting line method and parabolic-affecting line method. Based on the universal tooth profile equation, the meshing stiffness of gear pairs with different crack levels was calculated, and their effects on the vibration response are studied by fault detection indicators [2]. Considering three possible contact situations, the meshing stiffness under different torques is calculated and compared with the stiffness obtained by ISO 6336. The results of internal and external meshing stiffness calculated by various models are compared [3]. Based on the potential energy principle, Chen deduced the stiffness calculation model considering the modification and calculated the stiffness of spur gears under different modification coefficients [4].

Kahraman established a lumped parameter dynamic model of a single-stage planetary gear train and obtained the natural modes and forced vibration response caused by static transmission error. The influence of mesh phase relation on the dynamic characteristics of four-planet system was studied [5]. Kahraman also established a nonlinear time-varying dynamic model of planetary transmission considering tooth profile error and time-varying meshing stiffness, and analyzed the influence of errors on dynamic load-sharing coefficient [6]. Based on the harmonic decomposition of meshing forces acting on the sun-planet and ring-planet mesh, Parker and Abbasha studied the suppression of mesh phase relations on different vibration types of planetary gears. The nonlinear dynamic characteristics of system considering mesh phase difference are studied by using analytical model and finite element model [7, 8]. The results show that the rules of suppressing the vibration caused by mesh phase difference are applicable to the nonlinearity caused by contact loss of gear teeth [9].

Parker defined the mesh tooth number function to describe the mesh phase relations, deduced the calculation method of the phase difference between the gear pairs, and calculated the phase difference of each mesh gear pair in the single-stage planetary gear train of OH58 helicopter main reducer when the planet gear rotates in the positive and negative directions [10]. Chen Yong studied the relationship between mesh phase difference and torsional vibration and carried out experimental research under different rotational speeds and loads. The results show that the design considering mesh phase difference can effectively reduce the vibration and noise levels of planetary gear system, and the tooth profile contact ratio has no significant effect on the vibration and noise levels in the presence of phase differences [11].

Guo and Parker [12] proposed a method to accurately define the phase difference of the compound planetary gear system. By introducing the concepts of relative phase and reference datum, the mesh phase relations between any two gear pairs in a compound planetary gear train were derived. Wang [13] studied the influence of mesh phase on the vibration of spur planetary ring gear by using the superposition principle and Fourier series method and analyzed the relationship between modal characteristics and mesh phase difference. Gawande measured the noise level of the planetary gears and found that the noise level of the planetary gears with mesh phase difference was significantly lower than that of the planetary gears without mesh phase difference [14, 15].

Chaari [16] established a dynamic model of a single-stage planetary gear train with consideration of the mesh phase. The dynamic responses of healthy planetary gear train and planetary gear train containing defective gear tooth were compared in time and frequency domains, and the validity of Wigner-Ville method for fault diagnosis of planetary gear train was verified. Li [17, 18] studied the mesh phase relations of two-stage planetary gears with meshed-planet gears. On this basis, the time-varying meshing stiffness was calculated, and the simulation analysis and experiment of gear tooth fault signal were compared. Hu added the mesh phase difference introduced by helical angle to that of spur gear to obtain the mesh phase difference of the helical gear, established a load sharing calculation model considering the mesh phase, and calculated the load sharing factors of the system under different phase relations. The trajectory of the sun gear [19] was obtained, which is in good agreement with the experimental results [20, 21].

In these literatures, pitch points are used as reference points of mesh phase difference, and researches are mainly for the single-stage planetary gear systems. There is no specific calculation method for the mesh phase difference of the compound planetary gears with stepped-planet gears. In this paper, the starting point of the meshing is taken as the reference point to make the calculation of phase difference more concise. The calculation method of phase difference of the compound planetary gear train with stepped-planet gears and duplicate central gears is derived. The gear modification is considered in the derivation process, which makes the conclusions in this paper more versatile.

2. Definition of Mesh Phase Difference

The research object of this paper is an encased differential planetary gear train shown in Figure 1 which can be regarded as a combination of two different subsystems: the fixed-shaft gear train with stepped-planet gears, that is, the encased stage, and the differential planetary gears is the differential stage. The input power is split into two paths as it is transmitted to the system: one is transmitted to the ring gear r1 through the encased stage and the other is transmitted to the carrier c2 and gear ring r2 through the differential stage. When the number of teeth of each gear in a planetary gear train conforms to a specific relationship, the output shafts cs and rs can be output in reverse direction at the same rotational speed. The number of stepped planets ab and planets p is M and N, respectively.

The gear pair sun gear s1 and stepped gear ai (i=1, 2, …, M) is defined as s1ai, and the definitions of other gear pairs in the system are similar. The angular velocity of the component h is (h=s1, s2, r1, r2, c2, a, b, p) and the counterclockwise is positive, and the mesh period of all gear pairs is determined by the following:

As shown in Figure 2, represents the mesh phase difference of the mesh B relative to the mesh A [10], where and is the contact ratio of the mesh g. Starting from the reference position, the gear pair g is in contact at the meshing starting point for the first time at , and the mesh phase difference can be obtained. where is a fractional function and is the nearest integer less than x.

3. Mesh Phase Relations of Encased Differential Gear Train

The calculation process of mesh phase difference of the encased differential planetary gear train is shown in Figure 3.

The initial assembly position of the compound planetary gears is shown in Figure 4. In the figure, MN is the theoretical line of action, A is the meshing starting point, P is the pitch point of the mesh, E is the point where the gear tooth exits contact, D is the intersection point of the reference tooth midline and the addendum circle, and the definitions of radius and angle are given in the Appendix.

The angle of the midline of the planet gear b1 reference tooth relative to that of the planet gear a1 reference tooth is . In order to meet the assembly conditions, the relative angles of gears b and a of each stepped-planet gear need to be the same. The angle of the midline of the sun gear s2 reference space relative to that of the sun gear s1 reference space is .

3.1. Calculation of Phase Difference between Adjacent Mesh

For the sun-planet mesh s1a1, the calculations follow from the geometry of Figure 4(a).herewhere is the actual center distance between the sun gear s1 and the stepped-planet gear a1.

It can be concluded that when the stepped-planet a1 rotates from the initial position by an angle , the meshing starting point is in contact for the first time at .

3.1.1. Calculation of Phase Difference

In order to meet the assembly requirements, it is advisable to design the installation position of planet gear p1 to be the same as that of the stepped-planet gear a1 circumferentially; that is, the center line coincides with the center line in the circumferential direction, as shown in Figure 4(b).

For the sun-planet mesh s2p1, geometrical relation can be derived from Figures 4(b) and 4(c). herewhere is the actual center distance between the sun gear s2 and the planet gear p1.

When the planet p1 rotates from the initial position by an angle , the meshing starting point is in contact for the first time at .

The mesh phase difference between gear pairs s2p1 and s1a1 is

3.1.2. Calculation of Phase Difference

As shown in Figure 4(a), is the starting point of the involute of the reference gear teeth on the base circle, and the curves and are calculated as follows:where

Equations (14)-(16) givewhere

From the initial position, the meshing starting point of gear pair r1b1 is in contact for the first time at .

The mesh phase difference between gear pair r1b1 and gear pair s1a1 is

3.1.3. Calculation of Phase Difference

Figure 4(b) shows the following for length of the curves:

According to the above equations (21)-(23), we have

The phase difference between internal gear pair and external gear pair of the differential stage iswhere the tooth thickness on the base circle .

3.1.4. Calculation of Phase Difference

The number of teeth meshed when the sun gear s1 rotates one revolution is . When the sun gear s1 rotates from the position of the stepped-planet gear a1 to the position of the stepped-planet gear ai, the angle of rotation and the number of teeth meshed are

From (26) and (27), the mesh phase difference between gear pairs s1ai and s1a1 is

3.1.5. Calculation of Phase Difference

Similarly, when the ring gear r1 rotates from the position of the planet b1 to the position of the planet gear bi, the angle of rotation and the number of teeth meshed are as follows:

The mesh phase difference between the gear pairs r1bi and r1b1 obtained by (29) and (30) is

3.1.6. Calculation of Phase Difference

It is known that the angle between the j-th planet gear and the 1st planet gear in the differential stage is

When the sun gear s2 rotates over the angle relative to the carrier c2, the j-th planet gear moves to the initial position of the 1st planet gear, and the time required is

Considering that the time difference caused by mesh phase difference is , the time can also be expressed aswhere is an integer. Available from (1), (33) and (34), the mesh phase difference between the gear pairs s2pj and s2p1 is

3.1.7. Calculation of Phase Difference

Similar to the above, when the ring gear rotates relative to the carrier c2, the j-th planet gear moves to the initial position of the 1st planet gear, and the time required is

Considering that the time difference caused by mesh phase difference is , we obtain the following:where is an integer. The mesh phase difference between the gear pairs r2pj and r2p1 can be obtained by introducing (1) into (37).

3.2. Calculation of Comprehensive Mesh Phase Difference

Through the above analysis, the mesh phase differences between adjacent gear pairs are obtained. In order to unify the time difference caused by mesh phase difference to the absolute time of the system, it is necessary to synthesize the adjacent phase difference to get the phase difference of each gear pair relative to the base referred mesh. In the following analysis, the mesh s1a1 is selected as the base referred mesh of the system.

3.2.1. Calculation of Phase Difference

After the contact of the meshing start point of sun-planet mesh s2p1, the time at which the meshing start point of ring-planet mesh r2p1 contact for the first time is

Relative to the base referred mesh s1a1, the mesh phase difference of the gear pair r2p1 is

3.2.2. Calculation of Phase Difference

After the mesh r1b1 contacts at meshing starting point, the first contact time of meshing starting point of the mesh r1bi is

The mesh phase difference of the gear pair r1bi relative to the base referred mesh s1a1 is

3.2.3. Calculation of Phase Difference

After the mesh s2p1 contacts at meshing starting point, the first contact time of meshing starting point of the mesh s2pj is

The mesh phase difference of the gear pair s2pj relative to the base referred mesh s1a1 is

3.2.4. Calculation of Phase Difference

After the mesh r2p1 contacts at meshing starting point, the first contact time of meshing starting point of the mesh r2pj is

The mesh phase difference of the gear pair r2pj relative to the base referred mesh s1a1 is

Since the starting point of the referred gear pair s1a1 contacts at , the time at which the meshing starting point of the gear pair g contacts for the first time is

Then the mesh tooth variation function in the system can be expressed as [12]where is a periodic time-varying mesh tooth variation function of the gear pair g, corresponds to the meshing starting point, and t is absolute time of the system.

4. Example Calculation of Phase Difference in Encased Differential Gear Train

Table 1 shows the parameters of the encased differential planetary gear train with M=6, N=6. For the convenience of processing and assembly, it is usually guaranteed that the midline of referred gear tooth or space of the duplicate gear is coincident during machining, that is, the offset angle , . The mesh phase difference of all gear pairs in the system relative to the referred gear pair s1a1 is calculated as shown in Table 2.

When the input shaft speed is (r/min), the mesh period of all gear meshes in the example system can be obtained.

Without loss of generality, let the starting time of the reference gear pair s1a1 be in contact. The gear pair s2pj mesh tooth variation function changes with time as shown in Figure 5, and the other gear pairs in the system can be similarly unified to the system time t.

The time-varying meshing stiffness of each gear pair in a meshing period calculated by theoretical method or finite element method combined with the meshing time difference caused by the phase difference can be applied to the dynamic model of the system, which can accurately describe the meshing stiffness of all gear pairs in the model at any time and carry out the dynamic analysis of the gear system.

5. Conclusions

The accurate solution of the mesh phase difference of all gear pairs is an extremely important part of the study of the dynamics, fault diagnosis, and load sharing characteristics of compound planetary gear train. In this paper, for the encased differential compound planetary gear train widely used in coaxial twin-rotor helicopters, taking the meshing starting point as the reference point, considering the gear modification, the calculation method of mesh phase difference is studied.

Firstly, the initial position of the system and the referred gear pair of the system are defined, and the time when the referred gear pair is in contact at the meshing starting point is calculated. According to the angle relations of the meshing starting points and the distance relations of the starting points on the meshing line, from the initial position, the time at which the gear pair adjacent to the base referred gear pair first contacted at the starting point is calculated, and the phase difference can be obtained by the meshing time difference between the gear pairs. The phase difference between planet gears is calculated by the relations between the number of meshing teeth and the meshing time caused by the symmetry of planetary gear train.

Based on the phase difference of adjacent gear pairs, the calculation method of the phase difference of each gear pair relative to the referred gear pair is derived. The derivation process considers the processing and assembly requirements of the duplicate gear, and a specific compound planetary gear train is calculated as an example, which lays a foundation for the accurate introduction of time-varying meshing stiffness into the dynamic analysis of compound planetary gear train.

Nomenclature

:Number of teeth of gear
:Pressure angle on addendum circle of gear
:Working pressure angle of gear pair where central gear is located
:Involute function
:Involute function
:Tooth thickness at reference circle of gear
:Reference circle radius of gear
:Base circle radius of gear
:Addendum circle radius of gear
:Base circle pitch of gear
:Tooth thickness at base circle of gear
:Modification coefficient of gear
:Modulus of the gear .

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China [grant number 51775265.