潛艇電化學腐蝕電場特性研究?

劉春陽 汪小娜 王向軍

（海軍工程大學電氣工程學院 武漢 430033）

潛艇電化學腐蝕電場特性研究?

劉春陽 汪小娜 王向軍

（海軍工程大學電氣工程學院 武漢 430033）

初步分析潛艇在海水中腐蝕電場分布和衰減特性，為潛艇防腐和電場隱身提供理論和數據基礎?；陔姶艌龌痉匠虒С隹諝狻Ｋ畠蓪幽Ｐ椭校瑵撏Цg電場電位控制方程和邊界條件，利用有限元法對深海中潛艇腐蝕電場進行建模仿真。仿真結果表明：潛艇腐蝕電場在潛艇附近海水中各分量都能達到mV∕m量級，在潛艇螺旋槳附近幅值最大且變化劇烈，潛艇腐蝕電場衰減速度隨海水深度增加而減緩。潛艇腐蝕電場具有明顯的分布特性和衰減特性。

海水；潛艇；電化學；腐蝕電場；有限元

1 引言

海水中含有大量可溶性鹽，具有強烈的腐蝕性［1］。而潛艇的結構部件主要是金屬材料，不同材料的金屬在海水中產生不同的腐蝕電位，不同電位的金屬發生電連接將會產生電流，從而產生電化學腐蝕電場。腐蝕電流是艇殼鋼板腐蝕的主要原因，也是艦艇電場的主要來源［2～3］。潛艇腐蝕不僅使潛艇結構強度下降，縮短潛艇使用壽命，其產生的腐蝕電場在海水中傳播還會增加潛艇被探測、發現的可能性，成為水中兵器的信號源，嚴重影響潛艇隱身性能［3～5］。

采用電偶極子對潛艇腐蝕電場建模，簡單實用，不需要大量計算［6］，但是無法直接反映潛艇結構參數對電場分布的影響。利用有限元分析軟件對潛艇腐蝕電場建模，求解精度高，且可以直觀反映潛艇周圍海水中電場分布，可為潛艇電場防腐和隱身提供理論基礎。

2 腐蝕電場控制方程和邊界條件

當金屬或合金與電解質溶液接觸時，通常會因為原電池作用發生電化學腐蝕。潛艇的艇體和螺旋槳由不同合金材料制成，潛艇不同腐蝕部位與海水電解液接觸界面上的雙電層電荷密度不同，從而產生不同的電極電位。鋼質艇體和銅質螺旋槳由大軸連接，構成螺旋槳→連接軸→內電阻→艇體→海水→螺旋槳回路［7］。

當潛艇運動狀態和周圍海洋環境無劇烈變化時，潛艇電化學腐蝕速度幾乎不變，可認為電化學產生的腐蝕電流是恒定的。恒定電流將會產生恒定電場，海水可近似為各項同性導電媒質，海水中任意一點的電流密度矢量J和電場強度E滿足本構關系［7～8］：

恒定電場與靜電場一樣旋度為零，可用一個標量位函數的負梯度表示恒定電場的大?。?/p>

只要求出電勢φ在空間中的分布，就可由式（1）求解海水空間電場強度。海水中電位分布滿足拉普拉斯方程：

因此電化學腐蝕電場的計算和求解可以歸結為求解電位控制方程（2）。滿足拉普拉斯方程的解有無限個，因此必須給出電位滿足的邊界方程，才能求出海水中電位的唯一解，進而求出海水中電場分布。

恒定電場與靜電場一樣，電場強度切向分量和電流密度的法向分量在邊界面上連續。用電位函數描述邊界條件則為

深海區可以忽略海床影響，求解區域等效為空氣—海水兩層模型，在海水和空氣界面，空氣可視為是理想電介質，其電導率σ0=0，因此在空氣海水—分界面上法向電流為0，即海水中的電流沒有從分界面處流入空氣，電流只在海水內部流動。海水-空氣界面上電位函數滿足自然邊界條件：

因此求解海水中電場分布時無需對空氣部分建模。在發生電化學腐蝕的艇體和螺旋槳表面其電流密度和電勢關系滿足極化曲線，可由實驗測出。在離潛艇無窮遠處的海水區域電位為零，在實際計算中可取足夠遠為邊界，減少計算量［9～11］。

3 潛艇腐蝕電場有限元建模仿真

理論上根據潛艇周圍海水中電位控制方程和電位滿足的邊界條件可以求出潛艇周圍海水中的電位分布，進一步利用電磁場基本方程求出海水中電場強度分布。但是計算過程非常復雜，而且不能直觀得到海水中電場分布和潛艇腐蝕的關系。利用有限元仿真軟件對潛艇電化學腐蝕產生的恒流電場進行建模仿真，可以求解復雜潛艇模型的電場分布，而且求解精度高，可以得到海水中電場分布圖，直接反應腐蝕電場與艦船結構及參數變化的關系［12～13］。

本文利用Ansoft有限元仿真軟件建立潛艇仿真模型，為簡化分析過程，減少計算量，對潛艇模型適當簡化［9，11］：

1）用多個規則幾何體連接近似表示艇體光滑過度表面，艇體和螺旋槳之間由艦軸連接，艇體和螺旋槳都視為等勢體。

2）指揮臺和尾舵簡化為規則幾何體，且不考慮潛艇內部結構和電路系統對腐蝕電場的影響。

3）艇殼表面的防護涂層看成是有限導電薄膜，螺旋槳和艇體之間連接的主軸簡化為一個規則的圓柱體，單位電導率為40。

4）理論上無窮遠處海水電勢為0，為了減少計算量，取一個足夠大長方體作為求解域。

5）基本的電磁參數參考海水中的標準值，由于潛艇模型較為復雜，不考慮潛艇表面金屬的極化作用。

建立直角坐標系：以艇艏端點為坐標原點，X軸正方向為指向左舷，由艇艏指向艇尾為Y軸正方向，Z軸正方向垂直向上。潛艇仿真模型全長80 m，型寬8 m，潛艇仿真模型如圖1所示。

圖1 潛艇仿真模型

潛艇周圍海水等效為一個足夠大的長方體。在海水—空氣界面滿足自然邊界件，無須對其強加邊界條件，也不用對海面以上的空氣建立模型。潛艇外包海水X方向由300 m～-300 m，Y方向-500 m～500 m，Z負方向設為-250 m，Z正方向為潛艇下潛深度，海水外包模型如圖2所示。

圖2 海水外包模型

4 潛艇腐蝕電場特性分析

假設海水深度250m，潛艇下潛20m，潛艇艏部防腐涂層破損，其余部位防護涂層良好，即只有艇艏發生腐蝕。艇艏腐蝕電極電位為-0.64 V，銅質螺旋槳腐蝕電位-0.32 V。

應用Ansoft恒流電場求解器對潛艇電場進行求解，可以得到整個海水空間中腐蝕恒流電場分布。圖3為潛艇下方6 m與海平面平行的平面電場強度幅值分布云圖。

圖3 潛艇下方6 m與海平面平行的平面各方向電場分布圖

由圖3（a）可知，潛艇下方海水中電場X分量關于ZOY平面反對稱分布，在潛艇下方附近幅值較大，遠離潛艇區域幅值小。螺旋槳左右附近出現正負峰，螺旋槳附近峰值可達到0.84 mV∕m，艇體附近也有較大電場值，約為 0.44 mV∕m。由圖 3（b）可知，電場Y分量關于ZOY平面對稱分布，在螺旋槳前后附近出現正負峰值，正峰可達0.73 mV∕m，負峰可達到-1.68 mV∕m。在艇艏下方附近也有一個較小的正峰，幅值可達0.15 mV∕m。由圖3（c）可知電場Z分量關于ZOY平面對稱分布，有明顯的負峰，且負峰出現在螺旋槳后方，可達2.3 mV∕m，螺旋槳和艇體下方附近電場變化劇烈，遠離潛艇下方海水區域變化平緩。

圖4為距潛艇左舷6 m，與艏艉線平行不同深度求解路徑上電場各分量曲線圖。圖中實線為潛艇下方6 m曲線圖，短劃線為潛艇下方26 m曲線圖，長劃線為潛艇下方46 m曲線圖。x=100 m為艇艏頂端正下方，x=180 m為螺旋槳正下方。由圖4可知，電場各分量在不同深度變化趨勢相同，在螺旋槳下發附近出現峰值，且變化劇烈。隨著深度增加，潛艇腐蝕電場各分量衰減速度減慢，且X分量衰減最快，Z分量在潛艇下方附近下降較慢。深度為46m時，在潛艇下方X分量幾乎為0。

圖4 腐蝕電場不同深度時的分布曲線

5 結語

潛艇結構復雜、材料性質變化多樣，且不同海域海洋環境差異較大，因此很難用解析法準確求解潛艇電化學腐蝕電場在海水中分布。合理簡化潛艇模型和海水情況，對海水-空氣兩層模型中潛艇電化學腐蝕電場進行有限元仿真分析能對潛艇腐蝕電場進行初步分析。仿真結果表明，潛艇腐蝕電場具有明顯的分布和衰減特性，因此對潛艇腐蝕電場特性進行研究能為潛艇防腐和水下武器提供一定理論依據。

［1］江炎蘭，陳菊娜，吳世永.艦船的腐蝕與涂層保護技術［J］. 腐蝕與防護，2012，33（2）：139-143.

［2］鄭軍林，陳新剛，鄭春軍，等.艦船電場隱身技術［J］.中國艦船研究，2006，4（1）：48-51.

［3］牟蘭.國外艦船電場特性研究及其在水雷戰上的應用［J］. 艦船科學技術，2012，34（9）：138-142.

［4］DAGMAR M，PRAHA K.Continuous extendibility of solutions of the third problem for the laplace equation［J］.Czechoslovak Mathematical，2003，35（128）：669-688.

［5］ZOLOTAREVSKII Y M.Methods of measuring the lowfrequency electric and magnetic fields of ships［J］.Measurement Techniques，2005，48（11）：1140-1144.

［6］刁愛民，楊慶超，王杏青.海水中電偶極子電場分布有限元分析［J］. 四川兵工學報，2013，34（1）：136-138.

［7］王向軍，嵇斗，劉德紅等.艦船電場防護原理［M］.武漢：海軍工程大學，2016：72-80.

［8］何小祥，丁衛平，劉建霞.工程電磁場［M］.北京：電子工業出版社，2011：76-78.

［9］卞強，張民，柳懿等.一種基于ANSYS的艦船靜電場分析方法［J］. 海軍工程大學學報，2010，22（6）：65-69.

［10］岳瑞永，臧燕華，呂俊軍.基于邊界元理論模型的艦船腐蝕相關電場衰減規律［J］.艦船科學技術，2012，34（1）：11-14.

［11］曹寓，嵇斗，朱武兵.艦船水下腐蝕靜電場有限元仿真分析［J］. 艦船科學技術，2015，33（7）：69-72.

［12］孫建紅，鄭煒，王曉鵬.水面艦艇船體防腐和陰極保護的優化設計方法［J］. 中國艦船研究，2007，2（4）：60-64.

［13］汪小娜，肖昌漢，王向軍等.運動潛艇的感應電場分布［J］. 海軍工程大學學報，2012，24（4）：7-10.

Research on Electric Field Produced by Electrochemical Corrosion of Submarine

LIU ChunyangWANG XiaonaWANG Xiangjun
（College of Electrical Engineering ，Naval University of Engineering，Wuhan 430033）

A preliminary analysis of the submarine in seawater corrosion electric field distribution and attenuation characteristics provides theoretical and data basis for submarines anticorrosive and electric fields stealth.Electric field controllable equation and boundary condition in air-sea two layer structures are analyzed，based on the basic equation of electromagnetic.Corrosion related electric field of the submarine is modeled by finite element method.The simulation result shows that the electric field distribution of submarine has obvious symmetry，and the components are of magnitude of millivolts per meter.Corrosion electric field of submarine has obvious distribution and attenuation characteristics.

seawater，submarine，electrochemical，corrosion electric field，finite element

U674.76

10.3969∕j.issn.1672-9730.2017.10.004

Class Number U674.76

2017年4月7日，

2017年5月26日

劉春陽，男，碩士研究生，研究方向：電磁環境與防護。汪小娜，女，副教授，碩士生導師，研究方向：艦船電場防護。王向軍，男，教授，博士生導師，研究方向：艦船電場防護。