楊曉剛，林立，白振國，李兵

（1中國艦船研究院，北京100192；2中國船舶科學(xué)研究中心，江蘇無錫214082）

徑向激勵作用下有限長雙層圓柱殼聲振傳遞特性研究

楊曉剛1，林立1，白振國2，李兵1

（1中國艦船研究院，北京100192；2中國船舶科學(xué)研究中心，江蘇無錫214082）

文章利用Donnel殼體方程描述雙層圓柱殼振動，以附加阻抗的形式表示環(huán)肋、實肋板以及水介質(zhì)等對殼體的力作用，結(jié)合內(nèi)外殼振動位移模態(tài)展開形式，建立了雙層圓柱殼內(nèi)外殼之間的聲振傳遞矩陣。同時開展了雙層圓柱殼外殼均方速度響應(yīng)仿真分析，驗證了雙層圓柱殼內(nèi)外殼振動響應(yīng)之間存在固定的傳遞關(guān)系，不受激勵源特性影響，為雙層圓柱殼輻射噪聲預(yù)估提供了一定的理論指導(dǎo)。

聲振傳遞；傳遞矩陣；雙層圓柱殼

0 引言

機械噪聲是水下航行器低速工況下的主要噪聲源，因此機械設(shè)備激勵引起的結(jié)構(gòu)振動和聲輻射問題一直備受關(guān)注。水下航行器在航行過程中，機械設(shè)備通過隔振系統(tǒng)激勵基座，并沿基座傳遞到內(nèi)殼上，或者通過管路系統(tǒng)激起內(nèi)殼振動；然后部分振動能量經(jīng)雙層殼體之間的水介質(zhì)以及實肋板等支撐結(jié)構(gòu)傳遞到外殼體，并向外場輻射噪聲。

目前針對水下航行器輻射噪聲預(yù)報，中低頻段主要采用有限元法和邊界元法，中高頻段主要采用統(tǒng)計能量法，這些方法主要以設(shè)備激勵力或輸入功率為初始輸入?yún)?shù)，但是目前缺乏客觀合理的設(shè)備激勵力或其它有效的激勵輸入?yún)?shù)等[1]，難以由設(shè)備激勵或基座響應(yīng)準(zhǔn)確預(yù)報輻射噪聲。同時，由于外場水介質(zhì)作用，外殼表面難以布置大量振動傳感器[2]直接由外殼振動數(shù)據(jù)預(yù)報輻射噪聲。因此基于振動能量傳遞過程，利用內(nèi)殼振動數(shù)據(jù)預(yù)估外殼振動響應(yīng)，進而預(yù)估輻射噪聲成為比較有效的途徑，雙層圓柱殼聲振傳遞特性研究對輻射噪聲預(yù)估具有重要的意義。

國內(nèi)外學(xué)者針對雙層圓柱殼聲振傳遞特性開展了大量研究，主要是通過解析法和數(shù)值法建立雙層圓柱殼聲振耦合模型，獲取內(nèi)外殼聲振傳遞關(guān)系。白振國[3]利用模態(tài)展開法，建立了有限長雙層加肋圓柱殼聲振耦合模型，詳細(xì)分析了實肋板、環(huán)肋以及激勵特性等參數(shù)對雙層殼間聲振耦合特性的影響；陳美霞等[4-6]利用Flügge殼體理論，把環(huán)肋、縱骨、托板或?qū)嵗甙逡暈閯臃戳ψ饔迷跉んw上，建立了雙殼體結(jié)構(gòu)—聲—流體耦合方程，對比了連接形式等對雙層圓柱殼聲振特性的影響；金廣文[2]基于頻響函數(shù)求逆法，建立了內(nèi)外殼體之間的振動傳遞矩陣，該傳遞矩陣受激勵源、邊界條件等影響；文獻(xiàn)[7]利用有限元法和邊界元法建立了雙層殼體流固耦合數(shù)值模型，并對雙層圓柱殼的水下受激振動和聲輻射進行了數(shù)值分析。

考慮到實際水下航行器的復(fù)雜性，難以建立解析和數(shù)值簡化模型，同時存在計算量過大和精度不理想的問題。因此根據(jù)工程實際，可以利用模型或?qū)嶋H結(jié)構(gòu)的振動和噪聲測試數(shù)據(jù)，回歸獲取內(nèi)外殼體的聲振耦合傳遞關(guān)系，此技術(shù)途徑存在的一個關(guān)鍵性問題是首先要驗證內(nèi)外殼振動傳遞和聲場耦合是否存在與激勵源無關(guān)的傳遞關(guān)系，這樣才能以此為基礎(chǔ)，針對不同激勵力工況，由內(nèi)殼振動預(yù)估外殼振動，進一步預(yù)估輻射噪聲。為此，本文以有限長雙層圓柱殼為例，理論推導(dǎo)建立了內(nèi)殼振動到外殼振動的聲振傳遞矩陣，并數(shù)值仿真分析了激勵源特性對內(nèi)外殼聲振傳遞特性的影響。

1 雙層圓柱殼聲振傳遞理論

本文研究的模型為有限長雙層圓柱殼：內(nèi)殼體表面環(huán)向分布周期性環(huán)肋；內(nèi)外殼之間由實肋板連接并且充滿水介質(zhì)；外殼外部為無限大理想水介質(zhì)。雙層圓柱殼物理模型如圖1所示。且假設(shè)雙層圓柱殼兩端簡支，并有無限大圓柱形聲障板；殼間實肋板為聲透明，不影響殼間聲傳播。

雙層圓柱殼振動采用Donnel殼體理論描述。

圖1 有限長雙層圓柱殼模型Fig.1 Model of double finite length cylindrical shells

其中：ui、vi、wi分別為內(nèi)外殼的軸向、周向和徑向位移；ai、hi、分別為內(nèi)外殼半徑、厚度以及厚徑比；cp為殼體縱波聲速；E、μ分別為殼體材料楊氏模量和泊松比分別為環(huán)肋對內(nèi)外殼作用力分別為殼間實肋板對內(nèi)外殼作用力；分別為殼間以及外場水介質(zhì)對內(nèi)外殼的作用力為作用在圓柱殼上的激勵力。

僅考慮徑向?qū)ΨQ激勵，則雙層圓柱殼軸向、周向和徑向位移可表示為：

將方程式（2）代入到雙層圓柱殼振動方程（1）中，可獲得雙層圓柱殼振動矩陣方程，同時將環(huán)肋、實肋板等結(jié)構(gòu)對內(nèi)外殼體作用視為離散的反作用力，結(jié)合邊界條件，推導(dǎo)殼體附加阻抗的表達(dá)式，進而組裝形成雙層圓柱殼聲振耦合方程組（3）式。

其中雙層圓柱殼模態(tài)阻抗矩陣中子矩陣元素Zij以及環(huán)肋、實肋板、水介質(zhì)等的附加作用阻抗詳細(xì)求解過程參見文獻(xiàn)[3]，F(xiàn)為徑向模態(tài)作用力。

水下航行器上機械設(shè)備等產(chǎn)生的激勵力主要為徑向作用力，且殼體聲輻射也主要由徑向振動激發(fā)，因此本文僅考慮雙層圓柱殼內(nèi)外殼體徑向聲振耦合特性，為了便于分析，移動（3）式模態(tài)阻抗矩陣中Z31行元素至Z51行元素之后，同時相應(yīng)調(diào)整（3）式中模態(tài)分量矩陣和模態(tài)作用力矩陣中元素的位置，將內(nèi)外殼體的軸向模態(tài)分量U1、U2和周向模態(tài)分量V1、V2組合為一個新變量M，得到雙層圓柱殼模態(tài)運動方程，如（4）式所示。

為了研究內(nèi)殼和外殼徑向聲振傳遞特性，需要建立內(nèi)殼的徑向位移模態(tài)分量W1和外殼的徑向位移模態(tài)分量W2之間關(guān)系，將方程組（4）分解，可得：

由（5）式可得：

進而將（6）式代入（2）式，利用內(nèi)殼的徑向模態(tài)分量W1表示外殼的徑向振動位移w2，可得：

將（8）式代入（7）式中，可以獲得內(nèi)外殼體聲振傳遞關(guān)系：

根據(jù)內(nèi)外殼聲振傳遞矩陣表達(dá)式可知，聲振傳遞矩陣僅與雙層圓柱殼模態(tài)阻抗、內(nèi)外殼振動測點位置及數(shù)量相關(guān)，與激勵源特性無關(guān)，而模態(tài)阻抗由雙層圓柱殼結(jié)構(gòu)決定，即雙層圓柱殼結(jié)構(gòu)和振動測點位置固定時，內(nèi)外殼之間存在固定的聲振傳遞矩陣。后文對雙層圓柱殼內(nèi)外殼振動響應(yīng)進行了數(shù)值仿真，分別從激勵源數(shù)量、位置和形式三個方面研究激勵源特性對內(nèi)外殼聲振傳遞關(guān)系影響，以驗證內(nèi)外殼聲振傳遞關(guān)系與激勵源特性無關(guān)。

2 仿真計算和分析

本節(jié)選用的有限長雙層圓柱殼模型參數(shù)：內(nèi)殼半徑a1=3.5 m，外殼半徑a2=4.3 m，內(nèi)殼厚度h1=28 mm，外殼厚度h2=10 mm，殼體長度L=9.6 m，內(nèi)殼上均布4個T型環(huán)肋，內(nèi)外殼之間由3個實肋板連接，實肋板厚度5 mm，殼體材料為鋼，楊氏模量E=2.1×1011N/m2，泊松比μ=0.3，內(nèi)外殼之間以及外殼外部充滿理想水介質(zhì)。

經(jīng)雙層圓柱殼聲振耦合理論收斂性驗證，周向波數(shù)N=41，軸向波數(shù)M=60時，雙層圓柱殼輻射噪聲有效計算頻段可達(dá)1-500 Hz。由雙層圓柱殼內(nèi)外殼聲振傳遞矩陣?yán)碚撏茖?dǎo)過程可知，內(nèi)殼應(yīng)布置M×N個振動測點，因此本文內(nèi)殼上周向均布41個測點，軸向均布60環(huán)測點，同時為有效反映外殼振動狀態(tài)，外殼軸向和周向各均布11個測點，從而建立雙層圓柱殼內(nèi)外殼聲振傳遞矩陣。

本文根據(jù)（3）式可求解出雙層圓柱殼內(nèi)外殼徑向振動位移模態(tài)分量，進而由（2）式可直接計算出內(nèi)外殼振動響應(yīng)；同時利用計算得到的內(nèi)殼振動響應(yīng)和內(nèi)外殼聲振傳遞矩陣，由（13）式可間接求解外殼振動響應(yīng)，分別對比激勵源數(shù)量、位置、形式發(fā)生變化時外殼均方速度響應(yīng)間接求解和直接求解結(jié)果，分析激勵源特性對內(nèi)外殼聲振傳遞關(guān)系的影響。振動速度級基準(zhǔn)為：v0=1×10-6m/s。

圖2 內(nèi)外殼均方速度級Fig.2 Comparison of mean square velocity level both the inner and outer shell

2.1 激勵數(shù)量對雙層圓柱殼聲振傳遞關(guān)系影響

圖3和圖4分別給出了單位點激勵力作用下外殼均方速度響應(yīng)和外殼測點1振動響應(yīng)。激勵作用位置1：外殼測點1位置：xw1=4.7 m，φw1=0。圖5為雙層圓柱殼軸向均布10個、20個單位激勵作用下外殼均方速度響應(yīng)。

圖3和圖4可知，單位激勵力作用下，外殼速度響應(yīng)的直接求解結(jié)果和間接求解結(jié)果完全吻合，即利用內(nèi)外殼聲振傳遞關(guān)系可準(zhǔn)確預(yù)估點激勵作用下外殼響應(yīng)。

由圖5可知，增加激勵源數(shù)量，雙層圓柱殼外殼均方速度響應(yīng)變大，即雙層圓柱殼振動特性發(fā)生變化，但利用內(nèi)殼振動響應(yīng)和聲振傳遞矩陣間接求解的外殼振動響應(yīng)與直接求解結(jié)果完全吻合，這表明雙層圓柱殼聲振傳遞關(guān)系不受激勵源數(shù)量影響。

圖3 單位激勵作用下外殼均方速度級Fig.3 Mean square velocity level of the outer shell under unite force excitation

圖4 單位激勵作用下測點1速度級Fig.4 The velocity level of measured point 1 under unite force excitation

圖5 不同數(shù)量激勵作用下外殼均方速度級Fig.5 Mean square velocity level of the outer shell under different amount of force excitations

2.2 激勵位置對雙層圓柱殼聲振傳遞關(guān)系影響

圖6給出了點激勵作用于不同位置處外殼均方速度響應(yīng)。激勵作用位置2：激勵位置3：兩激勵位置關(guān)于（2.1）中激勵作用位置1對稱。

由圖6可知，在1-500 Hz頻段內(nèi)，利用文中直接求解方法獲得的外殼平均振動響應(yīng)和利用雙層圓柱殼聲振傳遞矩陣間接獲得的外殼平均振動響應(yīng)曲線完全吻合，這表明激勵源位置不會對雙層圓柱殼內(nèi)外殼聲振傳遞關(guān)系產(chǎn)生影響。

圖6 激勵作用于不同位置時外殼均方速度級Fig.6 Mean square velocity level of the outer shell under force excitation to different position

2.3激勵形式對雙層圓柱殼聲振傳遞關(guān)系影響

圖7為周向線激勵力、軸向線激勵力作用下外殼均方速度響應(yīng)。軸向和周向線激勵力都以（2.1）中激勵作用位置1為中心，作用長度為2 m，各個力與單位點激勵力在積分意義上相等。

由圖7可知，將點激勵更換為線激勵，激勵源形式發(fā)生變化，但外殼平均速度響應(yīng)的直接求解結(jié)果和間接求解結(jié)果仍完全吻合，即改變激勵力的形式，雙層圓柱殼內(nèi)外殼之間依然存在固定的聲振傳遞關(guān)系。

圖8為長方形面激勵、環(huán)形面激勵作用下外殼均方速度響應(yīng)。兩種形式面力都以為中心，長方形面激勵軸向長度為9.6 m，周向?qū)挾葹?.2 m；環(huán)形面激勵寬度為0.2 m，兩個面力與單位點激勵力在積分意義上相等。

圖7 線力作用下外殼均方速度級Fig.7 Mean square velocity level of the outer shell under different line force excitation

圖8 面力作用下外殼均方速度級Fig.8 Mean square velocity level of the outer shell under different surface force excitation

由圖8可知，兩種不同形式面激勵作用下，外殼振動響應(yīng)發(fā)生變化，尤其是200 Hz以下中低頻段，環(huán)形面力作用下，外殼振動響應(yīng)小于長方形面力作用下外殼振動響應(yīng)，但是利用內(nèi)殼響應(yīng)和內(nèi)外殼聲振傳遞矩陣依然能夠準(zhǔn)確預(yù)估外殼響應(yīng)，進一步表明文中獲得的聲振傳遞矩陣與激勵源形式無關(guān)。

2.4 內(nèi)殼測點位置偏差對外殼響應(yīng)預(yù)估的影響

考慮到實際水下航行器結(jié)構(gòu)復(fù)雜，難以獲得內(nèi)殼振動傳感器的精確位置，從而導(dǎo)致實際傳感器布置位置與聲振傳遞矩陣求解過程使用的內(nèi)殼測點位置有一定的偏差。因此，本節(jié)在前文研究基礎(chǔ)上，在一定范圍內(nèi)改變內(nèi)殼測點位置，利用單位激勵下，測點位置變化后內(nèi)殼振動數(shù)據(jù)和前文建立的聲振傳遞矩陣預(yù)估外殼振動響應(yīng)，并與利用測點位置未發(fā)生變化時內(nèi)殼振動數(shù)據(jù)預(yù)估的外殼響應(yīng)結(jié)果對比，從而分析內(nèi)殼測點位置偏差對外殼響應(yīng)預(yù)估精度的影響。圖9為單位激勵下，內(nèi)殼測點位置在軸向、周向以及在軸向和周向都發(fā)生偏差前后利用內(nèi)外殼聲振傳遞矩陣間接求解的外殼均方速度響應(yīng)對比曲線，測點位置軸向和周向的偏差范圍為5 cm。

由圖9可知，內(nèi)殼測點位置在軸向發(fā)生偏差時，40 Hz以下，150 Hz以上頻段內(nèi)殼測點位置發(fā)生偏差前后外殼響應(yīng)預(yù)估結(jié)果基本一致，40-150 Hz頻段內(nèi)內(nèi)殼測點位置發(fā)生偏差前后外殼響應(yīng)預(yù)估結(jié)果趨勢相同，幅值差別比較大；內(nèi)殼測點位置在周向發(fā)生偏差時，利用準(zhǔn)確的內(nèi)殼測點位置和利用發(fā)生偏差的內(nèi)殼位置處振動數(shù)據(jù)預(yù)估的外殼響應(yīng)曲線基本吻合，表明外殼響應(yīng)預(yù)估結(jié)果準(zhǔn)確性主要受軸向位置偏差影響，這是由于雙層圓柱殼實肋板和環(huán)肋周向連續(xù)分布，軸向離散分布導(dǎo)致的。

圖9 內(nèi)殼測點位置偏差時外殼均方速度級Fig.9 Mean square velocity level of the outer shell during measuring point position of the inner shell changed (a)Position deviation in circumference;(b)Position deviation in axis;(c)Position deviation both in axis and circumference

3 結(jié)論

本文對有限長雙層圓柱殼內(nèi)外殼之間聲振傳遞關(guān)系進行了研究，基于聲振耦合機理，理論推導(dǎo)出內(nèi)外殼之間的聲振傳遞矩陣僅與雙層圓柱殼結(jié)構(gòu)、測點位置和數(shù)量相關(guān)，與激勵源特性無關(guān)。雙層圓柱殼外殼響應(yīng)仿真計算證實激勵源數(shù)量、位置、形式等特性的改變不會對內(nèi)外殼振動傳遞特性產(chǎn)生響應(yīng)，仿真結(jié)果同時表明利用內(nèi)外殼聲振傳遞關(guān)系預(yù)估外殼振動響應(yīng)時，內(nèi)殼測點軸向位置偏差影響大于周向位置偏差影響。該結(jié)論可為后續(xù)利用水下航行器實測振動數(shù)據(jù)，回歸內(nèi)外殼體聲振耦合傳遞關(guān)系，實現(xiàn)基于內(nèi)殼振動響應(yīng)預(yù)估外殼振動響應(yīng)提供理論支撐。

[1]趙智偉.有限長圓柱殼聲與振動的快速預(yù)報方法研究[D].哈爾濱:哈爾濱工程大學(xué),2007. Zhao Zhiwei.Research on fast predicting method of sound and vibration from a finite cylindrical shell[D].Harbin:Harbin Engineering University,2007.(in Chinese)

[2]金廣文,何琳,姜榮俊.基于速度場重構(gòu)的雙層圓柱殼振動傳遞特性試驗研究[J].振動與沖擊,2007,10(26):180-183. Jin Guangwen,He Lin,Jiang Rongjun.Vibration of transmissibility of a ribbed cylindrical double-shell based on reconstructing velocity field[J].Journal of Vibration and Shock,2007,26(10):180-183.(in Chinese)

[3]白振國.雙層圓柱殼舷間聲振耦合特性及控制技術(shù)[D].無錫:中國船舶科學(xué)研究中心,2104. Bai Zhenguo.Vibration and sound coupling behavior and control method of double-layer cylindrical shells[D].Wuxi:China Ship Scientific Research Center,2014.(in Chinese)

[4]陳美霞,金家坤,彭旭等.內(nèi)、外殼對水下雙層圓柱殼聲振性能影響分析[J].船舶力學(xué),2009,4(13):628-634.Chen Meixia,Jin Jiakun,Peng Xu,et al.Influence by the inner shell and the outer shell on sound-vibration characteristics from double cylindrical shells in water[J].Journal of Ship Mechanics,2009,13(4):628-634.(in Chinese)

[5]高菊,陳美霞,陳清坤等.不同流場中雙層圓柱殼層間聲振傳遞特性研究[J].中國艦船研究,2010,5(10):34-39. Gao Ju,Chen Meixia,Chen Qingkun,et al.Sound and vibration transmission characteristics between inner and outer double cylindrical shells in different flow fields[J].Chinese Journal of Ship Research,2010,5(10):34-39.(in Chinese)

[6]夏齊強,陳志堅.殼間聲橋剛度對雙層圓柱殼聲振耦合特性的影響[J].海軍工程大學(xué)學(xué)報,2013,1(25):90-95. Xia Qiqiang,Chen Zhijian.Effect of stiffness of sound bridge between shells on vibro-acoustic coupling characteri-sticcs of double cylindrical shell[J].Journal of Naval University of Engineering,2013,1(25):90-95.(in Chinese)

[7]商德江,何祚鏞.加肋雙層圓柱殼振動聲輻射數(shù)值計算分析[J].聲學(xué)學(xué)報,2001,26(3):193-201. Shang Dejiang,He Zuoyong.The numerical analysis of sound and vibration from a ring-stiffened cylindrical double shell by FEM and BEM[J].Chinese Journal of Acoustics,2001,6(3):193-201.(in Chinese)

Sound and vibration transmission characteristic of finite length double cylindrical shells under radial excitation

YANG Xiao-gang1,LIN Li1,BAI Zhen-guo2,LI Bing1
(1.China Ship Research Development Academy,Beijing 100192,China;2.China Ship Scientific Research Center, Wuxi 214082,China)

The Donnell equation of thin shell was applied to describe double cylindrical shells motion.The effect of the rings,annular pales and the fluid was treated as forces and moments on shells,which were presented in terms of the additional impedance.By using the mode expanding expression of shell’s displacement,the sound and vibration transfer matrix between the inner and outer shell was constituted.The numerical analysis of mean square velocity of the outer shell was carried out.The result verified that the transfer matrix of the inner and outer shell was fixed,which was not influenced by excitation property.It offers the theoretical basis to predict the noise of double cylindrical shells.

sound and vibration transmission;transfer matrix;double cylindrical shell

TB532

A

10.3969/j.issn.1007-7294.2015.05.015

1007-7294（2015）05-0600-09

2015-03-05

楊曉剛（1989-），男，碩士研究生，E-mail：yxg2008300@126.com；

林立（1965-），女，研究員。