CN110262225A - 受约束的空间航天器轨道交会系统的切换控制器设计方法 - Google Patents
受约束的空间航天器轨道交会系统的切换控制器设计方法 Download PDFInfo
- Publication number
- CN110262225A CN110262225A CN201810974364.8A CN201810974364A CN110262225A CN 110262225 A CN110262225 A CN 110262225A CN 201810974364 A CN201810974364 A CN 201810974364A CN 110262225 A CN110262225 A CN 110262225A
- Authority
- CN
- China
- Prior art keywords
- spacecraft
- control system
- sat
- saturation
- switch 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 18
- 230000001133 acceleration Effects 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 4
- BTKMJKKKZATLBU-UHFFFAOYSA-N [2-(1,3-benzothiazol-2-yl)-1,3-benzothiazol-6-yl] dihydrogen phosphate Chemical compound C1=CC=C2SC(C3=NC4=CC=C(C=C4S3)OP(O)(=O)O)=NC2=C1 BTKMJKKKZATLBU-UHFFFAOYSA-N 0.000 claims description 3
- 238000013178 mathematical model Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 3
- 229920006395 saturated elastomer Polymers 0.000 abstract description 3
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
Classifications
-
- 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)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (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
本发明公开了一种受约束的空间航天器轨道交会的切换控制器设计方法,本发明通过空间航天器交会系统的数学建模、非对称饱和函数的转换、切换控制器的驻留时间方法、参量Riccati方程方法等手段,设计了一种基于驻留时间的切换控制器。本发明将具有非对称输入饱和的控制系统等价地表示为具有多个对称输入饱和的控制系统。实际上,得到的等价对称输入饱和系统本质上是一个切换系统。所设计控制器只需要离线求解参量Riccati方程,简单易操作,且计算量低。所提控制方法可以有效提高控制系统状态的收敛速度,满足航天领域对于控制系统快速收敛性能的要求。
Description
技术领域
本发明属于航天领域,具体涉及一种受有限推理约束的空间航天器轨道交会的切换控制器设计方法。
背景技术
空间航天器轨道交会最终要完成目标航天器和追踪航天器之间位置的对准,也就是说要求目标航天器相对于追踪航天器的相对距离为零。空间航天器轨道交会是航天器领域的一项重要的空间技术,是完成一些高级航天任务,如航天器间的拦截、在轨检修、空间垃圾处理等的先决条件。
由于设备的物理限制,执行器出现饱和是不可避免的。执行器饱和会影响控制系统的性能甚至会导致控制系统不稳定。因此,设计控制系统时需要考虑执行器的饱和问题。如空间航天器轨道交会系统,推力器所能提供的推力是有限的。如果根据计算得到的加速度超出了推力器所能供给的最大加速度则执行器发生饱和。这时控制系统不是按照设计的控制器来工作,将会导致空间航天器交会任务的失败。对称执行器饱和是非对称执行器饱和的一个特例,因此研究非对称执行器饱和更具有实际意义。
发明内容
本发明针对现有技术的不足,提出了一种受约束的空间航天器轨道交会的切换控制器设计方法。
本发明目的是设计一种空间航天器交会的切换控制器,来满足控制系统快速收敛性的要求。该设计方法首先通过将具有非对称输入饱和的控制系统等价地表示为具有多个对称输入饱和的控制系统。实际上,得到的等价对称输入饱和系统本质上是一个切换系统。然后基于驻留时间及参量 Riccati方程方法给出了切换控制器的设计。该方法不仅能够解决非对称输入饱和控制问题,同时能够满足控制系统状态快速收敛的要求并能够保证闭环系统的指数稳定性。
本发明的技术方案是通过空间航天器交会系统的数学建模、非对称饱和函数的转换、切换控制器的驻留时间方法、参量Riccati方程方法等手段,设计了一种基于驻留时间的切换控制器,所设计控制器可以有效的提高控制系统状态的收敛速度并保证系统的稳定性。
本发明方法的步骤包括:
步骤(1)建立空间航天器轨道交会系统的数学模型
根据牛顿运动理论及目标轨道坐标系,建立了航天器轨道交会系统的非线性微分方程,如下所示:
其中,x,y,z,分别表示追踪航天器相对于目标航天器在X、Y、Z轴的相对距离和相对速度,aX,aY和aZ分别表示在X、Y、Z轴的加速度,ω为目标航天器的轨道角速度,ω>0,η为引力常数,R为目标轨道半径;其中非线性项
利用泰勒级数展开,将非线性项进行线性化并保留到一阶项,同时选择状态向量得到如下状态空间表达式
其中Sat表示非对称饱和函数,
Sat(u)=[Sat(u1) Sat(u2) Sat(u3)]T,
u=[aX,aY,aZ]T,
并且υi max>0,υi min>0,υi max≠υi min;
重新给出执行器饱和水平
其中,k=4θ1+2θ2+θ3+1,k∈{1,2,3,4,5,6,7,8},并且
υmax=[υ1max υ2max υ3max]T,
υmin=[υ1min υ2min υ3min]T,
计算出的新的执行器饱和水平为:
基于新给出的执行器饱和水平,非对称饱和控制输入转换为8个对称的执行器饱和控制输入,控制系统重新表示为如下的切换系统
其中,
进而,把控制系统表示为具有单位输入饱和的切换系统
其中δ(t)=4θ1+2θ2+θ3+1,
步骤(2)设计切换控制器
(a)基于上述模型设计椭球集合
(b)P(ιj)通过求解下面的参量Riccati方程得到
P(ιj)A+ATP(ιj)-P(ιj)BBTP(ιj)=-ιjP(ιj)
(c)确定切换控制器驻留时间
切换控制器的最大驻留时间为
其中,符号λmax表示矩阵的最大特征值,R6表示欧几里得空间;
(d)构造切换控制器u
当系统的状态属于集合Ψ时,(Ψ=ESα+1-ESα),
u=-BTP(ια)X,α=0,1,2,...,N-1.
当系统的状态属于集合ESN时,
u=-BTP(ιN)X。
本发明与现有技术相比具有以下效果:
本发明将具有非对称执行器饱和的航天器交会系统控制问题转换为具有单位对称执行器饱和的切换系统进行研究。所设计控制器只需要离线求解参量Riccati方程,简单易操作,且计算量低。所提控制方法可以有效提高控制系统状态的收敛速度,满足航天领域对于控制系统快速收敛性能的要求。
具体实施方式
受约束的空间航天器轨道交会的切换控制器设计方法,该控制器的设计具体包括以下步骤:
a)建立了航天器轨道交会系统的非线性微分方程,如下所示:
其中,x,y,z,分别表示追踪航天器相对于目标航天器在X、Y、Z轴的相对距离和相对速度,aX,aY和aZ分别表示在X、Y、Z轴的加速度,ω>0 (η为引力常数,R为目标轨道半径)为目标航天器的轨道角速度。其中非线性项
(b)利用泰勒级数展开,将非线性项进行线性化并保留到一阶项,同时选择状态向量可以得到如下状态空间表达式
其中Sat表示非对称饱和函数,
Sat(u)=[Sat(u1) Sat(u2) Sat(u3)]T,
u=[aX,aY,aZ]T,
并且υi max>0,υi min>0,υi max≠υi min.
(c)根据下式计算新的执行器饱和水平
其中,k=4θ1+2θ2+θ3+1,并且
υmax=[υ1max υ2max υ3max]T,
υmin=[υ1min υ2min υ3min]T,
根据θi的不同取值,k可取1,2,3,4,5,6,7,8。
计算所得新的执行器饱和水平为:
基于上面计算出的新的执行器饱和水平,非对称饱和控制输入可以转换为8个对称的执行器饱和控制输入,控制系统可以重新表示为
其中,
进而,我们可以把控制系统表示为具有单位输入饱和的切换系统
其中δ(t)=4θ1+2θ2+θ3+1,
(d)根据上面计算出的新的执行器饱和水平,可以进而计算的值。
(e)选取满足条件0<ι0<ι1<ι2<…<ιN的参数ιj,j=0,1,2,...,N。并根据上述计算出的的值来构造如下的椭球集合:
(f)求解下面的参量Riccati方程得到矩阵P(ιj)
P(ιj)A+ATP(ιj)-P(ιj)BBTP(ιj)=-ιjP(ιj)
(g)计算切换控制器驻留时间
切换控制器的最大驻留时间为
其中,符号λmax表示矩阵的最大特征值。
(h)基于以上得到的参数,根据下面的公式计算所设计的切换控制器当系统的状态属于集合Ψ(Ψ=ESα+1-ESα)时,
u=-BTP(ια)X,α=0,1,2,...,N-1.
当系统的状态属于集合ESN时,
u=-BTP(ιN)X。
Claims (1)
1.受约束的空间航天器轨道交会系统的切换控制器设计方法,其特征在于,该控制器设计方法包括以下步骤:
步骤(1)建立空间航天器轨道交会系统的数学模型
根据牛顿运动理论及目标轨道坐标系,建立了航天器轨道交会系统的非线性微分方程,如下所示:
其中,x,y,z,分别表示追踪航天器相对于目标航天器在X、Y、Z轴的相对距离和相对速度,aX,aY和aZ分别表示在X、Y、Z轴的加速度,ω为目标航天器的轨道角速度,ω>0,η为引力常数,R为目标轨道半径;其中非线性项
利用泰勒级数展开,将非线性项进行线性化并保留到一阶项,同时选择状态向量得到如下状态空间表达式
其中Sat表示非对称饱和函数,
Sat(u)=[Sat(u1) Sat(u2) Sat(u3)]T,
u=[aX,aY,aZ]T,
并且υimax>0,υimin>0,υimax≠υimin;
重新给出执行器饱和水平
其中,k=4θ1+2θ2+θ3+1,k∈{1,2,3,4,5,6,7,8},并且
υmax=[υ1max υ2max υ3max]T,
υmin=[υ1min υ2min υ3min]T,
计算出的新的执行器饱和水平为:
基于新给出的执行器饱和水平,非对称饱和控制输入转换为8个对称的执行器饱和控制输入,控制系统重新表示为如下的切换系统
其中,
进而,把控制系统表示为具有单位输入饱和的切换系统
其中δ(t)=4θ1+2θ2+θ3+1,
步骤(2)设计切换控制器
(a)基于上述模型设计椭球集合
(b)P(ιj)通过求解下面的参量Riccati方程得到
P(ιj)A+ATP(ιj)-P(ιj)BBTP(ιj)=-ιjP(ιj)
(c)确定切换控制器驻留时间
切换控制器的最大驻留时间为
其中,符号λmax表示矩阵的最大特征值,R6表示欧几里得空间;
(d)构造切换控制器u
当系统的状态属于集合Ψ时,(Ψ=ESα+1-ESα),
u=-BTP(ια)X,α=0,1,2,...,N-1.
当系统的状态属于集合ESN时,
u=-BTP(ιN)X。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810974364.8A CN110262225B (zh) | 2018-08-24 | 2018-08-24 | 受约束的空间航天器轨道交会系统的切换控制器设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810974364.8A CN110262225B (zh) | 2018-08-24 | 2018-08-24 | 受约束的空间航天器轨道交会系统的切换控制器设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110262225A true CN110262225A (zh) | 2019-09-20 |
CN110262225B CN110262225B (zh) | 2022-01-25 |
Family
ID=67911604
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810974364.8A Active CN110262225B (zh) | 2018-08-24 | 2018-08-24 | 受约束的空间航天器轨道交会系统的切换控制器设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110262225B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112000132A (zh) * | 2020-07-08 | 2020-11-27 | 中国人民解放军军事科学院国防科技创新研究院 | 基于椭球体描述的航天器避障控制方法 |
CN112346472A (zh) * | 2020-11-20 | 2021-02-09 | 杭州电子科技大学 | 航天器交会系统的保成本鲁棒增益调度控制器设计方法 |
CN113110548A (zh) * | 2021-04-21 | 2021-07-13 | 北京控制工程研究所 | 一种航天器椭球集合演化的设计方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103728980A (zh) * | 2014-01-08 | 2014-04-16 | 哈尔滨工业大学 | 航天器相对轨道的控制方法 |
CN104076818A (zh) * | 2014-07-02 | 2014-10-01 | 哈尔滨工业大学 | 考虑线性化误差的空间交会系统的增益调度控制方法 |
CN104656447A (zh) * | 2015-01-16 | 2015-05-27 | 西北工业大学 | 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法 |
CN105159304A (zh) * | 2015-06-26 | 2015-12-16 | 哈尔滨工业大学 | 接近并跟踪空间非合作目标的有限时间容错控制方法 |
CN105843077A (zh) * | 2016-04-07 | 2016-08-10 | 杭州电子科技大学 | 一种航天器交会系统的增益切换方法及最大吸引域估计 |
CN106275509A (zh) * | 2016-08-19 | 2017-01-04 | 浙江大学 | 一种空间悬浮绳系组合体面内摆动控制方法 |
CN106407619A (zh) * | 2016-11-16 | 2017-02-15 | 哈尔滨工业大学 | 一种控制受限航天器交会控制系统的线性反馈全局镇定方法 |
CN106557659A (zh) * | 2016-11-21 | 2017-04-05 | 中国人民解放军装备学院 | 考虑j2摄动的大规模空间目标交会对象及次序搜索方法 |
CN106773679A (zh) * | 2016-12-01 | 2017-05-31 | 北京航空航天大学 | 一种基于角速度观测器的航天器容错控制方法 |
CN106970530A (zh) * | 2017-04-28 | 2017-07-21 | 西北工业大学 | 空间非合作目标自主视线交会的无模型预设性能控制方法 |
-
2018
- 2018-08-24 CN CN201810974364.8A patent/CN110262225B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103728980A (zh) * | 2014-01-08 | 2014-04-16 | 哈尔滨工业大学 | 航天器相对轨道的控制方法 |
CN104076818A (zh) * | 2014-07-02 | 2014-10-01 | 哈尔滨工业大学 | 考虑线性化误差的空间交会系统的增益调度控制方法 |
CN104656447A (zh) * | 2015-01-16 | 2015-05-27 | 西北工业大学 | 一种航天器抗干扰姿态跟踪的微分几何非线性控制方法 |
CN105159304A (zh) * | 2015-06-26 | 2015-12-16 | 哈尔滨工业大学 | 接近并跟踪空间非合作目标的有限时间容错控制方法 |
CN105843077A (zh) * | 2016-04-07 | 2016-08-10 | 杭州电子科技大学 | 一种航天器交会系统的增益切换方法及最大吸引域估计 |
CN106275509A (zh) * | 2016-08-19 | 2017-01-04 | 浙江大学 | 一种空间悬浮绳系组合体面内摆动控制方法 |
CN106407619A (zh) * | 2016-11-16 | 2017-02-15 | 哈尔滨工业大学 | 一种控制受限航天器交会控制系统的线性反馈全局镇定方法 |
CN106557659A (zh) * | 2016-11-21 | 2017-04-05 | 中国人民解放军装备学院 | 考虑j2摄动的大规模空间目标交会对象及次序搜索方法 |
CN106773679A (zh) * | 2016-12-01 | 2017-05-31 | 北京航空航天大学 | 一种基于角速度观测器的航天器容错控制方法 |
CN106970530A (zh) * | 2017-04-28 | 2017-07-21 | 西北工业大学 | 空间非合作目标自主视线交会的无模型预设性能控制方法 |
Non-Patent Citations (3)
Title |
---|
QIAN WANG等: "Robust Control for Spacecraft Rendezvous System with Actuator Unsymmetrical Saturation: A Gain Scheduling Approach", 《INTERNATIONAL JOURNAL OF CONTROL》 * |
WANG QIAN等: "Reliable Gain Scheduling Output Tracking Control for Spacecraft Rendezvous", 《INTERNATIONAL JOURNAL OF CONTROL》 * |
王茜等: "输入饱和系统的离散增益调度控制及其在在轨交会中的应用", 《自动化学报》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112000132A (zh) * | 2020-07-08 | 2020-11-27 | 中国人民解放军军事科学院国防科技创新研究院 | 基于椭球体描述的航天器避障控制方法 |
CN112346472A (zh) * | 2020-11-20 | 2021-02-09 | 杭州电子科技大学 | 航天器交会系统的保成本鲁棒增益调度控制器设计方法 |
CN112346472B (zh) * | 2020-11-20 | 2022-05-20 | 杭州电子科技大学 | 航天器交会系统的保成本鲁棒增益调度控制器设计方法 |
CN113110548A (zh) * | 2021-04-21 | 2021-07-13 | 北京控制工程研究所 | 一种航天器椭球集合演化的设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110262225B (zh) | 2022-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110262225B (zh) | 受约束的空间航天器轨道交会系统的切换控制器设计方法 | |
Zhang et al. | Integrated translational and rotational finite-time maneuver of a rigid spacecraft with actuator misalignment | |
CN112445234B (zh) | 一种航天器的姿态控制方法和装置 | |
Li et al. | Discrete-time pure-tension sliding mode predictive control for the deployment of space tethered satellite with input saturation | |
CN109164822B (zh) | 一种基于具有混合执行机构的航天器姿态控制方法 | |
Jiang et al. | Novel integral sliding mode control for small-scale unmanned helicopters | |
CN105843077B (zh) | 一种航天器交会系统的增益切换方法及最大吸引域估计 | |
Yamasaki et al. | Robust trajectory-tracking method for UAV guidance using proportional navigation | |
Guo et al. | Task space control of free-floating space robots using constrained adaptive RBF-NTSM | |
Li et al. | Design of asymptotic second-order sliding mode control for satellite formation flying | |
CN109625332B (zh) | 一种平动点轨道交会无需初始误差符号的预设性能控制方法 | |
Raj et al. | Robust higher-order sliding mode control systems for roll-coupled maneuvers of aircraft using output feedback | |
Yao et al. | Trajectory tracking control of a stratospheric airship with fuzzy sliding mode control | |
CN112835372A (zh) | 四旋翼无人机的固定时间控制方法 | |
Lu et al. | Prescribed-time attitude control during solid thrusting orbital maneuver | |
Xu et al. | Large angle attitude control of cubesat based on multi-objective optimization | |
Yin et al. | Attitude stability control of capturing spacecraft without angular velocity | |
Zhang et al. | Adaptive spacecraft relative position tracking control for non-cooperative target with thrust saturation | |
Liu et al. | Finite-time adaptive terminal sliding mode controller design for reusable launch vehicle in reentry phase | |
Esfahanian et al. | Time-varying dynamic sliding mode control of a nonlinear aerospace launch vehicle | |
Wang et al. | Finite-time adaptive sliding mode control of reentry Reusable Launch Vehicle | |
Lin et al. | Backstepping sliding mode control for flight attitude of quadrotor UAVs | |
Kemao et al. | Improved backstepping control for large angle maneuvers of spacecraft | |
Zhang et al. | Adaptive tracking control for spacecraft hovering to noncooperative target with trajectory generator | |
Wang et al. | A novel global sliding mode control strategy for attitude control of reusable launch vehicles in the reentry phase |
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 |