破冰船冲撞式破冰结构强度数值仿真方法研究

冯卫永，李 辉

(1. 中国人民解放军91404部队，河北 秦皇岛 066000；2. 哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨 150001)

破冰船冲撞式破冰结构强度数值仿真方法研究

冯卫永1，李 辉2

(1. 中国人民解放军91404部队，河北 秦皇岛 066000；2. 哈尔滨工程大学 船舶工程学院，黑龙江 哈尔滨 150001)

本文基于瞬态动力学分析法对破冰船冲撞式破冰过程进行有限元数值仿真研究。针对该方法中的船冰撞击分析、海冰浮力模拟等一些关键问题进行研究，根据上述方法对某破冰船冲撞式破冰过程进行数值模拟，并对应力计算结果进行校核。总之，本文建立的一套评价破冰船冲撞式破冰过程船体应力方法，对于中国老龄破冰船的强度评估具有重要意义。

破冰船；结构强度；冲击；数值模拟

0 引 言

我国黄渤海水域冬季冰期期间，破冰船需为在此海域通行作业的船舶或海洋装备开辟航道以保障其作业安全性。我国现役的破冰船大多存在服役年限超出、船体结构腐蚀严重等问题。对于老龄破冰船，其推进系统的衰减会导致破冰船的推力降低，使得破冰船的整体破冰能力下降，因此对于较薄厚度的冰层需要采用冲撞式进行破冰。船冰撞击过程中会对船体结构的安全性带来威胁，因此校核冲撞式破冰过程中船体的结构强度有重要的意义。目前国内外学者对破冰船连续式破冰数值仿真方法进行了研究，并取得了较大的进展[1 – 5]，但对于冲撞式破冰的研究相对较少。

本文采用有限元数值仿真的方式对破冰船在撞击海冰过程中的船体结构安全性进行研究。首先探讨模拟海冰脆性破坏方法，建立船体钢材弹塑性模型，提出海冰与船体接触模拟方法和海冰所受浮力模拟方法，在此基础上建立破冰船碰撞式破冰的仿真模型，并以某型破冰船进行结果验证。

1 材料失效模型的建立

1.1 海冰脆性破坏的模拟

海冰的脆性断裂强度需要通过单向拉伸和压缩实验测得[6]。海冰材料的参数如表1所示。海冰是拉压异性材料，其压缩强度要远远大于拉伸强度[7]，因此两者的破坏机理不同。

通常为了模拟材料的破坏需要通过实验来确定材料的多轴应力强度[8]，因此根据材料力学中的第一强度理论来表示海冰破坏的临界状态就可以通过有限元来模拟海冰的失效。总之，对于海冰的失效模拟可采用各向异性的弹性材料模拟海冰失效前的行为，通过定义海冰的最大失效压力作为海冰的临界失效准则。

表1 海冰材料模型的参数Tab.1 Material parameters of ice model

1.2 船体钢材的弹塑性模型

结构材料的本构关系是碰撞问题中的一项重要研究内容。为了更真实地反映材料的特性，采用线性强化弹塑性模型[9]。材料的屈服应力为：

式中：σ0为板材屈服应力，T型材取3.55×108MPa，球扁钢取2.35×108MPa；E=2.06×105MPa为弹性模量；Eh=1.18×103MPa为硬化模量。

船冰碰撞是一个动态响应过程，材料的动力特性不能忽略。由于船用低碳钢的塑性对应变率高度敏感，其屈服应力与拉伸强度极限随应变率的增加而增加，故本文在材料模型中引入应变率敏感性的影响。在材料应变率敏感性的众多本构方程中，采用与实验数据符合较好的Cowper-Symonds方程：

式中：σ0′为塑性应变率ε时的动屈服应力；σ0为相应的静屈服应力；D=40.4和q=5为钢材的常数。

一般来说，材料的失效过程较为复杂，本文采用MSC.DYTRAN程序中给出的最大塑性应变失效准则。即当结构单元的等效塑性应变达到定义的单元最大塑性失效应变时单元失效，失效单元将不再参与计算。

2 撞击过程中海冰受力分析

2.1 碰撞模型仿真

通过有限元接触面模型可以模拟刚体与刚体、变形体与刚体、变形体与变形体间的相互作用。为模拟撞击须定义相互接触的2个面及其方向，对于破冰过程的模拟定义船首为主面、海冰为从面。此外，碰撞过程中摩擦力也需要考虑。

本文采用罚函数法计算主从接触面间的的接触力。假定接触力与穿透量gn成正比，穿透量与罚函数Kn相关，则法向的接触力fn可由下式得到：

在使用罚函数方法计算接触力时，接触面间不允许存在初始穿透，但在计算过程中允许存在穿透，这样通过穿透量与Kn可以求出接触力fn，在采用罚函数法计算接触力时首先需要选取合适的罚参数值来保证数值求解过程的收敛。在船冰碰撞的计算过程中需要根据多次的试算来确定合适的罚参数Kn。

为计及摩擦力影响，位移偏量u由弹性位移ue和滑动位移us：

在考虑摩擦影响的接触模型中，摩擦力的与法向力fn成正比：

式中：μ为摩擦系数；fn可由式（3）求得。当相邻2点发生接触时，摩擦力可由下式计算：

式中：Kt为切向刚度；Fμ为临界摩擦力。

2.2 海冰浮力的模拟

在破冰过程中，海冰由于船首的作用会产生垂向运动，所以需考虑海冰浮力的变化。采用非线性弹簧元模型来计算浮力变化，弹簧单元对海冰的作用模型如图1所示。

连续的非线性弹簧如图2(a)所示，离散弹簧元如图2(b)所示。图中Kf为刚度系数；ρi为海冰密度；ρw为海水密度；h为冰厚；w为垂向位移；Aieff为有效面积。

3 瞬态动力学方程的求解

在瞬态动力学分析中，一般采用中心差分法进行时间推进：

在非线性显式求解方法中，时间步长由最小的网格单元决定，以保持计算的稳定：

4 算例分析

为了验证上述所建立的冲撞式破冰模型的合理性，采用某破冰船进行验证。

4.1 有限元模型建立

本算例中所采用破冰船型的轮廓图如图3所示，相应的主尺度及相关参数如表2所示。计算中的海冰的相关参数如表3所示。

在有限元模型建立中，破冰船的结构用板单元建立，海冰用六面体实体单元建立，周围的海水用欧拉网格模拟，海冰所受浮力用非线性弹簧单元模拟。所建立的破冰船、海水、海冰组成的整体模型如图4所示。

表2 破冰船主尺度及相关参数Tab.2 Main particulars of the icebreaker

表3 海冰相关参数Tab.3 Parameters of ice

本文所选取的算例工况中，碰撞过程中破冰船以10 kn的初始速度破冰，设计破冰厚度为0.6 m。破冰船在冲撞式破冰过程中，船体首部会与海冰发生剧烈碰撞[10]，并产生瞬时的高应力。为了准确模拟这一现象，需将船首结构网格进行细化以保证求解精度，如图5所示，而对于其他位置采用刚体模拟即可。

4.2 船冰碰撞过程受力分析

通过软件模拟得到的破冰船以10 kn的速度撞击冰域后的纵向碰撞力和前进速度的变化曲线分别如图6和图7所示，撞击某一瞬时的海冰破碎形状和应力分布如图8所示。

由模拟结果可看出，船冰在撞击过程中，冲击力随着破冰船进入海冰而逐渐增加，纵向碰撞力的最大值为7.6 MN。碰撞过程中碰撞力会出现瞬间的变化，这是由于海冰单元出现失效导致的碰撞力瞬间下降，直到破冰船再次与海冰单元接触撞击力变大，如此循环反复。在整个碰撞过程中，船体航行方向的加速度逐渐增大，直到速度降为0。

4.3 应力结果分析

在进行强度校核时，需把船体结构分成若干个构件分别进行强度校核，主要包括甲板、船体外板、横舱壁、横框架、纵舱壁、纵桁等。撞击过程中，各部分构件发生的最大应力情况如表4所示。

表4 各构件强度校核结果Tab.4 Strength assessment results

从以上的应力计算结果可知，当破冰船以10 kn 速度冲撞时，部分船体构件出现结构强度问题。船体最大应力发生在外壳板上，其应力分布云图如图9所示，最大应力356 MPa 已超出材料的屈服极限355 MPa。根据其他工况的计算还可以发现，当破冰船以6 kn 和8 kn 初速度撞击海冰时，船体受到的最大应力分别为253 MPa和216 MPa，此时船体结构较为安全。但当初始速度为10 kn 时，船体结构已发生危险。由此可知，出与船体结构安全因素考虑，冲撞式破冰船的破冰速度选取宜适中。

5 结 语

本文针对冲撞式破冰船提出了结构强度数值分析方法。通过非线性有限元方法进行船冰碰撞的数值模拟，本文对一些关键技术进行研究，其中包括船冰相互接触碰撞算法、接触面间摩擦系数在碰撞力计算中的影响、非线性弹簧元的使用以及显式求解方法等。

根据某破冰船冲撞式破冰过程的算例分析表明，本文所建立的模拟方法可行。并分析得到了撞击过程中船体的应力分布及危险构件强度校核。

[ 1 ]BIAO S, RISKA K, MOAN T. A numerical method for the prediction of ship performance in level ice [J]. Cold Regions Science and Technology, 2011(60): 177–188.

[ 2 ]BIAO S, RISKA K, MOAN T. Numerical simulation of local ice loads in uniform and randomly varying ice conditions[J]. Cold Regions Science and Technology, 2011, 65(2): 145–159.

[ 3 ]BIAO S, RISKA K, MOAN T. Numerical study of ice-induced loads on ship hulls[J]. Marine Structures, 2011, 4(65): 1–21.

[ 4 ]XIANG T. A six-degrees-of-freedom numerical model for level ice-ship interaction[J]. Cold Regions Science and Technology, 2013, 92(8): 1–16.

[ 5 ]XIANG T, RISKA K, MOAN T. Effect of dynamic bending of level ice on ship’s continuous-mode ice breaking[J]. Cold Regions Science and Technology, 2014, 106(10): 82–95.

[ 6 ]钢质船舶入级与建造规范[S]. 中国船级社(CCS), 2012

[ 7 ]沈梧, 林树枝. 渤海东部一年冰单轴压缩强度的应变速率敏感性[J]. 大连工学院学报, 1984, 23(4): 45–50.

[ 8 ]孟广琳, 张明远, 李志军. 渤海平整冰单轴抗压强度的研究[J]. 冰川冻土, 1987, 9(4): 329–338.

[ 9 ]季顺迎, 王安良, 苏洁. 环渤海海冰弯曲强度的试验测试及特性分析[J]. 水科学进展, 2011, 22(2): 266–272.

[10]季顺迎, 岳前进, 工程海冰数值模型及应用[M]. 北京: 科学出版社, 2011.

Numerical simulation of structural strength of icebreaker in the crash-breaking process

FENG Wei-yong1, LI Hui2
(1. No. 91404 Unit of PLA, Qinhuangdao 066000, China; 2. College of Shipbuilding Engineering , Harbin Engineering University, Harbin 150001, China)

Numerical simulation of collision of icebreaker with ice is investigated with FEM method in this paper. Key issues during the simulation are undertaken, such as impact procedure, calculation of ice buoyancy etc. Simulation of an icebreaker during collision with ice is conducted as an illustrative example and structural strength is accessed. In summarize, the method proposed in this paper is of great significance for strength evaluation of icebreakers.

icebreaker；structural；crash-breaking；numerical simulation

U674.77

A

1672 – 7619(2017)06 – 0028 – 05

10.3404/j.issn.1672 – 7619.2017.06.006

2016 – 12 – 05

冯卫永(1977 – )，男，工程师，主要从事船舶定型试验工作。