0 引 言

系泊系统对于保证船舶、船员和码头的安全性而言至关重要，国内外针对系泊系统所建的物理模型和所做数值研究较多。物理模型试验成本高、通用性差；计算机数值仿真方法可弥补该不足，但需在计算前定义船舶水动力性能特征等较多的参数，不仅费时费力，而且会耗费较多计算资源，计算时间较长[1]。利用相关船级社和国际组织颁布的规范来计算船舶系泊力可较为合理地解决船舶系泊力计算的经济性、科学性和安全性问题。

石油公司国际海事论坛（OCIMF）已通过大量的试验和模拟总结出相应的系泊力计算经验公式，通过给出船型、球艏形状、液罐形状、装载工况和深度吃水比等综合因素来确定船舶的安全系泊方案，可较为合理地描述环境载荷，相对于舾装数计算，该方法能对系泊力进行较为准确的估算，模拟船舶的真实受力情况，准确性和合理性要高[2]。本文简要介绍 OCIMF系泊力计算的方法和原理，并通过软件编程实现计算机辅助系泊力计算，供国内同行参考和借鉴。

1 系泊力计算原理

1.1 标准系泊模型

码头系泊的模式需根据船舶和码头的具体情况来确定。船舶行业针对码头系泊设定了一个标准的系泊模型，并考虑标准的环境载荷工况[3]。通常，艏部和艉部布置有数根横缆（主要用来限制船舶的横向运动），舯部靠艏艉附近布置有数根倒缆（主要用来限制船舶的纵向运动）。图1为码头标准系泊模式示意。

图1 码头标准系泊模式示意

1.2 标准环境载荷

作用在系泊状态的船舶上的环境载荷主要有风、海流、潮汐、过往船舶的船行波、波浪及冰载荷等。通常情况下，系缆若能抵抗最大的风和海流作用力，则其他较缓和的外界作用力也被认为是可承受的[4]。因此，规范对环境载荷的考虑主要涉及风和流作用力。

考虑任意方向60kn的风，同时考虑：0°或180°方向3kn的海流；10°或170°方向2kn的海流；90°方向0.75kn的海流。

风速选取海平面以上10m高度处30s的平均风速。若所给的风速不是10m高度处的平均风速，则需根据规范给出的公式换算成10m高度处的平均风速。若计算过程中给出的海流流速不是平均流速，则需根据规范给出的公式换算成等效平均流速。

1.3 系泊力计算方法

1.3.1 力学模型

OCIMF系泊力计算所考虑的相关环境载荷主要是风载荷和流载荷。在力和力矩的全局坐标系下，船舶受到的环境载荷可分解成x轴的纵向力Fx和y轴的横向力Fy，并外加一个分解到z轴的艏摇力矩Mxy。通常情况下，为方便计算，将全局坐标系下y轴的横向力Fy和z轴的艏摇力矩Mxy分解成艏柱处的横向力FyF及艉柱处的横向力FyA（见图2）。

图2 系泊力计算力学模型示意

1.3.2 风载荷计算

影响风载荷计算的主要因素有船舶纵向受风面积、横向受风面积、船长、吃水、风向角及球艏形状等，风载荷具体计算过程为

1) 求解出纵向风载荷Fxw，横向风载荷Fyw及艏摇风载荷Mxyw，即

式(1)～式(3)中：Cxw为纵向风载荷曳力系数；Cyw为横向风载荷曳力系数；Cxyw为艏摇风载荷曳力系数；wρ为空气密度；Vw为风速；AL为纵向受风面积；AT为横向受风面积；LBP为垂线间长。

2) 将横向风载荷和艏摇风载荷分解成艏柱处的横向风载荷FyFw及艉柱处的横向风载荷FyAw，即

1.3.3 海流载荷计算

影响海流载荷计算的主要因素有船舶深度吃水比、船长、吃水、流向角及球艏形状等，海流载荷具体计算过程与风载荷计算过程类似，这里不再赘述。

2 系泊力计算软件开发

2.1 编程语言介绍

Visual 即VB.Net，属于Basic系语言，其以极具亲和力的英文单词为基础标识，与自然语言有极为相近的逻辑表达，是Basic系语言中最强大的编程语言之一[5]。

2.2 软件主要组成

根据OCIMF系泊力计算原理和流程，使用Visual 在Visual Studio 2008平台上开发OCIMF系泊力计算软件，软件程序计算流程见图3。软件分为曲线数据模块、环境参数模块、船型参数模块、系泊模型模块、计算应用模块、系统设置模块和软件帮助模块等7个模块，各模块的功能结构见图4。

图3 OCIMF系泊力计算软件程序计算流程

图4 OCIMF系泊力计算软件各模块功能结构

2.2.1 曲线数据模块

文章来源：《电脑编程技巧与维护》 网址: http://www.dnbcjqywh.cn/qikandaodu/2020/1118/636.html