Abstract

A new dynamic model for a two-input two-path split torque transmission system which considers meshing error, time-varying meshing stiffness, and meshing-in impact is proposed. Time-varying meshing stiffness and meshing-in impact of each gear pair are accurately calculated based on tooth contact analysis and loaded tooth contact analysis. Equivalent displacements of eccentricity error and installation error along the meshing line of second- and third-stages gears are derived. The modified tooth surface of a third-stage double-helical gear is obtained by optimizing the amplitude of static loaded transmission error and meshing-in impact via nondominated sorting genetic algorithm-II (NSGA-II). Influence of modification on load sharing and dynamic load characteristics of split torque transmission system is investigated. The results indicate that the system’s dynamic meshing force increases when meshing-in impact is accounted for, which is unfavorable for the transmission. Following the modification of a double-helical gear, the dynamic load characteristics of the split torque transmission system are significantly improved, while its load sharing characteristics are improved to a certain extent.

1. Introduction

Representing a new and advanced configuration designed to replace the planetary transmission, the split torque transmission system (STTS) is increasingly employed in helicopter main reducer [15]. STTS is characterized by a driving pinion which simultaneously meshes with two gears, and each path only carries an equal half of the original load under ideal conditions. Due to the power closed-loop features of the STTS and the inevitable manufacturing and installation errors of the gears, a problem of uneven load distribution in each path of the STTS is inevitably caused. The unequal load sharing of the STTS increases the vibration and noise, while simultaneously affecting the STTS’ reliability and durability.

Many research works have conducted numerous investigations on natural and load sharing characteristics of planetary and star transmission systems. Kahraman [6, 7] established a dynamic model of a planetary transmission. The author predicted natural modes of the system and investigated the influence of manufacturing and installation error on dynamic load coefficient of the system. Sondkar et al. [8] investigated the influence of double-helical (DH) gear stagger angle on dynamic response of a double-helical planetary gear system and found that the left and right sides of the DH gear can bear different dynamic load due to the effect of the stagger angle. Liu et al. [9] analyzed the influence of the profile error before and after the gear tooth separation on the vibration and dynamic meshing force of the herringbone planetary gears. Li et al. [10] established a model to predict the reliability of a helicopter planetary gear train under partial load and concluded that when the twin-engine helicopters are running on only one engine, the likelihood of teeth fatigue fracture is increased. Ren et al. [11] established a generalized dynamic model for a herringbone planetary gear train which considers the manufacturing eccentric errors of components and tooth profiles errors. Mo et al. [12, 13] investigated the influence mechanism of multicoupling errors, flexible support stiffness, and floating components on load sharing characteristics of a herringbone planetary transmission system. Mo et al. [14] proposed a new dynamic model for double-helical star gearing system and obtained the vibration model, natural characteristics, and dynamic response of the system. Wang et al. [15] established a dynamic model of a star herringbone gear transmission system and discovered that the system has four typical vibration modes. Wang et al. [16] investigated the influence of bearing’s position and permutations on the load sharing performance of star gearing system. Bechhoefer et al. [17] processed the vibration signatures of a split torque gearbox through a number of gear analysis algorithms to quantify the gear fault performance. Gui et al. [18] analyzed the backlash influence on load sharing coefficient of the two-input cylindrical gear STTS. Dong et al. [19] established a dynamic model of the power-split transmission and obtained the frequency domain and time domain responses of the system. Zhao et al. [20] proposed a universal mathematical design method for a torque-split gear transmission and verified the correctness of the method through different types of gears. Hu et al. [21] established the gear-shaft-bearing dynamic model of two-path STTS and analyzed the influence of factors such as shaft angle and DH gear stagger angle on the system’s dynamic characteristics and load sharing characteristics. Jin et al. [22] analyzed the influence of law and weight of backlash and center distance error on STTS load sharing characteristics.

Meshing error [1114], time-varying meshing stiffness (TVMS) [2326], and meshing-in impact [2730] are important internal excitation of gear transmission system. Periodic TVMS excitation is the most obvious feature of gear system vibration that distinguished it from most mechanical systems. The finite element method and numerical analysis method are usually used to calculate TVMS. The former has high calculation accuracy but low calculation efficiency, while the latter has high calculation efficiency but low calculation accuracy. In this paper, the loaded tooth contact analysis (LTCA) method is used to calculate the TVMS of gear pair under static condition [3133]. Gear modification has proven to be an effective method for reducing vibration and noise of the gear system. The modification is mainly carried out around the three elements of modification length, modification amount, and modification curve [3437]. DH gear is generally used as the final-stage transmission in STTS because of its smooth transmission, compact structure, large carrying capacity, and self-balancing axial force. The final-stage transmission reduction ratio of approximately 10 : 1 to 14 : 1 can be achieved. Furthermore, transmission performance of the final-stage transmission has a crucial impact on the performance of the entire STTS. The existing research results show that the dynamic characteristics of a pair of DH gear transmission system can be significantly improved through the modification design of DH gear [3840]. However, there are few studies on the influence of DH gear modification on load sharing and dynamic load characteristics of STTS.

In this paper, TVMS and meshing-in impact of each gear pair are calculated accurately via tooth contact analysis (TCA) and LTCA. Dynamic model of a two-input two-path STTS accounting for meshing error, TVMS, and meshing-in impact is established. The dynamic model can be employed to analyze the load sharing and dynamic load characteristics of STTS before and after modification under any working conditions. The influence of gear eccentricity and installation error on load sharing characteristics of the system is studied. By optimizing the amplitude of static loaded transmission error and meshing-in impact, the modification design of DH gear is carried out, and the topological modification tooth surface of the DH gear is constructed. Finally, the modification influence on load sharing and dynamic load characteristics of two-input two-path STTS is investigated. This research has important theoretical guidance and significance for the load sharing design of STTS and application of DH gear modification technology in helicopter main reducer.

2. Internal Excitation Analysis of Two-Input Two-Path STTS

2.1. Physical Model of Two-Input Two-Path STTS

Figure 1 shows the two-input two-path STTS for helicopter main reducer. The helicopter main reducer has two engines which transmit power to the input gear, one on the left and one on the right side. The input power on both sides is collected in the output gear, and the power is transmitted to the main rotor and tail stabilator through the rotor shaft. The two-input two-path STTS adopts three--stage transmission. The first stage adopts a spiral bevel gear pair to achieve the reversal of the transmission system. The second stage adopts a spur gear pair; i.e., one spur pinion simultaneously meshes with two large spur gears to achieve the power split. The third stage adopts a DH gear pair; i.e., two small DH pinions simultaneously mesh with one large DH gear to achieve power convergence. For the convenience of the following description, the gears of the two-input two-path STTS are numbered in sequence from 1 to 8. The gears on the left and right sides of the first stage, second stage, and third stage are represented by zij (i = L, R, j = 1…7). A large center output DH gear is designated as the gear number 8, i.e., z8.

As shown in Figure 2, coordinates of two-input two-path STTS are established. Generalized STTS coordinate system is O-XYZ, while the axis direction of gear 8 is defined as the Z-axis direction. Parameters σ1 and σ2 represent system installation angles, and ψ represents the system layout angle. The angle between the center line of each gear pair and positive X-axis direction is indicated by γ and the corresponding subscript, i.e., γL34 and γR34. The angle between the meshing line of each gear pair and the positive direction of the X-axis is denoted as η and the corresponding subscript, i.e., ηL34 and ηR34. The angular displacement is taken counterclockwise from the positive X-axis as “+.” The remaining angles are not marked in Figure 2, and the local coordinate systems of the first-stage gear pair are omitted.

2.2. Meshing Error of Second- and Third-Stage Gears Analysis

Since the load sharing of two-input two-path STTS is independent of the first-stage gear pair error, only equivalent displacement formula of the eccentricity and installation error along the meshing line of second- and third-stage gears is studied in this paper.

Second-stage spur gear pair (denoted as 34) of left-input two-path STTS is considered as an example. The projection relationship between the eccentricity error EL3, EL4 and installation error AL3, AL4 of spur pinion zL3 and the spur gear zL4, and the meshing line MLL34 is shown in Figure 3. Parameters λL3, λL4 and μL3, μL4 represent phase angles of the eccentricity and installation errors, respectively. The eccentricity and installation errors of the spur pinion zL3 and the spur gear zL4 are projected to the meshing line to obtain their respective equivalent meshing errors. Then, the equivalent meshing errors of two gears are combined to obtain the equivalent cumulative meshing error on the meshing line of gear pair 34.

Therefore, the equivalent cumulative meshing error ei34 (i = L, R) on the meshing line of the second-stage spur gear pair 34 can be expressed as

Similarly, the equivalent cumulative meshing error on the meshing line of the second-stage spur gear pair 35 and the third-stage DH gear pair 68 and 78 can be expressed aswhere ωi3, ωi4, and ωi5 are the angular velocities of spur gear zi3, zi4, and zi5, respectively. ωi6, ωi7, and ω8 are the angular velocities of DH gear zi6, zi7, and z8. Parameters αts and αth are the transverse pressure angle of spur gear and DH gear, respectively.

2.3. Calculation of Time-Varying Meshing Stiffness of Gears

The approach for LTCA model is used to solve the comprehensive meshing stiffness of gear pairs. The LTCA technology organically combines the geometric analysis and mechanical analysis of the gear. It has the following advantages [31]: (1) the whole tooth profile of large and small gear including tooth root transition surface is obtained by cutter tools generation, and the tooth profile is not simplified. At the same time, the flange structure is considered, and the finite element mesh is automatically generated by programming. The shape of the element is relatively regular, so the gear model itself has high accuracy. (2) The finite element method only needs one calculation to solve the nodal flexibility matrix of working tooth surface, while the normal flexibility matrix of contact line discrete points at different meshing positions in the meshing period can be obtained by reinterpolation. Compared with commercial finite element software, it takes less time and has higher efficiency. (3) It can reflect the micron-level (tooth modification) geometric feature changes of the tooth surface. (4) It can reflect the influence of the micron-level geometric features of the tooth surface on the mechanical properties (strength and vibration).

The finite element models of a driving gear for spiral bevel gear pair, spur gear pair, and DH gear pair are shown in Figure 4.

Taking DH gear as an example, the LTCA model of a DH gear is shown in Figure 5, and the LTCA models of a spiral bevel gear and a spur gear are detailed in [32, 33]. The LTCA model is to simulate and analyze the gear pair under static conditions and then obtain the static transmission error (STE) and static loaded transmission error (SLTE) of the gear pair. A mesh period of DH gear is divided into different mesh positions. Under the action of the total normal load P, the normal deformation Z of different mesh positions can be obtained. The normal deformation Z is actually the SLTE expressed by the linear displacement. Converting it into the form of angular displacement, it can be expressed asHere, rb2 and βb are the base cylinder radius and base helix angle of the DH gear, respectively.

Then, the amplitude of SLTE can be expressed as

Finally, the comprehensive meshing stiffness k can be expressed aswhere Fn is normal force acting upon the tooth flank.

According to (3)∼(5), the SLTE and comprehensive meshing stiffness are two different description methods of the same nature of excitation. Thus, the two can be transformed into each other. Therefore, the fluctuation amplitude of the comprehensive meshing stiffness can be minimized by optimizing the amplitude of SLTE value in the subsequent optimization design.

The LTCA model can be used to accurately calculate the comprehensive meshing stiffness of each gear pair, and the time-varying meshing stiffness of each gear pair can be obtained by expressing it in the form of Fourier series, which can be directly introduced into dynamic equations of two-input two-path STTS.

2.4. Calculation of Gear Meshing-In Impact

In the actual transmission process of gear pair, due to the manufacturing error, installation error, and elastic deformation of gear teeth, the normal pitches of the driving pinion pb1 and the driven gear pb2 along the line of action are no longer equal, resulting in a relative base pitch difference. The essence of meshing-in impact is that the tooth pair that enters meshing-in advance produces a relative velocity difference along the meshing line at the meshing point. At this time, the position of the driven gear that is about to enter the mesh can be regarded as a slight angle of retraction based on the theoretical meshing position, which is defined as the reverse angle.

The key to calculate the meshing-in impact is to accurately obtain the position of the meshing point on the tooth surface. The two-input two-path STTS includes spiral bevel gear, spur gear, and double-helical gear. According to the geometric position relationship shown in Figure 6, for spur and DH gears with standard involute tooth surface, the reverse angle of the driven gear can be obtained by employing the reversal method. Then, the position of the meshing-in point of the standard tooth surface can be obtained [28]. In Figure 6, Δφkg is the reverse angle of the driven gear, and Δφkp is the reverse angle of the driving pinion. Point D is the actual meshing point, and the standard tooth surface contains meshing interference at this point. Point E is the normal engagement point, and E′ is the reversal point of the normal engagement point. rp and are pitch circle radii, and rbp and rbg are base circle radii of pinion and gear, respectively. rag is the addendum circle radius of the gear.

The exact position of the initial meshing-in point of engagement, D, can be accurately determined as follows:where φk = Z/rag,  = π/2−αt−∠PEOg, and αt is the transverse pressure angle.

Due to the complex tooth surface, it is difficult to determine the position of the initial meshing point through geometric analysis methods for spiral bevel gear and modified DH gear. In this paper, the reversal angle of the above driven gear is calculated based on TCA and LTCA technology. Then, the initial meshing point position is determined. As shown in Figure 7, point A is the intersection of STE and SLTE curves, which represents the theoretical meshing point. LTE2 is the SLTE value of the current meshing tooth pair (tooth pair 2) at the meshing point. TE1 is the STE value of the previous meshing tooth pair (tooth pair 1) to the corresponding meshing point. Parameter Δφ is the backward angle of the driven gear in theoretical engagement position, which can be expressed as

The rotation angle of the driven gear at the actual meshing point can be expressed aswhere φ2 is the theoretical rotation angle of the driven gear at the engagement position.

By employing a new angle, , the TCA of the tooth surface is performed again. Then, the position vector and normal vector of actual meshing-in impact point on the tooth surface of the driving pinion and driven gear can be obtained.

Based on the position of the meshing point, the relative normal velocity of two teeth surfaces at the meshing point is calculated, which can be expressed aswhere and are normal velocities of pinion and gear at the meshing point, respectively.

Then, the meshing-in impact of gear pair is calculated via impact mechanics method. The maximum meshing-in impact can be expressed as [29]where Jp and are the moments of inertia of pinion and gear, respectively. N is the deformation coefficient under the static state, and ks is single meshing stiffness at the engagement point.

The meshing-in impact of each gear pair is expressed as the sum of Fourier series and can be directly introduced into dynamic equation of two-input two-path STTS.

3. Dynamic Model of Two-Input Two-Path STTS

Figure 8 shows the dynamic model of the two-input two-path STTS. In Figure 8, mij and Iij are the mass and moment of inertia of each gear, respectively. m8 and I8 are the mass and moment of inertia of the central output gear Z8. TL1 and TR1 are left- and right-input torques. Tout is the output torque on the central output gear Z8. ki12, ki34, ki35, ki68, and ki78 are the comprehensive meshing stiffness of each gear pair. ci12, ci34, ci35, ci68, and ci78 are the normal meshing damping of each gear pair. kit23, kit46, and kit57 are the torsional stiffness of dual-gear shaft. cit23, cit46, and cit57 are the torsional damping of dual-gear shaft. φij represents torsional freedom of each gear about its Z-axis. φ8 represents torsional freedom of central output gear Z8 about its Z-axis. kix1, kix2, kix4, and kix5 denote the support stiffness of each gear in the X-direction. cix1, cix2, cix4, and cix5 denote the support damping of each gear in the X-direction. kiy1, kiy2, kiy4, and kiy5 indicate the support stiffness of each gear in the Y-direction. ciy1, ciy2, ciy4, and ciy5 indicate the support damping of each gear in the Y-direction. The rest of support stiffness and support damping are not shown in the figure.

Dynamic model of two-input two-path STTS represents a vibration system with 49 degrees of freedom. Its generalized displacement matrix X can be expressed aswhere x, y, and z are the microdisplacements for two-input two-path STTS in each direction and φ is the microrotation angle for two-input two-path STTS in Z-direction.

Normal relative displacement along the meshing line of each gear pair can be expressed aswhere pi12 (t) is the equivalent meshing line displacement of first-stage spiral bevel gear pair, pi34 (t) and pi35 (t) are the equivalent meshing line displacement of second-stage spur gear pair, and pi68 (t) and pi78 (t) are the equivalent meshing line displacement of third-stage DH gear pair. Parameters aix, aiy, and aiz are spiral bevel gear calculation coefficients.

The meshing force between the teeth of each gear pair can be expressed as

3.1. Dynamic Equations of First-Stage Spiral Bevel Gear

Dynamic equations for the left- and right-input first-stage spiral bevel gear can be expressed aswhere Ti1 is the input torque of a spiral bevel gear, Fis12 is the meshing-in impact of spiral bevel gear, and rib1 and rib2 are the equivalent pitch radii of spiral bevel gear, respectively. Parameters Fi12x = aix(Fin12 + Fis12), Fi12y = aiy (Fin12 + Fis12), and Fi12z = aiz (Fin12 + Fis12) are the component forces in each coordinate axis direction of a spiral bevel gear.

3.2. Dynamic Equations of Second-Stage Spur Gear

Dynamic equations for the left- and right-input second-stage spur gear can be expressed aswhere rib3, rib4, and rib5 are the base circle radii of second-stage spur gear, respectively. Fis34 and Fis35 are the meshing-in impact of spur gear pair.

3.3. Dynamic Equations of Third-Stage DH Gear

Dynamic equations for the left- and right-input third-stage DH gear can be expressed aswhere rib6, rib7, and rb8 are the base circle radii of DH gear, respectively. Tout is the output torque, while Fis68 and Fis78 are the meshing-in impact of DH gear pair, respectively.

3.4. Calculation of the Load Sharing Coefficient

The essence of load sharing coefficient (LSC) is to characterize the difference of power split in the two-path STTS, that is, the difference in the transmission load of each gear caused by various factors (manufacturing error, installation error, and vibration). A Runge–Kutta method is used to solve the dynamic equations of the system. It should be noted that in the process of using coordinate transformation to eliminate the rigid body displacement of the system, for the two-input two-path STTS, there is a rigid body displacement in the left engine closed-loop system and the right engine closed-loop system, which need to be eliminated together. Finally, a dynamic model with 48 degrees of freedom can be obtained.

The instantaneous LSC in each tooth frequency cycle of each path of the left- and right-input second stage and third stage is defined as follows:where k1 = 1,…,m1 and k2 = 1,…,m2. Parameters m1 and m2 are meshing gear frequency cycles of second stage and third stage within the system period, respectively. and are maximum meshing forces in the k1 tooth frequency cycle of the second stage, while and are the maximum meshing forces in the k2 tooth frequency cycle of the third stage.

The LSC of the system period of the left- and right-input second stage and third stage is defined as follows:

Therefore, the LSC of the left- and right-input two-path STTS can be expressed as

4. DH Gear Modification Optimization Design

Topological modification includes tooth profile modification and longitudinal modification, both of which are composed of two quadratic parabolas and a straight line, as shown in Figure 9. Parameters y1 and y3 are the maximum modification amount and modification length of the tooth root, respectively. Parameters y2 and y4 are the maximum modification amount and modification length of the tooth top, respectively. y5 is the maximum modification amount at both tooth ends. y6 is the length of the nonmodification area in the tooth direction. H and L are the effective tooth height and tooth length, respectively. s1, s2, s3, and s4 are the location coordinates of the modification points. For topological modification of DH gear (driving gear), left and right helical gear modification methods are the same.

Topological modification curve has to be defined prior to TCA simulation. The general form of the modification curve is as follows:where ξ is the modification amount, s is the modification point coordinate, s∈[s1, s4], Ca and Cq are modification values, and c is the index of the modification curve.

According to the modification curve, the grid node modification value ξ is calculated. Then, smooth modification amount surface is obtained by fitting the node data through the cubic B-spline. The relationship between the rotating projection surface and the theoretical tooth surface is as follows:where r1u, , and are the coordinate components of DH pinion theoretical tooth surface position vector.

The modified tooth surface is constructed by superimposing the theoretical tooth surface and the modification amount surface. Its position vector and normal vector can be expressed aswhere u1 and l1 are the tooth surface parameters, and r1 (u1, l1) and n1 (u1, l1) are the position and normal vectors of the DH pinion theoretical tooth surface, respectively.

The amplitude of SLTE and meshing-in impact are the main factors affecting the vibration and noise of the DH gear. In this paper, the optimization objective is to minimize the synthesis of the two. Parameters (y1, y2, y3, y4, y5, and y6) are employed for optimization. Modification optimization design model is represented as follows:where f10 and f20 are the amplitude of SLTE and meshing-in impact, respectively, of a nonmodified DH gear. f1 and f2 are the amplitude of SLTE and meshing-in impact, respectively, of a modified DH gear. Parameters qmin and qmax represent the lower and the upper bound, respectively, for y1, y2, and y5. hmin and hmax are the lower and the upper bound, respectively, for y3 and y4. lmin and lmax represent the lower and the upper bound, respectively, for y6.

NSGA-II employs fast nondominated sorting algorithm and a crowded distance comparison operator while introducing an elite strategy, which has low computational complexity and high computational efficiency. Hence, it is used in this paper to solve the aforementioned optimization model. The flowchart of topological modification optimization of DH gear is shown in Figure 10.

5. Numerical Analysis

Some specific parameters of two-input two-path STTS are shown in Table 1.

TVMS and meshing-in impact of first-stage spiral bevel gear are shown in Figure 11. The maximum meshing-in impact is equal to 4080.284 N.

TVMS and meshing-in impact of second-stage spur gear are shown in Figure 12. The maximum meshing-in impact is equal to 8046.040 N.

Optimized DH pinion topological modification tooth surface is shown in Figure 13. The optimized modification parameters are y1 = 0.023 mm, y2 = 0.035 mm, y3 = 2.30 mm, y4 = 2.01 mm, y5 = 0.013 mm, and y6 = 8.19 mm.

TVMS and meshing-in impact of third-stage DH gear before and after modification are shown in Figure 14. The maximum meshing-in impact is equal to 3909.481 N before modification and 1084.720 N after modification. The meshing-in impact is reduced by 72.25% after modification. In Figure 14 and the following figures, “non-mod” represents the results of a DH pinion tooth surface with nonmodification, while “mod” represents the results of a DH pinion tooth surface with the topological modification.

5.1. Dynamic Analysis of Two-Input Two-Path STTS

In two-input two-path STTS, the meshing period of the first-stage gear pair is the shortest, while that of the third-stage gear pair is the longest. In order to consider the influence of meshing-in impact of all gear pair levels, the following simulation calculation uses 1/100 of the meshing period of the first-stage gear pair as the step length. The simulation result employs the calculation results of the last 1500 meshing periods of the third-stage gear pair.

When the DH pinion is not modified, the influence of the installation error and eccentricity error of the second-stage spur gear and third-stage DH gear on the left-input two-path STTS LSC is shown in Figure 15. The law of the right-input two-path STTS LSC is similar to the law of the left-input two-path STTS. When a certain error emerges, the remaining errors are nullified and remain unaltered.

According to Figure 15, the LSC of the left-input two-path STTS linearly increases with an increase of the installation error (eccentricity error). In other words, the greater the installation error (eccentricity error), the greater the system’s LSC, and the worse the load sharing performance. The influence of third-stage DH gears installation error (eccentricity error) on LSC is greater than that of the second-stage spur gears. Under the action of gear comprehensive installation error, LSC may be less than or greater than that of gear installation error alone. However, under the action of gear comprehensive eccentricity error, LSC must be greater than that under the action of the gear eccentricity error alone. To summarize, the influence of the eccentricity error on STTS load sharing performance is greater than that of the installation error. In gear design, the eccentricity error of the third-stage DH gears must be strictly controlled.

In the subsequent analysis, errors of second and third-stage gears are selected according to the gear accuracy grade of 5, as shown in Table 2 [41].

Dynamic meshing force of first-stage gear pair of the left-input two-path STTS with and without considering meshing-in impact is shown in Figure 16. Dynamic force of the first-stage gear pair slightly increases when meshing-in impact is accounted for. Since this paper only analyzes the load sharing and dynamic load of the second and third-stage gear pairs of the system, the dynamic characteristics of the first-stage gear pair are not discussed too much.

The maximum meshing force in the tooth frequency cycle of a second and third-stage gear pair of the left-input two-path STTS with and without considering meshing-in impact is shown in Figure 17.

According to the ISO standard for gears, the dynamic load factor is defined as the ratio of dynamic load to the static load. The dynamic load factor of the second-stage spur gear pair of the STTS is increased from 1.41596 without considering meshing-in impact to 1.42443 with considering meshing-in impact, and the dynamic load performance is reduced by 2.00%. The dynamic load factor of the third-stage DH gear pair of the STTS is increased from 1.20645 without considering meshing-in impact to 1.20669 with considering meshing-in impact, and the dynamic load performance is reduced by 0.12%.

The instantaneous LSC of the second-stage spur gear pair and third-stage DH gear pair with and without the consideration of meshing-in impact is shown in Figure 18. According to the LSC definition, its value cannot be lower than 1. In other words, the average system load without error is the load that the path transmission must bear. The LSC of the second-stage spur gear pair decreases from 1.06387 without considering meshing-in impact to 1.06332 with considering meshing-in impact, thus increasing the load sharing performance by 0.86%. The LSC of the third-stage DH gear pair decreases from 1.09669 without considering meshing-in impact to 1.09603 with considering meshing-in impact, thereby increasing the load sharing performance by 0.68%. The LSC of the left-input two-path STTS depends on the LSC of the third-stage, which is equal to LSCL = 1.09603 with consideration of meshing-in impact. When accounting for meshing-in impact, although the dynamic meshing force of the system increases, the increase of denominator in (17) is larger than that of the numerator, which ultimately leads to a slight decrease in the LSC of the STTS. It should be noted that if the increase of the numerator in (17) is greater than the increase of the denominator after the meshing-in impact is considered, the LSC of the system will increase. However, after considering the meshing-in impact, the dynamic meshing force of the system will inevitably increase, and the dynamic load characteristics will inevitably decrease, which is unfavorable to the transmission characteristics of the entire system.

5.2. Influence of Modification on Load Sharing and Dynamic Load Characteristics

In this paper, the design load of DH gear pair is equal to 17500 Nm, and the theoretical contact ratio of DH gear pair is 3.42. Figure 19 shows the waveform change of SLTE under multiple loads for DH gear pair before and after modification. SLTE is illustrated three times by shifting the angular pitch. The rightmost column of Figures 19(a) and 19(b) shows the amplitude of SLTE of DH gear pair.

Figure 19(a) shows the case of the DH gear before modification. For DH gear with a standard tooth surface, the TE is approximately zero. In the process from entering meshing to exiting meshing, the meshing of quadruple tooth pairs and the meshing of triple tooth pairs alternately change. The SLTE of triple tooth pairs region is greater than that of quadruple tooth pairs. Because there is no STE to compensate for the difference of SLTE in different meshing zones, as the load increases, the change of the SLTE waveform gradually increases. In other words, the amplitude of SLTE increases linearly with the increase of the load, as shown by the black curve in Figure 19(c).

Figure 19(b) shows the case of DH gear after topological modification. When load is increased from the initial value of 1500 Nm to the design load, the amplitude of SLTE increases first and then decreases. This is due to the reduction of the contact ratio of the gear pair caused by the modification. When the load is small, there is an alternation of triple tooth pairs meshing and double tooth pairs meshing in the process of the tooth from meshing-in to meshing-out. When the load exceeds the design load, the amplitude of SLTE continues to gradually increase. This is because with the increase of the load, the actual contact ratio of the gear pair gradually increases and finally reaches the theoretical contact ratio. In the process of tooth entering meshing to exiting meshing, the quadruple tooth pairs meshing and triple tooth pairs meshing alternate, which is the same as the meshing situation of the DH gear before modification. However, because the STE after modification compensates the difference of SLTE in different meshing areas, the amplitude of SLTE after modification is still much smaller than that before modification. Before and after modification, the amplitude of SLTE changes from 2.495″ to 0.220″, respectively; i.e., it decreased by 91.18%, as shown in Figure 19(c).

The comprehensive meshing stiffness for DH gear pair under multiple loads before and after modification is shown in Figure 20. Similarly, comprehensive meshing stiffness is illustrated three times by shifting the angular pitch. The rightmost column in Figures 20(a) and 20(b) indicates the average value of the meshing stiffness of DH gear pair.

Corresponding to Figures 19(a), and 20(a) represents the comprehensive meshing stiffness of the DH gear before modification. Since the DH gear has a standard involute tooth surface, the two tooth surfaces of the pinion and gear are completely conjugated under load. Furthermore, each meshing moment is characterized by a linear contact. Thus, with an increase in load, the average meshing stiffness increases slightly.

Corresponding to Figures 19(b), and 20(b) represents the comprehensive meshing stiffness of the DH gear after modification. The modification destroys the conjugate characteristics of the standard tooth surface, while pinion and gear tooth surface changes from a line contact to the point contact. The length of the contact line and the area of the contact ellipse of the modified tooth surface gradually increase with an increase in load. When load is increased from 1500 Nm to the design load, the average meshing stiffness of the modified DH gear pair rapidly increases. When subjected to the design load, the average meshing stiffness decreases from 1.24269 × 109 N/m before modification to 1.11100 × 109 N/m after modification; i.e., it decreases by 10.60%. When the applied load exceeds the design load, the actual contact ratio of the gear pair is close to the theoretical contact ratio. The average meshing stiffness slowly increases, and it approaches the average meshing stiffness before modification, as shown in Figure 20(c).

The research results in Figures 19 and 20 show that, on the one hand, modification reduces the amplitude of SLTE and the fluctuation amplitude of comprehensive meshing stiffness of the DH gear pair. On the other hand, modification leads to the decrease of the average meshing stiffness. The reduction of the fluctuation amplitude of comprehensive meshing stiffness is beneficial to the STTS’s dynamic load characteristics, while the reduction of the average meshing stiffness will have an impact on the STTS’s load sharing characteristics, which will be analyzed below.

The maximum meshing force in the tooth frequency cycle of second and third-stage gear pair of left-input two-path STTS before and after modification is shown in Figure 21. Only the fractional part of the dynamic load coefficient greater than 1 (indicating vibration) can be reduced via gear modification and dynamic design to improve the dynamic performance of gear transmission and reduce vibration level. Dynamic load factor of the STTS second-stage is reduced from 1.42443 before modification to 1.40194 after modification. Thus, dynamic load performance is improved by 5.30%. However, third-stage dynamic load factor is significantly altered. The dynamic load factor is reduced from 1.20669 before modification to 1.13454 after modification, thereby increasing the dynamic load performance by 34.91%. The above research shows that the modification has a good effect on improving the dynamic load characteristics of the STTS third stage.

The instantaneous LSC of the second-stage spur gear pair and third-stage DH gear pair before and after modification is shown in Figure 22. Only the fractional part of the LSC greater than 1 (indicating vibration) can be reduced to improve the load sharing performance of the gear transmission through modification. The LSC of the second-stage spur gear pair decreased from 1.06332 before modification to 1.06183 after modification, hence increasing the load sharing performance by 2.35%. The LSC of the third-stage DH gear pair decreased from 1.09603 before modification to 1.09203 after modification, therefore increasing the load sharing performance by 4.17%. The LSC of the left-input two-path STTS depends on the LSC of the third stage. After modification, it is equal to LSCL = 1.09203. The results indicate that DH gear modification has a certain effect on improving the system’s load sharing performance.

6. Conclusions

In this paper, the reversal method is used to calculate the meshing-in impact of standard involute spur gears and double-helical gears, and the meshing-in impact of spiral bevel gears and modified double-helical gears is accurately calculated based on tooth contact analysis and loaded tooth contact analysis. A dynamic model of two-input two-path split torque transmission system considering meshing error excitation, time-varying meshing stiffness excitation, and meshing-in impact excitation is proposed. Influence of double-helical gear modification on load sharing and dynamic load characteristics of two-input two-path split torque transmission system is investigated. The main conclusions are listed as follows:(1)After considering the meshing-in impact, the dynamic meshing force of the gear pair increases, which is unfavorable to the transmission of the split torque transmission system. However, due to an increase in the average value of the maximum meshing force in each meshing period of the two paths, the load sharing coefficient of the system slightly decreases.(2)After double-helical gear modification, the amplitude of static loaded transmission error, the fluctuation amplitude of comprehensive meshing stiffness, and the meshing-in impact are all reduced. Dynamic load characteristics of the third-stage gear pair are significantly improved, but dynamic load characteristics of the second-stage gear pair are less affected by double-helical gear modification.(3)After double-helical gear modification, the average meshing stiffness is reduced. Simultaneously, due to an improvement of the split torque transmission system dynamic load characteristics, the load sharing characteristics of the system are also improved to a certain extent.

Abbreviations

zij (i = L, R, j = 1…7):Gears
z8:Center output double-helical gear
σ1, σ2:System installation angle
ψ:System layout angle
γ:The angle between the center line of each gear pair and positive X-axis direction
η:The angle between the meshing line of each gear pair and the positive direction of the X-axis
EL3, EL4:Eccentricity error of the spur pinion zL3 and the spur gear zL4
AL3, AL4:Installation error of the spur pinion zL3 and the spur gear zL4
λL3, λL4, μL3, μL4:Phase angle
MLL34:The meshing line
ei34, ei35:Equivalent cumulative meshing error on the meshing line of spur gear pair 34 and 35
ei68, ei78:Equivalent cumulative meshing error on the meshing line of double-helical gear pair 68 and 78
ωi3, ωi4, ωi5:Angular velocities of spur gear zi3, zi4, and zi5
ωi6, ωi7, ω8:Angular velocities of double-helical gear zi6, zi7, and z8
αts, αth:Transverse pressure angle of spur gear and double-helical gear
P:Total normal load
Z:Normal deformation
LTE:Static loaded transmission error
ALTE:Amplitude of static loaded transmission error
k:Comprehensive meshing stiffness
Δφkg, Δφkp:Reverse angles of the driven gear and driving pinion
rp, :Pitch circle radii of pinion and gear
rbp,rbg:Base circle radii of pinion and gear
rag:Addendum circle radius of the gear
Δφ:Backward angle of the driven gear
:Normal velocities of pinion and gear
:Relative normal velocity
Fs:Meshing-in impact
:Moments of inertia for pinion and gear
n:Deformation coefficient
ks:Single meshing stiffness at the engagement point
mij, Iij:Mass and moment of inertia of gear zij
m8, I8:Mass and moment of inertia of the central output gear Z8
TL1, TR1:Left- and right-input torques
Tout:Output torque on the central output gear Z8
ki12, ki34, ki35, ki68, ki78:Comprehensive meshing stiffness of each gear pair
ci12, ci34, ci35, ci68, ci78:Normal meshing damping of each gear pair
kit23, kit46, kit57:Torsional stiffness of dual-gear shaft
cit23, cit46, cit57:Torsional damping of dual-gear shaft
φij:Torsional freedom of each gear about their Z-axis
φ8:Torsional freedom of central output gear Z8 about its Z-axis
kix1, kix2, kix4, kix5:Support stiffness of each gear in the X-direction
cix1, cix2, cix4, cix5:Support damping of each gear in the X-direction
kiy1, kiy2, kiy4, kiy5:Support stiffness of each gear in the Y-direction
ciy1, ciy2, ciy4,ciy5:Support damping of each gear in the Y-direction
pi12 (t):Equivalent meshing line displacement of spiral bevel gear pair
pi34 (t), pi35 (t):Equivalent meshing line displacement of spur gear pair
pi68 (t), pi78 (t):Equivalent meshing line displacement of double-helical gear pair
aix, aiy, aiz:Spiral bevel gear calculation coefficients
Fin12, Fin34, Fin35, Fin68, Fin78:Meshing force of each gear pair
Fis12:Meshing-in impact of spiral bevel gear
rib1, rib2:Equivalent pitch radii of spiral bevel gear
rib3, rib4, rib5:Base circle radii of second-stage spur gear
Fis34, Fis35:Meshing-in impact of spur gear pair
rib6, rib7, rb8:Base circle radii of double-helical gear
Fis68, Fis78:Meshing-in impact of double-helical gear pair
y1, y2, y3, y4, y5, y6:Modification parameters
s1, s2, s3, s4:Location coordinates of the modification points
Ca, Cq:Modification values
u1, l1:Tooth surface parameters
r1 (u1, l1), n1 (u1, l1):Position and normal vectors of the double-helical pinion theoretical tooth surface
qmin, qmax:Lower and upper bound for y1, y2, and y5
hmin, hmax:Lower and upper bound for y3 and y4
lmin, lmax:Lower and upper bound for y6
NSGA-II:Nondominated sorting genetic algorithm-II
STTS:Split torque transmission system
DH gear:Double-helical gear
TVMS:Time-varying meshing stiffness
TCA:Tooth contact analysis
LTCA:Loaded tooth contact analysis
STE:Static transmission error
SLTE:Static loaded transmission error
LSC:Load sharing coefficient
non-mod:Nonmodification
mod:Modification.

Data Availability

All the data supporting this paper are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was financially supported by the National Natural Science Foundation of China (project no. 51675424).