趙 陽
(淮南師范學院 機械與電氣工程學院,安徽 淮南 232038)
我國鐵路運量大、消耗多,必須時刻監(jiān)測各信號設(shè)備狀態(tài)和各項性能指標。S700K道岔轉(zhuǎn)轍機從國外引進以來,經(jīng)不斷改進,以其維護少,功率輸出大等優(yōu)點迅速在我國鐵路提速區(qū)段得到應用,成為高速鐵路上控制列車運行方向的關(guān)鍵設(shè)備之一。道岔的轉(zhuǎn)換方向能和聯(lián)鎖控制命令相一致,道岔能完全轉(zhuǎn)換,是列車安全運行的基礎(chǔ),是鐵路運營安全理念的重要體現(xiàn)之一。
目前國內(nèi)外學者對S700K道岔轉(zhuǎn)轍機的研究主要集中在原理介紹[1-2]、電路優(yōu)化[3-4]、三維建模[5]、調(diào)整維護[6-7]等,通過智能算法研究交流電動轉(zhuǎn)轍機的故障診斷方案已有提出。例如,將智能算法包括神經(jīng)網(wǎng)絡(luò)算法[8]、貝葉斯網(wǎng)絡(luò)[9]、灰色關(guān)聯(lián)理論[10]、模糊聚類[11-13]、PSO 優(yōu)化HSMM[14-16]、支持向量機算法[17]等用于轉(zhuǎn)轍機的故障分類。以上智能算法在對S700K道岔轉(zhuǎn)轍機進行故障診斷時,沒有結(jié)合現(xiàn)有的微機監(jiān)測系統(tǒng),且故障診斷過程較為復雜,難以應用到實際現(xiàn)場,存在當診斷網(wǎng)絡(luò)擴大時可能出現(xiàn)指數(shù)爆炸的現(xiàn)象。
弗雷歇距離算法重點研究的是空間時序中兩曲線的路徑距離與兩曲線的相似度關(guān)系。本文提出的弗雷歇距離S700K道岔故障診斷算法省去了對大量樣本數(shù)據(jù)的訓練過程,無需對S700K道岔轉(zhuǎn)轍機動作時的電流、功率或其他工作曲線進行特征提取,診斷過程簡捷,能夠滿足電務(wù)實際對S700K轉(zhuǎn)轍機典型故障的診斷需求。
現(xiàn)階段我國主要采用微機監(jiān)測系統(tǒng)來對各信號設(shè)備的運行狀態(tài)進行遠程監(jiān)測和采集,通過觀測各信號設(shè)備運行時的電流、電壓、功率或其他指標變化情況來判斷其狀態(tài)是否良好。S700K道岔轉(zhuǎn)轍機在轉(zhuǎn)換過程中對道岔尖軌尖端的拉力和輸出功率存在如下關(guān)系:
式中,Re為轉(zhuǎn)轍機傳動系統(tǒng)等效力臂;n為電動機轉(zhuǎn)速;η為電動機效率。
由此可以看出,S700K道岔轉(zhuǎn)轍機的拉力大小能夠直接反映在其運行功率上。因此,S700K運行時的功率曲線為診斷其是否發(fā)生故障以及故障類型提供了判斷條件。
道岔的正常工作過程可分為5個階段:啟動-解鎖-轉(zhuǎn)換-鎖閉-構(gòu)通。正常轉(zhuǎn)換過程中輸出功率隨時間變化的曲線如圖1所示。在該轉(zhuǎn)換過程中,啟動階段是功率需求最大的階段,能夠在2 s內(nèi)使轉(zhuǎn)轍機的輸出功率急速增長到約1.5 kW,如圖1中的a階段。在道岔解鎖、轉(zhuǎn)換和鎖閉過程中,S700K道岔轉(zhuǎn)轍機的輸出功率會迅速下降至約0.5 kW,并在該功率水平保持一段時間,如b階段。道岔鎖閉完成后動作電路會自動被切斷,其功率輸出用于完成道岔位置的表示,此時S700K道岔轉(zhuǎn)轍機的功率降至0.2 kW,如c階段。
當?shù)啦砦恢猛瓿杀硎竞螅浔硎倦娐繁磺袛啵?出功率為0,如d階段。
經(jīng)過對上海鐵路局某段實地調(diào)研并結(jié)合S700K故障檢測的相關(guān)研究成果,將S700K典型故障分成8類,如表1所示,其對應的故障功率曲線如圖2所示。


圖2 S700K典型故障狀態(tài)功率曲線圖

表1 S700K道岔轉(zhuǎn)轍機典型故障列表
進行故障功率曲線與故障原因比對時,設(shè)計步驟如下:
①收集典型的故障曲線并進行聚類:從運營現(xiàn)場獲取S700K轉(zhuǎn)轍機動作功率曲線大數(shù)據(jù)并進行篩選分類;
②構(gòu)建各類故障曲線模板:根據(jù)①的故障曲線建立各類故障的曲線模板,用于對后續(xù)待測曲線進行相似度比對;
③計算弗雷歇距離:收集新的待測故障曲線計算其與各類故障模板之間的弗雷歇距離;
④相似度比較:匹配與待測故障曲線距離最近的故障模板曲線,將其對應的故障類型作為輸出。
Fr é chet distance(弗雷歇距離)是法國數(shù)學家Maurice René Fr é chet在1906年提出的一種路徑空間相似形描述,該理論為描述空間路徑相似度提供了新的思路[18]。弗雷歇距離的提出為描述空間路徑相似度提供了新的思路,使得空間中兩曲線之間相似度的評價效率得到提高。弗雷歇距離的定義如下:
設(shè)L1和L2是度量空間L上的兩條連續(xù)曲線,各滿足式(1)和式(2)。

又設(shè)α和β是單位區(qū)間[0,1]上的兩個重參數(shù)化函數(shù),即式(3)和式(4)。

曲線L1和L2之間的弗雷歇距離F(L1,L2)為式(5)。

其中,α,β∈[0,1]t∈[0,1],d是L上的度量函數(shù)。
基于弗雷歇將空間中兩條曲線之間距離具體化的思想,本文構(gòu)建的兩曲線相似性判斷流程如圖3所示。

圖3 兩曲線相似性判斷流程圖
待測曲線L1和模板曲線L2坐標化是曲線相似性判斷的基礎(chǔ)。如式(6)和式(7)分別為L1和L2各采樣點的坐標集合。

其中,p(n)=(xn,yn),n為曲線L1上采樣點序號,滿足n∈ [1 ,N],xn和yn分別是采樣點的橫坐標和縱坐標。

其中,q(m)=(x*m,y*m),m為曲線L2上采樣點序號,滿足m∈ [1 ,M],mx∈*m[1和,My*]m分別是采樣點的橫坐標和縱坐標。
曲線L1和曲線L2之間的距離矩陣D1為L1和L2上各采樣點之間距離形成的矩陣,如式(8)所示。

L1和L2上兩點間最大距離dmax=max(D1),最小距離為dmin=min(D1)。目標距離f設(shè)置為dmin。并設(shè)置目標距離f的循環(huán)增長間隔r為式(9)。

二值矩陣D2從矩陣D1變換而來。將矩陣D1中小于目標距離f的元素設(shè)置為1,將矩陣D1中大于目標距離f的元素設(shè)置為0即可得到二值矩陣D2,如式(10)。

目標間隔路徑R是否存在決定了在目標距離f下能否找到曲線L1和L2之間的弗雷歇距離。R需滿足:
R 的起點為 d′11,終點為 d′MN;
路徑在通過點d′mn后,其下一個通過點只能為中的一個;
R中所有點的值都必須為1。
弗雷歇距離算法部分代碼實現(xiàn)如圖4所示。

圖4 弗雷歇距離算法部分代碼
弗雷歇算法為兩曲線之間的相似度計算開辟了新研究思路。在對S700K故障曲線分類比較過程中,本文重點關(guān)注相似度的計算,定義相似度Z如式(11)。

其中,F(xiàn)為弗雷歇距離。
通過式(11)可看出,相似度和弗雷歇距離呈現(xiàn)反比例關(guān)系,弗雷歇距離越小,兩曲線間的相似度越高。
在上述故障功率曲線分類比較的過程中,設(shè)定的故障模板曲線分類不同會引起計算的弗雷歇距離和相似度不同。對S700K故障曲線進行分類就是對所得到的弗雷歇距離和相似度進行比較的過程,以此將待測的S700K故障匹配為相似度最高的模板曲線所屬的故障模式。
本實例驗證將表1中的故障類型及圖2中S700K的故障功率曲線作為模板曲線用于弗雷歇距離的計算。
從上海鐵路局某段的微機監(jiān)測系統(tǒng)中選取S700K屬于T5的功率故障曲線作為待測曲線,如圖5所示。使用本文提出的弗雷歇距離算法進行相似度計算,結(jié)果見表2。

圖5 S700K待測故障功率曲線

表2 T5類S700K故障曲線診斷結(jié)果
根據(jù)表2結(jié)果推斷待測故障曲線與T8類故障對應的功率曲線相似度最低,與T5類故障對應的功率曲線相似度最高,從而可以判定此次診斷的S700K故障為T5類故障,與預判故障類型設(shè)定相一致,說明采用弗雷歇距離算法能夠?qū)700K故障進行正確診斷。
在表1中,每類故障曲線各10條,共80條,對本文提出的S700K故障進行批量驗證,診斷結(jié)果如表3所示。

表3 各類型S700K故障曲線診斷結(jié)果
由表3得出,本文提出的采用弗雷歇距離算法進行S700K故障診斷的方法能夠?qū)?類典型故障進行正確診斷。
本文以S700K型轉(zhuǎn)轍機故障為研究對象總結(jié)出8種典型的故障類型并對其對應的功率曲線進行圖表分析。針對S700K轉(zhuǎn)轍機故障診斷準確性欠佳,提出采用計算弗雷歇距離判定兩曲線相似性的方法融入故障診斷過程。最后通過實例驗證,本文提出的故障診斷方法能夠正確識別出故障所屬類型,滿足實際運營維護中對S700K典型故障進行有效診斷識別的需求。不足之處在于受微機監(jiān)測系統(tǒng)自身存儲數(shù)據(jù)的局限性,收集到的典型故障曲線數(shù)量較少,不能夠?qū)θ抗收锨€進行有效匹配診斷。