石 权，龚雅萍，孙 峰，王帅军，张彦鑫，吴 姝，张 迪，陈 辉

1. 浙江海洋大学 a. 海洋工程装备学院，b. 船舶与海运学院，浙江 舟山 316022

2. 浙江舟山海王星蓝海开发有限公司，浙江 舟山 316100

随着现代海洋牧场建设的快速推进，发展海水养殖业已成为践行大食物观的重要举措[1-2]。中国是世界第一大水产养殖国，2022 年中国海水养殖产量达2 275.7 万t，其中深水网箱养殖产量约39.3 万t，相较2011 年，增长率高达600%[3-5]。网箱养殖已逐渐成为我国海水鱼类养殖的主要方式[6-7]。自1998 年海南省首次从挪威引进柔性抗风浪网箱开始，到2002 年我国第一套拥有自主知识产权的高密度聚乙烯 (High density polyethylene,HDPE) 升降式深水网箱成功下水，以及2018 年“德海1 号”成功抵御17 级超强台风“山竹”，再到集水产养殖、智慧渔业、休闲渔业、科技研发、科普教育等功能于一体的新型海洋牧场“耕海1 号”成功运营，提高网箱装备化、智能化、生产集约化、管理信息化、环境可持续化，推进网箱养殖朝着深远海方向发展是未来海洋渔业转型升级的必然趋势[8-12]。

在深水网箱结构受力计算方面，Yao 等[13]基于纳维-斯托克斯方程，提出了一种用于求解流固耦合的混合体积方法，并分析了刚性和柔性网箱在不同流速下承受海流负载的情况。Li 等[14-16]基于莫里森方程和浪的衍射和辐射理论，耦合浮架单元梁结构并模拟网衣桁架结构，分析了网箱在海流和波浪下的受力情况。Park 等[17]通过3 种不同直径的HDPE 管材的浸没程度，推导出网箱浮架的阻力系数，计算了网箱浮架在流动载荷作用下的变形和应力分布。杨新华等[18]基于绕射理论和Morison 方程，计算了波浪力并建立运动方程，对圆柱形浮式网箱进行受力分析并提出相应的数值计算方法。郭帅等[19]基于势流理论和非线性运动控制方程，建立了一种重力式抗风浪网箱的数值模型，对比分析了其在不同台风工况下与传统网格式网箱系泊缆张力、浮架运动和网箱变形的情况。

为缓解近岸养殖压力、拓宽深海水域养殖，基于中国海域养殖条件，本研究设计了一种可以应对恶劣海况条件的新型复式抗风浪深海养殖网箱，与传统HDPE 重力式网箱抗风浪结构不同，所设计的网箱创新采用了超静定抗风浪结构，以双浮筒结构取代传统的浮圈和浮架结构，并对极端海洋条件下网箱可能承受的风、浪、流等环境载荷进行计算和分析，为发展大型抗风浪养殖网箱提供了设计参考。

1 新型复式抗风浪网箱系统结构设计

1.1 主体结构设计

所设计的新型复式抗风浪养殖网箱，包括浮筒、立体浮框结构、系泊系统、网衣系统 (图1)。与传统单中心的养殖区网箱不同，网箱养殖区域分布在一体式立体浮框两侧，通过安装在养殖区域上方的浮筒提供浮力。浮筒与立体浮框结构之间通过连接结构形成刚性连接，当网箱受到风浪冲击时，两个浮筒随着海浪波动起伏，刚性立体浮框和浮筒相连保持浮筒在海上只做上下运动，使网箱整体在浮动中平衡稳定，减少风浪对网箱的影响。对称式双网箱的设计可以有效减少因刚性连接而发生的中拱和中垂现象，防止网箱框架受损破裂。

图1 新型复式网箱设计图Fig. 1 Design drawing of new type of composite aquaculture net cage

1.2 浮筒设计

国内常见的船型、组合式网架型网箱均以浮架提供浮力[20-22]，而对依靠浮筒为网箱整体提供浮力的设计并不多见，特别是双浮筒一体式抗风浪网箱，在国内尚未见类似的报道。图2 为筒形结构的浮筒，浮筒不仅为网箱提供漂浮海面的浮力，也是网箱各功能舱室设置的重要地点，从上往下依次为生活层(保证养殖工作人员的饮食起居)、机电设备区和辅助舱室(燃油舱、淡水舱)、饲料储存舱、投喂设备舱和压排载舱。浮筒采用双层壳体结构设计：外层设有金属防腐涂层，以防海水、海风侵蚀，并设有防撞层，以免遭受船舶撞击而发生破损，导致整个网箱出现结构性毁坏。底部压载区用于盛放海水，以达到调节网箱重心位置，使其始终在浮心以下，防止其倾覆。浮筒可根据作业海域的风浪状况调节距离位置，大于一个海域内出现概率多的波长，通过调整间距和压载重量分布，使得纵摇固有周期增加，避开出现概率最高的波浪周期。

图2 浮筒功能舱设计图Fig. 2 Design of buoy function tank

1.3 立体浮框结构设计

如图3 所示，立体浮框结构由底框、立柱和顶框组成，顶框与底框之间通过立柱连接后形成一体。底框结构由纵梁、横梁及斜撑组成 (图4)，3 根横梁和2 根纵梁一起围合出两片养殖区，双养殖区设计既可实现多元化养殖，亦可避免网衣破损后，单一养殖区鱼群全部逃逸的情况。同时，一体式立体浮框结构区别于需要连接结构的传统组合网箱，后者的连接点受浪、流载荷较大，极易发生断裂和破坏[23]，而一体式立体浮框结构则通过整体结构设计，使每根柱面受力而非点位受力，将受到的浪、流载荷分散开，从而实现抗风浪功能。

图3 立体浮框结构侧视图Fig. 3 Side view of three-dimensional floating frame structure

图4 立体浮框结构的仰视图Fig. 4 Upward view of three-dimensional floating frame structure

底框的材质选用HDPE 浮管，HDPE 材料不仅可以提供一定的浮力，有一定的抗拉性，且经济适用性好。顶框在结构设计上和底框相似，由纵梁I、横梁I 组成 (图5)。顶框主要用于连接定位，顶框钢管应选择高密度、高强度的合金钢管，特别需要在应力集中部位设置加强结构，整体增加中央钢管的抗剪、抗弯强度[24]。

图5 立体浮框结构的俯视图Fig. 5 Top view of three-dimensional floating frame structure

立体浮框结构是网箱的主体结构，也是网箱设计的关键，其结构强度关系到网箱的结构安全。所设计的立体浮框结构与浮筒组合共同为网箱提供浮力，基础框架和配重一起产生向下的拉力，保证网箱在海上垂直方向上的受力平衡。

1.4 网衣设计

网衣的主要作用是固定养殖空间，限制养殖鱼群的活动范围，防止鱼群逃逸，同时也可阻挡鲨鱼等捕食。网衣系统由顶网、侧网、底网、基础框架和配重5 个部分组成 (图6)。网衣系统根据立体浮框形状来设计，由侧网和底网共同围合成养殖区域，侧网和底网利用基础框架相连接，通过自身和配重的重力在水下张开成固定形状。

图6 网衣及配重设计图Fig. 6 Design drawing of net and additional weight

网衣材料采用超高分子量聚乙烯纤维 (UHMWPE) 材料，其分子量是普通聚乙烯(PE)的30～50 倍，断裂强度约为尼龙复丝的4 倍[25]。网衣网目形状为菱形，大小为35 mm×35 mm。由于网衣在自身重力和其他外界施压力的情况下极易变形，造成网衣破损，导致养殖鱼群逃逸，并且在流速较大的情况下，如在深远海，配重系统对网箱网衣变形程度的影响较大[26]，故需要增加配重系统，本网箱设计配重为2 342 kg。

1.5 网箱抗风浪结构设计

本网箱抗风浪功能实现的关键在于垂荡体和立体浮框结构连接形成的超静定结构。超静定结构是指通过浮筒下部锤挡板和空心柱组合形成的垂荡体(图7) 和立体浮框结构组成“十字形”稳定刚性连接，为网箱垂荡及纵向摇摆运动提供较大阻尼，实现网箱运动幅度与波浪运动形成相对“静止”和稳定。

图7 垂荡体和立体浮框的连接结构Fig. 7 Connection structure of dangling body and three-dimensional floating frame

垂荡体由浮筒、空心柱和锤挡板组成，上半部分柱状浮筒和空心柱相连部分设计了锥形波浪缓冲区，上半部分空心柱与缓冲、减摇作用的锤挡板直接连接，得益于锤挡板提供的较大垂荡以及纵向摇摆阻尼，因此在海浪较大的情况下，网箱运动幅度也能有效降低。

本网箱纵梁II 结构过渡区主要有两种作用：连接浮筒，给浮筒提供一个斜向支撑力，使其在风浪条件下做上下运动；做波浪缓冲带，能够使波浪作用在锥面上缓慢爬升释放波浪能量。与直壁结构不同的是，波浪在直壁结构上爬升得较高 (图8)，这样会增加网箱承受的波浪载荷，而斜向处理可以引导波浪沿斜面破碎，从而显著降低网箱承受的波浪载荷。

图8 直壁结构波浪爬升情况图Fig. 8 Wave run-up of straight wall structure

2 抗风浪网箱受海洋载荷的受力计算

计算网箱环境载荷，首先要考虑网箱放置海域的自然条件状况，如海风、海浪、海流等，特别是风暴潮、台风等极端天气，易破坏网箱结构，造成养殖鱼群逃逸、网箱框架损坏等严重后果。因此，设计抗风浪网箱时，对网箱结构所承受的环境载荷进行计算十分必要。

2.1 海流环境载荷

重力式抗风浪网箱的网架结构并不会因海流作用而产生惯性力，而其网衣受水下海流阻力的影响，因此分析海流载荷时，主要是对网衣形变进行计算。鉴于网衣表面结构不规律，故不考虑水通过网目时的绕流作用，而对网线进行分割计算，进而对水阻力进行评估分析。首先明确养殖海域的水流条件：海水密度ρ取1.025×103kg·m—3，水流速度v取1 m·s—1，网衣的受流面积A为287 m2。

根据莫里森公式、流体动力的一般公式，网衣在水流中的阻力为：

式中：FD为阻力 (N)；CD为水阻力系数。

网衣受阻计算的关键在于网衣水阻力系数的选取，本文参考李玉成和桂福坤[27]针对网衣垂直状态下水阻力系数与雷诺数关系的实验结果，在网衣形变程度一定的情况下，总体水动力估算采用正常状态下水阻力系数，网衣水阻系数与雷诺系数的关系可简化成：

式中：黏性系数μ=1.01×10−3，网线直径d1=0.001 8 m，则Re=1 826.73，可得水阻力系数CD=0.908 8。网箱受到的水阻力：FD=131 108 (N)=131.108 (kN)。

2.2 波浪环境载荷

网箱结构受波浪载荷分为网衣和浮架结构两部分，本文主要针对网箱主体浮架结构，利用莫里森公式分别计算网箱横、纵向框架柱的受力。选取蒲福氏风级表中8 级大风下的波浪参数：波高H=8 m，水深D=50 m，周期T=8 s，波长L=99.5 m。故波数k=2π·L—1，波频ω=2π·T—1，阻力系数CD=1.2，惯性力系数Cm=2，水流速度v=1 m·s—1。

2.2.1 纵向框架柱的受力计算

考虑到单个纵向框架柱所受波浪载荷包括速度力和惯性力，先进行计算单位长度下的速度力极值：

积分得：

根据现有标准和设计尺寸：附加质量系数γ=1.2，管直径d2=0.6 m，则单个垂直框架圆柱上受到的速度力极值为：

单位长度下的惯性力极值：

积分得：

则单个垂直框架圆柱上受到惯性力极值为：

则单个垂直框架圆柱所受的波浪力：

所有垂直框架圆柱所受波浪力：

2.2.2 横向框架柱的受力计算

采用修改后的莫里森公式计算水平圆管上的波浪力，即作用在单位长度上的水平圆管上的力：

引用线性波浪理论，时间t处位置 (x,z) 的波浪方程：

速度势：

水平速度：

所有水平框架圆柱所受波浪力：

则浮架总受力为：

2.3 海风环境载荷

因所设计的抗风浪网箱主体结构均在水下，仅浮筒和穿梭走廊受到的海上风载荷，考虑到我国频发台风海域的台风风级以及目前材料、制造水平，极限风速取12 级台风最大值37 m·s—1。

对作用在网箱表面的风载荷采用伯努利方程，计算海风作用在浮筒上的动压：

式中：Wp为风压 (k N·m—2)；ru为空气密度(kg·m—3)。

将空气密度 (ru) 替换成重度 (r)，其表达公式为：r=ru·g，代入式 (14)，得：

设在海面标准大气压101 325 Pa、温度15 ℃的标准状态下，空气重度，重力加速度g取9.8 m·s—2，代入式 (15) 得：

由此可得网箱所受风载荷为：

式中：A表示垂直风速方向作用下网箱的受风面积(m2)，而网箱阻力风向主要来源于水平风向，故得出网箱结构的受风面积A=309.3 (m2)。

将深水网箱作业环境极值风速与网箱受风面积代入式 (17)，得到风荷载F=264.645 (kN)。

2.4 环境载荷计算对比分析

根据计算结果并对比传统HDPE 重力式网箱，以黄小华等[28]针对沿海地区常用的40～80 m周长的深水网箱计算结果作为参考标准，在波高5 m、周期8 s、流速0.75 m·s−1条件下，周长80 m网箱所受的波浪力最大为85 kN，远小于新型复式抗风浪网箱设计承受的波浪载荷。海流载荷计算选取聂雪军等[29]和黄小华等[30]对重力式网箱在水流作用下受力计算的结果作为参考标准，二者在计算时均只考虑水流对网衣的作用，在流速1 m·s−1、网衣受流向面积210 m2的条件下，网箱受流阻力均值为109.4 kN，与新型复式抗风浪网箱设计承受载荷相近。风载荷的计算则以吴常文等[31]在舟山海区进行的HDPE 深水网箱抗风浪试验结果作为参考标准，当风速大于35 m·s−1时，周长48 m 网箱所受风载荷最大为21.923 kN，远小于新型复式抗风浪网箱设计承受的风载荷。基于以上对比，可见本研究所设计的网箱相较于传统HDPE 网箱能承受更大的风浪载荷。

3 小结

针对传统网箱结构简单，无法承受海洋中风、浪、流等环境载荷造成的压力冲击情况[32-33]， 本研究基于海上实际作业工况，设计了一款面向深海养殖的复式抗风浪网箱，分别采用莫里森公式和伯努利方程对网箱承受风、浪、流等环境载荷进行受力计算并与传统HDPE 网箱进行对比，总结如下：

1) 新型复式抗风浪网箱创新采用双浮筒的设计，能为网箱提供更大的浮力，并与立体浮框结构连接组成超静定抗风浪结构。该设计区别于传统HDPE 网箱抗风浪结构，网箱上部承受波浪载荷的变化不明显，当网箱上部载荷变化时，整体不会出现较大倾斜而影响功能性的使用，故网箱上部位置可以储备大量物资，增大了甲板使用面积。

2) 新型复式抗风浪网箱可通过调节两浮筒的间距来改变浮筒下方垂荡体的距离，使其大于工作海域极限波长，通过垂荡体的垂荡和减摇作用，给网箱提供较大的阻尼，从而实现控制网箱运动幅度，保证网箱的结构安全。

3) 立体浮框结构是承受波浪载荷的主要载体，其特征尺度简化为横杆和纵杆，受到速度力和惯性力的叠加，在海上8 级大风条件下所承受的极值负载约为344.82 kN。网箱因采用双浮筒结构，受风面积大，且选取12 级大风极值风速，其极值负载约为264.65 kN。网箱承受海流载荷的计算是通过网衣水下受流阻力计算而得，其极值负载约为131.11 kN。在近似计算条件下，所设计的网箱承受海洋环境载荷均高于传统HDPE 网箱，为后续研究开发与生产提供了参考。