Next Article in Journal
Improved Block-Pulse Functions for Numerical Solution of Mixed Volterra-Fredholm Integral Equations
Previous Article in Journal
Closure System and Its Semantics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Random Walk Analysis in a Reliability System under Constant Degradation and Random Shocks

by
Jewgeni H. Dshalalow
*,† and
Ryan T. White
Department of Mathematical Sciences, Florida Institute of Technology, College of Engineering and Science, Melbourne, FL 32940, USA
*
Author to whom correspondence should be addressed.
Current address: Florida Institute of Technology, 150 W. University Blvd., Melbourne, FL 32940, USA.
Axioms 2021, 10(3), 199; https://doi.org/10.3390/axioms10030199
Submission received: 4 July 2021 / Revised: 11 August 2021 / Accepted: 13 August 2021 / Published: 23 August 2021

Abstract

:
In this paper, we study a reliability system subject to occasional random shocks hitting an underlying device in accordance with a general marked point process with position dependent marking. In addition, the system ages according to a linear path that eventually fails even without any external shocks that accelerate the total failure. The approach for obtaining the distribution of the failure time falls into the area of random walk analysis. The results obtained are in closed form. A special case of a marked Poisson process with exponentially distributed marks is discussed that supports our claim of analytical tractability. The example is further confirmed by simulation. We also provide a classification of the literature pertaining to various reliability systems with degradation and shocks.

1. Introduction

1.1. Background

Most infrastructures, industrial systems, electronic devices, and manufactured products, experience multiple failures in their life times. Not only do they all deteriorate (as an inevitable wear) but they are also subject to random shocks of various natures. Through deterioration alone, sooner or later they die. The shocks further accelerate the decay, while regular maintenance prolongs their lives. A typical modeling of a basic system includes a formalism of deterioration, also referred to as degradation or fatigue, and the presence of shocks which are almost exclusively marked Poisson. A threshold is also set so that if underlying system’s operational capacity crosses that threshold when going downhill, the system fails, and a key target is the prediction of the failure time.
The literature on reliability with degradation and shocks is very rich with myriads of very interesting and practical models. To the best of our knowledge, they establish formulas (most often of total probability type) for the distribution of the failure time that can be computed numerically. The degradation process is commonly with stationary and independent increments, such as the gamma process (almost exclusively used reliability), reverse Gaussian process, and even Brownian motion with drift (called there a Wiener process).
The combination of a continuous (internal) deterioration, say Y ( t ) , with abrupt (external) shocks cause a so-called soft failure of the system, meaning that the cumulative wear and a series of instantaneous damages exceed the sustainability of the system and make it non-operational at some point of time. The sustainability threshold M of the system will be crossed or exceeded dependent on whether it happens by gradual wear or upon one of the shocks occurring at t 1 , t 2 , . These shocks, say X 1 , X 2 , , are accumulated over the time, so that at one point in time t ν , the sum k = 1 ν X k exceeds M (unless it happens earlier upon continuous crossing of M).
In addition, the same shocks that accelerate the wear of the system can knock down one of the internal components (such as a capacitor) and make the system instantaneously non-operational. This phenomenon is referred to as a hard failure of the system. A hard failure can take place upon t 1 , t 2 , and it is specified by another control threshold D and a sequence of single damages such as W 1 , W 2 , at t 1 , t 2 , , so that a single one of them, say W ν exceeds D at time t ν . With the presence of such additional damages, the global failure of the system occurs if at some point t ν when k = 1 ν X k exceeds M or W ν exceeds D or Y ( t ) + k = 1 ν 1 X k = M for sup { t : t < t ν } , whichever of the three comes first.

1.2. Pertinent Literature

Among most pertinent work related to our reliability systems are those with shocks and degradation. There are models entirely based on degradation (as in Tsai et al. [1]) and on shocks and both.

1.2.1. Shock Models

Most shock models fall into one of the five classes: cumulative shock models, extreme shock models, δ -shock models, run shock models, and mixed shock models. A mixed shock model must be a combination of at least two of the first four types.

1.2.2. Cumulative Shock Models

A cumulative shock model can be specified by a marked point process S = k = 1 X n ε t n , where the marks X n represent magnitude of shocks. S generally uses position dependent marking. The random index ν x   = inf n : S n = k = 1 n X k > x gives the total number of shocks until the failure. The failure time is t ν x . Gut [2] showed that t ν x / x E δ E X a.s. and in L 1 where δ t n t n 1 and X X n . See also the work of Gut and Hüsler [3] and Sumita and Shanthikumar [4].
Abolnikov and Dshalalow [5] studied marked point process S = n = 0 X n ε t n with position dependent marking (that is, when X n depends on δ n = t n t n 1 , but not on any other components) we assume a delay at X 0 , t 0 . More specifically, the sequence X 0 , t 0 , X 1 , t 1 , is a delayed renewal process. Additionally, X n N 0 .
With A = [ 0 , M ) , assuming M N , of interest is the time and position of S upon its escape from A , that is the failure time and the total damage to the system upon the failure. Thus, we have: ν = inf m : S m = X 0 + + X m A c giving the number of shocks that cause the system to fail, t ν —the failure time, S ν —the total amount of damage to the system on t ν (or excess value of M).
The transform
Φ ν = Φ ν ( ξ , u , v , ϑ , θ ) = E ξ ν u S ν 1 v S ν e ϑ t ν 1 θ t ν ,
with ξ , u , v B ¯ ( 0 , 1 ) , Re ( ϑ ) 0 , and Re ( θ ) 0 , is the joint probability-generating function (PGF) of ν and S ν and Laplace Stieltjes Transform of t ν , and includes two more useful pre-failure quantities S ν 1 and t ν 1 (pre-failure time), representing the damage to the system and time the total damage seen in set A before the failure. Φ ν was expressed in a closed form through the means of the D -operator, defined as follows
D x k φ ( x , y ) = lim x 0 1 k ! k x k 1 1 x φ ( x , y ) , k 0 0 , k < 0
for x B ¯ ( 0 , 1 ) and Re ( y ) 0 . Here, φ is a function, analytic at zero in the first variable. Suppose the joint transforms
γ 0 z , θ = E z X 0 e t 0 θ
γ z , θ = E z X 1 e t 1 t 0 θ ,
assuming z B ¯ ( 0 , 1 ) and Re ( θ ) 0 , are known.
Theorem 1.
The following formula holds:
Φ ν = D x M 1 γ 0 ( v , θ ) γ 0 ( x v , θ ) + ξ γ 0 ( x u v , ϑ + θ ) 1 ξ γ ( x u v , ϑ + θ ) γ ( v , θ ) γ ( x v , θ ,
where x B ( 0 , 1 ) (the open unit ball in C ) and the rest of domains are specified as in (1).
In the work of Dshalalow and Huang [6] and Dshalalow and Robinson [7], the authors studied the above marked point process S = n = 0 X n ε t n with position dependent marking, but now with X n real-valued. With  A = [ 0 , L ) , assuming L > 0 , define
μ = inf m : S m = X 0 + + X m A C
as the total number of shocks that knock down the system, t μ —the failure time, S μ , the total damage to the system the failure time at t μ (or excess value of L). Additionally, t μ 1 is the time of the μ 1 st shock preceding the failure time and S μ 1 is the total damage to the system at time t μ 1 .
The transform
Φ μ = Φ μ ( ξ , u , v , ϑ , θ ) = E ξ μ e u S μ 1 e v S μ e ϑ t μ 1 θ t μ ,
assuming ξ B ¯ ( 0 , 1 ) , Re ( u ) 0 , Re ( v ) 0 , Re ( ϑ ) 0 , and Re ( θ ) 0 , is the joint PGF of μ and LST of S μ , t μ , and two more useful pre-failure quantities S ν 1 and t ν 1 .
Theorem 2
(Dshalalow and Huang [6]). The functional Φ μ satisfies the formula
Φ μ ξ , u , v , ϑ , θ = L s 1 1 s G 0 1 G 0 + ξ g 0 1 ξ g G 1 G L ,
where
L i s t h e L a p l a c e t r a n s f o r m
g 0 = g 0 u + v + s , ϑ + θ g = g u + v + s , ϑ + θ
G 0 = g 0 v + s , θ G = g v + s , θ
G 0 1 = g 0 v , θ = E e v X 0 e t 0 θ G 1 = g v , θ = E e v X 1 e t 1 t 0 θ
In the work of Dshalalow and Robinson [7], under the special assumptions made on S , such as t 0 = 0 and X 0 (the initial damage) a constant, S with position independent marking, X j   Exp η , δ j = t j t j 1   Exp τ , the marginal joint transform Φ 1 , 0 , v , 0 , θ   = E e v S μ θ t μ of the damage to the system on the failure and the failure time satisfies the formula
E e v S μ θ t μ = η τ η + v τ + θ e η L X 0 e v η τ τ + θ e X 0 η τ τ + θ
that further implies that the marginal pdf of S μ is
f S μ x = η e η x L 1 L , x
and the marginal pdf of t μ is
f t μ t = τ e η L X 0 e τ t I 0 2 η τ L X 0 t 1 R + t ,
where I 0 x is the modified Bessel function of order zero.
Furthermore, the PDF F t μ 1 t   = P t μ 1 t of the pre-failure time is
F t μ 1 t = e η L X 0 e τ t I 0 2 η τ L X 0 t + τ z = 0 t e τ z I 0 2 η τ L X 0 z d z

1.2.3. Extreme Shock Models

Gut and Hüsler [3,8] studied a shock model with the input of extreme shocks represented by a marked point process S = j = 1 X j ε t j with position independent marking. Under  S , set ξ x   = inf n : X n > x , with t ξ x being the hard failure time. Obviously, ξ x     Geo p x , where p x   =   P X 1 > x . Furthermore, the authors established E t ξ x = E δ 1 p x and the asymptotic behavior t ξ x ξ x E δ 1 , as  x sup u : F X u < 1 .

1.2.4. δ -Shock Models

Li and Kong [9] studied a system that fails when the time lag between two consecutive shocks decreases than some positive δ . The δ -shock policy is often implemented whenever shock damages are hard to observe. Here, external shocks arrive according to a Poisson point process of rate λ . Let M N 0 , N t be the associated Poisson counting process, and  m t   = min { δ j = t j t j 1 : M < j N t } . Furthermore, the first M shocks are not fatal and the M + j th shock is fatal if δ M + j < δ , j = 1 , 2 , . In other words, the system is considered with regard to the confined σ -algebra
F M = F   X 1 δ , , X M δ
Let ξ = inf n : X n < δ . Then, P ξ = h X 1 δ , , X M δ   = p q h M 1 , q = e λ δ . The authors give the PDF of t ν on F M and the asymptotic behavior of t ν / E t ν   Exp 1 as δ 0 .
Eryilmaz and Bayramoglu [10] studied a δ -shock system where external shocks arrive according to a renewal process with interrenewal times uniformly distributed. Earlier, Eryılmaz [11] considered a system with shocks arriving according to a marked Poisson process with position independent marking, under the generalized model proposed by Eryilmaz in which the failure occurs when the respective lengths of k consecutive interarrival times of the shocks is less than δ . That is, the system’s failure occurs at t ν , where
ν = inf { n : t n k + 1 t n k < δ , , t n t n 1 < δ }
The authors found explicit distributions of ν (geometric distribution of order k) and of failure time t ν .
Eryilmaz further embellished this model by introducing
μ = inf n : X n k + 1 > M , , X n > M
and setting ρ = ν μ , so that the system’s failure occurs at time t ρ . In the combination of two types of failures (mixed model), the formulas for ρ and failure time t ρ are not of closed form except for special asymptotic cases.

1.2.5. Run Shock Models

Input of shocks is specified by a marked point process S = k = 1 X n ε t n , with  δ k being interarrival times of the shocks. Such a system was introduced and studied by [12,13], with the failure of the system defined as follows. Given a critical region R q R , let q = 1 p = P X j R q and
ν q , k = min n : X n , X n 1 , , X n k + 1 R q
being a critical run in a string of k critical shocks that cause the failure of the system, with the failure time t ν q , k . The authors studied asymptotic behavior of ν and t ν , showing that
q k ν q , k e weakly as q 0
p q k ν q , k e weakly as k
q k t ν q , k exp E δ 1 weakly as q 0
p q k t ν q , k exp E δ 1 weakly as k
Furthermore, [14] studied asymptotics in a mixed run and cumulative shock model.
Ozkut and Eryilmaz [15] studied a so-called Marshall–Olkin run shock model that stemmed from Marshall and Olkin [16,17]. This system consists of two components that are subject to shocks from three different sources such that the shocks from the first and second source damage the first and second component, respectively, while the shocks from the third source affect either component. The produced shocks are classified as critical or non-critical. A component fails if it is subject to k consecutive critical shocks from the same source. The authors assume that interarrival times of the shocks follow a phase-type distribution and derive explicitly the distribution of the survival function of the two components.

1.2.6. Mixed Shock Models

Eryilmaz and Tekin [18] considered a system with input of shocks being a marked point process with and without position dependence. They introduce two thresholds, d 1 < d 2 , such that
μ = inf n : X n k + 1 > d 1 , , X n > d 1 ( run shock )
ν = inf n : X n > d 2 ( extreme shock )
ρ = μ ν
The failure time is t ρ . No assumption is made on the nature of that point process. They found a closed formula for the probability distribution for ρ having a phase type distribution, and for the special case of δ n = t n t n 1 being of the phase type and the point process with position independent marking, the failure time t ρ was proved to be of a phase type.
Zhao et al. [19] dealt with a system where shocks arrive according to a Poisson process of rate λ . Each shock can be one of the two types, valid (with probability p) or invalid (with probability q = 1 p ). A valid shock causes some random damage to the system, while invalid shock does not. An invalid shock is called δ -invalid if the time elapsed from the preceding shock is greater than δ . A so-called self-healing mechanism kicks in when i valid shocks are followed by k invalid shocks, one after another. In this case, the last of the i damages is healed. If the total number of valid shocks hits d , then a self-healing ability is discontinued (end of stage 1), after which all valid shocks are accumulated until its number reaches D > d and this is when the system fails (end of stage 2). The reliability of this system is analyzed by a finite Markov chain. μ gives the total number of shocks that end stage 1, ν is the total number of shocks until the system fails. The authors provide computational formulas for both. However, they used a simulation to derive the probability distribution of the failure t ν .

1.2.7. Shock and Degradation Models

In the above models, the interest was more on shocks, with very little to no attention paid to the degradation process. There are numerous articles on reliability systems entirely focused on soft failures merely due to degradation and no shocks. Additionally, there are some notable works with the failures due to both degradation and shocks and, in addition, mixed shocks.
Peng et al. [20] studied a system under the following assumptions.
(a)
Soft failures are caused by continuous wear degradation and vocational non-fatal and extreme shocks that arrive in accordance with a marked Poisson process S = j = 1 X j , W j ε t j of rate λ for its support counting measure j = 1 ε t j with position independent marking and with marks X j ’s and W j ’s independent of each other. A soft failure occurs when the cumulative degradation process (due to wear and periodic shocks X j ’s) crosses a fixed threshold H . Additionally, N t   = j = 1 ε t j 0 , t ;
(b)
Hard failures caused by extreme shocks W j (also referred to as shock loads) affect a different unit of the system and that causes a catastrophic failure. Only in this case does the system fail if any such shock exceeds some level D.
So, it is a mixed shock model of cumulative and extreme shocks—whichever comes first on top of a continuous degradation process Y t being affine with Y t   = φ + β t , where the initial value φ is a.s. a constant and the degradation rate β   N μ , σ 2 . The non-fatal shocks X j are Gaussian with parameters μ X , σ X 2 . As usual, no assumption is needed on the distribution of extreme shocks W j .
The total degradation accrued by the system, Y s t , is the sum of the degradation due to continual wear Y t and the instantaneous damages due to non-fatal and extreme shocks. The cumulative damage size due to non-fatal shocks alone until time t, S t , is given by S t   = X 1 + + X N t . This implies a cumulative shock model. The overall degradation of the system, considering both wear degradation and non-fatal shock damages (but excluding the extreme shocks), is Y s t   = Y t + S t .
Therefore, the system reliability at time t is
R t = P Y t   < H , N t   = 0   + j = 1 [ P W j < D ] j P Y t   + i = 1 j X i < H N t   = j
that presents a simple computational formula
R ( t ) = Φ H μ t φ σ t e λ t + j = 1 P W j < D j Φ H μ t + φ + i μ X σ 2 t 2 + j σ X 2 e λ t λ t j / j !
with a rightly truncating series.
Hao et al. [21] considered a similar system with the degradation process, with Y t   = μ t + σ B t being Brownian motion with drift. Non-fatal shocks arrive according to a marked Poisson process S = j = 1 X j ε t j of rate λ for its support counting measure j = 1 ε t j with position independent marking and with the absence of extreme shocks. The non-fatal shock damages X 1 , X 2 , are Gaussian with parameters μ X , σ X 2 . The reliability R t obeys a similar formula as in the work of Peng et al. [20].
In another mixed shock model with degradation, Rafiee et al. [22] proposed a system where the degradation process is affine Y t   = φ + β t . All (non-fatal and extreme) shocks arrive according to a marked Poisson process S = j = 1 X j , W j ε t j of rate λ for its support counting measure j = 1 ε t j with position independent marking and with the marks X j and W j being independent of each other. Non-fatal shocks X 1 , X 2 , not only accelerate degradation but now make the slopes β (degradation rate) vary as to β 1 , β 2 , . Extreme shocks form W 1 , W 2 , and they become fatal single shocks according to variable thresholds D 1 , D 2 , . Furthermore, β j + 1 = β j + η j ( η j are iid r.v.’s) while D j + 1 = D j θ j ( θ j are iid r.v.). The model also involves the δ -shock policy. The total failure occurs when one of the three conditions are met: (a) a W-shock is above a fixed D e , (b) a time lag between two consecutive shocks is less than δ , or (c) the degradation process Y ( t ) and non-fatal shocks combined cross H.
In the work of Hao et al. [23], the authors propose a mixed model under the degradation process Y ( t ) , a linear gamma process, that is Y t Y 0   [ G a α t , β ] , and non-fatal and fatal shocks arrive according to a marked Poisson process S = j = 1 X j , W j ε t j of rate λ 0 for its support counting measure j = 1 ε t j with position independent marking, with each X j   N μ X , σ X 2 . A soft failure occurs when Y t   + S   t   > H and a hard failure occurs when W N t > D t   = a ( Y t   +   S t ) + b . Note that a stochastic process Y t is Gamma with shape parameter α t and scale parameter β if Y t Y 0 has independent and stationary increments and its single-dimensional pdf is g x ; α t , β   = β α t x α t 1 e β x / Γ α t . Thus, in this case, the shape parameter α t   = α t is linear.
Yousefi et al. [24] studied a series system with n components, each subject to degradation and periodic external shocks. The magnitudes of shocks exerted on each particular component are iid Gaussian r.v. The shocks arrive simultaneously at all components according to a marked Poisson process S = j = 1 X j 1 , , X j n ε t j of rate λ for its support counting measure j = 1 ε t j with position independent marking. The system ages according to the gamma process Y t with shape parameter α t and scale parameter β .
Cao et al. [25] studied a system with the degradation process being gamma with shape parameter α t   = k 1 t k 2 and scale parameter β . Shocks arrive in accordance with a marked Poisson process of rate λ , with position independent marking, and the marks with truncated Gaussian distributions. A soft failure occurs whenever the overall degradation level exceeds a fixed threshold H.
Oliveira et al. [26] studied a run shock two-component series system with two Marshall–Olkin shock structures [16,17]. According to Marshall and Olkin [17], there are three independent shock sources. A fatal shock from source 1 destroys component 1, which occurs at time W 1 ; a shock from source 2 destroys component 2, which occurs at time W 2 ; a shock from source 3 destroys both components, which occurs at time W 3 . In this case, the random lifetime of component 1 is X 1 = W 1 W 3 , while the random lifetime of component 2 is X 2 = W 2 W 3 .
In the other variant of Marshall and Olkin of [16], each component is subjected to independent stresses, say U 1 and U 2 , and the system has overall stress U 3 independently transmitted equally to both components. The observed shocks to either of the two components are X 1 = U 1 U 3 and X 2 = U 2 U 3 . In both structures, the authors assumed uniparametric Lindley distributions for W j and U j ( j = 1 , 2) with parameters β j ; for the latent random variables ( W 3 or U 3 ), they assumed two probability distributions: an exponential distribution and a Lindley distribution, both with parameter β 3 . Moreover, the dependence structure for the random variables X 1 and X 2 is related to the common source of shock or stress 3. Lifetime for component i is T i . The reliability function of the system is R t   = P T 1 T 2 > t . Note that a r.v. X is said to have a Lindley distribution with parameter β if its pdf is
f X x = β 2 1 + β 1 + x e β x 1 0 , x
for β > 0 .

1.3. Our Model

In the present paper, our degradation process is linear or affine with a constant slope (referred to as degradation rate), albeit under no assumptions upon the nature of shocks other than being represented by a general marked random measure with position dependent marking. (Of course, any marked Poisson process with position independent marking fits into this category.) Additionally, we are after a closed form expression for the failure time distribution and the overall damage upon the failure. From this context, with our focus on an analytical solution, the reliability model is very different from the mainstream literature. Our claim for analytical (rather than computational) tractability is due to a unique approach that we believe offers further embellishments and modifications applied to more general reliability systems such as those with multiple components (in series or parallel).
So, our system is driven by a deterministic degradation real-valued process Y ( t ) with a constant slope a that describes the well-being of the system. Under no further assumption on the system, the process Y ( t ) crosses a given threshold M at time 1 a [ M Y ( 0 ) ] . However, under random shocks occurring at random times, the crossing of M will happen at a random and possibly earlier time.
To be more specific, consider a piece of equipment that ages as it wears, meaning that it eventually deteriorates and becomes useless when its operating power Y ( t ) = Y ( 0 ) + a t crosses some positive sustainability threshold M. However, in a real-world system, the equipment deteriorates much faster due to occasional shocks occurring at random times t 1 , t 2 , of random magnitudes X 1 , X 2 , that altogether form a marked point process S = k = 1 X k ε t k ( ε b is the point mass) with position dependent marking.
Without aging, the crossing would be a more familiar task previously analyzed by the authors of [6,7], but along with aging it becomes a more challenging problem. The question is when the (soft) failure takes place and to what detriment. The latter may be of independent interest if the condition of an underlying device is beyond any repair or if it can still be fixed.
Denote
A n = k = 1 n X k + a δ k , where δ k = t k t k 1
and
ν = inf { n N : A n M }
The time dependent degradation of the system can be formalized by the marked random measure
A = k = 1 X k + a δ k ε t k , with A ( t ) = A [ 0 , t ] ,
where A ( t ) is the associated continuous time parameter process. No assumption is made on the nature of the point process { t k } and marks X k ’s and δ k ’s except that they are all mutually dependent.
From the Figure 1, the failure can take place at the crossing of M or by an excess.
Namely, if  A ν 1 + a δ ν M , then the crossing takes place in interval ( t ν 1 , t ν ) with crossing at time τ ν = t ν 1 + ( M A ν 1 ) / a , so that A ( τ ν ) = M sharp. If  A ν 1 + a δ ν < M , then the crossing of the horizontal line of height M takes place exactly at time t ν (the ν th shock) exceeding M and valued at A ν . Thus, A ν = A ( τ ν ) is the excess level of cumulative damage above threshold M. The location of A ν can be arbitrarily higher than M (in reality, process A ( t ) runs downhill and it can hit zero at worst).

1.3.1. Our Methodology

Our techniques fall into the category of fluctuation theory in the context of random walk processes. However, our model is far from a classical random walk where a generic particle or walker moves along a d-dimensional rectangular deterministic integer grid (also known as a lattice). The walker moves along all possible directions with equal probability and only one step at the time. It is placed in a bounded closed subset B of Z d and the associated problem is to find its escape location from that set B. In our present model, the walker moves along a randomly generated grid that is not even rectangular. (See more on different variants of random walks in Dshalalow and White [27].) Here, the walker moves in R 2 within a random rectangular set B = ( 0 , τ ν ] × ( 0 , M ] R + 2 and it escapes at the random point ( τ n u , M ) or ( t ν , A n u ) . (See the walker’s path in Figure 1).
The model as set is not very common in the context of random walks, although as far as grids, some walks allow looser configurations. Yet, this is a new identification of reliability models as random walks. Secondly, we apply and embellish the theory of fluctuations to arrive at analytically closed formulas and establish the main functional for the joint probability distribution of the first passage time (called so in fluctuation theory) and the position of walker or escape location associated with the failure time and the extent of the overall damage, respectively.
Our approach allows further generalizations by including other types of shocks. We earlier worked on multidimensional walks (cf. Dshalalow and White [27]), only this time the walker’s movements are not rectangular. The present paper points to ways of their modifications.

1.3.2. Paper’s Layout

The paper is organized as follows. In Section 2, we formalize the system with a linear degradation process combined with cumulative shock model of non-fatal shocks arriving according to a general marked point process S with position dependent marking that we introduced in the beginning of this section. We use an approach previously developed by the authors in [5,6,7] for a pure cumulative shock model but then expressed in terms of a random walk using fluctuation analysis. In this case, the model embellished by the presence of a degradation process required a quite different analysis, but still within the framework of the theory of fluctuations. We obtained a closed form functional of the joint distribution of the time failure and the detriment of the damage to the system, among other useful characteristics, like the status of the system on the shock preceding the fatal one or before the real crossing of the critical threshold. In Section 3, we support our claim of analytical tractability of the results by discussing a special case of the S of a marked Poisson process with exponentially distributed marks. The results are confirmed through a comparison with simulation in Section 4.

2. Fluctuation Analysis of the Linear Degradation Process with Shocks

Consider a piece of equipment that ages and degrades at a deterministic rate a of damage per unit time, eventually degrading to a point that it becomes unuseful when its accumulated damage crosses a positive threshold M. In this idealized scenario, this occurs precisely at time M / a . Practically, however, the piece of equipment deteriorates much faster due to occasional irreversible shocks of random magnitudes X 1 , X 2 , occurring at random times t 1 , t 2 , , which together form a marked point process. Estimating the crossing time of M and excess level of a marked point process crossing M by a piecewise constant jump process is a familiar task in the literature, but combining this with deterministic degradation poses a more challenging problem. The question is when the failure takes place and to what detriment. The latter may be of interest if the condition of the device is degraded beyond any repair.
In the context of the shock random measure S = k = 1 X k ε t k (introduced in subsection Our Model of Section 1), we assume that the support counting measure N = k = 1 ε t k represents a renewal process. That is, the times between consecutive shocks
δ k = t k t k 1
are i.i.d. for k 1 and assuming t 0 = 0 . Then, the cumulative damage due to continual degradation and periodic shocks is
A = k = 1 X k + a δ k ε t k ,
with A ( t ) = A [ 0 , t ] being the associated continuous time parameter counting process describing the status of damage at any time t. Note that A is not assumed to be with position independent marking, implying that the marks X k and δ k are not independent for any fixed k but the vectors ( X k , δ k ) , k = 1 , 2 , are i.i.d. with a common joint Laplace–Stieltjes transform (LST)
γ ( α , θ ) = E e α ( X k + a δ k ) θ δ k ,
which is assumed to be known or at least not difficult to estimate. With 
A n = k = 1 n X k + A δ k
define
ν = inf { n : A n M }
the minimum number of shocks that along with degradation leads to the failure.
The failure time is τ ν and it can be one of the two values. It either equals t ν if there no failure in interval ( t ν 1 , t ν ) or it is t ν 1 + ( M A ν 1 ) / a if the failure takes place in ( t ν 1 , t ν ) . Thus,
τ ν = t ν , if A ν 1 + a δ ν < M t ν 1 + ( M A ν 1 ) / a , if A ν 1 + a δ ν M
We denote by S ν the failure damage of S upon the failure time τ ν . It obviously is
S ν = A ν , if A ν 1 + a δ ν < M M , if A ν 1 + a δ ν M
We target the joint LST Φ ν of the joint distribution of the failure time τ ν , the total damage to the system S ν upon its failure (that is the crossing level), the pre-failure time (the time of the ( ν 1 ) st shock preceding the failure), and the total damage to the system A ν 1 brought by the ν th shock. Note that A ν 1 is not the biggest damage prior to the failure unless τ ν = t ν . The last two parameters can be of independent interest. Thus,
Φ ν ( α , β , ϑ , θ ) = E e α A ν 1 β S ν ϑ t ν 1 θ t ν
Theorem 3.
The joint LST Φ ν ( α , β , ϑ , θ ) of the degradation upon the shock before the failure, the damage upon failure, the time of the shock before the failure, and the failure time satisfies
Φ ν ( α , β , ϑ , θ ) = L x 1 1 x + β + θ a F d ( x ) ( M ) + L x 1 1 x F s ( x ) ( M ) ,
where
F d ( x ) = 1 γ ( 0 , a β + a x + θ ) 1 γ ( α + β + x , ϑ + θ )
F s ( x ) = γ ( β , θ + a x ) γ ( β + x , θ ) 1 γ ( α + β + x , ϑ + θ ) .
Proof. 
First, we introduce a set of failure indices { ν ( p ) = inf { n : A n > p } : p > 0 } so that ν = ν ( M ) , which creates a set of functionals { Φ ν ( p ) : p > 0 } . We will derive an expression for Φ ν ( p ) and then use operational calculus to find a formula for Φ ν .
The functional Φ ν ( p ) can be computed as a convenient partition of the sample space. The partition includes two subsets. First is the set of events where A j 1 + a δ j p A j 1 and ν ( p ) = j is given, i.e., the failure occurs between the ( j 1 ) th and jth shocks as a result of the constant degradation, which we call a degradation failure. Formally,
{ { ν ( p ) = j } { A j 1 + a δ j p A j 1 } : j = 1 , 2 , }
Second is the set of events where A j 1 + a δ j < p A j 1 and ν ( p ) = j is given, i.e., the failure occurs upon the jth shock, which we call a shock failure. Formally,
{ { ν ( p ) = j } { A j 1 + a δ j < p A j } : j = 1 , 2 , }
Therefore, we have
Φ ν ( p ) = j = 1 E e α A ν ( p ) 1 β S ν ( p ) ϑ t ν ( p ) 1 θ τ ν ( p ) 1 { ν ( p ) = j } 1 { A j 1 + a δ j p A j 1 } + j = 1 E e α A ν ( p ) 1 β S ν ( p ) ϑ t ν ( p ) 1 θ τ ν ( p ) 1 { ν ( p ) = j } 1 { A j 1 + a δ j < p A j } = j = 1 Φ j d p + Φ j s p .
The first expression, corresponding to the degradation failure, simplifies to
Φ j d p = E e α A j 1 β S j ϑ t j 1 θ τ j 1 { A j 1 + a δ j p A j 1 } = E e α A j 1 β p ϑ t j 1 θ t j 1 + 1 a ( p A j 1 ) 1 { A j 1 + a δ j p A j 1 } = e β + θ a p E e α A j 1 ϑ t j 1 θ t j 1 1 a A j 1 1 { A j 1 + a δ j p A j 1 } = e β + θ a p E e α θ a A j 1 ( ϑ + θ ) t j 1 1 { A j 1 + a δ j p A j 1 } .
Next, we apply a modified Laplace transform to Φ j d p acting on variable p, which bypasses all terms except those that depend on p,
x + β + θ a L p e β + θ a p 1 { A j 1 + a δ j p A j 1 } ( x ) = x + β + θ a p = A j 1 A j 1 + a δ j e β + θ a p e x p d p = e x + β + θ a A j 1 e x + β + θ a A j 1 + a δ j = e x + β + θ a A j 1 1 e a x + a β + θ δ j .
Therefore, we have
F j d ( x ) = L p Φ j d p ( x ) = E e α θ a A j 1 ( ϑ + θ ) t j 1 e x + β + θ a A j 1 1 e a x + a β + θ δ j = E e α + β + x A j 1 ( ϑ + θ ) t j 1 1 e a β + a x + θ δ j = γ j 1 ( α + β + x , ϑ + θ ) 1 γ ( 0 , a β + a x + θ ) .
Then,
F d ( x ) = j = 1 γ j 1 ( x + α + β , ϑ + θ ) 1 γ ( 0 , a x + a β + θ ) = 1 γ ( 0 , a x + a β + θ ) 1 γ ( x + α + β , ϑ + θ ) ,
which completes the proof for Formula (43).
The second expression, corresponding to the shock failure, simplifies as
Φ j s p = E e α A j 1 β S j ϑ t j 1 θ τ j 1 { A j 1 + a δ j < p A j } = E e α A j 1 β A j ϑ t j 1 θ t j 1 { A j 1 + a δ j < p A j } .
Next, we apply a different modified Laplace transform to Φ j s p acting on variable p, which bypasses all terms except those that depend on p,
x L p 1 { A j 1 + a δ j < p A j } ( x ) = x p = A j 1 + a δ j A j e x p d p = e x ( A j 1 + a δ j ) e x A j = e x ( A j 1 + a δ j ) 1 e x X j .
Therefore, we have
F j s ( x ) = L p Φ j s p ( x ) = E e α A j 1 β A j ϑ t j 1 θ t j e x ( A j 1 + a δ j ) 1 e x X j = E e ( α + β + x ) A j 1 β ( X j + a δ j ) ( ϑ + θ ) t j 1 θ δ j e a x δ j 1 e x X j = E e ( α + β + x ) A j 1 ( ϑ + θ ) t j 1 e β ( X j + a δ j ) ( a x + θ ) δ j 1 e x X j = E e ( α + β + x ) A j 1 ( ϑ + θ ) t j 1 e β ( X j + a δ j ) ( a x + θ ) δ j E e ( α + β + x ) A j 1 ( ϑ + θ ) t j 1 e ( β + x ) ( X j + a δ j ) ( a x + θ ) δ j = γ j 1 ( α + β + x , ϑ + θ ) γ ( β , θ + a x ) γ ( β + x , θ ) .
Then,
F s ( x ) = j = 1 γ j 1 ( α + β + x , ϑ + θ ) γ ( β , θ + a x ) γ ( β + x , θ ) = γ ( β , θ + a x ) γ ( β + x , θ ) 1 γ ( α + β + x , ϑ + θ ) ,
as a geometric series, which establishes Formula (44).
Therefore, applying the inverse operators evaluated at the true threshold M yields
Φ ν ( α , β , ϑ , θ ) = L x 1 1 x + β + θ a F d ( x ) ( M ) + L x 1 1 x F s ( x ) ( M )
as required.    □
In the following section, we will make some assumptions about the deterioration process to demonstrate the previous theorem can be applied to realistic special cases to find explicit probabilistic results about the deterioration.

3. Results for a Dual-Exponential Shock Process

To justify our claim for analytical tractability, we consider a special case under the following assumptions. Suppose the times between each shock δ k are exponentially distributed with parameter λ and damage due to each shock X k being independent of each other. Using the terminology of random measures, the marked point process S = k = 1 X k + a δ k ε t k is with position independent marking. Then,
γ ( α , θ ) = E e α ( X k + a δ k ) θ δ k = E e α X k ( a α + θ ) δ k = E e α X k E e ( a α + θ ) δ k .
Suppose further that the times δ k are exponentially distributed with parameter λ and the damage due to the shocks X k are exponentially distributed with parameter μ . Then,
γ ( α , θ ) = λ λ + α μ μ + a α + θ .
In this special case, we establish an explicit formula for Φ ν ( α , β , ϑ , θ ) .
Proposition 1.
For the marked point process S of the evolution of deterioration, let S have position independent marking. Furthermore, if the times between shocks δ k are exponentially distributed with parameter λ and the impacts of the shocks X k ’s are exponentially distributed with parameter μ, then the joint LST Φ ν ( α , β , ϑ , θ ) of the deterioration upon the shock before the failure, the deterioration upon failure, the time of the shock before the failure, and the failure time satisfies
Φ ν ( α , β , ϑ , θ ) = Φ d ( α , β , ϑ , θ ) + Φ s ( α , β , ϑ , θ ) ,
where
Φ d ( α , β , ϑ , θ ) = f ( α , ϑ , θ ) f ( α , ϑ , θ ) a λ μ e β + μ + θ a M + a λ μ 2 f ( α , ϑ , θ ) a λ μ e ( α + β ) M × 2 a α + a λ μ + ϑ θ D ( ϑ , θ ) e r ( ϑ , θ ) M e s ( ϑ , θ ) M e r ( ϑ , θ ) M + e s ( ϑ , θ ) M
and
Φ s ( α , β , ϑ , θ ) = α λ μ λ + β a λ μ θ g ( α , ϑ , θ ) α g ( α , ϑ , θ ) λ μ e ( λ + β ) M + λ μ λ + β a λ μ θ f ( α , ϑ , θ ) f ( α , ϑ , θ ) a λ μ e β + μ + θ a M λ 2 μ 2 2 ( λ + β ) 1 f ( α , ϑ , θ ) a λ μ 1 α g ( α , ϑ , θ ) λ μ e ( α + β ) M × [ 2 a 2 α 2 + 2 a λ μ + 2 a α ϑ a λ ϑ + μ ϑ + ϑ 2 + ϑ θ D ( ϑ , θ ) e r ( ϑ , θ ) M e s ( ϑ , θ ) M 2 a α + ϑ e r ( ϑ , θ ) M + e s ( ϑ , θ ) M ]
under the notation
f ( α , ϑ , θ ) = a 2 α 2 + a 2 α λ a α μ + a α ϑ + a λ ϑ μ ϑ a α θ ϑ θ
g ( α , ϑ , θ ) = a α a λ + μ + ϑ + θ ,
where r ( ϑ , θ ) and s ( ϑ , θ ) are the roots of the polynomial a x 2 + ( a λ + μ + ϑ + θ ) x + λ ( ϑ + θ ) and D ( ϑ , θ ) is the discriminant of the polynomial, with 
r ( ϑ , θ ) = ( a λ + μ + ϑ + θ ) D ( ϑ , θ ) 2 a
and
s ( ϑ , θ ) = ( a λ + μ + ϑ + θ ) + D ( ϑ , θ ) 2 a .
Proof. 
For the first term from Theorem 3, notice
1 x + β + θ a F d ( x ) = 1 x + β + θ a 1 γ ( 0 , a β + a x + θ ) 1 γ ( α + β + x , ϑ + θ ) = 1 x + β + θ a 1 μ μ + a β + a x + θ 1 λ λ + α + β + x μ μ + a α + a β + a x + ϑ + θ .
As we see, this is simply a rational function with degree 2 (in x) in the numerator and degree 3 in the denominator, which means the inverse Laplace transform can be computed easily with a partial fraction decomposition to find
Φ d ( α , β , ϑ , θ ) = L x 1 1 x + β + θ a F d ( x ) ( M ) = f ( α , ϑ , θ ) f ( α , ϑ , θ ) a λ μ e β + μ + θ a M + a λ μ 2 f ( α , ϑ , θ ) a λ μ e ( α + β ) M × 2 a α + a λ μ + ϑ θ D ( ϑ , θ ) e r ( ϑ , θ ) M e s ( ϑ , θ ) M e r ( ϑ , θ ) M + e s ( ϑ , θ ) M .
Similarly,
1 x F s ( x ) = 1 x γ ( β , θ + a x ) γ ( β + x , θ ) 1 γ ( α + β + x , ϑ + θ ) = 1 x λ λ + β μ μ + a β + a x + θ λ λ + β + x μ μ + a β + a x + θ 1 λ λ + α + β + x μ μ + a α + a β + a x + ϑ + θ .
This is, again, a rational expression—this time with degree 4 in the denominator, so its inverse Laplace transform can be computed easily with a partial fraction decomposition to find
Φ s ( α , β , ϑ , θ ) = L x 1 1 x F s ( x ) ( M ) = α λ μ λ + β a λ μ θ g ( α , ϑ , θ ) α g ( α , ϑ , θ ) λ μ e ( λ + β ) M + λ μ λ + β a λ μ θ f ( α , ϑ , θ ) f ( α , ϑ , θ ) a λ μ e β + μ + θ a M λ 2 μ 2 2 ( λ + β ) 1 f ( α , ϑ , θ ) a λ μ 1 α g ( α , ϑ , θ ) λ μ e ( α + β ) M × [ 2 a 2 α 2 + 2 a λ μ + 2 a α ϑ a λ ϑ + μ ϑ + ϑ 2 + ϑ θ D ( ϑ , θ ) e r ( ϑ , θ ) M e s ( ϑ , θ ) M 2 a α + ϑ e r ( ϑ , θ ) M + e s ( ϑ , θ ) M ] .
The sum of the two inverse Laplace transforms evaluated at M, i.e.,  Φ ν ( α , β , ϑ , θ ) , is the sum of (67) and (69) by Theorem 3.   □
The result above is a very information-rich functional in this special case and can give us expressions for many interesting probabilistic results. We outline a few such results below. First, the probability of failure occurring between shocks and probability of failure upon a shock.
Corollary 1.
Under the assumptions of Proposition 1, the probabilities of degradation failure and shock failure are
P ( d e g r a d a t i o n f a i l u r e ) = a λ a λ + μ + μ a λ + μ e λ + μ a M
and
P ( s h o c k f a i l u r e ) = μ a λ + μ μ a λ + μ e λ + μ a M .
Proof. 
For the probability of a degradation failure, simply let α = β = ϑ = θ = 0 in Φ d ( α , β , ϑ , θ ) to find
Φ d ( 0 , 0 , 0 , 0 ) = P ( degradation failure ) = f ( 0 , 0 , 0 ) f ( 0 , 0 , 0 ) a λ μ e μ a M + a λ μ 2 f ( 0 , 0 , 0 ) a λ μ × a λ μ D ( 0 , 0 ) e r ( 0 , 0 ) M e s ( 0 , 0 ) M e r ( 0 , 0 ) M + e s ( 0 , 0 ) M = a λ a λ + μ + μ a λ + μ e a λ + μ a M .
For the probability of shock failure, simply let α = β = ϑ = θ = 0 in Φ s ( α , β , ϑ , θ ) to find
Φ s ( 0 , 0 , 0 , 0 ) = μ a λ μ f ( 0 , 0 , 0 ) f ( 0 , 0 , 0 ) a λ μ e μ a M + λ μ 2 1 1 f ( 0 , 0 , 0 ) a λ μ a D ( 0 , 0 ) e r ( 0 , 0 ) M e s ( 0 , 0 ) M
= μ a λ + μ μ a λ + μ e λ + μ a M .
   □
Interestingly, as the threshold M grows, the probabilities approach constants a λ a λ + μ and μ a λ + μ . Of course, the sum of the probabilities is 1 since exactly one of these disjoint events will almost surely occur.
Corollary 2.
Under the assumptions of Proposition 1, the marginal LST of the failure time τ ν is
Φ ν ( 0 , 0 , 0 , θ ) = E e θ τ ν = 1 2 e r ( 0 , θ ) M + e s ( 0 , θ ) M + a λ + μ θ D ( 0 , θ ) e s ( 0 , θ ) M e r ( 0 , θ ) M .
Proof. 
Let α = β = ϑ = 0 , then
Φ ν ( 0 , 0 , 0 , θ ) = Φ d ( 0 , 0 , 0 , θ ) + Φ s ( 0 , 0 , 0 , θ ) ,
where
Φ d ( 0 , 0 , 0 , θ ) = f ( 0 , 0 , θ ) f ( 0 , 0 , θ ) a λ μ e μ + θ a M + a λ μ 2 f ( 0 , 0 , θ ) a λ μ × a λ μ θ D ( 0 , θ ) e r ( 0 , θ ) M e s ( 0 , θ ) M e r ( 0 , θ ) M + e s ( 0 , θ ) M = 1 2 a λ μ θ D ( 0 , θ ) e r ( 0 , θ ) M e s ( 0 , θ ) M e r ( 0 , θ ) M + e s ( 0 , θ ) M
and
Φ s ( 0 , 0 , 0 , θ ) = μ a λ μ θ f ( 0 , 0 , 0 ) f ( 0 , 0 , 0 ) a λ μ e μ + θ a M + μ f ( 0 , 0 , 0 ) a λ μ a λ μ D ( 0 , θ ) e r ( 0 , θ ) M e s ( 0 , θ ) M = μ D ( 0 , θ ) e s ( 0 , θ ) M e r ( 0 , θ ) M .
Summing expressions (76) and (77) above, we get E e θ τ ν .   □
This LST can easily yield means and moments, as we see below.
Proposition 2.
Under the assumptions of Proposition 1, the mean and variance of the failure time τ ν are
E τ ν = λ M a λ + μ + μ ( a λ + μ ) 2 1 e λ + μ a M = λ M + P ( s h o c k f a i l u r e ) a λ + μ
and
V a r ( τ ν ) = 2 λ μ M ( a λ + μ ) 3 + μ μ 4 a λ ( a λ + μ ) 4 + 2 μ 2 a 2 λ + a 2 λ 2 M 2 μ 2 M a ( a λ + μ ) 4 e λ + μ a M 2 μ 2 ( a λ + μ ) 4 e 2 λ + μ a M .
Proof. 
The first two moments of τ ν can readily be found by taking the first two derivatives of E e θ τ ν , managing the signs appropriately and taking limits, yielding the results above.   □
Let us now consider the failure damage S ν of the process.
Corollary 3.
Under the assumptions of Proposition 1, the marginal LST of the failure damage S ν is
E e β S ν = Φ ν ( 0 , β , 0 , 0 ) = e β M a λ a λ + μ + μ a λ + μ e λ + μ a M + λ λ + β μ a λ + μ 1 e λ + μ a M .
Proof. 
Let α = ϑ = θ = 0 , then
Φ ν ( 0 , β , 0 , 0 ) = Φ d ( 0 , β , 0 , 0 ) + Φ s ( 0 , β , 0 , 0 )
The first term simplifies to
Φ d ( 0 , β , 0 , 0 ) = f ( 0 , 0 , 0 ) f ( 0 , 0 , 0 ) a λ μ e β + μ a M + a λ μ 2 f ( 0 , 0 , 0 ) a λ μ e β M × a λ μ D ( 0 , 0 ) e r ( 0 , 0 ) M e s ( 0 , 0 ) M e r ( 0 , 0 ) M + e s ( 0 , 0 ) M = e β M a λ a λ + μ + μ a λ + μ e λ + μ a M
and the second term simplifies to
Φ s ( 0 , β , 0 , 0 ) = 1 ( λ + β ) e β M λ μ D ( 0 , 0 ) e r ( 0 , 0 ) M e s ( 0 , 0 ) M = λ λ + β e β M μ a λ + μ 1 e λ + μ a M .
Summing expressions (82) and (83), we obtain Formula (80) for the marginal LST of the failure damage S ν .   □
Next, we use the LST of S ν to find simple explicit formulas of its mean and variance.
Corollary 4.
Under the assumptions of Proposition 1, the mean and variance of the failure damage S ν are
E S ν = M + 1 λ μ a λ + μ 1 e λ + μ a M = M + E δ 1 P ( s h o c k f a i l u r e )
and
V a r ( S ν ) = 2 a μ λ ( a λ + μ ) 2 1 e λ + μ a M + μ 2 λ 2 ( a λ + μ ) 2 1 e 2 λ + μ a M .
Proof. 
The first two moments of S ν can readily be found by taking the first two derivatives of E e β S ν , managing the sign appropriately and taking limits, yielding the results above.    □
These results reveal a number of intuitive facts. The mean of S ν makes clear the obvious fact that the failure damage will always be at least M, but more interestingly shows that the excess of the failure damage is precisely E δ 1 P (shock failure), the mean time between shocks multiplied by the probability of a shock failure, meaning it cannot surpass the mean time between shocks. In addition, as M grows, the mean approaches M as shocks become small relative to M and the variance approaches the constant 2 a μ λ + μ 2 λ 2 ( a λ + μ ) 2 .

4. Comparison with Stochastic Simulation

In this section, predictions from the formulas for probabilities, means, and variances derived for the exponential special case in the previous section will be shown to agree with Monte Carlo simulation of the process under numerical assumptions on the parameters: the degradation rate a, the parameter of the exponentially distributed time between shocks λ , the parameter of the exponentially distributed shock damage μ , and the threshold M.
Simulation code is written in the Python language with the NumPy library for numerical computation. The code for simulating paths of the process has been included in Appendix A along with diagrams and a description of the code. The experiments run in this section have been made available publicly by the authors online: see the Data Availability Statement of this article.
Suppose λ = 2 and μ = 1 . Figure 2 below focuses on the probability of a degradation failure (rather than a shock failure). The curves correspond to the predictions for probability of degradation failure from Corollary 1, while the dots represent empirical probabilities corresponding to simulations of 10,000 paths of the process for each pair of parameters ( a , M ) { 0.5 , 1 , 5 } × { 0.5 , 0.75 , 1 , , 3 } .
Notice the degradation failure probabilities have an inverse relationship with the threshold M as a quick degradation failure become less probable and shocks become increasingly probable causes of failures. Further, we plotted the limiting probabilities as dotted horizontal lines of height a λ a λ + μ and note the probabilities all approach these, as predicted. Further, probabilities are inversely related the degradation rate a. Clearly, this confirms the intuitive fact that a higher degradation rate should result in more degradation failures.
Suppose λ = 2 and μ = 3 . Figure 3 below focuses on the mean and variance of the failure time τ ν . The curves correspond to the predictions for the mean and standard deviation from Proposition 2, while the dots represent empirical means and standard deviations corresponding to simulations of 10,000 paths of the process for each pair of parameters ( a , M ) { 0 , 0.01 , , 0.5 } × { 1 , 5 , 10 , 25 } .
Note that the mean failure time shrinks as the degradation rate a grows and failures occur more quickly while the threshold M and mean failure time have a positive relationship as failures occur more slowly. We see similar trends with the standard deviation: variability decreases as the degradation rate grows and the deterministic part of the degradation is more influential while variability grows as the threshold grows and more shocks will typically occur in paths.
Suppose that λ = 2 and μ = 3 . Figure 4 below focuses on the failure damage S ν . The curves correspond to the predictions for the mean and standard deviation from Proposition 4, while the dots represent empirical means and standard deviations corresponding to simulations of 1,000,000 paths of the process for each pair of parameters ( a , M ) { 0 , 0.01 , , 5 } × { 1 , 5 } . (Note we generated more paths in this case because the empirical standard deviations took longer to converge to the true predicted values.)
Here, as degradation rate grows, we saw that shock failures become less probable previously, which causes the mean failure damage to decrease toward M, as degradation failures with positions equal to M sharp are more likely. The standard deviation of failure damage has an inverse relationship with a and a positive relationship with M for reasons analogous to the same trends in the mean failure time.

5. Conclusions

In this paper, we studied a reliability system with linear degradation and external shocks (cumulative shocks model with degradation) that further accelerate its deterioration and lead to an inevitable failure (soft failure) when the overall degradation crosses or exceeds an M, its sustainability threshold. The gradual degradation and random damages to the system can be modeled by a random walk type process, although not in the classical sense where a walker moves along a rectangular grid. In a more recent work, random walk is considered on a variety of configurations, more general than rectangular, cf. Dshalalow and White [27]. Now the random walk modeling is not just an interpretation, but it suggested fluctuation analysis as a primary tool for dealing with this class of models, with the benefit of obtaining a closed form distribution of targeted parameters such as the failure time of the system and the total damage to the system at the failure time. The latter is noteworthy information in the event the targeted component can be repaired or completely discarded.
The time dependent status of degradation of the system can be formalized by the marked point random measure A = k = 1 X k + a δ k ε t k with position dependent marking. Here, t 1 , t 2 , is the arrival process of non-fatal shocks of respective magnitudes X 1 , X 2 , . The degradation is further enhanced by the random increment a δ k in interval ( t k 1 , t k ) where a is the constant degradation rate and δ k is the length of the time interval ( t k 1 , t k ) . We allow X k to depend on δ k as well, although in some special cases we can lift this assumption. No assumption is made about the nature of t k , X k and thus δ k . Note that position dependent marking refers to the sequence { ( X k , δ k = t k t k 1 ) } which is comprised of independent vectors but with mutually dependent components. The system fails at some random time τ ν when the continuous time parameter degradation process A ( t ) = A [ 0 , t ] crosses or exceeds some M. Note that the crossing of M can occur upon one of the arriving shocks at t 1 , t 2 , , say at t ν or within the interval ( t ν 1 , t ν ) that we associate with a time τ ν . The overall damage to the system upon its failure is measured by the quantity S ν being M sharp or A ( τ ν ) dependent on whether τ ν ( t ν 1 , t ν ) or τ ν = t ν .
We targeted the joint transform Φ ν ( α , β , ϑ , θ ) = E e α A ν 1 β S ν ϑ t ν 1 θ t ν that carries two more quantities, A ν 1 and t ν 1 , giving us information about the system upon a shock preceding the time of failure. We arrived (Section 2) at the formula
Φ ν ( α , β , ϑ , θ ) = L x 1 1 x + β + θ a 1 γ ( 0 , a β + a x + θ ) 1 γ ( α + β + x , ϑ + θ ) ( M ) + L x 1 1 x γ ( β , θ + a x ) γ ( β + x , θ ) 1 γ ( α + β + x , ϑ + θ ) ( M ) ,
where
γ ( α , θ ) = E e α ( X k + a δ k ) θ δ k ,
is the joint transform of X k and δ k , which is supposed to be known once we specify the { t k } and { X k } , and  L x 1 is the inverse Laplace transform in variable x.
To support our claim that the above formula is closed form, we considered in Section 3 a special case of position independent marking and the assumptions on δ k and X k as being exponential. Under these assumptions, the respective transforms were inverted in a rather compact formula, further extended by marginal distributions of the failure time τ ν and the failure damage S ν , and their means and variances. In Section 4, we presented numerical results for several concrete numbers and compared them with simulated results which all agreed.
The overall results are explicit, general, and use fluctuation techniques different from the mainstream literature. The approach allows us to further embellish the current system making it a series system of two or more components applying and embellishing previously developed techniques of stochastic games  (Dshalalow [28]).

Author Contributions

Conceptualization, J.H.D.; methodology, J.H.D. and R.T.W.; validation, R.T.W.; formal analysis, J.H.D. and R.T.W.; investigation, J.H.D. and R.T.W.; writing—original draft preparation, J.H.D. and R.T.W.; writing—review and editing, J.H.D. and R.T.W.; visualization, J.H.D. and R.T.W. Both authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The authors have made all code and simulated datasets used in the paper can be found in the the GitHub repository https://github.com/rtwhite1546/ShockDegradationModel, accessed on 18 August 2021.

Acknowledgments

The authors are indebted to anonymous referees whose insightful comments and suggestions contributed a great deal of improvement to the clarity and presentation of our work.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

   The following abbreviations are used in this manuscript:
General Acronyms
i.i.d.independent and identically distributed
LSTLaplace-Stieltjes transform
PDFprobability distribution function
pdfprobability density function
PGFprobability-generating function
r.v.random variable
Notation in the Literature Review
S = k = 1 X n ε t n soft shocks marked point process in cumulative shock model
N = k = 1 ε t n associated support counting measure of shocks’ arrivals
S ( t ) = S [ 0 , t ] the number of shocks in time interval [ 0 , t ]
ε a Dirac point mass (unity measure)
X n magnitudes of soft shocks
t n the time of the nth shock
δ n = t n t n 1
S n = k = 1 n X n
Y ( t ) pure degradation process
W n magnitudes of hard shocks
W = j = 1 W j ε t j marked point measure of hard shocks
δ lower threshold in δ -shock policy models
xfailure threshold
Msoft failure threshold
Dhard failure threshold
ν ( x ) number of soft shocks until soft failure with respect to threshold x
ξ ( x ) number of hard shocks until hard failure with respect to threshold x
A = [ 0 , M ) the set that process S escapes
Notation in Our Model
ν number of the first shock where degradation exceeds M
t ν 1 pre-failure time
t ν soft failure time (the escape from set A)
S ν 1 pre-failure cumulative damage
S ν cumulative damage to the system at the failure
Φ ν = E ξ ν u S ν 1 v S ν e ϑ t ν 1 θ t ν
D k D operator
γ 0 ( z , θ ) = E z X 0 e t 0 θ
γ ( z , θ ) = E z X 1 e ( t 1 t 0 ) θ
I 0 modified Bessel function of order zero
A = k = 1 X k + a δ k ε t k soft shocks and degradation marked point process
A [ 0 , t ] = A ( t ) cumulative degradation process including shocks
A n = k = 1 n X k + a δ k
δ k = t k t k 1
ν = inf { n N : A n M }
Mfailure threshold
τ ν failure time if it occurs upon constant degradation
t ν failure time if it occurs upon a shock
γ ( α , θ ) = E e α ( X k + a δ k ) θ δ k
Φ ν = E e α A ν 1 β S ν ϑ t ν 1 θ t ν
L Laplace transform
L 1 inverse Laplace transform
pfailure threshold ( p > 0 )
Φ ( x ) standard Gaussian PDF
[ E x p ( λ ) ] equivalence class of all exponential distributions with parameter λ
[ G e o ( p ( x ) ) ] equivalence classs of all geometric random variables with parameter p ( x )
[ G a ( α , β ) ] equivalence class of all gamma distributions with parameters ( α , β )
[ N ( μ , σ 2 ) ] equivalence class of all Gaussian distributions with parameters ( μ , σ 2 )

Appendix A

Appendix A.1. Simulation Code

Simulation code was written in Python using the NumPy library. First, we created a function simulatePath that takes the parameters of the process as inputs, simulates the process until a failure occurs, and returns the four terms in the functional Φ ν , i.e., the failure index ν , damage upon the shock before the failure A ν 1 , the failure damage S ν , the time of the shock before the failure t ν 1 , and the failure time τ ν as well as a Boolean flag indicating a degradation failure (if TRUE) or shock failure (if FALSE).
Figure A1 is a flow chart describing how the simulation proceeds.
Figure A1. The simulatePath function accepts numerical inputs for each parameter of the process: a, λ , μ , and M. It then simulates one path of the process and returns sampled values of ν , A ν 1 , S ν , t ν 1 and τ ν along with a flag indicating whether a degradation failure or shock failure occurred.
Figure A1. The simulatePath function accepts numerical inputs for each parameter of the process: a, λ , μ , and M. It then simulates one path of the process and returns sampled values of ν , A ν 1 , S ν , t ν 1 and τ ν along with a flag indicating whether a degradation failure or shock failure occurred.
Axioms 10 00199 g0a1
Next is the full Python code for the simulatePath function.
  • import numpy as~np
  •  
  • def simulatePath(a, lam, mu, M):
  •     # initialize outputs
  •     failureTime = 0
  •     failureDamage = 0
  •     failureIndex = 0
  •     degradationFailure = False
  •  
  •     # simulate the process
  •     while failureDamage < M:
  •         # save A_j-1
  •         oldDamage = failureDamage
  •  
  •         # save t_j-1
  •         oldTime = failureTime
  •  
  •         # compute waiting time before the next shock
  •         waitingTime = np.random.exponential(1/mu)
  •  
  •         # add degradation between shocks
  •         failureDamage += a ∗ waitingTime
  •  
  •         # if degradation causes damage to reach M...
  •         if failureDamage >= M:
  •             # compute tau_nu
  •             failureTime += (M - oldDamage)/a
  •  
  •             # set S_nu (total damage) to M
  •             failureDamage = M
  •  
  •             # mark degradation as the cause of the failure
  •             degradationFailure = True
  •  
  •             # exit the loop
  •             break
  •  
  •         # else, add the shock
  •         else:
  •             # add the waiting time
  •             failureTime += waitingTime
  •  
  •             # add the shock damage
  •             failureDamage += np.random.exponential(1/lam)
  •  
  •             # add 1 to the shock counter
  •             failureIndex += 1
  •  
  •     # gather the output values into a tuple
  •     outputs = (failureIndex, oldDamage, failureDamage,
  •                oldTime, failureTime, degradationFailure)
  •  
  •     # return values nu, A_nu-1, S_nu, t_nu-1, tau_nu, flag for
  •     # failure type
  •     return outputs
Then, a single path of the process can be simulated with the command
  • simulatePath(a, lam, mu, M)
for any parameters one chooses. All empirical results in Section 4 are created by generating many paths of the deterioration process with this function and finding means and variances of the outputs across generated paths.
The formulas derived in Corollary 1, Proposition 2, and Proposition 4 were implemented as the following Python functions to compute predicted results.
  • def degradationProbability(a, lam, mu, M):
  •     term = a ∗ lam + mu ∗ np.exp(-(lam + mu / a) ∗ M)
  •     return term / (a ∗ lam + mu)
  •  
  • def failureTimeMean(a, lam, mu, M):
  •     term = lam ∗ M + 1 - degradationProbability(a, lam, mu, M)
  •     return term / (a ∗ lam + mu)
  •  
  • def failureTimeVariance(a, lam, mu, M):
  •     dTerm = a ∗ lam + mu
  •     eTerm = np.exp(-(lam + mu / a) ∗ M)
  •     term1 = 2 ∗ lam ∗ mu ∗ M / dTerm ∗∗ 3
  •     term2 = mu ∗ (mu - 4 ∗ a ∗ lam) / dTerm ∗∗ 4
  •     term3 = 2 ∗ mu
  •     term3 ∗= 2 ∗ a ∗∗ 2 ∗ lam + (a ∗ lam) ∗∗ 2 ∗ M - 2 ∗ mu ∗∗ 2 ∗ M
  •     term3 ∗= eTerm / (a ∗ dTerm ∗∗ 4)
  •     term4 = 2 ∗ mu ∗∗ 2 ∗ eTerm ∗∗ 2 / dTerm ∗∗ 4
  •     return term1 + term2 + term3 - term4
  •  
  • def failureDamageMean(a, lam, mu, M):
  •     return M + (1 - degradationProbability(a, lam, mu, M)) / lam
  •  
  • def failureDamageVariance(a, lam, mu, M):
  •     dTerm = a ∗ lam + mu
  •     eTerm = np.exp(-(lam + mu / a) ∗ M)
  •     term1 = 2 ∗ a ∗ mu ∗ (1 - eTerm) / (lam ∗ dTerm ∗∗ 2)
  •     term2 = mu ∗∗ 2 ∗ (1 - eTerm ∗∗ 2) / (lam ∗ dTerm) ∗∗ 2
  •     return term1 + term2

References

  1. Tsai, C.C.; Tseng, S.T.; Balakrishnan, N. Mis-specification analyses of gamma and Wiener degradation processes. J. Stat. Plan. Inference 2011, 141, 3725–3735. [Google Scholar] [CrossRef]
  2. Gut, A. Cumulative shock models. Adv. Appl. Probab. 1990, 22, 504–507. [Google Scholar] [CrossRef] [Green Version]
  3. Gut, A.; Hüsler, J. Shock models. In Advances in Degradation Modeling: Applications to Reliability, Survival Analysis, and Finance; Nikulin, M.S., Limnios, N., Balakrishnan, N., Kahle, W., Huber-Carol, C., Eds.; Springer: Berlin/Heidelberg, Germany, 2010; Chapter 5; pp. 56–124. [Google Scholar]
  4. Sumita, U.; Shanthikumar, J.G. A class of correlated cumulative shock models. Adv. Appl. Probab. 1985, 17, 347–366. [Google Scholar] [CrossRef]
  5. Abolnikov, L.; Dshalalow, J.H. A first passage problem and its applications to the analysis of a class of stochastic models. J. Appl. Math. Stoch. Anal. 1992, 5, 83–97. [Google Scholar] [CrossRef] [Green Version]
  6. Dshalalow, J.H.; Huang, W. On noncooperative hybrid stochastic games. Nonlinear Anal. Hybrid Syst. 2008, 2, 803–811. [Google Scholar] [CrossRef]
  7. Dshalalow, J.H.; Robinson, R. On one-sided stochastic games and their applications to finance. Stoch. Model. 2012, 28, 1–14. [Google Scholar] [CrossRef]
  8. Gut, A.; Hüsler, J. Extreme shock models. Extremes 1999, 2, 295–307. [Google Scholar] [CrossRef]
  9. Li, Z.; Kong, X. Life behavior of δ-shock model. Stat. Probab. Lett. 2007, 77, 577–587. [Google Scholar] [CrossRef]
  10. Eryilmaz, S.; Bayramoglu, K. Life behavior of δ-shock models for uniformly distributed interarrival times. Stat. Pap. 2013, 55, 841–852. [Google Scholar] [CrossRef]
  11. Eryılmaz, S. Generalized δ-shock model via runs. Stat. Probab. Lett. 2012, 82, 326–331. [Google Scholar] [CrossRef]
  12. Mallor, F.; Omey, E. Shocks, runs and random sums. J. Appl. Probab. 2001, 38, 438–448. [Google Scholar] [CrossRef]
  13. Mallor, F.; Santos, J. Reliability of systems subject to shocks with a stochastic dependence for the damages. Test 2003, 12, 427–444. [Google Scholar] [CrossRef]
  14. Mallor, F.; Omey, E.; Santos, J. Asymptotic results for a run and cumulative mixed shock model. J. Math. Sci. 2006, 138, 5410–5414. [Google Scholar] [CrossRef]
  15. Ozkut, M.; Eryilmaz, S. Reliability analysis under Marshall–Olkin run shock model. J. Comput. Appl. Math. 2019, 349, 52–59. [Google Scholar] [CrossRef]
  16. Marshall, A.W.; Olkin, I. A multivariate exponential distribution. J. Am. Stat. Assoc. 1967, 62, 30–44. [Google Scholar] [CrossRef]
  17. Marshall, A.W.; Olkin, I. A generalized bivariate exponential distribution. J. Appl. Probab. 1967, 4, 291–302. [Google Scholar] [CrossRef]
  18. Eryilmaz, S.; Tekin, M. Reliability evaluation of a system under a mixed shock model. J. Comput. Appl. Math. 2019, 352, 255–261. [Google Scholar] [CrossRef]
  19. Zhao, X.; Guo, X.; Wang, X. Reliability and maintenance policies for a two-stage shock model with self-healing mechanism. Reliab. Eng. Syst. Saf. 2018, 172, 185–194. [Google Scholar] [CrossRef]
  20. Peng, H.; Feng, Q.; Coit, D.W. Reliability and maintenance modeling for systems subject to multiple dependent competing failure processes. IIE Trans. 2010, 43, 12–22. [Google Scholar] [CrossRef]
  21. Hao, H.; Su, C.; Qu, Z. Reliability analysis for mechanical components subject to degradation process and random shock with Wiener process. In The 19th International Conference on Industrial Engineering and Engineering Management; Springer: Berlin/Heidelberg, Germany, 2013; pp. 531–543. [Google Scholar] [CrossRef]
  22. Rafiee, K.; Feng, Q.; Coit, D.W. Reliability assessment of competing risks with generalized mixed shock models. Reliab. Eng. Syst. Saf. 2017, 159, 1–11. [Google Scholar] [CrossRef] [Green Version]
  23. Hao, S.; Yang, J.; Ma, X.; Zhao, Y. Reliability modeling for mutually dependent competing failure processes due to degradation and random shocks. Appl. Math. Model. 2017, 51, 232–249. [Google Scholar] [CrossRef]
  24. Yousefi, N.; Coit, D.W.; Song, S.; Feng, Q. Optimization of on-condition thresholds for a system of degrading components with competing dependent failure processes. Reliab. Eng. Syst. Saf. 2019, 192, 106547. [Google Scholar] [CrossRef] [Green Version]
  25. Cao, Y.; Liu, S.; Fang, Z.; Dong, W. Modeling ageing effects in the context of continuous degradation and random shock. Comput. Ind. Eng. 2020, 145, 106539. [Google Scholar] [CrossRef]
  26. Oliveira, R.P.; Achcar, J.A.; Mazucheli, J.; Bertoli, W. A new class of bivariate Lindley distributions based on stress and shock models and some of their reliability properties. Reliab. Eng. Syst. Saf. 2021, 211, 107528. [Google Scholar] [CrossRef]
  27. Dshalalow, J.H.; White, R.T. Current trends in random walks on random lattices. Mathematics 2021, 9, 1148. [Google Scholar] [CrossRef]
  28. Dshalalow, J.H. On the level crossing of multi-dimensional delayed renewal processes. J. Appl. Math. Stoch. Anal. 1997, 10, 355–361. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The linear degradation process with shocks.
Figure 1. The linear degradation process with shocks.
Axioms 10 00199 g001
Figure 2. Predicted and empirical probabilities of degradation failures with limiting probabilities.
Figure 2. Predicted and empirical probabilities of degradation failures with limiting probabilities.
Axioms 10 00199 g002
Figure 3. Predicted and empirical values for the mean and standard deviation of the failure time.
Figure 3. Predicted and empirical values for the mean and standard deviation of the failure time.
Axioms 10 00199 g003
Figure 4. Predicted and empirical values for the mean and standard deviation of the failure damage.
Figure 4. Predicted and empirical values for the mean and standard deviation of the failure damage.
Axioms 10 00199 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Dshalalow, J.H.; White, R.T. Random Walk Analysis in a Reliability System under Constant Degradation and Random Shocks. Axioms 2021, 10, 199. https://doi.org/10.3390/axioms10030199

AMA Style

Dshalalow JH, White RT. Random Walk Analysis in a Reliability System under Constant Degradation and Random Shocks. Axioms. 2021; 10(3):199. https://doi.org/10.3390/axioms10030199

Chicago/Turabian Style

Dshalalow, Jewgeni H., and Ryan T. White. 2021. "Random Walk Analysis in a Reliability System under Constant Degradation and Random Shocks" Axioms 10, no. 3: 199. https://doi.org/10.3390/axioms10030199

APA Style

Dshalalow, J. H., & White, R. T. (2021). Random Walk Analysis in a Reliability System under Constant Degradation and Random Shocks. Axioms, 10(3), 199. https://doi.org/10.3390/axioms10030199

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop