基于星下点机动的再入飞行器离轨规划

针对约束再入点地理位置的再入飞行器离轨问题,提出了一种基于星下点机动的离轨规划方法。再入飞行器的离轨轨道设计受到飞行器当前轨道状态和再入点参数的约束。首先,基于轨道飞行原理,建立了一般椭圆轨道冲量模型下离轨制动参数和再入点参数的关系,分析了最优离轨的推力施加原则;其次,在考虑地球自转的前提下,提出了直接离轨必要条件,针对约束再入点经纬度的问题,完善了利用非线性规划优化方法确定有限推力模型下离轨点位置的策略,同时给出了符合燃料最优目标的离轨制动参数;最后,探讨了一般情况下初始轨道不满足直接离轨必要条件时,为满足星下点约束而进行的轨道机动施加策略。      

气动辅助异面变轨在多碎片清除中的应用分析

碎片清除飞行器异面变轨需要消耗大量燃料.从气动辅助异面变轨优化设计及被清除碎片轨道高度差值、倾角差值等参数对变轨性能的影响出发,比较分析了优化气动辅助异面变轨与双脉冲霍曼轨道转移的燃料节约量,研究了不同轨道高度差对于实施气动辅助变轨燃料节约量的影响.当地球静止轨道（GEO）与低地轨道（LEO）间气动辅助变轨优化速度增量约为1.55km·s-1、质量面积比172kg·m-2、比冲310s、轨道倾角变化16°时,燃料节约率约为45%.对比研究了不同轨道高度差LEO轨道间实施气动辅助变轨的燃料节约情况.结果表明:随着轨道高度的增加,气动辅助优化效率逐渐降低;在相同高度轨道间实施异面变轨,随着轨道倾角的增加,气动辅助变轨燃料节约率先增大后减小,倾角改变量约为20°时,燃料节约率最大;当轨道倾角为5°时,采用气动辅助变轨和双脉冲变轨的燃料消耗量相同.      

空间碎片碰撞预警研究

介绍了空间碎片的分布和危害，阐明了航天器进行碎片预警规避的必要性，讨论了碎片碰撞预警技术中的关键问题，用已有的碰撞和规避实例对自行研制的软件进行了计算检验，结果证明了软件的正确性，也说明了国内进行碎片预警工作的可行性．     

基于驻留时间比的近地空间碎片环境建模

为分析近地空间碎片的分布规律,提出了一种以碎片在空间网格内驻留时间为基础的碎片环境统计建模方法.该方法利用多项式拟合和求根方法统计碎片在空间网格内的停留时间,获取模型基础数据,并据此采用多项式预测、插值和时间序列分析等技术,综合分析空间碎片的分布与演化规律.给出了一个基于双行根数(TLE,Two Line Elements)数据的建模实例,该实例通过了ORDEM2000模型的对比验证,并获得了一些更精细的近地空间碎片环境特征.所得建模方法和分析结论可为长期运行的近地航天器轨道设计、碰撞风险评估及防护等提供技术支撑.     

针对空间碎片捕获的绕飞轨道设计

地球轨道上日益增长的碎片云已引起各个航天国家的担忧,地球轨道上可编目的空间物体数量30多年内增长了2倍多,若不实施主动清除,碎片的数量将在未来200年内快速增长,给空间系统的安全带来极大的威胁。自然椭圆绕飞轨道可在目标附近长时间绕飞,可保证碎片捕获系统具有长时间的观测、捕获时间。文章提出通过设计绕飞轨道来实现捕获碎片的方案,介绍并分别推导了基于C-W方程和轨道根数两种方式绕飞轨道设计的方法。针对假想的捕获目标,基于轨道根数方法设计了5种脉冲变轨的轨道方案,并进行了相应的轨道算例仿真。仿真结果表明:该方案可适用于任意的初始相位差,具备一定的工程实现意义。     

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.     

Instability of the present LEO satellite populations

Several studies conducted during 1991–2001 demonstrated, with some assumed launch rates, the future unintended growth potential of the Earth satellite population, resulting from random, accidental collisions among resident space objects. In some low Earth orbit (LEO) altitude regimes where the number density of satellites is above a critical spatial density, the production rate of new breakup debris due to collisions would exceed the loss of objects due to orbital decay.     

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.     

基于多项式近似的空间碎片群体轨道预报算法

快速准确地分析空间碎片群轨道演化行为对于其他在轨航天器碰撞规避至关重要。在各摄动力的作用下,空间碎片群演化运动呈现出复杂的非线性特征。空间碎片群体个体数量巨大,如果通过对空间碎片群中每个空间碎片进行轨道积分来分析群体预报的方法会导致计算量过大。针对该问题,提出一种基于多项式近似的轨道快速预报分析方法。该方法将空间碎片群分为少量的标称碎片和其他大量关联碎片。针对标称碎片的轨道预报采用数值积分求解保证预报精度;而针对其他大量的关联碎片轨道预报问题,采用多项式泰勒展开半解析方法求解,从而在保证预报精度的前提下有效减少空间碎片群轨道预报的计算量。为了验证方法的有效性,对不同空间碎片群进行了轨道预报仿真。仿真结果表明,当轨道预报精度设定在1m范围内时,多项式近似算法的计算量较蒙特卡洛方法计算效率提高了2.2~17.2倍,验证了所提出方法的有效性。     

超高速撞击航天器二次碎片云能量特性分析

为了评估空间碎片超高速撞击航天器的碎片云破坏能力,挖掘超高速撞击数值模 拟结果数据的应用价值,基于9.53 mm铝球以6.64 km/s速度对2.2 mm铝靶撞 击的Ls-Dyna/SPH(Smoothed Particle Hydrodynamic)数值模拟研究结果,对靶后碎片云的 粒子动能进行求和统计,建立了碎片云比动能概念和函数形式;碎片云比动能综合考虑了靶 后所有碎片云粒子的动能,反映了一定距离处垂直于撞击方向平面上单位面积上的碎片云粒 子所蕴含的撞击能量;应用碎片云比动能概念,揭示出随着演化距离的增加,碎片云能量的 衰减规律;通过不同速度条件下的SPH计算,得到了碎片云的比动能函数的曲线形式随撞击 速度的变化规律;最后对采用2种材料模型进行数值模拟所对应的结果误差进行碎片云比动 能函数的曲线比较,反映出数值模拟中不同材料模型引起的差异.      

Identification of a debris cloud from the nuclear powered SNAPSHOT satellite with haystack radar measurements

Data from the Massachusetts Institute of Technology Lincoln Laboratory Long Range Imaging Radar (known as the Haystack radar) have been used in the past to examine families of objects from individual satellite breakups or families of orbiting objects that can be isolated in altitude and inclination. This is possible because, for some time after a breakup, the debris cloud of particles can remain grouped together in similar orbit planes. This cloud will be visible to the radar, in fixed staring mode, for a short time twice each day, as the orbit plane moves through the field of view. There should be a unique three-dimensional pattern in observation time, range, and range rate which can identify the cloud. Eventually, through slightly differing precession rates of the right ascension of ascending node of the debris cloud, the observation time becomes distributed so that event identification becomes much more difficult.     

Characterization of the cataloged Fengyun-1C fragments and their long-term effect on the LEO environment

The intentional breakup of Fengyun-1C on 11 January 2007 created the most severe orbital debris cloud in history. The altitude where the event occurred was probably the worst location for a major breakup in the low Earth orbit (LEO) region, since it was already highly populated with operational satellites and debris generated from previous breakups. The addition of so many fragments not only poses a realistic threat to operational satellites in the region, but also increases the instability (i.e., collision cascade effect) of the debris population there.     

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.     

月地转移轨道快速设计与特性分析

对采用直接大气再入方式的月地转移轨道,考虑大气再入界面参数和地面落点位置约束,提出了一种基于双二体模型的快速设计方法。该方法分为内外两层迭代循环,内层循环使月心段轨道和地心段轨道在月球影响球边界处连续,并采用Lambert问题与Newton-Raphson法相集合的方法求解满足再入角约束的地心段轨道参数;外层循环通过调整地心段轨道倾角和轨道置入时间使月地转移轨道满足地面落点位置约束。分析表明,存在四种类型的月地转移轨道满足大气再入界面约束,分别为降 降型、降 升型、升 降型和升 升型。在此基础上,对四种类型月地转移轨道的近地点地心距、置入分布点、再入点分布等特性进行了分析。仿真结果验证了所提出方法的有效性。     

The impact of alternative mission models on the future orbital debris environment

In reviewing discussions of future directions for space activity, it becomes obvious that there are a large number of groups formulating a wide diversity of plans for the future use of space. These plan alternatives are being made to account for user needs, technology development constraints, economic constraints, and launch support, and each of the plans will have direct or indirect effects on the orbital debris environment in terms of mass to orbit, deposition of operational debris, and control of accidental breakups. Thus it is important to develop the ability to project future debris states for a range of possible space traffic scenarios. The impact that these possible traffic environments would have on space operations forms the basis for studies of alternative options for the usage of space. In this paper, the effects on the orbital debris environment of a base-line mission model and two alternatives are investigated, using a numerical debris environment simulation code under development at JSC.     

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.     

空间碎片云演变过程的阶段划分

根据碎片云从破碎点开始向空间扩散过程中碎片密度和形状的变化规律,以几何形状和起主要作用的因素为特征,定义了球形、椭球形、绳形、螺旋线形、全方位弥漫直至球壳形六个演变阶段.论述了在各个阶段的主要特征和对演变过程起主要作用的因素.总结了与演变过程相关的轨道运动理论和研究方法,分析了各个阶段演变的动力学原理.在球形阶段起主要作用的是分离速度;椭球形阶段可以利用线性化相对运动方程进行分析;绳形与螺旋线形在几何上有质变,但都有结点和结线,并可以利用速度增量理论分析和解释其存在的原因.轨道摄动力消除了结点和结线,导致碎片云的全方位弥漫,并最终使碎片云趋于球壳形.推导和罗列了各阶段转换标志点时刻的计算公式,利用计算机仿真的方法,给出了近地轨道各个阶段碎片云分布示意图,验证了演变过程阶段划分的合理性.     

A simple and valid analysis method for orbit anomaly detection

A simple analysis method for orbit anomaly detection, called semi-major axis change method (SACM) was presented by using a relationship between the change of orbit parameters and velocity increments. In this method, the mean value and standard deviation of the semi-major axis change in different time intervals were first calculated according to historical data. Then, these two parameters, the mean value and standard deviation of the semi-major axis change, are chosen as basis variables and combined as an anomalous criterion. For orbit objects with different characteristic, anomalous thresholds were given in different time intervals for identifying the anomalies of the orbital objects. Finally, this method is used for low earth orbit (LEO) satellites and American–Russian breakup debris. By adopting this method, the characteristics of the orbit change were given. The accuracy rate of anomaly analysis for LEO satellites and American–Russian breakup debris can reach to 100%, which demonstrates that the method was rapid and valid.     

Space debris mitigation strategies and practices in geosynchronous transfer orbits

Missions to geosynchronous orbits remain one of the most important elements of space launch traffic, accounting for 40% of all missions to Earth orbit and beyond during the four-year period 2000–2003. The vast majority of these missions leave one or more objects in geosynchronous transfer orbits (GTOs), contributing on a short-term or long-term basis to the space debris population. National and international space debris mitigation guidelines seek to curtail the accumulation of debris in orbits which penetrate the regions of low Earth orbit and of geosynchronous orbit. The orbital lifetime of objects in GTO can be greatly influenced by the initial values of perigee, inclination, and right ascension of the orbital plane, leading to orbital lifetimes of from less than one month to more than 100 years. An examination of the characteristic GTOs employed by launch vehicles from around the world has been conducted. The consequences of using perigees above 300 km and super-synchronous apogees, typically above 40,000 km, have been identified. In addition, the differences in orbital behavior of launch vehicle stages and mission-related debris in GTOs have been investigated. Greater coordination and cooperation between space launch service providers and spacecraft designers and owners could significantly improve overall compliance with guidelines to mitigate the accumulation of debris in Earth orbit.     

空间碎片自然衰减和太阳活动

分析研究了空间碎片数随太阳辐射流量Ｆ１０．７的变化;给出预报Ｆ１０．７长期变化的计算方法和预测空间碎片数的数学模型。结果显示：①强太阳活动造成空间碎片年增长率下降;②空间碎片数与太阳活动１１年变化密切相关,相关数为０．９;③空间碎片增长率约为发射率的两倍;④若发射率保持不变,则到２０２０年,大于１０ｃｍ的碎片数将达到１４５００;⑤若小碎片的增长为大碎片增长的两倍,则到２０２０年,大于１ｃｍ的碎片数可达１２５０００。