The pursuit of power transmission systems that combine high efficiency, compactness, and high load capacity has been a constant endeavor in mechanical engineering. Among the various configurations, the toroidal planetary worm gear drive stands out as a particularly promising solution. This innovative mechanism ingeniously integrates rolling contact with planetary worm drive principles. It offers notable advantages such as smooth operation, a wide range of transmission ratios, high load-bearing capability, and a remarkably compact structure, positioning it as a superior form of mechanical transmission. However, the very features that grant its advantages—the complex spatial assembly and the intricate geometry of its key components, such as the toroidal stator and the enveloping worm—also pose significant challenges in manufacturing and precise assembly. Practical experience with physical prototypes has revealed unexpected levels of vibration and noise during operation, adversely affecting transmission smoothness and reliability. These undesirable dynamic characteristics underscore the critical necessity for in-depth research into the dynamics of the toroidal planetary worm gear drive.
Over the past two decades, foundational research has been conducted on the dynamics of this system, covering areas like model establishment, free vibration, modal analysis, and dynamic response. A common limitation in many previous studies is the treatment of the system as purely rigid-body and the simplification of the rolling elements (spherical rollers) as integral parts of the planet gears. This simplification overlooks the fundamental characteristic of the drive: the replacement of sliding meshing with rolling contact. While some researchers have begun to consider the independent motion of rollers, the analyses often remain at the level of complex theoretical derivations. Recently, the integration of multi-body dynamics simulation software like Adams with finite element analysis tools like ANSYS has provided a powerful avenue for investigating complex gear dynamics, including rigid-flexible coupled behaviors. This study leverages this virtual prototyping approach to construct and analyze both a rigid-body and a rigid-flexible coupled dynamics model of the toroidal planetary worm gear drive, with a particular focus on the influence of internal excitations.
Dynamic Modeling Methodology
Governing Dynamic Equations
The dynamic analysis is performed within a multi-body dynamics framework. The system’s equations of motion, considering internal excitations, can be generally expressed as:
$$ M\ddot{q} + C\dot{q} + Kq = F_{ext} + F_{int} $$
where \( M \), \( C \), and \( K \) are the global mass, damping, and stiffness matrices, respectively. \( q \), \( \dot{q} \), and \( \ddot{q} \) are the displacement, velocity, and acceleration vectors. \( F_{ext} \) represents external forcing terms (e.g., input torque, load torque), and \( F_{int} \) encapsulates the internal excitation forces arising from the meshing process.
The internal excitation force \( F_{int} \) is primarily composed of two key components: the time-varying mesh stiffness excitation and the transmission error excitation. It can be formulated as:
$$ F_{int} = -\Delta K(t) e(t) + S(t) $$
where \( \Delta K(t) \) is the time-varying component of the mesh stiffness, \( e(t) \) is the kinematic transmission error vector, and \( S(t) \) is an impact vector resulting from the interplay between stiffness variation and error, often expressed as \( S(t) = -\Delta K(t)(q – q_s) – K e(t) \), with \( q_s \) being the static displacement.
Virtual Prototype Development
A three-dimensional model of the toroidal planetary worm gear drive is developed. The primary design parameters for the key components are summarized in Table 1.
| Component | Number of Teeth/Threads | Pitch Diameter at Symmetry Section (mm) |
|---|---|---|
| Worm (Input) | 1 (Single-Start Thread) | 113 |
| Planet Gear | 8 | 125 |
| Toroidal Stator (Fixed Ring) | 23 | 363 |
| Spherical Roller | 16 (per planet, 4 planets total) | – |
This model is imported into Adams/View. All material properties, constraints, and contacts must be redefined within the simulation environment. The following kinematic joints are applied to reflect the actual working principle of the worm gear drive:
- A revolute joint between the worm shaft and the ground (housing).
- A revolute joint between the planet carrier (output) and the ground.
- Revolute joints between each planet gear and the planet carrier.
- A fixed joint between the toroidal stator and the ground.
- Spherical joints between each spherical roller and its corresponding planet gear (to allow for the rollers’ independent rotational degrees of freedom).
Contacts are defined between the worm and the rollers, and between the rollers and the toroidal stator’s internal helical raceway. A penalty-based contact method with an impact function is typically employed. A constant rotational speed is applied to the worm input, ramped up over 0.1 seconds to avoid step-function shocks. A constant resistive torque is applied to the planet carrier output, similarly ramped up.
Bearing Support Modeling
Bearings are critical flexible elements in any gear system. Rather than simplifying them as ideal joints or single-direction springs, a more realistic representation is used. The support provided by the bearings at the worm shaft and planet carrier is modeled using 6×6 stiffness and damping matrices, which relate forces and moments to displacements and rotations. For simplicity, and because coupling between directions is often secondary, the matrices are treated as diagonal:
$$ K_{bearing} = diag(K_x, K_y, K_z, K_{\theta_x}, K_{\theta_y}, K_{\theta_z}) $$
$$ C_{bearing} = diag(C_x, C_y, C_z, C_{\theta_x}, C_{\theta_y}, C_{\theta_z}) $$
The force-displacement relationship at a bearing location is then:
$$ F_{body} = K_{bearing} \cdot q_{body} + C_{bearing} \cdot \dot{q}_{body} $$
where \( q_{body} = [x, y, z, \theta_x, \theta_y, \theta_z]^T \) is the displacement vector of the body at the bearing connection point relative to a reference frame.
Flexibilization of the Toroidal Stator
A significant advancement over pure rigid-body models is the creation of a rigid-flexible coupled model. The toroidal stator, often a thin-walled component, is a prime candidate for flexibilization due to its potential for elastic deformation under load. The process involves:
- Exporting the stator geometry to ANSYS or a similar finite element (FE) software.
- Defining material properties (e.g., structural steel), meshing the component with suitable solid or shell elements, and defining connection points (nodes) where it interfaces with other parts (e.g., mounting points to ground).
- Exporting the meshed model as an MNF (Modal Neutral File), which contains information about its mass, stiffness, and mode shapes.
- Importing the MNF file into Adams and replacing the rigid stator body with this flexible component. The software automatically transfers previously defined constraints and forces to the appropriate nodes on the flexible body.
This integrated model now captures the dynamic interaction between the rigid components (worm, planet gears, rollers, carrier) and the flexible stator, allowing for a more accurate simulation of system vibration and dynamic mesh forces.
Analysis of Internal Excitations
The dynamic response of the worm gear drive is largely driven by internal excitations generated during operation. The two primary sources are quantitatively analyzed below.
Time-Varying Mesh Stiffness Excitation
The mesh stiffness in a toroidal planetary worm gear drive is not constant. It fluctuates because the number of rolling elements in simultaneous contact changes as the planet gears rotate. The instantaneous mesh stiffness \( K_m(t) \) for a given meshing pair (worm-roller or stator-roller) can be approximated as:
$$ K_m(t) = n(t) \cdot K_r $$
where \( K_r \) is the contact stiffness of a single roller-raceway pair (which itself can be calculated using Hertzian contact theory), and \( n(t) \) is the number of rollers sharing the load at time \( t \).
The number of contacting rollers is governed by the contact ratio. For a worm gear drive with a multi-threaded enveloping worm and a multi-toothed stator, the contact ratios for the worm-planet and stator-planet meshes are given by:
$$ \varepsilon_w = \frac{\psi_w Z_p}{360} $$
$$ \varepsilon_s = \frac{\psi_s Z_p}{360} $$
where \( \varepsilon_w \) and \( \varepsilon_s \) are the contact ratios for the worm and stator meshes, respectively. \( Z_p \) is the number of rollers on a planet gear (analogous to teeth). \( \psi_w \) and \( \psi_s \) are the angular spans of the worm thread and stator tooth that are in contact with the planet’s rollers, respectively. \( \psi_p = 360/Z_p \) is the angular pitch between rollers.
For the example system in Table 1, calculations yield \( \varepsilon_w \approx 2.2 \) and \( \varepsilon_s \approx 2.4 \). This implies that the number of rollers in contact with the worm fluctuates between 2 and 3, and with the stator between 2 and 3 as well. The periodic change in \( n(t) \) from 2 to 3 and back again directly causes a time-varying mesh stiffness \( \Delta K(t) \), which acts as a parametric excitation on the system. The stiffness variation for one planet gear’s meshes over time is conceptually shown in Figure 1 (Note: Actual simulation data would be used to plot this).
Transmission Error Excitation
Transmission error (TE) is defined as the difference between the actual position of the output and the position it would occupy if the drive were perfectly rigid and geometrically ideal. It is a primary source of vibration and noise in gear systems. For the toroidal planetary worm gear drive, TE arises from several factors:
- Manufacturing Errors: Deviations in the lead of the worm, pitch errors on the stator, and sphericity errors of the rollers.
- Assembly Errors: Eccentricity of the worm shaft or planet carrier, misalignments.
- Elastic Deformations: Deflections of the worm, stator, and supports under load.
The composite static transmission error \( e(t) \) along the line of action for each mesh can be modeled as a displacement excitation. It often contains harmonic components at the mesh frequency and its multiples. A simplified representation for a single mesh might be:
$$ e(t) = E_0 + \sum_{m=1}^{N} E_m \sin(m \omega_m t + \phi_m) $$
where \( E_0 \) is a constant error, \( E_m \) is the amplitude of the m-th harmonic, \( \omega_m \) is the mesh frequency (\( \omega_m = Z_s \cdot \omega_{carrier} \), where \( Z_s \) is stator tooth number and \( \omega_{carrier} \) is carrier speed), and \( \phi_m \) is the phase angle.
In the rigid-flexible coupled model, the elastic deformation of the stator, calculated by the FE component, dynamically contributes to \( e(t) \), making it more severe than in the purely rigid model, especially when operational clearances are present.
Dynamic Characteristics and Response Comparison
Simulations are conducted for both the rigid-body model (RBM) and the rigid-flexible coupled model (RFCM) under identical input conditions (worm speed = 157.08 rad/s ~ 1500 RPM, output torque applied to carrier). The dynamic responses are analyzed in the time and frequency domains.
Kinematic Response
The angular velocities of key components are extracted. The worm input speed is constant. The output carrier speed and the planet gear speeds oscillate around their theoretical mean values due to internal excitations.
$$
\omega_{carrier, theory} = \frac{\omega_{worm}}{(Z_s / Z_p) + 1} \approx 6.54 \text{ rad/s}
$$
$$
\omega_{planet, theory} = \omega_{carrier} \cdot (Z_s / Z_p) \approx 18.8 \text{ rad/s}
$$
Simulation results show the mean values for the RBM are close to these theoretical values. However, for the RFCM, the mean carrier and planet speeds are slightly higher (e.g., 6.72 rad/s and 21.3 rad/s). More importantly, the amplitude of oscillation is significantly larger in the RFCM. This is attributed to the elastic deformation of the stator, which effectively increases the system’s “backlash” or compliance, allowing for momentary loss of contact (jumping) and sliding between rollers and raceways. These nonlinear events manifest as large spikes in the angular velocity plots of the planets in the RFCM, which are absent in the smoother but still oscillatory RBM response.
Roller Behavior and Contact Forces
Modeling the rollers as independent bodies reveals their complex motion. A roller’s spin velocity is not constant. It experiences sudden acceleration when it engages with the high-speed worm, maintains a relatively high speed during non-contact phases (assuming low friction in the spherical joint), and then decelerates when it engages with the slow-moving stator raceway. This cyclic, discontinuous acceleration of all 32 rollers in the system acts as a significant internal disturbance, contributing to overall vibration.
The contact forces between the rollers and the raceways are highly discontinuous in both models, appearing as a series of pulses. This is a direct consequence of manufacturing/assembly errors and clearances modeled in the system, which prevent perfect continuous contact. The contact forces in the RFCM are generally larger in magnitude and have longer contact durations compared to the RBM. The flexible stator deforms to maintain contact over a slightly longer arc, but the forces required to cause this deformation are higher. The impulsive nature of these contact forces is a key driver of vibration and potential fatigue failure in the worm gear drive components.
Vibration Analysis in Time and Frequency Domains
The angular acceleration of the planet carrier is a direct indicator of output vibration and torsional oscillation. Time-domain plots show that the acceleration fluctuates with large amplitude around zero. The peak-to-peak amplitude of these fluctuations is an order of magnitude larger for the RFCM than for the RBM, confirming that the inclusion of structural flexibility drastically increases the perceived vibration level.
Frequency domain analysis (FFT of the carrier acceleration) provides further insight. The fundamental mesh frequency \( f_m \) for this drive can be calculated as:
$$ f_m = \frac{\omega_{carrier} \cdot Z_s}{2\pi} $$
With \( \omega_{carrier} \approx 6.7 \) rad/s and \( Z_s = 23 \), \( f_m \approx 24.5 \) Hz. The frequency spectra for both models show dominant peaks at the mesh frequency and its harmonics (\( 2f_m, 3f_m, \) etc.). However, the model in which the largest resonance peak occurs is different:
- RBM: The highest resonance peak typically occurs at a lower harmonic (e.g., \( 2f_m \)).
- RFCM: The highest resonance peak shifts to a higher harmonic (e.g., \( 6f_m \)). This indicates that the natural frequencies of the system increase when the stator’s local flexibility is considered, as the flexible body has its own set of higher-frequency modal characteristics.
Furthermore, the spectral peaks in both models are surrounded by sidebands. These sidebands are created by the modulation of the mesh frequency vibrations by other periodic phenomena, such as the rotation frequency of the planet carrier (\( f_c \)) and the passing frequency of the planets (\( N_p \cdot f_c \), where \( N_p \) is the number of planets). The presence of these sidebands is a classic signature of internal excitation in planetary gear systems and confirms that the time-varying stiffness and error are actively modulating the dynamic response of the worm gear drive.
Conclusion
This investigation into the rigid-flexible coupled dynamics of the toroidal planetary worm gear drive yields several critical insights for its design and application:
- Severity of Vibration: Both rigid-body and coupled models exhibit significant vibration primarily driven by internal excitations—time-varying mesh stiffness and transmission error. The vibration level in the rigid-flexible coupled model is substantially higher than in the rigid-body model. This increase is attributed to the elastic deformation of the stator, which amplifies effective errors and can lead to nonlinear events like tooth jumping and sliding, severely impacting transmission smoothness.
- Importance of Roller Dynamics: Modeling the spherical rollers as independent bodies is crucial for an accurate dynamic representation. Their motion cycle, characterized by sudden acceleration/deceleration and high-speed free spin during non-contact phases, constitutes a major source of internal disturbance that exacerbates system-wide vibration.
- Nature of Contact Forces: Due to inevitable errors and clearances, the contact forces in the worm gear drive are highly impulsive and discontinuous. This impulsive loading is a key concern for durability, noise generation, and component fatigue life, especially in high-power applications.
- Spectral Characteristics: The rigid-flexible coupled model exhibits resonance at higher frequencies compared to the rigid-body model, reflecting the influence of the flexible component’s modal properties. The frequency spectra for both models are rich in content, with mesh harmonics and modulation sidebands clearly visible, confirming that internal excitations profoundly shape the dynamic signature of the drive.
The integrated approach using multi-body dynamics and finite element analysis proves to be a powerful and necessary tool for understanding the complex behavior of advanced worm gear drive systems. It moves beyond simplified theoretical models to capture critical nonlinear interactions and flexibility effects. Future work should focus on investigating the influence of external load fluctuations, optimizing the flexible component’s design (e.g., stator ribbing) to mitigate vibration, and conducting experimental validation to correlate and refine these simulation models. This research provides a foundational reference for the dynamic optimization and reliable design of high-performance toroidal planetary worm gear drives.