Spiral bevel gears are used where a compact intersecting-shaft drive must carry high torque smoothly. Helicopter main reducers, vehicle differentials, and marine transmissions rely on their curved flanks and localized contact pattern. The flanks are not fully conjugate: a controlled mismatch is introduced so that the contact ellipse grows to a useful size under load instead of concentrating at the tooth edge. This is the basis of quiet operation, but it also makes the mesh sensitive to the relative position of the two axes. Manufacturing scatter, assembly tolerance, thermal growth, and elastic support deflection all appear as small changes in shaft position and angle. The contact pattern can migrate toward the toe or heel, the transmission-error curve can become asymmetric, and edge contact, noise, and fatigue can follow.
The design problem therefore has two objectives that are not independent. The mesh should retain a long effective contact path, which means a high contact ratio, and it should change little when the support moves, which means low sensitivity. These objectives often conflict. Shortening the contact path can make the pattern more stable, but it reduces the number of teeth sharing load. A practical robust design must quantify both effects and expose the trade-off rather than optimizing one and checking the other afterward.
Second-order contact parameters as design levers
In the Gleason local-conjugate method, the gear is generated first and the pinion machine settings are derived from a reference point on the gear flank. The designer prescribes three second-order quantities at that point: the tangent direction of the contact path, the first derivative of the transmission-ratio function, and the semi-major axis of the instantaneous contact ellipse. These are denoted here as \(x_1\), \(x_2\), and \(x_3\). They control the length, orientation, and curvature of the contact pattern. Changing them does not alter the macro geometry of the gear pair, but it changes how the pattern responds to a small displacement of the axes.
This is why the second-order parameters are natural optimization variables. They are set on the machine, they have a direct geometric meaning, and they can be adjusted without redesigning the blank. The difficulty is that their effect on contact ratio and on misalignment sensitivity is nonlinear. A parameter set that produces a long contact path may also produce a pattern that moves strongly when the pinion or gear shifts axially.
Equivalent misalignment in tooth contact analysis
To include support movement in the contact analysis, four equivalent misalignments are defined in a common assembly frame: pinion axial offset \(\Delta A_1\), gear axial offset \(\Delta A_2\), shaft-center distance offset \(\Delta E\), and shaft-angle offset \(\Delta\gamma\). The shaft-angle error is converted to an equivalent linear displacement at the mid-face cone distance so that all four quantities can be compared in length units. This is the basis of the error tooth contact analysis (ETCA) used here.
At each pinion roll angle, the two flanks are transformed into the assembly frame. Contact requires equal position vectors and equal unit normals:
$$
\mathbf{r}_h^{(1)}=\mathbf{r}_h^{(2)},\qquad
\mathbf{n}_h^{(1)}=\mathbf{n}_h^{(2)}.
$$
The position and normal vectors depend on the surface coordinates of the two flanks and on the two roll angles. With the pinion roll angle prescribed, the remaining unknowns are solved from the contact equations and the two meshing equations. Sweeping the pinion roll angle gives the contact path. The instantaneous contact ellipses along that path form the tooth contact pattern.
The transmission error is defined from the deviation between the actual gear roll angle and the roll angle implied by the nominal ratio:
$$
\delta(\phi_1)=
\left(\phi_2-\phi_2^{(0)}\right)
-\frac{N_1}{N_2}
\left(\phi_1-\phi_1^{(0)}\right),
$$
where \(\phi_1\) and \(\phi_2\) are the current pinion and gear roll angles, \(\phi_1^{(0)}\) and \(\phi_2^{(0)}\) are their values at the reference contact point, and \(N_1\) and \(N_2\) are the tooth numbers. A single branch of the curve is repeated with the pitch period to form the complete transmission-error diagram.
Quantifying contact-pattern sensitivity
A contact pattern is not a single number. It has area, position, and orientation. The pattern area \(S\) is approximated by connecting the major-axis endpoints of successive contact ellipses into a polygon and summing the component triangles. The pattern position \((X_c,Y_c)\) is taken as the mean of the polygon vertices. The pattern direction \(\alpha\) is the angle between the line joining the first and last contact centers and the gear blank axis.
These four metrics can be differentiated with respect to the four equivalent misalignments. The resulting sensitivity matrix has the form
$$
\mathbf{D}=
\begin{bmatrix}
\dfrac{\partial X_c}{\partial \Delta A_1} &
\dfrac{\partial X_c}{\partial \Delta A_2} &
\dfrac{\partial X_c}{\partial \Delta E} &
\dfrac{\partial X_c}{\partial \Delta\gamma}\\[6pt]
\dfrac{\partial Y_c}{\partial \Delta A_1} &
\dfrac{\partial Y_c}{\partial \Delta A_2} &
\dfrac{\partial Y_c}{\partial \Delta E} &
\dfrac{\partial Y_c}{\partial \Delta\gamma}\\[6pt]
\dfrac{\partial S}{\partial \Delta A_1} &
\dfrac{\partial S}{\partial \Delta A_2} &
\dfrac{\partial S}{\partial \Delta E} &
\dfrac{\partial S}{\partial \Delta\gamma}\\[6pt]
\dfrac{\partial \alpha}{\partial \Delta A_1} &
\dfrac{\partial \alpha}{\partial \Delta A_2} &
\dfrac{\partial \alpha}{\partial \Delta E} &
\dfrac{\partial \alpha}{\partial \Delta\gamma}
\end{bmatrix}.
$$
The matrix shows which misalignment dominates which pattern feature. For optimization, a scalar index is used:
$$
D=\sum_{i,j}\left|D_{ij}\right|.
$$
This index is an aggregate measure of pattern change, not a physical displacement. A lower \(D\) means that the pattern stays closer to its aligned position, size, and orientation when the support moves. It is useful for ranking candidate designs, but it should not be read as a tolerance in millimetres.
Contact ratio from the transmission-error curve
The contact ratio \(C\) is the ratio of the mesh duration of one tooth pair to the pitch duration:
$$
C=\frac{T_C}{T_Z},\qquad T_Z=\frac{360^\circ}{N_1},
$$
where \(T_C\) is the angle through which the gear rotates from tooth entry to exit and \(T_Z\) is the pinion single-tooth pitch angle. The design contact ratio is the upper bound reached when the tooth tip just participates without edge contact. The actual contact ratio under load is different. At near-zero load, the actual contact ratio of a locally conjugate spiral bevel gear tends to 1.0 because only the central part of the path carries load. As load increases, elastic deformation closes the gap between adjacent teeth and the actual contact ratio rises toward the design value. Edge contact can increase it slightly further, but that condition should be avoided because it concentrates stress and excites vibration.
For a locally conjugate flank, the root-sum-square formula based on transverse and overlap contact ratios is reliable mainly at full load. A more general method is to read the design contact ratio directly from the transmission-error curve. The width \(L’\) of one transmission-error branch is compared with the pitch \(L\) between adjacent branches:
$$
C=\frac{L’}{L}.
$$
This definition follows from the same contact-analysis model that produces the pattern, so it can be applied to loaded and unloaded cases. It also shows the geometric route to a higher contact ratio: reducing the angle between the contact path and the root cone lengthens the path and increases \(C\). Representative values from the same design family are given below.
| Angle between contact path and root cone (°) | Design contact ratio |
|---|---|
| 71.3 | 1.5 |
| 43.5 | 1.9 |
| 25.1 | 2.3 |
Bi-objective optimization formulation
The optimization seeks a high contact ratio and a low sensitivity index by adjusting the second-order contact parameters. The design vector is \(\mathbf{x}=(x_1,x_2,x_3)\). The problem can be written as
$$
\max_{\mathbf{x}} C(\mathbf{x}),\qquad
\min_{\mathbf{x}} D(\mathbf{x}),
$$
subject to the pattern metrics remaining within an allowable range \(U\) and the second-order parameters remaining within their design range \(S\). The two objectives are normalized and combined with variable weights:
$$
\min_{\mathbf{x}} f(\mathbf{x})
=\omega_1 f_1(\mathbf{x})+\omega_2 f_2(\mathbf{x}),
\qquad \omega_1+\omega_2=1,
$$
where \(f_1\) is the normalized inverse contact ratio and \(f_2\) is the normalized sensitivity. A fixed weighted sum produces one solution that depends strongly on the chosen weights. To trace neighboring non-dominated solutions, the weights are changed in steps:
$$
\omega_1^{j+1}=\omega_1^{j}+\Delta,\qquad
\omega_2^{j+1}=\omega_2^{j}-\Delta,\qquad
\Delta=\frac{1}{n}.
$$
For each parameter combination, ETCA gives the contact pattern and the transmission-error curve, from which \(C\) and \(D\) are computed. Response functions are fitted to these samples, and the weighted optimization is performed on the fitted functions. The result is a Pareto trajectory: a set of second-order parameter choices that trade contact ratio against sensitivity in a controlled way.
Example gear pair and misalignment response
The procedure is illustrated on a 90° spiral bevel gear pair with 38 pinion teeth and 51 gear teeth. The main data are listed below.
| Parameter | Pinion | Gear |
|---|---|---|
| Shaft angle (°) | 90 | 90 |
| Pressure angle (°) | 20 | 20 |
| Number of teeth | 38 | 51 |
| Module (mm) | 3.5 | 3.5 |
| Spiral angle (°) | 35 | 35 |
| Face width (mm) | 22 | 22 |
| Radial modification coefficient | +0.1737 | −0.1737 |
| Hand | Right | Left |
| Working flank | Concave | Convex |
| Normal backlash (mm) | 0.29–0.40 | 0.29–0.40 |
The four equivalent misalignments were applied one at a time while the other three were set to zero. The transmission-error curves and gear contact patterns respond differently in each case.
| Offset and test range | Transmission-error response | Contact-pattern response |
|---|---|---|
| Pinion axial \(\Delta A_1\): TE at ±0.05 mm; pattern at ±0.3 mm | Curve becomes asymmetric; peak transmission error rises. | Pattern shifts toward the toe (small end); contact-trace angle to the root cone decreases. |
| Gear axial \(\Delta A_2\): ±0.05 mm | As the offset goes positive, the left branch moves away from the parabola vertex and the right branch moves toward it; peak transmission error rises. | Pattern shifts toward the heel (large end); contact-trace angle to the root cone increases. |
| Shaft-center distance \(\Delta E\): ±0.05 mm | A positive offset makes the curve more symmetric and can reduce the peak transmission error relative to the aligned case. | Pattern shifts toward the heel; contact-ellipse major axis shortens; angle to the root cone increases; pattern area decreases. |
| Shaft angle \(\Delta\gamma\): ±0.002 rad | Amplitude increases, especially at \(-0.002\) rad. | Pattern shifts toward the toe; contact-trace angle to the root cone decreases. |
These trends are not additive in a simple way. A positive pinion axial offset moves the pattern toward the toe, while a positive gear axial offset moves it toward the heel. The shaft-center distance and shaft-angle errors can reinforce or cancel those shifts. The actual pattern centroid, area, and orientation are therefore the combined result of all four offsets. A design that is robust to one offset may still be weak to another, so the sensitivity matrix must include all four rather than using a single axial displacement as a proxy.
Optimization results and the contact-ratio–sensitivity trade-off
The optimization used a representative combined offset of \(\Delta A_1=-0.12\) mm, \(\Delta A_2=0.05\) mm, \(\Delta E=0.1\) mm, and \(\Delta\gamma=-0.00025\) rad. Three weight settings produced the local optima below.
| Contact-ratio weight \(\omega_1\) | \(x_1\) (°) | \(x_2\) | \(x_3\) | Contact ratio \(C\) | Sensitivity \(D\) |
|---|---|---|---|---|---|
| 0.2 | 72.763 | 0 | 0.15 | 1.458 | 32.127 |
| 0.5 | 70.126 | 0.008 | 0.15 | 1.750 | 33.948 |
| 0.8 | 70.209 | 0.008 | 0.15 | 1.875 | 131.044 |
The table makes the conflict explicit. At \(\omega_1=0.2\), the sensitivity is slightly lower than at the balanced setting, but the contact ratio falls to 1.458. At \(\omega_1=0.8\), the contact ratio rises to 1.875, but the sensitivity jumps to 131.044. The balanced setting at \(\omega_1=0.5\) keeps the contact ratio at 1.750 while holding the sensitivity to 33.948.
The initial design had a contact ratio of 1.750 and a sensitivity index of 70.740. After optimization at the balanced weight, the contact ratio is unchanged at 1.750, but the sensitivity index falls to 33.948, a reduction of 52.0%. Under the same prescribed misalignment, the optimized pattern shifts less, and the computed contact path is less curved. Both features favor smoother meshing and a more stable pattern.
| Design | Contact ratio \(C\) | Sensitivity \(D\) | Change in \(D\) |
|---|---|---|---|
| Initial | 1.750 | 70.740 | — |
| Optimized, \(\omega_1=0.5\) | 1.750 | 33.948 | −52.0% |
Experimental check with contact-pattern coloring
The predicted pattern shifts were checked with coloring tests on a high-speed spiral bevel gear rig at 600 r/min without external load. A 5:1 mixture of red lead powder and machine oil was applied to the pinion. With a pinion axial offset of \(\Delta A_1=+0.25\) mm, the pattern moved toward the toe. With a gear axial offset of \(\Delta A_2=+0.25\) mm, it moved toward the heel. These directions agree with the ETCA predictions.
The mixture has a low viscosity and is thrown outward by centrifugal action, so the painted region is larger than the true contact area. The measurement therefore used the largest continuous visible painted area for area, the angle of the line joining the endpoints of the area midline for direction, and the midpoint of that midline for the centroid. The comparison used a combined offset of \(\Delta A_1=0.25\) mm, \(\Delta A_2=0.25\) mm, \(\Delta E=0\), and \(\Delta\gamma=0\). The changes are relative to the aligned case.
| Change relative to aligned case | Before optimization | After optimization |
|---|---|---|
| Pattern area change (%) | 12.0 | 8.9 |
| Pattern direction change (°) | 1.7 | 0.8 |
| Centroid \(x\) change (cm) | 0.60 | 0.50 |
| Centroid \(y\) change (cm) | 0.12 | 0.05 |
All four change indicators decreased after optimization. The area change fell from 12.0% to 8.9%, the direction change from 1.7° to 0.8°, the horizontal centroid change from 0.60 cm to 0.50 cm, and the vertical centroid change from 0.12 cm to 0.05 cm. This supports the predicted improvement in geometric robustness.
Engineering implications and scope
For helicopter main reducers and accessory drives, the design question is not simply “minimize sensitivity” or “maximize contact ratio.” The two objectives compete. A weighted Pareto approach gives the designer a set of admissible second-order parameter choices, and the final choice can reflect the expected support stiffness, thermal range, and allowable pattern migration. Because the four equivalent offsets enter the same sensitivity matrix, the method can also be used to compare bearing arrangements or tolerance allocations before hardware is cut.
The reported sensitivity is a geometric index derived from contact-pattern metrics. It does not by itself quantify dynamic excitation, lubricant film behaviour, or tooth-root fatigue. The optimization also depends on the prescribed misalignment set; if the actual offsets are random, the same pattern metrics can be embedded in a probabilistic or six-sigma framework. The coloring test is a no-load geometric check, not a loaded endurance test. Within these limits, the method provides a clear way to design a spiral bevel gear pair that keeps a high contact ratio while reducing the pattern shift caused by support misalignment.