Mathematical Modeling for Performance Analysis and Inference of k-out of-n Repairable System Integrating Human Error and System Failure
1Department of Mathematics, School of Science & Technology, University of Fiji, Saweni, Fiji2Vision Institute of Technology Aligarh, G.B. Technical University, India
3. Notations and Description of the System
4. Evaluation of State Probabilities of the System
5. Solution of the Mathematical Model
Abstract
Present paper demonstrates mathematical modeling and evaluation of performance measures of k-out of-n repairable system wherein the most influencing constraints of human error and common-cause failure have been taken into consideration. Firstly the mathematical modeling is developed for performance analysis of k-out of-n repairable system with standby units involving human and common-cause failure. Then, a successful attempt has been made to evaluate various important performance measures such as availability of system, steady state availability and mean time of system failure (MTSF), mean operational time(MOT), expected busy period (EBP) and steady state busy period etc. Using the supplementary variable technique, Laplace transforms of various state probabilities are explored. Moreover, a particular case when repair rate follows exponential distribution has also been discussed. In addition, numerical illustration has also been presented in order to enable a better mode for understanding and testing the outcomes explored herein. Finally, tables and graphs for investigated results are displayed for drawing some significant conclusive observations for testing their validity and consistency.
At a glance: Figures
Keywords: Mathematical modeling, repairable system, human error, common-cause failure, steady state availability, Laplace transform, supplementary variable technique, performance analysis
American Journal of Modeling and Optimization, 2014 2 (1),
pp 16-24.
DOI: 10.12691/ajmo-2-1-3
Received January 20, 2014; Revised February 12, 2014; Accepted February 27, 2014
Copyright © 2013 Science and Education Publishing. All Rights Reserved.Cite this article:
- Maurya, Vishwa Nath. "Mathematical Modeling for Performance Analysis and Inference of k-out of-n Repairable System Integrating Human Error and System Failure." American Journal of Modeling and Optimization 2.1 (2014): 16-24.
- Maurya, V. N. (2014). Mathematical Modeling for Performance Analysis and Inference of k-out of-n Repairable System Integrating Human Error and System Failure. American Journal of Modeling and Optimization, 2(1), 16-24.
- Maurya, Vishwa Nath. "Mathematical Modeling for Performance Analysis and Inference of k-out of-n Repairable System Integrating Human Error and System Failure." American Journal of Modeling and Optimization 2, no. 1 (2014): 16-24.
Import into BibTeX | Import into EndNote | Import into RefMan | Import into RefWorks |
1. Introduction
Literature shows that performance analysis of networks with queues has received considerable attention by a large number of previous noteworthy researchers and it has occupied a prominent place in Operational Research in recent years. Repairable system incorporating human error and system failure is a particular type of queueing system. In series systems, the failures of one or more units result in the system failure. However, there exist systems that are not considered failed until at least k units or components have failed. Such systems are known as k-out of-n systems. Examples of such systems are: large airplanes usually have three or four engines, but two engines may be the minimum number required to provide a safe journey. Similarly, in many power-generation systems that have two or more generators, one generator may be sufficient to provide the power requirements. Performance analysis of wide range of queueing models in different frame works has been confined by several noteworthy researchers, (e.g. [3, 4, 5, 6, 7, 8, 11, 13]) and references therein. Among them, recently Maurya [7, 8] paid keen attention concerning with performance analysis of MX/(G1,G2)/1 queueing model with second phase optional service and Bernoulli vacation schedule and succeeded to investigate significant performance measures.
It is relevant to mention here that mathematical modeling pay a key role for performance analysis of system repair problems like other mathematical techniques. In this direction, some previous researchers (e.g. [3, 9, 10, 16]) and references therein are worth mentioning. Although, reliability and/or performance analysis of system repair problems have been attempted by some previous researchers, (e.g. [2, 14, 15, 16, 17]) but human error and common-cause failure have not been yet taken into consideration in their studies. However, it has been keenly observed that human error and common-cause failure are very much influencing constraints in performance and/or reliability analysis of system repair problems. Few earlier research workers confined the reliability behavior of k-out of-n systems taking into account the common-cause failure, for example, Hughes (1987) presented a new approach to common-cause failure in analyzing system repair problems and Who Kee Chang (1987) attempted for reliability analysis of a repairable parallel system with standby involving human failure and common-cause failures yet no attention has been paid to the performance evaluation of k-out of-n repairable system due to human error and system failure constraints. Jain et al (2002) studied a k-out of-n system with dependent failure and standby support and Moustafa (1997) also studied k-out of-n repairable system with dependent failure and imperfect coverage.
In this paper, our main objective is to investigate performance measures and valuable inferences of k-out of-n repairable system with standby units involving human and common-cause failure by using mathematical modeling, particularly to evaluate availability, steady state availability, mean time of system failure (MTSF) etc. The supplementary variable technique has been used here to evaluate Laplace transforms of various state probabilities of the repair model. Numerical computations and graphs for availability and MTSF have also been drawn. Following three particular cases are also developed:
(i). The case when there is repair due to hardware, human and common-cause failures
(ii). The case when there is repair due to hardware and human failures
(iii). The case with no repair
2. Hypotheses of the System
Here, we consider k-out of-n repairable system with standby units involving human and common-cause failure. Following are some assumptions taken into our present consideration for the repair model.
(i). The system consists of n identical main units and s standby units
(ii). Unit failure, human error and common-cause failure are constant
(iii). Repair rate from failed states due to unit failure, human error and common-cause failure are generally distributed
(iv). Initially n units are operating and s unit are kept as cold standby
(v). The entire system working if at least k out of n units or components are operating
(vi). The system is said to be in one of the failed if k+1 unit have failed due to unit failure, human error and common-cause failure
(vii). When any of the operating units fails, it is replaced immediately by standby unit
(viii). If all the standbys are consumed, the system works as degraded system until k-units works
(ix). No repair will be undertaken until the system has failed due to hardware and human error (i.e. until k+1 unit have failed) or due to common-cause failure
(x). We assume that a repaired system is as good as new
(xi). A perfect switch is used to switch-on the standby units and switch-over time is negligible
3. Notations and Description of the System
In order to analyze our present repair model in the direction as indicated earlier, we use underlying notations;
Pi,j(t) Probability that the system is in state (i, j) at time t. State (i, j) is the state of the system when i units failed due to hardware failure and j units due to human error, i, j = 0, 1, 2, …, k+1. State (0, 0) is the initial state at t = 0 and states (i, k+ 1 –i, i = 0, 1,……., k +1) are the failed states of the system. State cc is when the system has failed due to common-cause failure.
λ constant hardware failure rate of a unit.
h constant human error rate of a unit.
λc common-cause failure rate
Λ= n λ, H= n h
ρ repair rate when the system is in state (i, j)
ρc repair rate when the system is in state cc
Pi,j(y,t) Probability density function that the failed system is in state (i, j) and has an elapsed repair time of y at time t; j = k +1 –i, i= 0, 1, 2, ……….k +1
Pcc(y,t) Probability density function that the failed system is in state cc and has an elapsed repair time of y at time t.
s. Laplace transform variable, Laplace transform of the probability density function f(.)
4. Evaluation of State Probabilities of the System
In view of our assumptions and notations as in sections 2-3, the set of governing differential equations for the model can be expressed as following;
(4.1) |
(4.2) |
(4.3) |
(4.4) |
(4.5) |
(4.6) |
(4.7) |
(4.8) |
5. Solution of the Mathematical Model
By taking the Laplace transform of equations (4.1)-(4.7), one may get a set of following integral and differential equations:
(5.1) |
(5.2) |
(5.3) |
(5.4) |
(5.5) |
(5.6) |
(5.7) |
By integrating equations (5.3), (5.4), we get:
(5.8) |
(5.9) |
Again integrating equations (5.8), (5.9), using equations (5.5)-(5.7), we get
(5.10) |
(5.11) |
Using equations (5.5)-(5.7), we have following equations also from equations (5.8)-(5.9);
(5.12) |
(5.13) |
Lastly by substituting equations (5.11), (5.12) in equation (5.1), we can easily get
The Laplace transform of the reliability of the system and steady state availability of the system are:
6. Particular Case
When repair rate follows exponential distribution. Setting
(6.1) |
The Laplace transform of the probabilities of the system are:
(6.2) |
(6.3) |
(6.4) |
(6.5) |
Consider the case when n=2, s=1, k=1, , are the repair rates when system is in state (i, j), and in cc states respectively. In this case we have:
Moreover, from equations (6.2)-(6.5), one may get
(6.1.1) |
(6.1.2) |
(6.1.3) |
(6.1.4) |
(6.1.5) |
(6.1.6) |
(6.1.7) |
where :
Therefore the Laplace transform of probability of the system is in up-state is:
(6.1.8) |
where
Taking inverse Laplace transform of equation (6.1.8), one may get:
(6.1.9) |
where
and are the roots of the equation
6.1.1. Mean Operational Time (MOT)
From equation (6.1.9) one may get:
(6.1.10) |
6.1.2. Busy Period
(6.1.11) |
Taking inverse Laplace transform of equation (6.1.11), it is fairly easy to get busy period:
(6.1.12) |
where
6.1.3. Expected Busy Period (EBP)
From equation (6.1.12) one may get the expected busy period as following;
(6.1.13) |
6.1.4. Steady State Availability of the System
(6.1.14) |
where:
6.1.5. Steady State Busy Period
(6.1.15) |
6.1.6. Mean Time of System Failure (MTSF)
(6.1.16) |
where :
We consider the case when n=2, s=1, k=1, are the repair occurs when system is in state (i, j). In this case, we have: and from equations (6.2)-(6.5), one may get
(6.2.1) |
(6.2.2) |
(6.2.3) |
(6.2.4) |
(6.2.5) |
(6.2.6) |
(6.2.7) |
Therefore the Laplace transform of probability of the system is in up-state is:
(6.2.8) |
where
Taking inverse Laplace transform of equation (6.1.8), one may get after a little simplification:
(6.2.9) |
where
and are the roots of the equation
6.2.1. Expected Operational Time
From equation (6.2.9), one may get the expected operational time as below;
(6.2.10) |
6.2.2. Busy Period
(6.2.11) |
where
Taking inverse Laplace transform of equation (6.2.11) and after a little simplification, one may get:
(6.2.12) |
where
6.2.3. Expected Busy Period
From equation (6.2.12) , one may get
(6.2.13) |
6.2.4. Steady State Availability of the System
(6.2.14) |
6.2.5. Steady State Busy Period
(6.2.15) |
6.2.6. Mean Time of System Failure (MTSF)
(6.2.16) |
Consider the case when n=2, s=1, k=1, and no repair, we have:
And from equations (6.2)-(6.5), one may get
(6.3.1) |
(6.3.2) |
(6.3.3) |
(6.3.4) |
(6.3.5) |
(6.3.6) |
(6.3.7) |
Therefore the Laplace transform of probability of the system is in upstate is:
(6.3.8) |
where
Taking inverse Laplace transform of equation (6.1.8) one may get:
(6.3.9) |
where
and are the roots of the equation
6.3.1.Mean Operational Time (MOT)
From equation (6.3.9), one may get:
(6.3.10) |
6.3.2. Steady State Availability of the System
(6.3.11) |
6.3.3. Mean Time of System Failure (MTSF)
(6.3.12) |
7. Numerical Illustration
In this section, our endeavor is to test the validity of various performance measures explored herein by way of numerical illustration. Setting , one can compute results case wise as following;
7.1. Case Iwhen n=2, s=1, k=1, , are the repair rates when system is in state (i, j), and in cc states respectively
7.1.1. Availability Analysis
(7.1.1) |
7.1.2. Expected Operational Time
(7.1.2) |
7.1.3. Busy Period
(7.1.3) |
7.1.4. Expected Busy Period
(7.1.4) |
when n=2, s=1, k=1, are the repair occurs when system is in state (i, j)
7.2.1. Availability Analysis
(7.2.1) |
(7.2.2) |
7.2.2. Busy Period
(7.2.3) |
7.2.3. Expected Busy Period
(7.2.4) |
when n=2, s=1, k=1, and no repair
7.3.1. Availability Analysis
(7.3.1) |
(7.3.2) |
Setting , in equations (7.1.1), (7.2.1) and (7.3.1), one gets Table 1. Variation of availability w.r.t. time in three cases is shown in Figure 1.
Setting , in equations (7.1.2), (7.2.2) and (7.3.2), one gets Table 2. Variation of expected operational time w.r.t. time in three cases is shown in Figure 2.
Setting , in equations (6.1.14), (6.2.14), one gets Table 3. Variation of steady state availability for different values of human error in three cases is shown in Figure 3.
Setting , in equations (6.1.16), (6.2.16) and (6.3.12), one gets Table 4. Variation of MTSF for different values of human error in three cases is shown in Figure 4.
MTSF of the system with ρ and ρc = ∞, MTSF of the system with ρ = 500
8. Sensitivity Analysis and Conclusive Observations
Table 1 computes availability of the system at any time. Figure 1 shows availability of the system decreases in the interval (0, t] for the three cases. By comparing the availability with respect to time for three cases with and without repair graphically, it is observed that:
• The availability of system decreases with respect to time. We conclude that the system availability with ρ and ρc is greater than the system with ρ, and availability with ρ is greater than without repair.
Table 2 computes expected operational time for the system at any time. Figure 2 shows expected operational time increases in the interval (0, t] for the three cases. By comparing the expected operational time with respect to time t for three cases with and without repair graphically, we observe that:
• The system expected operational time increases with respect to time. And it is also noticeable that the system expected operational time with ρ and ρc is greater than the system with ρ and with ρ repair is greater than the system without repair.
Table 3 compute the relationship between steady state availability and human error for three cases. Figure 3 shows variation of steady state availability with respect to human failure for the two cases of repair. By comparing the steady state availability with respect to human failure for the systems with ρ and ρc graphically, it is remarked that:
• The human failure rate h increases, however the steady state availability of the system decreases at constant .
• The steady state availability of system with repair only is greater than the system availability with ρ and ρc (i.e., system include common-cause failure). The availability of the system in the steady-state is equal to zero.
Table 4 computes variation of mean time of system failure (MTSF) with respect to human failure. Figure (4) displays the variation of MTSF with respect to human failure for the second and third cases. By comparing the MTSF with respect to human failure for the system with ρ and without repair graphically, we observe that:
• The mean time of system failure (MTSF) with ρ and ρc = ∞
• MTSF with ρ = 500 (constant). However, in case of MTSF without repairing, as for as the value of human failure rate h increases at constant , the MTTF of the system decreases.
Finally with passing above remarks, we have succeeded to investigate significant performance measures of k-out of-n repairable system involving human and common-cause failures. It is highly expected that our present contribution for performance analysis and inference of k-out of-n repairable system involving human and common-cause failures will be useful for researchers, statisticians, mathematicians and management professionals for their future research and development in this direction, e.g. very recently Maurya (2013 c) proposed a computational approach to cost and profit analysis of k-out of-n repairable system integrating human error and system failure constraints.
Acknowledgements
The present paper is an integral part of the post-doctoral dissertation entitled “Performance Analysis And Inference of Mixed Poisson Queueing Models” published in Scholars Press Publishing, Germany. Prof. V.N. Maurya; author of the paper would like to express his heartiest gratitude to his D.Sc. advisor Prof. (Dr.) R.B. Misra. Ex-Vice Chancellor, Dr. Ram Manohar Avadh University Faizabad, India; Hon’ble Vice Chancellor, University of Fiji, Saweni and Dean, School of Science & Technology, University of Fiji for their encouragement and sincere support to provide all essential facilities.
References
[1] | Hughes, R.P., 1987. A new approach to common-cause failure. Reliability Engineering and System Safety, Vol. 17, pp. 211-236. | ||
In article | CrossRef | ||
[2] | Jain Madhu, Sharma G.C. and Alok Kumar, 2002. K-out of-n system with dependent failure and standby support. JKAU: Engineering Science, Vol. 14, No. 2, pp. 73-83 | ||
In article | |||
[3] | Maurya Avadhesh Kumar and Maurya V.N., 2013. A novel algorithm for optimum balancing energy consumption LEACH protocol using numerical simulation technique, International Journal of Electronics Communication and Electrical Engineering, Algeria, Vol. 3, No. 4, pp. 1-19. | ||
In article | |||
[4] | Maurya Avadhesh Kumar, Maurya V.N. and Arora D. Kaur, 2013. Linear regression and coverage rate performance analysis for optimization of received signal strength in antenna beam-tilt cellular mobile environment, International Journal of Electronics Communication and Electrical Engineering, Algeria, Vol. 3, No. 7, pp. 1-14. | ||
In article | |||
[5] | Maurya V.N., 2010. On busy period of an interdependent M/M/1:(∞;GD) queueing model with bivariate Poisson process and controllable arrival rates, IEEE Transactions, pp. 243-246 | ||
In article | |||
[6] | Maurya V.N., 2012. Determination of expected busy periods in faster and slower arrival rates of an interdependent M/M/1:(∞; GD) queueing model with controllable arrival rates, International Journal of Engineering Research and Technology, Engineering Science & Research Support Academy (ESRSA) Publication, Vadodara, India, Vol.1, No. 5, pp. 1-5. | ||
In article | |||
[7] | Maurya V. N., 2013 a. Sensitivity analysis on significant performance measures of bulk arrival retrial queueing MX/(G1,G2)/1 model with second phase optional service and Bernoulli vacation schedule, International Open Journal of Operations Research, Academic and Scientific Publishing, New York, USA, Vol. 1, No. 1, pp. 1-15. | ||
In article | |||
[8] | Maurya V.N., 2013 b. Performance analysis of MX/(G1,G2)/1 queueing model with second phase optional service and Bernoulli vacation schedule, American Journal of Modeling and Optimization, Science and Education Publishing, New York, USA, Vol. 2, No.1, pp 01-07. | ||
In article | |||
[9] | Maurya V.N., Arora Diwinder Kaur, Maurya Avadhesh Kumar and Gautam Ram Asrey, 2013 a. Numerical simulation and design parameters in solar photovoltaic water pumping systems, American Journal of Engineering Technology, Academic & Scientific Publishing, New York, USA, Vol.1, No. 1, pp. 1-09. | ||
In article | |||
[10] | Maurya V.N., Arora Diwinder Kaur, Maurya Avadhesh Kumar and Gautam R.A., 2013 b. Exact modeling of annual maximum rainfall with Gumbel and Frechet distributions using parameter estimation techniques, World of Sciences Journal, Engineers Press Publishing, Vienna, Austria, Vol. 1, No. 2, pp.11-26. | ||
In article | |||
[11] | Maurya V.N., 2000. A study of use of stochastic processes in some queueing models, Ph.D. Thesis, Department of Mathematics & Statistics, Dr. R.M.L. Avadh University, Faizabad, India. | ||
In article | |||
[12] | Maurya V.N., 2013 c. Computational approach to cost and profit analysis of k-out of-n repairable system integrating human error and system failure constraints, Physical Sciences Research International, Net Journals, Abuja, Nigeria, Vol. 1, No. 4, pp. 133-140. | ||
In article | |||
[13] | Mishra S.S. and Yadav D.K., 2010. Computational approach to cost and profit analysis of clocked queueing networks, Contemporary Engineering Sciences, Vol.3, No. 8, pp. 365-370. | ||
In article | |||
[14] | Moustafa S.M., 1997. Reliability analysis of k-out of-n systems with dependent failures and imperfect coverage. Reliability Engineering and System Safety, Vol. 58, pp. 15-17. | ||
In article | CrossRef | ||
[15] | Pham H., Suprasad A., and Misra R.B., 1997. Availability and mean life time predication of multistage degraded system with partial repair, Reliability Engineering and System Safety, Vol. 56, pp. 169-173. | ||
In article | CrossRef | ||
[16] | Shao J. and Lamberson L.R., 1991. Modeling a shared-load k-out of-n: G systems, IEEE Trans. Reliab. Vol. 40, No. 2, pp. 202-208 | ||
In article | CrossRef | ||
[17] | Wang K.H. and Ming Y., 1997. Profit analysis of M/EK/1 machine repair problem with a non-reliable service station, Computers and Industrial Engineering, Vol. 32, pp. 587-594. | ||
In article | CrossRef | ||
[18] | Who Kee Chang, 1987. Reliability analysis of a repairable parallel system with standby involving human failure and common-cause failures, Micro-electron. Reliab. Vol. 27, No 2, pp. 269-271. | ||
In article | CrossRef | ||