仿生蝠鲼胸鳍摆动推进机构设计与水动力分析
Design and hydrodynamic analysis of pectoral fin oscillation propulsion mechanism of bionic manta ray
收稿日期: 2018-03-30
Received: 2018-03-30
作者简介 About authors
云忠(1971—),男,博导,从事水下无人推进系统研究.orcid.org/0000-0003-4455-5967.E-mail:
为了提高水下无人航行器推进系统效率,设计仿生蝠鲼胸鳍摆动推进机构及其摆动形式. 在研究蝠鲼胸鳍运动机理的基础上,设计摆动幅值与急回系数相互解耦的曲轴联合曲柄摆杆机构,以实现仿生胸鳍的特定摆动形式. 推导并建立仿生胸鳍的鳍面运动方程,利用数值仿真分析单侧胸鳍在该摆动形式下的水动力特性,根据设定的推进性能评价指标说明该摆动形式相对正弦摆动形式的优越性. 通过实验验证装置的前进、转向等运动性能,分析胸鳍摆动幅值、频率与运动速度的关系. 结果表明,在特定的摆动形式下,该仿生胸鳍摆动推进机构能够实现水下无人航行器的前进、转向等预定运动,满足基本推进要求.
关键词:
A bionic pectoral fin oscillation propulsion mechanism of manta ray and its swing mode were designed to improve the propulsive efficiency of unmanned undersea vehicle (UUV). A crankshaft joint crank rocker mechanism, which achieved mutual decoupling between swing amplitude and jerk factor, was designed to achieve the specific swing mode of the bionic pectoral fin, based on the study of the pectoral fin movement mechanism of manta rays. The motion equation of the bionic pectoral fin was derived and established, and the hydrodynamic performance of unilateral pectoral fin in the swing mode was analyzed by numerical simulation, and the superiority of this swing mode relative to the sine swing mode was verified according to the set propulsion performance evaluation index. The advancement and steering performance of the device was verified by the experiments, and the relationship between the amplitude and frequency of the pectoral fin swing as well as the speed of the device was analyzed. Experimental results showed that the bionic pectoral fin oscillation propulsion mechanism can realize the predetermined movement of UUV, such as advancement and steering, and meet the basic propulsion requirements.
Keywords:
本文引用格式
云忠, 温猛, 蒋毅, 陈龙, 冯龙飞.
YUN Zhong, WEN Meng, JIANG Yi, CHEN Long, FENG Long-fei.
目前国内外学者对仿生蝠鲼UUV展开了相关研究,并取得了一定的成果. 李吉等[6]根据蝠鲼胸鳍的波动规律研制了具有双柔性胸鳍的仿生鱼BH-RAY3,该装置利用刚性胸鳍前缘的拍动和柔性胸鳍的被动变形产生运动,最大航速为450 mm/s. Wang等[7-9]设计了形状记忆合金丝驱动的仿生蝠鲼机器鱼,利用记忆合金丝的变形拉动柔性鳍面运动,最大直线游动速度为79 mm/s,最小转弯半径为118 mm. Chen等[10-11]利用离子聚合物-金属复合材料作为人造胸鳍的致动器以模拟蝠鲼的拍动过程,设计了2种带有人造胸鳍的UUV,航速分别为0.055 BL/s(每秒0.055倍体长)和0.067 BL/s. Zhou等[12-14]结合蝠鲼胸鳍拍动和滑翔原理设计并改进仿生UUV RoMan-Ⅱ,最大航速为400 mm/s. Gao等[15]在分析蝠鲼胸鳍运动特征和骨骼结构的基础上,设计了带有柔性硅胶胸鳍的机器鱼,最大航速为1.4 BL/s. Chew等[16]分析蝠鲼胸鳍运动机理并设计了仿蝠鲼摆动推进的机器鱼,最大航速为1.783 BL/s. Wang等[17]基于蝠鲼的结构外形,设计了由浮力驱动的仿生蝠鲼水下滑翔机,并利用仿真软件分析装置浮力等随重心位置的变化情况. 现有的研究成果大多仅涉及仿生蝠鲼胸鳍的运动机理、机构设计等,较少涉及到装置在直线行驶和转弯等状态下的稳定性分析,而UUV在各状态下的稳定性恰恰是限制航行器性能的重要因素. 此外,现有仿生装置柔性胸鳍的运动大多仅由刚性胸鳍前缘拍动来带动,无法实现对胸鳍弦向方向柔性部位局部性状的精确控制,装置运动的灵活性受限.
本研究设计由多根鳍条驱动的仿生胸鳍摆动推进机构,实现机构摆动幅度与急回系数的解耦,建立该机构的运动学方程,并通过软件仿真和实验验证分析该胸鳍在仿生装置各运动状态下的水动力学性能.
1. 仿生蝠鲼水下航器摆动推进机构设计
1.1. 仿生原型的生物学特征分析
蝠鲼是生活在热带和亚热带海底的软骨鱼类,身体扁平,呈菱形,拥有肥厚如翼状的胸鳍和细长的尾鳍,可以通过胸鳍的摆动实现远程洄游、凌空飞跃等运动,行动敏捷. 如图1所示,蝠鲼胸鳍的摆动包括外展、内收和前缘叠拢等一系列运动过程,可以分解为弦向和展向2个方向的波动运动,其中弦向波动是胸鳍推进力的主要来源,而展向波动主要用于提高系统的运动稳定性. 蝠鲼胸鳍的摆动属于典型的升力运动模式,相对划动运动等阻力运动模式在高速运动时更为高效. 研究发现,在蝠鲼胸鳍扇动过程中弦向波数约为0.4,且在下拍过程中存在急回运动,因而其弦向波形并非单纯的正弦波.
图 1
以胸鳍前沿延长线与中线交点为坐标原点o,弦向运动方向为x轴正方向,展向运动方向为y轴正方向建立坐标系. 为了简化模型,假设胸鳍为柔性薄板,且鳍面运动波形为正弦波,则蝠鲼胸鳍摆动的表达式[18]为
式中:Ac(x)、As(y)分别为弦向和展向波幅,
1.2. 仿生蝠鲼胸鳍推进机构设计
图 2
如图3所示,推进机构模仿蝠鲼以鳍条带动鳍面运动的形式,采用多根刚性鳍条以一定初相位差带动柔性硅胶板鳍面运动,从而实现鳍面特定的波动运动. 为了实现鳍条特定规律的摆动,对比分析常见的凸轮式、曲柄式和正弦式等能将转动运动转换为摆动运动的机构方案. 凸轮式机构可实现多种摆动规律,但摆动幅度较小;曲柄式相对凸轮式耗材较少,且不存在较大幅度下的干涉现象;正弦式通常由曲柄滑块和摆杆等机构组合而成,机构尺寸较大,且难以满足急回特性需求. 综上,曲柄式转化机构相对较好. 常用的曲柄式转换机构有曲柄导杆机构、曲柄滑块机构和曲柄摆杆机构等. 曲柄导杆机构的摆杆初端尺寸较大,对流体特性有较大影响;当曲柄滑块机构的摆杆做对称摆动时,其摆动幅度与急回系数相互耦合,无法针对特定幅度设计特定急回系数;曲柄摆杆机构能通过改变连杆长度设计出相同摆动幅度、不同急回系数的转换机构,实现机构摆动幅度与急回系数的解耦. 综合考虑推进能力与机构尺寸等因素,认为单侧单电机驱动曲轴联合曲柄摆杆机构是相对较优的选择.
图 3
图 3 仿生蝠鲼胸鳍摆动推进机构简图
Fig.3 Schematic diagram of pectoral fin oscillation propulsion mechanism of bionic manta ray
2. 水动力学分析
2.1. 机构运动学建模
建立胸鳍机构的运动学模型,确定电机输出参数与胸鳍运动参数之间的数学关系. 以右侧胸鳍为例,基本运动关系如图4所示. 图中,l1为曲柄长度;l2为连杆长度;l3为鳍条长度;l4为基座间距;
图 4
由三角函数关系可得
由图4中关系可得
式中:P、Q、K、H为与结构尺寸相关的常量,用于简化公式结构.
当曲柄AB匀速转动时,
式中:
假设
鳍条的摆动方程可以写为
式中:
鳍面运动由多根不同初始相位角的鳍条连续变化组成. 假设鳍面运动的波长为
设
2.2. 水动力数值仿真
采用较为流行的计算流体动力学(computational fluid dynamics,CFD)分析软件Fluent作为仿真软件,并利用动网格技术,通过用户自定义函数(user defined functions, UDF)编程控制胸鳍的摆动运动,以获得所需的翼展面流场信息,从而研究翼展面形状对胸鳍推进性能的影响.
2.2.1. 网格划分及参数设定
因胸鳍厚度相对翼展面积对推力的影响可忽略不计,假设胸鳍为无厚度曲面,胸鳍鳍基线为静止基准线. 根据式(4)、(9),取胸鳍摆动频率为2 Hz,波长为0.75 m,则鳍面运动周期T的运动序列如图5所示. 为了使仿真流场充分发展,应在来流方向取相对较长的长度,其余方向取较小的长度以减少计算量,本研究计算流域尺寸为6 000 mm×4 000 mm×2 000 mm;考虑到计算精度和计算时间的平衡,将计算流域分为内域和外域,内域采用非结构化四面体网格以适应复杂的模型特征,外域采用结构化网格以节省计算时间,网格划分结果如图6所示,网格规模约为480万. 进口边界采用速度边界,出口边界采用压力出口,壁面采用无滑移壁面. 采用有限体积法求解三维流域的雷诺平均N-S方程,湍流模型采用对复杂曲面运动有较好适应性的Realizablek-
图 5
图 6
图 6 胸鳍摆动仿真所采用的计算域网格配置
Fig.6 Calculate domain mesh configuration for pectoral fin swing simulation
2.2.2. 仿真结果及分析
通过水动力学仿真胸鳍在静水中的运动,得到在周期T内,在如图5所示运动序列下的鳍面压力云图序列,如图7所示. 图中,p为胸鳍拍动时上表面压力. 由图7可知,在胸鳍拍动过程中,鳍面上交替出现高压中心和低压中心. 由图7(a)~(c)可见,高压中心逐渐向鳍尖移动并趋于消失,同时低压中心开始出现;由图7(d)~(f)可见,低压中心产生并逐渐扩散至胸鳍中部;由图7(f)~(g)可见,低压中心向鳍尖收缩并逐渐消失;由图7(h)~(a)可见,高压中心开始产生并逐渐扩散至胸鳍中部. 综上所述,在胸鳍拍动过程中,高压中心和低压中心均从胸鳍前端开始出现并逐步移向鳍尖,鳍条能够较好地带动鳍面整体运动;此外,胸鳍尖部压力变化较大,根部压力变化较小,因而鳍根处的摆动对推进性能的影响较小.
图 7
图 7 胸鳍拍动时上表面不同时刻压力云图
Fig.7 Pressure contours of upper surface with pectoral fin flapping at different times
3. 推进机构性能分析
3.1. 性能评价指标
当仿生蝠鲼胸鳍在水下运动时,右侧鳍面受力示意图如图8所示. 设鳍面上每一单元节点相对于流体的速度为
图 8
图 8 水下运动时右侧胸鳍受力示意图
Fig.8 Force diagram of right pectoral fin during underwater movement
式中:Mkx、Mky、Mkz分别为相对坐标原点的滚转力矩、俯仰力矩和偏航力矩.
设胸鳍拍动时来流速度为U,则单侧胸鳍的推进效率为
式中:Fxm为胸鳍水平推力平均值,θkx为单元k在x方向上的摆角. 则胸鳍推进效率与滚转力矩的绝对值直接相关,且在推进力均值和推进速度均值一定时,与绝对扭转力矩均值
式中:
3.2. 胸鳍摆动水动力曲线
综合考虑推进力与效率,根据推力系数均值、升力系数均值分别与急回系数的关系,仿真得到最佳急回系数为0.56. 取仿生胸鳍的摆动频率为2 Hz,波长为0.75 m,分别对所设计的胸鳍采用式(9)的波动形式和正弦波动形式进行仿真计算,得到胸鳍在静水中各水动力参数变化曲线,如图9所示.
图 9
由图9可知,当采用式(9)的摆动形式时,
在相同条件下取不同的摆动频率,进行胸鳍在静水中的水动力学分析,仿真得到相应的水动力学参数随摆动频率f变化的关系曲线如图10所示. 由图可知,胸鳍的
图 10
图 10 胸鳍水动力参数随摆动频率的变化曲线
Fig.10 Variation curve of hydrodynamic parameters of pectoral fin versus swing frequency
4. 实验研究
设计实验样机,以验证仿生结构和胸鳍摆动方案的可行性. 主要验证仿生蝠鲼推进装置的推进能力、转向能力以及胸鳍摆动参数对装置游动能力的影响. 实验现场为人工湖,如图11所示. 利用高清摄像机对装置的运动姿态进行拍摄,然后通过图像处理得到装置的位姿数据,进而计算得到装置的运动速度等参数.
图 11
图 11 仿生蝠鲼推进装置实验现场情况
Fig.11 Experimental scene of propulsion device of bionic manta ray
实验样机是对设计方案的比例缩放,总长为600 mm,翼展为1 060 mm,曲轴摆杆机构传动比约为1∶4,输出摆幅A约为24°,则曲轴所需输入力矩和功率表示式为
式中:
图 12
图 12 理论输入转矩/功率与曲轴转速的关系
Fig.12 Relationship between theoretical input torque as well as power and crankshaft rotate speed
为了验证装置的推进性能,分别进行直线运动、转向运动、俯仰运动实验,结果如下. 1)直线运动. 如图13所示为当仿生装置胸鳍鳍条摆幅为15°、频率为15 Hz时的直线运动序列图. 在直线运动试验中,发现仿生装置存在左右偏转现象. 导致装置方向偏转的因素主要有胸鳍摆动波长、频率、幅度和初始相位等,其中波长、幅度、频率均可通过预先控制而保证较小误差,但各鳍条之间的初始相位较难保证. 分析对比仿真结果可知,单侧胸鳍在仿真实验中基本全程产生正向推力,因此想要保证直线运动的精度须使用精度较高的电机以控制鳍条的初始相位. 2)转向运动. 如图14所示为当鳍条最大摆幅为24°、频率为2 Hz时的转弯运动序列图. 结果表明装置转速约为35 °/s,不能实现原地转向,只能够实现较小半径的转向运动,分析认为是胸鳍向前和向后的推进力大小不同所致. 随着胸鳍摆动频率的增大,装置转向速度逐渐增大,而仿真结果表明胸鳍摆动频率与推进力成正相关,两者结果吻合较好. 因此,可以通过调整胸鳍的摆动频率来调整装置的机动能力. 3)俯仰运动. 在俯仰运动试验中,观察装置在直线运动时的升降情况和前端俯仰情况,发现其在运动过程中产生较小程度的上升、下降运动和前端俯仰运动. 对比仿真结果发现在胸鳍摆动过程中存在一定的升力分量和俯仰力矩,两侧胸鳍升力分量的叠加产生了装置的升降运动,而俯仰力矩的叠加产生了俯仰运动. 升降运动可以通过添加重心调节装置加以控制,而俯仰运动须通过改进胸鳍摆动方式来尽量消除,从而提高装置整体运动的稳定性. 实验进一步验证了摆动幅度和频率对胸鳍推进力的影响. 在实验中分别设置鳍条摆幅为6.0°、10.5°、15.0°、19.5°、24.0°这5个对照组,取一定梯度的鳍条摆动频率,通过测量装置的运动速度来间接测量胸鳍的推进力,装置的运动速度V与胸鳍摆动幅度、频率之间的关系如图15所示. 由图可以看出,胸鳍的推进力与摆动频率成正相关,且在一定范围内胸鳍推进力与摆幅成正相关.
图 13
图 14
图 15
图 15 装置速度随胸鳍摆动频率变化曲线
Fig.15 Variation curve of device speed versus swing frequency of pectoral fin
5. 结 论
(1)在仿生装置运动过程中,摆动机构能够通过多根鳍条驱动柔性鳍面整体按照预定的规律运动,但须采用较精密的控制方式对装置的运动过程进行控制,以提高装置的运动稳定性.
(2)在仿生装置前进运动过程中,仿生胸鳍摆动方式几乎能够全程产生正向推力,且相对正弦摆动方式有较大的推进力和较高的推进效率.
(3)在一定范围内,仿生胸鳍推进力的大小随摆动频率和摆动幅度的增加而逐渐增大.
(4)所设计的仿生蝠鲼胸鳍摆动形式的推进效率相对正弦摆动形式有一定的提高,但仍然和自然鱼类有一定差距. 下一步的工作将针对胸鳍的更多摆动方式和运动过程中胸鳍柔性部分具体状态对装置运动性能的影响进行进一步分析,以改善装置整体的水动力性能.
参考文献
仿生机器鱼研究进展及发展趋势
[J].
Research development and tendency of biomimetic robot fish
[J].
Pectoral fin locomotion in batoid fishes: undulation versus oscillation
[J].
On the relationship between efficiency and wake structure of a batoid-inspired oscillating fin
[J].
Performance limits of labriform propulsion and correlates with fin shape and motion
[J].
仿生蝠鲼机器鱼BH-RAY3的研制及水力实验
[J].
Development and hydrodynamics experiments of robotic manta ray BH-RAY3
[J].
形状记忆合金驱动仿生蝠鲼机器鱼的设计
[J].
Development of a biomimetic manta ray robot fish actuated by shape memory alloy
[J].
仿生机器魟鱼研制与游动性能实验研究
[J].
Development and swimming experimental research on bionic stingray robot
[J].
Bio-inspired robotic manta ray powered by ionic polymermetal composite artificial muscles
[J].DOI:10.1080/19475411.2012.686458 [本文引用: 1]
Better endurance and load capacity: an improved design of manta ray robot: RoMan-Ⅱ
[J].
Design and locomotion control of a biomimetic underwater vehicle with fin propulsion
[J].DOI:10.1109/TMECH.2011.2175004
Analysis of motion of batoid fins for thrust generation by using fluid-structure interaction method
[J].DOI:10.3795/KSME-A.2010.34.11.1575 [本文引用: 2]
/
〈 |
|
〉 |
