On the true circular orbit of a satellite

The osculating orbit of a planetary satellite moving in the equatorial plane of the central body under the influence of a rotational symmetric perturbation force is elliptical in first order approximation even if the true orbit is always circular. The satellite motion is influenced by a resonance effect due to this perturbing force. An inclined true satellite orbit cannot be circular.     

Some orbital characteristics of lunar artificial satellites

In this paper we present an analytical theory with numerical simulations to study the orbital motion of lunar artificial satellites. We consider the problem of an artificial satellite perturbed by the non-uniform distribution of mass of the Moon and by a third-body in elliptical orbit (Earth is considered). Legendre polynomials are expanded in powers of the eccentricity up to the degree four and are used for the disturbing potential due to the third-body. We show a new approximated equation to compute the critical semi-major axis for the orbit of the satellite. Lie-Hori perturbation method up to the second-order is applied to eliminate the terms of short-period of the disturbing potential. Coupling terms are analyzed. Emphasis is given to the case of frozen orbits and critical inclination. Numerical simulations for hypothetical lunar artificial satellites are performed, considering that the perturbations are acting together or one at a time.     

月球探测器过渡轨道的短弧定轨方法

论述的短弧定轨,是指在无先验信息情况下又避开多变元迭代的初轨计算方法,它需要相应的动力学问题有一能反映短弧内达到一定精度的近似分析解．探测器进入月球引力作用范围后接近月球时可以处理成相对月球的受摄二体问题,而在地球附近,则可处理成相对地球的受摄二体问题,但在整个过渡段的力模型只能处理成一个受摄的限制性三体问题．而限制性三体问题无分析解,即使在月球引力作用范围外,对于大推力脉冲式的过渡方式,相对地球的变化椭圆轨道的偏心率很大(超过Laplace极限),在考虑月球引力摄动时亦无法构造摄动分析解．就此问题,考虑在地球非球形引力(只包含J2项)和月球引力共同作用下,构造了探测器飞抵月球过渡轨道段的时间幂级数解,在此基础上给出一种受摄二体问题意义下的初轨计算方法,经数值验证,定轨方法有效,可供地面测控系统参考．     

基于STK的寿命末期GEO卫星轨道演化分析

地球同步轨道上的GEO卫星接近寿命末期时,星上的有效载荷和电子元器件等仍能正常应用,此时如果只做东西方向的位置保持,可以大幅度延长卫星的在轨工作寿命。简要介绍了引起GEO卫星轨道变化的主要摄动力,利用STK软件分析寿命末期GEO卫星的运动规律,并给出这类GEO卫星轨道演化的仿真分析流程,用以指导寿命末期GEO卫星的开发应用工作。     

Tilting Styx and Nix but not Uranus with a Spin-Precession-Mean-motion resonance

A Hamiltonian model is constructed for the spin axis of a planet perturbed by a nearby planet with both planets in orbit about a star. We expand the planet–planet gravitational potential perturbation to first order in orbital inclinations and eccentricities, finding terms describing spin resonances involving the spin precession rate and the two planetary mean motions. Convergent planetary migration allows the spinning planet to be captured into spin resonance. With initial obliquity near zero, the spin resonance can lift the planet’s obliquity to near 90\(^\circ \) or 180\(^\circ \) depending upon whether the spin resonance is first or zeroth order in inclination. Past capture of Uranus into such a spin resonance could give an alternative non-collisional scenario accounting for Uranus’s high obliquity. However, we find that the time spent in spin resonance must be so long that this scenario cannot be responsible for Uranus’s high obliquity. Our model can be used to study spin resonance in satellite systems. Our Hamiltonian model explains how Styx and Nix can be tilted to high obliquity via outward migration of Charon, a phenomenon previously seen in numerical simulations.     

The theory of asynchronous relative motion I: time transformations and nonlinear corrections

Using alternative independent variables in lieu of time has important advantages when propagating the partial derivatives of the trajectory. This paper focuses on spacecraft relative motion, but the concepts presented here can be extended to any problem involving the variational equations of orbital motion. A usual approach for modeling the relative dynamics is to evaluate how the reference orbit changes when modifying the initial conditions slightly. But when the time is a mere dependent variable, changes in the initial conditions will result in changes in time as well: a time delay between the reference and the neighbor solution will appear. The theory of asynchronous relative motion shows how the time delay can be corrected to recover the physical sense of the solution and, more importantly, how this correction can be used to improve significantly the accuracy of the linear solutions to relative motion found in the literature. As an example, an improved version of the Clohessy-Wiltshire (CW) solution is presented explicitly. The correcting terms are extremely compact, and the solution proves more accurate than the second and even third order CW equations for long propagations. The application to the elliptic case is also discussed. The theory is not restricted to Keplerian orbits, as it holds under any perturbation. To prove this statement, two examples of realistic trajectories are presented: a pair of spacecraft orbiting the Earth and perturbed by a realistic force model; and two probes describing a quasi-periodic orbit in the Jupiter-Europa system subject to third-body perturbations. The numerical examples show that the new theory yields reductions in the propagation error of several orders of magnitude, both in position and velocity, when compared to the linear approach.     

Analytical direct planetary perturbations on the orbital motion of satellites

An analytical expansion of the disturbing function arising from direct planetary perturbations on the motion of satellites is derived. As a Fourier series, it allows the investigation of the secular effects of these direct perturbations, as well as of every argument present in the perturbation. In particular, we construct an analytical model describing the evection resonance between the longitude of pericenter of the satellite orbit and the longitude of a planet, and study briefly its dynamic. The expansion developed in this paper is valid in the case of planar and circular planetary orbits, but not limited in eccentricity or inclination of the satellite orbit.     

Determination of an intermediate perturbed orbit from three position vectors

We suggest a new approach to solving the problem of finding the orbit of a celestial body from its three spatial position vectors and the corresponding times. It allows most of the perturbations in the motion of a celestial body to be taken into account. The approach is based on the theory of intermediate orbits that we developed previously. We construct the orbit the motion along which is a combination of two motions: the motion of a fictitious attracting center whose mass varies according to Mestschersky’s first law and the motion relative to the fictitious center. The first motion is generally parabolic, while the second motion is described by the equations of the Gylden-Mestschersky problem. The constructed orbit has such parameters that their limiting values at any reference epoch define a superosculating intermediate orbit with a fourth-order tangency. We have performed a numerical analysis to estimate the accuracy of approximating the perturbed motion of two minor planets, 145 Adeona and 4179 Toutatis, by the orbits computed using two-position procedures (the classical Gauss method and the method that we suggested previously), a three-position procedure based on the Herrick-Gibbs equation, and the new method. Comparison of the results obtained suggests that the latter method has an advantage.     

A set modified equinoctial orbit elements

Modified equinoctial elements are introduced which are suitable for perturbation analysis of all kinds of orbit. Equations of motion in Lagrangian and Gaussian forms are derived. Identities connecting the partial derivatives of the disturbing function with respect to equinoctial elements are established. Numerical comparisons of the evolution of a perturbed, highly eccentric, elliptic orbit analysed in equinoctial elements and by Cowell's method show satisfactory agreement.     

北斗三号卫星精密定轨中的光压模型研究

对于在轨运行的BDS (BeiDou Navigation Satellite System)卫星, 太阳光压是作用在卫星上主要的非引力摄动. 受多种因素的影响, 太阳光压摄动力难以精确建模, 是BDS卫星精密定轨和轨道预报过程中重要的误差来源. 由于ECOMC (Empirical CODE Orbit Model 1 and 2 Combined)模型兼顾了ECOM1 (Empirical CODE Orbit Model 1)和ECOM2 (Empirical CODE Orbit Model 2)模型的特点, 在模型中引入了较多的待估参数, 使得参数之间存在强相关性. 针对ECOMC模型的这一缺陷, 文中收集了2019年1月至2022年4月武汉大学分析中心提供的BDS-3卫星精密星历, 采用动力学轨道拟合方法得到了ECOMC模型的13个光压参数. 通过对该模型的光压参数进行时间序列分析, 分别给出了BDS-3 IGSO (Inclined Geosynchronous Orbit)和MEO (Medium Earth Orbit)卫星光压模型的参数选择策略. 并利用轨道拟合和轨道预报试验, 验证了光压模型参数选择策略的合理性. 结果表明, 采用改进型ECOMC模型进行BDS-3 IGSO和MEO卫星轨道拟合的效果最佳, 同时, 也能够提升BDS-3 IGSO和MEO卫星中长期轨道预报的精度.     

Satellite orbits perturbed by direct solar radiation pressure: General expansion of the disturbing function

An expression is derived for the solar radiation pressure disturbing function on an Earth satellite orbit which takes into account the variation of the solar radiation flux with distance from the Sun's centre and the absorption of radiation by the satellite. This expression is then expanded in terms of the Keplerian elements of the satellite and solar orbits using Kaula's method. The Kaula inclination functions are replaced by an equivalent set of modified Allan inclination functions.The resulting expression reduces to the form commonly used in solar radiation pressure perturbation studies (e.g. Aksnes, 1976), when certain terms are neglected. If, as happens quite often in practice, a satellite's orbit is in near-resonsnce with certain of these neglected terms, these near-resonant terms can cause changes in the satellite's orbital elements comparable to those produced by the largest term in Aksnes's expression. A new expression for the solar radiation pressure disturbing function expansion is suggested for use in future studies of satellite orbits perturbed by solar radiation pressure.     

Integrable approximation of J 2-perturbed relative orbits

Most existing satellite relative motion theories utilize mean elements, and therefore cannot be used for calculating long-term bounded perturbed relative orbits. The goal of the current paper is to find an integrable approximation for the relative motion problem under the J ₂ perturbation, which is adequate for long-term prediction of bounded relative orbits with arbitrary inclinations. To that end, a radial intermediary Hamiltonian is utilized. The intermediary Hamiltonian retains the original structure of the full J ₂ Hamiltonian, excluding the latitude dependence. This formalism provides integrability via separation, a fact that is utilized for finding periodic relative orbits in a local-vertical local-horizontal frame and determine an initialization scheme that yields long-term boundedness of the relative distance. Numerical experiments show that the intermediary-based computation of orbits provides long-term bounded orbits in the full J ₂ problem for various inclinations. In addition, a test case is shown in which the radial intermediary-based initial conditions of the chief and deputy satellites yield bounded relative distance in a high-precision orbit propagator.     

Approximate analysis for relative motion of satellite formation flying in elliptical orbits

This paper studies the relative motion of satellite formation flying in arbitrary elliptical orbits with no perturbation. The trajectories of the leader and follower satellites are projected onto the celestial sphere. These two projections and celestial equator intersect each other to form a spherical triangle, in which the vertex angles and arc-distances are used to describe the relative motion equations. This method is entitled the reference orbital element approach. Here the dimensionless distance is defined as the ratio of the maximal distance between the leader and follower satellites to the semi-major axis of the leader satellite. In close formations, this dimensionless distance, as well as some vertex angles and arc-distances of this spherical triangle, and the orbital element differences are small quantities. A series of order-of-magnitude analyses about these quantities are conducted. Consequently, the relative motion equations are approximated by expansions truncated to the second order, i.e. square of the dimensionless distance. In order to study the problem of periodicity of relative motion, the semi-major axis of the follower is expanded as Taylor series around that of the leader, by regarding relative position and velocity as small quantities. Using this expansion, it is proved that the periodicity condition derived from Lawden’s equations is equivalent to the condition that the Taylor series of order one is zero. The first-order relative motion equations, simplified from the second-order ones, possess the same forms as the periodic solutions of Lawden’s equations. It is presented that the latter are further first-order approximations to the former; and moreover, compared with the latter more suitable to research spacecraft rendezvous and docking, the former are more suitable to research relative orbit configurations. The first-order relative motion equations are expanded as trigonometric series with eccentric anomaly as the angle variable. Except the terms of order one, the trigonometric series’ amplitudes are geometric series, and corresponding phases are constant both in the radial and in-track directions. When the trajectory of the in-plane relative motion is similar to an ellipse, a method to seek this ellipse is presented. The advantage of this method is shown by an example.     

Exact Delaunay normalization of the perturbed Keplerian Hamiltonian with tesseral harmonics

A novel approach for the exact Delaunay normalization of the perturbed Keplerian Hamiltonian with tesseral and sectorial spherical harmonics is presented in this work. It is shown that the exact solution for the Delaunay normalization can be reduced to quadratures by the application of Deprit’s Lie-transform-based perturbation method. Two different series representations of the quadratures, one in powers of the eccentricity and the other in powers of the ratio of the Earth’s angular velocity to the satellite’s mean motion, are derived. The latter series representation produces expressions for the short-period variations that are similar to those obtained from the conventional method of relegation. Alternatively, the quadratures can be evaluated numerically, resulting in more compact expressions for the short-period variations that are valid for an elliptic orbit with an arbitrary value of the eccentricity. Using the proposed methodology for the Delaunay normalization, generalized expressions for the short-period variations of the equinoctial orbital elements, valid for an arbitrary tesseral or sectorial harmonic, are derived. The result is a compact unified artificial satellite theory for the sub-synchronous and super-synchronous orbit regimes, which is nonsingular for the resonant orbits, and is closed-form in the eccentricity as well. The accuracy of the proposed theory is validated by comparison with numerical orbit propagations.     

Minimum control for spacecraft formations in a J2 perturbed environment

Large ΔV amounts are often required to maintain the relative geometry which is needed to implement a formation flying concept. A wise use of the orbital environment makes the orbit keeping phase easier, reducing the overall propellant consumption. A first important step in this direction is the selection of formation configurations and orbits which, while still satisfying the mission requirements, are less subject to perturbations resulting naturally in closed relative motion. Within this frame, a number of studies have been recently carried out in order to identify possible sets of invariant relative orbits under the effects of the Earth oblateness, a conservative force commonly referred to as J2 which is also the most important perturbation for Low Earth Orbit. These efforts clearly marked the difficulties connected with achieving genuine periodic relative motion under J2 effect, but they also showed the existence of a set of conditions on the orbital parameters which allow for quasi-periodic J2 trajectories. This paper presents these particular trajectories, by means of deeper theoretical explanations, showing the dependency of the shape of the relative configurations on the orbital inclination. Since the quasi-periodic trajectories cannot be written analytically, and moreover, they are very sensitive with respect to the initial conditions, difficulties arise when trying to exploit these paths as reference for the control of a formation. This paper proposes a novel approach to find, from the actual quasi periodic natural trajectories, minimal control periodic reference trajectories. Next, it evaluates quantitatively the amount of propellant which is needed to control a spacecraft formation along these paths. The choice of Hill’s classical circular projected configuration as a nominal trajectory is considered as a comparison, showing the clear advantages of the proposed guidance design, which assumes low-perturbed periodic reference orbits as nominal trajectories.     

Kepler discretization in regular celestial mechanics

In this paper, a special extrapolation method for the numerical integration of perturbed Kepler problems (given in KS-formulation) is worked out and analyzed in detail. The underlying so-called Kepler discretization isexact for the pure (elliptic) Kepler motion. A numerically stable realization is presented together with a backward error analysis: this analysis shows that the effect of the arising rounding errors can be regarded as a small perturbation inferior to the physical perturbation. For test purposes, a well-known example describing the motion of an artificial Earth satellite in an equator plane subject to the oblateness perturbation is used to demonstrate the efficiency of the new extrapolation method.Proceedings of the Sixth Conference on Mathematical Methods in Celestial Mechanics held at Oberwolfach (West Germany) from 14 to 19 August, 1978.     

Density models for the upper atmosphere

Modeling the effects of atmospheric drag is one of the more important problems associated with the determination of the orbit of a near-earth satellite. Errors in the drag model can lead to significant errors in the determination and prediction of the satellite motion. The uncertainty in the drag acceleration can be attributed to three separate effects: (a) errors in the atmospheric density model, (b) errors in the ballistic coefficient, and (c) errors in the satellite relative velocity. In a number of contemporary satellite missions, the requirements for performing the orbit determination and predictions in near real-time has placed an emphasis on density model computation time as well as the model accuracy. In this investigation, a comparison is made of three contemporary atmospheric density models which are candidates for meeting the current orbit computation requirements. The models considered are the analytic Jacchia-Roberts model, the modified Harris-Priester model, and the USSR Cosmos satellite derived density model. The computational characteristics of each of the models are compared and a modification to the modified Harris-Priester model is proposed which improves its ability to represent the diurnal variation in the atmospheric density.This investigation was supported by the NASA Goddard Spaceflight Center under contract NAS5-20946 and Contract NSG 5154.     

A Perturbative Treatment of The Co-Orbital Motion

We develop a formalism of the non-singular evaluation of the disturbing function and its derivatives with respect to the canonical variables. We apply this formalism to the case of the perturbed motion of a massless body orbiting the central body (Sun) with a period equal to that of the perturbing (planetary) body. This situation is known as the co-orbital motion, or equivalently, as the 1/1 mean motion commensurability. Jupiter's Trojan asteroids, Earth's co-orbital asteroids (e.g., (3753) Cruithne, (3362) Khufu), Mars' co-orbital asteroids (e.g., (5261) Eureka), and some Jupiter-family comets are examples of the co-orbital bodies in our solar system. Other examples are known in the satellite systems of the giant planets. Unlike the classical expansions of the disturbing function, our formalism is valid for any values of eccentricities and inclinations of the perturbed and perturbing body. The perturbation theory is used to compute the main features of the co-orbital dynamics in three approximations of the general three-body model: the planar-circular, planar-elliptic, and spatial-circular models. We develop a new perturbation scheme, which allows us to treat cases where the classical perturbation treatment fails. We show how the families of the tadpole, horseshoe, retrograde satellite and compound orbits vary with the eccentricity and inclination of the small body, and compute them also for the eccentricity of the perturbing body corresponding to a largely eccentric exoplanet's orbit.This revised version was published online in October 2005 with corrections to the Cover Date.     

Lagrange’s planetary equations for the motion of electrostatically charged spacecraft

A spacecraft that generates an electrostatic charge on its surface in a planetary magnetic field will be subject to a perturbative Lorentz force. Active modulation of the surface charge can take advantage of this electromagnetic perturbation to modify or to do work on the spacecraft’s orbit. Lagrange’s planetary equations are derived using the Lorentz force as the perturbation on a Keplerian orbit, incorporating orbital inclination and true anomaly for the first time for an electrostatically charged vehicle. The planetary equations reveal that orbital inclination is a second-order effect on the perturbation, explaining results found in earlier studies through numerical integration. All of the orbital elements are coupled, but the coupling notably does not depend on the magnitude of the electrostatic charge or on the strength of the magnetic field. Analytical expressions that characterize this coupling are tested with a propellantless escape example at Jupiter. A closed-form solution exists that constrains the set of equatorial orbits for which planetary escape is possible, and a sufficient condition is identified for escape from inclined orbits. The analytical solutions agree with results from the numerically integrated equations of motion to within a fraction of a percent.     

Obliquity variations of terrestrial planets in habitable zones

We have investigated obliquity variations of possible terrestrial planets in habitable zones (HZs) perturbed by a giant planet(s) in extrasolar planetary systems. All the extrasolar planets so far discovered are inferred to be jovian-type gas giants. However, terrestrial planets could also exist in extrasolar planetary systems. In order for life, in particular for land-based life, to evolve and survive on a possible terrestrial planet in an HZ, small obliquity variations of the planet may be required in addition to its orbital stability, because large obliquity variations would cause significant climate change. It is known that large obliquity variations are caused by spin-orbit resonances where the precession frequency of the planet's spin nearly coincides with one of the precession frequencies of the ascending node of the planet's orbit. Using analytical expressions, we evaluated the obliquity variations of terrestrial planets with prograde spins in HZs. We found that the obliquity of terrestrial planets suffers large variations when the giant planet's orbit is separated by several Hill radii from an edge of the HZ, in which the orbits of the terrestrial planets in the HZ are marginally stable. Applying these results to the known extrasolar planetary systems, we found that about half of these systems can have terrestrial planets with small obliquity variations (smaller than 10°) over their entire HZs. However, the systems with both small obliquity variations and stable orbits in their HZs are only 1/5 of known systems. Most such systems are comprised of short-period giant planets. If additional planets are found in the known planetary systems, they generally tend to enhance the obliquity variations. On the other hand, if a large/close satellite exists, it significantly enhances the precession rate of the spin axis of a terrestrial planet and is likely to reduce the obliquity variations of the planet. Moreover, if a terrestrial planet is in a retrograde spin state, the spin-orbit resonance does not occur. Retrograde spin, or a large/close satellite might be essential for land-based life to survive on a terrestrial planet in an HZ.