Worm gear and worm pair drives are widely used in industrial applications due to their high reduction ratio and compact design. However, the trial machining and testing of worm gear teeth are costly, making it highly valuable to accurately assess the feasibility of a design during the development phase. In this study, I present a comprehensive rigid-flexible coupling meshing analysis of a worm gear and worm pair using advanced simulation tools. The goal is to establish a reliable modeling and analysis methodology that can reduce development costs and time. I first construct a precise geometric model of the worm gear and worm pair in CATIA, then perform finite element analysis (FEA) in ANSYS to evaluate contact stresses. Next, I conduct rigid-body dynamics and rigid-flexible coupling dynamics analyses in ADAMS. The results from the rigid-flexible coupling analysis are compared with those from the finite element method, demonstrating good agreement and validating the proposed approach.
1. Introduction
Worm gear mechanisms are essential in many mechanical systems, yet their complex geometry and contact behavior make traditional design methods challenging. The ability to simulate the meshing process numerically before physical prototyping can significantly shorten the design cycle. In my work, I adopt a parametric modeling approach using CATIA, which allows for easy editing and re-meshing for subsequent analyses. The model is then exported to ANSYS for static contact analysis and to ADAMS for dynamic simulation. For the dynamic analysis, I employ a rigid-flexible coupling technique where the worm is treated as a flexible body to capture its vibrational effects, while the worm gear remains rigid. This approach provides a more realistic representation of the actual operating conditions.
2. Accurate Modeling of Worm Gear and Worm Pair
The worm gear and worm pair under study originate from an industrial application. The key parameters are summarized in Table 1. I used these parameters to construct the geometry in CATIA. The worm material is 42CrMo and the worm gear material is QA19-4, with properties listed in Table 2.
| Parameter | Symbol | Value |
|---|---|---|
| Module | m / mm | 1.25 |
| Number of worm threads | Z1 | 1 |
| Number of worm gear teeth | Z2 | 42 |
| Torque on worm gear | T2 / Nm | 13.923 |
| Center distance | a / mm | 36.50 |
| Lead angle of worm | γ / (°) | 3.49 |
| Lead of worm | P / mm | 5.65 |
| Pressure angle | α / (°) | 20.00 |
| Pitch diameter of worm | d1 / mm | 20.50 |
| Pitch diameter of worm gear | d2 / mm | 52.50 |
| Component | Material | Density / kg·m-3 | Poisson’s ratio | Young’s modulus / GPa |
|---|---|---|---|---|
| Worm | 42CrMo | 7 850 | 0.28 | 212 |
| Worm gear | QA19-4 | 7 500 | 0.33 | 116 |
Based on the mechanical design handbook, the allowable contact stress for the worm gear under rated motor power is calculated as:
$$ \sigma_H = Z_E \sqrt{ \frac{9\,400\,T_2}{d_1 d_2^2} K_A K_V K_\beta } = 286\ \text{MPa} $$
where \(Z_E = 157\ \sqrt{\text{MPa}}\) is the elastic coefficient, \(K_A = 1.5\) is the application factor, \(K_V = 1\) is the dynamic factor, and \(K_\beta = 1.2\) is the load distribution factor. The bending stress is:
$$ \sigma_F = \frac{666\,T_2 K_A K_V K_\beta}{d_1 d_2 m} Y_{FS} Y_\beta = 49\ \text{MPa} $$
with the helix angle factor \(Y_\beta = 1 – \gamma / 120 = 0.97\) and the compound tooth form factor \(Y_{FS} = 3.8\). These values serve as reference for subsequent numerical simulations.
3. Finite Element Analysis of Worm Gear and Worm Pair
To evaluate the contact stress distribution, I performed a three-dimensional finite element analysis using ANSYS. The geometric model from CATIA was imported and cleaned. The material properties from Table 2 were assigned. For boundary conditions, I fixed the worm gear by constraining its rotation about the Z-axis with a torque of 13.923 Nm acting opposite to the rotation direction, while the worm was allowed to rotate about the X-axis by an angle of 5° with all other degrees of freedom constrained. Tetrahedral solid elements were used for meshing, resulting in 243 741 elements. The contact between the worm gear and worm was defined as frictional with a friction coefficient of 0.06 for dynamic and 0.09 for static conditions.
The simulation revealed that the maximum contact stress on the worm gear tooth surface is 266 MPa. This stress occurs at the tooth root as a localized peak, while the contact band exhibits stresses in the range of 100–250 MPa. The stress distribution shows a strip-like pattern along the tooth flank, which is consistent with theoretical expectations for worm gear meshing.
4. Rigid-Body Dynamics Analysis in ADAMS
For the dynamic analysis, I modeled the worm gear pair as rigid bodies in ADAMS. The contact is defined using a nonlinear impact model based on the Hertz contact theory. According to Hertz, the contact force \(F\) is related to the indentation depth \(\delta\) by:
$$ \delta = \frac{a^2}{R} = \sqrt[3]{\frac{9 F^2}{16 R E^{*2}}} $$
$$ F = K \delta^{1.5} $$
where the contact stiffness \(K\) is given by:
$$ K = \frac{4}{3} \sqrt{R} E^{*} $$
with the equivalent radius \(R\) and equivalent Young’s modulus \(E^{*}\) defined as:
$$ \frac{1}{R} = \frac{1}{R_1} + \frac{1}{R_2} \quad \text{and} \quad \frac{1}{E^{*}} = \frac{1-\mu_1^2}{E_1} + \frac{1-\mu_2^2}{E_2} $$
Here \(R_1\) and \(R_2\) are the pitch circle radii of the worm and worm gear, respectively. Using the material data from Table 2, I computed \(K = 3.005 \times 10^{11}\ \text{N/mm}^{3/2}\). A damping coefficient \(C = 40\ \text{N·s/mm}\) was chosen after trial simulations. The dynamic friction coefficient was set to 0.06, and the static friction coefficient to 0.09. The force exponent was 1.5 and the maximum damping penetration depth was 0.1 mm. A torque of 13.923 Nm was applied to the worm gear using the step function:
$$ \text{step}( \text{time}, 0, 0, 0.5, 338.8 ) $$
to ensure smooth load application. The resulting meshing force over time is shown in the rigid-body simulation. The force rises rapidly to 565 N within 0.01 s, then gradually increases to a peak of 587.49 N at around 0.19 s, after which it stabilizes around 570 N with small oscillations due to the initial impact.
5. Rigid-Flexible Coupling Dynamics Analysis
To capture the structural flexibility of the worm, I performed a rigid-flexible coupling analysis. First, a modal analysis of the worm was conducted in HYPERMESH. I filtered out modes below 1 Hz and computed 16 modes, with the first natural frequency at 725 Hz, which I considered effective. A modal neutral file (MNF) was generated by defining the material, element properties, rigid elements, load collectors, constraints, and load steps in HYPERMESH, then solving with OptiStruct. The MNF file was imported into ADAMS using the Flex module. The worm was replaced by the flexible body, while the worm gear remained rigid. Revolute joints were added between the worm, worm gear, and ground. A contact pair of type “Flexible Body to Solid” was defined. The worm was driven by a velocity motion given by:
$$ \text{step}( \text{time}, 0, 0, 0.2, 3388\ \text{deg/s} ) $$
and a torque of 13.923 Nm was applied to the worm gear in the opposite direction. The simulation ran for 0.35 s with a step size of 20 000. The maximum contact stress on the worm gear tooth root was found to be 251 MPa. This value is only 5.6% lower than the 266 MPa obtained from the static finite element analysis, indicating good agreement. The meshing force from the flexible body simulation rises to 573 N in the first 0.001 s, peaks at 795.49 N at 0.138 s, and then fluctuates around 585 N. The larger fluctuations compared to the rigid-body case are due to the realistic inclusion of flexibility, friction, and inertial effects, which cause more dynamic variations.
6. Conclusion
In this work, I successfully developed a complete simulation chain for the worm gear and worm pair meshing analysis. The geometric model built in CATIA was used for both finite element static contact analysis and multibody dynamics simulation. The finite element analysis predicted a maximum contact stress of 266 MPa on the tooth root, while the rigid-flexible coupling analysis in ADAMS gave a stress of 251 MPa, with a deviation of only 5.6%. This close correlation validates the accuracy of the modeling method and the simulation approach. The rigid-flexible coupling technique provides a more realistic representation of the worm gear behavior under dynamic loading, including the effects of tooth flexibility and impact. The presented methodology can be directly applied to industrial product development, reducing the need for costly physical prototypes and accelerating the design process for worm gear transmissions.