# Reliability and Sensitivity Analysis Method for a Multistate System with Common Cause Failure

Authors: Jinlei Qin; Zheng Li
Journal: Complexity
 (2019)

---

## Abstract

With the increasing complexity of industrial products and systems, some intermediate states, other than the traditional two states, are often encountered during reliability assessments. A system with more than two states is called a multistate system (MSS) which has already become a general phenomenon in the components and/or systems. Moreover, common cause failure (CCF) often plays a very important role in the assessment of system reliability. A method is proposed to assess the reliability and sensitivity of an MSS with CCF. Some components are not only in a failure state that can cause failure itself, but also in a state that can cause the failure of other components with a certain probability. The components that are affected by one type of CCF make up some sets which can overlap on some components. Using the technology of a universal generating function (UGF), the CCF of a component can be incorporated in the expression of its UGF. Consequently, indices of reliability can be calculated based on the UGF expression of an MSS. Sensitivity analysis can help engineers to judge which type of CCF should be eliminated first under various resource limitations. Examples illustrate and validate this method.


---

## Body

## 1. Introduction

Common cause failures (CCFs) are the dependent failures of multiple components originating from a common cause or single occurrence or condition. CCFs can be important contributors to system unavailability or accident risk. For example, lightning events can cause outages of unprotected electronic equipment and voltage surges caused by inappropriate switching can lead to multiple failures of components in a power system. Recognition of the fact that a CCF will increase the joint-failure probability and then reduce system reliability has inspired many researchers to model and estimate the reliability or availability of a system with CCF. Some basic concepts and theories have been investigated and developed [1].Two fundamental kinds of methods, termed implicit and explicit approaches, are used to incorporate CCFs in system reliability analysis [2]. The implicit approach proceeds as follows. The CCF is first ignored and the system logic is modeled with basic single-failure (component-level) events. The system reliability is given by algebraic probability expressions of the basic components, which are quantified to include the contribution of the CCF such that the system reliability or availability with the CCF can be expressed. This method has been studied in both binary- and multistate system (MSS) reliability [3]. In the explicit approach, the CCF is modeled as a basic cause-event in the system reliability block diagram or system fault tree, appearing as the repeated input to all elements or gates affected by the CCF. As the CCF can occur at random times in a redundant standby safety system, the expressions for the basic event probabilities of the explicit fault tree model have been developed [4].In the context of imperfect fault coverage (IPC), an uncovered failure may often lead to complete system failure [5]. Aiming at this situation, other methods such as ordered binary decision diagrams and generalized reliability block diagrams have been suggested to assess system reliability [6]. In the phased-mission system with IPC and CCF, a new binary decision diagram approach to reliability analysis was suggested [7]. Myers demonstrated how the coverage effect can be computed using combinational and recursive technologies for four coverage models [8]. In the reliability redundancy allocation problem three nonlinear optimization models mixing components with the inclusion of CCF events were addressed [9]. The reliability of an IPC system with functional or performance dependence was also investigated [10]. An extended object-oriented Petri net model was proposed for mission reliability simulation of a repairable phased-mission system with both external and internal CCFs [11]. Using the reliability block diagram method, the quantitative study of probability of failing safely (PFS) for a safety instrumented system showed that the CCF increases the PFS [12]. Under the assumption that a single-failure event may lead to simultaneous failure of multiple components, a nonparametric predictive inference for system reliability following the CCF of components was presented [13].In many cases, the system and/or its components can function in some different states characterized by several performance levels. Such systems are often referred to as MSSs. The theory of the MSS was investigated by Murchland in 1975 [14]. Many researchers have analyzed the reliability which is often seen as an important measure of the MSS to provide the desired performance level [15]. MSS can also be subjected to the CCF which can lead to the failure of an entire system or subsystem [16].Many studies have investigated the reliability and modeling of the MSS with CCF. For instance, an algorithm was proposed to evaluate the reliability of a complex nonrepairable series-parallel MSS with CCF under the assumption that the failure propagation time is a random value with a given distribution [17]. An MSS reliability analysis method incorporated a CCF into a Bayesian network [18] which is based on a well-defined theory of probabilistic reasoning and the ability to express the complex dependence between random variables. To reduce the computational complexity caused by the MSS, where the number of states will increase rapidly with the number of elements, a universal generating function (UGF) was adopted by Ushakov [19]. Further developments and applications have been studied in detail [20]. Recently, the reliability analysis method based on the UGF has been employed in many environments. A new reliability-evaluation methodology based on the UGF and a recursive algorithm was applied to a multistate weighted k-out-of-n system [21]. To analyze the reliability of the generalized linear multistate consecutively connected system, a UGF-based method was suggested [22]. Using UGF, the evaluation method of reliability characteristics for a nonrepair complex system has been presented [23].Although much interest has been focused on methods of reliability analysis for systems with CCF, less attention has been paid to the MSS with CCF, which only occurs in several independent elements. The reliability of a transmission MSS considering the CCF effects has been analyzed based on Bayesian network model [18]. This model can clearly express the influence of CCF on system reliability without the computation of minimum cut sets or the determination of algebraic expression of unreliability. However, these parameter values used in β-factor model reduce its accuracy because these values are based on engineering experiences or published statistics of CCF. An implicit two-stage procedure is proposed to evaluate reliability of MSS with CCF [3]. This procedure can obtain the system reliability function in both analytic and numerical solution. But the procedure can be applied only to series-parallel MSS in which each component can belong only to one common cause group.Aiming at the research limitations mentioned above, such as the inaccuracies of parameter values and nonoverlap of common cause group, an easily programmed recursive method is suggested for this situation, with the benefit of reduced computational complexity for reliability and performance distribution assessment. Using the sensitivity analysis method, it is easy to locate the optimal solution to eliminate the CCF and will help engineers to find the bottlenecks of system reliability design. This method is based on the application of UGF.This paper is organized as follows. Section2 formulates the MSS with CCF. In Section 3, based on the expression of the CCF, the algorithm to calculate the system reliability is presented. Several illustrative examples are shown in Section 4. Section 5 discusses the conclusions.


## 2. System Description

### 2.1. Formulations of the MSS

Under the framework of multistate modeling, an MSS consists ofn components connected in series-parallel. A characteristic of the MSS is that any component i can be in one of mi+1 performance levels corresponding to states s∈{0,1,2,…,mi},i=1,2,…,n. State 0 is for the perfect function, state mi for complete failure and other intermediate states are for partial failure or degradation. Without loss of generality, the performance level at each state of component i will take values from the set(1)hi=hi0,hi1,…,himi.Then the space of all combinations of performance levels for all components can be given by(2)Hn=h1×h2×⋯×hn,where the operator × denotes the Cartesian product. The event of every performance level for component i at any time is also assumed to a discrete random variable Hi∈hi. The probability set of Hi can be expressed by(3)pi=pi0,pi1,…,pimi.Obviously∑0mipis=1, because that the mi+1 states constitute a group of mutually exclusive events, i.e., the component i can be in one and only one state from s at any time. The relationship between the variable Hi and performance gi can be given by(4)pis=Pr⁡Hi=his.Another characteristic of the MSS is that it usually will consist of more than two states. Here, some external factors causing the failures or degradations of a system are not considered. Under this assumption, the performance levels of the MSS will be unambiguously ascertained by performance levels of its components, whose states will consequently determine the states of the MSS.Suppose that the MSS hasM+1 states and the variable wj is the performance level corresponding to system state j∈{0,1,2,…,M}. The MSS performance level can be seen as a discrete random variable W taking values from the set L={w0,w1,…,wM}. Using (2), the system function can be written as(5)W=ϕH1,H2,…,Hn:Hn→L.The above formula is the map of the performance level space from all components to the system. Similarly, the probability of system performance can be expressed as(6)qj=PrW=wj,j=0,1,…,M.From the analysis above, the keys of the MSS model are given as follows:(i)
The MSS consists ofn components connected in series-parallel.(ii)
Any componenti has mi+1 states. Each state’s performance level and probability of a corresponding performance level can be expressed as (1) and (3), respectively.(iii)
The MSS will haveM+1 states deduced from the states of the n components as defined in (5) and its performance level will take a value from L={w0,w1,…,wM}.(iv)
The probability of the MSS in a corresponding state can be expressed by (6).


### 2.2. Indices of Reliability

According to the origin of the CCF, there are two types of cause: external and internal cause. The CCF caused by the latter is often called a propagated failure because it will affect other system components. In an MSS where the CCF has an internal cause, the failed components will propagate to other components which may be mutually independent or even overlap with one or more components. Then the system performance level will be undermined with an induced decrease in reliability. When the system performance level is reduced to a limit valueθ, which is often called system demand, the system will be seen as unacceptable. The reliability of the MSS with CCF can be defined as the probability that the system satisfies the value of θ. From (6), one can obtain(7)Rθ=∑j=0MqjΔwj,θ,where(8)Δwj,θ=1,wj≥θΔwj,θ=0,otherwise.Another important measure is the conditional expected performance Ω(θ). It expresses the system’s expected performance under the condition that the MSS is in an acceptable state. Having the system reliability R(θ), this measure can be calculated by(9)Ωθ=∑j=0MqjwjΔwj,θRθ.To calculate these measures, the performance level distribution of a system should be first obtained according to (6). UGF has been proved an effective method for the reliability assessment of different types of MSS. In particular, the series-parallel system is likely to adopt it by the recursive method.


## 2.1. Formulations of the MSS

Under the framework of multistate modeling, an MSS consists ofn components connected in series-parallel. A characteristic of the MSS is that any component i can be in one of mi+1 performance levels corresponding to states s∈{0,1,2,…,mi},i=1,2,…,n. State 0 is for the perfect function, state mi for complete failure and other intermediate states are for partial failure or degradation. Without loss of generality, the performance level at each state of component i will take values from the set(1)hi=hi0,hi1,…,himi.Then the space of all combinations of performance levels for all components can be given by(2)Hn=h1×h2×⋯×hn,where the operator × denotes the Cartesian product. The event of every performance level for component i at any time is also assumed to a discrete random variable Hi∈hi. The probability set of Hi can be expressed by(3)pi=pi0,pi1,…,pimi.Obviously∑0mipis=1, because that the mi+1 states constitute a group of mutually exclusive events, i.e., the component i can be in one and only one state from s at any time. The relationship between the variable Hi and performance gi can be given by(4)pis=Pr⁡Hi=his.Another characteristic of the MSS is that it usually will consist of more than two states. Here, some external factors causing the failures or degradations of a system are not considered. Under this assumption, the performance levels of the MSS will be unambiguously ascertained by performance levels of its components, whose states will consequently determine the states of the MSS.Suppose that the MSS hasM+1 states and the variable wj is the performance level corresponding to system state j∈{0,1,2,…,M}. The MSS performance level can be seen as a discrete random variable W taking values from the set L={w0,w1,…,wM}. Using (2), the system function can be written as(5)W=ϕH1,H2,…,Hn:Hn→L.The above formula is the map of the performance level space from all components to the system. Similarly, the probability of system performance can be expressed as(6)qj=PrW=wj,j=0,1,…,M.From the analysis above, the keys of the MSS model are given as follows:(i)
The MSS consists ofn components connected in series-parallel.(ii)
Any componenti has mi+1 states. Each state’s performance level and probability of a corresponding performance level can be expressed as (1) and (3), respectively.(iii)
The MSS will haveM+1 states deduced from the states of the n components as defined in (5) and its performance level will take a value from L={w0,w1,…,wM}.(iv)
The probability of the MSS in a corresponding state can be expressed by (6).


## 2.2. Indices of Reliability

According to the origin of the CCF, there are two types of cause: external and internal cause. The CCF caused by the latter is often called a propagated failure because it will affect other system components. In an MSS where the CCF has an internal cause, the failed components will propagate to other components which may be mutually independent or even overlap with one or more components. Then the system performance level will be undermined with an induced decrease in reliability. When the system performance level is reduced to a limit valueθ, which is often called system demand, the system will be seen as unacceptable. The reliability of the MSS with CCF can be defined as the probability that the system satisfies the value of θ. From (6), one can obtain(7)Rθ=∑j=0MqjΔwj,θ,where(8)Δwj,θ=1,wj≥θΔwj,θ=0,otherwise.Another important measure is the conditional expected performance Ω(θ). It expresses the system’s expected performance under the condition that the MSS is in an acceptable state. Having the system reliability R(θ), this measure can be calculated by(9)Ωθ=∑j=0MqjwjΔwj,θRθ.To calculate these measures, the performance level distribution of a system should be first obtained according to (6). UGF has been proved an effective method for the reliability assessment of different types of MSS. In particular, the series-parallel system is likely to adopt it by the recursive method.


## 3. Analysis Method

The CCF caused by one component corresponds to the state0′ with the performance level hi0′=c. The performance value can coincide with component performance in the state of local failure and can usually be seen as zero. When the component i is in state 0′, all the system components affected by this CCF will also be in the failure state with performance c. If one component cannot arouse a CCF, the probability mass function (PMF) of the corresponding state should be zeroed: pi0′=0. The UGF of component i incorporating the CCF can be rewritten as(10)uiz=∑ki=0mipiki∗zhiki+pi0′∗zc.The conditional PMF of any componenti that cannot fail due to the CCF can be represented by(11)ui′z=∑ki=0mipiki1-pi0′∗zhiki,where pi0′ is the probability of causing the CCF.Suppose there areb≤n components that can independently and simultaneously cause the CCF. For component vector β={β(1),β(2),…,β(r),…,β(b)}, r∈[1,b], the probability that the CCF originated from one component is expressed as pβr0′. Since these CCFs can be combined independently, the number of combination is 2b.For any combinationδ,(0≤δ≤2b-1), the CCF originates from component β(r) if(12)εr=mod2δ2r-1=1,1≤r≤b.After evaluating the above formula for δ=0 to δ=2b-1, the set of components corresponding to the combination δ can be obtained. For every combination δ, the probability of the corresponding CCF is(13)ϕδ=∏r=1bpβr0′εr1-pβr0′1-εr,0≤δ≤2b-1.According to the properties of the CCF, other components failed because componentβ(r) can be denoted by Sβ(r). For every combination δ, the set of components affected by the CCF and aroused the CCF can be obtained by(14)μδ=⋃r=1bSβr′,where (15)Sβr′=Sβr∪βr,εr=1⌀,εr=0.Given that one component fails due to the CCF, the conditional PMF of its performance level can be represented byψ(z)=zc, which indicates the component can be only in the failure state with performance level c.When the CCF corresponding to combinationδ occurs, all components in the set μδ will come in the CCF mode, and their UGFs will be replaced by ψ(z). However, the UGFs of components not belonging to the set μδ must be represented by (11). Then the conditional PMF of the overall system performance will be obtained as uδ(z) in the form of a UGF. From (13), the system UGF can be calculated as(16)Usysz=∑δ=02b-1uδz∗ϕδ.Based on the UGF of the whole system, the relevant indices of reliability can be assessed by combining (7) and (9).Combined with the above analysis, the following steps should be adopted to realize the assessment:S1. Incorporate the CCF into each component’s UGF according to (10).S2. Obtain the conditional UGF ui′(z) of component i by (11).S3. Determine all possible combinations of components aroused by the CCF.S4. For every combination δ,(0≤δ≤2b-1), fix its corresponding components which can inspire the CCF according to (12).S5. Calculate the probability of every combination δ using (13).S6. Determine the set of components affected by the CCF given one combination δ according to (14).S7. Replace the UGF of all components belonging to set μδ with ψ(z), and use the UGF ui′(z) of the components not belonging to it.S8. Based on the physical structure of the series-parallel MSS, represent the conditional UGF uδ(z).S9. Obtain the UGF of the entire system Usys(z), according to (16).AfterUsys(z) has been obtained, the reliability indices can be further calculated according to (7) and (9). The following section will illustrate the above method.


## 4. Application Examples

### 4.1. Basic Application

A numerical example illustrates the above method for the reliability assessment of the MSS with CCF. A type of MSS with a series-parallel structure, as shown in Figure1, is a typical configuration. In a flow transition system, matter such as oil, water, or steam is often transmitted from terminal A to terminal B by interconnected components.Figure 1
Series-parallel structure MSS.For example, in the feeding water system of a power plant, components (pumps) are grouped into two subsystems; sub1 consists of C11, C12, and C13 connected in parallel, and sub2 consists of C21 and C22. Their parameters of performance distribution are listed in Table1. As shown in this table, all components can fail with a given probability. But only components C11 and C13 can arouse the CCF. C11 can destroy component C21 owing to its nearby location and then constitute the CCF group CG1 with component C21. Component C13 will stimulate the failure of components C21 and C22 at the same time, to form the second CCF group CG2. The reasons for the CCF can be events such as fire or lightning.Table 1
Performance distribution of components.
Component Perfect state Local failure Failure to cause CCF p h p p C11 0.88 3 0.07 0.05 C12 0.89 5 0.11 - C13 0.79 3 0.15 0.06 C21 0.86 6 0.14 - C22 0.85 2 0.15 -According to the above steps, the UGF of every component can be obtained as follows:(17)uC11z=0.88∗z3+0.07∗z0+0.05∗zcuC12z=0.89∗z3+0.11∗z0+0.0∗zcuC13z=0.79∗z5+0.15∗z0+0.06∗zcuC21z=0.86∗z6+0.14∗z0+0.0∗zcuC22z=0.85∗z2+0.15∗z0+0.0∗zc.From (11), the conditional UGFs of components take the following forms:(18)uC11′z=0.9263∗z3+0.0737∗z0uC12′z=0.89∗z3+0.11∗z0uC13′z=0.8404∗z5+0.1596∗z0uC21′z=0.86∗z6+0.14∗z0uC22′z=0.85∗z2+0.15∗z0.The conditional UGFs of component C12, C21, and C22 can be obtained by removing 0.0∗zc directly. Two components, β(1)=C11,β(2)=C13, cause the CCF, and then Sβ(1)={C21}, Sβ(2)={C21,C22}. The other variables are listed in Table 2.Table 2
Value ofε(r), ϕδ, and μδ for any combination δ.
δ ε ( r ) ϕ δ μ δ r = 1 r = 2 0 0 0 =(1-0.05)∗(1-0.06)=0.893 ⌀ 1 1 0 =0.05∗(1-0.06)=0.047 {C11,C21} 2 0 1 =(1-0.05)∗0.06=0.057 { C 13 , C 21 , C 22 } 3 1 1 =0.05∗0.06=0.003 { C 11 , C 21 } ∪ { C 13 , C 21 , C 22 }According to the physical nature of the components’ interconnection, the function within subsystems sub1 and sub2 should take the form of the sum, and the function between two subsystems should take the form of the minimum. Then the conditional UGF of the entire system can be denoted as(19)uδz=usub1z⊗minusub2z=uC11′z⊗sumuC12′z⊗sumuC13′z⊗minuC21′z⊗sumuC22′z.Forδ=0 and μδ=⌀, no UGF should be replaced with ψ(z). Equation (19) can be expressed as(20)u0z=0.9263∗z3+0.0737∗z0⊗sum0.89∗z3+0.11∗z0⊗sum0.8404∗z5+0.1596∗z0⊗min0.86∗z6+0.14∗z0⊗sum0.85∗z2+0.15∗z0=0.6094∗z8+0.2207∗z6+0.0059∗z5+0.0230∗z3+0.1188∗z2+0.0222∗z0.Forδ=1 and μδ={C11,C21}, only uC11′(z) and uC21′(z) will be substituted by ψ(z)=zc. In the case of failure, the performance level is zero. Equation (19) will be rewritten as(21)u1z=z0⊗sum0.89∗z3+0.11∗z0⊗sum0.8404∗z5+0.1596∗z0⊗minz0⊗sum0.85∗z2+0.15∗z0=0.835∗z2+0.165∗z0.Forδ=2 and μδ={C13,C21,C22}, the UGF of components C13,C21,C22 must be replaced by ψ(z) in (19). So one can obtain the following formula:(22)u2z=0.9263∗z3+0.0737∗z0⊗sum0.89∗z3+0.11∗z0⊗sumz0⊗minz0⊗sumz0=z0.Whenδ=3, the UGF of components C11,C13,C21,C22 should be substituted with ψ(z). Then (19) will be written as(23)u3z=z0⊗sum0.89∗z3+0.11∗z0⊗sumz0⊗minz0⊗sumz0=z0.The whole system’s UGF can be obtained by(24)Usysz=∑δ=03uδz∗ϕδ=0.5442∗z8+0.1971∗z6+0.0053∗z5+0.0205∗z3+0.1453∗z2+0.0876∗z0.Now the reliability of the entire system with CCF can be depicted as in Figure2.Figure 2
Reliability of the system.According to (7), the reliability of the system with performance θ=4 can be calculated as follows. Point A in Figure 2 is the reflection of the reliability value:(25)R4=∑j=08qjΔwj,4=0.5442+0.1971+0.0053=0.7466.Using (9), the conditional expected system performance level at reliability R(4) is(26)Ω4=8∗0.5442+6∗0.1971+5∗0.00530.7466=7.4507.


### 4.2. Advanced Application

A more practical example will be depicted and the sensitivity of the CCF can be further analyzed. Suppose an MSS production system consists of nine components connected in series-parallel, as shown in Figure3.Figure 3
An MSS production system.The system can be divided into the subsystems. Components C22, C31, and C33 have only binary states. All the others can be seen as multistate components with two operating states: perfect and degraded performance.The CCF can be inspired by four componentsβ={C11,C12,C32,C34}. The set of components affected by the CCF are Sβ(1)={C21}, Sβ(2)={C22,C23}, Sβ(3)={C21}, and Sβ(4)={C23,C32}. These components constitute four CCF groups, CG1, CG2, CG3, and CG4. The performance distributions of system components for operating states and failing states are listed in Table 3.Table 3
Performance distribution of components.
Component Perfection Degradation Local failure Failure of causing CCF p h p h p p C11 0.8 9 0.10 7 0.05 0.05 C12 0.75 12 0.09 8 0.05 0.11 C21 0.85 8 0.12 6 0.03 - C22 0.90 11 - - 0.10 - C23 0.82 7 0.13 5 0.05 - C31 0.93 13 - - 0.07 - C32 0.82 11 0.05 7 0.11 0.02 C33 0.89 9 - - 0.11 - C34 0.82 11 0.10 7 0.05 0.03According to the analysis method suggested above, the reliability curve at different performance levels can be represented as in Figure4.Figure 4
Curve of reliability at different performance levels.There are 21 performance levels corresponding to this MSS, as shown in the results. For example, the data points A, B, C, and D are the reliabilities at performance:R(5)=0.9782, R(10)=0.8154, R(15)=0.7228, and R(20)=0.4711, respectively.Under the limitation of cost, it is necessary to find the most valuable component than can improve the reliability of the system to the greatest extent when the CCF originating from it is eliminated. In other words, the sensitivity of components can be further analyzed by the suggested method. When the CCF of componenti is localized, its UGF will be rewritten as follows, according to (10)(27)uiz=∑ki=1mipiki∗zhiki+pi0+pi0′∗z0.The CCF can be incorporated in the failure state without arousing other components’ failures, and the performance level is also zeroed. Then the system reliability can be reassessed when the CCF in one single component is eliminated. The improvements shown in Table4 have been compared with the case where no CCF is removed.Table 4
Comparison of reliability improvement after the CCF is removed from one single component.
θ R(θ) Improvement after specified component’s CCF is eliminated β ( 1 ) = C 11 β ( 2 ) = C 12 β ( 3 ) = C 32 β ( 4 ) = C 34 8 0.943819 0.42% 1.50% 0.33% 0.05% 9 0.856588 0.38% 10.14% 0.17% 0.28% 10 0.815374 0.40% 0.00% 0.17% 0.43% 11 0.815374 0.40% 0.00% 0.17% 0.43% 12 0.808013 0.70% 0.00% 0.30% 0.63% 13 0.740249 0.00% 0.00% 0.31% 0.54% 14 0.722797 0.00% 0.00% 0.27% 0.72% 15 0.722768 0.00% 0.00% 0.27% 0.72%In Table4, the performance level θ is chosen from 8 to 15 and the reliability R(θ) corresponds to cases where no CCF is excluded. For the four components that can cause the CCF, the improved reliability is expressed in percentage form after the corresponding CCF has been eliminated.Moreover, for a specified performance levelθ, the best component can be chosen to localize its CCF. For θ=9, the component β(2)=C12 should be the best one to be improved because its improving proportion can reach 10.14%. At the same time, when one component has been chosen, the best performance level can be determined according to the extent of improvement. For example, because of some limitations only β(1)=C11 can be maintained to remove its CCF, and the best performance level corresponding to the biggest reliability improvement, 0.7%, is θ=12. The analysis of sensitivity to CCF localization will help engineers to find bottlenecks in the system reliability design and to properly direct investment in reliability enhancement.


## 4.1. Basic Application

A numerical example illustrates the above method for the reliability assessment of the MSS with CCF. A type of MSS with a series-parallel structure, as shown in Figure1, is a typical configuration. In a flow transition system, matter such as oil, water, or steam is often transmitted from terminal A to terminal B by interconnected components.Figure 1
Series-parallel structure MSS.For example, in the feeding water system of a power plant, components (pumps) are grouped into two subsystems; sub1 consists of C11, C12, and C13 connected in parallel, and sub2 consists of C21 and C22. Their parameters of performance distribution are listed in Table1. As shown in this table, all components can fail with a given probability. But only components C11 and C13 can arouse the CCF. C11 can destroy component C21 owing to its nearby location and then constitute the CCF group CG1 with component C21. Component C13 will stimulate the failure of components C21 and C22 at the same time, to form the second CCF group CG2. The reasons for the CCF can be events such as fire or lightning.Table 1
Performance distribution of components.
Component Perfect state Local failure Failure to cause CCF p h p p C11 0.88 3 0.07 0.05 C12 0.89 5 0.11 - C13 0.79 3 0.15 0.06 C21 0.86 6 0.14 - C22 0.85 2 0.15 -According to the above steps, the UGF of every component can be obtained as follows:(17)uC11z=0.88∗z3+0.07∗z0+0.05∗zcuC12z=0.89∗z3+0.11∗z0+0.0∗zcuC13z=0.79∗z5+0.15∗z0+0.06∗zcuC21z=0.86∗z6+0.14∗z0+0.0∗zcuC22z=0.85∗z2+0.15∗z0+0.0∗zc.From (11), the conditional UGFs of components take the following forms:(18)uC11′z=0.9263∗z3+0.0737∗z0uC12′z=0.89∗z3+0.11∗z0uC13′z=0.8404∗z5+0.1596∗z0uC21′z=0.86∗z6+0.14∗z0uC22′z=0.85∗z2+0.15∗z0.The conditional UGFs of component C12, C21, and C22 can be obtained by removing 0.0∗zc directly. Two components, β(1)=C11,β(2)=C13, cause the CCF, and then Sβ(1)={C21}, Sβ(2)={C21,C22}. The other variables are listed in Table 2.Table 2
Value ofε(r), ϕδ, and μδ for any combination δ.
δ ε ( r ) ϕ δ μ δ r = 1 r = 2 0 0 0 =(1-0.05)∗(1-0.06)=0.893 ⌀ 1 1 0 =0.05∗(1-0.06)=0.047 {C11,C21} 2 0 1 =(1-0.05)∗0.06=0.057 { C 13 , C 21 , C 22 } 3 1 1 =0.05∗0.06=0.003 { C 11 , C 21 } ∪ { C 13 , C 21 , C 22 }According to the physical nature of the components’ interconnection, the function within subsystems sub1 and sub2 should take the form of the sum, and the function between two subsystems should take the form of the minimum. Then the conditional UGF of the entire system can be denoted as(19)uδz=usub1z⊗minusub2z=uC11′z⊗sumuC12′z⊗sumuC13′z⊗minuC21′z⊗sumuC22′z.Forδ=0 and μδ=⌀, no UGF should be replaced with ψ(z). Equation (19) can be expressed as(20)u0z=0.9263∗z3+0.0737∗z0⊗sum0.89∗z3+0.11∗z0⊗sum0.8404∗z5+0.1596∗z0⊗min0.86∗z6+0.14∗z0⊗sum0.85∗z2+0.15∗z0=0.6094∗z8+0.2207∗z6+0.0059∗z5+0.0230∗z3+0.1188∗z2+0.0222∗z0.Forδ=1 and μδ={C11,C21}, only uC11′(z) and uC21′(z) will be substituted by ψ(z)=zc. In the case of failure, the performance level is zero. Equation (19) will be rewritten as(21)u1z=z0⊗sum0.89∗z3+0.11∗z0⊗sum0.8404∗z5+0.1596∗z0⊗minz0⊗sum0.85∗z2+0.15∗z0=0.835∗z2+0.165∗z0.Forδ=2 and μδ={C13,C21,C22}, the UGF of components C13,C21,C22 must be replaced by ψ(z) in (19). So one can obtain the following formula:(22)u2z=0.9263∗z3+0.0737∗z0⊗sum0.89∗z3+0.11∗z0⊗sumz0⊗minz0⊗sumz0=z0.Whenδ=3, the UGF of components C11,C13,C21,C22 should be substituted with ψ(z). Then (19) will be written as(23)u3z=z0⊗sum0.89∗z3+0.11∗z0⊗sumz0⊗minz0⊗sumz0=z0.The whole system’s UGF can be obtained by(24)Usysz=∑δ=03uδz∗ϕδ=0.5442∗z8+0.1971∗z6+0.0053∗z5+0.0205∗z3+0.1453∗z2+0.0876∗z0.Now the reliability of the entire system with CCF can be depicted as in Figure2.Figure 2
Reliability of the system.According to (7), the reliability of the system with performance θ=4 can be calculated as follows. Point A in Figure 2 is the reflection of the reliability value:(25)R4=∑j=08qjΔwj,4=0.5442+0.1971+0.0053=0.7466.Using (9), the conditional expected system performance level at reliability R(4) is(26)Ω4=8∗0.5442+6∗0.1971+5∗0.00530.7466=7.4507.


## 4.2. Advanced Application

A more practical example will be depicted and the sensitivity of the CCF can be further analyzed. Suppose an MSS production system consists of nine components connected in series-parallel, as shown in Figure3.Figure 3
An MSS production system.The system can be divided into the subsystems. Components C22, C31, and C33 have only binary states. All the others can be seen as multistate components with two operating states: perfect and degraded performance.The CCF can be inspired by four componentsβ={C11,C12,C32,C34}. The set of components affected by the CCF are Sβ(1)={C21}, Sβ(2)={C22,C23}, Sβ(3)={C21}, and Sβ(4)={C23,C32}. These components constitute four CCF groups, CG1, CG2, CG3, and CG4. The performance distributions of system components for operating states and failing states are listed in Table 3.Table 3
Performance distribution of components.
Component Perfection Degradation Local failure Failure of causing CCF p h p h p p C11 0.8 9 0.10 7 0.05 0.05 C12 0.75 12 0.09 8 0.05 0.11 C21 0.85 8 0.12 6 0.03 - C22 0.90 11 - - 0.10 - C23 0.82 7 0.13 5 0.05 - C31 0.93 13 - - 0.07 - C32 0.82 11 0.05 7 0.11 0.02 C33 0.89 9 - - 0.11 - C34 0.82 11 0.10 7 0.05 0.03According to the analysis method suggested above, the reliability curve at different performance levels can be represented as in Figure4.Figure 4
Curve of reliability at different performance levels.There are 21 performance levels corresponding to this MSS, as shown in the results. For example, the data points A, B, C, and D are the reliabilities at performance:R(5)=0.9782, R(10)=0.8154, R(15)=0.7228, and R(20)=0.4711, respectively.Under the limitation of cost, it is necessary to find the most valuable component than can improve the reliability of the system to the greatest extent when the CCF originating from it is eliminated. In other words, the sensitivity of components can be further analyzed by the suggested method. When the CCF of componenti is localized, its UGF will be rewritten as follows, according to (10)(27)uiz=∑ki=1mipiki∗zhiki+pi0+pi0′∗z0.The CCF can be incorporated in the failure state without arousing other components’ failures, and the performance level is also zeroed. Then the system reliability can be reassessed when the CCF in one single component is eliminated. The improvements shown in Table4 have been compared with the case where no CCF is removed.Table 4
Comparison of reliability improvement after the CCF is removed from one single component.
θ R(θ) Improvement after specified component’s CCF is eliminated β ( 1 ) = C 11 β ( 2 ) = C 12 β ( 3 ) = C 32 β ( 4 ) = C 34 8 0.943819 0.42% 1.50% 0.33% 0.05% 9 0.856588 0.38% 10.14% 0.17% 0.28% 10 0.815374 0.40% 0.00% 0.17% 0.43% 11 0.815374 0.40% 0.00% 0.17% 0.43% 12 0.808013 0.70% 0.00% 0.30% 0.63% 13 0.740249 0.00% 0.00% 0.31% 0.54% 14 0.722797 0.00% 0.00% 0.27% 0.72% 15 0.722768 0.00% 0.00% 0.27% 0.72%In Table4, the performance level θ is chosen from 8 to 15 and the reliability R(θ) corresponds to cases where no CCF is excluded. For the four components that can cause the CCF, the improved reliability is expressed in percentage form after the corresponding CCF has been eliminated.Moreover, for a specified performance levelθ, the best component can be chosen to localize its CCF. For θ=9, the component β(2)=C12 should be the best one to be improved because its improving proportion can reach 10.14%. At the same time, when one component has been chosen, the best performance level can be determined according to the extent of improvement. For example, because of some limitations only β(1)=C11 can be maintained to remove its CCF, and the best performance level corresponding to the biggest reliability improvement, 0.7%, is θ=12. The analysis of sensitivity to CCF localization will help engineers to find bottlenecks in the system reliability design and to properly direct investment in reliability enhancement.


## 5. Concluding Remarks

In this paper a reliability assessment method is suggested for the series-parallel MSS with CCF. The components that can arouse the CCF may cause failures of different subsets of system components. These subsets can even overlap on some components. This method is based on the UGF application, and the factor of the CCF is incorporated in the expression of the UGF. Because of its repetition for every combination, the method is conveniently programmed by iteration. The indices of reliability and expected conditional performance level can be obtained. The sensitivity of components can be further analyzed. This method will help for reliability engineers to decide which component is the most valuable to be invested in, so as to eliminate its CCF. According to the assumptions in this method, the structure of the MSS is applied only to the series-parallel system. The more complex MSS topologies such as bridge and G:(k/n) structures need to be considered in the future.



---

*Source: 6535726-2019-05-28.xml*