Life Behavior of a System under Discrete Shock Model
Serkan Eryilmaz1
Academic Editor: M. De la Sen
Received21 Jun 2012
Revised31 Jul 2012
Accepted06 Aug 2012
Published11 Sept 2012
Abstract
We study the life behavior of a system which is subjected
to shocks of random magnitudes over discrete time periods. We obtain the
survival function and mean time to failure of the system assuming that the
sizes of the shocks follow a discrete probability distribution under cumulative
and mixed shock models.
1. Introduction
There are various engineering systems which are subjected to shocks of random magnitudes at random times. The shock models can be classified in different ways. According to the cumulative shock model, the system breaks down because of a cumulative effect of shocks, while in an extreme shock model the system fails because of one single shock with large magnitude. See, for example, [1β9] for various problems on shock models.
Most of the studies on shock models focus on the evaluation of system failure time in a continuous setup, that is, the shocks arrive according to a renewal process, and the times between successive shocks have a continuous probability distribution. Some results on discrete case are in [3, 7, 10].
Consider a system which is subjected to periodic random shocks. A shock occurs with probability in each period . The period should be understood as hour, day, and so forth. The magnitude of the shock which occurs in period is a random variable denoted by . Assume that such a system fails if and only if the sum of the magnitudes of cumulative shocks exceed, the level for . Let be a binary random variable representing the shock occurrences that is, if a shock occurs in period and , otherwise. For , define
where the random variables and are independent in each time period. The random variable is strictly positive, and is a sequence of independent identically distributed (i.i.d.) random variables with cumulative distribution function (c.d.f.) and probability mass function (p.m.f.) .
Thus, under the cumulative shock model, the failure time of the system can be defined by the following waiting time random variable:
for .
In the case of a mixed shock model, a system fails if either a single shock with a large magnitude occurs or the sum of cumulative shocks exceeds the critical level. Thus, in this case the time to failure of the system is defined by the following compound waiting time random variable.
where
where for .
Such models can also be applied to insurance, replacing shock with claim and magnitude of the shock with claim amount. In this case, a period can be seen as a week, month, and so forth, and the random variable represents the waiting time until the cumulative sum of claim amounts exceeds the level . Similarly, the random variable is the waiting time until the first extreme claim size falls above the level .
The present paper is organized as follows. In Section 2, we derive recurrence formulae for the survival function and the mean time to failure (MTTF) of the system under the cumulative shock model. We also study two related characteristics and which represent, respectively, the number of shocks and the total shock that the system is subjected up to time when the system fails. Section 3 includes the results for mixed shock model.
2. Cumulative Shock Model
In the following, we derive two popular reliability characteristics: survival function and mean time to failure of the system under the cumulative shock model.
It is clear that
where , and denotes the -fold convolution of with itself, . By conditioning on the claim occurrence, one obtains
where .
Theorem 2.1. For ,
and , where is the endpoint of the support of .
Proof. From (2.1), it follows that
Thus, the proof is immediate from (2.2).
Proposition 2.2. For , the MTTF of the system can be computed from
with .
Proof. Using (2.1),
Thus, the proof follows from (2.2) since
Example 2.3. Let have a geometric distribution with pmf . Then under the conditions of Proposition 2.2,
with .
2.1. Related Characteristics
For , we define new random variables as follows:
It is clear that the random variables and represent, respectively, the number of shocks and the total shock that the system is subjected up to time when the system fails. These two characteristics might be useful for improvement purposes and can be effectively used in optimal system design.
Theorem 2.4. For ,
where , and and can be computed recursively from
for and for , and
for and for .
Proof. By conditioning on ,
The probability can be written as follows:
Thus, we need to get recurrences for and . By conditioning on the values of ,
On the other hand,
for . Thus, the proof is completed.
Before proceeding with the distribution of , it should be noted that the random variable denotes the total shock up to time and
for and .
Theorem 2.5. For ,
where
Proof. By the definition of ,
For ,
The proof follows by conditioning on and noting that for .
The following result readily follows from the definitions of and and Waldβs equation.
Proposition 2.6. For ,
In Table 1 we compute , , and whenever the shock size random variable has a geometric distribution with mean . From Table 1 we observe that an increase in leads to an increase in MTTF of the system. If the probability of observing a shock in a period increases, then the MTTF decreases. We also observe that MTTF is proportional to . Therefore, for the same shock size distribution the expected number of shocks and expected total shock remain the same for different values of .
3. Mixed Shock Model
For , the mixed shock model is same as the cumulative shock model. Thus we assume that . The following is a recursive equation for the survival probability of the system under mixed shock model.
Theorem 3.1. For and ,
and , where if .
Proof. For ,
By conditioning on the values of ,
By conditioning on ,
Thus, the proof is completed.
The following result can be proved similar to Proposition 2.2, and hence its proof is omitted.
Proposition 3.2. For , the MTTF of the system under mixed shock model can be computed from
with , where if .
In Table 2, using Proposition 3.2, we compute the MTTF of the system under mixed shock model when the shock size random variable has a geometric distribution with mean .
Theorem 3.3. For ,
where
Proof. By conditioning on ,
It is clear that
By the definition of ,
On the other hand,
Thus the proof is completed.
Before the derivation of the distribution of , we note the following recursion which will be useful in the sequel:
Theorem 3.4. For ,
for ,
Proof. By the definition of ,
For ,
Similarly, for ,
Thus, the proof is completed.
4. Summary and Conclusions
In this paper, we studied the life behavior of a system under discrete time cumulative and mixed shock models. The probability of getting a shock in any period is , and the shock occurrences are assumed to be independent over the periods. The size of the shock occuring in a period follows a discrete probability distribution and the systemβs lifetime coincides with the waiting time random variable which represents the time until the cumulative sum of shocks exceeds a specified level (cumulative shock model). We derived recurrence formulae for the survival function and the MTTF of the system. We also obtained recurrences for the distributions and expected values of the two related quantities which represent the number of shocks and the total shock that the system is subjected until failure. The results were illustrated for the case when the shock size distribution is geometric. We have also obtained a recurrence for the survival function of the system under a mixed shock model. The assumption of discrete shock size distribution enables us to obtain recursive formulae. However, the consideration of continuous shock size distribution might be of special interest in some applications. Therefore, a possible future work can be on discrete time shock models with a continuous shock size distribution.
In the model that was studied in the paper shock occurrence indicators are assumed to be independent and identical with a constant probability . As a future work, the case in which the shock occurrence indicators form a Markov chain can also be considered.
Acknowledgment
The author thanks referees for their very useful comments and suggestions, which improved the paper.
References
J. M. Bai, Z. H. Li, and X. B. Kong, βGeneralized shock models based on a cluster point process,β IEEE Transactions on Reliability, vol. 55, no. 3, pp. 542β550, 2006.
S. Eryilmaz, βOn the lifetime behavior of discrete time shock model,β Journal of Computational and Applied Mathematics, vol. 237, no. 1, pp. 384β388, 2013.
M. Finkelstein and F. Marais, βOn terminating Poisson processes in some shock models,β Reliability Engineering and System Safety, vol. 95, pp. 874β879, 2010.
U. Sumita and J. G. Shanthikumar, βA class of correlated cumulative shock models,β Advances in Applied Probability, vol. 17, no. 2, pp. 347β366, 1985.
T. Aven and S. Gaarder, βOptimal replacement in a shock model: discrete time,β Journal of Applied Probability, vol. 24, no. 1, pp. 281β287, 1987.