Estimation on Reliability Models of Bearing Failure Data

1Mechatronics Engineering College, Henan University of Science and Technology, Luoyang 471003, China 2Collaborative Innovation Center of Machinery Equipment Advanced Manufacturing of Henan Province, Henan University of Science and Technology, Luoyang 471003, China 3National Center for Material Service Safety, University of Science and Technology Beijing, Beijing 10083, China 4Luoyang Bearing Research Institute Co. Ltd., Luoyang 471039, China


Introduction
In machinery products and engineering projects, bearings are the joints and wearing parts in the whole transmission system.Their operational reliability is the basis to establish optimization and improvement strategies and implement failure factor analysis, which directly relates to the operation security of product during service time.In bearing reliability estimation, the selection of failure distribution model is of great importance, because it directly relates to the precision of reliability prediction and has a huge influence on the usability of bearing.If the predicted reliability value is too high and the product performance exceeds a fatigue limit of normal operation, particularly for aerospace, high-speed rail, nuclear reactor, precision meter, and such systems, it will result in major vicious accident or even affect the national security [1,2].If the predicted reliability value is too low, the product function cannot be fully exploited, making the product lose its environmental adaption and leading to huge waste of conditional resources.Therefore, in order to guarantee the safe and stable operation of product system, it is essential to implement effective monitoring and diagnosing and precise reliability model estimation for bearings [3][4][5][6][7].
Reliability analysis is aimed at searching for the failure distribution information which can exactly reflect that the failure mechanism of product components accords with the analysis results of failure data.After fitting the fault or failure data into certain distribution form, the reliability estimation and prediction will be carried out, in which the distribution function of product failure time is the basis to study reliability.Since the failure state of bearings may be affected by different operation conditions, such as structural composition, material, load, lubrication, and numerous uncertainty factors, the failure life of actual situation is random, accompanied by multiple failure modes among which each mode can be mutually affected, acted, and dynamically varied [8].
There are some difficulties how to quickly and effectively utilize failure data for precise reliability analysis and model selection of production information.Recently, there is a lack of exact plan to describe their gradual change process during operation, and no perfect theoretical system is formed yet [9].Traditional reliability estimation theory is established on the basis of a large number of failure data.However, in many cases, the probability distribution of problems met in engineering and experiment are not normal ones.For example, the aircraft bearing, due to its high cost, very few failure data, and extremely high requirements for precision and reliability, has a harsh demand for test equipment.This makes it difficult to implement large sample life test or obtain failure data within limited testing time.So, when the product's life probability distribution is unknown, and we only have small sample data for reference, it is impossible for us to use existing reliability theory to accurately describe its failure evolving law.
At present, the researches of bearing failure data are mostly about Weibull distribution, log-normal distribution, gamma distribution, and binomial distribution reliability analysis methods.In particular, the Weibull distribution and log-normal distribution are widely used in reliability theoretical analysis.Though this has achieved certain results, they show large error and low precision in the process of product reliability estimation [10][11][12][13][14][15].The accuracy of reliability models increasingly got the attention of scholars and experts.Rodriguez-Picon et al. [16] considered a gamma process to marginally model the degradation of a performance characteristic through two degradation test phases performed sequentially and obtained a robust model to get reliability estimates considering the effect of two serial degradation tests.Reuben et al. [17,18] proposed a reliability evaluation method by using Weibull equation and made reliability estimation on product failure data of gearbox bearing and ceramic material, respectively, where results showed very small discrepancy between its fitting curve and the points of failure data.Through the reliability modeling disposing of failure process of large-scale and complicated machinery equipment, Pulcini [19] declared that its failure strength is not so monotonous, and on this basis, he proposed the reliability analysis model of nonhomogeneous Poisson process.To solve the real-time online reliability problems, Hong and Meeker [20] proposed an intelligent reliability estimation method based on dynamic state information changes, which brought much convenience to timely judge the dynamic running state of workpiece.Khaleghei and Makis [21] proposed a new competing risk model to calculate the conditional mean residual life and conditional reliability function of a system subject to two dependent failure modes, namely, degradation failure and catastrophic failure.In order to ensure that the classifier can correctly inspect the system failure information, Hwang and Lee [22] presented a new approach to overcome class imbalance problem and human factor influence by using classification technique, thus speedily and effectively implementing system reliability estimation.Zhang et al. [23] applied ANSYS/PDS module to make simulated analysis on the reliability of agricultural machinery chassis drive axle housing and probed into the influence of random variables such as geometric dimension, load, and material strength on drive axle housing.By exploring discrete random variable and analyzing the expectation interval and information capacity of certain entropy, Aviyente et al. [24,25] successfully solved the time frequency distribution problem and interval forecast problem of entropy.Xia [26] proposed a grey bootstrap method based on poor information theory, which conducted a reliability analysis of zero-failure data when the probability distribution information is known or unknown in life test, thus providing a strong theoretical reference to the reliability of poor information of zero-failure data.
Based on this, lots of topics and literature sources are mentioned associated with bearing capacity, distribution types, reliability, lifetime, and so on, because there is a very close relationship between them.Firstly, bearing capacity, lubrication condition, rotational speed, and other working conditions are important determinants affecting the bearing lifetime, and the set of the same batch bearing lifetime makes up a number of failure data under the above working conditions.Secondly, the distribution types of a number of failure data can be obtained according to statistical theory, and then its probability density function can be acquired easily.As we all know, the probability density function is the hub of data analysis and solution, and then, according to the probability density function and the given integral interval, the failure probability of bearings can be obtained during their service.Finally, using the unit one to subtract the failure probability, the reliability of bearing failure data is acquired.Therefore, these topics on bearing capacity, distribution types, reliability, and lifetime have a very close coherent interlocking, and all of them have an evidently direct or indirect relationship with the calculation of the reliability of bearing failure data.
This article used the failure data obtained from simulated test and bearing life failure test and made comparative analysis via log-normal distribution, Weibull distribution, and improved maximum entropy distribution, so as to select the optimum and precise reliability analysis model.First, the reliability empirical value calculated by Johnson [27] method was taken as standard.In specific analysis, twoparameter log-normal distribution was compared with twoparameter Weibull distribution, and three-parameter lognormal distribution was compared with three-parameter Weibull distribution, while in parameter estimation process of three-parameter log-normal distribution, the integral transformation moment method, linear moment method, and probability weighted moment method were used for comparative analysis, respectively.The research indicates that Weibull distribution is not applicable to all bearing failure conditions, and sometimes, the log-normal distribution has a smaller standard deviation and lower relative life error in reliability analysis.Then, the reliability estimation method for improved maximum entropy probability distribution was put forward to make reliability analysis on failure data, and this novel method has a high fitting degree and can be applied to all failure cases; when comparatively analyzing with Weibull distribution model and log-normal distribution model, both of the standard deviation and relative life error between its

Classical Reliability Empirical
Value.Suppose X 0 is a failure data series group of research object, and each failure data is unequal and nonredundant, which is denoted as where X 0 is the vector composed of such failure data group,  0 is the th failure data of this data series,  is the serial number of th failure data, and  is the number of failure data.
In case that the probability distribution or distribution parameter of failure data is unknown, the reliability of life failure data of research object can be nonparametric estimated using Johnson's median rank empirical value formula.The reliability empirical formula [28,29] of such method can be expressed by vector as where R0 refers to reliability empirical value vector.The formula to calculate reliability median rank empirical value is where  is the th failure data and  is the number of failure data.

Two-Parameter Log-Normal Distribution and Weibull
Distribution.Both distributions are common reliability models in engineering applications, especially the Weibull distribution which is widely used in analyzing bearing failure data and has achieved good research results.The probability density function of two-parameter lognormal distribution is The reliability function is where  is the random variable of life,  is the proportional parameter,  is the shape parameter,  > 0,  > 0, and  > 0.
The probability density function of two-parameter Weibull distribution is The reliability function is where  is the random variable of life,  is the proportional parameter,  is the shape parameter,  > 0,  > 0, and  > 0.

Parameter Estimation.
Maximum likelihood method [30,31] is widely used in the parameter estimation of all kinds of reliability models, which is one of the frequently used parameter estimation methods.For two-parameter lognormal distribution and two-parameter Weibull distribution, the maximum likelihood method is used for parameter estimation of these two models, respectively.When the maximum likelihood method is used to estimate two-parameter log-normal distribution, the likelihood equation set is obtained as below: where (,  2 ) is the likelihood function of log-normal distribution.
When the maximum likelihood method is used to estimate two-parameter Weibull distribution, the likelihood equation set is obtained as below: where (, ) is the likelihood function of Weibull distribution.Using the iterative method to solve the equation set, we can acquire the estimated values of two parameters in twoparameter Weibull distribution.

Three-Parameter Log-Normal Distribution and Weibull
Distribution.The probability density function of threeparameter log-normal distribution is The reliability function is where  is the random variable of life, (, ) is the parameter of log-normal distribution:  is the proportional parameter and  is the shape parameter, and  is the location parameter. > 0,  > 0, and  > 0.
The probability density function of three-parameter Weibull distribution is The reliability function is where t is the random variable of life,  is the proportional parameter,  is the shape parameter, and  is the location parameter. > 0,  > 0, and  > 0.

The Parameter Estimation of Three-Parameter Log-Normal Distribution.
In the process of parameter estimation of three-parameter log-normal distribution, the integral transformation moment method [32], linear moment method [33], and probability weighted moment method [34] were used for comparative analysis, respectively.These parameter estimation methods are mature and widely used in many fields, and the details are as follows.
(1) Integral Transformation Moment Method.In the following formulas,  is the location parameter, which can be determined by mean value , coefficient of variation  V , and coefficient of skew   , that is, 2 ) where with where Φ stands for standard normal distribution.
(3) Probability Weighted Moment Method where  and  are function of   , both of which cannot be expressed by explicit formulation;  V and   are coefficient of variation and coefficient of skew of variable , respectively. 0 ,  1 , and  2 are the probability weighted moments [35] of zero-order, first-order, and two-order, respectively. with ) , where V and  are transition variables to solve function value for  and   .

The Parameter Estimation of Three-Parameter Weibull
Distribution.The -order exceeding probability weighted moment [36] equation of Weibull distribution is where Γ is gamma function; for convenience,  is valued as  = 0, 1, 3, and we can obtain  1,0,0 ,  1,0,1 , and  1,0,3 , and the three parameters of three-parameter Weibull distribution are The exceeding probability weighted moments of observed sample are

Improved Maximum Entropy Reliability Model.
Improved maximum entropy method can make an optimal estimation with minimum subjective bias on unknown probability distribution.Firstly, according to reliability empirical formula, the reliability empirical vector R0 can be obtained for failure data.Secondly, using the empirical value of vector R0 to adversely deduce a frequency vector for discrete failure, a statistical histogram is acquired, which is convenient to be used to calculate Lagrangian multipliers, and it is different for the traditional maximum entropy to use amount of sample data to solve the Lagrangian multipliers.
Then, based on an internal mapped method, probability density function () for the improved maximum entropy can be obtained.Finally, reliability function for estimate true value is then acquired by integrating the function of ().

Discrete Failure Frequency
Vector.According to statistic theory, from the reliability empirical vector R0 in (2), we can get the discrete cumulative failure probability vector F 1 : where ( 0 ) is the reliability of the th failure data of initial data X 0 in (1) and  is the number of initial data.Suppose the corresponding discrete failure probability of each failure life data is  0 .For the first data, that is, when  = 1, let its failure probability be  01 =  1 .So, from the second data, that is, when  = 2, 3, . . ., , the corresponding failure probability of each failure life can be obtained by cumulatively subtracting the elements in vector F 1 successively as follows:  0 =   −  −1 ,  = 2, 3, . . ., .
So the discrete failure frequency vector of its failure life data is Let (30) correspond to statistic histogram, in which the abscissa is discrete failure life data  0 , and ordinate is the frequency   =  0(−1) ,  = 2, . . .,  + 1 that corresponds to class mid-value   =  0(−1) of each group.Normally, the histogram can be expanded to  + 2 group, that is,  = 1, 2, . . .,  + 2, and let Here, the processing of histogram is beneficial for utilizing Newton's method to solve Lagrangian multipliers in maximum entropy probability density function in the following.

Maximum Entropy Probability Distribution Density
Function.Suppose the probability distribution density function with maximum entropy is where  is the random variable of life;  is origin moment order, generally let  = 3∼8, and commonly  = 5;   is the th Lagrangian multiplier,  = 0, 1, . . ., , totally  + 1.
The first Lagrangian multiplier  0 is Other  Lagrangian multipliers shall satisfy Newton iteration method can be used to solve Lagrangian multiplier vector c.

Improved Maximum Entropy Probability Distribution
Numerical Solution.There are some difficulties to obtain the solution procedure of probability distribution by improved maximum entropy method.To achieve a quick numerical solution with good convergence, this article adopted the internal mapped Newton iteration method.First, the failure data series were mapped onto dimensionless interval [−, ],  = 2.718282.Then, the mapping data were sorted from small to large into Q-2 groups, and histogram was drawn, to obtain mid-value   and frequency   of each group.Later, the histogram is extended into Q groups; namely,  1 =   = 0,  = 1, 2, . . ., .
The value of  order origin moment   changes to The integration variable  turns into mapped variable , and integrating range  is mapped onto [−, ].The improved maximum entropy probability distribution density function changes to where  and  are interval mapping parameters.Integrate the improved maximum entropy distribution density function () in interval  = [ 01 ,  0 ].The obtained cumulative failure probability function  is Therefore, the improved maximum entropy reliability estimation truth function can be expressed as     distribution function curve, which accurately describes the distribution law of failure data.In four simulation examples, the standard deviations between empirical vector R0 and known distribution R 1 are all 0.0079, which indicates there is tiny standard deviation between the known reliability of failure data and that obtained by the empirical method.The research shows that median ranks empirical model can well estimate the reliability of product failure data with known probability distribution, without estimating the parameters, which avoids possible errors in parameter estimation and greatly improves the accuracy of estimation.So, in the following reliability model estimation, empirical vector will be taken as criterion to comparatively analyze the fitting degree of other models and empirical value, so as to determine whether reliability model is good or bad.

Computer Simulation Verification
Though this empirical formula does not need parameter estimation, is easy to use, and can precisely compute reliability, the estimation results of reliability truth are discrete, fluctuant, and uncertain, making it hard to conduct continuous estimation.Moreover, when failure data repeatedly appear, this formula cannot estimate the failure probability precisely.In order to get more accurate and continuous reliability estimation model, it is necessary to make further study on fitting curve of failure data (researched later).

Two-Parameter Model and Improved Maximum Entropy
Model Verification.If using two-parameter log-normal distribution and two-parameter Weibull distribution as reliability model, it is required to verify the accuracy of parameter estimation on two-parameter Weibull distribution and two-parameter log-normal distribution by maximum likelihood method.The improved maximum entropy method does not take function distribution into account, but using simulated data of two-parameter Weibull distribution and two-parameter log-normal distribution to verify the practicability of improved maximum entropy method.
First, generate a group of two-parameter log-normal distributed random numbers by computer, with parameter setting as (, ) = (4, 0.4),  = 30, and  is the number of simulated data.Use maximum likelihood method's equation ( 8) for twoparameter log-normal (TWPLN) to implement parameter estimation on these 30 random numbers: Use internal mapped method's equation (35) for improved maximum entropy (ME) to implement probability density estimation on these 30 random numbers.
Substitute the estimated parameters into ( 5) and (37) to obtain relevant reliability function, and substitute 1 into (3) to obtain empirical point R0.The results are shown in Figure 5.
From Figure 5, it is observed that reliability curve of twoparameter log-normal distribution conforms to the distribution of empirical points.Let significance level  = 0.05.We can get the -s test values of log-normal distribution fitting model and improved maximum entropy model which are 0.07 and 0.09, respectively, and the critical value is  = 0.2417, which shows that the test values are all less than critical values.Furthermore, this proves that using maximum likelihood method to make parameter estimation on twoparameter log-normal distribution can get better effect, and improved maximum entropy probability distribution model can be well applied in two-parameter log-normal distribution.
Generate a group of random numbers on two-parameter Weibull distribution by computer, with parameter setting as (, ) = (60, 2.5),  = 30.Use the reliability empirical equation (3) for empirical value (R0) to implement point estimation on these 30 random numbers.
Use maximum likelihood method's equation ( 9) for twoparameter Weibull (TWPW) to implement parameter estimation on these 30 random numbers: Use internal mapped method's equation (35) for improved maximum entropy (ME) to implement probability density estimation on these 30 random numbers.
Substitute the estimated parameters into ( 7) and (37) to obtain relevant reliability function, and substitute 2 into (3) to obtain empirical point R0.The results are shown in Figure 6.
From Figure 6, it is observed that reliability curve of twoparameter Weibull distribution conforms to the distribution of empirical points.Use -s test method, and let significance level  = 0.05.With estimated parameters of maximum likelihood method, we can get the -s test values of twoparameter Weibull distribution fitting model and improved maximum entropy fitting model which are 0.06 and 0.07, respectively, and the critical value is  = 0.2417, which shows that the test values are all less than the critical values.Furthermore, this proves that using maximum likelihood method to make parameter estimation on two-parameter Weibull distribution is feasible, and improved maximum entropy probability distribution model can be well applied in two-parameter Weibull distribution.
To sum up, from the results of parameter estimation on random numbers of two-parameter log-normal distribution and two-parameter Weibull distribution, it is known that taking maximum likelihood method as parameter estimation method is feasible for these two models with highly accurate estimation results; improved maximum entropy reliability function basically completely coincides with empirical value vector, which proves that this model is suitable for the above two distributions and has small error and high precision.

Three-Parameter Model and Improved Maximum
Entropy Model Verification
Three-parameter log-normal uses integral transformation moment (ITM) method's equation ( 14) for parameter estimation on this random numbers group: (48) Three-parameter log-normal uses probability weighted moment (PWM) method's equations ( 18), (19), and (20) for parameter estimation on this random numbers group: Improved maximum entropy (ME) uses internal mapped method's equation (35) for probability density estimation on this random numbers group.
Substitute the estimated parameters into (11) and (37) to obtain relevant reliability function, and substitute 3 into (3) to obtain empirical point R0, as shown in Figure 7.
It is observed in Figure 7 that the reliability curve basically conforms to the distribution of empirical points, and threeparameter estimation methods and improved maximum entropy probability distribution demonstrate perfect fitting degree.Use -s test method and let significance level  = 0.05 to make hypothesis testing on the results, as shown in Table 1.
In reliability function image, it can be found that three estimation methods have good fitting degree.And it is known from -s test that the -s test values of threeparameter estimation methods are less than critical values; thus, three methods are all suitable for three-parameter log-normal distribution.So, when we conduct parameter estimation of three-parameter log-normal distribution on test failure data, the above three estimation methods can be used for estimation.Meanwhile, -s test value of improved maximum entropy is also less than critical value, indicating that improved maximum entropy probability distribution model can be perfectly applied in threeparameter log-normal distribution with perfect estimation effect.

Three-Parameter Weibull Distribution and Improved
Maximum Entropy Distribution.In order to verify the feasibility of -order exceeding probability weighted moment based parameter estimation method to three-parameter Weibull distribution, let Weibull parameter (, , ) = (60, 2.5, 10), and simulate three-parameter Weibull distribution and life data by computer system ( = 25) as below:   Empirical value (R0) uses the reliability empirical equation (3) for point estimation on this group of random numbers.
(51) Improved maximum entropy (ME) uses internal mapped method's equation (35) for probability density estimation on this group of random numbers.
Substitute the estimated parameters into ( 13) and (37) to obtain relevant reliability function, and substitute 4 into (3) to obtain empirical point R0, as shown in Figure 8.
From Figure 8, it is observed that reliability curves of three-parameter Weibull distribution and improved maximum entropy probability distribution basically conform to the distribution of empirical points.Use -s test method, and let significance level  = 0.05.We can get the -s test value of three-parameter Weibull distribution model which is 0.01, the values of improved maximum entropy fitting model are 0.02, and the critical value is  = 0.2641.As the test values are all less than critical values, it is feasible to use -order exceeding probability weighted moment method in parameter estimation on three-parameter Weibull distribution, and improved maximum entropy probability distribution model can be well applied in three-parameter Weibull distribution.
To sum up, it can be known from the simulation results of random numbers obeying three-parameter Weibull distribution, three-parameter log-normal distribution, and improved maximum entropy probability distribution that the reliability estimation methods of the above three models are all feasible with highly accurate estimation results.Improved maximum entropy reliability function basically completely coincides with empirical value vector and is suitable for the above two distribution models with small error and high precision.Meanwhile, it also proves that the distribution of sample data can be ignored in improved maximum entropy application process, and its major feature is that it applies to the poor information issues with unknown probability distribution, trends, or prior information.This is because the improved maximum entropy model does not need parameter estimation or consider any distribution, but objectively processing experimental data.That is to say, there is no ideal model artificially presupposed before data processing, which overcomes the influence produced by subjective factor and parameter estimation error, so that the certainty rule in data change can be directly reflected.
At the same time, median rank empirical model can effectively assess the reliability of product failure data with high precision, so, in the following experimental research section, this empirical value should be taken as a criterion to decide whether the reliability model is good or bad; in twoparameter and three-parameter reliability model verification, multiple parameter estimation methods are feasible and effective with high precision and small error; therefore, they are safe and reliable in the following practical case applications; the novel improved maximum entropy model is suitable for all the above conditions, which has laid a foundation to search for quasi-ideal model in the following sections.

Experimental Research
4.1.Experimental Facility and Conditions.This bearing life reliability research adopted NTN bearing life experimental facility and material samples.The experimental facility is Φ12 point-contact life testing machine, and the material samples are cylindrical roller of Φ12 mm × 22 mm obtained by processing different batches of material under the same heat treatment conditions.The experiment was conducted at motor speed of 3900 r/min, indoor temperature of 26 ∘ C, and humidity of 53%.By applied load from compression spring, the stress imposed on roller and steel ball is 2.55 KN, and the experimental applied maximum contact stress is 5.88 GPa.In the timing from the start of experiment, when surface peeling occurs on contact region of cylindrical roller and steel ball, it will make cylindrical roller vibrate more.And if the vibration value reaches certain amplitude, the inductor switch will automatically jump up, the motor will stop to run, and the experiment ends.At this time, read the test run time of testing machine, which is the failure life of this roller material.

Test Data.
There are 3 groups of failure data test in total.The failure time recorded by experimental facility is initial data with unit in minute.For easy calculation, failure data is conversed into data with unit in hour and sorted from small to large into a group of vector. In Substitute failure data 1 into improved maximum entropy (ME) model's equation (35) for probability density estimation.
Their reliability curves are shown in Figure 9.
Substitute failure data 2 into improved maximum entropy model for probability density estimation.
Their reliability curves are shown in Figure 10.
Substitute failure data 3 into improved maximum entropy model for probability density estimation.
Their reliability curves are shown in Figure 11.
From Figures 9 to 11, it can be discovered that three reliability curves basically conform to distribution of empirical points and show good fitting degree.Use -s test method and let significance level  = 0.05, to make hypothesis testing of parameter estimation results; the results are shown in Table 2.The research indicates that -s test values of three failure data distribution models in each example are all less than the critical values, so three curves can describe the distribution rule of failure data.
Calculate the standard deviation of reliability function value and empirical points and substitute failure data into empirical equation (3), two-parameter log-normal distribution reliability equation ( 5), two-parameter Weibull distribution reliability equation (7), and improved maximum entropy probability distribution reliability equation (37).Then we can obtain reliability empirical value vector R0 and reliability estimated truth vectors R 1 , R 2 , and R 3 and substitute them into standard deviation equation: where  is the number of failure data of each group and  = 1, 2, 3.
According to (61), we can get the standard deviation when each failure data group takes two-parameter log-normal (TWPLN) distribution, two-parameter Weibull (TWPW) distribution, and improved maximum entropy (ME) probability distribution as reliability models.At the same time, figure out the life value of three distributions when their life failure probability is 10% and 50%, respectively.That is to say, when reliability function () = 0.9, the value of life  is 1, and when () = 0.5,  value is 2.Comparing the life values under these two failure probabilities, and integrating standard deviation, we can determine which reliability model has small error and high precision; the results are shown in Table 2.This is the main problem with this article to be revealed; the difference of the standard deviation for each model may be small, but the difference of life value  under life reliability of 90% and 50% may be extremely large.Thus, an accurate reliability model of bearing failure data is vital and necessary.
It can be known from Table 2 that the -s test values of the three reliability models are all less than critical values, and all models conform to these failure data groups and thus can be taken as their reliability model.There is very small difference between the calculated standard deviations of each data group by three models, all at 0.01 orders of magnitude (except 0.1008), which indicates that three models have high precision in reliability estimation.However, in comparison with P1 or P2 value of each example, it is not difficult to discover that there is large difference between P1 and P2 values for different distribution models.That is also the greatest distinction in numerical solution by three models.In other words, under equivalent reliability, different estimation model will cause different bearing failure life.
In addition, the life value t under life reliability of 90% and 50% is an important indicator to study reliability.We took the distribution of small standard deviation, that is, the improved maximum entropy probability distribution, as datum to calculate the relative error of other two distributions under failure probability of 10% and 50%.When reliability is 90%, the relative error of 1 log-normal distribution is Similarly, the relative errors of each data group under different distributions can be obtained as shown in Table 3.
From Table 3, though there is very small difference between the standard deviations of three models, under reliability of 90% and 50%, their life values can be far different and the maximum relative error reaches 314.62%.Hence, the selection of distribution model will directly affect the precision of predicted reliability, which should not be blindly chosen in practical applications.That is to say, in reliability function image, provided a small change in ordinate value, it may bring very large change in abscissa value.To reduce life estimation error, this article judged the fitting degree of reliability model by whether standard deviation is large or small.
In the table, the standard deviations of reliability estimation truth-value vector and reliability empirical value vector for 1 and 3 two-parameter log-normal distribution groups are less than those of two-parameter Weibull distribution.But 2 groups are the opposite.This illustrates that on premise of that both distributions satisfy this failure data group, sometimes, the two-parameter log-normal distribution being taken as a reliability model may have higher fitting and smaller error than two-parameter Weibull distribution.So, in daily reliability analysis on bearing life, people should not merely use two-parameter Weibull distribution as reliability model for analysis, in order to prevent large error in bearing reliability predication and the occurrence of vicious accident.Regarding the above three examples, all the improved maximum entropy curves participate in the fitting, and fitting curves almost coincide with empirical values.What is more important is that the standard deviations of this model are all the smallest in reliability estimation on 1, 2, and 3 failure data groups.In summary, compared to two-parameter Three failure data groups were estimated by threeparameter log-normal distribution, three-parameter Weibull distribution, and improved maximum entropy probability distribution.Their reliability images are shown in Figures 12-14, in which, three-parameter log-normal distribution used integral transformation moment, linear moment, and probability weighted moment methods for parameter estimation, three-parameter Weibull distribution used -order exceeding probability weighted moment methods for parameter estimation, and the novel improved maximum entropy directly made probability fitting operation without considering failure data distribution condition.Then, the -s test method is used with significance level  = 0.05 to implement hypothesis testing on each method and results; finally, it figured out standard deviations according to (61), so as to easily analyze and compare the precision of reliability models; the results are shown in Tables 4-6.
In Tables 4-6, it can be known from hypothesis test results that, in three methods on parameter estimation of three-parameter log-normal distribution, only -s test values of integral transformation moment method are less than critical values: namely, only integral transformation moment method is suitable for three-parameter log-normal distribution model.But such "Suitable" only satisfies hypothesis test of data verification; in practical applying, the location parameter  of three-parameter log-normal distribution means the minimum failure data; that is,  should be larger than zero.While the location parameter  obtained by using integral transformation moment method is all less than zero in three groups of failure life tests, which goes against the practical significance of this formula in life test application, so integral transformation moment method does not apply to these failure data groups.That is to say, none of the above three-parameter estimations can be taken as an estimation approach for three-parameter log-normal distribution in failure data reliability model application.As a result, at present stage, it is considered that three-parameter lognormal distribution cannot be perfectly used in reliability estimation of bearing life, and its parameter estimation approach needs further exploration.
The -s test values of three-parameter Weibull distribution are less than the critical values with good fitting; namely, -order exceeding probability weighted moment method performs well when it is used for parameter estimation of three-parameter Weibull distribution.But for 2 failure data group, the location parameter of three-parameter Weibull distribution is  = −1.0515< 0, which goes against real meaning of life.This also indicates that three-parameter Weibull distribution performs better on bearing life failure than three-parameter log-normal distribution does, but it does not be the same with all situations.However, in engineering practices, it is widely believed that three-parameter Weibull distribution can better describe the distribution rule of bearing performance failure data, because it has more parameters and comprehensive information mining so that we can consider the general characteristics of research object from multiple aspects.But this research shows that threeparameter Weibull distribution has a large misunderstanding during its wide application in engineering practices, and it may not accurately recognize product's performance reliability in all cases.
Improved maximum entropy reliability curves participate in fitting in three examples, and the fitting curves almost coincide with empirical values.Its standard deviation is the minimum in three groups of life test compared to previous two models, declaring its highest fitting degree.It is observed in Figures 12-14 that the novel improved maximum entropy reliability curves trend to be deck chair shape, directly showing the change trend of initial failure data.The obtained curve does not require parameter estimation, thus preventing calculation error that may be brought by traditional mathematical hypothesis definition and condition limitation.In summary, relative to three-parameter log-normal distribution and Weibull distribution, the improved maximum entropy probability distribution can make reliability estimation with the highest precision and minimum error and directly reflect the general change rule of failure data.No matter there is truly distribution model (quasi-ideal distribution model) to calculate bearing failure data reliability or not, at least, the novel improved maximum entropy probability distribution is a precise simulation of quasi-ideal failure data distribution model.

Improved Maximum Entropy Reliability Test. Through calculation in test examples
, it can be obtained that the reliability estimation method of improved maximum entropy probability distribution is suitable for all situations, and it is observed that improved maximum entropy reliability model estimation approach has a perfect effect on reliability estimation on bearing failure data.In reliability model estimation on three failure data groups, the two-parameter log-normal distribution of 1 and 3 is superior to two-parameter Weibull distribution, and the two-parameter log-normal distribution of 2 is inferior to two-parameter Weibull distribution; the three-parameter Weibull distribution of 1 and 3 is superior to three-parameter log-normal distribution, and both distributions of 2 group are not satisfactory.So we compared the two-parameter log-normal (TWPLN) distribution of 1 and 3, two-parameter Weibull (TWPW) distribution of 2, and three-parameter Weibull (THPW) distribution of 1 and 3 with improved maximum entropy (ME) method, respectively, and the results are shown in Tables 7 and 8.
The results display that, in comparison with three reliability models, the standard deviation between estimated truthvalue of improved maximum entropy reliability model and empirical value vector is the smallest.Taking the improved maximum entropy model as datum, other two models possess a larger relative life error under failure probability of 10% and 50%, with maximum relative error of 314.58%, which means the log-normal distribution and Weibull distribution possess low precision and large error in reliability prediction of bearing failure data.Once again, this proves that the novel improved maximum entropy probability distribution being taken as reliability model estimation method possesses the best effect and lowest error in estimating bearing failure life data.
In order to verify that the novel proposed model of maximum entropy distribution can be applied to poor information problem with small sample and unknown probability distribution, another 4 group test was carried out, with data number  = 5.
Initial data series are 4 = [4.36 14.68 With the help of improved maximum entropy method, the reliability estimation results of this bearing failure data group are shown in Figure 15.The reliability empirical value vector of this example is R0 = [0.87040.6852 0.5000 0.3148 0.1296], and the reliability estimation truth-value vector obtained by improved maximum entropy method is R3 = [0.94250.7509 0.5486 0.2800 0.0895].
Depending on reliability empirical value vector R0 and reliability estimation truth-value vector R3 of improved maximum entropy method, it can be figured out that the standard deviation between R3 and R0 is very small at 0.0542, and the maximum difference value between both reliabilities is only 0.0721.The research results show that it is effective and feasible to use the reliability estimation truth-value function obtained by improved maximum entropy method to assess the reliability of small sample failure data with unknown probability distribution.In the condition of having failure data but without probability distribution or any prior information, the improved maximum entropy method can perfectly estimate the reliability function.And, at the same given life, the difference value between the reliability obtained by improved maximum entropy method and that obtained by empirical value is very small.This method also remedies the defect of existing method to only solve reliability estimation with known probability distribution.
To sum up, in simulation test, three reliability models display good fitting.In other words, in theoretical state, log-normal distribution, Weibull distribution, and improved maximum entropy model can be applied in reliability analysis on product performance failure problems.But in experiment part, we made comparative analysis of three reliability estimation methods according to actual bearing failure data, with the results showing that both standard deviation and relative error are the smallest between reliability empirical value vector R0 and reliability estimation truth-value vector R3 of improved maximum entropy method.Besides, this method can solve the poor information issue with small sample and unknown distribution that cannot be assessed by classical statistics.Therefore, based on all failure data series in this article, maximum reliability estimation method is the optimal reliability estimation approach due to its best fitting and highest precision.

Conclusions
Though the reliability estimation of bearing performance failure data can be realized by the above models, their prediction precision is far different.In practical application, we cannot rush to a conclusion by reliability obtained from single model.
The standard deviations of two-parameter log-normal distribution for 1 and 3 failure data groups are 0.0802 and 0.0950, respectively, whose performance is superior to two-parameter Weibull distribution.Although the threeparameter Weibull distribution taken as reliability model in example test of 1 and 3 groups performs better than threeparameter log-normal distribution, it cannot be applied to all analysis cases.So, in daily reliability analysis of bearing life, we should not only use Weibull distribution as reliability model for analysis.
For a novel improved maximum entropy model compared to Weibull distribution and log-normal distribution, the relative life error and standard deviation of its truthvalue vector and empirical value vector are the smallest, in which, the maximum relative error of log-normal distribution reaches 314.58%, and Weibull distribution reaches 189.37%.
Whether there are really other distribution models (quasi-ideal distribution model) to calculate bearing failure performance reliability or not, at least, the improved maximum entropy probability distribution is a precise simulation to such quasi-ideal distribution model.
The novel proposed model of maximum entropy probability distribution does not take parameter distribution into account and allows poor information issues with unknown probability distribution, unknown prior information, or trends.This provides important theoretical reference to many uncertain information and poor information issues in engineering and even aerospace field and remedies the deficiency of classic statistics.

Figure 1 :Figure 2 :
Figure 1: Verification of empirical model by two-parameter Weibull model.

Figure 3 :Figure 4 :
Figure 3: Verification of empirical model by two-parameter log-normal model.

Figure 5 :
Figure 5: Reliability model simulation of two-parameter log-normal distribution and improved maximum entropy distribution.

Figure 6 :
Figure 6: Reliability model simulation of two-parameter Weibull distribution and improved maximum entropy distribution.

Figure 7 :
Figure 7: Reliability model simulation of three-parameter log-normal distribution and improved maximum entropy model.

Figure 8 :
Figure 8: Reliability simulation of three-parameter Weibull distribution and improved maximum entropy model.

Figure 15 :
Figure 15: Improved maximum entropy reliability function image of small sample data 4.

Table 1 :
Comparison of three-parameter log-normal distribution and improved maximum entropy distribution.

Table 2 :
Comparative results of three reliability models of failure data.

Table 3 :
Relative errors of each group data under different distributions.

Table 6 :
Comparison results of failure data 3 by three models.

Table 7 :
Compare the standard deviations of three groups of experimental feasible reliability models.

Table 8 :
Compare the relative life errors of three groups of experimental feasible models.