趙 凡,齊 琛,李偉斌,馬洪林,王躍軍
(中國空氣動力研究與發(fā)展中心,四川綿陽621000)
高速列車明線交會時相對速度很大,在交會的短時間內(nèi)兩車之間氣流產(chǎn)生劇烈波動,并產(chǎn)生較大的壓力波和側(cè)向力。影響明線交會壓力波的主要因素有行駛速度、線間距、列車頭型等,但車速一直是最主要的、也是調(diào)節(jié)余地最大的影響因素。
國內(nèi)外對列車明線交會進行了大量研究,田紅旗等[1]開展了靜止交會、等速和不等速交會等3種工況的研究,并得到了一系列回歸關(guān)系式。LIU等[2-3]研究了明線會車壓力波幅系數(shù)與鐵路線間距、交會速度間的關(guān)系,給出了任意車速會車時3者間統(tǒng)一的關(guān)系式。楊明智[4]、梁習(xí)峰[5]等做了列車等速交會的試驗研究,擬合出了交會壓力波幅值與速度的關(guān)系式。郗艷紅等[6]開展了高速列車明線會車壓力波波幅研究,修正了Steinheur的列車表面壓力波波幅經(jīng)驗公式。何德華等[7]利用列車空氣動力學(xué)模擬和輪軌動力學(xué)相結(jié)合的方法研究了動車組明線交會氣動力對動力學(xué)的影響。ZHANG等開展了高速列車等速明線交會時繞流流場的數(shù)值模擬研究[8]。Reinhardt等采用面元法模擬了無黏流下列車明線交會問題,給出了列車簡化模型的氣動力系數(shù)[9]。BI等研究了磁浮車不同線間距時明線交會的壓力波幅值,并依據(jù)計算結(jié)果給出了壓力波幅值允許的最小線間距[10]。
目前許多研究中雖計算了不同速度下壓力波的擬合關(guān)系式,但并未揭示不等速交會與等速交會壓力波之間的關(guān)系。且研究中多采用簡化模型,列車模型編組較少,不能反映列車真實外形對壓力波的影響??紤]到影響交會壓力波幅值的因素較多且關(guān)系復(fù)雜,研究時給定了線間距5 m并給定列車頭型。采用8車編組的真實復(fù)雜外形,重點研究速度對壓力波幅值的影響,以及等速交會、不等速交會壓力波幅值之間的關(guān)系。基于求解低速流動問題的SIMPLE算法開展列車交會非定常狀態(tài)的數(shù)值模擬,對于兩車的相對運動,采用滑移網(wǎng)格技術(shù)。首先開展了多組速度下等速交會研究,找出等速交會時相同監(jiān)測點在不同速度下壓力波幅值變化的關(guān)系;隨后對不等速交會工況和一車靜止工況開展仿真研究,分析比較與等速交會壓力波之間的關(guān)系,并得出相應(yīng)的關(guān)系式,為兩車交會時車速的調(diào)節(jié)控制提供參考。
高速列車的運行速度為250~350 km/h,兩車會車時相對速度達到了500~700 km/h,需求解三維可壓縮RANS方程。算法仍采用求解低速流動的SIMPLE算法,湍流模型采用標(biāo)準(zhǔn)k-ε模型。
雷諾平均的Navier-Stokes方程的通式可表示為[11]式(1):
對于連續(xù)方程為式(2):
對于x方向的動量方程為式(3):
y、z方向的動量方程具有類似的表示。
能量方程為式(4):
湍動能方程為式(5):
湍流耗散率方程為式(6):
理想氣體方程為式(7):
由于兩列車交會時有相對運動,需采用動網(wǎng)格技術(shù)。動網(wǎng)格技術(shù)主要包括滑移網(wǎng)格法、重疊網(wǎng)格法及網(wǎng)格變形技術(shù)。兩節(jié)列車的運動是平動,采用滑移網(wǎng)格法較為合適且計算量最小?;凭W(wǎng)格法是動網(wǎng)格方法的特殊情況,即網(wǎng)格塊和網(wǎng)格單元在運動的過程中形狀不變,在兩塊網(wǎng)格的交界面處滑移運動,變量在交界面處進行插值。對于滑移網(wǎng)格,守恒方程應(yīng)寫成式(8):
式中:V→g為網(wǎng)格速度。
對于滑移面的插值算法,有守恒型插值算法和非守恒型插值算法兩類。非守恒型插值算法直接用插值公式進行插值,計算量小,但無法保證交界面上通量守恒,對采用SIMPLE算法的低速求解器,非守恒型插值算法計算難以收斂。文中的計算采用守恒型插值算法,主要思路是將網(wǎng)格重疊部分的投影面積作為通量插值的權(quán)重。
如圖1所示,假設(shè)單元a,b,c,d所在面為主面,單元E所在面為副面,將主面的通量傳到副面。
圖1 滑移網(wǎng)格守恒型插值示意圖
式(9)中:Si為主面上的單元i與副面上單元E重疊的面積。式(10)中:ψi為主面上的單元i流場物理量,ψ*為副面上單元E的流場物理量。
為精細(xì)模擬交會時壓力波變化,計算時保留完整的列車轉(zhuǎn)向架,并保留路基、鐵軌等。列車頭型為ROCKET,列車編組采用實際運行中的8車編組。明線交會時兩車在相對運動,因此做出兩塊網(wǎng)格。圖2(a)給出了交會時兩列車的車體表面網(wǎng)格和地面網(wǎng)格,網(wǎng)格總量1900萬。圖2(b)、圖2(d)為頭車和尾車表面壓力的檢測點示意圖,圖2(c)給出了中間6節(jié)車上的6個測點。
圖2 列車網(wǎng)格和各節(jié)車表面測點位置示意
計算的坐標(biāo)系定義為x軸為列車運動方向,y軸為側(cè)力方向,z軸正向為升力方向。
開展了250,300,350 km/h 3種速度的等速明線交會計算,并給出了各節(jié)車表面壓力波測點的壓力幅值如圖3所示。從圖3可以看出,列車明線交會時會出現(xiàn)2個壓力波峰(谷),分別是對面車的車頭和車尾經(jīng)過本車監(jiān)測點時所產(chǎn)生。由于對面車車頭經(jīng)過時引起的壓力波要大于車尾經(jīng)過時引起的壓力波,因此文中只給出了車頭經(jīng)過時引起的壓力波幅值(即第1個壓力波幅值)。
圖3 等速明線交會頭車區(qū)域測點壓力波(v=250 km?h-1)
從表1給出的壓力波幅值可以看出,相同測點在不同速度下的壓力波幅值與速度平方成正比關(guān)系。以某速度下等速明線交會的計算結(jié)果作為基準(zhǔn)值,可以推導(dǎo)出其他速度下相同測點的壓力幅值為式(11):
表1 等速交會測點壓力波幅值 單位:Pa
式中:v0、v1為列車行駛速度;Δp0為以速度v0行駛時某測點的壓力波幅值;Δp1為以速度v1行駛時某測點的壓力波幅值。
列車實際運行中,經(jīng)常出現(xiàn)兩車交會時速度不相等的情況。此時由于兩車對空氣的擾動強度不同,上面式(11)的關(guān)系式已經(jīng)不再適用。為研究不等速明線交會時車體表面壓力波幅值特性規(guī)律,取一輛車(慢車)速度為100 km/h,另一列車(快車)速度分別為250,350,400 km/h。另外開展了其中一輛車靜止情況的研究,取一輛車速度為0(靜止),另一列車速度分別為250,350,500 km/h。由于慢車表面壓力波幅值要明顯大于快車,慢車部分測點壓力波幅值見表2。慢車速度為100 km/h、快車速度為250 km/h時,慢車頭車部分測點的壓力隨時間變化的曲線如圖4所示。
圖4 不等速明線交會慢車頭車區(qū)域測點壓力波及表面壓力云圖
從表2數(shù)據(jù)可以看出,每個工況中間車測點壓力波幅值比較接近,且各工況中間車測點之間的規(guī)律性較好。
表2 不等速交會慢車測點壓力波幅值 單位:Pa
根據(jù)計算數(shù)據(jù),列車交會的壓力波幅值可以由兩列車速度的關(guān)系式構(gòu)成,為式(12):
式中:v0為對面車速度,v1為自身車速;Δp1為監(jiān)測點的壓力波幅值;k為經(jīng)驗系數(shù),對于不同頭型或不同線間距,系數(shù)k的取值會有所差異。在文中的計算工況下,對于中間車測點,取系數(shù)k≈0.009 Pa·h2/km2。
式(12)適用于等速或不等速明線交會,可以較好的預(yù)測交會速度較小時的壓力波幅值,但在速度較大的交會工況,預(yù)測值略有偏差。
對計算的認(rèn)識可以追溯到古希臘時代對“數(shù)”的認(rèn)識。畢達哥拉斯學(xué)派對“數(shù)”的探索最為“癡迷”。他們認(rèn)為,世界的本源并非原初物質(zhì),而是其結(jié)構(gòu)形式,可以用“數(shù)”來衡量,用“數(shù)”來描述世界萬物,由此奠定了其“萬物皆數(shù)”的世界觀。但是,畢達哥拉斯學(xué)派所說的“數(shù)”是一種關(guān)于世界本體的哲學(xué)思考,沒有指向人類心理活動的探索。在那個時代,計算和心智彼此孤立、尚未結(jié)合。
為找到不等速交會與等速交會壓力波幅值的關(guān)系,通過式(13)定義一個無量綱化的不等速交會壓力波幅值C:
式中:Δp1為不等速交會時,以速度v1行駛的慢車中間車某測點的壓力波幅值,其對面快車的車速為v0;Δp0為以速度v0等速交會時車體表面中間車對應(yīng)測點的壓力波幅值。C為慢車監(jiān)測點上的無量綱壓力波幅值,代表慢車壓力波幅值和等速交會壓力波幅值之比(等速交會車速為不等速交會時快車車速)。
對于不等速交會時壓力波幅值C的取值,采用線性回歸法計算。通過線性回歸計算的不同速度下交會時的壓力波幅值C的值見表3。
表3 不同車速比下無量綱壓力波幅值C的線性回歸值
C值的計算為式(14):
式(14)給出了C與車速比m之間的擬合公式。無量綱壓力波幅值C與車速比m之間的擬合曲線如圖5所示,R2趨近于1表示該擬合關(guān)系下C與m強相關(guān)??梢钥闯鯟在車速比m小于0.5時其數(shù)值在0.78~0.8之間,變化范圍很小,說明此時慢車自身的車速還不足以對壓力波幅值產(chǎn)生較大影響。當(dāng)車速比m大于0.5后,壓力波幅值C的值明顯增大,直至速度比m為1時(即等速交會),C=1。
圖5 無量綱壓力波幅值C的值與速度比m之間的擬合關(guān)系
對于頭尾車測點壓力波幅值,當(dāng)速度比m為0~1時,Δp1/Δp0的范圍在0.6~1之間(鼻尖處測點除外)。這是因為頭車和尾車周圍流動變化劇烈,且交會時頭型上測點與對面車車頭的距離較遠,對面快車的影響權(quán)重下降,自身速度影響的權(quán)重上升,因此慢車速度變化對壓力波幅值影響更大。無論等速或不等速交會,最大壓力幅值監(jiān)測點均為尾車測點m14,該測點處于尾車車頭與車身的過渡位置,不等速明線交會時壓力波仍可用式(13)、式(14)進行計算。
(1)列車等速明線交會時,相同位置測點壓力波幅值與速度平方成正比。
(2)明線交會壓力波最大幅值通常出現(xiàn)在頭尾車側(cè)面位置,在頭型與頭尾車車廂的變曲面過渡區(qū);明線交會車體表面壓力波大小主要由對面車決定,對面車速度越大,壓力波幅值越大。
(3)列車不等速交會時,其壓力波幅值與等速交會的關(guān)系主要由慢車與快車的速度比m決定,具體關(guān)系式見式(12)、式(13)。在速度比m小于0.5時,不等速與等速交會壓力波幅值之比在0.78~0.8之間,說明此時慢車的車速對自身壓力波幅值影響較小。當(dāng)速度比m大于0.5后,壓力波幅值C的值明顯增大。
(4)根據(jù)結(jié)論(3),兩車交會時如速度過高需降低車速而保證行車安全性,慢車車速降到快車的0.5倍,即可達到較好的效果。
(5)對于同一車型在相同線間距下,可根據(jù)某一速度下等速交會的壓力波幅值結(jié)果,用式(11)估算其他任意速度等速交會時壓力波幅值。對于不等速交會,中間車測點的壓力波幅值可以用式(11)、式(13)、式(14)估算。該系列公式具有較好的預(yù)測精度,且可以節(jié)省大量計算量和計算時間。