arxiv: v1 [math.pr] 2 Sep 2014

Transcription

1 A random shock model wh mxed effec, ncludng compeng sof and sudden falures, and dependence Sophe MERCIER & Ha Ha PHAM arxv: v1 [mah.pr] 2 Sep 214 Preprn acceped for publcaon n Mehodology and Compung n Appled Probably, 31 Augus 214, hp://dx.do.org/1.17/s Absrac A sysem s consdered, whch s subjec o exernal and possbly faal shocks, wh dependence beween he faaly of a shock and he sysem age. Apar from hese shocks, he sysem suffers from compeng sof and sudden falures, where sof falures refer o he reachng of a gven hreshold for he degradaon level, and sudden falures o accdenal falures, characerzed by a falure rae. A non-faal shock ncreases boh degradaon level and falure rae of a random amoun, wh possble dependence beween he wo ncremens. The sysem relably s calculaed by four dfferen mehods. Condons under whch he sysem lfeme s New Beer han Used are proposed. The nfluence of varous parameers of he shocks envronmen on he sysem lfeme s suded. Keywords: Relably; Bvarae non homogeneous compound Posson process; Hazard rae process; Posson random measure; Sochasc order; Ageng properes; Two-componen seres sysem. AMS MSC: 6K1 & 6G51 1 Inroducon Ths paper s devoed o he survval analyss of a sysem subjec o compeng falure modes whn an exernal sressng envronmen. The exernal envronmen s assumed o sress he sysem a random and solaed mes accordng o a random shock model. Such a model can represen exernal demands e.g., whch pu some sress on he sysem a her arrvals. Shock models have been he subjec of an exensve leraure. Followng Mallor and Sanos 23a, shock models may be classfed no dfferen caegores, accordng o wheher arrval mes and shock magnudes are correlaed Mallor and Omey, 21; Mallor and Sanos, 23b, or ndependen Cha and Fnkelsen, 29, and Laboraore de Mahémaques e de leurs Applcaons Pau UMR CNRS 5142, Unversé de Pau e des Pays de l Adour, Bâmen IPRA, Avenue de l Unversé, F-6413 PAU cedex, FRANCE; sophe.mercer@unv-pau.fr Research and Developmen Insue, Duyan Unversy, K7/25 Quang Trung, Da Nang, VIETNAM; phamhaha9@gmal.com 1

2 accordng o he assumpon pu on he shocks arrval process: homogeneous or non-homogeneous Posson process A-Hameed and Proschan, 1973; Cha and Fnkelsen, 29; Cha and M, 27, 211; Esary e al., 1973; Qan e al., 1999, renewal process Skoulaks, 2 or non-saonary pure brh process A-Hameed and Proschan, Accordng o he nfluence of shocks on he sysem, shock models may be furher classfed no hree ypes: exreme shock models a shock can cause he sysem mmedae falure, see Gu and Hüsler 1999, cumulave shock models a shock ncreases some nrnsc characersc of he sysem such as s deeroraon, falure rae, age, number of already endured shocks,..., see Cha and M 27; Qan e al and mxed shock models a shock can eher cause he sysem mmedae falure or ncrease some nrnsc characersc, see Cha and M 211; Gu 21. More references can be found n Fnkelsen and Cha 213; Mallor and Sanos 23a; Nakagawa 27; Sngpurwalla In hs paper, a shock model wh mxed effecs s consdered, where he occurrence of shocks s classcally modelled hrough a non-homogeneous Posson process and where each shock may resul n he sysem mmedae falure hrough a Bernoull ral, ndependen of he sysem nrnsc behavour. Apar from shocks, he sysem suffers from compeng sof and sudden falures, where sof falures refer o he reachng of some gven hreshold for he degradaon level, and sudden falures o accdenal falures, characerzed by a falure rae. A possble represenaon corresponds o a wo-componen seres sysem, where he frs componen s subjec o sof falures and he second one o sudden falure. The sysem may hence fal hrough hree dfferen compeng modes: raumac falure due o a faal shock; sof falure; sudden falure. Each non faal shock nduces some ncrease of boh deeroraon level and falure, wh possble dependence n-beween. In he oldes leraure, mos models consdered one sngle possble ype of falures for he sysem: e.g. sof falures n Marshall and Shaked 1979 or raumac falures due o shocks n Savs 1988, boh n he mulvarae seng. More recenly, dfferen models have been developed, whch consder wo dfferen ypes of falures. For nsance, compeon beween sof and sudden falures s suded n Zhu e al. 21. Several applcaon cases are proposed n he paper see also references heren. An ndusral example s also provded n Wang and Gao 214, whch sudes he relably of an arcraf engne. Noe ha Zhu e al. 21 assumes sof and sudden falures o be ndependen whereas he sressng envronmen of he presen paper makes hem dependen. Compeon beween sof falures and occurrence of a raumac even are consdered n Degradaon-Threshold-Shock models DTS-models, denomnaon of Lehmann 26 whch have been proposed by Lemone and Wenocur 1985 and furher suded by Lehmann 26, 29. A case sudy s provded n Hao e al. 213 for he analyss of fague crack growh, where he effecs of shocks on he degradaon are pu no evdence. Noe ha, conrary o he presen paper, DTS-models consder some possble nfluence of he deeroraon level of he sysem on he shocks arrval rae. However, he frs shock of a DTS-model s always faal o he sysem leadng o a sngle shock possbly endured by he sysem, whereas successve non faal shocks are here envsoned. We could no fnd any paper whch akes no accoun hree compeng falure modes, as n he presen paper. However, based on he case sudes from he prevous leraure, one can hnk ha our model can reflec los of sysems subjec o compeng 2

3 sof and hard falures, whn a sressng envronmen. For nsance, one can hnk of an aeral cabn hooked o a cable: he supporng cable s deerorang due o corroson and fague; he lnkng pulley s subjec o sudden falures and maybe also o deeroraon; boh cable and pulley endure shocks a each cabn ravel, whch ncrease jonly her respecve deeroraon level and falure rae. The presen paper consder several knds of dependence beween he hree compeng falure modes: a each non faal shock, he ncrease of deeroraon and falure rae s smulaneous. Ths nduces a frs ype of dependence beween sof and sudden falures. Each shock may nduce a falure, eher because he shock s faal, or because he deeroraon s suddenly ncreased beyond he hreshold level. Ths nduces a second ype of dependence beween sof and raumac falures, whch may be smulaneous. Also, some possble dependence s envsoned beween he ncremens of falure rae and of deeroraon a each non faal shock. Ths nduces a hrd ype of dependence beween sof and sudden falures. Fnally, followng Cha and Fnkelsen 29; Cha and M 211, he probably for a shock o be faal depends on he shock arrval me, whch nduces a las ype of dependence. Up o our knowledge, all hese knds of dependence have no been ye consdered alogeher and, as wll be seen all along he ex, hs model enlarges several ones from he prevous leraure. The paper s organzed as follows: he model s specfed n Secon 2. The sysem relably s compued hrough dfferen mehods n Secon 3. Suffcen condons are provded n Secon 4 for he sysem lfeme o be New Beer han Used. The nfluence of varous parameers of he shock envronmen on he lfeme s suded n Secon 5. Numercal expermens are proposed n Secon 6 and concludng remarks end he paper n Secon 7. 2 The model To make he model clear, a wo-componen seres sysem s consdered, where he frs componen s subjec o sudden falure and he second one o sof falure. Ths wo-un sysem s jus a represenaon for he compeng sof and sudden falure modes, wh no resrcon. In he deal condon for example n a laboraory envronmen, he lfeme of he frs componen s characerzed by s nrnsc hazard rae h, whle he second one s subjec o some accumulave deeroraon modeled by an ncreasng sochasc process G e.g. a gamma process. The second componen fals once s deeroraon level exceeds a falure hreshold L. The lfemes of he wo componens are made dependen by her common sressng envronmen. Ths envronmen s modelled by a random shock process, where he shocks arrve accordng o a non-homogeneous Posson process N wh nensy dλ = λd or cumulaed nensy Λ. To avod useless echncal deals, we assume ha Λ > for all >. More generally, one mgh consder ha Λ > only for greaer han some >, whch would mean ha he shocks would arrve only afer me. The pons of he Posson process are denoed by T 1,..., T n,... wh T = < T 1 <... < T n <... almos surely. A shock a me may cause he sysem mmedae falure faal shock wh probably p [, 1], whch depends on he age of he sysem a he shock arrval. A shock a me s non faal wh probably q = 1 p. A non faal shock a me T ncreases he deeroraon of boh componens n a dfferen way: 3

4 for he frs componen, s hazard rae s ncreased of a non negave random amoun V 1, for he second one, s accumulaed deeroraon s ncreased of a non negave random amoun V 2. The random vecors V = V 1, V 2, = 1, 2,... are assumed o be ndependen and dencally dsrbued..d. wh common dsrbuon µdv 1, dv 2, and ndependen of he shocks arrval mes T n n 1 and hence ndependen of he Posson process N. A each shock, he ncremens V 1 and V 2 are possbly dependen. When subscrp s unnecessary, we drop and se V = V 1, V 2 o be a generc copy of V = V 1, V 2. For j = 1, 2, he dsrbuon of V j s denoed by µ j dv j. We se A o be he bvarae compound Posson process defned by N N A = V 1 = A 1, A 2 1 wh k=1 k, k=1 A 1 = A 2 = V 2 k N k=1 N k=1 V 1 k, V 2 k, where k=1 =. The processes A and G are assumed o be ndependen. Provded ha he sysem s funconng up o me, he random varables A 1 and A 2 sand for he cumulaed ncremens on [, ] of he falure rae of he frs componen and of he deeroraon of he second componen due o he exernal envronmen, respecvely. Seng F = σ A s, s o be he σ feld generaed by A s s and provded ha he sysem s sll up a me, he condonal hazard rae of he frs componen gven F s X 1 = h + A 1 and he condonal deeroraon of he second componen gven F s X 2 = G + A 2. To make clearer, we nroduce τ, = 1, 2 o be he lfeme of he h componen under he exernal envronmen, whou akng no accoun he possbly of faal shocks for he sysem. To smplfy he wrng, we denoe by P B F he condonal expecaon E 1 B F for any measurable se B, where 1 B sands for he ndcaor funcon 1 B ω = 1 f ω B, elsewhere. We hen have: P τ 1 > F = e X1 s ds = e H e A1 P τ 2 > F = P X 2 L F = P s ds, 2 G L A 2 F = F G L A 2 3 4

5 where H = hs ds 4 s he cumulaed nrnsc falure rae of he frs componen and where F G sands for he cumulave dsrbuon funcon c.d.f. of G. Recall ha G s ndependen of F. We now le τ 3 o be he me o he frs faal shock for he sysem wh τ 3 = nf n 1 : he shock a me T n s faal and we assume ha he Bernoull rals faal shocks or no whch happen a each shock arrval are ndependen one wh each oher, and ha hey depend on F only hrough he q T n s, ha s: P τ 3 > F = N q T 1 [,] T = q T 5 where = 1. The sysem falure s nduced eher by a faal shock or by a componen falure sof or sudden falure, whaever arrves frs. The lfeme of he sysem hence s τ = mn τ 1, τ 2, τ 3. 6 We fnally make he addonal assumpon ha τ 1, τ 2 and τ 3 are condonally ndependen gven F : P τ 1 >, τ 2 >, τ 3 > F = P τ 1 > F P τ 2 > F P τ 3 > F. To sum up, he whole model s specfed by: V 1, V 2 : he generc random ncremens n falure rae frs componen, V 1 and deeroraon second componen, V 2, λ x dx : he nensy of he non-homogeneous Posson process, h x : he nrnsc falure rae of he frs componen sudden falure, G : he nrnsc deeroraon of he second componen sof falure, p : he probably for a shock a me o be faal a he sysem level, plus some ndependence assumpons. By akng specal cases for hese fve ngredens, we can see ha our model exens some well-known models from he leraure. For nsance, akng p = no faal shocks, V 1 = consan, V 2 = and G = all no second componen, one ges he sochasc falure model n random envronmen from Cha and M 27. Takng V 2 = and G = all one sngle componen, one ges he sochasc survval model for a sysem under randomly varable envronmen from Cha and M

6 Table 1: A few parcular models from he leraure Brown-Proschan model from Brown and Proschan, 1983 h = V 1 = V 2 =, G =, Deermnsc boundary n Cha and Fnkelsen, 29 V 1 = h =, G = V 2 exponenally dsrbued Cha and M 27 q = 1, V 1 s a consan, V 2 =, G =, Cha and M 211 V 2 =, G =, Cumulave damage hreshold models q = 1, V 2 =, λ s a consan n Secon 4 of Esary e al., 1973 h =, G =, Qan e al q s a consan, h = V 1 =, G =, Takng G = all, V 1 = V 2 =, he model resumes o a classcal exreme shock model one sngle componen, where sysem falures are only due o shocks arrvng accordng o a non-homogeneous Posson process, wh probably p for a shock o be faal and q o be harmless. Ths model s nerpreed as he Brown-Proschan model by Cha and Fnkelsen 29; see also Brown and Proschan 1983 where varous properes of he model are explored. Takng p = no faal shocks, V 1 = h = one sngle componen, λ x = λ homogeneous Posson process, G = all no nrnsc deeroraon for he second - and sngle - componen, one ges he cumulave damage hreshold model from Esary e al., 1973, Secon 4, case of..d. damage ncremens. Takng p = p consan, V 1 = h = one sngle componen, G = all, one ges he cumulave damage model wh wo knds of shocks from Qan e al., 1999, case of..d. damage ncremens. Takng V 1 = h =, V 2 exponenally dsrbued, G =, all, one ges he model from Subsecon 3.b n Cha and Fnkelsen 29. All hese models are summed up n Table 1. Noe ha we do no preend a any exhausbly and our model wll nclude los of oher prevous models whch are no provded here. 3 Calculaon of he sysem relably The objecve of hs secon s o calculae he relably R L of he sysem a me, wh R L = Pτ >, all, where we recall ha he sysem lfeme τ s defned by 6. A frs way o compue R L s o use classcal Mone-Carlo smulaons and o smulae a large number of ndependen hsores for he sysem up o me Mehod 1. Ths mehod wll serve as a comparson ool n he numercal 6

7 expermenaons n Secon 6. Ths mehod requres he smulaon of a random varable wh condonal hazard rae h + N k=1 V 1 k see Algorhm 12 n Secon 6 for deals and may mply long compuaonal mes for he sysem relably. We provde below a few alernae mehods whch may be qucker and also easer o mplemen. Proposon 1 Mehod 2 The relably s gven by R L = Pτ > = e H φ L, 7 where H s provded by 4 and where φ L = E F G L A 2 e A1 s N = E F G L V 2 q T N ds e N V 1 T N qt Proof. Due o he condonal ndependence of τ 1, τ 2 and τ 3 gven F, we have: Usng 2, 3, 5, we ge: R L = Pτ 1 >, τ 2 >, τ 3 > R L = E = E P τ 1 >, τ 2 >, τ 3 > F = E P τ 1 > F P τ 2 > F P τ 3 > F. e H e A1 s ds F G L A 2 whch provdes 8 and nex 9, due o 1 and A 1 s ds = = V 1 1 {T s} ds = V 1 T 1 {T } = N V 1 q T, N 1 {T s} ds 8 9 V 1 T. 1 Based on he prevous resul, he only pon o ge he relably R L s o compue φ L. The remanng of he secon s hence devoed o he compuaon of φ L. Sarng from 8 or 9, a possbly s o compue φ L hrough Mone-Carlo smulaons of N and A, whch s smpler and qucker han smulang rajecores of he sysem accordng o he nal model. Ths mehod s called Mehod 2 n Secon 6. Followng A-Hameed and Proschan 1973 and Esary e al. 1973, one may also use a seres expanson of φ L, as provded by he followng proposon. Proposon 2 Mehod 3, general case φ L = e Λ n= P n, L Λn n! 11 7

8 where P n, L = E F G L n V 2 n qz e ZV 1 12 and Z 1 are..d. random varables wh probably densy funcon p.d.f. λx Λ 1 [,]x and ndependen of V 1. Proof. Condonng on he Posson process N, we have: φ L = E f N = + n= f n Λn e Λ n! wh N f n = E F G L V 2 N qt e N 1 TV N = n. Now, gven ha N = n, he condonal jon dsrbuon of T 1,... T n s he same as he jon dsrbuon of he order sascs Z 1,..., Z n of n..d. random varables Z 1,..., Z n wh p.d.f. λx Λ 1 [,]x see Cocozza-Thven 1998 e.g.. Usng he fac ha V = V 1, V 2 s ndependen of N, we ge: f n = E F G L n V 2 n qz e n Z V 1 Nong ha he expresson whn he expecaon s nvaran hrough permuaon of he Z s, we derve ha: n n f n = E F G L qz e n 1 ZV, whch provdes he resul. V 2 Remark 3 Based on he prevous resul, one can see ha our model s equvalen o a classcal shock model, where he shocks arrve accordng o a nonhomogeneous Posson process wh nensy dλx and wh condonal probably of survval a me equal o P n, L, gven ha here has been n shocks up o me. Corollary 4 Mehod 3, ndependen case In he specal case where V 1 and V 2 are ndependen, we ge:. wh φ L = e Λ Q n, L = E n= F G Q n, L an n! L n V

9 and a = µ 1 qλ = µ 1 zqλ z dz, 15 where sands for he convoluon operaor, qλx = qxλx, all x, and µ 1 sands for he Laplace ransform of he dsrbuon µ 1 of V 1, wh µ 1 s = + e xs µ 1 dx, all s. Proof. Sarng from 12 and usng he ndependence of all V 1 s, V 2 s and Z s, and he dencal dsrbuons of all V 1 s and of all Z s, we ge: n P n, L = E F G L E qz 1 e Z1V 1 n 16 wh V 2 E qz 1 e Z1V 1 = 1 Λ = 1 Λ λzq z E e zv 1 dz λq z µ 1 z dz = a Λ. Subsung hs expresson no Eq. 16 and nex no Eq. 12 provdes he resul. Example 5 Le V 2 be an exponenally dsrbued wh mean 1/θ and V 1 an ndependen random varable. In ha case, n V 2 s Gamma dsrbued wh parameer n, θ. Ths provdes Q n, L = L F G L x θn x n 1 n 1! e θx dx, all n 1. Usng Q, L = F G L and Eq. 13, we ge: R L = e H Λ F G L + a n n=1 n! L F G L x θn x n 1 n 1! e θx dx. Takng G = for all, V 1 = and h = as a specal case, we ge and, for L : R L = e Λ 1 + = e Λ 1 + F G L = 1 { L}, n=1 n=1 a = qλ s ds qλ s L dsn θ n x n 1 n! n 1! e θx dx qλ s dsn + θ L k e θl n! k! k=n usng successve negraons by pars for he las negral. Ths las expresson s he resul of Theorem 2 n Cha and Fnkelsen 29, whch hence appears as a specal case of he prevous resuls. 9

10 From Proposon 2 and Corollary 4, one may derve he followng approxmaon for φ L. Corollary 6 Mehod 3, approxmaon For N, le φ N L be defned by n he general case, and by φ N L = e Λ φ N L = e Λ N n= N n= P n, L Λn n! Q n, L an n! n case V 1 and V 2 are ndependen, where P n, L and Q n, L are provded by 12 and 14, respecvely. Then, for all, he sequence φ N L N 1 ncreases o he lm φ L when N and for all N, we have where ɛ N = e Λ φ N L φ L φ N L + ɛ N n=n+1 a n = e Λ a PY > N, n! a s defned by 15 and Y s Posson dsrbued wh mean a. Proof. We jus look a he general case. From Proposon 2, we have φ L φ N L = e Λ Due o F G L n V 2 1, we ge: n P n, L E qz e n=n+1 ZV 1 based on he proof of Corollary 4. Ths provdes: φ L φ N L e Λ n=n+1 P n, L Λn. n! = n a Λ a n = e Λ a n! n=n+1 a a n e n! and he resul. The prevous proposon provdes numercal bounds for φ L, whch may be adjused as gh as necessary, akng N large enough. Also, he requred number of erms n he runcaed seres s gven, o ge a specfed precson. Ths mehod s que adaped as soon as s possble o compue he P n, L s o be known n full form, whch s he case e.g. when he V 2 s are consan or Gamma dsrbued. An example s provded n Example 5, n he specal case of an exponenal dsrbuon. In he mos general case, he compuaon of he P n, L s or of he Q n, L s may be as dffcul as he nal problem of compung φ L, so ha he prevous mehod s no always adaped. We fnally presen anoher mehod based on Laplace ransform, whch does no suffer from he same resrcon. or he Q n, L s. Ths mosly requres he dsrbuon of n V 2 1

11 Theorem 7 Mehod 4 We have: or equvalenly where ν s provded by s Laplace ransform φ s = F G s ν s, all s 17 φ L = F G ν L, all L, 18 ν s = e Λ+qλ µ,s, 19 wh µ he bvarae Laplace ransform of he dsrbuon µ of V 1, V 2 : µu, s = e uv1 sv2 µdv 1, dv 2, all u, s, R 2 + and µ, s : u µu, s, all s. In he specal case where V 1 and V 2 are ndependen, ν s may be smplfed no: ν s = e Λ+a µ2s where a s provded by 15 and where µ 2 s s he unvarae Laplace ransform of µ 2. Proof. Rememberng ha A 2 [ φ s = e sl E F G = N L A 2 [ = E e sl F G L A 2 V 2 e N V 1, we ge from 9 ha: N T dl e N V 1 ] qt dl N T ] qt 2 wh e sl F G L A 2 dl = A 2 e sl F G L A 2 dl because G s non negave. Seng w = L A 2 A 2 e sl F G L A 2 dl = e sa2 Subsung hs expresson no 2 provdes:, we oban e sw F G w dw = e sa2 FG s. φ s = F G s θ s 21 wh θ s = E [ e s N V 2 e N V 1 T ] N qt [ = E e ] N V 1 T +sv 2 ln qt. 11

12 Because N s equvalen o T, we ge: θ s = E e 1 ψs,v,v 2,T wh ψ s, v 1, v 2, w = wv 1 + sv 2 ln qw 1 {w }. 22 Nong ha he sequence V n 1, V n 2, T n are he pons of a Posson random measure M wh nensy ν dv 1, dv 2, dw = µdv 1, dv 2 λwdw, he funcon θ s may be nerpreed as a Laplace funconal wh respec of M : n θ s = E e Mψs,. The formula for Laplace funconals of Posson random measures Çnlar, 211, Theorem 2.9 nex provdes: θ s = exp 1 e ψ s, dν 23 wh R e ψ s, dν = Subsung ψ s, by s expresson 22, we ge: 1 e ψ s, dν R 3 + R e ψs,v1,v2,w µdv 1, dv 2 λwdw. = = R 3 + = Λ R q w 1 e wv1+sv2 ln qw µdv 1, dv 2 λw dw R 2 + e wv1+sv2 µdv 1, dv 2 qλ w µ w, s dw = Λ [qλ µ, s]. λw dw Subsung hs expresson no 23 and nex no 21 provdes φ s = F G s ν s, wh ν s gven by 19. Equaon 18 s a drec consequence. Fnally, n case V 1 and V 2 1 are ndependen, we have: and µw, s = µ 1 w µ 2 s ν s = e Λ+qλ µ1 µ2s = e Λ+a µ2s, 24 whch ends hs proof. Based on he prevous resul, one can compue φ L by nverng s Laplace ransform wh respec of L. Lookng a Equaon 18, he key pon s he nverson of he Laplace ransform ν s. We nex provde an example where he nverson s possble n full form. In he mos general case, hs can be done numercally usng some Laplace nverson sofware. 12

13 Example 8 Le h =, G = all, λ and q consan, V 1 = V 2 dencally exponenally dsrbued wh mean 1/θ so ha V 1 and V 2 are compleely dependen. Then: µu, s = R 2 + e uv1 sv2 θe θv1 dv 1 δ v1 dv 2 = θ u + s + θ, where δ v1 sands for he Drac mass a v 1. We easly ge: θ + s + θ qλ µ, s = qλ du = qλθ ln u + s + θ s + θ and For < θ, we have where ν s = e λ 1 + qλθ. s + θ ν s = e λ 1 + qλθ = n n= n+1 qλθ n + 1 s + θ qλθqλθ 1... qλθ n + 1. n! Inverng he Laplace ransform ν s, we oban: qλθ ν dx = e δ λ n+1 dx + e θx x n dx. n + 1 n! n= As F G = 1, we ge he followng full form for he relably: R L = 1 ν L n+1 qλθ = e 1 λ + F n+1,θl, n + 1 θ n= where F n+1,θ s he cumulave dsrbuon funcon of a gamma dsrbued random varable wh parameer n + 1, θ. In he specal case where V 1 and V 2 are ndependen, he Laplace nverson of ν s s reduced o nverng e a µ2s, or equvalenly o nverng e C µ2s, where C s a consan. Ths s hence easer han n he mos general case of correlaed V 1 and V 2. To sum up he secon, we have a our dsposal four dfferen mehods for compung he relably: Mehod 1 Drec MC smulaons The man drawbacks of hs mehod are ha suffers from long compuaon mes and ha s mplemenaon s less drec han for he oher mehods. Mehod 2 Compung φ L hrough formula 9 and MC smulaons Ths mehod s much qucker and much easer o mplemen han Mehod 1. Besdes, s always possble o use. However, Mehod 3 when possble and Mehod 4 are qucker. 13

14 Mehod 3 Truncaed seres expanson + conrol of he runcaon error hrough Corollary 6 Ths mehod provdes very good resuls as soon as he P n, L s or he Q n, L s are avalable n full form. Mehod 4 Laplace ransform nverson Ths mehod s he bes when s possble o nverse he Laplace ransform ν s n full form. Numercal Laplace nverson also provdes que good resuls. 4 An ageng propery for he sysem lfeme Le us recall ha a random varable Z or F Z s New Beer han Used NBU f PZ > s + PZ > spz >, 25 all s,. We here provde suffcen condons under whch τ s NBU. Theorem 9 Assume ha he nrnsc lfemes of boh componens are NBU, whch means ha: e Hs+ e Hs e H, all s,, 26 F G+s l F G l F Gs l, all l, s,, 27 where he second condon s rue as soon as G s a unvarae non negave Lévy process. Then, τ s NBU f one among he wo followng condons s sasfed: 1. q s non ncreasng and λ s consan, 2. q s consan and Λ s super-addve Λx+y Λx+Λy, all x, y. Proof. Le us frs noe ha, n case G s a unvarae non negave Lévy process, we have: F G+s u = P G + G +s G u P G u, G +s G u = P G u P G +s G u = F G u F Gs u due o he ndependen and homogenous ncremens of G for he hrd lne. Assumpon 27 s hence rue. Sarng agan from Pτ > = e H φ L, and based on he NBU assumpon 26, s suffcen o show ha φ L s NBU. Now: φ s+ L = E E F G+s L F G L N +s N +s V 2 V 2 e N +s V 1 +s T F Gs L N +s V 2 N +s qt e N +s V 1 +s T N +s qt 14

15 due o he second NBU assumpon 27. Under each of he wo provded condons, q s non ncreasng so ha qt qt, all N + 1. Usng N V 1 N + s T N +s L N +s L V 2 V 2 N L L V 1 T, N +s V 2, =N +1 V 2, and splng he exponenal and he produc no wo pars, one ges: N N φ s+ L E [F G L T qt Seng hen T n F Gs L n 1 N +s =N +1 T V 2 V 2 e N V 1 e N +s =N +1 V 1 +s T N +s =N +1 = T N+, all 1, are pons of he Posson process N s = N +s N adms λ + x dx for nensy. Equaon 28 now wres: N φ s+ L E [F G L e N N V 1 T qt or equvalenly: As φ s+ L N s F Gs L + n= s V 2 N V 2 j+n j=1 N E [1 {N=n}F G L F Gs L s s ndependen on N u u and as he V s are..d. and ndependen on N u u, one ges: s N V 2 j+n j=1 e N s j=1 V 1 j+n s T j e N s V 2 j=1 V 1 j+n s T j N s qt j j=1 e N V 1 T N s qt j j=1 qt. s, N qt. 28 wh φ s+ L 15 + n=1 a n b n

16 wh N a n = E [1 {N=n}F G L N s b n = E F Gs L V 2 j+n j=1 V 2 e N s e N V 1 T j=1 V 1 j+n s T j N s ] N qt qt j j=1. Nong ha b n s ndependen on n b n = b, all n and ha + n=1 a n = φ L, we fnally have where L = b = E F Gs φ s φ s+ L φ L φ s L N s L j=1 V 2 j e N s j=1 V 1 j s T j N s qt j j=1 The pon now s o prove ha φ s L φ s L under he wo dfferen assumpons. 1. If λ s a consan, hen We hence have and he resul s clear. 2. If q s consan, hen: N s s φ s L = E where f n n = E F Gs L. s dencally dsrbued as N u u. φ s L = φ s L n j=1 [f V 2 j T 1 {j s} ], e n j=1 V 1 j s j1 {j s} q n j=1 1 { j s} for n N { }. Moreover, he respecve cumulaed nenses of N u u and are Λx and Λx + Λ, wh N u u Λx + Λ Λx, all, x due o he super-addvy of Λx. We derve from Shaked and Shanhkumar, 26, Theorem 6.B.4, Example 6.B.41 ha n T n so T for all n 1, where so sand for he sandard sochasc order. As f n s non decreasng wh respec o each, we ge ha: [ n ] E [f n T n ] E f n T 16

17 for each n N. Seng n +, we derve by Lebesgue s domnaed convergence heorem ha [ n ] lm E [f n T n n + ] = φ s L lm E f n T n + = φ s L, whch acheves he proof. The condons of he prevous heorem means ha 1. he probably for a shock o be non faal decreases wh me and λ s consan, 2. he cumulaed rae of shocks arrvals s larger a me han a me = and q s consan. Such condons hence mean ha he envronmen s more and more sressng, or ha s more sressng afer a whle han a he begnnng. Such condons are que naural. 5 Influence of he dependence nduced by he sressng envronmen on he sysem lfeme We here sudy he nfluence on he lfeme τ of dfferen parameers of he sressng envronmen: we sudy he nfluence of probably q, of he dependence beween V 1 and V 2 and of he cumulaed nensy funcon Λ. We hence look a he nfluence on he lfeme τ of all characerscs of he sressng envronmen whch make he componens dependen. The nfluence of q s sraghforward. We menon for sake of compleeness. 5.1 Influence of q on he lfeme τ Le us consder wo dfferen sysems wh dencal parameers excep from q frs sysem and q second sysem and such ha qw qw for all w. Then, addng a lde o any quany referrng o he second sysem, we drecly ge from 8 ha R L R L, all, or equvalenly ha τ s smaller han τ n he sense of he sandard sochasc order τ s τ : P τ > P τ >, all. 29 As expeced, he lfeme τ s hence ncreasng wh he probably q for a shock o be non faal. 5.2 Influence of he dependence beween V 1 and V 2 on he lfeme τ We here sudy he nfluence of he dependence beween he wo margnal ncremens V 1 and V 2 on he lfeme τ. To measure he dependence level 17

18 beween V 1 and V 2, we use he lower or upper orhan order, where we recall ha V = V 1, V 2 s sad o be smaller han Ṽ Ṽ = 1, Ṽ 2 n he lower orhan order V lo Ṽ f P V 1 x 1, V 2 x 2 P Ṽ 1 x 1, Ṽ 2 x 2, all x 1, x 2 R, 3 or equvalenly f P V 1 > x 1, V 2 > x 2 P Ṽ 1 > x 1, Ṽ 2 > x 2, all x 1, x 2 R. 31 Proposon 1 Le us consder wo dfferen sysems, wh dencal parameers excep from V 1, V 2 frs sysem and Ṽ 1, Ṽ 2 second sysem. As prevously, a lde s added o any quany referrng o he second sysem. Assume ha V 1, V 2 lo Ṽ 1, Ṽ 2. Then τ s smaller han τ n he sense of he sandard sochasc order τ s τ. Proof. The pon s o show ha φ L φ L. Sarng from 9 and condonng on σ N, we have where Θn, 1,..., n N = E F G L = E = E F G k 1 n L n V 2 V 2 V 2 φ L = E ΘN, T 1,..., T N N n k 2 n qt e N V 1 T N = n, T 1 = 1,..., T n = n q e n V 1 V 1 for all n N, all 1... n, wh k 1 v 1 = F G L v 1, n k 2 v 2 = q e v2. Based on he ndependence beween V 1, V 2 and N, s suffcen o show ha Θn, 1,..., n Θn, 1,..., n, all n N, all 1... n o ge φ L φ L. As x 1, x 2 s non decreasng n x 1 and x 2, we frs know from Shaked and Shanhkumar, 26, Theorem 6.G.3. ha V 1, V 2 lo Ṽ 1, Ṽ 2, 18

19 all 1 n. As V 1, V 2 and Ṽ 1, Ṽ 2 are wo,n,n sequences of ndependen random vecors, we derve from he same heorem ha n V 1 T, n V 2 lo n Ṽ 1 T, n Ṽ 2 As boh funcons k 1 and k 2 are non decreasng, we derve same heorem ha Θn, 1,..., n Θn, 1,..., n, whch acheves hs proof. The prevous resul shows ha he more dependen V 1 and V 2 are, he larger he sysem lfeme s. 5.3 Influence of he cumulaed nensy funcon Λ on he lfeme τ We fnally sudy how he lfeme of he sysem depends on he frequency of shocks. Proposon 11 Le us consder wo dfferen sysems, wh dencal parameers excep from Λ frs sysem and Λ second sysem. As prevously, a lde s added o any quany referrng o he second sysem. Assume ha Λ Λ, and ha q s non decreasng. Then τ s smaller han τ n he sense of he sandard sochasc order τ s τ. Proof. Usng a smlar mehod as for he proof of Theorem 9, we can wre φ L as φ L = E g T where g n n = E F G L n V 2 1 { } n lnq V e 1 1 { } s non decreasng wh respec o each, all 1 n. As Λ Λ, we derve n he same way ha n T n so T for all n 1, whch allows o conclude. Ths resul s very naural. The more frequen he shocks occur, he shorer he lfeme s. 6 Numercal expermens 6.1 Valdaon of he resuls As already menoned n Subsecon 3, a frs possbly for compung he sysem relably R L s o use classcal Mone-Carlo MC smulaons Mehod. 19

20 1 and smulae a large number of ndependen hsores for he sysem up o me. We here provde he algorhm ha we have used, consderng he case where h =, all and where V 1 s almos surely posve P V 1 > = 1 no essenal assumpons. Algorhm 12 Repea M mes wh M large enough: 1. Smulae G wh gven dsrbuon. 2. Smulae N accordng o he Posson dsrbuon wh parameer Λ. 3. Smulae N..d. random varables W 1,..., W n wh p.d.f. The shock arrval mes are gven by λx Λ 1 [,]x. T 1,..., T N = W 1,..., W N, where W 1,..., W N s he order sascs of W 1,..., W N. 4. Smulae Z as he resul of a Bernoull ral beween a faal and a non faal shock 1 a me T, wh probably qt for a shock o be non faal, all = 1,..., N. 5. Smulae N..d. random vecors V 1, V 2, = 1,..., N accordng o dsrbuon µ. 6. Smulae he lfeme Y of he frs componen wh condonal hazard rae gven F = σ A s, s. Wh ha am, seng A 1 κ = A 1 s ds = N T V 1 see 1, κ s a one-o-one ncreasng funcon from [T 1, + no [, + and, for κ T j u < κ T j+1 wh j 1, we have: κ 1 u = u + j T V 1 j V. 1 I s hen known ha, f U s unformly dsrbued on [, 1], hen κ 1 ln U s dencally dsrbued as Y, see Cocozza-Thven, 1998, Proposon Compue w j = 1 {Y >} 1 {G+ N N V Z 2 L} where j refers o he j h MC hsory, wh 1 j M. A he end of he algorhm, symbol w j sands for he realzaon of a Bernoull ral W beween an up 1 or down sysem a me, wh probably R L for he sysem o be up. The relably R L s hen classcally approxmaed by he emprcal mean m W of he w j s and a 95% asympoc confdence nerval s compued. 2

21 Table 2: Valdaon of he resuls Inpu V 1, V 2 Mehod R1 95 % CI qx = e x [ ] V 1 E1, V 2 E1 Independence [ ] 3, qx = [ ] V 1 = V 2 E1 Complee dependence [ ] qx = e x,v 2 = V 1 + W [ ] V 1, W 2 ndependen Dependence [ ] V 1 E1, W 2 E1 Table 3: Parameers for he examples L h G q λ V 1 V 2 V 1, V 2 Ex E1 E1 Independen Ex E1 E1 - Ex e x - 1 E1 Independen Mehod 2 s based on MC smulaons of rajecores of A and of N see Secon 3. For boh mehods 1 & 2, MC smulaons are based on N = 1 5 hsores. Mehods 3 & 4 are descrbed n Secon 3. The four mehods are compared on a few specfc examples. In all hese examples, we compare he relably a me = 1 R1 and we suppose ha he shock are due o a homogeneous Posson process wh parameer λ = 1, L = 2, h = and G s a null process. All oher parameers are provded n Table 2, where T E1 means ha he random varable T s exponenally dsrbued wh mean 1. Mehods 1 and 2 may be used n any case. Mehod 2 s more effecve han Mehod 1 shorer c.p.u. me and gher 95% confdence nerval - IC -. Mehod 3 s more praccal when V 1 and V 2 are ndependen wh some specfc dsrbuon. Mehod 4 s adaped o he case where V 1 and V 2 are dependen. 6.2 Examples We here llusrae several properes from a numercal pon of vew on a few examples. Examples parameers are provded n Table 3. Example 13 Ths example llusraes he NBU propery of he lfeme when λ s consan and qx = e x s non ncreasng, see Theorem 9. Takng L = 2, Fg. 1 ndeed shows ha he remanng lfeme of a sysem wh age = 1 s sochascally smaller han he lfeme of a new sysem. On he conrary, when 21

22 1.8 Pτ> Pτ>+ /Pτ> Relably Fgure 1: Example 13, NBU case λ s sll consan bu qx = 1 e 3x s non decreasng, he remanng lfeme of a sysem wh age = 1 s no comparable wh ha of a new sysem, see Fg. 2 wh L = 4. So he NBU propery does no hold anymore n ha case. Example 14 We here consder wo exreme cases for he dependency beween V 1 and V 2 : ndependen or compleely dependen here V 1 = V 2. The relably n he compleely dependen case s always greaer han n he ndependen case Fg. 3. Ths resul s coheren wh Proposon 1. Example 15 Ths example shows he monoony of he relably wh respec o he nensy of he Posson process λ x when q s ncreasng, see Proposon 11. The more frequenly he shocks occur, he lower he relably s Fg Concludng remarks We here proposed a random shock model wh compeng falure modes, whch enlarges several models from he prevous leraure. The model akes no accoun dfferen ypes of dependence beween compeng falures modes, where he dependence s nduced by a common exernal shock envronmen. The relably has been calculaed by several dfferen mehods and condons have been provded under whch he sysem lfeme s New Beer han Used. Due o hs ageng propery, mgh be of neres o propose and sudy some manenance polcy o enlarge he sysem lfeme. Several versons mgh be proposed, accordng o he avalable nformaon. Also, he nfluence of he characerscs of he sressng envronmen on he lfeme τ has been suded. As expeced, we saw ha he lfeme was sochascally ncreasng wh he probably q for a shock o be non faal. Besdes, and ha resul was no necessarly so clear a frs sgh, we saw ha he lfeme was also sochascally ncreasng wh he dependence beween he wo margnal shock szes. Fnally, n case of a non decreasng funcon q, we saw ha he lfeme was sochascally decreasng wh he cumulaed frequency 22

23 1.8 Pτ> Pτ>+ /Pτ> Relably Fgure 2: Example 13, no NBU case 1.8 complee dependence ndependence Relably Fgure 3: Comparson of relably for wo dfferen ypes of dependence beween V 1 and V 2, Example 14 23

24 1.8 λx=x λx=2x Pτ> Fgure 4: Comparson of relably for wo dfferen λ x, Example 15 of shocks. Ths means ha he more frequen he shocks occur, he shorer he lfeme s. Ths resul s naural bu our proof s lmed o he case of a non decreasng funcon q. In he specal case of Cha and M 211, he survval funcon of τ s however gven by Pτ > = e H 1 µ1 wqwλw dw and s easy o check ha f λ λ sronger assumpon han Λ Λ hen τ s τ, whou any specal condon on q. So, he sochasc monooncy of he lfeme wh respec of he cumulaed? frequency of shocks mgh be rue under a more general seng han n he presen paper, whou assumng any monooncy condon on q. We however have no been able o conclude on hs pon, and wheher s rue or no remans an open queson. Acknowledgemen 16 Boh auhors hank he referees for her carefull readng of he paper and her consrucve remarks, whch lead o a beer nroducon and jusfcaon of he model, and o a clearer paper. Ths work has been naed durng Ha Ha PHAM s PhD sudes n Pau France, and has been suppored by he Consel Régonal d Aquane France. Ths work has also been suppored for boh auhors by he French Naonal Research Agency AMMSI projec, ref. ANR 211 BS1-21. References M. S. A-Hameed and F. Proschan. Nonsaonary shock models. Sochasc Processes and her Applcaons, 14:383 44, M. S. A-Hameed and F. Proschan. Shock models wh underlyng brh process. Journal of Appled Probably, 121:18 28, M. Brown and F. Proschan. Imperfec repar. Journal of Appled Probably, 24: ,

25 J. H. Cha and M. Fnkelsen. On a ermnang shock process wh ndependen wear ncremens. Journal of Appled Probably, 462: , 29. J. H. Cha and J. M. Sudy of a sochasc falure model n a random envronmen. Journal of Appled Probably, 441: , 27. J. H. Cha and J. M. On a sochasc survval model for a sysem under randomly varable envronmen. Mehodology and Compung n Appled Probably, 13 3: , 211. Erhan Çnlar. Probably and sochascs, volume 261 of Graduae exs n Mahemacs. Sprnger Scence + Busness Meda, 211. C. Cocozza-Thven. Processus sochasques e fablé des sysèmes, volume 28 of Mahémaques e Applcaons. Sprnger, J. D. Esary, A. W. Marshall, and F. Proschan. Shock models and wear processes. The Annals of Probably, 14: , M. Fnkelsen and J. H. Cha. Sochasc modelng for relably: Shocks, burnn and heerogeneous populaons. Sprnger Seres n Relably Engneerng. Sprnger, London, 213. A. Gu. Mxed shock models. Bernoull, 73: , 21. A. Gu and J. Hüsler. Exreme shock models. Exremes, 23:295 37, H.-B. Hao, C. Su, and Z.-Z. Qu. Relably analyss for mechancal componens subjec o degradaon process and random shock wh wener process. In 19h Inernaonal Conference on Indusral Engneerng and Engneerng Managemen, pages , 213. A. Lehmann. Degradaon-hreshold-shock models. Probably, Sascs and Modellng n Publc Healh, pages , 26. A. Lehmann. Jon modelng of degradaon and falure me daa. Journal of Sascal Plannng and Inference, 1395: , 29. Specal Issue on Degradaon, Damage, Fague and Acceleraed Lfe Models n Relably Tesng. A. J. Lemone and M. L. Wenocur. On falure modelng. Naval Research Logscs Quarerly, 323:497 58, F. Mallor and E. Omey. Shocks, runs and random sums. Journal of Appled Probably, 382: , 21. F. Mallor and J. Sanos. Classfcaon of shock models n sysem relably. Monografías del Semn. Maem. García de Galdeano, 27:45 412, 23a. F. Mallor and J. Sanos. Relably of sysems subjec o shocks wh a sochasc dependence for he damages. Tes, 122: , 23b. A.W. Marshall and M. Shaked. Mulvarae shock models for dsrbuons wh ncreasng hazard rae average. The Annals of Probably, 72: ,

26 T. Nakagawa. Shock and damage models n relably heory. Sprnger Seres n Relably Engneerng. Sprnger, London, 27. C. Qan, S. Nakamura, and T. Nakagawa. Cumulave damage model wh wo knds of shocks and s applcaon o he backup polcy. Journal of he Operaons Research Socey of Japan-Kee Kagaku, 424:51 511, T. H. Savs. Some mulvarae dsrbuons derved from a non-faal shock model. Journal of Appled Probably, 252:383 39, M. Shaked and J. G. Shanhkumar. Sochasc orders. Sprnger Seres n Sascs. Sprnger, 26. N. D. Sngpurwalla. Survval n dynamc envronmens. Sascal Scence, 1 1:86 13, G. Skoulaks. A general shock model for a relably sysem. Journal of Appled Probably, 374: , 2. H.W. Wang and J. Gao. A relably evaluaon sudy based on compeng falures for arcraf engnes. Eksploaacja Nezawodnosc: Manenance and Relably, 162: , 214. Y. Zhu, E.A. Elsayed, H. Lao, and L.Y. Chan. Avalably opmzaon of sysems subjec o compeng rsk. European Journal of Operaonal Research, 223: ,