基于FEM—MODENA的加筋板聲—固耦合分析

于士甲　張鵬　李彥斌　吳邵慶　費(fèi)慶國(guó)　韓曉林

摘要：模態(tài)能量分析（MODENA）是近年發(fā)展的研究聲-固耦合問(wèn)題的新方法。結(jié)合有限元法（FEM），形成基于FEM-MODENA的適用于復(fù)雜結(jié)構(gòu)的聲一固耦合動(dòng)響應(yīng)分析方法。以規(guī)則的平板-長(zhǎng)方體聲腔耦合系統(tǒng)為研究對(duì)象驗(yàn)證FEM-MODENA方法的準(zhǔn)確性和有效性，然后將該方法應(yīng)用于復(fù)雜加筋板-聲腔耦合系統(tǒng)的動(dòng)響應(yīng)分析中。結(jié)果表明：FEM-MODENA對(duì)復(fù)雜聲-固耦合系統(tǒng)的振動(dòng)能量響應(yīng)預(yù)示具有較好的精度，是一種可靠的動(dòng)響應(yīng)分析手段；對(duì)于系統(tǒng)在局部分析頻段內(nèi)的振動(dòng)響應(yīng)預(yù)示，F(xiàn)EM-MODENA仍可給出較為精確的結(jié)果。

關(guān)鍵詞：聲-固耦合；加筋板；模態(tài)能量分析；有限元法；振動(dòng)能量

引言

在航空、航天、機(jī)械、船舶等眾多工程領(lǐng)域中，廣泛存在噪聲環(huán)境引起的振動(dòng)問(wèn)題。隨著動(dòng)力機(jī)械設(shè)備功率的增大，結(jié)構(gòu)在振動(dòng)和噪聲方面的問(wèn)題日益嚴(yán)重。尤其是對(duì)于典型薄壁結(jié)構(gòu)，結(jié)構(gòu)和噪聲之問(wèn)易形成耦合效應(yīng)，即噪聲導(dǎo)致結(jié)構(gòu)發(fā)生振動(dòng)，結(jié)構(gòu)振動(dòng)反過(guò)來(lái)影響聲場(chǎng)的分布，進(jìn)而引起結(jié)構(gòu)和設(shè)備元件的振動(dòng)，造成儀器設(shè)備故障。因此，在設(shè)計(jì)階段很有必要考慮聲-固耦合效應(yīng)對(duì)結(jié)構(gòu)動(dòng)響應(yīng)的影響。

聲-固耦合響應(yīng)分析方法主要有試驗(yàn)分析方法和數(shù)值分析方法。試驗(yàn)分析結(jié)果真實(shí)可信，但巨大的試驗(yàn)耗費(fèi)、較長(zhǎng)的試驗(yàn)周期以及較為有限的試驗(yàn)條件和工況限制了其進(jìn)一步的發(fā)展與應(yīng)用。近年來(lái)，隨著計(jì)算機(jī)性能的提高，數(shù)值分析方法成為一種重要的分析手段，可為試驗(yàn)提供必要的補(bǔ)充及參考。有限元一有限元法（FEM-FEM）、有限元一邊界元法（FEM-BEM）和統(tǒng)計(jì)能量法（SEA）為目前較為常用的三種聲-固耦合響應(yīng)數(shù)值分析方法。FEM-FEM和FEM-BEM廣泛應(yīng)用于中低頻聲一固耦合分析，可適應(yīng)任意形狀的結(jié)構(gòu)及激勵(lì)形式，在結(jié)構(gòu)單元尺寸滿足一定條件后，能給出相當(dāng)精確的結(jié)果。但是FEM-FEM和FEM-BEM無(wú)法用于高頻動(dòng)響應(yīng)分析，主要原因有：當(dāng)分析頻率升高時(shí)，若要保證FEM和BEM有足夠的精度，須對(duì)單元進(jìn)行細(xì)化，這將導(dǎo)致計(jì)算規(guī)模呈指數(shù)級(jí)增長(zhǎng)；在高頻段，結(jié)構(gòu)響應(yīng)對(duì)結(jié)構(gòu)參數(shù)十分敏感，即使單元尺寸再小，F(xiàn)EM結(jié)果已不具有參考性。SEA理論通過(guò)引入均一化假設(shè)，只關(guān)心不同子系統(tǒng)落在分析頻帶內(nèi)的模態(tài)群之問(wèn)統(tǒng)計(jì)意義上的功率傳輸特性，即耦合損耗因子，以平均振動(dòng)能量為基本未知量建立系統(tǒng)功率流平衡方程。SEA方法對(duì)結(jié)構(gòu)細(xì)節(jié)要求不嚴(yán)格，采用“能量”的觀點(diǎn)，關(guān)注時(shí)域、頻域和空問(wèn)上的統(tǒng)計(jì)平均值，適用于高頻段內(nèi)密集模態(tài)的復(fù)雜結(jié)構(gòu)。但實(shí)踐中SEA的各項(xiàng)假設(shè)難以得到完全滿足，且SEA結(jié)果的精度也難以評(píng)估。

Maxit等基于結(jié)構(gòu)一聲場(chǎng)耦合系統(tǒng)的雙模態(tài)方程（Dual Modal Formulation），推導(dǎo)獲得了子系統(tǒng)模態(tài)振動(dòng)能量平衡方程，提出統(tǒng)計(jì)模態(tài)能量分布分析（Statistical modal Energy distribution Analysis，SInEdA）理論，該方法為研究聲一固耦合問(wèn)題提供了新的思路。SEA理論在均一化假設(shè)下以分析頻帶內(nèi)平均的各子系統(tǒng)總振動(dòng)能量為基本未知量，而SInEdA理論摒棄了均一化假設(shè)，以分析頻帶內(nèi)平均的各子系統(tǒng)各模態(tài)的振動(dòng)能量為基本未知量，可以考慮分析頻帶內(nèi)模態(tài)稀疏的情形。Totaro等在SInEdA理論的框架下，推導(dǎo)獲得了單一頻率點(diǎn)下子系統(tǒng)模態(tài)振動(dòng)能量平衡方程，提出了模態(tài)能量分析（MODal ENergy Analysis，MODENA）理論。相比于SmEdA理論只可獲取分析頻帶內(nèi)平均的模態(tài)振動(dòng)能量響應(yīng)，MODENA理論可以獲取任意頻率點(diǎn)處的模態(tài)振動(dòng)能量響應(yīng)，在計(jì)算薄壁構(gòu)件聲傳遞損耗時(shí)，后者可內(nèi)在地考慮非共振傳輸（Non-resonant Transmission）。各子系統(tǒng)的模態(tài)信息為MODENA理論的重要輸入條件之一，用于計(jì)算各模態(tài)問(wèn)的耦合系數(shù)、各模態(tài)的阻尼耗散功率及各模態(tài)上的載荷輸入功率。規(guī)則簡(jiǎn)單聲一固耦合系統(tǒng)的模態(tài)信息可由解析方法求得。對(duì)于復(fù)雜聲一固耦合系統(tǒng)，可由FEM方法獲取各子系統(tǒng)模態(tài)信息，進(jìn)而結(jié)合MODENA理論，分析獲得各子系統(tǒng)振動(dòng)能量響應(yīng)。相比于FEM分析方法，MODENA理論以統(tǒng)一的未知量（能量）建立系統(tǒng)動(dòng)力學(xué)方程，僅需獲取各單獨(dú)子系統(tǒng)模態(tài)信息，分析效率高。并且MODENA理論可以通過(guò)能量這一變量與其他適用于中高頻的能量方法相結(jié)合，可形成適用于寬頻激勵(lì)的復(fù)雜結(jié)構(gòu)的響應(yīng)分析。

本文結(jié)合FEM方法與MODENA理論，形成了基于FEM-MODENA的適用于復(fù)雜聲一固耦合系統(tǒng)的動(dòng)響應(yīng)分析方法。首先以規(guī)則的平板一長(zhǎng)方體聲腔耦合系統(tǒng)為研究對(duì)象，驗(yàn)證FEM-MODENA方法的準(zhǔn)確性和有效性；然后以一加筋板-聲腔耦合系統(tǒng)為研究對(duì)象，驗(yàn)證FEM-MODENA方法對(duì)復(fù)雜耦合系統(tǒng)的適用性；最后研究了特定分析頻帶內(nèi)的模態(tài)對(duì)FEM-MODENA方法在該頻帶內(nèi)分析結(jié)果的貢獻(xiàn)量。