不同風(fēng)擋內(nèi)傾角度對(duì)高速列車氣動(dòng)性能的影響*

韓運(yùn)動(dòng)， 董 巖， 李傳迎， 汪海燕

(1 中車青島四方機(jī)車車輛股份有限公司 國(guó)家高速動(dòng)車組總成工程技術(shù)研究中心， 山東青島 266111；2 日照鋼鐵控股集團(tuán)有限公司 稽核室， 山東日照 276806；3 中車青島四方機(jī)車車輛股份有限公司 技術(shù)中心， 山東青島 266111)

2017年9月21日京滬高鐵率先以速度350 km/h開跑，未來中國(guó)高鐵的持續(xù)運(yùn)營(yíng)速度最高可達(dá)350 km/h。隨著我國(guó)高速列車運(yùn)行速度的提高，列車氣動(dòng)特性也日益凸顯。風(fēng)擋是列車之間的柔性部分，可以吸收車輛之間的所有相對(duì)運(yùn)動(dòng)并使旅客能安全容易的通過，同時(shí)為了保證客室環(huán)境的舒適度需滿足列車的空氣動(dòng)力學(xué)及聲學(xué)要求[1]。高速動(dòng)車組的風(fēng)擋分為外風(fēng)擋和內(nèi)風(fēng)擋兩部分，共同承擔(dān)車輛之間的擠壓和相對(duì)運(yùn)動(dòng)，并且內(nèi)風(fēng)擋主要作為旅客安全通道，外風(fēng)擋用來降低車輛運(yùn)行時(shí)產(chǎn)生的噪聲，同時(shí)還使兩個(gè)車輛之間弧面聯(lián)接平滑一致，以降低高速動(dòng)車組運(yùn)行時(shí)的空氣阻力。為了適應(yīng)車外氣壓波的急劇變化，滿足氣動(dòng)載荷下的強(qiáng)度要求，風(fēng)擋要有足夠的強(qiáng)度。風(fēng)擋作為高速列車車端連接的重要部件，其拓?fù)浣Y(jié)構(gòu)及連接方式對(duì)列車的氣動(dòng)性能、運(yùn)行平穩(wěn)性及旅客的乘坐舒適性有重要的影響[2]。列車在高速運(yùn)行時(shí)，其風(fēng)擋區(qū)域由于自身外形結(jié)構(gòu)因素會(huì)形成強(qiáng)烈的具有顯著非定常特性的湍流流場(chǎng)，由此產(chǎn)生的非定常氣動(dòng)載荷作用于風(fēng)擋部件上，引起風(fēng)擋結(jié)構(gòu)的振動(dòng)。長(zhǎng)期的氣動(dòng)振動(dòng)勢(shì)必造成風(fēng)擋結(jié)構(gòu)的疲勞損壞，影響列車運(yùn)行的安全性、穩(wěn)定性以及舒適度。尤其當(dāng)氣動(dòng)載荷與風(fēng)擋結(jié)構(gòu)產(chǎn)生共振時(shí)，風(fēng)擋結(jié)構(gòu)振動(dòng)幅度加大，并伴有低沉的噪聲，嚴(yán)重影響運(yùn)行安全性及乘客舒適度。

目前，國(guó)內(nèi)外對(duì)于高速列車風(fēng)擋外形結(jié)構(gòu)對(duì)列車氣動(dòng)性能影響的研究相對(duì)較少，對(duì)風(fēng)擋內(nèi)傾角度的研究更是寥寥無幾。楊加壽等[3]針對(duì)CRH380A型動(dòng)車組3節(jié)編組簡(jiǎn)化模型，在速度分別為360 km/h、420 km/h、500 km/h的條件下，通過數(shù)值計(jì)算研究了典型風(fēng)擋、全封閉外風(fēng)擋、無外風(fēng)擋3種情況以及基于典型風(fēng)擋構(gòu)建的4種其他風(fēng)擋形式對(duì)高速列車氣動(dòng)性能的影響；黃志祥等人[4]進(jìn)行了3節(jié)編組高速列車縮比模型的風(fēng)洞試驗(yàn)，研究了半封閉外風(fēng)擋和全封閉外風(fēng)擋結(jié)構(gòu)形式對(duì)列車空氣阻力的影響；BELL等[5]利用低噪聲大型風(fēng)洞研究了高速列車的氣動(dòng)噪聲，其研究結(jié)果表明車輛連接處的間隙是高速列車上主要的氣動(dòng)噪聲聲源部位，通過合理設(shè)計(jì)車體連接處的外形，根據(jù)不同的外形設(shè)計(jì)安裝平滑罩結(jié)構(gòu)能夠有效減小此處帶來的氣動(dòng)噪聲。

本文采用三維數(shù)值計(jì)算的方法，研究?jī)蓚?cè)外風(fēng)擋的內(nèi)傾角度變化對(duì)列車氣動(dòng)性能的影響，數(shù)值計(jì)算的準(zhǔn)確性通過線路實(shí)車試驗(yàn)進(jìn)行了驗(yàn)證。研究得到了用于指導(dǎo)高速列車風(fēng)擋氣動(dòng)設(shè)計(jì)的重要結(jié)論，為解決高速列車風(fēng)擋振動(dòng)問題指明了方向。

1 數(shù)值計(jì)算方法

1.1 計(jì)算域、邊界條件設(shè)置

原風(fēng)擋結(jié)構(gòu)(風(fēng)擋內(nèi)傾角度定義為0°，下同。)三維模型見圖1所示。風(fēng)擋內(nèi)傾角度變化定義說明見圖2(a)所示，外側(cè)圓角在保證風(fēng)擋總高度和厚度以及下部寬度不變的情況下，以風(fēng)擋安裝孔為旋轉(zhuǎn)基點(diǎn)，分別旋轉(zhuǎn)2°、4°、6°、8°獲得不同風(fēng)擋內(nèi)傾角度，外風(fēng)擋內(nèi)傾0°、4°、8°對(duì)比示意見圖2(b)所示。

圖1 原結(jié)構(gòu)三維模型

圖2 風(fēng)擋內(nèi)傾角度變化說明

為了減少計(jì)算量，采用頭車+中間車向后延伸的計(jì)算編組形式，對(duì)列車幾何模型進(jìn)行適當(dāng)簡(jiǎn)化，頭車計(jì)算模型包含轉(zhuǎn)向架等重要結(jié)構(gòu)。在滿足計(jì)算要求的前提下，計(jì)算區(qū)域應(yīng)選取適當(dāng)?shù)拇笮。瑸榱瞬挥绊戯L(fēng)擋區(qū)域的流場(chǎng)發(fā)展，將車體向后延伸。最終搭建的計(jì)算域見圖3所示。

圖3 計(jì)算域

1.2 網(wǎng)格劃分

采用流體仿真軟件STAR-CCM+5.04對(duì)空間區(qū)域的網(wǎng)格進(jìn)行劃分，空間網(wǎng)格類型為Trimmer網(wǎng)格，車體表面布置邊界層網(wǎng)格，第1層高度為1 mm，增長(zhǎng)比為1.5，共6層邊界層網(wǎng)格；考慮到地面邊界層的影響，在地面同樣布置邊界層網(wǎng)格，第1層高度為4 mm，增長(zhǎng)比為1.5，共6層邊界層網(wǎng)格；為更好的模擬列車周圍的流場(chǎng)結(jié)構(gòu)，在幾何突變區(qū)域及尾流場(chǎng)和背風(fēng)側(cè)進(jìn)行網(wǎng)格加密，如圖4所示對(duì)風(fēng)擋區(qū)域設(shè)置網(wǎng)格加密域，加密域的尺寸為50 mm，風(fēng)擋結(jié)構(gòu)的面網(wǎng)格最小尺寸為2 mm，整體網(wǎng)格量約為108。風(fēng)擋結(jié)構(gòu)的體網(wǎng)格劃分結(jié)果如圖5所示。

圖4 頭車車體表面及風(fēng)擋區(qū)域周圍網(wǎng)格

圖5 風(fēng)擋網(wǎng)格劃分

1.3 計(jì)算方法

本文分析的列車運(yùn)行速度為350 km/h，來流馬赫數(shù)小于0.3，因此，計(jì)算模型采用不可壓計(jì)算方法；給定速度入口邊界，來流速度為運(yùn)行速度350 km/h，給定壓力出口邊界，車體為無滑移固壁邊界條件；為模擬地面效應(yīng)的影響，將地面設(shè)置為移動(dòng)壁面，移動(dòng)速度與來流速度相等。

計(jì)算方法采用基于κ-ωSST兩方程湍流模型的RANS方法。控制方程為定常N-S方程，方程的離散采用有限體積法。

2 方法驗(yàn)證

圖6給出了前期線路實(shí)車測(cè)試風(fēng)擋內(nèi)壓力測(cè)點(diǎn)以及振動(dòng)加速度測(cè)點(diǎn)的布置圖。圖7給出了仿真計(jì)算的壓力測(cè)點(diǎn)布置圖。

圖8給出了速度300 km/h速度下線路實(shí)車測(cè)試測(cè)點(diǎn)的壓力和振動(dòng)加速度的頻譜分析結(jié)果。圖9給出了同一速度下、相同位置測(cè)點(diǎn)仿真計(jì)算的壓力頻譜分析結(jié)果。由圖8可以看出，線路測(cè)試壓力頻譜分析結(jié)果和振動(dòng)加速度頻譜分析結(jié)果一致，存在一個(gè)8.118 Hz的主頻。由圖9可以看出，仿真計(jì)算壓力的頻譜分析結(jié)果同樣存在一個(gè)8.161 6 Hz的主頻，這與線路實(shí)車測(cè)試結(jié)果吻合，驗(yàn)證了仿真計(jì)算的準(zhǔn)確性和可靠性。

圖6 實(shí)車測(cè)試車間風(fēng)擋測(cè)點(diǎn)布置示意圖

圖7 仿真計(jì)算風(fēng)擋壓力測(cè)點(diǎn)布置圖

3 計(jì)算結(jié)果與分析

3.1 分析說明

風(fēng)擋結(jié)構(gòu)編號(hào)說明如圖10所示，原風(fēng)擋背風(fēng)側(cè)編號(hào)分別命名為gap1_01～gap1_06，對(duì)應(yīng)的迎風(fēng)側(cè)編號(hào)分別為gap2_01～gap2_06。風(fēng)擋區(qū)域的壓力測(cè)點(diǎn)布置如圖11所示。

表1 實(shí)車測(cè)試主要設(shè)備

圖8 實(shí)車測(cè)試壓力與加速度頻譜分析結(jié)果

圖9 仿真計(jì)算壓力頻譜分析結(jié)果

圖10 原風(fēng)擋編號(hào)說明示意圖

3.2 風(fēng)擋壓力速度云圖

圖12～圖13分別給出了迎風(fēng)側(cè)與背風(fēng)側(cè)不同風(fēng)擋內(nèi)傾角度下表面壓力云圖。由不同結(jié)構(gòu)的壓力云圖可以看出，風(fēng)擋間隙周邊存在相對(duì)較大正壓區(qū)，這是由于氣流由間隙加入流入所致。原風(fēng)擋方案，下風(fēng)擋和角風(fēng)擋存在負(fù)壓區(qū)，其他均為正壓區(qū)。風(fēng)擋內(nèi)傾2°、4°、6°、8°方案，整個(gè)風(fēng)擋區(qū)域基本處于正壓區(qū)，內(nèi)傾位置所受壓力值較大。

圖11 風(fēng)擋區(qū)域壓力測(cè)點(diǎn)布置圖

圖14給出了風(fēng)擋內(nèi)傾4°情況下風(fēng)擋區(qū)域的速度流線圖。由圖14可以看出，風(fēng)擋區(qū)域的流動(dòng)極為復(fù)雜，高速列車風(fēng)擋由內(nèi)外風(fēng)擋組成，兩者之間形成近似環(huán)狀的空腔，氣流由風(fēng)擋間隙流入空腔形成空腔流動(dòng)，空腔內(nèi)部的流動(dòng)為復(fù)雜的湍流，伴隨有各種渦結(jié)構(gòu)的脫落和破碎，導(dǎo)致內(nèi)部氣流及內(nèi)外風(fēng)擋表面壓力產(chǎn)生不規(guī)則的脈動(dòng)變化，引起風(fēng)擋的變形。其變形量的變化主要受內(nèi)部氣流的擾動(dòng)程度影響，空腔變小，流進(jìn)的氣流減少，影響空腔內(nèi)的渦結(jié)構(gòu)，流速變小，氣流壓力脈動(dòng)減小，一定程度上緩解了外風(fēng)擋的受力情況。

圖12 gap1背風(fēng)側(cè)不同內(nèi)傾角度下風(fēng)擋結(jié)構(gòu)表面的壓力云圖

圖13 gap2迎風(fēng)側(cè)不同內(nèi)傾角度下風(fēng)擋結(jié)構(gòu)表面的壓力云圖

圖14 風(fēng)擋內(nèi)傾4°風(fēng)擋區(qū)域的流線圖

3.3 風(fēng)擋測(cè)點(diǎn)壓力分析

表2給出了風(fēng)擋區(qū)域不同測(cè)點(diǎn)的壓力值。由表2可以看出，相對(duì)于原風(fēng)擋，風(fēng)擋內(nèi)傾8°各測(cè)點(diǎn)的壓力值均增大；而內(nèi)傾4°方案，除個(gè)別測(cè)點(diǎn)外，壓力基本增大。內(nèi)傾8°方案風(fēng)擋區(qū)域各測(cè)點(diǎn)的壓力值均為正壓。

3.4 風(fēng)擋氣動(dòng)力分析

下面主要從氣動(dòng)力的角度分析風(fēng)擋的變形。計(jì)算得到了不同風(fēng)擋結(jié)構(gòu)、不同風(fēng)擋塊在3個(gè)方向所受的阻力、升力、側(cè)向力，計(jì)算結(jié)果見表3所示。對(duì)于兩側(cè)風(fēng)擋主要分析其所受的側(cè)向力情況，對(duì)于上下風(fēng)擋主要分析其所受的升力情況。分背風(fēng)側(cè)(命名為gap1)和迎風(fēng)側(cè)(命名為gap2)。

由表3可以看出，對(duì)于原風(fēng)擋結(jié)構(gòu)，兩側(cè)風(fēng)擋gap1_01、gap1_02、gap2_01、gap2_02所受側(cè)向力方向均指向外側(cè)，背風(fēng)側(cè)兩側(cè)風(fēng)擋gap1_01、gap1_02(272 N、-383 N)所受側(cè)向力明顯大于迎風(fēng)側(cè)兩側(cè)風(fēng)擋gap2_01、gap2_02(81 N、-96 N)；迎風(fēng)側(cè)下風(fēng)擋與上風(fēng)擋gap2_03、gap2_04所受的升力較小，分別為-45 N、41 N，方向均指向外側(cè)。

風(fēng)擋內(nèi)傾4°方案：對(duì)比原風(fēng)擋，背風(fēng)側(cè)兩側(cè)風(fēng)擋gap1_01、gap1_02所受的側(cè)向力均變大，分別為452 N和-468 N，分別增大約21.5%和22.2%，方向不變，均指向外側(cè)。迎風(fēng)側(cè)兩側(cè)風(fēng)擋gap2_01、gap2_02所受側(cè)向力均變大，分別為236 N和-250 N，分別增大約191%和160%，受力方向不發(fā)生變化。背風(fēng)側(cè)與迎風(fēng)側(cè)上、下風(fēng)擋gap1_03、gap1_04、 gap2_03、gap2_04所受的升力均顯著增大，方向均指向外側(cè)。

風(fēng)擋內(nèi)傾8°方案：對(duì)比原風(fēng)擋，背風(fēng)側(cè)兩側(cè)風(fēng)擋gap1_01、gap1_02所受的側(cè)向力均變小，分別為340 N和-355 N，分別減小約8.6%和7.3%，方向不變，均指向外側(cè)。迎風(fēng)側(cè)兩側(cè)風(fēng)擋gap2_01、gap2_02所受側(cè)向力均變小，分別為-31 N和33 N，方向發(fā)生改變，均指向內(nèi)側(cè)。背風(fēng)側(cè)與迎風(fēng)側(cè)上、下風(fēng)擋gap1_03、gap1_04、 gap2_03、gap2_04所受的升力均變小，迎風(fēng)側(cè)上、下風(fēng)擋升力明顯小于背風(fēng)側(cè)上、下風(fēng)擋，方向發(fā)生改變，均指向內(nèi)側(cè)。

表2 不同測(cè)點(diǎn)壓力計(jì)算結(jié)果

表3 不同風(fēng)擋結(jié)構(gòu)側(cè)向力計(jì)算結(jié)果

圖15給出了不同風(fēng)擋內(nèi)傾角度情況下，兩側(cè)風(fēng)擋所受側(cè)向力對(duì)比圖。由圖可以看出，兩側(cè)風(fēng)擋所受側(cè)向力對(duì)稱性較好。不同內(nèi)傾角度工況下，背風(fēng)側(cè)兩側(cè)風(fēng)擋所受側(cè)向力方向均指向外側(cè)，呈現(xiàn)“外推”狀態(tài)。對(duì)于迎風(fēng)側(cè)兩側(cè)風(fēng)擋，0°、2°、4°工況下所受側(cè)向力方向指向外側(cè)，呈現(xiàn)“外推”狀態(tài)，而6°、8°工況下所受側(cè)向力方向指向內(nèi)側(cè)，呈現(xiàn)“內(nèi)壓”狀態(tài)。可見，風(fēng)擋區(qū)域復(fù)雜的流動(dòng)導(dǎo)致兩側(cè)風(fēng)擋所受側(cè)向力與內(nèi)傾角度并不是線性關(guān)系。

圖15 不同風(fēng)擋內(nèi)傾角度兩側(cè)風(fēng)擋側(cè)向力對(duì)比

4 結(jié) 論

通過研究不同風(fēng)擋內(nèi)傾角度的變化對(duì)列車氣動(dòng)性能的影響，可以獲得如下結(jié)論：

(1)相對(duì)于原風(fēng)擋，除個(gè)別測(cè)點(diǎn)外，風(fēng)擋內(nèi)傾2°、4°、6°、8°各測(cè)點(diǎn)的壓力值均增大；內(nèi)傾6°、8°風(fēng)擋區(qū)域各測(cè)點(diǎn)的壓力值均為正壓。

(2)兩側(cè)風(fēng)擋所受側(cè)向力對(duì)稱性較好，不同內(nèi)傾角度，背風(fēng)側(cè)兩側(cè)風(fēng)擋所受側(cè)向力方向均指向外側(cè)，呈現(xiàn)“外推”狀態(tài)；迎風(fēng)側(cè)兩側(cè)風(fēng)擋，0°、2°、4°所受側(cè)向力方向指向外側(cè)，呈現(xiàn)“外推”狀態(tài)，而6°、8°所受側(cè)向力方向指向內(nèi)側(cè)，呈現(xiàn)“內(nèi)壓”狀態(tài)。

(3)兩側(cè)風(fēng)擋所受側(cè)向力從小到大順序?yàn)?°<8°<0°<4°<2°，風(fēng)擋區(qū)域復(fù)雜的流動(dòng)導(dǎo)致兩側(cè)風(fēng)擋所受側(cè)向力與內(nèi)傾角度并不是線性關(guān)系。

研究結(jié)果為高速列車風(fēng)擋結(jié)構(gòu)的氣動(dòng)設(shè)計(jì)提供了技術(shù)支撐，為解決高速列車風(fēng)擋振動(dòng)變形問題指明了新的方向。