Spiral bevel gears transmit motion between intersecting shafts through teeth whose trace curves obliquely across the face width. Because contact starts at one end of the tooth and rolls progressively toward the other, the load transfer is gradual rather than instantaneous, which is why these gears dominate automotive final drives, differentials and aircraft accessory drives. The same oblique geometry, however, creates a design problem with several coupled outputs: the mesh must be smooth enough to keep transmission error — and therefore vibration and noise — small, while the tooth surface must stay cool and lightly stressed enough to avoid scuffing, pitting and early fatigue. Adjusting one geometric parameter rarely improves all of these responses at once.
Where the design conflicts arise
Three parameters control most of the trade-off. The helix angle β sets how obliquely the tooth crosses the face width; a larger β lengthens the path of contact and raises longitudinal overlap εr, so more tooth pairs share the load simultaneously and the mesh stiffness variation per cycle shrinks. The price is a larger axial thrust component, which must be absorbed by the bearings and the housing. The pressure angle α governs the direction and magnitude of the normal force at the pitch point: increasing it thickens the tooth root and lowers the contact stress, but raises the radial separating force. Face width b directly enlarges the contact area available for load sharing, but it also lengthens the bending span and adds mass.
These dependencies mean the usual single-objective approach — minimizing, for example, root stress alone — tends to push the remaining responses in unfavorable directions. A multi-objective formulation that treats transmission error and longitudinal overlap as competing objectives, with contact temperature and Hertzian stress checked as outcomes, gives a more complete picture of the design space.
Parametric geometry model
Building the geometry from a small set of basic inputs requires the standard cone relations. With pitch diameter d, module m, tooth number z, outer cone distance R and pitch cone angle δ1 for the pinion:
$$R = \frac{d}{\sin\delta_1}, \qquad d = m z$$
Addendum and dedendum follow from the addendum coefficient ha*, clearance coefficient c* and profile shift x:
$$h_a = \left(h_a^{*} + x\right)m, \qquad h_f = \left(h_a^{*} + c^{*} – x\right)m$$
The corresponding addendum and dedendum angles, and the pitch cone angles of the two members of a right-angle pair, are
$$\theta_a = \arctan\frac{h_a}{R}, \qquad \theta_f = \arctan\frac{h_f}{R}$$
$$\delta_1 = \arctan\frac{z_1}{z_2}, \qquad \delta_2 = 90^{\circ} – \delta_1$$
Feeding z, m, α and β into a commercial gear design environment (Kiss soft) generates the solid model directly from these relations, so that every geometric change propagates consistently to the tooth flank and the model can be passed to the contact solver without manual reconstruction. The baseline pair used throughout the study is a strongly unequal 8 × 37 reduction, typical of a low-pinion-count automotive drive, and is summarized in Table 1.
| Parameter | Pinion | Wheel |
|---|---|---|
| Number of teeth | 8 | 37 |
| Module / mm | 8.98 | 8.98 |
| Face width / mm | 70 | 70 |
| Pressure angle / (°) | 20 | 20 |
| Helix angle / (°) | 30 | 30 |
| Shaft angle / (°) | 90 | 90 |
| Pitch cone angle / (°) | 12.20 | 77.80 |
| Pitch diameter / mm | 71.87 | 332.39 |
| Addendum / mm | 10.86 | 4.07 |
| Dedendum / mm | 5.79 | 12.58 |
| Addendum angle / (°) | 4.85 | 1.82 |
| Dedendum angle / (°) | 1.82 | 4.85 |
| Outer cone distance / mm | 170.04 | 170.04 |
| Hand of spiral | Left | Right |
Load case and baseline contact behavior
The contact analysis is run at a pinion speed of 1 000 r/min, a transmitted power of 50 kW, a torque of 477.46 N·m and a design life of 30 000 h. Under these conditions the baseline pair shows the three symptoms that motivate the optimization: a transmission error amplitude of 6.61 μm (with the error curve reaching an absolute excursion of 32.94 μm), a peak contact temperature of 221.79 °C, and a peak Hertzian contact stress of 2 676.08 MPa. The temperature and stress fields are both concentrated toward the inner end of the tooth, where the circumferential speed and the sliding velocity are lowest but the tooth is stiffest — the classic location for scoring and pitting initiation in a low-tooth-count pinion.
Multi-objective formulation
The three geometric variables most directly connected to the conflicting responses are the helix angle, the pressure angle and the face width:
$$X = \left[\beta,\ \alpha,\ b\right]^{T}$$
Transmission error e is to be suppressed, while longitudinal overlap εr is to be raised; the second objective is written as its reciprocal so that both functions are minimized:
$$\min \begin{cases} f_1(\beta,\alpha,b) = e \\ f_2(\beta,\alpha,b) = 1/\varepsilon_r \end{cases}$$
Bounds are chosen to stay within manufacturable and assembly-compatible practice for this class of gear:
$$\text{s.t.} \begin{cases} 20^{\circ} \le \beta \le 35^{\circ} \\ 15^{\circ} \le \alpha \le 25^{\circ} \\ 70\ \text{mm} \le b \le 75\ \text{mm} \end{cases}$$
NSGA-II search procedure
Because the two objectives are in tension, a Pareto-based evolutionary search is appropriate. The second-generation non-dominated sorting genetic algorithm (NSGA-II) ranks a combined parent–offspring population by non-domination, so that the current best front F1 is carried forward intact, and then uses a crowding-distance measure to preserve diversity along that front when the new population would otherwise be truncated. This combination of elitist selection and diversity preservation is what allows a well-spread approximation of the Pareto front to emerge rather than a cluster of near-identical solutions.
The search used a population of 150 individuals, a crossover probability of 0.80, a mutation probability of 0.20 and 200 generations. The resulting Pareto set traces the achievable combinations of transmission error and longitudinal overlap, and spans roughly 1.3 to 1.8 in overlap against −36 to −26 μm in transmission error, confirming that the two responses cannot be improved independently. Selecting a single design from that front requires a weighting decision; weights of 0.7 on longitudinal overlap and 0.3 on transmission error were applied, reflecting the priority given to smooth, continuous tooth engagement in this application.
Optimization results
The selected design moves to the upper bound of the helix angle, a slightly increased pressure angle and the maximum face width, as shown in Table 2. All three variables migrate to or near their constraint limits, which is expected when the objectives are monotonic in those variables over the admissible range; the limits, not the optimizer, are what ultimately cap the improvement.
| Variable | Baseline | Optimized |
|---|---|---|
| Helix angle β / (°) | 30 | 35 |
| Pressure angle α / (°) | 20 | 22 |
| Face width b / mm | 70 | 75 |
The objective functions and the two contact responses evaluated on the optimized geometry are compared with the baseline in Table 3.
| Quantity | Baseline | Optimized | Change / % |
|---|---|---|---|
| Transmission error e / μm | 6.61 | 3.51 | −46.90 |
| Longitudinal overlap εr | 1.46 | 1.78 | +21.92 |
| Peak contact temperature / °C | 221.79 | 169.25 | −23.69 |
| Peak Hertzian stress / MPa | 2 676.08 | 1 814.66 | −32.19 |
| Maximum absolute error excursion / μm | 32.94 | 27.50 | −16.51 |
The transmission error curve retains its periodic form after optimization, but its amplitude is roughly halved, and the maximum absolute excursion of the curve falls from 32.94 μm to 27.50 μm. Since transmission error is the dominant excitation source for gear whine, a reduction of this magnitude in the mesh-frequency component translates directly into lower dynamic mesh force and lower radiated noise. The peak contact temperature drops by more than 50 °C, which matters because surface distress mechanisms such as scuffing are thermally activated and highly sensitive to small temperature reductions at the contact. The Hertzian stress falls by about a third, and its distribution along the face width becomes broader and more central, indicating that a larger fraction of the tooth width participates in carrying load rather than a narrow band near the inner end.
How the geometry changes produce the improvement
The three parameter changes act through different mechanisms and reinforce one another. Raising the helix angle from 30° to 35° lengthens the contact line and increases the number of tooth pairs in simultaneous mesh, which is the direct cause of the increase in longitudinal overlap from 1.46 to 1.78; with more pairs sharing the load, the per-pair elastic deflection and hence the mesh stiffness fluctuation fall, which is why transmission error drops at the same time. Increasing the face width from 70 mm to 75 mm adds contact area, lowering the mean Hertzian pressure for the same torque. Raising the pressure angle from 20° to 22° increases the radius of curvature at the pitch point and thickens the root section, both of which reduce contact and bending stress.
The trade-off is not eliminated, only relocated. A 35° helix angle increases the axial thrust component in proportion to tan β, so the bearings and housing of the final assembly must be sized for a higher axial load than the 30° baseline; the gain in acoustic and thermal performance is paid for at the support structure. Similarly, the pressure angle increase raises the radial separating force. These consequences fall outside the two-objective formulation and would normally be handled by adding bearing life or housing stiffness as a third objective or as an explicit constraint.
Scope of the results
The numbers above are tied to one load case — 1 000 r/min, 50 kW, 477.46 N·m — and to the contact model of the design software used to generate and evaluate the geometry. Transmission error, contact temperature and contact stress are all load-dependent, and the Pareto front would shift if the torque or speed range changed, so the selected weighting and the resulting parameter set should be re-derived for a different operating point. The outcome is best read as a demonstration that a two-objective formulation coupling transmission error with longitudinal overlap is sufficient to drive simultaneous, substantial reductions in the thermal and stress responses of a spiral bevel pair, and that parametric solid modeling linked to an evolutionary search can deliver that improvement without leaving the standard design environment.