CC BY
16
Optimal System for Five Units Serial Systems under Partial
and Complete Failure
Ibrahim Yusuf*1 and Abdullahi Sanusi2 •
1 Department of Mathematical Sciences, Bayero University, Kano, NIGERIA 2School of Continuing Education, Bayero University, Kano, NIGERIA
[email protected] [email protected]
Abstract
The present paper studies and compared some reliability characteristics of series-parallel systems containing five units each under partial and complete failure. Four different system configurations are considered in this paper. It is assumed that both the failure and repair rates of each system configuration follow exponential distribution. The steady-state availability, busy period of repairman due to partial and complete failure, profit function, mean time to failure (MTTF) have been derived, examined and compared. The system configurations are compared analytically in terms of their availability and mean time to failure (MTTF). Cost-benefit measure has been evaluated for all the system configurations. The computed results are presented in tables and figures. From the analysis, system configuration II is observed to be the optimal configuration.
Keywords: Reliability, availability, standby, partial, complete failure, configuration, cost benefit.
I. Introduction
In reliability analysis, the performance evaluation of repairable systems is a matter of great importance. Maintaining the reliability of the system is indispensable. System performance can be measure through some reliability characteristics such as availability, mean time to failure, profit and benefit-cost analysis. The system availability of some engineering systems depends on the system structure, preventive maintenance, redundancy and also on the component availability. To affirm system failure and high system performance of complex systems, it is necessary to have a system component of higher availability. Generally, increasing redundant units or using units with high availability can also enhance system performance. Performance of systems increases significantly through redundancy optimization, using components with high availability, system's structural design and maintenance through repair and preventive maintenance.
To achieve high system reliability and availability, the system must be maintained at the highest order. To achieve this end, numerous researchers have designed different types of mathematical models to study and compare their reliability, availability and mean time to failure. For instance, Singh and Abdul Kareem [8] discussed the cost assessment of complex repairable systems consisting two subsystems in series configuration using Gumbel Hougaard family copula. Berk et al [2] have discussed the reliability assessment of safety-critical sensor information. Sanusi et al [12] have recently studied the performance evaluation of an industrial configured as series-parallel system. Wang et al [16] have presented the reliability analysis of two-dissimilar unit warm standby repairable system with priority in use. Singh and Ayagi [15] discussed the study of availability of standby complex system under waiting repair and human failure using Gumbel-Hougaard family copula distribution.
Ibrahim Yusuf, Abdullahi Sanusi
OPTIMAL SYSTEM FOR FIVE UNITS SERIAL SYSTEMS RT&A, № 2 (62) UNDER PARTIAL AND COMPLETE FAILURE_Volume 16, June 2021
Harish Garg [4] discussed the study of the multi objective non-linear programming problem for reliability (GSA) and the results have been compared with the results computed by practice swarm optimization (PSO) methodology. Malik and Tewari [10] analyzed the performance of a system and maintenance priorities decision for the water flow system of a coal-based thermal power plant. Kumar et al [1] have recently studied the reliability analysis of a redundant system with 'FCFS' repair policy subject to weather conditions. Niwas and Garg [11] presented an approach for analyzing the reliability and profit of an industrial system based on the cost-free warranty repair policy. More recently, Sanusi and Yusuf [13] have presented the study of cost analysis of 2-out-of-4 system connected to two-unit parallel supporting device for operation. Mortazavi et al. [9] have evaluated the MTBF and other reliability parameters for a 2-out-of-3: G redundant repairable systems with common cause failures considering fuzzy rates for failures and repair via a case study of a centrifugal water pumping system. Saini et al [14] have investigated microprocessor systems using RAMD approach. Yang et al [19] discussed the reliability assessment of system with inconsistent priors and multi-level data. Gahlot et al. [5] investigated the performance assessment of serial system with different types of failure and repair policy. Zhang [23] dealt with the reliability analysis of computer networks based on intelligent cloud computing methods. Zhao et al [24] have discussed the reliability analysis of aero-engine compressor rotor system considering cruise characteristics. Ibrahim et al [6] have studied the reliability assessment of complex system consisting two subsystems connected in series configuration using Gumbel-Hougaard family copula distribution. Kakkar et al [7] have examined availability analysis of two parallel unit system under the provision of maintenance. Yusuf et al [22] have analyzed some reliability characteristics of a linear consecutive 2-out-of-4system connected to 2-out-of-4 supporting device for operation.
Some research works in the field of reliability and performance analysis of systems with standby components/units have shown that optimality among the systems under considerations is not unique and it depends on the value of some parameters. Some studies such as Wang and Learn
[17], Wang et al. [18], EL-Sherbeny [3] and Yen and Wang [20] did take into the effect of cost benefit such as cost/availability and cost/mean time to failure on system reliability. Wang and Learn [17], Wang et al. [18], EL-Sherbeny [3] and Yen and Wang [20] have studied cost benefit analysis of various standby systems in which the optimality among configurations in the study depends only on particular parameter using cost/MTTF and depends on the other parameter using cost/availability. The present paper is motivated by the work of Wang and Learn [17], Wang et al.
[18], EL-Sherbeny [3] and Yen and Wang [20] to study reliability of four 30MW power plant systems consisting of five units each arranged in series parallel and to determine the unique optimal system among the systems under study.
The contributions of this paper are as follow:
(i) To develop the explicit expressions for availability, busy period of repairman due partial and complete failure, mean time to failure and profit function.
(ii) To perform analytical comparison between the systems in order to rank them in terms of their availability and mean time to failure.
(iii) To study and compare the four systems in terms of their profit and cost benefit.
(iv) To determine the optimal system among the systems with cold standby.
The structure of this paper is organized as follows. Section 2 gives the description of the systems considered and their reliability block diagrams. Section 3 deals with the formulation of the models. Comparison between the systems analytically in terms of their availability and mean time to failure and numerically in terms of their profit and cost benefit are presented in Section 4. Conclusions are given in Sections 5.
II. Description of the System Configurations
The present study considers 3 0Af 11/ power plant arranged in the following four series parallel configurations as shown in Figures 1-4 below: System configuration I is a series parallel system
Ibrahim Yusuf, Abdullahi Sanusi
OPTIMAL SYSTEM FOR FIVE UNITS SERIAL SYSTEMS RT&A, № 2 (62)
UNDER PARTIAL AND COMPLETE FAILURE_Volume 16, June 2021
which has three lOAflV primary units and two 10Mil/ cold standby components. System configuration II is a series parallel system having two 15Af 11/ primary components and three 1SAf 11/ cold standby components. System configuration III is a series parallel system with three 1QMW primary units and two 10MW cold standby units. System configuration IV is a series parallel system which consists of three subsystems with two subsystems arranged in parallel and serial to the other subsystem with two 15MW primary components and three 10MW cold standby components. It is assumed that all switchover time are instantaneous and switching is perfect. It is also assumed that the switch from standby to operation is perfect. Each of the primary units fails independently of the state of the others and has an exponential failure time with parameter P0 and is replace with cold standby unit if available while the failed unit is immediately sent for repairs and the time to repair is exponential with parameter a0. All failures are assumed to be repairable. System failure occurs
when all units in the same subsystem have failed. A failure is partial if the system has not failed completely otherwise the failure is complete (system failure).
10 10
10
10 10
Figure 1: Reliability block diagram of System Configuration I
Figure 2: Reliability block diagram of System Configuration II
Figure 3: Reliability block diagram of System Configuration III
Figure 4: Reliability block diagram of Configuration IV
III. Reliability Models Formulation Formulation of System Configuration I
Mean time to failure (MTTF) Analysis of System Configuration I
Let the probability that the system is in state i at time t be ht (t) and H(t) = [h(t),h(t),---,hi(t)] be the probability row vector timet with initial conditions
11, k = 0
h (0) = \
[0, k = 1,2,3,...,11
The differential-difference equations derived from system configuration I are given by:
d h0 (t) = -3Poh0 (t) + «0h (t) + a0h2 (t) + «0^0 (t)
d_ dt d_ dt d_ dt
h(t) = - (3ß0 + «0 ) h (t) + ß0h0 (t) + «0h3 (t) + «0h5 (t) + «0h9 (t)
— h2 (t) = - (3ß0 + «0 ) h2 (t) + ß0h0 (t) + «0h3 (t ) + «0h4 (t ) + «0^1 (t )
— h (t) = - (3ß0 + « ) h3 (t) + M (t) + ß0h2 (t) + «0h6 (t ) + «0h7 (t ) + «0h8 (t )
dhi(t) = -«0h (t) + M (t)
i = 4,5,6,...,11 and j = 0,1,2,3
Equation (1) can be written in matrix form as:
H'(t)= MH (t)
Where
M =
f-3ßc «0 0 0 0 0 0 0 0 «0 0 1
ßc -¿0 0 «0 0 «0 0 0 0 «0 0 0
ßc 0 -¿0 «0 «0 0 0 0 0 0 0 «0
0 ß0 ß0 ¿1 0 0 «0 «0 «0 0 0 0
0 0 ß0 0 -« 0 0 0 0 0 0 0
0 ßo 0 0 0 -« 0 0 0 0 0 0
0 0 0 ßo 0 0 -« 0 0 0 0 0
0 0 0 ßo 0 0 0 -« 0 0 0 0
0 0 0 ßo 0 0 0 0 -« 0 0 0
0 ß0 0 0 0 0 0 0 0 -« 0 0
ß0 0 0 0 0 0 0 0 0 0 -« 0
I 0 0 ß0 0 0 0 0 0 0 0 0 -«0 )
(1)
(2)
To compute the MTTF of system configuration I, the procedure requires deleting rows and columns of absorbing states of matrix M and take the transpose to produce a new matrix, Q as
adopted in Wang and Kuo [17], Wang et al [16] and Wang et al [18]. The time expected to reach
the absorbing state is calculated from:
MTTF = H ( 0) (-Q-1 ) (1,1,1,1)"
Thus, the MTTF expression for system configuration I is:
mttf =
2«0 +11«0/0 + 17/
/ (2« +15^0^0 + 27/ )
Where H (0) = [1,0,0,0] and Q =
-3&0 &0 &0 0 1
«0 -^0 0 &0
«0 0 &0
0 «0 «0 -¿1 )
(3)
(4)
Availability and Busy period of System Configuration I
To compute the availability of system configuration I, the differential difference equation given in (2) are expressed in the form:
f K01 f «0 «0 0 0 0 0 0 0 0 «0 0 1 fK1
K /0 -¿0 0 «0 0 «0 0 0 0 «0 0 0 K
K /0 0 -¿0 «0 «0 0 0 0 0 0 0 «0 h2
K' 0 /0 /0 -¿1 0 0 «0 «0 «0 0 0 0 K
K4 0 0 /0 0 -«0 0 0 0 0 0 0 0 K4
K5 0 / 0 0 0 -«o 0 0 0 0 0 0 K
K6 0 0 0 / 0 0 -«o 0 0 0 0 0 K
Kj 0 0 0 /0 0 0 0 -«0 0 0 0 0 K
K 0 0 0 / 0 0 0 0 -«o 0 0 0 K
K 0 / 0 0 0 0 0 0 0 -«0 0 0 K
K / 0 0 0 0 0 0 0 0 0 -«0 0 K
V K 1) I 0 0 /0 0 0 0 0 0 0 0 0 -«0 ) V )
The steady state availability (i.e. the sum of the probabilities of all the operational states), busy period due to partial failure and complete failure are respectively given by:
AFi (m) = h0 (M) + K,(M) + K (M) + h (M) (5)
bi=h H+k (»)+k (m) (6)
Bh2 = h4 (M) + K (M) + h6 (M) + k (M) + K (M) + K (M) + ko (M) + hi (M) (7)
All the derivatives of state probabilities are set equal to zero in the steady state, therefore equation (2) becomes:
MH (t )T = 0 (8)
Which is
—;A «0 «0 0 0 0 0 0 0 0 «0 0 1 f K 1 f 01
A -¿0 0 «0 0 «0 0 0 0 «0 0 0 h 0
A 0 -¿0 «0 «0 0 0 0 0 0 0 «0 K2 0
0 A0 A0 -¿1 0 0 «0 «0 «0 0 0 0 K 0
0 0 A0 0 -« 0 0 0 0 0 0 0 K 0
0 A 0 0 0 -« 0 0 0 0 0 0 K 0
0 0 0 A 0 0 -« 0 0 0 0 0 K 0
0 0 0 Ao 0 0 0 -« 0 0 0 0 K 0
0 0 0 A 0 0 0 0 -« 0 0 0 K 0
0 A0 0 0 0 0 0 0 0 -« 0 0 K 0
a 0 0 0 0 0 0 0 0 0 -« 0 K 0
0 0 A0 0 0 0 0 0 0 0 0 —«0 0 v Kn0 V 0 0
Using the following normalizing condition:
ii
S K W = 1 (9)
n=0
To compute the state probabilities kt(t) i = 0,1,2, ...,11, (9) is substituted in the last of (8) to give:
-;A0 «0 «0 0 0 0 0 0 0 0 «0 01 f ¿0 1 f 01
A0 —¿0 0 «0 0 «0 0 0 0 «0 0 0 K 0
A0 0 —¿0 «0 «0 0 0 0 0 0 0 «0 h2 0
0 A0 A0 —¿1 0 0 «0 «0 «0 0 0 0 K 0
0 0 A0 0 -« 0 0 0 0 0 0 0 K 0
0 A0 0 0 0 -« 0 0 0 0 0 0 K 0
0 0 0 A0 0 0 -« 0 0 0 0 0 K 0
0 0 0 A0 0 0 0 -« 0 0 0 0 K 0
0 0 0 A 0 0 0 0 -« 0 0 0 K 0
0 A0 0 0 0 0 0 0 0 -« 0 0 K 0
a 0 0 0 0 0 0 0 0 0 -« 0 K 0
0 0 A0 0 0 0 0 0 0 0 0 —«0 0 V v10
(10)
Solving (10) using MATLAB package to obtain h (t) , the expressions for the steady-state availability, busy period due to partial failure and complete failure given in (5) to (7) are respectively given by:
A (x) = + 2^oA +^<A2 F1 ( ) a\ + 3«oA + 5«o^o2 + 3$3
Bu („) = . 2^o2A +M?
«03 + «0A0 + 5^0^02 +
, sa;+«02A+4«0A
0 0 0 2
«0;+;«02a+5^0 A2+A
(11)
(12)
Profit Analysis of System Configuration I
The units are exposed to corrective maintenance due to partial and complete failure, while the repairman is busy performing maintenance action to the failed units. Let K0, K and K2 be the revenue generated when the system is in working state and no income when in failed state, cost of each repair due to partial and complete failure respectively. The expected total profit of system configuration I per unit time incurred to the system in the steady-state is given by: Profit =total revenue generated - cost incurred by the repair man due to partial failure - cost incurred due to complete failure.
PFi = K! -(KBK + KBh2 ) (14)
Formulation of System Configuration II
Mean time to failure (MTTF) Analysis of System Configuration II
Let H(t) =[h0(t),hl(t),h2(t\...,hl0(t)] be the probability row vector at time t. The initial
11, j = 0 condition is given by: h, (0) = <
[0, j = 1,2,3,...,10
The corresponding set of differential-difference equations for system configuration II are expressed as:
Where:
M„
( t ) = M2 H ( t ) (15)
-2A, a0 a0 0 0 0 0 0 0 0 0 N
A0 -( ¥0 +a0 ) 0 0 a0 0 0 a0 0 0 0
A0 0 -( ¥0 +a0 ) a0 a0 0 0 0 0 0 0
0 0 A> -( 2 A) +a0 ) 0 a0 a0 0 0 0 0
0 A0 A0 0 -( 2A0 +2a0 ) a0 0 0 a0 0 0
0 0 0 A0 A> -( 2A0 + 2a0 ) 0 0 0 a0 a0
0 0 0 A 0 0 a 0 0 0 0
0 A 0 0 0 0 0 -a0 0 0 0
0 0 0 0 0 0 0 0 -a0 0 0
0 0 0 0 0 A 0 0 0 -a0 0
0 0 0 0 0 A 0 0 0 0 a
Using similar procedure presented in subsection 3.1.1, the expression for the mean time to failure of system configuration II is obtained through:
mttf2 = h ( 0)(-Q-1 )(1,1,1,1,1,1)T (16)
Thus,
mttf2 =
4a05 + 29alP0 + 97a03/ +119a0A03 + 211a0/ +100/
/ ( 4a0 + 25a0 A + 76a2 / +112^/ + 64,A )
Where H (0) = [1,0,0,0,0,0] and
(17)
Ö2 =
a ß0 0 0 0
«0 -( 2^0 +«0 ) 0 0 ß0 0
«0 0 -( 2ß +«0 ) ß0 ß0 0
0 0 «0 -( 2ß0 +«0 ) 0 ß0
0 «0 «0 0 -( 2ß0 + 2«0 ) ß0
0 0 0 «0 «0 -( 2ß0 + 2«0 ),
Availability and Busy period Analysis of System Configuration II To compute the availability of system configuration II, the differential-difference equations given in (14) are expressed in the form:
K '-2ßc «0 «0 0 0 0
K ßc -( 2ß0 +«0 ) 0 0 «0 0
K ßc 0 -( 2ß0 +«0 ) «0 «0 0
K 0 0 ß0 -( 2ß0 +«0 ) 0 «0
K 0 ß0 ß0 0 -( 2ß0 +2«0 ) «0
K = 0 0 0 ß0 ß0 -( 2ß0 +
K 0 0 0 ß0 0 0
K 0 ß0 0 0 0 0
K 0 0 0 0 0 0
K 0 0 0 0 0 ß0
k ; ,0 0 0 0 0 ß0
0 0
0 «
0 0
0
0
0 0 0 0
0 -«
00
00
0 0 N fK
0 0 K
0 0 K
0 0 K
0 0 K
«0 «0 K
0 0 h
0 0 K
0 0 K
-«0 0 K
0 -«0 ; , K,
The steady state availability (the proportion of the time the system is functioning or equivalently the sum of the probabilities of operational state), busy period due to partial failure and complete failure are given by:
av1 (ro) = h (ro)+h (ro)+h (ro)+h (ro)+K (ro)+h (ro) (18)
Bh3 = K (ro) + h (ro)+ h (ro) + h4 (ro) + h5 (ro) (19)
Bh 4 (ro) = h (ro)+h (ro)+h (ro)+h (ro)+ho (ro) (20)
In the steady state, the derivatives of states probabilities become zero and therefore (15) becomes:
m h (t )r = 0
In matrix form, we have:
(21)
(j
"2ßö «0 «0 0 0 0 0 0 0 0 0 N f K '01
ß0 -(2ß0 +«0 ) 0 0 «0 0 0 «0 0 0 0 K 0
ß 0 -(2ß0 +«0 ) «0 «0 0 0 0 0 0 0 K 0
0 0 ß0 -( 2ß0 +«0 ) 0 «0 «0 0 0 0 0 k 0
0 ß0 ß0 0 -( 2ß +2«0 ) «0 0 0 «0 0 0 K 0
0 0 0 ß0 ß0 -( 2ß0 +2«0 ) 0 0 0 «0 «0 h = 0
0 0 0 ß 0 0 -« 0 0 0 0 K 0
0 ß 0 0 0 0 0 -« 0 0 0 K 0
0 0 0 0 0 0 0 0 -«0 0 0 h 0
0 0 0 0 0 ß 0 0 0 -« 0 K 0
0 0 0 0 0 ß 0 0 0 0 -«0; v ^ ;
Using the following normalizing condition
a
0
0
0
a
0
a
V
10
S hn H = 1
(22)
To obtain the state probabilities h (t) i = 0,1,2,.. ,10, we substitute (22) in (21) to get
'-2A0 a0 a0 0 0 0 0 0 0 0 0 1 ' h0 '0 '
A0 -( 2A +a0 ) 0 0 a0 0 0 a0 0 0 0 K 0
A0 0 -( 2A0 +a0 ) a0 a0 0 0 0 0 0 0 h2 0
0 0 A0 -(2 A +a0 ) 0 a0 a0 0 0 0 0 h3 0
0 A0 A0 0 -( 2A +2a0 ) a0 0 0 a0 0 0 h4 0
0 0 0 A A -( 2A +2a0 ) 0 0 0 a0 a0 h5 = 0
0 0 0 A 0 0 -a0 0 0 0 0 h6 0
0 A 0 0 0 0 0 -a0 0 0 0 h 0
0 0 0 0 0 0 0 0 a0 0 0 h 0
0 0 0 0 0 A 0 0 0 -a0 0 h 0
v 0 0 0 0 0 A 0 0 0 0 v h10 ; J,
(20)
Solving (23) using MATLAB package to obtain h (t) , the explicit expressions for the steady-state availability, busy period due to partial failure and complete failure are given by:
AV 2 H =
Bh0 H =
Bh4 (®) =
a04 + 2a°A + 2a02Ac2 + a A
a04 + 2a00Aз + OalA + 3a0 A + 2 A
2alA + 2a,A + aA0 a\ + 2a0 A + 3a02A2 + 3aoA0 + 2 A
2 A4 +a02A02+a0A000 a04 + 2a3A + 3a02A02 + 3a0 A + 2A4
(24)
(25)
(26)
Profit Analysis of System Configuration II
Using similar procedure presented in subsection 3.1.3, the explicit expression for profit function of system configuration II is given by:
PF2 = K0 Av2 -(KBh3 + KBh4 ) (27)
Formulation of System Configuration III
Mean time to failure of Analysis System Configuration III
Let H(t) = [/z0 (/ ), /Z| (/ ),. _ _, /zx (7 )] be the probability row vector at time t with initial conditions 11, n = 0
h (0) H
n [0, n = 1,2,3,...,8
The corresponding set of differential-difference equations for system configuration III are expressed as:
H '( t )= M H (t ) Where
n=0
M =
Г —0 —0 0 0 0 0 0 0 1
P0 -( 3p0 +-0 ) 0 —0 —0 0 0 0 0
2P0 0 -( 3p0 +-0 ) 0 0 —0 —0 0 0
0 P0 0 -(3P0 + —0 ) 0 0 0 —0 —0
0 2P0 0 0 ——0 0 0 0 0
0 0 P0 0 0 ——0 0 0 0
0 0 2P0 0 0 0 0 0
0 0 0 2P0 0 0 0 -a 0 0
V 0 0 0 P0 0 0 0 0 —— —0 )
Using similar procedure presented in subsection 3.1.1, the explicit expression for the mean time to failure (MTTFJ of System Configuration III is obtained through:
MTTF3 = H (0) (-ÇC1 ) (1,1,1, if
Thus,
MTTF = - +1 + + 57^°3
P02 (8-02 + 45-Д + 81Д2 )
Where H ( 0) = [l, 0,0,0] and Q =
-3p0 p0 2p0 0
—0 — (3P0 0 ) 0 p0
—0 0 — ( 3P0 0 ) 0
0 —0 0 — ( 3P0 0 ),
(29)
(30)
Availability and Busy period Analysis of System Configuration III
To compute the availability of system configuration III, the differential difference equations given
in (28) are expressed in the form:
Г—3P0 —0 —0 0 0 0 0 0 0 1 ГК1
К P0 —У 0 0 —0 —0 0 0 0 0 К
h2 2P0 0 —У 0 0 0 —0 —0 0 0 К2
h^ 0 P0 0 —У0 0 0 0 —0 —0 К
К = 0 2P0 0 0 0 0 0 0 0 К4
К 0 0 P0 0 0 ——0 0 0 0 К5
К 0 0 2P0 0 0 0 ——0 0 0 К
К, 0 0 0 2P0 0 0 0 ——0 0 К
К) V 0 0 0 P0 0 0 0 0 —— —0 ) V К )
The steady state availability (the proportion of the time the system is functioning), busy period due to partial failure and complete failure are given by:
= н0 +h (да)+h (да)+h (да) (31)
Bh5 = h (да) + h2 (да) + к (да) (32)
Вне (да) = к4 (да) + к (да) + h6 (да) + к (да) + к (да) (33)
In the steady state, the derivatives of states probabilities become zero and therefore (28) becomes: MH (t )T = 0 (34)
In matrix form, we have:
-3Po «o «o 0 0 0 0 0 o 1 r Ko 1 r o 1
Po 0 «0 «0 0 0 0 o K o
2Po 0 0 0 «0 «0 0 o K o
0 P0 0 0 0 0 «0 «o K o
0 2P 0 0 -«0 0 0 0 o K = o
0 0 P 0 0 -« 0 0 o K o
0 0 2P 0 0 0 -« 0 o he o
0 0 0 2P 0 0 0 -« o K o
0 0 0 P 0 0 0 0 -« «o J V K J V o J
Using the following normalizing condition
8
E h (®)=i
n=0
To obtain the state probabilities h (t) i = 0,1,2,..., 8, (35) is substituted in (34) to give:
r-3P0 «o «0 0 0 0 0 0 o 1 r v 1 r o 1
Po - ^o 0 «0 «0 0 0 0 o Pi o
2Po 0 - ^o 0 0 «0 «0 0 o P2 o
0 P0 0 - ^o 0 0 0 «0 «o P3 o
0 2P 0 0 -« 0 0 0 o P4 = o
0 0 P 0 0 -« 0 0 o P5 o
0 0 2P0 0 0 0 -«0 0 o Pe o
0 0 0 2P0 0 0 0 -«0 o Pi o
V 1 1 1 1 1 1 1 1 1J V P8 J V 1 J
(35)
(36)
Solving (36) using MATLAB package to obtain h (t) ,the explicit expressions for the steady-state availability, busy period due to partial failure and complete failure are given by: al + 3«o2Ao + aoPo2
AV3 H =
Bh5 H = Bh6 («) =
al + 3«02Ac + 9«oAo2 + 3Ac3
« Po +^olo2
«o3 + 3«o2Po + 9«oP2 + 3Po3 3Po3 + 8^oPo2
«o3 + 3^o2Po + 9^oPo2 + 3Po3
(37)
(38)
(39)
Profit analysis of System Configuration III
Using similar procedure presented in subsection 3.1.3, the explicit expression for profit function of configuration III is given by:
pf3 = k0Av3 -(kxbh5 + K2bh6 ) (40)
Formulation of System Configuration IV
Mean time to failure of Analysis System Configuration IV
Define H (t ) = [K0 (t ), K (t ), K2 (t ),•••, K (t )]
to be the probability row vector at time t with initial
conditions:
|1, n = 0
|0, n = 1,2,3,...,8.
The corresponding set of differential-difference equations for system configuration IV is expressed as:
K (0) = ■
H '(t ) = M4 H (t ) Where:
(41)
M 4 =
f g0 g0 0 0 0 0 0 0
ß0 -( 2ß0 +a0 ) 0 g0 g0 0 0 0 0
ß0 0 -g0 0 0 0 0 0 0
0 ß0 0 -( 2ß0 +a0 ) 0 g0 g0 0 0
0 ß0 0 0 -g0 0 0 0 0
0 0 0 ß0 0 -( 2ß0 +a0 ) 0 g0 g0
0 0 0 ßo 0 0 -g0 0 0
0 0 0 0 0 ßo 0 -g0 0
v 0 0 0 0 0 ß0 0 0 -g0 J
Using similar procedure presented in subsection 3.1.1, the explicit expression for the mean time to failure (MTTFJ of System Configuration IV is obtained through:
mttf4=h(o)(-q-%ii\)t
Thus,
mttfa =
gQ + 5a2^ßo + 12^olo2 +15&3 al + 5a02ß0 +12a0ß02 + 16ß03 <~2ß«
(42)
(43)
Where H ( 0) = [1,0,0,0] and Ö4 =
go 0
0
ß0
-( 2 ß0 +a0
g0 0
)
0
ß0
-( 2 ß0 +a0 ) g
0 0
ß0
-( 2ß0 +a0 )
Availability and Busy period Analysis of System Configuration IV To compute the availability of system configuration IV, the differential difference equations given in (41) are expressed in the form:
f ^ f-2ß0 a0 a0 0 0 0 0 0 0 1 f h I
K ß0 -( 2ß0 +a0 ) 0 a0 a0 0 0 0 0 K
h'2 ß0 0 -a0 0 0 0 0 0 0 k
K 0 ß0 0 -( 2ß0 +a0 ) 0 a0 a0 0 0 hi
h4 = 0 ß0 0 0 -a0 0 0 0 0 h4
h5 0 0 0 ß0 0 -( 2ß0 +a0 ) 0 a0 a0 h5
h6 0 0 0 ß0 0 0 -a0 0 0 h6
h7' 0 0 0 0 0 ß0 0 -a0 0 K
V K2j v 0 0 0 0 0 ß0 0 0 -g0 j V h8 ,
The steady state availability (the proportion of the time the system is functioning or equivalently the sum of the probabilities of operational states), busy period due to partial failure and complete failure are given by:
AV4 (») = h0 (») + + h + h5 (») (44)
Bhl = h (<») + h3 H + h5 H (45)
Bh8 H = h2 H + h4 H + h6 H + h7 H + h8 H (46)
In the steady state, the derivatives of states probabilities become zero and therefore (41) becomes:
MA H (t )T = 0 (47)
These can be expressed in matrix form as:
i _
2ßo ao ao 0 0 0 0 0 0 N r ho ' 0 '
ßo -( 2ßo +ao ) 0 a0 ao 0 0 0 0 h, 0
ßo 0 -a 0 0 0 0 0 0 0 h 0
0 ßo 0 -(2ßo +ao ) 0 a0 ao 0 0 h 0
o ßo 0 0 -a 0 0 0 0 0 hl = 0
0 0 0 ß0 0 -( 2ßo +ao ) 0 a0 ao h 0
0 0 0 ßo 0 0 0 0 0 h6 0
0 0 0 0 0 ßo 0 -a 0 0 h 0
0 0 0 0 0 ßo 0 0 -a0 y V h8 > V 0 y
Using the following normalizing condition
8
E hn (®) = 1
(48)
obtained the state probabilities h (t) i = 0,1,2, ..,8 using (48) in the last row of (47) we gel
"2ßo a 0 0 0 0 0 0 N r ho ' 0 N
ßo -( 2ßo +«0 ) 0 ao ao 0 0 0 0 h, 0
ßo 0 -a 0 0 0 0 0 0 0 h2 0
0 ß0 0 -( 2ßo ) 0 ao ao 0 0 h 0
0 ßo 0 0 0 0 0 0 0 hi = 0
0 0 0 ß0 0 -( 2ßo ) 0 ao ao h5 0
0 0 0 ßo 0 0 -a o 0 0 h6 0
0 0 0 0 0 ßo 0 -a 0 0 h7 0
1 1 1 1 1 1 1 1 -a0 y V 1 y
(49)
Explicit expressions for the steady-state availability, busy period due to partial failure and complete failure are given by:
AV4 -
Bh7 =
al + cclßo + a02p02 + cc0ß0
a0 + 2a0ß0 + 2a ß0 + 2a0ß0 + 2ß a03ß0 + 2a02ß02
a0 + 2alß0 + 2a ß0 + 2a0ß0 + 2ß
Bh8 =■
aß +a\ß3 +aoßo3 + 2ßo4
a04 + 2a03p0 + + 2a0fâ + 2$ Profit analysis of System Configuration IV
(50)
(51)
n=0
Ibrahim Yusuf, Abdullahi Sanusi
OPTIMAL SYSTEM FOR FIVE UNITS SERIAL SYSTEMS RT&A, № 2 (62) UNDER PARTIAL AND COMPLETE FAILURE_Volume 16, June 2021
Using similar procedure presented in subsection 3.1.3, the explicit expression for profit function of
system configuration IV is given by:
pf4 = K0av4 -(kbh7 + Kbh8) (53)
IV. Discussion I. Comparison between the System Configurations
Analytical Comparisons
Here, the configurations are compared analytically in terms of their availability and mean time to failure to determine the optimal configuration by taking the difference between the configurations a0, /30 > 0 using MAPLE software package.
aoPo K + ^iPo + 5«o2AO2 + 3«o$ +Ao4 ) AV1 - Ar/1 = 7---77---7 > 0 (54)
V 2 V1 / -J - 9 „ „ „9 - „Q \ / A _ -J „ „9 ^9 „ ~ ,n4\ v '
v 2 V1 ^ + 3a02ß0 + 5a0ß02 + 3ß03 )(a04 + 2a03ß0 + 3a2ß2 + 3a0ß03 + 2ßf
^ Av2 > AV1 Va0,ß0 > 0
aß ( 7a0 + 14a23ß2 + Haß2 + 6aß(0 + ß ) Ar/9 - = 7---^---r- > 0 (55)
V 2 V 3 / A _ -J „ - 9 „9 - ^ v '
(a0 + 2a03ß0 + 3a2ß2 + 3aoßo3 + 2ß4 )(a03 + 3aX + aßl + 3ß03 )
^ Av2 > Av3 Vao,ßo > 0
aßo (ao5 + 2aoßo + 2ao3ß2 + 2a] ß\ + 2aoßo4 + ß5 )
V2 V4 = (ao4 + 2ao3ßo + 3ao2ßo2 + 3aoßo3 + 2ßo4)(ao4 + 2ao3ßo + 2ao2ß2 + 2aoßo3 + 2ßo4) >
(56)
^ Av 2 > Av4 Vao, ßo > 0
A _A = aoßo (ao5 - 4ao4ßo - 6a3j0 - 5ao2ß3 - 4aoßo4 - ßj )
V3 V4 = (ao3 + 3ao2ßo + 9a0ß2 + 3ß3)(ao4 + 2aoßo + 2ao2ß2 + 2aoß3 + 2ßo4) > 57
^ AK3 (œ) > AF4 (œ) if and only if a05 > (4a4ß0 + 6a\ßl + 5alßl + 4a0ß04 + ß05 ) for some
ao >ßo
aoßo2 (ao2 - 2aoßo -ßo2 )
- Ar/1 = 7-^-77-'--7 > 0 (58)
(ao3 + 3a2ßo + 9aoßo2 + 3ß3 ) (ao3 + 3a^ß() + 5aoß2 + 3ß^> )
^ AK3 (œ) > AF1 (œ) if and only if a2 > (2a0ß0 + ß2 ) for some a0 > ß0
aoßo2 ( 2ao3 + 2ao2ßo +a0ß0^ +ß3 )
A - A =_-_-_> 0 (59)
v 4 V1 (ao4 + 2ao3ßo + 2a2ß2 + 2a{)ß + 2ß4 )(a2 + 2aß0 + 3ß2 )
^ Av4 > Av 1 Vao, ßo > 0
Using availability models of all the system configurations, it is clear from (54) - (59) that
AV 2 > AV 3 > AV 4 > AV1
Thus, the optimal system configuration is configuration II
MTTF -MTTF = 8ao7 +110g°ßo + 643go5ßo + 2125go4ßo3 + 4421g2ßo4 + 5870a2^ß^ + 4589^ + 1617ßo7 >Q 2 1 ßo2 ( 2go2 + 15goßo + 27ß ) ( 4ao4 + 25a^iß0 + 76alßl +112goßo3 + 64ß4 )
^mttf2 >mttf va,ß>o
MTTF -MTTF = 28«7 + 343a06l0 +1890^ + 6041^ +12303^ + 16138^ +12585^ + 4452# ; Q (61) 2 3 Pi (4a4 + 25a\P + 16a0P02 + 112a0P03 + 64P04)(8a02 + 45a0P0 + 81P02)
^ mttf2 > mttf v a, P > 0
MTTF -MTTF - 4a8 + 45^0 + 245^ + 839^ +1942^ + 3088^ + 3284^ + 2128^ + 640P8 >Q 2 4 p20 (a3 + 5ap0 + 12^ + 16$ )(4a4 + 25ap + 16a2p2 +112^$ + 64$)
^ mttf2 > mttf v a, P > 0
mttf -mttf = 2a0 + 2ag4P + 96a03P2 + 233a02P02 + 306^ + 162P5 > q 3 1 P2 (2a02 + 15a0P0 + 21P2 )(8a02 + 45^ + 81P2)
^ mttf3 > mttf v a, P > 0
mttf -mttf = a06 + aP. + 23aQ4P2 + 8a03P03 -- 301a0P05 - 305/ > q
3 4 P02 (a03 + 5a02P0 + 12a0P02 + 16P3 )(8a2 + 45a0P0 + 81P2)
^ MTTF3 > mttf if and only if
(a06 + 8a05P + 23a4P2 + 8a0P03) > (112a02P04 + 301a0P05 + 305P06) for some a0, P > 0
MTTF -MTTF = 4a°4 + 30a3ß. + 96ao2Po2 + 8a3p°3 + + 133p4 > o
4 1 (a°3 + 5a°2p° + 12a°p2 + 16p3 )( 2a°2 + 15a°p + 27p°2 )
^ MTTF > MTTF Va, p > 0
Similarly, using mean time to failure (MTTFJ models of all the system configurations, it is clear from (60) - (65) that
mttf > mttf > mttf > mttfx
Thus, the optimal system configuration is configuration II
II. Numerical examples The purpose of this section is to rank the system configurations in terms of their availability and mean time to failure using MATLAB software package. The results are summarized in tables below
Table 1: Ranking between the systems configurations in terms of their availability and mean time to failure.
Case Parameter Range Results Constant Value
1 0 <p < 0.3 0.3 <p < 0.6 0.6 <P < 0.9 0.9 <P < 1.2 A (œ) > A4 (œ) > A (œ) > A (œ) A (œ) > A (œ) > A (œ) > A (œ) A (œ) > A (œ) > A (œ) > A (œ) A (œ) > A (œ) > Aj (œ) > A3 (œ) a = 0.6 c0 =100000 C = 500 C = 1QQQ
2 0 <p < 0.3 0.3 <p < 0.6 0.6 <p <0.9 0.9 <p < 1.2 mttf > mttf > mttf > mttf mttf > mttf > mttf > mttf mttf > mttf > mttf > mttf mttf > mttf > mttf > mttf
3 0 <a < 0.3 0.3 <a < Q.6 A (œ) > A3 (œ) > A4 (œ) > Aj (œ) A (œ) > A3 (œ) > A4 (œ) > Aj (œ)
0.6 <a0 < 0.9 0.9 <a0 < 1.2 A (œ) > A (œ) > A (œ) > A (œ) A (œ) > A3 (œ) > A4 (œ) > A1 (œ) ßo = 0.01 c0 =100000 C = 500 c = 1000
4 0 <a0< 0.3 0.3 <a0 < 0.6 0.6 <a0 < 0.9 0.9 <a0 < 1.2 mttf > mttf > mttf = mttf mttf > mttf > mttf = mttf mttf > mttf > mttf = mttf mttf > mttf > mttf = mttf
Profit Comparison
In this section, numerically comparison with respect to the profit functions for all configurations are discussed. For consistency, we fix the following set of parameters values throughout the simulations: a0= 0.6, p0 = 0.01, k0 = 100,000, k = 500 and k2 = 1,000
100000
80000
£ 60000 CD
£ 40000 20000 0
0,1 0,2 0,3
0,4 0,5 0,6 0,7
a
0,8 0,9
I Configuration I «Configuration II «Configuration III «Configuration IV
Figure 5: Profit Comparison for all configuration using a0
100000
80000
60000
o
Q. 40000
20000
0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
A
I Configuration I «Configuration II «Configuration III «Configuration IV
Figure 6: Profit Comparison for all configuration using ß
Ibrahim Yusuf, Abdullahi Sanusi
OPTIMAL SYSTEM FOR FIVE UNITS SERIAL SYSTEMS RT&A, № 2 (62) UNDER PARTIAL AND COMPLETE FAILURE_Volume 16, June 2021
Figures 5 and 6 depict the trends of profit for all the system configurations against the repair and
failure rates a0 and /30 respectively. In both figures, it is seen that as repair rate a0 increases, the
profit increases, while with increase in failure rate P0, the profit decreases. This means that preventive and major maintenance is significant in maximizing the system profit. It is also evident from these Figures that Configuration II has the highest profit as compared to the other three configurations.
Cost Benefit Comparison
In this section, the system configurations are compared based on their cost benefit, where the benefit is either availability or mean time to failure. Numerical values of Wang el al. (2006) parameter values are used to compare the configurations.
C = 48,000,000, c2 = 39,000,000, C3 = 42,000,000 and C4 = 39,000,000
I Configuration I «Configuration II «Configuration III «Configuration IV
Figure 7: C / A Comparison for all configuration using f3Q
80000000
60000000
c / mttfk
40000000 20000000 0
0 0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
A
■ Configuration I «Configuration II «Configuration III «Configuration IV
Figure 8: C / mttfk Comparison for all configuration using f3{)
250000000 200000000 ck / Ak150000000 100000000 50000000 0
I h lUUHLlLiLiLiLi
PVPPPTTTr
0,1 0,2 0,3 0,4 0,5 0,6 0,7 0,8 0,9
a0
■ Configuration I «Configuration II «Configuration III «Configuration IV Figure 9: / A Comparison for all configuration using <XQ
I Configuration I «Configuration II «Configuration III «Configuration IV
Figure 10: Ck / MTTFk Comparison for all configuration using <XQ
Figures 7-10 present the trends of cost benefit / a and C / MTTFk for all configurations against the repair and failure rates a0 and /30 respectively. It is observed from Figures 7 and 8 that ck / a and C / MTTFk increases as /30 increases for any system configuration. It is also observed from these figures that the optimal system configuration is configuration II. On the other hand, Figures 9 and 10 display the effects of / a and C / MTTFk for all the system configurations against the repair rate a0. These figures revealed that / a and C / MTTFk decreases as a0 increases for any system configuration. Also, from these figures, it is clear that the optimal configuration is configuration II.
V. CONCLUSION
In this paper, we have constructed four different standby serial systems each consisting of five units. The expressions for the system characteristics such as system availability, busy period of repairman due to partial and complete failure as well as profit functions for all the configurations have been obtained and validated by performing numerical experiments. Analysis of the effect of various system parameters on profit function and availability was performed. These are the main contributions of this study. On the basis of the numerical results obtained in Figures 5 - 10 and Tables 1-4 for a particular case, it is evident that the optimal system configuration is configuration II. This is supported from analytical comparison presented in terms of the availability and mean time to failure models obtained in which configuration II is the optimal configuration for all a0, P > 0
contrary to some studies where the optimality among the system configuration is not uniform as it depends on some system parameters.
References
[1] Ashok Kumar, D. Pawar1 & S. C. Malik. (2020). Reliability Analysis of a Redundant System with 'FCFS' Repair Policy Subject to Weather Conditions. International Journal of Advanced Science and Technology, 29(3),7568 - 7578.
[2] Berk, M., Schubert, O and Kroll, H,-M.. Buschardt, B and Straub, D. (2019). Reliability assessment of safety-critical sensor information. IEEE Transactions on Reliability. 68(4), 1227 -1241.
[3] EL-Sherbeny, M.S. (2010). The optimal system for series systems with mixed standby components. Journal of Quality in Maintenance Engineering, Vol. 16 No. 3, pp. 319-334.
[4] Garg, H. (2018). Multi objective non-linear programming problem for reliability optimization in intuitionistic fuzzy environment. Frontiers in Information Systems, Vol. 2, pp. 197-229.
[5] Gahlot, M., Singh, V.V., Ayagi, H. I. and Goel, C.K. (2018). Performance assessment of repairable system in series configuration under different types of failure and repair policies using copula linguistics. International Journal of Reliability and Safety, Vol.12 No.4, pp.348 - 363. DOI: 10.1504/IJRS.2018.096059.
[6] Ibrahim, K.H., Singh, V.V. and Lado, A. (2017). Reliability assessment of complex system consisting two subsystems connected in series configuration using Gumbel-Hougaard family copula distribution. Journal of Applied Mathematics and Bioinformatics, Vol. 7 No. 2, pp. 1-27.
[7] Kakkar, M.K., Bhardwaj, N., Sharma, N and Bhatti, J. (2017). Availability analysis of two parallel unit system under the provision of maintenance. International Journal of Recent Research Aspects, 4(2), 155-160.
[8] Lado A., and Singh, V. V. (2019). Cost assessment of complex repairable systems in series configuration using Gumbel Hougaard family copula. International Journal of Quality and Reliability Management, 36(10), pp. 1683-1698.
[9] Mortazavi, S.M., Mohamadi, M., Jouzdani, J. (2018). MTBF evaluation for 2-out-of-3 redundant repairable systems with a common cause and cascade failures considering fuzzy rates for failures and repair: a case study of a centrifugal water pumping system. Journal of Industrial Engineering Int, 14:281. https://doi.org/10.1007/s40092-017-0226-6.
[10] Malik, S., Tewari, P.C. (2018). Performance modeling and maintenance priorities decision for water flow system of a coal-based thermal power plant. International Journal of Quality and Reliability Management, Vol. 35(4), pp.996-1010.
[11] Niwas, R. and Garg, H. (2018). An approach for analyzing the reliability and profit of an industrial system based on the cost-free warranty policy. Journal of the Brazilian Society of Mechanical Sciences and Engineering, 40: 265.
Ibrahim Yusuf, Abdullahi Sanusi
OPTIMAL SYSTEM FOR FIVE UNITS SERIAL SYSTEMS RT&A, № 2 (62) UNDER PARTIAL AND COMPLETE FAILURE_Volume 16, June 2021
[12] Sanusi, A., Yusuf, I., Mamuda. B.Y. (2020). Performance evaluation of an industrial configured as series-parallel system. Journal of Mathematical and Computational Science, 10, pp. 692-712.
[13] Sanusi, A. and Yusuf, I. (2020). Cost analysis of 2-out-of-4 system connected to two-unit parallel supporting device for operation. Life Cycle Reliability and Safety Engineering, https://doi.org/10.1007/s41872-020-00134-z.
[14] Saini, M., Kumar, A. & Shankar, V.G. (2020). A study of microprocessor systems using RAMD approach. Life Cycle Reliab Saf Eng, https://doi.org/10.1007/s41872-020-00114-3.
[15] Singh, V.V. and Ayagi, H.I. (2018). Stochastic analysis of a complex system under preemptive resume repair policy using Gumbel-Hougaard family of copula. International Journal of Mathematics in Operational Research, 12(2). 273-292.
[16] Wang, J., Xie, N., and Yang. N. (2019). Reliability analysis of a two dissimilar unit warm standby repairable system with priority in use. Commun Stat Theory Methods, https://doi.org/10.1080/03610926.2019.1642488.
[17] Wang, K.-H. and Pearn, W. L. (2003). Cost benefit analysis of series systems with warm standby components. Mathematical Methods of Operations Research, 58, 247-258.
[18] Wang, K., Hsieh, C and Liou, C. (2006). Cost benefit analysis of series systems with cold standby components and a repairable service station. Journal of quality technology and quantitative management, 3(1), 77-92.
[19] Yang, L., Guo, Y and Wang, Q. (2019).Reliability assessment of a hierarchical system subjected to inconsistent priors and multilevel data. IEEE Transactions on Reliability, 68(4), 277-292.
[20] Yusuf, I., Yusuf, B., Babagana, M., Sani, B and Lawan, M.A. (2018). Some reliability characteristics of a linear consecutive 2-out-of-4system connected to 2-out-of-4 supporting device for operation. Int J Eng. Technol. 7(1):135-139.
[21] Yen, T, -S and Wang, K. -H (2018). Cost benefit analysis of three systems with imperfect coverage and standby switching failures. Int. Journal of Mathematics in Operational Research, 12 (2), 253-272.
[22] Yusuf, I., Yusuf, B., Babagana, M., Sani, B and Lawan, M.A. (2018). Some reliability characteristics of a linear consecutive 2-out-of-4system connected to 2-out-of-4 supporting device for operation. Int J Eng. Technol, 7(1):135-139.
[23] Zhang, F. (2019). Research on reliability analysis of computer networks based on intelligent cloud computing method. International Journal of Computers and Applications, 41(4), 283-288.
[24] Zhao, B., Xie, L., Li, H., Zhang, S., Wang, B., and Li, C. (2020). Reliability analysis of aero-engine compressor rotor system considering cruise characteristics. IEEE Trans Reliab, 68(4):245-259.
