劉傳波，劉 鈞

(武漢理工大學(xué) 機(jī)電工程學(xué)院，湖北 武漢 430070)

風(fēng)阻靈敏度是指在一定范圍內(nèi)，汽車造型的改變，如后視鏡、門把手、底盤離地高度、接近角等優(yōu)化;或者某些零部件的調(diào)整，如進(jìn)氣格柵、前風(fēng)窗傾角、擾流板等參數(shù)變化對風(fēng)阻系數(shù)的影響。前者相對于后者，對風(fēng)阻系數(shù)的影響較小。常用風(fēng)阻系數(shù)優(yōu)化主要依靠最佳實(shí)踐及工程經(jīng)驗(yàn)對主體外形特征參數(shù)進(jìn)行空氣動力學(xué)的分析和評價(jià)。

1991年，Takuya等[1]建立了一個(gè)基于CFD的關(guān)于汽車空氣動力學(xué)的計(jì)算分析系統(tǒng)，該系統(tǒng)包括數(shù)學(xué)計(jì)算和計(jì)算網(wǎng)格生成，運(yùn)用多塊法來解決復(fù)雜的問題。最終取得的試驗(yàn)?zāi)M數(shù)據(jù)和風(fēng)洞試驗(yàn)得到的數(shù)據(jù)相比很理想。2013年，法雷奧公司散熱系統(tǒng)的研究人員Samer等[2]通過風(fēng)洞試驗(yàn)，提出了一種通過減少發(fā)動機(jī)艙開口面積來達(dá)到降低風(fēng)阻的方案。黃永輝等[3]通過對某型車的研究發(fā)現(xiàn)，在模型車上增加擾流板和改變后風(fēng)窗高度等參數(shù)，車身側(cè)面的正壓區(qū)會明顯減小，而且增加后擾流板，對降低整車風(fēng)阻系數(shù)而言是一個(gè)非常有效的方案。孫琪等[4]先采取調(diào)查問卷的方式，確定對風(fēng)阻系數(shù)有影響的特征參數(shù)，然后經(jīng)過試驗(yàn)仿真，進(jìn)一步篩選特征參數(shù)，最后直接對其進(jìn)行正交試驗(yàn)分析，得到最終優(yōu)化模型。谷正氣等[5]通過在汽車尾部加裝噴射裝置，破壞其尾渦，從而來減小汽車尾部的負(fù)壓區(qū)域，降低整車風(fēng)阻系數(shù)。

以往在降低汽車風(fēng)阻的研究過程中，主要是基于風(fēng)洞試驗(yàn)改進(jìn)和設(shè)計(jì)整車造型，來找到對風(fēng)阻作用比較大的特征參數(shù)，然后再從中不斷地優(yōu)化，這樣得到的結(jié)果比較精確可信，缺點(diǎn)是資源要求比較大，而且大多數(shù)是單因素分析方式，筆者的研究主要采用計(jì)算流體動力學(xué)和正交試驗(yàn)相結(jié)合的方式，先通過整車流場分析，確定特征參數(shù)，進(jìn)行單因素分析并進(jìn)行局部尋優(yōu)并驗(yàn)證特征參數(shù)的重要性，最后進(jìn)行正交試驗(yàn)組合分析尋找全局最優(yōu)，以期得到耦合分析的最佳風(fēng)阻方案。

1 特征參數(shù)獲取

1.1 原車數(shù)學(xué)模型

研究對象為英國汽車工程協(xié)會MIRA(Motor Industry Research Association)快背式模型，該模型在乘用車上大量應(yīng)用，具有較大的現(xiàn)實(shí)意義。研究采用ANSYS Fluent仿真軟件對整車外流場進(jìn)行仿真計(jì)算，考慮到計(jì)算機(jī)硬件配置和計(jì)算機(jī)資源的合理利用，采用1∶5建立MIRA快背式模型，在相同的初始條件和邊界條件下能夠接近1∶1模型仿真計(jì)算結(jié)果[6]。

在ANSA中劃分網(wǎng)格，并建立計(jì)算域模型[7]：車身長度為L，高度為H，寬度為W，計(jì)算域幾何參數(shù)設(shè)置為：車前部3L，車后部5L，車頂部3H，車兩側(cè)各2W。同時(shí)為了仿真計(jì)算更加精確，在計(jì)算域中劃分一部分區(qū)域進(jìn)行加密：車前部1L，車頂部1H，車兩側(cè)各1W，車后部2L，如圖1所示。最后生成351萬體網(wǎng)格。

圖1 計(jì)算域

汽車造型任何參數(shù)的改變都會對整車風(fēng)阻系數(shù)造成一定的影響，門把手、后視鏡、雨刷和保險(xiǎn)桿等部件對整車風(fēng)阻系數(shù)影響較小，在本次建模分析中，將其忽略。選擇k-e湍流模型，定義空氣為不可壓縮氣體[8]。

計(jì)算域入口為空氣流入口，即速度入口，根據(jù)MIRA數(shù)據(jù)庫中的仿真分析條件，參數(shù)保持一致，速度設(shè)置為30 m/s。計(jì)算域出口為空氣流出口，即壓力出口，設(shè)置大氣壓為0 Pa,因?yàn)樵诜抡婺M分析中汽車和壁面是固定不動的，因此設(shè)置為壁面，同時(shí)為了更加精確地模擬車身底盤下的空氣流場，因?yàn)樵诒敬文M仿真中，車身是固定不動的，為了模擬相對運(yùn)動，設(shè)置內(nèi)部加密區(qū)域的下地面為移動壁面，同時(shí)給定速度為30 m/s,大計(jì)算域的下底面為固定面，設(shè)置為壁面，劃分的加密區(qū)域，在仿真計(jì)算時(shí)，空氣是要流經(jīng)該區(qū)域的，故設(shè)置為室內(nèi),邊界條件如表1所示。

1.2 仿真結(jié)果分析

根據(jù)MIRA所提供的快背式風(fēng)阻系數(shù)，對所建立模型進(jìn)行對比分析，并提取出對風(fēng)阻影響較大的特征參數(shù)。快背式模型的風(fēng)阻系數(shù)為0.27，本次仿真計(jì)算所得結(jié)果Cd=0.277 131，與給出的參考數(shù)據(jù)極為接近。模型的壓力云圖、對稱面壓力云圖和對稱面速度流線圖如圖2～圖4所示。

圖2 壓力云圖

圖3 對稱面壓力云圖

圖4 對稱面速度流線圖

從圖2可以看出，汽車車頭、發(fā)動機(jī)罩和前風(fēng)窗之間的壓力比較大，同時(shí)氣流在前風(fēng)窗和車頂之間的部分還存在一定的負(fù)壓區(qū)域；根據(jù)圖4也可以發(fā)現(xiàn)在模型車尾部分也存在著一定范圍的負(fù)壓區(qū)域。

從整車外流場仿真分析結(jié)果中，可選擇3個(gè)對風(fēng)阻影響較大的特征參數(shù)：發(fā)動機(jī)罩與水平面夾角(10°)、前風(fēng)窗與水平面夾角(45°)和離去角(10°)。

2 特征參數(shù)單因素分析

筆者的研究主要是獲得MIRA快背式車型每個(gè)特征參數(shù)對整車風(fēng)阻系數(shù)的影響，以及特征參數(shù)風(fēng)阻靈敏度的最佳組合。先對所提取出的3個(gè)特征參數(shù)進(jìn)行仿真研究，根據(jù)相關(guān)研究結(jié)果和實(shí)際情況，最終確定離去角仿真計(jì)算范圍為5°～15°，發(fā)動機(jī)罩與水平面夾角仿真計(jì)算范圍為5°～15°，前風(fēng)窗與水平面夾角仿真計(jì)算范圍為30°～50°。分析三個(gè)特征參數(shù)在單一情況下的風(fēng)阻靈敏度，驗(yàn)證其幾何參數(shù)改變對整車風(fēng)阻的重要性；并為最終耦合分析提供數(shù)據(jù)基礎(chǔ)。

2.1 離去角變化對風(fēng)阻系數(shù)影響

保持前風(fēng)窗與水平面夾角和發(fā)動機(jī)罩與水平面夾角幾何參數(shù)不變，將離去角上下變動5°，取5°和15°進(jìn)行建模仿真分析，仿真結(jié)果如表2所示。

表2 離去角變化仿真數(shù)據(jù)

在離去角為5°時(shí)，風(fēng)阻系數(shù)降低了1.02%，當(dāng)離去角增大到15°時(shí)，風(fēng)阻系數(shù)明顯變大，增加了5.35%，離去角的增大使得汽車后部負(fù)壓區(qū)域變大，同時(shí)壓差阻力也顯著增大。

2.2 發(fā)動機(jī)罩與水平面夾角變化對風(fēng)阻系數(shù)影響

保持離去角和前風(fēng)窗水平面夾角不變，將發(fā)動機(jī)罩與水平面夾角依次取5°、15°和20°進(jìn)行仿真建模分析，得到其在不同幾何參數(shù)下的風(fēng)阻系數(shù)，如表3所示。

表3 發(fā)動機(jī)罩與水平面夾角變化的仿真數(shù)據(jù)

通過仿真結(jié)果可以發(fā)現(xiàn)，隨著發(fā)動機(jī)罩與水平面夾角不斷變大，所研究車型前部的面積也不斷減小，風(fēng)阻系數(shù)降低了4.11%。

2.3 前風(fēng)窗與水平面夾角變化對風(fēng)阻系數(shù)影響

保持離去角和發(fā)動機(jī)罩與水平面夾角不變，將前風(fēng)窗與水平面夾角依次取30°、35°、40°、50°進(jìn)行仿真建模分析，仿真結(jié)果如表4所示。

表4 前風(fēng)窗水平面夾角變化仿真數(shù)據(jù)

通過仿真計(jì)算結(jié)果可以發(fā)現(xiàn)，在前風(fēng)窗與水平面夾角為30°的情況下，整車的風(fēng)阻系數(shù)最小，降低了1.5%，在50°的情況下，整車風(fēng)阻系數(shù)明顯變大，增加了17.04%，大大增大了整車的油耗。

2.4 綜合分析

根據(jù)上述仿真計(jì)算結(jié)果可知，通過依次改變所提取的特征參數(shù)，離去角、前風(fēng)窗與水平面夾角和發(fā)動機(jī)罩與水平面夾角均能在一定程度上降低整車的風(fēng)阻系數(shù)；尤其是降低發(fā)動機(jī)罩與水平面夾角，使得整車風(fēng)阻系數(shù)下降了4.11%。選取的3個(gè)特征參數(shù)的變動對整車風(fēng)阻系數(shù)數(shù)值改變的折線圖如圖5所示。

圖5 各特征參數(shù)風(fēng)阻系數(shù)變化圖

2.5 優(yōu)化模型分析

根據(jù)原車型仿真分析和提取出的特征參數(shù)仿真分析，選取離去角5°，前風(fēng)窗與水平面夾角30°，發(fā)動機(jī)罩與水平面夾角20°，重新建模，劃分網(wǎng)格并進(jìn)行仿真計(jì)算。邊界條件保持和初始仿真計(jì)算模型一致，導(dǎo)入到Fluent中進(jìn)行計(jì)算求解，得到仿真結(jié)果如圖6～圖8所示。

圖6 優(yōu)化模型壓力云圖

圖7 優(yōu)化模型對稱面壓力云圖

圖8 優(yōu)化模型對稱面速度流線圖

優(yōu)化模型的風(fēng)阻系數(shù)仿真值為0.259 921，比初始MIRA模型的風(fēng)阻系數(shù)減小了6.2%。對比初始模型可以發(fā)現(xiàn)：①車頂部和尾部的負(fù)壓區(qū)域明顯減小，尤其是前風(fēng)窗和發(fā)動機(jī)罩兩者之前的范圍。②車頭壓力區(qū)域也得到改善，優(yōu)化后的模型較好地減小了整車風(fēng)阻系數(shù)。③該3個(gè)特征參數(shù)經(jīng)過局部最優(yōu)后，通過壓力云圖發(fā)現(xiàn)其壓力值仍大于其他參數(shù)部分的壓力值，進(jìn)一步證明所選取的參數(shù)為特征參數(shù)。

3 耦合分析

上述研究是提取特征參數(shù)的單一因素變化對整車風(fēng)阻的影響，并沒有考慮到相互組合時(shí)對風(fēng)阻系數(shù)的作用，而正交試驗(yàn)?zāi)軌蛲ㄟ^選取一定量有代表性的試驗(yàn)，再加上合理科學(xué)的統(tǒng)計(jì)分析方法，高效地找到低風(fēng)阻情況下特征參數(shù)的組合和該組合情況下的風(fēng)阻靈敏度，而且可以研究分析每一個(gè)因素對整車風(fēng)阻系數(shù)影響的大小。

筆者研究的主要目的是降低整車的風(fēng)阻系數(shù)，得到不同特征參數(shù)在不同幾何參數(shù)下的風(fēng)阻靈敏度和其相應(yīng)的風(fēng)阻系數(shù)，因此將風(fēng)阻系數(shù)Cd的大小作為正交試驗(yàn)所優(yōu)化的目標(biāo)。在單因素分析中，選取了3個(gè)特征參數(shù)：①離去角在5°時(shí)風(fēng)阻系數(shù)降低最明顯，因此選取5°±2°為因素A(離去角角度)。②前風(fēng)窗與水平面夾角在30°時(shí)風(fēng)阻系數(shù)最低，因此選取30°±3°為因素B(前風(fēng)窗與水平面夾角)。③發(fā)動機(jī)罩與水平面夾角在20°時(shí)風(fēng)阻系數(shù)最低，因此選取20°±3°為因素C(發(fā)動機(jī)罩與水平面夾角)。以此為基礎(chǔ)，設(shè)計(jì)的正交試驗(yàn)水平和因素如表5所示。

根據(jù)正交試驗(yàn)設(shè)計(jì)水平和因素可知，正交試驗(yàn)為3因素3水平試驗(yàn)，根據(jù)經(jīng)驗(yàn)選取L9(34)4因素3水平正交試驗(yàn)表[9]，多出的一列用作誤差分析列，所設(shè)計(jì)的正交試驗(yàn)表頭如表6所示。

3.1 正交試驗(yàn)數(shù)據(jù)分析

表5 正交試驗(yàn)各因素取值表

表6 正交試驗(yàn)表頭設(shè)計(jì)

從表7可知，在A因素的3個(gè)水平中指標(biāo)K值最小的是A2，B因素的3個(gè)水平中指標(biāo)K值最小的是B1，C因素的3個(gè)水平中,指標(biāo)K最小的是C2。同時(shí)通過極差分析可以發(fā)現(xiàn)，在所選取的3個(gè)特征參數(shù)中B因素對整車風(fēng)阻系數(shù)影響最大,A因素和C因素對試驗(yàn)仿真結(jié)果影響較為接近。根據(jù)上述分析結(jié)果可知,降低風(fēng)阻系數(shù)的最佳特征參數(shù)組合為B1A2C2,即前風(fēng)窗與水平面夾角為30°,離去角為7°,發(fā)動機(jī)罩與水平面夾角為23°。

通過極差分析可知，B因素是對風(fēng)阻影響最大的因子，但尚無法確定B因素是否為顯著性因子，需要進(jìn)行顯著性檢驗(yàn)。在極差分析的基礎(chǔ)上進(jìn)一步進(jìn)行方差分析，離差平方和S表明了該因素水平變化所造成仿真試驗(yàn)結(jié)果的差異,方差分析結(jié)果如表8所示。

從表8可以看出，SA、SC和誤差列Se相差不大，因此相對于因素B而言,可以認(rèn)為因素A和因素C對試驗(yàn)指標(biāo)影響較小，可以將SA、SC、Se合并來估算誤差影響的大??；而且誤差項(xiàng)自由度越大，顯著性檢驗(yàn)也會更加精確。因素B的顯著性檢驗(yàn)公式為：

表7 正交試驗(yàn)結(jié)果的極差分析

表8 正交試驗(yàn)方差分析

(1)

式中：VB為因素B的均方值；V＇誤為調(diào)整后的誤差均方值。

V＇誤誤差均方值計(jì)算公式為：

(2)

式中：S為離差平方和；f為列自由度。計(jì)算顯著性檢查結(jié)果FB(2，6)=5.4，通過查臨界值表，F(xiàn)0.05(2，6)=5.14

3.2 最優(yōu)方案驗(yàn)證

選擇B1A2C2為最終優(yōu)化方案：前風(fēng)窗與水平面夾角為30°,離去角為7°,發(fā)動機(jī)罩與水平面夾角為23°。根據(jù)最佳特征參數(shù)組合方案進(jìn)行仿真計(jì)算[10]，最終模型仿真結(jié)果的對稱面壓力云圖如圖9所示。

圖9 最佳方案模型對稱面壓力云圖

從圖9可知，相對于初始的MIRA快背式車型，在選取的3個(gè)特征參數(shù)對應(yīng)的部分區(qū)域，壓力得到了明顯改善，最后得到風(fēng)阻系數(shù)Cd=0.252 031。相比初始MIRA快背式模型的風(fēng)阻系數(shù)降低了9.06%。

4 總結(jié)

以MIRA快背式車型為研究實(shí)例，先根據(jù)外流場分析，提取了3個(gè)特征參數(shù)；然后采用CFD仿真和正交試驗(yàn)分析相結(jié)合的方法，完成了該車3個(gè)特征參數(shù)的風(fēng)阻靈敏度分析和最后的優(yōu)化模型計(jì)算，得出以下結(jié)論：

(1)綜合分析了MIRA快背式車型離去角、前風(fēng)窗與水平面夾角、發(fā)動機(jī)罩與水平面夾角對整車風(fēng)阻系數(shù)的影響。研究結(jié)果表明，該3個(gè)特征參數(shù)角度的改變，均能在一定程度上改善了整車風(fēng)阻系數(shù)，其中發(fā)動機(jī)罩與水平面夾角和前風(fēng)窗與水平面夾角改變時(shí)，整車的風(fēng)阻系數(shù)減小較為明顯；離去角對整車風(fēng)阻系數(shù)的減小影響較低。

(2)通過單因素分析和正交試驗(yàn)分析，得到了最佳的3個(gè)特征參數(shù)的低風(fēng)阻組合方案，優(yōu)化后的MIAR快背式車型整車風(fēng)阻系數(shù)比初始模型降低了9.06%。

(3)通過正交試驗(yàn)結(jié)果分析中的極差分析和顯著性檢驗(yàn)可以發(fā)現(xiàn)，因素B(前風(fēng)窗與水平面夾角)對整車風(fēng)阻系數(shù)的降低效果最為顯著，在乘用車風(fēng)阻分析中可以重點(diǎn)研究前風(fēng)窗與水平面夾角的最佳風(fēng)阻幾何造型。

[1] Takuya Kataoka, Kunio Nakagawa, Masahiro Yoshida. Numerical Simulation of Road Vehicle Aerodynamics and Effect of Aerodynamic Devices[J].SAE Technical Paper,910579.

[2] Samer Saab, Alain Maiboom. Impact of the Underhood Opening Area on the Drag Coefficience and the Thermal Perfomance of a Vehicle[J].SAE International, 2013(1):869-879.

[3] 黃永輝.基于CFD的汽車外流場數(shù)值模擬及車身造型優(yōu)化分析[D].長沙：湖南大學(xué)，2011.

[4] 孫琪.汽車車身外流場空氣動力模擬及造型優(yōu)化[D].上海：上海工程技術(shù)大學(xué)，2014.

[5] 谷正氣，李學(xué)武，何億斌.汽車減阻新方法[J].汽車工程，2008,30(5)：441-443.

[6] 傅立敏.汽車設(shè)計(jì)與空氣動力學(xué)[M].北京：機(jī)械工業(yè)出版社，2006.

[7] Kisun Song, Kyung Hak Choo .A Spine-based Modeling Algorithm for Application to Aerodynamic Shape Optimization Based on CFD Analysis [J].SAE International, 2017(1):1510-1523.

[8] Lopes A M G, Carvalheira P. On the Application of Numerical Methods for the Calculation of the External Aerodynamics of a Streamlined Car Body[J]. SAE International,2003(1):1249-1256.

[9] 劉瑞江，張業(yè)旺.正交試驗(yàn)設(shè)計(jì)和分析方法研究[J].實(shí)驗(yàn)技術(shù)與管理，2010，27(9)：52-55.

[10] 莫易敏，雷志丹，楊宇.某MPV車型發(fā)型發(fā)動機(jī)艙散熱性能仿真與優(yōu)化研究[J].數(shù)字制造科學(xué)，2017，15(1-2)：1-6.