MPE Mathematical Problems in Engineering 1563-5147 1024-123X Hindawi 10.1155/2019/8947905 8947905 Research Article Reliability Assessment for Very Few Failure Data and Weibull Distribution https://orcid.org/0000-0001-6875-1538 Zhang Lulu Jin Guang You Yang Yang Mijia College of Systems Engineering National University of Defense Technology Changsha 410073 China nudt.edu.cn 2019 2592019 2019 03 07 2019 16 08 2019 05 09 2019 2592019 2019 Copyright © 2019 Lulu Zhang 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.

Only very few failure data can be obtained for the time censored test of high-reliability and long-life products. For very few failure data, the current methods fail to obtain both the point estimation and confidence interval of reliability parameters. If the point estimation and confidence interval of reliability parameters are obtained based on different methods, the results tend to be unreliable. In this study, based on the existing research, a Bayesian reliability evaluation method for very few failure data under the Weibull distribution was proposed. First, the range of failure probability was limited based on the convexity and self-features of the Weibull distribution function. Second, based on the background of the sample with very few failure data, the pretest distribution function and parameters were set and solved. The point estimation and confidence interval model of failure probability based on the Bayesian formula was established. The improved match distribution curve method was used to compute both the point estimation and confidence interval of reliability parameters. Furthermore, by comparing the results of numerical examples, the calculation results obtained by the proposed method were verified as being very reasonable. Finally, taking wet friction plates as an example, the results showed the effectiveness of this method in engineering practice.

National Natural Science Foundation of China 71371183 71071158
1. Introduction

In product design, research and development, production, improvement, and other aspects, it is usually necessary to verify the reliability of the product through reliability tests. For current products, the reliability is high and the life is longer and longer. The failure data collected by the reliability test have become less and less, and even sometimes the failure data cannot be collected. How to study the reliability of the product under the condition of very few failure data has become an urgent problem that needs to be solved.

The Bayesian method is an effective method to deal with the reliability problem of products with few failure information. It has been widely used in reliability parameter point estimation [1, 2]. Martz and Waller  studied the reliability evaluation problem under the condition of no failure data for the first time. After decades of development, research on this issue has achieved a lot of findings [4, 5]. Kurz et al.  studied the early life failure of a product and established the advanced Bayesian estimation models under the Weibull distribution. Jiang et al.  based on the different values of the shape parameters in the Weibull distribution found that the value of the probability of failure was determined by the characteristics of the concavity and distribution function, and the point estimation value of the failure probability was calculated. Yi et al.  proposed an E-Bayesian reliability evaluation method based on the condition of exponential distribution and zero-failure data, and the method was applied to the analysis of the seekers failure rate and reliability. In order to solve the problem that the Weibull parameter prior distribution is difficult to determine, Kan et al.  introduced the expert judgment process and established the multiple numerical control machine tools (NCMTs) Bayesian failure-free data reliability model and evaluation method. The above method studies the reliability parameter point estimation of products without failure data but did not discuss the confidence interval of reliability parameters.

With regard to the confidence interval problem of reliability parameters under the condition of no failure data, the optimal confidence limit method was proposed  and the point estimation of reliability parameters and the one-side confidence lower limit of parameters were obtained. Li et al.  proposed a new model to solve the consistency problem, which improved the accuracy of reliability interval estimation based on the reliability evaluation. The product reliability problem under unequal time censored failure data was studied . The time censored test was extended, and the solution method of the lower confidence limit under the Weibull distribution and exponential distribution was given. Xu and chen  proposed the method of the two-sided modified Bayesian (M-Bayesian) credible limit and used this method to study the failure probability and reliability confidence interval evaluation method under the condition of no failure data. It can be seen from the current research that the use of the distribution curve method is often limited to the point estimation of the parameters. It is rarely used in the parameter confidence interval solving process. In addition, the reliability of the product under the condition of no failure data has been studied more. However, it is rare to find a reliability assessment method for very few failure data conditions. The traditional reliability assessment method is generally based on obtaining more failure data, and the traditional method is often used in the case of very few failure data, resulting in a large deviation.

In the Weibull distribution, the match distribution curve method was improved for the characteristics of very few failure data, and the range of failure probability was limited according to the value of the shape parameters m and thus further improved the accuracy of the evaluation range. This method can analyze the reliability evaluation of products under the condition of few failure data, including the point estimation of parameters and the confidence interval estimation of parameters and fundamentally solve “disjointed” problems between parameter point estimation and parameter interval estimation.

2. Weibull Distribution

In the reliability evaluation of products, the first step is to assume the probability distribution form of product life compliance, on the basis of which the sample information is used to evaluate the distribution parameters, and further calculate the reliability indicators. Life distribution is the basis of reliability research. After years of development, the commonly used life distributions include exponential distribution, gamma distribution, Weibull distribution, and lognormal distribution, where the Weibull distribution is derived from the weakest link model. This model is like a chain in which many rings are connected in series. When both ends are subjected to tension, any one of the rings can break, the chain will fail, and the chain failure will occur in the most vulnerable link. Therefore, the life of units, components, devices, and equipment which are not working properly due to partial failure can be regarded as approximating the Weibull distribution. Most electronic products, mechanical products, and electrical products (such as bearings, generators, hydraulic pumps, and materials) are subjected to the Weibull distribution .

The Weibull distribution is defined as(1)ft=mηtηm1exptηm,t>0,where TWm,η, m is shape parameters, and η is scale parameters.

The cumulative distribution function (CDF) of Weibull distribution is defined as(2)Ft=PTt=1exptηm.

For Rt=1Ft,(3)Rt=exptηm.

Combined with the reliability formula λt=ft/Rt, the failure rate function under the Weibull distribution is λt=m/ηt/ηm1. When the shape parameter m<1, the failure probability density function ft and the failure rate function λt are both decreasing functions, which describe sudden failure, which is equivalent to the early failure of the product. When m=1, the Weibull distribution is the exponential distribution. When m=2, the Weibull distribution is the Rayleigh distribution. When m>3, the Weibull failure probability density function has single-peak symmetry, which approximates a normal distribution and describes the product gradual failure. At the same time, it indicates the combination of sudden failure and gradual failure, in which m can be adjusted according to different failure types.

3. Model

In the reliability test of products with higher reliability, test results often show a situation in which very few test samples fail during the entire reliability test or even no sample failures, and the number of such failed samples is less than 10% (including when the number of failed samples is 0), and this is called the minimum data failure test. Reliability research work based on this is important for predicting the life of highly reliable products.

In the product life comprehensive test verification method, the truncation life test is divided into two categories, one is the fixed censored test and the other is the time censored test. Of them, the time censored test needs to collect more test data in order to evaluate the reliability of the product with time and is used for products with a short life cycle and a low test cost. In reality, the problem of collecting only very few failure data occurs mainly in the time censored test. This study evaluates product reliability based on the minimum failure data obtained from the time censored test.

Assume there are N samples in the time censored test, and a total of kkN time censored tests are performed. The end of each time censored test τii=1,2,,k is indicated. There are ni samples observed at each censored test time τii=1,2,,k and i=1kni=N. If there are failed samples in the period of τi1,τi, the failure time tijj=1,2,,ri is recorded corresponding to the number of failed samples ri. The samples are numbered after each time, from 1 to N. After all the time censored tests are completed, the life data tii=1,2,,N are collected corresponding to the sample number. The special case is that the failure data are not collected at this time, they are a special case where there are very little failure data.

In order to facilitate the subsequent derivation, let si=j=iknj, si indicates the number of samples still participating at time τii=1,2,,k in the time censored test plus the number of internal τi1,τi failed samples ri. In addition, the probability of failure at each truncation moment is recorded pi=Fτi=Ptτi,i=1,2,,k, and it is known that p1p2pk based on the test results.

This study discussed the solution method of Rt based on the match distribution curve method. First, it is necessary to estimate the failure probability pi of each truncation moment τi and fit the distribution curve based on τi,pi, and then further estimate other required distribution parameters, wherein the failure probability is the key to the method.

4. Failure Probability Estimation

For the estimation of pi, the classical estimation method used in the past  was pi=ri+0.5/si+1, if there is no sample failure within the time period τi1,τi, pi=0.5/si+1. This failure probability estimation method is too rough, and the calculation accuracy is poor. In order to overcome the problem of poor calculation accuracy of the classical method, pi is set as a random variable and the Bayesian method [16, 17] is used to estimate the poster distribution of pi(4)πpit1,,tN=πpiLt1,,tNpipiπpiLt1,,tNpidpi,where t1,,tN is the sample information of the samples, Lt1,,tNpi is the likelihood function of the sample, and πpi is the prior distribution. It can be seen from the above formula that the Bayesian method can fuse a variety of information, which can further improve the estimation accuracy of the failure probability pi. There are two key problems that need to be solved by the Bayesian formula. First, it is necessary to determine the range of values of the failure probability pi. Second, it is necessary to determine the distribution form (i.e., distribution function) and distribution parameters of the prior distribution.

4.1. Range of Failure Probability

The value range of the failure probability pi is li,ui, and generally the value range of the failure probability pi is assumed to be 0,1. This assumption of the failure probability is too conservative, and the calculation accuracy is poor. In order to improve the accuracy of the range of failure probability, for the characteristics of the failure probability range of the product with very few failure data, this study used the concavity and function of the distribution function to set the range of failure probability.

By taking the second derivative of the Weibull distribution function (2), we have(5)d2Ftdt2=mtm2expt/ηmηmm1mtmηm=m2η2tηm2m1mtmηmexptηm.

When m1, it is easy find from (5) that d2Ft/dt2<0 and the distribution function Ft is a convex function, and according to the nature of the convex function, we have(6)p1t1>p2t2>>pi1ti1>piti>>pktk.

Combining p1p2pk, the interval of pi is(7)pi1pi<titi1pi1,i2.

First, it is necessary to define the value range 0<p1<plim of p1, where plim is the upper limit of p1. The value setting of plim is based on the actual situation of the sample and expert experience. The value is generally not more than 0.5 in actual engineering.

When m>1, it is difficult to determine the concavity and convexity of the distribution function. Based on the function characteristics, (2) is transformed into the following equation:(8)lnln11p=mlntmlnη.

Putting pi1 and pi into (8), it results in the following equation:(9)lnln1pi1ln1pi=mlnti1ti<lnti1ti.

Then, (9) is transformed into(10)ln1pi1ln1pi=mti1ti>ti1ti.

Then, pi>11pi1ti/ti1, and the lower limit of pi is 11pi1ti/ti1. Let the upper limit of pi be pupp, and the value of pupp is based on the actual situation. Based on the distribution function and the function characteristics, the value interval of pi is defined as(11)lipi1,m1,11pi1ti/ti1,m>1,uititi1pi1,m1,pupp,m>1.

Since the exact value of the failure probability pi cannot be obtained during the calculation, the estimated value p^i is used to replace failure probability pi.

4.2. Prior Distribution Form and Distribution Parameters

From the sample data and historical data, the prior distribution is a nonuniform distribution, which is closer to the beta distribution. In addition, the beta distribution Beta (a, b) is also the conjugate prior distribution of the parameters, and the conjugate prior distribution can be easily integrated from the historical information. It provides a reasonable premise for the analysis of future experimental data, and some parameters of its posterior distribution can also be well explained ; this study used Beta (a, b) as the prior distribution. According to the change of failure probability value range, the value range of Beta (a, b) is adjusted accordingly:(12)πpia,b=pia11pib1Bia,b,where Bia,b=liuixa11xb1dx.

Since the reliability and data collected in the time censored test are very small, the failure probability pi is small and the possibility of a large value is small. This is consistent with the characteristics of the subtraction function, and the prior distribution density function πpia,b is derived:(13)dπpia,bdpi=pia21pib2Bia,bab1pipib10.

Equation (13) is only true when a0,1 and b1. According to the nature of the probability density of the Bates distribution function, when a0,1, the larger the b, the finer the tail of the probability density function. From the perspective of Bayesian estimation of robustness, a finer right tail of the prior distribution makes the Bayesian estimation less robust, so the value must have an upper limit. It is c. In practical applications, c should be selected in [2, 8], and the arbitrary value in the interval [2, 8] has little influence on the reliability, failure probability, and other indicators, which can be ignored . It is very difficult to further determine the exact values of a and b. Rather than using expert experience to determine a uniquely determined value, as it is not as uniform as the hypothesis and within the respective range of a and b, we have(14)πia=U0,1,πib=U0,c.

The prior distribution of pi is(15)πpi=1c01pia11pib1c1Ba,bdadb.

4.3. Bayesian Failure Probability Estimation 4.3.1. Failure Probability Estimation without Zero-Failure Data

The failure probability acquisition under the condition of zero-failure data was studied . In the time-censored test, there are si test samples at time τi. The likelihood function of pi is Lsipi=1pisi. Let 1pi2 denote he core of the prior distribution. To meet the requirement of prior distribution pA1pi2dp=1, the prior distribution function is πpi=31pi2/1pl31pu3. Based on the Bayesian formula (4), the poster distribution can be found:(16)πpisi=si+31pisi+21plsi+31pusi+3.

Under the square loss assumption, the expectation of pi is(17)p^i=si+31plsi+31pusi+3qpuqpl,qx=1xsi+4si+41xsi+3si+3.

4.3.2. Failure Probability Estimation Based on Very Few Failure Data

From the time censored test, there are si samples remaining in the test at censoring time τi, where ri samples have failed. The likelihood function of pi is(18)Lsi,ripi=Csiripiri1pisiri.

Putting (15) and (18) into Bayesian formula (4), we can obtain the poster distribution:(19)πpisi,ri=1c01pia+ri11pisiri+b1/Bia,bdadb1c01liuipia+ri11pisiri+b1/Bia,bdpidadb,where pili,ui is the value obtained from (11). Under the square loss assumption, the poster distribution expected value of Bayesian poster point estimation p^i is(20)p^i=liuipiπpisi,ridpi=1c01liuipia+ri1pisiri+b1/Bia,bdpidadb1c01liuipia+ri11pisiri+b1/Bia,bdpidadb,where the integral calculation can be performed by the numerical integration method.

5. Reliability Estimation Based on Very Few Failure Data

After obtaining the point estimation p^i of failure probability pi, the distribution curve method can be used to fit the points into a distribution curve. In this study, the weighted least squares estimation method was used to fit the distribution curve.

Let y=lnln1p, x=lnt, b=mlnη, and then (2) is transformed into(21)y=mxb.

The weighted least squares estimation method is used to fit xi,y^i, obtaining the point estimation of parameters m and η, and then the reliability is estimated. By minimizing the square of fitting errors, we have(22)Q=i=1kwiy^imxi+mlnη2,where wi=niti/j=1knjtj,k2.

Let A=i=1kwixi, B=i=1kwixi2, C=i=1kwiyi, and D=i=1kwixiyi, and with the derivative calculation of (22), we obtain(23)m^=DACBA2,η^=expADBCDAC.

The estimated distribution parameters obtained from (23) are put into (3), getting the estimation R^t based on the Weibull distribution:(24)R^t=Rt;m^,η^=exptη^m^.

6. Confidence Interval Estimation Based on Very Few Failure Data

For the Weibull confidence interval, Han  proposed the calculation method of the optimal confidence lower bound under the condition of no failure data and then obtained the lower confidence limit of the Weibull distribution under the confidence level 1α:(25)RZLt=exptmlnαi=1kniτim,where the shape parameter m is the real value, but this is difficult to obtain in the actual use process, and therefore, the point estimation value m^ is generally used instead, which will bring a certain error into the final calculation result. In addition, the confidence interval estimation is different from the method used for point estimation, and it is easy to have a “disjointed” problem. If the acquired information is used, the Weibull partitioning interval is obtained based on the distribution curve method, and the secondary error can be further avoided.

Under confidence level 1α, first the upper bound p^iu of the failure probability pi is estimated. Based on the poster distribution in (19), Gx,s,r is defined as(26)Gx,s,r=lx1c01pa+r11psr+b1/Bia,bdadb1c011upa+r11psr+b1/Bia,bdpdadbdp.

According to the definition of the upper confidence limit, we have(27)Gp^iu,si,ri=1α.

It is known that Gli,si,ri=0, Gui,si,ri=1. Equation (27) is solved by dichotomy, and the upper confidence limit p^iu of failure probability pi is obtained. Based on the linearization process of (2), we have y^iu=lnln1p^iu and xi=lnti. According to the weighted least squares method, the upper limit curve of the failure probability is obtained by fitting points xi,y^iu. From (22), we have(28)Qu=i=1kwiy^iumxi+mlnη2.

Because the value of m is limited and the value has little effect on the reliability index value, the confidence interval of m is not assumed here. The point estimation value of m is m^. After obtaining the point estimation m^, the lower confidence limit ηL of the estimation value η is obtained. Putting m^ into (28), we have(29)Qu=i=1kwiy^ium^xi+m^lnη2.

The lower confidence limit ηL can be obtained when the value of Qu is the smallest:(30)ηL=expi=1kwixiyium^.

The point estimation value m^ and the lower confidence limit ηL are added to the Weibull reliability calculation formula, and then the lower confidence limit R of the reliability at the confidence level 1α can be obtained at any time t:(31)RLt=exptηLm^.

Similarly, the upper confidence limit for the probability of failure is(32)Fut=1RLt=1exptηLm^.

The upper confidence limit of the failure rate function is(33)λut=futRLt=m^ηLtηLm^1.

7. Simulation Verification

In order to verify the accuracy of the proposed method, experiments were used to generate simulation samples. Based on the generated simulation samples, the classical estimation method , the zero-failure estimation method , and the very few failure data evaluation method proposed in this study were used to evaluate the reliability of the samples and compare the evaluation result with the true value. The simulation experiment was divided into two cases, one was the simulation experiment without failure data and the other was the simulation experiment under the condition of very few failure data. The generation process of the simulation sample was described by Zhang et al. .

Generation of simulation samples with zero-failure data:

Set the true value of the shape parameter m and the scale parameter η in the Weibull distribution, and generate the failure time tii=1,2,,N through MATLAB, that is, generate a total of N experimental samples in the simulation experiment

Generate N uniformly distributed random numbers ui, ui0,1

ti is replaced by uiti to indicate no-failure truncation time ti

According to the value range of no-failure time ti, ti is divided into k group censored experiments. The ending data nj of each group is taken as the ending data, and the minimum value in the truncation data is taken as the truncation time, which constitutes the sample size tj,nj,j=1,2,,k

The generation of simulation samples under very few failure data conditions:

Set the true value of the shape parameter m and the scale parameter η in the Weibull distribution, and generate the failure time tii=1,2,,N through MATLAB, that is, generate a total of N experimental samples in the simulation experiment

Generate N uniformly distributed random numbers ui, ui0,1

Determine the censored level (CL), and on this basis, arbitrarily select rr=N×CL and the expiration time tfii=1,2,,r from all, denoted as tfii=1,2,,r, and the remaining expiration time tsii=1,2,,Nr is used

ti is replaced by uiti to indicate no-failure truncation time ti

According to the value of tfii=1,2,,r and tsj, the truncation time is grouped and the minimum truncation time is taken as the truncation time of the censored experiment of the group

It was assumed that 10 samples will be selected for 5 censored experiments. The shape parameters m=2 and scale parameters η=2500 in the Weibull distribution are known. The simulation samples were generated under the condition of zero-failure data, as shown in Table 1.

Simulation sample data.

Number Censoring time (h) Number of samples
1 435.357 1
2 938.124 1
3 1263.982 1
4 1517.114 2
5 1735.365 5

According to the time censored test data in Table 1, Figure 1 shows the point estimation value of the failure probability solved by the classical estimation method , the zero-failure estimation method , and the very few failure data reliability evaluation method proposed in this study. The absolute error with the true value is shown in Figure 2.

Comparison of the failure probability point estimation.

Comparison of the absolute error of the failure probability.

In the figure, Method 1 is the classical estimation method , Method 2 is the no-failure estimation method , and Method 3 is the minimum data failure reliability evaluation method proposed in this study. It can be seen from Figures 1 and 2 that the minimum data failure reliability assessment method can be used to estimate the failure probability of nonfailure data sample products, and the estimation result is closest to the true value. The classical estimation method  does not consider the truncation time versus failure in the time-censored test. The influence of probability makes the probability of failure “over-underestimated,” while the method of zero-failure estimation  has certain randomness in the process of determining the kernel of the prior distribution of failure probability, making the estimation probability of failure probability “large.”

As Jia et al.  and Jiang et al.  did not further analyze the confidence interval of the reliability parameter based on the estimation of the failure probability point, in order to facilitate the comparison, the confidence level was obtained by the optimal confidence limit method of equation (25). Under confidence level 0.9, the corresponding reliability lower confidence limit (LFM), the comparison results are shown in Figure 3.

Reliability confidence lower limit comparison.

Comparing Figure 1 with Figure 3, it can be seen that the lower confidence limit of reliability is calculated by using the optimal confidence limit method on the basis of the estimation results of failure probability points obtained by the classical estimation methods, namely, Method 1  and Method 2 . Compared with the estimation results of parameter points, there is a distinct “disconnection” problem, and the calculation results of the lower confidence limit are not credible. Method 3 proposed in this study effectively overcomes the problem of “disconnection” in the process of obtaining the lower confidence limit, and the calculation results are more reasonable and more reliable.

8. Case Study

Clutches in integrated transmissions are mostly hydraulically controlled wet multiplate clutches. This clutch allows the shifting process to be gentle, making it easy to shift power and automate the shifting process. The wet multiplate shifting clutch is composed of many annular double-sided steel sheets and friction plates, and the dual steel sheets and the friction plates are arranged between the two and can be used as an active sheet or a passive sheet. The friction sheet is generally used as an active sheet. The active piece is connected to the spline on the input shaft by an internal spline, and the passive piece is circumferentially fixed by the external spline and the internal spline of the clutch. The primary and passive sheets are in a cooling lubricant enclosed within the housing. When the clutch is in operation, the input shaft and the active piece connected thereto rotate, and the main and passive parts are pressed against each other by the hydraulic pressure to generate a friction torque, and the two synchronously rotate to realize power transmission. A section of a wet multiplate clutch is shown in Figure 4.

Sectional view of the integrated transmission clutch.

In the wet multiplate clutch, the normal state of the friction plate determines the key to whether the clutch can achieve its function. This study used the wet friction plate in the hydraulic control wet multiplate clutch as the research object and made a time censored test under normal working conditions. The test data are shown in Table 2.

Wet friction plate test data.

Subgroup number i 1 2 3 4
Censoring time τi/t 500 1000 1500 2000
Test samples ni 3 3 3 3
Remaining work samples si 12 9 5 3
Failure samples ri 0 0 1 0

From the number of failed samples in Table 2 and the actual project, the pi upper limit pupp0.5, so let pupp=0.5. The combination of (11) and (20) was used to calculate the each censoring time estimation of failure probability as p^ii=1,2,,4, where c=5. The calculation results are shown in Table 3.

Each censoring time ti point estimation value of failure probability point p^i.

Subgroup number i 1 2 3 4
Failure probability point estimation p^i 0.0256127 0.108449 0.258006 0.391425

Putting the data in Tables 2 and 3 into (23), m^=2.0985 and η^=2683.4 were calculated; the calculation result was substituted into (24) and the point estimation curve of reliability under the Weibull distribution was obtained. It is shown in Figure 5.

Reliability estimation curve.

This study sets the reliability R at the confidence level 1α=0.7, using (27), to obtain the point estimation of the upper bound p^iu of the failure probability pi. They are shown in Table 4.

Failure probability pi confidence upper limit point estimation.

Subgroup number i 1 2 3 4
Failure probability confidence upper bound point estimation p^iu 0.0259 0.1205 0.2890 0.4145

By y^iu=lnln1p^iu and xi=lnti, y^iu and xi, were obtained. They are shown in Table 5.

Transformed values y^iu.

Subgroup number i 1 2 3 4
y ^ i u −3.6404 −2.0526 −1.0756 −0.6249
x i 6.2146 6.9078 7.3132 7.6009

With (30) and (31), the lower confidence limit at any time t can be obtained when reliability R is at a confidence level of 0.7. It is shown in Figure 6.

Reliability confidence lower limit estimation curve.

With (33), the upper confidence limit of the failure rate function λut is obtained. It is shown in Figure 7.

Failure rate confidence upper limit point estimation curve.

9. Conclusion

In the case of the Weibull distribution, based on the time-censored test with very little failure data, the idea of the match distribution curve method combined with Bayesian theory and subtraction function method was used to study the reliability evaluation method, derived reliability parameter point estimation, and confidence interval estimation methods.

According to the different values of m in the Weibull distribution, the concave and convex properties of the distribution function and the function characteristics were used to determine the range of the failure probability. Combined with the reduction function of the failure probability, the parameters of the pretest distribution function were obtained. The value range further used the Bayesian method to obtain the point estimation failure probability pi based on the Bayesian poster distribution.

In view of the problem that the results of the reliability parameter estimation and the confidence interval estimation are “disjointed” in the current research, this study used the distribution curve method to unify the estimation method and avoid the “disjoint” phenomenon. The generation of the results and the comparison of the calculation results in the example further enhance the credibility of the results.

Applying the method proposed in this study to the data analysis of a time-censored test of the wet friction plate in a certain type of hydraulic control wet multiplate clutch, the applicability of the proposed method was further verified.

Data Availability

The data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This research was supported by the projects of the National Natural Science Foundation of China (grant nos. 71371183 and 71071158).

Kwon Y. Design of Bayesian zero-failure reliability demonstration test for products with Weibull lifetime distribution Journal of Applied Reliability 2014 14 4 220 224 Li H. Chen F. Yang Z. Bayesian Reliability Assessment Method for Single NC Machine Tool under Zero Failures. Applications and Techniques in Information Security 2015 Heidelberg, Germany Springer Martz H. F. Waller R. A. A bayesian zero-failure (BAZE) reliability demonstration testing procedure Journal of Quality Technology 1979 11 3 128 138 10.1080/00224065.1979.11980894 John Q. Estimating M. R. The probability of rare events: addressing zero failure date Risk Analysis 2011 31 7 1120 1132 Jia X. Jiang P. Guo B. The estimation to failure probability for zero-failure data in Weibull distribution Journal of Mechanical Strength 2015 37 2 288 294 Kurz D. Lewitschnig H. Pilz J. Advanced bauesian estimation of Weibull early life failure distributions Quality and Reliability Engineering International 2014 30 3 363 373 10.1002/qre.15772-s2.0-84896401028 Jiang P. Xing Y. Jia X. Guo B. Weibull failure probability estimation based on zero-failure data Mathematical Problems in Engineering 2015 2015 8 681232 10.1155/2015/6812322-s2.0-84930624678 Yin Y.-C. Hong H.-Z. Peng W. Li Y.-F. Mi J. An E-bayesian method for reliability analysis of exponentially distributed products with zero-failure date Maintenance and Reliability 2016 18 3 445 449 10.17531/ein.2016.3.172-s2.0-84976439831 Kan Y.-N. Yang Z.-J. Li G.-F. He J.-L. Wang Y.-K. Li H.-Z. Bayesian zero-failure reliability modeling and assessment method for multiple numerical control (NC) machine tools Journal of Central South University 2016 23 11 2858 2866 10.1007/s11771-016-3349-92-s2.0-85002642856 Chen J. D. Sun W. L. Li B. X. On the confidence limits in the case of no failure data Acta Mathematicae Applicatae Sinica 1995 18 1 90 100 Li H. Xie L. Li M. Ren J. Zhao B. Zhang S. Reliability assessment of high-quality and long-life products based on zero-failure data Quality and Reliability Engineering International 2019 35 1 470 482 10.1002/qre.23982-s2.0-85054180515 Jia X. The Estimation of Reliability for Satellite Platform with Multiply Type-І Censored Data 2017 Changsha, China National University of Defence Technology Xu T.-Q. Chen Y.-P. Two-sided M-bayesian credible limits of reliability parameters in the case of zero-failure data for exponential distribution Applied Mathematical Modelling 2014 38 9-10 2586 2600 10.1016/j.apm.2013.10.0542-s2.0-84898540016 Cheng K. Life Distribution Class and Reliability Mathematical Theory 1999 Beijing, China Science Press Jia X. Wang X. Guo B. Reliability assessment for very few failure data and zero-failure data Journal of Mechanical Engineering 2016 52 2 182 188 10.3901/jme.2016.02.1822-s2.0-84960367093 Lin Y.-J. Lio Y. L. Bayesian inference under progressive type-I interval censoring Journal of Applied Statistics 2012 39 8 1811 1824 10.1080/02664763.2012.6831702-s2.0-84862610380 Zhang J. K. Tang X. M. Bayes Methods 1989 Changsha, China National Defense Science and Technology University Press Mao S. S. Xia J. F. Guan W. Q. The reliability analysis of zero-failure data for bearing life test Chinese Journey of Applied Probability and Statistics 1993 9 3 326 331 Han M. Confidence limits of reliability parameters for Weibull distribution Journal of Mechanical Strength 2009 31 1 59 62 Zhang L. F. Xie M. Tang L. C. A study of two estimation approaches for parameter of Weibull distribution based on WPP Reliability Engineering and System Safety 2007 92 3 360 368 10.1016/j.ress.2006.04.0082-s2.0-33750980135