关于地球转移轨道碎片的轨道寿命问题

低轨卫星的轨道寿命主要取决于大气的耗散作用,其轨道在不断变小(即高度降低)变圆的状态下进入地球稠密大气层中陨落。但地球转移轨道(GTO)碎片的运行轨道是一个近地点高度为200km,远地点高度达36000km的大偏心率(e=0.73)椭圆轨道,其轨道寿命主要由第三体(日、月)引力摄动所决定,而且还与其轨道的初始状态有密切关系。本文将根据地球卫星轨道变化规律进行理论分析,阐明力学机制,并给出相应的数值验证。     

大气密度模型用于近地卫星定轨预报的比较

大气阻力是低轨卫星主要的摄动力,与高层大气密度的变化密切相关。由于目前对高层大气密度变化的机制尚未完全掌握,所使用的各种大气密度模型多属于半经验公式。在这些模型中并没有一种在任何情况下都是最好的,因此,对于特定轨道选择合适的大气密度模型对提高定轨预报的精度是非常重要的。通过对资源2号卫星实测GPS数据的分析计算,比较了常用的8种大气密度模型的定轨预报精度,探讨了预报24h应采用的定轨数据长度和大气密度模型。     

Modelling the Torques Affecting the Hipparcos Satellite

Reconstructed attitude data for the Hipparcos mission as obtained in the final stages of the data analysis for the published catalogue is used to derive detailed information on the dynamics of the satellite. Most elements of the inertia tensor of the satellite could be calibrated from the observed acceleration data, which are also used to reconstruct torques due to solar radiation and gravity gradient, and the magnetic moment of the satellite and it's interaction with the magnetic field surounding the Earth. The effects of the oblateness of the Earth on the gravity gradient are evaluated and shown to be negligable. The magnetic field model includes both the `main' and the `disturbance' fields. The remaining systematic effects in residual torques are most likely attributed to variations in the magnetic field that are local and are beyond the models used to describe it. The angular momentum vector for one of the gyros was reconstructed from the torque it asserted on the satellite while it was running in redundant mode. This revised version was published online in August 2006 with corrections to the Cover Date.     

椭圆轨道星座构型自稳定设计研究

椭圆轨道星座构型稳定性要求同时实现升交点赤经、近地点幅角和平近点角的稳定。通过初始偏差对椭圆轨道卫星的长期影响分析可知,轨道半长轴、倾角和偏心率的初始偏差对升交点赤经、近地点幅角和平近点角的长期摄动变化是线性的,因此通过主动偏置轨道半长轴、偏心率和倾角能够实现椭圆轨道星座构型的自稳定设计,从而提高其构型稳定性。对实例星座的设计结果表明,自稳定设计方法对提高椭圆轨道星座的构型稳定性是有效的。     

The Operational Environment of the Hipparcos Mission

We investigate links between the observational environment as experienced by the Hipparcos satellite and the performance of the spacecraft and payload instrumentation, with particular emphasis on finding out whether some of these effects may have been inadequately represented in instrument calibrations and could thus have affected the scientific results of the mission. Scan-coverage and radiation effects are primarily random effects with only some long-term systematics. However, long- (days to weeks) and short-term (hours) temperature variations reflected in the performance of some of the spacecraft instrumentation. It is shown that only a small sign of some long-term thermal variations could be detected in the payload instrumentation. These findings further limit the scope left for the occurrence of large-scale correlated errors in the Hipparcos astrometric data. On the other hand, a number of great circles were identified which showed a highly significant drift of the basic angle, which had not been detected in the preparation of the published data. The data from these circles may have, in some cases, led to, very localised, slightly anomalous results, in particular where stars are accidentally affected by two or more of such circles. This revised version was published online in August 2006 with corrections to the Cover Date.     

An Overview of the Fast Auroral SnapshoT (FAST) Satellite

The FAST satellite is a highly sophisticated scientific satellite designed to carry out in situ measurements of acceleration physics and related plasma processes associated with the Earth's aurora. Initiated and conceptualized by scientists at the University of California at Berkeley, this satellite is the second of NASA's Small Explorer Satellite program designed to carry out small, highly focused, scientific investigations. FAST was launched on August 21, 1996 into a high inclination (83°) elliptical orbit with apogee and perigee altitudes of 4175 km and 350 km, respectively. The spacecraft design was tailored to take high-resolution data samples (or `snapshots') only while it crosses the auroral zones, which are latitudinally narrow sectors that encircle the polar regions of the Earth. The scientific instruments include energetic electron and ion electrostatic analyzers, an energetic ion instrument that distinguishes ion mass, and vector DC and wave electric and magnetic field instruments. A state-of-the-art flight computer (or instrument data processing unit) includes programmable processors that trigger the burst data collection when interesting physical phenomena are encountered and stores these data in a 1 Gbit solid-state memory for telemetry to the Earth at later times. The spacecraft incorporates a light, efficient, and highly innovative design, which blends proven sub-system concepts with the overall scientific instrument and mission requirements. The result is a new breed of space physics mission that gathers unprecedented fields and particles observations that are continuous and uninterrupted by spin effects. In this and other ways, the FAST mission represents a dramatic advance over previous auroral satellites. This paper describes the overall FAST mission, including a discussion of the spacecraft design parameters and philosophy, the FAST orbit, instrument and data acquisition systems, and mission operations.     

Investigating suitable orbits for the Swarm constellation mission – The frozen orbit

This paper proposes a suitable orbit design for the lower pair of ESA's Swarm constellation mission, flying side-by-side in near-polar and circular orbits with a separation of only 1.4° at ascending node. Both orbits are suggested to be frozen orbits to minimize the evolution, and an along-track separation strategy is applied to avoid collision risk. The characteristics of the proposed orbit type are examined through numerical techniques including high-fidelity perturbation models. The prime change from the initial configuration is an along-track separation. The perturbations causing the along-track drift are analyzed by switching on/off certain perturbations. The results indicate that the tesseral harmonics and the atmospheric drag yield dominant effects. The atmospheric drag effect shows a dependence on the local time of the ascending node. From two months of orbit propagation for the altitude 300 km the maximum along-track drift we obtain is about 80 km, which is still within the measurement requirement range. Several maneuver strategies for maintaining the proposed orbit design are suggested. The results analyzed for the proposed orbit design show that collision risk can be avoided by along-track separation within the frozen orbit design. Consequently, this combination is considered as a suitable approach for Swarm's lower pair.     

Introduction to a Further Examination of the Hipparcos Data

We present an introduction to four papers on further analysis of the raw Hipparcos data. This analysis has lead to the recognition of how orbit, radiation and temperature conditions did or didn't affect the scientific results of the Hipparcos mission. It also led the way to a new reduction of the scientific data that shows from its initial results a real potential for a substantial improvement of the astrometric parameters for stars brighter than V=8.5. This short paper serves as an introduction to the four main papers, and provides some general references. This revised version was published online in August 2006 with corrections to the Cover Date.     

Estimation of ballistic coefficients of space debris using the ratios between different objects

This paper proposes a new method to estimate the ballistic coefficient (BC) of low earth orbit space debris.The data sources are the historical two-line elements (TLEs).Since the secular variation of semi-major axes is mainly caused by the drag perturbation for space objects with perigee altitude below 600 km,the ballistic coefficients are estimated based on variation of the mean semi-major axes derived from the TLEs.However,the approximate parameters used in the calcu lation have error,especially when the upper atmosphere densities are difficult to obtain and always estimated by empirical model.The proportional errors of the approximate parameters are cancelled out in the form of ratios,greatly mitigating the effects of model error.This method has been also been validated for space objects with perigee altitude higher than 600 km.The relative errors of esti mated BC values from the new method are significantly smaller than those from the direct estimation methods used in numerical experiments.The estimated BC values are used for the prediction of the semi-major axes,and good performance is obtained.This process is also a feasible method for prediction over a long period of time without an orbital propagator model.     

Tidal Models in a New Era of Satellite Gravimetry

The high precision gravity measurements to be made by recently launched (and recently approved) satellites place new demands on models of Earth, atmospheric, and oceanic tides. The latter is the most problematic. The ocean tides induce variations in the Earth's geoid by amounts that far exceed the new satellite sensitivities, and tidal models must be used to correct for this. Two methods are used here to determine the standard errors in current ocean tide models. At long wavelengths these errors exceed the sensitivity of the GRACE mission. Tidal errors will not prevent the new satellite missions from improving our knowledge of the geopotential by orders of magnitude, but the errors may well contaminate GRACE estimates of temporal variations in gravity. Solar tides are especially problematic because of their long alias periods. The satellite data may be used to improve tidal models once a sufficiently long time series is obtained. Improvements in the long-wavelength components of lunar tides are especially promising. This revised version was published online in August 2006 with corrections to the Cover Date.     

摄动对编队飞行星座相对构型的影响分析

在近圆轨道编队飞行的假设条件下，根据动力学关系推导出了环绕卫星相对参考卫星的运动学简化模型，并以此简化模型为基础，分别研究了大气阻力摄动、日月引力摄动、太阳光压摄动和地球扁率摄动对编队飞行星座的构型影响，着重分析了地球扁率摄动周期项和长期项对构型的影响，并以此对编队轨道设计提出建议。     

Dynamic Modelling of the Hipparcos Attitude

We present a new method for a high-accuracy reconstruction of the attitude for a slowly spinning satellite. This method, referred to as the fully-dynamic approach, explores the possibility to describe the satellite's attitude as that of a rigid body subject to continuous external torques. The method is tried out on the Hipparcos data and is shown to reduce the noise for the along-scan attitude reconstruction for that mission by about a factor two to three. The dynamic modelling is expected to give a more accurate representation of the satellite's attitude than was obtained with a pure mathematical modelling. As such, it decreases the degrees of freedom in the a posteriori reconstruction. Some of the decrease is obtained through accumulating and subsequently implementing information on high frequency components in the solar radiation torques, which show to be systematic and predictable. This could be expected, as they are primarily linked to the external geometry and optical properties of the satellite. In the context of an astrometric mission, the methods presented here can only be applied as a final iteration step: the star positions that are used to reconstruct the attitude are also part of the scientific objectives of the mission. An estimate for the potential of a re-reduction of the Hipparcos data using the fully-dynamic model for the attitude reconstruction was obtained from test reductions of the first 24 months of mission data. Improvement of the accuracies of the astrometric parameters for all stars brighter than Hp=9.0 appears possible. The noise on the astrometric parameters for these stars was affected significantly by the along-scan attitude noise, which dominated for stars brighter than Hp=4.5. The possible improvement for stars brighter than about Hp=4.5 may, after iterations, be as much as a factor three. The reduced noise levels also allow a more accurate calibration and monitoring of instrument parameters, leading potentially to a better understanding of the instrument and the scientific data obtained with it. This revised version was published online in August 2006 with corrections to the Cover Date.     

Minimax Approach to Trajectory Optimization of Multistage Launch Vehicles

An optimum pitch steering program is developed, using a minimax technique, for a multistage launch vehicle in the presence of large parameter uncertainties. The uncertainties are characterized by deterministic bounds on the respective parameters. The pitch steering program is obtained by maximizing two independent scalar performance ance indexes. 1) the coast apogee velocity for a specified altitude; 2) the perigee of the satellite orbit. The product of dynamic pressure and angle of attack is constrained so as to minimize the structural loads during the atmospheric flight. The values of the uncertain parameters are determined by minimizing the same performance indexes in order to achieve a worst case design. The existence of saddle point solution to this class of problems is shown using the techniques of differential game theory. The conjugate gradient algorithm has been used for computer er aided design of the minimax technique.     

磁暴下的LEO卫星精密轨道确定

对低轨卫星（LEO）,大气阻尼摄动是主要的定轨误差源.尤其在发生磁暴时,求解一个大气阻尼因子的定轨方法已不能充分吸收大气密度计算不准所造成的定轨误差,因而在标校统一S波段（USB）的测量系统差和随机差时往往计算失真.本文提出了一种求解折线型Cd因子的新方法,克服了动力学模型不准所带来的定轨误差,通过与独立的GPS数据比较,定轨精度有明显提高,同时给出的测量系统差和随机差更加真实可信.     

Fuel-optimal trajectories of aeroassisted orbital transfer withplane change

The fuel-optimal control problem arising in noncoplanar orbital transfer employing aeroassist technology is addressed. The mission involves the transfer from high Earth orbit to low Earth orbit with plane change. The complete maneuver consists of a deorbit impulse to inject a vehicle from a circular orbit to an elliptic orbit for atmospheric entry, a boost impulse at the exit from the atmosphere for the vehicle to attain a desired orbital altitude, and a reorbit impulse to circularize the path of the vehicle. In order to minimize the total fuel consumption, a performance index is chosen as the sum of the deorbit, boost, and reorbit impulses. The application of optimization principles leads to a nonlinear, two-point, boundary value problem, which is solved by a multiple shooting method     

Galileo trajectory design

The Galileo spacecraft was launched by the Space Shuttle Atlantis on October 18, 1989. A two-stage Inertial Upper Stage propelled Galileo out of Earth parking orbit to begin its 6-year interplanetary transfer to Jupiter. Galileo has already received two gravity assists: from Venus on February 10, 1990 and from Earth on December 8, 1990. After a second gravity-assist flyby of Earth on December 8, 1992, Galileo will have achieved the energy necessary to reach Jupiter. Galileo's interplanetary trajectory includes a close flyby of asteroid 951-Gaspra on October 29, 1991, and, depending on propellant availability and other factors, there may be a second asteroid flyby of 243-Ida on August 28, 1993. Upon arrival at Jupiter on December 7, 1995, the Galileo Orbiter will relay data back to Earth from an atmospheric Probe which is released five months earlier. For about 75 min, data is transmitted to the Orbiter from the Probe as it descends on a parachute to a pressure depth of 20–30 bars in the Jovian atmosphere. Shortly after the end of Probe relay, the Orbiter ignites its rocket motor to insert into orbit about Jupiter. The orbital phase of the mission, referred to as the satellite tour, lasts nearly two years, during which time Galileo will complete 10 orbits about Jupiter. On each of these orbits, there will be a close encounter with one of the three outermost Galilean satellites (Europa, Ganymede, and Callisto). The gravity assist from each satellite is designed to target the spacecraft to the next encounter with minimal expenditure of propellant. The nominal mission is scheduled to end in October 1997 when the Orbiter enters Jupiter's magnetotail.List of Acronyms ASI Atmospheric Structure Instrument - EPI Energetic Particles Instrument - HGA High Gain Antenna - IUS Inertial Upper Stage - JOI Jupiter Orbit Insertion - JPL Jet Propulsion Laboratory - LRD Lightning and Radio Emissions Detector - NASA National Aeronautics and Space Administration - NEP Nephelometer - NIMS Near-Infrared Mapping Spectrometer - ODM Orbit Deflection Maneuver - OTM Orbit Trim Maneuver - PJR Perijove Raise Maneuver - PM Propellant Margin - PDT Pacific Daylight Time - PST Pacific Standard Time - RPM Retropropulsion Module - RRA Radio Relay Antenna - SSI Solid State Imaging - TCM Trajectory Correction Maneuver - UTC Universal Time Coordinated - UVS Ultraviolet Spectrometer - VEEGA Venus-Earth-Earth Gravity Assist     

Spacecraft orbit determination using GPS navigation solutions

The orbit determination using the GPS navigation solutions for the KOMPSAT-1 spacecraft has been studied. The Cowell method of special perturbation theories was employed to develop a precision orbit propagation, and the perturbations due to geopotential, the gravity of the Sun and the Moon, solid Earth tides, ocean tides, the Earth's dynamic polar motion, solar radiation pressure, and atmospheric drag were modeled. Specifically, the satellite box-wing macro model was applied to minimize the drag errors at low altitude. The estimation scheme consisted of an extended Kalman filter and Bayesian least square method. To investigate the applicability of the method to the KOMPSAT-1 spacecraft, the orbit determination was accomplished using the GPS navigation solutions for the TOPEX/POSEIDON and TAOS satellites. The orbit determination results were compared with NASA POE generated by global laser tracking. The position and velocity accuracy was estimated about 16∼7 m and 0.0157∼0.0074 m·s⁻¹ RMS, respectively, for the two satellites in the presence of SA. These results verify that an orbit determination scheme using GPS navigation solutions can provide the static orbit information and reduce conspicuously the position and velocity errors of navigation solutions. It can be suggested that the sequential and batch orbit determination using the GPS navigation solutions be the most appropriate method in the KOMPSAT-1 type mission.