Gear root bending strength: statistical treatment of Single Tooth Bending Fatigue tests results

Source: Politecnico di Milano, Italy
Gear tooth failure due to tooth root bending fatigue is one of the most dangerous failure modes in gears. Therefore, the precise definition of gear bending fatigue strength is a key aspect. Indeed, STBF tests have been developed in order to estimate the gear load carrying capacity by testing directly the spur gear teeth. Here, with a specific focus on symmetric STBF test, a methodology based on two different statistical tools (i.e. Maximum Likelihood and statistic of extreme) is proposed with a view to estimate the gear tooth root bending load carrying capacity.
Tooth root bending fatigue is considered as one of the most dangerous gear failure mode because, typically, it results in the interruption of the power flow within the gearbox. Standards, such as ISO 6336-3 [1] and ANSI/AGMA 2001 [2], provide designers with an analytical framework, to assess a gear train, in respect to this failure mode, by comparing the maximum tooth root stress with the limit value of the gear itself. Data of several typical materials are present within the standards (e.g. [2], [3]), defined at 1% gear failure probability (or 99% reliability) for tooth root bending fatigue. Different reliability level can be reached using literature coefficients (e.g. [2], [4]).
However, when specific material data are not available, extensive gear testing is needed, especially if a certain reliability level, after a given service life, must be assured.
Within this context, a gear SN estimation approach, based two different statistical tools (i.e. maximum likelihood method and statistic of extreme) will be shortly presented here. More details on the proposed statistical methodology will be presented at the International Conference on Gears 2022.
Short remarks about STBF gear testing
Tooth root bending fatigue strength can be evaluated by means of running gears tests (e.g. [7], [8]), Single Tooth Bending Fatigue (STBF) tests (e.g. [9]–[18]), and tests on notched specimens (e.g. ISO 6336-3 annex A [1], [19]–[21]). Amongst them, due to its higher efficiency and capability to directly work with gears, STBF tests are the most adopted [22].
STBF test are performed in two configurations: asymmetric and symmetric. The first, mainly adopted in the U.S. [8], [12], is described in the SAE recommended practice J1619 [15]. In a asymmetric configuration, one tooth is directly loaded by the machine, while a second tooth and a pin work as supports. Here, when the loaded tooth breaks, the test is stopped.

Figure 1: Example of symmetric STBF test rig of an aerospace grade gear.
The symmetric STBF test configuration is the one mainly adopted in Europe [9], [10], [13]. Here, both tooth roots are subject to the same nominal toot root stress. In a symmetric STBF test configuration, each test is stopped when one of the two teeth breaks, as it is physically impossible to load again the same tooth roots. An example of symmetric STBF test rig is shown in Figure 1.
In the last two decades, the authors have estimated the tooth root bending fatigue load carrying capacity via a symmetric STBF approach for different kind of gear specimen. Some examples: gears for aeronautical applications [24], [25] , for small planetary gearboxes [26], for wind turbine gearboxes [27], [28] also additively manufactured gears [23], [29]. All of them have been performed on a Schenck pulsator, adopting R=0.1. For more details concerning the adopted testing procedure, the interested reader is referred to the aforementioned references
The classical approach to elaborate STBF test data
Unfortunately, STBF tests results cannot be directly applied to assess a gear pair. Indeed, the loading condition of a STBF test rig is not representative of the real case. The literature provides two different approaches for the treatment of STBF test data. Therefore, they can be used within the existing calculation methods. They are proposed within FVA report no. 304 [10] and within the works of Rao and Mc Pherson [8], [12].
Both authors highlight two main issues that must be considered in order to treat STBF test data:
- The testing and service failures are different. In STBF test, a failed test corresponds to the breakage of some tested teeth, while in the meshing case, the failure of the specimen (i.e. the gear itself) corresponds to the failure of the weakest tooth.
- The testing and service loading conditions are different [30]. In the real case, the force varies both, the application point, and its magnitude, while in STBF the force is a sinusoidal load whit a fixed amplitude and a given application point.
Stahl and Rao faced the problem with two different approaches. On the one hand, Stahl, K. [10] proposed several coefficients, based on experimental evidence: three different coefficients deal with the statistical difference: f_(p→Z) and f_(Z→ZR) translates STBF results to the case of a meshing gear while f_(1%ZR) is adopted to reach the 99% reliability level. A further coefficient f_korr based on Rettig experimental evidence [9] is used to adress the different loading conditions (the same is adopted here).
On the other hand, Rao and Mc Pherson [8], [12] define a simple statistical relation between the teeth and the gear, and by means of a ‘fit by eye’, they translate STBF test data, corrected to the real case load, to the desired reliability level. Concerning the different loading conditions, Allowable Stress Range (ASR) diagrams are used.
The proposed approach: teeth SN curve estimation by means of maximum likelihood

Figure 2: Typical gear SN curve
In the gear field, both ISO 6336-3 and ANSI-AGMA 2001-D04 reports classical gear SN curve, suggesting that gears do not present a fatigue limit (see Figure 2). However, the standards do not suggest a typical slope for the long-life region. Anyhow, in the gear field, the estimation of such secondary slope is becoming relevant if the ass
essment is based on a load spectra calculation, thus requiring a proper estimation of both slopes.
Therefore, the developed methodology proposes a model based on a curve shape compatible with the ones proposed by the standard, i.e. a curve presenting two slope, one for the limited life region, one for the long life one. This kind of curve has been firstly proposed by Spindel et. al. [31]:

where N is the number of cycles, S is the applied load, N_e and S_e represent the knee position, k is the slope of the limited life region and k_1 is the slope of the long-life region.
A first attempt to estimate this curve can be done by simply considering all the data as a failure and then fitting this curve. However, through a deep analysis of symmetric STBF test data, it is possible to notice that each experimental point “hides” the presence of another loaded tooth, which is survived. In other words, in a symmetric STBF test rig, each test interrupted by the failure of one of the two teeth is actually a failure and a survival. On the other hand, the case of a runout represents two runouts. This important presence of survived teeth can be properly handled only adopting MLE as a tool for the estimation of the teeth SN curve.
Indeed, MLE is an estimation technique typically applied in the determination of the S-N curve, as it has great advantages to consider, at the same time, different kind of data. Within MLE, exact data (e.g. failures) and censored points (e.g. runouts and survivals) are considered, with a different statistical meaning, within the same calculation procedure [6], [32], [33]. MLE has been applied for the estimation of the S-N curve of specimens (e.g. [31], [34]–[36]), as well as in other components, e.g. welded structures (e.g. [37], [38]).

Figure 3: Different estimated teeth SN curves.
The proposed methodology applies MLE in order to estimate the experimental S-N curve, that is the S-N curve of the teeth, and its associated dispersion.
The proposed approach: statistic of extreme to estimate the gear SN curve
Once the S-N curve of the teeth has been estimated, STBF result elaboration is necessary to make it representative of the gear. In the real case scenario, where gears are meshing and rotating, gear failure due to tooth root bending fatigue is ruled by failure of its weakest gear tooth [8], [9], [12].Therefore, the gear S-N curve can be defined as the one describing the weakest tooth amongst the z gear teeth. In other words, the CDF (Cumulative Density Function) of the teeth has to be elaborated in order to define the gear CDF, that is the CDF of the weakest tooth among the z gear teeth.
By means of a mathematical passage [6], statistic of extremes allows to define the CDF of the minimum value over n extractions of X. Focusing on the case of gear tooth root bending fatigue, it can be used to define the distribution of the tooth with the smallest resistance over z gear teeth.
The article header image explains graphically the procedure: the STBF CDF is elaborated, by means of statistics of extremes, to define the gear CDF.

Figure 4: Estimated gear SN curves
For gear failure due to tooth root bending fatigue, ISO 6336-5 [3] proposes a typical failure probability of 1% (i.e. a reliability of 99%), hence it is possible to define the curve at the required reliability level. In other words, estimating the load carrying capacity with the 99% percent probability of being exceeded for the weakest tooth over z teeth. However, different SN curve, at different reliability levels can be obtained by calculating different percentiles.
[1] ISO 6336-3:2019 - Calculation of load capacity of spur and helical gears — Part 3: Calculation of tooth bending strength. .
[2] ANSI~Standard, “Fundamental Rating Factors and Calculation Methods for Involute Spur and Helical Gear Teeth,” Am. Gear Manuf. Assoc. Alexandria, 2004.
[3] ISO~6336-5, “Calculation of Load capacity of Spur and Helical Gears, Part 5: Strength and quality of materials,” Int. Stand. Organ. Geneva, 2006.
[4] M. Hein, M. Geitner, T. Tobie, K. Stahl, and B. Pinnekamp, “Reliability of gears - Determination of statistically validated material strength numbers,” 2018.
[5] ISO 6336-6:2019, “Calculation of Load capacity of Spur and Helical Gears, Part 6: Calculation of service life under variable load,” Geneva, CH, 2019.
[6] S. Beretta, Affidabilità delle costruzioni meccaniche: Strumenti e metodi per l’affidabilità di un progetto. Springer Science & Business Media, 2010.
[7] I. J. Hong, A. Kahraman, and N. Anderson, “A rotating gear test methodology for evaluation of high-cycle tooth bending fatigue lives under fully reversed and fully released loading conditions,” Int. J. Fatigue, vol. 133, p. 105432, 2020.
[8] S. B. Rao and D. R. McPherson, “Experimental characterization of bending fatigue strength in gear teeth.,” Gear Technol., vol. 20, no. 1, pp. 25–32, 2003.
[9] H. Rettig, “Ermittlung von Zahnfußfestigkeits-Kennwerten auf Verspannungsprüfständen und Pulsatoren--Vergleich der Prüfverfahren und gewonnenen Kennwerte,” Antriebstechnik, vol. 26, pp. 51–55, 1987.
[10] K. Stahl, “Lebensdauerstatistik : Abschlussbericht : Forschungsvorhaben Nr. 304,” FVA, 1999.
[11] M. Benedetti, V. Fontanari, B.-R. Höhn, P. Oster, and T. Tobie, “Influence of shot peening on bending tooth fatigue limit of case hardened gears,” Int. J. Fatigue, vol. 24, no. 11, pp. 1127–1136, 2002.
[12] D. R. McPherson and S. B. Rao, “Methodology for translating single-tooth bending fatigue data to be comparable to running gear data,” Gear Technol., 2008.
[13] T. Tobie and P. Matt, “Empfehlungen zur Vereinheitlichung von Tragfähigkeitsversuchen an vergüteten und gehärteten Zylinderrädern (Technical Report FVA-Richtlinie 563-I),” 2012.
[14] F. Dobler, T. Tobie, and K. Stahl, “Influence of low temperatures on material properties and tooth root bending strength of case-hardened gears,” in International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, 2015, vol. 57205, p. V010T11A007.
[15] S. SAE, “J1619, Single Tooth Gear Bending Fatigue Test.” SAE International, 2017.
[16] C. Güntner, T. Tobie, and K. Stahl, “Influences of the Residual Stress Condition on the Load-Carrying Capacity of Case-Hardened Gears,” Am. Gear Manuf. Assoc. AGMA Tech. Pap. 17FTM20, 2017.
[17] K. J. Winkler, S. Schurer, T. Tobie, and K. Stahl, “Investigations on the tooth root bending strength and the fatigue fracture characteristics of case-carburized and shot-peened gears of different sizes,” Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci., vol. 233, no. 21–22, pp. 7338–7349, 2019.
[18] J. Koenig, S. Hoja, T. Tobie, F. Hoffmann, and K. Stahl, “Increasing the load carrying capacity of highly loaded gears by nitriding,” in MATEC Web of Conferences, 2019, vol. 287, p. 2001.
[19] D. J. Medlin, B. E. Cornelissen, D. K. Matlock, G. Krauss, and R. J. Filar, “Effect of thermal treatments and carbon potential on bending fatigue performance of SAE 4320 gear steel,” SAE Trans., pp. 547–556, 1999.
[20] J. J. Spice, D. K. Matlock, and G. Fett, “Optimized carburized steel fatigue performance as assessed with gear and modified Brugger fatigue tests,” SAE Trans., pp. 589–597, 2002.
[21] L. Vilela Costa, D. de Oliveira, D. Wallace, V. Lelong, and K. O. Findley, “Bending Fatigue in Low-Pressure Carbonitriding of Steel Alloys with Boron and Niobium Additions,” J. Mater. Eng. Perform., pp. 1–10, 2020.
[22] D. R. McPherson and S. B. Rao, “Mechanical Testing of Gears.,” Mater. Park. OH ASM Int. 2000., pp. 861–872, 2000.
[23] F. Concli et al., “Bending fatigue behavior of 17-4 ph gears produced by additive manufacturing,” Appl. Sci., vol. 11, no. 7, p. 3019, 2021, doi: 10.3390/app11073019.
[24] G. Gasparini, U. Mariani, C. Gorla, M. Filippini, and F. Rosa, “Bending fatigue tests of helicopter case carburized gears: Influence of material, design and manufacturing parameters,” in American Gear Manufacturers Association (AGMA) Fall Technical Meeting, 2008, pp. 131–142.
[25] C. Gorla, F. Rosa, E. Conrado, and F. Concli, “Bending fatigue strength of case carburized and nitrided gear steels for aeronautical applications,” Int. J. Appl. Eng. Res., vol. 12, no. 21, pp. 11306–11322, 2017.
[26] C. Gorla, E. Conrado, F. Rosa, and F. Concli, “Contact and bending fatigue behaviour of austempered ductile iron gears,” Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci., vol. 232, no. 6, pp. 998–1008, 2018.
[27] C. Gorla, F. Rosa, E. Conrado, and H. Albertini, “Bending and contact fatigue strength of innovative steels for large gears,” Proc. Inst. Mech. Eng. Part C J. Mech. Eng. Sci., vol. 228, no. 14, pp. 2469–2482, 2014.
[28] C. Gorla, F. Rosa, F. Concli, and H. Albertini, “Bending fatigue strength of innovative gear materials for wind turbines gearboxes: Effect of surface coatings,” in ASME 2012 International Mechanical Engineering Congress and Exposition, 2012, pp. 3141–3147.
[29] L. Bonaiti, F. Concli, C. Gorla, and F. Rosa, “Bending fatigue behaviour of 17-4 PH gears produced via selective laser melting,” in Procedia Structural Integrity, 2019, vol. 24, doi: 10.1016/j.prostr.2020.02.068.
[30] L. Bonaiti, A. B. M. Bayoumi, F. Concli, F. Rosa, and C. Gorla, “Gear root bending strength: a comparison between Single Tooth Bending Fatigue Tests and meshing gears,” J. Mech. Des., pp. 1–17, 2021, doi: 10.1115/1.4050560.
[31] J. E. Spindel and E. Haibach, “Some considerations in the statistical determination of the shape of SN curves,” in Statistical analysis of fatigue data, ASTM International, 1981.
[32] W. B. Nelson, Applied life data analysis, vol. 521. John Wiley & Sons, 2003.
[33] W. B. Nelson, Accelerated testing: statistical models, test plans, and data analysis, vol. 344. John Wiley & Sons, 2009.
[34] S. Beretta, P. Clerici, and S. Matteazzi, “The effect of sample size on the confidence of endurance fatigue tests,” Fatigue Fract. Eng. Mater. Struct., vol. 18, no. 1, pp. 129–139, 1995.
[35] S. Lorén, “Fatigue limit estimated using finite lives,” Fatigue Fract. Eng. Mater. Struct., vol. 26, no. 9, pp. 757–766, 2003.
[36] F. G. Pascual and W. Q. Meeker, “Estimating fatigue curves with the random fatigue-limit model,” Technometrics, vol. 41, no. 4, pp. 277–289, 1999.
[37] K. Wallin, “The probability of success using deterministic reliability,” in European Structural Integrity Society, vol. 23, Elsevier, 1999, pp. 39–50.
[38] G. Marquis and T. Mikkola, “Analysis of welded structures with failed and non-failed welds based on maximum likelihood,” Weld. World, vol. 46, no. 1, pp. 15–22, 2002.
Author
Luca Bonaiti, M. Sc., Department of Mechanical Engineering, Politecnico di Milano, Italy