The Influences of Gradual Wears and Bearing Clearance of Gear Transmission on Dynamic Responses

Abstract

Gears are important components of the transmission system. Tooth wear and bearing clearance are significant factors affecting the dynamics of the gear system. In order to reveal the effects of gradual wears and bearing clearance on the gear system dynamics, a six-degrees-of-freedom bending-torsion coupled model of gear-rotor-bearing which considers surface wear, bearing clearance and backlash is established. The Rung-Kutta method is used to solve the nonlinear dynamic system, and the dynamic responses of the system are obtained. The results show that the time-varying mesh stiffness decreases with the tooth surface from the unworn phase to severe wear phase. At the same time, the change of the mesh stiffness in the double-tooth mesh area and single-tooth area are different. Moreover, the amplitude of the X-displacement, Y-displacement and relative gear mesh displacement will be enlarged slightly with the increase of wear depth or bearing clearance. By analyzing variation tendency in the frequency domain, the different order harmonics show the different change characteristic with the variation of the wear phases or bearing clearances. This study provides a theoretical basis for improving the transmission performance and the selection of the bearing clearances in the gear system.

1. Introduction

Geared rotor bearing systems are commonly employed as power transmission mechanisms, and gear system dynamics has received significant attention in recent years. The multi-clearances, including bearing clearance and backlash, is an important factor affecting the dynamics of the gear system. In the process of gear meshing, due to the need for lubrication and manufacture errors, there will be a backlash between the conjugate teeth pair. The presence of the backlash causes an impact between the conjugate teeth pair, which affects the stability of the gear transmission [1]. The bearing clearance affects the number of rolling elements which can bear the load. Meanwhile, due to the presence of the bearing clearance, the gear center distance and pressure angle will change accordingly. Thus, the bearing clearance has an important impact on the dynamic characteristics of the system. In addition, tooth wear is inevitable, resulting in the tooth profile deviating from the ideal tooth profile. The gradual wears affect the backlash and the meshing stiffness [2,3,4], and also affect the friction between conjugate teeth pair. Gear tooth wear occurs over relatively long periods of service, severe wear can lead to an occurrence of other types of gear failures such as pitting and broken teeth [5]. It can be seen that consideration of these factors, which incorporate tooth wear and multi-clearances, has important implications for the dynamics of the gear system.

In the process of gear system dynamics research, scholars have gradually noticed the significant influence of backlash and bearing clearance on system dynamics. Researchers have made a lot of meaningful achievements in the study of system dynamics of bearing clearance and backlash from theory and experiment. The backlash and the time-varying stiffness gear model is proposed by Kahraman [6], and then the dynamic response characteristics of the gear system are analyzed. Chen et al. [7] derived the connection between the gear central distance error and backlash, and defined the dynamic backlash. After that, they established a multi-degree of freedom gear system model, and the influence of dynamic backlash on the system was studied. In order to study the influence of backlash and time-varying mesh stiffness on the dynamics of gear system, a set of test equipment with adjustable backlash was established by Cui et al. [8]. The gear system model is improved by Gou et al. [9], and the influence of the backlash, bearing clearance, flash temperature and mesh stiffness on the nonlinear dynamics is obtained. Xu et al. [10] explored the influence of the backlash and bearing clearance on the characteristics of the gear transmission system by establishing an eight-degrees-of-freedom dynamic model.

Tooth wear changes the profile of gears, which will inevitably affect the system dynamics, Zhang et al. [11] restructured the wear tooth profile, and obtained the dynamic meshing force, the result indicated tooth wear enlarge the peak of dynamic meshing force. The elastic dynamics model of gear teeth with two degrees of freedom was established by Wojnarowski et al. [12]. Later, they analyzed and studied the effects of deformation and wear on the dynamics of spur gears. Ding et al. [13] established a dynamic model considering the factors such as tooth surface wear and comprehensive errors, and studied the influence of surface wear on the dynamics of the gear system. For the ideal involute spur gear, Feng et al. [14] studied the calculation method of the gear meshing stiffness when the tooth surface is worn. Wang et al. [15] used the Archard model and Weber-Banaschek model to calculate the cumulative wear of the tooth surface and the time-varying meshing stiffness of the wear gear pair, and they proved that the gear system has complex dynamic characteristics due to the surface wear. Liu et al. [16] replaced the effects of wear with a modulated mesh impact and introduced them into the system dynamics model to derive the effect of wear on the dynamics of the gear system.

Scholars have done a lot of research on the dynamics of gear systems with factors such as multi-clearances and wear, and have achieved many meaningful results. During the gear transmission process, the tooth wear increases continuously. The wear affects mesh stiffness continuously, thereby changing the dynamic characteristics of the gear system. The current research has not comprehensively considered the effects of wear and multi-clearances on the dynamics of the gear system. At the same time, there is little research on the influence of bearing clearance on gear center distance. In view of this, firstly, the dynamic model of six-degrees-of-freedom single stage spur gear transmission system considering wear and multi-clearances is established in this paper. Then, the time-varying mesh stiffness corresponding to different wear depths are obtained, and thus the tooth wear is introduced into the system dynamics model. Therefore, a more realistic system model is obtained. Based on this model, the effects of gradual wears and bearing clearance on the dynamics of the gear system are analyzed by numerical simulation.

2. Model Description

2.1. Gear Dynamic Model

Figure 1 shows the six-degrees-of-freedom bending-torsion coupled model of the gear-rotor-bearing. The specific vibration form of the drive shaft is not considered. y_p, x_p and θ_p represent the vertical displacement, transverse displacement and torsional displacement of the driving gear, y_g, x_g and θ_g represent the vertical displacement, transverse displacement and torsional displacement of the driven gear. The subscripts p and g represent driving gear and driven gear, respectively. Other parameters are listed in

Table 1. The parameters of the dynamic model.

Nomenclature	Description
mp, mg	gear masses
zp, zg	number of tooth
α	the pressure angle
m	gear modulus
Rp, Rg	base circle radius
kh	time-varying mesh stiffness
ch	viscous meshing damping
kjk (j = p,g; k = x,y)	bearing stiffness
cjk (j = p,g; k = x,y)	bearing viscous damping coefficients
Ip, Ig	gear inertias
bjk (j = p,g; k = x,y)	bearing clearance
b	backlash
Tp, Tg	input/output torque

. The gear mesh is described by a non-linear displacement function and viscous damping. The friction force at the mesh point is assumed to be negligible. Further, it is assumed that input torque and output torque are constant.

2.2. Tooth Wear Model

In the process of gear meshing, the wear of the tooth surface makes the tooth profile deviate from the original tooth profile and causes the impact. As mentioned before, gear wear causes the change of gear mesh stiffness, and affects the dynamic characteristics of the gear system.

According to the literature [15] and Equation (1), the wear depth at each point of the gear tooth profile can be obtained.

$$h=4I_h\mu nt\sqrt{\frac{2\times 9.55\times {10}^6{T}_p}{\pi \psi m^2{z}_{p}n\cos \alpha }\frac{\frac{1-{\nu }_1^2}{E_1}+\frac{1-{\nu }_2^2}{E_2}}{\frac{1}{{\rho }_1}+\frac{1}{{\rho }_2}}}$$

(1)

where μ is the gear slide coefficient; I_h is the coefficient of wear characteristic, E₁ and E₂ are modulus of elasticity of driving gear and driven gear. ν₁ and ν₂ are the Poisson ratio of the driving gear and driven gear, m is the gear module, α is the pressure angle in pitch circle, T_p is the input torque, ψ is the width coefficient, t is the gear working time, n is the rotation speed, the parameters of gear and bearing are listed in

Table 2. Parameters of gear and bearing.

Parameter	Value	Parameter	Value
$z_p$	$30$	$z_g$	$30$
$m \left(\mathrm{mm}\right)$	$5$	$R_{b1 }\left(\mathrm{mm}\right)$	$70.47$
$\alpha$	$20°$	$R_{b2} \left(\mathrm{mm}\right)$	$70.47$
$m_{p }\left(\mathrm{Kg}\right)$	$3.75$	$m_g \left(\mathrm{Kg}\right)$	$3.75$
${\nu }_1$,${\nu }_2$	$0.3$	$n/\left(\mathrm{r}/\min \right)$	$1200$
$k_{px}$,$k_{py}$(N/m)	2.6 × 108	$k_{gx}$,$k_{gy}$(N/m)	2.6 × 108

.

Effect of Wear on Mesh Stiffness

The contact ratio of the gear is greater than 1; therefore, there are different logarithmic teeth alternately meshing during the gear meshing process, resulting in periodic variation in tooth meshing stiffness according to the Weber stiffness model, as shown in Figure 2. Assume that the gear tooth is the cantilever beam, the effective length of the cantilever beam is L_e, the gear tooth is divided into a series of infinitesimal elements along the X axis. The infinitesimal element is represented by i.

The gear deformation is divided into three parts: (a): bending deformation and shear deformation of gear tooth; (b): deformation produced by the tooth root (c): local deformation caused by contact stress.

δ₁ is the tooth deformation, it can be obtained by Equation (2).

$${\delta }_1=\frac{F_N}{E}\left|\begin{array}{l}{\cos }^2{w}_x\left[\frac{L_j^3+3L_j^2{S}_{ij}+3L_j{S}_{ij}^2}{3I_i}\right]-\cos w_x\sin w_x\left[\frac{L_j^2{Y}_j+2L_j{S}_{ij}}{2I_i}\right]\\ +{\cos }^2{w}_x\left[\frac{12\left(1+\nu \right)L_j}{5A_i}\right]+{\sin }^2{w}_x\left(\frac{L_j}{A_i}\right)\end{array}\right|$$

(2)

δ₂ is deformation produced by the tooth root, it can be obtained by Equation (3).

$${\delta }_2=\frac{F_N{\cos }^2{w}_x}{Eb}\left(1-{\nu }^2\right)\left|5.306{\left(\frac{L_f}{2Y_M}\right)}^2+2\left(\frac{1-\nu -2{\nu }^2}{1-{\nu }^2}\right)\frac{L_f}{2Y_M}+1.534\left(1+\frac{0.4167{\tan }^2{w}_x}{1+\nu }\right)\right|$$

(3)

$$L_f=X_j-X_M-Y_j\tan w_x$$

(4)

δ₃ is the local deformation caused by contact stress, it can be obtained by Equation (5).

$${\delta }_3=\frac{1.275F_N^{0.9}}{E^{0.9}{b}^{0.8}}$$

(5)

The meshing stiffness k_h can be obtained by Equation (6).

$$k_h=\frac{F_N}{{\left({\delta }_1+{\delta }_2\right)}_1+{\delta }_3+{\left({\delta }_1+{\delta }_2\right)}_2}$$

(6)

The ${\left({\delta }_1+{\delta }_2\right)}_1$ is the deformation of the driving gear and ${\left({\delta }_1+{\delta }_2\right)}_2$ is the deformation of the driven gear.

Here, F_N is the load. The L_j is the thickness of infinitesimal element, ω_x is the load angle, and A_i is the area of the infinitesimal element, S_ij is the distance from the infinitesimal element to the loading points. More detailed parameters are given in the reference [17].

The gradual wear causes the change of the involute gear tooth profile, introduce the wear depth on the tooth surface into Equations (2), (3) and (5). Finally, the gear time-varying meshing stiffness corresponding to different wear depths can be obtained by the Equation (6).

2.3. Bearing Clearance Model

2.3.1. Bearing Clearance Function

The bearing clearance is the gap between the bearing rolling element with the inner and outer ring housing of the bearing. The size of the bearing clearance affects the number of rolling elements which bear the load and will have a great influence on the response of the whole gear system.

According to Figure 3, the bearing clearance is processed into a non-linear displacement function as follows.

$$f\left(y_p\right)=\{\begin{array}{c}{y}_p-b_{py}\\ 0\\ y_p+b_{py}\end{array}\begin{array}{c}{y}_p>b_{py}\\ \left|y_p\right|\le b_{py}\\ y_p<-b_{py}\end{array}$$

(7)

$$f\left(x_p\right)=\{\begin{array}{c}{x}_p-b_{px}\\ 0\\ x_p+b_{px}\end{array}\begin{array}{c}{x}_p>b_{px}\\ \left|x_p\right|\le b_{px}\\ x_p<-b_{px}\end{array}$$

(8)

$$f\left(y_g\right)=\{\begin{array}{c}{y}_g-b_{gy}\\ 0\\ y_g+b_{gy}\end{array}\begin{array}{c}{y}_g>b_{gy}\\ \left|y_g\right|\le b_{gy}\\ y_g<-b_{gy}\end{array}$$

(9)

$$f\left(x_g\right)=\{\begin{array}{c}{x}_g-b_{gx}\\ 0\\ x_g+b_{gx}\end{array}\begin{array}{c}{x}_g>b_{gx}\\ \left|x_g\right|\le b_{gx}\\ x_g<-b_{gx}\end{array}$$

(10)

2.3.2. Influence of Bearing Clearance on Center Distance

Due to the presence of the bearing radial clearance, radial runout occurs during the operation of the gear system, which causes the change of the center distance. Thus, it is necessary to establish the relationship between the bearing clearance and the center distance. Figure 4 is the model of the bearing clearance and center distance.

BC indicates bearing clearance. C is the theoretical center distance and C′ is the center distance which considering the influence of the bearing clearance. F_x and F_y represent the contact force in X-direction and Y-direction on driving gear, respectively. F′_x and F′_y represent the contact force in X-direction and Y-direction on driven gear, respectively.

We can obtain the relationship between bearing clearance and center distance as follows:

$$C^{\prime }=\sqrt{{\left(C+2F_x/k_b+BC\right)}^2+{\left(2F_y/k_b+BC\right)}^2}$$

(11)

where k_b is the rolling bearing stiffness, and it is assumed to be constant.

2.3.3. The Model of the Center Distance and Backlash

The backlash is used to describe the gap of the teeth along the mesh line. Theoretically, the backlash is 0. However, due to the need of lubrication and manufacture errors, there will be a backlash between the conjugate teeth pair, we define this gap as the initial backlash b₀. Actually, the change in the gear center distance will cause a corresponding change in the backlash. As we all know, the backlash is a nonlinear factor affecting the dynamic behavior of the system. So, the relationship between the center distance and the backlash is particularly important. Figure 5 is the model of the center distance and backlash. The P is the mesh point. R₁ and R₂ are pitch circle.

In the process of gear mesh, due to presence of the initial backlash b₀, the point N does not coincide with the point M. According to the geometric relationship, the actual backlash can be expressed as:

$$b=\overline{PN}-\overline{PM}$$

(12)

A formula was derived to calculate the backlash b:

$$b=(\frac{2\pi R_1}{z_p}-s_p-s_g)\mathrm{c}\mathrm{o}\mathrm{s}(\alpha {}^{\prime })+b_0$$

(13)

where s_p and s_g are pitch tooth thickness of driving and driven gear. α${}^{\prime }$ is the actual pressure angle.

$${\alpha }^{\prime }=\mathrm{a}\mathrm{r}\mathrm{c}\mathrm{c}\mathrm{o}\mathrm{s}(\frac{C}{C^{\prime }}\mathrm{c}\mathrm{o}\mathrm{s}(\alpha ))$$

(14)

According to the Equations (13) and (14), the relationship between the backlash and the actual center distance can be expressed as:

$$b=2C^{\prime }\cos (\alpha )\left[inv({\alpha }^{\prime })-inv(\alpha )\right]+b_0$$

(15)

$$inv(\alpha )=\tan (\alpha )-\alpha$$

(16)

2.4. Equations of Motion

According to the gear dynamic model in the Section 2.1, we can obtain the equation of motion as follows:

$$\{\begin{array}{l}{m}_p{\ddot{y}}_p+c_{py}{\dot{y}}_p+k_{py}f\left(y_p\right)=-F_m\cos \alpha \\ m_g{\ddot{y}}_g+c_{gy}{\dot{y}}_g+k_{gy}f\left(y_g\right)=F_m\cos \alpha \\ m_p{\ddot{x}}_p+c_{px}{\dot{x}}_p+k_{px}f\left(x_p\right)=-F_m\sin \alpha \\ m_g{\ddot{x}}_g+c_{gx}{\dot{x}}_g+k_{gx}f\left(x_g\right)=F_m\sin \alpha \\ J_p{\ddot{\theta }}_p=T_p-R_p{F}_m\\ J_g{\ddot{\theta }}_g=-T_p+R_p{F}_m\end{array}$$

(17)

The static transfer error is not considered, so the relative gear mesh displacement can be defined as:

$$\delta =R_p{\theta }_p-R_g{\theta }_g+x_p\sin \alpha +y_p\cos \alpha -x_g\sin \alpha -y_g\cos \alpha$$

(18)

where T_P/T_g is the input/ output torque of the gear dynamic system, F_m is the dynamic gear force, it can be obtained by Equation (19)

$$F_m=c_h\dot{\delta }+k_{h}f\left(\delta \right)$$

(19)

f(δ) is the backlash function of the displacement, it can be rewritten as:

$$f(\delta )=\{\begin{array}{c}\delta -b\\ 0\\ \delta +b\end{array}\begin{array}{c}\delta >b\\ \left|\delta \right|\le b\\ \delta <b\end{array}$$

(20)

3. Numerical Simulation and Discussion

3.1. Influence of Gradual Wears on System Dynamics

This section is to investigate the influence of gradual wears on gear dynamics. The following sections will detail the response of the gear system under the effect of gradual wears.

3.1.1. Gear Mesh Stiffness Under Different Wear Phases

Figure 6 shows the time-varying mesh stiffness curves under the unworn, mild wear phase and severe wear phase, respectively.

The horizontal axis represents time and the vertical axis represents the mesh stiffness value. Figure 6 shows the meshing stiffness of gears in a cycle under different bearing clearances. In Figure 6, A and C represent the double tooth mesh area respectively, and B represents the single tooth mesh area. Overall, the increase of the tooth surface wear depth reduces the mesh stiffness. At the beginning of the A-segment (double-tooth mesh area near the tooth tip), the change of the tooth mesh stiffness is most obvious under the different wear depths. From the A-segment to the B-segment (single-tooth mesh area near the pitch point), the change of tooth mesh stiffness gradually becomes weak under the different wear depths. From the B-segment to the C-segment (double-tooth mesh area near the tooth root), the change of the tooth mesh stiffness gradually becomes apparent under the different wear depths, which is consistent with the reference [18]. The reason for this may be because the wear depth of the single-tooth mesh area is tiny during the gear mesh process, which is much smaller than the wear depth of the double-tooth mesh area, so the stiffness does not change significantly in the single-tooth mesh area.

3.1.2. Translational Displacement

Figure 7a depicts the Y-displacement under different wear phases. The amplitude of the displacement will be enlarged slightly with the increase of wear depth which is mainly caused by the deterioration of gear surface. There is a change in Y-displacement induced by the wear phases. The change in the Y-displacement is shown in Figure 7b. To investigate the variation of the Y-displacement, the first four harmonics of the meshing frequency are illustrated in Figure 8, and fm = 600 Hz (below is the same). As can be seen, the first order, second order and fourth order harmonics are enlarged with the increase of the wear depth and change linearly. However, the third harmonic is descending.

It can be obtained the same trend in X-displacement from the Figure 9 and Figure 10.

3.1.3. Relative Displacement

Figure 11a depicted the relative displacement under different wear phases. The amplitude of the relative displacement will be enlarged slightly with the increase of wear depth. The spectrum of relative displacement is shown in Figure 11b. To investigate the variation of the relative displacement, the first four harmonics of the meshing frequency are illustrated in Figure 12. As can be seen, the first order, second order and fourth order harmonics are enlarged with the increase of the wear depth. However, the third harmonic is descending for the wear phases from unworn to mild wear, and then increasing for the wear phases from mild wear to severe wear.

3.1.4. Gear Meshing Force

Figure 13a illustrated the time-varying gear meshing force in time domain. It can be obtained from Figure 13b that the amplitudes of eighth to ninth harmonics are much higher than that of other harmonics in the frequency band in which the signals are modulated by the resonance frequency of the gears. To investigate the variation of the gear meshing force, the first four harmonics of the meshing frequency are illustrated in Figure 14. As can be seen, the first order, third order and fourth order harmonics are descending for the wear phases from unworn to mild wear, and then increasing for the wear phases from mild wear to severe wear. However, the third harmonic is enlarged with the increase of the wear depth.

3.2. Influence of Bearing Clearances on System Dynamics

3.2.1. Gear Parameters and Gear Mesh Stiffness under Different Bearing Clearances

Figure 15 depicted the variation of gear parameters under different bearing clearances. With the increase of the bearing clearance, the actual gear center distance will increase linearly, as shown in Figure 15a. Due to the change of the gear center distance, the gear pitch circle will be increased; however, the base circle and addendum circle will remain unchanged, as these circle are machined dimensions. The pressure angle is related to the base circle and pitch circle, and the contact line length is related to the pitch circle and addendum circle. Figure 15b,c shown the variation of the pressure angle and contact line length under different bearing clearances. The reduction of the contact line length causes the decrease in contact ratio as shown in Figure 15d. The greater the contact ratio, the more stable the gear system runs.

Gear time-varying meshing stiffness is an important factor affecting the dynamics of gear system. The change of the bearing clearance leads to the shift of gear meshing position, which further leads to the change of gear meshing stiffness. Figure 16 shows the meshing stiffness of gears in a cycle under different bearing clearances. In Figure 17, A and C represent the double tooth mesh area respectively, and B represents the single tooth mesh area. The time-varying meshing stiffness will increase gradually with the increase of bearing clearance. At the same time, in the zoomed-in figure we can see that the double tooth mesh area was reduced, and the single tooth mesh area increased. The reason for this is that the increase of bearing clearance leads to the decrease of contact ratio, which means that the double tooth mesh area is reduced. In order to explain the change of meshing stiffness of gears with bearing clearance more clearly, the average meshing stiffness of gears in a period with different clearance is calculated as shown in Figure 17. It can be seen that the average meshing stiffness decreases linearly with the increase of bearing clearance.

3.2.2. Translational Displacement

Figure 18a depicted the Y-displacement under different bearing clearances. The larger the bearing clearance is, the more shift is generated. More detail on displacement in time domain is depicted in the zoomed-in figure. The spectrum of Y-displacement is shown in Figure 18b. The first four harmonics of the Y-displacement are illustrated in Figure 19. As can be seen, the first order harmonic and second order harmonic are enlarged with the increase of the wear depth and change linearly. The third order harmonic is linearly descending. However, the fourth harmonic presents a nonlinear tendency with the increase of bearing clearance.

The same trend can be obtained in X-displacement from Figure 20 and Figure 21.

3.2.3. Relative Displacement

Figure 22a depicted the relative displacement under different bearing clearances. It can be seen that the displacement curve moves upwards as a whole. More detail about displacement in the time domain is depicted in the zoomed-in figure. The spectrum of relative displacement is shown in Figure 22b. The first four harmonics of the Y-displacement are illustrated in Figure 23. As can be seen, the first order harmonic and third order harmonic are enlarged with the increase of the bearing clearance. The second order harmonic is descending. However, the fourth harmonic presents a nonlinear tendency with the increase of bearing clearance.

3.2.4. Gear Meshing Force

Figure 24a illustrated the time-varying gear meshing force in time domain. It can be obtained from Figure 24b that the amplitudes of eighth to tenth harmonics are much higher than that of other harmonics in the frequency band in which the signals are modulated by the resonance frequency of the gears. To investigate the variation of the gear meshing force, the first ten harmonics of the meshing frequency are illustrated in Figure 25. As can be seen, the amplitude variations of the first four harmonics are nonlinear, which is caused by the nonlinearity effect of the bearing dynamics as presented in many researches. The first order harmonics descend when the bearing clearance increases, and then they increase. The second order, fifth order, eighth order, ninth order and tenth harmonics descend when the bearing clearance increases. However, the third order, fourth order and sixth order harmonics are enlarged with the increase of the bearing clearance.

4. Conclusions

(1)

The change of bearing clearance leads to the change of the actual center distance of gears, and the change of the center distance further affects the backlash. The relationship of the bearing clearance and backlash is deduced. It can be obtained that the backlash will increase linearly under static conditions.

(2)

The time-varying meshing stiffness of gear meshing under different wear phases is calculated. The time-varying meshing stiffness will increase gradually with the increase of wear depth. However, the mesh stiffness does not change significantly in the single-tooth mesh area. The translational displacement, relative displacement and dynamic meshing force will increase with the increase of the wear depth. This shows that the increase of the wear depth has negative effects on the vibration of the gear system.

(3)

The time-varying meshing stiffness of gear meshing under different bearing clearances is calculated. The time-varying meshing stiffness will increase gradually with the increase of bearing clearance, the double-tooth mesh area will reduce and single-tooth mesh area will increase. Also, the gear mesh parameters, like contact ratio and pressure angle, will change as the bearing clearance increases. In the meantime, the change of the bearing clearance affects the dynamic response of the system. The translational displacement, relative displacement and dynamic meshing force will increase with the increase of the bearing clearance. This shows that the increase of the bearing clearance has negative effects on the vibration of the gear system.