In modern precision machinery, gear drives are widely employed due to their stable transmission and ease of control, especially in unidirectional motion. However, the inherent backlash in gear pairs introduces return errors and degrades positioning accuracy. To mitigate these adverse effects, various anti-backlash mechanisms have been developed. Traditional solutions, such as dual-lead worm gears and dual-servo motor systems, offer high transmission ratios and smooth operation but suffer from manufacturing complexity, high cost, and maintenance difficulties. In response to these challenges, we propose a novel anti-backlash drive mechanism based on double worm gears and opposing preloaded nuts. This paper presents a comprehensive study of the nonlinear dynamic behavior of this mechanism, focusing on the role of worm gears in eliminating backlash and the influence of meshing damping and stiffness on system stability.
1. Description of the Double Worm Gear Anti-Backlash Mechanism
The proposed mechanism is specifically designed for rotary tables in machining centers, where high precision and zero backlash are critical. As illustrated below, the mechanism consists of a large gear rigidly connected to a rotary table, which meshes with two identical small gears. Each small gear is driven by a worm gear pair (worm and worm wheel) mounted on the same shaft. Two electric motors drive the worms independently. A pair of preloading nuts is threaded onto the worm shafts in opposite directions to generate an axial preload that forces the two worms apart. This preload ensures that the worm–worm wheel interfaces and the large–small gear interfaces maintain single-sided contact at all times, effectively eliminating backlash in both transmission chains.
When the rotary table requires clockwise rotation, the left tooth flank of small gear (driven by worm–worm wheel set A) contacts the right tooth flank of the large gear, transmitting torque from motor A. Conversely, counterclockwise rotation is achieved through worm–worm wheel set B, which drives the other small gear. This configuration allows the table to reverse direction without any lost motion, as the preloaded worm gears guarantee that one side of each gear pair is always engaged. The use of worm gears in this context is advantageous because worm gears inherently provide high reduction ratios, self-locking capabilities, and smooth meshing. However, the nonlinearities introduced by the worm gear geometry, such as time-varying mesh stiffness and damping, require careful dynamic analysis.
2. Dynamic Model of the Mechanism
2.1 Lumped‑Parameter Model
To analyze the nonlinear dynamic characteristics, we establish a lumped‑parameter model of the transmission system. The following assumptions are made:
- All components are considered as rigid bodies with concentrated inertia.
- Elastic deformations occur only at the gear meshes, represented by springs and dampers.
- Bearings and shafts are assumed to be sufficiently stiff, so their flexibility is neglected.
- The two small gears are identical, and the large gear is symmetric.
The model includes three rotational degrees of freedom corresponding to the angular displacements of the large gear (\(\theta_2\)) and the two small gears (\(\theta_1\) and \(\theta_3\)). The base radii are \(r_1\), \(r_2\), \(r_3\), and the moments of inertia are \(I_1\), \(I_2\), \(I_3\). The meshing damping coefficient is \(c\), and the time‑varying mesh stiffness is \(K(t)\). Static transmission error is denoted by \(\varepsilon(t)\). External torques applied to the large gear and the two small gears are \(T_2\), \(T_1\), and \(T_3\) respectively.
Applying Newton’s second law to each rotating mass yields three coupled equations:
$$
\begin{aligned}
T_1 &= I_1 \ddot{\theta}_1 + c r_1 \left[ r_1 \dot{\theta}_1 – r_2 \dot{\theta}_2 – \dot{\varepsilon}(t) \right] + r_1 K(t) \, f\!\left[ r_1\theta_1 – r_2\theta_2 – \varepsilon(t) \right] \\
– T_2 &= I_2 \ddot{\theta}_2 – c r_2 \left[ r_1 \dot{\theta}_1 – r_2 \dot{\theta}_2 – \dot{\varepsilon}(t) \right] – r_2 K(t) \, f\!\left[ r_1\theta_1 – r_2\theta_2 – \varepsilon(t) \right] \\
T_3 &= I_3 \ddot{\theta}_3 + c r_3 \left[ r_3 \dot{\theta}_3 – r_2 \dot{\theta}_2 – \dot{\varepsilon}(t) \right] + r_3 K(t) \, f\!\left[ r_3\theta_3 – r_2\theta_2 – \varepsilon(t) \right]
\end{aligned}
$$
where \(f(\cdot)\) is the gap nonlinear function that describes backlash. For conventional gears with backlash \(2b\), the function is piecewise linear:
$$
f(x) =
\begin{cases}
x – b, & x > b \\
0, & -b \le x \le b \\
x + b, & x < -b
\end{cases}
$$
Because our mechanism eliminates backlash via preloading, we set \(b = 0\) in this analysis. The dynamic transmission error is defined as \(u_d = r_1\theta_1 – r_2\theta_2\). The net relative displacement along the line of action is \(u(t) = u_d – \varepsilon(t)\). Combining the equations and introducing the equivalent mass \(m_e\) and average load \(F_m\) leads to a single second‑order differential equation:
$$
m_e \ddot{u}(t) + c \dot{u}(t) + K(t) f[u(t)] = F_m – m_e \ddot{\varepsilon}(t)
$$
where
$$
m_e = \frac{I_1 I_2}{I_1 r_2^2 + I_2 r_1^2}, \qquad F_m = \frac{T_1}{r_1} = \frac{T_2}{r_2}
$$
2.2 Time‑Varying Parameters and Excitation
The time‑varying mesh stiffness \(K(t)\) is periodic due to alternating single‑tooth and double‑tooth contact. It can be expressed as the sum of a mean value \(k_m\) and a fluctuating component approximated by a cosine function:
$$
K(t) = k_m + k_n \cos(\omega_e t + \phi_n)
$$
where \(k_n\) is the amplitude of the stiffness variation, \(\omega_e\) is the mesh frequency, and \(\phi_n\) is the phase angle (taken as zero). The static transmission error \(\varepsilon(t)\) is also periodic and can be expanded in a Fourier series:
$$
\varepsilon(t) = \sum_{n=1}^{\infty} \varepsilon_n \sin(n\omega_e t + \varphi_n)
$$
with \(\varepsilon_n\) and \(\varphi_n\) being the harmonic coefficients and phases.
2.3 Dimensionless Form
To facilitate numerical integration, we perform a dimensionless normalization. Introduce the natural frequency \(\omega_0 = \sqrt{k_m / m_e}\), a characteristic length \(l\), and the dimensionless time \(\tau = \omega_0 t\). Define the dimensionless displacement \(x(\tau) = u(t)/l\). Then the equation becomes:
$$
\ddot{x}(\tau) + 2\zeta \dot{x}(\tau) + \left[ 1 + k_a \cos(\Omega \tau) \right] f[x(\tau)] = f_0 + f_a \Omega^2 \sum_{n=1}^{\infty} \sin(n\Omega \tau + \varphi_n)
$$
where the dimensionless parameters are:
| Symbol | Definition | Meaning |
|---|---|---|
| \(\zeta\) | \(c / (2 m_e \omega_0)\) | Mesh damping ratio |
| \(k_a\) | \(k_n / k_m\) | Stiffness fluctuation amplitude ratio |
| \(\Omega\) | \(\omega_e / \omega_0\) | Dimensionless mesh frequency |
| \(f_0\) | \(F_m / (m_e l \omega_0^2)\) | Dimensionless mean load |
| \(f_a\) | \(\varepsilon_1 / l\) | Dimensionless excitation amplitude |
3. Numerical Solution and Dynamic Response
Equation (10) is a second‑order nonlinear ordinary differential equation. We solve it using the fourth‑order Runge–Kutta method implemented in MATLAB. To apply the method, we first convert it into a system of first‑order equations:
$$
\begin{aligned}
\dot{x}_1 &= x_2 \\
\dot{x}_2 &= -2\zeta x_2(\tau) – \left[ 1 + k_a \cos(\Omega\tau) \right] f[x_1(\tau)] + f_0 + f_a \Omega^2 \sum_{n=1}^{\infty} \sin(n\Omega\tau + \varphi_n)
\end{aligned}
$$
where \(x_1 = x\) and \(x_2 = \dot{x}\). The initial conditions are set to \(x(0)=0.5\) and \(\dot{x}(0)=0.5\). The gear parameters used in the simulation are listed in Table 2.
| Parameter | Value |
|---|---|
| Module | 6 mm |
| Number of teeth (small gear / large gear) | 24 / 192 |
| Face width | 80 mm |
| Mass (small gear / large gear) | 10.70 kg / 137.75 kg |
| Moment of inertia (small gear \(I_1\), \(I_3\)) | 3.255 × 10⁻² kg·m² |
| Moment of inertia (large gear \(I_2\)) | 46.101 kg·m² |
| Input torque (active gear) | 400 N·m |
| Input speed | 24 r/min |
| Equivalent mass \(m_e\) | 5.150 kg |
| Mean force \(F_m\) | 5128 N |
| Natural frequency \(\omega_0\) | 1393.466 rad/s |
| Mesh frequency \(\omega_e\) | 60.319 rad/s |
From these values, the dimensionless parameters become: \(\zeta = 0.05\), \(k_a = 0.1\), \(\Omega = 0.043\), \(f_0 = 512.8\), \(f_a = 1538.4\), and \(\varphi_n = 0\). The backlash width \(b\) is zero because the worm gears are preloaded.
The numerical results for the displacement and velocity time histories are shown in Fig. 4 (not reproduced here in HTML, but described). The displacement response initially exhibits a pronounced oscillatory transient. After about \(\tau = 100\), the transient decays and the response settles into a steady‑state harmonic oscillation with constant amplitude and period. The velocity response also shows a decreasing amplitude over time, approaching a constant amplitude. This behavior confirms that the anti‑backlash mechanism, by eliminating the gap nonlinearity through the preloaded worm gears, achieves stable and smooth transmission with low vibration. The worm gears are therefore essential in providing both zero backlash and favorable nonlinear dynamics.
4. Influence of Meshing Damping and Stiffness
4.1 Effect of Mesh Damping Ratio \(\zeta\)
To investigate how damping affects the transient and steady‑state behavior, we vary \(\zeta\) from 0.02 to 0.15 while keeping all other parameters constant. Figure 5 (described) shows the displacement time histories for four damping ratios.
| \(\zeta\) | Transient duration (approx. \(\tau\)) | Steady‑state amplitude |
|---|---|---|
| 0.02 | ~180 | Moderate |
| 0.06 | ~120 | Slightly reduced |
| 0.10 | ~80 | Further reduced |
| 0.15 | ~60 | Significantly reduced |
As the damping ratio increases, the initial oscillation decays more rapidly, and the steady‑state vibration amplitude decreases. This indicates that higher mesh damping, which can be achieved by proper lubrication or material selection in the worm gears, helps the system reach stable operation faster and reduces residual vibrations. For worm gears, the sliding contact naturally provides higher damping compared to spur gears, which is beneficial for dynamic stability.
4.2 Effect of Mesh Stiffness Variation \(k_a\)
The stiffness fluctuation amplitude ratio \(k_a\) is another critical parameter. We examine cases with \(k_a = 0\) (constant stiffness), 0.05, 0.2, and 0.4. Figure 6 (described) illustrates the displacement responses.
| \(k_a\) | Transient amplitude | Steady‑state behavior |
|---|---|---|
| 0 | Large initial oscillation | Almost zero amplitude after transient; very smooth |
| 0.05 | Moderate initial oscillation | Small harmonic amplitude; stable |
| 0.2 | Small initial oscillation | Noticeable harmonic amplitude; still periodic |
| 0.4 | Very small initial oscillation | Response becomes quasi‑periodic or chaotic; large amplitude fluctuations |
The results indicate that as \(k_a\) increases, the transient overshoot decreases, but the steady‑state vibration amplitude increases. When \(k_a\) surpasses a critical threshold (around 0.3–0.4 in this system), the response loses periodicity and enters more complex nonlinear regimes. Therefore, minimizing \(k_a\) is desirable for stable operation. Since \(k_a = k_n/k_m\), reducing the stiffness variation amplitude \(k_n\) or increasing the mean stiffness \(k_m\) will lower \(k_a\). In worm gears, the mean mesh stiffness can be enhanced by using larger module, wider face width, or stiffer materials, while the variation \(k_n\) depends on the tooth profile and contact ratio.
5. Discussion and Practical Implications
The double worm gear anti‑backlash mechanism effectively eliminates the backlash nonlinearity, resulting in a system that behaves as a linear oscillator under moderate parameter ranges. The presence of preloaded worm gears ensures that the system operates in the single‑sided contact regime, avoiding the discontinuous piecewise‑linear stiffness that causes jump phenomena and chaos in conventional gear drives. Our numerical analysis confirms that both mesh damping and stiffness play vital roles in shaping the dynamic response.
In practical design, engineers should aim for high mesh damping (e.g., through viscous lubricants or elastomeric inserts in worm gears) to accelerate transient decay and reduce vibration amplitudes. Simultaneously, increasing the average mesh stiffness (by optimizing worm gear geometry or using hardened materials) lowers the stiffness fluctuation ratio and promotes a more stable steady‑state response. The worm gears themselves can be designed with a high contact ratio to smooth the stiffness variation. For example, a multi‑tooth contact worm gear can reduce \(k_n\) significantly.
Furthermore, the preload applied by the opposing nuts must be carefully chosen. Too little preload may not eliminate backlash under all load conditions, while excessive preload increases friction and wear, potentially altering the damping and stiffness characteristics. A trade‑off exists between backlash elimination and dynamic performance, and the optimal preload should be determined experimentally.
6. Conclusion
We have presented a comprehensive study of the nonlinear dynamic characteristics of a novel double worm gear anti‑backlash drive mechanism. The main conclusions are:
- The mechanism successfully eliminates backlash by using preloaded opposing worm gears, transforming the system into a continuous nonlinearity‑free model.
- Numerical simulations based on a lumped‑parameter model show that the response transitions from transient oscillation to steady‑state harmonic motion, confirming stable and smooth transmission.
- Increasing the mesh damping ratio reduces both the transient duration and the steady‑state vibration amplitude, enhancing dynamic performance.
- Reducing the stiffness fluctuation amplitude ratio (by raising the mean stiffness or lowering the variation) suppresses steady‑state vibrations and prevents the onset of chaotic behavior.
- The use of worm gears in this anti‑backlash configuration offers inherent advantages in damping and stiffness controllability, making it a promising solution for high‑precision rotary tables and similar applications.
Future work will focus on experimental validation of the model and optimization of worm gear parameters for minimal vibration and maximum reliability.
References
- J. H. Baek, Y. K. Kwak, and S. H. Kim, “On the frequency bandwidth change of a servo system with a gear reducer due to backlash and motor input voltage,” Archive of Applied Mechanics, vol. 73, no. 5/6, pp. 367–376, 2003.
- R. Liu, G. Zhang, L. Gong, et al., “Variable inertial damper and its application in zero‑drive gear hobbing machine,” Chinese Journal of Mechanical Engineering, vol. 19, no. 3, pp. 439–441, 2006.
- Z. Yang, J. Shang, Z. Luo, et al., “Comprehensive mesh stiffness of anti‑backlash mechanism with torsion spring loaded double gears,” Journal of Mechanical Engineering, vol. 49, no. 1, pp. 23–30, 2013.
- Y. Li, W. Chen, Q. Chen, et al., “Influence of backlash on nonlinear vibration characteristics of gear pair,” Mechanical Transmission, vol. 37, no. 5, pp. 1–5, 2013.
- X. Wang, S. Wu, and X. Zhou, “Bifurcation and chaos analysis of gear transmission system with backlash nonlinearity,” Journal of Vibration and Shock, vol. 27, no. 1, pp. 53–56, 2008.
- T. Kiekbusch and I. Howard, “A common formula for the combined torsional mesh stiffness of spur gears,” in Proc. 5th Australasian Congress on Applied Mechanics, Brisbane, Australia, Dec. 2007.
- M. Hamed and S. Hassan, “Analysis of nonlinear oscillations in spur gear pairs with approximated modeling of backlash nonlinearity,” Mechanism and Machine Theory, vol. 51, pp. 14–31, 2012.
- Y. Li, T. Chen, X. Wang, et al., “Nonlinear dynamic characteristics of gear system under external dynamic excitation,” Journal of Xi’an Jiaotong University, vol. 48, no. 1, pp. 101–105, 2014.
- S. Chen, J. Tang, and C. Luo, “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.
- Y. Wu and D. Zhao, CNC Machining Equipment and Technology for Special Stone, Beijing: Science Press, 2011.
- L. Rao, “Dynamic characteristics of HTM50200 special stone machining center,” M.S. thesis, Shenyang Jianzhu University, 2012.
- X. Mu, K. Zhang, and D. Zhao, “The topological optimization design of cross beams structure of shaped stone turning‑milling CNC (HTM50200),” Advanced Materials Research, vol. 694–697, pp. 2725–2728, 2013.
- L. Wang, R. Li, T. Lin, et al., “Research on nonlinear vibration characteristics of gear system with time‑varying stiffness and backlash,” China Mechanical Engineering, vol. 14, no. 13, pp. 1143–1146, 2003.