Abstract
The failure of the key parts, such as bearings and gears, in cutting drum long chain gear driving system of longwall shearer is a crucial factor for service life of gear systems. A short driving system is driven by multiple electric motors for cutting drum has been proposed as an alternative to long chain gear driving system. An electromechanical dynamic model based on lumpedparameter method of short driving system is established considering the effect of intrinsic mechanical characteristic of electric motor. By taking into account timevarying stiffness, damping, electric motor synchronization error, manufacturing eccentric error and transmission error of gears, the loadsharing characteristic and dynamic characteristic of the multiple motor driving unit (MMDU) of short driving system are analyzed under shock load. The results of simulation show that the electric motor synchronization error has significant effect on the loadsharing characteristic of MMDU. The loadsharing coefficient evidently increases as the electric motor synchronization error go up, decreasing the electric motor synchronization error can improve the loadsharing characteristic of MMDU; and the synchronization error must be controlled fewer than 5 % to ensure the loadsharing coefficient of the MMDU under 1.1. The imbalance of dynamic meshing force distribution of pinions of the MMDU can be obviously decreased and the loadsharing characteristic of the MMDU can be greatly improved by controlling the synchronization error of electric motor, but the effect is unobvious by controlling the manufacturing eccentric error or transmission error. It provides an effective method to choose a suitable control ranges under different synchronization errors of electric motor according to different loadsharing requirements for the MMDU.
1. Introduction
Coal mining is a very complex system and process. The harsh working conditions, hazardous environment and the high international demand for this critical and valuable mining product are the most important factors that affect the coal mining process. Longwall mining, as the most important and applied coal mining methods has many differences from other mining methods because of its high mechanization and productivity [1] and commonly used longwall mining equipment is made up of midbody, haulage unit, hydraulic system and cutting unit as shown in Fig. 1. The cutting unit is the major component element of longwall mining equipment, which taking on the task of cutting. The cutting unit is composed of the cutting motor, the gear transmission system, and cutting drum, and the cutting drum interact with coal seam directly. The transmission chain (as shown in Fig. 2) of the cutting unit is so long with several parallelaxis gears and one or two stage of planetary gears. Moreover, the coal seam is usually inhomogenous with hard inclusion and rock intercalation, which induces the following features including heavy load, strong impact and large fluctuation to the cutting drum load. The cutting drum load transfer to the gear transmission system by cutting drum, which will lead to damage of the gears and bearings in transmission system as shown in Fig. 3, and make the coal wining machine stop working. So the gear transmission chain is one of weakest part of the cutting unit. In order to solve this problem, and improve the operational reliability of the coal wining machine, a short driving system (as shown in Fig. 4) in cutting unit was proposed in this paper.
Fig. 1The longwall shearer
Fig. 2The transmission chain of cutting unit
As shown in Fig. 4, the short driving system type of arrangement is also sometimes called a ‘power confluence design’, which is consisting of a multiple motor driving unit (MMDU) and a herringbone planetary gear unit (HPGU). The MMDU is composed of electric motor and parallelaxis transmission gear set and is driven by three or more electric motors via parallelaxis transmission gear set. Then, the cutting drum is driven by the herringbone planetary gear set, which connect with parallelaxis transmission gear set (including three pinions and a bull wheel). Compared with the gear transmission system used in the existing longwall shearer cutting unit, the short driving system can offer the following advantages:
a) A reduced number of gear stages;
b) Lower energy losses;
c) Increased reliability because of separate drive paths;
d) Fewer gears and bearings;
e) Lower noise levels from gear meshes;
f) Lower overall drive system weight.
In MMDU, the input torque applied should theoretically be shared equally by each transmission path, i.e. each path should carry 1/n of the total input torque [2, 3]. On the one hand, a certain manufacturing errors and assembly variations of components are present; on the other hand, traditional control theory is difficult to meet the performance requirements of control system, especially multiple equipment coordinated system [4]. A synchronous control system is usually used for the MMDU and guarantees the pinions which connect with motor have the same velocity or torque. In practice, the synchronous control system has control accuracy errors existence (namely, the electric motor has synchronization error existence between each other), so it’s hard to make the torque or speed of each electric motor complete identical. By this token, the load sharing performance is very important for the dynamic behavior and failure life of a gear system [5]. The coupling effect between the synchronous error of electric motor and the gears manufacturing eccentricity error or transfer errors etc. extrinsic and intrinsic factor may lead to wear and failure of a gear system because of bad load sharing performance. As a consequence, study the load sharing characteristic of the MMDU is necessary.
Fig. 3The damage bearing and gear
a) The damage bearing
b) The damage gear
The study of load sharing characteristic of gear transmission system for many years and improve load sharing characteristic of gear transmission system is beneficial to slow the wear and tear of gear and enhance the life of gear. The multimotorgear system can be considered as the extension of multipleparalleltransmissions, which have been studied with different approaches, based mainly on gear dynamics [6]. Jing Wei et al. [79] developed a dynamic model of multigear drive in a tunnel boring machine and the effects of inertia on the driving mechanism and cutter head are considered. The study demonstrates that the loadsharing coefficients at different rotating cutter head speeds and input torques are determined and are found to reach up to 1.21.3. A load sharing index was proposed according to the minimum and maximum of the Floquet multipliers and the simplified dynamic model of an example gear system driven by two pinions simultaneously is presented by Haidong Yu [5] et al., in which mechanical parameters such as meshing frequencies, bearing stiffness of pinions, mounting locations of pinions, and gear configuration coefficients are involved. The correlations of the load sharing behavior of every mesh gear pair and the mechanical parameters are discussed. Shaoping Bao [10] et al. developed a nonlinear dynamic model of the multimotorgear system and in which nonlinearity of gear meshing stiffness, gear geometric and assembly errors, and the motor characteristics are duely considered. A companion study to develop a method to investigate the effects of machining errors, assembly, and tooth modifications on loading capacity, loadsharing ratio, and transmission error of a pair of spur gears by using specially developed finite element method software by Li [11]. Singh and N. Alfred [1214] have studied static loadsharing performance of planetary transmission with finite element method, then Singh has provided a physical explanation for the basic mechanism causing the unequal loadsharing phenomenon. Ligata et al. [15] have presented a simplified discrete model to predict load sharing among the planets of a planetary gear set with planet carrier position errors and proposed a translational representation of the torsional system that includes any number of planets positioned at any spacing configuration. Ajit Bodas and Ahmet Kahraman [2] have mainly considered the effect of manufacturing errors on the static loadsharing behavior of planetary gear sets and proposed three parameters of the loadsharing coefficient and static loadsharing coefficient describing the loadsharing behavior of the planetary gear trains. Abousleiman and Velex [16] have presented a model which enables the simulation of the threedimensional dynamic behavior of planetary/epicyclic spur and helical gears. Dynamic loadsharing behavior and loadsharing coefficient of star gear trains with effect of each levels connection stiffness and star gear eccentric errors have been analyzed by Fang et al. [17].
However, current research on the loadsharing characteristic of MMDU takes only the nonlinear coupling factors between the bull wheel and multiple pinions into account, less research has considered the influences of the electric motor intrinsic mechanical property and synchronization error. In this study, firstly, a translationtorsional dynamic model for the parallelaxis transmission gear set of the MMDU is proposed, which is connected to the electric motor model for loadsharing characteristic analysis by the angular displacement. Next, an electromechanical dynamic model for the MMDU is constructed including the electric motor, the parallelaxis transmission gear set and the cutting drum. Then, the loadsharing characteristics are simulated when the shock loads are acted on the cutting drum and the MMDU. By taking into account the electric motor intrinsic mechanical property and synchronization error, the dynamic loadsharing characteristics of the MMDU are analyzed based on the electromechanical dynamic simulation, the influence of manufacturing eccentricity error and transfer error of gears on the loadsharing characteristics are also investigated in order to compare with the electric motor synchronization error. It provides the date support in the gear design and synchronization control of the MMDU. At last, the study contents are reviewed and conclusions are briefly concluded.
Fig. 4The short driving system
2. Model of the electric motor
The steadystate model of electric motor is used in the electromechanical dynamic analysis. The steadystate model of driving motor, namely, the intrinsic mechanical characteristic of electric motor reflects the relationship of the output rotation speed and the output torque of electric motor. The intrinsic mechanical characteristic of electric motor [18] is derived by steadystate equivalent electric circuit of electric motor as shown in Fig. 5. The intrinsic mechanical characteristic of electric motor is derived as Eqs. (1) and (2):
Combined with Eqs. (1) and (2), electromagnetic torque can be found as Eq. (3):
where ${T}_{e}$ is the electromagnetic torque of electric motor; ${P}_{e}$ is the power of electric motor; ${I}_{2}^{\text{'}}$ is equivalent rotor current value; ${m}_{1}$ is phase number; ${\mathrm{\Omega}}_{s}$ is synchronous angular speed; ${U}_{\varphi}$ is phase voltage; ${R}_{1}$ is stator resistance; ${R}_{F}$ is equivalent rotor resistance; $j{X}_{1}$ is stator leakage reactance; $j{X}_{F}$ is equivalent rotor leakage reactance; $j{X}_{M}$ is the magnetizing inductance; $s$ is slip ratio.
Fig. 5The equivalent circuit of this induction motor
3. Dynamic model of parallelaxis transmission gear set
The dynamic mathematic model of the parallelaxis transmission gear set is established by using centralized parameter method. The translationtorsional equivalent dynamic model of the parallelaxis transmission gear set is shown in Fig. 6.
Fig. 6Equivalent dynamic model for the translationtorsion of parallelaxis transmission gear set
Fig. 6 illustrated the equivalent dynamic model of the parallelaxis transmission gear set; each component of the parallelaxis transmission gear set is regarded as a rigid body. The elastic deformation of support shafts, bearings and gears are represented by equivalent spring stiffness and damping among which ${k}_{pi}\left(t\right)$ represents the timevarying mesh stiffness and ${c}_{pi}\left(t\right)$ ($i=$1,…, $N$ and $N$ is the number of pinions) represents the timevarying mesh damping between pinions and bull wheel. ${k}_{xpi}$, ${k}_{xbw}$ and ${k}_{ypi}$, ${k}_{ybw}$ ($i=$1,…, $N$ and $N$ is the number of pinions) represent support equivalent spring stiffness of pinions and bull wheel in $x$direction and $y$direction, respectively, ${c}_{xpi}$, ${c}_{xbw}$ and ${c}_{ypi}$, ${c}_{ybw}$ ($i=$1,…, $N$ and $N$ is the number of pinions) represent support equivalent spring damping of pinions and bull wheel in $x$direction and $y$direction, respectively. ${r}_{bpi}$ and ${r}_{bbw}$ respectively represent the base circle radius of the pinions and bull wheel. ${\phi}_{n}$ is the position angle of the theoretical center of $n$th pinion, where ${\phi}_{n}=2\pi \left(n1\right)/N$.
The system in Fig. 6 has (3$N+$3) degree of freedom, where $N$ is the number of pinions. The generalized vibration displacement is $X=\left({x}_{pi},{y}_{pi},{\theta}_{pi,}{x}_{bw},{y}_{bw},{\theta}_{bw}\right)$, where ${x}_{pi}$, ${x}_{bw}$, ${y}_{pi}$, ${y}_{bw}$ and ${\theta}_{pi}$, ${\theta}_{bw}$ donate $x$direction, $y$direction translational displacement and angular displacement of pinions and bull wheel, respectively.
Supposing the relative dynamic displacements along the meshing line in each pinion and bull wheel are ${\delta}_{i}$ ($i=$1,…, $N$ and $N$ is the number of pinions), the gear timevarying stiffness ${k}_{pi}\left(t\right)$ and the dynamic meshing forces ${F}_{pi}$ can be calculated according to [19, 20]. The differential equation of equivalent mathematic model is:
where ${m}_{pi}$, ${m}_{bw}$ and ${J}_{pi}$, ${J}_{bw}$ are the mass and moments of inertia of the pinions and bull wheel, respectively.
The load distribution between pinions and bull wheel are affected by the errors, which are manufacturing eccentric error of pinions ${E}_{pi}$ and manufacturing eccentric error of bull wheel ${E}_{bw}$, transmission error of pinions ${A}_{pi}$ and transmission error of bull wheel ${A}_{bw}$. The equivalent accumulated meshing error caused by the errors which are mentioned in the above, along with the meshing line can be defined as follows according to [8, 21]:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{A}_{pi}\mathrm{s}\mathrm{i}\mathrm{n}\left({\omega}_{pi}t+{\beta}_{pi}\right)+{A}_{bw}\mathrm{s}\mathrm{i}\mathrm{n}\left({\omega}_{bw}t+{\beta}_{bw}\right),$
where ${\alpha}_{t}$ is the meshing angle and ${\phi}_{i}$ is the azimuth angle of each pinion. In this paper, the MMDU is driven by three electric motors and the motors are evenly distributed, ${\phi}_{i}=2\pi \left(i1\right)/3$. $\gamma $ and $\beta $ are the phase angles of manufacturing eccentric error and transmission error. ${\omega}_{pi}$ and ${\omega}_{bw}$ represent the angular velocity of pinion and bull wheel, respectively.
4. Electromechanical dynamic model of the short driving system
The electromechanical dynamic model of the short driving system as shown in Fig. 7, which is composed of electric motors, the parallelaxis transmission gear set, herringbone planetary gear set and the cutting drum. In this paper, the loadsharing characteristic of the MMDU of the short driving system will be studied, so an electromechanical dynamic model of the MMDU is constructed including the electric motors, the parallelaxis transmission gear set and the cutting drum, where the planetary set is simplified to a rotary inertia ${J}_{p}$ and a transmission ratio ${i}_{p}$ due to the loadsharing coefficient of first stage planetary system has no effect on the loadsharing coefficient of the other stage for multistage planetary system [9]. In Fig. 7, ${T}_{mi}$, ${J}_{mi}$ and ${\theta}_{mi}$ ($i=$1, 2, 3) are the torque acting on the electric motor rotor, the rotary inertia and angular displacement of the electric motor rotor. ${T}_{d}$, ${J}_{d}$ and ${\theta}_{d}$ are the torque acting on the cutting drum, the rotary inertia and angular displacement of the cutting drum. ${K}_{mi}$ and ${K}_{d}$ are the connecting stiffness between the electric motors and the pinions, ${C}_{mi}$ and ${C}_{d}$ are the connecting damping between the electric motors and the pinions, respectively.
Fig. 7Electromechanical dynamic model of short driving system
The torque acting on the electric motor rotor (${T}_{mi}$) can be calculated by Eq. (3), the torque acting on each pinion (${T}_{pi}$) and bull wheel (${T}_{bw}$) in Eq. (4) can be obtained by Eqs. (6) and (7). The differential equation of the cutting drum is derived as Eq. (8). The differential equation of the electric motor is derived as Eq. (9). The electromechanical model can be constructed with Eqs. (4), (8) and (9):
5. Solution of loadsharing coefficient
In what follows, the loadsharing characteristic of the MMDU will be studied which use the developed electromechanical dynamic model, with interests in the influences of the electric motor synchronization error, manufacturing eccentric error and transmission error. Here, the loadsharing coefficient is a measurement index for load uneven level in each pinion at a certain time. The loadsharing coefficient in each tooth frequency cycle can be calculated as [8]:
where ${\lambda}_{pij}$ is the loadsharing coefficient in each tooth frequency cycle; $i=1N$ and $j=1n$ ($n$ is the number of tooth frequency cycle in one system period). The dynamic meshing forces ${F}_{pi}$ can be calculated according to [19, 20].
Define the loadsharing coefficient in one system period as:
where ${\kappa}_{pi}$ is the loadsharing coefficient in one system period.
6. Simulation and analysis
The coal seam is usually inhomogenous with hard inclusion and rock intercalation, which induces the shock load on the driving system, so the loadsharing characteristics of system are simulated under shock load with different electric motor synchronization error, different manufacturing eccentric error and different transmission error. The shock load is a transient process when coal winning machine is working, so in order to investigate the effect of each error which mentions above on the loadsharing characteristic of MMDU under the action of the shock load, a square signal is used to simulate the load mutation process which is acting on the cutting drum, as shown in Fig. 8. The shock load starts at 0.5 s and ends at 0.6 s with the amplitude is 70 KN. the parameters of the multiple motor driving system are shown in Table 1.
Table 1Parameters of the multiple motor driving unit
Parameters  
Electric motor  ${P}_{e}=$100 Kw, ${\mathrm{\Omega}}_{s}=$ 154.4616 rad/s, ${U}_{\mathrm{\Phi}}=$ 658.1793 V, ${R}_{1}=$0.1749 Ω; $j{X}_{1}=$1.503 Ω, $j{X}_{F}=$1.503 Ω, ${R}_{F}=$0.0408 Ω, ${m}_{1}=$ 3, ${J}_{mi}=$7.23 Kg∙m^{2}, the electric motor is connected in $Y$ 
Parallelaxis gear set  The teeth number of pinions and bull wheel are 16 and 99, respectively; the module is 4 mm 
Supporting and torsional stiffness  ${k}_{xpi}={k}_{ypi}=$10^{9} N/m, ${k}_{mi}=$ 10^{8} N/m, ${k}_{d}=$10^{9} N/m 
Fig. 8Torque acting on the cutting drum
The supporting damping, torsional damping and meshing damping can be obtained as follows:
where ${\zeta}_{S}$ is the flexural damping ratio of material, the value is 0.030.1; ${\zeta}_{T}$ is the torsional damping ratio of material, the value is 0.0050.075; ${\zeta}_{V}$ is the meshing damping ratio, the value is 0.030.17.
6.1. Effect of synchronization error on loadsharing coefficient and dynamic load
The loadsharing coefficients are not same under different synchronization error of electric motor when master slave control method is used in MMDU of the short driving system, which imply that the torque acting on each pinion is different. Synchronization error is defined as the ratio of the difference value of the output torque of master motor and slave motor to the output torque of the master motor. A value of 0 represents the output torque of slave motor equal to the output torque of master motor, a value greater or less than 0 represents the output torque of slave motor less or higher than master motor. The loadsharing coefficients of MMDU can be obtained using RungeKutta integration method based on electromechanical dynamic model (as shown in Fig. 7).
Fig. 9Loadsharing coefficient and dynamic meshing force with synchronization error of 10 %
a) Loadsharing coefficient
b) Dynamic meshing force
Loadsharing coefficient and dynamic meshing force with synchronization error of 10 % is shown as Fig. 9 when the steadystate electric motor model is used and the initial rotating speed of the electric motor is set as the rated speed. According to Fig. 9, the loadsharing coefficient of pinion 1 is greater than the loadsharing coefficient of pinion 2 and 3, $\delta $ represents the difference value of loadsharing coefficient between pinion 1 and other two pinions; but the loadsharing coefficient of pinion 2 and 3 is nearly equal, as shown in Fig. 9(a). This is to say, the dynamic meshing force acting on pinion 1 is greater than the dynamic meshing force acting on pinion 2 and 3, and the dynamic meshing force acting on pinion 2 and 3 is almost equal, which can be proved by calculating the dynamic meshing force ( as shown in Fig. 9(b)) of each pinion. The reason for this phenomenon is that the output torque of the electric motor which connects with pinions are different as shown in Fig. 10; it is obvious that the output torque of electric motor 1 (master motor) is greater than electric motor 2 and 3 (salve motor), and the output torque of electric motor 2 and 3 is almost equal.
Fig. 10The output torque of electric motor with synchronization error of 10 %
As shown in Fig. 9(a), 0.50.6 s shows the loadsharing coefficient under the shock load, the loadsharing coefficient of each pinion has a nonlinear fluctuation process before stable. This is because the difference value of the output torque between electric motor 1 and electric motor 2 and 3 has a decrease and increase process (as the area A, B, C and D shown in Fig. 10) during 0.50.6 s. The maximum loadsharing coefficient of pinions is 1.453 at 0.6 s, because the difference value of the output torque between the electric motor 1 and the electric motor 2 and 3 reach the maximum value. And a similar process can be observed after the shock load; this is also because the difference value of the output torque between electric motor 1 and electric motor 2 and 3 has a decrease and increase process during 0.60.7 s. As shown in Fig. 9(b), 0.50.6 s shows the dynamic meshing force variable process under the shock load. The dynamic meshing force begins to increase immediately at the beginning of the shock load, and does not reach the steady state until the end of the shock load. The dynamic meshing force begins to decrease promptly and reaches the steady state quickly at the end of the shock load. The similar rules also can be obtained when the synchronization error of electric motor changes from 0 to 50 %, and need not be repeated here.
Fig. 11The influence of manufacturing eccentric error on loadsharing coefficient and dynamic meshing force as shown in a) and b), separately, and the influence of transmission error on loadsharing coefficient and dynamic meshing force as shown in c) and d), respectively
a) Loadsharing coefficient
b) Dynamic meshing force
c) Loadsharing coefficient
d) Dynamic meshing force
6.2. Effect of manufacturing eccentric error/transmission error on loadsharing coefficient and dynamic load
In order to compare the influence of manufacturing eccentric error/transmission error on loadsharing coefficient and dynamic load with the influence of electric motor synchronization error on loadsharing coefficient and dynamic load, only the manufacturing eccentric error/transmission error of pinions is assumed to change in the case where the manufacturing eccentric errors/transmission errors for other components are assumed to be constant as 0 μm and the electric motor synchronization error is 0, as given in Eq. (5). Fig. 11 represents the variation of the MMDU loadsharing coefficients and dynamic meshing force with the manufacturing eccentric error/transmission error of pinions is 20 μm, that case will be analyzed as follow.
It can be seen from Fig. 11(a) the loadsharing coefficient appears cyclical change under the function of manufacturing eccentric error, and the maximum loadsharing coefficient of pinion 1 is slightly greater than other two pinions under the steady state, which means the dynamic meshing force of the meshing pair 1 is approximately above the dynamic meshing force of other two meshing pairs as shown in Fig. 11(b). The fluctuation amplitude of the loadsharing coefficient decrease quickly at the beginning of the shock load; and the fluctuation amplitude of the loadsharing coefficient increase immediately at the end of the shock load. The variation of loadsharing coefficient is mainly caused by the sudden increase of dynamic meshing force according to Eq. (10). The dynamic meshing force shows the opposite variation rules compare with the loadsharing coefficient, as shown in Fig. 11(b). Namely, the dynamic meshing force begins to increase immediately at the beginning of the shock load, and does not reach the steady state until the end of the shock load. The dynamic meshing force begins to decrease promptly and reaches the steady state quickly at the end of the shock load.
The similar phenomenon also can be seen in Fig. 11(c) and (d) under the function of transmission error, but the maximum loadsharing coefficient of pinion 3 is lower than other two pinions before the shock load and greater than other two pinions after the shock load. This is because the output torque of the electric motor 3 is greater than the output torque of the electric motor 1 and 2 after the shock load.
Fig. 12The relationship between the RMS of loadsharing coefficient and different errors
a) The relationship between the RMS of loadsharing coefficient and synchronization error
b) The relationship between the RMS of loadsharing coefficient and eccentric error
c) The relationship between the RMS of loadsharing coefficient and transmission error
Fig. 12 shows the relationship between the root mean square (RMS) of loadsharing coefficient of system and the three errors (the electric motor synchronization error, manufacturing eccentric error and transmission error) as discussed above. Here, the RMS of loadsharing coefficient will be used to investigate the effect degree of the electric motor synchronization error, manufacturing eccentric error and transmission error of pinions on the loadsharing characteristic of the MMDU.
From Fig. 12(a), it can be observed that the RMS of loadsharing coefficient of system obviously increases with the synchronization error of electric motor goes up. In other words, as the synchronization error of electric motor increases, the load distribution situation among different pinions will become worse. From Fig. 12(b) and (c) respectively, it can be found that the RMS of loadsharing coefficient of system increases with the manufacturing eccentric error and transmission error go up, but the amplification of RMS is not obvious compared with the synchronization error of electric motor. Meanwhile, comparing Fig. 12(a), (b) and (c) can be found that the influence of the electric motor synchronization error on the RMS of loadsharing coefficient is greater than the influence of the other two errors (manufacturing eccentric error and transmission error) on the RMS of loadsharing coefficient, and the influence of the transmission error on the RMS of loadsharing coefficient is greater than the influence of the manufacturing eccentric error on the RMS of loadsharing coefficient. So it is concluded that synchronization error of electric motor has an significant effect on the loadsharing characteristic of system, decreasing the synchronization error of electric motor can improve the loadsharing characteristic of system as shown in Fig. 12(a); and the synchronization error must be controlled under 5 % to ensure the loadsharing coefficient of the MMDU under 1.1. Decreasing the manufacturing eccentric error and transmission error also can improve the loadsharing characteristic of system as shown in Fig. 12(b) and (c), but the effect is unobvious compared with reducing the synchronization error of electric motor.
7. Conclusions
An electromechanical dynamic model for MMDU has been established to study the effect of the synchronization error of electric motor, manufacturing eccentric error and transmission error of pinions on loadsharing performance and dynamic load. Then the main conclusions are:
1) The shock load has less effect on the loadsharing characteristic of MMDU with only considering the synchronization error of electric motor, or manufacturing eccentric error of pinions, or transmission error of pinions. But it has significant effect on dynamic meshing force of gear meshing pair. The dynamic meshing force begins to increase immediately at the beginning of the shock load, and does not reach the steady state until the end of the shock load. The dynamic meshing force begins to decrease promptly and reaches the steady state quickly at the end of the shock load.
2) The RMS of loadsharing coefficient increases with the errors go up, the load distribution situation among different pinions will become worse. And the influence of the electric motor synchronization error on the loadsharing characteristic of MMDU is greater than the influence of the transmission error on the loadsharing characteristic of MMDU, and the influence of the manufacturing eccentric error on the loadsharing characteristic of MMDU is lower than the influence of the transmission error on the loadsharing characteristic of MMDU.
3) The electric motor synchronization error has significant effect on the loadsharing characteristic of MMDU. The loadsharing coefficient evidently increases with the electric motor synchronization error go up, decreasing the electric motor synchronization error can improve the loadsharing characteristic of MMDU, and the synchronization error must be controlled fewer than 5 % to ensure the loadsharing coefficient of the MMDU under 1.1.
References

Seyed Hadi Hoseinie, Mohammda Ataie, Reza Khalookakaei, Uday Kumar Reliability modeling of hydraulic system of drum shearer machine. Journal of Coal Science and Engineering. Vol. 17, Issue 4, 2011, p. 450456.

Bodas A., Kahraman A. Influence of carrier and gear manufacturing errors on the static load sharing behavior of planetary gear sets. Jsme International Journal Series CMechanical Systems Machine Elements and Manufacturing, Vol. 47, Issue 3, 2004, p. 908915.

Singh A., Kahraman A., Ligata H. Internal gear strains and load sharing in planetary transmissions: Model and experiments. Journal of Mechanical Design, Vol. 130, 2008, p. 0726027.

Tan C., Si L., Zhou X., Wang Z., Wang K. A cooperative control method for fully mechanized mining machines based on fuzzy logic theory and neural networks. Advances in Mechanical Engineering, 2014.

Yu H., Eberhard P., Zhao Y., Wang H. Sharing behavior of load transmission on gear pair systems actuated by parallel arrangements of multiple pinions. Journal of Mechanism and Machine Theory, Vol. 65, 2013, p. 5870.

Özgüven H. Nevzat, Houser D. R. Mathematical models used in gear dynamics – a review. Journal of Sound and Vibration, Vol. 121, Issue 3, 1988, p. 383411.

Wei J., Sun Q., Sun W., Cai J., Zeng J. Dynamic analysis and loadsharing characteristic of multiple pinion drives in tunnel boring machine. Journal of Mechanical Science and Technology, Vol. 27, Issue 5, 2013, p. 13851392.

Wei J., Sun Q., Sun W., Ding X., Tu W., Wang Q. Loadsharing characteristic of multiple pinions driving in tunneling boring machine. Chinese Journal of Mechanical Engineering, Vol. 26, Issue 3, 2013, p. 532540.

Sun W., Li X., Wei J., Zhang A., Ding X., Hu X. A study on loadsharing structure of multistage planetary transmission system. Journal of Mechanical Science and Technology, Vol. 29, Issue 4, 2015, p. 15011511.

Shaoping Bai, Stephance Caro A nonlinear motorgear model and its application to shareloading analysis of wind turbine yawing mechanisms. Mutibody Dynamics, Eccomas Thematic Conference, 2011, p. 110.

Li S. Effects of machining errors, assembly errors and tooth modifications on loading capacity, loadsharing ratio and transmission error of a pair of spur gears. Journal of Mechanism and Machine Theory, Vol. 42, Issue 6, 2007, p. 698726.

Singh A. Load sharing behavior in epicyclic gears: Physical explanation and generalized formulation. Journal of Mechanism and Machine Theory, Vol. 45, Issue 3, 2010, p. 511530.

Singh A. Application of a system level model to study the planetary load sharing behavior. Journal of Mechanical Design, Vol. 127, Issue 3, 2005, p. 469476.

Montestruc A. N. Influence of planet pin stiffness on load sharing in planetary gear drives. Journal of Mechanical Design, Vol. 133, 2011, p. 0145011.

Ligata H., Kahraman A., Singh A. A closedform planet load sharing formulation for planetary gear sets using a translational analogy. Journal of Mechanical Design, Vol. 131, 2009, p. 0210072.

Abousleiman V., Velex P. A hybrid 3D finite element/lumped parameter model for quasistatic and dynamic analyses of planetary/epicyclic gear sets. Journal of Mechanism and Machine Theory, Vol. 41, Issue 6, 2006, p. 725748.

Fang Z. D., Shen Y. W., Huang Z. D. The dynamic behavior of star gearing with three branches. Journal of Aerospace, Vol. 17, Issue 7, 1990, p. 725748.

Chapman S. Electric Machinery Fundamentals. Tata McGrawHill Education, 2005.

Liu C., Qin D., Liao Y. Dynamic model of variable speed process for herringbone gears including friction calculated by variable friction coefficient. Journal of Mechanical Design, Vol. 136, 2014, p. 0410064.

Liu C., Qin D., Lim T. C., Liao Y. Dynamic characteristics of the herringbone planetary gear set during the variable speed process. Journal of Sound and Vibration, Vol. 333, Issue 24, 2014, p. 64986515.

Parker R. G. A physical explanation for the effectiveness of planet phasing to suppress planetary gear vibration. Journal of Sound and Vibration, Vol. 236, Issue 4, 2000, p. 561573.
About this article
The authors would like to acknowledge the support and contribution from the State Key Lab of Mechanical Transmission, Chongqing University, China. This research was funded by the National Major Basic Research Program of China (973 Program, Grant No. 2014CB046304) and the Fundamental Research Funds for Central Universities (No. 106112015CDJXY110001).