吳梓峰

摘? ? 要：為研究結(jié)構(gòu)聲強(qiáng)技術(shù)在船舶振動分析中的應(yīng)用，以一艘航行于珠江水系的雙體游覽船作為研究對象。首先基于結(jié)構(gòu)聲強(qiáng)理論，通過編寫結(jié)構(gòu)聲強(qiáng)計算模塊對該游覽船進(jìn)行結(jié)構(gòu)聲強(qiáng)分析及其可視化;然后分析振動能量在船體結(jié)構(gòu)中的主要流動路徑和耗散方式;最后與實船艙室噪聲測試結(jié)果作比較，為該類船舶減振降噪問題的研究提供參考。

關(guān)鍵詞：結(jié)構(gòu)聲強(qiáng);振動;流向;游覽船

中圖分類號：U674.11 ?? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識碼：A

Abstract： This paper introduces the application of the structural intensity technology in vibration analysis of a traveling cruiser sailing in the Pearl River system. Firstly， based on the structural intensity theory， the structural intensity calculation and visualization of the cruiser are carried out by the calculation module. Then the main flow paths and the dissipation mechanism of vibration energy in hull structure are analyzed. Finally， the calculation results are compared with those of real ship cabin noise test， which can provide reference for the research on vibration and noise reduction of this kind of ship.

Key words： Structural intensity; Vibration; Flow direction; Traveling cruiser

1? ? ?前言

考慮到游覽船的耐久性以及旅客對舒適游覽環(huán)境的要求，關(guān)于船舶振動控制的標(biāo)準(zhǔn)越來越嚴(yán)格。目前研究船舶振動問題的主要手段包括：船體梁固有頻率分析和總振動響應(yīng)計算;上層建筑固有頻率分析和船體-尾部-上層建筑耦合振動計算;局部結(jié)構(gòu)振動分析等。然而在實際的工程應(yīng)用中，仍難以據(jù)此有效和合理地修改結(jié)構(gòu)設(shè)計以減少船舶振動帶來的問題。

結(jié)構(gòu)聲強(qiáng)法同時考慮結(jié)構(gòu)的振動響應(yīng)和激勵引起的結(jié)構(gòu)內(nèi)力，可給出激勵源的位置信息，揭示振動能量在復(fù)雜結(jié)構(gòu)中的流向特性和耗散方式。目前，運(yùn)用結(jié)構(gòu)聲強(qiáng)法研究薄板和加筋板的振動特性已有不少：李凱[2]利用功率流理論研究加筋板振動能量的傳輸、分布和耗散特性;梁日興[3]基于結(jié)構(gòu)聲強(qiáng)法研究了平板、加筋板及箱體結(jié)構(gòu)等的振動特性;Xu X，D[4]研究了不同加強(qiáng)筋布置對復(fù)合材料層板結(jié)構(gòu)聲強(qiáng)分布的影響;何鵬[5]研究了不同激勵下機(jī)艙雙層底的振動傳遞特性。從以上文獻(xiàn)來看，學(xué)者們主要以簡單的板架等結(jié)構(gòu)作為研究對象，得出不少富有參考價值的研究成果，但關(guān)于結(jié)構(gòu)聲強(qiáng)法在全船結(jié)構(gòu)振動分析中的應(yīng)用卻較少。本文提出一種適用于船體結(jié)構(gòu)聲強(qiáng)分析及可視化的數(shù)值模擬方法，可辨識振動能量從激發(fā)源到上部結(jié)構(gòu)的主要流動路徑和耗散方式，為船舶振動控制的研究提供新思路。

2? ? 船體結(jié)構(gòu)聲強(qiáng)計算公式

對于船體結(jié)構(gòu)的有限元分析，通常采用shell單元離散船殼外板、艙壁板、甲板板、上層建筑圍壁和船體骨架的高腹板等;采用beam單元離散加強(qiáng)筋、扶強(qiáng)材、支柱和船體骨架的面板等。若要采用有限元的方法進(jìn)行船體結(jié)構(gòu)的結(jié)構(gòu)聲強(qiáng)分析，首先要建立起beam單元和shell單元的內(nèi)力、彎矩、位移和轉(zhuǎn)角等參數(shù)與一維梁和二維板殼的結(jié)構(gòu)聲強(qiáng)之間的函數(shù)關(guān)系。

3? ?基于PATRAN結(jié)果后處理的結(jié)構(gòu)聲強(qiáng)可視化

本文的研究思路是：首先在MSC/Patran中建立全船三維有限元模型，然后調(diào)用求解器MSC/Natran進(jìn)行全船結(jié)構(gòu)的頻率響應(yīng)分析，并輸出所有beam單元和shell單元的內(nèi)力、彎矩、位移、轉(zhuǎn)角等數(shù)據(jù)（存儲于filename.f06中）;其次，采用Matlab編寫結(jié)構(gòu)聲強(qiáng)計算模塊，讀取和處理filename.f06中的數(shù)據(jù)，并輸出為MSC/Patran后處理模塊中可讀的結(jié)構(gòu)聲強(qiáng)矢量結(jié)果文件（filename.els）;最后，在MSC/Patran的后處理模塊中查看全船的結(jié)構(gòu)聲強(qiáng)矢量圖、幅值圖等，從而實現(xiàn)結(jié)構(gòu)聲強(qiáng)的可視化?；赑atran結(jié)果后處理的結(jié)構(gòu)聲強(qiáng)可視化流程圖，如圖1所示。

在復(fù)雜的三維全船結(jié)構(gòu)有限元模型中，諸如彎曲的船殼外板上的shell單元，不同單元之間的法向量是各不相同的，意味著每個單元自身的局部坐標(biāo)系是完全不同的。而有限元求解器MSC/Natran輸出的結(jié)果文件中，關(guān)于位移、轉(zhuǎn)角、速度、角速度的信息是基于絕對坐標(biāo)系的，且是以節(jié)點的形式輸出的;關(guān)于軸力、剪力、彎矩、扭矩的信息是基于每個單元自身的局部坐標(biāo)系的，且是以單元的形式輸出的。因此，利用上述（1）、（2a）、（2b）公式計算單元的結(jié)構(gòu)聲強(qiáng)矢量之前，首先要將單元各節(jié)點的變形值均分到單元上，即轉(zhuǎn)化為單元的變形值;其次，還需要將基于絕對坐標(biāo)系下的單元變形值，轉(zhuǎn)換為每個單元自身局部坐標(biāo)系下的單元變形值。

4? ? 某游覽船的結(jié)構(gòu)聲強(qiáng)分析