Dynamics Modeling and Experimental Investigation of Gear Mechanism with Coupled Small Clearances

Internal gear mechanism is widely used in micro-nano satellites due to its compact structure and high precision transmission. However, the vibration coupling caused by the small clearance coupling is more obvious and cannot be ignored under low speed, light load and zero gravity conditions. Based on the geometric relationship between radial clearance and backlash, a coupled model between dynamic backlash and radial clearance of internal meshing gear is established. Based on the conformal contact theory, the radial collision force model of the gear shaft and shaft sleeve considering the small clearances is established. Additionally, a multi-clearance gear rotor system test device is built to measure the vibration acceleration of the internal gear rotor system by an acceleration sensor and transmitted to the industrial computer through a signal collector for data processing. Through the comparison of simulation and experiment, the accuracy of the gear dynamics model is verified. The analysis results show that, compared with the traditional model, the calculation results of the gear mechanism model considering the small clearance coupling is closer to the experimental data.


Introduction
With the rapid development of micro-satellites, the transmission mechanism as the core component has become the focus of attention. The internal gear mechanism is widely used in the field of aerospace due to its compact structure and high precision. However, due to manufacturing and assembly errors, fatigue and wear, the clearance in the internal gear is inevitable. Under the conditions of low speed, light load and zero gravity, the vibration of the gear with small clearance coupling is more obvious and cannot be ignored. In addition, these clearances will lead to contact and impact between gear groups, thus affecting the control accuracy. With the development of high-speed and high-precision mechanical engineering, the motion accuracy and operational stability of mechanical systems are more demanding than ever before, and the effect of clearance on the dynamic characteristics of the system has gradually attracted the attention of scholars [1][2][3]. The previous models are too much simplified, and the object of study is relatively simple; therefore, these models are basically studied for the single clearance problem. However, there are multiple clearances in the configurations actually used in engineering. For example, there are four rotary hinge clearances in the crank-rocker mechanism; three rotary hinge clearances and one sliding clearance exist in the crank slider mechanism and the radial clearance of the bearing and the backlash of the gear are included in the gear system. Therefore, the study of the multi-clearance mechanism has great application value in mechanical structure design.
The dynamic research of gear mechanism considering multi-clearance started from Kahraman's nonlinear model of gear system which considered both backlash and radial clearance [4][5][6][7]. Then, Byrtus et al. analyzed the influence of nonlinear factors such as clearance and time-varying mesh stiffness on the vibration characteristics of the system in detail [8]. R Farère et al. believed that it is necessary to deal with bearing clearance with the establishment of a dynamic model of gear mechanism [9]. Farshidianfar and Anooshirvan analyzed the nonlinear characteristics of multi-clearance gear mechanism [10]. Subsequently, Fernandez and others established more comprehensive multi-clearance models considering factors on this basis [11][12][13]. However, in the existing models, the radial clearance model and the backlash model of the gear are independent of each other, ignoring the influence of the radial run-out of the gear on the position of the gear tooth meshing point and the backlash, unable to describe the coupling between radial clearance and the backlash. In fact, when the radial clearance of the bearing causes the gear to run out in the gear transmission process, the actual center distance between the gear teeth changes, and the size of the meshing point and the backlash changes accordingly. In some special application backgrounds, such as aerospace agencies, low speed, light load and zero gravity will make the coupled vibration between multiple clearances more obvious and difficult to suppress. However, gear mechanisms are widely used in aerospace structures.
Therefore, a dynamic model of internal meshing gear rotor system considering multiclearance coupling is established to further study the influence of multi-clearance coupling on the dynamic characteristics of internal meshing gear mechanism [14]. The test platform of the multi-clearance inner meshing gear rotor system was established to carry out a dynamics verification test study. Through the comparative analysis of the test data, the difference in the dynamic response between the traditional dynamic model of the multi-clearance gear mechanism and the multi-clearance model proposed in this paper is obtained.

Dynamic Model of Multi-Clearance Internal Gear Mechanism System
An internal gear mechanism is widely used in aerospace field due to its compact structure and high precision transmission. However, in the internal meshing gear mechanism, clearances exist not only in the bearing, but also in the tooth side, which will lead to radial run out of the gear and dynamic change of the backlash, and the torsional vibration at the gear place will influence the radial vibration. Therefore, in order to reflect the phenomenon of vibration caused by coupling of multiple clearances, a nonlinear dynamic model of internal gear mechanism with multiple clearances is established, as shown in Figure 1. The dynamic research of gear mechanism considering multi-clearance started from Kahraman's nonlinear model of gear system which considered both backlash and radial clearance [4][5][6][7]. Then, Byrtus et al. analyzed the influence of nonlinear factors such as clearance and time-varying mesh stiffness on the vibration characteristics of the system in detail [8]. R Farère et al. believed that it is necessary to deal with bearing clearance with the establishment of a dynamic model of gear mechanism [9]. Farshidianfar and Anooshirvan analyzed the nonlinear characteristics of multi-clearance gear mechanism [10]. Subsequently, Fernandez and others established more comprehensive multi-clearance models considering factors on this basis [11][12][13]. However, in the existing models, the radial clearance model and the backlash model of the gear are independent of each other, ignoring the influence of the radial run-out of the gear on the position of the gear tooth meshing point and the backlash, unable to describe the coupling between radial clearance and the backlash. In fact, when the radial clearance of the bearing causes the gear to run out in the gear transmission process, the actual center distance between the gear teeth changes, and the size of the meshing point and the backlash changes accordingly. In some special application backgrounds, such as aerospace agencies, low speed, light load and zero gravity will make the coupled vibration between multiple clearances more obvious and difficult to suppress. However, gear mechanisms are widely used in aerospace structures.
Therefore, a dynamic model of internal meshing gear rotor system considering multiclearance coupling is established to further study the influence of multi-clearance coupling on the dynamic characteristics of internal meshing gear mechanism [14]. The test platform of the multi-clearance inner meshing gear rotor system was established to carry out a dynamics verification test study. Through the comparative analysis of the test data, the difference in the dynamic response between the traditional dynamic model of the multi-clearance gear mechanism and the multi-clearance model proposed in this paper is obtained.

Dynamic Model of Multi-Clearance Internal Gear Mechanism System
An internal gear mechanism is widely used in aerospace field due to its compact structure and high precision transmission. However, in the internal meshing gear mechanism, clearances exist not only in the bearing, but also in the tooth side, which will lead to radial run out of the gear and dynamic change of the backlash, and the torsional vibration at the gear place will influence the radial vibration. Therefore, in order to reflect the phenomenon of vibration caused by coupling of multiple clearances, a nonlinear dynamic model of internal gear mechanism with multiple clearances is established, as shown in Figure 1.   respectively, e brj is the clearance vector of the gears, e rj and e bj are the rotor center vector and bearing center vector of the gears, respectively, R j is the pitch radius of the gears, I j and θ j are moment of inertia and angular displacement of the gears, respectively. Here, when j is p, it represents the driving gear, when j is g, it represents the driven gear. In the dynamic modeling, the dynamic gear backlash model is first established, and then the gear torsional vibration model and the radial vibration model are given. Finally, the dynamic model of the multi-clearance gear mechanism considering the radial clearance and dynamic gear backlash is given.

Multi-Clearance Coupled Model
In order to establish the radial vibration model of the gear, the bearing can be simplified as a planar shaft and shaft sleeve, and the radius difference between them is the radial clearance C ri . In Figure 1, it can be represented by the radial clearance vector e bri , which is computed with the following Equation: Further, the maximum embedding depth between the shaft to the collar is described by the clearance vector, that is, the maximum embedding amount δ imax , which can be calculated using the following relationship: Additionally, the gear rotor will collide with the shaft sleeve radially, causing the gear to run out radially, which will cause the gear backlash to dynamically change during the gear transmission process. The actual center distance A R of the gear can be expressed by the center distance vector e pg in Figure 1. The actual center distance A R can be calculated as: As shown in the formula, the actual center distance A R changes with the change of the clearance vector, and the change of A R will further cause the change of the size of the backlash. According to the geometric relationship of gear meshing, when the meshing position changes with the actual center distance, the dynamic backlash can be expressed as: where R p is the index circle radius of the driving gear, b 0 is the initial backlash of gear, Z p is the number of teeth of the driving gear, s p is the driving gear pitch tooth thickness, s g is the thickness of the driven gear pitch, α is the actual meshing angle when the main and driven gear teeth mesh, as given by: According to Equations (1), (3) and (4), the functional expression of the dynamic backlash b t is given as It can be seen from Equations (5) and (6) that in the dynamic verification test, adjusting the center distance A R , the actual meshing angle α of the gear changes and the backlash b t also changes. Additionally, the actual meshing angle α is only affected by the radial run out of the bearing, and the radial run-out amplitude can be achieved by controlling the radial radius difference between the shaft and the sleeve in the experiment.

Radial Collision Force Model and Gear Meshing Force Model
With reference to the Hertz model theory, the shaft is regarded as a rigid wedge, and the shaft sleeve is established with reference to the Winkler elastic foundation model, as shown in Figure 2, in which there is no interaction force between the springs.
According to Equations (1), (3) and (4), the functional expression of the dynamic backlash bt is given as: It can be seen from Equations (5) and (6) that in the dynamic verification test, adjusting the center distance AR, the actual meshing angle α′ of the gear changes and the backlash bt also changes. Additionally, the actual meshing angle α′ is only affected by the radial run out of the bearing, and the radial run-out amplitude can be achieved by controlling the radial radius difference between the shaft and the sleeve in the experiment.

Radial Collision Force Model and Gear Meshing Force Model
With reference to the Hertz model theory, the shaft is regarded as a rigid wedge, and the shaft sleeve is established with reference to the Winkler elastic foundation model, as shown in Figure 2, in which there is no interaction force between the springs. The Hertz model is used to establish the radial impact force model of gear bearing at a contact point, the expression is: where Fri(t) is the bearing collision force Frp(t) and Frg(t) of the driving and driven gears, Kri is the contact stiffness coefficient, n is the force index,  i is the amount of embedding at that point,  i is the relative velocity in the radial direction and Cri is the damping coefficient. According to the Equation (7), the radial force at point A can be obtained as: In the Equation (8), the contact stiffness Kri can be expressed as: where υbi and Ebi are the Poisson's ratio and elastic modulus of the sleeve, respectively, υji and Eji are the Poisson's ratio and elastic modulus of the shaft, respectively.
The damping coefficient Cri of the radial impact force can be given by: The Hertz model is used to establish the radial impact force model of gear bearing at a contact point, the expression is: where F ri (t) is the bearing collision force F rp (t) and F rg (t) of the driving and driven gears, K ri is the contact stiffness coefficient, n is the force index, δ i is the amount of embedding at that point, δ i is the relative velocity in the radial direction and C ri is the damping coefficient. According to the Equation (7), the radial force at point A can be obtained as: In the Equation (8), the contact stiffness K ri can be expressed as: where υ bi and E bi are the Poisson's ratio and elastic modulus of the sleeve, respectively, υ ji and E ji are the Poisson's ratio and elastic modulus of the shaft, respectively. The damping coefficient C ri of the radial impact force can be given by: where c ei is the coefficient of restitution, which can be used to describe the energy loss during a collision, is the relative speed at the critical collision. The radial force in y direction of j at any contact point between the shaft and the sleeve can be written as: F riy (t) = K ri δ n ij + C ri Insertion amount of any contact point between shaft and sleeve can be represented as follows: where R is the radius of the shaft. The angle between the contact boundary point of shaft and sleeve and the vertical direction of shaft axis can be expressed as: The resultant force of radial force can be calculated using: The dynamic meshing force Ft(t) between the driving and driven gears can be written as follows: where the right side of the equation is the elastic meshing force and damping force of the gear, respectively. In elastic meshing force, K t (t) is the time-varying meshing stiffness, which is a periodic function used to express the influence of the change of the number of coincident teeth on the meshing stiffness when the gear teeth are meshed. The time-varying meshing stiffness K t (t) caused by Kahraman [5] can be computed by: where k m is the average meshing stiffness, k a is the stiffness amplitude, ϕ is the initial phase and ω m is the meshing frequency, which is a function related to the number of teeth and speed.
The main parameter of the meshing damping force in Equation (15) is the nonlinear damping coefficient C t (t), which is used to describe the energy loss. The nonlinear damping coefficient C t (t) can be obtained from: where C m is the maximum damping coefficient, d is the maximum embedding amount and β is the embedding ratio, which can be expressed as: β = f g (t)/d. In Figure 1, the master and driven gears move in three directions in the coordinate system, namely translation along the x direction, translation along the y direction and rotation around the z axis. The components of the clearance vector in the x and y directions can be expressed as e brpx , e brpy , e brgx and e brgy . When considering the radial movement of the gear, the relative meshing displacement can be written as: Through the clearance function of Kahraman [15], the relative meshing relationship between the driving and driven gears can be obtained as: Entropy 2021, 23, 834 6 of 14 when the relative displacement is in the state of g t ≥ 0, the main and driven gears mesh normally, when the relative displacement is in the state of −b t < g t < 0, the main and driven gears are in a state of disengagement, when the relative displacement is in the state of g t ≤ −b t , tooth back collision occurs between the driving and driven gears.

Gear System Dynamic Mode
According to the characteristics of the multi-clearance gear mechanism, a dynamic model of the gear mechanism system considering multi-clearance coupling can be established, which can be calculated as follows: where .. θ rp and m p are the angular acceleration and mass of the driving gear, respectively, .. θ rg and mg are the angular acceleration and mass of the driven gear, respectively, T p is the driving torque and T g is the load torque, which is not added in the dynamic verification, only the dynamic characteristics of the driven gear under free motion are considered. g is the acceleration of gravity (it can also represent the acceleration under microgravity conditions). When the organization is placed horizontally, δ py are the component of the translational acceleration of the driving gear and the driven gear on the x and y axes in the global coordinate system, respectively. F rpx (t), F rpy (t), F rgx (t) and F rgy (t) are the component of the bearing collision force, respectively, F tpx (t), F tpy (t), F tgx (t) and F tgy (t) are the component of gear meshing force, respectively.

Analysis of Vibration Characteristics of Gear Rotor System
In this section, through the combination of simulation and experimental research, the non-coupled model, the coupled model and the experimental data are compared to verify the correctness of the multi-clearance coupled theoretical model. Study on vibration characteristics of gear rotor system under different working conditions was done by changing rotational speed and clearance size.
The specific structure of the test device, as shown in Figure 3, mainly includes driving gear, driving gear shaft, driven internal gear, internal gear frame, support block, base, pulley, etc. Among them, the driving gear and the driven internal gear are all involute gear, the support block is placed in the groove of the base, the internal gear and the internal gear frame are fixed together and the center hole of the internal gear frame is assembly with the cylinder on the support block. In the test, the center distance and shaft hole clearance of the main and driven gears can be adjusted by changing the support blocks with different sizes and specifications, so as to change the backlash and radial size. In the test, the driven gear was not loaded, so as to be able to simulate low speed, light load and zero gravity, and then to study the phenomenon of radial collision and tooth side collision in a free state. In the aerospace structure, the execution end (except the control moment gyro) generally has a speed below 1 rpm. Considering the deceleration capacity of the reducer, the input speed generally does not exceed 180 rpm. The horizontal arrangement of the gear shaft in most gear tests is not adopted in the test device, the arrangement of the gear perpendicular to the horizontal plane is adopted in order to equivalent the space microgravity environment, so as to ensure that the rotation plane of the gear is not affected by gravity, and the radial collision is freer. a free state. In the aerospace structure, the execution end (except the control moment gyro) generally has a speed below 1 rpm. Considering the deceleration capacity of the reducer, the input speed generally does not exceed 180 rpm. The horizontal arrangement of the gear shaft in most gear tests is not adopted in the test device, the arrangement of the gear perpendicular to the horizontal plane is adopted in order to equivalent the space microgravity environment, so as to ensure that the rotation plane of the gear is not affected by gravity, and the radial collision is freer. After the schematic design of the test device is completed, the actual design, processing and assembly problems should be considered. In addition, there are several factors to consider. First of all, it is the selection of gears, because the adjustment of gear backlash is realized by changing the actual center distance of the two gears. According to Equation (6), the proportional relationship between the increase of center distance that can be written as À A A A     and gear backlash is about 1:1. Therefore, in order to ensure that the backlash can be adjusted in a larger range, a gear with a larger modulus must be selected, so that the tooth height is longer, and the center distance will not be too large to be unable to mesh. Finally, a gear with a gear module of 4 is selected. Second, the experimental device must be able to achieve adjustable speed to study the influence of different speeds on the dynamic response of the system, so a DC servo motor is selected, and the adjustable range of speed is 0~300 rpm. Finally, the power is transmitted through the synchronous belt, and the motor is installed separately from the base to avoid the influence of the vibration of the motor itself on the experimental data collection. Considering the gear inertia and the size of the entire device, a spur gear with 20 teeth and an internal gear with 60 teeth are selected. Through repeated design and verification, a multiclearance gear rotor system test device was finally developed. The actual product is shown in Figure 4. After the schematic design of the test device is completed, the actual design, processing and assembly problems should be considered. In addition, there are several factors to consider. First of all, it is the selection of gears, because the adjustment of gear backlash is realized by changing the actual center distance of the two gears. According to Equation (6), the proportional relationship between the increase of center distance that can be written as ∆A∆A = A − A and gear backlash is about 1:1. Therefore, in order to ensure that the backlash can be adjusted in a larger range, a gear with a larger modulus must be selected, so that the tooth height is longer, and the center distance will not be too large to be unable to mesh. Finally, a gear with a gear module of 4 is selected. Second, the experimental device must be able to achieve adjustable speed to study the influence of different speeds on the dynamic response of the system, so a DC servo motor is selected, and the adjustable range of speed is 0~300 rpm. Finally, the power is transmitted through the synchronous belt, and the motor is installed separately from the base to avoid the influence of the vibration of the motor itself on the experimental data collection. Considering the gear inertia and the size of the entire device, a spur gear with 20 teeth and an internal gear with 60 teeth are selected. Through repeated design and verification, a multi-clearance gear rotor system test device was finally developed. The actual product is shown in Figure 4.  On this basis, the multi-gear rotor system test platform was built, including the multigear rotor system test device, motor driver, three-axis acceleration sensor, signal collector, host computer, etc. Among them, the motor driver is used to control the motor speed to realize the adjustable speed of the driving gear. The three-axis acceleration sensor is used to measure the vibration of the system, which is installed on the shaft of the driven gear, and the acceleration along the center line is measured emphatically. In this direction, not only the vibration frequency of the contact and collision between the shaft and the gear can be observed, but also the component of the gear meshing vibration in this direction can be observed. The acceleration sensor transmits the acceleration signal collected from the driven wheel to the upper computer through the signal collector, and the upper computer is used to further analyze the spectral characteristics. The characteristics of the test system are that it can easily adjust the size of the bearing radial clearance and backlash, and can control the driving speed, and then study the vibration characteristics of the gear rotor system under multi-clearance coupling. The main parameters of the equipment used On this basis, the multi-gear rotor system test platform was built, including the multigear rotor system test device, motor driver, three-axis acceleration sensor, signal collector, host computer, etc. Among them, the motor driver is used to control the motor speed to realize the adjustable speed of the driving gear. The three-axis acceleration sensor is used to measure the vibration of the system, which is installed on the shaft of the driven gear, and the acceleration along the center line is measured emphatically. In this direction, not only the vibration frequency of the contact and collision between the shaft and the gear can be observed, but also the component of the gear meshing vibration in this direction can be observed. The acceleration sensor transmits the acceleration signal collected from the driven wheel to the upper computer through the signal collector, and the upper computer is used to further analyze the spectral characteristics. The characteristics of the test system are that it can easily adjust the size of the bearing radial clearance and backlash, and can control the driving speed, and then study the vibration characteristics of the gear rotor system under multi-clearance coupling. The main parameters of the equipment used in the test platform are shown in Table 1. The dynamic verification test is carried out on the test platform of multi-clearance gear rotor system, and the difference of dynamic characteristics between the multi-clearance coupled model proposed in this paper and the traditional multi-clearance model is compared and analyzed based on the test data. Among them, the comparison model based on Kahraman's multi-clearance gear system dynamics model [4][5][6][7], in the radial clearance, considering the nonlinear contact stiffness and nonlinear damping coefficient, and in the gear meshing, considering the time-varying meshing stiffness and damping characteristics. However, it is considered that the backlash is a constant, which is mainly related to the machining and assembly, precision without considering the dynamic change of backlash caused by the motion deviation of the radial backlash. That is to say, the dynamic equation of the system is basically the same as that in Equation (7), but the backlash is a constant. The multi-clearance coupled model and multi-clearance non-coupled model are used to compare with the test data in the test verification process. The specific structural parameters of the gear are shown in Table 2. According to the above parameters, the corresponding numerical calculation and test are carried out, and the numerical calculation results are compared with the test results. In the test, the acceleration of the active gear is mainly controlled by the motor and has a negligible effect on the experimental results, so the vibration acceleration of the driven gear shaft is mainly measured. Therefore, the calculated data are also obtained in the numerical calculation, and then the frequency spectrum analysis is performed to facilitate comparison with the test data.

Study on Dynamic Characteristics of Gear Considering Multi-Clearance Coupling at Different Rotation Speeds
First, we studied the influence of different speeds on the vibration characteristics of the system. The bearing radial clearance c r and the initial backlash b 0 were considered to be 100 µm, the driving speed ω p of the driving gear was 60 rpm, 120 rpm and 180 rpm, respectively. The comparison between numerical calculation data and experimental data is shown in Figures 5-7. tween radial clearance and dynamic lateral clearance, which can reflect the nonlinear phenomenon caused by multi-clearance coupling in the test. Compared with the multi-clearance uncoupled model, it is more practical. Furthermore, the difference between the numerical calculation and the experimental data is analyzed and compared mainly in terms of frequency and amplitude, as shown in Tables 3 and 4.      nomenon caused by multi-clearance coupling in the test. Compared with the multi-clearance uncoupled model, it is more practical. Furthermore, the difference between the numerical calculation and the experimental data is analyzed and compared mainly in terms of frequency and amplitude, as shown in Tables 3 and 4.     ance uncoupled model, it is more practical. Furthermore, the difference between the numerical calculation and the experimental data is analyzed and compared mainly in terms of frequency and amplitude, as shown in Tables 3 and 4.  Acceleration data obtained in the experiment are transformed into frequency spectrum by Fourier transform. Vibration spectrums of driven gears at different rotational speeds are shown in Figures 5-7, respectively. Therefore, the ordinate uses the acceleration amplitude to represent the vibration characteristics of the system. According to the rotational speed and gear structure parameters, the corresponding frequencies of the three vibration peaks are, respectively; the radial vibration frequency f r between the inner gear support frame and the support block, which is mainly caused by the continuous contact between the inner hole of the inner gear support frame and the cylindrical part of the support block, gear meshing frequency f m is caused by meshing in and out of gear teeth, the gear meshing frequency is 2 f m , which is mainly caused by the existence of nonlinear factors such as backlash and radial clearance. In the numerical calculation, if the above three frequencies can be reflected in the frequency spectrum of multi-clearance coupled model and noncoupled model, which can tell the two models are correct in calculating the basic vibration law of gear. Judging from the changing trend of amplitude with speed in the above three figures, when the speed is 120 rpm in the test data, the amplitude of radial vibration and the meshing vibration and its multiplier are significantly higher than those at the other two speeds, but the frequency characteristics are still consistent with the other two sets of data, so the influence of test errors can be excluded. This phenomenon is the vibration caused by the coupling of radial clearance and the backlash, which makes the amplitude suddenly increase. From the frequency spectrum of the coupled model and the uncoupled model, it can be seen that the amplitude of the coupled model also appears a similar phenomenon, when the speed is 120 rpm, the amplitude is larger but does not appear in the uncoupled model. The amplitude of the uncoupled model increases with the increase of rotational speed, and this trend is also quite different from the experimental data.
It can be seen from the above qualitative analysis that due to the simultaneous interaction of radial clearance and backlash in the test, the system vibration caused by the coupling of radial clearance and the backlash, leads to a certain degree of small fluctuations in the frequency spectrum. From the perspective of experiment, it is shown that there is indeed a coupling phenomenon between radial clearance and backlash in the gear mechanism, which will adversely affect the vibration characteristics of the system. The multi-clearance coupled model proposed in this paper considers the coupling relationship between radial clearance and dynamic lateral clearance, which can reflect the nonlinear phenomenon caused by multi-clearance coupling in the test. Compared with the multi-clearance uncoupled model, it is more practical. Furthermore, the difference between the numerical calculation and the experimental data is analyzed and compared mainly in terms of frequency and amplitude, as shown in Tables 3 and 4.  It can be seen from Tables 3 and 4 that the error between the coupled model and the uncoupled model in the calculation of radial vibration frequency and engagement vibration frequency is smaller than that of the test data, but there is a big difference between the uncoupled model and the test in the calculation of the amplitude, especially when the rotation speed is 120 rpm, the calculation error is very large, indicating that the uncoupled model cannot accurately reflect the coupling vibration that may occur in the multi-clearance gear mechanism. The amplitude error of the coupled model is relatively small, and it can reflect the vibration caused by the coupling of radial clearance and the backlash at 120 rpm, which indicates that the coupled model can reflect the vibration characteristics of multi-clearance gear mechanism more accurately due to the consideration of the coupling relationship between radial clearance and dynamic lateral clearance and is more in line with the reality than the traditional multi-clearance model.

Study on Gear Dynamic Characteristics Considering Multi-Clearance Coupling under Different Clearance Sizes
The driving speed was set as 120 rpm to further verify the vibration characteristics of a gear rotor system with different radial clearance and backlash. By changing the size and position of the cylindrical part of the supporting block at the driven gear, the experimental data of three different clearance combinations are obtained, namely c r = 200 µm, b 0 = 100 µm; c r = 100 µm, b 0 = 200 µm; c r = 200 µm, b 0 = 200 µm. The numerical calculation values and test results of its vibration frequency spectrum are shown in Figures 8-10

Study on Gear Dynamic Characteristics Considering Multi-Clearance Coupling under Different Clearance Sizes
The driving speed was set as 120 rpm to further verify the vibration characteristics of a gear rotor system with different radial clearance and backlash. By changing the size and position of the cylindrical part of the supporting block at the driven gear, the experimental data of three different clearance combinations are obtained, namely cr = 200 μm, b0 = 100 μm; cr = 100 μm, b0 = 200 μm; cr = 200 μm, b0 = 200 μm. The numerical calculation values and test results of its vibration frequency spectrum are shown in Figures 8-10, respectively.   Figures 8-10 shows the vibration acceleration spectrum of driving gear with different clearances, which reflects the vibration characteristics of the system with different combinations of radial clearance and backlash. It can be seen from the test data of different clearance combinations that the amplitude of gear is obviously higher than that of the other two clearance combinations when the radial clearance cr = 200 μm and the initial backlash b0 = 100 μm, which indicates that the coupling vibration of the system is caused by the

Study on Gear Dynamic Characteristics Considering Multi-Clearance Coupling under Different Clearance Sizes
The driving speed was set as 120 rpm to further verify the vibration characteristics of a gear rotor system with different radial clearance and backlash. By changing the size and position of the cylindrical part of the supporting block at the driven gear, the experimental data of three different clearance combinations are obtained, namely cr = 200 μm, b0 = 100 μm; cr = 100 μm, b0 = 200 μm; cr = 200 μm, b0 = 200 μm. The numerical calculation values and test results of its vibration frequency spectrum are shown in Figures 8-10, respectively.   Figures 8-10 shows the vibration acceleration spectrum of driving gear with different clearances, which reflects the vibration characteristics of the system with different combinations of radial clearance and backlash. It can be seen from the test data of different clearance combinations that the amplitude of gear is obviously higher than that of the other two clearance combinations when the radial clearance cr = 200 μm and the initial backlash b0 = 100 μm, which indicates that the coupling vibration of the system is caused by the

Study on Gear Dynamic Characteristics Considering Multi-Clearance Coupling under Different Clearance Sizes
The driving speed was set as 120 rpm to further verify the vibration characteristics of a gear rotor system with different radial clearance and backlash. By changing the size and position of the cylindrical part of the supporting block at the driven gear, the experimental data of three different clearance combinations are obtained, namely cr = 200 μm, b0 = 100 μm; cr = 100 μm, b0 = 200 μm; cr = 200 μm, b0 = 200 μm. The numerical calculation values and test results of its vibration frequency spectrum are shown in Figures 8-10, respectively.    Figures 8-10 shows the vibration acceleration spectrum of driving gear with different clearances, which reflects the vibration characteristics of the system with different combinations of radial clearance and backlash. It can be seen from the test data of different clearance combinations that the amplitude of gear is obviously higher than that of the other two clearance combinations when the radial clearance cr = 200 μm and the initial backlash b0 = 100 μm, which indicates that the coupling vibration of the system is caused by the  Figures 8-10 shows the vibration acceleration spectrum of driving gear with different clearances, which reflects the vibration characteristics of the system with different combinations of radial clearance and backlash. It can be seen from the test data of different clearance combinations that the amplitude of gear is obviously higher than that of the other two clearance combinations when the radial clearance c r = 200 µm and the initial backlash b 0 = 100 µm, which indicates that the coupling vibration of the system is caused by the multi-clearance collision under this clearance combination. However, the amplitude of the uncoupled model varies from the experiment to some extent. The uncoupled model reflects that the vibration amplitude of the system is the largest when c r = 200 µm and b 0 = 200 µm, there is no obvious coupling vibration in the system. It can be seen that the multi-clearance coupled model is in good agreement with the experimental data in the calculation of the influence of the clearance change on the vibration characteristics of the system, which can reflect the coupling vibration phenomenon of the system when c r = 200 µm and b 0 = 100 µm.
Furthermore, the difference between numerical calculation and experimental data is analyzed quantitatively, mainly from the frequency and amplitude, as shown in Tables 5  and 6. It can be seen that the error between the calculation of radial vibration frequency and meshing vibration frequency of multi-clearance coupled model and non-coupled model and the test data is small, but the error between the calculation of amplitude of non-coupled model and the test is large, especially when the clearance size is c r = 200 µm and b 0 = 100 µm, the error is large, indicating that the non-coupled model cannot accurately reflect the coupling vibration phenomenon of multi-clearance gear mechanism. The amplitude error of the coupled model is relatively small and can reflect the vibration phenomenon caused by the coupling of radial clearance and the backlash when c r = 200 µm and b 0 = 100 µm. In summary, in the test of gear mechanism considering multiple clearances, a coupling vibration phenomenon was observed and is caused by the interaction of the radial bearing clearance and the dynamic backlash on the gear. Through the analysis of test data, it is proved that the coupled model proposed in this paper can more accurately reflect this coupling vibration phenomenon. Additionally, it is also more accurate in the calculation of frequency and amplitude. However, the coupling relationship between radial clearance and backlash was not considered in the non-coupled model, and the dynamic change of backlash was ignored, resulting in a large clearance with the test results.
Through the experimental analysis, it is found that there is still a certain error between the theoretical calculation and the experimental data. The analysis shows that the error source is mainly caused by the following reasons: (1) In the numerical calculation, the selection of contact stiffness and damping is calculated by the empirical formula, and there is a certain error with the actual device, which directly leads to the deviation between the calculation of collision force and the actual collision force that ultimately leads to the inaccurate calculation of vibration amplitude.
(2) Due to the limitation of the machining accuracy of the driven wheel shaft and sleeve, the surface roughness and the overall morphology are different from the theoretical model, which is one of the reasons for the large clearance between the vibration amplitude in the test and the theoretical calculation. Thus, machining precision of shaft can be further improved to reduce test error.

Conclusions
In this paper, dynamical model of gear mechanism with multiple clearances coupling is investigated by use of experimental test and theoretical modeling. The following conclusions can be summarized as: (1) A new dynamic model of the gear-rotor system considering small clearance coupling is established. This model considers nonlinear problems such as bearing radial clearance, dynamic backlash and time-varying meshing stiffness, which can more comprehensively reflect the influence of multi-clearance coupling effect on system dynamic characteristics.
(2) Through the dynamic test of multi-clearance gear mechanism, it is found that there exists a coupling vibration phenomenon in multi-clearance gear mechanism with bearing radial clearance and backlash, which will adversely affect the vibration characteristics of the system.
(3) By comparing the theoretical data with the experimental data, it is found that the multi-clearance coupled model proposed in this paper can well reflect the coupling vibration characteristics of the system. Compared with the traditional multi-clearance model, the calculation of vibration frequency and amplitude is in better agreement with the experimental data.
Influence of multi-clearance coupling on multi-stage gears transmission will be further studied by experiment and simulation analysis.