A Retrospective Geometric Analysis of the Long-Term Evolution of Orbits and the Ballistic Lifetime for Satellites of the Prognoz Series

We consider the applied aspects of the geometrical analysis of solutions in a restricted circular double-averaged three-body problem that are concerned with the design of high-apogee satellite orbits. Based on the analysis of the long-term evolution and the ballistic lifetime for orbit families of the Prognoz satellites launched into their orbits in the period 1972–1995, we suggest some practical suggestions for choosing long-lived high-apogee orbits with taking into account various requirements for the domain of evolution of the orbital elements.     

On the features of long-term evolution of a high-apogee orbit of the SPECTR-R spacecraft

The practical tasks related to qualitative investigation of long-term evolution of high-apogee orbits of artificial Earth satellites (AES), for which the main perturbing factors are gravitational perturbations from the Moon and the Sun, are considered. Attention is given to the problem of the ballistic lifetime of similar orbits, and the issues associated with possibilities of the correction of orbits for ensuring the required duration of their ballistic lifetime are considered. The orbit of the SPECTR-R spacecraft launched in July of 2011 is considered as an example.     

Choice of high-apogee AES orbits on the basis of the qualitative methods of the theory of perturbations and situational analysis. Part I. Situational studies based on orbital tori

The paper discusses the problems of the choice of high-apogee orbits of artificial Earth satellites (AES), proceeding from the tasks of space experiments aimed at studying near-earth space and taking into account the features of the orbital evolution and ballistic lifetime. The suggested methods of the choice of orbits consist of two components. The first is based on the use of mathematical models of studied regions of near-earth space and various techniques of situation analysis, among which the annual and daily orbital tori developed by the author about 35 years ago are key. The second component is based on qualitative methods of the theory of perturbations of high-apogee AES orbits developed by M.L. Lidov more than 50 years ago.     

同步双小行星系统共振轨道设计

研究了同步双小行星系统中共振轨道的设计方法及演化规律。首先，基于双椭球模型建立探测器运动方程，并给出共振轨道初值选取方法。然后，利用改进并行打靶法，提出一种双小行星系统平面共振轨道两步修正方法。同时结合稳定性理论及分岔理论，给出双小行星系统三维共振轨道生成和延拓方法；最后，以双小行星系统1999KW4为例，设计了共振比为1∶1,1∶2,1∶3,1∶4,2∶3的平面和空间共振轨道族，并分析了共振轨道的特性及轨道周期和轨道能量的变化规律。给出的双小行星系统中共振轨道的设计方法具有普适性，对未来双小行星系统探测任务中的轨道设计具有一定的参考意义与借鉴价值。     

Pervasive orbital eccentricities dictate the habitability of extrasolar earths

The long-term habitability of Earth-like planets requires low orbital eccentricities. A secular perturbation from a distant stellar companion is a very important mechanism in exciting planetary eccentricities, as many of the extrasolar planetary systems are associated with stellar companions. Although the orbital evolution of an Earth-like planet in a stellar binary system is well understood, the effect of a binary perturbation on a more realistic system containing additional gas-giant planets has been very little studied. Here, we provide analytic criteria confirmed by a large ensemble of numerical integrations that identify the initial orbital parameters leading to eccentric orbits. We show that an extrasolar earth is likely to experience a broad range of orbital evolution dictated by the location of a gas-giant planet, which necessitates more focused studies on the effect of eccentricity on the potential for life.     

太阳同步回归轨道的长期演变与控制

近地轨道的遥感卫星绝大部分都采用太阳同步回归轨道。这类轨道由于受到大气阻力的影响,半长轴将不断地衰变并导致地面轨迹的东漂,为保持回归特性需周期性地对半长轴进行调整。另一类长期变化是太阳引力引起的倾角变化,这是太阳同步轨道特有的。倾角长期的变化又进一步导致回归轨道的标称半长轴和降交点地方时的相应变化。文章给出了这些变化的解析模型以及轨道控制的策略。     

Continuum Model for Space Debris Evolution with Account of Collisions and Orbital Breakups

The paper discusses the mathematical modeling of long-term orbital debris evolution taking into account mutual collisions of space debris particles of different sizes. Investigations and long-term forecasts of orbital debris environment evolution in low Earth orbits are essential for future space mission hazard evaluation and for adopting rational space policies and mitigation measures. The paper introduces a new approach to space debris evolution mathematical modeling based on continuum mechanics incorporating partial differential equations. This is an alternative to the traditional approaches of celestial mechanics incorporating ordinary differential equations to model fragments evolution. The continuum approach to orbital debris evolution modeling has essential advantages for describing the evolution of a large number of particles, because it replaces the traditional tracking of space objects by modeling the evolution of their density of distribution.     

月球空间站停泊轨道选择分析

系统梳理了月球附近可用于部署空间站的停泊轨道类型,给出了不同类型轨道的定义及主要参数,通过仿真计算和对比国内外轨道研究结论,分析了不同类型轨道的能量需求、登月任务支持性、空间环境等特点,在此基础上,分析了美国"深空之门"空间站部署轨道选择的主要考虑和依据,并提出了各类轨道的应用建议,可为未来部署月球空间站的停泊轨道选择...     

Solution of the flyby problem for large space debris at sun-synchronous orbits

the paper considers the flyby problem related to large space debris (LSD) objects at low earth orbits. The data on the overall dimensions of known last and upper stages of launch vehicles makes it possible to single out five compact groups of such objects from the NORAD catalog in the 500–2000 km altitude interval. The orbits of objects of each group have approximately the same inclinations. The features of the mutual distribution of the orbital planes of LSD objects in the group are shown in a portrait of the evolution of deviations of the right ascension of ascending nodes (RAAN). In the case of the first three groups (inclinations of 71°, 74°, and 81°), the straight lines of relative RAAN deviations of object orbits barely intersect each other. The fourth (83°) and fifth (97°–100°) LSD groups include a considerable number of objects whose orbits are described by straight lines (diagonals), which intersect other lines many times. The use of diagonals makes it possible to significantly reduce the temporal and total characteristic velocity expenditures required for object flybys, but it complicates determination of the flyby sequence. Diagonal solutions can be obtained using elements of graph theory. A solution to the flyby problem is presented for the case of group 5, formed of LSD objects at sun-synchronous orbits.     

State transition matrix approximation with geometry preservation for general perturbed orbits

This paper presents a method for approximating the state transition matrix for orbits around a primary body and subject to arbitrary perturbations. The primary objective of this method is to provide an accurate state transition matrix for orbits with realistic perturbations, which has a sufficiently simple form for implementation onboard spacecraft. The averaging method is employed to isolate the high and low frequency spectra of the perturbation terms, and construct a functional form of the approximate state transition matrix composed only of elementary analytic functions. In addition to the methodology of the approximation, it is shown that the symplectic property, which is a fundamental mathematical structure of Hamiltonian systems, can be incorporated into this method. This not only reduces the number of parameters required for approximations but also makes it possible to preserve the physically true structure of the state transition matrix. The resulting state transition matrix approximation is valid for tens of orbital revolutions without having to update the parameters. Numerical simulations show that this method is valid for arbitrary eccentricity orbits with semimajor axis ranging from LEO up to around 10 Earth radii when applied to Earth orbits.     

Synergistic use of debris environment models for flexible, very-long term projections

Long-term debris environment projections are of great importance for assessing the necessity and effectiveness of debris mitigation measures. Two types of models have been developed to predict these environments. Environment evolution models like the EVOLVE code are using detailed mission model data to input spacecraft, upper stages, and operational debris into specific orbits at specific times; debris from fragmentations are placed in orbits defined by the state vector of the fragmenting object(s) and the breakup model. The second type, typified by the CHAIN program, uses a particle-in-box model that bins the environment in size and altitude rather than following the orbit evolution of individual debris fragments. A 3-Step approach using both the EVOLVE and CHAIN model in a synergistic way was used to increase the reliability of long term environment projections. EVOLVE historical projections 1957–1995 could be validated by comparison to measurements. The comparison of 100 year projection runs of EVOLVE and CHAIN for different traffic scenarios showed a good agreement. In this paper, for the first time, CHAIN projections up to 10,000 years, based on validated boundary conditions derived by EVOLVE are presented, indicating clearly the need of early implementation of effective mitigation measures to prevent exponential population growth by collisional cascading effects.     

Control of satellite imaging formations in multi-body regimes

Libration point orbits may be ideal locations for satellite imaging formations. Therefore, control of these arrays in multi-body regimes is critical. A continuous feedback control algorithm is developed that maintains a formation of satellites in motion that is bounded relative to a halo orbit. This algorithm is derived based on the dynamic characteristics of the phase space near periodic orbits in the circular restricted three-body problem (CR3BP). By adjusting parameters of the control algorithm appropriately, satellites in the formation follow trajectories that are particularly advantageous to imaging arrays. Image reconstruction and coverage of the (u, v) plane are simulated for interferometric satellite configurations, demonstrating potential applications of the algorithm and the resulting motion.     

Optimization of space orbits design for Earth orbiting missions

In Earth orbiting space missions, the orbit selection dictates the mission parameters like the ground resolution, the area coverage, and the frequency of coverage parameters. To achieve desired mission parameters, usually Earth regions of interest are identified and the spacecraft is maneuvered continuously to visit only these regions. This method is expensive, it requires a propulsion system onboard the spacecraft, working throughout the mission lifetime. It also requires a longer time to cover all the regions of interest, due to the very weak thrust forces compared to that of the Earth's gravitational field. This paper presents a methodology to design natural orbits, in which the regions of interest are visited without the use of propulsion systems, depending only on the gravitational forces. The problem is formulated as an optimization problem. A genetic algorithm along with a second order gradient method is implemented for optimization. The design process takes into consideration the gravitational second zonal harmonic, and hence allows for the design of repeated Sun-synchronous orbits. The field of view of the payload is also taken into consideration in the optimization process. Numerical results are presented that demonstrates the efficiency of the proposed method.     

太阳帆绕地球周期轨道研究

  地球同步和太阳同步卫星在各个领域有着广泛的应用。静止轨道是一种特殊的地球同步轨道,轨道资源有限。利用化学推进或电推进可以实现轨道高度不同的同步轨道,如悬挂轨道,但需要消耗较多的燃料,工程上无法承受。本文考虑利用太阳帆实现地球同步和太阳同步轨道。太阳光压力在轨道平面内沿拱线方向,选择光压力与平面的夹角使得轨道平面的旋转速率与太阳光同步。研究表明,设计合适的半长轴和偏心率可以使得轨道旋转速率与地球自转速率一致。假设太阳光与赤道平面平行,可以得到准静止轨道,太阳帆将在传统静止轨道的附近运动,星下点的经度将在一个固定值附近振动。实际上太阳光是与黄道面平行,黄道面与赤道面之间存在夹角。考虑黄赤交角的情况下,太阳帆将在一定纬度和经度范围内运动。适合于对某个区域进行长期观测任务。     

基于傅立叶平均法下的连续小推力动力学分析

通过将小推力展开为偏近点角的傅立叶级数,并对高斯摄动方程在一个轨道周期上的平均,将原方程的推力转化为仅由14个傅立叶系数表示的控制变量。仿真计算表明,平均化后的高斯方程使计算量与牛顿积分相比显著减少,且对小推力而言有足够的精度。对利用平均化后的高斯方程计算轨道根数时产生误差的原因进行了研究,并进一步分析小推力的范围和小推力近似表达式对上述误差的影响,为今后小推力下非开普勒轨道动力学分析提供了理论依据和参数。     

State propagation in an uncertain asteroid gravity field

Various spacecraft have been and will be sent to asteroids to characterize them. Generally, an asteroid's gravity field is very irregular and not accurately known when compared to the gravity field of a major planet, Earth in particular. It has been well studied that the irregularity significantly affects the trajectory of an orbiting spacecraft, and causes it to impact or to escape from the asteroid. Complementary to that, this paper focuses on the influence of the limited knowledge of this gravity field on the evolution of the spacecraft's orbit. It develops a general method by which this influence can be quantified. This method comprises specific Monte Carlo simulations with a discrete set of low-altitude orbits, taking into account the uncertainties in the gravity-field parameters. For illustration purposes, it is applied to two different asteroids. Already after three revolutions, the gravity-field uncertainties propagate to significant position uncertainties; this specifically holds for prograde orbits, and around the smaller asteroid. Applying this robust and accurate method helps mission designers and planners to assess the risk posed by gravity uncertainties, and take appropriate measures such as choosing the most favorable orbital geometries and/or lowering the orbit more slowly.     

Particle swarm optimization applied to impulsive orbital transfers

The particle swarm optimization (PSO) technique is a population-based stochastic method developed in recent years and successfully applied in several fields of research. It mimics the unpredictable motion of bird flocks while searching for food, with the intent of determining the optimal values of the unknown parameters of the problem under consideration. At the end of the process, the best particle (i.e. the best solution with reference to the objective function) is expected to contain the globally optimal values of the unknown parameters. The central idea underlying the method is contained in the formula for velocity updating. This formula includes three terms with stochastic weights. This research applies the particle swarm optimization algorithm to the problem of optimizing impulsive orbital transfers. More specifically, the following problems are considered and solved with the PSO algorithm: (i) determination of the globally optimal two- and three-impulse transfer trajectories between two coplanar circular orbits; (ii) determination of the optimal transfer between two coplanar, elliptic orbits with arbitrary orientation; (iii) determination of the optimal two-impulse transfer between two circular, non-coplanar orbits; (iv) determination of the globally optimal two-impulse transfer between two non-coplanar elliptic orbits. Despite its intuitiveness and simplicity, the particle swarm optimization method proves to be capable of effectively solving the orbital transfer problems of interest with great numerical accuracy.     

Elegant propositions of the boundary-value problem

The fundamental ellipse of the two-body, two-point orbital boundary-value problem of orbital mechanics is introduced and its key role in the basic geometry of this famous problem is developed. Analytic formulas for the orbital parameters of conjugate orbits are derived together with further interesting properties of the mean point-locus.     

An analysis of peculiarities of the satellite autonomous navigation in sighting unknown landmarks on the planet surface

An analysis of characteristics of autonomous navigation by using optical sighting of unknown landmarks was carried out for artificial satellites in circular orbits around the Earth with altitudes from 300 to 900 km. It is shown that navigational errors vary approximately in inverse proportion to duration of measurements. Also they decrease with an increasing number of sighted landmarks. The errors poorly depend on the orbit altitude if the duration of measurement for each landmark is chosen in proportion to the orbit altitude, or a maximal zenith angle of sighting the landmark remains constant.Some characteristics of the navigational algorithm are given. Admissible deviations of initial orbital elements from their exact values are investigated. For a properly chosen measurement procedure the algorithm proves to be stable even for rather crudely given initial parameters.     

火星卫星的冻结轨道研究

冻结轨道是一种稳定的轨道,地球、火星、月球的卫星因引力场的南北不对称,都存在冻结轨道.由于主星体引力场的不同,它们卫星的冻结轨道也有不同的特性.地球卫星的冻结执道的偏心率非常小,对卫星遥感非常有利,国内外已有相当多的近地遥感卫星采用这种轨道.月球卫星的冻结轨道偏心率随轨道倾角的不同有很大的变化,对月球卫星冻结轨道的研究...