文潤旭，袁麗蕓，潘 鶴，王林波，韓小曼

（1.中國航發(fā)貴州紅林航空動力控制科技有限公司 產(chǎn)品設(shè)計所，貴州 貴陽 550002；2.廣西科技大學(xué) 機械與汽車工程學(xué)院，廣西 柳州 545006；3.空軍裝備部成都局駐貴陽地區(qū)第二軍事代表室，貴州 貴陽 550002）

0 引言

傳統(tǒng)有限元法具有計算精度不高、收斂性較差、“過硬”等的現(xiàn)象，因此，國內(nèi)外研究人員對該領(lǐng)域進行了深入研究，Liu[1-2]對傳統(tǒng)有限元法中的應(yīng)變矩陣進行了平滑梯度運算，從而避免了傳統(tǒng)有限元法的局限性，建立了一種計算精度較高的數(shù)值方法是光滑有限元法。Li 等[3]針對結(jié)構(gòu)-聲學(xué)耦合模型聲振特性模型，用光滑有限元法建立結(jié)構(gòu)-聲學(xué)耦合模型，并且通過數(shù)值算例，光滑有限元法比有限元法精度更高，具有對有限元“過硬”情況軟化的效果。柴應(yīng)彬[4]把光滑有限元法推廣應(yīng)用到船舶聲固耦合聲振分析中，與有限元法相比，光滑有限元表現(xiàn)出較高精確性。劉嘉明[5]基于FEM/ESFEM 耦合建立簡單矩形聲振模型，通過數(shù)值算例證明了此方法的正確性，并有較好的耦合性，以上可得出：彈性空腔應(yīng)用ESFEM 建立聲場域能夠軟化剛度矩陣，并且具有較高的計算精度，應(yīng)用性較好。

Albadr[6]提出了基于自然選擇理論的遺傳算法，與傳統(tǒng)遺傳算法相比，有較好的搜索性能，應(yīng)用不廣泛。陶猛[7]使用遺傳算法對圓柱空腔敷設(shè)吸聲材料的參數(shù)進行優(yōu)化研究，優(yōu)化后可得吸聲材料有較好的吸聲性能。陳文炯[8]建立空腔吸聲覆蓋層優(yōu)化模型，討論了材料參數(shù)對覆蓋層吸聲性能影響，基于遺傳算法搜索最優(yōu)化參數(shù)，提高覆蓋層的吸聲性能。高處[9]研究船舶艙室聲學(xué)優(yōu)化問題，基于遺傳算法對船舶艙室敷設(shè)基體材料、吸聲材料、阻尼材料的組合優(yōu)化設(shè)計，得出了降噪材料最優(yōu)布置，完成預(yù)定的降噪效果。

綜上，遺傳算法在敷設(shè)吸聲材料的各種空腔聲學(xué)優(yōu)化中，有較高的準(zhǔn)確性，搜索最優(yōu)結(jié)果效率高，魯棒性較強。因此，本文對吸聲材料的敷設(shè)方式優(yōu)化，采用穩(wěn)定性較好的遺傳算法。

1 彈性空腔聲振特性優(yōu)化模型

1.1 彈性空腔FEM/ESFEM 聲振耦合模型

彈性空腔結(jié)構(gòu)是由彈性梁與聲場域構(gòu)成。先針對彈性梁結(jié)構(gòu)采用有限元方法建立動力學(xué)方程，彈性層采用兩節(jié)點兩自由度的梁單元構(gòu)造插值函數(shù)，然后應(yīng)用變分原理，對梁單元集成，得出彈性梁的運動控制方程：

式中，Kb為梁結(jié)構(gòu)的剛度矩陣，Mb為梁結(jié)構(gòu)的質(zhì)量矩陣，ub為位移向量，F(xiàn)b為外激勵向量，F(xiàn)fb為聲壓耦合力。

結(jié)合二維聲場整體的光滑剛度矩陣Kf，并利用節(jié)點聲壓變分δP的任意性，可以得到二維聲場域光滑有限元動力學(xué)方程：

式中，Kf為二維聲場光滑剛度矩陣，Mf為二維聲場質(zhì)量矩陣，Cf為吸聲材料邊界處產(chǎn)生的阻尼矩陣，F(xiàn)f為二維空腔在邊界處的振動所引起的激勵力向量。

結(jié)合式（1）和（2），可獲得彈性空腔FEM/ESFEM聲振耦合運動控制方程：

式中ρω2HTu耦合力

FEM/ESFEM 耦合算法的二維彈性空腔系統(tǒng)方程采用MTALAB 軟件進行編程求解，為了直接有效的評價復(fù)雜彈性空腔內(nèi)部某一點的聲學(xué)性能，計算出的最大n個聲壓級（后文將具體闡述n的合理取值），再取平均值，定義為平均聲壓級為：

1.2 基于遺傳算法的彈性空腔聲振特性優(yōu)化模型

針對本研究的模型進行聲振特性研究，分析彈性空腔內(nèi)部某一參考點的聲壓級隨吸聲材料敷設(shè)方式變化的最優(yōu)解，選取設(shè)計變量為可敷設(shè)吸聲材料位置處的單元“生死”判斷變量為設(shè)計變量，將敷設(shè)吸聲材料的復(fù)雜彈性空腔聲學(xué)優(yōu)化模型描述如下：

其中，X為設(shè)計向量，表示吸聲材料在彈性空腔上可敷設(shè)吸聲材料處單元的生死變量；n為該類單元總數(shù)；ASPL表示為優(yōu)化目標(biāo)函數(shù)，即一定頻率范圍內(nèi)，復(fù)雜彈性空腔內(nèi)部聲壓級頻響曲線中最大的k個峰值的平均值，用來評價彈性空腔某一參考點的聲壓水平。

針對（5）式給出的優(yōu)化模型，采用遺傳算法。其編碼為：初始種群設(shè)計，適應(yīng)度計算，遺傳操作（包括選擇、交叉、變異）逐一解釋說明。

（1）個體編碼：擬用應(yīng)用性較好的二進制編碼對本文吸聲材料在彈性空腔上的位置上是否敷設(shè)情況進行表示，即符號0 代表不敷設(shè)吸聲材料，符號1 代表敷設(shè)吸聲材料。若有n個可敷設(shè)單元，則個體編碼長度為n，每個個體均由n個長度為1 的二進制編碼構(gòu)成。

（2）初始群體的產(chǎn)生：表示起始搜索點的初始群體。即（5）式中設(shè)計變量X的選擇，為了后文遺傳算法的進行，需選擇M個初始變量構(gòu)成種群，其中這M個初始變量可根據(jù)不同敷設(shè)方式隨意選取。后續(xù)算例中，為了遺傳算法的簡單實施，我們選擇M= 4。

（3）計算適應(yīng)度：遺傳算法中將以適應(yīng)度值的大小評價個體優(yōu)劣程度。針對上一步驟選取種群，我們采用公式（5）分別計算M個種群個體的聲壓級函數(shù)（ASPLi），作為適應(yīng)度相關(guān)函數(shù)，用以評判不同種群的適應(yīng)度，ASPLi越小，其適應(yīng)能力越強。

（4）選擇運算：選擇的目的是按照一定的評判標(biāo)準(zhǔn)將適應(yīng)能力強的個體遺傳到下一代。一般要求適應(yīng)度較高的個體將有更多的機會遺傳到下一代群體中。本文中，由于要求ASPLi越小個體適應(yīng)能力越強，這跟常用的遺傳算法要求適應(yīng)度函數(shù)越高適應(yīng)能力越強有所不同。在通用的遺傳算法中，直接采用適應(yīng)度函數(shù)在所有個體適應(yīng)度函數(shù)的總和中所占比例，來計算其適應(yīng)度。這里，每個個體的適應(yīng)度我們采用式（6）進行計算：

這樣，可以使得所有個體的適應(yīng)度總和為1，而且ASPLi越小其適應(yīng)度pi越大，滿足選擇需求，從而實現(xiàn)選擇，其具體操作過程是：

1）計算出種群中所有個體的ASPLi（i= 1.2，…，M）；

2）計算出每個個體的適應(yīng)度大小pi，作為它被遺傳到下一代種群中的概率；

3）每個概率值組成一個區(qū)域，全部pi之和為1；

4）進行（M-1）次隨機運算，生成一個0 到1 之間的數(shù)，采用輪盤賭方式，依據(jù)該隨機數(shù)出現(xiàn)在上述哪一個概率區(qū)域內(nèi)來確定各個個體被選中的次數(shù)，從而選擇性能較優(yōu)的個體進行保留，性能不佳的個體進行剔除，然后進行后續(xù)交叉運算。

（5）交叉運算：交叉操作是從種群中隨機選擇的兩個個體進行交換重組，產(chǎn)生優(yōu)秀的個體，采用應(yīng)用最為廣泛、便于理解的是單點交叉。

（6）變異運算：變異的種類較多，本研究采用二進制編碼，運用單點變異稱位點變異，通過位點變異之后，產(chǎn)生的新一個個體就會出現(xiàn)0 變?yōu)?、相反1 變?yōu)?。

基于遺傳算法的彈性空腔聲學(xué)振動優(yōu)化過程，設(shè)計了遺傳算法的機理圖，如圖1 所示。

圖1 遺傳算法機理圖

2 數(shù)值算例

2.1 某汽車駕駛室FEM/ESFEM 聲振耦合分析

使用彈性梁有限元法建模、聲腔域光滑有限元法建模，根據(jù)（3）式得出FEM/ESFEM 耦合算法的二維復(fù)雜彈性空腔的運動學(xué)控制方程，采用MATLAB 軟件編程實現(xiàn)復(fù)雜彈性空腔的聲振特性分析再進行新耦合方式合理性驗證。以某汽車駕駛室空腔（圖2）作為具體的算例來分析該類彈性空腔的聲振耦合特性。汽車駕駛室空腔尺寸參數(shù)不變，底盤視為彈性梁，彈性梁的厚度為hb= 0.003 mm、彈性模量為Eb=20.6 × 1010Pa、泊松比為ub= 0.3、密度為ρb= 7800 kg/m3，假定流體域為空氣，即密度為ρ= 1.29 kg/m3，聲速為c= 340 m/s 對此進行分析的FEM/ESFEM 耦合算法的二維駕駛室空腔模態(tài)見表1。

表1 FEM/ESFEM 耦合算法的二維駕駛室空腔前二十階固有頻率對比

圖2 某汽車二維駕駛室空腔示意圖

從表1 中可知，針對二維駕駛室聲固耦合模型分析得出的前二十階固有頻率結(jié)果，分析可得，本文方法和網(wǎng)格背景下的ANSYS 所得固有頻率結(jié)果與參考解都比較相近，固有頻率誤差值都小于3%，符合誤差要求，基于同一種網(wǎng)格的三種計算結(jié)果相對于參考解而言，在1～14 階頻率處精度不相上下，說明了本文使用的耦合算法具有正確性；本文方法誤差值小于ANSYS 和FEM 值，進一步得出了ESFEM 方法比傳統(tǒng)有限元法精度高。從15 階頻率結(jié)果開始，可以明顯看出，ESFEM 方法所得固有頻率結(jié)果更接近精細(xì)網(wǎng)格所得參考解。說明ESFEM 方法即是在處理耦合問題方面，也比通常有限元法精度高。而聲固耦合問題中，F(xiàn)EM 法由于采用低階形函數(shù)插值，很難在高頻問題中得到精確的解，在本文中可以采用ESFEM 方法，得到一定程度的彌補。

2.2 某汽車駕駛室空腔聲振特性優(yōu)化分析

上一節(jié)汽車駕駛室空腔聲振特性分析數(shù)值算例驗證了FEM/ESFEM 耦合算法的正確性。針對汽車駕駛室空腔內(nèi)部的噪聲問題，采取在駕駛室空腔頂蓋處敷設(shè)吸聲材料進行減振降噪處理，但為了減少吸聲材料的使用量和制造成本，通過MATLAB 編程實現(xiàn)遺傳算法對吸聲材料最優(yōu)布局的設(shè)計，自動搜索出了聲壓級的最優(yōu)解。

根據(jù)圖2 并結(jié)合底盤的材料屬性為鋁材，鋁材的彈性模量為Ea= 8 × 1010Pa、泊松比為ua= 0.35、密度為ρ= 2700 kg/m3，頂蓋處全敷設(shè)吸聲材料，擋風(fēng)玻璃處和其他部分進行剛性化處理，并且在發(fā)動機正下方邊界處施加法向振動速度vn= 0.001 m/s。考慮吸聲材料的位置布局設(shè)計進行優(yōu)化，通過設(shè)置30 個吸聲材料單元的死活來表示是否敷設(shè)，其敷設(shè)吸聲材料為活單位表示為1，未敷設(shè)吸聲材料為死單位表示為0，設(shè)定駕駛室員耳旁聲壓級最小值為目標(biāo)函數(shù)，以遺傳算法進行聲壓級最優(yōu)值搜索，初始種群取4 個?？紤]本文采用算例，我們通過迭代80 次，花費時間28 h，繪制出如圖3 遺傳進化曲線。

圖3 80 代的遺傳進化曲線

由圖3 可以得出，經(jīng)過遺傳進化次數(shù)為50 代時，ASPL 逐漸收斂，當(dāng)種群接近80 代時，最小目標(biāo)函數(shù)值已經(jīng)收斂。因此我們選擇第80 代的結(jié)果作為最優(yōu)解。與初始種群對比，目標(biāo)函數(shù)值A(chǔ)SPL 從初始種群中最小個體適應(yīng)度函數(shù)值112.7dB 下降到110.6 dB以下；其中優(yōu)化種群個體1 的ASPL 為108.54 dB、個體2 的ASPL 為109.77 dB、個體3 的ASPL 為110.61 dB，個體4 由于性能相對這三個個體而言較差，這里略去不表述。

圖4 給出第80 代的最終3 個選擇保留個體的SPL 頻譜曲線，通過相互對比可知，個體1 在0 ～1000 Hz 頻段內(nèi)的整體減振降噪效果最佳，其ASPL值也為最低。

圖4 80 代的最優(yōu)種群

針對初始種群設(shè)定計算出的聲壓級為優(yōu)化前，并且與遺傳搜索出來的種群優(yōu)化個體1 獲得的聲壓級為優(yōu)化后進行比較，如圖5 所示，優(yōu)化后的ASPL 值降低了1.97 dB。在頂蓋處吸聲材料全敷設(shè)的情況下獲得聲壓級曲線與優(yōu)化后的聲壓級曲線相比，全敷設(shè)和優(yōu)化后的曲線變化趨勢基本一致，ASPL 值降低了0.35 dB，并且全敷設(shè)為30 個活單元，個體1 為21 個活單元，其具體敷設(shè)方式如圖6 所示。從單元吸聲材料使用數(shù)量上來看，降低為9 個敷設(shè)單元；全敷設(shè)和優(yōu)化后的聲壓級曲線對比，如圖7 所示，曲線變化趨勢基本相似，相對于全敷設(shè)ASPL 值，優(yōu)化后的ASPL值降低了0.35 dB。經(jīng)過遺傳算法對吸聲材料位置布局優(yōu)化設(shè)計后，聲壓級得以下降，吸聲材料的使用量下降，具有良好的降噪效果，制造成本降低。

圖5 優(yōu)化前后駕駛室耳旁聲壓級曲線對比

圖6 吸聲材料的最優(yōu)布局

圖7 全敷設(shè)與優(yōu)化后聲壓級曲線對比

3 結(jié)語

采用耦合性能較好，計算精度較高的FEM/ESFEM 法建立彈性空腔聲振特性模型，再進一步考慮提高彈性空腔的減振降噪性能和敷設(shè)吸聲材料的使用量，利用收斂性強，搜索速度快的遺傳算法對吸聲材料的布局進行優(yōu)化。通過二維駕駛室空腔數(shù)值算例分析得出以下結(jié)論。

（1）仿真分析彈性空腔FEM/ESFEM 聲振耦合模型，計算精度較高、耦合性能較好。

（2）通過遺傳算法優(yōu)化出二維汽車駕駛室空腔頂蓋處吸聲材料最優(yōu)布局，減振降噪效果明顯，吸聲材料使用量降低。