Research Article  Open Access
Dynamic Characteristics of Spur Gear Pair with Dynamic Center Distance and Backlash
Abstract
Effect of dynamic backlash and rotational speed is investigated on the sixdegreeoffreedom model of the gearbearing system with the timevarying meshing stiffness. The relationship between dynamic backlash and center distance can be defined clearly. The nonlinear differential equations of the model are solved by the Newmark method. The results show that system amplitude increases in the wake of increasing rotational speed. After reaching a certain rotational speed, the system jumps from periodic motion to chaos motion, and the effective amplitude is changed violently. Comparing the dynamic backlash with fixed backlash, the amplitude of the dynamic backlash is augmented and the frequency components are diversified. The vibration displacement is enlarged by the dynamic backlash and the chaotic behavior of the system becomes complex with increasing rotational speed. The numerical results provide a useful reference source for engineers to select rotational speed section for steady running.
1. Introduction
Gear transmission system is called as one of the most important mechanisms used to transmit power and motion in the modern machinery industry, whose mechanism and complex nonlinear phenomenon have been the focus. The lifetime of gear transmission system is reduced by a series of nonlinear factors and nonlinear excitation, including dynamic backlash, timevarying meshing stiffness, dynamic transmission error, and gear eccentricity. Therefore, it is necessarily analyzed for the dynamic behaviors of the gear transmission system and effect of the parameters on the system. In order to greatly enhance performance, reliability, and dynamic stability of the gear system, the effects of dynamic parameters and excitation have been the focus on research.
Plentiful research achievements have been accomplished by many outstanding scientists. S Theodossiades and S Natsiavas [1] researched dynamic characteristics of gear pair system, including backlash and timeperiodic stiffness. Their theory is used to guide later scientists to investigate dynamic characteristics of the gear pair system. L Walha et al. [2] presented a twostage spur gear system, in which periodically changing stiffness and fixed backlash were both taken into consideration. They indicated that amplitude jumps were generated by fixed backlash. Hamed Moradi and Hassan Salarieh [3] investigated dynamic response of gear system with nonlinear backlash, which is more realistic than constant backlash. The nonlinear backlash function was expressed by a smooth thirdorder polynomial and oscillatory different equations of the system were solved using the multiple scale method. Wei S et al. [4] proposed a singlemesh gear system to investigate effect of periodic meshing stiffness and backlash nonlinearity based on improved IHBM method.
There scholars have made contributions to research of gear transmission system with nonlinear factors. However, in order to be more tally with the actual situation, more and more scholars developed dynamic factors and the timevarying parameters on gearbearing system or gearrotorbearing system. Cao Z et al. [5] studied dynamic responses in the planetary gear system, in which gear eccentricity was taken into consideration and observed meshing stiffness changes due to the timevarying center distance and gear eccentricity. Lu JW et al. [6] established a dynamic model of the gear transmission system with the stochastic perturbation of parameters. They indicated that motion of the system became more complex due to increasing perturbation intensity. But the unstable motion of the system transformed the steady periodic motion when the perturbation intensity is considerably great. The extended research of Lu JW et al. [7] established onestage gear system model, in which stochastic backlash was included. The relationship between stability of system and random backlash was discussed. The dynamic characteristics of the onestage gear system were influenced by random backlash, and chaotic behaviors of the system could be enlarged due to random backlash. Researching the perturbation of the gear system parameters was made among scholars, and the foundation is provided for further research.
However, some facts are ignored. Some scientists studied considering more factors. Chen Q [8] researched the gear system with fractal theory. The fractal theory was employed to express the dynamic instead of constant backlash. It was clearly generated that the fractal backlash was more reasonable form than fixed backlash and random backlash. The research of the gear system with the fractal characteristics is used to guide the design, manufacture, and choice of the accuracy grade of the gear tooth surface. To further explore dynamic backlash, the timevarying mesh stiffness of a spur gear pair was changed with changing the center distance error, which was developed by Yong L et al. [9]. Chen SY et al. [10] raised a 6DOF gear meshing model, in which the timevarying mesh stiffness, dynamic backlash, and effect of friction were included. The variation of dynamic meshing force became complex with increasing rotational speed. Xiang L et al. [11] studied the gear transmission model with dynamic friction and dynamic backlash based on the RungeKutta numerical method and observed effect of backlash and gear eccentricity under different rotational speeds.
Although some dynamic factors of the gear system are discussed, few researches have investigated the effect of the rotational speed on a dynamic model of the gearbearing transmission system with dynamic center distance and backlash. In this paper, a 6DOF dynamic model of the gearbearing transmission system is established. The coupled relationship between the dynamic backlash and the vibration responses of gear system is researched by analyzing effect of the rotational speed and dynamic backlash variation. The differential equations of motion are solved using the Newmarkβ method and the simulation results including 3D frequency diagram and bifurcation figure are described. The dynamic backlash on the system with the center distance error is clearly defined to further develop the gear system running mechanism. And we will focus on effect of the rotational frequency on the multidegreeoffreedom dynamic gear system. The research provides a reference for selecting appropriate rotational frequency interval with gear system. It is practically significant to ensure the stability of the gear transmission.
2. Dynamic Model of the GearBearing System
2.1. The Lumped Mass of the System
Figure 1 is set to express simplified 6DOF lumped mass model in the paper. Figure 1 presents schematic diagram of the dynamic system with gears and bearings.
The nonlinear model of the gear system with gear inertia and , gear masses and , and base circle radius and as shown in Figure 1 is considered here. The symbol (j=p,g) is the dynamic angular displacements of the gears. The bearings are modeled by elements with viscous damping coefficients c_{j} (j=p, g) and linear springs k_{j} (j=p, g). The timevarying meshing stiffness is expressed as [12]where is the averaged tooth stiffness and is variant part of the meshing stiffness. The stiffness is shown in Figure 2. Tm is the gear meshing period. The highfrequency internal excitation caused by the static transmission error [13] is written aswhere (t) is the static transmission error, is averaged static transmission error, and is the variant part of the static transmission error. The relative displacement between the pinion and the gear is The dynamic meshing force acting on the meshing points iswhere is the meshing viscous damping coefficient. is the backlash function which will be defined later in the paper.
2.2. The Dynamic Backlash
Figure 3 presents simplified illustration of the dynamic center distance with variation of gear eccentricity. Real circles , and dotted circles , are, respectively, original position of the gears and the position of the gear system of generating vibration. The new gear axle center coordinates are, respectively, and . The actual center distance and the theoretical center distance are denoted by and a_{0} as described in Figure 3. The instantaneous distance between the pinion and the gear is expressed as [14]The dynamic pressure angle is defined as [15]where and are the base radius of driving gear and driven gear. The position angle is defined asTheoretically, the backlash should be zero because the gears are designed without the real situation of the gears running. Actually, in order to avoid seizing after the frictional heating expansion of the gear system and make sure to form lubricating film between the tooth profiles, a backlash must be considered in the gear design. If the center distance is changed, instantaneous backlash will be altered. The graphic about the tooth flanks with backlash is shown in Figure 4. The dynamic backlash is expressed aswhere is the dynamic backlash and is the backlash. The backlash function of the nonlinear dynamic system can be written as [13]
2.3. Equation of Motion
The equations of motion of the gearbearing system can be expressed aswhere (j=p,g) is the base radius of the pinion and gear, respectively. A dimensionless form of Equations (10) is obtained by letting ), , , , , , , , , , , , )/b_{0}, /b_{0}, and , where is the meshing frequency. The 5DOF dimensionless equations of the gear motion are concluded aswhere Δ is the dimensionless relative displacement expressed as
3. Results and Discussion
The parameters of gear transmission are listed in Table 1. In this section, in order to explore the nonlinear responses of the gearbearing transmission system, Eq. (11) with dimensionless equations can be solved using the Newmark method. Meshing frequency is regarded as the key parameter, which is controlled. In the following analysis, the displacement and amplitude from dimensionless equation will be depicted as shown in 3D frequency diagram, bifurcation spectrum, Poincaré map, and so on.

3.1. Dynamic Backlash and Dynamic Meshing Angle
For the gear system with dynamic backlash, it is first for dynamic backlash and dynamic meshing angle to make comparison with different as shown in Figures 5 and 6 because the frequency is very significant and a lot of brilliant scientists were devoted to investigating the nonlinear dynamic characteristics of the gear transmission system.
(a)
(b)
(c)
(a)
(b)
(c)
In Figures 5(a), 5(b), and 5(c), the dynamic backlash is steadily on the system with the frequency . The dynamic backlash presents serious fluctuation when the frequency is 2. With the frequency reaching 3, this fluctuation is decreased. The amplitude is dimensionless.
The dynamic meshing angle of the gear system is shown in Figures 6(a), 6(b), and 6(c) at different frequency severally, including , , and . The characteristics of dynamic meshing angle are similar to the dynamic backlash. To further explore the dynamic backlash and the dynamic meshing angle, we will research the global characteristics of the backlash and angle by using the amplitudefrequency response curve of the dynamic backlash and dynamic meshing angle with respect to the frequency as shown in Figures 7(a) and 7(b), respectively. The amplitude of the backlash and the degree of the meshing angle are both steady values in the range of . The increasing and fluctuation present in the range of and values of the backlash and the angle are both shown as fluctuation decrease and small value in the range of .
(a)
(b)
3.2. Nonlinear Characteristics of Gear System with Constant Backlash
Firstly, for the nonlinear system with constant backlash, meshing frequency is controlled and the frequency is determined as , , and . The phase diagrams and Poincaré maps are exhibited in Figure 8. It is explicitly illuminated that the motion trail of phase plane becomes more complex with increasing the frequency . The amplitude variation of the system is showed as the tendency of inverted parabolic. As shown in Figure 8, Tperiodic motion, multiperiodic motion, and chaos motion are described with increasingly dispersing the points in Figures 8(d), 8(e), and 8(f). The phase diagrams and Poincaré maps only express a few parts of the nonlinear characteristics of the gear transmission system. In order to explore more detailed nonlinear responses, bifurcation diagram and 3D frequency spectrum are showed in the following.
(a)
(b)
(c)
(d)
(e)
(f)
3.3. Nonlinear Characteristics of Gear System with Dynamic Backlash
Firstly, for the sake of developing the nonlinear response of the multidegreeoffreedom model, comparing constant backlash model with dynamic backlash model, the different characteristics are described by bifurcation map and 3D frequency spectrum as shown in Figures 9, 10, and 11. The frequency as control parameter is the most important, which is key to explore the dynamic response of the multidegreeoffreedom system. Detailedly, the same system parameters and numerical method are applied to two different models involving constant backlash model and dynamic backlash model.
(a)
(b)
(a)
(b)
(a)
(b)
Here, the amplitude maps of vibration are showed in Figure 9. The amplitude is dimensionless value. The is regarded as varying in the range of . The amplitude value in torsional direction is expressed as the difference between the highest and the lowest amplitude values, following divided by 2. With the increase of the frequency , in Figure 9(a), it can be clearly surveyed that amplitude goes from increasing to the decreasing within the limits of , as the characteristic of amplitude of the constant backlash model at the slow meshing frequency. Subsequently, the tendency of amplitude value verges on stabilization in the range of . In addition, the amplitude variation presents an approximate linear characteristic.
Afterword, the tendency causes explosive increase in the range of , which is discrete distribution. At , the amplitude further becomes feeble discreteness, and the situation is stability in this range. In Figure 9(b), the nonlinear responses of the system are described. It is obvious that steady amplitude value is shown in the range of . Then the distribution of points is more complex and discrete than the constant backlash model in the range of . According to the above comprehensive analysis, it can be gotten for a matching rotational speed to steadily run. For the stability of the system, the appropriate speed of the system with dynamic backlash is in the range of .
In the dynamic model, 3D frequency spectrum diagrams are depicted by Figure 10. The comparison between Figures 10(a) and 10(b) shows that the major frequency amplitude increases with the increase of the frequency on the system with constant backlash, and superharmonic and subharmonic frequencies are kind of vague. On the contrary, the major frequency amplitude is more complex and irregular and the frequency component presents diversity in Figure 10(b).
The bifurcation diagram of the system of constant backlash is shown in Figure 11(a). Complex bifurcation and chaos phenomenon about dimensionless relative displacement on meshing point emerged with increasing the frequency . As the dimensionless frequency is low, the system motion is Tperiodic motion in the range of . Then, the system suddenly changes to chaos motion within the limit of . Successively, the system undergoes from periodic motion at to quasiperiodic motion at and to periodic motion in the range of . Conclusively, through 2periodic motion in the range of , the system is going to back to chaos motion again.
The bifurcation diagram of the system of dynamic backlash is described in Figure 11(b). The bifurcation behavior of the dimensionless relative displacement of mesh point is obviously reduced. The system motion is from the Tperiodic motion, through chaos motion at , to the quasiperiodic in the range of . Finally, the system behavior is chaos motion again, by 2periodic bifurcation in the range of . Thus, it can be seen that the type of system motion is simplified by dynamic backlash, involving 1periodic motion and chaos motion, but chaos motion is more complex and unpredictable than nTperiodic motion. In other words, the system behavior becomes more complex and dynamic backlash can defer the beginning of the chaos motion.
4. Conclusions
In this paper, the nonlinear dynamic model of gearbearing transmission system is established as the key to explore the nonlinear responses. In detail, the effects of the constant and dynamic backlash on gearbearing transmission system have been developed. The relationship between gear center distance and the dynamic backlash is defined considering the gear and bearing factors. To adopt by bifurcation diagram, 3D frequency spectrum, and phase diagram, the influence of the frequency on the system is explored. The results show the following:
With increasing the frequency , the amplitude of the system with dynamic backlash is heightened. In other words, it is observed that vibrational energy between gears becomes bigger. Effect of dynamic backlash on the system amplitude showed that the energy changes and becomes more and more complex.
In 3D frequency spectrum, the amplitude of the system with constant backlash is gradually rising in the wake of increase of the frequency . Nevertheless, the amplitude of the system with dynamic backlash becomes very complicated and changeable. Superharmonic and subharmonic of the response are a lot fresher on the system of dynamic backlash than the system with constant backlash.
By means of numerical simulation, bifurcation diagram, phase diagram, and Poincaré map are obtained. With the dynamic backlash, the gear system undergoes from Tperiodic motion, through chaos motion, to qusiaperiodic motion and back to periodic motion following varying the frequency . After synthetically considering conclusions and , the most appropriate rotational speed can be identified to ensure that the gear transmission system operates normally.
Data Availability
All the data is in the manuscript. If the researchers are interested in obtaining the numerical solution files, please contact email address: 1274186512@qq.com.
Conflicts of Interest
The authors declared no potential conflicts of interest with respect to the research, authorship, and publication of this article.
Acknowledgments
The authors gratefully acknowledge the financial support provided by Natural Science Foundation of China (No. 51675350).
References
 S. Theodossiades and S. Natsiavas, “Periodic and chaotic dynamics of motordriven gearpair systems with backlash,” Chaos, Solitons & Fractals, vol. 12, no. 13, pp. 2427–2440, 2001. View at: Publisher Site  Google Scholar
 L. Walha, T. Fakhfakh, and M. Haddar, “Backlash effect on dynamic analysis of a twostage spur gear system,” Journal of Failure Analysis and Prevention (JFAP), vol. 6, no. 3, pp. 60–68, 2006. View at: Publisher Site  Google Scholar
 H. Moradi and H. Salarieh, “Analisis of nonlinear oscillations in spur pairs with approximated modelling of backlash nonlinearity,” Mechanism and Machine Theory, vol. 51, pp. 14–31, 2012. View at: Google Scholar
 S. Wei, Q. Han K, and X. Dong J, “Dynami response of a singlemesh gear syatem with periodic mesh stiffness and backlash nonlinearity under uncertainty,” Nonlinear Dynamics, vol. 89, no. 1, pp. 49–60, 2017. View at: Publisher Site  Google Scholar
 Z. Cao, Y. Shao, M. Rao, and W. Yu, “Effects of the gear eccentricities on the dynamic performance of a planetary gear set,” Nonlinear Dynamics, vol. 91, pp. 1–15, 2018. View at: Google Scholar
 J. Lu, H. Chen, F. Zeng, A. F. Vakakis, and L. A. Bergman, “Influence of system parameters on dynamic behavior of gear pair with stochastic backlash,” Meccanica, vol. 49, no. 2, pp. 429–440, 2014. View at: Publisher Site  Google Scholar
 J Lu, F Zeng, J Xin et al., “Influence of Stochastic Perturbation of Parameters on Dynamic Behavior of Gear System,” Journal of Mechanical Science and Technology, vol. 25, no. 7, pp. 1667–1673, 2011. View at: Google Scholar
 Q. Chen, Y Ma, S Huang et al., “Research on gears’ dynamic performance influence by gear backlash based on fractal theory,” Applied Surface Science, vol. 313, pp. 325–332, 2014. View at: Google Scholar
 L. Yong, Natalie. Baddour, and L. Ming, “Effect of Gear Center Distance Variation on Time Varying Mesh Stiffness of a Spur Gear Pair,” Engineering Failure Analysis, vol. 75, pp. 37–53, 2017. View at: Publisher Site  Google Scholar
 Y. S. Chen, Y. J. Tang, W. C. Lou et al., “Nonlinear dynamic characteristics of geared rotor bearing system with dynamic backlash and friction,” Mechanism and Machine Theory, vol. 46, pp. 466–478, 2011. View at: Google Scholar
 L. Xiang and N. Gao, “Coupled torsion–bending dynamic analysis of gearrotorbearing system with eccentricity fluctuation,” Applied Mathematical Modelling, vol. 50, pp. 569–584, 2017. View at: Publisher Site  Google Scholar
 J. H. Kuang and Y. T. Yang, “An estimate of mesh stiffness and load sharing ratio of a spur gear pair,” in Proceedings of the International Power Transmission and Gear Conference, vol. 1, pp. 1–9, 1992. View at: Google Scholar
 A. Kahraman and R. Singh, “Interactions between timevarying mesh stiffness and clearance nonlinearities in a geared system,” Journal of Sound and Vibration, vol. 146, no. 1, pp. 135–156, 1991. View at: Publisher Site  Google Scholar
 C. Siyu, T. Jinyuan, L. Caiwang, and W. Qibo, “Nonlinear dynamic characteristics of geared rotor bearing systems with dynamic backlash and friction,” Mechanism and Machine Theory, vol. 46, no. 4, pp. 466–478, 2011. View at: Publisher Site  Google Scholar
 Y. Luo, N. Baddour, and M. Liang, “Effects of gear center distance variation on time varying mesh stiffness of a spur gear pair,” Engineering Failure Analysis, vol. 75, pp. 37–53, 2017. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2019 Jie Liu 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.