袁志群，吉斯?jié)h，谷正氣，劉顯貴，3，許西安

前言

隨著車速的提高，汽車空氣動(dòng)力學(xué)性能備受重視，風(fēng)洞試驗(yàn)是評(píng)估汽車空氣動(dòng)力學(xué)性能的重要手段。因此，世界各大汽車公司都不惜投入巨資建設(shè)專有的汽車風(fēng)洞，國(guó)內(nèi)外高校也投入大量資金建設(shè)和改造現(xiàn)有的風(fēng)洞，在高速試驗(yàn)段加裝天平測(cè)力系統(tǒng)和邊界層控制系統(tǒng)，借助先進(jìn)的設(shè)備和科學(xué)的方法開(kāi)展汽車空氣動(dòng)力學(xué)性能研究。

從現(xiàn)有研究可知，目前國(guó)內(nèi)的高校都投入了大量資金購(gòu)買先進(jìn)的邊界層控制系統(tǒng)和模型支撐系統(tǒng)，用來(lái)建設(shè)汽車風(fēng)洞高速試驗(yàn)段。吉林大學(xué)汽車風(fēng)洞試驗(yàn)室采用移動(dòng)帶地面效應(yīng)模擬系統(tǒng)[1]，其洞體結(jié)構(gòu)為回流開(kāi)式[2]，帶有寬帶式移動(dòng)地板和前部垂直抽吸系統(tǒng)，試驗(yàn)段氣流穩(wěn)定，速度梯度??；湖南大學(xué)汽車風(fēng)洞試驗(yàn)室采用抽吸固定地板的方式消除地板邊界層的影響，根據(jù)抽吸參數(shù)變化對(duì)流場(chǎng)的影響，對(duì)抽吸參數(shù)進(jìn)行了優(yōu)化，提高了試驗(yàn)精度[3]，試驗(yàn)?zāi)Ｐ筒捎盟妮喼蔚男问脚c六分力浮框式應(yīng)變天平連接[4]，設(shè)計(jì)制造了一套固定鎖緊裝置，防止模型在試驗(yàn)過(guò)程移動(dòng)，影響測(cè)量結(jié)果[5]；同濟(jì)大學(xué)上海地面交通工具風(fēng)洞中心是國(guó)內(nèi)唯一的全尺寸汽車風(fēng)洞試驗(yàn)室，試驗(yàn)段也是回流開(kāi)口式，地面效應(yīng)采用五帶移動(dòng)路面系統(tǒng)和轉(zhuǎn)盤(pán)，試驗(yàn)?zāi)Ｐ团c天平轉(zhuǎn)盤(pán)采用輪下連接的固定方式[6]，采用固定模式的錨定方式，試驗(yàn)方案更安全，測(cè)量結(jié)果更精確[7]。試驗(yàn)段噴口處加裝了渦流發(fā)生器，使流場(chǎng)品質(zhì)顯著提高[8]。

廈門理工學(xué)院風(fēng)工程試驗(yàn)室于2017年竣工驗(yàn)收，是福建省首座風(fēng)洞試驗(yàn)室，可用于汽車模型風(fēng)洞試驗(yàn)，本文中采用計(jì)算機(jī)輔助設(shè)計(jì)與分析的方法，在該風(fēng)洞高速試驗(yàn)段設(shè)計(jì)了一種經(jīng)濟(jì)簡(jiǎn)便、滿足工程精度要求的汽車模型風(fēng)洞試驗(yàn)臺(tái)架，建立了某汽車風(fēng)洞試驗(yàn)的數(shù)值分析模型，根據(jù)該模型的風(fēng)洞試驗(yàn)結(jié)果，分析了影響試驗(yàn)精度的關(guān)鍵參數(shù)，并對(duì)試驗(yàn)臺(tái)架的干擾進(jìn)行了分析與修正。

1 試驗(yàn)臺(tái)架結(jié)構(gòu)與原理

試驗(yàn)臺(tái)架位于風(fēng)洞高速試驗(yàn)段，由底座、地板、動(dòng)力系統(tǒng)和試驗(yàn)?zāi)Ｐ?部分組成，如圖1所示。

圖1 試驗(yàn)臺(tái)架結(jié)構(gòu)

底座與風(fēng)洞地板采用螺栓連接，動(dòng)力系統(tǒng)安裝在底座的中間立柱上，四周的6根立柱起輔助支撐作用；地板位于底座之上，抬高的地板可以消除地面邊界層的影響?；啺惭b在底座四周的6根立柱上，地板可繞底座的中心軸旋轉(zhuǎn)，以模擬不同橫擺角的工況；動(dòng)力系統(tǒng)動(dòng)力源為交流伺服電機(jī)，通過(guò)轉(zhuǎn)化電信號(hào)輸出動(dòng)力和運(yùn)動(dòng)，可實(shí)現(xiàn)精準(zhǔn)的速度和轉(zhuǎn)角控制，滿足設(shè)計(jì)功能要求。動(dòng)力通過(guò)立式電機(jī)輸出，經(jīng)過(guò)1級(jí)斜齒輪減速后傳遞給地板。試驗(yàn)?zāi)Ｐ凸潭ㄔ诘匕迳?，氣?dòng)力天平安裝在模型之內(nèi)，實(shí)現(xiàn)氣動(dòng)力測(cè)量。

2 試驗(yàn)臺(tái)架數(shù)學(xué)模型的建立

為驗(yàn)證臺(tái)架結(jié)構(gòu)的試驗(yàn)精度，利用計(jì)算流體動(dòng)力學(xué)方法建立數(shù)學(xué)分析模型。在三維建模軟件Siemens NX 10.0中，構(gòu)建試驗(yàn)臺(tái)架的幾何模型，如圖1所示，保留主體結(jié)構(gòu)，并對(duì)局部細(xì)節(jié)進(jìn)行適當(dāng)簡(jiǎn)化；在前處理軟件ICEM-CFD中，采用OCTREE方法對(duì)求解域進(jìn)行離散，劃分5層三棱柱網(wǎng)格精確捕捉車身表面附面層，在流動(dòng)分離區(qū)域?qū)W(wǎng)格進(jìn)行細(xì)化，共劃分非結(jié)構(gòu)化網(wǎng)格540萬(wàn)個(gè)，計(jì)算模型網(wǎng)格如圖2所示，定義該計(jì)算模型為模型1。

圖2 模型1計(jì)算網(wǎng)格

汽車屬非流線型物體，車身表面曲率變化大，流動(dòng)分離劇烈，屬于三維不可壓、等溫、定常流場(chǎng)[9-10]。求解計(jì)算在大型商業(yè)軟件Fluent中完成，根據(jù)現(xiàn)有研究現(xiàn)狀可知，Realizable k-e湍流模型廣泛應(yīng)用于汽車外流場(chǎng)分析，其計(jì)算精度比其它湍流模型更高[11-13]。車身附面層采用標(biāo)準(zhǔn)壁面函數(shù)求解，空間離散格式采用2階迎風(fēng)格式，計(jì)算方法運(yùn)用SIMPLE算法，側(cè)風(fēng)的計(jì)算方法參見(jiàn)文獻(xiàn)[14]和文獻(xiàn)[15]，其它邊界設(shè)置見(jiàn)表1。

為驗(yàn)證數(shù)學(xué)分析模型的準(zhǔn)確性和試驗(yàn)臺(tái)架支架對(duì)汽車流場(chǎng)的干擾影響，另建立汽車外流場(chǎng)分析的常規(guī)模型，建模方法詳見(jiàn)文獻(xiàn)[14]，定義該計(jì)算模型為模型2。為使之與模型1和風(fēng)洞試驗(yàn)的地板邊界層數(shù)據(jù)接近，使計(jì)算結(jié)果具有可比性，將地板一分為二，如圖3所示。模型2網(wǎng)格劃分和其它邊界條件設(shè)置與模型1完全相同。模型2風(fēng)洞試驗(yàn)在湖南大學(xué)HD-2風(fēng)洞完成，如圖4所示，流場(chǎng)測(cè)量采用粒子圖像測(cè)速技術(shù)(PIV)。

表1 邊界條件設(shè)置

圖3 模型2地板邊界條件

圖4 模型2風(fēng)洞試驗(yàn)

3 結(jié)果分析與討論

通過(guò)大量的數(shù)值分析得知，臺(tái)架地板直徑、臺(tái)架高度和臺(tái)架形狀是影響模型氣動(dòng)特性數(shù)據(jù)變化的主要因素，其中臺(tái)架地板直徑影響最大，它會(huì)影響地板邊界層的速度分布和汽車周圍流場(chǎng)，氣動(dòng)阻力和氣動(dòng)升力會(huì)隨之改變，試驗(yàn)結(jié)果如圖5～圖12所示，圖中L為車長(zhǎng)。

圖5為模型1在橫擺角為0°、風(fēng)速為30m/s工況下氣動(dòng)力系數(shù)隨地板直徑變化的關(guān)系曲線。由圖5可見(jiàn)，隨著地板直徑的增加，氣動(dòng)阻力系數(shù)先減小后稍有增加，而氣動(dòng)升力系數(shù)則一直呈增加趨勢(shì)，但變化越來(lái)越緩慢。地板直徑較小時(shí)，地板前緣的氣流分離和再附著、地板后緣的氣流分離對(duì)汽車近地面流場(chǎng)產(chǎn)生嚴(yán)重干擾。地板前緣的再附著氣流使接近汽車前端的氣流速度方向發(fā)生改變，不再是水平來(lái)流方向，而是與水平線成一定夾角，如圖6中箭頭所示。地板直徑越小，夾角越大。夾角增大時(shí)，X軸方向氣流速度減小，氣動(dòng)阻力減小，Z軸方向氣流速度增加，氣動(dòng)升力增加。圖7和圖8分別為模型1縱對(duì)稱面的車身表面壓力云圖和車底速度。由圖可見(jiàn)，當(dāng)?shù)匕逯睆捷^大時(shí)，汽車底部邊界層增厚，車底氣流阻塞效應(yīng)增強(qiáng)，車底流速減小，壓力升高，因此氣動(dòng)阻力和氣動(dòng)升力都會(huì)增加。

圖5 氣動(dòng)力系數(shù)隨地板直徑變化關(guān)系

圖6 模型1縱對(duì)稱面流線圖

由圖8可見(jiàn):當(dāng)來(lái)流遇到臺(tái)架地板時(shí)，由于地板前緣的影響，氣流在此處加速，并改變運(yùn)動(dòng)方向，速度會(huì)出現(xiàn)第一個(gè)峰值；當(dāng)氣流遇到車頭時(shí)，由于車頭的阻滯作用，速度急劇下降，在車頭位置降到最低，然后氣流在車頭下緣加速，出現(xiàn)另一個(gè)速度峰值；由于汽車底部的阻塞效應(yīng)和黏性的影響，車底氣流速度逐漸減小。地板直徑越大，車身底部同樣位置的邊界層越厚，速度越低，當(dāng)?shù)匕逯睆綇?.8L增加到2.8L時(shí)，車底平均速度降低了6.2%。

圖7 模型1縱對(duì)稱面和車身壓力云圖

圖8 模型1縱對(duì)稱面車底速度分布

圖9 模型1尾部湍流動(dòng)能云圖

圖9 為不同地板直徑下距離車尾0.55m處湍流動(dòng)能分布云圖。由圖可見(jiàn):當(dāng)?shù)匕逯睆捷^小時(shí)，臺(tái)架支架的尾流將會(huì)嚴(yán)重影響汽車尾部氣流，汽車尾渦與臺(tái)架尾渦相互作用。從圖6也可明顯看出車尾氣流的相互干涉，這些現(xiàn)象都會(huì)改變汽車尾部流場(chǎng)，進(jìn)而導(dǎo)致氣動(dòng)力的變化。

通過(guò)以上分析可知，地板直徑大小的變化會(huì)影響汽車前端氣流方向、車底流速分布和尾渦形態(tài)。結(jié)合圖5數(shù)據(jù)分析可知，存在一個(gè)合理的地板直徑，風(fēng)洞臺(tái)架對(duì)模型周圍流場(chǎng)干擾最小，汽車周圍流場(chǎng)與實(shí)際情況相符，氣動(dòng)力系數(shù)接近實(shí)際值。以模型2風(fēng)洞試驗(yàn)數(shù)據(jù)為參照，當(dāng)?shù)匕逯睆郊s為2.45m(2.6L)時(shí)，氣動(dòng)力系數(shù)計(jì)算值與風(fēng)洞試驗(yàn)值偏差最小，如表2所示。

表2 氣動(dòng)力系數(shù)結(jié)果對(duì)比

誤差產(chǎn)生的原因主要包括兩方面:一方面是由于數(shù)值計(jì)算本身帶來(lái)的誤差，包括模型一致性、網(wǎng)格劃分和邊界條件設(shè)置等，這些誤差不可避免，根據(jù)以往的工程經(jīng)驗(yàn)可知，該部分誤差可控制在5%以內(nèi)[14]；另一方面是臺(tái)架支架的干擾，雖然某一尺寸的地板直徑能將誤差控制到最低，但或多或少都會(huì)影響汽車近地面的流場(chǎng)。為研究支架干擾產(chǎn)生的誤差量，對(duì)比了模型1和模型2數(shù)值計(jì)算結(jié)果，如圖10和圖11所示。

圖10 氣動(dòng)力系數(shù)隨風(fēng)速變化關(guān)系

圖11 氣動(dòng)力系數(shù)隨橫擺角變化關(guān)系

圖10 為橫擺角為0°時(shí)，模型1和模型2的氣動(dòng)力系數(shù)隨風(fēng)速而變化的關(guān)系，模型1臺(tái)架地板直徑取2.6L。由圖10可見(jiàn)，模型1和模型2氣動(dòng)力系數(shù)變化趨勢(shì)相同，且都比較穩(wěn)定。隨著風(fēng)速的升高，氣動(dòng)力系數(shù)稍有降低。只是由于支架對(duì)流場(chǎng)的干擾影響，模型1氣動(dòng)力系數(shù)均高于模型2。通過(guò)對(duì)比模型2的計(jì)算分析結(jié)果與風(fēng)洞試驗(yàn)結(jié)果可知，兩者比較接近，以試驗(yàn)值為基準(zhǔn)，氣動(dòng)阻力系數(shù)平均計(jì)算誤差為3.1%，氣動(dòng)升力系數(shù)平均計(jì)算誤差為4.8%，證明了本文數(shù)值計(jì)算方法可行，結(jié)果可信。

圖11為風(fēng)速為30m/s時(shí)模型1和模型2的氣動(dòng)力系數(shù)隨橫擺角而變化的關(guān)系。由圖可見(jiàn)，氣動(dòng)力系數(shù)隨著橫擺角的增加而增加，氣動(dòng)側(cè)向力系數(shù)增加明顯，其次為氣動(dòng)升力，而氣動(dòng)阻力變化比較平緩，與前人的研究結(jié)論吻合[16-17]。另外，由于支架的影響，模型1的氣動(dòng)力系數(shù)均高于模型2，且差值隨著橫擺角的增加而增大。當(dāng)橫擺角為15°時(shí)，支架的干擾使氣動(dòng)阻力系數(shù)增大了0.014，即4.2%，使氣動(dòng)側(cè)力增大了0.029，即3.9%，使氣動(dòng)升力增大了0.024 6，即9.7%。從百分比看，支架干擾對(duì)氣動(dòng)升力影響較大，對(duì)氣動(dòng)阻力和氣動(dòng)側(cè)力的影響較小。這也與上面分析的結(jié)論吻合，因?yàn)榕_(tái)架的干擾會(huì)影響汽車前端流場(chǎng)和邊界層分布，是引起升力變化的主要原因。通過(guò)對(duì)比分析模型1和模型2的氣動(dòng)力數(shù)據(jù)，可為臺(tái)架的干擾修正提供理論依據(jù)。

雖然臺(tái)架會(huì)對(duì)流場(chǎng)產(chǎn)生干擾，使氣動(dòng)力數(shù)據(jù)不同程度增加，但合適的臺(tái)架地板直徑可使干擾量降到最低，氣動(dòng)阻力系數(shù)和氣動(dòng)升力系數(shù)誤差在8%以內(nèi)，說(shuō)明該風(fēng)洞臺(tái)架的氣動(dòng)力試驗(yàn)精度較好，能滿足工程實(shí)際的要求。為了進(jìn)一步驗(yàn)證風(fēng)洞臺(tái)架的試驗(yàn)精度，接下來(lái)將從流場(chǎng)測(cè)量數(shù)據(jù)進(jìn)行對(duì)比，如圖12～圖14所示。

圖12 縱對(duì)稱面速度云圖和流線圖(β＝0°)

圖12 為模型1和模型2縱對(duì)稱面速度云圖和流線圖。由圖可見(jiàn)，車尾存在兩個(gè)旋向相反的橢圓形尾渦，順時(shí)針渦A是由尾部行李箱的上緣氣流拖曳形成，逆時(shí)針渦B是由車底部氣流翻卷至尾流中形成，且渦A的尺寸比渦B大，占據(jù)尾渦空間的主要部分。兩個(gè)模型的氣動(dòng)數(shù)據(jù)計(jì)算結(jié)果比較接近，車頭和車尾速度分布與氣流流向、尾渦形態(tài)與位置都非常相似，但在車底氣流速度分布存在一定的差別。對(duì)模型2車尾主渦A進(jìn)行PIV測(cè)量，結(jié)果如圖13所示。由圖可見(jiàn)，主渦A的形態(tài)、渦核的位置計(jì)算準(zhǔn)確，局部區(qū)域速度大小存在一定差別，進(jìn)一步證明了風(fēng)洞臺(tái)架流場(chǎng)干擾小，試驗(yàn)精度高。

圖13 模型2車尾主渦A風(fēng)洞試驗(yàn)測(cè)量值(β＝0°)

圖14 給出了模型1和模型2在橫擺角為0°、風(fēng)速為30m/s時(shí)，縱對(duì)稱面車身上表面壓力系數(shù)分布。從圖可見(jiàn)，壓力系數(shù)大小、分布形態(tài)基本相同，說(shuō)明臺(tái)架對(duì)車身上表面流場(chǎng)干擾較小。

圖15為模型1和模型2在橫擺角為15°、風(fēng)速為30m/s時(shí)車身迎風(fēng)側(cè)壓力云圖。從圖可見(jiàn)，壓力大小、正壓和負(fù)壓分布區(qū)域及壓力變化基本一致。

圖14 縱對(duì)稱面車身上表面壓力系數(shù)分布(β＝0°)

圖15 車身迎風(fēng)側(cè)壓力云圖(β＝15°)

通過(guò)對(duì)比模型1和模型2的各主要?dú)鈩?dòng)特性數(shù)據(jù)可知，當(dāng)試驗(yàn)臺(tái)架地板具有合適的直徑時(shí)，可將支架干擾量降到最低，汽車模型氣動(dòng)力試驗(yàn)誤差最小，流場(chǎng)測(cè)量數(shù)據(jù)與風(fēng)洞試驗(yàn)結(jié)果接近。

4 結(jié)論

本文中設(shè)計(jì)了一種簡(jiǎn)易的汽車模型風(fēng)洞試驗(yàn)臺(tái)架，利用計(jì)算流體動(dòng)力學(xué)方法建立了其數(shù)學(xué)模型，總結(jié)了影響該裝置試驗(yàn)精度的關(guān)鍵參數(shù)，分析了其影響規(guī)律，得到了該模型在不同風(fēng)速、不同橫擺角下的氣動(dòng)力數(shù)據(jù)，并與常規(guī)數(shù)學(xué)模型建模方法得到的數(shù)據(jù)和風(fēng)洞試驗(yàn)數(shù)據(jù)進(jìn)行了對(duì)比，得出以下結(jié)論。

(1)臺(tái)架地板直徑是影響試驗(yàn)精度的重要參數(shù)。氣動(dòng)阻力隨著直徑的增加先減小后增加，氣動(dòng)升力隨著直徑的增加而增大；地板直徑較小時(shí)，主要影響車身前端氣流和尾渦形態(tài)，地板直徑較大時(shí)，主要影響地板附面層。

(2)試驗(yàn)臺(tái)架支架干擾不可忽略。對(duì)氣動(dòng)升力影響較大，對(duì)氣動(dòng)阻力和氣動(dòng)側(cè)力的影響較小。當(dāng)橫擺角為15°時(shí)，臺(tái)架支架的干擾使氣動(dòng)阻力增大4.2%，氣動(dòng)側(cè)力增大 3.9%，而使氣動(dòng)升力增大9.7%。

(3)臺(tái)架裝置的試驗(yàn)精度較好。當(dāng)?shù)匕逯睆綖?.6L時(shí)，支架干擾量最小，氣動(dòng)阻力計(jì)算精度為6.7%，氣動(dòng)升力計(jì)算精度為7.5%，且車身流場(chǎng)與風(fēng)洞試驗(yàn)結(jié)果接近，試驗(yàn)精度達(dá)到工程要求。

[1] 張英朝，詹大鵬，趙婧，等.汽車風(fēng)洞地面效應(yīng)模擬系統(tǒng)影響的研究[J].汽車工程，2016，38(12):1434-1439.

[2] 張英朝，張喆，王婧宇，等.汽車模型表面壓力系數(shù)的測(cè)試與數(shù)值模擬[J].汽車工程，2011，33(11):951-954.

[3] 楊易，蔡圣康，徐永康，等.抽吸系統(tǒng)對(duì)汽車風(fēng)洞試驗(yàn)的影響及優(yōu)化[J].中國(guó)公路學(xué)報(bào)，2015，28(2):119-126.

[4] 張勇，谷正氣，劉長(zhǎng)水，等.基于風(fēng)洞試驗(yàn)與CFD分析的汽車前后軸氣動(dòng)升力計(jì)算的研究[J].汽車工程，2017，39(3):275-280.

[5] 張勇，谷正氣，劉長(zhǎng)水.車身姿勢(shì)對(duì)風(fēng)洞試驗(yàn)氣動(dòng)升力測(cè)量影響研究[J].汽車工程，2015，37(3):295-299.

[6] 周華，朱暉，楊志剛.整車風(fēng)洞力/壓同測(cè)實(shí)驗(yàn)方法研究[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2016，44(4):614-619.

[7] 陳力，劉曉輝，龐加斌，等.風(fēng)洞試驗(yàn)中車輛錨定方式對(duì)氣動(dòng)力測(cè)量的影響分析[J].汽車工程，2014，36(11):1369-1371.

[8] 賈青，吳桐，楊志剛.渦流發(fā)生器對(duì)3/4汽車模型風(fēng)洞流場(chǎng)品質(zhì)的影響[J].同濟(jì)大學(xué)學(xué)報(bào)(自然科學(xué)版)，2014，42(6):956-962.

[9] KIM Jeong Jae， LEE Sangseung， KIM Myeongkyun， et al.Salient drag reduction of a heavy vehicle using modified cab-roof fairings[J].Journal of Wind Engineering and Industrial Aerodynamics，2017，164:138-151.

[10] TIENPHUC D，ZHENGQI G，ZHEN C.Numerical simulation of the flow filed around generic formula one[J].Journal of Applied Fluid Mechanics，2016，9(1):443-450.

[11] 楊易，鄭萌，黃劍鋒，等.基于非光滑表面與渦流干擾的車身氣動(dòng)減阻方法[J].中國(guó)機(jī)械工程，2016，27(7):982-988.

[12] 谷正氣，陳陣，黃泰明，等.基于改進(jìn)LRN k-e模型的汽車氣動(dòng)特性研究[J].中國(guó)機(jī)械工程，2015，26(18):2550-2555.

[13] 汪怡平，王濤，黎帥.基于自由變形技術(shù)的汽車氣動(dòng)減阻優(yōu)化[J].機(jī)械工程學(xué)報(bào)，2017，53(9):135-143.

[14] 袁志群，谷正氣，楊明智，等.底部結(jié)構(gòu)對(duì)轎車側(cè)風(fēng)氣動(dòng)特性的影響分析[J].汽車工程，2017，39(1):28-34.

[15] 黃泰明，谷正氣，豐成杰.不同波長(zhǎng)瞬態(tài)側(cè)風(fēng)對(duì)汽車氣動(dòng)性能影響分析[J].中國(guó)機(jī)械工程，2016，27(22):3111-3117.

[16] KIM Se-Jin，YOO Chul-Hwan， KIM Ho-Kyung.Vulnerability assessment for the hazards of crosswinds when vehicles cross a bridge deck[J].Journal of Wind Engineering and Industrial Aerodynamics，2016，156:62-71.

[17] TSUBOKURA M，NAKASHIMA T，KITAYAMA M，et al.Large eddy simulation on the unsteady aerodynamic response of a road vehicle in transient crosswinds[J].International Journal of Heat and Fluid Flow，2010，31:1075-1086.