CN113885508B - 水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法 - Google Patents
水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法 Download PDFInfo
- Publication number
- CN113885508B CN113885508B CN202111224379.0A CN202111224379A CN113885508B CN 113885508 B CN113885508 B CN 113885508B CN 202111224379 A CN202111224379 A CN 202111224379A CN 113885508 B CN113885508 B CN 113885508B
- Authority
- CN
- China
- Prior art keywords
- underwater robot
- moment
- coordinate system
- force
- body coordinate
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/02—Control of position or course in two dimensions
- G05D1/0206—Control of position or course in two dimensions specially adapted to water vehicles
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Abstract
本发明提供了一种水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法,其根据受到的流体惯性力、流体位置力、流体阻尼力、浮力、重力、推进器推力、刀锋腿运动产生的水动力及崎岖面对支撑腿的反作用力,建立混合驱动机器人的动力学模型,构建支撑腿不打滑的不等式约束条件。选取推进器的驱动力和力矩作为六维状态量,并在六维空间内确定若干个探索方向,利用最优化方法探寻在保证水下机器人在崎岖面稳定攀爬下的各方向许用的最大驱动力和力矩。此外,该发明利用反向传播神经网络来训练不同运动状态下各方向的最大许用驱动力和力矩,大幅度减小探索水下机器人稳定空间的时间。最后,获得最佳的推进器推力,以增强混合驱动水下机器人的稳定性。
Description
技术领域
本发明涉及水下机器人技术领域,具体为一种多推进器和六条半圆形刀锋腿混合驱动水下机器人在崎岖面攀爬的稳定空间生成及稳定性增强方法。
背景技术
水下机器人是海洋资源开发、海洋监测以及海洋生态保护的重要装备之一。通过搭载不同类型的传感器和执行器,水下机器人能够有效地实现对海洋探索、开发、监测以及侦查等多项任务。
当水下机器人在不平整工作面执行攀爬或行走等作业任务时,就需要确保水下机器人具有一定的越障能力。腿式机器人作为一种具有卓越越障能力的机器人受到国内外学者的广泛关注。因此,研究一类由多推进器和六条半圆形刀锋腿混合驱动且能在不平整工作面攀爬或行走的水下机器人具有重要的经济和社会效益。
目前,多足机器人常见的稳定性判据方法有基于零力矩点(ZMP)的稳定性判据方法,其本质是多足机器人不发生欠驱动翻转运动的约束条件。但该方法仅能用于表征在平坦地面行走的多足机器人稳定性,且需要假设支撑腿不与地面打滑。Kensuke Harada等人在《IEEE TRANSACTIONS ON ROBOTICS》(2006年第3期)发表文章《Dynamics and Balanceof a Humanoid Robot During Manipulation Tasks》,在ZMP稳定性判据的基础上,提出了一种普适性更好的稳定性判据方法,通过分析带操作臂的人形机器人所受的重力、惯性力以及重心处的力矩是否同时在由多足摩擦力所组成的多面体内来判断机器人是否稳定。这种方法考虑了支撑腿与壁面之间的摩擦约束,且能用于判断陆上多足机器人在崎岖壁面攀爬的稳定性。CNRS-AIST联合实验室的Herve Audren等在《IEEE TRANSACTIONS ONROBOTICS》(2018年第2期)发表文章《3-D Robust Stability Polyhedron inMulticontact》研究了当腿式机器人静态支撑于崎岖面时稳定空间的构建问题。在分析机器人本体静态力学约束和支撑腿摩擦力约束的基础上,该文章利用搜索策略,提出了一种多足机器人稳定质心域的生成方法。
发明内容
(1)技术问题
现有方法仅适用于生成由单一驱动机构(腿)驱动的陆上机器人在崎岖面攀爬的稳定空间,不能用于生成混合驱动水下机器人的稳定空间。此外,现有基于搜索策略的腿式机器人稳定空间生成方法仅适用于静态支撑的腿式机器人。最后,现有方法还存在难以兼顾质心域生成准确度高和生成时耗少的问题。
不同于陆上多足机器人,本发明所研究的混合驱动水下机器人由6条刀锋腿和8个推进器共同驱动。这种新型驱动方式会产生新的稳定空间生成方法和稳定性增强方法。具体地,在生成这类机器人稳定空间并增强其稳定性方面,主要存在以下难点:1、如何融入崎岖轮廓的信息和腿滑移的约束表征这类混合驱动水下机器人在崎岖面攀爬的稳定空间;2、如何设计合适的策略以提高稳定空间生成的准确性和减少稳定空间生成的时耗;3、如何设计合适的推进器推力大小,来增强这类水下机器人攀爬崎岖面的稳定性。
(2)技术方案
为了现有技术的不足,本发明提出了一种水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法,包括如下步骤:
步骤1:在分析水下机器人本体受到的流体惯性力、流体位置力、流体阻尼力、浮力、重力、推进器推力、刀锋腿运动产生的水动力及崎岖面对支撑腿的反作用力的基础上,建立水下机器人的动力学模型,同时,构建支撑腿不打滑的不等式约束条件;
步骤2:选取推进器的驱动力和力矩作为六维状态量,并在六维空间内确定若干个搜索方向,利用最优化方法探寻水下机器人在保证稳定攀爬下的各方向许用的最大驱动力和力矩,进而获取水下机器人在崎岖面攀爬的稳定空间;
步骤3:利用BP神经网络训练不同运动状态下各方向的最大许用驱动力和驱动力矩,从而生成水下机器人的稳定空间;
步骤4:根据水下机器人当前的驱动力、驱动力矩和构建的稳定空间,设计最佳的推进器推力,以增强混合驱动水下机器人在崎岖面攀爬的稳定性。
根据本发明的实施例,步骤1中建立混合驱动水下机器人在体坐标系下的动力学模型,其可表示为
其中,为体坐标系下水下机器人浮心的速度,/>为体坐标系下水下机器人的旋转角速度,/>为体坐标系下作用于水下机器人上的全部外力,/>为体坐标系下作用于水下机器人上的全部外力矩,水下机器人的动量/>和动量矩/>可表示为
其中,为在体坐标系下水下机器人质心的速度,m为水下机器人的质量,且
vc=vo+ω×rc
其中,为水下机器人质心至浮心的矢径,它在体坐标系中三个轴上的分量为xc,yc、zc,/>为水下机器人相对浮心的转动惯量矩阵,且定义为
同时,作用于水下机器人的全部外力F和外力矩M可分别表示为
其中,为体坐标系下水下机器人所受的流体惯性力,/>为体坐标系下水下机器人所受的流体位置力,/>为体坐标系下水下机器人所受的阻尼力,和/>分别为体坐标系下水下机器人所受的浮力和重力,/>为体坐标系下水下机器人受到的推进器推力,/>为体坐标系下水下机器人腿运动所产生的力,为体坐标系下水下机器人受到的支撑壁面对支撑腿的反作用力;/>为体坐标系下水下机器人受到的流体惯性力矩,/>为体坐标系下水下机器人受到的流体位置力矩,/>为体坐标系下水下机器人受到的阻尼力矩,/>为体坐标系下水下机器人受到的重力矩,/>为体坐标系下水下机器人受到的推进器推力所形成的力矩,为体坐标系下腿运动所产生的力矩,/>为体坐标系下水下机器人受到的支撑壁面对支撑腿的反作用力所形成的力矩。
根据本发明的实施例,步骤1中支撑腿不打滑约束可以表示为
其中,μ为刀锋腿与崎岖壁面之间的摩擦系数,和/>分别为第i(i=1,2,…6)个支撑腿沿崎岖壁面前向、垂向和横向的力。
根据本发明的实施例,步骤2中设定基于推进器驱动力F3和驱动力矩M3的六维探索状态量,同时,设置N个探索方向其中,l=1,…,N;
在上述动力学模型等式约束和支撑腿不打滑的不等式约束下,利用最优化方法最大化其中,/>可获取在当前运动状态下第l个方向上的最大许用驱动力和驱动力矩;
通过解算N个探索方向上的最大许用驱动力和力矩,便可获取当前运动状态下混合驱动水下机器人在崎岖面攀爬的稳定空间。
根据本发明的实施例,步骤3中通过收集大量不同运动状态下的稳定空间数据和状态数据,来训练运动状态量与稳定控制之间的映射关系;具体地,采用BP神经网络来表征这种映射关系;当参数训练完成时,输入水下机器人的运动状态,即可快速生成混合驱动水下机器人的稳定空间。
根据本发明的实施例,步骤4中水下机器人第j(j=1,2,…,M)组运动状态下的第l个方向上的最大许用驱动力和力矩矢量可定义为:
其中,M为水下机器人运动状态的总量,为第j组运动状态下的第l个方向上的驱动力和力矩,经过矢量/>末端且与矢量/>垂直的6维超平面方程可定义为:
其中,表示水下机器人六维状态空间的第k维状态量,k=1,…,6;
稳定空间内的任一驱动力和力矩到六维超平面的距离/>可表示为
其中,||·||为矢量的范数符号;
利用最优化方法寻找出第j组运动状态下Dj最大的最优驱动力F3和力矩M3,以使水下机器人的稳定裕度最大;
根据最优驱动力F3和力矩M3求解8个推进器的推力值。
(3)有益效果
(1)提供了一种在崎岖面攀爬的混合驱动水下机器人稳定空间生成和稳定性增强方法,其可以大幅度减小探索水下机器人稳定空间的时间,并增强混合驱动水下机器人的稳定性;
(2)为解释混合驱动水下机器人的复杂运动提供了新工具,也为混合驱动水下机器人稳定控制器的设计和攀爬步态的规划提供了一种新思路。
附图说明
图1为一种在崎岖面攀爬的混合驱动水下机器人稳定空间生成和稳定性增强方法示意图。
图2为混合驱动水下机器人定义的地面坐标系和体坐标系。
图3为混合驱动水下机器人六条腿髋关节位置示意图。
图4为混合驱动水下机器人八个推进器安装位置示意图。
图5为第i条刀锋腿在崎岖工作面上攀爬的示意图。
图6为第i条刀锋腿绕髋关节转动的示意图。
图7为利用正八边形近似代替圆形的示意图。
具体实施方式
现结合实施例、附图对本发明的具体实施例进行详细描述:
首先,对附图进行阐述和说明:
图1为一种混合驱动水下机器人在崎岖面攀爬的稳定空间生成和稳定性增强方法示意图。本发明在获取机器人本体的水动力(矩)、重浮力(矩)、惯性力(矩)、推进器推力(矩)、刀锋腿运动所产生的水动力(矩)和刀锋腿与崎岖面作用力(矩)的基础上,考虑支撑腿不打滑约束,选取推进器的驱动力和驱动力矩为待搜索的状态量,确定搜索方向,探寻水下机器人的稳定空间。同时,针对这类稳定空间搜索方法耗时过长的问题,本发明利用BP神经网络训练不同运动状态下稳定空间,其能大幅度减小稳定空间的搜索时间,提升了该稳定性判据方法的简便性。最后,设计合适的推进器推力大小来增加水下机器人稳定性。
图2为混合驱动水下机器人的体坐标系Oxyz和地面坐标系O0x0y0z0示意图。图中,体坐标系Oxyz与机器人本体固连,即相对机器人本体静止不动的坐标系。体坐标系的原点选在机器人本体的浮心;Ox轴沿机器人本体纵轴并指向前方;Oy轴垂直于Ox轴并指向上方;Oz轴垂直于轴Ox和Oy,其方向满足右手坐标系定义法则。
地面坐标系O0x0y0z0与大地固连。地坐标系原点O0选在地面的某处,O0y0轴铅直向上,O0x0轴在水平面内指向任一方向,O0z0轴垂直于轴O0x0和O0y0,其方向使该坐标系成为右手坐标系。
图3为混合驱动水下机器人六条腿髋关节位置示意图。图中,六条腿均匀布置于机器人本体的两侧,且腿3与腿4的连线经过机器人本体的中点并与体坐标系中的Ox轴垂直。同时,d1为腿1与腿2(腿3与腿4、腿5与腿6)之间的距离,d2为腿1与腿3(腿3与腿5、腿2与腿4、腿4与腿6)之间的距离。
图4为混合驱动水下机器人八个推进器安装位置示意图。图中,4个垂直推进器垂直布置,4个水平推进器以45°角矢量布置,这8个推进器皆内嵌于机器人本体中。l1为垂直推进器1与垂直推进器2(垂直推进器3与垂直推进器4)之间的距离,l2为垂直推进器1与垂直推进器4(垂直推进器2与垂直推进器3)之间的距离,l3为水平推进器5与水平推进器6(水平推进器7与水平推进器8)之间的距离,l4为水平推进器5与水平推进器8(水平推进器6与水平推进器7)之间的距离。本发明所指的水下机器人即指混合驱动水下机器人,或者具有多推进器和多条半圆形刀锋腿混合驱动水下机器人。
图5为第i条刀锋腿在崎岖工作面上攀爬的示意图。如图所示,θ为机器人本体的俯仰角,Ci和Di分别为第i条刀锋腿的髋关节和触地点,qi为第i条刀锋腿的旋转角度,r为半圆形刀锋腿的半径,γi为第i条刀锋腿的支撑角,和/>分别为第i个支撑腿沿壁面前向、垂向和横向的力;腿坐标系Cixiyizi与第i条刀锋腿固连,其坐标系原点为第i条腿髋关节Ci,Cixi轴指向第i条腿的另一个末端;Cixi轴与Ciyi轴构成的平面与刀锋腿所构成的平面平行,且Cixi轴逆时针旋转90°的方向为轴Ciyi的方向;Cizi轴垂直于轴Cixi和Ciyi,其方向满足右手坐标系定义法则。
以下将结合附图详细介绍本发明的具体实施方法。
步骤1:在分析混合驱动水下机器人本体受到的流体惯性力(矩)、流体位置力(矩)、流体阻尼力(矩)、浮力(矩)、重力(矩)、推进器推力(矩)、刀锋腿运动产生的水动力(矩)及崎岖面对支撑腿的反作用力(矩)的基础上,建立了混合驱动机器人的动力学模型,同时构建了支撑腿不打滑的不等式约束条件;
步骤2:选取推进器的驱动力和力矩作为六维状态量,并在六维空间内确定若干个搜索方向,利用最优化方法探寻混合驱动水下机器人在崎岖面稳定攀爬下的各方向许用的最大驱动力和力矩,进而生成混合驱动水下机器人在崎岖面攀爬的稳定空间;
步骤3:利用BP神经网络训练不同运动状态下各方向的最大许用驱动力和驱动力矩,即生成混合驱动水下机器人的稳定空间。
步骤4:根据水下机器人当前的驱动力、驱动力矩和稳定空间,设计最佳的推进器推力,以增强混合驱动水下机器人在崎岖面攀爬的稳定性。
混合驱动水下机器人在体坐标系下的动力学模型可表示为
其中,为体坐标系下水下机器人浮心的速度,/>为体坐标系下水下机器人的旋转角速度,/>为体坐标系下作用于水下机器人上的全部外力,/>为体坐标系下作用于水下机器人上的全部外力矩,水下机器人的动量/>和动量矩/>可表示为
其中,为在体坐标系下水下机器人质心的速度,m为水下机器人的质量,且
vc=vo+ω×rc (3)
其中,为水下机器人质心至浮心的矢径,它在体坐标系中三个轴上的分量为xc,yc、zc。/>为水下机器人相对浮心的转动惯量矩阵,其定义为
此外,作用于水下机器人的全部外力F和外力矩M可分别表示为
其中,为体坐标系下水下机器人所受的流体惯性力,/>为体坐标系下水下机器人所受的流体位置力,/>为体坐标系下水下机器人所受的阻尼力,和/>分别为体坐标系下水下机器人所受的浮力和重力,/>为体坐标系下水下机器人受到的推进器推力,/>为体坐标系下水下机器人腿运动所产生的力,为体坐标系下水下机器人受到的支撑壁面对支撑腿的反作用力;/>为体坐标系下水下机器人受到的流体惯性力矩,/>为体坐标系下水下机器人受到的流体位置力矩,/>为体坐标系下水下机器人受到的阻尼力矩,/>为体坐标系下水下机器人受到的重力矩,/>为体坐标系下水下机器人受到的推进器推力所形成的力矩,为体坐标系下腿运动所产生的力矩,/>为体坐标系下水下机器人受到的支撑壁面对支撑腿的反作用力所形成的力矩。
将式(2)、(3)和(4)代入式(1)中,可得
接下来,分别介绍作用于机器人的各类外力与外力矩。
①体坐标系下的流体惯性力位置力/>和阻尼力/>为
其中,Xuu、Yvv和Zww分别为相对vox、voy和voz的阻力因数、升力因数和侧力因数,/>为机器人升力因数对角速度wz的旋转导数,/>为机器人侧力因数对角速度wy的旋转导数。
②体坐标系下浮力和重力/>为
其中,B为机器人的浮力,G=mg为机器人的重力,g为重力加速度,θ为水下机器人本体的俯仰角,为水下机器人本体的横滚角。
③体坐标系下八个推进器的推力F3为
F3=C1T
其中,为推力矢量,Ti为第i个推进器推力大小,C1为8个推进器推力分量至体坐标系三个轴方向力的矩阵,其可表示为
④体坐标系下刀锋腿运动(相对地坐标系平动+相对髋关节转动)所产生的水动力F4为
其中,和/>分别为第i条刀锋腿运动所产生的水动力在体坐标系Ox和Oy方向的分量。
如图6所示,任取第i条刀锋腿上的一小段圆弧S。该小段圆弧的速度可表示为
其中,αi为直线Cixi绕Ci旋转至直线CiSi的旋转角度,qi为第i条刀锋腿的旋转角度,为第i条刀锋腿的旋转角速度,r为刀锋腿的半径。
速度沿腿坐标系Cixi和Ciyi方向的分量为
髋关节固连于机器人本体,且在三维空间内运动。当刀锋腿绕髋关节旋动时,需考虑髋关节运动对刀锋腿内任意位置处的速度的影响。刀锋腿任意微小圆弧S受到的水阻力dfxi(沿腿坐标系Cixi方向)和dfyi(沿腿坐标系Ciyi方向)可表示为
其中,W为刀锋腿的宽度,ρ为水的密度,μc为刀锋腿的阻力系数,为第i个髋关节Ci沿腿坐标系Cixi和Ciyi方向的速度,其可表示为
其中,和/>为第i条刀锋腿髋关节Ci沿体坐标系Ox和Oy方向的速度,为第i条刀锋腿髋关节相对体坐标系的速度矢量,/> 为第i条腿髋关节至水下机器人浮心的矢径。
将式(6)代入至和/>中,可得
由于式(8)难以用于直接判断的正负号,将导致不能通过直接积分的方式获得刀锋腿的水动力的合力。考虑到刀锋腿在绝大多数情况下是正转,因此,本发明假设在此基础上,令/>可得/>
当时,表明/>始终小于零。当/>时,/>可能大于零也可能小于零。为了描述在这种情况下/>的正负号,我们定义了两个角度当/>时,表明/>小于零;反之,当/>或/>时,表明/>大于零。当/>时,表明/>始终大于零。
在上述分析的基础上,可推导出刀锋腿沿腿坐标系Cixi方向所受的水动力fxi为
其中,
将式(8)代入至式(9),求积分,可得
同理,在定义角度的基础上,可推导出刀锋腿沿腿坐标系Ciyi方向所受的水动力fyi为
其中,
将式(8)代入至式(11),求积分,可得
第i条刀锋腿所受的水动力在体坐标系Oy和Oy方向的分量和/>为
⑤体坐标系下壁面对刀锋腿的反作用力F5为
其中,壁面对第i条刀锋腿的反作用力在体坐标系中的各分量为
其中,和/>分别为第i个腿沿壁面前向、垂向和横向的力,γi∈[0,π]为第i条腿的触地角,其具体定义如图5所示。
和/>分别为第i个支撑腿沿壁面前向、垂向和横向的力,ai为第i条腿的触地状态。当ai=0时,第i条半圆形刀锋腿未触地;当ai=1时,第i条半圆形刀锋腿触地。
接下来,介绍触地状态ai和触地角γi的获取方法。利用二元函数y=f(x,z)来表征崎岖地面的轮廓,即在O0x0z0平面内的任一坐标位置(x,z)处,其崎岖轮廓的高度为y。第i条半圆形刀锋腿上的任一坐标值可定义为将该坐标值代入轮廓函数/>当存在γi∈[0,π]使/>时,表明第i条半圆形刀锋腿与崎岖地形接触,即第i条刀锋腿为支撑腿,同时,第i条半圆形刀锋腿的触地状态ai=1;当不存在γi∈[0,π]使时,表明第i条半圆形刀锋腿未与崎岖地形接触,即这条刀锋腿为摆动腿,同时,第i条半圆形刀锋腿的触地状态ai=0。
值得注意的是,由于本发明假设崎岖工作面的曲率半径始终大于刀锋腿半径r,因此,第i条腿的触地角γi不存在多解。
⑥体坐标系下水下机器人所受的流体惯性力矩位置力矩/>和阻尼力矩/>分别为
其中,λ44、λ55和λ66为绕体坐标系Ox、Oy和Oz方向的附加转动惯量,和/>分别为水下机器人的横滚力矩、偏航力矩和俯仰力矩因数对速度voz、voz和vox的位置导数,和/>分别为水下机器人的横滚力矩、偏航力矩和俯仰力矩因数对角速度wx、wy和wz的旋转导数。
⑦体坐标系下重力矩M2和推进器推力所形成的力矩M3可分别表示为
其中,C2为8个推机器你生成体坐标系三轴方向的合力矩的矩阵,其可表示为
其中,l3为推进器5和6(或7和8)之间的距离,l4为推进器6和7(或5和8)之间的距离。
⑧体坐标下由于刀锋腿运动所产生的力矩M4可表示为
其中,C4和C5分别为力矩分配矩阵其可表示为
其中,为第i条绕髋关节转动时所产生的力矩。
将第i条刀锋腿上的部分圆弧S处所受到力元投影至Cixi和Ciyi方向,分别乘以相应的力臂,积分即可得到绕体坐标系Oz方向的力矩为
/>
对式(13)求积分,可得
由于刀锋腿转动时不会产生绕体坐标系Ox和Oy方向的力矩,即
综上,第i条刀锋腿绕髋关节转动时所产生的力矩可表示为
⑨体坐标系下壁面对刀锋腿的反作用力所形成的力矩M5为
支撑腿不打滑约束可以表示为
设定基于驱动力F3和驱动力矩M3的六维探索状态量,同时,设置N个探索方向其中,l=1,…,N。在上述动力学模型等式约束和支撑腿不打滑的不等式约束下,利用诸如机器人操作系统(ROS)中的scipy.opt imize库中最优化方法最大化/>其中,可获取在当前运动状态下第l个向量上驱动力和驱动力矩的最大许用值/>其中,下标j表示机器人的第j组运动状态。求N个探索向量上的驱动力和驱动力矩的最大许用值,即可获取当前运动状态下混合驱动水下机器人的稳定空间。
需要说明的是,水下机器人真实的稳定空间的求解代价是极其昂贵的。这将导致无法实时地获取机器人的真实稳定空间,进而导致无法用于实际的机器人系统。因此,本发明沿着有限个探索方向来探索水下机器人的最大许用驱动力和驱动力矩,进而生成水下机器人的近似稳定空间。
如图7所示,本发明利用二维平面内的正八边形近似代替圆形区域的方式来形象地描述上述近似方法,其中,正八边形代表近似的稳定空间,圆形表示真实的稳定空间。具体地,在圆形区域内均匀地设置八个探索方向,其中,相邻的两个探索方向之间的夹角为45°。在此基础上,以圆心为起点,沿着八个探索方向探索与圆形区域(真实稳定空间)边界的八个交点(最大许用驱动力和驱动力矩)。这八个交点的切线可形成一个正八边形(近似稳定空间),其可用于近似代替圆形区域,即近似替代真实的稳定空间。
如图7所示,设置的八个探索方向分别定义为
其中,4个探索方向与两个坐标轴的正负方向保持平行,另外的4个探索方向均布于坐标轴之间。
受二维平面内探索方向设置方法的启发,针对本发明中的六维稳定空间,设置了76个探索方向,其中,6*2=12个探索方向分别与6条坐标轴正负方向保持平行,2^6=64个探索方向均布于6个坐标轴之间,其分别定义为
值得注意的是这76个探索方向均设置为单位向量。
针对无法实时地获取稳定空间的问题,本发明提出通过利用BP神经网络离线表征运动状态与稳定空间之间的映射关系的方法。通过多次的水下机器人攀爬崎岖面的实验,收集M组运动状态数据集ΩD={S1,…,SM},其中,第j组运动状态定义为M组运动状态下的76个稳定空间探索向量上的最大驱动力和驱动力矩许用值数据集/>其中,/>表示第j(j=1,…,M)组运动状态下76个稳定空间探索方向上的最大许用驱动力和驱动力矩。接下来,利用BP神经网络来离线训练数据集ΩD和数据集Ωζ之间的映射关系,即获取运动状态与稳定空间之间的映射关系。当参数训练完成时,输入水下机器人的任一组运动状态Sj,即可快速生成混合驱动水下机器人的稳定空间/>
混合驱动水下机器人第j组运动状态下的第l个方向上的最大许用驱动力和力矩矢量可定义为经过矢量/>末端且与矢量/>垂直的6维超平面方程可定义为
/>
其中,表示水下机器人六维状态空间第k(k=1,…,6)维的状态量。
稳定空间内的任一驱动力和力矩到式(6)中的六维超平面的距离/>可表示为
其中,||·||为矢量的范数符号。
利用诸如机器人操作系统(ROS)中的scipy.optimize库中最优化方法寻找出第j组运动状态下使Dj最大的最优驱动力F3和力矩M3,其能使混合驱动水下机器人的稳定裕度最大。根据最优驱动力F3和力矩M3可解算出8个推进器的推力值。
具体地,8个推进器的推力T可表示为
其中,为推力矢量,Ti为第i个推进器推力大小,pinv(C)表示矩阵C的伪逆矩阵,C为将8个推进器推力映射至体坐标系三个轴方向力和力矩的矩阵,其可表示为
其中,l1、l2、l3和l4的定义如图4所示。
这些推力值能大幅增强混合驱动水下机器人攀爬崎岖面的稳定性。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,虽然本发明已以较佳实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本发明的技术实质对以上实施例所作任何的简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。
Claims (5)
1.一种水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法,水下机器人由8个推进器、6条刀锋腿和机器人本体构成;其中,4个垂向推进器垂直布置,4个水平推进器以45°角矢量布置,六条半圆形刀锋腿的髋关节与机器人本体左右两侧的六个驱动电机连接,其特征在于,所述稳定空间生成和稳定性增强方法包括如下步骤:
步骤1:在分析水下机器人本体受到的流体惯性力、流体位置力、流体阻尼力、浮力、重力、推进器推力、刀锋腿运动产生的水动力及崎岖面对支撑腿的反作用力的基础上,建立水下机器人的动力学模型,同时,构建支撑腿不打滑的不等式约束条件;
步骤2:选取推进器的驱动力和力矩作为六维状态量,并在六维空间内确定若干个搜索方向,利用最优化方法探寻水下机器人在保证稳定攀爬下的各方向许用的最大驱动力和力矩,进而获取水下机器人在崎岖面攀爬的稳定空间;
步骤3:利用BP神经网络训练不同运动状态下各方向的最大许用驱动力和驱动力矩,从而生成水下机器人的稳定空间;
步骤4:根据水下机器人当前的驱动力、驱动力矩和构建的稳定空间,得到推进器推力,以增强混合驱动水下机器人在崎岖面攀爬的稳定性;
步骤1中建立水下机器人的动力学模型包括:
水下机器人在体坐标系下的动力学模型可表示为
其中,为体坐标系下水下机器人浮心的速度,/>为体坐标系下水下机器人的旋转角速度,/>为体坐标系下作用于水下机器人上的全部外力,为体坐标系下作用于水下机器人上的全部外力矩,水下机器人的动量/>和动量矩/>可表示为
其中,为在体坐标系下水下机器人质心的速度,m为水下机器人的质量,且
vc=vo+ω×rc
其中,为水下机器人质心至浮心的矢径,其在体坐标系中三个轴上的分量为xc,yc、zc,/>为水下机器人相对浮心的转动惯量矩阵;
同时,作用于水下机器人上的全部外力和外力矩/>可分别表示为
其中,为体坐标系下水下机器人所受的流体惯性力,/>为体坐标系下水下机器人所受的流体位置力,/>为体坐标系下水下机器人所受的阻尼力,/>和分别为体坐标系下水下机器人所受的浮力和重力,/>为体坐标系下水下机器人受到的推进器推力,/>为体坐标系下水下机器人腿运动所产生的力,/>为体坐标系下水下机器人受到的支撑壁面对支撑腿的反作用力;/>为体坐标系下水下机器人受到的流体惯性力矩,/>为体坐标系下水下机器人受到的流体位置力矩,/>为体坐标系下水下机器人受到的阻尼力矩,/>为体坐标系下水下机器人受到的重力矩,/>为体坐标系下水下机器人受到的推进器推力所形成的力矩,/>为体坐标系下腿运动所产生的力矩,/>为体坐标系下水下机器人受到的支撑壁面对支撑腿的反作用力所形成的力矩。
2.根据权利要求1所述的稳定空间生成和稳定性增强方法,其特征在于,步骤1中支撑腿不打滑的不等式约束条件可以表示为:
其中,μ为刀锋腿与崎岖壁面之间的摩擦系数,和/>分别为第i个支撑腿沿崎岖壁面前向、垂向和横向的力,其中,i=1,2,3,4,5,6。
3.根据权利要求1所述的稳定空间生成和稳定性增强方法,其特征在于,步骤2包括:
设定基于驱动力和驱动力矩/>的六维探索状态量,同时,设置N个探索方向其中,l=1,…,N;
基于上述动力学模型和支撑腿不打滑的不等式约束条件,利用最优化方法最大化其中,/>获取在当前运动状态下第l个方向上的最大许用驱动力和力矩;
通过解算N个探索方向上的最大许用驱动力和力矩,便可获取当前运动状态下水下机器人在崎岖面攀爬的稳定空间。
4.根据权利要求1所述的稳定空间生成和稳定性增强方法,其特征在于,步骤3包括:
通过收集水下机器人在不同运动状态下的稳定空间数据和状态数据,来训练运动状态量与稳定空间之间的映射关系,映射关系采用BP神经网络来表征;当映射关系训练完成时,输入水下机器人的运动状态,从而生成水下机器人的稳定空间。
5.根据权利要求1所述的稳定空间生成和稳定性增强方法,其特征在于,步骤4包括:
水下机器人第j,j=1,2,…,M组运动状态下的第l个方向上的最大许用驱动力和力矩矢量可定义为:
其中,M为水下机器人运动状态的总量,为第j组运动状态下的第l个方向上的驱动力和力矩,经过矢量/>末端且与矢量/>垂直的6维超平面方程可定义为:
其中,表示水下机器人六维状态空间的第k维状态量,k=1,…,6;
稳定空间内的任一驱动力和力矩到六维超平面的距离/>可表示为
其中,||·||为矢量的范数符号;
利用最优化方法寻找出第j组运动状态下Dj最大的最优驱动力F3和力矩M3,以使水下机器人的稳定裕度最大;
根据最优驱动力F3和力矩M3求解8个推进器的推力值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111224379.0A CN113885508B (zh) | 2021-10-20 | 2021-10-20 | 水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111224379.0A CN113885508B (zh) | 2021-10-20 | 2021-10-20 | 水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113885508A CN113885508A (zh) | 2022-01-04 |
CN113885508B true CN113885508B (zh) | 2023-07-21 |
Family
ID=79003886
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111224379.0A Active CN113885508B (zh) | 2021-10-20 | 2021-10-20 | 水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113885508B (zh) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004050966A (ja) * | 2002-07-19 | 2004-02-19 | Niigata Tlo:Kk | 脚型走行方式とその装置 |
CN102923207A (zh) * | 2012-11-13 | 2013-02-13 | 中国计量学院 | 可适应崎岖路面爬行的仿蜘蛛机器人的躯干装置 |
WO2014185791A1 (en) * | 2013-05-16 | 2014-11-20 | Høgskolen I Ålesund | Underwater vehicle |
CN112306076A (zh) * | 2020-11-06 | 2021-02-02 | 西北工业大学 | 一种混合驱动水下机器人动力学模型建立方法 |
CN113353217A (zh) * | 2021-07-23 | 2021-09-07 | 西北工业大学 | 一种水下机器人定向和稳定行走方法 |
-
2021
- 2021-10-20 CN CN202111224379.0A patent/CN113885508B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004050966A (ja) * | 2002-07-19 | 2004-02-19 | Niigata Tlo:Kk | 脚型走行方式とその装置 |
CN102923207A (zh) * | 2012-11-13 | 2013-02-13 | 中国计量学院 | 可适应崎岖路面爬行的仿蜘蛛机器人的躯干装置 |
WO2014185791A1 (en) * | 2013-05-16 | 2014-11-20 | Høgskolen I Ålesund | Underwater vehicle |
CN112306076A (zh) * | 2020-11-06 | 2021-02-02 | 西北工业大学 | 一种混合驱动水下机器人动力学模型建立方法 |
CN113353217A (zh) * | 2021-07-23 | 2021-09-07 | 西北工业大学 | 一种水下机器人定向和稳定行走方法 |
Non-Patent Citations (7)
Title |
---|
Adaptive sliding-mode attitude control for autonomous underwater vehicles with input nonlinearities;崔荣鑫;《Ocean Engineering》;全文 * |
Hydrodynamic design of an underwater hull cleaning robot and its evaluation;Man Hyung Lee;《International Journal of Naval Architecture and Ocean Engineering》;全文 * |
Neural Learning-Based Integrated Guidance and Control Algorithm of Multiple Underactuated AUVs;陈乐鹏;《2018 IEEE 8th International Conference on Underwater System Technology: Theory and Applications (USYS)》;全文 * |
Optimization of structural parameters of jet end in the underwater intelligent pigging robot;Shangjunnan Liu;《Ocean Engineering》;全文 * |
基于模型参考的多自主水下航行器自适应覆盖控制;严卫生;《系统工程与电子技术》;全文 * |
水下探伤机器人水动力系数与运动控制研究;褚洪贵;《中国优秀硕士论文 工程科技Ⅱ辑》;全文 * |
矢量推进自主水下航行器动力学建模及仿真;王玉;《天津大学学报(自然科学与工程技术版)》;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113885508A (zh) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103942383B (zh) | 一种深海作业型水下机器人的动力学和运动学估计方法 | |
CN104076689B (zh) | 一种全驱动式自主水下航行器协同控制方法 | |
Liu et al. | Investigation of the vectored thruster AUVs based on 3SPS-S parallel manipulator | |
CN107991872B (zh) | 欠驱动auv的虚拟锚泊水平面区域镇定控制方法及实现方法 | |
CN111857165B (zh) | 一种水下航行器的轨迹跟踪控制方法 | |
CN109540575B (zh) | 深海自主姿态调节取样系统及其姿态调节方法 | |
CN106840143B (zh) | 一种判别水下机器人姿态稳定的方法 | |
Bai et al. | Hydrodynamics of a flexible flipper for an underwater vehicle-manipulator system | |
CN112925332B (zh) | 无人艇与水下无人潜器联合系统的协同交会对接控制方法 | |
CN113885508B (zh) | 水下机器人崎岖面攀爬的稳定空间生成和稳定性增强方法 | |
Shi et al. | A fuzzy PID control method for the underwater spherical robot | |
CN108363400B (zh) | 一种基于虚拟锚泊的欠驱动auv三维控制区域镇定控制方法 | |
CN116736721B (zh) | 一种变结构水下机器人游走模式运动控制方法和系统 | |
CN107697248A (zh) | 一种深海作业型机器人推进器的误差及参数修正方法及方法 | |
Wang et al. | ABA*–Adaptive Bidirectional A* Algorithm for Aerial Robot Path Planning | |
CN108227723A (zh) | 一种水下机器人及其稳定性分析及结构优化的应用方法 | |
Kiselev et al. | Selected Features of Autonomous Underwater Robot Dynamics under Near-Bottom Equidistant Motion Control | |
Li et al. | Dynamics modeling and experiments of wave driven robot | |
CN113419545B (zh) | 一种具有滑移管制的水下机器人优化控制方法 | |
CN116679693A (zh) | 融合螺旋桨水动力的无人艇自抗扰控制方法 | |
CN113671977B (zh) | 一种海上作业船状态同步稳定鲁棒控制方法 | |
Tang et al. | A novel underwater snake-like robot with gliding gait | |
CN115755939A (zh) | 一种四旋翼水下航行器前进运动状态估计方法 | |
Guo et al. | Guidance and control of a biomimetic autonomous underwater vehicle using body-fin propulsion | |
Hou et al. | The impulse output control of multiple thrusters AUV based on the S-plane algorithm |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |