基于分段常值推力的水滴悬停构型控制策略

2019-04-08白晟州王慧疆韩潮张斯航

北京航空航天大学学报 2019年3期

白晟州，王慧疆，韩潮，张斯航

(北京航空航天大学宇航学院，北京 100083)

随着空间领域的研究、开发以及应用的不断提高，航天器功能与结构日趋复杂，航天器在轨服务技术可以有效地保证航天器在复杂的空间环境中持久、稳定、高质量地在轨运行，因而成为当前空间技术研究的热门[1-6]。航天器在轨服务技术主要包含在轨检查、交会对接和编队飞行等，其中涉及的一个核心问题是航天器的绕飞问题，即通过对任务航天器施加脉冲或推力，使其绕目标航天器作近距离周期运动。根据任务航天器绕目标航天器运动的周期与目标航天器本身运动的周期的比值，可以将绕飞分为“快速绕飞”与“慢速绕飞”。从更广的意义上说，通过对任务航天器施加脉冲或推力，使其与目标航天器距离始终保持在一定范围内，也可以称作绕飞。

作为经典相对运动方程，C-W(Clohessy-Wiltshire)方程得到了广泛的关注与应用。赵书阁和张景瑞[7]基于C-W方程研究了航天器共面圆型快速绕飞问题并相应分析了2种典型方法；林来兴[8]研究了绕飞轨道动力学和稳定性；师鹏等[9]基于线性动力学模型，分析了有限推力下的航天器绕飞性质；潘屹[10]对C-W方程的状态转移矩阵进行了推导，给出了选择悬停轨道的方法。这些理论进一步完善和发展了C-W方程理论，但工程应用价值有待提高。为了实现任务航天器对目标航天器的悬停或绕飞，设计基于脉冲控制和小推力控制的绕飞构型策略成为热门研究。Straight等[11]提出了以圆形轨道参考卫星为中心的特定环形区域内的特定点的脉冲制导方案； Hope和Trask[12]提出一种脉冲控制下的水滴悬停构型；Lovell和Tollefson[13]进一步给出了水滴悬停构型的参数表示方法；饶殷睿等[14]基于相对轨道要素描述了水滴悬停构型；王功波等[15]在圆参考轨道和连续小推力条件下，推导了快速绕飞策略；罗建军等[16]分析了快速受控绕飞；朱小龙等[17]提出了一种参数延拓方法，实现了有限推力bang-bang控制下的绕飞轨迹优化问题；张冉等[18]推导了4种受迫绕飞构型的解析表达式和脉冲控制策略，完善了受迫绕飞构型设计理论。

上述方法均可实现悬停或受迫绕飞，但模型构建存在一定的局限性。大部分研究都是基于圆轨道假设，即参考轨道是圆轨道，无法有效地解决非圆参考轨道下的受迫绕飞问题；对于悬停构型，大多采用脉冲推力或连续小推力实现构型的保持，工程应用难度较大。针对上述情况，本文提出了多段常值推力控制实现水滴悬停构型的打靶方程，分析了近距离相对运动条件下两段常值推力控制的可行性，数值仿真显示分段常值小推力可以实现水滴悬停相对运动，与脉冲推力或连续小推力控制相比，更加符合工程实际。

1 水滴悬停构型

水滴悬停构型是航天器悬停构型中一种典型构型[19]，可同时满足悬停和高精度要求。将构型建立在质心非惯性坐标系中，如图1所示。

图1 水滴悬停构型三维示意图Fig.1 Schematic diagram of 3D teardrop hovering configuration

1.1 水滴悬停打靶方程

给出经典轨道摄动方程：

(1)

根据水滴悬停构型，可以得到施加脉冲前后任意时刻对应的轨道要素，不妨设施加脉冲前某时刻为t0，施加脉冲后的某时刻为tf，则对应的轨道要素分别为为X(t0)和X(tf)，若在X(t0)处施加多段常值推力，使其在tf时刻轨道要素刚好为X(tf)，则实现了用多段常值推力代替脉冲实现水滴悬停构型。

假定任务航天器在一个水滴周期内通过N段常值推力维持悬停构型，则常值小推力f可表示为

f=f(t1,t2,…,tN-1,Ft1,Fn1,Fh1,…，FtN,FnN,

FhN)=f(t1,t2,…,tN-1,μ1,μ2，…，μN)

(2)

根据上述分析，可建立相应的打靶方程。

(3)

1.2 最小二乘法求解常值推力

因为摄动方程是高度非线性的微分方程组，无法直接求解常值推力解析解，同时每段推力的作用时间也未知，也是待求量，所以将原方程问题转化为极值问题求解，考虑到构造极值问题：

by optf,t1,t2,…,tN-1

subject to:

(4)

显然，优变量随着推力段数的增加而增加，相应的计算时间也会大量增加，计算精度不能得到保证，且推力段数过多导致控制复杂，不利于工程实际应用。应当考虑采用尽可能少的推力段数实现常值推力控制。

2 单段常值推力与多段常值推力

首先分析只采用一段常值推力下的情况：

(5)

式中：ft,fn,fh为t时刻μ1在3个方向上的投影分量。将μ1看作参数，根据常微分方程解对参数的连续依赖性，X(μ1,tf)为关于μ1的可微连续函数，即X(μ1,tf)局部上是一个三维微分流形，而X(μ1,tf)的值本身是六维空间的元素。显然，三维流形无法覆盖六维空间，也就是说，很多情况下，不存在满足条件的μ1，使得

X(tf)=X(μ1,tf)=X(μ1,t0)+

(6)

其次，考虑两段常值推力，假设两段常值推力作用时间相同，则代入最小二乘优化函数里，优化变量只有两段常值推力，即六维变量。小邻域定理保证了解的局部存在性。

小邻域定理：近距离运动假设下，固定时间内，对于参考轨道，存在一个小邻域对其中任意一点，一定存在两段常值推力解。

数学描述如下：

设X1(μ1,μ2,tf)是两段常值推力轨道，在t0时刻与参考轨道X0(t0)相等，在[t0,(t0+tf)/2]内推力为μ1；在[(t0+tf)/2,tf]内推力为μ2；对于tf时刻参考轨道瞬时要素X(tf)，一定存在以其为圆心的一个邻域B(X(tf),δ)(见图2)，对其中任意一点Xtf，一定存在μ1,μ2，使得X1(μ1,μ2,tf)在tf时刻的轨道要素恰好是Xtf。

证明根据小推力线化方程[20]，有

(7)

式中：E0、Eh和Ef为参考轨道X0(t)分别在t0、(t0+tf)/2和tf时的偏近点角；A(E0,Eh)与A(Eh,Ef)均表示线化矩阵，具体表示详见文献[20]在近距离运动假设下，小推力线化方程的系数矩阵C可以近似代替雅可比阵。detC≠0，所以雅可比阵行列式不为零，根据逆映射定理，一定存在一个小邻域满足一一映射。

证毕

根据小邻域定理，只要设计的推力曲线终端六要素位于以参考轨道终端时刻的六要素为圆心的小球内，常值推力一定有解，因此，若设计的推力曲线终端六要素与参考轨道终端时刻的六要素满足近距离相对运动假设，很大可能存在一对常值推力满足：

图2 小邻域定理示意图Fig.2 Schematic diagram of small neighborhood theorem

(8)

其中:推力f满足：

(9)

3 修正公式

采用最小二乘法得到的双段常值推力解往往精度较差，考虑采用迭代方法对解进行修正，提高精度准确性。

设映射函数F:Rn→Rn。

若存在根x0，使得F(x0)=0，且F(x0)的Jacobi矩阵JF(x0) 满秩，若矩阵A≈JF(x0)。

证明对任意x,y∈B(x0,δ)，定义：

g(x)=x-A-1F(x)，g(y)=y-A-1F(y)

g(x)-g(y)=(x-y)-A-1(F(x)-F(y))=

(x-y)-A-1JF(y)(x-y)-A-1L(x-y)=

(I-A-1JF(y))(x-y)-A-1L(x-y)

这是因为当x和y与x0接近，A-1JF(y)≈I,又有，若xk∈δ1，xk+1=xk-A-1F(xk)，则

xk+1-x0=xk-x0-A-1F(xk)-

A-1F(x0)=(I-A-1JF(x0))(xk-

x0)+m(xk-x0)xk-x0→0,m→0

由公式：

得

(10)

故xk+1∈B(x0,δ)。

证毕

(11)

可以提高最小二乘解的精度。

4 数值仿真

设参考星K的轨道要素为X=(a,e,i,Ω,ω,M)=(7 555 000 m，0，π/4，0，0，0)，其中半长轴单位为m。任务航天器W绕K作水滴悬停运动，水滴构型几何参数为(xhover,yhover,zhover,du)=(1 000 m，0，2 000 m，π/3)，其中du表示水滴构型维持一周时参考星K绕地球转过的弧度，用以表征水滴构型的周期。

首先，根据方程式(4)，采用最小二乘法尝试寻找一组多段推力解(μ1,μ2,μ3,μ4,μ5)。常值推力时间Δt为参考星K绕地球转过0.124 rad所需时间。代入方程式(4)使用最小二乘法求解，并进行了计算机仿真，仿真结果如图3所示。

在(x,z)=(930,275 0) m时，任务航天器开始施加常值推力，直至(x,z)=(1 100,2 750) m时停止施加推力，其中一共采用了5段常值推力，推力结束后，任务航天器进入自由段，不再施加推力，直到下一个周期/往复循环，表明确实存在多段常值推力控制可以代替脉冲推力使目标航天器绕参考星进行周期性水滴悬停运动，同时也说明最小二乘法可以较好地求解多段常值推力。

设参考星K的轨道要素仍然为X=(a,e,i,Ω,ω,M)=(7 555 000 m，0，π/4，0，0，0)，半长轴单位为m，任务航天器W绕K作水滴悬停运动，水滴构型几何参数为(xhover,yhover,zhover,du)=(1 000 m，0，1 000 m，2π/3)。常值推力时间Δt为参考星K绕地球转过π/9所需时间。

尽管5段常值推力实现了周期性水滴悬停运动，但计算时间较长，控制较复杂。为了减少求解时间，得到更简单的推力控制策略，考虑两段常值推力实现悬停构型，即在N=2的情况下进行了仿真，首先通过最小二乘法得到双推力初值，得到的结果精度差于5段常值推力控制，然后使用精度法对两段推力解进行精度修正，修正后的解精度较好，最后仿真结果如表1和表2所示。

图4为仅用最小二乘法得到的两段推力解。绿色轨迹代表施加脉冲推力的仿真效果，蓝色段表示为无推力段，红色段代表施加两段常值推力作用后的结果。可以看出，两段常值推力基本实现了水滴悬停构型控制，但控制精度较差，与原先构型重合度较差。

如图5所示，绿色轨迹代表施加脉冲推力的仿真效果，蓝色段表示为无推力段，红色段代表施加两段常值推力作用后的结果。在(x,z)=(-300,3 500) m时，任务航天器开始施加常值推力，直至(x,z)=(2 250,3 300)m时停止施加推力，推力结束后，任务航天器进入自由段，不再施加推力，直到下一个周期往复循环。可以看到，红色段与绿色段重合度较高，说明双段常值推力解的精度较高，可以代替脉冲推力进行水滴悬停构型控制。