CN111722645B - 模型参数已知的带负载四旋翼无人机速度控制方法 - Google Patents
模型参数已知的带负载四旋翼无人机速度控制方法 Download PDFInfo
- Publication number
- CN111722645B CN111722645B CN202010628316.0A CN202010628316A CN111722645B CN 111722645 B CN111722645 B CN 111722645B CN 202010628316 A CN202010628316 A CN 202010628316A CN 111722645 B CN111722645 B CN 111722645B
- Authority
- CN
- China
- Prior art keywords
- unmanned aerial
- aerial vehicle
- load
- quad
- controller
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000013461 design Methods 0.000 claims abstract description 28
- 239000011159 matrix material Substances 0.000 claims abstract description 23
- 238000013519 translation Methods 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 7
- 239000013598 vector Substances 0.000 claims description 31
- 239000000725 suspension Substances 0.000 claims description 15
- 230000004888 barrier function Effects 0.000 claims description 4
- 125000004432 carbon atom Chemical group C* 0.000 claims description 4
- 230000000694 effects Effects 0.000 claims description 4
- 230000009471 action Effects 0.000 claims description 3
- 239000003795 chemical substances by application Substances 0.000 claims description 3
- 230000008878 coupling Effects 0.000 claims description 3
- 238000010168 coupling process Methods 0.000 claims description 3
- 238000005859 coupling reaction Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 230000008859 change Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000005484 gravity Effects 0.000 description 4
- 238000002474 experimental method Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 238000011217 control strategy Methods 0.000 description 2
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/10—Simultaneous control of position or course in three dimensions
- G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种模型参数已知的带负载四旋翼无人机速度控制方法,包括步骤1)建立带负载四旋翼无人机的动力学模型,2)设定对无人机的速度跟踪误差ev=v‑v*,ew=w‑w*,3)在无人机上加装传感器测量飞行过程中的摆动角α和β,计算时变增益矩阵B,4)设计模型参数已知情况下带负载四旋翼无人机的平移子系统和旋转子系统的控制器。本发明考虑了带负载无人机在时变增益矩阵、执行器饱和与执行器故障等实际条件下的飞行速度控制问题,解决了带负载后增益矩阵时变下的控制器设计问题,所提出的控制器不仅可以保证闭环系统中所有内部信号有界且连续,并且确保最终误差收敛于一个小的紧集,而且可以在取得良好的轨迹跟踪性能过程中确保不违反状态约束条件。
Description
技术领域
本发明涉及无人机控制技术领域,特别涉及一种带负载的四旋翼无人机的速度控制方法。
背景技术
旋翼式无人机的动力系统是由多个旋翼提供推力的。例如四旋翼无人机,当四个带电机的螺旋桨产生的推力大于无人机自身的重力和空气阻力的时候,无人机就能够起飞。通过控制系统反馈状态信息,调整每个电机的转速来改变每个旋翼的转速来进行无人机的姿态飞行控制。
在实际应用中,四旋翼无人机不可避免地需要携带任务设备或相机来执行特定任务。通常,四旋翼无人机一般采用夹持和悬挂两种方法。例如,航拍摄影的时候,就必须采用夹持方式,这是为了防止四旋翼无人机因为大风等外部环境原因飞行不稳定导致相机产生强烈的抖动,影响拍摄效果;当检测地雷时,又必须采用悬挂方法,只有绳索足够长才能确保四旋翼无人机的安全和完成既定任务。因此,研究四旋翼无人机携带载荷具有重要的现实意义。
目前,许多的研究集中在四旋翼无人机的建模和控制上。一般来说,四旋翼无人机的6自由度动力学模型可以通过使用拉格朗日方法或欧拉方程来推导。但是,大多数研究都采用了经典四旋翼无人机模型,其中假设无人机质心是固定的并且是已知的,即质心与无人机机身几何中心处于同一位置。实际上,在携带载荷以后,整个四旋翼无人机运载系统的质心通常都会发生一定的偏移,导致质心与机身几何中心不在同一位置,同时转动惯量也会发生变化,这些不确定性因素都会影响四旋翼无人机控制系统的稳定性。因此,给无人机的数学建模工作和控制策略设计带来艰巨性、复杂性和挑战性。
由于四旋翼无人机始终需要在动态环境中保证连续稳定飞行,控制器所需求的控制量可能超过无人机执行器所提供的最大控制量,这就会造成无人机控制系统性能急速下降,无法正常完成指定任务,甚至引发控制系统的崩溃。同时,无人机也可能发生执行器、状态、传感器等故障,也将导致性能显著下降,甚至发生坠机等灾难性后果。此外,四旋翼无人机在实际飞行中存在许多物理条件约束问题,例如位置约束、速度约束、性能、安全规格等。
因此,在四旋翼无人机的控制系统设计中,同时考虑带负载、无人机质心时变、执行器饱和、执行器故障和状态约束条件更加具有实际价值。
发明内容
有鉴于此,本发明的目的是提供一种模型参数已知的带负载四旋翼无人机速度控制方法,以解决在同时考虑带负载、受负载摆动影响情况下系统模型中增益矩阵时变、执行器饱和、执行器故障和状态约束条件下的四旋翼无人机速度控制问题。
本发明模型参数已知的带负载四旋翼无人机速度控制方法,其包括以下步骤:
1)建立带负载四旋翼无人机的动力学模型,具体表示如下:
上式中,fa(·)和fb(·)代表模型的耦合项,A=diag{m+m0;m+m0;m+m0}和表示增益矩阵,m表示无人机的质量,m0表示负载的质量,M1表示惯性矩阵,m1=m0zh(zh+lcosα),m2=m0lzh sinαsinβ,m3=-m0lzhsinαcosβ,其中l是无人机悬挂负载的缆绳的长度,α是缆绳长度l与Zb负方向之间的角度,β是l投影到机体坐标系XbObYb平面的投影线和YbObZb平面之间的角度,zh是负载悬挂点到机体坐标系XbObYb原点的距离;J为在机体坐标系下四旋翼无人机的转动惯量,J对于每一台无人机来说是固定的,可通过查产品手册直接计算得到、或可直接在设计制造无人机时将其转动惯量给定在产品手册中、或者可以根据实验测得无人机的转动惯量;ut是速度控制量,ur是角速度控制输入;Fp1是负载施加在无人机上的额外力Fp中不含状态量和v的剩余项;Mp1是负载施加于无人机的额外力矩Mp中不含有状态变量和w的剩余项;da(·)和db(·)代表系统模型的不确定项;向量v=[vx,vy,vz]T表示在机体坐标系下的线速度,为v的一阶导数,为带负载四旋翼无人机平移子系统;向量w=[wx,wy,wz]T表示在机体坐标系下的角速度,是w的一阶导数,为带负载的四旋翼无人机的旋转子系统;机体坐标系表示为Ob=(Xb,Yb,Zb),机体坐标系的原点Ob取在四旋翼的质心位置上,Xb轴在四旋翼对称轴内指向机头方向,Zb轴在无人机对称平面内,Zb轴垂直于Xb轴指向上方,按右手定则确定Yb轴;
考虑无人机的物理结构和性能指标限制,无人机的速度和角速度都需要满足如下限制:|v|≤kv,|w|≤kw,其中kv是无人机的速度所能达到的最大极限值,kw是无人机的角速度极所能达到的最大极限值;
考虑输入饱和影响,由于转速限制,ui=[ut,ur],i=1,...,6;不再是设计的控制输入,系统模型表示如下:
此处,p(ui)表示受非对称非光滑的饱和非线性影响的控制输入,定义为
那么,p(ui)可以表示为其中Ψ(ui)是p(ui)和之差;因函数和饱和度函数p(ui)的有界性质,因此函数Ψ(ui)是有界的,即|Ψ(ui)|≤Ψm,其中Ψm是正定未知的常数;为方便起见,对函数采用均值定理,变成其中且0<λ<1;通过选择ui0=0并利用事实p(0)=0,得到由此引入以下变量:对于所有是不变的,并且存在正定常数gmax,使得因此,有0<gi≤gmax<∞;最后非线性系统(2)变成:
上式中,Lt1(·)=AΨ(ut)+Fp1/(m+m0)+da(·),Gt1=diag{g1,g2,g3},Lr1(·)=BΨ(ur)+M1 -1 Mp1+db(·),Gr1=diag{g4,g5,g6}。
考虑执行器故障,此时的控制输入不再是ui而是uai=ρi(t)ui+δi(t)(i=1,...,6),其中ρi(t)∈(0,1]表示执行器效率因子,δi(t)表示由部分控制动作产生的时变且不可测量的矢量函数;将式(4)的模型变成如下等式:
上式中,Gt=Gt1ρt,ρt=diag{ρ1,ρ2,ρ3},Lt2(·)=AGt1δt+Lt1(·),Gr=Gr1ρr,ρr=diag{ρ4,ρ5,ρ6},Lr2(·)=BGr1δr+Lr1(·);
2)设定对无人机的速度跟踪误差,具体包括:
设定平移跟踪误差为ev=v-v*,旋转跟踪误差为ew=w-w*,定义期望参考为和期望参考v*(t)和w*(t)是已知和有界的,即满足对任意的t≥0,有|v*|≤Av<kv和|w*|≤Aw<kw,其中kv和kw是无人机的速度和角速度极限值,Av=[Av1,Av2,Av3]T和Aw=[Aw1,Aw2,Aw3]T是小于极限值的正常数已知向量;
将式(5)转换成跟踪误差的动力学模型如下等式:
3)在无人机上加装传感器测量飞行过程中的摆动角α和β,计算时变增益矩阵B;
4)设计模型参数已知情况下带负载四旋翼无人机的平移子系统和旋转子系统的控制器,并通过设计的控制器控制无人机系统,具体包括:
第一步,考虑式(6)里面的平移子系统,设计基于模型参数已知的受状态约束平移控制器,并实现所需的跟踪控制目标:
为了确保不违反v的约束,即|v|<kv,根据BLF的性质,将李雅普诺夫函数的第一部分定义为
其中kb为正定常数向量,与此同时,定义一个紧凑集合Ωev={ev:|ev|<kb},满足V1在紧凑集合Ωev中是有效的;为了使不等式|v|<kv成立,kb被选定成:
kb=kv-Av (8)
由于ev=v-v*和|v*|≤Av,所以|v|≤|ev|+|v*|<kb+Av=kv-Av+Av=kv;
控制器ut的定义为
其中,σ1>0是控制器设计参数;
第二步,考虑式(6)误差跟踪动态的旋转子系统,设计基于模型参数已知的受状态约束姿态控制器,并实现所需的跟踪控制目标:
为了确保不违反角速度w的约束,即|w|<kw,根据势垒函数BLF的性质,将李雅普诺夫函数的第一部分定义为
其中,kc是控制器设计参数;同时,定义紧凑集合Ωew={ew:|ew|<kc},其满足V2在紧凑集合Ωew中有效,为了满足不等式|w|<kw,选择kc为
kc=kw-Aw (12)
ew=w-w*和|w*|≤Aw,有|w|≤|ew|+|w*|<kc+Aw=kw-Aw+Aw=kw;
控制器uri和更新律定义为:
本发明的有益效果:
1、本发明模型参数已知的带负载四旋翼无人机速度控制方法,以带负载的四旋翼无人机为研究对象,在同时考虑时变增益矩阵、执行器饱和与执行器故障等实际应用条件下,建立了带负载四旋翼无人机的6自由度动力学模型;并针对受状态约束的情况,设计了两种自适应容错控制器。理论分析和仿真结果表明,所提出的控制器不仅可以保证闭环系统中所有内部信号有界且连续,并且确保最终误差收敛于一个小的紧集,而且可以在取得良好的轨迹跟踪性能过程中确保不违反状态约束条件。
2、本发明模型参数已知的带负载四旋翼无人机速度控制方法,针对现有关于带负载四旋翼无人机的动力学模型研究默认无人机质心是固定的并且是已知的,即质心与无人机机身几何中心处于同一位置,忽略了负载摆动会使质心时变,进而系统模型中增益矩阵B也是时变的情况。本发明考虑了增益矩阵B时变的问题,可以更准确地反映无人机系统的动态过程,所设计的控制器能更好地满足无人机控制的实际需求。
3、本发明模型参数已知的带负载四旋翼无人机速度控制方法,针对实际飞行中经常遇到在执行器饱和情况时控制性能下降的问题,本发明通过定义良好的平滑函数,无需事先知道输入饱和度的界限即可处理具有未知不对称性和非平滑度的执行器输入饱和故障;针对实际飞行中受物理条件限制情况下的系统状态约束问题,本发明通过引入李雅普诺夫势垒函数和反步法,所设计的控制器能自适应容错和满足状态约束要求。
附图说明
图1为悬挂了载荷的四旋翼无人机动力学示意图。
图2为带状态约束的速度跟踪及误差vx。
图3为带状态约束的速度跟踪及误差vy。
图4为带状态约束的速度跟踪及误差vz。
图5为带状态约束的角速度跟踪及误差wx。
图6为带状态约束的角速度跟踪及误差wy。
图7为带状态约束的角速度跟踪及误差wz。
图8为带状态约束下四旋翼无人机控制输入。
图9为带状态约束下悬挂负载的摆动角。
具体实施方式
下面结合附图和实施例对本发明作进一步描述。
本实施例中模型参数已知的带负载四旋翼无人机速度控制方法,包括以下步骤:
1)建立带负载四旋翼无人机的动力学模型,包括:
为了精确的描述悬挂负载的四旋翼无人机的姿态和位置信息,定义如下坐标系:
①地球固联坐标系(e系)
也称地面坐标系,表示为Oe=(Xe,Ye,Ze)。将地心作为坐标原点Oe,Xe轴在水平面内指向某一方向,Ze轴垂直于地面向上。然后,按右手定则确定Ye轴。
②机体坐标系(b系)
也称体轴系,表示为Ob=(Xb,Yb,Zb)。与四旋翼无人机固连,其原点Ob取在四旋翼的质心位置上。Xb轴在四旋翼对称轴内指向机头方向。Zb轴在无人机对称平面内,Zb轴垂直于Xb轴指向上方。然后,按右手定则确定Yb轴。
③悬挂点坐标系(h系)
绳子的悬挂点与无人机质心位置是不一致的,为了准确描述悬挂负载的运动,需要建立悬挂点坐标系,表示为Oh=(Xh,Yh,Zh)。其原点Oh取四旋翼机体悬挂点位置。其它三轴的方向分别与四旋翼无人机机体坐标系的方向矢量平行。
定义ξ=[x,y,z]T表示四旋翼无人机在e系下的位置向量,η=[φ,θ,ψ]T表示四旋翼无人机在e系下的姿态向量,其中三个欧拉角分别表示横滚角φ,俯仰角θ和偏航角ψ。向量v=[vx,vy,vz]T表示在b系下的线速度,向量ve=[vex,vey,vez]T表示Ob在e系下的线速度,有向量w=[wx,wy,wz]T表示在b系下的角速度。位置矢量和线速度之间的转换关系为:
类似地,b系下的角速度与e系下的欧拉角之间的转换关系可以写成
假设机体坐标系在地面坐标系中以角速度w旋转,对于任意向量A,满足
上式中,m为四旋翼无人机的质量,向量F=[Fx,Fy,Fz]T表示作用在四旋翼无人机上xyz三个方向的分力。由力学分析表示,四旋翼无人机受到的合外力包括四旋翼无人机总升力和自身重力。因此,在机体坐标系下,四旋翼无人机的总升力FT在xyz三个方向的分力如下:
其中κ>0为升力系数,Ωi为第i(i=1,2,3,4)个螺旋桨的转速。在机体坐标系下,无人机的重力可表示为
因此,四旋翼无人机所受的合外力为:
F=FT+FG (9)
再由角动量定理,物体角动量的变化等于力矩的冲量,可得
Mdt=d(Jw) (10)
其中,向量M=[uφ,uθ,uψ]T是指输入转矩矢量,J为在机体坐标系下四旋翼无人机的转动惯量,且假定为常数,如下表示
将(式5)和(6)结合起来,可得到
上式(12)为不考虑悬挂负载的四旋翼无人机平移子系统的动力学模型。
通常,可以假设将转动惯量矩阵J的非对角项视为一个小数项,当有两个小数项相乘时,可以被忽略。此假设比普通无人机模型更仔细地考虑了四旋翼无人机的质量不对称性。所以,将式(5)和(10)结合起来,可得到
上式中,将J看作一个常数矩阵,使之计算较为简便。
结合式(12)和(13),四旋翼无人机的6-DOF动力学模型可以如下描述
上式中的系数矩阵以及其详细表达如下所示:
图1表示带有悬挂负载的四旋翼无人机动力学示意图。如果不考虑附加的有效载荷,无人机质心将始终与几何中心(即原点Ob)重合。但是,在实际飞行中,由于悬挂有效载荷施加摆动的影响,无人机质心将随时间变化,不再与原点Ob一致。将变化的质心表示为C点。由于将额外的有效负载视为质点P,因此可以忽略悬挂负载的空气阻力。将无人机和负载分开单独分析,先分析得出负载的运动模型,将负载施加在无人机的额外力和额外力矩视作与无人机本体独立的附加分量,直接导入四旋翼无人机数学模型中,就可以得到带负载的四旋翼无人机动力学模型。坐标系框架h的单位矢量始终保持与框架b的单位矢量平行。负载的位置信息由两个摆动角α和β描述,其中α是缆绳长度l与Zb负方向之间的角度,β是l投影到XbObYb平面的投影线和YbObZb平面之间的角度。根据刚体力学定理,负载在悬挂坐标系(h系)中的位置矢量由下式给出:
rp=[xp,yp,zp]T=[lsinαsinβ,-lsinαcosβ,-lcosα]T (15)
定义悬挂点Oh在机体坐标系上的位置矢量为rh=[0,0,-zh]T。则负载在机体坐标系下的绝对速度vp为:
其中r=rp+rh是负载在机体坐标系下的位置矢量。悬挂负载在机体坐标系下的绝对加速度ap为:
此外,负载在机体坐标系下所受的重力矢量gp表示如下:
其中m0表示负载的质量。忽略负载所受的空气阻力,根据牛顿第二运动定律,负载通过绳子作用在四旋翼无人机悬挂点的拉力Fp为
Fp=-m0ap+gp (19)
负载通过绳子作用在四旋翼无人机悬挂点的力矩Mp为
Mp=rh×Fp (20)
上式中,矢量Fp=[Fpx,Fpy,Fpz]T和Mp=[Mpx,Mpy,Mpz]T都是α和β的函数。负载的运动受到悬挂点处力矩平衡的限制,可得:
-rp×(-m0ap+gp)=0 (21)
将上式展开,可得α和β的表达式,简化后为
其中,c1,c2,c3表达式如下所示:
首先,通过将额外力Fp施加到机体平移的动力学模型(12)中,得到带负载的四旋翼无人机的平移子系统:
然后,将Mp视为负载带来的额外力矩施加到式(13)中,得到带负载的四旋翼无人机的旋转子系统:
最后,结合式(23)和(24),得到带有负载的四旋翼无人机的动力学模型,如下所示:
m1=m0zh(zh+lcosα),m2=m0lzhsinαsinβ,m3=-m0lzh sinαcosβ
c12=wzm0zh(zh+lcosα)+wxm0zhlsinαsinβ-Jyzwy-(Jz-Jy)wz
c22=m0zhwylsinαsinβ+wxzhm0lsinαcosβ+Jyzwx-Jxywz
c31=Jxywx+(Jx-Jy)wy,c32=Jxywy+Jxzwz,c33=Jyzwx
为了方便控制器设计,可以将式(25)分解为以下两个子系统的动力学模型:
平移动力学:
旋转动力学:
其中,和代表模型的耦合项,A=diag{m+m0;m+m0;m+m0}和表示增益矩阵。与式(14)表示的四旋翼无人机模型相比,式(27)表示的旋转子系统的控制增益不再是已知增益矩阵B1。而J对于每一台无人机来说是固定的,可通过查产品手册直接计算得到、或可直接在设计制造无人机时将其转动惯量给定在产品手册中、或者可以根据实验测得无人机的转动惯量;又由于M1与时变的摆动角α和β相关,因此旋转子系统的控制增益矩阵B成为时变的,但只要实时测量得到摆动角α和β,即可计算得到增益矩阵B。选择控制量ut=[Fx/(m+m0),Fy/(m+m0),Fz/(m+m0)]T和ur=[uφ,uθ,uψ]T。da(·)和db(·)代表系统模型的不确定项。
考虑无人机的物理结构和性能指标限制,无人机的速度和角速度都需要满足如下限制:|v|≤kv,|w|≤kw,其中kv和kw是无人机的速度和角速度极所能达到的最大极限值;
考虑输入饱和影响对无人机系统的影响,由于转速限制,ui=[ut,ur](i=1,...,6)不再是设计的控制输入。系统模型(26)和(27)表示如下
此处,p(ui)表示受非对称非光滑的饱和非线性影响的控制输入,定义为
那么,p(ui)可以表示为其中Ψ(ui)是p(ui)和之差。因函数和饱和度函数p(ui)的有界性质,可知函数Ψ(ui)是有界的。即|Ψ(ui)≤Ψm,其中Ψm是正定未知的常数。为方便起见,对函数采用均值定理,变成其中且0<λ<1。通过选择ui0=0并利用事实p(0)=0,得到由此引入以下变量:对于所有是不变的,并且存在正定常数gmax,使得因此,有0<gi≤gmax<∞。最后,非线性系统(2.28)和(2.29)变成
上式中,Lt1(·)=AΨ(ut)+Fp1/(m+m0)+da(·),Gt1=diag{g1,g2,g3},Lr1(·)=BΨ(ur)+M1 -1 Mp1+db(·),Gr1=diag{g4,g5,g6}。
此外,在实际飞行中,四旋翼无人机存在发生执行器故障的情况。只要执行器故障产生,系统控制输入不再是ui而是uai=ρi(t)ui+δi(t)(i=1,...,6),其中ρi(t)∈(0,1]表示执行器效率因子,δi(t)表示由部分控制动作产生的时变且不可测量的矢量函数。当ρi=1,δi=0时,意味着四旋翼无人机的执行器是完美的并且可以正常工作;ρi≠0,δi≠0时,意味着四旋翼无人机发生了部分执行器故障;ρi=0意味着四旋翼无人机中对应的执行器完全无效,即完全失控。
由此,考虑到执行器故障,将式(31)和(32)分别改写为:
平移动力学:
旋转动力学:
上式中,Gt=Gt1ρt,ρt=diag{ρ1,ρ2,ρ3},Lt2(·)=AGt1δt+Lt1(·),Gr=Gr1ρr,ρr=diag{ρ4,ρ5,ρ6},Lr2(·)=BGr1δr+Lr1(·)。
2)设定对无人机的速度跟踪误差,具体包括:
为了便于控制器设计和稳定性分析,定义期望参考为和期望参考v*(t)和w*(t)是已知和有界的,即满足对任意的t≥0,有|v*|≤Av<kv和|w*|≤Aw<kw,其中kv和kw是无人机的速度和角速度极限值,Av=[Av1,Av2,Av3]T和Aw=[Aw1,Aw2,Aw3]T是小于极限值的正常数已知向量。
因此,将平移跟踪误差和旋转跟踪误差分别定义为
ev=v-v* (35)
ew=w-w* (36)
然后,从式(33)到(36)可知,跟踪误差的动力学模型可以写成
3)在无人机上加装传感器测量飞行过程中的摆动角α和β,计算时变增益矩阵B。
4)设计模型参数已知情况下带负载四旋翼无人机的平移子系统和旋转子系统的控制器,并通过设计的控制器控制无人机系统,具体包括:
第一步,考虑式(6)里面的平移子系统,设计基于模型参数已知的受状态约束平移控制器,并实现所需的跟踪控制目标:
为了确保不违反v的约束,即|v|<kv,根据BLF的性质,将李雅普诺夫函数的第一部分定义为
其中kb为正定常数向量,与此同时,定义一个紧凑集合Ωev={ev:|ev|<kb},满足V1在紧凑集合Ωev中是有效的;为了使不等式|v|<kv成立,kb被选定成:
kb=kv-Av (40)
由于ev=v-v*和|v*|≤Av,所以|v|≤|ev|+|v*|<kb+Av=kv-Av+Av=kv;
控制器ut的定义为
其中,σ1>0是控制器设计参数;
第二步,考虑式(6)误差跟踪动态的旋转子系统,设计基于模型参数已知的受状态约束姿态控制器,并实现所需的跟踪控制目标:
为了确保不违反角速度w的约束,即|w|<kw,根据势垒函数BLF的性质,将李雅普诺夫函数的第一部分定义为
其中,kc是控制器设计参数;同时,定义紧凑集合Ωew={ew:|ew|<kc},其满足V2在紧凑集合Ωew中有效,为了满足不等式|w|<kw,选择kc为
kc=kw-Aw (44)
ew=w-w*和|w*|≤Aw,有|w|≤|ew|+|w*|<kc+Aw=kw-Aw+Aw=kw;
控制器uri和更新律定义为:
下面采用仿真实验来验证本实施例中控制方法的有效性。
实验中,给定的期望值轨迹为v*=[0.2*cos(t);0.1*sin(t);0.3*cos(t)],w*=[0.1*cos(t);0.1*sin(t);0.2*cos(t)]。初值条件v(0)=[0.1;0.1;0.1],w(0)=[0.1;0.1;0.05]。此外,选择系统状态约束分别为kv=[1.8;1.8;1.8]和kw=[1.8;1.8;1.8],以及系统状态的期望值约束分别为Av=[0.1;0.2;1.2]和Aw=[0.6;0.6;0.2]。
四旋翼无人机的速度和角速度轨迹跟踪结果,分别如图2-7所示。仿真结果表明,本实施例中设计的控制器在整个跟踪控制过程中都没有违反状态约束的限制,并且具有相当好的跟踪性能,跟踪误差ev和ew始终保持在很小的范围内,明显优于经典的PID控制算法。同时,从局部放大图可以看出,无人机在4秒时出现故障,误差会突然变大,但在本实施例设计的控制器的调节下仍能精确跟踪期望轨迹。
图8显示的是受状态约束下,考虑存在执行器故障与时变转动惯量,所设计的自适应容错控制器的控制输入。图9显示了跟踪过程中吊挂载荷的摆角大小。从中可以看出,所提出的控制策略能使四旋翼跟随参考轨迹移动,并且时间延迟很小,具有较好的飞行稳定性。
最后说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的宗旨和范围,其均应涵盖在本发明的权利要求范围当中。
Claims (1)
1.一种模型参数已知的带负载四旋翼无人机速度控制方法,其特征在于:包括以下步骤:
1)建立带负载四旋翼无人机的动力学模型,具体表示如下:
上式中,fa(·)和fb(·)代表模型的耦合项,A=diag{m+m0;m+m0;m+m0}和表示增益矩阵,m表示无人机的质量,m0表示负载的质量,M1表示惯性矩阵,m1=m0zh(zh+l cosα),m2=m0lzh sinα sinβ,m3=-m0lzh sinαcosβ,其中l是无人机悬挂负载的缆绳的长度,α是缆绳长度l与Zb负方向之间的角度,β是l投影到机体坐标系XbObYb平面的投影线和YbObZb平面之间的角度,zh是负载悬挂点到机体坐标系XbObYb原点的距离,J为在机体坐标系下四旋翼无人机的转动惯量,ut是速度控制量,ur是角速度控制输入;Fp1是负载施加在无人机上的额外力Fp中不含状态量和v的剩余项;Mp1是负载施加于无人机的额外力矩Mp中不含有状态变量和w的剩余项;da(·)和db(·)代表系统模型的不确定项;向量v=[vx,vy,vz]T表示在机体坐标系下的线速度,为v的一阶导数,为带负载四旋翼无人机的平移子系统;向量w=[wx,wy,wz]T表示在机体坐标系下的角速度,是w的一阶导数,为带负载四旋翼无人机的旋转子系统;机体坐标系表示为Ob=(Xb,Yb,Zb),机体坐标系的原点Ob取在四旋翼的质心位置上,Xb轴在四旋翼对称轴内指向机头方向,Zb轴在无人机对称平面内,Zb轴垂直于Xb轴指向上方,按右手定则确定Yb轴;
考虑无人机的物理结构和性能指标限制,无人机的速度和角速度都需要满足如下限制:|v|≤kv,|w|≤kw,其中kv是无人机的速度所能达到的最大极限值,kw是无人机的角速度极所能达到的最大极限值;
考虑输入饱和影响,由于转速限制,ui=[ut,ur],i=1,...,6;不再是设计的控制输入,系统模型表示如下:
此处,p(ui)表示受非对称非光滑的饱和非线性影响的控制输入,定义为
那么,p(ui)表示为其中Ψ(ui)是p(ui)和之差;因函数和饱和度函数p(ui)的有界性质,因此函数Ψ(ui)是有界的,即|Ψ(ui)|≤Ψm,其中Ψm是正定未知的常数;为方便起见,对函数采用均值定理,变成其中且0<λ<1;通过选择ui0=0并利用事实p(0)=0,得到由此引入变量:对于所有 是不变的,并且存在正定常数gmax,使得因此,有0<gi≤gmax<∞;最后非线性系统(2)变成:
考虑执行器故障,此时的控制输入不再是ui而是uai=ρi(t)ui+δi(t)(i=1,...,6),其中ρi(t)∈(0,1]表示执行器效率因子,δi(t)表示由部分控制动作产生的时变且不可测量的矢量函数;将式(4)的模型变成如下等式:
上式中,Gt=Gt1ρt,ρt=diag{ρ1,ρ2,ρ3},Lt2(·)=AGt1δt+Lt1(·),Gr=Gr1ρr,ρr=diag{ρ4,ρ5,ρ6},Lr2(·)=BGr1δr+Lr1(·);
2)设定对无人机的速度跟踪误差,具体包括:
设定平移跟踪误差为ev=v-v*,旋转跟踪误差为ew=w-w*,定义期望参考为和期望参考v*(t)和w*(t)是已知和有界的,即满足对任意的t≥0,有|v*|≤Av<kv和|w*|≤Aw<kw,其中kv和kw是无人机的速度和角速度极限值,Av=[Av1,Av2,Av3]T和Aw=[Aw1,Aw2,Aw3]T是小于极限值的正常数已知向量;
将式(5)转换成跟踪误差的动力学模型如下等式:
3)在无人机上加装传感器测量飞行过程中的摆动角α和β,计算时变增益矩阵B;
4)设计模型参数已知情况下带负载四旋翼无人机的平移子系统和旋转子系统的控制器,并通过设计的控制器控制无人机系统,具体包括:
第一步,考虑式(6)里面的平移子系统,设计基于模型参数已知的受状态约束平移控制器,并实现所需的跟踪控制目标:
为了确保不违反v的约束,即|v|<kv,根据BLF的性质,将李雅普诺夫函数的第一部分定义为
其中kb为正定常数向量,与此同时,定义一个紧凑集合Ωev={ev:|ev|<kb},满足V1在紧凑集合Ωev中是有效的;为了使不等式|v|<kv成立,kb被选定成:
kb=kv-Av (8)
由于ev=v-v*和|v*|≤Av,所以|v|≤|ev|+|v*|<kb+Av=kv-Av+Av=kv;
控制器ut的定义为
其中,σ1>0是控制器设计参数;
第二步,考虑式(6)误差跟踪动态的旋转子系统,设计基于模型参数已知的受状态约束姿态控制器,并实现所需的跟踪控制目标:
为了确保不违反角速度w的约束,即|w|<kw,根据势垒函数BLF的性质,将李雅普诺夫函数的第一部分定义为
其中,kc是控制器设计参数;同时,定义紧凑集合Ωew={ew:|ew|<kc},其满足V2在紧凑集合Ωew中有效,为了满足不等式|w|<kw,选择kc为
kc=kw-Aw (12)
ew=w-w*和|w*|≤Aw,有|w|≤|ew|+|w*|<kc+Aw=kw-Aw+Aw=kw;
控制器uri和更新律定义为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010628316.0A CN111722645B (zh) | 2020-07-01 | 2020-07-01 | 模型参数已知的带负载四旋翼无人机速度控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010628316.0A CN111722645B (zh) | 2020-07-01 | 2020-07-01 | 模型参数已知的带负载四旋翼无人机速度控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111722645A CN111722645A (zh) | 2020-09-29 |
CN111722645B true CN111722645B (zh) | 2022-09-16 |
Family
ID=72571282
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010628316.0A Active CN111722645B (zh) | 2020-07-01 | 2020-07-01 | 模型参数已知的带负载四旋翼无人机速度控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111722645B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112462798B (zh) * | 2020-12-04 | 2021-05-28 | 三生万物(北京)人工智能技术有限公司 | 一种无人机及提高无人机航线飞行性能的方法 |
CN112865897B (zh) * | 2021-01-13 | 2022-08-02 | 山东师范大学 | 一种无人机对地面场景的非平稳信道仿真方法及系统 |
CN113741502B (zh) * | 2021-09-15 | 2023-08-25 | 重庆大学 | 输出约束下带负载四旋翼无人机的神经网络控制方法 |
CN115129089B (zh) * | 2022-08-29 | 2022-12-02 | 国网湖北省电力有限公司技术培训中心 | 无人机拖挂横幅飞行轨迹容错控制方法及设备 |
CN116819976B (zh) * | 2023-08-31 | 2023-11-10 | 中国人民解放军空军工程大学 | 控制输入受约束动力学系统的预定时间容错控制设计方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102193557B (zh) * | 2010-03-04 | 2013-01-02 | 南京航空航天大学 | 无人机的鲁棒受限飞行控制方法 |
US9146557B1 (en) * | 2014-04-23 | 2015-09-29 | King Fahd University Of Petroleum And Minerals | Adaptive control method for unmanned vehicle with slung load |
CN106774373B (zh) * | 2017-01-12 | 2019-11-05 | 哈尔滨工业大学 | 一种四旋翼无人机有限时间姿态跟踪控制方法 |
CN108563124B (zh) * | 2018-04-26 | 2021-01-05 | 宋永端 | 基于api可用户设置轨迹的旋翼无人机自适应控制方法 |
CN108681247B (zh) * | 2018-05-14 | 2021-03-30 | 中国人民解放军火箭军工程大学 | 一种完全分布式保性能通信容错编队控制方法 |
US11693373B2 (en) * | 2018-12-10 | 2023-07-04 | California Institute Of Technology | Systems and methods for robust learning-based control during forward and landing flight under uncertain conditions |
CN109856972B (zh) * | 2019-02-21 | 2020-03-17 | 南京航空航天大学 | 一种无人直升机鲁棒容错跟踪控制方法 |
-
2020
- 2020-07-01 CN CN202010628316.0A patent/CN111722645B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111722645A (zh) | 2020-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111722645B (zh) | 模型参数已知的带负载四旋翼无人机速度控制方法 | |
Bisheban et al. | Geometric adaptive control with neural networks for a quadrotor in wind fields | |
CN111610721B (zh) | 模型参数完全未知的带负载四旋翼无人机速度控制方法 | |
CN110990947B (zh) | 一种火箭助推无人机发射过程多场耦合仿真分析方法 | |
Xu et al. | Sliding mode control of a quadrotor helicopter | |
Das et al. | Dynamic inversion with zero-dynamics stabilisation for quadrotor control | |
Bisgaard et al. | Modeling of generic slung load system | |
Dhadekar et al. | Robust control of quadrotor using uncertainty and disturbance estimation | |
Bonna et al. | Trajectory tracking control of a quadrotor using feedback linearization | |
CN111538255B (zh) | 一种反蜂群无人机的飞行器控制方法及系统 | |
Bulka et al. | Autonomous control of agile fixed-wing UAVs performing aerobatic maneuvers | |
Sridhar et al. | Tilt-rotor quadcopter xplored: Hardware based dynamics, smart sliding mode controller, attitude hold & wind disturbance scenarios | |
CN111459188A (zh) | 一种基于四元数的多旋翼非线性飞行控制方法 | |
Kobilarov | Trajectory tracking of a class of underactuated systems with external disturbances | |
CN113848990B (zh) | 一种考虑结构非线性的无人机稳定飞行控制方法 | |
Ansari et al. | Quadrotor control using generalized dynamic inversion and terminal sliding mode | |
Prach et al. | Development of a state dependent riccati equation based tracking flight controller for an unmanned aircraft | |
Yu et al. | Compensating aerodynamics of over-actuated multi-rotor aerial platform with data-driven iterative learning control | |
Setiawan et al. | Development of real-time flight simulator for quadrotor | |
CN112327629A (zh) | 基于动态补偿的小型无人直升机自适应容错控制方法 | |
Neff et al. | Velocity control for a quad-rotor uav fly-by-camera interface | |
An et al. | Relative position control design of receiver UAV in flying-boom aerial refueling phase | |
CN113238572B (zh) | 基于预设性能控制的预设时间四旋翼无人机姿态跟踪方法 | |
Lee et al. | Design, modeling and control of t 3-multirotor: a tilting thruster type multirotor | |
CN109212969B (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 |