Dynamic Load Sharing Behaviours of Planetary Gear Trains and Parameter Study through Perturbation Analysis

In this paper, an analytical solution on the dynamic mesh forces of planetary gear trains (PGTs) is proposed by investigating a lumped-parameter model. By using themethod of multiple-scales (MMS), closed-form expressions of mesh force under the effects of manufacturing and assembly errors are obtained. From these expressions, the effects of several key factors such as the tooth thickness error, pin position error, applied torque, support stiffness of sun gear, and tooth profile modifications (TPM) on dynamic load sharing behaviours are explored. Numerical integration is carried out to verify the validation of the proposed method, and the developed expressions are also validated by comparing the results with previously published predictions. *e results for several examined PGTsystems show that the key factors abovementioned affect the dynamic load sharing behaviours as both static and dynamic factors. An important new conclusion obtained by this work is that proper tooth profile modifications keep the dynamic load sharing factors almost equal to the results obtained under static conditions. *is conclusion provides the possibility to simplify the dynamic analysis to the static analysis on the dynamic load sharing problems.


Introduction
In the applications of high-power mechanical transmission areas, to guarantee the load on each component within the safe ranges, split-torque transmission systems are widely used. Planetary gear train (PGT) is one of the most popular split-torque transmission systems, because of its advantages such as high-power density, high transmission ratio, low bearing load, and compactness.
eoretically, the input torque should be shared evenly by all the planets. In practical applications, however, each path of a PGT will carry uneven load due to the presence of errors, such as tooth thickness error and pin position error. e load sharing behaviours are affected by a variety of factors such as the gravity, support stiffness of central components, bearing clearance, backlashes, applied torque, and component flexibility [1]. Uneven load sharing leads one or several planet gears to carry load which is more than the nominal value, which shortens the lifetime of PGT systems and increases potential damaging risk.
In the past 20 years, nonlinearities in geared systems caused by factors including backlash [2], time-varying mesh stiffness [3], friction [4] between the mating gear teeth, and the error of the gear transmission have been the main concerns of researchers. Neglecting the effects of dynamic vibrations, the quasistatic load sharing behaviours of PGTs have been extensively studied [5]. Previous quasistatic analysis showed that the load sharing behaviours of PGTs are sensitive to the manufacturing errors [6], and this sensitivity increases with the increasing of the number of planet gears [7]. e results obtained from finite element models showed that load sharing is dramatically affected by the applied torque [8]. Floating one or several central components can improve the load sharing and decrease critical tooth stress [9]. ese aforementioned conclusions obtained by quasistatic analysis have been confirmed by experimental results [10][11][12]. Ligata et al. [13] presented planet load sharing formulas under quasistatic conditions, and these formulas revealed the influence of various parameters and errors on the load sharing. Singh [14,15] gave a physical explanation for the basic mechanism which causes the unequal load sharing phenomenon. Quasistatic analysis shows the mechanism of the influence of various parameters and errors on the load sharing as static factors and provides researchers a basic understanding on the loading sharing problems. However, these studies cannot explore the effects of the dynamic factors which are unavoidable in practical PGTs under the actual working conditions.
In order to further predict the load sharing behaviours under practical dynamic working conditions, some researchers investigated the dynamic load sharing by using lumped-parameter models. Kahraman [16] established a relationship between dynamic load and static load by defining the dynamic load factor to indicate the effects of dynamic factors on mesh load. Considering the planet pin position errors [17,18] and eccentricities [19], Gu and Velex investigated the influence of centrifugal forces and pin flexibility on system with rotating carriers. Results indicate the centrifugal forces of rotating carrier may change the dynamic load sharing behaviours. It is believed that the manufacturing errors may significantly increase the dynamic vibrations [20][21][22], and the dynamic behaviours have different sensitivity to different kinds of errors [23]. Mo et al. [24,25] conducted analytical investigation of load sharing characteristics for herringbone planetary gear train with flexible support and multipower face gear split flow system. e impacts of manufacturing errors, assembly errors, manufacturing error phases, assembly error phases, meshing damping, support stiffness, and the input power on the load sharing coefficients were analyzed. e load sharing may be improved by modifying the tooth profile [26] and floating one central component [27]. However, the floating central component will lead large vibrations [28]. e aforementioned dynamic analysis was proposed by numerical methods. ese numerical methods may provide accurate solutions; however, they were not able to reveal the mechanism of the influence of the system parameters and the variety of errors on the dynamic load sharing.
In this effort, an analytical solution on the dynamic load sharing behaviours of PGTs is proposed by investigating a lumped-parameter model. Based on the previous works [29,30] and the authors' efforts [31,32], the method of multiple-scales (MMS) is applied to analyse the dynamic load sharing behaviours under the effects of tooth thickness errors and the pin position errors for the first time. Closedform expressions of mesh force versus mesh frequency are derived by MMS. rough these expressions, the effects of several key parameters, tooth thickness errors, pin position errors, and tooth profile modifications (TPM) on the dynamic load sharing behaviours are illustrated. Numerical integration is employed to verify the proposed method. is study provides guidance for the selections of the key parameters of the PGT systems. e rest of this article is organized as follows. In the second section, the modelling and problem definition are presented. In the third section, closed-form expressions of mesh forces are driven through the method of multiple-scales (MMS). Analytical and numerical results and discussions are presented in the fourth section. Finally, conclusions are collected together and presented.

Modelling with Various Kinds of Errors.
A lumped-parameter model of PGTs with a sun gear, a fixed carrier, a ring gear, and N planet gears is shown in Figure 1. In this model, only the rotational motions of each component and two translational motions of the sun gear are considered. e number of degree of freedom is N + 5. In Figure 1, k bs is the translational stiffness of the sun gear. k sp and k rp are the mesh stiffness of the pth s-p and r-p mesh, respectively. k su , k ru , and k cu are the torsional stiffness of sun, ring, and carrier, respectively. u denotes the rotational motions. Figures 2 and 3 are two diagrams to describe the r-p and s-p mesh, respectively. In these two figures, e and h are equivalent tooth thickness errors and the gap induced by tooth profile modifications, respectively. δ rp and δ sp are respectively the relative mesh deformation on r-p and s-p meshes caused by the motions of the components. Considering the mesh force between the ring gear and planet gears, one can write the equation of ring gear as where δ rp � − u c cos α r + u r − u p , where f erp and f mrp are respectively the additional mesh forces induced by errors and TMP. c denotes the damping term, and α r is the pressure angle of r-p mesh. J is the moment and r is the base radius. e rp and h rp are respectively the tooth thickness errors and TPM functions on r-p mesh. e tooth separation is modelled by Θ(δ), where δ is the deformation of the gear mesh. Figure 4 shows the effects of the backlash on the mesh force of the tooth pair. When δ > 0, where δ is the deformation of the gear mesh, the tooth pair maintains contact; when − b < δ < 0, where b is the width of the backlash, contact loss occurs, and when δ < − b, backside tooth contact occurs. Backside tooth contact is not normally observed in practice because of the gear preload and backlash, so it is neglected. e tooth separation functions Θ(δ) are e sun gear links with N planet gears and the equation of motion could be expressed as 2 Shock and Vibration  where δ sp � − u c cos α s − x s sin ψ sp + y s cos ψ sp + u s + u p , where x s and y s are the translational motions of sun gear in the x and y directions, respectively, α s is the pressure angle of s-p mesh, Ψ p is the angle between x-axis and the line oo p and Ψ sp � Ψ p − α, T s is the transmitted torque applying on sun gear, and e sp and h sp are respectively the tooth thickness errors and TPM functions on s-p mesh. From equations (1) and (4), the equation of motion associated with the pth planet gear is As shown in Figure 5, all the planet gears are fixed on the carrier and the equation of carrier rotational motion is where r c is the length of oo 1 , as shown in Figure 1. As shown in Figure 6(a), the tooth thickness error is the deviation between the actual tooth thickness and the design involute tooth thickness. When the actual tooth is thicker, the additional mesh force is positive and vice versa. e pin position error is shown in Figure 6  is not shown. Converting all the errors to the gear pair meshing lines, the equivalent tooth thickness errors obtained through superposition could be expressed as In equation (8), E pr , E ps , and E pp represent the tooth thickness errors of ring, sun gear, and the pth planet gear, respectively. E cc , E cr , E cs , and E cp represent the pin position errors of the carrier, ring, sun gear, and the pth planet gear, respectively. E is the time-invariant component and E is the amplitude of the time-varying component. E sp and E rp are the amplitudes of the tooth thickness errors of s-p and r-p mesh, respectively. φ is the angle between the directions of the corresponding pin position error and the mesh line, and c is the initial phase of the time-varying errors. Ω is the rotational speed and t denotes time. ω m is the mesh frequency, which is determined by the gear tooth and the power flow path as a function of rotational speed, as given in the following equation [16]: Tooth thickness error A ct ua l to ot h D es ig n to ot h where Z s and Z r are the tooth number of sun gear and ring gear, respectively. en, the additional mesh force introduced by equivalent tooth thickness errors on the pth s-p mesh and r-p mesh is obtained, ere are several methods to modify gear tooth surfaces, including crowning, tip relief, and root relief with linear or parabolic variations with roll angle. e TPM curves, the magnitude of the relief, and the modification length are the three key factors that determine the effects of the TPM on vibration reduction. Without loss of generality, linear relief is applied to double tooth pair contact areas about the tooth tip and root in this study. As shown in Figure 8, for the spur PGTs, TPM is applied on the double teeth contact area about the tooth tip and root. e additional mesh force introduced by TPM on the pth s-p and r-p mesh could be expressed as From equations (1)-(7), the system equation in matrix form is x � u c , u r , x s , y s , u s , u 1 , . . . , u N T .
M is the mass matrix, K b is the support stiffness matrix between the PGT and the fixture. K mv and K m0 are the varying part and mean part of stiffness matrix, respectively.  Figure 7: Contributions of different errors to the equivalent tooth thickness errors. 6 Shock and Vibration where k and k are the average and varying components of the time-varying mesh stiffness, respectively. Considering the tooth separation, one can write the mesh stiffness as In this study, the rectangle waves [30,33] are applied to approximate the time-varying mesh stiffness. As shown in Figure 9, mesh stiffness varies as the number of contact tooth pair changes. K sp is the nondimensional mesh stiffness matrix and can be written as K rp has a similar form as K sp . C is the damping matrix in the form of F t is the external load vector. F d and F m are respectively the inner force vectors introduced by the equivalent tooth thickness errors and TPM. Shock and Vibration 7 where D sp and D rp are the component vectors of f esp (x, t) and f erp (x, t) on each degree of freedom. e time-varying mesh stiffnesses and the additional mesh forces induced by errors and TPM can be expressed in terms of Fourier series as where c.c. stands for the complex conjugate of the preceding terms. e eigenvalue problem associated with the linear free vibration of the PGTs is where v i are the eigenvectors and ω i are the natural frequencies. Applying the modal transformation, x � Vz, and let where where s ⟶ r represent the corresponding terms of r-p meshes. From the results of modal analysis, one can find that the 1 st mode is rigid mode, the 2 nd and (N + 2) th are rotational modes with distinct natural frequency, and the 3 rd to (N + 1) th modes have equal natural frequency. It has been proved that, the vibration of PGTs system is mainly the superposition of these two rotational modes [34]. Figure 10, the time histories of s-p mesh forces of a 5-planet example system with given parameters in Table 1 are presented. Without any error, the mesh forces for all the s-p mesh pairs are equal. As shown in Figure 10(a), the value of mesh force at point A is the maximum. Let f δsp and f δrp denote the peak value of s-p and r-p mesh forces, respectively. With tooth thickness error on the 1 st gear mesh E p1 � 10 μm, the corresponding time histories of s-p mesh forces are shown in Figure 10(b). It is clear that the 1 st s-p mesh force f δs1 is dramatically increased and f δs2 � f δs5 , f δs3 � f δs4 . e ideal load condition for a PGT system is that each path will carry an equal load and the dynamic mesh force values are low. To illustrate the dynamic load sharing, the load sharing coefficients of the s-p mesh L sp and r-p mesh L rp are defined as

Problem Definition. In
e s-p load sharing factor L s and the r-p load sharing factor L r are defined as In order to indicate the maximum mesh force, the dynamic load factor Γ is defined as In the case of the PGTs manufacture perfectly without error, as shown in Figure 10(a), L s � 1 and Γ � 1.055. In the case of E p1 � 10 μm, as shown in Figure 10(b), L s � 1.13 and Γ � 1.19. It is clear that the error causes uneven mesh forces and the increasement of the maximum mesh forces.

Perturbation Analysis on Dynamic Mesh Forces
To investigate the dynamic load sharing and load factor of PGTs, an approximate solution for the mesh force is sought by using the MMS [35]. e tooth separation function can be approximated as e approximated expression of Θ rp is similar to that of Θ sp . e phases ϕ sp and ϕ rp will be chosen subsequently such that the tooth separation is in-phase with the mesh deflection. Similar to the procedure as presented in references [24,27], an expression of the response amplitude a i and mesh frequency ω m can be obtained as Shock and Vibration 9 sp G spii + s ⟶ r, where f tw , f dw , and f mw are the wth element of the vectorsF vt F vd , and F vm , respectively. G spiw is the (i, w) elements of G sp and R spi is the ith elements of R sp . ζ i is the damping ratio. e dynamic mesh forces include the static components and timevarying components. e peak values of the time-varying s-p and r-p mesh forces can be expressed as where v δsp and v δrp are two coefficients associated with the modal vectors V and can be expressed as Substituting equation (29) into the definition of the dynamic load sharing coefficients in equation (24), Let Υ s represent the equal parts among all the s-p meshes: Substituting equation (32) into equation (31), equation (31) could be simplified as From equation (33), one can see that the load sharing coefficients are determined by Υ s and e sp . Υ s includes the effects of the transmit torque, tooth profile modification, and dynamic vibrations, and the equivalent tooth thickness error e sp includes the effects of both tooth thickness errors and the pin position errors.

Results and Discussions
In this investigation, the effects of the errors, the support stiffness of sun gear, the applied torque, and tooth profile modification on the load sharing factors and the dynamic load factors will be discussed following. e main parameters of the example systems with different number of planet gear are listed in Table 1. As shown in Table 1, the only difference among the given several PGT systems is the number of planet gears. e two natural frequencies associated with the first two rotational vibration modes, ω r1 and ω r2 , are listed in Table 2.

Validation of the Proposed Method.
To verify the proposed method, assume the 1 st planet gear has the timeinvariant tooth profile errors, e r1 � e s1 � 50 μm, and the other components are perfectly manufactured and assembled for the examined PGTs. It should be noted that these equivalent tooth thickness errors on the 1 st planet may come from either the tooth thickness error or the pin position errors of the 1 st planet gear. e comparisons of the load coefficients obtained by MMS and NI with varying number of planet gear are shown in Figures 11 and 12. e loads on r-p meshes are mainly concerned withω m ≈ ω r1 , while the loads on s-p meshes are mainly concerned with ω m ≈ ω r2 .
With sun gear fixed, as shown in Figure 11, the 1 st planet gear takes the heaviest load due to the errors. Because of the symmetry positions to the 1 st planet gear, the 2 nd and N th planet gears carry equal loads; the 3 rd and (N − 1) th planet gears carry equal loads. With sun gear floated, as shown in Figure 12, all the planet gears in 3-planet PGT systems carry almost equal load and the opposite planet gears in 4-planet gears carry equal load. With the increase of the number of planet gears, the load sharing coefficients become more sensitive to the errors. ese conclusions agree well with the conclusions obtained by previous static analysis [6,7,10]. In Figures 11 and 12, all the curves obtained by MMS agree well with the NI results.
One can find the dynamic load sharing coefficients vary with the changing of mesh frequency. is is because the amplitudes of mesh forces are the functions of the modal vibration amplitudes a i , which is associated with the mesh frequency, as expressed in equation (28). e interesting phenomenon is that, with ω m ≈ ω r1 , ω r2 , the load sharing conditions are improved.
is is because the vibration amplitudes are increased, and the dynamic mesh forces caused by the vibration are the dominant factors in the primary resonance ranges. Consequently, the effects of errors are relatively decreased, and the load sharing coefficient associated with the 1 st planet gear is decreased. However, with this improvement of the load sharing, the dynamic load conditions get worse, as show in Figure 13. e dynamic load factors dramatically increase near the primary resonance due to the large vibrations. Under dynamic working conditions, both of the load sharing factors and dynamic load factors are considered to evaluate the load conditions among all the planet gears.
It should be noted that, the softening nonlinearity and vibration jump phenomenon appear in the primary resonance ranges. e frequency-load sharing coefficient curves have three branches near the primary resonance, and the middle branch is unstable.

Effects of the Errors in Planet Gears.
It is believed that both of the planet pin position errors and tooth thickness errors are parts of the vibration excitations of PGTs [31,32]. In order to explore the effects of tooth thickness errors and pin position errors of planet gear on the planet load sharing, especially on the dynamic load factor, the effects of the equivalent tooth thickness error on the load sharing factors and dynamic load factors are discussed here. Figure 14 shows the effects of e r1 and e s1 for the 5-planet system versus mesh frequency. With the increasing of e r1 and e s1 , the load sharing factors increase and the frequency ranges with contact loss expanded. It is worth noting that the increase of the load sharing factors is almost proportional to the increase of e r1 and e s1 , as shown in Figure 15. e results in Figure 15(a) have similar trend with the results presented in references [9,15]. According to Bodas and Kahraman [6], the load sharing is proportional to the absolute magnitude of carrier and gear manufacturing errors. From Figure 15(b), one can see that, in dynamic conditions, the load sharing is also proportional to the absolute magnitude of equivalent tooth thickness errors. However, with same amount of equivalent tooth thickness error, the load sharing in dynamic conditions is better.
It is believed that the error induces vibration and then leads the increase of the dynamic load factor, as shown in Figure 16. In addition, by this figure, one can find the dynamic load factor is proportional to the equivalent tooth thickness error. Selecting high-precision gears will certainly reduce the dynamic mesh load and improve load sharing performance.

Effects of Pin Position Error of Central Components.
Pin position error of central components is a critical factor affecting the dynamic load sharing. In Figure 17 load sharing, as shown in Figures 17(c) and 17(d). e load of the 1 st planet is minimum. is is because the 1 st planet is furthest to the ring gear with this pin position error. ese results also agree well with the conclusions obtained from static analysis in reference [19]. e effects of the pin position errors on dynamic load sharing factors versus mesh frequency are shown in Figure 18. e load sharing factors are proportional to the pin position errors of sun gear and ring gear, and the frequency ranges with tooth separations are increased by the increasing of these errors. e curves in Figure 19 indicate that there exist little effects of the pin position errors on the dynamic factors. It suggests that the pin position errors of sun gear and ring gear affect the load sharing behaviours mainly as static factors.

Effects of Applied Torque.
Applied torque is another key factor that affects the dynamic characteristic of PGT system. Figure 20 shows the effects of the applied torque on the dynamic factors versus mesh frequency with e r1 � e s1 � 50 μm. From Figure 20, one can see that heavy applied torque effectively decreases the dynamic load sharing factor. From equations (32) and (33), a reasonable explanation is that large torque increases the value of Υ s , and the effects of the errors on the load sharing factors are suppressed. Heavy applied torque may decrease the load sharing factors but increase the loads on all the meshes. e load sharing coefficients of 5-planet system with varying applied torque in both quasistatic and dynamic conditions are shown in Figure 21. e curves in Figure 21(a) agree well with the results obtained under quasistatic conditions in reference [8]. Once again, the dynamic load sharing is better than the corresponding results in quasistatic, because of the effects of vibrations. Considering the static strength of gear tooth, increasing the applied torque to suppress the effects of equivalent tooth thickness error is not the first option.

Effects of the Support Stiffness of Sun Gear.
In order to explore the mechanism of the float sun gear on the improvement of the dynamic load sharing, the comparisons of the translational displacements of sun gear of a five-planet system with/without float sun gear are shown in Figure 22.
With sun gear floated, the equilibrium position of sun gear is put away from the origin point, and the vibration amplitudes of the translational displacements are increased. From equation (29), one can find that float sun gear with low support stiffness allows the equilibrium position to change and induces the translational motions, which compensates the uneven effects of the equivalent tooth thickness errors.
In practical PGTs, the sun gear is not absolutely float or fixed. To further explore the effects of the sun gear support stiffness, the influences of the support stiffness on the load sharing factors and the dynamic load factor of sun gear for the 5-planet system are shown in Figures 23 and 24, respectively. e load sharing factor decreases with the support stiffness of sun gear decreasing, and the load sharing is improved. e dynamic load factor slightly decreases with the decreasing of the support stiffness of sun gear.
In the case of low rotational speed, the load sharing condition is close to those in quasistatic condition. As shown in Figure 25(a), the load sharing coefficients are almost constant whenk bs ∈ [10 6 , 10 8 ] N/m. While k bs increases up to [10 8 , 10 10 ] N/m, the load sharing coefficient of 1 st planet gear increases. ese curves in Figure 25 agree well with the results presented in reference [9] in the case of quasistatic conditions. In dynamic conditions, as shown in Figure 25(b), the varying of load sharing coefficients with the support stiffness of sun gear has similar trend with the results in quasistatic conditions. While in dynamic load sharing, under the influence of system vibration, the 1 st planet gear load sharing coefficient is decreased and the load sharing is improved. As discussed above, floating sun gear decreases the dynamic mesh load and optimizes the load sharing, which is still one of the effective methods with high priority in dynamic conditions despite of larger sun vibrations.

Effects of Tooth Profile
Modifications. TPM has been proved as an effective method to decrease the vibrations of PGT system [30][31][32]. Without TPM, tooth separations occur when ω m ≈ ω r1 , ω r2 and the "soften" phenomenon appears, as shown in Figures 26(a) and 26(b). Aiming at eliminating the tooth separations of both s-p and r-p meshes with ω m ≈ ω r1 , ω r2 , simultaneously, the proper magnitudes of tooth profile modifications for the s-p and r-p meshes are 28μ m and 15μ m, respectively. ese optimal TMP amounts could be obtained by using the method presented in Ref. [31]. With proper TPM, the vibration amplitudes are dramatically decreased, and tooth separations are eliminated, as shown in Figures 26(c) and 26(d). Proper TPM suppresses the contributions of vibration to the mesh forces.
In Figures 27(a) and 27(b), the load sharing factors L r and L s are shown for the system with/without TPMs. By comparing the results obtained for the system without TPMs with the results obtained for the system with TPM, one can say TPM can help to eliminate the fluctuation of the load sharing factors with the change of mesh frequency, as shown in Figures 27(a) and 27(b). at means, with proper TPM, it is very reasonable to approximate the dynamical load sharing factors with the result obtained in quasistatic conditions. is conclusion provides the possibility to simplify  the dynamic analysis to the static analysis on the dynamic load sharing problems.
Since the TPM is effective to decrease the system vibration amplitudes, TPM must be an effective method to decrease the dynamic load factors. As shown in Figures 27(c) and 27(d), TPM dramatically decreases the dynamic load factors in the primary resonance frequency ranges. e static and quasistatic results of the load sharing factors agree well with dynamic analysis, and the static and quasistatic analysis has advantages in analysis efficiency. However, dynamic analysis can offer a better understanding of the dynamic load factors. So, selection of static and dynamic analysis depends on the main focus on the PGTs.

Conclusion
In this study, a simplified discrete model is presented to investigate the load sharing among the planet meshes of PGTs with several type errors. Both of the cases of fixed and float sun gear are investigated to study the effects of the support stiffness of sun gear on the load sharing. Timevarying mesh stiffness and tooth separations are also considered. e method of multiple-scales (MMS) is used to obtain the response and closed-form expressions of mesh force are derived over the important mesh frequency ranges. From these expressions, the effects of several key factors such as the tooth thickness and pin position errors, applied torque, support stiffness of sun gear, and tooth profile modifications on dynamic load sharing behaviours are explored. e validation of MMS is obtained by the results of numerical integration and previously published predictions. Several conclusions are obtained: (1) e amplitudes of dynamic mesh force are the function of vibration amplitude which is associated with the mesh frequency. at means the load sharing factors and the dynamic load factors are the functions of mesh frequency.
(2) For PGTs with different number of planet gears, the load sharing coefficients have similar trend versus mesh frequency. With equivalent tooth thickness error on the 1 st planet gear, the dynamic load sharing factor is proportional to the absolute magnitude of equivalent tooth thickness errors. While with the increasing of the planet gear number, the load sharing factors and dynamic load factors became more sensitive to the equivalent tooth thickness errors. (3) With equivalent tooth thickness error on the 1 st planet gear, for 3-planet gear PGTs, the 2 nd and 3 rd planet gears carry equal load because of the geometric symmetry of position. For the same reason, the 2 nd and 4 th planet gears for 4-planet PGTs carry equal load. For 5-and 6-planet systems, the 2 nd and N th and 3 rd and (N − 1) th planets carry equal load, respectively. (4) Floating sun gear decreases the dynamic mesh load while optimizing the load sharing although the vibration amplitude of sun gear increases. at means floating sun gear is one of the effective methods to improve the dynamic load sharing with high priority.   Dynamic load factors Θ: Tooth separation function Ψ p : Position angle of planet gear α s , α r : Pressure angles φ: Angle between pin position errors and the mesh line Υ: Initial phase of errors ω m : Mesh frequency ω i : Natural frequency δ: Mesh deflection ζ: Damping ratio.

Data Availability
e MATLAB data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that they have no conflicts of interest.