呂國(guó)策

長(zhǎng)春職業(yè)技術(shù)學(xué)院，吉林長(zhǎng)春 130000

基于邊界元法聲場(chǎng)輻射預(yù)測(cè)研究

呂國(guó)策

長(zhǎng)春職業(yè)技術(shù)學(xué)院，吉林長(zhǎng)春 130000

采用邊界元法預(yù)測(cè)結(jié)構(gòu)振動(dòng)的聲場(chǎng)聲壓，已經(jīng)得到廣泛研究，然而邊界積分方程建立的過(guò)程中遇到奇異積分，雖然此時(shí)也采用高斯積分，但計(jì)算收斂性降低。本文針對(duì)奇異積分問(wèn)題采取廣義極坐標(biāo)變換，對(duì)比了兩種結(jié)構(gòu)網(wǎng)格對(duì)邊界元法計(jì)算結(jié)果的影響，結(jié)果表明，采用廣義極坐標(biāo)變換后，邊界元法可以預(yù)測(cè)振動(dòng)輻射聲場(chǎng)聲壓，當(dāng)頻率升高時(shí)，某些頻率結(jié)果誤差增加，為此，在提高網(wǎng)格質(zhì)量后，在低頻情況兩種結(jié)果與解析解極為接近，高頻段較好的網(wǎng)格質(zhì)量度數(shù)值計(jì)算結(jié)果的精度有所提高，但計(jì)算時(shí)間乘倍增加，不利于工程實(shí)際，對(duì)于大型問(wèn)題要綜合考慮，在不影響計(jì)算精度的同時(shí)，可以適當(dāng)降低網(wǎng)格數(shù)量，以降低計(jì)算時(shí)間。

邊界元；奇異積分；廣義極坐標(biāo)變換

由于機(jī)械振動(dòng)產(chǎn)生強(qiáng)烈的噪聲，工程中振動(dòng)產(chǎn)生的噪聲越來(lái)越受到人們的關(guān)注，而噪聲對(duì)現(xiàn)場(chǎng)工作人員產(chǎn)生極大心理和勝利危害，為解決此問(wèn)題，各國(guó)都相應(yīng)出臺(tái)了相應(yīng)的措施限制現(xiàn)場(chǎng)噪聲，機(jī)械制造商相應(yīng)的也對(duì)其產(chǎn)生的噪聲采取了相應(yīng)的控制措施。對(duì)生產(chǎn)制造商而言，在進(jìn)行制造前提前預(yù)知機(jī)械產(chǎn)生的噪聲強(qiáng)度對(duì)后續(xù)工作有指導(dǎo)意義，所以需要一種有效地計(jì)算輻射聲場(chǎng)聲壓方法[1-3]。

1 邊界元方程的建立

聲在介質(zhì)中傳播，假設(shè)介質(zhì)為理想流體，流體不存在粘滯性，沒(méi)有擾動(dòng)，宏觀上靜止，且介質(zhì)絕熱，并且聲波是小振幅聲波，聲壓遠(yuǎn)小于介質(zhì)靜態(tài)聲壓，由以上假設(shè)可以將非線性問(wèn)題轉(zhuǎn)化為線性聲輻射問(wèn)題，而我們關(guān)注的多數(shù)問(wèn)題都是脈動(dòng)振源產(chǎn)生的穩(wěn)定的簡(jiǎn)諧聲波，所以有聲壓輻射公式[4-5]

根據(jù)理想流體小振幅聲波輻射的波動(dòng)方程，可以得到helmholtz等式。

其中k為波數(shù)， /k cω= ，其伴隨方程為

伴隨方程的基本解為

我們僅考慮輻射聲壓，所以有B=0，有伴隨方程的基本解

利用第二格林公式建立邊界積分方程

假設(shè)基本解和聲壓滿足此方程，得到場(chǎng)點(diǎn)聲壓與邊界信息之間的關(guān)系

2 邊界奇異積分的處理

對(duì)于給定的結(jié)構(gòu)，可能僅了解邊界表面振動(dòng)部分信息，然而通過(guò)以上邊界積分方程可知，要得到場(chǎng)點(diǎn)聲壓，必須要有所有邊界信息，包括表面振動(dòng)速度和聲壓，當(dāng)給定邊界條件時(shí)，可以將場(chǎng)點(diǎn)p配置到邊界表面上，然而基本解這時(shí)存在奇異性，對(duì)于奇異積分的處理，本文采用極坐標(biāo)變換法引進(jìn)雅克比消除奇異積分的方法。

如圖1所示，將三角主單元轉(zhuǎn)化為正方形單位單元，由此的到轉(zhuǎn)化關(guān)系[6]

由上式帶入圖坐標(biāo)點(diǎn)可得

變換后的得到的雅克比

圖1 廣義極坐標(biāo)變換

3 脈動(dòng)求輻射聲場(chǎng)驗(yàn)證

脈動(dòng)球源是指隨時(shí)間在進(jìn)行均勻振動(dòng)的球體，表面做簡(jiǎn)諧振動(dòng)，脈動(dòng)求向它存在的聲場(chǎng)輻射聲壓，滿足本文假設(shè)條件，輻射聲場(chǎng)聲壓可以通過(guò)解析公式計(jì)算解析解，作為驗(yàn)證邊界元法的一種有效手段，與邊界元法計(jì)算得到的數(shù)值解進(jìn)行對(duì)比，驗(yàn)證邊界元法的有效性。

對(duì)脈動(dòng)球半徑r0=0.5(m)，介質(zhì)密度ρ0=1.21kg/m3，脈動(dòng)求表面振動(dòng)速度ua=1m/s建兩種單元模型，脈動(dòng)球聲場(chǎng)輻射聲壓解析解計(jì)算公式有

對(duì)脈動(dòng)球邊界離散采用四邊形網(wǎng)格，選用兩種網(wǎng)格尺寸，得到場(chǎng)點(diǎn)網(wǎng)格，第一種為四邊形網(wǎng)格單元，單元尺寸為50mm，節(jié)點(diǎn)數(shù)為1269，單元數(shù)為1267，第二種為單元尺寸30mm，節(jié)點(diǎn)數(shù)和網(wǎng)格數(shù)分別為3164、3162，對(duì)比網(wǎng)格質(zhì)量提高后對(duì)計(jì)算聲場(chǎng)聲壓精度的影響。

由于Fortran語(yǔ)言計(jì)算速度較快，適合大型程序開(kāi)發(fā)與維護(hù)，針對(duì)大型計(jì)算問(wèn)題，可以將模型分配幾種情況單獨(dú)計(jì)算，這樣可以縮減計(jì)算時(shí)間，本文邊界元積分方程程序采Fortran語(yǔ)言編寫(xiě)。

對(duì)于脈動(dòng)球場(chǎng)聲壓本文分別計(jì)算不同場(chǎng)點(diǎn)和不同頻率下聲壓進(jìn)行對(duì)比，表1為場(chǎng)點(diǎn)網(wǎng)格1對(duì)應(yīng)的不同頻率下場(chǎng)點(diǎn)聲壓解析解和數(shù)值解對(duì)比結(jié)果，從表可以的得出，在低頻段，解析解與數(shù)值解吻合較好，當(dāng)頻率升高后，數(shù)值解相應(yīng)誤差增加。

表2為場(chǎng)點(diǎn)網(wǎng)格2對(duì)應(yīng)解析解與數(shù)值解對(duì)別結(jié)果，與表3變化趨勢(shì)類似，隨頻率增加后聲壓幅值誤差有所增加，對(duì)比表1，在頻率700Hz以上，計(jì)算效果明顯改善，分析原因由于脈動(dòng)球建模劃分網(wǎng)格質(zhì)量原因造成表4結(jié)果更近解析解。

以上結(jié)果分別對(duì)比了不同網(wǎng)格不同頻率對(duì)比結(jié)果，表3為網(wǎng)格1計(jì)算的到不同場(chǎng)點(diǎn)聲壓對(duì)比結(jié)果，頻率為50Hz，不同場(chǎng)點(diǎn)得到的聲壓在低頻情況邊界元法計(jì)算得到的數(shù)值結(jié)果與解析結(jié)果誤差較小，并且場(chǎng)點(diǎn)與脈動(dòng)球源距離變化后誤差并沒(méi)有明顯變

表4為在1500Hz頻率下計(jì)算場(chǎng)點(diǎn)聲壓結(jié)果對(duì)比，當(dāng)頻率升高后，場(chǎng)點(diǎn)聲壓解析解與數(shù)值解仍相差較小，誤差在1%以下，滿足工程應(yīng)用需要。

對(duì)于兩種網(wǎng)格，計(jì)算時(shí)間也有所區(qū)別，網(wǎng)格1節(jié)點(diǎn)較少，單元較少，而網(wǎng)格2節(jié)點(diǎn)數(shù)和單元數(shù)為網(wǎng)格1的2.5倍，兩種網(wǎng)格計(jì)算時(shí)間分別為277.479 （s）、1568.574（s），網(wǎng)格2的計(jì)算時(shí)間是1的5.65倍，嚴(yán)重增加了計(jì)算時(shí)間。

4 結(jié)論

場(chǎng)點(diǎn)聲壓隨頻率變化數(shù)值解與解析解吻合較好，當(dāng)頻率增加聲壓略微產(chǎn)生波動(dòng)。場(chǎng)點(diǎn)與脈動(dòng)球距離變化不大時(shí)，聲壓誤差沒(méi)有明顯變化。

網(wǎng)格質(zhì)量提高后對(duì)聲場(chǎng)聲壓預(yù)測(cè)低頻精度沒(méi)有明顯提高，高頻誤差減小，但計(jì)算時(shí)間數(shù)倍增長(zhǎng)。

存在誤差可能原因是在離散邊界是認(rèn)為邊界單元為平面單元，造成誤差。在計(jì)算過(guò)程中考慮到實(shí)際問(wèn)題不可能將邊界離散很多單元，所以本文離散邊界單元較大，也造成誤差。但這不影響對(duì)實(shí)際問(wèn)題的分析。

表1 網(wǎng)格1不同頻率對(duì)比結(jié)果

表2 網(wǎng)格2不同頻率對(duì)比結(jié)果

表3 網(wǎng)格1場(chǎng)點(diǎn)聲壓結(jié)果

表4 網(wǎng)格2場(chǎng)點(diǎn)聲壓結(jié)果

[1]高貴兵，張道兵，岳文輝，等.一種具有唯一解的邊界元法計(jì)算三維聲散射問(wèn)題[J].機(jī)械科學(xué)與技術(shù)，2013，32（8）：1172-1175.

[2]郭偉，蔡明，馬俊.聲輻射邊界元數(shù)值計(jì)算[J].船海工程2011，40（1）：101-103.

[3]黎勝，趙德有.用邊界元法計(jì)算結(jié)構(gòu)振動(dòng)輻射聲場(chǎng)[J]，大連理工大學(xué)學(xué)報(bào)，2000，40（4）：391-394.

[4]杜功煥，朱哲民，龔秀芬.聲學(xué)基礎(chǔ)[M].南京：南京大學(xué)出版社.

[5]Ursell F. On the exterior problems of acoustics[C]//Mathematical Proceedings of the Cambridge Philosophical Society. Cambridge University Press， 1973， 74(01): 117-125.

[6] Boundary element acoustics: Fundamentals and computer codes[M]. Wit Pr/Computational Mechanics，2000.

O29

A

1674-6708（2015）145-0113-02