CN113867352A - 一种全垫升气垫船路径跟踪方法 - Google Patents

一种全垫升气垫船路径跟踪方法 Download PDF

Info

Publication number
CN113867352A
CN113867352A CN202111151049.3A CN202111151049A CN113867352A CN 113867352 A CN113867352 A CN 113867352A CN 202111151049 A CN202111151049 A CN 202111151049A CN 113867352 A CN113867352 A CN 113867352A
Authority
CN
China
Prior art keywords
hovercraft
path
zero
track
angle
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.)
Granted
Application number
CN202111151049.3A
Other languages
English (en)
Other versions
CN113867352B (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 CN202111151049.3A priority Critical patent/CN113867352B/zh
Publication of CN113867352A publication Critical patent/CN113867352A/zh
Priority to GB2210207.3A priority patent/GB2618860A/en
Priority to PCT/CN2022/070350 priority patent/WO2023050636A1/zh
Application granted granted Critical
Publication of CN113867352B publication Critical patent/CN113867352B/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/02Control of position or course in two dimensions
    • G05D1/0206Control of position or course in two dimensions specially adapted to 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)

Abstract

本发明公开了一种全垫升气垫船路径跟踪方法,包括:建立全垫升气垫船航迹向控制数学模型和纵向速度控制数学模型;设计SFLOS的曲线路径导引算法,得到路径参数更新率
Figure DDA0003287118610000011
和期望航迹角ψwd;分别设计扩张状态观测器观测外界不确定性干扰;运用积分型障碍李亚普诺夫函数设计转艏力矩控制律对回转率进行约束,使得回转率在安全限界内,全垫升气垫船的实际航迹角跟踪期望航迹角ψwd;运用对称型对数型障碍李亚普诺夫函数设计纵向推力控制律对侧滑角进行约束,使得侧滑角在安全限界内,全垫升气垫船的实际纵向速度跟踪期望纵向速度。本发明实现了全垫升气垫船在跟踪上期望路径的前提下保证回转率和侧滑角分别在各自的安全限界以内。

Description

一种全垫升气垫船路径跟踪方法
技术领域
本发明涉及一种全垫升气垫船路径跟踪方法,特别是带侧滑角以及回转率约束的全垫升 气垫船路径跟踪方法,属于控制技术领域。
背景技术
全垫升气垫船是一种依靠垫升系统来垫态航行于水面上的特种船舶。全垫升气垫船在高 速航行过程中,特别是在回转过程中,如果操纵不当容易发生侧滑、甩尾乃至失稳翻船的安 全事故。因此,在对全垫升气垫船的运动控制进行研究时,非常有必要考虑到全垫升气垫船 的航行安全性,在对其路径跟踪进行研究中同样如此。
对于全垫升气垫船路径跟踪而言,影响其航行安全性的主要因素是回转率和侧滑角。 高速航行的全垫升气垫船在回转过程中的回转率容易超出其安全限界,进一步由过大的回转 率而引起全垫升气垫船在回转过程中的侧滑角超出其安全限界,因此需选用控制方法将回转 率和侧滑角分别限制在各自的安全限界以内。
发明内容
针对上述现有技术,本发明要解决的技术问题是提供一种基于SFLOS导引算法的带侧滑 角以及回转率约束的全垫升气垫船路径跟踪方法,既能跟踪直线路径又能跟踪曲线路径且通 过约束控制方法将回转率和侧滑角分别限制在各自的安全限界以内。
为解决上述技术问题,本发明的一种全垫升气垫船路径跟踪方法,包括以下步骤:
步骤1:建立全垫升气垫船四自由度运动数学模型,并在所建模型的基础上得到全垫升 气垫船航迹向控制数学模型和纵向速度控制数学模型;
步骤2:设计SFLOS的曲线路径导引算法,得到路径参数更新率
Figure BDA0003287118590000011
和期望航迹角ψwd
步骤3:依据步骤1所得的航迹向控制数学模型和纵向速度控制数学模型分别设计扩张 状态观测器观测外界干扰的不确定性;
步骤4:运用积分型障碍李亚普诺夫函数设计转艏力矩控制律对路径跟踪控制中的回转 率进行约束,使得路径跟踪过程中的回转率在安全限界内,全垫升气垫船的实际航迹角跟踪 期望航迹角ψwd
步骤5:运用对称型对数型障碍李亚普诺夫函数设计纵向推力控制律对路径跟踪控制中 的侧滑角进行约束,使得路径跟踪过程中的侧滑角在安全限界内,全垫升气垫船的实际纵向 速度跟踪期望纵向速度。
本发明还包括:
1.步骤1中全垫升气垫船航迹向控制数学模型和纵向速度控制数学模型具体为:
Figure BDA0003287118590000021
Figure BDA0003287118590000022
其中,ψw为实际航迹角,β为侧滑角且假设为可测的量,d3为外界不确定性干扰在转艏 自由度上的分量,Dr为转艏自由度上的动力学不确定性,d1为外界不确定性干扰在纵荡自由 度上的分量,Du为纵荡自由度上的动力学不确定性;r为船体坐标系下的全垫升气垫船的转 艏角速度/回转率;φ表示北东坐标系下的全垫升气垫船的横倾角;τR为空气舵产生的转艏力 矩;τP为空气螺旋桨产生的纵向推进力;v表示船体坐标系下的全垫升气垫船的横向速度;MzD表示艏摇自由度中除去τP和τR的合力,FxD表示纵荡自由度中除去τP和τR的合力。
2.步骤2中路径参数更新率
Figure BDA0003287118590000023
和期望航迹角ψwd满足:
Figure BDA0003287118590000024
其中,ψθ为北东坐标系向SF坐标系转换时的旋转角度,且ψθ=arctan2(y′θ,x′θ),xθ、yθ 为期望路径P(θ)上的任意路径点的坐标,满足:
Figure BDA0003287118590000025
其中,
Figure BDA0003287118590000026
kx为大于零 的设计参数;Δ为P(θ)处沿设定路径切向的大于零的前视距离参数,Δ的具体形式如下:
Figure BDA0003287118590000027
其中,μ、Δmax、Δmin均为大于零的设计参数,且Δmax>Δmin,可通过调节μ来调节Δ的 收敛速率。
3.步骤3中扩张状态观测器具体为:
针对航迹向控制数学模型,扩张状态观测器具体为:
Figure BDA0003287118590000031
其中,
Figure BDA0003287118590000032
为ψw的观测值,
Figure BDA0003287118590000033
为ψw的观测误差,
Figure BDA0003287118590000034
为r的观测值,
Figure BDA0003287118590000035
为Dr的观测值,λ11、 λ12、λ13均为大于零的设计参数,fal(e,α,δ)为具有如下形式的连续非线性幂次函数:
Figure BDA0003287118590000036
e即为
Figure BDA0003287118590000037
δ为大于零的设计参数,表示fal(e,α,δ)在原点附近线性段的长度,α为大 于零且小于1的设计参数,sgn(·)为符号函数;
针对纵向速度控制数学模型,设计扩张状态观测器具体为:
Figure BDA0003287118590000038
其中,
Figure BDA0003287118590000039
为u的观测值,eu为u的观测误差,
Figure BDA00032871185900000310
为Du的观测值,λ21、λ22均为大于零的设计参数,e即为eu
4.步骤4中转艏力矩控制律具体为:
Figure BDA00032871185900000311
其中,τR为空气舵产生的转艏力矩,MzD表示艏摇自由度中除去τR的合力,Iz为z轴的转动惯量,Dr为转艏自由度上的动力学不确定性,
Figure BDA00032871185900000312
为Dr的观测值,
Figure BDA00032871185900000313
K3均为大于零 的常数,re为回转率偏差且满足:re=r-rd,rd为期望回转率,r为转艏角速度/回转率;rmax为回转率r的安全限界,ρ1(re,rd)满足:
Figure BDA0003287118590000041
5.步骤5中纵向推力控制律具体为:
Figure BDA0003287118590000042
其中,τP为空气螺旋桨产生的纵向推进力,FxD表示纵荡自由度中除去τP的合力,Du为 纵荡自由度上的动力学不确定性,
Figure BDA0003287118590000043
为Du的观测值,
Figure BDA0003287118590000044
和η2均为大于零的常数,m为全垫 升气垫船的质量,v表示船体坐标系下的全垫升气垫船的横向速度,r表示船体坐标系下的全 垫升气垫船的转艏角速度/回转率;
Figure BDA0003287118590000045
为期望纵向速度ud的一阶导数,su1为选取的线性滑模 面,满足:su1=λ2ue,其中,λ2为滑模面增益且大于零,|ue|=|u-ud|,
Figure BDA0003287118590000046
Figure BDA0003287118590000047
umin为最小速度。
本发明的有益效果:与现有技术相比,本发明采用的SFLOS导引算法既能跟踪直线路径 又能跟踪曲线路径,解决了常规LOS导引算法只能跟踪直线路径的局限性;本发明将基于 SFLOS导引算法的带回转率约束的转艏力矩控制器以及带侧滑角约束的纵向推力控制器应用 到全垫升气垫船这种比较特殊的对象模型的路径跟踪过程当中,不但解决了全垫升气垫船在 路径跟踪回转过程中容易发生的回转率以及侧滑角分别超出其各自的安全限界的问题;而且, 所提出的基于SFLOS导引算法的带回转率和侧滑角约束的路径跟踪控制器能够使全垫升气 垫船路径跟踪过程中的纵向位置误差以及横向位置误差分别渐进收敛于零。
附图说明
图1带状态约束的全垫升气垫船路径跟踪原理图;
图2全垫升气垫船跟踪期望路径以及实际运动路径对比图;
图3全垫升气垫船路径跟踪的纵向位置误差与横向位置误差对比图;
图4βmax=5.5°时的全垫升气垫船运动的侧滑角对比图;
图5全垫升气垫船运动的实际纵向速度对比图;
图6全垫升气垫船运动的航迹角对比图;
图7rmax=1.2°/s时的全垫升气垫船运动的实际回转率对比图;
图8外界干扰的不确定性的实际值及其观测值。
具体实施方式
下面结合说明书附图和具体实施方式对本发明做进一步说明。
结合图1,本发明实施例具体包括:
步骤1:建立包括全垫升气垫船四自由度运动学数学模型和全垫升气垫船四自由度动力 学数学模型的全垫升气垫船四自由度运动数学模型,并在所建模型的基础上得到全垫升气垫 船航迹向控制数学模型以及纵向速度控制数学模型,为之后推导纵向推力控制律以及转艏力 矩控制律做准备。
步骤1.1:建立形如微分方程组形式的包括纵荡、横荡、横倾、艏摇四个自由度的全垫升 气垫船四自由度运动数学模型为:
Figure BDA0003287118590000051
其中,x、y、φ、ψ分别表示北东坐标系下的全垫升气垫船的北向位置、东向位置、横倾角、艏向角,u、v、p、r分别表示船体坐标系下的全垫升气垫船的纵向速度、横向速度、 横倾角速度、转艏角速度/回转率;其中,
Figure BDA0003287118590000052
Figure BDA0003287118590000053
其中,m为全垫升气垫船的质量,Fx、Fy、Mx、Mz分别为全垫升气垫船的纵向合力、 横向合力、横倾合力矩、转艏合力矩,Ix、Iz分别为全垫升气垫船绕x、z轴的转动惯量,τP为空气螺旋桨产生的纵向推进力,也是后面要设计的纵向推力控制律,τR为空气舵产生的转 艏力矩,也是后面要设计的转艏力矩控制律;各个分力(矩)当中包含的下角标a、m、h、 R分别表示空气动力(矩)、空气动量力(矩)、水动力(矩)、空气舵力(矩)在纵荡、横荡、 横倾、艏摇四个自由度上的分力(矩);FxD、FyD、MxD、MzD分别为纵荡、横荡、横倾、 艏摇四个自由度中除去τP和τR(如果有)的其它合力;侧滑角β与纵向、横向速度之间存在 的如下的关系:
Figure BDA0003287118590000061
步骤1.2:在上述所建模型的基础上并考虑到全垫升气垫船所受到的外界不确定性干扰, 得到全垫升气垫船航迹向控制数学模型以及纵向速度控制数学模型分别如下:
Figure BDA0003287118590000062
Figure BDA0003287118590000063
其中,ψw为实际航迹角,β为侧滑角且假设为可测的量,d3为外界不确定性干扰在转艏 自由度上的分量,Dr为转艏自由度上的动力学不确定性,d1为外界不确定性干扰在纵荡自由 度上的分量,Du为纵荡自由度上的动力学不确定性。
步骤2:设计SFLOS的曲线路径导引算法,给出路径参数更新率
Figure BDA0003287118590000064
和期望航迹角ψwd,SFLOS导引算法的具体设计流程如下:
步骤2.1:跟踪误差转换:定义全垫升气垫船期望路径P(θ)为:[xθ,yθ]T,其中,θ为路 径参数。SF坐标系下位置误差向量为:[xe,ye]T,北东坐标系下位置误差向量为: [x-xθ,y-yθ]T。两者满足如下等式:
Figure BDA0003287118590000065
其中,ψθ为北东坐标系向SF坐标系转换时的旋转角度,且ψθ=arctan2(y′θ,x′θ)。
若全垫升气垫船跟踪上期望路径,只需设计导引律保证[xe,ye]T渐进收敛于零即可。
步骤2.2:对SF坐标系下的跟踪误差xe和ye求导可得:
Figure BDA0003287118590000071
步骤2.3:构造李雅普诺夫函数如下:
Figure BDA0003287118590000072
对Vε求导可得:
Figure BDA0003287118590000073
步骤2.4:为使xe能渐进收敛于零,根据步骤2.3中的
Figure BDA0003287118590000074
可设计路径参数更新率
Figure BDA0003287118590000075
如下:
Figure BDA0003287118590000076
其中,kx为大于零的设计参数。将设计的
Figure BDA0003287118590000077
带入
Figure BDA0003287118590000078
可得:
Figure BDA0003287118590000079
步骤2.5:为使ye能渐进收敛于零,将全垫升气垫船的航迹角ψw看作虚拟控制量,设计 虚拟控制律ψwd如下:
Figure BDA00032871185900000710
其中,Δ为大于零的目标跟踪点P(θ)处沿设定路径切向的前视距离参数。
步骤2.6:设计动态前视距离参数如下:
Figure BDA00032871185900000711
其中,μ、Δmax、Δmin均为大于零的设计参数,且Δmax>Δmin;可通过调节μ来调节Δ的 收敛速率;选取Δ为时变量,Δ能自适应变化,当全垫升气垫船与期望路径相距较远时,ye比 较大,Δ趋近于Δmin,Δ取值较小,全垫升气垫船能快速接近期望路径,ye能快速缩小;当 全垫升气垫船航行至期望路径附近时,ye比较小,Δ趋近于Δmax,Δ取值较大,ye缩小速度 较慢,有助于避免ye超调。
步骤2.7:镇定航迹角误差。定义虚拟控制量ψwd的跟踪误差如下:
Figure BDA0003287118590000081
则有:
Figure BDA0003287118590000082
步骤2.8:将ψw代入
Figure RE-GDA0003385510030000083
中(3-19)可得:
Figure BDA0003287118590000084
若设计控制律使
Figure BDA0003287118590000085
能渐进收敛于零,即ψwwd能渐进收敛于零,即可得到:
Figure BDA0003287118590000086
综上可知,SFLOS导引算法的具体表达式为:
Figure BDA0003287118590000087
其中,
Figure BDA0003287118590000088
为路径参数更新率,ψwd为期望航迹角,U为实际合速度,ψθ为北东坐标系向SF坐标系转换时的旋转角度,且ψθ=arctan2(y′θ,x′θ);xθ、yθ为期望路径P(θ)上的任意路径点的坐标,满足下式:
Figure BDA0003287118590000089
其中,
Figure BDA00032871185900000810
kx为大于零的设计参数;Δ为P(θ)处沿设定路径切向的大于零的前视 距离参数;
步骤3:依据步骤1.2的航迹向控制数学模型和纵向速度控制数学模型,分别设计扩张状 态观测器来观测外界不确定性干扰,并分别在后续转艏力矩控制律和纵向推力控制律的设计 过程中对外界干扰的不确定性进行补偿。。
步骤3.1:针对带有Dr的全垫升气垫船航迹向控制数学模型,设计扩张状态观测器如下:
Figure BDA0003287118590000091
其中,
Figure BDA0003287118590000092
为ψw的观测值,
Figure BDA0003287118590000093
为ψw的观测误差,
Figure BDA0003287118590000094
为r的观测值,
Figure BDA0003287118590000095
为Dr的观测值,λ11、 λ12、λ13均为大于零的设计参数,fal(e,α,δ)为具有如下形式的连续非线性幂次函数:
Figure BDA0003287118590000096
此处的e即为
Figure BDA0003287118590000097
δ为大于零的设计参数,表示fal(e,α,δ)在原点附近线性段的长度,α 为大于零且小于一的设计参数,sgn(·)为符号函数。
步骤3.2:针对带有Du的全垫升气垫船纵向速度控制数学模型,设计扩张状态观测器如 下:
Figure BDA0003287118590000098
其中,
Figure BDA0003287118590000099
为u的观测值,eu为u的观测误差,
Figure BDA00032871185900000910
为Du的观测值,λ21、λ22均为大于零的设计参数,此处的e即为eu
步骤4:基于积分型障碍李亚普诺夫函数IBLF对路径跟踪过程中的回转率进行约束来设 计转艏力矩控制律,不但能使全垫升气垫船的实际航迹角跟踪上由SFLOS导引律生成的期望 航迹角,进而使全垫升气垫船的横向位置误差能渐进收敛于零;而且,保证了全垫升气垫船 路径跟踪过程中的回转率在回转率安全限界以内。
步骤4.1:定义航迹角偏差如下:
Figure BDA00032871185900000911
其中,ψwd为期望航迹角;基于全垫升气垫船航迹向控制模型可得:
Figure BDA0003287118590000101
步骤4.2:取二阶滑模面如下:
Figure BDA0003287118590000102
其中,
Figure BDA0003287118590000103
为滑模面增益且为大于零的常数,则
Figure BDA0003287118590000104
步骤4.3:取等速趋近律如下:
Figure BDA0003287118590000105
其中,K2为常数且大于零,sgn(·)为符号函数,则
Figure BDA0003287118590000106
可设计虚拟控制律rd如下:
Figure BDA0003287118590000107
步骤4.4:由于
r=re+rd
将rd带入到r中,再将r带入到
Figure BDA0003287118590000108
中可得:
Figure BDA0003287118590000109
步骤4.5:定义回转率偏差如下:
re=r-rd
其中,rd为期望回转率;对上式等号两边进行求导并基于步骤1.3中的全垫升气垫船航 迹向控制模型可得:
Figure BDA00032871185900001010
选取IBLF为:
Figure BDA00032871185900001011
其中,rmax为回转率r的安全限界,即
|r|<rmax
步骤4.6:对V3求导可得:
Figure BDA0003287118590000111
其中,选取γ为积分变量,且根据分部积分法及用σ=γ·re进行替换可知:
Figure BDA0003287118590000112
可设计转艏力矩控制律如下:
Figure BDA0003287118590000113
其中,
Figure BDA0003287118590000114
为外界干扰Dr的观测值,
Figure BDA0003287118590000115
K3均为大于零的常数。
步骤4.7进行稳定性分析,取Lapunov函数如下:
Figure BDA0003287118590000116
对V2等号两边进行求导可得:
Figure BDA0003287118590000117
由上式可以看出,如果re
Figure BDA0003287118590000118
均能渐进收敛于零,或者re
Figure BDA0003287118590000119
均能收敛于零的一个很小 的邻域内,则
Figure BDA00032871185900001110
Figure BDA00032871185900001111
将渐进收敛于零,故ψwe
Figure BDA00032871185900001112
将渐进收敛于零;
将τR代入
Figure BDA00032871185900001113
中可得:
Figure BDA0003287118590000121
其中,
Figure BDA0003287118590000122
为外界干扰的观测误差,
Figure BDA0003287118590000123
Figure BDA0003287118590000124
的界值,即
Figure BDA0003287118590000125
因此,设计增益
Figure BDA0003287118590000126
可使
Figure BDA0003287118590000127
故re将渐进收敛于零,且τR的设计过程中运用到指数趋近律的思想,故re的收敛于零的过程中是平滑且基本无抖振的,则
Figure BDA0003287118590000128
将收敛于零的一个很小的邻域内;进一步,
Figure BDA0003287118590000129
恒成立,则
Figure BDA00032871185900001210
将渐进收敛于零,故ψwe
Figure BDA00032871185900001211
将渐进收敛于零。
步骤5:运用对称型对数型障碍李亚普诺夫函数LBLF对路径跟踪过程中的侧滑角进行 约束来设计纵向推力控制律,不但使全垫升气垫船的实际纵向速度跟踪上期望纵向速度,而 且可以保证全垫升气垫船路径跟踪过程中的侧滑角在侧滑角安全限界以内。同时所设计的纵 向推力控制律既能够保证全垫升气垫船避开阻力峰速附近的航速范围内,以防止全垫升气垫 船出现不稳定的航行状态,而且也能够防止全垫升气垫船在回转过程中出现失速现象。
步骤5.1:将侧滑角约束转化为纵向速度约束。根据侧滑角β与纵向、横向速度之间存在 的关系可求出:
Figure BDA00032871185900001212
设侧滑角满足如下条件:
Figure BDA00032871185900001213
其中,βmax为侧滑角安全限界。
步骤5.2:由
Figure BDA00032871185900001214
可得:
Figure BDA0003287118590000131
步骤5.3:基于纵向速度约束和期望速度进行纵向速度规划。假设初始纵向速度u(0)满足 下式:
umin(0)<u(0)<2ud(0)-umin(0)
其中,umin(0)为航行所要达到的最小速度umin的初值,umin(0)的表达形式如下:
Figure BDA0003287118590000132
umin(0)为航行所要达到的最小速度umin的初值;
初始期望纵向速度ud(0)满足下式:
ud(0)>umin(0)
步骤5.4:设计纵向速度的规划函数如下:
Figure BDA0003287118590000133
其中,
Figure BDA0003287118590000134
ku、δ、εu min均为大于零的常数。选取这种时变的ud,再通过选取合适的参数,可使ud>umin恒成立。
步骤5.5:基于LBLF设计纵向推力控制律。定义纵向速度偏差的绝对值如下:
|ue|=|u-ud|
定义LBLF边界函数的相关变量如下:
Figure BDA0003287118590000135
设ue的绝对值满足如下条件:
Figure BDA0003287118590000136
步骤5.6取线性滑模面如下:
su1=λ2ue
其中,λ2为滑模面增益且大于零,则
Figure BDA0003287118590000141
步骤5.7:选取LBLF如下:
Figure BDA0003287118590000142
Figure BDA0003287118590000143
则上式可以改写为:
Figure BDA0003287118590000144
对上式求导并结合步骤1.3中的纵向速度控制数学模型可得:
Figure BDA0003287118590000145
由上式可设计纵向推力控制律如下:
Figure BDA0003287118590000146
其中,
Figure BDA0003287118590000147
为Du的观测值,
Figure BDA0003287118590000148
和η2均为大于零的常数。
步骤5.8:进行稳定性分析,将τP代入到
Figure BDA0003287118590000149
中可得:
Figure BDA00032871185900001410
其中,
Figure BDA00032871185900001411
Figure BDA00032871185900001412
的界值,即
Figure BDA00032871185900001413
取增益
Figure BDA00032871185900001414
则可使
Figure BDA00032871185900001415
则ue将渐进收敛 于零。
给出与基于SFLOS导引算法的带回转率和侧滑角约束的路径跟踪控制器的作对比的基 于SFLOS导引算法的不带状态约束的路径跟踪控制器的表达式,其具体形式如下:
Figure BDA0003287118590000151
Figure BDA0003287118590000152
Figure BDA0003287118590000153
其中,
Figure BDA0003287118590000154
sr、su均为二阶滑模面;rd为虚拟控制律,τR为转艏力矩控制律,τp为纵向推力控制律;
Figure BDA0003287118590000155
K2
Figure BDA0003287118590000156
K3,λ1,K1均为大于零的常数。
附图中的仿真图部分主要是针对基于SFLOS导引算法的带回转率以及侧滑角约束的全 垫升气垫船路径跟踪仿真图和基于SFLOS导引算法的不带状态约束的全垫升气垫船路径跟 踪仿真图之间的对比。
从图2可以看出,带回转率以及侧滑角约束的全垫升气垫船路径跟踪控制器和不带状态 约束的全垫升气垫船路径跟踪控制器(后面简称“带状态约束的和不带状态约束的”)均能够 实现对期望正弦曲线路径的跟踪。从图3可以看出,在371秒左右,不带状态约束的纵向位 置偏差几乎为零,带状态约束的纵向位置偏差的绝对值约为0.2米;从350秒至450秒,不 带状态约束的横向位置偏差几乎为零,带状态约束的横向位置偏差有比较大的波动,说明不 带状态约束的跟踪正弦曲线路径的效果更好。从图4、图5、图7可以看出,带状态约束的全 垫升气垫船的侧滑角、实际纵向速度、回转率均分别没超出各自的界限值βmax=5.5°、umin、 rmax=1.2°/s,不带状态约束的全垫升气垫船的侧滑角、实际纵向速度、回转率均在70秒左右 分别超出各自的界限值,说明了所设计的带状态约束的控制器的有效性。从图6可以看出, 不带状态约束和带状态约束的航迹角控制器均能够分别实现对SFLOS导引算法所生成的期 望航迹角的跟踪。从图8可以看出,针对外界不确定性干扰在纵荡自由度以及转艏自由度上 的分量,所设计的扩张状态观测器均能够较好的观测外扰动。
图5中,带侧滑角和回转率约束与带侧滑角和回转率约束的umin是一组,不带状态约束 与不带状态约束的umin是一组;且从图5可以看出,带侧滑角和回转率约束的umin始终小于带 侧滑角和回转率约束的全垫升气垫船的实际纵向速度,不带状态约束的umin在100秒左右远 超过不带状态约束的全垫升气垫船的实际纵向速度。

Claims (6)

1.一种全垫升气垫船路径跟踪方法,其特征在于,包括以下步骤:
步骤1:建立全垫升气垫船四自由度运动数学模型,并在所建模型的基础上得到全垫升气垫船航迹向控制数学模型和纵向速度控制数学模型;
步骤2:设计SFLOS的曲线路径导引算法,得到路径参数更新率
Figure FDA0003287118580000011
和期望航迹角ψwd
步骤3:依据步骤1所得的航迹向控制数学模型和纵向速度控制数学模型分别设计扩张状态观测器观测外界干扰的不确定性;
步骤4:运用积分型障碍李亚普诺夫函数设计转艏力矩控制律对路径跟踪控制中的回转率进行约束,使得路径跟踪过程中的回转率在安全限界内,全垫升气垫船的实际航迹角跟踪期望航迹角ψwd
步骤5:运用对称型对数型障碍李亚普诺夫函数设计纵向推力控制律对路径跟踪控制中的侧滑角进行约束,使得路径跟踪过程中的侧滑角在安全限界内,全垫升气垫船的实际纵向速度跟踪期望纵向速度。
2.根据权利要求1所述的一种全垫升气垫船路径跟踪安全控制方法,其特征在于:步骤1所述全垫升气垫船航迹向控制数学模型和纵向速度控制数学模型具体为:
Figure FDA0003287118580000012
Figure FDA0003287118580000013
其中,ψw为实际航迹角,β为侧滑角且假设为可测的量,d3为外界不确定性干扰在转艏自由度上的分量,Dr为转艏自由度上的动力学不确定性,d1为外界不确定性干扰在纵荡自由度上的分量,Du为纵荡自由度上的动力学不确定性;r为船体坐标系下的全垫升气垫船的转艏角速度/回转率;φ表示北东坐标系下的全垫升气垫船的横倾角;τR为空气舵产生的转艏力矩;τP为空气螺旋桨产生的纵向推进力;v表示船体坐标系下的全垫升气垫船的横向速度;MzD表示艏摇自由度中除去τP和τR的合力,FxD表示纵荡自由度中除去τP和τR的合力。
3.根据权利要求1或2所述的一种全垫升气垫船路径跟踪方法,其特征在于:步骤2所述路径参数更新率
Figure FDA0003287118580000021
和期望航迹角ψwd满足:
Figure FDA0003287118580000022
其中,ψθ为北东坐标系向SF坐标系转换时的旋转角度,且ψθ=arctan2(y′θ,x′θ),xθ、yθ为期望路径P(θ)上的任意路径点的坐标,满足:
Figure FDA0003287118580000023
其中,
Figure FDA0003287118580000024
kx为大于零的设计参数;Δ为P(θ)处沿设定路径切向的大于零的前视距离参数,Δ的具体形式如下:
Figure FDA0003287118580000025
其中,μ、Δmax、Δmin均为大于零的设计参数,且Δmax>Δmin,可通过调节μ来调节Δ的收敛速率。
4.根据权利要求1或2所述的一种全垫升气垫船路径跟踪方法,其特征在于:步骤3所述扩张状态观测器具体为:
针对航迹向控制数学模型,扩张状态观测器具体为:
Figure FDA0003287118580000026
其中,
Figure FDA0003287118580000027
为ψw的观测值,
Figure FDA0003287118580000028
为ψw的观测误差,
Figure FDA0003287118580000029
为r的观测值,
Figure FDA00032871185800000210
为Dr的观测值,λ11、λ12、λ13均为大于零的设计参数,fal(e,α,δ)为具有如下形式的连续非线性幂次函数:
Figure FDA00032871185800000211
e即为
Figure FDA00032871185800000212
δ为大于零的设计参数,表示fal(e,α,δ)在原点附近线性段的长度,α为大于零且小于1的设计参数,sgn(·)为符号函数;
针对纵向速度控制数学模型,设计扩张状态观测器具体为:
Figure FDA0003287118580000031
其中,
Figure FDA0003287118580000032
为u的观测值,eu为u的观测误差,
Figure FDA0003287118580000033
为Du的观测值,λ21、λ22均为大于零的设计参数,e即为eu
5.根据权利要求1或2所述的一种全垫升气垫船路径跟踪方法,其特征在于:步骤4所述转艏力矩控制律具体为:
Figure FDA0003287118580000034
其中,τR为空气舵产生的转艏力矩,MzD表示艏摇自由度中除去τR的合力,Iz为z轴的转动惯量,Dr为转艏自由度上的动力学不确定性,
Figure FDA0003287118580000035
为Dr的观测值,
Figure FDA0003287118580000036
K3均为大于零的常数,re为回转率偏差且满足:re=r-rd,rd为期望回转率,r为转艏角速度/回转率;rmax为回转率r的安全限界,ρ1(re,rd)满足:
Figure FDA0003287118580000037
6.根据权利要求1或2所述的一种全垫升气垫船路径跟踪方法,其特征在于:步骤5所述纵向推力控制律具体为:
Figure FDA0003287118580000038
其中,τP为空气螺旋桨产生的纵向推进力,FxD表示纵荡自由度中除去τP的合力,Du为纵荡自由度上的动力学不确定性,
Figure FDA0003287118580000039
为Du的观测值,
Figure FDA00032871185800000310
和η2均为大于零的常数,m为全垫升气垫船的质量,v表示船体坐标系下的全垫升气垫船的横向速度,r表示船体坐标系下的全垫升气垫船的转艏角速度/回转率;
Figure FDA00032871185800000311
为期望纵向速度ud的一阶导数,su1为选取的线性滑模面,满足:su1=λ2ue,其中,λ2为滑模面增益且大于零,|ue|=|u-ud|,
Figure FDA0003287118580000041
Figure FDA0003287118580000042
umin为最小速度。
CN202111151049.3A 2021-09-29 2021-09-29 一种全垫升气垫船路径跟踪方法 Active CN113867352B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202111151049.3A CN113867352B (zh) 2021-09-29 2021-09-29 一种全垫升气垫船路径跟踪方法
GB2210207.3A GB2618860A (en) 2021-09-29 2022-01-14 Path tracking method for air cushion vehicle
PCT/CN2022/070350 WO2023050636A1 (zh) 2021-09-29 2022-01-14 一种全垫升气垫船路径跟踪方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111151049.3A CN113867352B (zh) 2021-09-29 2021-09-29 一种全垫升气垫船路径跟踪方法

Publications (2)

Publication Number Publication Date
CN113867352A true CN113867352A (zh) 2021-12-31
CN113867352B CN113867352B (zh) 2024-09-06

Family

ID=78992747

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111151049.3A Active CN113867352B (zh) 2021-09-29 2021-09-29 一种全垫升气垫船路径跟踪方法

Country Status (2)

Country Link
CN (1) CN113867352B (zh)
WO (1) WO2023050636A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023050636A1 (zh) * 2021-09-29 2023-04-06 哈尔滨工程大学 一种全垫升气垫船路径跟踪方法
GB2618860A (en) * 2021-09-29 2023-11-22 Univ Harbin Eng Path tracking method for air cushion vehicle

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116753961B (zh) * 2023-08-16 2023-10-31 中国船舶集团有限公司第七〇七研究所 基于状态观测的动力定位船舶高速循迹导航方法
CN117539267B (zh) * 2024-01-09 2024-04-23 哈尔滨工程大学三亚南海创新发展基地 一种气垫船的运动控制方法及系统

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618611A (zh) * 2019-10-23 2019-12-27 哈尔滨工程大学 一种基于回转率约束的无人艇轨迹跟踪安全控制方法
CN113156965A (zh) * 2021-04-30 2021-07-23 哈尔滨工程大学 一种基于纵向速度规划的气垫船高速回转控制方法
WO2023050636A1 (zh) * 2021-09-29 2023-04-06 哈尔滨工程大学 一种全垫升气垫船路径跟踪方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE404682B (sv) * 1974-12-11 1978-10-23 Kockums Automation Forfarande och anordning for vesentligen kinematisk styrning av ett fartyg
WO2008110519A1 (de) * 2007-03-09 2008-09-18 Continental Teves Ag & Co. Ohg Automatische stabilisierungseinheit für wasserfahrzeuge
CN105825714B (zh) * 2016-05-24 2018-04-10 哈尔滨工程大学 鲁棒航迹引导律的确定方法及装置
CN107450318B (zh) * 2017-08-21 2020-07-28 哈尔滨工程大学 一种基于二阶滑模控制的气垫船路径跟踪控制方法
CN108037662B (zh) * 2017-12-06 2019-11-08 浙江工业大学 一种基于积分滑模障碍李雅普诺夫函数的四旋翼飞行器输出受限反步控制方法
CN108594639B (zh) * 2018-03-28 2020-12-22 哈尔滨工程大学 一种基于强化学习的全垫升气垫船航迹跟踪控制方法
CN112346465B (zh) * 2020-11-27 2022-09-02 哈尔滨工程大学 基于ialos导引律的欠驱无人船自适应模糊控制方法
CN113189979B (zh) * 2021-04-02 2023-12-01 大连海事大学 一种无人船的分布式队列有限时间控制方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110618611A (zh) * 2019-10-23 2019-12-27 哈尔滨工程大学 一种基于回转率约束的无人艇轨迹跟踪安全控制方法
CN113156965A (zh) * 2021-04-30 2021-07-23 哈尔滨工程大学 一种基于纵向速度规划的气垫船高速回转控制方法
WO2023050636A1 (zh) * 2021-09-29 2023-04-06 哈尔滨工程大学 一种全垫升气垫船路径跟踪方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王豫洛: "无人艇路径跟踪自抗扰控制方法研究", 《中国优秀硕士论文全文数据库(工程科技Ⅱ辑)》, 15 May 2021 (2021-05-15) *
陈光祥;李仲彬;刘伟;: "基于反步法的全垫升气垫船航向增稳控制", 自动化技术与应用, no. 11, 25 November 2018 (2018-11-25) *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023050636A1 (zh) * 2021-09-29 2023-04-06 哈尔滨工程大学 一种全垫升气垫船路径跟踪方法
GB2618860A (en) * 2021-09-29 2023-11-22 Univ Harbin Eng Path tracking method for air cushion vehicle

Also Published As

Publication number Publication date
CN113867352B (zh) 2024-09-06
WO2023050636A1 (zh) 2023-04-06

Similar Documents

Publication Publication Date Title
CN113867352A (zh) 一种全垫升气垫船路径跟踪方法
CN111580523B (zh) 一种基于侧滑角补偿的无人艇路径跟踪自抗扰控制方法
CN109189071B (zh) 基于模糊观测器的鲁棒自适应无人船路径跟踪控制方法
CN111123967B (zh) 一种基于自适应动态逆的固定翼无人机着舰控制方法
CN106292287B (zh) 一种基于自适应滑模控制的uuv路径跟踪方法
CN106950974B (zh) 欠驱动自主水下航行器的对三维路径进行理解及跟踪控制方法
Zheng et al. Adaptive sliding mode relative motion control for autonomous carrier landing of fixed-wing unmanned aerial vehicles
CA2569227C (en) Systems and methods for controlling dynamic systems
CN109782785B (zh) 基于直接力控制的飞机自动着舰控制方法
CN113341953B (zh) 一种吊舱式无人艇航向控制方法
CN106444822B (zh) 一种基于空间矢量场制导的平流层飞艇路径跟踪控制方法
CN110618611A (zh) 一种基于回转率约束的无人艇轨迹跟踪安全控制方法
CN107844123A (zh) 一种非线性飞行器航迹控制方法
CN106200665A (zh) 携带不确定负载的四轴飞行器的建模与自适应控制方法
Wang et al. Horizontal tracking control for AUV based on nonlinear sliding mode
CN114488803A (zh) 基于事件触发的无人船/机鲁棒自适应神经协同控制方法
CN112527018A (zh) 一种欠驱动自主水下航行器的三维镇定控制方法
CN105334854A (zh) 应用于气垫船航向控制与横倾控制的解耦控制装置及方法
Zhang et al. Line-of-sight path following control on UAV with sideslip estimation and compensation
CN115113524A (zh) 一种基于干预lvs制导的asv多端口事件触发路径跟踪控制方法
CN113110527B (zh) 一种自主水下航行器有限时间路径跟踪的级联控制方法
CN117234086A (zh) 一种全垫升气垫船航速航向指定时间快速控制方法
CN116500890A (zh) 基于rleso的气垫船滤波反步轨迹跟踪控制方法
CN116679693A (zh) 融合螺旋桨水动力的无人艇自抗扰控制方法
CN116203840A (zh) 可重复使用运载器自适应增益调度控制方法

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