#### Abstract

This paper deals with the Bayes prediction of the future failures of a deteriorating repairable mechanical system subject to minimal repairs and periodic overhauls. To model the effect of overhauls on the reliability of the system a proportional age reduction model is assumed and the 2-parameter Engelhardt-Bain process (2-EBP) is used to model the failure process between two successive overhauls. 2-EBP has an advantage over Power Law Process (PLP) models. It is found that the failure intensity of deteriorating repairable systems attains a finite bound when repeated minimal repair actions are combined with some overhauls. If such a data is analyzed through models with unbounded increasing failure intensity, such as the PLP, then pessimistic estimates of the system reliability will arise and incorrect preventive maintenance policy may be defined. On the basis of the observed data and of a number of suitable prior densities reflecting varied degrees of belief on the failure/repair process and effectiveness of overhauls, the prediction of the future failure times and the number of failures in a future time interval is found. Finally, a numerical application is used to illustrate the advantages from overhauls and sensitivity analysis of the improvement parameter carried out.

#### 1. Introduction

A repairable system is a system that, after failing to perform one or more of its functions satisfactorily, can be restored to satisfactory performance.

Most repairable mechanical systems are subjected to degradation phenomena with operating time, so that the failures become increasingly frequent with time. Such systems often undergo a maintenance policy. Maintenance extends system's lifetime or at least the mean time to failure, and an effective maintenance policy can reduce the frequency of failures and the undesirable consequences of such failures. Maintenance can be categorized into two classes: corrective and preventive actions. Corrective maintenance, called repair, is all actions performed to restore the system to functioning condition when it fails. Preventive maintenance is all actions performed to prevent failures when the system is operating. Corrective and preventive maintenance actions are generally classified in terms of their effect on the operating conditions of the system. Pham and Wang [1] classified them as perfect maintenance, minimal maintenance, imperfect maintenance, and worse maintenance. At one extreme is the assumption of perfect maintenance, that a system is restored to good-as-new condition after maintenance. At the other extreme is the bad-as-old assumption that the failure rate of a system is not enhanced by maintenance. In real-world situations, maintenance generally enhances the condition of the equipment at a level between these two extremes, that is, imperfect maintenance.

We consider a system that deteriorates with age and receives two kinds of maintenance actions: minimal repair and overhaul. When a failure occurs, minimal repair is carried out. The minimal repair is a corrective maintenance action that brings the repaired equipment to the conditions it was just before the failure occurrence (bad-as-old). Hence, the reliability of the system decreases with operating time until it reaches unacceptable values. When it reaches unacceptable values or at prefixed epochs, preventive maintenance action (overhaul) is performed so as to improve the system condition and hence reduce the probability of failure occurrence in the following interval. However, overhaul cannot return the system to “good-as-new”, and thus it can be treated as imperfect repair. When the overhaul is effective, the reliability of the system improves significantly. An overhaul usually consists of a set of preventive maintenance actions such as oil change, cleaning, greasing, and replacing some worn components of the system.

Many imperfect repair models have already been proposed [1]. Malik [2] proposed a general approach to model the improvement effect of maintenance, where each maintenance reduces the age of the unit in the view of the rate of occurrences of failures. Each maintenance is assumed to reduce proportionally the operating time elapsed from the previous maintenance. Malik's proportional age reduction (PAR) model for imperfect maintenance is a generalization of good-as-new and bad-as-old. On the basis of this general model, Shin et al. [3] have proposed a PAR model which assumes that each major overhaul reduces proportionally the age of the equipment by a fraction of the epoch of the overhaul. Shin's paper proposes a method of data analysis to estimate the parameters of the failure process and the maintenance effect for a repairable unit. Two parametric families of intensity functions are considered, power law [4–11] and log linear [12]. To model different effects of rejuvenation of the preventive maintenance, two classes of models have been proposed by Jack [13].

In contrast to the classical approach used by them, Bayes approach has been used by several authors as it helps in incorporating prior information and/or technical knowledge on the failure mechanism and on the overhaul effectiveness into the inferential procedure. Pulcini [14] deals with the statistical analysis, from a Bayes viewpoint, of the failure data of repairable mechanical units subjected to minimal repairs and periodic overhauls. The effect of overhauls on the reliability of the system is modeled by PAR model and the power law process (PLP) is used to model the failure process between two successive overhauls. Pulcini provided Bayes point and interval estimation of model parameters (and functions thereof), as well as testing procedures on the effectiveness of the performed overhauls. Further, Pulcini [15] deals with the prediction, from a Bayes viewpoint, of future failures of repairable mechanical units subjected to minimal repairs and periodic overhauls. In [15], also PLP is used to model the failure process between two successive overhaul epochs, and the effect of overhauls on the reliability of the system is modeled by proportional age reduction model.

In PLP models, the increasing failure intensity tends to infinity as the system age increases. However, it is noted that the failure intensity of deteriorating repairable systems attains a finite bound when beginning from a given system age, repeated minimal repair actions are combined with some overhauls performed in order to oppose the growth of failure intensity with the operating time. The average behavior of the intensity function due to the consecutive steps with increasing intensity between two subsequent overhauls results in globally constant asymptotic intensity. If such data is analyzed through models with unbounded increasing failure intensity, such as the PLP, then pessimistic estimates of the system reliability will arise, and incorrect preventive maintenance policy may be defined.

Engelhardt and Bain [16] were the first to suggest NHPP with an increasing bounded intensity, called Engelhardt-Bain process (EBP). But the mathematical simplicity of their 1-parameter model makes it inadequate to analyze failure data in many cases. Most recently, Pulcini [17] proposed a 2-parameter NHPP called the bounded intensity process (BIP) to fit the failure data of deteriorating repairable systems showing bounded intensity function. Attardi and Pulcini [18] proposed a new bounded intensity process called 2-parameter Engelhardt-Bain process (2-EBP) whose failure intensity function is given by The failure intensity in (1) tends to its asymptote more slowly than the BIP intensity and is in a way a compromise between the PLP and BIP models. When , the intensity in (1) degenerates into the EBP intensity.

This paper deals with the Bayes prediction of the future failures of a deteriorating repairable mechanical system subjected to minimal repairs and periodic overhauls. To model the effect of overhauls on the reliability of the system, a proportional age reduction model is assumed, and the 2-parameter Engelhardt-Bain process (2-EBP) is used to model the failure process between two successive overhauls. On the basis of the observed data and of a number of suitable prior densities reflecting varied degrees of belief on the failure/repair process and effectiveness of overhauls, the prediction of the future failure times and the number of failures in a future time interval is found. Finally, a numerical application is used to illustrate the advantages of overhauls.

#### 2. Basic Assumptions

(i)Failure rate of the system is an increasing function of time that attains a finite bound as tends to infinity. (ii)System is subjected to two kinds of maintenance actions: minimal repair and overhaul.(iii)The times to perform maintenance actions are ignored. (iv)Minimal repair will restore the failure rate only to bad-as-old condition. But overhauls will improve the system to a condition between bad-as-old and good-as-new. (v)The failure density function is not changed by overhauls.(vi)The quality of an overhaul is dependent on improvement factor .(vii)The improvement parameter has a uniform value over all the overhaul actions. (viii)The th overhaul reduces the age of the system by a fraction of the epoch of the overhaul.(ix)The effect of overhauls on the reliability of the system is modeled by proportional age reduction model, and the 2-parameter Engelhardt-Bain process (2-EBP) is used to model the failure process between two successive overhauls, say .(x) are the overhaul epochs, which may coincide with failure times.#### 3. Model Formulation

Let denote the failure times of the repairable system observed till . If , that is, the process is observed till th failure, then it is failure-truncated sampling, and is a random variable. If , then it is time-truncated sampling, where is a prefixed quantity, and is a random variable.

The 2-EBP model is an NHPP whose failure intensity is of the form which is an increasing bounded intensity function with the operating time , equal to 0 at and approaching an asymptote of as increases. So, the parameter is the asymptotic value of the intensity function: , and the parameter is a measure of the initial increasing rate of : the smaller is, the faster the failure intensity increases until it approaches .

Thus, the initial failure intensity, that is, the intensity function till the first overhaul epoch , is The conditional intensity function at a generic time in the interval , given is, The expected number of failures between two successive overhaul epochs is The cumulative value of expected number of failures from 0 up to a generic time in the interval is

#### 4. Likelihood Function

The likelihood function based on observed data: , is

The conditional probability that the unit of age will fail before is The conditional reliability just before the th failure will be and the conditional pdf is .

If the data is time truncated, the joint pdf of the failure times is where is the truncation time.

As are the overhaul epochs, they provide disjoint fixed intervals . Since the overhaul epochs are prescheduled, they are not random. Let the time interval be the th period.

From (4) and (9a), for , the joint pdf of failure times in is Since the overhaul epochs are prescheduled, the failure times in disjoint epoch periods are independent, and therefore, the joint pdf of the failure times is simply the product form of (9b), that is, where .

#### 5. Bayesian Procedure

##### 5.1. Informative Prior

Suppose that the analyst is able to anticipate the following: (a)the asymptotic value of the intensity function, (b)the time at which the failure intensity is times the asymptotic value , that is, Now, we assume a Gamma density for and , both for the mathematical tractability and flexibility. Therefore, where the gamma parameters are related to the prior mean and standard deviation by The analyst can use the fact that a value closer to 0 (to 1) indicates less (more) effectiveness of the overhaul actions to formulate an informative prior on . We assume that the analyst is able to anticipate both a prior mean and a standard deviation on , say , then we formalize this information through a beta density where the values of the parameters related to the prior mean and standard deviation are

##### 5.2. Joint Prior Density

Assuming the prior independence of the parameters , and using (11) to find density of , the joint prior density is given by

##### 5.3. Joint Posterior Density on Model Parameters

By combining the likelihood (10) with the joint prior density (16), the joint posterior density on the model parameters is given by where

#### 6. Posterior Inference on the Expected Number of Failures in a Future Time Interval

Let be the expected number of failures in the future time interval , where If no overhaul is performed after , the conditional distribution of is Poisson with mean equal to the expected number of failures in the time interval , where The posterior mean of is given by the posterior expectation of , After the transformation of variables, we have the following: A conservative Bayes upper credibility limit for is obtained by finding the smallest value of which satisfies

#### 7. Posterior Inference on the Future Failure Times

Given the observed data, we are interested in predicting the th future failure time , under the assumption that no overhaul is performed after , where is gamma distributed with shape parameter m and unit scale parameter. Thus, the conditional pdf of , with , is given by where The predictive density function of is given by A Bayes prediction interval of a given probability content is defined as the interval such that

#### 8. Sensitivity Analysis

Sensitivity analysis is carried out with respect to the prior information on , namely the prior mean and the standard deviation , by evaluating : expected number of failures in the future time interval when overhaul is performed at and not performed at , over a reasonable range of values for and . If the slight change in values of and does not have much effect on then the proposed prior is robust.

#### 9. Numerical Application

Consider the following hypothetical data for illustrative purpose.

The failure times and overhaul epochs are given in Table 1, major overhauls marked with*. We have assumed it to be a time-truncated sample with failures observed for 1500 units. Four major overhauls are assumed to be performed at times different from failure times.

The proportional age reduction 2-EBP model is adequate for this hypothetical data set in contrast to pre-existing proportional age reduction power law process model as its log-likelihood value obtained using (10) is −95.4606 which is greater than the log-likelihood value −96.1537 of power law process model.

Suppose that analyst is able to anticipate a prior mean and .

In addition, from previous experiences, the analyst possesses a vague belief that the failure intensity, at the time units, is nearly half its asymptotic value: Then, he formalizes his prior knowledge on t_{r} through the exponential density having mean , so that and . As , therefore . Hence, .

Again, the analyst possesses a vague belief that the overhaul actions are quite effective, and then he chooses the beta density for the improvement parameter with prior mean and standard deviation .

Now, suppose that the above failure process is time truncated at units. On the basis of the observed data, we predict the future failure times under two hypotheses:(a)major overhaul is performed at , (b)no overhaul is performed at .

All further repair actions are assumed to be minimal repairs. The remaining observed values, that is, have been used for comparison with the predicted values.

In Table 2, we compare the occurred failure times with the Bayes prediction of the th future failure times under both the above hypotheses; in particular, the lower and upper limits of the 0.90 Bayes equal-tail credibility interval are shown. It clearly shows the advantages of performing overhauls as the future failure times are delayed considerably by overhaul.

In Table 3, we compare the observed values of with the lower and upper limits of the Bayes 0.90 equal-tail credibility interval and with the posterior median obtained using (30) and (31). This allows the analyst to assess the adequacy of the selected model and the correctness of prior information. Table 3 shows that all occurred failures fall well within the Bayes prediction interval and are quite close to the respective posterior median, thus indicating the adequacy of the 2-EBP to model the given data set.

Table 4 clearly shows that the, 95% upper credibility limit for the number of failures in future time interval (1200, 1500) when major overhaul is performed at is equal to 11 failures and the posterior mean = 5.327 failures, very close to the number of failures, namely 5, actually occurred in (1200, 1500).

On the contrary in the case of no overhaul 95% upper credibility limit for the number of failures in future time interval (1200, 1500) is equal to 14 failures which is more than that in major overhaul case, and the posterior mean is 7.612 which is again more than the value obtained in the major overhaul case. These results show the considerable advantages arising from performing a major overhaul at .

Tables 5 and 6 show that the change in values of and does not have much effect on , so the proposed prior is robust.

#### 10. Conclusion

In this paper, the prediction of future failures of a deteriorating repairable mechanical system subject to minimal repairs and periodic overhauls has been done using Bayesian approach. The effect of overhauls on the reliability of the system has been modeled using a proportional age reduction model, and the failure process between two successive overhauls has been modeled using 2-parameter Engelhardt-Bain process (2-EBP). The prediction of the future failure times and the number of failures in a future time interval has been done on the basis of the observed data and of a number of suitable prior densities reflecting varied degrees of belief on the failure/repair process and effectiveness of overhauls. The advantages of overhauls have been highlighted using a numerical application and sensitivity analysis of the improvement parameter carried out.

#### Nomenclature

: | Total no. of failures |

: | Total no. of overhauls |

: | th overhaul epoch, |

: | Parameters of 2-EBP |

: | Improvement factor |

: | Time at which the failure intensity is times the asymptotic value |

: | Initial failure intensity |

: | Conditional intensity function at a generic time in the interval |

: | Expected number of failures between two successive overhaul epochs |

: | Prior pdf on |

: | Prior pdf on |

: | Prior pdf on |

: | Prior gamma parameters of |

: | Prior gamma parameters of |

: | Prior beta parameters of |

: | Prior mean and variance of |

: | Prior mean and variance of |

: | Prior mean and variance of |

: | Joint prior pdf of |

: | Joint posterior pdf of |

: | Expected number of failures in the future time interval |

: | th future failure time |

NHPP: | Nonhomogeneous poisson process |

2-EBP: | 2-parameter Engelhardt-Bain process |

PLP: | Power law process. |

#### Acknowledgment

The authors are grateful to the referees for their valuable comments and suggestions that improved this paper by making it more informative.