高陽，李啟良，王毅剛

(1.中車長春軌道客車股份有限公司 工程實驗室，吉林 長春 130062； 2.同濟大學(xué) 上海地面交通工具風(fēng)洞中心，上海 201804)*

全尺寸高速列車頭型氣動噪聲分析方法

高陽1，李啟良2，王毅剛2

(1.中車長春軌道客車股份有限公司 工程實驗室，吉林 長春 130062； 2.同濟大學(xué) 上海地面交通工具風(fēng)洞中心，上海 201804)*

建立高速列車頭型氣動噪聲分析方法有利于了解頭型與空氣相互作用產(chǎn)生的氣動噪聲特性及其對車內(nèi)外的影響.為此，先后建立了兩個頭型的1∶8縮比三節(jié)編組氣動噪聲仿真模型，并開展氣動噪聲仿真計算，得到外場測點平均總聲壓級.通過與風(fēng)洞試驗結(jié)果相比較，兩者量值相差小于3 dBA，且頭型1均小于頭型2，驗證數(shù)值仿真結(jié)果.為了實現(xiàn)全尺寸高速列車頭型氣動噪聲數(shù)值仿真，提出在三節(jié)編組的計算域中截取一部分—子域法.子域法和整車得到頭型部位的氣動特性一致性間接表明子域法的合理性.利用子域法開展了全尺寸頭型1和頭型2氣動噪聲仿真計算，得到頭型表面聲功率、表面和外場總聲壓級，可為頭型選型和優(yōu)化提供依據(jù)，從而建立了基于數(shù)值仿真的全尺寸高速列車頭型氣動噪聲分析方法，解決了以往無法通過風(fēng)洞試驗和數(shù)值仿真進行全尺寸高速列車頭型氣動噪聲分析.

高速列車；頭型；氣動噪聲；子域法

0 引言

高速列車產(chǎn)生的氣動噪聲對周圍居民的生產(chǎn)和生活環(huán)境造成很大的影響[1- 2]，也影響乘員乘坐的舒適性[3]，為此受到研究人員廣泛重視.高速列車氣動噪聲的研究方法主要有實車路試、聲學(xué)風(fēng)洞試驗和數(shù)值計算三種.相比于路試和風(fēng)洞試驗，數(shù)值計算的方法能夠為列車開發(fā)階段的氣動噪聲評估和降噪提供更多的數(shù)據(jù)支持.高速列車氣動噪聲的數(shù)值計算均是基于Lighthill聲類比方法，早期研究集中在門窗、絕緣子、受電弓桿件部件[4- 7].隨著計算能力的提升，學(xué)者們開始采用數(shù)值方法評估三車編組列車的氣動噪聲分布及聲源特性[8- 11]，并針對轉(zhuǎn)向架[12- 13]開展降噪研究.

頭車是高速列車運行中重要的氣動噪聲源，頭型的參數(shù)影響高速列車氣動噪聲等各種性能，合理的頭型氣動噪聲分析方法有助于了解頭車與空氣相互作用和產(chǎn)生氣動噪聲的機理.肖友剛[14- 15]分析了頭車曲面外形對列車車內(nèi)及遠場氣動噪聲的影響.劉加利[16]計算了頭尾車組成的簡化模型遠場氣動噪聲特性，但都缺少與聲學(xué)風(fēng)洞試驗的直接或間接對比驗證.

為此本文首先開展了1∶8縮比三節(jié)編組氣動噪聲仿真分析和風(fēng)洞試驗研究，擬通過對比外場測點總聲壓級來驗證數(shù)值仿真準(zhǔn)確性.在此基礎(chǔ)上提出在三節(jié)編組的計算域中截取一部分——子域法，希望通過子域法來實現(xiàn)全尺寸高速列車頭型氣動噪聲數(shù)值仿真，從而彌補以往無法通過風(fēng)洞試驗和數(shù)值仿真進行全尺寸高速列車頭型氣動噪聲分析.

1 數(shù)值仿真方法的驗證

1.1 風(fēng)洞試驗方法與數(shù)值仿真方法

為了進行頭型的選型，受風(fēng)洞大小的限制.目前國內(nèi)外通常只能進行1∶8或更小縮比的風(fēng)洞試驗.為了驗證數(shù)值仿真方法，選取了2個頭型在綿陽二九基地的5.5 m×4.0 m的開口聲學(xué)風(fēng)洞進行1∶8縮比三節(jié)編組的氣動噪聲試驗.試驗?zāi)Ｐ桶惭b在試驗段的支撐平臺上，并采用地板裝置支撐.頭車模型的鼻尖距地板斜坡為1.0 m，如圖1所示.模型安裝到位后，對模型的姿態(tài)角進行標(biāo)定，確保試驗?zāi)Ｐ妥藨B(tài)角準(zhǔn)確.

圖1 1∶8縮比三節(jié)編組試驗?zāi)Ｐ?/p>

車外的聲場環(huán)境是評估頭型噪聲的主要指標(biāo)之一，也是頭型選型的主要依據(jù)，為此在試驗?zāi)Ｐ妥髠?cè)布置共計30個遠場麥克風(fēng).遠場麥克風(fēng)距車身中心截面6.25 m，分3排布置.首排測點高度距地板0.4 m，第二排和第三排分別距地板0.8 m和1.2 m.在流動方向上，兩個麥克風(fēng)相距0.8 m，如圖2所示.

圖2 測點示意圖

試驗?zāi)Ｐ秃吐晫W(xué)測量設(shè)備安裝完成后，采用聲學(xué)校準(zhǔn)器校準(zhǔn)所有麥克風(fēng)，并使用激光器確保模型處于0°.當(dāng)試驗風(fēng)速穩(wěn)定在250 km/h時，設(shè)定采樣頻率和采樣時間分別為51.2 kHz和15 s，由動態(tài)數(shù)據(jù)采集系統(tǒng)進行數(shù)據(jù)采集.采集到時域信號通過快速傅里葉變化得到測點頻譜和總聲壓級.

建立與風(fēng)洞試驗完全一致的數(shù)值仿真模型，包括試驗地板、軌道、轉(zhuǎn)向架、受電弓、車廂風(fēng)擋等所有細節(jié)部件.為了更好模擬列車流場和聲場，創(chuàng)建如圖3所示的計算域.計算域的長、寬和高分別為18.5 m、5.5 m和2.5 m.計算域進口距離地板前緣2.5 m，設(shè)定與風(fēng)洞試驗相同的速度和湍流度，即250 km/h和0.2%.計算域出口距離地板后緣5.5 m，設(shè)定相對壓力為0 Pa的壓力出口.計算域兩側(cè)距模型中心2.75 m，頂部距試驗段地面2.5 m.它們均設(shè)定為對稱邊界.所有固定不動的壁面均設(shè)定為無滑移壁面.

圖3 計算域

依次對1∶8縮比的三節(jié)編組進行面網(wǎng)格劃分.頭車的鼻錐、排障器、裙板和雨刮凹槽的面網(wǎng)格大小為2.5 mm，其它區(qū)域的面網(wǎng)格大小為5.0 mm，生成約24萬個三角形面網(wǎng)格.拖車轉(zhuǎn)向架和動力轉(zhuǎn)向架的面網(wǎng)格大小分別為2 mm和3 mm，分別生成約12萬個和31萬個三角形面網(wǎng)格.受電弓保留受電弓頭、弓架、基座和支撐絕緣子，它們的面網(wǎng)格大小為0.5 mm，生成24萬個三角形面網(wǎng)格.保留車廂連接部位風(fēng)擋的具體結(jié)構(gòu)，采用3.0 mm的面網(wǎng)格大小，生成3.5萬個三角形面網(wǎng)格.直流段面網(wǎng)格大小為5 mm.軌道面網(wǎng)格為2 mm.靠近軌道和遠離軌道的地板面網(wǎng)格分別為10 mm和30 mm.整車面網(wǎng)格數(shù)量約為500萬個.2個頭型對應(yīng)整車面網(wǎng)格數(shù)相差4萬個，差異不到1%.

使用STARCCM+的Trim+layer策略生成計算域的體網(wǎng)格.在距車身0.3 m創(chuàng)建包絡(luò)面，用于實現(xiàn)將列車、軌道和地板表面產(chǎn)生的偶極子和空間渦流產(chǎn)生四極子噪聲輻射到外場.包絡(luò)面里側(cè)和外側(cè)體網(wǎng)格大小分別為20 mm和40 mm.基于可壓縮的DES湍流模型對近壁區(qū)網(wǎng)格要求，生成了第一層為0.5 mm，增長率為1.15，共計7層邊界層網(wǎng)格.整個計算域共計生成約3 200萬個體網(wǎng)格，滿足網(wǎng)格質(zhì)量要求.不同頭型對應(yīng)體網(wǎng)格差異小于1%.

氣動噪聲計算基于STARCCM+框架中的可壓縮DES，其近壁面采用k-ω湍流模型，遠離壁面采用LES.計算首先使用SSTk-ω模型得到定常解.然后基于該定常解，使用可壓縮DES進行氣動噪聲計算.設(shè)定非定常計算的時間步長和每時間步長迭代次數(shù)分別為0.000 25s和25.考慮計算精度和計算效率，共計算物理時間0.375s，最終的0.25s用于數(shù)據(jù)分析.應(yīng)該指出的是，每時間步長迭代步數(shù)都是經(jīng)過分析殘差、監(jiān)測受電弓表面壓力系數(shù)等參數(shù)來確定.所設(shè)定的時間步長能保證流動準(zhǔn)周期內(nèi)有足夠樣本點.所有計算均在同濟大學(xué)的512核高性能計算機集群下完成.每個工況在120核的刀片服務(wù)器下，花費7天.

1.2 總聲壓級對比

圖4給出3排測點通過風(fēng)洞實驗和數(shù)值仿真得到測點總聲壓級.測點1～10為第1排，測點11～20為第2排，測點21～30為第3排.從圖中可以看出，所有測點的總聲壓級處于76～85dBA之間.不同離地高度的各排測點總聲壓級稍有差異，但變化不大.

圖4 測點總聲壓級

對于第一排測點，頭型2對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為80.62dBA和78.62dBA，兩者相差2.00dBA.頭型1對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為82.78dBA和80.25dBA，兩者相差2.52dBA.對于第二排測點，頭型2對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為80.69dBA和78.68dBA，兩者相差2.01dBA.頭型1對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為82.83dBA和80.41dBA，兩者相差2.41dBA.對于第三排測點，頭型2對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為80.71dBA和76.96dBA，兩者相差3.75dBA.頭型1對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為82.81dBA和78.85dBA，兩者相差3.96dBA.

基于三排測點，頭型2對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為80.67dBA和78.09dBA，兩者相差2.59dBA.頭型1對應(yīng)測點平均總聲壓級的仿真和試驗結(jié)果分別為82.81dBA和79.84dBA，兩者相差2.97dBA.由此可見，無論是仿真還是試驗，均表明在三節(jié)編組情況下，頭型2優(yōu)于頭型1.數(shù)值仿真方法可信.

2 實車頭型對比

2.1 仿真方法

受整車風(fēng)洞建設(shè)和運行成本限制，相當(dāng)長時間都無法實現(xiàn)實車三節(jié)編組的氣動噪聲風(fēng)洞試驗.為此有必要尋找實車頭型氣動噪聲的仿真方法，來解決模型縮比細節(jié)缺失和雷諾數(shù)低的問題.

為了滿足1∶8縮比三節(jié)編組氣動噪聲要求，創(chuàng)建了高達3 000萬個體網(wǎng)格，若進行實車三節(jié)編組氣動噪聲計算，保守估計網(wǎng)格數(shù)也將超過3億.這對于絕大多數(shù)研究人員都是無法滿足.為了解決網(wǎng)格數(shù)過大、計算時間過長，本文提出用于解決實車頭型氣動噪聲計算的方法-子域法.所謂子域法是只選取整車計算的一部分，并在此基礎(chǔ)上創(chuàng)建較小的計算域，具體如圖5所示.

圖5 子域法

子域進口距頭型鼻尖約6m，出口在頭車直流段的一半；兩側(cè)各距車廂中截面8.4m，約2.2倍車廂寬度；底部為地面，頂部距地面12m.整個子域的長、寬和高分別為25.3m，16.8m和12.0m.為了便于對比，直接將子域固壁表面的面網(wǎng)格放大8倍.子域四周邊界面網(wǎng)格為80mm.使用STARCCM+的Trim+layer策略生成計算域的體網(wǎng)格.在距車身1.2m創(chuàng)建包絡(luò)面，用于實現(xiàn)將固壁表面產(chǎn)生的偶極子和空間渦流產(chǎn)生四極子噪聲輻射到外場.包絡(luò)面里側(cè)和外側(cè)體網(wǎng)格大小分別為20mm和80mm.基于可壓縮的DES湍流模型對近壁區(qū)網(wǎng)格要求，生成了第一層為0.5mm，增長率為1.15，共計7層邊界層網(wǎng)格.整個子域共計生成約3 000萬個體網(wǎng)格，滿足網(wǎng)格質(zhì)量要求.不同頭型對應(yīng)體網(wǎng)格差異小于1%.

使用子域法進行氣動噪聲計算，不能按照通常在計算域進出口指定均勻的邊界條件，而應(yīng)利用實車三節(jié)編組得到特定截面的速度、湍流結(jié)果作為邊界條件.由于目前無法進行實車三節(jié)編組氣動噪聲計算，但可以進行氣動計算.因此可通過實車三節(jié)編組氣動計算得到子域兩側(cè)、頂部和前部的速度分量、湍流強度和耗散率，以及后部的靜壓、湍流強度和耗散率.子域法的其它計算設(shè)置與1∶8縮比三節(jié)編組氣動噪聲計算設(shè)置完全相同，在此不再敘述.

2.2 子域法合理性

由于缺少實車三節(jié)編組氣動噪聲結(jié)果，因此評估子域法的合理性不能直接拿兩者氣動噪聲結(jié)果進行對比.將通過兩者氣動力、表面壓力和中截面速度等氣動信息來間接說明子域法的合理性.

實車三節(jié)編組氣動計算的面網(wǎng)格來自于縮比模型的直接放大.邊界層層數(shù)由原來的7層減少到4層、體網(wǎng)格大小由原來的20增加到100，使之最終網(wǎng)格數(shù)為3 000萬左右.其它計算設(shè)置與縮比模型氣動噪聲計算中定常部分完全相同.通過實車三節(jié)編組氣動計算可以得到頭型1和頭型2、直流段和轉(zhuǎn)向架的氣動阻力系數(shù)分別為0.087和0.095，而子域法得到氣動阻力系數(shù)分別為0.090和0.093.對于相同頭型，兩者相差小于3%，而且均表明頭型1小于頭型2.

圖6給出兩者計算得到表面靜壓和中截面速度分布.兩者一致性表明子域法得到的氣動仿真結(jié)果能夠代表真實情況，具有一定合理性.

(a)頭型1

(b)頭型2

2.3 頭型對比

2.3.1 外場測點總聲壓級

與縮比模型相同，仍通過外場測點平均總聲壓級來進行頭型的選型.對于實車外場噪聲的評估，通常選擇距車身中截面25 m.首排測點位于距地面高度0.5 m處，測點間隔1 m，序號標(biāo)定為1～12.其中，序號12正對頭型末端與直流段連接處.第二排和第三排測點分別距地面1.5 m和2.5 m，序號標(biāo)定分別為13～24和25～36.

從圖7的測點總聲壓級分布可以看到，所有測點的總聲壓級處于80～88 dBA，流動方向后端比前端要大，最大差異約5 dBA.不同離地高度的測點，總聲壓級稍有差異，但變化不大.頭型1和頭型2對應(yīng)測點平均總聲壓級分別為82.4 dBA和85.6 dBA.兩者排序與1∶8縮比三節(jié)編組得到的排序相反.其主要原因是1∶8縮比三節(jié)編組外場測點給出的總聲壓級是整車、軌道、地板和地面等多個聲源共同結(jié)果，而實車頭型給出的僅是單個頭型、直流段和地面的聲源，可見聲源類型和個數(shù)遠小于前者，更能反映出不同頭型差異導(dǎo)致氣動噪聲的不同.

圖7 測點總聲壓級對比

圖8給出外場36個測點在各自頻率的能量平均聲壓級.從圖中可以看到，頭型1輻射到外場測點的聲壓級幾乎在所有頻率都低于頭型2.在250 Hz以下和500 Hz以上差異較為明顯.

圖8 測點聲壓級頻譜對比

2.3.2 表面聲功率與總聲壓級

評估頭型好壞，除了考慮外場噪聲外，還需考慮作用在車身鈑金件的激勵，并最終導(dǎo)致車內(nèi)噪聲的大小.表面聲功率的大小反映了作用在表面激勵的聲能量.為此分別計算了鼻尖、排障器、轉(zhuǎn)向架、曲線部位和直線部位的聲功率，并由此得到它們的總聲功率.通過計算得到頭型1表面聲功率比頭型2小1 dB，仍能表明頭型1優(yōu)于頭型2.應(yīng)該指出的是，分析1∶8縮比三節(jié)編組在相同部位的表面聲功率可以看到，頭型1對應(yīng)的表面聲功率也小于頭型2.

圖9分別給出2個頭型表面總聲壓級分布.從圖中可以看到轉(zhuǎn)向架是總聲壓級最大的區(qū)域，表明它是主要聲源，占有最大的聲能量.排障器兩側(cè)的總聲壓級低于轉(zhuǎn)向架，是不能忽略的部件.曲線部位的總聲壓級也不小，且面積大.由于雨刮凹槽區(qū)域正對著司機室.在密封較好的情況下，該區(qū)域總聲壓級不能忽略.頭型1上部總聲壓級明顯低于鳳頭型2，但底部卻不如頭型1，特別是排障器區(qū)域.

(a)頭型1

(b)頭型2

圖9 總聲壓級對比

3 結(jié)論

為了建立全尺寸高速列車頭型氣動噪聲分析方法，開展了1∶8縮比三節(jié)編組氣動噪聲風(fēng)洞試驗和數(shù)值仿真，得到2個頭型在外場測點的總聲壓級.在此基礎(chǔ)上，建立了子域法，并通過子域法開展了實車頭型氣動噪聲研究.通過上述研究得到結(jié)論如下：

(1)建立了與風(fēng)洞試驗相一致的1∶8縮比三節(jié)編組氣動噪聲仿真模型，得到了與風(fēng)洞試驗相一致的氣動噪聲結(jié)果：基于外場測點平均總聲壓級，數(shù)值仿真和風(fēng)洞試驗得到外場側(cè)窗平均總聲壓級相差小于3 dBA，均發(fā)現(xiàn)頭型1優(yōu)于頭型2;

(2)將子域法引入到全尺寸高速列車頭型氣動噪聲仿真中.子域法和整車得到氣動阻力相差小于3%，靜壓系數(shù)和無量綱速度分布一致性間接表明子域法的合理性;

(3)利用子域法能夠獲得2個頭型的外場總聲壓級、表面聲功率級和總聲壓級分布，可為頭型選型和優(yōu)化提供依據(jù)，從而建立了全尺寸高速列車頭型氣動噪聲仿真分析方法.

[1]繆炳榮，張衛(wèi)華，鄧永權(quán)，等.新一代中國高速鐵路動車組面臨的技術(shù)挑戰(zhàn)與策略研究[J].中國工程科學(xué)，2015(4)：98- 112.

[2]田紅旗.中國高速軌道交通空氣動力學(xué)研究進展及發(fā)展思考[J].中國工程科學(xué)，2015(4)：30- 41.

[3]范蓉平，孫旭，孟光,等.高速列車車內(nèi)噪聲特性研究[J].振動工程學(xué)報,2004，17(2)：1097- 1110.

[4]TALOTTE C,GAUTIER P E,THOMPSON D J,et al.Identification,modeling and reduction potential of railway noise sources: a critical suvey[J].Journal of Sound and Vibration,2003,267(3):447- 468.

[5]TALOTTE C.Aerodynamic noise:a critical survey[J].Journal of Sound and Vibration,2000,231(3):549- 662.

[6]SASSA T,SATO T,YATSHI S.Numerical analysis of aerodynamic noise radiated from a high speed train surface[J].Journal of Sound and Vibration,2001,47(3):407- 416.

[7]YU H,LI J,ZHANG H.On aerodynamic noises radiated by the pantograph system of high-speed trains[J].Acta Mechanica Sinica,2013,29(3):399- 410.

[8]孫振旭，王一偉，安亦然.高速列車氣動噪聲的計算研究[J].水動力學(xué)研究與進展，2010，25(5):660- 667.

[9]鄭拯宇，李人憲.高速列車表面氣動噪聲偶極子聲源分布數(shù)值分析[J].西南交通大學(xué)學(xué)報(自然科學(xué)版)，2011，46(6)：996- 1002.

[10]張軍，孫幫成，郭濤，等.高速列車整車氣動噪聲及分布規(guī)律研究[J].鐵道學(xué)報，2015，37(2)：10- 17.

[11]張軍，郭濤，孫幫成，等.高速列車氣動噪聲源特性研究[J].鐵道學(xué)報，2015，37(6)：10- 18.

[12]黃莎，楊明智，李志偉.高速列車轉(zhuǎn)向架部位氣動噪聲數(shù)值模擬及降噪研究[J].中南大學(xué)學(xué)報(自然科學(xué)版)，2011，42(12):3899- 3904.

[13]ZHU J Y,HU Z H,THOMPSON D J.Flow simulation and aerodynamic noise prediction for a high-speed train wheelset[J].International Journal of Aeroacoustics,2014,13:533- 552.

[14]肖友剛，田紅旗，張洪.高速列車司機室內(nèi)氣動噪聲預(yù)測[J].交通運輸工程學(xué)報，2008，8(3):10- 14.

[15]肖友剛，康志成.高速列車車頭曲面氣動噪聲的數(shù)值預(yù)測[J].中南大學(xué)學(xué)報(自然科學(xué)版)，2008，39(6):1267- 1272.

[16]劉加利，張繼業(yè)，張衛(wèi)華.高速列車車頭氣動噪聲數(shù)值分析[J].鐵道學(xué)報，2011，33(9):19- 26.

Analysis Method of Aerodynamic Noise of Full Scale High Speed Train Head Shape

GAO Yan1，LI Qiliang2，WANG Yigang2

(1.CRRC Changchun Railway Vehicles Limited Company Engineering Laboratory,Changchun 130062,China； 2.Tongji University-Shanghai Automotive Wind Tunnel Center,Shanghai 201804,China)

The establishment of an aerodynamic noise analysis method of high-speed trains head type is conducted.Two head shapes with three marshalling 1∶8 scaled were modeled for aerodynamic noise simulation,and the average total sound pressure level of the far field measured points was calculated.By comparison with the results of wind tunnel tests, the two values differ by less than 3 dBA,and that of the head shape 1 is less.In order to achieve aerodynamic noise numerical simulation of a full-size high-speed train head,the subdomain method is proposed,and part of the computational domain of the three marshalling is intercepted.The aerodynamic characteristics obtained from the subdomain and the whole train are consistent,which shows the rationality of the subdomain method indirectly.The subdomain method is used on the full-size head shapes 1 and 2 to perform aerodynamic noise simulation,and surface acoustic power of the head,total sound pressure level of the surface and the field were obtained,which can provide the basis for selection and optimization of the head.An aerodynamic noise analysis method is established based on numerical simulation of the full size high-speed train and solved the problem that wind tunnel tests and numerical simulation analysis could not be carried out on full-size high-speed trains.

high-speed train；head shape；aerodynamic noise；subdomain

1673- 9590(2017)03- 0030- 06

2016- 10- 12基金項目:國家重點研發(fā)計劃資助項目(2016YFB1200503)

高陽(1980-),男,高級工程師, 碩士,從事高速列車噪聲與振動控制的研究E-mail:ck_gy@126.com.

A