王遠(yuǎn)霏， 孫海榮， 裴春興， 陳永春， 劉先升

(中車唐山機(jī)車車輛有限責(zé)任公司 產(chǎn)品技術(shù)研究中心， 河北唐山 063035)

軸承故障診斷技術(shù)是高速列車安全運(yùn)行的關(guān)鍵技術(shù)之一。軸承作為動(dòng)車中重要的旋轉(zhuǎn)零件，也是主要故障源之一[2]。隨著列車?yán)锍痰脑黾樱S箱軸承也不斷受到磨損，特別是由于安裝質(zhì)量不佳，潤(rùn)滑油品質(zhì)不好，外界異物進(jìn)入時(shí)，軸承內(nèi)部會(huì)出現(xiàn)剝離、卡傷、裂紋、滾子碎裂等現(xiàn)象[3]。對(duì)于高速行駛的列車，做好軸承運(yùn)轉(zhuǎn)狀態(tài)的監(jiān)測(cè)與故障診斷能有效提高高速鐵路的運(yùn)行管理水平，并具有顯著的經(jīng)濟(jì)效益。

國(guó)內(nèi)外有關(guān)動(dòng)車組軸箱軸承故障的研究有許多。代宏偉[4]針對(duì)快速鐵路動(dòng)車軸承故障檢測(cè)的不確定性，提出了引入權(quán)重分析優(yōu)化的快速鐵路動(dòng)車軸承故障檢測(cè)方法，以人工神經(jīng)網(wǎng)絡(luò)和專家系統(tǒng)相結(jié)合，引入權(quán)重分析優(yōu)化的概念對(duì)動(dòng)車軸承參數(shù)權(quán)重進(jìn)行優(yōu)化。李佳睿等[5]針對(duì)動(dòng)車組走行部軸箱振動(dòng)信號(hào)非平穩(wěn)、非線性且背景噪聲頻率分布范圍廣的復(fù)雜特性，設(shè)計(jì)了改進(jìn)的希爾伯特-黃變換算法(Hilbert-Huang Transform, HHT)，綜合利用了HHT方法有效平穩(wěn)化信號(hào)及共振解調(diào)方法有效提高信噪比的優(yōu)點(diǎn)，對(duì)軸箱軸承故障信號(hào)進(jìn)行分析。李永華等[6]綜合考慮影響軸箱軸承壽命的模糊因素，運(yùn)用集對(duì)分析理論和模糊集理論對(duì)軸箱軸承進(jìn)行可靠性壽命評(píng)估。楊天蘊(yùn)[7]建立了故障軸承動(dòng)力學(xué)模型，采用數(shù)值仿真對(duì)不同故障特征的軸承振動(dòng)特性進(jìn)行分析，并模擬軸承的實(shí)際工作環(huán)境設(shè)計(jì)了試驗(yàn)，為轉(zhuǎn)向架軸承的故障檢測(cè)提供理論依據(jù)。郝燁江[8]研究了軸承接觸疲勞壽命的評(píng)估方法，設(shè)計(jì)了軸箱軸承可靠性試驗(yàn)方法。韓威[9]研究了基于EEMD和小波包的動(dòng)車組軸箱軸承故障診斷系統(tǒng)。汪志君[10]研究了基于EEMD的共振解調(diào)技術(shù)在CRH380BL動(dòng)車組軸箱軸承故障診斷中的應(yīng)用，并通過(guò)仿真證明了該方法適用于動(dòng)車組軸箱軸承的早期故障診斷。

首先對(duì)重點(diǎn)故障進(jìn)行關(guān)聯(lián)分析和故障診斷，挖掘故障間的關(guān)聯(lián)關(guān)系，探尋故障與外界因素的關(guān)系，結(jié)合重點(diǎn)故障類型，建立軸箱軸承故障模式庫(kù)、故障特征庫(kù)、故障特性經(jīng)驗(yàn)庫(kù)等知識(shí)庫(kù)，同時(shí)研究了軸箱軸承故障演變規(guī)律和故障間隙，建立安全性能狀態(tài)的置信區(qū)間和安全健康指標(biāo)閾值。根據(jù)滾動(dòng)軸承的發(fā)熱機(jī)理、傳熱過(guò)程及軸承溫度場(chǎng)分布等利用RBF 神經(jīng)網(wǎng)絡(luò)對(duì)軸箱軸承的溫度變化趨勢(shì)進(jìn)行建模。最終結(jié)合平臺(tái)數(shù)據(jù)倉(cāng)庫(kù)中存儲(chǔ)的歷史數(shù)據(jù)和實(shí)時(shí)監(jiān)控?cái)?shù)據(jù)，采用深度機(jī)器學(xué)習(xí)、模式匹配、神經(jīng)網(wǎng)絡(luò)等技術(shù)，結(jié)合已有的專家系統(tǒng)知識(shí)庫(kù)和故障樹，挖掘潛在故障關(guān)聯(lián)關(guān)系，預(yù)測(cè)潛在故障，指導(dǎo)作業(yè)人員提前進(jìn)行維修。

1 故障分析

經(jīng)過(guò)分析2011—2017年某鐵路局CRH380B(L)型動(dòng)車重點(diǎn)故障數(shù)據(jù)庫(kù)，可知2011年7月6日至2017年5月7日，軸箱軸承故障共發(fā)生了335起，其中由于溫度升高引起的故障有295起，軸箱軸承滲油故障29起，軸箱軸承異聲/TADS報(bào)警11起。由溫度升高引起的295起故障中，屬于1級(jí)故障(影響運(yùn)行安全)的有1起；屬于2級(jí)故障(影響車組運(yùn)行)的有11起；屬于3級(jí)故障(影響系統(tǒng)運(yùn)行)的有64起；屬于4級(jí)故障(一般常見故障)的有50起；屬于嚴(yán)重故障的有4起。

從以上分析可以看出，軸箱軸承的故障主要是溫度異常，因此，軸箱軸承故障數(shù)據(jù)挖掘模型的關(guān)鍵是刻畫軸箱軸承溫度變化規(guī)律和影響因素的關(guān)系，在軸箱軸承溫度預(yù)測(cè)的基礎(chǔ)上進(jìn)行故障預(yù)測(cè)。

目前針對(duì)動(dòng)車組軸溫升高的處理機(jī)制如下：

運(yùn)行途中發(fā)生軸承異常溫升的動(dòng)車組入庫(kù)時(shí)需根據(jù)動(dòng)車組軸承溫度(t)或者軸承溫升(Δt)進(jìn)行重點(diǎn)檢查[4]。

1.1 符合以下條件之一的，需入庫(kù)開蓋檢查軸承

(1)換輪后試運(yùn)行動(dòng)車組：47℃≤Δt<65℃或110℃≤t<120℃。

(2)換輪后運(yùn)行5萬(wàn)km以內(nèi)動(dòng)車組：40℃≤Δt<47℃或7日內(nèi)有3日同一軸承出現(xiàn)110℃≤t<115℃。

(3)換輪后運(yùn)行超過(guò)5萬(wàn)km動(dòng)車組：105℃≤t<115℃或30℃≤Δt<45℃(裝用美爾森接地裝置的輪對(duì)在每年的2月至6月期間按照25℃≤Δt<45℃執(zhí)行)。

(4)若檢查無(wú)異常，安排再次試運(yùn)行或安排局管內(nèi)短交路或?？空据^多交路運(yùn)行至Δt<45℃，溫度有下降趨勢(shì)，可安排長(zhǎng)大交路運(yùn)行。

1.2 符合以下條件之一的，需更換輪對(duì)

(1)換輪后試運(yùn)行動(dòng)車組：Δt≥65℃或t≥120℃。

(2)試運(yùn)行1萬(wàn)km內(nèi)動(dòng)車組：Δt≥50℃或t≥120℃。

(3)1萬(wàn)km≤動(dòng)車組試運(yùn)行里程<5萬(wàn)km：Δt≥47℃或t≥115℃。

換輪后運(yùn)行超過(guò)5萬(wàn)km動(dòng)車組：Δt≥45℃或t≥115℃或一周(7天)內(nèi)有2天同一軸箱軸承溫度t≥110℃。

這套規(guī)則是根據(jù)現(xiàn)行的軸溫預(yù)警規(guī)則，結(jié)合動(dòng)車使用經(jīng)驗(yàn)總結(jié)得出。但該規(guī)則只能在軸溫突破閾值(比如A級(jí)預(yù)警的120 ℃)后才能讓動(dòng)車使用方做出反應(yīng)，如果列車在行進(jìn)中軸溫再上升到達(dá)120 ℃, 列車必須限速運(yùn)行(250 km/h)，從而打亂列車運(yùn)行。通過(guò)研究軸溫與外溫、速度、輪對(duì)里程、軸承位置等特征之間的關(guān)系，分析產(chǎn)生軸溫過(guò)高的原因，從而實(shí)現(xiàn)對(duì)軸溫的智能預(yù)測(cè)，進(jìn)而提前指導(dǎo)維修。

2 數(shù)據(jù)分析

首先觀察單個(gè)車組在一天之內(nèi)的行駛和軸溫情況。圖1，圖2記錄了兩天之內(nèi)3517車組1車廂3軸左軸承溫度在兩天之內(nèi)與列車行駛速度和室外溫度的關(guān)系。我們能看到軸溫在一天之內(nèi)與列車行駛速度有關(guān)。車組于2016年8月11日早上出發(fā)，在13點(diǎn)左右到達(dá)G12車次的終點(diǎn)站北京南, 此時(shí)軸溫達(dá)到一天當(dāng)中的最高點(diǎn), 我們能夠看到列車保持高速時(shí)軸溫上升。在北京南休息了約50 min左右，下午車組又以G3車次返回上海虹橋。即便回程停站少，高速行駛時(shí)間比上午要長(zhǎng)，軸溫在下午并未超過(guò)中午的最高點(diǎn)。同樣的，在其他軸承上也發(fā)現(xiàn)了類似的規(guī)律： 即軸溫在上午的長(zhǎng)里程車次中升至一天最高點(diǎn)后，休息一段時(shí)間，下午同樣里程的車次很難使軸溫達(dá)到上午的最高溫。同時(shí)，從圖中可以看出軸溫與室外溫度也有一定的正相關(guān)，軸溫的升降和室外氣溫的升降幾乎是同步的。

圖1 3517車組1車廂3軸左軸溫與列車行駛速度

圖2 3517車組1車廂3軸左軸溫與室外溫度

2.1 軸溫與室外溫度

從圖3觀察可得，室外溫度與軸溫呈明顯的正相關(guān)，且軸溫突破最低閾值都發(fā)生在室外溫度高于25 ℃時(shí)。但從圖4可以看到，和行駛速度相比，室外溫度并不是軸溫的決定性因素。只要列車保持高速運(yùn)行，即便室外溫度降低至一天最低點(diǎn)，軸溫也會(huì)持續(xù)升高。

2.2 軸溫與走行速度

圖5反映了軸溫與列車運(yùn)行速度關(guān)系。3536車組在2016年8月19日運(yùn)行G230/G231車次，以300 km/h左右的最高速度在京滬高鐵上從上海行駛至濟(jì)南，然后通過(guò)膠濟(jì)線以200 km/h左右的最高速度行駛至終點(diǎn)站青島，下午再通過(guò)G236/G233原路返回。我們能看到軸承大幅升溫發(fā)生在速度300 km/h時(shí)。當(dāng)速度在200 km/h時(shí)，軸溫能升至一定程度并保持，但升溫幅度明顯降低。我們?cè)谄渌嚱M的軸承上也發(fā)現(xiàn)了類似情況。由此看來(lái)，軸溫升溫幅度與速度呈很明顯的正相關(guān)。在建立模型樣本時(shí)，我們有必要將列車運(yùn)行時(shí)間按不同的速度區(qū)間進(jìn)行劃分。

圖3 室外溫度與軸溫

圖4 室外溫度、速度與軸溫

圖5 軸溫與列車行駛速度的關(guān)系

根據(jù)數(shù)據(jù)分布情況我們創(chuàng)建了4個(gè)速度區(qū)間：320～370 km/h, 270～320 km/h, 220～270 km/h和170～220 km/h。圖6中從左至右的速度區(qū)間依次為：270～320 km/h, 220～270 km/h和170～220 km/h。從圖6能看到只有速度區(qū)間270～320 km/h所占比重與軸溫有持續(xù)正相關(guān)(分析樣本中沒(méi)有速度320 km/h以上的行駛記錄)。列車只有長(zhǎng)時(shí)間高速運(yùn)行時(shí)，軸溫才有可能不斷上升突破閾值。

2.3 軸溫與輪對(duì)里程

由趨勢(shì)圖分析知，不同的走行里程對(duì)于軸溫變化也是有影響的。圖7描述三級(jí)修后的輪對(duì)里程與每日最高軸溫和平均軸溫的關(guān)系。樣本中有兩個(gè)車組在2016年1月至今經(jīng)歷過(guò)三級(jí)修。由圖可以看出，在5萬(wàn)km以下時(shí)，每日平均軸溫和最高軸溫都還處在一個(gè)較為離散的狀態(tài)，突破最低閾值也發(fā)生在這個(gè)里程區(qū)間內(nèi)。在三級(jí)修后的10萬(wàn)km內(nèi)這兩個(gè)車組的軸溫與里程沒(méi)有明顯的正相關(guān)或負(fù)相關(guān)。

圖6 不同速度區(qū)間運(yùn)行時(shí)間占比與軸溫

圖8描述了四級(jí)修后的輪對(duì)里程與軸溫的關(guān)系。樣本中只有一個(gè)車組在2016年1月至今經(jīng)歷了四級(jí)修。在2萬(wàn)km內(nèi)，每日平均軸溫與最高軸溫都與里程呈正相關(guān)。但這之后至5萬(wàn)km時(shí)，軸溫緩慢下滑。軸溫與里程沒(méi)有呈線性關(guān)系。

圖7 輪對(duì)里程與軸溫(三級(jí)修后)

圖8 輪對(duì)里程與軸溫(四級(jí)修后)

圖9描述的是臨時(shí)換上的軸承溫度與里程的關(guān)系。小樣本中共有72個(gè)軸承在2016年1月后被臨時(shí)更換。換上的軸承在里程15萬(wàn)km以下時(shí)，每日平均軸溫和最高軸溫都相對(duì)較為離散。度過(guò)這個(gè)磨合期之后，每日軸溫在里程20萬(wàn)km時(shí)降至最低點(diǎn)，之后再緩慢回升。這些臨時(shí)換上的新輪突破最低軸溫閾值甚至都發(fā)生在里程不到1萬(wàn)km時(shí)。

圖9 輪對(duì)里程與軸溫(臨時(shí)換輪后)

綜上所述，輪對(duì)里程對(duì)軸溫是有影響的，但不是簡(jiǎn)單的正相關(guān)或負(fù)相關(guān)，與修程和走行里程的多少相關(guān)聯(lián)。

2.4 軸溫與軸承位置

由圖10可以看到第1和第16車廂出現(xiàn)的軸溫突破閾值的次數(shù)較多，第2、第6、第11和第12車廂的軸溫突破閾值的次數(shù)次之；軸位分布中1、2為動(dòng)軸，3、4為拖軸，其中第3軸發(fā)生的軸溫突破閾值次數(shù)最多，其他軸的軸溫突破閾值次數(shù)較為均勻。

圖10 突破溫度閾值的軸承在車廂和軸位的分布

為了進(jìn)一步檢驗(yàn)軸溫與軸位之間是否有顯著關(guān)系，對(duì)樣本數(shù)據(jù)進(jìn)行了T 檢驗(yàn)，結(jié)果證實(shí)動(dòng)軸的日最高和日平均軸溫都比拖軸高，其中平均每日最高軸溫66.5 ℃>64.25 ℃,每日平均軸溫50.1 ℃>49.1 ℃。但是，在0.05的顯著性水平下，T 檢驗(yàn)的P值大于0.05，說(shuō)明動(dòng)軸軸承的平均溫度與拖軸平均溫度并無(wú)顯著差異。因此，軸承位置中的軸位因素未加入模型中。

此外，由于文中樣本數(shù)據(jù)選取的是同一車廂的數(shù)據(jù)，軸承位置中的車廂變量為控制變量，因此模型中的輸入變量未考慮軸承位置。

3 RBF神經(jīng)網(wǎng)絡(luò)模型

RBF神經(jīng)網(wǎng)絡(luò)是一種三層前向網(wǎng)絡(luò)，通過(guò)輸入層空間到隱含層空間的非線性變換以及隱含層空間到輸出層空間的線性變換，實(shí)現(xiàn)輸入層空間到輸出層空間的映射。這兩個(gè)層間變換參數(shù)的學(xué)習(xí)可以分別進(jìn)行，使得RBF神經(jīng)網(wǎng)絡(luò)的學(xué)習(xí)速度較快且可避免局部極小問(wèn)題。

RBF(Radial Basis Function，徑向基函數(shù))是某種沿徑向?qū)ΨQ的標(biāo)量函數(shù)，通常定義為空間中點(diǎn) 到某一中心 之間歐氏距離的單調(diào)函數(shù)，最常用的徑向基函數(shù)是高斯函數(shù)，形式為式(1)

(1)

其中Cj為函數(shù)中心向量，σj為寬度向量。高斯函數(shù)的作用域表現(xiàn)出局部性，即當(dāng)x遠(yuǎn)離σj時(shí)函數(shù)取值較小。

如圖11所示，本次建模選擇的RBF神經(jīng)網(wǎng)絡(luò)的結(jié)構(gòu)從左至右分為3層，依次是輸入層、隱含層和輸出層：

圖11 RBF神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)圖

網(wǎng)絡(luò)的輸出如式(2)所示。

(2)

計(jì)算流程如圖12。

圖12 基于RBF神經(jīng)網(wǎng)絡(luò)的軸承溫度預(yù)測(cè)模型計(jì)算過(guò)程

4 模型結(jié)果

以預(yù)測(cè)2016年2月1日至20日的3537車01-1L軸箱軸承最高溫度為例，模型預(yù)測(cè)結(jié)果如圖13和表1。

圖13 2016年2月1日至20日3537車01—1L軸箱軸承最高溫度預(yù)測(cè)

日期最高軸溫預(yù)測(cè)值最高軸溫實(shí)際值絕對(duì)誤差相對(duì)誤差2月1日54.977 357.82.822 70.048 82月2日61.124 953.927.204 90.133 62月3日61.102 168.046.937 90.1022月4日56.224 551.74.524 50.087 52月5日60.5761.090.520.008 52月6日62.050 163.491.439 90.022 72月7日64.614 463.491.124 40.017 72月8日66.064 268.041.975 80.0292月9日72.838 472.590.248 40.003 42月10日69.028 469.180.151 60.002 22月11日46.514 246.430.084 20.001 82月12日64.178 163.490.688 10.010 82月13日60.593 458.941.653 40.028 12月14日65.184 565.360.175 50.002 72月15日58.97961.222.2410.036 62月16日62.454 762.350.104 70.001 72月17日67.081 666.90.181 60.002 72月18日69.610 269.180.430 20.006 22月19日49.791 749.840.048 30.0012月20日65.828 365.760.068 30.001

根據(jù)預(yù)測(cè)結(jié)果對(duì)軸箱軸承溫度進(jìn)行實(shí)時(shí)預(yù)警，預(yù)警的閾值如下：

一級(jí)預(yù)警：

(1)預(yù)測(cè)軸箱軸承溫度t≥110℃或115℃或≥120℃或t≥140℃；

(2)預(yù)測(cè)與同車同側(cè)其他3套軸承溫度平均值Δt≥50℃；

(3)預(yù)測(cè)與同車同側(cè)其他3套軸承溫度平均值Δt≥65℃。

二級(jí)預(yù)警：

(1)軸箱軸承溫度跳變：|軸箱軸承溫度(非空)-前1 min軸箱軸承(非空)|>12℃；

(2)10 min內(nèi)發(fā)生兩次軸箱軸承溫度跳變且速度>0，判定軸箱軸承溫度傳感器故障；

(3)預(yù)測(cè)軸箱軸承溫度t≥110℃或115℃，且據(jù)上次該故障發(fā)生小于等于7天(當(dāng)日多次發(fā)生不算)；

(4)預(yù)測(cè)與同車同側(cè)其他3套軸承溫度平均值Δt≥45℃。

三級(jí)預(yù)警：

預(yù)測(cè)軸承溫度≥80℃，且相對(duì)溫度*同側(cè)溫差≥2 600。其中，相對(duì)溫度=軸溫-環(huán)境溫度，同側(cè)溫差=軸溫-同側(cè)其他三軸平均值。

5 結(jié) 論

通過(guò)對(duì)CRH380動(dòng)車組軸箱軸承的故障數(shù)據(jù)進(jìn)行分析，總結(jié)了軸箱軸承故障發(fā)生的原因及占比，證明了溫度升高對(duì)于軸箱軸承的影響。為了研究軸箱軸承溫度升高的相關(guān)因素，分別比較了軸箱軸承溫度與外溫、運(yùn)行速度、輪對(duì)里程、軸承位置等特征之間的關(guān)系。最終選取了外界溫度、運(yùn)行速度和走行里程3個(gè)特征，利用RBF神經(jīng)網(wǎng)絡(luò)算法對(duì)軸箱軸承的溫度進(jìn)行建模預(yù)測(cè)。實(shí)驗(yàn)結(jié)果表明該模型具有較高的準(zhǔn)確率，可以為今后的動(dòng)車組維修保養(yǎng)工作提供參考性指導(dǎo)。另外，文中還存在一定的局限性，只考慮了軸溫過(guò)高(t≥120℃)情況下，軸溫與軸箱軸承故障的預(yù)警作用，對(duì)于軸溫小于120°的情況還有待繼續(xù)研究。