General Examples & Case Studies: Difference between revisions

From ReliaWiki
Jump to navigation Jump to search
Line 132: Line 132:
<br>
<br>


==Arrhenius Parameter Interpretation Example==
==Arrhenius Parameters Interpretation Example==
A mechanical component was put into an accelerated test with temperature as the accelerated stress. The following times-to-failure were observed.
A mechanical component was put into an accelerated test with temperature as the accelerated stress. The following times-to-failure were observed.



Revision as of 17:19, 21 March 2012

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/accelerated_life_testing_data_analysis

Chapter 12: General Examples & Case Studies


ALTAbox.png

Chapter 12  
General Examples & Case Studies  

Synthesis-icon.png

Available Software:
ALTA

Examples icon.png

More Resources:
ALTA Examples

In this chapter we present some application examples utilizing ReliaSoft's ALTA. We assume that you have previously consulted the ALTA User's Guide and familiarized yourself with the software.


Paper Clip Example

To illustrate the principles behind accelerated testing, consider the following simple example which involves a paper clip and can be easily and independently performed by the reader. The objective was to determine the mean number of cycles-to-failure of a given paper clip. The use cycles were assumed to be at a 45º bend. The acceleration stress was determined to be the angle to which we bend the clips, thus two accelerated bend stresses of 90º and 180º were used. The paper clips were tested using the following procedure for the 90º bend. A similar procedure was also used for the 180º and 45º test.


1) Open the Paper Clip.

Openclip.gif


2) With one hand, hold the clip by the longer, outer loop.

3) With the thumb and forefinger of the other hand, grasp the smaller, inner loop.

4) Pull the smaller, inner loop out and down 90 degrees so that a right angle is formed as shown.

5) Close the Paper Clip.

Closeclip.gif


6) With one hand, continue to hold the clip by the longer, outer loop.

7) With the thumb and forefinger of the other hand, grasp the smaller, inner loop.

8) Push the smaller inner loop up and in 90 degrees so that the smaller loop is returned to the original upright position in line with the larger, outer loop as shown.

9) This completes one cycle.

10) Repeat until the paper clip breaks. Count and record the cycles-to-failure for each clip.


At this point, the reader must note that the paper clips used in this example were Jumbo paper clips capable of repeated bending. Different paper clips will yield different results. Additionally, and so that no other stresses are imposed, caution must be taken to assure that the rate at which the paper clips are cycled remains the same across the experiment.


For the experiment, a sample of six paper clips was tested to failure at both 90º and 180º bends. A base test sample of six paper clips was tested at a 45º bend (the assumed use stress level) to confirm the analysis. The cycles-to-failure are given next.


Cycles-to-Failure at [math]\displaystyle{ 90^\circ }[/math] 90degree.gif


16, 17, 18, 21, 22, 23 cycles.


Cycles-to-Failure at [math]\displaystyle{ 180^\circ }[/math] 180degree.gif


4, 4, 5, 6, 6, 8 cycles.


Cycles-to-Failure at [math]\displaystyle{ 45^\circ }[/math] 45degree.gif


58, 63, 65, 72, 78, 86 cycles.


The accelerated test data set was then analyzed in ALTA, assuming a lognormal life distribution (fatigue) and an inverse power law relationship (non-thermal) for the life-stress model. The analysis and some of the results are shown in the next figures.


The accelerated test data analyzed in ALTA.



Lognormal probabiliyt plot of both stress levels from ALTA.



The resulting Acceleration Factor vs. Stress plot from ALTA.


The MTTF at 45 degrees is found to be 71.375467 cycles as shown next.


The Mean Life at 45 degrees bend.


The next two figures show the analysis of the base data in Weibull++ and the base MTTF estimate. In this case, our accelerated test correctly predicted the MTTF as verified by our base test.


The base [math]\displaystyle{ 45^\circ }[/math] data analyzed in ReliaSoft's Weibull++, utilizing a lognormal distribution shown on lognormal probability paper.


Calculated Mean Life from base test.


Temperature-Voltage Test Example


Twelve electronic devices were put into a continuous accelerated test. The accelerated stresses were temperature and voltage, with use level conditions of 328K and 2V respectively. The data set obtained is shown in the table below:

11 2ex.gif


Do the following:

1) Using the T-NT Weibull model, analyze the data in ALTA and determine the MTTF and B(10) life for these devices at use level. Determine the upper and lower 90% 2-sided confidence intervals on the results.

2) Examine the effects of each stress on life.


Solution

1. The data set was analyzed in ALTA and the following MTTF and [math]\displaystyle{ B(10) }[/math] life results were obtained: (The results are shown in ALTA's QCP in the following figures.)


Pv ex2 1.gif


Pv ex2 2.gif




2. The next two figures examine the effects of each stress on life.


Specifically, the next figure shows the Life vs. Voltage plot with temperature held constant at 328K.


The effects of voltage on life, with temperature held constant at 328K.


The next figure shows the Life vs. Temperature plot with voltage held constant at 2V.


The effects of temperature on life, with voltage held constant at 2V.



The next figure shows a 3D plot of reliability (at a constant time) versus both stresses. One can observe that reliability declines slightly faster with the change in temperature than with the change in voltage. Note that in the following plot, Stress 1 refers to temperature and Stress 2 refers to voltage.


The combined effects of voltage and temperature on the reliability, as plotted in ALTA.



Arrhenius Parameters Interpretation Example

A mechanical component was put into an accelerated test with temperature as the accelerated stress. The following times-to-failure were observed.


343 K 363 K 383 K
266.66 618.54 351.12
430.09 666.72 355.1
570.45 724.4 672.69
890.42 950.89 923.35
1046.65 1148.4 948.22
1158.14 1202.94 1277.04
1396.01 1492.56 1538.81
1918.38 1619.59 2020.34
2028.86 2592.29 2099.03
2785.58 3596.85 2173.04


1) Determine the parameters of the Arrhenius-Weibull model.

2) What is your observation?



Solution

The parameters of the Arrhenius Weibull model were estimated using ALTA with the following results:


[math]\displaystyle{ \beta =1.771460,\text{ }B=86.183591,\text{ }C=1170.423770. }[/math]


A small value for [math]\displaystyle{ B }[/math] was estimated. The following observations can then be made:


• Life is not accelerated with temperature, or
• the stress increments were not sufficient, or
• the test stresses were well within the "specification limits" for the product (see discussion here).


A small value for [math]\displaystyle{ B }[/math] is not the only indication for this behavior. One can also observe from the data that at all three stress levels, the times-to-failure are within the same range. Another way to observe this is by looking at the Arrhenius plot. The scale parameter, [math]\displaystyle{ \eta }[/math] , and the mean life are plotted next.

Eta and Mean Life vs. Stress.


It can be seen that life ( [math]\displaystyle{ \eta }[/math] , and the mean life) is almost invariant with stress.

IPL-Lognormal Example

A tensile component of a landing gear was put through an accelerated reliability test to determine whether the life goal would be achieved under the designed-in load. Fifteen units, N=15, were tested at three different shock loads. The component was designed for a peak shock load of 50 kips with an estimated return of 10% of the population by 10,000 landings. Using the inverse power law lognormal model, determine whether the design life was met.

Failure, in Cycles Shock Load, kips
1855 73
2158 73
2425 73
2736 73
3206 73
614 95
770 95
781 95
890 95
1055 95
152 128
162 128
173 128
209 128
216 128


Solution

The data were entered into ALTA, and the following estimates were obtained for the parameters of the IPL-lognormal model:


[math]\displaystyle{ \begin{align} Std=\ & 0.179080 \\ K=\ & 8.953055E-13 \\ n=\ & 4.638882 \end{align} }[/math]


The probability plot at each test stress level (73 kips, 95 kips and 123 kips) was then obtained, as shown next.


Probability plots at different shock loads.


In this plot, it can be seen that there is a good agreement between the data and the fitted model. The probability plot at a use stress level of 50 kips is shown next.


The probability plot at a use stress level of 50 kips.


Using the Use Level Probability plot, the cycles-to failure for a 10% probability of failure can be estimated. From the plot, this time is estimated to be [math]\displaystyle{ T(Q=0.10=10%)\approxeq 12,000 }[/math] landings. Another way to obtain life information is by using a Life vs. Stress plot. Once a Life vs. Stress plot is selected, the 10% unreliability (probability of failure) line must be plotted. To do this, click Life Lines... in the plot panel and and select the 10% Unreliability box (or enter 10 in one of the boxes and select the box) in the Specify Life Lines window, as shown next.


Ipl example specify life.gif


In the following plot, the 10% unreliability line is plotted (the first line from the left). For a stress of 50 psi (X-axis) and for the 10% unreliability line, the cycles-to-failure can be obtained by reading the value on the Y-axis. Again, [math]\displaystyle{ T(0.1)\approxeq 12,000 }[/math] landings.


Median Life line (right) and 10% unreliability line (left) vs. stress.


However, a more accurate way to obtain such information is by using the Quick Calculation Pad (QCP) in ALTA. Using the QCP, the life for a 10% probability of failure at 50 kips is estimated to be 11,669.1 cycles (or landings), as shown next.


Ipl example2 qcp.gif


The design criterion of 10,000 landings is met.

IPL-Lognormal Confidence Bounds Example

In the previous example, it was found that the 10,000 landings life criterion for a 90% reliability was met. However, the estimate of 11,669.1 cycles (or landings) is very close to the requirement. Additionally, this estimate was obtained at the 50% confidence level. In other words, 50% of the time life will be greater than 11,669.1 landings, and 50% of the time life will be less. A lower confidence level is thus crucial before any decisions are made. Repeat the calculation of the previous example, and this time include a 90% lower 1-sided confidence bound on the estimation.

Solution


Using ALTA's QCP, the 90% lower 1-sided confidence bound on the 11,669.1 cycles-to-failure can be obtained.

Ipl example2 qcp.gif



The 90% lower 1-sided confidence bound is estimated to be 9,680.84 landings. Thus the 10,000 landings criterion is not quite met.


IPL-Weibull Example

ACME manufacturing has implemented an accelerated testing program for their new design. A total of 40 units were tested at four different pressure levels. The operating stress level is 170 psi.


Stress Level, psi 220 psi 230 psi 240 psi 250 psi
165 93 72 26
177 106 73 44
238 156 99 63
290 170 124 68
Times-to-failure, hr 320 185 134 69
340 214 150 72
341 220 182 77
380 236 186 96
449 252 190 131
544 288 228 140


Do the following:


1) Determine the parameters of the inverse power law Weibull model.


2) Obtain the use level probability plot with 90% 2-sided confidence bounds on time.


3) Obtain the Life vs. Stress plot with 90% 2-sided confidence bounds.


Solution

1) The parameters of the IPL-Weibull model are estimated to be:


[math]\displaystyle{ \begin{align} \beta =\ & 3.009236 \\ K=\ & 3.267923E-27 \\ n=\ & 10.218104 \end{align} }[/math]


2) The use level probability plot is shown next.

The probability plot at a use stress level.


The confidence bounds on time are the Type 1 bounds (Time Bounds) in ALTA.


The probability plot at a use stress level with 90% Type I confidence bounds.


3) Similarly, the Life vs. Stress plot with the confidence bounds can be obtained.


Arrhenius-Weibull Example

Twenty-seven circuit boards were tested with temperature as the accelerated stress (or stimuli). The goal was to estimate the reliability after 10 years at 100ºC. The circuit boards were tested and the following data obtained:


Test Chamber Temperature 463K
# of units Time, hr Observation
1 2,403 Failure
1 2,668 Unrelated Failure
1 3,669 Test Fixture Failure
1 3,863 Failure
2 4,400 Test Fixture Failure
2 4,767 Failure
1 5,219 Failure
4 5,276 Power Surge Failure
2 7,517 Failure
1 7,840 Failure
2 8,025 Failure
1 8,571 Removed - Test Terminated


Test Chamber Temperature 488K
# of units Time, hr Observation
2 346 Power Regulator Failure
1 1,416 Failure
1 2,197 Unrelated Failure
1 2,533 Failure
1 2,630 Failure
1 2,701 Test Fixture Failure
1 3,000 Failure
1 3,489 Failure
1 6,720 Removed - Test Terminated



The data were entered in ALTA and analyzed utilizing an Arrhenius-Weibull model. Note that all units that did not fail in the test but failed for unrelated reasons were entered as suspensions. The reliability after 10 years at 100ºC is 45.67% as shown in the next figure.

Reliability plot at use stress level


Arrhenius-Weibull Example

An electronic component was redesigned and was tested to failure at three different temperatures. Six units were tested at each stress level. At the 406K stress level however, a unit was removed from the test due to a test equipment failure which led to a failure of the component. A warranty time of one year is to be given, with an expected return of 10% of the population. The times-to-failure and test temperatures are given next:



11 8ex.gif



The operating temperature is 356K. Using the Arrhenius-Weibull model, determine the following:


1) Should the first failure at 406K be included in the analysis?


2) Determine the warranty time for a 90% reliability.


3) Determine the 90% lower confidence limit on the warranty time.


4) Is the warranty requirement met? If not, what steps should be taken?


5) Repeat the analysis with the unrelated failure included. Is there any difference?


6) If the unrelated failure occurred at 500 hr, should it be included in the analysis?


Solution

1) Since the failure occurred at the very beginning of the test and for an unrelated reason, it can be omitted from the analysis. If it is included it should be treated as a suspension and not as a failure.


2) The first failure at 406K was neglected and the data were analyzed using ALTA. The following parameters were obtained:


[math]\displaystyle{ \begin{align} \beta =\ & 2.965820 \\ B=\ & 10679.567542 \\ C=\ & 2.396615\cdot {{10}^{-9}} \end{align} }[/math]


The use level probability plot (at 356K) can then be obtained. The warranty time for a reliability of 90% (or an unreliability of 10%) can be estimated from this plot as shown next.


The probability plot at use stress level.


This estimate can also be obtained from the Arrhenius plot (a Life vs. Stress plot). The 10th percentile (time for a reliability of 90%) is plotted versus stress. This type of plot is useful because a time for a given reliability can be determined for different stress levels.


10% unreliability vs stress.


A more accurate way to determine the warranty time would be to use ALTA's Quick Calculation Pad (QCP). By selecting the Reliable Life option from in the QCP and entering 356 for the temperature and 0.9 for the required reliability, a reliable life of 11,977.792719 hr can be determined, as shown next:


Ex1 qcp.gif


3) The warranty time for a 90% reliability was estimated to be approximately 12,000 hr. This is above the 1 year (8,760 hr) requirement. However, this is an estimate at the 50% confidence level. In other words, 50% of the time life will be greater than 12,000 hr and 50% of the time life will be less. A known confidence level is therefore crucial before any decisions are made. Using ALTA, confidence bounds can be plotted on both probability and Arrhenius plots. In the following use level probability plot, the 90% lower confidence level (LCL) is plotted. Note that percentile bounds are Type 1 confidence bounds in ALTA.


The use level probability plot with the Type I 90% lower bound.


An estimated 4,437 hr warranty time at a 90% lower confidence level was obtained from the use level probability plot. This means that 90% of the time, life will be greater than this value. In other words, a life of 4,437 hr is a bounding value for the warranty.


The Arrhenius plot with the 90% lower confidence level is shown next.


Life vs Stress plot with the 90% lower confidence level.


Using the QCP and specifying a 90% lower confidence level, a warranty time of 4436.5 hr is estimated, as shown next.

Ex1 qcp2.gif


4) The warranty time for this component is estimated to be 4,436.5 hr at a 90% lower confidence bound. This is much less than the 1 year warranty time required (almost 6 months). Thus the desired warranty is not met. In this case, the following four options are available:


a) Redesign.
b) Reduce the confidence level.
c) Change the warranty policy.
d) Test additional units at stress levels closer to the use level.


5) Including the unrelated failure of 0.3 hr at 406 K (by treating it as a suspension at that time), the following results are obtained:


[math]\displaystyle{ \begin{align} \beta =\ & 2.965804 \\ B=\ & 10679.566625 \\ C=\ & 2.396619\cdot {{10}^{-9}} \end{align} }[/math]


These results are essentially identical to the ones with the unrelated failure excluded. A small difference can be seen only if more significant digits are considered. The warranty time with the 90% lower 1-sided confidence bound was estimated to be:


[math]\displaystyle{ \begin{align} T=\ & 11977.729145\text{ hr} \\ {{T}_{L}}=\ & 4436.462084\text{ hr} \end{align} }[/math]


Again, the difference is negligible. This is due to the very early time at which this unrelated failure occurred.


6) The analysis is repeated treating the unrelated failure at 500 hr as a suspension, with the following results:


[math]\displaystyle{ \begin{align} \beta =\ & 3.022714 \\ B=\ & 10959.523204 \\ C=\ & 1.23808\cdot {{10}^{-9}} \end{align} }[/math]


In this case, the results are very different. The warranty time with the 90% lower 1-sided confidence bound is estimated to be:


[math]\displaystyle{ \begin{align} T=\ & 13780.208409\text{ hr} \\ {{T}_{L}}=\ & 5303.672050\text{ hr} \end{align} }[/math]


It can be seen that in this case, it would be a mistake to neglect the unrelated failure. By neglecting this failure, we would actually underestimate the warranty time. The important observation in this example is that every piece of life information is crucial. In other words, unrelated failures also provide information about the life of the product. An unrelated failure occurring at 500 hr indicates that the product has survived for that period of time under the particular stress level, thus neglecting it would be a mistake. On the other hand it would also be a mistake to treat this data point as a failure, since the failure was caused by a test equipment failure.

Example 9

An electronic component was subjected to a voltage stress, starting at 2V (use stress level) and increased to 7V in stepwise increments. The following steps, in hours, were used to apply stress to the products under test: 0 to 250, 2V; 250 to 350, 3V; 350 to 370, 4V; 370 to 380, 5V; 380 to 390, 6V; and 390 to 400, 7V. This profile is represented graphically next.

ALTA14.20.gif


The objective of this test was to determine the [math]\displaystyle{ B10 }[/math] life of these components at the normal use stress level of 2V.
In this experiment, the overall test time was 385 hrs. If the test had been performed at use conditions, one would expect the test duration to be approximately 1700 hrs if the test were run until all units failed.
Eleven units were available for the test. All eleven units were tested using the same stress profile. Units that failed were removed from the test and their total time on test recorded. The following times-to-failure were observed in the test, in hours: 280, 310, 330, 352, 360, 366, 371, 374, 378, 381 and 385.


Solution to Example 9

After creating a Standard Folio for a single stress (voltage), the analyst added a Stress Profile under Tools in the Project Manager to create the stress profile for this experiment. In ALTA's stress profile, you can define any stress profile in segments. The stress applied during these segments can be a constant value (as is the case in this step example) or defined as a function of a time variable (t). The stress profile for this analysis is displayed next.

Chp14-65.gif


The voltage test profile as in ALTA.


The times-to-failure data are entered in the ALTA Standard Folio with the following analysis options selected: cumulative damage power life-stress relationship and Weibull for the underlying life distribution. The ALTA Standard Folio with the data entered and the estimated parameters is shown next. Notice that the voltage step-stress profile has been assigned to each data point.


Chp14-66.gif


The [math]\displaystyle{ B10 }[/math] life at the 2V use stress level, can be calculated with ALTA's Quick Calculation Pad (QCP), as shown next.


Ch12-Ex9-4.gif


CD-Weibull Example

Consider a test in which multiple stresses were applied simultaneously to a particular automotive part in order to precipitate failures more quickly than they would occur under normal use conditions. The engineers responsible for the test were able to quantify the combination of applied stresses in terms of a "percentage stress" as compared to typical stress levels (or assumed field conditions). In this scenario, the typical stress (field or use stress) was defined as 100% and any combination of the test stresses was quantified as a percentage over the typical stress. For example, if the combination of stresses on test was determined to be two times higher than typical conditions, then the stress on test was said to be at 200%.


The test was set up and run as a step-stress test (i.e. the stresses were increased in a stepwise fashion) and the time on test was measured in hours. The step-stress profile used was as follows: until 200 hours, the equivalent applied stress was 125%; from 200 to 300 hrs, it was 175%; from 300 to 350 hrs, it was 200% and from 350 to 375 hrs, it was 250%. The test was terminated after 375 hours and any units that were still running after that point were right-censored (suspended).


Additionally, and based on prior analysis/knowledge, the engineers also stated that each hour on test under normal use conditions (i.e., at 100% stress measure) was equivalent to approximately 100 miles of normal driving. The test was conducted and the following times-to-failure and times-to-suspension under the stated step-stress profile were observed (note that + indicates a non-failed unit, i.e., suspension): 252, 280, 320, 328, 335, 354, 361, 362, 368, 375+, 375+, 375+.


After performing failure analysis on the failed parts, it was determined that the failure that occurred at 328 hrs was due to mechanisms other than the ones considered. That data point was therefore identified as a suspension in the currently analysis. The modified data set for this analysis was: 252, 280, 320, 328+, 335, 354, 361, 362, 368, 375+, 375+, 375+.


The test objective was to estimate the [math]\displaystyle{ B1 }[/math] life for the part (i.e. time at which reliability is equal to 99%) at the typical operating conditions (i.e., Stress=100%), in miles.


Solution

Utilizing ALTA, the analyst first created a Standard Folio for non-grouped time-to-failure and time-to-suspension data and added a Stress Profile under Tools in the Project Manager to define the stress profile, as shown next.

Autoexampleprofile.gif


Once the profile was defined, the analyst entered the observed times, their state (i.e., failed F or non-failed S) and a reference to the profile used in the ALTA Standard Folio. The analyst selected the cumulative damage life-stress relationship (to use a time-varying stress) based on a power model (since the effect of the stress was deemed to be mechanical and more appropriately modeled by a power function). The Weibull distribution was selected as the underlying life distribution. Additionally, note that the use stress was set to 100 and all results were then extrapolated to the typical stress level. The Standard Folio and the selected model is shown next.

Autoexamplefolio.gif



The estimated parameters are shown next.

Autoexampleresults.gif



The last part remaining was to determine the [math]\displaystyle{ B1 }[/math] life at the part's use stress level. Using the QCP, the [math]\displaystyle{ B1 }[/math] life was found to be 658 hours, as shown next.

Autoexampleqcp.gif


Based on the given multiplier, the [math]\displaystyle{ B1 }[/math] life in miles would then be 658 test-hr*100 (miles/test-hr)= 65,800 miles.


The [math]\displaystyle{ B1 }[/math] life can also be obtained from the use level probability plot, as shown next.

Autoexampleplot.gif


Arrhenius-Weibull Interval Data Example

A sample of and electronic device was tested in two testing chambers set at different temperatures and monitored on a weekly basis. The test was performed continuously for 15 weeks. 50 units were put in each chamber. The goal of the test is to estimate the 90% lower bound on the [math]\displaystyle{ B10 }[/math] life (in hours) of the electronic device. The following figure represents the test results, it shows the test results entered in ALTA in an interval data format. The normal use condition is assumed to be 300K.


Chapter14 80.gif


Solution

Assuming that the failure times of the devices follow a Weibull distribution and that the Arrhenius model is an adequate life-stress relationship, we obtain the following estimated model parameters:


[math]\displaystyle{ \begin{align} \beta =\ & 3.425937 \\ B=\ & 1488.412639 \\ C=\ & 24.059577 \end{align} }[/math]


The next figure show the probability lines at the different stress levels along with the extrapolated use level probability line.


Probability lines at the different test stress levels with the extrapolated use level probability line.


The following is the life versus stress plot showing the [math]\displaystyle{ B10 }[/math] life line.


Life vs. Stress with the B10 life line (bottom).


The 90% lower bound on the [math]\displaystyle{ B10 }[/math] life can be estimated as follows.

Intervalb10qcp.gif


[math]\displaystyle{ }[/math]