STAR-CCM+系泊问题：漂浮式海洋牧场养殖装置系泊系统设计

1.1 时域计算方法

对于研究浮体在波浪下的运动问题，首先需要求解流场速度势[10]。通过线型叠加入射势、辐射势及绕射势，以此来表示浮体周围流场的总速度势：

式中：Φ₁(x,y,z,t)为入射波速度势；Φ_R(x,y,z,t)为辐射速度势；Φ_D(x,y,z,t)为绕射速度势；(x,y,z)为流场中的位置坐标；t为时间。

结合运动学边界条件、动力学条件、远场边界条件以及物体表面边界条件，最终得到总速度势定解条件，见式（2），

其中：n为物面法向量；i为虚数单位；g为重力加速度；k为波数，满足k=ω²/g；v_n为浮体湿表面S₀的法向速度；ω为波浪频率。

式（2）为速度势的定解条件，入射波的速度势已知后，即可求得辐射速度势和绕射速度势。将浮体所受波浪力分成3种力，即入射、绕射和辐射，而入射波浪力与绕射波浪力构成波浪激励力。通过所得辐射势和绕射势，求解浮体的波浪力，进而建立浮体的频域运动方程，见式（3）。

式（3）中：m为浮体质量；λ_a为额外阻尼系数；μij为附加质量；λ_ij为辐射阻尼；K为静水回复系数；K_a为额外刚度系数；x_i为浮体的位移；F为波浪激励力。

之后将浮体的频域运动方程转换为时域运动方程。当浮体进行六自由度微幅运动时，可以把运动看成一系列脉冲运动。将浮体周围的速度势进行叠加求解。浮体运动下，流体辐射势可以表示为

式中：x_j(t) 为浮体在j方向运动的速度或者是角速度；ψ_j为浮体单位脉冲运动产生的速度势，且在j方向上；χ_j为浮体做单位脉冲运动后t时刻的速度势。

通过伯努利方程求得浮体表面动压力并进行湿表面积分，得到浮体辐射水动力和力矩：

式中：m_ij为附加质量矩阵；K_ij为延迟函数；X_j(t)为浮体在j方向运动的加速度或者是角加速度。进而推导出时域运动方程表达式为

式中：X_i(t)为第i个自由度上的波浪激励力；C_ij为浮体静水回复力矩阵；M_ij惯性矩阵。通过对式（6）换元得到：

在求解时域运动方程前课通过面元法，获得浮体的频域水动力系数。在通过频域水动力系数求解时域运动方程的水动力系数。同时波浪激励力也需要进行频域到时域的转换，这种方法称为间接时域法。在建立浮体的时域运动方程后，利用数值方法4阶龙格-库塔求解方程，最终得到浮体各自由度时历运动数据。

1.2 计算模型

本文中的漂浮式海洋牧场养殖装置的母型是由中国船舶集团有限公司建造的“海洋渔场1号”养殖网箱，整个平台基础结构式为半潜式，主要分为四大构件：上层平台、立柱、网箱结构和下浮体。表1为海洋牧场主尺度；

图1为漂浮式海洋牧场养殖装置的总体结构图；图2为海洋牧场水动力计算模型，其中，立柱和下浮体基于面元法计算，而横撑由于属于小尺度构建，采用莫里森单元进行计算。

2.1 环境参数

假定本文海洋牧场所处位置为黄海，该海域接近中国重要港口城市，具有运输方便的优点。综合水文环境资料，作业工况和生存工况环境参数见表2。

2.2风力机载荷

本文的风力机使用沈阳工业大学开发的 200 kW 风力机，不同风速下风力机所受推力（见图 3）和侧向力（见图 4）通过查阅文献所得[11]。根据所得数据将风力机在相应风速下所得载荷按照定常力和力矩的形式加入数值模型中。

2.3 水轮机载荷

水轮机选择翼型为NACA0018、3叶片的 5kW 潮流能水轮机。为了确定水轮机载荷，通过STAR-CCM+ 数值模拟得到本文水轮机侧向力（见图 5）和推力曲线（见图 6）。根据所得数据将水轮机在相应流速下所得载荷按照定常力和力矩的形式加入数值模型中。

2.4 系泊设计方案

整个海洋牧场的系泊系统设计为悬链线形式。整个系泊系统的设计不仅要考虑浮体的运动，还需要考虑系泊材料的属性等性能。本文将对比纯锚链方案和锚链-钢缆-锚链组合式系泊方案，系泊材料参数见表3。

根据海洋牧场形式，浮体的固定系泊点设置在结构物浮箱底部，每个系泊点连接两条锚链并与海底锚点连接。锚链对称分布，并且两两平行。海域水深为100m。设海洋牧场水平面中心为坐标原点，系泊线布置方式及系泊点数据具体坐标见图7和表4。

2.5 系泊设计要求

系泊系统的设计不仅要满足浮体所规定的位移要求，还需要满足系泊线的强度要求。在作业工况和生存工况下需要保证浮体不会发生倾覆，系泊系统不会断裂等情况。中国船社级的《海上移动平台入级规范》对锚链不同状态下张力的安全系数要求见表5。

由于海洋牧场上并没有长期居住的作业人员，也没有对运动响应有严格要求的装置，因此运动响应只需要保证在作业工况下装置的正常工作，生存工况下系泊达到规定的安全系数即可。而对于合理的位移偏移量，也没有具体的设计规范，因此通过查数据以及经验，要求作业工况下，平台的纵荡和横荡运动限定在水深的10%左右。由于风力机具有自动偏航功能，横摇和艏摇对发电效率影响小，仅考虑纵摇影响。限定作业工况下纵摇角度小于5°。

3.1 作业工况结果分析

3.1.1 纯锚链方案

设计流向与风向为0°，设置不同的浪向角，研究海洋牧场在纯锚链系泊方案下3h内的六自由度运动情况以及系泊系统张力变化，统计结果见表6。

观察表6数据可知：作业工况下，平台在浪向角为0°时的纵荡最大，为5.68m，符合浮体位移要求。而平台在其他自由度的运动数值较小，这是由于浪向角在0°时，浪向与风流同向，对平台纵荡方向的影响较大。

作业工况下，纯锚链系泊在各浪向角下的六自由度最大摇荡的时历曲线见图8，系泊系统的最大张力时历变化曲线见图9。由图8和图9可知：纯锚链系泊方案中，系泊锚链张力最大为665kN，最小安全系数为7.750，锚泊系统符合规范要求。在作业工况下，海洋牧场上的风力机与水轮机受到风速和流速的作用，在风速流速方向产生较大位移。

3.1.2 锚链-钢缆-锚链方案

在使用锚链-钢缆-锚链进行组合系泊时，对海洋牧场运动响应及系泊受力进行分析，具体结果见表7。

由表7可知：在各个浪向角下，纵荡均有较大的运动，最大达到了11.78m，略微超过平台位移设计要求。对比2种系泊方式在作业工况下的数据，见表8。

对比作业工况下的2种系泊方案的数据可知：组合系泊方案相比纯锚链系泊方案，在纵荡方向的位移更大，同时系泊的安全系数也更小，而在其他方向上的位移相差不大。综合对比后发现，纯锚链系泊方案更加符合设计要求。

3.2 生存工况结果分析

3.2.1 纯锚链方案

生存工况下，海上牧场搭载的风力机和水轮机停止工作，只考虑海况对平台和系泊系统的影响。海上牧场在纯锚链系泊方案下，六自由度运动最大值统计和锚链最大张力计算结果见表9。

由表9可知，生存工况下，平台除纵荡以外的自由度变化较小，纵荡方向的最大值为12.47m，横荡方向的最大值为3.1m，这是由于在纵荡方向平台受到波浪作用以及风力机水轮机受到风速流速作用共同影响；锚链最大张力为2037kN，安全系数为2.53，大于最小安全系数1.67，满足要求。

生存工况下纯锚链系泊方案在不同浪向时，平台六自由度最大摇荡时历曲线见图10。系泊系统发生最大张力的时历曲线见图11。

3.2.2 锚链-钢缆-锚链方案

在使用锚链-钢缆-锚链进行组合系泊时，对海洋牧场运动响应及系泊受力进行分析，结果见表10。

由表10可发现，在各个浪向角下，纵荡运动依然最大，最大达到了12.39m，而其他自由度的运动极值变化不大。对比两种系泊方式在设计工况下的数据，见表11。

由表10和表11数据可知，在生存工况下，2种系泊方案的系泊安全系数均达到规范要求。

从浮体位移的角度来看，组合式系泊方案下的海上牧场在各浪向角时的纵荡位移要略大于纯锚链系泊方案，其余自由度运动变化幅值上相差不大；而从系泊系统角度来看，组合式系泊方案的系泊张力安全系数要小于纯锚链方案。综合考虑，本文中的海上牧场建议使用纯锚链系泊方案。

3.3 系泊系统部分失效结果分析

在海上牧场在海上作业及生存时，可能会由于磨损或海底撞击等原因导致系泊系统部分失效，为了提高平台的生存性能，设计系泊系统在各浪向角时受力最大的系泊线断开，研究其他锚链是否可以支撑整个系统的生存。纯锚链系泊方案部分失效状态系泊张力结果见表12。

通过研究发现，在6号锚链破损断裂情况下，系泊系统锚链张力最大值为2478kN，最小安全系数为2.08，大于所需的安全系数，验证了系泊部分失效后，在作业工况和生存工况下，剩余的锚链系泊系统仍能够保证海洋牧场的安全。

本文以半潜式海上牧场为时域计算目标，对海洋结构物在波浪中的时域计算方法进行分析。在考虑风力机、网衣和水轮机载荷的情况下，结合本文中海上牧场的结构特点和有关文献研究，设计了六点式系泊方案，设计了纯锚链系泊方案和组合式系泊方案2种系泊形式并进行了性能对比分析。通过对比分析2种系泊方案，得到以下结论：

1）在作业工况和生存工况下，两种系泊方案均符合系泊系统设计的安全要求，但组合式系泊方案下的海上牧场在各浪向角时的纵荡位移要略大于纯锚链系泊方案，系泊张力安全系数要小于纯锚链方案。综合考虑，选择采取纯锚链作为海上牧场的系泊系统。

2）分析系泊系统部分失效状态下的性能，结果表明：受到最大张力的系泊线因意外断裂时，其余锚链在生存工况下安全系数仍大于所规定的系泊安全系数，说明该系泊方案设计具有足够的安全性。