余 鵬，張曉宇，陳 杰，杜兆偉，殷 洪

(武漢第二船舶設(shè)計(jì)研究所，湖北 武漢 430205)

0 引 言

在水下航行器的初步設(shè)計(jì)階段，對(duì)其主要組成的圓柱殼結(jié)構(gòu)進(jìn)行聲振響應(yīng)預(yù)報(bào)，可以為結(jié)構(gòu)的聲學(xué)優(yōu)化設(shè)計(jì)提供理論基礎(chǔ)[1-2]。

在實(shí)際工程中，為了對(duì)結(jié)構(gòu)振動(dòng)數(shù)據(jù)進(jìn)行采樣，需要在結(jié)構(gòu)表面布置相關(guān)傳感器，但在整個(gè)結(jié)構(gòu)上均布置傳感器實(shí)施難度大、不現(xiàn)實(shí)。所以，通過一種好的測(cè)點(diǎn)布置方法，實(shí)現(xiàn)較少數(shù)目的傳感器對(duì)結(jié)構(gòu)聲振響應(yīng)的準(zhǔn)確預(yù)報(bào)就顯得非常重要。目前，國(guó)內(nèi)研究人員基本是采用均勻布置方法[3-7]，這種方法簡(jiǎn)單易行，但測(cè)點(diǎn)數(shù)目及位置的選擇存在較大的主觀性和盲目性，且需借助大量的數(shù)值分析算例。

幾十年來(lái)，學(xué)者們根據(jù)自然界中的生物進(jìn)化或物理演變過程，提出了很多基于概率變換來(lái)實(shí)現(xiàn)最優(yōu)組合的隨機(jī)類算法，如遺傳算法等[8]。這類算法可在全域內(nèi)進(jìn)行搜索，不易陷入局部最優(yōu)解，比較適合解決組合優(yōu)化類的問題，但也存在收斂速度慢、迭代次數(shù)多、計(jì)算效率有待改善等缺點(diǎn)。

針對(duì)圓柱殼結(jié)構(gòu)的測(cè)點(diǎn)布置問題，本文提出一種改進(jìn)方法—均勻子集隨機(jī)選取測(cè)點(diǎn)方法：將所有待選測(cè)點(diǎn)分成相互獨(dú)立的均勻子集，子集數(shù)量與測(cè)點(diǎn)數(shù)量保持一致，測(cè)點(diǎn)從各個(gè)子集中隨機(jī)選取，再利用設(shè)定的目標(biāo)函數(shù)對(duì)隨機(jī)選取的測(cè)點(diǎn)組進(jìn)行評(píng)價(jià)，從而選擇出好的測(cè)點(diǎn)布置方案。這種測(cè)點(diǎn)布置改進(jìn)方法不僅具有均勻布置方法中取樣全面、操作簡(jiǎn)單等優(yōu)點(diǎn)，同時(shí)通過劃分子集進(jìn)行隨機(jī)選擇，可避免選擇范圍廣、收斂速度慢、迭代次數(shù)多等缺點(diǎn)。

根據(jù)模態(tài)疊加原理，可將結(jié)構(gòu)振速展開為模態(tài)疊加的形式，再利用均勻子集隨機(jī)選取方法得到的測(cè)點(diǎn)布置方案，與均勻布置方法對(duì)比。通過上述不同測(cè)點(diǎn)布置方案得到的水下圓柱殼振動(dòng)聲輻射預(yù)報(bào)結(jié)果對(duì)比表明，均勻子集隨機(jī)選取測(cè)點(diǎn)方法是一種優(yōu)秀的測(cè)點(diǎn)布置方法。

1 基本理論

1.1 模態(tài)疊加理論

水下彈性結(jié)構(gòu)振動(dòng)響應(yīng)方程如下：

其中：M為質(zhì)量矩陣；C為阻尼矩陣；K為剛度矩陣；p(t)為 流體與結(jié)構(gòu)交界面上的流體載荷；f(t)為除了p(t)以 外的外界激勵(lì)；x(t)為結(jié)構(gòu)各節(jié)點(diǎn)的位移向量。

經(jīng)拉普拉斯變換可得：

將結(jié)構(gòu)振速展開為如下形式：

式中：V為結(jié)構(gòu)速度矢量； ?r為第r階模態(tài)振型；qr為第r階模態(tài)坐標(biāo)。

速度與位移的關(guān)系如下：

結(jié)合式（2）～式（4）可得：

式中：ms=?Ts M?s為模態(tài)質(zhì)量，cs=?TsC?s為模態(tài)阻尼，ks=?TsK?s為模態(tài)剛度。

由此，可得第S階模態(tài)坐標(biāo)為：

當(dāng)阻尼及外部激勵(lì)均為0 時(shí)，式（5）可表示如下：

式中： ωr為結(jié)構(gòu)在空氣中的固有頻率。對(duì)式（7）進(jìn)行質(zhì)量歸一，即ms=?sTM?s=1，可得：

根據(jù)式（7），當(dāng)結(jié)構(gòu)處于流體介質(zhì)中時(shí)，其對(duì)結(jié)構(gòu)振動(dòng)影響可表示成附加質(zhì)量的形式：

式中： ?ms為模態(tài)附連水質(zhì)量； ωs為結(jié)構(gòu)水中模態(tài)的固有頻率。

結(jié)合式（7）～式（9），可得：

即可通過式（10），對(duì)結(jié)構(gòu)空氣和水中的振動(dòng)固有頻率進(jìn)行轉(zhuǎn)換。

1.2 均勻子集隨機(jī)選取測(cè)點(diǎn)方法

均勻子集隨機(jī)選取測(cè)點(diǎn)包含以下3 個(gè)方面內(nèi)容：

1）子集劃分

將所有待選測(cè)點(diǎn)劃分成相互獨(dú)立的均勻子集，子集數(shù)量與測(cè)點(diǎn)數(shù)量保持一致。并且可根據(jù)結(jié)構(gòu)形狀選擇合適的子集劃分方式，針對(duì)本文中的單層圓柱殼結(jié)構(gòu)，可沿軸向和周向進(jìn)行劃分，圖1 所示將圓柱殼沿軸向分成了4 個(gè)子集，圖2 所示將圓柱殼沿周向分成了8 個(gè)子集，子集的表示方式為AiBj(i=1,···,4;j=1,···,8)， 子集A1和A2交 界處的節(jié)點(diǎn)視為子集A1中的節(jié)點(diǎn)，子集B1和B2交 界處的節(jié)點(diǎn)視為子集B2中的節(jié)點(diǎn)。

圖1 沿軸向劃分子集示意圖Fig. 1 The axial subset

2）測(cè)點(diǎn)選取

從每個(gè)子集中隨機(jī)選擇一個(gè)測(cè)點(diǎn)形成測(cè)點(diǎn)組合，并且引入循環(huán)策略，循環(huán)選擇的次數(shù)應(yīng)足夠多，以保證能夠大概率的選到子集中每一個(gè)待選測(cè)點(diǎn)。

3）目標(biāo)函數(shù)設(shè)定

圖2 沿周向劃分子集示意圖Fig. 2 The circumferential subset

由于每個(gè)測(cè)點(diǎn)均是從對(duì)應(yīng)子集中隨機(jī)選取，所以通過設(shè)定目標(biāo)函數(shù)，對(duì)選取的測(cè)點(diǎn)組合進(jìn)行評(píng)判?？筛鶕?jù)實(shí)際情況設(shè)定合適的目標(biāo)函數(shù)，本文將測(cè)點(diǎn)組模態(tài)矩陣的MAC 值[9-10]作為目標(biāo)函數(shù)，如下式：

式 中： ?i為第i階模態(tài)向量； ?j為第j階模態(tài)向量；MAC值為模態(tài)向量的余弦值，其值越小，模態(tài)向量的空間夾角就越大，對(duì)應(yīng)振型越易識(shí)別。

本文通過計(jì)算不同測(cè)點(diǎn)組的模態(tài)矩陣MAC值，從而選擇MAC值最小的測(cè)點(diǎn)組，再對(duì)水下單層圓柱殼結(jié)構(gòu)的振動(dòng)聲輻射進(jìn)行預(yù)報(bào)，以驗(yàn)證本文方法的準(zhǔn)確性。

1.3 振動(dòng)聲輻射預(yù)報(bào)基本理論

假設(shè)模態(tài)階數(shù)為N，節(jié)點(diǎn)數(shù)為n，以矩陣形式將振速在模態(tài)空間中展開如下：

根據(jù)式（12），若測(cè)點(diǎn)數(shù)為M個(gè)，則測(cè)點(diǎn)振速可表示如下：

已知測(cè)點(diǎn)振速 [V]M×1， 和前N階模態(tài)向量 [Φ]M×N，根據(jù)式（13）可求解得到模態(tài)坐標(biāo)，進(jìn)而通過式（12）可得到結(jié)構(gòu)表面所有節(jié)點(diǎn)振速 [V]n×1，再結(jié)合邊界元法得到的聲傳遞矢量ATV，即可求解輻射聲壓如下：

2 數(shù)值分析

2.1 有限元建模

本文中的單層圓柱殼結(jié)構(gòu)相關(guān)參數(shù)如下：

總 長(zhǎng)L=1.2 m， 半 徑R=0.4 m，殼 體 厚t=0.01 m ，密 度 ρ=7 800 kg/m3，彈 性 模 量E=210 GPa，泊 松 比μ=0.3，阻尼比ξ =0.005，尺寸如圖3 所示。

圖3 結(jié)構(gòu)尺寸示意圖Fig. 3 The structure size

對(duì)結(jié)構(gòu)兩端面上的節(jié)點(diǎn)約束橫向和垂向位移，在結(jié)構(gòu)中間正下方作用一垂向單位力。同時(shí)在距離結(jié)構(gòu)100 m 處，與激勵(lì)力相同方向和垂直方向分別布置1 個(gè)聲壓場(chǎng)點(diǎn)，如圖4 所示。

圖4 聲壓場(chǎng)點(diǎn)位置示意圖Fig. 4 The sound pressure field point

采用Ansys 軟件對(duì)結(jié)構(gòu)進(jìn)行有限元建模，在軸向和周向分別劃分24 和48 個(gè)網(wǎng)格，頻率分析范圍由2～1000 Hz，間隔2 Hz。有限元模型如圖5 所示。

圖5 結(jié)構(gòu)有限元模型Fig. 5 The structural finite element model

2.2 模態(tài)分析

根據(jù)上述理論分析，結(jié)合有限元軟件，可得到結(jié)構(gòu)的模態(tài)分析結(jié)果如表1 所示。

根據(jù)上表分析結(jié)果，第47 階空氣振動(dòng)模態(tài)對(duì)應(yīng)的水中頻率已達(dá)1000 Hz 以上。為了避免頻率范圍內(nèi)的模態(tài)遺漏，本文選取了空氣中前100 階模態(tài)進(jìn)行分析。

表1 固有頻率換算表Tab. 1 Natural frequency conversion table

2.3 振動(dòng)聲輻射預(yù)報(bào)

為驗(yàn)證本文方法的準(zhǔn)確性，開展對(duì)比分析。

1）均勻測(cè)點(diǎn)布置

目前，均勻測(cè)點(diǎn)布置方法在圓柱殼結(jié)構(gòu)振動(dòng)采樣中應(yīng)用最為廣泛，所以，本文先針對(duì)這種方法預(yù)報(bào)結(jié)構(gòu)振動(dòng)聲輻射開展分析研究。

根據(jù)上述建模情況，模型總的節(jié)點(diǎn)數(shù)目為1200，去除邊界節(jié)點(diǎn)，其余均可作為待選測(cè)點(diǎn)，共計(jì)1104 個(gè)待選測(cè)點(diǎn)。按照測(cè)點(diǎn)布置的疏密程度不同，本文設(shè)計(jì)了以下3 種方案：方案1 的測(cè)點(diǎn)總數(shù)為264，其中軸向11 排，周向24 列，測(cè)點(diǎn)布置較為密集；方案2 相比方案1，將軸向測(cè)點(diǎn)減為6 排，周向測(cè)點(diǎn)不變，測(cè)點(diǎn)總數(shù)144；方案3 相比方案2，軸向測(cè)點(diǎn)不變，將周向測(cè)點(diǎn)減為12 列，測(cè)點(diǎn)總數(shù)72。

表2 三種均勻布置方案對(duì)比Tab. 2 Comparison of three uniform layout schemes

圖6 均勻布置方案1Fig. 6 The first uniform layout scheme

圖7 均勻布置方案2Fig. 7 The second uniform layout scheme

圖8 均勻布置方案3Fig. 8 The third uniform layout scheme

根據(jù)截?cái)嗟那?00 階空氣中模態(tài)、及上述3 種均勻測(cè)點(diǎn)布置方案，通過式（13）對(duì)結(jié)構(gòu)表面振速場(chǎng)進(jìn)行重構(gòu)。其中方案1 和方案2 的測(cè)點(diǎn)數(shù)目M大于模態(tài)數(shù)目N，采用最小二乘法求解；方案3 的測(cè)點(diǎn)數(shù)目M小于模態(tài)數(shù)目N，采用欠定分離法解決解的非唯一性問題。根據(jù)重構(gòu)的結(jié)構(gòu)表面振速，結(jié)合邊界元法求得的聲傳遞矢量ATV，通過式（14）便可得到各個(gè)場(chǎng)點(diǎn)的輻射聲壓，如圖9～圖11 所示。

圖9 均方振速對(duì)比結(jié)果Fig. 9 The change curve of square vibration velocity

圖10 1 號(hào)場(chǎng)點(diǎn)輻射聲壓對(duì)比結(jié)果Fig. 10 The radiated sound pressure of field point No.1

通過對(duì)比上述結(jié)果可知，在分析頻段內(nèi)，方案1和方案2 的預(yù)報(bào)結(jié)果精度更高，且方案2 的預(yù)報(bào)結(jié)果已基本達(dá)到收斂，但方案3 的預(yù)報(bào)結(jié)果精度明顯降低。

2）均勻子集隨機(jī)選取測(cè)點(diǎn)布置

圖11 2 號(hào)場(chǎng)點(diǎn)輻射聲壓對(duì)比結(jié)果Fig. 11 The radiated sound pressure of field point No.2

根據(jù)上述理論分析，該方法將所有待選測(cè)點(diǎn)分成相互獨(dú)立的均勻子集，子集數(shù)量與測(cè)點(diǎn)數(shù)量保持一致，從每個(gè)子集中隨機(jī)選擇一個(gè)測(cè)點(diǎn)形成測(cè)點(diǎn)組合，再通過目標(biāo)函數(shù)對(duì)每種測(cè)點(diǎn)組合進(jìn)行評(píng)價(jià)，從而選擇出好的測(cè)點(diǎn)布置方案。

以上述方案3 為參照，對(duì)應(yīng)其測(cè)點(diǎn)布置方案，本文將結(jié)構(gòu)沿周向劃分了12 個(gè)子集，軸向劃分了6 個(gè)子集，劃分的子集總數(shù)為72 個(gè)，如圖12 和圖13 所示。在循環(huán)策略中，為了盡可能選到每一個(gè)待選測(cè)點(diǎn)，本文將循環(huán)次數(shù)定為單個(gè)子集中待選測(cè)點(diǎn)數(shù)目的10 倍，考慮單個(gè)子集中最多有16 個(gè)待選測(cè)點(diǎn)，相應(yīng)需進(jìn)行160 次循環(huán)操作。再根據(jù)式（11）中設(shè)定的目標(biāo)函數(shù)，當(dāng)目標(biāo)函數(shù)MAC 值為0.1461 時(shí)達(dá)到最小，其對(duì)應(yīng)的測(cè)點(diǎn)布置方案如圖12 和圖13 所示。利用該組測(cè)點(diǎn)對(duì)結(jié)構(gòu)振動(dòng)聲輻射進(jìn)行預(yù)報(bào)，得到的預(yù)報(bào)結(jié)果與均勻布置方法對(duì)比曲線如圖14～圖16 所示。

圖12 軸向測(cè)點(diǎn)分布示意圖Fig. 12 The axial measuring point distribution

圖13 周向測(cè)點(diǎn)分布示意圖Fig. 13 The circumferential measuring point distribution

圖14 均方振速對(duì)比結(jié)果Fig. 14 The change curve of square vibration velocity

圖15 1 號(hào)場(chǎng)點(diǎn)輻射聲壓對(duì)比結(jié)果Fig. 15 The radiated sound pressure of field point No.1

圖16 2 號(hào)場(chǎng)點(diǎn)輻射聲壓對(duì)比結(jié)果Fig. 16 The radiated sound pressure of field point No.2

由圖12 和圖13 可知，本文通過劃分子集的方式選取測(cè)點(diǎn)，使得測(cè)點(diǎn)能夠覆蓋曲面且分布相對(duì)均勻。通過設(shè)定目標(biāo)函數(shù)MAC 值，可保證能夠選取到模態(tài)向量空間夾角小的測(cè)點(diǎn)組合，對(duì)應(yīng)的模態(tài)振型更易辨識(shí)。最后通過圖14～圖16 中的振動(dòng)聲輻射預(yù)報(bào)結(jié)果對(duì)比曲線可知，當(dāng)其他因素（如分析頻段、測(cè)點(diǎn)數(shù)目、模態(tài)截?cái)鄶?shù)目等）相同時(shí)，相比均勻測(cè)點(diǎn)布置方法，均勻子集隨機(jī)選取測(cè)點(diǎn)方法能夠更準(zhǔn)確地進(jìn)行振動(dòng)聲輻射預(yù)報(bào)，且能與參考結(jié)果吻合較好。

3 結(jié) 語(yǔ)

通過上述分析，可得到以下結(jié)論：

1）均勻子集隨機(jī)選取測(cè)點(diǎn)方法測(cè)點(diǎn)覆蓋曲面、分布均勻，并且測(cè)點(diǎn)選取不存在主觀性和盲目性。

2）本文以測(cè)點(diǎn)組模態(tài)矩陣的MAC 值作為目標(biāo)函數(shù)，可保證測(cè)點(diǎn)組的模態(tài)向量間的空間夾角較小，其模態(tài)振型相應(yīng)較易識(shí)別。

3）通過數(shù)值分析，當(dāng)分析頻段、測(cè)點(diǎn)數(shù)目、模態(tài)截?cái)鄶?shù)目等相同時(shí)，該方法對(duì)結(jié)構(gòu)振動(dòng)聲輻射預(yù)報(bào)結(jié)果的精度比均勻布置方法更高。