考慮質量相關性需求與完全預防性維修的EPQ決策

盧 震，徐 健，楊允鋒

(東北大學工商管理學院，遼寧 沈陽 110819)

1 引言

經(jīng)濟生產(chǎn)批量(Economic Production Quantity，EPQ)模型最早由Taft[1]提出，并作為控制庫存的有力工具，其模型被廣泛運用于工業(yè)領域。但傳統(tǒng)的EPQ模型與實際生產(chǎn)并不是十分符合，許多學者結合實際問題對模型進行擴展修正。Massonnet等[2]將需求率設定成以時間為自變量的函數(shù)，進而構建最優(yōu)生產(chǎn)批量模型；在此基礎上，Tayal等[3]和Das等[4]考慮在一個設定的周期內(nèi)，產(chǎn)品的需求率是一個隨時間變化的線性函數(shù)，在此條件下擴展了EPQ模型；Andriolo等[5]進一步考慮了缺貨情況于以上的模型當中；Mandal和Pal[6]構建了基于斜坡型需求率的EPQ模型。Hsu等[7]以及Taleizadeh和Mohammadi[8]對經(jīng)濟生產(chǎn)批量問題的擴展研究進行了相關的文獻綜述；國內(nèi)學者謝金星等[9]針對能力受限的批量問題的數(shù)學模型與算法，介紹了一些主要的研究問題和方法。

由于傳統(tǒng)的EPQ模型對于其所提出的假設過于嚴格，越來越多的學者也開辟了更多新的視角，Ghare和Schrader[10]在構建EPQ模型時假設了一個固定的退化率；Geetha和Udayakumar[11]、Li Shoude[12]等學者均對Ghare和Schrader所提出的模型進行了擴展研究。對于其他方面的考慮，Wee等[13]針對庫存失效的經(jīng)濟生產(chǎn)批量模型研究了部分缺貨的情形；Goyal和Gunasekaran[14]將營銷策略和產(chǎn)品失效引入經(jīng)濟生產(chǎn)批量模型；Teng和Chang[15]構建了一個包括產(chǎn)品價格以及與存貨相關需求的考慮產(chǎn)品失效的經(jīng)濟生產(chǎn)批量模型。國內(nèi)的學者也在EPQ模型的擴展上做了許多相關的研究，魯奎[16]在考慮運輸成本和多供應商的情況下研究了動態(tài)經(jīng)濟批量問題。

此外，為了有效改善生產(chǎn)過程中設備的質量退化現(xiàn)象，許多學者在EPQ問題研究中引入了預防性維修策略，并依據(jù)相鄰兩次預防性維修的間隔期是否固定以及通過預防性維修是否將設備修復如新將預防性維修策略分為周期性完全、周期性不完全、非周期性完全以及非周期性不完全四類預防性維修策略。其中，Aghezza等[17]和Nourelfath等[18]構建了周期性完全預防性維修和生產(chǎn)批量計劃的集成優(yōu)化模型;Fitouhi和Nourelfath[19]和Hui和Gede[20]研究了非周期性完全預防性維修策略下的EPQ問題。Mustapha等[21]在EPQ決策中考慮了周期性不完全預防性策略；Mabrouka等[22]則進一步研究了非周期性不完全預防性策略下的EPQ問題。Charabory和Giri[23]提出了有檢查策略的EPQ模型，在設備生產(chǎn)過程中存在“受控”和“非受控”兩種狀態(tài)，通過定期或不定期的檢查以決定維修或預防性維修來提高設備的可靠性。

由于在大部分擴展EPQ問題均為整數(shù)非線性規(guī)劃問題(INLP)，考慮到問題求解的復雜性，許多啟發(fā)式算法被用于此類問題的求解，遺傳算法(GA)由于其簡單易用性，而成為一種較為常見的求解擴展EPQ問題的啟發(fā)式算法[24-26]。

綜上所述可以看出，已有的EPQ問題的研究中，對于可變需求模式的研究主要集中在庫存依賴需求、時間相關需求等方面，缺乏對產(chǎn)品質量相關需求的考慮，尤其是對于在預防性維修策略下同時考慮質量相關性需求及其與產(chǎn)品質量差異相關聯(lián)的研究并不多見。在實際生產(chǎn)經(jīng)營過程中，產(chǎn)品需求率的大小在影響企業(yè)經(jīng)濟生產(chǎn)批量決策的同時，還會受到產(chǎn)品質量水平的影響；而適當?shù)脑O備維護計劃可以在有效改善設備工作狀態(tài)同時，提高產(chǎn)品質量水平，進而提升產(chǎn)品需求率。因此，在EPQ問題中綜合考慮預防性維修策略與質量相關性需求有助于更為實際的解決EPQ問題。本文在已有研究的基礎上提出單機環(huán)境下考慮產(chǎn)品質量相關性需求的EPQ與維護計劃的集成優(yōu)化模型。同時，為有效降低問題研究的復雜性，本文對設備的維護采用周期性完全預防性維修策略[17,18]。

2 假設與符號說明

2.1 研究問題與相關假設

考慮設備在生產(chǎn)過程中隨著運行時間的增加而發(fā)生退化,設備的運行狀態(tài)隨著時間的推移，會有可能從原來的“受控”狀態(tài)轉為“非受控”狀態(tài)(“非受控”狀態(tài)下設備仍能運轉，但生產(chǎn)的產(chǎn)品質量會大大降低)，設備狀態(tài)轉化的概率隨設備運行時間增加而增加。設備在生產(chǎn)過程中會產(chǎn)出三種質量水平的產(chǎn)品：高質量合格品、低質量合格品和不合格品。其中不合格品只在設備“非受控”狀態(tài)下以一定比率產(chǎn)生，且能被及時檢測并進行修復。設備在“受控”和“非受控”狀態(tài)下均會生產(chǎn)一定比例的低質量合格品，“非受控”狀態(tài)下的低質量合格品產(chǎn)出率高于“受控”狀態(tài)。低質量合格品和高質量合格品雖然均為合格產(chǎn)品，但在未來用戶的使用當中會有不同的質量表現(xiàn)，因此，從長期來看，市場的需求是受產(chǎn)品質量(產(chǎn)品中低質量合格品率)的影響。在這種情況下，通過對設備進行周期性的完全預防性維修可使設備的生產(chǎn)狀態(tài)得到恢復，從而提高設備可靠性，降低設備生產(chǎn)狀態(tài)轉換的概率，減少產(chǎn)品中低質量合格品比率。針對以上問題，本文以單位時間的利潤最大化為目標，綜合考慮質量相關性需求與完全預防性維修下的經(jīng)濟生產(chǎn)批量策略的確定。

為了便于展開研究，結合企業(yè)生產(chǎn)與維護的實際情況，本文給出以下相關研究假設：

1)本文的研究假設由一臺機器生產(chǎn)一種產(chǎn)品，生產(chǎn)率已知且恒定；

2)在生產(chǎn)過程中，設備不發(fā)生導致設備停機的故障；

3)在生產(chǎn)過程中對設備進行檢查，檢查時如果設備處于“非受控”狀態(tài)，則進行完全的恢復，否則進行完全的預防性維修；

4)檢查時間、恢復時間、預防性維修時間均忽略不計；(此項假設主要是為在突出本文研究重點的同時，簡化問題研究的復雜性，此類假設可參見文[17-18])

5)不合格品修復后都是低質量合格品，且其修復過程不占用正常的生產(chǎn)時間；

6)產(chǎn)品的價格由市場決定，為外生變量，產(chǎn)品需求只受產(chǎn)品的質量影響，且在一個生產(chǎn)周期內(nèi)保持不變，產(chǎn)品的質量由產(chǎn)品中的低質量合格品比率來表示。

2.2 符號說明

p單位時間的生產(chǎn)率，已知且恒定；

d單位時間的需求率，是低質量合格品的函數(shù)

Pr單位產(chǎn)品的售價；

n檢查次數(shù)，決策量；

T檢查周期，決策量；

c0單位產(chǎn)品生產(chǎn)成本；

c1生產(chǎn)周期設置成本；

cs單次檢查成本；

ch單位產(chǎn)品單位時間庫存持有成本；

F(t) 生產(chǎn)過程中設備從“受控”到“非受控”發(fā)生狀態(tài)轉換的時間節(jié)點t的概率分布函數(shù)；

R(t) 設備恢復費用，R(t)=a+b(T-t)，其中a>0b>0

cp單次完全預防性維修的費用；

c單位不合格品修復成本；

α“非受控”狀態(tài)下所生產(chǎn)的不合格品比率；

θ1“受控”狀態(tài)下所生產(chǎn)的低質量合格品比率；

θ2“非受控”狀態(tài)下所生產(chǎn)的低質量合格品比率，θ2>θ1；

Q經(jīng)濟生產(chǎn)批量，Q=nTp；

μ需求率調(diào)解系數(shù)。

3 模型構建

針對綜合考慮質量相關性需求與完全預防性維修的EPQ策略問題，本文采用企業(yè)單位時間的利潤最大化為系統(tǒng)最優(yōu)化目標，其中生產(chǎn)總成本包括庫存成本、產(chǎn)品生產(chǎn)成本、檢查成本、周期設置成本、不合格品修復成本、設備預防性維修成本和恢復成本。

3.1 需求率函數(shù)

基于產(chǎn)品的需求率d與其質量有關，當所有產(chǎn)品均為高質量合格品時，產(chǎn)品的需求率達到其最大值dmax。有需求率函數(shù)如下所示：

d=dmax(1-μ·ρ)

(1)

其中，μ為調(diào)解系數(shù)，0<μ≤1；ρ為一次生產(chǎn)周期所生產(chǎn)的低質量合格品的總比率。

一個生產(chǎn)周期中產(chǎn)生的低質量合格品由三部分組成，分別為“受控”狀態(tài)下所生產(chǎn)的低質量合格品、“非受控”狀態(tài)所生產(chǎn)的低質量合格品、不合格品修復后的低質量合格品。

(1)“受控”狀態(tài)生產(chǎn)的低質量合格品數(shù)：

(2)

(2)“非受控”狀態(tài)生產(chǎn)的低質量合格品數(shù)：

(3)

(3)不合格品修復后的低質量合格品的數(shù)量：

(4)

則一個生產(chǎn)周期所產(chǎn)生的低質量合格品總數(shù)為：

A=A1+A2+A3

(5)

一個生產(chǎn)周期所生產(chǎn)的產(chǎn)品總數(shù)為npT，則低質量合格品比率ρ如式所示：

(6)

則與質量相關的需求率函數(shù)如下：

d=dmax(1-μ·ρ)

(7)

3.2 預防性維修成本

設備狀態(tài)轉換時間的概率分布函數(shù)設定為F(t)，其概率密度函數(shù)為f(t)。因為只有在檢查時設備處于“受控”狀態(tài)時才進行預防性維修，因此，當檢查周期為T，每個生產(chǎn)周期的檢查次數(shù)為n，每次完全預防性維修的費用為cp，則可知周期期望預防性維修總成本為：

(8)

3.3 恢復成本

只有在檢查時設備處于“非受控”狀態(tài)時才進行恢復操作。在對設備進行恢復時，單次的恢復費用為R(t)=a+b(T-t)，其中a>0、b>0，T-t為設備發(fā)生狀態(tài)轉換到檢查時的時間間隔，則周期期望恢復總成本為：

(9)

3.4 庫存持有成本

在下面圖1中描述了產(chǎn)品在整個生產(chǎn)過程中的庫存變化，其中T1到Tn為n次預防性維修的時間節(jié)點，其間隔期為T，此期間庫存增加率為p-d，Tn時刻設備停止生產(chǎn)，Tn到Te期間為單純的庫存消耗階段，其消耗率為需求率d，到Te時刻庫存消耗到零，并由此刻開始下一個生產(chǎn)批次。在整個生產(chǎn)過程中，庫存大于零，沒有缺貨發(fā)生。

圖1 產(chǎn)品庫存情況

令單位產(chǎn)品單位時間的庫存持有成本為ch，一個生產(chǎn)周期的期望庫存可以用圖1中三角形面積表示。則期望周期庫存持有總成本如下：

(10)

3.5 不合格品修復成本

當單位不合格品修復成本為c時，由式(4)可知期望周期不合格品修復總成本如下：

(11)

3.6 其他成本

(1)生產(chǎn)周期設置成本

在生產(chǎn)周期開始前，通常發(fā)生各種生產(chǎn)事項的準備費用，設定生產(chǎn)設置成本為一固定的數(shù)值c1，則有：

E(TC)=c1

(12)

(2)周期產(chǎn)品生產(chǎn)成本

令單件產(chǎn)品的生產(chǎn)成本為c0，而一個生產(chǎn)周期的生產(chǎn)批量為npT，則有一個生產(chǎn)周期的周期產(chǎn)品生產(chǎn)成本為：

E(PC)=c0nTp

(13)

(3)周期檢查成本

一個生產(chǎn)周期的檢查次數(shù)為n，設定一次檢查所需要的成本為cs，則有周期檢查成本為：

E(SC)=ncs

(14)

3.7 最優(yōu)EPQ模型的構建

一個生產(chǎn)周期總成本包括設置成本、生產(chǎn)成本、檢查成本、庫存持有成本、預防性維修成本、恢復成本和不合格品修復成本，可以表示為式(15)，其中d如式(7)所示。

(15)

一個生產(chǎn)周期的期望時間為：

(16)

以企業(yè)單位時間的利潤最大化為目標來構建最優(yōu)的EPQ模型，進而求出最優(yōu)的檢查次數(shù)n、檢查周期T，由檢查次數(shù)和檢查周期可確定經(jīng)濟生產(chǎn)批量。企業(yè)一個生產(chǎn)周期的總利潤為總收入減去總成本，總收入為Pr˙npT，總成本如式(15)，單位時間的期望利潤就是周期的期望利潤除以周期的期望時間，由此可得出企業(yè)的單位時間的期望利潤為：

(17)

模型的約束條件如下所示：

(18)

上述約束條件分別表示:設備預防性維修成本小于設備的恢復成本；需求有上限；生產(chǎn)率大于需求率；n是正整數(shù)；單位產(chǎn)品的售價大于單位產(chǎn)品的生產(chǎn)成本。問題的主要決策變量為檢查周期T、檢查或預防性維護的次數(shù)n。此問題是典型的整數(shù)非線性規(guī)劃問題，考慮目標函數(shù)的復雜性，本文采用遺傳算法對模型進行數(shù)值求解。

4 利用遺傳算法求解模型

遺傳算法(GA)是一種模仿自然選擇和生物進化過程的搜索技術。遺傳算法已被廣泛應用于各種應用中，特別是在整數(shù)非線性規(guī)劃問題(INLP)問題中，并被證明能夠在合理的時間內(nèi)提供接近最優(yōu)的解。

結合本文的最優(yōu)EPQ策略問題，具體的算法步驟如下：

1)解表示和種群初始化：種群中的一個解決方案由對應于EPQ策略中的完全預防性維修的間隔周期T和預防性維修次數(shù)n組成的向量來表示。初始種群是通過將隨機生成的T和n的值分配給每個種群個體而產(chǎn)生的。種群的規(guī)模選擇為50。

2)適應度函數(shù)值的計算：由于標準遺傳算法中適應度值是取最小，結合本文研究的目標為利潤最大化，因此，令適應度函數(shù)為公式(17)的結果取負值。并將式(18)中的前三個約束作為懲罰函數(shù)并入到適應度函數(shù)中，形式分別為M*[cp-R(t)]+,M*[d-dmax]+,M*[d-p]+，其中M為懲罰因子，取10000。則適應度函數(shù)的最終表達式為

F=-П+M*[cp-R(t)]++M*[d-dmax]++M*[d-p]+

依據(jù)適應度函數(shù)公式計算種群中每個個體的適應度函數(shù)值。

3)選擇運算：依據(jù)步驟2)所得的50個個體的適應度值，采用輪盤賭選擇法進行50個新個體的選擇。

4)交叉與變異：設定交叉概率和變異概率分別為0.8和0.05。依據(jù)設定的交叉概率和變異概率對步驟3)所得的新個體基因進行交叉和變異操作，得到新一代的子代種群。

5)終止條件判斷：依據(jù)設定的終止代數(shù)50代，判斷是否滿足終止條件？如果不滿足，則重復步驟2到5，直至終止條件滿足。

5 數(shù)值分析

5.1 考慮質量相關需求的EPQ最優(yōu)策略

使用MATLAB遺傳算法工具箱對問題求解，可得最優(yōu)預防性維修次數(shù)n*=3，最優(yōu)的檢查周期T*=0.23月，以及最優(yōu)單位時間利潤П(n*,T*)= 16695.61元/月，遺傳算法最優(yōu)適應度值的迭代過程如圖2所示：

圖2 遺傳算法適應度值迭代過程

在此基礎上，進一步分析不同預防性維修成本cp與不合格品比率θ2對EPQ策略、企業(yè)收益和市場需求的影響，計算結果分別如表1、2所示。

由表1可以看出，當設備處于“非受控”狀態(tài)所產(chǎn)生的不合格品率θ2增加時，市場的質量相關需求保持基本穩(wěn)定，即產(chǎn)品的質量水平基本保持不變，但是最優(yōu)的預防性維修次數(shù)n*增加，維修間隔期T*縮短，單位時間的利潤∏*降低。這是因為當“非受控”狀態(tài)下不合格率增加時，為減少低質量產(chǎn)品的產(chǎn)生，保持穩(wěn)定的市場需求，就需要降低生產(chǎn)過程中由受控狀態(tài)轉為非受控狀態(tài)的可能性，因此所投入的預防性維修次數(shù)增加，且間隔期同時縮短，相應的預防性維修的成本隨之增加，從而導致單位時間利潤有所下降。

表1 不同不合格頻率θ2對最優(yōu)EPQ策略的影響

表2 不同預防性維修成本cp對最優(yōu)EPQ策略的影響

從表2可以看出，當隨著預防性維修的單次成本cp增加時，最優(yōu)的預防性維修次數(shù)n*減少，維修間隔周期T*增大，市場需求d小幅減少，單位時間利潤∏*降低。這是因為當單次預防性維修成本增加時，為有效控制預防性維修總成本，以保證合理的企業(yè)收益，就需要適當減少預防性維修次數(shù)以及加大維修間隔期。這一維修策略的調(diào)整必然在一定程度上增加生產(chǎn)過程中“非受控”狀態(tài)的發(fā)生可能，產(chǎn)品中低質量產(chǎn)品比率因此增加，從而導致市場需求產(chǎn)生小幅降低，單位時間收益相應減少。

由此可見，在市場需求受產(chǎn)品質量影響時，加強對生產(chǎn)過程控制、加強設備維修管理，可以減少生產(chǎn)系統(tǒng)失效狀態(tài)發(fā)生的可能，有效提高產(chǎn)品質量的同時增強產(chǎn)品質量的穩(wěn)定性，為企業(yè)獲取最大化收益。

5.2 與忽略質量相關需求的EPQ策略的對比分析

為進一步驗證考慮質量相關需求EPQ策略的有效性，下面將分析當市場需求受產(chǎn)品質量影響時，采用傳統(tǒng)靜態(tài)需求(即忽略需求的質量相關性)下的EPQ決策與應用考慮質量相關需求的EPQ決策，企業(yè)所獲單位時間利潤之間的關系。首先，在構建基于完全預防性維修的EPQ模型時，若忽略質量相關性需求，則需求為不受產(chǎn)品質量影響的定值，即需求率調(diào)節(jié)系數(shù)μ=0，則d=dmax=900件/月。原問題式(17)在除μ以外其他參數(shù)不變情況下，使用MATLAB遺傳算法工具箱對其進行求解，可得其EPQ最優(yōu)決策為nd*=1,Td*=0.768月，單位時間利潤為∏d(nd*,Td*)= 17451.54元/月(其中上標d表示忽略需求質量相關性下的EPQ決策)。若將此決策應用到實際生產(chǎn)經(jīng)營過程中，即將決策nd*=1，Td*代入到具有質量相關性需求的EPQ模型式(17)中，則產(chǎn)生的實際的單位時間利潤為∏(nd*,Td*)=15654.89元/月，小于原最優(yōu)單位時間利潤∏(n*,T*)=16695.61元/月，其利潤差為1040.72元/月。表3展示了當市場中的需求率調(diào)節(jié)系數(shù)μ在0.1～0.9區(qū)間變化時，企業(yè)考慮質量相關需求所得到的最優(yōu)EPQ決策及產(chǎn)生的單位時間利潤，與忽略質量相關需求下的最優(yōu)EPQ決策及其產(chǎn)生的實際單位時間利潤之間的關系。

表3 不同的需求率調(diào)節(jié)系數(shù)下質量相關性需求EPQ決策與靜態(tài)需求EPQ決策下的利潤比較

首先，由表3左面四列數(shù)據(jù)可以看出，當需求率調(diào)節(jié)系數(shù)μ增加時，產(chǎn)品的需求率受產(chǎn)品質量的影響增大，從保證企業(yè)收益的角度，有必要加大設備預防性維修力度，從而在增加預防性維修次數(shù)n*的同時進一步縮短維修間隔期T*。其次，由表3右面四列數(shù)據(jù)可以看出，忽略需求質量相關性下的EPQ決策不受需求率調(diào)節(jié)系數(shù)變化的影響，而若將此決策應用到實際生產(chǎn)經(jīng)營中，其所產(chǎn)生的企業(yè)收益將小于考慮質量相關性需求的EPQ決策為企業(yè)帶來的收益，且隨著市場需求受產(chǎn)品質量影響程度的增加，這一利潤的差值將明顯增加，由此可以看出本文所給出的EPQ策略會為企業(yè)帶來更多的收益。

6 結語

本文通過在EPQ決策中，綜合考慮質量相關性需求和完全預防性維修策略，對經(jīng)典的最優(yōu)經(jīng)濟生產(chǎn)批量問題進行了擴展性研究。在建模階段，依據(jù)產(chǎn)品需求受產(chǎn)品質量(低質量產(chǎn)品比率)影響這一客觀背景，構建產(chǎn)品的需求率函數(shù)，在此基礎上綜合考慮單批次生產(chǎn)周期內(nèi)發(fā)生的各項成本，以單位時間利潤最大化為目標，構建最優(yōu)經(jīng)濟生產(chǎn)批量模型。考慮到模型的復雜性，采用遺傳算法對預防性維修次數(shù)、維修間隔周期以及經(jīng)濟生產(chǎn)批量的最優(yōu)決策進行數(shù)值求解及分析，并得出以下幾點結論：(1)隨著不合格品率θ2的增加，市場的質量相關性需求基本穩(wěn)定，預防性維修次數(shù)n*增多，維修間隔期T*縮短，單位時間的利潤∏*降低；(2)隨著預防性維修的單次成本cp增加，預防性維修次數(shù)n*減少，維修間隔周期T*增大，市場需求d小幅減少，單位時間利潤∏*降低；(3)隨著市場需求受產(chǎn)品質量影響程度的增加，預防性維修次數(shù)n*增加，維修間隔期T*縮短，單位時間利潤∏*降低；(4)與忽略質量相關性需求的EPQ策略相比，考慮質量相關性需求的EPQ策略會為企業(yè)帶來更多的收益，且隨著市場需求受產(chǎn)品質量影響程度的增加，這一收益的增加變得更為明顯。由于本文采取的是定周期的完全預防性維修策略，故在接下來的研究中可進一步考慮不完全預防性維修且變周期的情況；對于影響需求因素的多樣性方面，可以進一步考慮產(chǎn)品售價、保修期等因素對產(chǎn)品需求的影響。