摘要
针对当前渔业养殖监测装备续航能力低、成本高、难以普遍应用,提出并设计了一种渔业移动与定点监测自主水下机器人(Autonomous underwater vehicle,AUV)。该AUV采用三体模块化结构设计,3个Myring型舱体呈“品”字型分布以易于坐底稳定,并采用T形支架进行连接以减少阻力。考虑三体结构3个舱体间阻力干扰作用,基于计算流体力学方法进行不同舱体间距阻力计算和设计。以直航阻力最小为优化目标,建立3个舱体和T形支架的参数化模型,通过改变舱体艏部和艉部曲线方程与T字形支架横截面曲线方程的参数来实现对AUV外形的控制,进行对AUV外形优化设计。结果表明,T形支架相比于传统的三角形支架直航阻力值降低了26.28%,AUV外形优化后的直航阻力值降低了24.05%,提高了续航力。研究结果为渔业养殖监测领域提供了一种新装备,同时也为其他三体结构AUV优化设计提供参考。
我国是世界第一渔业养殖大国,目前正着力建设以生态健康养殖为主的现代设施渔业,对渔业养殖监测装备的机械化、自动化、智能化提出了更高的要
减少航行阻力,提高续航能力是渔业监测AUV的重要性能指标之
本文结合渔业养殖的低成本、低能耗、立体移动和定点监测需求,提出并设计一种三体模块化结构的渔业监测AUV,进而开展三体结构及其连接框架的水动力性能分析与外形优化设计。
水产动物对食物的需求不同,栖息活动主要分布在水体中下层,越接近水底,水产动物密度越大,大多数水产动物喜欢群居且喜欢安

图1 渔业监测AUV定点和运动作业模式
Fig.1 Fishery monitoring AUV fixed-point and moving operation mode
AUV常用多推进器组合实现悬停来进行定点作业,但推进器产生较大的噪声,容易引起水产动物应激反应,影响其生长。同时,AUV靠自身携带的电池供电,推进器持续工作耗能大,难以满足长期作业需求。为了保证在水底安静地定点观测生物的生长状况和对环境的监测,提出并设计一种易于稳定坐底的三体结构渔业监测AUV,即AUV由3个回转体形的舱体构成“品”字形结构,即1个上舱体和2个下舱体,如

图2 渔业监测AUV三体模块化结构
Fig.2 Fishery monitoring AUV three-body modular structure
渔业监测AUV的上舱体主要负责水面、水中和坐底观测、电控与通信、浮力调节等,由水面通信浮标模块、视频监测模块、照明模块、浮力调节模块、电控模块和上舱体艉部模块组成。传统的AUV大多采用水下惯性导航系统和水声通信进行导航定位和通信,成本高,传输速率低。为此,本渔业监测AUV设计了一种内置4G、LoRa 和 GNSS 天线的低成本水面通信浮标模块,实现AUV定位导航及与岸基实时通信等功能。为了减少推进器对渔业环境的影响和降低能耗,浮力调节模块主要是通过控制外置油囊注油调节排水体积,实现无动力下潜上浮垂直运动。另外,浮力调节模块还包括一个辅助垂直推进器,当油囊浮力调节机构产生的浮力不足以克服水底黏土的吸附力或浮力调节机构出现故障时,提供额外上浮动力。
为了保证渔业监测AUV正常工作和安全作业,通过压载重量及其分布的衡重设计,保证AUV在水中呈水平姿态,且浮心和重心之间垂直距离大,上轻下重,即保证足够的稳心高度,以保持 AUV 航行姿态平稳。为此,渔业监测AUV设计2个相同的下舱体,并将沉重的电池布放在下舱体。另外,下舱体还包括艏部模块和艉部推进器模块。
三体结构AUV水下航行过程中,阻力的影响因素较为复杂。不同于单体结构AUV,三体结构AUV的阻力分析,不仅需要研究舱体的外形,还需要研究各舱体间流场扰动影响所产生的航行阻力,进而设计合适的舱体间距;另外,三体结构AUV为保持结构稳定性,需要通过连接支架进行固定连接。不同结构的连接支架会对AUV的水动力性能产生影响,设计易于搭载仪器装置和低阻力的连接支架很重要。
采用CFD进行数值模拟时,复杂的模型结构使得网格数量和计算资源需求增加,从而影响网格质量和计算结果的精
由于3个舱体左右对称,为减少计算量,取一半模型建立半径为1.5Lc的圆柱体计算域,其中Lc为舱体整体纵向长度,如

图3 舱体间距水动力分析
Fig.3 Hydrodynamic analysis of cabin spacing
由于舱体间呈正三角形,间距相等。以间距w为自变量,控制其余变量不变,分析不同间距对直航阻力F的影响。为了提高计算效率,将w进行参数化设置,利用Isight实现舱体间距的变更、网格划分、数值模拟和数据读取,同时确保计算结果的一致性。考虑到本AUV在水下主要进行水平直航运动,则计算中攻角设为0,考虑到实际监测中对航速的需求,航速设为1 m/s,直航阻力F计算结果如
由
传统的三体结构AUV连接方式,即连接支架形状大多采用三角形支架,其抗风浪性能较好,但安装拆卸较麻烦,不利于升级维
为了分析连接支架的水动力性能,将T形支架与三角形支架对比。按照上文设计,舱体间距一定的情况下,支架连接在相同位置且横截面形状都为水动力性能较好的椭圆形,建立对应的简易模型,计算两种连接支架在航速分别为0.25、0.5、0.75和1 m/s时的直航阻力F。两种连接支架模型在不同航速下的直航阻力值如

图4 连接支架水动力分析
Fig.4 Hydrodynamic analysis of connecting bracket
由
由前述分析与设计,通过结构设计、电控系统设计、各模块组装和系统联调等,研制的渔业监测AUV样机,如

图5 渔业监测AUV样机
Fig.5 Fishery monitoring AUV prototype
以渔业监测AUV直航阻力F最小为优化目标,采用参数化建模方法,在Isight上通过试验设计生成多组样本点来构造近似模型,进而使用全局优化算法基于近似模型对AUV的外形进行优化设计。
本文渔业监测三体结构AUV相比于单体AUV,外形难以由明确的公式直接表达,考虑优化设计时需要计算不同参数下外形的直航阻力,普通建模方式需要反复建模,计算效率低。因此,为提高计算效率,对渔业监测AUV进行参数化建模。渔业监测AUV的外形主要包括3个回转体舱体和T字形连接支架,如

图6 AUV外形示意图
Fig.6 AUV shape diagram
a1为上舱体艏段长度;b1为上舱体平行中体长度;c1为上舱体为艉部长度;Lq为上下舱体端面距离;Lt为T字形支架长度;h为上下舱体距离;a2为下舱体艏段长度;b2为下舱体平行中体长度;c2为下舱体艉部长度;dt为支架厚度;at为支架半椭圆长度。
a1 is the length of the bow section of the upper hull;b1 is the length of the parallel center body of the upper compartment; c1is the length of the upper hull for the transom;Lq is the distance between the upper and lower hull end faces;Lt is the T-bracket length;h is the upper and lower hull distance;a2 is the length of the bow section of the lower hull;b2 is the length of the parallel center body of the lower hull;c2 is the lower hull transom length;dt is the thickness of the bracket;at is the half-ellipse length of the bracket.
渔业监测AUV舱体选用Myring型回转体,Myring型回转体包括艏部、平行中体和艉部3段曲线,其曲线方程表达式参见文献[
艏部曲线方程:
(1) |
平行中体方程:
(2) |
艉部曲线方程:
(3) |
式中:x为长轴上点到艏部顶点的距离;a为艏段长度;b为平行中体长度;c为艉部长度;d为平行中体直径;n为艏部锐度因子;θ为艉部离去角;L为舱体总长。
T形支架迎流面由其横截面椭圆曲线方程控制,椭圆曲线方程表达式为
(4) |
式中:为支架截面椭圆锐度因子。
考虑渔业监测AUV各模块的结构和所需搭载器件和仪器对内部空间的要求因素,本文渔业监测AUV首先设计并确定的几何尺寸如下
参数名称 Parameter name | 参数尺寸 Parameter size/mm |
---|---|
D | 130 |
dt | 20 |
b1 | 820 |
b2 | 620 |
w | 520 |
h | 450.33 |
Lq | 100 |
Lt | 540 |
同时由
(5) |
(6) |
结合公式(
上舱体:
(7) |
T字形支架:
(8) |
下舱体:
(9) |
根据AUV的曲线方程组公式(
设计变量 Design variable | 初始值 Initialization value | 取值范围 Value ranges |
---|---|---|
a1 | 97.5 mm | (65 mm,130 mm) |
c1 | 195 mm | (130 mm,220 mm) |
n1 | 2 | (1,3) |
θ1 | 20° | (10°,30°) |
n2 | 2 | (1,3) |
θ2 | 20° | (10°,30°) |
at | 20 mm | (10 mm,35 mm) |
nt | 2 | (1,3) |
在Solidworks中设置设计变量为全局变量参数,以渔业监测AUV外形曲线方程为驱动方程组,并在拟合曲线上建立图形约束,完成参数化模型的建立,即由全局变量参数值直接控制模型的变化。
为保证近似模型的拟合精度以及提高计算效率,需要安排合理的试验设计方法。理论上试验设计生成的样本点越多精度会越高,但样本点数量过大,计算量激增导致优化流程时间过长。考虑最优拉丁超立方设计(Optimal latin hypercube design,OLHD)生成的样本点在设计空间内具有更好的均衡性和填充性,能够使所有样本点均匀覆盖在可行域范围
为减小计算量并且缩短计算周期,采用二阶多项式响应面模型(Response surface model,RSM)构造近似模型, 二阶多项式响应面模型数学表达式
(10) |
式中:m0、ei、fii、gij分别为常数项、一次项、二次项和交叉项的待定系数;xi为设计变量;k为设计变量的个数。
利用二阶多项式回归技术对OLHD得到的80组样本点和响应阻力进行最小二乘拟合得到二阶多项式拟合模型,拟合曲线系数如
常数项系数 Constant term coefficient | m0 | |||||||
---|---|---|---|---|---|---|---|---|
31.0740 | ||||||||
一次项系数 Monomial coefficient | e1 | e2 | e3 | e4 | e5 | e6 | e7 | e8 |
-0.048 6 | -0.062 6 | -0.459 9 | -0.126 1 | -0.537 7 | -0.079 3 | -0.475 7 | -2.708 3 | |
二次项系数 Quadratic coefficient | f11 | f22 | f33 | f44 | f55 | f66 | f77 | f88 |
0.000 1 | 0.000 1 | 0.086 0 | 0.000 5 | 0.094 9 | 0.000 7 | 0.004 7 | 0.417 4 | |
交叉项系数 Cross term coefficient | g12 | g13 | g14 | g15 | g16 | g17 | g18 | |
7.453 4 | -0.001 2 | -0.000 1 | 0.001 2 | 0.000 2 | 0.000 2 | 0.002 1 | ||
g23 | g24 | g25 | g26 | g27 | g28 | g34 | ||
0.001 0 | 0.000 4 | -0.000 6 | -1.047 7 | 7.305 7 | -0.001 1 | 0.000 7 | ||
g35 | g36 | g37 | g38 | g45 | g46 | g47 | ||
-0.015 3 | 0.006 3 | -0.003 3 | 0.007 0 | 0.005 3 | 0.000 3 | 0.000 6 | ||
g48 | g56 | g57 | g58 | g67 | g68 | g78 | ||
0.005 9 | -0.003 2 | -0.001 4 | 0.094 1 | 0.000 5 | 0.002 0 | 0.017 9 |
由生成的二阶多项式响应面模型,重新选取20个样本点进行误差分析,得到阻力的复相关系数
自适应模拟退火算法(Adaptive security algorithm,ASA)具有良好的全局搜索能力,计算效率高,广泛应用于工程领
(11) |
在Isight中采用ASA,设定初温为1,最大迭代次数为10 000,当迭代到6 460步时得到优化结果Fmin=7.379 N。相比初始值计算得到的阻力,阻力降低了24.05%。
为了验证优化结果的准确性,将优化结果中的数据重新进行CFD仿真分析,验证结果如
参数名称 Parameter name | a1/mm | c1/mm | n1 | θ1/(°) | n2 | θ2/(°) | at/mm | nt/mm | F/N |
---|---|---|---|---|---|---|---|---|---|
初始Initial | 97.50 | 195.00 | 2.00 | 20 | 2.00 | 20.00 | 20 | 2.00 | 9.716 |
优化Optimized | 65.25 | 218.92 | 1.62 | 10 | 2.59 | 27.69 | 35 | 2.19 | 7.379 |
CFD | 65.25 | 218.92 | 1.62 | 10 | 2.59 | 27.69 | 35 | 2.19 | 7.435 |
优化前后,渔业监测AUV的数值仿真结果如

图7 优化前后仿真结果
Fig.7 Simulation results before and after optimization
本文结合渔业养殖监测领域的现状和需求,提出并设计研制了一台三体稳定结构的渔业监测AUV。AUV设有定点和运动两种作业模式,同时采用模块化结构设计方法,设计各个功能模块,其中浮力调节模块和水面通信浮标模块降低了能耗和成本。基于计算流体力学方法,分析了三体结构舱体间距的干扰作用,确定了舱体间距为520 mm,进而针对渔业监测AUV的三体结构连接提出并设计了一种T形支架,通过与传统的三角形支架对比分析,T形支架直航阻力相比于三角形支架降低了26.28%,表明T形支架水动力性能更好。为了进一步减小直航阻力,对AUV的3个舱体和T形支架进行了参数化建模,选取舱体艏部和艉部与T形支架的曲线方程参数为设计变量,在Isight中采用OLHD生成样本点,并通过二阶多项式响应面模型和ASA得到了优化后的3个舱体的艏艉部和T形支架外形,使相同工况下阻力降低了24.05%,提高了AUV的续航能力。
参考文献
毛晓雅, 祖爽. 现代设施农业如何走得稳?[N]. 农民日报, 2023-11-10(008). [百度学术]
MAO X Y, ZU S. How can modern facility agriculture walk steadily?[N]. Farmers' Daily, 2023-11-10(008). [百度学术]
吴诗昊, 周悦, 敖琪, 等. 渔业养殖监测用AUV及其控制系统[J]. 上海海洋大学学报, 2023, 32(5): 953-966. [百度学术]
WU S H, ZHOU Y, AO Q, et al. AUV for aquaculture monitoring and its control system[J]. Journal of Shanghai Ocean University, 2023, 32(5): 953-966. [百度学术]
BAO H M, ZHANG Y, SONG M P, et al. A review of underwater vehicle motion stability[J]. Ocean Engineering, 2023, 287: 115735. [百度学术]
康帅, 俞建成, 张进. 微小型自主水下机器人研究现状[J]. 机器人, 2023, 45(2): 218-237. [百度学术]
KANG S, YU J C, ZHANG J. Research status of micro autonomous underwater vehicle[J]. Robot, 2023, 45(2): 218-237. [百度学术]
许裕良, 杜江辉, 雷泽宇, 等. 水下机器人在渔业中的应用现状与关键技术综述[J]. 机器人, 2023, 45(1): 110-128. [百度学术]
XU Y L, DU J H, LEI Z Y, et al. Review: applications status and key technologies of underwater robots in fishery[J]. Robot, 2023, 45(1): 110-128. [百度学术]
天津瀚海蓝帆海洋科技有限公司. 智帆AUV-ZF-01[EB/OL]. (2024-03-28)[2024-01-15]. http://tjhhlf.com/h-col-111.html. [百度学术]
Tianjin Hanhai Lanfan Marine Technology Co., Ltd. Zhifan AUV-ZF-01[EB/OL]. (2024-03-28)[2024-01-15]. http://tjhhlf.com/h-col-111.html. [百度学术]
郑荣, 胡志强, 朱兴华, 等. 一种三体构型的长期定点观测型水下机器人: 中国, 105644743B[P]. 2017-06-16. [百度学术]
ZHENG R, HU Z Q, ZHU X H, et al. A three-body configuration of long-term fixed observation underwater robot: CN, 105644743B[P]. 2017-06-16. [百度学术]
黄海, 周浩, 许锦宇, 等. 一种三体流线型自主作业水下机器人平台: 中国, 106477008B[P]. 2018-07-31. [百度学术]
HUANG H, ZHOU H, XU J Y, et al. A three-body streamlined autonomous operational underwater robotic platform: CN, 106477008B[P]. 2018-07-31. [百度学术]
杨睿, 冯旭琛, 于辉, 等. 三体水下机器人: 中国, 114670996B[P]. 2022-08-23. [百度学术]
YANG R, FENG X C, YU H, et al. Three-body underwater robot: CN, 114670996B[P]. 2022-08-23. [百度学术]
赵月. 水下航行器阻力及流噪声特性数值预报[D]. 大连: 大连理工大学, 2021. [百度学术]
ZHAO Y. Numerical prediction of resistance and flow noise characteristics of underwater vehicles[D]. Dalian: Dalian University of Technology, 2021. [百度学术]
李宏源, 吕鹏宇, 杜增智, 等. 水下滑移边界减阻技术研究综述[J]. 舰船科学技术, 2022, 44(9): 1-6. [百度学术]
LI H Y, LYU P Y, DU Z Z, et al. A review of drag reduction technology for underwater slip boundary[J]. Ship Science and Technology, 2022, 44(9): 1-6. [百度学术]
张贝, 张代雨, 王志东, 等. 基于自由变形的自主水下航行器壳体优化设计方法研究[J]. 船舶力学, 2022, 26(9): 1315-1325. [百度学术]
ZHANG B, ZHANG D Y, WANG Z D, et al. A free-form deformation optimization design approach for AUV hull shape[J]. Journal of Ship Mechanics, 2022, 26(9): 1315-1325. [百度学术]
HU F, HUANG Y, XIE Z B, et al. Conceptual design of a long-range autonomous underwater vehicle based on multidisciplinary optimization framework[J]. Ocean Engineering, 2022, 248: 110684. [百度学术]
VARDHAN H, SZTIPANOVITS J. Search for universal minimum drag resistance underwater vehicle hull using CFD[C]//International Conference on Computational and Experimental Simulations in Engineering. Cham: Springer, 2024: 1297-1303. [百度学术]
MILLER L M, NJAKA T, BRIZZOLARA S, et al. Drag reduction and power optimization due to an innovative, toroidal hull form of an AUV[J]. IOP Conference Series: Materials Science and Engineering, 2023, 1288(1): 012039. [百度学术]
戴鹏. 基于参数化的三体组合式自主水下航行器多目标优化设计研究[D]. 广州: 华南理工大学, 2020. [百度学术]
DAI P. Research on multi-objective optimal design of three-body combined autonomous underwater vehicle based on parameterization[D]. Guangzhou: South China University of Technology, 2020. [百度学术]
李为, 荣宽, 覃丽蓉, 等. 水产动物行为及其在渔业中的应用研究进展[J]. 水生生物学报, 2021, 45(5): 1171-1180. [百度学术]
LI W, RONG K, QIN L R, et al. Research progress of aquatic animal behavior and its application in fisheries[J]. Acta Hydrobiologica Sinica, 2021, 45(5): 1171-1180. [百度学术]
丑晴. 湍流环境下垂直轴双风机气动性能优化研究[D]. 哈尔滨: 哈尔滨工业大学, 2021. [百度学术]
CHOU Q. Optimization analysis of aerodynamic performance of twin vertical-axis wind turbines in turbulent environments[D]. Harbin: Harbin Institute of Technology, 2021. [百度学术]
姜宜辰, 赵月, 熊济时, 等. 水下航行器艇体形状对阻力及流噪声综合影响[J]. 哈尔滨工程大学学报, 2022, 43(1): 76-82, 138. [百度学术]
JIANG Y C, ZHAO Y, XIONG J S, et al. Comprehensive influence of underwater vehicle hull shape on resistance and flow noise[J]. Journal of Harbin Engineering University, 2022, 43(1): 76-82, 138. [百度学术]
郭斌, 杨雷, 倪庆, 等. 新一代载人飞船气动外形设计与优化[J]. 宇航学报, 2023, 44(9): 1337-1349. [百度学术]
GUO B, YANG L, NI Q, et al. Aerodynamic shape design and optimization of China’s next-generation manned spacecraft[J]. Journal of Astronautics, 2023, 44(9): 1337-1349. [百度学术]
何隆, 张亚, 李世中, 等. 基于响应面方法的AUV水动力外形优化[J]. 兵器装备工程学报, 2022, 43(12): 43-50. [百度学术]
HE L, ZHANG Y, LI S Z, et al. Hydrodynamic profile optimization of AUV based on response surface method[J]. Journal of Ordnance Equipment Engineering, 2022, 43(12): 43-50. [百度学术]
LUO W L, GUO X M, DAI J W, et al. Hull optimization of an underwater vehicle based on dynamic surrogate model[J]. Ocean Engineering, 2021, 230: 109050. [百度学术]