徐明林，張利華，皮君濤，葛澤昊，劉聰睿，周 波

（1.大連理工大學(xué) 工程裝備結(jié)構(gòu)分析國家重點(diǎn)實驗室 運(yùn)載工程與力學(xué)學(xué)部船舶工程學(xué)院，遼寧 大連 116024；2.煙臺中集來福士海洋工程有限公司，山東 煙臺 64000）

隨著國民生活水平的提升，人們對環(huán)境舒適性的要求也越來越高，其中振動與噪聲是最常被考慮到的重要指標(biāo)。板類結(jié)構(gòu)應(yīng)用廣泛，是構(gòu)成汽車、高速列車、船舶、飛行器等交通工具的艙體結(jié)構(gòu)以及海洋平臺等建筑物艙壁結(jié)構(gòu)的重要組成部分[1]。從聲源到聲音接收端，噪聲的傳播路徑主要有結(jié)構(gòu)聲傳播和空氣聲傳播兩種。板殼結(jié)構(gòu)作為傳播路徑上的重要環(huán)節(jié)，對艙室內(nèi)的噪聲環(huán)境產(chǎn)生直接的影響。因此，對開孔板結(jié)構(gòu)的振動和傳聲特性進(jìn)行研究，并在此基礎(chǔ)上抑制結(jié)構(gòu)的振動和噪聲，可為工程結(jié)構(gòu)設(shè)計與技術(shù)應(yīng)用提供參考，具有重要意義。

工程中為了控制結(jié)構(gòu)物重量與實現(xiàn)相應(yīng)應(yīng)用功能比如通風(fēng)換氣、通過線路等，常常會在板上開孔，比如海洋平臺艙壁結(jié)構(gòu)、船舶內(nèi)底板架腹板等。開孔板結(jié)構(gòu)作為工程中的熱點(diǎn)研究對象，國內(nèi)外不僅進(jìn)行了應(yīng)力集中、變形等方面的研究，還對此類結(jié)構(gòu)的聲振特性開展了深入的研究。Kwak等[2]運(yùn)用獨(dú)立坐標(biāo)耦合法研究了矩形開孔板結(jié)構(gòu)振動特性。Putra等[3]以孔徑為毫米級的多孔板為研究對象，建立了多孔板聲輻射模型，進(jìn)行了多孔板的輻射實驗并將數(shù)值結(jié)果與實驗結(jié)果進(jìn)行了比較，推導(dǎo)出了射孔效應(yīng)的近似公式。Chang 等[4]運(yùn)用有限元法研究了內(nèi)部開孔薄板結(jié)構(gòu)的振動特性。陳美霞等[5]等利用SYSNOISE 軟件采用有限元法+邊界元法計算了四邊簡支平板在空氣中和水中的振動與聲輻射特性，并將結(jié)果與有關(guān)解析解進(jìn)行了比較分析，驗證了方法的合理性和準(zhǔn)確性。王剛[6]從能量的角度出發(fā)，建立了“聲腔-開孔板-聲腔”耦合模型，對三角形開孔板和復(fù)雜開孔板的聲場進(jìn)行了預(yù)測，并分析了開孔板的聲場特性。邱昌林等[7]采用有限元與邊界元相結(jié)合的方法，對開有圓孔、四邊簡支、無障板鋼制平板的水下振動及聲輻射特性進(jìn)行了研究。

綜上，現(xiàn)有文獻(xiàn)對開孔板的聲輻射研究主要集中在對板結(jié)構(gòu)的振動特性以及傳播介質(zhì)的研究上，并沒有深入探索影響板結(jié)構(gòu)聲輻射的主要參數(shù)以及聲輻射的指向性特性。因此，本文基于間接邊界元法，應(yīng)用商用軟件Patran 和Virtual.Lab Acoustics 對四邊剛固的開孔板進(jìn)行了振動和聲學(xué)特性的仿真計算，并對開孔板結(jié)構(gòu)的輻射聲功率以及指向性進(jìn)行分析，系統(tǒng)地探討了開孔形狀、開孔大小以及加筋板對結(jié)構(gòu)聲輻射特性的影響。

1 基本理論

間接邊界元法是由直接邊界元法推導(dǎo)得來，即在邊界的兩側(cè)同時建立直接邊界元的Helmholtz 積分方程，然后將兩積分方程相減，計算得到任意觀測點(diǎn)的聲壓[8]，如式（1）所示：

式中：σ(rq)為表面兩側(cè)的法向壓力梯度差（單層勢），μ(rq) 為表面兩側(cè)的聲壓差（雙層勢），當(dāng)一側(cè)的聲壓p(rq2)=0時，μ(rq)=p(rq1)。將場點(diǎn)定義在邊界上，其位置向量為rq，則可以建立邊界條件與未知主要變量的關(guān)系：

式中：υ(rp)為邊界模型表面上rp處的振速。由式（2）可知，只要以σ(rp)、μ(rq)中的一個為未知量即可。然后，對邊界進(jìn)行離散，則邊界元上的主要未知變量可以用模型上節(jié)點(diǎn)的未知變量及形函數(shù)表示，得到一般形式的方程為：

式中：[D]為對稱矩陣，由邊界幾何形狀和介質(zhì)決定，通過格林函數(shù)及格林函數(shù)的偏導(dǎo)積分求得；{x}為邊界單元模型上的未知向量，即σ(rp)、μ(rq)中的一個或兩者的組合，且需滿足未知量的總個數(shù)等于邊界上的節(jié)點(diǎn)數(shù)；{Fa}為作用在流體上的載荷向量。

2 數(shù)值仿真

2.1 模型的建立

考慮到本文將探討加筋板對板結(jié)構(gòu)的聲輻射的影響，這里選用了參考文獻(xiàn)[9]中船舶艙室的典型加筋板尺寸，材質(zhì)屬性選為普通鋼材，具體的模型參數(shù)如表1所示。

表1 模型參數(shù)

在Virtual.Lab Acoustics 中，建立聲學(xué)仿真計算模型一般包括確定結(jié)構(gòu)網(wǎng)格、聲學(xué)網(wǎng)格、接受數(shù)據(jù)的場點(diǎn)網(wǎng)格以及邊界條件等。本文用到的結(jié)構(gòu)網(wǎng)格和聲學(xué)網(wǎng)格均通過有限元軟件Patran生成，如圖1至圖2所示，分別為開孔平板與帶有加筋板的開孔板的有限元網(wǎng)格。開孔板結(jié)構(gòu)的面板和加筋板均采用2維殼單元，面板與加筋板的連接通過合并節(jié)點(diǎn)（Equivalence）實現(xiàn)。在賦予板結(jié)構(gòu)材料屬性后，在板的四周施加剛固約束，求解結(jié)構(gòu)模態(tài)，生成了“bdf”和“op2”兩種格式的文件。

圖1 開孔平板結(jié)構(gòu)網(wǎng)格

圖2 帶有加筋板的開孔板結(jié)構(gòu)網(wǎng)格

本文研究的對象是板結(jié)構(gòu)，對于此類結(jié)構(gòu)采用間接邊界元法進(jìn)行聲振耦合計算時，要求聲學(xué)邊界網(wǎng)格與結(jié)構(gòu)網(wǎng)格相重合。故以“bdf”格式文件作為聲學(xué)邊界網(wǎng)格，以含有模態(tài)求解結(jié)果的“op2”格式文件作為結(jié)構(gòu)網(wǎng)格。對于聲學(xué)邊界網(wǎng)格與結(jié)構(gòu)網(wǎng)格中的重名節(jié)點(diǎn)和單元，可以通過Virtual.Lab Acoustics中的“Check Id Conflicts”進(jìn)行重新編號。對于聲學(xué)有限元和邊界元模型，一般要求在一個波長內(nèi)至少有6 個單元，即單元長度應(yīng)滿足：L≤c/6fmax[10]，c為聲音在相應(yīng)介質(zhì)中的傳播速度，fmax為最大求解頻率，本文研究的最大頻率為1 000 Hz，求得最大單元長度為56.7 mm，為了保證開孔板結(jié)構(gòu)模態(tài)具有一定的精度，這里取25 mm。為了得到板結(jié)構(gòu)輻射聲功率以及指向性聲功率，在距板邊緣50 mm 處施加大小為10 N、方向垂直平板的力，并通過符合ISO(International Organization for Standardization)聲功率測試標(biāo)準(zhǔn)的場點(diǎn)網(wǎng)格與半徑為3000mm的圓形指向性場點(diǎn)網(wǎng)格進(jìn)行聲功率的采集。另外，為了消除聲音沿夾層板邊緣的衍射，可以運(yùn)用Virtual.Lab Acoustics 軟件提供的障板功能。具體仿真模型如圖3所示。

圖3 聲學(xué)仿真模型

2.2 場點(diǎn)輻射聲功率的分布

模型的準(zhǔn)確性直接影響數(shù)值仿真結(jié)果的準(zhǔn)確性與可靠性[11]，故建立與文獻(xiàn)[12]相同的模型，運(yùn)用Patran計算模型的前3階模態(tài)，并對計算結(jié)果進(jìn)行了比較分析。如圖4所示，圓形開孔板結(jié)構(gòu)的前3階模態(tài)振型相似；如表2所示，本文與文獻(xiàn)中計算的固有頻率相差較小，最大誤差在5%以內(nèi)，驗證了模型的準(zhǔn)確性。

圖4 開圓孔板的模態(tài)振型與文獻(xiàn)[12]比較

表2 開孔板的振動頻率比較/Hz

為了保證仿真計算結(jié)果的準(zhǔn)確性，本文選取了文獻(xiàn)[13]中的四邊固支板算例進(jìn)行驗證。矩形板的物理參數(shù)為：a=1m，r=a/b=1.2，h=6 mm，ρ=7 800 kg/m3，E=2.07×1011Pa，μ=0.3，ρ0=1.2 kg/m3，c=340 m/s。如圖5所示，四邊固支板在中心點(diǎn)受單位力的輻射聲功率，數(shù)值仿真結(jié)果和理論計算基本吻合，差值均在3 dB以內(nèi)。

圖5 四邊固支矩形板輻射聲功率與文獻(xiàn)[13]對比

基于上述模型，對開孔板結(jié)構(gòu)的聲輻射進(jìn)行數(shù)值計算。圓形開孔板受到垂直方向大小為10 N力，頻率為315 Hz，輻射聲功率分布情況如圖6所示。在激勵側(cè)與激勵側(cè)正對面，開孔板輻射聲功率最大；聲功率由激勵處和正對激勵側(cè)向左右兩側(cè)衰減，并在側(cè)面的斜對角上方出現(xiàn)極小值。根據(jù)場點(diǎn)網(wǎng)格數(shù)據(jù)，得到開孔板隨頻率變化的輻射聲功率曲線。如圖7所示，開孔板的輻射聲功率隨著頻率增加而增加。且開孔板的輻射聲功率在31.5 Hz、100 Hz、400 Hz 附近出現(xiàn)峰值，這與結(jié)構(gòu)固有頻率31.15 Hz（2階），104.89 Hz（7階），396.94 Hz（28階）相對應(yīng)，可認(rèn)為是由結(jié)構(gòu)共振引起。

圖6 開圓孔板輻射聲功率分布圖

圖7 圓形開孔板輻射聲功率曲線

圓形開孔板在各個頻率下的指向性輻射聲功率分布如圖8所示。開孔板結(jié)構(gòu)的輻射聲功率整體上沿著孔的形狀逐漸向外擴(kuò)散，在部分方向出現(xiàn)峰谷值。其中，63 Hz 頻率下輻射聲功率擴(kuò)散均勻，各個方向的輻射聲功率相差較?。?1.5 Hz頻率下輻射聲功率在板±90度方向出現(xiàn)低谷；在500 Hz和1 000 Hz下輻射聲功率在±135度方向出現(xiàn)低谷。

圖8 圓形開孔板指向性輻射聲功率圖

3 結(jié)構(gòu)參數(shù)影響研究

3.1 形狀對聲輻射的影響

為了提高板材的利用效率，人們針對不同的情況常在板上開不同形狀的孔。本文參考文獻(xiàn)[14]，選取了7 種開孔形狀，如圖9所示，通過改變孔的形狀，探討不同開孔形狀的板材的輻射聲功率。

圖9 不同開孔形狀

如圖10所示，當(dāng)頻率在200 Hz 以下時，板結(jié)構(gòu)的輻射聲功率出現(xiàn)振蕩現(xiàn)象；當(dāng)頻率在200 Hz～400 Hz范圍內(nèi)，相較于未開孔板，開孔板的輻射聲功率均大于未開孔板的，且出現(xiàn)了峰值；在400 Hz～800 Hz 范圍內(nèi)，開孔板的輻射聲功率基本均小于未開孔板，其中矩形、六邊形和圓形的開孔影響較為顯著。

圖10 不同開孔形狀的開孔板輻射聲功率曲線圖

不同開孔形狀下的指向性輻射聲功率如圖11所示。在31.5 Hz頻率下，未開孔板和開孔板的輻射聲功率均呈兩個半圓形，在±90度方向出現(xiàn)低谷；在63 Hz 頻率下，板的輻射聲功率場呈橢圓形或圓形，均勻向外擴(kuò)散；在頻率為125 Hz 下，板的輻射聲功率場變成“十”字形，在4個正方向出現(xiàn)峰值，在斜對角方向出現(xiàn)谷值。當(dāng)頻率為250 Hz時，矩形孔板的輻射聲功率場與未開孔板相似，但是聲功率數(shù)值較大，橢圓形孔板的輻射聲功率場呈長條形，其他開孔形狀板呈“十”字形。當(dāng)頻率為500 Hz時，板的輻射聲功率場近似為橢圓形。在頻率為1 000 Hz 時，輻射聲功率場變成近似的圓形，即向各個方向的輻射聲功率近似相等。綜上可以認(rèn)為在125 Hz及以下，板的開孔形狀主要影響輻射聲功率的大小，在250 Hz～500 Hz，開孔形狀既影響聲輻射的大小也影響擴(kuò)散方向；在頻率大于500 Hz 時，開孔形狀只在局部位置影響輻射聲功率。

圖11 不同形狀開孔板在不同頻率下的指向性輻射聲功率

3.2 開孔尺寸對聲輻射的影響

矩形是對輻射聲功率影響最大的形狀，改變矩形開孔的尺寸大小，研究不同開孔大小對板材輻射聲功率的影響。矩形開孔大小為0.2a、0.3a、0.4a時的輻射聲功率曲線如圖12所示，開孔板的輻射聲功率整體變化趨勢一致，只有峰谷位置隨著開孔大小變化而不同。開孔的大小會影響板材的固有頻率，可認(rèn)為開孔的大小主要影響板結(jié)構(gòu)的固有頻率，進(jìn)而影響聲輻射。

圖12 不同開孔大小的開孔板輻射聲功率曲線

3.3 加筋板對聲輻射的影響

加筋板是工程中常用來起輔助加強(qiáng)作用的結(jié)構(gòu)，探討采用不同加筋方式的加筋板對開孔板結(jié)構(gòu)輻射聲功率的影響。采用不同加筋形式的加筋板輻射聲功率曲線如圖13所示，在200 Hz 以下，輻射聲功率隨著頻率增加而快速增大；在200 Hz～400 Hz，加筋板降低了開孔板結(jié)構(gòu)的輻射聲功率；在400 Hz～800 Hz，加筋板的輻射聲功率大于平板，且單向加筋和雙向加筋的影響趨勢一致。單向加筋和雙向加筋板的指向性輻射聲功率如圖14至圖15所示，加筋減小了不同方向上輻射聲功率分布的差異性，使得聲輻射在各個方向更均勻。

圖13 不同加筋形式的加筋板輻射聲功率曲線

圖14 單向加筋板指向性輻射聲功率圖

圖15 雙向加筋板指向性輻射聲功率圖

4 結(jié)語

板類結(jié)構(gòu)是結(jié)構(gòu)噪聲傳播路徑上的重要環(huán)節(jié)，研究板結(jié)構(gòu)的振動和聲輻射特性具有重要意義。本文運(yùn)用聲學(xué)仿真軟件Virtual.Lab Acoustics，基于間接邊界元法對開孔板的輻射聲功率以及聲輻射的指向性進(jìn)行了研究。并在此基礎(chǔ)上，探討了不同開孔形狀、開孔大小以及加筋板對板結(jié)構(gòu)聲輻射的影響，并得出以下結(jié)論：

（1）開孔能夠減小板結(jié)構(gòu)的聲輻射，且板材的固有頻率是影響聲輻射的重要因素。

（2）在不同的頻率范圍內(nèi)，開孔形狀對聲輻射的影響不一樣。在低頻段，開孔形狀主要影響板輻射聲功率的幅值大小；在中頻段，開孔形狀會影響聲輻射的大小和方向；在高頻段，板的輻射場趨于固定形狀，開孔形狀主要在局部影響聲輻射大小。

（3）板結(jié)構(gòu)的開孔大小主要通過改變結(jié)構(gòu)固有頻率來影響聲輻射特性。

（4）單項加筋與雙向加筋對板材的聲輻射影響一致，均使之向高頻移動。在指向性上，加筋板減小了輻射聲功率的差異性，使得聲輻射在各個方向更均勻。