In the realm of power transmission, screw gears, particularly those involving worm drives, represent a critical class of mechanisms prized for their high reduction ratios, compact design, and self-locking potential. Among various innovative designs, the inclined double-roller enveloping hourglass worm drive presents a sophisticated evolution, aiming to address common drawbacks found in traditional designs, such as sensitivity to alignment errors, inability to adjust backlash, and weak load-bearing capacity at the worm root. The performance and longevity of these precision screw gears are intrinsically tied to the lubrication condition between the meshing tooth surfaces. In practical engineering, perfectly smooth surfaces are unattainable; machining processes like grinding inevitably leave a certain degree of surface roughness. When the magnitude of this roughness is comparable to, or even exceeds, the thickness of the lubricant film—which in elastohydrodynamic lubrication (EHL) regimes is often on the order of micrometers or sub-micrometers—its influence can no longer be neglected. Therefore, a comprehensive analysis that incorporates surface roughness effects is paramount for accurately predicting the lubrication performance and ensuring the reliable operation of these advanced screw gears.
The core of this analysis hinges on elastohydrodynamic lubrication theory, which describes the formation of a protective fluid film between two elastic, non-conforming surfaces in relative motion, where the fluid’s viscosity increases dramatically with pressure. For the inclined double-roller design, the contact between the worm thread and the roller surface can be mathematically simplified to a line contact problem. This involves modeling the interaction as that between an equivalent elastic cylinder and a rigid plane. The primary governing equation is the Reynolds equation, which describes pressure generation within the fluid film. For a transient, line-contact EHL problem considering roughness under isothermal conditions, the dimensionless form of this equation is central:
$$
\frac{d}{dX}\left( \varepsilon \frac{dP}{dX} \right) = \frac{d(\rho^* H)}{dX} + \frac{d(\rho^* H)}{dT}
$$
Here, $X = x/b$ is the dimensionless coordinate, $P = p/p_h$ is the dimensionless pressure, $H = h/R$ is the dimensionless film thickness, and $T$ is the dimensionless time. The parameter $\varepsilon$ is defined as $\varepsilon = \frac{\rho^* H^3}{\eta^* \lambda}$, where $\rho^*$ and $\eta^*$ are the dimensionless density and viscosity, respectively, and $\lambda$ is a dimensionless speed parameter. The right-hand side includes both the entrainment term ($d(\rho^* H)/dX$) and the squeeze film term ($d(\rho^* H)/dT$), which is crucial for capturing transient effects as the contact point moves along the tooth flank.
The film thickness equation must account for both the macroscopic geometry (Hertzian contact deformation) and the microscopic surface topography. Considering roughness on both the worm and roller surfaces, modeled here as transverse cosine waves, the dimensionless film thickness equation becomes:
$$
H(X) = H_0 + \frac{X^2}{2} – \frac{1}{\pi} \int_{X_{in}}^{X_{out}} \ln|X – X’| P(X’) dX’ – S(X, T)
$$
In this equation, $H_0$ is the dimensionless central rigid film thickness, the integral term represents the elastic deformation of the surfaces under pressure $P$, and $S(X, T)$ is the combined dimensionless roughness function for the two surfaces. If we denote the roughness amplitudes as $A_w$ and $A_r$, and the wavelengths as $l_w$ and $l_r$, for the worm and roller respectively, with surface velocities $v_w$ and $v_r$, the function can be expressed as:
$$
S(X, T) = A_w \cos\left[\frac{2\pi}{l_w}(X – U_w T)\right] + A_r \cos\left[\frac{2\pi}{l_r}(X – U_r T)\right]
$$
where $U_w$ and $U_r$ are the dimensionless surface speeds. The pressure-viscosity and pressure-density relationships are described by established empirical models. The Barus-like equation, often used in a more accurate exponential form, governs the viscosity change:
$$
\eta^* = \exp\left\{ (\ln \eta_0 + 9.67)\left[ (1 + 5.1 \times 10^{-9} p)^z – 1 \right] \right\}
$$
with $z = \alpha / [5.1 \times 10^{-9} (\ln \eta_0 + 9.67)]$, and $\alpha$ being the pressure-viscosity coefficient. The density variation is given by:
$$
\rho^* = \frac{1 + 0.6 \times 10^{-9}p}{1 + 1.7 \times 10^{-9}p}
$$
Finally, the solution must satisfy the force balance equation, ensuring the integrated pressure supports the applied load $W$:
$$
W = \int_{X_{in}}^{X_{out}} P(X) dX = \frac{\pi}{2}
$$
The geometrical and kinematic inputs for these equations, specific to the inclined double-roller enveloping hourglass screw gears, are derived from its meshing theory. Key parameters include the equivalent radius of curvature $R$, the entrainment velocity $v_{jx}$, and the load per unit length $w_i$ at the $i$-th meshing tooth pair. These parameters vary significantly as a tooth pair moves from the entry to the exit of the mesh zone.
| Parameter | Symbol | Value |
|---|---|---|
| Number of Worm Threads | $Z_1$ | 1 |
| Number of Worm Wheel Teeth | $Z_2$ | 25 |
| Center Distance | $A$ | 125 mm |
| Roller Radius | $R_k$ | 6.5 mm |
| Roller Offset | $c_2$ | 7 mm |
| Throat Diameter Coefficient | $k_1$ | 0.4 |
| Inclination Angle | $\gamma$ | 6° |
| Ambient Lubricant Viscosity | $\eta_0$ | 0.028 Pa·s |
| Ambient Lubricant Density | $\rho_0$ | 870 kg/m³ |
| Pressure-Viscosity Coefficient | $\alpha$ | 2.2 × 10⁻⁸ m²/N |
Solving the coupled system of equations (Reynolds, film thickness, force balance, and constitutive relations) is a computationally intensive task requiring advanced numerical techniques. The approach employed here utilizes the multi-grid method, which dramatically accelerates convergence by solving the problem on a hierarchy of grids with different discretization levels. The pressure is solved using the multi-grid solver, while the elastic deformation integral in the film thickness equation is efficiently calculated using the multi-grid integration method. The numerical procedure follows a rigorous cycle:
- Initialization: Start with an initial guess for pressure $P$ and film thickness $H$ on the coarsest grid.
- Relaxation: Use the Gauss-Seidel iterative method to solve the discretized Reynolds equation on the current grid level.
- Fine-to-Coarse (Restriction): Transfer the residual of the solution to a coarser grid.
- Coarse Grid Correction: Solve the error equation on the coarser grid.
- Coarse-to-Fine (Prolongation): Interpolate the correction back to the finer grid and update the solution.
- Cycle: Repeat steps 2-5 using a W-cycle pattern across 6 grid levels until convergence criteria for pressure ($E_p < 0.001$) and load balance ($E_w < 0.001$) are met on the finest grid (961 nodes).
- Transient Marching: The converged solution for one meshing instant serves as the initial guess for the next, allowing the simulation to track a tooth pair through the entire engagement.
The results reveal profound insights into the lubrication behavior of these screw gears. Comparing the “smooth” solution (ignoring roughness) with the “rough” solution (including roughness) at various meshing instants—single-tooth, double-tooth, triple-tooth, and quadruple-tooth contact phases—shows clear and significant differences. While the overall trends for both cases are similar (pressure peak increases then decreases, film thickness decreases then increases as the tooth moves through the mesh), the rough surface solution exhibits pronounced fluctuations superimposed on these trends. Each asperity passing through the high-pressure contact zone generates a local pressure spike and a corresponding dimple in the film thickness profile. Crucially, these fluctuations mean that the maximum pressure peak in the rough contact is substantially higher than in the smooth case, and the minimum film thickness is lower. This directly implies that surface roughness detrimentally affects the lubrication performance, increasing the risk of surface fatigue and wear, and potentially leading to a breakdown of the lubricant film in severe cases. The most critical phase occurs during triple-tooth engagement near the worm throat, where the load per unit length is highest and the entrainment velocity is lowest, resulting in the thinnest film.
Further analysis investigates the influence of key design parameters of the screw gears on the EHL performance under rough surface conditions. The parameters varied include the roller radius ($R_k$), the throat diameter coefficient ($k_1$), the roller offset distance ($c_2$), and the roller inclination angle ($\gamma$). The effects, evaluated at the critical triple-tooth meshing instant, are summarized qualitatively and quantitatively below.
| Design Parameter | Trend (Increase) | Effect on Pressure Fluctuation | Effect on Max Pressure | Effect on Min Film Thickness | Relative Influence |
|---|---|---|---|---|---|
| Roller Radius ($R_k$) | Increase | More Severe | Increases | Decreases | Medium |
| Throat Diameter Coeff. ($k_1$) | Increase | Less Severe | Decreases | Increases | Highest |
| Roller Offset ($c_2$) | Increase | More Severe | Increases | Decreases | Medium |
| Inclination Angle ($\gamma$) | Increase | Slightly More Severe | Slightly Increases | Slightly Decreases | Lowest |
The underlying mechanics can be interpreted as follows. Increasing the roller radius, offset, or inclination angle generally alters the contact geometry and kinematics in a way that reduces the equivalent radius of curvature or adversely affects the entrainment motion, making it harder to generate a thick hydrodynamic film. Consequently, the film becomes thinner, pressure rises, and the relative impact of surface roughness becomes more pronounced, leading to stronger fluctuations. A small throat diameter coefficient ($k_1$) means a smaller worm throat diameter relative to the center distance, which sharply increases the contact load density and severely impairs film formation, hence its dominant influence.
In conclusion, this detailed elastohydrodynamic lubrication analysis underscores critical aspects for the design and application of high-performance screw gears like the inclined double-roller enveloping hourglass worm drive. The explicit inclusion of surface roughness in the EHL model is not merely an academic refinement but a necessary step for realistic performance prediction. The results conclusively demonstrate that surface roughness induces significant fluctuations in pressure and film thickness, elevating peak pressures and reducing minimum film thickness, thereby posing a threat to the durability and efficiency of the gear set. For optimal lubrication performance, design parameters must be chosen judiciously: the roller radius, roller offset, and inclination angle should not be excessively large, and the throat diameter coefficient should not be too small. This work provides a robust numerical framework and practical guidance for enhancing the reliability of such sophisticated power transmission systems, ensuring that these advanced screw gears can operate under demanding conditions with improved longevity and reduced risk of failure.