Low-thrust transfer to Backflip orbits

The aim of the work is to design a low-thrust transfer from a Low Earth Orbit to a “useful” periodic orbit in the Earth–Moon Circular Restricted Three Body Model (CR3BP). A useful periodic orbit is here intended as one that moves both in the Earth–Moon plane and out of this plane without any requirements of propellant mass. This is achieved by exploiting a particular class of periodic orbits named Backflip orbits, enabled by the CR3BP. The unique characteristics of this class of periodic solutions allow the design of an almost planar transfer from a geocentric orbit and the use of the Backflip intrinsic characteristics to explore the geospace out of the Earth–Moon plane. The main advantage of this approach is that periodic plane changes can be obtained by performing an almost planar transfer. In order to save propellant mass, so as to increase the scientific payload of the mission, a low-powered transfer is considered. This foresees a thrusting phase to gain energy from a departing circular geocentric orbit and a second thrusting phase to match the state of the target Backflip orbit, separated by an intermediate ballistic phase. This results in a combined application of a low-thrust manoeuvre and of a periodical solution in the CR3BP to realize a new class of missions to explore the Earth–Moon neighbourhoods in a quite inexpensive way. In addition, a low-thrust transit between two different Backflip orbits is analyzed and considered as a possible extension of the proposed mission. Thus, also a Backflip-to-Backflip transfer is addressed where a low-powered probe is able to experience periodic excursions above and below the Earth–Moon plane only performing almost planar and very short transfers.     

一种适用于平动点周期轨道初值计算的简化路径搜索修正法

在限制性三体问题中，路径搜索修正法是一种基于平动点周期轨道垂直穿越Poincare截面的几何对称性计算平面及空间平动点周期轨道近似初值的方法.采用路径搜索修正法的一种简化形式，在圆形限制性三体模型中，对地月系中几种典型的平面及空间周期轨道近似初值进行了计算.结果表明，该简化方法得到的周期轨道近似初值不唯一，由近似初值经微分修正得到的精确结果中通常同时存在Halo轨道和大幅值逆行轨道（DRO）.进一步分析表明，在某些临界初值下，精确结果中Halo轨道将消失，同时可能出现平面Lyapunov轨道及Vertical轨道.上述计算中，搜索初值与结果中轨道类型的对应关系值得进一步研究.      

基于火卫二视线测量的火星接近段自主导航算法

针对火星接近段导航通信时延大、存在通信盲区、自主导航可用观测信息有限等问题,提出了一种基于一组火卫二相对探测器视线矢量测量的天文自主导航算法,每个导航周期测量一组火卫二视线矢量可得到其中某一时刻探测器的完整轨道信息的估计值,该方法不依赖于轨道渐近线方向等先验信息。考虑到存在火卫二和火星在同一视场的情形,此时结合火星中心视线矢量方向以及火卫二的星历信息可估计出精度较高的探测器轨道半径,作为第一种方法的补充观测量,提高导航精度。最后给出仿真校验,验证了该方法的导航精度和可行性,表明该方案能够满足未来火星探测接近段自主导航需求。     

GEO螺旋巡游轨道的相对运动分析和螺旋环设计

摘要： GEO螺旋巡游轨道采用螺旋巡游方式，以不同的构型“上下浮动”在GEO轨道附近，可实现对该轨道上空间目标和空间环境的高精度探测.本文在分析GEO轨道航天器运动规律的基础上，应用小偏差理论分析螺旋巡游轨道与GEO目标之间的相对运动，给出平面螺旋环和三维螺旋环的设计方法，为GEO螺旋巡游轨道的设计奠定基础.     

Near real-time BDS GEO satellite orbit determination and maneuver analysis with reversed point positioning

The Geostationary Earth Orbit (GEO) satellite is a crucial part of the BeiDou Navigation Satellite System (BDS) constellation. However, due to various perturbation forces acting on the GEO satellite, it drifts gradually over time. Thus, frequent orbit maneuvers are required to maintain the satellite at its designed position. During the orbit maneuver and recovery periods, the orbit quality of the maneuvered satellite computed with broadcast navigation ephemeris will be significantly degraded. Furthermore, the conventional dynamic Precise Orbit Determination (POD) approach may not work well, because of a lack of publicly available satellite information for modeling the thrust forces. In this paper, a near real-time approach free of thrust forces modeling is proposed for BDS GEO satellite orbit determination and maneuver analysis based on the Reversed Point Positioning (RPP). First, the station coordinates and receiver clock offsets are estimated by GPS/BDS combined Single Point Positioning (SPP) with single-frequency phase-smoothed pseudorange observations. Then, with the fixed station coordinates and receiver clock offsets, the RPP method can be conducted to determine the GEO satellite orbits. When no orbit maneuvers occur, the proposed method can obtain orbit accuracies of 0.92, 2.74, and 8.30?m in the radial, along-track, and cross-track directions, respectively. The average orbit-only Signal-In-Space Range Error (SISRE) is 1.23?m, which is slightly poorer than that of the broadcast navigation ephemeris. Using four days of GEO maneuvered datasets, it is further demonstrated that the derived orbits can be employed to characterize the behaviors of GEO satellite maneuvers, such as the time span of the maneuver as well as the satellite thrusting accelerations. These results prove the efficiency of the proposed method for near real-time GEO satellite orbit determination during maneuvers.     

Solar sails for perturbation relief: Application to asteroids

A major cause of spacecraft orbital variation comes from natural perturbations, which, in close proximity of a body, are dominated by its non-spherical nature. For small bodies, such as asteroids, these effects can be considerable, given their uneven (and uncertain) mass distribution. Solar sail technology is proposed to reduce or eliminate the net secular effects of the irregular gravity field on the orbit. Initially, a sensitivity analysis will be carried out on the system which will show high sensitivity to changes in initial conditions. This presents a challenge for optimisation methods which require an initial guess of the solution. As such, the Genetic Algorithm (GA) is proposed as the preferred optimisation method as this requires no initial guess from the user. A multi-objective optimisation is performed which aims to achieve a periodic orbit whilst also minimising the effort required by the sail to do so. Given the system sensitivity, the control law for one orbit is not necessarily applicable for any subsequent orbit. Therefore, a new method of updating the control law for subsequent orbits is presented, based on linearisation and use of a Control Transition Matrix (CTM). The techniques will later find application in a multiple asteroid rendezvous mission with a solar sail as the primary propulsion system.     

Orbit determination by genetic algorithm and application to GEO observation

This paper demonstrates an initial orbit determination method that solves the problem by a genetic algorithm using two well-known solutions for the Lambert’s problem: universal variable method and Battin method. This paper also suggests an intuitive error evaluation method in terms of rotational angle and orbit shape by separating orbit elements into two groups. As reference orbit, mean orbit elements (original two-lines elements) and osculating orbit elements considering the J2 effect are adopted and compared. Our proposed orbit determination method has been tested with actual optical observations of a geosynchronous spacecraft. It should be noted that this demonstration of the orbit determination is limited to one test case. This observation was conducted during approximately 70 min on 2013/05/15 UT. Our method was compared with the orbit elements propagated by SGP4 using the TLE of the spacecraft. The result indicates that our proposed method had a slightly better performance on estimating orbit shape than Gauss’s methods and Escobal’s method by 120 km. In addition, the result of the rotational angle is closer to the osculating orbit elements than the mean orbit elements by 0.02°, which supports that the estimated orbit is valid.     

基于分段常值的全电推进GEO卫星制导策略

电推进技术因其比冲高的技术特点在GEO轨道转移中应用可大大减少燃料质量，提高有效载荷质量比，延长任务寿命等。针对全电推进GEO卫星入轨的轨迹优化和制导问题，首先利用间接法获得小推力燃料最优GEO轨道转移的数值解，提出一种多项式曲线拟合最优轨迹的方法，多项式曲线形式简单，可作为参考轨道在星上存储和使用。在多项式参考轨道的基础上，建立了一种分段常值推力跟踪参考轨道的闭环制导策略，在常值推力条件下，轨道要素控制量与控制力有解析关系，简化了制导律设计；将多圈轨道转移问题分解为多个单圈轨道优化问题。结果显示，本文提出的分段常值跟踪制导策略跟踪精度高，和最优轨道相比多消耗7%的燃料。本制导策略控制结构简单，易于工程实施。     

电磁航天器编队悬停鲁棒协同控制方法

针对电磁航天器编队近地轨道悬停问题,提出一种在缺少参考轨道准确信息时的协同控制方法。用TH方程描述航天器间的相对运动,选择与参考轨道同周期的圆轨道为标称轨道。将参考轨道相对于标称圆轨道的偏差、地球非球形引力、大气阻力及其他天体引力等参数单独归类,视其为不确定量,构成不确定系统。通过引入一致性理论,在电磁作用模型和动力学方程均存在不确定性的条件下,针对航天器编队悬停的目标设计了鲁棒协同控制律。考虑能量消耗最优和均衡以及轨道姿态解耦,给出了通过优化进行磁矩配置的方案。仿真结果表明,所设计的鲁棒协同控制律能够实现编队电磁航天器高精度悬停,所给出的磁矩配置方案能够实现磁矩的合理分配。      

基于差分进化算法的导航星座在轨重构构型设计

针对有卫星失效下的导航星座,提出了对在轨卫星进行相位机动的方法来对星座的空间构型进行重构,以实现修复和改善星座性能的目的。首先,利用共面高轨和共面低轨变相的方式来对在轨卫星的相位进行调整,建立共面轨道相位机动的数学模型;其次,提出了除重构时间和重构能量外的其他重构指标,包括重构能量均衡度和重构构型恢复性,并建立了各重构指标的数学模型;然后,建立了重构构型的优化模型,并对优化问题中个体的评价手段与编码方式进行设计;最后,以北斗中MEO卫星失效为例,利用差分进化算法对优化模型进行求解,得到以不同重构指标为目标函数下的Pareto前沿。从结果中可以看出该重构方法对星座性能的提升最大可以达到41.2%,同时Pareto前沿中对应的所有重构策略中,重构能量、重构构型的均衡度、重构构型的健壮性的数量级均维持在-1、-3和0的水平。     

基于GPS的地球资源卫星精密星历外推新方法

随着技术的发展,通过星载GPS接收机直接确定卫星星历成为卫星定位的一个重要手段.GPS接收机获取的卫星星历数据是某一时刻的瞬时状态,要获取连续的卫星星历数据还需要进一步处理.常用的处理方法有几何法与动力学法.在GPS接收机给定瞬时星历频率较低的情况下,几何法的计算误差比较大,特别是只有一组瞬时星历时,无法用几何法进行轨道的外推.在分析地球资源卫星轨道特点的基础上,提出一种新的轨道缩减动力模型,该模型将卫星运动在直角坐标系中分解为简谐运动,利用模型实现了轨道外推的算法.通过试验验证,该算法可以达到较高的精度.      

Orbit determination of BeiDou navigation satellite system maneuvered satellites based on instantaneous velocity pulses parameters

The BeiDou navigation satellite system (BDS) comprises geostationary earth orbit (GEO) satellites as well as inclined geosynchronous orbit (IGSO) and medium earth orbit (MEO) satellites. Owing to their special orbital characteristics, GEO satellites require frequent orbital maneuvers to ensure that they operate in a specific orbital window. The availability of the entire system is affected during the maneuver period because service cannot be provided before the ephemeris is restored. In this study, based on the conventional dynamic orbit determination method for navigation satellites, multiple sets of instantaneous velocity pulses parameters which belong to one of pseudo-stochastic parameters were used to simulate the orbital maneuver process in the orbital maneuver arc and establish the observed and predicted orbits of the maneuvered and non-maneuvered satellites of BeiDou regional navigation satellite system (BDS-2) and BeiDou global navigation satellite system (BDS-3). Finally, the single point positioning (SPP) technology was used to verify the accuracy of the observed and predicted orbits. The orbit determination accuracy of maneuvered satellites can be greatly improved by using the orbit determination method proposed in this paper. The overlapping orbit determination accuracy of maneuvered GEO satellites of BDS-2 and BDS-3 can improve 2–3 orders of magnitude. Among them, the radial orbit determination accuracy of each maneuvered satellite is basically better than 1 m. simultaneously, the combined orbit determination of the maneuvered and non-maneuvered satellites does not have a great impact on the orbit determination accuracy of the non-maneuvered satellites. Compared with the multi GNSS products (indicated by GBM) from the German Research Centre for Geosciences (GFZ), the impact of adding the maneuvered satellites on the orbit determination accuracy of BDS-2 satellites is less than 9 %. Furthermore, the orbital recovery time and the service availability period are significantly improved. When the node of the predicted orbit is traversed approximately 3 h after the maneuver, the accuracy of the predicted orbit of the maneuvered satellite can reach that of the observed orbit. The SPP results for the BDS reached a normal level when the node of the predicted orbit was 2 h after the maneuver.     

火卫一周期准卫星轨道及入轨分析

围绕火卫一的准卫星轨道（QSOs）因其具有良好的稳定性，是火卫一探测任务最为实用的轨道。在平面圆型限制性三体问题模型下，利用庞加莱截面和KAM环迭代方法探究了准卫星轨道的周期轨道族，并给出不同能量准卫星周期轨道的初始条件。针对火卫一周期准卫星轨道入轨，提出一种转移轨道设计方法:对准卫星周期轨道调整速度后进行反向积分，直至离开火卫一邻近区域，从而得到由火星环绕轨道向火卫一周期准卫星轨道的转移轨道，并调整转移轨道参数对燃料与时间消耗进行优化。研究结果表明，当周期准卫星轨道能量处于特定区间时，存在特定速度脉冲区间，可利用火卫一引力实现较少燃料消耗的轨道转移；在该速度脉冲区间中，通过选取较小的速度脉冲，可缩短转移时间。      

双臂空间机器人关节运动的一种增广鲁棒控制方法

讨论了载体位置与姿态均不受控制情况下,漂浮基双臂空间机器人系统的控制问题.结合系统动量守恒关系进行的运动学、动力学分析表明,可以得到一组与适当选择的惯性参数呈线性函数关系的、欠驱动形式的系统动力学方程.以此为基础,并采用增广变量的思想,克服了通常情况下,空间机器人系统动力学方程关于系统惯性参数呈非线性关系的难点,针对双臂空间机器人末端爪手所持载荷参数不确定,但误差范围可确定的情况,设计了漂浮基双臂空间机器人关节运动的变结构鲁棒控制方案.该控制方案的优点在于,不需要反馈、测量漂浮基的位置、移动速度及移动加速度;与自适应控制方案相比,化积分运算为简单四则运算,计算量大为减少,有利于实时应用.通过对一个平面双臂空间机器人系统的数值仿真,证实了算法的有效性.      

地月空间的远距离逆行轨道族及其分岔研究

地月系统中存在着一类绕月逆行、高度稳定的轨道族,称为远距离逆行轨道族(DRO)。以圆型限制性三体问题(CR3BP)为动力学模型研究了DRO轨道族周边的动力系统结构。利用Broucke稳定性图寻找分叉点,判断分叉类型,基于数值延拓计算分岔后产生的一系列新轨道分支。分叉类型主要有切分叉与多倍周期分叉(从3倍周期开始),轨道维度包含平面轨道族与三维轨道族。计算新轨道族的特征,包括形状、周期、能量、稳定性、双曲流形结构等。探讨周期轨道的轨道周期与能量的关系,以几何化的方式展现分叉结构、多周期轨道的双曲流形结构等。该动力结构将为基于DRO轨道族的地月空间任务提供重要的理论支持。      

基于引力场不对称性的三体系统轨道自主导航

以月球背面的中继通信为背景,提出了基于三体系统引力场不对称特性的星–星测距自主定轨方案。该方案以环月极轨卫星和地–月L2点Halo轨道卫星组成中继通信网,以实现对月球两极和背面的覆盖。通过采集极轨卫星与Halo轨道卫星的测距信息,结合卡尔曼滤波在日–地–月动力学模型下获得两颗卫星的绝对轨道。数值仿真结果表明:本文方法能将导航的位置精度和速度精度分别提高到百米和厘米/秒量级。该自主导航方法还可以扩展到不规则引力场小天体附近星群运动的自主导航。     

基于经验加速度的低轨卫星轨道预报新方法

研究将定轨过程中的经验加速度应用于地球低轨卫星轨道预报的新方法. 利用GPS伪距观测数据和简化动力学最小二乘批处理方法对地球低轨卫星定 轨, 其中卫星位置、速度及大气阻力系数和辐射光压系数可以直接用于轨道预报. 作为简化动力学最重要特征的经验加速度呈现准周期、余弦曲线特点, 可通过 傅里叶级数拟合建模. 确定性动力学模型与补偿大气阻力模型误差的切向经验 加速度级数拟合模型组成增强型动力学模型用于提高轨道预报精度. 应用 GRACE-A星载GPS伪距观测数据和IGS超快星历定轨并进行轨道预报, 结果表明 轨道预报初值位置精度达到0.2m, 速度精度达到1×10-4m·s-1, 预报3天位置精度优于60m, 比只利用确定性动力学模型进行预报精度平 均提高2.3倍. 先定轨后预报的模式可用在星上自主精确导航系统中.      

Dynamics of a spacecraft and normalization around Lagrangian points in the Neptune–Triton system

The problem of a spacecraft orbiting the Neptune–Triton system is presented. The new ingredients in this restricted three body problem are the Neptune oblateness and the high inclined and retrograde motion of Triton. First we present some interesting simulations showing the role played by the oblateness on a Neptune’s satellite, disturbed by Triton. We also give an extensive numerical exploration in the case when the spacecraft orbits Triton, considering Sun, Neptune and its planetary oblateness as disturbers. In the plane a × I (a = semi-major axis, I = inclination), we give a plot of the stable regions where the massless body can survive for thousand of years. Retrograde and direct orbits were considered and as usual, the region of stability is much more significant for the case of direct orbit of the spacecraft (Triton’s orbit is retrograde). Next we explore the dynamics in a vicinity of the Lagrangian points. The Birkhoff normalization is constructed around L₂, followed by its reduction to the center manifold. In this reduced dynamics, a convenient Poincaré section shows the interplay of the Lyapunov and halo periodic orbits, Lissajous and quasi-halo tori as well as the stable and unstable manifolds of the planar Lyapunov orbit. To show the effect of the oblateness, the planar Lyapunov family emanating from the Lagrangian points and three-dimensional halo orbits are obtained by the numerical continuation method.     

The influence of orbital maneuver on autonomous orbit determination of an extended satellite navigation constellation

The right ascension of the ascending node is unobservable if only the inter-satellite ranging is used for autonomous orbit determination (AOD) of an Earth navigation constellation. However, if an Earth-Moon libration point satellite is added to the Earth navigation constellation to construct an extended navigation constellation, all the orbital elements can be determined with only the inter-satellite ranging. Furthermore, the extended navigation constellation can provide navigation information for interplanetary probes. For such an extended navigation constellation, orbital control needs to be considered due to the instability of the libration-point satellite orbit. This study concerns the influence of satellite orbital maneuver on the AOD of the extended navigation constellation. An AOD method under orbital maneuver is proposed. A low thrust controller is designed to achieve libration point satellite autonomous orbit maintenance by using AOD results. A navigation constellation consisting of three GPS satellites and one libration point satellite are designed for simulation. The simulation results show that libration point satellites can achieve autonomous navigation and autonomous orbit maintenance by only using inter-satellite ranging information. The rotation drift error of the Earth navigation constellation is also suppressed.     

基于多模型最优融合的双星定位系统一体化精密定轨方法

针对卫星动力学模型复杂且不准,考虑到卫星定轨中待估参数在时间和空间上的相关性,提出了一种基于多模型最优融合的双星定位系统一体化参数建模的近地卫星精密定轨新方法.利用节点自由分布B样条描述卫星运动,实现了对卫星粗略动力学模型的抑制作用;同时结合双星观测模型,使该方法转化为关于求解卫星轨道样条表示参数和定轨系统误差的多模型融合的非线性优化问题;通过引入模型结构确定最优融合权值的选取准则,在最小二乘准则下,采用非线性最优化方法搜寻样条的最优节点分布,得到了待估参数的最优估计,完成了近地卫星的精密定轨.理论分析和仿真计算表明,该方法确实有效,不仅提高了卫星的定轨精度,而且使状态估计的结构更加稳定.