Mathematical Problems in Engineering

Volume 2015, Article ID 681232, 8 pages

http://dx.doi.org/10.1155/2015/681232

## Weibull Failure Probability Estimation Based on Zero-Failure Data

^{1}College of Information System and Management, National University of Defense Technology, Changsha 410073, China^{2}College of Continuing Education, National University of Defense Technology, Changsha 410073, China

Received 27 May 2014; Revised 4 September 2014; Accepted 15 October 2014

Academic Editor: Wenbin Wang

Copyright © 2015 Ping Jiang 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.

#### Abstract

Reliability testing is often carried out with small sample sizes and short duration because of increasing costs and the restriction of development time. Therefore, for highly reliable products, zero-failure data are often collected in such tests, which could not be used to evaluate reliability by traditional methods. To cope with this problem, the match distribution curve method was proposed by some researchers. The key step needed to exercise this method is to estimate the failure probability, which has yet to be solved in the case of the Weibull distribution. This paper presents a method to estimate the intervals of failure probability for the Weibull distribution by using the concavity or convexity and property of the distribution function. Furthermore, to use the method in practice, this paper proposes using the approximate value of the shape parameter determined by either engineering experience or by hypothesis testing through a *p* value. The estimation of the failure probability is thus calculated using a Bayesian approach. A numerical example is presented to validate the effectiveness and robustness of the method.

#### 1. Introduction

Reliability testing is usually required in product development to evaluate product reliability. Product lifetimes are becoming longer than in previous decades because of the improvement of reliability. Therefore, with the restrictions of increasing costs and short development times, reliability testing is often carried out with small sample sizes and short duration, which will often lead to zero-failure data [1]. Accordingly, it is desirable to estimate product reliability using the zero-failure data.

The zero-failure scenario is described as follows.

Let denote the lifetime distribution of a product, where is the parameter of the distribution and is the parameter space. The reliability test is composed of truncated tests with corresponding censored time , which satisfies . The sample size for the th test is . As no failure is observed in the tests, the zero-failure data is expressed by , . Let denote the failure probability of the product at time . Then satisfies the following statements:(1), when ;(2).

Let denote the number of samples at time ; that is, there are samples in the reliability test at time .

Estimating product reliability based on zero-failure data is challenging. Welker and Lipow [1] first raised this problem and, since then, some researchers have made progress on the topic [2, 3]. For the binomial distribution, Bailey [4] proposes a model to predict failure probability from zero-failure data; however, his model requires a large sample size. Based on zero-failure data, Wang and Langanke [5] present an approach that compares the reliability index—mean time to failure (MTTF)—between a newly designed product and the old product, but their approach is based on the assumption that the shape parameters for the two products are the same, which limits the application in practice. Miller et al. [6] study the estimation of failure probability for software when no errors are observed in testing. Jiang et al. [7] construct the shrinkage preliminary test estimator (SPTE) to estimate the reliability of a product following a Weibull distribution, when a prior estimate is available. Chen et al. [8] introduce the optimal confidence limit method to obtain the optimal lower confidence limit of some parameters of reliability distributions, in the case of zero-failure data. However, their method only considers one sample in each test and could not be applied to the case of multiple samples in each test. As no failure is observed in the tests, the maximum likelihood estimation (MLE) approach could not be applied in such a scenario. To solve this problem, L. Wang and B. Wang [9] propose the modified maximum likelihood estimation (MMLE), introducing a new parameter, , to modify the results obtained by applying MLE on the zero-failure data. The key to MMLE is the value of parameter . However, some researchers have found that the parameter, , will often cause the overestimation of parameters. The overestimation usually results from the improper use of the zero-failure data, which has raised debates among researchers for a long time. One solution is to introduce failure information when dealing with zero-failure data [10]. The key to this method is to acquire the failure information (i.e., failure time), whereas it is often estimated by the zero-failure data. This will result in an inaccurate estimate when the test duration is far less than the real lifetime of a product [11]. Mao and Luo [12] present the match distribution curve (MDC) method to solve the evaluation on zero-failure data: first, estimating the failure probability at censoring time is carried out; then, the data pairs are used to construct a distribution curve to estimate the parameters of the distribution; finally, the reliability can be evaluated based on the distribution. The MDC method is widely used in various distributions with zero-failure data; however, it is not discussed in the case of a Weibull distribution because of the computational complexity of the distribution. Motivated by this problem, we focus our research on the failure probability estimation method in a Weibull distribution.

#### 2. Weibull Distribution

When evaluating reliability using test data, we often assume that the product life follows some kind of distribution. Then, the parameters of the distribution are estimated based on the data, and reliability is evaluated by using the determined distribution. The Weibull distribution is one of the most commonly used distributions in reliability evaluation because of its ability to take on various forms by adjusting its parameters [13]. The two-parameter Weibull distribution is defined aswhere is the shape parameter and is the scale parameter. The CDF of the Weibull distribution is defined as

The failure rate function is an increasing function, when , which describes the character of various products. The Weibull distribution will become another type of distribution when the value of shape parameter varies. For example, it becomes an exponential distribution when ; it becomes a Rayleigh distribution when ; and it approximates a normal distribution when . Ning [10] summarizes that a Weibull distribution can be used to describe traumatic failures, when , and it can be applied to describe degradation failures, when . Therefore, it can be adopted to describe the combined effects of traumatic failures and degradation failures, when , where the ratio of the two failures is determined by the value of . Due to the flexibility of Weibull distribution, it is widely used in reliability evaluation in practice, even in the cases of zero-failure data.

As the key step in the MDC method is failure probability estimation, which is yet to be solved in the case of a Weibull distribution, we present this estimation procedure in the following section.

#### 3. Failure Probability Estimation

Let , , and ; then (2) is transformed into

The linear regression method could be applied to calculate the parameters of the Weibull distribution. Here, the weighted least square estimation method is used to obtain the estimates of parameters and , which minimizes

As it is already proven that the weighted least square estimates, and , are unbiased, we propose the use of the method and briefly introduce it here.

Denote the weight , where is failure time. Let , , , and . The weighted least square estimates, and , are obtained by referring to

After and are acquired, given the test time , it is easy to calculate the failure probability at time for the Weibull distribution by

The above procedure is typical to obtain the failure probability at time for the Weibull distribution given that is the corresponding failure time. To calculate the failure probability at time for the Weibull distribution with zero-failure data, we have to consider an alternative solution. This is discussed in the next subsection for Weibull distribution.

##### 3.1. Failure Probability Estimation with Zero-Failure Data

In cases of zero-failure data, Ning [10] proposes the following equation to estimate at censoring time :

Equation (7) is designed to calculate the mean value of the upper limit and lower limit 0, which is too simple and arbitrary, and often results in inaccurate estimation.

Bayesian theory is the most popular method to estimate as it combines prior information with test data, making the estimation more accurate than if only prior information or test data was used alone. In Bayesian theory, distribution parameter is regarded as a random variable, and the prior distribution is determined by using historical data, experts’ judgments, or data from similar products to obtain the posterior distribution of bywhere is the PDF of the population , or it can be replaced by likelihood function , and is the prior distribution of . The point estimate of is , which is the expectation of the posterior distribution under the square loss.

In the case of zero-failure data, to estimate , the likelihood function is usually defined as . Therefore, in terms of Bayes’ theorem, the key to estimating is to choose the prior distribution of . To achieve this, two important answers should be sought here: the right choice of the prior distribution and the determination of the interval of in the prior distribution.

*(1) The Right Choice of the Prior Distribution of *. For the choice of the prior distribution of , one may simply assume that follows the uniform distribution. However, since the data exhibits zero failure, the value of failure probability is more likely to be small. Therefore, Han and Li [14] present the idea of utilizing the decreasing function to construct the prior distribution. We also adopt this idea and use it to construct the prior distribution in the next section.

*(2) The Interval of *. To determine the interval of for a normal distribution, Zhang [15] proves that the CDF of is a concave function in , when . Then, based on the properties of the concave function, Zhang concludes that the following inequalities stand:Then the interval of is , where is the upper limit preset by expert or engineering experience.

Inspired by this theory, and given the distribution type, researchers tried to use the properties of the distribution to determine the interval of if there is no other prior information: in the case of an exponential distribution, Ning [10] obtains the interval of as , and in the case of an extreme distribution, Li [11] obtains the interval of as . However, in the case of a Weibull distribution, there is no reference to determine the interval of , which is why we choose to study this problem according to the convex and concave properties of the Weibull distribution.

##### 3.2. Weibull Failure Probability Estimation

Based on the accumulated test data from Weibull-distributed products in worldwide practice, Han [16] states that the shape parameter, , for Weibull-distributed products is usually within the interval of . So, the interval for is set to be .

To obtain the convex and concave properties of the Weibull distribution, by taking the second derivative of (2), we have

From (10), it is easy to find that , when , and when , let the inflexion point ; then , when ; , when . So, according to the convex and concave criteria, we have the following: is convex, when , and is concave, when and . This conclusion provides us with the criteria for determining the interval of .

When , is convex, so we have

As we have assumed that , the following inequality stands:

Therefore, the interval of is , which could be calculated by an iterative method beginning from ().

The calculation begins with , which can be estimated by the Bayesian method. Let the interval of be , where is the upper limit set by experts. In practice, is often used. However, this value of is conservative here, as no failure occurs; we can assume that , so could be set at less than 0.5.

When , is not strictly convex or concave, which makes it difficult to analyze. Therefore, we present our solution as follows.

We can acquire the median time by referring to (1), when . The inflexion point (time) is obtained in a similar manner, which is denoted by .

As the data reveals no failure, we assume that . Therefore, to find the interval, within which exhibits convexity or concavity, we will determine the interval that satisfies .

Let and , where . Then it is easy to prove that both and are strictly increasing in and , when , as shown in Figure 1. Therefore, when , we have , and is concave.