Criticality Analysis: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
No edit summary
 
(30 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{Banner Xfmea Articles}}{{Navigation box}}
=== Introduction ===
Failure mode criticality, as used in FMECA, can be generically described as the joint probability of three events:  
Failure mode criticality, as used in FMECA, can be generically described as the joint probability of three events:  


#The probability of failure of the component;
#The probability of failure of the component.
#The probability that the mode under consideration was the culprit; and
#The probability that the mode under consideration was the culprit.
#The probability that this mode would result in a system failure (which for the case of non-redundancy or if in series =1).
#The probability that this mode would result in a system failure (which for the case of non-redundancy or if in series =1).


This number is then used to provide relative ranking for the different failure modes. This general approach is described in Kececioglu and expanded further here.  
This number is then used to provide relative rankings for the different failure modes. This general approach is described in Kececioglu and expanded further here.  


=== General Criticality Number Computation  ===
=== General Criticality Number Computation  ===
   
   
In the case that probability of occurrence of each mode is defined, and the modes are assumed to be in series, then a better way to obtain a criticality value  is:
::<math>{{C}_{{e}_{i}}}={{\beta }_{\bmod {{e}_{i}}}}\cdot \frac{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}}</math>


Mathematically, and for a specific failure mode,  
Mathematically, and for a specific failure mode,  


::<math>{{C}_{mode}}={{\beta }_{\bmod e}}\cdot {{\alpha }_{\bmod e}}\cdot Q{{(t)}_{component}} \,\!</math>
::<math>{{C}_{mode}}={{\beta }_{mode}}\cdot {{\alpha }_{mode}}\cdot Q{{(t)}_{component}} \,\!</math>
 
Where
* <math>{{C}_{component}}=</math>Criticality number for the failure mode.
* <math>{{\beta }_{\bmod e}}=</math>Conditional probability of mission loss.
* <math>{{\alpha }_{\bmod e}}=</math>Failure mode ratio.
* <math>Q{{(t)}_{\bmod e}}=</math>Unreliability or probability of failure due to that mode at time t.


In the case that the probability of occurrence of each mode is defined, and the modes are assumed to be in series, then a better way to obtain a criticality value is:
where:
* <math>{{C}_{mode}}= \,\!</math> Criticality number for the failure mode.
* <math>{{\beta }_{mode}}=  \,\!</math> Conditional probability of mission loss.
* <math>{{\alpha }_{mode}}= \,\!</math> Failure mode ratio.
* <math>Q{{(t)}_{mode}}= \,\!</math> Unreliability or probability of failure due to that mode at time t.
* <math>Q{{(t)}_{component}}= \,\!</math> Unreliability or probability of failure of the component at time t.


::<math>{{C}_{\bmod {{e}_{i}}}}={{\beta }_{\bmod {{e}_{i}}}}\cdot \frac{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}}</math>
In the case that the probability of occurrence of each mode is defined, and the modes are assumed to be in a series, then a better way to obtain a criticality value is:  


Note that this differs from the specialized definition given in Task 102 of MIL-STD-1629A, Section 3.2.1.6.
::<math>{{C}_{{{mode}_{i}}}}={{\beta }_{{{mode}_{i}}}}\cdot \frac{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}}\,\!</math>


Specifically a procedure is given for computing the “Failure mode criticality number” when an assumption of a constant failure mode (i.e., Exponential distribution) is used and for small values of lambda.


This is given by:
Note that this differs from the specialized definition given in Task 102 of MIL-STD-1629A, Section 3.2.1.6 (described next).


::<span class="texhtml" sab="496">''C''<sub sab="498">''m''</sub> = βαλ<sub sab="500">''P''</sub>''t''</span>
=== Task 102 MIL-STD-1629A ===


<br>where:  
Specifically, a procedure is given for computing the “Failure mode criticality number” when an assumption of a constant failure mode (i.e., an exponential distribution) is used and for small values of lambda:


*<span class="texhtml" sab="507">''C''<sub sab="509">''m''</sub> = </span>Criticality number for a failure mode.
This is given by:
*<span class="texhtml" sab="512">β = </span>Conditional probability of mission loss.
*<span class="texhtml" sab="514">α = </span>Failure mode ratio.
::<math>{{C}_{m}}=\beta \alpha {{\lambda }_{P}}t\,\!</math>
*<span class="texhtml" sab="516">λ<sub sab="517">''P''</sub> = </span>Part failure rate.
*<span class="texhtml" sab="520">''t'' = </span>Mission duration.


In addition to the constant failure rate assumption, it is important to note that this formulation assumes that the part failure rate values for a mode are relatively small and thus the following simplification is made:
where:
* <math>{{C}_{m}}=\,\!</math> Criticality number for a failure mode.
* <math>\beta =\,\!</math> Conditional probability of mission loss.
* <math>\alpha =\,\!</math> Failure mode ratio.
* <math>{{\lambda }_{P}}=\,\!</math> Part failure rate.
* <math>t=\,\!</math> Mission Duration.


<br>
In addition to the constant failure rate assumption, it is important to note that this formulation assumes that the part failure rate values for a mode are relatively small and thus the following simplification is made:


::<math>Q(t)=1-{{e}^{-{{\lambda }_{P}}\cdot t}}\approx {{\lambda }_{P}}\cdot t \,\!</math>
::<math>Q(t)=1-{{e}^{-{{\lambda }_{P}}\cdot t}}\approx {{\lambda }_{P}}\cdot t \,\!</math>


This approximation holds true when using an exponential distribution and small values of lambda (&lt;1E-4), but not for all values. ReliaSoft chose to use the general formulation instead of the specialized and limited 1629A formulation.
This approximation holds true when using an exponential distribution and small values of lambda (<1E-4), but not for all values.  
ReliaSoft chose to use the general formulation instead of the specialized and limited 1629A formulation.


=== Example  ===
===Example===
   
To illustrate this assume that a component has three failure modes (FM1, FM2 and FM3) where any mode can cause system failure (i.e., a series system).
The unreliability is defined by the failure model of each mode at t=1,000 hours. (In this case FM1 is Weibull
(<math>\beta=3\,\!</math>,<math>\eta=3000\,\!</math> hrs) and FM2 and 3 are exponential with MTTF=10,000 hours.) 


To illustrate this assume&nbsp;that a component has three failure modes, FM1, FM2 and FM3, and that any mode can cause system failure (i.e., a&nbsp;series system). The unreliability is defined by the failure model of each mode at t=1,000 hrs. (In this case FM1 is Weibull (<span class="texhtml" sab="451">β = 3</span>,<span class="texhtml" sab="452">η = 3000</span> hrs) and FM2 and 3 are Exponential with MTTF=10,000 hrs.)


<br>Then  
Then:


::<math>\begin{align}
::<math>\begin{align}
   & \sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}=\frac{0.036360+0.095163+0.095163}{0.211038} \\  
   & \sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}=\frac{0.036360+0.095163+0.095163}{0.211038} \\  
  & =\frac{0.226686}{0.211038} \\  
  & =\frac{0.226686}{0.211038} \\  
  & =1.074148   
  & =1.074148   
\end{align}</math>
\end{align}</math>


<br>
 


::<math>\begin{align}
::<math>\begin{align}
   & {{C}_{FM1}}={{\beta }_{FM1}}\cdot \frac{\frac{Q{{(t)}_{FM1}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}} \\  
   & {{C}_{FM1}}={{\beta }_{FM1}}\cdot \frac{\frac{Q{{(t)}_{FM1}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}} \\  
  & =1\cdot \frac{\frac{0.36360}{0.211038}}{1.074148} \\  
  & =1\cdot \frac{\frac{0.36360}{0.211038}}{1.074148} \\  
  & =0.160398   
  & =0.160398   
\end{align}</math>
\end{align}</math>


<br>


::<math>\begin{align}
::<math>\begin{align}
   & {{C}_{FM2}}={{\beta }_{FM2}}\cdot \frac{\frac{Q{{(t)}_{FM2}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}} \\  
   & {{C}_{FM2}}={{\beta }_{FM2}}\cdot \frac{\frac{Q{{(t)}_{FM2}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}} \\  
  & =1\cdot \frac{\frac{0.095163}{0.211038}}{1.074148} \\  
  & =1\cdot \frac{\frac{0.095163}{0.211038}}{1.074148} \\  
  & =0.419801   
  & =0.419801   
\end{align}</math>
\end{align}</math>


and  
and:


::<math>\begin{align}
::<math>\begin{align}
   & {{C}_{FM3}}={{\beta }_{FM3}}\cdot \frac{\frac{Q{{(t)}_{FM3}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{\bmod {{e}_{i}}}}}{Q{{(t)}_{Component}}}}} \\  
   & {{C}_{FM3}}={{\beta }_{FM3}}\cdot \frac{\frac{Q{{(t)}_{FM3}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}} \\  
  & =1\cdot \frac{\frac{0.095163}{0.211038}}{1.074148} \\  
  & =1\cdot \frac{\frac{0.095163}{0.211038}}{1.074148} \\  
  & =0.419801   
  & =0.419801   
\end{align}</math>
\end{align}</math>


=== Derivation ===
===Derivation===
 
Given that&nbsp;the component failed before time T, the probability that it is caused by mode B is:<br>


<br>
Given that the component failed before time T, the probability that it is caused by mode B is:


::<math>\begin{align}
::<math>\begin{align}
Line 109: Line 101:
\end{align}</math>
\end{align}</math>


<br>
The upper bound of P1 is:
 
<br>The upper bound on P<sub>1</sub> is:  


::<math>{{P}_{1}}=\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{{{Q}_{component}}}<\frac{\int_{0}^{T}{{{f}_{B}}(t)dt}}{{{Q}_{component}}}=\frac{{{Q}_{B}}}{{{Q}_{component}}}=\frac{{{Q}_{B}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{AB}}}</math>.
::<math>{{P}_{1}}=\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{{{Q}_{component}}}<\frac{\int_{0}^{T}{{{f}_{B}}(t)dt}}{{{Q}_{component}}}=\frac{{{Q}_{B}}}{{{Q}_{component}}}=\frac{{{Q}_{B}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{AB}}}</math>.


<br>The lower bound of P<sub>1</sub> is:  
The lower bound of P1 is:


::<math>\begin{align}
::<math>\begin{align}
Line 122: Line 112:
\end{align}</math>
\end{align}</math>


Thus, the probability that the failure is caused by mode B is  
Thus, the probability that the failure is caused by mode B is:
 
::<math>\frac{{{Q}_{B}}-{{Q}_{B}}{{Q}_{A}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{A}}{{Q}_{B}}}<{{P}_{1}}<\frac{{{Q}_{B}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{AB}}}</math>
::<math>\frac{{{Q}_{B}}-{{Q}_{B}}{{Q}_{A}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{A}}{{Q}_{B}}}<{{P}_{1}}<\frac{{{Q}_{B}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{AB}}}</math>


<br>Thus the approximation  
Thus the approximation:
 
::<math>{{\tilde{P}}_{1}}=\frac{{{Q}_{B}}}{{{Q}_{B}}+{{Q}_{A}}}</math>
::<math>{{\tilde{P}}_{1}}=\frac{{{Q}_{B}}}{{{Q}_{B}}+{{Q}_{A}}}</math>
is within the bounds.  
is within the bounds.  


<br>Similarly, given that the component failed before time T, the probability that it is caused by mode A is:  
Similarly, given that the component failed before time T, the probability that it is caused by mode A is:
 
<br>


::<math>\begin{align}
::<math>\begin{align}
Line 141: Line 126:
  & =\frac{\int_{0}^{T}{{{f}_{A}}(t){{R}_{B}}(t)dt}}{1-{{R}_{A}}(T){{R}_{B}}(T)}   
  & =\frac{\int_{0}^{T}{{{f}_{A}}(t){{R}_{B}}(t)dt}}{1-{{R}_{A}}(T){{R}_{B}}(T)}   
\end{align}</math>
\end{align}</math>
 
<br>Then the sum of P<sub>1</sub> and P<sub>2</sub> is:  
Then the sum of P1 and P2 is:


::<math>\begin{align}
::<math>\begin{align}
Line 152: Line 137:
  & =\frac{{{F}_{A}}(T)+{{F}_{B}}(T)-{{F}_{A}}(T){{F}_{B}}(T)}{1-{{R}_{A}}(T){{R}_{B}}(T)} \\  
  & =\frac{{{F}_{A}}(T)+{{F}_{B}}(T)-{{F}_{A}}(T){{F}_{B}}(T)}{1-{{R}_{A}}(T){{R}_{B}}(T)} \\  
  & =1 \\  
  & =1 \\  
\end{align}</math>
\end{align}</math>  
 
Thus the approximation


Thus the approximation:
::<math>{{\tilde{P}}_{2}}=\frac{{{Q}_{A}}}{{{Q}_{B}}+{{Q}_{A}}}</math>
::<math>{{\tilde{P}}_{2}}=\frac{{{Q}_{A}}}{{{Q}_{B}}+{{Q}_{A}}}</math>
is within these bounds.


is within these bounds.
Also note:
 
::<math>{{\tilde{P}}_{1}}+{{\tilde{P}}_{2}}=1</math>
Also note that
 
::<math>{{\tilde{P}}_{1}}+{{\tilde{P}}_{2}}=1</math>.

Latest revision as of 05:19, 23 August 2012

Xfmea Articles Banner.png


Introduction

Failure mode criticality, as used in FMECA, can be generically described as the joint probability of three events:

  1. The probability of failure of the component.
  2. The probability that the mode under consideration was the culprit.
  3. The probability that this mode would result in a system failure (which for the case of non-redundancy or if in series =1).

This number is then used to provide relative rankings for the different failure modes. This general approach is described in Kececioglu and expanded further here.

General Criticality Number Computation

Mathematically, and for a specific failure mode,

[math]\displaystyle{ {{C}_{mode}}={{\beta }_{mode}}\cdot {{\alpha }_{mode}}\cdot Q{{(t)}_{component}} \,\! }[/math]

where:

  • [math]\displaystyle{ {{C}_{mode}}= \,\! }[/math] Criticality number for the failure mode.
  • [math]\displaystyle{ {{\beta }_{mode}}= \,\! }[/math] Conditional probability of mission loss.
  • [math]\displaystyle{ {{\alpha }_{mode}}= \,\! }[/math] Failure mode ratio.
  • [math]\displaystyle{ Q{{(t)}_{mode}}= \,\! }[/math] Unreliability or probability of failure due to that mode at time t.
  • [math]\displaystyle{ Q{{(t)}_{component}}= \,\! }[/math] Unreliability or probability of failure of the component at time t.

In the case that the probability of occurrence of each mode is defined, and the modes are assumed to be in a series, then a better way to obtain a criticality value is:

[math]\displaystyle{ {{C}_{{{mode}_{i}}}}={{\beta }_{{{mode}_{i}}}}\cdot \frac{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}}\,\! }[/math]


Note that this differs from the specialized definition given in Task 102 of MIL-STD-1629A, Section 3.2.1.6 (described next).

Task 102 MIL-STD-1629A

Specifically, a procedure is given for computing the “Failure mode criticality number” when an assumption of a constant failure mode (i.e., an exponential distribution) is used and for small values of lambda:

This is given by:

[math]\displaystyle{ {{C}_{m}}=\beta \alpha {{\lambda }_{P}}t\,\! }[/math]

where:

  • [math]\displaystyle{ {{C}_{m}}=\,\! }[/math] Criticality number for a failure mode.
  • [math]\displaystyle{ \beta =\,\! }[/math] Conditional probability of mission loss.
  • [math]\displaystyle{ \alpha =\,\! }[/math] Failure mode ratio.
  • [math]\displaystyle{ {{\lambda }_{P}}=\,\! }[/math] Part failure rate.
  • [math]\displaystyle{ t=\,\! }[/math] Mission Duration.

In addition to the constant failure rate assumption, it is important to note that this formulation assumes that the part failure rate values for a mode are relatively small and thus the following simplification is made:

[math]\displaystyle{ Q(t)=1-{{e}^{-{{\lambda }_{P}}\cdot t}}\approx {{\lambda }_{P}}\cdot t \,\! }[/math]

This approximation holds true when using an exponential distribution and small values of lambda (<1E-4), but not for all values. ReliaSoft chose to use the general formulation instead of the specialized and limited 1629A formulation.

Example

To illustrate this assume that a component has three failure modes (FM1, FM2 and FM3) where any mode can cause system failure (i.e., a series system). The unreliability is defined by the failure model of each mode at t=1,000 hours. (In this case FM1 is Weibull ([math]\displaystyle{ \beta=3\,\! }[/math],[math]\displaystyle{ \eta=3000\,\! }[/math] hrs) and FM2 and 3 are exponential with MTTF=10,000 hours.)


Then:

[math]\displaystyle{ \begin{align} & \sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}=\frac{0.036360+0.095163+0.095163}{0.211038} \\ & =\frac{0.226686}{0.211038} \\ & =1.074148 \end{align} }[/math]


[math]\displaystyle{ \begin{align} & {{C}_{FM1}}={{\beta }_{FM1}}\cdot \frac{\frac{Q{{(t)}_{FM1}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}} \\ & =1\cdot \frac{\frac{0.36360}{0.211038}}{1.074148} \\ & =0.160398 \end{align} }[/math]


[math]\displaystyle{ \begin{align} & {{C}_{FM2}}={{\beta }_{FM2}}\cdot \frac{\frac{Q{{(t)}_{FM2}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}} \\ & =1\cdot \frac{\frac{0.095163}{0.211038}}{1.074148} \\ & =0.419801 \end{align} }[/math]

and:

[math]\displaystyle{ \begin{align} & {{C}_{FM3}}={{\beta }_{FM3}}\cdot \frac{\frac{Q{{(t)}_{FM3}}}{Q{{(t)}_{Component}}}}{\sum\limits_{i=1}^{M}{\frac{Q{{(t)}_{{{mode}_{i}}}}}{Q{{(t)}_{Component}}}}} \\ & =1\cdot \frac{\frac{0.095163}{0.211038}}{1.074148} \\ & =0.419801 \end{align} }[/math]

Derivation

Given that the component failed before time T, the probability that it is caused by mode B is:

[math]\displaystyle{ \begin{align} & {{P}_{1}}=\operatorname{P}(A\gt B|\text{At least one of either A or B is less than T)} \\ & \text{=}\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{P(\text{At least one of A and B is less than T})} \\ & =\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{1-{{R}_{A}}(T){{R}_{B}}(T)}=\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{{{Q}_{component}}} \end{align} }[/math]

The upper bound of P1 is:

[math]\displaystyle{ {{P}_{1}}=\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{{{Q}_{component}}}\lt \frac{\int_{0}^{T}{{{f}_{B}}(t)dt}}{{{Q}_{component}}}=\frac{{{Q}_{B}}}{{{Q}_{component}}}=\frac{{{Q}_{B}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{AB}}} }[/math].

The lower bound of P1 is:

[math]\displaystyle{ \begin{align} & {{P}_{1}}=\frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{{{Q}_{component}}}\gt \frac{\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(T)dt}}{{{Q}_{component}}}=\frac{{{Q}_{B}}{{R}_{A}}(T)}{{{Q}_{component}}}=\frac{{{Q}_{B}}(1-{{Q}_{A}})}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{A}}{{Q}_{B}}} \\ & =\frac{{{Q}_{B}}-{{Q}_{B}}{{Q}_{A}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{A}}{{Q}_{B}}} \\ \end{align} }[/math]

Thus, the probability that the failure is caused by mode B is:

[math]\displaystyle{ \frac{{{Q}_{B}}-{{Q}_{B}}{{Q}_{A}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{A}}{{Q}_{B}}}\lt {{P}_{1}}\lt \frac{{{Q}_{B}}}{{{Q}_{A}}+{{Q}_{B}}-{{Q}_{AB}}} }[/math]

Thus the approximation:

[math]\displaystyle{ {{\tilde{P}}_{1}}=\frac{{{Q}_{B}}}{{{Q}_{B}}+{{Q}_{A}}} }[/math]

is within the bounds.

Similarly, given that the component failed before time T, the probability that it is caused by mode A is:

[math]\displaystyle{ \begin{align} & {{P}_{2}}=\operatorname{P}(B\gt A|\text{At least one of A and B is less than T)} \\ & \text{=}\frac{\int_{0}^{T}{{{f}_{A}}(t){{R}_{B}}(t)dt}}{P(\text{At least one of A and B is less than T})} \\ & =\frac{\int_{0}^{T}{{{f}_{A}}(t){{R}_{B}}(t)dt}}{1-{{R}_{A}}(T){{R}_{B}}(T)} \end{align} }[/math]

Then the sum of P1 and P2 is:

[math]\displaystyle{ \begin{align} & {{P}_{1}}+{{P}_{2}}=\frac{\int_{0}^{T}{{{f}_{A}}(t){{R}_{B}}(t)dt}+\int_{0}^{T}{{{f}_{B}}(t){{R}_{A}}(t)dt}}{1-{{R}_{A}}(T){{R}_{B}}(T)} \\ & =\frac{{{F}_{A}}(T)-\int_{0}^{T}{{{f}_{A}}(t){{F}_{B}}(t)dt}+{{F}_{B}}(T)-\left[ {{F}_{A}}(T){{F}_{B}}(T)\left| \begin{align} & T \\ & 0 \\ \end{align} \right.-\int_{0}^{T}{{{F}_{B}}(T){{f}_{A}}(T)dt} \right]}{1-{{R}_{A}}(T){{R}_{B}}(T)} \\ & =\frac{{{F}_{A}}(T)+{{F}_{B}}(T)-{{F}_{A}}(T){{F}_{B}}(T)}{1-{{R}_{A}}(T){{R}_{B}}(T)} \\ & =1 \\ \end{align} }[/math]

Thus the approximation:

[math]\displaystyle{ {{\tilde{P}}_{2}}=\frac{{{Q}_{A}}}{{{Q}_{B}}+{{Q}_{A}}} }[/math]

is within these bounds.

Also note:

[math]\displaystyle{ {{\tilde{P}}_{1}}+{{\tilde{P}}_{2}}=1 }[/math]