吳 寧,董孝卿,林鳳濤,文 彬,王悅明
(1 中國鐵道科學(xué)研究院 機車車輛研究所,北京100081;2 華東交通大學(xué) 載運工具與裝備省部共建教育部重點實驗室,江西南昌330013)
輪軌幾何接觸是解釋鐵道車輛動力學(xué)性能的基礎(chǔ),等效錐度是輪軌幾何接觸中的重要參數(shù)。當鐵道車輛運行在直線線路或大半徑曲線線路上時,等效錐度決定著輪軌之間的匹配程度[1-2]。
對我國多條線路上的動車組振動狀態(tài)及車輪磨耗跟蹤研究的結(jié)果表明,輪軌匹配等效錐度的大小與動車組運行過程中的動態(tài)響應(yīng)密切相關(guān)。如等效錐度過小導(dǎo)致動車組運用過程中出現(xiàn)晃車現(xiàn)象[3-4],等效錐度過大引起動車組車輛構(gòu)架橫向振動報警[5]。因此,定期獲取運用動車組車輪廓形并準確計算等效錐度是很必要的。
在關(guān)于輪軌關(guān)系的研究領(lǐng)域內(nèi),等效錐度是一個重要參數(shù)。等效錐度與鐵道車輛的動力學(xué)性能有著密切的聯(lián)系,歐洲標準UIC518[6]將等效錐度作為車輛型式試驗的關(guān)鍵參數(shù)并規(guī)定了試驗過程中等效錐度的范圍。國內(nèi)鐵道車輛領(lǐng)域內(nèi)很多研究工作探討了等效錐度和車輛動力學(xué)性能的關(guān)系[7-11]。文獻[12]比較了幾種等效錐度計算方法的差別。
對于確定的輪軌廓形和參數(shù),等效錐度的計算結(jié)果應(yīng)該是確定的。等效錐度算法并不局限于一種,且計算結(jié)果受輪軌廓形平滑插值方法影響較大。因此需要有相關(guān)驗證標準對等效錐度算法進行驗證,一種等效錐度算法只有在完全通過相關(guān)標準驗證后其計算結(jié)果才能被認可。歐洲標準EN 15302用于等效錐度算法的驗證[1],UIC519給出了等效錐度的計算方法[2]。
錐形踏面輪對在線路上具有相對固定的正弦運動軌跡,Klingel理論指出這一運動軌跡的波長取決于輪對踏面的錐度角以及左右輪軌接觸斑之間的距離,即Klingel公式:
式中λ為輪對運動軌跡的波長,m;e為左右輪軌接觸斑之間的距離,即為跨距,m;r0為輪對處在對中位置時的車輪滾動圓半徑,m;γ為車輪踏面的錐度角,弧度。
實際運用中的車輪踏面具有變化的錐度角,可以采用對左右車輪滾動圓半徑差(取決于輪對橫移量)進行積分的方法獲得輪對的運動波長,再將該波長與Klingel理論中的相應(yīng)波長對照,得出輪對等效錐度。即:
式中tanγe為等效錐度。
等效錐度定義為給定輪對運動波長相同的錐形車輪輪對的車輪錐度角正切值。如圖1所示。
圖1 等效錐度定義
等效錐度的實質(zhì)是非錐形踏面輪對在每一橫移幅值下均等效于一個錐形踏面輪對。其等效的依據(jù)是兩者的運動軌跡波長相等。
等效錐度計算的流程見圖2。主要分為以下幾個步驟:(1)對輪軌廓形進行平滑和插值處理;(2)根據(jù)輪軌廓形和位置參數(shù)計算得到輪軌接觸點;(3)根據(jù)輪軌接觸點得到與輪對橫移量對應(yīng)的左右車輪滾動圓半徑差;(4)基于左右車輪滾動圓半徑差計算結(jié)果,計算等效錐度。
圖2 等效錐度計算流程
以上4個步驟中,每一步均可采用不同的計算方法。步驟(1)可選擇采用通用的數(shù)學(xué)平滑插值處理方法,步驟(2)和步驟(3)采用的計算方法均為一般意義上的解析幾何方法。以下簡要由滾動圓半徑差計算等效錐度步驟(4)的2種方法(圖3)。
方法1:半徑差函數(shù)線性回歸算法。根據(jù)滾動圓半徑差積分確定輪對橫移幅值和橫移量的對應(yīng)關(guān)系,再由橫移幅值對應(yīng)的橫移量確定滾動圓半徑差的線性回歸范圍,對該范圍的滾動圓半徑差進行線性擬合,擬合后的斜率之半即為對應(yīng)輪對橫移幅值下的等效錐度。
方法2:非線性微分方程積分算法。根據(jù)滾動圓半徑差積分確定輪對橫移幅值和橫移量的對應(yīng)關(guān)系,再由邊界條件(即橫移幅值對應(yīng)的最大橫移量和最小橫移量對應(yīng)的搖頭角Ψ均為0)計算求解搖頭角Ψ時的常數(shù)C,得到常數(shù)C后積分確定搖頭角Ψ與橫移量y的函數(shù),之后再進行積分,確定出輪對橫移幅值對應(yīng)的輪對運動波長λ,最后由公式(2)得到等效錐度tanγ。
圖3 由滾動圓半徑差計算等效錐度的2種算法
根據(jù)上面介紹的等效錐度計算方法,采用方法1實現(xiàn)了等效錐度計算的算法編程,以適應(yīng)批量計算實測輪軌廓形等效錐度的需要。
(1)輪軌廓形的插值平滑處理
輪軌廓形的插值平滑處理主要根據(jù)輪軌廓形坐標數(shù)據(jù)及輪軌幾何參數(shù)(軌距,輪對內(nèi)側(cè)距,軌底坡等)實現(xiàn)輪軌廓形坐標的插值、平滑及輪軌位置定位。采用兩次插值處理,一次平滑處理實現(xiàn)輪軌廓形坐標插值平滑處理:先對初始輪軌廓形坐標值進行插值處理,接著對插值處理后的坐標進行平滑處理,最后再對平滑處理后的坐標進行插值。其中平滑處理采用濾波的方法實現(xiàn)。插值平滑處理完成后根據(jù)最終的輪軌廓形坐標及輪軌幾何參數(shù)輸出用于計算輪軌接觸點的左右側(cè)輪軌廓形坐標。
(2)計算輪軌接觸點
采用左右側(cè)輪軌垂向距離最小極值的方法計算輪軌接觸點,具體的計算流程見圖4。輪對橫移量一般從-20mm計算到20mm。左右側(cè)輪軌最小垂向距離之差絕對值的控制變量ε是一個很小的正數(shù),其值越小,輪軌接觸點的計算精度越高。但若ε過小,則可能得到錯誤的輪對側(cè)滾角,從而得到錯誤的輪軌接觸點。ε的選擇與輪軌插值點橫向坐標間距有關(guān),如輪軌插值點橫向坐標間距取0.1mm時,ε可取0.01mm。圖5為計算得到的輪軌接觸點示意圖。
(3)計算左右車輪滾動圓半徑差
在得到輪軌接觸點后,根據(jù)輪對橫移量對應(yīng)的左右車輪上接觸點可直接得到左右車輪滾動圓半徑差。
(4)根據(jù)左右車輪滾動圓半徑差計算等效錐度
等效錐度的實質(zhì)是非錐形踏面輪對在每一橫移幅值下均等效于一個錐形踏面輪對,等效錐度是輪對橫移幅值的函數(shù),因此首先要確定輪對橫移量和輪對橫移幅值的對應(yīng)關(guān)系。
確定半徑差為0時對應(yīng)的橫移量yem,該橫移量對應(yīng)的輪對位置即為輪對相對于軌道的平衡位置。對半徑差進行積分,積分S(y)=-∫Δr·dy分別從yem向橫移量正方向y+和負方向y-各進行一次,確定出在輪對橫移振幅ym下輪對的最小橫移量yemin和最大橫移量yemax。對最小橫移量yemin和最大橫移量yemax之間的半徑差進行直線擬合,擬合得到的斜率之半即為輪對橫移振幅ym下的等效錐度。
圖4 輪軌接觸點及計算流程
圖5 輪軌接觸點示意圖
如果Δr=f(y)特征在區(qū)間[yemin,yemax]內(nèi)關(guān)于原點對稱性較差,應(yīng)對每部分(正和負)分別進行回歸得到兩個等效錐度值tanγe,p和tanγe,n,等效錐度可通過下式得到:
由于計算等效錐度可采用不同的計算方法,且計算等效錐度的每一步均可選取不同的數(shù)學(xué)算法。對于確定的輪軌廓形及輪軌參數(shù),等效錐度計算結(jié)果應(yīng)是確定的、穩(wěn)定的。因此必須采用相關(guān)標準對整個等效錐度算法進行驗證。一個具體的等效錐度算法只有完全通過相關(guān)標準驗證后才能用于實測輪軌廓形等效錐度的計算。
歐洲標準EN 15302:2008規(guī)定了驗證等效錐度算法的過程、驗證采用的基準輪軌廓形及采用基準輪軌廓形計算得到的等效錐度基準結(jié)果和容許誤差。采用EN 15302:2008對整個等效錐度算法程序進行驗證。整個驗證過程共分為3步。
(1)采用EN 15302:2008附錄中給出的基準輪軌廓形數(shù)據(jù)及匹配參數(shù)作為輸入,考核被驗證的算法程序的計算結(jié)果是否在規(guī)定的容許誤差界定范圍內(nèi)。其目的是考核計算程序能否準確地計算出等效錐度,即驗證計算程序的準確性。圖6給出了一個驗證結(jié)果示例,本例中左側(cè)車輪廓形為標準附錄中的車輪B廓形,右側(cè)車輪為標準附錄中的車輪A廓形,軌道廓形為標準附錄中的軌道A廓形。圖中同時給出了ADAMS軟件輪軌接觸計算模塊的計算結(jié)果。
圖6 等效錐度算法程序驗證第1步圖例
(2)將EN 15302:2008附錄中定義的隨機誤差添加到基準輪軌廓形數(shù)據(jù)中,考核被驗證算法程序的計算結(jié)果是否在規(guī)定的容許誤差界定范圍內(nèi)。其目的是當考核輪軌廓形有一定誤差時計算程序的計算結(jié)果是否穩(wěn)定,即驗證計算程序的穩(wěn)定性。圖7給出了一個驗證結(jié)果示例,本例中車輪廓形為標準附錄中的車輪A廓形,軌道廓形為標準附錄中的軌道A廓形。
圖7 等效錐度算法程序驗證第2步圖例
(3)將EN 15302:2008附錄中定義的柵格誤差添加到基準輪軌廓形數(shù)據(jù)中,考核被驗證的算法程序計算結(jié)果是否在規(guī)定的容許誤差界定范圍內(nèi)。其目的是考核輪軌廓形數(shù)字化過程中引入柵格誤差后計算程序的計算結(jié)果是否穩(wěn)定,即驗證計算程序能否用于實測輪軌廓形等效錐度的計算。圖7給出了一個驗證結(jié)果示例,本例中車輪廓形為標準附錄中的車輪A廓形,軌道廓形為標準附錄中的軌道A廓形。
圖8 等效錐度算法程序驗證第3步圖例
詳細地介紹了等效錐度的定義、計算方法和算法驗證過程。以期對等效錐度算法編程實現(xiàn)及算法程序驗證有很好的參考作用。
[1]EN15302:2008,Railway applications-Method for determining the equivalent conicity[S].
[2]UIC 519;2004;Method for determining the equivalent conicity[S].
[3]董孝卿,王悅明,倪純雙,吳 寧,等.京津城際鐵路動車組車輪運用情況深入研究,TY字第3041號[R].北京:中國鐵道科學(xué)研究院機車車輛研究所,2010.
[4]董孝卿,王悅明,王林棟,倪純雙,等.京津城際鐵路CRH3C型動車組車體晃動專題研究報告,TY字 第3056號[R].北京:中國鐵道科學(xué)研究院機車車輛研究所,2011.
[5]康 熊,董孝卿,王悅明,文 彬.車輪外形磨耗及動車組振動狀態(tài)跟蹤研究,TY字第3217號[R]北京:中國鐵道科學(xué)研究院機車車輛研究所,2011.
[6]UIC518;2005;Test and approval of railway vehicles from the point of view of their dynamic behaviour-Safety-Track fatigue-Ride quality[S].
[7]李 霞,溫澤峰,張劍,等.軌底坡對輪軌滾動接觸行為的影響[J].機械強度,2009,(03):130-135.
[8]羅世輝.軌距對機車車輛穩(wěn)定性影響的研究[J].中國鐵道科學(xué),2010,(2),56-60.
[9]張 劍,金學(xué)松,等.基于CRH5型高速動車組車輛的輪對動態(tài)特性與等效錐度關(guān)系初探[J].鐵道學(xué)報,2010,32(3):20-27.
[10]樸明偉,樊令舉,等.基于輪軌匹配的車輛橫向穩(wěn)定性分析[J].機械工程學(xué)報,2008,44(3):21-28.
[11]陳厚嫦,黃體忠,等.輪對內(nèi)側(cè)距對機車車輛動力學(xué)性能影響的試驗研究[J].中國鐵道科學(xué),2006,27(5):99-103.
[12]任利惠,謝 綱,沈 鋼,等.踏面等效斜率計算方法的比較[J].城市軌道交通研究,2011,(6):8-13.