高 暢,張繼業(yè),李 田,孫 瑤

(西南交通大學牽引動力國家重點實驗室,成都 610031)

隨著高速列車的不斷提速，列車空氣動力學問題日益凸顯[1-2]，速度的增加使列車阻力大幅增加，同時，列車在側(cè)風、隧道內(nèi)等環(huán)境下運行，會引發(fā)出危及列車運行安全的問題[3-7]。但是通過利用列車高速運行時的空氣阻力轉(zhuǎn)化為制動力，服務(wù)于列車制動[8]。有效可靠的列車制動在列車安全運行中扮演著重要的角色。風阻制動作為一種非黏著制動技術(shù)，具有減少輪軌損耗、結(jié)構(gòu)簡單、綠色環(huán)保等特點，風阻制動是通過增加空氣阻力來提供制動力，列車速度越高，制動效果越明顯。日本最早開展了高速列車風阻板的研究，通過實驗開展了山梨超導磁懸浮試驗線風組制動的研究，風阻制動很好地彌補了高速制動工況下黏著制動的缺陷[9-10]。隨著風阻制動的研究深入，日本普遍認為速度350 km/h以上的高速列車采用風阻制動作為緊急情況的輔助制動方式，于是在Fastech360z型列車設(shè)有風阻制動裝置。韓國也積極開展了風阻制動的研究，采用四節(jié)編組列車，開展了頭車制動板的角度對制動性能影響的研究[11]。我國在風阻制動的研究相對日本起步較晚，隨著我國高速列車速度的快速攀升，高速列車的風阻制動的需求日益凸顯。國內(nèi)學者的研究也在逐步深入，高立強等研究了制動板之間的氣動干擾效應(yīng)。制動板之間的軸向距離越大干擾越小，制動效果越好[12-13]，同時，采用數(shù)值仿真的方法對制動板的板型及其對流場結(jié)構(gòu)和氣動特性的影響進行了探究[14]。田春等采用固定位置的八套風阻制動板，研究了列車制動板縱向布置規(guī)律[15-17]。

目前國內(nèi)外學者沒有對制動板的具體最佳位置進行探討。所謂的最佳位置，相同大小的制動板相對于其他位置產(chǎn)生的空氣阻力最大。高速列車周圍的流場結(jié)構(gòu)復雜，制動板的位置不同產(chǎn)生的空氣阻力也大不相同。本文將對制動板的位置進行探究，優(yōu)化風阻制動效果，為制動板的分布設(shè)計提供參考依據(jù)。

1 控制方程數(shù)值仿真模型

1.1 流體控制方程

當高速列車高速運行時，周圍流場考慮為三維黏性定常的湍流流場。當列車運行速度為400 km/h，風速為111.11 m/s>98.40 m/s，對應(yīng)的馬赫數(shù)大于0.3，列車附近的流場可近似處理為可壓縮流場，湍流模型采用k-ε標準模型，指出其控制方程形式為[18-19]

(1)

式中，t為時間；ρ為空氣密度；Γ為廣義擴散系數(shù)；u為流場速度矢量，u=[u v ω]；ut為網(wǎng)格移動速度矢量，ut=[utvω]；φ為流場通量；S為廣義源項。

1.2 幾何模型

風阻制動板采用整體式制動板，采用同空調(diào)機組的安裝結(jié)構(gòu)密封方式，因此制動板周圍車體表面可以考慮為光滑表面。制動板沿著列車縱向中間對稱面對稱分布。制動板的開合角度為75°最優(yōu)[20]，本文采用風阻板與列車運行方向的角度為75°，考慮到列車限界以及風組制動裝置的開合結(jié)構(gòu)，制動板提供1 280×273.35 mm2的迎風面積，板的厚度為25 mm。風阻制動板的外形尺寸見圖1。

列車模型采用國內(nèi)某型號高速動車組全尺寸模型為研究對象，采用頭車+兩節(jié)中間車+尾車的四車編組方式，忽略車體外部某些復雜的細部結(jié)構(gòu)，如門把手，風擋，受電弓等部件。列車模型如圖2所示。列車基本參數(shù)見表1。

圖1 風阻制動板模型

圖2 列車及擋風板幾何模型

表1 列車基本參數(shù)

計算區(qū)域見圖3，長度、寬度和高度分別為424，80 m和40 m。列車沿著x正方向運行，列車正前方速度入口邊界距頭車鼻尖處1倍車長，長度106 m，列車正后方壓力出口邊界距尾車鼻尖處212 m，頂面邊界距地面40 m，列車左右側(cè)邊界距離車身相同，距離為40 m，列車底板與地面之間的距離為0.376 m。

圖3 計算區(qū)域(單位：m)

列車以400 km/h的速度沿著x正方向運行。列車表面設(shè)置為無滑移壁面邊界條件，為了模擬地面效應(yīng)，地面設(shè)置為滑移壁面，其滑移速度為列車運行速度。列車運行的正前方為速度入口，速度入口風速為列車運行速度，風速為111.11 m/s，列車運行正后方為壓力出口，頂面與側(cè)面均設(shè)為對稱邊界。

本文主要是研究前排風阻板的位置對列車風阻制動的影響，需要對風阻板的位置進行有效標記。首排風阻板即沿著列車運行方向的第一塊制動板，位于頭車。頭車分為前端流線型車身和后端平順型車身?？紤]到實際安裝所需空間，風阻板制動板安裝于頭車的平順型車身。為便于描述風阻板位置，流線型尾端，即流線型與平順型車身交界處，記為D0，距離流線型尾端為x的位置記為Dx，比如在平順型車身上的制動板，距離流線型尾端6 m處的位置記為D6，具體見圖4。

圖4 制動板位置示意

1.3 數(shù)學模型

采用基于有限體積法商業(yè)軟件FLUENT對列車風阻制動板的氣動特性進行分析，列車空氣動力學的數(shù)值模擬采用三維、定常、不可壓縮Navier-Stoke方程和k-ε兩方程湍流模型，CFD模擬分析所采用的高階有限差分格式見表2。

表2 CFD模擬所采用的高階差分格式

2 網(wǎng)格獨立性驗證

本文采用ICEM-CFD軟件對計算區(qū)域進行網(wǎng)格劃分，網(wǎng)格形式采用四面體網(wǎng)格。車體網(wǎng)格最小尺寸為40 mm。為保證計算精度，在制動板和列車周圍處設(shè)置加密區(qū)，同時在列車表面劃分為邊界層。為避免網(wǎng)格差異對計算結(jié)果的影響，制動板位于D0處，劃分了5套不同尺度的網(wǎng)格并分別進行氣動性能計算。比較網(wǎng)格變化對列車頭車阻力的影響，對比結(jié)果見表3。

由表3可以看出：第1套網(wǎng)格的頭車阻力較第2套大92 N，較第3套大71 N，同時在局部加密網(wǎng)格后(第4套網(wǎng)格和第5套網(wǎng)格)，頭車阻力值變化很小，可見第3套網(wǎng)格滿足網(wǎng)格獨立性要求。因此，采用第3套網(wǎng)格進行數(shù)值模擬，邊界層第1層厚度0.1 mm，邊界層層數(shù)為10層，增長率為1.2。網(wǎng)格數(shù)量約3 235萬。圖5為局部網(wǎng)格示意。

表3 網(wǎng)格獨立性檢驗對比

圖5 局部網(wǎng)格示意

3 前排風阻制動板對列車氣動性能的影響

3.1 風阻制動板制動力的分布規(guī)律

為了了解風阻制動板制動力的總體分布規(guī)律，在列車上設(shè)置了8塊風阻制動板，首排風阻板位于頭車流線型尾端位置，即D0處，然后沿著x負方向每隔12 m設(shè)置同樣的制動板。第1塊制動板記為ban1，第2塊板記為ban2，以此類推，最后一塊板記為ban8。列車以400 km/h的速度運行，各個制動板的制動力見圖6。

圖6 各制動板產(chǎn)生的制動力

由圖6可知，各個風阻制動板產(chǎn)生的制動力大小不一，首排風阻板即第1塊風阻板，產(chǎn)生的制動力最大，隨后的制動板產(chǎn)生的制動力越來越小。這是因為前排風阻板改變了氣流的流場結(jié)構(gòu)，氣流撞擊前排制動板能量損耗，產(chǎn)生速度相對較慢的繞流流場，因此后排制動板的制動效果相對較差。具體流場的壓力云圖和速度云圖見圖7。

圖7 列車流場結(jié)構(gòu)

3.2 首排風阻制動板制動力的研究

由圖6可知，首排風阻制動板產(chǎn)生的制動力最大，為總制動力的22%，因此對首排制動板制動力的優(yōu)化尤為重要。在頭車平順型車身上的不同位置設(shè)置風阻板，探究首排峰值制動板的最佳位置，以獲取最大的制動力。列車以400 km/h的速度運行，制動板產(chǎn)生的制動力見圖8。

圖8 首排制動板不同位置的制動力曲線

由圖8可知，整體來看，首排風阻制動板的位置越往后，產(chǎn)生的制動力越低。但是在距流線型尾端2～5 m范圍內(nèi)，制動力大小基本不變。為了更好的解釋現(xiàn)象，以不加制動板的工況作為對比。截取列車中間縱向?qū)ΨQ面的流場結(jié)構(gòu)為研究對象，不加制動板的頭車周圍的流場見圖9。

圖9 無制動板列車周圍的流場結(jié)構(gòu)(單位：Pa)

在對稱面上，沿著z正方向，距離列車車體上表面的0.2 m位置的x方向速度變化和壓力變化見圖10。記D0位置為原點，沿著與列車相反的方向為正方向。

圖10 距離車頂表面0.2 m處的速度與壓力曲線

風阻制動板的制動力來源主要來自于制動板前后的壓差。制動板迎風面產(chǎn)生正壓區(qū)，背風面由于氣流分離產(chǎn)生負壓區(qū)，從而產(chǎn)生制動力。從圖9可知，在頭車前端，氣流在頭車鼻尖處產(chǎn)生分離，由于空氣壓力作用在鼻尖處產(chǎn)生了正壓區(qū)；一部分氣流沿著流線型車頭高速流動，尤其在頭車頭部與頭車車身過渡區(qū)域，由于曲率較大，氣流流速較快，產(chǎn)生負壓區(qū)；在頭車車身位置，速度慢慢降低，壓力開始回升。在距離流線型尾端2 m處，即D2位置處，此時流線型對氣流的影響消失，但是此時依然處于負壓；D2位置后，由于平穩(wěn)型車身對氣流的黏性阻力起主要作用，氣流速度慢慢降低，壓力緩慢回升，在D3.7位置附近，氣流的速度等于列車運行速度，壓力值也在大氣壓值附近，這就導致在D3.7位置附近風阻制動板前后的壓力差值變化不大，見圖11。因此在距離尾端2～5 m處的制動力變化不大；由于黏性阻力的影響，車體表面附近的流場速度逐步降低至其穩(wěn)定值，壓力也慢慢升到其穩(wěn)定值，由于氣流分離出氣旋渦，使壓力和速度曲線呈現(xiàn)鋸齒狀。因此制動板的位置越往后，制動力就越小，最后制動力趨于穩(wěn)定。

圖11 不同位置制動板壓力云圖

3.3 首排風阻制動板對列車氣動性能的影響

頭車設(shè)置風阻制動板，改變了頭車以及整個列車的流場結(jié)構(gòu)，首排制動板不同位置對列車氣動性能的影響見圖12。

圖12 不同制動板位置的列車氣動性能

從圖12可以看出，首排風阻制動板的布置位置越靠后，頭車阻力越小。列車阻力主要有壓差阻力和黏性阻力構(gòu)成，從圖12可以看出，制動板布置位置靠前，頭車的壓差阻力較大，黏性阻力越小。因為氣體具有黏性，制動板會破壞列車表面附件氣流結(jié)構(gòu)，制動板的位置越往前制動板的影響范圍就越大，黏性阻力越小。中間車和尾車的氣動性能沒有改變。由圖12(b)可知，列車的總阻力變化不同于制動板制動力的變化趨勢，列車總阻力逐漸減小。這是由于制動板雖然在D2～D5段的制動力沒有發(fā)生變化，但是制動板迫使頭車阻力減小，從而導致總阻力下降。隨著制動板的后移，制動板的制動力和列車總阻力趨于平穩(wěn)。制動板的布置會使列車車身的阻力減小，但是變化量較小，制動板的制動力起主導作用。

4 結(jié)論

(1)列車采用風阻制動時，提供最大制動力的是首排風阻制動板，即位于頭車的第一套風阻制動板。

(2)首排風阻板的不同位置，制動效果差異較大。位于頭車流線型尾端，制動板產(chǎn)生的制動力最大，緊接著開始下降，最后，制動力保持不變。

(3)制動板的設(shè)置會改變車身附近的流場結(jié)構(gòu)。隨著制動板頭車的推后，頭車壓差阻力降低，黏性阻力上升，列車總阻力下降。