彭卓凱，鄧志強(qiáng)，黃振邦

（大長(zhǎng)江集團(tuán) 研發(fā)中心，廣東 江門(mén) 529030）

空濾器結(jié)構(gòu)優(yōu)化和輻射噪聲分析

彭卓凱，鄧志強(qiáng)，黃振邦

（大長(zhǎng)江集團(tuán) 研發(fā)中心，廣東 江門(mén) 529030）

針對(duì)試驗(yàn)中發(fā)現(xiàn)的空濾器殼體輻射噪聲大問(wèn)題，綜合運(yùn)用模態(tài)計(jì)算、進(jìn)氣壓力頻譜計(jì)算、拓?fù)鋬?yōu)化、形貌優(yōu)化和聲學(xué)有限元法解決了該問(wèn)題。首先，通過(guò)模態(tài)計(jì)算，得出了空濾器外殼的各階振型和模態(tài)頻率。接著分析了發(fā)動(dòng)機(jī)的進(jìn)氣壓力波頻譜，找出了進(jìn)氣壓力較大的頻段，由此可知空濾器輻射噪聲大是由于第1階頻率較低引起的。然后，通過(guò)拓?fù)鋬?yōu)化和形貌優(yōu)化，找出了筋的最佳布置位置，使空濾器外殼的第1階頻率得到大幅提高。最后，通過(guò)輻射噪聲的計(jì)算，確認(rèn)了改進(jìn)后的空濾器外殼的總聲功率級(jí)降低了13.8 dB，噪聲改進(jìn)效果非常明顯。

聲學(xué)；空濾器；進(jìn)氣壓力；優(yōu)化；輻射噪聲；有限元法

摩托車(chē)空濾器是發(fā)動(dòng)機(jī)進(jìn)氣系統(tǒng)的一個(gè)重要組成部分，其作用是預(yù)先對(duì)進(jìn)入氣缸的空氣過(guò)濾，防止灰塵或沙粒進(jìn)入氣缸，以減少氣缸、活塞、活塞環(huán)的磨損，延長(zhǎng)發(fā)動(dòng)機(jī)的工作壽命。在現(xiàn)代摩托車(chē)的空濾器結(jié)構(gòu)設(shè)計(jì)中，還注重進(jìn)氣消聲；另外，空濾器還具有優(yōu)化匹配發(fā)動(dòng)機(jī)性能的作用[1,2]。

空濾器大體由兩部分組成：空濾器殼體和濾芯。為了達(dá)到較好的濾清效果和減少過(guò)濾時(shí)對(duì)空氣的阻力以及增大儲(chǔ)灰能力，空濾器應(yīng)設(shè)計(jì)得盡量大一些。但當(dāng)空濾器較大時(shí)，在進(jìn)氣壓力波的激勵(lì)下，空濾器殼體容易發(fā)生共振而發(fā)出較大的噪聲。某空濾器在試驗(yàn)過(guò)程中被發(fā)現(xiàn)其外殼體輻射噪聲較大，此時(shí)，發(fā)動(dòng)機(jī)轉(zhuǎn)速在6 000 r/min左右?，F(xiàn)本文綜合運(yùn)用模態(tài)計(jì)算、進(jìn)氣壓力頻譜計(jì)算、拓?fù)鋬?yōu)化、形貌優(yōu)化和聲學(xué)有限元法來(lái)解決該問(wèn)題。

1 空濾器模態(tài)計(jì)算

試驗(yàn)時(shí)出現(xiàn)共振的空濾器外殼體如圖1、2所示，該件通過(guò)6個(gè)安裝點(diǎn)與空濾器內(nèi)殼體組裝在一起。為解決該共振問(wèn)題，首先對(duì)空濾器外殼體進(jìn)行模態(tài)分析。將在NX軟件中建好的3D模型導(dǎo)入有限元分析軟件Nastran中，進(jìn)行網(wǎng)格劃分，采用2 mm二階四面體單元。接著定義部件的材料屬性，材料為PP塑料，彈性模量為1 000 MPa，泊松比為0.41，密度為890 kg/m3。再定義邊界條件，將6個(gè)安裝孔約束住。最后定義模態(tài)提取分析步，提取1 000 Hz以前的所有模態(tài)。計(jì)算出的原始空濾器外殼的模態(tài)頻率見(jiàn)表1，第1階模態(tài)振型見(jiàn)圖3，從該振型可知，中部的圓形區(qū)域振幅最大。從計(jì)算出的模態(tài)頻率和振型看，第1階模態(tài)頻率較低是造成輻射噪聲大的主要原因，因?yàn)樵谶M(jìn)氣壓力波的激勵(lì)下，空濾器外殼會(huì)發(fā)生共振。為進(jìn)一步驗(yàn)證該推斷，下面進(jìn)行進(jìn)氣壓力頻譜計(jì)算。

圖1 原始空濾器外殼（里側(cè)）

圖2 原始空濾器外殼（外側(cè)）

圖3 原始空濾器外殼第1階振型

表1 改進(jìn)前、后的空濾器外殼模態(tài)頻率

2 進(jìn)氣壓力頻譜計(jì)算

該發(fā)動(dòng)機(jī)為125 ml排量的單缸自然風(fēng)冷發(fā)動(dòng)機(jī)。在發(fā)動(dòng)機(jī)模擬分析軟件boost中，對(duì)該發(fā)動(dòng)機(jī)進(jìn)行建模，設(shè)置相應(yīng)的進(jìn)氣系統(tǒng)、燃燒系統(tǒng)、排氣系統(tǒng)等參數(shù)，在空濾器內(nèi)設(shè)置壓力監(jiān)測(cè)點(diǎn)，測(cè)量發(fā)動(dòng)機(jī)轉(zhuǎn)速在6 000 r/min時(shí)空濾器內(nèi)壓力隨曲軸轉(zhuǎn)角的變化關(guān)系，并將此壓力在時(shí)域內(nèi)的信號(hào)通過(guò)快速傅里葉變換（FFT）轉(zhuǎn)換為頻域信號(hào)[3]。在6 000 r/min時(shí)，空濾器內(nèi)壓力的頻譜如圖4所示。

圖4 空濾器內(nèi)壓力在6 000 r/min時(shí)的頻譜

從該圖可知，在300 Hz和低于300 Hz的區(qū)域，空濾器內(nèi)壓力較大，在高于300 Hz的區(qū)域，空濾器內(nèi)壓力隨頻率增加而逐漸減小。從上面的模態(tài)計(jì)算可知，空濾器外殼體第1階模態(tài)頻率為312 Hz，這個(gè)頻率位于進(jìn)氣壓力較大的頻域范圍附近，故而會(huì)出現(xiàn)共振而產(chǎn)生較大的輻射噪聲。

3 拓?fù)鋬?yōu)化和形貌優(yōu)化

優(yōu)化設(shè)計(jì)有三要素，即設(shè)計(jì)變量、目標(biāo)函數(shù)、約束條件。設(shè)計(jì)變量是在優(yōu)化過(guò)程中發(fā)生改變從而提高性能的一組參數(shù)。目標(biāo)函數(shù)就是要求的最優(yōu)性能，是關(guān)于設(shè)計(jì)變量的函數(shù)。約束條件是對(duì)設(shè)計(jì)的限制。目標(biāo)函數(shù)、約束函數(shù)是從有限元分析中獲得的結(jié)構(gòu)響應(yīng)。設(shè)計(jì)變量的選擇依賴于優(yōu)化類(lèi)型。在拓?fù)鋬?yōu)化中，設(shè)計(jì)變量是單元的密度，在形貌優(yōu)化中，設(shè)計(jì)變量為形狀擾動(dòng)的線性組合因子[4]。拓?fù)鋬?yōu)化能在給定的設(shè)計(jì)空間內(nèi)尋求最佳的材料分布。形貌優(yōu)化是一種形狀最佳化的方法，即在板形結(jié)構(gòu)中尋找最優(yōu)的凹凸型加強(qiáng)筋的分布，與拓?fù)鋬?yōu)化不同的是，形貌優(yōu)化不刪除材料。

從以上的模態(tài)和進(jìn)氣壓力頻譜分析的結(jié)果可知，空濾器外殼體第1階模態(tài)頻率低，從而出現(xiàn)了共振問(wèn)題。故提高第1階模態(tài)頻率可解決該問(wèn)題。從原始空濾器外殼體的3D模型看（見(jiàn)圖1、2），其上面布置了大量的網(wǎng)狀加強(qiáng)筋，且還帶有3個(gè)四邊形的凹凸筋。從提高模態(tài)頻率的角度考慮，這些加強(qiáng)筋和凹凸筋的形狀和位置未必是合理的。為了找到最佳的加強(qiáng)筋和凹凸筋的布置位置，現(xiàn)對(duì)空濾器外殼體進(jìn)行拓?fù)鋬?yōu)化和形貌優(yōu)化，以提高其第1階模態(tài)頻率。其中，拓?fù)鋬?yōu)化用于確定加強(qiáng)筋（增加材料的筋）的最佳布置，形貌優(yōu)化用于確定凹凸筋（不增加材料的筋）的最佳布置。

將原始的加強(qiáng)筋和凹凸筋都去除，再用優(yōu)化計(jì)算重新布置這些筋。優(yōu)化計(jì)算的模型采用1 mm四邊形殼單元，如圖5所示，其中灰色區(qū)域?yàn)樵O(shè)計(jì)區(qū)域（可優(yōu)化區(qū)域），邊緣黑色區(qū)域?yàn)榉窃O(shè)計(jì)區(qū)域。優(yōu)化目標(biāo)為第1階模態(tài)頻率最大。在拓?fù)鋬?yōu)化計(jì)算中，加強(qiáng)筋高度設(shè)為3.5 mm。在形貌優(yōu)化計(jì)算中，凹凸筋的深度設(shè)為3 mm，起筋角度設(shè)為60°。

圖5 空濾器外殼優(yōu)化計(jì)算模型

形貌優(yōu)化計(jì)算結(jié)果見(jiàn)圖6，其中，除黑色以外的白色區(qū)域和灰色區(qū)域?yàn)樽罴训陌纪菇畈贾梦恢?。從圖可知，凹凸筋應(yīng)布置在外殼的中部，呈近似橢圓形。而原始空濾器外殼的凹凸筋為間隔一定距離的3個(gè)四邊形，與優(yōu)化結(jié)果有較大差異，這表明，原始凹凸筋的設(shè)計(jì)不盡合理。

圖6 空濾器外殼形貌優(yōu)化結(jié)果

拓?fù)鋬?yōu)化計(jì)算結(jié)果見(jiàn)圖7，其中，除黑色以外的灰色區(qū)域和接近灰色的白色區(qū)域?yàn)樽罴训募訌?qiáng)筋布置位置。從圖可知，加強(qiáng)筋應(yīng)布置在外殼的邊緣，中間的大部分區(qū)域都無(wú)需布置加強(qiáng)筋。而原始空濾器外殼的加強(qiáng)筋呈網(wǎng)狀密布，與優(yōu)化結(jié)果差異大，這表明，原始加強(qiáng)筋的設(shè)計(jì)既浪費(fèi)材料也不利于提高模態(tài)頻率。

圖7 空濾器外殼拓?fù)鋬?yōu)化結(jié)果

現(xiàn)根據(jù)拓?fù)鋬?yōu)化和形貌優(yōu)化計(jì)算結(jié)果，對(duì)空濾器外殼的加強(qiáng)筋和凹凸筋重新設(shè)計(jì)，重新設(shè)計(jì)后的3D模型如圖8所示。

圖8 重新設(shè)計(jì)后的空濾器外殼

4 空濾器輻射噪聲計(jì)算

4.1 模態(tài)和動(dòng)態(tài)響應(yīng)計(jì)算

為進(jìn)一步確認(rèn)改進(jìn)前、后的空濾器外殼的聲學(xué)響應(yīng)特性，現(xiàn)進(jìn)行輻射噪聲的分析?？諡V器外殼在進(jìn)氣壓力波的激勵(lì)下，表面會(huì)產(chǎn)生振動(dòng)，從而發(fā)出噪聲。故首先須計(jì)算空濾器外殼表面的振動(dòng)響應(yīng)。對(duì)于帶阻尼的多自由度系統(tǒng)，其受迫振動(dòng)的微分方程為

解這個(gè)運(yùn)動(dòng)方程一般有兩種方法，一種是直接積分法，就是按時(shí)間歷程對(duì)上述微分方程直接進(jìn)行數(shù)值積分；另一種解法是模態(tài)疊加法[5,6]。對(duì)于計(jì)算系統(tǒng)在諧波激勵(lì)下的線性響應(yīng)，后一種基于模態(tài)的計(jì)算方法比直接積分的方法計(jì)算速度更快。故本文采用模態(tài)疊加法計(jì)算空濾器外殼表面的振動(dòng)速度。首先須計(jì)算空濾器外殼的模態(tài)，再基于模態(tài)結(jié)果計(jì)算空濾器外殼在進(jìn)氣壓力波激勵(lì)下的振動(dòng)速度。改進(jìn)后空濾器外殼的結(jié)構(gòu)有限元模型見(jiàn)圖9，采用2.5 mm二階四面體單元。改進(jìn)前、后的空濾器外殼模態(tài)頻率見(jiàn)表1，從該表可知，在不增加空濾器外殼基本壁厚的前提下，僅僅改變筋的布置，其第1階模態(tài)頻率已經(jīng)由312 Hz大幅提高到了487 Hz。

圖9 改進(jìn)后空濾器外殼的結(jié)構(gòu)有限元模型

4.2 輻射噪聲計(jì)算

為求解流場(chǎng)中任意位置（x,y,z）處的聲壓p，可通過(guò)流體的連續(xù)方程、運(yùn)動(dòng)方程、能量方程、物態(tài)方程推導(dǎo)出Helmholtz聲學(xué)波動(dòng)方程，其表達(dá)式為[7]

?2p(x,y,z)-k2p(x,y,z)=-jρ0ωq(x,y,z)

式中k為波數(shù)，ω為角頻率，ρ0為定常流動(dòng)下的密度，q體積速度。用數(shù)值方法求解Helmholtz方程的常用方法有聲學(xué)有限元法（FEM）和聲學(xué)邊界元法（BEM）。當(dāng)模型不是特別大時(shí)，有限元法在計(jì)算內(nèi)聲場(chǎng)和輻射噪聲方面的速度更快，故本文采用有限元法。用有限元法求解Helmholtz方程時(shí)，首先需要把計(jì)算的聲場(chǎng)離散成一定數(shù)量的小聲場(chǎng)，每個(gè)小聲場(chǎng)稱(chēng)為單元，單元之間用節(jié)點(diǎn)連接。對(duì)于三維聲場(chǎng)，常用的單元有四面體單元和六面體單元。

本文采用3 mm四面體單元對(duì)空濾器外殼周?chē)穆晥?chǎng)進(jìn)行離散，聲學(xué)網(wǎng)格模型的截面圖見(jiàn)圖10。

并以空濾器外殼為中心布置半徑為1 m的球面場(chǎng)點(diǎn)（見(jiàn)圖11），以計(jì)算該球面場(chǎng)點(diǎn)上的輻射聲功率。將上一步計(jì)算的空濾器外殼的表面振動(dòng)速度導(dǎo)入聲學(xué)計(jì)算軟件LMS Virtual Lab，將表面振動(dòng)速度映射到聲學(xué)網(wǎng)格上作為輻射噪聲計(jì)算的邊界條件，求得改進(jìn)前、后的空濾器外殼的輻射聲功率如圖12所示。從計(jì)算結(jié)果可知，原始空濾器外殼在312 Hz和739 Hz附近出現(xiàn)了聲功率峰，312 Hz處的最大聲功率級(jí)為127.2 dB。新設(shè)計(jì)的空濾器外殼在492 Hz和882 Hz附近出現(xiàn)了聲功率峰，492 Hz處的最大聲功率級(jí)為112.5 dB。對(duì)比模態(tài)計(jì)算結(jié)果可知，改進(jìn)前、后的空濾器外殼的聲功率峰都出現(xiàn)在第1階和第4階模態(tài)頻率附近。原始空濾器外殼在整個(gè)頻段的總聲功率級(jí)為133 dB，新設(shè)計(jì)的空濾器外殼在整個(gè)頻段的總聲功率級(jí)為119.2 dB?？偮暪β始?jí)降低了13.8 dB；最大聲功率級(jí)降低了14.7 dB。從以上輻射噪聲的計(jì)算結(jié)果可知，噪聲改進(jìn)效果非常明顯。

圖11 球面場(chǎng)點(diǎn)

圖12 改進(jìn)前、后的空濾器外殼的輻射聲功率

5 結(jié)語(yǔ)

（1）通過(guò)模態(tài)計(jì)算和進(jìn)氣壓力波計(jì)算，得出了空濾器外殼的各階振型和模態(tài)頻率以及進(jìn)氣壓力較大的頻段，找出了空濾器外殼輻射噪聲大的原因；

（2）通過(guò)拓?fù)鋬?yōu)化和形貌優(yōu)化，找出了筋的最佳布置位置，使空濾器外殼的第1階頻率得到大幅提高；

（3）通過(guò)輻射噪聲的計(jì)算，進(jìn)一步確認(rèn)了空濾器外殼的改進(jìn)效果：總聲功率級(jí)降低了13.8 dB。噪聲改進(jìn)效果非常明顯。

以上分析流程及相關(guān)的模擬計(jì)算方法為快速、準(zhǔn)確地解決輻射噪聲大問(wèn)題提供了解決方案，減少了費(fèi)時(shí)費(fèi)力的試驗(yàn)，也節(jié)約了部件的材料成本（相對(duì)將部件簡(jiǎn)單的加厚的對(duì)策方案），對(duì)解決類(lèi)似的車(chē)輛部件輻射噪聲大問(wèn)題具有一定的參考價(jià)值。

[1]馮超，鄔惠樂(lè).汽車(chē)工程手冊(cè)摩托車(chē)篇[M].北京：人民交通出版社，2001.330-331.

[2]曾憲錚，吳振根.中外流行摩托車(chē)構(gòu)造使用維修[M].長(zhǎng)沙：湖南科學(xué)技術(shù)出版社，1996.164-165.

[3]朱訪君，吳堅(jiān).內(nèi)燃機(jī)工作過(guò)程數(shù)值計(jì)算及其優(yōu)化[M].北京：國(guó)防工業(yè)出版社，1997.89-143.

[4]張勝蘭，鄭冬黎.基于Hyper Works的結(jié)構(gòu)優(yōu)化設(shè)計(jì)技術(shù)[M].北京：機(jī)械工業(yè)出版社，2007.159-163.

[5]許本文，焦群英.機(jī)械振動(dòng)與模態(tài)分析基礎(chǔ)[M].北京：機(jī)械工業(yè)出版社，1998.70-78

[6]劉延柱，陳文良，陳立群.振動(dòng)力學(xué)[M].北京：高等教育出版社，1998.85-90.

[7]李增剛，詹福良.Virtual.Lab Acoustics聲學(xué)仿真計(jì)算高級(jí)應(yīng)用實(shí)例[M].北京：國(guó)防工業(yè)出版社，2010.4-17.

Structure Optimization and Radiation NoiseAnalysis of Air Cleaners

PENG Zhuo-kai,DENG Zhi-qiang,HUANG Zhen-bang
（Research&Development Center,Dachangjiang Group,Jiangmen 529030,Guangdong China）

Modal calculation,frequency spectrum calculation of intake pressure,topology analysis,topography analysis and acoustic FEM were employed to solve the serious radiation noise problem of an air cleaner’s shell.Firstly,vibration modes and eigen-frequencies of the air cleaner shell were obtained by modal computation.Secondly,the intake pressure spectrum was analyzed and the frequency zone with large intake pressure was determined.It was found that the large radiation noise of the air cleaner was resulted from the low 1st order frequency.Finally,the optimal position of the ribs was discovered by topology and topography analysis,so that the 1st order frequency of the air cleaner could be enhanced greatly. The result of the radiation noise computation show that the total sound power level of the air cleaner can be reduced by 13.8 dB after the modification and the improvement effect is very obvious.

acoustics;air cleaner;intake pressure;optimization;radiation noise;finite element method

TB53

A

10.3969/j.issn.1006-1335.2015.02.049

1006-1355(2015)02-0222-04

2014－10－17

彭卓凱（1979－），男，碩士，主要從事摩托車(chē)及其發(fā)動(dòng)機(jī)的CAE和NVH分析。E-mail:pengzhuokai@126.com