Time-Dependent System Reliability (Analytical): Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
m (→‎Conditional Reliability: changed R(T,t) to R(t|T))
 
(461 intermediate revisions by 9 users not shown)
Line 1: Line 1:
{{Template:bsbook|5}}
{{Template:bsbook|4}}
In the [[RBDs and Analytical System Reliability]] chapter, different system configuration types were examined, as well as different methods for obtaining the system's reliability function analytically.  Because the reliabilities in the problems presented were treated as probabilities (e.g., <math>P(A)\,\!</math>, <math>{{R}_{i}}\,\!</math> ), the reliability values and equations presented were referred to as ''static'' (not time-dependent). Thus, in the prior chapter, the life distributions of the components were not incorporated in the process of calculating the system reliability. In this chapter, time dependency in the reliability function will be introduced. We will develop the models necessary to observe the reliability over the life of the system, instead of at just one point in time. In addition, performance measures such as failure rate, MTTF and warranty time will be estimated for the entire system. The methods of obtaining the reliability function analytically remain identical to the ones presented in the previous chapter, with the exception that the reliabilities will be functions of time.  In other words, instead of dealing with <math>{{R}_{i}}\,\!</math>, we will use <math>{{R}_{i}}(t)\,\!</math>.  All examples in this chapter assume that no repairs are performed on the components. Repairable systems analysis will be introduced in a [[Introduction_to_Repairable_Systems|subsequent chapter]]. 


==Analytical Life Predictions==
The analytical approach presented in the prior chapter involved the determination of a mathematical expression that describes the reliability of the system, expressed in terms of the reliabilities of its components.  So far we have estimated only static system reliability (at a fixed time).  For example, in the case of a system with three components in series, the system's reliability equation was given by:
<br>
::<math>{{R}_{s}}={{R}_{1}}\cdot {{R}_{2}}\cdot {{R}_{3}}  \ \,\!</math>


In the previous chapter, different system configuration types were examined, as well as different methods for obtaining the system's reliability function analytically.  Because the reliabilities in the problems presented were treated as probabilities (e.g.  <math>P(A)</math> ,  <math>{{R}_{i}}</math> ), the reliability values and equations presented were referred to as static (not time-dependent).  Thus, in the prior chapter, the life distributions of the components were not incorporated in the process of calculating the system reliability.  In this chapter, time dependency in the reliability function will be introduced.  We will develop the models necessary to observe the reliability over the life of the system, instead of at just one point in time.  In addition, performance measures such as failure rate, MTTF and warranty time will be estimated for the entire system.  The methods of obtaining the reliability function analytically remain identical to the ones presented in the previous chapter, with the exception that the reliabilities will be functions of time.  In other words, instead of dealing with  <math>{{R}_{i}}</math> , we will use  <math>{{R}_{i}}(t)</math> .  All examples in this chapter assume that no repairs are performed on the components.   
The values of <math>{{R}_{1}}\,\!</math>, <math>{{R}_{2}}\,\!</math> and <math>{{R}_{3}}\,\!</math> were given for a common time and the reliability of the system was estimated for that timeHowever, since the component failure characteristics can be described by distributions, the system reliability is actually time-dependentIn this case, the equation above can be rewritten as:


=Analytical Life Predictions=
::<math>{{R}_{s}}(t)={{R}_{1}}(t)\cdot {{R}_{2}}(t)\cdot {{R}_{3}}(t)\,\!</math>
The analytical approach presented in the prior chapter involved the determination of a mathematical expression that describes the reliability of the system, expressed in terms of the reliabilities of its components.  So far we have estimated only static system reliability (at a fixed time).  For example, in the case of a system with three components in series, the system's reliability equation was given by:
<br>


::<math>{{R}_{s}}={{R}_{1}}\cdot {{R}_{2}}\cdot {{R}_{3}}</math>
The reliability of the system for any mission time can now be estimated.  Assuming a Weibull life distribution for each component, the first equation above can now be expressed in terms of each component's reliability function, or:
<br>


The values of  <math>{{R}_{1}}</math> ,  <math>{{R}_{2}}</math>  and  <math>{{R}_{3}}</math>  were given for a common time and the reliability of the system was estimated for that time.  However, since the component failure characteristics can be described by distributions, the system reliability is actually time-dependent.  In this case, Eqn. (1) can be rewritten as:
<br>


::<math>{{R}_{s}}(t)={{R}_{1}}(t)\cdot {{R}_{2}}(t)\cdot {{R}_{3}}(t)</math>
::<math>{{R}_{s}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{2}}} \right)}^{{{\beta }_{2}}}}}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{3}}} \right)}^{{{\beta }_{3}}}}}}\,\!</math>


<br>
In the same manner, any life distribution can be substituted into the system reliability equation.  Suppose that the times-to-failure of the first component are described with a Weibull distribution, the times-to-failure of the second component with an exponential distribution and the times-to-failure of the third component with a normal distribution.  Then the first equation above can be written as:


The reliability of the system for any mission time can now be estimated.  Assuming a Weibull life distribution for each component, Eqn. (1) can now be expressed in terms of each component's reliability function, or:
::<math>{{R}_{s}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}}\cdot {{e}^{-{{\lambda }_{2}}t}}\cdot \left[ 1-\Phi \left( \frac{t-{{\mu }_{3}}}{{{\sigma }_{3}}} \right) \right]\,\!</math>


<br>
It can be seen that the biggest challenge is in obtaining the system's reliability function in terms of component reliabilities, which has already been discussed in depth.  Once this has been achieved, calculating the reliability of the system for any mission duration is just a matter of substituting the corresponding component reliability functions into the system reliability equation.


::<math>{{R}_{s}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{2}}} \right)}^{{{\beta }_{2}}}}}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{3}}} \right)}^{{{\beta }_{3}}}}}}</math>


<br>
===Advantages and Disadvantages===
The primary advantage of the analytical solution is that it produces a mathematical expression that describes the reliability of the system.  Once the system's reliability function has been determined, other calculations can then be performed to obtain metrics of interest for the system. Such calculations include: 


In the same manner, any life distribution can be substituted into the system reliability equation. Suppose that the times-to-failure of the first component are described with a Weibull distribution, the times-to-failure of the second component with an exponential distribution and the times-to-failure of the third component with a normal distribution.  Then Eqn. (1) can be written as:
:*Determination of the system's ''pdf''.
<br>


:*Determination of warranty periods.


::<math>{{R}_{s}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}}\cdot {{e}^{-{{\lambda }_{2}}t}}\cdot \left[ 1-\Phi \left( \frac{t-{{\mu }_{3}}}{{{\sigma }_{3}}} \right) \right]</math>
:*Determination of the system's failure rate.


<br>
:*Determination of the system's MTTF.


It can be seen that the biggest challenge is in obtaining the system's reliability function in terms of component reliabilities, which has already been discussed in depth.  Once this has been achieved, calculating the reliability of the system for any mission duration is just a matter of substituting the corresponding component reliability functions into the system reliability equation.
===Advantages of the Analytical Method===
The primary advantage of the analytical solution is that it produces a mathematical expression that describes the reliability of the system.  Once the system's reliability function has been determined, other calculations can then be performed to obtain metrics of interest for the system. Such calculations include: 
<br>
:• Determination of the system's  <math>pdf.</math><br>
:• Determination of warranty periods.<br>
:• Determination of the system's failure rate.<br>
:• Determination of the system's MTTF.<br>
<br>
In addition, optimization and reliability allocation techniques can be used to aid engineers in their design improvement efforts.  Another advantage in using analytical techniques is the ability to perform static calculations and analyze systems with a mixture of static and time-dependent components.  Finally, the reliability importance of components over time can be calculated with this methodology.
In addition, optimization and reliability allocation techniques can be used to aid engineers in their design improvement efforts.  Another advantage in using analytical techniques is the ability to perform static calculations and analyze systems with a mixture of static and time-dependent components.  Finally, the reliability importance of components over time can be calculated with this methodology.
===Disadvantages of the Analytical Method===
 
The biggest disadvantage of the analytical method is that formulations can become very complicated.  The more complicated a system is, the larger and more difficult it will be to analytically formulate an expression for the system's reliability.  For particularly detailed systems this process can be quite time-consuming, even with the use of computers.  Furthermore, when the maintainability of the system or some of its components must be taken into consideration, analytical solutions become intractable.  In these situations, the use of simulation methods may be more advantageous than attempting to develop a solution analytically.  Simulation methods are presented in later chapters.
The biggest disadvantage of the analytical method is that formulations can become very complicated.  The more complicated a system is, the larger and more difficult it will be to analytically formulate an expression for the system's reliability.  For particularly detailed systems this process can be quite time-consuming, even with the use of computers.  Furthermore, when the maintainability of the system or some of its components must be taken into consideration, analytical solutions become intractable.  In these situations, the use of simulation methods may be more advantageous than attempting to develop a solution analytically.  Simulation methods are presented in later chapters.
===Looking at a Simple ``Complex'' System Analytically===


The complexity involved in an analytical solution can be best illustrated by looking at the simple ``complex'' system with 15 components, as shown in Figure Ch5fig1.
==Looking at a Simple "Complex" System Analytically==
<br>
The complexity involved in an analytical solution can be best illustrated by looking at the simple ''complex'' system with 15 components, as shown below.
[[File:5-1.png|center]]
 
<center>Figure 5.1: An RBD of a complex system.</center>
[[Image:5-1.png|center|600px|<div align="center"> An RBD of a complex system.</div>|link=]]
<br>
 
The system reliability for this system (computed using BlockSim) is shown next.  The first solution is provided using BlockSim's symbolic solution.  In symbolic mode, BlockSim breaks the equation into segments, identified by tokens, that need to be substituted into the final system equation for a complete solution.  This creates algebraic solutions that are more compact than if the substitutions were made.
The system reliability for this system (computed using BlockSim) is shown next.  The first solution is provided using BlockSim's symbolic solution.  In symbolic mode, BlockSim breaks the equation into segments, identified by tokens, that need to be substituted into the final system equation for a complete solution.  This creates algebraic solutions that are more compact than if the substitutions were made.
<br>


::<math>\begin{align}
::<math>\begin{align}
Line 67: Line 57:
D1 = & +{{R}_{D}}\cdot ID \\  
D1 = & +{{R}_{D}}\cdot ID \\  
ID = & -{{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}}   
ID = & -{{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}}   
\end{align}</math>
\end{align}\,\!</math>
<br>
 
Substituting the terms yields:  
Substituting the terms yields:  
<br>


::<math>\begin{align}
::<math>\begin{align}
Line 80: Line 69:
& -{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}} \\  
& -{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}} \\  
& +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}}\}   
& +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}}\}   
\end{align}</math>
\end{align}\,\!</math>


<br>
BlockSim's automatic algebraic simplification would yield the following format for the above solution:  
BlockSim's automatic algebraic simplification would yield the following format for the above solution:  
<br>


::<math>\begin{align}
::<math>\begin{align}
Line 94: Line 81:
& -{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\  
& -{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\  
& +RI\cdot {{R}_{O}}\cdot {{R}_{F}} \\  
& +RI\cdot {{R}_{O}}\cdot {{R}_{F}} \\  
& +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}})){{R}_{L}})   
& +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}})){{R}_{L}})  \
\end{align}</math>
\end{align}\,\!</math>
 
In this equation, each <math>{{R}_{i}}\,\!</math> represents the reliability function of a block.  For example, if <math>{{R}_{A}}\,\!</math> has a Weibull distribution, then each <math>{{R}_{A}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{A}}} \right)}^{{{\beta }_{A}}}}}}\,\!</math> and so forth.  Substitution of each component's reliability function in the last <math>{{R}_{System}}\,\!</math> equation above will result in an analytical expression for the system reliability as a function of time, or <math>{{R}_{s}}(t)\,\!</math>, which is the same as <math>(1-cd{{f}_{System}}).\,\!</math>


<br>
==Obtaining Other Functions of Interest==
In this equation, each  <math>{{R}_{i}}</math>  represents the reliability function of a block.  For example, if  <math>{{R}_{A}}</math>  has a Weibull distribution, then each  <math>{{R}_{A}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{A}}} \right)}^{{{\beta }_{A}}}}}}</math>  and so forth.  Substitution of each component's reliability function in Eqn. (rseq) will result in an analytical expression for the system reliability as a function of time, or <math>{{R}_{s}}(t)</math> , which is the same as  <math>(1-cd{{f}_{System}}).</math>
Once the system reliability equation (or the cumulative density function, ''cdf'') has been determined, other functions and metrics of interest can be derived.   


===Obtaining Other Functions of Interest===
Once the system reliability equation (or the cumulative density function,  <math>cdf</math> ) has been determined, other functions and metrics of interest can be derived. 
<br>
Consider the following simple system:
Consider the following simple system:
<br>
[[File:5-2.png|center]]
<br>
<center>Figure 5.2: Simple two-component system.</center>
<br>
Furthermore, assume that component 1 follows an exponential distribution with a mean of 10,000 (<math>\mu =10,000,</math>  <math>\lambda =1/10,000)</math>  and component 2 follows a Weibull distribution with  <math>\beta =6</math>  and  <math>\eta =10,000</math> .  The reliability equation of this system is:


<br>
[[Image:5-2.png|center|200px|<div align="center"> Simple two-component system. </div>|link=]]
 
Furthermore, assume that component 1 follows an exponential distribution with a mean of 10,000 (<math>\mu =10,000,\,\!</math> <math>\lambda =1/10,000)\,\!</math> and component 2 follows a Weibull distribution with <math>\beta =6\,\!</math> and <math>\eta =10,000\,\!</math>.  The reliability equation of this system is:
 
::<math>\begin{align}
::<math>\begin{align}
{{R}_{S}}(t)= & {{R}_{1}}(t)\cdot {{R}_{2}}(t) \\  
{{R}_{S}}(t)= & {{R}_{1}}(t)\cdot {{R}_{2}}(t) \\  
= & {{e}^{-\lambda t}}\cdot {{e}^{-{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \\  
= & {{e}^{-\lambda t}}\cdot {{e}^{-{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \\  
= & {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}   
= & {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \  
\end{align}</math>
\end{align}\,\!</math>


<br>
The system ''cdf'' is:
The system <math>cdf</math>  is:


<br>
::<math>\begin{align}
::<math>\begin{align}
{{F}_{S}}(t)= & 1-({{R}_{1}}(t)\cdot {{R}_{2}}(t)) \\  
{{F}_{S}}(t)= & 1-({{R}_{1}}(t)\cdot {{R}_{2}}(t)) \\  
= & 1-\left( {{e}^{-\lambda t}}\cdot {{e}^{-{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \right) \\  
= & 1-\left( {{e}^{-\lambda t}}\cdot {{e}^{-{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \right) \\  
= & 1-\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)   
= & 1-\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)   
\end{align}</math>
\end{align}\,\!</math>
 
<br>
====System  <math>pdf</math>====
Once the equation for the reliability of the system has been obtained, the system's  <math>pdf</math>  can be determined. The  <math>pdf</math>  is the derivative of the reliability function with respect to time or:
<br>


::<math>{{f}_{s}}(t)=-\frac{d({{R}_{s}}(t))}{dt}</math>
===System ''pdf''===
<br>
Once the equation for the reliability of the system has been obtained, the system's ''pdf'' can be determined. The  ''pdf'' is the derivative of the reliability function with respect to time or:


<br>
::<math>{{f}_{s}}(t)=-\frac{d({{R}_{s}}(t))}{dt} \ \,\!</math>
For the system shown in Figure Ch5fig2, this is:


<br>
For the system shown above, this is:


::<math>\begin{align}
::<math>\begin{align}
Line 146: Line 121:
& +{{e}^{-\tfrac{1}{10,000}t}}\left[ -\frac{d}{dt}\left( {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right) \right] \\  
& +{{e}^{-\tfrac{1}{10,000}t}}\left[ -\frac{d}{dt}\left( {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right) \right] \\  
= & {{f}_{1}}(t)\cdot {{R}_{2}}(t)+{{f}_{2}}(t)\cdot {{R}_{1}}(t)   
= & {{f}_{1}}(t)\cdot {{R}_{2}}(t)+{{f}_{2}}(t)\cdot {{R}_{1}}(t)   
\end{align}</math>
\end{align}\,\!</math>
<br>
 
The next figure shows a plot of the ''pdf'' equation.


Figure Ch5fig3 shows a plot of Eqn. (fig2pdf).
[[Image:BS5.3.png|center|650px|pdf plot of the two-component system|link=]]
<br>


===Conditional Reliability===
====Conditional Reliability====
Conditional reliability is the probability of a system successfully completing another mission following the successful completion of a previous mission.  The time of the previous mission and the time for the mission to be undertaken must be taken into account for conditional reliability calculations.  The system's conditional reliability function is given by:
Conditional reliability is the probability of a system successfully completing another mission following the successful completion of a previous mission.  The time of the previous mission and the time for the mission to be undertaken must be taken into account for conditional reliability calculations.  The system's conditional reliability function is given by:


::<math>R(t|T)=\frac{R(T+t)}{R(T)}\,\!</math>


::<math>R(T,t)=\frac{R(T+t)}{R(T)}</math>
Equation above gives the reliability for a new mission of duration <math>t\,\!</math> having already accumulated <math>T\,\!</math> hours of operation up to the start of this new mission. The system is evaluated to assure that it will start the next mission successfully.
 


Eqn. (conditional) gives the reliability for a new mission of duration  <math>t</math>  having already accumulated  <math>T</math>  hours of operation up to the start of this new mission. The system is evaluated to assure that it will start the next mission successfully.
For the simple two-component system, the reliability for mission of <math>t\,\!</math> = 1000 hours, having an age of <math>T\,\!</math> = 500 hours, is:
<br>
For the system in Figure Ch5fig2, the reliability for mission of <math>t=1,000</math> hours, having an age of <math>T=500</math> hours, is:
<br>


::<math>\begin{align}
::<math>\begin{align}
{{R}_{S}}(T=500,t=1000)= & \frac{R(T+t)}{R(T)} \\  
{{R}_{S}}(t=1000|T=500)= & \frac{R(T+t)}{R(T)} \\  
= & \frac{R(1500)}{R(500)} \\  
= & \frac{R(1500)}{R(500)} \\  
= & \frac{{{e}^{-\tfrac{1500}{10,000}}}\cdot {{e}^{-{{\left( \tfrac{1500}{10,000} \right)}^{6}}}}}{{{e}^{-\tfrac{500}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{500}{10,000} \right)}^{6}}}}} \\  
= & \frac{{{e}^{-\tfrac{1500}{10,000}}}\cdot {{e}^{-{{\left( \tfrac{1500}{10,000} \right)}^{6}}}}}{{{e}^{-\tfrac{500}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{500}{10,000} \right)}^{6}}}}} \\  
= & 0.9048=90.48%   
= & 0.9048=90.48%   
\end{align}</math>
\end{align}\,\!</math>
[[File:BS5.3.png|center]]
<br>
<br>
[[File:BS5.4.png|center]]
<br>


====Conditional Reliability for Components====
===Conditional Reliability for Components===
Now in this formulation, it was assumed that the accumulated age was equivalent for both units. That is, both started life at zero and aged to 500.  It is possible to consider an individual component that has already accumulated some age (used component) in the same formulation.  To illustrate this, assume that component 2 started life with an age of T = 100.  Then the reliability equation of the system, as given in <math>{{R}_{S}}(t)= {{R}_{1}}(t)\cdot {{R}_{2}}(t)\,\!</math>, would need to be modified to include a conditional term for 2, or:


Now in this formulation, it was assumed that the accumulated age was equivalent for both units. That is, both started life at zero and aged to 500.  It is possible to consider an individual component that has already accumulated some age (used component) in the same formulation.  To illustrate this, assume that component 2 started life with an age of 100.  Then the reliability equation of the system, as given in Eqn. (reqgsyst), would need to be modified to include a conditional term for 2, or:
::<math>{{R}_{S}}(t)={{R}_{1}}(t)\cdot \frac{{{R}_{2}}({{T}_{2}}+t)}{{{R}_{2}}({{T}_{2}})} \ \,\!</math>
<br>
 
::<math>{{R}_{S}}(t)={{R}_{1}}(t)\cdot \frac{{{R}_{2}}({{T}_{2}}+t)}{{{R}_{2}}({{T}_{2}})}</math>
 
<br>


In BlockSim, the start age input box may be used to specify a starting age greater than zero.
In BlockSim, the start age input box may be used to specify a starting age greater than zero.
====System Failure Rate====
Once the distribution of the system has been determined, the failure rate can also be obtained by dividing the  <math>pdf</math>  by the reliability function:


<br>
===System Failure Rate===
 
Once the distribution of the system has been determined, the failure rate can also be obtained by dividing the ''pdf'' by the reliability function:
::<math>{{\lambda }_{s}}\left( t \right)=\frac{{{f}_{s}}\left( t \right)}{{{R}_{s}}\left( t \right)}</math>
<br>


::<math>{{\lambda }_{s}}\left( t \right)=\frac{{{f}_{s}}\left( t \right)}{{{R}_{s}}\left( t \right)}  \ \,\!</math>


For the system in Figure Ch5fig2:  
For the simple two-component system:  
<br>


::<math>\begin{align}
::<math>\begin{align}
Line 203: Line 161:
= & \frac{-\tfrac{d}{dt}\left( {{e}^{-\tfrac{1}{10,000}t}} \right)}{{{e}^{-\tfrac{1}{10,000}t}}}+\frac{-\tfrac{d}{dt}\left( {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)}{{{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}} \\  
= & \frac{-\tfrac{d}{dt}\left( {{e}^{-\tfrac{1}{10,000}t}} \right)}{{{e}^{-\tfrac{1}{10,000}t}}}+\frac{-\tfrac{d}{dt}\left( {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)}{{{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}} \\  
= & \frac{{{f}_{1}}}{{{R}_{1}}}+\frac{{{f}_{2}}}{{{R}_{2}}} \\  
= & \frac{{{f}_{1}}}{{{R}_{1}}}+\frac{{{f}_{2}}}{{{R}_{2}}} \\  
= & {{\lambda }_{1}}+{{\lambda }_{2}}
= & {{\lambda }_{1}}+{{\lambda }_{2}}   \
\end{align}</math>
\end{align}\,\!</math>
<br>
 
The following figure shows a plot of the equation.


::Figure 5.4 shows a plot of Eqn. (9).
[[Image:BS5.4.png|center|650px|Failure rate function plot of the two component system|link=]]
<br>


BlockSim uses numerical methods to estimate the failure rate.  It should be pointed out that as <math>t\to \infty </math> , numerical evaluation of Eqn. (system failure rate) is constrained by machine numerical precision. That is, there are limits as to how large <math>t</math> can get before floating point problems arise.  For example, at <math>t=5,000,000</math> both numerator and denominator will tend to zero (e.g. <math>{{e}^{-\tfrac{5,000,000}{10,000}}}=7.1245\times {{10}^{-218}}</math> ).  As these numbers become very small they will start looking like a zero to the computer, or cause a floating point error, resulting in a <math>\tfrac{0}{0}</math> or <math>\tfrac{X}{0}</math> operation.  In these cases, BlockSim will return a value of " <math>N/A</math> " for the result.  Obviously, this does not create any practical constraints.
BlockSim uses numerical methods to estimate the failure rate.  It should be pointed out that as <math>t\to \infty \,\!</math>, numerical evaluation of the first equation above is constrained by machine numerical precision. That is, there are limits as to how large <math>t\,\!</math> can get before floating point problems arise.  For example, at <math>t=5,000,000\,\!</math> both numerator and denominator will tend to zero (e.g., <math>{{e}^{-\tfrac{5,000,000}{10,000}}}=7.1245\times {{10}^{-218}}\,\!</math> ).  As these numbers become very small they will start looking like a zero to the computer, or cause a floating point error, resulting in a <math>\tfrac{0}{0}\,\!</math> or <math>\tfrac{X}{0}\,\!</math> operation.  In these cases, BlockSim will return a value of "<math>N/A\,\!</math>" for the result.  Obviously, this does not create any practical constraints.


====System Mean Life (Mean Time To Failure)====
===System Mean Life (Mean Time To Failure)===
The mean life (or mean time to failure, MTTF) can be obtained by integrating the system reliability function from zero to infinity:  
The mean life (or mean time to failure, MTTF) can be obtained by integrating the system reliability function from zero to infinity:  
<br>


::<math>MTTF=\mathop{}_{0}^{\infty }{{R}_{s}}\left( t \right)dt</math>
::<math>MTTF=\int_{0}^{\infty }{{R}_{s}}\left( t \right)dt   \ \,\!</math>


<br>
The mean time is a performance index and does not provide any information about the behavior of the failure distribution of the system.


The mean time is a performance index and does not provide any information about the behavior of the failure distribution of the system.
For the simple two-component system:  
<br>
For the system in Figure Ch5fig2:  
<br>


::<math>\begin{align}
::<math>\begin{align}
MTTF= & \mathop{}_{0}^{\infty }\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)dt \\  
MTTF= & \int_{0}^{\infty }\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)dt \\  
= & 5978.9   
= & 5978.9   
\end{align}</math>
\end{align}\,\!</math>


<br>
===Warranty Period and BX Life===
====Warranty Period and BX Life====
Sometimes it is desirable to know the time value associated with a certain reliability.  Warranty periods are often calculated by determining what percentage of the failure population can be covered financially and estimating the time at which this portion of the population will fail.  Similarly, engineering specifications may call for a certain BX life, which also represents a time period during which a certain proportion of the population will fail.  For example, the B10 life is the time in which 10% of the population will fail.   
Sometimes it is desirable to know the time value associated with a certain reliability.  Warranty periods are often calculated by determining what percentage of the failure population can be covered financially and estimating the time at which this portion of the population will fail.  Similarly, engineering specifications may call for a certain BX life, which also represents a time period during which a certain proportion of the population will fail.  For example, the B10 life is the time in which 10% of the population will fail.   
This is obtained by setting <math>{{R}_{S}}(t)</math>  to the desired value and solving for  <math>t.</math>
This is obtained by setting <math>{{R}_{S}}(t)\,\!</math> to the desired value and solving for <math>t.\,\!</math>  
<br>
For the system in Figure Ch5fig2:
<br>
 
::<math>{{R}_{s}}\left( t \right)={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}</math>
<br>
 
To compute the time by which reliability would be equal to 90%, Eqn. (warranty) is recast as follows and solved for  <math>t.</math>
 
<br>
::<math>0.90={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}</math>
 
<br>
 
In this case,  <math>t=1053.59</math> .  Equivalently, the B10 life for this system is also  <math>1053.59</math> .
Except for some trivial cases, a closed form solution for  <math>t</math>  cannot be obtained.  Thus, it is necessary to solve for  <math>t</math>  using numerical methods.  BlockSim uses numerical methods.
 
===Example 1===
Consider a system consisting of three exponential units in series with the following failure rates (in failures per hour):  <math>{{\lambda }_{1}}</math>  = 0.0002,  <math>{{\lambda }_{2}}</math>  = 0.0005 and <math>{{\lambda }_{3}}</math>  = 0.0001.
<br>
:• Obtain the reliability equation for the system.<br>
:• What is the reliability of the system after 150 hours of operation?<br>
:• Obtain the system's  <math>pdf.</math> <br>
:• Obtain the system's failure rate equation.<br>
:• What is the MTTF for the system?<br>
:• What should the warranty period be for a 90% reliability?<br>
====Solution to Example 1====
<br>
:• The analytical expression for the reliability of the system is given by:
 
 
::<math>\begin{align}
{{R}_{s}}(t)= & {{R}_{1}}(t)\cdot {{R}_{2}}(t)\cdot {{R}_{3}}(t) \\
= & {{e}^{-{{\lambda }_{1}}t}}\cdot {{e}^{-{{\lambda }_{2}}t}}\cdot {{e}^{-{{\lambda }_{1}}t}} \\
= & {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}} 
\end{align}</math>
 
<br>
 
At 150 hours of operation, the reliability of the system is:
<br>
 
 
::<math>\begin{align}
{{R}_{s}}(t)= & {{e}^{-(0.0002+0.0005+0.0001)150}} \\
= & 0.8869\text{ or }88.69% 
\end{align}</math>
 
<br>
 
:• In order to obtain the system's  <math>pdf</math> , the derivative of the reliability equation given in Eqn. (exp units reliability) is taken with respect to time, or:
<br>
::<math>\begin{align}
{{f}_{s}}(t)= & -\frac{d[{{R}_{s}}(t)]}{dt} \\
= & -\frac{d\left[ {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}} \right]}{dt} \\
= & ({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})\cdot {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}} 
\end{align}</math>
:• The system's failure rate can now be obtained simply by dividing the system's  <math>pdf</math>  given in Eqn. (exp units pdf) by the system's reliability function given in Eqn. (exp units reliability), or:
 
 
::<math>\begin{align}
{{\lambda }_{s}}\left( t \right)= & \frac{{{f}_{s}}\left( t \right)}{{{R}_{s}}\left( t \right)} \\
= & \frac{({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})\cdot {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}}}{{{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}}} \\
= & ({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}}) \\
= & 0.0008\text{ }fr/hr 
\end{align}</math>
 
<br>
 
:• Combining Eqn. (system MTTF) and Eqn. (exp units reliability), the system's MTTF can be obtained:
 
<br>
 
<math>\begin{align}
MTTF= & \mathop{}_{0}^{\infty }{{R}_{s}}\left( t \right)dt \\
= & \mathop{}_{0}^{\infty }{{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}}dt \\
= & \frac{1}{({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})} \\
= & 1250\text{ }hr 
\end{align}</math>
<br>
 
 
:• Solving Eqn. (exp units reliability) with respect to time will yield the corresponding warranty period for a 90% reliability.  In this case, the system reliability equation is simple and a closed form solution exists.  The warranty time can now be found by solving:
 
<br>
 
::<math>\begin{align}
t= & -\frac{\ln (R)}{{{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}}} \\
= & -\frac{\ln (0.9)}{0.0008} \\
= & 131.7\text{ }hr 
\end{align}</math>
<br>
 
 
Thus, the warranty period should be 132 hours.
===Example 2===
Consider the system shown in Figure fig4a.
<br>
[[File:BS5.5.png|center]]
<br>
<center>Figure 5.5: Complex bridge system in Example 2.</center>
<br>
 
Components  <math>A</math>  through  <math>E</math>  are Weibull distributed with  <math>\beta =1.2</math>  and  <math>\eta =1230</math>  hours.  The starting and ending blocks cannot fail.  <br>
Determine the following:
<br>
:• The reliability equation for the system and its corresponding plot.<br>
:• The system's  <math>pdf</math>  and its corresponding plot.<br>
:• The system's failure rate equation and the corresponding plot.<br>
:• The  <math>MTTF</math> .<br>
:• The warranty time for a 90% reliability.<br>
:• The reliability for a 200-hour mission, if it is known that the system has already successfully operated for 200 hours.<br>
====Solution====
<br>
The first step is to obtain the reliability function for the system.  The methods described in the previous chapter can be employed, such as the event space or path-tracing methods.  Using BlockSim, the following reliability equation is obtained:
<br>
 
::<math>\begin{align}
{{R}_{s}}(t)= & ({{R}_{Start}}\cdot {{R}_{End}}(2{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\
& -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}-{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\
& -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{B}}\cdot {{R}_{E}}-{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\
& -{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}}+{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\
& +{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}+{{R}_{A}}\cdot {{R}_{D}}+{{R}_{B}}\cdot {{R}_{E}})) 
\end{align}</math>
<br>
 
Note that since the starting and ending blocks cannot fail,  <math>{{R}_{Start}}=1</math>  and  <math>{{R}_{End}}=1,</math>  Eqn. (brequation) can be reduced to:
 
<br>
 
::<math>\begin{align}
{{R}_{s}}(t)= & 2\cdot {{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\
& -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}-{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\
& -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{B}}\cdot {{R}_{E}}-{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\
& -{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}}+{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\
& +{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}+{{R}_{A}}\cdot {{R}_{D}}+{{R}_{B}}\cdot {{R}_{E}} 
\end{align}</math>
 
<br>
 
Where  <math>{{R}_{A}}</math>  is the reliability equation for Component A, or:
 
<br>
 
::<math>{{R}_{A}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{A}}} \right)}^{{{\beta }_{A}}}}}}</math>
 
<br>
 
::<math>{{R}_{B}}</math>  is the reliability equation for Component  <math>B</math> , etc.
<br>


Since the components in this example are identical, the system reliability equation can be further reduced to:
For the simple two-component system:  


<br>
::<math>{{R}_{s}}\left( t \right)={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}\,\!</math>


::<math>{{R}_{s}}(t)=2R{{(t)}^{2}}+2R{{(t)}^{3}}-5R{{(t)}^{4}}+2R{{(t)}^{5}}</math>
To compute the time by which reliability would be equal to 90%, equation above is recast as follows and solved for <math>t.\,\!</math>  


<br>
::<math>0.90={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}\,\!</math>


Or, in terms of the failure distribution:
In this case, <math>t=1053.59\,\!</math>.  Equivalently, the B10 life for this system is also 1053.59.
Except for some trivial cases, a closed form solution for <math>t\,\!</math> cannot be obtained.  Thus, it is necessary to solve for <math>t\,\!</math> using numerical methods.  BlockSim uses numerical methods.


<br>


::<math>{{R}_{s}}(t)=2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}</math>
===Examples===<!-- THIS SECTION HEADER IS LINKED TO THE PAGES THAT ARE TRANSCLUDED IN THIS SECTION. IF YOU RENAME THE SECTION, YOU MUST UPDATE THE LINK. -->
'''Components in Series'''<hr>
{{:Time-Dependent System Reliability for Components in Series}}


<br>


The corresponding plot is given in Figure 5-6.
'''Components in Parallel'''<hr>
<br>
{{:Time-Dependent System Reliability for Components in Parallel}}
[[File:BS5.6.png|center]]
<br>
<center>Figure 5.6: Reliability plot for the system in Figure 5.5.</center>
<br>
<br>
In order to obtain the system's  <math>pdf</math> , the derivative of the reliability equation given in Eqn. (system eqn 3) is taken with respect to time, resulting in:
<br>


::<math>\begin{align}
==Approximating the System cdf==<!-- THIS SECTION HEADER IS LINKED TO: Example Using a Distribution to Approximate the CDF. IF YOU RENAME THE SECTION, YOU MUST UPDATE THE LINK. -->
{{f}_{s}}(t)= & 4\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+6\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \\
In many cases, it is valuable to fit a distribution that represents the system's times-to-failure.  This can be useful when the system is part of a larger assembly and may be used for repeated calculations or in calculations for other systems.  In cases such as this, it can be useful to characterize the system's behavior by fitting a distribution to the overall system and calculating parameters for this distribution.  This is equivalent to fitting a single distribution to describe <math>{{R}_{S}}(t)\,\!</math>.  In essence, it is like reducing the entire system to a component in order to simplify calculations. 
& -20\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+10\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} 
\end{align}</math>


<br>
For the system shown below:
[[Image:5-2.png|center|250px|<div align="center">Simple two-component system. </div>|link=]]


The  <math>pdf</math>  can now be plotted for different time values,  <math>t</math> , as shown in Figure Ch5fig6.
::<math>{{R}_{S}}(t)={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}\,\!</math>
<br>
The system's failure rate can now be obtained by dividing the system's  <math>pdf</math>  given in Eqn. (bridge pdf) by the system's reliability function given in Eqn. (system eqn 3), or:
 
<br>
 
::<math>\begin{align}
{{\lambda }_{s}}(t)= & \frac{4\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+6\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}}{2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}} \\
& +\frac{-20\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+10\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}}{2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}} 
\end{align}</math>
 
<br>
 
The corresponding plot is given in Figure 5.7.
<br>
[[File:BS5.7.gif|center]]
<br>
<center>Figure 5.7: ''pdf'' plot for the system in Figure 5.5.</center>
<br>
The  <math>MTTF</math>  of the system is obtained by integrating the system's reliability function given by Eqn. (system eqn 3) from time zero to infinity, as given by Eqn. (system MTTF).  Using BlockSim's Analytical QCP, an  <math>MTTF</math>  of 1007.8 hours is calculated, as shown in Figure Ch5fig8.
<br>
 
The warranty time can be obtained by solving Eqn. (system eqn 3) with respect to time for a system reliability  <math>{{R}_{s}}=0.9</math> .  Using the Analytical QCP and selecting the <br>
Warranty Time option, a time of 372.72 hours is obtained, as shown in Figure Ch5fig9.
<br>
[[File:BS5.8.gif|center]]
<br>
<center>Figure 5.8: Failure rate for the system in Figure 5.5.</center>
<br>
[[File:BS5.9.png|center]]
<br>
<center>Figure 5.9: MTTF of the system in Figure 5.5.</center>
<br>
[[File:BS5.10.png|center]]
<br>
<center>Figure 5.10: Time at which ''R''=0.9 or 90% for the system in Figure 5.5.</center>
<br>
[[File:BS5.11.png|center]]
<br>
<center>Figure 5.11: Conditional reliability calculation for the system in Figure 5.5.</center>
<br>
 
Lastly, the conditional reliability can be obtained using Eqn. (conditional) and Eqn. (system eqn 3), or:
<br>
 
::<math>\begin{align}
R(200,200)= & \frac{R(400)}{R(200)} \\
= & \frac{0.883825}{0.975321} \\
= & 0.906189 
\end{align}</math>
 
<br>
 
This can be calculated using BlockSim's Analytical QCP, as shown in Figure Ch5fig10.
<br>
 
=Approximating the System CDF=
 
In many cases, it is valuable to fit a distribution that represents the system's times-to-failure.  This can be useful when the system is part of a larger assembly and may be used for repeated calculations or in calculations for other systems.  In cases such as this, it can be useful to characterize the system's behavior by fitting a distribution to the overall system and calculating parameters for this distribution.  This is equivalent to fitting a single distribution to describe  <math>{{R}_{S}}(t</math> ).  In essence, it is like reducing the entire system to a component in order to simplify calculations. 
<br>
For the system in Figure Ch5fig2:
<br>
 
::<math>{{R}_{S}}(t)={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}</math>
 
<br>
To compute an approximate reliability function for this system, <math>{{R}_{A}}(t)\simeq {{R}_{S}}(t)</math> , one would compute  <math>n</math>  pairs of reliability and time values and then fit a single distribution to the data, or:


<br>
<br>
To compute an approximate reliability function for this system, <math>{{R}_{A}}(t)\simeq {{R}_{S}}(t)\,\!</math>, one would compute <math>n\,\!</math> pairs of reliability and time values and then fit a single distribution to the data, or:


::<math>\begin{align}
::<math>\begin{align}
Line 497: Line 224:
& ... \\  
& ... \\  
{{R}_{S}}(t= & 1,053.6)=90%   
{{R}_{S}}(t= & 1,053.6)=90%   
\end{align}</math>
\end{align}\,\!</math>


<br>
A single distribution, <math>{{R}_{A}}(t)\,\!</math>, that approximates <math>{{R}_{S}}(t)\,\!</math> can now be computed from these pairs using life data analysis methods.  If using the Weibull++ software, one would enter the values as free-form data.


A single distribution,  <math>{{R}_{A}}(t)</math> , that approximates  <math>{{R}_{S}}(t)</math>  can now be computed from these pairs using life data analysis methods.  If using the Weibull++ software, one would enter the values as free form data.
===Example===
{{:Example Using a Distribution to Approximate the CDF}}


===Example 3===
==Duty Cycle==
Compute a single Weibull distribution approximation for the system in Example 2.
Components of a system may not operate continuously during a system's mission, or may be subjected to loads greater or lesser than the rated loads during system operation. To model this, a factor called the Duty Cycle ( <math>{{d}_{c}}\,\!</math> ) is usedThe duty cycle may also be used to account for changes in environmental stress, such as temperature changes, that may effect the operation of a component. The duty cycle is a positive value, with a default value of 1 representing continuous operation at rated load, and any values other than 1 representing other load values with respect to the rated load value (or total operating time).   A duty cycle value higher than 1 indicates a load in excess of the rated value. A duty cycle value lower than 1 indicates that the component is operating at a load lower than the rated load or not operating continuously during the system's mission.  For instance, a duty cycle of 0.5 may be used for a component that operates only half of the time during the system's mission.
====Solution to Example 3====
The system in the previous example, shown in Figure fig4a, can be approximated by use of a 2-parameter Weibull distribution with  <math>\beta =2.02109</math>  and  <math>\eta =1123.51</math> .  In BlockSim, this is accomplished by representing the entire system as a block by going to the Subdiagram tab of the Block Properties window and selecting the Represent <br>
Diagram as Distribution option (Figure Ch5fig11).  
<br>
<br>
[[File:5.12.gif|center]]
<br>
<math></math>
<br>
<center>Figure 5.12: Representing a system with a distribution.</center>
<br>
[[File:BS5.13.png|center]]
<br>
<center>Figure 5.13: Distribution Settings window.</center>
<br>
Click the Approximate Failure Distribution Settings button, the Distribution Settings window will appear (Figure Ch5fig12).
<br>
In this window you can select a distribution to represent the data. BlockSim will then generate a number of system failure times based on the system's reliability function. The system's reliability function can be used to solve for a time value associated with that unreliability value. The distribution of the generated time values can then be fitted to a probability distribution function.
<br>


Consider a value of  <math>F(t)=0.11</math> .  Using the system's reliability equation and solving for time, the corresponding time-to-failure for a 0.11 unreliability can be calculated.  <br>
The reliability metrics for a component with a duty cycle are calculated as follows. Let <math>{{d}_{c}}\,\!</math> represent the duty cycle during a particular mission of the component, <math>t\,\!</math> represent the mission time and <math>{t}'\,\!</math> represent the accumulated age. Then:
For the system of Example 2, the time for a 0.11 unreliability is 389.786 hours. 
<br>
When enough points have been generated, the selected distribution will be fitted to this data set and the distribution's parameters will be returned.  In addition, if ReliaSoft's Weibull++ 7 is installed, the generated data can be viewed/analyzed using a Weibull++ instance, as shown in Figure Ch5fig13. 
<br>
[[File:BS5.14.png|center]]
<br>
<center>Figure 5.14: using Weibull++ to calculate distribution parameters.</center>
<br>
It is recommended that the analyst examine the fit to ascertain the applicability of the approximation.


=Duty Cycle=
::<math>{t}'={{d}_{c}}\times t\,\!</math>
 
Components of a system may not operate continuously during a system's mission, or may be subjected to loads greater or lesser than the rated loads during system operation.  To model this, a factor called the Duty Cycle ( <math>{{d}_{c}}</math> ) is used.  The duty cycle may also be used to account for changes in environmental stress, such as temperature changes, that may effect the operation of a component.  The duty cycle is a positive value, with a default value of 1 representing continuous operation at rated load, and any values other than 1 representing other load values with respect to the rated load value (or total operating time).  A duty cycle value higher than 1 indicates a load in excess of the rated value.  A duty cycle value lower than 1 indicates that the component is operating at a load lower than the rated load or not operating continuously during the system's mission.  For instance, a duty cycle of 0.5 may be used for a component that operates only half of the time during the system's mission.
<br>
The reliability metrics for a component with a duty cycle are calculated as follows. Let  <math>{{d}_{c}}</math>  represent the duty cycle during a particular mission of the component,  <math>t</math>  represent the mission time and  <math>{t}'</math>  represent the accumulated age. Then:
<br>
 
::<math>{t}'={{d}_{c}}\times t</math>
 
<br>


The reliability equation for the component is:
The reliability equation for the component is:
<br>
::<math>R({t}')=R({{d}_{c}}\times t)</math>
<br>


::<math>R({t}')=R({{d}_{c}}\times t)\,\!</math>


The component ''pdf'' is:
The component ''pdf'' is:
<br>
::<math>f({t}')=-\frac{d(R({t}'))}{dt}=-\frac{d(R({{d}_{c}}\times t))}{dt}={{d}_{c}}f({{d}_{c}}\times t)</math>
<br>


::<math>f({t}')=-\frac{d(R({t}'))}{dt}=-\frac{d(R({{d}_{c}}\times t))}{dt}={{d}_{c}}f({{d}_{c}}\times t)\,\!</math>


The failure rate of the component is:
The failure rate of the component is:
<br>
::<math>\lambda ({t}')=\frac{f({t}')}{R({t}')}=\frac{{{d}_{c}}f({{d}_{c}}\times t)}{R({{d}_{c}}\times t)}={{d}_{c}}\lambda ({{d}_{c}}\times t)</math>
===Example 4===
Consider a computer system with three components: a processor, a hard drive and a CD drive in series as shown next.  Assume that all three components follow a Weibull failure distribution with the parameters  <math>{{\beta }_{1}}=1.5</math>  and  <math>{{\eta }_{1}}=5000</math>  for the processor,  <math>{{\beta }_{2}}=2.5</math>  and  <math>{{\eta }_{2}}=3000</math>  for the hard drive, and  <math>{{\beta }_{3}}=2</math>  and  <math>{{\eta }_{3}}=4000</math>  for the CD drive.  Determine the reliability of the computer system after one year (365 days) of operation, assuming that the CD drive is used only 30% of the time.
<br>
[[File:BS5ex4.png|center]]
<br>
<br>
====Solution to Example 4====
<br>
The reliability of the processor after 365 days of operation is given by:
::<math>\begin{align}
  {{R}_{processor}}(365)= & {{e}^{-{{\left( \tfrac{365}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}} \\
  = & {{e}^{-{{\left( \tfrac{365}{5000} \right)}^{1.5}}}} \\
  = & 0.9805\text{ or }98.05% 
\end{align}</math>
The reliability of the hard drive after 365 days of operation is given by:
::<math>\begin{align}
  {{R}_{harddrive}}(365)= & {{e}^{-{{\left( \tfrac{365}{{{\eta }_{2}}} \right)}^{{{\beta }_{2}}}}}} \\
  = & {{e}^{-{{\left( \tfrac{365}{3000} \right)}^{2.5}}}} \\
  = & 0.9948\text{ or }99.48% 
\end{align}</math>
The reliability of the CD drive after 365 days of operation (taking into account the 30% operation using a duty cycle of 0.3) is given by:
::<math>\begin{align}
  {{R}_{CDdrive}}(365)= & {{e}^{-{{\left( \tfrac{{{d}_{c}}\times 365}{{{\eta }_{3}}} \right)}^{{{\beta }_{3}}}}}} \\
  = & {{e}^{-{{\left( \tfrac{0.3\times 365}{4000} \right)}^{2}}}} \\
  = & 0.9993\text{ or }99.93% 
\end{align}</math>
<br>
[[File:5-15.png|center]]
<br>
<center>Figure 5.15: Result for the computer system reliability.</center>
<br>
Thus the reliability of the computer system after 365 days of operation is:
::<math>\begin{align}
  {{R}_{s}}(365)= & {{R}_{processor}}(365)\cdot {{R}_{harddrive}}(365)\cdot {{R}_{CDdrive}}(365) \\
  = & 0.9805\cdot 0.9948\cdot 0.9993 \\
  = & 0.9747\text{ or }97.47% 
\end{align}</math>
<br>
This result can be obtained in BlockSim as shown in Figure DutyCycleExResults.
<br>


=Load Sharing=
::<math>\lambda ({t}')=\frac{f({t}')}{R({t}')}=\frac{{{d}_{c}}f({{d}_{c}}\times t)}{R({{d}_{c}}\times t)}={{d}_{c}}\lambda ({{d}_{c}}\times t)\,\!</math>
As presented in earlier chapters, a reliability block diagram (RBD) allows you to graphically represent how the components within a system are reliability-wise connected.  In most cases, independence is assumed across the components within the system.  For example, the failure of component A does not affect the failure of component B.  However, if a system consists of components that are sharing a load, then the assumption of independence no longer holds true.
<br>


If one component fails, then the component(s) that are still operating will have to assume the failed unit's portion of the load. Therefore, the reliabilities of the surviving unit(s) will change. Calculating the system reliability is no longer an easy proposition.  In the case of load sharing components, the change of the failure distributions of the surviving components must be known in order to determine the system's reliability.
===Example===<!-- THIS SECTION HEADER IS LINKED TO The PAGE TRANSCLUDED IN THIS SECTION. IF YOU RENAME THE SECTION, YOU MUST UPDATE THE LINK. -->
<br>
{{:Example Calculating System Reliability with Duty Cycles}}


To illustrate this, consider the a system of two units connected reliability-wise in parallel (Figure ch5loadparallel).
==Load Sharing==
<br>
As presented in earlier chapters, a reliability block diagram (RBD) allows you to graphically represent how the components within a system are reliability-wise connected. In most cases, independence is assumed across the components within the system. For example, the failure of component A does not affect the failure of component B. However, if a system consists of components that are sharing a load, then the assumption of independence no longer holds true.
 
[[File:5-16.png|center]]
<center>Figure 5.16: Two units connected reliability-wise in parallel.</center>
<br>
Assume that the units must supply an output of 8 volts and that if both units are operational, each unit is to supply 50% of the total output.  If one of the units fails, then the surviving unit supplies 100%. Furthermore, assume that having to supply the entire load has a negative impact on the reliability characteristics of the surviving unit. Since the reliability characteristics of the unit change based on whether both or only one is operating, a life distribution along with a life-stress relationship (as discussed in Chapter 3) will be needed to model each component.
<br>


To illustrate the steps needed, we will create the model starting from raw data. Assume that a total of 20 units were tested to failure at 7, 10 and 15 volts.  The test data set is presented in the next table.
If one component fails, then the component(s) that are still operating will have to assume the failed unit's portion of the load. Therefore, the reliabilities of the surviving unit(s) will change. Calculating the system reliability is no longer an easy proposition. In the case of load sharing components, the change of the failure distributions of the surviving components must be known in order to determine the system's reliability.
<br>


[[File:5-17.png|center]]
To illustrate this, consider a system of two units connected reliability-wise in parallel as shown below.
<br>
[[Image:5-16.png|center|400px|<div align="center"> Two units connected reliability-wise in parallel.</div>|link=]]
Assume that the units must supply an output of 8 volts and that if both units are operational, each unit is to supply 50% of the total output. If one of the units fails, then the surviving unit supplies 100%. Furthermore, assume that having to supply the entire load has a negative impact on the reliability characteristics of the surviving unit.  
For this example, Units 1 and 2 are the same component. Therefore, only one set of data was collected. However, it is possible that the load sharing components in a system may not be the same.  If that were the case, data would need to be collected for each component.
<br>


The data set in Table 1 was analyzed using ReliaSoft's ALTA software (as shown in Figure fig14) with the Inverse Power Law as the underlying life-stress relationship and Weibull as the life distribution.
Because the reliability characteristics of the unit change based on the load it is sharing, a method that can model the effect of the load on life should be used. One way to do this is to use a life distribution along with a life-stress relationship (as discussed in [[Statistical_Background#A_Brief_Introduction_to_Life-Stress_Relationships|A Brief Introduction to Life-Stress Relationships]]) for each component. The detailed discussion for this method can be found at [[Additional Information on Load Sharing]]. Another simple way is to use the concept of acceleration factors and assume that the load has a linear effect on the failure time. If the load is doubled, then the life of the component will be shortened by half.  
<br>


::<math></math>
For the above load sharing system, the reliability of each component is a function of time and load. For example, for Unit 1, the reliability and the probability density function are:
<br>


The estimated model parameters,  <math>\beta </math> , <math>K</math> and <math>n</math> , are shown next.
::<math>{{R}_{1}}(t,{{S}_{1}})\,\!\,\!</math> and <math>{{f}_{1}}(t,{{S}_{1}})\,\!\,\!</math>
<br>


where <math>{{S}_{1}}\,\!\,\!</math> is the load shared by Unit 1 at time ''t'' and the total load of the system is <math>S={{S}_{1}}+{{S}_{2}}\,\!\,\!</math>. At the beginning, both units are working. Assume that Unit 1 fails at time ''x'' and Unit 2 takes over the entire load. The reliability for Unit 2 at time ''x'' is:


::<math>\begin{align}
::<math>\begin{align}
\beta = & 1.9239 \\  
  & {{R}_{2}}(x,{{S}_{2}})={{R}_{2}}({{t}_{2e}},S) \\  
K= & 3.2387\times {{10}^{-7}} \\  
& {{t}_{2e}}=\frac{{{S}_{2}}}{S}x \\  
n= & 3.4226 
\end{align}\,\!\,\!</math>
\end{align}</math>
<br>


Or:  
<math>{{t}_{2e}}\,\!\,\!</math> is the equivalent time for Unit 2 at time ''x'' if it is operated with load ''S''. The equivalent time concept is illustrated in the following plot.
<br>
[[Image:BS5.19.png|center|600px| Illustrating <math>{{t}_{e}}\,\!</math>|link=]]
The system reliability at time ''t'' is:


::<math>\begin{align}
::<math>\begin{align}
{{R}_{1}}(t,{{S}_{1}})= & {{e}^{-{{\left( KS_{1}^{n}t \right)}^{\beta }}}} \\  
  & R(t)={{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}}) \\
= & {{e}^{-{{\left( 3.2387\times {{10}^{-7}}S_{1}^{3.4226}t \right)}^{1.9239}}}}   
& +\int_{0}^{t}{{{f}_{1}}(x,{{S}_{1}})\cdot {{R}_{2}}({{t}_{2e}}+(t-x),S)dx} \\  
\end{align}</math>
& +\int_{0}^{t}{{{f}_{2}}(x,{{S}_{2}})\cdot {{R}_{1}}({{t}_{1e}}+(t-x),S)dx}   
\end{align}\,\!\,\!</math>


<br>
In BlockSim, the failure time distribution for each component is defined at the load of ''S''. The reliability function for a component at a given load is calculated as:


::<math>{{f}_{1}}(t,{{S}_{1}})=\beta KS_{1}^{n}{{\left( KS_{1}^{n}t \right)}^{\beta -1}}{{e}^{-{{\left( KS_{1}^{n}t \right)}^{\beta }}}}</math>
::<math>{{R}_{i}}(t,{{S}_{i}})={{R}_{i}}(t\times \frac{{{S}_{i}}}{S},S)\,\!\,\!</math>
::<math>{{f}_{i}}(t,{{S}_{i}})=\frac{{{S}_{i}}}{S}{{f}_{i}}\left( t\times \frac{{{S}_{i}}}{S},S \right)\,\!\,\!</math>
From the above equation, it can be seen that the concept used in the calculation for load sharing is the same as the concept used in the calculation for [[Time-Dependent_System_Reliability_(Analytical)#Duty_Cycle|duty cycle]].


<br>
===Example===
In the following load sharing system, Block 1 follows a Weibull failure distribution with<math>{{\beta }_{1}}=1.5\,\!</math>, and<math>{{\eta }_{1}}=1,000\,\!</math>
. Block 2 follows a Weibull failure distribution with<math>{{\beta }_{2}}=2\,\!</math>, and<math>{{\eta }_{2}}=2,000\,\!</math>. The load for Block 1 is 1 unit, and for Block 2 it is 3 units. Calculate the system reliability at time 1,500.


And for this case:
[[Image:loadsharingconfig.png|center|345px| Load Sharing System|link=]]
Block 1 shares 25% (P1) of the entire load, and Block 2 shares 75% (P2) of it. Therefore, we have the following equations for calculating the system reliability:


::<math>{{R}_{i}}(t,{{S}_{i}})={{R}_{i}}({{P}_{i}}\times t)\,\!</math>,
::<math>{{f}_{i}}(x,{{S}_{i}})={{P}_{i}}{{f}_{i}}({{P}_{i}}\times x)\,\!</math>


::<math>\begin{align}
and:
{{R}_{1}}(t,{{S}_{1}})= & {{R}_{2}}(t,{{S}_{2}}) \\
{{f}_{1}}(t,{{S}_{1}})= & {{f}_{2}}(t,{{S}_{2}}) 
\end{align}</math>
<br>
<br>
[[File:5.17.gif|center]]
<br>
<center>Figure 5.17: Using ALTA to calculate component parameters.</center>
<br>
 
Figure fig15 shows a plot of Eqn. (iplweib).
<br>


Now that the failure properties have been determined using the test data, the reliability of the system at some time,  <math>t</math> , can be calculated using the following equation:
::<math>{{R}_{1}}({{t}_{1e}}+(t-x))={{R}_{1}}({{P}_{1}}x+t-x)={{R}_{1}}(t-{{P}_{2}}x)\,\!</math>
::<math>{{R}_{2}}({{t}_{2e}}+(t-x))={{R}_{2}}({{P}_{2}}x+t-x)={{R}_{2}}(t-{{P}_{1}}x)\,\!</math>


<br>
Using the above equations in the system reliability function, we get:


::<math>\begin{align}
::<math>\begin{align}
R(t,S)= & {{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}}) \\
  & R(t)={{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}})+\int_{0}^{t}{{{f}_{1}}(x,{{S}_{1}})\cdot {{R}_{2}}({{t}_{2e}}+(t-x),S)dx} \\
& +\underset{o}{\overset{t}{\mathop \int }}\,{{f}_{1}}\left( x,{{S}_{1}} \right)\cdot {{R}_{2}}(x,{{S}_{2}})\cdot \left( \frac{{{R}_{2}}({{t}_{1e}}+(t-x),S)}{{{R}_{2}}({{t}_{1e}},S)} \right)dx \\
& +\int_{0}^{t}{{{f}_{2}}(x,{{S}_{2}})\cdot {{R}_{1}}({{t}_{1e}}+(t-x),S)dx} \\
& +\underset{o}{\overset{t}{\mathop \int }}\,{{f}_{2}}\left( x,{{S}_{2}} \right)\cdot {{R}_{1}}(x,{{S}_{1}})\cdot \left( \frac{{{R}_{1}}({{t}_{2e}}+(t-x),S)}{{{R}_{1}}({{t}_{2e}},S)} \right)dx  
& ={{R}_{1}}({{P}_{1}}\times t)\cdot {{R}_{2}}({{P}_{2}}\times t)+\int_{0}^{t}{{{P}_{1}}{{f}_{1}}({{P}_{1}}x)\cdot {{R}_{2}}(t-{{P}_{1}}x)dx} \\  
\end{align}</math>
& +\int_{0}^{t}{{{P}_{2}}{{f}_{2}}({{P}_{2}}x)\cdot {{R}_{1}}(t-{{P}_{2}}x)dx}   
\end{align}\,\!</math>


<br>
The calculated system reliability at time 1,500 is 0.8569, as given below.
[[Image:loadsharingresults.png|center|450px|Calculated System Reliability at Time 1,500|link=]]


Where:  
==Standby Components== <!-- THIS SECTION HEADER IS LINKED TO: Standby_Configuration_Example. IF YOU RENAME THE SECTION, YOU MUST UPDATE THE LINK. -->
<br>
In the previous section, the case of a system with load sharing components was presented.  This is a form of redundancy with dependent components. That is, the failure of one component affects the failure of the other(s).  This section presents another form of redundancy: standby redundancy.  In standby redundancy the redundant components are set to be under a lighter load condition (or no load) while not needed and under the operating load when they are activated.


::<math>\begin{align}
In standby redundancy the components are set to have two states: an active state and a standby state.  Components in standby redundancy have two failure distributions, one for each state. When in the standby state, they have a quiescent (or dormant) failure distribution and when operating, they have an active failure distribution.
{{S}_{1}}= & {{P}_{1}}S \\
{{S}_{2}}= & {{P}_{2}}S  
\end{align}</math>


<br>
In the case that both quiescent and active failure distributions are the same, the units are in a simple parallel configuration (also called a hot standby configuration). When the rate of failure of the standby component is lower in quiescent mode than in active mode, that is called a warm standby configuration. When the rate of failure of the standby component is zero in quiescent mode (i.e., the component cannot fail when in standby), that is called a cold standby configuration.
[[File:5-18.png|center]]
<br>
<center>Figure 5.18: Reliability curves for different voltage output conditions.</center>
<br>
[[File:BS5.19.png|center]]
<br>
<br>


And:
===Simple Standby Configuration===
<br>
Consider two components in a standby configuration.  Component 1 is the active component with a Weibull failure distribution with parameters <math>\beta = 1.5\,\!</math> and <math>\eta  = 1,000 \,\!</math>.  Component 2 is the standby component.  When Component 2 is operating, it also has a Weibull failure distribution with <math>\beta  = 1.5\,\!</math> and <math>\eta  = 1,000 \,\!</math>.  Furthermore, assume the following cases for the quiescent distribution.


:• <math>S</math>  is the total load (or required output).<br>
:*Case 1:  The quiescent distribution is the same as the active distribution (hot standby).
:• <math>{{P}_{1}}</math>  and  <math>{{P}_{2}}</math>  are the portion of the total load that each unit supports when both units are operational.  In this case,  <math>{{P}_{1}}={{P}_{2}}=0.5=50%.</math><br>
:• <math>{{S}_{1}}</math>  and  <math>{{S}_{2}}</math>  represent the portions of the load that Unit 1 and Unit 2 must support when both units are operational.<br>
:• <math>{{t}_{{{1}_{e}}}}</math> is the equivalent operating time for Unit 1 if it had been operating at  <math>S</math>  instead of  <math>{{S}_{1}}</math> .  A graphical representation of the equivalent time is shown in Figure fig16, where the curve marked by L represents the low stress (load) and the curve marked by H represents the high stress (load).<br>
<br>
::<math>{{t}_{1e}}</math>  can be calculated by:
<br>


:*Case 2:  The quiescent distribution is a Weibull  distribution with <math>\beta = 1.5\,\!</math> and <math>\eta = 2000\,\!</math> (warm standby).


::<math>\begin{align}
:*Case 3: The component cannot fail in quiescent mode (cold standby).
{{R}_{1}}(t)= & {{R}_{1}}({{t}_{1e}}) \\
{{e}^{-{{(tKS_{1}^{n})}^{\beta }}}}= & {{e}^{-{{({{t}_{1e}}K{{S}^{n}})}^{\beta }}}} \\
tS_{1}^{n}= & {{t}_{1e}}{{S}^{n}} \\
{{t}_{1e}}= & t{{\left( \frac{{{S}_{1}}}{S} \right)}^{n}},\text{    }{{S}_{1}}={{P}_{1}}S \\
\therefore  & {{t}_{1e}}=tP_{1}^{n} 
\end{align}</math>


<br>
In this case, the reliability of the system at some time, <math>t\,\!</math>, can be obtained using the following equation:


::<math>{{t}_{2e}}</math> can be calculated the same way, or:
::<math>R(t)={{R}_{1}}(t)+\underset{0}{\overset{t}{\mathop \int }}\,{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x)\cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}dx  \ \,\!</math>


<br>
where:


::<math>\begin{align}
:*<math>{{R}_{1}}\,\!</math> is the reliability of the active component.
{{R}_{2}}(t)= & {{R}_{2}}({{t}_{2e}}) \\  
\therefore  & {{t}_{2e}}=tP_{2}^{n} 
\end{align}</math>
<br>


:*<math>{{f}_{1}}\,\!</math> is the ''pdf'' of the active component.


In this example, the reliability equations for Unit 1 and Unit 2 are the same since they are the same type of component and demonstrate the same failure properties.  In addition, the total output is divided equally between the two units (when both units are operating), so  <math>{{t}_{1e}}</math> and  <math>{{t}_{2e}}</math>  will also be the same.
:*<math>{{R}_{2;SB}}\,\!</math> is the reliability of the standby component when in standby mode (quiescent reliability).
<br>
The next step is to determine the reliability of the system after 8,760 hours,  <math>R(t=8,760)</math> .  Using Eqn. (loadshare3) the system reliability is found to be:
<br>


:*<math>{{R}_{2;A}}\,\!</math> is the reliability of the standby component when in active mode.


:*<math>{{t}_{e}}\,\!</math> is the equivalent operating time for the standby unit if it had been operating at an active mode, such that:
::<math>\begin{align}
::<math>\begin{align}
R(t=8760)= & 0.8567 \\
{{R}_{2;SB}}(x)={{R}_{2;A}}({{t}_{e}})   
= & 85.67%  
\end{align}\,\!</math>
\end{align}</math>
<br>


The second equation above can be solved for <math>{{t}_{e}}\,\!</math> and substituted into the first equation above.
The following figure illustrates the example as entered in BlockSim using a standby container.


===Load Sharing in BlockSim===
[[Image:5_24_new.png|center|200px|<div align="center"> Standby container.</div>|link=]]


BlockSim uses this formulation when computing reliabilities of units in a load sharing configurationWhen using the System Reliability Equation window, BlockSim returns a single token for the reliability of units in a load sharing configuration (as well as in the case of standby redundancy, discussed in the next section)As an example, consider the following RBD with Unit 1 in series with a container that includes two load sharing units.
The active and standby blocks are within a container, which is used to specify standby redundancySince the standby component has two distributions (active and quiescent), the Block Properties window of the standby block has two pages for specifying each oneThe following figures illustrate these pages.
<br>
[[Image: Fig 5.25.PNG|center|600px|<div align="center"> Defining the active failure distribution </div>|link=]]
[[File:BS5.19.2.png|center]]
<br>  
BlockSim will return the system equation as:
<br>


::<math>{{R}_{System}}=+{{R}_{LS}}\cdot {{R}_{1}}</math>
<br>


Where  <math>{{R}_{LS}}</math>  implies a form similar to Eqn. (loadshare3). BlockSim allows for  <math>k</math> -out-of- <math>n</mathunits in a load sharing configuration.
[[Image:Fig 5.26.PNG|center|600px|<div align="center"> Defining the quiescent failure distribution </div>|link=]]
   
The system reliability results for 1000 hours are given in the following table:


===Example 5===
[[Image:5-24.png|center|250px|link=|]]
A component has five possible failure modes,  <math>A</math> ,  <math>{{B}_{A}}</math> ,  <math>{{B}_{B}}</math> ,  <math>{{B}_{C}}</math>  and  <math>C</math> , and the  <math>B</math>  modes are interdependent.  The system will fail if mode  <math>A</math>  occurs, mode  <math>C</math>  occurs or two out of the three  <math>B</math>  modes occur.
<br>
 
Modes  <math>A</math>  and  <math>C</math>  each have a Weibull distribution, with a  <math>\beta =2</math>  and  <math>\eta =10,000</math>  and 15,000 respectively.  Events  <math>{{B}_{A}}</math> ,  <math>{{B}_{B}}</math>  and  <math>{{B}_{C}}</math>  each have an exponential distribution with a mean of 10,000 hours.
<br>
 
If any  <math>B</math>  event occurs (i.e.  <math>{{B}_{A}}</math> ,  <math>{{B}_{B}}</math>  or  <math>{{B}_{C}}</math> ), the remaining  <math>B</math>  events are more likely to occur.  Specifically, the mean times of the remaining  <math>B</math>  events are halved.  Determine the reliability at 1000 hours for this component.
====Solution to Example 5====
The first step is to create the RBD.  Modes  <math>A</math> ,  <math>C</math>  and a load sharing container with the  <math>{{B}_{i}}</math>  modes must be drawn in series, as illustrated next.
<br>
[[File:BS5.19.3.png|center]]
<br>
   
   
The next step is to define the properties for each block, including those for the container.  Setting the failure distributions for modes  <math>A</math> and  <math>C</math>  is simple.  The more difficult part is setting the properties for the container and the contained blocks. Based on the problem statement, the <math>B</math>  modes are in a 2-out-of-3 load sharing redundancy.  When all three are working (i.e. when no  <math>B</math>  mode has occurred), each block has an exponential distribution with  .. If one  <math>B</math>  mode occurs, then the two surviving units have an exponential distribution with  <math>\mu =5,000.</math>
Note that even though the <math>\beta \,\!</math> value for the quiescent distribution is the same as in the active distribution, it is possible that the two can be different. That is, the failure modes present during the quiescent mode could be different from the modes present during the active mode.  In that sense, the two distribution types can be different as well (e.g., lognormal when quiescent and Weibull when active).
<br>


Assume a Power Life-Stress relationship for the componentsThen:
In many cases when considering standby systems, a switching device may also be present that switches from the failed active component to the standby componentThe reliability of the switch can also be incorporated into
<br>
<math>R(t)={{R}_{1}}(t)+\underset{0}{\overset{t}{\mathop \int }}\,{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x)\cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}dx  \ \,\!</math>


::<math>\begin{align}
as presented in the next section.
{{\mu }_{1}}= & \frac{1}{KV_{1}^{n}} \\
{{\mu }_{2}}= & \frac{1}{KV_{2}^{n}} 
\end{align}</math>
<br>


BlockSim's System Reliability Equation window returns a single token for the reliability of units in a standby configuration.  This is the same as the load sharing case presented in the previous section.


Substituting  <math>{{\mu }_{1}}=10,000</math>  and  <math>{{V}_{1}}=1</math> in Eqn. (miu1) and casting it in terms of  <math>K</math>  yields:
===Reliability of Standby Systems with a Switching Device===
<br>
In many cases when dealing with standby systems, a switching device is present that will switch to the standby component when the active component fails. Therefore, the failure properties of the switch must also be included in the analysis.


[[Image:BS5_26_2_new.png|center|250px|link=]]


::<math>\begin{align}
In most cases when the reliability of a switch is to be included in the analysis, two probabilities can be considered.  The first and most common one is the probability of the switch performing the action (i.e., switching) when requested to do so.  This is called Switch Probability per Request in BlockSim and is expressed as a static probability (e.g., 90%).  The second probability is the quiescent reliability of the switchThis is the reliability of the switch as it ages (e.g., the switch might wear out with age due to corrosion, material degradation, etc.). Thus it is possible for the switch to fail before the active component fails.  However, a switch failure does not cause the system to fail, but rather causes the system to fail only if the switch is needed and the switch has failed.  For example, if the active component does not fail until the mission end time and the switch fails, then the system does not fail.  However, if the active component fails and the switch has also failed, then the system cannot be switched to the standby component and it therefore fails.
10,000= & \frac{1}{K} \\
K = & \frac{1}{10,000}=0.0001  
\end{align}</math>


<br>
In analyzing standby components with a switching device, either or both failure probabilities (during the switching or while waiting to switch) can be considered for the switch, since each probability can represent different failure modes.  For example, the switch probability per request may represent software-related issues or the probability of detecting the failure of an active component, and the quiescent probability may represent wear-out type failures of the switch.


Substituting <math>{{\mu }_{2}}=5,000</math> , <math>{{V}_{2}}=1.5</math> (because if one fails, then each survivor takes on an additional 0.5 units of load) and Eqn. (Kval) for  <math>K</math> in Eqn. (mieu2) yields:
To illustrate the formulation, consider the previous example that assumes perfect switching. To examine the effects of including an imperfect switch, assume that when the active component fails there is a 90% probability that the switch will switch from the active component to the standby component.  In addition, assume that the switch can also fail due to a wear-out failure mode described by a Weibull distribution with <math>\beta = 1.7\,\!</math> and <math>\eta = 5000\,\!</math>.


<br>
Therefore, the reliability of the system at some time, <math>t\,\!</math>, is given by the following equation.


::<math>\begin{align}
::<math>\begin{align}
5,000= & \frac{1}{0.0001\cdot {{(1.5)}^{n}}} \\  
R(t)= & {{R}_{1}}(t) +\underset{0}{\overset{t}{\mathop \int }}\,\{{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x) \\  
0.5= & {{(1.5)}^{-n}} \\
& \cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}\cdot {{R}_{SW;Q}}(x)\cdot {{R}_{SW;REQ}}(x)\}dx  
\ln (0.5)= & -n\ln (1.5) \\
\end{align}\,\!</math>
n = & 1.7095  
\end{align}</math>
<br>


where:


This also could have been computed in ALTA, as shown in Figure fig17, or with the Load & Life Parameter Experimenter in BlockSim, as shown in Figure Ch5fig18.
:*<math>{{R}_{1}}\,\!</math> is the reliability of the active component.
<br>


::<math></math>
:*<math>{{f}_{1}}\,\!</math> is the ''pdf'' of the active component.
<br>


At this point, the parameters for the load sharing units have been computed and can be set, as shown in Figure Ch5fig19.
:*<math>{{R}_{2;SB}}\,\!</math> is the reliability of the standby component when in standby mode (quiescent reliability).
<br>


::<math></math>
:*<math>{{R}_{2;A}}\,\!</math> is the reliability of the standby component when in active mode.
<br>


The next step is to set the weight proportionality factor.  This factor defines the portion of the load that the particular item carries while operating, as well as the load that shifts to the remaining units upon failure of the item.  To illustrate, assume three units (1, 2 and 3) are in a load sharing container with weight proportionality factors of 1, 2 and 3 respectively (and a 1-out-of-3 requirement).
:*<math>{{R}_{SW;Q}}\,\!</math> is the quiescent reliability of the switch.
<br>


:• Unit 1 carries  <math>\left( \tfrac{1}{1+2+3} \right)=0.166</math>  or 16.6% of the total load.
:*<math>{{R}_{SW;REQ}}\,\!</math> is the switch probability per request.
:• Unit 2 carries  <math>\left( \tfrac{2}{1+2+3} \right)=0.333</math>  or 33.3% of the total load.
:• Unit 3 carries  <math>\left( \tfrac{3}{1+2+3} \right)=0.50</math> or 50% of the total load.
<br>


The actual load on each unit then becomes the product of the entire load defined for the container times the portion carried by that unit.  For example, if the container load is 100 lbs, then the portion assigned to Unit 1 will be  <math>100\cdot 0.166=16.6</math> lbs.
:*<math>{{t}_{e}}\,\!</math> is the equivalent operating time for the standby unit if it had been operating at an active mode.
<br>


In the current example, all units share the same load and thus have equal weight proportionality factors.  Because these factors are relative, if the same number is used for all three items then the results will be the same.  For simplicity, factor is set equal to 1 for each item.
This problem can be solved in BlockSim by including these probabilities in the container's properties, as shown in the figures below.  In BlockSim, the standby container is acting as the switch.
<br>
 
[[File:BS5.20.png|center]]
<br>
<center>Figure 5.20: Calculation performed in ALTA.</center>
<br>
<br>
 
[[File:5.21.gif|center]]
<br>
<center>Figure 5.21: blockSim's Load & Life Parameter Experimenter.</center>
<br>
<br>
[[File:5.22.gif|center]]
<center>Figure 5.22: Load sharing parameters.</center>
<br>
 
<br>
[[File:BS5.23.png|center]]
<br>
<center>Figure 5.23: Defining total load for load sharing units.</center>
<br>
 
The last properties that need to be defined are the total load and the 2-out-of-3 redundancy.  The total load is dependent on how the parameters were computed.  In this case, total load was assumed to be 3 when the parameters were computed (i.e. the load per item was 1 when all worked and 1.5 when two worked).  This is defined at the container level, as shown in Figure fig20.
   
   
When all of the parameters have been specified in BlockSim, the reliability at 1,000 hours can be determined.  From the Analytical QCP, this is found to be 93.87%.
[[Image:Fig 5.28_2.PNG|center|600px|<div align="center"> Standby container (switch) failure probabilities while attempting to switch </div>|link=]]
 
=Standby Components=
 
In the previous section, the case of a system with load sharing components was presented.  This is a form of redundancy with dependent components. That is, the failure of one component affects the failure of the other(s).  This section presents another form of redundancy: standby redundancy.  In standby redundancy the redundant components are set to be under a lighter load condition (or no load) while not needed and under the operating load when they are activated.
<br>
 
In standby redundancy the components are set to have two states: an active state and a standby state.  Components in standby redundancy have two failure distributions, one for each state.  When in the standby state, they have a quiescent (or dormant) failure distribution and when operating, they have an active failure distribution.
<br>
 
In the case that both quiescent and active failure distributions are the same, the units are in a simple parallel configuration (also called a hot standby configuration).  When the rate of failure of the standby component is lower in quiescent mode than in active mode, that is called a warm standby configuration.  When the rate of failure of the standby component is zero in quiescent mode (i.e. the component cannot fail when in standby), that is called a cold standby configuration. 
 
===Simple Standby Configuration===
 
Consider two components in a standby configuration.  Component 1 is the active component with a Weibull failure distribution with parameters  <math>\beta </math>  = 1.5 and  <math>\eta </math>  = 1,000.  Component 2 is the standby component.  When Component 2 is operating, it also has a Weibull failure distribution with  <math>\beta </math>  = 1.5 and  <math>\eta </math>  = 1,000.  Furthermore, assume the following cases for the quiescent distribution.
<br>
:• Case 1:  The quiescent distribution is the same as the active distribution (hot standby).<br>
:• Case 2:  The quiescent distribution is a Weibull  distribution with  <math>\beta </math>  = 1.5 and  <math>\eta </math>  = 2000 (warm standby).<br>
:• Case 3: The component cannot fail in quiescent mode (cold standby).<br>
<br>
In this case, the reliability of the system at some time,  <math>t</math> , can be obtained using the following equation:
<br>
 
 
::<math>R(t)={{R}_{1}}(t)+\underset{0}{\overset{t}{\mathop \int }}\,{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x)\cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}dx</math>
<br>
 
Where:
<br>
:• <math>{{R}_{1}}</math>  is the reliability of the active component.<br>
:• <math>{{f}_{1}}</math>  is the  <math>pdf</math>  of the active component.<br>
:• <math>{{R}_{2;SB}}</math>  is the reliability of the standby component when in standby mode (quiescent reliability).<br>
:• <math>{{R}_{2;A}}</math>  is the reliability of the standby component when in active mode.<br>
:• <math>{{t}_{e}}</math>  is the equivalent operating time for the standby unit if it had been operating at an active mode, such that: <br>
:• <math>{{R}_{2;SB}}(x)={{R}_{2;A}}({{t}_{e}})</math><br>
<br>
Eqn. (stbte) can be solved for  <math>{{t}_{e}}</math>  and substituted into Eqn. (stb1).
Figure ch5StandbyCont illustrates the example as entered in BlockSim using a standby container.
<br>
<br>
[[File:5.24.gif|center]]
<br>
<center>Figure 5.24: Standby container.</center>
<br>
<math></math>
<br>
 
The active and standby blocks are within a container, which is used to specify standby redundancy.  Since the standby component has two distributions (active and quiescent), the Block Properties window of the standby block has two pages for specifying each one.  Figures fig21 and fig22 illustrate these pages.
<br>
The system reliability results for 1000 hours are given in the following table:
<br>
[[File:5-24.png|center]]
<br>
Note that even though the  <math>\beta </math>  value for the quiescent distribution is the same as in the active distribution, it is possible that the two can be different. That is, the failure modes present during the quiescent mode could be different from the modes present during the active mode.  In that sense, the two distribution types can be different as well (e.g. lognormal when quiescent and Weibull when active).
<br>
 
In many cases when considering standby systems, a switching device may also be present that switches from the failed active component to the standby component.  The reliability of the switch can also be incorporated into Eqn. (stb1), as presented in the next section.
<br>
 
BlockSim's System Reliability Equation window returns a single token for the reliability of units in a standby configuration.  This is the same as the load sharing case presented in the previous section.
<br>
 
[[File:BS5.25.png|center]]
<br>
<br>
 
[[File:BS5.26.png|center]]
<br>
 
===Reliability of Standby Systems with a Switching Device===
 
In many cases when dealing with standby systems, a switching device is present that will switch to the standby component when the active component fails.  Therefore, the failure properties of the switch must also be included in the analysis.
<br>
 
[[File:BS5.26.2.png|center]]  
<br>


In most cases when the reliability of a switch is to be included in the analysis, two probabilities can be considered.  The first and most common one is the probability of the switch performing the action (i.e. switching) when requested to do so.  This is called Switch Probability per Request in BlockSim and is expressed as a static probability (e.g. 90%).  The second probability is the quiescent reliability of the switch.  This is the reliability of the switch as it ages (e.g. the switch might wear out with age due to corrosion, material degradation, etc.). Thus it is possible for the switch to fail before the active component fails.  However, a switch failure does not cause the system to fail, but rather causes the system to fail only if the switch is needed and the switch has failed.  For example, if the active component does not fail until the mission end time and the switch fails, then the system does not fail.  However, if the active component fails and the switch has also failed, then the system cannot be switched to the standby component and it therefore fails.
<br>


In analyzing standby components with a switching device, either or both failure probabilities (during the switching or while waiting to switch) can be considered for the switch, since each probability can represent different failure modes.  For example, the switch probability per request may represent software-related issues or the probability of detecting the failure of an active component, and the quiescent probability may represent wear-out type failures of the switch.
[[Image:Fig 5.27.PNG|center|600px|<div align="center"> Standby container (switch) failure distribution while waiting to switch </div>|link=]]
<br>


To illustrate the formulation, consider the previous example that assumes perfect switchingTo examine the effects of including an imperfect switch, assume that when the active component fails there is a 90% probability that the switch will switch from the active component to the standby component.  In addition, assume that the switch can also fail due to a wear-out failure mode described by a Weibull distribution with  <math>\beta </math> = 1.7 and <math>\eta </math>  = 5000.
Note that there are additional properties that can be specified in BlockSim for a switch, such as Switch Restart Probability, No. of Restarts and Switch Delay TimeIn many applications, the switch is re-tested (or re-cycled) if it fails to switch the first time.  In these cases, it might be possible that it switches in the second or third, or <math>{{n}^{th}}\,\!</math> attempt.   
<br>


Therefore, the reliability of the system at some time, <math>t</math> , is given by the following equation.
The Switch Restart Probability specifies each additional attempt's probability of successfully switching and the Finite Restarts specifies the total number of attempts.  Note that the Switch Restart Probability specifies the probability of success of each trial (or attempt). The probability of success of <math>n\,\!</math> consecutive trials is calculated by BlockSim using the binomial distribution and this probability is then incorporated into equation above.  The Switch Delay Time property is related to repairable systems and is considered in BlockSim only when using simulation.  When using the analytical solution (i.e., for a non-repairable system), this property is ignored.


Solving the analytical solution (as given by the above equation), the following results are obtained.


::<math>\begin{align}
[[Image:5-30.png|center|600px]]
R(t)= & {{R}_{1}}(t) \\
& +\underset{0}{\overset{t}{\mathop \int }}\,\{{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x) \\
& \cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}\cdot {{R}_{SW;Q}}(x)\cdot {{R}_{SW;REQ}}(x)\}dx 
\end{align}</math>
 
Where:
<br>
:• <math>{{R}_{1}}</math>  is the reliability of the active component.<br>
:• <math>{{f}_{1}}</math>  is the  <math>pdf</math>  of the active component.<br>
:• <math>{{R}_{2;SB}}</math>  is the reliability of the standby component when in standby mode (quiescent reliability).<br>
:• <math>{{R}_{2;A}}</math>  is the reliability of the standby component when in active mode.<br>
:• <math>{{R}_{SW;Q}}</math>  is the quiescent reliability of the switch.<br>
:• <math>{{R}_{SW;REQ}}</math>  is the switch probability per request.<br>
:• <math>{{t}_{e}}</math>  is the equivalent operating time for the standby unit if it had been operating at an active mode.<br>
This problem can be solved in BlockSim by including these probabilities in the container's properties, as shown in Figures fig23 and fig24.  In BlockSim, the standby container is acting as the switch.
   
   
<br>
From the table above, it can be seen that the presence of a switching device has a significant effect on the reliability of a standby system.  It is therefore important when modeling standby redundancy to incorporate the switching device reliability properties.  It should be noted that this methodology is not the same as treating the switching device as another series component with the standby subsystem.  This would be valid only if the failure of the switch resulted in the failure of system (e.g., switch failing open).  In equation above, the Switch Probability per Request and quiescent probability are present only in the second term of the equation.  Treating these two failure modes as a series configuration with the standby subsystem would imply that they are also present when the active component is functioning (i.e., first term of equation above).  This is invalid and would result in the underestimation of the reliability of the system.  In other words, these two failure modes become significant only when the active component fails.
 
[[File:BS5.27.png|center]]
<br>
<br>
 
[[File:BS5.28.png|center]]
<br>
<br>
 
Note that there are additional properties that can be specified in BlockSim for a switch, such as Switch Restart Probability, Finite Restarts and Switch Delay Time.  In many applications, the switch is re-tested (or re-cycled) if it fails to switch the first time.  In these cases, it might be possible that it switches in the second or third, or  <math>{{n}^{th}}</math> attempt. 
<br>
The Switch Restart Probability specifies each additional attempt's probability of successfully switching and the Finite Restarts specifies the total number of attempts.  Note that the Switch Restart Probability specifies the probability of success of each trial (or attempt).  The probability of success of  <math>n</math>  consecutive trials is calculated by BlockSim using the binomial distribution and this probability is then incorporated into Eqn. (stb2a).  The Switch Delay Time property is related to repairable systems and is considered in BlockSim only when using simulation.  When using the analytical solution (i.e. for a non-repairable system), this property is ignored.
<br>
 
Solving the analytical solution (as given by Eqn. stb2a), the following results are obtained.
<br>
 
[[File:5-30.png|center]]
<br>
From the table above, it can be seen that the presence of a switching device has a significant effect on the reliability of a standby system.  It is therefore important when modeling standby redundancy to incorporate the switching device reliability properties.  It should be noted that this methodology is not the same as treating the switching device as another series component with the standby subsystem.  This would be valid only if the failure of the switch resulted in the failure of system (e.g. switch failing open).  In Eqn. (stb2a), the Switch Probability per Request and quiescent probability are present only in the second term of the equation.  Treating these two failure modes as a series configuration with the standby subsystem would imply that they are also present when the active component is functioning (i.e. first term of Eqn. stb2a).  This is invalid and would result in the underestimation of the reliability of the system.  In other words, these two failure modes become significant only when the active component fails.
<br>


As an example, and if we consider the warm standby case, the reliability of the system without the switch is 70.57% at 1000 hours.  If the system was modeled so that the switching device was in series with the warm standby subsystem, the result would have been:
As an example, and if we consider the warm standby case, the reliability of the system without the switch is 70.57% at 1000 hours.  If the system was modeled so that the switching device was in series with the warm standby subsystem, the result would have been:
<br>


::<math>\begin{align}
::<math>\begin{align}
Line 1,028: Line 435:
= & 0.7057\cdot 0.9372\cdot 0.9 \\  
= & 0.7057\cdot 0.9372\cdot 0.9 \\  
= & 0.5952   
= & 0.5952   
\end{align}</math>
\end{align}\,\!</math>
 
<br>


In the case where a switch failure mode causes the standby subsystem to fail, then this mode can be modeled as an individual block in series with the standby subsystem.
In the case where a switch failure mode causes the standby subsystem to fail, then this mode can be modeled as an individual block in series with the standby subsystem.


===Example 6===
===Example===
Consider a car with four new tires and a full-size spare.  Assume the following failure characteristics:
{{:Standby Configuration Example}}
<br>
:• The tires follow a Weibull distribution with a  ..  and an  <math>\eta =</math>  40,000 miles while on the car due to wear.<br>
:• The tires also have a probability of failing due to puncture or other causes.  For this, assume a constant rate for this occurrence with a probability of 1 every 50,000 miles.<br>
:• When not on the car (i.e. is a spare), a tire's probability of failing also has a Weibull distribution with a  <math>\beta =</math>  2 and  <math>\eta =</math>  120,000 miles.<br>
<br>
Assume a mission of 1,000 miles.  If a tire fails during this trip, it will be replaced with the spare.  However, the spare will not be repaired during the trip.  In other words, the trip will continue with the spare on the car and if the spare fails the system will fail.  Determine the probability of system failure.
====Solution to Example 6====
Active failure distribution for tires:
<br>
:• Due to wear-out, Weibull  <math>\beta =4</math>  and  <math>\eta =40,000</math>  miles.<br>
:• Due to random puncture, exponential  <math>\mu =50,000.</math> <br>
:• The quiescent failure distribution is a Weibull distribution with <math>\beta =2</math>  and  <math>\eta =120,000</math>  miles.<br>
<br>
The block diagram for each tire has two blocks in series, one block representing the wear-out mode and the other the random puncture mode, as shown next:
<br>
 
<math></math>
<br>
<br>
[[File:small5.gif|center]]
<br>
 
There are five tires, four active and one standby (represented in the diagram by a standby container with a 4-out-of-5 requirement), as shown next:
<br>
[[File:BStirecontainer.png|center]]
For the standby Wear block, set the active failure and the quiescent distributions, but for the Puncture block, set only the active puncture distribution (because the tire cannot fail due to puncture while stored).  Using BlockSim, the probability of system failure is found to be 0.003 or 0.3%.
<br>


=Note Regarding Numerical Integration Solutions=
<div class="noprint">
{{Examples Box|BlockSim_Examples|<p>More examples on load sharing and standby configurations are available! See also:</p>
{{Examples Link External|http://www.reliasoft.com/BlockSim/examples/rc4/index.htm|Modeling Failure Modes}}<nowiki/>
{{Examples Link|Load_Sharing_Configuration_Example|Load Sharing Configuration Example}}}}
</div>


Load sharing and standby solutions in BlockSim are performed using numerical integration routines. As with any numerical analysis routine, the solution error depends on the number of iterations performed, the step size chosen and related factors, plus the behavior of the underlying function. By default, BlockSim uses a certain set of preset factors. In general, these defaults are sufficient for most problems. If a higher precision or verification of the precision for a specific problem is required, BlockSim's preset options can be modified and/or the integration error can be assessed using the Integration Parameters... option for each container. For more details, you can refer to the documentation on the Algorithm Setup window in the BlockSim 7 User's Guide.
==Note Regarding Numerical Integration Solutions==
Load sharing and standby solutions in BlockSim are performed using numerical integration routines. As with any numerical analysis routine, the solution error depends on the number of iterations performed, the step size chosen and related factors, plus the behavior of the underlying function. By default, BlockSim uses a certain set of preset factors. In general, these defaults are sufficient for most problems. If a higher precision or verification of the precision for a specific problem is required, BlockSim's preset options can be modified and/or the integration error can be assessed using the '''Integration Parameters''' option for each container. For more details, you can refer to the documentation on the Algorithm Setup window in the BlockSim help file.

Latest revision as of 20:09, 8 February 2017

New format available! This reference is now available in a new format that offers faster page load, improved display for calculations and images, more targeted search and the latest content available as a PDF. As of September 2023, this Reliawiki page will not continue to be updated. Please update all links and bookmarks to the latest reference at help.reliasoft.com/reference/system_analysis

Chapter 4: Time-Dependent System Reliability (Analytical)


BlockSimbox.png

Chapter 4  
Time-Dependent System Reliability (Analytical)  

Synthesis-icon.png

Available Software:
BlockSim

Examples icon.png

More Resources:
BlockSim examples

In the RBDs and Analytical System Reliability chapter, different system configuration types were examined, as well as different methods for obtaining the system's reliability function analytically. Because the reliabilities in the problems presented were treated as probabilities (e.g., [math]\displaystyle{ P(A)\,\! }[/math], [math]\displaystyle{ {{R}_{i}}\,\! }[/math] ), the reliability values and equations presented were referred to as static (not time-dependent). Thus, in the prior chapter, the life distributions of the components were not incorporated in the process of calculating the system reliability. In this chapter, time dependency in the reliability function will be introduced. We will develop the models necessary to observe the reliability over the life of the system, instead of at just one point in time. In addition, performance measures such as failure rate, MTTF and warranty time will be estimated for the entire system. The methods of obtaining the reliability function analytically remain identical to the ones presented in the previous chapter, with the exception that the reliabilities will be functions of time. In other words, instead of dealing with [math]\displaystyle{ {{R}_{i}}\,\! }[/math], we will use [math]\displaystyle{ {{R}_{i}}(t)\,\! }[/math]. All examples in this chapter assume that no repairs are performed on the components. Repairable systems analysis will be introduced in a subsequent chapter.

Analytical Life Predictions

The analytical approach presented in the prior chapter involved the determination of a mathematical expression that describes the reliability of the system, expressed in terms of the reliabilities of its components. So far we have estimated only static system reliability (at a fixed time). For example, in the case of a system with three components in series, the system's reliability equation was given by:

[math]\displaystyle{ {{R}_{s}}={{R}_{1}}\cdot {{R}_{2}}\cdot {{R}_{3}} \ \,\! }[/math]

The values of [math]\displaystyle{ {{R}_{1}}\,\! }[/math], [math]\displaystyle{ {{R}_{2}}\,\! }[/math] and [math]\displaystyle{ {{R}_{3}}\,\! }[/math] were given for a common time and the reliability of the system was estimated for that time. However, since the component failure characteristics can be described by distributions, the system reliability is actually time-dependent. In this case, the equation above can be rewritten as:

[math]\displaystyle{ {{R}_{s}}(t)={{R}_{1}}(t)\cdot {{R}_{2}}(t)\cdot {{R}_{3}}(t)\,\! }[/math]

The reliability of the system for any mission time can now be estimated. Assuming a Weibull life distribution for each component, the first equation above can now be expressed in terms of each component's reliability function, or:


[math]\displaystyle{ {{R}_{s}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{2}}} \right)}^{{{\beta }_{2}}}}}}\cdot {{e}^{-{{\left( \tfrac{t}{{{\eta }_{3}}} \right)}^{{{\beta }_{3}}}}}}\,\! }[/math]

In the same manner, any life distribution can be substituted into the system reliability equation. Suppose that the times-to-failure of the first component are described with a Weibull distribution, the times-to-failure of the second component with an exponential distribution and the times-to-failure of the third component with a normal distribution. Then the first equation above can be written as:

[math]\displaystyle{ {{R}_{s}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}}\cdot {{e}^{-{{\lambda }_{2}}t}}\cdot \left[ 1-\Phi \left( \frac{t-{{\mu }_{3}}}{{{\sigma }_{3}}} \right) \right]\,\! }[/math]

It can be seen that the biggest challenge is in obtaining the system's reliability function in terms of component reliabilities, which has already been discussed in depth. Once this has been achieved, calculating the reliability of the system for any mission duration is just a matter of substituting the corresponding component reliability functions into the system reliability equation.


Advantages and Disadvantages

The primary advantage of the analytical solution is that it produces a mathematical expression that describes the reliability of the system. Once the system's reliability function has been determined, other calculations can then be performed to obtain metrics of interest for the system. Such calculations include:

  • Determination of the system's pdf.
  • Determination of warranty periods.
  • Determination of the system's failure rate.
  • Determination of the system's MTTF.

In addition, optimization and reliability allocation techniques can be used to aid engineers in their design improvement efforts. Another advantage in using analytical techniques is the ability to perform static calculations and analyze systems with a mixture of static and time-dependent components. Finally, the reliability importance of components over time can be calculated with this methodology.

The biggest disadvantage of the analytical method is that formulations can become very complicated. The more complicated a system is, the larger and more difficult it will be to analytically formulate an expression for the system's reliability. For particularly detailed systems this process can be quite time-consuming, even with the use of computers. Furthermore, when the maintainability of the system or some of its components must be taken into consideration, analytical solutions become intractable. In these situations, the use of simulation methods may be more advantageous than attempting to develop a solution analytically. Simulation methods are presented in later chapters.

Looking at a Simple "Complex" System Analytically

The complexity involved in an analytical solution can be best illustrated by looking at the simple complex system with 15 components, as shown below.

An RBD of a complex system.

The system reliability for this system (computed using BlockSim) is shown next. The first solution is provided using BlockSim's symbolic solution. In symbolic mode, BlockSim breaks the equation into segments, identified by tokens, that need to be substituted into the final system equation for a complete solution. This creates algebraic solutions that are more compact than if the substitutions were made.

[math]\displaystyle{ \begin{align} {{R}_{System}}= & D2\cdot D3\cdot {{R}_{L}} \\ D3= & +{{R}_{K}}\cdot IK \\ IK= & +{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{O}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{O}}\cdot {{R}_{G}}\cdot {{R}_{F}} \\ & -{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\ & -{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}} \\ & +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}} \\ D2 = & +{{R}_{A}}\cdot {{R}_{E}}\cdot IE \\ IE = & -D1\cdot {{R}_{M}}\cdot {{R}_{N}}+{{R}_{M}}\cdot {{R}_{N}}+D1 \\ D1 = & +{{R}_{D}}\cdot ID \\ ID = & -{{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}} \end{align}\,\! }[/math]

Substituting the terms yields:

[math]\displaystyle{ \begin{align} {{R}_{System}}= & {{R}_{A}}\cdot {{R}_{E}}\cdot {{R}_{L}}\cdot {{R}_{K}} \\ & \cdot \{({{R}_{D}}\cdot {{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}})\cdot {{R}_{M}}\cdot {{R}_{N}} \\ & +{{R}_{M}}\cdot {{R}_{N}}-{{R}_{D}}\cdot {{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}}\} \\ & \cdot \{{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{O}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{O}}\cdot {{R}_{G}}\cdot {{R}_{F}} \\ & -{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\ & -{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}} \\ & +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}}\} \end{align}\,\! }[/math]

BlockSim's automatic algebraic simplification would yield the following format for the above solution:

[math]\displaystyle{ \begin{align} {{R}_{System}}= & (({{R}_{A}}\cdot {{R}_{E}}(-({{R}_{D}}(-{{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}})){{R}_{M}}\cdot {{R}_{N}} \\ & +{{R}_{M}}\cdot {{R}_{N}} \\ & +({{R}_{D}}(-{{R}_{B}}\cdot {{R}_{C}}+{{R}_{B}}+{{R}_{C}})))) \\ & ({{R}_{K}}({{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{O}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\ & -{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{O}}\cdot {{R}_{G}}\cdot {{R}_{F}}-{{R}_{I}}\cdot {{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\ & -{{R}_{I}}\cdot {{R}_{O}}\cdot {{R}_{F}}\cdot {{R}_{H}}-{{R}_{J}}\cdot {{R}_{G}}\cdot {{R}_{F}}\cdot {{R}_{H}} \\ & +RI\cdot {{R}_{O}}\cdot {{R}_{F}} \\ & +{{R}_{I}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{F}}\cdot {{R}_{H}}+{{R}_{J}}\cdot {{R}_{G}})){{R}_{L}}) \ \end{align}\,\! }[/math]

In this equation, each [math]\displaystyle{ {{R}_{i}}\,\! }[/math] represents the reliability function of a block. For example, if [math]\displaystyle{ {{R}_{A}}\,\! }[/math] has a Weibull distribution, then each [math]\displaystyle{ {{R}_{A}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{A}}} \right)}^{{{\beta }_{A}}}}}}\,\! }[/math] and so forth. Substitution of each component's reliability function in the last [math]\displaystyle{ {{R}_{System}}\,\! }[/math] equation above will result in an analytical expression for the system reliability as a function of time, or [math]\displaystyle{ {{R}_{s}}(t)\,\! }[/math], which is the same as [math]\displaystyle{ (1-cd{{f}_{System}}).\,\! }[/math]

Obtaining Other Functions of Interest

Once the system reliability equation (or the cumulative density function, cdf) has been determined, other functions and metrics of interest can be derived.

Consider the following simple system:

Simple two-component system.

Furthermore, assume that component 1 follows an exponential distribution with a mean of 10,000 ([math]\displaystyle{ \mu =10,000,\,\! }[/math] [math]\displaystyle{ \lambda =1/10,000)\,\! }[/math] and component 2 follows a Weibull distribution with [math]\displaystyle{ \beta =6\,\! }[/math] and [math]\displaystyle{ \eta =10,000\,\! }[/math]. The reliability equation of this system is:

[math]\displaystyle{ \begin{align} {{R}_{S}}(t)= & {{R}_{1}}(t)\cdot {{R}_{2}}(t) \\ = & {{e}^{-\lambda t}}\cdot {{e}^{-{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \\ = & {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \ \end{align}\,\! }[/math]

The system cdf is:

[math]\displaystyle{ \begin{align} {{F}_{S}}(t)= & 1-({{R}_{1}}(t)\cdot {{R}_{2}}(t)) \\ = & 1-\left( {{e}^{-\lambda t}}\cdot {{e}^{-{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \right) \\ = & 1-\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right) \end{align}\,\! }[/math]

System pdf

Once the equation for the reliability of the system has been obtained, the system's pdf can be determined. The pdf is the derivative of the reliability function with respect to time or:

[math]\displaystyle{ {{f}_{s}}(t)=-\frac{d({{R}_{s}}(t))}{dt} \ \,\! }[/math]

For the system shown above, this is:

[math]\displaystyle{ \begin{align} {{f}_{s}}(t)= & -\frac{d}{dt}\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right) \\ = & -\frac{d}{dt}\left( {{e}^{-\tfrac{1}{10,000}t}} \right)\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \\ & +{{e}^{-\tfrac{1}{10,000}t}}\left[ -\frac{d}{dt}\left( {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right) \right] \\ = & {{f}_{1}}(t)\cdot {{R}_{2}}(t)+{{f}_{2}}(t)\cdot {{R}_{1}}(t) \end{align}\,\! }[/math]

The next figure shows a plot of the pdf equation.

pdf plot of the two-component system

Conditional Reliability

Conditional reliability is the probability of a system successfully completing another mission following the successful completion of a previous mission. The time of the previous mission and the time for the mission to be undertaken must be taken into account for conditional reliability calculations. The system's conditional reliability function is given by:

[math]\displaystyle{ R(t|T)=\frac{R(T+t)}{R(T)}\,\! }[/math]

Equation above gives the reliability for a new mission of duration [math]\displaystyle{ t\,\! }[/math] having already accumulated [math]\displaystyle{ T\,\! }[/math] hours of operation up to the start of this new mission. The system is evaluated to assure that it will start the next mission successfully.

For the simple two-component system, the reliability for mission of [math]\displaystyle{ t\,\! }[/math] = 1000 hours, having an age of [math]\displaystyle{ T\,\! }[/math] = 500 hours, is:

[math]\displaystyle{ \begin{align} {{R}_{S}}(t=1000|T=500)= & \frac{R(T+t)}{R(T)} \\ = & \frac{R(1500)}{R(500)} \\ = & \frac{{{e}^{-\tfrac{1500}{10,000}}}\cdot {{e}^{-{{\left( \tfrac{1500}{10,000} \right)}^{6}}}}}{{{e}^{-\tfrac{500}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{500}{10,000} \right)}^{6}}}}} \\ = & 0.9048=90.48% \end{align}\,\! }[/math]

Conditional Reliability for Components

Now in this formulation, it was assumed that the accumulated age was equivalent for both units. That is, both started life at zero and aged to 500. It is possible to consider an individual component that has already accumulated some age (used component) in the same formulation. To illustrate this, assume that component 2 started life with an age of T = 100. Then the reliability equation of the system, as given in [math]\displaystyle{ {{R}_{S}}(t)= {{R}_{1}}(t)\cdot {{R}_{2}}(t)\,\! }[/math], would need to be modified to include a conditional term for 2, or:

[math]\displaystyle{ {{R}_{S}}(t)={{R}_{1}}(t)\cdot \frac{{{R}_{2}}({{T}_{2}}+t)}{{{R}_{2}}({{T}_{2}})} \ \,\! }[/math]

In BlockSim, the start age input box may be used to specify a starting age greater than zero.

System Failure Rate

Once the distribution of the system has been determined, the failure rate can also be obtained by dividing the pdf by the reliability function:

[math]\displaystyle{ {{\lambda }_{s}}\left( t \right)=\frac{{{f}_{s}}\left( t \right)}{{{R}_{s}}\left( t \right)} \ \,\! }[/math]

For the simple two-component system:

[math]\displaystyle{ \begin{align} {{\lambda }_{s}}\left( t \right)= & \frac{-\tfrac{d}{dt}\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)}{{{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}} \\ = & \frac{-\tfrac{d}{dt}\left( {{e}^{-\tfrac{1}{10,000}t}} \right)}{{{e}^{-\tfrac{1}{10,000}t}}}+\frac{-\tfrac{d}{dt}\left( {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)}{{{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}} \\ = & \frac{{{f}_{1}}}{{{R}_{1}}}+\frac{{{f}_{2}}}{{{R}_{2}}} \\ = & {{\lambda }_{1}}+{{\lambda }_{2}} \ \end{align}\,\! }[/math]

The following figure shows a plot of the equation.

Failure rate function plot of the two component system

BlockSim uses numerical methods to estimate the failure rate. It should be pointed out that as [math]\displaystyle{ t\to \infty \,\! }[/math], numerical evaluation of the first equation above is constrained by machine numerical precision. That is, there are limits as to how large [math]\displaystyle{ t\,\! }[/math] can get before floating point problems arise. For example, at [math]\displaystyle{ t=5,000,000\,\! }[/math] both numerator and denominator will tend to zero (e.g., [math]\displaystyle{ {{e}^{-\tfrac{5,000,000}{10,000}}}=7.1245\times {{10}^{-218}}\,\! }[/math] ). As these numbers become very small they will start looking like a zero to the computer, or cause a floating point error, resulting in a [math]\displaystyle{ \tfrac{0}{0}\,\! }[/math] or [math]\displaystyle{ \tfrac{X}{0}\,\! }[/math] operation. In these cases, BlockSim will return a value of "[math]\displaystyle{ N/A\,\! }[/math]" for the result. Obviously, this does not create any practical constraints.

System Mean Life (Mean Time To Failure)

The mean life (or mean time to failure, MTTF) can be obtained by integrating the system reliability function from zero to infinity:

[math]\displaystyle{ MTTF=\int_{0}^{\infty }{{R}_{s}}\left( t \right)dt \ \,\! }[/math]

The mean time is a performance index and does not provide any information about the behavior of the failure distribution of the system.

For the simple two-component system:

[math]\displaystyle{ \begin{align} MTTF= & \int_{0}^{\infty }\left( {{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}} \right)dt \\ = & 5978.9 \end{align}\,\! }[/math]

Warranty Period and BX Life

Sometimes it is desirable to know the time value associated with a certain reliability. Warranty periods are often calculated by determining what percentage of the failure population can be covered financially and estimating the time at which this portion of the population will fail. Similarly, engineering specifications may call for a certain BX life, which also represents a time period during which a certain proportion of the population will fail. For example, the B10 life is the time in which 10% of the population will fail. This is obtained by setting [math]\displaystyle{ {{R}_{S}}(t)\,\! }[/math] to the desired value and solving for [math]\displaystyle{ t.\,\! }[/math]

For the simple two-component system:

[math]\displaystyle{ {{R}_{s}}\left( t \right)={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}\,\! }[/math]

To compute the time by which reliability would be equal to 90%, equation above is recast as follows and solved for [math]\displaystyle{ t.\,\! }[/math]

[math]\displaystyle{ 0.90={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}\,\! }[/math]

In this case, [math]\displaystyle{ t=1053.59\,\! }[/math]. Equivalently, the B10 life for this system is also 1053.59. Except for some trivial cases, a closed form solution for [math]\displaystyle{ t\,\! }[/math] cannot be obtained. Thus, it is necessary to solve for [math]\displaystyle{ t\,\! }[/math] using numerical methods. BlockSim uses numerical methods.


Examples

Components in Series


Consider a system consisting of 3 exponential units, connected in series, with the following failure rates (in failures per hour): [math]\displaystyle{ {{\lambda }_{1}} = 0.0002\,\! }[/math], [math]\displaystyle{ {{\lambda }_{2}} = 0.0005\,\! }[/math] and [math]\displaystyle{ {{\lambda }_{3}} = 0.0001\,\! }[/math].

  • Obtain the reliability equation for the system.
  • What is the reliability of the system after 150 hours of operation?
  • Obtain the system's pdf.
  • Obtain the system's failure rate equation.
  • What is the MTTF for the system?
  • What should the warranty period be for a 90% reliability?


Solution

The analytical expression for the reliability of the system is given by:

[math]\displaystyle{ \begin{align} {{R}_{s}}(t)= & {{R}_{1}}(t)\cdot {{R}_{2}}(t)\cdot {{R}_{3}}(t) \\ = & {{e}^{-{{\lambda }_{1}}t}}\cdot {{e}^{-{{\lambda }_{2}}t}}\cdot {{e}^{-{{\lambda }_{1}}t}} \\ = & {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}} \end{align}\,\! }[/math]

At 150 hours of operation, the reliability of the system is:

[math]\displaystyle{ \begin{align} {{R}_{s}}(t)= & {{e}^{-(0.0002+0.0005+0.0001)150}} \\ = & 0.8869\text{ or }88.69% \end{align}\,\! }[/math]

In order to obtain the system's pdf, the derivative of the reliability equation given in the first equation above is taken with respect to time, or:

[math]\displaystyle{ \begin{align} {{f}_{s}}(t)= & -\frac{d[{{R}_{s}}(t)]}{dt} \\ = & -\frac{d\left[ {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}} \right]}{dt} \\ = & ({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})\cdot {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}} \end{align}\,\! }[/math]


The system's failure rate can now be obtained simply by dividing the system's pdf given in the equation above by the system's reliability function given in the first equation above, and:

[math]\displaystyle{ \begin{align} {{\lambda }_{s}}\left( t \right)= & \frac{{{f}_{s}}\left( t \right)}{{{R}_{s}}\left( t \right)} \\ = & \frac{({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})\cdot {{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}}}{{{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}}} \\ = & ({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}}) \\ = & 0.0008\text{ }fr/hr \end{align}\,\! }[/math]


Combining [math]\displaystyle{ MTTF=\int_{0}^{\infty }{{R}_{s}}\left( t \right)dt \,\! }[/math] and the first equation above, the system's MTTF can be obtained:

[math]\displaystyle{ \begin{align} MTTF= & \int_{0}^{\infty }{{R}_{s}}\left( t \right)dt \\ = & \int_{0}^{\infty }{{e}^{-({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})t}}dt \\ = & \frac{1}{({{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}})} \\ = & 1250\text{ }hr \end{align}\,\! }[/math]


Solving the first equation above with respect to time will yield the corresponding warranty period for a 90% reliability. In this case, the system reliability equation is simple and a closed form solution exists. The warranty time can now be found by solving:

[math]\displaystyle{ \begin{align} t= & -\frac{\ln (R)}{{{\lambda }_{1}}+{{\lambda }_{2}}+{{\lambda }_{3}}} \\ = & -\frac{\ln (0.9)}{0.0008} \\ = & 131.7\text{ }hr \end{align}\,\! }[/math]


Thus, the warranty period should be 132 hours.


Components in Parallel


Consider the system shown next.

Complex bridge system in Example 2.

Components [math]\displaystyle{ A\,\! }[/math] through [math]\displaystyle{ E\,\! }[/math] are Weibull distributed with [math]\displaystyle{ \beta =1.2\,\! }[/math] and [math]\displaystyle{ \eta =1230\,\! }[/math] hours. The starting and ending blocks cannot fail.

Determine the following:

  • The reliability equation for the system and its corresponding plot.
  • The system's pdf and its corresponding plot.
  • The system's failure rate equation and the corresponding plot.
  • The MTTF.
  • The warranty time for a 90% reliability.
  • The reliability for a 200-hour mission, if it is known that the system has already successfully operated for 200 hours.


Solution

The first step is to obtain the reliability function for the system. The methods described in the RBDs and Analytical System Reliability chapter can be employed, such as the event space or path-tracing methods. Using BlockSim, the following reliability equation is obtained:

[math]\displaystyle{ \begin{align} {{R}_{s}}(t)= & ({{R}_{Start}}\cdot {{R}_{End}}(2{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\ & -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}-{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\ & -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{B}}\cdot {{R}_{E}}-{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\ & -{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}}+{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\ & +{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}+{{R}_{A}}\cdot {{R}_{D}}+{{R}_{B}}\cdot {{R}_{E}})) \end{align}\,\! }[/math]


Note that since the starting and ending blocks cannot fail, [math]\displaystyle{ {{R}_{Start}}=1\,\! }[/math] and [math]\displaystyle{ {{R}_{End}}=1,\,\! }[/math] the equation above can be reduced to:

[math]\displaystyle{ \begin{align} {{R}_{s}}(t)= & 2\cdot {{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\ & -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}-{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\ & -{{R}_{A}}\cdot {{R}_{D}}\cdot {{R}_{B}}\cdot {{R}_{E}}-{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}} \\ & -{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}\cdot {{R}_{E}}+{{R}_{A}}\cdot {{R}_{C}}\cdot {{R}_{E}} \\ & +{{R}_{D}}\cdot {{R}_{C}}\cdot {{R}_{B}}+{{R}_{A}}\cdot {{R}_{D}}+{{R}_{B}}\cdot {{R}_{E}} \end{align}\,\! }[/math]


where [math]\displaystyle{ {{R}_{A}}\,\! }[/math] is the reliability equation for Component A, or:

[math]\displaystyle{ {{R}_{A}}(t)={{e}^{-{{\left( \tfrac{t}{{{\eta }_{A}}} \right)}^{{{\beta }_{A}}}}}}\,\! }[/math]


[math]\displaystyle{ {{R}_{B}}\,\! }[/math] is the reliability equation for Component [math]\displaystyle{ B\,\! }[/math], etc.


Since the components in this example are identical, the system reliability equation can be further reduced to:

[math]\displaystyle{ \begin{align} {{R}_{s}}(t)=2R{{(t)}^{2}}+2R{{(t)}^{3}}-5R{{(t)}^{4}}+2R{{(t)}^{5}} \end{align}\,\! }[/math]

Or, in terms of the failure distribution:

[math]\displaystyle{ {{R}_{s}}(t)=2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}\,\! }[/math]

The corresponding plot is given in the following figure.

Reliability plot for the system.

In order to obtain the system's pdf, the derivative of the reliability equation given above is taken with respect to time, resulting in:

[math]\displaystyle{ \begin{align} {{f}_{s}}(t)= & 4\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+6\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \\ & -20\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+10\cdot \frac{\beta }{\eta }{{\left( \frac{t}{\eta } \right)}^{\beta -1}}{{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}} \end{align}\,\! }[/math]

The pdf can now be plotted for different time values, [math]\displaystyle{ t\,\! }[/math], as shown in the following figure.

pdf plot for the system.

The system's failure rate can be obtained by dividing the system's pdf, given in equation above, by the system's reliability function given in

[math]\displaystyle{ \begin{align} {{R}_{s}}(t)= & 2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}\end{align}\,\! }[/math], or:
[math]\displaystyle{ \begin{align} {{\lambda }_{s}}(t)= & \frac{4\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+6\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}}{2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}} \\ & +\frac{-20\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+10\cdot \tfrac{\beta }{\eta }{{\left( \tfrac{t}{\eta } \right)}^{\beta -1}}{{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}}{2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}} \end{align}\,\! }[/math]

The corresponding plot is given below.

Failure rate for the system.

The [math]\displaystyle{ MTTF\,\! }[/math] of the system is obtained by integrating the system's reliability function given by [math]\displaystyle{ {{R}_{s}}(t)=2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}\,\! }[/math] from time zero to infinity, as given by [math]\displaystyle{ MTTF=\int_{0}^{\infty }{{R}_{s}}\left( t \right)dt \ \,\! }[/math]. Using BlockSim's Analytical QCP, an [math]\displaystyle{ MTTF\,\! }[/math] of 1007.8 hours is calculated, as shown in the figure below.

MTTF of the system.


The warranty time can be obtained by solving [math]\displaystyle{ {{R}_{s}}(t)\,\! }[/math] with respect to time for a system reliability [math]\displaystyle{ {{R}_{s}}=0.9\,\! }[/math]. Using the Analytical QCP and selecting the Reliable Life option, a time of 372.72 hours is obtained, as shown in the following figure.

Time at which R=0.9 or 90% for the system.


Lastly, the conditional reliability can be obtained using [math]\displaystyle{ R(t|T)=\frac{R(T+t)}{R(T)}\,\! }[/math] and [math]\displaystyle{ {{R}_{s}}(t)=2\cdot {{e}^{-2{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-3{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}-5\cdot {{e}^{-4{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}+2\cdot {{e}^{-5{{\left( \tfrac{t}{\eta } \right)}^{\beta }}}}\,\! }[/math], or:

[math]\displaystyle{ \begin{align} R(200|200)= & \frac{R(400)}{R(200)} \\ = & \frac{0.883825}{0.975321} \\ = & 0.906189 \end{align}\,\! }[/math]

This can be calculated using BlockSim's Analytical QCP, as shown below.

Conditional reliability calculation for the system.

Approximating the System cdf

In many cases, it is valuable to fit a distribution that represents the system's times-to-failure. This can be useful when the system is part of a larger assembly and may be used for repeated calculations or in calculations for other systems. In cases such as this, it can be useful to characterize the system's behavior by fitting a distribution to the overall system and calculating parameters for this distribution. This is equivalent to fitting a single distribution to describe [math]\displaystyle{ {{R}_{S}}(t)\,\! }[/math]. In essence, it is like reducing the entire system to a component in order to simplify calculations.

For the system shown below:

Simple two-component system.
[math]\displaystyle{ {{R}_{S}}(t)={{e}^{-\tfrac{1}{10,000}t}}\cdot {{e}^{-{{\left( \tfrac{t}{10,000} \right)}^{6}}}}\,\! }[/math]


To compute an approximate reliability function for this system, [math]\displaystyle{ {{R}_{A}}(t)\simeq {{R}_{S}}(t)\,\! }[/math], one would compute [math]\displaystyle{ n\,\! }[/math] pairs of reliability and time values and then fit a single distribution to the data, or:

[math]\displaystyle{ \begin{align} {{R}_{S}}(t= & 10,396.7)=10% \\ {{R}_{S}}(t= & 9,361.9)=20% \\ & ... \\ {{R}_{S}}(t= & 1,053.6)=90% \end{align}\,\! }[/math]

A single distribution, [math]\displaystyle{ {{R}_{A}}(t)\,\! }[/math], that approximates [math]\displaystyle{ {{R}_{S}}(t)\,\! }[/math] can now be computed from these pairs using life data analysis methods. If using the Weibull++ software, one would enter the values as free-form data.

Example

Compute a single Weibull distribution approximation for the parallel system in the previous example.

Solution

Assume that the system can be approximated by a 2-parameter Weibull distribution with [math]\displaystyle{ \beta =2.02109\,\! }[/math] and [math]\displaystyle{ \eta =1123.51\,\! }[/math]. In BlockSim, this is accomplished by representing the entire system as one distribution. To do this, click the Distribution fit area on the diagram's control panel, as shown next.


Representing a system with a distribution.


This opens the Distribution Estimator window, which allows you to select a distribution to represent the data.

Distribution Estimator.

When you analyze the diagram, BlockSim will generate a number of system failure times based on the system's reliability function. The system's reliability function can be used to solve for a time value associated with an unreliability value. For example, consider an unreliability value of [math]\displaystyle{ F(t)=0.11\,\! }[/math]. Using the system's reliability equation, the corresponding time-to-failure for a 0.11 unreliability is 389.786 hours. The distribution of the generated time values can then be fitted to a probability distribution function.

When enough points have been generated, the selected distribution will be fitted to the generated data set and the distribution's parameters will be calculated. In addition, if ReliaSoft's Weibull++ is installed, the fit of the distribution can be analyzed using a Weibull++ standard folio, as shown in the next figure. It is recommended that the analyst examine the fit to ascertain the applicability of the approximation.


Using Weibull++ to calculate distribution parameters.

Duty Cycle

Components of a system may not operate continuously during a system's mission, or may be subjected to loads greater or lesser than the rated loads during system operation. To model this, a factor called the Duty Cycle ( [math]\displaystyle{ {{d}_{c}}\,\! }[/math] ) is used. The duty cycle may also be used to account for changes in environmental stress, such as temperature changes, that may effect the operation of a component. The duty cycle is a positive value, with a default value of 1 representing continuous operation at rated load, and any values other than 1 representing other load values with respect to the rated load value (or total operating time). A duty cycle value higher than 1 indicates a load in excess of the rated value. A duty cycle value lower than 1 indicates that the component is operating at a load lower than the rated load or not operating continuously during the system's mission. For instance, a duty cycle of 0.5 may be used for a component that operates only half of the time during the system's mission.

The reliability metrics for a component with a duty cycle are calculated as follows. Let [math]\displaystyle{ {{d}_{c}}\,\! }[/math] represent the duty cycle during a particular mission of the component, [math]\displaystyle{ t\,\! }[/math] represent the mission time and [math]\displaystyle{ {t}'\,\! }[/math] represent the accumulated age. Then:

[math]\displaystyle{ {t}'={{d}_{c}}\times t\,\! }[/math]

The reliability equation for the component is:

[math]\displaystyle{ R({t}')=R({{d}_{c}}\times t)\,\! }[/math]

The component pdf is:

[math]\displaystyle{ f({t}')=-\frac{d(R({t}'))}{dt}=-\frac{d(R({{d}_{c}}\times t))}{dt}={{d}_{c}}f({{d}_{c}}\times t)\,\! }[/math]

The failure rate of the component is:

[math]\displaystyle{ \lambda ({t}')=\frac{f({t}')}{R({t}')}=\frac{{{d}_{c}}f({{d}_{c}}\times t)}{R({{d}_{c}}\times t)}={{d}_{c}}\lambda ({{d}_{c}}\times t)\,\! }[/math]

Example

Consider a computer system with three components: a processor, a hard drive and a CD drive in series as shown next. Assume that all three components follow a Weibull failure distribution. The processor has the following parameters, [math]\displaystyle{ {{\beta }_{1}}=1.5\,\! }[/math] and [math]\displaystyle{ {{\eta }_{1}}=5000\,\! }[/math]. For the hard drive, the parameters are [math]\displaystyle{ {{\beta }_{2}}=2.5\,\! }[/math] and [math]\displaystyle{ {{\eta }_{2}}=3000\,\! }[/math], and for the CD drive they are [math]\displaystyle{ {{\beta }_{3}}=2\,\! }[/math] and [math]\displaystyle{ {{\eta }_{3}}=4000\,\! }[/math]. Determine the reliability of the computer system after one year (365 days) of operation, assuming that the CD drive is used only 30% of the time.


BS5ex4 new.png

Solution

The reliability of the processor after 365 days of operation is given by:


[math]\displaystyle{ \begin{align} {{R}_{processor}}(365)= & {{e}^{-{{\left( \tfrac{365}{{{\eta }_{1}}} \right)}^{{{\beta }_{1}}}}}} \\ = & {{e}^{-{{\left( \tfrac{365}{5000} \right)}^{1.5}}}} \\ = & 0.9805\text{ or }98.05% \end{align}\,\! }[/math]


The reliability of the hard drive after 365 days of operation is given by:


[math]\displaystyle{ \begin{align} {{R}_{harddrive}}(365)= & {{e}^{-{{\left( \tfrac{365}{{{\eta }_{2}}} \right)}^{{{\beta }_{2}}}}}} \\ = & {{e}^{-{{\left( \tfrac{365}{3000} \right)}^{2.5}}}} \\ = & 0.9948\text{ or }99.48% \end{align}\,\! }[/math]


The reliability of the CD drive after 365 days of operation (taking into account the 30% operation using a duty cycle of 0.3) is given by:


[math]\displaystyle{ \begin{align} {{R}_{CDdrive}}(365)= & {{e}^{-{{\left( \tfrac{{{d}_{c}}\times 365}{{{\eta }_{3}}} \right)}^{{{\beta }_{3}}}}}} \\ = & {{e}^{-{{\left( \tfrac{0.3\times 365}{4000} \right)}^{2}}}} \\ = & 0.9993\text{ or }99.93% \end{align}\,\! }[/math]


Thus the reliability of the computer system after 365 days of operation is:


[math]\displaystyle{ \begin{align} {{R}_{s}}(365)= & {{R}_{processor}}(365)\cdot {{R}_{harddrive}}(365)\cdot {{R}_{CDdrive}}(365) \\ = & 0.9805\cdot 0.9948\cdot 0.9993 \\ = & 0.9747\text{ or }97.47% \end{align}\,\! }[/math]


The result can also be obtained in BlockSim by creating an RBD of the system and using the Quick Calculation Pad (QCP) to calculate the reliability, as shown in the following figure.


Result for the computer system reliability.

Load Sharing

As presented in earlier chapters, a reliability block diagram (RBD) allows you to graphically represent how the components within a system are reliability-wise connected. In most cases, independence is assumed across the components within the system. For example, the failure of component A does not affect the failure of component B. However, if a system consists of components that are sharing a load, then the assumption of independence no longer holds true.

If one component fails, then the component(s) that are still operating will have to assume the failed unit's portion of the load. Therefore, the reliabilities of the surviving unit(s) will change. Calculating the system reliability is no longer an easy proposition. In the case of load sharing components, the change of the failure distributions of the surviving components must be known in order to determine the system's reliability.

To illustrate this, consider a system of two units connected reliability-wise in parallel as shown below.

Two units connected reliability-wise in parallel.

Assume that the units must supply an output of 8 volts and that if both units are operational, each unit is to supply 50% of the total output. If one of the units fails, then the surviving unit supplies 100%. Furthermore, assume that having to supply the entire load has a negative impact on the reliability characteristics of the surviving unit.

Because the reliability characteristics of the unit change based on the load it is sharing, a method that can model the effect of the load on life should be used. One way to do this is to use a life distribution along with a life-stress relationship (as discussed in A Brief Introduction to Life-Stress Relationships) for each component. The detailed discussion for this method can be found at Additional Information on Load Sharing. Another simple way is to use the concept of acceleration factors and assume that the load has a linear effect on the failure time. If the load is doubled, then the life of the component will be shortened by half.

For the above load sharing system, the reliability of each component is a function of time and load. For example, for Unit 1, the reliability and the probability density function are:

[math]\displaystyle{ {{R}_{1}}(t,{{S}_{1}})\,\!\,\! }[/math] and [math]\displaystyle{ {{f}_{1}}(t,{{S}_{1}})\,\!\,\! }[/math]

where [math]\displaystyle{ {{S}_{1}}\,\!\,\! }[/math] is the load shared by Unit 1 at time t and the total load of the system is [math]\displaystyle{ S={{S}_{1}}+{{S}_{2}}\,\!\,\! }[/math]. At the beginning, both units are working. Assume that Unit 1 fails at time x and Unit 2 takes over the entire load. The reliability for Unit 2 at time x is:

[math]\displaystyle{ \begin{align} & {{R}_{2}}(x,{{S}_{2}})={{R}_{2}}({{t}_{2e}},S) \\ & {{t}_{2e}}=\frac{{{S}_{2}}}{S}x \\ \end{align}\,\!\,\! }[/math]

[math]\displaystyle{ {{t}_{2e}}\,\!\,\! }[/math] is the equivalent time for Unit 2 at time x if it is operated with load S. The equivalent time concept is illustrated in the following plot.

Illustrating [math]\displaystyle{ {{t}_{e}}\,\! }[/math]

The system reliability at time t is:

[math]\displaystyle{ \begin{align} & R(t)={{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}}) \\ & +\int_{0}^{t}{{{f}_{1}}(x,{{S}_{1}})\cdot {{R}_{2}}({{t}_{2e}}+(t-x),S)dx} \\ & +\int_{0}^{t}{{{f}_{2}}(x,{{S}_{2}})\cdot {{R}_{1}}({{t}_{1e}}+(t-x),S)dx} \end{align}\,\!\,\! }[/math]

In BlockSim, the failure time distribution for each component is defined at the load of S. The reliability function for a component at a given load is calculated as:

[math]\displaystyle{ {{R}_{i}}(t,{{S}_{i}})={{R}_{i}}(t\times \frac{{{S}_{i}}}{S},S)\,\!\,\! }[/math]
[math]\displaystyle{ {{f}_{i}}(t,{{S}_{i}})=\frac{{{S}_{i}}}{S}{{f}_{i}}\left( t\times \frac{{{S}_{i}}}{S},S \right)\,\!\,\! }[/math]

From the above equation, it can be seen that the concept used in the calculation for load sharing is the same as the concept used in the calculation for duty cycle.

Example

In the following load sharing system, Block 1 follows a Weibull failure distribution with[math]\displaystyle{ {{\beta }_{1}}=1.5\,\! }[/math], and[math]\displaystyle{ {{\eta }_{1}}=1,000\,\! }[/math] . Block 2 follows a Weibull failure distribution with[math]\displaystyle{ {{\beta }_{2}}=2\,\! }[/math], and[math]\displaystyle{ {{\eta }_{2}}=2,000\,\! }[/math]. The load for Block 1 is 1 unit, and for Block 2 it is 3 units. Calculate the system reliability at time 1,500.

Load Sharing System

Block 1 shares 25% (P1) of the entire load, and Block 2 shares 75% (P2) of it. Therefore, we have the following equations for calculating the system reliability:

[math]\displaystyle{ {{R}_{i}}(t,{{S}_{i}})={{R}_{i}}({{P}_{i}}\times t)\,\! }[/math],
[math]\displaystyle{ {{f}_{i}}(x,{{S}_{i}})={{P}_{i}}{{f}_{i}}({{P}_{i}}\times x)\,\! }[/math]

and:

[math]\displaystyle{ {{R}_{1}}({{t}_{1e}}+(t-x))={{R}_{1}}({{P}_{1}}x+t-x)={{R}_{1}}(t-{{P}_{2}}x)\,\! }[/math]
[math]\displaystyle{ {{R}_{2}}({{t}_{2e}}+(t-x))={{R}_{2}}({{P}_{2}}x+t-x)={{R}_{2}}(t-{{P}_{1}}x)\,\! }[/math]

Using the above equations in the system reliability function, we get:

[math]\displaystyle{ \begin{align} & R(t)={{R}_{1}}(t,{{S}_{1}})\cdot {{R}_{2}}(t,{{S}_{2}})+\int_{0}^{t}{{{f}_{1}}(x,{{S}_{1}})\cdot {{R}_{2}}({{t}_{2e}}+(t-x),S)dx} \\ & +\int_{0}^{t}{{{f}_{2}}(x,{{S}_{2}})\cdot {{R}_{1}}({{t}_{1e}}+(t-x),S)dx} \\ & ={{R}_{1}}({{P}_{1}}\times t)\cdot {{R}_{2}}({{P}_{2}}\times t)+\int_{0}^{t}{{{P}_{1}}{{f}_{1}}({{P}_{1}}x)\cdot {{R}_{2}}(t-{{P}_{1}}x)dx} \\ & +\int_{0}^{t}{{{P}_{2}}{{f}_{2}}({{P}_{2}}x)\cdot {{R}_{1}}(t-{{P}_{2}}x)dx} \end{align}\,\! }[/math]

The calculated system reliability at time 1,500 is 0.8569, as given below.

Calculated System Reliability at Time 1,500

Standby Components

In the previous section, the case of a system with load sharing components was presented. This is a form of redundancy with dependent components. That is, the failure of one component affects the failure of the other(s). This section presents another form of redundancy: standby redundancy. In standby redundancy the redundant components are set to be under a lighter load condition (or no load) while not needed and under the operating load when they are activated.

In standby redundancy the components are set to have two states: an active state and a standby state. Components in standby redundancy have two failure distributions, one for each state. When in the standby state, they have a quiescent (or dormant) failure distribution and when operating, they have an active failure distribution.

In the case that both quiescent and active failure distributions are the same, the units are in a simple parallel configuration (also called a hot standby configuration). When the rate of failure of the standby component is lower in quiescent mode than in active mode, that is called a warm standby configuration. When the rate of failure of the standby component is zero in quiescent mode (i.e., the component cannot fail when in standby), that is called a cold standby configuration.

Simple Standby Configuration

Consider two components in a standby configuration. Component 1 is the active component with a Weibull failure distribution with parameters [math]\displaystyle{ \beta = 1.5\,\! }[/math] and [math]\displaystyle{ \eta = 1,000 \,\! }[/math]. Component 2 is the standby component. When Component 2 is operating, it also has a Weibull failure distribution with [math]\displaystyle{ \beta = 1.5\,\! }[/math] and [math]\displaystyle{ \eta = 1,000 \,\! }[/math]. Furthermore, assume the following cases for the quiescent distribution.

  • Case 1: The quiescent distribution is the same as the active distribution (hot standby).
  • Case 2: The quiescent distribution is a Weibull distribution with [math]\displaystyle{ \beta = 1.5\,\! }[/math] and [math]\displaystyle{ \eta = 2000\,\! }[/math] (warm standby).
  • Case 3: The component cannot fail in quiescent mode (cold standby).

In this case, the reliability of the system at some time, [math]\displaystyle{ t\,\! }[/math], can be obtained using the following equation:

[math]\displaystyle{ R(t)={{R}_{1}}(t)+\underset{0}{\overset{t}{\mathop \int }}\,{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x)\cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}dx \ \,\! }[/math]

where:

  • [math]\displaystyle{ {{R}_{1}}\,\! }[/math] is the reliability of the active component.
  • [math]\displaystyle{ {{f}_{1}}\,\! }[/math] is the pdf of the active component.
  • [math]\displaystyle{ {{R}_{2;SB}}\,\! }[/math] is the reliability of the standby component when in standby mode (quiescent reliability).
  • [math]\displaystyle{ {{R}_{2;A}}\,\! }[/math] is the reliability of the standby component when in active mode.
  • [math]\displaystyle{ {{t}_{e}}\,\! }[/math] is the equivalent operating time for the standby unit if it had been operating at an active mode, such that:
[math]\displaystyle{ \begin{align} {{R}_{2;SB}}(x)={{R}_{2;A}}({{t}_{e}}) \end{align}\,\! }[/math]

The second equation above can be solved for [math]\displaystyle{ {{t}_{e}}\,\! }[/math] and substituted into the first equation above. The following figure illustrates the example as entered in BlockSim using a standby container.

Standby container.

The active and standby blocks are within a container, which is used to specify standby redundancy. Since the standby component has two distributions (active and quiescent), the Block Properties window of the standby block has two pages for specifying each one. The following figures illustrate these pages.

Defining the active failure distribution


Defining the quiescent failure distribution

The system reliability results for 1000 hours are given in the following table:

5-24.png

Note that even though the [math]\displaystyle{ \beta \,\! }[/math] value for the quiescent distribution is the same as in the active distribution, it is possible that the two can be different. That is, the failure modes present during the quiescent mode could be different from the modes present during the active mode. In that sense, the two distribution types can be different as well (e.g., lognormal when quiescent and Weibull when active).

In many cases when considering standby systems, a switching device may also be present that switches from the failed active component to the standby component. The reliability of the switch can also be incorporated into [math]\displaystyle{ R(t)={{R}_{1}}(t)+\underset{0}{\overset{t}{\mathop \int }}\,{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x)\cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}dx \ \,\! }[/math]

as presented in the next section.

BlockSim's System Reliability Equation window returns a single token for the reliability of units in a standby configuration. This is the same as the load sharing case presented in the previous section.

Reliability of Standby Systems with a Switching Device

In many cases when dealing with standby systems, a switching device is present that will switch to the standby component when the active component fails. Therefore, the failure properties of the switch must also be included in the analysis.

BS5 26 2 new.png

In most cases when the reliability of a switch is to be included in the analysis, two probabilities can be considered. The first and most common one is the probability of the switch performing the action (i.e., switching) when requested to do so. This is called Switch Probability per Request in BlockSim and is expressed as a static probability (e.g., 90%). The second probability is the quiescent reliability of the switch. This is the reliability of the switch as it ages (e.g., the switch might wear out with age due to corrosion, material degradation, etc.). Thus it is possible for the switch to fail before the active component fails. However, a switch failure does not cause the system to fail, but rather causes the system to fail only if the switch is needed and the switch has failed. For example, if the active component does not fail until the mission end time and the switch fails, then the system does not fail. However, if the active component fails and the switch has also failed, then the system cannot be switched to the standby component and it therefore fails.

In analyzing standby components with a switching device, either or both failure probabilities (during the switching or while waiting to switch) can be considered for the switch, since each probability can represent different failure modes. For example, the switch probability per request may represent software-related issues or the probability of detecting the failure of an active component, and the quiescent probability may represent wear-out type failures of the switch.

To illustrate the formulation, consider the previous example that assumes perfect switching. To examine the effects of including an imperfect switch, assume that when the active component fails there is a 90% probability that the switch will switch from the active component to the standby component. In addition, assume that the switch can also fail due to a wear-out failure mode described by a Weibull distribution with [math]\displaystyle{ \beta = 1.7\,\! }[/math] and [math]\displaystyle{ \eta = 5000\,\! }[/math].

Therefore, the reliability of the system at some time, [math]\displaystyle{ t\,\! }[/math], is given by the following equation.

[math]\displaystyle{ \begin{align} R(t)= & {{R}_{1}}(t) +\underset{0}{\overset{t}{\mathop \int }}\,\{{{f}_{1}}(x)\cdot {{R}_{2;SB}}(x) \\ & \cdot \frac{{{R}_{2;A}}({{t}_{e}}+t-x)}{{{R}_{2;A}}({{t}_{e}})}\cdot {{R}_{SW;Q}}(x)\cdot {{R}_{SW;REQ}}(x)\}dx \end{align}\,\! }[/math]

where:

  • [math]\displaystyle{ {{R}_{1}}\,\! }[/math] is the reliability of the active component.
  • [math]\displaystyle{ {{f}_{1}}\,\! }[/math] is the pdf of the active component.
  • [math]\displaystyle{ {{R}_{2;SB}}\,\! }[/math] is the reliability of the standby component when in standby mode (quiescent reliability).
  • [math]\displaystyle{ {{R}_{2;A}}\,\! }[/math] is the reliability of the standby component when in active mode.
  • [math]\displaystyle{ {{R}_{SW;Q}}\,\! }[/math] is the quiescent reliability of the switch.
  • [math]\displaystyle{ {{R}_{SW;REQ}}\,\! }[/math] is the switch probability per request.
  • [math]\displaystyle{ {{t}_{e}}\,\! }[/math] is the equivalent operating time for the standby unit if it had been operating at an active mode.

This problem can be solved in BlockSim by including these probabilities in the container's properties, as shown in the figures below. In BlockSim, the standby container is acting as the switch.

Standby container (switch) failure probabilities while attempting to switch


Standby container (switch) failure distribution while waiting to switch

Note that there are additional properties that can be specified in BlockSim for a switch, such as Switch Restart Probability, No. of Restarts and Switch Delay Time. In many applications, the switch is re-tested (or re-cycled) if it fails to switch the first time. In these cases, it might be possible that it switches in the second or third, or [math]\displaystyle{ {{n}^{th}}\,\! }[/math] attempt.

The Switch Restart Probability specifies each additional attempt's probability of successfully switching and the Finite Restarts specifies the total number of attempts. Note that the Switch Restart Probability specifies the probability of success of each trial (or attempt). The probability of success of [math]\displaystyle{ n\,\! }[/math] consecutive trials is calculated by BlockSim using the binomial distribution and this probability is then incorporated into equation above. The Switch Delay Time property is related to repairable systems and is considered in BlockSim only when using simulation. When using the analytical solution (i.e., for a non-repairable system), this property is ignored.

Solving the analytical solution (as given by the above equation), the following results are obtained.

5-30.png

From the table above, it can be seen that the presence of a switching device has a significant effect on the reliability of a standby system. It is therefore important when modeling standby redundancy to incorporate the switching device reliability properties. It should be noted that this methodology is not the same as treating the switching device as another series component with the standby subsystem. This would be valid only if the failure of the switch resulted in the failure of system (e.g., switch failing open). In equation above, the Switch Probability per Request and quiescent probability are present only in the second term of the equation. Treating these two failure modes as a series configuration with the standby subsystem would imply that they are also present when the active component is functioning (i.e., first term of equation above). This is invalid and would result in the underestimation of the reliability of the system. In other words, these two failure modes become significant only when the active component fails.

As an example, and if we consider the warm standby case, the reliability of the system without the switch is 70.57% at 1000 hours. If the system was modeled so that the switching device was in series with the warm standby subsystem, the result would have been:

[math]\displaystyle{ \begin{align} {{R}_{S}}(1000)= & {{R}_{Standby}}(1000)\cdot {{R}_{sw,Q(1000)}}\cdot {{R}_{sw,req}} \\ = & 0.7057\cdot 0.9372\cdot 0.9 \\ = & 0.5952 \end{align}\,\! }[/math]

In the case where a switch failure mode causes the standby subsystem to fail, then this mode can be modeled as an individual block in series with the standby subsystem.

Example

Consider a car with four new tires and a full-size spare. Assume the following failure characteristics:

  • The tires follow a Weibull distribution with a [math]\displaystyle{ \beta=4\,\! }[/math] and an [math]\displaystyle{ \eta = 40,000\,\! }[/math] miles while on the car due to wear.
  • The tires also have a probability of failing due to puncture or other causes. For this, assume a constant rate for this occurrence with a probability of 1 every 50,000 miles.
  • When not on the car (i.e., is a spare), a tire's probability of failing also has a Weibull distribution with a [math]\displaystyle{ \beta = 2\,\! }[/math] and [math]\displaystyle{ \eta = 120,000\,\! }[/math] miles.

Assume a mission of 1,000 miles. If a tire fails during this trip, it will be replaced with the spare. However, the spare will not be repaired during the trip. In other words, the trip will continue with the spare on the car and, if the spare fails, the system will fail. Determine the probability of system failure.

Solution

The active failure distribution for the tires are:

  • Due to wearout, Weibull [math]\displaystyle{ \beta =4\,\! }[/math] and [math]\displaystyle{ \eta =40,000\,\! }[/math] miles.
  • Due to random puncture, exponential [math]\displaystyle{ \mu =50,000.\,\! }[/math]
  • The quiescent failure distribution is a Weibull distribution with [math]\displaystyle{ \beta =2\,\! }[/math] and [math]\displaystyle{ \eta =120,000\,\! }[/math] miles.

The block diagram for each tire has two blocks in series, one block representing the wearout mode and the other the random puncture mode, as shown next:


Small5.png


There are five tires, four active and one standby (represented in the diagram by a standby container with a 4-out-of-5 requirement), as shown next:


BStirecontainer.png


For the standby Wear block, set the active failure and the quiescent distributions, but for the Puncture block, set only the active puncture distribution (because the tire cannot fail due to puncture while stored). Using BlockSim, the probability of system failure is found to be 0.003 or 0.3%.

Examples heading.png

More examples on load sharing and standby configurations are available! See also:

Examples link.png Modeling Failure Modes

Examples link.png Load Sharing Configuration Example


Note Regarding Numerical Integration Solutions

Load sharing and standby solutions in BlockSim are performed using numerical integration routines. As with any numerical analysis routine, the solution error depends on the number of iterations performed, the step size chosen and related factors, plus the behavior of the underlying function. By default, BlockSim uses a certain set of preset factors. In general, these defaults are sufficient for most problems. If a higher precision or verification of the precision for a specific problem is required, BlockSim's preset options can be modified and/or the integration error can be assessed using the Integration Parameters option for each container. For more details, you can refer to the documentation on the Algorithm Setup window in the BlockSim help file.