A spiral bevel gear pair is often the first stage of a helicopter transmission because its curved teeth combine smooth direction-changing motion with high load capacity in a compact layout. That same curved tooth trace makes the lubricant film direction-dependent: the contact point moves across the face width, the entrainment velocity changes direction, and the contact ellipse rotates during mesh. Lubrication therefore does more than reduce friction. It changes the mesh stiffness, the mesh damping, and the dynamic transmission error that the gear pair must accommodate.
A curved-tooth spiral bevel gear pair. The contact point travels across the face width, so the entrainment velocity and contact-ellipse orientation vary during mesh.
Why the Lubricant Film Belongs in the Dynamic Model
Dry-contact gear dynamics usually treats the mesh as a time-varying spring and damper. That representation captures the periodic change from single-tooth to double-tooth contact, but it omits the compliance and viscous dissipation of the oil film. In a spiral bevel gear running at high speed, the film is not a passive surface layer. It adds a contact-scale spring in series with the tooth-body compliance, and it dissipates energy as the film is squeezed and sheared. Ignoring those effects can shift the predicted mesh force, transmission error, and vibration amplitude.
The coupling used here treats the film as part of the mesh constitutive behaviour. The dry contact stiffness is built from Hertz contact, tooth bending, shear, axial compression, and foundation compliance. In the lubricated case, the oil-film stiffness replaces the dry Hertz contact term at the interface, while the tooth-body terms remain. The mesh damping is formed from a structural component proportional to the instantaneous mesh stiffness and an oil-film component obtained from the energy dissipated in a loaded EHL contact. The dynamic transmission error is then modified by the central film thickness, so the film becomes a kinematic term as well as a stiffness and damping term.
Gear Geometry and Operating Point
The analysis uses a Gleason-type spiral bevel gear pair with a 23-tooth pinion and a 115-tooth wheel. The large ratio and the 90° shaft angle are representative of a helicopter input stage. The main geometry and material inputs are listed in Table 1.
| Parameter | Value |
|---|---|
| Teeth, pinion/wheel | 23/115 |
| Module | 6 mm |
| Pressure angle | 25° |
| Midpoint helix angle | 35° |
| Shaft angle | 90° |
| Face width | 60 mm |
| Pitch diameter, pinion/wheel | 138 mm/690 mm |
| Cone distance | 351.832 mm |
| Mass, pinion/wheel | 5.148 kg/165.65 kg |
| Elastic modulus | 207 GPa |
| Poisson ratio | 0.25 |
The contact is treated as a point contact between two ellipsoids. Because the tooth trace is curved, the instantaneous curvature radius along the mesh path is not constant. It is small near the start of mesh, rises toward the pitch region, and falls again near the end of mesh. The entrainment velocity, formed from the vector sum of the two surface velocities projected onto the tangent plane, changes continuously along the same path. These two kinematic features are the main reasons why a quasi-static, constant-curvature EHL model is insufficient for spiral bevel gears.
Time-Varying Mesh Stiffness and Damping With a Film
The dry mesh stiffness is assembled from five compliances: Hertz contact, bending, shear, axial compression, and gear-body foundation. The foundation term matters for thin-walled aerospace gears, where the web and rim deform under load. For a gear pair, the tooth-body compliances of the pinion and wheel act in series with the contact compliance. The resulting single-tooth stiffness rises as the contact moves from the tooth tip or root toward the fully engaged region, then falls as contact area is lost.
In the lubricated case, the oil film introduces an additional compliance at the interface. The film stiffness is derived from the EHL film-thickness solution, using the derivative of the central or minimum film thickness with respect to load. The gear-pair damping is scaled from the instantaneous mesh stiffness through a damping ratio of 0.07 and the equivalent mass of the pair. The film damping is obtained from the hysteresis loop of a harmonically loaded EHL contact: the area enclosed by the film-thickness versus film-force curve gives the energy dissipated per cycle, and that energy defines an equivalent viscous damping coefficient.
When the gear rotates, the periodic exchange between single-tooth and double-tooth contact produces a stepwise stiffness curve. Double-tooth contact raises the effective stiffness; single-tooth contact lowers it. In the calculated cases, adding the film raises both the effective mesh stiffness and the effective mesh damping relative to dry contact. The film also buffers the load transition between teeth, so the mesh-force fluctuation is smaller even though the effective stiffness is higher. This combination is the central dynamic effect: the film does not simply soften the mesh, it changes the amplitude and phase of the mesh response.
Point-Contact EHL Along a Curved Tooth Trace
The lubricated contact is solved as a transient point-contact EHL problem. The governing set contains the Reynolds equation, the film-thickness equation with elastic deformation, a viscosity–pressure relation, a density–pressure relation, and a load-balance equation. The equations are made dimensionless and discretized on a rectangular grid in the contact plane. The Reynolds equation is solved by finite-difference iteration with the boundary condition that pressure vanishes at the inlet and outlet and cannot fall below zero in the cavitated region.
The EHL solution gives the film thickness and pressure distribution at each mesh position. The film is asymmetric because the entrainment vector is not aligned with the contact-ellipse axes. At low speed, the main pressure peak is the global maximum. As speed or viscosity rises, the secondary pressure peak near the outlet can become the global maximum. That switch is a characteristic EHL feature and is visible in the calculated pressure fields.
Two operating variables were swept: input speed at fixed viscosity, and lubricant viscosity at fixed speed. The relative changes in minimum film thickness and maximum film pressure are summarized in Table 2. The percentages are relative to the preceding case in each sweep.
| Condition | Change in minimum film thickness | Change in maximum film pressure |
|---|---|---|
| 500 to 1000 r/min at 0.02135 Pa·s | +39.016% | +0.584% |
| 1000 to 1500 r/min at 0.02135 Pa·s | +27.980% | +5.229% |
| 0.00439 to 0.02135 Pa·s at 1500 r/min | +235.101% | +4.652% |
| 0.02135 to 0.03190 Pa·s at 1500 r/min | +22.046% | +6.231% |
The film thickness is much more sensitive to speed and viscosity than the peak pressure. The first viscosity increase produced the largest film-thickness change, which is consistent with the strong nonlinearity of the viscosity–pressure relation at low inlet viscosity. The pressure rise is modest by comparison, but it is still large enough to shift the location of the maximum pressure and to change the contact-load distribution.
Coupled Dynamic Response
The dynamic model uses eight degrees of freedom: three translations and one rotation for each gear. The equations include the time-varying mesh stiffness and damping, backlash, static transmission error, bearing stiffness and damping, and external torque. The mesh force is
$$F_m = k_t \lambda_n + c_t \dot{\lambda}_n$$
where \(k_t\) is the comprehensive time-varying mesh stiffness, \(c_t\) is the comprehensive time-varying mesh damping, and \(\lambda_n\) is the dynamic transmission error along the normal direction. The transmission error includes the geometric and static-error terms and the central film thickness, so the lubricant film directly changes the relative displacement that the mesh must accommodate.
The dynamic equations and the EHL equations are solved in an iterative loop. The Runge–Kutta method advances the lumped-parameter dynamic model; the finite-difference EHL solver updates the film thickness, film pressure, film stiffness, and film damping. The loop continues until the relative changes in the summed stiffness and damping fall below 0.001. The response calculation was run at 1500 r/min with a 100 N·m input torque.
Compared with dry contact at the same nominal operating point, the lubricated model predicts a larger effective mesh stiffness and damping, a smaller dynamic transmission error, and a smaller mesh-force swing. The peak-to-peak values are compared in Table 3.
| Response quantity | Dry contact | EHL model | Reduction |
|---|---|---|---|
| Dynamic transmission error | 1.85299 × 10−5 m | 1.51471 × 10−5 m | 18.256% |
| Single-tooth mesh force | 1928.0189 N | 1658.6190 N | 13.973% |
| Displacement | 5.12421 × 10−5 m | 4.11541 × 10−5 m | 19.687% |
| Velocity | 1.378908 × 10−3 m/s | 1.128098 × 10−3 m/s | 18.189% |
| Acceleration | 45.0901 m/s2 | 36.9326 m/s2 | 18.091% |
The reduction in transmission error and mesh force is not a simple consequence of a softer contact. The film adds compliance, but it also adds viscous damping and changes the phase of the mesh force. In the model, the net effect is a lower dynamic transmission error and a lower force amplitude. The displacement, velocity, and acceleration amplitudes follow the same trend because they are driven by the same mesh force and transmission-error histories.
Experimental Check
Validation used a spiral bevel gearbox test rig with a drive motor, the test gearbox, a planetary speed increaser, and an eddy-current brake. Two IEPE piezoelectric triaxial accelerometers were mounted on the housing near the pinion shaft and the wheel shaft. A DH5916 data acquisition system recorded the signals at a 10 kHz sampling rate. The lubricant was a 4106 synthetic aviation oil. Its viscosity was set by oil temperature: 100 °C gave 0.00439 Pa·s, 40 °C gave 0.02135 Pa·s, and 20 °C gave 0.03190 Pa·s. Tests were run at 500, 1000, 1500, and 2000 r/min with a 100 N·m output load. Dry-contact runs were kept short because of the risk of scuffing.
With increasing input speed, the measured X-direction acceleration at the pinion housing rose strongly. In the dry case, the RMS acceleration increased from 1.82099 m/s² at 500 r/min to 8.32360 m/s² at 2000 r/min, a factor of about 4.6. The EHL model reproduced the same upward trend. The largest difference between the model and the test RMS for the speed sweep was 19.305%, at 2000 r/min.
The viscosity sweep was less monotonic. At 500 and 2000 r/min, the highest viscosity gave the lowest RMS acceleration. Relative to dry contact, 0.03190 Pa·s reduced the RMS by 42.67% at 500 r/min and by 6.79% at 2000 r/min. Changing from 0.00439 to 0.03190 Pa·s reduced the RMS by 40.64% at 500 r/min and by 3.06% at 2000 r/min. At 1000 and 1500 r/min, however, the measured RMS at the highest viscosity was close to the dry value, within about 4%, and the trend was not clearly monotonic. The model predicted a monotonic decrease with viscosity, while the test did not. The largest model–test difference in this sweep was 28.229%, in the dry case.
| Input speed | Dry RMS acceleration | RMS at 0.03190 Pa·s | Reduction |
|---|---|---|---|
| 500 r/min | 1.82099 m/s² | 1.04400 m/s² | 42.67% |
| 1000 r/min | 2.84958 m/s² | 2.73984 m/s² | 3.85% |
| 1500 r/min | 5.51987 m/s² | 5.48040 m/s² | 0.72% |
| 2000 r/min | 8.32360 m/s² | 7.75883 m/s² | 6.79% |
The disagreement at intermediate speeds is physically plausible. The test changed viscosity by changing oil temperature, so viscosity and thermal effects were coupled. The EHL model is isothermal and assumes smooth surfaces. If the contact at 1000 or 1500 r/min did not develop a fully separating film, the measured vibration would be controlled by mixed or boundary lubrication rather than by the fully flooded EHL trend. Housing and bearing dynamics also contribute to the measured acceleration and are not represented in the eight-degree-of-freedom gear model.
Design Implications and Limits
The results point to a practical conclusion: for spiral bevel gears in high-speed, high-load service, the lubricant film should be treated as part of the dynamic design problem, not only as a wear-control measure. A film that raises effective mesh stiffness and damping can reduce dynamic transmission error and mesh-force swing. The benefit is operating-point dependent. Higher viscosity does not always produce lower housing vibration, especially when the contact is not fully separating or when speed-related excitation dominates.
The model has clear limits. The EHL solution is isothermal, Newtonian, and smooth-surface; it does not resolve roughness asperities or mixed lubrication. The dynamic model is a lumped-parameter gear-pair model and does not include the rotor, bearings, or housing as coupled flexible components. The experiments covered speeds up to 2000 r/min, and viscosity was changed through temperature rather than by using different oil formulations. Those limits matter for extrapolation to higher-speed helicopter transmissions.
Within those limits, the coupled EHL–dynamic model gives a consistent explanation of the measured trends. Speed is the dominant driver of vibration in the test range. Lubrication reduces vibration most clearly at the lower speed and at the higher speed where the film can develop, while the intermediate-speed behaviour shows that film formation and thermal effects cannot be ignored. The central film thickness should therefore be carried into the transmission-error calculation, and the film stiffness and damping should be updated with the operating point rather than held constant.