可靠性是反映產(chǎn)品性能的重要指標(biāo)之一, 對于為客戶提供貨物倉儲服務(wù)的物流倉儲企業(yè)而言, 合理評價其服務(wù)過程的可靠性是提升服務(wù)質(zhì)量與客戶滿意度的重要保障。1966年美國軍用標(biāo)準(zhǔn)MIL-STD-721B將可靠性定義為:“產(chǎn)品在規(guī)定的條件下和規(guī)定的時間內(nèi)完成規(guī)定功能的能力?!?span style="font-size:12px;line-height:0;vertical-align:baseline;">[1]可靠性評價就是根據(jù)產(chǎn)品的既有信息, 利用概率統(tǒng)計方法對產(chǎn)品的可靠性特征量進(jìn)行統(tǒng)計與推斷, 包括對產(chǎn)品可靠度的點估計、區(qū)間估計和假設(shè)檢驗等。可靠度 (Reliability) 是用來反映和衡量可靠性高低的尺度, 通常指產(chǎn)品在規(guī)定條件下和規(guī)定時間內(nèi), 完成規(guī)定功能的概率。如果設(shè)λ (t) 為產(chǎn)品在t時刻的失效速率, 則產(chǎn)品的可靠度R (t) 為:
在系統(tǒng)可靠性研究領(lǐng)域, 國內(nèi)外學(xué)者都取得了豐富的研究成果。Lindsay等 (2015) 應(yīng)用可靠性指數(shù)概率分布來衡量電力系統(tǒng)的可靠性。[2]Coolen等 (2015) 就組件在常規(guī)性原因失效條件下, 系統(tǒng)可靠性的非參數(shù)預(yù)測問題展開了深入研究。[3]Chang等 (2014) 建立了一種基于模糊理論的可靠性評價模型, 來評估勞動密集型制造企業(yè)的可靠性。[4]在研究方法上, Billinton等 (2006) 應(yīng)用序貫蒙特卡羅仿真技術(shù)來評價大電力系統(tǒng)的可靠性。[5]Chen等 (2014) 給出了基于Markov狀態(tài)空間和功能分解的分層抽樣蒙特卡羅方法來評價系統(tǒng)的可靠性。[6]李翔等 (2013) 則構(gòu)建了基于Markov鏈的電傳操作系統(tǒng)軟件可靠性評價模型。[7]Majdara等 (2009) 提出了一種基于模塊化建模方法的故障樹自動生成技術(shù), 其研究思路為本文應(yīng)用Petri網(wǎng)模型來進(jìn)行故障樹的仿真求解, 提供了很好的參考與借鑒。[8]
綜合考慮各種研究方法的應(yīng)用特點及論文研究問題的特性, 本文采用具有更好通用性和實用性的故障樹分析法來處理物流倉儲系統(tǒng)的可靠性評價問題, 并依托Petri網(wǎng)模型強(qiáng)大的動態(tài)系統(tǒng)描述能力來進(jìn)行故障樹的仿真求解, 以實現(xiàn)系統(tǒng)可靠度計算的模型化和自動化, 并有效提高系統(tǒng)可靠性評價的科學(xué)性、準(zhǔn)確性。
倉儲系統(tǒng) (Warehousing System, WS) 是指產(chǎn)品分揀或儲存接收過程中使用的設(shè)施、設(shè)備和運作策略的有機(jī)組合, 是物流系統(tǒng)的核心子系統(tǒng)之一。要準(zhǔn)確評價倉儲系統(tǒng)運行的可靠性, 首先應(yīng)全面分析和把握影響其可靠性的各類因素。 (1) 基礎(chǔ)設(shè)施、設(shè)備的可靠性。基礎(chǔ)設(shè)施、設(shè)備是倉儲系統(tǒng)日常生產(chǎn)和運行的基礎(chǔ), 系統(tǒng)功能的發(fā)揮離不開倉庫、貨場、裝卸機(jī)具和運載工具等物理資源的支撐。 (2) 信息系統(tǒng)的可靠性。以倉儲管理系統(tǒng) (Warehousing Management System, WMS) 為代表的計算機(jī)軟件系統(tǒng), 是支撐倉儲系統(tǒng)高效、準(zhǔn)確運行的又一重要保障。WMS可以在充分保證精確度的前提下, 實現(xiàn)倉儲系統(tǒng)有效產(chǎn)出的最大化。 (3) 組織管理的可靠性。主要包括:管理人員的可靠性、現(xiàn)場工作人員的可靠性、管理信息和指令傳達(dá)的可靠性等。 (4) 服務(wù)水平的可靠性。一般體現(xiàn)在倉儲過程的安全性、貨物存取的便捷性, 以及服務(wù)人員的素質(zhì)和態(tài)度等客戶的直接體驗與切身感受, 其主要反映倉儲系統(tǒng)的時間效能和服務(wù)能力。
故障樹分析法 (FTA) 是以一個不希望發(fā)生的事件 (頂事件) 為焦點, 通過自上而下的逐層分析, 找出導(dǎo)致該事件發(fā)生的全部直接原因和間接原因 (基本事件和中間事件) , 建立其邏輯關(guān)系 (門) , 畫出樹狀圖, 并輔以定性分析和定量計算的可靠性評價方法。在FTA中, 若一個集合中的基本事件同時發(fā)生時頂事件必然發(fā)生, 則稱其為割集。在所有割集中, 如果其中的任一基本事件不發(fā)生時頂事件不發(fā)生, 則該割集稱為最小割集。最小割集法是求解故障樹頂事件發(fā)生概率的常用方法。設(shè)Ci (i=1, 2, …, n) 為故障樹的第i個最小割集, 則頂事件X可以表示為X=∪Ci;頂事件發(fā)生的概率為P (X) =P (∪Ci) 。若已知最小割集Ci中各個基本事件x1, x2, …, xk發(fā)生的概率, 則最小割集發(fā)生的概率為P (Ci) =P (∪xj) 。因而, 在已知最小割集發(fā)生概率情況下, 即可求得頂事件發(fā)生概率:
在實際應(yīng)用中, 由于基本事件發(fā)生的概率較小, 只需要保留前三項即可。
通常情況下, 最小割集的階數(shù) (最小割集所含基本事件的數(shù)目) 越低重要性越大。而基本事件的出現(xiàn)次數(shù)越多, 其重要度越高。實際操作過程中, 可以使用下行法 (Fussell-Vesely) 來求取故障樹的最小割集。
以“倉儲系統(tǒng)失效”為頂事件進(jìn)行分析, 其中間事件包括貨物狀態(tài)改變和貨物不能及時出入庫, 中間事件的發(fā)生原因主要有:人為因素、自然災(zāi)害和設(shè)施設(shè)備因素等。造成“貨物狀態(tài)改變”的原因通常有:貨物損壞 (操作不當(dāng)、管理疏失) 、外部環(huán)境因素、設(shè)施設(shè)備因素 (收發(fā)設(shè)施設(shè)備故障、存儲設(shè)備故障、搬運和運輸設(shè)備故障) 和偷盜損失等。貨物不能及時出入庫分為信息傳遞不及時和倉儲能力不足兩種情況, 可能來源于人工操作緩慢、軟硬件故障、儲存空間不足等原因?;谏鲜龇治? 建立倉儲系統(tǒng)可靠性評價的故障樹, 如圖1所示。其中, Ai (i=1, 2, 3, …) 為中間事件, Xj (j=1, 2, 3…) 為基本事件, 各基本事件故障率取值為:X1=0.04%;X2=0.50%;X3=0.30%;X4=0.02%;X5=0.09%;X6=1.18%;X7=0.05%;X8=0.02%;X9=0.04%;X10=1.50%;X11=0.90%;X12=1.5%;X13=1.00%;X14=0.40%;X15=0.03%;X16=0.03%。
對倉儲系統(tǒng)的故障樹采用下行法逐級展開, 具體操作過程見表1。
表1中的最后一列為不能繼續(xù)分解的割集, 但其存在重復(fù), 根據(jù)吸收率可得最小割集為:{X2}、{X3}、{X8}、{X4, X5}、{X4, X6}、{X1, X15}、{X1, X16}、{X9, X10}、{X9, X11}和{X12, X13, X14}。依據(jù)最小割集中的元素情況可知各基本事件的結(jié)構(gòu)重要度次序為:X2=X3=X8>X1=X4=X9>X5=X6=X10=X11=X15=X16>X12=X13=X14>X7??梢? 操作不當(dāng)、管理疏失和盜竊丟失對倉儲系統(tǒng)的可靠性影響較大, 同時故障樹的結(jié)構(gòu) (或門較多) 反映出影響倉儲系統(tǒng)的不穩(wěn)定性因素很多, 因此在產(chǎn)品運行過程中, 應(yīng)加強(qiáng)對此類因素的觀測和檢查, 從而提高系統(tǒng)運行的可靠性。
表1 倉儲系統(tǒng)故障樹下行法展開過程 下載原表
在大型復(fù)雜系統(tǒng)可靠性評價問題中, 使用傳統(tǒng)的FTA來求取最小割集, 計算量極為龐大, 甚至?xí)霈F(xiàn)NP難題。此時, Petri網(wǎng)作為一種特殊的有向網(wǎng), 憑借其優(yōu)越的對象表述能力, 可以用來取代FTA。應(yīng)用Petri網(wǎng)模型進(jìn)行系統(tǒng)可靠性分析的實質(zhì)是:將頂事件作為頂庫所, 逐級向下找出導(dǎo)致頂事件發(fā)生的所有可能因素作為中間庫所和底庫所;將故障樹中的各事件用Petri網(wǎng)中的庫所表示, 故障樹中的各種邏輯門用Petri網(wǎng)中的庫所、變遷和有向弧表示, 即用Petri網(wǎng)的基本元素庫所和變遷的不同連接替代故障樹中的邏輯關(guān)系。當(dāng)故障樹中不含有重復(fù)事件時, 按圖2規(guī)則進(jìn)行轉(zhuǎn)化;若故障樹中含有重復(fù)事件, 可將重復(fù)事件用同一庫所表示, 如圖3所示。
關(guān)聯(lián)矩陣是通過矩陣的形式描述各個點和每條邊之間的關(guān)系。對于無向圖G (V, E) , V表示頂點集合, E表示邊集合, bv, e表示在關(guān)聯(lián)矩陣中的點v和邊e的關(guān)系。若點v和邊e相連, bv, e=1, 否則bv, e=0。如圖4所示, 關(guān)聯(lián)矩陣的第一行{1 1 10}表示頂點v1與邊e1、e2、e3相連, 與邊e4不相連。對于有向圖, 若bv, e=1, 表示邊e由點v發(fā)出, 若bv, e=-1, 表示邊e進(jìn)入點v, 若bv, e=0, 則邊e和點v不相連。
利用關(guān)聯(lián)矩陣, 可將復(fù)雜系統(tǒng)問題的評價過程順利轉(zhuǎn)化為數(shù)學(xué)問題。以矩陣形式來表示Petri網(wǎng)的步驟如下:假設(shè)從庫所P到變遷T的輸入函數(shù)值為非負(fù)整數(shù)w, 記作I (P, T) =w, 并用從P到T的一有向弧和旁注w表示;假設(shè)從變遷T到庫所P的輸出函數(shù)值為非負(fù)整數(shù)w, 記作O (P, T) =w, 則用從T到P的一有向弧和旁注w表示。若w=1, 一般不標(biāo)注;若I (P, T) =0或O (P, T) =0, 則不必畫弧。I與O均可表示為n×m非負(fù)整數(shù)矩陣, O與I的差A(yù)=O-I叫做關(guān)聯(lián)矩陣。對于規(guī)范網(wǎng), w=1。以圖2的Petri網(wǎng)為例, 其相關(guān)矩陣I、O和A分別為:
在關(guān)聯(lián)矩陣中, -1表示有向弧由庫所指向變遷, 1表示由變遷指向庫所。以關(guān)聯(lián)矩陣A為依據(jù), 求取最小割集的Matlab仿真程序設(shè)計如圖5所示。
按照上述方法求得最小割集為:{P1, P4}、{P2, P4}、{P5, P6}、{P7}??煽慷萈10為:
由此可見, Petri網(wǎng)模型的應(yīng)用不僅使圖形簡單明了, 而且算法簡便, 有效的避免了相同符號的出現(xiàn), 還可以借助計算機(jī)仿真軟件來予以實現(xiàn)。
依據(jù)上述轉(zhuǎn)換方法, 將圖1的故障樹模型轉(zhuǎn)化為Petri網(wǎng)結(jié)構(gòu)模型, 對重復(fù)事件和中間事件進(jìn)行等價簡化, 所得結(jié)果如圖6所示。
進(jìn)而得到倉儲系統(tǒng)可靠性評價Petri網(wǎng)模型的輸入矩陣I、輸出矩陣O和關(guān)聯(lián)矩陣A, 即可求得物流倉儲系統(tǒng)的最小割集為:{P2}、{P3}、{P8}、{P4, P5}、{P4, P6}、{P1, P15}、{P1, P16}、{P9, P10}、{P9, P11}和{P12, P13, P14}與故障樹求得最小割集的結(jié)果相同, 說明利用Petri網(wǎng)模型求取最小割集非常準(zhǔn)確, 同時求解過程更為簡便快捷。進(jìn)一步結(jié)合各基本事件的故障率, 最終得到物流倉儲系統(tǒng)的可靠度為0.991 8。
通過構(gòu)建物流倉儲系統(tǒng)故障樹模型, 不僅可以直觀展現(xiàn)產(chǎn)品故障的發(fā)生和發(fā)展過程, 而且有利于在故障發(fā)展的不同階段, 采取相應(yīng)措施和手段阻斷其發(fā)展, 從而達(dá)到預(yù)防事故發(fā)生的目的, 對于保障系統(tǒng)安全可靠運行具有重要意義。論文在全面分析物流倉儲系統(tǒng)可靠性影響因素的基礎(chǔ)上, 建立了系統(tǒng)可靠性評價故障樹, 并利用Petri網(wǎng)仿真技術(shù)對系統(tǒng)的故障樹進(jìn)行仿真建模, 得到系統(tǒng)運行的可靠度, 以此作為物流倉儲系統(tǒng)可靠性預(yù)測和故障診斷及預(yù)防的科學(xué)依據(jù)。研究結(jié)果表明, Petri網(wǎng)仿真模型具有強(qiáng)大的動態(tài)系統(tǒng)描述能力, 將其應(yīng)用于復(fù)雜系統(tǒng)的可靠性研究, 能夠有效提高可靠性評價問題的求解效率。同時論文已經(jīng)證明, Petri網(wǎng)模型與故障樹模型具有很好的等價性, 在處理可靠性評價問題時科學(xué)準(zhǔn)確、簡單易行。
權(quán)所有©:上海陽合儲運
專業(yè)承接上海倉庫租賃、上海倉儲配送物流、上海電商倉儲企業(yè)服務(wù)與微笑同在"的先進(jìn)理念不斷發(fā)展壯大。