CN112327892B - 一种auv误差受限的抗扰控制方法 - Google Patents

一种auv误差受限的抗扰控制方法 Download PDF

Info

Publication number
CN112327892B
CN112327892B CN202011287662.3A CN202011287662A CN112327892B CN 112327892 B CN112327892 B CN 112327892B CN 202011287662 A CN202011287662 A CN 202011287662A CN 112327892 B CN112327892 B CN 112327892B
Authority
CN
China
Prior art keywords
auv
underwater
mechanical arm
interference
error
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
Application number
CN202011287662.3A
Other languages
English (en)
Other versions
CN112327892A (zh
Inventor
魏延辉
于恬
付敬恒
邵红
王玥玥
樊雅妮
郭彤旭
谢吉顺
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN202011287662.3A priority Critical patent/CN112327892B/zh
Publication of CN112327892A publication Critical patent/CN112327892A/zh
Application granted granted Critical
Publication of CN112327892B publication Critical patent/CN112327892B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05DSYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
    • G05D1/04Control of altitude or depth
    • G05D1/06Rate of change of altitude or depth
    • G05D1/0692Rate of change of altitude or depth specially adapted for under-water vehicles

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明涉及一种AUV误差受限的抗扰控制方法,基于AUV在执行水下定点任务所存在的问题,如水下机械臂耦合干扰与AUV系统自身的结构性不确定性和海流干扰等因素所引起的非结构性不确定性的影响,设计了一种AUV误差受限的抗扰控制方法,可以在补偿水下机械臂耦合干扰,抵抗AUV结构和非结构性不确定性的同时,使得AUV在三维空间的轨迹和姿态跟踪误差可以按照设定的性能在有限时间内收敛于期望轨迹,满足水下定点作业的要求。本发明将AUV和水下机械臂的控制器进行有机整合,整体控制性能好,为AUV水下定点作业的抗扰动提供重要帮助,具有重要工程价值和实际意义。

Description

一种AUV误差受限的抗扰控制方法
技术领域
本发明涉及一种AUV误差受限的抗扰控制方法,属于水下AUV运载水下机械臂定点作业技术领域。
背景技术
AUV发展的关键技术包括提高续航能力技术、水下高速通信技术、精确导航定位技术、高级自主决策技术和运动控制技术等。在诸多关键技术中,运动控制技术尤为重要,只有具备了精准度高,鲁棒性强的运动控制方法AUV才能在远洋、深海、恶劣的海况环境下自主地完成任务。水下自主机器人难以控制的主要原因有:机器人高度的非线性、时变的动态特性,水动力系数的不确定性,各自由度运动间的耦合性,存在未知水流的扰动等。
欠驱动AUV无法进行完成悬停状态下的水下定点作业,因此,在对UVMS的运动过程进行规划时,不仅要考虑AUV携带机械臂从当前运动到目标定点处进行作业,而且,在AUV在离目标物的距离尚有余量时,水下机械臂的末端执行器的必须运动到UVMS任务规划层所指定的位置且保持不变,剩余距离由AUV携带机械臂运动到目标定点处完成作业任务。在整个水下定点作业的运动过程当中,UVMS因为AUV本体和水下机械臂存在着动力学耦合,水下机械臂的运动状态会对AUV的位置和姿态产生影响,影响系统整体的作业精度,UVMS运动控制系统的精确性、快速性和鲁棒性是其是否能够完成水下定点作业的关键。
发明内容
针对上述现有技术,本发明要解决的技术问题是提供一种AUV误差受限的抗扰控制方法,可以针对海流干扰,海底环境复杂,AUV自身结构不确定性等问题对AUV的姿态与位置进行控制,补偿水下机械臂的耦合干扰,抵抗AUV结构和非结构性不确定性,使得AUV在三维空间的姿态与轨迹跟踪误差可以按照设定的性能在有限时间内收敛于期望轨迹。
为解决上述技术问题,本发明的一种AUV误差受限的抗扰控制方法,包括以下步骤:
步骤1:对AUV的模型进行简化得到AUV五自由度的运动学与动力学模型;
步骤2:建立基于视线法和设定性能函数的三维空间误差动力学模型,再通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰,并通过直接补偿的方法抵消水下机械臂对于AUV本体的影响;
步骤3:根据AUV的误差动力学模型设计AUV抗扰控制器,通过神经网络控制补偿AUV的结构性的不确定性,通过自适应控制补偿AUV的非结构性的不确定性和水下机械臂的估算误差,通过误差增益控制项减少三维空间轨迹和姿态跟踪过程中的误差;
步骤4:完成水下定点作业过程中的AUV的运动控制任务。
本发明还包括:
1.步骤1中AUV五自由度的运动学模型具体为:
Figure GDA0004015391610000021
其中,[x y z]T和[φ θ ψ]T是AUV在惯性坐标系下的位置和姿态信息,[u v w]T是AUV载体坐标系原点线速度在xb、yb和zb轴上投影,[p q r]T是xb、yb和zb轴上是AUV载体坐标系原点的角速度在xb、yb和zb轴上投影;
AUV的动力学模型具体为:
Figure GDA0004015391610000022
Figure GDA0004015391610000023
其中,MD与MP分别为有驱动和无驱动条件下AUV本体的惯性矩阵,包括刚体质量矩阵和附加质量矩阵;DD(υ)与DP(υ)分别为有驱动和无驱动条件下AUV本体的阻尼阵;CD(υ)与CP(υ)分别为有驱动和无驱动条件下AUV本体的哥氏力阵和向心力阵;τ是AUV由推进器推力、舵角输入共同构成的向量形式,GD(θ)是AUV所受重力和浮力矩阵组成的恢复力矩阵,FD(vD)与FP(vP)分别为有驱动和无驱动条件下AUV动力学方程的未建模动态和高阶非线性项,τw1和τw2为海流引起的有界干扰信号,τm1与τm2为机械臂对AUV本体的干扰,vD=[u q r]T为驱动方向上的速度向量,vP=[v w]为无驱动方向上的速度向量。
2.步骤2中通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰具体为:
Figure GDA0004015391610000024
其中,
Figure GDA0004015391610000025
3.步骤3中AUV抗扰控制器具体为:
Figure GDA0004015391610000031
其中,
Figure GDA0004015391610000032
和/>
Figure GDA0004015391610000033
为正定对角增益矩阵,τnc是神经网络控制项,用于补偿转换后的AUV误差动力学方程的结构性不确定项,τac是自适应控制项,用于补偿转换后的AUV误差动力学方程中的由海流干扰引起的不确定性项,/>
Figure GDA0004015391610000034
为机械臂耦合干扰的前馈补偿项。
本发明的有益效果:本发明是是一种当AUV控制器在有限时间内跟踪UVMS任务规划层给出的期望轨迹时的一种抗干扰避免碰撞的控制系统。针对在水下机械臂耦合干扰影响下的AUV三维空间的运动控制问题,提出了一种AUV误差受限抗扰控制方法,对AUV的三维空间中的位置和姿态进行控制。基于视线法和设定性能函数得到AUV的误差动力学方程,基于此方程设计AUV抗扰控制器,并通过Newton-Eular方程法估算法机械臂对于AUV本体的耦合干扰并进行实时补偿,采用神经网络控制补偿AUV系统的结构性不确定性,自适应控制补偿神经网络估计误差,非结构性不确定性和机械臂耦合干扰的补偿误差。具有较好的抑制能力,可以确保水下定点作业过程中AUV的位置和姿态控制,达到良好的作业效果。
(1)本发明提供了一种误差受限抗扰控制方法,基于视线法和设定函数建立了三维空间的AUV误差动力学方程,并基于此方程设计了AUV抗扰控制器。该控制器采用牛顿-欧拉法实现了对水下机械臂对于AUV本体的耦合干扰的实时估计与补偿。
(2)本发明通过神经网络控制方法减少AUV结构性不确定对其运动带来的影响,通过自适应控制方法抵消AUV非结构性不确定性和水下机械臂的补偿误差的影响,具有较好的抗扰能力和鲁棒性能。并通过水下机械臂静止状态下的仿真实验验证了该控制器可以有效实现欠驱动AUV在设定性能下的三维空间中的轨迹跟踪,并且对水下机械臂耦合干扰,结构性不确定和海流干扰等其他非结构性的不确定性具有较强的抑制能力。
(3)本发明提出的抗扰控制方法,可以在机械臂耦合干扰、外加海流干扰,未建模动态存在的情况下,设定性能参数,调整控制参数,不仅可以使其稳定跟踪误差,还可以满足其跟踪误差的设定性能。其稳定性得到验证,控制效果良好,解决了AUV本体和水下机械臂因为存在动力学耦合,导致水下机械臂的运动状态会对AUV的位置和姿态产生影响,从而影响系统整体的作业效率与作业效果的问题。
(4)本发明将AUV与水下机械臂有机融合,整个水下定点作业过程中,AUV的三维空间轨迹仍按设定性能以较高精度趋近于期望轨迹,AUV抗扰控制器可以较好抑制水下机械臂的耦合干扰、AUV的结构性不确定性与非结构性不确定性的影响,具有较强的鲁棒性。
因为与AUV的运动学耦合,水下机械臂末端执行器的跟踪误差有所增大,但跟踪精度仍可以达到2.5cm左右,可以满足某些水下定点作业的需求,例如,海底热液和冷泉处的生物取样。
附图说明
图1是运动控制流程图;
图2是UVMS系统的整体示意图;
图3是AUV三维空间轨迹跟踪框图;
图4是AUV抗扰动控制器框图;
图5是水下机械臂任务空间的控制框图;
图6是水下机械臂关节空间的控制框图;
图7是AUV在三维空间的距离跟踪误差;
图8是AUV在三维空间的姿态跟踪误差;
图9是水下机械臂的任务空间轨迹。
具体实施方式
下面结合附图对本发明具体实施方式做进一步说明。
本发明针对AUV在执行水下定点任务所存在的问题,如水下机械臂耦合干扰与AUV系统自身的结构性不确定性和海流干扰等因素所引起的非结构性不确定性的影响,设计了一种AUV误差受限的抗扰控制方法,可以在补偿水下机械臂耦合干扰,抵抗AUV结构和非结构性不确定性的同时,使得AUV在三维空间的轨迹和姿态跟踪误差可以按照设定的性能在有限时间内收敛于期望轨迹,满足水下定点作业的要求。本发明将AUV和水下机械臂的控制器进行有机整合,整体控制性能好,为AUV水下定点作业的抗扰动提供重要帮助,具有重要工程价值和实际意义。
考虑到欠驱动AUV无法实现悬停状态下的定点作业、减少整个系统能量消耗和提高系统整体的控制效率,UVMS系统的任务规划层应分三个阶段对水下机械臂运动轨迹进行规划:AUV从当前位置运动到定点位置的初期阶段,水下机械臂的各个关节保持静止状态;在接近水下定点目标的阶段,水下机械臂按照UVMS任务规划层规划好的轨迹进行展开,在距离目标定点还有一段既定的缓冲距离,水下机械臂末端执行器到达既定的位置并保持不变;最后阶段AUV按照任务规划层给出的期望轨迹,调节自身的位置和姿态,携带机械臂运动到定点位置协助其末端完成作业任务。
水下机械臂系统的控制目标为在有限时间内末端执行器跟踪上期望轨迹且对水下环境具有较强的鲁棒性。水下机械臂的任务空间控制主要针对水下定点作业过程中其末端执行器的闭环控制任务,要求其末端可以在设定时间内快速到达期望位置。基于其水下作业的任务需求,采用一种自适应迭代滑模控制方法来实现水下机械臂任务空间轨迹的跟踪控制,通过引入初始状态的修正因子来克服初始重置误差带来的影响
本发明技术方案如下:
步骤1:通过对AUV的模型的简化,得到五自由度的AUV运动学与动力学模型。
步骤2:建立基于视线法和设定性能的三维空间的误差动力学模型,再通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰。
步骤3:基于步骤2设计AUV抗扰控制器,通过直接补偿的方法抵消水下机械臂对于AUV本体的影响,通过神经网络控制补偿AUV的结构性的不确定性,通过自适应控制补偿AUV的非结构性的不确定性和水下机械臂的估算误差,通过误差增益控制项减少三维空间轨迹和姿态跟踪过程中的误差,完成水下定点作业过程中的AUV的运动控制任务。
步骤4:采用一种自适应迭代滑模控制方法来实现水下机械臂任务空间轨迹的跟踪控制,再结合AUV运动学模型、动力学模型以及AUV抗扰控制算法可以建立控制系统模型,并基于该模型进行仿真,验证了该控制方法具有较好的抑制能力,可以确保水下定点作业过程中AUV的位置和姿态控制。
结合图1,本发明的方法包括如下步骤:
步骤1:通过对AUV的模型的简化,得到五自由度的AUV运动学与动力学模型。
步骤2:建立基于视线法和设定性能的三维空间的误差动力学模型,再通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰。
步骤3:基于步骤2设计AUV抗扰控制器,通过直接补偿的方法抵消水下机械臂对于AUV本体的影响,通过神经网络控制补偿AUV的结构性的不确定性,通过自适应控制补偿AUV的非结构性的不确定性和水下机械臂的估算误差,通过误差增益控制项减少三维空间轨迹和姿态跟踪过程中的误差,完成水下定点作业过程中的AUV的运动控制任务。
步骤4:采用一种自适应迭代滑模控制方法来实现水下机械臂任务空间轨迹的跟踪控制,再结合AUV运动学模型、动力学模型以及AUV抗扰控制算法可以建立控制系统模型,并基于该模型进行仿真,验证了该控制方法具有较好的抑制能力,可以确保水下定点作业过程中AUV的位置和姿态控制。
实施2:如附图2所示,UVMS系统主要由AUV和水下机械臂两个部分组成:
AUV本体:AUV整体呈鱼雷形,可以有效减少其水中运动的阻力,采用推进器和十字舵的共同作用去控制AUV的位置和姿态,因为系统控制变量的个数为3,但是其运动的自由度为6,所以该AUV为欠驱动系统。
水下机械臂:UVMS系统选用带有末端执行器的四自由度机械臂作为作业工具,水下机械臂由腰部和腕部的回转关节、大臂和小臂的摆动关节四个部分组成,通过四个关节的组合旋转运动,可以完成水下机械臂末端执行器的位置与姿态的控制。
由AUV本体与水下机械臂所组成的UVMS系统是本发明所针对的对象。
下面通过对AUV的模型的简化,得到步骤1中五自由度的AUV运动学与动力学模型。
AUV相关参数说明如下:
Figure GDA0004015391610000061
其中,[X Y Z]T为AUV在xb、yb和zb轴上所受到的力,[K M N]T是AUV在各轴上所受到的力矩,[u v w]T是AUV载体坐标系原点线速度在xb、yb和zb轴上投影,为[p q r]T是xb、yb和zb轴上是AUV载体坐标系原点的角速度在xb、yb和zb轴上投影。[x y z]T和[φ θ ψ]T是AUV在惯性坐标系下的位置和姿态信息。
忽略横摇对AUV运动的影响,对AUV的数学模型进行简化,AUV五自由度的运动学可以写为:
Figure GDA0004015391610000062
根据有驱动和无驱动的方向,当欠驱动AUV的位姿信息可以传感器实时测得并进行反馈时,可以AUV的动力学方程分为两个部分:
Figure GDA0004015391610000063
Figure GDA0004015391610000064
其中,MD与MP分别为有驱动和无驱动条件下,AUV本体的惯性矩阵,包括刚体质量矩阵和附加质量矩阵;同理DD(υ)与DP(υ)分别为有驱动和无驱动条件下,AUV本体的阻尼阵;CD(υ)与CP(υ)分别为有驱动和无驱动条件下,AUV本体的哥氏力阵和向心力阵;τ是AUV由推进器推力、舵角输入共同构成的向量形式。GD(θ)是AUV所受重力和浮力矩阵组成的恢复力矩阵,FD(vD)与FP(vP)分别为有驱动和无驱动条件下,AUV动力学方程的未建模动态和高阶非线性项,τw1和τw2为海流引起的有界干扰信号,τm1与τm2为机械臂对AUV本体的干扰,vD=[u qr]T为驱动方向上的速度向量,vP=[v w]为无驱动方向上的速度向量。
实施3:跟据步骤2,建立基于视线法和设定性能的三维空间的误差动力学模型,再通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰。
AUV相关参数说明如下:
Figure GDA0004015391610000071
如附图3所示,AUV跟踪的位置误差(xe、ye和ze)与AUV基于虚拟向导的跟踪误差(ρe、γ和φ)存在如下的几何关系:
xe=ρecosφcosγ
ye=ρecosφsinγ
ze=ρesinφ
为了方便设计,实际AUV和虚拟向导在惯性坐标下的距离和姿态误差可以转换到AUV的载体坐标系上:
Figure GDA0004015391610000072
对其进行求导:
Figure GDA0004015391610000073
Figure GDA0004015391610000074
Figure GDA0004015391610000075
为了方便表示,定义α1,α2和α3
Figure GDA0004015391610000081
Figure GDA0004015391610000082
Figure GDA0004015391610000083
在AUV的动力学模型存在机械臂耦合干扰,AUV结构性不确定性和海流干扰等非结构性不确定性的情况下,使得AUV的三维空间的轨迹跟踪误差收敛到原点,且满足预先设定动态性能。根据视线法,定义基于虚拟向导的跟踪误差如下式所示:
Figure GDA0004015391610000084
Figure GDA0004015391610000085
/>
Figure GDA0004015391610000086
对式(10)、式(11)和式(12)进行求导,将式(4)、式(5)、式(6)、式(7)、式(8)和式(9)代入可得到跟踪误差如下:
Figure GDA0004015391610000087
Figure GDA0004015391610000088
Figure GDA0004015391610000089
为使跟踪误差按照设定的动态和静态性能收敛,对跟踪误差进行设定:
a1ξ1(t)<ρe(t)<b1ξ2(t) (16)
-a2ξ2(t)<φ(t)<b2ξ2(t) (17)
-a3ξ3(t)<γ(t)<b3ξ3(t) (18)其中,b1>a1,a2,b2,a3和b3为正常数,ξi(t)(i=1,2,3)函数设计为如下形式:
ξi(t)=(ξi0i∞)e-λiti∞ (19)
函数ξi(t)为单调递减函数,ξi0为ξi(t)的初始状态,ξi∞为任意小的常数且为ξi(t)的极限,λi与跟踪误差的逼近速率相关,ξi0>ξi∞>0,ξi∞>0。
为了实现控制目和确保跟踪误差的设定的动态和静态性能,设计三个严格递增函数
Figure GDA00040153916100000912
y2φ)和y3γ),/>
Figure GDA00040153916100000913
εφγ∈L,且满足如下性质:
Figure GDA00040153916100000914
Figure GDA0004015391610000091
Figure GDA0004015391610000092
y1(0)=a11,y2(0)=0,y3(0)=0 (23)
因此,函数
Figure GDA00040153916100000915
y2φ)和y3γ)应为:
Figure GDA0004015391610000093
Figure GDA0004015391610000094
Figure GDA0004015391610000095
其中,υ1=0.5ln(ε1/(b1-a11)),υ2=0.5ln(a1/b1)和υ3=0.5ln(a3/b3),且pi>0,i=1,2,3为常数。
Figure GDA00040153916100000916
φ=ξ2y2φ)和γ=ξ3y3γ),因为ξ1(t),ξ2(t),ξ3(t),/>
Figure GDA00040153916100000917
y2φ)和y3γ)为单调递增函数,且ξ123≠0,所以,
Figure GDA0004015391610000096
Figure GDA0004015391610000097
存在,具体形式如下:
Figure GDA0004015391610000098
Figure GDA0004015391610000099
Figure GDA00040153916100000910
为了方便表示,令
Figure GDA00040153916100000911
Figure GDA00040153916100000918
εφ和εγ是ρe,φ和γ进行变换后的误差。对上式再进行求导,得到变换误差的动态方程。
Figure GDA0004015391610000101
定义新的误差向量
Figure GDA0004015391610000109
将式(13)、式(14)和式(15)代入上式可得:
Figure GDA0004015391610000102
将上式写成紧凑的形式:
Figure GDA0004015391610000103
由式(32)可以得到vD及其
Figure GDA0004015391610000104
Figure GDA0004015391610000105
将上式代入到AUV动力学模型中,两边乘以R-Te),得到转换后的误差动力学方程,即建立基于视线法和设定性能的三维空间的误差动力学模型:
Figure GDA0004015391610000106
/>
其中,
Figure GDA0004015391610000107
表示由海流干扰引起的不确定性项等/>
Figure GDA0004015391610000108
机械臂耦合干扰。
水下机械臂在进行运动时,会对AUV本体的姿态产生干扰,使得AUV偏离期望的轨迹,为了抵消水下机械臂对AUV本体位置和姿态的干扰,需要对水下机械臂对于AUV本体进行估计并进行实时补偿。因为本发明的AUV为欠驱动水下机器人,在横摇方向上没有驱动,如果水下机械臂的腰部回转关节在运动时,会对AUV在横摇方向产生干扰,使AUV失去控制,所以,为了保证AUV系统的可控性,保持AUV腰部关节为静止状态,大臂和小臂关节只可以在AUV的纵倾面上运动,此时的水下机械臂可以退化成纵倾面上的二自由度机械臂。
牛顿-欧拉公式是建立在机械臂连杆之间的所有力平衡关系的基础之上,得出的方程组,其结构使得解具有递归形式,前向递归用于连杆速度和加速度传递,后向递归用于力的传递。
对于转动关节存在如下关系:
(1)前向递归
Figure GDA0004015391610000111
Figure GDA0004015391610000112
Figure GDA0004015391610000113
Figure GDA0004015391610000114
其中,
Figure GDA0004015391610000115
和/>
Figure GDA0004015391610000116
是坐标系i的角速度和角加速度,/>
Figure GDA0004015391610000117
和/>
Figure GDA0004015391610000118
是坐标系i原点的线速度和线加速度,/>
Figure GDA0004015391610000119
是连杆i的质心的线速度,/>
Figure GDA00040153916100001110
和/>
Figure GDA00040153916100001111
分别表示坐标系i的z轴与在坐标系i下的坐标系i-1的原点到坐标系i的距离,i-1Ri∈R3×3i-1Pi∈R3×3分别代表{i}系在{i-1}系下的姿态矩阵和位置向量,左上角的角标i表示以第i连杆为基准建立的坐标系,下文同理。
得到合力与力矩:
Figure GDA00040153916100001112
Figure GDA00040153916100001113
(2)后向递归
Figure GDA00040153916100001114
Figure GDA00040153916100001115
Figure GDA00040153916100001116
对于移动关节,有如下关系
(1)前向递归
Figure GDA00040153916100001117
Figure GDA00040153916100001118
Figure GDA00040153916100001119
/>
Figure GDA00040153916100001120
得到合力与力矩:
Figure GDA0004015391610000121
Figure GDA0004015391610000122
(2)后向递归
Figure GDA0004015391610000123
Figure GDA0004015391610000124
Figure GDA0004015391610000125
下面就水下机械臂对AUV本体的耦合扰动作具体的分析,为了方便起见,假设水下机械臂连杆1、2的质量分别是m1和m2,质心位于连杆中心位置,令l1,l2为两连杆质心到关节轴之间的距离,L1和L2为两连杆的长度,连杆1、2的惯性可以表示为
Figure GDA0004015391610000126
和/>
Figure GDA0004015391610000127
连杆1、2的重心及其末端的可以表示向量形式/>
Figure GDA0004015391610000128
和/>
Figure GDA0004015391610000129
通过编码器测得机械臂两个关节的转角为/>
Figure GDA00040153916100001210
和/>
Figure GDA00040153916100001211
向前进行递推,对于连杆1的速度和加速度信息如下:
Figure GDA00040153916100001212
Figure GDA00040153916100001213
Figure GDA00040153916100001214
Figure GDA00040153916100001215
Figure GDA0004015391610000131
/>
Figure GDA0004015391610000132
对于连杆2的速度和加速度信息如下:
Figure GDA0004015391610000133
Figure GDA0004015391610000134
Figure GDA0004015391610000135
Figure GDA0004015391610000136
Figure GDA0004015391610000137
Figure GDA0004015391610000138
向后进行递推,对于连杆2的力和力矩信息如下:
Figure GDA0004015391610000141
Figure GDA0004015391610000142
因此可以得到对于连杆1的力和力矩信息如下:
Figure GDA0004015391610000143
将力和力矩1f和1μ1转换到坐标系0上,有:
Figure GDA0004015391610000144
其中,s1,s2,s3和s4分别表示sinq1,sinq2,sinq3和sinq4,c1,c2,c3和c4分别表示cosq1,cosq2,cosq3和cosq4
综上,可以得到水下机械臂末端不夹持物体和夹持物体对机械臂基座的作用力,为了方便分析水下机械臂对AUV本体的干扰和简化计算,将水下机械臂的基座坐标系放在AUV载体坐标系的正下方,水下机械臂的基座坐标系的原点在AUV载体坐标系下表示为0rB=[0pBx 0pBy 0pBz]T,其中0pBx=0和0pBy=0,则水下机械臂对AUV本体的耦合干扰为:
Figure GDA0004015391610000145
Figure GDA0004015391610000146
因此,通过牛顿-欧拉法估算出的水下机械臂对AUV本体的耦合干扰的作用为
Figure GDA0004015391610000147
实施4:AUV抗扰控制器的框图如附图4示。根据AUV在水下定点作业过程中的任务需求,基于实施3由视线法和设定函数建立的三维空间的AUV误差动力学方程以及牛顿-欧拉(Newton-Eular)公式推导得出的水下机械臂对于AUV本体的扰动力矩,在AUV抗扰控制器中直接进行补偿。AUV抗扰控制器由误差控制项,神经网络控制项,自适应控制项和机械臂耦合干扰估计项组成。误差控制项来减少AUV运动过程中的跟踪误差,机械臂耦合干扰估计项实时补偿机械臂对于AUV本体的干扰,神经网络控制项补偿AUV系统的结构性的不确定项,自适应控制项补偿神经网络估计误差、机械臂耦合干扰补偿误差和外部海流干扰引起的非结构性误差。
为了抵消水下机械臂对AUV本体位姿的影响,采用前馈补偿的方式抵消水下机械臂对AUV的干扰,令
Figure GDA0004015391610000151
在AUV抗扰控制器中直接进行补偿,因为本发明并没有考虑水动力,传感器测量误差及信息传输延迟和欠驱动补偿的不对称性等因素,机械臂耦合干扰估计存在误差,并不能完全消除水下机械臂对AUV的影响,估计误差需要鲁棒控制算法在系统闭环中消除。
基于转换后的AUV误差动力学方程设计抗扰控制器,将实施三得到的式(34)转换为:
Figure GDA0004015391610000152
定义新的误差变量
Figure GDA0004015391610000153
其中,
Figure GDA00040153916100001511
为正对角增益矩阵。
将式(72)代入式(71),可得一阶的AUV误差动力学方程:
Figure GDA0004015391610000154
其中,
Figure GDA0004015391610000155
Θ有界,即|Θ||≤c5||xe||+c6||xe||2,/>
Figure GDA0004015391610000156
c5和c6为未知正常数。
则得到AUV的抗扰控制器设计为:
Figure GDA0004015391610000157
其中,
Figure GDA0004015391610000158
和/>
Figure GDA0004015391610000159
为正定对角增益矩阵,τnc是神经网络控制项,用于补偿转换后的AUV误差动力学方程的结构性不确定项,τac是自适应控制项,用于补偿转换后的AUV误差动力学方程中的由海流干扰引起的不确定性项,/>
Figure GDA00040153916100001510
为机械臂耦合干扰的前馈补偿项。
该控制器可以有效实现欠驱动AUV在设定性能下的三维空间中的轨迹跟踪,并且对水下机械臂耦合干扰,结构性不确定和海流干扰等其他非结构性的不确定性具有较强的抑制能力。通过有机整合AUV和水下机械臂的控制算法实现了水下定点过程中的AUV可以在设定的时间内完成携带水下机械臂从当前位置到目标位置的运动控制任务,水下机械臂的末端可以快速跟踪上期望轨迹并在仿真设定的时间到达期望位置,且末端的作业精度可以定点作业任务的需求,同时也确保了分离式UVMS运动控制方案在水下定点作业任务中的可行性和有效性。
实施5:如附图5、附图6所示,因水下机械臂通常工作在复杂的水下的环境,机械臂关节存在机械摩擦、外加海流干扰等多种不确定因素,所以传统的切换控制方式,因为不连续的控制会造成严重的滑模抖振。使用边界层的方式来进行抗扰,将切换控制连续化,可以有效降低抖振,但是系统的鲁棒性会受到影响,同时系统由有限时间收敛变为渐进收敛。为了解决这个问题,采用一种自适应边界层的控制方式,根据外界环境的变化,控制器自发调节边界层的参数,可以有效提高水下机械臂对于复杂环境的适应能力和作业的精度。
水下机械臂关节空间控制主要应用在UVMS运动过程中遭遇高危地形(如礁群和峡湾)时的水下机械臂关节的快速调整问题,因此,要求所设计的关节空间控制器可以在有限时间内实现快速收敛,且在复杂的水下环境下具有较强的鲁棒性。采用新型的非奇异快速终端滑模自适应控制方法,可以确保其关节在有限时间内收敛到期望位置,并且通过引入非线性分段函数来设计非奇异快速终端滑模面,来改善其远离滑模面的收敛过慢问题,通过自适应调整边界层,进一步增强在复杂海况下的水下机械臂关节空间控制的鲁棒性和关节空间控制精度的同时,降低滑模抖振。
有效改善了UVMS系统在进行水下定点作业时,AUV和水下机械臂在运动状态下,水下环境无法精确建模和外界扰动给定点作业带来的不良影响;为UVMS系统找到了一种在水下定点作业状况下的控制其姿态和位置的有效方法;确保了UVMS运动控制系统的精确性、快速性和鲁棒性。
如附图7、附图8、附图9所示,图7中位于最上方的虚线曲线是b1ξ1,位于最下方的虚线曲线是a1ξ1,图8中,上图中位于最上方的虚线曲线是b2ξ2,位于最下方的虚线曲线是-a2ξ2,下图中位于最上方的虚线曲线是b3ξ3,位于最下方的虚线曲线是-a3ξ3。跟据仿真结果,本发明整合的UVMS控制器,可以较好满足水下定点作业过程中的AUV和水下机械臂的任务需求,AUV可以在运动过程中有效抵抗水下机械臂的耦合干扰和其他的不确定性,并且以设定的性能跟踪上期望的三维空间轨迹,水下机械臂可以在有限时间内快速跟踪上任务空间的轨迹,UVMS系统的作业精度可以满足水下定点作业的需求。

Claims (1)

1.一种AUV误差受限的抗扰控制方法,其特征在于,包括以下步骤:
步骤1:对AUV的模型进行简化得到AUV五自由度的运动学与动力学模型;
步骤2:建立基于视线法和设定性能函数的三维空间误差动力学模型,再通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰,并通过直接补偿的方法抵消水下机械臂对于AUV本体的影响;
步骤3:根据AUV的误差动力学模型设计AUV抗扰控制器,通过神经网络控制补偿AUV的结构性的不确定性,通过自适应控制补偿AUV的非结构性的不确定性和水下机械臂的估算误差,通过误差增益控制项减少三维空间轨迹和姿态跟踪过程中的误差;
步骤4:完成水下定点作业过程中的AUV的运动控制任务;
步骤1所述AUV五自由度的运动学模型具体为:
Figure FDA0004015391600000011
Figure FDA0004015391600000012
Figure FDA0004015391600000013
Figure FDA0004015391600000014
Figure FDA0004015391600000015
其中,[x y z]T和[φ θ ψ]T是AUV在惯性坐标系下的位置和姿态信息,[u v w]T是AUV载体坐标系原点线速度在xb、yb和zb轴上投影,[p q r]T是xb、yb和zb轴上是AUV载体坐标系原点的角速度在xb、yb和zb轴上投影;
AUV的动力学模型具体为:
Figure FDA0004015391600000016
Figure FDA0004015391600000017
其中,MD与MP分别为有驱动和无驱动条件下AUV本体的惯性矩阵,包括刚体质量矩阵和附加质量矩阵;DD(vD)与DP(vD)分别为有驱动和无驱动条件下AUV本体的阻尼阵;CD(vD)与CP(vD)分别为有驱动和无驱动条件下AUV本体的哥氏力阵和向心力阵;τ是AUV由推进器推力、舵角输入共同构成的向量形式,GD(θ)是AUV所受重力和浮力矩阵组成的恢复力矩阵,FD(vD)与FP(vP)分别为有驱动和无驱动条件下AUV动力学方程的未建模动态和高阶非线性项,τw1(t)和τw2(t)为有驱动方向与无驱动方向下海流引起的有界干扰信号,τm1与τm2为有驱动与无驱动条件下机械臂对AUV本体的干扰,vD=[u q r]T为驱动方向上的速度向量,vP=[v w]为无驱动方向上的速度向量;
步骤2所述通过牛顿-欧拉法估算出水下机械臂对于AUV本体的耦合干扰具体为:
Figure FDA0004015391600000021
其中,
Figure FDA0004015391600000022
步骤3所述AUV抗扰控制器具体为:
Figure FDA0004015391600000023
其中,
Figure FDA0004015391600000024
Figure FDA0004015391600000025
为正定对角增益矩阵,τnc是神经网络控制项,用于补偿转换后的AUV误差动力学方程的结构性不确定项,τac是自适应控制项,用于补偿转换后的AUV误差动力学方程中的由海流干扰引起的不确定性项,
Figure FDA0004015391600000026
为机械臂耦合干扰的前馈补偿项,
Figure FDA0004015391600000027
CN202011287662.3A 2020-11-17 2020-11-17 一种auv误差受限的抗扰控制方法 Active CN112327892B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011287662.3A CN112327892B (zh) 2020-11-17 2020-11-17 一种auv误差受限的抗扰控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011287662.3A CN112327892B (zh) 2020-11-17 2020-11-17 一种auv误差受限的抗扰控制方法

Publications (2)

Publication Number Publication Date
CN112327892A CN112327892A (zh) 2021-02-05
CN112327892B true CN112327892B (zh) 2023-03-24

Family

ID=74322401

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011287662.3A Active CN112327892B (zh) 2020-11-17 2020-11-17 一种auv误差受限的抗扰控制方法

Country Status (1)

Country Link
CN (1) CN112327892B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113093771B (zh) * 2021-04-07 2022-07-26 安徽大学 基于神经网络的水下机器人-机械手系统的建模方法及系统
CN113110527B (zh) * 2021-04-12 2023-12-01 大连海事大学 一种自主水下航行器有限时间路径跟踪的级联控制方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102768539A (zh) * 2012-06-26 2012-11-07 哈尔滨工程大学 基于迭代的自主水下航行器三维曲线路径跟踪控制方法
US8468108B1 (en) * 2009-06-30 2013-06-18 Zytek Communications Corporation Modeling efficiency over a range of velocities in underwater vehicles
KR101370649B1 (ko) * 2012-09-04 2014-03-10 주식회사 한화 무인잠수정의 경로 제어방법
CN104898688A (zh) * 2015-03-27 2015-09-09 哈尔滨工程大学 Uuv四自由度动力定位自适应抗扰滑模控制系统及控制方法
CN105843224A (zh) * 2016-03-25 2016-08-10 哈尔滨工程大学 基于神经动态模型反步法的auv水平面路径跟踪控制方法
CN107168312A (zh) * 2017-05-17 2017-09-15 哈尔滨工程大学 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法
CN111650948A (zh) * 2020-06-10 2020-09-11 哈尔滨工程大学 一种可底栖式auv的水平面轨迹快速跟踪控制方法
CN111880409A (zh) * 2020-08-06 2020-11-03 天津大学 一种基于数据驱动的自主水下航行器轨迹跟踪控制方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8468108B1 (en) * 2009-06-30 2013-06-18 Zytek Communications Corporation Modeling efficiency over a range of velocities in underwater vehicles
CN102768539A (zh) * 2012-06-26 2012-11-07 哈尔滨工程大学 基于迭代的自主水下航行器三维曲线路径跟踪控制方法
KR101370649B1 (ko) * 2012-09-04 2014-03-10 주식회사 한화 무인잠수정의 경로 제어방법
CN104898688A (zh) * 2015-03-27 2015-09-09 哈尔滨工程大学 Uuv四自由度动力定位自适应抗扰滑模控制系统及控制方法
CN105843224A (zh) * 2016-03-25 2016-08-10 哈尔滨工程大学 基于神经动态模型反步法的auv水平面路径跟踪控制方法
CN107168312A (zh) * 2017-05-17 2017-09-15 哈尔滨工程大学 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法
CN111650948A (zh) * 2020-06-10 2020-09-11 哈尔滨工程大学 一种可底栖式auv的水平面轨迹快速跟踪控制方法
CN111880409A (zh) * 2020-08-06 2020-11-03 天津大学 一种基于数据驱动的自主水下航行器轨迹跟踪控制方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
Modified Line-of-Sight Guidance Law With Adaptive Neural Network Control of Underactuated Marine Vehicles With State and Input Constraints;Raja Rout,等;《IEEE Transactions on Control Systems Technology》;20200610;第28卷(第5期);1902-1914 *
Research on Active Disturbance Rejection Control of Underwater Manipulator Based on Terminal Sliding Mode;Jialin He等;《2020 IEEE International Conference on Mechatronics and Automation (ICMA)》;20201026;856-861 *
基于串级变结构神经网络的机器人补偿控制;马毅潇等;《上海交通大学学报》;19990730(第07期);80-84 *
基于扰动观测器的机械臂自适应反演滑模控制;司彦娜等;《计算机测量与控制》;20180825(第08期);108-111+124 *
欠驱动AUV机动目标跟踪控制研究;焦文龙;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20200315(第03(2020)期);C036-287 *
欠驱动自主水下航行器空间曲线路径跟踪控制研究;苗建明等;《兵工学报》;20170915(第09期);125-135 *
水下机器人的轨迹跟踪控制方法研究;曾德伟;《中国优秀硕士学位论文全文数据库 信息科技辑》;20180715(第07(2018)期);I140-267 *

Also Published As

Publication number Publication date
CN112327892A (zh) 2021-02-05

Similar Documents

Publication Publication Date Title
CN109634307B (zh) 一种无人水下航行器复合航迹跟踪控制方法
CN107168312B (zh) 一种补偿uuv运动学和动力学干扰的空间轨迹跟踪控制方法
CN111650948B (zh) 一种可底栖式auv的水平面轨迹快速跟踪控制方法
Lee et al. A robust neural controller for underwater robot manipulators
CN112091976B (zh) 一种水下机械臂任务空间控制方法
CN112965371B (zh) 基于固定时间观测器的水面无人艇轨迹快速跟踪控制方法
CN112327892B (zh) 一种auv误差受限的抗扰控制方法
CN111857165B (zh) 一种水下航行器的轨迹跟踪控制方法
Han et al. Modeling and fuzzy decoupling control of an underwater vehicle-manipulator system
CN112068440B (zh) 基于模型预测控制的auv回收对接动力定位控制方法
Cui et al. A unified force control approach to autonomous underwater manipulation
CN111506095B (zh) 一种双刚体特征点间饱和固定时间相对位姿跟踪控制方法
Shojaei Three-dimensional tracking control of autonomous underwater vehicles with limited torque and without velocity sensors
Sarkar et al. Adaptive control of underwater vehicle-manipulator systems subject to joint limits
CN114967714A (zh) 一种自主式水下机器人抗扰运动控制方法及系统
CN108363400B (zh) 一种基于虚拟锚泊的欠驱动auv三维控制区域镇定控制方法
Cui et al. Impedance control of underwater vehicle-manipulator systems (UVMS)
CN113110512B (zh) 一种减弱未知干扰与抖振影响的可底栖式auv自适应轨迹跟踪控制方法
Santhakumar et al. Coupled, non-linear control system design for autonomous underwater vehicle (AUV)
CN114527772A (zh) 一种auv轨迹跟踪控制器设计方法与系统
An et al. Underwater motion control of a bio-inspired father-son robot based on hydrodynamics analysis
Morton et al. Generalized trajectory control for tree-structured aerial manipulators
Zhou et al. Multi-UUV formation coordination control based on virtual navigator
CN113305831B (zh) 空间机械臂自适应协调控制的位置观测器设计方法及系统
Sugiyama et al. A nonlinear disturbance observer using delayed estimates-its application to motion control of an underwater vehicle-manipulator system

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