CN115390587A - 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 - Google Patents
一种非预设参数的欠迹向驱动多航天器编队重构控制方法 Download PDFInfo
- Publication number
- CN115390587A CN115390587A CN202211012265.4A CN202211012265A CN115390587A CN 115390587 A CN115390587 A CN 115390587A CN 202211012265 A CN202211012265 A CN 202211012265A CN 115390587 A CN115390587 A CN 115390587A
- Authority
- CN
- China
- Prior art keywords
- sliding mode
- under
- linear transformation
- tracking
- dimension
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 230000015572 biosynthetic process Effects 0.000 title claims abstract description 22
- 230000009466 transformation Effects 0.000 claims abstract description 34
- 239000011159 matrix material Substances 0.000 claims abstract description 21
- 230000014509 gene expression Effects 0.000 claims abstract description 18
- 239000000126 substance Substances 0.000 claims description 9
- 238000009795 derivation Methods 0.000 claims description 7
- 238000013461 design Methods 0.000 claims description 4
- 238000004088 simulation Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000012360 testing method Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- JJWKPURADFRFRB-UHFFFAOYSA-N carbonyl sulfide Chemical compound O=C=S JJWKPURADFRFRB-UHFFFAOYSA-N 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000012546 transfer Methods 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
- G05D1/104—Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying
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
本发明公开了一种非预设参数的欠迹向驱动多航天器编队重构控制方法,所方法包括将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的非预设参数线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的非预设参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
Description
技术领域
本发明涉及欠驱动航天器控制技术领域,具体涉及一种非预设参数的欠迹向驱动多航天器编队重构控制方法。
背景技术
多颗航天器被发射到太空中,首先需要进行构型初始化,接着进行队形的保持,然后根据任务的调整再进行队形重新构造,这些步骤涉及航天器在运行过程中的变轨技术,整体形态发生变化(数量、形状等)。无论是编队飞行、星座构建、集群任务、在轨服务、甚至航天器群博弈,都要求每颗航天器从当前位置转移到指定目标位置,因此重构是多航天器任务执行层面。航天器在运行过程中,由于太空垃圾、宇宙粒子、机械故障等原因可能导致个别推机器失效,此时全驱动控制算法无法应对欠驱动情况,因此研究欠驱动航天器控制算法具有实际意义。
欠驱动情形可分为两种:欠径向驱动和欠迹向驱动,线性化的动力学方程在没有沿轨迹输入的情况下存在一个不可控的特征值,根据当前欠迹向驱动的研究结论指出欠迹向情形不可控状态变量等于圆轨道编队自然构型的保持条件且恒为零,所以不影响系统的能空性。欠迹向驱动系统仅包含5维变量,控制输入只有径向和法向2个自由度,因此设计稳定可控的欠迹向驱动控制方法是一大难点。
进一步的说,当前的面向输入饱和的迹向欠驱动航天器编队重构的控制方法忽视了一个问题,即为了简化控制器的推导,预先设定了欠驱动线性变换矩阵中的某些参数的值。这种操作不仅缺乏理论依据,还可能会限制系统状态之间的耦合关系,预先决定系统的稳定重构精度,从而影响所提出的控制器的普适性。本发明提供的是一种非预设参数的控制器,系统状态之间的耦合关系不像在以前的方法中那样受预设参数的限制,系统收敛性的推导过程及相关表达式更为简洁,同时可以预估制控制器在三个通道位置和速度的控制精度,对实际工程更具借鉴意义。
发明内容
针对现有技术的不足,本发明旨在提供一种非预设参数的欠迹向驱动多航天器编队重构控制方法,通过本发明的方法可以在仅包含径向和法向两个自由度推力时完成卫星编队重构任务。
为了实现上述目的,本发明采用如下技术方案:
一种非预设参数的欠迹向驱动多航天器编队重构控制方法,所方法包括:
S1将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的非预设参数线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;
S2利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;
S3分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
需要说明的是,所述步骤S1包括:
S1.1假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T;
其中,a3,b3,f3均为常数;
S2.1根据滑模面设计规则选择如下滑模面:
其中,c2>0为常数。
S2.2等效控制律ui1可通过对滑模面求导得到参数形式的等效控制律:
其中,
选择开关控制为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2;
需要说明的是,所述步骤S3包括:
S3.2当线性变换矩阵P21中的非预设参数满足
a21=b3/a3>0;
可以的得到相对位置跟踪误差在x,y,z三个通道的收敛误差为:
相对位置和相对运动速度跟踪误差在迹向的收敛误差为:
相对位置和相对运动速度跟踪误差在法向的收敛误差为:
设置滑模控制参数及满足线性变换矩阵P21中的参数值,根据三个通道的相对位置和相对速度收敛误差的表达式,预估控制器(11)式在x,y,z三个通道的相对位置和相对速度的控制精度。
与现有技术相比,本发明提供的是一种非预设参数的控制器,系统状态之间的耦合关系不像以前的方法中那样受预设参数的限制,系统收敛性的推导过程及相关表达式更为简洁,同时可以预估制控制器在三个通道位置和速度的控制精度,对实际工程更具借鉴意义。
附图说明
图1为本发明的基于线性变换的欠迹向驱动多航天器编队重构控制方法的流程图;
图2为本发明仿真试验中跟随者航天器的控制输入示意图,其中,图2a为跟随者航天器1的控制输入、图2b为跟随者航天器2的控制输入,图2c跟随者航天器3的控制输入;
图3为本发明仿真试验中跟随者航天器的相对于领航者航天器的速度示意图,其中,图3a为跟随者航天器1的相对于领航者航天器的速度、图3b为跟随者航天器1的相对于领航者航天器的速度、图3c为跟随者航天器1的相对于领航者航天器的速度;
图4为本发明仿真试验中3颗跟随者航天器在xy平面、xz平面、yz平面和三维空间的轨迹示意图,其中,图4a为3颗跟随者航天器在xy平面的轨迹、图4b为3颗跟随者航天器在xz平面的轨迹、图4c为3颗跟随者航天器在yz平面的轨迹、图4d为3颗跟随者航天器在三维空间轨迹。
具体实施方式
以下将对本发明作进一步的描述,需要说明的是,以下实施例以本技术方案为前提,给出了详细的实施方式和具体的操作过程,但本发明的保护范围并不限于本实施例。
如图1所示,本发明为一种非预设参数的欠迹向驱动多航天器编队重构控制方法,所方法包括:
S1将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;
S2利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;
S3分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
进一步的,本发明的所述步骤S1包括:
S1.1假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T;
其中,a3,b3,f3均为常数;
S2.1根据滑模面设计规则选择如下滑模面:
其中,c2>0为常数。
S2.2等效控制律ui1可通过对滑模面求导得到参数形式的等效控制律:
其中,
选择开关控制为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2;
进一步的,本发明的所述步骤S3包括:
S3.2当线性变换矩阵P21中的参数满足
a21=b3/a3>0;
可以的得到相对位置跟踪误差在x,y,z三个通道的收敛误差为:
相对位置和相对运动速度跟踪误差在迹向的收敛误差为:
|ey|≤Δey=ζin /a 21;
相对位置和相对运动速度跟踪误差在法向的收敛误差为:
设置滑模控制参数及满足线性变换矩阵P21中的参数值,根据三个通道的相对位置和相对速度收敛误差的表达式,预估控制器(11)式在x,y,z三个通道的相对位置和相对速度的控制精度。
实施例1
步骤1.1:假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:其中eiρ=[ex ey ez]T为位置误差,为速度误差。已知欠迹向驱动动力学方程为其中,
将其改写为由一个3维变量和一个2维变量的形式:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T (4)
其中,a3,b3,f3为常数。
其中,c2>0为常数,并且可知:
步骤2.2:对滑模面求导:
则可以得到参数形式的效控制律ui1:
其中,
开关控制选择为ui2=-k1si-k2sigγ(si),其中k1>0,k2>0和0<γ<1为滑模参数,于是欠迹向推力控制器ui可以表达为:
ui=ui1+ui2,(11)
其中,δi=(c2P21+P21A21)di2u+di2a,||δi||≤δi2m=(c2||P12||+||P12A21||+1)dim,而(12)式可以转化为:
上式可以进一步改写为如下两种形式:
其中,ζip=-P21dj2u,且有ζip≤-||P21||dim。进一步,将上式改写为如下两种形式:
对(6)式第二项求导,可得法向的相对速度收敛误差:
根据(6)式第一项可得一个一阶系统:
a21=b3/a3>0,(25)
此时式(24)可以写为如下两种表达式:
于是在径向的相对位置收敛误差为:
径向的相对速度收敛误差为:
步骤3.3:线性变换矩阵P21和滑模控制法的控制器参数设置为:a3=-454.5,b3=-0.4,f3=1,c3=0.003γ=0.5,k1=0.003,k2=0.000001,c2=0.003。假设扰动的边界值为δi2m=1×10-6m/s2,根据式(16)、(20)-(23)、(27)-(28)和(32)-(33)可以预估滑模面、线性化后的相对位置误差、线性化后相对速度误差及径向、迹向和法向的相对位置和相对速度的收敛误差上届依次为: 即在实际的仿真中,重构后期的稳定重构精度小于上述数值。
仿真试验
本发明的实验结果如下:
选择一般圆轨道为当前队形轨道,投影圆轨道为期望构型轨道:
航天器期望的相对运动位置矢量为当前相对位置矢量一般圆轨道和投影圆轨道的半径分别设置为rig=500m,rip=1000m,3颗跟随者航天器的初始队形相角分别为而期望构型的相角为领航者航天器的轨道要素选择为:半长轴6878.137km、偏心率为0、轨道倾角42deg、升交点赤经-60deg、纬度幅角30deg。
图2(a)-(c)为3颗跟随者航天器在径向和法向的控制输入,数值大小为10-3量级。图3(a)-(c)表示3颗跟随者航天器相对于领航者的速度变化情况,相对速度的大小影响着控制器参数的设置,在本实施案例中,在10-3量级控制输入对应于100量级相对运动速度。图4(a)-(c)中三颗跟随者航天器在xy、xz、yz平面上的轨迹;图4(d)为3D重构轨迹,每颗跟随者在图2所示径向和法向两个方向的作用下转移至目标点,最终在投影圆轨道上形成了期望的编队构型,如果控制过大则重构轨迹将变成直线。
对于本领域的技术人员来说,可以根据以上的技术方案和构思,给出各种相应的改变和变形,而所有的这些改变和变形,都应该包括在本发明权利要求的保护范围之内。
Claims (4)
1.一种非预设参数的欠迹向驱动多航天器编队重构控制方法,其特征在于,所方法包括:
S1将欠迹向驱动动力学方程分为3维和2维两个状态变量,设计了一个2×3阶的非预设参数线性变换矩阵将3维变量降为2维,于是欠驱动系统由两个2维变量组成;
S2利用经线性化变换后的2维变量设计滑模面,使用滑模控制法推导出仅包含径向和法向的欠驱动重构控制器;
S3分析当闭环系统收敛至期望构型时线性变换矩阵和滑模控制法需满足的非预设未知参数条件,并推导相对位置和相对速度在径向、迹向和法向的收敛误差表达式,接着设置控制参数,由此预估制控制器在三个通道位置和速度的控制精度。
2.根据权利要求1所述的基于线性变换的欠迹向驱动多航天器编队重构控制方法,其特征在于,所述步骤S1包括:
S1.1假设领航者航天器运行在圆形轨道上,其纬度幅角为n0,定义第i颗跟随者卫星相对领航者卫星的状态误差为:
其中,
ui=[ux uz]T,di2u=[0 diy 0]T,di2a=[dix diz]T;
其中,a3,b3,f3均为常数;
4.根据权利要求1所述的基于线性变换的欠迹向驱动多航天器编队重构控制方法,其特征在于,所述步骤S3包括:
S3.1当k2-δim||si||-γ>0或k1-δim||si||-1>0成立时,闭环系统将在有限时间收敛到期望轨迹,滑模面的收敛误差精度△si可以被描述为:
S3.2当线性变换矩阵P21中的参数满足
a21=b3/a3>0;
可以的得到相对位置跟踪误差在x,y,z三个通道的收敛误差为:
相对位置和相对运动速度跟踪误差在迹向的收敛误差为:
|ey|≤△ey=ζin/a21;
相对位置和相对运动速度跟踪误差在法向的收敛误差为:
设置滑模控制参数及满足线性变换矩阵P21中的参数值,根据三个通道的相对位置和相对速度收敛误差的表达式,预估控制器(11)式在x,y,z三个通道的相对位置和相对速度的控制精度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211012265.4A CN115390587B (zh) | 2022-08-23 | 2022-08-23 | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211012265.4A CN115390587B (zh) | 2022-08-23 | 2022-08-23 | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115390587A true CN115390587A (zh) | 2022-11-25 |
CN115390587B CN115390587B (zh) | 2023-12-15 |
Family
ID=84120418
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211012265.4A Active CN115390587B (zh) | 2022-08-23 | 2022-08-23 | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115390587B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105739511A (zh) * | 2016-01-31 | 2016-07-06 | 中国人民解放军国防科学技术大学 | 一种缺失迹向控制的欠驱动航天器悬停渐近控制方法 |
CN105786013A (zh) * | 2016-03-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 面向输入饱和的迹向欠驱动航天器编队重构控制方法 |
CN106773713A (zh) * | 2017-01-17 | 2017-05-31 | 北京航空航天大学 | 针对欠驱动海洋航行器的高精度非线性路径跟踪控制方法 |
CN114237272A (zh) * | 2021-12-15 | 2022-03-25 | 中国空间技术研究院 | 一种多卫星同步编队重构控制方法 |
-
2022
- 2022-08-23 CN CN202211012265.4A patent/CN115390587B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105739511A (zh) * | 2016-01-31 | 2016-07-06 | 中国人民解放军国防科学技术大学 | 一种缺失迹向控制的欠驱动航天器悬停渐近控制方法 |
CN105786013A (zh) * | 2016-03-29 | 2016-07-20 | 中国人民解放军国防科学技术大学 | 面向输入饱和的迹向欠驱动航天器编队重构控制方法 |
CN106773713A (zh) * | 2017-01-17 | 2017-05-31 | 北京航空航天大学 | 针对欠驱动海洋航行器的高精度非线性路径跟踪控制方法 |
CN114237272A (zh) * | 2021-12-15 | 2022-03-25 | 中国空间技术研究院 | 一种多卫星同步编队重构控制方法 |
Non-Patent Citations (2)
Title |
---|
JIANG SHAO 等: ""Formation Reconfiguration Hovering Technique for Underactuated Satellite"", PROCEEDINGS OF THE 2021 5TH CHINESE CONFERENCE ON SWARM INTELLIGENCE AND COOPERATIVE CONTROL * |
于金龙 等: ""欠驱动航天器编队重构的切换神经网络控制"", 宇航学报 * |
Also Published As
Publication number | Publication date |
---|---|
CN115390587B (zh) | 2023-12-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Horri et al. | Practical implementation of attitude-control algorithms for an underactuated satellite | |
CN109375648B (zh) | 一种多约束条件下椭圆轨道卫星编队构形初始化方法 | |
US8346410B2 (en) | Method for improving maneuverability and controllability by simultaneously applying both reaction wheel-based attitude controller and thruster-based attitude controller | |
CN110989644B (zh) | 一种考虑目标点多终端约束的飞行器轨迹规划方法 | |
CN112550770B (zh) | 一种基于凸优化的火箭软着陆轨迹规划方法 | |
Bevilacqua et al. | Multiple spacecraft rendezvous maneuvers by differential drag and low thrust engines | |
CN113602535B (zh) | 一种微纳卫星在轨自主交会控制的方法及计算机设备 | |
CN112395689A (zh) | 基于凸优化的火箭故障后在线重构方法 | |
Shirazi | Analysis of a hybrid genetic simulated annealing strategy applied in multi-objective optimization of orbital maneuvers | |
CN115309174A (zh) | 一种近地轨道卫星编队联合控制的方法及电子设备 | |
CN109582039A (zh) | 一种采用相对导航信息的j2摄动下最优队形重构方法 | |
CN113343442B (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及系统 | |
CN114237272A (zh) | 一种多卫星同步编队重构控制方法 | |
Liang et al. | Decoupling trajectory tracking for gliding reentry vehicles | |
CN112629339B (zh) | 一种基于直接法的火箭软着陆轨迹规划方法 | |
CN113569391A (zh) | 地月转移轨道的参数确定方法、装置、设备及介质 | |
CN115390587A (zh) | 一种非预设参数的欠迹向驱动多航天器编队重构控制方法 | |
Naasz et al. | Application of several control techniques for the ionospheric observation nanosatellite formation. | |
Yoshimura et al. | Global trajectory design for position and attitude control of an underactuated satellite | |
CN113485396A (zh) | 一种航天器交会对接最后逼近段相对轨道与姿态跟踪控制方法 | |
Benosman et al. | Output trajectory tracking for a switched nonlinear non-minimum phase system: The vstol aircraft | |
CN115390576B (zh) | 一种非预设参数的欠径向推力多卫星编队重构控制方法 | |
Wehbeh et al. | Geometric MPC techniques for reduced attitude control on quadrotors with bidirectional thrust | |
CN111674570A (zh) | 一种航天器三超控制可变包络操纵律设计方法 | |
Shen et al. | Constrained attitude control of agile spacecraft using CMGs |
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 |