某轎車車身外流場空氣動(dòng)力學(xué)研究

胡樹山，劉榮娥

（鄂爾多斯應(yīng)用技術(shù)學(xué)院，內(nèi)蒙古 鄂爾多斯 017000）

引言

汽車空氣動(dòng)力性能是汽車重要特性之一，它直接影響了汽車的動(dòng)力性能、經(jīng)濟(jì)型、操縱穩(wěn)定性和冷卻通風(fēng)、風(fēng)噪聲以及安全性[1]。因此，汽車制造行業(yè)一直致力于研究降低汽車的空氣阻力系數(shù) Cd值，以降低燃油消耗率，從而達(dá)到節(jié)省燃油的目的[2]。

降低空氣阻力系數(shù)對(duì)改善燃油經(jīng)濟(jì)性極其可觀，轎車氣動(dòng)阻力系數(shù)由0.42降為0.3時(shí)，混合循環(huán)工況下燃油經(jīng)濟(jì)性可提高約9%；而當(dāng)車速達(dá)到150km／h時(shí)，燃油經(jīng)濟(jì)性改善可達(dá) 25%之多[3]。為了減少汽車燃料的消耗，現(xiàn)代汽車的外形設(shè)計(jì)越來越趨于流線型，使得轎車的風(fēng)阻系數(shù)越來越小[4]。

1 模型建立

本文通過某款轎車四視圖及外形尺寸，利用 CATIA V5 R20中影像草圖功能模塊，導(dǎo)入視圖，再通過自由曲面和創(chuàng)成式外形設(shè)計(jì)功能模塊進(jìn)行正向開發(fā)。將上述視圖文件中的主要外形尺寸進(jìn)行整理，受計(jì)算機(jī)CPU運(yùn)算能力的限制，按1：5比例計(jì)算得到模型外形參數(shù)尺寸，如表1所示。

對(duì)實(shí)車做了如下簡化：忽略車身外部突起物如后視鏡、刮雨器等部分；將輪胎簡化為無輻板式實(shí)體輪胎；且對(duì)車身底部做了簡化，如圖1所示。

表1 模型參數(shù)

圖1 轎車整車建模

CATIA V5 R20軟件提供了豐富的曲面分析工具，能夠精確控制曲面質(zhì)量，在建模過程中，所有曲線都須做質(zhì)量分析，特別是給車身覆蓋件做鋪面的曲線，必須達(dá)到要求。本模型車身表面的曲面都達(dá)到A級(jí)曲面的標(biāo)準(zhǔn)。

通過車身光順性分析[5]可知，從每個(gè)曲面來看，走勢順暢、非常連續(xù)，沒有出現(xiàn)突變的部位；過渡曲面與其相鄰的曲面過渡流暢，沒有間斷，較為美觀；從整體效果來看，高亮線排列比較均勻且都是連續(xù)的，紋理過渡平滑自然，曲面光滑平整，造型整體效果較好，整車的表面質(zhì)量很高。

圖2 等照度線映射分析

圖3 環(huán)境對(duì)映分析

2 控制方程

2.1 質(zhì)量守恒方程

任何流動(dòng)問題都必須滿足質(zhì)量守恒定律，又稱連續(xù)方程。按這一規(guī)律質(zhì)量守恒方程如下：

式中：ρ為密度；t為時(shí)間，u、v和w分別為速度矢量在x、y和z方向的分量。若流體不可壓，ρ為常數(shù)，則方程變?yōu)椋?/p>

2.2 動(dòng)量守恒方程

動(dòng)量守恒方程也是任何流動(dòng)問題都必須滿足的守恒定律，即牛頓第二定律。按這一規(guī)律，在x、y和z方向上的動(dòng)量守恒方程如下：

式中：p為流體微元上的壓力；τxx、τxy和τxz等分別為因分子黏性作用在微元體上的黏性應(yīng)力τ的分量；Fx、Fy和Fz為微元體上的體力。

2.3 能量守恒方程

能量守恒定律是含熱交換流動(dòng)系統(tǒng)必須滿足的基本定律，即熱力學(xué)第一定律。方程如下：

式中：cp為比熱容；T為溫度；k為流體的傳熱系數(shù)；ST為流體的內(nèi)熱源及由于黏性作用的流體機(jī)械能轉(zhuǎn)化為熱能的部分，有時(shí)簡稱ST為粘性耗散項(xiàng)。

對(duì)于不可壓縮流體，若熱交換量很小甚至可忽略，可以不考慮能量守恒方程，這樣，只需聯(lián)立求解連續(xù)方程和運(yùn)動(dòng)方程即可。

2.4 湍流控制方程

一般認(rèn)為，無論湍流運(yùn)動(dòng)多么復(fù)雜，非穩(wěn)態(tài)的連續(xù)方程和動(dòng)量方程對(duì)于湍流的瞬時(shí)運(yùn)動(dòng)仍然是適用的。因此，不考慮壓縮流動(dòng)，使用笛卡爾坐標(biāo)，速度矢量在x、y、z上的分量為u、v和w，湍流瞬時(shí)控制方程如下：

3 車身外流場數(shù)值模擬計(jì)算

3.1 計(jì)算空間建立

計(jì)算域尺寸的確定其實(shí)就是確定“模擬風(fēng)洞”的大小。其在建立計(jì)算域尺寸大小的時(shí)候最主要的考慮就是阻塞比φ對(duì)汽車氣動(dòng)特性的影響，阻塞比φ的計(jì)算公式如下：

式中：Aref表示轎車迎風(fēng)面積；WC表示計(jì)算域的寬度；HC表示計(jì)算域的高度。

關(guān)于計(jì)算域尺寸的確定，參考了相關(guān)文獻(xiàn)，設(shè)車身長度為L，車身寬度為W，車總高為H，為了盡量解除數(shù)值風(fēng)洞的有限體積對(duì)車身外流場數(shù)值計(jì)算的影響，計(jì)算域設(shè)定車前部為3L，車后部為6L，長方體高為6H，寬為6W。具體計(jì)算域尺寸參數(shù)設(shè)置，見表2。

此方案的阻塞比φ大小為2.3%，滿足SAEJ1252推薦的模型風(fēng)洞的阻塞比小于5%的要求[6]。

表2 計(jì)算域相關(guān)尺寸數(shù)據(jù)

圖4 CATIA計(jì)算域半剖模型

3.2 網(wǎng)格劃分

圖5 計(jì)算域的網(wǎng)格劃分

圖6 對(duì)局部加密處理

圖7 車身表面的網(wǎng)格處理

靠近車身以及車身與地面之間的計(jì)算區(qū)域由于流動(dòng)比較復(fù)雜，也是研究的重點(diǎn)區(qū)域，采用較密的網(wǎng)格劃分，這樣有利于流場的分析。而在遠(yuǎn)離車身的計(jì)算區(qū)域，希望網(wǎng)格生成盡量簡便，節(jié)省機(jī)時(shí)，采用較稀疏的網(wǎng)格。最后，網(wǎng)格總數(shù)為1209433，節(jié)點(diǎn)數(shù)為218887。

3.3 軟件參數(shù)設(shè)置

1）在本文計(jì)算當(dāng)中，邊界條件設(shè)置可見圖8，具體參數(shù)如下：

（1）入口邊界：取車身前遠(yuǎn)端面為入口邊界，計(jì)算時(shí)給定入口速度條件為25m/s。

（2）出口邊界：取車身后遠(yuǎn)端面為出口邊界，邊界條件設(shè)置壓力出口。

（3）固定壁面條件：將車輛模型外表面設(shè)置為固定壁面，滿足固定壁面無滑移定律。

（4）由于計(jì)算機(jī)配置有限，無法進(jìn)行過于復(fù)雜的計(jì)算，故不考慮地面效應(yīng)。

圖8 邊界條件設(shè)置

2）設(shè)定介質(zhì)物理參數(shù)

因?yàn)槠嚨男旭偹俣韧ǔ?-300km/h之間，流體介質(zhì)只為空氣。馬赫數(shù)對(duì)于普通轎車來講遠(yuǎn)小于 1，因此可以認(rèn)為空氣是不可壓縮的，且不考慮空氣物理參數(shù)值的變化，即空氣的溫度、粘性和壓強(qiáng)不變[7]，如表3所示。

表3 空氣物理參數(shù)

本文選擇標(biāo)準(zhǔn)k-ε模型作為計(jì)算的湍流模型，迭代方式采用SIMPLEC算法，離散格式為二階迎風(fēng)格式[8]。

4 結(jié)果分析

4.1 氣動(dòng)阻力

根據(jù)建模的坐標(biāo)系方向，y方向?yàn)闅鈩?dòng)阻力方向，選擇y方向（0，1，0），顯示出車身模型在數(shù)值模擬計(jì)算中所受到的氣動(dòng)阻力值，如圖 9。同理，選擇z方向（0，0，1），即可顯示車身模型在外流場模擬中所受的的氣動(dòng)升力值，車身模型在外流場模擬中所受的氣動(dòng)升力值，如圖10。

圖9 計(jì)算結(jié)果（一）

圖10 計(jì)算結(jié)果（二）

4.2 流場圖形分析

（1）速度云圖

根據(jù)圖中的箭頭疏密程度和顏色可以看出，氣流從車頭到車尾的空氣流動(dòng)形態(tài)。氣流在撞擊車頭后分為兩股，一股順車身模型的引擎蓋流向擋風(fēng)玻璃，另外一股向下沿車身從底部流走。

圖11 對(duì)稱面速度矢量分布

在車身上部分的氣流在引擎蓋端緣開始加速，由于受到前風(fēng)窗的阻擋，氣流迅速減緩，并形成一個(gè)小的阻滯區(qū)，然后速度開始增加，在前風(fēng)窗與車頂相接處達(dá)到最大值。之后，在車身上面的氣流速度略有減小，但是數(shù)值仍然較大。當(dāng)氣流到達(dá)車頂與后風(fēng)窗相接處的時(shí)候，氣流速度減小，并發(fā)生分離現(xiàn)象，在車身尾部形成了一個(gè)渦流區(qū)。

在車身下部分的氣流則從車身底部和路面之間強(qiáng)制通過，氣流并未在車身尾部直接流走，而是沿汽車尾部與車身上表面的氣流匯合，并形成渦流區(qū)域，因此產(chǎn)生很大的阻力。

圖12 對(duì)稱面上速率分布圖

圖12是對(duì)稱面速率分布圖，從圖中可以看出，在尾部有負(fù)的速率分布區(qū)，圖中尾部藍(lán)色區(qū)域，即回流區(qū)。由于擾動(dòng)過大，車身底部表面邊界層很快發(fā)展成湍流邊界，且迅速增大甚至充滿整個(gè)車身模型與地面間的空間。從圖中還可以明顯看見車身表面氣流速度為零，這是附著在車身表面的一層邊界層，而且從車身上部流出的氣體速度明顯高于車身底部流出的氣體速度，這兩種氣流速度的不同，導(dǎo)致了車身產(chǎn)生了向上的升力。

（2）壓力云圖

根據(jù)兩個(gè)云圖可以看出，在汽車頭部及四個(gè)輪胎的迎風(fēng)面呈深紅色，可見氣流受到垂直方向的阻值降為零，全部動(dòng)壓轉(zhuǎn)換為靜壓。流向車身上方的這部分氣流在流過前風(fēng)窗下邊緣時(shí)，受到阻力，轉(zhuǎn)化為靜壓，可見前風(fēng)窗角影響重大。所以在車身設(shè)計(jì)時(shí)，前風(fēng)窗應(yīng)盡量平躺，采用大傾角設(shè)計(jì)。

圖13 對(duì)稱面上的壓力分布圖

圖14 車身表面靜壓力分布圖

在有一定斜度引擎蓋上流過的氣流速度仍然較快，所以在引擎蓋起緣處壓力仍為負(fù)值，導(dǎo)致該處負(fù)壓區(qū)的產(chǎn)生。當(dāng)氣流到達(dá)引擎蓋和前擋風(fēng)玻璃相交處，擋風(fēng)玻璃的存在致使氣流速度降低，形成一個(gè)滯區(qū)。而且在滯區(qū)中有內(nèi)部渦流的產(chǎn)生，所以該區(qū)具有正壓力。后風(fēng)窗處由于頂部的氣流受到行李廂的阻礙，所以受到阻力增大，形成局部正壓區(qū)。

（3）頭部氣流

由矢量圖可以看出，氣流在受到車頭阻擋后，速度下降很快，幾乎到接近于零，并導(dǎo)致阻滯區(qū)的形成。速度急速下降后的氣流在此處分為兩股，一股經(jīng)引擎蓋向上流動(dòng)，并在前風(fēng)窗與引擎蓋的相交處發(fā)生分離，阻滯區(qū)再次形成。由于上翹角影響到車身底部與車頭前端的過渡區(qū)域的氣流，致使車身底部氣流通過速度較小。

圖15 車頭對(duì)稱面速度分布圖

（4）尾部氣流

在車身尾部存在有兩個(gè)旋渦，在車身上面形成的大旋渦按順時(shí)針方向旋轉(zhuǎn)，在車身下方形成的小旋渦按逆時(shí)針方向旋轉(zhuǎn)。這是因?yàn)楫?dāng)遠(yuǎn)方的氣流流經(jīng)整個(gè)車身后，兩股氣流中的大部分氣流會(huì)在身后很遠(yuǎn)處交匯成一股氣流，所以在交匯之前有一段很長的負(fù)壓區(qū)。受負(fù)壓區(qū)的壓力影響，從車身頂部和車身底部流出氣流中的各一部分氣，會(huì)被吸進(jìn)負(fù)壓區(qū)。從而產(chǎn)生了這兩個(gè)按不同方向旋轉(zhuǎn)的旋渦，這也是氣動(dòng)阻力的最主要來源。

圖16 車尾對(duì)稱面速度分布圖

由圖17可以看出，湍流動(dòng)能基本集中于汽車的尾部分離區(qū)域，雖然流速較低，但由于受到分離擾動(dòng)的影響，氣流流動(dòng)紊亂，從而湍流動(dòng)能比較大，致使湍流耗散率也相對(duì)較大，湍流動(dòng)能最后都轉(zhuǎn)化為熱能耗散掉，這也使汽車克服空氣阻力多做功。因此要改善汽車的氣動(dòng)性，就應(yīng)該控制尾渦系數(shù)的強(qiáng)度，減弱它的湍流程度，降低氣流的湍流能量消耗。

圖17 汽車尾部對(duì)稱面湍流動(dòng)能分布圖

4.3 結(jié)果計(jì)算

由圖9、圖10，得到原始模型的結(jié)果數(shù)據(jù)為：

利用CATIA軟件得到在y方向的投影面積，如圖18。即迎風(fēng)面積為：

圖18 車身模型在y方向的投影面

阻力系數(shù)：

根據(jù)企業(yè)報(bào)道，該款實(shí)車阻力系數(shù)為0.321。

誤差為：（0.321-0.317）/0.321=1.2%。滿足工程分析需要。

5 優(yōu)化方案

對(duì)汽車空氣動(dòng)力學(xué)性能影響因素有：前部上翹角（A），前風(fēng)窗角（B），后風(fēng)窗角（C），尾部上翹角（D）。其中前部上翹角主要影響汽車車身底部的氣流分布；前風(fēng)窗角主要影響在車身上部的氣流分離點(diǎn)；后風(fēng)窗角主要影響汽車車身的尾部氣流分布，尾部上翹角主要影響氣動(dòng)阻力大小。對(duì)A、B、C、D四因子分別提取大、中、小三種水平，生成車身三維模型。構(gòu)造正交實(shí)驗(yàn)表，如表4。

表4 正交實(shí)驗(yàn)表

借鑒了相關(guān)研究經(jīng)驗(yàn)[9]，在此范圍分為大、中、小三段，分別取代表值。數(shù)據(jù)的參數(shù)分類如下：

A：4-8，9-11，12-16。取7、10、14計(jì)算。

B：25-28，29-32，33-36。取27、30、34計(jì)算。

C：21-24，25-27，28-31。取22、26、30計(jì)算。

D：5-9，10-12，13-17。取7、11、15計(jì)算。

由此得到的參數(shù)數(shù)據(jù)如表5所示。

表5 模型計(jì)算數(shù)據(jù)表

根據(jù)表5的影響車身氣動(dòng)性的相關(guān)參數(shù)建立車身模型。圖19為選定第1組數(shù)據(jù)建立的參數(shù)化模型。

圖19 第1組參數(shù)化模型

將模型導(dǎo)入數(shù)值仿真軟件，得到各個(gè)參數(shù)模型所受到的氣動(dòng)阻力和氣動(dòng)升力整理如表6所示。由公式計(jì)算得到車身模型阻力系數(shù)和升力系數(shù)，整理后如表7所示。

為了盡可能的降低燃油消耗和增加附著力，所以要選擇阻力系數(shù)和升力系數(shù)都相對(duì)較小的參數(shù)組合。對(duì)比可知第 4組和第7組系數(shù)值較小。

表6 參數(shù)化車身模型計(jì)算結(jié)果

表7 參數(shù)化車身所對(duì)應(yīng)的阻力和升力系數(shù)值

可知，四個(gè)關(guān)鍵參數(shù)的取值范圍為：A是 10～14°，B是27°，C是26～30°，D是11～15°。

選定A為12°，B為27°，C為28°，D為13°。對(duì)該優(yōu)化后的模型進(jìn)行數(shù)值計(jì)算。選定上述參數(shù)后的優(yōu)化車身曲面如圖20所示。

圖20 優(yōu)化后車身曲面模型

數(shù)值計(jì)算結(jié)果顯示，優(yōu)化后模型車身外流場氣動(dòng)性能得到了較大改善。

優(yōu)化后結(jié)果數(shù)據(jù)為：

升力系數(shù)：

6 結(jié)論

（1）在車身上部分的氣流在引擎蓋端緣開始加速。在前風(fēng)窗與車頂相接處達(dá)到最大值。當(dāng)氣流到達(dá)車頂與后風(fēng)窗相接處的時(shí)候，氣流速度減小，并發(fā)生分離現(xiàn)象，在車身尾部形成了一個(gè)渦流區(qū)。

（2）在車身下部分的氣流則從車身底部和路面之間強(qiáng)制通過，氣流并未在車身尾部直接流走，而是沿汽車尾部與車身上表面的氣流匯合，并形成渦流區(qū)域，因此產(chǎn)生很大的阻力。

（3）由壓力云圖顯示，前風(fēng)窗角影響重大?？刹捎么髢A角設(shè)計(jì)。

（4）要改善汽車的氣動(dòng)性，就應(yīng)該控制尾渦系數(shù)的強(qiáng)度，減弱它的湍流程度，降低氣流的湍流能量消耗。

（5）通過正交實(shí)驗(yàn)，得到優(yōu)化模型。即：A為12°，B為27°，C為28°，D為13°。

參考文獻(xiàn)

[1] 張國忠,賴征海.汽車空氣動(dòng)力學(xué)與車身造型研究最新進(jìn)展[J].沈陽大學(xué)學(xué)報(bào),2005（6）:40-43.

[2] 陳振明,尹華鑫.改善汽車空氣動(dòng)力學(xué)性能的措施[J].公路與汽運(yùn).2007（5）:4-6.

[3] 徐永康.汽車氣動(dòng)升力的實(shí)驗(yàn)與仿真研究[D].長沙:湖南大學(xué),2014.

[4] 徐勛,許濤.汽車外流場的數(shù)值模擬[J].內(nèi)燃機(jī),2012（3）:26-28.

[5] 李學(xué)志,李若松.CATIA實(shí)用教程[M].北京:清華大學(xué)出版社,2011.

[6] Guilmineau E. Numerical simulations of flow around a realistic gene-ric car model[J].SAE International Journal of Passenger Cars Me-chanical Systems,2014（2）:646-653.

[7] 張偉.EQ1118GA 運(yùn)輸車外流場的數(shù)值模擬與分析[D].杭州:浙江大學(xué),2006.

[8] 安震,李旭,李居莉.某轎車外流場數(shù)值模擬與分析[J].農(nóng)業(yè)裝備與車輛工程,2012（12）:48-52.

[9] 龍鋼.某車型車身造型設(shè)計(jì)及空氣動(dòng)力學(xué)仿真研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2010.