基于改進板塊元法的潛艇目標強度預報仿真

孫乃葳， 李建辰， 萬亞民， 趙 罡， 呂 維， 范若楠

（1. 中國船舶重工集團公司 第705研究所， 陜西 西安， 710077； 2. 水下信息與控制重點實驗室， 陜西 西安，710077）

基于改進板塊元法的潛艇目標強度預報仿真

孫乃葳1，2， 李建辰1， 萬亞民1， 趙 罡1， 呂 維1， 范若楠1

（1. 中國船舶重工集團公司 第705研究所， 陜西 西安， 710077； 2. 水下信息與控制重點實驗室， 陜西 西安，710077）

針對板塊元法計算潛艇目標強度時積分分母可能為零， 從而導致計算結果不穩定的問題， 文中將Gordon積分算法用于潛艇目標強度的預報仿真， 并針對復雜目標簡化面元遮擋判斷流程， 首先對剛性球體進行建模計算，其次對Benchmark潛艇進行3D建模， 運用改進板塊元法對潛艇目標強度進行預報。仿真結果表明， 改進模型的計算結果更加穩定， 能夠較好反映潛艇目標強度的特征。

潛艇； 板塊元； Gordon積分； 遮擋判斷； 目標強度預報

0 引言

有效獲取潛艇目標強度特性是主動聲吶探測的關鍵技術之一， 準確掌握目標在不同態勢下的目標特性， 對魚雷等水中兵器的目標檢測、參數估計具有十分重要的意義［1］。當今， 魚雷和聲吶技術的發展， 要求目標回波特性預報的精度更高。目標強度預報關鍵在于求解 3D流體空間中的目標受聲波激勵產生的滿足表面邊界條件、波動方程和輻射條件的散射聲場［2］。

分離變量法雖然是求解波動方程的重要方法， 但它只能應用于那些表面能用正交曲線坐標表示的規則形狀物體， 從而限制了它的應用范圍。Helmholtz積分方法可以在遠場條件下， 求解形狀較復雜物體的散射聲場。Kirchhoff近似法是為了避免 Helmholtz積分方程求解的復雜性而提出的， 其適用于高頻。板塊元方法就是基于Kirchhoff近似法， 用一組平面板塊元近似目標曲面， 將所有面元的散射聲場疊加近似得到總散射聲場。由于板塊元方法把散射聲場的積分運算轉化成代數運算， 大大提高了計算速度。馬忠誠［3］依據艇體線型值對表面進行插值計算， 完成曲面擬合， 采用Kirchhoff近似積分方法進行目標強度的計算， 該方法計算過程較為復雜， 不適應快速預報的需要； 范軍［4］等人建立了面向實際工程應用的潛艇回波特性快速預報模型， 采用Kirchhoff近似方法， 將潛艇表面分成多組三角形面元， 借助傅里葉積分運算， 將面積分轉成圍線積分， 從而獲得一種精確、快速的目標回聲預報模型， 但該方法在板塊元空間坐標劃分時會出現積分分母為零的情況， 因而由面元計算得出的目標強度可能存在奇異值， 從而導致計算不穩定。

文中應用將 Gordon積分法引入潛艇目標強度板塊元預報方法， 并改進面元遮擋判斷流程，運用改進模型對潛艇目標強度進行預報， 預報結果表明， 該改進模型避免了計算奇異值的出現，可有效用于潛艇目標強度的快速預報。

1 水下目標強度預報方法

1.1 板塊元方法

根據 Kirchhoff近似， 遠場條件下（ 1）kr? 得到遠場條件下的目標強度［5］

目標表面所有面元通過坐標變換統一變換到某個確定的平面上， 應用傅里葉積分變換［6］，

圖1 網格劃分及坐標轉換Fig. 1 Meshing and coordinates conversion

可以將面元積分轉化為圍線積分， 對所有面元的散射聲場求和就得到目標散射聲場的近似值

1.2 改進算法建模

用 L表示式（3）中的圍線積分部分， 將式（4）帶入得

其中

根據羅比達法則， 下式成立

此時積分有界， 該方法具有穩定解。

將 Gordon積分算法帶入， 可得基于改進板塊元算法的目標強度預報模型［8］

1.3 遮擋算法改進

根據物理聲學法的假設， 入射聲波照射 3D物體時， 陰影區域的面元對散射沒有貢獻， 因此陰影區的三角形面元不參與目標強度的計算。定義法向朝向聲波入射方向的面元稱之為亮區面元，否則稱之為暗區面元。因此式（3）中， 參與計算的G H× 個面元應為亮區面元。板塊元法本質上是將 3D實體模型表面離散成三角形面元， 總的目標強度近似為各個有散射貢獻的面元目標強度之和。因此對于復雜形狀目標， 不能忽略面元遮擋的影響， 必須對各個頂點及面元進行判斷， 去除遮擋面元對聲波反射的影響。

該方法需要對目標表面所有面元進行兩兩判斷， 才能在計算目標的回波強度時， 把被遮擋的板塊排除在外。

為避免對目標表面所有的面元兩兩進行遮擋判斷導致的計算量大， 效率低的問題， 這里考慮應用模型都是凸表面幾何體， 對亮區面元的遮擋問題簡化為僅與亮區面元而不對所有面元進行兩兩遮擋判斷［9］。參與目標強度計算的面元篩選流程見圖2。

圖2 面元篩選流程Fig. 2 Selection process of planar elements

圖3 面元遮擋示意圖Fig. 3 Schematic of planner occlusion

第1步: 判別入射波是否首先到達面元1所在平面。若入射波首先到達面元 2所在平面， 有可能遮擋面元1。

第2步: 判別交點P是否在面元2內部。若在， 則判斷面元1被面元2遮擋。

這里引入一個參數α， 令

由于

將式（10）代入式（9）， 整理得

第3步: 判別交點P是否在面元2內部。

2 仿真驗證與結果分析

2.1 算法驗證

對于任意的曲面， 都可以近似地由1組三角形表示。文中仿真計算將模型表面劃分三角形板塊元， 即 N= 3。

可見目標強度值在8 dB左右起伏。為得到剛球體的目標強度， 考慮球體的對稱結構， 因此可將各個角度的計算目標強度進行算術平均， 得到統計意義上的平均目標強度值

如圖5所示。

圖4 剛性球體目標強度Fig. 4 Target strength of rigid sphere

圖5 改進板塊元法與分離變量法對比Fig. 5 Comparison between improved planar element method and variable separation method

由圖5可知， 改進板塊元法得到的結果同解析解計算結果的變化逐漸趨向一致， 二者計算結果都在8 dB左右， 差值極小。改進板塊元法與分離變量法直接得到的剛性球體目標強度值的差值較小， 從而驗證了改進板塊元算法的正確性。

2.2 Benchmark潛艇建模仿真

用UG對國際通用標準的Benchmark潛艇（以下簡稱潛艇）進行3D建模， 艇長62 m， 艇身圍殼直徑7.5 m， 如圖6所示。

圖6 Benchmark潛艇3D模型Fig. 6 3D model of Benchmark submarine

為便于將模型文件導入 ANSYS分析， 需要將文件格式轉為.igs格式文件。為方便計算， 網格劃分類型選擇Tri（三角形）。由于僅需要潛艇表面殼層的節點與拓撲信息， 因此 ANSYS中選擇shell劃分網格。網格屬性對計算結果沒有影響，文中定義網格屬性為3D 4node 181。并將潛艇表面劃分成35 854個三角形面元， 提取18 306個表面坐標， 如圖7所示。

圖7 ANSYS潛艇模型網格劃分Fig. 7 Submarine model meshing by ANSYS

考慮面元遮擋影響， 用改進板塊元方法對潛艇0～180°°范圍的目標強度進行預報仿真。考慮算法應用條件， 設置入射頻率為高頻（high frequency， HF）的仿真條件， 對潛艇的目標強度進行預報，結果如圖8所示。

圖8 高頻條件下潛艇目標強度Fig. 8 Submarine target strength at high frequency

從改進板塊元方法預報結果可以看出， 潛艇艇艏和艇艉目標強度明顯小于艇舯， 符合潛艇目標強度的蝶形圖特征［10］。預報結果能夠反映潛艇目標強度的一般規律。

圖9 高頻條件下2種方法對比Fig. 9 Comparison between two methods at high frequency

3 結束語

文中針對板塊元方法對潛艇目標強度預報時， 由于板塊元積分分母可能為零， 從而導致計算結果出現不穩定的情況， 將 Gordon積分方法用于潛艇目標強度的板塊元計算， 考慮并優化了模型中面元遮擋計算流程， 建立了基于改進板塊元方法的潛艇目標強度快速預報模型。仿真結果表明， 模型計算結果能夠反映潛艇目標強度的特征和規律， 可對潛艇目標強度進行快速預報， 從而改進模型的預報結果， 具有更好的穩定性。針對改進方法在計算耗時的不足， 下一步考慮將板塊元方法和改進板塊元方法相結合， 提出更適用于工程快速計算的優化算法。

［1］ 劉伯勝， 雷家煜. 水聲學原理［M］. 哈爾濱: 哈爾濱工程大學出版社， 2009.

［2］ 范軍. 水下復雜目標回聲特性研究［D］. 上海: 上海交通大學， 2001.

［3］ 馬忠誠， 范景理. 中國船舶科技報告: 水下模型散射聲場模型［R］. 大連: 第760研究所， 1995.

［4］ 范軍， 萬琳， 湯渭霖. 沉底水雷目標強度預報［J］. 聲學技術， 2001， 20（4）: 145-148. Fan Jun， Wan Lin， Tang Wei-lin. Prediction Target Strength of Mines Lying on the Seabed［J］. Technical Acoustics，2001， 20（4）: 145-148.

［5］ 王汝夯， 李冰. 基于板塊元法的水面艦船遠場聲目標強度預報方法［J］. 聲學技術， 2015， 34（1）: 6-10. Wang Ru-hang， Li Bin. A Method for Predicting Far-field Target Strength of Surface Ship Based on Planar Elements Theory［J］. Techcical Acoustics， 2015， 34（1）: 6-10.

［6］ Houshmand B， Chew W C， LEE S W. Fourier Transform of Linear Applications in Electromagnetic［J］. IEEE Transactions on Antennas and Propagation， 1991， 39（2）：252-254.

［7］ Gordon W B. Far Field Approximations of the Kirchhoff-Helmholtz Representations of Scattered fields［J］. IEEE Transactions on Antennas and Propagation， 1975， 23（5）: 590-592.

［8］ 劉成元， 張明敏， 程廣利， 等. 一種改進的板塊元目標回聲計算方法［J］. 海軍工程大學學報， 2008， 20（1）: 25-27. Liu Cheng-yuan， Zhang Ming-min， Cheng Guang-li， et al. Improved Planar Element Method for Computing Target Echo［J］. Journal of Naval University of Engineering， 2008，20（1）: 25-27.

［9］ 姬金祖， 劉戰合. 基于面元分組的電磁遮擋算法及其優化［J］. 北京航空航天大學學報， 2009， 35（4）: 453-456. Ji Jin-zu， Liu Zhan-he. Electrom Agnetic Occultation Algorithm Based on Facets Group and Optimization［J］. Joumal of Beijing University of Aeronautics and Astronautics， 2009， 35（4）: 453-456.

［10］ 何心怡， 張鵬， 王光宇， 等. 大俯仰角時的潛艇目標強度研究［J］. 艦船科學技術，2007， 29（4）: 37-39. He Xin-yi， Zhang Peng， Wang Guang-yu， et al. Target Strength Research of Submarine under Large Pitching Angle［J］. Ship Science and Technology， 2007， 29（4）: 37- 39.

（責任編輯: 楊力軍）

Simulation of Submarine Target Strength Forecast Based on Improved Planar Element Method

SUN Nai-wei1，2， LI Jian-chen1， WAN Ya-min1， ZHAO Gang1， Lü Wei1， FAN Ruo-nan1
（1. The 705 Research Institute， China Shipbuilding Industry Corporation， Xi′an 710077， China； 2. Science and Technology on Underwater Information and Control Laboratory， Xi′an 710077， China）

Aiming at the instability in computing submarine target strength with the planar element method in case the integral denominator is zero， Gordon integral algorithm is used for submarine target strength forecast， and the process of planar occlusion judgment for complicated target is simplified. Firstly， a rigid sphere is modeled. Then， a 3D model of Benchmark submarine and a forecast model of submarine target strength based on the improved planar element method are built to forecast submarine target strength. Simulation results show that the proposed forecast model gains more stable computation， and the obtained data can well reflect the features of submarine target strength.

submarine； planar element； Gordon integral； occlusion judgment； target strength forecast

U674.941； TB566

A

1673-1948（2016）04-0254-06

10.11993/j.issn.1673-1948.2016.04.003

2016-06-01；

2016-07-05.

孫乃葳（1993-）， 男， 在讀碩士， 研究方向為自導仿真技術.