基于大幅值Lyapunov轨道的地月转移轨道设计研究

基于平面圆形限制性三体问题模型,利用与绕月轨道相切的大幅值Lyapunov周期轨道,提出了一种新的地月转移轨道设计方法。根据Poincaré截面与限制性三体问题动力学系统对称性计算得到的大幅值Lyapunov轨道,通过与绕月轨道拼接,将地月转移问题转化为地球到大幅值Lyapunov轨道的转移问题。为保证探测器能够从近地轨道（LEO）切向逃逸到达大幅值Lyapunov轨道,通过计算其稳定流形,采用最近点作为Poincaré截面的终止条件求解探测器的初始状态,并根据初始状态完成地月轨道的设计。仿真结果表明,该地月转移策略相比于Hohmann转移,在同样只需要两次速度增量的前提下,约节约100 m/s的速度增量,该研究为地月转移轨道的设计提供了一种新思路。     

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     

Satellite Skin-Force Modelling for Atmospheric Drag Calculations

Satellites in low Earth orbits are influenced by the Earth’s atmosphere. The interactions between the molecules and the spacecraft cause the highest non-gravitational force, which in magnitude is comparable to planetary disturbances. Therefore the modelling of atmospheric drag effects is important for many missions with a scientific background like STEP (Satellite Test of Equivalence Principle). With the STEP mission variations between gravitational and inertial mass shall be measured with an accuracy of 10^?18. The results are of great interest for cosmological and gravitational theories. To achieve the aimed accuracy, a precise model of external disturbances is necessary. In this article the method of Ray-Tracing is used to quantify the atmospheric drag forces and torques for spacecrafts of arbitrary shape.     

Simultaneous spacecraft attitude and orbit estimation using magnetic field vector measurements

Based on magnetometer measurements only, three-axis attitude, rate, and orbit estimation are successfully achieved. A single Augmented Dynamics Extended Kalman Filter (ADEKF) is configured by combining the spacecraft nonlinear attitude dynamics and quaternion kinematics with orbital mechanics. The filter design is adopted for three-axis stabilized spacecraft in low Earth orbits where the aerodynamic drag is the dominant source of disturbances in addition to the spacecraft magnetic residuals. To reduce the computational burden, another Interlaced Extended Kalman Filter (IEKF) is developed to uncouple the attitude/rate from the orbit dynamics. Both filters are implemented using the magnetometer measurements and their corresponding time derivatives. As a part of EgyptSat-1 flight scenario, detumbling and standby modes are used for performance testing of the ADEKF. The concept of local observability is applied to the basic filter and the stability is investigated by incorporating extensive Monte Carlo simulations with uniformly distributed initial conditions. The filter shows the capability of estimating the attitude better than 5 deg and rate of order 0.03 deg/s in each axis. In orbit estimation, the filter is capable of estimating the position with accuracy less than 8 km and velocity upto 5 m/s in each axis.     

Earth-moon Trajectory Optimization Using Solar Electric Propulsion

The optimization of the Earth-moon trajectory using solar electric propulsion is presented. A feasible method is proposed to optimize the transfer trajectory starting from a low Earth circular orbit (500 km altitude) to a low lunar circular orbit (200 km altitude). Due to the use of low-thrust solar electric propulsion, the entire transfer trajectory consists of hundreds or even thousands of orbital revolutions around the Earth and the moon. The Earth-orbit ascending (from low Earth orbit to high Earth orbit) and lunar descending (from high lunar orbit to low lunar orbit) trajectories in the presence of J2 perturbations and shadowing effect are computed by an analytic orbital averaging technique. A direct/indirect method is used to optimize the control steering for the trans-lunar trajectory segment, a segment from a high Earth orbit to a high lunar orbit, with a fixed thrust-coast-thrust engine sequence. For the trans-lunar trajectory segment, the equations of motion are expressed in the inertial coordinates about the Earth and the moon using a set of nonsingular equinoctial elements inclusive of the gravitational forces of the sun, the Earth, and the moon. By way of the analytic orbital averaging technique and the direct/indirect method, the Earth-moon transfer problem is converted to a parameter optimization problem, and the entire transfer trajectory is formulated and optimized in the form of a single nonlinear optimization problem with a small number of variables and constraints. Finally, an example of an Earth-moon transfer trajectory using solar electric propulsion is demonstrated.     

Two-maneuver indirect contingency return from a low lunar orbit

The problem of contingency return from the low lunar orbit is studied. A novel two-maneuver indirect return strategy is proposed. By effectively using the Earth’s gravity to change the orbital plane of the transfer orbit, the second maneuver in the well-known three-maneuver return strategy can be removed, so the total delta-v is reduced. Compared with the single-maneuver direct return, our strategy has the advantage in that the re-entry epoch for the minimum delta-v cost can be advanced in time, with a minimum delta-v value similar to that of the direct return. The most obvious difference between our strategy and the traditional single- or multiple- maneuver strategies is that the complete transfer orbit is a patch between a two-body conic orbit and a three-body orbit instead of two conic orbits. Our strategy can serve as a useful option for contingency return from a low lunar orbit, especially when the delta-v constraint is stringent for a direct return and the contingency epoch is far away from the return window.     

MARC-a system for simulation and visualization of space missionscenarios

MARC (modeling, animation, rendering, and compositing), a system using advanced computer graphics and animation techniques for spacecraft mission simulation, is described. The MARC system provides capabilities for generating complex models of both man-made and natural phenomena. The system models orbital dynamics of terrestrial satellites, supports solid models for the Earth, Sun, and Moon, and simulates the dynamics of terrestrial satellites for arbitrary elliptical orbits. A stellar background including magnitudes and spectral types is generated. The elements of the MARC system, including object modeling tools, orbital animation techniques, the rendering system used to compute individual frames, and the compositing techniques used, are discussed. The software architecture of the MARC system and the hardware used to support the system are described     

基于STK的小卫星轨道交会设计研究

利用球面三角形的几何关系进行了基于STK(卫星工具包)的小卫星轨道交会的规划和设计。采用一种较为方便的轨道交会设计方法,并使用功能强大的专业STK软件来仿真和演示,取得了良好的仿真效果,能够较好地满足航天器轨道交会设计的要求,对于其它类型的轨道规划有一定的借鉴作用,为各种航天器的仿真研究提供了一种新的方法。     

A Detailed Analysis of the Operational Orbit of the Hipparcos Satellite

An analysis of the orbital evolution of the ESA's Hipparcos satellite is presented. Hipparcos operated between August 1989 and March 1993 in a highly elliptical orbit: a geostationary transfer orbit with increased perigee height. The requirements of the scientific mission included high accuracy knowledge of the position and velocity vectors of the spacecraft as a function of time. Through a study of the variations in the total orbital energy, the loss of energy during the mission as a result of non-conservative forces is recovered. These are explained as largely due to atmospheric drag during perigee passages. Apparent variations in the drag coefficient are in agreement with orientation variations of the satellite during those perigee passages. Two different models used for calculating the atmospheric drag give significantly different results, confirming earlier findings by other users of those models. This revised version was published online in August 2006 with corrections to the Cover Date.     

Science Objectives and Rationale for the Radiation Belt Storm Probes Mission

The NASA Radiation Belt Storm Probes (RBSP) mission addresses how populations of high energy charged particles are created, vary, and evolve in space environments, and specifically within Earth’s magnetically trapped radiation belts. RBSP, with a nominal launch date of August 2012, comprises two spacecraft making in situ measurements for at least 2 years in nearly the same highly elliptical, low inclination orbits (1.1×5.8 RE, 10^°). The orbits are slightly different so that 1 spacecraft laps the other spacecraft about every 2.5 months, allowing separation of spatial from temporal effects over spatial scales ranging from ～0.1 to 5 RE. The uniquely comprehensive suite of instruments, identical on the two spacecraft, measures all of the particle (electrons, ions, ion composition), fields (E and B), and wave distributions (d E and d B) that are needed to resolve the most critical science questions. Here we summarize the high level science objectives for the RBSP mission, provide historical background on studies of Earth and planetary radiation belts, present examples of the most compelling scientific mysteries of the radiation belts, present the mission design of the RBSP mission that targets these mysteries and objectives, present the observation and measurement requirements for the mission, and introduce the instrumentation that will deliver these measurements. This paper references and is followed by a number of companion papers that describe the details of the RBSP mission, spacecraft, and instruments.     

The STEREO Mission: An Introduction

The twin STEREO spacecraft were launched on October 26, 2006, at 00:52 UT from Kennedy Space Center aboard a Delta 7925 launch vehicle. After a series of highly eccentric Earth orbits with apogees beyond the moon, each spacecraft used close flybys of the moon to escape into orbits about the Sun near 1 AU. Once in heliospheric orbit, one spacecraft trails Earth while the other leads. As viewed from the Sun, the two spacecraft separate at approximately 44 to 45 degrees per year. The purposes of the STEREO Mission are to understand the causes and mechanisms of coronal mass ejection (CME) initiation and to follow the propagation of CMEs through the inner heliosphere to Earth. Researchers will use STEREO measurements to study the mechanisms and sites of energetic particle acceleration and to develop three-dimensional (3-D) time-dependent models of the magnetic topology, temperature, density and velocity of the solar wind between the Sun and Earth. To accomplish these goals, each STEREO spacecraft is equipped with an almost identical set of optical, radio and in situ particles and fields instruments provided by U.S. and European investigators. The SECCHI suite of instruments includes two white light coronagraphs, an extreme ultraviolet imager and two heliospheric white light imagers which track CMEs out to 1 AU. The IMPACT suite of instruments measures in situ solar wind electrons, energetic electrons, protons and heavier ions. IMPACT also includes a magnetometer to measure the in situ magnetic field strength and direction. The PLASTIC instrument measures the composition of heavy ions in the ambient plasma as well as protons and alpha particles. The S/WAVES instrument uses radio waves to track the location of CME-driven shocks and the 3-D topology of open field lines along which flow particles produced by solar flares. Each of the four instrument packages produce a small real-time stream of selected data for purposes of predicting space weather events at Earth. NOAA forecasters at the Space Environment Center and others will use these data in their space weather forecasting and their resultant products will be widely used throughout the world. In addition to the four instrument teams, there is substantial participation by modeling and theory oriented teams. All STEREO data are freely available through individual Web sites at the four Principal Investigator institutions as well as at the STEREO Science Center located at NASA Goddard Space Flight Center.     

IBEX Backgrounds and Signal-to-Noise Ratio

The Interstellar Boundary Explorer (IBEX) mission will provide maps of energetic neutral atoms (ENAs) originating from the boundary region of our heliosphere. On IBEX there are two sensors, IBEX-Lo and IBEX-Hi, covering the energy ranges from 10 to 2000 eV and from 300 to 6000 eV, respectively. The expected ENA signals at 1 AU are low, therefore both sensors feature large geometric factors. In addition, special attention has to be paid to the various sources of background that may interfere with our measurement. Because IBEX orbits the Earth, ion, electron, and ENA populations of the Earth’s magnetosphere are prime background sources. Another potential background source is the magnetosheath and the solar wind plasma when the spacecraft is outside the magnetosphere. UV light from the night sky and the geocorona have to be considered as background sources as well. Finally background sources within each of the sensors must be examined.     

一种基于轨道根数约束的最优制导方法

针对航天器空间变轨的制导问题,研究了一种基于轨道根数约束的最优制导方法。在地心惯性坐标系下直接建立航天器的最优控制模型,给出了位置速度表达式和协态变量初值之间的关系;进一步,在不约束真近点角的前提下,推导了5个轨道根数的约束方程,并通过对轨道根数和最优控制理论中协态方程的特性分析,获得了另外两个约束方程。协态变量初值可直接通过求解7个完整约束方程组获得,进而得到最优推力方向。仿真验证了所提制导方法的有效性。     

Gravitation and celestial mechanics investigations with Galileo

The gravitation and celestial mechanics investigations during the cruise phase and Orbiter phase of the Galileo mission depend on Doppler and ranging measurements generated by the Deep Space Network (DSN) at its three spacecraft tracking sites in California, Australia, and Spain. Other investigations which also rely on DSN data, and which like ours fall under the general discipline of spacecraft radio science, are described in a companion paper by Howard et al. (1992). We group our investigations into four broad categories as follows: (1) the determination of the gravity fields of Jupiter and its four major satellites during the orbital tour, (2) a search for gravitational radiation as evidenced by perturbations to the coherent Doppler link between the spacecraft and Earth, (3) the mathematical modeling, and by implication tests, of general relativistic effects on the Doppler and ranging data during both cruise and orbiter phases, and (4) an improvement in the ephemeris of Jupiter by means of spacecraft ranging during the Orbiter phase. The gravity fields are accessible because of their effects on the spacecraft motion, determined primarily from the Doppler data. For the Galilean satellites we will determine second degree and order gravity harmonics that will yield new information on the central condensation and likely composition of material within these giant satellites (Hubbard and Anderson, 1978). The search for gravitational radiation is being conducted in cruise for periods of 40 days centered around solar opposition. During these times the radio link is least affected by scintillations introduced by solar plasma. Our sensitivity to the amplitude of sinusoidal signals approaches 10^-15 in a band of gravitational frequencies between 10^-4 and 10^-3 Hz, by far the best sensitivity obtained in this band to date. In addition to the primary objectives of our investigations, we discuss two secondary objectives: the determination of a range fix on Venus during the flyby on 10 February, 1990, and the determination of the Earth's mass (GM) from the two Earth gravity assists, EGA1 in December 1990 and EGA2 in December 1992.     

Conceptual design for an orbiting debris detection radar

There exists a large amount of man-made debris in low Earth orbit. The quantity of this debris is growing every year as a result of on-going activities in space. Much of the debris consists of particles which are too small to be tracked from the ground, but nevertheless pose a threat to spacecraft. This paper examines the possibility of an orbiting radar which, over a period of time, could measure debris positions and velocities and thereby build up a database of debris orbits. Spacecraft could use this database for advance warning of impending collisions and maneuver themselves out of the collision path, thereby mitigating the long-term risk of collision damage     

The contribution of observations of satellites to the determination of the Earth's gravitational potential

The history of the determination of the external gravitational potential of the Earth is sketched briefly. A discussion of the principles by which the potential may be derived from the observations of changes in the orbits of artificial satellites is followed by outlines of the principal theories and by detailed consideration of the formal differences between them that arise from differences in the ways that the orbits are described and it is shown that those formulae on which most of the numerical results depend are equivalent in the principal theories. The usual methods of treatment break down in certain special conditions and the analysis of these cases is also considered although they are not of great practical importance in the derivation of the potential; similarly a short account is given of the behaviour of a satellite having an orbital angular velocity commensurable with the spin angular velocity of the Earth. Methods by which satellites are observed are mentioned and the main numerical results on the external potential of the Earth are discussed critically. Finally the results are compared with those derived from observations of gravity on the surface of the Earth and the application of the results to problems of geodesy, of the physical state of the Earth and of the motion of the Moon are described.     

航天器多普勒测速平均误差的分析与修正

对航天器多普勒测速平均误差进行了分析,详细推导了圆轨道和椭圆轨道时该误差的计算公式,计算了不同采样周期和轨道高度时的误差大小。经过分析指出,该项误差对于高轨航天器影响较小,对于低轨航天器可以通过缩短采样周期或利用3个点或多个点的连续测量数据进行修正。     

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     

连续地月转移系统动力学研究与能量分析

为了研究新型连续地月转移系统的动力学及能量需求,采用Lagrange方法,在系绳为刚性杆假设的前提下,同时忽略第三体引力、地球扁率和系绳轴向变形等扰动因素的影响,建立了驱动型动量交换绳系卫星(MMET)系统的三维刚性动力学模型。对所建立的动力学模型进行了数值仿真及对比分析,仿真结果验证了所建模型的正确性。研究表明,外力矩对系统轨道运动参数影响甚小,对姿态运动参数影响明显。采用MMET方式进行载荷转移,推导出了实现载荷地月轨道转移所需的入口速度条件以及时间周期条件,并求解出了载荷在2次任务之间的时间间隔。给定初始条件下,当MMET系统以0.231 6 rad/s的旋转角速度绕其质心旋转1 448.5圈,其绕地心刚好运行5圈时,载荷可顺利进入地月转移轨道。最后,对连续地月转移系统实现载荷的地月转移进行了能量对比分析,结果表明,相同条件下,MMET载荷转移方式相比于传统脉冲变轨方式在载荷转移过程中消耗更少的能量。     

绕飞慢自旋小天体的航天器运动分析

针对可以用三轴椭球体近似建模的小天体,给出了非球形引力势函数,建立了航天器绕飞小天体的轨道动力学方程。利用Jacobi积分常数绘制了探测器在小天体周围的零速度曲线,并分析了探测器的可能运动区域,给出了航天器不碰撞小天体的边界条件。针对绕飞慢自旋小天体的情况,基于平均轨道根数的近似解分析了小天体扁率和椭率的摄动影响,并给出了几条冻结轨道及其稳定条件。