楊　楊，陸森林

(江蘇大學(xué) 汽車與交通工程學(xué)院，江蘇 鎮(zhèn)江　212013)

?

車內(nèi)噪聲的數(shù)值化分析與降噪設(shè)計(jì)

楊楊，陸森林

(江蘇大學(xué) 汽車與交通工程學(xué)院，江蘇 鎮(zhèn)江212013)

針對(duì)一般腔體結(jié)構(gòu)的內(nèi)部聲場(chǎng)公式，利用有限元法得到聲壓響應(yīng)幅度判定參數(shù)的計(jì)算公式。以某客車車身結(jié)構(gòu)和聲腔有限元模型為例，計(jì)算得到低頻結(jié)構(gòu)車身頂棚的聲壓響應(yīng)幅度判定參數(shù)。為降低車身頂棚的聲壓響應(yīng)幅度判定參數(shù)，提出車身頂棚的改進(jìn)方案。利用聲學(xué)軟件計(jì)算改進(jìn)前后的車身車內(nèi)噪聲，并分析駕駛員和前3排乘客場(chǎng)點(diǎn)處的聲壓級(jí)曲線。結(jié)果表明：通過計(jì)算頂棚聲壓響應(yīng)幅度判定參數(shù)以降低車內(nèi)噪聲的方法是可行的。

客車車身；聲壓響應(yīng)幅度判定參數(shù)；有限元方法；聲壓級(jí)曲線

目前，車內(nèi)的噪聲水平?jīng)Q定了駕駛員和乘客的乘坐舒適性，成為評(píng)價(jià)車輛乘坐舒適性的重要指標(biāo)之一[1-4]。目前，國(guó)內(nèi)外對(duì)車內(nèi)噪聲的研究方法包括實(shí)驗(yàn)法和CAE法兩種。實(shí)驗(yàn)法通常采取鋪設(shè)隔聲吸聲材料、涂裝阻尼、添加結(jié)構(gòu)加強(qiáng)筋等措施以降低車室內(nèi)的低頻結(jié)構(gòu)噪聲。CAE法的特點(diǎn)是：在產(chǎn)品成型或制造前，通過相關(guān)軟件分析出車身結(jié)構(gòu)的薄弱環(huán)節(jié)，然后進(jìn)行結(jié)構(gòu)優(yōu)化。該方法可以縮短產(chǎn)品上市周期，節(jié)約大量生產(chǎn)成本。

有研究指出，對(duì)于一般腔體的無聲源內(nèi)部聲場(chǎng)表達(dá)式，可以通過對(duì)整個(gè)車身的聲壓響應(yīng)幅度判定參數(shù)幅值進(jìn)行分析，以達(dá)到降低車內(nèi)噪聲的目的[5]。這種方法的優(yōu)點(diǎn)在于將降噪設(shè)計(jì)和有限元法相結(jié)合，將車身的噪聲問題轉(zhuǎn)化為有限元模型節(jié)點(diǎn)的計(jì)算問題，但其計(jì)算量太大。若能對(duì)此判定參數(shù)進(jìn)行分塊處理，則可以大大減少工作量，即通過降低車身某部分的聲壓響應(yīng)幅度判定參數(shù)，以及聲壓貢獻(xiàn)值，達(dá)到整體車內(nèi)噪聲下降的效果。因此，以計(jì)算車身頂棚聲壓響應(yīng)幅度判定參數(shù)為例，提出車身頂棚改進(jìn)策略，并分析車身改進(jìn)前后7個(gè)場(chǎng)點(diǎn)處的聲壓級(jí)曲線，證實(shí)這種采用局部計(jì)算聲壓響應(yīng)幅度判定參數(shù)以達(dá)到降低車內(nèi)噪聲的思路是比較合理的。

1　車身有限元模型

在CATIA中建立客車白車身模型。模型應(yīng)能反映結(jié)構(gòu)力學(xué)特性，以及車身結(jié)構(gòu)的實(shí)際形狀和相互間的搭接關(guān)系。模型忽略尺寸較小的圓角、倒角、圓孔以及一些非承載件等[6]。將裝配完成的模型以“.IGS”格式導(dǎo)出，并由Hyper mesh軟件導(dǎo)入，對(duì)白車身進(jìn)行幾何清理。采用40 mm的shell63單元?jiǎng)澐志W(wǎng)格。車身主體的單元厚度為2 mm，如車架、頂棚加強(qiáng)筋等，其他部件單元厚度為1 mm。賦予車身材料和屬性，運(yùn)用rigid命令建立點(diǎn)焊連接，單元數(shù)量為37 177個(gè)，節(jié)點(diǎn)數(shù)量為36 286個(gè)。

將結(jié)構(gòu)有限元模型導(dǎo)入Hyper mesh中，運(yùn)用ruled命令對(duì)白車身單元進(jìn)行封閉處理，再利用Hyper mesh中的NVH模塊提取車身結(jié)構(gòu)內(nèi)壁與空氣接觸的表面構(gòu)成封閉的聲學(xué)空腔[7]，快速生成白車身的聲學(xué)模型。在劃分網(wǎng)格時(shí)，采用六面體和四面體混合單元網(wǎng)格劃法，單元的大小要求不大于計(jì)算頻率最短波長(zhǎng)的1/6[8]，即單元長(zhǎng)度L滿足式(1)。

(1)

考慮局部結(jié)構(gòu)和聲腔節(jié)點(diǎn)的對(duì)應(yīng)關(guān)系，聲腔模型的單元網(wǎng)格亦取40 mm,劃分單元16 4681個(gè)。所得結(jié)構(gòu)和聲腔有限元模型如圖1所示。

圖1　車身有限元模型

2　一般腔體內(nèi)部聲場(chǎng)公式分析

2.1一般腔體內(nèi)部聲場(chǎng)解析式

(2)

其中：p(r)為極坐標(biāo)下系統(tǒng)內(nèi)部點(diǎn)r處的聲壓；ρ0為常溫空氣密度；c0為常溫空氣聲速；S為結(jié)構(gòu)聲腔接觸的表面積；F為外界激勵(lì)力；n和p分別為聲腔和結(jié)構(gòu)的模態(tài)階數(shù)；Cnp為聲腔第n階模態(tài)和結(jié)構(gòu)第p階模態(tài)的振型耦合系數(shù)；Fnp為聲腔第n階模態(tài)和結(jié)構(gòu)第p階模態(tài)的頻率重疊系數(shù)；ψn(r)聲腔振型在響應(yīng)點(diǎn)處的分量；φp(ρ)為結(jié)構(gòu)振型在激勵(lì)點(diǎn)處的分量；Λn為聲腔第n階模態(tài)質(zhì)量參數(shù)；Λp為結(jié)構(gòu)第p階模態(tài)質(zhì)量。

將式(2)由極坐標(biāo)轉(zhuǎn)換到笛卡爾坐標(biāo)系中，有振型耦合系數(shù)和頻率重疊系數(shù)表達(dá)式如式(3)、(4)所示。

(3)

(4)

將振型耦合系數(shù)、頻率重疊系數(shù)以及激勵(lì)點(diǎn)處的結(jié)構(gòu)振型分量三者的乘積幅值作為車內(nèi)噪聲聲壓響應(yīng)幅度判定參數(shù)[9]。激勵(lì)點(diǎn)處的結(jié)構(gòu)振型分量由模態(tài)振型可得，記為phis。

2.2振型耦合系數(shù)分析

由于研究的車身形狀為不規(guī)則形，運(yùn)用有限元的方法，將振型耦合系數(shù)轉(zhuǎn)化為關(guān)于有限元模型節(jié)點(diǎn)的函數(shù)，得到振型耦合系數(shù)近似解[10]。假設(shè)式(3)中的dx=lsx/q1,dy=lsy/q2，其中q1和q2分別為兩條積分邊lsx和lsy上的總段數(shù)，并將式(3)轉(zhuǎn)化為離散的求和表達(dá)式：

(5)

利用有限元自由模態(tài)分析結(jié)果，可以得到耦合面上聲腔和結(jié)構(gòu)模態(tài)振型在各個(gè)節(jié)點(diǎn)上的分量，因此，將結(jié)構(gòu)和聲腔的模態(tài)振型函數(shù)轉(zhuǎn)換為以節(jié)點(diǎn)編號(hào)為自變量的函數(shù)：

(6)

我的身體太搖晃了，除了盡力躲開他我什么都做不了，在場(chǎng)上允許的范圍里，能躲多遠(yuǎn)就多遠(yuǎn)。他沖了過來，狠狠踢向我的腹部，這一腳把我肺里的空氣都逼了出來，很疼。也許是太疼了，也許是因?yàn)楸惶吡艘荒_，我感到無法呼吸，到底是什么原因說不清，只是倒在地上。

2.3頻率重疊系數(shù)分析

3　聲壓響應(yīng)幅度判定參數(shù)的計(jì)算

3.1振型耦合系數(shù)計(jì)算

在ANSYS中分別對(duì)車身結(jié)構(gòu)和聲腔模型進(jìn)行自由模態(tài)分析。結(jié)構(gòu)自由模態(tài)分析采用蘭索法(Block Lanczos)[11]，聲腔自由模態(tài)分析采用非對(duì)稱法。對(duì)車身結(jié)構(gòu)模型求解,導(dǎo)出頂棚耦合面上各個(gè)節(jié)點(diǎn)的前200階模態(tài)的第p階固有頻率和第p階固有振型。對(duì)車身聲腔模型求解,導(dǎo)出頂棚耦合面上各個(gè)節(jié)點(diǎn)前20階模態(tài)的第n階固有頻率和第n階固有振型在XY平面內(nèi)的振型分量。

利用式(5)，通過Matlab編程計(jì)算出振型耦合系數(shù)，其數(shù)值有正有負(fù)，正負(fù)號(hào)僅代表方向。從計(jì)算結(jié)果發(fā)現(xiàn)，對(duì)于頂棚模態(tài)振型很小的階數(shù)，其振型耦合系數(shù)也很小。圖2中僅列出40階頂棚模態(tài)振型較大的振型耦合系數(shù)幅值結(jié)果。

3.2頻率重疊系數(shù)計(jì)算

計(jì)算頻率重疊系數(shù)時(shí)，需要激勵(lì)力的頻率數(shù)值。車身在運(yùn)動(dòng)過程中承受著來自發(fā)動(dòng)機(jī)和路面等激勵(lì)作用。以發(fā)動(dòng)機(jī)激勵(lì)為例，通過分析車身頂棚在發(fā)動(dòng)機(jī)激勵(lì)下的受迫振動(dòng)響應(yīng)曲線，找到具有代表性的發(fā)動(dòng)機(jī)激勵(lì)力圓頻率點(diǎn)，此圓頻率點(diǎn)應(yīng)是引起較大車內(nèi)聲壓值的發(fā)動(dòng)機(jī)激勵(lì)頻率。在車身發(fā)動(dòng)機(jī)懸置點(diǎn)定義發(fā)動(dòng)機(jī)激勵(lì)，響應(yīng)點(diǎn)定義為車身頂棚前部和后部中點(diǎn)，運(yùn)用Virtual.Lab軟件分析得到頻率區(qū)間在20～100 Hz下的頂棚受迫振動(dòng)響應(yīng)曲線，如圖3所示。

圖2　改進(jìn)前振型耦合系數(shù)

圖3　頂棚受迫振動(dòng)響應(yīng)曲線

由圖3可見：在發(fā)動(dòng)機(jī)激勵(lì)頻率為36，78 Hz等時(shí)頂棚有較大的振動(dòng)。限于篇幅，取發(fā)動(dòng)機(jī)激勵(lì)力頻率為36 Hz，通過式(4)計(jì)算頻率重疊系數(shù)，并對(duì)這些數(shù)據(jù)取幅值，結(jié)果如圖4所示。頻率重疊系數(shù)的峰值出現(xiàn)在結(jié)構(gòu)第41階模態(tài)、聲腔第2階模態(tài)處，此處結(jié)構(gòu)模態(tài)固有頻率、聲腔模態(tài)固有頻率與發(fā)動(dòng)機(jī)激勵(lì)頻率最接近。

3.3聲壓響應(yīng)幅度判定參數(shù)計(jì)算結(jié)果

已知振型耦合系數(shù)結(jié)果、頻率重疊系數(shù)結(jié)果和結(jié)構(gòu)在激勵(lì)力點(diǎn)處的振型分量，最終得到頂棚聲壓響應(yīng)幅度判定參數(shù)，如圖5所示。從圖5中發(fā)現(xiàn)：頻率重疊系數(shù)對(duì)聲壓響應(yīng)幅度判定參數(shù)有很大影響，決定了判定參數(shù)的整體基本變化趨勢(shì)，但在判定參數(shù)峰值區(qū)域，振型耦合系數(shù)起到了關(guān)鍵作用，使得判定參數(shù)趨于無規(guī)律的變化。

圖4　改進(jìn)前頻率重疊系數(shù)

圖5　改進(jìn)前聲壓響應(yīng)幅度判定參數(shù)

4　車身頂棚改進(jìn)方案

聲壓響應(yīng)幅度判定參數(shù)由結(jié)構(gòu)振型函數(shù)、聲腔振型函數(shù)、結(jié)構(gòu)固有頻率、聲腔固有頻率以及激勵(lì)力頻率等因素決定。結(jié)構(gòu)固有頻率、聲腔固有頻率以及激勵(lì)力頻率在改進(jìn)前后的變化不大。因此，可以通過減小結(jié)構(gòu)振型函數(shù)和聲腔振型函數(shù)來降低聲壓響應(yīng)幅度判定參數(shù)的數(shù)值。考慮到車身結(jié)構(gòu)已經(jīng)成形，通過改變聲腔形狀降低聲腔振型函數(shù)的方法不容易實(shí)現(xiàn)。因此采取提高車身結(jié)構(gòu)剛度的方法來降低車內(nèi)噪聲，而改變結(jié)構(gòu)剛度可以通過在結(jié)構(gòu)壁板上增加加強(qiáng)筋來實(shí)現(xiàn)。因此，考慮改變結(jié)構(gòu)剛度以降低結(jié)構(gòu)振型函數(shù)，從而達(dá)到降低聲壓響應(yīng)幅度判定參數(shù)的目的。

具體方法：對(duì)頂棚上橫梁的數(shù)量和分布情況進(jìn)行調(diào)整。如圖6所示，將頂棚左端起第2根梁前移，第5根梁后移，并在第2、3根梁中間增加1根梁，第5、6根梁中間增加1根梁，使得改進(jìn)后的前4根梁間距、后4根梁間距相等。

圖6　改進(jìn)前后的車身頂棚

對(duì)改進(jìn)后結(jié)構(gòu)進(jìn)行模態(tài)分析，得到結(jié)構(gòu)前200階模態(tài)固有頻率和振型。同樣計(jì)算得到改進(jìn)后40階頂棚振型較大的振型耦合系數(shù)值、發(fā)動(dòng)機(jī)激勵(lì)頻率為36 Hz時(shí)的頻率重疊系數(shù)以及聲壓響應(yīng)幅度判定參數(shù)，見圖7。

通過與改進(jìn)前各參數(shù)的比較，改進(jìn)后頂棚的振型耦合系數(shù)有所下降，而頻率重疊系數(shù)有所變化，因此聲壓響應(yīng)幅度判定參數(shù)也有所下降。

5　改進(jìn)前后的噪聲分析

在聲學(xué)軟件Virtual.Lab中進(jìn)行發(fā)動(dòng)機(jī)激勵(lì)下的耦合車內(nèi)噪聲分析，車身結(jié)構(gòu)和聲腔模型應(yīng)含有模態(tài)信息。為降低車身內(nèi)部聲場(chǎng)，在車身結(jié)構(gòu)上定義阻尼層以及多孔材料參數(shù)，在聲腔包絡(luò)面上定義吸聲系數(shù)。場(chǎng)點(diǎn)采用駕駛員和6位乘客的右耳耳旁網(wǎng)格，聲壓級(jí)采用A計(jì)權(quán)標(biāo)準(zhǔn)。分析得到了頂棚改進(jìn)前和改進(jìn)后各場(chǎng)點(diǎn)的聲壓級(jí)曲線,如圖8所示。

從聲壓級(jí)曲線進(jìn)行分析，得到各頻率下聲壓級(jí)的變化：從聲壓級(jí)曲線的局部進(jìn)行分析，改進(jìn)后各場(chǎng)點(diǎn)處的峰值都有下降，其中在前排乘客右側(cè)場(chǎng)點(diǎn)處最大下降了6.16 dB；但在80Hz附近，中排左側(cè)乘客場(chǎng)點(diǎn)處的聲壓值略有增加；車內(nèi)聲壓較大值多出現(xiàn)在頻率40 Hz附近和頻率80 Hz附近，與圖2中頂棚的受迫振動(dòng)響應(yīng)曲線在這2個(gè)頻率點(diǎn)附近有較大振幅的情況相吻合。從聲壓級(jí)曲線總體上分析，在聲壓值較小的頻率區(qū)間里，改進(jìn)前后的聲壓級(jí)數(shù)值有升高也有降低，這是因?yàn)楦倪M(jìn)前后車身的固有頻率發(fā)生了變化，聲壓級(jí)曲線的增減趨勢(shì)也發(fā)生了改變，但改進(jìn)后聲壓級(jí)曲線多處于改進(jìn)前曲線下方。

圖7　改進(jìn)后各參數(shù)結(jié)果

圖8　改進(jìn)前后場(chǎng)點(diǎn)處聲壓級(jí)曲線

6　結(jié)束語

對(duì)于一般腔體結(jié)構(gòu)，可以用聲壓響應(yīng)幅度判定參數(shù)來表征其內(nèi)部聲場(chǎng)情況，計(jì)算聲壓響應(yīng)幅度判定參數(shù)可以達(dá)到對(duì)內(nèi)部聲場(chǎng)的數(shù)值化分析效果，直觀地反映出車內(nèi)聲場(chǎng)與聲固耦合模型模態(tài)之間的內(nèi)在聯(lián)系。但計(jì)算整車的聲壓響應(yīng)幅度判定參數(shù)計(jì)算量過大，故考慮只計(jì)算車身某部的判定參數(shù)。

通過車身局部部件的聲壓響應(yīng)幅度判定參數(shù)有限元計(jì)算分析車身聲固耦合作用，減小了計(jì)算的工作量。通過改進(jìn)頂棚結(jié)構(gòu)，最終降低了頂棚的聲壓響應(yīng)幅度判定參數(shù)幅值。

從聲壓級(jí)曲線的分析中可以發(fā)現(xiàn)：大部分聲壓級(jí)都有所降低，個(gè)別頻率點(diǎn)上的聲壓級(jí)略有升高。說明對(duì)車身頂棚的改進(jìn)方案可行，也說明在本例中，通過計(jì)算車身局部部件的聲壓響應(yīng)幅度判定參數(shù)分析車內(nèi)噪聲情況具有一定的理論指導(dǎo)作用。

[1]王彥博，陳劍.汽車內(nèi)部聲場(chǎng)分析[J].噪聲與振動(dòng)控制，2009(2)：84-87.

[2]胡遠(yuǎn)志, 雷雨, 鄭光澤,等. 車身密封性對(duì)車內(nèi)噪聲的影響[J]. 重慶理工大學(xué)學(xué)報(bào)(自然科學(xué)), 2014(9):1-6.

[3]夏勇, 鄒波, 賀才春,等. 某城市客車車內(nèi)噪聲試驗(yàn)與治理[J]. 客車技術(shù)與研究, 2015(3):59-62.

[4]康強(qiáng), 左曙光, 韓惠君. 汽車空調(diào)系統(tǒng)離心風(fēng)機(jī)氣動(dòng)噪聲數(shù)值計(jì)算[J]. 江蘇大學(xué)學(xué)報(bào)(自然科學(xué)版), 2013, 34(1):1-6.

[5]高書娜.轎車車身低頻結(jié)構(gòu)聲耦合機(jī)理及降噪設(shè)計(jì)研究[D].重慶：重慶大學(xué)，2011.

[6]唐維偉，何仁，游專，等.純電動(dòng)客車車身優(yōu)化設(shè)計(jì)[J].機(jī)械設(shè)計(jì)與制造，2013(12)：20-23.

[7]馮蘭芳，夏兆義，王宏曉，等.微型客車車身結(jié)構(gòu)模態(tài)與聲腔模態(tài)分析[J].機(jī)械研究與應(yīng)用，2012(4)：27-28.

[8]劉成武，黃鼎鍵，鐘勇.基于NVH的車內(nèi)聲腔模態(tài)分析[J].福建工程學(xué)院學(xué)報(bào)，2009(4)：347-350.

[9]高書娜，鄧兆祥.轎車車內(nèi)低頻噪聲的判定參數(shù)探討[J].汽車工程，2013(8):744-748.

[10]鄧兆祥，高書娜.結(jié)構(gòu)聲腔耦合系統(tǒng)的振型耦合特性分析及應(yīng)用[J].吉林大學(xué)學(xué)報(bào)(工學(xué)版)，2012(5)：1120-1126.

[11]張學(xué)榮.NJ6400白車身結(jié)構(gòu)動(dòng)力學(xué)分析及優(yōu)化設(shè)計(jì)[J].機(jī)械研究與應(yīng)用，2002(3)：38-40.

(責(zé)任編輯劉舸)

Numerical Analysis and Design to Lower Vehicle Interior Noise

YANG Yang, LU Sen-lin

(School of Automotive and Traffic Engineering, Jiangsu University, Zhenjiang 212013, China)

In view of the general formula for the internal sound field of a ruled structure, calculation formula for determinant parameter of sound pressure response level was derived using FEM (finite element method). Taking the structure and acoustic infinite model of a passenger car body as an example, the determinant parameter of the roof part in the low frequency body was calculated. In order to decrease the determinant parameter, an improvement program of the car body was provided. The interior noise of both the original body and the improved body was calculated using acoustic software. The sound pressure level curves for driver and passengers in first three rows were analyzed. It was confirmed that the solution to achieve interior noise reduction based on calculation of determinant parameter of sound pressure level of the roof part was feasible.

car body; determinant parameter of sound pressure response level; finite element method; sound pressure level curve

2015-05-18

江蘇省自然科學(xué)基金資助項(xiàng)目(BK2009212)

楊楊(1990—)，男，江蘇徐州人，碩士，主要從事車輛振動(dòng)噪聲控制方面的研究，E-mail:15862227665@163.com; 陸森林(1957—)，男，江蘇江陰人，博士，教授，主要從事車輛振動(dòng)噪聲控制方面的研究。

format：YANG Yang, LU Sen-lin.Numerical Analysis and Design to Lower Vehicle Interior Noise[J].Journal of Chongqing University of Technology(Natural Science)，2016(9):26-32.

10.3969/j.issn.1674-8425(z).2016.09.004

U461;TB532

A

1674-8425(2016)09-0026-07

引用格式：楊楊，陸森林.車內(nèi)噪聲的數(shù)值化分析與降噪設(shè)計(jì)[J].重慶理工大學(xué)學(xué)報(bào)(自然科學(xué))，2016(9):26-32.