Artificial satellites orbits in 2:1 resonance: GPS constellation

In this work, the resonance problem in the artificial satellites motion is studied. The development of the geopotential includes the zonal harmonics J₂₀ and J₄₀ and the tesseral harmonics J₂₂ and J₄₂. Through an averaging procedure and successive Mathieu transformations, the order of dynamical system is reduced and the final system is solved by numerical integration. In the simplified dynamical model, three critical angles are studied. The half-width of the separatrix is calculated through a linearized model which describes the behavior of the dynamical system in a neighborhood of each critical angle. Through the resonance overlap criterion the possible regular and irregular motions are investigated by the time behavior of the semi-major axis, argument of perigee and eccentricity. The largest Lyapunov exponent is used as tool to verify the chaotic motion.     

Orbit prediction and shuttle rendezvous strategy for eureca

The European Retrievable Carrier (EURECA) is a platform to be launched, deployed and retrieved in low Earth orbit by the Space Shuttle.A newly developed analytical orbit prediction method is described which meets the severe requirements for EURECA's orbit propagation. It is based on an averaging procedure including the Earth's zonal harmonics J₂, J₃ and J₄ and a refined treatment of the air drag perturbation where EURECA's large solar panels are taken into account. Some orbit prediction results are included.In order to offer more flexibility for the Shuttle retrieval of EURECA, it is proposed to execute a part of the rendezvous manoeuvres by EURECA. A corresponding strategy is described.     

Lorentz-force-perturbed orbits with application to J2-invariant formation

The Lorentz force acting on an electrostatically charged spacecraft in the Earth's magnetic field provides a new propellantless means for controlling a spacecraft's orbit. Assuming that the Lorentz force is much smaller than the gravitational force, the perturbation of a charged spacecraft's orbit by the Lorentz force in the Earth's magnetic field, which is simplified as a titled rotating dipole, is studied in this article. Our research starts with the derivation of the equations of motion in geocentric equatorial inertial Cartesian coordinates using Lagrange mechanics, and then derives the Gauss variational equations involving Lorentz-force perturbation using a set of nodal inertial coordinates as an intermediate step. Subsequently, the approximate averaged changes in classical orbital elements, including single-orbit-averaged and one-day-averaged changes, are obtained by employing orbital averaging. We have found that the approximate analytic one-day-averaged changes in semi-major axis, eccentricity, and inclination are nearly zero, and those in the other three angular orbital elements are affected by J₂ and Lorentz-force perturbations. This characteristic is applied to model bounded relative orbital motion in the presence of the Lorentz force, which is termed Lorentz-augmented J₂-invariant formation. The necessary condition for J₂-invariant formation is derived when the chief spacecraft's reference orbit is either circular or elliptical. It is shown that J₂-invariant formation is easier to implement if the deputy spacecraft is capable of establishing electric charge. All conclusions drawn from the approximate analytic solutions are verified by numerical simulation.     

Nonlinear relative position control of precise formation flying using polynomial eigenstructure assignment

A nonlinear relative position control algorithm is designed for spacecraft precise formation flying. Taking into account the effect of J₂ gravitational perturbations and atmospheric drag, the relative motion dynamic equation of the formation flying is developed in a quasi-linear parameter-varying (QLPV) form without approximation. Base on this QLPV model, polynomial eigenstructure assignment (PEA) is applied to design the controller. The resulting PEA controller is a function of system state and parameters, and produces a closed-loop system with invariant performance over a wide range of conditions. Numerical simulation results show that the performance can fulfill precise formation flying requirements.     

Lagrangian coherent structures in various map representations for application to multi-body gravitational regimes

The Finite-Time Lyapunov Exponent (FTLE) has been demonstrated as an effective metric for revealing distinct, bounded regions within a flow. The dynamical differential equations derived in multi-body gravitational environments model a flow that governs the motion of a spacecraft. Specific features emerge in an FTLE map, denoted Lagrangian Coherent Structures (LCS), that define the extent of regions that bound qualitatively different types of behavior. Consequently, LCS supply effective barriers to transport in a generic system, similar to the notion of invariant manifolds in autonomous systems. Unlike traditional invariant manifolds associated with solutions in an autonomous system, LCS evolve with the flow in time-dependent systems while continuing to bound distinct regions of behavior. Moreover, in general, FTLE values supply information describing the relative sensitivity in the neighborhood of a trajectory. Here, different models and variable representations are used to generate maps of FTLE, and the resulting structures are applied to design and analysis within an astrodynamical context. Application of FTLE and LCS to transfers from LEO to the L₁ region in the Earth–Moon system are presented and discussed. In an additional example, an FTLE analysis is offered of a few stationkeeping maneuvers from the Earth–Moon mission ARTEMIS (Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun).     

On Saturn’s rotation relative to a center of mass under the action of the gravitational moments of the Sun and Jupiter

Saturn’s rotation relative to a center of mass is considered within an elliptic restricted three-body problem. It is assumed that Saturn is a solid under the action of gravity of the Sun and Jupiter. The motions of Saturn and Jupiter are considered elliptic with small eccentricities e_S and e_J, respectively; the mean motion of Jupiter n_J is also small. We obtain the averaged Hamiltonian function for a small parameter of ε = n_J and integrals of evolution equations. The main effects of the influence of Jupiter on Saturn’s rotation are described: (α) the evolution of the constant parameters of regular precession for the angular momentum vector I₂; (β) the occurrence of new libration zones of oscillations I₂ near the plane of the celestial equator parallel to the plane of the Jupiter’s orbit; (γ) the occurrence of additional unstable equilibria of vector I₂ at the points of the north and south poles of the celestial sphere and, as a result, the existence of homoclinic trajectories; and (δ) the existence of periodic trajectories with arbitrarily large periods near the homoclinic trajectory. It is shown that the effects of (β), (γ), and (δ) are caused by the eccentricity e of the Jupiter’s orbit and are practically independent of Jupiter’s mass (within satellite approximation).     

Acceleration of solar cosmic rays and the fine spectral structure of type II radio bursts

On the basis of data, obtained by means of the ground-based solar service RSTN (Radio Solar Telescope Network) and the geostationary satellite system GOES, the relationship between the solar cosmic rays (SCR) intensity I _p with the proton energy E _p > 1 MeV and parameters of meter-decameter type II radio bursts in the frequency range of 25–180 MHz is studied. The process of proton acceleration by shock waves was characterized by the frequency drift velocity of radio bursts V _mII and the relative difference between radio emission frequencies at the first two harmonics b. It is shown that the coefficient of correlation between I _p and b increases with E _p growing from 0.40 to 0.70, while a similar coefficient between I _p and V _mII does not exceed 0.30. Indications in favor of the two-stage SCR acceleration model are obtained.     

Equilibrium attitude and stability of a spacecraft on a stationary orbit around an asteroid

The stationary orbits around an asteroid, if exist, can be used for communication and navigation purposes just as around the Earth. The equilibrium attitude and stability of a rigid spacecraft on a stationary orbit around a uniformly-rotating asteroid are studied. The linearized equations of attitude motion are obtained under the small motion assumption. Then, the equilibrium attitude is determined in both cases of a general and a symmetrical spacecraft. Due to the higher-order inertia integrals of the spacecraft, the equilibrium attitude is slightly away from zero Euler angles. Then necessary conditions of stability of this conservative system are analyzed based on the linearized equations of motion. The effects of different parameters, including the harmonic coefficients C₂₀ and C₂₂ of the asteroid and higher-order inertia integrals of the spacecraft, on the stability are assessed and compared. Due to the significantly non-spherical shape and rapid rotation of the asteroid, the effects of the harmonic coefficients C₂₀ and C₂₂ are very significant, while effects of the third- and fourth-order inertia integrals of the spacecraft can be neglected. Considering a spacecraft on a stationary orbit around an example asteroid, we show that the classical stability domain predicted by the Beletskii–DeBra–Delp method on a circular orbit in a central gravity field is modified due to the non-spherical mass distribution of the asteroid. Our results are confirmed by a numerical simulation.     

Luni-solar perturbations in the extended phase space representation of the Vinti problem

Perturbation theory is applied to the Vinti problem—motion about an oblate spheroid—to include the gravitational effects of the sun and moon. The problem is formulated using the extended phase space method which introduces a new independent variable similar to the true anomaly. The disturbing Hamiltonian H₁ for third bodies is of order J₂² (second order) and the final goal is a theory including second order short and long period terms and third order secular terms. The current paper however carries the development only to the second order in the secular terms and the first order in the periodic terms. Problems of including the higher orders are discussed. Therefore, in the development of H₁ all terms of order 10^?9 or larger are retained. The lunar emphemeris retains terms to e′² in the lunar eccentricity. The perturbation analysis is carried out by means of Lie series and is developed through the first order only which is consistent with the final accuracy desired. The generating function W₁ is obtained and separated into the long period, short period and secular terms. From W₁ the coordinates are defined from the Lie series by means of a transformation equation. These coordinates are non-singular for small eccentricity and inclination. Because of the complexity of the equations all algebraic computations were accomplished by means of a computerized Poisson series manipulator developed at the Naval Research Laboratory.     

航天器共振轨道研究

建立航天器非开普勒运动理论是航天技术发展的必然。提出了一类非开普勒轨道——共振轨 道。共振是自然界的一种普遍现象,当发生共振时,很小的输入可以使系统的状态产生较大 变化。研究表明航天器在推力作用下的非开普勒运动在参数平面内可以视为一种受迫振 动,也会发生共振现象。因此,可以利用共振原理来研究航天器的运动,称这样一类非开普 勒轨道为共振轨道。首先通过合理地选择轨道描述参数、时间尺度和推力描述方式建立 航天器共振轨道的动力学模型。然后讨论航天器在推力作用下轨道运动的振动规律,并给出 共振轨道的概念及轨道方程。最后提出基于共振轨道的机动轨道设计方法。
     

Mapping orbits with low station keeping costs for constellations of satellites based on the integral over the time of the perturbing forces

The present paper is concerned with the search for orbits that have potential to require low fuel consumption for station-keeping maneuvers for constellations of satellites. The method used to study this problem is based on the integral over the time of the undesired perturbing forces. This integral measures the change of velocity caused by the perturbation forces acting on the satellite, so mapping orbits that are less perturbed, which generates good candidates for orbits that requires low fuel consumption for station-keeping maneuvers. The integral over the time depends only on the orbit of the spacecraft and the dynamical system considered. The type of engine and the control technique applied to the spacecraft are not considered to search for those orbits. It can be a good strategy to be applied for a first mapping of orbits. For this search, it is analyzed the integral of orbits with different values of the Keplerian elements in order to find the best ones with respect to this criterion. The perturbations considered are the ones caused by the third body, which includes the Sun and the Moon, and the J₂ term of the geopotential. The results presented here show numerical simulations to obtain the integral of those perturbing forces for different orbits. The GPS and the Molniya constellations are used as examples for those calculations.     

含板类柔性附件的充液航天器动力学研究

对含有板类柔性附件和曲壁轴对称充液储腔的复杂航天器系统进行动力学建模和耦合机理研究。首先,采用Kirchhoff-Love薄板理论对航天器的板类柔性附件进行研究,通过D’Alembert原理得到柔性附件的振动方程,运用模态假设法将混合方程转换为常微分方程。其次,通过推导充液航天器储腔内任意点的运动,得到储腔液体的牵连速度势函数,采用Gauss超几何级数得到液体相对速度势函数的解析形式,通过Hamilton变分原理推导液体晃动的运动方程,以及液体速度势函数模态系数的控制方程。最后采用准坐标Lagrange方程得到耦合航天器系统的状态方程,通过数值仿真校验系统动力学模型的有效性。研究结果表明,刚性平台、液体、柔性附件的相互耦合效应使得航天器系统存在复杂动力学行为,在复杂航天器系统动力学建模过程中需要充分考虑液体晃动和柔性附件振动的影响,柔性附件的安装位置对于耦合航天器系统的动力学行为也有着重要影响。     

Variations of the magnetopause position versus the level of geomagnetic activity (according to data of the INTERBALL-1 Satellite for 1995–1997)

The variations in the deviation of the observed position of the magnetosphere boundary from its mean position predicted by the Shue at al., 1997 (Sh97) model [7] are studied as a function of the substorm activity level (the AE-index value) and magnetic storm intensity (the value of the corrected D _st ^* index). The results obtained make it possible to state that the amplitude of motion of the magnetospheric boundary on the dayside and in the low-latitude tail is small. It is likely that the position of the boundary is either independent of the AE and D _st ^* indices or this dependence is weak. At the same time, the boundary of the high-latitude tail shifts inward on the average by 1.5R _E with an increase of the AE-index in the case of absence of magnetic storms (contraction of the magnetospheric tail). On the contrary, in the presence of magnetic storms, this boundary shifts outward by up to 3R _E with an increase of the AE-index (inflation of the magnetospheric tail). It is also shown that the boundary of the high-latitude tail moves outward with an increase of the D _st ^* index, both at low substorm activity and in periods of high substorm activity. The amplitude of the outward motion of the high-latitude tail of the magnetosphere is by a factor of two higher for moderate magnetic storms with strong substorms than for moderate magnetic storms with weak substorms.     

Halo orbits in the vicinity of the L 2 libration point in the Sun-Earth system

Methods are proposed for constructing the orbits of spacecraft remaining for long periods of time in the vicinity of the L ₂ libration point in the Sun-Earth system (so-called halo orbits), and the trajectories of uncontrolled flights from low near-Earth orbits to halo orbits. Halo orbits and flight trajectories are constructed in two stages: A suitable solution to a circular restricted three-body problem is first constructed and then transformed into the solution for a restricted four-body problem in view of the real motions of the Sun, Earth, and Moon. For a halo orbit, its prototype in the first stage is a combination of a periodic Lyapunov solution in the vicinity of the L ₂ point and lying in the plane of large-body motion, with the solution for the linear second-order system describing small deviations of the spacecraft from this plane along the periodic solution. The desired orbit is found as the solution to the three-body problem best approximating the prototype in the mean square. The constructed orbit serves as a similar prototype in the second stage. In both stages, the approximating solution is constructed by continuation along a parameter that is the length of the approximation interval. Flight trajectories are constructed in a similar manner. The prototype orbit in the first stage is a combination of a solution lying in the plane of large-body motion and a solution for a linear second-order system describing small deviations of the spacecraft from this plane. The planar solution begins near the Earth and over time tends toward the Lyapunov solution existing in the vicinity of the L ₂ point. The initial conditions of both prototypes and the approximating solutions correspond to the spacecraft’s departure from a low near-Earth orbit at a given distance, perigee, and inclination.     

Bifurcation Regime of a Synchronous Resonance in the Translational–Rotational Motion of Nonspherical Natural Satellites of Planets

Period-doubling bifurcations of the synchronous spin-orbit resonance in the motion of a nonspherical natural planetary satellite along the elliptic orbit are studied. The satellite spin axis is assumed to coincide with the axis of its largest principal moment of inertia and is perpendicular to the orbital plane. The period-doubling bifurcations take place when the value of satellite's dynamical asymmetry parameter falls in the parametric resonance domain. Theoretical dependences of the amplitude of the bifurcation oscillations of a satellite at the pericenter of its orbit upon the eccentricity and dynamical asymmetry parameter are investigated. Three different methods of calculating the amplitude of bifurcation oscillations are presented and compared. These theoretical estimates can be used to predict the opportunity to observe the bifurcation regime. The possibility of the occurrence of the bifurcation regime in the motion of natural planetary satellites is studied. It is concluded that the bifurcation regime is possible in the motion of Deimos, Epimetheus, Helen, Pandora, and Phobos. Phobos is the most probable candidate for finding the bifurcation regime of a synchronous rotation. The identification of such a regime would allow one to impose stringent constraints on the values of the inertial parameters of the satellite observed.     

Applications of multi-body dynamical environments: The ARTEMIS transfer trajectory design

The application of forces in multi-body dynamical environments to permit the transfer of spacecraft from Earth orbit to Sun–Earth weak stability regions and then return to the Earth–Moon libration (L₁ and L₂) orbits has been successfully accomplished for the first time. This demonstrated that transfer is a positive step in the realization of a design process that can be used to transfer spacecraft with minimal Delta-V expenditures. Initialized using gravity assists to overcome fuel constraints; the ARTEMIS trajectory design has successfully placed two spacecrafts into Earth–Moon libration orbits by means of these applications.     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.     

A small parameter method in problems of maneuvering space vehicles with aerodynamic efficiency

The motion of a satellite with aerodynamic efficiency along a low near-circular orbit is considered in the paper. The controls of bank angle γ and lift coefficient C_y are used as control functions. The introduction of a small parameter ($\text{? = (}\text{?}{}_0\text{· S · g}{}_0\text{2G}\text{)}$) makes it possible to integrate an adjoint system of equations and to obtain an approximate solution to the complete problem in the class of piecewise-constant control functions. Maximum values for the coordinates of heading angle η and lateral derivation from the plane of a reference orbit ?, which are connected with orbit plane angle by the relation cos i = cos ? · cos h, are used as criteria of maneuvering capability for a satellite with aerodynamic efficiency. Optimal programs for bank angle and incidence variation are derived and the influence of lift-to-drag ratio on the vehicle maneuvering capabilities has been estimated.It is shown that the process of the optimal motion is a special kind of gravitational skipping similar to the Keplerian motion but with continuous descent.     

绳系空间太阳能电站动力学响应分析

利用绝对节点坐标方法研究绳系空间太阳能电站在轨飞行的太阳能电池板动力响应。通过勒让德变换引入广义动量,在约束哈密尔顿体系下建立轨道、姿态和弹性振动耦合的动力学方程。基于祖冲之类方法的思想,结合辛龙格-库塔方法对微分-代数方程进行数值求解。数值算例说明本文建模方法和数值算法都是有效的,能很好地保持系统约束和能量。最后分析了绳长、平台系统的质量、轨道高度对于梁中点挠度和轴向平均应变的影响。