Reliability-Engineering - PDF Safe
Reliability-Engineering - PDF Safe
������������
����������������������������������� ��
������������������������
��������������������������������������������������������������������������������������������������������
�
2
��
�
3 Reliabilty Engineering
�����������������������������������������������������������������������������������������������
F O R E W O R D ...................................................... 11
Nomenclature .............................................................. 13
References ................................................449
T A B L E S .........................................455
�����������������������������������������
�
�
4
�
5 Reliabilty Engineering
C O N T E N T S
F O R E W O R D ...................................................... 11
Nomenclature .............................................................. 13
�
7 Reliabilty Engineering
�
9 Reliabilty Engineering
�
11 Reliabilty Engineering
�����������������
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
������������� �������������
Symbols �� Description ��
H,L,N,O R
active �
component system( both
active)
ℎ( )� Instantaneous failure rate
������ The reliability function of
function the switch
���� Lower Specification Level System reliability function
�������
�
13 Reliabilty Engineering
������������� �������������
� ) The reliability function of The failure rate of
ith unit component a
� The failure rate of
component b�
S,T �
In load sharing:the failure
rate for the component
S� Stress, Load
under half load�
Smax� Maximum of stress(load),
upper extreme value of S �
In load sharing:The failure
rate for the component
under full load�
Smean Average stress value
The failure rate of the
�
Smin Minimum of stress(load),
lower extreme value of S
switch
The failure rate of the
t* Duration of experiment
� standby component
while in standby mode
T� Total operating time of all
items under the lifetest The mean of a population
� Average load
U,V,W,X Average stength
The standard devaition of
U� The upper specification
level of a dimention of a
� apopulation
product Standard devation of load
sistribution
� − ) The difference of U al L
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
� � Coefficient of variation of measurement in a part
safety factor(SF)
In AGREE Method: the standard devaition of a
�
The importance factor
of the ith subsystem
sum � measurement
in the assembly
α, β, … �� composition operator
acting over a function
(in UGF Technique)
� The probability of Type I
End��� Eexample
error
End of theorem ; End of
The probability of Type II ������� proof
error
�����
The strength of the ith link
δ� (component)
� �
14
� �
�
15 Reliabilty Engineering
Chapter 1
Introduction
&
Basic
Concepts
Chap 1 Introduction and Basic Concepts 16
1
Introduction and Basic Concepts
1.1 Introduction
�
17 Reliabilty Engineering
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)
�
21 Reliabilty Engineering
MTBF=MTTR+MTTF (1-1-2)
WT
M T�BF (1-2)
N
∞
���� = ���X > ) = ���X ≥ t) = ∫ f(x)�� − � − �(x) (1-3)
where
� �� ( > ) , k = 1, 2, . . .. (1-4-1)
( ) = ∑∞ ( ), (1-4-2)
( ) = (1 − ) , k = 1, 2, . . .. (1-5)
�
23 Reliabilty Engineering
� )
( )≅ . (1-6)
∞
( )=∫ ∞
( ) . (1-7)�
∞
( ) = ∫ [1 − ( )] −∫ ∞
( ) (1-8)
∞
( ) = ∫ [1 − ( )] (1-9)
or from :
( )=∫ ( ) . (1-10)
Example 1-1
Solution
( )= = , ( )= ≅ 0.��
Example 1-2
�
25 Reliabilty Engineering
Solution
Let
Min (X,T)+Max(X,T)=X+T⟹
�(X + T) = E(X) + T
then
λ λ λ
�{����X, T)}=(�� + λ) × e +T 1−e = T+ λe
Chap 1 Introduction and Basic Concepts 26
= � � − (T + e λ � ⇒ �E(Y) = (1 − e λ �
λ λ λ
Pr(t1 X t2 ) � ( )− ( )= ( ) − ( ).
( )− ( )
( < < | > ) ( ) ( )− ( )
= =
− − ( − ) ( )
�
27 Reliabilty Engineering
� � ( ∆ ) ( ) ( )
ℎ( ) = ���∆ → (∆ )× ( )
= ( )
− = ( )
, (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) � �k)
h(k) = � . (�� − 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:
�
31 Reliabilty Engineering
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)
a
f (x )dx
Chap 1 Introduction and Basic Concepts 32
�
� �� , ≥ 0, >0 = > 0 (1-19-1)
∞
� )= ( > )=∫ � � ≥ 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
x 2
f (x )
1
e 2 2 , x (1-20-1)
2
The rate function is:
1
1 t 2 t
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.
ℎ�
��
�� � �
t=TTF
Fig. 1.1 shows the rate function of 2 normal distributions.
Example 1.3
����� − �����
�������� = ��� > ������ = �� � > �
����
Pr(Z>0.5)=1-Pr(Z<-0.5)=0.69146=69.15%
�
35 Reliabilty Engineering
- t-200002
2
1
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);
�
37 Reliabilty Engin
neering
( )= −1 (1-22-3
3)
median= (1-22-4
4)
Figure 1.3 shows
s the pdf of 2 sample lognormal
distribution
ns.
( )
= �� ( )
� � ��������������� (1-22
2-5)
= �� ( ) − . (1-22
2-6)
�
Fig 1.4 Normal and lognormal distributions�
To calculate the probabilities in this distribution proceed as
follows:
FT (t ) Pr(T t )
lnT ln t ln t
Pr(lnT ln t ) Pr( ) Pr(Z )
ln t
Then FT (t ) Z ( ) (1-22-7)
where �is the CDF of standard normal.
The reliability function or survivor function is given by :
�
39 Reliabilty Engineering
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
( )= ( = ���) ≅ 0.����
√ ,
�
41 Reliabilty Engineering
( ) /
ℎ(���� � ( )
= = 0.���� 1 =
.
�� �������� ����� .
( )= ≤ ≤
(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��������� 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
f t e B
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 ��� ���2 ��� � (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
(1-24-2)
E (X )
(1-24-3)
V ar ( X )
2
(1-24-4)
M GF (t ) ( )
t
������������������������������
Γ( + ) (1-25)
� + > 0�
Γ( )
�
49 Reliabilty Engin
neering
Fig
g. 1-10 pdf, h(t),R(t) of some distributions
Chap 1 Introduction and Basic Concepts 50
�
51 Reliabilty Engin
neering
�
Fig 1-11-2
2 Some variations of bathtub curve(Kuo&Zuo,2
2003)
Fig 1-11
1-4 Another variation of Bathtub failure rate functtion
(Nahmias, 2004, Fig.12-4)
Figure 1-11-5
5 shows different forms of bathtub failure curve due
d to
different leveels of stress on some mechanical devices.
�
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):
�
55 Reliabilty Engineering
∫ ( ) 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
�
57 Reliabilty Engineering
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 n
n 1 x x n n 1
x 1
E ( X ) x p x 1 p n p 1 p np p 1 p
n x n x n x
x 0 x x 1 x 1 x 1 x 1
n 1 n 1
y
x 1 y E (X ) np p 1 p
n 1 y
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�������0�-0.0256������
Pr(X=0)= =1.8���-6
0!
Distribution h(x )
( = ) ( ≤ )
Discrete uniform 1 x 1 x
� ∈ ��, 1 , … , n}
n n n x
Binomial
( − )
(1 − p) ( − )
∑ ( − )
∈ � � �, … , }
Poisson
x
� ∈ {0, 1 , … } e
x x k
e x 1 k
x! k 1 k ! x !(1 )
k 0 k!
Geometric
( − ) −� − ) �� p
� ∈ {0, 1 , … }
Hyper-geometric
−
max 0, m+ n -N −
− −
− −
x ∑ ( , ) −
−
min n, m
Chap 1 Introduction and Basic Concepts 62
�
63 Reliabilty Engineering
Grouped Observations
���
Chap 1 Introduction and Basic Concepts 64
.
interval (a� b� ) (a� b� ) (ai bi ) (an bn ) sum
... ..
.
Frequ. o
.. ..
and let
N (bi )
R� (bi ) (1-34)
N�
�
65 Reliabilty Engineering
ℎ( ) = (1-35-2)
N (a i ) ( t ) i
Example 1-8
Solution
b ( ) ℎ� �
i � � multiply multiply
��� � � ���� �
i ai ai ( t ) i ��� �
�
�b ) by:
��
by:
��
=
� ai t b i
�
67 Reliabilty Engineering
� �� (������ �� − �� ��
b1 �� → � ������� = �
��
�
��
� 0.���
N� N�
N ( a i ) − N ( bi )
ℎ( ) = →
N ( ai ) bi ai
h� ( t ) � �� − ��
= 0.��� � �� � 0 t 20000
����(������
ai − bi
� �� � →��
N� � (b a )
i i
20000<t 40000 � ( )=
( )( )
= 0.�� � �� �
Example 1-9
Solution
N (b i )
� )= �� b1 �� ,�ℎ� � � �
N � × (b a ) N� N ( a i ) (b a )
i i
i i
�� � � N � (b�i ai ) � � R� (b i )
�� � � ��� �
�
� � (b i ai )
�
The followin
ng figures shows the functions of Example 1--9:
�
69 Reliabilty Engineering
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
f� � � i �� , R� (bi ) NN(b�i ) , i
h� � � ��
N (b a ) N (a )(b a )
� i i i i i
Chap 1 Introduction and Basic Concepts 70
i ����� � � �� ��� � � ��
failures
���
1 0 3 185 185/(800�3)=0.0771 185/(800�3)=0.0771
���
���
2 3 6 42 42/(800���������� 42/(615����������
���
���
3 6 9 36 36/(800��������� 36/(573����������
���
���
4 9 12 30 30/(800���������� 30/(537����������
���
���
5 12 15 17 17/(800���������� 17/(507����������
���
���
6 15 18 8 8/(800���������� 8/(490����������
���
���
7 18 21 14 14/(800���������� 14/(482����������
���
���
8 21 24 9 9/(800����������� 9/(468����������
���
���
9 24 27 6 6/(800���������� 6/(459����������
���
���
10 27 30 3 3/(800���������� 3/(453����������
���
�
71 Reliabilty Engineering
Cumulative Distr .
function and
Reliability function�
Example 1-11
Solution
��(t)٬ f(t) , R(t) were calculated using Eqs. 1-35-1, 1-36-1 and
1-37. The following table shows the results:
�
73 Reliabilty Engineering
20 18 16 14 13 12 11 10 9 8 7 6
= 2 × ��
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
�
F(t) t �in 20 38 54 68 81 93 104 114 123 131 138 144
1000 hours
200 200 200 200 200 200 200 200 200 200 200 200
� �
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
�
75 Reliabilty Engineering
Solution
.
� � =0.0024 � =1− = R (195)=0.9167
( )( . ) .
.
� � =0.0050 � = 1− = R(300)=0.5595
( )( . ) .
.
� = =0.0082 � = 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
� − 0.3
1 �
�� � 0.4
1 5 (�� − 5)(�� � 0.4) 0.9327 ��� − 5)(�� − 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
10 117.5 Cannot be computed 0.9327 0.0673 Cannot be computed
End of Example
�����������������������������������������������������������
1
Prepared by:Mr M Morrdi former student of Kerman University
�
77 Reliabilty Engineering
�� 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:
�
79 Reliabilty Engineering
( )=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≤ ≤
�� ( )
(y) = 1 − [1 − F (y)] �
� �
�1− 1−
� )= 1− 1− 0≤ ≤ �
� ��� � �� 1− ≥0
→
=1 ( ) , =2 ( ) ,… , = ( )
�
81 Reliabilty Engineering
( ) = ��� � ��1−
→
( )
� )= = ( )= ≥ 0.End of proof
� = = ( )
Derivation of ( )
i.e.�the CDF of sample minimum
or � ) when →∞
� ) = ⋯ in terms of
� ( )
(y) = ���limiting � � � )
���U< ⋯)=�� .
Step 5 Calculate � ( )
(y), considering step 4 and
Example 1-15
�
83 Reliabilty Engineering
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 − )
1−
Let �− ; then
1 1 1
� ) =− − − − −⋯
2 3
Step 3
( ) → = × .
Chap 1 Introduction and Basic Concepts 86
Step 4
� ( )
(y) = ��� ( ) < ) = ≤ = ( ≤ )
⟹ ( )
( )= ( )
Step 5
( )
( )= ( )=1− , ≥0 ⟹
( )
( )=1− , ≥0
nλy
(y) = e y≥0
f ( )
0 ������
End of Example
�
87 Reliabilty Engineering
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
�
89 Reliabilty Engineering
� (x) = , � ≤ 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− ( )=
Example 1-18
Solution
�
91 Reliabilty Engineering
∫
Solving the equation = in MATLAB:
landa=127.64972
= �� � , = �� � , = , = �����, ( ) = 0.9�
Yields � ≅ ������.
�����������������������������������������������������������
�
���������������������������������
�
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(� − , 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
�
95 Reliabilty Engineering
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 20
B 19.34
20+1
1+
6 20
χ ,
= ��������0.������ = �������
��
χ ,
= ��������0.������ = �������
�
fall outside χ ,
χ ,
. Therefore the
�
∑ . .
� = ≅ ������� ( )= × = .
.
>>Data=[�
50.1
�
10.7];
>> H=kstest(Data, [Data expcdf(Data,mean(Data))] , 0.1)
H=0
�
97 Reliabilty Engineering
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
i a
F (i ) (� � � )
n b
�.�� �.��
where a �.�� b �.��
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)
n x (i ) x (i ) n x� ( i ) 2 x� ( i )
2 2 2
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. , � 11.7287
c4
or � ) = (0.������������������� = ��������
�
101 Reliabilty Engineering
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
∫ ��� + ≤ | = ) ( ) =∫ 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 − y) ( )
�
105 Reliabilty Engineering
�a − y) ≠ 0 0≤ − ≤ 1 ���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 − y) ( )
�a − y) ≠ 0 for � − y ≥ 0 ≤ .
� )
f (a) = ∫ f (a − y)f (y)�� � ∫ � e �� � �� � a ≥ 0
Solution No.2
( )=∫ (a − x) ( ) .
�
107 Reliabilty Engineering
� )
( )= = a≥0
Solution No.3
Example 1-22
Solution
� �= (a − y) ( )
�a) = (a − y) × 1
1
��) = (t)(− ) = →
√2
�
109 Reliabilty Engineering
and also
� ∗ F��� ∗ ������������������������������������������
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
�������������������������������������������������
���������������������������������
�
111 Reliabilty Engineering
− ≥ ≤ , >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
��� − < < + )> 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��
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
�����������������������������������������������������������
�
��������� �����������������������������������
�
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
�
119 Reliabilty Engineering
Solution
k
n
L Px i Px i p x i q n x i
i � xi
k
n n n x k n x k
L p x i q n x i p x �q n x � ... p q
i � x i x� xk
k k
k
n xi ( n x i )
L p i �
q i �
i � x i
k k
k
n xi (n x i )
ln(L ) ln( ) ln( p i � ) ln(q i � )
i � x i
k
n k k
ln(L ) ln x i ln p (n x i ) ln(� p ),
i � x i i � i �
k k k k
ln(L ) x i (n x i ) � p (n x i ) nk x i
� i �
i �
� i �
i �
p p � p p k k
x
i �
i x
i �
i
� nk xi
� k � MLE ( p ) i �
p kn
xi i �
End of Example
x
k � MLE ( p )
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
�
121 Reliabilty Engineering
Example 1-30
2- Methods of Moments(MOM)
i �
x k
i
Mk
n .
Methods of moments for estimating the parameters of the
distribution of random sample X is based on the fact the g the
kth moment of X could be estimated by the kth sample moment
i.e. � �= , k=1, 2��
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, ����
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
� X i n
E (X ) M � �
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
X i
E X M � n
X
E � M � �
� � X i
n
X i
� X
n
� ^ �
X i
� X i �
�
X
�
n
n
Chap 1 Introduction and Basic Concepts 124
x i
�
(x i X )� x i
�
nX �
S� i �
i �
is a good estimate for the
n � n �
variance of X because of being unbiased i.e. � )= .
�
125 Reliabilty Engineering
B-Parameter estimation
[muhat,sigmahat] = normfit(Data)
BhatandChat= wblfit(data)
���������������������������������
H = kstest(Data,CDF, )�
where
�
127 Reliabilty Engineering
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 � ���� .
>>�X=[���������
pd=makedist(distribution name ); e.g.
pd=makedist ('exponential')
pd=makedist ('Gamma')
qqplot(X,pd)
�
129 Reliabilty Engineering
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?
���� = b − c t, 0≤t≤t �
h(t) = b − c t − c (t − t ), t < ≤t �
�
133 Reliabilty Engineering
- 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
= ( ∩ …∩ ).
Assuming
g the operation of the subsystems are indepen
ndent of
each other, we
w have:
= ( )… ( )= ×…× ⇒
=∏ . (2-1)
n
Since � � �� therefore ≤ min � � ; that is the reliab
bility of
i�
Therefore
( − 1)
�1 − ) = � � (− ) + �− � �� � � ��− �
2
Ignoring hig
gher order terms, could be approximated
d as
follows if < 1:
�
139 Reliabilty Engineering
: = (1 − ) ≅ 1 − , < 1.
n n
R sys � q i q i �
i � i �
Example 2-1
Solution
=1− = 0.9949
Using approximation:
�1− = .���
Chap. 2 Static Reliability Models 140
= ( ∩ ∩ …∩ ∩ ),
Using the chain rule for factorization, this joint probability can
be rewritten as follows:
= 1 − (1 − )(1 − ) … (1 − )=1− (1 − ) ( − )
�
141 Reliabilty Engineering
End of proof
Example 2-2
Solution
1
a- = 1 − (1 − )6 ≥ 0.��� ⇒ ≥ 1 − (1 − 0.8)6 =
������� = �������
b- = 1 − (0.5) = 0.���� ⇒ = ��
�
143 Reliabilty Engineering
Solution
�
145 Reliabilty Engineering
1-Active redundant
2-Standby redundant
���������������������������������������������������������
�������������������
�
147 Reliabilty Engiineering
=
�����������
������ � ��� A ������ � ���( A ��������������� × ��� B ������
�
�� + (1 − )= + −
Therefore th
he reliability of a two-component standby sy
ystem is
the same as to
t that of a 2-component active system as given below
b :
2
R 1 (1 R 1 )(1 R 2 ) Ps�1 (2-3-1
1)
syys
�������������
������������������������������������������
������
��������������������
Fig
g. 2-3-1 RDB of n-component standby system
m
(https:///egyankosh.ac.in/bitstream/123456789/35169/1/Unit-1
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 − )
�
151 Reliabilty Engineering
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 � � �� � ��� �� �� � � � �� �� � �� � � �
a b c
a b c
= 1 − [(1 − )]
�
153 Reliabilty Engineering
��� − (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
�
155 Reliabilty Engineering
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
� � + − )
=[ + −( )( )]( )
− =( )( − )
�
157 Reliabilty Engineering
������������������������������
( = )= (1 − )
The n-k system will function if there are no more than n-k
failures. Thus the reliability is as follows:
( ≤ − )= (1 − )
�
159 Reliabilty Engineering
ℎ :
���� � ��������� − �, n, 1 − R) (2-5-1)
The system failure probabilty or system unreliabilty is calculated
from:
�− =1− ( − , ,1 − )
Alternatively,
we could say that system work as far as ��� 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:
���� � � − ��������� − �, 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−1 (
=3 3 ∫0 1 − )5−3 ���������������
R=0.9 ⟹ ���������
�
161 Reliabilty Engineering
k−1 3−1
� 5
���� � � − 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.����
Eq. 2-6-1:
=1− (2,5,0.9) = 1 − 0.�����=0.99144
End of Example
Example 2-13
The system shown in the following figure has only 4
components ��A , A , A ���� 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,
���� 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 ��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 � &2 work,
If �� &3 work,
If �� &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
� = (1 − R) (R)
x
Then
=1− (1 − ) (R)
(1 − ) (R) ( )
=1−
− +1
(1 − ) (R) ( )
− −⋯
− +2
Then
≤1− (1 − ) (2 − 7)
− +1
�
165 Reliabilty Engineering
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.��� � 0.���
= 1− ( − 1, , )
1− (2,6,0.���� = 0.��������������
Notice that:
In Reliability literature related to k-out-of-n configuration
sometimes the binary variables and are defined as:
1
� �= �
0 ℎ
⎧1 ≥
⎪
� �
⎨
⎪0 < .
⎩
����������������������������
�
167 Reliabilty Engineering
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 − ) ��
�
171 Reliabilty Engineering
����������������������������������������������������
����������������������������� (Grosh,1989 Page137)�����
Suppose a delta configuration is completely known and given
in Fig 2-7 and we would like to find a star configuration which
has the same reliability. The reliability equivalence of the 2
configurations in Fig. 2-7 requires that:
1)The reliability of the section between vortices V1 and V 2 in
star configuration ( R1 in series with R3 ) must be equivalent to
the reliability of same section in delta configuration
(subsystem� r in parallel with" in series with "); i.e.
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
R� (2-10-1)
B
[� (� r� )(� r�r� )][� (� r� )(� r�r� )][� (� r� )(� r� r� )]
1 R�
� (� r� )(� r�r� )
ABC
2 R � (2-10-2)
C
ABC
3 R � (2-10-3)
A
A � (� r r )(� r ) B � (� r r )(� r )
� � � � � �
C � (� r� r� )(� r� )
Example 2-18
�
173 Reliabilty Engineering
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:
�
175 Reliabilty Engineering
Example 2-19
Solution
End of Example
Example 2-19
Solution No.1
N : star-delta approach
The vortice
ces 1,2 and 3 in the given RBD constitute a deelta
configuration
n. The equivalent star configuration has 3
components with the following reliabilities:
where
A = B = C = l-[l-(0.9)(0.9)](l-0.9)=0.981��������
3
R =R =R = 0.981 =0.9904
1 2 3 0.981
Replacing th
he delta configuration with this star configurat
ation
�
177 Reliabilty Engiineering
The reliabilit
ity of the system based on the new RBD is cal
alculated
as follows:
RSys = 0.9904
4 ���-(1-0.9904 ��������-0.9904 ������= 0.9787.
Solution No.2
N Conditional Probability Approach
RSP =1-(1-0.9
9������������������������
Chap. 2 Staticc Reliability Models 17
78
Suppose Eleement No. 3 in the original RBD works all thee time.
Replacing it with a line would result the following RBD with
w a
R PS =[1-(1-0
0.9)(1-0.9)][1-(1-0.9)(1-0.9)]=0.9801
Rsys =RSP������3���PS��3�����������������������������785 �
���������
���������������������������������������
������������������������������������������������
�
179 Reliabilty Engiineering
Before perfo
orming the calculations, some definitions are
introduced:
Definition of
o minimal cut set
This set is
i a minimal set of components which by faiiling
guarantee th
he failure of the system(Grosh,1989 page125);; in
other words the failure of its components cause interruptiion of
all paths from
m the input to the output of the system.
Definition of
o minimal path
A minimal path
p is a minimal set of components by functiioning
ensures the system
s operation(Grosh,1989 page125).
Chap. 2 Staticc Reliability Models 18
80
Notice that
Example 2-21
2 (Grosh,1989 page 125)
Find the miinimal path and cut sets for the following RB
BD.
Solution
The minimal
al cut sets are:
2 C4C5
3 C1C3C5
4 C2C3C4
The minimal
al path are:
�
181 Reliabilty Engiineering
mple
End of Exam
Example 2-22
2 (O'connor,2003 page175)
The RBD
D of a system is given below. Show the miniimal cut
Solution
The minim
mal cut sets are (1-3 2-3 and 4 )as shown below:
The minim
mal path sets are (1-2-4 3-4 ) shown below::
Chap. 2 Staticc Reliability Models 18
82
End of exam
mple
Example 2-23
2 (����������������������
What are thee minimal cut and path sets between source a and
target b giveen below?
Solution
If all co
omponents fail but Component 1,the con
nnection
nd b does not interrupts then the set �1} is a tie set
between a an
If all compon
nents fail but Components 2& 5, the connectiion
nd b does not interrupts then the set �2,5} is a tie
between a an
set
�
Functioning of the set with minimum elements 6,4 ensur
ures the
system operaation , 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
or Rsys ,
of network N2 ( R N 2 ),we have the following inequality fo
the reliability
ty of the original system (based on Grosh, 1989 p 125 -6)
RN 2 Rsys RN 1 (2-13)
For somee reasons such as dependence of the subsysstems of
Network N1, The calculated upper bound( R N 1 ) is usuallly over-
estimated by
y this method. If R N 1 exceeds 1 set R N 1 1
Example 2-24
2
Draw th
he auxiliary networks for calculating the up
pper and
lower bound
ds for RBD of Example 2-21.
Solution a)Auxiliary network N1
Network N1 which is a parallel configuration of all the
h elements in series is shown below:
minimal path
�
185 Reliabilty Engiineering
Network N2
The reliabilit
ity 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 thee auxiliary networks for calculating the upp
per and
lower bound
ds for the following RBD. The number in thee box is
the reliability
ty 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
K 1 i B K
Ri
Solution
The minimal cut sets and the minimal paths are given in the
following table.
�
189 Reliabilty Engineering
�
�
Cut Set No. j� components�
Path No. k� components�
j=1� 1-3
k=1� 1-2-4
j=2� 2-3
k=2� 3-4
� j=3� 4
����������� ������������
��������������������������������� �����������������������������������
�
�
The reliability of the system( Rsys )lies between:
and
�����������������������������������������������
Pr( Hi ) Pr( B / Hi )
Pr( Hi | B) k (2-15)
Pr(H ) Pr(B H )
i 1
i i
�
191 Reliabilty Engineering
B2 = the event that the mean life is 30,000 kilometers, Pr(B� ) �.�
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���������������� ��
(50000)(0.10)
=Pr(Z�����������
40000-30000
Pr(A B2 )=Pr(Z> )=0.00045
30000�����
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%
�
193 Reliabilty Engineering
Exercises1
�����������������������������������������������������������
�
���������������������
Chap. 2 Static Reliability Models 194
�
195 Reliabilty Engiineering
of losing a wheel
w bolt are 0.00001, what is the reliability
y of this
bolt system??
6- -(Problemnn 7 Page 70 K&L) The system diagram given beelow
describes thee circuitry for a neutral start switch on a manu
ual
automobile transmission.
t According to the service manuaal, in
order to startt this vehicle the clutch pedal must be fully
depressed an
nd the ignition switch must be in the start posi
sition.
(a) Define reliability
r as it relates to this system
(b) Draw an
n appropriate reliability block diagram(RBD)
(c) Assumin
ng that each functional block in your RBD has a
0.0001 chanc
nce of failing, calculate the system reliability.
7- A manufa
facturer wishes to know the reliability of a skiid
protection sy
ystem Io be used on military tractor trailers. The
T
system consiists of:
(a) Two batt
ttery or generator powered sensors per wheel..
(b) One logiic unit per sensor to predict wheel skid.
(c) A comm
mand unit, which operates an electric or an en
ngine
vacuum soleenoid.
Chap. 2 Staticc Reliability Models 19
96
8) A DC bat
attery has a time to failure that is normally disstributed
with a mean
n of 30 hours and a standard deviation of 30 ho
ours,
(a) What is the 25-hour reliability?
(b) When sh
hould a battery be replaced to ensure, that theere is no
more than a 10% chance of failure prior to replacement?
(c) Two baatteries are connected in parallel to power a light.
Assuming, that
t the light does not fail, what is the 35- hour
reliability fo
or the power source?
(d) A particcular battery has been in continuous use for 30
0 hours.
What is thee probability that this battery will last an
nother 4
hours?
9)Calculate the reliability of the following two systemss, where
each compon
nent has the indicated reliability
�
197 Reliabilty Engineering
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
�
199 Reliabilty Engineering
Chapter 3
Reliability
Considerations
in Design
+UGF
Technique
Chap. 3 Reliability in Design + UGF Technique� 200
3
Reliability Considerations in Design + UGF Technique
�
201 Reliabilty Engineering
1 2 3 n
�
203 Reliabilty Engineering
�
205 Reliabilty Engineering
1
R=0.9
0.95
R=0.6
0.9 R=0.7
R=component reliability in
������������������������
System Reliability
0.85
series configuration
0.8
0.75
0.7
0.65
0.6
1 2 3 4 5 6 7 8 9 10
no. of components
.
�
207 Reliabilty Engineering
Reliability of LL redundancy
The reliability of LL configuration in Fig 3.7 is given by:
= [1 − (1 − ) ] � (3-1)
Reliability of HL redundancy
The reliability of HL configuration in Fig 3.9 is:
� = 1 − [1 − ] ، (3-2)
where
�
209 Reliabilty Engineering
�
211 Reliabilty Engineering
�
213 Reliabilty Engineering
k
� ∑
( )= � �.
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
∑
� )= � ������������������������������� − ��
i 1
4. If in the definition of the MGF � is replaced by z, then one
gets z-transform of the random variable.
Example 3.3
End of Example
�
215 Reliabilty Engineering
x= �, … , p= �, … ,
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 − ����
U (z) =� .
f ( x i j1 ,..., x n j n )
n
⊗ P Z =
pij i z
… ��� (3 − ���
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 � � ��
� � = 0.6 + 0.�� + 0.�� + 0.��
� (1, 2, 4, ���
�������������� = = (0.6, 0. ��� 0.��� 0.���
End of example
�
221 Reliabilty Engineering
( ) ( )
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
� × 3 × 3 × 2 × 2 = ��; however applying UGF technique as
performed below reduces this amount to only 26 (Livitin,2010
page11).
� ( , ) = ( , ) = +
( , ) ( , )
���� + 0.�� +
� � . � � .
( , ) ( , )
0.�� ���
� × . � � .
( , ) ( , )
+0.�� + 0.�� ⟹
( ) = 0.�� � �.�� + 0.1
( )= ( )⊗ ( )
= (0.1 + 0.5 + 0.4 ) ⊗ (0.6 + 0.4 )=
( , ) ( , ) ( , [) ( , )
0.�� +0.�� +0.3 0.2
( , ) ( , )
+0.�� + 0.�� ⟹
( ) = 0.�� � 0.��
� )= ( )⊗ ( )=
End of Example
�
223 Reliabilty Engineering
Example 3-7
The universal generating function of a system is:
�
� ) = 0.1 + 0.�� + 0.4 + 0.�� ,
Find 20-hour reliability of the system.
Solution
= ( > ��� = 0.�� � 0.�� � �.���
End of Example
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 � is the state of jth subsystem(with 2 values :either
� =1=working or � �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 � 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)
�
225 Reliabilty Engineering
= = 1− = 0 ��� = 1 = 1,2,3
( = ∩ = ∩ = )=
(1 − ) (1 − ) (1 − )
�
v�������
, , Pr( X x ) p X (x)
� 0 d���� = 1 working
� =( + − ) = [1 − (1 − )(1 − )].�
R3*(1-(1-R1)*(1-R2))= 0.8458
�
227 Reliabilty Engineering
� )= ( )= ( )=
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
����������
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����������
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:
E ( ) �S
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)
f ( y s )f (s )ds y 0
s
s
y
(4-6)
R Pr(Y 0)
y 0
fY (y ) 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
1 s-15
R= f s (s) 1-Fδ (s) ds=
� ��
1- ds
s=0 s=20 25-20 25-15
1 25
50 s=20
R= (25-s)ds=0.25
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
�����������������������������������������������������������
�
237 Reliabilty Engineering
�
Fig 4-2-1Normally distributed load and strength: interfere���
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
S
2
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
.999968�
.9999999
.999999
.97725
.99865
.99977
.5987
.6915
.7734
.8413
.8944
.9332
.9599
.9937
R=
.50
ϕ ( )
�
239 Reliabilty Engineering
�
��������� ���������������������������������������������������������������
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)
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 R=
λ Exp(λδ ) Exp(λ )
(4-9)
λ λδ
Ref:K&L page 157
2 σ λ Normal Exp( λ)
R = 1 − ��� (−μ × λ + ) (μ� σ)
(4-10)
������������������������������������������ Ref:K&L page 157
3 �� =
μδ μ
�
(4-11) N(σδﻭμδ ) �(σ ﻭμ )
σδ σ K&L page 126
� = ϕ (SM)=normcdf(SM)
Calculable in MATLAB
4 � (z)=normcdf(z)���������������������������� �����σδ ﻭμδ ) �����σ ﻭμ )
�
� �ﻭ ��� �ﻭare the parameters not
the mean and standard deviation
5
s
(a b )
x
s Gamma Gamma
R
(a )(b ) x a � (� x )b �dx , (a,�, δ � � (b,�, )
(4-13)
x �
Ref:K&L page 141
betacdf ( s s , a, b ) : MATLAB
Calculable in MATLAB
6 � = ��� > ) =
−
+ ×
Weibull
(A,B,C)
Normal
(μ� σ)
√2
Ref:K&L page
(4-14)
∞
× ∫ �����− − 0.5 + ] ,
=
142& Appendix III
Calculable in Matlab , Maple���
7 ��1−
∞
Bδ
e ���� �−( y
B
δ �
�δ − �
�δ
) ]��
����� ,
, )
����� ,
, )
In K&L is seen in
C
= d the expression which
seems to be a typo
Chap. 4 Structural Reliability Analysis 242
�
243 Reliabilty Engineering
Solution
F = 1 − R = 1 − ��(Z < ),
μ −μ δ−μ −
�� = = =
σ +σ σ +0
( )
�( ) ( )
����� �����������(4-15-1)
=1− e ���� �− � ) ]
[ )]
=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 ����� Kpa and
standard deviation of������ Kpa is applied to a structure which
Chap. 4 Structural Reliability Analysis 246
( ) �����
� �� + 1 = ���( � 1) =
( ) �����
(������2/������2 + 1) = 0.�����
�� = �� � ) − = log(60000)-0.5*.1054= 10.9494
( ) ����� �����
� �� + 1 = ���( � 1) = � 1 = 0.��
( ) ������ ������
� �� � ) − =log(100000)-0.5*.01= 11.5079
� (z)=normcdf(z) z�
�
. .
�� = = 1.�����
√ . .
�
247 Reliabilty Engineering
0.����
� �= = ������� � = ������
�
( )= −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.���) fits the
strength. Since the exponential distribution could be considered
Gamma(b=1, = 0.��), therefore according to Eq. 4-13 the
reliability of the structure is:
λs
x= λ
s �λ δ
�����
R=
������
x=0
x a-1 (1-x) b-1dx=
λ
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
������ ��� 0
x=0
x=0
or by MATLAB
R betacdf ( 0.01+0.018
0.01
,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.
�����������������������������������������������������������
�
����������������������������������������������������������������������
�
249 Reliabilty Engineering
Solution
From Eq. 4-9:
�= = = 0.��� ������������
= = , = ⟹d = d ⟹
1
=1− e
>>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
�������������������
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
�
253 Reliabilty Engineering
answer in part a.
Solution
a)
���� − ����
� � ≤ = ⎛� ≤ ⎞ = 0/������
��� � ���
⎝ ⎠
b)
L=norminv(1-0.99598,5000,400)=�3939.85N
�
255 Reliabilty Engineering
���������������������������������������������������
���������������������������
=F δ � � ) = =1− �
0.8 = 1 − ⟹ = 0.2�
�
257 Reliabilty Engineering
� Pr S � ( ) = 1 − Fδ ��) (4 − ���
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 � (X ) 1]
using MATLAB
����log(4002 /50002 +1))����
2
μ�������� σ ����5000�� ���� �������
2 2
R=Pr ������ �
Pr(ln������������������� 8.294-8.5122
0.1 =-2.182��������
or using MATLAB:
End of Example
�
259 Reliabilty Engineering
���������������������������������������������
������������������������������ ��������������
�������������������������������������
×
�≥1− ( )
(4-21)
where
� The mean of SF
� = V = ������
n
Chap. 4 Structural Reliability Analysis 260
S�
������������������ n � � ��������������(4-22-1 )
S S
�����
�
� � �
�� �
S
S
(4-22-2 )
S S S
SF
Vn(%)�
Fig 4.13 Plot of � � versus Vn for six levels of R
�
261 Reliabilty Engineering
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.����
� = � 0.��
μ ������
. ∗ .
≥1− ( . )
= 0.��.
. ∗ .
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
�
263 Reliabilty Engineering
as the upper limit for S and the lower limit for respectively
forming the interference interval [ ��� � � ��� ]. min and S max are
Step 1
from:
s max s max
R=Pr(s>δ� δ min
�δ ���� � �����
min
�� FS �u ��f �u �du (4-25)
Step 2
Step 3
�
265 Reliabilty Engineering
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
�
267 Reliabilty Engineering
n
i n n
������� 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
p� 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
�
271 Reliabilty Engineering
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 �
i Pr(s min ) Pr( Z
3
) 0.5
n
50-30 2
i �
q i Pr( s max ) � exp[ (
60
) ]=0.105.
n
0.02 0.02
0.5-( ) 0.5 p i 0.5+( ) 0.5
2 i � 2
n
0.02 0.02
0.105-( ) 0.105 q i 0.105+( ) 0.105
2 i � 2
Chap. 4 Structural Reliability Analysis 272
�
273 Reliabilty Engineering
s .t .
� � �
�� �������� ���� �������� � �� �������� ���� ��������
� � �
�� �������� ���� �������� � �� �������� ���� ��������
n
ap pi b p � .5-(0.02/2)*0.5 p1 ... p10 0.5+(0.02/2)*0.5
i �
n
aq q i b q 0.1055-(0.02/2)*0.105 q� ... q�� 0.105+(0.02/2)*0.105
i �
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:
�
275 Reliabilty Engineering
End of Example
��������������������������������������
������������������������������������������
�����������
�� � − 1. (p-1)
Therefore f mw = and:
SF
�� = − 1� (����� or
×
�� = ���������������������������
� + )×
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.
�
277 Reliabilty Engineering
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:
�
279 Reliabilty Engineering
Find R.
s 2 s2 s 2 s2
R Z ( ) exp( s ) [1 Z ( )]
s 2 s
������������������������������������������������������������������������������))�
�
281 Reliabilty Engineering
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
������������������
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
�
285 Reliabilty Engineering
����������������������������������������������
�����������������
ℎ( ) = × [ ( )] (5-2)
where
ℎ( ) = × ( )+ × ( ) (5-3)
Example 5.1
Random variable X is normally distributed with density
2
x
function(pdf) 1 and � . Find the pdf of Y.
f (x ) e 2 2
2
Solution
� � )= ,�
�������Y = ℎ( ) = × [ ( )] , ( )= ( )�
1
= ⟹ = = ( )�������� = ��
2
ln y
1 2 �
f (k ( y )) e 2
2
2
ℎ( ) = × [ ( )]=
1
× e
1
ln y
;
2 2
2
�
287 Reliabilty Engineering
Solution
4A 4A
= ⟹ �1 ��2 ;
4
ℎ( ) = × ( )+ × ( )
4A 4A
ℎ( ) = × � × �
1 4 4 1 1 1 dD 1
D A 2 | dD | 2 A 2 | | �
dA
dA A
2 2
4A 4A
[exp( ) exp(
)] ���
1 1
h (A )
A 2 2 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
��������������
( . )
� )≅ ( )+ "( ) = (2 ) + (2 ) =������ �
From Eq. 5-9:
End of Example
Chap. 5 Combinations of random variables + tolerance 290
� � ≅ � ,…, )+ ∑ ⃒ ( )����������������������
Or in vector form:
|
⎛ . ⎞
� ) ≅ ( ,…, )+ ( )…. ( ) ⎜ . ⎟����������������
⎜ . ⎟
⎝ | ⎠
⎧ ⎫
⎪ ( ,…, ) ⎪
( )≅∑ × ⃒ �����������������
⎨ ⎬
⎪ ⋮ ⎪
⎩ ⎭
�
291 Reliabilty Engineering
Solution
= = ( , ),
⃒
1 ⎛ ⎞
� )≅ ( , )+ ( )⎜ ⎟�
2
⃒
⎝ ⎠
2 2
2 �����
=− , = . ⃒ = = ����
�
1 �����
= , =0. ( , )= ⟹�
�
2 0
����� 1(0.4 0)
��� ⟹ �
( )≅ +
� �
( ) ≅ ����� = ������ = �� �
�����������������������
Chap. 5 Combinations of random variables + tolerance 292
( )≅ ( ⃒ ) + ( ⃒ )
1 −�����
= (����� + (0.4) �
5 �
End of Example
approximately :
X X Y
�
E � �����������������
Y Y Y
1 Y ������������������������
2
�� 1
E 1
Y Y Y
�
�
X X
� �
�
293 Reliabilty Engineering
Example 5.5
Solution
P A
�
����� 0.4
�� ) = ≅ � � ∗ 1+ = �����8 N / cm 2
A A � 5
⟹ E( ) = ����� �
�
�
P P
� �
( ) = Var P A A
A A P A A
����� ���� 0.4 0.4
� ∗ � − = ������
� ����� 5 5
�����������������������������
tolerance is 0.12
100 2.500= 0.003.
Sum
2
nPart
2
(5-13-2)
�
295 Reliabilty Engineering
part
sum
2
(5-13-3)
n
(U L ) part .
n
�
( ) ( )
Let��∆� �and � then:
∆
���������� � �������������������������
√
��������∆� √ ������������������������������������������
would be ∑ a ± ∆′ where:
∆′ = 2
. (5 − ����
=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.��� . End of Example
√ √
Example 5-8
�
297 Reliabilty Engineering
Solution
The mean off the distributions are equal to their nominal values.
v
Solution
Taylor seeries 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
������
�
299 Reliabilty Engineering
1
E 0 E 1*N *K *
1!
E 1 E 1* N *K * N N * E 1* K * K K * E 1* N *
*
E 0 N *K *E 1 E 1 K * N E 1*N *K 2E 1*N * K *
0
E E N *K *E E *K *N E *N *K 2N *E *K *
1 1 1 1
Since the mean of the distributions of � � equal ��* , � * , � *
* *
E N * K * EE E 1* K * NN E 1* K * N * 2N * K * E 1*
0
1
1
or Eo �� �* � * � *
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
σ = � � ���������������������������
therefore:
�
��� ����� ������
E �
�
� N � �
� ������ � k � �����
� � �
Substituting numerical values in the right hand side of Eq. 5-18
* * *
gives E ����� . Eo E� N K 40� 12 ����� is the nominal value
�
�
301 Reliabilty Engineering
Δ Δ Δ Δ
�.��� � = = � ⟹ Δ = 1.���
. 0.����� 3
�����
� �� : �� � 1.����� ���or�since�� ������ E 0 � �������� ��
��
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����� 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.� �
������������
In the following problems assume all dimensions are normally distributed and the tolerance
range is 6 sigma(��σ�
�����������������������������������������������������������
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
�
305 Reliabilty Engineering
specification of V=XYZ.
Variable
mean 4 2 1
��stand. devia 0.4 0.2 0.1
�
307 Reliabilty Engiineering
Determine th
he value which the random variable smean willl exceed
only 1.3% off the time. Also determine the value which th
he
random variaable samp will not exceed 90% of the time.
Hint: For lar
argish n a gamma distribution (n, ) could be
√
approximateed with a normal distribution ( � � � �;
R 1 100 , R 1 10 , R 2 200 , R 2 15 ,
�
309 Reliabilty Engiineering
Anyo
one pursuing his goals honestlly
does not slip up
and iff he does, he can seek a way out
o
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 310
Chapter 6
Estimation of
Mean Life &
Reliability,
Related
Experiments
&Tests
�
311 Reliabilty Engineering
6
����������������������������������������������
����������������������������
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
This equation could be used for any product with any lifetime
distribution.
6.3 ��������� 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 ������������� ������������ ��� ������������
.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
X=[�����������������������xponential', mean(X));qqplot(X,pd).
Calculation of T for = = =
���������������������������������������������������
T nt *
� ������������������
r r
∗
Where � ≥ 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
� ���������������������
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
Cycles at
No
No failure
3100, 8100
which
6700
4600
5400
2600
4700
failure
failure
occurred
Solution
� n � � � ����� � �������
T ������
M T�B F � θ = � � � ������
r ��
Solution
T 24������������� ������
*
� T nt 120000 =8571hr �
MTTF
r r 14 End of Example
������������������������������������������������������
�
319 Reliabilty Engineering
where
Example 6-4:
Solution
�������������������������������������������������������
� ) ( ) ∑
= = �����
��
�� ith value in the sample containing failure times
� ( ) The time of the rth failure
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)�����������������������������5
i �
Then according to Eq. 6-6
����������������������������������������������������
nx
( r ) �������������������������������������
r
T nx ( r )
MT�BF � ��������������������
r r
�
323 Reliabilty Engineering
(Type-I Replacement ( n r )t * x ( i )
T
censoring) � i 1 �
r r
Without 6-5 ( − ) ( ) +∑
Replacement �
Failure- Without 6-6 nx
terminated Replacement (r ) �
(Type-II r
censoring) With 6-8 T nt *
Replacement � �
r r
�
�����������������������������������������������
At the end of the subject of experiments, it is worth
mentioning that there are some life tests called accelerated life
tests for quick obtaining of lifetime data.
Conventional way for preparing the lifetime data of a product
is to place some units of it on life test under normal use
�������������������������������������������
������������������������������
�
325 Reliabilty Engineering
2T 2T (6-11)
Pr 2 2 1
2 r ,2 2 r ,1
2
Example 6.7
Solution
�
327 Reliabilty Engineering
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
∑ � �
= = = = �����
2T 2T
2
θ 2 .
χ α χ α
2r, 2r,1-
2 2
2
χ16,0.975 ���������������� �����
2
χ16,0/025 ����������������� �����
2T 2T 2������� ��������
2
θ 2 θ
χ α χ α ����� ����
2r, 2r,1-
2 2
Example 6.8:
2T 2T
Pr 2 θ 2 ���α ��� ��θ��������� ��α����
χ α χ
2r, 2 2r1 ,(1- )
α
2
�
329 Reliabilty Engineering
Solution
The test is of Type II without replacement, therefore:
r
T (n r ) x( r ) xi , r =1,x (1) =15.5 ,
i �
r
T=(n-r)x (r) + x i =(n-1)x (1) +x1 =(n-1)x (1) +x (1) =nx (1)
i=1
�����������������������������������������������������������
�
�������������������������������������������������������������
Chap 6 Estimation of Lifetime &Reliability Experiments &Tests 330
1-���� �����
MATLAB
χ 22,0.975 �������������������������
χ 22,0/975 ���������������������������
2������� 2���15.5
n=5,x (1) =15.5 θ �� θ ����
7.3778 0.0506
End of Example
�
331 Reliabilty Engineering
estimate for : � × ��
where
� is when (the time or the kilometer or the temperature
or������������������ 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.
������������������������������������������
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
�����������������������������������������������������������
�
�������������������������������������������������������������������
�����������������������������������������������������������������������������
�������������������������������������������
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 − ����
� �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.��), and to reject 95% of the
lots with mean life θ � ������� (β = 0.���� 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 �
� � � � � �
335 3 5 �� �� 3Reliabilty
5 Engineering
�� ��
/
r �n n ��n �n n n n ��n
1 0 r
� �.��� � ﻭ.�� � ��� � ﻭ.0.��
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.��ﻭβ � 0.�� α � 0.��ﻭβ � 0.��
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.��ﻭ � �.1� = 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.��ﻭ � �.��� � �.��� � ﻭ.��
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�
� �.1 = ﻭ0.��� � �.��� � ﻭ.��
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�
α � ��1 = ﻭ0.�� α � 0.��� � ﻭ.��
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�
α � ��1 = ﻭ0.�� α = 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.��� α � 0.��� � ﻭ.��
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
�����������������������������������������������������������
������������������������������
�
337 Reliabilty Engineering
= ��� ≤ ) = ∑ (1 − ) (6-16)
b)the lifetim
me is exponentially distributed then for T=100
0⟹
�����������
��������������������������������������������
������
������������
���������������������
Example 6-12
6 (OC curve:Pa plotted versus mean lifetime)
(based on Gran
ant&Leevenworth,1988 page 585)
Consideer the following acceptance sampling plan
n to be
taken from a largish lot and plot its OC curves using various
MTTF. The plan is as follows:
�
339 Reliabilty Engineering
�
341 Reliabilty Engiineering
To know
w how Pa is calculated , suppose 0.0003 ; since
n=11000 theen np ' 11000 3.3 . The approximate value
va for
Pa from Poiisson CDF table : Pa Pr( X 2) 0.359 0.360 � .
The exact vaalue of Pa is calculable from MATLAB :
= ����≤ 2) = ���������������� �.����� = 0.���
���
The followin
ng MATLAB commands plots the OC curve.
p=1/17000:.0
00001:1/1000;Pa=binocdf(2,11000,p);plot(1../p,Pa)
�
Fig. 6-2 OC curve for Example 14.15
(acceptance probabilityversus mean life).
End of exam
mple
���������������������������������������������
����
�������
�����������
�
343 Reliabilty Engineering
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) )
(15-10)(2099199-190437)
θ� i=1
� �
r-1 10-1
=1908762
(190437-190437)+(245595-190437)+....+( 2099199-190437+ )
�������
10-1
nx (1) 15�������
F0 = = =1.03
�
θ� 2767421
H 0 : 0
H1 : 0
�
345 Reliabilty Engineering
8. Reject Ho if
2r
������������������
2
2,2r �������������������������������(6-18���������������
0
0
Example 6-14
Solution
2r
2
r 10, n 15, 5%
0
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
%5,20
2
31.41 . 02 %5,20
2
H � : θ0 is strongly rejected, &
�
347 Reliabilty Engineering
S = ( ), ( ), … , ( )
. For simplicity � will be
where
x ,x
2(1) 1(1)
The occurrence time of first failure in Design 1 & 2 respectively
C r2
( n2 r2 )( x2( r2 ) x2(1) ) ( x2( j ) x2(1) )
j 1
C r1
( n1 r1 )( x1( r1 ) x1(1) ) ( x1( j ) x1(1) )
i 1
Rejection Criteria:
Reject � ��� is outside
For H1 :1 2
F α F ,2 r2 2,2 r1 2
1 2 ,2r2 -2,2r1 -2 2 �
�
349 Reliabilty Engineering
�����������
�����������������������������������������������������������
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
�
351 Reliabilty Engineering
�
353 Reliabilty Engineering
1
10. Plot 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
����������������������������������������������
�������������
( )=∏ ( ) (7-1)
Proof No.1
' m m
h sys (t)=
-R 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
'
R sys (t) j 1
h sys (t)= j i
m �
R sys (t)
R i (t )
t 1
f 1 (t ) f m (t ) m
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 ) ,
0
h (t )
dt
ln R (t )
d d m
ln R sys (t ) ln R i (t ) hi (t )
dt dt i 1
d m
hsys (t ) ln R sys (t ) h i (t ) hsys (t ) hi (t )
dt i 1
End of proof
∫ ( )
� )= =
( ) = ��� − ∑ − ∑ (7-2-1)
Let ∗ = ∑ , ∗
=∑ , = ∑ = ∗
∑ ∗
��� →∞ ( )= = = (7-2-2)
Chap 7 Dynamic Models+Availability, Markov Chain 358
Example 7-1
n n
h sys (t)= h i (t)= λ i � ��������
i=1 i=1
�
359 Reliability Engineering
-t λ
������������������������ R sys (t)=e i ��������������������(7-4-1)
1
1 1 n 1
M T BF (7-4-2)
i=1 i
n n
1
λ
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.
�
361 Reliability Engineering
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
FX (1) ( y ) 1 1 FX ( y ) (7-8)
Chap 7 Dynamic Models+Availability, Markov Chain 362
G ( y ) F (1) ( y ) 1 1 F ( y ) �
n
Where
1-FC (y)=1- 1-1-Fδ (y)
n
we could conclude that the reliability
of the n-link chain which equals ������ ��� is obtainable from:
RC f S (s ) 1 F (s ) ds
n
0 (7-9)
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
10
1 500 600
s s
R sys e e ds
0 500
56 S
1 3000 1 3000 6
R sys
500 s 0
e ds
500 56
56
End of Example
Chap 7 Dynamic Models+Availability, Markov Chain 364
�
365 Reliability Engineering
�����������������������������������������
Proof
( ≤ )= [ ( ,…, )≤ ]= ( ≤ ,…, ≤ )
( ≤ ,…, ≤ )= ( ≤ )… ( ≤ )�
( ≤ )= ( ≤ )… ( ≤ )
⇒
�− ( ) = [1 − ( )][1 − ( )] … [1 − ( )] ⇒
( ) = 1 − [1 − ( )][1 − ( )] … [1 − ( )]
��������������������������������������������������������������������
���������������������������������
n
R sys (t ) 1 (1 e t i ) t 0 ������������������(7-11)�
i 1
�
367 Reliability Engineering
( ) = 1 − (1 − ( ))(1 − ( )) ��
( )=1− 1− 1− =
1t
� e e 2t et ( 1 2 ) t 0 ��������������
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
�
1 1 1 1 1 1 1 ��� − 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
k 1
Mn .
3 1 1 11 11
n=3 MTBFsys = -3� � � � θ .
λ �λ �λ �λ �
�
369 Reliability Engineering
R s y s ( t ) = 1 -( 1 -e -λ � ) � ����������������(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 )
������������������������������������������������
1�
�
� 2
� Switch
n
�
371 Reliability Engineering
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
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-λ1� + λ1��λ2 � ���λ 2 �λ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
2 2
R 91%, R 84.5%
standby active
2 2
for t 100 , 100
5
,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
t
2
3
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: � �� , = , λ = �
= , = 0.�� ���
�
377 Reliability Engineering
t
202
2
R (t ) e t 1 t
20 2
R (20) e (1 25 25 ) 0.932
20
25
2 2
The entire system reliability is currently:
n 1
( t )
i
n
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).
���������������������������������������������������
�������������������������������������������
�
379 Reliability Engineering
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)
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 -λ� ) � = 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− = +
( )= ( − )
�e t �e t
� �
Proof:
Rsys (t)=R 1 (t)+ R 2 (t-t ) f1 (t ) dt =e 1t e 2 (t t ') 1e 1t ' dt
t t
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
�
385 Reliability Engineering
n n
1
Remember as proved earlier MTTFsys MTTFi
i 1 l 1 i
gives
where
�
387 Reliability Engineering
����������������������������������������������
���������������������������������������
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
t
(t ) e t 1
2
R 1 e s t 0 (7-����
sys s
2
1e s t dt
MTTFsys R (t )dt e t 1
0
sys
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 Eng
gineering
2
-Substituting P s 1 in the above relationships yields R (t ) and
sys
given
om Eq. 7-43 i.e. R (t ) e t 1 st
2
Calculated fro 1 e
sys
s
λs as the sw
witch failure rate
2
ated from Eq. 7-39 i.e. R (t ) e t (1 Ps t ) if thee
Or is calculat
sys
reliability of the
t switch is given as a constant value Ps .
7 1
Example 7-6
A system
m is composed of 1 active unit A, 1 red
dundant
standby unit
it B and a switch S. The lifetime of the units is
exponentiallly distributed with 1 failure/1000hr. The keey has a
constant ratee of failure per 10000 hours. Find 500-hr reeliability
of the system
m and its mean lifetime.
�����������������������������������������������������������
1
From Dr Esh hargh's pamphlet, Faculty member of Sharif University
y of Tech
, Tehran. Figurre 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:
�
��� � e−λt (1 + PS × λt)�
R���
�
393 Reliability Engineering
�
������ = e−0.5 (1 + 0.�� � ��� � ��−6 × �������������
R���
� λ .
(����� = e (1 + � � 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
sys
t
1 (1 p )t e t
(1 Ps t ) t 0
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 dt t
1 e s
switch average reliability Ps 0
(7-46)
t 0 s t
R (t ) e t 1 PS (1 e t ) (7-47)
sys
�
395 Reliability Engineering
t t
t t
MTTFsys Rsys (t )dt e
0 0 1 (1 p ) (1 e
) dt �
Ps 1 p
1 1
1 1
MTTFsys �(7-48)�
Where
Example 7-8
t
R (t ) e t 1 1 e s
2
t0
sys
s
Chap 7 Dynamic Models+Availability, Markov Chain 396
Answer :a)
t
1 e s
Ps
s 1 e t
1 1
f x f a x a f a x a f a ...
2
1! 2!
Therefore e
t is expanded as follows:
�
397 Reliability Engineering
×
=1+ !
(− )+ ⋯�
1 e t t (7-49)
t t t s t
Ps
1e s
1e s
0
e dt
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 ��� 1 e s �����λs ≅ �� ∗ ��−6 �������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
t
e t 1 1 e s
2
R 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
λ + � � �� � λ :
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).
Where
�
401 Reliability Engineering
-λ f � 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.���� ℎ �����λ �= 4 × �� � ℎ �
According to Eq.7-51:
2 g t 2g 2 g t
R sys (t ) e (e f t e )
2g f
� (����� =
× . ×
2 × 0.��� . × × . ×
e + �e −e )
� × 0.��� − 0.���
Using MATLAB�
�exp(-2*0.001*1000)+2*0.001*(exp(-0.004*1000)-exp(-2*0.001*1000))/(2*.001-0.004)
�����������������������������������
Serviceability�
�
403 Reliability Engineering
Maintainability
Maintainability function
t
M (t ) m (x )dx ����������������������
0
Solution
t
M (t ) 1 e t 1 e
End of Example��� ����
�
405 Reliability Engineering
��������������������������
(MTTR )
i i n
MTTR i 1
i
(MTTRi ) ��������������������7-53 )�������
i
i 1 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.
����������������������(OR)
�
"The term operational readiness(O.R.) is defined as the
probability that either a system is operating or can operate
Chap 7 Dynamic Models+Availability, Markov Chain 406
Solution
�
407 Reliability Engineering
This mean that the machine is ready to perform its function 70%
of the time. End of Example��� ����
�����������������
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.r.t.= total down time- administrative& logistic times=6-2=4.
8
� = = 0.�� � ���
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���
�
409 Reliability Engineering
����������������������������������������������������������
�������������������������������������������
( )t
A (t ) e ������� ��������������(7-56)
Chap 7 Dynamic Models+Availability, Markov Chain 410
or � = λ
�� (7-58-2)
�
μ λ
�
411 Reliability Engineering
process (see page 287 of Ross,1985) that A(t) in the long range
approaches to = .
�
μ
�
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 �� working hours. The average of its
actual repair time is 89.3 hr. The other subsystem failure rate is
4.78 failures per �� 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 � ��−6 ⟹
θ MTBF1=1/80.37*10^-6=12442.4
� )=exp(-(t/ ))
(26280)=exp (-26280*80.37*(10^-6))=0.1210
.
� . . .
= . .
= 0.9929
�
415 Reliability Engineering
1 T
T 0
A * (T ) A (t )dt (7-60)
1 T
A * (T )
T
0
R (t )dt (7-61)
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 )
T 0
R (t )dt TTF ~ exp( ) R (t ) e s t
1 T 1 e sT
Ps
T
0
e s t dt
sT
Since the Taylor expansion of f(x) about a is;
1 1
f x f a x a f a x a f a ...
2
1! 2!
T
Then the expansion of e about a=0 :
- T
T-0 (T-0) 2
e s =e-λs ��� + (-λ s ���λs � �T=0 � ��λ s ��λ s ���λs� ��T=0 ����
1! 2!
1 1
e- sT 1-λ s �� �λ s ��2 � �λ s ��3 �����
2 6
1 -λ �
Then for sT 1 , approximately e s 1-λs�� �λs�� and:
2
�
417 Reliability Engineering
1
1-1+λ s �� �λ s �� 2
1-e-λs � 2 1
Ps = 1- λ s � A * �T �
λs� λs � �
∗� � = 0.�� then
b)
1 �
0.95 1- λs� T T ��� MTTF
2 10
s
End of Example
������ � � , … , ) �������������������(7-63)
where
Example 7-18
The lifetime and downtimes of n independent components,
= ( ,…, )= ×. . .×
�
419 Reliability Engineering
= 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
�
421 Reliability Engineering
������������������������������������������
���������������������
�
423 Reliability Engineering
������������
�
425 Reliability Engineering
Chapter 8
Enhancement,
Optimization
&Allocation
of Reliability
Chap 8 Enhancement, Optimization & Allocation of Reliability 426
8
�������������������������������������������
������������
���������������������������������������
������������
�
427 Reliability Engineering
Example 8.1
The series system shown in Fig. 8.1 is composed of k=3
identical components with a reliability of 0.80 each.
Solution
a) R=�.8 =0.512.
�
429 Reliability Engineering
�����������������������������������������������������������
�
����������������������������������������������������������������������������� �
Chap 8 Enhancement, Optimization & Allocation of Reliability 430
������������������������������
/ ���� ( ) ��
. . ( )≤0 = 1,2, … ,
ℎ( )=0 = 1,2, … ,
����������� ≥ ��
�
431 Reliability Eng
gineering
As an illustration
i consider the following series--parallel
system show
wn in Fig 8-5.
k
Total num
mber of components in the system is n
i 1
i
. If we
k
R sys 1 (1 R i ) ni �����������������������
i 1
where
k numbber of subsystems
Ri The reliability
r of each component in Subsystem i
ni Numb ber of components in Subsystem i
Now let
k = Maxiimum budget available
Ci = The cost
c of each component in Subsystem i
Then :
Chap 8 Enhancement, Optimization & Allocation of Reliability 432
k
�������� 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
M ax R 1 (1 R ) i
n
i 1
i
s .t .
k
a
i 1
ij ni b j i 1, 2,..., k j 1, 2,.., �
0 Ri 1
n 0 n i =Integer
i
where
�
433 Reliability Engineering
� Number of materials
th
ni � Number of components in i subsystem
k
Min Z ni C i �
i 1
�����
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 (1 Ri )
ni
i 1
R 0
0<R <1 i=1,2,...,k
i
n 0
i
n = Integer
i
�
435 Reliability Engineering
��������������������������������
�����������������������������������������������������������
1
The refrence of this chapter is mainly K&L, Chap 14.
Chap 8 Enhancement, Optimization & Allocation of Reliability 436
Where
�
437 Reliability Engineering
�..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
� � (1 − R) (R)
x
�������������������������������������
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
Example 8-3
Solution
Example 8-4
Solution
� (1 − ) ( ) �
then: �
�
(1 − ) ( ) − 0.����� � 0
�
441 Reliability Engineering
∗
The ARINC method tries to choose such that(K&L p 407):
∗
+. . . . + ∗
≤ �
i 1
i
∗
= ��������������(8-7)�
Solution
�.���
�� � 0.����� w = 0.����� �
�.��� � 0.��� � 0.���
w = 0.����
To find the required failure rate() for this exponentially-
distributed-lifetime system we could write:
�
443 Reliability Engineering
∗
= (0.�����0.������ = 0.�����
∗
= (0.�����0.������ = 0.������
∗
= (0.�����0.������ = 0.�������
∗ (��� ∗ ( . )
=e =e = 0.��������������
∗ ���� ∗
=e = 0.���, R∗ (��� = 0.������ �
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).
( )( )( )
���� � [ ( )]
= 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
= �
[− ( )]
��� ∗ � ∗ ��
= � �����4 �
�����− 0.���
1 1
= = = 0.�������� � ����� � ��
�����4
��� ∗ �.�� ∗ �
= � �����2 = 0.�������
(����− 0.���
��� ∗ � ∗ ��
= � �������� = 0.������
����− 0.���
�
447 Reliability Engineering
��� ∗ �.9 ∗ 8
= � ����� = 0.������
(���(− 0.���
∗ 1−[ ( )]
=1− �
∗
1 − (0.���
=1− = 0.������
1
∗ ∗ ( . )( )
��� = e =e = 0.�����
similarly
∗
1 − (0.���
��− � �.������
0.��
∗ ∗ ( . )( )
or � e =e = 0.�����
∗ (t)
� = 0.������,
∗ (t)
�� = 0.�����
∗ ∗ ∗ ∗
� × × = 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 ����� are less than unity.End of Example
����������
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
�����������
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.,�����
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
�
453 Reliability Engineering
�����������������������������������
�
�
455 Reliability Engineering
�
�����������������
�
�
�������������������������������������������������������( 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 10
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.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 12
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 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
λ ������������������������������������������������������� 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
66 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
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
λ ������������������������������������������������������� 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
λ ������������������������������������������������������� 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����Z < )) e.g. ��� < −3.��� = 0.������ �
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
�
Table C continued e.g.������ < 1.��� = 0.������
465 ���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
2.3 0.01072 0.01044 0.01017 0.0099 0.00964 0.00939� 0.00914 0.00889 0.00866 0.00842
2.4 0.0082 0.00798 0.00776 0.00755 0.00734 0.00714 0.00695 0.00676 0.00657 0.00639
2.5 0.00621 0.00604 0.00587 0.0057 0.00554 0.00539 0.00523 0.00508 0.00494 0.0048
2.6 0.00466 0.00453 0.0044 0.00427 0.00415 0.00402 0.00391 0.00379� 0.00368 0.00357
2.7 0.00347 0.00336 0.00326 0.00317 0.00307 0.00298 0.00289 0.0028 0.00272 0.00264
2.8 0.00256 0.00248 0.0024 0.00233 0.00226 0.00219 0.00212 0.00205 0.00199 0.00193
2.9 0.00187 0.00181 0.00175 0.00169 0.00164 0.00159 0.00154� 0.00149 0.00144 0.00139�
3 0.00135 0.00131 0.00126 0.00122 0.00118 0.00114 0.00111 0.00107 0.00104 0.001
�
469 ���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
��������������������������������������������������������������
pdf MGF Mean Variance
1 ,a x b (t ) e e
tb ta
(a b ) (b a ) 2
Uniform (a,b)�
b a t (b a ) 2 12 �
Exponential with 1 1
Parmeter
x , x 0 � t 2
�
ln t
2
2
1 2
exp
2μ�σ2 σ2
Lognormal (μ��σ)� e 2 t 0 e (e -1)
t 2 2 �
�
�
473 ���Reliability Engineering
C x A C 1 ( x A )C
A B (1 1 )
x A C
(1 )
2 ( )
B (1 ) B
2 2 2 1 ( ) B B e B Weibull�
C C
C � 1 e B x A� x A
1
xA C
1 (1 C x-A ) -1- C1 valid on
1 (1 C ) valid on
B B B
� A B
A x w hen C 0 or A x when C 0 or GPD
1 C � A x A B w hen C 0 A x A B when C 0
C � C
B
0 x A C & C0
0 x A B & C0
C
1
1 x
1
A C 1 xA
x A (1 C ) e(1 C B ) C
(1 C ) C B
� � e B
for
x A B &C 0 GEV�
B C
B
x A &C 0
1
x A & C0
C C
1 x A B &C 0
C
�
�
2B2 A γB x A xA
6 � 0.05772� e e B
�
1
B
exp( x A e
B
B )
�
FT1
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
������������������������������������������������������������
Probability function p(x) MGF� Z-transform mean Variance�
n x n x
Binomial (n,p),0≤�≤� p (1 p ) x 0,1, , n
[ pe t (1 p )]n [ pz (1 p )]n np n p (1 p ) �
x
�
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, 1 (1 p ) z p
0 p 1� 1 (1 p )e t p2 �
p
Geometric for failure , p 1 p 1 p
with parameter p(1 p) x
x 0,1, 2,� 1 (1 p )e t
0 p 1� t ln(1 p )
1 (1 p )z p � p2 �
�
�
�����������������1����������������������������������
��
�� 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(μ��σ����� 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(μ��σ����� 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
���������������������
���� ������� ��������� � ���� ����� ��� ����������� ������������ ����� ����� � ��
��������������������������������������������������������������
��� ��� ����� ����������� �������������������� ���� ��� ������ � ��� ����
������������������������������������������������������������
�����������������������������������������������������������������
���������� ��� ����� ���� ���� ������ �������� ��� ��������� ����
��������������������������������������������������������������
��������������������������������
���������������������������������������������������������������������
������������ ���� ��������� ���� ����������� �������� ���������
�������� ��� ���� �������� ��� ����� ����� ���� ���� ��� �� ��������
������� ��� �� ����������� ��� ���� ��������� �������� �������
��������� ��� ��� ���� ��� ������������ ���� ������ ����
����������������� ��� ���� �������� ��� ������������ ��� ��������
�������������������������������������
�����������������������������������������
������������������������������������������
����������������
������������������������������������� �