1,090
Views
3
CrossRef citations to date
0
Altmetric
Articles

Time-dependent reliability analysis for repairable consecutive-k-out-of-n:F system

ORCID Icon &
Pages 139-147 | Received 27 Jan 2020, Accepted 08 Aug 2021, Published online: 13 Sep 2021

Abstract

In a repairable consecutive C(k,n:F) system, after the system operates for a certain time, some components may fail, some failed components may be repaired and the state of the system may change. The models developed in the existing literature usually assume that the state of the system varies over time depending on the values of n and k and the state of the system is known. Since the system reliability will vary over time, it is of great interest to analyse the time-dependent system reliability. In this paper, we develop a novel and simple method that utilizes the eigenvalues of the transition rate matrix of the system for the computation of time-dependent system reliability when the system state is known. In addition, the transition performance probabilities of the system from a known state to the possible states are also analysed. Computational results are presented to illustrate the applicability and accuracy of the proposed method.

1. Introduction

In many engineering applications, consecutive-k-out-of-n systems are encountered and the evaluation of the system reliability is of interest. For example, consecutive-k-out-of-n systems appear in telecommunications, microwave relay stations, oil pipeline systems, vacuum systems in accelerators, computer ring networks and spacecraft relay stations. These systems can be classified with respect to the logical or physical connections among components as either linear or circular and the functioning principle as either failed (F) or good (G). A consecutive-k-out-of-n: F system, or a C(k,n:F) system for short, consists of an ordered sequence of n components such that the system fails if and only if at least k consecutive components fail. Another special type of system related to the C(k,n:F) system is the consecutive-k-out-of-n: G system, or a C(k,n:G) system for short. A C(k,n:G) system consists of an ordered sequence of n components such that the system works if and only if at least k consecutive components work. Both C(k,n:F) system and C(k,n:G) system can be either linear or circular. The first report for these systems was presented by Kontoleon (Citation1980). Then, much research work has been done on consecutive-k-out-of-n systems in the literature since Kontoleon (Citation1980). For further reference, see, for example, Chiang and Niu (Citation1981), Bollinger and Salvia (Citation1982), Derman et al. (Citation1982), Zuo and Kuo (Citation1990), Chang et al. (Citation2000), Kuo and Zuo (Citation2003), Eryilmaz (Citation2010) and Gökdere and Gurcan (Citation2016). However, in the aforementioned work, it was assumed that the components in the system are non-repairable. If we consider that components are repairable, then such a system is called a repairable system.

Repairable C(k,n:F) and C(k,n:G) systems have also been studied extensively in the literature (see, for example, Cheng & Zhang, Citation2001; Eryilmaz, Citation2014; Guan & Wu, Citation2006; Hongda et al., Citation2019; Lam & Ng, Citation2001; Lam & Zhang, Citation1999Citation2000; Wang et al., Citation2021; Yam et al., Citation2003; Yuan & Cui, Citation2013; Zhang & Lam, Citation1998; Zhang & Wang, Citation1996; Zhang et al., Citation1998Citation2000). In most of these research papers, it is assumed that the working time and the repair time of the components are both exponentially distributed and the repair completely restores all the properties of failed components. Moreover, when the system is in a failed state, each failed component in the system is classified as either a key component or an ordinary component. A failed component is called a key component if repairing this component will return the system to a working state (Zhang et al., Citation2000). The key components have a higher priority in repair than ordinary components. Then, by using the definition of generalized transition probability, for both linear and circular systems, the state transition rate matrix of the system can be derived. Usually, the Laplace transform method or Runge-Kutta method is used to obtain the reliability of a consecutive-k-out-of-n system.

In most of the models for repairable C(k,n:F) and C(k,n:G) systems in the existing literatures, the system reliability is calculated by assuming that all components are working in the system within a certain period of time. However, in a repairable system, the state of the system may change over time after the system starts operating. In this paper, we study the repairable linear and circular C(k,n:F) systems by assuming that the state of the system is known. We aim to provide a simple method to determine the time-dependent system reliability and analyse the transition performance probabilities of the system when the system state is known. We have first modified the transition rate matrix considered in Cheng and Zhang (Citation2001) and Yam et al. (Citation2003) and then proposed a novel method to compute the time-dependent system reliability using the eigenvalues of the transition rate matrix. The proposed method can be served as a simple alternative method to the Laplace transform and the Runge-Kutta method for the determination of time-dependent reliability and transition performance probabilities of the system. Also, it can be applied to both repairable and irreparable systems. In Section 4.1 and 4.2, to show the accuracy of our method we consider the case that the system is irreparable.

This paper is organized as follows. In Section 2, the preliminary assumptions, the state transition probabilities and the transition rate matrix are presented for the repairable C(k,n:F) system with n linearly and circularly arranged components. In Section 3, we present the proposed method for evaluation of the system reliability and transition performance probabilities when the state of the system is known. In Section 4, some numerical examples are provided to illustrate the usefulness of the proposed method for different linear and circular repairable C(k,n:F) systems. In Section 5, some concluding remarks are provided.

2. Repairable C(k,n:F) systems and model assumptions

Here is the list of the model assumptions used throughout this paper:

Assumption 2.1

The system is a linear or circular repairable C(k,n:F) system with identical components.

Assumption 2.2

At time t = 0, all components are new and working.

Assumption 2.3

There is a single repairman in the system and the repair completely restores all properties of failed components.

Assumption 2.4

The working time and the repair time of the components are both exponentially distributed with the parameters λ>0 and μ>0, respectively.

Assumption 2.5

While the system is in the failed state, the components that have not failed will not fail.

Assumption 2.6

The probability that two or more than two components fail or complete their repairs simultaneously in a very small time interval is negligible.

Based on these model assumptions, the state of the system at time t, denoted by N(t), varies depending on the values of n and k. For linear and circular repairable C(k,n:F) systems, the possible states can be expressed as follows: N(t)={0,if at time t,all components work thesystem works,1,if at time t,one component fails thesystem works,k,if at time t,k components fail thesystem works,d,if at time t,d components fail thesystem works,1,if at time t,l components fail thesystem fails, where l=k,k+1,,d+1 and d takes different values depending on whether the system is linear or circular. If the system is linear, d=n[n/k] and [x] is the largest integer less than or equal to value x. If the system is circular, d=nn/k and x is the smallest integer greater than or equal to value x.

In our proposed method, all states (k,k+1,,d+1) that caused the system to fail were considered as a single state. Then, {N(t),t0} is a continuous-time homogeneous Markov process with state space Ω={0,1,,d,1}. The set of working state is W={0,1,,d} and the set of failed states is F={1}.

According to N(t), the generalized transition probability from state i to state j in time Δt is defined as (1) pij(Δt)=P{N(t+Δt)=jN(t)=i}.(1) In order to obtain the transition rate matrix Q, we need to calculate the state transition probabilities of the system. These state transition probabilities are presented as follows (for the detailed derivations, one can refer to Cheng & Zhang, Citation2001; Yam et al., Citation2003):

For i = 0 and jΩ in (1), we have (2) p0j(Δt)={1nλΔt+o(Δt),j=0,nλΔt+o(Δt),j=1,o(Δt),for all other j values.(2) For 1id and jΩ in (1), we have (3) pij(Δt)={μΔt+o(Δt),j=(i1),(i+1)M(i+1)MiλΔt+o(Δt),j=(i+1),1((ni)λ+μ)Δt+o(Δt),j=i,o(Δt),for all otherj values.(3) For an explain (Equation3) with j=(i+1), refer to Section 4 of Guan and Wu (Citation2006). In Equation (Equation3), Mi is the number of different possible cases which the system under consideration in state i. We need to find a formula of Mi for iW and i0.

In linear C(k,n:F) systems Mi=MiL=r=0m(1)r(ni+1r)(nkrni), where m=min(ni+1,[i/k]).

In circular C(k,n:F) systems Mi=MiC=MiLl=kmin(i,2(k1))(2k1l)×r=0ml(1)r(ni1r)(nlkr2ilkr), where ml=min(ni+1,[(il)/k]).

For 1id and j = 1 in (1), we have (4) pi1(Δt)=((ni)(i+1)M(i+1)Mi)λΔt+o(Δt).(4) For i = 1 and jΩ in (Equation1), we have (5) p1j(Δt)=o(Δt).(5) Using the state transition probabilities in the system in Equations (Equation2)–(Equation5), we can obtain the transition rate matrix Q as Q=(qij), where qij=limΔt0pij(Δt)Δt, and qj=limΔt01pjj(Δt)Δt, for i,jΩ and ij. Also, qjj=qj, jΩ. In the following section, we present a novel and simple method to compute the time-dependent system reliability by using the eigenvalues of the matrix Q.

3. Proposed method

Let Υ(t)=(Ψij(t);i,jΩ and t0) be a transition matrix function defining a continuous-time homogeneous Markov process, where Ψij(t) represents the time-dependent performance probability of the system between known state i and possible state j at time t. If Q was a finite matrix the solution of ddtΥ(t)=Υ(t)Q could be written down at once in the form Υ(t)=eQt=k=0(Qt)kk!,Q0=I, where I=(Iij) is the unit matrix for i,jΩ.

The abovementioned matrix Q is a ξ×ξ real square matrix, where ξ is equal to the number of elements of the state space Ω. Then, there exists a singular value decomposition of Q of the form Q=UΣVT, where Σ=diag(σ1,σ2,,σξ1,σξ) has diagonal elements arranged in ascending order of magnitude and U=(u1,u2,,uξ) and V=(v1,v2,,vξ) are orthogonal. Now, consider the linear transformation of ξ-dimensional vectors defined by matrix Q. In this case, Ψij(t) is found to be (6) Ψij(t)=A0ij+A1ijeσ1t+A2ijeσ2t++Aξ1ijeσξ1t,(6) where A0ij,A1ij,,Aξ1ij are constant coefficients for i,jΩ and σ1,σ2,,σξ1 are the non-zero eigenvalues of Q. To obtain the above constant coefficients, the following system of differential equations is solved: (7) ddtΥ(k)(t=0)=Qk,(7) for k=0,1,2,,ξ1. We can express (Equation7) as A0ij+A1ij+A2ij++Aξ1ij=(Iij),(σ1)A1ij+(σ2)A2ij++(σξ1)Aξ1ij=(qij),(σ1)2A1ij+(σ2)2A2ij++(σξ1)2Aξ1ij=(qij)2,(σ1)ξ1A1ij+(σ2)ξ1A2ij++(σξ1)ξ1Aξ1ij=(qij)ξ1. According to the Equation (Equation6) for j = 1, we can obtain the time-dependent system reliability at time t when the state of the system is known as (8) Ri(t)=1Ψi1(t),(8) where iΩ.

In the proposed method, we need to obtain the eigenvalues of the matrix Q. When μ=0, the eigenvalues of the Q can be easily calculated based on the parameter λ. However, obtaining the eigenvalues of the matrix Q can be difficult and complex in general based on the parameters λ and μ. In such a case, we overcome this difficulty by providing the numerical values of the parameters to simplify the computation. It is noteworthy that the existing methods in the literature deal with the assumption that all components in the system are working, whereas such an assumption is not needed in our proposed method. Hence, the proposed method can be applied to a wider range of situations in system engineering.

4. Numerical illustrations

In the following, numerical examples in the existing literature are used to illustrate the proposed method developed in Section 3.

4.1. A linear repairable C(2,5:F) system

In this section, we analyse a linear repairable C(2,5:F) system. First, we focus on obtaining the time-dependent transition performance probabilities and the system reliability under the condition that there is a single failed component in the system at time t. Then, to show the accuracy of the proposed method, the situation where μ=0 is analysed. Based on the results presented in Section 2, the state space Ω={0,1,2,3,1}, the set of working state W={0,1,2,3} and the set of failed states F={1} of the system. The Q-matrix can be obtained as Q=(5λ5λ0μ(μ+4λ)125λ0μ(μ+3λ)00μ00000085λ12λ52λ(μ+2λ)2λ00).Now consider the Q-matrix for λ=0.5 and μ=1.5. For i = −1 and jΩ in (6) we have (9) Ψ1j(t)=A01j+A11jeσ1t+A21jeσ2t+A31jeσ3t+A41jeσ4t,(9) where A01j, A11j, A21j, A31j and A41j are constant coefficients and σ1=5.5045, σ2=3.2078, σ3=2.1437 and σ4=0.6439 are the non-zero eigenvalues of Q. To determine the above constant coefficients, Equation (Equation7) has to be solved for jΩ and k = 0, 1, 2, 3, 4. Finally, we have the following results for Equation (Equation9).

For j = 0, (10) Ψ10(t)=0.2859e5.5045t0.0545e3.2078t+0.0295e2.1437t+0.3209e0.6439t.(10) For j = −1, (11) Ψ11(t)=0.5726e5.5045t+0.0257e3.2078t+0.0047e2.1437t+0.397e0.6439t.(11) For j = −2, (12) Ψ12(t)=0.2887e5.5045t+0.096e3.2078t+0.0284e2.1437t+0.2212e0.6439t.(12) For j = −3, (13) Ψ13(t)=0.024e5.5045t0.0339e3.2078t0.0199e2.1437t+0.0298e0.6439t.(13) For j = 1, (14) Ψ11(t)=10.022e5.5045t0.0332e3.2078t+0.0241e2.1437t0.9689e0.6439t.(14) Using Equations (Equation10)–(Equation14), for t=0.5,1.0,,5, we can obtain the time-dependent transition performance probabilities for the linear repairable C(2,5:F) system (see, Table ). Moreover, using Equation (Equation14) for i = −1 in Equation (Equation8), we can obtain the time-dependent reliability of the system when a single component is in the failed state at time t.

Table 1. Time-dependent transition performance probabilities for the linear repairable C(2,5:F) system when one component fails.

In order to demonstrate the accuracy of our method, we consider the case that μ=0 and the system is linear irreparable C(2,5:F) system. The Q-matrix in this situation becomes Q˘=(5λ5λ00004λ125λ085λ003λ12λ52λ0002λ2λ00000). Assume that at time t all the five components are working. To obtain the time-dependent system reliability, let i = 0 and j = 1 in Equation (Equation6), then (15) Ψ01(t)=A001+A101eσ1t+A201eσ2t+A301eσ3t+A401eσ4t,(15) where A001, A101, A201, A301 and A401 are constant coefficients and σ1=5λ, σ2=4λ, σ3=3λ and σ4=2λ are the non-zero eigenvalues of Q˘. To determine the constant coefficients, we have to solve the following system of differential equations: A001+A101+A201+A301+A401=0,(5λ)A101+(4λ)A201+(3λ)A301+(2λ)A401=0,(5λ)2A101+(4λ)2A201+(3λ)2A301+(2λ)2A401=8λ2,(5λ)3A101+(4λ)3A201+(3λ)3A301+(2λ)3A401=42λ3,(5λ)4A101+(4λ)4A201+(3λ)4A301+(2λ)4A401=140λ4. The solutions of the above equations are A001=1, A101=1, A201=4, A301=3, A401=1. Finally, Equation (Equation15) can be obtained as (16) Ψ01(t)=1e5λt+4e4λt3e3λte2λt.(16) Using Equation (Equation16) for i = 0 in Equation (Equation8), we can obtain the time-dependent reliability of the system as (17) R0(t)=1Ψ01(t)=e2λt+3e3λt4e4λt+e5λt.(17) Note that Equation (Equation17) is equivalent to Equation (Equation12) presented in Zhang and Wang (Citation1996).

4.2. A circular repairable C(2,6:F) system

In this section the circular repairable C(2,6:F) system is considered. The state space is Ω={0,1,2,3,1}, the set of working state is W={0,1,2,3} and the set of failed states is F={1}. Based on the results in Section 2, the Q-matrix is Q=(6λ6λ0μ(μ+5λ)3λ0μ(μ+4λ)00μ0000002λ23λ103λ(μ+3λ)3λ00).Assume that, at time t, two components fail. Let, λ=0.5 and μ=1.5, for i = −2 and jΩ in Equation (Equation6), and then we have (18) Ψ2j(t)=A02j+A12jeσ1t+A22jeσ2t+A32jeσ3t+A42jeσ4t,(18) where A02j, A12j, A22j, A32j and A42j are constant coefficients and σ1=6.2507, σ2=3.7738, σ3=2.5737 and σ4=0.9018 are the non-zero eigenvalues of Q. By applying the proposed method in Section 3 in determining the constant coefficients, we can obtain the following results for Equation (Equation18):

For j = 0, (19) Ψ20(t)=0.1501e6.2507t0.2039e3.7738t0.13e2.5737t+0.1838e0.9018t.(19)

For j = −1, (20) Ψ21(t)=0.3254e6.2507t+0.1052e3.7738t0.037e2.5737t+0.2572e0.9018t.(20)

For j = −2, (21) Ψ22(t)=0.1879e6.2507t+0.4237e3.7738t+0.225e2.5737t+0.1634e0.9018t.(21)

For j = −3, (22) Ψ23(t)=0.0193e6.2507t0.1825e3.7738t+0.1758e2.5737t+0.026e0.9018t.(22)

For j = 1, (23) Ψ21(t)=1+0.0066e6.2507t0.1425e3.7738t0.2337e2.5737t0.6305e0.9018t.(23) Using Equations (Equation19)–(Equation23), for t=0.5,1.0,,5, the time-dependent transition performance probabilities for the circular repairable C(2,6:F) system can be obtained (see, Table ). In addition, using Equation (Equation23) for i = −2 in Equation (Equation8), we can obtain the time-dependent reliability of the system when two components failed at time t.

Table 2. Time-dependent transition performance probabilities for the circular repairable C(2,6:F) system when two components fail.

For μ=0, we have the circular irreparable C(2,6:F) system and the Q-matrix becomes Q˘=(6λ6λ00005λ3λ02λ004λ23λ103λ0003λ3λ00000). Suppose that all the five components are working at time t, for i = 0 and j = 1 in Equation (Equation6), we have (24) Ψ01(t)=A001+A101eσ1t+A201eσ2t+A301eσ3t+A401eσ4t,(24) where A001, A101, A201, A301 and A401 are constant coefficients and σ1=6λ, σ2=5λ, σ3=4λ and σ4=3λ are the non-zero eigenvalues of Q˘. To determine the constant coefficients, we have to solve the following system of differential equations: A001+A101+A201+A301+A401=0,(6λ)A101+(5λ)A201+(4λ)A301+(3λ)A401=0,(6λ)2A101+(5λ)2A201+(4λ)2A301+(3λ)2A401=12λ2,(6λ)3A101+(5λ)3A201+(4λ)3A301+(3λ)3A401=72λ3,(6λ)4A101+(5λ)4A201+(4λ)4A301+(3λ)4A401=228λ4. The solutions of the above equations are A001=1, A101=2, A201=6, A301=3, A401=2. Then, Equation (Equation24) can be obtained as (25) Ψ01(t)=12e6λt+6e5λt3e4λt2e3λt.(25) Finally, using Equation (Equation25) for i = 0 in Equation (Equation8), the time-dependent reliability of the system can be obtained as (26) R0(t)=1Ψ01(t)=2e3λt+3e4λt6e5λt+2e6λt.(26) Note that Equation (Equation26) is equivalent to Equation (42) in Zhang et al. (Citation2000).

4.3. A linear repairable C(3,8:F) system

As another example, suppose we have a linear repairable C(3,8:F) system. In this case the state space is Ω={0,1,2,3,4,5,6,1}, the set of working state is W={0,1,2,3,4,5,6} and the set of failed states is F={1}. The Q-matrix is Q=(8λ8λ00μ(μ+7λ)7λ00μ(μ+6λ)7514λ00μ(μ+5λ)000μ00000000000000000000000914λ185λ0075λ(μ+4λ)169λ0209λμ(μ+3λ)38λ218λ0μ(μ+2λ)2λ0000). Suppose that all the eight components are working at time t, if λ and μ are given, we can obtain the time-dependent reliability of the system by using the method proposed in Section 2. Here, we consider three sets of parameters: (i) λ=0.5,μ=1.5; (ii) λ=0.5,μ=3.0; and (iii) λ=0.25,μ=1.5.

Let, i = 0 in Equation (Equation8), we have (27) R0(t)=1Ψ01(t)=1(A001+A101eσ1t+A201eσ2t++A701eσ7t),(27) where A001, A101,…, A701 are constant coefficients and σ1, σ2,…, σ7 are the non-zero eigenvalues of Q. Using the above Q-matrix for each parameter sets respectively, the non-zero eigenvalues of Q are computed and they are presented in Table .

Table 3. The non-zero eigenvalues of Q for three parameter sets.

Using eigenvalues as seen in Table  and Equation (Equation7) for k=0,1,2,,7, we can obtain the constant coefficients (see, Table ). Using the constant coefficients in Tables  and , we plotted the reliability R0(t) against time t for the linear repairable C(3,8:F) system in Figure  for three different parameter sets.

Figure 1. The plot of the reliability R0(t) of the linear repairable C(3,8:F) system against time t.

Figure 1. The plot of the reliability R0(t) of the linear repairable C(3,8:F) system against time t.

Table 4. The constant coefficients in Equation (Equation27) for the three parameter sets.

4.4. A circular repairable C(3,8:F) system

In this section, we study a circular repairable C(3,8:F) system. The state space is Ω={0,1,2,3,4,5,1}, the set of working state is W={0,1,2,3,4,5} and the set of failed states is F={1}. The Q-matrix is Q=(8λ8λ00μ(μ+7λ)7λ00μ(μ+6λ)7514λ00μ(μ+5λ)000μ00000000000000000000914λ185λ075λ(μ+4λ)169λ209λμ(μ+3λ)218λ0μ2λ000). Suppose that all the eight components are working at time t, to evaluate the accuracy of the proposed method, we used the three parameter sets described in Section 4.3. Let i = 0 in Equation (Equation8), we have (28) R0(t)=1Ψ01(t)=1(A001+A101eσ1t+A201eσ2t++A601eσ6t),(28) where A001, A101,…, A601 are constant coefficients and σ1, σ2,…, σ6 are the non-zero eigenvalues of Q. Using the Q-matrix for each of the three parameter sets, the non-zero eigenvalues of Q are obtained and presented in Table .

Table 5. The non-zero eigenvalues of Q for three parameter sets.

Using the eigenvalues presented in Table  for k=0,1,2,,6, we can obtain the constant coefficients in Equation (Equation28) (see, Table ). When the values in Tables  and  are used in Equation (Equation28), we plotted the reliability R0(t) against time t for the circular repairable C(3,8:F) system in Figure  based on the three parameter sets.

Figure 2. The plot of the reliability R0(t) of the circular repairable C(3,8:F) system against time t.

Figure 2. The plot of the reliability R0(t) of the circular repairable C(3,8:F) system against time t.

Table 6. The constant coefficients in Equation (Equation27) for the three parameter sets.

5. Concluding remarks

In this paper, a repairable C(k,n:F) system with n linearly and circularly arranged components is studied. A novel and simple method for determining the time-dependent transition performance probability and time-dependent reliability of the system when the state of the system is known within a certain of time is developed. Some numerical examples are used to illustrate the usefulness of the proposed method. Computer code written in R Core Team (Citation2019) to compute the time-dependent transition performance probability can be obtained from the authors upon request. The proposed method can serve as a simple alternative to the Laplace transform and the Runge-Kutta method which are commonly used to determine the time-dependent reliability of the repairable systems.

Disclosure statement

No potential conflict of interest was reported by the author(s).

Additional information

Funding

H. K. T. Ng's work was also supported by a grant from the Simons Foundation [Grant Number 709773].

Notes on contributors

Gökhan Gökdere

Gökhan Gökdere is currently an Associate Professor of Statistical Science with the Fırat University, Elazıg, TURKEY. His research interests include reliability and statistical inference.

Hon Keung Tony Ng

Hon Keung Tony Ng is currently a Professor of Statistical Science with the Southern Methodist University, Dallas, TX, USA. His research interests include reliability, censoring methodology, ordered data analysis, nonparametric methods, and statistical inference. He is a Fellow of the American Statistical Association, an elected senior member of IEEE, and an elected member of the International Statistical Institute.

References