# Optimization Problems of Costs in a MultiState System with Imperfect Maintenance

Authors: Soufiane Gasmi; Nizar Mannai
Journal: Journal of Quality and Reliability Engineering
 (2014)

---

## Abstract

We consider optimal replacement policies with periodic imperfect maintenance actions and minimal repairs. The multistate system is minimally repaired at failure and imperfect maintenance actions are regularly carried out for preventive maintenance. The discrete modified Weibull distribution is introduced and some cost functions applied to this distribution are defined in order to be minimized. Moreover, we assume that the costs of preventive maintenance depend on the degree of repair via a Kijima type 2 model. For illustrative purpose, the obtained results are applied on sets of simulated data.


---

## Body

## 1. Introduction

In reliability theory, discrete failure data arise in several common situations. It is sometimes impossible to measure the life length of items with continuous scale. A discrete distribution is appropriate when we have a multistate system and the number of states prior to preventive maintenance is observed. A few results on discrete life distributions are introduced in the literature (Salvia and Bollinger [1], Kalbfleisch and Prentice [2], Shmueli et al. [3], and Nakagawa and Osaki [4]).Recently, in reliability, much attention has been paid to the evaluation and application of multistate systems (MSS). We note that the MMS model provides a more flexible tool for modeling engineering systems in real life. This model is at first introduced in Barlow and Wu [5] and El-Neveihi et al. [6]. Some recent contributions on the subject are the works of (Zuo and Tian [7]; Li and Zuo [8]; Lisnianski and Ding [9]; Lisnianski and Frenkel [10]).The most commonly used models for the failure process of a repairable system are known as perfect repair or as good as new and as minimal repair or as bad as old.In the perfect repair case, each repair is perfect and leaves the system as if it were new; hence we obtain a renewal processes (RP). In the minimal repair case, each repair leaves the system in the same state as it was before failure, so we obtain nonhomogeneous Poisson processes (NHPP). It is well known in practice that the reality is between these two extreme cases. The repair may not yield a functioning item which is as good as new and the minimal repair assumption seems to be too pessimistic in repair strategies. From this it is seen that the imperfect repair is of great signification in practice.Furthermore, the imperfect maintenance is a broader subject, in this way many authors have considered other aspects associated with the imperfect maintenance that corresponds to errors when doing the maintenance activities. In fact, when the maintainer makes a mistake it may correspond not only in terms of cost, but also in terms of reliability.Currently, Berrade et al. [11] and Berrade et al. [12] have supposed that inspections are subject to error, and false positives or false negatives are possible. A false positive occurs when the inspection says the system is defective when in fact it is good. Consequently, the system is unnecessarily perfectly replaced and this implies an additional cost to the system owner. A false negative occurs when the inspection says the system is good when in fact it is defective. In this case, the system reliability will be reduced.In this paper, we are concerned with the modeling of an imperfect maintenance model; that is, the impact of preventive maintenance is between the two boundary cases: (a) minimal (as bad as old) and (b) perfect (as good as new).Further, it is assumed that costs of preventive maintenance are not constant but depend on the degree of repair via a virtual age process.A (MSS) withN states is considered. On failure, corrective maintenance actions are minimal; that means that the item is minimally repaired at each failure and it will have the same state as just before failure. Therefore, we consider the inspection as a possible action from the maintenance. In this way we generalize the case that after preventive maintenance actions, only perfect repairs are considered.In order to model the optimal maintenance actions, a new discrete modified Weibull distribution is introduced. At last some possible cost functions are proposed and optimal maintenance strategies are discussed.


## 2. Modeling the Multistate System

Most operating units (items) are repaired or replaced when they have failed. However, it may need much time and high cost to repair a failed item, so it is essential to maintain it in order to prevent failures and to reduce costs.After failure the maintenance is called corrective maintenance (CM) and before failure it is called preventive maintenance (PM). Some PM policies and their optimization problems for shock and damage models are summarized in Nakagawa [13].In practice, PM is used to lengthen the useful lifetime of items and to decrease average running cost by minimizing the occurrence of failures. Here we can note that most of the models concerning the modeling of repairable systems identify the minimal repair and the imperfect repair actions. Naturally, this popular assumption is a very unreal one.In the following, a multistate system withN states is considered. To simplify the notation and to use results from the imperfect maintenance models with discrete lifetime distribution, a simple time scale is suitable (Kahle [14]). Therefore, a time scale is introduced, and it will be assumed that in time 1 the item is in state one, in time 2 the item is in state two, and so forth. We remark that this time scale becomes more simple and can be generalized to any time scale. We studied the case that PM actions are regularly carried out. Furthermore, failures between PM actions are minimally repaired. The PM actions reset the failure rate proportional to the virtual age of the item. We assume that the failure rate at state, k=1,2,…,N, is increasing in time t. In the case of multistate system, the effect of the PM actions is to reset the item in a lower state that means a younger state. We denote the state after PM actions by state s.It is supposed that the reduction of the failure rate is done by using the concept of virtual age reduction introduced by Kijima [15] and Kijima et al. [16]. We will particularly use the Kijima’s type 2 imperfect repair model to describe the PM actions. We note that the virtual age of the item at time tk, following the repair, is defined by Vk=ζk(Vk-1+tk-tk-1), k≥1, where tk is time of the kth sojourn and ζk(0≤ζk≤1)is the degree of repair at that time. We assume that t0=0 and V0=0. After repairing the failure intensity during the (k+1)th sojourn is determined by λk+1(t)=λ(t-tk+Vk), tk≤t<tk+1, k≥0.The process defined by(t,ζk,k=1,…,N)=t-tk+Vk,tk≤t<tk+1, k≥0, is called the virtual age process. We can note that in Gasmi et al. [17] a generalized Kijima type 2 model was considered, where repairs give an additional impact made by using a major repair.


## 3. The Discrete Modified Weibull Distribution

In some situations, discrete lifetime distributions are appropriate to model lifetimes. The first discrete reliability distribution has been defined by Nakagawa and Osaki [4]. Furthermore, Kahle [14] studied some discrete lifetime distributions like the discrete uniform distribution, the Poisson distribution, and the truncated discrete Weibull distribution. In this paper the discrete modified Weibull distribution (DMWD) will be introduced.We note that the modified Weibull distribution (MWD) is recently proposed by Sarhan and Zaindin [18]. This distribution is a general form for some most commonly used distributions in survival analysis such as exponential, Rayleigh, linear failure rate, and Weibull distribution. We will use the notation MWD (α,β,γ) to denote the modified Weibull distribution with the scale parameter α and the two shape parameters β and γ.LetX be MWD (α,β,γ) with density function:
(1)f(x,θ)=(α+βγxγ-1)exp⁡(-αx-βxγ),x>0,
where θ=(α,β,γ), γ>0, α,β≥0, such that α+β>0.For the discretization of this distribution, we will put the probability mass of the interval(k-1,k] into the point k; that is,
(2)P(X=k)=∫k-1kf(u)du=exp⁡(-α(k-1)-β(k-1)γ)-exp⁡(-αk-βkγ),k=1,2,….
The reliability function is given by the probability that the system is still alive at time k:
(3)R(k)=P(X>k)=∑i=k+1∞P(X=i);
therefore we have P(X=k)=R(k-1)-R(k).In the case of the DMWD we obtainR(k)=exp⁡(-αk-βkγ), k=1,2,….The failure rate of a discrete distribution has been defined by Barlow et al. [19]:
(4)λ(k)=R(k-1)-R(k)R(k-1),k=1,2,….
The definition of the failure rate as in (4) has many disadvantages. We cite that in the discrete case, the failure rate is bounded by 1 while most of continuous reliability models have unbounded failure rates. Moreover, λ(k) is not a convex function. This restricts the interpretation of the failure rate in praxis. Besides, λ(k) can never grow linearly or exponentially. Xie et al. [20] have shown that the cumulative hazard function, denoted by H(k), verified the relation H(k)=-ln⁡R(k)≠∑i=1kλ(i).They proposed a new redefined failure rate function for discrete distributions, denoted by(5)r(k)=-ln⁡R(k)R(k-1)=ln⁡R(k-1)R(k).
Consequently, the above mentioned problems are resolved.Both failure ratesλ(k) and r(k) are linked by r(k)=-ln⁡(1-λ(k)). On the other hand, we can see that the two failure rate definitions have the same monotonicity property.Xie et al. [20] proved that if one is increasing/decreasing, the other is also increasing/decreasing. Therefore, the redefined failure rate function is more appropriate and is used for a much useful discrete reliability functions.For the DMWD(α,β,γ), we have then
(6)r(k)=(αk+βkγ)-(α(k-1)+β(k-1)γ),hhhhhhhhhhhhhhhhhhhhhhlhk=1,2,….
The cumulative hazard function is defined by
(7)H(k)=∑i=1kr(i).
We obtain then the same form as for the continuous MWD.                         (8)H(k)=(αk+βkγ),k=1,2,….
The rth moment of the DMWD is given by E(Xr)=∑i=1∞ir[R(i-1)-R(i)]. If we put r equal to 1, we obtain E(X)=∑i=0∞exp⁡(-αi-βiγ) and the mean time to failure MTTF of the DMWD has then the form MTTF=∑i=0∞exp⁡(-αi-βiγ). We can remark that in the case of β=0 we obtain the discrete exponential distribution with parameter α. By putting α=0 and γ=2 we obtain the discrete Rayleigh distribution with parameter β. If γ=2, a discrete linear failure rate with parameters α and β is obtained. Finally if α=0 the Weibull distribution with parameters β and γ will be derived.Figures1 and 2 show, respectively, different patterns of the discrete failure rate r(k) and the cumulative hazard function H(k) in the case of linear failure rate, Rayleigh distribution, Weibull distribution, and the modified Weibull distribution. We note that the parameters for all these distributions are chosen so that the expectation is equal to 5 and N=20.Figure 1
Discrete failure ratesr(k).Figure 2
Cumulative hazard functionsH(k).Assumption 1.
We suppose that the shape parameterγ≥2, this implies that λ(k) is increasing and H(k) is convex.


## 4. The Maintenance Policy

A sequential imperfect preventive maintenance policy was proposed in Nakagawa [21] and Chen et al. [22]. This policy supposes that PM is undertaken at fixed intervals.We consider sequential PM policies for the system, where the distance between PM actions is always equal tov. Furthermore, all CM actions are minimal, that means that the item is minimally repaired at each failure and after repair the item has the same state as just before failure. It is assumed that the duration for CM and PM is negligible.Following the strategy of Kahle [14], we assume that the multistate system has N states. Let s be the state of the item after such PM actions and the next PM executed v states later. Every time when the item reaches the state (s+v) the PM action resets it to the state s. In this case the following pairs of v and s are possible:
(9)v=N-1,s=1,v=N-2,s=1ors=2,…⋮v=1,s=1,2,…,ors=N-1.
Let ζ=s/(s+v) be the degree of repair following the Kijima type 2 model. We have assumed that 0<ζ<1 and we eliminated the two cases:(i)
ζ=0 which correspond to perfect repair action that involves that the item is reset to state s=0;(ii)
ζ=1 which correspond to minimal repair action and engender v=0.Now, our purpose is to obtain the pairs v and s minimizing the costs of PM and CM per time unit. In numerical examples, we investigate the effect of model parameters to the expected cost.


## 5. The Optimal Expected Cost Function

We consider maintenance policies for repairable systems in which the economical approach is our main objective. In our model each PM action reduces the age of the system tos (state s) and every CM is minimal. The idea in some papers (Nakagawa [23]) assuming that each PM action has the same costs is not realistic. In Kahle [14] a new cost function which describes the cost of the PM action according to the degree of repair is introduced, where PM actions are undertaken at fixed times v, 2v, ….LetCCM and CPM be, respectively, the cost of a corrective maintenance action and the cost of preventive maintenance action. The cost function per time unit is defined as follows:
(10)Ctot(s,v)=CPM(s,v)+CCM(H(s+v)-H(s))v.
Our main objective is to minimize this function with respect to s and v. Now, we will introduce new defined cost functions applied to the DMWD.


### 5.1. The Costs of PM Action Proportional to the State before Repair

We assume that the cost functionCPM depends on the state just before the preventive maintenance action:
(11)CPM(s,v)=Cs(1s+v)δ,
where δ∈IR+ has the role to make this cost function more flexible.Let the ratio of costsCCM/Cs=2, δ=0.5, and N=20. The cost function per time unit is plotted for the DMWD in Figure 3. We remark that the points will be connected with lines to make more visibility and that for different values of v. Our cost function is convex; it follows then the uniqueness of the minimum.Figure 3
Cost function for DMWD (α=0.03, β=6.1 10−3, and γ=3).In Table1, we will present for different values of δ and different ratio costs the optimal states before the PM (s+v) and after the PM (s).Table 1
Optimal values fors and v and corresponding optimal costs Ctot.
DLFRD(α, β)
DRD(β)
DWD(β, γ)
DMWD(α, β, γ)
α = 0.01 β = 0.0366
β = 0.0388
β = 0.0078 γ = 3
α = 0.03   β = 6.1 10−3  γ = 3
C
CM
/
C
s
=
3
δ = 0.5
s = 1v = 3  Ctot = 0.7457
s = 1v = 3Ctot = 0.7485
s = 1v = 2Ctot = 0.5937
s = 1v = 2Ctot = 0.6182
δ = 1
s = 1v = 2Ctot = 0.6359
s = 1 v = 2Ctot = 0.6321
s = 1v = 2Ctot = 0.4717
s = 1v = 2Ctot = 0.4962
δ = 2
s = 1v = 2Ctot = 0.5248
s = 1v = 2Ctot = 0.5210
s = 1v = 2Ctot = 0.3605
s = 1v = 2Ctot = 0.3851
δ = 1/8
s = 1v = 3Ctot = 0.8593
s = 1v = 3Ctot = 0.8621
s = 1v = 2Ctot = 0.7408
s = 1v = 3Ctot = 0.7572
δ = 8
s = 1v = 1Ctot = 0.3633
s = 1v = 1Ctot = 0.3530
s = 1v = 1Ctot = 0.1681
s = 1v = 1Ctot = 0.2229
C
CM
/
C
s
=
2
δ = 0.5
s = 1v = 3Ctot = 0.5527
s = 1v = 3Ctot = 0.5545
s = 1v = 2Ctot = 0.4920
s = 1v = 3Ctot = 0.4846
δ = 1
s = 1v = 3Ctot = 0.4694
s = 1v = 3Ctot = 0.4712
s = 1v = 2Ctot = 0.3700
s = 1v = 2Ctot = 0.3863
δ = 2
s = 1v = 2Ctot = 0.3684
s = 1v = 2Ctot = 0.3659
s = 1v = 2Ctot = 0.2589
s = 1v = 2Ctot = 0.2752
δ = 1/8
s = 1v = 4Ctot = 0.6637
s = 1v = 3Ctot = 0.6682
s = 1v = 3Ctot = 0.6087
s = 1v = 3Ctot = 0.5982
δ = 8
s = 1v = 1Ctot = 0.2435
s = 1v = 1Ctot = 0.2366
s = 1v = 1Ctot = 0.1134
s = 1v = 1Ctot = 0.1499
C
CM
/
C
s
=
1
δ = 0.5
s = 1v = 4Ctot = 0.3414
s = 1v = 4Ctot = 0.3445
s = 1v = 3Ctot = 0.3309
s = 1v = 3Ctot = 0.3256
δ = 1
s = 1v = 3Ctot = 0.2763
s = 1v = 3Ctot = 0.2773
s = 1v = 3Ctot = 0.2476
s = 1v = 3Ctot = 0.2423
δ = 2
s = 1v = 2Ctot = 0.2120
s = 1v = 2Ctot = 0.2107
s = 1v = 2Ctot = 0.1572
s = 1v = 2Ctot = 0.1654
δ = 1/8
s = 1v = 5Ctot = 0.4261
s = 1v = 5Ctot = 0.4314
s = 1v = 3Ctot = 0.4445
s = 1v = 4Ctot = 0.4248
δ = 8
s = 1v = 1Ctot = 0.1237
s = 1v = 1Ctot = 0.1203
s = 1v = 1Ctot = 0.0586
s = 1v = 1Ctot = 0.0769We have used the following discrete distributions: linear failure rate distribution (DLFRD), Rayleigh distribution (DRD), Weibull distribution (DWD), and the modified Weibull distribution (DMWD).We also compared the cost function for these distributions and we noted that the parameters were chosen so that the expectation was the same for all distributions. Let this expectation be 5.


### 5.2. The Costs of PM Action Proportional to the Impact of Repair

In this case, we will assume that the cost functionCPM depends only on the state after the PM action:
(12)CPM(s,v)=CI·(1s)δ,
where δ∈IR+·Let the ratio of costsCPM/CI=2, δ=2, and N=20. In Figure 4, each line corresponds to a different value of v. This figure illustrates the cost function per time unit for the discrete modified Weibull distribution. As earlier, to make more visibility, the points are connected with lines. The states (s+v) and (s), corresponding to the unique optimum for the cost function, are given in Table 2.Table 2
Optimal values fors and v and corresponding optimal costs Ctot.
DLFRD(α, β)
DRD(β)
DWD(β, γ)
DMWD(α, β, γ)
α  =  0.01 β  =  0.0075
β  = 0.0086
β  = 8.31 10−4  γ = 3
α = 0.03    β  = 4.76  10−4  γ  =  3
C
PM
/
C
I
=
3
δ = 0.5
s  =  2v  = 6Ctot = 0.3730
s  =  1v  = 6Ctot = 0.3727
s  =  2v  = 4Ctot = 0.3065
s  =  2v  = 5Ctot = 0.3273
δ = 1
s  =  2v  = 5Ctot = 0.3326
s  =  2v  = 4Ctot = 0.3310
s  =  2v  = 4Ctot = 0.2547
s  =  2v  = 5Ctot = 0.2858
δ = 2
s  =  3v  = 2Ctot = 0.2657
s  =  3v  = 2Ctot = 0.2615
s  =  3v  = 2Ctot = 0.1778
s  =  3v  = 2Ctot = 0.2157
δ = 1/8
s  =  1v  = 7Ctot = 0.3755
s  =  1v  = 6Ctot = 0.3727
s  =  1v  = 5Ctot = 0.3073
s  =  1v  = 7Ctot = 0.3373
δ = 8
s  =  2v  = 1Ctot = 0.1465
s  =  2v  = 1Ctot = 0.1327
s  =  2v  = 1Ctot = 0.0513
s  =  2v  = 1Ctot = 0.1211
C
PM
/
C
I
=
2
δ = 0.5
s  =  2v  = 7Ctot = 0.2861
s  =  2v  = 6Ctot = 0.2895
s  =  2v  = 5Ctot = 0.2529
s  =  2v  = 6Ctot = 0.2580
δ = 1
s  =  3v  = 5Ctot = 0.2518
s  =  2v  = 5Ctot = 0.2545
s  =  2v  = 5Ctot = 0.2114
s  =  3v  = 4Ctot = 0.2187
δ = 2
s  =  3v  = 3Ctot = 0.1921
s  =  3v  = 3Ctot = 0.1915
s  =  3v  = 2Ctot = 0.1371
s  =  3v  = 3Ctot = 0.1571
δ = 1/8
s  =  1v  = 8Ctot = 0.2951
s  =  1v  = 8Ctot = 0.2967
s  =  1v  = 6Ctot = 0.2615
s  =  1v  = 8Ctot = 0.2718
δ = 8
s  =  2v  = 1Ctot = 0.0990
s  =  2v  = 1Ctot = 0.0897
s  =  2v  = 1Ctot = 0.0355
s  =  2v  = 1Ctot = 0.0820
C
PM
/
C
I
=
1
δ = 0.5
s  =  2v  = 10Ctot = 0.1858
s  = 2 v  = 9Ctot = 0.1901
s  =  2v  = 7Ctot = 0.1867
s  =  2v  =  8Ctot = 0.1775
δ = 1
s  =  3v  = 7Ctot = 0.1552
s  =  3v  = 6Ctot = 0.1586
s  =  3v  = 5Ctot = 0.1473
s  =  3v  = 6Ctot = 0.1413
δ = 2
s  =  3v  = 4Ctot = 0.1128
s  =  3v  = 4Ctot = 0.1136
s  =  3v  = 3Ctot = 0.0894
s  =  4v  = 3Ctot = 0.0952
δ = 1/8
s  =  1v  = 12Ctot = 0.1984
s  =  1v  = 11Ctot = 0.2025
s  =  1v  = 8Ctot = 0.2007
s  =  1v  = 10Ctot = 0.1934
δ = 8
s  =  2v  = 1Ctot = 0.0514
s  =  2v  = 1Ctot = 0.0468
s  =  2v  = 1Ctot = 0.0197
s  =  2v  = 1Ctot = 0.0430Figure 4
Cost function for DMWD (α=0.03, β=4.76 10−4, and γ=3).The cost functions for distributions used in Section5.1 are compared. We note that all distributions have the same expectation 10.


### 5.3. The Costs of PM Action Proportional to the Degree of Repair

Now, we consider a new cost function per time unit which depends on the degree of repairζ.We know that for the Kijima type 2 model, the degree of repairζ corresponds to ζ=s/(s+v).LetCPM(s,v)=CPR·(1-(s/(s+v))exp⁡(s/(s+v)-1))δ be the new cost function where CPR are the perfect repair costs. This function can be well interpreted.Remark 2.
(a) For perfect repair,ζ=0, implying that (s=0), we obtain CPM(s,v)=CPR.
(b) For minimal repair,ζ=1, involving that (v=0), we obtain CPM(s,v)=0.
(c) Ifδ>1, the cost function CPM(s,v) is convex and with increasing δ the costs of large repair are relatively small.
(d) Ifδ<1, the cost function CPM(s,v) is concave and with decreasing δ, the costs of relatively small repairs are close to that of perfect repair CPR.Remark 3.
(i) Letζ=0.1; Table 3 shows the costs of large repair for different values of δ.
(ii) Letζ=0.9; Table 4 tells us about the costs of small repair for different values of δ.Table 3
Costs of large repairs.
δ
2
3
4
8
20
C
PM
(
s
,
v
)
0.9203
*
C
PR
0.8829
*
C
PR
0.8470
*
C
PR
0.7174
*
C
PR
0.4359
*
C
PRTable 4
Costs of small repairs.
δ
1/2
1/3
1/4
1/8
1/20
C
PM
(
s
,
v
)
0.4308
*
C
PR
0.5704
*
C
PR
0.6564
*
C
PR
0.8101
*
C
PR
0.9192
*
C
PRThe resulting cost function per time unit has then the following form:(13)Ctot(s,v)=(CPR·(1-(ss+v)exp⁡(ss+v-1))δnnl(1-(ss+v)exp⁡(ss+v-1))δi+CCM(H(s+v)-H(s)))×(v)-1.
This function should be minimized with respect to s and v.Similar to Sections5.1 and 5.2, Figure 5 illustrates the cost function defined in (13) for the DMWD with N states, N=20, CCM/CI=3, and δ=1/2.Figure 5
Cost function for DMWD (α=0.03, β=4.76 10−4, and γ=3); δ=0.5.Ifδ=2 and CCM/CI=3, we obtain Figure 6. For each value of v, it can be seen that CPM(s,v) has a unique minimum. The obtained results are given in Table 5.Table 5
Optimal values fors and v and corresponding optimal costs Ctot.
DLFRD(α, β)
DRD(β)
DWD(β, γ)
DMWD(α, β, γ)
α = 0.01 β  =  5.4  10−3
β  =  8.6  10−3
β  = 8.39  10−4  γ = 3
α = 0.03   β  = 4.76  10−4  γ = 3
C
CM
/
C
I
=
3
δ = 0.5
s  =  1v  = 8Ctot = 0.3755
s  =  1v  = 6Ctot = 0.3675
s  =  1v  = 5Ctot = 0.3009
s  =  1v  = 6Ctot = 0.3331
δ = 1
s  =  1v  = 7Ctot = 0.3724
s  =  1v  = 6Ctot = 0.3626
s  =  1v  = 5Ctot = 0.2938
s  =  1v  = 6Ctot = 0.3281
δ = 2
s  =  4v  = 1Ctot = 0.3461
s  =  4v  = 1Ctot = 0.3508
s  =  3v  = 1Ctot = 0.2661
s  =  4v  = 1Ctot = 0.2963
δ = 1/8
s  =  1v  = 8Ctot = 0.3776
s  =  1v  = 6Ctot = 0.3714
s  =  1v  = 5Ctot = 0.3064
s  =  1v  = 7Ctot = 0.3363
δ = 8
s  =  2v  = 1Ctot = 0.1772
s  =  1v  = 1Ctot = 0.1328
s  =  2v  = 1Ctot = 0.0534
s  =  2v  = 1Ctot = 0.1227
C
CM
/
C
I
=
2
δ = 0.5
s  = 1 v  = 9Ctot = 0.2886
s  =  1v  = 7Ctot = 0.2936
s  =  1v  = 6Ctot = 0.2572
s  =  1v  = 7Ctot = 0.2687
δ = 1
s  =  1v  = 9Ctot = 0.2864
s  =  1v  = 7Ctot = 0.2899
s  =  1v  = 6Ctot = 0.2522
s  =  1v  = 7Ctot = 0.2650
δ = 2
s  =  5v  = 1Ctot = 0.2666
s  =  4v  = 1Ctot = 0.2735
s  = 4 v  = 1Ctot = 0.2214
s  =  5v  = 1Ctot = 0.2336
δ = 1/8
s  =  1v  = 10Ctot = 0.2902
s  =  1v  = 8Ctot = 0.2959
s  =  1v  = 6Ctot = 0.2610
s  =  1v  = 8Ctot = 0.2711
δ = 8
s  =  2v  = 1Ctot = 0.1200
s  =  2v  = 1Ctot = 0.0914
s  =  2v  = 1Ctot = 0.0374
s  =  2v  = 1Ctot = 0.0837
C
CM
/
C
I
=
1
δ = 0.5
s  =  1v  = 13Ctot = 0.1875
s  =  1v  = 11Ctot = 0.2010
s  =  1v  = 8Ctot = 0.1985
s  =  1v  = 10Ctot = 0.1916
δ = 1
s  =  1v  = 13Ctot = 0.1864
s  =  1v  = 10Ctot = 0.1993
s  =  1v  = 8Ctot = 0.1957
s  =  1v  = 9Ctot = 0.1895
δ = 2
s  =  7v  = 1Ctot = 0.1636
s  =  6v  = 1Ctot = 0.1776
s  =  5v  = 1Ctot = 0.1632
s  =  6v  = 1Ctot = 0.1566
δ = 1/8
s  =  1v  = 14Ctot = 0.1883
s  =  1v  = 11Ctot = 0.2021
s  =  1v  = 8Ctot = 0.2006
s  =  1v  = 10Ctot = 0.1930
δ = 8
s  =  2v  = 1Ctot = 0.0628
s  =  2v  = 1Ctot = 0.0485
s  =  2v  = 1Ctot = 0.0215
s  =  2v  = 1Ctot = 0.0446Figure 6
Cost function for DMWD (α=0.03, β=4.76 10−4, and γ=3); δ=2.Once again, the choice of parameters of all used distributions in this table verified expectation 10.The summary of computation results is presented in Table5. Hence, we obtain the following remark.Remark 4.
(a) Ifδ≤1, then the state s=1 is always the optimal state after PM actions. This result is expected and can be interpreted as follows: the costs of PM actions are close to the costs of perfect repair CPR; this implies that the optimal policy will be a perfect repair strategy.
(b) Ifδ>1, then v is always equal to 1 and the optimal distance between maintenance is short. This result can be interpreted as follows: the optimum for the cost function is obtained for relatively small repairs.
(c) If the failure rate is DMWD or DWD the obtained results are quite similar.
(d) If the failure rate is that of DLFRD, then the optimal distance between PM actions is larger than that for the rest failure rates.
(e) Ifδ=8, then the optimal cost per time unit is very small.
(f) In the case of DMWD with parametersα=0.03, β=4.76 10−4, and γ=3, If CCM/CI=3 and δ=1/2, we obtain then the unique optimum for both s=1 and v=6. This implies the repair degree ζ=1/7 and the total costs are Ctot=0.3331.
(g) If we use the same parameters as in (f) and after replacingδ by 8, we obtain then the unique optimum in s=2 and v=1. It follows then that ζ=2/3 and Ctot=0.1227.


## 5.1. The Costs of PM Action Proportional to the State before Repair

We assume that the cost functionCPM depends on the state just before the preventive maintenance action:
(11)CPM(s,v)=Cs(1s+v)δ,
where δ∈IR+ has the role to make this cost function more flexible.Let the ratio of costsCCM/Cs=2, δ=0.5, and N=20. The cost function per time unit is plotted for the DMWD in Figure 3. We remark that the points will be connected with lines to make more visibility and that for different values of v. Our cost function is convex; it follows then the uniqueness of the minimum.Figure 3
Cost function for DMWD (α=0.03, β=6.1 10−3, and γ=3).In Table1, we will present for different values of δ and different ratio costs the optimal states before the PM (s+v) and after the PM (s).Table 1
Optimal values fors and v and corresponding optimal costs Ctot.
DLFRD(α, β)
DRD(β)
DWD(β, γ)
DMWD(α, β, γ)
α = 0.01 β = 0.0366
β = 0.0388
β = 0.0078 γ = 3
α = 0.03   β = 6.1 10−3  γ = 3
C
CM
/
C
s
=
3
δ = 0.5
s = 1v = 3  Ctot = 0.7457
s = 1v = 3Ctot = 0.7485
s = 1v = 2Ctot = 0.5937
s = 1v = 2Ctot = 0.6182
δ = 1
s = 1v = 2Ctot = 0.6359
s = 1 v = 2Ctot = 0.6321
s = 1v = 2Ctot = 0.4717
s = 1v = 2Ctot = 0.4962
δ = 2
s = 1v = 2Ctot = 0.5248
s = 1v = 2Ctot = 0.5210
s = 1v = 2Ctot = 0.3605
s = 1v = 2Ctot = 0.3851
δ = 1/8
s = 1v = 3Ctot = 0.8593
s = 1v = 3Ctot = 0.8621
s = 1v = 2Ctot = 0.7408
s = 1v = 3Ctot = 0.7572
δ = 8
s = 1v = 1Ctot = 0.3633
s = 1v = 1Ctot = 0.3530
s = 1v = 1Ctot = 0.1681
s = 1v = 1Ctot = 0.2229
C
CM
/
C
s
=
2
δ = 0.5
s = 1v = 3Ctot = 0.5527
s = 1v = 3Ctot = 0.5545
s = 1v = 2Ctot = 0.4920
s = 1v = 3Ctot = 0.4846
δ = 1
s = 1v = 3Ctot = 0.4694
s = 1v = 3Ctot = 0.4712
s = 1v = 2Ctot = 0.3700
s = 1v = 2Ctot = 0.3863
δ = 2
s = 1v = 2Ctot = 0.3684
s = 1v = 2Ctot = 0.3659
s = 1v = 2Ctot = 0.2589
s = 1v = 2Ctot = 0.2752
δ = 1/8
s = 1v = 4Ctot = 0.6637
s = 1v = 3Ctot = 0.6682
s = 1v = 3Ctot = 0.6087
s = 1v = 3Ctot = 0.5982
δ = 8
s = 1v = 1Ctot = 0.2435
s = 1v = 1Ctot = 0.2366
s = 1v = 1Ctot = 0.1134
s = 1v = 1Ctot = 0.1499
C
CM
/
C
s
=
1
δ = 0.5
s = 1v = 4Ctot = 0.3414
s = 1v = 4Ctot = 0.3445
s = 1v = 3Ctot = 0.3309
s = 1v = 3Ctot = 0.3256
δ = 1
s = 1v = 3Ctot = 0.2763
s = 1v = 3Ctot = 0.2773
s = 1v = 3Ctot = 0.2476
s = 1v = 3Ctot = 0.2423
δ = 2
s = 1v = 2Ctot = 0.2120
s = 1v = 2Ctot = 0.2107
s = 1v = 2Ctot = 0.1572
s = 1v = 2Ctot = 0.1654
δ = 1/8
s = 1v = 5Ctot = 0.4261
s = 1v = 5Ctot = 0.4314
s = 1v = 3Ctot = 0.4445
s = 1v = 4Ctot = 0.4248
δ = 8
s = 1v = 1Ctot = 0.1237
s = 1v = 1Ctot = 0.1203
s = 1v = 1Ctot = 0.0586
s = 1v = 1Ctot = 0.0769We have used the following discrete distributions: linear failure rate distribution (DLFRD), Rayleigh distribution (DRD), Weibull distribution (DWD), and the modified Weibull distribution (DMWD).We also compared the cost function for these distributions and we noted that the parameters were chosen so that the expectation was the same for all distributions. Let this expectation be 5.


## 5.2. The Costs of PM Action Proportional to the Impact of Repair

In this case, we will assume that the cost functionCPM depends only on the state after the PM action:
(12)CPM(s,v)=CI·(1s)δ,
where δ∈IR+·Let the ratio of costsCPM/CI=2, δ=2, and N=20. In Figure 4, each line corresponds to a different value of v. This figure illustrates the cost function per time unit for the discrete modified Weibull distribution. As earlier, to make more visibility, the points are connected with lines. The states (s+v) and (s), corresponding to the unique optimum for the cost function, are given in Table 2.Table 2
Optimal values fors and v and corresponding optimal costs Ctot.
DLFRD(α, β)
DRD(β)
DWD(β, γ)
DMWD(α, β, γ)
α  =  0.01 β  =  0.0075
β  = 0.0086
β  = 8.31 10−4  γ = 3
α = 0.03    β  = 4.76  10−4  γ  =  3
C
PM
/
C
I
=
3
δ = 0.5
s  =  2v  = 6Ctot = 0.3730
s  =  1v  = 6Ctot = 0.3727
s  =  2v  = 4Ctot = 0.3065
s  =  2v  = 5Ctot = 0.3273
δ = 1
s  =  2v  = 5Ctot = 0.3326
s  =  2v  = 4Ctot = 0.3310
s  =  2v  = 4Ctot = 0.2547
s  =  2v  = 5Ctot = 0.2858
δ = 2
s  =  3v  = 2Ctot = 0.2657
s  =  3v  = 2Ctot = 0.2615
s  =  3v  = 2Ctot = 0.1778
s  =  3v  = 2Ctot = 0.2157
δ = 1/8
s  =  1v  = 7Ctot = 0.3755
s  =  1v  = 6Ctot = 0.3727
s  =  1v  = 5Ctot = 0.3073
s  =  1v  = 7Ctot = 0.3373
δ = 8
s  =  2v  = 1Ctot = 0.1465
s  =  2v  = 1Ctot = 0.1327
s  =  2v  = 1Ctot = 0.0513
s  =  2v  = 1Ctot = 0.1211
C
PM
/
C
I
=
2
δ = 0.5
s  =  2v  = 7Ctot = 0.2861
s  =  2v  = 6Ctot = 0.2895
s  =  2v  = 5Ctot = 0.2529
s  =  2v  = 6Ctot = 0.2580
δ = 1
s  =  3v  = 5Ctot = 0.2518
s  =  2v  = 5Ctot = 0.2545
s  =  2v  = 5Ctot = 0.2114
s  =  3v  = 4Ctot = 0.2187
δ = 2
s  =  3v  = 3Ctot = 0.1921
s  =  3v  = 3Ctot = 0.1915
s  =  3v  = 2Ctot = 0.1371
s  =  3v  = 3Ctot = 0.1571
δ = 1/8
s  =  1v  = 8Ctot = 0.2951
s  =  1v  = 8Ctot = 0.2967
s  =  1v  = 6Ctot = 0.2615
s  =  1v  = 8Ctot = 0.2718
δ = 8
s  =  2v  = 1Ctot = 0.0990
s  =  2v  = 1Ctot = 0.0897
s  =  2v  = 1Ctot = 0.0355
s  =  2v  = 1Ctot = 0.0820
C
PM
/
C
I
=
1
δ = 0.5
s  =  2v  = 10Ctot = 0.1858
s  = 2 v  = 9Ctot = 0.1901
s  =  2v  = 7Ctot = 0.1867
s  =  2v  =  8Ctot = 0.1775
δ = 1
s  =  3v  = 7Ctot = 0.1552
s  =  3v  = 6Ctot = 0.1586
s  =  3v  = 5Ctot = 0.1473
s  =  3v  = 6Ctot = 0.1413
δ = 2
s  =  3v  = 4Ctot = 0.1128
s  =  3v  = 4Ctot = 0.1136
s  =  3v  = 3Ctot = 0.0894
s  =  4v  = 3Ctot = 0.0952
δ = 1/8
s  =  1v  = 12Ctot = 0.1984
s  =  1v  = 11Ctot = 0.2025
s  =  1v  = 8Ctot = 0.2007
s  =  1v  = 10Ctot = 0.1934
δ = 8
s  =  2v  = 1Ctot = 0.0514
s  =  2v  = 1Ctot = 0.0468
s  =  2v  = 1Ctot = 0.0197
s  =  2v  = 1Ctot = 0.0430Figure 4
Cost function for DMWD (α=0.03, β=4.76 10−4, and γ=3).The cost functions for distributions used in Section5.1 are compared. We note that all distributions have the same expectation 10.


## 5.3. The Costs of PM Action Proportional to the Degree of Repair

Now, we consider a new cost function per time unit which depends on the degree of repairζ.We know that for the Kijima type 2 model, the degree of repairζ corresponds to ζ=s/(s+v).LetCPM(s,v)=CPR·(1-(s/(s+v))exp⁡(s/(s+v)-1))δ be the new cost function where CPR are the perfect repair costs. This function can be well interpreted.Remark 2.
(a) For perfect repair,ζ=0, implying that (s=0), we obtain CPM(s,v)=CPR.
(b) For minimal repair,ζ=1, involving that (v=0), we obtain CPM(s,v)=0.
(c) Ifδ>1, the cost function CPM(s,v) is convex and with increasing δ the costs of large repair are relatively small.
(d) Ifδ<1, the cost function CPM(s,v) is concave and with decreasing δ, the costs of relatively small repairs are close to that of perfect repair CPR.Remark 3.
(i) Letζ=0.1; Table 3 shows the costs of large repair for different values of δ.
(ii) Letζ=0.9; Table 4 tells us about the costs of small repair for different values of δ.Table 3
Costs of large repairs.
δ
2
3
4
8
20
C
PM
(
s
,
v
)
0.9203
*
C
PR
0.8829
*
C
PR
0.8470
*
C
PR
0.7174
*
C
PR
0.4359
*
C
PRTable 4
Costs of small repairs.
δ
1/2
1/3
1/4
1/8
1/20
C
PM
(
s
,
v
)
0.4308
*
C
PR
0.5704
*
C
PR
0.6564
*
C
PR
0.8101
*
C
PR
0.9192
*
C
PRThe resulting cost function per time unit has then the following form:(13)Ctot(s,v)=(CPR·(1-(ss+v)exp⁡(ss+v-1))δnnl(1-(ss+v)exp⁡(ss+v-1))δi+CCM(H(s+v)-H(s)))×(v)-1.
This function should be minimized with respect to s and v.Similar to Sections5.1 and 5.2, Figure 5 illustrates the cost function defined in (13) for the DMWD with N states, N=20, CCM/CI=3, and δ=1/2.Figure 5
Cost function for DMWD (α=0.03, β=4.76 10−4, and γ=3); δ=0.5.Ifδ=2 and CCM/CI=3, we obtain Figure 6. For each value of v, it can be seen that CPM(s,v) has a unique minimum. The obtained results are given in Table 5.Table 5
Optimal values fors and v and corresponding optimal costs Ctot.
DLFRD(α, β)
DRD(β)
DWD(β, γ)
DMWD(α, β, γ)
α = 0.01 β  =  5.4  10−3
β  =  8.6  10−3
β  = 8.39  10−4  γ = 3
α = 0.03   β  = 4.76  10−4  γ = 3
C
CM
/
C
I
=
3
δ = 0.5
s  =  1v  = 8Ctot = 0.3755
s  =  1v  = 6Ctot = 0.3675
s  =  1v  = 5Ctot = 0.3009
s  =  1v  = 6Ctot = 0.3331
δ = 1
s  =  1v  = 7Ctot = 0.3724
s  =  1v  = 6Ctot = 0.3626
s  =  1v  = 5Ctot = 0.2938
s  =  1v  = 6Ctot = 0.3281
δ = 2
s  =  4v  = 1Ctot = 0.3461
s  =  4v  = 1Ctot = 0.3508
s  =  3v  = 1Ctot = 0.2661
s  =  4v  = 1Ctot = 0.2963
δ = 1/8
s  =  1v  = 8Ctot = 0.3776
s  =  1v  = 6Ctot = 0.3714
s  =  1v  = 5Ctot = 0.3064
s  =  1v  = 7Ctot = 0.3363
δ = 8
s  =  2v  = 1Ctot = 0.1772
s  =  1v  = 1Ctot = 0.1328
s  =  2v  = 1Ctot = 0.0534
s  =  2v  = 1Ctot = 0.1227
C
CM
/
C
I
=
2
δ = 0.5
s  = 1 v  = 9Ctot = 0.2886
s  =  1v  = 7Ctot = 0.2936
s  =  1v  = 6Ctot = 0.2572
s  =  1v  = 7Ctot = 0.2687
δ = 1
s  =  1v  = 9Ctot = 0.2864
s  =  1v  = 7Ctot = 0.2899
s  =  1v  = 6Ctot = 0.2522
s  =  1v  = 7Ctot = 0.2650
δ = 2
s  =  5v  = 1Ctot = 0.2666
s  =  4v  = 1Ctot = 0.2735
s  = 4 v  = 1Ctot = 0.2214
s  =  5v  = 1Ctot = 0.2336
δ = 1/8
s  =  1v  = 10Ctot = 0.2902
s  =  1v  = 8Ctot = 0.2959
s  =  1v  = 6Ctot = 0.2610
s  =  1v  = 8Ctot = 0.2711
δ = 8
s  =  2v  = 1Ctot = 0.1200
s  =  2v  = 1Ctot = 0.0914
s  =  2v  = 1Ctot = 0.0374
s  =  2v  = 1Ctot = 0.0837
C
CM
/
C
I
=
1
δ = 0.5
s  =  1v  = 13Ctot = 0.1875
s  =  1v  = 11Ctot = 0.2010
s  =  1v  = 8Ctot = 0.1985
s  =  1v  = 10Ctot = 0.1916
δ = 1
s  =  1v  = 13Ctot = 0.1864
s  =  1v  = 10Ctot = 0.1993
s  =  1v  = 8Ctot = 0.1957
s  =  1v  = 9Ctot = 0.1895
δ = 2
s  =  7v  = 1Ctot = 0.1636
s  =  6v  = 1Ctot = 0.1776
s  =  5v  = 1Ctot = 0.1632
s  =  6v  = 1Ctot = 0.1566
δ = 1/8
s  =  1v  = 14Ctot = 0.1883
s  =  1v  = 11Ctot = 0.2021
s  =  1v  = 8Ctot = 0.2006
s  =  1v  = 10Ctot = 0.1930
δ = 8
s  =  2v  = 1Ctot = 0.0628
s  =  2v  = 1Ctot = 0.0485
s  =  2v  = 1Ctot = 0.0215
s  =  2v  = 1Ctot = 0.0446Figure 6
Cost function for DMWD (α=0.03, β=4.76 10−4, and γ=3); δ=2.Once again, the choice of parameters of all used distributions in this table verified expectation 10.The summary of computation results is presented in Table5. Hence, we obtain the following remark.Remark 4.
(a) Ifδ≤1, then the state s=1 is always the optimal state after PM actions. This result is expected and can be interpreted as follows: the costs of PM actions are close to the costs of perfect repair CPR; this implies that the optimal policy will be a perfect repair strategy.
(b) Ifδ>1, then v is always equal to 1 and the optimal distance between maintenance is short. This result can be interpreted as follows: the optimum for the cost function is obtained for relatively small repairs.
(c) If the failure rate is DMWD or DWD the obtained results are quite similar.
(d) If the failure rate is that of DLFRD, then the optimal distance between PM actions is larger than that for the rest failure rates.
(e) Ifδ=8, then the optimal cost per time unit is very small.
(f) In the case of DMWD with parametersα=0.03, β=4.76 10−4, and γ=3, If CCM/CI=3 and δ=1/2, we obtain then the unique optimum for both s=1 and v=6. This implies the repair degree ζ=1/7 and the total costs are Ctot=0.3331.
(g) If we use the same parameters as in (f) and after replacingδ by 8, we obtain then the unique optimum in s=2 and v=1. It follows then that ζ=2/3 and Ctot=0.1227.


## 6. Conclusions

We considered an incomplete maintenance model; that is, the impact of PM is not minimal (ABAO) and not perfect (AGAN) but lies between these boundary cases. The PM actions reset the failure rate of the item proportional to the virtual age. We assumed that the reduction of the failure rate is done by using the Kijima’s type 2 virtual age model. Further, all CM actions are minimal, which means that the item is minimally repaired at each failure.In this research, the costs of PM actions are not constant but depend on (a) the state of the item just before the PM actions, (b) on the impact of repair, and (c) on the degree of repair. We note that the last case is very realistic and the obtained function can be well interpreted in praxis. Furthermore a periodic replacement policy is studied. We discussed the discrete failure data and the situations where the use of a continuous scale in order to measure the life length is quite artificial, so a DMWD is considered to describe the failure of the item under study.The expected total optimal maintenance cost under discrete modified Weibull distribution (DMWD) is obtained. It is shown that the optimum for the cost function is unique.Furthermore, we can consider models with imperfect repairs at failures. We expect that PM models will be applied to real systems and some results to these subjects will be studied in near future.



---

*Source: 592527-2014-02-27.xml*