Why the Contact Pattern Is a Design Variable
Spiral bevel gears transmit rotation between intersecting axes and are among the most heavily loaded components in automotive axles, aircraft engines, helicopter transmissions, and ship drives. Their teeth are cut on a face-milling or face-hobbing machine with a rotating cutter head, so the flank geometry is not specified directly by a drawing: it emerges from the interplay of blank dimensions, cutter geometry, and machine settings such as radial and angular cutter position, workhead tilt, vertical offset, axial setting, sliding base, and roll ratio.
Because of that indirect relationship, the meshing quality of a spiral bevel gear pair is judged by two computed or measured quantities rather than by the flank coordinates themselves. The first is the contact pattern, the trace left on the tooth surface when the pair is rolled under light load; the second is the transmission error, the deviation between the actual angular position of the driven gear and its theoretical value. Both are highly sensitive to the settings used for the pinion, and both can be predicted before any metal is cut.
The classical Gleason approach to this problem is a local-conjugate procedure. The wheel (the larger member) is cut first and its flank is treated as given; a reference point is chosen on that flank, and the pinion surface is derived so that the two surfaces are conjugate at that point. The difficulty is that the second-order properties of the conjugate surface at the reference point, meaning its principal curvatures and the orientation of the contact ellipse, are not known in advance. The pinion curvature must therefore be corrected by trial and error, which is slow and offers no guarantee of reaching a favorable contact pattern.
Two developments changed the situation. In 1961, M. L. Baxter of the Gleason Works published the formulation of tooth contact analysis (TCA), which made it possible to compute the contact trace and the transmission error of a gear pair numerically, without cutting it. In 1981, Krenzer extended the idea to loaded conditions through loaded tooth contact analysis (LTCA). The idea of prescribing the second-order contact conditions rather than discovering them afterwards is associated with F. L. Litvin and co-workers, whose local synthesis method forms the basis of the design procedure discussed here. For orientation, the spiral bevel concept itself is much older; T. Olivier described an enveloping-surface construction for skew-axis gearing in 1820, and the first Gleason machine for cutting teeth with a circular arc longitudinal form was built in 1913.
Kinematics Behind Local Synthesis
Consider two tooth surfaces, $\Sigma_1$ on the pinion and $\Sigma_2$ on the wheel, in continuous tangency. At a common contact point they share a common normal, and the velocity components along that normal must be equal. Writing the relative velocity of the two surfaces as $\mathbf{v}^{(12)}$ and the unit normal as $\mathbf{n}^{(2)}$, the condition of continuous tangency is the equation of meshing
$$\mathbf{n}^{(2)}\cdot\mathbf{v}^{(12)}=0 .$$
Differentiating along the path of contact gives the second fundamental relation
$$\dot{\mathbf{n}}^{(2)}\cdot\mathbf{v}^{(12)}+\mathbf{n}^{(2)}\cdot\dot{\mathbf{v}}^{(12)}=0 ,$$
which, together with Rodrigues’ formula $\dot{\mathbf{n}}^{(i)}=-K^{(i)}\mathbf{v}_r^{(i)}$ for each surface, converts a purely kinematic statement into a statement about curvature. When both surfaces are projected onto their common tangent plane, the relation can be written as a linear system in the two components of the relative velocity along the principal directions of one surface. Two cases follow from that system.
In the first case, the surfaces are in line contact. This is the situation during generation, when the cutter head surface and the workpiece surface are conjugate along a line. The relative velocity direction is undetermined, the linear system is singular in a weaker sense, and its augmented matrix has rank less than two. Three algebraic conditions then yield the principal curvatures and principal directions of the generated surface directly from those of the cutter surface, with no free parameter.
In the second case, the two finished members are in point contact. Here the relative velocity direction is determined, the coefficient matrix has rank exactly two, and the vanishing of its determinant produces a single scalar relation among the curvatures of the two surfaces, their relative orientation, and the derivative of the transmission ratio. That single relation contains one more unknown than the kinematics supply, and the gap is filled by prescribing three parameters at the reference point.
Three Presets That Shape the Contact
Local synthesis is therefore the act of choosing, at a reference point $M$ on the wheel flank, three second-order quantities that define the desired contact behavior in the neighborhood of that point.
Slope of the transmission ratio function, $m’_{21}$. For a truly conjugate pair the ratio $N_1/N_2$ is constant and the transmission error vanishes identically. Real spiral bevel gears are deliberately non-conjugate, and the deviation is written as
$$\delta\varphi_1=\left(\varphi_2-\varphi_2^{\,0}\right)-\frac{N_1}{N_2}\left(\varphi_1-\varphi_1^{\,0}\right),$$
which, expanded to second order about the reference position, becomes
$$\delta\varphi_1=\tfrac{1}{2}\,m’_{21}\left(\varphi_1-\varphi_1^{\,0}\right)^2 .$$
The preset $m’_{21}$ therefore sets the amplitude of a parabolic error curve. A negative value produces a crown-shaped curve in which the wheel lags its theoretical position in the middle of the mesh and leads at the ends, a form that spreads the load away from the tooth edges. Practical ranges used in design are roughly $-0.008$ to $0.008$, with negative values on the working flank and positive values on the coast flank.
Direction of the contact path, $\eta_2$. This angle, measured in the tangent plane of the wheel from a reference principal direction, fixes the orientation of the contact trace and hence how the contact migrates across the face width as the gears rotate. For a left-hand pinion meshing with a right-hand wheel, the working flank is assigned values between $0^\circ$ and $90^\circ$ and the coast flank between $90^\circ$ and $180^\circ$. Once $\eta_2$ is fixed, the corresponding direction $\eta_1$ on the pinion and the components of relative velocity follow from the basic equations.
Semi-major axis of the instantaneous contact ellipse, $a$. Under load the theoretical point contact spreads into an ellipse whose axes are governed by the combined curvature of the two surfaces and by the elastic approach $\varepsilon$, conventionally taken as $0.00635\ \mathrm{mm}$:
$$a=\sqrt{\frac{\varepsilon}{A}},\qquad b=\sqrt{\frac{\varepsilon}{B}},$$
where $A$ and $B$ are built from the curvature sums $K_\Sigma=k_1+k_2$, the curvature differences $g=k_1-k_2$ of the two surfaces, and the angle $\sigma^{(12)}$ between their principal directions. Prescribing $a$ between $0.15$ and $0.20$ of the face width effectively fixes the contact ratio and the nominal size of the loaded footprint; the pinion curvatures and the orientation angle then follow in closed form. Larger values give a broader, more forgiving pattern at the cost of a longer contact path and greater sensitivity to misalignment.
From Presets to Machine Settings
The computational chain that turns these presets into settings an operator can dial in runs as follows. Blank data, namely tooth counts, module, face width, pressure angle, pitch, face and root cone angles, spiral angle, addendum and whole depth, are established first, typically following the Gleason calculation card. Wheel settings are then obtained from the blank: radial and angular cutter positions from the mean cone distance and mean spiral angle, sliding base from the root apex distance, and roll ratio from the pitch cone angle.
The wheel flank is described analytically by transforming the cutter-head cone from the cutter frame through the cradle frame and the workhead frame into the wheel frame, and then eliminating the surface parameter with the equation of meshing. The result is a two-parameter surface $\mathbf{r}_2(\theta_g,\phi_g)$ with an accompanying normal $\mathbf{n}_2(\theta_g,\phi_g)$.
The reference point $M$ is located on that surface by intersecting it with two conditions: a prescribed position along the face width, expressed through an axial coordinate $XL$, and a prescribed position along the tooth height, expressed through a radial coordinate $RL$. The resulting nonlinear system is solved iteratively, with initial guesses of $\theta_g^{(0)}=270^\circ+\beta$ and $\phi_g^{(0)}=0^\circ$ for the left-hand pinion and right-hand wheel combination. Substituting the solution back gives the position vector and normal at $M$, from which the principal curvatures and directions of the wheel flank follow by applying the line-contact relations to the cutter surface.
With the wheel characterized at $M$, the two surfaces are brought into mesh in a fixed frame, with the wheel axis inclined at the shaft angle. The meshing condition at $M$, evaluated with the theoretical ratio, determines the wheel rotation angle at which contact occurs. The three presets are then applied through the point-contact relations to yield the principal curvatures $k_I^{(1)},k_{II}^{(1)}$ and directions of the pinion at the same point.
The final step reverses the direction of the calculation. The pinion is generated by a cutter, and cutter and workpiece are in line contact, so the line-contact conditions now run backwards: the known pinion curvatures determine the cutter tip radius, and the known relative velocity determines the roll ratio. The remaining settings follow from the geometry of the pinion generation frame, including the vertical offset, the axial setting, the sliding base, and the radial and angular cutter positions.
Case: A 9 × 39 Axle Gear Set
A representative design case illustrates the procedure. The gear pair has 9 pinion teeth and 39 wheel teeth, a module of 5.69 mm, a face width of 35 mm, and a 22.5° pressure angle, giving a ratio of 4.33 and a geometry typical of a light-vehicle drive axle. The wheel is cut by the generating method with a double-sided cutter, and the pinion by the modified-roll method, with separate settings for the concave and convex flanks.
| Parameter | Pinion | Wheel |
|---|---|---|
| Number of teeth | 9 | 39 |
| Module (mm) | 5.69 | |
| Face width (mm) | 35 | |
| Pressure angle (°) | 22.5 | |
| Pitch cone angle (°) | 12.99 | 77.01 |
| Face cone angle (°) | 16.57 | 79.72 |
| Root cone angle (°) | 10.28 | 73.43 |
| Mean spiral angle (°) | 29.43 | |
| Whole depth (mm) | 11.38 | |
| Addendum (mm) | 5.98 | 4.27 |
| Wheel machining parameter | Value | |
|---|---|---|
| Cutter diameter (mm) | 190.50 | |
| Cutter top distance (mm) | 1.50 | |
| Profile angle, convex flank (°) | 22 | |
| Profile angle, concave flank (°) | −18 | |
| Blank mounting angle (°) | 73.43 | |
| Sliding base (mm) | 0 | |
| Radial cutter position (mm) | 96.64 | |
| Angular cutter position (°) | 59.21 | |
| Pinion machining parameter | Concave flank | Convex flank |
|---|---|---|
| Cutter tip diameter (mm) | 188.46 | 187.90 |
| Profile angle (°) | 18 | −21 |
| Blank mounting angle (°) | 10.28 | 10.28 |
| Axial setting (mm) | −7.89 | 11.57 |
| Sliding base (mm) | 1.41 | −2.28 |
| Vertical offset (mm) | −9.44 | 2.26 |
| Radial cutter position (mm) | 83.80 | 142.41 |
| Angular cutter position (°) | 60.56 | 68.23 |
| Roll ratio | 5.46 | 3.28 |
The asymmetry between the two pinion flanks is a direct consequence of the local synthesis formulation. Each flank is treated as a separate surface with its own preset set, so the cutter tip radius, the radial cutter position, and the roll ratio differ substantially between the working and coast sides. The wheel, by contrast, is cut with one double-sided cutter, and its two flanks share a common set of machine settings.
Tooth Contact Analysis as a Feedback Filter
TCA converts the computed settings into an observable contact pattern. In a fixed frame attached to the housing, the two surfaces must have a common point and a common normal at the instantaneous contact, giving six scalar equations in eight unknowns: the two surface parameters and the rotation angle of each member. Fixing the pinion rotation angle reduces the system to five equations in five unknowns, which is solved numerically and repeated for a sequence of pinion angles to trace the path of contact.
Applying this to the 9 × 39 pair gave a contact trace that is nearly straight and slightly inclined across the face width, with contact ellipses of comparable size at the entry, middle, and exit of the mesh. The transmission error curve is approximately symmetric about the mid-mesh position, continuous through the handover from one tooth pair to the next, and free of the sharp upturn that signals edge contact. This last feature matters most: when the error curve turns upward at the end of the mesh, the contact has run off the tooth end, and the resulting edge loading is a common precursor of tooth breakage.
The presets are not final after one pass. If the computed trace is too long, too short, or too close to the toe or heel, $a$, $\eta_2$, and $m’_{21}$ are adjusted and the pinion settings recomputed. This loop, in which TCA acts as the acceptance criterion for a candidate setting, is what replaces the trial-and-error curvature correction of the purely local-conjugate method. It reduces the number of physical cutting trials, but it does not eliminate iteration, and the choice of presets remains a matter of engineering judgment informed by experience with similar gear sets.
Solid Modeling and Virtual Assembly
A geometry that exists only as a pair of two-parameter surface equations is of limited use for downstream analysis, so the flanks are converted into solid models through a rotating projection scheme. Each tooth surface is projected onto its axial section, where the available region is bounded by four straight segments corresponding to the face cone, root cone, front cone, and back cone. The four corner points follow from the cone angles, the outer cone distance, the addendum, and the whole depth; a grid of, for example, 5 points across the tooth height and 9 along the face width is then generated by interpolating the boundary lines and intersecting the two families of lines.
Each planar grid node $(X,Y)$ corresponds to a unique point on the three-dimensional flank through the relations
$$x=X,\qquad y^2+z^2=Y^2,$$
which, for the wheel, form a nonlinear system in the two surface parameters $\theta_g$ and $\phi_g$. For the pinion the equation of meshing must be appended, giving three equations in the surface parameters and the cradle rotation. Solutions are obtained numerically and exported as point clouds, which are then fitted with splines in a CAD environment, converted to surfaces, trimmed against the blank, patterned around the axis, and assembled.
The assembly itself carries one non-obvious detail. The two members are positioned with a small negative offset of $-0.00635\ \mathrm{mm}$ between mating flanks, equal to the assumed elastic approach used in the contact ellipse calculation. This offset avoids the interference that would otherwise appear when the theoretical surfaces are brought into tangency, and it produces a contact patch in the motion simulation whose shape and location can be compared directly against the TCA prediction. Agreement between the two is a useful check that the surface equations, the discretization, and the assembly constraints are all consistent.
Contact Stress by Finite Element Analysis
The contact stress analysis was performed on a three-tooth segment of the assembled pair, a standard reduction that exploits the symmetry of the gear and the fact that no more than three tooth pairs share load at any instant. The material is 20CrMnTi case-hardening steel, with a density of $7800\ \mathrm{kg/m^3}$, a Poisson ratio of 0.25, and an elastic modulus of $2.1\times10^{5}\ \mathrm{MPa}$.
Meshing used free tetrahedral elements with a global size of 0.6 mm and a minimum size of 0.013 mm, producing 1,341,875 elements and 1,953,178 nodes. The contact pair was defined between the pinion flank as the target surface and the wheel flank as the contact surface, with frictional contact at a coefficient of 0.1 and the augmented Lagrangian algorithm. Two numerical controls required calibration rather than assumption. The normal stiffness factor $FKN$ was raised progressively from small values at a 15 N·m torque; beyond about 3 the maximum contact stress became insensitive to further increase, so $FKN=3$ was adopted. The penetration tolerance was reduced in the same empirical fashion and fixed at 0.02.
With the wheel fully constrained at its bore and sides, a drive torque of 200 N·m applied to the pinion produced a maximum contact stress of 1227.9 MPa. The predicted stress field is an irregular ellipse that decays outward from its peak, located in the region where the wheel addendum meets the pinion dedendum. Comparing this result with the allowable stress obtained from Hertzian theory gives the following check.
| Item | Value |
|---|---|
| Drive torque (N·m) | 200 |
| Maximum contact stress (MPa) | 1227.9 |
| Allowable contact stress (MPa) | 1280 |
| Deviation (%) | 4.07 |
A 4.07% margin is thin but positive, and it is of the order expected given the differences between a Hertzian line-contact idealization and a finite element solution of a curved, compliant, frictionally constrained contact. The result should be read as a confirmation that the design is not obviously understrength, not as a precise safety factor. Because bending stress at the pinion root reduces the effective contact pressure there, the peak stress on the pinion flank came out lower than on the wheel flank, which is the favorable distribution for avoiding concentrated damage.
A transient analysis extended the picture across the full mesh cycle. The pinion was driven at 1000 r/min through a rotational joint and the wheel loaded with a 200 N·m resisting torque. The contact stress migrates progressively from the heel toward the toe and from the tip toward the root during engagement, consistent with the contact path found in the motion simulation, and the peak value fluctuates rather than remaining constant. This cyclic variation, rather than the static peak, is what governs surface durability and is the quantity that a fatigue assessment would need to use.
Cutting and Roll Test
The final validation is physical. The gear pair was cut on a GH-35 face-milling machine of Japanese origin that had been retrofitted with servo drives in place of its mechanical feed train, improving positioning accuracy and repeatability. Blank mounting angle, eccentric angle, horizontal and vertical work offsets, roll ratio, and cradle angle were set from the computed values, and the teeth were cut on both members.
The cut pair was then mounted on a rolling tester, coated with red lead oxide, and run under light load so that contact between the flanks would transfer the marking compound and leave a visible footprint. On the wheel’s convex flank the pattern fell in the middle region of the tooth surface, with clearance at the toe, heel, and root, which is the intended outcome for a gear set of this type. The corresponding pattern on the pinion was of comparable size and position. Both matched the TCA prediction closely enough to support the conclusion that the settings derived through local synthesis were correctly transferred to the machine.
Two caveats apply to this verification. A roll test at light load shows the geometric contact pattern, not the loaded one, and the pattern under service torque will be larger and shifted. The test also confirms the settings collectively rather than isolating the effect of each preset, so it validates the procedure as a whole rather than any individual parameter choice.
What the Method Does and Does Not Control
The strength of local synthesis lies in making the second-order contact conditions explicit design inputs. Instead of accepting whatever curvature the conjugate construction produces at the reference point, the designer states the desired length and orientation of the contact path, the size of the contact ellipse, and the slope of the transmission error parabola, and the machine settings follow. Combined with TCA as an acceptance filter, this shortens the loop between a blank specification and a cuttable set of settings, and it makes the outcome predictable enough that the finite element and cutting stages become confirmations rather than discoveries.
The limitation is equally clear. Local synthesis controls behavior at and immediately around one reference point. It says nothing about third- and higher-order terms, so the shape of the transmission error curve away from the reference point, and the ease-off topography of the whole flank, remain unconstrained. This is precisely the gap that later work has addressed: H. J. Stadtfeld proposed replacing the parabolic transmission error with higher-order functions to improve dynamic behavior, and A. Fuentes and F. L. Litvin designed near-straight contact patterns and optimized second- and third-order modification coefficients to reduce sensitivity to misalignment and to lower mesh noise. Finite element studies such as those of S. Barone and co-workers quantified how misalignment and profile modification shift the loaded contact and the contact ratio.
Three practical extensions follow from the present results. A sensitivity study of the three presets, quantifying how $m’_{21}$, $\eta_2$, and $a$ each move the contact pattern and the error curve, would replace judgment with data at the start of the loop. Closing the chain from parameter calculation through solid modeling to finite element analysis in a single automated workflow would remove the manual file transfers that currently interrupt it. Finally, experimental measurement of root and flank strain on the cut pair would let the computed stress distribution be checked directly, rather than only through the indirect evidence of a satisfactory contact pattern.
For spiral bevel gears produced on existing machine tools, the significance of the procedure is that it does not require new hardware. It requires a reliable description of the generating process, a numerical evaluation of the meshing conditions, and a willingness to iterate on a small number of well-understood parameters until the predicted contact pattern and transmission error meet the application’s requirements.