海洋用A3钢阴极保护电场的三维有限元建模.pdfVIP

海洋用A3钢阴极保护电场的三维有限元建模 孙虎X,o王在峰①② 孙立娟① (①中国科学院海洋研究所；②中国科学院研究生院) 摘要根据有限元法能够预测复杂体系中阴极保护电化学场的分布原理。为了解介质中 的LAPLACE方程，使用8节点立方体和合适的形状函数对构筑物表层和周围场域进行了离散。 通过对各单元的分析和整个求解区域的综合，建立r三维的有限元数学模型，根据此模型编写 了相应的C语言程序。通过对模型的计算得出了电位电流密度分布图，发现规律与经典电场原 理相符台。 关键词 有限元 三维模型 阴极保护 电化学场Seidel迭代 对于复杂的实海受阴极保护体，阴极保护电位电流的现场原位测鼍和监测比较困难，使用 有限元计算可以求出阴极保护体表层和周围场域中的电化学场分布。有限元法在设计阴擞保 护体系和目I极保护监控方面已体现出了其巨大优点。Raymond，Fu，邱枫等人都对限元法计算 阴极保护电化学场进行了研究，作者也曾建立了：维的有限元模型进行了一些实例分析。但 是以前的研究大都借助NASTRAN，MARC等力学、热学有限元程序来分析电化学场；由于计算 复杂，研究多以二维为主，并且对阴极表面非线性边界条件的处理也4i多见。作者建立了三维 模刊的有限元模型并剐一个模拟实例进行r计算，对非线性边界条件解决提出了一种新方法。 l三维数学模型的建立 阴极保护电场中不考虑时间的变化影响，设各向电导率为常数，金属表面和介质中的电位 分布可用IAPI^cE方程和边界条件表示： 域内 (嚣)2+(器)2+(器)2：o E=Eo 边界厂． 塑一Ⅱ盟 边界工12 an—o E为电位，d腐蚀介质电导率，i(E)金属的极化曲线作边界条件，n为边界内法线 方向。 用有限元法解上而偏微分方程需化成求其等价泛函的极小值： ，(E)：Ⅲ吉【(器)2+(等)2+(差)2]av+盯掣Eas ‘7 岛 =min E=Eo 孙虎元．研究员，博士，主要从事海洋腐蚀监测研究。E-mail：sun@msqdioac．cn。 223 为求厂(E)的极小值，将整个域划分成M个单元，一个顶点，厂(鳓为多元函数要使其达 到极小值需满足$互掣：o(z：j，2，…，¨在每个单元内构造如下插值函数： F， E| E=＼N．～j·一～。 M 巨，， 将泛函对(E。 簟… E。)分别求导 坌丛啦一盯『垫驾 挚dxdydz+ aE： 一圳【孤ax ■‘吖B 盯f垫蚋…丛1 警dxdYdz+ 上u L ay J ay ay 皿[嚣嚣 E弓M‘鹾弓M 萼dXdYdz E 同理掣…哿也可表示成上面的形式． 将各单元的分析结果扩展到整个区