#### Abstract

In order to solve the dynamic vibration characteristics of the power-split transmission system, the system of the dynamic mechanical model is established. Firstly, according to the theoretical analysis method of the tooth contact analysis (TCA) and loaded tooth contact analysis (LTCA), the actual meshing process of each gear pair is simulated, and the time-varying mesh stiffness excitation is obtained, which can improve the numerical precision. Next, by using the lumped mass method, the bending-torsional coupling three-dimensional dynamical model of the power-split transmission is established. The identical dimensionless equations are deduced by eliminating the effect of rigid displacement and the method of dimensional normalization. Next, the frequency domain and time domain responses of this system are obtained. The dynamic load change characteristics of each gear pair are analyzed. The results show that establishment, solution, and analysis of the system dynamics model could provide a basis for the dynamic design and have an important significance for the dynamic efficiency analysis and dynamic performance optimization design of the power-split transmission. Through theoretical data compared with the experimental data, we verified the correctness of the method proposed.

#### 1. Introduction

A power-split transmission adopts the power-split technology; it can realize the power split in four ways. Due to compact structure of power-split transmission, it is designed to meet the operating conditions with high speed and overloaded operation. It has a broad application prospect in the aerospace and marine industry. The vibration and dynamic of power-split transmission system are of important in the application; it is directly related to the transmission system reliability and security.

Many researchers [1–4] had already done lots of analysis on the load-sharing systems with elastic supports at abroad and home. Kish [5, 6] has reported on the development and testing of a power-split transmission that featured a torsional elastomeric load-sharing device. A companion study to develop a method to analyze and optimize the load sharing of split-path gear boxes has also been completed; the results of that study are reported separately by Krantz et al. [7–9]. The effect of time-varying mesh stiffness, friction static transmission error, and elastic supporting has been considered. Dynamic load-sharing behavior, the load-sharing coefficient of star gear trains with effect of each level of connection stiffness, and star gear eccentric errors have been analyzed by Guo et al. [10]. Some other scholars [11–19] have also conducted researches on the power-split transmission.

However, much of recent research on power-split transmission systems minimized the influence of gear surface tooth contact; most of these researches just adopted equivalent average mesh stiffness of a gear pair to express progress of load research and could not accurately reflect a real meshing process. In this paper, a real meshing process of gear pairs is dispersed into a limited meshing point, according to the theoretical analysis method of tooth contact analysis (TCA) and loaded tooth contact analysis (LTCA). Statics characteristic of each meshing position is analyzed. The mechanical properties are obtained. This method can improve accuracy of the calculation.

In this paper, in order to obtain good dynamic performance of double road power-split transmission systems, the lumped parameter method is used to establish the bend-torsion coupling dynamic model. Then, the differential equation is solved to obtain the frequency domain and time domain response characteristics and the dynamic load coefficient of the system. Through the analysis and solution, it will help to further understand the dynamic characteristics of the system. It will provide a theoretical basis for the design and manufacturing of dual power-split transmission. It will have a great significance to the reasonable design and successful use for the dual power-split transmission system.

#### 2. Dynamics Model Building

Figure 1 is a schematic diagram of the three-dimensional structure of power-split transmission system. The system adopts power-split technology. The power is from the input shaft 1. Level I adopts helical gears. Level II adopts spur gear. Two torsional axes are used to link the two levels of transmission. The two gears 5 of level II simultaneously mesh big gear 6 of level II. This transmission system can further effectively solve the drive big torque and speed ratio.

Lumped parameter method is adopted to establish the dynamic model of bend-torsion coupling. Dynamic model can be expressed in Figure 2. Here, is input torque; is output torque; is time-varying mesh stiffness; is meshing damping; and are torsion angle; is integrated error; is moment of inertia concentrated mass; and are torsional rigidity; and are torsional damper.

The meshing force can be represented by

The damping force can be represented bywhere and are displacement deformation along the -axis; and are displacement deformation along the -axis; and are the amplitudes of errors along the -axis; and are the amplitudes of errors along the -axis, respectively, for the pinion and gear ; is actual operating pressure positive angle; is the base radius.

The dynamics differential equation of lateral-bending vibration can be represented bywhere and are the equivalent supporting rigidity; and are the equivalent bearing damping; , and are the transverse vibration acceleration, velocity, and displacement along -direction; , and are the transverse vibration acceleration, velocity, and displacement along -direction.

The angular displacement of generalized coordinates is translated into the line displacement. The differential equations of torsional direction can be represented bywhere is equivalent mass. The variable-step four-order Runge-Kutta method is used to solve (2) and (3). is the definition of dimensional time. is displacement nominal scale.

The equations (2) and (3) are processed by the dimensional normalization.

Due to the influence of the initial value, the initial transient response cycles are deleted. Then, vibration responses of displacement, velocity, and acceleration are calculated. By using the fast Fourier transform (FFT) method, the frequency responses of displacement, velocity, and acceleration are derived. Finally, the dynamic load coefficient of the system is calculated.

The results of dynamic response characteristics are generated into (1) to get dynamic load. Dynamic load coefficient can be expressed by

The dynamic load coefficient directly reflects the dynamic response of the system. The vibration of the system is more obvious as the dynamic load coefficient get bigger.

#### 3. Time-Varying Mesh Stiffness Based on the LTCA

Under the external load, the torsional angle of tooth will be changed. The change of torsional angle consists of geometry transmission errors, bending deformations, and contact deformations of tooth [15]. Functional relations between and are expressed by where , , and are the deformations of torsional angle; is the torque of the th meshing position in a meshing cycle; , , and are constants.

The LTCA (load tooth contact analysis) models are shown in Figure 3, where the two pairs of teeth which contacted each other at a specific moment in the meshing cycle are denoted by I and II. The tooth surface curve is vertical along the relative principal direction in the normal plane, which is shown in Figure 3. is the contact point; is a point along the relative principal direction [20].

Under the load , the driving gear goes through an approach . Due to the tooth deformation, contact load will become distributed. After contact deformation, the state can be described by the following equation for the tooth pair :where ; ; ; ; is the normal flexibility coefficient matrix of the gear pair ; is the contact load supported at point of the tooth pair ; is the final tooth clearance at point ; is the tooth approach. Apparently, the contact force must satisfy the following conditions: where if , then ; if , then . , and are known. The contact forces , final tooth clearance , and tooth approach are unknown.

The known parameters and unknown parameters constitute a nonlinear program model. According to the tooth approach , the following objective function is established by

Equations (8) and (9) represent a constrained nonlinear programming problem, which is solved by the modified simplex method.

The objective function (9) forms a nonlinear programming model with functions (7) and (8) as constraint conditions [14]: where is the artificial variables; of each element is equal to 1.

The corresponding angular error is determined byHere, is the helix angle.

The load distribution on the contact lines of the tooth surface is shown in Figure 4.

**(a)**

**(b)**

The tooth approach solved from the nonlinear programming problem for each contact position is actually the loaded tooth transmission errors as the amount of linear displacement error of the driven gear along the contact normal (the line of action). The corresponding angular transmission error under load for the contact position is determined by reversing (5). The column vector is solved from the programming problem that represents the discrete distribution of the contact load along the contact line that coincides with the relative principal direction.

By solving (5), we can obtain the coefficient of , , and . Then, functional relations between loaded transmission errors and some nominal load of may be proposed. The calculation curves are supplied in a meshing cycle. The time-varying mesh stiffness is represented by

The mesh stiffness could reflect real meshing elastic properties at the meshing position more directly. The discrete value of meshing stiffness is fitted by the polynomial and through the Fourier series transformation to spread out into a periodic function.

#### 4. Spline Clearance Floating

In order to improve the uniform load distribution of the power-split transmission system and to solve the problem that elastic torsion shaft cannot completely satisfy the load-sharing characteristics, a structure of level pinion floating is proposed.

Level floating pinion is installed on one end of input shaft with high speed characteristics. It is connected with output components through a short spline. The spline can transmit the torque. However, floating pinion cannot completely float freely under the constraint of spline coupling. The supporting rigidity of floating pinion can be described in Figure 5.

In the process of the spline transmit torque, friction will be produced between internal and external spline. is friction; here, is the positive pressure between internal and external spline and is friction coefficient. The floating quantum can be represented byHere, and are the floating quantum along the -direction and -direction, respectively; is the iterations.

The floating pinion is affected by both of the engaging force of the two associated gears and support reaction of spline coupling. It can be seen from Figure 5 that when the support reaction is less than the friction, the internal and external spline cannot produce a slippage. Here, the bending deflection of input shaft will adapt to the change of position of floating pinion (0 to ). When the support reaction is greater than the friction, the internal and external spline will produce a slippage. Here, the slippage will adapt to the change of positions of floating pinion ( to ). However, if the slippage is beyond , namely, radial clearance between internal and external spline is eliminated, the bending deflection of input shaft will again adapt to the change of position of floating pinion. represent the radial clearance between the internal and external spline.

and represent the support reaction of floating pinion projected on the -axis and -axis, respectively:Here, is the flexural rigidity of spline shaft; is a direction angle of vector of .

and represent the damping force of floating pinion projected on the -axis and -axis, respectively:Here, is the damping of spline shaft bending and is frictional damping.

The support equilibrium conditions of floating pinion can be represented by

Equations (16) will be combined with other dynamics differential equations to establish clearance nonlinear mathematical model. Then, through solving this nonlinear mathematical model, the transmission torque of each gear pair is obtained. Finally, the load-sharing coefficient will be obtained.

#### 5. Examples

To verify the correctness of the method proposed in this paper, the initial parameters of the numerical examples will be completely based on the initial data in [8] of the NASA Research Institutions. Finally, the theoretical analysis results proposed in this paper and experimental results of [8] will be compared.

Gear parameters are shown in Table 1. Input power and input speed r/min.

The loaded transmission errors (LTE) of each gear pair of the system related to Table 1 are shown in Figure 6.

**(a)**

**(b)**

Figure 6 shows the load transmission error (LTE) of different meshing positions on the tooth face. Figure 6(a) is the LTE of level I helical gear pairs. Here, the LTE fluctuation is ()/′′. Figure 6(b) is the LTE of level II gear pair. Here, the LTE fluctuation is ()/′′. Considering the above analysis, the range of fluctuation of the spur gear pairs LTE is larger than the one of the helical gear pairs.

The time-varying mesh stiffness incentive of each gear pair abased on the LTCA is shown in Figure 7; meanwhile, the time-varying mesh stiffness incentive based on the empirical formula (IS06336.1996) [19] is shown.

**(a)**

**(b)**

Figure 7 shows the mechanical properties of different meshing positions on the tooth face. Figure 7(a) is the time-varying mesh stiffness incentive of the first-stage gear pair. Here, the meshing stiffness fluctuation is ()/N·mm^{−1}. Figure 7(b) is the time-varying mesh stiffness incentive of the second-stage gear pair. Here, the meshing stiffness fluctuation is ()/N·mm^{−1}.

Differential equation is expressed in formula of . The natural characteristics of problem can be converted into solving the equations of , is system mass matrix, is stiffness matrix, and the time-varying meshing stiffness is replaced with average mesh stiffness. 30 order natural frequency results are calculated and shown in Table 2.

As seen from Table 2, the fundamental frequency of system is 142.802 Hz. Then, due to the symmetric structure of dual power-split transmission system, the natural frequency of the equal value appeared. Due to the complexity and particularity of the system, the quadruple frequency phenomenon appeared.

The natural frequency corresponding to the various degrees of freedom vibration model is shown in Figure 8. It can be seen from the vibration mode curve, the node between 23 and 29 of the fundamental frequency is the largest relative amplitude.

From the above analysis, the torsional-vibration natural frequency of transmission system has a reasonable distribution. Meanwhile, it can meet the requirements of design.

The dynamic load reflects the vibration amplitude. Figure 9 is time-history dynamic load and FFT spectrum diagram of the gears at all levels of system.

**(a) The I-stage dynamic load and the corresponding FFT spectrum**

**(b) Level II dynamic load and the corresponding FFT spectrum diagram**

It can be seen from Figure 9 that level I of meshing frequency is 4682.67 Hz and level II of meshing frequency is 1019.13 Hz. After a certain calculation, dynamic load coefficient of the gear pair at all levels of system is obtained in the standard condition without errors, respectively, which is 1.152 and 1.065. At the same time, it can be seen that a strong vibration of the whole system happened in meshing frequency of the gear pair.

Vibration acceleration reflects the impact force. The system vibration acceleration time history is shown in Figure 10 with the standard conditions.

**(a) Level I acceleration and the corresponding FFT spectrum**

**(b) Level II acceleration and the corresponding FFT spectrum**

It can be seen from Figure 10 that the vibration acceleration RMS (m/s^{2}) of the gear pair at all levels of system is 4.989 and 2.163, respectively. It can be seen that a strong vibration of the vibration acceleration happened in meshing frequency of the gear pair at all levels.

When the input speed was given 4390 r/min, the vibration acceleration and FFT spectrum of the II grade gear pair are shown in Figure 11.

It can be seen from Figure 11 that when the input speed is reduced to 4390 r/min, the vibration acceleration RMS is 1.290 and 1.290, respectively. Here, the biggest vibration amplitude is in the fourth harmonic.

With the change of rotating speed, dynamic load coefficient is shown in Figure 12; here, is meshing damping ratio.

**(a) Level I dynamic load coefficient**

**(b) Level II dynamic load coefficient**

It can be seen from Figure 12 that the dynamic load coefficient is changed with the change of rotating speed. By the above analysis of system dynamic characteristic, transmission system must stay away from the resonance speed according to the reasonable working speed requirements. With the increase of damping ratio, the dynamic load of near critical speed can be reduced.

The center distance installation errors influence the load-sharing characteristics of transmission system. Such as when mm is given, the load-sharing coefficient is 1.0983. If these errors individually affect the load-sharing characteristics of transmission system, the result is shown in Figure 13.

**(a)**Gear 2 error of

**(b)**Gear 4 error of

**(c)**Gear 6 error ofFigure 13 shows that the torque is cyclically fluctuating at each meshing position in different errors of , and , which reflect the load distribution at different engagement position in the tooth surface. With the influence of , and , the load-sharing coefficient is, respectively, 1.0207, 1.0783, and 1.0641.

Load-sharing coefficient with a single influence of the center distance installation error is shown in Figure 14.

Figure 14 shows that level II pinion plays a vital role in the load-sharing coefficient. Thus, in the process of the system installation, level II pinion errors should be mainly considered.

The influence of the floating pinion based on spline clearance floating is shown in Figure 15.

**(a) Spline clearance is 0 mm**

**(b) Spline clearance is 0.85 mm**

**(c) Spline clearance is 1.231 mm**

Figure 15 shows that load-sharing coefficient is 1.0307, 1.0133, and 1.0044 under the influence of the spline clearance.

Under the different error conditions, the curves of load-sharing coefficient changed with the spline clearance are shown in Figure 16.

**(a)**Clearance floating coefficient changes with

**(b)**Clearance floating coefficient changes withIt can be seen from Figure 16 that the load-sharing coefficient is decreased with the increase of spline clearance. In other words, load-sharing performance of system will be better.

#### 6. Dynamics Model Validation

In this section, the theoretical calculation based on the method proposed in this paper and the analysis results of [8] will be compared. Dual power-split transmission system mentioned in [8] is a power transmission device used in aviation helicopters. Experimental equipment of dual power-split transmission system is given in [8]; this is shown in Figure 17(a). Dual power-split transmission device is composed of a power flow of closed loop system. The principle of system is shown in Figure 17(b).

**(a) Experimental split torque gearbox test facility**

**(b) Power closed circuit principle diagram**

Theoretical analysis results of the time-history dynamic load of the gear pair at all levels of system are compared with the calculation results of [8]. Then, the results are shown in Figure 18. Here, the condition is the standard no errors conditions; the data is the time domain steady state response data. At the same time, in the process of transmission torque, the static loads of the various gears are given, respectively, which are 4.05 × 10^{5} N·mm, 7.86 × 10^{5} N·mm, and 1.02 × 10^{7} N·mm. In order to make a better comparison, the unit of inch-pounds (in.-lb) of [8] was transformed to the N·mm.

**(a) Level I dynamic load from [8] analysis**

**(b) Level I dynamic load from the theory analysis**

**(c) Level II dynamic load from [8] analysis**

**(d) Level II dynamic load from the theory analysis**

**(e) Level III dynamic load from [8] analysis**

**(f) Level III dynamic load from the theory analysis**

It can be seen from Figure 18(a) that dynamic load of level I gear pair is fluctuated 4.05 × 10^{5} N·mm from [8]. Figure 18(b) shows that the dynamic load of level I gear pair is fluctuated 4.11 × 10^{5} N·mm from theoretical calculation and analysis. Here, dynamic load coefficient is 1.152 and 1.117, respectively. It can be seen from Figure 18(c) that dynamic load of level II gear pair is fluctuated 7.86 × 10^{5} N·mm from [8]. Figure 18(d) shows that the dynamic load of level II gear pair is fluctuated 7.92 × 10^{5} N·mm from theoretical calculation and analysis. Here, dynamic load coefficient is 1.065 and 1.112, respectively. It can be seen from Figure 18(e) that dynamic load of level III gear pair is fluctuated 1.02 × 10^{6} N·mm from [8]. Figure 18(f) shows that the dynamic load of level III gear pair is fluctuated 1.14 × 10^{6} N·mm from theoretical calculation and analysis. Here, dynamic load coefficient is 1.062 and 1.077, respectively. Through the above analysis, the theoretical analysis and the calculation data of [8] are consistent. The dynamic load coefficient is also consistent. It is effective to prove the feasibility of dynamic calculation method proposed in this paper.

Figure 19(a) is the FFT spectrum of level II gear pair dynamic load from [8]; Figure 19(b) is the FFT spectrum of level II gear dynamic load from the theoretical analysis.

**(a) FFT spectrum diagram from [8] analysis**

**(b) FFT spectrum diagram from the theory analysis**

It can be seen from Figure 19(a) that the meshing fundamental frequency is 1034.67 Hz from [8]. Figure 19(b) shows that the meshing fundamental frequency is 1019.13 Hz from theoretical calculation and analysis. Through the above analysis, trend of the theoretical analysis and the calculation data of [8] are consistent.

The influence of installation error (0.0040 in. (0.1016 mm)) is considered in [8]. The results of two branches of time-history dynamic load are shown in Figure 20(a). Then, by using the calculation method proposed in this paper, the theoretical calculation results are given in Figure 20(b).

**(a) Dynamic load from [8] analysis**

**(b) Dynamic load from the theory analysis**

The power allocation is, respectively, 56.59% and 55.91% with the analysis of [8] and theoretical analysis in this paper. The loading-sharing coefficient is, namely, 1.132 and 1.118. The results of [8] and theoretical analysis are consistent.

The equivalent average meshing stiffness is adopted in [8], and it cannot reflect the real tooth surface meshing process. The mathematical model putted forward in this paper is based on the tooth contact analysis (TCA) and load tooth contact analysis (LTCA). This method can be dispersed in the whole meshing process of each gear pair to limited meshing point. Each meshing position will be analyzed by the mechanical properties.

Time-varying mesh stiffness using the above analysis method has a higher accuracy than the traditional stiffness excitation using average mesh stiffness. This mesh stiffness more accurately reflects the dynamic characteristics of the whole system. Mathematical model of dual power-split transmission system proposed in this paper will have some theoretical value; at the same time, this mathematical model could be applied in the planetary and star gear transmission system.

#### 7. Conclusions

Through the analysis and comparison, we can get the following main conclusions.(1)It can get more accurate mesh stiffness using TCA and LTCA technology. The dynamic model is established through an effective integration of power-split transmission system.(2)Time-varying mesh stiffness based on LTCA was analyzed by using the Fourier series transfer; the periodic function was used to improve the calculation accuracy. The dynamics model of power-split transmission system can provide a basis for solving and analyzing the dynamic design.(3)The installation errors accumulatively influence the load-sharing characteristics. The installation errors of level II components should be paid more attention. The floating pinion can effectively improve the load-sharing characteristics of system. The quantity of spline clearance should not be excessive. Too much clearance will lead the system to produce serious vibration and shock.

#### Conflict of Interests

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

#### Acknowledgments

The authors would like to acknowledge the support provided to the project by the National Natural Science Foundation of China (51375384) and the Research Foundation of Education Bureau of Shaanxi Province, China (12jk0692).