A joint optimal policy of inspection and age based replacement based on a threestage failure process
Ruifeng Yang^{1} , Jianshe Kang^{2} , Huijuan Wang^{3}
^{1, 2}Mechanical Engineering College, Shijiazhuang, 050003, China
^{3}Zhejiang University of Finance and Economics, Hangzhou, 310000, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 17, Issue 6, 2015, p. 30033022.
Received 17 December 2014; received in revised form 2 September 2015; accepted 10 September 2015; published 30 September 2015
JVE Conferences
Preventive maintenance (PM) and conditionbased maintenance (CBM) are two dominant maintenance policies in industrial applications. Inspection activities are the foundation of PM and CBM policies as to provide the operating information of system through processing the collected vibration data. Age based replacement is one of the most used preventive maintenance policy aiming at avoiding unplanned downtime and higher failure loss. This paper proposes a joint optimal policy of inspection and age based replacement based on a threestage failure process for a single component system. The threestage failure process, which is closer to reality, divides the failure process of system into three stages: namely normal, minor defective and severe defective. When the severe defective stage is identified, maintenance action is carried out immediately. The system is replaced once it reaches certain age. However, two potential actions are considered and analyzed in this paper when the minor defective stage is identified: halving the subsequent inspection interval or replacing the item immediately. As inspection may not be perfect because of the complexity of plant items, both perfect and imperfect inspection cases are considered. Finally, a case study is presented to demonstrate the efficiency of the proposed models.
Keywords: preventive maintenance, conditionbased maintenance, inspection, delay time, threestage failure process, age based replacement.
1. Introduction
Preventive maintenance (PM) and conditionbased maintenance (CBM) are two popular maintenance policies in practice. Preventive maintenance, by which maintenance is implemented with a planned interval aiming at preventing potential failures from occurring, is still the most widely used policy in industry due to its easy implementation [1]. Nowadays CBM attracts more interest because of the efficiency of deriving the actual condition information [2].
Vibration data can express the condition of system exactly when operating. Through collecting and processing the vibration data, maintenance operators could make maintenance decisions on the basis of the actual condition obtained from inspection. Therefore, inspection is the foundation of implementing the PM and CBM policies as it can provide the information on the status of system checked to facilitate the determination and execution of repair and replacement decisions [3, 4]. Nowadays inspection based on vibration information has been widely applied in industry. Of course, inspection can be online or offline. As to the offline inspection, maintenance operators often inspect the system using handheld inspection tools, such as SPM (Shock pulse method) instrument.
SPM has been widely used as a quantitative method in determining bearing health condition. Unlike vibration analysis that monitors a broad vibration band and then tries to isolate unique frequencies, SPM has developed a means to only “look” at the high frequency signals of rotating bearings [5, 6]. Through monitoring and analyzing high frequency shock waves generated by a bearing while rotating, SPM can provide a direct shock value indicating the bearing condition. The major benefit of SPM is providing a direct indication of bearing condition on a GreenYellowRed scale. Green means a good bearing, Yellow is a bearing with early damage and Red is more severe damage [7]. This is very important when monitoring installed bearings on machines for which there are no trends or comparable readings. Through years of application, SPM has been proven successfully as a diagnostic tool for bearings.
As maintenance operators inspect the bearings using SPM with a planned interval, the determination of inspection intervals is one of the key decisions. There are numerous PM models to determine the optimal inspection intervals [8, 9]. The delaytime concept, which assumes the state of system before failure to be either normal or defective, has obvious advantage for optimizing the inspection intervals because it can capture the relationship between inspection interval and number of failures [10]. Most importantly, numerous successful case studies have been published with actual applications in industry using delaytime concept [1113].
However, the state of system may by described more than two before failure as three color scheme is mostly used to divide the state before failure into green (normal), yellow (need attention), and red (need immediate attention) in industrial applications [14]. For example, SPM divides the lifetime of bearings into three stages, namely normal, minor defective and severe defective, expressed by green, yellow and red respectively. Therefore, Wang firstly extended the twostage delay time to a threestage failure process, where the traditional failure delay time stage is further divided into two stages: minor defective stage and severe defective stage [15]. Finally a case study is presented to show that delaying maintenance is better than maintenance immediately and validate the efficiency of the proposed threestage failure process. A preventive maintenance model with a two level inspection policy based on a threestage failure process is proposed by Wang et al. [16]. Minor inspection and major inspection consist of the twolevel inspection policy. Moreover, a planned PM is considered as maintenance window to prevent the occurrence of failure. Yang et al. [17, 18] developed imperfect inspection maintenance models based on a threestage failure process. Monte Carlo simulation is used to present the efficiency of the threestage failure process. An inspection optimization model based on a threestage failure process is presented to optimize the shortening proportion of inspection interval when the minor defective stage is identified [19]. But age based replacement policy is not considered in [1519].
Age based replacement policy is one of the most commonly used preventive maintenance and replacement policies. According to this policy, the system is replaced upon failure or at fixed age, whichever occurs first [20]. Preventive replacement is carried out to avoid unplanned downtime and higher failure loss, thus age based replacement policy is widely used to determine optimum preventive replacement time. The joint policy of regular inspections and age based preventive replacement has been widely used in industrial applications. Most researches have considered the age based preventive replacement policy [2123], but joint optimization for inspection and replacement intervals based on the threestage failure process has not been developed.
In this paper, a joint optimal policy of inspection and age based replacement based on a threestage failure process is introduced. Both cases of perfect and imperfect inspections are considered. The system is replaced when it fails or reaches the certain age. Once the severe defective stage is identified, repair will be carried out. However, when the minor defective stage is identified, two options are assumed here for comparison: halving the subsequent inspection interval or repair it immediately. The longrun availability is used to jointly optimize the inspection and replacement intervals.
The remaining part of the paper is organized as follows. Section 2 introduces the modeling assumptions and notation. Section 3 provides the proposed availability models with imperfect inspection and perfect inspection case is modeled in Section 4. Section 5 presents a case study and Section 6 concludes the paper.
2. Modelling assumptions and notation
The following modeling assumptions and notation are presented for model building.
1) A single component system is considered, which only subject to a single failure mode.
2) The failure process of system is divided into three stages: normal stage ${X}_{1}$, minor defective stage ${X}_{2}$ and severe defective stage ${X}_{3}$. These three stages are assumed to be independent.
3) The system is subjected to inspection interval $t$ and age based replacement interval $T$. Age based replacement interval is a multiple of inspection interval as $T=Nt$ for simplicity.
4) Both perfect and imperfect inspection cases are considered. Perfect inspection can always reveal the defective stage no matter it is minor and severe. However, imperfect inspection may miss the minor defective stage with a probability of ($1r$), but always can identify the severe defective stage.
5) If the system is found in the normal stage, do nothing.
6) On condition that the system is identified in the minor defective stage by an inspection, two options are assumed here, the first one is to shorten the subsequent inspection interval to be half of the current interval, and the other one is to repair immediately.
7) System is found to be in the severe defective stage, it is always repaired immediately.
8) Failure can be observed immediately and replacement is always carried out at once.
9) Replacement renews the system once the age of system reaches certain age $T$, which is an age limit of system.
10) Repair or replacement is regarded as renewing the system, though it may be the only option for a single component system.
3. Availability models with imperfect inspection
The availability models are formulated in this section to jointly optimize the inspection and age replacement intervals. The imperfect inspection case mainly means that an inspection may miss the minor defective stage with a probability of ($1r$) but can always identify the severe defective stage. Fig. 1 presents the joint policy of imperfect inspection and age based replacement.
Fig. 1. The system fails in $\left(\left(i1\right)t,it\right)$ before identifying any defective stages
Maximizing the availability is a main object of maintenance optimization. The decision objective is to derive the maximum longrun availability. Therefore, the long run expected renewal cycle downtime and length need to be derived firstly so that [24]:
According to the different decisions when the minor defective stage of system is identified depending on the assumption Eq. (6), two models are proposed.
3.1. Model 1
The system is renewed at failure, at the time of age based replacement and that a severe defect was first found prior to failure. However, if minor defective stage is identified, subsequent inspection intervals are halved.
3.1.1. The expected renewal cycle downtime
There exist three different kinds of renewal scenarios, namely failure renewal, inspection renewal caused by the identification of severe defective stage and age based replacement renewal. The downtime caused by these renewals in a renewal cycle consist the expected renewal cycle downtime, which will be introduced respectively.
3.1.1.1. The failure renewal scenarios
Failure renewal means the system fails in an inspection interval before identifying the severe defective stage or reaching the certain age $T$. However, according to the case whether the minor defective stage is identified by an inspection, there are two different failure renewals.
1) The system fails in ${T}_{f}$, ${T}_{f}\in \left(\right(i1)t,it)$ before any defective stages are identified, as shown in Fig. 2. The minor defective stage starts in $\left(\left(j1\right)t,jt\right)$ and ends in $\left(\left(i1\right)t,it\right)$. All the inspections within $\left(jt,\left(i1\right)t\right)$ miss the minor defective stage because of the imperfect inspection. The severe defective stage and failure happen within the same interval $\left(\left(i1\right)t,it\right)$. The probability of such an event is given as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=P\left(\right(j1)t{X}_{1}jt,(i1)t{X}_{1}{X}_{2}it{X}_{1},0{X}_{3}it{X}_{1}{X}_{2})$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{\int}_{\lambda \left(\right(i1)tx)}^{itx}(1r{)}^{ij}{f}_{{X}_{1}}(x\left){f}_{{X}_{2}}\right(y\left){F}_{{X}_{3}}\right(itxy)dydx,$
where $i=$1, 2,…, $N$. Here we define $\lambda \left(x\right)=x$ if $x\ge 0$, otherwise $\lambda \left(x\right)=0$.
Fig. 2. The system fails in $\left(\left(i1\right)t,it\right)$ before identifying any defective stages
2) The system fails in ${T}_{f}$, ${T}_{f}\in (it+(k1)t/2,it+kt/2)$ before identifying the severe defective stage by an inspection. The system degradates into minor defective stage within $\left(\left(j1\right)t,jt\right)$. However, the minor defective stage is identified at $it$ and the subsequent inspection interval is halved, as depicted in Fig. 3. Similarly, the severe defective stage and failure happen within the same interval $\left(it+\left(k1\right)t/2,it+kt/2\right)$. The probability of such a renewal is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=P\left(\begin{array}{c}\left(j1\right)t{X}_{1}jt,\\ it+\left(k1\right)\frac{t}{2}{X}_{1}{X}_{2}it+\frac{kt}{2}{X}_{1},\\ 0{X}_{3}it+\frac{kt}{2}{X}_{1}{X}_{2}\end{array}\right)$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{\left(j1\right)t}^{jt}{f}_{{X}_{1}}\left(x\right){\int}_{it+\left(k1\right)\frac{t}{2}x}^{it+\frac{kt}{2}x}{f}_{{X}_{2}}\left(y\right)(1r{)}^{ij}r{F}_{{X}_{3}}\left(it+\frac{kt}{2}xy\right)dydx.$
Fig. 3. The system fails in $\left(it+\left(k1\right)t/2,it+kt/2\right)$ after identifying minor defective stage at $it$
Therefore, the expected renewal cycle downtime due to failure renewal can be derived from Eqs. (2) and (3):
$+\sum _{i=1}^{N1}\sum _{k=1}^{2(Ni)}\left({D}_{f}+\left(i+k1\right){D}_{s}\right)P\left(it+\left(k1\right)\frac{t}{2}{T}_{f}it+\frac{kt}{2}\right).$
3.1.1.2. The inspection renewal scenarios
Inspection renewal occurs once the defective stages are identified. However, when the minor defective stage is identified, the subsequent inspection interval is halved. So the inspection renewal here only means the renewal caused by the identification of severe defective stage. Similarly, there are also two different renewals depending on whether the minor defective stage is identified.
1) Fig. 4 shows the scenario that the system is renewed when the severe defective stage is found at $it$ before any minor defective stage is identified. The inspections miss the minor defective stage $\left(ij\right)$ times. So the probability of such an inspection renewal is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right){\int}_{\lambda \left(\right(i1)tx)}^{itx}{f}_{{X}_{2}}\left(y\right)(1r{)}^{ij}(1{F}_{{X}_{3}}(itxy))dydx.$
Fig. 4. The system renews at $it$ before any minor defective stage is identified
2) Fig. 5 presents the scenario that inspection renewal occurs at $it+kt/2$ when the minor and severe defective stages are identified at $it$ and $it+kt/2$ respectively. The inspections within $\left(jt,\left(i1\right)t\right)$ miss the minor defective stage but identify it at $it$ and then the subsequent inspection interval is halved:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=P\left(\begin{array}{c}\left(j1\right)t{X}_{1}jt,\\ it+\frac{\left(k1\right)t}{2}{X}_{1}{X}_{2}it+\frac{kt}{2}{X}_{1},\\ {X}_{3}it+\frac{kt}{2}{X}_{1}{X}_{2}\end{array}\right)$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{\left(j1\right)t}^{jt}{f}_{{X}_{1}}\left(x\right){\int}_{it+\frac{\left(k1\right)t}{2}x}^{it+\frac{kt}{2}x}{f}_{{X}_{2}}\left(y\right)(1r{)}^{ij}r\left(1{F}_{{X}_{3}}\left(it+\frac{kt}{2}xy\right)\right)dydx.$
Fig. 5. The system renews at $it+kt/2$ after the minor defective stage is identified at $it$
Then the expected renewal cycle downtime of an inspection renewal can be derived using Eqs. (5) and (6) as:
3.1.1.3. The age based replacement renewal scenarios
The system is renewed once it reaches certain age $T$ according to the assumption Eq. (9). There are also two scenarios.
1) The age based replacement renewal is happened before any defective stages are identified. However, there may exist three possible scenarios, see Fig. 6. Firstly, the normal stage ${X}_{1}$ exceeds the age $T$. Secondly, the normal stage ends within $\left(\left(j1\right)t,jt\right)$. The minor defective stage is longer than $T{X}_{1}$ but missed by inspections. Finally, the minor defective stage starts within $\left(\left(j1\right)t,jt\right)$, ends in $\left(\left(N1\right)t,T\right)$ but all missed. The severe defective stage is longer than $T{X}_{1}{X}_{2}$. The probability of such a renewal is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+P\left(\left(j1\right)t{X}_{1}jt,\left(N1\right)t{X}_{1}{X}_{2}Nt{X}_{1},{X}_{3}Nt{X}_{1}{X}_{2}\right)$
$={\int}_{Nt}^{\mathrm{\infty}}{f}_{{X}_{1}}\left(x\right)dx+\sum _{j=1}^{N}{\int}_{\left(j1\right)t}^{jt}(1r{)}^{Nj}{f}_{{X}_{1}}(x\left)\right(1{F}_{{X}_{2}}(Ntx))dx$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{j=1}^{N}{\int}_{\left(j1\right)t}^{jt}(1r{)}^{Nj}{f}_{{X}_{1}}(x){\int}_{\left(N1\right)tx}^{Ntx}{f}_{{X}_{2}}\left(y\right)\left(1{F}_{{X}_{3}}\left(Ntxy\right)\right)dydx.$
Fig. 6. The system renews at age $T$ before any defective stages are identified
2) The system renews at age $T$ after the minor defective stage is found at $it$. There also exist two different scenarios, see Fig. 7. Firstly, The normal stage ends in $\left(\left(j1\right)t,jt\right)$ but the minor defective stage is longer than $T{X}_{1}\text{.}$ The second is the minor defective stage ends in $\left(\left(N1\right)t+t/2,T\right)$ and the severe defective stage exceeds the length of $T{X}_{1}{X}_{2}$. All the inspections between $jt$ and $\left(i1\right)t$ miss the minor defective stage in both scenarios:
$+P\left(\right(j1)t{X}_{1}jt,(N1)t+\frac{t}{2}{X}_{1}{X}_{2}Nt{X}_{1},{X}_{3}Nt{X}_{1}{X}_{2}$
$=\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right)(1r{)}^{ij}r(1{F}_{{X}_{2}}(Ntx))dx$
$+\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right)(1r{)}^{ij}r{\int}_{(N1)t+\frac{t}{2}x}^{Ntx}{f}_{{X}_{2}}\left(y\right)(1{F}_{{X}_{3}}(Ntxy\left)\right)dydx.$
Fig. 7. The system renews at age $T$ after identifying the minor defective stage at $it$
The expected renewal cycle downtime of age based replacement renewal can be derived from Eqs. (8) and (9) as:
3.1.2. The expected renewal cycle length
These three kinds of expected renewal cycle length can be derived according to the three different renewal scenarios.
3.1.2.1. The failure renewal scenarios
In order to derive the expected cycle length of failure renewal, the pdf of failure needs to be formulated. Firstly, if the system fails in $\left(\left(i1\right)t,it\right)\text{,}$ the pdf of failure at ${T}_{f}\text{,}$${T}_{f}\in \left(\right(i1)t+z,(i1)t+z+dz)$, $z\in (0,t)$ can be derived from Eq. (2) as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{\int}_{\lambda \left(\right(i1)tx)}^{(i1)t+zx}(1r{)}^{ij}{f}_{{X}_{1}}(x\left){f}_{{X}_{2}}\right(y\left){f}_{{X}_{3}}\right((i1)t+zxy)dydxdz.$
Similarly, the pdf of failure at ${T}_{f}$, ${T}_{f}\in (it+(k1)t/2+z,it+(k1)t/2+z+dz)$, $z\in (0,t/2)$ can be derived from Eq. (3) as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right){\int}_{it+(k1)\frac{t}{2}x}^{it+(k1)\frac{t}{2}+zx}{f}_{{X}_{2}}\left(y\right)(1r{)}^{ij}r{f}_{{X}_{3}}\left(it+(k1)\frac{t}{2}\right.$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+zxy)dydxdz.$
Accordingly, the expected renewal cycle length caused by failure renewal is expressed as:
$+\sum _{i=1}^{N1}\sum _{k=1}^{2\left(Ni\right)}{\int}_{0}^{\frac{t}{2}}\left(it+\left(k1\right)\frac{t}{2}\mathrm{}+z\right)P\left(it+\left(k1\right)\frac{t}{2}\mathrm{}+z{T}_{f}\right.it$
$\left.\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\left(k1\right)\frac{t}{2}+z+dz\right).$
3.1.2.2. The inspection renewal scenarios
The expected renewal cycle length caused by an inspection renewal is given as:
3.1.2.3. The age replacement renewal scenarios
The expected renewal cycle length of age based replacement renewal is:
3.1.3. The longrun availability
Based on the different expected renewal cycle downtime and length, the longrun availability is given using the renewal theorem [25], as shown in Eq. (16):
3.2. Model 2
The system is replaced at the time of defect identification no matter it is minor or major prior to failure. The system is also renewed at failure and at the time of age based replacement.
In model 2, there also exist three kinds of renewals, namely failure renewal, age based replacement renewal and inspection renewal. However, inspection renewal may be caused by the identification of minor or severe defective stage, which is different from model 1. The expected renewal cycle downtime and length are derived as follows.
3.2.1. The expected renewal cycle downtime
3.2.1.1. The failure renewal scenarios
As the system is renewed at the time of the identification of two defective stages, failure only happens in the same interval before any defective stages are found. Therefore, the probability of failure renewal is same to the Eq. (2). The expected downtime of failure renewal can be derived as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+(i1){D}_{s}){\int}_{\left(j1\right)t}^{jt}{\int}_{\lambda \left(\left(i1\right)tx\right)}^{itx}(1r{)}^{ij}{f}_{{X}_{1}}(x\left){f}_{{X}_{2}}\right(y\left){F}_{{X}_{3}}\right(itxy)dydx.$
3.2.1.2. The inspection renewal scenarios
The inspection renewal may be carried out once the defective stages are identified no matter it is minor or severe. Therefore, two different renewal scenarios are considered.
1) The system is renewed at it when the severe defective stage is found firstly before the minor defective stage is identified. The probability of this scenario is same to the Eq. (5). Thus the expected renewal cycle downtime of this scenario is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{i=1}^{N1}\sum _{j=1}^{i}({D}_{ps}+i{D}_{s}){\int}_{\left(j1\right)t}^{jt}{f}_{{X}_{1}}\left(x\right){\int}_{\lambda \left(\left(i1\right)tx\right)}^{itx}{f}_{{X}_{2}}\left(y\right)(1r{)}^{ij}(1{F}_{{X}_{3}}(it$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}xy\left)\right)dydx.$
2) The system is renewed at it when the minor defective stage is found firstly prior to failure, see Fig. 8. All the inspections between $\left(\left(j1\right)t,it\right)$ miss the minor defective stage. The probability of such a scenario is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{j=1}^{i}{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right)(1r{)}^{ij}r(1{F}_{{X}_{2}}(itx))dx.$
Fig. 8. The system renews at it when the minor defective stage is identified
Accordingly, the expected renewal cycle downtime of this renewal is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{i=1}^{N1}\sum _{j=1}^{i}({D}_{pm}+i{D}_{s}){\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right)(1r{)}^{ij}r(1{F}_{{X}_{2}}(itx))dx.$
3.2.1.3. The age replacement renewal scenarios
The probability of this scenario is the same to the Eq. (8). Thus the expected renewal cycle downtime of this event is given as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{j=1}^{N}{\int}_{\left(j1\right)t}^{jt}(1r{)}^{Nj}{f}_{{X}_{1}}(x\left)\right(1{F}_{{X}_{2}}(Ntx))dx$
$\left.\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{j=1}^{N}{\int}_{(j1)t}^{jt}(1r{)}^{Nj}{f}_{{X}_{1}}(x){\int}_{(N1)tx}^{Ntx}{f}_{{X}_{2}}\left(y\right)(1{F}_{{X}_{3}}(Ntxy\left)\right)dydx\right).$
3.2.2. The expected renewal cycle length
Similarly, three different renewal scenarios should be considered in order to derive the expected renewal cycle length.
3.2.2.1. The failure renewal scenarios
The system fails in an inspection interval before any defective stage is identified. The pdf of such a renewal is similar to Eq. (11). So the expected cycle length of this failure renewal is given by:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{i=1}^{N}\sum _{j=1}^{i}{\int}_{0}^{t}{\int}_{(j1)t}^{jt}{\int}_{\lambda \left(\right(i1)tx)}^{(i1)t+zx}\left(\right(i1)t+z)(1r{)}^{ij}{f}_{{X}_{1}}(x\left){f}_{{X}_{2}}\right(y\left){f}_{{X}_{3}}\right((i1)t$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+zxy)dydxdz$.
3.2.2.2. The inspection renewal scenarios
The expected cycle length of an inspection renewal caused by the identification of severe defective stage is:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{i=1}^{N1}\sum _{j=1}^{i}it{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right){\int}_{\lambda \left(\right(i1)tx)}^{itx}{f}_{{X}_{2}}\left(y\right)(1r{)}^{ij}(1{F}_{{X}_{3}}(itxy))dydx.$
Moreover, the expected cycle length of an inspection renewal caused by the identification of minor defective stage is given as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}=\sum _{i=1}^{N1}\sum _{j=1}^{i}it{\int}_{(j1)t}^{jt}{f}_{{X}_{1}}\left(x\right)(1r{)}^{ij}r(1{F}_{{X}_{2}}(itx))dx.$
3.2.2.3. The age replacement renewal scenarios
The expected renewal cycle length of age based replacement renewal is given by:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{j=1}^{N}{\int}_{(j1)t}^{jt}(1r{)}^{Nj}{f}_{{X}_{1}}(x\left)\right(1{F}_{{X}_{2}}(Ntx))dx$
$\left.\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{j=1}^{N}{\int}_{(j1)t}^{jt}(1r{)}^{Nj}{f}_{{X}_{1}}(x){\int}_{(N1)tx}^{Ntx}{f}_{{X}_{2}}\left(y\right)(1{F}_{{X}_{3}}(Ntxy\left)\right)dydx\right).$
3.2.3. The longrun availability
In the similar way as Eq. (16), the longrun availability of model 2 is expressed as:
4. Availability models with perfect inspection
The case of imperfect inspection has been presented in Section 3. Perfect inspection means inspection can always reveal the state of system. The defective stage will be found immediately once it occurs no matter it is minor or severe. Similarly, three different renewal scenarios are considered, namely failure renewal, age based replacement renewal and inspection renewal with the identification of minor or severe defective stage. The availability model with perfect inspection follows the same principle as the imperfect inspection case. Therefore, the derivation of renewal probabilities with perfect inspection will be left to Appendix.
5. Case study
A case study is presented in this section to demonstrate the efficiency of the proposed models and derive the optimal inspection and replacement intervals.
5.1. Wind turbine gearbox experiment
The experimental system used in this paper mainly concentrates on the bearings condition inspection of wind turbines. The gearbox of industrial wind turbine is inspected aperiodically by SPM instrument, as shown in Fig. 9. The system includes two gearboxes, one is speeddown gearbox, other is speedup gearbox. Four motors are used, a 1000 kw for driving the gearboxes and three motors for loading with power 355 kw, 500 kw and 655 kw respectively. A speed and torque sensor is to measure the rotating speed. The speedup gearbox is the test gearbox. In the test gearbox, the ring gear is stationary, a sun gear rotates around a ﬁxed center, and planet gears not only rotate around their own centers but also revolve around the center of the sun gear. The planet gears mesh simultaneously with both the sun gear and the ring gear. The inner structure of the test gearbox and measuring points of SPM sensor are depicted as Fig. 10.
Fig. 9. The gearbox test rig of FL600 wind turbine
The maintenace data we collected include threefolds: the time that the bearings starts to operate; the time that the bearings are replaced because of the identified defect, age replacement or failure; the time that last inspection before renewal or replacement.
After a period of experiment, the maintenance data of bearings are collected, as shown in Table 1.
Table 1. The renewal data of rolling bearings
Start time

Happen time

Minor defect

Severe defect

Age replacement

Failure

Length

Inspection time

0

219

√

219

194


325

370

√

45

21


398

413

√

15

0


459

497

√

38

17


518

717

√

199

149


727

785

√

58

29


801

878

√

77

54


883

924

√

41

26


939

952

√

13

0


965

1022

√

57

28


1037

1109

√

72

42


1112

1407

Normal when the inspection is ended

283

Fig. 10. The test gearbox and measuring points location
5.2. Maximum likelihood method
The maximum likelihood method are used to estimate the parameters as this method only needs to konw the probability density function of each renewal event. The parameters are estimated through maximizing the multiplication of each density function.
The inspection and age replacement intervals are assumed as $t$ and $T$ respectively. Maintenance will be implemented once the minor defective stage is identified. Otherwise, replacement will be carried out.
After the inspection for a period of time, we can derive some information as follows:
1) Total ${N}_{m}$ times of minor defect renewal are occurred at the time (${t}_{1}$, ${t}_{2}$,…, ${t}_{j}$,…, ${t}_{Nm}$) respectively;
2) Total ${N}_{s}$ times of severe defect renewal are occurred at the time (${t}_{1}$, ${t}_{2}$,…, ${t}_{n}$,…, ${t}_{Ns}$) respectively;
3) Total ${N}_{v}$ times of age replacement renewal are occurred at the time (${T}_{1}$, ${T}_{2}$,…, ${T}_{u}$,…, ${T}_{Nv}$) respectively;
4) Total ${N}_{f}$ times of failure renewal are occurred at the time (${x}_{1}\text{,}$${x}_{2}$,…, ${x}_{q}\text{,}$…, ${x}_{Nf}$) respectively;
5) The bearings are normal at the ending time $t$.
The likehood function can be derive by multipling the probabilty of each renewal event as:
where ${P}_{m}\left({t}_{j}\right)$: The probability of $j$th minor defect renewal at ${t}_{j}$; ${P}_{s}\left({t}_{n}\right)$: The probability of $n$th severe defect renewal at ${t}_{n}$; ${p}_{f}\left({x}_{q}\right)$: The probability of $q$th failure renewal at ${x}_{q}$; ${P}_{AR}\left({T}_{u}\right)$: The probability of $u$th age replacement renewal at ${T}_{u}$; ${P}_{n}\left(t\right)$: The probability of no renewal event at the ending time $t$.
The Eq. (27) can be changed as:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\mathrm{l}\mathrm{o}\mathrm{g}\left[{P}_{n}\left(t\right)\right].$
From the Eq. (28), the log likehood function can be derived according to the probability of each renewal event. The estimated values of parameters are achieved by maximizing the log likehood function. Here, the probability of each renewal event can be derived according to the Sections 3 and 4.
5.3. Parameter estimation
Before estimating the parameters, the distribution type of these three stages should be determined, which can be viewed as the foundation of maintenance decision. As the Weibull distribution is one of the most commonly used distributions in reliability [26], these three stages are assumed to follow Weibull distributions and the pdf of the $n$th ($n=$1, 2 and 3) stage is:
where ${a}_{n}$, ${b}_{n}$ are scale parameter and shape parameter respectively.
The useful information can be derived from the Table 1:
1) Total 5 times of minor defect renewal are occurred at the time ${t}_{j}=$ (219, 45, 199, 77, 57) respectively and the last inspection time before renewal is (194, 21, 149, 54, 28) respectively;
2) Total 3 times of severe defect renewal are occurred at the time ${t}_{n}=$ (15, 58, 13) respectively and the last inspection time before renewal is (0, 29, 0) respectively;
3) Total 2 times of age replacement renewal are occurred at the time ${T}_{u}=$(38, 41) respectively and the last inspection time before renewal is (17, 26) respectively;
4) Total 1 time of failure renewal are occurred at the time ${x}_{q}=$72 respectively and the last inspection time before renewal is 42 respectively;
5) The bearings are normal at the ending time 283.
The log likehood function can be derived according to above information:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\mathrm{l}\mathrm{o}\mathrm{g}\left[{P}_{n}\left(283\right)\right].$
Through maximizing the log $L$, the distribution parameters of these three stages can be estimated as Table 2.
On the other hand, the downtime parameters are given by the experimence of maintenance operators in the industry, as shown in Table 3.
Table 2. The estimation value of lifetime distribution parameters
${a}_{1}$

${b}_{1}$

${a}_{2}$

${b}_{2}$

${a}_{3}$

${b}_{3}$

0.0154

1.156

0.0174

1.758

0.0182

2.973

Table 3. Downtime parameters
${D}_{s}$

${D}_{r}$

${D}_{pm}$

${D}_{ps}$

${D}_{f}$

$r$

1

3

5

10

50

0.6

5.4. Maintenance decisionmaking
Depending on the given distribution and downtime parameters, we can determine the optimal decision variables and analyze the outputs of the proposed models.
As for the availability models with imperfect inspection in Section 3, Figs. 11 and 12 show the results of the two proposed models in terms of the inspection and replacement intervals. In model 1 (halving the subsequent inspection intervals when the minor defect is identified with imperfet inspection), the optimal intervals are ${t}^{*}=$37, ${N}^{*}=$2 with maximum longrun availability ${A}^{*}$(37, 2)$=$0.9269. However, the optimal intervals are ${t}^{*}=$39, ${N}^{*}=$ 2 with maximum longrun availability ${A}^{*}$(39, 2) $=$0.9236 in model 2 (repair immediately when the minor defect is identified with imperfet inspection).
On the other hand, the results of the availability models with perfect inspection case in terms of the two decision variables are presented in Figs. 13 and 14. The optimal intervals are ${t}^{*}=$40, ${N}^{*}=$ 2 with maximum long run availability ${A}^{*}$(40, 2)$=$0.9248 in model 3 (halving the subsequent inspection intervals when the minor defect is identified with perfet inspection). But the optimal solutions are ${t}^{*}=$48, ${N}^{*}=$ 2 with maximum long run availability ${A}^{*}$(48, 2) $=$0.9208 in model 4 (repair immediately when the minor defect is identified with perfet inspection).
Fig. 11. The longrun availability of the model 1 with imperfect inspection
Fig. 12. The longrun availability of the model 2 with imperfect inspection
Fig. 13. The longrun availability of the model 3 with perfect inspection
Fig. 15 presents the results of the four proposed models for comparison. As all of the optimal $N$ of four models are 2, Fig. 15 gives the availability values with constant $N=$2 and different inspection intervals $t$.
From the results, we can see that the optimal availability value of model 1 is larger than model 2 under the imperfect inspection case. The same results can be derived under the perfect inspection case. It shows that no matter the perfect inspection or not, the models 1 and 3 is optimal compared with models 2 and 4 with the given parameters from the numerical examples. The results mean that shortening the inspection interval is better than replacing immediately regarding the maximum longrun availability under the joint policy of inspection and age replacement. The results also shows the efficiency and importance of delaying the maintenance action based on the threestage failure process. Therefore, when the minor defective stage is identified, the optimal action is to shortening the subsequent inspection interval for these systems subject to threestage failure process. This conclusion can help the managers make scientific maintenance decisions.
Fig. 14. The longrun availability of the model 4 with perfect inspection
Fig. 15. The comparison of the four models
6. Conclusions
Regular inspections and age based preventive replacement are common activities in industrial maintenance actions. In this paper, a joint optimal policy of inspection and age based replacement for a single component based on a threestage failure process is proposed aiming at jointly optimizing the inspection and replacement intervals. The objective is the maximum longrun availability. Three states of system are divided before failure by the concept of threestage failure process, namely normal, minor defective and severe defective. Age based replacement policy is also considered due to its wide applications in industry. The system is replaced when it fails or reaches certain age. Once the severe defective stage is identified, maintenance action is carried out. Two different models are proposed according to the different measures of identifying the minor defective stage. A case study is presented to demonstrate the efficiency of the proposed models. The results show that delaying the maintenance action is better than replacing immediately based on the threestage failure process with this joint maintenance policy.
Further research with the concept of threestage failure process can be developed such as: 1) The complex or multicomponent system which are more commom in practice can be considered, 2) The combination of conditionbased maintenance and the concept of threestage failure process should be considered as the identification of the three stages will be easier with conditionbased maintenace policy. These issues will be researched in the future.
References
 Zhang Z. Q., Wu S., Li B. F., Lee S. Optimal maintenance policy for multicomponent systems under Markovian environment changes. Expert Systems with Applications, Vol. 40, 2013, p. 73917399. [Search CrossRef]
 Adriaan V. H., Liliane P. A dynamic predictive maintenance policy for complex multicomponent systems. Reliability Engineering and System Safety, Vol. 120, 2013, p. 3950. [Search CrossRef]
 Wang W. An overview of the recent advances in delaytimebased maintenance modelling. Reliability Engineering and System Safety, Vol. 106, 2012, p. 165178. [Search CrossRef]
 Koochaki J., Bokhorst J., Wortmann H., Klingenberg W. Evaluating condition based maintenance effectiveness for two processes in series. Journal of Quality in Maintenance Engineering, Vol. 17, 2011, p. 398414. [Search CrossRef]
 Li Z., He Z. J., Zi Y. Y., Chen X. F. Bearing condition monitoring based on shock pulse method and improved redundant lifting scheme. Mathematics and Computers in Simulation, Vol. 79, 2008, p. 318338. [Search CrossRef]
 Tandon N., Yadava G. S., Ramakrishna K. M. A comparison of some condition monitoring techniques for the detection of defect in induction motor ball bearings. Mechanical Systems and Signal Processing, Vol. 21, Issue 1, 2007, p. 244256. [Search CrossRef]
 Yang R. F., Kang J. S., Zhang X. H., Teng H. Z., Li H. P. Research on operating condition effect on the shock pulse method. Telkomnika Journal of Electrical Engineering, Vol. 12, Issue 7, 2014, p. 53925398. [Search CrossRef]
 Bartholomew M., Christianson B., Zuo M. J. Optimizing preventive maintenance models. Computational Optimization and Applications, Vol. 35, 2006, p. 261279. [Search CrossRef]
 Nicola R. P., Dekker R. Optimal Maintenance of MultiComponent Systems: a Review. Murthy DNP, Kobbacy AKS, Editors, Complex System Maintenance Handbook. Springer, Amsterdam, 2008. [Search CrossRef]
 Wang W. Models of inspection, routine service, and replacement for a serviceable onecomponent system. Reliability Engineering and System Safety, Vol. 116, 2013, p. 5763. [Search CrossRef]
 Cavalcante C. A. V., Scarf P. A., Almeida A. T. A study of a twophase inspection policy for a preparedness system with a defective state and heterogeneous lifetime. Reliability Engineering and System Safety, Vol. 96, Issue 6, 2011, p. 627635. [Search CrossRef]
 Cunningham A., Wang W., Zio W., Allanson E., Wall D., Wang J. A. Application of delay time analysis via Monte Carlo simulation. Journal of Marine Engineering and Technology, Vol. 10, Issue 3, 2011, p. 5772. [Search CrossRef]
 Lu W. Y., Wang W. Modelling Preventive maintenance based on the delay time concept in the context of a case study. Journal Maintenance and Reliability, Vol. 3, 2011, p. 410. [Search CrossRef]
 Wang W., Scarf P. A., Smith M. A. J. On the application of a model of conditionbased maintenance. Journal of the Operational Research Society, Vol. 51, 2000, p. 12181227. [Search CrossRef]
 Wang W. An inspection model based on a threestage failure process. Reliability Engineering and System Safety, Vol. 96, 2011, p. 838848. [Search CrossRef]
 Wang W., Zhao F., Peng R. A preventive maintenance model with a twolevel inspection policy based on a threestage failure process. Reliability Engineering and System Safety, Vol. 121, 2014, p. 207220. [Search CrossRef]
 Yang R. F., Yan Z. W., Kang J. S. An inspection maintenance model based on a threestage failure process with imperfect maintenance via Monte Carlo simulation. International Journal of Systems Assurance Engineering and Management, 2014. [Search CrossRef]
 Yang R. F., Zhao F., Kang J. S., Li H. P., Teng H. Z. Inspection optimization model with imperfect maintenance based on a threestage failure process. Maintenance and Reliability, Vol. 17, Issue 2, 2015, p. 165173. [Search CrossRef]
 Yang R. F., Zhao F., Kang J. S., Zhang X. H. An inspection optimization model based on a threestage failure process. International Journal of Performability Engineering, Vol. 10, Issue 7, 2014, p. 775779. [Search CrossRef]
 Scarf P. A., Cavalcante C. A. V., Dwight R. A., Gordon P. An agebased inspection and replacement policy for heterogeneous components. IEEE Transactions on Reliability, Vol. 58, Issue 4, 2009, p. 641648. [Search CrossRef]
 Sharareh T., Dragan B. Optimal inspection of a complex system subject to periodic and opportunistic inspections and preventive replacements. European Journal of Operational Research, Vol. 220, 2012, p. 649660. [Search CrossRef]
 Huynh K. T., Castro I. T., Barros A., Bérenguer C. Modeling agebased maintenance strategies with minimal repairs for systems subject to competing failure modes due to degradation and shocks. European Journal of Operational Research, Vol. 218, Issue 1, 2012, p. 140151. [Search CrossRef]
 Scarf P. A., Cavalcante C. A. V. Modelling quality in replacement and inspection maintenance. International Journal of Production Economics, Vol. 135, Issue 1, 2012, p. 372381. [Search CrossRef]
 Wang W. Delay Time Modelling. Murthy DNP, Kobbacy AKS, Editors. Complex System Maintenance Handbook. Springer, Amsterdam, 2008, p. 345370. [Search CrossRef]
 Ross S. M. Introduction to Probability Models. Eighth Edition, Academic Press, 2003. [Search CrossRef]
 Murthy D. N. P., Xie M., Jiang R. Weibull Models. Willey, New Jersey, 2004. [Search CrossRef]