基于改進(jìn)的Pareto遺傳算法的車身氣動多目標(biāo)優(yōu)化*

韋 甘，楊志剛，李啟良

(同濟(jì)大學(xué)上海地面交通工具風(fēng)洞中心，上海 201804)

前言

降低風(fēng)阻系數(shù)是汽車車身氣動優(yōu)化最重要的任務(wù)之一，通過風(fēng)洞試驗等方法可以設(shè)計出風(fēng)阻系數(shù)很低的車身。但傳統(tǒng)低阻車身都沒有多大的實用價值，因為這些低阻車身往往不符合傳統(tǒng)整車布置、機(jī)件布置、乘員舒適性等空間要求。車內(nèi)空間大小也是車身重要的性能之一，它和車身氣動性能往往發(fā)生沖突。例如，快背車型比直背車型的氣動性能優(yōu)越，但室內(nèi)空間就明顯減小。如何找到氣動性能和空間性能的平衡點(diǎn)，是設(shè)計出實用性較高的低阻車身的關(guān)鍵。

本文中把氣動性能和空間性能同時作為優(yōu)化目標(biāo)，把車身氣動優(yōu)化問題作為多目標(biāo)優(yōu)化問題?；赑areto最優(yōu)解的定義，參考Pareto遺傳算法在飛機(jī)翼型多目標(biāo)優(yōu)化問題上的應(yīng)用方法[1-3]，求解不同車型的Pareto波陣面。

Pareto多目標(biāo)遺傳算法常用方法有等級分類法、小生境技術(shù)、自適應(yīng)變異率等[3]。本文中根據(jù)車身多目標(biāo)優(yōu)化問題的特點(diǎn)，對常用的Pareto遺傳算法進(jìn)行改進(jìn)，以提高優(yōu)化效率和效果。

1 多目標(biāo)優(yōu)化方法

1.1 車身多目標(biāo)優(yōu)化問題

車內(nèi)空間性能是車身的重要性能之一。把車身氣動性能和空間性能作為優(yōu)化目標(biāo)進(jìn)行多目標(biāo)優(yōu)化。車身的氣動性能指標(biāo)采用風(fēng)阻系數(shù)Cd值。用二維車身(即車身側(cè)視圖)的剩余面積百分?jǐn)?shù)Ra和三維車身的剩余體積百分?jǐn)?shù)Rv作為二維車身和三維車身的空間性能指標(biāo)。Ra和Rv的定義如下：

Ra=(1-A/((H-D)×L))×100

Rv=(1-V/((H-D)×L×W))×100

式中：A為側(cè)視圖面積；V為三維車身體積；L、W、H和D分別為整車的長、寬、高和離地間隙。該多目標(biāo)優(yōu)化問題是雙最小值問題，得到的Pareto波陣面是一條向左下方凸出的弧線。

1.2 氣動性能的評估

圖1和圖2分別示出兩種典型造型車型的側(cè)視圖、俯視圖和橫截面圖。其中，普通車型是指至今最常見的、包括兩廂和三廂式的前部造型；流線型車型是指20世紀(jì)80年代前后才出現(xiàn)的、采用更強(qiáng)流線型的多功能車(MPV)之類汽車的造型，國內(nèi)俗稱“子彈頭”車，國外稱為“蛋形造型”。圖中用參數(shù)控制曲線關(guān)鍵點(diǎn)(圖中圓點(diǎn))的位置和部分曲線的端點(diǎn)斜率，用直線和保凸曲線構(gòu)成車身的側(cè)視圖、俯視圖和橫截面圖。其中曲線關(guān)鍵點(diǎn)主要由車身最高點(diǎn)、最低點(diǎn)、最前后端點(diǎn)、發(fā)動機(jī)罩與風(fēng)窗的交點(diǎn)等特征點(diǎn)及其之間的過渡點(diǎn)組成。結(jié)合半車身[4]和翼型彎度線[5]的概念，實現(xiàn)三維形體的構(gòu)建。圖3為三維車身的縮比樹脂模型。

把流線型車型側(cè)視圖的12個參數(shù)和普通車型側(cè)視圖的18個參數(shù)作為各自的優(yōu)化對象。俯視圖和橫截面圖的參數(shù)組合根據(jù)乘員舒適性等空間要求進(jìn)行設(shè)計并固定為如圖2所示的形狀。作為優(yōu)化約束條件之一的整車的長、寬、高和離地間隙設(shè)計在A級車尺寸范圍內(nèi)，其中長高比為3.0。

用Fluent商業(yè)軟件對車身外流場進(jìn)行模擬計算以求解其Cd值。來流速度設(shè)為30m/s，湍流模型選用可實現(xiàn)k-ε模型。在噴口面積為0.123m2的模型風(fēng)洞中對1∶12縮比模型進(jìn)行測力試驗(圖3)，支撐方式包括中央單支撐柱和四支撐柱兩種。車身Cd值的CFD計算結(jié)果和試驗結(jié)果誤差小于7%。

1.3 優(yōu)化策略

使用基于Pareto最優(yōu)解的Pareto多目標(biāo)遺傳算法作為優(yōu)化方法[1-3]。直接對三維車身進(jìn)行氣動優(yōu)化耗時較長。本文中把二維優(yōu)化和三維優(yōu)化相結(jié)合，先利用耗時較短的二維算例做全局搜索找到目標(biāo)Pareto波陣面的附近位置，然后將優(yōu)化得到的二維車身Pareto波陣面作為三維車身多目標(biāo)優(yōu)化階段的初代，繼續(xù)做進(jìn)一步優(yōu)化，局部逼近和加密Pareto波陣面。二維優(yōu)化階段的初代是用優(yōu)化拉丁方的方法設(shè)計的較大規(guī)模均勻分布群體。二維和三維優(yōu)化階段的群體大小都為60，進(jìn)化代數(shù)都為15。

2 改進(jìn)的Pareto遺傳算法

單目標(biāo)優(yōu)化結(jié)果是一個點(diǎn)，多目標(biāo)優(yōu)化問題得到的是Pareto波陣面，后者可以看做是前者的組合。固定Rv或Ra為定值，優(yōu)化得到的Cd最小點(diǎn)即為 Pareto波陣面上的其中一個點(diǎn)，許多這樣的點(diǎn)組成了整個Pareto波陣面。當(dāng)求解波陣面上的某一個點(diǎn)時，按照單目標(biāo)優(yōu)化的方法應(yīng)該是就近選擇附近區(qū)域內(nèi)的點(diǎn)逐步逼近該點(diǎn)最優(yōu)解。常用的Pareto遺傳算法為了保證Pareto波陣面的均勻性，往往強(qiáng)調(diào)保持物種多樣性[6]，例如遠(yuǎn)親交配制度選用距離較大的個體作為子代[3]。這些方法忽略了鄰近個體之間存在的近似線性關(guān)系。尤其在作為局部逼近階段的三維優(yōu)化階段，鄰近個體能夠提供重要信息指導(dǎo)進(jìn)化方向。

改進(jìn)的遺傳算法將總庫中的現(xiàn)有樣本按角度均勻劃分到多個小家族中，然后以小家族為單位進(jìn)行進(jìn)化。這樣既可以保證樣本的多樣性和均勻度，又可以最大限度地利用鄰近樣本點(diǎn)的信息，提高優(yōu)化效率。所有已經(jīng)計算過的樣本點(diǎn)都納入二維或三維車身總庫中，為下一代的進(jìn)化盡量提供更多的信息。下面以耗時較少的二維車身的單代進(jìn)化為例，講解和驗證改進(jìn)的Pareto遺傳算法。

2.1 單目標(biāo)預(yù)測

提前預(yù)測無須通過CFD計算求解的Rv或Ra值，可以提供更精確的進(jìn)化方向，提高優(yōu)化效率。Ra值可以通過Matlab程序快速獲得，而Rv值的單算例計算時間長達(dá)半個小時。流線型車型和普通車型的Rv值和Ra值用線性方程擬合后得到的決定系數(shù)分別為0.978和0.989，線性關(guān)系明顯，這說明用Ra的變化關(guān)系預(yù)測Rv的變化關(guān)系是可行的。為了在保證精度的同時盡量減少預(yù)測時間，預(yù)測Ra時所取橫截面間隔為三維建模時的4倍。

2.2 鄰點(diǎn)交叉法

(1)劃分小家族 使用如圖4(a)所示的角度劃分方法1進(jìn)行小家族劃分。圖中圓點(diǎn)為總庫中的現(xiàn)有樣本點(diǎn)，圓圈內(nèi)點(diǎn)X1、Xj、Xm為上一代Pareto最優(yōu)解。用個體的目標(biāo)值除以最大值進(jìn)行歸一化，以相鄰兩點(diǎn)與中心點(diǎn)o′的連線為各個小家族的分界線，根據(jù)夾角φj+1占90°的比例大小分配各個小家族將要產(chǎn)生的子代個數(shù)。

(2)交叉 分配好小家族的子代個數(shù)后，以小家族的兩個端點(diǎn)作為親代進(jìn)行部分交叉然后進(jìn)行變異。將初步得到的子代用單目標(biāo)預(yù)測功能進(jìn)行篩選，使最終確定的子代盡量落在綜合性能較優(yōu)的區(qū)域內(nèi)。例如，只保留空間性能優(yōu)于兩個親代的個體作為子代，如果一次變異得到的符合要求的子代個數(shù)不夠，繼續(xù)對交叉得到的個體進(jìn)行變異操作，直到滿足該小家族的子代數(shù)目為止。

2.3 分象限外推法

(1)劃分小家族 在角度劃分法1的基礎(chǔ)上稍作改動，以上一代Pareto最優(yōu)解的夾角中分線作為每個小家族的分界線，得到如圖4(b)所示的角度劃分法2。按歸一化后的樣本點(diǎn)的角度所占比例來分配子代個數(shù)的角度劃分法可以最大限度地保證波陣面上各個區(qū)域的點(diǎn)均勻分布。

(2)劃分象限 以任何個體為原點(diǎn)都可以把其余個體按相對角度劃分到4個象限中。例如，對于點(diǎn)Xj，所有Ra>Ra(Xj)且Cd>Ra(Xj)的點(diǎn)都被劃分到第一象限中。

(3)外推變異 選擇每個小家族中的Pareto最優(yōu)解為原點(diǎn)，選擇特定象限中與最優(yōu)解距離較近的點(diǎn)為推點(diǎn)，按下式進(jìn)行變異：

式中：i表示基因位；Xoi為子代；Xai為原點(diǎn)；Xbi為推點(diǎn)，d為變異范圍，隨機(jī)數(shù)r∈[0,1]。

外推變異以原點(diǎn)為唯一親代，搜索方向從推點(diǎn)指向原點(diǎn)，目的是讓子代盡量落在推點(diǎn)所在象限的相隔象限。單目標(biāo)預(yù)測功能可以使子代落在目標(biāo)象限的概率增大。按照每個小家族中推點(diǎn)的數(shù)目分為單點(diǎn)外推法和多點(diǎn)外推法。

對于Pareto波陣面上的中間點(diǎn)，使用第1象限點(diǎn)為推點(diǎn)，目標(biāo)是讓子代落在第3象限，將Pareto波陣面向左下角外推。為了使Pareto波陣面向左上方和右下方延伸和加密，對于上端點(diǎn)(靠近0°水平線的端點(diǎn))，使用第4象限點(diǎn)為推點(diǎn)；對于下端點(diǎn)，使用第2象限點(diǎn)為推點(diǎn)。

2.4 進(jìn)化性能評估指標(biāo)

(1)分布性能評估指標(biāo) 引用小生境技術(shù)中的小生境數(shù)的概念來衡量子代個體沿周向的擁擠程度。小生境技術(shù)中的共享函數(shù)定義為

式中：σshare為共享半徑；dij為i和j兩個個體之間的歐拉距離。把衡量個體間擁擠程度的小生境數(shù)mi定義為共享函數(shù)之和，所有個體的平均擁擠程度定義為群體中所有個體的小生境數(shù)之和除以群體數(shù)，即：

式中：M為群體的平均小生境數(shù)。為了消除不同目標(biāo)函數(shù)值尺度和單位的影響，先根據(jù)圖4所示的方法計算個體的角度值，再根據(jù)角度值計算平均小生境數(shù)。

(2)進(jìn)化效果評估指標(biāo) 把圖4(a)所示的所有個體與中心點(diǎn)o′的連線長度Li的平均值L作為評估子代進(jìn)化效果的指標(biāo)，稱L為平均半徑，其定義為

2.5 二維車身的多目標(biāo)優(yōu)化測評

下面用PT、LJ、WT分別表示普通遺傳算法、鄰點(diǎn)交叉法和分象限外推法的數(shù)值實驗算例。PT3、4，LJ1-4和WT1-3使用了單目標(biāo)預(yù)測算子。PT3、4，LJ2、4、6和WT3、5、7的變異步長為參數(shù)變化范圍的20%～40%，其余算例的變異步長是它們的1/2。PT1-4只使用了變異算子，PT5只使用了交叉算子。LJ3、4和LJ1、2的單目標(biāo)預(yù)測算子的目標(biāo)區(qū)域不同。WT1、4、5使用的是單點(diǎn)外推法，其余WT算例使用的是多點(diǎn)外推法。不同進(jìn)化參數(shù)和進(jìn)化策略可以檢驗優(yōu)化方法的魯棒性。圖5為不同優(yōu)化方法的單代進(jìn)化實例的子代平均小生境數(shù)M和平均半徑L。

LJ系列和WT系列的子代的平均小生境數(shù)M值普遍低于PT系列的M值，L值普遍高于PT系列的L值，根據(jù)平均小生境數(shù)和平均半徑的定義，LJ系列和WT系列的子代擁擠程度明顯比PT系列的低，子代平均進(jìn)化性能明顯比PT系列高，這說明改進(jìn)后的鄰點(diǎn)交叉法和分象限外推法的優(yōu)化效率和優(yōu)化效果比普通的優(yōu)化方法有明顯提高。

PT系列的M值和L值的變化范圍分別約為1.6和0.06，M值和L值的變化范圍都較大，說明普通進(jìn)化方法對進(jìn)化參數(shù)和進(jìn)化方案較敏感。LJ系列的M值和L值的變化范圍分別約為0.6和0.02；WT系列的M值和L值的變化范圍分別約為0.3和0.04。和PT系列相比，LJ系列和WT系列的M值和L值的變化范圍明顯較小，說明改進(jìn)后得到的鄰點(diǎn)交叉法和分象限外推法對進(jìn)化參數(shù)的敏感性降低，魯棒性能有所增強(qiáng)。

圖6列出了各系列中進(jìn)化效果最優(yōu)的部分算例，其中標(biāo)題帶(Ra)的算例使用了單目標(biāo)預(yù)測功能。圖中圓圈為上一代Pareto波陣面，點(diǎn)為子代。

LJ系列和WT系列與PT系列相比，前兩者的子代有更多的樣本點(diǎn)分布在上一代波陣面的左下方，且在Rv>30%和Cd>0.3的區(qū)域有最多的個體分布，子代分布更狹長。其中WT系列的子代比LJ系列的子代分布更均勻，性能更優(yōu)秀。使用了單目標(biāo)預(yù)測算子的PT3、LJ3、WT3與其余算例相比，前者的子代大部分都落在初代波陣面外部而后者的子代大部分落在初代波陣面內(nèi)部，可見單目標(biāo)預(yù)測算子與鄰點(diǎn)交叉法、分象限外推法的配合可以明顯提高進(jìn)化效率。

3 優(yōu)化結(jié)果分析

通過二維和三維兩個優(yōu)化階段最終優(yōu)化得到的Pareto波陣面如圖7所示。兩種車型的Pareto波陣面基本都均勻分布，且覆蓋范圍較廣，Cd變化范圍約在0.07～0.14之間，Rv變化范圍約在20%～40%之間，其中普通車型落在兩端位置的最優(yōu)解更多，Pareto波陣面分布范圍更廣。

三維無輪車身Cd值在0.07～0.09范圍內(nèi)屬于低阻區(qū)域[7]。在Cd≤0.11的中低阻區(qū)，流線型車型的Pareto最優(yōu)解比普通車型的更向左下角突出，綜合性能更優(yōu)。在中低阻區(qū)，偏低的尾部是車身的最大特點(diǎn)，尾部空間較小使得頭部空間對Rv的影響較顯著。而流線型車型的頭部空間調(diào)節(jié)幅度比普通車型更大，相同的Cd值下Rv可以達(dá)到更小值，所以Pareto波陣面有下移的趨勢。在Cd>0.11的高阻區(qū)，情況正好相反，普通車型的Pareto最優(yōu)解比流線型車型的更向左下角突出，綜合性能更優(yōu)秀。在高阻區(qū)，偏高的尾部是車身的最大特點(diǎn)。較大的尾部空間減弱了頭部空間對Rv的影響。普通車型因為在側(cè)視圖的頭部使用了更多的控制點(diǎn)，外形變化更自由，通過調(diào)整頭部附近氣流來降低車身氣動阻力的能力更優(yōu)，同樣的Rv值下Cd可以達(dá)到更小值，所以Pareto波陣面有左移的趨勢。

圖8為圖7中字母標(biāo)出點(diǎn)所對應(yīng)的車身表面油流圖。左上端點(diǎn)的Cd值最小，Rv最大，形體特點(diǎn)是尾部低而且長；右下端點(diǎn)的Cd值最大，Rv最小，形體特點(diǎn)是車尾高，接近矩形。普通車型的Pareto最優(yōu)解的頭部凹陷均不明顯，尤其是b、c點(diǎn)，且所選4個點(diǎn)與對應(yīng)的流線型車型的4個點(diǎn)的背部、尾部、底部后翹的造型非常相似，這是因為尾部造型是決定Cd和Rv的最重要因素。

流線型車型車身A的尾部兩側(cè)拖曳渦形成和分離的痕跡較明顯，而車身C、D的尾部兩側(cè)流線平直。普通車型的示例車身a-d有類似的現(xiàn)象。這說明當(dāng)車身尾部較低長時，外流場的三維效應(yīng)較明顯，尾部兩側(cè)拖曳渦較明顯，誘導(dǎo)阻力占風(fēng)阻的分量上升。當(dāng)車身尾部較高時，外流場的三維效應(yīng)減弱，兩側(cè)拖曳渦不明顯，甚至沒有生成拖曳渦，氣流在截尾處突然分離，壓差阻力占風(fēng)阻的分量上升。

4 結(jié)論

(1)鄰點(diǎn)交叉法和分象限外推法和單目標(biāo)預(yù)測法配合使用，可以有效降低平均子代的平均生境數(shù)，增加平均半徑，提高子代的分布性能和優(yōu)化效率，提高優(yōu)化方法的魯棒性。這種改進(jìn)的Pareto遺傳算法可以解決車身氣動多目標(biāo)優(yōu)化中計算量過大的問題。

(2)把車身氣動性能和空間性能作為優(yōu)化目標(biāo)，用改進(jìn)的Pareto遺傳算法求解得到基本均勻分布的Pareto波陣面，其中Cd變化范圍約在0.07～0.14之間。在中低阻區(qū)，流線型車型的Pareto最優(yōu)解更向左下角突出，綜合性能更優(yōu)；在高阻區(qū)，普通車型的Pareto最優(yōu)解更向左下角突出，綜合性能更優(yōu)。

[1] Franco Mastroddi, Stefania Gemma. Analysis of Pareto Frontiers for Multidisciplinary Design Optimization of Aircraft[J]. Aerospace Science and Technology, 2013,28(1):40-55.

[2] Lee Sang-Hwan, Lee Juhee. Aerodynamic Analysis and Multi-objective Optimization of Wings in Ground Effect[J]. Ocean Engineering，2013， 68(8)：1-13.

[3] 王曉鵬.遺傳算法及其在氣動優(yōu)化設(shè)計中的應(yīng)用研究[D].西安:西北工業(yè)大學(xué)，2001.

[4] 傅立敏. 汽車空氣動力學(xué)[M]. 北京:機(jī)械工業(yè)出版社, 2006,71-72.

[5] 何憶斌,谷正氣,李偉平,等. 汽車?yán)硐霘鈩有误w數(shù)字化模型構(gòu)建及氣動性能試驗[J]. 航空動力學(xué)報，2011, 25(6): 1031-1035.

[6] Jaroslav Hájeka, András Sz?ll?sa,Jakubístekb. A New Mechanism for Maintaining Diversity of Pareto Archive in Multi-objective Optimization[J]. Advances in Engineering Software, 2010,41(7-8):1031-1057.

[7] Buchheim R, Deutenbah K R, Luckoff H J. Necessity and Premises for Reducing the Aerodynamic Drag of Future Passenger Cars[C]. SAE Technical Paper, 1982:758-771.