潜艇电化学腐蚀电场特性研究∗

刘春阳 汪小娜 王向军

（海军工程大学电气工程学院 武汉 430033）

潜艇电化学腐蚀电场特性研究∗

刘春阳 汪小娜 王向军

（海军工程大学电气工程学院 武汉 430033）

初步分析潜艇在海水中腐蚀电场分布和衰减特性，为潜艇防腐和电场隐身提供理论和数据基础。基于电磁场基本方程导出空气—海水两层模型中，潜艇腐蚀电场电位控制方程和边界条件，利用有限元法对深海中潜艇腐蚀电场进行建模仿真。仿真结果表明：潜艇腐蚀电场在潜艇附近海水中各分量都能达到mV∕m量级，在潜艇螺旋桨附近幅值最大且变化剧烈，潜艇腐蚀电场衰减速度随海水深度增加而减缓。潜艇腐蚀电场具有明显的分布特性和衰减特性。

海水；潜艇；电化学；腐蚀电场；有限元

1 引言

海水中含有大量可溶性盐，具有强烈的腐蚀性［1］。而潜艇的结构部件主要是金属材料，不同材料的金属在海水中产生不同的腐蚀电位，不同电位的金属发生电连接将会产生电流，从而产生电化学腐蚀电场。腐蚀电流是艇壳钢板腐蚀的主要原因，也是舰艇电场的主要来源［2～3］。潜艇腐蚀不仅使潜艇结构强度下降，缩短潜艇使用寿命，其产生的腐蚀电场在海水中传播还会增加潜艇被探测、发现的可能性，成为水中兵器的信号源，严重影响潜艇隐身性能［3～5］。

采用电偶极子对潜艇腐蚀电场建模，简单实用，不需要大量计算［6］，但是无法直接反映潜艇结构参数对电场分布的影响。利用有限元分析软件对潜艇腐蚀电场建模，求解精度高，且可以直观反映潜艇周围海水中电场分布，可为潜艇电场防腐和隐身提供理论基础。

2 腐蚀电场控制方程和边界条件

当金属或合金与电解质溶液接触时，通常会因为原电池作用发生电化学腐蚀。潜艇的艇体和螺旋桨由不同合金材料制成，潜艇不同腐蚀部位与海水电解液接触界面上的双电层电荷密度不同，从而产生不同的电极电位。钢质艇体和铜质螺旋桨由大轴连接，构成螺旋桨→连接轴→内电阻→艇体→海水→螺旋桨回路［7］。

当潜艇运动状态和周围海洋环境无剧烈变化时，潜艇电化学腐蚀速度几乎不变，可认为电化学产生的腐蚀电流是恒定的。恒定电流将会产生恒定电场，海水可近似为各项同性导电媒质，海水中任意一点的电流密度矢量J和电场强度E满足本构关系［7～8］：

恒定电场与静电场一样旋度为零，可用一个标量位函数的负梯度表示恒定电场的大小：

只要求出电势φ在空间中的分布，就可由式（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日

刘春阳，男，硕士研究生，研究方向：电磁环境与防护。汪小娜，女，副教授，硕士生导师，研究方向：舰船电场防护。王向军，男，教授，博士生导师，研究方向：舰船电场防护。