Template:Crow extended model rga: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
 
(7 intermediate revisions by 2 users not shown)
Line 1: Line 1:
==Crow Extended Model==
#REDIRECT [[Reliability_Growth_Planning#Crow_Extended_Model]]
===Introduction===
The Crow Extended model for reliability growth planning is a revised and improved version of the MIL-HDBK-189 growth curve [13]. MIL-HDBK-189 can be considered as the growth curve based on the Crow-AMSAA (NHPP) model. Using MIL-HDBK-189 for reliability growth planning assumes that the corrective actions for the observed failure modes are incorporated during the test and at the specific time of failure. However, in actual practice, fixes may be delayed until after the completion of the test or some fixes may be implemented during the test while others are delayed and some are not fixed at all. The Crow Extended model for reliability growth planning provides additional input to be able to account for a specific management strategy and delayed fixes with specified effectiveness factors.
<br>
<br>
===Definitions===
<br>
{{management strategy ratio crow extended}}
 
{{effectiveness factor crow extended}}
 
====MTBF Goal====
<br>
When putting together a reliability growth plan, a goal MTBF  <math>{{M}_{G}}</math>  (or goal failure intensity  <math>{{\lambda }_{G}}</math> ) is defined as the requirement or target for the product at the end of the growth program.
<br>
====Growth Potential====
The failure intensity remaining in the system at the end of the test will depend on the management strategy given by the classification of the Type A and Type B failure modes. The engineering effort applied to the corrective actions determines the effectiveness factors. In addition, the failure intensity depends on  <math>h(t)</math> , which is the rate at which problem failure modes are being discovered during testing. The rate of discovery drives the opportunity to take corrective actions based on the seen failure modes and it is an important factor in the overall reliability growth rate. The reliability growth potential is the limiting value of the failure intensity as time  <math>T</math>  increases. This limit is the maximum MTBF that can be attained with the current management strategy. The maximum MTBF will be attained when all Type B modes have been observed and fixed.
<br>
If all seen Type B modes are corrected by time  <math>T</math> , that is, no deferred corrective actions at time  <math>T</math> , then the growth potential is the maximum attainable with the Type B designation of the failure modes and the corresponding assigned effectiveness factors. This is called the nominal growth potential. In other words, the nominal growth potential is the maximum attainable growth potential assuming corrective actions are implemented for every mode that is planned to be fixed. In reality, some fixes to modes might be implemented at a later time due to schedule, budget, engineering, etc.
<br>
If some seen Type B modes are not corrected at the end of the current test phase then the prevailing growth potential is below the maximum attainable with the Type B designation of the failure modes and the corresponding assigned effectiveness factors.
If all Type B failure modes are seen and corrected with an average effectiveness factor,  <math>d</math> , then the maximum reduction in the initial system failure intensity is the growth potential failure intensity:
 
 
::<math>{{\lambda }_{GP}}={{\lambda }_{A}}+\left( 1-d \right){{\lambda }_{B}}</math>
 
 
The growth potential MTBF is:
 
 
::<math>{{M}_{GP}}=\frac{1}{{{\lambda }_{GP}}}</math>
 
 
Note that based Eqns. (lambda GP), (lambda initial) and (msr), the initial failure intensity is equal to:
 
 
::<math>{{\lambda }_{I}}=\frac{{{\lambda }_{GP}}}{1-d\cdot msr}</math>
 
 
====Growth Potential Design Margin====
<br>
The Growth Potential Design Margin ( <math>GPDM</math> ) can be considered as a safety margin when setting target MTBF values for the reliability growth plan. It is common for systems to degrade in terms of reliability when a prototype product is going into full manufacturing due to variation in material, processes, etc. Furthermore, the in-house reliability growth testing usually overestimates the actual product reliability, since the field usage conditions may not be perfectly simulated during growth testing. Typical values for the  <math>GPDM</math>  are around 1.2. Higher values yield less risk for the program, but require a more rigorous reliability growth test plan. Lower values imply higher program risk, with less safety margin.
 
During the planning stage, the growth potential MTBF,  <math>{{M}_{GP}},</math>  can be calculated based on the goal MTBF,  <math>{{M}_{G}},</math>  and the growth potential design margin,  <math>GPDM</math> :
 
 
::<math>{{M}_{GP}}=GPDM\cdot {{M}_{G}}</math>
 
 
or in terms of failure intensity:
 
 
::<math>{{\lambda }_{GP}}=\frac{{{\lambda }_{G}}}{GPDM}</math>
 
 
===Nominal Idealized Growth Curve===
<br>
During developmental testing, management should expect that certain levels of reliability will be attained at various points in the program in order to have assurance that reliability growth is progressing at a sufficient rate to meet the product reliability requirement. The idealized curve portrays an overall characteristic pattern, which is used to determine and evaluate intermediate levels of reliability and construct the program planned growth curve. Note that growth profiles on previously developed, similar systems provide significant insight into the reliability growth process and are valuable in the construction of idealized growth curves.
The nominal idealized growth curve portrays a general profile for reliability growth throughout system testing. The idealized curve has the baseline value  <math>{{\lambda }_{I}}</math>  until an initialization time,  <math>{{t}_{0}},</math>  when reliability growth occurs. From that time and until the end of testing, which can be one or, most commonly, multiple test phases, the idealized curve increases steadily according to a learning curve pattern until it reaches the final reliability requirement,  <math>{{M}_{F}}</math> . The slope of this curve on a log-log plot is the growth rate of the Crow Extended model [13].
 
====Nominal Failure Intensity Function====
The nominal idealized growth curve failure intensity as a function of test time  <math>t</math>  is:
 
::<math>{{r}_{NI}}(t)={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta {{t}^{\left( \beta -1 \right)}}\text{ for }t\ge {{t}_{0}}</math>
 
:and:
 
::<math>{{r}_{NI}}(t)={{\lambda }_{I}}\text{ for }t\le {{t}_{0}}</math>
 
where  <math>{{\lambda }_{I}}</math>  is the initial system failure intensity,  <math>t</math>  is test time and  <math>{{t}_{0}}</math>  is the initialization time, which is discussed in the next section.
<br>
It can be seen that Eqn. (Nominal FI) is the failure intensity equation of the Crow Extended model.
<br>
<br>
 
====Initialization Time====
<br>
Reliability growth can only begin after a Type B failure mode occurs, which cannot be at a time equal to zero. Therefore, there is a need for an initialization time, different than zero, to be defined. The nominal idealized growth curve failure intensity is initially set equal to the initial failure intensity,  <math>{{\lambda }_{I}},</math>  until the initialization time,  <math>{{t}_{0}}</math> :
 
 
::<math>{{r}_{NI}}({{t}_{0}})={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta t_{0}^{(\beta -1)}</math>
 
:Therefore:
 
::<math>{{\lambda }_{I}}={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta t_{0}^{(\beta -1)}</math>
 
:Then:
 
::<math>{{t}_{0}}={{\left[ \frac{{{\lambda }_{I}}-{{\lambda }_{A}}-(1-d){{\lambda }_{B}}}{d\lambda \beta } \right]}^{\tfrac{1}{\beta -1}}}</math>
 
Using Eqn. (lambda initial) to substitute  <math>{{\lambda }_{I}}</math>  we have:
 
::<math>{{t}_{0}}={{\left[ \frac{{{\lambda }_{A}}+{{\lambda }_{B}}-{{\lambda }_{A}}-(1-d){{\lambda }_{B}}}{d\cdot \lambda \cdot \beta } \right]}^{\tfrac{1}{\beta -1}}}</math>
 
:Then:
 
::<math>{{t}_{0}}={{\left( \frac{{{\lambda }_{B}}}{\lambda \cdot \beta } \right)}^{\tfrac{1}{\beta -1}}}</math>
 
The initialization time,  <math>{{t}_{0}},</math>  allows for growth to start after a Type B failure mode has occurred.
<br>
<br>
 
====Nominal Time to Reach Goal====
<br>
Assuming that we have a target MTBF or failure intensity goal, we can solve Eqn. (Nominal FI) to find out how much test time,  , is required, (based on the Crow Extended model and the nominal idealized growth curve) to reach that goal:
 
 
::<math>{{t}_{N,G}}={{\left[ \frac{{{r}_{G}}-{{\lambda }_{A}}-(1-d){{\lambda }_{B}}}{d\cdot \lambda \cdot \beta } \right]}^{\tfrac{1}{\beta -1}}}</math>
 
 
Note that when  <math>{{\lambda }_{I}}<{{r}_{G}}</math>  or, in other words, the initial failure intensity is lower than the goal failure intensity, there is no need to solve for the nominal time to reach the goal, because the goal is already met. In this case, no further reliability growth testing is needed.
====Growth Rate for Nominal Idealized Curve====
<br>
The growth rate for the nominal idealized curve is defined in the same context as the growth rate for the Duane Postulate [8]. The nominal idealized curve has the same functional form for the growth rate as the Duane Postulate and the Crow-AMSAA (NHPP) model.
For both the Duane Postulate and the Crow-AMSAA (NHPP) model, the average failure intensity is given by:
 
 
::<math>C(t)=\frac{\lambda {{t}^{\beta }}}{t}=\lambda {{t}^{(\beta -1)}}</math>
 
 
Also, for both the Duane Postulate and the Crow-AMSAA (NHPP) model, the instantaneous failure intensity is given by:
 
 
::<math>r(t)=\lambda \beta {{t}^{(\beta -1)}}</math>
 
 
Taking the difference,  <math>D(t),</math>  between the average failure intensity,  <math>C(t)</math> , and the instantaneous failure intensity,  <math>r(t)</math> , yields:
 
 
::<math>D(t)=\lambda {{t}^{(\beta -1)}}-\lambda \beta {{t}^{(\beta -1)}}</math>
 
 
:Then:
 
 
::<math>D(t)=\lambda {{t}^{(\beta -1)}}[1-\beta ]</math>
 
 
For reliability growth to occur,  <math>D(t)</math>  must be decreasing.
The growth rate for both the Duane Postulate and the Crow-AMSAA (NHPP) model is the negative of the slope of  <math>\log (D(t))</math>  as a function of  <math>\log (t)</math> :
 
 
::<math>{{\log }_{e}}(D(t))=\text{constant}-(1-\beta ){{\log }_{e}}(t)</math>
 
 
 
The slope is negative under reliability growth and equals:
 
 
::<math>\text{slope}=-(1-\beta )</math>
 
 
The growth rate for both the Duane Postulate and the Crow-AMSAA (NHPP) model is equal to the negative of this slope:
 
 
::<math>\text{Growth Rate}=(1-\beta )</math>
 
 
The instantaneous failure intensity for the nominal idealized curve is:
 
 
::<math>{{r}_{NI}}(t)={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta {{(t)}^{(\beta -1)}}</math>
 
 
The cumulative failure intensity for the nominal idealized curve is:
 
 
::<math>{{C}_{NI}}(t)={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda {{(t)}^{(\beta -1)}}</math>
 
 
:Therefore:
 
 
::<math>{{D}_{NI}}(t)=[{{C}_{NI}}(t)-{{r}_{NI}}(t)]=\lambda {{t}^{(\beta -1)}}[1-\beta ]</math>
 
 
:and:
 
 
::<math>{{\log }_{e}}({{D}_{NI}}(t))=\text{constant}-(1-\beta ){{\log }_{e}}(t)</math>
 
 
Therefore, in accordance with the Duane Postulate and the Crow-AMSAA (NHPP) model,  <math>a=1-\beta </math>  is the growth rate for the reliability growth plan.
<br>
 
====Lambda - Beta Parameter Relationship====
<br>
Under the Crow-AMSAA (NHPP) model, the time to first failure is a Weibull random variable. The MTTF of a Weibull distributed random variable with parameters  <math>\beta </math>  and  <math>\eta </math>  is:
 
 
::<math>MTTF=\eta \cdot \Gamma \left( 1+\frac{1}{\beta } \right)</math>
 
 
The parameter lambda is defined as:
 
 
::<math>\lambda =\frac{1}{{{\eta }^{\beta }}}</math>
 
 
Using Eqn. (lambda eta relationship), the MTTF relationship shown in Eqn. (Weibull MTTF) becomes:
 
 
::<math>MTB{{F}_{B}}=\frac{\Gamma \left( 1+\tfrac{1}{\beta } \right)}{{{\lambda }^{\left( \tfrac{1}{\beta } \right)}}}</math>
 
 
Or, in terms of failure intensity:
 
 
::<math>{{\lambda }_{B}}=\frac{{{\lambda }^{\left( \tfrac{1}{\beta } \right)}}}{\Gamma \left( 1+\tfrac{1}{\beta } \right)}</math>
 
 
===Actual Idealized Growth Curve===
<br>
The actual idealized growth curve differs from the nominal idealized curve in that it takes into account the average fix delay that might occur in each test phase. The actual idealized growth curve is continuous and goes through each of the test phase target MTBFs.
 
====Fix Delays and Test Phase Target MTBF====
<br>
Fix delays reflect how long it takes from the time a problem failure mode is discovered in testing, to the time the corrective action is incorporated into the system and reliability growth is realized. The consideration of the fix delay is often in terms of how much calendar time it takes to incorporate a corrective action fix after the problem is first seem. However, the impact of the delay on reliability growth is reflected in the average test time it takes between finding a problem failure mode and incorporating a corrective action. The fix delay is reflected in the actual idealized growth curve in terms of test time.
In other words, the average fix delay is calendar time converted to test hours. For example, say that we expect an average fix delay of two weeks: if in two weeks the total test time is 1000 hours, the average fix delay is 1000 hours. If in the same two weeks the total test time is 2000 hours (maybe there are more units available or more shifts) then the average fix delay is 2000 hours.
<br>
There can be a constant fix delay across all test phases or, as a practical matter, each test phase can have a different fix delay time. In practice, the fix delay will generally be constant over a particular test phase.  <math>{{L}_{i}}</math>  denotes the fix delay for phase  <math>i=1,...,P,</math>  where  <math>P</math>  is the total number of phases in the test. RGA 7 allows for a maximum of seven test phases.
<br>
<br>
====Actual Failure Intensity Function====
<br>
Consider a test plan consisting of  <math>i</math>  phases. Taking into account the fix delay within each phase, we expect the actual failure intensity to be different (i.e. shifted) from the nominal failure intensity. This is because fixes are not incorporated instantaneously, thus growth is realized at a later time compared to the nominal case.
Specifically, the actual failure intensity will be estimated as follows:
 
'''Test Phase 1'''
<br>
For the first phase of a test plan, the actual idealized curve failure intensity,  <math>{{r}_{AI}}(t)</math> , is <math>:</math>
 
 
::<math>{{r}_{AI}}(t)={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta {{\left[ \left( \frac{{{T}_{1}}-{{L}_{1}}}{{{T}_{1}}} \right)t \right]}^{(\beta -1)}}\text{ for }0<t\le {{T}_{1}}</math>
 
 
Note that the end time of Phase 1,  <math>{{T}_{1}},</math>  must be greater than  <math>{{L}_{1}}+{{t}_{0}}</math> . That is,  <math>{{T}_{1}}>{{L}_{1}}+{{t}_{0}}</math> .
The actual idealized curve initialization time for Phase 1,  <math>T_{0}^{AIC},</math>  is calculated from:
 
 
::<math>{{r}_{AI}}(T_{0}^{AIC})={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta {{\left[ \left( \frac{{{T}_{1}}-{{L}_{1}}}{{{T}_{1}}} \right)T_{0}^{AIC} \right]}^{(\beta -1)}}</math>
 
 
Where  <math>{{r}_{AI}}(T_{0}^{AIC})={{r}_{I}}.</math>
Therefore, using Eqn. (lambdaqualsominal):
 
 
::<math>{{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta {{\left[ \left( \frac{{{T}_{1}}-{{L}_{1}}}{{{T}_{1}}} \right)T_{0}^{AIC} \right]}^{(\beta -1)}}={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta t_{0}^{(\beta -1)}</math>
 
 
Solving Eqn. (ActIhase1) for  <math>T_{0}^{AIC}</math>  we get:
 
 
::<math>T_{0}^{AIC}=\frac{{{t}_{0}}}{\left( \tfrac{{{T}_{1}}-{{L}_{1}}}{{{T}_{1}}} \right)}</math>
 
 
 
'''Test Phase  <math>i</math>'''
<br>
For any test phase  <math>i</math> , the actual idealized curve failure intensity is given by:
 
 
::<math>{{r}_{AI}}(t)={{\lambda }_{A}}+(1-d){{\lambda }_{B}}+d\lambda \beta {{\left[ {{T}_{i-1}}-{{L}_{i-1}}+\left( \frac{{{T}_{i}}-{{L}_{i}}-{{T}_{i-1}}+{{L}_{i-1}}}{{{T}_{i}}-{{T}_{i-1}}} \right)(t-{{T}_{i-1}}) \right]}^{(\beta -1)}}</math>
 
 
where  <math>{{T}_{i-1}}\le t\le {{T}_{i}}</math>  and  <math>{{T}_{i}}</math>  is the test time of each corresponding test phase.
<br>
The actual idealized curve MTBF is:
 
 
::<math>{{M}_{AI}}=\frac{1}{{{r}_{AI}}(t)}</math>
 
 
====Actual Time to Reach Goal====
<br>
The actual time to reach the target MTBF or failure intensity goal,  <math>{{t}_{AC,G}},</math>  can be found by solving Eqn. (ActIhase):
 
 
::<math>\begin{align}
  & {{r}_{AI}}({{t}_{AC,G}})= & {{\lambda }_{A}}+(1-d){{\lambda }_{B}} \\
&  & +d\lambda \beta {{\left[ {{T}_{i-1}}-{{L}_{i-1}}+\left( \frac{{{T}_{i}}-{{L}_{i}}-{{T}_{i-1}}+{{L}_{i-1}}}{{{T}_{i}}-{{T}_{i-1}}} \right)({{t}_{AC,G}}-{{T}_{i-1}}) \right]}^{(\beta -1)}} 
\end{align}</math>
 
 
Since the actual idealized growth curve depends on the phase durations and average fix delays, there are three different cases that need to be treated differently in order to determine the actual time to reach the MTBF goal. The cases depend on when the actual MTBF that can be reached within the specific phase durations and fix delays becomes equal to the MTBF goal. This can be determined by solving Eqn. (ActIhase) for phases  <math>1</math>  through  <math>i</math> , then solving in terms of actual MTBF using Eqn. (MTBFctualq) for each phase and finding the phase during which the actual MTBF becomes equal to the goal MTBF. The three cases are presented next.
 
''Case 1: MTBF goal is met during the last phase''
<br>
If  <math>{{T}_{F}}</math>  indicates the cumulative end phase time for the last phase and  <math>{{L}_{F}}</math>  indicates the fix delay for the last phase, then we have:
 
 
::<math>\begin{align}
  & {{r}_{G}}= & {{\lambda }_{A}}+(1-d){{\lambda }_{B}} \\
&  & +d\lambda \beta {{\left[ {{T}_{F-1}}-{{L}_{F-1}}+\left( \frac{{{T}_{F}}-{{L}_{F}}-{{T}_{F-1}}+{{L}_{F-1}}}{{{T}_{F}}-{{T}_{F-1}}} \right)({{t}_{AC,G}}-{{T}_{F-1}}) \right]}^{(\beta -1)}} 
\end{align}</math>
 
 
Starting to solve for  <math>{{t}_{AC,G}}</math>  yields:
 
 
::<math>{{\left[ \frac{{{r}_{G}}-{{\lambda }_{A}}-(1-d){{\lambda }_{B}}}{d\lambda \beta } \right]}^{\tfrac{1}{\beta -1}}}={{T}_{F-1}}-{{L}_{F-1}}+\left( \frac{{{T}_{F}}-{{L}_{i}}-{{T}_{F-1}}+{{L}_{F-1}}}{{{T}_{F}}-{{T}_{F-1}}} \right)({{t}_{AC,G}}-{{T}_{F-1}})</math>
 
 
We can substitute the left term by using Eqn. (tgoalominal), thus we have:
 
 
::<math>{{t}_{N,G}}={{T}_{F-1}}-{{L}_{F-1}}+\left( \frac{{{T}_{F}}-{{L}_{F}}-{{T}_{F-1}}+{{L}_{F-1}}}{{{T}_{F}}-{{T}_{F-1}}} \right)({{t}_{AC,G}}-{{T}_{i-1}})</math>
 
 
:Therefore:
 
 
::<math>{{t}_{AC,G}}=\frac{{{t}_{N,G}}-{{T}_{F-1}}+{{L}_{F-1}}}{\left( \tfrac{{{T}_{F}}-{{L}_{F}}-{{T}_{F-1}}+{{L}_{F-1}}}{{{T}_{F}}-{{T}_{F-1}}} \right)}+{{T}_{F-1}}</math>
 
 
''Case 2: MTBF goal is met before the last phase''
<br>
Eqn. (toalctual) still applies, but in this case  <math>{{T}_{F}}</math>  and  <math>{{L}_{F}}</math>  are the time and fix delay of the phase during which the goal is met.
<br>
<br>
''Case 3: MTBF goal is met after the final phase''
<br>
If the goal MTBF,  <math>{{M}_{G}},</math>  is met after the final test phase, then the actual time to reach the goal is not calculated, since additional phases have to be added with specific duration and fix delays. The reliability growth program needs to be re-evaluated with the following options:
<br>
<br>
:• Add more phase(s) to the program.
:• Re-examine the phase duration of the existing phases.
:• Investigate whether there are potential process improvements in the program that can <br>
reduce the average fix delay for the phases.
<br>
Other alternative routes for consideration would be to investigate the rest of the inputs in the model:
<br>
<br>
:• Change the management strategy.
:• Consider if further program risk can be acceptable, and if so, reduce the growth potential design margin.
:• Consider if it is feasible to increase the effectiveness factors of the delayed fixes by using more robust engineering redesign methods.
<br>
Note that each change of input variables into the model can significantly influence the results. <br>
With that in mind, any alteration in the input parameters should be justified by actionable decisions that will influence the reliability growth program. For example, increasing the average effectiveness factor value should be done only when there is proof that the program will pursue a different, more effective path in terms of addressing fixes.
<br>
<br>

Latest revision as of 01:40, 27 August 2012

Redirect to:

Retrieved from "https://www.reliawiki.com/index.php?title=Template:Crow_extended_model_rga&oldid=33933"