Reliability Engineering
Reliability Engineering
Engineering
Therory and Practice
Hamid Bazargan_Harandi
January 2023
2
To my parents
&
C H A P T E R T I T L E S
F O R E W O R D ...................................................... 11
Nomenclature .............................................................. 13
References ................................................449
T A B L E S .........................................455
About the author …………………………………………….476
4
5 Reliabilty Engineering
C O N T E N T S
F O R E W O R D ...................................................... 11
Nomenclature .............................................................. 13
FO REWORD
This book, whose Persian version written by the same author
is published by our univerisity, is the outcome of teaching a
course on reliability for several years to graduate students
using many books especially the book written by Dr Kapur and
Lamberson (abbreviated by K&L throughout the book) . The
main prerequisite for understanding the materials of the book is
probability. As evident from the chapter titles, the book
introduces readers with reliability and availability of products. I
would like to the students who helped the author in some phases
of editing.
At the end of this foreword a software and a symposium are
introduced.
Software
ReliaSoft's reliability software tools facilitate a set of
reliability engineering modeling and analysis techniques.
(downloable from http://www.reliasoft.com)
A symposium: RAMS associated to IEEE
Reliability and Maintainability www.rams.org
The proceedings are available from
http://ieeexplore.ieee.org/xpl/mostRecentIssue.jsp?punumber=6516162
It is suggested to the readers, especially those working in
industry, to read books on design for reliability after reading this
book.
Thanks God for making me successful to present this work
which I hope to be useful in both academic and industrial
environments.
The author would be pleased if the readers write him about
any kind of deficiencies in the book.
Hamid Bazargan
Jan 2023
College of Engineering,
Shahid Bahonar University of Kerman, Iran
[email protected]
The wise is one who puts everything
in its right place
12
Nomenclature Nomenclature
Symbols Description
Nomenclature Nomenclature
( ) The reliability function of The failure rate of
ith unit component a
The failure rate of
component b
Average load
U,V,W,X Average stength
The upper specification The standard devaition of
U level of a dimention of a apopulation
product Standard devation of load
The difference of U al L sistribution
( − ) standard devation of
in a part
strength distribution
δ lower bound of strength
( − ) The difference of U al L in
the assembled unit
σ part the standard devaition of a
V Coefficient of variation of measurement in a part
safety factor(SF)
In AGREE Method: the standard devaition of a
The importance factor measurement
of the ith subsystem
sum in the assembly
Chapter 1
Introduction
&
Basic
Concepts
Chap 1 Introduction and Basic Concepts 16
1
Introduction and Basic Concepts
1.1 Introduction
1-2 Failure
1-3 Reliability
1
From (https://www.techslang.com/definition/what-is-reliability-engineering/)
19 Reliabilty Engineering
1-3-1 abbreviations
Note that TBF is equal to the sum of TTR and TTF (Fig.1.1).
TBF=TTR+TTF (1-1-1)
MTBF=MTTR+MTTF (1-1-2)
WT
M TˆBF (1-2)
N
∞
R(t) = Pr(X > ) = Pr(X ≥ t) = ∫ f(x)dx − 1 − F(x) (1-3)
where
( )= ( > ) , k = 1, 2, . . .. (1-4-1)
( ) = ∑∞ ( ), (1-4-2)
( ) = (1 − ) , k = 1, 2, . . .. (1-5)
( )
( )≅ . (1-6)
∞
( )=∫ ∞
( ) . (1-7)
∞
( ) = ∫ [1 − ( )] −∫ ∞
( ) (1-8)
∞
( ) = ∫ [1 − ( )] (1-9)
or from :
( )=∫ ( ) . (1-10)
Example 1-1
Solution
( )= = , ( )= ≅ 0.37
Example 1-2
Solution
Let
Min (X,T)+Max(X,T)=X+T⟹
E(X + T) = E(X) + T
then
λ λ λ
E{Max(X, T)}=( T + λ) × e +T 1−e = T+ λe
Chap 1 Introduction and Basic Concepts 26
= + T − (T + e λ ) ⇒ E(Y) = (1 − e λ )
λ λ λ
Pr(t1 X t2 ) ( )− ( )= ( ) − ( ).
( )− ( )
( < < | > ) ( ) ( )− ( )
= =
− − ( − ) ( )
( ) ( ∆ ) ( ) ( )
ℎ( ) = lim∆ → (∆ )× ( )
= ( )
− = ( )
, (1-11)
Then for devices with continuous lifetime having pdf f(t) and
reliability(survivor) function R(t), the hazard function is defined
as the ratio of the probability density function to the survivor
function.
( ) ′( )
ℎ( ) = ( )
= ( )
. (1-12)
ℎ( ) ≥ 0 ≥0 (1-13-1)
∫ ℎ( ) = ∞. (1-13-2)
29 Reliabilty Engineering
P (k) P (k)
h(k) = = . (14 − 1)
∑∞ P (k) R(k)
( ) ( )
h(k)= =∑ ( )
= .
( ) ! ∑∞
!
( , )
in MATLAB: h(k)= .
( , )
Chap 1 Introduction and Basic Concepts 30
h(t) is known:
∫ ( )
( )= (1-14-1)
( ) = ℎ( ) ( ) (1-14-2)
( )=1− ( ) (1-14-3)
f(t) is known:
F(t) is known:
( )
( )= ( ) ( )=1− ( ) ℎ( ) = ( )
R(t) known
-R'(t)
h(t)= ( ) = ℎ( ) ( ) ( )=1− ( )
R(t)
1
f a b (x ) b
f (x ) (1-18)
f (x )dx
a
Chap 1 Introduction and Basic Concepts 32
∞
( )= ( > )=∫ = ≥ 0 (1-19-2)
f (t ) 1
h (t ) (1-19-3)
R (t )
1
some softwares such as ARENA could determine the best distributions that
fit a data set(e.g in ARENA tools-input analyzer- new-file data file- use
existing- fit all)
33 Reliabilty Engineering
0< < .
x 2
f (x )
1
e 2 2 , x (1-20-1)
2
The rate function is:
1 t 2 t
1
h (t ) exp 1 Z (1-20-2)
2 2 2
Chap 1 Introduction and Basic Concepts 34
where
Z is the cumulative distribution function(CDF) of standard normal.
ℎ(
=1
=2
= 0.2
= 0.3
t=TTF
Fig. 1.1 shows the rate function of 2 normal distributions.
Example 1.3
19000 − 20000
R(19000) = Pr( > 19000) = Pr Z > =
2000
Pr(Z>0.5)=1-Pr(Z<-0.5)=0.69146=69.15%
35 Reliabilty Engineering
1 - t-200002
f(t)= e 2(2000 ) f(19000)=0.000176
2000 2π
f(19000) 0.000176
h(19000)= = =0.000245 failures/cycle
R(19000) 0.69146
μ
f (t ) = e σ , t ≥ 0 , σ > 0 . < μ < ∞ (1-21-1)
σ√ π
where
∞
1 − −
= = > =1−Φ ( )
√2
Note that
-The probabilities of this distribution is not calculated in the
same manner which is done in classical normal distributions.
- The mean and variance of this distribution does not equal
and 2 . The mean is
Chap 1 Introduction and Basic Concepts 36
2
1 μ
-
E X e 2 σ (1-21-2)
a
This truncated distribution has an increasing rate function.
Figure 1-2 shows the function for typical one plotted with the
following MATLAB commands:
mu=3;sigma=0.1;t=0:.01:10;f=(normpdf(t,mu,sigma))/(normcdf(mu/si
gma))./(1-normcdf(t,mu,sigma));plot(t,f);
( )= −1 (1-22-3)
3)
median= (1-22-4)
4)
Figure 1.3 shows the pdf of 2 sample lognormal
distribution
distributions.
1-8-4-1
1 Calculation of the parameters(μ,σ)
parameters( of lognormal
distribution from the mean and variance
( )
= ln ( )
+1 , (1-22
22-5)
= ln ( ) − . (1-22
22-6)
1-8-4-2
2 The relationship between lognormal (μ,σ)
( ) and normal
(μ,σ)) distributions
Chap 1 Introduction and Basic Concepts 38
Fig 1-5 The hazard function for some lognormal variable with μ = 1
(the horizontal axis is the time and the vertical is the failure function)
Example 1-4
Solution
( )= ( = 150) ≅ 0.0027
√ ,
41 Reliabilty Engineering
( ) /
ℎ(150) = ( )
= = 0.0053 1 =
.
53 failures 10000 .
≤ ≤
( )= (1-23-1)
0 .
In this distribution:
( )= (1-23-2)
( )= (1-23-3)
( )
ℎ( ) = ( )
= ≤ ≤ (1-23-4)
Chap 1 Introduction and Basic Concepts 42
Figures 1-6 and 1-7 shows the density function f (x) and the
hazard function ℎ( ).
Fig 1.7 The hazard function of Fig 1.7 The density function of
a uniform distribution a uniform distribution
C
C 1 tA
C t A
B
f t e tA (1-24-1)
B B
( )=1− (1-24-2)
43 Reliabilty Engineering
( )= (1-24-3)
C 1
( ) C t A
ℎ( ) = = (1-24-4)
( ) B B
1
E (T ) A B (1 ) (1-24-5)
C
2 1
Var(T) B2 Γ(1+ )-B2 Γ(1+ ) (1-24-6)
C C
1
median A B (ln 2) C (1-24-7)
In weibull distribution:
Fig. 1-8 Plot of the Weibull distribution for scale parameter B=1 and
five values of shape parameter (extracted from Grant&Leavenworth ,1988
page605)
∑ ∑
= ∑
− (1-23-8)
∑
= (1-23-9)
45 Reliabilty Engineering
Example 1-5
Write a MALAB code to estimate the parameters B and C
of a Weibull distribution from which the following random
sample is at hand:
X=[113.0634 49.5432 53.4872 93.7147 74.0594
114.3216 97.1033 61.5069 74.7216 52.8807];
Furthermore estimate the parameters with wblfit MATLAB
command.
Solution
%Sample X=[X(1)......X(n)]
X=[113.0634 49.5432 53.4872 93.7147 74.0594
114.3216 97.1033 61.5069 74.7216 52.8807];
for C=.01:0.001:40
for I=1:length(X)
LNX(I)=log(X(I));
XIC(I)=X(I)^C;XICLNX(I)=XIC(I)*LNX(I);
end
A= C-(sum(XICLNX)/sum(XIC)-sum(LNX)/length(X))^(-1);
if ( abs(A)<= 0.001 ) C1=C; disp(sprintf('C= %6.4f ', C1))
end
end
B=(sum(X.^C1)/(length(X)))^(1/C1);
disp(sprintf('B= %6.4f ', B))
with MATLAB command wblfit:
Chap 1 Introduction and Basic Concepts 46
>>wblfit(X)
ans = 87.1543 3.7149
x 1e x x 0
f x (1-24-1)
0 x 0
Furthermore:
1 Extreme Value
2
Fisher-Tippet I or Gumbel distribution
47 Reliabilty Engineering
E (X )
(1-24-2)
V ar ( X )
(1-24-3)
2
M GF (t ) (
)
(1-24-4)
t
=
Γ( + )
+ >0
(1-25)
Γ( )
Fig. 1-10
1 pdf, h(t),R(t) of some distributions
Chap 1 Introduction and Basic Concepts 50
Fig 1-11-2
2 Some variations of bathtub curve(Kuo&Zuo,2003)
cu ve(Kuo&Zuo,2003)
Fig 1-11-3
3 A variation of bathtub curve for some mechanical devices
(Ireson,1995 page18-2)
Chap 1 Introduction and Basic Concepts 52
Fig 1-11
11-4 Another variation of Bathtub failure rate function
(Nahmias, 2004, Fig.12-4)
Figure 1-11-5
5 shows different forms of bathtub failure curve due to
different levels of stress on some mechanical devices.
Fig1-11--5
5 Effects of stress levels on mechanical failure rates
(Ireson,199
n,1995 Fig18-1, Stamatis,2010 Fig. 6-3
3 page163 )
53 Reliabilty Engineering
For each part of the bath tub curve different hazard function
is appropriate. So next some kinds of failure rate function are
considered(K&L page 28) .
a- Constant hazard( failure rate) function
If the failure rate function is ℎ( ) = i.e. is constant and
does not depend on time ,according to the relationship
between the hazard function and the density function of
lifetime(TTF);
∫ ( ) ∫
( ) = ℎ( ) = = >0
( )
( )= = =
ℎ( )
∫ ( )
( )= = (1-26)
Chap 1 Introduction and Basic Concepts 54
∫ ( )
( ) = ℎ( ) ⟹
∫
( )= = ≥0 (1-27-1)
( )= (1-27-2)
For the case that the filature rate of a device is like a bathtub except
the first and last part are linear, the hazard function is as follows(K&L
page 29):
∫ ( ) C t
f (t ) = h (t )e = e C, B > 0
B B
R(t) = e
Which corresponds to a 2-parameter Weibull distribution
∫ ( )
( )= =
( )= (1-29-1)
( )= (1-29-2)
Chap 1 Introduction and Basic Concepts 56
and
1 1-p
E(X)= 1-31-2 Var(X)= 1-31-3
p p2
( )= (1 − ) , = 0,1, … (1-32-1)
E(X)=np. (1-32-2)
The variance is
Var(X)=np(1-p). (1-32-3)
Example 1-6
Find the mean of a binomial distribution .
59 Reliabilty Engineering
Solution
E(X)=∑ (1 − )
Since the value of the first term( i.e. for x=0) is zero and
k k 1 then
i k
i i 1
n n x n n 1
n n x n 1 x n x x 1 n x
E ( X ) x p x 1 p n p 1 p np p 1 p
x 0 x x 1 x 1 x 1 x 1
n 1 n 1
y n 1 y
x 1 y E (X ) np p 1 p
y o y
n
n
(a b ) n a i b n i E ( X ) np ( p 1 p ) n 1 np
i 1 i
End of Example
and
Var(X)=E(X)= (1-33-2)
( )
( )= = 0,1,2, … (1-33-3)
!
and
Var(X)=E(X)= (1-33-4)
Solution
=0.0256
61 Reliabilty Engineering
(0.0256×515.8)0e-0.0256×515.8
Pr(X=0)= =1.8×10-6
0!
Distribution h(x )
( = ) ( ≤ )
Discrete uniform 1 x 1 x
x ∈ {0, 1 , … , n}
n n n x
Binomial
( − )
(1 − p) ( − )
∑ ( − )
∈ { , ,…, }
Poisson
x k x
x ∈ {0, 1 , … } e x x 1 k
e
x! k 1 k ! x !(1 )
k 0 k!
Geometric
( − ) −( − ) p
x ∈ {0, 1 , … }
Hyper-geometric
−
max 0, m+ n -N −
− −
− −
x ∑ ( , ) −
−
min n, m
Chap 1 Introduction and Basic Concepts 62
Grouped Observations
.
Chap 1 Introduction and Basic Concepts 64
.
interval (a1 b1 ) (a2 b2 ) (ai bi ) (an bn ) sum
... ..
.
Frequ. o
.. ..
and let
N (bi )
Rˆ (bi ) (1-34)
N0
65 Reliabilty Engineering
ℎ( ) = (1-35-2)
N (a i ) ( t ) i
1
N (bi )
( )=1− ( )= (1-37)
N0
Example 1-8
Solution
b ( ) ℎ( )
i ( ) multiply multiply
N(a i ) N(bi )
ai N(b ) by: by:
i i
ai ( t ) i (b )
10 10
=
N ai t b i
b1 =
( 1)
→ (20000) =
(20000)
=
46 − 19 27
46
=
46
= 0.587
N0 N0
N ( a i ) − N ( bi )
ℎ( ) = →
N ( ai ) bi ai
hˆ ( t ) = 46 − 27
= 0.207 × 10 0 t 20000
(46)(20000)
ai − bi
( )= →
N0 × (b a )
i i
20000<t 40000 ( )=
( )( )
= 0.12 × 10
Example 1-9
Solution
N (b i )
( )= , b1 = , ℎ( ) = .
N 0 × (b a ) N0 N ( a i ) (b a )
i i i i
fˆ t = N 0 (b1i ai ) Rˆ (b i )
hˆ t = N(a )1(b
i ai )
a i t bi for ai t bi
i
N (bi )
(× 0.01) N0
b
(× 0.01)
9
1 0 5 1/(10*5)= 2 1/(10*5)= 2
10
8
2 5+ 10 1/(10*5)=2 1/(9*5)= 2.22
10
7
3 10+ 17.5 1/(10*7.5)= 1.33 1/(8*7.5)= 1.67
10
6
4 17.5+ 30 1/(10*12.5)=0. 1/(7*12.5)= 1.14
10
5 30+ 40 1/(10*10)= 1 5
1/(6*10)= 1.67
10
4
6 40+ 55 1/(10*15)=0. 67 1/(5*15)= 1.33
10
3
7 55+ 67.5 1/(10*12.5)=0. 8 1/(4*12.5)= 2
10
2
8 67.5+ 82.5 1/(10*15)=0. 67 1/(3*15)= 2.22
10
1
9 82.5+ 100 1/(10*17.5)=0. 1/(2*17.5)= 2.86
10
10 100+ 117.5 1/(10*17.5)=0. 0 1/(1*17.5)= 5.7
Example 1-10
800 units of a product were placed on the life test and every 3
hours the number of failures were recorded (see table below).
Estimate and plot the density function, the reliability function &
hazard function for the time intervals.(Example. 1.2 Grosh,1989 )
Solution
f
fˆ t =
f
hˆ t =
i , Rˆ (bi ) NN(b0i ) , i
N (b a ) N (a )(b a )
0 i i i i i
Chap 1 Introduction and Basic Concepts 70
a i t bi a i t bi
i ai bi
failures
1 0 3 185 185/(800×3)=0.0771 615 185/(800×3)=0.0771
800
2 3 6 42 42/(800×3)=0.0175 573 42/(615×3)=0.0227
800
3 6 9 36 36/(800×3)=0.015 537 36/(573×3)=0.0209
800
4 9 12 30 30/(800×3)=0.0125 507 30/(537×3)=0.0175
800
490
5 12 15 17 17/(800×3)=0.0071 17/(507×3)=0.0112
800
482
6 15 18 8 8/(800×3)=0.0033 8/(490×3)=0.0054
800
468
7 18 21 14 14/(800×3)=0.0058 14/(482×3)=0.0097
800
8 21 24 9 9/(800×3)=0.00375 459 9/(468×3)=0.0064
800
9 24 27 6 6/(800×3)=0.0025 453 6/(459×3)=0.0044
800
10 27 30 3 3/(800×3)=0.0013 450 3/(453×3)=0.0022
800
71 Reliabilty Engineering
Cumulative Distr .
function and
Reliability function
Example 1-11
Solution
h(t)٬ f(t) , R(t) were calculated using Eqs. 1-35-1, 1-36-1 and
1-37. The following table shows the results:
20 18 16 14 13 12 11 10 9 8 7 6
= 2 × 10
D D D D D D D D D D D D
f̂(t)
t in hour
f̂(t) 20 18 16 14 13 12 11 10 9 8 7 6
t in 1000
200 200 200 200 200 200 200 200 200 200 200 200
hours
ˆ ˆ
R(t)=1-F(t) 0.8
0.90 0.73 0.66 0.6 .54 .48 0.43 0.39 .36 .48 0.28
1
End of Example
Chap 1 Introduction and Basic Concepts 74
( ) < () <⋯< ( ).
.
() = (1-38)
.
.
() =1− (1-39)
.
() − ( )
ℎ () = ⟹
( )− () ()
ℎ () = ( . )
(1-40)
( ) ()
( )− ( )
() = ⟹
( )− ()
() = (
(1-41)
( ) () . )
Example 1-12
Solution
.
h = =0.0024 R =1− = R (195)=0.9167
( )( . ) .
.
h = =0.0050 R = 1− = R(300)=0.5595
( )( . ) .
.
h = =0.0082 R = 1− = R(325)=0.3214
( )( . ) .
Chap 1 Introduction and Basic Concepts 76
Example 1-13
Use Eqs. 1-38 through 1-41 to estimate F (t), R (t), f (t) and h (t)
Solution
1 − 0.3
1 1
10 + 0.4
1 5 (10 − 5)(10 + 0.4) 0.9327 (10 − 5)(10 − 1 + 0.7)
=
= 0.0192 = 0.0206
0.0673
2 10 0.0128 0.1635 0.8365 0.0153
3 17. 5 0.0077 0.2596 0.7404 0.0104
4 30 0.0096 0.3558 0.6442 0.0149
5 40 0.0064 0.4519 0.5481 0.0117
6 55 0.0077 0.5481 0.4519 0.0170
7 67.5 0.0064 0.6442 0.3558 0.0180
8 82. 5 0.0055 0.7404 0.2596 0.0212
9 100 0.0055 0.8365 0.1635 0.0336
117.5 Cannot be computed 0.9327 0.0673 Cannot be computed
10
End of Example
1
Prepared by:Mr M Morrdi former student of Kerman University
77 Reliabilty Engineering
Pr X ( ) > = [1 − F (y)] or
( )
( ) = 1 − [1 − ( )] −∞ < < ∞ (1-42)
( )
( )= ( )
( ) −∞< <∞ (1-43)
Example 1-14
( )= ≥ 0, ( )= 1− ≥0
Solution
1− ≥0
( )
( )= 1− 1−1+ =
0 <0
therefore:
( )=1− 1− (1-44)
Proof:
( )= ( ≤ )= ( ) ≤ ⇒
( )= ( ) ≤ = ( ) ≤ ⇒
1
This assumption in the fracture of structures is logical because the number
of their defects are large.
Chap 1 Introduction and Basic Concepts 80
( )= ( )
0≤ ≤
F ( )
(y) = 1 − [1 − F (y)]
( )
=1− 1−
( )= 1− 1− 0≤ ≤
lim ( )= 1− ≥0
→
=1 ( ) , =2 ( ) ,… , = ( )
81 Reliabilty Engineering
( ) = lim ( )=1−
→
( )
( )= = ( )= ≥ 0.End of proof
= = = ( )
Derivation of ( )
i.e. the CDF of sample minimum
or ( ) when →∞
( ) = ⋯ in terms of
F ( )
(y) = Pr(limiting ( ) < )
Pr(U< ⋯)= .
Step 5 Calculate F ( )
(y), considering step 4 and
Eq. 1-45 i. e. Pr ( ≤ ) = 1 − .
Example 1-15
1
( )= 0≤ ≤
0 .
Step 1
−0 ( )
( )= ⇒ ( ) =
Step 2
( )
= ( ) ⇒ = ⇒ ( ) =
Step 3
( ) → = ( )
Step 4
( )
( )= ( ) ≤ = ≤ = ≤
⟹ ( )
( )= ≤ =
Chap 1 Introduction and Basic Concepts 84
Step 5
( )
( )= =1− , ≥0 ⟹
( )
( )=1− , ≥0
( )= ≥ 0 .End of Example
( )
0 ℎ
Example 1-16
( )= ≥0
0 .
85 Reliabilty Engineering
Step 1
( )=1− =1− ( )
( )
Step 2
= ⇒ = 1− ( ) ⟹
( )
1 1 1
( ) = ln = − ln(1 − )
1−
Let =− ; then
1 1 1
( ) =− − − − −⋯
2 3
Step 3
( ) → = × .
Chap 1 Introduction and Basic Concepts 86
Step 4
F ( )
(y) = Pr( ( ) < ) = ≤ = ( ≤ )
⟹ ( )
( )= ( )
Step 5
( )
( )= ( )=1− , ≥0 ⟹
( )
( )=1− , ≥0
nλy
(y) = e y≥0
f ( )
0 others
End of Example
Example 1-17
An applications of the extreme value distribution is to the study the failures of car exhauhs caused by
corrosion.
Solution
Symbols
( )= = 0≤ ≤
(0≤ ≤ ) 1−
( )= ( ≤ )= ( − )≤ = ( − ≤ ) ⟹
( )= − ≤ .
( )= − ≤ ≤ = ( )− −
where
F (x) = , 0 ≤ x ≤ D, i = 1,2, . . . , N
Therefore
( )
( )= ( )− − = - =
.
= ⇒
−1
( )= 0≤ ≤ , i = 1,2, . . . , N
−1
N
= m in ( ) and its CDF is:
i 1
( ) = 1− ( > )… ( > )= 1− 1− ( ) … 1− ( )
( )= ( < ) = 1− 1− ( )
( )=1− ⇒ ( )=1− ( )=
Example 1-18
Solution
∫
Solving the equation = in MATLAB:
landa=127.64972
= 10 , = 10 , = , = 128 , ( ) = 0.9
1
Extreme value 1=Fisher Tippet 1
93 Reliabilty Engineering
1 r 1 r
2r ln i 1t i i 1 ln t i
r r
B
(1-46)
r 1
1
6r
χ ,
is read from Table E or calculated in MATLAB from
chi2cdf(1 − , r − 1)
χ ,
is read from Table E or calculated in MATLAB
Example 1-19
50.1 20.9 31.1 96.5 36.3 99.1 42.6 84.9 6.2 32.0
30.4 87.7 14.2 4.6 2.5 1.8 11.5 84.6 88.6 10.7
Solution
1 r 1 r
2r ln i 1t i i 1 ln t i
r r
B
r 1
1
6r
20
i 1 i
t 50.1 20.9 ... 88.6 10.7 836.3
20
i 1
Ln (t i ) Ln 50.1 Ln 20.9 ... Ln 88.6 Ln10.7 63.9385
836.3 63.94
2×20 Ln -
20 20
B 19.34
20+1
1+
6 20
χ ,
= chi2inv(0.05,19) = 10.1170
χ ,
= chi2inv(0.95,19) = 30.1415
fall outside χ ,
χ ,
. Therefore the
∑ . .
= = ≅ 41.82, ( )= × = .
.
>>Data=[…
50.1
…
10.7];
>> H=kstest(Data, [Data expcdf(Data,mean(Data))] , 0.1)
H=0
i
F (i ) , i 1,...,n
n 1 (1-47)
(Besterfield,1990 page52):
i 0.5
F (i ) (1-47-1)
n
or (Goda,2000 page 287):
i 0.375
F (i ) . (1-47-2)
n 0.25
Chap 1 Introduction and Basic Concepts 98
(1 4 8 )
i a
F (i )
n b
0.27 0.23
where a 0.20 b 0.20
C C
i 0.47
F (i ) (1-48-1)
n 0.43
and fit a line to the points. The closer this line to the
angle bisector of the first quarter of the plane, the better
fits the distribution to the dataset. It is worth knowing
that the better the distribution fits the data set the closer
the correlation coefficient of x(i) & xˆ ( i ) to 1; but the vice
99 Reliabilty Engineering
n x ( i ) xˆ ( i ) x ( i ) xˆ ( i )
R (1-49)
2 2 2
n x (i ) x (i ) n xˆ ( i ) 2 xˆ ( i )
Example 1-20
Solution
To answer, a Q-Q plot is drawn. The mean and variance of
the distribution is estimated as follows:
s
ˆ X 54.81 , ˆ 11.7287
c4
or ( ) = (0.0294,54.81,11.8751) = 32.3698.
QQplot
80
pp=(i-0.375)/(n+0.25)
70
60
Predicted
50
40
30
30 40 50 60 70 80
Observed
F(i )=
i-0.375 .
with F(i )= i-0.5 .
n+0.25 n
Since the points are near to the fitted line and the line is close to
the angle bisector of the first quarter of the X-Y coordinate
plane, it is concluded that the normal distribution fits the dataset.
r=corrcoef(X,Xhat);r=R(1,2)
Chap 1 Introduction and Basic Concepts 102
where
X is the vector containing x(i), i,=1,2,3..
1-15 Convolution
( )= ∗ ( )=∫ F (a − y) ( ) (1-50)
103 Reliabilty Engineering
∫ Pr( + ≤ | = ) ( ) =∫ F (a − y) ( ) .
where
( )= F (a − y) ( ) = F (a − y) ( ) ⇒
( )= (a − y) ( ) .
Notice that
-If the distribution of X+Y is known for us, there is no need for
the above integrations
(a) = (a − y) ( )
(a − y) ≠ 0 0≤ − ≤ 1 or a − 1 ≤ ≤a
Therefore
(a) = (a − y) ( )
∞
⎧
1×1 0≤a≤1
⎪
⎪
= ⇒
⎨
⎪
⎪ 1×1 1≤a≤2
⎩
Chap 1 Introduction and Basic Concepts 106
0≤a≤1
(a) =
2− 1≤a≤2
Example 1-21
Solution No.1
(a) = (a − y) ( )
(a − y) ≠ 0 for a − y ≥ 0 ≤ .
( )
f (a) = ∫ f (a − y)f (y)dy = ∫ e e dy = ae , a ≥ 0
Solution No.2
( )=∫ (a − x) ( ) .
107 Reliabilty Engineering
( )
( )= = a≥0
Solution No.3
Example 1-22
Solution
( )= (a − y) ( )
(a − y) ≠ 0 for −∞ ≤ a − y ≤ ∞ −∞≤y≤∞.
(a) = (a − y) × 1
1
(a) = (t)(− ) = →
√2
and also
F ∗ F= F ∗ (1-54)
Example 1-23
Find the probability density function of the sum of the
exponentially distributed lifetimes of 3 independent
components with parameter .
Solution
Using convolution:
( ) = [ ( )]( )∗ = ( ) ∗ [ ( )]( )∗
= ∗ ( )∗
[f(t)]( )∗ = [ f (t )] ∗
[f (a)]2* f (a t )f (t )dt a 0
a a
[f (a)]2* e (at )e (t )dt 2e a dt 2ae a
0 0
2 t
Continue with convolving te and λe to reach the
− ≥ ≤ , >1 (1-56-1)
or
− ≤ ≤ + >1− , > 1 (1-56-2)
where
and are the mean an standard deviation of Y.
c) A table containing some intervals and frequencies
In this case if the frequency distribution shows the
distribution of Y has only one mode and the mode is the same as
the arithmetic mean and the frequencies decline continuously on
Chap 1 Introduction and Basic Concepts 112
Or
Pr( − < < + )>1− , ≥1 (1-57-2)
Example 1-24
Solution
Normal distribution
a b
Pr a X b Pr( Z )
34-40 46-40
Pr( <Z< )=Pr -3<Z<3
2 2
From Table C:
Pr -3<Z<3 =0.99865 - 0.00135=0.9973
113 Reliabilty Engineering
With MATLAB:
More than 1- 1 4
=1- 2 of the product fall within (34 46) ;
2.25k 2 9k
or
1 1
Pr 34<X<46 >1- =1- =%95.06
2 2.25×9
2.25k
L f ( x1 xn ) f x1 ( x1 ) f xn ( xn )
(1-58-1)
1
Refer to page 25.33 Handbook of Reliability by Irenson et al,1996
115 Reliabilty Engineering
L ( e x 1 )... ( e x n ) n e x i
(1-58-3)
for binomial distribution:
=∏ ∑
× (1 − )∑ ( )
(1-58-4)
1
A function is monotonic if its first derivative is always positive or negative
Chap 1 Introduction and Basic Concepts 116
Example 1-25
a)Given a random sample ,…, from an exponential
distribution, use MLE method to estimate .
b)If the sample is (1.1, 0.9, 1.21, 0.8)calculate the value of .
Solution
a)Since the sample is random, 's are independent.
L n e x i lnL ln n e x i nln x i
d lnL n 1
0 MLE ( ) ˆ
d xi X
b) ( )= = =∑ = =1
. . . .
Lemma1
1
Based on page 290 Bowker & Liberman (1972).
117 Reliabilty Engineering
End of lemma
Example 1-26
If ,…, is a random sample of size n, taken from an
Solution
Since θ (λ) = <0 therefore θ has a unique inverse and
Example 1-27
Given a random sample …., from a Weibull
distribution with location parameter A=0, use MLE method to
derive the relations for calculating the scale and shape
parameters B and C.
Answer
∑ ∑
= ∑
− (1-59-1)
∑
= (1-59-2)
Chap 1 Introduction and Basic Concepts 118
Example 1-29
Solution
k
n
L Px i Px i p x i q n x i
i 1 xi
k
n n n x k n x k
L p x i q n x i p x 1q n x 1 ... p q
i 1 x i x1 xk
k k
k
n xi
i 1
( n x i )
i 1
L p q
i 1 x i
k k
k
n xi (n x i )
ln(L ) ln( ) ln( p i 1 ) ln(q i 1 )
i 1 x i
k
n k k
ln(L ) ln x i ln p (n x i ) ln(1 p ),
i 1 x i i 1 i 1
k k k k
x i (n x i ) 1 p (n x i ) nk x i
0 0
ln(L ) i 1 i 1 i 1 i 1
1 p
k
k
p p p
x
i 1
i x
i 1
i
1 xi
1 k 1 MLE ( p )
nk i 1
p kn
xi i 1
End of Example
k 1 MLE ( p )
x
n
Example 1-30
The following sample shows the lifetime(in year) of some
units randomly taken from a batch of a device having an
exponentially distributed lifetime. Estimate the parameter of the
distribution.
0.04 0.15 0.04 0.09 0.03 0.01 0.04 0.06 0.01 0.15
Solution
X=[ 0.04 0.15 0.04 0.09 0.03 0.01 0.04 0.06 0.01 0.15 ] ; expfit(X)
failures.
Example 1-30
The annual number of failures of a device has a Poisson
distribution; given the following sample estimate the distribution
parameter.
17 13 19 8 17 17 12 19 18 19
Solution
Example 1-30
2- Methods of Moments(MOM)
Steps of MOM
To estimate parameters ,…, of the statistical
distribution of random variable X ,
i-Compute ( ), = , … , in terms of ,…, .
Notice that always the first moment of a distribution is its
mean and the second moment is equal to the variance of the
distribution plus the squared mean.
ii-Form the k equations
( )= , j=1, ….,k
iii-Solve the equations for the parameters ,…, .
The resulting values are called method of moments
estimators for the parameters.
123 Reliabilty Engineering
Example 1-32
Let ,…, be a random sample taken from an
exponential distribution with parameter , Estimate the
parameter by the method of moments.
Solution
1 X i n
E (X ) M 1 ˆ
n X i
End of example
Example 1-33
Let ,…, be a random sample taken from an normal
distribution with parameters , estimate the
parameters by the method of moments.
Solution
Xi
E X M 1
n
2 2
E
X M 2 2 2 X i
n
X i
ˆ X
2
n
2 ^
2
Xi
2
2 X i 2 X
n
n
Chap 1 Introduction and Basic Concepts 124
(x i X )2 x i
2
nX 2
S2 i 1 i 1
is a good estimate for the
n 1 n 1
B-Parameter estimation
thetahat=expfit(Data)
[muhat,sigmahat] = normfit(Data)
BhatandChat= wblfit(data)
H = kstest(Data,CDF, )
where
if omitted or
Example 1-33
[110, 520, 645, 680, 330, 75, 95, 480, 360, 575, 1065, 170,
415, 15, 20,1275, 270, 90, 1500, 1923, 715, 1523,427, 730,
1120, 390, 240, 40, 220, 673, 2397, 1032, 315]
Solution
Entering the data as a column vector:
>>Data=[…
110
520
…
1032
315];
Giving the command:
Chap 1 Introduction and Basic Concepts 128
The answer for H is 0; i.e. it is not rejected that the data belongs
to an exponential distribution with significance level = 0.05 .
>> X=[….data];
pd=makedist(distribution name ); e.g.
pd=makedist ('exponential')
pd=makedist ('Gamma')
qqplot(X,pd)
The more the points near to the line and the line near to the
bisector of the first quarter, the better the distribution fits the
data.
D Calculation of Reliability
answer :
>>x=expinv(0.3935,100)
x=50
answer
>> p=1-expcdf(50,100)
p = 0.6065
c)Find the lifetime value for which the reliability of the device
is 0.6065.
>>x= expinv(1-0.6065,100)
x = 50.0051
131 Reliabilty Engineering
∫
which satisfies the equation = , the following
( )
landa=127.65
Exercises
(b) For a 15 hour guaranteed life what should the choice be?
h(t) = b − c t, 0≤t≤t
h(t) = b − c t − c (t − t ), t < ≤t
- 2 t
i)h(t)=10-6t ii)h(t)=10-4 e 100 t km
that h '( x ) h 2 ( x ) .
( ) ( ) ( ) ( ) ( )
Solution:ℎ( ) = ( )
⟹ℎ( )= ( )
( ) ( ) ( ) ( ) × ( ) ( ) ( ) ( ) ( )
⟹ℎ( )= = ⟹ℎ( )= =
( ) ( ) ( )
h '( x ) h 2 ( x )
( ) ( )
×
( ) ( )
135 Reliabilty Engineering
Chapter 2
Static Models
in
Reliability
Chap. 2 Static Reliability Models 136
2
Static Reliability Models
Here, the word static means that the time coordinate is not
presented in the calculations. In modeling a system from a
reliability stand point using static models, the component or
subsystem reliabilities are considered to be constants; thus some
base time period is implied(K&L page 55). Before dealing with
some conventional component configurations, a graph is
introduced below.
137 Reliabilty Engineering
1 2 3 n
= ( ∩ …∩ ).
= ( )… ( )= ×…× ⇒
=∏ . (2-1)
n
Since 0 < <1 therefore ≤ min { } ; that is the reliability of
i1
Therefore
( − 1)
(1 − ) = 1 + (− ) + (− ) +. . . +(− )
2
: = (1 − ) ≅ 1 − , < 1.
n n
R sys 1 q i q i 1
i 1 i 1
Example 2-1
Solution
=1− = 0.9949
Using approximation:
=1− = .995
Chap. 2 Static Reliability Models 140
= ( ∩ ∩ …∩ ∩ ),
Using the chain rule for factorization, this joint probability can
be rewritten as follows:
=Pr(nth component is on|other components are on)× … × Pr(1st component on) (2-2)
= 1 − (1 − )(1 − ) … (1 − )=1− (1 − ) ( − )
141 Reliabilty Engineering
= 1 − (1 − )(1 − ) … (1 − )=1− (1 − ). (2 = 3 − 2)
End of proof
Example 2-2
Solution
1
a- = 1 − (1 − )6 ≥ 0.80 ⇒ ≥ 1 − (1 − 0.8)6 =
02353 = 23.53%
b- = 1 − (0.5) = 0.999 ⇒ = 10
Solution
1-Active redundant
2-Standby redundant
=
Pr(system works)
works = Pr( A works) + Pr ( A does not work ) × Pr( B works)
works
= + ((1 − )= + −
Therefore the reliability of a two-component
component standby system is
the same as to that of a 2-component
component active system as given below:
2
R 1 (1 R 1 )(1 R 2 ) Ps=1 (2-3-1
1)
sys
2-4-1-2 n-component
component system with 1 active and n-1
1 standby-
standby
Perfect switching
Consider an n-component
component standby system with one normally
operating subsystem and n-1 in standby status(Fig 2-3--1). The
system is functional until n failures occur.. The reliability of this
system is:
n
standby R 1 (1 R 1 ).....(1 R n ) (2--3-2)
sys
Fig. 2-3-1
2 RDB of n-component
component standby system
(https://egyankosh.ac.in/bitstream/123456789/35169/1/Unit 15.pdf)
(https://egyankosh.ac.in/bitstream/123456789/35169/1/Unit-15.pdf)
Chap. 2 Static Reliability Models 148
Stand by 3
Reliabilty1-active 1 standby
R 1 1 [1 (1 R1 )(1 R 2 )]1 R 3
sys
Standby 3
R (1 R1 )(1 R 2 )(1 R 3 ) Ps 1
sys
Example 2-5
Solution
2
R =1-(1-R1 )(1-R 2 ) =1-(1-0.9)(1-0.8)=0.98
sys
149 Reliabilty Engineering
Example 2-6
2 4
1 2
B
3 4
Solution
= 1 − (1 − )(1 − )
Example 2-8
In example 2-7 let R 1 0.95 , R 2 0.85 , R3 0.75, R 4 0.8 and
R B = 1 -(1 -0 .9 5 × 0 .8 5 )(1 -0 .7 5 × 0 .8 )= 0 .9 2 3 0
a b c
a b c
= 1 − [(1 − )]
1 − (1 − )(1 − ) = 2 −
= (2 − ) 2 − (2 − )
If = = then
− =6 (1 − ) > 0 ⇒ > .
Example 2-9
1 2 3 n
1 2 3 n
1 2 3 n
1 1 1
2 2 2
3 3 3
n n n
Fig.2-5-2 The RBD of part b of Example 2-9
Solution
The figure shows that the system has m parallel subsystems each
having n components. The reliability of this LL redundancy
system is given by:
= [1 − ∏ (1 − )] (2-4-2 )
1 2
4
3
Chap. 2 Static Reliability Models 156
Example 2-11
In the following 2-component system, the reliability of the
components are respectively R , R . If the reliability of this
system is not enough for us, and we have the following two
options to use, which one has more reliability? A or B?
Solution
= ( + − )
=[ + −( )( )]( )
− =( )( − )
( = )= (1 − )
The n-k system will function if there are no more than n-k
failures. Thus the reliability is as follows:
( ≤ − )= (1 − )
ℎ :
Rsys = binocdf(n − k, n, 1 − R) (2-5-1)
The system failure probabilty or system unreliabilty is calculated
from:
1− =1− ( − , ,1 − )
Alternatively,
we could say that system work as far as n-k components out
of n component fail. Then the probability of system
failure is given by:
( > − )= (1 − )
−
=1− (1 − )
= − +1
(1 − ) −
=
=
= ∫ (1 − ) .
In summary, if the reliability of components in a k-n system is
denoted by R then:
=∑ (1 − ) = ∫ (1 − ) . (2-6)
where
R component reliability
system reliability
n Total number of components in the system
i no. of components that work
n-i no. of components that fail
using MATLAB:
Rsys = 1 − binocdf(k − 1, n, R) (2-6-1)
Example 2-12
A system has a 3 out of 5 active redundancy configuration.
The reliability of each component is R=0.9. Calculate the
reliability of the system.
Solution
Using integral of Eq. 2-6:
5
=3 3 ∫0 3−1 (
1 − )5−3 =6R5-15R4+10R3
R=0.9 ⟹ =0.99144
k−1 3−1
n 5
Rsys = 1 − Ri (1 − R)n−i = 1 − 0.9i (1 − 0.9)5−i
i i
i=0 i=0
With Matlab:
Eq. 6-2
3*nchoosek(5,3)*int(R^2*(1-R)^2)= 6R5-15R4+10R3
Eq. 2-5-1
( − , , 1 − ), = (5 − 3,5,1 − .9) =
0.9914
Eq. 2-6-1:
=1− (2,5,0.9) = 1 − 0.00856=0.99144
End of Example
Example 2-13
The system shown in the following figure has only 4
components A , A , A and A . Each component works,
independently of all the other components .Determine the
configuration of the system.
Chap. 2 Static Reliability Models 162
Solution
It is evident from the figure that the system works in the
following conditions:
If A , A & work,
If A , A & A work,
If A , A & work,
If A , A & work,
If A ,A , A & work.
Therefore the system works if at least 3 components work; i.e.
the system has a k=3/n=4 configuration. End of Example
Example 2-14
The system shown in the following figure has only 3
components 1,2 3 . Each component works, independently
of the two others .Determine the configuration of the system.
1 2
2 3
1 3
Solution
It is evident from the figure that the system works in the
following conditions:
If 1 &2 work,
If 1 &3 work,
If 2 &3 work,
163 Reliabilty Engineering
= (1 − )
substituting k=n:
= (1 − ) = (1 − ) =
k=n⟹ =
substituting k=1:
= (1 − )
= (1 − ) − (1 − ) ⟹
0
= (1 − ) − (1 − )
( + ) = ( )
then
(1 − ) = (1 − + ) =1 =1
therefore
= (1 − ) − (1 − ) = 1 − (1 − )
n
R = (1 − R) (R)
x
Then
=1− (1 − ) (R)
(1 − ) (R) ( )
=1−
− +1
(1 − ) (R) ( )
− −⋯
− +2
Then
≤1− (1 − ) (2 − 7)
− +1
Example 2-16
A pressure vessel is equipped with 6 relief values that work
independent of each other. Three values are enough for the safe
operation of the vessel. The failure probability of each valve is
0.1%. Calculate the probability of safe working of the vessel.
Solution
The reliability of this system, which has a 3-out-of-6
reliability configuration, is given by Eq. 2-6-1:
= 1 − 0.001 = 0.999
= 1− ( − 1, , )
1− (2,6,0.999) = 0.99999999998502
Notice that:
In Reliability literature related to k-out-of-n configuration
sometimes the binary variables and are defined as:
1
=
0 ℎ
⎧1 ≥
⎪
=
⎨
⎪0 < .
⎩
Solution
Component 2a is chosen as the key element and system
operation is conditioned on:
i) 2a works all the time
ii) 2a does not work at all
With the following symbols :
Y The event that the original system works successfully
x The event that Component 2a fails
169 Reliabilty Engineering
{ }= ( | ) ( )+ ( | ) ( ) (2-8)
Let { }= { }=
Then { }= 1− and and Eq. 2-8 is could be written as:
= ( )= (1 − )+ (2-9)
= ( | )=
Chap. 2 Static Reliability Models 170
= ( | ) = [1 − (1 − ) ][1 − (1 − ) ]
R R 1 (1 r r )(1 r ) B r r r r r r
2 3 23 1 1 23 123
result in:
ABC
R1 (2-10-1)
B
[1 (1 r3 )(1 r2r1 )][1 (1 r1 )(1 r2r3 )][1 (1 r2 )(1 r3 r1 )]
1 R1
1 (1 r1 )(1 r2r3 )
ABC
2 R 2 (2-10-2)
C
ABC
3 R 3 (2-10-3)
A
A 1 (1 r r )(1 r ) B 1 (1 r r )(1 r )
1 2 3 2 3 1
C 1 (1 r3 r1 )(1 r2 )
Example 2-18
Solution
0.900497
ABC=0.900497 R1 = =0.9831
0.916
0.900497 0.900497
R2 = = 0.9725 R 3 = =.9419
0.926 0.956
%% calculations
options=optimset('Display','Off');
Eq=fsolve(@Delta,[0.5,0.5,0.5],options);
Eq1=Eq(1);
Eq2=Eq(2);
Eq3=Eq(3);
display(['r1= ' num2str(Eq1), ' r2= '
num2str(Eq2) ,' r3= ' num2str(Eq3)])
The sub-code Delta.m used above is as follows:
function W=Delta(r)
global A B C;
W=[1-(1-r(1)*r(2))*(1-r(3))-A;
1-(1-r(3)*r(2))*(1-r(1))-B;
1-(1-r(1)*r(3))*(1-r(2))-C];
Performing star2Delta using the data in Example 2-11
>>star2Delta
Please Insert R1 value: 0.9831
Please Insert R2 value: 0.9725
Please Insert R3 value: 0.9419
Results: r1= 0.69991 r2= 0.79993 r3= 0.90017
Special case: identical components
If a delta configuration consists of 3 identical component
with a reliability of r ; the equivalent star configuration must
have the following components:
Example 2-19
Solution
End of Example
Example 2-19
where
A = B = C = l-[l-(0.9)(0.9)](l-0.9)=0.981
l
3
R =R =R = 0.981 =0.9904
1 2 3 0.981
Suppose Element No. 3 in the original RBD works all the time.
Replacing it with a line would result the following RBD with a
R PS =[1-(1-0.9)(1-0.9)][1-(1-0.9)(1-0.9)]=0.9801
=[1-(1-0.9)(1-0.9)][1-(1-0.9)(1-0.9)]=0
2-8
8 Calculation of upper and lower bounds for
complex system using cut and tie set
set analysis
Notice that
-Each
Each complex system has usually several minimal path and cut
sets. A set could be both cut and path.
Example 2-21
2 (Grosh,1989 page 125)
Find the minimal path and cut sets for the following RBD.
Solution
2 C4C5
3 C1C3C5
4 C2C3C4
End of Example
xample
Example 2-22
2 (O'connor,2003 page175)
Solution
End of example
Example 2-23
2 (Shooman,2002 page285)
What are the minimal cut and path sets between source a and
target b given below?
Solution
If all components fail but Component 1,the connection
between a and b does not interrupts then the set {1} is a tie set
If all components fail but Components 2& 5, the connection
between a and b does not interrupts then the set {2,5} is a tie
set
…
Functioning of the set with minimum elements 6,4 ensures
ensure the
system operation , then Set {6,4} is a minimal path.
…..
183 Reliabilty Engineering
The other cut and tie sets are determined with similar reasoning.
The following table shows the minimal cut sets and minimal
paths.
End of Example
the reliability of the original system (based on Grosh, 1989 p 125 -6)
RN 2 Rsys RN 1 (2-13)
For some reasons such as dependence of the subsystems of
Network N1, The calculated upper bound( R N 1 ) is usually over-
over
Example 2-24
2
Draw the auxiliary networks for calculating the upper and
lower bounds for RBD of Example 2-21.
2
Solution a)Auxiliary network N1
Network N1 which is a parallel configuration of all the
minimal path elements in series is shown below:
Network N2
The reliability of the above RBD given by the following
relationship is the lower bound for the original system
reliability:
RN2 =
1-(1-r1)(1-r2 ) 1-(1-r4 )(1-r5 ) 1-(1-r1)(1-r3)(1-r5 ) 1-(1-r2 )(1-r3)(1-r4 )
Then RN 2 Rsys RN 1
Example 2-25
2
Draw the auxiliary networks for calculating the upper and
lower bounds for the following RBD. The number in the box is
the reliability of the component.
Chap. 2 Static Reliability Models 186
Solution
The minimal paths are:
1-3-6-7 1-3-5-7 1-2-4-7 1-2-5-7
The minimal cut set are:
1 7 2-3 4-5-6 2-5-6 3-4-5
Network N1
187 Reliabilty Engineering
Network N2
C
1 (1 R i ) =
j 1 i A j
Ri
K 1 i B K
Solution
The minimal cut sets and the minimal paths are given in the
following table.
189 Reliabilty Engineering
Network N1 Network N2
T=2, B1={1, 2, 4}, B2={ 3, 4} C=3, A1={1, 3}, A2={2, 3}, A3={4}
and
Pr( Hi ) Pr( B / Hi )
Pr( Hi | B) k (2-15)
Pr(H ) Pr(B H )
i 1
i i
B2 = the event that the mean life is 30,000 kilometers, Pr(B2 ) 0.2
SS=B1 UB2
Pr(Bi )Pr(A Bi )
Pr(Bi A)= i=1,2
Pr(B1 )Pr(A B1 )+Pr(B2 )Pr(A B2 )
Pr(B1 )=0.8
40000-50000
Pr(A B1 )=Pr(X³40000,μ=50000)=Pr Z>
(50000)(0.10)
=Pr(Z³-2)=0.9772
40000-30000
Pr(A B2 )=Pr(Z> )=0.00045
30000×0.10
Pr(B1 )Pr(A B1 )
Pr(B1|A)= =0.9999
Pr(B1 )Pr(A B1 )+Pr(B2 )Pr(A B2 )
Pr(B2 |A)=1-Pr(B1|A)=1-0.9999
Suppose later, after production stage, one of these devices
was selected for a mission of 35000-km- experiment. Assuming
a normally distributed life time with mean 50000 km and
=5000 km the reliability is given by
35000-50000
R=Pr(X>35000)=Pr(Z> =Pr(Z>-3)=0.99865
5000
We are 99.99% confident about the mean used for this
calculation because Pr(B1 | A) 0.9999. Therefore we are 99.99%
Exercises1
1
From K&L chapter 3
Chap. 2 Static Reliability Models 194
Ans. RI RII .
11)To have a 6-component series system of at least 95%
reliability, how many components do you suggest?
12) A system consists of several components with 94%
reliability To have a system with 95% reliability, what
configuration do you suggest and many components?
14) In the following RBD, each component has the indicated
reliability. B2 is a standby component which replaces B1 by a
switch. The failure probability of the switch and the standby
component when they are needed is negligible. Calculate the
reliability the system with the given RBD and compare it with
the case as if there is no redundant standby component.
Chap. 2 Static Reliability Models 198
Chapter 3
Reliability
Considerations
in Design
+UGF
Technique
Chap. 3 Reliability in Design + UGF Technique 200
3
Reliability Considerations in Design + UGF Technique
1 2 3 n
1
R=0.9
0.95
R=0.6
0.9 R=0.7
R=component reliability in
System Reliability
0.8
0.75
0.7
0.65
0.6
1 2 3 4 5 6 7 8 9 10
no. of components
.
Reliability of LL redundancy
= [1 − (1 − ) ] (3-1)
Reliability of HL redundancy
The reliability of HL configuration in Fig 3.9 is:
= 1 − [1 − ] ، (3-2)
where
1
Jinhua, Mi, et al, 2015
Belief Universal Generating Function Analysis of Multi-State Systems
Under Epistemic Uncertainty and Common Cause Failures IEEE
Transactions on Reliability Vol. 64 No.4
Chap. 3 Reliability in Design + UGF Technique 212
∑
( )= ( ).
i 1
(k )
d
3. E (X k ) ( )|
dt k
i.e. the kth derivative of ( )with respect to t gives the
th
value of k moment of the distribution about the origin,
for t=0.
4. The additive property of some distributions such as
Poisson and normal distributions could be verified by
MGF.
3.2-2 Z-transform or probability generating function of
discrete random variables
The probability generating function of a discrete random
variable is defined as follows:
( )= ( )= (3 − 3)
End of Example
( )= ( )= ( )= (3 − 5)
=1 ⟹ ( )=∑ = ( )
3.The z-transform of the sum of some independent rv's is equal
to the product of the z-transforms of the rv's:
k
∑
( )= ( ) (3 − 6)
i 1
4. If in the definition of the MGF e is replaced by z, then one
gets z-transform of the random variable.
Example 3.3
End of Example
215 Reliabilty Engineering
x= 1, … , p= 1, … ,
that the function f takes on some value is equal to the sum of the
combination producing this value(Levitin,2010 page 6). Let
be a set of combinations producing the value . If the total
number of different realization of the function ( … ) is H,
then the probability function of f is:
Y f h :1 h H ,
n (3 − 7)
q pij i :1 h H
( x 1 j1 ,..., x 1 j n )A h i 1`
Example 3-4
Solution
All possible combinations of X 1 , X 2 and the probability
function of Y is given in the following table:
X1
= X2 ( )
1 1 1 0.5
0.6 0.1=0.06
2 4 0.5 2 0.4 0.1=0.04
3 1 1 1 0.6 0.6=0.36
4 4 1 4 0.4 0.6=0.24
5 1 2 1 0.6 0.3=0.18
6 4 2 16 0.4 0.3=0.12
As the table shows some combinations have the same value .
Since all combinations are mutually independent then the
217 Reliabilty Engineering
( )= (3 − 8)
UGF Technique
( ) =⊗ ( ), ( )] = [ ( )⊗ ( ) (3 − 9)
(z) ( ) =⊗ ( ( ), … , ( )) (3 − 10)
U (z) = .
n
f ( x i j1 ,..., x n j n )
⊗ P Z =
pij i z
… (3 − 11)
i 1
The technique based on using z-transform and composition
operators ⊗ is named universal z-transform or universal
(moment) generating function (UGF) technique (Livitin, 2010
page 8 ). UGF technique has applications such as finding the
probability function of an arbitrary function of several
independent random variables and finding the reliability of
complicated systems. For other applications refer to reverences
such as chapter 2& 3 in Levitin(2010).
219 Reliabilty Engineering
U Y ( z ) U 1 ( z ) power U 2 ( z )
U Y ( z ) (0.6 z 1 0.4z 4 ) power (0.1z 0.5 0.6z 1 0.3z 2 )
1 2 4 16
( ) = 0.6 + 0.04 + 0.24 + 0.12
= (1, 2, 4, 16)
probabilities = = (0.6, 0. 04, 0.24, 0.12)
End of example
( ) ( )
x 1,0 =8 p1,0 =0.7 5 0.6 0.1 0 0.6 1 0.5
x 1,1 =10 p1,1 = 0.3 8 0.3 0.5 2 0.4 1.5 0.5
12 0.1 0.4
= ( … )=[ ( , )+ ( , )]( ).
Solution
The total number of term multiplication procedures that one has
to perform using enumerative approach is
2 × 3 × 3 × 2 × 2 = 72; however applying UGF technique as
performed below reduces this amount to only 26 (Livitin,2010
page11).
= ( , ) = ( , ) = +
( )= ( )⊗ ( )=
= (0.7 + 0.3 ) ⊗ (0.6 + 0.3 + 0.1 )=
Chap. 3 Reliability in Design + UGF Technique 222
( , )
0.42 ( , )
+ 0.21 +
. × . . × .
( , ) ( , )
0.07 018
. × . . × .
( , ) ( , )
+0.09 + 0.03 ⟹
( ) = 0.63 + 0.27 + 0.1
( )= ( )⊗ ( )
= (0.1 + 0.5 + 0.4 ) ⊗ (0.6 + 0.4 )=
( , ) ( , ) ( , [) ( , )
0.06 +0.04 +0.3 0.2
( , ) ( , )
+0.24 + 0.16 ⟹
( ) = 0.64 + 0.36
( )= ( )⊗ ( )=
( )= ( ) ⊗× ( )
= (0.4032 + 0.3996 + 0.01612
.
+ 0.036 ) ⊗× (0.5 + 0.5 )
After necessary calculations and simplification, the final answer
for U (z) is:
End of Example
Example 3-7
The universal generating function of a system is:
Symbols
System Reliability
The reliability of jth subsystem
The state of jth subsystem(either 1=working or 0=down)
X The system structure function
The UGF method is very effective for the reliability analysis of
multistate systems; however it could be used for binary-state
systems, though that effective as compared to conventional
methods (see Kuo &Zuo,2003).
= [ , ( , )]
where X is the state of jth subsystem(with 2 values :either
x =1=working or x =0=down).
Let
Rj=Pj be the reliability of jth subsystem for a fixed mission
time, the probability that is on working conditions during the
mission time
and 1-Pj be the probability that the jth subsystem is down.
Then the expected value of is:
=0 1− +1 = =
where R is jth subsystem reliability.
Therefore for a fixed mission time the system reliability
equals the expected value of .
Similarly the reliability of the system for a fixed mission
equals the expected value of the system structure function X:
= ( ), (3-12)
Where
= ( , ,…, )
The state of jth subsystem(either 1=working or
0=failed)
= = 1− = 0 or = 1 = 1,2,3
( = ∩ = ∩ = )=
(1 − ) (1 − ) (1 − )
value of
, , Pr( X x ) p X (x)
= 0 down, = 1 working
0,1,0 0 (1 − ) (1 − )
0,1,1 1 (1 − )
1,0,0 0 (1 − )(1 − )
1,0,1 1 (1 − )
1,1,0 0 (1 − )
1,1,1 1
According to Eq. 3-12 for a binary system :
= ( )= ( ) = [(1 − )(1 − )(1 − ) × 0] + ⋯ + [ × 1]
=( + − ) = [1 − (1 − )(1 − )].
R3*(1-(1-R1)*(1-R2))= 0.8458
( )= ( )= ( )=
Then this series-parallel system reliability is:
( ) = { ( )] = 3( )[ 1( )+ 2( )− 1( ) 2( ))]=
( )
+ − ,
( )= ( )[1 − (1- ( ))(1 − ( ))]= [1 − (1- )(1 − )]
End of Example
( )
E(X)= (1) = = |z=1 (3-14)
where U (z ) [U (z ),...,U n (z )] .
1
Exercises
0.3and 0.5 . Calculate the UGF or U(z) of this system and the 10-daay
reliability of the system.
Keep in mind
that you are
never absent
from
God’s sight,
so keep looking
how you are
acting
229 Reliabilty Engineering
Chapter 4
Structural
reliability
Analysis
Chap. 4 Structural Reliability Analysis 230
4.1 Introduction
Designers of systems such as structures take many factors
into considerations including the reliability. Strength(capacity)
and load(stress) are 2 variables that affect the reliability of
structures(dams ,bridges; communication networks and
antenna). To be reliable, structures require to withstand ultimate
loads without failure.
There are 2 approaches for this purpose: deterministic and
probabilistic.
231 Reliabilty Engineering
SF= (4-1)
1
When the strength and load are independent random
variables, the average SF,E(SF), is approximately:
2S
1 2
E ( )
E SF
2
E ( s ) E (s )
MIL-HDBK-17-3E, Working Draft page 6-7
https://www.gla.ac.uk/external/asranet/Resources/milhdbk.pdf
233 Reliabilty Engineering
R P r( S ) (4-2-1)
R P r( S 0 ) (4-2-2)
(4-2-3)
R Pr( 1)
S
R P r( S F 1) (4-2-4)
Let = − then
= ( > 0) (4-2-5)
R Pr S f ,S , S d dS , (4 3)
S
where f ,S , S is the joint probability density function of
R f ( ) f s (s )ds d f s (s ) f ( )d ds
0 s 0 s 0 s
R f s (s ) 1 F (s )ds (4 4)
s 0
where
f ( y s )f s (s )ds y 0
0
f Y ( y ) f ( y s )f s (s )ds (4-5)
s f ( y s )f (s )ds y 0
s
y
R Pr(Y 0) fY (y ) (4-6)
y 0
f ( y s )f
y 0 s 0
s (s )dsdy
Example 4-1
The stress and the strength distributions for a component are
uniform over the interval :
Strength: [15 25]
Stress : [20 25
235 Reliabilty Engineering
Solution
¥ 25 1 s-15
R= f s (s) 1-Fδ (s) ds= 1- ds
s=0 s=20 25-20 25-15
1 25
R= (25-s)ds=0.25
50 s=20
100(1-R)=75% break.
End of Example
where
are the means of the strength and load,
are the variances of the strength and load.
In a structure, if < , the more s far from the less failure
probability and the more reliability. Then the more the
denominator the more the reliability; the less variation of the
load and strength (or the less the denominator), the more we are
confident. Therefore the greater SM>0, the more reliable the
Chap. 4 Structural Reliability Analysis 236
μ
structure. It is worth noticing that actually SM equals where
1
In a exceptional case where the distribution is Weibul(A,B,C)
with shape parameter C=3.44, the distribution could be
approximated with a normal distribution with parameters
(Carter,1986 as refrence by O'Connor, 2003 )
= + Γ(1 + ) ≅ A + 0.9B , = Γ(1 + Γ( )) ≅ 0.3
237 Reliabilty Engineering
Suppose in a structure
the load is normally distributes with parameters and .
the strength is normal with parameters and
the strength and the load are independent.
Then:
= − ~ ( S , 2 s2 )
= ( > 0)
S S
R Pr( Z ) Pr( Z )
2 2
S S2 2
Since = then
R (SM ) (4-8)
Z
Therefore when the load and the strength are independent and
normally distributed, the reliability(R) is calculated from Eq. 4-
8. The more SM the greater R. More specifically on the average
the more the difference ((strength-load) or the less the variances
of load and strength the more R. Moreover
0.5
1.5
0.25
0.75
1.25
1.75
2.75
3.5
3
4
5
6
1
.9999999
.999968
.999999
.97725
.99865
.99977
.5987
.6915
.7734
.8413
.8944
.9332
.9599
.9937
R=
.50
ϕ ( )
239 Reliabilty Engineering
Fig. 4-3 Plot of Reliability versus Safety Margin( K&L page 80 Redrawn)
0
log(1-R)(log of failure probability per application of load)
-2
-4
-6
-8
-10
-12
-14
-16
0 2 4 6 8 10
safety margin
Fig 4-4 Plot of Logarithm of unreliability versus safety margin(SM)
=0:.25:10;R=normcdf(SM); F=log10(1-R);plot(SM , F)
Table 4.1 gives the reader an idea about the variability in R
related to different magnitudes of variability in normal y
distributed strength and stress random variables(K&L page 79)
Strength Stress SF R=
SM
( )
S S / S
End of Example
Figure 5.4 is a sample plot of log(1-R) versus SM per
application of load. The figure shows if the SM of a design
lies in the third region (i.e. if SM is greater than a threshold), the
logarithm of failure probability is very small and the failure
probability becomes infinitesimal and the design is said to be
intrinsically reliable .
241 Reliabilty Engineering
Eq No.
Dist .of Dist .of
Eq. for reliability Strength( δ) Load(s)
1 λ Exp(λδ ) Exp(λ )
(4-9)
R=
λ λδ
Ref:K&L page 157
2 σ λ Normal Exp( λ)
(4-10)
R = 1 − exp (−μ × λ + ) (μ, σ)
Ref:K&L page 157
3 SM =
μδ μ (4-11) N(σδﻭμδ ) N(σ ﻭμ )
σδ σ K&L page 126
R = ϕ (SM)=normcdf(SM)
Calculable in MATLAB
4 = (z)=normcdf(z) logN(σδ ﻭμδ ) logN(σ ﻭμ )
z=
(4-12)
5 x
s Gamma Gamma
x a 1 (1 x )b 1dx ,
(a b )
R (a, , δ ) (b, , )
(a )(b )
(4-13)
x 0
Ref:K&L page 141
betacdf ( s s , a, b ) : MATLAB
Calculable in MATLAB
6 − Weibull Normal
R = Pr( > ) = + × (A,B,C) (μ, σ)
√2
Ref:K&L page
(4-14)
∞
× ∫ exp [− − 0.5 + ] ,
=
142& Appendix III
Calculable in Matlab , Maple…
7 ∞
Bδ Aδ − A Weib( , Weib( ,
R=1− e exp [−( y δ + ) ]dy , ) , )
B Bδ
In K&L is seen in
C
= d the expression which
seems to be a typo
Chap. 4 Structural Reliability Analysis 242
Solution
F = 1 − R = 1 − Pr(Z < ),
μ −μ δ−μ −
SM = = =
σ +σ σ +0
( )
=( ) ( )
(4-15-1)
=1− e exp [− ( ) ]
[ )]
=1− e ⟹
245 Reliabilty Engineering
c
=1− e , = d ⟹
c [ ]( )
=1− e d
1 c
=1− 1+ e d ⟹
1+
∞ C
− 1+ − 1+ (0)C
e −e ( )C
=1− =1−
( )C + ( )C
1+
( )
⟹R=( ) ( )
End of proof
Notice that
-if C=1 i.e. the load and strength are both exponentially
distributed with parameters λ = and λδ = respectively,
-if the load and strength are both Rayleigh distributed with
scale parameters and respectively, then the reliability is
calculated from Eq. 4-15-1 for C=1:
( )
=( ) ( )
. ( 4-15-2 )
Example 4-3-1
A lognormal distributed load with mean of 60000 Kpa and
standard deviation of 20000 Kpa is applied to a structure which
Chap. 4 Structural Reliability Analysis 246
( ) 20000
= ln + 1 = ln ( + 1) =
( ) 60000
(20000^2/60000^2 + 1) = 0.1054
= ln ( ) − = log(60000)-0.5*.1054= 10.9494
( ) 10000 10000
= ln + 1 = ln ( + 1) = + 1 = 0.01
( ) 100000 100000
= ln ( ) − =log(100000)-0.5*.01= 11.5079
= (z)=normcdf(z) z=
. .
z= = = 1.6441
√ . .
247 Reliabilty Engineering
0.1054
( )= = 10.9494 + = 60000
2
( )= −1 .
approximations.
Example 4-3-2
The following random sample is from the load random
variable applied to a structure. The values are in KPa.
284.9188 104.1661 20.6819 461.9137 197.4067
159.5707 161.5850 50.4525 130.1263 161.5384
418.1608 29.1977 80.9464 7.0621 76.3582
87.0721 16.6974 64.4067 159.6288 39.7292
Solution
Using softwares such as ARENA1 or goodness-of-fit tests or
Q-Q plot help us to consider the load is exponentially distributed
with mean 100KPa and know that a Rayleigh distribution with
mean 111 or equivalently a Gamma(a=2, = 0.018) fits the
strength. Since the exponential distribution could be considered
Gamma(b=1, = 0.01), therefore according to Eq. 4-13 the
reliability of the structure is:
λs
x= λ
s +λ δ
Γ(a+b)
R= x a-1 (1-x) b-1dx=
Γ(a)Γ(b) x=0
λ
betacdf( λs +λs δ ,a,b) MATLAB
x= 0.01 0.01
0.01
x= 0.01+0.018 0.028
x 2-1 (1-x)1-1dx= 2 x(1-x) 0dx=x 2 ] 0.028 =0.1276
Γ(2+1)
Γ(2)Γ(1)
x=0
1×1 0
x=0
or by MATLAB
0.01
R betacdf ( 0.01+0.018 ,2,1) =0.1276
Example 4-3-3
The strength( ) of a component and the stress (S) applied to
it are exponentially distributed with means 150 and 100 psi
respectively. Find the reliability( R) of the component using
Eqs. 4.9 &4-15.
1
ARENAtools-input analyzer- new-file data file- use existing- fit all
249 Reliabilty Engineering
Solution
From Eq. 4-9:
R= = = 0.60
= = , = ⟹d = d ⟹
1
=1− e
instructions results in R= .
>>bd=150;bs=100;syms x;W=[(1/bd)*exp(-(x/bd)-x/bs)]
>>R=1- int(W,x,0,inf)
S
LR
2 2
(4-16)
S
where
1
Wu,Y. Xi,L. 2010 Load-roughness impact on reliability considering
dependent failure ,Proceeding 16th ISSAT conference on Reliability and
Quality in Design
251 Reliabilty Engineering
answer in part a.
Solution
a)
5000 − 3500
= ≤ = ⎛ ≤ ⎞ = 0/99598
400 + 400
⎝ ⎠
b)
L=norminv(1-0.99598,5000,400)= 3939.85N
=F δ = ( ) = =1−
0.8 = 1 − ⟹ = 0.2
= Pr S = ( ) = 1 − Fδ (s) (4 − 19)
Example 4-7
A component is subject to the fixed load of 4000N. The
strength is log-normally distributed with mean of 5000N and
Chap. 4 Structural Reliability Analysis 258
Solution
ln[var(X ) / E 2 (X ) 1]
using MATLAB
σ=sqrt(log(4002 /50002 +1))=0.1
2
μ=lnE(X)- σ =ln(5000)- 0.01 =8.5122
2 2
R=Pr δ>S=4000 =
Pr(lnδ>ln4000=8.294)=Pr(Z> 8.294-8.5122
0.1 =-2.182)=0.9854
or using MATLAB:
End of Example
259 Reliabilty Engineering
×
R≥1− ( )
(4-21)
where
n The mean of SF
V = V =
n
Chap. 4 Structural Reliability Analysis 260
S2
1 2
n (4-22-1 )
S S
2
2 2 4
2 S S
SF
S (4-22-2 )
S S
Vn(%)
Fig 4.13 Plot of n = versus Vn for six levels of R
n S
=Φ ( ) (4-23)
2 2
s
Solution
From Fig. 4-21 the minimum of the reliability is 0.95.
using Inequality 4-21:
σ 0.8357
= = = 0.18
μ 4.6426
= 0.95.
. ∗ .
≥1− ( . )
. ∗ .
End of Example
Chap. 4 Structural Reliability Analysis 262
Example 4-8-1
The reliability of a structure should be at least 0.95 . the
coefficient of variation is 18%. How much is the safety factor
on average?
Solution
End of Example
R Pr S f ,S , S d dS
S
or equivalently, in the case of independence of stress and
strength, the failure probability cannot be computed from:
R Pr(S ) F (s )f s (s )ds [1 Fs ( )]f ( )d (4-24)
Where
as the upper limit for S and the lower limit for respectively
forming the interference interval [ min , S max ]. min and S max are
Step 1
from:
s max s max
R=Pr(s>δ) Fδ (u)f s (u)du [1 FS (u )]f (u )du (4-25)
δ min min
Step 2
Step 3
p i P r(a i 1 S a i ) i 1, . .. , n (4-25)
q i P r (a i 1 a i ) i 1, . . . , n (4-26)
Now it can be shown that Eq. 4-25 could be approximated
by(K&L page89):
n
i n n
R pi q k qi pk (4-27)
i 1 k 1 i 1 k i
Step 4
Step 5
n
Pr(s min ) p i ( 4-29-2)
i 1
n n
ap pi b p , aq q i bq
i 1 i 1
Step 6
n
i n n
Max/Min R p i q k q i p k
i 1 k 1 i 1 k i
s .t .
L p ,i p i U p ,i i 1,..., n
L q ,i q i U q ,i i 1,..., n
n
ap p i b p ,
i 1
n
aq q i bq
i 1
pi 0 qi 0
End of algorithm.
Find the lower and upper bound for the reliability of this
structure.
Solution
32-30
p1 Pr(30 s 32) Pr(0 Z ) 0.2475
3
269 Reliabilty Engineering
Step 4
Upi=pi+(0.02)*pi ; Lpi=pi-(0.02)*pi
Uqi=qi+(0.02)*qi ; Lqi=qi-(0.02)*qi
-10 -10
Up10= 9.7350 10 +(0.02)* 9.7350 10 ;
-10 -10
Lp10= 9.7350 10 - (0.02)* 9.7350 10 ;
As sample calculation:
q1=0.0011;Uq1=0.0011+(0.02)*0.0011;Lq1=0.0011-(0.02)* 0.0011
n
30-30
p
i 1
i Pr(s min ) Pr( Z
3
) 0.5
n
q i Pr( s max ) 1 exp[ (
50-30 2
) ]=0.105.
i 1 60
n
0.02 0.02
0.5-( ) 0.5 p i 0.5+( ) 0.5
2 i 1 2
n
0.02 0.02
0.105-( ) 0.105 q i 0.105+( ) 0.105
2 i 1 2
Chap. 4 Structural Reliability Analysis 272
s .t .
n
ap pi b p 0 .5-(0.02/2)*0.5 p1 ... p10 0.5+(0.02/2)*0.5
i 1
n
aq q i b q 0.1055-(0.02/2)*0.105 q1 ... q10 0.105+(0.02/2)*0.105
i 1
0 pi 1 0 pi 1
Softwares such as Lingo or GAMS give the following results for the
maximization and minimization :
The objective function has the optimal values:
End of Example
SM = − 1. (p-1)
Therefore f mw = and:
SF
SM = −1 (p-2) or
×
SF = ( p-3)
( + )×
where
SM The safety margin or margin of safety
SF The safety factor
The strength of the material
The maximum allowable working stress
The stress applied to the structure.
Exercises
In the following problems the stress and the load(stress) are independent.
1-(Problem 5 Page159 K&L) The strength and the stress S for the
design of a component are logomrally distributed with the
following infonnation on and S:
Find R.
s 2 s2 s 2 s2
R Z ( ) exp( s ) [1 Z ( )]
s 2 s
Chapter 5
on the
Combination of
Random Variables
in Design
& a Glance at the
Tolerance Concept
Chap. 5 Combinations of random variables + tolerance 282
5
On the combinations of random variables in design ;
A glance at the tolerance concept
5.1 Introduction
M c
S (5-1-1)
I
283 Reliabilty Engineering
where
S = the maximum stress at the farthest surface from the neutral
axis (it can be at top or at bottom), kPa
ℎ( ) = × [ ( )] (5-2)
where
ℎ( ) = × ( )+ × ( ) (5-3)
Example 5.1
Random variable X is normally distributed with density
2
x
function(pdf) 1 2 2 and = . Find the pdf of Y.
f (x ) e
2
Solution
= ( )= ,
pdf of Y = ℎ( ) = × [ ( )] , ( )= ( )
1
= ⟹ = = ( ) =
2
ln y
1 2
f (k ( y )) e 2
2
2
1 1
ln y
ℎ( ) = × [ ( )]= × e
2
;
2 2
Solution
4A 4A
= ⟹ D1 ,D2 ;
4
ℎ( ) = × ( )+ × ( )
4A 4A
ℎ( ) = × + ×
1 4 4 1 1 1 dD 1
D A 2 | dD
dA
| 2 A 2 | |
dA A
2 2
4A 4A
1 1
h (A ) [exp( 2
) exp(
2
)] .
A 2 2 2
End of Example
Chap. 5 Combinations of random variables + tolerance 288
E[X+Y]=E[X]+E[Y] or +
= 1
+ 2
(5-4)
End of Proof
= + −2 ( )= + −2 . (5-6-2)
289 Reliabilty Engineering
= + . (5-7)
( ) ≅ [ ′( )] . (5-9)
Solution
The cross section of the is calculated from:
= ( )= ′( ) = 2 "( ) = 2
From Eq. 5-8:
( . )
( )≅ ( )+ "( ) = (2 ) + (2 ) = 4.01
From Eq. 5-9:
End of Example
Chap. 5 Combinations of random variables + tolerance 290
( ) ≅ ( ,…, )+ ∑ ⃒ ( ) (5-10-1)
Or in vector form:
|
⎛ . ⎞
( ) ≅ ( ,…, )+ ( )…. ( ) ⎜ . ⎟ (5-10-2)
⎜ . ⎟
⎝ | ⎠
⎧ ⎫
⎪ ( ,…, ) ⎪
( )≅∑ × ⃒ (5-11)
⎨ ⎬
⎪ ⋮ ⎪
⎩ ⎭
Solution
= = ( , ),
⃒
1 ⎛ ⎞
( )≅ ( , )+ ( )⎜ ⎟
2
⃒
⎝ ⎠
2 2
2 20000
=− , = . ⃒ = = 160
5
1 10000
= , =0. ( , )= ⟹
5
2 0
10000 1(0.4 0)
160 ⟹
( )≅ +
5 2
( ) ≅ 2000 = 20000 = 20
( )≅ ( ⃒ ) + ( ⃒ )
1 −10000
= (1000) + (0.4)
5 5
End of Example
approximately :
2
X X Y
E 1 (5-12-1)
Y Y Y
2
1 1 Y (5-12-2)
E 1
Y Y Y
4
2 2 2
X X
x Y Y (5-12-3)
V ar
Y Y X Y Y
293 Reliabilty Engineering
Example 5.5
Solution
2
P A
1
10000 0.4
E( ) = ≅ = ∗ 1+ = 2012.8 N / cm 2
5 5
A A
⟹ E( ) = 20128
4
2 2 2
P P
P A A
( ) = Var
A A P A A
10000 1000 0.4 0.4
= ∗ + − = 65436
5 10000 5 5
tolerance is 0.12
100 2.500= 0.003.
part
2
sum (5-13-3)
n
( ) ( )
Let ∆= and = then:
∆
= (5-15-1)
√
∆= √ (5-15-2)
∆′ = 2
. (5 − 16)
=1
Chap. 5 Combinations of random variables + tolerance 296
Example 5-6
Solution
Example 5-7
Solution
.
The specification of each of the 10 parts must be ±
∆ .
where = = = 0.012 . End of Example
√ √
Example 5-8
297 Reliabilty Engineering
Solution
Solution
Taylor series of a 3-variable function up to the first order is
as follows:
1
f x , y , z f a,b,c x a y b z c f x , y , z ...
1! x y z x a, y b, z c
(5-17)
1
E 0 E 1*N *K * E 1 E 1* N *K * N N * E 1* K * K K * E 1* N *
1!
*
E 0 N *K *E 1 E 1 K * N E 1*N *K 2E 1* N * K *
E E N *K *E E *K *N E *N *K 2N *E *K *
0 1 1 1 1
* *
E N * K * EE E 1* K * NN E 1* K * N * 2N * K * E 1*
0
1
1
or Eo =E 1* N * K *
E N *K *E E *K *N E *N *K 2N *E *K *
0 1 1 1 1
2 2 2
2 * * 2 * * 2
E N K E E1 K E1*N * K2
(5-18)
0 1 N
σ = (5-19)
where z is the critical value related to normal distribution
given in Table D or by a software such as MATLAB; e.g. for
=0.27% Z z0.00135 norminv(1-0.000135) 3 and
2
therefore:
0.5 1
×0.01 3×0.02
E = , N = 2
= 0.0017 , k = =0.02
1
3 3 3
Substituting numerical values in the right hand side of Eq. 5-18
gives E 0.512 . * * *
Eo E1 N K 40× 12 ×3=60 is the nominal value
0
Δ Δ Δ Δ
0.512 = = = = ⟹ Δ = 1.536
. 0.00135 3
1.536
=2.56% E 0 : 60±2.56% .
60
: 60 ± 1.536 or since
R Pr(S ) Z (z ) (5-20)
where
= and
( . )
The diameter of the bar has a mean of r 0.12635 inches and its
specification is r p % . The load and the strength are normally
distributed. To know how the variations in the rod diameter affect the
rod reliability, conduct a sensitivity analysis of the rod reliability
with respect to the rod radius.
Solution
1
For proof see K&L page 165
303 Reliabilty Engineering
p z R Z (z ) normcdf z
0 3.760 0.999915
0.5 3.756 0.999914
1.0 3.74 0.999908
1.5 3.72 0.999900
3.0 3.61 0.999847
5.0 3.37 0.999624
7.0 3.10 0.999032
If for example the specification of the radius of the rod is
0/12635±1.5% i.e. p=1.5, z turns to be z=3.72 and the rod
reliability would be 99.99% as calculated below:
p=1.5;
z=(10^5-
(4000/(pi*(0.12635)^2)))/sqrt(5000^2+((100^2+((4*(0.01*p)^2)
/9)*4000^2))/(pi^2*0.12635^4))
R= normcdf(z)
End of Example.
Exercises1
In the following problems assume all dimensions are normally distributed and the tolerance
range is 6 sigma(±3σ)
1
Problems 1 through 7 are from Chap. 5 page 113 problems(1,3,5,7,9,11,13)
Chap. 5 Combinations of random variables + tolerance 304
specification of V=XYZ.
Variable
mean 4 2 1
stand. devia 0.4 0.2 0.1
Determine the value which the random variable smean will exceed
only 1.3% of the time. Also determine the value which the
random variable samp will not exceed 90% of the time.
Hint: For largish n a gamma distribution (n, ) could be
√
approximated with a normal distribution ( = , = );
The four forces F1, F2, F3, and F4 are random variables, the
distributions of which are given in the following table
Chap. 5 Combinations of random variables + tolerance 308
R 1 100 , R 1 10 , R 2 200 , R 2 15 ,
Chapter 6
Estimation of
Mean Life &
Reliability,
Related
Experiments
&Tests
311 Reliabilty Engineering
6
Estimation of Mean Lifetime &Reliability and
Related Experiments &Tests
6.1 Introduction
The problem of estimation of the lifetime and the reliability
of products is a common problem in the control of products
quality. When we have the lifetimes of a random sample of the
product, one obvious way to estimate the mean lifetime is
calculating the sample mean. Another way is performing special
life tests on the sample and then calculating the mean. A third
way of life testing is called accelerated life testing which
involves the acceleration of failures to quantify the life
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 312
( 6-1)
∑
θ= =
This equation could be used for any product with any lifetime
distribution.
6.3 Tests for Estimating Mean Life
Consider a life testing where n items are simultaneously
placed on test. The purpose of the life tests here is to calculate a
point estimate and sometimes interval estimates for the product
mean life. It often occurs that we need to discontinue the life test
before all the elements in the sample fail. In such cases, we say
that the test has been “suspended,” “censored,” or “truncated”
.Censoring schemes employed during the life test make the
inspection as a cost effective one. Time censoring (Type I),
failure censoring (Type II) are 2 common types of the censoring
313 Reliabilty Engineering
ℎ
⎧ :
⎪ ℎ
⎨ ℎ
⎪ :
⎩ ℎ
Estimate of
θ= (6-2)
where
T = the total operation time of all items placed on test including
those failed,
r = total number of failures occurred during the life test
315 Reliabilty Engineering
Calculation of T for = = =
∗
In type-I censoring at a predetermined time, say , the life
test is terminated. The test could be performed without or with
replacement.
T nt *
ˆ (6-3)
r r
∗
Where r ≥ 1 is the number of failures dung the test time .
An application of this equation is , for example, when we have
where we have n test stands, and we cycle each
test stand for cycles. As items fail they are replaced. Where a
truncation time is specified this is called Type I censoring. Here
we have(K&L page 252)
317 Reliabilty Engineering
n
ˆ (6-4)
r
Example 6-2(Example 10-12 K&L page 252 based on Example 10-3 page 241)
Nine stands are used for testing the life of a kind of switch.
As items fail they are replaced. Each stand was cycled 20,000
times, and a counter recorded the cycle number at which failures
occurred. The following table contains the data. Estimate MTBF.
Stand no. 1 2 3 4 5 6 7 8 9
4100, 18100, 18950
No
Cycles at No failure
3100, 8100
which
6700
4600
5400
2600
4700
failure
failure
occurred
Solution
= n = 9 × 20000 = 180000
T 180000
M TˆB F = θ = ˆ = = 18000
r 10
Solution
T 24×5000=120,000 hr
*
ˆ T nt 120000 =8571hr .
MTTF
r r 14 End of Example
where
Example 6-4:
Solution
( ) ( ) ∑
= = (6-6)
It is obvious that:
321 Reliabilty Engineering
1. If the time of the life test is such that all n items fail then
∑
r=n and Eq. 6.1 i.e. θ = is obtained.
Example 6-5
Solution
r
T (n r ) x( r ) xi (20-10)×41.2+27.4+...+14.4+41.2=624.5
i 1
Then according to Eq. 6-6
nx
(r ) (6-8)
r
T nx ( r )
MTˆBF ˆ (6-9)
r r
(Type-I Replacement ( n r )t * x ( i )
T i 1
censoring) ˆ
r r
Without 6-5 ( − ) +∑
( )
Replacement =
Failure- Without 6-6 nx
terminated Replacement (r )
(Type-II r
censoring) With 6-8
Replacement
T nt *
ˆ
r r
2T 2T (6-11)
Pr 2 2 1
2 r ,2 2 r ,1
2
Example 6.7
Solution
r 8
T=(n-r)x (r) + x i =(8-8)x (8) + x i
i=1 i=1
8 8
T x i = x (i) =669320
i=1 i=1
∑ ( )
= = = = 83665
2T 2T
2
θ 2 .
χ α χ α
2r, 2r,1-
2 2
2
χ16,0.975 chi2inv(.025,16) 6.91,
2
χ16,0/025 chi2inv(0.975,16) 28.85
2T 2T 2×669320 2×669320
2
θ 2 θ
χ α χ α 28.85 6.91
2r, 2r,1-
2 2
Example 6.8:
2T 2T
Pr 2 θ 2 =1-α ,T= r×θ=20×5000, 1-α=0.9
χ α χ α
2r, 2 2r1 ,(1- )
2
Solution
The test is of Type II without replacement, therefore:
r
T (n r ) x( r ) xi , r =1,x (1) =15.5 ,
i 1
r
T=(n-r)x (r) + x i =(n-1)x (1) +x1 =(n-1)x (1) +x (1) =nx (1)
i=1
1
from ttps://web.cortland.edu/matresearch/CensorDatSTART.pdf
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 330
1-α=0.95 α=0.05,
MATLAB
χ 22,0.975 =chi2inv(0.025,2)=0.0506,
χ 22,0/975 = chi2inv(1-0.025,2)=7.3778
2×5×15.5 2×5×15.5
n=5,x (1) =15.5 θ 21 θ 3063
7.3778 0.0506
End of Example
estimate for : = × ln
where
= is when (the time or the kilometer or the temperature
or…) that fraction p of the working product fail or fraction R=1-p
of them do not.
Last part of problem 6 of this chapter exercises uses Sec.
6.6.2.
6-7 Reliability Acceptance Sampling Plans
During the past years several researches have been done on
the subject of sampling from a lot of products to accept or reject
it based on product lifetime. Single, double and multiple
sampling plans have been developed in this regard which are
called Reliability Acceptance Sampling Plans (RASP) ,utilized
to inspect the quality of a lot for acceptance.
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 332
1
H108., Quality Control and Reliability Handbook (Interim) Sampling
Procedures and Tables for Life and Reliability Testing (Based on Exponential
Distribution),, in (Supply and Logistics)
2
https://opentextbc.ca/oerdiscipline/chapter/industrial-engineering or
https://archive.org/details/statistical-methods-august-2020-bazargan or
https://opentextbc.ca/oerdiscipline/chapter/statistics/
333 Reliabilty Engineering
=1− (1 − ) (6 − 14)
= (1 − ) (6 − 15)
Example 6.10
Design a plan whose test is of type "life tests terminated at
pre-assigned time" in such a way the test time does not exceed
T=500 hr . The plan is wanted to accept 90% the lots having
mean life θ =10000 hr (α = 0.10), and to reject 95% of the
lots with mean life θ = 2000hr (β = 0.05). The life is
assumed to be exponentially distributed.
Solution
= = , = =
Table 6-2 Acceptance Sampling Plans for some , , & with a test
terminated at pre-assigned time with replacement (Table 2C-4 in H108 Handbook)
T / 0 T / 0
1 1 1 1 1 1 1 1
335 3 5 10 20 3Reliabilty
5 Engineering
10 20
/
r n n n n n n n n
1 0 r
= 0.01 = ﻭ0.01 = 05 = ﻭ0.0.01
2thirds 136 3 551 1103 220 9 23 397 795 1591
1half 46 9 158 317 634 3 72 120 241 483
1third 19 3 51 103 206 1 25 38 76 153
1 fifth 9 1 17 35 70 7 9 16 32 65
1tenth 5 4 6 12 25 4 4 6 13 27
α = 0.01ﻭβ = 0.05 α = 0.05ﻭβ = 0.05
2thirds 101 2 395 790 158 6 16 270 541 1082
1half 35 6 113 227 454 2 47 78 157 314
1third 15 2 37 74 149 1 16 27 54 108
1 fifth 8 8 14 29 58 5 6 10 19 39
1tenth 4 3 4 8 16 3 3 4 8 16
= 0.01ﻭ = 0.1 = 0.05 = ﻭ0.1
2thirds 83 1 316 632 126 5 13 216 433 867
1half 30 5 93 187 374 1 37 62 124 248
1third 13 1 30 60 121 8 11 19 39 79
1 fifth 7 7 11 23 46 4 4 7 13 27
1tenth 4 2 4 8 16 3 3 4 8 16
= 0.01ﻭ = 0.25 = 0.05 = ﻭ0.25
2thirds 60 1 217 434 869 3 77 129 258 517
1half 22 3 62 125 251 1 23 38 76 153
1third 10 1 20 41 82 6 7 13 26 52
1 fifth 5 4 7 13 25 3 3 4 8 16
1tenth 3 2 2 4 8 2 1 2 3 7
= 0.1 = ﻭ0.01 = 0.25 = ﻭ0.01
2thirds 77 1 329 659 131 5 14 234 469 939
1half 26 5 98 197 394 1 42 70 140 281
1third 11 2 35 70 140 7 15 25 50 101
1 fifth 5 7 12 24 48 3 5 8 17 34
1tenth 3 3 5 11 22 2 2 4 9 19
α = 0.1 = ﻭ0.05 α = 0.25 = ﻭ0.05
2thirds 52 1 214 429 859 3 84 140 280 560
1half 18 3 64 128 256 1 25 43 86 172
1third 8 1 23 46 93 5 10 16 33 67
1 fifth 4 5 8 17 34 2 3 5 10 19
1tenth 2 2 3 5 10 2 2 4 9 19
α = 0.1 = ﻭ0.1 α = 0.25 = ﻭ0.1
2thirds 41 9 165 330 660 2 58 98 196 392
1half 15 3 51 102 205 8 17 29 59 119
1third 6 9 15 31 63 4 7 12 25 50
1 fifth 3 4 6 11 22 2 3 4 9 19
1tenth 2 2 2 5 10 1 1 2 3 5
α = 0.1 = ﻭ0.25 α = 0.25 = ﻭ0.25
2thirds 25 5 94 188 376 1 28 47 95 190
1half 9 1 27 54 108 5 10 16 33 67
1third 4 5 8 17 34 2 2 4 9 19
1 fifth 3 3 5 11 22 1 1 2 3 6
1tenth 2 1 2 5 10 1 1 1 2 5
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 336
= Pr( ≤ ) = ∑ (1 − ) (6-16)
b)the lifetime
time is exponentially distributed then for T=100⟹
T=100
6-7-3-2
2 Operating Characteristic curve for single sampling
plans (Pa versus mean lifetime)
Example 6-12
6 (OC curve:Pa plotted versus mean lifetime)
End of example
6-8
8 statistical hypothesis testing on mean and
minimum lifetime
n x(1)
Ho is rejected if F0 F ,2,2r 2 (6 17)
ˆ '
Cycles to
190437
245593
277761
432298
530100
626300
1043307
1055528
1221393
2099199
failure
x(i)
Solution
H0 : 0
H1 : 0
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 344
10
(n-r)(x (r) -x (1) )+ (x (i) -x (1) )
i=1 (15-10)(2099199-190437)
θ= = +
r-1 10-1
=1908762
(190437-190437)+(245595-190437)+....+( 2099199-190437+ )
θ=2767421
10-1
nx (1) 15×190437
F0 = = =1.03
ˆ
θ' 2767421
H 0 : 0
H1 : 0
8. Reject Ho if
2 2r
0 2,2r (6-18 )
0
Example 6-14
Solution
2 2r
0 r 10, n 15, 5%
0
The test is of non-replacement failure-terminated type, therefore
according to Eq. 6.6
(n r )x ( r ) x i (2)(10)(2871391)
2871391, 02 57
r 106
S = ( ), ( ), … , ( )
. For simplicity S will be
where
x ,x
2(1) 1(1) The occurrence time of first failure in Design 1 & 2 respectively
Rejection Criteria:
Reject 0 if is outside
For H1 :1 2
F α F ,2 r2 2,2 r1 2
1 2 ,2r2 -2,2r1 -2 2
For H1 :1 2 Reject 0 if F0 >F ,2r1 -2,2r2 -2
Exercises1
1
Exercises 1 through 8 are from Chap. 9 K&L page 269 problems
1,2,3,4,6,8,9,10
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 350
10. Plot1 the OC curve (Pa versus MTBF) for the following
sampling plan:
Twenty units are randomly taken from a largish sample, and
simultaneously placed on lifetest. Whenever an item fails, it
is replaced with another item selected at random from the
lot. If the test continues for 500 hr with not more than 2
failures accept the lot. If 3 failures occur before the 500 hr of
testing, reject the lot and terminate the test.
1
Solution on Page 586 Grant &Leavnworth(1988)
Chap 7 Dynamic Models+Availability, Markov Chain 354
Chapter 7
Dynamic
Models,
Availability,
Application of
Markov Chain to
Reliability
355 Reliability Engineering
7
Dynamic Models+ Availability, Application of
Markov Chain
( )=∏ ( ) (7-1)
Proof No.1
' m m
-R sys (t)
h sys (t)= R sys (t ) R i (t ) [1 Fi (t )]
R sys (t) t 1 t 1
m m
m
m
R (t)= f i (t ) R j (t ) f i (t ) R j (t )
'
sys
i 1 j 1
i 1 j 1
j i
j i
Therefore
m m
f (t ) R (t )
'
i
i 1 j 1
j
R sys (t) j i
h sys (t)=
m
R sys (t)
R i (t )
t 1
m
f 1 (t ) f m (t ) f i (t ) m
hsys (t)= +...+ = hi (t )
R1 (t ) R m (t ) i 1 R i (t ) i 1
357 Reliability Engineering
Proof No.2
h ( ) d
d
R (t ) e 0
h ( )d ln R (t ) , h (t ) ln R (t )
0
dt
R sys (t ) R 1 (t ) ... R n (t ) ln R sys (t ) ln R i (t )
d d m
ln R sys (t ) ln R i (t ) hi (t )
dt dt i 1
m
d
hsys (t ) ln R sys (t ) h i (t ) hsys (t ) hi (t )
dt i 1
End of proof
∫ ( )
( )= =
( ) = exp − ∑ − ∑ (7-2-1)
Let ∗ = ∑ , ∗
=∑ , = ∑ = ∗
∑ ∗
lim →∞ ( )= = = (7-2-2)
Chap 7 Dynamic Models+Availability, Markov Chain 358
Example 7-1
n n
h sys (t)= h i (t)= λ i = constant
i=1 i=1
359 Reliability Engineering
-t λ
R sys (t)=e i (7-4-1)
1
1 1 n 1
M T BF n
n
(7-4-2)
1 i=1 i
λ
i=1
i
i=1 M T BFi
1
where MTBFi i . End of Example
λi
-n t
R sys (t)=e (7-4-3)
1 MTBFpart
MTBF (7-4-4)
n n
= ( ,1, … , ) ( 7-5)
where Rj is the reliability of the ith component and describes the
component's resistance to failure from thermal stresses.
Ri s 0 s s f ( )d ds
f ( s ) f s (s ) 1 F (s ) ds (7-7)
0
Fig 7-1 The region where the strength is more than stress
n
C min i
i 0
n (7-8)
FX (1) ( y ) 1 1 FX ( y )
Chap 7 Dynamic Models+Availability, Markov Chain 362
Where
n
1-FC (y)=1- 1-1-Fδ (y)
we could conclude that the reliability
of the n-link chain which equals RC=Pr( >s) is obtainable from:
n
RC f S (s ) 1 F (s ) ds (7-9)
0
Where
RC The reliability of an n-link chain
Note the similarity of Eq. 7-7 with Eq. 7-9 for n=1.
363 Reliability Engineering
Example 7-2
Solution
t t
1 600
f part (t ) e , F part (t ) 1 e 600 ,
600
s
1 500
The pdf of the stress applied to the system is f e
s 500
According to Eq. 7-9:
n
R sys f (s ) 1 F part (s ) ds n 10
0 s
s s 10
1 500 600
R sys e e ds
0 500
56 S
1 3000 1 3000 6
R sys e ds
500 s 0 500 56 56
End of Example
Chap 7 Dynamic Models+Availability, Markov Chain 364
Proof
1
From: , Lewis(1994) page 263, Grosh(1989 )page169,Li(2016)and
https://www.weibull.com/hotwire/issue21/
Chap 7 Dynamic Models+Availability, Markov Chain 366
( ≤ )= [ ( ,…, )≤ ]= ( ≤ ,…, ≤ )
( ≤ ,…, ≤ )= ( ≤ )… ( ≤ )
( ≤ )= ( ≤ )… ( ≤ )
⇒
1− ( ) = [1 − ( )][1 − ( )] … [1 − ( )] ⇒
( ) = 1 − [1 − ( )][1 − ( )] … [1 − ( )]
n
R sys (t ) 1 (1 e t i ) t 0 (7-11)
i 1
( ) = 1 − (1 − ( ))(1 − ( ))
( )=1− 1− 1− =
1 1 1
(MTBF ) sys Rsys (t )dt (7-13)
0 1 2 1 2
( ) = 1 − (1 − ( ))(1 − ( ))(1 − ( ))
( ) = 1− 1− 1− 1− (7 − 14 − 1)
MTBFSys R sys (t )dt
0
1 1 1 1 1 1 1 (7 − 14 − 2)
+ + - - - +
λ1 λ 2 λ 3 λ1 +λ 2 λ1 +λ 3 λ 2 +λ 3 λ1 +λ 2 +λ 3
n
1
MTBFSysActive (7-15)
k 1 k
n
1
Table 7-1 gives the values of k Mn .
k 1
3 1 1 11 11
n=3 MTBFsys = -3× + = = θ .
λ 2λ 3λ 6λ 6
369 Reliability Engineering
R s y s ( t ) = 1 -( 1 -e -λ t ) n (7-16)
It is seen from Table 7-1 that "for a system to have double the
mean life of a single component, it must consist of 4
components. To triple the mean life the system must have 11
components. Theoretically, there is no limit to how much the
system mean life can be extended but the cost of extending life
through mere redundancy is usually prohibitive. Redesign
should be performed to excessive redundancy"(Grosh,1989 p 135).
Given the reliability function of an active parallel system,
R active (t ) , the failure rate function is obtained from:
Chap 7 Dynamic Models+Availability, Markov Chain 370
d
Ractive (t )
hactive (t ) dt (7-17)
Ractive (t )
Switch
n
2 t d
R (t)=R 1 (t)- R 2 (t-t ) R 1 (t ) dt =
standby 0
dt (7-18)
t
=R1 (t)+ R 2 (t-t )f1 (t )dt
0
Where
1
Pradip Kundu & Asok Nanda, Redundancy Allocation in a System: A Brief Review
Chap 7 Dynamic Models+Availability, Markov Chain 372
Example 7-3-1
active unit fails a perfect switch replaces it with the other unit.
Calculate the reliability of the system.
Solution
2 t d
R (t)=R1 (t)- R 2 (t-t ) R 1 (t ) dt =
standby 0
dt
t (t-t ) d '
e 1t e 2 e 1t dt '
0
dt
t (t-t ) t
=e 1t e 2 1e 1t ' dt ' e 1t 1e 2t e( 2 1 )t 'dt '
0 0
2
1
R (t) = e-λ1t + λ1e-λ2 t (e(λ 2 -λ1 )t -1)
standby λ 2 - λ1
373 Reliability Engineering
1 251 2 101
2
R (20)=exp(-20/10)+(1/10)*exp(-20/25)*(exp(20/(1/25-1/10))-1)/(1/25-1/10) 0.8842
standby
End of Example
Assume that the main and standby units are identical, each
with a constant failure rate and the switch is perfect. Then:
2
R Sys (t ) e t (1 t ) t 0 (7-19)
Proof
Since R1(t)= R2(t)= e t then according to Eq. 7-18:
2 t d
R (t)=R1 (t)- R 2 (t-t ) R 1 (t ) dt =
standby 0
dt
t
e t e (t t ') -e t ' dt =e t te t e t (1 t )
0
End of Proof .
Example 7-3-2
Calculate and compare the reliability of two configurations
of a two-unit system (active parallel and standby parallel). The
failure rate of the unit is 5%.
Solution
Chap 7 Dynamic Models+Availability, Markov Chain 374
2
R (t) et (1t) t 0
standby
2
R sys (t ) 1 (1 e t i ) t 0 1 2
i 1
5
Then for t 10 , 100
2 2
R 91%, R 84.5%
standby active
2 2
5
for t 100 , 100 ,R 4%, R 1.3%
standby active
where
375 Reliability Engineering
3
R The reliability of 3-unit standby system
standby
2
R The reliability of 2-unit standby system
standby
2
3
t
t
R (t ) e 1 t t 0 (7-21)
sys
2
Example 7-4
The reliability of the water supply system of a city is a
concern of the city council. The council would like to ensure
the inhabitants that the system will work for 20 years with a
reliability of 95%. At the time being the water is supplied by a
reservoir and a river ( with mean lifetime of 25 and 10 years
respectively) in parallel. The water is then disinfected in a
building that has an active unit for disinfection and 2 standby
backup units. Each of the disinfection units is designed for a
useful life of 25 years. After disinfection the water goes to
distribution subsystem, which is 99% reliable. The council is to
decide whether to allow the municipality to add a deep well to
Chap 7 Dynamic Models+Availability, Markov Chain 376
Solution
Fig 7.3 The RBD of the wateter supply system of Example 7-4
Data: = 20 , = , λ =
= , = 0.99
2 202
t 20 20 2
R (t ) e t 1 t R (20) e 25 (1 25 25 ) 0.932
2 2
The entire system reliability is currently:
n 1 i
n ( t )
R (t ) e t i 0 (7-22)
standby
i!
n
- R standby '(t)
h standby (t)= n
(7-23)
R standby (t)
Chap 7 Dynamic Models+Availability, Markov Chain 378
Figure 7.4 shows "both the reliability and the failure rate for
[2-component] active and standby parallel systems, along with
the results for a system consisting of a single unit with the
assumption that the lifetimes are exponentially distributed. At
intermediate times the failure rate for the standby system is
smaller than for the active parallel system. This is reflected in a
larger reliability for the standby system(Lewis,1994 page 256).
1 dR active (t )
hactive (t )
R active (t ) dt
1 e t
hactive (t ) t
(7-24)
1 0.5e
t
hstandby (t ) (7-25)
1 t
MTTF R (t )dt . (7-26)
0
n
MTBFstandby = =n MTBFpart (7-27-4)
381 Reliability Engineering
Proof
n
E (X standby ) E ( X i ) E (X i ) MTBFi
i 1
1
MTBFsys (7-27-5)
i
n
If 1 ... n then MTBFstandby = =n MTBFpart and the
proof is complete.
Example 7-5
Notice that if both units were active the according to Eq. 7-10:
2
R (100)=1-(1-e -λt ) 2 = 0.9909 .
active
t
f sys (t ) f 1 (z )f 2 (t z )dz (7-28).
z 0
where
383 Reliability Engineering
( )= ( − ) (7-29)
t
f sys (t ) f 1 (z )f 2 (t z )dz
z 0
( ) ( )( )
( )= =
( )
= 1− = + ⟹
( )= ( − )
2 1
1e t 2e t
R sys (t) (7-30)
1 2 2 1
Chap 7 Dynamic Models+Availability, Markov Chain 384
Proof:
t t
Rsys (t)=R 1 (t)+ R 2 (t-t ) f1 (t ) dt =e 1t e 2 (t t ') 1e 1t ' dt
0 0
2 t 1t
1e e
Rsys (t) 2
1 2 2 1
Notice that for = = use Eq. 7-19 i.e. R sys (t ) e t (1 t )
(7-33)
MTTFsys xf sys (x )dx Rsys (t )dt
x 0 0
where
n n
1
Remember as proved earlier MTTFsys MTTFi gives
i 1 l 1 i
where
n
Let R (t ) denote the reliability function of a system having 1
sys
2 t
R (t)=R1(t)+Ps f1(t')R2 (t-t)dt (7-37)
sys 0
where
2
R (t ) Reliability function of 2-component standby system
sys
2
R (t)=[1+(1-p)t]et et (1 Ps t) (7-39)
sys
where
2 t
n 2 R (t) R1(t) f1(t) Rs (t ') R2 (t t)dt (7-40)
sys 0
where
389 Reliability Engineering
2
R (t ) Reliability function of 2-component standby system
sys
Let failure rates of the active unit, the cold redundant standby
unit and the switch is constant and equal 2 , 1 & then:
s
t 2 1 1 1
MTTFsys R (t )dt
0 sys 1 (1 2 )(1 ) (1 2 )2
s s s
or
1 1 1 1
MTTFsys ( ) (7-41)
1 1 2 1 2
s s
Furthermore if 1 2 then :
1 1 1
MTTFsys-standby ( )
s s
1 1
MT TFsys-standby (7-42)
s ( )
s s
where
2 t
R (t ) e t 1 1 e s t 0 (7-43)
sys s
2
MTTFsys R (t )dt e t 1 1e s t dt
sys
0 0 s
1 1
MTTFsys which is Eq.7-42.
( )
s s s
2
1+P s
M TTFsys R sys (t )dt (1+P s ) M T T F part (7-45)
0
Notice:
391 Reliability Engineering
2
-Substituting
Substituting P s 1 in the above relationships yields R (t ) and
sys
-The
The reliability function of a two-component
two component system having similar
active and standby units with the same failure rate λ is :
2
Calculated from Eq. 7-43
7
i.e. R (t ) e t 1 st given
sys
1 e
s
λs as the switch failure rate
2
Or is calculated from Eq. 7-39
7 i.e. R (t ) e t (1 Ps t ) if the
sys
7 1
Example 7-6
1
From Dr Eshargh's pamphlet, Faculty member of Sharif University of Tech
, Tehran. Figure from K & L p 59.
Chap 7 Dynamic Models+Availability, Markov Chain 392
Solution:
1 1
MTTFsys 1909.9 hr End of Example
s s (s )
Solution:
2
(t) = e−λt (1 + PS × λt)
Rsys
393 Reliability Engineering
2
(1000) = e−0.5 (1 + 0.97 × 500 × 10−6 × 1000)=0.9007
Rsys
2 λ .
(1000) = e (1 + λt) = e (1 + 0.5)=0..9098
R
2 e s t
1)R (t ) e t 1 1 t 0
sys s
2
2) R (t ) e t 1 (1 p )t e t (1 Ps t ) t 0
sys
t e t 1 (1 p )t e t (1 P t )
e t 1 1 e s s
s
t
1e s
Ps
s t
function ( e s t )of constant- failure- rate switch over period(0 t):
t
s t
e
0
dt t
1 e s
switch average reliability Ps (7-46)
t 0 s t
R (t ) e t
1 PS (1 e t ) (7-47)
sys
t t
t t
MTTFsys Rsys (t )dt e 1 (1 p ) (1 e ) dt
0 0
Ps 1 p
1 1
1 1
MTTFsys (7-48)
Where
Example 7-8
2 t
R (t ) e t 1 1 e s t0
sys
s
Chap 7 Dynamic Models+Availability, Markov Chain 396
2
R (t ) e t
1 Ps (1 e t ) .
sys
a)Specify a condition under which these two expressions
identical if .
Answer :a)
t
1 e s
Ps
s 1 e t
1 1 2
f x f a x a f a x a f a ...
1! 2!
Therefore e
t is expanded as follows:
397 Reliability Engineering
×
=1+ (− )+ ⋯
!
1 e t t (7-49)
t t t s t
1e s 1e s e dt
0
Ps
s t s t t
That is under the condition that the value given for the switch
reliability( Ps ) equals the mean of reliability function of
exponential- distributed –lifetime switch over period(0 t) the 2
expression give the same results.
Answer: b)
Since the specified equals in Example 7-7, therefore if
t 1000
1e s or 1 e s or λs ≅ 61 ∗ 10−6 the 2 relationship
Ps 0.97
st s 1000
2
R (t ) e t
1 (1 p ) (1 e t )
sys
Chap 7 Dynamic Models+Availability, Markov Chain 398
2 t
R e t 1 1 e s t 0
sys
s
2
R (t ) e t 1 (1 p )t e t (1 Ps t ) t 0
sys
Give the same results if the failure rate of the redundant unit
while in standby mode equals the failure rate of the active unit
( ) and the given reliability value for the switch(Ps) equals
the value obtained from Eq. 7-46.
Example 7-9( Lewis,1994 page 262)
An engineer designs a standby system with two identical
units to have an idealized MTTF of 1000 days. To be
conservative, she then assumes a switching failure probability of
10%( p=0.10 or PS =0.9 ) and the failure rate of the unit in
Solution
active. If the standby unit were cold, according to Eq. 7-27-4 for the
idealized case:
399 Reliability Engineering
λ + = 0 .1 0 λ :
1 p 1 0.1
1 1
1 1 0.1 909days
MTTFsys
0.002
End of Example
The relationships for 3-component standby systems could be
studied in K&L pages 221-222. The interested readers in the
reliability function and the MTTF of the general case n=n, could
refer to Niaki & Yaghoubi(2020).
1
K&L page 222, Lewis(11996)page 260&chapter11
Chap 7 Dynamic Models+Availability, Markov Chain 400
Where
-λ f t 2 g t 2g 2 g t
R f (t)=e R sys (t ) e (e f t e ) (7-51)
2g f
Where
Example 7-10
λ = 0.001 ℎ , λ = 4 × 10 ℎ
According to Eq.7-51:
2 g t 2g 2 g t
R sys (t ) e (e f t e )
2g f
R (1000) =
× . ×
2 × 0.001 . × × . ×
e + (e −e )
2 × 0.001 − 0.004
Using MATLAB
=exp(-2*0.001*1000)+2*0.001*(exp(-0.004*1000)-exp(-2*0.001*1000))/(2*.001-0.004)
Serviceability
Maintainability
Maintainability function
t
M (t ) m (x )dx (7-52-2)
0
Solution
t
M (t ) 1 e t 1 e
End of Example
(MTTR )
i i n
MTTR i 1
i
(MTTRi ) (7-53 )
i 1
i i
where
i Failure rate of the ith repairable component
MTTR i Mean time to repair ith repairable unit
n number components in the system
i A fraction of failures per unit time related to
i
ith unit
It might be useful for some readers to know that some
references such as the manual of MIL-HDBK-472 standard
deal with MTTR in details. This manual is comprehensive
design tool for maintainability prediction analysis including
calculating MTTR.
Operational readiness(OR)
Solution
This mean that the machine is ready to perform its function 70%
of the time. End of Example
7-3 Availability
operating time
A (7-54-1)
operating time down time
Example 7-12:
The down time in one complete cycle of a machine is 6 time
units and its operating time is 8 units(See Fig 7-5). Find the
machine availability (A).
Solution
operating time 8
A A =0.57 57%
operating time down time 6+8
End of Example
Solution
A =
. .
a.r.t.= total down time- administrative& logistic times=6-2=4.
8
A = = 0.69 = 66%
8+4
Examples 7-13&14 show that by eliminating the administrative
and logistics time in the repair cycle, the current availability of
0.57 can be increased in the limit to the intrinsic availability of
0.66. There is a potential for a 9% improvement in availability.
End of Example
( )t
A (t ) e (7-56)
Chap 7 Dynamic Models+Availability, Markov Chain 410
λ
or = (7-58-2)
μ λ
process (see page 287 of Ross,1985) that A(t) in the long range
approaches to = .
μ
Example 7-15
It is desired that a machine which has an exponentially
distributed lifetime with mean 3000 hours to possess a steady-
413 Reliability Engineering
MTBF 3000
A 0.9995 MTTR 1.6 hr
MTBF MTTR 3000+MTTR
Example 7-16
A series system has 2 subsystems. One the is a compressor
with 80.37 failures per 10 working hours. The average of its
actual repair time is 89.3 hr. The other subsystem failure rate is
4.78 failures per 10 working hours and the actual repair time
on the average is 890.3. Calculate 26280-hr reliability of the
each subsystem , their steady state intrinsic availability
Solution
Constant failure rate λ implies that the life time distribution
1
is exponential with mean θ = and reliability function .
λ
therefore
For compressor:
λ =80.37 × 10−6 ⟹
θ MTBF1=1/80.37*10^-6=12442.4
( )=exp(-(t/ ))
(26280)=exp (-26280*80.37*(10^-6))=0.1210
.
= = = 0.9929
. . . . .
1 T
A * (T ) A (t )dt (7-60)
T 0
1 T
A * (T ) R (t )dt (7-61)
T 0
Example 7-17
The lifetime of an non-repairable switch with known MTTF
is exponentially distributed with parameter λ .
a)Calculate the parametric average of the reliability function of
the switch.
b)(Lewis ,1994 p301)
The system mission availability must be 0.95. Find the
maximum design life that can be tolerated in terms of the
MTTF.
Solution
a)For this switch which is non-repairable, A * i.e. the mission
availability equals the average of the reliability function(Ps) .
1 T
Ps A * (T ) R (t )dt TTF ~ exp( ) R (t ) e s t
T 0
1 T 1 e sT
Ps e s t dt
T 0 sT
Since the Taylor expansion of f(x) about a is;
1 1 2
f x f a x a f a x a f a ...
1! 2!
T
Then the expansion of e about a=0 :
- T T-0 (T-0) 2
e s =e-λs (0) + (-λ s )e-λs T |T=0 + [-λ s (-λ s )e-λsT ]|T=0 +...
1! 2!
1 1
e- sT 1-λ s T+ (λ s T)2 + (λ s T)3 +....
2 6
1 -λ T
Then for sT 1 , approximately e s 1-λsT+ (λsT) and:
2
2
417 Reliability Engineering
1
1-1+λ s T- (λ s T) 2
1-e-λs T 2 1
Ps = 1- λ s T A * (T )
λsT λs T 2
∗(
b) ) = 0.95 then
1 1
0.95 1- λsT T T 0.1 MTTF
2 10
s
End of Example
= ( ,…, ) (7-63)
where
Example 7-18
The lifetime and downtimes of n independent components,
= ( ,…, )= ×. . .×
= 2
( 1, … , ) = 1 − (1 − 1) ×. . .× (1 − )
= ( ,…, ) =1 − (1 − ). . . (1 − )
Example 7-19
A 2-unit system fails when either of its units fail. The units
have the steady state availability of 0.9958 and 0.9929.
Calculate the steady state availability of the system.
Solution
= ( 1, 2) = 1 × 2
= ( , )= × A=0.9929*0.9958=0.9887
availability and costs. In the steady state, i.e. after any transient
behavior has settled down and assuming that maintenance
actions occur at a constant rate" (O'connor,2003 page 402):
C
= × (7-64 )
+ C+T
where
C= Preventive maintenance cycle [ e.g. every 1000 hr]
T=Total time required to perform preventive maintenance tasks
Fig 7-8
8 shows a diagram of the state space of the critical
system.. The possible states are denoted by symbols such as
(Computer A is active) , (Computer B is in standby mode)
Fig 7-8
7 The 9 possible states of a two-unit
unit system
(Barlow and Proschan , 1996p120)
passes around the perimeter of the square…. States 4,6 & 8 are
unfavorable. (Barlow and Proschan , 1996page121).
The user of this system may be interested in such information
as the mean system down time during a specified time interval,
the probability that the system is down more than x minutes at
any one time. Under certain reasonable assumptions on time to
failure(TTF) , the time to perform repair(TTR), etc. The
operation of the system can be described by a semi-Markov
process to get the desired information. Chapter 5 of Barlow&
Proschan(1996) deal with this system in detail.
As another example consider a system having 3 components
or units a, b, c. To use a Markov chain the states of this system
are defined as combinations of operating and failed components.
As the following table shows the system, depending on the
operation or failure of the components(o=operating X= failed),
has 8 states(Lewis,1994page 326):
unit State
1 2 3 4 5 6 7 8
a o X o o X X o X
b o o X o X o X X
c o o o X o X X X
o=operating X= failed
Exercises
system reliability.
Chapter 8
Enhancement,
Optimization
&Allocation
of Reliability
Chap 8 Enhancement, Optimization & Allocation of Reliability 426
8
Enhancement, Optimization & Allocation of
Reliability
Example 8.1
Solution
a) R=0.8 =0.512.
1
https://reliabilityanalyticstoolkit.appspot.com/standby_redundancy_integrate
Chap 8 Enhancement, Optimization & Allocation of Reliability 430
/ Max ( )
. . ( )≤0 = 1,2, … ,
ℎ( )=0 = 1,2, … ,
≥0
Fig.
Fig.8-5 A series-parallel system(Faghih,1996 p102)
k
Total number of components in the system is n i
. If we
i 1
k
R sys 1 (1 R i ) ni (8-1)
i 1
where
k number of subsystems
Ri The reliability of each component in Subsystem i
ni Number of components in Subsystem i
Now let
k = Maximum budget available
Ci = The cost of each component in Subsystem i
Then :
Chap 8 Enhancement, Optimization & Allocation of Reliability 432
n C
i 1
i i C ; (8-2)
k
M ax R sys f ( n 1 ,..., n k ) 1 (1 R i ) n i
i 1
s .t .
k
n C
i 1
i i C
0 R i 1 i 1,..., k
n 0 n i =Integer
i
k
n
M ax R 1 (1 R ) i
i 1
i
s .t .
k
a ij ni b j i 1, 2,..., k j 1, 2,.., m
i 1
0 Ri 1
n 0 n i =Integer
i
where
Number of materials
ni Number of components in ith subsystem
k
Min Z ni C i
i 1
s.t.
k
R sys 1 (1 R i )ni R 0
i 1
0 R i 1 i 1,..., k
n 0
i
n = Integer
i
1
The refrence of this chapter is mainly K&L, Chap 14.
Chap 8 Enhancement, Optimization & Allocation of Reliability 436
Where
…..e ≥ (8-4)
or
+ λ +. . . . +λ ≤ (8-5)
Theoretically, the above equation has an infinite number of
solutions, assuming no restrictions on the allocation. The
problem is to establish a procedure that yields a unique or
limited number of solutions by which consistent and reasonable
reliabilities may be allocated1. Some of these procedures are:
1
From: http://reliabilityanalytics.com/blog/2011/10/09/reliability-allocation/
Chap 8 Enhancement, Optimization & Allocation of Reliability 438
8. Dynamic Programming
= × × …× ⟹ =
= = i = 1, … , (8-5-1)
=R i = 1, … , n ⟹
= 1 − (1 − R) = 1- ⟹1-R= (1 − )⟹
=1− (1 − ) (8-5-2)
439 Reliability Engineering
n
R = (1 − R) (R)
x
(∑ (1 − R) (R) − R )=0.
Example 8-2
The reliability requirement for a 3-component series system
is = 0.8573. Find the reliability of each component using
equal apportionment technique.
Solution
= √0.8573 = 0.95
Example 8-3
Solution
Example 8-4
Solution
= (1 − ) ( )
then:
5
(1 − ) ( ) − 0.99144 = 0
∗
The ARINC method tries to choose such that(K&L p 407):
∗ ∗
+. . . . + ≤ 0
i 1
i
∗
= (8-7)
Solution
0.005
= = 0.555, w = 0.333 ,
0.005 + 0.003 + 0.001
w = 0.111
To find the required failure rate() for this exponentially-
distributed-lifetime system we could write:
∗
= (0.555)(0.00256) = 0.00142
∗
= (0.333)(0.00256) = 0.000582
∗
= (0.111)(0.00256) = 0.000284
∗ (20) ∗ ( . )
=e =e = 0.97
∗ (20) ∗
=e = 0.98 , R∗ (20) = 0.99
which are taken out of the system before the end of the mission
time considered for the system. Notice that the total system
under consideration is not truly a series system unless all of the
importance indices (wi's) equal unity and the mission time of all
subsystems ((ti's) are equal(Grosh,1989 p150).
To reach a target MTBF for the system, this method uses
Equation 8-8 which calculates an approximate value for the
MTBF of each subsystem(Grosh,1989 p150).
( )( )( )
MTBF = [ ( )]
= 1,2, … (8-8)
This is equivalent to
the following failure rate for the subsystem(K&L p409):
[ ∗( )]
= = 1,2, … (8-9-1)
or
the following reliability for the subsystem
∗
( )= (8-9-2)
where
Number of components of in the ith subsystem
N Total number of components in the system: N= ∑
failure rate of ith subsystem
( ) The required system reliability for a mission time t
∗
The reliability allocated to ith subsystem
System mission time
The mission time for ith subsystem; the time period
required for the ith subsystem to operate from the
445 Reliability Engineering
or
∗ [ ( )]
= 1− (8-10-2)
1 10 1 15
2 9 0.95 25
3 10 1 100
4 8 0.9 70
sum N=210
Solution
=
[− ( )]
210 ∗ 1 ∗ 10
= = 2729.4
(15)(− 0.95)
1 1
= = = 0.00036638 = 36638 × 10
2729.4
210 ∗ 0.95 ∗ 9
= = 1400.2 = 0.0007142
(25)(− 0.95)
210 ∗ 1 ∗ 10
= = 409.41 = 0.002442
100(− 0.95)
447 Reliability Engineering
210 ∗ 0.9 ∗ 8
= = 421.1 = 0.002374
(70)(− 0.95)
∗ 1−[ ( )]
=1−
∗
1 − (0.95)
=1− = 0.99634
1
∗ ∗ ( . )( )
or = e =e = 0.99634
similarly
∗
1 − (0.95)
=1− = 0.99359
0.95
∗ ∗ ( . )( )
or = e =e = 0.99359
∗ (t)
= 0.975870,
∗ (t)
= 0.98116
∗ ∗ ∗ ∗
× × × = 0.94788=94.79 % for the system which is
slightly less than the system reliability requirement 0.95. This is
a result of the approximate nature of the AGREE formula and
that and are less than unity.End of Example
Exercises
Number of Operating
Subsystem No.( ) subsystems time
ni wi ti
1 25 1.00 10
2 80 0.97 9
3 45 1.00 10
4 60 0.93 7
5 70 1.00 10
449 Reliability Engineering
References
American National Standard ,1989
Bazargan, Hamid,
Statistical methods in Quality control
Downloadable from:
https://opentextbc.ca/oerdiscipline/chapter/industrial-engineering or
https://archive.org/details/statistical-methods-august-2020-bazargan or
https://opentextbc.ca/oerdiscipline/chapter/statistics/
Bazovsky, Igor, 2004
Reliability Theory and Practice
Dover Publication Inc
Barlow,E.B.& Proschan F,L1996
Math. Theory of Reliability,
S.I.A.M
Billinton, R, Allan, R.,1992
Reliability Evaluation of Engineering Systems: Concepts and Techniques
Plenum Publishers
Bowker, A., H. and. Leiberman, G. J., 1972
Engineering Statistics
Prentice Hall, Inc., Englewood Cliffs, N.J
Cabarbaye,A.2019
Implementation of accelerated life testing
Cab Innovation, France(www.cabinovation.com)
Carter,A.D.S., 1986
Mechanical Reliability
Macmillan
Chen, Qiming, 2004
"The probabibity, identification, and prevention of rare events in
power systems "
Iowa State univ. Retrospective Theses andDissertations.Paper 1149.
Coles , Stuart, 2001
An Introduction to Statistical Modeling of Extreme Values
Springer Verlag, London.
Dhillon, B.S.2006
Maintainability, Maintenance, and Reliability for Engineers
CRC Press
Taylor & Francis Group
Dao-Thein, M.&Massoud,M.1974
On the Relation Between the Factor of Safety and Reliability
Jr of Engineering for Industry , 96(3): 853–857.
doi: https://doi.org/10.1115/1.3438452
https://citeseerx.ist.psu.edu/viewdoc/download?doi=
10.1.1.1058.9057&rep=rep1&type=pdf
Ebeling, Charles E.1997
Introduction to Reliability and Maintainability Engineering
McGraw-Hill
451 Reliability Engineering
T A B L E S
Table A Crtical values of F distribution Fm,n ( Mood et al,1974 ) Example . , , = .
1 2 3 4 5 6 7 8 9 10 12 15 20 30 60 120
.10 39.9 49.5 53.6 55.8 57.2 58.2 58.9 59.4 59.9 60.2 60.7 61.2 61.7 62.3 62.8 63.1 63.3
.05 161 200 216 225 230 234 237 239 241 242 244 246 248 250 252 253 254
.02 5 648 800 864 900 922 937 948 957 963 969 977 985 993 1000 1010 1010 1020
.01 1 4050 5000 5400 5620 5760 5860 5930 5980 6020 6060 6110 6160 6210 6260 6310 6340 6370
.00 5 16200 20000 21600 22500 23100 23400 23700 23900 24100 24200 24400 24600 24800 25000 25200 25400 25500
.10 8.53 9.00 9.16 9.24 9.29 9.33 9.35 9.37 9.38 9.39 9.41 9.42 9.44 9.46 9.47 9.48 9.49
.05 18.5 19.0 19.2 19.2 19.3 19.3 19.4 19.4 19.4 19.4 19.4 19.4 19.5 19.5 19.5 19.5 19.5
.02 5 38.5 39.0 39.2 39.2 39.3 39.3 39.4 39.4 39.4 39.4 39.4 39.4 39.4 39.5 39.5 39.5 39.5
2
.01 98.5 99.0 99.2 99.2 99.3 99.3 99.4 99.4 99.4 99.4 99.4 99.4 99.4 99.5 99.5 99.5 99.5
.00 5 199 199 199 199 199 199 199 199 199 199 199 199 199 199 199 199 199
.10 5.54 5.46 5.39 5.34 5.31 5.28 5.27 5.25 5.24 5.23 5.22 5.20 5.18 5.17 5.15 5.14 5.13
.05 10.1 9.55 9.28 9.12 9.01 8.94 8.89 8.85 8.81 8.79 8.74 8.70 8.66 8.62 8.57 8.55 8.53
.02 5 17.4 16.0 15.4 15.1 14.9 14.7 14.6 14.5 14.5 14.4 14.3 14.3 14.2 14.1 14.0 13.9 13.9
3
.01 34 .1 30.8 29.5 28.7 28.2 27.9 27.7 27.5 27.3 27.2 27.1 26.9 26.7 26.5 26.3 26.2 26.1
.00 5 55.6 49.8 47.5 46.2 45.4 44.8 44.4 44.1 43.9 43.7 43.4 43.1 42.8 42.5 42.1 42.0 41.8
.10 4.54 4.32 4.19 4.11 4.05 4.01 3.98 3.95 3.93 3.92 3.90 3.87 3.84 3.82 3.79 3.78 3.76
.05 7.71 6.94 6.59 6.39 6.26 6.16 6.09 6.04 6.00 5.96 5.91 5.86 5.80 5.75 5.69 5.66 5.63
.02 5 12.2 10.6 9.98 9.60 9.36 9.20 9.07 8.98 8.90 8.84 8.75 8.66 8.56 8.46 8.3& 8.31 8.26
.01 4 21.2 18.0 16.7 16.0 15.5 15.2 15.0 14.8 14.7 14.5 14.4 14.2 14.0 13.8 13.7 13.6 13.5
.00 5 31.3 26.3 24.3 23.2 22.5 22.0 21.6 21.4 21.1 21.0 20.7 20.4 20.2 19.9 19.6 19.5 19.3
457 Reliability Engineering
Table A -continued
.10 4.06 3.78 3.62 3.52 3.45 3.40 3.37 3.34 3.32 3.30 3.27 3.24 3.21 3.17 3.14 3.12 3.11
.05 6.61 5.79 5.41 5.19 5.05 4.95 4.88 4.82 4.77 4.74 4.68 4.62 4.56 4.50 4.43 4.40 4.37
.02 5 10.0 8.43 7.76 7.39 7.15 6.98 6.85 6.76 6.68 6.62 6.52 6.43 6.33 6.23 6.12 6.07 6.02
5
.01 16.3 13.3 12.1 11.4 11.0 10.7 10.5 10.3 10.2 10.1 9.89 9.72 9.55 9.38 9.20 9.11 9.02
.00 5 22.8 18.3 16.5 15.6 14.9 14.5 14.2 14.0 13.8 13.6 13.4 13.1 12.9 12.7 12.4 12.3 12.1
.10 3.78 3.46 3.29 3.18 3.11 3.05 3.01 2.98 2.96 2.94 2.90 2.87 2.84 2.80 2.76 2.74 2.72
.05 5.99 5.14 4.76 4.53 4.39 4.28 4.21 4.15 4.10 4.06 4.00 3.94 3.87 3.81 3.74 3.70 3.67
.02 5 8.81 7.26 6.60 6.23 5.99 5.82 5.70 5.60 5.52 5.46 5.37 5.27 5.17 5.07 4.96 4.90 4.85
6
.01 13.7 10.9 9.78 9.15 8.75 8.47 8.26 8.10 7.98 7.87 7.72 7.56 7.40 7.23 7.06 6.97 6.88
.00 5 18.6 14.5 12.9 12.0 11.5 11.1 10.8 10.6 10.4 10.2 10.0 9.81 9.59 9.36 9.12 9.00 8.88
.10 3.59 3.26 3.07 2.96 2.88 2.83 2.78 2.75 2.72 2.70 2.67 2 . 63 2.59 2.56 2.51 2.49 2.47
.05 5.59 4.74 4.35 4.12 3.97 3.87 3.79 3.73 3.68 3.64 3.57 3.51 3.44 3.38 3.30 3.27 3.23
.02 5 8.07 6.54 5.89 5.52 5.29 5.12 4.99 4.90 4.82 4.76 4.67 4.57 4.47 4.36 4.25 4.20 4.14
.01 7 12.2 9.55 8.45 7.85 7.46 7.19 6.99 6.84 6.72 6.62 6.47 6.31 6.16 5.99 5.82 5.74 5.65
.00 5 16.2 12.4 10.9 10.1 9.52 9.16 8.89 8.68 8.51 8.38 8.18 7.97 7.75 7.53 7.31 7.19 7.08
.10 3.46 3.11 2.92 2.81 2.73 2.67 2.62 2.59 2.56 2.54 2.50 2.46 2.42 2.38 2.34 2.31 2.29
.05 5.32 4.46 4.07 3.84 3.69 3.58 3.50 3.44 3.39 3.35 3.28 3.22 3.15 3.08 3.01 2.97 2.93
.02 5 7.57 6.06 5.42 5.05 4.82 4.65 4.53 4.43 4.36 4.30 4.20 4.10 4.00 3.89 3.78 3.73 3.67
.01 8 11.3 8.65 7.59 7.01 6.63 6.37 6.18 6.03 5.91 5.81 5.67 5.52 5.36 5.20 5.03 4.95 4.86
.00 5 14.7 11.0 9.60 8.81 8.30 7.95 7.69 7.50 7.34 7.21 7.01 6.81 6.61 6.40 6.18 6.06 5.95
458
1 2 3 4 5 6 7 8 9 10 12 15 20 30 60 120 ∞
.10 3.36 3.01 2.81 2.69 2.61 2.55 2.51 2.47 2.44 2.42 2.38 2.34 2.30 2.25 2.21 2.18 2.16
.05 5.12 4.26 3.86 3.63 3.48 3.37 3.29 3.23 3.18 3.14 3.07 3.01 2.94 2.86 2.79 2.75 2.71
.02 5 7.21 5.71 5.08 4.72 4.48 4.32 4.20 4.10 4.03 3.96 3.87 3.77 3.67 3.56 3.45 3.39 3.33
.01 9 10 .6 8.02 6.99 6.42 6.06 5.80 5.61 5.47 5.35 5.26 5.11 4.96 4.81 4.65 4.48 4.4O 4.31
.005 13.6 10.1 8.72 7.96 7.47 7.13 6.88 6.69 6.54 6.42 6.23 6.03 5.83 5.62 5.41 5.30 5.19
.10 3.29 2.92 2.73 2.61 2.52 2.46 2.41 2.38 2.35 2.32 2.28 2.24 2.20 2.15 2.11 2.08 2.06
.05 4.96 4.10 3.71 3.48 3.33 3.22 3.14 3.07 3.02 2.98 2.91 2.84 2.77 2.70 2.62 2.58 2.54
.025 6.94 5.46 4.83 4.47 4.24 4.07 3.95 3.85 3.78 3.72 3.62 3.52 3.42 3.31 3.20 3.14 3.08
.01 10 10.0 7.56 6.55 5.99 5.64 5.39 5.20 5.06 4.94 4.85 4.71 4.56 4.41 4.25 4.08 4.00 3.91
.005 12.8 9.43 8.08 7.34 6.87 6.54 6.30 6.12 5.97 5.85 5.66 5.47 5.27 5.07 4.86 4.75 4.64
.10 3.18 2 81 2.61 2.48 2 39 2.33 2.28 2.24 2.21 2.19 2.15 2.10 2.06 2 .01 1.96 1.93 1.90
.05 4.75 3.89 3.49 3.26 3.11 3.00 2.91 2.85 2.80 2.75 2.69 2.62 2.54 2.47 2.38 2.34 2.30
.025 6.55 5.10 4.47 4.12 3.89 3.73 3.61 3.51 3.44 3.37 3.28 3.18 3.07 2.96 2.85 2.79 2.72
.01 12 9.33 6.93 5.95 5.41 5.06 4.82 4.64 4.50 4.39 4.30 4.16 4.01 3.86 3.70 3.54 3.45 3.36
.005 11.8 8.51 7.23 6.52 6.07 5.76 5.52 5.35 5.20 5.09 4.91 4.72 4.53 4.33 4.12 4.01 3.90
.10 3.07 2.70 2.49 2.36 2.27 2.21 2:16 2.12 2.09 2.06 2.02 1.97 1.92 1.87 1.82 1.79 1.76
.05 4.54 3.68 3.29 3.06 2.90 2.79 2.71 2.64 2.59 2.54 2.48 2.40 2.33 2.25 2.16 2.11 2.07
.025 6.20 4.77 4.15 3.80 3.58 3.41 3.29 3.20 3.12 3.06 2.96 2.86 2.76 2.64 2.52 2.46 2.40
.01 15 8.68 6.36 5.42 4.89 4.56 4.32 4.14 4.00 3.89 3.80 3.67 3.52 3.37 3.21 3.05 2.96 2.87
.005 10.8 7.70 6.48 5.80 5.37 5.07 4.85 4.67 4.54 4.42 4.25 4.07 3.88 3.69 3.48 3.37 3.26
459 Reliability Engineering
Table A -continued
.10 2.97 2.59 2.38 2.25 2.16 2.09 2.04 2.00 1.96 1.94 1.89 1.84 1.79 1.74 1.68 1.64 1.61
.05 4.35 3.49 3.10 2.87 2.71 2.60 2.51 2.45 2.39 2.35 2.28 2.20 2.12 2.04 1.95 1.90 1.84
.025 5.87 4.46 3.86 3.51 3.29 3.13 3.01 2.91 2.84 2.77 2.68 2.57 2.46 2.35 2.22 2.16 2.09
.01 20 8.10 5.85 4.94 4.43 4.10 3.87 3.70 3.56 3.46 3.37 3.23 3.09 2.94 2.78 2.61 2.52 2.42
.005 9.94 6.99 5.82 5.17 4.76 4.47 4.26 4.09 3.96 3.85 3.68 3.50 3.32 3.12 2.92 2.81 2.69
.10 2.88 2.49 2.28 2.14 2.05 1.98 1.93 1.88 1.85 1.82 1.77 1.72 1.67 1.61 1.54 1.50 1.46
.05 4.17 3.32 2.92 2.69 2.53 2.42 2.33 2.27 2.21 2.16 2.09 2.01 1.93 1.84 1.74 1.68 1.62
.025 5.57 4.18 3.59 3.25 3.03 2.87 2.75 2.65 2.57 2.51 2.41 2.31 2.20 2.07 1.94 1.87 1.79
.01 30 7.56 5.39 4.51 4.02 3.70 3.47 3.30 3.17 3.07 2.98 2.84 2.70 2.55 2.39 2.21 2.11 2.01
.005 9.18 6.35 5.24 4.62 4.23 3.95 3.74 3.58 3.45 3.34 3.18 3.01 2.82 2.63 2.42 2.30 2.18
.10 2.79 2.39 2.18 2.04 1.95 1.87 1.82 1.77 1.74 1.71 1.66 1.60 1.54 1.48 1.40 1.35 1.29
.05 4.00 3.15 2.76 2.53 2.37 2.25 2.17 2.10 2.04 1.99 1.92 1.84 1.75 1.65 1.53 1.47 1.39
.025 5.29 3.93 3.34 3.01 2.79 2.63 2.51 2.41 2.33 2.27 2.17 2.06 1.94 1.82 1.67 1.58 1.48
.01 60 7.08 4.98 4.13 3.65 3.34 3.12 2.95 2.82 2.72 2.63 2.50 2.35 2.20 2.03 1.84 1.73 1.60
.005 8.49 5.80 4.73 4.14 3.76 3.49 3.29 3.13 3.01 2.90 2.74 2.57 2.39 2.19 1.96 1.83 1.69
.10 2.75 2.35 2.13 1.99 1.90 1.82 1.77 1.72 1.68 1.65 1.60 1.54 1.48 1.41 1.32 1.26 1.19
.05 3.92 3.07 2.68 2.45 2.29 2.18 2.09 2.02 1.96 1.91 1.83 1.75 1.66 1.55 1.43 1.35 1.25
.025 5.15 3.80 3.23 2.89 2.67 2.52 2.39 2.30 2.22 2.16 2.05 1.94 1.82 1.69 1.63 1.43 1.31
.01 120 6.85 4.79 3.95 3.48 3.17 2.96 2.79 2.66 2.56 2.47 2.34 2.19 2.03 1.86 1.66 1.53 1.38
.005 8.18 5.54 4.50 3.92 3.55 3.28 3.09 2.93 2.81 2.71 2.54 2.37 2.19 1.98 1.75 1.61 1.43
.10 2.71 2.30 2.08 1.94 1.85 1.77 1.72 1.67 1.63 1.60 1.55 1.49 1.42 1.34 1.24 1.17 1.00
.05 3.84 3.00 2. 60 2.37 2. 21 2.10 2.01 1.94 1.88 1.83 1.75 1.67 1.57 1.46 1.32 1.22 1.00
.025 5.02 3.69 3.12 2.79 2.57 2.41 2.29 2.19 2.11 2.05 1.94 1.83 1.71 1.57 1.39 1.27 1.00
.01 6.63 4.61 3.78 3.32 3.02 2.80 2.64 2.51 2.41 2.32 2.18 2.04 1.88 1.70 1.47 1.32 1.00
.005 7.88 5.30 4.28 3.72 3.35 3.09 2.90 2.74 2.62 2.52 2.36 2.19 2.00 1.79 1.53 1.36 1.00
460
λ Table B Some values of CDF of Poisson Distribution Pr(X k )
or k
np 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
0.01 0.990 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.10 0.905 0.995 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.20 0.819 0.982 0.999 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.30 0.741 0.963 0.996 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.40 0.670 0.938 0.992 0.999 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
0.80 0.449 0.809 0.953 0.991 0.999 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
1.00 0.368 0.736 0.920 0.981 0.996 0.999 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000 1.000
2.00 0.135 0.406 0.677 0.857 0.947 0.983 0.995 0.999 1.000 1.000 1.000 1.000 1.000 1.000 1.000
3.00 0.050 0.199 0.423 0.647 0.815 0.916 0.966 0.988 0.996 0.999 1.000 1.000 1.000 1.000 1.000
4.00 0.018 0.092 0.238 0.433 0.629 0.785 0.889 0.949 0.979 0.992 0.997 0.999 1.000 1.000 1.000
5.00 0.007 0.040 0.125 0.265 0.440 0.616 0.762 0.867 0.932 0.968 0.990 0.995 0.998 0.999 1.000
6.00 0.002 0.017 0.062 0.151 0.285 0.446 0.606 0.744 0.847 0.916 0.960 0.980 0.991 0.996 0.999
6.20 0.002 0.015 0.054 0.134 0.259 0.414 0.574 0.716 0.826 0.902 0.950 0.975 0.989 0.995 0.998
6.40 0.002 0.012 0.046 0.119 0.235 0.384 0.542 0.687 0.803 0.886 0.940 0.969 0.986 0.994 0.997
6.60 0.001 0.010 0.040 0.105 0.213 0.355 0.511 0.658 0.780 0.869 0.930 0.963 0.982 0.992 0.997
461 Reliability Engineering
λ Table B Some values of CDF of Poisson Distribution Pr(X k )
or k
np 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
7.00 0.001 0.007 0.030 0.082 0.173 0.301 0.450 0.599 0.729 0.830 0.900 0.947 0.973 0.987 0.994
8.00 0.000 0.003 0.014 0.042 0.100 0.191 0.313 0.453 0.593 0.717 0.820 0.888 0.936 0.966 0.983
9.00 0.000 0.001 0.006 0.021 0.055 0.116 0.207 0.324 0.456 0.587 0.710 0.803 0.876 0.926 0.959
10.00 0.000 0.000 0.003 0.010 0.029 0.067 0.130 0.220 0.333 0.458 0.580 0.697 0.792 0.864 0.917
11.00 0.000 0.000 0.001 0.005 0.015 0.038 0.079 0.143 0.232 0.341 0.460 0.579 0.689 0.781 0.854
12.00 0.000 0.000 0.001 0.002 0.008 0.020 0.046 0.090 0.155 0.242 0.350 0.462 0.576 0.682 0.772
12.50 0.000 0.000 0.000 0.002 0.005 0.015 0.035 0.070 0.125 0.201 0.300 0.406 0.519 0.628 0.725
13.00 0.000 0.000 0.000 0.001 0.004 0.011 0.026 0.054 0.100 0.166 0.250 0.353 0.463 0.573 0.675
13.50 0.000 0.000 0.000 0.001 0.003 0.008 0.019 0.041 0.079 0.135 0.210 0.304 0.409 0.518 0.623
14.00 0.000 0.000 0.000 0.000 0.002 0.006 0.014 0.032 0.062 0.109 0.180 0.260 0.358 0.464 0.570
14.50 0.000 0.000 0.000 0.000 0.001 0.004 0.010 0.024 0.048 0.088 0.140 0.220 0.311 0.413 0.518
15.00 0.000 0.000 0.000 0.000 0.001 0.003 0.008 0.018 0.037 0.070 0.120 0.185 0.268 0.363 0.466
15.50 0.000 0.000 0.000 0.000 0.001 0.002 0.006 0.013 0.029 0.055 0.100 0.154 0.228 0.317 0.415
16.00 0.000 0.000 0.000 0.000 0.000 0.001 0.004 0.010 0.022 0.043 0.080 0.127 0.193 0.275 0.368
16.50 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.007 0.017 0.034 0.060 0.104 0.162 0.236 0.323
17.00 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.005 0.013 0.026 0.050 0.085 0.135 0.201 0.281
17.50 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.004 0.009 0.020 0.040 0.068 0.112 0.170 0.243
18.00 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.007 0.015 0.030 0.055 0.092 0.143 0.208
462
λ Table B Some values of CDF of Poisson Distribution Pr(X k )
or k
np 0 1 2 3 4 5 6 7 8 9 10 11 12 13 14
18.50 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.005 0.012 0.020 0.044 0.075 0.119 0.177
19.00 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.004 0.009 0.020 0.035 0.061 0.098 0.150
19.50 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.007 0.010 0.027 0.049 0.081 0.126
20.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.005 0.010 0.021 0.039 0.066 0.105
20.50 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.004 0.010 0.017 0.031 0.054 0.087
21.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.010 0.013 0.025 0.043 0.072
21.50 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.005 0.010 0.019 0.035 0.059
22.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.004 0.008 0.015 0.028 0.048
22.50 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.006 0.012 0.022 0.039
23.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.004 0.009 0.017 0.031
23.50 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.003 0.007 0.014 0.025
24.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.003 0.005 0.011 0.020
24.50 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.001 0.002 0.004 0.008 0.016
463 Reliability Engineering
Table C Area under standard normal curve(Pr(Z < )) e.g. Pr( < −3.00) = 0.0013
x
z
0.09 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0.00
-3.5 0.00017 0.00017 0.00018 0.00019 0.00019 0.0002 0.00021 0.00022 0.00022 0.00023
-3.4 0.00024 0.00025 0.00026 0.00027 0.00028 0.00029 0.0003 0.00031 0.00032 0.00034
-3.3 0.00035 0.00036 0.00038 0.00039 0.0004 0.00042 0.00043 0.00045 0.00047 0.00048
-3.2 0.0005 0.00052 0.00054 0.00056 0.00058 0.0006 0.00062 0.00064 0.00066 0.00069
-3.1 0.00071 0.00074 0.00076 0.00079 0.00082 0.00084 0.00087 0.0009 0.00094 0.00097
-3 0.001 0.00104 0.00107 0.00111 0.00114 0.00118 0.00122 0.00126 0.00131 0.00135
-2.9 0.00139 0.00144 0.00149 0.00154 0.00159 0.00164 0.00169 0.00175 0.00181 0.00187
-2.8 0.00193 0.00199 0.00205 0.00212 0.00219 0.00226 0.00233 0.0024 0.00248 0.00256
-2.7 0.00264 0.00272 0.0028 0.00289 0.00298 0.00307 0.00317 0.00326 0.00336 0.00347
-2.6 0.00357 0.00368 0.00379 0.00391 0.00402 0.00415 0.00427 0.0044 0.00453 0.00466
-2.5 0.0048 0.00494 0.00508 0.00523 0.00539 0.00554 0.0057 0.00587 0.00604 0.00621
-2.4 0.00639 0.00657 0.00676 0.00695 0.00714 0.00734 0.00755 0.00776 0.00798 0.0082
-2.3 0.00842 0.00866 0.00889 0.00914 0.00939 0.00964 0.0099 0.01017 0.01044 0.01072
-2.2 0.01101 0.01130 0.0116 0.01191 0.01222 0.01255 0.01287 0.01321 0.01355 0.01390
-2.1 0.01426 0.01463 0.015 0.01539 0.01578 0.01618 0.01659 0.01700 0.01743 0.01786
-2 0.01831 0.01876 0.01923 0.0197 0.02018 0.02068 0.02118 0.02169 0.02222 0.02275
-1.9 0.0233 0.02385 0.02442 0.025 0.02559 0.02619 0.0268 0.02743 0.02807 0.02872
-1.8 0.02938 0.03005 0.03074 0.03144 0.03216 0.03288 0.03362 0.03438 0.03515 0.03593
-1.7 0.03673 0.03754 0.03836 0.0392 0.04006 0.04093 0.04182 0.04272 0.04363 0.04457
464
Z 0.09 0.08 0.07 0.06 0.05 0.04 0.03 0.02 0.01 0.00
-1.6 0.04551 0.04648 0.04746 0.04846 0.04947 0.0505 0.05155 0.05262 0.0537 0.0548
-1.5 0.05592 0.05705 0.05821 0.05938 0.06057 0.06178 0.06301 0.06426 0.06552 0.06681
-1.4 0.06811 0.06944 0.07078 0.07215 0.07353 0.07493 0.07636 0.0778 0.07927 0.08076
-1.3 0.08226 0.08379 0.08534 0.08691 0.08851 0.09012 0.09176 0.09342 0.0951 0.0968
-1.2 0.09853 0.10027 0.10204 0.10383 0.10565 0.10749 0.10935 0.11123 0.11314 0.11507
-1.1 0.11702 0.119 0.121 0.12302 0.12507 0.12714 0.12924 0.13136 0.1335 0.13567
-1 0.13786 0.14007 0.14231 0.14457 0.14686 0.14917 0.15151 0.15386 0.15625 0.15866
-0.9 0.16109 0.16354 0.16602 0.16853 0.17106 0.17361 0.17619 0.17879 0.18141 0.18406
-0.8 0.18673 0.18943 0.19215 0.19489 0.19766 0.20045 0.20327 0.20611 0.20897 0.21186
-0.7 0.21476 0.2177 0.22065 0.22363 0.22663 0.22965 0.2327 0.23576 0.23885 0.24196
-0.6 0.2451 0.24825 0.25143 0.25463 0.25785 0.26109 0.26435 0.26763 0.27093 0.27425
-0.5 0.2776 0.28096 0.28434 0.28774 0.29116 0.2946 0.29806 0.30153 0.30503 0.30854
-0.4 0.31207 0.31561 0.31918 0.32276 0.32636 0.32997 0.3336 0.33724 0.3409 0.34458
-0.3 0.34827 0.35197 0.35569 0.35942 0.36317 0.36693 0.3707 0.37448 0.37828 0.38209
-0.2 0.38591 0.38974 0.39358 0.39743 0.40129 0.40517 0.40905 0.41294 0.41683 0.42074
-0.1 0.42465 0.42858 0.43251 0.43644 0.44038 0.44433 0.44828 0.45224 0.4562 0.46017
0 0.46414 0.46812 0.4721 0.47608 0.48006 0.48405 0.48803 0.49202 0.49601 0.5
465
Table C continued e.g. Pr( < 1.04) = 0.14917 Reliability Engineering
0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
0 0.5 0.50399 0.50798 0.51197 0.51595 0.51994 0.52392 0.5279 0.53188 0.53586
0.1 0.53983 0.5438 0.54776 0.55172 0.55567 0.55962 0.56356 0.56749 0.57142 0.57535
0.2 0.57926 0.58317 0.58706 0.59095 0.59483 0.59871 0.60257 0.60642 0.61026 0.61409
0.3 0.61791 0.62172 0.62552 0.6293 0.63307 0.63683 0.64058 0.64431 0.64803 0.65173
0.4 0.65542 0.6591 0.66276 0.6664 0.67003 0.67364 0.67724 0.68082 0.68439 0.68793
0.5 0.69146 0.69497 0.69847 0.70194 0.7054 0.70884 0.71226 0.71566 0.71904 0.7224
0.6 0.72575 0.72907 0.73237 0.73565 0.73891 0.74215 0.74537 0.74857 0.75175 0.7549
0.7 0.75804 0.76115 0.76424 0.7673 0.77035 0.77337 0.77637 0.77935 0.7823 0.78524
0.8 0.78814 0.79103 0.79389 0.79673 0.79955 0.80234 0.80511 0.80785 0.81057 0.81327
0.9 0.81594 0.81859 0.82121 0.82381 0.82639 0.82894 0.83147 0.83398 0.83646 0.83891
1 0.84134 0.84375 0.84614 0.84849 0.85083 0.85314 0.85543 0.85769 0.85993 0.86214
1.1 0.86433 0.8665 0.86864 0.87076 0.87286 0.87493 0.87698 0.879 0.881 0.88298
1.2 0.88493 0.88686 0.88877 0.89065 0.89251 0.89435 0.89617 0.89796 0.89973 0.90147
1.3 0.9032 0.9049 0.90658 0.90824 0.90988 0.91149 0.91309 0.91466 0.91621 0.91774
1.4 0.91924 0.92073 0.9222 0.92364 0.92507 0.92647 0.92785 0.92922 0.93056 0.93189
1.5 0.93319 0.93448 0.93574 0.93699 0.93822 0.93943 0.94062 0.94179 0.94295 0.94408
1.6 0.9452 0.9463 0.94738 0.94845 0.9495 0.95053 0.95154 0.95254 0.95352 0.95449
1.7 0.95543 0.95637 0.95728 0.95818 0.95907 0.95994 0.9608 0.96164 0.96246 0.96327
1.8 0.96407 0.96485 0.96562 0.96638 0.96712 0.96784 0.96856 0.96926 0.96995 0.97062
466
Z 0.00 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09
1.9 0.97128 0.97193 0.97257 0.9732 0.97381 0.97441 0.975 0.97558 0.97615 0.9767
2 0.97725 0.97778 0.97831 0.97882 0.97932 0.97982 0.9803 0.98077 0.98124 0.98169
2.1 0.98214 0.98257 0.983 0.98341 0.98382 0.98422 0.98461 0.985 0.98537 0.98574
2.2 0.9861 0.98645 0.98679 0.98713 0.98745 0.98778 0.98809 0.9884 0.9887 0.98899
2.3 0.98928 0.98956 0.98983 0.9901 0.99036 0.99061 0.99086 0.99111 0.99134 0.99158
2.4 0.9918 0.99202 0.99224 0.99245 0.99266 0.99286 0.99305 0.99324 0.99343 0.99361
2.5 0.99379 0.99396 0.99413 0.9943 0.99446 0.99461 0.99477 0.99492 0.99506 0.9952
2.6 0.99534 0.99547 0.9956 0.99573 0.99585 0.99598 0.99609 0.99621 0.99632 0.99643
2.7 0.99653 0.99664 0.99674 0.99683 0.99693 0.99702 0.99711 0.9972 0.99728 0.99736
2.8 0.99744 0.99752 0.9976 0.99767 0.99774 0.99781 0.99788 0.99795 0.99801 0.99807
2.9 0.99813 0.99819 0.99825 0.99831 0.99836 0.99841 0.99846 0.99851 0.99856 0.99861
3 0.99865 0.99869 0.99874 0.99878 0.99882 0.99886 0.99889 0.99893 0.99896 0.999
3.1 0.99903 0.99906 0.9991 0.99913 0.99916 0.99918 0.99921 0.99924 0.99926 0.99929
3.2 0.99931 0.99934 0.99936 0.99938 0.9994 0.99942 0.99944 0.99946 0.99948 0.9995
3.3 0.99952 0.99953 0.99955 0.99957 0.99958 0.9996 0.99961 0.99962 0.99964 0.99965
3.4 0.99966 0.99968 0.99969 0.9997 0.99971 0.99972 0.99973 0.99974 0.99975 0.99976
3.5 0.99977 0.99978 0.99978 0.99979 0.9998 0.99981 0.99981 0.99982 0.99983 0.99983
467 Reliability Engineering
.001 .005 .010 .025 .050 .100 .900 .950 .975 .990 .995 .999
1 1.83 7.88 6.63 5.02 3.84 2.71 0.02 0.00 0.00 0.00 0.00 0.00
2 13.82 10.60 9.21 7.38 5.99 4.61 0.21 0.10 0.05 0.02 0.01 0.00
3 16.27 12.84 11.34 9.35 7.81 6.25 0.58 0.35 0.22 0.11 0.07 0.02
4 18.47 14.86 13.28 11.14 9.49 7.78 1.06 0.71 0.48 0.30 0.21 0.09
5 20.52 16.75 15.09 12.83 11.07 9.24 1.61 1.15 0.83 0.55 0.41 0.21
6 22.46 18.55 16.81 14.45 12.59 10.64 2.20 1.64 1.24 0.87 0.68 0.38
7 24.32 20.28 18.48 16.01 14.07 12.02 2.83 2.17 1.69 1.24 0.99 0.60
8 26.13 21.95 20.09 17.53 15.51 13.36 3.49 2.73 2.18 1.65 1.34 0.86
9 27.88 23.59 21.67 19.02 16.92 14.68 4.17 3.33 2.70 2.09 1.73 1.15
10 29.59 25.19 23.21 20.48 18.31 15.99 4.87 3.94 3.25 2.56 2.16 1.48
11 31.26 26.76 24.72 21.92 19.68 17.28 5.58 4.57 3.82 3.05 2.60 1.83
12 32.91 28.30 26.22 23.34 21.03 18.55 6.30 5.23 4.40 3.57 3.07 2.21
13 34.53 29.82 27.69 24.74 22.36 19.81 7.04 5.89 5.01 4.11 3.57 2.62
14 36.12 31.32 29.14 26.12 23.68 21.06 7.79 6.57 5.63 4.66 4.07 3.04
15 37.70 32.80 30.58 27.49 25.00 22.31 8.55 7.26 6.26 5.23 4.60 3.48
16 39.25 34.27 32.00 28.85 26.30 23.54 9.31 7.96 6.91 5.81 5.14 3.94
470
Table E continud
α .001
.005 .010 .025 .050 .100 .900 .950 .975 .990 .995 .999
17 4.79 35.72 33.41 30.19 27.59 24.77 10.09 8.67 7.56 6.41 5.70 4.42
18 42.31 37.16 34.81 31.53 28.87 25.99 10.86 9.39 8.23 7.01 6.26 4.91
19 43.82 38.58 36.19 32.85 30.14 27.20 11.65 10.12 8.91 7.63 6.84 5.41
20 45.32 40.00 37.57 34.17 31.41 28.41 2.44 10.85 9.59 8.26 7.43 5.92
21 46.80 41.40 38.93 35.48 32.67 29.62 13.24 11.59 10.28 8.90 8.03 6.45
22 48.27 42.80 40.29 36.78 33.92 30.81 14.04 12.34 10.98 9.54 8.64 6.98
23 49.73 44.18 41.64 38.08 35.17 32.01 14.85 13.09 11.69 10.20 9.26 7.53
24 51.18 45.56 42.98 39.36 36.42 33.20 15.66 13.85 12.40 1.86 9.89 8.09
25 52.62 46.93 44.31 40.65 37.65 34.38 16.47 14.61 13.12 11.52 10.52 8.65
26 54.05 48.29 45.64 41.92 38.89 35.56 17.29 15.38 13.84 12.20 11.16 9.22
27 55.47 49.64 46.96 43.19 40.11 36.74 18.11 16.15 14.57 12.88 11.81 9.80
28 56.89 50.99 48.28 44.46 41.34 37.92 18.94 16.93 15.31 13.56 12.46 10.39
29 58.30 52.34 49.59 45.72 42.56 39.09 19.77 17.71 16.05 14.26 13.12 10.99
30 59.70 53.67 50.89 46.98 43.77 40.26 20.60 18.49 16.79 14.95 13.79 11.59
40 73.40 66.77 63.69 59.34 55.76 51.81 29.05 26.51 24.43 22.16 20.71 17.92
50 86.67 79.49 76.15 71.42 67.50 63.17 37.69 34.76 32.36 29.71 27.99 24.67
60 99.61 91.95 88.38 83.30 79.08 74.40 46.46 43.19 40.48 37.48 35.53 31.74
70 112.32 104.21 100.43 95.02 90.53 85.53 55.33 51.74 48.76 45.44 43.28 39.04
80 124.84 116.32 112.33 106.63 101.88 96.58 64.28 60.39 57.15 53.54 51.17 46.52
90 137.20 128.30 124.12 118.14 113.15 107.57 73.29 69.13 65.65 61.75 59.20 54.16
100 149.45 140.17 135.81 129.56 124.34 118.50 82.36 77.93 74.22 70.06 67.33 61.92
471 Reliability Engineering
A Weibull distribution with A= 0 & C=1 is an exponential distribution. A Weibull distribution with A= 0 & C=2 is Rayleigh distribution
A GPD distribution with A= 0 & C=0 is an exponential distribution
474
Poisson with x
parameter e
x 0,1,2,... exp[(et 1)] e ( z 1)
x!
Geometric for success
x 1 pe t pz 1 1 p
with parameter p(1 p) x 1,2,
0 p 1 1 (1 p )e t 1 (1 p ) z p p2
Geometric for failure p
x
, p 1 p 1 p
with parameter p(1 p) x 0,1, 2, 1 (1 p )e t
0 p 1 t ln(1 p )
1 (1 p )z p p2
Table H MATLAB1 Commands related to Distributions
Parameter Density/probabality
estimator
Random numbers Inverse of CDF CDF(F(x)) function
Distribution
Beta betafit(X) betarnd(A,B,m,n,o,...) betainv(P,A,B) betacdf(x,A,B) betapdf(x,A,B)
Poisson poissfit(X) poissrnd(λ,m,n) poissinv(P, λ) poisscdf(x, λ) poisspdf(x, λ)
Binomial binofit(X,n binornd(N,P,m,n) binoinv(Y,N,P) binocdf(x,N,P) binopdf(x,N,P)
Neg Bino. nbinfit(X) nbinrnd(R,P,m,n) nbininv(Y,R,P) nbincdf(x,R,P) nbinpdf(x,R,P)
Hyp. Geo. hygernd(M,K,N,m,n) hygeinv(P,M,K,N) hygecdf(x,M,K,N) hygepdf(x,M,K,N
Gamma gamfit(X) gamrnd(n, λ,m,n) gaminv(P, n, λ) gamcdf(x, n, λ) gampdf(x,n, λ)
Lognormal lognfit(X) lognrnd(μ, σ,m,n) logninv(P,μ, σ) logncdf(x,μ, σ) lognpdf(x,μ, σ)
Chi-Squa.. chi2rnd(V,m,n) chi2inv(P,V) chi2cdf(x,V) chi2pdf(x,V)
Normal normfit(X) normrnd(μ, σ,m,n) norminv(P,μ, σ) normcdf(x μ, σ) normpdf(x, μ, σ)
Exponential expfit(X) exprnd(mu,m,n) expinv(P,mu) expcdf(x, mu) exppdf(x, mu)
Geometry geornd(P,m,n) geoinv(Y,P) geocdf(x,P) geopdf(x,P)
Weibull wblfit(X) wblrnd(B,C,m,n) wblinv(P, B,C) wblcdf(x, B,C) wblpdf(x-A, B,C)
Uniform unifit(X) unifrnd(A,B,m,n) unifinv(P,A,B) unifcdf(x,A,B) unifpdf(x,A,B)
F frnd(V1,V2,m,n) finv(P, V1,V2) fcdf(x, V1,V2) fpdf(x, V1,V2)
GEV gevfit(X) gevrnd(C,B,A) gevinv (P,C,B,A) gevcdf(x,C,B,A) gevpdf(C,B,A)
GPD gpfit gprnd gpinv gpcdf gppdf
Rayleigh raylfit(X) raylrnd(B,m,n) raylinv(P,B) raylcdf(x,B) raylpdf(x,B)
t trnd(V,m,n) tinv(P,V) tcdf(x,V) tpdf(x,V)
1
Prepared by by: Mr Mohsen Abyar : A graduate of Shahid Bahonar University of Kerman, Iran
ABOUT THE AUTHOR
The author received his B.S. in Industrial Engineering (IE) from a
University of Technology in Tehran, in 1976 and his MS degree
in IE from University of Pittsburgh(Pitt) ,PA in 1978. He was
employed as a faculty member in Kerman, Iran in 1979. He
started to continue his studies for PhD at Pitt in 1985; after 2
semesters he left USA for home; however he received PhD
from Brunel University of London in July 2006. He has taught
some courses for over 30 years