Spiral bevel gears transmit torque between intersecting shafts in vehicle axles, marine drives, and aerospace transmissions. Their curved teeth engage gradually, giving the mesh high load capacity and smoother motion than straight bevel gears. Smoothness under load, however, depends on more than nominal tooth geometry. Torque bends and compresses the teeth, the contact ellipse migrates across the flank, and the driven gear no longer follows its theoretical kinematic angle. The resulting angular error—the loaded transmission error (LTE)—acts as a displacement excitation at the mesh. It generates dynamic mesh force and, through the shaft, bearing, and housing structure, produces vibration and noise.
A representative 16/27-tooth spiral bevel gear pair illustrates the mechanism. The pinion has 16 teeth and the gear 27, with a 4.25 mm module, 35° mid-face spiral angle, 20° normal pressure angle, and 90° shaft angle. The pinion is left-hand and the gear right-hand. These proportions are typical of a final drive: the low pinion tooth count produces a large ratio, while the spiral angle provides overlapping contact and a more favorable contact pattern.
| Parameter | Pinion | Gear |
|---|---|---|
| Number of teeth | 16 | 27 |
| Module (mm) | 4.25 | 4.25 |
| Mid-face spiral angle (°) | 35 | 35 |
| Normal pressure angle (°) | 20 | 20 |
| Shaft angle (°) | 90 | 90 |
| Hand | Left | Right |
| Pitch cone angle (°) | 30.650 7 | 59.349 3 |
| Addendum (mm) | 4.682 | 2.542 |
| Dedendum (mm) | 3.342 | 5.482 |
| Clearance (mm) | 0.8 | 0.8 |
Building the finite element model from the tooth surface itself
The flank of a spiral bevel gear is not a simple analytic surface. It is the envelope produced by a rotating cutter as the workpiece and cutter move through a generating motion. Machine-tool settings define the cutter radius, radial and angular cutter positions, roll ratio, and blank mounting angle. Evaluating the tooth-surface equations on a grid gives the position and normal vector of each flank point. In a rotation-projection plane, the grid nodes satisfy
$$
x_2(u_{ij},v_{ij})=L_{ij},\qquad
\sqrt{y_2^2(u_{ij},v_{ij})+z_2^2(u_{ij},v_{ij})}=R_{ij},
$$
for \(i=1,\ldots,n\) and \(j=1,\ldots,m\), where \(L_{ij}\) is the coordinate along the projection plane, \(R_{ij}\) is the radial coordinate, and \(u\) and \(v\) are the two surface parameters. The indices \(n\) and \(m\) denote the number of grid divisions along tooth height and tooth length.
Using these calculated points directly as finite element nodes preserves the manufactured tooth topology. A CAD solid is not generated and then re-meshed, so the micron-scale flank modifications that control contact are not lost in geometry cleanup. Hexahedral elements are used because the contact problem is geometrically nonlinear and because the teeth must resist distortion under high contact pressure. The root transition is refined with six node layers, since bending fatigue is highly sensitive to the local stress concentration where the tooth fillet meets the root.
Tooth contact analysis (TCA) supplies the second ingredient: the actual meshing position. TCA computes the contact path and the geometric transmission error for the unloaded pair. From these results, the instantaneous angular positions of pinion and gear are obtained and transformed into a common coordinate system. A five-tooth segment is then assembled, with the neighboring teeth retained to represent the stiffness of the surrounding gear body without the cost of a full-wheel model.
The contact model uses alloy steel with Young’s modulus \(2.1\times 10^5\) MPa, Poisson’s ratio 0.3, and density \(7.9\times 10^{-9}\) t/mm³. The gear flank is defined as the slave surface and the pinion flank as the master surface to limit penetration during contact iteration. The pinion is fully constrained, while the gear is free to rotate about its axis and receives the applied torque. Two implicit static steps are used: a 1 s step closes the backlash and establishes contact, and a 3 s step applies the load. This sequencing avoids an initial rigid-body mismatch and improves convergence of the contact solution.
Contact stress and load sharing
At 30 N·m, the contact ellipse moves from the heel and addendum region toward the toe and dedendum as the mesh progresses. The maximum principal stress occurs near the center of the contact ellipse. When only one tooth pair carries the load, the contact stress reaches its highest value; when two pairs share the load, the stress in each contact falls. This single-tooth/double-tooth alternation is the basic source of mesh stiffness variation and thus of dynamic excitation.
The contact pattern remains inside the flank at 30 N·m. At the exit of the middle tooth, the maximum stress sits just below the addendum line rather than at the tooth edge. That condition is favorable: an interior contact ellipse distributes pressure over the flank and avoids the stress concentration that appears when the contact runs off the tooth boundary. At higher torque the situation changes, because the loaded deformation shifts the contact pattern toward the edge.
Loaded transmission error is not monotonic with torque
The LTE is the angular difference between the actual driven-gear position and the theoretical position under load. It is produced by bending of the teeth, local contact deformation, and deflection of the gear body. At a pinion rotation angle \(\varphi\), the LTE can be written as
$$
E_L(\varphi)=\delta\varphi(\varphi)+E_T(\varphi),
$$
where \(E_T\) is the geometric transmission error obtained from TCA and \(\delta\varphi\) is the load-induced angular deformation extracted from the finite element contact solution. The first term describes the unloaded kinematic irregularity designed into the flank; the second describes how the mesh compliance changes that irregularity under torque.
For the 16/27-tooth pair, the peak-to-peak fluctuation of the LTE is 18.64 μrad at 10 N·m, 11.18 μrad at 30 N·m, and 31.24 μrad at 50 N·m. The trend is non-monotonic: the fluctuation first decreases and then increases. At 50 N·m, the LTE reaches a maximum of 259.15 μrad, and the loaded curve lies above the right-hand endpoint of the geometric transmission error curve. That condition marks the threshold of edge contact. Once the contact pattern reaches the tooth boundary, the maximum stress moves toward the tooth tip or root, the mesh stiffness variation increases, and the transmission becomes less stable.
| Torque (N·m) | LTE fluctuation (μrad) | Contact and vibration observation |
|---|---|---|
| 10 | 18.64 | Interior contact; first-order vibration higher than at 30 N·m and lower than at 50 N·m |
| 30 | 11.18 | Minimum fluctuation; measured axial acceleration 0.1097 m/s² |
| 50 | 31.24 | LTE peak 259.15 μrad; near edge contact; measured axial acceleration 0.1802 m/s² |
The 30 N·m case is therefore the best operating point within this load series. The result does not mean that lower torque is always smoother. Reducing torque changes the contact area and the share of load between neighboring teeth, so the LTE fluctuation can rise again at light load. The design target is not minimum load or maximum load in isolation, but the load at which the contact pattern and mesh stiffness variation produce the smallest angular error fluctuation.
From transmission error to mesh force and vibration
The time-varying LTE generates a dynamic mesh force. To evaluate the system response, the LTE curve is used as the excitation in a gearbox model built in MASTA. The response is examined in an order-domain waterfall covering the first three mesh orders. The first mesh order dominates the acceleration response and therefore contributes most to the overall vibration signature.
At 30 N·m, the first-order response is the smallest of the three load cases. The 10 N·m and 50 N·m cases produce successively larger maxima. The important feature is not simply the overall level but where the response peaks. At 30 N·m, the first-order mesh response shows clear peaks when the excitation aligns with the natural frequencies at 1.66 kHz and 2.19 kHz. The corresponding accelerations are 131.91 m/s² and 128.93 m/s². At 2.10 kHz, the mesh excitation also coincides with a natural frequency, but it does not produce severe vibration.
| Frequency (kHz) | Modal compliance condition | Phase relation | Dynamic mesh force | First-order response |
|---|---|---|---|---|
| 1.66 | Equal magnitudes | Opposite | Peak | 131.91 m/s² |
| 2.10 | Unequal magnitudes | Not opposite | No strong net peak | No severe vibration |
| 2.19 | Equal magnitudes | Opposite | Peak | 128.93 m/s² |
The distinction is mechanically significant. Frequency coincidence identifies a potential resonance region, but it does not by itself guarantee severe vibration. In this model, a strong dynamic mesh-force peak appears only when the two modal compliance magnitudes are equal and their phases are opposite. At 2.10 kHz, the compliance magnitudes are unequal and the phase difference is not 180°, so the modal contributions do not add into a large net mesh force. Resonance assessment based only on a natural frequency matching the mesh frequency can therefore overpredict risk; the compliance and phase relation between mesh-force contributions must also be considered.
The load ranking is consistent with the LTE result. The 30 N·m case has the smallest LTE fluctuation and the smallest first-order vibration maximum. The 50 N·m case has the largest fluctuation and the largest vibration maximum. Reducing the peak-to-peak LTE therefore reduces the excitation that drives the gearbox response.
Experimental trend check
Validation was carried out on a closed power-flow bevel gear test rig using an M+P dynamic signal analyzer. An accelerometer was mounted on the pinion shaft close to the meshing position to measure axial vibration. Tests were run at an input speed of 540 r/min with gear loads of 10, 30, and 50 N·m. At 540 r/min, the 16-tooth pinion gives a fundamental mesh frequency of 144 Hz.
The measured spectra contain the mesh frequency and its harmonics, with the first mesh order producing the largest acceleration. This agrees with the order-domain waterfall result. Under the same speed, the first-order acceleration is smallest at 30 N·m, 0.1097 m/s², and largest at 50 N·m, 0.1802 m/s². The 10 N·m case lies between them. The acceleration therefore decreases from 10 to 30 N·m and then increases from 30 to 50 N·m, following the same non-monotonic trend as the LTE fluctuation.
The comparison is one of trend rather than absolute amplitude. The test measures the axial response of the shaft structure, whereas the vibration model uses the LTE as a mesh excitation and evaluates the resulting system acceleration. Absolute levels depend on the measurement point, mounting stiffness, and structural damping. The load ranking, however, is a meaningful check: the operating point with the smallest loaded transmission error also produces the smallest measured vibration.
Design consequences
The first design consequence concerns the geometric transmission error. Flank modifications are often evaluated in the unloaded state, where the TCA curve is clean and the contact pattern is easy to read. Under torque, the same flank deflects, the contact pattern moves, and the effective transmission error changes. A geometric transmission error that is well behaved at light load can produce a large LTE fluctuation at service load. The design target should therefore be the loaded curve at the dominant operating torque, not only the unloaded TCA curve.
The second consequence concerns edge contact. The 50 N·m case shows how a load increase can push the contact pattern toward the tooth boundary. Once edge contact begins, the contact stress concentrates at the tip or root, the mesh stiffness variation grows, and the LTE fluctuation rises sharply. Maintaining an interior contact ellipse under the highest expected torque is more effective than compensating for edge contact after it appears.
The third consequence concerns resonance assessment. A natural frequency coinciding with a mesh harmonic is a necessary warning, but not a sufficient prediction of severe vibration. The modal compliance magnitudes and phases determine whether the dynamic mesh force adds or cancels. In this spiral bevel gear pair, the 2.10 kHz coincidence does not produce a strong response, while the 1.66 kHz and 2.19 kHz coincidences do. Design iterations should therefore examine the mesh-force summation condition, not just the frequency coincidence.
Because LTE is load-dependent, a gear set optimized at one torque is not necessarily optimal across a duty cycle. The 30 N·m optimum in the 10–50 N·m series is a useful illustration: it is neither the lightest nor the heaviest load, but the condition at which the contact pattern, mesh stiffness, and angular error fluctuation combine most favorably. For automotive, marine, and aerospace applications, the practical target is to weight the geometric transmission error design toward the torque range that dominates service, while checking that light-load and high-load edges of the duty cycle remain free of unstable contact and excessive vibration.