水陆两栖机器人水下动力学建模及操纵性研究

陈龙; 童保成; 祖斯羽; 张豪; 胡桥; CHEN Long; TONG Baocheng; ZU Siyu; ZHANG Hao; HU Qiao

引 en

水陆两栖机器人水下动力学建模及操纵性研究

陈龙¹

机构：

1. 西安交通大学机械工程学院，陕西西安 710049

×
，童保成¹

机构：

1. 西安交通大学机械工程学院，陕西西安 710049

×
，祖斯羽¹

机构：

1. 西安交通大学机械工程学院，陕西西安 710049

×
，张豪¹

机构：

1. 西安交通大学机械工程学院，陕西西安 710049

×
，胡桥^1,2

机构：

1. 西安交通大学机械工程学院，陕西西安 710049

2. 西安交通大学陕西省智能机器人重点实验室，陕西西安 710049

×

1. 西安交通大学机械工程学院，陕西西安 710049；
2. 西安交通大学陕西省智能机器人重点实验室，陕西西安 710049；

Research on Underwater Dynamics Modeling and Maneuverability of Amphibious Robots

CHEN Long¹

Affiliation：

1. School of Mechanical Engineering，Xi′an Jiaotong University，Xi′an 710049 ，China

×
， TONG Baocheng¹

Affiliation：

1. School of Mechanical Engineering，Xi′an Jiaotong University，Xi′an 710049 ，China

×
， ZU Siyu¹

Affiliation：

1. School of Mechanical Engineering，Xi′an Jiaotong University，Xi′an 710049 ，China

×
， ZHANG Hao¹

Affiliation：

1. School of Mechanical Engineering，Xi′an Jiaotong University，Xi′an 710049 ，China

×
， HU Qiao^1,2

Affiliation：

1. School of Mechanical Engineering，Xi′an Jiaotong University，Xi′an 710049 ，China

2. Shanxi Key Laboratory of Intelligent Robots，Xi′an Jiaotong University，Xi′an 710049 ，China

×

1. School of Mechanical Engineering，Xi′an Jiaotong University，Xi′an 710049 ，China；
2. Shanxi Key Laboratory of Intelligent Robots，Xi′an Jiaotong University，Xi′an 710049 ，China；

作者简介:

陈龙（1996-），男，硕士生，主要从事水下仿生机器人结构设计与水动力仿真、水陆两栖机器人等研究。

通讯作者:

胡桥（1977-），男，博士，教授，主要从事海洋智能感知、海洋仿生机器人等研究。

中图分类号:TP242.6

文献标识码:A

文章编号:2096-5753(2023)04-0487-09

DOI:10.19838/j.issn.2096-5753.2023.04.012

全文
评论
参考文献
出版信息

参考文献 1

WU Z Y，QI J，ZHANG S.Amphibious robots：a review[J].Applied Mechanics and Materials，2014，494-495：1036-1041.

查找原文

参考文献 2

RAFEEQ M，TOHA S F，AHMAD S，et al.Locomotion strategies for amphibious robots-a review[J].IEEE Access，2021，9：26323-26342.

查找原文

参考文献 3

FRANCESCHI A，PIAGGIO B，TONELLI R，et al.Assessment of the manoeuvrability characteristics of a twin shaft naval vessel using an open-source CFD code[J].Journal of Marine Science and Engineering，2021，9（6）：665.

查找原文

参考文献 4

GO G，AHN H T.Hydrodynamic derivative determination based on CFD and motion simulation for a tow-fish[J].Applied Ocean Research，2019，82：191-209.

查找原文

参考文献 5

DUBBIOSO G，BROGLIA R，ZAGHI S.CFD analysis of turning abilities of a submarine model[J].Ocean engineering，2017，129：459-479.

查找原文

参考文献 6

SHEN Z R，WAN D C，CARRICA P M.Dynamic overset grids in OpenFOAM with application to KCS self-propulsion and maneuvering[J].Ocean Engineering，2015，108：287-306.

查找原文

参考文献 7

尹盛林.仿生波动鳍水陆两栖机器人设计与控制技术研究[D].西安：西安交通大学，2021.

查找原文

参考文献 8

王雷.水下机器人运动控制研究[D].合肥：中国科学技术大学，2018.

查找原文

参考文献 9

胡芳琳，冯康佳，马向能.动态PMM实验的数值模拟及线性导数获取方法[J].水动力学研究与进展：A 辑，2021，36（2）：9.

查找原文

参考文献 10

施生达.潜艇操纵性[M].北京：国防工业出版社，1995.

查找原文

参考文献 11

FOSSEN T I.Guidance and control of ocean vehicles[M].New York：John Willey & Sons，1994.

查找原文

目录contents

摘要Abstract
关键词Keywords
0 引言
1 机器人坐标系建立及数值求解方法
1.1 坐标系建立
1.2 计算模型
1.3 网格划分
1.4 仿真参数设定
2 水动力系数求解与动力学建模
2.1 纯纵荡运动
2.2 水动力系数求解汇总
2.3 五自由度动力学及运动学建模
3 机器人水下操纵性运动特性研究
3.1 运动仿真数值计算方法
3.2 机器人运动仿真结果及分析
3.2.1 机器人直航运动仿真
3.2.2 机器人水平面回转运动
4 机器人水下操纵性运动实验
4.1 直航运动实验
4.2 水平面回转运动
5 结束语
参考文献

摘要

针对多模态水陆两栖机器人作业环境复杂使得水下运动状态难以预报等难题，基于 CFD 方法求解的水动力系数，构建了机器人水下运动的五自由度动力学和运动学模型。基于机器人水下动力学模型，采用四阶经典龙格库塔法，开展了机器人直航运动及水平面回转运动数值仿真研究，并进行了水池试验验证。试验结果数据与数值仿真结果误差均不超过 10%，验证了机器人水动力系数及五自由度动力学模型的准确性，为水陆两栖机器人研制提供了理论与技术支撑。

Abstract

In view of the difficulties in predicting the underwater motion state of a multimodal amphibious robot caused by complex working environment，a five-degree-of-freedom dynamic and kinematics model for underwater motion of a robot was constructed based on the hydrodynamic coefficient solved by CFD method. Based on the underwater dynamics model of the robot，the fourth order classical Runge-Kutta methods is used to conduct numerical simulation research on the direct sailing motion and horizontal rotation motion of the robot，and pool experiments are carried out for verification. The error between the experimental data and the numerical simulation results is less than 10%，which verifies the accuracy of hydrodynamic coefficient and five-degree-of-freedom dynamic model of the robot，providing reliable theoretical guidance for robot development.

关键词

水陆两栖机器人； CFD ；动力学模型；数值仿真

Keywords

amphibious robots ； CFD ； dynamic model ； numerical simulation

0 引言
性差和稳定性弱等问题，本文将仿生推进机构与传统推进机构结合设计出一种高性能组合推进型水陆两栖机器人。机器人的水下操纵性影响其工作状态，因此研究机器人的水动力系数以及水下操纵性有重要意义。性差和稳定性弱等问题，本文将仿生推进机构与传统推进机构结合设计出一种高性能组合推进型水陆两栖机器人。机器人的水下操纵性影响其工作状态，因此研究机器人的水动力系数以及水下操纵性有重要意义。
通过 CFD 方法建立机器人水下动力学模型的方法已通过标模试验数据的验证，得到了国内外学者的广泛应用。2021 年，FRANCESCHI 等人利用开源 CFD 代码 OpenFoam 求解了船体本地和附加质量的水动力系数，完成了对双轴舰船的操纵性质量评估^[3]，研究表明，应用 CFD 流体分析方法可以得到较为准确的结果。2019 年，GO 等人提出了一种利用计算流体动力学技术求解拖鱼水下航行器水动力导数的新方法^[4]，并进行了不同速度下拖鱼在 3 种运动场景的六自由度仿真，最后通过虚拟仿真结果与理论对比分析，验证了方法的有效性^[3]。2017 年，意大利 DUBBIOSO 等人利用非定常 CFD 求解器详细分析了潜艇的机动特性，针对十字形和 X 形舵进行水平面的三自由度运动分析，并通过实验进行验证，研究结果表明 X 形舵具备更加优异的转弯能力^[5]。上海交通大学万教授团队基于 OpenFOAM 仿真平台，开发了一款船舶自航行求解器 naoe-FOAM-SJTU^[6]，实现了船舶航行的直接数值模拟。结合上述研究现况分析可知，基于 CFD 的间接求解方法计算效率和精度高，且易于参数化，便于进行动力学模型的修正和机器人优化设计。
本文以一种利用波动鳍推进的轮–鳍–桨机器人为研究对象，通过 CFD 技术针对机器人水下多种运动模态进行数值求解，分析得到机器人不同运动模态下的水动力特性，求解水动力参数，建立机器人的动力学模型。构建水下运动仿真平台，通过机器人直航、回转等运动模态研究其水下操纵性，并验证动力学模型的准确性和有效性^[6]，为进一步研究多模态水陆两栖机器人操纵特性奠定了理论与技术基础。
1 机器人坐标系建立及数值求解方法
1.1 坐标系建立
如下图1所示，本文采用了国际水池会议（ITTC）推荐的体系，建立水陆两栖机器人的坐标系：
大地坐标系 $O_{e} x_{e} y_{e} z_{e}$ –固定坐标系。该坐标系 O_e 为空间中任一点，且 $O_{e} z_{e}$ 正向指向地心， $O_{e} x_{e}$ 和 $O_{e} y_{e}$ 位于水平面内且互相垂直。
随体坐标系 $O_{b} x_{b} y_{b} z_{b}$ –运动坐标系。原点与机器人浮心重合， $O_{b} x_{b}$ 轴沿机器人前进方向为正， $O_{b} y_{b}$ 轴沿机器人右侧方向为正， $O_{b} z_{b}$ 轴垂直于 $O_{b} x_{b} y_{b}$ 平面指向下为正。
图1 水陆两栖机器人坐标系示意
Fig.1 Schematic diagram of the coordinate system for amphibious robots
鳍面坐标系 $O_{f} x_{f} y_{f} z_{f}$ –鳍面坐标系。各坐标轴正向与波动鳍动力学建模分析中坐标系相同。
1.2 计算模型
本文以机器人三维模型为基础构建仿真计算模型，计算模型将其内部填充为实体，保留整体的外形轮廓；同时针对局部细节进行简化处理，提高计算速率。最终构建的机器人水下仿真模型如下图2 所示，其尺寸为 1 150 mm×607 mm×280 mm（长× 宽×高）。
图2 机器人流体仿真计算模型
Fig.2 Computational model for robot fluid simulation
1.3 网格划分
根据机器人尺寸划分外流场和密度盒内外 2 部分计算域。其中：外流场尺寸为 8L×6L×6L，密度盒尺寸为 1.3L×0.8L×0.5L，L 为机器人体长。计算模型放置于密度盒中，距离流体域速度入口的距离为 3L，尾部距离压力出口 4L，如图3 所示。
图3 流场范围及边界条件
Fig.3 Flow field range and boundary conditions
本文采用 ICEM 软件进行网格划分。外流场计算域尺寸较大，且无计算模型，故采用结构六面体网格进行划分；由于计算模型在部分工况中非静止，故密度盒内网格划分采用适应性强的四面体非结构网格，该网格支持 smooth 和 remesh 共 2 种动网格节点更新方法，在机器人计算模型表面进行加密，保证计算精度；2 个计算域间采用 interface 边界条件连接，并在接触面处进行加密，如图4 所示。
图4 计算域网格划分
Fig.4 Computational domain grid partitioning
1.4 仿真参数设定
机器人运动的外流场通常呈现大雷诺数湍流，因此采用可准确描述流场演化规律的 N-S 方程进行求解，笛卡尔坐标系下 RANS 方程如下：

\{\begin{matrix} \frac{\partial ρ}{\partial t} + \frac{\partial}{\partial x_{i}} (ρ u_{i}) = 0 \\ \frac{\partial}{\partial t} (ρ u_{i}) + \frac{\partial}{\partial x_{j}} (ρ u_{i} u_{j}) = - \frac{\partial ρ}{\partial x_{i}} + \\ \frac{\partial}{\partial x_{j}} [μ (\frac{\partial u_{i}}{\partial x_{j}} + \frac{\partial u_{j}}{\partial x_{i}} - \frac{2}{3} δ_{i j} \frac{\partial u_{i}}{\partial x_{j}})] + \frac{\partial}{\partial x_{j}} (- ρ_{j} \bar{u_{i}^{'} u^{'}}) \end{matrix}

(1)

式中： $ρ_{j} \bar{u_{i}^{'} u^{'}}$ 为雷诺应力； $μ$ 为动力粘性系数； $δ_{i j}$ 为克罗内克函数。
经无关性验证，确定密度盒内流场网格数量为 2640 000 和时间步长为 1/40 T 时满足精度要求，同时计算效率较高。
机器人仿真计算参数设定如表1 所示。
2 水动力系数求解与动力学建模
根据受力分析可知，若要建立机器人水下动力学模型，首先需要求解机器人在水下运动过程中受到的水动力^[6]。求解机器人水动力系数成为建立动力学模型的关键。出于研究周期和成本考虑，常采用 CFD 技术求解机器人不同运动状态下的水动力系数。
表1 机器人运动仿真参数
Table1 Robot motion simulation parameters
2.1 纯纵荡运动
进行纯纵荡运动仿真可求解机器人沿 x 轴方向速度 u 的惯性水动力系数 $X_{u}^{'}$ ，其运动方程如式（2）所示，根据运动方程编写流体仿真 UDF 文件，实现计算模型在流场内的运动。

\{\begin{matrix} U = 0, ω = 2 π f \\ x = A s i n (ω t) \\ u = \dot{x} = A ω c o s (ω t) \\ \dot{u} = - A ω^{2} s i n (ω t) \end{matrix}

(2)

式中：U 为流场流速，m/s；A 为机器人正弦运动的幅值，设置 A=40 mm；f 为机器人纯纵荡运动频率，取 f=0.3 Hz、f=0.4 Hz、f=0.5 Hz 共 3 组工况进行求解。
纯纵荡运动中机器人所受力的表达式为
$X = X_{\dot{u}} \dot{u} + X_{u} u + X_{0} =$
$- A ω^{2} X_{\dot{u}} s i n (ω t) + A ω X_{u} c o s (ω t) + X_{0} =$

X_{a} s i n (ω t) + X_{b} c o s (ω t) + X_{0}

(3)

计算收敛之后，机器人纯纵荡运动中受到的 x 轴方向的阻力 X 如图5 所示，由于数值仿真前期数据波动较大，可靠性低，因此本文选择 10~20 s 仿真时间内的数据进行拟合，此时计算已稳定。
通过曲线拟合，得到公式系数如表2 所示。
图5 纯纵荡运动中 x 轴受力曲线图
Fig.5 X-axis force curves in pure surge motion
由于机器人在水下运动过程中很少进行反向运动，若采用纯纵荡运动求解粘性水动力系数会造成较大误差，因此该工况仅用于求解惯性水动力系数，而粘性水动力系数由稳态直航运动进行求解。公式系数与水动力系数的关系式为

X_{a} = - A ω^{2} X_{\dot{u}} = - \frac{1}{2} ρ L^{3} A ω^{2} X_{\dot{u}}^{'}

(4)

通过建立公式系数 X_a和 $- 1 / 2 ρ L^{3} A ω^{2}$ 之间的线性关系求解纵荡运动的水动力系数，曲线的斜率即为水动力系数，而从小幅度振荡的线性理论可知，由多工况绘制的曲线应该具备过原点特征^[9]，因此建立的过原点曲线如图6 所示。
表2 机器人纯纵荡运动阻力曲线拟合结果
Table2 Fitting results of resistance curves for robot pure pitch motion
图6 过原点拟合 X 轴加速度系数曲线
Fig.6 A curve of fitting X-axis acceleration coefficient through the origin
由图6 可以看出，公式系数 X_a和 $- 1 / 2 ρ L^{3} A ω^{2}$ 之间整体呈线性关系，且数据拟合效果较好，证明了动态运动仿真求解机器人水动力系数的有效性，求得机器人纯纵荡运动无因次化水动力系数：

X_{\dot{u}}^{'} = - 0.02259

(5)

2.2 水动力系数求解汇总
以纯纵荡运动求解过程为例，可以求解与无因次化处理其它水动力系数。具体可参考文献^[10]。机器人水动力系数计算结果汇总如表3 所示。
表3 水动力系数计算结果汇总
Table3 Summary of hydrodynamic coefficient calculation results
2.3 五自由度动力学及运动学建模
结合机器人受力分析及水动力系数求解结果，可建立其动力学模型，忽略机器人横摇小扰动的影响^[7]，将机器人运动方程简化为五自由度动力学方程：

\{\begin{matrix} (m - X_{\dot{u}}) \dot{u} = m (v r - w q) + X_{u u} u^{2} + P_{1} + P_{2} \\ (m - Y_{\dot{v}}) \dot{v} - Y_{\dot{r}} \dot{r} = - m u r + Y_{v} v + Y_{r} r + Y_{v | v |} v | v | \\ (m - Z_{\dot{w}}) \dot{w} - Z_{\dot{q}} \dot{q} = m u q + Z_{w} w + Z_{w | w |} w | w | + Z_{q} q + Z_{δ} δ \\ (I_{y} - M_{\dot{q}}) \dot{q} - M_{\dot{w}} \dot{w} = M_{\dot{w}} w + M_{w | w |} w | w | + M_{q} q + \\ m g x_{g} c o s ϕ c o s θ + m g z_{g} s i n θ + z_{1} P_{1} + z_{2} P_{2} + M_{δ} δ \\ (I_{z} - N_{\dot{r}}) \dot{r} - N_{\dot{v}} \dot{v} = N_{v} v + N_{r} r + N_{v | v |} v | v | - \\ m g x_{g} s i n ϕ c o s θ - m g y_{g} s i n θ + y_{1} P_{1} - y_{2} P_{2} \end{matrix}

(6)

同理，可将机器人运动学方程进一步简化，得到其平动方程：

[\begin{matrix} x \\ y \\ z \end{matrix}] = [\begin{matrix} c o s θ c o s ψ & - s i n ψ & s i n θ c o s ψ \\ c o s θ s i n ψ & c o s ψ & s i n θ s i n ψ \\ - s i n θ & 0 & c o s θ \end{matrix}] [\begin{matrix} u \\ v \\ w \end{matrix}]

(7)

转动方程为

[\begin{matrix} \dot{θ} \\ \dot{ψ} \end{matrix}] = [\begin{matrix} 1 & 0 \\ 0 & 1 / c o s θ \end{matrix}] [\begin{matrix} q \\ r \end{matrix}]

(8)

基于数值仿真结果建立了机器人五自由度运动学和动力学模型，为机器人水下操纵性分析研究奠定了基础。
3 机器人水下操纵性运动特性研究
机器人水下操纵性运动特性分析是基于机器人运动仿真方程，通过设定初始运动参数及机器人推力，以数值积分的方式求解机器人下一时刻的运动状态，可以直观地看到水动力对机器人运动的影响。可以预报机器人运动，验证运动模型的准确性；同时可以总结机器人运动规律，为机器人的优化和研制提供指导。
3.1 运动仿真数值计算方法
为了保证仿真结果的准确性，机器人的仿真方程计算过程中，本文采用精度高、易收敛的四阶经典龙格库塔法进行积分。计算方法如下：

y_{i + 1} = y_{i} + \frac{h}{6} (k_{1} + 2 k_{2} + 2 k_{3} + k_{4})

(9)

其中：

\{\begin{matrix} k_{1} = f (t_{i}, y_{i}) \\ k_{2} = f (t_{i} + \frac{h}{2}, y_{i} + \frac{h}{2} k_{1}) \\ k_{3} = f (t_{i} + \frac{h}{2}, y_{i} + \frac{h}{2} k_{2}) \\ k_{4} = f (t_{i} + h, y_{i} + h k_{3}) \end{matrix}

(10)

3.2 机器人运动仿真结果及分析
基于操纵性仿真运动模型，本文针对机器人直航运动及水平面回转运动开展仿真研究，分析机器人运动状态，并对水动力系数和动力学模型的准确性进行验证。
3.2.1 机器人直航运动仿真
初始状态下，机器人的初始运动速度和姿态角均为 0，首先设置推进器推力为 10 N，接着每隔 30 N 进行一次运动仿真，求解不同推力下机器人航速和最大推力下的极限速度，结果如图7 所示。
图7 不同推力下机器人航速曲线
Fig.7 Robot speed curves under different thrusts
由图7 中机器人运动航速曲线可知：
1）当推进器推力 P=10 N 时，机器人航速 u=0.36 m/s；当 P=40 N 时， u =0.72 m/s； P = 70 N， u = 0.95 m/s；而推进器最大合力 106 N 左右，机器人最大航速约为 1.17 m/s。
2）机器人在螺旋桨推进力的作用下，由静止开始加速，且加速时间随着推进力的增大而缩短，最终速度趋于稳定。
3）当水动力与推力平衡时，机器人进入匀速行驶阶段。随着推力增大，机器人航速越大，阻力逐渐增大，而航速增长速率逐渐减小，符合机器人水下直航的一般规律。
3.2.2 机器人水平面回转运动
设置机器人左右推进器推力 P = 10 N，摇艏力矩 $N = - 3 N \cdot m$ ，机器人由静止开始运动，仿真时间设置为 70 s。其运动轨迹、运动过程中机器人在随体坐标系中的纵向速度 u、横向速度 v、垂向速度 w 的变化曲线如图8 所示。
图8 机器人水平面回转运动
Fig.8 Robot horizontal rotation motion
由图8 可知：
1）机器人在推力 P 和摇艏力矩 N 的作用下，由静止开始加速，在 10 s 之后开始进入匀速回转阶段。
2）当机器人开始匀速回转时，此时纵向速度 u = 0.32 m/s，横向速度 v = 0.035 m/s，回转半径 R 约为 3.095 m。
3）机器人在推进器推力和摇艏力矩的作用下可实现水下回转运动，符合水下机器人水平回转的一般规律。
结果表明，机器人水下运动基本符合一般水下机器人的运动规律，验证了水动力系数求解的准确性和五自由度运动方程的有效性，该方程可通过一系列仿真预报机器人运动状态，总结运动规律，为机器人闭环运动控制研究奠定了基础。
4 机器人水下操纵性运动实验
为了进一步研究机器人水下运动操纵性，本文开展水下操纵性实验，与数值仿真结果对比分析，从而验证机器人水动力系数和动力学模型的准确性。
4.1 直航运动实验
为了验证机器人直航运动仿真结果的准确性，本实验通过改变推进器推力进行机器人不同航速测试。机器人在推进器作用下从静止开始加速，加速完成后通过 10 m 标准航段，监测所需时长，从而计算出机器人运动速度，实验结果如下图9 所示。
图9 不同推力下机器人水下航行速度
Fig.9 Underwater navigation speed of robots under different thrusts
将实验结果与仿真数据对比分析，其结果如下表4 所示。
表4 不同推力下机器人航行速度
Table4 Robot navigation speed under different thrusts
由实验结果可知，机器人水下航行速度随着推力增大而不断增大，而阻力也显著提高，因此机器人速度增长速率逐渐减小，符合水下机器人的一般运动规律；且速度误差小于 6%，可见实验结果与仿真数据吻合较好，初步验证了机器人水动力系数和动力学模型的准确性，说明了利用动力学模型研究机器人水下运动特性的可行性。
4.2 水平面回转运动
为了进一步验证水动力系数求解和动力学模型的准确性，进行机器人水平面回转运动实验。实验通过 UWB（Ultra Wide Band，超宽带）高精度导航定位系统采集机器人的位置信息，如图10 所示。UWB 定位系统由基站、标签和控制台 3 部分组成，基站的 4 个 Link Track P-B 模块分别位于矩形实验区域四角；标签同样采用 Link Track P-B 模块，安装在机器人上，标签可测量到基站的距离并解算出坐标信息；控制台可与标签、基站通信，在上位机软件中实时显示机器人位置信息。UWB 定位系统的 2 维定位精度为 10 cm；当定位频率为 200 Hz 时，延迟为 5 ms；最远通信距离为 1 500 m。
实验过程中，初始条件设置机器人推进力合力为 10 N，接着基于初始条件分别设置两侧的推进器推力，在满足推力要求的同时使机器人摇艏力矩为 3 N·m。实验过程如图11 所示，机器人从静止加速进行水平面回转运动。
实验数据由于噪声干扰出现小范围波动，因此将实验数据通过局部二次回归进行平滑处理，机器人水平面回转运动的仿真数据、实验数据及拟合数据如图12 所示。
由于游泳池水流循环系统的影响，实验场地并非静水环境，存在随机水流的干扰。由运动轨迹可知，在第 2 圈回转运动过程中，机器人于 1/4 航程处偏航角发生了明显变化，导致后续运动轨迹发生了偏转，因实验环境与仿真环境一致性不高，误差较大。故以第 1 圈机器人轨迹进行数据分析，在推进力和摇艏力矩作用下机器人进行水平面回转运动，运动轨迹近似为圆，与仿真数据重合度较高，规律一致，最大偏移量为 0.15 m，偏差小于 5%，满足精度要求。实验结果表明机器人水平面回转运动特性与仿真结果相符，进一步验证了动力学模型的准确性。
图10 UWB 定位系统原理示意图
Fig.10 Schematic diagram of UWB positioning system principle
图11 机器人水平面回转运动实验
Fig.11 Experiments on robot horizontal rotation motion
图12 机器人水平面回转运动轨迹
Fig.12 Robot horizontal rotation motion trajectories
5 结束语
仿生波动鳍水陆两栖机器人凭借隐蔽性强、机动性好等优点成为海洋设备的研究热点，然而由于现阶段仿生机器人仿生程度低、推进性能差、作业环境复杂等客观因素，使其在执行任务中面临着推进速度慢、效率低等困境。针对上述问题，本文围绕新型水陆两栖机器人水下动力学建模、水下操纵性数值仿真及实验验证等方面开展研究，仿真结果与实验结果误差均不超过 10%，验证了动力学建模的准确性。相关成果为进一步研究多模态水陆两栖机器人运动性能奠定了理论与技术基础。
参考文献
- [1] WU Z Y，QI J，ZHANG S.Amphibious robots：a review[J].Applied Mechanics and Materials，2014，494-495：1036-1041.
- [2] RAFEEQ M，TOHA S F，AHMAD S，et al.Locomotion strategies for amphibious robots-a review[J].IEEE Access，2021，9：26323-26342.
- [3] FRANCESCHI A，PIAGGIO B，TONELLI R，et al.Assessment of the manoeuvrability characteristics of a twin shaft naval vessel using an open-source CFD code[J].Journal of Marine Science and Engineering，2021，9（6）：665.
- [4] GO G，AHN H T.Hydrodynamic derivative determination based on CFD and motion simulation for a tow-fish[J].Applied Ocean Research，2019，82：191-209.
- [5] DUBBIOSO G，BROGLIA R，ZAGHI S.CFD analysis of turning abilities of a submarine model[J].Ocean engineering，2017，129：459-479.
- [6] SHEN Z R，WAN D C，CARRICA P M.Dynamic overset grids in OpenFOAM with application to KCS self-propulsion and maneuvering[J].Ocean Engineering，2015，108：287-306.
- [7] 尹盛林.仿生波动鳍水陆两栖机器人设计与控制技术研究[D].西安：西安交通大学，2021.
- [8] 王雷.水下机器人运动控制研究[D].合肥：中国科学技术大学，2018.
- [9] 胡芳琳，冯康佳，马向能.动态PMM实验的数值模拟及线性导数获取方法[J].水动力学研究与进展：A 辑，2021，36（2）：9.
- [10] 施生达.潜艇操纵性[M].北京：国防工业出版社，1995.
- [11] FOSSEN T I.Guidance and control of ocean vehicles[M].New York：John Willey & Sons，1994.

基本信息

中图分类号: TP242.6
文献标识码: A
DOI: 10.19838/j.issn.2096-5753.2023.04.012
文章编号: 2096-5753(2023)04-0487-09

基金信息

成果转化项目“水陆两栖轮鳍式机器人技术”（2016330ZD01000208）；

引用信息

陈龙，童保成，祖斯羽，等. 水陆两栖机器人水下动力学建模及操纵性研究[J]. 数字海洋与水下攻防，2023， 6（4）：487-495.

稿件历史

收稿日期: 2023-05-25

参考文献

[1] WU Z Y，QI J，ZHANG S.Amphibious robots：a review[J].Applied Mechanics and Materials，2014，494-495：1036-1041.
[2] RAFEEQ M，TOHA S F，AHMAD S，et al.Locomotion strategies for amphibious robots-a review[J].IEEE Access，2021，9：26323-26342.
[3] FRANCESCHI A，PIAGGIO B，TONELLI R，et al.Assessment of the manoeuvrability characteristics of a twin shaft naval vessel using an open-source CFD code[J].Journal of Marine Science and Engineering，2021，9（6）：665.
[4] GO G，AHN H T.Hydrodynamic derivative determination based on CFD and motion simulation for a tow-fish[J].Applied Ocean Research，2019，82：191-209.
[5] DUBBIOSO G，BROGLIA R，ZAGHI S.CFD analysis of turning abilities of a submarine model[J].Ocean engineering，2017，129：459-479.
[6] SHEN Z R，WAN D C，CARRICA P M.Dynamic overset grids in OpenFOAM with application to KCS self-propulsion and maneuvering[J].Ocean Engineering，2015，108：287-306.
[7] 尹盛林.仿生波动鳍水陆两栖机器人设计与控制技术研究[D].西安：西安交通大学，2021.
[8] 王雷.水下机器人运动控制研究[D].合肥：中国科学技术大学，2018.
[9] 胡芳琳，冯康佳，马向能.动态PMM实验的数值模拟及线性导数获取方法[J].水动力学研究与进展：A 辑，2021，36（2）：9.
[10] 施生达.潜艇操纵性[M].北京：国防工业出版社，1995.
[11] FOSSEN T I.Guidance and control of ocean vehicles[M].New York：John Willey & Sons，1994.

首页

本刊简介

编委会

收录情况

期刊订阅

征稿简则

联系我们

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

水陆两栖机器人水下动力学建模及操纵性研究

Research on Underwater Dynamics Modeling and Maneuverability of Amphibious Robots

摘要

Abstract

关键词

Keywords

0 引言

1 机器人坐标系建立及数值求解方法

1.1 坐标系建立

1.2 计算模型

1.3 网格划分

1.4 仿真参数设定

(1)

2 水动力系数求解与动力学建模

2.1 纯纵荡运动

(2)

(3)

(4)

(5)

2.2 水动力系数求解汇总

2.3 五自由度动力学及运动学建模

(6)

(7)

(8)

3 机器人水下操纵性运动特性研究

3.1 运动仿真数值计算方法

(9)

(10)

3.2 机器人运动仿真结果及分析

3.2.1 机器人直航运动仿真

3.2.2 机器人水平面回转运动

4 机器人水下操纵性运动实验

4.1 直航运动实验

4.2 水平面回转运动

5 结束语

参考文献

基本信息

基金信息

引用信息

稿件历史

参考文献