基于邊界元的艦船水下靜態電場推算方法

何　芳，王向軍，高俊吉

(海軍工程大學電氣工程學院,湖北 武漢　430033)

?

基于邊界元的艦船水下靜態電場推算方法

何芳，王向軍，高俊吉

(海軍工程大學電氣工程學院,湖北 武漢430033)

0引言

在艦船靜態電場建模方面，目前主要有邊界元法(BEM)[1-2]、有限元法(FEM)[3-4]、等效電偶極子法[5-6]，艦船電場的深度換算是電場建模的重要研究內容。艦船靜態電場的深度換算目前國外未見有報道，國內可見報道很少，可查資料顯示只有海軍工程大學在該方面做了初步工作，主要采用點電源法[7-8]、電偶極子法[9-11]、平面迭代法[12]進行深度換算。文獻[8]采用點電源產生的靜電場模擬船舶，利用泰勒展開式的收斂半徑對場源的個數和位置進行了準確的計算，建立了基于點電源的靜態模型，再反演得到模擬源的模型，最后利用此模型計算目標深度上的靜電場。此方法對建模要求很高，模型的準確度對換算結果影響較大，計算過程較為復雜。文獻[9-11]主要是基于電偶極子進行無限大平面的深度換算，解決了電場某分量向下進行較大深度自換算和互換算的問題，文獻[12]利用微分原理，根據靜態電場各分量與其垂直分量的偏導數關系，建立各無限大平面之間的數學關系，采用迭代算法解決海水中靜態電場向較淺深度換算問題。但是基于無限大平面理論的深度換算在轉化為數值積分后，對換算深度、測量點距及測量平面尺寸都有要求，若換算深度小于測量點距或測量平面尺寸接近換算深度時都會產生較大誤差，同時在工程實際應用中，無法在無限大平面上測量，只能在有限平面內取值，造成了換算精度的下降，不利用指導工程實踐。本文針對此問題，提出了基于邊界元的艦船水下靜態電場推算方法。

1水下靜態電場的邊界元模型

建立水下靜態電場邊界元模型的思路是：先建立半無限域任意形狀的閉合曲面邊界面積分模型；然后通過對邊界面進行正方形面元剖分，將邊界積分方程離散為代數方程組，解出邊界源點近似解，并對邊界奇異9+++點進行數值處理；最后根據域內場點與邊界源點近似解的關系式，求解半無限域內的靜態電場分布。

1.1任意形狀的閉合曲面邊界面積分模型

艦船靜態電場是由海水中的穩恒電荷分布產生的，基于媒質內部均勻且各向同性的假設,穩恒電荷只分布在空氣和海水、海水和艦船以及海水和海床的分界面上, 海水內部無凈電荷存在，電荷面密度與分界面兩邊媒質的電容率有關。若在近場區取一任意形狀閉合曲面，海水無限遠處邊界面視為S∞,S0-S∞構成一無限域V,半無限域V的內邊界面S0的外法線方向指向無限域內部，如圖1所示。半無限域V內的靜態電場則滿足Laplace方程。

(1)

設邊界面S0上有一源點M0,M為無限域V中的任意一場點，r為源點到場點的距離。半無限域內場值與內邊界面S0上函數值關系，可以通過Green第三公式[13]推導得出：

(2)

圖1　艦船水下電場半無限區域平面圖Fig.1　Schematic of vessel electric field in the semiinfinite domain

1.2邊界積分方程離散及奇異點處理

設邊界面S0上的有限測量點為M0，根據Green第二公式[13]推導出邊界積分公式為：

2πu(M0)=

(3)

(4)

式中，i=1,2,…,k。ri,j為面積元dSi的中心與面積元dSj的中心之間的距離。

(5)

當i=j時 ，場點和源點在同一積分面元上，邊界面積分出現奇異點。用正方形小面元(面元邊長為a)，正方形面元坐標設置如圖2所示，計算得

(6)

圖2　正方形面元的坐標設置Fig.2　setting of coordinate in a element

用矩陣方程描述即為Ak×kU1×k=Bk×k,系數矩陣A，B均為k×k的滿秩方陣，U1×k即邊界面上的整個電場分布。為了便于說明，在笛卡兒坐標系下對代數方程組系數A，B進行確定。以平行于潛艇縱向方向由艦尾向艦首為X軸正方向，以平行于潛艇橫向方向由左舷向右舷為Y軸正方向，以垂直于水面指向地心為Z軸正方向。以潛艇為中心取一長方體表面作為邊界閉合面，六個面上的法線方向都指向長方體內部，如圖3所示。

圖3　長方體邊界面及坐標系Fig.3　Boundary surface and its coordinate system

ri,j為面積元dSi的中心與面積元dSj的中心之間的距離：

(7)

只要邊界面的形狀和具體剖分方法確定了，矩陣方程Ak×k=Bk×k的系數矩陣A，B就能惟一確定。解此方程組,可以直接計算出U1×k，得到邊界面上各面元中心處的電勢u，即邊界面上源點的近似解。

1.3基于邊界元法的水下靜態電場分布

將邊界面上各面元中心處的電勢代入式(2)，得到半無限場域內電場強度的函數解析解。

(8)

(9)

2基于邊界元的艦船水下靜態電場推算方法

為了接近實際情況(方便潛艇通過)，不考慮邊界面S1，S2兩個側面的電場分布，驗證采用長方體四個面近似閉合曲面的水下靜態電場邊界元模型是否準確可靠，采用MATLAB進行仿真。將一單位電偶極子放置于水深5 m處，設計邊界面為一長方體外表面，如圖3所示，長寬高分別為x∈(-6,6)m,y∈(-1.5,1.5)m,z∈(-1.5,1.5)m。忽略長方體S1，S2兩個側面，將其余四個側面每個面都剖分為648個正方形面元，共計剖分2 592個面元，每個面元的中心點電場強度作為整個面元電場強度的平均值。 應用式(9)，換算得到長方體表面積外的半無限區域內任意位置的電場三分量，以此作為水下靜態電場邊界元模型。

2.1利用邊界元模型外推直線處電場

利用水下靜態電場邊界元模型外推直線z=20 m,y=0,x∈(-50,50) m處的電場三分量，以及電偶極子在相同位置直接產生的電場三分量對照情況如圖4所示，與文獻[5]的電場分布解析表達式仿真結果(理論值)完全吻合。

2.2利用邊界元模型外推平面處電場

利用水下靜態電場邊界元模型外推平面z=20 m,y∈(-50,50) m,x∈(-50,50) m處的電場三分量，如圖5所示。仿真結果與文獻[9]采用無限大平面法換算得到的相同平面處電場三分量在數量級上是完全吻合的。

仿真結果證明了采用長方體四個面近似閉合曲面的邊界元模型是可行的，且驗證了基于該邊界元模型的水下靜態電場推算方法的正確性和可靠性。

圖4　電場三分量仿真推算值和理論值對比圖Fig.4　The results of the calculated and the theoretical for the components of electric field intensity

圖5　邊界元換算平面電場分布圖Fig.5　Result of BEM on plane

3實驗驗證

為了驗證基于邊界元水下靜態電場推算方法的有效性及換算精度，在實驗室采用潛艇模型進行實驗驗證。

3.1實驗條件和實驗步驟

實驗硬件設施包括無磁性實驗水池、實驗潛艇模型、電場三分量傳感器及矩形測量架。無磁性實驗水池，尺寸為8m×5m×1.5m,水深0.8m，將工業用鹽倒入池水中，充分溶解并混合均勻，模擬淡鹽水，測得所配制的海水電導率為3.96S/m。潛艇模型長1.5m,寬0.13m，高0.15m，按實艇縮小比例1∶100制得。艇模船殼鐵制，螺旋槳青銅制，船模的不同金屬結構在模擬淡鹽水中因電化學作用產生腐蝕電流，為防止船殼腐蝕，艇模上裝載了陰極保護系統(ICCP)，而防腐系統的開啟也會在模擬淡鹽水中產生電流。這些腐蝕和防腐電流在模擬海水(池水)中形成了靜電場。

矩形測量架寬0.4m,高0.4m。完全置于水中固定，測量架每條邊上均勻布放3個電場三分量傳感器。

實驗步驟如下：

1)將艇模移至距測量架較遠處，電場三分量傳感器采集背景電場。

2)以測量架為中心建立笛卡兒坐標系，潛艇沿X軸勻速穿過測量架，如圖6所示(潛艇由紙內向紙面運動)，電場三分量傳感器可以采集到艇模在各個位置時的電場值，即可得到一個除前后兩個端面外的矩形包絡面。用矩形包絡面上各測量點處的電場減去背景電場即可得到艇模在矩形包絡面上的電場分布。

圖6　潛艇與傳感器位置示意圖Fig.6　Schematic of position of submarine and sensor

3)根據所建立的靜態電場邊界元模型推算距潛艇左舷3m,水深0.2m,即直線 z=0.2m,y=3m,x∈(-8,8)m處的電場。

圖7　電場強度測量值與邊界元推算值的比較Fig.7　The comparision of electric field intensity betweenthe measured and the computed

3.2實測數據結果分析

根據等效偶極子理論[5-6]，實驗中船模可以看做電場源點，z=0.2m,y=3m,x∈(-8,8)m直線處各點作為場點。以船模中心為對稱軸，在偏離源點較近處電場模值出現大幅度衰減，偏離源點較遠處衰減變弱。在該實驗建立的笛卡兒坐標系中，船模中心坐標(0，0，0)，場點坐標(0，3，0.2)及(6，3，0.2)代入文獻[5]中靜電場分布解析公式驗證，可以發現在船長約4倍處，靜態電場強度模量衰減約75%，4倍以內衰減較快，4倍以外衰減變緩。該實驗實測數據與文獻[5]中直流電偶極子等效理論公式在相同源點與場點處計算結果一致，證明該實驗是有效的。

3.3邊界元模型換算數據分析

根據式(9)可知，電場強度與場點和源點之間的距離成反比。當場點與源點距離較近時，電場強度衰減越快；場點與源點距離較遠時，電場強度衰減越慢。定性上來看，靜態電場邊界元模型推算的電場強度模值是符合靜電場在海水中的衰減規律的。根據邊界元模型換算數據顯示，在船長約4倍處，靜態電場強度模量衰減約75%，與實驗實測數據衰減率接近，可以把船長4倍處看做電場強度模值衰減特征點。

4結論

本文提出了基于邊界元的艦船水下靜態電場推算方法。該方法不局限于無限大平面，可對半無限域任意形狀的閉合曲面建立邊界面積分模型，并進一步采用長方體四個面近似閉合曲面建立水下靜態電場邊界元模型，基于該邊界元模型進行水下靜態電場的推算。實驗驗證和結果分析證明了基于邊界元的水下靜態電場推算方法的正確性和可靠性，其外推電場值與實測電場值相對誤差為5%左右，能達到較高的換算精度，相比于無限大平面法的深度換算，更適宜指導工程實踐。

基于邊界元的艦船水下靜態電場推算方法可以在以下方面做進一步改進：1)基于實際測量的原因，本文忽略了包絡面的前、后兩個端面，來近似封閉面，由于前后兩端面面積較小，仿真結果表明誤差不大，但仍會產生一定的計算誤差。2)在對邊界面進行剖分時，將每個面元中心點上的值作為該面元上的平均值，這會產生離散誤差，因此小面元剖分得越細，面元越多，面元的面積越小，邊界元推算的精度越高。3)當場點趨于源點時，邊界面積分涉及到對奇異點的數值處理，邊界面積分精度很大程度上依賴于奇異積分的計算精度，可以在對奇異積分處理方面對文中的邊界元模型進行改進。

參考文獻:

[1]ZAMANING.Boundaryelementsimulationofthecathodicprotectionsysteminaprototypeship[J].AppliedMathematicsandComputation, 1988, 26:119-123.

[2]IWATAM,HUANGY,FUJIMOTOY.ApplicationofBEMtodesignoftheimpressedcurrentcathodicprotectionsystemforshiphull[J].JoftheSocietyofNavalArchitectsofJapan, 1992, 171: 377-380.

[3]DOIGP,FLFWITTPEJ.Afinitedifferencenumericalanalysisofgalvaniccorrosion[J].Electrochem,1979,126(12):2057-2063.

[4]KASPERRG,APRILMG.Electrogalvanicfiniteelementanalysisofpartiallyprotectedmarinestructures[J].Corrosion, 1983, 39(5): 181-183.

[5]劉勝道.艦船水下電場的測試技術與電偶極子模型研究[D].武漢：海軍工程大學, 2002：43-45.

[6] 劉文寶，王向軍，嵇斗.基于電偶極子模型的艦船靜電場深度換算[J].空軍雷達學院學報，2010(6)：435-438.

[7] 劉忠樂，龔沈光.海水中穩恒電流電場的點電極計算模型[J].海軍工程大學學報，2004，16(1)：35-39.

[8] 姜潤翔，胡英娣，龔沈光.基于點電源的船舶靜態電場深度換算方法研究[J].電波科學學報，2014，29(4)：685-693.

[9] 陳聰，李定國，龔沈光.基于拉氏方程的艦船靜態電場深度換算[J].電子學報，2010，38(9)：2025-2030.

[10]李定國，陳聰，劉華波.艦船水下靜態電場深度換算的格林函數法[J].海軍工程大學學報，2012，24(3)：1-5.

[11] 陳聰，龔沈光，李定國.艦船靜態電場深度換算方法[J].哈爾濱工程大學學報，2009，30(6)：719-722.

[12]胡英娣，龔沈光，閆永貴.一種新的船舶靜態電場深度換算方法[J].海軍工程大學學報，2013，10(25)：17-20.

[13]姚端正，梁家寶.數學物理方法[M].武漢：武漢大學出版社，2010:250-252.

摘要:針對基于無限大平面理論的深度換算在轉化為數值積分后，對換算深度、測量點距及測量平面尺寸都有要求，若換算深度小于測量點距或測量平面尺寸接近換算深度時都會產生較大誤差等問題，提出了基于邊界元的艦船水下靜態電場推算方法。該方法不局限于無限大平面，可對半無限域任意形狀的閉合曲面建立邊界面積分模型，進一步采用長方體四個面近似閉合曲面建立水下靜態電場邊界元模型，并利用該邊界元模型進行了水下靜態電場的推算。實驗驗證和結果分析證明了此方法的有效性，能達到較高的換算精度。

關鍵詞:邊界元；靜態電場；半無限區域

Extrapolation of Vessel Static Electric Field Based on Boundary ElementHE Fang，WANG Xiangjun，GAO Junji

(Electrical Engineering College，Naval University of Engineering，Wuhan 430033，China)

Abstract:When the depth is less than?measuring distance, or the size of measuring plane is close to the depth, extrapolation based on the infinite plane will have a greater error. If extrapolation is based on the finite plane, a decrease in accuracy is caused. To solve this problem, a method of depth extrapolation of static electric field was proposed based on boundary element method in the semi infinite domain. The method was not limited to the infinite plane, but also could establish the boundary integral model for arbitrary closed surface in semi infinite domain. Static electric field boundary element model under water was built up based on four surface of rectangular approximated the closed surface, and the underwater static electric field could be calculated using this boundary element model. The submarine model experiment testified that this method was of higher precision and stronger operability.

Key words:boundary element; static electric field; semi infinite domain

中圖分類號:TM46

文獻標志碼:A

文章編號：1008-1194(2015)06-0055-06

作者簡介:何芳(1981—), 女,湖北荊門人，講師, 博士研究生，研究方向：電場防護與隱身技術。E-mail:wskay1103@163.com。

基金項目:國家自然科學基金項目資助(41476153)

*收稿日期:2015-03-27