水陆两栖球形机器人辅助鳍的水动力特性

第52卷第5期2021年5月
中南大学学报(自然科学版)
Journal of Central South University (Science and Technology)
V ol.52No.5
May 2021
水陆两栖球形机器人辅助鳍的水动力特性
迟兴，战强
(北京航空航天大学机械工程及自动化学院，北京，100191)
摘要：为将传统球形机器人的应用领域从陆地拓展到水域，通过在球壳上黏附辅助鳍，增大球形机器人水面运动的推进力，设计一种紧凑高效的水陆两栖球形机器人。

分析特殊平面形状的单片辅助鳍在运动周期中入水面积和形心的变化规律，基于机翼理论推导液体作用在辅助鳍上的升力和阻力计算公式，给出水平分力和竖直分力表达式。

利用傅里叶级数对水平分力和竖直分力进行拟合，针对球壳上不同数量辅助鳍配置方式，给出傅里叶级数表示的水平分力合力和竖直分力合力计算公式。

研究结果表明，增大球壳旋转角速度和辅助鳍入水面积可提高水平分力；水平分力合力呈现周期性变化规律；增加配置的辅助鳍数量可增大水平分力合力的平均值，减小合力的波动幅度；辅助鳍结构能够给球形机器人水面运动提供持续稳定的推进力。

关键词：球形机器人；水陆两栖；辅助鳍；升力；阻力中图分类号：TP242.6
文献标志码：A
文章编号：1672-7207（2021）05-1456-09
Hydrodynamic performance of assistant fins of an amphibious
spherical robot
CHI Xing,ZHAN Qiang
(School of Mechanical Engineering and Automation,Beihang University,Beijing 100191,China)
Abstract:To expand the application of a traditional spherical robot from the terrestrial environment to the aquatic environment,a compact and efficacious amphibious spherical robot design was proposed by attaching assistant fins onto the outer spherical shell to increase the robot's propulsion force.The wetted surface area of one assistant fin and its centroid were calculated,then equations of the lift force and the drag force were derived based on the airfoil theory,and afterwards equations of the horizontal force and the vertical force were presented.The horizontal force and the vertical force were fitted separately using Fourier series,and after that,the equations of the resultant force exerted on all those fins submerged into the water were given.The results show that increasing the rotational speed of the spherical shell and the wetted surface area of the fin increase the horizontal force.
DOI:10.11817/j.issn.1672-7207.2021.05.007
收稿日期：2020−05−30；修回日期：2020−07−06
基金项目(Foundation item)：国家自然科学基金资助项目(50705003)(Project(50705003)supported by the National Natural Science
Foundation of China)
通信作者：战强，博士，教授，从事机器人机构设计与控制、球形机器人和仿生灵巧手研究；E-mail:**************.cn
引用格式：迟兴,战强.水陆两栖球形机器人辅助鳍的水动力特性[J].中南大学学报(自然科学版),2021,52(5):1456−1464.Citation:CHI Xing,ZHAN Qiang.Hydrodynamic performance of assistant fins of an amphibious spherical robot[J].Journal of Central South University(Science and Technology),2021,52(5):1456−
1464.
第5期迟兴，等：水陆两栖球形机器人辅助鳍的水动力特性
Attaching more fins onto the spherical shell increases the average of the resultant horizontal force,which changes periodically,and decreases its undulation.The assistant fins provide steady propulsion force for the spherical robot.
Key words:spherical robot;amphibious;assistant fin;lift force;drag force
球形机器人将驱动机构、控制电路和感知系统等重要部件封闭在球形外壳的内部，对其形成密闭保护。

该结构特点使得球形机器人容易适应复杂严苛的外部环境，在环境探测、安保防爆、军事侦察等领域具有广阔应用前景[1]。

传统的球形机器人应用局限于陆地环境，通过配备额外的驱动能力可以拓展球形机器人的应用环境到水域，使其具有水陆两栖特性，增强其环境适应性。

为此，需要解决传统球形机器人在水域中运动时驱动能力不足这一关键问题。

一些研究者通过给传统球形机器人配备适合于水中运动的驱动装置，研制了不同类型的水陆两栖球形机器人。

GUO等[2−4]介绍了一种由复合驱动系统提供动力的水陆两栖球形机器人，该复合驱动系统由4个相同的驱动单元组成，每个驱动单元使用1个喷水推进器和2个垂直布置的伺服电机，有四足步行模式和喷水推进模式2种运动模式，分别用于机器人的陆地运动和水中运动。

KAZNOV等[5]研发了水陆两栖球形机器人GroundBot，该机器人利用重心偏移原理实现陆地运动和水面运动，球形外壳上设计了点状凸起带，以提高机器人的水面移动速度。

李艳生等[6−7]设计了一种利用螺旋桨推进的水陆两栖球形机器人，机器人中心位置处安装了1个内部包含推进螺旋桨的柱形涵道，该机器人在陆地上利用重心偏移实现滚动运动，在水中时则利用螺旋桨推进的方式运动。

通过在球壳两极地带黏附辅助鳍，提高水面运动推进力，可实现利用1套驱动装置完成球形机器人两栖运动的目标。

已有一些研究参考机翼理论对辅助鳍进行了力学特性分析[8]。

HEALEY等[9]给出了AUV叶片整周旋转中的升力和阻力计算公式，GEORGIADES等[10]通过实验探究了叶片受力和往复运动振幅和周期的关系，BEHBAHANI 等[11−12]研究了仿生鱼胸鳍划水动作时的水动力特性，READ等[13−14]给出了叶片俯仰和升沉运动的推进力，刘鹏等[15−16]分析了运动参数对水翼推进性能的影响。

YANG等[17−18]利用CFD方法给出了水轮机叶片之间的速度分布和压力分布，张万超等[19−20]分析了叶轮受浪流扰动产生横摇运动时的水动力性能。

在常规的叶片水动力分析中，叶片通常整体浸入液体中，其湿面积保持不变，压力中心也不变。

本文提出利用辅助鳍改善球形机器人水面运动性能的设计方案，通过在球壳上黏附辅助鳍，设计一种紧凑高效的水陆两栖球形机器人。

分析特殊平面形状的单片辅助鳍在运动周期中入水面积和形心的变化规律，基于机翼理论推导液体作用在辅助鳍上的升力和阻力计算公式，给出水平分力和竖直分力表达式。

利用傅里叶级数对水平分力和竖直分力进行拟合，针对球壳上不同数量辅助鳍配置方式，给出傅里叶级数表示的水平分力合力和竖直分力合力计算公式。

1水陆两栖球形机器人系统简介
附加辅助鳍的水陆两栖球形机器人由内部驱动单元和球壳2部分组成，内部驱动单元包括机器人的动力系统和控制电路，球壳由左、右两半球壳组成，作为运动执行部件，同时对内部驱动单元形成密闭保护，如图1所示。

黏附在球壳外表面两极地带的辅助鳍数量可灵活配置，一种典型的配置方式是在左右两半球壳上各安装8
片辅助鳍，且位于半球壳上的辅助鳍
图1水陆两栖球形机器人
Fig.1Amphibious spherical robot
1457
第52卷中南大学学报(自然科学版)关于中心轴呈放射性分布，其示意图如图2所示。

图中O −XYZ 为固定在球壳上的随体坐标系，Y 轴与中心轴重合，R 为球壳半径；1~8代表8个辅助鳍，相邻辅助鳍之间的相位差为45°。

辅助鳍ABC 区域的横截面形状近似为钝角三角形，在球壳表面的安装位置由ζ1和ζ2确定，ζ1，ζ2和顶角ζt 共同决定了辅助鳍的平面形状。

当机器人执行陆地运动时，位于底部的辅助鳍高于地面，即使机器人产生沿X 轴方向的滚转扰动，辅助鳍也不易与地面发生碰撞；当机器人执行水面运动时，位于底部的辅助鳍部分浸没到液体中，液体会在辅助鳍上产生额外的推进力。

2单片辅助鳍力学特性
2.1
单片辅助鳍升力和阻力
辅助鳍受力分析如图3所示，在辅助鳍上建立惯性坐标系P −XYZ ，其中原点P 为辅助鳍的压力中心，坐标轴Z 垂直纸面向外，倾角φ为辅助鳍平面从X 轴正向逆时针转过的角度。

液体作用在辅助鳍上的力可以分解为升力L 和阻力D ，升力方向垂直于来流方向，阻力方向与来流方向一致。

考虑机器人在静止液体表面的运动，当球壳
绕中心轴旋转时，机器人整体沿X 轴负向以速度V s 移动。

辅助鳍随球壳一起平移和旋转，其相对液体的速度V f 是牵连速度V s 和相对速度V n 的矢量和，即
V f =V s +V n
(1)
相对速度V n 方向垂直于辅助鳍平面，幅值为
|V n
|=ωs l
PY。

ωs 为球壳相对中心轴的旋转角速度，
l PY 为压力中心P 到中心轴间的距离。

运用运动转换原理，即认为机器人是固定的，
观察液体相对机器人的运动，可以认为液体以一定速度冲击到辅助鳍上，从而对辅助鳍施加力的作用。

相对辅助鳍的来流速度为V w ，V w 与V f 幅值相等而方向相反，即
V w =-V s -V n
(2)
记速度矢量V s 和-V n 之间的夹角为β，β∈[0,π/2]，则β与倾角φ之间的关系如下：
β=
{
π/2-φ,φ-π/2,φ≤π/2
φ>π/2(3)
来流速度V w 的幅值为
|V w
|=|
|V s
2
+||V n 2
-2||V s ||V n cos β
(4)
V w 和V n 之间的夹角γ的取值范围为[0,π/2]，γ
可通过下式计算：
γ=arcsin (||V s sin β/|
|V w )
(5)
攻角αe 为来流速度V w 与辅助鳍平面之间的夹角，攻角αe 和γ之间的关系如下：
αe =
{
π/2-γ,π/2+γ,φ≤π/2
φ>π/2
(6)
(a)侧视图；(b)正视图图2辅助鳍安装示意图
Fig.2
Installation diagram of assistant
fins
图3辅助鳍受力分析Fig.3
Force analysis of a fin
1458
第5期
迟兴，等：水陆两栖球形机器人辅助鳍的水动力特性当液体以速度V w 冲击到辅助鳍上时，液体对辅助鳍施加的合力作用于压力中心，合力可以分解为升力L 和阻力D [10]：
L =1
2
ρwater |V w |2SC L max sin (2αe )(7)
D =1
2
ρwater |V w |2SC D max (1-cos (2αe ))
(8)
式中：ρwater 为水的密度，S 为浸入水中的辅助鳍平
面的面积；C L max 为最大升力系数，C D max 为最大阻力系数。

由式(7)和(8)可知，升力和阻力与来流速度幅值|V w |和面积S 正相关。

在辅助鳍随球壳运动的过程中，|V w |和S 均为辅助鳍倾角φ的函数，这使得升力和阻力的计算变得复杂。

将作用于辅助鳍上的升力和阻力在惯性坐标系P-XYZ 中分解为
F X =-L sin (φ+αe )+D cos (φ+αe )(9)F Y =L cos (φ+αe )+D sin (φ+αe )
(10)
式中：F X 和F Y 分别为液体对单片辅助鳍的作用力在水平方向和竖直方向的分量。

水平分力F X 的方向与机器人运动方向相同，表现为机器人的推进力。

竖直分力F Y 对机器人运动影响较小，仅进行简单讨论。

2.2
工作周期中的关键角度
辅助鳍ABC 工作周期中临界状态如图4所示。

辅助鳍前缘A 和中心轴之间的距离为l AY ，l AY =R sin ζ2；浸入液面下的球壳高度为h 。

当前缘A 刚
好和液面接触，即将开始浸入液体的临界状态，倾角φe 为
φe =arcsin (()R -h /l AY )
(11)
在辅助鳍逆时针转过一定角度，前缘A 即将离
开液面的临界状态，倾角φd 为
φd =π-φe
(12)
当辅助鳍顶点C 刚好位于液面上时，倾角φc 为
φc =arcsin
R -h
l AY -R ()cos ζ1-cos ζ2tan ()
ζt -π/2(13)
在机器人水面运动过程中，当某辅助鳍倾角φe ≤φ≤φd 时，该辅助鳍部分浸没到液体中，液体
能够对其施加力的作用；否则，该辅助鳍位于液面上方，液体对其无力的作用。

2.3
辅助鳍入水面积和形心
辅助鳍的升力和阻力与其浸入液体中部分的面积和形心有关。

以前缘A 点为原点，建立局部坐标系A −uv ，其中uv 平面与辅助鳍平面平行，u 轴与中心轴平行，如图5所示。

原点A 在随体坐标系O −XYZ 中的坐标为
(R sin ζ
2
cos φ,R cos ζ2,R sin ζ2sin φ)。

在局部坐标系
A −uv 中，直线AC 和BC 的方程分别为
v =tan (ζt -π/2)u
u =R (cos ζ1-cos ζ2),
圆弧
AB 的方程为
(u +R cos ζ2)
2
+(v -R sin ζ2)2
=R 2。

当φe ≤φ≤φd 时，浸入液体中的辅助鳍长
度l im 为
l im =l AY -()R -h /sin φ
(14)
由于辅助鳍平面形状近似为三角形，针对顶点C 高于液面或低于液面的情况，计算浸入液体中的面积和形心。

当顶点C 高于液面时，即φe ≤φ
≤
图4辅助鳍工作周期中的临界状态Fig.4Critical states in a work cycle of
fin
(a)φe ≤φ≤φc ；(b)φc <φ≤π/2图5
辅助鳍入水面积示意图
Fig.5
Wetted surface area of a fin
1459
第52卷
中南大学学报(自然科学版)φc ，浸入液体中的面积S im1为
S im1=
∫0l im
∫
R 2
-()
v -R sin ζ2
2
-R cos ζ2
v /tan ()
ζt
-π2d u d v
(15)
形心在局部坐标系A −uv 中的坐标(u P 1,v P 1)可通过下式计算：
u P 1=M V 1/S im1(16)v P 1=M U 1/S im1
(17)
式中：M V 1和M U 1为浸入液体中辅助鳍平面形状的一阶矩：
M V 1=∫0l im ∫
R 2-(
)
v -R sin ζ22
-R cos ζ2
v /tan ()ζt -π/2
u d u d v
(18)M U 1=
∫0
l im ∫
R 2-(
)
v -R sin ζ2
2
-R cos ζ2
v /tan ()ζt -π/2
v d u d v
(19)
当顶点C 低于液面时，即φc <φ≤π/2，记L C =R (cos ζ1-cos ζ2)tan (ζt -π/2)，浸入液体中的
辅助鳍面积为
S im2=
∫0L C ∫
R 2-(
)v -R sin ζ22
-R cos ζ2
v /tan ()ζt -π/2
d u d v +
∫L C
l im ∫
R 2-(
)
v -R sin ζ2
2
-R cos ζ2
R ()
cos ζ1-cos ζ2
d u d v
(20)
形心的坐标(u P 2,v P 2)为：
u P 2=M V 2/S im2(21)v P 2=M U 2/S im2
(22)
式中：M V 2和M U 2为浸入液体中辅助鳍平面形状的一阶矩：
M V 2=
∫0L C ∫
R 2-(
)
v -R sin ζ2
2
-R cos ζ2
v /tan ()ζt -π/2
u d u d v +
∫L C
l im ∫
R 2-(
)
v -R sin ζ22
-R cos ζ2
R ()cos ζ1-cos ζ2
u d u d v
(23)
M U 2=
∫0
L C ∫
R 2
-(
)
v -R sin ζ2
2
-R cos ζ2
v /tan ()ζt -π/2
v d u d v +
∫L C
l im ∫
R 2-(
)
v -R sin ζ2
2
-R cos ζ2
R ()cos ζ1-cos ζ2
v d u d v
(24)
忽略不同深度压强不同引起的压力分布微小差异，假设液体作用在辅助鳍上的力在浸入液体中的部分均匀分布，则辅助鳍的压力中心P 与浸入液体中的部分的形心重合。

基于形心的计算公式，可以得到压力中心P 到中心轴之间的距离l PY 的表达式：
l PY =
{
l AY -v P 1,l AY -v P 2,φe ≤φ≤φc φc <φ≤π/2
(25)
2.4
数值分析结果与讨论
水陆两栖球形机器人的参数为R =175mm ，ζ1=
13.7°，ζ2=55.2°，ζt =115°，h =116.4mm ，代入式(11)~(13)计算可得关键角度φe =24.1°，φc =31.9°；φd =155.9°。

取C Lmax =0.92；C Dmax =1.12，进行数值仿真分析。

浸入水中的辅助鳍面积S 和l PY 随倾角φ的变化曲线如图6所示。

由图6可见：面积变化曲线和l PY 变化曲线关于φ=π/2对称。

在辅助鳍刚浸入液面和即将离开液体的阶段，S 和l PY 随倾角φ变化较为显著；而在辅助鳍接近与液面垂直的阶段，面积和l PY 随倾角φ变化比较平缓。

球形机器人运动稳定后，即球壳以额定转速ωs =2πrad/s 旋转，机器人以额定平移速度|V s |=
0.4m/s 做直线运动时，来流速度幅值|V w |和攻角αe 随倾角φ变化的曲线如图7所示。

由图7可知：倾角φ从φe 增大到π/2的过程中，来流速度幅值|V w |快速减小，攻角先减小再逐渐增大；|V w |的曲线关于φ=π/2对称，攻角的曲线关于点(π/2,π/2)中心对称，该点同时是攻角曲线的拐点；|V w |持续减小的主要影响因素为V s 和-V n 之间的夹角β减小。

运动稳定后，作用在单片辅助鳍上的升力和阻力变化规律如图8所示。

从图8可见：在倾角φ从φe 增大到π/2的过程中，由于辅助鳍入水面积增大，阻力先从0增大到最大值，然后，由于来流速度幅值|V w |下降成为主导，阻力缓慢下降；在φ=π/2即辅助鳍平面与液面垂直时，阻力取得局部极
小值，局部极小值约为阻力最大值的一半，此时，
图6
面积S 和l PY 变化曲线
Fig.6
Changing curves of area S and length l PY
1460
第5期
迟兴，等：水陆两栖球形机器人辅助鳍的水动力特性
|V w
|也取到极小值；阻力变化曲线关于φ=π/2对
称。

当φe ≤φ≤π/2时，升力的变化趋势与阻力一
致，但是在φ>π/2后，升力的符号发生了变化。

原因在于，以从X 轴正向逆时针转过的角度为升力和阻力的相位，当φ<π/2时，升力的相位超前阻力的相位，当φ>π/2后，升力的相位滞后阻力的相位，φ=π/2时升力为0，是相位变化的临界点。

升力符号的改变反映了升力和阻力间相位的变化。

水平分力和竖直分力随倾角φ变化的规律如图9所示。

从图9可见：水平分力和阻力有相似的变化规律，但两者的符号相反，竖直分力和升力有相似的变化规律。

根据水平分力计算公式，运动周期中，当倾角φ接近φe 或φd 时，在水平分力的成分中，升力的水平分量占比较大，阻力的水平分量增大了该阶段水平分力的幅值。

在倾角φ增大到π/2的过程中，升力的水平分量占比减小，阻力的
水平分量占比逐渐增大；当倾角φ接近π/2时，阻力的水平分量在水平分力中起主导作用，该阶段升力趋向于0，升力对水平分力的贡献很小。

在运动周期的中间阶段，水平分力整体变化较为平缓，辅助鳍能给机器人提供持续稳定的推进力。

参考升力和阻力计算公式可知，增大辅助鳍入水面积和球壳旋转速度可提高中间阶段水平分力的幅值。

3液体作用在辅助鳍上的合力
3.1水平分力和竖直分力计算公式
水平分力和竖直分力关于倾角φ的函数表达式
比较复杂，为利于进一步计算和分析，有必要对其形式进行简化。

由于水平分力和竖直分力呈现明显周期性变化规律，可用傅里叶方法进行拟合，得到水平分力和竖直分力的傅里叶级数近似表达。

m 阶傅里叶级数的方程为
F (m ,φ)=a 0+∑i =1
m a i cos ()iw φ+b i sin (iw φ)(26)
将上述方程中的基写为列向量形式：
B F (m ,φ)=[1,cos (w φ),sin (w φ),cos (2w φ),]
sin (2w φ),⋯,cos ()mw φ,sin ()
mw φT
(27)
相应的系数写为行向量形式A F (m )=[a 0
a 1
b 1
a 2
⋯a m b m ](28)
则m 阶傅里叶级数可以表示为系数向量与基向量的积：
F (m ,φ)=A F (m )B F (m ,φ)
(29)
图7来流速度幅值和攻角变化曲线
Fig.7
Changing curves of amplitude of inlet water and
angle of
attack
图9水平分力和竖直分力变化曲线
Fig.9
Changing curves of horizontal force and
vertical
force
图8升力和阻力变化曲线
Fig .8Changing curves of lift force and drag force
1461
第52卷
中南大学学报(自然科学版)液体作用在辅助鳍上的水平分力可以用傅里
叶级数表示的分段函数来描述：
F X (φ)=
{
F ()m ,φ,φe ≤φ≤φd
0,0≤φ<φe ,φd <φ≤2π
(30)
定义带通滤波函数
D (m ,φ)=
{
I 2m +1,φe ≤φ≤φd
0,0≤φ<φe ,φd <φ≤2π
(31)
则上述水平分力又可以表示为
F X (φ)=A F X (m )D (m ,φ)B F X (m ,φ)
(32)
类似地，竖直分力可以表示为
F Y (φ)=A F Y (m )D (m ,φ)B F Y (m ,φ)
(33)
其中：A F X (m )和A F Y (m )分别为水平分力和竖直分力对应的系数向量，基向量B F X (m ,φ)和B F Y (m ,φ)的区别在于各自的频率w 不同。

3.2
傅里叶级数表示的合力计算公式
当半球壳上配置的辅助鳍数量为n ∈N 时，任意相邻的辅助鳍之间的相位差为2π/n 。

以辅助鳍1作为参考基准，其倾角为φ时，辅助鳍k 的相位为
φk =
{
φ-2()k -1π/n ,k ≤()n +2/2φ+2()n -k +1π/n ,k >()n +2/2
(34)
当φe ≤φ≤φd 且φe ≤φk ≤φd ，辅助鳍1和辅
助鳍k 同时位于液面下方，液体作用在辅助鳍k 上的力计入合力。

辅助鳍k 相应的带通滤波函数为D (m ,k ,φ)=ìíîï
ï
I 2m +1,φe +2()k -1π/n ≤φ≤φd ,φe ≤φ≤φd -2()n +k -1π/n 0,其他
(35)
液体作用于辅助鳍k 的水平分力为
F kX (φ)=F X (φk )=A F X (m )D (m ,k ,φ)B F X (m ,φk )
(36)
当F kX (φ)≠0时，将F kX (φ)的表达式用三角函数和差公式展开，分离出其中的常值系数项，表示为系数向量和基向量的积：
F kX (φ)=A F X (m )D (m ,k ,φ)W kX (m ,k )B F X (m ,φ)
(37)
式中：W kX (m ,k )为系数矩阵，
W kX ()m ,k =
éëêêêêêêêêêêùû
úú
úú
úúúúú
ú1000⋯000cos ϕk sin ϕk 0⋯00
0-sin ϕk cos ϕk 0⋯00
00
cos (2ϕk )⋯
00
⋮⋮⋮⋮⋮⋮⋮0000⋯cos (m ϕk )sin (m ϕk )0
⋯
-sin (m ϕk )cos (m ϕk )(38)
ϕk =
{
2()k -1w π/n ,k ≤()n +2/2
-2()n -k +1w π/n ,k >()n +2/2
(39)
对作用在半球壳上浸入液面下的所有辅助鳍上的水平分力求和：
F t X (φ)=∑k =1
n F kX ()φ=
∑k =1
n A
F X
()m D ()m ,k ,φW kX ()m ,k B F X ()m ,φ(40)
记
W t X (m ,φ)=∑k =1
n
D (m ,k ,φ)W kX ()
m ,k (41)
则水平分力的合力可以简化为
F t X (φ)=ΑF X (m )W t X (m ,φ)B F X (m ,φ)
(42)类似地，竖直分力的合力可简化为
F t Y (φ)=A F Y (m )W t Y (m ,φ)B F Y (m ,φ)
(43)
3.3数值分析结果与讨论
在给定的机器人参数下，利用傅里叶级数对
单片辅助鳍水平分力和竖直分力进行拟合，结果如图10所示。

从图10可见：在φe ≤φ≤φd 范围内，当傅里叶级数的阶数m =6时，水平分力拟合
相应的和方差SSE 为0.0002432，相关系数R 2为0.9998；竖直分力拟合相应的SSE 为0.004267，R 2为0.9995
，显示拟合曲线与原始曲线具有较高的图10水平分力和竖直分力拟合曲线Fig.10
Fitting curves of horizontal force and
vertical force
1462
第5期迟兴，等：水陆两栖球形机器人辅助鳍的水动力特性
匹配度。

将水平分力拟合得到的频率和系数向量代入水平分力合力表达式中，计算辅助鳍数量n 分别为4，8，12时的水平分力合力，结果如图11所示。

从图11可以看出：对不同的辅助鳍数量配置，水平分力合力均呈现周期性变化规律。

增加配置的辅助鳍数量可以增大合力的平均幅值，平抑合力的波动，从而减小振幅，同时，合力的振荡频率也会增加。

其中，当辅助鳍数量为12时，水平分力合力平均幅值约为0.72N ，振幅约为平均幅值的3.5%，是较为合适的配置方式。

两半球壳上的辅助鳍相互对称配置，故作用在球壳整体的合力等于作用于半球壳的合力的2倍。

4结论
1)利用辅助鳍改善球形机器人水面运动性能，将传统球形机器人的应用环境从陆地拓展到水面，增强了其环境适应性。

2)基于机翼理论推导了液体作用在辅助鳍上的升力和阻力计算公式，给出了水平分力和竖直分力表达式，并利用傅里叶级数对水平分力和竖直分力进行拟合，给出了傅里叶级数表示的水平分力合力和竖直分力合力计算公式。

3)辅助鳍在入水和出水临界状态之间的中间阶段，液体作用在辅助鳍上的水平分力变化平缓，辅助鳍能给机器人提供持续稳定的推进力；提高
球壳的旋转角速度和增大辅助鳍的入水面积都可以增大中间阶段作用在辅助鳍上的水平分力，从而提升机器人的推进力；作用在半球壳上辅助鳍上的水平分力合力呈现周期性变化规律，增加辅助鳍数量，可以增大合力的平均值，减小振幅。

参考文献：
[1]
战强,李伟.球形移动机器人的研究进展与发展趋势[J].机械工程学报,2019,55(9):1−17.
ZHAN Qiang,LI Wei.Research progress and development trend of spherical mobile robots[J].Journal of Mechanical Engineering,2019,55(9):1−17.[2]
GUO Shuxiang,LI Maoxun,YUE Chunfeng.Performance evaluation on land of an amphibious spherical mother robot in different terrains[C]//IEEE International Conference on Mechatronics and Automation.Takamatsu,Japan:IEEE,
2013:1173−1178.
[3]
LI Maoxun,GUO Shuxiang,HIRATA H,et al.Design and performance evaluation of an amphibious spherical robot[J].Robotics and Autonomous Systems,2015,64:21−34.[4]
LI Maoxun,GUO Shuxiang,HIRATA H,et al.A roller-skating/walking mode-based amphibious robot[J].Robotics and Computer-Integrated Manufacturing,2017,44:17−29.[5]
KAZNOV V ,SEEMAN M.Outdoor navigation with a spherical amphibious robot[C]//IEEE/RSJ International Conference on Intelligent Robots and Systems.Taipei,Taiwan,China:IEEE,2010:5113−5118.[6]
李艳生,孙汉旭,贾庆轩,等.一种摆式球形水下机器人水底滚动特性分析[J].中南大学学报(自然科学版),2016,47(11):3664−3669.
LI Yansheng,SUN Hanxu,JIA Qingxuan,et al.Analysis of rolling characteristics at water bottom for spherical robot with heavy pendulums[J].Journal of Central South University(Science and Technology),2016,47(11):3664−3669.[7]
LI Yansheng,YANG Meimei,SUN Hanxu,et al.A novel amphibious
spherical
robot
equipped
with
flywheel,
pendulum,and propeller[J].Journal of Intelligent &Robotic Systems,2018,89(3/4):485−501.[8]
盛振邦.船舶原理[M].2版.上海:上海交通大学出版社,2017:23−25.
SHENG Zhenbang.Naval principal[M].2nd ed.Shanghai:Shanghai Jiao Tong University Press,2017:23−25.[9]
HEALEY A J,ROCK S M,CODY S,et al.Toward an improved understanding of thruster dynamics for underwater vehicles[J].IEEE Journal of Oceanic Engineering,1995,20(4):
354-361.
图11不同辅助鳍数量对应的水平分力合力Fig.11
Resultant horizontal force of different
numbers of fin
1463
第52卷中南大学学报(自然科学版)
[10]GEORGIADES C,NAHON M,BUEHLER M.Simulation
of an underwater hexapod robot[J].Ocean Engineering, 2009,36(1):39−47.
[11]BEHBAHANI S B,TAN Xiaobo.Design and modeling of
flexible passive rowing joint for robotic fish pectoral fins[J].
IEEE Transactions on Robotics,2016,32(5):1119−1132. [12]MORGANSEN K A,TRIPLETT B I,KLEIN D J.
Geometric methods for modeling and control of free-swimming fin-actuated underwater vehicles[J].IEEE Transactions on Robotics,2007,23(6):1184−1199.
[13]READ D A,HOVER F S,TRIANTAFYLLOU M S.Forces
on oscillating foils for propulsion and maneuvering[J].
Journal of Fluids and Structures,2003,17(1):163−183. [14]BAIK Y S,BERNAL L P.Experimental study of pitching
and plunging airfoils at low Reynolds numbers[J].
Experiments in Fluids,2012,53(6):1979−1992.
[15]刘鹏,苏玉民,刘焕兴,等.串列异步拍动翼推进性能分析
[J].上海交通大学学报,2014,48(4):457−463.
LIU Peng,SU Yumin,LIU Huanxing,et al.Propulsive performance analysis of tandem asynchronous flapping foil [J].Journal of Shanghai Jiao Tong University,2014,48(4): 457−463.
[16]刘鹏,苏玉民,李宏伟.运动参数对近波面拍动翼推进性能
的影响[J].中南大学学报(自然科学版),2016,47(12): 4062−4069.
LIU Peng,SU Yumin,LI Hongwei.Effect of kinematic parameters on propulsion performance of flapping foil[J].
Journal of Central South University(Science and Technology),2016,47(12):4062−4069.
[17]WANG Tao,WANG Chuan,KONG Fanyu,et al.
Theoretical,experimental,and numerical study of special impeller used in turbine mode of centrifugal pump as turbine [J].Energy,2017,130:473−485.
[18]YANG Sunsheng,KONG Fanyu,CHEN Hao,et al.Effects
of blade wrap angle influencing a pump as turbine[J].
Journal of Fluids Engineering,2012,134(6):061102. [19]张万超,王凯,王树齐,等.受迫横摇运动浮式立轴叶轮的
水动力性能分析[J].华中科技大学学报(自然科学版), 2018,46(8):18−23.
ZHANG Wanchao,WANG Kai,WANG Shuqi,et al.The hydrodynamic performance analysis of floating vertical-axis tidal turbine considering forced rolling motion[J].Journal of Huazhong University of Science and Technology(Natural Science Edition),2018,46(8):18−23.
[20]王凯,张万超,王树齐,等.不同横摇幅值下浮式立轴叶轮
的水动力性能分析[J].哈尔滨工程大学学报,2019,40(9): 1569−1575.
WANG Kai,ZHANG Wanchao,WANG Shuqi,et al.
Hydrodynamic performance analysis of floating vertical-axis tidal turbine with different rolling motion amplitudes[J].
Journal of Harbin Engineering University,2019,40(9): 1569−1575.
(编辑赵俊)
1464。