(海軍工程大學(xué)兵器工程系 武漢 430033)
備件配置優(yōu)化,即為裝備確定備件品種和數(shù)量最合理的配置,以在有限的資源約束條件下實(shí)現(xiàn)最大的備件供應(yīng)效能,或者以最小的資源實(shí)現(xiàn)要求的供應(yīng)效能水平[1~2]。對于工程中由多部件組成的系統(tǒng)而言,為了確保其連續(xù)工作能力,通常會對其配置一定數(shù)量的備件。在配置過程中,為了使系統(tǒng)具備較高的連續(xù)工作能力,往往需要配置較多的備件,這就產(chǎn)生了高額的保障費(fèi)用(如采購和庫存費(fèi)用)。而配置備件較少,又會頻繁出現(xiàn)備件短缺的情況,導(dǎo)致系統(tǒng)不可用。因此,在有限經(jīng)費(fèi)條件下,如何合理地配置備件以滿足系統(tǒng)的使用要求,具有重要的經(jīng)濟(jì)意義和使用價值。
針對由不同壽命分布類型不可修部件組成的系統(tǒng),國內(nèi)外學(xué)者開展了大量的相關(guān)研究。通過綜合考慮系統(tǒng)的戰(zhàn)備完好性或任務(wù)可靠度,以及保障費(fèi)用等條件建立了相應(yīng)的備件配置模型,并且通過采用邊際效益、動態(tài)規(guī)劃或遺傳算法等優(yōu)化算法對備件配置方案進(jìn)行了優(yōu)化。如文獻(xiàn)[1]針對部件服從任意壽命分布的多部件系統(tǒng),提出了基于蒙特卡羅仿真-邊際效益分析的備件配置優(yōu)化算法。文獻(xiàn)[3]綜合考慮備件的費(fèi)用、體積和質(zhì)量等因素,提出了系統(tǒng)保障度分配模型并采用動態(tài)規(guī)劃的優(yōu)化方法得到備件的配置方案;文獻(xiàn)[4]利用邊際效益方法和仿真優(yōu)化方法解決備件配置優(yōu)化問題的思路;文獻(xiàn)[5]針對Weibull分布型多部件備件組成的串聯(lián)系統(tǒng),提出了備件需求量的改進(jìn)算法。盡管上述方法成功地解決了部分問題,但存在一定的適用性。比如大部分文獻(xiàn)[6~7]僅研究了指數(shù)壽命分布組成的系統(tǒng)的備件配置問題,對于實(shí)際中存在的不同壽命分布類型組成的系統(tǒng)備件,由于其自身的復(fù)雜性(壽命規(guī)律不同),導(dǎo)致備件配置模型十分復(fù)雜,往往僅能通過仿真手段進(jìn)行求解。特別是對于并聯(lián)系統(tǒng)而言,由于結(jié)構(gòu)的特殊性,導(dǎo)致較少文獻(xiàn)對其進(jìn)行研究。由此可見,針對多壽命分布類型組成的并聯(lián)系統(tǒng),如何建立其備件配置優(yōu)化模型,從而可以準(zhǔn)確、直觀地獲得備件優(yōu)化解析解尚需進(jìn)一步研究。
針對由不同壽命分布類型部件組成的并聯(lián)系統(tǒng)的備件配置問題,首先依據(jù)累積失效相等原則,將不同壽命分布部件等效為指數(shù)型部件;然后在此基礎(chǔ)上建立以系統(tǒng)保障度為約束條件,以保障費(fèi)用為目標(biāo)的數(shù)學(xué)模型,并利用邊際效應(yīng)算法,對不同壽命分布類型的部件的備件配置數(shù)量進(jìn)行了優(yōu)化求解;最后通過具體案例的應(yīng)用,對計算方法進(jìn)行了進(jìn)一步的說明和檢驗(yàn)。
衡量一種分布的基本參數(shù)包括失效率、平均壽命、可靠度等。已有的文獻(xiàn)研究表明,盡管平均壽命等效方法對于非指數(shù)型備件的數(shù)量確定具有適應(yīng)性,但具有保守性。但如果按照這種方法進(jìn)行備件需求預(yù)測,方便的同時又會造成備件的浪費(fèi),有時甚至是過度浪費(fèi)。通過數(shù)值分析表明了,累積失效相等方法由于具有較高的近似精度而更加適合艦船裝備備件需求數(shù)量的確定[8~9]。
所謂失效率是指部件在工作到給定時間后,單位時間內(nèi)發(fā)生故障的概率。若部件的壽命分布函數(shù)為F(t),其失效率函數(shù)為λ(t),平均壽命為μ,則其在(0,t]時間內(nèi)平均故障累積函數(shù)(Cumulative Hazard Function,CHF)為
平均故障累積函數(shù)反映了備件在時間t內(nèi)的平均故障次數(shù)。當(dāng)該分布為指數(shù)型分布時,有H′其中λ′為常數(shù)。若采用平均故障累積函數(shù)相等原則,將非指數(shù)分布等效為指數(shù)分布時,即要求H(t)=H′(t)=λ′t,得到
稱上述等效原則為累積失效相等原則,簡稱λ-等效原則[10]。通過該方法,可將非指數(shù)型分布等效成失效率為常數(shù)的指數(shù)分布,表明在相等時間內(nèi)平均累積故障次數(shù)相等。
根據(jù)文獻(xiàn)[5]的推導(dǎo)計算,分別有威布爾分布、正態(tài)分布、伽馬分布型備件的λ-等效方法如下:
1)威布爾型備件的λ-等效方法
對于服從參數(shù)為(β,λ)的威布爾型分布備件,當(dāng)形狀參數(shù)β>1時,其平均壽命為(1+1/β),失效率可等效為λ1=λtβ-1。
當(dāng)保障時間T≤μ時:
當(dāng)保障時間T>μ時:
2)正態(tài)型備件的λ-等效方法
對于服從(μ,σ)的正態(tài)型備件,其失效率可等效為
當(dāng)保障時間T≤μ時:
當(dāng)保障時間T>μ時:
3)伽馬型備件的λ-等效方法
對于服從(α,λ)的伽馬型備件,其平均壽命μ=α/λ,失效率可等效為
當(dāng)保障時間T≤μ時:
當(dāng)保障時間T>μ時:
為簡化計算,本文建立的并聯(lián)系統(tǒng)備件需求預(yù)測模型假設(shè)如下:1)所有部件均為不可修件,且初始時各部件為全新件;2)不考慮環(huán)境、氣候、人為因素等對部件壽命的影響;3)忽略換件時間。
假設(shè)并聯(lián)系統(tǒng)由m個部件組成,系統(tǒng)各部件服從失效率為λ′i(t)(i=1,2,…,m)的不同壽命分布,通過λ-等效方法將各部件服從不同壽命分布的失效率λ′i(t)分別等效為服從λi(i=1,2,…,m)的指數(shù)分布。不同的備件配置費(fèi)用不同,記為Si(i=1,2,…,m)。
在實(shí)際使用中,裝備總是優(yōu)先考慮完成任務(wù)的能力,在此基礎(chǔ)上使費(fèi)用消耗最低。因此,本文建模以并聯(lián)系統(tǒng)的總費(fèi)用S為目標(biāo)函數(shù),以保障度P為約束條件,尋求最優(yōu)的備件配置數(shù)量,使并聯(lián)系統(tǒng)在保障度限制條件下達(dá)到最小的費(fèi)用消耗,則并聯(lián)系統(tǒng)的數(shù)學(xué)模型為
其中,ni為各部件對應(yīng)的備件配置數(shù)量,λk為各部件等效為指數(shù)分布后的失效率,α為并聯(lián)系統(tǒng)保障度。
使用邊際效益方法對并聯(lián)系統(tǒng)的備件配置進(jìn)行優(yōu)化,在尋求系統(tǒng)保障度與備件費(fèi)用之間的最佳平衡時,重復(fù)以下兩個步驟直到滿足保障度要求和費(fèi)用限制:計算給定庫存水平下的系統(tǒng)的保障度Ps;使用邊際效益方法確定各部件對應(yīng)備件的購置順序,選擇購置效益最高的備件作為下一步的購置對象,計算備件購置后的系統(tǒng)保障度。具體計算步驟如下:
步驟1 計算給定庫存水平下的系統(tǒng)的保障度Ps。假設(shè)并聯(lián)系統(tǒng)由m個指數(shù)型部件組成,其第i個部件的備件庫存為ni,對于壽命服從參數(shù)為λi的指數(shù)分布的部件,其系統(tǒng)保障度由式(9)求得。
步驟2 使用邊際效益分析方法。經(jīng)過步驟1,可以算出庫存水平為(n1,n2,…,nm)時的系統(tǒng)保障度Ps,若此時的Ps未達(dá)到目標(biāo)保障度P,就在原庫存水平下選擇某種部件增加一個備件,則有m可能,分別為(n1+1,n2,…,nm),(n1,n2+1,…,nm),…,(n1,n2,…,nm+1)。定義增長因子Ai為第i個部件庫存量從ni增加到ni+1時,系統(tǒng)的新舊保障度與該部件單個費(fèi)用的比值,其計算表達(dá)式如下
其中,Pi+1表示第i個部件增加一個備件后的系統(tǒng)保障度,Si表示第i個部件的單個費(fèi)用。根據(jù)式(10)計算出所有部件對應(yīng)的增長因子Ai,并比較出值最大的Ai,則值最大的Ai就是要增加庫存的部件,其庫存量增加1。把增加后的庫存作為新的庫存水平,重復(fù)以上步驟,直到系統(tǒng)保障度達(dá)到目標(biāo)保障度。
步驟3 得出備件配置數(shù)量。在步驟2中,系統(tǒng)保障度達(dá)到目標(biāo)保障度時庫存水平即為備件配置數(shù)量。
假設(shè)并聯(lián)系統(tǒng)分別由如表1所示的三個部件組成,其壽命分布特征參數(shù)、單價在表中列出,初始庫存為0,要求該系統(tǒng)在工作時間1000小時內(nèi),系統(tǒng)保障度達(dá)到0.97。采用近似方法和邊際效益優(yōu)化后得到的備件配置數(shù)量如表1所示,總費(fèi)用為2.22,系統(tǒng)保障度為0.9972。
表1 并聯(lián)系統(tǒng)各部件相關(guān)數(shù)據(jù)及配置數(shù)量
為驗(yàn)證該配置數(shù)量是否有效,利用Monte-Carlo方法進(jìn)行仿真驗(yàn)證,仿真流程圖如圖1所示。通過10000次的仿真驗(yàn)算,在此配置方案下,該并聯(lián)系統(tǒng)的系統(tǒng)保障度為0.9868,與近似計算得出的系統(tǒng)保障度比較接近,且達(dá)到系統(tǒng)保障度要求,說明近似計算方法是有效的。
本文介紹了將非指數(shù)型部件等效為指數(shù)型部件的λ-等效方法,建立了并聯(lián)系統(tǒng)以保障度為約束條件,以保障費(fèi)用為目標(biāo)的數(shù)學(xué)模型;針對由不同壽命分布類型不可修部件組成并聯(lián)系統(tǒng)的備件配置問題,通過運(yùn)用λ-等效原則,結(jié)合邊際效應(yīng)方法近似計算出系統(tǒng)各部件的備件配置數(shù)量,并通過Monte-Carlo方法進(jìn)行仿真驗(yàn)證,證明該方法是有效的,為并聯(lián)系統(tǒng)的配置優(yōu)化問題提供了一種新思路。
[1]龍軍,康銳,馬麟,等.任意壽命分布的多部件系統(tǒng)備件配置優(yōu)化算法[J].北京航空航天大學(xué)學(xué)報,2007,33(6):698-700.
[2]Willem V J,Rommert D.Spare parts stock control for redundant systems using reliability centered maintenance data[J].Reliability Engineering &System Safety,2011,96(11):1576-1586.
[3]丁定浩.備件保障性的綜合與分配模型[J].電子產(chǎn)品可靠性與環(huán)境試驗(yàn),2006,24(2):1-6.
[4]劉喜春,鄭華,等.備件配置優(yōu)化問題研究[J].系統(tǒng)工程與電子技術(shù),2008,30(10):1934-1937.
[5]劉天華,張志華,等.weibull型多部件備件需求量的改進(jìn)算法[J].海軍工程大學(xué)學(xué)報,2012,24(4):41-45.
[6]陶小創(chuàng),郭霖瀚,肖波平,等.基于備件保障概率分配的備件需求量預(yù)測模型[J].兵工學(xué)報,2012,33(8):975-979.
[7]Karin S S D,Nicole P E,Ana I B,et al.A spare parts model with cold-standby redundancy on system level[J].Computers & Operations Research,2011,38(7):985-991.
[8]劉天華,張志華,梁勝杰,等.一種Weibull型備件需求量的改進(jìn)算法[J].系統(tǒng)工程理論與實(shí)踐,2012,32(5):1124-1128.
[9]馮玉姣.串聯(lián)系統(tǒng)的備件配置方法[J].電子產(chǎn)品可靠性與環(huán)境試驗(yàn),2012,30(5):28-33.
[10]劉天華.艦船裝備備件需求確定與配置優(yōu)化方法研究[D].武漢:海軍工程大學(xué),2012.