金家善,蔡芝明,陳硯橋

(1. 海軍工程大學(xué) 艦船動力工程軍隊重點實驗室， 湖北 武漢　430033；2. 海軍工程大學(xué) 動力工程學(xué)院， 湖北 武漢　430033)

?

基于CML模型的編隊器材攜帶方案優(yōu)化方法*

金家善1,2,蔡芝明1,2,陳硯橋1,2

(1. 海軍工程大學(xué) 艦船動力工程軍隊重點實驗室， 湖北 武漢430033；2. 海軍工程大學(xué) 動力工程學(xué)院， 湖北 武漢430033)

摘要：隨船備件配置優(yōu)化問題研究目前主要集中在只有定量或定性約束上，而對既有定性約束又有定量約束的混合約束問題未見文章提及，同時研究對象多為單船，而對艦艇編隊研究比較少。針對此問題，以艦艇編隊出海任務(wù)準(zhǔn)備階段為背景，以保障費用、艦艇倉庫空間、艦艇最大排水量及艦員維修能力作為模型的約束條件，以編隊備件保障概率為目標(biāo)函數(shù)，采用正態(tài)逆向云模型、邊際效應(yīng)法及拉格朗日乘子法(稱為CML模型)對此類問題進(jìn)行求解，給出系統(tǒng)資源因子求解步驟及分析方法，并對本模型進(jìn)行動態(tài)調(diào)整和優(yōu)化。實例分析表明本模型求解方法和步驟可以為多個定量和定性約束下艦船備件攜帶問題提供新的參考。

關(guān)鍵詞：混合約束；編隊；備件配置；云模型；邊際效應(yīng)法；拉格朗日乘子法

艦艇編隊是現(xiàn)在及未來海軍的主流作戰(zhàn)樣式，其備件攜帶方案將直接影響任務(wù)成功率及戰(zhàn)備完好性，在其離碼頭去執(zhí)行上級賦予的使命任務(wù)之前，必須對其備件攜帶方案進(jìn)行優(yōu)化，達(dá)到在滿足多個定量或定性約束下保障方案最優(yōu)及在滿足任務(wù)要求的前提下最大限度地提高編隊系統(tǒng)可用度，同時達(dá)到使后勤保障經(jīng)費降低的目的。

文獻(xiàn)[1-2]依據(jù)備件領(lǐng)域經(jīng)典的METRIC理論[3]建立了艦船裝備多層級模型；劉勇等[4]以可靠性、質(zhì)量及經(jīng)費為約束條件，建立了艦船遠(yuǎn)航備件存儲決策模型；魏曙寰等[5]依據(jù)艦艇海上維修器材保障模式，建立了備件攜帶方案在保障經(jīng)費及倉庫空間約束下的優(yōu)化模型；費廣玉等[6]對隨船備件影響因素及指標(biāo)進(jìn)行了詳細(xì)說明，選擇滿足率作為評價指標(biāo)；文獻(xiàn)[7-8]主要針對航空類裝備在經(jīng)費和重量約束條件下的可用度進(jìn)行了研究，并建立了相應(yīng)的模型；阮旻智等[9-10]以經(jīng)費、備件滿足率等為約束條件，對艦載裝備維修等級和備件庫存模型進(jìn)行了優(yōu)化；王乃超等[11]建立了備件保障概率在兩個約束因素下的模型，并對庫存指標(biāo)模型進(jìn)行了優(yōu)化；文獻(xiàn)[12-15]以維修保障能力、虛警、運輸時間及保障任務(wù)等作為定性約束條件，給出了相應(yīng)約束條件下的備件最優(yōu)攜帶方案和求解步驟。綜合分析上述研究成果可以看出：首先，艦船攜帶備件方案在多個約束作用下，建立相應(yīng)模型是有一定難度的，故其是研究的熱點，同時也必然是難點；其次，目前研究主要集中在只有定性或定量約束問題上，而對既有定量約束又有定性約束的混合約束問題未見文章提及，同時可以看出定量約束問題的解決方法不能應(yīng)用于定性約束，反之亦不成立，由此可見，尋找一種能夠解決混合約束問題的方法的重要性；再次，單艦隨船備件攜帶方案是之前國內(nèi)外專家研究的主要對象，而把編隊作為目標(biāo)來研究的比較少，未來海軍主要作戰(zhàn)樣式是以艦艇編隊為單位進(jìn)行的，其備件配置方案將直接影響任務(wù)成功率及戰(zhàn)備完好性；最后，從技術(shù)和科學(xué)這方面來說，之前文獻(xiàn)建立了許多編隊備件攜帶優(yōu)化模型，但因邊際效應(yīng)法中的邊際成本只能為一個定量變量，所以這些模型大部分只針對變量為一個的情況，當(dāng)變量為多個時，現(xiàn)有技術(shù)、模型和方法無法解決此類工程問題。因此，提出用正太逆向云模型將定性約束條件轉(zhuǎn)換為定量約束條件，通過引入約束因素因子，再用拉格朗日乘子法將4個約束因素轉(zhuǎn)換成1個拉格朗日乘子，最后用邊際效應(yīng)法求出最優(yōu)方案。

艦艇編隊備件攜帶方案需考慮多個定性及定量約束綜合權(quán)衡得到，常用到的將定性約束條件轉(zhuǎn)換成定量約束條件的方法有：專家打分法[16]、量化加權(quán)[17]、層次分析[18]等，這些方法利用不同模型對不確定性進(jìn)行了研究，優(yōu)點是：目標(biāo)明確，邊界條件清晰，在研究時分別從模糊性和不確定性兩個方面建立模型求解，但二者往往是不能完全分開的。因而，李德毅等建立了云模型[19]，通過相互映射的方法將二者聯(lián)系到了一起，同時該分析方法已被廣泛應(yīng)用到天氣預(yù)測[20]、信息技術(shù)[21]等領(lǐng)域中。

艦船出海不能無限制地攜帶維修需要的所有備件，限制因素主要有保障經(jīng)費、排水量、艦員維修能力和倉庫空間，這些因素對艦艇備件攜帶都具有重要的影響，在建立模型求解的過程中，不能顧此失彼，必須同時考慮，已有模型和方法針對此問題研究比較少，本文以備件費用、質(zhì)量、體積及維修工時作為約束條件，以備件保障概率作為目標(biāo)函數(shù)，給出了基于CML模型的艦艇編隊備件攜帶方案的優(yōu)化程序和步驟，最后通過對單個約束條件下的備件攜帶方案和多個定量及定性約束條件下的備件攜帶方案進(jìn)行對比，驗證了模型和計算步驟的合理性。

1模型建立

1.1目標(biāo)函數(shù)

艦船備件攜帶方案和陸地裝備主要以后勤保障經(jīng)費為約束，不同的是，前者還要考慮排水量、倉庫空間大小及艦員維修能力，在滿足多個約束條件下使其最優(yōu)，文獻(xiàn)[6, 22-25]建立了不同目標(biāo)函數(shù)下的最優(yōu)備件攜帶方案，但備件保障概率很少被作為優(yōu)化對象，上述目標(biāo)在艦艇編隊實際出海備件攜帶方案制定過程中，管理部門是難以掌控和把握的，因而，本文以編隊備件保障概率Ps為目標(biāo)函數(shù)。船上電子設(shè)備累計概率密度函數(shù)都服從指數(shù)分布，計算每項備件保障概率Pi是一樣的，即：

(1)

式中：i(i=0,1,2,…,I)為備件編號；j(j=0,1,…,J)為第i項元器件的備件需求量；Ni為第i項電子元件單機安裝數(shù)；λi為第i項元器件故障率；ti為第i項備件供應(yīng)周期。

艦艇裝備主要由兩個層級組成，即艦船現(xiàn)場可更換或拆卸單元(Line-Replaceable Unit， LRU)和車間可更換或拆卸單元(Shop-Replaceable Unit， SRU)，LRUi年平均故障率λi為：

(2)

式中：DCi為占空比；RIPi為故障件可以在故障現(xiàn)場修理的概率；T0為裝備年平均工作時間；L為隨船裝備配置數(shù)量；MTBFi為部件i的平均故障間隔時間；RtOKi為部件i的重測完好率。依據(jù)故障樹分析原理，母體的需求率是由其各個分部件SRUk(k∈sub(i))的需求率所引起的，若裝備LRUi發(fā)生故障是由其組成部件引起的概率為Pik，則SRUk的需求率(故障率)λk為：

λk=λi·Pik

(3)

(4)

把式(4)代入式(3)得：

(5)

第y(y=1,2,…,Y)條船備件保障概率psy為：

psy=pscy·psdy

(6)

式中：pscy為修理所需求的備件放在故障現(xiàn)場的概率；psdy為假設(shè)在事先設(shè)定的約束條件下，需求的備件放在故障現(xiàn)場且可以獲得的概率。其計算方法為：

(7)

(8)

式中，λTy為第y船總故障率，且

(9)

根據(jù)式(6)～(9)，按上述給定的計算步驟和方法可得到第y船的保障概率Psy，即：

(10)

編隊艦艇在海上彼此之間的距離遠(yuǎn)比離岸上保障基地距離要近的多，因而可將編隊看成一個整體系統(tǒng)，艦船修理車間不進(jìn)行串件拼修的前提下，艦艇編隊備件保障概率依據(jù)式(10)并考慮一些其他修理因素綜合權(quán)衡分析可得到，即：

(11)

其中：

(12)

此時，目標(biāo)函數(shù)為:

maxPs

(13)

1.2約束條件

艦船備件攜帶問題可以理解為在多個定量和定性影響因素下，綜合考慮各種維修因素，通過攜帶不同種類和數(shù)量的備件，使艦艇編隊Ps達(dá)到最優(yōu)。

艦艇倉庫攜帶的所有備件體積之和必須小于艦艇倉庫空間所能承受的上限，即:

(14)

式中:viy為備件LRU，SRU的體積；Vt為艦艇編隊倉庫空間所能承受的上限值。

雖然與倉庫體積及排水量相比，保障經(jīng)費并非剛性約束，即以費用最小化和保障概率最大化共同作為優(yōu)化目標(biāo)更為實際，但為簡化問題，本文將其作為約束因素，編隊出海購買所有備件的費用總和必須小于保障經(jīng)費的上限，即：

(15)

式中：ciy為單件LRU,SRU的費用；Ct為決策者所能給予的保障經(jīng)費的最大值。

編隊攜帶備件所有重量之和必須小于艦船載荷所能承受最大值的上限，即：

(16)

式中：miy為單件LRU，SRU的質(zhì)量；Mt為事先給定的倉庫質(zhì)量指標(biāo)。

維修能力對于艦船的重要性是顯而易見的[6,14-15,23]。標(biāo)準(zhǔn)工時[26]：在規(guī)定的工作環(huán)境里，按規(guī)定的流程，以標(biāo)準(zhǔn)工作者正常速度完成裝備修理所必要的時間。編隊艦員更換及調(diào)試裝備所需的維修工時之和必須小于所有維修人員所能承受的最大值，即：

(17)

式中：wiy為第y條船的第i個備件從拆卸開始到測試結(jié)束所需要的維修工時；Wt為編隊艦員所能承受維修工時的最大值；w′y為第y船上單位時間內(nèi)所具有的平均維修保障能力；Ty為任務(wù)周期。

1.3模型的建立

編隊離碼頭執(zhí)行上級下達(dá)的某項具體任務(wù)之前，決策者要綜合影響備件攜帶的所有因素，在滿足各種因素所能承受值的最大值的前提下，求出編隊備件保障概率Ps最優(yōu)值，所建模型如下：

(18)

2模型求解方法

2.1正態(tài)云模型的引入

目前將定性轉(zhuǎn)換成定量的方法有很多，如專家打分法、模糊層次分析法等[16-18]，其包含著一些數(shù)學(xué)方法和計算模型，但上述方法并沒有將研究對象的模糊性和隨機性很好地統(tǒng)一起來[19]，因而李德毅建立了定性條件和定量條件轉(zhuǎn)換模型，并給出了其求解步驟和程序，通過構(gòu)建二者之間的映射關(guān)系，將定性和定量之間的信息進(jìn)行處理，達(dá)到相互轉(zhuǎn)化的目的。

定義1[27]正態(tài)分布通常用均值和方差兩個特征值來表示；模糊理論中用的最多的是鐘形隸屬度函數(shù)，通常用

(19)

來定義，正態(tài)云模型是以兩者為基礎(chǔ)發(fā)展而來的一種新的模型。

定義2[27]設(shè)U是一個由很多精確數(shù)值確定的定量論域，C屬于U的定性概念表現(xiàn)，若確定的數(shù)量值x∈U，且x是定性概念C的一次隨機實現(xiàn)，若x滿足：

x～N(Ex,En′2)

(20)

其中

En′～N(En,He2)

(21)

且x對C的確定度滿足：

(22)

則x在論域U上的分布稱為正態(tài)云。

定義3設(shè)U是維修能力T用精確數(shù)值給定的定量論域，X?U，T是U空間上的關(guān)于維修能力的定性概念，若對于維修工時x(x∈X)都存在一個穩(wěn)定傾向的隨機數(shù)CT(x)∈[0,1]，稱為維修工時x對維修能力T的隸屬度，即：

CT(x):U→[0,1],?x∈X(X?U),x→CT(x)

(23)

無論在自然科學(xué)還是其他學(xué)科(如社會科學(xué))都已經(jīng)論證了李德毅所給出的正態(tài)云模型的科學(xué)性、合理性及普適性[27]。

一個云滴就是云模型中一次定性的隨機實現(xiàn)，云模型中云滴是無限的，但其取值并不是絕對或連續(xù)的無限，文中處理的維修工時就是一次維修能力的定性實現(xiàn)，當(dāng)?shù)玫蕉鄠€維修工時云滴，便得到了裝備的維修能力；常規(guī)意義上的定性概念(如高、中、低等)是有限取值和具體離散的值。從以上兩個方面可以看出，二者相互聯(lián)系但也有客觀上的區(qū)別，這些不同決定了云模型中的定性概念和常規(guī)意義上的定性概念有本質(zhì)區(qū)別。

2.2基于正態(tài)逆向云模型的定性與定量變量轉(zhuǎn)換模型

從文獻(xiàn)[27]給出的正態(tài)云模型可以看出，用該方法來還原云模型中的數(shù)字特征，更加符合艦船備件攜帶工程方面的要求，易于推廣到高維逆向云模型中去，該方法是以統(tǒng)計原理為基礎(chǔ)的，在實際計算過程中，不同的樣本值有可能引起少量的誤差，主要體現(xiàn)在計算超熵He時，容易出現(xiàn)超熵值為虛數(shù)的情況，若出現(xiàn)超熵值為虛數(shù)，則說明本次求解過程失敗，需要選擇新的云模型重新計算，為此，對算法進(jìn)行了改進(jìn)，對差值取絕對值并加上了控制參數(shù)，提高了計算的魯棒性，減小了計算誤差。算法具體計算步驟如下：

步驟1：輸入第i個備件R個云滴樣本的數(shù)字值xir(r=1,2,…,R);

步驟2：依據(jù)R個云滴數(shù)字值xir可求解出均值大小：

(24)

步驟3：計算樣本期望：

(25)

步驟4：計算樣本熵：

(26)

步驟5：計算R個樣本云滴方差：

(27)

步驟6：計算樣本超熵：

(28)

式中，Ci為控制參數(shù)，一般情況下取值為1，可對超熵進(jìn)行適當(dāng)調(diào)節(jié)，使其處于一個管理者可接受的合理范圍。

步驟7：重復(fù)步驟1～6共I次，得到所有備件云模型的特征參數(shù)。

2.3模型的求解方法

由于涉及的約束變量為定性1個和定量3個，而邊際效應(yīng)法求解的原理是用邊際效益除以邊際成本得到邊際效益值，通過逐次比較邊際效益值的大小，從而確定編隊攜帶備件的種類和數(shù)量，但邊際效應(yīng)法要求邊際成本必須為一個定量變量，因而，一方面需要尋找一種可以將定性變量變成定量變量的方法，另一方面需要尋找將4個定量變量變成1個定量變量的方法，最后，用邊際分析法求解得到編隊庫存?zhèn)浼詈侠淼臄y帶方案。

用云模型將定性變量轉(zhuǎn)換為定量變量，具體見2.1和2.2節(jié)，此方法求得的定量值可以作為拉格朗日乘子輸入值；單純拉格朗日乘子法主要思路是引入新的參數(shù) 進(jìn)而求出極值解，本文主要是借用該數(shù)學(xué)方法的思想，引入4個變量因子，用拉格朗日乘子法將4個定量變量轉(zhuǎn)換為1個定量變量(一個備件對應(yīng)一個定量值)，文獻(xiàn)[9,11]已驗證了本文所用方法的正確性；最后，將備件定量值作為邊際成本變量輸入到邊際效應(yīng)方法中，求出編隊最優(yōu)備件攜帶方案，下面是該方法的具體求解步驟。

步驟1：從式(11)和式(18)可以看出，艦艇編隊備件攜帶問題的影響因素有很多，需要將這些因素(艦艇載荷大小、備件體積、經(jīng)費及艦員維修能力)統(tǒng)一轉(zhuǎn)化為zi，即：

zi=λmmi+λvvi+λcci+λwwi

(29)

式中：zi為單元i的系統(tǒng)保障資源約束因子；λm為質(zhì)量因子；λv為體積因子；λc為費用因子；λw為艦員維修能力因子。

步驟2：將備件的3個定量因素和1個定性約束因素，依據(jù)式(29)轉(zhuǎn)化為系統(tǒng)資源約束Zi后，在其基礎(chǔ)之上運用適當(dāng)?shù)臄?shù)學(xué)變換，即可求出編隊Ps，即：

(30)

步驟3：依據(jù)凸函數(shù)疊加之后仍然是凸函數(shù)的原理，依據(jù)式(29)及文中給的程序和步驟重復(fù)計算I次，可得到當(dāng)前系統(tǒng)約束組合為：

zt=(z1,z2,…，zi，…,zI)

(31)

步驟4：對式(31)做一階差分，即：

ΔPzi(si)=Pzi(si+1)-Pzi(si)

(32)

步驟5：將式(31)除以式(32)，即：

對Δ矩陣按邊際效應(yīng)分析法進(jìn)行處理，逐次選擇最大值，如此進(jìn)行多次循環(huán)，當(dāng)所有影響因素不再滿足指標(biāo)上限時停止。

步驟6：循環(huán)步驟4～5，直到攜帶備件種類和數(shù)量不能再增加時停止。

3初始約束因子的確定及動態(tài)調(diào)整

在求解得到備件綜合資源約束Zi之前，必須先求出艦船備件配置初始方案，從而確定艦船載荷因子λm、艦船倉庫空間因子λv、費用因子λc、維修因子λw的初值，記初始艦船載荷因子為λm0，初始倉庫空間因子為λv0，初始維修因子為λw0，模型求解方法和優(yōu)化步驟如下：

步驟1：當(dāng)艦艇倉庫空間、艦艇載荷及艦員維修能力不被作為約束條件時，通過費效曲線計算得到一組初始備件攜帶方案，即：

sc0=(s01j,s02j,…,s0ij,…,s0Ij)

(34)

步驟2：艦艇載荷、維修工時、保障經(jīng)費及艦艇倉庫空間都可在式(34)的基礎(chǔ)上，運用文中模型給的步驟和程序求出，然后用兩兩之間的相互比值λcm0，λcv0，λcw0作為初始因子，即：

(35)

(36)

(37)

步驟3：將式(35)～(37)值代入式(29)，即：

zc0ij=cij+λcv0vij+λcm0mij+λcw0wij

(38)

步驟4：依據(jù)式(33)和式(38)按文中給定步驟再次求解備件攜帶方案，即：

sc1=(sc11j,sc12j,…,sc1ij,…,sc1Ij)

(39)

步驟5：在計算過程中可能會出現(xiàn)一種情況，計算得到總的備件質(zhì)量、倉庫體積和保障費用超過了指標(biāo)最大值上限，需要對初始各種約束條件影響因子進(jìn)行適當(dāng)變換，即：

(40)

(41)

(42)

步驟6：依據(jù)式(40)～(42)對約束條件因子進(jìn)行適當(dāng)處理之后，絕大部分時候都會滿足約束條件指標(biāo)要求，但還有可能出現(xiàn)一種情況：無論怎么調(diào)整備件因子值的大小，求出的備件攜帶方案都不能滿足所有約束指標(biāo)的要求，總會有1個或多個指標(biāo)不滿足要求。此時，需要調(diào)整約束指標(biāo)上限值，可增加艦艇載荷所能承受的指標(biāo)M、艦船倉庫所能承受指標(biāo)V、維修人員能承受最大工作工時指標(biāo)W或者艦船保障經(jīng)費指標(biāo)C。

依據(jù)第2節(jié)和本節(jié)內(nèi)容，可以得出CML模型計算的流程如圖1所示。

圖1　CML模型流程圖Fig.1　Flow chart of CML model

4案例分析

3艘最新型艦船準(zhǔn)備出海去完成上級賦予的使命任務(wù)，離碼頭之前需要綜合權(quán)衡保障經(jīng)費、維修能力、倉庫空間及排水量等約束條件，使保障方案在滿足所有約束條件前提下最優(yōu)。由于艦船裝備種類和數(shù)量繁多不能一一列舉，同時本文重點在于建立編隊備件方案求解模型，并給出模型優(yōu)化步驟，因此，表1只給出了會給艦艇航行帶來致命性影響的一些重要裝備，并對屬性值進(jìn)行了適當(dāng)處理，其他裝備求解

過程和計算流程是一樣的。

按照文獻(xiàn)[27]中關(guān)于維修工時的標(biāo)準(zhǔn)，去該類艦船進(jìn)行調(diào)研，依據(jù)裝備維修工作記錄本等資料獲取這10個備件維修所需工時，具體見表2。

按文中2.2節(jié)給定的計算步驟，得到10個備件的三個特征參數(shù)，并據(jù)此得出10個備件的云分布圖，云滴分布符合決策人員要求的標(biāo)準(zhǔn)，具體見表3。依據(jù)表1～3的數(shù)據(jù)，按2.3節(jié)給定的模型求解方法及第3節(jié)模型優(yōu)化方法，求解得到不同影響因素下對應(yīng)的保障方案，具體見表4。

表1　備件屬性表

表2　備件維修工時表

表3　備件云模型數(shù)字特征

表4　不同約束條件下編隊隨船備件攜帶方案

艦艇編隊給定的保障經(jīng)費最大值為25.2萬元，載重最大值為307 kg，倉庫最大體積為13.5 m3，維修工時指標(biāo)為65.1 h。依據(jù)表4的結(jié)果可以計算推導(dǎo)出表5，從表5的5個可行性方案可以看出，只有方案3和方案5的各項約束值均滿足給定的影響因素所能承受最大值上限的要求，方案1、方案2和方案4分別有指標(biāo)超過了給定指標(biāo)，均不滿足要求。

表5　不同約束下編隊隨船備件攜帶方案結(jié)果

從表5可以看出：

1)方案1，2和4因其都是將一個影響因素作為約束條件，求出的備件最優(yōu)配置方案都不滿足所有約束因素的要求，從上面可以看出隨船備件配置必須考慮所有的影響因素，此時求出的才是符合所有約束因素要求的最優(yōu)方案。

2) 一方面，綜合考慮編隊所有影響因素，按文中模型計算得到的備件攜帶方案，滿足所有約束條件且是最優(yōu)方案；另一方面，只考慮編隊一個影響因素，按邊際效應(yīng)法計算得到的備件攜帶方案，是否能滿足所有約束條件是偶然的，因此，這也從另一角度證明了本文研究的必要性。

3) 已有研究方法主要側(cè)重對單船或陸地裝備進(jìn)行研究，取得了很多成果，但這些方法無法科學(xué)、合理地解決艦艇編隊備件攜帶方案的制定問題，故提出CML模型和方法，通過和已有方法將編隊中的艦艇分開來研究的情況進(jìn)行對比，證明了本文方法在提高裝備戰(zhàn)備完好性的同時，可以有效降低保障經(jīng)費。

不同約束下的艦艇編隊最優(yōu)保障曲線如圖2～6所示。

圖2　費用約束下的最優(yōu)Ps曲線Fig.2　Optimal curve of support probability under cost constraints

圖3　質(zhì)量約束下的最優(yōu)Ps曲線Fig.3　Optimal curve of support probability under mass constraints

圖4　維修能力約束下的最優(yōu)Ps曲線Fig.4　Optimal curve of support probability under maintenance per hour constraints

圖5　體積約束下的最優(yōu)Ps曲線Fig.5　Optimal curve of support probability under volume constraints

圖6　混合約束下的最優(yōu)Ps曲線Fig.6　Optimal curve of support probability under mix-constraints

綜合分析約束因素下的最優(yōu)Ps曲線，可以得到：

1)隨船備件配置方案不僅受備件種類和數(shù)量的影響，同時也受其各種指標(biāo)條件約束。根據(jù)不同因素條件，運用文中給定的模型求解得出不同優(yōu)化曲線；

2)方案制定者依據(jù)最優(yōu)Ps曲線，可以制定編隊出海之前在相應(yīng)影響因素下的最優(yōu)備件攜帶種類和數(shù)量。例如，以圖5中點(9.68，0.817 4)為例進(jìn)行說明，當(dāng)編隊Ps達(dá)到81.74%時，所有備件體積之和為9.68 m3，因此，設(shè)定的體積最小指標(biāo)值必須達(dá)到Vt≥9.68 m3，若給定的體積指標(biāo)小于9.68 m3時，相應(yīng)地應(yīng)該通過適當(dāng)降低Ps來達(dá)到滿足約束條件上限要求；

3)上述影響因素和最優(yōu)Ps在約束條件和目標(biāo)函數(shù)之間是可以相互轉(zhuǎn)換的，可根據(jù)具體工程決策者的要求具體分析。

5結(jié)論

主要針對編隊隨船備件配置方案中既有定量約束又有定性約束的混合問題進(jìn)行了研究。可以為解決多個定量及定性約束的混合問題提供借鑒，當(dāng)約束條件增加時，只要對模型系統(tǒng)資源約束因子適當(dāng)調(diào)整即可，計算流程不變，模型的適用性比較強。模型可以有效克服主觀方法的不足，同時兼顧復(fù)雜系統(tǒng)的隨機性和模糊性。建立的模型、求解方法及優(yōu)化程序?qū)ε灤瑐浼y帶有重要參考價值，而且對于其他復(fù)雜系統(tǒng)領(lǐng)域裝備等同樣具有一定的參考意義。一方面，編隊的組成形式主要取決于任務(wù)樣式，任務(wù)不同組成形式不同，但求解過程和用到的模型及方法是相似的；另一方面，若編隊(有2艘或以上艦船即可組成編隊)中艦艇數(shù)量隨著實際情況有變化(假設(shè)有Z艘且Z>3)，艦船數(shù)量的增加將主要是增加計算工作量，而用到的模型和方法都是相同的。

隨船備件配置方案不僅與備件種類和數(shù)量有關(guān)，還與任務(wù)強度及倉庫地理位置等相關(guān)，同時本文并未考慮艦艇編隊內(nèi)部備件的串件拼修，這部分內(nèi)容將是接下來研究的重點。

參考文獻(xiàn)(References)

[1]Sleptchenko A, Van der Heijden M C, Van Harten A. Effects of finite repair capacity in multi-echelon, multi-indenture service part supply systems[J]. International Journal of Production Economics, 2002, 79(3): 209-230.

[2]Rustenburg W D, Van Houtum G J, Zijm W H M. Spare parts management at complex technology-based organizations: an agenda for research[J]. International Journal of Production Economics, 2001, 71(1/2/3): 177-193.

[3]Sherbrooke C C. METRIC: a multi-echelon technique for recoverable item control[J]. Operations Research, 1968, 16(1): 122-141.

[4]劉勇, 蓋強, 趙翀, 等. 多約束下艦船遠(yuǎn)航備件儲量決策[J]. 艦船科學(xué)技術(shù), 2013, 35(11): 144-147.

LIU Yong, GAI Qiang, ZHAO Chong, et al. Research on spare parts storage decision for ships on sailing in multi-restrain[J].Ship Science and Technology, 2013, 35(11): 144-147. (in Chinese)

[5]魏曙寰, 陳硯橋, 金家善. 空間和費用約束條件下的隨船備件配置優(yōu)化方法[J]. 系統(tǒng)工程與電子技術(shù), 2013, 35(12): 2540-2544.

WEI Shuhuan, CHEN Yanqiao, JIN Jiashan. Warship spare parts allotment method under space and cost constraints[J]. Systems Engineering and Electronics, 2013, 35(12): 2540-2544. (in Chinese)

[6]費廣玉, 張志華, 劉軍, 等. 艦船隨艦備件配置方法[J]. 指揮控制與仿真, 2014, 36(3): 133-136.

FEI Guangyu, ZHANG Zhihua, LIU Jun, et al. Configuration method of spares on warship[J]. Command Control & Simulation, 2014, 36(3): 133-136. (in Chinese)

[7]Lee L H, Chew E P, Teng S, et al. Multi-objective simulation based evolutionary algorithm for an aircraft spare parts allocation problem[J]. European Journal of Operational Research, 2008, 189(2): 476-491.

[8]Kline R C, Bachman T C. Estimating spare parts requirements with commonality and redundancy[J]. Journal of Spacecraft and Rockets, 2007, 44(4): 977-984.

[9]阮旻智, 李慶民, 張光宇, 等. 多約束下艦船裝備攜行備件保障方案優(yōu)化方法[J]. 兵工學(xué)報, 2013, 34(9): 1144-1149.

RUAN Minzhi, LI Qingmin, ZHANG Guangyu, et al.Optimization method of carrying spare parts warship equipment under multi-constraints[J]. Acta Armamentarii, 2013, 34(9): 1144-1149. (in Chinese)

[10]阮旻智, 李慶民, 彭英武, 等. 多指標(biāo)約束下艦載裝備維修級別建模與優(yōu)化[J]. 系統(tǒng)工程與電子技術(shù), 2012, 34(5): 955-960.

RUAN Minzhi, LI Qingmin, PENG Yingwu, et al. Modeling and optimization for repair level of shipborne equipment under multi-constraints[J]. Systems Engineering and Electronics, 2012, 34(5): 955-960. (in Chinese)

[11]王乃超, 康銳. 多約束條件下備件庫存優(yōu)化模型及分解算法[J]. 兵工學(xué)報, 2009, 30(2): 247-251.

WANG Naichao， KANG Rui. An optimization model for inventory spares under multi-constraints and its decomposition algorithm[J]. Acta Armamentarii, 2009, 30(2): 247-251.(in Chinese)

[12]卞潔輝, 郭霖瀚, 康銳, 等. 虛警影響下的備件庫存優(yōu)化[J]. 北京航空航天大學(xué)學(xué)報, 2014, 40(3): 1-4.

BIAN Jiehui, GUO Linhan, KANG Rui, et al. Spares optimization under the influence of false alarm[J]. Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(3): 1-4. (in Chinese)

[13]徐立, 李慶民, 阮旻智. 具備有限維修能力的艦船編隊保障方案優(yōu)化 [J]. 系統(tǒng)工程與電子技術(shù), 2014, 36(11): 2226-2232.

XU Li, LI Qingmin, RUAN Minzhi. Support project optimization for warship formation with finite repair capacity[J]. Systems Engineering and Electronics, 2014, 36(11): 2226-2232. (in Chinese)

[14]王正元, 曹繼平, 朱昱, 等. 考慮維修能力的戰(zhàn)時備件資源配置方法研究[J]. 兵工學(xué)報, 2014, 35(5): 719-724.

WANG Zhengyuan,CAO Jiping,ZHU Yu, et al. An optimization model of wartime spare parts resource allocation on consideration of maintainability[J]. Acta Armamentarii, 2014, 35(5): 719-724. (in Chinese)

[15]王正元, 曹繼平, 朱亞紅, 等. 可移動備件倉庫設(shè)置優(yōu)化方法[J]. 系統(tǒng)工程與電子技術(shù), 2013, 35(11): 2352-2355.

WANG Zhengyuan,CAO Jiping, ZHU Yahong, et al. Optimization model on removable spare parts warehouse deployment[J]. Systems Engineering and Electronics, 2013, 35(11): 2352-2355. (in Chinese)

[16]郭海鵬, 黃勝, 王超. 基于改進(jìn)德爾菲法的艦船總體方案群決策方法[J]. 上海交通大學(xué)學(xué)報, 2014, 48(4): 13.

GUO Haipeng, HUANG Sheng, WANG Chao. Group decision-making method of warship overall scheme based on improved Delphi[J]. Journal of Shanghai Jiaotong University, 2014, 48(4): 13. (in Chinese)

[17]蔡坦, 劉衛(wèi)寧, 劉波. 一種新的基于直覺模糊集的制造云服務(wù)優(yōu)選方法[J]. 中國機械工程, 2014, 25(3): 352-356.

CAI Tan，LIU Weining，LIU Bo. A new method of cloud manufacturing service optimal-selection based on intuitionistic fuzzy set[J]. Chinese Journal of Mechanical Engineering, 2014, 25(3): 352-356. (in Chinese)

[18]楊兆軍, 王繼利, 李國發(fā), 等. 沖壓機床可靠性增長的模糊層次分析預(yù)測方法[J]. 吉林大學(xué)學(xué)報: 工學(xué)版, 2014，44(3): 686-691.

YANG Zhaojun, WANG Jili, LI Guofa, et al. Reliability growth prediction based on fuzzy analytical hierarchy process for punching machines[J]. Journal of University (Engineering and Technology Edition), 2014， 44(3): 686-691. (in Chinese)

[19]李德毅, 孟海軍. 隸屬云和隸屬云發(fā)生器[J]. 計算機研究與發(fā)展, 1995, 32(6): 15-20.

LI Deyi, MENG Haijun. Membership cloud and membership cloud generators[J]. Journal of Computer Research and Development, 1995, 32(6): 15-20. (in Chinese)

[20]Farley R D, Nguyen P, Orville H D. Numerical simulation of cloud seeding using a three-dimensional cloud model[J]. The Journal of Weather Modification, 2014, 26(1): 113-124.

[21]Rezaei R, Chiew T K, Lee S P, et al. Interoperability evaluation models: a systematic review[J]. Computers in Industry, 2014, 65(1): 1-23.

[22]Molenaers A, Baets H, Pintelon L, et al. Criticality classification of spare parts: a case study[J]. International Journal of Production Economics, 2012, 140(2): 570-578.

[23]Romeijnders W, Teunter R, Van Jaarsveld W. A two-step method for forecasting spare parts demand using information on component repairs [J]. European Journal of Operational Research, 2012, 220(2): 386-393.

[24]Costantino F, Di Gravio G, Tronci M. Multi-echelon, multi-indenture spare parts inventory control subject to system availability and budget constraints[J]. Reliability Engineering & System Safety, 2013, 119: 95-101.

[25]Caglar D, Li C L, Simchi-levi D. Two-echelon spare parts inventory system subject to a service constraint[J]. IIE Transactions, 2004, 36(7): 655-666.

[26]張峰, 殷秀清. 基于產(chǎn)品批量生產(chǎn)的標(biāo)準(zhǔn)工時測量[J]. 山東理工大學(xué)學(xué)報 (自然科學(xué)版), 2014, 28(3): 56-59.

ZHANG Feng, YIN Xiuqing. The standard time measurement based on product batch production[J]. Journal of Shangdong University of Technology(Natural Science Edition), 2014, 28(3): 56-59. (in Chinese)

[27]李德毅, 劉常昱. 論正態(tài)云模型的普適性[J]. 中國工程科學(xué), 2004, 6(8): 28-34.

LI Deyi, LIU Changyu. Study on the universality of the normal cloud model[J]. Engineering Science, 2004, 6(8): 28-34. (in Chinese)

doi:10.11887/j.cn.201602022

*收稿日期：2015-03-12

基金項目：國家部委基金資助項目(51319060103)；中國博士后科學(xué)基金資助項目(2013T60921)

作者簡介：金家善(1962—)，男，山東青島人，教授，博士，博士生導(dǎo)師，E-mail：jinjiashan401@163.com

中圖分類號：TJ83; U662.1

文獻(xiàn)標(biāo)志碼：A

文章編號：1001-2486(2016)02-128-10

Optimization method of carrying maintenance material support project for warship formation based on CML model

JIN Jiashan1,2, CAI Zhiming1,2, CHEN Yanqiao1,2

(1. Military Key Laboratory for Naval Ship Power Engineering, Naval University of Engineering, Wuhan 430033, China；2. College of Power Engineering, Naval University of Engineering, Wuhan 430033, China)

Abstract：At present, studies in spares allotment optimization area mainly focus on the conditions that there is only qualitative or quantitative constraint. Nevertheless, the mix-conditions have not been reported to the best of our knowledge from the literatures. The researches usually concern the single warship, but fewer show concern about the warship formation. In view of this problem, the mission preparation of warship formation before sailing was taken as a research background; the warship formation support cost, warship storage space, warship carrying ability and crew maintenance ability were adopted as the mixed-constraints conditions; the spare parts support probability was served as the objective function. A normal reverse distribution cloud model, the marginal effect and the Lagrange multiplier methods (CML model) were used to solve such warship formation spare parts allotment problem. A method was proposed for the determination and dynamical updating of system resource constraint factors, meanwhile, the selecting methods and calculating process of the model was optimized. The reliability of the method which can provide a new trial for the problem of warship formation spare parts allotment under the mixed constraints condition was proved by an actual case.

Key words：mixed-conditions; warship formations; spares allotment; cloud model; marginal effect method; Lagrange multiplier method

http://journal.nudt.edu.cn