潘顯俊， 張 煒， 趙 田， 蘇金波

(1. 裝備學(xué)院 光電裝備系， 北京 101416； 2. 裝備學(xué)院 研究生管理大隊， 北京 101416； 3. 裝備學(xué)院 裝備指揮系， 北京 101416)

基于系統(tǒng)動力學(xué)模型的備件供應(yīng)系統(tǒng)計算實驗研究

潘顯俊1,2， 張 煒3， 趙 田2， 蘇金波2

(1. 裝備學(xué)院 光電裝備系， 北京 101416； 2. 裝備學(xué)院 研究生管理大隊， 北京 101416； 3. 裝備學(xué)院 裝備指揮系， 北京 101416)

為了提高備件供應(yīng)系統(tǒng)效率，降低備件供應(yīng)成本，實現(xiàn)備件供應(yīng)“精確化”，以某型武器“基層級-基地級”兩級保障模式為研究對象，分析影響基層級備件供應(yīng)系統(tǒng)備件保障能力的主要因素，建立基層級備件供應(yīng)系統(tǒng)動力學(xué)模型和系統(tǒng)動力學(xué)方程，進(jìn)行了模型校驗以保證模型的有效性。以備件入庫延遲時間、出庫延遲時間、供需差額反應(yīng)時間和供需差額調(diào)整強(qiáng)度為實驗因素，以平均庫存量和供需差額絕對值均值為考察對象進(jìn)行計算實驗。實驗結(jié)果表明：備件供應(yīng)系統(tǒng)中較大的延遲時間和備件供需差額調(diào)整強(qiáng)度將導(dǎo)致備件庫存水平較大波動，較大出庫延遲時間將顯著增加備件庫存水平。

裝備保障；備件供應(yīng)；系統(tǒng)動力學(xué)；計算實驗

優(yōu)良的戰(zhàn)備完好率是裝備維修保障和備件供應(yīng)的重要目標(biāo)。如某型裝備的維修保障模式主要采用“基層級-基地級”兩級維修保障模式，其中基層級維修保障能力直接影響裝備的戰(zhàn)備完好率?；鶎蛹壘S修主要方式是進(jìn)行直接換件維修，不進(jìn)行部件的離位維修。在通常情況下，換件維修耗時短，裝備的不可用時間主要取決于故障檢測時間和備件等待時間，備件供應(yīng)水平成為影響裝備使用可用度和完好率的主要因素[1]。在基層級維修，主要是由備件短缺導(dǎo)致系統(tǒng)的可用度降低，直接影響裝備的戰(zhàn)備完好率，從而影響任務(wù)完成率。因此，科學(xué)合理的基層級備件存儲和供應(yīng)策略是提高裝備完好率的重要保證。

目前，初始備件庫存模型，主要以VARI-METRIC模型為基礎(chǔ)，假設(shè)裝備部件壽命分布為指數(shù)分布，備件需求數(shù)量為泊松分布[2-3]。也有學(xué)者在此模型基礎(chǔ)上，對備件的庫存優(yōu)化方法進(jìn)行了研究[4-5]。但是VARI-METRIC模型沒有考慮在備件供應(yīng)系統(tǒng)中的不同環(huán)節(jié)的延遲動態(tài)變化問題，備件供應(yīng)系統(tǒng)是一個動態(tài)系統(tǒng)，VARI-METRIC模型無法描述其動態(tài)特性。該模型通常適用于裝備論證階段保障系統(tǒng)備件保障方案的初步設(shè)計和裝備部署時初始備件方案設(shè)計。

備件短缺率增加，系統(tǒng)可用度降低，尤其是在基層級維修普遍采用換件維修的模式，備件短缺成為影響裝備可用度的主要因素。備件的供應(yīng)包括訂貨、接收入庫、供應(yīng)出庫和備件使用等環(huán)節(jié)。備件供應(yīng)系統(tǒng)的目標(biāo)就是滿足備件供應(yīng)需求，同時降低備件庫存成本。裝備故障或者計劃維修保養(yǎng)進(jìn)行維修，需要消耗一定的備品備件，從而產(chǎn)生備件需求。在基層級，當(dāng)一次故障或者一次計劃維修需要進(jìn)行換件時，即產(chǎn)生一次備件需求，因此將備件需求率作為備件供應(yīng)系統(tǒng)的輸入。利用系統(tǒng)動力學(xué)模型分析備件供應(yīng)系統(tǒng)，分析備件供應(yīng)系統(tǒng)要素和系統(tǒng)結(jié)構(gòu)，建立備件供應(yīng)系統(tǒng)中的物質(zhì)流和信息流運行機(jī)制，探索備件庫存控制的優(yōu)化策略。

1 兩級備件供應(yīng)系統(tǒng)描述

典型的兩級備件供應(yīng)系統(tǒng)運行過程如圖1所示。LRU(Line-Replaceable Unit)為外場更換單元，SRU(Shop-Replaceable Unit)為內(nèi)場更換單元?；鶎蛹壉Ｕ蠙C(jī)構(gòu)進(jìn)行現(xiàn)場換件維修，不存儲備件。裝備LRU故障時，基層級保障機(jī)構(gòu)向基地級保障機(jī)構(gòu)進(jìn)行備件申請，基地級保障機(jī)構(gòu)供應(yīng)LRU備件；備件送達(dá)時，基層級保障機(jī)構(gòu)進(jìn)行換件維修，并將故障LRU后送到基地級維修保障機(jī)構(gòu)?；鼐S修機(jī)構(gòu)負(fù)責(zé)故障LRU的修理，修理完好后的LRU送到備件倉庫存儲，同時負(fù)責(zé)LRU和SRU訂購。

圖1 典型的兩級備件供應(yīng)系統(tǒng)示意圖

2 備件供應(yīng)系統(tǒng)SD原理

2.1 備件供應(yīng)系統(tǒng)因果分析

備件供應(yīng)系統(tǒng)的功能滿足備件需求，備件供應(yīng)系統(tǒng)根據(jù)備件供應(yīng)完成情況，調(diào)整系統(tǒng)控制參數(shù)，對供應(yīng)活動加以控制，是一個典型反饋控制系統(tǒng)[6-7]。通常描述備件供應(yīng)系統(tǒng)狀態(tài)的要素有：庫存、訂貨、備件出庫、備件故障、備件更換、備件需求、備件供需差額等。備件庫存隨著訂購備件的入庫和備件的出庫發(fā)生變化，出庫使庫存減少，入庫使庫存增加；換件維修時，備件故障需要進(jìn)行更換的速率受制于備件供應(yīng)速率，當(dāng)備件供應(yīng)充足時，故障備件可以得到及時的更換，減少了備件需求的積累；若備件供應(yīng)不足，則故障備件無法得到及時更換，從而增加需求量的積累，此時需要提高備件出庫速率，以滿足備件更換需求，即備件故障需求率的增加將引起備件出庫率的增加；備件供需差額等于備件需求量減去備件庫存量，即備件需求增加引起備件差額增加，庫存增加引起備件差額降低。

將庫存的出庫率作為制定訂貨率的參考依據(jù)之一，同時考慮需求量和供給量的差額信息，與出庫率信息一起作為控制訂貨速率的綜合依據(jù)。庫存的出庫速率增大時，需要增加訂貨率，庫存的出庫率減小時，相應(yīng)的減少訂貨率，從而使庫存水平保持穩(wěn)定。同時，考慮供需矛盾，當(dāng)需求量大于庫存量時，需求供應(yīng)差值大于零，表明供不應(yīng)求，應(yīng)在上一次訂貨率的基礎(chǔ)上增加備件的訂貨率；反之，減小訂貨率，最終達(dá)到供需平衡，即達(dá)到備件的入庫率等于備件的出庫率。

通過以上分析建立備件供應(yīng)系統(tǒng)因果關(guān)系圖，如圖2所示。這是典型反饋控制，利用誤差來消除誤差，從而達(dá)到對系統(tǒng)的控制。根據(jù)反饋控制原理，在利用誤差信息控制系統(tǒng)時，可以通過提高開環(huán)增益的方式，提高控制的精度，但是這樣就有可能引起系統(tǒng)振蕩導(dǎo)致系統(tǒng)不穩(wěn)定。

圖2 備件供應(yīng)系統(tǒng)因果關(guān)系圖

2.2 備件供應(yīng)系統(tǒng)SD模型

圖3 備件供應(yīng)系統(tǒng)SD模型

根據(jù)備件供應(yīng)系統(tǒng)因果關(guān)系圖，利用Anylogic建模仿真工具建立備件供應(yīng)SD流圖，如圖3所示。SD流圖可以更加清晰的刻畫圖2所示的因果關(guān)系，通過建立存量和流速率的數(shù)學(xué)關(guān)系模型，更加具體地反映物質(zhì)流的結(jié)構(gòu)和反饋的量化關(guān)系。同時，加入備件供應(yīng)環(huán)節(jié)的延遲信息，可以更加準(zhǔn)確的描述備件供應(yīng)系統(tǒng)的真實情況。

為了建立各個變量間的數(shù)學(xué)模型方程，定義變量如表1所示。

表1 備件供應(yīng)系統(tǒng)SD模型參數(shù)定義

下面針對各個存量建立系統(tǒng)動力學(xué)方程。

對于訂貨量O(t)，反映在做出訂貨決策后，備件在運輸過程中尚未入庫的備件數(shù)量，滿足

(1)

對于庫存量S(t) ，滿足

(2)

對需求量D(t)，滿足

(3)

考慮入庫延遲情況，訂貨決策到備件入庫之間存在物流延遲。所以有

r(t)=o(t-TRK)

(4)

式(4)表示t時刻的入庫率等于TRK前的訂貨率。

初始時刻，裝備無故障，正常運行，此時無備件需求，備件需求量D(t)=0，備件需求率d(t)=0，備件出庫率c(t)=0，庫存S(t)為初始庫存狀態(tài)。當(dāng)系統(tǒng)發(fā)生故障，產(chǎn)生備件需求，即d(t)>0，產(chǎn)生出庫率c(t)>0。如果考慮管理等原因產(chǎn)生出庫延遲，則有

c(t)=d(t-TCK)

(5)

對換件維修，不考慮換件時間，認(rèn)為備件出庫后直接安裝到系統(tǒng)上，則有

u(t)=c(t)

(6)

至此，得到了備件供應(yīng)系統(tǒng)的基本存量和流速率之間的關(guān)系方程。

作為決策者，庫存策略和訂貨策略應(yīng)該滿足一定的備件滿足率。科學(xué)的思考是，以庫存的出庫率c(t)為反饋信息決定訂貨率，以補(bǔ)充庫存，保證一定的庫存水平。同時，應(yīng)考察需求量D(t)和庫存量S(t)之間的差，即供需差額E(t) ：

E(t)=D(t)-S(t)

(7)

綜合考慮出庫率c(t) 和供需差額E(t) 信息共同決定是否調(diào)整訂貨率o(t) 。合理的方式是，當(dāng)E(t)>0時，即需求量大于供給量，供不應(yīng)求，訂貨率應(yīng)在上一周期的出庫率基礎(chǔ)上有所增加；當(dāng)E(t)<0，即供過于求，則訂貨率應(yīng)在上一周期的出庫率基礎(chǔ)上有所減少。增加或減少的幅度就是調(diào)整強(qiáng)度系數(shù)g，調(diào)整的方向由E(t)的符號決定，g的作用類似反饋系統(tǒng)的開環(huán)增益。即訂貨率為

o(t)=c(t-1)+g·E(t-TE)

(8)

為方便理論分析備件供應(yīng)系統(tǒng)性能，不考慮各項延遲時間。令TE=0、TRK=0、TCK=0，由式(4)～式(6)和式(8)得到

r(t)=o(t)=u(t)=c(t)=d(t)

(9)

o(t)=c(t-1)+g·E(t)

(10)

將c(t)=d(t)代入式(2)得到

(11)

將u(t)=d(t)代入式(3)得到

D(t+1)=D(t)

(12)

將r(t)=o(t)、u(t)=c(t)和式(7)代入式(10)，并保留含有r(t)和d(t)項，得到

r(t+1)=

(13)

將式(11)～式(13)寫成矩陣形式，得到

(14)

即得到以備件需求率d(t)作為系統(tǒng)輸入，以備件需求量D(t)、入庫率r(t)、庫存量S(t)為狀態(tài)空間變量的系統(tǒng)狀態(tài)空間方程。

該離散系統(tǒng)特征方程為

2.3 模型的檢驗

系統(tǒng)動力學(xué)模型的檢驗包括一致性檢驗和有效性，主要方法是量綱一致性檢驗、有效性檢驗和方程式極端條件檢驗[8]。

量綱一致性檢驗——通過分析方程兩端變量量綱進(jìn)行。模型涉及的量綱有存量和速率2類，這里的存量主要是備件的數(shù)量，速率主要是單位時間內(nèi)備件的數(shù)量；通過檢查，量綱一致。

圖4 需求率=2時系統(tǒng)行為

方程式極端條件檢驗——假設(shè)備件需求率為0，初始庫存量為8，通過仿真得到：庫存保持不變，備件訂購率為0，與實際情況一致。假設(shè)庫存為0，備件需求率為4，供應(yīng)差額調(diào)整強(qiáng)度為0，即不考慮供需差額的問題，其他參數(shù)均為0，得到仿真結(jié)果訂購率為4，等于備件的需求率。這與實際情況相符。同樣可以逐一檢驗其他參數(shù)的極端條件，均與實際相符。

通過改變備件的需求率，來驗證系統(tǒng)的有效性。分別令需求率為2，4，初始庫存為8，備件差額調(diào)整強(qiáng)度為1(即差多少就在前次基礎(chǔ)上多訂購多少)，其他參數(shù)為0，進(jìn)行系統(tǒng)仿真。系統(tǒng)開始時備件充足，訂購率為0；隨著備件的消耗，當(dāng)庫存逐漸減少供不應(yīng)求時，訂購率增加，逐漸達(dá)到備件的需求率，如圖4和圖5所示，這與實際庫存情況相符。如果訂購策略是按照供需差額的10倍進(jìn)行，系統(tǒng)仿真結(jié)果見如圖6所示，這與庫存控制模型的理論分析一致[9]，模型通過有效性檢驗。

圖5 需求率=4時系統(tǒng)行為

圖6 訂購策略g=10的庫存變化情況

3 基于備件供應(yīng)SD模型的計算實驗

實驗設(shè)計(DesignofExperiments，DOE)是一種實驗策略，通過有目的的改變一個系統(tǒng)的輸入來考察輸出改變的情況。結(jié)合統(tǒng)計學(xué)方法進(jìn)行實驗設(shè)計，對實驗結(jié)果運用統(tǒng)計分析工具加以分析，最大限度的解釋分析數(shù)據(jù)，發(fā)現(xiàn)各個因素之間的相互關(guān)系，揭示內(nèi)在規(guī)律。在多因素、多水平實驗中，如果考察每個因素的每個水平都相互搭配進(jìn)行全面實驗，則實驗次數(shù)很多，需要耗費大量的實驗成本。通常采用正交實驗，使用正交表Ln(mk)。式中：n是表的行數(shù)，就是實驗方案數(shù)(實驗次數(shù))；k是表的列數(shù)，表示實驗因素數(shù)；m是各個因素的水平數(shù)。

3.1 實驗設(shè)計

假設(shè)初始庫存量S(t)=6，平均備件需求速率為d(t)=4(件/供應(yīng)周期)。對于穩(wěn)定備件供應(yīng)系統(tǒng)，初始備件庫存不影響系統(tǒng)穩(wěn)態(tài)響應(yīng)，即系統(tǒng)經(jīng)過一個時期的自我調(diào)整，會自發(fā)趨于一個穩(wěn)定的狀態(tài)，這是具有負(fù)反饋閉環(huán)系統(tǒng)的特點之一。

針對供需差額調(diào)整系數(shù)g、供需差額反應(yīng)時間TE、入庫延遲時間TRK和出庫延遲時間TCK4個因素，每個因素水平數(shù)為3，則使用正交表L9(34)進(jìn)行。

以供需差額絕對值均值和平均庫存量作為考察指標(biāo)。供需差額的絕對值反映了備件供應(yīng)系統(tǒng)精確性，值越小表示保障越精確。庫存量可以體現(xiàn)庫存成本，庫存越小庫存成本越低，體現(xiàn)了系統(tǒng)的經(jīng)濟(jì)性。

平均庫存量

式中，N為仿真實驗記錄的離散時間點數(shù)。

供需差額絕對值平均值表示系統(tǒng)絕對誤差，它表征了系統(tǒng)精度，是體現(xiàn)“精確保障”的指標(biāo)。

供需差額絕對值平均值

供需差額反應(yīng)時間TE、入庫延遲時間TRK和出庫延遲時間TCK水平選取{0，1，2}，g因素水平取g={0.5，1.0，1.5}。實驗方案與系統(tǒng)仿真結(jié)果如表2和圖7～圖10所示。

表2 備件供應(yīng)系統(tǒng)計算實驗方案及結(jié)果

圖8 方案2仿真結(jié)果

圖9 方案4仿真結(jié)果

圖10 方案6仿真結(jié)果

3.2 實驗結(jié)果與結(jié)論

1) 當(dāng)備件供應(yīng)系統(tǒng)某些環(huán)節(jié)(存在較大延遲時，如出庫延遲時間、入庫延遲時間、供需差額反應(yīng)時間較大等)，庫存水平一直處于波動狀態(tài)，系統(tǒng)存在不穩(wěn)定情況(方案3、6、7、8、9)。如供需差額反應(yīng)時間較大，系統(tǒng)不穩(wěn)定，說明實際利用的供需關(guān)系信息不及時，不能準(zhǔn)確反映供需矛盾，導(dǎo)致備件訂購策略失效。這表明提高庫存管理信息化水平，及時掌握備件供需矛盾，是提高備件供應(yīng)系統(tǒng)能力的有效手段。

2) 供需差額調(diào)整強(qiáng)度系數(shù)過大時，系統(tǒng)可能不穩(wěn)定。說明過分關(guān)注備件短缺的歷史信息，加大備件訂貨力度，矯枉過正，導(dǎo)致庫存積累(方案3、8)。如某些單位根據(jù)出現(xiàn)一次備件短缺的事件，就不顧備件長期以來使用的規(guī)律特點，盲目提高備件訂購數(shù)量，導(dǎo)致訂購該備件數(shù)量過多，庫存積壓。

3) 當(dāng)備件出庫沒有延遲時，供應(yīng)系統(tǒng)經(jīng)過反饋調(diào)節(jié)，使長期穩(wěn)定庫存水平趨于零，實現(xiàn)零庫存供應(yīng)保障，降低庫存成本，提高備件供應(yīng)效費比。反之，備件出庫存在延遲時，會導(dǎo)致備件長期庫存水平處于一個穩(wěn)定值(方案2、4、5)。這個穩(wěn)定庫存水平與出庫延遲時間和入庫延遲時間之間的差成正比(方案4、5)，出庫延遲時間比入庫延遲時間大，導(dǎo)致長期庫存水平增加(因為總有一部分備件需求遲滯在倉庫中，而沒有及時發(fā)放給備件需求單位)。

4) 由于備件供應(yīng)環(huán)節(jié)存在延遲，導(dǎo)致系統(tǒng)不穩(wěn)定，長期庫存水平出現(xiàn)較大波動，且波動幅度逐漸增加。可以通過降低備件供需差額調(diào)整強(qiáng)度，降低庫存水平波動的幅度，甚至可以使庫存水平穩(wěn)定。說明當(dāng)物流延遲嚴(yán)重時，可以通過延長訂貨周期，給基層維修單位一點調(diào)整時間，降低庫存波動水平。

4 結(jié) 束 語

通過備件供應(yīng)系統(tǒng)SD模型的計算實驗，對不同的入庫延遲時間、出庫延遲時間、供需差額反應(yīng)時間和供需差額調(diào)整強(qiáng)度水平，進(jìn)行了仿真分析。本文只初步得出了備件供應(yīng)系統(tǒng)中各個因素對保障效能的影響關(guān)系，給出了相關(guān)定性結(jié)論。在此基礎(chǔ)上，還可以進(jìn)一步在以下方面開展深入研究：(1) 在計算實驗基礎(chǔ)上，應(yīng)用定量分析技術(shù)進(jìn)一步分析各影響因素與保障效能之間的相關(guān)關(guān)系，進(jìn)行靈敏度分析，探索改進(jìn)保障能力的途徑；(2) 以模型為基礎(chǔ)，以備件短缺率和庫存量最小化為優(yōu)化目標(biāo)，研究備件供應(yīng)策略優(yōu)化問題。

References)

[1]丁定浩，陸軍，王斌.新型戰(zhàn)備完好率模型[J].電子產(chǎn)品可靠性與環(huán)境試驗.2010,28(3):1-5.

[2]SHERBROOKE C C.Optimal inventory modeling of systems multi-echelon techniques[M].2nd ed.Boston:Kluwer Academic Publisher,2004:20-31.

[3]李陽.基于裝備戰(zhàn)備完好性的備件配置優(yōu)化研究[D].成都：電子科技大學(xué),2009:32-40.

[4]楊平律，包磊.不確定需求下多級備件庫存系統(tǒng)優(yōu)化模型[J].指揮控制與仿真.2011,33(6):45-49.

[5]楊平律，包磊.基于METRIC的備件庫存優(yōu)化模型[J].艦船電子工程.2010,30(12):161-164.

[6]田瑾，趙廷弟，劉錚.備件供應(yīng)系統(tǒng)的動力學(xué)原理及動態(tài)性研究[J].航空學(xué)報.2007,28(5):1104-1109.

[7]宋揚(yáng).復(fù)雜裝備可維修備件供應(yīng)策略仿真優(yōu)化研究[D].北京：北京科技大學(xué),2015:55-72.

[8]鐘永光，李旭，賈曉菁.系統(tǒng)動力學(xué)[M].北京:科學(xué)出版社,2009:11-20.

[9]尹大力，顏雁，戚龍.庫存控制模型研究[J].長春理工大學(xué)學(xué)報(自然科學(xué)版).2007,30(4):122-124.

(編輯：李江濤)

Research on Computational Experiments of Spare Parts Supply System Based on System Dynamic Model

PAN Xianjun1,2， ZHANG Wei3， ZHAO Tian2， SU Jinbo2

(1. Department of Optical and Electronic Equipment， Equipment Academy, Beijing 101416, China； 2. Department of Graduate Management， Equipment Academy， Beijing 101416， China 3. Department of Equipment Command， Equipment Academy， Beijing 101416， China)

To improve the spare parts supply system efficiency, reduce the spare parts supply cost, and achieve the spare parts supply accuracy, this paper, with the “operation-level to deport-level” support model as its research object, analyses the influence factors of spare parts support ability level of spare parts supply system. It also establishes the operation-level spare parts supply system dynamic model and system dynamic equation, and checks the model to ensure its validity. For the computational experiments, it takes the storage delay time, delivery delay time, supply and demand balance reaction time, supply and demand balance adjustment strength as the experimental factors, and the average inventory and average absolute value of the difference between supply and demand as the object. Experiments show shows that a larger delay time and a larger adjustment intensity of spare parts supply and demand will lead to a larger fluctuation of spare parts inventory level, and a larger time delay will increase the spare parts inventory level significantly.

equipment support; spare parts supply; system dynamics (SD); computational experiments

2017-01-17

國家社會科學(xué)基金資助項目(15GJ005-039);全軍軍事類研究生資助項目(2013JY383)

潘顯俊(1980—)，男，講師，主要研究方向為裝備保障管理。panjun_mail@163.com

E92

2095-3828(2017)03-0051-06

A DOI 10.3783/j.issn.2095-3828.2017.03.009