李少偉 林家祥 周茜琳 胡俊 孫乙城

摘 要：汽車發(fā)動機啟停以及在加速過程中產(chǎn)生的振動激勵會通過前儀表板對汽車駕駛室輻射噪聲，針對此類情況，對汽車駕駛室聲輻射進行研究。通過聲學(xué)軟件ACTRAN對駕駛室前儀表板分別在不加裝飾、添加20 mm泡沫吸音層以及添加20 mm泡沫吸音層和3 mm隔音涂層的情況，在隨機激勵下的70～250 Hz低頻段進行仿真研究。運用有限元-無限元的方法對殼體的固定點位添加激勵，選取駕駛員右耳位作為監(jiān)測點，通過對比監(jiān)測點在3種情況下的振動響應(yīng)曲線，發(fā)現(xiàn)添加吸音、隔音材料后，振動響應(yīng)曲線振動幅值范圍減小。最后在監(jiān)測點位選取頻率進行研究，導(dǎo)出相關(guān)頻率下前儀表板在3種不同情況下的汽車駕駛室聲壓云圖。通過對選取頻率下3種方案的研究分析，為吸音、隔音材料復(fù)合情況下的降噪水平研究提供參考。

關(guān)鍵詞：振動噪聲；聲固耦合；振動特性；聲輻射

中圖分類號：TB523；U467.493 DOI：10.16375/j.cnki.cn45-1395/t.2023.01.002

0 引言

汽車NVH指汽車的噪聲（noise）、振動（vibration）和舒適性（harshness）。NVH作為汽車舒適性的重要指標，反映了汽車制造商的工藝水平，已成為衡量汽車質(zhì)量不可或缺的指標之一[1]。目前汽車制造技術(shù)不斷改進提升，使各汽車制造商在同等價位汽車的動力性能以及安全性能方面的差距越來越小，因此，用戶在購買私家車時，會著重考慮車輛的駕駛舒適性[2]。噪聲和振動作為影響汽車駕駛舒適性的主要因素被各汽車制造商不斷改進優(yōu)化[3]。

國內(nèi)外眾多學(xué)者對聲固耦合做了大量研究，已積累了豐富的研究成果和實踐經(jīng)驗。1976年，Wolf Jr等[4]針對汽車聲學(xué)結(jié)構(gòu)的耦合問題建立了二維車輛模型，并使用NASTRAN軟件在20～200 Hz范圍內(nèi)進行仿真分析，針對駕駛室的板件振動和受迫響應(yīng)問題采用模態(tài)法進行分析。20世紀80年代中期有限元方法迅速普及，美國的通用公司率先引進該方法并建立精確的模型進行仿真分析。1984年，Nefske等[5]對車內(nèi)結(jié)構(gòu)進行三維建模，運用有限元法建立了駕駛室結(jié)構(gòu)和聲學(xué)耦合模型并進行了噪聲仿真分析。郭榮等[6]介紹了傳遞路徑試驗（TPA）原理、傳遞函數(shù)以及激勵力的測量方法，并使用該試驗對某乘用車型進行試驗優(yōu)化。馬天飛等[7]對某商用車模型進行簡化后，通過模態(tài)相似原則建立駕駛室的有限元模型，并通過有限元模型和聲腔耦合建立駕駛室聲固耦合模型，對駕駛室聲固耦合模型進行試驗，驗證了模型的正確性，通過對實車道路試驗得到的振動噪聲響應(yīng)仿真結(jié)果與試驗結(jié)果進行對比分析，利用模型預(yù)測車內(nèi)振動噪聲水平。

目前在汽車NVH聲學(xué)研究中，普遍使用HyperWorks等有限元軟件進行處理分析，但是在處理聲學(xué)問題中，有限元軟件操作繁瑣，影響效率，有一定局限性。聲學(xué)軟件ACTRAN具有完備的結(jié)構(gòu)單元庫、材料庫、邊界條件以及求解方法等，在處理聲學(xué)問題上有著得天獨厚的優(yōu)勢[8]。本文基于聲學(xué)軟件ACTRAN在不同聲學(xué)材料情況下對汽車前儀表板進行仿真分析及優(yōu)化，為分析復(fù)雜的結(jié)構(gòu)振動和聲輻射以及為后續(xù)ACTRAN結(jié)合實際案例仿真優(yōu)化提供參考。

1 聲學(xué)理論基礎(chǔ)

結(jié)構(gòu)振動的聲學(xué)問題就是結(jié)構(gòu)振動與聲之間的相互關(guān)系。結(jié)構(gòu)振動會產(chǎn)生聲，聲也會引起結(jié)構(gòu)振動，在工程上是兩個不同的問題，但是在力學(xué)表達上是一個問題，僅是力學(xué)系統(tǒng)的輸入不同[8]。在有限元模型中，結(jié)構(gòu)控制方程和聲場控制方程通過運動和壓力進行耦合[9]。

結(jié)構(gòu)振動的有限元動力學(xué)方程可以表示為：

[Msu+Ksu=Fs+Fas]. （1）

其中：[u]為位移；[Ms]為質(zhì)量矩陣；[Ks]為剛度矩陣；[Fs]為結(jié)構(gòu)載荷；[Fas]為聲場對結(jié)構(gòu)的作用。

[Fas=Feas=ΓeNTspndΓe]. （2）

其中：[Γ]為邊界；[Ns]為結(jié)構(gòu)單元形函數(shù)；[p]為聲壓；[n]為方向余弦。根據(jù)聲學(xué)流體單元的插值關(guān)系，有：

[p=Nape] . （3）

其中，[Na]為聲學(xué)流體單元的形函數(shù)。

[Feas=ΓeNTspndΓe=ATepe] ， （4）

[Fas=Feas=ATepe=ATp]. （5）

把式（5）代入方程（1）得聲學(xué)有限元方程式，得：

[Msu+Ksu-ATp=Fs]. （6）

式（6）對應(yīng)的頻域形式為：

[Ks-ω2Msu-ATp=Fs]. （7）

聲學(xué)有限元方程式為：

[-Ka-iωCa+ω2MaΨ=F]. （8）

其中：i為虛數(shù)；[Ψ]為單元節(jié)點勢向量。為了簡化推導(dǎo)過程，因此不考慮流體速度、吸聲邊界條件，但存在假設(shè)聲源，可得：

[p=-iωρΦa=-iωΨ]. （9）

式（8）右端可表示為：

[F=Fa+Fsa]. （10）

其中：[Fa]是聲載荷；[Fsa]是由結(jié)構(gòu)振動產(chǎn)生的。得：

[Fesa=-iωΓeNTaundΓe]. （11）

其中：[u]為結(jié)構(gòu)振動位移；[n]為法線方向。根據(jù)結(jié)構(gòu)單元插值關(guān)系，有：

[u=Nsue]， （12）

[Fesa=-iωATu]. （13）

把式（9）、式（10）、式（13）代入式（8），忽略阻尼影響，得：

[Ka-ω2Map-ω2Au=Fa]. （14）

由式（7）和式（14）得：

[Ks-ω2Ms-AT-ω2AKa-ω2Maup=FsFa] . （15）

從而得出聲固耦合矩陣方程，方程中的矩陣[A]為耦合矩陣。通過對耦合矩陣方程進行分析可知，當固體在載荷作用下振動時會在周圍聲腔內(nèi)產(chǎn)生聲波，聲腔內(nèi)的聲源輻射聲波時也會對結(jié)構(gòu)進行作用，引起固體的振動[10]。