常值径向推力下的飞行器运动轨迹

研究了初始停泊轨道为椭圆时,空间飞行器在常值径向推力下运动的有界性和周期性.首先建立了飞行器运动的动力学方程,并通过能量积分和角动量积分进行了简化.然后将有界性的研究转化为一个一元三次不等式的求解,并在此基础上针对不同的初始真近点角分别进行了研究,得到了运动的边界和有界性条件.接下来利用椭圆积分研究了运动的周期性,分别研究了径向运动、极角转动以及整体运动的周期性.最后用数值算法得到了运动的周期轨道.     

Mapping orbits around the asteroid 2001SN263

The present paper has the goal of mapping orbits, with respect to the perturbations, for a spacecraft traveling around the asteroid 2001SN₂₆₃. This asteroid is a triple system, which center of mass is in an elliptic orbit around the Sun. The perturbations considered in the present model are the ones due to the oblateness of the central body, the gravity field of the two satellite bodies (Beta and Gamma), the Sun, the Moon, the asteroids Vesta, Pallas and Ceres and all the planets of the Solar System. This mapping is important, because it shows the relative importance of each force for a given orbit for the spacecraft, helping to make a decision about which forces need to be included in the model for a given accuracy and nominal orbit. Another important application of this type of mapping is to find orbits that are less perturbed, since it is expected that those orbits have good potential to require a smaller number of station-keeping maneuvers. Simulations under different conditions are made to find those orbits. The main reason to study those trajectories is that, currently, there are several institutions in Brazil studying the possibility to make a mission to send a spacecraft to this asteroid (the so-called ASTER mission), because there are many important scientific studies that can be performed in that system. The results showed that Gamma is the main perturbing body, followed by Beta (10 times smaller) and the group Sun–Mars-oblateness of Alpha, with perturbations 1000 times weaker than the effects of Gamma. The other bodies have perturbations 10⁷ times smaller. The results also showed that circular and polar orbits are less perturbed, when compared to elliptical and equatorial orbits. Regarding the semi-major axis, an internal orbit is the best choice, followed by a larger external orbit. The inclination of the orbit plays an important role, and there are values for the inclination where the perturbations show minimum and maximum values, so it is important to make a good decision on those values.     

零动量卫星轮控系统可重构性设计研究

针对卫星姿态控制系统设计中常见的动量轮配置问题,总结了星上常用的5种轮控系统构型,首先对各种构型从角动量包络、干扰力矩下动量轮的饱和情况、系统功耗和可靠性等方面进行了可重构性分析和比较,在此基础上提炼了动量轮可重构性设计准则,可为轮控系统的设计提供参考依据。     

Shape optimization of a three-dimensional membrane-structured solar sail using an angular momentum unloading strategy

A shape of the satellite’s solar sail membrane is essential for unloading angular momentum in the three-axis stabilized attitude control system because the three-dimensional solar sail can receive solar radiation pressure from arbitrary directions. In this paper, the objective is the shape optimization of a three-dimensional membrane-structured solar sail using the angular momentum unloading strategy. We modelled and simulated the solar radiation pressure torque, for unloading angular momentum. Using the simulation system, since the unloading angular momentum rate is maximized, the shape of the three-dimensional solar sail was optimized using a Genetic algorithm and Sequential Quadratic Programming. The unloading velocity in the optimized shaped solar sail was greatly improved with respect to a conventional flat or pyramid solar sail.     

Jitter reduction of a reaction wheel by management of angular momentum using magnetic torquers in nano- and micro-satellites

Nowadays, nano- and micro-satellites, which are smaller than conventional large satellites, provide access to space to many satellite developers, and they are attracting interest as an application of space development because development is possible over shorter time period at a lower cost. In most of these nano- and micro-satellite missions, the satellites generally must meet strict attitude requirements for obtaining scientific data under strict constraints of power consumption, space, and weight. In many satellite missions, the jitter of a reaction wheel degrades the performance of the mission detectors and attitude sensors; therefore, jitter should be controlled or isolated to reduce its effect on sensor devices. In conventional standard-sized satellites, tip-tilt mirrors (TTMs) and isolators are used for controlling or isolating the vibrations from reaction wheels; however, it is difficult to use these devices for nano- and micro-satellite missions under the strict power, space, and mass constraints. In this research, the jitter of reaction wheels is reduced by using accurate sensors, small reaction wheels, and slow rotation frequency reaction wheel instead of TTMs and isolators. The objective of a reaction wheel in many satellite missions is the management of the satellite’s angular momentum, which increases because of attitude disturbances. If the magnitude of the disturbance is reduced in orbit or on the ground, the magnitude of the angular momentum that the reaction wheels gain from attitude disturbances in orbit becomes smaller; therefore, satellites can stabilize their attitude using only smaller reaction wheels or slow rotation speed, which cause relatively smaller vibration. In nano- and micro-satellite missions, the dominant attitude disturbance is a magnetic torque, which can be cancelled by using magnetic actuators. With the magnetic compensation, the satellite reduces the angular momentum that the reaction wheels gain, and therefore, satellites do not require large reaction wheels and higher rotation speed, which cause jitter. As a result, the satellite can reduce the effect of jitter without using conventional isolators and TTMs. Hence, the satellites can achieve precise attitude control under low power, space, and mass constraints using this proposed method. Through the example of an astronomical observation mission using nano- and micro-satellites, it is demonstrated that the jitter reduction using small reaction wheels is feasible in nano- and micro-satellites.     

Solar particle dynamics during magnetic storms of July 23–27, 2004

It is a case study of a chain of three magnetic storms with a special attention to the particle dynamics based on CORONAS-F and SERVIS-1 low altitude satellite measurements. Solar proton penetration inside the polar cap and inner magnetosphere and dynamics at different phases of the magnetic storms was studied. We found, that solar protons were captured to the inner radiation belt at the recovery phase of the first and the second magnetic storms and additionally accelerated during the last one. No evidence of sudden commencement (SC) particle injection was found. Enhanced solar proton belt intensity with small pitch angles decreased slowly during satellite orbits for 30 days until the next magnetic storm. Then in 20–30 h we registered strong precipitation of these protons followed by the trapped proton flux dropout. Intensity decrease was more pronounced at lower altitudes and higher particle energies.     

Estimating the spin axis orientation of the Echostar-2 box-wing geosynchronous satellite

For the first time, the spin axis orientation of an inactive box-wing geosynchronous satellite has been estimated from ground-based optical photometric observations of Echostar-2’s specular reflections. Recent photometric light curves obtained of Echostar-2 over four years suggest that unusually bright and brief specular reflections were occurring twice within an observed spin period. These bright and brief specular reflections suggested two satellite surfaces with surface normals separated by approximately 180°. The geometry between the satellite, the Sun, and the observing location at the time of each of the brightest observed reflections, was used to estimate Echostar-2’s equatorial spin axis orientation coordinates. When considering prograde and retrograde rotation, Echostar-2’s spin axis orientation was estimated to have been located within 30° of either equatorial coordinate pole. Echostar-2’s spin axis was observed to have moved approximately 180° in right ascension, within a time span of six months, suggesting a roughly one year spin axis precession period about the satellite’s angular momentum vector.     

A study of reaction wheel configurations for a 3-axis satellite attitude control

The satellite reaction wheel’s configuration plays also an important role in providing the attitude control torques. Several configurations based on three or four reaction wheels are investigated in order to identify the most suitable orientation that consumes a minimum power. Such information in a coherent form is not summarized in any publication; and therefore, an extensive literature search is required to obtain these results. In addition, most of the available results are from different test conditions; hence, making them difficult for comparison purposes. In this work, the standard reaction wheel control and angular momentum unloading schemes are adopted for all the reaction wheel configurations. The schemes will be presented together with their governing equations, making them fully amenable to numerical treatments. Numerical simulations are then performed for all the possible reaction wheel configurations with respect to an identical reference mission. All the configurations are analyzed in terms of their torques, momentums and attitude control performances. Based on the simulations, the reaction wheel configuration that has a minimum total control torque level is identified, which also corresponds to the configuration with minimum power consumption.     

太阳同步轨道卫星电源系统设计计算方法研究

给出了太阳同步轨道卫星电源系统的两种基本结构框架 ;分析了这两种结构的电源系统的设计计算方法 ;推导出了基于单圈能量平衡和多圈能量平衡的计算公式。以一个算例验证了该算法 ,并计算出所需太阳电池阵的面积和蓄电池组的容量 ,最后进行了能量平衡分析 ,得到了单圈或多圈能量平衡的结果。     

Near real-time BDS GEO satellite orbit determination and maneuver analysis with reversed point positioning

The Geostationary Earth Orbit (GEO) satellite is a crucial part of the BeiDou Navigation Satellite System (BDS) constellation. However, due to various perturbation forces acting on the GEO satellite, it drifts gradually over time. Thus, frequent orbit maneuvers are required to maintain the satellite at its designed position. During the orbit maneuver and recovery periods, the orbit quality of the maneuvered satellite computed with broadcast navigation ephemeris will be significantly degraded. Furthermore, the conventional dynamic Precise Orbit Determination (POD) approach may not work well, because of a lack of publicly available satellite information for modeling the thrust forces. In this paper, a near real-time approach free of thrust forces modeling is proposed for BDS GEO satellite orbit determination and maneuver analysis based on the Reversed Point Positioning (RPP). First, the station coordinates and receiver clock offsets are estimated by GPS/BDS combined Single Point Positioning (SPP) with single-frequency phase-smoothed pseudorange observations. Then, with the fixed station coordinates and receiver clock offsets, the RPP method can be conducted to determine the GEO satellite orbits. When no orbit maneuvers occur, the proposed method can obtain orbit accuracies of 0.92, 2.74, and 8.30?m in the radial, along-track, and cross-track directions, respectively. The average orbit-only Signal-In-Space Range Error (SISRE) is 1.23?m, which is slightly poorer than that of the broadcast navigation ephemeris. Using four days of GEO maneuvered datasets, it is further demonstrated that the derived orbits can be employed to characterize the behaviors of GEO satellite maneuvers, such as the time span of the maneuver as well as the satellite thrusting accelerations. These results prove the efficiency of the proposed method for near real-time GEO satellite orbit determination during maneuvers.     

空间绳系拖拽系统摆动特性与平稳控制

考虑了任务星与废星的姿态运动以及系统组合体的面内外姿态运动,建立了绳系拖拽离轨系统动力学与控制模型,以切向常值推力下绳系拖拽轨道转移为任务过程,分析了任务星在喷气和零动量轮的限制姿态反馈控制条件下飞行时,废星姿态摆动、系统组合体面内外摆动和任务星姿态运动的规律及相互影响关系。采用留位和阻尼控制相结合的系绳张力复合控制方法,并结合任务星姿态控制,确保绳系拖拽转移安全平稳进行。仿真结果表明:常值推力下绳系拖拽轨道转移时,牵挂点偏置诱发的废星姿态周期性摆动会激发绳系组合体的面内外同频率高阶摆动,星体姿态运动是任务星姿态扰动力矩产生的主要因素;采用张力复合控制可有效消除废星姿态摆动并保持星间相对距离,结合任务星姿态控制,可实现离轨过程的平稳与安全,大幅减少任务星的姿控能耗。      

Old tail lobes effect on the solar-wind – Magnetosphere energy transport for the 27 August 2001 substorm

The magnetic flux of tail lobes Ψ is divided in two parts of comparable values Ψ₁ and Ψ₀₂, with the first that appears during substorm and the second, observed before substorm start. The first was named “new magnetic flux”, the second – “old magnetic flux”. The first, Ψ₁, is known to play a definitive role in the energy transport from the solar wind into the magnetosphere-ionosphere-atmosphere system, but the role of Ψ₀₂ in this transport is not well known. From the 27 August 2001 substorm data we study the involvement in the above transport process of the old flux Ψ₀₂. This involvement is observed in the polar cap (PC) area, which existed prior to the substorm and is called respectively “the old PC”. In this study, as distinct from earlier works, we use the balance equation of the energy stored in magnetosphere and energy consumed. Activation of the old PC magnetic flux Ψ₀₂ was found to increase the energy input by ∼85% in the event under consideration.     

Characteristics of cosmic ray cutoffs for a satellite orbiting at 400 km; The effect of the solid earth

Studies of galactic cosmic ray intensity and composition by means of earth orbiting space vehicles often rely on the earth's magnetic field as a momentum analyzer or threshold energy indicator. We have determined, by the trajectory-tracing process, cosmic ray cutoff characteristics for a satellite at 400 km altitude. These calculations indicate that cosmic rays have direct access to the satellite from the magnetic west down to angles of 135° from the zenith. For spacecraft observations it is necessary to employ additional definitions of cutoff beyond the classical Stormer definition because there are regions with allowed orbits where the cutoff, in the classical sense, does not exist.     

A BRIEF INTRODUCTION TO GEOSPACE DOUBLE STAR PROGRAM

The Geospace Double Star Exploration Project (DSP) contains two small satel lites operating in the near-earth equatorial and polar regions respectively. The tasks of DSP are: (1) to provide high-resolution field, particle and wave mea surements in several important near-earth active regions which have not been covered by existing ISTP missions, such as the near-earth plasma sheet and its boundary layer, the ring current, the radiation belts, the dayside magnetopause boundary layer, and the polar region; (2) to investigate he trigger mechanisms of magnetic storms, magnetospheric substorms, and magnetospheric particle storms, as well as the responses of geospace storms to solar activities and in terplanetary disturbances; (3) to set up the models describing the spatial and temporal variations of the near-earth space environment.To complete the mission, there are eight instruments on board the equatorial satellite and the polar satellite, respectively. The orbit of the equatorial satellite is proposed with a perigee at 550km and an apogee at 60 000km, and the inclination is about 28.5°; while the orbit of the polar satellite with a perigee at 700 km and an apogee at 40 000 km, as well as an inclination about 90°. The equatorial and polar satellites are planed to be launched into orbits in June 2003 and December 2003 respectively to take coordinating measurements with Cluster Ⅱ and other missions.     

H-reversal trajectory: A new mission application for high-performance solar sails

The aim of this paper is to quantify the performance of a flat solar sail to perform a double angular momentum reversal maneuver and produce a new class of two-dimensional, non-Keplerian orbits in the ecliptic plane. For a given pair of orbital parameters, the orbital period and the perihelion distance, it is possible to find the minimum solar sail characteristic acceleration required to fulfil a double angular momentum reversal trajectory. This problem is addressed using an optimal formulation and is solved through an indirect approach. The new trajectories are symmetrical with respect to the sun-perihelion line and exhibit a bean-like shape. Two main difficulties must be properly taken into account. On one side the sail is required to perform a rapid reorientation maneuver when it approaches the perihelion. Suitable simulations have shown that such a maneuver is feasible. In the second place the new trajectories require the use of high performance solar sails. For example, assuming an orbital period equal to 5 years, the required solar sail characteristic acceleration is greater than 3.4 mm/s². Such a value, although beyond the currently available sail performance, is comparable to what is required by the original concept of H-reversal maneuvers introduced by Vulpetti in 1996.     

On the structure and dynamics of the thermosphere

Thermospheric temperature, composition and wind measurements from the Dynamics Explorer satellite (DE-2) are interpreted using a three dimensional, multiconstituent spectral model. The analysis accounts for tides driven by the absorbed solar radiation as well as energy and momentum coupling involving the magnetosphere and lower atmosphere. We discuss phenomena associated with the annual tide, polar circulation, magnetic storms and substorms.     

Radiative forces and LAGEOS' orbit

The orbit analysis of LAGEOS satellite has resulted into the discovery and/or reassessment of several very small sources of perturbation on satellite orbits. The analysis of orbital arcs of duration ranging from one month to several years has revealed that perturbative effects are present, having unpredicted long-term or secular components down to the 10⁻¹² m/s² acceleration level. It was soon realized that those perturbations have a non-gravitational origin.

In recent years, we have devoted some effort to the physical modelling of radiative perturbations, caused by momentum exchanges with an appropriate radiative field, and have considered their potential role in the analysis of LAGEOS orbit residuals. These perturbations include: (i) direct solar radiation pressure; (ii) radiation pressure from the Earth's emitted/reflected/diffused radiation flux; (iii) the so-called thermal thrust force.

The main results of this work are reviewed, discussing its relationships with models developed by other research groups. In particular, we present a list of the physical processes which still appear to need more detailed and realistic modelling to reach a better understanding of LAGEOS dynamics at the 10⁻¹² m/s² level.     

小卫星动平衡研究--质心不对中对姿态的影响

首先定性分析质心不对中对星体动力学特性的影响，然后由拉格朗日第二类方程导出引力场下星体运动方程的矩阵表达式，最后以50kg级长方体3轴稳定小卫星为例，分析不同初角速度和不同时间历程情况下质心不对中对星体运动姿态的影响，并分析其对星体运动稳定性的影响．本文研究为同类品种小卫星的构型设计、姿态控制系统设计以及制定卫星动平衡精度提供依据。     

一种用于SGCMG奇异规避的CVP轨迹规划

采用单框架控制力矩陀螺(SGCMG)作为执行机构的小型敏捷卫星在姿态机动过程中存在着奇异问题.本文从SGCMG姿态控制系统整体出发，将奇异问题转化为状态约束的动态控制问题，基于控制变量参数化(CVP)方法，设计了一种用于SGCMG奇异规避的轨迹规划.该算法在实现小型敏捷卫星大角度姿态机动过程无奇异的基础上，将SGCMG框架角转速的最优轨迹通过CVP方法进行分段线性规划.这种规划策略对框架伺服系统的算法设计无复杂要求，仅需要简单的加减速控制，从而节约了星上资源.在轨迹规划实现过程中，考虑了工程实际中的约束条件，可以按照姿态机动任务要求规划出一条综合考虑能量资源和目标精度的最优轨迹.仿真结果表明：该算法实现了姿态参数轨迹和星体角速度轨迹的平缓变化，目标误差在1×10-3量级，星体在机动过程中运行稳定，SGCMG不会出现奇异现象.