# Dynamic Vibration Characteristic Analysis for the Power-Split Transmission System Based on Loaded Tooth Contact Analysis

Authors: Hao Dong; Yan Cao; Zhou Fang
Journal: Shock and Vibration
 (2015)

---

## 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.


---

## Body

## 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

Figure1 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.Figure 1
Example of a power-split transmission design with power paths. Illustration of(1) input shaft, (2) I-stage helical pinion, (3) I-stage helical gear, (4) torsion shaft, (5) II-stage spur pinion, (6) II-stage spur gear, and (7) II-stage spur pinion.Lumped parameter method is adopted to establish the dynamic model of bend-torsion coupling. Dynamic model can be expressed in Figure2. Here, Tin is input torque; Tout is output torque; Kij(i=1,4,5;j=2,3,6) is time-varying mesh stiffness; Cij is meshing damping; φi and φj are torsion angle; etij is integrated error; Mi is moment of inertia concentrated mass; k24 and k35 are torsional rigidity; c24 and c35 are torsional damper.Figure 2
Schematic of mechanical structure model.The meshing force can be represented by(1)Pij=Kijrbiφi-rbjφj+xi-xjcosζij+yi-yjsinζijkkkk+ΔAxi-ΔAxjcosζijkkkk+ΔAyi-ΔAyjsinζij-eijt.The damping force can be represented by(2)Dij=Cijrbiφ˙i-rbjφ˙j+x˙i-x˙jcosζij+y˙i-y˙jsinζijkkkkkk+ΔAx˙i-ΔAx˙jcosζijkkkkkk+ΔAy˙i-ΔAy˙jsinζij-e˙ijt,where xi and xj are displacement deformation along the x-axis; yi and yj are displacement deformation along the y-axis; ΔAxi and ΔAxj are the amplitudes of errors along the x-axis; ΔAyi and ΔAyj are the amplitudes of errors along the y-axis, respectively, for the pinion i and gear j; ζij is actual operating pressure positive angle; rbi is the base radius.The dynamics differential equation of lateral-bending vibration can be represented by(3)∑Pij+Dij·cosζij+Mix¨i+Cxix˙i+Kxixi=0,∑Pij+Dij·sinζij+Miy¨i+Cyiy˙i+Kyiyi=0,where Kxi and Kyi are the equivalent supporting rigidity; Cxi and Cyi are the equivalent bearing damping; x¨i,x˙i, and xi are the transverse vibration acceleration, velocity, and displacement along x-direction; y¨i,y˙i, and yi are the transverse vibration acceleration, velocity, and displacement along y-direction.The angular displacement(ui=riφi) of generalized coordinates is translated into the line displacement. The differential equations of torsional direction can be represented by(4)meq,1u¨1+P12+D12+P13+D13=Tinr1cos⁡an,meq,2u¨2-P12+D12+c24rb2u˙2rb2-u˙4rb4+k24rb2u2rb2-u4rb4=0,meq,3u¨3-P13+D13+c35rb3u˙3rb3-u˙5rb5+k35rb3u3rb3-u5rb5=0,meq,4u¨4+P46+D46-c24rb4u˙2rb2-u˙4rb4-k24rb4u2rb2-u4rb4=0,meq,5u¨5+P56+D56-c35rb5u˙3rb3-u˙5rb5-k35rb5u3rb3-u5rb5=0,meq,6u¨6-P46+D46-P56+D56=-Toutrb6,where meq,i(meq,i=Ii/rbi2) is equivalent mass. The variable-step four-order Runge-Kutta method is used to solve (2) and (3). t=τwn is the definition of dimensional time. bc 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 (KD) can be expressed by(5)KD46KD56T=maxP46maxP56TTin/2i12.The dynamic load coefficient(KD) 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 Δφij(Tij(k)) and Tij(k) are expressed by (6)ΔφijTijk=δ1Tk+δ2Tk+δ3Tk=a+bTijk+cTijk2/3,where δ1, δ2, and δ3 are the deformations of torsional angle; Tij(k)(k=1,2,…,5) is the torque of the kth meshing position in a meshing cycle; a, b, and c are constants.The LTCA (load tooth contact analysis) models are shown in Figure3, 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. ik(k=I,II) is the contact point; jk is a point along the relative principal direction [20].Figure 3
Load tooth contact analysis model.Under the loadP, the driving gear goes through an approach Z. 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 k:(7)fkpjk+wk=Z+djk,where k=I,II; [pj]k=[p1,p2,…,pjk,…,pnk]T; [dj]k=[d1,d2,…,djk,…,dnk]T; [Z]=Zk[1,1,…,1,…,1]T; [f]k is the normal flexibility coefficient matrix of the gear pair k; pj(j=1,2,…,nk) is the contact load supported at point j of the tooth pair k; dj(j=1,2,…,nk) is the final tooth clearance at point j; Zk is the tooth approach. Apparently, the contact force pj(j=1,2,…,nk) must satisfy the following conditions: (8)∑j=1npjΙ+∑j=1npjΙΙ=P,where if [pj]k>0, then [d]jk=0; if [pj]k=0, then djk>0. f,P, and w are known. The contact forces [pj]k, final tooth clearance [d]jk, and tooth approach [Z] are unknown.The known parameters(f,P,w) and unknown parameters (p,d,Z) constitute a nonlinear program model. According to the tooth approach Zk, the following objective function is established by(9)MinimizepTFp+w-Z.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]: (10)Min∑j=12n+1Xjkkkkk-Fp+Z+d+X=wS.t.keTp+X2n+1=Pkkkkkkpj≥0,dj≥0,Z≥0,Xj≥0;kkkkkkpj=0ordj=0,where Xjj=1,2,…,2n+1 is the artificial variables; e of each element is equal to 1.The corresponding angular errorΔφij(Tij(k)) is determined by(11)ΔφijTijk=ΔδijTijkrbi·cos⁡ξi.Here, ξi is the helix angle.The load distribution on the contact lines of the tooth surface is shown in Figure4.Figure 4
Load distribution on the tooth surface. Illustration of (a) I-stage helical gear pairs and (b) II-stage spur gear pairs.
(a)
(b)The tooth approachZ solved from the nonlinear programming problem for each contact position is actually the loaded tooth transmission errors as the amount of linear displacement error (Δδij) of the driven gear along the contact normal (the line of action). The corresponding angular transmission error (Δφij(Tij(k))) under load for the contact position is determined by reversing (5). The column vector [p] 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 a, b, and c. Then, functional relations between loaded transmission errors and some nominal load of Tij(k) may be proposed. The calculation curves are supplied in a meshing cycle. The time-varying mesh stiffness is represented by(12)Kijk=Tijkricos⁡αn·1rbiΔφijTijk.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.Figure 5
The schematic of supporting rigidity.In the process of the spline transmit torque, friction will be produced between internal and external spline.Fm=τFm is friction; here, Fm is the positive pressure between internal and external spline and τ is friction coefficient. The floating quantum can be represented by(13)D-ri=X-in2+Y-in2.Here, X-i(n) and Y-i(n) are the floating quantum along the x-direction and y-direction, respectively; n 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 Figure5 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 S1). 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 (S1 to S2). However, if the slippage is beyond S2, 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. S1-S2 represent the radial clearance between the internal and external spline.Δ
F
x
1
(
n
) and ΔFy1(n) represent the support reaction of floating pinion projected on the x-axis and y-axis, respectively:(14)ΔPx1n=-Kw·D-r1·cosξ1n0≤D-r1≤S1-Fm·cosξ1nS1<D-r1≤S2-Fm+Kw·D-r1-S2·cosξ1nD-r1>S2,ΔPy1n=-Kw·D-r1·sinξ1n0≤D-r1≤S1-Fm·sinξ1nS1<D-r1≤S2-Fm+Kw·D-r1-S2·sinξ1nD-r1>S2.Here, Kw is the flexural rigidity of spline shaft; ζn is a direction angle of vector of (x1n,y1n).Δ
D
x
1
(
n
) and ΔDy1(n) represent the damping force of floating pinion projected on the x-axis and y-axis, respectively:(15)ΔDx1n=-Cw·D-r1·cosξ1n0≤D-r1≤S1-CFm·cosξ1nS1<D-r1≤S2-CFm+Cw·D-r1-S2·cosξ1nD-r1>S2,ΔDy1n=-Cw·D-r1·sinξ1n0≤D-r1≤S1-CFm·sinξ1nS1<D-r1≤S2-CFm+Cw·D-r1-S2·sinξ1nD-r1>S2.Here, Cw is the damping of spline shaft bending and CFm is frictional damping.The support equilibrium conditions of floating pinion can be represented by(16)M1x¨1+ΔPx1n+ΔDx1n+P12+D12cosζ12+P13+D13cosζ13=0,M1y¨1+ΔPy1n+ΔDy1n+P12+D12sinζ12+P13+D13sinζ13=0,meq,1u¨1+P12+D12+P13+D13=0.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 Table1. Input power P=373Kw and input speed n1=8780 r/min.Table 1
Gear parameters.
Gear
Teethz
Pitch diameterd/(mm)
Tooth widthb/(mm)
Pressure angleα/(°)
Helical angleβ/(°)
1
32
51.1
44.5
20
6
2, 3
124
197.9
38.1
20
6
4, 5
27
68.6
66
20
0
6
176
447.0
59.9
20
0The loaded transmission errors (LTE) of each gear pair of the system related to Table1 are shown in Figure 6.Figure 6
Transmission errors. Illustration of (a) I-stage helical gear pairs and (b) II-stage spur gear pairs.
(a)
(b)Figure6 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 (-7.21~0.12)/′′. Figure 6(b) is the LTE of level II gear pair. Here, the LTE fluctuation is (-11.5~0.11)/′′. 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 Figure7; meanwhile, the time-varying mesh stiffness incentive based on the empirical formula (IS06336.1996) [19] is shown.Figure 7
The curves of time-varying mesh stiffness incentive. Illustration of (a) level-I helical gear pairs and (b) level-II spur gear pairs.
(a)
(b)Figure7 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 (7.29×105~1.17×106)/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 (1.09×105~1.27×106)/N·mm−1.Differential equation is expressed in formula of[M]{X′′}+[C]{X′}+([K(t)]-e(t)){X}={P0}. The natural characteristics of problem can be converted into solving the equations of ωi2[M]=[K-], [M] is system mass matrix, [K-] 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.Table 2
The natural frequency analysis results.
Order
Natural  frequency
Order
Natural  frequency
Order
Natural  frequency
Order
Natural  frequency
Order
Natural  frequency
1
0
7
199.877
13
468.116
19
1018.471
25
2084.417
2
142.802
8
273.048
14
523.136
20
1018.471
26
2084.418
3
152.281
9
322.915
15
824.536
21
1748.084
27
3774.596
4
153.447
10
354.337
16
824.536
22
1794.751
28
3774.596
5
191.011
11
386.456
17
824.544
23
2084.417
29
11694.802
6
191.011
12
387.002
18
824.550
24
2084.417
30
11694.802As seen from Table2, 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 Figure8. It can be seen from the vibration mode curve, the node between 23 and 29 of the fundamental frequency is the largest relative amplitude.Figure 8
Each natural principal mode of system.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. Figure9 is time-history dynamic load and FFT spectrum diagram of the gears at all levels of system.Figure 9
Dynamic load time history and FFT spectrum under standard conditions.
(a)
The I-stage dynamic load and the corresponding FFT spectrum
(b)
Level II dynamic load and the corresponding FFT spectrum diagramIt can be seen from Figure9 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 Figure10 with the standard conditions.Figure 10
Vibration acceleration time history and FFT spectrum under standard conditions.
(a)
Level I acceleration and the corresponding FFT spectrum
(b)
Level II acceleration and the corresponding FFT spectrumIt can be seen from Figure10 that the vibration acceleration RMS (m/s2) 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 Figure11.Figure 11
II-stage time-history vibration acceleration and FFT spectrum.It can be seen from Figure11 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 Figure12; here, ξg is meshing damping ratio.Figure 12
The relationship between the dynamic factor and speed.
(a)
Level I dynamic load coefficient
(b)
Level II dynamic load coefficientIt can be seen from Figure12 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ΔAx2=ΔAx4=ΔAx6=0.05 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.Figure 13
Load-sharing characteristics changed with installation error.
(a)
Gear 2 error ofΔAx2
(b)
Gear 4 error ofΔAx4
(c)
Gear 6 error ofΔAx6Figure13 shows that the torque is cyclically fluctuating at each meshing position in different errors of ΔAx2,ΔAx4, and ΔAx6, which reflect the load distribution at different engagement position in the tooth surface. With the influence of ΔAx2,ΔAx4, and ΔAx6, 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 Figure14.Figure 14
Load-sharing coefficient curves with the change of installation errors.Figure14 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 Figure15.Figure 15
Dynamic load under the influence of the spline clearance.
(a)
Spline clearance is 0 mm
(b)
Spline clearance is 0.85 mm
(c)
Spline clearance is 1.231 mmFigure15 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 Figure16.Figure 16
Load-sharing coefficient changed with spline clearance.
(a)
Clearance floating coefficient changes withΔAx2
(b)
Clearance floating coefficient changes withΔAx4It can be seen from Figure16 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).Figure 17
Experimental equipment and the principle diagram.
(a)
Experimental split torque gearbox test facility
(b)
Power closed circuit principle diagramTheoretical 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 × 105 N·mm, 7.86 × 105 N·mm, and 1.02 × 107 N·mm. In order to make a better comparison, the unit of inch-pounds (in.-lb) of [8] was transformed to the N·mm.Figure 18
Contrast between theory analysis and [8] calculation of dynamic load.
(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 analysisIt can be seen from Figure18(a) that dynamic load of level I gear pair is fluctuated 4.05 × 105 N·mm from [8]. Figure 18(b) shows that the dynamic load of level I gear pair is fluctuated 4.11 × 105 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 × 105 N·mm from [8]. Figure 18(d) shows that the dynamic load of level II gear pair is fluctuated 7.92 × 105 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 × 106 N·mm from [8]. Figure 18(f) shows that the dynamic load of level III gear pair is fluctuated 1.14 × 106 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.Figure19(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.Figure 19
FFT spectrum of dynamic load contrast between theory analysis and [8] calculation.
(a)
FFT spectrum diagram from [8] analysis
(b)
FFT spectrum diagram from the theory analysisIt can be seen from Figure19(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).Figure 20
Contrast between theory analysis and [8] calculation of dynamic load with installation errors.
(a)
Dynamic load from [8] analysis
(b)
Dynamic load from the theory analysisThe 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.



---

*Source: 871894-2015-10-19.xml*