葛超,賀德強(qiáng),譚文舉
?
基于圖參數(shù)估計(jì)法和IPSO的地鐵列車關(guān)鍵設(shè)備檢修決策模型研究
葛超1,賀德強(qiáng)1,譚文舉2
(1. 廣西大學(xué) 機(jī)械工程學(xué)院,廣西 南寧 530004; 2. 南寧軌道交通集團(tuán)有限責(zé)任公司廣西 南寧 530029)
為了提高地鐵列車關(guān)鍵設(shè)備檢修決策能力,建立服從三參數(shù)威布爾分布,利用圖參數(shù)估計(jì)法和IPSO算法求解可靠度函數(shù),從而達(dá)到優(yōu)化修規(guī)程目的的檢修決策模型。通過(guò)圖參數(shù)估計(jì)法求解得到故障數(shù)據(jù)的三參數(shù)威布爾分布初始估計(jì)值和搜索空間;采用IPSO算法進(jìn)行迭代計(jì)算獲得最優(yōu)參數(shù),并進(jìn)行K-S檢驗(yàn);結(jié)合可靠度衍變規(guī)律和故障數(shù)據(jù)對(duì)關(guān)鍵設(shè)備的檢修決策進(jìn)行優(yōu)化。研究結(jié)果表明:提出的決策模型能夠優(yōu)化關(guān)鍵設(shè)備的修規(guī)程,且IPSO算法適用于解決三參數(shù)威布爾參數(shù)估計(jì)問(wèn)題,為地鐵列車關(guān)鍵設(shè)備的檢修提供新的方法。
地鐵列車;檢修;圖參數(shù)估計(jì)法;IPSO算法;可靠性;三參數(shù)威布爾分布
地鐵列車是城市軌道交通的重要裝備,為了安全完成運(yùn)輸任務(wù),必須使列車設(shè)備處于良好的技術(shù)狀態(tài)。地鐵列車每運(yùn)行一段時(shí)間或相應(yīng)里程后會(huì)執(zhí)行檢修工作,以恢復(fù)故障設(shè)備的功能,保證正線運(yùn)營(yíng)。以可靠性為中心的維修[1](RCM)提高了維修工作的針對(duì)性和效率,受到學(xué)術(shù)界的關(guān)注。馬亮等[2]提出了一種地鐵車輛預(yù)防性檢修多目標(biāo)優(yōu)化模型和求解算法,提高了地鐵車輛檢修工作的效率;王紅等[3]提出了一種基于維護(hù)方式多樣性的地鐵列車設(shè)備預(yù)防性維護(hù)決策模型,該模型能夠制定合理的維護(hù)計(jì)劃,降低維護(hù)成本;呂德峰等[4]建立了地鐵車輛系統(tǒng)檢查周期變化的檢查優(yōu)化模型,車輛可用度和檢修經(jīng)濟(jì)性得以提高;孫超等[5]針對(duì)城市軌道交通車輛設(shè)備維修模式的不足,提出了基于可靠度的最佳維修周期確定方法,該方法未使用優(yōu)化算法求解分布參數(shù)。夏軍等[6]通過(guò)可靠性分析確定了地鐵車門系統(tǒng)危害較大的故障模式,但未提出具體的車門檢修優(yōu)化方法。王瑩等[7]基于列生成算法,根據(jù)動(dòng)車組運(yùn)用維修規(guī)程特點(diǎn),建立了動(dòng)車組檢修計(jì)劃優(yōu)化模型;王靈芝等[8]研究了基于RCM的多部件成組預(yù)防性維修策略,建立了對(duì)應(yīng)的維修費(fèi)用優(yōu)化模型;Doostparast等[9]研究了系統(tǒng)衰退設(shè)備的預(yù)防性維護(hù)問(wèn)題,以保持可靠性及降低維修費(fèi)用為目標(biāo)提出了維修優(yōu)化模型;Tsai等[10]研究了基于可用度的多部分系統(tǒng)的預(yù)防性維護(hù),提出了可靠度模型,該模型基于可靠度最大優(yōu)化維修時(shí)間間隔。目前,關(guān)于地鐵列車RCM的研究相對(duì)較少,針對(duì)列車關(guān)鍵設(shè)備檢修的研究更是罕見(jiàn)。為此,本文提出一種新的檢修決策模型,該模型假設(shè)設(shè)備故障規(guī)律服從三參數(shù)威布爾分布,利用圖參數(shù)估計(jì)法和IPSO算法求解最優(yōu)參數(shù)并進(jìn)行K-S擬合優(yōu)度檢驗(yàn),從而得到故障設(shè)備的可靠度模型,結(jié)合可靠度衍變規(guī)律和故障檢修數(shù)據(jù)對(duì)地鐵列車關(guān)鍵設(shè)備的檢修進(jìn)行優(yōu)化,為地鐵列車檢修科學(xué)決策提供依據(jù)。
由式(1)可得樣本的似然函數(shù)為
對(duì)式(4)取自然對(duì)數(shù)后可得:
圖參數(shù)估計(jì)法是三參數(shù)Weibull分布參數(shù)估計(jì)常用的方法,具有簡(jiǎn)單、直觀、方便等優(yōu)點(diǎn)。本文利用圖參數(shù)估計(jì)法為算法提供解的初始估計(jì)值以及搜索空間。
對(duì)式(3)進(jìn)行2次自然對(duì)數(shù)變換后轉(zhuǎn)換為:
三參數(shù)威布爾分布利用圖參數(shù)估計(jì)法估計(jì)參數(shù)的步驟如下:
求方程組(7)可得基于圖參數(shù)估計(jì)法的威布爾分布3個(gè)參數(shù)的估計(jì)值分別為:
基本PSO(Particle Swarm Optimization)算法是一類隨機(jī)優(yōu)化進(jìn)化算法,其核心思想是模擬鳥群覓食的社會(huì)行為,通過(guò)鳥之間的競(jìng)爭(zhēng)與協(xié)作在飛行空間中尋得食物,即尋找最優(yōu)解的過(guò)程[12?13]。該算法由于具有較少的可調(diào)參數(shù),易編程實(shí)現(xiàn)等特點(diǎn),得到了深入的研究和廣泛應(yīng)用。
IPSO(Improved Particle Swarm Optimization)算法在基本PSO算法的基礎(chǔ)上引入先進(jìn)的理論,極大地提高了算法的收斂性、全局尋優(yōu)性及收斂速度等性能。帶時(shí)變慣性權(quán)重的IPSO算法能夠控制算法的探測(cè)和開發(fā)能力,影響到全局和局部尋優(yōu)性。
本文采取的IPSO算法進(jìn)化方程為:
式中:1和2為加速常數(shù);1和2為服從[0,1]區(qū)間均勻分布的隨機(jī)數(shù);為粒子序數(shù);為當(dāng)前迭代步數(shù);為目標(biāo)搜索空間的維數(shù);x()為第步時(shí)第個(gè)粒子所處的位置;v()為第步時(shí)第個(gè)粒子的進(jìn)化速度;p為第粒子歷史最優(yōu)適應(yīng)度位置;p為全體粒子歷史最優(yōu)適應(yīng)度位置;為慣性權(quán)重。
地鐵列車一般被分成多個(gè)子系統(tǒng),各個(gè)子系統(tǒng)下包含一系列設(shè)備,其中故障頻率較高和故障后果嚴(yán)重的為“關(guān)鍵設(shè)備”,應(yīng)當(dāng)是檢修工作的重點(diǎn)。本文以南寧軌道交通1號(hào)線全部列車故障檢修數(shù)據(jù)為研究對(duì)象。1號(hào)線運(yùn)營(yíng)工作由30列B2型車承擔(dān),該列車共有13個(gè)子系統(tǒng),即:供風(fēng)與制動(dòng)、轉(zhuǎn)向架、車門、車鉤、貫通道、空調(diào)、車體及內(nèi)裝、牽引/制動(dòng)、乘客信息、列車通信控制、主電路、照明子、輔助系統(tǒng)。對(duì)2016?03~2017?03(含試運(yùn)營(yíng))30列列車的故障數(shù)據(jù)進(jìn)行整理和統(tǒng)計(jì),統(tǒng)計(jì)結(jié)果如圖1所示,發(fā)現(xiàn)在全部子系統(tǒng)中,車門系統(tǒng)的故障次數(shù)占比(子系統(tǒng)故障次數(shù)/所有子系統(tǒng)故障次數(shù)之和)最高達(dá)17%,說(shuō)明車門系統(tǒng)具有較高的故障率。線上營(yíng)運(yùn)列車車門出現(xiàn)故障很可能造成列車延誤、掉線甚至嚴(yán)重的安全事故,影響乘車人的舒適度和人身安全。
圖1 各子系統(tǒng)故障數(shù)據(jù)統(tǒng)計(jì)
本文選擇以客室車門關(guān)鍵設(shè)備的運(yùn)營(yíng)故障記錄為故障樣本,進(jìn)行可靠性分析。根據(jù)故障數(shù)據(jù),選取故障次數(shù)較多的設(shè)備進(jìn)行統(tǒng)計(jì),結(jié)果如圖2 所示。
由圖2可見(jiàn),平衡壓輪和門控器發(fā)生故障最為頻繁,選擇其為客室車門的關(guān)鍵設(shè)備,下面以這2種設(shè)備為例進(jìn)行分析,檢驗(yàn)上述決策模型的有效性。對(duì)平衡壓輪壓力和門控器的故障數(shù)據(jù)進(jìn)行圖參數(shù)估計(jì)初步計(jì)算的結(jié)果如表1和表2所示。
圖2 客室車門故障統(tǒng)計(jì)
表1 平衡壓輪故障數(shù)據(jù)
Table 1 Fault data of balance pressure wheels
序號(hào)ti /(d)R(ti)xiyi 1250.979 03.218 9?3.852 6 2290.949 13.367 3?2.951 9 3310.919 23.434 0?2.473 9 4420.889 23.737 7?2.141 9 5440.859 33.784 2?1.886 3 6480.829 33.8712?1.675 7 7560.799 44.025 4?1.496 6 8570.769 54.043 1?1.339 4 9580.739 54.060 4?1.198 1 10600.709 64.094 3?1.069 9 11690.679 64.234 1?0.951 3 12740.649 74.304 1?0.841 1 13800.619 84.382 0?0.737 4 14940.589 84.543 3?0.638 7 15990.559 94.595 1?0.544 7 16990.529 94.595 1?0.454 0 171010.500 04.615 1?0.366 5 181010.470 14.615 1?0.281 3 191060.440 14.663 4?0.197 5
201180.410 24.770 7?0.115 3 211270.380 24.844 2?0.033 5 221280.350 34.852 00.047 8 231280.320 44.852 00.129 4 241380.290 44.927 30.212 3 251420.260 54.955 80.296 5 261450.230 54.976 70.383 6 271480.200 64.997 20.474 0 281680.170 75.124 00.569 8 291720.140 75.147 50.673 5 301770.110 85.176 10.788 5 312050.080 85.323 00.922 6 322860.050 95.656 01.091 2 333220.021 05.774 61.351 5
表2 門控器故障數(shù)據(jù)
根據(jù)表1和表2,繪制故障數(shù)據(jù)的散點(diǎn)圖,如圖3和圖4所示。
圖3 平衡壓輪的故障數(shù)據(jù)WPP及漸近線圖
圖4 門控器的故障數(shù)據(jù)WPP及漸近線圖
分析圖3和圖4,發(fā)現(xiàn)散點(diǎn)并非分布于同一條直線,而是呈曲線波動(dòng)狀,且存在明顯的拐點(diǎn),即存在顯著的凹凸性變化,這與三參數(shù)威布爾分布的特征相契合。由此可以假定故障數(shù)據(jù)符合三參數(shù)威布爾分布。在圖3和4中作直線ASY1和ASY2,分別為散點(diǎn)的鉛垂?jié)u近線和斜漸近線。通過(guò)擬合計(jì)算得到平衡壓輪和門控器的漸進(jìn)線方程:
方程組(13)為平衡壓輪的漸近線,方程組(14)為門控器的漸近線。
平衡壓輪故障數(shù)據(jù)的初始估計(jì)值及搜索空間為:
門控器故障數(shù)據(jù)的初始估計(jì)值及搜索空間為:
根據(jù)IPSO算法編寫程序?qū)δ繕?biāo)函數(shù)進(jìn)行迭代尋優(yōu)計(jì)算,得到平衡壓輪和門控器的迭代過(guò)程,如圖5和6所示。
圖5 平衡壓輪的迭代尋優(yōu)結(jié)果
IPSO算法的具體求解流程如圖7所示。
算法中各參數(shù)取值如下:加速常數(shù)1和2都取2;目標(biāo)搜索空間的維數(shù)取3;迭代次數(shù)取100;慣性權(quán)重采用線性遞減方式,最大值和最小值分別為0.9和0.6;粒子數(shù)為40。
圖6 門控器的迭代尋優(yōu)結(jié)果
圖7 IPSO算法具體流程
平衡壓輪K-S檢驗(yàn)結(jié)果:
門控器K-S檢驗(yàn)結(jié)果:
從式(15)及(16)可以看出,檢驗(yàn)條件完全成立,因此可以判定平衡壓輪和門控器的故障時(shí)間間隔服從三參數(shù)威布爾分布。
由式(3)得:
平衡壓輪可靠度函數(shù)如下:
門控器可靠度函數(shù)如下:
根據(jù)可靠度函數(shù)和故障數(shù)據(jù),在同一張圖中繪制可靠度衍變曲線,如圖8~9所示。
圖8 平衡壓輪故障間隔時(shí)間可靠度函數(shù)圖
分析圖8~9,發(fā)現(xiàn)實(shí)際故障數(shù)據(jù)均勻分布于可靠性函數(shù)曲線的兩側(cè),且變化趨勢(shì)相同。門控器擬合圖中的某些故障點(diǎn)的擬合效果較差,與人為檢修誤差影響了數(shù)據(jù)精度有關(guān)??傮w來(lái)看,這進(jìn)一步證明了地鐵列車關(guān)鍵設(shè)備的故障時(shí)間間隔分布規(guī)律與三參數(shù)威布爾分布相吻合,同時(shí)也表明本文提出的決策模型有效。
可靠度是衡量設(shè)備在規(guī)定條件下實(shí)現(xiàn)預(yù)定功能能力的重要指標(biāo),并且二者之間呈正比例關(guān)系。保證設(shè)備維持在高可靠度重要手段就是進(jìn)行高效的預(yù)防性定期檢修,因此,必須考量檢修策略的合理性。
圖9 門控器故障間隔時(shí)間可靠度函數(shù)圖
本文從平衡壓輪和門控器的可靠度衍變規(guī)律和故障數(shù)據(jù)記錄入手,結(jié)合當(dāng)前的檢修修程和規(guī)程,優(yōu)化地鐵列車關(guān)鍵設(shè)備的檢修策略。平衡壓輪和門控器檢修內(nèi)容如表3所示。
表3 平衡壓輪和門控器檢修內(nèi)容
為了證明IPSO算法對(duì)于解決本文參數(shù)估計(jì)問(wèn)題的適用性,選擇使用另一種群優(yōu)化算法IABC (Improved Artificial Bee Colony)算法對(duì)目標(biāo)函數(shù)進(jìn)行迭代尋優(yōu),所得的計(jì)算結(jié)果如圖10~11所示。
圖10 平衡壓輪故障數(shù)據(jù)優(yōu)化迭代結(jié)果
圖11 門控器故障數(shù)據(jù)優(yōu)化迭代結(jié)果
把2種算法的迭代過(guò)程和結(jié)果進(jìn)行對(duì)比,結(jié)果如表4所示。從表中可以看出,2種算法的求解精度基本相同,但是在收斂速度上迥異。IPSO算法的收斂速度明顯快于IABC算法,就平衡壓輪的最快收斂迭代次數(shù)來(lái)看,二者相差34倍。同時(shí),IPSO算法的收斂速度穩(wěn)定,通過(guò)25~30次迭代即可得到最優(yōu)解,而IABC算法收斂速度極不穩(wěn)定,最快收斂迭代次數(shù)相差750次。較快的求解速度和穩(wěn)定性得益于ISPO算法本身特點(diǎn)及算法參數(shù)的合理選擇。因此,本文提出的IPSO算法用于解決本文的參數(shù)估計(jì)問(wèn)題更具適用性。
表4 IPSO和IABC的迭代效果
1) 三參數(shù)威布爾分布在設(shè)備壽命可靠性分析中具有重要意義,對(duì)于故障數(shù)據(jù)具有強(qiáng)大的擬合能力。結(jié)合圖參數(shù)估計(jì)法和IPSO算法求解威布爾分布三個(gè)最優(yōu)參數(shù),進(jìn)而得到故障設(shè)備的可靠度函數(shù),建立起檢修決策模型。
2) 分析實(shí)際地鐵列車的故障數(shù)據(jù),選取關(guān)鍵設(shè)備,并應(yīng)用決策模型進(jìn)行解算,驗(yàn)證了決策模型的有效性。通過(guò)算法對(duì)比,證明了IPSO算法適用于解決三參數(shù)威布爾參數(shù)估計(jì)問(wèn)題。
3) 根據(jù)關(guān)鍵設(shè)備的可靠度模型和故障數(shù)據(jù),對(duì)列車檢修進(jìn)行優(yōu)化,證明了決策模型具有一定的實(shí)際應(yīng)用價(jià)值,為地鐵列車檢修決策提供了新方法。
[1] 賈希勝. 以可靠性為中心的維修決策模型[M]. 北京:國(guó)防工業(yè)出版社, 2007. JIA Xisheng. The deision models for reliability centred maintenance[M]. Beijing: National Defence Industry Press, 2007.
[2] 馬亮, 郭進(jìn), 張曉霞. 地鐵車輛預(yù)防性檢修計(jì)劃優(yōu)化模型與算法[J]. 鐵道學(xué)報(bào), 2016, 38(7): 8?14. MA Liang, GUO Jin, ZHANG Xiaoxia. Optimation model and algorithm of preventive maintenance schedule for metro rolling stock[J]. Journal of the China Railway Society, 2016, 38(7): 8?14.
[3] 王紅, 劉志龍, 杜維鑫, 等. 基于多種維護(hù)方式的地鐵車輛設(shè)備預(yù)防性維護(hù)策略及維護(hù)模型[J]. 中國(guó)鐵道科學(xué), 2016, 37(4): 108?114. WANG Hong, LIU Zhilong, DU Weixin, et al. Preventive maintenance strategy and model for metro vehicle equipment based on various maintenance methods[J]. China Railway Science, 2016, 37(4): 108?114.
[4] 呂德峰, 過(guò)秀成, 孔哲, 等. 基于可用度的地鐵車輛維修檢查間隔優(yōu)化方法[J]. 東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2011, 41(4): 877?881. Lü Defeng, GUO Xiucheng, KONG Zhe, et al. Maintence inspection interval of metro vehicle based on availability[J]. Journal of Southeast University (Natural Science Edition), 2011, 41(4): 877?881.
[5] 孫超, 徐永能. 城市軌道交通車輛走行部維修周期優(yōu)化模型[J]. 城市軌道交通研究, 2015, 31(1): 71?73. SUN Chao, XU Yongneng. Optimum maintenance circle model for rail transit vehicle’s running gear[J]. Urban Mass Transit, 2015, 31(1): 71?73.
[6] 夏軍, 郭翔, 邢宗義, 等. 地鐵車門系統(tǒng)故障模式可靠性分析[J]. 組合機(jī)床與自動(dòng)化加工技術(shù), 2014, 55(3): 49?52.XIA Jun, GUO Xiang, XING Zongyi, et al. Failure mode critically analysis of metro door system[J]. Modular Machine Tool & Automatic Manufacturing Technique, 2014, 55(3): 49?52.
[7] 王瑩, 劉軍, 苗建瑞. 基于列生成算法的動(dòng)車組檢修計(jì)劃優(yōu)化[J]. 中國(guó)鐵道科學(xué), 2010, 31(2): 115?120.WANG Ying, LIU Jun, MIAO Jianrui. Column generation algorithms based optimation method for maintenance scheduling of muitiple units[J]. China Railway Science, 2010, 31(2): 115?120.
[8] 王靈芝, 徐宇工, 張家棟. 以可靠性為中心的多部件系統(tǒng)預(yù)防性維修優(yōu)化研究[J]. 計(jì)算機(jī)集成制造系統(tǒng), 2008, 126(10): 2042?2046. WANG Lingzhi, XU Yugong, ZHANG Jiadong. Optimization of reliability-centered preventive maintenance for multi-component systems[J]. Computer Integrated Manufacturing Systems, 2008, 126(10): 2042?2046.
[9] Doostparast M, Kolahan F, Doostparast M. A reliability-based approach to optimize preventive maintenance scheduling for coherent systems[J]. Reliability Engineering & System Safety, 2014, 126(126): 98?106.
[10] Tsai Y T, Wang K S, Tsai L C. A study of availa- blity-centered preventive maintenance for multi- componet systems[J]. Reliability Engineering & System Safety 2004, 84(3): 261?270.
[11] 楊謀存, 聶宏. 三參數(shù)Weibull分布參數(shù)的極大似然估計(jì)數(shù)值解法[J]. 南京航空航天大學(xué)學(xué)報(bào), 2007, 39(1): 22?25. YANG Moucun, NIE Hong. Advanced algorithm for maximum likelihood estimation of three parameter weibull distribution.[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2007, 39(1): 22?25.
[12] Kennedy J, Eberhart R. Particle swarm optimization[C]// IEEE International Conference on Neural Networks, 1995. Proceedings. IEEE, 2002: 1942?1948.
[13] Eberhart R, Kennedy J. A new optimizer using particle swarm theory[C]// International Symposium on MICRO Machine and Human Science. IEEE, 2002: 39?43.
[14] 孫妍婷. 動(dòng)車組關(guān)鍵零部件維修決策模型的研究[D]. 北京: 北京交通大學(xué), 2012. SUN Yanting. The research on maintenance decision model of EMU key components[D]. Beijing: Beijing Jiaotong University, 2012.
[15] 麻曉敏, 張士杰, 胡麗琴, 等. 可靠性數(shù)據(jù)威布爾分析中秩評(píng)定算法改進(jìn)研究[J]. 核科學(xué)與工程, 2007, 27(2): 152?155. MA Xiaomin, ZHANG Shijie, HU Liqin et al. An improved rank assessment method for Weibull analysis of reliability data[J]. Chinese Journal of Nuclear Science and Engineering, 2007, 27(2): 152?155.
(編輯 蔣學(xué)東)
Research on overhaul decision model of key equipment of metro train based on graph parameter estimation method and IPSO algorithm
GE Chao1, HE Deqiang1, TAN Wenju2
(1. College of Mechanical Engineering, Guangxi University, Nanning 530004, China; 2. Nanning Rail Transit Co., Ltd, Nanning 530029, China)
For the improvement of key equipment overhaul decision ability of metro trains, the graph parameter estimation method and IPSO algorithm are used to solve the reliability function, and the maintenance decision model which is subject to the three-parameter Weibull distribution is established, achieving the purpose of optimizing repair procedure. Graph parameter estimation method is first used to obtain the initial estimated value and the search space of three parameter Weibull distribution of fault data, then IPSO algorithm is used to get optimum parameters by iterative computation, and then K-S test is employed, at last, key equipment overhaul is optimized by combining evolution law of reliability with fault data. Simulation analysis of actual metro train fault data shows that the decision model can optimize the class of repairs and procedures of key equipment, IPSO algorithm is suitable to solve three parameter Weibull distribution parameter estimation problem and a new method is provided for the overhaul of metro train key equipment.
metro trains; overhaul; graph parameter estimation method; IPSO algorithm; reliability; three parameter Weibull distribution
10.19713/j.cnki.43?1423/u.2018.10.028
TP391.9;U279
A
1672 ? 7029(2018)10 ? 2661 ? 09
2017?07?27
國(guó)家自然科學(xué)基金資助項(xiàng)目(51765006);廣西科技攻關(guān)資助項(xiàng)目(桂科攻1598009-6);南寧市科技攻關(guān)資助項(xiàng)目(20151021);廣西制造系統(tǒng)與先進(jìn)制造技術(shù)重點(diǎn)實(shí)驗(yàn)主任資助項(xiàng)目(15-140-30S003)
賀德強(qiáng)(1973?),男,湖南桃江人,教授,博士,從事機(jī)車車輛故障診斷與智能維護(hù)、列車網(wǎng)絡(luò)與控制、網(wǎng)絡(luò)化制造研究;E?mail:hdqianglqy@126.com