Spiral bevel gears carry the final torque in vehicle drivelines, armored-vehicle transmissions and aviation power trains, and a tooth-root fracture in these components usually ends the life of the whole assembly. Design against bending fatigue therefore depends on knowing the stress-life relationship of the actual tooth root, not only the allowable stress of the material. For cylindrical gears this relationship is supported by standardized pulsator testing; for spiral bevel gears the same evidence base is thin, test rigs are costly, and a full-scale fatigue program for every variant in a product family is rarely affordable. The practical consequence is that bevel gear designs are often derived by analogy with existing parts, with limited quantitative confirmation of life.
A workable route around this bottleneck combines three elements: a material S-N curve obtained from specimens of the same heat and heat treatment as the gears, a transient three-dimensional contact model that resolves the stress history along the curved tooth root, and an accelerated rotating-bench test in which a short-life, high-load result is used to infer the long-life strength at service load. The evaluation summarized here applied this chain to a 16-tooth/27-tooth spiral bevel gear pair in 20CrNiMo steel, and cross-checked simulated lives against bench failures.
Test Gear Geometry and Manufacturing Quality
The pair was designed as a tapered-tooth (Gleason-type) spiral bevel set with a 90° shaft angle, following the geometric conventions used for cylindrical gear fatigue specimens. The pinion has 16 teeth (left hand) and the wheel 27 teeth (right hand), both with an outer transverse module of 4.25 mm and a face width of 17 mm. The remaining data are collected in Table 1.
| Parameter | Pinion | Wheel |
|---|---|---|
| Number of teeth, z | 16 | 27 |
| Hand of spiral | left | right |
| Outer transverse module, mt (mm) | 4.25 | |
| Face width, b (mm) | 17 | |
| Shaft angle (°) | 90 | |
| Normal pressure angle, αn (°) | 20 | |
| Mean spiral angle, β (°) | 35 | |
| Addendum coefficient | 0.85 | |
| Clearance coefficient | 0.188 | |
| Height modification coefficient | +0.252 | −0.252 |
| Tangential modification coefficient | +0.005 | −0.005 |
Machine-tool settings were computed with a loaded tooth contact analysis (LTCA) loop that treats the contact pattern and the transmission error as acceptance criteria: the wheel was cut by the duplex method and the pinion by the single-side method, with a cutter tip radius of 0.3 mm on both tools. Cutting parameters were adjusted until the predicted contact ellipse and the measured pattern coincided. After grinding, the teeth were inspected on a CNC gear measuring center and reached grade 6 of GB/T 11365—2019, and the rolled contact pattern reproduced the simulated one. That agreement matters for fatigue work because the position of the contact pattern governs the load distribution along the root fillet, and therefore which cross-section of the root carries the highest alternating stress.
Material Baseline: Tensile and Fatigue Strength of 20CrNiMo
Specimens were cut from the same batch of 20CrNiMo low-carbon alloy steel and underwent the same carburizing and hardening cycle as the gears, so that the specimen data could be transferred to the component analysis without a heat-treatment correction.
Five specimens were pulled to fracture in a universal testing machine. All failed in a brittle manner, with 6.5% elongation and less than 0.5% reduction of area, and a mean ultimate tensile strength of 1664.7 MPa. This places the steel at essentially the same strength level as AISI 9310, whose specification value is 1683.0 MPa — a reference point widely used for heavily loaded bevel gears in aviation applications.
Fatigue testing used a GPS200 high-frequency resonance machine at a stress ratio of 0.05 on 18 specimens, of which 16 produced valid results. Fracture surfaces showed negligible necking, faint dimples in the core and intergranular cracking in the carburized case, consistent with a case-controlled initiation mechanism. Regression of the valid points gave the 50% survival S-N curve, with 883.90 MPa at 3 × 106 cycles. Because a single tooth failure in a 16-tooth pinion corresponds to 15 of 16 teeth surviving, the component reliability implied by the test is 93.75%, so the material curve was shifted downward to a 93.75% confidence (survival) S-N curve. This shifted curve, not the mean curve, supplied the fatigue database used in the subsequent component simulation — a conservative choice that keeps the simulated life on the safe side of the specimen statistics.
Transient Contact Model of the Meshing Pair
The numerical tooth flanks were generated from the cutting kinematics rather than from an ideal involute-like approximation, and the root fillet was produced by sweeping the cutter tip radius, so the geometric stress concentration of the real ground fillet is represented in the model. The finite element model contained five teeth (two in mesh plus three following teeth), a choice justified by a contact ratio of 2.06, which allows up to three tooth pairs to share load simultaneously; the pinion concave flank was meshed against the wheel convex flank, the driving side under forward torque.
Mesh density was 0.7 mm on the tooth flanks and refined to 0.1 mm at the root fillet, where the bending stress gradient is steepest. Material inputs were a density of 7850 kg/m³, an elastic modulus of 1.8 × 1011 Pa, a Poisson ratio of 0.3 and a flank friction coefficient of 0.03. The operating point for the long-life case was a pinion speed of 900 r/min with a wheel torque of 499 N·m.
The transient solution returned the full bending stress history at the root. As a tooth moves from entering mesh at the heel to leaving mesh at the toe, the peak root stress migrates from the heel section through the mid-section to the toe section, so the root experiences a travelling load rather than a stationary one. This is the main reason a single nominal root stress from a standard formula cannot fully describe bevel gear bending damage, and the reason the fatigue assessment was carried out on the simulated stress-time history.
Fatigue Life Simulation at Two Load Levels
The stress history was processed in nCode with the 93.75% survival S-N curve assigned to the gear solid and Gerber correction applied for mean stress. Damage was accumulated with Miner’s rule, and the critical node — the node with the shortest life in the root region — was taken as the predicted failure site; contact-flank fatigue was excluded from the assessment. Because the target was a specified life rather than a specified load, the wheel torque was adjusted iteratively at 1 N·m resolution until the minimum root life matched the target.
Two life levels were computed. At T2 = 499 N·m the predicted life was 3.009 × 106 cycles against a 3.0 × 106 target, and at T2 = 570 N·m the predicted life was 1.0 × 106 cycles. The critical location lay at the mid-section of the root, slightly toward the heel, and moved further toward the heel as the load increased. Table 2 lists the associated stresses, together with the nCode stress and the tooth-root stress calculated according to ISO 10300-3 for the same load.
| Wheel torque T2 (N·m) | Life N (cycles) | Max. principal stress (MPa) | Max. equivalent stress (MPa) | nCode stress (MPa) | ISO 10300 stress (MPa) |
|---|---|---|---|---|---|
| 570 | 1.0 × 106 | 826.87 | 669.48 | 692.00 | 727.95 |
| 499 | 3.0 × 106 | 745.47 | 603.31 | 623.10 | 637.27 |
The nCode stress and the ISO 10300 result differ by 5.2% at the 106-cycle level and 2.3% at the 3 × 106-cycle level. Since the two are obtained by independent routes — one from the local stress history, the other from the standardized nominal stress with its own load-distribution factors — the proximity supports the credibility of the contact model. Plotting all four stress measures against life on logarithmic axes yields the simulated gear S-N lines, with slopes of −0.095465 (nCode) and −0.121097 (ISO).
Stress Acceleration Factor
Accelerated testing replaces a long-life, low-stress condition with a short-life, high-stress one, and the conversion depends entirely on the slope of the S-N line. In log-log coordinates the line is straight:
$$\lg \sigma = k \lg N + b$$
Taking two points on the line, with σx > σy corresponding to Nx < Ny, subtracting the two equations eliminates the intercept:
$$\frac{\sigma_x}{\sigma_y} = \left(\frac{N_x}{N_y}\right)^{k}$$
and the stress acceleration factor is defined as
$$a = \frac{\sigma_x}{\sigma_y}$$
For a factor of three in life (3 × 106 replaced by 1 × 106 cycles) the simulated slopes give a = 1.1106 with the nCode stress and a = 1.1423 with the ISO stress, corresponding to wheel torques of 570 N·m and 499 N·m respectively. Extending the same relation to a factor of ten in life with the ISO slope gives a = 1.3216. The factor is thus not a property of the gear but of the S-N slope, and it grows quickly as the extrapolation range widens; the reliability of a tenfold acceleration has not been demonstrated experimentally and should be treated as an extrapolation.
Bench, Lubrication and Running-In Protocol
A pulsator cannot load a spiral bevel tooth in a representative way, because the real tooth is loaded through a rolling-sliding motion with a migrating contact pattern. The tests therefore used a closed power flow (mechanically recirculating) durability bench developed for bevel gears, in which the load torque is set once and held by the loop stiffness. The stiffness of the housing and shaft system is far higher than that of the teeth, so misalignment of the mesh under load remains negligible and does not perturb the root stress.
Gears were installed to reproduce the design contact pattern. Lubrication was a GL-5 85W-90 heavy-duty gear oil, replaced once during the program, and the oil sump temperature was held between 45 °C and 55 °C. To limit the risk of pitting before the bending comparison, a stepped running-in schedule was applied: 1.0 × 105 cycles at quarter torque, 7.2 × 104 cycles at half torque and 3.6 × 104 cycles at three-quarter torque, a total of 2.0 × 105 cycles at 1200 r/min. The speed was then reduced to the test value of 900 r/min and the torque raised to the test level.
Test Results and Cross-Validation
The accelerated stage was run first at T2 = 570 N·m. Three pairs were mounted; the first showed pitting and did not fail by tooth breakage, so only the two bending failures were retained, giving lives of 6.46381 × 105 and 1.903183 × 106 cycles and a mean of 1.274782 × 106 cycles for a nominal 106-cycle target. The verification stage at T2 = 499 N·m accumulated more than 130 h of running on two pairs and produced lives of 1.653261 × 106 and 5.304999 × 106 cycles, a mean of 3.47913 × 106 cycles against a 3 × 106 target (Table 3).
| Wheel torque T2 (N·m) | Specimen | Life (105 cycles) | Mean life (105 cycles) |
|---|---|---|---|
| 570 | 1 | 6.46381 | 12.74782 |
| 2 | 19.03183 | ||
| 499 | 3 | 16.53261 | 34.7913 |
| 4 | 53.04999 |
The measured lives exceed the simulated targets by factors of 1.27 and 1.16, which corresponds to about 1.8% and 1.0% on a log-life basis; in other words the simulation is slightly conservative, as intended by the use of the 93.75% survival curve. Applying the simulated S-N slope to the 570 N·m mean life predicts a life of 3.824 × 106 cycles at 499 N·m, 1.63% above the 3 × 106 target in log terms, and the measured 3.479 × 106 cycles then falls within 0.6% of that prediction.
Post-test inspection showed that the teeth had reached a full-load contact condition, with the wear trace matching the preset pattern and no evident edge loading; all fractures initiated at the heel end of the pinion root, the location identified by the simulation. Converting the measured lives into allowable root bending stress and comparing with the simulation gives the deviations in Table 4.
| Life level | Test stress, nCode scale (MPa) | Simulated nCode stress (MPa) | Deviation | Test stress, ISO scale (MPa) | Simulated ISO stress (MPa) | Deviation |
|---|---|---|---|---|---|---|
| 1 × 106 | 709.77 | 692.00 | −2.5% | 751.75 | 727.95 | −3.2% |
| 3 × 106 | 632.82 | 623.10 | −1.5% | 649.91 | 637.27 | −1.9% |
| 3 × 106, inferred from the 106-cycle test | 639.10 | — | +1.0% vs. measured | 658.11 | — | +1.3% vs. measured |
Fitting the four valid test points produces the experimental gear S-N line, with slopes of −0.104454 (nCode scale) and −0.132503 (ISO scale), both steeper than the simulated lines. The corresponding threefold-life acceleration factors are 1.1216 and 1.1567, close to the simulated 1.1106 and 1.1423. Using the ISO 10300-3 evaluation, the measured tooth-root bending strength at 3 × 106 cycles is 649.91 MPa. The largest discrepancy between test and simulation across both stress scales and both life levels is 3.2%.
| Quantity | From simulated S-N line | From tested S-N line |
|---|---|---|
| Slope k, nCode stress | −0.095465 | −0.104454 |
| Slope k, ISO stress | −0.121097 | −0.132503 |
| Factor a, 3× life, nCode stress | 1.1106 | 1.1216 |
| Factor a, 3× life, ISO stress | 1.1423 | 1.1567 |
| Factor a, 10× life, ISO stress | 1.3216 | not established |
Interpretation and Limits of the Approach
The result that carries the most weight is not the 3.2% agreement itself but the consistency of the whole chain: a specimen-based survival S-N curve, transferred into a contact model that resolves the travelling root stress, predicts both the failure location and the life at two load levels, and the bench reproduces both. Because the simulation uses a 93.75% survival curve while the tests were run on only a few pairs, agreement at the level observed is partly a matter of statistics as well as mechanics.
Several constraints should be stated plainly. The evidence rests on one gear size and one material, with four valid bending failures in total, so the scatter of the underlying distribution is poorly characterized; the reported deviations describe this pair rather than a general bias for bevel gears. The 1 × 106 test level also produced a pitting failure on one pair, which indicates that at the accelerated load the contact-flank damage mode competes with bending and that the acceleration window is not unlimited. The threefold acceleration factor is supported by the test data; the tenfold factor is an extrapolation from the fitted slope and has no experimental confirmation. Finally, the comparison relies on mean stress correction (Gerber) and Miner’s rule, both of which introduce their own uncertainty when the load history is not purely constant-amplitude.
For design practice, the useful outcome is a procedure that can be executed in a reasonable time: establish the material curve on specimens from the production heat, calibrate a transient contact model against a measured contact pattern, iterate load to the required life, and then confirm with a short high-load bench run converted through the S-N slope. Where a full-scale program is impossible, this sequence gives a quantitative, traceable estimate of bevel gear bending strength instead of an analogy-based one, and it identifies the critical root section — at the heel side of the pinion in this case — that manufacturing and inspection effort should protect.