Research Article  Open Access
Reliability Modeling and Estimation of the Gear System
Abstract
Traditional reliability models of the gear system either simplify the gear system unreasonably or are too complicated to be applicable for some gear systems. Therefore, a reliability model of the gear system considering statistically dependent failure based on the theory of the order statistics is developed in this paper. Firstly, the gear bending fatigue test of small sample is implemented to obtain the PSN curves that are fitted by using the least square method and the linear regression method of the statistical parameters under the test stress levels. Then, according to stressstrength interference theory and the theory of the order statistics, the reliability models considering the number of load cycles from a gear to gear system are built, which indicates clear modeling process of the gear system. Finally, the proposed reliability model of gear system is validated by using the Monte Carlo Simulation (MCS) and the system reliability under special load history is analyzed qualitatively.
1. Introduction
Gear is an important component of the mechanical transmission equipment. Its reliability has the great influence on the system service life. Therefore, scholars have done lots of research works to estimate the gear system reliability. For instance, the simplification that a tooth denotes a gear is widely used for reliability estimation [1–5]. Based on the linear damage theory, a longterm fatigue damage analysis method for the gear bending fatigue of the wind turbine is developed [6, 7]. By the introduction of correlation coefficient, the reliability models of the gear system considering failuredependence are developed [8–10]. Recently, Xie [11] considered the gear system as multiconfiguration system and built the reliability models considering the toothpair meshing sequence based on the stressstrength interference theory.
Though lots of reliability methods of the gear system have been proposed, they either considered the gear as a component or built too complicated models. The methods introducing the correlation coefficient that is calculated experientially may cause the great error. Therefore, building an efficient and applicable reliability model of the gear system is necessary.
The life distribution under the special stress level can be determined by the test. The loadweighted average model for the fatigue reliability calculation is developed based on the life distribution [12]. Thus, for the random constant amplitude load history, the reliability can be calculated using this model.
Based on this idea, the reliability model of gear system considering statistically dependent failure is developed according to the theory of the order statistics in this paper.
2. Test of Gear Bending Fatigue Life
The gear bending fatigue life is tested by the gear test rig of closed power flow (Figure 1(a)). In the equipment, the failure of gear bending fatigue is monitored by the vibration signal detection sensor. The tested gear is lubricated by the spraying oil below 60°C. The oil is cooled by the water cooling system. The error of cycle counting is between 0.1% and 0.1% [13]. The revolution speed is 2000 rpm. The gear parameters are listed in Table 1. The tested gear and failure mode are shown in Figures 1(b) and 1(c).

(a)
(b)
(c)
The permissible bending stress is calculated in accordance with ISO 63363:2006 (5) as follows [14]. where the factors can be found in ISO gear standards, =520 MPa is estimated by the GB/T 384097 [15], and =1.5 is recommended by GB/T 384093 for the reliability purpose, whose failure probability is less than 1E4.
Therefore, 700 MPa is used for the maximum stress level of the gear bending test. Since =520 MPa, 550 MPa is used for the minimum stress level. 40% of total stress range [13] (i.e., 60 MPa for this test) is recommended for the stress interval between the maximum stress level and secondary stress level. Then, 50 MPa and 40 MPa are used for the secondary and third stress intervals, respectively. So, 700 MPa, 640 MPa, 590 MPa, and 550 MPa are used for the four stress levels.
Not less than 5 samples are recommended for each stress level in GB/T 14230. Thus, the testing method of small sample is used for the test of gear bending fatigue and 5 samples are used for each stress level. The gear bending fatigue lives are shown in Figure 2.
3. Statistical Analysis of Test Data
The distribution function of test data is unknown. In the test of gear bending fatigue, the normal distribution, Lognormal distribution, and twoparameter Weibull distribution are recommended for the test of the life distribution [13] as follows. where denotes the tested data; P denotes the fatigue life distribution; and denote the mean and standard deviation of tested data, respectively; and denote the mean and standard deviation of logarithmic data, respectively; λ denotes the scale parameter; η denotes the shape parameter.
For the test of fit goodness, the weighted least square method is used. The linear regression for the three distributions can be expressed uniformly as follows. where and denote, respectively, the intercept and slop of the linear regression function, which are calculated by the weighted least square method (WLSM) [16, 17]. where is the ordinal of tested data from small to large; denotes the weight factor; , , and for the three distributions are given in Table 2, whose linear functions are given in Appendix A.

The failure probability is estimated by the empirical failure probability as follows [18].where denotes the number of tested lives under each stress level.
The correlation coefficient is used for the comparisons of fit goodness between the three distributions, which is expressed as follows.
By virtue of the least square method, the statistical parameters of the fatigue lives under each stress level are shown in Table 3.

From Table 3, the correlation coefficient associated with the lognormal distribution is better than other two distributions. When the minimum value of linear correlation coefficient satisfies the three distributions ( represents the high correlation), the lognormal distribution with maximum correlation coefficient has the priority to be selected [13].
The studies [19–21] illustrate that the linear relationship exists between the stress level and the mean/std of lognormal distribution with parameters and . Then, the statistical parameters of lognormal distribution can be fitted with 95% confidence level by the least square method, which are shown in Figure 3. E() and D() denote, respectively, the mean and variance of lognormal distribution, which can be expressed as
Then, the parameters , and probability density function (pdf) can be expressed as follows.
The parameters and of each stress level calculated by using (7) and (8) are listed in Table 4.

The logarithmic fatigue life under each stress level with the failure probability can be calculated by (9) [22]. Then, PSN curves can be fitted by the least square method, which are shown in Figure 4. where is the onesided tolerance coefficient with a failure probability P, a confidence level of 1α, and degrees of freedom υ [22].
Shown in Figure 4 are the median SN curve and the PSN curves of gear bending fatigue with the confidence level of 95%, failure probability 10% that is fitted by the proposed regression method, and Ref. [23] method that shares the fatigue lives under the different stress levels and fits the PSN curves of small sample based on the principle of sample polymerization (detailed method is shown in Appendix B). Both methods provide the different life dispersions under different stress levels, which is different from ISO 12107:2012 method that has the same slope. The relative error of the slopes for the PSN curves of C95P10 obtained by using the proposed method and Ref. [23] method is 2.08%. This illustrates that the proposed method is applicable to the fitting of PSN curves.
4. Gear System Reliability Model
4.1. Reliability of a Gear Tooth
During the operation of the aircraft engine or wind turbine, the gear is subjected to the torque load, which can be described using a variable amplitude load history. For any meshing toothpair, the load can be considered as random variable.
For a simple case, under constant amplitude load history , the reliability of a gear tooth, whose bending fatigue life is greater than the special life , can be expressed as where denotes the pdf of gear bending fatigue life under the stress level s, which is the function of torque load .
If the occurrence probability of load is (γ=1,2,…m, ), then, the reliability of a gear tooth bending fatigue life can be expressed as
If the torque load is a random variable and can be described as the pdf: , then can be expressed as where is the load discrete interval.
Then, (11) can be expressed as
Let tend to zero and tend to infinity; (13) can be expressed as
4.2. Reliability of the Gear System
For the gear with teeth, the strengths of its teeth are independent and identically distributed (iid) random variables [11]. With respect to a gear, a tooth failure can cause the failure of the gear. Therefore, a gear can be considered as a series system composed of its teeth. Then, the reliability of a gear with teeth using the traditional method can be expressed as where denotes the pdf of strength under the special life .
The scholars [24, 25] have proved that the probability that the fatigue life is less than the special life under the stress equal to the probability that the strength is less than under the special life (i.e., ). Then, (15) can be expressed as
It can be seen that (16) is the square of (14). It considers that the failure of the gear tooth is statistically independent. Then, the reliability of gear system with δ gears can be expressed as where e=1, 2, 3, …, δ denotes the Ids of gears.
Eq. (17) is established from the viewpoint of statistically independent failure. Actually, a gear failure depends on the weakest component (tooth) in the viewpoint of probability. Then, for the gear with teeth, according to the definition of order statistics, the distribution of the fatigue life equals the distribution of the minimum order statistics of its teeth lives. Therefore, for the samples (teeth), the pdf and reliability of the gear fatigue life can be expressed as (18) and (19), respectively.
Due to the gear system subjected to the common input torque load, the failures of these gears are statistically dependent. Then, the reliability of the gear system with δ gears can be expressed as
Eqs. (17) and (20) are two kinds of gear reliability models which are developed from two viewpoints. The proposed model (i.e., (20)) considers the statistically dependent failure between gears.
5. Reliability Analysis of a Planetary Gear System under Variable Amplitude Load History
The topology of a planetary gear system is shown in Figure 5, which includes three planetary gears. The parameters of all gears are listed in Table 5.

The stress of tooth root is calculated as [14] where the variables and factors can be found in the ISO gear standards [14, 15].
The cycles of each gear are different. The relationship between them can be described using the gear ratio. Thus, cycles of each gear can be determined as where , , , , , , and denote, respectively, the cycles of the sun gear, ring gear, planetary carrier, gear1, gear2, gear3, and gear4. denotes the ratio of the tooth numbers between the ring gear and sun gear. denotes the ratio of the tooth numbers between the gear1 and gear2. denotes the ratio of the tooth numbers between the gear3 and gear4.
In order to validate the reliability model, MCS (Monte Carlo Simulation) is used for the reliability analysis. In the MCS, each tooth is given randomly an SN curve from the PSN curves of gear fatigue life. The MCS procedure is shown in Figure 6. It is assumed that the torque history of the input shaft follows the normal distribution, T~N(800, 50^{2}). Two kinds of methods of (17) and (20) are used for the reliability estimation. Eq. (17) is denoted by the traditional method that denotes the reliability model of independent failure between the gears and (20) is denoted by the proposed method. The system reliability is shown in Figure 7.
From Figure 7, the reliability obtained by the proposed method is in good agreement with that of the MCS. The traditional method underestimates the reliability of the gear system. It considers that the failures of each tooth and gear are statistically independent, which may cause the incorrect results.
The proposed dynamic reliability model (i.e., (20)) is built based on the PSN curves and load distribution. Because the PSN curves are determined by the test, the load dispersions also have the influences on the system reliability, which will be analyzed in the next section. Two cases are used for the influence analyses.
Case 1. The torque load distributions with different means are listed in Table 6. The estimation results with different means are shown in Figure 8.

Case 2. The torque distributions with different standard deviations are listed in Table 7. The estimation results with different standard deviations are shown in Figure 9.

From Figure 8, the fatigue life decreases with the increase of the mean of the torque distribution, which reflects the correct fatigue failure law. From Figure 9, at the early stage of the fatigue life, the reliability of the gear system with a large dispersion of the torque load is low. At the later stage of the fatigue life, a large dispersion of the torque load leads to a high reliability of the gear system. When the mean value is constant, a large dispersion of the torque load leads to a large dispersion of the fatigue life, which increases the system reliability at the later stage.
6. Conclusions
The proposed dynamic reliability model is established based on the PSN curves rather than the initial strength of component, which avoids the difficulties of determining the initial strength and strength degradation function [26]. The failure relationship between the teeth of a gear is expressed by the order statistics, which simplifies the complicated modeling process considering the statistically dependent failure [11] based on the stress intensity interference theory. The main conclusions are as follows.
(1) To obtain the accurate PNS curves, the fitting method of PSN curves based on the LSM and regression functions of the statistical parameters is proposed. The fitting results obtained by the proposed method and Ref. [23] method show the good agreement.
(2) The modeling method from a component (tooth) to the system (gear system) is used for the reliability modeling of the gear system, which indicates clear modeling process of the gear system. According to the principle of the system failure starting from the weakest component and the common cause failure, the reliability model of the gear system is built based on the theory of the order statistics.
(3) Through comparison with MCS, the reliability estimated by the proposed model is in good agreement with that of the MCS. By using the proposed model, the influences of the load dispersion on the system reliability at different service life stages are analyzed.
Appendix
A.
Linear regression functions can be expressed as follows:
(1) Normal distribution:
(2) Lognormal distribution:
(3) Weibull distribution:
B.
Taking the tested fatigue lives with 5555 samples, for example (4 stress levels with 5 samples of each stress level are used for the test), the procedures of estimating and are as follows.(1)The fatigue lives tested under the four stress levels should range from 10^{4} to 10^{7}.(2)The logarithmic life mean under each stress level can be estimated by (B.1) and the median SN curve can be fitted by the least square method. Then, the mean value of the fatigue lives under each stress level can be calculated by the median SN curve: , , , . (3)Let standard deviation of stress level (α=0.001, 0.002, …) and let =+Δk (Δk is the step length, such as 1e6). According to each , select all possible values for the estimation of the standard deviation of other stress level by (B.2). Then, convert the fatigue lives of the other several stress levels into the equivalent fatigue life under the stress level by (B.3). Thus, all the equivalent fatigue lives can be obtained. (4)Calculate the equivalent of all fatigue lives under stress level. Calculate the relative error between the equivalent and by (B.4). If (u=0,1,2,3…) satisfies (B.4), the right and can be determined. Otherwise, =+Δk (Δk is the step length, such as 1e6), return to the step and recalculate , until the right and are obtained. Then, the PSN curves can be fitted by the least square method based on the obtained parameters.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors acknowledge (1) Natural Science Foundation of China “Study on Reliability Modeling and Simulation of Mechanical Structure of Aero Engine” (Grant No. 51335003) and (2) The Collaborative Innovation Center of Major Machine Manufacturing in Liaoning. The first author wishes to thank other authors for the help of the test and discussions on the reliability modeling of the gear system.
References
 M. Ognjanović, M. Ristić, and P. Živković, “Reliability for design of planetary gear drive units,” Meccanica, vol. 49, no. 4, pp. 829–841, 2014. View at: Publisher Site  Google Scholar
 W. Dong, Y. Xing, T. Moan, and Z. Gao, “Time domainbased gear contact fatigue analysis of a wind turbine drivetrain under dynamic conditions,” International Journal of Fatigue, vol. 48, pp. 133–146, 2013. View at: Publisher Site  Google Scholar
 Q. J. Yang, “Fatigue test and reliability design of gears,” International Journal of Fatigue, vol. 18, no. 3, pp. 171–177, 1996. View at: Publisher Site  Google Scholar
 E.S. Aziz and C. Chassapis, “Comparative analysis of toothroot strength using stressstrength interference (SSI) theory with FEMbased verification,” International Journal on Interactive Design and Manufacturing, vol. 8, no. 3, pp. 159–170, 2014. View at: Publisher Site  Google Scholar
 D. Qin, Z. Zhou, J. Yang, and H. Chen, “Timedependent reliability analysis of gear transmission system of wind turbine under stochastic wind load,” Journal of Mechanical Engineering, vol. 48, no. 3, pp. 1–8, 2012 (Chinese). View at: Publisher Site  Google Scholar
 A. R. Nejad, Z. Gao, and T. Moan, “Longterm analysis of gear loads in fixed offshore wind turbines considering ultimate operational loadings,” Energy Procedia, vol. 35, no. 41, pp. 187–197, 2013. View at: Google Scholar
 A. R. Nejad, Z. Gao, and T. Moan, “On longterm fatigue damage and reliability analysis of gears under wind loads in offshore wind turbine drivetrains,” International Journal of Fatigue, vol. 61, pp. 116–128, 2014. View at: Publisher Site  Google Scholar
 P. Mou, F. Tao, C. Jia, and W. Ma, “A copulabased function model in fuzzy reliability analysis on the planetary steering gear,” in Proceedings of the International Conference on Quality, Reliability, Risk, Maintenance, and Safety Engineering (QR2MSE '13), pp. 375–378, July 2013. View at: Google Scholar
 Y.H. Liao, D.T. Qin, and C.Z. Liu, “Timedependent reliability analysis for a herringbone planetary gear set with failure dependency under random loads,” Journal of Vibroengineering, vol. 17, no. 6, pp. 2933–2946, 2015. View at: Google Scholar
 W. H. Chen, C. P. Zheng, and L. I. QiZhi, “Gear reliability analysis of 2.5MW wind turbine gearbox based on Copula function,” Chinese Journal of Engineering Design, 2015. View at: Google Scholar
 L. Xie, N. Wu, and W. Qian, “Time domain series system definition and gear set reliability modeling,” Reliability Engineering & System Safety, vol. 155, pp. 97–104, 2016. View at: Publisher Site  Google Scholar
 L. Xie, “A loadweighted average model for fatigue reliability calculation under stochastic constant amplitude cyclic load,” in Proceedings of the 12th ISSAT International Conference on Reliability and Quality in Design, pp. 274–278, August 2006. View at: Google Scholar
 “GB/T 1423093: Standard of test method for bending load capacity of gears”. View at: Google Scholar
 “ISO 63363: Calculation of load capacity of spur and helical gears  Part 3: Calculation of tooth bending strength,” International Organization for Standardization. Geneva. 2006. View at: Google Scholar
 “ISO 63362: Calculation of load capacity of spur and helical gears  Part 3: Calculation of surface durability (pitting),” International Organization for Standardization. Geneva. 2006. View at: Google Scholar
 K. E. Ahmad, “Modified weighted leastsquares estimators for the threeparameter Weibull distribution,” Applied Mathematics Letters, vol. 7, no. 5, pp. 53–56, 1994. View at: Publisher Site  Google Scholar  MathSciNet
 B. Bergman, “Estimation of Weibull parameters using a weight function,” Journal of Materials Science Letters, vol. 5, no. 6, pp. 611–614, 1986. View at: Publisher Site  Google Scholar
 Y.X. Zhao, Q. Gao, and J.N. Wang, “An approach for determining an appropriate assumed distribution of fatigue life under limited data,” Reliability Engineering & System Safety, vol. 67, no. 1, pp. 1–7, 2000. View at: Publisher Site  Google Scholar
 P. H. Wirsching, “Statistical summaries of fatigue data for design purposes,” NASA Contractor Reports, 1983. View at: Google Scholar
 P. H. Wirsching and Y.T. Wu, “Probabilistic and statistical methods of fatigue analysis and design,” Pressure Vessel and Piping Technology, pp. 793–819, 1985. View at: Google Scholar
 J. A. Collins, Failure of Materials in Mechanical Design, John Wiley & Sons, New York, NY,USA, 1981.
 “Metallic materials – Fatigue testing – statistical planning and analysis of data,” ISO 12107:2012. View at: Google Scholar
 L. Xie, “Principle of sample polymerization and method of PSN curve fitting,” Journal of Mechanical Engineering, vol. 49, no. 15, pp. 96–104, 2013. View at: Google Scholar
 Y. weiXing, Fatigue Life Prediction of Structures , National Defence Industry Press, Beijing, China, 2003 (Chinese).
 F. HuiMin, “Probability density function of fatigue strength,” Beijing University of Aeronautics and Astronautics, pp. 31–40, 1986 (Chinese). View at: Google Scholar
 P. Gao, S. Yan, L. Xie, and J. Wu, “Dynamic reliability analysis of mechanical components based on equivalent strength degradation paths,” Journal of Mechanical Engineering, vol. 59, no. 6, pp. 387–399, 2013. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2018 Enjun Bai et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.