8. Assessing Product Reliability
8.1. Introduction
8.1.9. How can you model reliability growth?

8.1.9.1.

NHPP power law

If the Power Law applies, Repair Rates improve over time according to the formula. The exponent  lies between 0 and 1 and is called the reliability growth slope

This repairable system model was described in  Section 8.1.7.2. The expected number of failures by time t has the form M(t) = at^b and the repair rate has the form m(t) = abt^b-1. This will model improvement when 0 < b < 1, with larger improvements coming when b is smaller. As we will see in the next section on Duane Plotting, it is convenient to define = 1 - b and  =  ab, and write the repair rate as 

m(t) = 

Again we have improvement when 0 <  < 1, with larger improvement coming from larger values of .  is known as the Duane Plot slope or the reliability improvement Growth Slope. 

In terms of the original parameters for M(t), we have 

Use of the Power Law model for reliability growth test data generally assumes the following:

1. While the test is ongoing, system improvements are introduced that produce continual improvements in the rate of system repair. 

2. Over a long enough period of time the effect of these improvements can be modeled adequately by the continuous polynomial repair rate improvement model .

When an improvement test ends, the MTBF stays constant at its last achieved value 

3. When the improvement test ends at test time T and no further improvement actions are ongoing, the repair rate has been reduced to . The repair rate remains constant from then on at this new (improved) level. 

Assumption 3 means that when the test ends, the HPP constant repair rate model takes over and the MTBF for the system from then on is the reciprocal of the final repair rate or . If we estimate the expected number of failures up to time T by the actual number observed, the estimated MTBF at the end of a reliability test (following the Power Law) is: 

with T denoting the test time, r  is the total number of test failures and  is the reliability growth slope. A formula for estimating  from system failure times is given in the Analysis Section for the Power Law model.

Simulated Data Example

Simulating NHPP Power Law Data

Step 1: User inputs the positive constants a and b. 

Step 2: Simulate a vector of n uniform (0,1) random numbers. Call these U₁, U₂, U₃, . . . U_n. 

Step 3: Calculate Y₁ = {-1/a * ln U₁} ** 1/b

Step i: Calculate Y_i = {(Y_i-1 ** b) -1/a * ln U_i}**1/b  for i = 2, . . ., n 

The n numbers Y₁, Y₂, . . ., Y_n are the desired repair times simulated from an NHPP Power Law process with parameters a, b (or  = 1 - b and  = ab). 

The Dataplot Macro powersim.dp will ask the user to input the number N of repair times desired and the parameters a and b. The program will output the N simulated repair times and a plot of these repair times. 

Example

Below powersim.dp is used to generate 13 random repair times from the NHPP Power Law process with a = .2 and b = .4.

CALL powersim.dp 

Enter number N of simulated repair times desired
13
Enter value for shape parameter a (a > 0)
.2
Enter value for shape parameter b (b > 0)
.4 

FAILNUM	FAILTIME
1	26
2	182
3	321
4	728
5	896
6	1268
7	1507
8	2325
9	3427
10	11871
11	11978
12	13562
13	15053