Spiral bevel gears are widely used in helicopter main reducers because they transmit high torque between intersecting shafts with smooth, gradual tooth engagement. In service, the main reducer is not only a power path but also a vibration source: time-varying mesh stiffness, backlash, alternate tooth engagement, and transmission error act together, and the resulting housing motion can propagate into the airframe and rotor system. The design target is therefore a contact pattern that carries load without stress concentrations and a mesh excitation spectrum that is as low as practical.
Spiral bevel gears are not fully conjugate. Their tooth surfaces are deliberately modified so that contact is localized and remains tolerant to misalignment. Local synthesis defines the contact path and transmission error near a reference point through a small set of tooth-surface parameters. For a non-orthogonal helicopter gear pair, the practical question is how those parameters should be selected when the objective is vibration reduction rather than contact pattern alone.
A representative case is a civil-helicopter main reducer with a 95.5° shaft angle, an 11-tooth right-hand pinion, and a 49-tooth left-hand gear. The pair operates at 106.647 kW input power and 2097 r/min pinion speed. Its non-orthogonal layout, high torque density, and strict vibration requirements make it a useful example of how tooth-surface microgeometry can be treated as a source-control variable.
Mesh Excitation and the Eight-Degree-of-Freedom Model
The 95.5° shaft angle is not a cosmetic detail. The shaft angle and the pitch-cone angles determine how the normal mesh force is resolved into the local \(x\), \(y\), and \(z\) directions of each gear. In the lumped model, the gear bodies are treated as rigid conical disks, shaft and bearing flexibility is lumped at the mid-face, and tooth friction and rocking are neglected. Each gear retains three translations and one rotation about its own axis, giving eight degrees of freedom:
$$\left\{X_p,Y_p,Z_p,\theta_p,X_g,Y_g,Z_g,\theta_g\right\}^{T}$$
At the mid-face, the normal mesh force is
$$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 along the mesh line. The backlash function is
$$f(\lambda_n)=\begin{cases}
\lambda_n-b_n, & \lambda_n>b_n,\\
0, & |\lambda_n|\le b_n,\\
\lambda_n+b_n, & \lambda_n<-b_n,
\end{cases}$$
with \(b_n\) equal to half the normal backlash. The relative displacement along the mesh line is composed of the translational and rotational displacements of the two gears:
$$\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 \(a_1,\ldots,a_6\) are direction cosines of the normal force in the two local coordinate frames, \(r_p\) and \(r_g\) are the mesh-point radii, and \(e_n(t)\) is the transmission error. The transmission error is represented by a Fourier series:
$$e_n(t)=e_0+\sum_{i=1}^{\infty}e_i\cos(2\pi i f_e t+\phi_i)$$
where \(e_0\) is the mean error, \(e_i\) is the amplitude of the \(i\)-th harmonic, and \(\phi_i\) is its phase. The mesh frequency \(f_e\) is set by the pinion speed and tooth count:
$$f_e=\frac{n_1z_1}{60}$$
The equations of motion combine these mesh forces with support stiffness and damping, input torque, and load torque. They are made dimensionless before integration because the stiffness, damping, and displacement terms differ by orders of magnitude. The steady-state response is obtained after discarding the initial transient.
Case Gear Set and Analysis Chain
The gear data for the case are summarized below.
| Parameter | Pinion | Gear |
|---|---|---|
| Number of teeth | 11 | 49 |
| Module | 6.285 mm | 6.285 mm |
| Pressure angle | 20° | 20° |
| Spiral angle | 30° | 30° |
| Face width | 46.584 mm | 46.584 mm |
| Shaft angle | 95.5° | |
| Hand | Right | Left |
The 11/49 pair gives a nominal reduction ratio of 4.45:1. At 2097 r/min, the mesh frequency is approximately 384.5 Hz. The analysis chain has four parts. Tooth contact analysis gives the contact ratio for a candidate tooth surface. A three-dimensional finite-element model in ANSYS gives discrete mesh stiffness values, which are fitted by a Fourier series to form \(k_h(t)\). The eight-degree-of-freedom equations are integrated with the Runge-Kutta method. The RMS of the steady-state vibration displacement is then used as the optimization objective. Because the tooth-surface parameters do not appear explicitly in the dynamic equations, they enter through contact ratio and single-tooth/composite mesh stiffness.
Tooth-Surface Design Variables and Optimization Setup
Local synthesis controls the meshing geometry near the reference point with three parameters used here as design variables:
- Contact-trace angle \(\eta\): the angle between the contact trace on the gear and the tooth-height direction, bounded by \(0^\circ<\eta<90^\circ\). It sets the orientation of the contact path across the tooth.
- Transmission-ratio first derivative \(m’_{pg}\): the local slope of the transmission ratio at the reference point. On the working flank it is negative, with the bound \(-0.05<m'_{pg}<0\). curve.
- Contact-ellipse ratio \(K_D\): the ratio of the contact-ellipse semi-major axis to the face width, bounded by \(0.15<k_d<0.25\). available="" contact-ellipse="" controls="" face="" it="" li="" proportion="" relative="" the="" to="" width.
</k_d
</m'_{pg}
The three variables are coupled. The optimization problem is
$$\begin{aligned}
\min_{\eta,m’_{pg},K_D}\quad & V(\eta,m’_{pg},K_D)\\
\text{subject to}\quad & 0^\circ<\eta<90^\circ,\\
& -0.05<m'_{pg}<0,\\ &="" 0.15<k_d<0.25,="" \end{aligned}$$
where \(V\) is the RMS value of the steady-state vibration displacement of the transmission system. A genetic algorithm is used to search the bounded space. The settings are listed below.
| Item | Value |
|---|---|
| Contact-trace angle \(\eta\) | 0°–90° |
| Transmission-ratio first derivative \(m’_{pg}\) | −0.05 to 0 |
| Contact-ellipse ratio \(K_D\) | 0.15–0.25 |
| Population size | 10 |
| Generations | 10 |
| Crossover probability | 0.8 |
| Mutation probability | 0.07 |
| Encoding | Binary |
| Crossover | Single-point |
Optimization Results
The first and final iterations are compared below.
| Iteration | \(\eta\) (°) | \(m’_{pg}\) | \(K_D\) | RMS displacement (μm) | Change |
|---|---|---|---|---|---|
| First | 1.7 | −0.0175 | 0.1738 | 153.86 | — |
| Final | 71 | −0.0195 | 0.1512 | 139.35 | −9.43% |
Relative to the first iteration, the RMS displacement falls from 153.86 to 139.35 μm, a reduction of 14.51 μm or 9.43%. The population-average objective changed little after the first eight iterations, and the ten-iteration run was treated as converged. The largest design change was the contact-trace angle, from 1.7° to 71°. The transmission-ratio derivative became slightly more negative, from −0.0175 to −0.0195, and the contact-ellipse ratio decreased from 0.1738 to 0.1512. All final values remain within the prescribed bounds.
The computed displacement-time histories show the same qualitative mesh oscillation for the two designs, with a smaller amplitude envelope for the final iteration. Within the model, the optimization therefore reduces mesh-induced vibration at the chosen operating point.
Interpretation and Limits
The 9.43% figure is a model-based reduction in RMS displacement, not a measured sound-pressure or fatigue-life improvement. It is useful as a source-excitation metric: if the mesh displacement falls, the force transmitted to the housing can fall, but the actual structural path and response determine the acoustic result.
The scope is narrow. The model uses rigid gear bodies, lumped shaft and bearing properties, a mid-face mesh, no tooth friction, and deterministic stiffness, backlash, and transmission error. The optimization is performed at one power and speed: 106.647 kW and 2097 r/min. The optimum is conditional on those assumptions and on the parameter bounds. A population of 10 and ten generations mean that the result is the best design found within the search, not a proof of global optimality. The final contact-ellipse ratio, 0.1512, is also close to the lower bound of 0.15, so the result is sensitive to that bound.
The dynamic objective excludes contact stress, edge contact, scoring, efficiency, and thermal behavior. A tooth-surface modification that lowers RMS displacement can shift the loaded contact pattern; the final design therefore needs loaded tooth contact analysis, stress checks, and experimental validation before it can be accepted for a flight-qualified reducer.
Design Implications
For helicopter main reducers, the case supports treating tooth-surface microgeometry as a source-control variable. The practical loop is: define the macro geometry, parameterize the local synthesis, pass candidate surfaces through tooth contact analysis and finite-element mesh-stiffness calculations, solve the dynamic model, and optimize the RMS response. This loop is complementary to isolation, damping, and tolerance control, which address the transmission path and manufacturing variability rather than the mesh excitation itself.
The result also shows why the three parameters should be optimized together. In this case the best individual required a large reorientation of the contact trace, from 1.7° to 71°, a modest change in transmission-error slope, and a smaller contact-ellipse ratio. A search restricted to contact-ellipse size alone would not have reached the selected design.
Within the stated model, the optimized surface lowers the computed RMS displacement by 9.43%, from 153.86 to 139.35 μm, at 106.647 kW and 2097 r/min. That is a meaningful reduction in mesh excitation for a main-reducer gear pair, but it is not by itself evidence of a quieter or longer-lived helicopter; those claims require system-level acoustic analysis and durability testing.
</m'_{pg}