An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

Periodic attitude motions along planar orbits in the elliptic restricted three-body problem

A way to improve the accuracy of the three-body problem model is taking into account the eccentricity of primary attractors. Elliptic Restricted Three-Body Problem (ER3BP) is a model for studying spacecraft trajectory within the three-body problem such that the orbital eccentricity of primaries is reflected in it. As the principal cause of perturbation in the employed dynamical model, the primaries eccentricity changes the structure of orbits compared to the ideal Circular Restricted Three-Body Problem (CR3BP). It also changes the attitude behavior of a spacecraft revolving along periodic orbits in this regime. In this paper, the coupled orbit-attitude dynamics of a spacecraft in the ER3BP are exploited to find precise periodic solutions as the spacecraft is considered to be in planar orbits around Lagrangian points and Distant Retrograde Orbits (DRO). Periodic solutions are repetitious behaviors in which spacecraft whole dynamics are repeated periodically, these periodic behaviors are the main interest of this study because they are beneficial for future mission designs and allow delineation of the system’s governing dynamics. Previous studies laid the foundation for spacecraft stability analysis or studying pitch motion of spacecraft in the ER3BP regime. While in this paper, at first, initial guesses for correction algorithms were derived through verified search methods, then correction algorithms were used to refine calculated orbit-attitude periodic behaviors. Periodic orbits and full periodic solutions are portrayed and compared to previous studies and simpler models. Natural periodic solutions are valuable information eventuate in the longer functional lifetime of spacecraft. Since the problem assumption considered in this paper is much closer to real mission conditions, these results may be the means to use natural bounded motions in the actual operational environment.     

椭圆轨道编队的构形变化控制方法

推导了椭圆轨道编队的一些典型构形模态,指出椭圆轨道卫星编队的模态不同于圆轨道的构形模态,是极为丰富的。基于T-H方程的解析解,推导了副星在沿航向常值推力作用下的椭圆轨道编队的构形变化控制方法。推导的构形变化控制方法,通过合理的选择控制量的作用时刻,可达到大量节约燃料的目的。最后,给出仿真算例。     

Formation flying for magnetic sails around artificial equilibrium points

This paper investigates the problem of magnetic sail-based spacecraft formation control around the artificial equilibrium points (AEPs), which can eliminate the requirement of the propellant. The thrusts are achieved by utilizing the interaction between the solar wind and the artificial magnetosphere generated by superconducting current coil onboard. The circular restricted three-body problem (CRTBP) of magnetic sail is discussed including the allowed regions and linear stability of AEPs, the locations of collinear AEPs and the possibility of existence of periodic orbits around the collinear AEPs. Next, the dynamical models of magnetic sail formation around the collinear AEPs are established. A novel fast fixed-time nonsingular terminal sliding mode controller (FFNTSM) based on fixed-time disturbance observer (FTDO) is developed to account for external disturbances. Several numerical simulations are conducted to substantiate that spacecraft formation can be precisely controlled by the proposed propellantless propulsion method in the presence of external disturbances.     

基于目标特性的航天器太阳光压面积计算方法

太阳光压是影响深空探测航天器轨道确定与预报精度最主要的摄动力.针对实际任务需求,采用了一种基于目标特性的光压面积建模与计算方法,根据航天器形状、尺寸、表面材料以及材料光学特性等信息,实现了分析型光压模型的建立与求解,提高了计算效率和精度,可快速计算目标在光照方向上的光压面积、投影面积以及光压比例因子等参数.通过长方体光压面积理论值与仿真值的对比,验证了该方法的准确性和有效性.针对复杂结构探测器开展了光压面积计算,可为深空探测航天器精密定轨中的光压模型解算、定轨及预报提供参考.      

Bound orbits of solar sails and general relativity

We study general relativistic effects on the bound orbits of solar sails. The combined effects of spacetime curvature and solar radiation pressure (SRP) lead to deviations from Kepler’s third law. Such kind of deviations also arise from frame dragging, the gravitational multipole moments of the sun, a net electric charge on the sun, and a positive cosmological constant. The SRP increases these deviations by several orders of magnitude, possibly rendering some of them detectable. We consider how the SRP modifies the perihelion shift of non-circular orbits, as well as the Lense-Thirring effect involving the precession of polar orbits. We investigate how the pitch angle for non-Keplerian orbits changes due to the partial absorption of light, general relativistic effects, and the oblateness of the sun. It is predicted that there is an analog of the Lense-Thirring effect for non-Keplerian orbits, in that the orbital plane precesses around the sun. We also consider the Poynting–Robertson effect and show that this effect can, in principle, be compensated for by an extremely small tilt of the solar sail.     

Dynamics and control of spacecraft hovering using the geomagnetic Lorentz force

To achieve hovering, a spacecraft thrusts continuously to induce an equilibrium state at a desired position. Due to the constraints on the quantity of propellant onboard, long-time hovering around low-Earth orbits (LEO) is hardly achievable using traditional chemical propulsion. The Lorentz force, acting on an electrostatically charged spacecraft as it moves through a planetary magnetic field, provides a new propellantless method for orbital maneuvers. This paper investigates the feasibility of using the induced Lorentz force as an auxiliary means of propulsion for spacecraft hovering. Assuming that the Earth’s magnetic field is a dipole that rotates with the Earth, a dynamical model that characterizes the relative motion of Lorentz spacecraft is derived to analyze the required open-loop control acceleration for hovering. Based on this dynamical model, we first present the hovering configurations that could achieve propellantless hovering and the corresponding required specific charge of a Lorentz spacecraft. For other configurations, optimal open-loop control laws that minimize the control energy consumption are designed. Likewise, the optimal trajectories of required specific charge and control acceleration are both presented. The effect of orbital inclination on the expenditure of control energy is also analyzed. Further, we also develop a closed-loop control approach for propellantless hovering. Numerical results prove the validity of proposed control methods for hovering and show that hovering around low-Earth orbits would be achievable if the required specific charge of a Lorentz spacecraft becomes feasible in the future. Typically, hovering radially several kilometers above a target in LEO requires specific charges on the order of 0.1 C/kg.     

Mapping orbits around the asteroid 2001SN263

The present paper has the goal of mapping orbits, with respect to the perturbations, for a spacecraft traveling around the asteroid 2001SN₂₆₃. This asteroid is a triple system, which center of mass is in an elliptic orbit around the Sun. The perturbations considered in the present model are the ones due to the oblateness of the central body, the gravity field of the two satellite bodies (Beta and Gamma), the Sun, the Moon, the asteroids Vesta, Pallas and Ceres and all the planets of the Solar System. This mapping is important, because it shows the relative importance of each force for a given orbit for the spacecraft, helping to make a decision about which forces need to be included in the model for a given accuracy and nominal orbit. Another important application of this type of mapping is to find orbits that are less perturbed, since it is expected that those orbits have good potential to require a smaller number of station-keeping maneuvers. Simulations under different conditions are made to find those orbits. The main reason to study those trajectories is that, currently, there are several institutions in Brazil studying the possibility to make a mission to send a spacecraft to this asteroid (the so-called ASTER mission), because there are many important scientific studies that can be performed in that system. The results showed that Gamma is the main perturbing body, followed by Beta (10 times smaller) and the group Sun–Mars-oblateness of Alpha, with perturbations 1000 times weaker than the effects of Gamma. The other bodies have perturbations 10⁷ times smaller. The results also showed that circular and polar orbits are less perturbed, when compared to elliptical and equatorial orbits. Regarding the semi-major axis, an internal orbit is the best choice, followed by a larger external orbit. The inclination of the orbit plays an important role, and there are values for the inclination where the perturbations show minimum and maximum values, so it is important to make a good decision on those values.     

Spacecraft fuel-optimal and balancing maneuvers for a class of formation reconfiguration problems

This paper presents fuel optimal and balancing methodologies for reconfiguring multiple spacecraft in formation subject to a Newtonian gravity field. For a kind of continuous-thrust propulsion system, a fuel-optimal control problem is formulated to minimize the integral squared control subject to the linearized Hill or Clohessy–Wiltshire dynamics of relative motion with respect to a circular reference orbit. Palmer’s analytical solution for general reconfiguration is adapted to maneuvers between projected circular orbits, resulting in the optimal fuel consumption index as a function of configuration parameters such as orbit radius, phase angle, and transfer time. Parametric analyses reveal unique characteristics of individual fuel optimality and gross fuel consumption: for an arbitrary selection of initial/terminal orbit radii, (i) there exist special transfer times such that individual fuel consumption is optimally uniform for all phase angles, and (ii) the total fuel expenditure for a group of three or more spacecraft is invariant for the relatively same configuration with respect to the departure phase. These results serve to effectively design fuel balancing strategies for formation reconfiguration of multiple spacecraft.     

考虑太阳摄动的小行星附近轨道动力学

本文研究了艳后星(216 Kleopatra)和爱神星(433 Eros)附近的周期轨道,在考虑太阳引力摄动的情况下,发现了以往所遗漏的216 Kleopatra轨道族和环绕433 Eros的12族周期轨道,并且给出了它们的特性。研究结果表明,太阳引力对小行星平衡点位置的影响很小,但是对平衡点上航天器运动的影响较大。同族不稳定轨道中,大Jacobi常数轨道更容易在摄动后保持轨道原来特性,这很好地解释了小行星卫星在较远轨道上长期存在的可能性。     

近地轨道集群航天器电磁编队飞行非线性反馈控制方法

针对近地轨道集群航天器电磁编队飞行的动力学和控制问题, 提出了一种非线性反馈控制方法. 基于电磁力模型和地磁场模型, 分析了地磁场对近地轨道电磁编队的影响; 建立了集群航天器电磁编队高精度相对轨道动力学模型; 基于Lyapunov稳定性理论设计了一种非线性反馈控制律, 利用该方法对两星电磁编队维持控制进行了仿真验证. 仿真结果表明, 地磁场引起的电磁干扰力可以忽略, 但是电磁干扰力矩的影响必须考虑; 近地轨道集群航天器电磁编队是可控的, 所设计的控制方法是可行的.      

Design of a formation of solar pumped lasers for asteroid deflection

This paper presents the design of a multi-spacecraft system for the deflection of asteroids. Each spacecraft is equipped with a fibre laser and a solar concentrator. The laser induces the sublimation of a portion of the surface of the asteroid, and the resultant jet of gas and debris thrusts the asteroid off its natural course. The main idea is to have a formation of spacecraft flying in the proximity of the asteroid with all the spacecraft beaming to the same location to achieve the required deflection thrust. The paper presents the design of the formation orbits and the multi-objective optimisation of the formation in order to minimise the total mass in space and maximise the deflection of the asteroid. The paper demonstrates how significant deflections can be obtained with relatively small sized, easy-to-control spacecraft.     

On the stabilizing effect of Solar Radiation Pressure in the Earth-Moon system

Solar sails change the natural dynamics of systems: Trajectories that are driven by gravitational forces can be displaced and changed because of the effect of Solar Radiation Pressure (SRP). Moreover, if the lightness number of the sail is large enough, the instability of certain orbits can be diminished and even removed. In this paper we modify two models for the motion of a probe in the Earth-Moon system that include the effect of Sun’s gravity to take also into account the effect of SRP. These models, the Bicircular Problem (BCP) and the Quasi-Bicircular Problem (QBCP), are periodic perturbations of the Earth-Moon Restricted Three Body Problem (RTBP). The models are modified to consider the effect of the SRP upon a solar sail. We provide examples of periodic orbits that are stabilized (or made less unstable) due to the effect of SRP.     

Tetrahedral satellite formation: Geomagnetic measurements exchange and interpolation

This note presents a study of a four-satellite tetrahedral formation to collect, process, and exchange multipoint measurements of geomagnetic field in a near-polar orbit. The study is conducted as a series of numerical experiments based on simulated spacecraft orbits and corresponding geomagnetic field models output. The four satellites are assumed to move in near-circular orbits specifically chosen to maintain the tetrahedron quality. The satellites exchange their simulated magnetometers readings and the collected multipoint measurements are processed on board of any of them thus creating an instantaneous interpolated map of the geomagnetic field in the interior of the tetrahedron. Interpolation is carried out with the use of Kriging algorithms, known in geostatistics for capturing spatial correlation of the data and taking into account statistical properties of the interpolated variables. We propose a concept of a servicing formation, and analyze interpolation accuracy for different formation sizes. It is then discussed how the processed multipoint measurements can be provided as a service to other nearby satellites. Finally, we show that using the existing COTS magnetometers it is possible to obtain real-time interpolation data, which are more precise at a given point and time than a conventional onboard magnetic field model, thus ensuring better attitude determination routines performance in the serviced spacecraft.     

Improving magnetosphere in situ observations using solar sails

Past and current magnetosphere missions employ conventional spacecraft formations for in situ observations of the geomagnetic tail. Conventional spacecraft flying in inertially fixed Keplerian orbits are only aligned with the geomagnetic tail once per year, since the geomagnetic tail is always aligned with the Earth-Sun line, and therefore, rotates annually. Solar sails are able to artificially create sun-synchronous orbits such that the orbit apse line remains aligned with the geomagnetic tail line throughout the entire year. This continuous presence in the geomagnetic tail can significantly increase the science phase for magnetosphere missions. In this paper, the problem of solar sail formation design is explored using nonlinear programming to design optimal two-craft, triangle, and tetrahedron solar sail formations, in terms of formation quality and formation stability. The designed formations are directly compared to the formations used in NASA’s Magnetospheric Multi-Scale mission.     

Spin motion determination of the Envisat satellite through laser ranging measurements from a single pass measured by a single station

The Satellite Laser Ranging (SLR) technology is used to accurately determine the position of space objects equipped with so-called retro-reflectors or retro-reflector arrays (RRA). This type of measurement allows to measure the range to the spacecraft with high precision, which leads to determination of very accurate orbits for these targets. Non-active spacecraft, which are not attitude controlled any longer, tend to start to spin or tumble under influence of the external and internal torques and forces.If the return signal is measured for a non-spherical non-active rotating object, the signal in the range residuals with respect to the reference orbit is more complex. For rotating objects the return signal shows an oscillating pattern or patterns caused by the RRA moving around the satellite’s centre of mass. This behaviour is projected onto the radial component measured by the SLR.In our work, we demonstrate how the SLR ranging technique from one sensor to a satellite equipped with a RRA can be used to precisely determine its spin motion during one passage. Multiple SLR measurements of one target over time allow to accurately monitor spin motion changes which can be further used for attitude predictions. We show our solutions of the spin motion determined for the non-active ESA satellite Envisat obtained from measurements acquired during years 2013–2015 by the Zimmerwald SLR station, Switzerland. All the necessary parameters are defined for our own so-called point-like model which describes the motion of a point in space around the satellite centre of mass.     

Attitude stability and periodic attitudes of rigid spacecrafts on the stationary orbits around asteroid 216 Kleopatra

In this work, equilibrium attitude configurations, attitude stability and periodic attitude families are investigated for rigid spacecrafts moving on stationary orbits around asteroid 216 Kleopatra. The polyhedral approach is adopted to formulate the equations of rotational motion. In this dynamical model, six equilibrium attitude configurations with non-zero Euler angles are identified for a spacecraft moving on each stationary orbit. Then the linearized equations of attitude motion at equilibrium attitudes are derived. Based on the linear system, the necessary conditions of stability of equilibrium attitudes are provided, and stability domains on the spacecraft’s characteristic plane are obtained. It is found that the stability domains are distributed in the first and third quadrants of the characteristic plane and the stability domain in the third quadrant is separated into two regions by an unstable belt. Subsequently, we present the linear solution around a stable equilibrium attitude point, indicating that there are three types of elemental periodic attitudes. By means of numerical approaches, three fundamental families of periodic solutions are determined in the full attitude model.     

航天器相对运动与Hill解的适用性分析

研究了空间飞行器编队中最具基础性的问题之一 ,即相对运动的解析表达及 Hill方程的适用条件。通过建立相对运动的通解公式 ,针对不同性质的初值深入地分析了其相对运动轨迹的本质特征 ,并给出了 Hill方程的适用条件。此外 ,文中还给出了一个新的编队设计简化公式     

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.