The purpose of this paper is to create an interval estimation of the fuzzy system reliability for the repairable multistate series–parallel system (RMSS). Two-sided fuzzy confidence interval for the fuzzy system reliability is constructed. The performance of fuzzy confidence interval is considered based on the coverage probability and the expected length. In order to obtain the fuzzy system reliability, the fuzzy sets theory is applied to the system reliability problem when dealing with uncertainties in the RMSS. The fuzzy number with a triangular membership function is used for constructing the fuzzy failure rate and the fuzzy repair rate in the fuzzy reliability for the RMSS. The result shows that the good interval estimator for the fuzzy confidence interval is the obtained coverage probabilities the expected confidence coefficient with the narrowest expected length. The model presented herein is an effective estimation method when the sample size is n≥100. In addition, the optimal α-cut for the narrowest lower expected length and the narrowest upper expected length are considered.
1. Introduction
Most researches on reliability theory involve traditional binary reliability models where each component in a system basically consists of two functional states, perfect functionality and complete failure. However, in the system reliability of multistate components, the entire system performance will be considered from different performance levels and several failure modes. The evolution of such a system is represented by a continuous-time discrete state stochastic process. The multistate system is widely used in various industrial areas such as power generation systems, computer systems, and transportation systems (Lisnianki and Levitin [1]). Compared with a binary system model, a multistate system (MSS) model provides a more flexible tool for representing engineering systems in real life as first introduced in [2–4]. Recent research has focused on reliability evaluation and optimization of MSS [5–7]. In conventional multistate theory, it is assumed that the exact probability of each component state is given. However, with the progress of modern industrial technologies, the product development cycles have become shorter, while the lifetime of products has become longer [8]. In many highly reliable applications, there may be only a few available observations of the system’s failures. Therefore, it may be difficult to obtain sufficient data to estimate the precise values of the probabilities and performance levels of these systems. Moreover, the inaccuracy of system models, caused by human errors, is difficult to quantify using conventional reliability theory alone [9]. In light of these significant challenges, new techniques are needed to solve these fundamental problems related to reliability.
More recently, fuzzy reliability theory has been developed on the basis of fuzzy theory (Zadeh [10]). Presently, the applications of fuzzy idea in reliability theory that deal with the problem of lacking of inaccuracy or fluctuation data can be seen in many areas. In reliability analysis, many theories and methods have been constructed to facilitate the multistate system reliability assessment such as the universal generating function (Levitin [11]) and the multistate weight system (Li and Zuo [12]). Ding and Lisnianski [13] proposed the fuzzy universal generating function method to derive the fuzzy probability distribution and fuzzy system availability of the overall system when the component’s performance rate and state probabilities take fuzzy values. Liu et al. [14] investigated the dynamic fuzzy system state probabilities, fuzzy availability, and fuzzy performance rewards of a multistate system under a continuous-time Markov model. Liu and Huang [15] proposed the fuzzy multistate system that has extended the multistate system model to the cases that the transition and performance rates of multistate elements are uncertain.
This research considers the problem of interval estimation of fuzzy system reliability when the parameter of interest is fuzzy and the data are observations from fuzzy random variables. A random sample of fuzzy data set is generated in constructing a fuzzy confidence interval. In some studies on fuzzy confidence interval, Corral and Gil [16] considered the problem of constructing confidence interval using fuzzy data without considering any fuzzy random variables. Geyer and Meeden [17] introduced a fuzzy confidence interval that is the optimality of UMP and UMPU test. Viertl [18] investigated statistical inference about an unknown parameter based on fuzzy observations and developed testing hypothesis for crisp parameter based on fuzzy data. Wu [19] proposed an approach based on fuzzy random variables for constructing a fuzzy confidence interval for an unknown fuzzy parameter. Škrjanc [20] introduced a method to define a fuzzy confidence interval that combines a fuzzy identification methodology with some idea from applied statistics in finding the confidence interval defined by the lower and upper fuzzy bounds. Chachi and Taheri [21] proposed a method to construct the one-sided and two-sided fuzzy confidence intervals for an unknown fuzzy parameter based on normal fuzzy random variable. Škrjanc [22] presented a new method of confidence interval identification for Takagi-Sugeno fuzzy models in the case of the data with regionally changeable variance. The method combines a fuzzy identification methodology with some ideas from applied statistics. Some studies showed that much research proposed an approach based on fuzzy data for constructing a fuzzy confidence interval, but without considering the performance analysis of interval estimator. Then the performance of fuzzy confidence interval will be assessed based on the coverage probability and the expected length in this research.
In this research, the fuzzy system reliability of RMSS is constructed, where the fuzzy failure rate and the fuzzy repair rate for each component are the triangular fuzzy number. An approach to construct interval estimation of the fuzzy system reliability of RMSS which subsequently will be used in estimation of the fuzzy confidence interval of fuzzy system reliability is developed. Finally, the analytic expression to find the coverage probability and the expected length that it is used to interpret the efficiency of fuzzy confidence interval are presented.
2. Materials and Methods2.1. Markov Model for Multistate Element
A Markov model has been used for evaluating the expected number of failures at an arbitrary time interval in many practical cases and can be described as a Poisson process (Lisnianki and Levitin [1]).
Definition 1 (see Ibe [23]).
A stochastic process {X(t)∣t≥0} is a continuous-time Markov chain if, for all t, Δt≥0 and nonnegative integers i, j, k,
(1)P[X(t+Δt)=j∣X(Δt)=i,X(u)=k,0≤u≤Δt]=P[X(t+Δt)=j∣X(Δt)=i].
This means that in a continuous-time Markov chain the conditional probability of the future state at time t+Δt given the present state at s and all part state depends only on the present state and is independent of the past. If in addition P[X(t+Δt)=j∣X(Δt)=i] is dependent on s, then the process {X(t)∣t≥0} is said to be time homogeneous or have the time homogeneity property. Time homogeneous Markov chains have stationary transition probabilities. Let Pij(t)=P[X(t+Δt)=j∣X(Δt)=i] be the probability that a Markov chain in state i will be in state j after an additional time t. Thus, the Pij(t) is the transition probability function that satisfies the condition 0≤Pij(t)≤1. Also, ∑jPij(t)=1.
2.2. Kolmogorov Differential Equations
Using a set of differential equations to find Pij(t) (Rausand and Hoyland [24]), it is start by considering the Chapman-Kolmogorov equations
(2)Pij(t+Δt)=∑k=0rPik(Δt)Pkj(t).
The interval of (0,t+Δt) is split into two parts. First, consider a transition from state i to state k in the small interval (0,Δt) and then a transition from state k to state j in the rest of the interval. It is seen that
(3)P˙ij(t)=∑k=0k≠irαikPkj(t)-αiPij(t)=∑k=0rαikPkj(t),
where αii=-αi and the following notation for the time derivative P˙ij(t)=(d/dx)Pij(t)is introduced. The differential equation (3) is known as the Kolmogorov backward equations. They are called backward equations because we start with a transition back by the start of the interval. The Kolmogorov backward equations may also be written in matrix format as
(4)P˙(t)=A·P(t).
Likewise, split the time interval (0,t+Δt) into two parts. Consider a transition from i to k in the interval (0,t) and then a transition from k to j in the small interval (t,t+Δt). It is seen that
(5)P˙ij(t)=∑k=0k≠jrαkjPik(t)-αjPij(t)=∑k=0rαkjPik(t),
where, as before, αjj=-αj. The differential equation (5) is known as the Kolmogorov forward equations. The interchange of the limit and the sum above does not hold in all cases but is always valid when the state space is finite.
Consider the following:
(6)P˙(t)=P(t)·A.
For the Markov processes, the backward and the forward equations have the same unique solutionP(t), where ∑j=0rPij(t)=1 for all i in X.
2.3. Repairable Multistate Elements
In this section, the repairable multistate element assumes that minor failures can happen and minor repairs can be done (Xie et al. [25]). A minor failure is a failure that causes the element transition from state i to state i-1 denoted by λi,i-1. On the other hand, a minor repair is a repair that causes the element transition from state i-1 to state i denoted by μi-1,i. It is actually a birth and death process as presented in Figure 1.
The state-space diagram of repairable multistate element with minor failures and minor repairs.
The Chapman-Kolmogorov equations for the general case are as follows:
(7)dP1(t)dt=-μ1,2P1(t)+λ2,1P2(t),dP2(t)dt=μ1,2P1(t)-(λ2,1+μ2,3)P2(t)+λ3,2P3(t),⋯=⋯,dPk(t)dt=μk-1,kPk-1(t)-λk,k-1Pk(t).
Assume that the initial state is in the state k with the best performance. Therefore, by solving (7) of differential equations under the initial condition Pk(0)=1, Pk-1(0)=Pk-2(0)=⋯=P1(0)=0. The unreliability function for the multistate element will be a sum of the probabilities of the unacceptable states 1,…,i. Therefore reliability function is given by
(8)Rrepair(t)=1-∑j=1iPj(t).
2.4. Basic Concept of Fuzzy Set Theory
The basic concepts which are used for analysing the fuzzy system reliability are discussed in this section. In classical set theory, an element x in a universe U may or may not be a membership of some crisp set A. This binary membership can be represented by the following indicator function:
(9)χA(x)={1ifx∈A0ifx∉A.
Zadeh’s [10] extended the notation of binary membership to accommodate various degrees of membership on the real continuous interval [0,1] and defined the fuzzy set A~ by the membership function ξA~(x)∈[0,1], given that ξA~(x) is the degree of membership of element x in fuzzy set A~={(x,ξA~(x));x∈X}. Consider a closed interval of real numbers [a,b]={x∈R∣a≤x≤b}∀a,b∈R. The following are formulas for four basic arithmetic operations on closed intervals of real numbers (Ross [26]):
(10)[a,b]+[c,d]=[a+c,b+d],[a,b]-[c,d]=[a-d,b-c],[a,b]·[c,d]=[min(ac,ad,bc,bd),max(ac,ad,bc,bd)],[a,b][c,d]=[a,b][1/c,1/d]when0∈[c,d].
The triangular membership function of fuzzy set A~ is given by
(11)ξA~={0;x<Lx-LM-L;L≤x≤MU-xU-M;M≤x≤U0;x>U.
Let R be a universal set of real numbers and A~ a fuzzy subset of R. A~ is referred to as triangular fuzzy number (TFN) as shown in Figure 2, denoted by A~=(L,M,U). Let A~α be the α-cut level, so we have
(12)A~α=[A~αL,A~αU]=[L+α(M-L),U-α(U-M)];∀α∈[0,1].
Fuzzy input membership function (the triangular membership function).
2.5. Repairable Fuzzy Multistate Elements
Based on the conventional multistate elements, the state space diagram of a repairable multistate system takes the form presented in Figure 1 where state k is the best state and state 1 is the worst state. The minor failures between states i and i-1 are determined by fuzzy value λ~i,i-1. And the minor repairs between states i-1 and i are determined by fuzzy value μ~i-1,i. With the fuzzy transition intensities, the state probability of elements at time t must also be a fuzzy value P~i(t), and then the Chapman-Kolmogorov equation with fuzzy transition intensities can be written as
(13)dP~1(t)dt=-μ~1,2P~1(t)+λ~2,1P~2(t),dP~2(t)dt=μ~1,2P~1(t)-(λ~2,1+μ~2,3)P~2(t)+λ~3,2P~3(t),⋯=⋯dP~k(t)dt=μ~k-1,kP~k-1(t)-λ~k,k-1P~k(t),
where the initial conditions are P~k(0)=1 and P~k-1(0)=P~k-2(0)=⋯=P~1(0)=0. The unreliability function for the multistate element will be a sum of the probabilities of the unacceptable states 1,…,i(∑j=1iPj(t)). Laplace transform is adopted to transform (13) into linear equation as follows:
(14)sP~1(s)=-μ~1,2P~1(s)+λ~2,1P~2(s),sP~2(s)=μ~1,2P~1(s)-(λ~2,1+μ~2,3)P~2(s)+λ~3,2P~3(s),⋯=⋯sP~k(s)-1=μ~k-1,kP~k-1(s)-λ~k,k-1P~k(s),
where sP~i(s)-P~i(0)=L[dP~i(t)/dt]. Given that P~i(s) is a function of λ~, μ~, and s, then the inverse Laplace transform is executed to get the P~i(t) in time domain:
(15)P~i(t)=L-1[P~i(s)]=P~i*(λ~,t),
where L-1[·] is an inverse Laplace operator and P~i(t) is a function in terms of fuzzy variables λ~={λ~k,k-1,…,λ~i,j,…,λ~2,1} and μ~={μ~k-1,k,…,μ~j,i,…,μ~1,2} at any time t. The fuzzy state probabilities can be obtained in the form of P~i(t) at the α-cut level:
(16)P~iα(t)=[P~iαL(t),P~iαU(t)];t≥0,0≤α≤1,
where P~Liα(t)=minP~i*(λ,t;μλi,j(λi,j)≥α) and P~iαU(t)=maxP~i*(λ,t;μλi,j(λi,j)≥α).
Then the fuzzy reliability of repairable multistate element is given by
(17)[R~αL(t),R~αU(t)]=[1-∑i=1nP~iαU(t),1-∑i=1nP~iαL(t)].
The fuzzy state probabilities can be obtained in the form that P~iα(t)=[P~iαL(t),P~iαU(t)] at the α-cut level where α∈[0,1]. Let F~iα(t) be fuzzy unreliability functions of repairable multistate element in each element (F~iα(t)=1-R~iα(t)).
2.6. Multistate Element under Series-Parallel System
In this section, analysis of the fuzzy system reliability using the repairable multistate element with series-parallel system will be demonstrated. Figure 3 represents a system containing m subsystem connection in series where each subsystem consists of k components in parallel. Let R~ij be a fuzzy reliability function and F~ij a fuzzy unreliability function of subsystem i, which is connected in series, and component j, which is connected in parallel (i=1,2,…,m and j=1,2,…,k), respectively. Let a failure rate and a repair rate be represented by λ~ij and μ~ij at time t, respectively.
Fuzzy series-parallel system.
Then the fuzzy system reliability of repairable multistate series-parallel system (RMSS) at time t is given by
(18)R~(sp)α=∏i=1m[(1-∏j=1k(1-(R~ij)αL)),(1-∏j=1k(1-(R~ij)αU))],kkkkkkkkkkkkkkkwikkwkkkkkkkkkkkkkkkkkkaaα∈[0,1].
2.7. Fuzzy Confidence Interval Probability for Fuzzy Parameter
The interval estimation of the fuzzy system reliability of RMSS is constructed by extending the concepts of two-sided 100(1-γ)% confidence interval of the real parameters to the case where both parameter and random variables are fuzzy as shown in Figure 4.
Two-sized 100(1-γ)% confidence intervals for fuzzy parameter.
According to Wu [19], R1,…,Rn are independent and identically distributed random variables. Let L(R) and U(R) be two statistics such that L(R)≤U(R), where R=(R1,…,Rn). If the random interval [L(R),U(R)] satisfies P[L(R)≤θ≤U(R)]=1-γ, then [L(r),U(r)] is a confidence interval for θ with confidence coefficient 1-γ, where r=(r1,…,rn) and each ri is an observed value of Ri when i=1,2,…,n. It can be applied to fuzzy confidence interval. Let R~1,…,R~n be independent and identically distributed fuzzy random variables with fuzzy parameter θ~. Let r~i be the observed value of R~i for i=1,2,…,n, where each r~i is a fuzzy number for i=1,2,…,n. Therefore, r~iαL and r~iαU are the observed values of R~iαL and R~iαU for α∈[0,1]. Then R~1αL,R~2αL,…,R~nαL and R~1αU,R~2αU,…,R~nαU are identically distributed fuzzy random variables. Suppose that the distribution of R~i is unknown for i=1,2,…,n. Then the approximate fuzzy confidence interval can be constructed using the central limit theorem when the sample size n is sufficiently large (Ross [27]). Let R~1,…,R~n be independent and identically distributed fuzzy random variables. There exists a unique fuzzy number E(R~)=(E(R~αL),E(R~αU)) where R~=(R~αL,R~αU) for all α∈[0,1]; then E(R~) is called the expectation of R~. Suppose that E(R~)=μ~R, E(R~αL)=(μ~R)αL, and E(R~αU)=(μ~R)αU, respectively. Therefore R~iαL and R~iαU have finite expectations (μ~R)αL and (μ~R)αU with unknown variances for i=1,2,…,n and α∈[0,1].
Let R~αL¯=(1/n)∑i=1nR~iαLandR~αU¯=(1/n)∑i=1nR~iαU and let SαL2=(1/n)∑i=1n(R~iαL-R~αL¯)2an SαU2=(1/n)∑i=1n(R~iαU-R~αU¯)2.
It is shown that SαL2→Var(R~iαL) and SαU2→Var(R~iαU) as →∞. The central limit theorem gives that n(R~αL¯-(μ~R)αL)/SαL as well as n(R~αU¯-(μ~R)αU)/SαU has approximately N(0,1) distribution. Therefore the approximate 100(1-γ)% fuzzy confidence interval for (μ~R)αL and (μ~R)αU is given by
(19)[R~αL--zγ/2SαLn,R~αL-+zγ/2SαLn]≡[L(R~αL),U(R~αL)],(20)[R~αU--zγ/2SαUn,R~αU-+zγ/2SαUn]≡[L(R~αU),U(R~αU)],
respectively.
We assign [L(R~αL),U(R~αL)] as lower fuzzy confidence interval and [L(R~αU),U(R~αU)] as upper fuzzy confidence interval.
2.8. Coverage Probability and Expected Length of Fuzzy Confidence Interval
For an interval estimator [L(R),U(R)] of a true parameter θ, the coverage probability of [L(R),U(R)] is the probability that the random interval [L(R),U(R)] covers θ. It is denoted by P(θ∈[L(R),U(R)]). In this research, the performance of fuzzy confidence interval is assessed based on the coverage probability and the expected length. The analytic expressions for the coverage probability of fuzzy confidence interval of true parameter μ~Rα=[(μ~R)αL,(μ~R)αU] are derived.
Let P((μ~R)αL∈[L(R~αL),U(R~αL)]) and P((μ~R)αU∈[L(R~αU),U(R~αU)]) be the coverage probability of fuzzy confidence interval [L(R~αL),U(R~αL)] and [L(R~αU),U(R~αU)], respectively. Then the lower coverage probability of fuzzy confidence interval of (μ~R)αL is given by
(21)P((μ~R)αL∈[L(R~αL),U(R~αL)])=P[L(R~αL)<(μ~R)αL<U(R~αL)]=P[A<Z<B]=E[I{A<Z<B}(τ)],
where
(22)I{A<Z<B}(τ)={1,ifτ∈{A<Z<B}0,otherwise,
where A=R~αL¯-U(R~αL)/(SαL/n), B=R~αL¯-L(R~αL)/(SαL/n), and Z is the standard normal distribution.
Likewise, the upper coverage probability of fuzzy confidence interval of (μ~R)αU is given by
(23)P((μ~R)αU∈[L(R~αU),U(R~αU)])=P[L(R~αU)<(μ~R)αU<U(R~αU)]=P[C<z<D]=E[I{C<Z<D}(τ)],
where
(24)I{C<Z<D}(τ)={1,ifτ∈{C<Z<D}0,otherwise,
where C=R~αU¯-U(R~αU)/(SαU/n), D=R~αU¯-L(R~αU)/(SαU/n), and Z is the standard normal distribution.
2.9. General Procedure for Investigating the Coverage Probabilities and Expected Length for Fuzzy Confidence Interval of Fuzzy System Reliability
In this section, the coverage probability and expected length for fuzzy system reliability of RMSS is investigated. Suppose that fuzzy parameter (population mean) is denoted by μ~Rα=[(μ~R)αL,(μ~R)αU]. A method to calculate the fuzzy coverage probability and expected length for μ~Rα is demonstrated as follows.
Step 1.
Generate fuzzy random variables from normal distribution with R~1αL,R~2αL,…,R~nαL~N((μ~R)αL,σ2) and R~1αU,R~2αU,…,R~nαU~N((μ~R)αU,σ2). The sample sizes of this example are n=30,50,100,150,200 and σ2=1 for α∈[0,1].
Step 2.
Compute fuzzy sample mean [R~αL¯,R~αU¯] and fuzzy sample variance [SαL2,SαU2] for α∈[0,1].
Step 3.
Compute fuzzy confidence intervals of μ~R=[(μ~R)αL,(μ~R)αU] for α∈[0,1], the lower fuzzy confidence interval [L(R~αL),U(R~αL)], and the upper fuzzy confidence interval [L(R~αU),U(R~αU)].
Step 4.
Compute the fuzzy coverage probability (CPα) and the expected length (ELα) for α∈[0,1],
(25)CPα=(Numberoftimesthatthefuzzyconfidencewwiintervalcoversthetrueparameterineachα-cut)×m-1,ELα=∑i=1mUi-Lim,
where Ui is the lower fuzzy confidence interval and Li is the upper fuzzy confidence interval at repeated time i, i=1,2,…,m.
Step 5.
Repeat Steps 2–4 for a given condition.
After receiving the coverage probabilities and the expected length of the fuzzy confidence interval, the next step is to consider the sample size (n) that gives the coverage probability higher than the expected confidence coefficient, where the number of repeatition is m=10,000. The Z-test statistic hypothesis testing is used in confirming the level of the coverage probability as follows:
the coverage probability is not less than the expected confidence coefficient (CP≥CP0),
the coverage probability is lower than the expected confidence coefficient (CP<CP0),
when the significance level is γ=0.95 and the test statistic
(26)Z=C^P-CP0CP0(1-CP0)/m.
Since the criterion in the test is Z<-Zγ, the hypothesis H0 cannot be rejected if
(27)C^P≥CP0-ZγCP0(1-CP0)m,
where CP is the coverage probability, C^P is the coverage probability estimated from this study, CP0 is the expected confidence coefficient, and m is the repeated time.
3. Numerical Example and Results
In this section, a RMSS which consists of 3 subsystems in series and 2 components in each subsystem in parallel is considered. In each element, it has repairable multistate with a fuzzy failure rate and a fuzzy repair rate as shown in Figure 5. Since a failure rate and a repair rate cannot be recorded precisely due to human errors, machine errors, or some unexpected situations, triangular fuzzy numbers are used to describe the fuzzy failure rate and fuzzy repair rate. The parameters of these functions are shown in Table 1.
Triangular fuzzy number of fuzzy failure rates and fuzzy repair rates (per year).
Subsystem 1
Fuzzy failure rate
Fuzzy repair rate
λ~2,111=(0.20,0.25,0.35)
μ~1,211=(0.015,0.020,0.025)
λ~2,112=(0.30,0.35,0.40)
μ~1,212=(0.010,0.015,0.020)
Subsystem 2
Fuzzy failure rate
Fuzzy repair rate
λ~2,121=(0.22,0.25,0.30)
μ~1,221=(0.025,0.030,0.035)
λ~2,122=(0.20,0.25,0.35)
μ~1,222=(0.022,0.025,0.030)
Subsystem 3
Fuzzy failure rate
Fuzzy repair rate
λ~2,113=(0.15,0.20,0.22)
μ~1,213=(0.012,0.015,0.022)
λ~2,123=(0.28,0.35,0.39)
μ~1,223=(0.018,0.020,0.025)
Multistate system and state diagrams for RMSS.
For Markov model, the Chapman-Kolmogorov equations of each element can be written as
(28)dP~1ij(t)dt=-μ~1,2ijP~1ij(t)+λ~2,1ijP~2ij(t)dP~2ij(t)dt=μ~1,2ijP~1ij(t)-λ~2,1ijP~2ij(t).fori=1,2,3,j=1,2,
where P~1ij(t) is probabilities of good function and P~2ij(t) is probabilities of fail function with fuzzy failure rate (λ~ij) and fuzzy repair rate (μ~ij) for each element in multistate model. Using the inverse Laplace transform, the fuzzy state probabilities are obtained as functions of time t in the form of
(29)P~1ij(t)=λ~ij2,1λ~ij2,1+μ~ij1,2(e-(λ~ij2,1+μ~ij1,2)·t-1)P~2ij(t)=μ~ij1,2λ~ij2,1+μ~ij1,2+λ~ij2,1λ~ij2,1+μ~ij1,2e-(λ~ij2,1+μ~ij1,2)·twwwwwwwwfori=1,2,3,j=1,2.
From Figure 5, fuzzy reliability functions and fuzzy unreliability functions of repairable for each element are given by
(30)R~ij(t)=P~2ij(t),F~ij(t)=P~1ij(t)fori=1,2,3,j=1,2,
where F~ij(t) is fuzzy unreliability function of repairable multistate system at subsystem i connected in series and component j connected in parallel, respectively. Suppose that λ~2,1ij and μ~2,1ijare triangular membership function in each α∈[0,1], as shown in Table 1. It can be written in the form of α-cut as follows:
(31)[(λ~2,1ij)αL,(λ~2,1ij)αU]=[aij+α(bij-aij),cij-α(cij-bij)],[(μ~2,1ij)αL,(μ~2,1ij)αU]=[eij+α(eij-dij),fij-α(fij-eij)].
From the example system in Figure 6, RMSS with fuzzy failure rate (λ~2,1ij)α and a fuzzy repair rate (μ~2,1ij)αare considered. The fuzzy system reliability in each α∈[0,1] is given by
(32)R~sα=[R~sαL(t),R~sαU(t)]=∏i=13[[1,1]α-[∏j=12(1-(R~ij)αU),∏j=12(1-R~ij)αL]].
Substituting fuzzy reliability function (R~ij(t)) in each element into (32), then the fuzzy system reliability of RMSS is shown in Table 2.
The fuzzy parameter of fuzzy system reliability for RMSS in each α∈[0,1].
α-cut
R~sαL(t)=(μ~R)αL
R~sαU(t)=(μ~R)αU
0
0.5357
0.9587
0.1
0.5825
0.9530
0.2
0.6257
0.9465
0.3
0.6654
0.9392
0.4
0.7016
0.9310
0.5
0.7346
0.9219
0.6
0.7644
0.9116
0.7
0.7913
0.9001
0.8
0.8155
0.8872
0.9
0.8371
0.8727
1
0.8565
0.8565
Fuzzy 95% confidence intervals for fuzzy parameter μ~Rα=[(μ~R)αL,(μ~R)αU].
Suppose that fuzzy system reliability for RMSS in Table 2 is the fuzzy parameter (population mean) denoted by μ~Rα=[(μ~R)αL,(μ~R)αU]. The confidence coefficient with γ=0.95 is defined and fuzzy random variables from normal distribution are generated by using MATLAB [28]. Then the fuzzy coverage probability and the expected length for μ~Rα are shown in Figure 6.
After receiving the coverage probabilities and the expected length of fuzzy confidence interval, the next step is to compare the values between fuzzy coverage probabilities and the expected confidence coefficient by the test of hypothesis. In this example, let m=10,000 and the criterion used in comparing the coverage probability at significance level 95% will be
(33)C^P≥0.95-1.6450.95(1-0.95)10,000=0.9464.
Considering a fuzzy confidence interval at α-cut which gives the coverage probability higher than 0.9464 at significant level 95% will be the coverage probability that is covered in the expected confidence coefficient. Then only the coverage probability that is covered in the expected confidence coefficient will be used in the most appropriate expected length estimation. In addition, the fuzzy confidence interval of the α-cut which gives the narrowest expected length is considered as the most appropriate expected length. Both lower and upper fuzzy coverage probabilities and expected lengths results are shown in Tables 3 and 4, respectively.
Lower fuzzy coverage probabilities and lower expected lengths for 95% confidence interval where m=10,000 and C^P≥0.9464.
α-cut
n=30
n=50
n=100
n=150
n=200
Lower CP
Lower EL
Lower CP
Lower EL
Lower CP
Lower EL
Lower CP
Lower EL
Lower CP
Lower EL
0
0.9402
0.7091
0.9456
0.5524
0.9485
0.3908
0.9532
0.3196
0.9500
0.2769
0.1
0.9377
0.7102
0.9426
0.5518
0.9472
0.3908
0.9494
0.3199
0.9464
0.2769
0.2
0.9390
0.7089*
0.9440
0.5516
0.9467
0.3909
0.9466
0.3196
0.9487
0.2768
0.3
0.9425
0.7090
0.9443
0.5518
0.9489
0.3908
0.9465
0.3196
0.9467
0.2768
0.4
0.9443
0.7100
0.9409
0.5521
0.9518
0.3914
0.9467
0.3201
0.9491
0.2772
0.5
0.9414
0.7093
0.9441
0.5513*
0.9493
0.3908
0.9465
0.3195
0.9468
0.2767
0.6
0.9411
0.7088
0.9409
0.5516
0.9466
0.3913
0.9504
0.3196
0.9502
0.2769
0.7
0.9387
0.7101
0.9465
0.5520
0.9480
0.3911
0.9476
0.3195
0.9467
0.2769
0.8
0.9460
0.7109
0.9396
0.5515
0.9469
0.3911
0.9482
0.3194
0.9477
0.2770
0.9
0.9417
0.7090
0.9477
0.5513*
0.9466
0.3907*
0.9481
0.3193*
0.9475
0.2766*
1
0.9408
0.7102
0.9442
0.5515
0.9510
0.3908
0.9511
0.3198
0.9555
0.2769
Coverage probabilities that are higher than the expected confidence coefficient are shown in bold numbers.
*The narrowest expected length.
Upper fuzzy coverage probabilities and upper expected lengths for 95% confidence interval where m=10,000 and C^P≥0.9464.
α-cut
n=30
n=50
n=100
n=150
n=200
Upper CP
Upper EL
Upper CP
Upper EL
Upper CP
Upper EL
Upper CP
Upper EL
Upper CP
Upper EL
0
0.9392
0.7092
0.9460
0.5514
0.9469
0.3908
0.9465
0.3197
0.9501
0.2769
0.1
0.9435
0.7112
0.9466
0.5511
0.9466
0.3910
0.9494
0.3194
0.9482
0.2768
0.2
0.9404
0.7117
0.9482
0.5516
0.9475
0.3909
0.9483
0.3196
0.9465
0.2768
0.3
0.9360
0.7098
0.9425
0.5516
0.9470
0.3914
0.9477
0.3197
0.9470
0.2771
0.4
0.9384
0.7092
0.9441
0.5506*
0.9484
0.3911
0.9475
0.3193
0.9470
0.2767
0.5
0.9411
0.7101
0.9444
0.5514
0.9465
0.3912
0.9474
0.3197
0.9477
0.2769
0.6
0.9420
0.7101
0.9410
0.5517
0.9509
0.3908
0.9465
0.3196
0.9467
0.2768
0.7
0.9426
0.7084
0.9421
0.5511
0.9477
0.3906*
0.9466
0.3191*
0.9474
0.2765*
0.8
0.9400
0.7096
0.9425
0.5514
0.9475
0.3910
0.9508
0.3195
0.9466
0.2768
0.9
0.9387
0.7075*
0.9436
0.5522
0.9469
0.3907
0.9509
0.3198
0.9475
0.2768
1
0.9453
0.7103
0.9412
0.5509
0.9466
0.3907
0.9489
0.3195
0.9525
0.2767
Coverage probabilities that are higher than the expected confidence coefficient are shown in bold numbers.
*The narrowest expected length.
4. Numerical Results
From the numerical example, the fuzzy reliability for the RMSS is calculated. Estimation of fuzzy confidence interval for fuzzy system reliability model (μ~R) in each α∈[0,1] revealed that the fuzzy confidence interval can be divided into 2 parts which are lower fuzzy confidence interval and upper fuzzy confidence interval as shown in Figure 6.
From Figure 6, the lower bound and upper bound of the fuzzy parameter are estimated at a 95% significance level. Next, the performance of the estimated parameter of the fuzzy confidence interval of fuzzy system reliability model is considered. The coverage probability and expected length is used in calculation with the above method with repeated time at m=10,000.
In Figure 7(a), it is seen that at the sample size n=30,50(n<100) there are some α-cut that give lower coverage probability covering in the expected confidence coefficient when C^P≥0.9464. However, if considering, at sample size n=100,150,200(n≥100), the lower coverage probability of every α∈[0,1] is covering in the expected confidence coefficient, likewise, Figure 7(b) shows that at the sample size n=30,50(n<100) there are only few α-cut numbers that contain upper coverage probability in the expected confidence coefficient at C^P≥0.9464, but, at the sample size =100,150,200(n≥100), the upper coverage probability of everyα∈[0,1] is in the expected confidence coefficient.
Boxplots of the lower and upper fuzzy coverage probability with sample size n=30,50,100,150, and 200 at significance level 95%.
Then only the coverage probability that is covered in the expected confidence coefficient will be used in estimation of the most appropriate expected length. The result shows that the fuzzy confidence interval at α-cut=0.9 gives the narrowest lower expected length. It is considered as the most appropriate lower expected length. It is also seen that the fuzzy confidence interval at α-cut=0.7 gives the narrowest upper expected length. It is considered as the most appropriate upper expected length as shown in Tables 3 and 4.
It seems that at the larger sample size (n) more α-cut numbers are obtained for the estimated parameter of the fuzzy confidence interval of fuzzy system reliability model. The good interval estimator for the fuzzy confidence interval is the obtained coverage probabilities that can cover the expected confidence coefficient C^P≥0.9464 with the narrowest expected length.
5. Conclusions and Discussions
This paper presents an innovative modeling approach when dealing with uncertainties in the RMSS. The Markov process for the RMSS with a fuzzy failure rate and a fuzzy repair rate is considered and the fuzzy system reliability is constructed. Interval estimation of the fuzzy system reliability model is constructed by extending the concepts of two-sided 100(1-γ)% confidence interval of the true parameters to the case where both parameter and random variables are fuzzy based on the central limit theorem. Recently, much research proposed an approach based on fuzzy data for constructing a fuzzy confidence interval, but without considering the performance analysis of interval estimator. In this research, the performance of fuzzy confidence interval is assessed based on the coverage probability and the expected length.
From the study, it is seen that estimation of the fuzzy confidence interval of fuzzy system reliability for RMSS will be effective when the sample size is n≥100. This results in lower coverage probability and upper coverage probability which covered in the expected confidence coefficient at α∈[0,1] and the narrowest lower expected length when α-cut=0.9 and the narrowest upper expected length when α-cut=0.7. Accordingly, we conclude that the model presented herein is an effective estimation method at α-cut=0.9 and α-cut=0.7. This study is suitable for the system reliability of multistate system where the accurate data are fuzzy values. In further work, fuzzy confidence intervals of system reliability for more complex systems are created. In particular, performance of interval estimator is also being considered based on the coverage probability and the expected length.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work is supported by the Science Achievement Scholarship of Thailand (SAST) and the Faculty of Science, KMUTT.
LisnianskiA.LevitinG.2003SingaporeWorld scientific PublishingBarlowR. E.WuA. S.Coherent systems with multi-state components1978342752812-s2.0-0018031576MurchlandJ.Fundamental concepts and relations for reliability analysis of multi-state systems and fault tree analysis1975Philadelphia, Pa, USASIAM581618El-NeveihiE.ProchanF.SetharamanJ.Multi-state coherent systems197815675688YehW.-C.The k-out-of-n acyclic multistate-node networks reliability evaluation using the universal generating function method20069178008082-s2.0-3364524381010.1016/j.ress.2005.08.002HuangJ.ZuoM. J.FangZ.Multi-state consecutive-k-out-of-n systems20033565275342-s2.0-003813665710.1080/07408170304418LevitinG.LisnianskiA.Ben-HaimH.ElmakisD.Redundancy optimization for series-parallel multi-state systems19984721651722-s2.0-003208365810.1109/24.722283HuangH.-Z.ZuoM. J.SunZ.-Q.Bayesian reliability analysis for fuzzy lifetime data200615712167416862-s2.0-3364626407110.1016/j.polymdegradstab.2005.12.004HuangH.-Z.TongX.ZuoM. J.Posbist fault tree analysis of coherent systems20048421411482-s2.0-214280877810.1016/j.ress.2003.11.002ZadehL. A.Fuzzy sets1965833383532-s2.0-34248666540LevitinG.2005London, UKSpringerLiW.ZuoM. J.Reliability evaluation of multi-state weighted k-out-of-n systems20089311611682-s2.0-3454857466610.1016/j.ress.2006.11.009DingY.LisnianskiA.Fuzzy universal generating functions for multi-state system reliability assessment200815933073242-s2.0-3604904612910.1016/j.fss.2007.06.004LiuY.HuangH.-Z.LevitinG.Reliability and performance assessment for fuzzy multistate element200822267568610.1243/1748006XJRR180LiuY.HuangH.-Z.Reliability assessment for fuzzy multi-state systems20104143653792-s2.0-7795109870110.1080/00207720903042939CorralN.GilM. A.A note on interval estimation with fuzzy data19882822092152-s2.0-38249028344GeyerC. J.MeedenG. D.Fuzzy and randomized confidence intervals and P-values20052043583662-s2.0-3114447254310.1214/088342305000000340ViertlR.Univariate statistical analysis with fuzzy data20065111331472-s2.0-3375030004210.1016/j.csda.2006.04.002WuH.-C.Statistical confidence intervals for fuzzy data2009362267026762-s2.0-5634910941310.1016/j.eswa.2008.01.022ŠkrjancI.Confidence interval of fuzzy models: an example using a waste-water treatment plant20099621821872-s2.0-6424908320610.1016/j.chemolab.2009.01.009ChachiJ.TaheriS. M.Fuzzy confidence intervals for mean of Gaussian fuzzy random variables2011385524052442-s2.0-7915147432610.1016/j.eswa.2010.10.034ŠkrjancI.Fuzzy confidence interval for pH titration curve2011358408340902-s2.0-7995511656910.1016/j.apm.2011.02.033IbeO. C.2009San Diego, Calif, USAElsevier Academic PressRausandM.HoylandA.2004New York, NY, USAJohn Wiley & SonsXieM.DaiY. S.PohK.-L.2004New York, NY, USAKluwer Academic/Plenum PublishersRossT. J.20103rdLondon, UKJohn Wiley & SonsRossG. G.1972New York, NY, USAAddison-WesleyThe Math WorksMATLAB, 7.6.0(R2009a)License Number 350306, February 2009