Research Article  Open Access
Shuai Mo, Ting Zhang, Guoguang Jin, Zhanyong Feng, Jiabei Gong, Shengping Zhu, "Dynamic Characteristics and Load Sharing of Herringbone Wind Power Gearbox", Mathematical Problems in Engineering, vol. 2018, Article ID 7251645, 24 pages, 2018. https://doi.org/10.1155/2018/7251645
Dynamic Characteristics and Load Sharing of Herringbone Wind Power Gearbox
Abstract
In this study, the dynamic model for the herringbone planetary gear transmission system is established by the lumped parameter method based on the system dynamics and the Lagrange equation, and the impact of the support stiffness and the torsional stiffness on dynamic characteristics is studied. The research results have a guiding significance for the design of the herringbone gear transmission system. In this model, the herringbone gear is treated as a special gear coupled by 2 opposite helical gears, where the stagger angle, comprehensive meshing error, support stiffness, support damping, and load inertia are considered in the analysis of dynamics. Moreover, the dynamic characteristic of the carrier is considered as well. By calculating the meshing force curve of the transmission system, the impact of the stagger angle, supporting stiffness, and the torsional stiffness on meshing force and load sharing coefficient is analyzed. The results show that the stagger angle has an obvious impact on load sharing coefficient while it has little impact on maximum meshing force. And the support stiffness has a more obvious impact on the dynamic characteristics of the system. The recommendary support stiffness of the system is that all of the support stiffness of the sun gear, planetary gear, ring gear, and carrier is 10^{7} N/m. The torsional stiffness has little impact on the dynamic characteristics of transmission system, except the torsional stiffness of planetary gear, and carrier has an obvious impact on load sharing coefficient. The commercial software ADAMS carried out dynamics analysis of the transmission system to verify the necessity validity of the theoretical analysis.
1. Introduction
Herringbone gear planetary transmission system, which is one of the most important parts of wind power generator, is wildly used in ship, aviation, and other transmission fields because of its distinctive advantages, such as small volume, light weight, large transmission ratio, strong bearing capacity, high transmission efficiency, and so on. The load distribution of planetary gears in the transmission directly affects the working life, stability, and reliability of gears. It is of great practical and engineering significance to carry out the dynamic load characteristics of herringbone gear transmission system.
Many scholars around the world have done many theoretical and experimental studies on the load sharing characteristics of planetary transmission. SHI Wankai et al. [1] studied the transmission dynamic characteristic of wind speed increase gearbox to study the effect on wind turbine gearbox dynamic load characteristics caused by stagger angle and leftright coupling torsional stiffness. ZHANG Linlin et al. [2] researched load sharing characteristic of herringbone gear transmission system influenced by meshing phase and conducted an experiment on different meshing phase planetary transmission system to study load sharing on singlestage herringbone planetary gear test stand. MO Shuai et al. [3–5] proposed a new dynamic model of herringbone star gear transmission system to study meshing error, tooth thickness error, base circle error, and assembly error influencing the load sharing characteristic of GTF aeroengine. Ming LI [6] established a model to predict the reliability of helicopter planetary gear train under the condition with partial load and analyzed structure and part load state of the transmission system in detail based on the study of load sharing characteristics of a helicopter gearbox. The impact of assembly error on load sharing characteristic of closed differential herringbone gear transmission system was researched by ZHU Zengbao [7]. LU Junhua [8] established a calculation model of 2KH transmission system to analyze load sharing characteristic of planetary transmission system from the point of view of dynamics. XU Xiangyang [9] studied impact of stiffness of flexible pin shaft on load sharing of star transmission system and verified the effectiveness of theoretical result. Kahraman A et al. [10–14] researched the impact of gear indexing error on dynamic response of a typical multigrid gear system and compared subharmonic resonance of response with numerical integration result directly to prove the accuracy of prediction and predicted other nonlinear phenomena by nonlinear model and studied the variation of the position errors of the planetary gear and the accompanying star wheel under the condition of light load in the small module planetary gear system as well. Howard I [15] proposed a new dynamic model of flexible supporting gear train to calculate the meshing stiffness with varied gear center distance and coupled it to dynamic model by iteration. Wang J et al. [16] developed a three free degrees torsionalvibration model of locomotive spur gear transmission system, where nonlinear gap, static transfer error, and time varying meshing stiffness and other factors were considered, and the influence of speed and stiffness of pinion on the dynamic characteristics of the transmission system is studied by numerical method. FernandezDelRincon A et al. [17] proposed a model which can simultaneously consider the internal excitation caused by variable meshing stiffness and the variable compliance results of bearing. The model combines with hybrid formula for calculating meshing force and nonlinear variable flexibility of bearing to prove a detail described for meshing force. Liang X et al. [18] focus on modeling, detection, and diagnosis of gearbox failure based on dynamics. The research result was limited on some key and basic aspects by literature review: gear meshing stiffness evaluation, gearbox damage modeling and fault diagnosis technology, gearbox transmission path modeling, and method verification. TAN Yuanqiang [19] researched dynamic and load sharing characteristics of closed differential planetary gearbox by comparative trial, which proves that the average load performance of the closed stage is better than that of the differential stage, and the greater the load, the better the load performance of the two stages. The transmission performance of the planetary system with the three kinds of load sharing methods is better than the single load sharing system. Singh [20] established an analysis model of planetary transmission load sharing characteristic to physically explain basic mechanism of the phenomenon which causes inhomogeneous load. Wink et al. [21] researched the calculation method of load distribution of tooth surface and proved three different methods to solve load distribution of gear teeth in calculation of static transmission error of gear pair.
The herringbone gear planetary transmission system of wind power train was taken as the research object in this paper. The support stiffness, support damping, meshing stiffness, meshing damping, and other factors of herringbone gear in the condition of stagger angle which is 0 are overall considered to study their impact on dynamic characteristic. Compared with other papers, this paper especially researched the impact of carrier on dynamic characteristic and deduced the formula of carrier as one part of system dynamics formula. The dynamic model of herringbone gear planetary transmission system is established to analyze the dynamic characteristics of the system. The theoretical dynamic characteristic of system is verified through simulation.
2. Dynamics Theories and Equations
2.1. Dynamic Model of Planetary Transmission System
The threedimensional model of herringbone planetary transmission system is shown in Figure 1.
(a) 3D model with carrier of system
(b) 3D model without carrier of system
The impact of the input and output shafts is not considered in this paper. Thus, the input torque is directly acting on carrier. Because each singlehelical gear, respectively, has the translational motion along the direction of x, y, and z axes and turning motion in the direction of the z axis, each herringbone gear, which is coupled by 2 helical gear with opposite direction, has 8 free degrees. The carrier has the translational motion along the direction of x, y, and z axes and turning motion in the direction of the z axis. Thus, the carrier has 4 free degrees. Because the herringbone gear is divided into left and right parts, subscript L means the left side gear while subscript R means the rightside gear. The transmission system has free degrees, and the generalized coordinate is
where x, y, z, respectively, are lateral, longitudinal, and axial displacement of a component and θ is the angle displacement of the component. is the ith planetary gear.
According to the 3D model of herringbone planetary transmission system, the dynamic model of herringbone planetary transmission system is established, as shown in Figure 2.
OXY is the fixed coordinate system, and oxy is the moving reference system.
To simplify the angle displacement as the equivalent line displacement on the corresponding base circle, where means the left side and right side of the component. means the jth planetary gear. r_{bs}, r_{bpj}, r_{br}, respectively, are the radius of base circle of sun gear, the jth planetary gear, and ring gear. r_{c} is the distance from the center of planetary gear to gyration center of carrier,where r_{s}, r_{p}, respectively, are the radius of pitch circle of sun gear and planetary gear.
2.2. Equivalent Displacement
Figure 3 shows the projection relation of each component of the planetary gear system on the gear face.
According to Figure 3, projection of relative displacement of each component in the meshing line or coordinate system can be calculated. The projection of relative displacement of sun gear and planetary gear along the direction of meshing line isThe projection of relative displacement of planetary gear and ring gear along the direction of meshing line isAnd the projection of relative displacement of planetary gear and carrier along the directions of x_{c}, y_{c}, z_{c}, and u_{c} iswhere Φ_{pi} is the position angle formed by the ith planetary gear relative to the first planetary gear, , α_{sp} is the end face meshing angle of sunplanetary gear meshing, and α_{rp} is the end face meshing angle of planetaryring gear meshing. e_{sp} is the comprehensive meshing error in the sunplanetary gear meshing line, e_{rp} is the comprehensive meshing error in the planetaryring gear meshing line, L is the left side, and R is the right side.
2.3. Equilibrium Differential Equation of System Dynamics
According to the Newton second law and system dynamics, the differential equation of the dynamics of sun gear can be written as (6) and (7). Equation (6) is the dynamics equation of the left side of sun gear and (7) is the dynamics equation of the right side of sun gear.where m_{s} is the equivalent mass of sun gear, K_{s} is the radial support stiffness of sun gear, C_{s} is the radial support damping of sun gear, K_{sa} is the axial support stiffness, C_{sa} is the axial support damping, K_{s12} is the axial stiffness between two helical gears, C_{s12} is the axial damping between two helical gears, K_{1t} is the leftright coupling torsional stiffness, C_{1t} is the leftright coupling torsional damping, and T_{in} is the input torque.
The axial support stiffness and support damping of sun gear are as follows:
The differential equation of the dynamics of planetary gear can be written as (10) and (11). Equation (10) is the dynamics equation of the left side of planetary gear and (11) is the dynamics equation of the right side of planetary gear.where m_{p} is the equivalent mass of planetary gear, K_{p} is the radial support stiffness of planetary gear, C_{p} is the radial support damping of planetary gear, K_{pa} is the axial support stiffness, C_{pa} is the axial support damping, K_{p12} is the axial stiffness between two helical gears, C_{p12} is the axial damping between two helical gears, K_{2t} is the leftright coupling torsional stiffness, and C_{2t} is the leftright coupling torsional damping.
The axial support stiffness and support damping of planetary gear are as follows:
The differential equation of the dynamics of ring gear can be written as (13) and (14). Equation (13) is the dynamics equation of the left side of ring gear and (14) is the dynamics equation of the right side of ring gear.where m_{r} is the equivalent mass of planetary gear, K_{r} is the radial support stiffness of planetary gear, C_{r} is the radial support damping of planetary gear, K_{ra} is the axial support stiffness, C_{ra} is the axial support damping, K_{r12} is the axial stiffness between two helical gears, C_{r12} is the axial damping between two helical gears, K_{3t} is the leftright coupling torsional stiffness, and C_{3t} is the leftright coupling torsional damping.
The axial support stiffness and support damping of ring gear are as follows:
Figure 4 is the dynamics model of carrier, which combined with the Newton second law can deduct the differential equation of the carrier dynamics, as shown in where m_{c} is the equivalent mass of carrier, K_{c} is the radial support stiffness of carrier, C_{c} is the radial support damping of carrier, K_{ca} is the axial support stiffness, C_{ca} is the axial support damping, K_{ct} is the torsional stiffness, and C_{ct} is the torsional damping.
The axial support stiffness and support damping of carrier is as follows:
All above equations can be arranged as matrix form as (18):where M is the generalized mass matrix, X is the generalized coordinate displacement matrix, C is the damping matrix, K is the stiffness matrix, and T is the external excitation column matrix.
2.4. Meshing Error Excitation
Shortperiod error is mainly composed by base circle error and involute tooth profile error. In the filed of gear transmission, the product of the rotation frequency of the shaft and the number of teeth, which is named meshing frequency, is noted as .where n_{s}, n_{p}, n_{r}, respectively, are the revolution speed of sun gear, planetary gear, and ring gear and z_{s}, z_{p}, z_{r}, respectively, are the teeth number of sun gear, planetary gear, and ring gear.
The circular frequency corresponding to the meshing frequency is .
There will be errors in the manufacture and installation of gears due to objective factors, which is the main factor to make vibration of gear system. The variation rule of short period error can be expressed by Fourier series, which is based on the fundamental frequency of the wave. There is a total of 4N periodic excitation.
Error excitation of inner meshing isAnd error excitation of outer meshing iswhere is the 1st order amplitude of the transfer error between the sun gear and the planetary gear, is the initial phase of the transfer error between the sun gear and the planetary gear, is the 1st order amplitude of the transfer error between the sun gear and the planetary gear, is the initial phase of the transfer error between the sun gear and the planetary gear, is the phase difference between the meshing of ring gear and planetary gear and the meshing of sun gear and planetary gear, and is the stagger angle of herringbone gear.
3. Impact of Stagger Angle on Dynamic Characteristic of System
There is a wind power gearbox whose main parameters are shown in Table 1. And the initial amplitude and phase of gear transmission error can be calculated in Kisssoft. The numerical value of first 3order harmonics amplitude, shown in Table 2, is needed because the numerical value of higher order harmonics amplitude is very small.
 
Note. ξ is an arbitrary value in the range of 0.070.13. ξ is 0.1 in this paper. 

Because of the special tooth shape of herringbone gear, there will be formed unavoidable stagger angle in the processing of gear manufacturing. As Figure 5 shown, left side is the plan sketch of stagger angle while the right side is the 3D model of it.
(a) Sketch of stagger angle is 0
(b) Sketch of stagger angle is
While stagger angle, respectively, is 0, π/4, π/2, 3π/4, and π, the corresponding maximum meshing force can be calculated. Figure 6 shows the maximum meshing force curve.
(a) Inner meshing
(b) Outer meshing
According to the equation of load sharing coefficient, the load sharing coefficient of a single gear pair can be defined aswhere n is the period number of meshing frequency, N is the number of planetary gear, and F is the corresponding meshing force of the meshing cycle.
Likewise, the load sharing coefficient of herringbone gear inner and outer meshing is
The load sharing coefficient of inner and outer meshing can be calculated by (23), (24), and Figure 6, as shown in Figure 7.
(a) Inner meshing
(b) Outer meshing
In summary, the change of stagger angle has obvious influence on the dynamic characteristics of the left side of the herringbone gear but has little effect on the right side. In the processing of stagger angle increase from 0 to π, both the left and right side maximum meshing force have an unobvious change. But the load sharing coefficient of left side decreases at first and increases at last, while the right side decreases all the time. The load sharing coefficient of inner meshing meets the minimum when stagger angle is π/2, and the outer meshing meets the minimum when stagger angle is 3π/4. Overall factors which are comprehensively considered, the best stagger angle of this system is 3π/4.
4. Impact of Support Stiffness on Dynamic Characteristic of System
4.1. Theoretical Dynamics Analysis on Initial Parameter
To analyze the impact of support stiffness on dynamic characteristics of the herringbone gear planetary transmission system from the view of practical application, the parameters researched in this paper were gotten from the wind power gearbox described. In the transmission system, all types of the initial support stiffness of sun gear, planetary gear, and ring gear are 1×10^{7}N/m, and the stagger angle is 0. All of these parameters are substituted in the dynamics program; the meshing force is calculated and shown in Figure 8.
(a) Leftside meshing force of inner mesh
(b) Rightside meshing force of inner mesh
(c) Leftside meshing force of outer mesh
(d) Rightside meshing force of outer mesh
4.2. Theoretical Dynamics Analysis on Changing Support Stiffness
Keep the stiffness of planetary gear and ring gear be a constant value which is 10^{7}N/m, and change the stiffness of sun gear from 10^{6}N/m to 10^{8}N/m. The stiffness, respectively, is 10^{6}, 5×10^{6}, 10^{7}, 5×10^{7}, and 10^{8}N/m, and the corresponding meshing force can be gotten and be coupled as 3D sketch as shown in Figures 9–14.
(a) 3D curve of leftside meshing force
(b) 3D curve of rightside meshing force
(a) Leftside meshing force 3D curve
(b) Rightside meshing force 3D curve
(a) 3D curve of leftside meshing force
(b) 3D curve of rightside meshing force
(a) 3D curve of leftside meshing force
(b) 3D curve of rightside meshing force
(a) 3D curve of leftside meshing force
(b) 3D curve of rightside meshing force
(a) 3D curve of leftside meshing force
(b) 3D curve of rightside meshing force
The maximum meshing force of corresponding support stiffness of sun gear can be coupled as shown in Figure 15. Therefore, the load sharing coefficient curve is shown in Figure 16.
(a) Inner meshing maximum meshing force curve
(b) Outer meshing maximum meshing force curve
(a) Inner meshing load sharing coefficient curve
(b) Outer meshing load sharing coefficient curve
Likewise, keep the stiffness of sun gear and ring gear be a constant value which are 10^{7}N/m, and change the stiffness of planetary gear from 10^{6}N/m to 10^{8}N/m. The stiffness is 10^{6}, 5×10^{6}, 10^{7}, 5×10^{7}, and 10^{8}N/m, respectively. And the corresponding maximum meshing force and load sharing coefficient can be gotten as shown in Figures 17 and 18. Keep the stiffness of sun gear and planetary gear be a constant value which are 10^{7}N/m, and change the stiffness of ring gear from 10^{6}N/m to 10^{8}N/m. The corresponding maximum meshing force and load sharing coefficient can be shown in Figures 19 and 20.
(a) Inner meshing maximum meshing force curve
(b) Outer meshing maximum meshing force curve
(a) Inner meshing load sharing coefficient curve
(b) Outer meshing load sharing coefficient curve
(a) Inner meshing maximum meshing force curve
(b) Outer meshing maximum meshing force curve
(a) Inner meshing load sharing coefficient curve
(b) Outer meshing load sharing coefficient curve
And change the support stiffness of carrier while the support stiffness of gears is a constant value, as shown in Figures 21 and 22. According to the research results, commendatory support stiffness of the studied system is as follows: support stiffness of sun gear 10^{7}N/m, support stiffness of planetary gear 10^{7}N/m, support stiffness of ring gear 10^{7}N/m, and support stiffness of carrier 10^{7}N/m. Curves of meshing fore and load sharing coefficient with commendatory support stiffness are shown in Figures 23 and 24, where load sharing coefficient is written as LSC.
(a) Inner meshing maximum meshing force curve
(b) Outer meshing maximum meshing force curve
(a) Inner meshing load sharing coefficient curve
(b) Outer meshing load sharing coefficient curve
(a) The leftside meshing force of inner mesh
(b) The rightside meshing force of inner mesh
(c) The leftside meshing force of outer mesh
(d) The rightside meshing force of outer mesh
(a) The leftside LSC of inner mesh
(b) The rightside LSC of inner mesh
(c) The leftside LSC of outer mesh
(d) The rightside LSC of outer mesh
In summary, firstly, with the increase of the support stiffness of sun gear, both of the maximum meshing force and load sharing coefficient of inner meshing increase while these of outer meshing decrease. Secondly, with the increase of support stiffness of planetary gear, both of the maximum meshing force and load sharing coefficient of inner meshing increase; the load sharing coefficient of outer meshing increases while the maximum meshing force of outer meshing fluctuates near the support stiffness of 10^{7}N/m. Thirdly, with the increase of support stiffness of ring gear increase, both of the maximum meshing force and load sharing coefficient of inner meshing decrease while these of outer meshing increase. Finally, there is mere change of the maximum meshing force and load sharing coefficient of inner meshing and outer meshing, no matter how much the support stiffness of carrier changes.
5. Impact of Torsional Stiffness on Dynamic Characteristic of System
5.1. Impact of Coupling Torsional Stiffness of Gears on the Dynamic Characteristics
Suppose the support stiffness of system is the commendatory value, the torsional stiffness of leftright helical gear is changed to study the impact of torsional stiffness on dynamic characteristic of system. The torsional stiffness of planetary gear and ring gear is a constant value that is 10^{8} Nm/rad, while changing it of sun gear which is 10^{7}, 5×10^{7}, 10^{8}, 5×10^{8}, and 10^{9} Nm/rad, respectively. The inner and outer meshing force and load sharing coefficient of system are shown in Figures 25 and 26.
(a) The maximum meshing force of inner meshing
(b) The maximum meshing force of outer meshing
(a) The load sharing coefficient of inner meshing
(b) The load sharing coefficient of outer meshing
Likewise, when the torsional stiffness of sun gear and ring gear is a constant value, the curves of maximum meshing force and load sharing coefficient with the changing of torsional stiffness of planetary gear from 10^{7} Nm/rad to 10^{9} Nm/rad are shown in Figures 27 and 28. Change the torsional stiffness of ring gear; the curves of maximum meshing force and load sharing coefficient with the changing of torsional stiffness of ring gear from 10^{7} Nm/rad to 10^{9} Nm/rad are shown in Figures 29 and 30.
(a) The maximum meshing force of inner meshing
(b) The maximum meshing force of outer meshing
(a) The load sharing coefficient of inner meshing
(b) The load sharing coefficient of outer meshing
(a) The maximum meshing force of inner meshing
(b) The maximum meshing force of outer meshing
(a) The load sharing coefficient of inner meshing
(b) The load sharing coefficient of outer meshing
5.2. Impact of Torsional Stiffness of Carrier on the Dynamic Characteristics
For researching impact of torsional stiffness of carrier on dynamic characteristic of system, the torsional stiffness of carrier is changed from 10^{7} Nm/rad to 10^{9}Nm/rad, while the leftright coupling torsional stiffness of sun gear, planetary gear, and ring gear is a constant value. The curves of maximum meshing force and load sharing coefficient are shown in Figures 31 and 32.
(a) The maximum meshing force of inner meshing
(b) The maximum meshing force of outer meshing
(a) The load sharing coefficient of inner meshing
(b) The load sharing coefficient of outer meshing
In summary, the leftright coupling torsional stiffness of sun gear and ring gear has seldom impact on dynamic characteristic when it changes in a reasonable range, and the leftright coupling torsional stiffness of planetary gear has obvious impact on load sharing coefficient but has little impact on maximum meshing force. Likewise, the torsional stiffness of carrier has obvious impact on load sharing coefficient but has little impact on maximum meshing force.
6. Dynamic Simulation Verifications
To verify the validity of theoretical calculation, 3D model of the researched system is built to be simulation analyzed in commercial software ADAMS. The similar simulation curves are fitted together, as shown in Figure 33, where simulation curve of meshing force is written as SCMF.
(a) SCMF of inner meshing of leftside
(b) SCMF of inner meshing of rightside
(c) SCMF of outer meshing of leftside
(d) SCMF of outer meshing of rightside
Meanwhile, according to the simulation meshing force, (23) and (24), load sharing coefficient can be calculated. The curve of load sharing coefficient is shown in Figure 34, where LSC means load sharing coefficient and ROM, RIM, LIM, and LOM mean rightside outer meshing, rightside inner meshing, leftside inner meshing, and lestside outer meshing, respectively.
(a) The Maximum LSC of LIM
(b) The Maximum LSC of RIM
(c) The Maximum LSC of LOM
(d) The Maximum LSC of ROM
In summary, the tendency and phase distribution of meshing force of simulation curves are similar to these curves of theoretical calculation. And the curves of load sharing coefficient of simulation are relatively consistent with these curves of theoretical calculation. Both of them verify the correctness of theoretical dynamic analysis method.
7. Conclusion
(1) Considering the impact of support stiffness and leftright coupling torsional stiffness on dynamics characteristic of herringbone gear transmission and the impact of carrier as one part of system, whose dynamic characteristic will influence the system dynamics, the dynamic equations of herringbone gear transmission system were established by lumped parameter method based on system dynamics and Lagrange equation.
(2) Study the impact of stagger angle of gear on dynamic characteristic of system. Stagger angle has obvious impact on load sharing coefficient but little on maximum meshing force. In the processing of stagger angle increase from 0 to π, both the left and rightside maximum meshing force have an unobvious change. But the load sharing coefficient of left side decreases at first and increases at last, while the right side decreases all the time. The load sharing coefficient of inner meshing meets the minimum when stagger angle is π/2, and the outer meshing meets the minimum when stagger angle is 3π/4, and the best stagger angle of this system is 3π/4.
(3) In the herringbone gear transmission system, the changing of support stiffness will have impact on dynamic characteristics. With the increase of the support stiffness of sun gear, both of the maximum meshing force and load sharing coefficient of inner meshing increase while these of outer meshing decrease. With the increase of support stiffness of planetary gear, both of the maximum meshing force and load sharing coefficient of inner meshing increase; the load sharing coefficient of outer meshing increases while the maximum meshing force of outer meshing fluctuates near the support stiffness of 10^{7}N/m. With the increase of support stiffness of ring gear increase, both of the maximum meshing force and load sharing coefficient of inner meshing decrease while these of outer meshing increase. And, there is mere change of the maximum meshing force and load sharing coefficient of inner meshing and outer meshing, no matter how much the support stiffness of carrier changes. The commendatory support stiffness of this paper is as follows: all of support stiffness of sun gear, planetary gear, ring gear, and carrier is 10^{7}N/m.
(4) The impact of torsional stiffness on system characteristic was considered. The meshing force and load sharing coefficient seldomly changed while the changing of leftright coupling torsional stiffness changed from 10^{7}N m/rad to 10^{9}Nm/rad except the change of coupling torsional stiffness of planetary gear. When the coupling torsional stiffness changed from 10^{7}N m/rad to 10^{9}Nm/rad, the maximum meshing force is nearly a constant value while the load sharing coefficient increases slowly. Moreover, with the change of torsional stiffness of carrier, the meshing force is merely changed and the load sharing coefficient changes a lot.
(5) The dynamic characteristic was simulated by commercial software ADAMS. And curves of meshing force and load sharing coefficient gotten from ADAMS simulation were compared with these of theoretical calculation. The result is shown that both the meshing force curves and the load sharing curves are in good agreement with the theoretical curves, which validates the necessity and correctness of theoretical analysis method.
Data Availability
All data generated or analyzed during this study are included in this published article, and other pieces of information are available from the corresponding author on reasonable request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
This research is supported by National Natural Science Foundation of China (no. 51805368), Natural Science Foundation of Tianjin (no. 17JCQNJC04300), Open Funding of the State Key Laboratory of Materials Processing and Die & Mould TechnologyHuazhong University of Science and Technology (no. P2019022), Fundamental Research Funds for the Tianjin Universities (no. 2017KJ083), Applied Basic Research Project of China Textile Industry Association (no. J201806), and the Program for Innovative Research Team in University of Tianjin (no. TD135037). Thanks also go to Professor Haruo Houjoh, Shigeki Matsumura of Tokyo Institute of Technology, Professor Izhak Bucher of Israel Institute of Technology, and Professor Yidu Zhang of Beihang University for their constant assistances.
References
 W. Shi, H. Qiu, Z. Han, and P. Qin, “Dynamics research of wind turbine gearbox with herringbone planetary gear,” Chongqing Daxue Xuebao/Journal of Chongqing University, vol. 38, no. 1, pp. 95–132, 2015. View at: Google Scholar
 L. Zhang, “Impact of Meshing Phase on Load Sharing for Herringbone Planetary Train,” Journal of Mechanical Engineering, vol. 54, no. 11, pp. 129–140, 2018. View at: Publisher Site  Google Scholar
 S. Mo, Y. Zhang, Q. Wu, F. Wang, S. Matsumura, and H. Houjoh, “Load Sharing Behavior of Star Gearing Reducer for Geared Turbofan Engine,” Chinese Journal of Mechanical Engineering, vol. 30, no. 4, pp. 796–803, 2017. View at: Publisher Site  Google Scholar
 S. Mo, Y. Zhang, Q. Wu, H. Houjoh, and S. Matsumura, “Research on natural characteristics of doublehelical star gearing system for GTF aeroengine,” Mechanism and Machine Theory, vol. 106, pp. 166–189, 2016. View at: Publisher Site  Google Scholar
 S. Mo, Y. Zhang, and Q. Wu, “Research on multiplesplit load sharing of twostage star gearing system in consideration of displacement compatibility,” Mechanism and Machine Theory, vol. 88, pp. 1–15, 2015. View at: Publisher Site  Google Scholar
 M. Li, L. Xie, and L. Ding, “Load sharing analysis and reliability prediction for planetary gear train of helicopter,” Mechanism and Machine Theory, vol. 115, pp. 97–113, 2017. View at: Publisher Site  Google Scholar
 Z. Zhu, R. Zhu, Y. Li, Y. Chen, and Z. Zhu, “Impact of installation error on dynamics load sharing characteristic for encased differential herringbone train,” Jixie Gongcheng Xuebao/Journal of Mechanical Engineering, vol. 48, no. 3, pp. 16–24, 2012. View at: Google Scholar
 J.H. Lu, R.P. Zhu, and G.H. Jin, “Basiccomponents' offset analysis of planetary gearing in need of loading sharing,” Zhongnan Daxue Xuebao (Ziran Kexue Ban)/Journal of Central South University (Science and Technology), vol. 39, no. 1, pp. 143–148, 2008. View at: Google Scholar
 C. Zhu, X. Xu, T. C. Lim, X. Du, and M. Liu, “Effect of flexible pin on the dynamic behaviors of wind turbine planetary gear drives,” Proceedings of the Institution of Mechanical Engineers, Part C: Journal of Mechanical Engineering Science, vol. 227, no. 1, pp. 74–86, 2013. View at: Publisher Site  Google Scholar
 A. AlShyyab and A. Kahraman, “Nonlinear dynamic analysis of a multimesh gear train using multiterm harmonic balance method: periodone motions,” Journal of Sound and Vibration, vol. 284, no. 12, pp. 151–172, 2005. View at: Publisher Site  Google Scholar
 M. Inalpolat, M. Handschuh, and A. Kahraman, “Influence of indexing errors on dynamic response of spur gear pairs,” Mechanical Systems and Signal Processing, vol. 60, pp. 391–405, 2015. View at: Publisher Site  Google Scholar
 A. Bodas and A. Kahraman, “Influence of carrier and gear manufacturing errors on the static load sharing behavior of planetary gear sets,” JSME International Journal Series C Mechanical Systems, Machine Elements and Manufacturing, vol. 47, no. 3, pp. 908–915, 2004. View at: Publisher Site  Google Scholar
 H. Ligata, A. Kahraman, and A. Singh, “An experimental study of the influence of manufacturing errors on the planetary gear stresses and planet load sharing,” Journal of Mechanical Design, vol. 130, no. 4, Article ID 041701, 2008. View at: Publisher Site  Google Scholar
 A. Kahraman, “Load sharing characteristics of planetary transmissions,” Mechanism and Machine Theory, vol. 29, no. 8, pp. 1151–1165, 1994. View at: Publisher Site  Google Scholar
 S. Xue and I. Howard, “Dynamic modelling of flexibly supported gears using iterative convergence of tooth mesh stiffness,” Mechanical Systems and Signal Processing, vol. 80, pp. 460–481, 2016. View at: Publisher Site  Google Scholar
 J. Wang, G. He, J. Zhang, Y. Zhao, and Y. Yao, “Nonlinear dynamics analysis of the spur gear system for railway locomotive,” Mechanical Systems and Signal Processing, vol. 85, pp. 41–55, 2017. View at: Publisher Site  Google Scholar
 A. FernandezdelRincon, P. Garcia, A. DiezIbarbia, A. deJuan, M. Iglesias, and F. Viadero, “Enhanced model of gear transmission dynamics for condition monitoring applications: Effects of torque, friction and bearing clearance,” Mechanical Systems and Signal Processing, vol. 85, pp. 445–467, 2017. View at: Publisher Site  Google Scholar
 X. Liang, M. J. Zuo, and Z. Feng, “Dynamic modeling of gearbox faults: A review,” Mechanical Systems and Signal Processing, vol. 98, pp. 852–876, 2018. View at: Publisher Site  Google Scholar
 Y. Tan, C. Hu, Y. Zhang, W. Wang, S. Xiong, and H. Zuo, “Testing research on dynamic load sharing performance of encased differential planetary gearbox,” Jixie Gongcheng Xuebao/Journal of Mechanical Engineering, vol. 52, no. 9, pp. 28–35, 2016. View at: Google Scholar
 A. Singh, “Load sharing behavior in epicyclic gears: physical explanation and generalized formulation,” Mechanism and Machine Theory, vol. 45, no. 3, pp. 511–530, 2010. View at: Publisher Site  Google Scholar
 C. H. Wink and A. L. Serpa, “Performance assessment of solution methods for load distribution problem of gear teeth,” Mechanism and Machine Theory, vol. 43, no. 1, pp. 80–94, 2008. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Shuai Mo et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.