支洋

（中國國家鐵路集團有限公司 工電部，100844）

0 引言

軌道結(jié)構(gòu)是整個高速鐵路基礎(chǔ)設(shè)施的重要組成部分，車輛運行安全性和舒適性要求線路具備高平順性和高可靠性。目前，我國已投入運營高鐵線路總里程數(shù)將近4萬km，隨著運營時間的增加，軌道在車輛動載荷反復(fù)作用下結(jié)構(gòu)狀態(tài)將產(chǎn)生變化，微小變形積累到一定程度會形成軌道不平順，如不及時整治將演變?yōu)檐壍啦『Α２ㄩL1 m以下的短波病害，如焊接接頭不良、擦傷、鋼軌波磨等，易引起軌道-車輛系統(tǒng)的劇烈振動，縮短車輛及軌道部件的使用壽命，增加鐵路養(yǎng)護維修費用，是危害行車安全的潛在隱患。

軌道短波病害主要分為2類，一類是單點沖擊短波病害，如焊接接頭不良、凹坑和軌面擦傷等；另一類是連續(xù)周期性短波病害，如打磨痕跡、鋼軌波磨等。軌道幾何檢測系統(tǒng)的有效檢測波長范圍通常在2～120 m，難以評判波長1 m以下的短波病害對高速鐵路軌道-車輛系統(tǒng)運行的安全性和舒適性的影響，需要研究其他診斷和評判方法。

軸箱直接與輪對相連，利用軸箱加速度評判軌道短波狀態(tài)，不但可直接反映短波不平順激發(fā)的外界力對車輛動力學(xué)的影響，同時可以突破軌道幾何檢測系統(tǒng)波長限制和人工地面檢測效率低下的瓶頸，為保障行車安全、實現(xiàn)高速鐵路的高效維修與管理提供強有力的理論與技術(shù)支撐。

1 短波狀態(tài)動態(tài)檢測技術(shù)框架

利用高速綜合檢測列車實測的軸箱加速度動態(tài)診斷軌道短波病害，技術(shù)框架見圖1，主要包括短波病害動力學(xué)仿真、特征提取和動態(tài)診斷方法。在短波病害動力學(xué)仿真方面，早期采用多體動力學(xué)方法建模，模擬短波病害幾何參數(shù)對車輛動態(tài)響應(yīng)的影響；后來隨著計算技術(shù)特別是并行技術(shù)的發(fā)展，能更真實模擬短波病害形狀和材料非線性的有限元仿真成了首選。利用時域、頻域、時頻方法從軸箱加速度中挖掘能表征短波病害的時域特征，確定帶通濾波的頻率范圍。在短波病害動力學(xué)仿真和特征提取的基礎(chǔ)上，利用軸箱加速度的特征指標(biāo)如軌道沖擊指數(shù)，動態(tài)評判軌道短波的平順狀態(tài)。

圖1 高速鐵路軌道短波狀態(tài)動態(tài)檢測技術(shù)框架

2 短波狀態(tài)仿真分析

早期Jenkins等［1］利用現(xiàn)場實測的輪軌力統(tǒng)計分析鋼軌焊接接頭不良對軌道-車輛系統(tǒng)動力學(xué)性能的影響，發(fā)現(xiàn)車輪經(jīng)過鋼軌低接頭過程中會出現(xiàn)2個輪軌力峰值，分別為P1和P2力。P1力持續(xù)時間較短，主要受鋼軌等效質(zhì)量、列車運行速度以及車輛系統(tǒng)簧上質(zhì)量影響，為高頻瞬態(tài)沖擊力；P2力持續(xù)時間較長，主要與軌道系統(tǒng)等效剛度相關(guān)，為中低頻響應(yīng)力，對軌道結(jié)構(gòu)的破壞較為明顯。

隨著計算機技術(shù)的快速發(fā)展，仿真分析方法包括多體系統(tǒng)動力學(xué)和有限元方法逐漸在軌道-車輛系統(tǒng)動力學(xué)分析中占重要地位。SUN等［2］利用自主開發(fā)的軌道-車輛系統(tǒng)耦合仿真分析軟件，探索了鋼軌焊接接頭處列車運行速度與輪軌作用力的關(guān)系。

與多體系統(tǒng)動力學(xué)相比較，有限元方法能更準(zhǔn)確模擬軌道短波病害、輪軌材料和接觸非線性對軌道-車輛系統(tǒng)動態(tài)特性的影響。Dong［3］利用垂向軌道-車輛系統(tǒng)有限元模型，得出了車輛通過鋼軌低接頭時P1力和P2力的波動規(guī)律。Molodova等［4］利用ANSYS建立三維輪軌有限元模型，探明了鋼軌凹痕與軸箱加速度的關(guān)聯(lián)關(guān)系。于淼等［5］應(yīng)用ABAQUS軟件建立三維高速輪軌非穩(wěn)態(tài)滾動接觸有限元模型，初步探明了復(fù)合短波病害的形成機理。

利用有限元方法模擬短波病害對軌道-車輛系統(tǒng)影響的基本思路是，先建立軌道和輪對的幾何模型，通過修改鋼軌表面節(jié)點坐標(biāo)對軌道短波病害建模，接下來計算軌道-車輛系統(tǒng)的靜力平衡，最后采用顯式有限元方法模擬車輛通過短波病害區(qū)域的動態(tài)響應(yīng)。

3 短波病害動態(tài)響應(yīng)特征提取方法

從軸箱加速度中提取軌道短波病害特征的方法包括時域、頻域和時頻分析等方法。信號的時域特征指標(biāo)主要包括峰值、均值、方差、標(biāo)準(zhǔn)差等；頻域分析主要采用傅里葉變換揭示頻率與能量的對應(yīng)關(guān)系。

由于輪軌作用是強非線性接觸問題，用來診斷軌道短波病害的軸箱加速度信號呈現(xiàn)高度非線性和非平穩(wěn)特性。單純的時域或頻域方法適合于平穩(wěn)信號，對于非平穩(wěn)信號，時頻分析方法聯(lián)合時間-頻率表示一維時間信號的特征，具有獨特的優(yōu)勢。經(jīng)典時頻分析包括STFT、Wigner-Ville分布、Chirplet變換、EEMD等。2010年Daubechies等［6］提出了同步壓縮小波變換（SST），根據(jù)瞬時頻率對時頻分布進行重排。徐曉迪等［7］對SSTFT進行了改進，優(yōu)化了脊線提取方式，將零碎曲線連接成完整的頻率曲線。

鋼軌擦傷區(qū)段在列車經(jīng)過時接觸應(yīng)力增加導(dǎo)致剪切應(yīng)力也增大，嚴(yán)重時會增加鋼軌表面疲勞，并引起鋼軌掉塊等病害。利用改進的SSTFT方法對不同深度的鋼軌擦傷區(qū)段的軸箱加速度進行時頻分析，結(jié)果見圖2。圖中能量較大的點為焊接接頭對應(yīng)位置，鋼軌擦傷區(qū)段的軸箱加速度在頻域內(nèi)呈現(xiàn)全頻段的白噪聲特性，這是由于列車經(jīng)過擦傷區(qū)域時輪軌接觸點恰好在引起掉塊的地方，從而產(chǎn)生了沖擊響應(yīng)。同時，在100～400 Hz之間能量比較集中，這是由于擦傷區(qū)域長度以及軌面粗糙度較大引起的中高頻響應(yīng)。對比圖2（a）、圖2（b）可知，擦傷深度為2.5 mm時產(chǎn)生的沖激響應(yīng)能量要比擦傷深度為1.5 mm時的大，全頻段的白噪聲特性也比較明顯。對比圖2（c）和圖2（d）可知當(dāng)綜合檢測列車以更高的速度經(jīng)過鋼軌擦傷區(qū)段時，速度越大，軸箱加速度的振動能量越大。

圖2 不同深度的鋼軌擦傷區(qū)段軸箱加速度時頻分布

4 短波病害動態(tài)診斷方法

如果將輪對近似地看成是剛性結(jié)構(gòu)，則軸箱的振動加速度可直接反映由軌道短波不平順?biāo)鸬耐饨缂_力對車輛動力學(xué)的影響。因此，許多學(xué)者采用軸箱加速度輔助評判軌道短波不平順對車輛動力學(xué)性能的影響。

2004年，日本學(xué)者Sunaga等［8］提出利用軸箱加速度的有效值診斷新干線上的軌道短波病害，發(fā)現(xiàn)軌枕松動、軌面粗糙、鋼軌波磨等短波病害對應(yīng)不同頻段的軸箱加速度，可通過對軸箱加速度設(shè)置不同的濾波范圍進行診斷。

2016年，Molodova等［9］利用軸箱加速度對焊接接頭的健康狀態(tài)進行診斷，并給出了定量的評判指標(biāo)，即尺度平均小波功率譜指標(biāo)（scale-averaged wavelet power，SAWP）。王林棟［10］在借鑒國內(nèi)外算法的基礎(chǔ)上提出了基于振動響應(yīng)的鋼軌波磨快速檢測方法；當(dāng)線路里程長且采樣頻率高于1 000 Hz時，常規(guī)的診斷方法會碰到計算量大滿足不了在線應(yīng)用要求的問題。劉金朝等［11-12］通過挖掘軸箱加速度能量指標(biāo)計算過程中的遞推映射關(guān)系，提出了快速診斷軌道短波病害的軌道沖擊指數(shù)方法，詳細(xì)過程可描述如下：

（1）根據(jù)軌道短波病害對應(yīng)的軸箱加速度的頻率分布特性，取通濾波頻率一般為 ［10，500］Hz，對實測的軸箱加速度進行帶通濾波；

（2）利用能量指標(biāo)的遞推映射關(guān)系，快速計算帶通濾波后的軸箱加速度的移動有效值S及其平均值；

（3）計算軸箱加速度的軌道沖擊指數(shù)TII，定義為軸箱加速度的移動有效值S與其平均值的比值；

（4）對軌道進行單元劃分，單元長度一般取為50 m；

（5）提取各單元TII的最大值，記為TIIe，同時記錄對應(yīng)的里程信息；

（6）結(jié)合閾值評判軌道短波狀態(tài)，即若TIIe大于閾值，則診斷軌道短波狀態(tài)不良。

利用軌道沖擊指數(shù)方法對某線路K1+523附近的軸箱加速度進行分析，得到軌道沖擊指數(shù)（見圖3）。軌道沖擊指數(shù)超過閾值6，診斷該處焊接接頭狀態(tài)不良?，F(xiàn)場復(fù)核發(fā)現(xiàn)鋼軌表面魚鱗傷嚴(yán)重、同時存在剝落掉塊（見圖4），驗證了軌道沖擊指數(shù)方法能有效診斷鋼軌短波病害。

圖3 某線路K1+523附近的軌道沖擊指數(shù)

圖4 鋼軌剝落掉塊

5 結(jié)論與展望

為了從演變機理上找到抑制病害發(fā)展的養(yǎng)護措施，軌道短波病害對軌道-車輛系統(tǒng)動力學(xué)性能影響的研究一般采用多體系統(tǒng)動力學(xué)和有限元方法，相比較而言，有限元方法由于能細(xì)化考慮材料、接觸和高速轉(zhuǎn)動的非線性，能更精確描述短波病害與車輛響應(yīng)的關(guān)聯(lián)關(guān)系。在挖掘提取軌道短波特性的時域、頻域、時頻分析方法中，時頻方法因為能同時提取短波病害發(fā)生的位置和頻率特性而顯得優(yōu)勢明顯。短波病害的診斷方法的發(fā)展路線是從時域、頻域、時頻，然后再回到時域即軌道沖擊指數(shù)方法。

雖然國內(nèi)外對軌道短波病害的仿真分析和診斷方法已開展了大量研究，但仍有一些深層次的問題有待解決。如何結(jié)合數(shù)值仿真方法和信號挖掘技術(shù)，找出高速鐵路鋼軌波磨的產(chǎn)生機理、發(fā)展規(guī)律從而找到針對性的抑制措施是有待攻克的難題；輪軌接觸共振是公認(rèn)的鋼軌波磨形成的理論，但共振發(fā)生的隨機性還有待揭示。此外，復(fù)合型短波病害的診斷、識別、評價等也是有待深入研究的問題，如何從車輛響應(yīng)特性上將單點沖擊短波病害和周期性短波病害分離是其識別的難點。