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.     

Evolution of the orbital elements for objects with high area-to-mass ratios in geostationary transfer orbits

A new population of uncatalogued objects in geosynchronous Earth orbits (GEO), with a mean motion of about 1 rev/day and eccentricities up to 0.6, has been identified recently. The first observations of this new type of objects were acquired in the framework of the European Space Agency’s (ESA) search for space debris in GEO and the geostationary transfer orbit (GTO) using the ESA 1-m telescope on Tenerife. Earlier studies have postulated that the perturbations due to the solar radiation pressure can lead to such large eccentricities for GEO objects with a high area-to-mass ratio (A/M). The simulations showed that the eccentricities of GEO objects with large A/M exhibit periodic variations with periods of about one year and amplitudes depending on the value of A/M. The findings of these studies could be confirmed by observations from the ESA 1-m telescope on Tenerife.     

The classical Laplace plane as a stable disposal orbit for geostationary satellites

The classical Laplace plane is a frozen orbit, or equilibrium solution for the averaged dynamics arising from Earth oblateness and lunisolar gravitational perturbations. The pole of the orbital plane of uncontrolled GEO satellites regress around the pole of the Laplace plane at nearly constant inclination and rate. In accordance with Friesen et al. (1993), we show how this stable plane can be used as a robust long-term disposal orbit. The current graveyard regions for end-of-life retirement of GEO payloads, which is several hundred kilometers above GEO depending on the spacecraft characteristics, cannot contain the newly discovered high area-to-mass ratio debris population. Such objects are highly susceptible to the effects of solar radiation pressure exhibiting dramatic variations in eccentricity and inclination over short periods of time. The Laplace plane graveyard, on the contrary, would trap this debris and would not allow these objects to rain down through GEO. Since placing a satellite in this inclined orbit can be expensive, we discuss some alternative disposal schemes that have acceptable cost-to-benefit ratios.     

Apparent rotation properties of space debris extracted from photometric measurements

Knowledge about the rotation properties of space debris objects is essential for the active debris removal missions, accurate re-entry predictions and to investigate the long-term effects of the space environment on the attitude motion change. Different orbital regions and object’s physical properties lead to different attitude states and their change over time.Since 2007 the Astronomical Institute of the University of Bern (AIUB) performs photometric measurements of space debris objects. To June 2016 almost 2000 light curves of more than 400 individual objects have been acquired and processed. These objects are situated in all orbital regions, from low Earth orbit (LEO), via global navigation systems orbits and high eccentricity orbit (HEO), to geosynchronous Earth orbit (GEO). All types of objects were observed including the non-functional spacecraft, rocket bodies, fragmentation debris and uncorrelated objects discovered during dedicated surveys. For data acquisition, we used the 1-meter Zimmerwald Laser and Astrometry Telescope (ZIMLAT) at the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald, Switzerland. We applied our own method of phase-diagram reconstruction to extract the apparent rotation period from the light curve. Presented is the AIUB’s light curve database and the obtained rotation properties of space debris as a function of object type and orbit.     

Dynamical evolution of high area-to-mass ratio objects in Molniya orbits

For special demands, some notable orbit types have been developed by human, including the Molniya orbits, which have a relatively high eccentricity up to about 0.7, and a period of 12 h. Considering that space debris with high area-to-mass ratio (A/M) has been discovered, such objects may also exist in Molniya orbits due to spacecraft and upper stages fragmentation events. However, there are not sufficient studies of the complex dynamical phenomena of such orbits. These studies can enrich the knowledge about the long-term evolution of these orbits, be helpful to propose uncatalogued objects observation and identification, and also set the protected region as well as active debris removal. In this paper, the characteristics of 2:1 resonance of Molniya satellite orbits are studied. A large set of numerical simulations, including all the relevant perturbations, is carried out to further investigate the main characteristics, and special attention is payed to the dynamical evolution of objects with high A/M, particularly affected by the direct solar radiation pressure. The long-term dynamical evolution of orbital elements, as well as the dependency of lifetime on the A/M value, is discussed.     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.     

基于主动约束层阻尼控制的大型挠性航天器动力学模型

以大型挠性航天器简化模型——中心刚体和固接板混合系统为对象,同时考虑板表面粘贴的主动约束层阻尼结构,根据拟坐标拉格朗日方程,建立了包含所有高阶小量的完整动力学模型。建模时,采用混合坐标描述中心刚体和主动约束阻尼板的运动,用Golla-Hughes-Mctavish模型描述阻尼材料的复模量。基于该动力学模型可进行空间柔性板振动的主动约束阻尼控制器设计和大型挠性航天器的姿态动力学研究。     

A deorbiter CubeSat for active orbital debris removal

This paper introduces a mission concept for active removal of orbital debris based on the utilization of the CubeSat form factor. The CubeSat is deployed from a carrier spacecraft, known as a mothership, and is equipped with orbital and attitude control actuators to attach to the target debris, stabilize its attitude, and subsequently move the debris to a lower orbit where atmospheric drag is high enough for the bodies to burn up. The mass and orbit altitude of debris objects that are within the realms of the CubeSat’s propulsion capabilities are identified. The attitude control schemes for the detumbling and deorbiting phases of the mission are specified. The objective of the deorbiting maneuver is to decrease the semi-major axis of the debris orbit, at the fastest rate, from its initial value to a final value of about 6471?km (i.e., 100?km above Earth considering a circular orbit) via a continuous low-thrust orbital transfer. Two case studies are investigated to verify the performance of the deorbiter CubeSat during the detumbling and deorbiting phases of the mission. The baseline target debris used in the study are the decommissioned KOMPSAT-1 satellite and the Pegasus rocket body. The results show that the deorbiting times for the target debris are reduced significantly, from several decades to one or two years.     

万有引力场中带挠性轴太阳帆板航天器的姿态稳定性

研究了万有引力场中带挠性轴太阳帆板航天器的姿态运动、导出了航天器系统的广义势能,利用Liapunov直接方法判断带挠性轴太阳帆板航天器在轨道坐标系内相对平衡的稳定性,得到航天器姿态稳定性的充分条件。讨论了挠性轴扭转刚度及太阳帆板质量几何等因素对航天器姿态稳定性的影响。      

大面质比空间碎片在太阳光压和引力作用下的轨道演化

分析了较高轨道(a > 10000km)大面质比空间碎片的轨道动力学演化问题. 重点讨论了位于地球同步轨道的空间碎片轨道演化问题, 并给出轨道偏心率 随时间演化的表达式. 通过进一步分析得出, 倾角大于63°26'的GTO轨 道空间碎片, 仅在J₂和第三体摄动影响下, 会出现轨道偏心率升高; 而对 于大面质比空间碎片, 在J₂项和太阳光压同时作用下, 当近地点指向的角 变率与太阳平黄经变化率接近时, 会出现长期共振现象, 导致轨道偏心率升 高, 近地点降低. 分析还得出, 轨道演化过程中, 偏心率的最大值与初始轨 道近地点的指向有关.      

Properties of the high area-to-mass ratio space debris population at high altitudes

In the framework of its space debris research activities ESA established an optical survey program to study the space debris environment at high altitudes, in particular in the geostationary ring and in the geostationary transfer orbit region. The Astronomical Institute of the University of Bern (AIUB) performs these surveys on behalf of ESA using ESA’s 1-m telescope in Tenerife. Regular observations were started in 1999 and are continued during about 120–140 nights per year. Results from these surveys revealed a substantial amount of space debris at high altitudes in the size range from 0.1 to 1 m. Several space debris populations with different dynamical properties were identified in the geostationary ring. During the searches for debris in the geostationary transfer orbit region a new population of objects in unexpected orbits, where no potential progenitors exist, was found. The orbital periods of these objects are clustered around one revolution per day; the eccentricities, however, are scattered between 0 and 0.6. By following-up some of these objects using the ESA telescope and AIUB’s 1-m telescope in Zimmerwald, Switzerland, it was possible to study the properties of this new population. One spectacular finding from monitoring the orbits over time spans of days to months is the fact that these objects must have extreme area-to-mass ratios, which are by several orders of magnitudes higher than for ‘normal-type’ debris. This in turn supports the hypothesis that the new population actually is debris generated in or near the geostationary ring and which is in orbits with periodically varying eccentricity and inclination due to perturbations by solar radiation pressure. In order to further study the nature of these debris, multi-color and temporal photometry (light curves) were acquired with the Zimmerwald telescope. The light curves show strong variations over short time intervals, including signals typical for specular reflections. Some objects exhibit distinct periodic variations with periods ranging from 10 to several 100 s. All this is indicative for objects with complicated shapes and some highly reflective surfaces.     

基于压电纤维复合材料的航天器动力学建模与振动抑制

压电纤维复合材料(MFC)在柔性航天器的振动主动抑制中具有很好的应用前景。利用哈密顿原理和压电驱动的载荷比拟方法,建立了带MFC压电驱动的离散形式的刚柔耦合动力学方程,采用线性二次型最优控制(LQR)算法进行主动控制。结果表明：在航天器的柔性体受到脉冲载荷激励条件下,使用MFC驱动器可以实现航天器挠性振动的快速抑制,并且同时保持中心刚体姿态的稳定性,即能够实现挠性振动与姿态运动的协同控制。基于MFC的主动控制方法对于高频响应也具有较好的控制效果。对于柔性占优的航天器,采用MFC的主动控制优于被动控制。本文方法在处理具有复杂柔性体的航天器时更具优势,更适合于工程应用。     

挠性卫星轨控期间动力学与姿态控制

卫星轨道控制期间,轨道控制推力会激振挠性太阳帆板,从而也影响卫星的姿态。用拉格朗日方程建立了带有大型太阳帆板的卫星动力学模型,分析了卫星质心运动、姿态运动与挠性振动的耦合关系。根据姿态控制推力器的输出特性,设计了轨道控制期间卫星姿态控制方案。通过数学仿真验证了轨道控制推力对挠性帆板与卫星姿态的影响,验证了轨控期间姿态控制方案的有效性。     

Shape optimization of a three-dimensional membrane-structured solar sail using an angular momentum unloading strategy

A shape of the satellite’s solar sail membrane is essential for unloading angular momentum in the three-axis stabilized attitude control system because the three-dimensional solar sail can receive solar radiation pressure from arbitrary directions. In this paper, the objective is the shape optimization of a three-dimensional membrane-structured solar sail using the angular momentum unloading strategy. We modelled and simulated the solar radiation pressure torque, for unloading angular momentum. Using the simulation system, since the unloading angular momentum rate is maximized, the shape of the three-dimensional solar sail was optimized using a Genetic algorithm and Sequential Quadratic Programming. The unloading velocity in the optimized shaped solar sail was greatly improved with respect to a conventional flat or pyramid solar sail.     

Deorbiter CubeSat mission design

This paper presents the mission design for a CubeSat-based active debris removal approach intended for transferring sizable debris objects from low-Earth orbit to a deorbit altitude of 100 km. The mission consists of a mothership spacecraft that carries and deploys several debris-removing nanosatellites, called Deorbiter CubeSats. Each Deorbiter is designed based on the utilization of an eight-unit CubeSat form factor and commercially-available components with significant flight heritage. The mothership spacecraft delivers Deorbiter CubeSats to the vicinity of a predetermined target debris, through performing a long-range rendezvous maneuver. Through a formation flying maneuver, the mothership then performs in-situ measurements of debris shape and orbital state. Upon release from the mothership, each Deorbiter CubeSat proceeds to performing a rendezvous and attachment maneuver with a debris object. Once attached to the debris, the CubeSat performs a detumbling maneuver, by which the residual angular momentum of the CubeSat-debris system is dumped using Deorbiter’s onboard reaction wheels. After stabilizing the attitude motion of the combined Deorbiter-debris system, the CubeSat proceeds to performing a deorbiting maneuver, i.e., reducing system’s altitude so much so that the bodies disintegrate and burn up due to atmospheric drag, typically at around 100 km above the Earth surface. The attitude and orbital maneuvers that are planned for the mission are described, both for the mothership and Deorbiter CubeSat. The performance of each spacecraft during their operations is investigated, using the actual performance specifications of the onboard components. The viability of the proposed debris removal approach is discussed in light of the results.     

太阳帆航天器编队维持和重构的方法研究

针对高面质比航天器可以利用太阳光压进行轨道控制的特点,本文提出一种太阳帆航天器编队构型维持和重构的方法.该方法通过控制主从航天器太阳帆姿态角和反射系数,调整主从航天器之间的光压差,产生抵消编队成员间相对运动受到摄动差或进行轨道机动时所需的连续小推力,从而实现编队构型的维持和重构.仿真结果表明,在主航天器太阳帆的姿态角和反射系数相对固定的条件下,对于太阳同步轨道上的高面质比太阳帆航天器编队,使用滑模控制方法,能够调整编队中从航天器太阳帆的姿态角和反射系数产生推力抵消摄动力影响,达到长期维持太阳帆航天器编队构型的目的;通过开环控制方法,能够调整编队中从航天器太阳帆的姿态角和反射系数产生连续小推力,在较长时间周期内实现编队重构.     

基于特征模型的挠性航天器姿态快速机动研究

刚体卫星的大角度姿态机动可以用常规的四元数反馈控制,当挠性帆板的振动和中心刚体的耦合系数很大时,大角度快速机动后姿态的稳定度较差．结合特征建模理论,设计一种卫星大角度机动的黄金分割控制算法,对三轴带挠性帆板的航天器姿态机动进行仿真,仿真结果表明,机动完成后的控制精度比四元数反馈控制方法的精度高一个数量级以上。     

Faster algorithm of debris cloud orbital character from spacecraft collision breakup

Space debris is polluting the space environment. Collision fragment is its important source. NASA standard breakup model, including size distributions, area-to-mass distributions, and delta velocity distributions, is a statistic experimental model used widely. The general algorithm based on the model is introduced. But this algorithm is difficult when debris quantity is more than hundreds or thousands. So a new faster algorithm for calculating debris cloud orbital lifetime and character from spacecraft collision breakup is presented first. For validating the faster algorithm, USA 193 satellite breakup event is simulated and compared with general algorithm. Contrast result indicates that calculation speed and efficiency of faster algorithm is very good. When debris size is in 0.01–0.05 m, the faster algorithm is almost a hundred times faster than general algorithm. And at the same time, its calculation precision is held well. The difference between corresponding orbital debris ratios from two algorithms is less than 1% generally.     

On-orbit verification of fuel-free attitude control system for spinning solar sail utilizing solar radiation pressure

This paper introduces a new attitude control system for a solar sail, which leverages solar radiation pressure. This novel system achieves completely fuel-free and oscillation-free attitude control of a flexible spinning solar sail. This system consists of thin-film-type devices that electrically control their optical parameters such as reflectivity to generate an imbalance in the solar radiation pressure applied to the edge of the sail. By using these devices, minute and continuous control torque can be applied to the sail to realize very stable and fuel-free attitude control of the large and flexible membrane. The control system was implemented as an optional attitude control system for small solar power sail demonstrator named IKAROS (Interplanetary Kite-craft Accelerated by Radiation Of the Sun). In-orbit attitude control experiments were conducted, and the performance of the controller was successfully verified in comparison with the ground-based analytical performance estimation.