EE5725 Spring 2010 Homework Set 2 Solutions EE5725 HW1

EE5725 HW1 SOLUTION
EE5725 Spring 2010 Homework Set 2 Solutions
Problem 1
For a given component: h(t) = λ
where t is the operating time in hours, and λ = 0.005 (hr)―1
1) What is the mean time to failure MTTF of such a component?
MTTF = 1/λ = 1/0.005 = 200 (hours)
2) Given a population of 1000 such components operating independently, what is the
expected number of them still operating after 100 hours.
Consider the probability of surviving to t time
R(t) = 1 ― Q(t) = Ns(t) / N0
Hence, Ns(t) = R(t) N0 .
Also, λ(t) =
1
1 dR(t )
F (t ) = −
i
.
R(t )
R (t ) dt
Let R(t) = x(t), it is shown that
λ(t) = −
1 dx(t )
dx(t )
i
, and hence,
= − ∫ λ (t )dt ;
x(t ) dt
x(t )
100
100
ln x(t) = − ∫0 λ (t )dt = − ∫0 0.005dt = −0.5 , x(t) = R(t) = exp(-0.5) = 0.6065.
Therefore, Ns(t) = R(t) N0 = 0.6065*1000 = 606.5 ≈ 606. (Remove the uncertain point
in the rear); It is estimated that at t = 100 hours, 606 units (that had started operating at
0 hours) are expected to survive.
Problem 2
A utility company has a transmission line which goes over a mountain and is exposed
to many storms, wind, etc. and is rather unreliable. This line has a capacity of 500 MW.
They wish to add an additional 500 MW capacity to this interconnection. Two plans
have been submitted for increasing the capacity of the interconnection:
Plan A) Add a second 500 MW capacity circuit to the line over the mountain. We
shall assume that the failure and repair rates for the second circuit are the same as for
the first.
1
The failure and repair statistics are:
λa = 1/1000, μa = 1/200
We shall assume that the failure and repair rates for line b are identical. We shall
assume that in any state (Plan A has four possible failure states) that the transition from
a state is either due to a failure or a repair of one single component. Further, the λij for a
state transition will be the value of λ for that component which failed, similarly for μ’s.
Plan B) Requires building a double circuit line, each circuit having 250 MW
capacity through a longer route around the mountain. The circuits are treated as
independent with respect to failure rates and repair rates. The circuit for plan B is as
follows:
The capacity and the failure and repair rate for line a will be the same as in plan A. The
failure and repair rate for line c are:
λc = 1/8000, μa = 1/1000
The failure and repair rates for line d are the same as line c. Plan B has 9 states and the
λij and μij are equal to the λ and μ for the component which failed or was repaired
respectively.
The utility company is very concerned that the capacity of the transmission
interconnection be equal to or greater than 500 MW if they are going to add the
capacity to the system.
FIND:
For each plan, find the probability of having 500 MW or more capacity.
Note: This will require you to find the probability of being in one or more states whose
capacity meets the 500 or more MW criteria. To do so, find the probability of each
possible “operate/failed” state for the interconnection, and then label those that meet
the criteria and add the probabilities. (that is, the criteria is met if we are in state1 OR
state2 OR state3 OR ... each meeting the criteria).
Sol:
For Plan A:
λa = λb = 1/1000, μa = μb = 1/200
The status table of transition is depicted as follows:
A
1
1
0
0
Table I. Status table of Plan A
B
power
1
1000
0
500
1
500
0
0
2
status
operate
operate
operate
fail
Let {pi} be the probability of being in State i (i = 1, 2, 3, 4)
p1=
μ a μb
λa μb
μa λb
μa λb
, p2=
, p3=
, p4 =
;
( μa + λa )( μb + λb )
( μa + λa )( μb + λb )
( μa + λa )( μb + λb )
( μa + λa )( μb + λb )
Taking λa = λb = 1/1000, μa = μb = 1/200 into account, it is easy to find
p1 =
25
5
5
1
, p2 =
, p3 = , p4 =
36
36
36
36
The operating Probability = p1 + p2 + p3 =
35
= 0.9722, which indicates the probability
36
of transmission capacity being equal to as greater than 500MW.
For Plan B:
λa = 1/1000, λc = 1/8000, λd = 1/8000 ;
μa = 1/200, μc = 1/100, μd = 1/100 .
A
1
1
0
0
0
0
1
1
Table II. Status table of Plan B
C
D
power
1
1
1000
1
0
750
1
0
250
1
1
500
0
1
250
0
0
0
0
0
500
0
1
750
status
operate
operate
fail
operate
fail
fail
operate
operate
Construct the 8*8 status matrix as follows:
λd
0
0
0
0
0
λc
⎡1− (λc + λd )
⎤
⎢
⎥
μd
1− (λa + μd )
λa
0
0
0
0
0
⎢
⎥
⎢
⎥
0
μa
1− (μa + μd )
μd
0
0
0
0
⎢
⎥
0
0
λd
1− (λc + λd )
λc
0
0
0
⎢
⎥
P=
⎢
⎥
0
0
0
μc
1− (λd + μc )
λd
0
0
⎢
⎥
0
0
0
0
μd
1− (μa + μd )
μa
0
⎢
⎥
⎢
⎥
0
0
0
0
0
λa
1− (λa + μd )
μd
⎢
⎥
0
0
0
0
0
0
1− (μc + λd )⎦⎥
λd
⎣⎢
3
Do the MATLAB Program as below:
% Plan B_Program for State Probability Calculation
% Failure rates
lam_a=1/1000;
lam_c=1/8000;
lam_d=1/8000;
% Repair rates
mu_a=1/200;
mu_c=1/100;
mu_d=1/100;
% Transition Matrix
P=[1-(lam_c+lam_d) lam_d 0 0 0 0 0 lam_c; mu_d 1-(lam_a+mu_d) lam_a
0 0 0 0 0; 0 mu_a 1-(mu_a+mu_d) mu_d 0 0 0 0; 0 0 lam_d 1(lam_c+lam_d) lam_c 0 0 0; 0 0 0 mu_c 1-(lam_d+mu_c) lam_d 0 0; 0 0
0 0 mu_d 1-(mu_d+mu_a) mu_a 0; 0 0 0 0 0 lam_a 1-(lam_a+mu_d) mu_d;
mu_c 0 0 0 0 0 lam_d 1-(mu_c+lam_d)]
I=eye(8);
A=P-I;
Trans_A=transpose(A);
prob_matrix=[1 0 0 0 0 0 0 0]';
Trans_A(1,:)=1; % Replace the first row of Trans_A with 1
% State Probability Matrix
state_prob=inv(Trans_A)*prob_matrix
% Operating State (>=500MW)
P_operate =
state_prob(1)+state_prob(2)+state_prob(4)+state_prob(7)+state_prob(
8)
% Failed State (<500MW)
P_failed = state_prob(3)+state_prob(5)+state_prob(6)
4
>>P =
0.9998
0.0001
0
0
0
0
0
0.0001
0.0100
0.9890
0.0010
0
0
0
0
0
0
0.0050
0.9850
0.0100
0
0
0
0
0
0
0.0001
0.9998
0.0001
0
0
0
0
0
0
0.0100
0.9899
0.0001
0
0
0
0
0
0
0.0100
0.9850
0.0050
0
0
0
0
0
0
0.0010
0.9890
0.0100
0.0100
0
0
0
0
0
0.0001
0.9899
state_prob =
0.8129
0.0102
0.0020
0.1626
0.0020
0.0000
0.0001
0.0102
P_operate =
0.9959
P_failed =
0.0041
Verify that: P_operate + P_failed = 1.
Compare the Probability of operating, P_operate(B) = 0.9959 > P_operate(A) = 0.9722, it
clearly shows that Plan B offers an improved reliability over Plan A.
5