Spiral bevel gears carry the input stage of most helicopter main transmissions because a single pair can turn the drive through 90° while keeping the gearbox short and the contact ratio high. Gleason-type teeth, with their curved longitudinal form and tapered tooth depth, combine a large face contact area with good resistance to pitting, which is why they dominate this position despite the cost of generating them. The operating condition that receives least attention in conventional dynamic models is also the one that governs surface survival: the lubricating film between the teeth.
In a spiral bevel pair the contact is a point contact that sweeps across the tooth surface. The contact ellipse is inclined at a changing angle to the sliding-velocity vector, the entrainment velocity rotates along the path of contact, and the composite radius of curvature passes through a maximum near the pitch line. The film therefore has a strongly asymmetric, time-dependent thickness distribution rather than the near-uniform film assumed in quasi-static lubrication checks. Because that film sits directly between the two elastic bodies, it modifies the local compliance of the mesh and dissipates energy through viscous shear. Both effects appear in the same coefficients that a gear dynamics model uses — mesh stiffness and mesh damping — and both feed back into the dynamic transmission error, which in turn changes the load that generates the film. Treating lubrication and dynamics as separate problems breaks this loop.
Geometry of the Analysis Case
The case examined here is a Gleason-generated pair with a 23-tooth left-hand pinion driving a 115-tooth right-hand wheel on a 90° shaft angle. The design follows the usual constraints for this tooth form: a mid-face helix angle between 25° and 40° so that the face contact ratio stays in the 1.25–1.75 band, a pressure angle of 25° to avoid undercut without sacrificing transverse contact ratio, and identical clearance at both ends of the tooth, which forces the pinion addendum angle to equal the wheel dedendum angle. Tooth depth tapers from the outer to the inner end, so every radial station has a different section modulus.
| Parameter | Pinion | Wheel |
|---|---|---|
| Number of teeth, Z | 23 | 115 |
| Module, m | 6 mm | 6 mm |
| Pressure angle, α | 25° | 25° |
| Mid-face helix angle, βm | 35° | 35° |
| Shaft angle, Σ | 90° | 90° |
| Face width, B | 60 mm | 60 mm |
| Hand | Left | Right |
| Pitch diameter, d | 138.00 mm | 690.00 mm |
| Base circle diameter, db | 129.677 mm | 648.388 mm |
| Tip diameter, da | 152.407 mm | 691.119 mm |
| Root diameter, df | 129.458 mm | 686.821 mm |
| Cone distance, Re | 351.832 mm | 351.832 mm |
| Pitch cone angle, δ | 12.683° | 77.317° |
| Face cone angle, γ | 13.392° | 78.636° |
| Root cone angle, δf | 10.670° | 79.330° |
| Addendum, ha | 6.972 mm | 3.228 mm |
| Mounting distance | 358.00 mm | 175.00 mm |
Time-Varying Mesh Stiffness of the Dry Contact
For a pair of this size the mesh stiffness is built from five compliances in series: Hertzian contact, tooth bending, tooth shear, axial compression, and the fillet-foundation rotation of the thin web that carries the tooth. The first is closed-form under the Hertz assumption; for smooth elastic contacts of this geometry the elliptical pressure distribution is normally taken as accurate to within about half a percent of the true distribution. The other four are obtained from strain energy. Each tooth is sliced into thin sections normal to the pitch cone, and each slice is treated as a varying-section cantilever loaded at the contact point, so that
$$U_b=\frac{F^2}{2k_b},\qquad U_s=\frac{F^2}{2k_s},\qquad U_a=\frac{F^2}{2k_a}$$
with the bending, shear and axial energies integrated along the tooth from the root circle to the load line. Because the sections are conical, the moment arm, the section height and the second moment of area all vary with the angular position αk of the contact point, and the resulting stiffnesses are functions of that angle rather than constants. The foundation term cannot be neglected for aerospace gears, whose web is deliberately thin: the deflection of the web under load is of the same order as the tooth deflection itself. It is written as a closed-form expression in the tooth-root geometry and the load angle, with the four dimensionless coefficients L*, M*, P*, Q* fitted as polynomials in the web slenderness ratio and the half-tooth angle.
Coupling the five compliances for both members gives the dry mesh stiffness
$$\frac{1}{k_g}=\sum_{i=z,c}\left(\frac{1}{k_{h,i}}+\frac{1}{k_{b,i}}+\frac{1}{k_{s,i}}+\frac{1}{k_{a,i}}+\frac{1}{k_{f,i}}\right)$$
where the sum runs over pinion (z) and wheel (c). As the contact travels from tip to root the single-tooth stiffness rises to a maximum near the pitch line and falls away at both ends, following an approximately parabolic profile. Because the contact ratio exceeds one, single- and double-tooth intervals alternate, and the total stiffness becomes a rectangular-wave-like staircase in which double-tooth intervals show the sum of two single-tooth contributions. This staircase is the dominant internal excitation of the system.
Adding the Film: Stiffness and Damping of the Lubricated Mesh
Once a film separates the surfaces, the near-surface compliance is no longer that of the Hertzian contact but that of the film itself. The series chain is therefore modified by replacing kh with a film stiffness ko and summing over the n tooth pairs instantaneously in mesh. The mesh damping follows a proportional rule: the dry pair damping is taken as proportional to the instantaneous stiffness through the damping ratio ζ and the equivalent mass of the pair,
$$c_g=2\zeta k_g\sqrt{\frac{m_e}{k_{\mathrm{avg}}}},\qquad m_e=\frac{m_z m_c}{m_z+m_c}$$
with ζ = 0.07 for the geared system and kavg the stiffness averaged over a complete mesh cycle. Under lubrication the film damping co is added to the pair damping. Since both components scale with the same kinematic quantities, the combined damping inherits the staircase shape of the stiffness curve, and lubrication raises both coefficients at every instant of the cycle.
Point-Contact Elastohydrodynamic Model
The film is resolved with the standard isothermal point-contact formulation, with the entrainment direction along the minor axis of the contact ellipse. The Reynolds equation in its thin-film form,
$$\frac{\partial}{\partial x}\left(\frac{\rho h^3}{12\eta}\frac{\partial p}{\partial x}\right)+\frac{\partial}{\partial y}\left(\frac{\rho h^3}{12\eta}\frac{\partial p}{\partial y}\right)=u_s\frac{\partial(\rho h)}{\partial x}$$
is solved together with the film-thickness equation
$$h(x,y)=h_0+\frac{x^2}{2R_X}+\frac{y^2}{2R_Y}+\frac{2}{\pi E_x}\iint\frac{p(x’,y’)\,dx’\,dy’}{\sqrt{(x-x’)^2+(y-y’)^2}}$$
the Roelands viscosity–pressure relation
$$\eta=\eta_0\exp\left\{(\ln\eta_0+9.67)\left[-1+\left(1+5.1\times10^{-9}p\right)^{z}\right]\right\}$$
the Dowson–Higginsong density–pressure relation
$$\rho=\rho_0\left(1+\frac{0.6\times10^{-9}p}{1+1.7\times10^{-9}p}\right)$$
and the load balance
$$w_m=\iint p(x,y)\,dx\,dy$$
The domain is normalized by the contact half-widths, the Hertzian peak pressure and the ambient viscosity and density, which removes the several-orders-of-magnitude spread between the pressure, film and elasticity terms. The normalized Reynolds equation is discretized on a uniform grid with central differences for the pressure-gradient terms and a forward difference for the wedge term, and solved by finite-difference iteration with the elastic deformation computed from the influence-coefficient matrix. The two inputs that vary along the path of contact are computed at each mesh position: the entrainment speed, obtained by projecting the pinion and wheel surface velocities onto the tangent plane at the contact point and averaging them,
$$u_s=\frac{1}{2}\sqrt{\left(u_{z\xi}+u_{c\xi}\right)^2+\left(u_{z\psi}+u_{c\psi}\right)^2}$$
and the composite radius of curvature, which follows from the instantaneous radii of the two teeth measured from the pitch point.
For the present pair the entrainment speed increases monotonically along the path of contact, while the composite radius of curvature peaks near the pitch line and falls toward both ends. These two trends act in opposition — higher speed thickens the film, smaller radius of curvature thins it — which is why the film thickness at entry and exit of mesh does not follow a simple monotonic law.
Film Response to Speed and Viscosity
With the lubricant held at 0.02135 Pa·s (24.40 mm²/s at 40 °C, density 870 kg/m³, pressure–viscosity coefficient 2.082×10⁻⁸ Pa⁻¹), raising the pinion speed from 500 to 1500 r/min thickens the film substantially but changes the peak pressure comparatively little. The regime also shifts: at 500 r/min the maximum pressure coincides with the primary Hertzian peak, whereas at 1000 and 1500 r/min the secondary pressure spike in the outlet region becomes the highest point in the distribution.
| Speed (r/min) | Change in minimum film thickness vs. previous step | Change in maximum film pressure vs. previous step | Location of maximum pressure |
|---|---|---|---|
| 500 | — | — | Primary peak |
| 1000 | +39.016% | +0.584% | Secondary peak |
| 1500 | +27.980% | +5.229% | Secondary peak |
Viscosity acts more strongly than speed on film thickness at fixed speed. At 1500 r/min, moving from 0.00439 to 0.03190 Pa·s more than triples the minimum film thickness across two steps, while the peak pressure rises by roughly a tenth.
| Dynamic viscosity (Pa·s) | Change in minimum film thickness vs. previous step | Change in maximum film pressure vs. previous step | Location of maximum pressure |
|---|---|---|---|
| 0.00439 | — | — | Primary peak |
| 0.02135 | +235.101% | +4.652% | Secondary peak |
| 0.03190 | +22.046% | +6.231% | Secondary peak |
The pattern is consistent with classical point-contact scaling: film thickness is far more sensitive to viscosity and entrainment speed than to load, while the pressure distribution stays close to Hertzian except in the outlet spike. That last feature matters for dynamics, because the spike is a local, load-dependent stress concentration that the smooth Hertzian model cannot represent.
Extracting Film Stiffness and Film Damping
Three routes are available for converting a film solution into a stiffness. The average-film method divides the load increment by the change in the mean film thickness over the Hertzian zone; it suffers from too few usable data points and from boundary contamination. The global method treats each grid node as a parallel spring and sums the nodal contributions; it is accurate but expensive. The empirical-derivative method differentiates a closed-form film-thickness relation with respect to load. The Hamrock–Dowson point-contact relations
$$h_{\min}=3.63\,u_s^{0.68}G^{0.49}w_m^{-0.073}R_X^{0.64}\left(1-e^{-0.68R_Y/R_X}\right)$$
$$h_{c}=2.69\,u_s^{0.67}G^{0.53}w_m^{-0.067}R_X^{0.75}\left(1-0.61e^{-0.73R_Y/R_X}\right)$$
give a film stiffness that varies smoothly along the path of contact and is cheap to evaluate at every dynamic time step, which is decisive when the stiffness must be updated inside an iterative dynamic solution. The price is that the local outlet spike is smoothed out.
Film damping is obtained from an energy argument rather than from a derivative. A dimensionless sinusoidal load is applied to the contact, the resulting contact force is integrated over the cycle, and the area of the resulting hysteresis loop in the force–film-thickness plane equals the energy dissipated per cycle,
$$E_C=\oint f\,dh_C=\pi C A^2\Omega$$
which yields the dimensionless damping C. The loop is nearly elliptical, confirming that the film behaves as a viscous element in series with the elastic bodies. For engineering use the numerical result is replaced by a fit in the Moes parameters,
$$C=\frac{4.3\left(L_o/M_o\right)^{0.85}}{\left[1+3\left(L_o/M_o\right)^{1.2}\right]^{2}}$$
$$L_o=\alpha_{\mathrm{oil}}E_x\left(\frac{2\eta_0u_s}{E_xR_X}\right)^{1/4},\qquad M_o=\frac{w_m}{E_xR_X}\left(\frac{2\eta_0u_s}{E_xR_X}\right)^{-1/2}$$
which reproduces the numerical damping across the load and speed range of interest. Under normal running the film damping is several orders of magnitude below the structural pair damping and therefore has little effect on the total. The exception is instructive: under high-speed, heavily loaded or boundary-lubricated conditions the two become comparable, and the viscous dissipation of the film then modulates the combined mesh damping appreciably. A model that omits the film is therefore not wrong everywhere — it is wrong precisely in the regimes where surface distress is most likely.
Eight-Degree-of-Freedom Dynamics with a Film-Modified Transmission Error
The gear pair is represented by a lumped-mass system with three translational coordinates and one rotational coordinate per member, giving the vector {X1, Y1, Z1, θ1, X2, Y2, Z2, θ2}, with bearing stiffness and damping in all three radial directions at each shaft and a torque applied at each end. The dynamic transmission error along the line of action is the projection of the relative translational and torsional displacements onto the normal, plus the static transmission error, which is represented by its first harmonic at the mesh frequency.
The distinguishing feature of the formulation is the treatment of the film in that projection. In a dry model the mating surfaces are in direct contact, so the relative approach of the two bodies is entirely elastic. With a film present, the two surfaces are separated by the central film thickness hcenter, and this separation adds directly to the transmission error:
$$\lambda_n=\lambda_n^{\mathrm{dry}}+h_{\mathrm{center}}$$
The mesh force then follows as
$$F_m=k_t\,f(\lambda_n)+c_t\,\dot{\lambda}_n$$
where f(·) is the backlash function that clips the transmission error at the mean clearance, and kt and ct are the combined stiffness and damping defined above. This closes the loop: the transmission error sets the mesh force, the mesh force and the kinematics set the film, and the film sets both the transmission error and the stiffness and damping coefficients that produce the force.
Numerically the loop is broken by alternation. The eight equations of motion are integrated in dimensionless form with a fourth-order Runge–Kutta scheme from quiescent initial conditions, giving the composite radius of curvature, entrainment speed and mesh force at each time step. Those three quantities drive the elastohydrodynamic solution, from which the film stiffness, film damping and central film thickness are extracted. The updated coefficients are returned to the dynamic model and the cycle repeats until the relative change in the cycle-summed stiffness and damping both fall below 0.001. Typical convergence requires only a few outer iterations because the film responds to the mesh force much faster than the mesh force responds to the film.
Predicted Dynamic Response With and Without the Film
Simulations were run at 1500 r/min with an input torque of 100 N·m, an equivalent mass of 5.148 kg for the pinion and 165.65 kg for the wheel, and radial bearing stiffness of 6.56×10⁷ N/m with 1.8×10⁵ N·s/m damping. Compared with dry contact, the lubricated mesh shows higher stiffness and higher damping at every instant, and the peak-to-peak amplitudes of every response quantity fall.
| Quantity | Dry contact | Lubricated | Reduction |
|---|---|---|---|
| Dynamic transmission error | 1.85299×10⁻⁵ m | 1.51471×10⁻⁵ m | 18.256% |
| Single-tooth mesh force | 1928.0189 N | 1658.6190 N | 13.973% |
| Displacement along line of action | 5.12421×10⁻⁵ m | 4.11541×10⁻⁵ m | 19.687% |
| Velocity along line of action | 1.378908×10⁻³ m/s | 1.128098×10⁻³ m/s | 18.189% |
| Acceleration along line of action | 45.0901 m/s² | 36.9326 m/s² | 18.091% |
The mechanism is straightforward once the coefficients are separated. The film raises the effective mesh stiffness, which shortens the interval over which the contact loads and unloads, and it adds viscous dissipation, which damps the oscillation that would otherwise persist across the mesh cycle. The two effects act in the same direction on the force amplitude, so the transmission error contracts and the displacement, velocity and acceleration responses contract with it. The force reduction (14%) is smaller than the transmission-error reduction (18%) because the force is the product of a larger stiffness and a smaller displacement, and the two partly offset.
Test-Stand Validation
Measurements were taken on a dedicated spiral bevel gearbox test stand built around the same 23/115 pair. A drive motor feeds the test gearbox through a planetary speed increaser that extends the achievable input speed, and an eddy-current brake loads the output. Two IEPE piezoelectric triaxial accelerometers are mounted on the housing adjacent to the wheel shaft and the pinion shaft, with axial sensitivities near 10 mV/(m/s²), frequency response from 0.5 Hz to 5 kHz on the X axis and to 7 kHz on the Y and Z axes, and an operating range of −40 to +120 °C. Signals were recorded on an eight-channel acquisition unit at a 10 kHz sampling rate, with channels 1–3 assigned to the pinion-side sensor and 5–7 to the wheel-side sensor; the remaining channels were disabled to raise the effective resolution.
The lubricant was a synthetic aviation oil (Great Wall 4106) with kinematic viscosity of 24.40 mm²/s at 40 °C, 5.02 mm²/s at 100 °C and 8652 mm²/s at −40 °C, a flash point of 250 °C and a pour point below −54 °C. Because no measured value exists at 20 °C, the viscosity there was extrapolated with a double-logarithmic temperature relation and, after conversion with a density of 875 kg/m³, gave 0.03190 Pa·s. Viscosity was varied by controlling oil-sump temperature rather than by changing oil grade, which keeps the additive package constant but introduces a temperature effect on the housing and bearings as well.
| Sump temperature (°C) | Dynamic viscosity (Pa·s) | Input speed (r/min) |
|---|---|---|
| 100 | 0.00439 | 500 / 1000 / 1500 / 2000 |
| 40 | 0.02135 | 500 / 1000 / 1500 / 2000 |
| 20 | 0.03190 | 500 / 1000 / 1500 / 2000 |
| Dry (no lubricant) | 0 | 500 / 1000 / 1500 / 2000 |
Speed is the dominant variable. The RMS of the pinion-side X acceleration rises by a factor of roughly 3.5 between 500 and 2000 r/min, and the model reproduces both the trend and the order of magnitude, with a maximum deviation of 19.305% at 2000 r/min. The corresponding comparison for viscosity is much less satisfactory: the model predicts a monotonic decrease in RMS with increasing viscosity, whereas the measurements scatter without a clear trend, the largest discrepancy being 28.229% in the dry case. A plausible explanation is that under these conditions the film is not fully developed over the whole cycle, so the assumption of a continuous lubricating film — on which both the film stiffness and the film damping depend — is only partly satisfied.
| Input speed (r/min) | Dry | 0.00439 Pa·s | 0.02135 Pa·s | 0.03190 Pa·s |
|---|---|---|---|---|
| 500 | 1.82099 | 1.75888 | 1.45713 | 1.04400 |
| 1000 | 2.84958 | 2.62270 | 2.32551 | 2.73984 |
| 1500 | 5.51987 | 5.35679 | 5.80646 | 5.48040 |
| 2000 | 8.32360 | 8.00403 | 7.81259 | 7.75883 |
Two observations survive the scatter. First, lubricated running is consistently gentler than dry running at the extremes of the speed range: at 500 and 2000 r/min the highest-viscosity oil lowers the RMS by 42.67% and 6.79% respectively relative to dry contact. Second, the sensitivity to viscosity is strongly speed dependent. At 500 r/min, moving from 0.00439 to 0.03190 Pa·s cuts the RMS by 40.64%; at 2000 r/min the same change yields only 3.06%. At 1000 and 1500 r/min no monotonic ordering appears at all. A consistent reading of these numbers is that viscosity matters most when the entrainment speed is low and the film is hardest to establish, and that once the film is adequate at high speed, adding more viscosity buys little in vibration terms while costing churning losses.
What the Evidence Supports and Where It Stops
The coupled formulation rests on a defensible chain: the film modifies the contact compliance and adds viscous dissipation; those changes raise the effective mesh stiffness and damping; the raised coefficients reduce the dynamic transmission error and the force and motion amplitudes that follow from it. The predicted reductions of 18% in transmission error and 14% in mesh force are of the size that would matter in a vibro-acoustic assessment of a helicopter input stage, and they rest on a lubrication solution that resolves the actual, time-varying entrainment speed and curvature of a spiral bevel contact rather than a quasi-static equivalent.
The chain also has visible weak points. The film stiffness used in the dynamic loop comes from a closed-form fit rather than from the full numerical solution, so the outlet pressure spike — the feature most sensitive to load and most relevant to surface distress — is smoothed away. The thermal problem is decoupled: viscosity is treated as a function of ambient pressure and a fixed inlet temperature, although in a high-speed aircraft gearbox the film temperature rise across the contact can be large and can reduce the film thickness appreciably. The experimental programme covers input speeds up to 2000 r/min only, well below the regime where the film damping becomes comparable with the structural damping; the most interesting part of the model, from the standpoint of the earlier argument, is therefore the least validated. Viscosity was varied by temperature control rather than by changing oil grade, which confounds the viscosity effect with thermal expansion and clearance changes.
The larger limitation is scope. A gearbox is a rotor–bearing–housing system, and the mesh is one excitation path among several. Bearing clearance nonlinearities, housing flexibility and load fluctuations from the rotor system all shape the measured acceleration spectrum, and none of them appears in an eight-degree-of-freedom model with rigid supports. The agreement of 19% between model and measurement at the highest speed is respectable for a model of this class, but it also marks the level at which the neglected dynamics begin to dominate. Closing that gap requires embedding the lubricated mesh, with its film-modified stiffness and damping, into a flexible housing and bearing model — a step that the present results justify but do not attempt.