Isothermal Elastohydrodynamic Lubrication Analysis of an Inclined Double-Roller Enveloping Hourglass Worm Gear Drive

In the field of power transmission, worm gear drives have long been valued for their high reduction ratios and compact design. However, traditional worm gear drives often suffer from significant sliding friction between the worm and worm wheel teeth, leading to reduced efficiency, increased wear, and a heightened risk of scuffing failure. To address these limitations, various forms of “live-rack” or “rolling-element” worm gear drives have been developed, where the teeth of the worm wheel are replaced by rolling elements such as balls, needles, or rollers. This fundamental change converts the predominant sliding friction into rolling friction, thereby promising substantial improvements in transmission efficiency and durability. Among these innovations, the inclined double-roller enveloping hourglass worm gear drive presents a particularly interesting configuration. In this worm gear drive, the worm wheel is split into two halves, each carrying circumferentially arranged rollers that can rotate about their own axes. The worm thread is generated by the envelope of these rollers, which are positioned at an inclination angle relative to the radial direction of the wheel. This design allows for adjustable backlash and maintains continuous contact. A critical aspect for the reliable operation of any worm gear drive, including this novel type, is its lubrication performance. The formation of a full elastohydrodynamic lubrication (EHL) film between the contacting surfaces is essential to prevent direct metal-to-metal contact, minimize wear, reduce operating temperatures, and enhance resistance to scuffing. Therefore, a detailed analysis of the EHL characteristics is paramount. This work focuses on establishing and solving an isothermal, steady-state, line-contact EHL model for the inclined double-roller enveloping hourglass worm gear drive. We aim to elucidate the pressure distribution and film thickness profile along the contact lines, analyze the lubrication performance at different meshing positions, and investigate the influence of operational parameters, specifically lubricant viscosity. The findings will provide a theoretical foundation for optimizing the design and lubrication of this efficient worm gear drive, ultimately contributing to its practical application and reliability.

The core of the EHL analysis lies in accurately modeling the contact conditions within the worm gear drive. The meshing between the worm and the rollers is a complex spatial conjugate action. At any given instant, multiple rollers are in contact with the worm thread, each forming a spatial contact line. To analyze the EHL performance, we simplify this complex line-contact problem based on Hertzian contact theory. The contact between the cylindrical roller (with radius $R_g$) and the worm thread surface (with local radius of curvature $R_w$ at the contact point) is equivalent to the contact between an equivalent elastic cylinder and a rigid plane. The equivalent radius $R$ and equivalent elastic modulus $E’$ are given by:

$$
R = \frac{R_w R_g}{R_w + R_g}
$$

$$
\frac{1}{E’} = \frac{1}{2} \left( \frac{1 – \mu_1^2}{E_1} + \frac{1 – \mu_2^2}{E_2} \right)
$$

where $\mu_1, E_1$ and $\mu_2, E_2$ are the Poisson’s ratio and elastic modulus of the worm and roller (worm wheel material) respectively. The worm thread radius $R_w$ varies along the contact path and can be determined from the meshing geometry. The simplified physical EHL model for our worm gear drive thus consists of a series of discrete equivalent cylindrical contacts distributed along each instantaneous spatial contact line. The load, surface velocities, and geometrical parameters change along this line. For the numerical solution, we consider a generic point $P$ on a contact line. The entrainment velocity $v_{jx}$, which is crucial for film formation, is the average of the surface velocities of the worm and roller projected onto the common normal direction at the contact point. For our specific worm gear drive geometry, these velocities are derived from its kinematic relations. If $v_{g1}$ and $v_{g2}$ are the normal velocities of the worm and roller surfaces at point P, then:

$$
v_{jx} = \frac{v_{g1} + v_{g2}}{2}
$$

The expressions for $v_{g1}$ and $v_{g2}$ involve parameters such as the rotational speeds, geometrical dimensions, the inclination angle $\gamma$, and the instantaneous position coordinates. The load per unit length $w_i$ on the $i$-th contacting roller is calculated considering the load distribution among simultaneously meshing teeth:

$$
w_i = \frac{F_{ni}}{L} = \frac{2 K_i T_1}{L d_1 \cos \alpha_n \cos \beta}
$$

where $F_{ni}$ is the normal force, $L$ is the contact line length, $T_1$ is the input torque on the worm, $d_1$ is the worm reference diameter, $\alpha_n$ is the normal pressure angle, $\beta$ is the lead angle, and $K_i$ is the load-sharing factor for the $i$-th contact pair. The factors $K_i$ are proportional to the contact stiffness, which relates to the instantaneous tooth thickness of the worm and wheel.

Having established the physical model, we now formulate the mathematical equations governing the isothermal steady-state line-contact EHL problem for this worm gear drive. The governing equations are the Reynolds equation, the film thickness equation, the viscosity-pressure relation, the density-pressure relation, and the force balance equation. We employ Hertzian contact parameters to non-dimensionalize these equations for numerical stability and generality. Let $h$ be film thickness, $p$ pressure, $x$ the coordinate along the rolling direction, $b$ the Hertzian contact half-width, $p_h$ the maximum Hertzian pressure, $\eta_0$ the ambient viscosity, $\rho_0$ the ambient density, and $v_{jx}$ the entrainment speed. The non-dimensional variables are defined as: $H = h/R$, $P = p/p_h$, $X = x/b$, $\eta^* = \eta/\eta_0$, $\rho^* = \rho/\rho_0$, $W = w/(E’R)$, $U = \eta_0 v_{jx} / (E’R)$. Here, $w$ is the load per unit length. The non-dimensional equations are as follows:

1. Reynolds Equation:

$$
\frac{d}{dX} \left( \varepsilon \frac{dP}{dX} \right) = \frac{d (\rho^* H)}{dX}
$$

where $\varepsilon = \frac{\rho^* H^3}{\eta^* \lambda}$ and $\lambda = \frac{12 \eta_0 v_{jx} R^2}{p_h b^2}$. The boundary conditions are $P(X_{in}) = 0$ at the inlet and $P(X_{out}) = dP(X_{out})/dX = 0$ at the outlet.

2. Film Thickness Equation:

$$
H(X) = H_0 + \frac{X^2}{2} – \frac{1}{\pi} \int_{X_{in}}^{X_{out}} \ln |X – S| \, P(S) \, dS
$$

where $H_0$ is the non-dimensional central film thickness.

3. Viscosity-Pressure Relation (Roelands equation):

$$
\eta^* = \exp \left\{ (\ln \eta_0 + 9.67) \left[ -1 + (1 + p_h P / p_0)^z \right] \right\}
$$

where $p_0$ is a reference pressure and $z$ is the viscosity-pressure index.

4. Density-Pressure Relation (Dowson-Higginson equation):

$$
\rho^* = 1 + \frac{0.6 \times 10^{-9} p_h P}{1 + 1.7 \times 10^{-9} p_h P}
$$

5. Load Balance Equation:

$$
\int_{X_{in}}^{X_{out}} P(X) \, dX = \frac{\pi}{2}
$$

These five equations constitute the complete isothermal line-contact EHL mathematical model for our analysis of the worm gear drive.

To obtain the pressure distribution $P(X)$ and film thickness profile $H(X)$, we solve the coupled integro-differential system numerically. We employ the Newton-Raphson method combined with a relaxation technique (the “descent method”). The Reynolds equation is discretized using finite differences. The film thickness equation is discretized using a standard influence coefficient technique. The resulting set of non-linear algebraic equations for the nodal pressures $P_i$, the central film thickness $H_0$, and the outlet boundary condition (often handled via the $\rho^*_{out} H_{out}$ term) is solved iteratively. The discretized forms are:

Discretized Reynolds Equation:

$$
f_i = \left( \frac{H_i^3}{\eta^*_i} \frac{dP}{dX} \bigg|_i \right) – A \left( \bar{\rho}_i H_i – \frac{\rho^*_{out} H_{out}}{\bar{\rho}_i} \right) = 0
$$

where $A = \frac{3\pi^2 U}{4 W^2}$, and overbars denote averaged values between nodes.

Discretized Film Thickness Equation:

$$
H_i = H_0 + \frac{X_i^2}{2} + \sum_{j=1}^{n} K_{ij} P_j
$$

where $K_{ij}$ are the influence coefficients derived from the elastic deformation integral.

Discretized Load Balance Equation:

$$
\sum_{i=1}^{n} P_i \Delta X = \frac{\pi}{2}
$$

In the iteration process, a damping factor $\lambda$ is introduced to ensure convergence: $P^{k+1} = P^k + \lambda \Delta P$, $H_0^{k+1} = H_0^k + \lambda \Delta H_0$, etc. This numerical framework allows us to analyze the EHL performance at any point along the contact lines in the worm gear drive.

We now present the results of our isothermal EHL analysis for a specific design of the inclined double-roller enveloping hourglass worm gear drive. The key parameters are: worm threads $z_1=1$, wheel teeth (rollers) $z_2=25$, center distance $A=125$ mm, throat diameter coefficient $k_1=0.4$, roller radius $R_g=6.5$ mm, roller offset $c_2=7$ mm, inclination angle $\gamma=6^\circ$. The lubricant has an ambient viscosity $\eta_0 = 0.028$ Pa·s. The material properties are: $E_1 = E_2 = 210$ GPa, $\mu_1 = \mu_2 = 0.3$. The Roelands parameters are $z=0.68$, $p_0$ is taken standard. The analysis focuses on the contact conditions at different positions: along a specific contact line (e.g., Contact Line 2), at the dedendum circle, reference circle, and addendum circle of the worm wheel, across different contact lines, and under varying lubricant viscosity.

General EHL Characteristics: For a typical point on Contact Line 2, the calculated non-dimensional pressure $P$ and film thickness $H$ profiles exhibit classical line-contact EHL features. The central region of the contact shows a nearly parallel film, with pressure closely following the Hertzian ellipse. In the exit region, a sharp secondary pressure spike is observed, immediately followed by a rapid drop to ambient pressure. Corresponding to this spike, the film profile shows a pronounced constriction or “necking,” and the minimum film thickness $H_{min}$ occurs at this constriction. This behavior is fundamental to the elastohydrodynamic lubrication in heavily loaded contacts like those in this worm gear drive.

Lubrication Characteristics Along the Path of Contact: We analyze the EHL performance at three characteristic circles on the worm wheel: the dedendum circle (root), the reference circle (pitch), and the addendum circle (tip). The results are summarized in the table below, showing key parameters like the secondary pressure spike magnitude $P_{spike}$, its location $X_{spike}$, and the minimum film thickness $H_{min}$.

Position Secondary Pressure Spike $P_{spike}$ Spike Location $X_{spike}$ Minimum Film Thickness $H_{min}$ Necking Onset
Dedendum Circle High Closer to exit Smallest Most delayed
Reference Circle Medium Intermediate Medium Intermediate
Addendum Circle High Closer to inlet Largest Earliest

The pressure spike magnitudes are relatively high at both root and tip but occur at different locations. The film thickness increases significantly from the dedendum to the addendum. The variation in film thickness is smaller between the reference circle and addendum circle compared to the variation between the dedendum and reference circles. This indicates that the contact region from the pitch circle to the tip experiences superior lubrication performance in this worm gear drive, with thicker films that better separate the surfaces.

Lubrication Characteristics Across Different Contact Lines: In a multi-tooth engagement, typically four contact lines exist simultaneously. We label them Contact Lines 1, 2, 3, and 4 in sequence. Their EHL characteristics differ due to variations in load sharing, entrainment velocity, and contact geometry. The following table compares these characteristics.

Contact Line Secondary Pressure Spike $P_{spike}$ Spike Location $X_{spike}$ Minimum Film Thickness $H_{min}$ Necking Behavior
Line 1 Medium-High Intermediate Medium Intermediate
Line 2 Medium Intermediate Medium-Low Intermediate
Line 3 Highest Closest to exit Smallest Most delayed
Line 4 Lowest Closest to inlet Largest Earliest

Contact Line 3 experiences the most severe conditions with the highest pressure spike and thinnest film. Contact Line 4, conversely, shows the most favorable conditions with the lowest pressure spike and thickest film. The necking phenomenon appears earliest for Contact Line 4 and latest for Contact Line 3. This variation must be considered in the design to ensure all meshing pairs in the worm gear drive have adequate lubrication.

Influence of Lubricant Viscosity: Lubricant properties significantly impact EHL performance. We analyze the effect of increasing the ambient dynamic viscosity $\eta_0$ while keeping other worm gear drive parameters constant. The results are summarized below.

Viscosity $\eta_0$ (Pa·s) Secondary Pressure Spike $P_{spike}$ Spike Location Shift Minimum Film Thickness $H_{min}$ Necking Onset
0.020 (Lower) Lower Towards exit Smaller Delayed
0.028 (Reference) Reference Reference Reference Reference
0.040 (Higher) Higher Towards inlet Larger Earlier

As viscosity increases, both the secondary pressure spike magnitude and the minimum film thickness increase. The location of the pressure spike moves towards the inlet region of the contact. Consequently, the film constriction or necking occurs earlier in the exit region. The enhancement of film thickness with higher viscosity is particularly pronounced. This underscores the importance of selecting a sufficiently high-viscosity lubricant for this worm gear drive to ensure robust elastohydrodynamic film formation and protection against wear and scuffing.

The comprehensive isothermal elastohydrodynamic lubrication analysis of the inclined double-roller enveloping hourglass worm gear drive leads to several important conclusions. First, the established physical and mathematical models, based on Hertzian contact theory and the line-contact EHL framework, are effective for analyzing the lubrication performance of this complex worm gear drive. The numerical solution using the Newton-Raphson method provides detailed pressure and film thickness profiles. Second, the lubrication characteristics are not uniform across the meshing zone. The region from the reference (pitch) circle to the addendum (tip) circle of the worm wheel generally exhibits superior EHL performance with thicker oil films compared to the dedendum (root) region. This insight is valuable for design and failure analysis. Third, among the simultaneously engaged contact lines, significant variation exists. Contact Line 4 consistently shows the most favorable conditions with the largest film thickness and smallest secondary pressure spike, while Contact Line 3 experiences the most severe contact. This non-uniformity in the worm gear drive’s lubrication must be accounted for in load distribution calculations and life predictions. Fourth, lubricant viscosity is a critical operational parameter. Higher viscosity directly leads to the formation of thicker elastohydrodynamic films, thereby improving the lubrication performance and likely enhancing the scuffing resistance and efficiency of the worm gear drive. These findings form a theoretical foundation for further investigations into thermal effects, transient conditions, and friction losses in this innovative worm gear drive. Ultimately, optimizing the EHL performance is key to realizing the full potential of this high-efficiency, low-backlash worm gear drive in demanding power transmission applications.

Scroll to Top