Spiral bevel gears transmit rotation between intersecting shafts through curved teeth that engage progressively across the face width. This gradual engagement supports high torque density, high ratios and comparatively low vibration and noise, which is why the topology is widely used in automotive final drives, aviation accessory drives and industrial reducers. The same curved geometry also creates a tightly coupled design problem. Transmission error, longitudinal overlap, contact temperature and Hertzian stress all respond to the same small set of macro-geometry variables, and a change that improves one metric can degrade another.
The optimization workflow examined here treats the spiral angle, pressure angle and face width as design variables, uses transmission error and longitudinal overlap as competing objectives, and applies the second-generation non-dominated sorting genetic algorithm (NSGA-II) to generate a Pareto set. A parametric model built in KISSsoft supplies the gear geometry and the contact calculation. After a compromise design is selected from the Pareto set, contact temperature and Hertzian stress are recalculated on the optimized tooth form. The resulting design raises longitudinal overlap from 1.46 to 1.78 while reducing transmission error from 6.61 to 3.51 μm; the same contact analysis gives lower peak temperature and stress.
Geometry and Parametric Model
Spiral bevel gear geometry is defined at the outer cone. For a pinion with \(z_1\) teeth meshing with a wheel with \(z_2\) teeth on a 90° shaft angle, the pitch cone angles are
$$ \delta_1 = \arctan\left(\frac{z_1}{z_2}\right), \qquad \delta_2 = 90^\circ – \delta_1 $$
The pitch diameter \(d\), outer cone distance \(R\), addendum \(h_a\) and dedendum \(h_f\) follow from the module \(m\), the addendum modification \(x\), the addendum coefficient \(h_a^*\) and the clearance coefficient \(c^*\):
$$ d = mz, \qquad R = \frac{d}{2\sin\delta_1} $$
$$ h_a = (h_a^* + x)m, \qquad h_f = (h_a^* + c^* – x)m $$
These relations show why the parameter set is compact but not independent. The cone angles are fixed by the tooth-count ratio, the pitch diameter follows from module and tooth count, and the outer cone distance follows from pitch diameter and cone angle. Once the coefficients and modification are chosen, the tooth flank, root fillet and contact pattern are generated as a coupled geometry rather than as independent CAD features.
| Parameter | Pinion | Wheel |
|---|---|---|
| Number of teeth | 8 | 37 |
| Module (mm) | 8.98 | 8.98 |
| Face width (mm) | 70 | 70 |
| Pressure angle (°) | 20 | 20 |
| Spiral 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 | Left | Right |
The model is built in the bevel gear module of KISSsoft. The inputs are the tooth counts, module, pressure angle, spiral angle, face width, shaft angle and hand of rotation. The software reconstructs the tooth flanks and the assembly from those parameters and exports the solid model for contact analysis. Because the flank is regenerated from the parameter set, a change in spiral angle or pressure angle propagates directly to the tooth surface and the contact pattern; there is no need to rebuild the solid from scratch or to edit individual flank patches. This is the practical advantage of a dedicated gear modeler over a general CAD kernel, which can hold the solid but does not by itself enforce gear-specific geometry.
The pinion has 8 teeth and the wheel 37, giving a ratio of 4.625:1. The shaft angle is 90°, and the initial pressure angle, spiral angle and face width are 20°, 30° and 70 mm, respectively. These values define the baseline design that is used for the contact calculation and the subsequent optimization.
Baseline Contact Analysis
The contact calculation is performed at 1 000 r/min and 50 kW, which corresponds to 477.46 N·m. A design life of 30 000 h is entered as an input condition; it is not a fatigue life predicted by the contact model. For the baseline geometry, the calculated transmission error is 6.61 μm, with a maximum peak-to-peak fluctuation of 32.94 μm. The maximum contact temperature is 221.79 °C and the maximum Hertzian stress is 2 676.08 MPa.
Transmission error is the difference between the rotation that an ideal conjugate pair would produce and the rotation delivered under load. It is not eliminated by good design; it is shaped. The peak-to-peak amplitude is the relevant excitation measure because it governs the dynamic force transmitted to the housing and bearings. A fluctuation of 32.94 μm at this load indicates a mesh-stiffness variation that can excite the driveline. The 221.79 °C contact temperature and 2 676.08 MPa peak Hertzian stress are the thermal and fatigue-side consequences of the same contact pattern.
Multi-Objective Formulation
Three variables are selected because they control the contact pattern and the variation of mesh stiffness through the cycle. The spiral angle \(\beta\) fixes the inclination of the tooth trace and therefore the face contact ratio. The pressure angle \(\alpha\) fixes the distribution of normal load across the profile and the shape of the contact ellipse. The face width \(b\) fixes the loaded width available for load sharing. The optimization vector is
$$ \mathbf{X} = [\beta,\alpha,b]^T $$
The two objectives are the transmission error \(e\) and the longitudinal overlap \(\varepsilon_r\). Because overlap is to be maximized, its reciprocal is minimized:
$$ \min f_1(\mathbf{X}) = e(\mathbf{X}), \qquad \min f_2(\mathbf{X}) = \frac{1}{\varepsilon_r(\mathbf{X})} $$
The search is restricted to the following box:
$$ 20^\circ \le \beta \le 35^\circ, \qquad 15^\circ \le \alpha \le 25^\circ, \qquad 70\,\text{mm} \le b \le 75\,\text{mm} $$
| Variable | Symbol | Initial value | Lower bound | Upper bound |
|---|---|---|---|---|
| Spiral angle | \(\beta\) | 30° | 20° | 35° |
| Pressure angle | \(\alpha\) | 20° | 15° | 25° |
| Face width | \(b\) | 70 mm | 70 mm | 75 mm |
NSGA-II is used because the two objectives conflict. Increasing longitudinal overlap generally requires a larger spiral angle or face width, which can raise contact stress or temperature, while reducing transmission error may call for a different pressure angle. A weighted single-objective search would collapse this trade-off before the designer sees it. NSGA-II instead ranks solutions by Pareto dominance, uses crowding distance to maintain spread along the front, and combines parent and offspring populations through elitist selection. The settings are a population of 150, crossover probability 0.80, mutation probability 0.20 and 200 generations. The search returns a set of non-dominated designs rather than a single optimum. A weighted compromise is then selected from the Pareto set with a weight of 0.7 on longitudinal overlap and 0.3 on transmission error.
Optimization Results
The selected design uses a 35° spiral angle, a 22° pressure angle and a 75 mm face width. Both the spiral angle and the face width lie on their upper bounds, while the pressure angle settles inside its range.
| Variable | Before | After |
|---|---|---|
| Spiral angle | 30° | 35° |
| Pressure angle | 20° | 22° |
| Face width | 70 mm | 75 mm |
| Metric | Before | After | Change |
|---|---|---|---|
| Transmission error (μm) | 6.61 | 3.51 | −46.90% |
| Longitudinal overlap | 1.46 | 1.78 | +21.92% |
| Maximum contact temperature (°C) | 221.79 | 169.25 | −23.69% |
| Maximum Hertzian stress (MPa) | 2 676.08 | 1 814.66 | −32.19% |
At the selected point, the longitudinal overlap is 1.78, a 21.92% increase over 1.46. The transmission error falls to 3.51 μm, a 46.90% reduction, and the maximum fluctuation falls from 32.94 to 27.50 μm. Re-running the contact analysis on the optimized geometry gives a maximum contact temperature of 169.25 °C and a maximum Hertzian stress of 1 814.66 MPa, reductions of 23.69% and 32.19%, respectively.
Two features of the solution are worth noting. First, the selected design is a boundary solution rather than an interior optimum. Within the stated box, the optimizer takes both the spiral angle and the face width to their upper limits; the pressure angle settles at 22°. This suggests that the stated ranges, rather than a local balance among the three variables, constrain the achievable overlap. Expanding the spiral angle or face width could raise overlap further, but would also increase axial thrust, mass and packaging envelope, so the bound is an engineering choice as much as a numerical one. Second, temperature and stress were not objective functions. Their reduction is a post-optimization contact result, not something NSGA-II directly minimized. It is consistent with a more evenly loaded contact pattern, but it should not be read as a guaranteed coupled improvement at every operating point.
Interpretation and Engineering Limits
The main mechanism is the shift in contact ratio and load sharing. A 35° spiral angle and a 75 mm face width lengthen the contact line and increase the number of tooth pairs carrying load at any instant. The result is a smaller variation in mesh stiffness through the cycle, which appears as a lower transmission-error amplitude. A 22° pressure angle modifies the contact ellipse and the normal-force distribution, contributing to the lower peak Hertzian stress. Lower peak stress and more gradual load transfer are consistent with the lower contact temperature, although the contact calculation does not separate sliding friction, lubricant-film behaviour and bulk tooth temperature.
The engineering appeal of the result is straightforward: a smoother mesh, lower vibration and noise excitation, lower thermal load at the contact and higher load capacity for the same envelope. The limits are equally important. The analysis represents one speed and one load case, with nominal mounting and a single gear pair. It does not by itself cover housing compliance, shaft misalignment, manufacturing scatter, lubricant condition or transient thermal effects. The 30 000 h figure is a design input, not a calculated service life. The weighted selection of 0.7 and 0.3 is a designer preference; a different weighting would select a different point on the same Pareto front.
Qualification of the design would require measured transmission error under load, contact-pattern inspection, temperature measurement and strain or stress verification on a gear test rig. A load and speed sweep would show whether the gains persist away from the nominal operating point, and a formulation that treats temperature and stress as explicit objectives or constraints would test whether the post-optimization reductions can be pushed further without sacrificing overlap.
Conclusion
The workflow links a gear-specific parametric model to a multi-objective evolutionary search. For the 8/37 tooth pair, NSGA-II identifies a compromise with a 35° spiral angle, a 22° pressure angle and a 75 mm face width. Longitudinal overlap increases from 1.46 to 1.78, a 21.92% gain, while transmission error decreases from 6.61 to 3.51 μm, a 46.90% reduction. The optimized contact analysis gives a maximum contact temperature of 169.25 °C and a maximum Hertzian stress of 1 814.66 MPa, reductions of 23.69% and 32.19% relative to the baseline. Within the studied load case, the design moves the mesh toward more gradual load transfer and lower thermal and stress peaks. Because the result is a single-point simulation with the optimum on two variable bounds, the next steps are experimental verification and extension to load and speed spectra, misalignment, and a formulation that treats temperature and stress as explicit objectives or constraints.