Research Article  Open Access
S. B. Munoli, Rohit R. Mutkekar, "Estimation of Reliability for a Two Component Survival StressStrength Model", Journal of Quality and Reliability Engineering, vol. 2011, Article ID 721962, 8 pages, 2011. https://doi.org/10.1155/2011/721962
Estimation of Reliability for a Two Component Survival StressStrength Model
Abstract
The reliability function for a parallel system of two identical components is derived from a stressstrength model, where failure of one component increases the stress on the surviving component of the system. The Maximum Likelihood Estimators of parameters and their asymptotic distribution are obtained. Further the Maximum Likelihood Estimator and Bayes Estimator of reliability function are obtained using the data from a lifetesting experiment. Computation of estimators is illustrated through simulation study.
1. Introduction
Several authors have considered estimation of system reliability based on stressstrength models. Here are a few references of contributions towards these models. Church and Harris [1] considered estimation of reliability from stressstrength relationship. Downton [2] considered the case of estimation of reliability for a stressstrength model under normal distribution. Wani and Kabe [3] have considered the problem of estimation of system reliability where life time of each component has gamma distribution. Constantine et al. [4] considered estimation of stressstrength relationship under the assumption of stressstrength random variables following gamma distribution with known shape parameters. Bhattacharya and Johnson [5] have considered estimation of reliability in a multicomponent stressstrength model. Kunchur and Munoli [6] have considered estimation of reliability for a multicomponent survival stressstrength model based on exponential distribution. Dan and Krausz [7] have obtained inference for a multistep stressstrength model of parallel system. Kunchur and Munoli [8] have considered estimation of reliability in Freund’s model for a two component system. Hanagal [9] has considered the problem of estimation of system reliability in a two component stressstrength model with cases on distribution of stress as exponential and gamma. Kundu and Gupta [10] have obtained inference of stressstrength relationship for generalized exponential scale family of distribution. Bhattacharya [11] has proposed Bayesian approach to lifetesting and reliability estimation. Draper and Guttman [12] have obtained the Bayes estimator of reliability in a multicomponent stressstrength model.
In the present study, we are considering a system of two components. The component survives as long as the stress on it is smaller than its strength. The system survives if at least one component functions (parallel system). Here the stress and strength associated with the components of the system are random variables. To carry out the inference, we assume certain probability distribution for these random variables. Let the strength of the two components be and , where , are independently and identically distributed gamma random variables with shape parameter and scale parameter . Let , be the stress on the two components, respectively. Initially , are independently and identically distributed as exponential random variables with parameter . Statistically dependent failures are typical for modern systems, which involve complicated interactions among component parts, and it is reasonable to assume that the failure of one component does change the stress on surviving components. The strength on the other hand is dependent on various inbuilt qualities of the component such as the technology by which it has been manufactured, the raw materials used, and so forth. Failure of one component of the system need not change these inbuilt qualities of the surviving component of the system. Hence, change in the stress of the surviving component is assumed on failure of other components of the system. That is, in a two component parallel system, the distribution of stress of the surviving component changes upon the failure of the other component. Then the exponential distribution with parameter follows. The system fails whenever both the components of the system fail.
A block diagram for the proposed model is given in Figure 1.
One can quote several examples for two component parallel systems such as pair of kidneys, eyes, hands, and legs, pair of elevators, and pair of engines in an aircraft. To illustrate the function of the proposed model let us consider the example of a pair of kidneys. Here the function of the kidneys is to purify the blood and thus help to maintain the body in healthy condition. Here the pair of kidneys perform the same function as per their natural builtup mechanism (strength). Failure of one kidney increases the purification work load on the surviving kidney (stress). Here the surviving kidney should carry out the purification function the same as it was when both the kidneys were functioning. Taking this scenario under consideration, we consider that failure of one component of the system changes only the stress and not the strength of the surviving component.
The reliability function is derived in Section 2. Lifetesting experiment is explained in Section 3. The Maximum Likelihood Estimators (MLEs) are obtained in the same section. Section 4 deals with Bayes estimation of reliability function. Computations of estimators of reliability function (MLE and Bayes) along with the findings of the study are discussed in Section 5. Some results that support the findings of this research paper are proved in the appendix.
2. Reliability Function
In order to find the reliability function of the model discussed in Section 1, let us consider that is the minimum and is the maximum and let .
Here follows exponential distribution with parameter and follows exponential distribution with parameter .
The reliability function of this system is given by The details of these derivations are given in Lemma A.1.
3. LifeTesting Experiment
In order to obtain the estimators of , suppose that “” systems whose life distribution is characterized by the reliability function derived in Section 2 are put on lifetesting experiment. Here, , are observed and , are independently and identically distributed gamma random variables with shape parameter and scale parameter . Also, the data of stress and , are obtained separately from a simulation of conditions of the operating environment. , are exponential random variables with parameters and , respectively.
Now, the joint probability density function of the random variables, , , is given by where The MLEs (the estimators that maximize the likelihood function) of parameters , , , are given in the following expressions, respectively, where denote
Using the invariance property of MLEs, the MLE of reliability function is obtained by substituting the MLEs of parameters , , , and in expression (1), that is,
The details of these findings are given in the form of Lemma A.2.
From the asymptotic properties of MLEs under regularity conditions and multivariate central limit theorem, we have , where are the corresponding elements of the inverse of Fisher Information Matrix “” as and its inverse is given by
4. Bayes Estimation of Reliability Function
In order to obtain Bayes estimator of reliability function, assume that parameters and are known (i. e., and and consider the prior distribution of parameters and [13] as
The Bayes estimator of reliability function is obtained as the posterior expectation of and is given by where with with Substituting values of expressions (13), and (14) in expression (12) we obtain the Bayes estimator of reliability function.
The details of these findings are given in the form of Lemma A.3.
5. Computation of Estimators
For the system, the random variables , (with respect to strength) and random variables (with respect to stress) are generated independently as follows.
Step 1. Initialize , , , , for the 1st and 2nd components of the system. Uniform random numbers , are generated from . Further for a given value of , exponential random variables , are obtained for the 1st and the 2nd components of the system, respectively. Now, is incremented by 1 and the process of generating the exponential random variables is repeated for both the components of the system. This repetition process is continued until and the subsequent exponential random variable values generated are noted for both the components of the system. Here and for. Here denote gamma random variables with shape parameter and scale parameter. We also compute the values and .
Step 2. The whole procedure in Step 1 is repeated for number of systems, and the statistics , , are computed.
Step 3. Initialize , , , . Uniform random numbers , are generated from . Further for a given value of , exponential random variables are obtained. The value of is incremented by 1, and the above process of generated exponential random variables is repeated. This repetition process is continued until and the subsequent exponential random variable values generated , for are computed. The statistics , are computed.
Step 4. With the help of the statistics , , , , and , the MLEs of parameters , , , of the model are obtained. Using these MLEs in the expression of reliability function, the MLE of reliability function is obtained.
For the parameter values , , , and the value of reliability function is also obtained.
The Bayes estimator of reliability function is obtained using the simulated values of , , , , and for given values of and .
Tables 1, 2, and 3 give the results of the above simulation experiment for different values of and .



6. Conclusion
The reliability function for the proposed model is evaluated in terms of stressstrength relationship rather than considering the time factor, as it is realistic to observe reliability of a system functioning under the influence of external factors (stress) as compared to longevity of working time associated with the system.
Though the expression for reliability function involves sum of infinite series, it is observed that for given values of the parameters the value of the reliability function stabilizes at the 15th value of the running variable involved in the sums.
Stress and strength associated with the components of the system possess different physical properties; thereby data on stress and strength in the lifetesting experiment is observed separately based on their corresponding operative environments.
The MLEs are sufficient, efficient and also maximizes the likelihood of the joint distribution function. Further, using the invariance property of MLE, it is easy to obtain the MLE of reliability function.
Bayes estimator is based on the prior information obtained through certain pilot study that helps in synthesizing the information to be generated for the system under function. Hence, Bayes estimator of reliability function is obtained by considering certain prior information for the parameters. But the process of obtaining Bayes estimator of reliability function is quite tedious as it involves lengthy numerical calculations.
From Tables 1, 2, and 3, it is clear that for greater values of “” and “” (large sample size) both MLE and Bayes estimators perform better. In the majority of the cases both the estimators overestimate the true value of reliability function “.” Here we observe that Bayes estimator is a better estimator in terms of bias for the given data set.
Appendix
Lemma A.1. The reliability function given in expression (1) is derived as follows. The reliability function for the system under study is given by
As , its distribution is given by
Here represents the distribution function and represents probability density function of random variable Z, and as Z follows gamma distribution with shape parameter and scale parameter , one has
Using the relationship between incomplete gamma distribution and Poisson sum, one has
Substituting the above integral value in expression (A.3), one gets
Now, substituting the value of from expression (A.5) in the expression for (expression (A.1)), one will solve the integrals associated with expression (A.1) as follows:
where,
Using the results of expressions (A.8) and (A.6) in expression (A.7), one has
Substituting the results of expressions (A.6) and (A.9) in the expression (A.1), one gets the reliability function given in expression (1).
Lemma A.2. The MLEs of parameters given in expressions (4), (5), (6), and (7) and further used to determine the MLE of reliability function in expression (8) are derived as follows.
The loglikelihood function of expression (2) in is
Now,
Similarly,
Now
Solving (A.11) and (A.13) simultaneously one obtains the MLEs of parameters of and given in expression (4) and (5), respectively.
Now again,
Similarly,
Solving (A.14) and (A.15) simultaneously one obtains the MLEs of parameters of and given in expression (6) and (7), respectively.
Lemma A.3. The posterior distribution of and , that is, used in expression (12), is derived as follows.
The joint probability density function of random variables, , , , , and is given by
Integrating with respect to and over their respective range one gets, the joint probability distribution of , , , as
Dividing in expression (A.16) by in expression (A.17) one gets the posterior distribution of and as
Acknowledgment
The authors thank the reviewers for many helpful comments and suggestions on an earlier version which substantially has improved this paper.
References
 J. D. Church and B. Harris, “The estimation of reliability from stressstrength relationship,” Technometrics, vol. 12, no. 1, pp. 49–54, 1970. View at: Google Scholar
 F. Downton, “The estimation of $P(Y<X)$ in the normal case,” Technometrics, vol. 15, pp. 551–558, 1973. View at: Google Scholar
 J. K. Wani and D. G. Kabe, “Point estimation of reliability of a system comprised of k elements from the same gamma model,” Technometrics, vol. 13, pp. 859–864, 1971. View at: Google Scholar
 K. Constantine, S. Tse, and M. Karson, “Estimation of $P(Y<X)$ in the gamma distribution case,” Communications in Statistics, vol. 15, no. 2, pp. 365–388, 1986. View at: Google Scholar
 G. K. Bhattacharya and R. A. Johnson, “Estimation of reliability in a multi component stressstrength model,” Journal of the American Statistical Association, vol. 69, no. 348, pp. 966–970, 1974. View at: Google Scholar
 S. H. Kunchur and S. B. Munoli, “Estimation of reliability for a multi component survival stressstrength model based on exponential distribution,” Communications in Statistics, vol. 22, no. 3, pp. 769–779, 1993. View at: Google Scholar
 S. Dan and A. S. Krausz, “A multistep stressstrength model of a parallel system,” IEEE Transactions on Reliability, vol. R35, no. 1, pp. 119–123, 1986. View at: Google Scholar
 S. H. Kunchur and S. B. Munoli, “Estimation of reliability in Freunds model for two component system,” Communications in Statistics, vol. 23, no. 11, pp. 3273–3283, 1994. View at: Google Scholar
 D. D. Hanagal, “Estimation of system reliability in a two component stressstrength models,” Economic Quality Control, vol. 11, pp. 145–154, 1996. View at: Google Scholar
 D. Kundu and R. D. Gupta, “Estimation of $P(Y<X)$ for generalized exponential distribution,” Metrika, vol. 61, no. 3, pp. 291–308, 2005. View at: Publisher Site  Google Scholar
 S. K. Bhattacharya, “Bayesian approach to life testing and reliability estimation,” Journal of the American Statistical Association, vol. 62, no. 317, pp. 48–62, 1967. View at: Google Scholar
 N. R. Draper and I. Guttman, “Bayesian analysis of reliability in multi component stressstrength models,” Communications in Statistics, vol. A7, pp. 441–451, 1978. View at: Google Scholar
 G. E. P. Box and G. C. Tiao, Bayesian Inference in Statistical Analysis, Addison Wesley, 1973.
Copyright
Copyright © 2011 S. B. Munoli and Rohit R. Mutkekar. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.