第28卷第3期 水下无人系统学报 Vol.
28
No.3
2020年6月
JOURNAL OF UNMANNED UNDERSEA SYSTEMS Jun. 2020
收稿日期: 2019-10-11; 修回日期: 2019-11-19.
基金项目: 国家重点研发计划重点专项(2017YFC0305902); 青岛海洋科学与技术国家实验室“问海计划”项目(2017WHZZB
0101); 天津市自然科学基金重点基金(18JCZDJC40100).
作者简介: 孙秀军(1981-), 男, 教授, 硕导, 主要研究方向为海洋移动观测平台技术.
252 Journal of Unmanned Undersea Systems www.yljszz
[引用格式] 孙秀军, 王力伟, 桑宏强. 波浪滑翔器水下牵引机滑翔动力分析[J]. 水下无人系统学报, 2020, 28(1): 252-258. 波浪滑翔器水下牵引机滑翔动力分析
实验室分析天平
孙秀军1,2,3, 王力伟2, 桑宏强4
(1. 河北工业大学 机械工程学院, 天津, 300130; 2. 中国海洋大学 物理海洋教育部重点实验室, 山东 青岛, 266100; 3. 青岛海洋与技术试点国家实验室, 山东 青岛, 266237; 4. 天津工业大学 机械工程学院, 天津, 300387) 摘 要: 波浪滑翔器的纵向速度与水下牵引机结构参数有很大关系, 其中水下牵引机自重和水翼最大 俯仰角对其影响尤为突出。文中以“海哨兵”波浪滑翔器为研究对象, 采用牛顿-欧拉方程建立水翼动力学模型, 对最优自重及最优角度进行计算流体力学仿真, 并将仿真结果与试验数据进行对比分析。分析结果可知, 受到多种因素的影响, 两者数据存在一定偏差, 但两者总体趋势相近, 数据显示20°俯仰角为水翼的最优俯仰角, 3级海况条件下, 45 kg 的自重使得水下牵引机在下滑时间和下滑位移方面更加匹配波浪特性。文中在部分参数确定的情况下, 研究自重和俯仰角对波浪滑翔器推进效率的影响, 可为波浪滑翔器后续的参数优化和设计提供参考。 关键词: 水下牵引机; 水翼; 波浪滑翔器; 自重; 俯仰角
中图分类号: U674.941; O353.2 文献标识码: A 文章编号: 2096-3920(2020)03-0252-07 DOI: 10.11993/j.issn.2096-3920.2020.03.002
Gliding Dynamics Analysis for Underwater Tractor of Wave Glider
美容喷雾器
SUN Xiu-jun 1,2,3, WANG Li-wei 1, SANG Hong-qiang 4
(1. School of Mechanical Engineering, Hebei University of Technology, Tianjin 300130, China; 2. Physical Oceanography Laboratory, Ocean University of China, Qingdao 266100, China; 3. Qingdao National Laboratory for Marine Science and Technology, Qingdao 266237, China; 4. School of Mechanical Engineering, Tianjin Polytechnic University, Tianjin 300387, China)
Abstract: The longitudinal velocity of wave glider is closely related to the structural parameters of the underwater trac-tor, especially, the self-weight of the underwater tractor and the maximum pitch angle of the hydrofoil are crucial. In this study, a dynamic model of the hydrofoil of the Sea Sentry wave glider was established by using the Newton-Euler equa-tions to simulate the optimum weight of the underwater tractor and the optimum pitch angle of the hydrofoil in compu-tational fluid dynamics(CFD) simulation software. Furthermore, the simulation results were compared with the test data. Due to the influence of various factors, there is a certain deviation between the two kinds o
f data, but the general trends are similar, and the data show that 20° pitch angle is the optimal pitch angle of hydrofoil. Under the level three sea con-dition, the self-weight of 45 kg can make the underwater tractor more suitable for wave characteristics in terms of glid-ing time and gliding displacement. The effects of the weight and pitch angle on propulsion efficiency were studied under the condition that some parameters were fixed, which provided reference for the subsequent parameter optimization and design of wave glider.
Keywords: underwater tractor; hydrofoil; wave glider; self-weight; pitch angle
2020年6月
孙秀军, 等: 波浪滑翔器水下牵引机滑翔动力分析 第3期
水下无人系统学报 www.yljszz 253
0 引言
波浪滑翔器是一种能够将波浪能转化为前向动能的水面无人航行器, 具有工作时间长、位置自主可控等优点, 可广泛应用于水文气象观测、水质生态调查、卫星遥感验证以及水面通信中继等领域[1-3]。当前, 如何提升波浪滑翔器的运动速度逐渐成为该领域的研究焦点之一。提升波浪滑翔器的运动速度也就是提高波浪滑翔器的波浪动力转化效率。波浪动力转化效率主要取决于波浪滑翔器的双体结构参数以及水下牵引机的结构参数等。对于水下牵引机来说, 其自重、翼板翼型和翼板弹簧机构等参数选择, 以及不同参数的组合对波浪转化效率都会产生不同程度的影响, 在各个参数相互耦合的关系中, 出波浪动力转化效率的优化路线以及完成耦合参数的优选, 对波浪滑翔器的运动速度提升至关重要。
刘鹏等
[4]
基于雷诺平均(Reynolds average
Navier-Stockes, RANS)方程利用二维面元法计算了摆幅角、翼间距及升沉幅度等因素对水下牵引机推进性能的影响, 并给出小翼间距可产生更大推力的结论。Yang 等[5]基于URANS 方程利用计算流体力学(computational fluid dynamics, CFD)软件FINE/Marine 与STAR-CCM+联合仿真了波浪滑翔器整个运动周期。Caiti 等[6]利用拉格朗日方法对波浪滑翔器纵垂面进行建模分析, 求解了总体位移和速度与时间的对应关系。贾立娟[7]利用Kane 方法对波浪滑翔器进行动力学建模, 分析了波浪滑翔器水下牵引机升阻比的影响因素。但上述研究主要集中在理论层面, 尚未在工程设计中给出水下牵引机细节机构的参数选择依据, 例如, 水下牵引机自重、弹簧劲度系数、水翼限位轴与转动轴的距离等与推进效率的对应关系。
文中以“海哨兵”波浪滑翔器的牵引机结构为研究对象, 在相关结构参数确定的情况下分析水下牵引机自重与波浪滑翔器滑翔效率的对应关系, 到已有结构下水下牵引机的最优自重, 进而确定最优俯仰角。最后, 对波浪滑翔器的实际滑翔效率进行了试验, 在趋势上验证了理论分
析、流体仿真和真实情况之间的一致性。
1 结构组成及工作原理
波浪滑翔器是由水面浮体船、柔性铠装缆、水下牵引机组成的双体结构系统。水面浮体船在波浪的作用下产生上下起伏运动, 通过柔性铠装缆拉拽水下牵引机上下运动, 受迫垂向运动的水下牵引机水翼与水质点发生相对滑翔运动, 水质点对水翼产生前向作用力, 水下牵引机产生向前运动, 反之, 水下牵引机通过柔性铠装缆带动水面浮体船前向运动, 其运动过程如图1所示。
图1 波浪滑翔器运动原理 Fig. 1 Motion principle of wave glider
为便于求解水下牵引机滑翔效率, 做出如下假设: 1) 波浪滑翔器的垂向运动完全符合波浪周期T , 且为周期性正弦运动; 2) 以6对水翼代表整个水下牵引机, 不考虑水下牵引机其余部件的水动力影响; 3) 水下牵引机向下滑落时水翼保持最大滑翔角的时间为整个波浪周期的1/5, 下落位移为波高H 的1/2, 此状态下水下牵引机在水中处于自由落体状态, 柔性铠装缆对水下牵引机不施加外力作用。
根据波浪滑翔器运动特性, 其水翼的运动规律如图2所示。
图2 水翼运动曲线
Fig. 2 Motion curves of hydrofoil
2 数值模型
2.1 动力学模型
水翼受力情况如图3和图4所示, 利用牛顿-
2020年6月
水下无人系统学报 第28卷
254 Journal of Unmanned Undersea Systems www.yljszz
折叠麻将桌欧拉方程[8]建立水翼下滑过程中纵向速度u 、垂向速度w 的动力学模型, 用以揭示水翼运动与自重和俯仰角的特征关系, 水翼动力学方程的显式表达形式为
3233cos cos (1)()d sin sin (1)()d cos sin cos sin k k k y x z
zz F F u k u k t m F F G
w k w k t m GL F L F L F L L δβ
δβθθθθΩ⎧-+=+
⋅⎪⎪
⎪+-+=+
⋅⎨⎪
---⎪'=⎪⎩
(1)
式中: m 为水翼质量; d t 为仿真时间步长;β为弹
簧与水平方向夹角; F 为水动力合力; δ为水动力合力与水平方向夹角; G 为水翼重力; Ω'为水翼绕轴旋转的角加速度; L zz 为水翼相对于质心z 方向的转动惯量; F x 、F y 分别为F 在x 、y 方向上的分力; F k 为弹簧拉力, 由于文中着重研究水翼处于最大俯仰角的状态, 此状态下, 弹簧与水翼的相对状态不变, 弹簧拉力为恒定值。
图3和图4中: α为水翼攻角; θ为水翼俯仰角; L 1为弹簧拉拽点与质心的水平距离; L 2为旋转轴与限位轴的距离; L 3为旋转轴与质心的距离。
图3 水翼下滑受力情况示意图
Fig. 3
Forces on hydrofoil in sliding state
图4 水翼上滑受力情况示意图
Fig. 4 Forces on hydrofoil in up-sliding state
文中由于水翼下滑时间t=t 1, 垂向位移y ≤y 1，因此, 水下牵引机水翼的最优方程可以表达为式
(2), 式中相关参数均由“海哨兵”波浪滑翔器已有数据测得, 其几何参数及动力学参数如表1所示。
1
110101
(cos cos )d (1)d (sin sin )d (1)d t k x t k y F δF t x u k t m F δF G t y w k t m y y ββ--⎧-⎡
⎤=-+⎪⎢⎥⎣⎦⎪⎪⎨+-⎡⎤=-+⎪⎢⎥⎣⎦
⎪⎪⎩∑∑
≤ (2)
表1 波浪滑翔器几何参数及动力学参数
Table 1 Geometric parameters and dynamic parame-ters of wave glider
参数描述 数值 m /kg m 1 G /N 9.8 m 1 F k /N 37.8 α/(º)
N/A
β/(º) 1.51 θ/(º) 20 δ/(º) N/A L 1/mm 123.5 L 2/mm 39.5 L 3/mm 24.04
2.2 仿真模型及参数
文中采用类NACA0012翼型[9]单对水翼进行仿真, 水翼特征弦长b =0.16 m, 质心位置由SolidWorks 确定并使之与原点重合, 利用CFD 计
算软件FLUENT 对二维水翼进行仿真, 仿真模型计算域为15c ×20c , c 为水翼弦长, 左边界和上边界距离翼型质心的距离均为5c , 全局网格采用非结构性网格, 边界层为10层, 并严格控制第1层网格高度Y +值小于60, 如图5所示。
图5 模型边界及水翼近壁网格
Fig. 5 The boundary of model and the grid nearby
hydrofoil
2020年6月
酸洗缓蚀剂孙秀军, 等: 波浪滑翔器水下牵引机滑翔动力分析 第3期
水下无人系统学报 www.yljszz 255
考虑到流体中的非稳状态, 水翼的动量方程和能量方程均较难显式表达, 因此文中在FLUENT 仿真中采用标准-k ε湍流模型, 利用非平衡壁面函数对水翼壁面进行修正[4]。 2.3 网格无关性验证
为验证网格的可行性, 分别采用Grid1、Grid2、Grid3和Grid4这4种尺寸网格下比较水翼在20︒攻角下的升力系数差异, 如表2所示, 其中设置流体密度ρ=998.2 kg/m, 截面平均流速v =1 m/s, 动力粘性u =1.003×10–3 Pa·s, 翼型特征长度d =0.06 m, 根据公式Re=ρvd/u 得出雷诺数为5.9×105。根据表2可知, Grid 3与Grid 4尺寸下水翼的升力系数差值变化较小, 网格尺寸对计算精度的影响减弱, 考虑到计算的准确性及计算资源的合理运用, 因此采用Grid 3网格尺寸进行下一步仿真计算。
表2 FLUENT 网格无关性验证
Table 2 Grid independence verification in Fluent
网格 数量 攻角/(º)流速/(m·s –1)
升力系数
Grid 1 80 406 20 1 0.309 51 Grid 2 179 270 20 1 0.417 84 Grid 3 716 458 20 1 0.535 45 Grid 4
1 99
2 009
20
1
0.558 16
3 数值计算仿真与分析
太平洋海况波谱如表3所示[10], 据表显示太平洋中3级海况最为常见, 由于海况的复杂性和不确定性增加了对海浪系统性仿真的难度, 因此拟采用太平洋中最常见的3级海况进行仿真以达到简化的效果。
表3 太平洋海况波谱
Table 3 State spectrum of Pacific Ocean
浪级
波高/m
周期/s
二维码打印设备浪级概率/%
0~1 <0.25 <2
0.70
2 0.50
3 6.80 3 1.00 5 27.53
4 2.00 7 24.00
5 3.00 10 20.61
6 5.00 12 13.15
≥7 >7.50 14~24 7.21
利用6自由度动网格方法对水翼20°俯仰角下5~9 kg 不同水翼质量进行仿真[11], 得到水翼纵向速度曲线见图6。从图中可知, 当水翼质量为2 kg 时, 在初始的一段时间内由于自身质量较轻, 水翼尾缘分流效果作用明显, 导致其存在非稳状态: 水翼纵向速度达到约0.7 m/s 后, 其运动状态逐渐稳定, 且在0.25 s 下滑时间内, 该质量水翼的纵向速度远超其余情况纵向速度; 当下滑时间较大时, 该质量水翼运动状态趋于稳定, 纵向加速度较小, 纵向速度小于其余质量水翼的纵向速度。5~9 kg 水翼的纵向运动较为相似, 纵向速度随着质量的增加而增加。
图6 不同质量水翼纵向速度随时间变化曲线
Fig. 6 Curves of longitudinal velocity of hydrofoil ver-sus time under different mass
图7给出了水翼的纵向位移曲线, 从图中可以看出, 7.5 kg 水翼更加吻合假定3)中关于下滑时间和下滑位移的设定, 因此将45 kg(7.5 kg×6)质量视为水下牵引机最优自重进行下一步仿真。
图7 不同质量水翼纵向位移随时间变化曲线
Fig. 7 Curves of longitudinal displacement of hydrofoil
versus time under different mass
在最优自重情况下, 分别对10°、15°、20°
和25°俯仰角的纵向位移和下滑时间进行仿真, 其结果如图8和图9所示。从图中可以看出, 当水翼俯仰角在10°~20°范围内时, 其纵向位移随
2020年6月
水下无人系统学报 第28卷
256 Journal of Unmanned Undersea Systems www.yljszz
着俯仰角的增大而增大; 在20°时达到峰值; 当水翼俯仰角为25°时, 其纵向位移在初始一段时间内较大, 但由于其俯仰角较大导致升力系数减小, 该俯仰角度的水翼在后续时间的纵向速度小于其他水翼。
图8 不同俯仰角水翼纵向位移随时间变化曲线
Fig. 8 Curves of longitudinal displacement of hydrofoil
versus time under different pitch angles
图9 不同俯仰角水翼纵向速度随时间变化曲线 Fig. 9 Curves of longitudinal velocity of hydrofoil ver-sus time under different pitch angles
根据仿真数据绘得1 s 时间内不同俯仰角的水翼纵向位移曲线, 如图10所示。可以看出, 相同下滑时间时, 20°俯仰角有较大纵向位移; 当俯仰角为25°时, 水翼的纵向位移逐渐减小, 因此认为20°为水下牵引机在最优自重情况下的最优俯仰角。
图10 1 s 滑翔时间内水翼纵向位移随俯仰角变化仿真曲线 Fig. 10 Simulation curve of longitudinal displacement of
hydrofoil versus pitch angle in one second of gliding time
4 水池试验与数据对比分析
4.1 水池试验
为验证仿真数据的准确性, 现利用一定机构进行水池试验验证, 其结构如图11所示。定滑轮将水下牵引机固定在一定高度, 利用单一变量法将水下牵引机的结构参数根据仿真条件进行相关调整, 具体试验结构参数如表4所示。将水下牵引机在水中自由释放, 观察记录水下牵引机垂向位移0.5 m 所需时间及其下滑的纵向位移。图12为水下牵引机水池试验现场照片。
图11 水下牵引机水池试验装置结构图
Fig. 11 Structure of pool test device of underwater tr-
oahactor
表4 水下牵引机试验结构参数
Table 4 Test structure parameters of underwater trac-tor
自重/kg
俯仰角/(°)
垂向位移/m
30~54 20 0.5 45 10~25 0.5
图12 水下牵引机水池试验
Fig. 12 Pool test of underwater tractor
4.2 数据对比分析
由于试验装置的局限性, 并未对 2 kg 和
20 kg 翼重对应的水下牵引机自重进行试验, 利用去最值平均法对每组参数测得的数据进行处理, 将处理后的数据与仿真数据进行对比, 其结果如图13~图16所示。
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议