航天器交会远程快速自主制导方法

针对近地近圆轨道航天器交会任务,设计了基于经典轨道要素的远程快速自主制导算法.对于任意初始纬度幅角偏差的远程导引,通过建立纬度幅角偏差与半长轴偏差的关系,将远程导引段分为初始轨道飞行、调相轨道飞行和调整轨道飞行3个阶段.初始轨道飞行进行轨道共面修正和调相机动;在调相轨道飞行期间,进行自然调相以及调相轨道到调整轨道的机动;调整轨道飞行阶段进行追踪航天器的远地点高度和近地点高度的修正,以及再次共面修正.所有变轨机动都以制导脉冲的形式给出,并都在轨道特殊点执行.精确轨道仿真验证了远程快速自主接近制导算法的可行性.     

Semi-analytical investigations of high area-to-mass ratio geosynchronous space debris including Earth’s shadowing effects

This paper investigates the long-term perturbations of the orbits of geosynchronous space debris influenced by direct radiation pressure including the Earth’s shadowing effects. For this purpose, we propose an extension of our homemade semi-analytical theory [Valk, S., Lemaître, A., Deleflie, F. Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence. Adv. Space Res., submitted for publication], based on the method developed by Aksnes [Aksnes, K. Short-period and long-period perturbations of a spherical satellite due to direct solar radiation. Celest. Mech. Dyn. Astron. 13, 89–104, 1976] and generalized into a more convenient non-singular formalism. The perturbations accounting for the direct radiation pressure with the Earth’s shadow are computed on a revolution-by-revolution basis, retaining the original osculating Hamiltonian disturbing function. In this framework, we compute the non-singular mean longitude at shadow entry and shadow exit at every orbital revolution in opposition to classical approaches where the singular eccentric anomalies at shadow entry and shadow exit are computed. This new algorithm is developed using non-singular variables. Consequently, it is particularly suitable for both near-circular and near-equatorial orbits as well as orbits which transit periodically around null eccentricities and null inclinations.The algorithm is tested by means of numerical integrations of the equations, averaged over the short periods, including radiation pressure, J₂, the combined Moon and Sun third body attraction as well as the long-term effects of the 1:1 resonance occurring for geosynchronous objects. As an extension of [Valk, S., Lemaître, A., Anselmo, L. Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios influenced by solar radiation pressure. Adv. Space Res., doi:10.1016/j.asr.2007.10.025, 2007b], we especially apply our analysis to space debris with area-to-mass as high as 20 m²/kg. This paper provides numerical and semi-analytical investigations leading to a deep understanding of the long-term evolution of the semi-major axis. Finally, these semi-analytical investigations are compared with accurate numerical integrations of the osculating equations of motion over time scales as high as 25 years.     

Configuration space and stability analysis of solar sail near-vertical earth-trailing orbits

Regions outside the reach of traditional propulsion systems or the ones that require significant propellant, may be reached by harnessing the solar radiation pressure and leveraging coupled dynamics to maneuver a sail-based spacecraft. Earth-trailing orbits have recently been investigated for getting a unique perspective of the Sun while maintaining the spacecraft in close proximity to Earth. Vertical orbits trailing the Earth exhibit the additional capability to view the Sun from above and below the ecliptic plane. In this work, families of sail-based orbits are explored for varying Earth-trailing angles and Z amplitudes in the Sun-Earth circular restricted three-body problem. Optimization is carried out to ensure that the non-traditional vertical orbits exhibit a constant pitch angle control history, as well as symmetry across the X-Y plane. The stability of the resulting orbit families is assessed using an extension of Flouquet theory to Differential Algebraic Equations. Results indicate that sail-based Earth-trailing vertical orbits can be more stable than traditional sub-L₁ sail-based vertical orbits.     

基于GLONASS星历的预报轨道的误差分析

推导了协议地球坐标系下的卫星运动方程.通过分析由GLONASS(Global Navigation Satellite System)广播星历参数确定的卫星预报轨道的拟合精度,指出了摄动力模型的简化、积分器的选择,以及忽略了极移影响等因素是引起拟合误差的主要因素,其中摄动力模型的简化起最主要的作用.通过对卫星轨道运动方程积分30?min,可知由摄动力模型的简化、积分器和忽略极移影响等因素引起的拟合误差分别为0.827?m,0.224?m和0.025?m.要提高预报轨道拟合的精度,关键是要对摄动力简化特别是地球引力摄动高阶项的截断以及日月引力场简化造成的轨道预报精度损失加以控制.     

Multiobjective genetic optimization of Earth–Moon trajectories in the restricted four-body problem

In this paper, optimal trajectories of a spacecraft traveling from Earth to Moon using impulsive maneuvers (ΔV maneuvers) are investigated. The total flight time and the summation of impulsive maneuvers ΔV are the objective functions to be minimized. The main celestial bodies influencing the motion of the spacecraft in this journey are Sun, Earth and Moon. Therefore, a three-dimensional restricted four-body problem (R4BP) model is utilized to represent the motion of the spacecraft in the gravitational field of these celestial bodies. The total ΔV of the maneuvers is minimized by eliminating the ΔV required for capturing the spacecraft by Moon. In this regard, only a mid-course impulsive maneuver is utilized for Moon ballistic capture. To achieve such trajectories, the optimization problem is parameterized with respect to the orbital elements of the ballistic capture orbits around Moon, the arrival date and a mid-course maneuver time. The equations of motion are solved backward in time with three impulsive maneuvers up to a specified low Earth parking orbit. The results show high potential and capability of this type of parameterization in finding several Pareto-optimal trajectories. Using the non-dominated sorting genetic algorithm with crowding distance sorting (NSGA-II) for the resulting multiobjective optimization problem, several trajectories are discovered. The resulting trajectories of the presented scheme permit alternative trade-off studies by designers incorporating higher level information and mission priorities.     

Analysis of the orbit lifetime of CubeSats in low Earth orbits including periodic variation in drag due to attitude motion

It is estimated that more than 22,300 human-made objects are in orbit around the Earth, with a total mass above 8,400,000 kg. Around 89% of these objects are non-operational and without control, which makes them to be considered orbital debris. These numbers consider only objects with dimensions larger than 10 cm. Besides those numbers, there are also about 2000 operational satellites in orbit nowadays. The space debris represents a hazard to operational satellites and to the space operations. A major concern is that this number is growing, due to new launches and particles generated by collisions. Another important point is that the development of CubeSats has increased exponentially in the last years, increasing the number of objects in space, mainly in the Low Earth Orbits (LEO). Due to the short operational time, CubeSats boost the debris population. One of the requirements for space debris mitigation in LEO is the limitation of the orbital lifetime of the satellites, which needs to be lower than 25 years. However, there are space debris with longer estimated decay time. In LEÓs, the influence of the atmospheric drag is the main orbital perturbation, and is used in maneuvers to increment the losses in the satellite orbital energy, to locate satellites in constellations and to accelerate the decay.The goal of the present research is to study the influence of aerodynamic rotational maneuver in the CubeSat?s orbital lifetime. The rotational axis is orthogonal to the orbital plane of the CubeSat, which generates variations in the ballistic coefficient along the trajectory. The maneuver is proposed to accelerate the decay and to mitigate orbital debris generated by non-operational CubeSats. The panel method is selected to determine the drag coefficient as a function of the flow incident angle and the spinning rate. The pressure distribution is integrated from the satellite faces at hypersonic rarefied flow to calculate the drag coefficient. The mathematical model considers the gravitational potential of the Earth and the deceleration due to drag. To analyze the effects of the rotation during the decay, multiple trajectories were propagated, comparing the results obtained assuming a constant drag coefficient with trajectories where the drag coefficient changes periodically. The initial perigees selected were lower than 400 km of altitude with eccentricities ranging from 0.00 to 0.02. Six values for the angular velocity were applied in the maneuver. The technique of rotating the spacecraft is an interesting solution to increase the orbit decay of a CubeSat without implementing additional de-orbit devices. Significant changes in the decay time are presented due to the increase of the mean drag coefficient calculated by the panel method, when the maneuver is applied, reducing the orbital lifetime, however the results are independent of the angular velocity of the satellite.     

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.     

椭圆轨道编队构型的初始化控制研究

研究了椭圆参考轨道下的编队构型初始化问题,给出轨道面内绕飞构型的三种初始化策略.基于状态转移矩阵研究了构型初始化的两脉冲机动方法;基于高斯摄动方程,以降低燃料消耗为原则,研究了构型初始化的三迹向脉冲和四迹向脉冲方法,推导了速度增量的求解方程;对上述三种构型初始化策略进行了仿真分析,利用序列二次规划和遗传算法对计算结果进行优化.数值仿真表明,三种初始化策略均能完成构型初始化,三迹向脉冲初始化比两脉冲初始化节省约41%的速度增量,四迹向脉冲初始化需要的速度增量与三迹向脉冲初始化基本相等.      

Dynamic orbit propagation in a gravitational field of an inhomogeneous attractive body using the Lagrange coefficients

The analytical methods have nearly been replaced by the numerical methods due to their higher accuracy and accessibility of computation facilities. The semi-analytical Lagrange method of orbit propagation using f and g series is a competitive alternative to the numerical integration technique if the Lagrange coefficients are derived in a full gravitational field. In this paper, a generalization of the Lagrange method of orbit propagation is introduced. In other words, we introduce a complete form of the Lagrange coefficients in all force fields developed in the spherical harmonics for example full gravitational field of the Earth. The method is numerically compared with the numerical integration technique. In order to show the numerical performance of the method, it has been implemented for orbit propagation of a GPS-like MEO and CHAMP-like LEO satellites. Discrepancy at centimeter level for CHAMP-like and sub-millimeter accuracy for GPS-like satellites shows relatively high performance of the developed algorithm. Compared to integration method, the proposed Lagrange method is nearly faster by a factor two for small Nmax and four for large Nmax.     

对地观测卫星有限推力快速轨道机动优化方法

研究了对地观测卫星在有限推力作用下,快速轨道机动到目标区域上空的优化算法.分析了任务需求并选取两次点火的变轨方法,使算法更具普适性.推导出考虑J2摄动的三维空间的有限推力动力学方程.使用Lambert方法确定有限推力时间变量的取值范围,并在此基础上使用遗传算法进行优化,最终得到时间最优意义的优化结果.通过一个算例对整个优化算法做了验证.结果表明该算法精度高、计算时间短,适用于快速轨道机动任务.      

Is the Earth’s orbital motion linked to the spin rotation of the Sun?

Time variations of the magnetic field of the Sun, seen as a star (the data 1968–2018, with more than 27 thousand daily measurements of the solar mean magnetic field), allowed to specify the rotation period of the gravitating solar mass: 27.027(6)?days, synodic. This indicates a presumably unknown physical connection between motions of the Sun and the Earth: in the course of a year our star accomplishes nearly 27 half-revolutions, while the planet itself performs an identical number of its spinnings during one complete axial revolution of the Sun. True origin of this strange Sun–Earth resonance is unknown, but it is supposed the phenomenon might be caused by slight coherent perturbations of gravity within the solar system.     

Properties of transit trajectory in the restricted three and four-body problem

In the restricted three-body problem if the Jacobi constant is just below the value corresponding to Lagrangian point only a little neck exists around the equilibrium point and capture trajectories are indicated as low-energy. Capture properties depend on the dynamics around these critical points and qualitative results can be obtained using linearized systems. In this paper, to study transit trajectory properties in the restricted three and four-body problem, the Earth–Moon–Sun–Satellite system is considered as example and studied using different models. In the circular restricted three-body problem (Earth–Moon–Satellite), transit, non transit and asymptotic trajectories, are easily identified by using the principal reference frame. Dynamics around Lagrangian point are then studied introducing the Moon eccentricity into the elliptical restricted three-body model. A preferential region for transit orbit is individuated and studied as a function of eigenvalue properties. To introduce the Sun effect, the bi-circular four-body model is considered and dynamics around Lagrangian points studied as a function of angular distance between Earth–Sun and Earth–Moon line. Finally, results obtained in the elliptical three-body model and bi-circular four-body model, are compared with numerical simulations using real Sun–Moon–Earth ephemeris.     

Halo轨道探测器的姿态描述与建模

日地系Halo轨道位于平动点附近,且不存在以地球为中心的轨道平面,这使得Halo轨道探测器的姿态描述和建模问题完全不同于近地卫星.针对三轴稳定方式的探测器,给出了2种姿态描述方案:对日定向方案(方案I,基于日、地敏感器)和旋转坐标系下定向方案(方案II,基于星敏感器).依据两者所定义的不同轨道坐标系,分别建立了姿态动力学模型.分析2类描述方案的相对误差,发现方案I的上限不超过5%,而方案II的上限不超过0.005%.并就方案I设计了姿态镇定控制器,以验证所提出的姿态描述方案的合理性.结果表明:具体工程应用中应以方案II为正常工作模式,并定期修正星历表的误差;方案I作为备份,保证探测器的正常运行.      

A simulation of Martian gravity field recovery by using a near equatorial orbiter

An improvement to the Martian gravity field may be achieved by means of future orbiting spacecraft with small eccentricity and low altitude exemplified through a newly proposed mission design that may be tested in upcoming reconnaissance of Mars. Here, the near equatorial orbital character (with an inclination approximating 10°, eccentricity as 0.01 and semi-major axis as 4000 km) is considered, and its contribution to Martian gravity field solution is analyzed by comparing it with a hypothetical polar circular orbiter. The solution models are evaluated in terms of the following viewpoints: power spectra of gravity field coefficients, correlations of low degree zonal coefficients, precise orbit determination, and error distribution of both Mars free air gravity anomaly and areoid. At the same time, the contributions of the near equatorial orbiters in low degree zonal coefficients time variations are also considered. The present results show that the near equatorial orbiter allows us to improve the accuracy of the Martian gravity field solution, decrease correlation of low degree zonal coefficients, retrieve much better time variable information of low degree zonal coefficients, improve precise orbit determination, and provide more accurate Mars free air gravity anomaly and areoid around the equatorial region.     

Low-energy near Earth asteroid capture using Earth flybys and aerobraking

Since the Sun-Earth libration points L₁ and L₂ are regarded as ideal locations for space science missions and candidate gateways for future crewed interplanetary missions, capturing near-Earth asteroids (NEAs) around the Sun-Earth L₁/L₂ points has generated significant interest. Therefore, this paper proposes the concept of coupling together a flyby of the Earth and then capturing small NEAs onto Sun–Earth L₁/L₂ periodic orbits. In this capture strategy, the Sun-Earth circular restricted three-body problem (CRTBP) is used to calculate target Lypaunov orbits and their invariant manifolds. A periapsis map is then employed to determine the required perigee of the Earth flyby. Moreover, depending on the perigee distance of the flyby, Earth flybys with and without aerobraking are investigated to design a transfer trajectory capturing a small NEA from its initial orbit to the stable manifolds associated with Sun-Earth L₁/L₂ periodic orbits. Finally, a global optimization is carried out, based on a detailed design procedure for NEA capture using an Earth flyby. Results show that the NEA capture strategies using an Earth flyby with and without aerobraking both have the potential to be of lower cost in terms of energy requirements than a direct NEA capture strategy without the Earth flyby. Moreover, NEA capture with an Earth flyby also has the potential for a shorter flight time compared to the NEA capture strategy without the Earth flyby.     

On the lunar node resonance of the orbital plane evolution of the Earth’s satellite orbits

This paper aims to investigate the effects of lunar node resonance on the circular medium Earth orbits (MEO). The dynamical model is established in classical Hamiltonian systems with the application of Lie transform to remove the non-resonant terms. Resonant condition, stability and phase structures are studied. The lunar node resonance occurs when the secular changing rates of the orbital node (with respect to the equator) and the lunar node (with respect to the ecliptic) form a simple integer ratio. The resonant conditions are satisfied for both inclined and equatorial orbits. The orbital plane would have long period (with typical timescales of several centuries) fluctuation due to the resonance.     

Yarkovsky-Schach effect on space debris motion

The Yarkovsky-Schach effect is a small perturbation affecting Earth satellites and space debris illuminated by the Sun. It was first applied to the orbit of LAGEOS satellites as an explanation of the residuals in orbital elements. In this work, we carry out several numerical integration tests taking into consideration various orbit and rotation parameters, in order to analyse this effect in a broader context. The semi-major axis variations remain small and depend on the spin axis attitude with respect to the Sun. We show that the force amplitude is maximised for orbits inclined with $i\approx$?20–30°. We also observe the influence on other orbital elements, notably on the orbit inclination. However, these effects are clearly observed only on long timescales; in our simulations, we propagated the orbits for 200?y. The Yarkovsky-Schach effect is thus confirmed to have a minuscule magnitude. It should be taken into account in studies requiring high-precision orbit determination, or on expanded timescales.     

大偏心率轨道星上快速计算方法

卫星星上轨道的计算方法通常有解析方法和数值积分方法等.采用数值积分方法计算量较大,对于大偏心率的探测器任务轨道,解析方法也不适用.通过引入精密轨道与二体轨道的位置差,利用位置差拟合切比雪夫多项式系数,大大提高拟合精度.该方法计算精度高、计算速度快,可以很好地解决大偏心率椭圆轨道的星上轨道计算问题.