趙 丹，周 巖，孟憲云，李海霞

（燕山大學(xué)理學(xué)院，河北秦皇島 066004）

壽命服從PH分布的冷貯備可修系統(tǒng)

趙 丹，周 巖，孟憲云，李海霞

（燕山大學(xué)理學(xué)院，河北秦皇島 066004）

研究了由兩個(gè)不同部件組成的冷貯備系統(tǒng)。假定部件的工作時(shí)間服從的是PH分布，且假定部件能夠修復(fù)如新，部件失效后的修理時(shí)間服從指數(shù)分布。利用引進(jìn)變量法、廣義的馬爾可夫過(guò)程和PH分布的性質(zhì)，得到了系統(tǒng)的瞬時(shí)概率分布、系統(tǒng)首次故障前的平均時(shí)間、系統(tǒng)穩(wěn)態(tài)概率分布等可靠性指標(biāo)的Laplace變換表達(dá)式。

PH分布；冷貯備；引進(jìn)變量法；廣義的馬爾可夫過(guò)程；可修系統(tǒng)；可靠性

0 引言

在可靠性分析中，常用的貯備系統(tǒng)有溫貯備系統(tǒng)和冷貯備系統(tǒng)，冷貯備系統(tǒng)應(yīng)用比較廣泛，因此受到學(xué)者們的關(guān)注。所謂冷貯備就是指部件在貯備過(guò)程中既不失效也不劣化，貯備期的長(zhǎng)短與部件日后使用時(shí)的工作壽命無(wú)關(guān)，貯備部件替換失效部件是通過(guò)轉(zhuǎn)換開(kāi)關(guān)來(lái)完成的。文獻(xiàn)［1］研究了壽命和修理時(shí)間服從一般分布的可修系統(tǒng)，文獻(xiàn)［2－3］分別研究了含有兩個(gè)同型和不同型部件的冷貯備可修系統(tǒng)，文獻(xiàn)［4］在此基礎(chǔ)上針對(duì)帶有優(yōu)先權(quán)的兩不同型部件冷貯備系統(tǒng)進(jìn)行了可靠性分析，文獻(xiàn)［5］研究了有優(yōu)先權(quán)且有兩不同修理設(shè)備的兩部件冷貯備可修系統(tǒng)，文獻(xiàn)［6］研究了修理時(shí)間服從PH分布、壽命服從指數(shù)分布的兩個(gè)不同型部件、一個(gè)修理工組成的并聯(lián)可修系統(tǒng)，更多關(guān)于帶有PH分布的可修系統(tǒng)的可靠性分析可參考文獻(xiàn)［7－9］。本文在以往文獻(xiàn)的研究基礎(chǔ)上，對(duì)于修理時(shí)間服從指數(shù)分布、壽命服從PH分布的兩個(gè)不同型部件、一個(gè)修理工組成的冷貯備系統(tǒng)進(jìn)行了討論。

1 模型假設(shè)

（Ⅰ）部件1的工作壽命服從的是m階PH分布，其分布函數(shù)為G1（t）＝1－αeTte，t≥0，過(guò)程的初始概率向量為α＝（α1，α2，…，αm），狀態(tài)轉(zhuǎn)移概率矩陣為T(mén)＝（Tij）m×m，且滿(mǎn)足T0＝－Te1＝（，…）′，其中e1為分量全為1的m階列向量。

部件2的工作壽命服從的是n階PH分布，其分布函數(shù)為G2（t）＝1－βeRte，t≥0，過(guò)程的初始概率向量為β＝（β1，β2，…，βn），狀態(tài)轉(zhuǎn)移概率矩陣為R＝（Rij）n×n，且滿(mǎn)足R0＝－Re2＝（，…，）′，其中e2為分量全為1的n階列向量。

（Ⅱ）部件i故障后的修理時(shí)間服從參數(shù)為μi的指數(shù)分布，其分布函數(shù)為Fi（t）＝1－e－μit，t≥0，（μi＞0），i＝1，2。

（Ⅲ）上述所有隨機(jī)變量是相互獨(dú)立的。

（Ⅳ）在t＝0時(shí)刻，部件1先開(kāi)始工作，部件2貯備，且兩個(gè)部件都是完好無(wú)損的。轉(zhuǎn)換開(kāi)關(guān)完全可靠且轉(zhuǎn)換是瞬時(shí)的。當(dāng)正在工作的部件癱瘓時(shí)，貯備部件立即去替換而進(jìn)入工作狀態(tài)，修理工立即修理失效部件；如果兩個(gè)部件都?jí)牧?，則遵循先壞先修的原則，假定修復(fù)如新。

2 主要結(jié)果及證明

X（t）＝i表示時(shí)刻t系統(tǒng)處于狀態(tài)i，i＝0，1，2，3，4，5。由于系統(tǒng)中所假設(shè)的部件壽命不服從指數(shù)分布，那么可以驗(yàn)證｛X（t），t≥0｝不是馬氏過(guò)程。現(xiàn)在引入補(bǔ)充變量：令Z（t）＝j(luò)表示正在工作中的部件在時(shí)刻t所處的位相，則｛X（t），Z（t），t≥0｝構(gòu)成了一個(gè)廣義的馬爾可夫過(guò)程。此時(shí)系統(tǒng)的狀態(tài)可有

證明由系統(tǒng)的模型假設(shè)可得系統(tǒng)的狀態(tài)為

狀態(tài)0：部件1開(kāi)始工作，部件2進(jìn)入貯備狀態(tài)。

狀態(tài)1：部件2開(kāi)始工作，部件1進(jìn)入貯備狀態(tài)。

狀態(tài)2：部件1修理，部件2開(kāi)始工作。

狀態(tài)3：部件2修理，部件1開(kāi)始工作。

狀態(tài)4：部件1正在修理，部件2處于待修。

狀態(tài)5：部件2正在修理，部件1處于待修。

證明由可用度的定義和定理1的結(jié)果即可得出結(jié)論。

定理2設(shè)Z為系統(tǒng)首次失效時(shí)間，則Z服從2（m＋n）階PH分布F（t），且可得系統(tǒng)的可靠度函數(shù)為R（t）＝P｛Z≥t｝＝γePte，e為分量全為1的2（m＋n）維列向量，且該函數(shù)的Lap lace變換為R*（s）＝

證明將上述馬氏過(guò)程中的故障狀態(tài)（4，5）歸并為一個(gè)吸收態(tài)N，則狀態(tài)空間為j＝｛（0，1），…，（0，m），（1，1），…，（1，n），（2，1），…，（2，n），（3，1），…，（3，m），N｝，對(duì)應(yīng)的轉(zhuǎn)移概率矩陣為：

劃去最后一行和最后一列的矩陣為P，其中，P0＝－Pe＝（0，0，R0，T0）T；γ＝（α，0，0，…，0）∈R2（m＋n）；e為分量全為1的2（m＋n）維列向量。則由PH分布的定義可知：系統(tǒng)的首次失效時(shí)間Z服從2（m＋n）階PH分布，再由可靠度的定義及Lap lace變換即可得到定理后半部的結(jié)論。

定理3設(shè)π＝（π0，π1，π2，π3，π4，π5）為系統(tǒng)的平穩(wěn)概率分布，其中πi是系統(tǒng)處于狀態(tài)i的穩(wěn)態(tài)概率，其中，π4，π5∈R＋；π1，π2是n維行向量；π0，π3是m維行向量，則有

證明該隨機(jī)過(guò)程是一個(gè)不可約有限狀態(tài)馬氏過(guò)程，所以其平穩(wěn)概率π＝（π0，π1，π2，π3，π4，π5）存在，且滿(mǎn)足方程組

將式（8）按方塊矩陣展開(kāi)為：

由式（13）和式（14）可得

將式（15）代入式（11）和式（12）化簡(jiǎn)得

將上式代入式（15）則有

由式（9）、式（10）和式（11）可得

推論2假設(shè)系統(tǒng)的穩(wěn)態(tài)可用度為A，則有

證明由可用度的定義和定理3的結(jié)果即可得出結(jié)論。

3 結(jié)論

本文討論了部件壽命服從PH分布、修理時(shí)間服從指數(shù)分布的兩不同型部件的冷貯備系統(tǒng)。PH分布是指數(shù)分布和幾何分布的矩陣類(lèi)比，它的參數(shù)矩陣具有明顯的概率意義，將之應(yīng)用于可靠性分析中，可以較易地得到系統(tǒng)的一些主要可靠性指標(biāo)，所得結(jié)果更具有一般性，因此更加具有實(shí)際意義。

［1］ 曹晉華，程侃.可靠性數(shù)學(xué)引論［M］.北京：科學(xué)出版社，2006.

［2］ 李才良，蒲冰遠(yuǎn)，唐應(yīng)輝，等.兩不同部件冷貯備系統(tǒng)的可靠性分析［J］.電子科技大學(xué)學(xué)報(bào)，2003，32（4）：447－450.

［3］ 周美秀，范偉峰，胡曉曉.兩相同部件冷貯備可修系統(tǒng)的定性分析［J］.數(shù)學(xué)的時(shí)間與認(rèn)識(shí)，2009，39（18）：50－53.

［4］ 張民悅，包林濤，段紅星.有優(yōu)先權(quán)的兩個(gè)不同型部件冷貯備可修系統(tǒng)［J］.蘭州理工大學(xué)學(xué)報(bào)，2007，33（2）：152－155.

［5］ 賈鵬茹，孟憲云，張曉爽，等.有優(yōu)先權(quán)且有兩不同修理設(shè)備的兩部件冷貯備可修系統(tǒng)的可靠性分析［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2010，4（10）：187－192.

［6］ 董兵，唐應(yīng)輝.修理時(shí)間服從PH分布的兩部件可靠性分析［J］.西華大學(xué)學(xué)報(bào)，2008，27（3）：69－71，76.

［7］ 岳德權(quán)，李蘭巧，祁洪娟.PH／M／c可修排隊(duì)系統(tǒng)［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2010，40（8）：132－137.

［8］ Delia M C，Rafael P O.A Deteriorating Two-system with Two Repair Modes and Sojourn Times Phase-type Distributed［J］.Reliability Engineering and System Safety，2006，91：1－9.

［9］ Juan E R C，Rafael P O，Gemma F V.Modelling a Reliability System Governed by Discrete Phase-type Distributions［J］. Reliability Engineering and System Safety，2008，93：1650－1657.

O213.2

A

1672－6871（2014）01－0092－05

全國(guó)統(tǒng)計(jì)科研計(jì)劃基金項(xiàng)目（2010LC33）；河北省教育廳計(jì)劃基金項(xiàng)目（2007323）

趙 丹（1986－），女，黑龍江鶴崗人，碩士生；周 巖（1964－），男，黑龍江齊齊哈爾人，副教授，碩士生導(dǎo)師，研究方向?yàn)橄到y(tǒng)的可靠性分析及優(yōu)化.

2012－04－20