CC BY
11
A FUZZY RELIABILITY MODEL FOR "SAFETY SYSTEM-PROTECTED OBJECT" COMPLEX
A. I. Pereguda, D. A. Timashov
Obninsk Institute for Nuclear Power Engineering, Obninsk, Russia e-mail: [email protected]
ABSTRACT
The paper presents a new fuzzy reliability model for automated "safety system-protected object" complex. It is supposed that parameters of reliability model and reliability indices are fuzzy variables. Scheduled periodic inspections of safety system are also taken into account. Asymptotic estimates of mean time to accident membership function are proposed.
1 INTRODUCTION
Many researchers applied the concept of fuzzy reliability on various systems (Cai et al 1991, Cai et al. 1993, Cai 1996, Chen & Mon 1993, Onisava & Kacprzyk 1995, Utkin & Gurov 1995, Verma et al. 2007). Those researches are based on possibility instead of probability assumption or fuzzy state instead of binary state assumption. This paper presents a slightly different concept of fuzzy reliability. It is supposed that parameters of the reliability model are fuzzy variables. According to the random fuzzy variables theory presented by Liu (Liu 2002) reliability indices in this case are also fuzzy variables. In order to develop a complete practical methodology of the fuzzy reliability assessment we also consider some aspects of the fuzzy parameter estimation and numerical methods of the fuzzy arithmetic.
In the present study we set out to analyze the reliability of the automated "safety system-protected object" complex. Systems of such kind are quite common in the nuclear power engineering. We follow Pereguda (Pereguda 2001) in assuming that the operation of the complex can be described using a superposition of alternating renewal processes. We also utilize the concept of a random fuzzy renewal process (Shen et al. 2008, Zhao et al. 2006). Our objective is to provide an asymptotic estimation for the mean time to accident membership function.
2 MODEL DESCRIPTION
Let us consider an automated complex of safety system and protected object. The safety system and the protected object are repairable. They are restored to an as-good-as-new state. It is assumed that failures of the safety system can be detected only during periodic inspections of the safety system. All failures are supposed to be independent.
Let {F (t;X(0)),0e&} be a family of probability distributions on the probability space
(Q, A, P) with a common fuzzy parameters vector X on the credibility space (0, n, Cr) which induces a joint membership function, ^(x), then x is a random fuzzy variable (Guo et al. 2007). For example, x is an exponentially distributed random fuzzy variable, x~EXP(A) if
[1 - e-A(0)t,ift > 0 I 0, otherwise.
Fx(t;m) =
If x is a random fuzzy variable defined on the credibility space (0, n, Cr) than the probability Pr(x e A) is a fuzzy variable for any Borel set A e R and the expected value E[x] is a fuzzy variable provided that E[x(0)] is finite for each 0 e & (Liu 2002).
By X, i = 1,2,... denote the time to the i-th protected object failure. Let x, i = 1,2,... be independent and identically distributed (i.i.d) random fuzzy variables (Li & Liu 2006) with CDF Fx(t; Xx(d)). By y, i = 1,2,. denote the time to the protected object repair after it's i-th failure. Let Y, i = 1,2,. be i.i.d. random fuzzy variables with CDF FY(t; ^(6*)). Suppose that moments of the protected object repair are renewal points of the operation process of the complex. By i = 1,2,. denote the time to the i-th failure of the safety system. Let i = 1,2,. be i.i.d. random fuzzy variables with CDF F^t; X,^(6>)). By n, i = 1,2,. denote the time to the safety system repair after it's i-th failure. Let n, i = 1,2,. be i.i.d. random fuzzy variables with CDF Fn(t; X,n(6>)). Suppose that moments of the safety system repair are renewal points of the operation process of the safety system. By T denote the period of scheduled inspections of the safety system. By S denote the duration of scheduled inspections of the safety system. The safety system is inactive during the inspection. By v denote the number of renewal intervals before the accident. Let v be an integer random fuzzy variable. By rn denote the time to accident. An accident takes place when the protected object fails during the period of the safety system inactivity. Our aim is to estimate the membership function fj.Mm(y) of the mean time to accident.
2 MAIN RESULTS
Since the operation process of the complex is a superposition of alternating renewal processes, it follows that
v — 1
® = Z (X + Y ) + Xv •
i=1
Taking into account the fact that E[®(0)], E[#(0)], E|>(0)], E[£(0)], and E[v(6>)] are crisp
variables for each fixed 0 e ® we obtain
E[a>{9)] = E
'v(e)-\
Kx (0))+ xv[e)(0)
i=\
for each fixed 0 e & . Since all random fuzzy variables of interest are independent it follows that
Pr| Z(x(0)+^(0)) + Xv(o)(0)>r
i=\
k-\
K Pr(v(0) = k)Prl K ( (0) + Y (0)) + Xk (0) > r
k=\
- i=\
Note that
<» f v(0)-\
E[ (0)] = J Prl K ( (0) + Yi (0)) + Xv(o) (0)
>r
dr
k-\
(0) + Yi (0)) + Xk (0)
= X Pr(v(0) = k )E
k=1 L i=1
Since all random fuzzy variables of interest are i.i.d. it follows that
Xi
E[®(0)]=YPr(v(0)= k)[kx(0)+ (k -1)(0)] = E[v(0)x(0)+ (y(0)~ )))
k=1
Therefore
E[®]= E[vx + (v - 1)7].
Taking into account the fact that all random fuzzy variables of interest are independent we obtain
00
i=\
0
E[o(e)]=-^fXe)] ^ eW)]
= ^ - x >{0)+^^ - X^)
for each fixed , where q(ff) is the probability of the accident during a renewal interval.
Let Q(0) be the set of intervals where the safety system is active and let Q (O) be the set of intervals where safety system is inactive. We obviously have
ro
q(0) = JPr( e Q~(0))dFz(f;X»)
0
for each fixed 0 e® . Note that
Pr (t e Q- (0)) = 1 - Pr (t e Q+(0)) = 1 - Pr + (t; 0). Applying the law of total probability we obtain
ro ro
Pr+(t;0)= {{Pr( e Q +(0)| ^(0)= x,n(0) = y) (;Xv(0))dF,(x;X,(0)).
0 0
Since the operation process of the safety system is an alternating renewal process, it follows that Pr + (t;0) = {{Pr(( e Q+(0) | ,(0) = x, n(0) = y) (y; Xv (0^ (x; X, (0))
Tss (x,y)<t
+ {{Pr( e Q+ (0) | ,(0) = x,r,(0) = y)(;Xn(0)(x;X,(0)) = I1 +12,
Tss (x,y)>t
where zss (,,n) = 1 ^ +1( + 8) + n is the length of the renewal interval of the safety system
V\T + 8/ J
operation process and <x> is an integer part of x. We see that
'2 = {{
T+8
Tss (x,y)>t
£ Jte[m(T +S),m(T+S)+T) + J VI x ( . m = 0 H + 8)
dFn(y; Xn(0))dF,(x; X,(0)),
where JA is an indicator function of the event A. It now follows that
I2 = (1 - F,(t;X,(0)))-X (1 - F,(m(T + 8);X,(0))\j(m^+s)+t<, - Jm(T+s),t)
m=1
= Fz(t; XM-F,(t; X,(0)),
where
Fz(;X,(0))= 1 - £ (1 - F,(m( + 8);X^))-1)(t+s)+t<t - Jm(T+8)t).
m =1
Note that
t
I1 = {{ Pr + (t(x,yddF^;Xn(0)]dF,(x;X,(0)) = {Pr + (t -z;0)dFTss (z;X»,X,(0)).
T.SS (x,y)<t 0
Finally,
Pr+(t; 0) = f (t; 0) + { Pr+(t - z; 0)dFTss (z; Xv (0), X, (0)),
0
where f (t;0) = Ff(t; X,(0))-F,(t; X,(0)) and F^ (z; X» X,(0)) = di^y ^ + 8) + n()< z^ The application of Laplace-Stieltjes transform and tauberian theorems yields
limpr+M)=e m]-e №)]
SPr E[(0)]
where
e[z(0)]=se
T + Si
= & k (f( + 1)(T + S); F,(k (F + S); k()),
k=1
X
E ,(0)]=JF - f,(f \(e))t,
E[rss(0)] = J(l - Fn(kK fa)))* + (T + S) 1 + XkF((k + 1)T + S);k()- F,(k(k + S);^(0))) I.
0 V k=1 )
Therefore
01 - ^^M]1 E [ss (0)]
Taking into account the definition of a function of fuzzy variables and Zadeh Extension Principle we obtain
ae[*]F) = suP min(A ,(x1)> (x2^ A, (x3)> (x4)) '
y=f(x1.
where
f (X1' x 2' X3. x 4 ) = -Tk1--1 J(1 - Fx(t; X1 ))dt + 1 "F(x3x x)4) J(1 - Fy(t; x 2)))
q\x.3,2.4)0 0
a(x x ^ 1 - E[,(x3)]-E[Z(x3)]
E[tSS (x3' x4 )]
E ZF3)] = SX k F (F + 1)(T + S); x 3) - F, (k (k + S); x 3)
k=1
X
E [,(x 3)] = J (1 - F,(t; x 3 ))t
it
J \ , \ y 3//
0
X fx \
E[s (x3, x 4 )]=J(1 - Fv(t; x 4 ))t + (T + S) 1 + £ k ((( + 1)(k + S); x3)-F,(k (T + S); x3 ))
0 V k=1 )
Therefore it is now possible to estimate the fuzzy mean time to accident. In order to perform defuzzyfication we use expected value operator suggested by Liu (Liu & Liu 2003):
EH = 2 JI sup Am* (y) +1 - sup Am* F) \dr ■
2 0 V y - r y < r )
It is clearly evident that the most difficult part in the proposed methodology is the evaluation of p.Mm(y) according to Zadeh Extension Principle. To overcome these difficulties one should use a suitable numerical method. We suggest to use General Transformation Method (Hanss 2005) for this task. This method consists of a decomposition of input fuzzy variables, a transformation of the input intervals, an evaluation of the model and a retransformation of the output array. Another important issue is the fuzzy parameter estimation. We suggest to use fuzzy estimators developed by Buckley (Buckley 2006). These estimators are based on confidence intervals and allow the estimation of the membership function of the distribution parameter from the sample data.
Consider now the following example. Suppose /~EXP(AZ), y~EXP(^), ^-EXP(A^), ^-EXP(^). Therefore
Ae [*](y )= sup min(Ax(*1 )'Ar(x2 )'Ai(x3)' An(x4))'
y = f F'x2'x3'x4 )
where
X
f(x x x x )=_1_l+1 - q(x3,x4) 1
J V 1' 2' 3' 4 / ( ) ^ ( )
^x3 5 x4 / xi ^x3 5 x4 7 x2
q(x x . - E[¿(x.)]-E[Z(x3)]
^x3J xJ ~ 1 E[T (x x )] 5
^ L' SS 3' A4 j\
E [C(x3 )] = S
- x3(T + S)
1 - e
-x3 (T + S)
E [fe )] =
E[TSS (x3, x4 )]= — + (T + S)
1 + •
e
- x3 (T + S) \
1 - e~x.(T+S)
v I H y
y 107h
1.5 2.0 2.5 3.0 3.5 4.0 Figure 1. Membership function of the mean time to accident.
Let all fuzzy parameters be triangular fuzzy variables: /uA (x) = a(1 x 10 6h :,1.5 x 10 6h ',2 x 10 6h 1), (x) = a(1h-1,1.5h-1, 2h-1), p (x) = a(1x10-4h_1,1.5 x10-4h"1, 2 x10-4h-), ^(x) = a^h"1,1.5h_1, 2h-), T = 500h, S = 0.1h. Membership function of the mean time to accident is shown on Figure 1.
3 CONCLUSIONS
The proposed model permits to assess the reliability of one specific class of technological systems with fuzzy parameters. In particular the suggested approach allows to evaluate the membership function of the mean time to accident for the "safety system-protected object" complex. The proposed approach allows to take into account the uncertainty of reliability model parameters and reliability indices. The solution obtained is useful for reliability assessment of nuclear power plants and similar dangerous technological objects.
e
4
REFERENCES
Buckley, J.J. 2006. Fuzzy Probability and Statistics. Springer-Verlag.
Cai, K.Y., Wen, C.Y. & Zhang, M.L. 1991. Fuzzy variables as a basis for a theory of fuzzy reliability in the possibility
context. Fuzzy Sets Syst. 42: 145-172. Cai, K.Y., Wen, C.Y. & Zhang, M.L. 1993. Fuzzy states as a basis for a theory of fuzzy reliability. Microelect. Reliabil. 33: 2253-2263.
Cai, K.Y. 1996. Introduction to Fuzzy Reliability. Norwell: Kluwer Academic Publishers.
Chen, C.-H. & Mon, L. 1993. Fuzzy system reliability analysis by interval of confidence. Fuzzy Sets Syst. 56: 29-35.
Guo, L., Zhao, R.Q., Guo, D. & Dunne, T. 2007. Random fuzzy variable modelling on repairable system. Journal of
Uncertain Systems 1, no. 3: 222-234. Hanss, M. 2005. Applied Fuzzy Arithmetic: An Introduction with Engineering Applications. Springer-Verlag. Li, X. & Liu, B. 2006. New independence definition of fuzzy random variable and random fuzzy variable. World
Journal of Modelling and Simulation 2, no. 5: 338-342. Liu, B. 2002. Theory and Practice of Uncertain Programming. Physica Verlag Heidelberg.
Liu, Y. & Liu, B. 2003. Expected value operator of random fuzzy variable and random fuzzy expected value models.
Int. J. Uncertainty, Fuzziness Knowl.-BasedSyst. 11: 195-215. Onisava, T. & Kacprzyk, J. 1995. Reliability and safety analysis under fuzziness. Physica Verlag Heidelberg. Pereguda, A.I. 2001. Calculation of the Reliability Indicators of the System Protected Object-Control and Protection
System. Atomic Energy 90: 460-468. Shen, Q., Zhao, R. & Tang, W. 2008. Random fuzzy alternating renewal processes. Soft Computing - A Fusion of
Foundations, Methodologies and Applications 13, no. 2: 139-147. Utkin, L.V. & Gurov, S.V. 1995. A general formal approach for fuzzy reliability analysis in the possibility context.
Fuzzy Sets Syst. 83: 203-213. Verma, A.K. et al. 2007. Fuzzy-Reliability Engineering: Concepts and Applications. Narosa Publishing House. Zhao, R., Tang, W. & Yun, H. 2006. Random fuzzy renewal process. European Journal of Operational Research 169: 189-201.
