楊力，馬小兵

(北京航空航天大學(xué)可靠性與系統(tǒng)工程學(xué)院，北京100191)

0 引言

貯存可用性是指產(chǎn)品在貯存狀態(tài)下，當(dāng)任務(wù)需要時，產(chǎn)品處于可使用狀態(tài)的程度。貯存可用度作為貯存可用性的概率度量，是一種特殊的工作狀態(tài)可用度。對于大部分貯存類產(chǎn)品而言，其大修期或報廢期一般較長。在該期間內(nèi)，由于貯存環(huán)境和維修管理等因素的影響，產(chǎn)品的貯存可用度會逐漸降低，直至一個不可接受的程度。目前，工程上對長期貯存產(chǎn)品一般采用定期檢測維修和定期更換這兩種策略，來保持和恢復(fù)其貯存可用度。

需要注意的是，檢測和更換間隔時間過長，會使產(chǎn)品長時間處于故障狀態(tài)，不可用導(dǎo)致的損失較多，而間隔過短則會導(dǎo)致過度維修，造成不必要的資金浪費。因此，有必要對基于可用度的貯存系統(tǒng)費用問題進(jìn)行研究，尋求最佳維修保障策略。

在貯存可靠性研究方面，文獻(xiàn)［1 -2］基于Bootstrap 對貯存系統(tǒng)的可靠度進(jìn)行了置信限評估。文獻(xiàn)［3 -4］構(gòu)建了由3 類部件組成的系統(tǒng)的貯存可靠度模型，旨在尋找實現(xiàn)使平均可靠度最大的檢修策略。

在工作系統(tǒng)維護(hù)費用優(yōu)化方面，文獻(xiàn)［5 -7］建立了包含不同壽命分布組件的系統(tǒng)混合檢測模型，并用仿真方法確定了長期運行的可用度和費用率。文獻(xiàn)［8 -9］提出了基于連續(xù)檢測的兩部件退化系統(tǒng)視情維修策略，并對其整體費用率進(jìn)行了優(yōu)化。文獻(xiàn)［10 -11］分別針對單組件和多組件系統(tǒng)建立了基于風(fēng)險比例模型的視情維修策略，通過費用率確定最優(yōu)預(yù)防性維修閾值。文獻(xiàn)［12 -15］對于定期檢測與機(jī)會檢測、隨機(jī)檢測結(jié)合的維修更換模型的費用進(jìn)行了對比分析。

然而，目前對于多部件貯存系統(tǒng)的可用度和維修費用缺少研究。該系統(tǒng)的維修保障特點與工作系統(tǒng)差異在于:第一，貯存系統(tǒng)的故障只有通過定期檢測發(fā)現(xiàn)并予以維修，或通過定期預(yù)防性更換消除;第二，貯存系統(tǒng)一般是由幾類維護(hù)策略和維護(hù)周期均不同的部件構(gòu)成。為此，本文開展了對由檢修類和更換類部件構(gòu)成的貯存系統(tǒng)的可用度建模和費用評估優(yōu)化工作。該系統(tǒng)中，第一類部件需進(jìn)行定期更換，更換耗時固定;第二類部件需進(jìn)行定期檢測，檢測有一定的概率無法發(fā)現(xiàn)故障，發(fā)現(xiàn)故障則立即維修如新，維修耗時固定。根據(jù)貯存期間的可靠性變化規(guī)律和具體維護(hù)情況，分別構(gòu)建了兩類部件的瞬時貯存可用度模型，并推導(dǎo)了系統(tǒng)的平均貯存可用度。在此基礎(chǔ)上，建立了該系統(tǒng)在固定貯存期限內(nèi)的維護(hù)費用模型，并基于更新理論對檢測周期以及更換周期進(jìn)行了綜合的優(yōu)化。最后，通過具體的數(shù)值和仿真算例對優(yōu)化結(jié)果進(jìn)行分析和驗證。

1 貯存可用度模型

1.1 模型假設(shè)

1)貯存系統(tǒng)由定期檢測維修和定期更換兩類相互獨立的部件串聯(lián)構(gòu)成。出于簡化模型的考慮，假設(shè)每一類均為單個部件。

2)兩類部件從同一時刻開始貯存，設(shè)為零時刻。

3)更換類部件的更換周期為檢修類部件檢測周期的整數(shù)倍。

4)對檢修類部件的檢測是不完善的，有一定的概率發(fā)生漏檢，而對該類部件的故障維修是完善的，能將部件恢復(fù)到全新的狀態(tài)。

5)檢修類部件的檢測時間忽略不計，而維修時間不可忽略，更換類部件的更換時間不可忽略。

6)更換和故障維修期間系統(tǒng)均為不可用狀態(tài)。

1.2 更換類部件可用度模型

第一類部件通過定期更換策略來保證其高可用度。無論是否發(fā)生故障，該類部件每隔規(guī)定時間均需被全新的部件替換。實際貯存中，更換新部件消耗資源較多，且費用較高，不宜頻繁進(jìn)行，故更換周期一般應(yīng)大于檢修類部件的檢測周期。從方便系統(tǒng)整體維護(hù)的角度出發(fā)，可將該類部件的更換周期定為檢修類部件檢測周期的整數(shù)倍，以避免由于系統(tǒng)處于不可用狀態(tài)而造成的損失。更換-檢測周期比不為整數(shù)以及兩周期相同的情況將在后續(xù)研究中討論，本文暫不考慮。更換策略如圖1所示。

圖1 定期更換類貯存部件的更換策略示意圖Fig.1 Schematic diagram of replacement strategy of replaceable components

設(shè)T'為更換周期，T 為檢測周期，則有T' =NT(N≥2). 并設(shè)更換耗時為Tr，t 時刻的故障累積分布函數(shù)和可靠度函數(shù)分別為Fr(t)和Rr(t)。該類部件在任意時刻t 的瞬時貯存可用度可表示為

式中:m 為正整數(shù)。

具體推導(dǎo)過程如下:

1)當(dāng)0≤t ＜NT 時，該類部件還未進(jìn)行過任何更換，故瞬時貯存可用度為Ar(t)=P(Xr(t)=1)=Rr(t)，其中Xr(t)為該類部件的貯存狀態(tài)隨機(jī)變量

2)當(dāng)mNT≤t ＜mNT +Tr時，該類部件正在更換過程中，即處在不可用狀態(tài)，故瞬時貯存可用度為Ar(t)=P(Xr(t)=1)=0.

3)當(dāng)mNT+Tr≤t ＜(m+1)NT 時，部件處于貯存狀態(tài)。注意到部件每一次更換都相當(dāng)于完成一次更新過程，并考慮到t 之前最近的更新時刻為mNT+Tr，故t 時刻的瞬時貯存可用度為Ar(t)=P(Xr(t)=1)=P(Xr(t -mNT -Tr)=1)=Rr(t -mNT-Tr).

1.3 檢修類部件可用度模型

該類部件每隔T 時刻進(jìn)行定期檢測，檢測有p的概率無法檢測到故障。若檢測結(jié)果為正常，則不對部件進(jìn)行任何操作，繼續(xù)貯存;若檢測結(jié)果為故障，則進(jìn)行故障維修，耗時Tc(0 ＜Tc＜T)且修復(fù)如新，并繼續(xù)貯存。該類部件的檢測和維修策略如圖2所示。

圖2 檢修類貯存部件的檢測維修策略示意圖Fig.2 Schematic diagram of inspection and repair strategy of repairable components

設(shè)t 時刻的故障累積分布函數(shù)和可靠度函數(shù)分別為Fc(t)和Rc(t). 則該類部件在任意時刻t 的瞬時貯存可用度Ac(t)可表示為

式中:k 為正整數(shù);Si=Rc(iT)-Rc((i+1)T).

設(shè)T*f為檢修類部件的首次故障發(fā)生時間，Tf為其首次維修開始時間。根據(jù)貯存時刻的不同，其瞬時貯存可用度可分為以下兩種情況:

1)當(dāng)0 ＜t≤T +Tc時，檢修類部件的瞬時貯存可用度可以表示為

式中:Xc(t)為該類部件的貯存狀態(tài)隨機(jī)變量，Xc(t)=)表示該類部件貯存到t 時刻未發(fā)生過故障的概率，故有P(Xc(t)= 1，t ＜T*

f)= Rc(t);P(Xc(t)=1，T*f≤t≤Tf+Tc)表示該類部件從故障發(fā)生到故障維修完成這段時間內(nèi)處于可用狀態(tài)的概率，顯然有P(Xc(t)=1，T*f≤t≤Tf+Tc)=0. 因此，當(dāng)0 ＜t≤T+Tc時，該類部件的瞬時貯存可用度為Ac(t)=Rc(t)+0 =Rc(t).

2)當(dāng)t ＞T+Tc時，該類部件的瞬時貯存可用度可表示為

結(jié)合情況1 中的推導(dǎo)結(jié)果，可化簡為Ac(t)=Rc(t)+P(Xc(t)=1，t ＞Tf+Tc)，式中:P(Xc(t)=1，t ＞Tf+Tc)表示部件在t 時刻之前已完成過故障的維修，且在t 時刻部件處于正常貯存狀態(tài)的概率。

設(shè)kT+Tc＜t≤(k+1)T+Tc，k =1，2，…. 若部件在t 時刻可用，則部件在t 之前的首次故障可能發(fā)生在［iT，(i+1)T］(i=0，1，…，k-1)內(nèi)。由于每一次定檢均有p 的概率無法發(fā)現(xiàn)故障，因此該故障可能在jT(j=i +1，i +2，…，k)時刻被檢測到，繼而在jT+Tc時刻完成首次故障維修。從零時刻開始至首次故障維修結(jié)束即為一次更新過程。因此:

故當(dāng)kT+Tc＜t≤(k+1)T+Tc時，部件的瞬時貯存可用度為

1.4 貯存系統(tǒng)可用度模型

貯存系統(tǒng)由1.2 節(jié)和1.3 節(jié)中的定期更換和定期檢修兩類部件串聯(lián)構(gòu)成，且兩類部件相互獨立，故障間無影響。因此，某一時刻任意一類部件處于不可用狀態(tài)時，貯存系統(tǒng)也將處于不可用狀態(tài)。綜上，該系統(tǒng)在t 時刻的瞬時貯存可用度為

假定更換耗時Tr大于維修耗時Tc，且N≥2，即不考慮更換周期等于檢測周期的情況。此時系統(tǒng)的瞬時貯存可用度可表示為

式中:m 為正整數(shù);nt=?t/T」;Si=Rc(iT)-Rc((i +1)T).

該貯存系統(tǒng)在時間區(qū)間［0，L］內(nèi)的平均貯存可用度可表示為

式中:L 為固定貯存期限。

2 貯存系統(tǒng)費用模型

實際工程中，大部分貯存產(chǎn)品具有固定的貯存期限，超過該期限一般對產(chǎn)品予以報廢處理。本節(jié)將研究在固定期限內(nèi)系統(tǒng)的貯存費用情況。關(guān)于無貯存期限的系統(tǒng)費用建模將在后續(xù)研究中進(jìn)行，本文暫不考慮。

系統(tǒng)貯存期內(nèi)產(chǎn)生的費用主要由以下四部分組成:更換類部件的更換費用、檢修類部件的檢測費用、檢修類部件的維修費用、貯存系統(tǒng)不可用導(dǎo)致的損失。費用參數(shù)設(shè)置如下:更換類部件單次更換費用為Cr;檢修類部件單次檢測費用為CI;單次維修費用為Cc，系統(tǒng)單位時間不可用導(dǎo)致的損失為Cd.

更換類部件的更換次數(shù)mr為

檢修類部件的檢測次數(shù)mI為

檢修類部件的期望維修次數(shù)E(Nc)為(檢修類部件在iT 時刻不可用)=

貯存系統(tǒng)在［0，L］上期望不可用時間E(Td)為

本模型的優(yōu)化目標(biāo)唯一，即為系統(tǒng)固定貯存期限內(nèi)的平均貯存費用該優(yōu)化目標(biāo)與2 個決策變量相關(guān):1)檢測周期T;2)更換周期與檢測周期長度比N. 由于更換周期NT 不能超過固定貯存期L，且N≥2，故檢測周期T 不能超過L/2. 此外，維修和更換時間一般遠(yuǎn)小于檢測周期，故T ＞max (Tr，Tc).

將(6)式～(9)式代入(10)式并化簡，可得到該費用優(yōu)化問題的數(shù)學(xué)描述:

式中:Ac(iT)為關(guān)于T 的可用度函數(shù)，可通過(2)式求得(0，L)為關(guān)于(N，T)的可用度函數(shù)，可通過(4)式、(5)式求得?？梢园l(fā)現(xiàn)，該貯存系統(tǒng)費用率與貯存部件及貯存系統(tǒng)的可用度均密切相關(guān)。

由于2 個優(yōu)化變量取值均為整數(shù)且數(shù)據(jù)量相對較小，因此擬采用枚舉法，從滿足約束條件(11)式的若干組決策變量(T，N)中搜索使(T，N)達(dá)到最小的最優(yōu)解(T*，N*)。

3 算例

3.1 貯存可用度分析

設(shè)系統(tǒng)的固定貯存期限L 為180 個月。更換類部件故障時間服從威布爾分布，其尺度參數(shù)ηr為550 個月，形狀參數(shù)mr= 1.1，部件更換耗時為0.2 個月;檢修類部件的故障時間服從威布爾分布，其尺度參數(shù)ηc為120 個月，形狀參數(shù)mI=1.7，每次檢測未能檢出故障的概率為0.05，故障維修耗時為0.1 個月。

更換類部件的更換周期T'為24 個月時其瞬時貯存可用度變化規(guī)律如圖3所示。

圖3 更換類部件瞬時貯存可用度變化曲線(更換周期為24 個月)Fig.3 The curres of transient storage availability of replaceable components (replacement cycle:24 moaths)

由圖3可見，首個更換周期內(nèi)貯存可用度不存在零點。從第2 個周期開始，在初始的更換階段可用度為0. 而更換完成時變?yōu)?，且各周期內(nèi)變化規(guī)律相同。該部件的貯存可用度一直維持在某一較高的水平，表明定期更換對部件可用度的保持和恢復(fù)作用較顯著。

當(dāng)檢測周期T 為12 個月時，檢修類部件的瞬時貯存可用度變化規(guī)律如圖4所示。

圖4 檢修類部件瞬時貯存可用度變化曲線(檢測周期為12 個月)Fig.4 The curves of transient storage availability of repairable components (inspestion cycle:24 months)

由圖4可見，瞬時貯存可用度周期性上下振蕩。隨時間增大，振蕩區(qū)間內(nèi)的最低可用度將收斂于某一固定值。此外，由于受維修期間不可用時間的影響，零時刻之后的可用度峰值無法達(dá)到1.

當(dāng)定檢周期為12 個月，更換周期為24 個月時，該系統(tǒng)平均貯存可用度變化規(guī)律如圖5所示。

圖5 系統(tǒng)平均貯存可用度變化規(guī)律(檢測和更換周期分別為12 個月、24 個月)Fig.5 The change rule of average availability of system (inspetion cycle:12 months;raplacement cycle:24 months)

從圖5可看到，系統(tǒng)平均貯存可用度整體上逐漸下降并趨近于某一定值。通過計算表明，該維修-更換串聯(lián)系統(tǒng)存在穩(wěn)態(tài)平均貯存可用度。此外，平均貯存可用度在局部范圍內(nèi)呈周期性的先增后減的變化規(guī)律，但這種變化逐漸趨于不明顯。

3.2 系統(tǒng)貯存費用研究

貯存期間的各項費用如下:更換類部件單次更換費用Cr=100 000 元，檢修類部件單次檢測費用CI=5 000 元，檢修類部件單次維修費用Cc=40 000 元，系統(tǒng)單位時間不可用導(dǎo)致?lián)p失Cd=180 000 元/月。

計算結(jié)果表明，當(dāng)T 為4 個月且N =6 時，即更換類部件更換周期為24 個月，檢修類部件檢測周期為4 個月時，系統(tǒng)在固定貯存期限180 個月內(nèi)的平均費用率最小，為19 026 元/月。平均費用率與兩決策變量的關(guān)系如圖6所示。

圖6 平均費用率與檢測周期及周期比的關(guān)系Fig.6 Average cost rate vs. inspection cycle and cycle radio

當(dāng)更換-檢測周期比N 分別為4、6、8 時，平均費用率與檢測周期的關(guān)系如圖7所示。

圖7 平均費用率與檢測周期的關(guān)系(N 分別為4、6、8)Fig.7 Average cost rate vs. inspection cycle (N=4，6，8)

當(dāng)檢測周期分別為2 個月、4 個月、6 個月時，期望費用率與更換-檢測周期比N 的關(guān)系如圖8所示。

由圖7可見，當(dāng)N 確定時，系統(tǒng)平均費用率隨檢測周期增大先單調(diào)減后單調(diào)增，存在唯一最小值。類似的，由圖8可見，當(dāng)T 確定時，隨著倍數(shù)N 的增大，系統(tǒng)期望費用率先減后增，存在唯一最小值。圖6～圖8綜合分析表明:若檢測和更換周期過短，部件維護(hù)過于頻繁，相應(yīng)費用較多，導(dǎo)致費用率偏高;而若周期過長，則系統(tǒng)不可用時間較多，不可用損失較大，同樣不利于降低費用率。因此，對檢測周期和更換周期進(jìn)行合理的規(guī)劃配比是節(jié)省系統(tǒng)貯存費用的重要途徑。

圖8 平均費用率與周期比N 的關(guān)系(T 為2、4、6 個月)Fig.8 Averge cost rate vs. replacement-inspestion cycle radio N(T=2，4，6 monts)

3.3 費用仿真驗證

費用蒙特卡洛仿真的各項故障分布、費用及維修相關(guān)參數(shù)值見3.1 節(jié)。需要生成的隨機(jī)變量包括:1)服從威布爾分布的更換類部件故障時間f1、f2;2)服從二項分布的檢修類部件故障檢測結(jié)果P(P 取1 或0，1 表示檢測結(jié)果為正常，0 表示檢測結(jié)果為故障)。仿真運行次數(shù)K=10 000. 仿真流程如圖9所示。

部分仿真和計算結(jié)果如表1所示。通過對比可以看到，貯存費用率的計算結(jié)果與仿真結(jié)果基本吻合，從而驗證了模型的準(zhǔn)確性。

表1 貯存費用的計算值和仿真值Tab.1 The calculated and simulated values of storage cost

4 結(jié)論

圖9 貯存費用仿真流程圖Fig.9 Flow chart of storage cost simulation

本文基于實際中長期貯存產(chǎn)品的特點，分別對定期檢修類和定期更換類貯存部件的維修保障策略及特點進(jìn)行了研究，并以此為基礎(chǔ)構(gòu)建了貯存系統(tǒng)可用度模型及貯存費用模型。通過對檢測周期以及更換周期的優(yōu)化選取，實現(xiàn)固定貯存期限內(nèi)平均費用率最低的目標(biāo)。計算結(jié)果表明，最優(yōu)系統(tǒng)維護(hù)策略存在且唯一，且不同的檢測-更換周期組合對費用率有較為明顯的影響。因此對系統(tǒng)各類周期的合理計劃安排是提升貯存經(jīng)濟(jì)性的重要途徑。最后，通過蒙特卡洛仿真驗證了該優(yōu)化計算結(jié)果的準(zhǔn)確性。

References)

［1］羅巍，張春華，譚源源，等. 基于Bootstrap 的可修系統(tǒng)貯存可用度近似置信下限評估方法［J］. 兵工學(xué)報，2010，31(3):391 -395.LUO Wei，ZHANG Chun-hua，TAN Yuan-yuan，et al. Estimating method of approximate confidence lower limits of repairable system storage availability based on bootstrap［J］. Acta Armamentarii，2010，31(3):391 -395.(in Chinese)

［2］羅巍，張春華，譚源源，等. 系統(tǒng)貯存可靠度近似置信下限的Boots trap 評估方法［J］. 宇航學(xué)報，2009，30(4):1725 -1730.LUO Wei，ZHANG Chun-hua，TAN Yuan-yuan，et al. Bootstrap estimate method of approximate confidence lower limits of system storage reliability［J］. Journal of Astronautics，2009，30 (4):1725 - 1730.(in Chinese)

［3］Ito K，Nakagawa T，Nishi K. Extended optimal inspection policies for a system in storage［J］. Mathematical and Computer Modelling，1995，22(10):83 -87.

［4］Ito K，Nakagawa T. Optimal inspection policies for a storage system with degradation at periodic tests［J］. Mathematical and Computer Modelling，2000，31(10):191 -195.

［5］Scarf P A，Cavalcante C A V. Hybrid block replacement and inspection policies for a multi-component system with heterogeneous component lives［J］. European Journal of Operational Research，2010，206(2):384 -394.

［6］Cavalcante C A V，Scarf P A，De Almeida A T. A study of a twophase inspection policy for a preparedness system with a defective state and heterogeneous lifetime［J］. Reliability Engineering ＆System Safety，2011，96(6):627 -635.

［7］Laggounea R，Chateauneuf A，Aissani D. Preventive maintenance scheduling for a multi-component system with non-negligible replacement time［J］. International Journal of Systems Science，2010，41(7):747 -761.

［8］Liu X，Li J，Al-Khalifa K N，et al. Condition-based maintenance for continuously monitored degrading systems with multiple failure modes［J］. IIE Transactions，2013，45(4):422 -435.

［9］Castanier B，Grall A，Bérenguer C. A condition-based maintenance policy with non-periodic inspections for a two-unit series system［J］. Reliability Engineering ＆ System Safety，2005，87(1):109 -120.

［10］Golmakani H R. Condition-based inspection scheme for condition-based maintenance［J］. International Journal of Production Research，2012，50(14):3920 -3935.

［11］Tian Z，Liao H. Condition based maintenance optimization for multicomponent systems using proportional hazards model［J］. Reliability Engineering ＆ System Safety，2011，96(5):581-589.

［12］Laggoune R，Chateauneuf A，Aissani D. Impact of few failure data on the opportunistic replacement policy for multi-component systems［J］. Reliability Engineering ＆ System Safety，2010，95(2):108 -119.

［13］Taghipour S，Banjevic D. Optimum inspection interval for a system under periodic and opportunistic inspections［J］. IIE Transactions，2012，44(11):932 -948.

［14］Taghipour S，Banjevic D. Optimal inspection of a complex system subject to periodic and opportunistic inspections and preventive replacements［J］. European Journal of Operational Research，2012，220(3):649 -660.

［15］Nakagawa T，Mizutani S，Chen M. A summary of periodic and random inspection policies［J］. Reliability Engineering and System Safety. 2010，95(8):906 -911.