響應(yīng)傳導(dǎo)函數(shù)在船舶機(jī)械噪聲預(yù)報(bào)中的應(yīng)用及實(shí)驗(yàn)驗(yàn)證

周軍偉，何 琳，徐榮武，崔立林

（1．海軍工程大學(xué)振動(dòng)與噪聲研究所，武漢 430033；2．海軍工程大學(xué)船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室，武漢 430033）

響應(yīng)傳導(dǎo)函數(shù)在船舶機(jī)械噪聲預(yù)報(bào)中的應(yīng)用及實(shí)驗(yàn)驗(yàn)證

周軍偉1，2，何 琳1，2，徐榮武1，2，崔立林1，2

（1．海軍工程大學(xué)振動(dòng)與噪聲研究所，武漢 430033；2．海軍工程大學(xué)船舶振動(dòng)噪聲重點(diǎn)實(shí)驗(yàn)室，武漢 430033）

據(jù)線性時(shí)不變系統(tǒng)自身具備的響應(yīng)傳導(dǎo)特性，將力－響應(yīng)傳遞函數(shù)模型轉(zhuǎn)化為響應(yīng)－響應(yīng)傳導(dǎo)函數(shù)模型用于船舶機(jī)械噪聲預(yù)報(bào)。理論推導(dǎo)傳導(dǎo)函數(shù)法用于響應(yīng)預(yù)報(bào)的約束條件，分析其用于船舶機(jī)械噪聲實(shí)時(shí)預(yù)報(bào)的可行性。給出傳導(dǎo)函數(shù)的求解方法－工況聯(lián)合法，并對(duì)其應(yīng)用條件進(jìn)行分析論證。雙層圓柱殼體模型試驗(yàn)結(jié)果表明，傳導(dǎo)函數(shù)法用于近、遠(yuǎn)場(chǎng)機(jī)械噪聲預(yù)報(bào)時(shí)均具有較高精確度。

機(jī)械噪聲；實(shí)時(shí)預(yù)報(bào)；頻響函數(shù)；傳導(dǎo)函數(shù)；激勵(lì)源；工況聯(lián)合

對(duì)船舶輻射噪聲的主要成分－機(jī)械噪聲進(jìn)行測(cè)量、評(píng)估意義重大［1］。水下結(jié)構(gòu)輻射聲的計(jì)算方法已有諸多研究［2－5］。該方法按原理可分三類，即解析法［3，6］、數(shù)值法［3，6］及實(shí)驗(yàn)測(cè)量法［3，7］。由于船舶結(jié)構(gòu)的復(fù)雜性，解析法及數(shù)值法用于輻射噪聲評(píng)估受到一定限制，而實(shí)驗(yàn)測(cè)量法基于船舶真實(shí)結(jié)構(gòu)響應(yīng)所得遠(yuǎn)場(chǎng)輻射噪聲結(jié)果則相對(duì)準(zhǔn)確。

實(shí)驗(yàn)測(cè)量法中傳遞函數(shù)法應(yīng)用通常有三步：①由人工激勵(lì)試驗(yàn)獲得傳遞函數(shù)，亦稱頻率響應(yīng)函數(shù)（FRF）；②識(shí)別實(shí)際工況設(shè)備作用于結(jié)構(gòu)的激勵(lì)力；③將前兩步所得傳遞函數(shù)與激勵(lì)力相乘得遠(yuǎn)場(chǎng)輻射噪聲。受測(cè)量手段及設(shè)備安裝條件限制，實(shí)際工況激勵(lì)力較難獲得。針對(duì)傳遞函數(shù)法面臨的激勵(lì)力識(shí)別難題，考慮避開(kāi)求解激勵(lì)力而直接由激勵(lì)力作用系統(tǒng)各點(diǎn)響應(yīng)間的傳遞關(guān)系實(shí)現(xiàn)對(duì)系統(tǒng)未知響應(yīng)的預(yù)報(bào)評(píng)估，將問(wèn)題轉(zhuǎn)化為在設(shè)備激勵(lì)下求解系統(tǒng)近場(chǎng)參考點(diǎn)與遠(yuǎn)場(chǎng)目標(biāo)點(diǎn)響應(yīng)間傳遞關(guān)系。由于該傳遞關(guān)系在一定約束條件下保持不變，因而可據(jù)此傳遞關(guān)系與實(shí)測(cè)所得參考點(diǎn)響應(yīng)對(duì)遠(yuǎn)場(chǎng)目標(biāo)點(diǎn)響應(yīng)實(shí)時(shí)預(yù)報(bào)。為將力與響應(yīng)間傳遞關(guān)系（用H表示，稱傳遞函數(shù)）與響應(yīng)－響應(yīng)間傳遞關(guān)系（用T表示）進(jìn)行區(qū)分，稱一定約束條件下系統(tǒng)響應(yīng)間傳遞關(guān)系為傳導(dǎo)函數(shù)，其輻射噪聲計(jì)算方法則稱傳導(dǎo)函數(shù)法。

本文首先建立物理模型對(duì)傳導(dǎo)函數(shù)原理說(shuō)明、推導(dǎo)，獲得傳導(dǎo)函數(shù)保持不變需滿足的系統(tǒng)約束條件，即本文方法的適用范圍。再通過(guò)數(shù)學(xué)分析論證多工況聯(lián)合法求解響應(yīng)傳導(dǎo)函數(shù)的可行性，給出傳導(dǎo)函數(shù)法用于船舶機(jī)械噪聲實(shí)時(shí)預(yù)報(bào)方案。并通過(guò)雙層圓柱殼體模型試驗(yàn)檢驗(yàn)傳導(dǎo)函數(shù)法的應(yīng)用效果。

1 傳導(dǎo)函數(shù)T推導(dǎo)

系統(tǒng)在設(shè)備激勵(lì)下的響應(yīng)傳導(dǎo)函數(shù)T由激勵(lì)力與結(jié)構(gòu)響應(yīng)間傳遞函數(shù)H推導(dǎo)獲得。本質(zhì)上傳遞函數(shù)H為基于力與響應(yīng)間因果關(guān)系，而傳導(dǎo)函數(shù)T則在系統(tǒng)及激勵(lì)條件下響應(yīng)與響應(yīng)間的共存關(guān)系，二者均為關(guān)于頻率ω的函數(shù)［8］。為表達(dá)簡(jiǎn)潔，將力、加速度、聲壓均默認(rèn)為頻率ω的函數(shù)。

所建設(shè)備激勵(lì)力作用的系統(tǒng)響應(yīng)物理模型，見(jiàn)圖1。用加速度代表響應(yīng)，亦適用其它響應(yīng)（如位移、聲壓等）［9－10］。將一組作用力FA施加于系統(tǒng)，其參考點(diǎn)加速度響應(yīng)集合記為AK，目標(biāo)點(diǎn)加速度響應(yīng)集合記為AU，由于力、參考點(diǎn)、目標(biāo)點(diǎn)數(shù)量均大于1，故FA，AK，AU均為列向量集合。

圖1 系統(tǒng)受力響應(yīng)模型Fig．1 System with force-response relationship

2 系統(tǒng)約束條件

由式（5）知，T由傳遞函數(shù)HUA及HKA的逆矩陣約束，使T保持不變需滿足約束條件：①系統(tǒng)為線性時(shí)不變系統(tǒng)即滿足響應(yīng)疊加原理，傳遞函數(shù)能反映系統(tǒng)固有特性，不隨激勵(lì)力幅值變化［10－12］。②激勵(lì)點(diǎn)、參考點(diǎn)、預(yù)報(bào)點(diǎn)三者數(shù)量及位置固定不變即滿足矩陣HKA，HUA中元素?cái)?shù)值不變。③參考點(diǎn)數(shù)量大于等于激勵(lì)力數(shù)量即滿足矩陣HKA的逆存在條件。

滿足系統(tǒng)約束條件即能用試驗(yàn)所測(cè)響應(yīng)傳導(dǎo)函數(shù)結(jié)合實(shí)測(cè)參考點(diǎn)響應(yīng)計(jì)算目標(biāo)點(diǎn)的聲響應(yīng)。將設(shè)備激勵(lì)的船舶與周圍空間聲傳播介質(zhì)作為研究系統(tǒng)（圖示橢圓包圍的區(qū)域），并近似為線性時(shí)不變系統(tǒng)，見(jiàn)圖2。由于系統(tǒng)中各激勵(lì)設(shè)備數(shù)量、位置均固定（用EA表示），只要選足夠數(shù)量參考點(diǎn)（用XK表示），將其與目標(biāo)點(diǎn)（用XU表示）位置固定，則兩點(diǎn)響應(yīng)間傳導(dǎo)函數(shù)T保持不變。通過(guò)試驗(yàn)測(cè)得T，當(dāng)船舶航行或停靠時(shí)無(wú)論處于何種設(shè)備激勵(lì)組合，用布置船體的參考點(diǎn)XK響應(yīng)數(shù)據(jù)與T聯(lián)合計(jì)算即可求出目標(biāo)點(diǎn)XU的聲響應(yīng)。

圖2 設(shè)備激勵(lì)下船舶噪聲輻射模型Fig．2 Ship mechanical noisemodel with equipmentworking

3 工況聯(lián)合法求解傳導(dǎo)函數(shù)

3.1 原理分析

因分析知，傳導(dǎo)函數(shù)T為頻率ω的函數(shù)，見(jiàn)圖3。對(duì)激勵(lì)頻率為ω0的作用力F（ω0），參考點(diǎn)與目標(biāo)點(diǎn)響應(yīng)的傳導(dǎo)函數(shù)T也是關(guān)于頻率ω0及倍頻的函數(shù)，即T（ω0），T（2ω0），…；而對(duì)設(shè)備的寬帶激勵(lì)如激勵(lì)頻率為ω1－ω2，則亦可將寬帶的傳導(dǎo)函數(shù)T分解為T（ω1），T（ω1＋1），T（ω1＋2），…，T（ω2）進(jìn)行逐個(gè)求解。

圖3 傳導(dǎo)函數(shù)頻域?qū)傩訤ig．3 Frequency dependency of TF

圖4 力－響應(yīng)及響應(yīng)－響應(yīng)關(guān)系模型Fig．4 Force-response vs．response-response

設(shè)三臺(tái)設(shè)備的工作頻率為ω0，其激勵(lì)力幅值分別為f1，f2，f3；三參考點(diǎn)位置響應(yīng)分別為a1，a2，a3；目標(biāo)點(diǎn)響應(yīng)為au，見(jiàn)圖4。由響應(yīng)－響應(yīng)傳導(dǎo)函數(shù)得

由于試驗(yàn)條件下式（6）中參考點(diǎn)響應(yīng)［a1a2

a3］與目標(biāo)點(diǎn)響應(yīng)au均可通過(guò)傳感器直接測(cè)得，但T1u，

T2u，T3u三未知數(shù)并不能由1個(gè)方程計(jì)算獲得。據(jù)傳導(dǎo)函數(shù)T保持不變時(shí)應(yīng)滿足的約束條件，可在約束條件下設(shè)計(jì)不同工況，記為m組。只要使每組工況各設(shè)備激勵(lì)力幅值互不相同，即能連續(xù)獲得一系列響應(yīng)傳導(dǎo)方程，即

式中：上標(biāo)（1）、（2）…（m）分別代表第1～m組工況所測(cè)響應(yīng)。將式（7）變形，得傳導(dǎo)函數(shù)T的解為

求解傳導(dǎo)函數(shù)T需式（7）本身為適定或超定，即滿足工況數(shù)m等于或大于參考點(diǎn)數(shù)k的條件。且對(duì)式（8）參考點(diǎn)響應(yīng)矩陣AK的求逆運(yùn)算盡可能精確，即矩陣AK趨于滿秩，其行、列向量間保持線性無(wú)關(guān)。由式（9）看出，參考點(diǎn)響應(yīng)由激勵(lì)力與傳遞函數(shù)共同作用結(jié)果。不同工況的激勵(lì)力幅值間相關(guān)性直接決定矩陣AK的行向量間相關(guān)性；傳遞函數(shù)矩陣各列之間相關(guān)性直接決定矩陣AK的列向量間相關(guān)性。因此，為獲得參考點(diǎn)矩陣AK滿秩結(jié)果，工況聯(lián)合法求解傳導(dǎo)函數(shù)T時(shí)應(yīng)滿足的約束條件為：①試驗(yàn)工況數(shù)大于或等于參考點(diǎn)數(shù)，即滿足式（7）適定或超定。②不同工況的激勵(lì)力幅值保持線性無(wú)關(guān)，即滿足矩陣AK行向量間線性無(wú)關(guān)。③參考點(diǎn)應(yīng)遵循分散、非對(duì)稱布置原則，且應(yīng)避開(kāi)低頻模態(tài)節(jié)點(diǎn)，即滿足矩陣AK列向量間線性無(wú)關(guān)。

3.2 工況設(shè)置

求解某頻率處傳導(dǎo)函數(shù)時(shí)，為獲得多組線性無(wú)關(guān)的激勵(lì)力，可通過(guò)將工作在該頻率的所有設(shè)備依次開(kāi)啟，或?qū)⒃O(shè)備遍歷組合開(kāi)啟，或在保持設(shè)備同時(shí)開(kāi)啟情況下分別獨(dú)立改變各設(shè)備激勵(lì)幅值實(shí)現(xiàn)。以三臺(tái)設(shè)備為例，工況設(shè)置方案分別對(duì)應(yīng)于表1～表3。表中數(shù)字0代表設(shè)備關(guān)閉，1，2，3分別代表設(shè)備開(kāi)啟的各作用力相對(duì)幅值。

表1 依次獨(dú)立開(kāi)啟Tab.1 Turn on separately in sequence

表2 遍歷組合開(kāi)啟Tab.2 Turn on in combination and include all ones

表3 同時(shí)開(kāi)啟獨(dú)立變幅值Tab.3 Tu rn on sim ultaneously but vary individually

用以上工況設(shè)置方案可實(shí)現(xiàn)由較少試驗(yàn)工況數(shù)求得傳導(dǎo)函數(shù)T，進(jìn)而用于設(shè)備任意組合工況的遠(yuǎn)場(chǎng)響應(yīng)評(píng)估，且該特性會(huì)隨設(shè)備數(shù)量增多呈越顯著優(yōu)勢(shì)。設(shè)在某頻率處施加激勵(lì)設(shè)備數(shù)量為10，則該10臺(tái)設(shè)備對(duì)應(yīng)的組合開(kāi)啟工況共N種，即

由式（10）求出N＝1 023。用工況聯(lián)合法在試驗(yàn)中只需將10臺(tái)設(shè)備分部開(kāi)啟建立10種線性無(wú)關(guān)工況即可求得傳導(dǎo)函數(shù)T，可用于1023種工況的遠(yuǎn)場(chǎng)響應(yīng)計(jì)算。

4 噪聲預(yù)報(bào)實(shí)施流程

對(duì)傳導(dǎo)函數(shù)T應(yīng)用條件及求解方法論證分析后，給出傳導(dǎo)函數(shù)用于船舶機(jī)械噪聲實(shí)時(shí)預(yù)報(bào)的實(shí)施流程為：①將船舶機(jī)械設(shè)備按工作頻率分組；②確定各工作頻率下機(jī)械設(shè)備數(shù)量；③確定各工作頻率下參考點(diǎn)數(shù)量及位置；④在遠(yuǎn)場(chǎng)目標(biāo)位置布放水聽(tīng)器，開(kāi)展多組線性無(wú)關(guān)的訓(xùn)練工況試驗(yàn)，通過(guò)所測(cè)得參考點(diǎn)及目標(biāo)點(diǎn)響應(yīng)求解各工作頻率的傳導(dǎo)函數(shù)存入數(shù)據(jù)庫(kù)；⑤在船舶航行時(shí)據(jù)參考點(diǎn)實(shí)時(shí)響應(yīng)數(shù)據(jù)由數(shù)據(jù)庫(kù)中調(diào)取對(duì)應(yīng)的傳導(dǎo)函數(shù)，二者相乘獲得船舶遠(yuǎn)場(chǎng)機(jī)械輻射噪聲。

5 試驗(yàn)驗(yàn)證

5.1 近場(chǎng)結(jié)構(gòu)響應(yīng)預(yù)報(bào)

試驗(yàn)所用雙層圓柱殼體模型外部尺寸及內(nèi)部設(shè)備布置見(jiàn)圖5、圖6。為檢驗(yàn)傳導(dǎo)函數(shù)法用于響應(yīng)預(yù)報(bào)的一般性，對(duì)模型在空氣中外殼加速度響應(yīng)進(jìn)行預(yù)報(bào)驗(yàn)證。試驗(yàn)場(chǎng)所為半消聲室，用于激勵(lì)模型殼體的三臺(tái)設(shè)備分別為電機(jī)、1＃激振器、2＃激振器。試驗(yàn)所用工況設(shè)置見(jiàn)表4，參考點(diǎn)及目標(biāo)點(diǎn)位置見(jiàn)表5。

圖5 雙層圓柱殼體外形尺寸（單位：mm）Fig．5 External dimension of cylindrical double-shellmodel

圖6 圓柱殼體內(nèi)部設(shè)備布Fig．6 Inner arrangement of cylindricalmodel

表4 陸上試驗(yàn)工況設(shè)置Tab.4 The training and validation conditions on land

表5 參考點(diǎn)及目標(biāo)點(diǎn)位置Tab.5 The locations of reference and target points

由于三臺(tái)設(shè)備激勵(lì)頻率不重合，理論上可只用1個(gè)參考點(diǎn)計(jì)算傳導(dǎo)函數(shù)進(jìn)行目標(biāo)點(diǎn)響應(yīng)預(yù)報(bào)。將參考點(diǎn)數(shù)量設(shè)為大于激勵(lì)源數(shù)量。由式（8）據(jù)訓(xùn)練工況下所測(cè)參考點(diǎn)及目標(biāo)點(diǎn)加速度響應(yīng)求出傳導(dǎo)函數(shù)T，將驗(yàn)證工況的參考點(diǎn)加速度與傳導(dǎo)函數(shù)T相乘求得目標(biāo)點(diǎn)預(yù)報(bào)響應(yīng)，并將其與驗(yàn)證工況實(shí)測(cè)目標(biāo)點(diǎn)響應(yīng)進(jìn)行比較，結(jié)果見(jiàn)圖7。由圖7看出，在設(shè)備激勵(lì)頻率處預(yù)報(bào)結(jié)果與實(shí)測(cè)值符合很好，各激勵(lì)頻率處預(yù)報(bào)誤差分別為：61 Hz／2．45 dB；180 Hz／0．11 dB；300 Hz／0．56 dB。其它頻率處因設(shè)備未進(jìn)行有效激勵(lì)，因而區(qū)域信號(hào)信噪比較低，預(yù)報(bào)結(jié)果意義不大。

5.2 水下輻射噪聲評(píng)估

水下試驗(yàn)的消聲水池見(jiàn)圖8。雙層圓柱殼體模型位于水下5 m，模型內(nèi)安裝3臺(tái)激振器，在同深度距模型中心10 m處布設(shè)水聽(tīng)器，用于測(cè)量雙層圓柱殼體遠(yuǎn)場(chǎng)輻射聲壓，同時(shí)作為輻射噪聲評(píng)估的目標(biāo)點(diǎn)。模型外殼分別布設(shè)20個(gè)加速度計(jì)及6個(gè)水聽(tīng)器用于近場(chǎng)參考點(diǎn)響應(yīng)選取，試驗(yàn)中各激振器工況設(shè)置見(jiàn)表6，各工況參考點(diǎn)與遠(yuǎn)場(chǎng)目標(biāo)點(diǎn)響應(yīng)同測(cè)得。

表6 水下試驗(yàn)工況設(shè)置Tab.6 The training and validation conditions in underwater test

圖7 結(jié)構(gòu)加速度響應(yīng)預(yù)報(bào)結(jié)果與實(shí)測(cè)值比較Fig．7 Comparison ofmeasured vs．estimated value on acceleration responses

圖8 雙層圓柱殼體水下試驗(yàn)方位布置Fig．8 The arrangement diagram of underwater experiment

圖9 外殼加速度作為參考點(diǎn)響應(yīng)的預(yù)報(bào)結(jié)果（單頻點(diǎn)）Fig．9 Comparison ofmeasured vs．estimated value using accelerometers measurement：single frequency point

圖10 外殼加速度作為參考點(diǎn)響應(yīng)預(yù)報(bào)結(jié)果（1／3倍頻程）Fig．10 Comparison ofmeasured vs．estimated value using accelerometers measurement：one-third octave

圖11 外殼水聽(tīng)器作為參考點(diǎn)響應(yīng)預(yù)報(bào)結(jié)果（單頻點(diǎn)）Fig．11 Comparison ofmeasured vs．estimated value using hydrophones measurement：single frequency point

圖12 外殼水聽(tīng)器作為參考點(diǎn)響應(yīng)的預(yù)報(bào)結(jié)果（1／3倍頻程）Fig．12 Comparison ofmeasured vs．estimated value using hydrophones measurement：one-third octave

用外殼加速度計(jì)作為參考點(diǎn)響應(yīng)所得遠(yuǎn)場(chǎng)輻射聲壓預(yù)報(bào)值與實(shí)測(cè)值比較結(jié)果見(jiàn)圖9、圖10。用外殼水聽(tīng)器作為參考點(diǎn)響應(yīng)所得遠(yuǎn)場(chǎng)輻射聲壓預(yù)報(bào)值與實(shí)測(cè)值比較結(jié)果見(jiàn)圖11、圖12。各試驗(yàn)結(jié)果分別以單頻點(diǎn)及1／3倍頻程形式表示。

5.3 試驗(yàn)結(jié)果分析

由雙層圓柱殼體模型試驗(yàn)結(jié)果，可知：

（1）傳導(dǎo)函數(shù)法用于響應(yīng)預(yù)報(bào)精度較高。在中高頻段（本試驗(yàn)?zāi)Ｐ蛯?duì)應(yīng)150 Hz以上頻段）的頻帶級(jí)預(yù)報(bào)誤差小于2 dB，低頻段（本試驗(yàn)?zāi)Ｐ?50 Hz以下頻段）的頻帶級(jí)預(yù)報(bào)誤差小于5 dB，中高頻段預(yù)報(bào)結(jié)果明顯優(yōu)于低頻段。其原因可能為：①模型結(jié)構(gòu)在低頻區(qū)模態(tài)密度小，結(jié)構(gòu)各參考點(diǎn)在單頻點(diǎn)響應(yīng)被有限模態(tài)主導(dǎo)，導(dǎo)致各頻點(diǎn)參考點(diǎn)響應(yīng)間相關(guān)性較高，使傳導(dǎo)函數(shù)T計(jì)算結(jié)果不準(zhǔn)確。隨頻率升高結(jié)構(gòu)模態(tài)會(huì)越密集，每個(gè)頻點(diǎn)均有較多模態(tài)對(duì)結(jié)構(gòu)振動(dòng)產(chǎn)生影響，各參考點(diǎn)響應(yīng)間能保持互相獨(dú)立，所求傳導(dǎo)函數(shù)T精確度高，預(yù)報(bào)結(jié)果準(zhǔn)確。由于模態(tài)密度與結(jié)構(gòu)尺度關(guān)聯(lián)，對(duì)如船舶大尺度結(jié)構(gòu)，其模態(tài)密度對(duì)傳導(dǎo)函數(shù)法預(yù)報(bào)結(jié)果影響尚需進(jìn)一步研究。②低頻段由于水中流固耦合影響，水在外殼體產(chǎn)生的等效附加質(zhì)量較大，導(dǎo)致阻抗變大，結(jié)構(gòu)不易被激勵(lì)，導(dǎo)致所測(cè)響應(yīng)信號(hào)信噪比差，影響計(jì)算及預(yù)報(bào)結(jié)果。

（2）1／3倍頻程的預(yù)報(bào)結(jié)果優(yōu)于單頻點(diǎn)。由圖9、圖11看出，預(yù)報(bào)結(jié)果主要在響應(yīng)極低的頻率點(diǎn)（反共振頻率）及響應(yīng)極高的頻率點(diǎn)（共振頻率）處誤差較大，由于反共振頻率處信號(hào)信噪比低，共振頻率處信號(hào)測(cè)量誤差被放大，二者皆會(huì)導(dǎo)致預(yù)報(bào)誤差變大。而用1／3倍頻程處理預(yù)報(bào)結(jié)果時(shí)能在一定程度上抵消各頻點(diǎn)處由測(cè)量隨機(jī)誤差導(dǎo)致的影響，使預(yù)報(bào)結(jié)果更準(zhǔn)確。

（3）用外殼水聽(tīng)器的數(shù)據(jù)作為參考點(diǎn)響應(yīng)所得預(yù)報(bào)結(jié)果優(yōu)于用外殼加速度計(jì)數(shù)據(jù)。原因?yàn)榭紤]加速度計(jì)易受安裝條件影響，當(dāng)加速度計(jì)軸線偏離結(jié)構(gòu)表面法向時(shí)會(huì)導(dǎo)致測(cè)量誤差。

（4）由參考點(diǎn)位置選取對(duì)試驗(yàn)結(jié)果影響看出，參考點(diǎn)選取原則是成立的，即參考點(diǎn)位置盡量分散，非對(duì)稱，避開(kāi)低階模態(tài)節(jié)點(diǎn)。

6 結(jié) 論

本文據(jù)線性時(shí)不變系統(tǒng)自身存在的響應(yīng)傳導(dǎo)特性，研究響應(yīng)傳導(dǎo)關(guān)系用于船舶機(jī)械噪聲評(píng)估的可行性，通過(guò)理論推導(dǎo)得出傳導(dǎo)函數(shù)法的應(yīng)用約束條件。給出傳導(dǎo)函數(shù)求解方法－工況聯(lián)合法，并對(duì)工況聯(lián)合法的應(yīng)用條件進(jìn)行分析論證。雙層圓柱殼體模型試驗(yàn)表明，傳導(dǎo)函數(shù)法預(yù)報(bào)結(jié)構(gòu)響應(yīng)及評(píng)估遠(yuǎn)場(chǎng)輻射噪聲效果均較好。

［1］王之程，陳宗岐．艦船噪聲測(cè)量與分析［M］．北京：國(guó)防工業(yè)出版社，2004：3－5．

［2］金廣文．雙層圓柱殼體水下噪聲實(shí)時(shí)預(yù)報(bào)方法研究［D］．武漢：海軍工程大學(xué)，2008．

［3］湯智胤．潛艇聲隱身態(tài)勢(shì)評(píng)估研究［D］．武漢：海軍工程大學(xué)，2011．

［4］王斌．基于表面振動(dòng)監(jiān)測(cè)的大型水下結(jié)構(gòu)輻射噪聲預(yù)報(bào)研究［D］．上海：上海交通大學(xué)，2008．

［5］彭旭，駱東平．船舶結(jié)構(gòu)建模及水下振動(dòng)和輻射噪聲預(yù)報(bào)［J］．噪聲與振動(dòng)控制，2003（6）：9－12．

PENG Xu，LUO Dong-ping．The modeling of the ship and numerical simulation of coupled vibro-acoustic behavior by FEM／BEM［J］．Noise and Vibration Control，2003（6）：9－12．

［6］孫雪榮，朱錫．船舶水下結(jié)構(gòu)噪聲的研究概況與趨勢(shì)［J］．振動(dòng)與沖擊，2005，24（1）：106－113．

SUN Xue-rong，ZHU Xi．Status and trend of underwater noise research on ship structure［J］．Journal of Vibration and Shock，2005，24（1）：106－113．［7］何元安，楊德森．水下結(jié)構(gòu)輻射噪聲估算方法實(shí)驗(yàn)研究［J］．船舶力學(xué)，1999，3（2）：58－61．

HE Yuan-an，YANG De-sen．Experimental study on underwater structure radiated noise evaluation［J］．Journal of Ship Mechanics，1999，3（2）：58－61．

［8］Maia N M，Silva JM．The transmissibility concept in multidegree-of freedom systems［J］．Mechanical Systems and Signal Processing，2001，15（1）：129－137．

［9］張磊，曹躍云，楊自春，等．船舶噪聲源的傳遞率矩陣識(shí)別方法及試驗(yàn)驗(yàn)證［J］．振動(dòng)工程學(xué)報(bào)，2013，26（2）：291－297．

ZHANG Lei，CAO Yue-yun，Yang Zi-chun，et al．Transmissibility matrix method for identification of ship noise source and corresponding test verification［J］．Journal of Vibration Engineering，2013，26（2）：291－297．

［10］張磊，曹躍云，楊自春，等．雙層圓柱殼體水下振動(dòng)噪聲結(jié)構(gòu)傳遞路徑分析［J］．振動(dòng)與沖擊，2012，31（20）：12－16．

ZHANG Lei，CAO Yue-yun，YANG Zi-chun，et al．Structural transfer path analysis for vibration and noise of a submerged cylindrical double-shell［J］．Journal of Vibration and Shock，2012，31（20）：12－16．

［11］傅志方，華宏星．模態(tài)分析理論與應(yīng)用［M］．上海：上海交通大學(xué)出版社，2000：1－20．

［12］曹躍云，張磊，楊自春，等．船舶振動(dòng)噪聲源傳遞路徑分析及試驗(yàn)驗(yàn)證［J］．振動(dòng)與沖擊，2013，32（22）：159－162．

CAO Yue-yun，ZHANG Lei，YANG Zi-chun，et al．A new OPA model for ship noise sources and test validation［J］．Journal of Vibration and Shock，2013，32（22）：159－162．

Practical app lication and experimental verification of transm issibility function in ship mechanical noise prediction

ZHOU Jun-wei1，2，HE Lin1，2，XU Rong-wu1，2，CUILi-lin1，2
（1．Institute of Noise＆vibration，Naval University of Engineering，Wuhan 430033，China；2．National Key Laboratory on ship Vibration＆Noise，Wuhan 430033，China）

The traditional response estimation model is based on force-response relationship．According to the transmissibility concept in linear time invariant system，a new response estimation method by the name of transmissibility function（TF）method was provided based on response-response relationship．The principle of themethod was interpreted in mathematic formula and the restriction condition was presented．Then，an engineering technique called operational conditions combination was presented for calculating the transmissibility function（TF）of vibration responses，and the guideline on how to devise the scheme of working conditions was discussed．The performance of the developed transmissibility function method was validated by a cylindrical double-shell vibration test．

mechanical noise；real-time prediction；frequency response function（FRF）；transmissibility function（TF）；exciting force；operational conditions combination

TB53

：A

10．13465／j．cnki．jvs．2014．22．014

國(guó)家自然科學(xué)基金（51209214）；國(guó)防科技重點(diǎn)實(shí)驗(yàn)室基金（SYSZC2012006）

2013－08－15 修改稿收到日期：2013－11－28

周軍偉男，博士生，1983年生

徐榮武男，博士，副研究員，1980年生