Same Example Modeled with RBDs or Fault Trees

This example also appears in the System analysis reference. The following example demonstrates how you can model the same analysis scenario using either RBDs or fault trees in BlockSim. The results will be the same with either approach. This discussion presents the RBD and fault tree solutions together so you can compare the methods. As an alternative, you could also review RBD Model and Fault Tree Model, which present the steps for each modeling approach separately.

Problem Statement
Assume that a component can fail due to six independent primary failure modes: A,  B,  C,  D,  E  and  F. Some of these primary modes can be broken down further into the events that can cause them, or sub-modes. Furthermore, assume that once a mode occurs, the event also occurs and the mode does not go away. Specifically:
 * The component fails if mode A,  B  or  C  occurs.
 * If mode D,  E  or  F  occurs alone, the component does not fail; however, the component will fail if any two (or more) of these modes occur (i.e.,  D  and  E ;  D  and  F ;  E  and  F).
 * Modes D,  E  and  F  have a constant rate of occurrence (exponential distribution) with mean times of occurrence of  200,000, 175,000  and 500,000 hours, respectively.
 * The rates of occurrence for modes A,  B  and  C  depend on their sub-modes.

Do the following:


 * 1) Determine the reliability of the component after 1 year (8,760 hours).
 * 2) Determine the B10 life of the component.
 * 3) Determine the mean time to failure (MTTF) of the component.
 * 4) Rank the modes in order of importance at 1 year.
 * 5) Recalculate results 1, 2 and 3 assuming mode B  is eliminated.

To begin the analysis, modes A,  B  and  C  can be broken down further based on specific events (sub-modes), as defined next.

Mode A
There are five independent events (sub-modes) associated with mode A : events  S1,  S2,  T1,  T2  and  Y. It is assumed that events  S1  and  S2  each have a constant rate of occurrence with a probability of occurrence in a single year (8,760 hours) of  1  in  10,000  and  1  in  20,000, respectively. Events T1  and  T2  are more likely to occur in an older component than a newer one (i.e., they have an increasing rate of occurrence) and have a probability of occurrence of  1  in  10,000  and  1  in  20,000, respectively, in a single year and  1  in  1,000  and  1  in  3,000, respectively, after two years. Event Y  also has a constant rate of occurrence with a probability of occurrence of  1  in  1,000  in a single year. There are three possible ways for mode A  to manifest itself:


 * 1) Events S1  and  S2  both occur.
 * 2) Event T1  or  T2  occurs.
 * 3) Event Y  and either event  S1  or event  S2  occur (i.e., events  Y  and  S1  or events  Y  and  S2  occur).

RBD Solution for Mode A

The RBD that satisfies the conditions for mode A  is shown in the figure below.



Each mode is identified in the RBD. Furthermore, two additional items are included: a starting block (NF) and an end node (2/2). The starting block and the end node are set so they cannot fail and, therefore, will not affect the results. The end node is used to define a 2-out-of-2 configuration (i.e., both paths leading into the node must work).

Fault Tree Solution for Mode A

The fault tree for mode A  is shown in the figure below.



Each mode is identified as an event in the fault tree. The following figure shows an alternative representation for mode A  using mirrored events for  S1  and  S2.



Mode A  Discussion

The system reliability equation for this configuration (regardless of how it is drawn) is:


 * $$\begin{align}

R(t)= & -2{{R}_{T2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{S1}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{S2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{T1}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{Y}} \\ & +{{R}_{T2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{S1}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{S2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{T1}} \\ & +{{R}_{T2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{S1}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{T1}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{Y}} \\ & +{{R}_{T2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{S2}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{T1}}\text{ }\!\!\cdot\!\!\text{ }{{R}_{Y}} \end{align}\,\!$$

Based on the given probabilities, distribution parameters are computed for each block (either RBD block or the fault tree event block). One way is to compute them using the Parameter Experimenter, as shown in the figure below. In this figure and for S1, the probability is  1  in  10,000  in one year (8,760  hours), thus the exponential failure rate is 1.1416e-8. This can be repeated for S2  and  Y.



Events T1  and  T2  need to be modeled using a life distribution that does not have a constant failure rate. Using BlockSim's Parameter Experimenter and selecting the Weibull distribution, the parameter values for events T1  and  T2  are shown in the figures below.



Mode B
There are three dependent events associated with mode B : events  BA,  BB  and  BC.


 * Two out of the three events must occur for mode B  to occur.
 * o	Events BA,  BB  and  BC  all have an exponential distribution with a mean of  50,000  hours.
 * o	The events are dependent (i.e., if BA,  BB  or  BC  occurs, then the remaining events are more likely to occur). Specifically, when one event occurs, the  MTTF  of the remaining events is halved.

This is basically a load sharing configuration. The reliability function for each block will change depending on the other events. Therefore, the reliability of each block is not only dependent on time, but also on the stress (load) that the block experiences.

RBD Solution for Mode B

The reliability block diagram for mode B  is shown in the figure below.

Fault Tree Solution for Mode B

The fault tree for mode B  is shown in the figure below. A Load Sharing gate is used.



Mode B  Discussion

To describe the dependency, a Load Sharing gate and dependent event blocks are used. Since the failure rate is assumed to be constant, an exponential distribution is used. Furthermore, for simplicity, an Arrhenius life-stress relationship is used with the parameters B=2.0794 and C=6250.

Mode C
There are two sequential events associated with mode C :  CA  and  CB.


 * Both events must occur for mode C  to occur.
 * Event CB  will only occur if event  CA  has occurred.
 * If event CA  has not occurred, then event  CB  will not occur.
 * Events CA  and  CB  both occur based on a Weibull distribution.
 * For event CA, $$\beta \,\!$$ = 2 and $$\eta \,\!$$ = 30,000 hours.
 * For event CB, $$\beta \,\!$$ = 2 and $$\eta \,\!$$ = 10,000 hours.

RBD Solution for Mode C

To model this, you can think of a scenario similar to standby redundancy. Basically, if CA  occurs then  CB  gets initiated. A Standby container can be used to model this, as shown in the figure below.



In this case, event CA  is set as the active component and  CB  as the standby. If event CA  occurs,  CB  will be initiated. For this analysis, a perfect switch is assumed. The properties are set in BlockSim as follows:

Contained Items


 * CA : Active failure distribution, Weibull distribution ($$\beta \,\!$$ = 2, $$\eta \,\!$$ = 30,000).
 * CA : Quiescent failure distribution: None, cannot fail or age in this mode.
 * CB : Active failure distribution, Weibull distribution ($$\beta \,\!$$ = 2, $$\eta \,\!$$ = 10,000).
 * CB : Quiescent failure distribution: None, cannot fail or age in this mode.

Switch


 * Active Switching: Always works (100% reliability) and instant switch (no delays).
 * Quiescent Switch failure distribution: None, cannot fail or age in this mode.

Fault Tree Solution for Mode C

The fault tree for mode C  is shown in the figure below. Note that the sequence is enforced by the Standby gate (used as a Sequence Enforcing gate).



Mode C  Discussion

The failure distribution settings for event CA  are shown in the figure below.



The failure distribution properties for event CB  are set in the same manner.

Modes D, E and F
Modes D,  E  and  F  can all be represented using the exponential distribution. The failure distribution properties for modes D,  E  and  F  are:


 * D : MTTF  = 200,000 hours.
 * E : MTTF  = 175,000 hours.
 * F : MTTF  = 500,000 hours.

The Entire Component
The last step is to set up the model for the component based on the primary modes (A, B,  C,  D,  E  and  F). Modes A,  B  and  C  can each be represented by single blocks that encapsulate the subdiagrams already created. The RBD in the first figure below represents the primary failure modes for the component while the fault tree in second figure below illustrates the same. The node represented by 2/3 in the RBD indicates a 2-out-of-3 configuration. The Voting OR gate in the fault tree accomplishes the same. Subdiagrams are used in both configurations for the sub-modes.





Once the diagrams have been created, the reliability equation for the system can be obtained, as follows:


 * $$\begin{align}

R{{(t)}_{System}}= & R{{(t)}_{A}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{B}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{F}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{D}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{C}} \\ & +R{{(t)}_{A}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{B}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{F}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{C}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{E}} \\ & +R{{(t)}_{A}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{B}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{D}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{C}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{E}} \\ & -2(R{{(t)}_{A}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{B}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{F}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{D}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{C}}\text{ }\!\!\cdot\!\!\text{ }R{{(t)}_{E}}) \end{align}\,\!$$

where $$R{{(t)}_{A}}\,\!$$, $$R{{(t)}_{B}}\,\!$$ and $$R{{(t)}_{C}}\,\!$$ are the reliability equations corresponding to the sub-modes.

Analysis and Discussion
The questions posed earlier can be answered using BlockSim. Regardless of the approach used (i.e., RBD or FTA), the answers are the same.


 * 1.	The reliability of the component at 1 year (8,760 hours) can be calculated using the Analytical Quick Calculation Pad (QCP) or by viewing the reliability vs. time plot, as displayed in the following figure. $$R(t=8760)=86.4975%\,\!$$.


 * 2.	Using the Analytical QCP, the B10 life of the component is estimated to be 7,373.94 hours.
 * 3.	Using the Analytical QCP, the mean life of the component is estimated to be 21,659.68 hours.
 * 4.	The ranking of the modes after 1 year can be shown via the Static Reliability Importance plot, as shown in the figure below.


 * 5.	Re-computing the results for 1, 2 and 3 assuming mode B  is removed:
 * a)	 R(t=8760) =98.72%.
 * b)	 B10 = 16,928.38 hours.
 * c)	 MTTF = 34,552.89 hours.

There are multiple options for modeling systems with fault trees and RBDs in BlockSim. The first figure below shows the complete fault tree for the component without using subdiagrams (Transfer gates) while the second figure below illustrates a hybrid analysis utilizing an RBD for the component and fault trees as the subdiagrams. The results are the same regardless of the option chosen.