基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算

2013-06-04  by:廣州有限元分析、培訓(xùn)中心-1CAE.COM  來源:仿真在線

本文運(yùn)用Sysnoise軟件中的邊界元方法對大型循環(huán)水槽中脈動球源的聲輻射進(jìn)行了數(shù)值計算,得到水槽試驗(yàn)段的自由聲場修正數(shù)據(jù),并與試驗(yàn)測試結(jié)果進(jìn)行了比較。然后,計算偶極子聲源在循環(huán)水槽試驗(yàn)段中的輻射,得到了偶極子聲源的水槽試驗(yàn)段聲場修正數(shù)據(jù)。

陳奕宏 孫紅星 劉竹青 韓用波 來源:LMS
關(guān)鍵字:CAE LMS 循環(huán)水槽 聲場 數(shù)值計算

1 引言
中國船舶科學(xué)研究中心的大型循環(huán)水槽是國內(nèi)開展船舶推進(jìn)器模型噪聲試驗(yàn)研究最為理想試驗(yàn)設(shè)施。就模型噪聲測量試驗(yàn)而言,其具有:低背景噪聲;可以開展整船模型帶推進(jìn)器的噪聲試驗(yàn),模擬伴流場更加接近于實(shí)艇真實(shí)流場;建有專用的水聲測量艙;水聽器與推進(jìn)器更遠(yuǎn),有效測試頻率范圍更寬等優(yōu)勢。該試驗(yàn)設(shè)備自建成以后,已先后完成多型水面船舶和水下航行體的推進(jìn)器模型聲學(xué)試驗(yàn),為我國船舶推進(jìn)器噪聲研究作出了重要貢獻(xiàn)。
開展推進(jìn)器模型噪聲試驗(yàn)的主要目的是預(yù)報自由聲場條件下推進(jìn)器的聲源譜級,而試驗(yàn)室的聲場卻不可能滿足自由場的條件。因此,需要研究循環(huán)水槽試驗(yàn)室聲場特性,進(jìn)而建立循環(huán)水槽試驗(yàn)室聲場與自由聲場之間的修正關(guān)系。早在1987年,國際拖曳水池大會(International Towing Tank Conference,ITTC)就對空泡水筒的自由聲場修正方法進(jìn)行了討論,并提出了自由場修正的試驗(yàn)方法。K.OHTA[2,3](1987年)則采用邊界元方法對日本船舶研究所(Ship Research Institute SRI)的空泡水筒進(jìn)行了聲場數(shù)值計算,計算結(jié)果表明空泡水筒的聲場與自由聲場差異較大,數(shù)值計算的結(jié)果與試驗(yàn)測量結(jié)果趨勢一致。國外的空泡水筒和循環(huán)水槽均采用試驗(yàn)方法建立各自的自由場修正關(guān)系,但由于涉及軍事秘密未能公開其具體過程與修正數(shù)據(jù)。
本文將采用Sysnoise軟件中提供的邊界元方法對循環(huán)水槽試驗(yàn)段中的聲場進(jìn)行數(shù)值計算,并與自由場中聲源的輻射場進(jìn)行對比,得到循環(huán)水槽試驗(yàn)段聲場修正量的數(shù)值計算結(jié)果。本文計算分為兩個部分:首先根據(jù)水槽試驗(yàn)段聲場修正的試驗(yàn)方法,計算脈動球作為單極子聲源的水槽試驗(yàn)段聲場修正結(jié)果,并與試驗(yàn)測量結(jié)果進(jìn)行比較。然后,根據(jù)螺旋槳噪聲偶極子特性,計算偶極子聲源的水槽試驗(yàn)段聲場修正結(jié)果。
本文第二部分將介紹循環(huán)水槽試驗(yàn)段以及水聲測量艙的邊界元模型,第三部分中給出單極子聲源的水槽試驗(yàn)段聲場修正結(jié)果及與試驗(yàn)測量結(jié)果的比較、偶極子聲源的水槽試驗(yàn)段聲場修正結(jié)果。
2 循環(huán)水槽試驗(yàn)段的計算模型
2.1 循環(huán)水槽試驗(yàn)段和水聲艙的模型
循環(huán)水槽試驗(yàn)段的尺寸為10.5m(長)×2.2m(寬)×2.0m(高)。試驗(yàn)段下方設(shè)有水聲測量艙,尺寸為9.5m(長)×2.2m(寬)×2.0m(高)。圖1、圖2分別為循環(huán)水槽試驗(yàn)段和水聲測量艙的照片。

基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+項目圖片圖片1  基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+項目圖片圖片2
圖1 大型循環(huán)水槽     圖2 循環(huán)水槽試驗(yàn)段下方的水聲艙


為了減少水流對噪聲測量的影響,水槽試驗(yàn)段和水聲艙采用有機(jī)玻璃進(jìn)行隔離。但有機(jī)玻璃的透聲性能較好,計算中忽略其對聲傳播的影響。為了減弱壁面反射對噪聲測量的影響,水聲艙的四壁布滿消聲尖劈。因此,水聲艙的計算模型尺寸取為:9.5m(長)×2.0m(寬)×2.0m(高)。圖3為水槽試驗(yàn)段和水聲艙的計算模型,網(wǎng)格尺寸為0.25m。

基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+項目圖片圖片3
圖3 循環(huán)水槽試驗(yàn)段和水聲艙的計算模型


2.2 邊界條件
模型邊界條件是本文計算中的關(guān)鍵。循環(huán)水槽中聲源輻射場計算時,試驗(yàn)段的四壁作為剛性邊界處理。水聲艙的四壁布有消聲尖劈,因此其邊界則作為消聲邊界處理??紤]到試驗(yàn)段前后方分別為水槽的收縮段和擴(kuò)散段,因此試驗(yàn)段的兩個端面作為阻抗邊界處理,其聲阻抗為:

zn0c0

其中,ρ0、c0分別為水的密度和聲速。
    自由聲場中聲源輻射場計算時,所有壁面均作消聲邊界處理,也即壁面無聲反射。
3 計算結(jié)果與分析
3.1 循環(huán)水槽中單極子球聲源輻射聲場計算結(jié)果
首先參照循環(huán)水槽試驗(yàn)段自由聲場修正的試驗(yàn)方法,計算了脈動球源作為單極子聲源的自由聲場和水槽試驗(yàn)段中的聲輻射。然后,對比自由場和水槽試驗(yàn)段中聲輻射計算結(jié)果,得到單極子聲源循環(huán)水槽試驗(yàn)段的聲場修正量。計算時,脈動球源設(shè)置于試驗(yàn)段橫截面的中心位置(距入口6.75m處),接收點(diǎn)取為水聲艙內(nèi)脈動球正下方1.25m處。
圖4和圖5是自由聲場和水槽試驗(yàn)段中脈動球源的聲輻射分布圖。

基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+項目圖片圖片4    基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+應(yīng)用技術(shù)圖片圖片5
圖4 自由聲場中球聲源輻射    圖5 水槽試驗(yàn)段中球聲源輻射


為了驗(yàn)證計算的可靠性,本文開展了循環(huán)水槽試驗(yàn)段自由聲場修正試驗(yàn)。試驗(yàn)中采用直徑為40mm的小球作為聲源,B&K8105水聽器作為接收水聽器,試驗(yàn)中采用1/3Oct分析儀處理噪聲測量結(jié)果。自由場中脈動球聲輻射測量試驗(yàn)在消聲水池中進(jìn)行。循環(huán)水槽中球源輻射噪聲測量時,接收水聽器位于水聲艙內(nèi)聲源正下方1.25m處。最后,根據(jù)試驗(yàn)測量結(jié)果對比,得到循環(huán)水槽試驗(yàn)段聲場修正量。圖6是循環(huán)水槽試驗(yàn)段聲場修正量的計算結(jié)果與試驗(yàn)結(jié)果比較。從圖中可知,在315Hz~415Hz范圍內(nèi)邊界元計算結(jié)果與試驗(yàn)結(jié)果差別在2dB以內(nèi),吻合較好。

基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+應(yīng)用技術(shù)圖片圖片6
圖6 計算結(jié)果與試驗(yàn)結(jié)果比較


3.2 循環(huán)水槽中偶極子聲源輻射聲場計算結(jié)果
眾所周知,螺旋槳噪聲有十分明顯的偶極子聲源特性。為了考察偶極子聲源對循環(huán)水槽自由場修正的影響,本文分別計算了自由聲場和循環(huán)水槽試驗(yàn)段中偶極子聲源的輻射噪聲,得出偶極子聲源循環(huán)水槽試驗(yàn)段的聲場修正量。偶極子聲源的位置與球聲源位置相同,接收點(diǎn)的位置也相同。偶極子聲源強(qiáng)度定義為:(D=10i+0j+10k)

基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+應(yīng)用技術(shù)圖片圖片7    基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+應(yīng)用技術(shù)圖片圖片8
圖7 自由聲場中偶極子聲源輻射 圖8 水槽試驗(yàn)段中偶極子聲源輻射

基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算+應(yīng)用技術(shù)圖片圖片9
圖9 偶極子聲源結(jié)果與單極子聲源結(jié)果比較


圖7和圖8是自由聲場和水槽試驗(yàn)段中偶極子聲源的聲輻射分布圖。圖9是偶極子和單極子聲源計算的循環(huán)水槽試驗(yàn)段聲場修正結(jié)果比較。從圖中可知,在315Hz~415Hz范圍內(nèi)偶極子和單極子聲源循環(huán)水槽試驗(yàn)段聲場修正數(shù)據(jù)隨頻率變化的趨勢一致,兩者在各個頻率點(diǎn)上的修正數(shù)據(jù)存在差異。
4 結(jié)論
本文采用Sysnoise軟件中提供的邊界元方法分別對自由聲場和循環(huán)水槽試驗(yàn)段中脈動球聲源的輻射聲場進(jìn)行數(shù)值計算,并得到了單極子的循環(huán)水槽試驗(yàn)段聲場修正數(shù)據(jù)。計算結(jié)果表明:數(shù)值計算得到的循環(huán)水槽試驗(yàn)段聲場修正數(shù)據(jù)與試驗(yàn)結(jié)果吻合較好。然后,本文對自由聲場和循環(huán)水槽試驗(yàn)段中偶極子聲源的輻射聲場進(jìn)行數(shù)值計算,得到了偶極子聲源的循環(huán)水槽試驗(yàn)段聲場修正數(shù)據(jù)。計算結(jié)果表明:不同的聲源在循環(huán)水槽試驗(yàn)段內(nèi)的輻射聲場不同,得到的循環(huán)水槽試驗(yàn)段聲場修正關(guān)系也不同。


開放分享:優(yōu)質(zhì)有限元技術(shù)文章,助你自學(xué)成才

相關(guān)標(biāo)簽搜索:基于Sysnoise軟件的循環(huán)水槽聲場數(shù)值計算 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn) Abaqus培訓(xùn) Autoform培訓(xùn) 有限元培訓(xùn) Solidworks培訓(xùn) UG模具培訓(xùn) PROE培訓(xùn) 運(yùn)動仿真 

編輯
在線報名:
  • 客服在線請直接聯(lián)系我們的客服,您也可以通過下面的方式進(jìn)行在線報名,我們會及時給您回復(fù)電話,謝謝!
驗(yàn)證碼

全國服務(wù)熱線

1358-032-9919

廣州公司:
廣州市環(huán)市中路306號金鷹大廈3800
電話:13580329919
          135-8032-9919
培訓(xùn)QQ咨詢:點(diǎn)擊咨詢 點(diǎn)擊咨詢
項目QQ咨詢:點(diǎn)擊咨詢
email:kf@1cae.com