#### Abstract

Effect of dynamic backlash and rotational speed is investigated on the six-degree-of-freedom model of the gear-bearing system with the time-varying 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, time-varying 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 time-periodic 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 two-stage 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 third-order polynomial and oscillatory different equations of the system were solved using the multiple scale method. Wei S et al. [4] proposed a single-mesh 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 time-varying parameters on gear-bearing system or gear-rotor-bearing 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 time-varying 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 one-stage 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 one-stage 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 time-varying 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 6-DOF gear meshing model, in which the time-varying 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 Runge-Kutta 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 gear-bearing transmission system with dynamic center distance and backlash. In this paper, a 6-DOF dynamic model of the gear-bearing 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 multi-degree-of-freedom 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 Gear-Bearing System

##### 2.1. The Lumped Mass of the System

Figure 1 is set to express simplified 6-DOF 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 time-varying 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 high-frequency 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 gear-bearing 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 5-DOF 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 gear-bearing 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 amplitude-frequency 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, T-periodic motion, multi-periodic 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 multi-degree-of-freedom 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 multi-degree-of-freedom 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 T-periodic 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 quasi-periodic motion at and to -periodic motion in the range of . Conclusively, through 2-periodic 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 T-periodic motion, through chaos motion at , to the quasi-periodic in the range of . Finally, the system behavior is chaos motion again, by 2-periodic bifurcation in the range of . Thus, it can be seen that the type of system motion is simplified by dynamic backlash, involving 1-periodic motion and chaos motion, but chaos motion is more complex and unpredictable than nT-periodic 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 gear-bearing transmission system is established as the key to explore the nonlinear responses. In detail, the effects of the constant and dynamic backlash on gear-bearing 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 T-periodic motion, through chaos motion, to qusia-periodic 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).