魯金明，姚玉南，李澤之

(武漢理工大學(xué)能源與動(dòng)力工程學(xué)院，湖北武漢 430063)

維修保養(yǎng)是維持設(shè)備正常使用的重要技術(shù)手段。船舶設(shè)備的維修保養(yǎng)體系經(jīng)歷了3個(gè)發(fā)展階段:事后維修、定期維修 (計(jì)劃預(yù)防維修)、視情維修。其中視情維修時(shí)以狀態(tài)檢測(cè)和預(yù)測(cè)技術(shù)為技術(shù)手段，針對(duì)設(shè)備的運(yùn)行狀態(tài)監(jiān)測(cè)數(shù)據(jù)以及預(yù)測(cè)結(jié)果來(lái)確定設(shè)備的最佳維修周期，即在故障發(fā)生前維修、更換可能發(fā)生故障的部件;但是目前存在需要對(duì)監(jiān)測(cè)投入較大，對(duì)技術(shù)要求比較高等問題?，F(xiàn)有船舶維修保養(yǎng)體系中還是以事后維修和定期維修為主。定期維修需要解決的關(guān)鍵問題是確定最佳的維修周期，維修周期太長(zhǎng)，導(dǎo)致設(shè)備的可靠性降低，增加了設(shè)備運(yùn)行的風(fēng)險(xiǎn);維修周期太短，又增加了設(shè)備的維修保養(yǎng)的成本，降低了設(shè)備的使用效率。所以確定一個(gè)合理的維修周期，對(duì)降低設(shè)備維修成本、提高設(shè)備使用效率具有重要現(xiàn)實(shí)意義。

國(guó)內(nèi)很多學(xué)者對(duì)維修周期展開了研究，取得了相當(dāng)大的進(jìn)展。宋之杰等［1］考慮了維修成本的動(dòng)態(tài)變化性，構(gòu)建了設(shè)備預(yù)防維修決策模型，該模型以設(shè)備單位周期內(nèi)的總維修成本為目標(biāo)函數(shù)，以設(shè)備可用度為約束條件的預(yù)防維修周期決策模型。用MATLAB仿真軟件對(duì)模型進(jìn)行求解，求解結(jié)果證明了模型的有效性。劉棟等［2］對(duì)現(xiàn)行三級(jí)維修機(jī)構(gòu)保障的復(fù)雜設(shè)備，以規(guī)定可用度為約束條件，以單位工作時(shí)間內(nèi)的平均維修費(fèi)用最低為目標(biāo)，建立了設(shè)備一個(gè)更新周期內(nèi)，維修周期與維修費(fèi)用關(guān)系模型。王靈芝等［3］人在基于以可靠性為中心的預(yù)防性維修計(jì)劃的基礎(chǔ)上分析系統(tǒng)維修費(fèi)用構(gòu)成和系統(tǒng)有效度，并以系統(tǒng)總體維修費(fèi)用最小化、系統(tǒng)有效度最大化為目標(biāo)，建立多部件成組預(yù)防修策略優(yōu)化模型，提出多部件設(shè)備非周期預(yù)防維修計(jì)劃的優(yōu)化方法。陳雅菊［4］以修理經(jīng)費(fèi)和時(shí)間作為目標(biāo)函數(shù)，討論了船舶維修計(jì)劃多目標(biāo)問題;利用了該井的非支配排序遺傳算法對(duì)優(yōu)化模型進(jìn)行了求解;為維修計(jì)劃的制定提供了更多選擇。

對(duì)于可修系統(tǒng)來(lái)說(shuō)，相比較于可靠性，可用性更能體現(xiàn)系統(tǒng)的可靠性和維修性?？尚尴到y(tǒng)可以停機(jī)維修，但出于經(jīng)濟(jì)性和任務(wù)連續(xù)性考慮，一般將嚴(yán)格限定停機(jī)時(shí)間長(zhǎng)短。建模這類系統(tǒng)的可用度的最恰當(dāng)模型是任務(wù)可用度模型。但有關(guān)任務(wù)可用度的研究尚不多見，所以將系統(tǒng)的任務(wù)可用度與維修成本優(yōu)化模型結(jié)合起來(lái)，以船舶主動(dòng)力裝置為研究對(duì)象，分析其任務(wù)可用度。建立以任務(wù)可用度、系統(tǒng)可靠度為約束條件，以維修周期為變量的維修費(fèi)用最低的數(shù)學(xué)模型，同時(shí)利用MATLAB求解模型，并且以典型船舶主動(dòng)力裝置為案例，驗(yàn)證該方法的有效性。

1 維修成本優(yōu)化模型

1.1 模型假設(shè)

設(shè)備的維修屬于一個(gè)目標(biāo)優(yōu)化問題，常見的優(yōu)化目標(biāo)函數(shù)有維修費(fèi)用最低、停機(jī)時(shí)間最短、設(shè)備的維修周期等。作為可修系統(tǒng)，還需要做出以下假設(shè)問題［5］。

1)修復(fù)假設(shè)。修復(fù)假設(shè)可分為修復(fù)如新、修復(fù)如舊以及介于兩者之間的不完全維修。修復(fù)如新表示修復(fù)后設(shè)備的故障率如新設(shè)備;修復(fù)如舊表示修復(fù)后部件的故障率與故障發(fā)生時(shí)的故障率相當(dāng)。

2)系統(tǒng)任務(wù)結(jié)果只有2種情況:失敗與成功。

3)不考慮部件失效間的差異，每個(gè)部件的失效間都相互獨(dú)立。4)不考慮外部環(huán)境、維修方式帶來(lái)的影響。5)系統(tǒng)修復(fù)時(shí)間不考慮故障診斷、備品備件的運(yùn)輸時(shí)間等間接等待時(shí)間。

所以本文構(gòu)建的維修模型要求系統(tǒng)由相互獨(dú)立的部件組成，且特征參數(shù)服從獨(dú)立的指數(shù)分布，以便利用指數(shù)分布的無(wú)記憶性求解模型。設(shè)備無(wú)故障運(yùn)行到預(yù)防性維修周期時(shí)間T，則對(duì)其進(jìn)行一次預(yù)防性維修。若當(dāng)該設(shè)備工作時(shí)間小于預(yù)防維修周期時(shí)間T時(shí)發(fā)生故障，則對(duì)其立刻進(jìn)行修復(fù)性維修，修復(fù)完成后設(shè)備繼續(xù)投入工作直至工作時(shí)間累積到預(yù)防性維修周期，再對(duì)其進(jìn)行一次預(yù)防性維修，如圖1任務(wù)周期所示。

1.2 任務(wù)可用度

圖1 任務(wù)周期

可用度是設(shè)備可靠性與維修性的綜合體現(xiàn)?？捎枚鹊亩x:在規(guī)定的時(shí)間內(nèi)，系統(tǒng)正處于能完成規(guī)定功能狀態(tài)的概率，可用度有多種表現(xiàn)形式，包括:瞬間可用度、穩(wěn)態(tài)可用度、平均可用度以及任務(wù)可用度。假定系統(tǒng)執(zhí)行某次任務(wù)，任務(wù)開始到任務(wù)結(jié)束之間的時(shí)間間隔為一個(gè)任務(wù)周期。系統(tǒng)的任務(wù)結(jié)果只有任務(wù)成功和任務(wù)失敗2種，那么任務(wù)可用度的定義為任務(wù)成功的概率。

國(guó)外學(xué)者針對(duì)任務(wù)可用度進(jìn)一步提出了可修系統(tǒng)的任務(wù)可用度的定義。對(duì)于可修系統(tǒng)，在單次獨(dú)立失效所需的修復(fù)時(shí)間不超過(guò)維修時(shí)限tf條件下，在計(jì)劃運(yùn)行時(shí)間T(維修周期)內(nèi)成功執(zhí)行任務(wù)的概率［6］，即在任務(wù)期間，發(fā)生m次故障，每次故障修復(fù)恢復(fù)時(shí)間都在維修時(shí)限tf內(nèi)的的概率公式如下:

式中:F(T)為系統(tǒng)的累計(jì)失效分布函數(shù); ［Fm(T) －Fm+1(T)］為系統(tǒng)在T內(nèi)發(fā)生m次失效的概率;［G(tf)］m為m次修復(fù)中每次修復(fù)時(shí)間的時(shí)間都小于tf的概率;G(t)為修復(fù)時(shí)間分布函數(shù)。

公式 (1)分為2個(gè)部分，1－F(T)表示在時(shí)間T內(nèi)不發(fā)生故障的概率;后面部分表示在時(shí)間T內(nèi)發(fā)生m次故障，每次故障都在tf時(shí)間內(nèi)修復(fù)完成的概率?？梢岳斫馊蝿?wù)周期中包含2種情形，一種情形是系統(tǒng)在執(zhí)行任務(wù)未發(fā)生故障，此時(shí)完成任務(wù)的概率可以描述為系統(tǒng)可靠性;另外一種情形是系統(tǒng)發(fā)生故障，但修復(fù)恢復(fù)時(shí)間未超過(guò)限定時(shí)間。

1.3 維修成本優(yōu)化模型

以任務(wù)可用度以及可靠度為約束條件，以維修周期為變量，建立總維修成本最低的優(yōu)化目標(biāo)函數(shù)，建立設(shè)備維修成本優(yōu)化模型如下:

式中:f(T)為總維修費(fèi)用;Cw為維修成本;Cl為設(shè)備停機(jī)損失;Cp為預(yù)防維修成本;Cc為故障損失成本;βd為設(shè)備單位時(shí)間停機(jī)成本;A(T，tf)代表維修周期為T，修復(fù)時(shí)間限制tf下的任務(wù)可用度;A0為系統(tǒng)的許用可用度;R為系統(tǒng)的可靠度;R0為許用可靠度。

2 模型求解

所建立的維修費(fèi)用最低數(shù)學(xué)模型屬于非線性規(guī)劃問題。非線性規(guī)劃問題的求解需要將模型轉(zhuǎn)變?yōu)榈臉?biāo)準(zhǔn)形式，多變量約束優(yōu)化問題的標(biāo)準(zhǔn)形式為:

式中:f(x)為目標(biāo)函數(shù);αx為線性不等式約束;ax≤b代表線性規(guī)劃問題的不等式約束條件，b為不等式約束的值;c(x)為非線性不等式約束;ceq(x)為非線性等式約束;Aeqx=beq為線性等式約束;lb和ub為變量x的上下限約束。

所以將所建立的數(shù)學(xué)模型變成標(biāo)準(zhǔn)型。

然后利用MATLAB中自帶的fmincon函數(shù)求解多變量有約束非線性函數(shù)的最小值;得到目標(biāo)函數(shù)的最小值以及最佳維修周期。MATLAB中Fmincon的調(diào)用格式如下:

式中:fun代表目標(biāo)函數(shù);nonlcon表示非線性不等式約束;C(x)≤x和等式Ceq(x) =0的函數(shù)句柄。

所以構(gòu)建相應(yīng)目標(biāo)函數(shù)以及約束條件函數(shù)句柄，利用MATLAB中fmincon即可求解。

3 案例分析

船舶主推進(jìn)系統(tǒng)是船舶的動(dòng)力來(lái)源，它的安全可靠是保證船舶正常航行的前提。根據(jù)船舶主推進(jìn)系統(tǒng)功能圖，得到其系統(tǒng)可靠性框圖;如圖2是一個(gè)典型的船舶動(dòng)力系統(tǒng)可靠性框圖，從圖2可以看出船舶動(dòng)力系統(tǒng)的結(jié)構(gòu)是比較簡(jiǎn)單的，都是一些部件的串、并聯(lián)組合。

圖2 船舶主動(dòng)力裝置可靠性框圖

由船舶主動(dòng)力裝置的可靠性框圖中編號(hào)A－E的部分可以看成并聯(lián)冗余系統(tǒng)，其他部分為串聯(lián)，所以系統(tǒng)的可靠度計(jì)算如下:

主推進(jìn)系統(tǒng)在計(jì)劃運(yùn)行時(shí)間 (0，T］內(nèi)發(fā)生m次故障的概率可由計(jì)數(shù)過(guò)程中泊松過(guò)程求得，以下公式表明在時(shí)間 (t，t+Δt)發(fā)生m次故障的概率:

式中:Λ(t)為累計(jì)故障強(qiáng)度函數(shù)。

假設(shè)系統(tǒng)的失效率服從指數(shù)分布，那么將公式(7)代入公式 (3)中得到任務(wù)可用度計(jì)算公式:

我們假設(shè)22個(gè)部件的可靠性均服從指數(shù)分布，且各部件的失效率查得如表1［7］，維修成本如表2所示。

表1 部件故障率 10－6

維修成本如表2所示。

表2 部件維修維護(hù)成本 103元

βd設(shè)備單位時(shí)間停機(jī)成本為750元/小時(shí);維修時(shí)限為60 h;復(fù)雜機(jī)械系統(tǒng)的許用可靠度、許用可用度一般取0.9;假設(shè)系統(tǒng)的故障率服從指數(shù)分布，失效率為1/4 000。利用MATLAB仿真軟件對(duì)維修成本優(yōu)化模型進(jìn)行求解，約束條件為主推進(jìn)系統(tǒng)的任務(wù)可用度以及可靠度，優(yōu)化變量為系統(tǒng)的計(jì)劃維修周期T。

將目標(biāo)函數(shù)轉(zhuǎn)變?yōu)榉蔷€性規(guī)劃問題的標(biāo)準(zhǔn)形式;利用MATLAB求解優(yōu)化模型得到T=3 488.4 h，維修費(fèi)用隨維修周期關(guān)系如下:

從圖3中，維修成本開始隨維修周期增加而降低;當(dāng)維修周期大于3 488.4 h后，維修成本開始上升;符合實(shí)際規(guī)律。船舶的故障維修保養(yǎng)和傳統(tǒng)機(jī)械維修維護(hù)存在很大的不同，和其他復(fù)雜機(jī)械可以隨時(shí)維護(hù)維修不同的是，船舶設(shè)備的維護(hù)維修需要考慮港口信息。船舶在航行過(guò)程中，拘于備件、技術(shù)條件等原因，很難在航行途中進(jìn)行維護(hù)管理，所以本文得出的最佳維修周期可以作為船舶調(diào)度中權(quán)衡的因素，在既定的航運(yùn)任務(wù)條件下，以最佳維修周期為參考變量，選擇合適的港口來(lái)作為船舶維修維護(hù)的?？奎c(diǎn)。

圖3 系統(tǒng)維修成本與維修周期的關(guān)系

4 結(jié)束語(yǔ)

本文針對(duì)典型的船舶主動(dòng)力裝置系統(tǒng)，采取故障后維修、定期預(yù)防性維修的維修策略，建立了以任務(wù)可用度以及可靠度為約束條件，維修周期與維修成本的關(guān)系模型，最后以算例驗(yàn)證模型及求解方法的適用性和靈敏性，為船舶主動(dòng)力裝置系統(tǒng)最佳維修周期確定提供參考。

［1］宋之杰，楊志秀，趙玉忠，等.可用度及動(dòng)態(tài)維修成本下的維修決策模型 ［J］.工業(yè)工程，2014(2):17－22.

［2］劉棟，李有為，葛陽(yáng).復(fù)雜設(shè)備最佳維修周期優(yōu)化決策模型研究 ［J］.數(shù)學(xué)的實(shí)踐與認(rèn)識(shí)，2013，43(20):177－182.

［3］王靈芝，徐宇工，張家棟.基于設(shè)備有效度和可靠度的預(yù)防修經(jīng)濟(jì)優(yōu)化模型 ［J］.機(jī)械工程學(xué)報(bào)，2010，46(4):163－168.

［4］陳雅菊.基于多目標(biāo)遺傳算法的船舶修理計(jì)劃優(yōu)化研究 ［J］.中國(guó)修船，2011，24(3):46－48.

［5］周瑜.可修系統(tǒng)維修效果評(píng)估及維修決策研究［D］.成都:電子科技大學(xué)，2014.

［6］Attila csenki.Mission Availability For Repairable Semi－Markov Systems:Analytical Results And Computational Implementation［J］.Statistics A Journal of Theoretical＆Applied Statistics，1995，26(1):75 －87.

［7］楊勇虎.船舶動(dòng)力裝置壽命預(yù)測(cè)與風(fēng)險(xiǎn)評(píng)估的研究及系統(tǒng)開發(fā) ［D］.武漢:武漢理工大學(xué)，2003.