孫 超，何元安 ，商德江 ，尚大晶，劉月嬋 ，張若愚

(1.哈爾濱工程大學(xué) 水聲技術(shù)重點實驗室，哈爾濱 150001;2.船舶系統(tǒng)工程部，北京 100036;3.北京神州普惠科技有限公司，北京 100085)

現(xiàn)代潛艇及其表面均覆蓋有消聲材料，同時也采取覆蓋隔聲材料，其結(jié)果使其輻射噪聲的水平大大降低。圓柱殼結(jié)構(gòu)作為潛艇的典型結(jié)構(gòu)，對于其振動輻射問題的研究在實際中具有重要意義。水下結(jié)構(gòu)振動輻射噪聲的分析方法主要有數(shù)值法和實驗法。前者主要有有限元法、邊界元法、無限元法及統(tǒng)計能量分析法等［1-3］。而實驗方法主要包括均方聲壓法、聲強法和聲全息法等［4-6］。何元安等［7］利用已知模型輻射效率和實驗測量的結(jié)構(gòu)表面均方加速度估算水下輻射聲功率。金廣文等［8-9］研究了雙層圓柱殼體速度場重構(gòu)理論和聲輻射因子的測量來實現(xiàn)水下結(jié)構(gòu)噪聲的預(yù)報，并給出了輻射因子的測量方法。陳明等［10］通過實驗驗證了多源激勵下的輻射聲場可以近似為單(源激勵下的非相干疊加。李冰茹等［11］運用波數(shù)分解法和波疊加法通過測量殼體表面振速來預(yù)報圓柱殼體近場聲壓。

為了尋找一種求解聲輻射問題或散射問題的簡單算法，Koopmann等［12-13］提出了波疊加方法。其主要思想是結(jié)構(gòu)的輻射聲場可以由置于該輻射體內(nèi)部若干個不同大小源強的虛擬源產(chǎn)生的聲場疊加代替，適用于分析任意形狀物體的聲輻射，國外學(xué)者對此已有較深入的研究。Miller等［14］將邊界元法和波疊加法分析剛性體和彈性殼體的散射場進行了比較。Gargouri［15］采用波疊加法在測量潛艇舷側(cè)陣基陣對應(yīng)的內(nèi)部振動后，得到了舷側(cè)陣部位的機械噪聲。Sarkissian［16-17］針對聲場局部全息重建問題開展了研究。在國內(nèi)，于飛等［18］提出了一種穩(wěn)健的全波數(shù)空間聲場重構(gòu)技術(shù)，適用于分析任意形狀物體的聲輻射。畢傳興等［19］利用等效源法對多源聲場和半自由聲場重建等問題進行了研究，并通過音箱實驗驗證了該方法的準(zhǔn)確性。李加慶等［20］則用數(shù)值仿真分析了各因素對波疊加法重建準(zhǔn)確性的影響。薛瑋飛等［21］針對測量工作量較大等問題提出了混合波疊加技術(shù)，主要應(yīng)用于機械噪聲特征提取及聲源識別等方面。目前關(guān)于波疊加法的研究大多處于驗證和發(fā)展階段，工程實際應(yīng)用在國外已經(jīng)開始涉及，國內(nèi)還較少應(yīng)用。本文利用波疊加方法通過圓柱殼體外水聽器陣列測得的復(fù)聲壓來預(yù)報空間其它場點聲壓，并將預(yù)報值和測量值進行對比。仿真分析說明了該方法的正確性。水池實驗和外場湖上實驗分別對柱形換能器聲源和受寬帶激勵源作用的雙層加肋圓柱殼在光殼和敷設(shè)聲學(xué)覆蓋層后的輻射聲場進行預(yù)報，并分析了測量面尺寸對預(yù)報精度的影響。針對工程中解決水下結(jié)構(gòu)聲學(xué)設(shè)備測量及聲場預(yù)報的問題很有研究意義。

1 理論及仿真分析

1.1 基于波疊加法的聲全息技術(shù)

假設(shè)輻射體置于均勻介質(zhì)中，如圖1所示。

圖中S為輻射體的閉合表面，外部區(qū)域為E，內(nèi)部區(qū)域為D，邊界條件為Neumann邊界條件?？紤]聲場隨時間簡諧變化，并且對于無限域中的外部聲輻射問題，還要滿足Sommerfeld輻射條件，以保證無窮遠處不存在反射波。則此類問題可表述為:

圖1 聲場中的輻射體Fig.1 A radiator in sound field

利用SVD，將矩陣G寫成:

其中:對角陣［λ］=diag(λ1，λ2，…，λN)，對角元素為奇異值 λi，且滿足 λ1≥λ2≥…≥λN≥0;［U］為M×M維酉陣，［V］為N×N維酉陣。U和V可寫為:U=(u1，…，uM)，V=(v1，…，vN)。其中 ui為M維向量，稱作左奇異矢量，vi為N維向量，稱作右奇異矢量，并且均具有正交性。得到源強后，空間聲場中任意點的復(fù)聲壓即可表示為:

以上就是基于波疊加法聲全息變換的基本公式，該方法在進行聲場預(yù)報過程中首先利用全息面上的復(fù)聲壓信息求出不同位置處虛擬源強度，然后通過已知的源強求出空間場點的聲學(xué)量實現(xiàn)聲場預(yù)報。

1.2 仿真分析

水下結(jié)構(gòu)輻射噪聲低頻段的主要成分是機械噪聲。因此，對結(jié)構(gòu)在受到力激勵時輻射聲場的預(yù)報問題進行分析是很有必要的。仿真模型是一個有限長圓柱殼體，如圖2所示。結(jié)構(gòu)長5.4 m，半徑1.05 m。模型被劃分為898個結(jié)點，896個四邊形單元。材料密度7 850 kg/m3，楊氏模量2.1 ×1011N/m2，泊松比為0.3。流體介質(zhì)密度1 000 kg/m3，聲速1 500 m/s。測量面為共形的圓柱面，長 10.8 m(-5.4 ～5.4 m)，半徑 1.5 m，周向測點間隔15°，軸向間隔0.49 m，共552個測量點。模型中心和測量面中心位于坐標(biāo)原點。預(yù)報場點位于R=3 m、6 m、12 m的圓周上，周向間隔15°共72個場點。

圖2 圓柱殼模型及預(yù)報場點Fig.2 The cylindrical shell model and prediction points

為更好地分析預(yù)報效果和精度，首先定義誤差:

通過邊界元軟件SYSNOISE計算得到測量面的聲壓值，采用基于波疊加法的全息技術(shù)得到了預(yù)測面72個場點的復(fù)聲壓，并將其與軟件計算得到預(yù)測面聲壓進行對比。仿真分析了兩種情況:① 只有F1激勵(F1=1 000 N);②F1-F4同時激勵 (F1=1 000 N，F(xiàn)2=F3=F4=2 000 N)。為了分析測量面尺寸對預(yù)報結(jié)果的影響，在原測量面中選取中間312個部分測點作為新測量面進行比較，部分測點(-2.95～2.95 m)軸向間隔0.49 m，周向測點間隔15°。其測量面孔徑和結(jié)構(gòu)尺寸相當(dāng)。仿真結(jié)果如圖3所示。激勵力頻率分別為250 Hz，500 Hz和1 000 Hz。以下只給出250 Hz時聲壓幅值的預(yù)報結(jié)果，其它頻率時的誤差見表1。

圖3 預(yù)報聲壓幅值對比Fig.3 Comparisons of prediction pressure amplitude

表1 聲壓預(yù)報誤差Tab.1 Pressure prediction error

從圖示及表格可以看出，應(yīng)用波疊加法進行聲場預(yù)報是一種穩(wěn)健的全波數(shù)算法。無論采用全部測點或部分測點，不同頻率下的預(yù)報結(jié)果與數(shù)值結(jié)果吻合較好，誤差均小于10%。所以在實際工程應(yīng)用中，類似局部近場聲全息，測量面尺寸不必過大，只要包含感興趣的聲源位置即可重建其對應(yīng)的聲場。這不同于基于空間聲場變換(STSF)的近場聲全息技術(shù)，后者要求測量孔徑至少大于被測對象尺寸的兩倍，可見本方法可以節(jié)省很大的工作量，提高了工作效率。

2 實驗驗證

2.1 柱形換能器聲源

為了說明該方法用于水下聲源聲場預(yù)報的有效性，首先在消聲水池進行了驗證性實驗，水池大小16 m×8 m×8 m，在水池壁面和水面上覆有吸聲尖劈，水池消聲頻率f≥2 kHz。實驗所用聲源為柱形換能器聲源，其諧振頻率為3 kHz。換能器發(fā)射面半徑0.15 m，高0.426 m。測量面采用圓柱面，掃描半徑0.635 m，軸向掃描區(qū)域0～1.8 m，采樣間隔0.06 m，周向掃描330°，采樣間隔6°，共1 705個空間測量點，測試示意圖見圖4。實驗分析頻率為3 kHz和5 kHz，測量之前對和水聽器的靈敏度進行校驗。為了比較測量面尺寸對預(yù)報結(jié)果的影響，分別采用了全部測點的聲壓數(shù)據(jù)和部分測點的聲壓數(shù)據(jù)進行預(yù)報。部分測點位于軸向0.54～1.26 m的區(qū)域內(nèi)共715個測量點。預(yù)報區(qū)域同樣為包圍結(jié)構(gòu)的圓柱面，半徑0.735 m，軸向區(qū)域0.6～1.2 m，間隔0.06 m，周向間隔72°，共55 個場點。聲場聲壓級定義為:L=20×log(P/Pr)，單位為dB。P為實測值或預(yù)報值，Pr取10-6Pa為基準(zhǔn)聲壓。圖5為預(yù)報值和測量值的聲壓級比較。

圖4 換能器實驗示意圖Fig.4 The schematic diagram of transducer experiment

由圖比較可知，聲場預(yù)報值和測量值有較好的一致性，整體趨勢也相近。說明該方法可以有效地預(yù)測空間聲場。當(dāng)采用部分測點進行預(yù)報時，整體效果較好，但是個別點誤差較大。其原因可能由于實驗過程水聽器發(fā)生的微小偏移以及采集通道的不一致性，同時該柱面換能器不是理想換能器，其輻射聲場也非沿軸均勻分布的柱面波，這些因素都會產(chǎn)生誤差從而影響預(yù)報精度。

2.2 水下加肋圓柱殼聲源

2.2.1 實驗系統(tǒng)及參數(shù)

外場實驗在千島湖水聲實驗站進行，實驗站水域開闊湖面平靜，附近漁船較少可近似為自由場。實驗?zāi)Ｐ蜑殡p層加肋圓柱殼。尺寸參數(shù)為:圓柱殼長5.4 m，外層殼體直徑2.1 m，外殼鋼板厚0.006 m，內(nèi)層殼體直徑1.5 m，內(nèi)殼鋼板厚0.004 m。測試示意圖參見圖6。

圖5 預(yù)報值測量值聲壓級對比Fig.5 The comparisons of predicted value and measured value

實驗采用寬帶激勵方式，激勵源位于模型中心正橫方向位置。模型通過連接桿吊入水中，連接桿另一端固定在測量轉(zhuǎn)臺上，測量轉(zhuǎn)臺連同模型可在360°內(nèi)旋轉(zhuǎn)。模型中心距離水面14 m，在距模型中心10 m正橫位置布放水聽器測量列陣，最上端1號水聽器入水深度為2 m，最下端19號水聽器入水深度為26 m。1號～4號，16號～19號水聽器間隔2 m，其它水聽器間隔1 m。距模型中心11～15 m位置布放20號～24號水聽器。模型與水聽器陣列相對旋轉(zhuǎn)一周形成柱形測量面。實驗時模型周向每隔15°測量一次。預(yù)測11～15 m處的場點聲壓級。

圖6 圓柱殼實驗示意圖Fig.6 The schematic diagram of cylindrical shell experiment

2.2.2 實驗結(jié)果及分析

預(yù)報圓柱殼體輻射的遠場聲壓級，首先驗證了聲源輻射聲場的球面波特性。取正橫方向20號～24號水聽器聲壓和按球面波傳播規(guī)律衰減的理論值進行比較，結(jié)果見圖7。分析頻率為250 Hz、500 Hz、1 kHz和2 kHz。

可以看到，遠場輻射聲場近似按球面波規(guī)律衰減，測量值和理論值偏差基本小于1 dB。圖8曲線為不同頻率下在未敷設(shè)聲學(xué)覆蓋層時，采用全部測量點預(yù)報的遠場聲壓級，并和實測值結(jié)果進行對比。可見遠場聲壓級預(yù)報結(jié)果和實測結(jié)果在各個頻率下都吻合較好，曲線的變化趨勢基本一致。同時也計算了只用5號～15號水聽器位置的聲場測量數(shù)據(jù)預(yù)報的遠場聲壓級。兩種測點情況下的預(yù)測聲壓級和實測聲壓級之差見圖9。

圖7 圓柱殼實驗示意圖Fig.7 The schematic diagram of cylindrical shell experiment

圖8 光殼時遠場聲壓級的預(yù)報值和實測值對比Fig.8 The comparisons of far field pressure level values

圖9 兩種測點時的預(yù)測值和實測值聲壓級差Fig.9 The difference of pressure level with different points

比較發(fā)現(xiàn)，不論采用全部測點或是部分測點，在不同頻率下的誤差基本小于3 dB，預(yù)報精度較高。說明該方法可以準(zhǔn)確地預(yù)報水下大型圓柱殼體輻射的遠場聲壓級，并且不需要較大的測量陣列。

將圓柱殼外表面敷設(shè)聲學(xué)覆蓋層后，采用同樣的激勵方式激勵殼體。分別對兩種測點情況下的預(yù)報結(jié)果與實測結(jié)果進行比較，其聲壓級差參見圖10?？梢钥闯?，該方法對于含有覆蓋層結(jié)構(gòu)的遠場聲壓級預(yù)報依然有效，整體誤差較小，除個別點稍大外，誤差一般不超過3 dB，預(yù)報精度仍然滿足工程需要?？梢姡捎貌ǒB加法的聲場預(yù)報方法不必使用大尺寸的測量陣列就可以達到較高的預(yù)報精度。總體上看，兩次采用全部測點進行遠場聲壓級預(yù)報的效果要稍好，但需要較多的水聽器和采集通道，工作量較大。因此，該方法在工程上有較強的研究意義。通過圓柱殼體敷設(shè)聲學(xué)覆蓋層后的遠場聲壓級預(yù)報獲得了材料的隔聲性能，這也為評價材料的降噪效果提供了一種新的方法。

圖10 敷設(shè)覆蓋層后兩種測點時的預(yù)測值和實測值聲壓級差Fig.10 The difference of pressure level with different points when the cylindrical shell coated with acoustic covering layer

3 結(jié)論

本文針對水下聲源輻射聲場的預(yù)報問題，采用了基于波疊加方法的聲全息技術(shù)預(yù)報聲源的三維空間聲場。通過水池實驗及外場湖試實驗，比較了柱形換能器聲源和水下圓柱殼體受激振動輻射時聲場的預(yù)報效果。結(jié)果表明:該方法是一種穩(wěn)健有效地聲場預(yù)報方法，預(yù)報值和實測值基本吻合，誤差一般在3 dB以內(nèi)，精度較高能夠滿足工程需要。并在此基礎(chǔ)上發(fā)現(xiàn)采用稍大于結(jié)構(gòu)尺寸的測量陣列即能夠夠較精確地預(yù)測聲場，提高了工作效率。表明其在水下結(jié)構(gòu)聲場預(yù)報和聲學(xué)測量方面有良好的應(yīng)用前景。

［1］商德江，何祚鏞.加肋雙層圓柱殼振動聲輻射數(shù)值計算分析［J］.聲學(xué)學(xué)報，2001，26(3):193-201.

SHANG De-jiang，HE Zuo-yong.The numerical analysis of sound and vibration from a ring-stiffened cylindrical doubleshell by FEM and BEM［J］.Acta Acustic，2001，26(3):193-201.

［2］商德江.復(fù)雜彈性殼體水下結(jié)構(gòu)振動和聲場特性研究［D］.哈爾濱:哈爾濱工程大學(xué)，2004.

［3］劉鐘文，李天勻，駱東平.某潛艇綜合聲納聲腔自噪聲的統(tǒng)計能量方法分析［J］.艦船科學(xué)技術(shù)，2000，28(3):61-64.

LIU Zhong-wen，LI Ttian-yun，LUO Dong-ping.The statistical energy analysis for self noise of submarine sonar［J］.Ship Science and Technology，2000，28(3):61-64.

［4］高 巖，李 兵，戴 江.圓柱殼結(jié)構(gòu)模型水下輻射聲功率測量方法研究［J］.船舶力學(xué)，2010，14(10):1189-1194.

GAO Yan，LI Bing，DAI Jiang.Investigation on underwater radiation sound power measurements from cylindrical shells model［J］.Journal of Ship Mechanics，2010，14(10):1189-1194.

［5］何元安，何祚鏞，姜 軍.基于聲強測量的近場聲全息及其在水下聲源聲輻射分析中的應(yīng)用［J］.聲學(xué)學(xué)報，1996，21(4):297-305.

HE Yua-nan，HE Zuo-yong，JIANG Jun.Near field acoustic holography based on intensity measurement and application in radiating analysis of underwater sound sources.［J］.Acta Acustic，1996，21(4):297-305.

［6］何元安.大型水下結(jié)構(gòu)近場聲全息的理論與實驗研究［D］.哈爾濱:哈爾濱工程大學(xué)，2000.

［7］何元安，楊德森.水下結(jié)構(gòu)輻射噪聲估算方法實驗研究［J］.船舶力學(xué)，1999，3(2):58-62.

HE Yua-nan，YANG De-sen.Experiment study of predicting method for radiating noise by submerged vibrated structure［J］.Journal of Ship Mechanics，1999，3(2):58-62.

［8］金廣文，章林柯，繆旭弘，等.水下加肋雙層圓柱殼體振動傳遞特性分析［J］.振動與沖擊，2011，30(5):218-221.

JIN Guan-gwen，ZhangLin-ke，MIAO Xu-hong，etal.Vibration transmissibility of a submerged cylindrical doubleshell based on reconstructing velocity field［J］.Journal of Vibration and Shock，2011，30(5):218-221.

［9］金廣文，何 琳.水下雙層圓柱殼體結(jié)構(gòu)輻射噪聲實時預(yù)報方法研究［J］.聲學(xué)學(xué)報，2010，35(4):427-433.

JIN Guangwen，HE Lin.Calculation of noises from submerged cylindrical shell structure by radiator factor［J］.Acta Acustic，2010，35(4):427-433.

［10］陳 明，王 斌，曹為午.水下復(fù)雜殼體結(jié)構(gòu)在多源激勵下的振動及聲輻射特性研究［J］.聲學(xué)學(xué)報，2009，34(6):498-505.

CHEN Ming，WANG Bin，CAO Weiwu.Research on vibration and acoustic radiation of the submerged complicated shell structure with multiple excitations［J］.Acta Acustic，2009，34(6):498-505.

［11］李冰茹，王宣銀，葛輝良.圓柱殼體近場輻射噪聲預(yù)報與實驗研究［J］.浙江大學(xué)學(xué)報，2010，44(3):563-568.

LI Bin-gru，WANG Xuan-yin，GE Hui-liang.Prediction and experiment for near-field sound radiation of cylindrical shells［J］.Journal of ZheJiang University(Engineering Science)，2010，44(3):563-568.

［12］Koopmann G H.A method for computing acoustic fields based on the principle of wave superposition［J］.Journal of the Acoustical Society of America，1989，86(5):2433-2438.

［13］ Song L.Numerical errors associated with the method of superposition for computing acoustic fields［J］.Journal of the Acoustical Society of America，1991，89(6):2625-2633..

［14］ Miller R D，Meyer E T.A comparison between the boundary Element method and the wave superposition approach for the analysis of the scattered fields from rigid bodies and elastic shells［J］.Journal of the Acoustical Society of America，1991，89(5):2185-2196.

［15］ Gargouri Y.A study of method of characterizing the effect of internal noise source on submarine flank arrays［J］.Applied Acoustics，1998，53(4):349-367.

［16］ Sarkissian A.Extension of measurement surface in near-field acoustic holography［J］.Journal of the Acoustical Society of America，2004，115(4):1593-1596.

［17］ Sarkissian A.Method of superposition applied to patch nearfield acoustic holography［J］.Journal of the Acoustical Society of America，2005，118(2):671-678.

［18］于 飛，陳 劍，李衛(wèi)兵.一種穩(wěn)健的全波數(shù)聲場重構(gòu)技術(shù)［J］.中國科學(xué)，E 輯，工程科學(xué)，2004，34(9):1069-1080.

YU Fei，CHEN Jian，LI Wei-bing. A robust sound field reconstrction technique in full wave number domain［J］.Science China Technological Sciences ，2004，34(9):1069-1080.

［19］畢傳興.基于等效源法的近場聲全息技術(shù)［J］.中國科學(xué)，E 輯，工程科學(xué)，2005，35(5):535-548.

BIChuan-xi.Near field holography based on wave superposition method ［J］.Science China Technological Sciences，2005，35(5):535-548.

［20］李加慶，陳 進.自由場波疊加噪聲源識別的仿真研究［J］.振動與沖擊，2006，25(4):58-60.

LI Jia-qing，CHEN Jin.Simulation of noice source identif ication based on the wave superpos ition method［J］.Journal of Vibration and Shock，2006，25(4):58-60.

［21］薛瑋飛，陳 進，李加慶.機械噪聲故障特征提取的波疊加法［J］.機械科學(xué)與技術(shù)，2006，25(9):1105-1108.

XUE Wei-fei，CHEN Jin，LI Jia-qing.Wave superposition method for extracting mechanical noise features［J］.Mechan Ical Science and Technology，2006，25(9):1105-1108.

［22］ Fahnline J，Koopmann G H.A numerical solution for the general radiation problem based on the combined methods of superposition and singular value decomposition［J］.Journal of the Acoustical Society of America，1991，90:2808-2818.