基于改进人工势函数的航天器近距离安全控制方法1)

doi:10.6052/0459-1879-20-112

基于改进人工势函数的航天器近距离安全控制方法¹⁾

许丹丹^,^*^,²⁾, 张进^,^†^,³⁾

^*空军工程大学航空机务士官学校航空机械工程系, 河南信阳 464000

^†国防科技大学空天科学学院, 长沙 410073

A COLLISION-AVOIDANCE CONTROL ALGORITHM FOR SPACECRAFT PROXIMITY OPERATIONS BASED ON IMPROVED ARTIFICIAL POTENTIAL FUNCTION ¹⁾

Xu Dandan^,^*^,²⁾, Zhang Jin^,^†^,³⁾

^*Department of Aeronautical Mechanical Engineering, Air Cadets School, Air Force Engineering University, Xinyang 464000, Henan, China

^†College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China

通讯作者: 2) 许丹丹, 助教, 主要研究方向: 航天动力学与控制. E-mail:1106124587@qq.com;3) 张进, 副教授, 主要研究方向: 航天飞行任务规划. E-mail:zhangjin@nudt.ed.cn

收稿日期: 2020-04-14 接受日期: 2020-06-28 网络出版日期: 2020-11-18

基金资助:

1) 国家自然科学基金资助项目. 11572345

Received: 2020-04-14 Accepted: 2020-06-28 Online: 2020-11-18

作者简介 About authors

摘要

针对航天器近距离操作的安全问题,提出了一种基于人工势函数改进的碰撞规避控制算法.根据航天器与目标、障碍物之间的实时状态, 利用人工势函数算法,计算航天器的实时加速度, 规划航天器的轨迹. 为改进人工势函数方法的适用性,提出三个方面的改进措施: 首先, 在人工势函数算法中, 为提高碰撞预警的准确性,减少额外机动, 碰撞预警采用碰撞概率代替相对距离. 其次, 为了提高对接安全,降低接近目标航天器的相对速度, 利用相对速度的安全接近走廊来计算目标排斥力.最后, 针对大多航天器不能提供任意连续变化推力的情况, 设置两种实用的推力形式,如bang-bang控制的推力形式和恒定变化率的推力形式, 代替连续变推力形式.通过对不同算例的比较,成功地揭示了主要任务参数(如碰撞预警方法、速度安全边界和实际加速度形式)对近距离操作安全的影响.结果表明, 该方法可以提高航天器近距离操作的安全性、效率性, 并且结构简单,实时性强.

关键词： 近距离操作 ; 碰撞规避 ; 人工势函数 ; 瞬时碰撞概率 ; 安全接近走廊

Abstract

Aiming at the proximity operation of spacecraft, a collision-avoidance control algorithm based on improved Artificial Potential Fields (APF) method is proposed. According to the real-time state between the spacecraft and the target and obstacles, the APF method is used to calculate the real-time acceleration of the spacecraft, and the trajectory of the spacecraft is planned. In order to improve the applicability of the artificial potential function method, three improvement measures are proposed. First, in order to improve the accuracy of collision warning and reduce additional maneuvers, the collision probability combined with the relative distance, instead of only the relative distance, is used to evaluate the collision. Second, in order to increase the docking safety and slow down the approaching relative velocity, the safety boundary and control margin of relative velocity are used to calculate the target repulsion force. Third, in view of the fact that most spacecraft cannot provide any continuously varying thrust, two practical thrust forms, including the thrust with upper limit and constant variation rate and the bang-bang thrust, are used to substitute for the continuously variable thrust form. Numerical simulations are executed to validate the proposed method. The effects of the major mission parameters, such as the collision warning method, the target repulsion acceleration and acceleration forms, are successfully revealed by the comparison between different examples. The results show that the proposed method can improve the safety and efficiency of the spacecraft proximity operations, and has simple structure and strong real-time performance.

Keywords： proximity operation ; collision avoidance ; artificial potential field ; instantaneous collision probability ; safely approaching Corridor

PDF (3952KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

许丹丹, 张进. 基于改进人工势函数的航天器近距离安全控制方法¹⁾. 力学学报[J], 2020, 52(6): 1581-1589 DOI:10.6052/0459-1879-20-112

Xu Dandan, Zhang Jin. A COLLISION-AVOIDANCE CONTROL ALGORITHM FOR SPACECRAFT PROXIMITY OPERATIONS BASED ON IMPROVED ARTIFICIAL POTENTIAL FUNCTION ¹⁾. Chinese Journal of Theoretical and Applied Mechanics[J], 2020, 52(6): 1581-1589 DOI:10.6052/0459-1879-20-112

引言

随着空间技术的进步, 世界主要大国或组织大力发展空间在轨近距离操作服务技术.对自主交会对接、空间目标监视、推进剂补给等单航天器操作技术进行了较为充分的验证^[1-3].大部分接触式在轨服务需要靠近目标进行操作, 而在近距离空间任务中,航天器不与其他空间目标相撞是轨道飞行的基本要求^[4-5]. 因此,航天器近距离接近技术首先需要解决安全逼近控制问题.

人工势场法(artificial potential fields, APF)是最常用的考虑碰撞约束的航天器安全制导与控制方法. APF的优势是:能够引导飞行器进行避撞机动, 并且能够利用势场函数的负梯度方向来设计避撞控制律. APF方法结构简单、计算量小,在实时避撞和平滑的轨迹控制方面得到广泛的应用. 在1986年, APF方法最早由Khatib在机械臂避撞任务中提出^[6], 1993年,McInnes将APF引入到航天器的自主接近机动控制^[7]. 目前, APF研究已经扩展到考虑分布式控制^[8]、具有固定障碍物回避的自主交会^[9]轨道组装的自主控制^[10]和星群编队控制^[11].McCamish等^[12]、Nastasi等^[13]提出结合线性二次型调节器(linear quadratic regulator, LQR)和APF用于解决圆轨道目标最优接近控制问题.在此基础上, Bevilacqua等提出沿墙跟踪方法^[14], Palacios等实现椭圆轨道目标的跟踪控制^[15], Wang等^[16]提出基于碰撞概率的目标交会控制.美国的Prox-1任务首次成功完成基于APF的自主接近操作^[17-18]. 此外,朱安等^[19]、艾海平等^[20]还提出了空间机器人双臂捕获卫星操作过程避免关节冲击破坏的避撞柔顺控制问题.

然而, 将APF控制算法直接用于航天器近距离操作任务, 还存在以下不足. 首先,基于最小距离进行碰撞预警的存在虚警率较高的问题^[21-22]. 其次,防止服务航天器与目标航天器的接近碰撞是最重要的,使用APF控制方法来控制接近目标的速度可能违背目标附近接近速度安全走廊要求.最后, 因为大多数航天器都只安装一组不同推力大小的常推力发动机,不能提供任意连续变化的推力,所以基于APF控制方法计算得到的任意连续变化加速形式的控制律在大多数实际任务中很难实现.

针对上述问题, 为了提高机动效率, 增大机动安全性,本文以航天器交会对接的近距离操作任务为背景,设计一种基于改进APF的防撞控制方法, 并利用数据仿真来验证改进算法的特性.

1 基于人工势函数的安全逼近

如图1所示, 在目标航天器轨道坐标系(LVLH)中, 服务航天器进行交会、逼近机动时, 不仅需要确保自身的轨迹安全, 即不与空间内的障碍物相碰撞, 还要保证与目标航天器安全接触. 因此, 当航天器进行轨迹规划时, 一定要考虑碰撞规避问题.

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 相对运动框架

Fig. 1 Relative motion Frame

采用人工势函数控制航天器时, 控制加速度表示为所有引力和斥力的和. 本文研究探讨了与位置误差、速度误差相关的人工势函数^[23]用于控制航天器.服务航天器的控制加速度表示为所有引力和斥力的和

(1)$a=a_{g} +a_{o} +a_{m}$

其中, $a_{g} $为目标航天器的引力加速度, $a_{o} $为障碍物的斥力加速度, $a_{m} $为目标航天器的斥力加速度.

1.1 引力加速度

引力加速度主要是指由终端位置引力场所产生的加速度, 导致服务航天器从初始位置运动到终端位置. 常见的引力场是二次型引力场,引力场函数$V_{g} $选为

(2)$V_{g} =\frac{1}{2}r_{cg}^{T} Pr_{cg}$

式中, $r_{cg} =r_{g} -r_{c} $为服务航天器到终端目标的位置矢量, $r_{g} $为终端位置矢量, $r_{c} $为初始位置, 取$P=1/\|r_{cg}\|$.

由于终端位置的吸引力, 服务航天器所需的速度$v_{g} $为

(3)$v_{g} =-k_{g} \cdot (r_{cg} /\| r_{cg}\|)$

其中, 吸引力速度影响函数$k_{g} $为

(4)$k_{g} =k\cdot v_{\max} \cdot (1-{e}^{-t})\cdot \left[ {1-{e}^{-\left( {b_{g} \cdot V_{g} } \right)}} \right]$

$k_{g} $的选择决定了控制算法的收敛性. 式(4)中, $(1-{e}^{-t})$可以避免在初始时刻,服务航天器的控制加速度过大, $k$和$b_{g} $为常数, $v_{\max}$为估计的最大相对运动速度.

由终端引力场所产生的航天器加速度$a_{g} $表示为

(5)$a_{g} =\Delta v/\Delta t=(v_{g} -v)/\Delta t$

其中, $v$为服务航天器的相对速度矢量, $\Delta t$为服务航天器的制导控制计算的时间步长.

1.2 斥力加速度

斥力加速度主要是由远离障碍物的斥力场提供的, 斥力场是距离的单调减函数: 航天器距离障碍物越近, 斥力势函数的值就越大. 障碍势能选高斯函数形式

(6)$V_{o} =\left({e}^{-\left( {\left\| {r_{co} } \right\|^2/(2\sigma ^2)} \right)}-{e}^{-\left( {D^2/(2\sigma ^2)} \right)} \right)\cdot\\ \left( {e}^{-\left( {L_{o}^2 /(2\sigma ^2)} \right)}-{e}^{-\left( {D^2/(2\sigma ^2)} \right)} \right)^{-1}$

其中, $r_{co}=r_{c} -r_{o} $表示服务航天器与障碍物的相对位置矢量; $r_{c} $为服务航天器的位置矢量; $r_{o} $为障碍物的位置矢量; 标准偏差$\sigma=(L_{o} +D)/3$, 当航天器相对速度接近零时障碍物表面略大于一个标准偏差, 围绕障碍物的安全区域更合理.

障碍物控制区域的球半径$D$为

(7)$D=d_0 (L_0 +D_{s} )$

式中, $d_0 $为常数, $L_0 $为障碍物与服务航天器的半径和, $D_{s} =\| v \|^2/(2a_{\max} )$是最小制动距离, $a_{\max}$为最大加速度数值.

当服务航天器与障碍物的相对位置$\| {r_{co} } \|\leqslant D$时, 由$n$个障碍物的斥力场所产生的服务航天器加速度$a_{o} $可表示为

(8)$a_{o} =(\sum_{obs=0}^n {v_{co} } -v)\Bigg/\Delta t$

式中, 由于避开障碍物所需的速度$v_{co} $为

(9)$v_{co} =k_{i} \cdot k_{g} \cdot V_{o} \cdot (r_{co} /\left\| {r_{co} } \right\|)=k_{o} \cdot (r_{co} /\left\| {r_{co} } \right\|)$

式中, $k_{i} $和$k_{o} $为常数. 斥力速度影响参数$k_{o} $的选择与吸引力速度影响参数$k_{g} $有关, 以确保航天器的速度能达到所需的临界阻尼.

1.3 目标斥力加速度

服务航天器进行交会、逼近机动时, 终端位置为目标航天器, 因此, 航天器在防止与其他物体碰撞的同时, 需要安全逼近到目标表面.尤其当服务航天器与目标航天器非常接近时, 应对防撞控制机动的范围需要进行一定限制, 因此, 安全对接距离约束表示为

(10)$\left\| {r_{\min} } \right\|\geqslant L$

其中, $r_{\min}$表示两航天器的最近距离, $L$表示服务航天器与目标航天器的半径和.

当相对距离$\| r_{cg}\|\leqslant D_{m} $ $(D_{m} =d_{m} +D_{s})$, $d_{m}$为常数)时, 使用对接安全参数$k_{s}$来修改服务航天器在接近目标时所需的排斥速度, 表达式为

(11)$v_{m} =k_{s} \cdot k_{o} \cdot (r_{cg} \Big/\left\| {r_{cg} } \right\|)$

其中, $k_{s} $通常大于或等于0. 如果$k_{s} =0$, 则$v_{m} $不受目标位置的影响.

由目标航天器的斥力场所产生的航天器加速度, 表示为

(12)$a_{m} =(v_{m} -v)/\Delta t$

2 改进的人工势函数

针对APF控制方法做出三点改进: 第一点是采用瞬时碰撞概率法代替最小距离法进行碰撞预警.第二点是基于速度安全接近走廊设计目标斥力来减慢服务航天器的接近速度.第三点是改进执行加速度实现形式, 设计了两种实现形式的控制执行加速度.

2.1 基于瞬时碰撞概率的人工势函数

为了提高决策准确度并减少机动次数, 设计一种基于瞬时碰撞概率的APF控制方法. 引力和斥力函数的形式与影响因素保持不变, 只改变碰撞预警方式, 即当航天器与障碍物的瞬时碰撞概率$P_{c} \leqslant P$ $(P$为触发碰撞预警的门限值)时,考虑由障碍物的斥力场所产生的服务航天器加速度.

将服务航天器简化为质点, 其位置分布服从相对轨迹偏差传播模型, 瞬时碰撞概率$P_{c}$^[24]的计算公式为

(13)$P_{c} =\iiint_{\varOmega _{c} } {\rho \mbox{d}x\mbox{d}y\mbox{d}z}=\frac{1}{(2\pi )^{3/2}\left\| C \right\|^{1/2}}\cdot\\ \iiint_{\varOmega _{c} } {\exp \left( {-\frac{1}{2}r_{co}^{T} C^{-1}r_{co} } \right)\mbox{d}x\mbox{d}y\mbox{d}z}$

式中, $\rho$为概率密度函数, $C$表示相对位置分布的协方差矩阵, 积分域$\varOmega _{c} $为障碍物的影响区域.

2.2 基于速度走廊设计目标航天器斥力

服务航天器执行终端位置为目标航天器的机动时, 无论过程中出现任何故障或偏差, 都必须要保证服务航天器与目标航天器不发生碰撞.以往应对目标防撞控制机动的范围限制是根据经验的最小距离约束, 没有考虑速度安全接近走廊约束. 因此, 本文考虑根据速度安全接近走廊设计目标航天器斥力.

速度安全接近走廊可表征为

(14)$v_d \leqslant \left\| v \right\|\leqslant \omega \cdot \left\| {r_{cg} } \right\|$

式中, $v_d $是一个常数, 本节取$v_d =5.0\times 10^{-4}$ m/s. 如果与标准状态向量的偏差达到了预先计算的安全界限, 也就是说$\left\| v \right\|>\omega \cdot \left\| {r_{cg} } \right\|$, 这个轨迹就会被判断为"不安全'', 从而发出防撞控制机动指令.

在此, 采用接近安全参数$k_{s} $来修改服务航天器接近目标时所需的排斥速度, 表达式为

(15)$v_{m} =k_{s} \cdot k_{o} \cdot (r_{cg} /\left\| {r_{cg} } \right\|)$

本文取$k_{s} =0.8[1-{e}^{-(V_{g} -L_{o} /2)}]$.

由目标航天器的斥力场所产生的航天器加速度形式仍如式(12).

2.3 设计控制加速度变化形式

本文考虑实际发动机运动机理, 将实际推力形式与人工势函数结合, 控制服务航天器安全机动. 实际飞行中的推力形式通常为以下两种形式, 具有上限和恒定变化率的加速度, 以及设置多个档位开关的bang-bang加速度. 假设在近距离操作的研究中忽略质量变化对控制的影响,则航天器被控加速度的设计直接对应于推进器的设计. 发动机在3个主轴的运动方向上均有一个最大加速度限制$a_{\max}$.以下加速度设计表示的是一个主轴方向上的加速度$a_{i}\ (i=x,y,z)$.

2.3.1 加速速率设计

假设航天器由安装在某方向的1台发动机提供推力, 最大加速度上限为$a_{\max } $, 在到达最大加速度前的加速度变化率为$\dot{{a}}$; 发动机关机时, 加速度由当前值逐渐减小到零的变化率也为$\dot{{a}}$, 则发动机工作过程中某时刻$t$的加速度满足均满足式(16)

(16)$a_{i(t)}= \begin{cases} a_{\max},& a_{i(t)}\geqslant a_{\max}\\ a_{i(t-1)}+\dot{a}\Delta t, & a_{i(t)}> a_{i(t-1)} \ \text{and}\ |a_{i(t)}|<a_{\max}\\ a_{i(t-1)}, & a_{i(t)} = a_{i(t-1)} \ \text{and}\ |a_{i(t)}|<a_{\max}\\ a_{i(t-1)}-\dot{a}\Delta t, & a_{i(t)}< a_{i(t-1)} \ \text{and}\ |a_{i(t)}|<a_{\max}\\ -a_{\max},&a_{i(t)}\leqslant a_{\max}\\ \end{cases}$

其中, $a_{t-1} $是前一时刻的加速度. 假设发动机推力由零增加到最大或由最大减小到零需要一定时间为d$t$, 加速速率$\dot{{a}}=a_{\max } /{d}t$.

2.3.2 bang-bang控制多档位设计

假设航天器某方向上的推力由该方向上安装的4个推力相同为$a_0 =a_{\max } /4$的发动机共同提供, 设置开关逻辑控制发动机工作,开关分档如式(17)所示

(17)$a_{i} ={\begin{cases} a_{\max },& a_{i} \geqslant a_{\max } \\ 0.75a_{\max },& 0.75a_{\max } <a_{i} \leqslant a_{\max } \\ 0.5a_{\max }, &0.5a_{\max } <a_{i} \leqslant 0.75a_{\max } \\ 0.25a_{\max }, &0.25a_{\max } <a_{i} \leqslant 0.5a_{\max } \\ 0,&-0.25a_{\max } <a_{i} \leqslant 0.25a_{\max } \\ -0.25a_{\max },& -0.5a_{\max } <a_{i} \leqslant -0.25a_{\max } \\ -0.5a_{\max },& -0.75a_{\max } <a_{i} \leqslant -0.5a_{\max } \\ -0.75a_{\max },& -a_{\max } <a_{i} \leqslant -0.75a_{\max } \\ -a_{\max },& a_{i} \leqslant -a_{\max } \\ \end{cases}}$

3 仿真算例

近距离机动中, 时间和推进剂消耗是评估控制算法性能的主要指标. 机动时间$t_{d}$单位为s, 推进剂消耗质量$\Delta m$单位为kg. 假设障碍物为固定障碍物, 其运动状态相对目标不发生变化, 在目标航天器质心相对坐标下, 服务航天器以及固定障碍物的仿真初始参数配置如表1所示.

表1 基于最小距离的APF仿真初始参数

Table 1 Initial parameters of APF simulation based on minimum distance

新窗口打开| 下载CSV

仿真中采用四阶Runge-Kutta方法积分C-W方程来计算相对轨道, 时间步长为$\Delta t=1$ s. 服务航天器的半径为1 m, 质量均为100 kg. 假定航天器的质心位于几何中心, 位置和距离传感器可以提供理想的信息; 通过三个主轴上各安装最大推力为1.0 N的推进器执行机动, 航天器的最大相对加速度不超过0.01 m/s$^{2}$, 最大相对速度不超过1 m/s. 卫星采用化学推进系统, 比冲$I_{sp} =300$ s, 脉冲施加次数没有限制.

3.1 基于碰撞概率的APF仿真分析

交会策略的目的是将补给星带进预定的范围内, 同时确保避免碰撞. 当补给航天器到达目标外表面的0.5 m范围内, 交会结束.加速度随着APF制导方法的计算结果自由变化, 没有限制. 假设3个固定障碍物放置在最差的情况下, 即追踪航天器运动的路径上, 直径为1.2 m. 服务航天器的状态标准差为球形高斯分布, 各个方向均相同. 初始状态导航偏差协方差为

$C_{p0}=\mbox{diag}\Big(1 \text{m}^2,1 \text{m}^2,1 \text{m}^2,0.1 \text{cm}^2/\text{s}^2,\\ 0.1 \text{cm}^2/\text{s}^2,0.1 \text{cm}^2/\text{s}^2\Big)$

3.1.1 静态障碍物

图2给出了采用两种不同APF控制方法的服务航天器轨迹图, 其中红色实线代表基于最小距离的APF控制方法、黑色虚线代表基于碰撞概率的APF控制方法. 从图2中可以看出, 两种轨迹差别较小.这一方面验证了基于瞬时碰撞概率的APF控制方法的可行性, 另一方面, 瞬时碰撞概率仅仅作为APF控制方法制导的门限, 对于静态障碍物的控制轨迹改进并不显著.

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 两种碰撞预警门限的航天器轨迹对比

Fig. 2 Comparison of spacecraft trajectory between two collision warning thresholds

3.1.2 动态障碍物

为进一步检验控制方法对复杂环境的有效性以及适应性, 将球形障碍改变为移动障碍: 其球心按照相对轨道运动模型进行运动, 球心初始位置矢量为($-$2.93 m, 14.89 m, 7.47 m), 速度矢量为(8.4 mm/s, 6.7 mm/s, 11 mm/s). 利用基于最小距离APF和基于瞬时碰撞概率的APF制导方法分别对服务航天器的交会任务进行控制, 同时统计机动过程的总推进剂消耗.

两种预警方法仿真得到的轨迹对比如图3所示, 图中红色粗实线为基于最小距离机动的服务航天器轨迹图, 推进剂消耗为0.049 kg;黑色细虚线为基于瞬时碰撞概率机动的服务航天器轨迹图, 推进剂消耗为0.023 kg. 从图3中可以明显看出, 基于最小距离机动的情况下出现虚假报警情况, 导致服务航天器进行了较多的轨道机动, 从而造成额外的推进剂消耗; 而基于碰撞概率的APF方法有效地减少了虚警率, 避免额外的机动, 减少了不必要的推进剂消耗.

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 服务航天器和动态障碍物相对目标航天器的运动轨迹

Fig. 3 Motion trajectories of service spacecraft and dynamic obstacles relative to target spacecraft

3.2 基于速度走廊设计目标航天器斥力的仿真分析

以增加目标斥力的APF方法进行仿真, 所获得的被控航天器轨迹如图4中红色细实线所示, 而黑色粗虚线为目标斥力不考虑速度安全接近走廊的服务航天器轨迹.

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 目标斥力考虑与不考虑速度安全接近走廊的航天器轨迹对比

Fig. 4 Comparison of spacecraft trajectory for Target repulsion with and without consider the speed safely approaching corridor

图5显示了有目标斥力与没有目标斥力的相对速度随时间的变化情况,图6显示了相对速度和相对距离的对应关系,其中红色细实线表示目标斥力考虑速度安全接近走廊的航天器接近速度,黑色粗虚线表示目标斥力不考虑速度安全接近走廊的航天器接近速度.从图5和图6中可以看出, 当被控航天器快接近目标航天器时,目标斥力考虑速度安全接近走廊可更有效减缓靠近目标的速度,满足速度安全接近走廊要求, 保证了接近的安全性. 同时,图4中目标斥力考虑比不考虑速度安全接近走廊的航天器任务总时间大大增加.

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 航天器接近速度随时间变化情况(目标斥力考虑与不考虑速度安全接近走廊)

Fig. 5 Spacecraft approaching speed vs time (target repulsion with and without consider the speed safely approaching corridor)

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 航天器接近速度随距离变化情况(目标斥力考虑与不考虑速度安全接近走廊)

Fig. 6 Spacecraft approaching velocity vs distance (target repulsion with and without consider the speed safely approaching corridor)

3.3 设计加速度的仿真分析

在本例中, 分别使用了具有上限和恒定变化率的加速度和多个档位开关的bang-bang加速度.设置初始任务参数与表1相同, 机动过程与交会机动相同. 获得的轨迹如图7所示.

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 不同加速度实现形式的航天器轨迹对比

Fig. 7 Comparison of spacecraft trajectories with different acceleration implementations

图8和图9显示了具有加速度上限和恒定变化率的模拟结果,图10和图11显示了多个档位开关的bang-bang加速度.从图8到图11可以看出加速度满足所设定的要求, 而且加速度在初始阶段迅速变化,因为引力加速度与距离有关. 在初始阶段,目标和受控制的航天器之间的距离相对较长, 对应引力加速度很大.由于受控制的航天器与障碍物之间的斥力增加, 机动过程的加速有几个显著的变化.图9和图11中速度的变化趋势与图5中所给出的历史变化趋势相似,都是先增大后逐渐趋近零, 验证了改进的APF方法的安全性和有效性. 然而,如图9所示, 当受控制的航天器接近目标时, 速度总是会略微振荡,对实际的任务操作是不利的. 与图9相比, 图11显示了更平滑的速度历史,这表明使用多个档位开关的bang-bang加速效果更好.

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 具有上限和常变化率加速度形式的加速度随时间变化

Fig. 8 Time vs the acceleration with upper limit and constant variation rate of acceleration

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 具有上限和常变化率加速度形式的速度随时间变化

Fig. 9 Time vs the speed with upper limit and constant variation rate of acceleration

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 多个档位开关bang-bang加速度形式的加速度随时间变化

Fig. 10 Time vs the acceleration with the bang-bang thrust with multi-speed switches of acceleration

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11 多个速度档位开关bang-bang加速度形式的速度随时间变化

Fig. 11 Time vs the speed with the bang-bang thrust with multi-speed switches of acceleration

针对两种加速度形式, 统计了交会过程中总任务时间和发动机开关切换次数,如表2所示. 根据表2, 这两种方法的总任务时间和发动机切换的次数是相似的,这表明使用多台小发动机来代替一台大发动机能带来更平滑的速度剖面,但不需要额外的任务时间或更多地发动机开关负荷.

表2 不同加速度形式机动过程的统计参数

Table 2 Statistical parameters of maneuvering processes with different acceleration forms

新窗口打开| 下载CSV

4 结论

本文提出了一种基于人工势函数控制方法改进的航天器近距离操作安全控制方法.首先, 将碰撞概率计算纳入了判断碰撞的标准, 避免虚假预警, 提高决策准确度;其次, 利用目标的速度安全裕度和控制边缘来制定目标排斥力,以减慢航天器逼近的速度. 第三, 设计两种实用的推力形式,包括带有上限和恒定变化率的推力, 以及多个档位开关的bangbang推力,用来代替已有研究中使用的连续变量推力形式.

结果表明, 改进的基于人工势函数的方法能够成功地帮助航天器安全完成近距离操作,采用碰撞概率的评判可以减少额外的轨道机动, 节省燃料;考虑目标的速度安全裕度和控制边缘设计目标斥力, 能保证安全的相对接近速度.与带有上限和恒定变化率的推力形式相比,多个档位开关的bang-bang推力形式可以提供一个较为平滑的速度.本文设计的安全控制算法方法可以提高航天器近距离操作的安全性和有效性.

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Schulte

, Spencer

Development of an integrated spacecraft guidance, navigation, and control subsystem for automated proximity operations

Acta Astronautica, 2016,118:168-186

DOI URL [本文引用: 1]

[2]

Zhang

, Chu

, Zhang

, et al.

Safe-trajectory optimization and tracking control in ultra-close proximity to a failed satellite

Acta Astronautica, 2018,144:339-352

DOI URL

[3]

Romano

, Friedman

, Shay

Laboratory experimentation of autonomous spacecraft approach and docking to a collaborative target

Journal of Spacecraft and Rockets, 2017,1(44):164-173

DOI URL [本文引用: 1]

[4]

罗刚桥

美俄碰撞分析与启示

国际太空, 2009(6):23-27