Spiral bevel gears transmit motion between intersecting shafts and are widely used in vehicle drivetrains, differentials, and final drives. Their curved teeth engage gradually, giving a high contact ratio, high load capacity, and relatively smooth power transfer. The same sliding–rolling contact that supports torque also makes the tooth surface vulnerable to wear. As material is removed, tooth thickness and profile change, mesh stiffness and backlash shift, and the contact pattern moves. The result is growing angular transmission error, vibration, noise, and loss of precision.
Much of the dynamic wear literature has focused on spur gears or on low-order lumped-parameter models. Those models do not capture the high contact ratio and load sharing of spiral bevel gears, where a worn tooth may be partly relieved by neighbouring teeth. For condition monitoring and tooth-profile design, the relevant question is how wear depth and operating load combine to change the dynamic transmission error.
Torsional Vibration Model and Wear Excitation
Transmission error can be represented as the relative displacement along the mesh line. For a pinion–gear pair,
$$
\delta = R_p\theta_p – R_g\theta_g
$$
where \(R_p\) and \(R_g\) are the reference pitch radii, and \(\theta_p\) and \(\theta_g\) are the torsional displacements of the pinion and gear about their shaft axes. The analysis uses a two-degree-of-freedom pure torsional model with rigid supports. Each gear body is represented by a concentrated inertia, and the mesh is represented by stiffness, damping, and backlash elements.
$$
I_p \ddot{\theta}_p + R_p c_m\left(\dot{\delta}-\dot{e}\right) + R_p k_m f\left(\delta-e\right) = T_p
$$
$$
I_g \ddot{\theta}_g + R_g c_m\left(\dot{\delta}-\dot{e}\right) + R_g k_m f\left(\delta-e\right) = T_g
$$
Here, \(I_p\) and \(I_g\) are the rotary inertias, \(T_p\) and \(T_g\) are the driving and load torques, \(e\) is the static transmission error, \(c_m\) is the mesh damping, \(k_m\) is the time-varying mesh stiffness, and \(f\) is the backlash function. Defining the dynamic mesh displacement as \(x=\delta-e\) gives the combined form
$$
m_e \ddot{x} + c_m \dot{x} + k_m f(x) =
m_e\left(\frac{R_p T_p}{I_p}+\frac{R_g T_g}{I_g}\right)-m_e\ddot{e}
$$
where \(m_e\) is the equivalent mass of the gear pair. The equation retains time-varying stiffness, backlash, and static transmission error, so the response is nonlinear. Tooth wear enters the model as a controlled reduction of tooth thickness. The pinion is treated as the fault component because, for the same material and service period, it accumulates more contact cycles than the gear. Single-tooth and full-tooth wear cases are simulated at wear depths \(\lambda\) of 25, 50, 75, and 100 μm.
Gear Geometry, Material, and Simulation Conditions
The gear pair has a 14-tooth pinion and a 41-tooth gear. Its main geometry is listed in Table 1.
| Parameter | Pinion | Gear |
|---|---|---|
| Number of teeth | 14 | 41 |
| Face width / mm | 55 | 55 |
| Pitch diameter / mm | 111.72 | 327.18 |
| Module at large end / mm | 7.98 | 7.98 |
| Spiral angle / ° | 35 | 35 |
| Pressure angle / ° | 20 | 20 |
The solid model uses 45Cr steel. Its elastic modulus is 209,000 MPa, Poisson’s ratio is 0.3, and density is 7,860 kg·m⁻³. Rayleigh damping is included with a mass-proportional coefficient \(\alpha=0.03\) and a stiffness-proportional coefficient \(\beta=3\times10^{-6}\). The mesh remains in frictional contact with a friction coefficient of 0.1. The pinion speed is held at 210 r/min, and the torque on the driven gear is varied from 500 to 3,000 N·m. The main comparisons use 500, 1,000, 2,000, and 3,000 N·m; the full-tooth fit also includes a 1,500 N·m case. Dynamic transmission error is extracted from the computed pinion and gear rotation angles over the first 0.13 s of response.
Single-Tooth Wear: Local Excitation, Peak Error, and Load Sharing
For a single worn tooth at 500 N·m, the response settles into steady state after a short transient. When the worn tooth enters mesh, the local stiffness and backlash conditions produce a stronger excitation. The instantaneous dynamic transmission error \(\gamma\) rises in magnitude to a peak \(\Gamma\), then falls as the worn tooth leaves the contact zone. Increasing the wear depth \(\lambda\) increases both the instantaneous error and the peak magnitude. The curves for 25, 50, 75, and 100 μm show this monotonic growth in the initial wear range.
The high contact ratio of the spiral bevel pair changes the trend at larger wear depths. Once \(\lambda\) exceeds a critical level, the worn tooth no longer carries its original share of load, and neighbouring teeth in the edge region take over part of the mesh force. The peak error then approaches a load-dependent plateau \(\gamma_1\). This saturation is a load-sharing effect, not an absence of damage: the worn tooth has lost its normal contact condition, but the remaining teeth limit the further rise in transmission error.
Load has two competing effects. A higher torque \(P\) produces larger elastic deformation, so for the same wear depth the instantaneous error \(\gamma\) and the plateau \(\gamma_1\) both increase. At the same time, the deformed contact can compensate part of the geometric disturbance, reducing fluctuation around the mean response, shortening the transient, and smoothing the curve. The peak \(\Gamma\) is also delayed relative to the low-load case.
At a fixed load, the peak magnitude \(\Gamma\) increases monotonically with wear depth \(\lambda\) in the initial range. The slope \(\varepsilon_1\) of the fitted peak-error curve then falls to zero at the critical wear level, marking the plateau \(\gamma_1\). Increasing the load increases the absolute values of both the slope \(\varepsilon_1\) and the intercept \(\eta_1\), and raises the plateau \(\gamma_1\). A load-dependent threshold is therefore needed when peak transmission error is used as a wear indicator.
Full-Tooth Wear: Mean Error and Empirical Fit
When all teeth wear equally, every mesh cycle is affected. At 500 N·m, increasing the wear depth from 25 to 100 μm increases the absolute dynamic transmission error. Higher load drives the response into steady state more quickly and reduces volatility, but it also produces a larger absolute error. Because no healthy tooth can take over the load of a worn tooth, the full-tooth fault produces a more stable error signal than the single-tooth case.
The mean dynamic transmission error \(E\) varies approximately linearly with wear depth \(\lambda\). The fitted lines for different loads are nearly parallel, which means that wear mainly shifts the mean error while load changes its offset. The load-specific fit parameters are listed in Table 3.
| Load \(P\) / N·m | Slope \(\varepsilon_2\) | Intercept \(\eta_2\) / μm |
|---|---|---|
| 500 | −1.332 | −10.62 |
| 1,000 | −1.463 | −28.33 |
| 1,500 | −1.472 | −49.26 |
| 2,000 | −1.475 | −60.27 |
The linear form is
$$
E = \varepsilon_2 \lambda + \eta_2
$$
A combined fit across the full-tooth cases gives
$$
E = -1.475\lambda – 0.021P – 3.025
$$
with \(E\) and \(\lambda\) in μm and \(P\) in N·m. This expression is useful as an engineering estimate of the mean error under full-tooth wear. When the load is known, the load-specific coefficients in Table 3 are preferable because they preserve the local fit more closely.
Design and Monitoring Implications
The two wear modes produce different signatures. Single-tooth wear creates a local, phase-locked peak in dynamic transmission error. Its growth saturates because the high contact ratio allows adjacent teeth to share load. A global RMS value may therefore understate severe local damage. Phase-resolved or angle-domain peak error is more sensitive to a single worn tooth, especially when the worn tooth can be identified in the mesh cycle.
Full-tooth wear shifts the mean transmission error in a nearly linear manner. The mean value is a robust indicator of distributed wear, but it must be normalized against load because the intercept moves with torque. The fitted slope becomes more negative as load increases, so a fixed alarm threshold will not be valid across the operating range.
- For local tooth damage, monitor peak and phase-resolved transmission error rather than only the mean.
- For distributed wear, use the mean error with load compensation and load-specific baseline coefficients.
- Tooth-profile design and profile modification can delay the critical wear level by maintaining favourable contact and load sharing.
- High contact ratio improves load sharing but can mask local wear in global metrics, so multi-feature monitoring is preferable.
The analysis is limited to a two-degree-of-freedom torsional formulation with rigid supports. It does not resolve shaft and bearing flexibility, thermal effects, or lubricant-film behaviour. Absolute error values and alarm thresholds for a specific axle should therefore be confirmed experimentally, while the load-dependent trends provide a basis for tooth-profile design and wear-sensitive monitoring of spiral bevel gears.