Spiral bevel gears transmit torque between intersecting or non-orthogonal shafts and are widely used in helicopter main reducers, automotive axles, and industrial drives. Their loaded contact pattern is sensitive to shaft and bearing deflection, so mesh stiffness, static transmission error, backlash, and tooth-surface friction are coupled. Friction is often treated as secondary because the friction force is only a fraction of the normal mesh force. However, its point of action moves across the tooth, its direction reverses at the pitch point, and it generates a moment about the gear axis. These features make friction an internal excitation rather than a simple damping term.
Eight-degree-of-freedom mesh model
The dynamic formulation treats each gear as a rigid disk with three translations and one rotation, giving eight degrees of freedom: \(X_p,Y_p,Z_p,\theta_p\) for the pinion and \(X_g,Y_g,Z_g,\theta_g\) for the wheel. Rocking motion is neglected. Shaft and bearing flexibility is represented by equivalent support stiffness and damping acting at the mid-face, and the mesh point is also taken at the mid-face. The wheel local frame is rotated by the shaft angle \(\Sigma\) about the pinion \(Y\)-axis, which allows non-orthogonal arrangements such as the \(75^\circ\) shaft angle used in helicopter main reducers.
The normal mesh force combines viscous mesh damping and a nonlinear elastic term:
$$F_n=c_m\dot{\lambda}_n+k_h(t)f(\lambda_n)$$
where \(c_m\) is the mesh damping, \(k_h(t)\) is the time-varying mesh stiffness, and \(\lambda_n\) is the normal relative displacement at the mesh point. The backlash function is
$$f(\lambda_n)=
\begin{cases}
\lambda_n-b_n, & \lambda_n \gt b_n,\\
0, & |\lambda_n|\le b_n,\\
\lambda_n+b_n, & \lambda_n \lt -b_n,
\end{cases}$$
with \(b_n\) denoting half of the normal backlash. Direction cosines \(a_i\) project the normal force onto the local coordinate axes. Friction is introduced as
$$F_f=\mu F_n,\qquad F_{fi}=n_{fi}F_f,\qquad M_{fi}=r_{fi}F_f,\qquad i=p,g$$
where \(\mu\) is the friction coefficient, \(n_{fi}\) is the friction-force direction cosine in the local frame, and \(r_{fi}\) is the corresponding moment arm. The pinion friction opposes the relative sliding velocity, while the wheel friction acts in the sliding direction. Because the sliding velocity changes sign at the pitch point, the friction force and its moment alternate during each mesh cycle.
The normal relative displacement includes translational and torsional contributions and the static transmission error:
$$\lambda_n=-a_1X_p-a_2Y_p-a_3Z_p-a_4X_g-a_5Y_g-a_6Z_g+a_3r_p\theta_p+a_6r_g\theta_g-e_n(t)$$
Here \(r_p\) and \(r_g\) are the mesh-point radii, and \(e_n(t)\) is the static transmission error. For the present case, the error is represented by a zero mean and a first-order harmonic, \(e_n(t)=e_1\cos(\omega_e t)\), where \(\omega_e\) is the mesh frequency. The translational and rotational equations take the general form
$$m_i\ddot{\mathbf{r}}_i+\mathbf{c}_i\dot{\mathbf{r}}_i+\mathbf{k}_i\mathbf{r}_i=\mathbf{F}_{ni}+\mathbf{F}_{fi},\qquad
J_i\ddot{\theta}_i+M_{ni}+M_{fi}=T_i$$
for \(i=p,g\). The two torsional coordinates are not independent because rigid-body rotation is present; \(\lambda_n\) is therefore used as an additional generalized coordinate to eliminate the rigid-body component.
Time-varying stiffness and mixed-lubrication friction
The mesh stiffness is obtained from a three-tooth loaded tooth contact analysis. Because the number of tooth pairs in contact changes periodically, the stiffness is a periodic function and is fitted by a Fourier series:
$$k_h(t)=k_m+\sum_{i=1}^{3}\left[k_{ai}\cos(i\omega_e t+\phi_{ki})+k_{bi}\sin(i\omega_e t+\phi_{ki})\right]$$
where \(k_m\) is the mean mesh stiffness, \(k_{ai}\) and \(k_{bi}\) are harmonic amplitudes, and three harmonics are retained to include sub-harmonic and super-harmonic effects.
For gear teeth in mixed lubrication, the friction coefficient is not a fixed material constant. It depends on the share of load carried by the lubricating film and by asperity contacts:
$$\mu_{mb}=f_{\rm l}^{1.2}\mu_{ml}+(1-f_{\rm l})\mu_{bl}$$
where \(\mu_{ml}\) is the full-film elastohydrodynamic friction coefficient, \(\mu_{bl}\) is the boundary-lubrication coefficient, typically in the range \(0.1\) to \(0.2\), and \(f_{\rm l}\) is the load-share fraction carried by the film. The full-film term depends on the normal load per unit face width, sliding speed, oil dynamic viscosity, composite surface roughness, pitch radius, and composite curvature radius. Under starved lubrication, the effective coefficient is written as a load-share combination of film and asperity contact:
$$\mu_f=k_c\mu_{ml}+k_b\mu_{bl},\qquad k_c+k_b=1$$
where \(k_c\) and \(k_b\) are the load fractions carried by the elastohydrodynamic film and by asperity peaks. This formulation captures the transition from normal mixed lubrication to oil-starved and dry running without treating friction as an arbitrary constant.
Baseline response of a helicopter main-reducer gear pair
The case considered is a spiral bevel gear pair from a helicopter main reducer, rated at \(1{,}200\ \mathrm{kW}\) input power and \(20{,}000\ \mathrm{r/min}\) input speed. The shaft angle is \(75^\circ\), so the axes are non-orthogonal. The basic geometry is summarized in Table 1.
| Parameter | Pinion | Wheel |
|---|---|---|
| Number of teeth | 29 | 76 |
| Module | \(3.75\ \mathrm{mm}\) | \(3.75\ \mathrm{mm}\) |
| Pressure angle | \(20^\circ\) | \(20^\circ\) |
| Spiral angle | \(30^\circ\) | \(30^\circ\) |
| Face width | \(40\ \mathrm{mm}\) | \(40\ \mathrm{mm}\) |
| Hand | Left | Right |
| Shaft angle | \(75^\circ\) | |
| Rated input power | \(1{,}200\ \mathrm{kW}\) | |
| Input speed | \(20{,}000\ \mathrm{r/min}\) | |
The eight equations are integrated with a variable-step Runge-Kutta method. The first several hundred mesh cycles are discarded to remove the initial transient and obtain the steady-state response. Under normal mixed lubrication, with \(\mu=0.0244\), the normal relative displacement fluctuates between \(80.0458\ \mu\mathrm{m}\) and \(85.4185\ \mu\mathrm{m}\). The dynamic normal mesh force indicates single-sided impact rather than repeated backlash crossing. The phase portrait is periodic, and the FFT spectrum is dominated by a narrow set of frequency components.
Effect of lubrication state and friction coefficient
Three lubrication conditions are compared: normal mixed lubrication, starved lubrication, and dry running. The corresponding friction coefficients and root-mean-square responses are listed in Table 2. The normal relative displacement range is also given to show the shift in the response band.
| Lubrication state | Friction coefficient \(\mu\) | RMS vibration displacement (\(\mu\mathrm{m}\)) | RMS normal mesh force (\(\mathrm{N}\)) | Normal relative displacement range (\(\mu\mathrm{m}\)) |
|---|---|---|---|---|
| Normal mixed lubrication | 0.0244 | 82.937 | 15,070 | 80.0458–85.4185 |
| Starved lubrication | 0.050 | 82.939 | 15,071 | 80.0461–85.4210 |
| Dry running | 0.52 | 82.977 | 15,092 | 80.0457–85.4722 |
As the oil supply decreases, the friction coefficient rises from \(0.0244\) to \(0.52\), and both the vibration displacement and the normal mesh force increase. The trend is monotonic but modest: the RMS displacement increases by only about \(0.05\%\), and the RMS normal mesh force increases by about \(22\ \mathrm{N}\), or approximately \(0.15\%\). The dry-running range extends to \(85.4722\ \mu\mathrm{m}\), compared with \(85.4185\ \mu\mathrm{m}\) under normal lubrication. In this model, friction enters the normal-direction equations through the directional terms \(a_i+\mu n_{fi}\), so its influence is bounded by the friction coefficient relative to the normal-force projection.
Effect of contact ratio through the contact trace angle
Contact ratio is controlled here without changing the gear macro-geometry. Active tooth-surface design is used to vary the angle between the wheel contact trace and the root cone. The allowed range is \(0^\circ\) to \(90^\circ\); three representative angles, \(15^\circ\), \(45^\circ\), and \(75^\circ\), give design contact ratios of \(2.125\), \(1.417\), and \(1.250\), respectively. As the contact trace angle increases, the contact ratio decreases and the mesh stiffness also decreases.
| Contact trace angle to root cone (\(^\circ\)) | Design contact ratio \(\varepsilon\) | RMS vibration displacement (\(\mu\mathrm{m}\)) | RMS normal mesh force (\(\mathrm{N}\)) | Normal relative displacement range (\(\mu\mathrm{m}\)) |
|---|---|---|---|---|
| 15 | 2.125 | 83.612 | 14,735 | 80.0112–87.2067 |
| 45 | 1.417 | 89.263 | 14,829 | 86.5481–91.8493 |
| 75 | 1.250 | 92.085 | 14,865 | 88.7469–95.3988 |
The contact-ratio effect is considerably stronger than the friction-coefficient effect in these computed cases. Reducing the design contact ratio from \(2.125\) to \(1.250\) raises the RMS vibration displacement from \(83.612\ \mu\mathrm{m}\) to \(92.085\ \mu\mathrm{m}\), an increase of about \(10\%\), while the RMS normal mesh force rises from \(14{,}735\ \mathrm{N}\) to \(14{,}865\ \mathrm{N}\), slightly less than \(1\%\). The displacement band also shifts upward, from \(80.0112\)–\(87.2067\ \mu\mathrm{m}\) at \(15^\circ\) to \(88.7469\)–\(95.3988\ \mu\mathrm{m}\) at \(75^\circ\). Lower contact ratio means fewer tooth pairs share the load and the mesh stiffness is lower and more abruptly transferred between teeth, which amplifies the dynamic response.
Design implications for spiral bevel gears
For high-speed, high-power spiral bevel gear transmissions, lubrication state and contact-path design affect the dynamic response through different mechanisms. Friction changes the tangential force and friction moment at the mesh point, and dry running raises both the friction coefficient and the normal mesh force. However, the computed sensitivity to friction is small compared with the sensitivity to contact ratio. Maintaining an adequate oil film is necessary for thermal and wear reasons, but it is not sufficient to suppress mesh vibration if the contact path produces a low contact ratio.
The contact trace angle provides a design lever: a smaller angle to the root cone preserves a higher contact ratio, lower mesh stiffness variation, and lower RMS vibration in the reported configuration. This benefit must be balanced against other loaded contact metrics, such as contact-pattern location, edge loading, and sensitivity to misalignment. The results also depend on the assumptions of the eight-degree-of-freedom model: rigid gear disks, lumped shaft and bearing flexibility, mid-face mesh representation, and numerical integration without housing or thermal coupling. Within those limits, the comparison shows that contact-ratio control is the stronger of the two factors examined, while tooth-surface friction remains a secondary but non-negligible internal excitation in the oil-starved and dry-running conditions.