于 博 文，张 维 英，陈 静，于 欣，于 洋，周 俊 秋

(大连海洋大学 航海与船舶工程学院，辽宁 大连 116023 )

0 引 言

据中国渔业互保协会统计，1994～2015年的渔船事故中，全损(无法修复)渔船2 996艘，部分(可修复)损失渔船76 454艘．事故原因方面，由人员操作不当造成的碰撞事故占比44.81%，触损事故占比39.12%，触礁事故占比49.23%，搁浅事故占比36.04%[1]．所以探讨渔船在破舱条件下的稳性，尤其是非对称进水过程的六自由度变化，可以判断渔船的分舱合理性和给出破舱应急操作的建议，以达到避免或减少渔船事故损失的目的，对渔船的安全性具有实用价值[2]．

不断发生的海难事件促使人们不断发展船舶破舱稳性的研究，在模型试验之外对于破舱进水的研究方法可以分为两种．准静态法：将整体分成破舱管涌、破损舱室内的水运动以及波浪运动下的船舶运动[3-9]．在这些研究中，基于伯努利方程水动力模型被用来确定通过开口的流速．假定涌入舱内的水是瞬间沉降到一个水位[3，5-7]，并作为块状物在特定路径表面上自由移动[8-9]，或者晃荡并且遵循浅水方程[4]．势流理论法：应用切片理论[3，4，6]或面元法[5，7-9]计算船-波相互作用动力学．计算机硬件的飞速发展，增强了数值模拟用于解决非线性问题的能力．破舱进水过程中水的涌入、冲击、飞溅等都是典型的非线性状况，通过数值模拟方法解决非线性问题的应用日趋成熟，对基于计算流体动力学方法的高保真度模拟的需求也在增加．

刘强等[10]基于Fluent软件用Navier-Stokes(N-S)方程……