嚴(yán)正峰, 劉　猛, 彭建剛

(1.合肥工業(yè)大學(xué) 1.機(jī)械與汽車工程學(xué)院； 2.汽車工程技術(shù)研究院，安徽 合肥 230009)

?

不確定環(huán)境下復(fù)雜機(jī)械裝配系統(tǒng)設(shè)備維修方法

嚴(yán)正峰1, 劉猛1, 彭建剛2

(1.合肥工業(yè)大學(xué) 1.機(jī)械與汽車工程學(xué)院； 2.汽車工程技術(shù)研究院，安徽 合肥 230009)

摘要：為保障不確定環(huán)境下復(fù)雜機(jī)械裝配系統(tǒng)的連續(xù)性并降低其維修成本，提出一種以馬爾可夫決策理論為基礎(chǔ)的設(shè)備維修策略動(dòng)態(tài)選擇方法。在綜合考慮系統(tǒng)運(yùn)行成本、緩沖庫(kù)存成本、設(shè)備維修成本及停機(jī)損失成本的基礎(chǔ)上，構(gòu)建了裝配系統(tǒng)可靠性成本模型。該模型以帶有中間緩沖區(qū)的二級(jí)裝配系統(tǒng)為研究對(duì)象，以設(shè)備狀態(tài)和緩沖庫(kù)存量為自變量，以可靠性成本為目標(biāo)函數(shù)。分析了裝配系統(tǒng)的不同運(yùn)行狀態(tài)，利用模擬退火算法和模糊非線性混合整數(shù)目標(biāo)規(guī)劃對(duì)可靠性成本模型求解，制定裝配系統(tǒng)最優(yōu)維修方法。該方法降低了裝配系統(tǒng)停機(jī)時(shí)間，減少了設(shè)備維修次數(shù)，可為生產(chǎn)線設(shè)計(jì)和維修計(jì)劃的制定提供依據(jù)。最后，通過(guò)算例分析驗(yàn)證了模型的有效性和可行性。

關(guān)鍵詞：裝配系統(tǒng); 庫(kù)存; 緩沖區(qū); 可靠性成本; 維修方法

在不確定環(huán)境下，各種不確定因素(如設(shè)備磨損、模具更換、緩沖區(qū)的堵塞與饑餓)會(huì)導(dǎo)致復(fù)雜機(jī)械裝配系統(tǒng)出現(xiàn)故障，尤其是關(guān)鍵設(shè)備故障會(huì)給高速運(yùn)行的裝配系統(tǒng)帶來(lái)重大的經(jīng)濟(jì)損失[1]。為保障裝配系統(tǒng)的連續(xù)運(yùn)行，減少單臺(tái)設(shè)備故障造成整個(gè)裝配系統(tǒng)停止運(yùn)行的狀況，可在工作站之間增設(shè)緩沖區(qū)將設(shè)備-設(shè)備之間的“剛性連接”轉(zhuǎn)變?yōu)樵O(shè)備-緩沖區(qū)-設(shè)備的“柔性連接”。以緩沖區(qū)庫(kù)存量和設(shè)備狀態(tài)為依據(jù)制定合理的維修計(jì)劃，能夠有效降低設(shè)備故障率并提高生產(chǎn)效益及裝配系統(tǒng)的穩(wěn)定性。

緩沖區(qū)庫(kù)存量和設(shè)備維修時(shí)機(jī)的動(dòng)態(tài)關(guān)系研究具有一定的理論難度和較高的實(shí)用價(jià)值，為此激發(fā)了一大批學(xué)者開(kāi)展了這一領(lǐng)域的探索，并取得了許多階段性成果。Federgruen和So[2]以設(shè)備長(zhǎng)期運(yùn)行的期望成本最小化為目標(biāo)函數(shù)，確定了單臺(tái)設(shè)備的最佳預(yù)防維修周期；Van der Duyn Schouten和Vanneste等[3]以串聯(lián)裝配系統(tǒng)中的單臺(tái)設(shè)備為研究對(duì)象，考慮了設(shè)備失效率的變化，用嵌入法制定了以設(shè)備服務(wù)時(shí)間和緩沖庫(kù)存量為基礎(chǔ)的預(yù)防維修方法，并證明了該方法的整體最優(yōu)性。Russell D Meller和David S Kim[4]建立了二級(jí)裝配系統(tǒng)模型，利用馬爾可夫鏈確定了最優(yōu)的緩沖庫(kù)存量，通過(guò)建立成本模型確定了設(shè)備的預(yù)防維修時(shí)間及維修后啟動(dòng)生產(chǎn)的時(shí)間。Kyriakidis, Dimitrakos和Karamatsoukis[5-8]等以設(shè)備狀態(tài)和緩沖庫(kù)存量為依據(jù)，分別考慮了設(shè)備維修方法和下級(jí)緩沖區(qū)、上級(jí)緩沖區(qū)庫(kù)存量的關(guān)系，以設(shè)備運(yùn)行成本最低為目標(biāo)函數(shù)確定了設(shè)備修復(fù)后最優(yōu)的等待時(shí)間。Yevkin和Krivtsov[9]以失效率為Weibull分布的設(shè)備為例，建立了設(shè)備的更新函數(shù)，用改進(jìn)的蒙特卡洛方法求解出維修量最小的維修方法。Borrero, Akhavan-Tabatabaei[10]對(duì)設(shè)備狀態(tài)的退化和緩沖庫(kù)存量進(jìn)行了量化，以設(shè)備長(zhǎng)期運(yùn)行的折扣成本最小化為目標(biāo)函數(shù)，建立了兩個(gè)不同的優(yōu)化模型，制定了單臺(tái)設(shè)備的最優(yōu)維修方法。孟吉偉、王少華等[11]研究了緩沖庫(kù)存最小安全量的計(jì)算，提出了非均勻節(jié)拍控制法，其可有效對(duì)緩沖庫(kù)存進(jìn)行補(bǔ)給。韓幫軍、范秀敏等[12-13]提出了等效役齡的概念，建立了預(yù)防維修周期間故障率的遞推關(guān)系式，以裝配系統(tǒng)中的單臺(tái)設(shè)備為對(duì)象建立可靠性成本模型，并用遺傳算法對(duì)模型進(jìn)行了優(yōu)化。付芳[14]使用離散變量構(gòu)建維修計(jì)劃模型，設(shè)計(jì)和聲搜索算法求解模型，對(duì)可行解生成過(guò)程和替換算子進(jìn)行改進(jìn)以保證解的多樣性。通過(guò)算例驗(yàn)證了算法的有效性，并分別對(duì)維修間隔期、故障水平進(jìn)行了敏感性分析以輔助決策。姚運(yùn)志、孟晨[15]引入了關(guān)聯(lián)系數(shù)以描述部件間的失效相關(guān)性，以兩部件系統(tǒng)為例構(gòu)建了系統(tǒng)的失效率函數(shù)，得到了復(fù)雜裝備預(yù)防維修模型的解析形式。綜上所述，緩沖庫(kù)存量的控制對(duì)上下級(jí)設(shè)備的維修方法的制定有重要影響，目前研究大多以單臺(tái)設(shè)備為對(duì)象，本文將同時(shí)考慮緩沖區(qū)和其上下級(jí)設(shè)備的約束關(guān)系，從裝配系統(tǒng)整體的角度考慮設(shè)備運(yùn)行成本和維護(hù)成本，制定具體維修方法。

本文以帶有緩沖區(qū)的二級(jí)裝配系統(tǒng)為研究對(duì)象，并將其抽象為可靠性框圖模型，構(gòu)建裝配系統(tǒng)可靠性成本函數(shù)，利用馬爾可夫決策理論確定裝配系統(tǒng)中設(shè)備狀態(tài)和緩沖庫(kù)存的最優(yōu)匹配。模型能夠根據(jù)設(shè)備工作狀態(tài)和緩沖庫(kù)存量確定最優(yōu)的維修方法并指導(dǎo)現(xiàn)場(chǎng)的維修調(diào)度。

1模型表達(dá)

1.1裝配系統(tǒng)基本模型和狀態(tài)分析

以帶有緩沖區(qū)的二級(jí)裝配系統(tǒng)為研究對(duì)象，不考慮裝配系統(tǒng)上級(jí)緩沖區(qū)的物料短缺和下級(jí)緩沖區(qū)的物料堵塞情況。裝配系統(tǒng)中設(shè)備M1以效率d1生產(chǎn)工件并運(yùn)送至B1儲(chǔ)存，M2以效率d2從緩沖區(qū)B1獲取工件(d1≥d2)，不考慮工件在M1和B1，B1和M2之間的運(yùn)送時(shí)間。詳見(jiàn)圖1。

圖1　裝配系統(tǒng)可靠性框圖模型

隨著設(shè)備M1、M2役齡的增加，其失效率也將增大，可以定期對(duì)設(shè)備狀態(tài)進(jìn)行檢查評(píng)估并將其分為m+2個(gè)狀態(tài)，以0代表設(shè)備的全新?tīng)顟B(tài)，無(wú)任何磨損。1,2,…,m代表設(shè)備磨損程度逐漸增加的過(guò)程，但仍可以正常工作。m+1代表設(shè)備處于失效狀態(tài)，必須進(jìn)行故障維修。設(shè)備由ξ時(shí)刻的狀態(tài)i轉(zhuǎn)移到ξ+1時(shí)刻的狀態(tài)i′的概率為pii′,其中i′∈(i,i+1,i+2,…,m+1), 0

對(duì)設(shè)備可以采取的維修方式u∈{0,1,2}，u=0代表不對(duì)設(shè)備進(jìn)行維修活動(dòng)，使其以現(xiàn)有狀態(tài)繼續(xù)工作，u=1代表對(duì)設(shè)備進(jìn)行預(yù)防維修，u=2代表對(duì)設(shè)備進(jìn)行故障維修。預(yù)防維修和故障維修的維修時(shí)間均服從幾何分布，且其修復(fù)率為定值。當(dāng)設(shè)備處于狀態(tài)1,2，…,m時(shí)，可以對(duì)設(shè)備采取維修方式u=0，也可以采取維修方式u=1，當(dāng)設(shè)備處于狀態(tài)m+1時(shí)，必須采取維修方式u=2即故障維修。

裝配系統(tǒng)的狀態(tài)由設(shè)備M1的狀態(tài)i,緩沖區(qū)B1的庫(kù)存量x和設(shè)備M2的狀態(tài)j決定，即

S=(i,x,j)。

其中，i∈{0,1,2，…,m,m+1,PM}；x∈{0,1,2,…,K}；j∈{0,1,2，…,m,m+1,PM}。

1.2可靠性成本函數(shù)

在構(gòu)建可靠性成本模型時(shí)，將裝配系統(tǒng)的可靠性成本歸納為任意狀態(tài)下的運(yùn)行成本Cij，緩沖區(qū)庫(kù)存成本h*x (單位庫(kù)存成本h和庫(kù)存量x的乘積)，預(yù)防維修成本Cp，故障維修成本Cf和停機(jī)造成的生產(chǎn)損失成本Cl的總和(可靠性成本函數(shù)可參考文獻(xiàn)[5]、[7]、[8]、[13]、[19])，所以可靠性成本函數(shù)為

Cs=Cij+hx+Cp+Cf+Cl。

(1)

當(dāng)裝配系統(tǒng)處于狀態(tài)S=(i,x,j)時(shí)，如對(duì)M1進(jìn)行維修需要考慮M1的停機(jī)損失成本及緩沖庫(kù)存B1不能滿足M2的需求而造成的生產(chǎn)損失，M2維修需要考慮M2的停機(jī)損失成本及緩沖庫(kù)B1達(dá)到最大容量K造成的生產(chǎn)損失。生產(chǎn)損失成本等于設(shè)備單位時(shí)間生產(chǎn)損失量與停機(jī)時(shí)間的乘積。即生產(chǎn)損失成本

Cl=d1t1+d2t2。

(2)

其中，t1表示設(shè)備M1的停機(jī)時(shí)間；t2表示設(shè)M2的停機(jī)時(shí)間；d1、d2分別表示設(shè)備M1、M2的生產(chǎn)率。將式(2)代入式(1)可得

Cs=Cij+hx+Cp+Cf+(d1t1+d2t2)。

(3)

設(shè)備預(yù)防維修和故障維修時(shí)間均服從幾何分布，假設(shè)設(shè)備修復(fù)率為a，則修復(fù)所用時(shí)間t的概率p(t)=a(1-a)(t-1)。

2維修方法分析

Cαn(i,x,j)=min{Cij+hx+α∑m+1i'=i∑m+1j'=jpit'pjj'Cαn(i',x,j')+Cl,Cαn((i,PM),x,(j,PM)}, (4)Cαn(m+1,x,j)=min{Cm+1,j+hx+Cf1+ααCαn[0,min(x-d2t,0),j]+α(1-α) Cαn[m+1,min(x-d2t,0),j]+Cl,Cαn(m+1,x,PM)}, (5)Cαn(i,x,m+1)=min{Ci,m+1+hx+Cf2+αbCαn[i,max(x+d1t,K),j]+α(1-b) Cαn[i,max(x+d1t,K),j]+Cl,Cαn(PM,x,m+1)}。 (6)ì?í????????????

針對(duì)不同維修方法，可對(duì)式(4)～(6)展開(kāi)：

結(jié)合式(4)～(6)，可得出引理1。

引理1：對(duì)任意n=0,1,2,…有

引理1的證明過(guò)程詳見(jiàn)文獻(xiàn)[7]中Lemma 1。

令Cα(i,x,j)為初始狀態(tài)折扣成本的最小期望值，其中(i,x,j)∈S，因?yàn)?/p>

(7)

且存在正數(shù)B滿足：對(duì)所有(i,x,j)∈S且α∈(0,1)有

|Cα(i,x,j)-Cα(i,x*,j)|≤B,x*∈(0,m)。

(8)

證明過(guò)程詳見(jiàn)文獻(xiàn)[6]。

C(i,x,j)=min{Cij+hx+α∑m+1i'=i∑m+1j'=jpii'pjj'C(i',x,j')+Cl-g,C((i,PM),x,(j,PM))}, (9)C(m+1,x,j)=min{Cm+1,j+hx+Cf1+αa'C(0,min(x-d2t,0),j)+α(1-a')+Cl-g+ C(m+1,min(x-d2t,0),j),C(m+1,x,PM)}, (10)C(i,x,m+1)=min{Ci,m+1+hx+Cf2+αb'C(i,max(x+d1t,K),j)+α(1-b')+Cl-g+ C(i,max(x+d1t,K),j),C(PM,x,m+1)}。 (11)ì?í????????

式(7)、(8)表明存在(i,x,j)∈S有常數(shù)g滿足：結(jié)合引理1中(1)～(4)和式(9)～(11)有以下結(jié)論。

即在圖1所示的裝配系統(tǒng)可靠性模型中，設(shè)備M1存在關(guān)鍵狀態(tài)值i*∈{0,1,2，…,m}，有(i*,x,j)∈S滿足此時(shí)對(duì)設(shè)備M1進(jìn)行預(yù)防維修時(shí)的可靠性成本小于裝配系統(tǒng)繼續(xù)維持運(yùn)行的可靠性成本。即(i,x,j)∈S滿足

C(PM,x,j)≤C(i,x,j),0≤i≤m,0≤x≤K,0≤j≤m，

C(PM,x,m+1)≤C(i,x,m+1),0≤i≤m,0≤x≤K。

對(duì)設(shè)備M2有相同結(jié)論，即

C(i,x,PM)≤C(i,x,j),0≤i≤m,0≤x≤K,0≤j≤m，

C(m+1,x,PM)≤C(m+1,x,j),0≤x≤K,0≤j≤m。

3模型的求解

該問(wèn)題是由設(shè)備狀態(tài)轉(zhuǎn)移和維修時(shí)機(jī)選擇這兩個(gè)子問(wèn)題組成的，關(guān)鍵問(wèn)題是確定設(shè)備狀態(tài)值i,j和緩沖庫(kù)存量x之間的最佳匹配問(wèn)題(即求解變量C(i,x,j))，一旦設(shè)備狀態(tài)轉(zhuǎn)移函數(shù)確定，那么維修時(shí)機(jī)選擇問(wèn)題即可轉(zhuǎn)化為三參數(shù)模型求極值問(wèn)題來(lái)解決。這里可以采用遺傳算法作為整體框架來(lái)解決，用模擬退火遺傳算法來(lái)解決Pij,而利用模糊非線性混合整數(shù)目標(biāo)規(guī)劃來(lái)求解可靠性成本的最優(yōu)值，整個(gè)模型的求解流程如圖2所示。

3.1設(shè)備狀態(tài)轉(zhuǎn)移問(wèn)題

從數(shù)學(xué)模型角度出發(fā)，可把模擬退火描述為：在給定領(lǐng)域結(jié)構(gòu)后，模擬退火過(guò)程是從一個(gè)狀態(tài)到另一個(gè)狀態(tài)不斷地隨機(jī)游動(dòng)。結(jié)合本文設(shè)備狀態(tài)轉(zhuǎn)移問(wèn)題，可以用馬爾科夫鏈描述這一過(guò)程。當(dāng)設(shè)備役齡為t時(shí)，設(shè)備狀態(tài)的轉(zhuǎn)移概率定義為

其中，|D|表示狀態(tài)集合中狀態(tài)的個(gè)數(shù); Gij(t)稱為從i到j(luò)的產(chǎn)生概率，其表示設(shè)備由t時(shí)刻的狀態(tài)i轉(zhuǎn)變?yōu)閠+1時(shí)刻狀態(tài)j的概率；Aij(t)為接受概率，Aij(t)表示設(shè)備從狀態(tài)i轉(zhuǎn)移到狀態(tài)j后，接受狀態(tài)j的概率。在某一溫度為Tk和該優(yōu)化問(wèn)題的一個(gè)解xk，可以生成x′，接受x′作為下一個(gè)新解x(k+1)的概率為

P(xk?x′)=

(12)

式(12)中，f(x′)為狀態(tài)x′的目標(biāo)值；f(xk)為狀態(tài)xk的目標(biāo)值。

起始溫度、退溫操作與停止溫度的確定如下。

圖2　模型求解流程

起始溫度的確定：根據(jù)初始解x0，采用2變換法隨機(jī)產(chǎn)生下一個(gè)解x1，基于Metropolis準(zhǔn)則的初始概率來(lái)確定，則起始溫度T0=(f(x1)-f(x0))/lnpo，其中f(x0)和f(x1)表示初始個(gè)體的第一代適應(yīng)度與下一代適應(yīng)度的值；po為初始接受概率，且要求po∈[0.9,0.95]。

退溫操作：Te+1=a×Te，其中，e是迭代步數(shù)；a為降溫系數(shù)，且a∈[0.90,0.95]。

停止溫度：采用基于不改變規(guī)則的控制法，即f(xi+1)-f(xi)=0，就可以認(rèn)為達(dá)到停止溫度。

3.2最小可靠性成本問(wèn)題

如何將問(wèn)題的解編碼成為染色體是遺傳算法使用中的關(guān)鍵問(wèn)題。Gen,Ida和Kim[16]提出了一種用于解決表示為f-nMIGP模型的串并聯(lián)系統(tǒng)可靠性優(yōu)化問(wèn)題的遺傳算法。本文將二級(jí)裝配系統(tǒng)看作f-nMIGP模型的串聯(lián)系統(tǒng)的一部分，在對(duì)遺傳算法進(jìn)行改進(jìn)的基礎(chǔ)上提出多維參數(shù)求極值的方法。具體解決步驟如下。

Step1：設(shè)置參數(shù)和初始化。

設(shè)置種群規(guī)模(pop_size)，變異率(PM)，雜交率(Pc)，最大代數(shù)(max_gen)和初始代數(shù)0。將(i,x,j)作為染色體，其表示設(shè)備狀態(tài)值和緩沖庫(kù)存量的匹配方案，最后選用矩陣的形式來(lái)表征。例如設(shè)備有6個(gè)狀態(tài)值，緩沖庫(kù)存量最大值為8，則可用矩陣M表示。

式中，i表示設(shè)備M1的狀態(tài)值；x表示緩沖區(qū)庫(kù)存量；j表示設(shè)備M2的狀態(tài)值。矩陣M表示設(shè)備M1狀態(tài)值為0時(shí)，緩沖區(qū)庫(kù)存量為1且設(shè)備M2的狀態(tài)值為0。同理，設(shè)備M1狀態(tài)值為1時(shí)，緩沖區(qū)庫(kù)存量為1且設(shè)備M2的狀態(tài)值為1，以此類推，可得出整個(gè)矩陣M。

Step2：適應(yīng)度函數(shù)。

對(duì)于某個(gè)個(gè)體所對(duì)應(yīng)的匹配方案，要判斷其優(yōu)劣，一是看其是否滿足設(shè)備狀態(tài)轉(zhuǎn)移的約束條件；二是計(jì)算其目標(biāo)函數(shù)值，即各個(gè)體所對(duì)應(yīng)的可靠性成本。根據(jù)該問(wèn)題的特點(diǎn)以及初始種群的編碼方式，隱含設(shè)備狀態(tài)值不會(huì)主動(dòng)減少以及緩沖區(qū)庫(kù)存量和上下級(jí)設(shè)備的生產(chǎn)效率保持線性關(guān)系的約束條件。為此，需要對(duì)每個(gè)個(gè)體進(jìn)行判斷，看其是否滿足約束條件，若不滿足，則個(gè)體為不可行解；再計(jì)算其目標(biāo)函數(shù)值。對(duì)于任一個(gè)體(i,x,j),其下一狀態(tài)中不可行解個(gè)數(shù)為Vn,目標(biāo)函數(shù)值為C則該個(gè)體的適應(yīng)度f(wàn)j為

公式中D為定值，由決策者制定，此處令D=10。

Step3：選擇操作。

選擇方法組合了輪盤賭選擇和最優(yōu)性方法。輪盤賭選擇根據(jù)每個(gè)適應(yīng)值占總適應(yīng)值的比例來(lái)產(chǎn)生下一代，最優(yōu)方法為下一代保持最優(yōu)染色體，克服了采樣帶來(lái)的隨機(jī)誤差。根據(jù)適應(yīng)度大小對(duì)種群個(gè)體進(jìn)行選擇。

Step4：交叉操作。

以PC為概率進(jìn)行交叉操作，采用類OX法進(jìn)行交叉，具體步驟如下。

1)隨機(jī)選擇一個(gè)父代交配區(qū)域如下。

2)將N中的交配區(qū)域中放在M的前面，同理也將M中的交配區(qū)域中放在N的前面， 得到M′和N′。

3)去除原有交配位置上的片段得到第二代個(gè)體。

Step5：變異操作。

保證種群的多樣性，對(duì)父代以概率PM進(jìn)行變異，變異概率PM定義了執(zhí)行變異操作的期望染色體數(shù)量PM×pop_size，具體變異方法如下。

隨機(jī)選取若干個(gè)狀態(tài)(如選取2個(gè)狀態(tài))，以最大狀態(tài)值5減去每個(gè)狀態(tài)的值，得到一次變異結(jié)果：

Step6：執(zhí)行終止試驗(yàn)。

如果t

4實(shí)例驗(yàn)證

為了驗(yàn)證以上結(jié)論的正確性，現(xiàn)給出以下實(shí)例進(jìn)行分析驗(yàn)證。以汽車離合器蓋及壓盤總成裝配線中綜合性能檢測(cè)機(jī)和平衡機(jī)構(gòu)成的裝配系統(tǒng)為例，綜合性能檢測(cè)機(jī)對(duì)離合器蓋及壓盤總成進(jìn)行綜合性能檢測(cè)，檢測(cè)完成后半成品通過(guò)托盤運(yùn)送至緩沖區(qū)B1,等待平衡機(jī)對(duì)其進(jìn)行平衡測(cè)試，詳見(jiàn)圖3。

圖3　二級(jí)裝配系統(tǒng)模型

裝配系統(tǒng)已知參數(shù)有

m=10,K=5,g=0,當(dāng)i,j∈(0,m)時(shí)，Cij=0.3(i+j);當(dāng)i=PM或i=m+1時(shí)，Cij=0.3j；當(dāng)j=PM或j=m+1時(shí)，Cij=0.3i;CP1=1.8,a=0.98,Cf1=2.1,a′=0.95,CP2=2.1,b=0.96,Cf2=3.0,b′=0.94。

具體采用的遺傳參數(shù)有：pop_size=20,max_gen=1 500.PC=0.4,PM=0.1, 把以上參數(shù)代入式(8)～(10),在Matlab7.0環(huán)境下對(duì)模型求解，最終計(jì)算結(jié)果如表1所示。

表1　最優(yōu)維修方法

5結(jié)論

針對(duì)傳統(tǒng)制造系統(tǒng)維修方法研究中未全面考慮系統(tǒng)連續(xù)性和可靠性成本的問(wèn)題，提出一種帶有緩沖區(qū)的復(fù)雜機(jī)械裝配系統(tǒng)設(shè)備維修研究方法?？紤]了設(shè)備狀態(tài)基于役齡的隨機(jī)變化，構(gòu)建了基于馬爾可夫理論的可靠性成本模型，并詳細(xì)推導(dǎo)出系統(tǒng)在不同狀態(tài)下的可靠性成本計(jì)算公式。算例分析部分以具體裝配系統(tǒng)為例，通過(guò)改進(jìn)的遺傳算法求出模型最優(yōu)解，即給出了具體裝配系統(tǒng)的最優(yōu)維修方法。驗(yàn)證結(jié)果表明裝配系統(tǒng)維修方法是狀態(tài)閾值問(wèn)題。實(shí)際上，設(shè)備的維修方法可能包含各種周期性大修、小修，但基于可靠性成本的維修模型給出了維修方法的選擇與緩沖庫(kù)存之間的變化關(guān)系。該研究對(duì)于帶有緩沖區(qū)的裝配系統(tǒng)設(shè)備維修計(jì)劃的制定具有一定的指導(dǎo)作用；可靠性成本模型對(duì)降低維修成本，提高裝配系統(tǒng)的生產(chǎn)率和連續(xù)性有實(shí)踐價(jià)值。在后續(xù)的研究工作中，將重點(diǎn)考慮設(shè)備失效率因子和修復(fù)率因子變化情況下的維修方法的制定。

參考文獻(xiàn)：

[1]QI H, ALZAABI R N, WOOD A S, et al. A fuzzy criticality assessment system of process equipment for optimised maintenance management[J]. International Journal of Computer Integrated Manufacturing, 2013, 28(1):112-125.

[2]FEDERGRUEN A, SO K C. Optimal maintenance policies for single-server queueing systems subject to breakdowns[J]. Operations Research, 1990,38(2):330-343.

[3]VAN DER DUYN SCHOUTEN F A, VANNESTE S G. Maintenance optimization of a production system with buffercapacity[J]. European Journal of Operational Research, 1991, 82(2):323-338.

[4]MELLER R D, KIM D S. The impact of preventive maintenance on system cost and buffer size[J]. European Journal of Operational Research, 1996, 95(3):577-591.

[5]DIMITRAKOS T D, KYRIAKIDIS E G. A semi-Markov decision algorithm for the maintenance of a production system with buffer capacity and continuous repair times[J]. General Information, 2008, 111(2):752-762.

[6]KARAMATSOUKIS C C, KYRIAKIDIS E G. Optimal maintenance of a production-inventory system with idle periods[J]. European Journal of Operational Research, 2009, 196(2):744-751.

[7]PAVITSOS A, KYRIAKIDIS E G. Markov decision models for the optimal maintenance of a production unit with an upstream buffer.[J]. Computers & Operations Research, 2009, 36(6):1993-2006.

[8]KARAMATSOUKIS C C, KYRIAKIDIS E G. Optimal preventive maintenance of a production-inventory system when the action of “Idling” is permissible[J]. Springer Optimization & Its Applications, 2014,91:217-218.

[9]YEVKIN O, KRIVTSOV V. Comparative Analysis of Optimal Maintenance Policies Under General Repair With Underlying Weibull Distributions[J]. IEEE Transactions on Reliability, 2013, 62(1):82 - 91.

[10]BORRERO J S, AKHAVAN-TABATABAEI R. Time and inventory dependent optimal maintenance policies for single machine workstations: An MDP approach[J]. European Journal of Operational Research, 2013, 228(3):545-555.

[11]孟吉偉, 王少華, 李小龍. 串行連續(xù)生產(chǎn)線的緩沖庫(kù)存控制與可靠性研究[J]. 機(jī)械設(shè)計(jì)與制造, 2012(5):6-8.

MENG Jiwei, WANG Shaohua, LI Xiaolong. Buffer stock control and reliability of serial continuous production line[J].Machinery Design & Manufacture 2012 (5): 6-8.

[12]韓幫軍, 范秀敏, 馬登哲. 基于可靠度約束的預(yù)防性維修策略的優(yōu)化研究[J]. 機(jī)械工程學(xué)報(bào), 2003,39(6):102-105.

HAN Bangjun, FAN Xiumin, MA Dengzhe. Research on the optimal preventive maintenance policy based on reliability constraints[J]. Chinese Journal of Mechanical Engineering, 2003,39(6):102-105.

[13]韓幫軍, 范秀敏, 馬登哲. 用遺傳算法優(yōu)化制造設(shè)備的預(yù)防性維修周期模型[J]. 計(jì)算機(jī)集成制造系統(tǒng). 2003,9(3):206-209.

HAN Bangjun, FAN Xiumin, MA Dengzhe. Optimization of equipment preventive maintenance cycle in production system based on genetic algorithm[J]. Computer Intergrated Manufactueing Systems. 2003, 9(3):206-209.

[14]付芳. 基于可靠性的機(jī)電設(shè)備維修計(jì)劃優(yōu)化[J]. 計(jì)算機(jī)集成制造系統(tǒng). 2013,19(12):2940-2946.

FU Fang. Maintenance scheduling optimization for electromechanical equipment based on reliability[J]. Computer Intergrated Manufactueing Systems, 2013, 19(12):2940-2946.

[15]姚運(yùn)志, 孟晨, 王成. 考慮失效相關(guān)的多部件系統(tǒng)最優(yōu)預(yù)防維修策略[J]. 計(jì)算機(jī)集成制造系統(tǒng), 2013,19(12):2976-2981.

YAO Yunzhi,MENG Chen,WANG Cheng. Optimal preventive maintenance policies for multi- unit system considering failure interactions[J]. Computer Intergrated Manufactueing Systems, 2013, 19(12):2976-2981.

[16]GEN M, IDA K, KIM J. System reliability optimization with fuzzy goals using genetic algorithm[J]. Japanese Journal of Fuzzy Theory anti Systems, 1998,10(2):356-365.

Equipment Maintenance Method of a Complex Mechanical Assembly System under Uncertain Environment

YAN Zhengfeng1, LIU Meng1, PENG Jiangang2

(1.School of Mechanical and Automotive Engineering, 2.Automobile Engineering & Technology Research Institute, Hefei University of Technology, Hefei 230009, China)

Abstract：To guarantee continuity and decrease maintenance cost of a complex mechanical assembly system under uncertain environment, a dynamic selection method of maintenance policy based on Markov decision is put forward. On the basis of the comprehensive consideration to the cost of operating, inventory, maintenance and the cost due to the lost production, a reliability-cost model is established. The model takes the secondary production system as the research object, equipment condition and the buffer inventory as independent variables, and reliability-cost as objective function. The best maintenance policy in different conditions of production system is derived by simulated annealing algorithm and fuzzy nonlinear programming. The method decreases downtime, reduces equipment maintenance frequency, and provides basis for designing the production line and scheduling the maintenance plan. The calculation analysis of samples proves the effectiveness of the model.

Key words：assembly system; inventory; buffer; reliability cost; maintenance policy

中圖分類號(hào)：TH186

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

文章編號(hào)：1007-7375(2016)01- 0024- 07

doi:10.3969/j.issn.1007- 7375.2016.01.004

作者簡(jiǎn)介：嚴(yán)正峰(1968-)，安徽省人，男，教授，主要研究方向?yàn)樯a(chǎn)系統(tǒng)可靠性.

基金項(xiàng)目：安徽省科技廳工程技術(shù)研究中心建設(shè)計(jì)劃(201106G01015)；安徽省重點(diǎn)實(shí)驗(yàn)室建設(shè)項(xiàng)目

收稿日期：2014- 11- 30