Home Exponential Stability and Availability Analysis of a Complex Standby System
Article Publicly Available

Exponential Stability and Availability Analysis of a Complex Standby System

  • Chunli Wang , Mingyan Teng and Fu Zheng EMAIL logo
Published/Copyright: June 25, 2014
Become an author with De Gruyter Brill

Abstract

In this paper, we first investigate the solution of the two correlated units redundant system with two types of failure. By using the method of functional analysis, especially, the c0 semigroup theory of bounded linear operators on Banach space, we prove the well-posedness and the existence of positive solution of the system. By analyzing the spectra distribution of the system operator, we prove that the dynamic solution of the system asymptotically converges to the nonnegative steady-state solution which is the eigenfunction corresponding to eigenvalue 0 of the system operator. Furthermore, we discuss the exponential stability of the system. Finally, we analyze the reliability the system with the help of our main results and present some reliability indices of the system at the end of the paper.

1 Introduction

Reliability is an important concept of various complex systems. The need of obtaining highly reliable system has became more and more important with the development of the modern technology. The high degree of reliability is usually achieved by introducing redundancy. The redundant system discussed here comes from [1]. It consists of two sub-systems A and B, connected in parallel having one regular and one expert repairman. This system has four states, viz. normal, degraded, failed and waiting. The sub-system B is in standby mode and consists of further two units, viz. B1 and B2, connected in parallel. It is assumed that when both the units of sub-system B start functioning simultaneously, then sub-system B works like the sub-system A. So, when some failure occurs in sub-system A, the two units of sub-system B will be in operable condition through perfect switching over device. If any one unit of subsystem B fails (after failing sub-system A), the system will work but in degraded stage. The system will be in down state only when both sub-system A and B fail or due to occurrence of environmental failure.

In [1], the well-posedness and asymptotic stability of the redundant system have been obtained. Because it is difficult to determine the convergence rate of the dynamic solution only by asymptotic stability, the study of the exponential stability is in demand, which is a more valuable property in application. We will continue to discuss the above system and obtain the exponential stability of the system. It should be pointed out that, in recent years, though there are several papers (e.g. see [3-6]) on the exponential stability of the repairable systems, the method of [2-5] is involved in quasi-compactness and is hard to be applied here. Moreover, as the applications of the main results of [1], we will calculate the stationary availability of the system with another method different from the traditional one [6-9].

The rest of this paper is organized as follows. In section 2, we introduce the mathematical model and some notations of the system. In the meantime, some results on the well-posedness and stability are introduced in this section. In section 3, the exponential of the system is presented. In the last section, we investigate the stationary availability of the system with another method different from the traditional one and give several numerical simulation examples.

2 Mathematical model formulation and some results

Most materials of this section come from [1] and for convenience, we do not restate and introduce them freely. Let a = λ0 + λ1 + λ2, p0(t) be the probability that the system is in operable state at time t, p1(t) probability that the system is in waiting state at time t, pi(t, x)dx probability that the system is in operable state at time and elapsed repair time lies between x and x + dx, where j = 8, 9, 10, 11, and pi(t, x)dx probability that the system is in failed state at time t and elapsed repair time lies between x and x + dx, where i = 7, 12. With some additional assumptions, the differential equations of the system are obtained,

(1)p0(t)=0p2(t,x)μ1(x)dx+0p3(t,x)μ2(x)dx+0p4(t,x)μ(x)dx+0η(x)[p5(t,x)+p6(t,x)]dx+0p7(t,x)μ3(x)dx+0p12(t,x)β(x)dx
(2)[ddt+r1]p1(t)=λ00p4(t,x)dx
(3)[t+x+μ1(x)+λ2]p2(t,x)=0
(4)[t+x+μ2(x)+λ1]p3(t,x)=0
(5)[t+x+μ(x)+λ]p4(t,x)=0
(6)[t+x+η(x)+λ1+λ2+r2]p5(t,x)=0
(7)[t+x+η(x)+λ1+λ2]p6(t,x)=0
(8)[t+x+μ3(x)]p7(t,x)=0
(9)[t+x+μ1(x)+r2]p8(t,x)=0
(10)[t+x+μ2(x)+r2]p9(t,x)=0
(11)[t+x+μ1(x)+λ2]p10(t,x)=0
(12)[t+x+μ2(x)+λ1]p11(t,x)=0
(13)[t+x+β(x)]p12(t,x)=0

and with the boundary conditions

(14)p2(t,0)=λ1p0(t),p3(t,0)=λ2p0(t)
(15)p4(t,0)=λ20p2(t,x)dx+λ10p3(t,x)dx
(16)p5(t,0)=λ0p0(t)+r1p1(t)+0μ1(x)p8(t,x)dx+0μ2(x)p9(t,x)dx
(17)p6(t,0)=r20p5(t,x)dx+0μ1(x)p10(t,x)dx+0μ2(x)p11(t,x)dx
(18)p7(t,0)=λ3p0(t),p8(t,0)=λ10p5(t,x)dx,p9(t,0)=λ10p5(t,x)dx
(19)p10(t,0)=λ10p6(t,x)dx+r20p8(t,x)dx
(20)p11(t,0)=λ20p6(t,x)dx+r20p9(t,x)dx
(21)p12(t,0)=λ20p10(t,x)dx+λ10p11(t,x)dx

and the initial conditions p0 = 1; p1 = p1(0, x) = · · · = p2(0, x) = 0. You can refer to [1] for the concrete meanings of the symbols, such as λ0 and r1 etc.

The equations (1)~(21) can be rewritten as an the abstract Cauchy problem (see [1112]) in Banach space X:

(ACP)dp(t)dt=Ap(t),t0,p(0)=(1,0,,0)T.

in which 𝓐 is defined in [4] and X = ℝ2×(L1(ℝ+))11 is Banach space. For p = (p0, p1, p2 (x), · · ·, p12(x)) ∈ X, The norm of p is given by

||p||=|p0|+|p1|+i=2120|pi(x)|dx.

Let µ4(x) = µ(x), µ5(x) = β(x), µ6(x)= η(x). The main results of this paper are obtained under the following assumptions.

General Assumptions

There exist positive constants H and c, such that for any t ≥ 0,

1tte0uμi(s)dsdu>H,

and

c=min{infxR+μi(x),i=1,2,,6}.

The physical meaning of the general assumption is the mean repair rate is not too small. By the general assumptions, it is easily to obtain the following identify

0ui(x)e0xui(τ)d(τ)dx=1.

The following theorems are main results of [4].

Theorem 1

𝓐 generates a positive C0 contractive semigroup T(t) on X and the unique time dependent positive solution to the equations (1)~(21) exists.

Theorem 2

Zero is an algebraically simple eigenvalue of 𝓐 and the system is asymtotically stable. The steady-state solution is the positive unit eigenfunction corresponding eigenvalue 0 of the generator 𝓐.

3 The exponential stability of the system

Let T(t) be the semigroups generated by the operator 𝓐 and p(0) the initial value, the asymptotic stability shows that

T(t)p(0)=(p(0),Q)p~+T1(t)p(0)=p~+T1(t)p(0),limtT1(t)p(0)=0,

in which is the positive unit eigenfunction corresponding to eigenvalue 0 of 𝓐 and Q = (1, · · · , 1) ∈ x*.

In general, the above stability does not implies the exponential stability, that is to say there exist positive ω and M such that

||T1(t)p(0)||Meωt,t0.

But in some special cases this implication holds. Now we will show the exponential stability of the system. First, we introduce some notations and results which come from references [3] and [14].

Let X be a Banach space and T(t) be a c0 semigroups on X and 𝓐 its generator, for λρ(𝓐), then we have the symbols

ω1(T)=limt+ln||T(t)R(λ,A)||t,s0(A)=sup{ω>s(A):supReλ=ω||R(λ,A)||<+},

in which ω1(T) is the growth bound of the semigroups T(t) and s0(A) is abscissa of uniform boundedness of the resolvent of 𝓐. Their relationship is given by the following lemma which comes from [3].

Lemma 1

Let T(t) be a c0 semigroup on Banach space X and 𝓐 its generator, then ω1(T) < s0(𝓐).

The following Lemma comes from [14].

Lemma 2

For each yX, b ∈ ℝ and ib ∈ ℂ, the unique solution to the equations (ibA)p = y exists and p is given by

(22)(a+ib)p0=0p2(x)μ1(x)dx+0p3(x)μ2(x)dx+0p4(x)μ(x)dx+0η(x)[p5(x)+p6(x)]dx+0p7(x)μ3(x)dx+y0
(23)(r1+ib)p1=λ00p4(x)dx+y1
(24)p2(x)=p2(0)e0x[ib+λ2+μ1(s)]ds+0eτx[ib+λ2+μ1(s)]dsy2(τ)d(τ)
(25)p3(x)=p3(0)e0x[ib+λ1+μ2(s)]ds+0eτx[ib+λ1+μ2(s)]dsy3(τ)d(τ)
(26)p4(x)=p4(0)e0x[ib+λ0+μ(s)]ds+0eτx[ib+λ0+μ(s)]dsy4(τ)d(τ)
(27)p5(x)=p5(0)e0x[ib+λ1+λ2+r2+η(s)]ds+0eτx[ib+λ1+λ2+r2+η(s)]dsy5(τ)d(τ)
(28)p6(x)=p6(0)e0x[ib+λ1+λ2+η(s)]ds+0eτx[ib+λ1+λ2+η(s)]dsy6(τ)d(τ)
(29)p7(x)=p7(0)e0x[ib+μ3(s)]ds+0eτx[ib+μ3(s)]dsy7(τ)d(τ)
(30)p8(x)=p8(0)e0x[ib+r2+μ1(s)]ds+0eτx[ib+r2+μ1(s)]dsy8(τ)d(τ)
(31)p9(x)=p9(0)e0x[ib+r2+μ2(s)]ds+0eτx[ib+r2+μ2(s)]dsy9(τ)d(τ)
(32)p10(x)=p10(0)e0x[ib+λ2+μ1(s)]ds+0eτx[ib+λ2+μ1(s)]dsy10(τ)d(τ)
(33)p11(x)=p11(0)e0x[ib+λ1+μ2(s)]ds+0eτx[ib+λ1+μ2(s)]dsy11(τ)d(τ)
(34)p12(x)=p12(0)e0x[ib+β(s)]ds+0eτx[ib+β(s)]dsy12(τ)d(τ)
(35)p2(0)=λ1p0,p3(0)=λ2p0
(36)p4(0)=λ20p2(x)dx+λ10p3(x)dx
(37)p5(0)=λ0p0+r1p1+0μ1(x)p8(x)dx+0μ2p9(x)dx
(38)p6(0)=r20p5(x)dx+0μ1p10(x)dx+0μ2p11(x)dx
(39)p7(0)=λ3p0,p8(0)=λ10p5(x)dx,p9(0)=λ20p5(x)dx
(40)p10(0)=λ10p6(x)dx+r20p8(x)dx
(41)p11(0)=λ20p6(x)dx+r20p9(x)dx
(42)p12(0)=λ20p10(x)dx+λ10p11(x)dx

Theorem 3

The system is exponentially stable.

Proof

Let one rank projection P be the associated spectral projection of the eigenvalue 0 since 0 is an algebraically simple eigenvalue of 𝓐. Thus X has the invariant decomposition

X=X0+Y,X0=PX,Y=(IP)X,

in which X0 is the eigenfunction subspace whose dimension is one. Y is also Banach space equipped with the norm of X. Let T1(t) and A1 be the restrictions of T(t) and A on Y respectively. Obviously, we have σ(A1) = σ(A) \ {0}. Now we will show that s0(A1) ≤ 0. However, 0 is a regular point of A1, so it is sufficient to show that

supγ=ib,|b|>δ>0,bR||R(λ,A1)||<+.

In fact, for r = ib, b ∈ ℝ and each given yX, from Theorem 2 the equation (rIA)p = y has solution pD(A) and pi(x) is given by (17) and (18). Thus we have

(43)|b|0|pi(x)|dx0|pi(0)re0x[r+γi+hiξ]dξ|dx+|b|0|0yi(τ)eτx[r+γi+hi(ξ)d(ξ)]|dx
(44)|b|0|p2+i(x)|dx0|p2+i(0)re0x[r+hiξ]dξ|dx+0|pi(0)r(e0x[r+hiξ]dξre0x[r+γi+hiξ]dξ)|dx+|b|c0y2+i(x)dx+|b|0|0xyi(τ)(eτx[r+hi(ξ)]dξeτx[r+γi+hiξ]dξ)dτ|dx

From (13), we have the following inequality

(45)(|b|+γ)p0y0+i=12|0pi(0)hi(x)e0(r+γi+hi(ξ))dξdx|+i=120|yi(x)|dx

Moreover, by the boundary conditions (14)~(16) we have

(46)|p1(0)|λ1|p0|+0|p2(0)(1re0x[r+h2(ξ)]dξh2(x)e0x[r+h2(ξ)]dξ)|dx+0|p4(0)(1re0x[r+h2(ξ)]dξ)|dx+20|y2(x)|dx+0|y4(x)|dx
(47)|p2(0)|λ2|p0|+0|p1(0)(1re0x[r+h1(ξ)]dξh1(x)e0x[r+h1(ξ)]dξ)|dx+0|p3(0)(1re0x[r+h1(ξ)]dξ)|dx+20|y1(x)|dx+0|y3(x)|dx

If adding (43)~(47) together, we would have the following estimate

(48)|b|(|p0|+i=140|pi(x)|dx)y0+(1+|b|c)i=140|yi(x)|dx+2i=120|yi(x)|dx

By the definition of the norm on the space X and (48), we obtain that

||R(rIA1)||1|b|(2+|b|c).

Therefore we know that,

limb+||R(rIA1)||=1c,

and

supr=ib,|b|>δ>0,bR||R(rIA1)||<+.

From the definition of s0(A) and Lemma 2, we obtain that s0(A1) ≤ 0 and ω1(T1) < 0 due to ω1(T1) < s0 (A1). However, by [14] we know semigroups T(t) is exponentially stable if and only if ω1(T) < 0. Thus we complete the proof of Theorem 3.

In fact, the exponential stability of the repairable system was firstly given in [14]. You can refer to the reference [14] in the modified version for more information. To show the effectiveness of the reliability of the repair system, we include this part again in this note. The method is mainly employed from [3].

4 Some steady-state reliability indices

According to [2],

p0(t)+i=260pi(x,t)dx+i=8110pi(x,t)dx

and

AV=limt(p0(t)+i=260pi(x,t)dx+i=8110pi(x,t)dx)

are the instantaneous availability and the stationary availability of the repairable system, respectively. The availability of the system is one of the most important reliability indices and engineers are especially interested in the steady-state availability. Now we will obtain the stationary availability of the system with another method different from the traditional one. For convenience, set

f2(x)=e0x[μ1(s)+λ2]ds,a2=0μ1(x)f2(x)dx,b2=0f2(x)dx,f3(x)=e0x[μ2(s)+λ1]ds,a3=0μ2(x)f3(x)dx,b2=0f3(x)x,f4(x)=e0x[μ(s)+λ0]ds,a4=0μ(x)f4(x)dx,b4=0f4(x)dx,f5(x)=e0x[η(s)+λ1+λ2+r2]ds,a5=0η(x)f5(x)dx,b5=0f5(x)dx,f6(x)=e0x[η(s)+λ1+λ2]ds,a6=0η(x)f6(x)dx,b6=0f6(x)dx,f7(x)=e0xμ3(s)ds,b7=0f7(x)dx,f8(x)=e0x[μ1(s)+r2]ds,a8=0μ1(x)f8(x)dx,b8=0f8(x)dx,f9(x)=e0x[μ2(s)+r2]ds,a9=0μ2(x)f9(x)dx,b9=0f9(x)dx,f12(x)=e0xβ(s)ds,b12=0f12(x)dx.

Theorem 4

The stationary availability of the system is

(49)AV=1+i=46bici+i=89bici+b2(c10+λ1)+b3(c11+λ2)1+c1+i=46bici+λ3b7+i=89bici+b2(c10+λ1)+b3(c11+λ2)+b12c12

in which

c1=r11λ0λ1λ2b4(b2+b3),c4=λ1λ2(b2+b3),c5=[1b5(λ1a8+λ2a9)]1[λ0+λ0λ1λ2b4(b2+b3)],c8=λ1b5[1b5(λ1a8+λ2a9)]1[λ0+λ0λ1λ2b4(b2+b3)],c9=λ2b5[1b5(λ1a8+λ2a9)]1[λ0+λ0λ1λ2b4(b2+b3)],c6=r2(b5c5+a2b8c8+a3b9c9)[1b6(λ1a2+λ2a3)]1,c10=λ1b6c6+r2b8c8,c11=λ2b6c6+r2b9c9,c12=λ2b2c10+λ1b3c11,

Proof

In order to compute the stationary availability of the discussed system, we need to obtain the formulation of the positive unit eigenfunction p corresponding to eigenvalue 0 of the operator 𝓐. Now we complete this job. In fact, 𝓐p = 0 is the following set of equations:

(50)p2(x)=[μ1(x)+λ2]p2(x),p3(x)=[μ2(x)+λ1]p3(x),p4(x)=[μ(x)+λ0]p4(x),p5(x)=[η(x)+λ1+λ2+r2]p5(x),p6(x)=[η(x)+λ1+λ2]p6(x),p7(x)=μ3(x)p7(x),p8(x)=[μ1(x)+r2]p8(x),p9(x)=[μ2(x)+r2]p9(x),p10(x)=[μ1(x)+λ2]p10(x),p11(x)=[μ2(x)+λ1]p11(x)

with the initial conditions:

(51)r1p1=λ00p4(x)dx,p2(0)=λ1p0,p0(0)=λ2p0
(52)p4(0)=λ20p2(x)dx+λ10p3(x)dx
(53)p5(0)=λ0p0+r1p1+0μ1(x)p8(x)dx+0μ2(x)p9(x)dx
(54)p6(0)=r20p5(x)dx+0μ1(x)p10(x)dx+0μ2(x)p11(x)dx
(55)p7(0)=λ3p0,p8(0)=λ10p5(x)dx,p8(0)=λ20p5(x)dx
(56)p10(0)=λ10p6(x)dx+r20p8(x)dx
(57)p11(0)=λ20p6(x)dx+r20p9(x)dx
(58)p12(0)=λ20p10(x)dx+λ10p11(x)dx

Solving equations (50), we have

(59)pi(x)=pi(0)fi(x),i=2,3,,9,p10(x)=p10(0)f2(x)p11(x)=p11(0)f2(x),p12(x)=p12(0)f12(x)

By the Theorem 4.1 of [1], because 0 is a simple eigenvalue of 𝓐 and there exists a corresponding positive eigenvector, we can set p0 = 1.With the help of

0ui(x)e0xui(τ)dτdx=1

and the notations of ci, it is easy to know that

p1=c1,p2(0)=λ1,p3(0)=λ2,pi(0)=ci(i=4,5,6),p7(0)=λ3,

and

pj(0)=cj(j=8,9,,12).

Thus,

P=(1,c1,λ1f2(x),λ2f3(x),c4f4(x),c5f5(x),c6f6(x),λ3f7(x),c8f8(x),c9f9(x),c10f10(x),c11f11(x),c12f12(x))

is the positive eigenfunction corresponding to eigenvalue 0 of the operator 𝓐 and = ||P||–1P. Therefore, we can get the expression AV = 1 – (c1 + b7c7 + c12c12)||P||–1 by the definition of the stationary availability, that is the formula (49) due to

||P||=1+c1+i=46bici+λ3b+i=89bici+b2(c10+λ1)+b3(c11+λ2)+b12c12.

We complete the proof of the Theorem 4. With the help of Theorem 4, we are able to calculate the mean up-time of the redundant system.

Theorem 5

The mean up-time (MUT) of the redundant system is equal to

MUT=1+i=46bici+λ3b+i=89bici+b2(c10+λ1)+b3(c11+λ2)λ0(1+b4c4)+λ2b2(c10+λ1)+λ1b3(c11+λ2)+(λ1+λ1)(b5c5+b6c6)+λ32b7.

Proof

Let M be the rate of failures of the system, then

M=||p||1[λ0(1+p4)+λ2(p2+p10)+λ1(p3+p11)+(λ1+λ1)(p5+p6)+λ3b7].

From the proof of the Theorem 4, we know p2 = λ1b2, p3 = λ2b3, p7 = λ3b2, pi = cibi, i = 4, 5, 6, p10 = c10b2, and p11 = c11b3. However, the mean up-time (MUT) of the redundant system is AVM. Computing directly leads to the formulation of MUT.

Corollary 1

The stationary availability and MUT of the system without standby are

AV=11+λ3b7+λ0b12,MUT=1λ0+λ32b7.

Supported by the National Natural Science Foundation of China (Grant No. 11201037, 11371070)


Acknowledgments

The authors would like to thank the referees for their helpful suggestions and comments.

References

[1] Itoi T, Kodama M, Nishida T. Behavior of a two correlated units redundant system with many types of failure. Microeletronics Reliability, 1978, 17: 269–271.10.1016/0026-2714(78)90972-1Search in Google Scholar

[2] Cao J H, Cheng K. Introduction to reliability mathematics. Science Press, Beijing, China, 1986.Search in Google Scholar

[3] Wang W L, Xu G Q. Well-posedness and stability of solution of a series system with repair facility. Applied Mathematics A Journal of Chinese Universities (Ser A), 2007, 22: 474–482.Search in Google Scholar

[4] Wang W L, Xu G Q. Stability analysis of a complex standby system with constant waiting and different repairman criteria incorporating environmental failure. Applied Mathematical Modelling, 2009, 33: 724–743.10.1016/j.apm.2007.12.008Search in Google Scholar

[5] Hu W W, Shen Z F, Xin Y H, et al. Exponential stability of a repairable system with imperfect switching mechanism. Asymptotic Analysis, 2007, 54: 93–102.Search in Google Scholar

[6] Shen Z F, Hu X X, Fan W F, et al. Exponential asymptotic property of a parallel repairable system with warm standby under common-cause failure. Journal of Mathematical Analysis and Applications, 2008, 341: 457–466.10.1016/j.jmaa.2007.10.016Search in Google Scholar

[7] Guo L N, Xu H B, Gao C, et al. Stability analysis of a new kind series system. IMA Journal of Applied Mathematics, 2010, 75: 439–460.10.1093/imamat/hxq013Search in Google Scholar

[8] Gupur G, Li X Z, Zhu G T. Functional analysis method in queueing theory. Hertfordshire, Research Information Ltd: United Kingdom, 2001.Search in Google Scholar

[9] Wang J T. An M/G/1 queue with second potional service and server breakdowns. Computers and Mathematics with Applications, 2004, 47: 1713–1723.10.1016/j.camwa.2004.06.024Search in Google Scholar

[10] Haji A, Radl A. A semigroup approach to queueing systems. Semigroup Forum, 2007, 75: 609–623.10.1007/s00233-007-0726-6Search in Google Scholar

[11] Pazy A. Semigroup of linear operators and application to partial differential equations. Spring-verlag, New York, 1983.10.1007/978-1-4612-5561-1Search in Google Scholar

[12] Engel K J, Nagel R. A short course on operator semigroups. Springer, 2006.Search in Google Scholar

[13] Dunford N, Schwartz J T. Linear operators I. Wiley, New York, 1958.Search in Google Scholar

[14] Zheng F, Yang J Y, Zhu G T. Exponential stability of a complex standby system. Advanced Mathematics, 2012, 41(6): 713–722.Search in Google Scholar

[15] Dyson J, Villella-Bressan R, Webb G F. Asynchronous exponential growth in an age structured population of proliferating and quiescent cells. Mathematical Biosciences, 2002, 177–178: 73–83.10.1016/S0025-5564(01)00097-9Search in Google Scholar

[16] Arino O, Bertuzzi A, Gandolfi A, et al. The asynchronous exponential growth property in a model for the kinetic heterogeneity of tumour cell populations. Journal of Mathematical Analysis and Applications, 2005, 302: 521–542.10.1016/j.jmaa.2004.08.024Search in Google Scholar

Received: 2013-10-24
Accepted: 2014-3-6
Published Online: 2014-6-25

© 2014 Walter de Gruyter GmbH, Berlin/Boston

Downloaded on 20.11.2025 from https://www.degruyterbrill.com/document/doi/10.1515/JSSI-2014-0279/html
Scroll to top button