Loss of lubrication is a credible failure mode in helicopter drivetrains. If the oil system is breached by hostile weather, ballistic damage, or a seal failure, the spiral bevel gears at the high-speed input of the main gearbox or in the intermediate and tail gearboxes must continue to transmit torque with progressively less oil. The contact begins in elastohydrodynamic lubrication (EHL), passes through mixed and boundary lubrication, and eventually reaches dry friction. Friction power raises tooth temperature, thermal growth reduces backlash, the film collapses, and scuffing or seizure can follow. Helicopter transmission designs commonly require at least 30 min of dry running after oil loss, while advanced designs target 60 min or more. Direct dry-running tests are expensive, can damage the test rig, and are difficult to instrument at the tooth contact; tooth-film thickness and local tooth temperature are especially hard to measure. A predictive model therefore has to connect micro-scale film behavior to the macro-scale thermal response of the whole gearbox.
A three-level coupling problem
The physics spans several orders of magnitude in both time and space. At the tooth surface, EHL film formation and pressure generation occur on time scales of microseconds and length scales of micrometres. At the tooth-meshing level, friction heat, windage, and surface convection occur on time scales of tenths of a millisecond and length scales of tenths of a millimetre. At the gear-system level, transient heat conduction and convection through gears, shafts, bearings, oil-air mixture, and housing occur on time scales of tenths of a second to seconds and length scales of decimetres. A single transient simulation covering all three levels is impractical, so the levels are coupled through selected parameters rather than resolved in one mesh.
| Level | Typical time scale | Typical length scale | Main outputs passed to other levels |
|---|---|---|---|
| Tooth-surface EHL | 10−6–10−5 s | 10−6 m | Film thickness, pressure, friction coefficient, friction heat source |
| Tooth meshing | 10−4–10−3 s | 10−4 m | Meshing friction power, windage power, surface heat-transfer coefficients |
| Gear-system heat transfer | 10−2–100 s | 10−1 m | Component temperatures, internal oil-air temperature, housing temperature |
The coupling loop is directional but iterative. The EHL model supplies a time-dependent friction coefficient to the meshing model; the meshing model converts friction and windage into heat sources and surface heat-transfer coefficients; the system model predicts the temperature of the internal oil-air mixture and the gearbox components. The internal air temperature is then returned to the meshing model as its environmental temperature and to the EHL model as the initial oil temperature. In this way, the model does not treat the gearbox interior as ambient room temperature, which would overestimate heat removal and underpredict tooth temperature.
Tooth-contact EHL under full lubrication
Spiral bevel gears operate with local conjugate point contact. The contact ellipse changes size and orientation during meshing, and the entrainment velocity is generally not aligned with either the major or minor axis of the ellipse. Load, sliding speed, entrainment speed, and surface curvature all vary with meshing position. For a representative Gleason spiral bevel pair with a 31-tooth pinion and a 46-tooth wheel, a module of 6.3 mm, a face width of 57.4 mm, a mean spiral angle of 35°, and a pinion speed of 5,000 r/min, the sampled contact points gave a minor axis of about 0.82–0.89 mm, a major axis of about 1.46–1.53 mm, and an entrainment velocity of about 4.30–6.98 m/s. The angle between the entrainment velocity and the minor axis was approximately 86°, so the contact is highly skewed relative to the rolling direction.
The local EHL model solves the Reynolds equation in the contact coordinate system:
$$ \frac{\partial}{\partial x}\left(\frac{\rho h^{3}}{\eta}\frac{\partial p}{\partial x}\right)+\frac{\partial}{\partial y}\left(\frac{\rho h^{3}}{\eta}\frac{\partial p}{\partial y}\right)=12u\frac{\partial}{\partial x}\left(\frac{\rho h}{\eta}\right)+12v\frac{\partial}{\partial y}\left(\frac{\rho h}{\eta}\right) $$
where x and y are the minor- and major-axis directions, p is the film pressure, h is the film thickness, ρ is the oil density, η is the effective viscosity, and u and v are the entrainment-velocity components. The film thickness includes the geometric gap and the elastic deformation of the contacting bodies:
$$ h=h_{0}+\frac{x^{2}}{2R_{x}}+\frac{y^{2}}{2R_{y}}+\delta(x,y) $$
$$ \delta(x,y)=\frac{2}{\pi E’}\iint_{\Omega}\frac{p(x’,y’)}{\sqrt{(x-x’)^{2}+(y-y’)^{2}}}\,dx’\,dy’ $$
Here h0 is the central film thickness, Rx and Ry are the effective radii of curvature, δ is the elastic deformation, E′ is the combined elastic modulus, and Ω is the contact area. The viscosity is described by a Roelands pressure–temperature relation, the density by a pressure–temperature relation, and the oil-film energy equation is solved together with the Reynolds equation because friction heating changes viscosity inside the contact.
The friction coefficient is obtained by integrating the wall shear stress over the contact and dividing by the normal load:
$$ f_{\mathrm{EHL}}=\frac{\iint_{\Omega}\tau(x,y)\,dx\,dy}{\iint_{\Omega}p(x,y)\,dx\,dy} $$
For the example gear pair, the full-film friction coefficient over the sampled meshing points was approximately 0.010–0.030. This low friction coefficient is characteristic of a fully developed EHL film under high contact pressure, where the film is of micrometre order and the pressure is of gigapascal order. The calculation requires a strongly coupled numerical solution; multigrid methods are used for the elastic deformation, and progressive mesh refinement is used for the Reynolds equation.
From full film to starved contact
During oil loss, the contact is no longer fully replenished. The film becomes thinner at each meshing cycle, and the lubrication state evolves from full EHL to mixed lubrication, boundary lubrication, and finally dry friction. The regime can be indexed by the film ratio
$$ \lambda=\frac{h_{\min}}{R_{a}} $$
where hmin is the minimum film thickness and Ra is the average surface roughness. In the full-film regime, the surfaces are separated by oil and the load is carried by the film. In the mixed regime, part of the load is carried by the film and part by asperities and boundary films. In the boundary regime, asperity contact and boundary films carry almost all the load. In dry friction, no continuous oil film remains and the metal–boundary-film contact carries the load.
A film-retention parameter is used to describe the decay of the central film with time:
$$ R(n)=\frac{h_{c}}{h_{\mathrm{EHL}}}\approx(\omega t)^{-1/\gamma} $$
where hc is the central film thickness after time t, hEHL is the fully flooded central film thickness, ω is the angular speed, and γ is a film-decay exponent. For the example, γ = 3. With an average roughness of Ra = 0.2 µm and a fully flooded central film of hEHL = 0.5 µm, the onset of boundary lubrication is estimated at hc ≈ 0.16 µm, corresponding to R(n) ≈ 0.32. The model therefore treats R(n) > 0.32 as the full-film regime, 0.01 < R(n) < 0.32 as the mixed/boundary regime, and R(n) < 0.01 as the dry-friction regime.
The load fraction carried by the film can be approximated by
$$ \epsilon=\frac{1.21\lambda^{0.64}}{1+0.37\lambda^{1.26}} $$
and the friction coefficient is blended between the full-film value, a boundary-lubrication value of about 0.11, and a dry-friction value of about 0.6. As the film fraction falls, the friction coefficient rises from the EHL range of 0.01–0.03 toward the boundary and dry values. This time-dependent friction coefficient is the key input for the friction-heat source used in the meshing and system models. The blending assumes that asperity and boundary-film contact is distributed uniformly over the nominal contact; that is a modeling approximation, not a direct measurement of every asperity interaction.
Meshing-scale heat generation and windage
At the tooth-meshing level, the dominant heat sources are sliding friction at the contact and windage from the rotating gear. Because the meshing time scale is very short, the tooth-level calculation is treated as a steady CFD problem rather than a transient one. The friction power is obtained by summing the product of friction coefficient, normal load, and sliding velocity over the contact points:
$$ P_{\mathrm{mesh}}=\sum_{i}f_{i}F_{i}v_{s,i} $$
where fi is the local friction coefficient, Fi is the normal load at point i, and vs,i is the local sliding velocity. To obtain windage, a thin fluid region is built around the teeth and a moving-reference-frame (MRF) CFD model is used. The dimensionless drag-torque coefficient is
$$ C_{m}=\frac{M}{0.5\rho_{\mathrm{ref}}v_{\mathrm{ref}}^{2}A_{\mathrm{ref}}L_{\mathrm{ref}}} $$
where M is the drag torque, and ρref, vref, Aref, and Lref are reference density, velocity, area, and length. The windage power is then
$$ P_{w}=M\omega $$
The CFD solution also provides surface heat-transfer coefficients for the tooth flanks, tooth sides, gear end faces, and nearby bearing surfaces. Windage is not negligible at high speed: the dimensionless drag coefficient depends on the tooth surface—convex flank, concave flank, and other surfaces—and increases with rotational speed. The meshing model passes the friction heat, windage heat, and surface heat-transfer coefficients to the system-level thermal model.
System-level transient thermal response
The gear-system model resolves transient heat conduction and convection over the gears, shaft, bearings, internal oil-air mixture, and housing. To keep the computational cost manageable, the two gears can be represented as equivalent disks at this level, while the heat sources and surface heat-transfer coefficients are inherited from the meshing-level CFD. The housing loses heat to the surrounding air by natural convection, approximated by
$$ h=0.3\,Re^{0.57}\frac{k}{L} $$
where Re is the Reynolds number based on the housing scale, k is the air thermal conductivity, and L is the housing diameter. Under fully lubricated steady conditions, the system-level calculation gives a peak temperature of about 91.6 °C. The gears are the hottest components, the pinion is hotter than the wheel, and the housing is hottest near the bearing outer race. The internal oil-air mixture is about 58.7 °C, which is significantly above room temperature. This value is returned to the meshing model as its environment and to the EHL model as the initial oil temperature.
Under oil loss, the system temperature no longer has a steady, low-temperature baseline. The predicted maximum surface temperatures of the pinion and wheel are listed in the table below. The pinion heats faster and reaches higher temperatures because it has the higher speed and more severe sliding conditions in this gear pair.
| Time after oil loss (s) | Pinion maximum surface temperature (°C) | Wheel maximum surface temperature (°C) |
|---|---|---|
| 600 | 108 | 95 |
| 1,200 | 141 | 122 |
| 1,800 | 176 | 152 |
| 2,400 | 209 | 180 |
| 3,000 | 239 | 205 |
| 3,600 | 266 | 227 |
| 4,200 | 289 | 246 |
| 4,800 | 309 | 263 |
Predicted dry-running time and failure mode
The coupled calculation shows a characteristic thermal history. Friction heat rises as the film becomes starved, but the gear-system temperature does not follow the friction curve instantaneously. The temperature first rises steeply, then more slowly, and finally approaches a plateau. The thermal lag comes from the heat capacity of the gears and housing and from the finite rate of heat conduction from the contact into the gear body. The pinion remains the critical component: at 1 h its maximum surface temperature is about 266 °C, and at 1.33 h it reaches about 309 °C.
Boundary lubricants provide a useful failure criterion. Adsorbed boundary films typically break down near 150 °C, while chemical reaction films can survive to approximately 300 °C. Once the pinion surface exceeds about 300 °C, the chemical film can no longer protect the contact, and the tooth pair enters a complete dry-friction state. After that point, wear and scuffing progress rapidly and the gear can fail. The coupled model therefore predicts failure at approximately 1.5 h after oil loss, with the pinion as the first component to lose effective lubrication.
This result should be read with the model assumptions in mind. The bearings are treated as fully lubricated in the calculation, and the housing wall is heated mainly by the internal air rather than by direct contact with the gears. In an actual helicopter transmission, an emergency oil reservoir may keep a small amount of oil on the teeth after the main supply is lost, so the real dry-running time can be longer than the conservative model prediction. The value of the coupled approach is not a single guaranteed life number, but a mechanism-based way to rank designs, identify the critical component, and set test conditions. It connects film-scale lubrication, tooth-scale heat generation, and gearbox-scale thermal response without pretending that one mesh can resolve all three.
Design implications for spiral bevel gears
For spiral bevel gears in high-speed helicopter transmissions, the analysis points to three design levers. First, the pinion is thermally critical, so surface finish, boundary-film chemistry, and cooling near the pinion should receive more attention than a simple bulk-temperature check would suggest. Second, the internal air temperature is a real boundary condition, not a cosmetic detail: using room temperature for the meshing CFD would misrepresent heat removal and shift the predicted failure time. Third, the thermal lag between friction heating and component temperature means that a short burst of high friction does not immediately produce peak metal temperature, while a sustained starved condition does. The coupled model captures that lag by using a steady meshing calculation inside a transient system calculation.
The same framework can be used to compare surface roughness values, alternative boundary lubricants, emergency-oil strategies, and housing ventilation concepts. Its main limitations are the assumed uniformity of asperity contact, the empirical film-decay exponent, and the uncertainty in convective boundary conditions. These are appropriate targets for validation with instrumented rig tests and for refinement with measured oil-air temperatures and post-test tooth-surface observations. Within those limits, the multi-level model provides a practical way to estimate the dry-running window of spiral bevel gears and to explain why the pinion, rather than the more heavily loaded wheel, can become the first point of failure.