基于FEM/BEM的浸水結(jié)構(gòu)振動(dòng)及聲輻射特性國(guó)內(nèi)研究綜述

王崢，洪明，b，劉城（大連理工大學(xué).船舶工程學(xué)院；b.工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連116024）

基于FEM/BEM的浸水結(jié)構(gòu)振動(dòng)及聲輻射特性國(guó)內(nèi)研究綜述

王崢a，洪明a，b，劉城a
（大連理工大學(xué)a.船舶工程學(xué)院；b.工業(yè)裝備結(jié)構(gòu)分析國(guó)家重點(diǎn)實(shí)驗(yàn)室，遼寧大連116024）

文章基于有限元（FEM）或有限元結(jié)合邊界元（FEM/BEM）分析方法，對(duì)浸水結(jié)構(gòu)振動(dòng)及聲輻射特性國(guó)內(nèi)公開發(fā)表的文獻(xiàn)進(jìn)行梳理，從低頻流固耦合的振動(dòng)特性、弱耦合與強(qiáng)耦合的聲輻射特性、直/間接邊界元的非唯一性與奇異性等著重流體數(shù)值模擬方面，對(duì)國(guó)內(nèi)相關(guān)研究進(jìn)行了綜述分析，就浸水結(jié)構(gòu)振動(dòng)及聲輻射特性的研究發(fā)展動(dòng)態(tài)做出了總結(jié)，為相關(guān)科研人員研究浸水結(jié)構(gòu)振動(dòng)特性及聲輻射提供了有益參考。

有限元；直接邊界元；間接邊界元；流固耦合；聲輻射；非唯一性；奇異性

1 引言

若結(jié)構(gòu)在流體中振動(dòng)，則流體與結(jié)構(gòu)之間存在相互作用，相應(yīng)稱為流固耦合問題。根據(jù)耦合的強(qiáng)弱可以采用不同的分析模型。結(jié)構(gòu)與低密度流體介質(zhì)是弱耦合，對(duì)于介質(zhì)中結(jié)構(gòu)的振動(dòng)特性及聲輻射問題可分別求解，先計(jì)算出結(jié)構(gòu)的振動(dòng)響應(yīng)，輻射聲場(chǎng)可以通過(guò)結(jié)構(gòu)的振動(dòng)響應(yīng)來(lái)確定。對(duì)于研究浸水結(jié)構(gòu)的振動(dòng)特性及聲輻射問題，一般在物理上視為結(jié)構(gòu)與重流體介質(zhì)耦合，對(duì)此強(qiáng)耦合系統(tǒng)振動(dòng)響應(yīng)和水中聲場(chǎng)必須耦合求解。以振動(dòng)板殼聲振耦合研究為例，耦合的強(qiáng)弱可用特征量λ來(lái)判別[1]：

其中：c是聲在介質(zhì)中傳播速度；ρ是結(jié)構(gòu)質(zhì)量密度；h是板厚度；ω是結(jié)構(gòu)與流體耦合系統(tǒng)振動(dòng)頻率。當(dāng)λ＜1時(shí)，認(rèn)為系統(tǒng)是弱耦合；當(dāng)λ＞1時(shí)，認(rèn)為系統(tǒng)是強(qiáng)耦合。

當(dāng)浸水結(jié)構(gòu)在低頻振動(dòng)時(shí)，流體對(duì)結(jié)構(gòu)的作用可看成是附加質(zhì)量效應(yīng)；而高頻時(shí)流體的阻抗相當(dāng)復(fù)雜，包括質(zhì)量與阻尼兩方面的影響。一般低頻時(shí)可用有限元（FEM）或有限元和邊界元（FEM/BEM）相結(jié)合來(lái)研究結(jié)構(gòu)在流固耦合作用下的振動(dòng)及聲輻射特性，而在高頻時(shí)則用統(tǒng)計(jì)能量分析方法（SEA）來(lái)研究。

考慮到篇幅綜述文獻(xiàn)限于結(jié)構(gòu)滿足線彈性基本假定；對(duì)流體假定為理想流體且無(wú)旋，在線性小擾動(dòng)的情況下的聲輻射，而且流體區(qū)域討論的范圍限定于無(wú)限或半無(wú)限三維空間。若流體假設(shè)為可壓縮的，則聲壓或速度勢(shì)滿足Helmholtz方程，若將流體假設(shè)為不可壓縮的，則聲壓或速度勢(shì)滿足Laplace方程。當(dāng)研究結(jié)構(gòu)聲輻射問題時(shí)，考慮到聲音在介質(zhì)中傳播的速度不能為無(wú)限大，所以假設(shè)流體是可壓縮的；僅對(duì)于流固耦合作用下結(jié)構(gòu)低頻振動(dòng)的計(jì)算，則可假設(shè)流體是不可壓縮的，分析模型簡(jiǎn)單且精度可以接受。

無(wú)論是求解流固耦合振動(dòng)特性問題還是聲輻射問題，由于結(jié)構(gòu)和流體的相互作用，所以要考慮到結(jié)構(gòu)運(yùn)動(dòng)方程與流體運(yùn)動(dòng)方程的耦合方式。如果將結(jié)構(gòu)利用有限元理論進(jìn)行處理，流體利用邊界元理論進(jìn)行處理，則其耦合方式有兩種：第一種是將流體的邊界積分方程耦合到結(jié)構(gòu)的運(yùn)動(dòng)方程中，進(jìn)而求解有限元方法離散得到的運(yùn)動(dòng)方程；第二種是將結(jié)構(gòu)的運(yùn)動(dòng)方程耦合到流體的Helmholtz微分方程或Laplace微分方程中后，將Helmholtz微分方程或Laplace微分方程轉(zhuǎn)化為積分方程來(lái)進(jìn)行求解未知變量。顯然，第一種方法是利用邊界元的基本原理對(duì)流體進(jìn)行降一維處理，而第二種方法由于方程耦合過(guò)程中導(dǎo)致了Helmholtz微分方程或Laplace微分方程中的主要變量發(fā)生變化，其后續(xù)處理的積分形式也發(fā)生變化，致使處理十分困難，所以一般利用結(jié)構(gòu)有限元結(jié)合流體邊界元理論對(duì)流固耦合問題進(jìn)行處理時(shí)選用第一種方式。如果結(jié)構(gòu)和流體都利用有限元理論來(lái)處理，考慮到結(jié)構(gòu)和流體的相互作用，其結(jié)構(gòu)運(yùn)動(dòng)方程與流體運(yùn)動(dòng)方程同樣存在類似的耦合方式，只是無(wú)論哪種耦合方式對(duì)于結(jié)構(gòu)還是流體涉及到的積分并不能進(jìn)行降維處理。

2 結(jié)構(gòu)有限元與流體邊界元理論梗概

2.1 有限元理論

對(duì)于彈性體結(jié)構(gòu)采用有限元法計(jì)算其振動(dòng)特性是非常方便的，相應(yīng)的振動(dòng)運(yùn)動(dòng)方程如下：

其中：M為系統(tǒng)的質(zhì)量陣；C為系統(tǒng)的阻尼矩陣；K為系統(tǒng)的剛度陣；Fs為作用在結(jié)構(gòu)上激勵(lì)力矢量。

若計(jì)算結(jié)構(gòu)的固有振動(dòng)特性，忽略阻尼，則方程變?yōu)椋?/p>

2.2 直/間接邊界元（D-BEM/I-BEM）理論

假設(shè)流場(chǎng)為無(wú)旋、可壓縮且做小擾動(dòng)運(yùn)動(dòng)，線性齊次聲學(xué)波動(dòng)方程滿足Helmholtz微分方程：

其中：p為流體動(dòng)壓力；k=ω/c為波數(shù)；ω為圓頻率；c為流體介質(zhì)中的聲速。

假設(shè)流場(chǎng)為無(wú)旋、不可壓縮且做小擾動(dòng)運(yùn)動(dòng)，則c→∞，Helmholtz微分方程簡(jiǎn)化為L(zhǎng)aplace微分方程：

對(duì)振動(dòng)結(jié)構(gòu)在無(wú)限流體域中聲輻射的數(shù)值模擬，通常是將流體域中控制微分方程（Helmholtz微分方程）轉(zhuǎn)化成邊界上的積分方程，然后采用直接或間接邊界元法進(jìn)行數(shù)值計(jì)算。直接邊界元方法理論的Helmholtz積分方程是基于封閉結(jié)構(gòu)建立，且內(nèi)場(chǎng)外場(chǎng)問題分別是不同的分析模型，如圖1所示。

其Helmholtz積分方程如下：內(nèi)場(chǎng)問題：

圖1 直接邊界元法主要變量示意圖Fig.1 Primary variables of D-BEM

其中：p為場(chǎng)點(diǎn)聲壓；rX為場(chǎng)點(diǎn)的位置矢量；G為Helmholtz方程的Green函數(shù)；GL為L(zhǎng)aplace方程的Green函數(shù)；C1（rX）和C2（rX）分別為內(nèi)外場(chǎng)問題的邊界積分系數(shù)；S為結(jié)構(gòu)表面，無(wú)論是內(nèi)場(chǎng)問題還是外場(chǎng)問題，當(dāng)P∈S時(shí)C1（rX）或C2（rX）的表達(dá)式適用于任意形狀的結(jié)構(gòu)表面；當(dāng)S為光滑表面時(shí)，C1（rX）或C2（rX）的值均為1/2。

直接邊界元理論的Helmholtz積分方程的求解是將邊界離散后，利用配點(diǎn)法來(lái)計(jì)算結(jié)構(gòu)表面未知變量，然后利用內(nèi)外場(chǎng)問題的公式（6）和（7）等來(lái)計(jì)算場(chǎng)點(diǎn)聲壓、輻射聲功率及聲輻射效率等物理量。

對(duì)于結(jié)構(gòu)向三維半空間內(nèi)進(jìn)行聲輻射，并且結(jié)構(gòu)可以簡(jiǎn)化為鑲嵌在無(wú)限大剛性障板中的問題，直接邊界元法理論的Helmholtz積分方程可簡(jiǎn)化為Rayleigh積分形式：

其中：p為半空間介質(zhì)中或結(jié)構(gòu)表面上的任意點(diǎn)聲壓；S為結(jié)構(gòu)表面。Rayleigh積分可看成是直接邊界元法理論的Helmholtz積分方程的特殊形式。

圖2 場(chǎng)點(diǎn)在結(jié)構(gòu)封閉域外Fig.2 Field point P in the exterior acoustic domain

圖3 場(chǎng)點(diǎn)在結(jié)構(gòu)邊界上Fig.3 Field point P on the boundary

圖4 場(chǎng)點(diǎn)在結(jié)構(gòu)封閉域內(nèi)Fig.4 Field point P in the interior acoustic domain

間接邊界元法理論是利用直接邊界元理論內(nèi)場(chǎng)和外場(chǎng)的邊界積分公式推導(dǎo)得到的。如圖2-4所示，結(jié)合直接邊界元法的理論計(jì)算公式可得出結(jié)論，無(wú)論場(chǎng)點(diǎn)P是在聲場(chǎng)域中任意位置都有：

由公式（6）和（7）相加結(jié)合公式（9）可得聲場(chǎng)域中任意點(diǎn)的邊界積分公式：

由于間接邊界元理論推導(dǎo)所涉及的主要變量為振動(dòng)結(jié)構(gòu)邊界表面的聲壓差和聲壓梯度差，由公式（10）可看出，求出結(jié)構(gòu)邊界表面的聲壓差和聲壓梯度差便能得到輻射場(chǎng)，所以必須將場(chǎng)點(diǎn)分別移動(dòng)到邊界的兩個(gè)表面上，進(jìn)而推導(dǎo)出結(jié)合了邊界條件的邊界積分公式，才能進(jìn)一步求得結(jié)構(gòu)邊界表面的聲壓差和聲壓梯度差。比如，假設(shè)聲場(chǎng)滿足Neumann邊界條件，則推導(dǎo)出的結(jié)合了邊界條件的間接邊界元理論的Helmholtz積分方程為：

其中：ρ為介質(zhì)密度；nX為聲場(chǎng)邊界場(chǎng)點(diǎn)的法線方向與nY方向規(guī)定相同；）為聲場(chǎng)邊界的速度邊界條件。

間接邊界元理論的Helmholtz積分方程離散后，邊界單元上未知的主要變量可以用形函數(shù)插值單元節(jié)點(diǎn)上的值來(lái)表示，但是其不能像直接邊界元那樣利用配點(diǎn)法來(lái)計(jì)算邊界表面的未知變量，而普遍的求解方法是利用變分原理來(lái)進(jìn)行邊界表面未知變量的求解。比如對(duì)于Neumann邊界條件，最終得到的泛函方程為：

其中：Q為系統(tǒng)矩陣，其為對(duì)稱陣；μ為結(jié)構(gòu)表面節(jié)點(diǎn)上未知變量的列向量；b為邊界上速度推導(dǎo)的列向量。

利用變分原理，公式（12）可得：

對(duì)于Neumann邊界條件，可通過(guò)求解方程（13）得到結(jié)構(gòu)邊界表面的聲壓差，從而求得場(chǎng)點(diǎn)聲壓、輻射聲功率以及聲輻射效率等物理量。

2.3 有限元方法或有限元與邊界元結(jié)合方法計(jì)算浸水結(jié)構(gòu)振動(dòng)及聲輻射特性運(yùn)動(dòng)方程

對(duì)于浸水結(jié)構(gòu)振動(dòng)及聲輻射特性的求解，無(wú)論是利用FEM或FEM/BEM結(jié)合理論，終歸于對(duì)耦合的運(yùn)動(dòng)方程進(jìn)行求解，耦合流體后的結(jié)構(gòu)運(yùn)動(dòng)方程為：

其中：Fa為流固交界面作用在結(jié)構(gòu)上流體動(dòng)壓力矢量，流體動(dòng)壓力Fa是一個(gè)復(fù)數(shù)，從物理上解釋體現(xiàn)為流體對(duì)結(jié)構(gòu)的作用力，其可以等效為對(duì)結(jié)構(gòu)的慣性力以及阻尼力的影響，其中的慣性力也就是對(duì)結(jié)構(gòu)的附加質(zhì)量效應(yīng)（對(duì)應(yīng)Fa的實(shí)部），阻尼力也就是對(duì)結(jié)構(gòu)的阻尼效應(yīng)（對(duì)應(yīng)Fa的虛部）。在激勵(lì)力頻率較低的低頻區(qū)域，流體對(duì)結(jié)構(gòu)的作用主要表現(xiàn)為附加質(zhì)量效應(yīng)，對(duì)結(jié)構(gòu)的阻尼的影響很小。也就是說(shuō)作用于結(jié)構(gòu)濕表面的流體動(dòng)壓力可以近似轉(zhuǎn)換為結(jié)構(gòu)動(dòng)力學(xué)方程的加速度項(xiàng)。這樣如果只考慮激勵(lì)力的頻率是在中低頻的范圍內(nèi)，那么最終水下結(jié)構(gòu)的動(dòng)力學(xué)方程可轉(zhuǎn)換為如下的形式：

其中：Ma為浸水結(jié)構(gòu)的附加質(zhì)量陣。當(dāng)研究進(jìn)水結(jié)構(gòu)固有振動(dòng)特性時(shí)Fs為零，忽略阻尼，上式對(duì)應(yīng)的廣義特征值問題為：

3 國(guó)內(nèi)浸水結(jié)構(gòu)振動(dòng)及聲輻射特性研究現(xiàn)狀

3.1 浸水結(jié)構(gòu)振動(dòng)特性研究現(xiàn)狀

對(duì)于浸水結(jié)構(gòu)振動(dòng)特性的研究，包括浸水結(jié)構(gòu)低頻振動(dòng)時(shí)的附加質(zhì)量陣、浸水結(jié)構(gòu)在流固耦合作用下的振動(dòng)特性計(jì)算以及相關(guān)實(shí)驗(yàn)的研究。

3.1.1 浸水結(jié)構(gòu)低頻振動(dòng)時(shí)的附加質(zhì)量的研究

對(duì)浸水結(jié)構(gòu)的附加質(zhì)量陣Ma的研究是求解浸水結(jié)構(gòu)振動(dòng)及聲輻射特性的基礎(chǔ)，其能體現(xiàn)出流固耦合的相互作用效應(yīng)，Ma的求解過(guò)程也是求解浸水結(jié)構(gòu)振動(dòng)及聲輻射特性的難點(diǎn)所在。對(duì)于研究附加質(zhì)量陣的計(jì)算方法，可以用有限元的理論直接將流體域進(jìn)行離散，也可以用邊界元的理論將三維流體域降一維再對(duì)邊界進(jìn)行離散計(jì)算。無(wú)論是哪種對(duì)于流體的處理方法，Ma的推導(dǎo)思路基本相同，首先是對(duì)Helmholtz微分方程公式（4）或Laplace微分方程公式（5）進(jìn)行三維或二維域內(nèi)積分并變換得到三維有限域積分方程或邊界積分方程，將流體動(dòng)壓力Fa用其積分后的式子來(lái)表示，此時(shí)流體動(dòng)壓力Fa是ω2的函數(shù)，這樣Fa就可以與（14）式中的質(zhì)量陣M合并，此時(shí)稱Fa中可以與質(zhì)量陣M合并的矩陣即為附加質(zhì)量矩陣Ma，它的意義是接觸面上結(jié)構(gòu)振動(dòng)加速度引起流體壓力的附加值。也就是說(shuō)對(duì)于附加質(zhì)量陣的研究，就是對(duì)于流體動(dòng)壓力Fa的研究。

將無(wú)限流體域截?cái)酁橛邢蘖黧w域，將有限域流體和結(jié)構(gòu)都利用有限元進(jìn)行數(shù)值模擬，即結(jié)構(gòu)有限元耦合流體有限元（FEM/FEM），這種方法基本思想是將公式（4）或公式（5）在有限體積域內(nèi)積分，將積分后的表達(dá)式進(jìn)而轉(zhuǎn)化成流體動(dòng)壓力Fa的形式，進(jìn)而求得附加質(zhì)量矩陣Ma，這是比較常用的處理流固耦合問題附加質(zhì)量的方法，其難點(diǎn)在于流體截?cái)嘤虻倪x擇、流體有限元的劃分以及流體單元類型的選擇等。汪庠寶和韓繼文等[2-3]對(duì)懸臂板和連續(xù)板的流固耦合振動(dòng)，給出了浸水結(jié)構(gòu)的附加質(zhì)量有限元法的計(jì)算格式，其對(duì)于板格在不同的邊界條件約束下的附加質(zhì)量系數(shù)進(jìn)行了研究。裴智勇和吳衛(wèi)國(guó)等[4]對(duì)高速船艙壁加筋板流固耦合振動(dòng)問題，將結(jié)構(gòu)有限元法和流體有限元法結(jié)合，根據(jù)Garlerkin法和Hamilton變分原理分別推導(dǎo)出離散的流體和結(jié)構(gòu)的運(yùn)動(dòng)方程，推導(dǎo)出流體附加質(zhì)量矩陣。文獻(xiàn)[5-6]同樣利用有限元方法離散流體的思想，早期研究人員多基于有限元方程編制相應(yīng)數(shù)值計(jì)算程序進(jìn)行計(jì)算。對(duì)于這類方法，可用于幾何形狀相對(duì)較簡(jiǎn)單結(jié)構(gòu)附加質(zhì)量陣的求解，其優(yōu)點(diǎn)是方程的理論推導(dǎo)及數(shù)值實(shí)現(xiàn)相對(duì)簡(jiǎn)潔，對(duì)一些有可能存在規(guī)律性結(jié)果的問題進(jìn)行研究比較有益；其缺點(diǎn)是一旦結(jié)構(gòu)比較復(fù)雜，其流體劃分有限元就比較復(fù)雜，單元和節(jié)點(diǎn)增多，導(dǎo)致計(jì)算量比較大。

對(duì)流體動(dòng)壓力Fa來(lái)說(shuō)，針對(duì)理想流體流場(chǎng)可用速度勢(shì)來(lái)表示，域中速度勢(shì)滿足Laplace方程，對(duì)于這類流場(chǎng)可使用源匯分布法（Green函數(shù)法）討論結(jié)構(gòu)物與周圍流體的穩(wěn)態(tài)相互作用問題。可求出流場(chǎng)中任一點(diǎn)的速度勢(shì)，進(jìn)而求出流體的動(dòng)壓力，從而求解附加質(zhì)量陣。王國(guó)強(qiáng)和楊建民[7]利用結(jié)構(gòu)有限元模型，結(jié)構(gòu)濕表面上布置源匯來(lái)模擬流體，采用線性三角形源匯法與結(jié)構(gòu)有限元模型相結(jié)合的辦法計(jì)算出了振動(dòng)結(jié)構(gòu)的附加質(zhì)量陣。鄭治國(guó)和趙德有等[8]對(duì)利用結(jié)構(gòu)有限元結(jié)合Green函數(shù)的邊界積分法計(jì)算結(jié)構(gòu)-流體的耦合作用，計(jì)算出非對(duì)稱的附加質(zhì)量矩陣，其利用的邊界積分方程是基于源匯法而得出的。錢勤等[9]借助Green第二恒等式，將浸水結(jié)構(gòu)的附加質(zhì)量陣計(jì)算通過(guò)將速度勢(shì)化為邊界上的積分方程來(lái)計(jì)算，即邊界積分公式。然后其用邊界元數(shù)值方法將邊界離散，從而計(jì)算出節(jié)點(diǎn)速度勢(shì)，然后將速度勢(shì)回帶計(jì)算出附加質(zhì)量陣。文獻(xiàn)[10-13]同樣利用源匯的思想，將流體運(yùn)用布置在結(jié)構(gòu)濕表面的源匯進(jìn)行模擬，其優(yōu)點(diǎn)是可以將流體從三維降至二維進(jìn)行處理；但是由于其在源匯的布置方面以及物體表面源匯強(qiáng)度方面不連續(xù)等缺點(diǎn)，其在數(shù)值方法求解的時(shí)候比較繁瑣，近年來(lái)使用這種方法解決附加質(zhì)量陣問題相對(duì)較少。

將三維無(wú)限流體域降維轉(zhuǎn)化為結(jié)構(gòu)邊界面上二維形式來(lái)處理，可以利用結(jié)構(gòu)有限元耦合流體直接邊界元的方法來(lái)求解附加質(zhì)量陣，這種方法基本思想是將公式（4）或公式（5）在結(jié)構(gòu)邊界內(nèi)做積分變換，將流體動(dòng)壓力Fa用（6）、（7）式或（8）式變換后的形式表示，進(jìn)而求得附加質(zhì)量矩陣Ma，其難點(diǎn)在于邊界積分方程的推導(dǎo)和變換以及邊界元所涉及的奇異積分處理問題。鄒元杰和趙德有[14]針對(duì)水下結(jié)構(gòu)聲固耦合振動(dòng)的特征值計(jì)算問題進(jìn)行了研究，其水下結(jié)構(gòu)附加質(zhì)量陣的計(jì)算是基于直接邊界元法推導(dǎo)的，而且指出假設(shè)鑲嵌在無(wú)限大剛性障板上的結(jié)構(gòu)，直接邊界元理論的積分公式可進(jìn)一步簡(jiǎn)化為Rayleigh積分來(lái)計(jì)算。付喜華和洪明等[15-17]對(duì)于可壓縮流體和不可壓縮流體的附加質(zhì)量陣進(jìn)行了研究，并利用無(wú)限大剛性障板的Rayleigh積分公式推導(dǎo)出了可壓縮流體的動(dòng)壓力，進(jìn)而推導(dǎo)出了其附加質(zhì)量矩陣，同時(shí)給出了不可壓縮流體的附加質(zhì)量矩陣。施衛(wèi)華[18]介紹了半無(wú)限流體域內(nèi)聲場(chǎng)的邊界元解法，從半無(wú)限流體域邊界積分方程出發(fā)，推導(dǎo)得到了板結(jié)構(gòu)鑲嵌在無(wú)限大剛性障板中時(shí)的Rayleigh積分方程，并求取了板式結(jié)構(gòu)在水中的附加水質(zhì)量。

魏建輝和陳美霞等[19]采用基于直接邊界元基本理論的Helmholtz積分方程，利用單層圓柱殼在空氣中的固有頻率和固有振型，在模態(tài)坐標(biāo)下推導(dǎo)出了結(jié)構(gòu)附加質(zhì)量的計(jì)算公式。安小同和洪明等[20]利用結(jié)構(gòu)有限元原理和流體直接邊界元原理對(duì)于浸水結(jié)構(gòu)的附加質(zhì)量陣進(jìn)行了研究，給出了適用于低頻作用下的基于直接邊界元的附加質(zhì)量陣的計(jì)算形式。文獻(xiàn)[21-23]同樣利用邊界元理論對(duì)流體附加質(zhì)量進(jìn)行處理。綜上可以看出Rayleigh積分公式可以用來(lái)處理假設(shè)鑲嵌在無(wú)限大剛性障板上的結(jié)構(gòu)（如局部加筋板），其優(yōu)點(diǎn)是較直接邊界元方法編程簡(jiǎn)單，其缺點(diǎn)是處理問題存在局限性；直接邊界元方法其優(yōu)點(diǎn)在于在理論上可以處理任意形狀的與介質(zhì)單面接觸的封閉結(jié)構(gòu)，但是其缺點(diǎn)在于對(duì)封閉結(jié)構(gòu)的特征頻率處解不唯一。無(wú)論是利用Rayleigh積分公式還是直接邊界元方法處理浸水結(jié)構(gòu)的附加質(zhì)量陣Ma，其優(yōu)點(diǎn)都在于將流體進(jìn)行了降維處理，可以有效減少計(jì)算量。

對(duì)于可壓縮流體和不可壓縮流體的附加質(zhì)量陣Ma，張升明[23]利用結(jié)構(gòu)有限元和流體直接邊界元方法，分別對(duì)流體為可壓縮和不可壓縮時(shí)，其對(duì)于結(jié)構(gòu)振動(dòng)特性的影響進(jìn)行了研究。其指出：將流體視為可壓縮時(shí)流體內(nèi)各質(zhì)點(diǎn)的運(yùn)動(dòng)相位是有差異的，運(yùn)動(dòng)相位差的存在說(shuō)明了阻尼的存在，從而耗散能量；將流體視為不可壓縮，流場(chǎng)內(nèi)各質(zhì)點(diǎn)運(yùn)動(dòng)是同相位的，因而它不產(chǎn)生阻尼。因此流體是否可壓，就相當(dāng)于是否考慮流體的阻尼。在低頻時(shí)，流體壓縮性對(duì)流體結(jié)構(gòu)的振動(dòng)固有頻率影響很小，可以忽略。因此在低頻時(shí)分析流體的附加質(zhì)量，可以把流體作為不可壓縮的介質(zhì)；而無(wú)論低頻或者高頻，在研究分析結(jié)構(gòu)在流體中的振動(dòng)響應(yīng)時(shí)，必須把流體作為可壓縮的介質(zhì)來(lái)處理。毛程亮和洪明等[24]針對(duì)浸水結(jié)構(gòu)動(dòng)力特性問題進(jìn)行了計(jì)算，表明在低頻范圍內(nèi)選取附加質(zhì)量，流體的可壓縮性對(duì)浸水結(jié)構(gòu)的固有頻率的影響非常小，在高頻范圍內(nèi)隨著計(jì)算頻率的逐漸增加，流體的可壓縮性對(duì)浸水結(jié)構(gòu)的固有頻率的影響將逐漸增大。

關(guān)于利用商業(yè)軟件計(jì)算浸水結(jié)構(gòu)的附加質(zhì)量陣Ma，吳紹亮和金咸定[25]利用流體直接邊界元理論得到流體的附加質(zhì)量陣，對(duì)大型通用軟件（MSC/NASTRAN）的流固耦合計(jì)算方法進(jìn)行了應(yīng)用開發(fā)。用MSC/PATRAN建立有限元模型，在MSC/NASTRAN平臺(tái)上主要對(duì)流體直接邊界元法、經(jīng)驗(yàn)公式法進(jìn)行結(jié)構(gòu)的流體附加質(zhì)量計(jì)算的對(duì)比研究和計(jì)算。這種方法的優(yōu)點(diǎn)在于可以對(duì)大型復(fù)雜結(jié)構(gòu)的附加質(zhì)量進(jìn)行求取，缺點(diǎn)在于有限元方法處理流體導(dǎo)致的單元以及節(jié)點(diǎn)增多，計(jì)算量增大。

3.1.2 浸水結(jié)構(gòu)在流固耦合作用下的振動(dòng)特性的研究

對(duì)于浸水結(jié)構(gòu)在流固耦合作用下的振動(dòng)特性的研究，普遍的研究是在解決附加質(zhì)量陣Ma的基礎(chǔ)上來(lái)進(jìn)行的，研究人員普遍關(guān)心的問題是結(jié)構(gòu)在流體中的振動(dòng)模態(tài)。

無(wú)論是利用何種方法處理流體動(dòng)壓力Fa，對(duì)于浸水結(jié)構(gòu)振動(dòng)特性的研究其最終都需要對(duì)廣義特征值方程（16）進(jìn)行求解，模態(tài)疊加法是求解廣義特征值方程的一種有效的方法。一般流固耦合系統(tǒng)的模態(tài)分析，根據(jù)計(jì)算過(guò)程的不同，可分為濕模態(tài)法和干模態(tài)法。

所謂干模態(tài)法是先求出結(jié)構(gòu)在真空中較低的若干階固有頻率和固有模態(tài)（稱為干模態(tài)），即先對(duì)運(yùn)動(dòng)方程（3）進(jìn)行求解，得到干模態(tài)；然后將干模態(tài)進(jìn)行模態(tài)坐標(biāo)變換考慮流體的影響，利用求得的干模態(tài)結(jié)合耦合運(yùn)動(dòng)方程（16）進(jìn)行求解，得到結(jié)構(gòu)在流體中的固有頻率和固有模態(tài)，其難點(diǎn)在于用何種方法如何求解方程（16）。汪庠寶等[2-3]給出了懸臂板和連續(xù)板的流固耦合振動(dòng)的有限元計(jì)算格式的結(jié)構(gòu)—流體系統(tǒng)耦合的無(wú)阻尼自由振動(dòng)方程式，利用干模態(tài)法計(jì)算了懸臂版和連續(xù)板的流固耦合自由振動(dòng)，給出了懸臂版在不同浸水深度下的固有頻率。楊建民和汪庠寶[27]提出了一種計(jì)及流體與結(jié)構(gòu)相互作用時(shí)船體結(jié)構(gòu)模態(tài)計(jì)算的方法，船體結(jié)構(gòu)用有限元法模擬，流體運(yùn)動(dòng)則用布置在船體表面的線性分布的三角形元源匯模擬。文中同時(shí)介紹了“濕模態(tài)法”以及“干模態(tài)法”計(jì)算流固耦合系統(tǒng)模態(tài)，并對(duì)兩種方法進(jìn)行了分析。徐張明和華宏星[42]利用FEM/I-BEM結(jié)合建立了結(jié)構(gòu)—聲耦合方程計(jì)算結(jié)構(gòu)表面的振動(dòng)特性，將結(jié)構(gòu)的物理坐標(biāo)用模態(tài)坐標(biāo)來(lái)表示，對(duì)耦合的結(jié)構(gòu)動(dòng)力學(xué)方程及聲場(chǎng)方程進(jìn)行解耦。并分別對(duì)浸沒于空氣中和水中的橢圓球殼進(jìn)行計(jì)算，其計(jì)算得出的結(jié)構(gòu)表面振動(dòng)速度表明，流體介質(zhì)改變了結(jié)構(gòu)的共振特性。文獻(xiàn)[28-30]也利用了干模態(tài)的方法求解了浸水結(jié)構(gòu)的固有振動(dòng)特性，可見，干模態(tài)法較之濕模態(tài)法有一個(gè)明顯的優(yōu)點(diǎn)，即可分別利用現(xiàn)成的結(jié)構(gòu)力學(xué)及流體力學(xué)領(lǐng)域中的最新研究成果來(lái)求解。同時(shí)，通過(guò)分析可知，較之濕模態(tài)法，干模態(tài)法中求解濕模態(tài)方程的自由度數(shù)目可大大縮小，即計(jì)算量大大減少；其缺點(diǎn)在于必須先求出結(jié)構(gòu)的干模態(tài)，才能進(jìn)一步求解結(jié)構(gòu)的濕模態(tài)，過(guò)程較為繁復(fù)，且干模態(tài)階數(shù)的選取在一定程度上影響求解精度。

所謂濕模態(tài)法是求解系統(tǒng)固有頻率和固有模態(tài)時(shí)直接計(jì)及周圍流體的動(dòng)力影響對(duì)耦合運(yùn)動(dòng)方程（16）進(jìn)行求解的方法。裴智勇和吳衛(wèi)國(guó)等[4]將結(jié)構(gòu)有限元法和流體有限元法結(jié)合起來(lái)，求得了高速船艙壁加筋板流固耦合系統(tǒng)的流體附加質(zhì)量矩陣，利用濕模態(tài)法對(duì)高速船艙壁加筋板流固耦合系統(tǒng)進(jìn)行了模態(tài)分析。鄭治國(guó)和趙德有等[8]采用20節(jié)點(diǎn)三維等參元離散結(jié)構(gòu)，在結(jié)構(gòu)表現(xiàn)上布置源匯，應(yīng)用Green函數(shù)的邊界積分法計(jì)算流體—結(jié)構(gòu)的相互影響，推導(dǎo)出系統(tǒng)的運(yùn)動(dòng)方程，用濕模態(tài)法對(duì)螺旋槳在水中的動(dòng)力特性進(jìn)行了計(jì)算。劉連海和洪明等[10，26]對(duì)流固耦合問題中結(jié)構(gòu)采用有限元和流體采用邊界元耦合的廣義特征方程求解方法進(jìn)行了研究，其指出由于流體附加質(zhì)量矩陣的存在，由此耦合而成系統(tǒng)的總質(zhì)量矩陣具有不對(duì)稱性（這里針對(duì)的是直接邊界元模型，對(duì)于間接邊界元模型來(lái)說(shuō)其是對(duì)稱的）且是滿陣，給特征值問題的求解帶來(lái)了困難，其利用的雙正交Lanczos法在一定程度上克服了上述困難。付喜華和洪明等[15-17]在求解流固耦合方程的特征值問題時(shí)，采用了迭代的Lanczos算法，即Lanczos算法與迭代算法結(jié)合的方法求解不對(duì)稱頻率相關(guān)特征值問題，安小同等[20]也利用這種方法對(duì)流固耦合方程的特征值問題進(jìn)行了求解。文獻(xiàn)[37-41，43-47]也利用了濕模態(tài)的方法求解了浸水結(jié)構(gòu)的固有振動(dòng)特性，濕模態(tài)方法較干模態(tài)方法優(yōu)點(diǎn)在于其可以直接對(duì)浸水結(jié)構(gòu)的固有模態(tài)進(jìn)行求解，可以直觀體現(xiàn)流體和結(jié)構(gòu)的相互作用；其缺點(diǎn)在于計(jì)算量往往較干模態(tài)大。但近年來(lái)，由于數(shù)值計(jì)算水平的發(fā)展研究者一般偏重直接求取浸水結(jié)構(gòu)的模態(tài)，即濕模態(tài)。

無(wú)論是利用干模態(tài)的方法還是濕模態(tài)的方法求解浸水結(jié)構(gòu)的固有振動(dòng)特性，終歸是要對(duì)方程（16）進(jìn)行求解，由有限元理論可知，方程（16）中質(zhì)量陣M是對(duì)稱矩陣，對(duì)稱矩陣在計(jì)算機(jī)數(shù)值計(jì)算時(shí)是比較節(jié)省存儲(chǔ)空間的，而附加質(zhì)量矩陣Ma的對(duì)稱性取決于所選用的對(duì)于流體的處理方法，文獻(xiàn)[2-6]所利用的有限元處理流體的方法，其得到的附加質(zhì)量矩陣Ma是非對(duì)稱的；文獻(xiàn)[19-23]所利用的直接邊界元處理流體的方法，其得到的附加質(zhì)量矩陣Ma同樣是非對(duì)稱的；文獻(xiàn)[42，46]所利用的間接邊界元處理流體的方法，其得到的附加質(zhì)量矩陣Ma是對(duì)稱的。可見利用間接邊界元方法處理流體較直接邊界元方法和有限元方法更能節(jié)省儲(chǔ)存空間。所以，求解浸水結(jié)構(gòu)的固有振動(dòng)特性，其難點(diǎn)之一便是如何解決附加質(zhì)量矩陣Ma的非對(duì)稱性問題。

朱曉芳和金咸定等[31]利用有限元軟件PATRAN/NASTRAN建立浸入水中圓柱體模型，并對(duì)其固有頻率進(jìn)行計(jì)算，分析比較了三種常用的浸水柱狀結(jié)構(gòu)固有頻率的計(jì)算方法：干模態(tài)的經(jīng)驗(yàn)修正系計(jì)算方法、流體力學(xué)兩維理論計(jì)算方法和流固耦合數(shù)值計(jì)算方法。文獻(xiàn)[32-36]同樣論述了利用有限元軟件PATRAN/NASTRAN或ANSYS來(lái)建立浸水結(jié)構(gòu)模型，并對(duì)其固有頻率進(jìn)行數(shù)值計(jì)算，該方法已用于單

層和雙層環(huán)肋增強(qiáng)圓柱殼的動(dòng)力特性計(jì)算。此方法優(yōu)點(diǎn)是可用于任意復(fù)雜浸水結(jié)構(gòu)動(dòng)力特性分析；缺點(diǎn)是計(jì)算節(jié)點(diǎn)多，計(jì)算量大。

3.1.3 相關(guān)實(shí)驗(yàn)研究

相對(duì)于數(shù)值模擬來(lái)說(shuō)，國(guó)內(nèi)對(duì)于實(shí)驗(yàn)驗(yàn)證數(shù)值模擬準(zhǔn)確性的研究以及試驗(yàn)的研究相對(duì)較少。

裴智勇和吳衛(wèi)國(guó)等[4]用濕模態(tài)法求得流體附加質(zhì)量矩陣，依據(jù)結(jié)構(gòu)有限元和流體有限元理論建立的加筋板結(jié)構(gòu)流固耦合動(dòng)力控制方程，編制了加筋板流固耦合振動(dòng)計(jì)算程序。其為了認(rèn)識(shí)高速船艙壁加筋板流固耦合動(dòng)力性能，建立了船舶的局部模擬艙段結(jié)構(gòu)，并進(jìn)行模態(tài)試驗(yàn)，研究了艙內(nèi)不同水深時(shí)艙壁加筋板結(jié)構(gòu)流固耦合振動(dòng)模態(tài)，同時(shí)與有限元計(jì)算結(jié)果進(jìn)行了比較。

王杰德和楊永謙等[29]利用其對(duì)于船體薄壁梁彎扭動(dòng)力分析的流固耦合數(shù)值計(jì)算結(jié)果與鋁制船模在空氣中和水中的實(shí)驗(yàn)得到的固有頻率進(jìn)行了對(duì)比，驗(yàn)證了其數(shù)值結(jié)果的可行性。謝志勇和周其斗等[41]采用鋼質(zhì)雙層柱殼結(jié)構(gòu)模擬潛艇結(jié)構(gòu)，利用有限元軟件建立數(shù)值分析模型并進(jìn)行了相應(yīng)的模態(tài)數(shù)值計(jì)算，之后分別對(duì)懸掛于空氣中的柱殼：殼間充滿空氣的模態(tài)和殼間充水的模態(tài)兩類情況進(jìn)行了模態(tài)試驗(yàn)測(cè)量，并將試驗(yàn)結(jié)果與有限元法所得的計(jì)算結(jié)果進(jìn)行了對(duì)比。

楊吉新和張可等[36]將用有限元軟件ANSYS建模計(jì)算的懸臂柱在不同水位下的模態(tài)計(jì)算結(jié)果和實(shí)驗(yàn)進(jìn)行了對(duì)比，實(shí)驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果的差值小于10%，說(shuō)明了這種數(shù)值方法是可行的。

童宗鵬和王國(guó)治[48]采用模態(tài)試驗(yàn)，利用單點(diǎn)激勵(lì)多點(diǎn)響應(yīng)的方式，對(duì)整個(gè)艦艇殼模型進(jìn)行了測(cè)量，通過(guò)對(duì)艦艇結(jié)構(gòu)傳遞函數(shù)的測(cè)量及試驗(yàn)數(shù)據(jù)的處理，得到了固有頻率和振型，并與有限元的計(jì)算結(jié)果進(jìn)行了對(duì)比。

胡選利等[49]對(duì)于聲振耦合目前求解的三種常用方法：聲彈性法、有限元法和邊界元法進(jìn)行了分析，對(duì)采用各種簡(jiǎn)化假設(shè)時(shí)三種方法的差異及統(tǒng)一進(jìn)行了討論。通過(guò)對(duì)矩形截面彈性薄壁箱體振動(dòng)進(jìn)行的實(shí)驗(yàn)測(cè)量，對(duì)各種理論模型所用簡(jiǎn)化假設(shè)的合理性進(jìn)行了驗(yàn)證。

姚熊亮等[50]從實(shí)驗(yàn)和數(shù)值計(jì)算兩個(gè)方面，對(duì)敷有隔聲去耦材料的水下加筋圓柱殼體振動(dòng)進(jìn)行討論分析。其實(shí)驗(yàn)地點(diǎn)選擇在大型消音水池進(jìn)行，在實(shí)驗(yàn)中測(cè)量了模型在水中的固有頻率（濕模態(tài)），采用的方法是尋找響應(yīng)曲線峰值點(diǎn)所對(duì)應(yīng)的頻率值，即共振頻率。

對(duì)于浸水結(jié)構(gòu)動(dòng)力特性的實(shí)驗(yàn)研究相對(duì)較少，主要是由于其操作起來(lái)很復(fù)雜，如水下的實(shí)驗(yàn)器械選取以及布置安裝等困難；而且其受到實(shí)驗(yàn)條件的影響極大，如用有限的水池模擬無(wú)限水域空間以及邊界條件的選取等難題。但是，基于數(shù)值算法的驗(yàn)證、修正以及問題的發(fā)現(xiàn)和方法的改進(jìn)，這類實(shí)驗(yàn)的開展還是十分有必要的。

3.2 浸水結(jié)構(gòu)聲輻射特性研究現(xiàn)狀

3.2.1 弱耦合作用下的聲輻射問題

對(duì)于結(jié)構(gòu)在空氣中（輕流體）的聲輻射問題，一般視之為弱耦合，可忽略耦合作用。但是空氣中結(jié)構(gòu)聲輻射問題的研究方法是研究浸水結(jié)構(gòu)聲輻射問題的基礎(chǔ)，其對(duì)于流固耦合聲輻射問題的研究，尤其是對(duì)于耦合方程中的質(zhì)量陣，是非常具有借鑒意義的。

利用直接邊界元方法來(lái)求解結(jié)構(gòu)在輕流體中的聲輻射問題是比較常用的方法，由于其不用考慮流固耦合的相互作用，可直接利用有限元方法對(duì)結(jié)構(gòu)運(yùn)動(dòng)方程（2）求解得出結(jié)構(gòu)在空氣中的響應(yīng)，將其作為求解聲輻射問題的已知條件，利用邊界積分公式（6）或（7）變換得到未知變量的方程組后，進(jìn)而利用配點(diǎn)法等對(duì)結(jié)構(gòu)聲輻射特性進(jìn)行求解。利用間接邊界元方法來(lái)求解結(jié)構(gòu)在輕流體中的聲輻射問題是近年來(lái)研究的熱點(diǎn)，若不用考慮流固耦合的相互作用，同樣可直接利用有限元方法對(duì)結(jié)構(gòu)運(yùn)動(dòng)方程（2）求解得出結(jié)構(gòu)在空氣中的響應(yīng)，將其作為求解聲輻射問題的已知條件，進(jìn)而利用邊界積分公式（10）進(jìn)行變換計(jì)算，若聲場(chǎng)滿足Neumann邊界條件則得到邊界積分公式（11），利用公式（11）得到泛函方程（12），應(yīng)用變分原理得到方程（13），進(jìn)而求解結(jié)構(gòu)聲輻射特性。無(wú)論是直接或間接邊界元方法其難點(diǎn)之一是對(duì)于各類奇異積分的處理。趙翔等[51]首先建立求解自由場(chǎng)中結(jié)構(gòu)體聲輻射的數(shù)學(xué)物理模型，從Helmholtz方程出發(fā)，用加權(quán)余量法推導(dǎo)Helmholtz方程的邊界積分方程，即直接邊界元的Helmholtz邊界積分方程，然后用線性三角形單元離散邊界，得到了邊界積分的數(shù)值計(jì)算公式。編制了用于計(jì)算任意形狀結(jié)構(gòu)體聲源的三維輻射聲場(chǎng)的程序。對(duì)旋轉(zhuǎn)式冰箱壓縮機(jī)的輻射噪聲場(chǎng)進(jìn)行計(jì)算，并與實(shí)測(cè)結(jié)果進(jìn)行了比較。黎勝等[52]采用結(jié)構(gòu)有限元方法和流體邊界元方法建立了輕流體介質(zhì)中復(fù)雜結(jié)構(gòu)振動(dòng)聲輻射的模態(tài)分析模型，文中以板和加筋板結(jié)構(gòu)的振動(dòng)聲輻射為例進(jìn)行了模態(tài)分析，并與用遠(yuǎn)場(chǎng)Rayleigh積分計(jì)算所得結(jié)果進(jìn)行了比較。研究了激勵(lì)力作用位置對(duì)板結(jié)構(gòu)聲輻射的影響和加筋板與板的模態(tài)輻射效率特點(diǎn)，并指出板加筋后會(huì)改變板的模態(tài)輻射效率規(guī)律。劉洪林[53]詳細(xì)地推導(dǎo)出了對(duì)應(yīng)于無(wú)限域和半無(wú)限域的振動(dòng)聲輻射問題的邊界積分計(jì)算公式，并且比較了聲輻射問題的直接邊界元法和間接邊界元法。利用直接邊界元理論的積分公式處理結(jié)構(gòu)在空氣中的聲輻射問題，其可以對(duì)任意形狀的封閉結(jié)構(gòu)進(jìn)行場(chǎng)點(diǎn)聲學(xué)特性的計(jì)算，其優(yōu)點(diǎn)也在于此；缺點(diǎn)是其只能解決封閉結(jié)構(gòu)問題，對(duì)于開口結(jié)構(gòu)問題不能處理，而且其還存在特征頻率解不唯一問題。文獻(xiàn)[53]涉及到了間接邊界元方法處理結(jié)構(gòu)在空氣中的聲輻射問題，間接邊界元方法優(yōu)點(diǎn)在于可以處理任意形狀封閉結(jié)構(gòu)以及開口結(jié)構(gòu)的聲輻射問題；其缺點(diǎn)是推導(dǎo)較直接邊界元法繁瑣且數(shù)值處理復(fù)雜。

對(duì)于可以假定鑲嵌在無(wú)限大障板上的結(jié)構(gòu)，直接邊界元的積分公式可以簡(jiǎn)化為Rayleigh積分公式（8），利用Rayleigh積分公式處理聲輻射問題的難點(diǎn)同樣是對(duì)奇異積分的處理。郭新毅和洪明[54-56]采用結(jié)構(gòu)有限元和流體邊界元（Rayleigh積分）相結(jié)合的方法對(duì)損傷結(jié)構(gòu)振動(dòng)及聲輻射特性建立計(jì)算分析數(shù)值模型，得到結(jié)構(gòu)的輻射聲功率和輻射指向性，考慮到不同損傷存在形式，計(jì)算分析了損傷對(duì)振動(dòng)頻率、模態(tài)以及輻射聲功率和指向性的影響。文獻(xiàn)[57-61]同樣利用Rayleigh積分公式來(lái)處理結(jié)構(gòu)在空氣中的聲輻射問題，這種方法優(yōu)點(diǎn)是其由直接邊界元積分方程簡(jiǎn)化而來(lái)，積分公式相對(duì)簡(jiǎn)單，數(shù)值處理也相對(duì)簡(jiǎn)潔；其缺點(diǎn)在于其只能對(duì)于少部分可以假定鑲嵌在無(wú)限大障板上的結(jié)構(gòu)進(jìn)行計(jì)算。

3.2.2 浸水結(jié)構(gòu)的聲輻射問題

對(duì)于結(jié)構(gòu)在流固耦合下的聲輻射，一般研究的是結(jié)構(gòu)在重流體中的強(qiáng)耦合問題，是與浸水結(jié)構(gòu)在流固耦合作用下的振動(dòng)問題相輔相成的。無(wú)論是用那種方法處理結(jié)構(gòu)和流體之間的相互關(guān)系，其最終都?xì)w結(jié)到對(duì)于結(jié)構(gòu)表面的聲壓或聲壓差的求解，進(jìn)而計(jì)算浸水結(jié)構(gòu)聲功率和聲輻射效率等聲輻射相關(guān)物理量。

對(duì)于浸水結(jié)構(gòu)流固耦合作用下的聲輻射問題，如果浸水結(jié)構(gòu)的附加質(zhì)量陣Ma是按照結(jié)構(gòu)有限元耦合流體有限元的理論進(jìn)行研究，其聲輻射的求解可以直接利用有限元原理進(jìn)行研究。不計(jì)阻尼影響，結(jié)構(gòu)有限元耦合流體有限元模型（FEM/FEM）的系統(tǒng)運(yùn)動(dòng)方程如下：

其中：KS為結(jié)構(gòu)模型的剛度矩陣；MS為結(jié)構(gòu)模型的質(zhì)量矩陣；KC為耦合剛度矩陣；KF為流體模型的剛度矩陣；MS為流體模型的質(zhì)量矩陣；FSi為作用在結(jié)構(gòu)上的激勵(lì)力矢量；ui為節(jié)點(diǎn)位移；pi為節(jié)點(diǎn)動(dòng)壓力，即結(jié)構(gòu)表面聲壓。由于這種方法在單元離散的過(guò)程中，隨著計(jì)算頻率的升高單元的密度將大大增加，因此科研人員往往利用另一種方式求解，即先求出浸水結(jié)構(gòu)在流體中的振動(dòng)響應(yīng)，然后再利用所求出的振動(dòng)響應(yīng)結(jié)合邊界元理論計(jì)算出浸水結(jié)構(gòu)的聲輻射特性。

如果浸水結(jié)構(gòu)的附加質(zhì)量陣Ma是按照結(jié)構(gòu)有限元耦合流體邊界元理論（直接或間接）進(jìn)行研究，其聲輻射特性一般是將流體動(dòng)壓力Fa表示成結(jié)構(gòu)表面聲壓或聲壓差的函數(shù)，即將流體的動(dòng)壓力Fa用邊界積分方程（6）、（7）、（8）或（10）的形式來(lái)表示，進(jìn)而利用結(jié)構(gòu)表面聲壓或聲壓差與結(jié)構(gòu)表面法向振動(dòng)速度的關(guān)系，以及法向振動(dòng)速度與結(jié)構(gòu)表面振速及位移的關(guān)系，對(duì)耦合方程（15）進(jìn)行處理后聯(lián)立Fa，求解聯(lián)立方程組（18）或（19）得到浸水結(jié)構(gòu)的聲輻射特性。

不計(jì)阻尼影響，結(jié)構(gòu)有限元耦合流體直接邊界元模型（FEM/D-BEM）的系統(tǒng)運(yùn)動(dòng)方程如下：

其中：LC為耦合矩陣；T為法向速度和位移之間的轉(zhuǎn)換矩陣；Aω（）和Bω（）為直接邊界元影響矩陣。

不計(jì)阻尼影響，且滿足單層勢(shì)σ=0，則結(jié)構(gòu)有限元耦合流體間接邊界元模型（FEM/I-BEM）的系統(tǒng)運(yùn)動(dòng)方程如下：

其中：Qω（）為間接邊界元影響矩陣；ui為節(jié)點(diǎn)雙層勢(shì)（壓力跳動(dòng)量），即結(jié)構(gòu)表面聲壓差。

對(duì)于可以假定鑲嵌在無(wú)限大障板上的浸水結(jié)構(gòu)，直接邊界元的積分公式（6）或（7）可以簡(jiǎn)化為Rayleigh積分公式（8），將簡(jiǎn)化后的Rayleigh積分公式變換得到流體動(dòng)壓力Fa，將Fa代入耦合流體后的結(jié)構(gòu)有限元運(yùn)動(dòng)方程（14）進(jìn)而求出方程（15）聯(lián)立Fa，并求解聯(lián)立方程組（18）便可得到浸水結(jié)構(gòu)的聲輻射特性，其難點(diǎn)是有限元和Rayleigh積分方程的耦合以及奇異積分的處理。王其峰[12]利用Hess-Smith源匯方法對(duì)浸水結(jié)構(gòu)的附加質(zhì)量陣進(jìn)行了推導(dǎo)，得出考慮流體加載效應(yīng)的結(jié)構(gòu)有限元方程并求解得出解決聲輻射問題的邊界條件。對(duì)于流固耦合作用下的聲輻射問題，給出了一般問題的直接邊界元法的Helmholtz積分方程以及鑲嵌在無(wú)限大剛性障板上結(jié)構(gòu)的Rayleigh積分方程，并概述了輻射聲功率的多種計(jì)算方法。毛程亮[13]對(duì)流固耦合聲輻射問題的研究方法及思路與王其峰類似，對(duì)半無(wú)限空間中的“軟”“硬”邊界進(jìn)行了介紹。施衛(wèi)華[19]利用Rayleigh積分方法對(duì)浸水結(jié)構(gòu)的附加質(zhì)量陣進(jìn)行了推導(dǎo)，得出了考慮流體加載效應(yīng)的結(jié)構(gòu)有限元方程，并求解得出四邊簡(jiǎn)支板的振動(dòng)響應(yīng)。對(duì)于流固耦合作用下的聲輻射問題，其將四邊簡(jiǎn)支板假設(shè)為鑲嵌在無(wú)限大剛性障板上的結(jié)構(gòu)，利用Rayleigh積分方程結(jié)合之前求解的振動(dòng)響應(yīng)得到了板在空氣中和可壓縮流體中的聲輻射特性。文獻(xiàn)[71-72，74，83-85]同樣利用Rayleigh積分公式來(lái)處理浸水結(jié)構(gòu)的聲輻射問題，這種方法優(yōu)點(diǎn)在于積分公式相對(duì)簡(jiǎn)單，數(shù)值處理也較直接邊界元方法簡(jiǎn)潔，缺點(diǎn)在于其只能對(duì)于少部分可以假定鑲嵌在無(wú)限大障板上的單面浸水結(jié)構(gòu)進(jìn)行計(jì)算。

對(duì)于單面浸水的封閉式結(jié)構(gòu)，將直接邊界元的積分公式（6）或（7）變換得到流體動(dòng)壓力Fa，將Fa代入耦合流體后的結(jié)構(gòu)有限元運(yùn)動(dòng)方程（14）得到方程（15）并聯(lián)立Fa求解聯(lián)立方程組（18）便可得到浸水結(jié)構(gòu)的聲輻射特性，其難點(diǎn)在于非對(duì)稱耦合矩陣以及奇異積分的處理。鄒元杰和趙德有[43]利用結(jié)構(gòu)有限元和流體直接邊界元理論，建立耦合振動(dòng)方程，并結(jié)合淺水域的聲學(xué)邊界條件，研究了水深對(duì)聲輻射的影響。其對(duì)方箱結(jié)構(gòu)聲輻射進(jìn)行數(shù)值計(jì)算后，得出的方箱結(jié)構(gòu)在淺水中的聲輻射功率等計(jì)算結(jié)果表明：淺水中的結(jié)構(gòu)聲輻射特性對(duì)水深有較強(qiáng)的依賴性。韓峰和閆桂榮等[45]利用LMS virtual lab中的Acoustic BEM模塊實(shí)現(xiàn)了圓錐殼結(jié)構(gòu)FEM/BEM聲振耦合模型的仿真計(jì)算，基于結(jié)構(gòu)聲振耦合響應(yīng)的計(jì)算求解結(jié)果，得到了所關(guān)心位置處的聲壓級(jí)，并將其與混響噪聲實(shí)驗(yàn)測(cè)量結(jié)果進(jìn)行比較，驗(yàn)證了結(jié)構(gòu)聲振耦合計(jì)算的正確性和合理性。沈順根等[63]利用結(jié)構(gòu)有限元耦合流體邊界元，將結(jié)構(gòu)有限元和修正的Helmholtz邊界積分方程結(jié)合，建立了一種求解結(jié)構(gòu)聲輻射及物面動(dòng)響應(yīng)的矩陣分析方法。其對(duì)于流體的處理主要是利用直接邊界元的基本原理來(lái)進(jìn)行的。俞孟薩等[64]根據(jù)彈性結(jié)構(gòu)振動(dòng)和聲輻射理論，采用結(jié)構(gòu)有限元和流體邊界元，系統(tǒng)地分析了水下運(yùn)動(dòng)體典型結(jié)構(gòu)—有限長(zhǎng)彈性加肋圓柱殼振動(dòng)、耦合振動(dòng)和聲輻射的相似性，在公式推導(dǎo)的過(guò)程中強(qiáng)調(diào)了結(jié)構(gòu)阻尼和流體負(fù)載的作用。崔宏武和趙德有等[73]對(duì)水下三維彈性結(jié)構(gòu)振動(dòng)的聲輻射問題進(jìn)行了流體—結(jié)構(gòu)聲耦合振動(dòng)的數(shù)值計(jì)算。對(duì)結(jié)構(gòu)部分采用有限元法（FEM），流體部分采用直接邊界元法（D-BEM），推導(dǎo)并給出了基于流體直接邊界元法和結(jié)構(gòu)有限元方法協(xié)調(diào)的流體結(jié)構(gòu)聲耦合振動(dòng)聯(lián)立方程式，并討論了其求解思路。文中著重研究了邊界元法中矩陣元素的數(shù)值計(jì)算，編制了相應(yīng)的計(jì)算機(jī)程序，利用脈動(dòng)球算例驗(yàn)證了本文方法的正確性。文獻(xiàn)[75-76，79-81]同樣利用直接邊界元理論的積分公式處理浸水結(jié)構(gòu)的聲輻射問題，其優(yōu)點(diǎn)在于其可以解決任意形狀單面浸水封閉結(jié)構(gòu)的聲輻射問題；其缺點(diǎn)在于求解運(yùn)動(dòng)方程（15）和Fa聯(lián)立方程組的時(shí)候無(wú)論是對(duì)于非對(duì)稱矩陣的存儲(chǔ)還是對(duì)非對(duì)稱矩陣方程的求解都存在較大的困難，而且其不能處理開口結(jié)構(gòu)問題。

對(duì)于雙面浸水的封閉式結(jié)構(gòu)或開口結(jié)構(gòu)，將間接邊界元的積分公式（10）變換得到流體動(dòng)壓力Fa，將Fa代入耦合流體后的結(jié)構(gòu)有限元運(yùn)動(dòng)方程（14）得到方程（15）聯(lián)立Fa，并求解聯(lián)立方程組（19）便可得到雙面浸水浸水結(jié)構(gòu)的聲輻射特性。其難點(diǎn)在于有限元和間接邊界積分方程的耦合以及超奇異積分或超奇異積分轉(zhuǎn)化后的雙層弱奇異積分的處理。徐張明和華宏星等[42]將有限元和間接邊界元結(jié)合得到聲振耦合方程，進(jìn)而計(jì)算相應(yīng)的輻射聲場(chǎng)。其對(duì)一水下受均布諧載荷的球殼進(jìn)行數(shù)值計(jì)算來(lái)驗(yàn)證其方法的正確性，分別得到了球殼在空氣中和水中的輻射聲壓、輻射聲功率和輻射效率。這種方法的最大優(yōu)點(diǎn)是可以同時(shí)計(jì)算結(jié)構(gòu)表面振動(dòng)和輻射聲場(chǎng)，而且適用于任意復(fù)雜結(jié)構(gòu)。魏建輝和陳美霞等[46]基于結(jié)構(gòu)有限元和流體間接邊界元方法（FEM/I-BEM）研究了雙層圓柱殼的振動(dòng)及聲輻射，并與解析法及目前常用的數(shù)值算法：結(jié)構(gòu)有限元耦合結(jié)構(gòu)有限元，結(jié)構(gòu)有限元耦合直接邊界元方法進(jìn)行了對(duì)比，驗(yàn)證了本方法的正確性。指出了結(jié)構(gòu)有限元和流體間接邊界元計(jì)算方法可以很好地用于解決水下雙層圓柱殼的振動(dòng)及聲輻射問題。利用間接邊界元方法處理浸水結(jié)構(gòu)的聲輻射問題，其優(yōu)點(diǎn)在于可以處理任意形狀雙面浸水的封閉結(jié)構(gòu)以及開口結(jié)構(gòu)的聲輻射問題，而且其得到的耦合矩陣是對(duì)稱陣，對(duì)于矩陣的存儲(chǔ)可以節(jié)省內(nèi)存空間，并且對(duì)稱矩陣方程的求解較非對(duì)稱矩陣方程容易求解；其缺點(diǎn)是推導(dǎo)過(guò)程繁瑣復(fù)雜，較直接邊界元在奇異積分等數(shù)值處理上更為困難。綜上可見，無(wú)論是Rayleigh積分公式、直接邊界元法還是間接邊界元法，其都可以有效地結(jié)合有限元方法，將結(jié)構(gòu)模擬為有限元處理，將流體降維利用邊界元來(lái)處理，進(jìn)而直接求解流固耦合運(yùn)動(dòng)方程。

姚熊亮等[32，50]對(duì)雙殼體結(jié)構(gòu)的聲輻射特性進(jìn)行了研究，其主要是將有限元軟件ANSYS建模并計(jì)算的浸水雙殼體的振動(dòng)響應(yīng)帶入到了聲學(xué)軟件SYSNOISE利用其直接邊界元計(jì)算模塊來(lái)計(jì)算浸水雙殼體的振動(dòng)特性。文獻(xiàn)[33-35，38，40，62，65-70，77-78，82]同樣論述了利用有限元軟件ANSYS結(jié)合SYSNOISE來(lái)建立浸水結(jié)構(gòu)模型，對(duì)其聲輻射特性進(jìn)行數(shù)值計(jì)算，該方法已用于單層、雙層環(huán)肋增強(qiáng)圓柱殼等浸水結(jié)構(gòu)的聲輻射特性計(jì)算，此方法優(yōu)點(diǎn)是可用于任意復(fù)雜浸水結(jié)構(gòu)的聲輻射特性分析；缺點(diǎn)是其有限元和邊界元的計(jì)算是分開進(jìn)行的，其存在將有限元的計(jì)算結(jié)果導(dǎo)入邊界元程序的復(fù)雜過(guò)程，而且其將流體截?cái)酁橛邢拊M(jìn)行處理，若要得到較精確結(jié)果其計(jì)算單元和節(jié)點(diǎn)會(huì)較多進(jìn)而計(jì)算量會(huì)較大。

3.2.3 關(guān)于直/間接邊界元方法解的非唯一性問題

聲輻射問題一般采用Helmholtz微分方程在邊界條件下形成的定解方程來(lái)描述，其解具有唯一性，但是一旦采用Helmholtz邊界積分方程對(duì)上述定解方程進(jìn)行求解，在某些波數(shù)（以結(jié)構(gòu)封閉表面為界的內(nèi)部Dirichlet問題的特征頻率）處就會(huì)出現(xiàn)非唯一性的現(xiàn)象，產(chǎn)生非唯一性的原因是由于在這些波數(shù)處，解的邊界積分表述與原問題之間不存在對(duì)等性。

對(duì)于特征頻率處解的非唯一性問題，有兩種最具有代表性的方法，一種1968年Schenck年提出的CHIEF方法（Combined Helmholtz Integral Equation Formulation），另一種是1971年Burton與Miller提出的Burton-Miller法。眾多其它的方法都是在這兩種方法的基礎(chǔ)上對(duì)其進(jìn)行改進(jìn)而進(jìn)行的。施衛(wèi)華[19]對(duì)于無(wú)限域自由場(chǎng)中的三維封閉彈性體結(jié)構(gòu)，采用直接邊界元法計(jì)算結(jié)構(gòu)的聲輻射。利用CHIEF法解決特征頻率處解不唯一的問題。文獻(xiàn)[67，69]同樣采用了CHIEF法解決Helmholtz積分方程解的非唯一性問題。文獻(xiàn)[89-92]同樣利用CHIEF方法處理特征頻率處解不唯一的問題，其優(yōu)點(diǎn)在于，不在非常高的激勵(lì)頻率下，CHIEF點(diǎn)的選取并不困難，只要有一個(gè)CHIEF點(diǎn)有效，CHIEF方法就有效；其缺點(diǎn)在于，在高頻情形下，由于此時(shí)振型十分密集，節(jié)點(diǎn)的密度大大增加，CHIEF點(diǎn)的選取在振型節(jié)點(diǎn)上的概率也大大增加，選取合適的CHIEF點(diǎn)也更加困難。文獻(xiàn)[86-88]所使用的內(nèi)部Helmholtz方程關(guān)于內(nèi)點(diǎn)坐標(biāo)取導(dǎo)后構(gòu)成補(bǔ)充方程式方法與CHIEF方法相比較，能更有效地處理特征頻率處解不唯一的問題。

蔣偉和何正耀等[93]應(yīng)用了一種無(wú)界聲媒質(zhì)空間中任意形狀輻射體輻射聲場(chǎng)的新方法—綜合H-elmholtz積分公式法（Synthetic Helmholtz Integral Equation Formulation,SHIEF）。其推導(dǎo)是基于直接邊界元理論的，主要解決了當(dāng)輻射體的振動(dòng)頻率等于或接近其特征頻率時(shí)，聲輻射出現(xiàn)奇異解的問題。該方法將關(guān)于輻射體的內(nèi)部、表面和外部Helmholtz積分方程有機(jī)地組合在一起，在已知輻射體表面法向振速的條件下，求出輻射體表面的聲壓分布，進(jìn)而求得其輻射聲場(chǎng)分布。該方法將外部、表面以及內(nèi)部Helmholtz積分方程有機(jī)地結(jié)合起來(lái)，減少了對(duì)輻射體內(nèi)部點(diǎn)的依賴，并且充分利用內(nèi)外輻射聲場(chǎng)的信息，克服了CHIEF法內(nèi)部點(diǎn)選取的困難和超定外部Helmholtz積分公式法在特征頻率處和高頻段收斂速度慢以及所需外點(diǎn)數(shù)多的缺點(diǎn)。

3.2.4 關(guān)于直/間接邊界元奇異積分的處理

不論是直接邊界元理論的Helmholtz積分方程還是間接邊界元理論的Helmholtz積分方程，亦或是直接邊界元理論的Helmholtz積分方程的特殊形式Rayleigh積分方程，都存在奇異積分的問題，對(duì)于奇異積分的處理是解決Helmholtz積分方程或Rayleigh積分方程的重中之重也是難點(diǎn)所在。奇異積分根據(jù)邊界歸化的途徑不同，可能有三種類型：弱奇異O）；Cauchy型奇異O）；超奇異。超奇異積分又稱為Hadamard型奇異積分，其具有比Cauchy型奇異積分更高階的奇異性，根據(jù)經(jīng)典微積分學(xué)的概念，Hadamard型奇異積分是發(fā)散且沒有意義的，無(wú)法用經(jīng)典的數(shù)值積分公式計(jì)算出其具有一定精度的近似值，它們只有在Cauchy主值意義下才存在，邊界元法中的Hadamard型奇異積分實(shí)際上是其有限部分。

對(duì)于弱奇異積分的研究，劉連海[10]對(duì)于Rayleigh積分中弱奇異積分的處理利用的是退化單元的方法。對(duì)于四邊形單元，其將四邊形單元分解為若干三角形單元，對(duì)含奇異點(diǎn)的三角形單元進(jìn)行適當(dāng)?shù)淖鴺?biāo)變換以消除奇異性。柯兵和譚林森等[61]對(duì)于直接邊界元法存在的奇異性問題，借助于引進(jìn)退化單元，將三角形單元再分成三個(gè)三角形單元進(jìn)行處理，解決了基本解的弱奇異性問題。楊勇[69]在處理基于表面的Rayleigh積分方程的弱奇異積分時(shí)，對(duì)于4節(jié)點(diǎn)線性單元，其引入了退化元消除了單元上的積分奇異性。黎勝等[90]對(duì)于直接邊界元Helmholtz積分方程的積分奇異性處理，采用三角形斜坐標(biāo)系、退化單元方法。對(duì)于線性三角形單元采用三角形斜坐標(biāo)系來(lái)消除積分奇異性；對(duì)于四邊形線性元和四邊形二次元采用退化元來(lái)消除積分奇異性。李宏偉和文立華等[100]提出了采用三角形或四邊形改進(jìn)的二次插值函數(shù)，并結(jié)合局部極坐標(biāo)變換或引入退化單元的方法，解決了任意三維非光滑結(jié)構(gòu)表面聲學(xué)計(jì)算邊界元法中的奇異積分問題。其采用的三角形單元為六節(jié)點(diǎn)三角形單元，四邊形單元為八節(jié)點(diǎn)四邊形單元，但是其插值函數(shù)并非一般常見的二次插值函數(shù)，而是改進(jìn)的二次差值函數(shù)。其所說(shuō)的退化單元，是將三角形或四邊形單元分為兩個(gè)子單元來(lái)進(jìn)行奇異積分的計(jì)算。利用退化單元的方法來(lái)解決奇異積分問題，其優(yōu)點(diǎn)是對(duì)于處理弱奇異積分是非常有效的；缺點(diǎn)是其不能對(duì)高階奇異積分進(jìn)行處理，且對(duì)于單元的形狀也有一定要求。

劉洪林[53]指出對(duì)于弱奇異積分的數(shù)值計(jì)算，目前最有效的處理方法是極坐標(biāo)變換法；對(duì)于強(qiáng)奇異積分，可以首先將其變換為主值積分和奇點(diǎn)自由項(xiàng)的組合，然后通過(guò)適當(dāng)?shù)臄?shù)學(xué)變換（如配項(xiàng)處理、坐標(biāo)變換等）將其降奇后再進(jìn)行求解，而且由于這種強(qiáng)奇性積分顯然都位于系數(shù)矩陣的對(duì)角線上，因而利用物理學(xué)中的均勻場(chǎng)的概念，可以間接地對(duì)強(qiáng)奇性積分進(jìn)行求解。李小明[67]是基于直接邊界元理論的Helmholtz積分方程來(lái)求解聲輻射問題，其采用極坐標(biāo)法消除四邊形線性元中的弱奇異積分。趙鍵和汪鴻振等[86]在解決直接邊界元的Helmholtz積分方程的奇異積分方面，提出了在奇異點(diǎn)附近的積分區(qū)域，利用極坐標(biāo)變換消除積分的奇異性，將奇異積分變?yōu)楸环e函數(shù)無(wú)奇異點(diǎn)的普通積分來(lái)處理，并給出了表面Helmholtz方程、CHIEF及該文方法處理與不處理奇異積分的比較結(jié)果，得出經(jīng)該文介紹的辦法處理奇異積分后，各方法的計(jì)算誤差減小了一個(gè)數(shù)量級(jí)。利用極坐標(biāo)變化的方法來(lái)解決奇異積分問題，這種方法較退化單元法在處理弱奇異積分在積分變換方面更有效，但是其缺點(diǎn)同樣是不能處理高階奇異積分問題。

程昊等[58]對(duì)于直接邊界元法所涉及的奇異積分處理，采用替換點(diǎn)法，其基本思路是在奇異點(diǎn)周圍，不連續(xù)常單元內(nèi)選擇適當(dāng)?shù)膸讉€(gè)點(diǎn)，分別將這些點(diǎn)看作新的源點(diǎn)，通過(guò)采用高斯積分方法計(jì)算這些新源點(diǎn)與原奇異點(diǎn)，即場(chǎng)點(diǎn)之間的系數(shù)值，用這些選定的奇異點(diǎn)周圍的點(diǎn)代替奇異點(diǎn)計(jì)算積分值，然后將所有替換點(diǎn)所求得的值取平均值用以代替原奇異點(diǎn)處的值。這種方法在處理弱奇異積分時(shí)，可以在一定程度上減少計(jì)算量。

對(duì)于Cauchy奇異積分的處理，朱春浩和譚林森[95]對(duì)于直接邊界元法應(yīng)用于計(jì)算流固耦合問題時(shí)的弱奇異積分以及Cauchy奇異積分問題進(jìn)行討論，尤其是對(duì)于二次奇異積分問題的處理進(jìn)行了詳細(xì)的推導(dǎo)和證明，并以八節(jié)點(diǎn)曲邊單元為例討論了二次奇異積分的處理方式。孫亮和滕斌等[97]對(duì)于弱奇異積分進(jìn)行了研究，對(duì)源點(diǎn)附近的單元即存在弱奇異積分的單元進(jìn)行細(xì)分，從而改變計(jì)算區(qū)域上高斯積分點(diǎn)的分布。就是說(shuō)在弱奇異積分點(diǎn)附近進(jìn)行積分點(diǎn)加密，對(duì)于遠(yuǎn)離弱奇異積分點(diǎn)的積分使積分點(diǎn)分布稀疏一些。張秀珍[98]對(duì)于彈性問題的三維邊界積分方程中的弱奇異積分和Cauchy奇異積分進(jìn)行了討論，提出利用線性位移函數(shù)消除弱奇異因子；用剛體位移原理消除Cauchy奇異因子，避免了強(qiáng)奇異因子向弱奇異因子轉(zhuǎn)化計(jì)算效率低的問題。趙志高和黃其柏[99]提出了一種非等參單元的四邊形坐標(biāo)變換，其將積分的曲面單元映射為另一四邊形單元，通過(guò)兩次坐標(biāo)變換引入的雅可比行列式可以消除Helmholtz聲學(xué)邊界積分方程中的弱奇異型積分，利用?R/?n以及坐標(biāo)變換可以同時(shí)消除坐標(biāo)變換無(wú)法消除的Cauchy型奇異積分。李小瑜和傅志方[101]對(duì)邊界積分方程和邊界元技術(shù)計(jì)算結(jié)構(gòu)的穩(wěn)態(tài)外輻射聲場(chǎng)的方法進(jìn)行了研究，對(duì)邊界元方法所固有的奇異數(shù)值積分給出了處理方法，對(duì)于弱奇異積分其通過(guò)線性變換來(lái)進(jìn)行處理，對(duì)于Cauchy奇異積分的處理則是給出了利用離散后得到的線性方程組的特點(diǎn)來(lái)進(jìn)行處理。

對(duì)于超奇異積分的處理，閆再友等[89，96]在利用復(fù)合Helmholtz積分方程時(shí)，由于引入了表面Helmholtz積分方程的法向?qū)?shù)方程，所以導(dǎo)致了其積分算子中存在超奇異積分，其利用正則化關(guān)系式將超奇異積分算子降為弱奇異積分，然后對(duì)弱奇異積分進(jìn)行處理，從而解決超奇異積分計(jì)算的困難。其優(yōu)點(diǎn)在于可以直接快速地求出超奇異積分算子所對(duì)應(yīng)的系數(shù)矩陣。董春迎和姚振漢等[94]對(duì)邊界積分方程中存在的超奇異積分的數(shù)值解法作了介紹，其分別介紹了分部積分法、Stokes理論、散度理論、Guiggiani法和Kutt法等9種處理超奇異積分方法，文中雖然只針對(duì)于彈性問題的奇異積分方程的奇異積分進(jìn)行具體描述，但是對(duì)間接邊界元的Helmholtz積分方程的超強(qiáng)奇異積分處理是可以借鑒的。

在一般情況下，利用間接邊界元方法解Helmholtz積分方程，其超奇異積分的處理都是轉(zhuǎn)化成為雙層弱奇異積分來(lái)進(jìn)行的，只是這種雙層弱奇異積分的處理較直接邊界元涉及到的弱奇異積分的處理更為困難，而且這種奇異積分問題處理是否得當(dāng)是決定間接邊界元方法計(jì)算是否準(zhǔn)確的核心因素之一，國(guó)內(nèi)對(duì)于這類問題的研究非常少。

3.3 浸水結(jié)構(gòu)振動(dòng)及聲輻射特性工程問題的一些處理方法

對(duì)于浸水結(jié)構(gòu)在流固耦合作用下振動(dòng)及聲輻射特性的實(shí)際工程問題來(lái)說(shuō)，其研究的結(jié)構(gòu)往往是復(fù)雜的而且相對(duì)尺寸是非常龐大的，現(xiàn)階段研究人員一般較為常用的處理方法是直接利用有限元商業(yè)軟件進(jìn)行結(jié)構(gòu)和流體的建模，基于圣維南原理將無(wú)限域流體簡(jiǎn)化為有限域流體，利用其內(nèi)部的求解器進(jìn)行浸水結(jié)構(gòu)振動(dòng)特性的求解。

童宗鵬和王國(guó)治[48]以艦艇結(jié)構(gòu)振動(dòng)與噪聲預(yù)報(bào)為研究背景，采用有限元方法利用有限元軟件ANSYS對(duì)艦艇的干模型進(jìn)行研究，通過(guò)有限元模態(tài)計(jì)算和模態(tài)試驗(yàn)對(duì)理論計(jì)算模型進(jìn)行修正，確定艦艇干模型的參數(shù)；在該工作的基礎(chǔ)上，建立艦艇的流固耦合模型，模擬激振力施加在發(fā)動(dòng)機(jī)上，進(jìn)行流固耦合計(jì)算，得到艦艇外層殼體表面節(jié)點(diǎn)的位移。徐張明和華宏星等[102]利用有限元軟件ANSYS計(jì)算船舶結(jié)構(gòu)與水下聲場(chǎng)的耦合振動(dòng)，其在研究的過(guò)程中利用ANSYS軟件建立了水下船舶振動(dòng)的有限元模型，并計(jì)算了模型的濕模態(tài)以及在動(dòng)力系統(tǒng)中模擬的激勵(lì)源作用下的船舶殼板振動(dòng)響應(yīng)，其船體結(jié)構(gòu)和流體都是利用有限元單元來(lái)進(jìn)行模擬的。基于對(duì)浸水結(jié)構(gòu)振動(dòng)特性的計(jì)算，將計(jì)算得到的外殼面上的振動(dòng)速度作為聲場(chǎng)計(jì)算的邊界條件，利用邊界元軟件SYSNOISE計(jì)算聲學(xué)物理量。該文的方法為解決大型復(fù)雜結(jié)構(gòu)的耦合聲振預(yù)報(bào)提供了一個(gè)典型的實(shí)例。文獻(xiàn)[103-105]同樣利用有限元軟件ANSYS等結(jié)合邊界元軟件SYSNOISE等計(jì)算實(shí)際工程浸水結(jié)構(gòu)流固耦合振動(dòng)及聲輻射特性問題，這種方法需要分別對(duì)流體和結(jié)構(gòu)進(jìn)行有限元的單元?jiǎng)澐郑摲椒ㄒ延盟鹿こ探Y(jié)構(gòu)物以及水面艦艇的動(dòng)力特性計(jì)算，此方法優(yōu)點(diǎn)是可用于任意復(fù)雜實(shí)際工程浸水結(jié)構(gòu)振動(dòng)及聲輻射特性分析；缺點(diǎn)是由于其流體區(qū)域的有限元?jiǎng)澐执嬖诮財(cái)嗾`差，若需要得到較為精確的計(jì)算結(jié)果，往往單元?jiǎng)澐值脮?huì)比較多且比較精細(xì)，進(jìn)而計(jì)算量會(huì)非常大。而且在計(jì)算聲輻射特性時(shí)需要將有限元的計(jì)算結(jié)果導(dǎo)入邊界元軟件中，操作較繁瑣。

4 結(jié)語(yǔ)

對(duì)國(guó)內(nèi)相關(guān)文獻(xiàn)進(jìn)行分類總結(jié)后，關(guān)于浸水結(jié)構(gòu)振動(dòng)及聲輻射特性研究現(xiàn)狀有如下觀點(diǎn)：

（1）基于FEM/FEM的浸水結(jié)構(gòu)流固耦合振動(dòng)特性的研究，相對(duì)于其他研究方法在一些領(lǐng)域較為成熟，但是其缺點(diǎn)也比較明顯，即對(duì)于無(wú)限流體域的截?cái)嗾`差和對(duì)于三維流體域的計(jì)算量較大等；

（2）基于FEM/D-BEM的浸水結(jié)構(gòu)流固耦合振動(dòng)特性的研究發(fā)展相對(duì)較FEM/I-BEM在某些領(lǐng)域的研究發(fā)展要成熟一些，但其本身仍存在一些局限。如基于直接邊界元理論推導(dǎo)得出的附加質(zhì)量矩陣是非對(duì)稱的滿陣，在求解廣義特征方程的過(guò)程中存在一定的困難；不適用于開口結(jié)構(gòu)的研究等；

（3）基于FEM/I-BEM的浸水結(jié)構(gòu)流固耦合振動(dòng)特性的研究，優(yōu)點(diǎn)是基于間接邊界元推導(dǎo)得出的附加質(zhì)量陣是對(duì)稱的滿陣，在解決廣義特征方程的過(guò)程中比D-BEM更有優(yōu)勢(shì)，適用于開口結(jié)構(gòu)的研究等。但是由于其存在超奇異積分等特點(diǎn)，同樣給問題的解決帶來(lái)很大困難；

（4）對(duì)于實(shí)際工程中重流體中低頻振動(dòng)結(jié)構(gòu)附加質(zhì)量的計(jì)算，尤其是船舶的附加質(zhì)量的公式估算，本文雖沒過(guò)多涉及，但其無(wú)疑是科研和工程設(shè)計(jì)人員感興趣的。Korotkin[106]對(duì)船舶及相關(guān)構(gòu)件附加質(zhì)量系統(tǒng)地給出了適合工程實(shí)際的簡(jiǎn)便算法，包括各種船舶結(jié)構(gòu)在多種條件下（如固壁附近和自由表面附近等）的附加質(zhì)量計(jì)算公式，為設(shè)計(jì)提供了簡(jiǎn)單實(shí)用方法；

（5）對(duì)于浸水結(jié)構(gòu)的聲輻射問題，利用FEM/D-BEM方法，無(wú)論是在解決其特征頻率處解的非唯一性問題還是非超奇異積分的問題已經(jīng)有一些相對(duì)有效的處理辦法；但對(duì)I-BEM在解決某些問題時(shí)其積分方程涉及到超奇異積分問題，解決的方法相對(duì)于直接邊界元解決非超奇異積分的方法還是少一些；

（6）無(wú)論是實(shí)驗(yàn)驗(yàn)證還是試驗(yàn)研究，對(duì)于結(jié)構(gòu)流固耦合振動(dòng)特性問題都相對(duì)較少；對(duì)于結(jié)構(gòu)在流固耦合作用下的聲輻射方面的實(shí)驗(yàn)驗(yàn)證及試驗(yàn)研究，由于其操作起來(lái)更加困難，所以這方面的工作要更少一些。故對(duì)于這兩方面的研究還有待深入發(fā)展；

（7）雖然對(duì)結(jié)構(gòu)聲輻射靈敏度問題目前有些研究[58-60]，但是對(duì)于浸水結(jié)構(gòu)聲輻射靈敏度問題還是鮮有研究；

（8）對(duì)于本文所涉及的相關(guān)流固耦合振動(dòng)及聲輻射特性問題的文章，按照處理結(jié)構(gòu)和流體耦合時(shí)的三種不同的方法：FEM/FEM、FEM/D-BEM和FEM/I-BEM進(jìn)行分類。從圖5可以看出對(duì)于FEM/IBEM的研究相對(duì)較少；

（9）對(duì)于本文所涉及的相關(guān)流固耦合振動(dòng)及聲輻射特性問題的文章，按照處理聲輻射問題的兩種不同的方法：D-BEM和I-BEM進(jìn)行分類。從圖6可以看出對(duì)于I-BEM的研究相對(duì)較少；

（10）對(duì)于本文所涉及的相關(guān)流固耦合振動(dòng)及聲輻射特性問題的文章，按照數(shù)值模擬方法的兩種不同途徑分類，圖7表明文章基于計(jì)算機(jī)編程研究要多于基于商業(yè)軟件計(jì)算；

（11）對(duì)于本文所涉及的相關(guān)流固耦合振動(dòng)及聲輻射特性問題的文章，覆蓋的年代范圍是1985年到2011年，其年份分布如圖8所示。

圖5 文章分類：按照處理結(jié)構(gòu)和流體耦合問題時(shí)的三種方法Fig.5 Paper categories:According to the threemethods of processing fluid and structure coupling problem

圖6 文章分類：按照處理聲輻射問題的兩種方法Fig.6 Paper categories:According to the twomethods of processing acoustic radiation problem

圖7 文章分類：按照數(shù)值模擬方法的兩種不同途徑Fig.7 Paper categories:According to the two ways of numerical simulationmethod

圖8 文章分類：按年份分布Fig.8 Paper categories:According to the chronology

[1]Atalla N.Review of numerical solutions for low-frequency structural-acoustic problems[J].Applied Acoustics,1994,43 (3):271-294.

[2]汪庠寶,韓繼文,陸鑫森.懸臂板和連續(xù)板的流固耦合振動(dòng)解[J].中國(guó)造船,1988(2):53-63.

[3]汪庠寶,韓繼文,陸鑫森.船底板流固耦合振動(dòng)的有限元解[J].振動(dòng)與沖擊,1985(4):8-15.

[4]裴智勇,吳衛(wèi)國(guó),翁長(zhǎng)儉.高速船艙壁加筋板流固耦合振動(dòng)分析[J].工程力學(xué),2003(2):159-162.

[5]童予靖,劉正興.流固耦合問題中的附連水質(zhì)量研究[J].上海力學(xué),1997(4):311-320.

[6]金占禮,王宗利等.結(jié)構(gòu)在無(wú)限流體域中振動(dòng)時(shí)附連水質(zhì)量的數(shù)值計(jì)算方法[J].上海交通大學(xué)學(xué)報(bào),2000(8):1078-10820.

[7]王國(guó)強(qiáng),楊建民.用線性三角形元源匯法計(jì)算振動(dòng)船體的附連水質(zhì)量[J].上海交通大學(xué)學(xué)報(bào),1988(3):99-103.

[8]鄭治國(guó),王鎬章,趙德有.螺旋槳的流固耦合動(dòng)力分析[J].大連理工大學(xué)學(xué)報(bào),1996(2):219-223.

[9]錢勤,黃玉盈,劉忠族.求附連水質(zhì)量的一種直接方法[J].力學(xué)與實(shí)踐,1996(5):20-22.

[10]劉連海.含損傷加筋結(jié)構(gòu)水中振動(dòng)與聲輻射特性研究[D].大連:大連理工大學(xué),2006.

[11]劉連海,洪明,郭新毅.含損傷加筋結(jié)構(gòu)流固耦合動(dòng)力特性與聲輻射特性研究[J].中國(guó)艦船研究,2007(2):9-14.

[12]王其峰.流體加載下含損傷板式結(jié)構(gòu)聲振特性研究[D].大連:大連理工大學(xué),2009.

13]毛程亮.船舶板式結(jié)構(gòu)振動(dòng)及聲輻射特性研究[D].大連:大連理工大學(xué),2009.

[14]鄒元杰,趙德有.水下結(jié)構(gòu)聲固耦合振動(dòng)的特征值計(jì)算[J].船舶力學(xué),2004,8(2):109-120.

[15]付喜華,洪明,劉連海.艦船加筋結(jié)構(gòu)聲輻射數(shù)值模擬方法研究[Z].中國(guó)陜西西安:2007:28-34.

[16]付喜華.艦船加筋結(jié)構(gòu)聲輻射阻尼研究[D].大連:大連理工大學(xué),2008.

[17]付喜華,洪明,郭新毅.流體中彈性板聲輻射阻尼研究[J].中國(guó)艦船研究,2008(1):42-47.

[18]施衛(wèi)華.浸水結(jié)構(gòu)動(dòng)力特性與聲輻射特性研究[D].大連:大連理工大學(xué),2010.

[19]魏建輝,陳美霞,和衛(wèi)平等.基于空氣中模態(tài)參數(shù)的水中結(jié)構(gòu)固有頻率估算方法研究[C]//第十二屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集.中國(guó)長(zhǎng)沙,2009:115-121.

[20]安小同,洪明,李艮田.基于BEM/FEM計(jì)算浸水結(jié)構(gòu)振動(dòng)特性[J].船海工程,2011(6):55-58.

[21]蘇海東,黃玉盈.求半無(wú)限域流場(chǎng)中物體附連水質(zhì)量的一種簡(jiǎn)便解法[J].華中科技大學(xué)學(xué)報(bào)(城市科學(xué)版),2003(4): 14-16.

[22]李華東,朱錫,羅忠等.附連水質(zhì)量的邊界元法求解[J].海軍工程大學(xué)學(xué)報(bào),2009(2):45-49.

[23]張升明.流體的可壓縮性對(duì)彈性結(jié)構(gòu)振動(dòng)的影響[J].水動(dòng)力學(xué)研究與進(jìn)展(A輯),1994(4):429-436.

[24]毛程亮,洪明.流體可壓縮性對(duì)板及加筋板聲振特性數(shù)值模擬的影響[C]//第四屆全國(guó)船舶與海洋工程學(xué)術(shù)會(huì)議論文集.中國(guó)鎮(zhèn)江,2009:23-28.

[25]吳紹亮,金咸定.流固耦合計(jì)算方法在船舶局部結(jié)構(gòu)中的應(yīng)用[J].振動(dòng)與沖擊,2003(4):28-30.

[26]劉連海,洪明等.結(jié)構(gòu)流固耦合動(dòng)力分析中非對(duì)稱廣義特征值問題Lanczos算法[C]//2006中國(guó)大連國(guó)際海事論壇論文集.中國(guó)大連,2006.

[27]楊建民,汪庠寶,王國(guó)強(qiáng).船體結(jié)構(gòu)流固耦合模態(tài)分析[J].上海交通大學(xué)學(xué)報(bào),1993(1):19-26.

[28]吳芳.板和加筋板附加水質(zhì)量算法研究[D].大連:大連理工大學(xué),2006.

[29]王杰德,楊永謙,毛政良.船體薄壁梁彎扭動(dòng)力分析的流固耦合邊界元法[J].武漢交通科技大學(xué)學(xué)報(bào),1995(4):345-352.

[30]沈惠明,趙德有,羅志雍.流固耦合振動(dòng)問題的特征值解法[J].大連理工大學(xué)學(xué)報(bào),1990(3):369-371.

[31]朱曉芳,杜志鵬,金咸定.浸入水中的船舶尾軸架固有頻率的計(jì)算方法探討[J].振動(dòng)與沖擊,2005(2):111-115.

[32]姚熊亮,楊娜娜,陶景橋.雙層殼體水下振動(dòng)和聲輻射的仿真分析[J].哈爾濱工程大學(xué)學(xué)報(bào),2004(2):136-140.

[33]陳樂佳,駱東平,關(guān)珊珊等.潛艇雙層殼體振動(dòng)特性與聲輻射性能FEM/BEM分析[C]//第十屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集.中國(guó)煙臺(tái),2005:104-110.

[34]高菊,陳美霞,段曄鑫.基于FEM/BEM的環(huán)肋圓柱殼振動(dòng)及聲輻射相似性數(shù)值分析[C]//第十二屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集.中國(guó)長(zhǎng)沙,2009:44-51.

[35]尹緒超,姚熊亮,金葉青.軸向激勵(lì)下加肋圓柱殼體振動(dòng)及聲輻射特性研究[C]//第十三屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集.中國(guó)鷹潭,2011:166-171.

[36]楊吉新,張可,黨慧慧.基于ANSYS的流固耦合動(dòng)力分析方法[J].船海工程,2008(6):86-89.

[37]朱克強(qiáng),李道根,姚震球等.海洋浮式結(jié)構(gòu)流固耦合動(dòng)力建模分析[J].華東船舶工業(yè)學(xué)院學(xué)報(bào)(自然科學(xué)版),2002(1): 1-6.

[38]孫玉東,俞孟薩,李東升等.帶液艙水面船舶振動(dòng)和聲輻射的計(jì)算與分析[C]//第十屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集.中國(guó)煙臺(tái),2005:111-116.

[39]孫士麗,王詩(shī)平,姚熊亮等.板和加筋板流固耦合振動(dòng)特性研究[C]//黑龍江省造船工程學(xué)會(huì)2008年學(xué)術(shù)年會(huì)論文集.中國(guó)哈爾濱,2008:142-146.

[40]陳美霞,杜磊,陳樂佳等.基于邊界元法的平板結(jié)構(gòu)聲振性能數(shù)值計(jì)算[J].武漢理工大學(xué)學(xué)報(bào),2009(6):1048-1051.

[41]謝志勇,周其斗,紀(jì)剛.雙層柱殼的流固耦合模態(tài)計(jì)算與試驗(yàn)研究[J].海軍工程大學(xué)學(xué)報(bào),2009(2):97-101.

[42]徐張明,沈榮瀛,華宏星.利用FEM/IBEM計(jì)算流體介質(zhì)中的殼體的結(jié)構(gòu)聲耦合問題[J].振動(dòng)工程學(xué)報(bào),2002(3): 119-123.

[43]鄒元杰,趙德有.結(jié)構(gòu)在淺水中的振動(dòng)和聲輻射特性研究[J].振動(dòng)工程學(xué)報(bào),2004(3):21-26.

[44]岳鵬,黎勝,趙德有.船體薄壁梁彎扭耦合振動(dòng)的流固耦合分析[J].船舶力學(xué),2009(2):290-297.

[45]韓峰,胡迪科,閆桂榮.圓錐殼結(jié)構(gòu)聲振耦合特性分析[J].噪聲與振動(dòng)控制,2009(5):30-33.

[46]魏建輝,陳美霞,牟彬杰等.基于IBEM的水下雙層圓柱殼聲振性能研究[J].艦船科學(xué)技術(shù),2011(7):9-13.

[47]白長(zhǎng)青,周進(jìn)雄,閆桂榮.聲振耦合對(duì)薄壁圓柱結(jié)構(gòu)動(dòng)力特性的影響[J].機(jī)械工程學(xué)報(bào),2011(5):78-84.

[48]童宗鵬,王國(guó)治.艦艇結(jié)構(gòu)水下振動(dòng)和聲輻射特性研究[J].華東船舶工業(yè)學(xué)院學(xué)報(bào)(自然科學(xué)版),2003(2):18-22.

[49]胡選利,陳花玲,趙建平.聲振耦合數(shù)值模型合理性分析和實(shí)驗(yàn)驗(yàn)證[J].西安交通大學(xué)學(xué)報(bào),1997(7):3-8.

[50]姚熊亮,劉慶杰,翁強(qiáng)等.水下加筋圓柱殼體的振動(dòng)與近場(chǎng)聲輻射研究[J].中國(guó)艦船研究,2006(2):13-19.

[51]趙翔,黃幼玲,謝壯寧.自由場(chǎng)結(jié)構(gòu)體聲輻射研究[J].聲學(xué)學(xué)報(bào),1994(1):22-31.

[52]黎勝,趙德有.用有限元/邊界元方法進(jìn)行結(jié)構(gòu)聲輻射的模態(tài)分析[J].聲學(xué)學(xué)報(bào),2001(2):174-179.

[53]劉洪林.船舶結(jié)構(gòu)振動(dòng)與聲輻射分析[D].上海:上海交通大學(xué),2003.

[54]郭新毅,洪明,李艮田.含損傷加筋板結(jié)構(gòu)輻射聲功率及指向性變異研究[J].船舶力學(xué),2005,9(4):119-129.

[55]郭新毅,洪明,李艮田.含損傷加筋板結(jié)構(gòu)聲輻射模態(tài)變異研究[J].船舶力學(xué),2005,9(2):115-123.

[56]郭新毅.含損傷加筋板結(jié)構(gòu)振動(dòng)與聲特性變異研究[D].大連:大連理工大學(xué),2005.

[57]趙志高,黃其柏,何锃.基于有限元邊界元方法的薄板聲輻射分析[J].噪聲與振動(dòng)控制,2008(1):39-43.

[58]程昊,高煜,張永斌等.振動(dòng)體聲學(xué)靈敏度分析的邊界元法[J].機(jī)械工程學(xué)報(bào),2008(7):45-51.

[59]陳劍,程昊,高煜等.基于有限元-邊界元的聲學(xué)構(gòu)形靈敏度分析[J].振動(dòng)工程學(xué)報(bào),2009(2):213-217.

[60]張波,沈火明,支偉.有限元/邊界元法求解變厚度板的聲輻射[J].科學(xué)技術(shù)與工程,2010(36):8963-8967.

[61]柯兵,譚林森.用有限元/邊界元方法計(jì)算結(jié)構(gòu)體振動(dòng)輻射聲場(chǎng)[J].船海工程,2001(S2):97-99.

[62]陳美霞,邱昌林,駱東平.基于FEM/BEM法的內(nèi)部聲激勵(lì)水下圓柱殼聲輻射計(jì)算[J].中國(guó)艦船研究,2007(6):50-54.

[63]沈順根,李琪華,王大云等.加肋旋轉(zhuǎn)殼結(jié)構(gòu)噪聲聲輻射水彈性研究[J].中國(guó)造船,1992(2):55-64.

[64]俞孟薩,史小軍,陳克勤.采用有限元和邊界元方法分析彈性加肋圓柱殼的聲學(xué)相似性[J].中國(guó)造船,1999(3):65-71.

[65]張阿漫,錢德進(jìn),姚熊亮.結(jié)構(gòu)型式對(duì)雙層殼聲輻射特性影響研究[J].中國(guó)艦船研究,2007(3):1-6.

[66]石煥文,盛美萍,孫進(jìn)才等.加縱肋平底圓柱殼振動(dòng)和聲輻射的FEM/BEM研究[J].振動(dòng)與沖擊,2006(2):88-92.

[67]李小明.加肋圓柱殼振動(dòng)聲輻射特性的數(shù)值計(jì)算研究[D].大連:大連理工大學(xué),2008.

[68]袁自強(qiáng).圓柱殼的振動(dòng)聲輻射的數(shù)值分析[D].武漢:華中科技大學(xué),2008.

[69]楊勇.周期加強(qiáng)薄板聲輻射數(shù)值方法研究[D].大連:大連理工大學(xué),2009.

[70]金葉青,姚熊亮,尹緒超等.復(fù)雜錐柱結(jié)構(gòu)聲輻射數(shù)值計(jì)算研究[C]//第十三屆船舶水下噪聲學(xué)術(shù)討論會(huì)論文集.中國(guó)鷹潭,2011:127-133.

[71]洪明,郭新毅,劉連海.含損傷加筋結(jié)構(gòu)的振動(dòng)與聲輻射特性研究[J].中國(guó)艦船研究,2007(6):35-41.

[72]張升明,潘旭初.板架結(jié)構(gòu)的振動(dòng)噪聲研究[J].噪聲與振動(dòng)控制,1995(5):9-13.

[73]崔宏武,趙德有,羅志雍等.結(jié)構(gòu)振動(dòng)的水中聲輻射計(jì)算[J].中國(guó)造船,1990(4):51-56.

[74]黎勝,趙德有.流體加載下加肋板結(jié)構(gòu)的聲輻射特性研究[J].應(yīng)用聲學(xué),2000(6):28-32.

[75]黎勝,趙德有.用耦合有限元/邊界元方法研究加筋板的聲傳輸[J].振動(dòng)工程學(xué)報(bào),2001(3):120-123.

[76]黎勝,趙德有.水下結(jié)構(gòu)振動(dòng)和聲輻射的模態(tài)分析及控制研究[C]//水下噪聲學(xué)術(shù)論文選集(1985-2005).中國(guó)煙臺(tái),2005:45-56.

[77]楊德慶,鄭靖明,王德禹等.基于SYSNOISE軟件的船舶振動(dòng)聲學(xué)數(shù)值計(jì)算[J].中國(guó)造船,2002(4):32-37.

[78]楊德慶,王德禹,劉洪林.艦艇振動(dòng)聲學(xué)特性數(shù)值分析[J].上海交通大學(xué)學(xué)報(bào),2002(11):1537-1539.

[79]郝承智.水下加肋殼體振動(dòng)與聲輻射特性研究[D].西北工業(yè)大學(xué),2006.

[80]洪明,王其峰,劉連海等.輕/重流體介質(zhì)中振動(dòng)板聲輻射模態(tài)下的聲輻射阻尼研究[J].船舶力學(xué),2009,13(5): 813-821.

[81]郭新毅,洪明.結(jié)構(gòu)振動(dòng)聲輻射阻尼研究[C]//2004年船舶與海洋工程學(xué)術(shù)研討會(huì)論文集.2004:189-196.

[82]程武山.ANSYS軟件與SYSNOISE軟件結(jié)合的聲輻射分析[C]//全國(guó)先進(jìn)制造技術(shù)高層論壇暨第九屆制造業(yè)自動(dòng)化與信息化技術(shù)研討會(huì)論文集.中國(guó)武夷山,2010:299-302.

[83]楊婧媛.水下結(jié)構(gòu)聲輻射的響應(yīng)變異性分析及代理模型研究[D].大連:大連理工大學(xué),2010.

[84]施衛(wèi)華,鄧海華,陳明.浸水板式結(jié)構(gòu)的動(dòng)力特性與聲輻射特性數(shù)值分析[J].艦船科學(xué)技術(shù),2011(8):29-37.

[85]安小同.基于有限元/邊界元的水下結(jié)構(gòu)振動(dòng)聲輻射研究[D].大連:大連理工大學(xué),2011.

[86]趙鍵,汪鴻振,朱物華.邊界元法計(jì)算已知振速封閉面的聲輻射[J].聲學(xué)學(xué)報(bào),1989(4):250-257.

[87]趙鍵,汪鴻振,朱物華.改進(jìn)計(jì)算對(duì)稱殼體聲輻射的邊界積分方程法[J].上海交通大學(xué)學(xué)報(bào),1989(4):61-69.

[88]趙鍵,汪鴻振.用邊界積分方程法計(jì)算封閉體的聲輻射[C]//水下噪聲學(xué)術(shù)論文選集(1985-2005).中國(guó)煙臺(tái),2005: 4-17.

[89]閆再友,姜楫,嚴(yán)明.利用邊界元法計(jì)算無(wú)界聲場(chǎng)中結(jié)構(gòu)體聲輻射[J].上海交通大學(xué)學(xué)報(bào),2000(4):520-523.

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

[91]白楊,汪鴻振.聲學(xué)-結(jié)構(gòu)設(shè)計(jì)靈敏度分析[J].振動(dòng)與沖擊,2003(3):45-47.

[92]趙志高,黃其柏.復(fù)雜結(jié)構(gòu)的聲輻射解耦及其聲輻射效率分析[J].振動(dòng)工程學(xué)報(bào),2004(3):78-83.

[93]蔣偉,何正耀,王冬海.輻射體聲場(chǎng)計(jì)算的綜合Helmholtz積分公式法[J].電聲技術(shù),2006(10):9-12.

[94]董春迎,杜慶華,謝志成,姚振漢.邊界積分方程中超奇異積分的解法[J].力學(xué)進(jìn)展,1995(3):424-429.

[95]朱春浩等.用邊界元法計(jì)算流固耦合問題時(shí)高次奇異積分的處理方法[J].武漢船舶職業(yè)技術(shù)學(xué)院學(xué)報(bào),2002(4):49-52.

[96]閆再友,姜楫,何友聲等.聲學(xué)邊界元方法中超奇異數(shù)值積分處理的新方法[J].聲學(xué)學(xué)報(bào),2001(3):282-286.

[97]孫亮,滕斌,寧德志.處理準(zhǔn)奇異積分的自適應(yīng)高斯積分法[J].大連理工大學(xué)學(xué)報(bào),2007(1):106-112.

[98]張秀珍.三維邊界積分方程中的奇異積分的處理[J].裝甲兵工程學(xué)院學(xué)報(bào),2003(3):17-19.

[99]趙志高,黃其柏.Helmholtz聲學(xué)邊界積分方程中奇異積分的計(jì)算[J].工程數(shù)學(xué)學(xué)報(bào),2004(5):779-784.

[100]李宏偉,文立華,陳克安.非光滑結(jié)構(gòu)表面聲學(xué)邊界元法中奇異積分計(jì)算方法[J].西北工業(yè)大學(xué)學(xué)報(bào),2006(4): 497-500.

[101]李小瑜,傅志方.結(jié)構(gòu)振動(dòng)輻射聲場(chǎng)的預(yù)估-邊界積分方程中奇異積分的間接處理[J].振動(dòng)工程學(xué)報(bào),1989(1):59-65.

[102]徐張明,汪玉,華宏星等.船舶結(jié)構(gòu)的建模及水下振動(dòng)和輻射噪聲的FEM/BEM計(jì)算[J].船舶力學(xué),2002,6(4): 89-95.

[103]徐張明,汪玉,華宏星等.雙層殼體的船舶動(dòng)力艙振動(dòng)與聲輻射的有限元結(jié)合邊界元數(shù)值計(jì)算[J].中國(guó)造船, 2002(4):39-44.

[104]彭旭,駱東平.船舶結(jié)構(gòu)建模及水下振動(dòng)和輻射噪聲預(yù)報(bào)[J].噪聲與振動(dòng)控制,2003(6):9-12.

[105]鄒春平,陳端石,華宏星.船舶水下輻射噪聲特性研究[J].船舶力學(xué),2004,8(1):113-124.

[106]Korotkin A I.Added Masses of Ship Structures[M].Springer,2007.

Domestic review of the submerged structure vibration and acoustic radiation characteristics based on FEM/BEM

WANG Zhenga,HONGMinga,b,LIU Chenga
(a.School of Naval Architecture Engineering;b.State Key Laboratory of Structure Analysis for Industrial Equipment, Dalian University of Technology,Dalian 116024,China)

This paper summarizes the papers about the submerged structure vibration and acoustic radiation characteristics based on FEM or FEM/BEM published in China,emphasizes the numerical simulation on fluid,reviews and analyses the domestic research about the vibration characteristics of fluid-structure coupling,the acoustic radiation ofweak and strong coupling and the non-uniqueness and singularity problems of D-BEM and I-BEM,summarizes the research and the least development about the submerged structure vibration and acoustic radiation characteristics,and provides some reference to other researcherswho study the submerged structure vibration and acoustic radiation characteristics.

FEM;DBEM;IBEM;fluid-structure coupling;acoustic radiation;non-uniqueness;singularity

O326

A

10.3969/j.issn.1007-7294.2014.11.015

1007-7294（2014）11-1397-18

2014-06-05

國(guó)家自然科學(xué)基金資助項(xiàng)目（51079027）

王崢（1987-），男，碩士研究生，E-mail：zhwang@yahoo.cn；洪明（1959-），男，大連理工大學(xué)船舶工程學(xué)院教授，E-mail：mhong@dlut.edu.cn。