深空探测器定位在共线平动点附近的线性策略

深空探测器需要定位在日-地(月)系的共线平动点L1或L2附近执行探测任务,但由于共线平动点的不稳定性,必须在运行期间进行轨控。对于条件周期轨道(如晕轨道)必须在控制过程中考虑高次项,控制条件复杂,技术上实现相对困难。而某些探测任务,探测器定位在共线平动点附近的条件拟周期轨道(对应L issajous轨道)上亦可以。这种类型的轨道可以离共线平动点较近,那么只需要在控制过程中考虑线性项即可,控制条件简单。以圆型限制性三体问题作为基本模型,采用预估-校正法逼近线性化模型下的目标轨道,给出在轨运行期间的轨控策略亦是可取的,这种控制措施相对而言较简单,容易实现。仿真计算结果表明是可行的,能够提供较高的位置精度。     

RCD太阳帆航天器人工平动点轨道保持与控制

针对航天器平动点轨道保持问题,研究了含有反射率控制设备(RCD)的太阳帆航天器在日地系共线人工平动点处的轨道保持与控制,同时降低因频繁改变航天器姿态所带来的振动问题。首先,基于太阳帆圆型限制性三体问题,计算了RCD型太阳帆人工平动点位置,给出了太阳帆共线人工平动点三阶Halo轨道,并将其作为参考轨道;然后,将太阳帆动力学方程线性化,采用跟踪控制输出的方法对线性模型进行控制;最后,通过合理选择控制变量矩阵,将控制律代入非线性模型中进行轨道保持控制。仿真结果表明,通过控制RCD太阳帆反射率设备参数及姿态角,实现了长时间的Halo轨道保持,同时大幅减小了太阳帆姿态角的改变,从而减小了帆面振动,为太阳帆航天器长期轨道任务的实现提供了良好的理论依据。     

低轨卫星对高轨卫星仅测角初轨计算方法

针对我国地面测站对高轨卫星监视能力缺乏的问题,提出一种低轨卫星对高轨卫星仅测角初轨计算方法。该算法引入天基跟踪坐标系,消除测距信息影响,建立仅测角观测方程;引入法向运动,增加摄动因素影响,建立扩展拉普拉斯动力学模型;推导分析观测模型和动力学模型的关系方程,将初轨计算问题转换为非线性方程求解问题,利用高斯全主元消去法完成方程求解。通过实战和仿真测角数据对方法进行检验,结果表明,该方法能利用仅测角数据对非合作目标进行初轨确定,精度在公里量级,可为我国地基监视系统提供补充参考。     

The ARTEMIS Mission

The Acceleration, Reconnection, Turbulence, and Electrodynamics of the Moon??s Interaction with the Sun (ARTEMIS) mission is a spin-off from NASA??s Medium-class Explorer (MIDEX) mission THEMIS, a five identical micro-satellite (hereafter termed ??probe??) constellation in high altitude Earth-orbit since 17 February 2007. By repositioning two of the five THEMIS probes (P1 and P2) in coordinated, lunar equatorial orbits, at distances of ??55?C65 R _E geocentric (??1.1?C12 R _L selenocentric), ARTEMIS will perform the first systematic, two-point observations of the distant magnetotail, the solar wind, and the lunar space and planetary environment. The primary heliophysics science objectives of the mission are to study from such unprecedented vantage points and inter-probe separations how particles are accelerated at reconnection sites and shocks, and how turbulence develops and evolves in Earth??s magnetotail and in the solar wind. Additionally, the mission will determine the structure, formation, refilling, and downstream evolution of the lunar wake and explore particle acceleration processes within it. ARTEMIS??s orbits and instrumentation will also address key lunar planetary science objectives: the evolution of lunar exospheric and sputtered ions, the origin of electric fields contributing to dust charging and circulation, the structure of the lunar interior as inferred by electromagnetic sounding, and the lunar surface properties as revealed by studies of crustal magnetism. ARTEMIS is synergistic with concurrent NASA missions LRO and LADEE and the anticipated deployment of the International Lunar Network. It is expected to be a key element in the NASA Heliophysics Great Observatory and to play an important role in international plans for lunar exploration.     

Finite time L1 Approach for Missile Overload Requirement Analysis in Terminal Guidance

This article analyzes the problem about the missile overload requirement in a homing terminal guidance under various engagement scenarios. An augmented proportional navigation guidance (APNG) model is introduced on the basis of linear kinematics. To analyze the peak-to-peak performance of the terminal guidance system, a new finite time L₁ performance measure for a linear time-varying (LTV) continuous system is proposed. Then, according to the idea of the adjoint system, a novel method for computing the L₁ norm of a linear continuous system is first derived. Within the finite time L₁ framework, the quantitative relation between the guidance loop dynamics and the maximum missile-target maneuver ratio is offered. This relation is expressed in the form of graphs and formulas that can be used to synthesize some of the major subsystem specifications for the missile guidance system. The illustrative examples show that a significant performance improvement is achieved with the proposed guidance loop dynamics.     

Artificial equilibrium points and bi-impulsive maneuvers to observe 243 Ida

In the restricted three-body problem, the traditional Lagrange points L₁ and L₂ are the only equilibrium points near the asteroid 243 Ida. The thrust generated by a solar sail over a spacecraft enables the existence of new artificial equilibrium points, which depend on the position of the spacecraft with respect to the asteroid and the attitude of the solar sail. Such equilibrium points generate new spots to observe the body from above or below the plane of motion. Such points are very good observational locations due to their stationary condition. This work provides a preliminary analysis to observe Ida through the use of artificial equilibrium points as spots combined with transfer maneuvers between them. Such combination can be used to observe the asteroid from more different points of view in comparison to fixed ones. The analyses are made for a spacecraft equipped with a solar sail and capable of performing bi-impulsive maneuvers. The solar radiation pressure is used both to maintain the equilibrium condition and to reduce the costs of the transfers and/or to create transfers with longer duration. This is a new aspect of the present research, because it combines the continuous thrust with initial and final small impulses, which are feasible for most of the spacecraft, because the magnitudes of the impulses are very low. These combined maneuvers may reduce the transfer times of the maneuvers in most of the cases, compared with the maneuvers based only on continuous thrust. Several options involved in these transfers are shown, like to minimize the fuel spent ($\mathrm{\Delta }\mathit{v}$) as a function of the transfer time or to extend the duration of the travel between the points. Extended transfer times can be useful when observations are required during the transfers.     

Nonlinear dynamics and station-keeping control of a rotating tethered satellite system in halo orbits

The dynamics of a rotating tethered satellite system （TSS） in the vicinity of libration points are highly nonlinear and inherently unstable. In order to fulfill the station-keep control of the rotating TSS along halo orbits, a nonlinear output tracking control scheme based on the θ- D technique is proposed. Compared with the popular time-variant linear quadratic regulator （LQR） controller, this approach overcomes some limitations such as on-line computations of the algebraic Riccati equation. Besides, the obtained nonlinear suboptimal controller is in a closed form and easy to implement. Numerical simulations show that the TTS trajectories track the periodic reference orbit with low energy consumption in the presence of both tether and initial injection errors. The axis of rotation can keep pointing to an inertial specific object to fulfill an observation mission. In addition, the thrusts required by the controller are in an acceptable range and can be implemented through some low-thrust propulsion devices.     

基于地月三角平动点的卫星自主定轨

近年提出了利用地月系平动点建立深空导航星座的设想。在受太阳摄动的真实力模型下,地月系平动点是不稳定的,从而会导致导航星座必须通过控制才能定点在特定区域。针对此问题,引入一种特殊的平动点轨道,即动力学替代轨道。平动点轨道卫星星座可利用星间测距数据自主定轨,由于动力学替代轨道具有长期稳定性,整个自主定轨过程不需要来自地面的测控支持,且定轨精度可达到观测精度。研究结果表明,观测资料的长短、导航卫星垂直白道面的运动分量都将影响到导航星间的自主定轨精度。该研究成果可以应用在以后的地月系导航星座中。     

描述弹塑性材料损伤演变行为的3种计算方法的研究与应用

借助于双向组合坐标系,就弹塑性材料的损伤演变速率和寿命提出了3种计算方法.它包括:单参数法、比值法、乘积法,同时本文给出了许多新的计算方程,还提供了一个同各个常用材料参数有着函数关系的综合性材料常数的新概念.这些常用参数有:疲劳强度系数σ'f、疲劳强度指数b'1、疲劳延性系数ε'f、疲劳延性指数c'1、平均应力、平均应变和临界加载时间等等.用汽车的一个零件做了实例计算,并对计算结果进行了一些比较.预计这些计算方法对避免过多的疲劳试验,对节省疲劳试验的时间、人力和资金,对方便工程应用在某种程度上将有着实际意义.     

混合小推力航天器日心悬浮轨道保持控制

针对太阳帆、太阳电混合推进航天器日心悬浮轨道保持控制问题进行了研究。为解决基于局部线性化模型设计轨道保持控制器时存在的控制精度不高、模型精确性过度依赖等问题,应用自抗扰控制(ADRC)技术设计了轨道保持控制器。首先,采用圆形限制性三体问题(CRTBP)模型推导了混合小推力航天器日心悬浮轨道动力学方程;然后,考虑系统模型不确定性和外部扰动,提出了一种基于扰动估计和补偿的轨道保持控制方法;最后,数值仿真表明存在系统模型不确定性、初始入轨误差及地球轨道偏心率扰动等因素的情况下,所设计的控制器仅需很小的速度增量即可实现高精度的日心悬浮轨道保持控制。     

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

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

Efficient aerodynamic shape optimization using variable-fidelity surrogate models and multilevel computational grids

A variable-fidelity method can remarkably improve the efficiency of a design optimization based on a high-fidelity and expensive numerical simulation, with assistance of lower-fidelity and cheaper simulation(s). However, most existing works only incorporate “two” levels of fidelity, and thus efficiency improvement is very limited. In order to reduce the number of high-fidelity simulations as many as possible, there is a strong need to extend it to three or more fidelities. This article proposes a novel variable-fidelity optimization approach with application to aerodynamic design. Its key ingredient is the theory and algorithm of a Multi-level Hierarchical Kriging (MHK), which is referred to as a surrogate model that can incorporate simulation data with arbitrary levels of fidelity. The high-fidelity model is defined as a CFD simulation using a fine grid and the lower-fidelity models are defined as the same CFD model but with coarser grids, which are determined through a grid convergence study. First, sampling shapes are selected for each level of fidelity via technique of Design of Experiments (DoE). Then, CFD simulations are conducted and the output data of varying fidelity is used to build initial MHK models for objective (e.g. C_D) and constraint (e.g. C_L, C_m) functions. Next, new samples are selected through infill-sampling criteria and the surrogate models are repetitively updated until a global optimum is found. The proposed method is validated by analytical test cases and applied to aerodynamic shape optimization of a NACA0012 airfoil and an ONERA M6 wing in transonic flows. The results confirm that the proposed method can significantly improve the optimization efficiency and apparently outperforms the existing single-fidelity or two-level-fidelity method.     

地月系共线平动点探测器的星上轨道预报问题

针对定点在地月系共线平动点附近运动的探测器，研究其轨道预报问题，并探讨了星上外推预报时的力模型简化问题。首先定量分析了共线平动点L1和L2附近探测器的受力情形；然后构造了真实力模型下这些点附近的目标轨道并研究了其误差传播规律，指出这类轨道相比一般卫星轨道预报的不同特征；最后结合星上有限的计算和存储资源限制探讨了力模型的简化问题。此项工作可为定点在平动点附近的探测器的星上预报提供参考。     

Nonlinear adaptive flight control system: Performance enhancement and validation

The problem of decreasing stability margins in L₁ adaptive control systems is discussed and an out-of-loop L₁ adaptive control scheme based on Lyapunov’s stability theorem is proposed. This scheme enhances the effectiveness of the adaptation, which ensures that the system has sufficient stability margins to achieve the desired performance under parametric uncertainty, additional delays, and actuator faults. The stability of the developed control system is demonstrated through a series of simulations. Compared with an existing control scheme, the constant adjustment of the stability margins by the proposed adaptive scheme allows their range to be extended by a factor of 4–5, bringing the stability margin close to that of variable gain PD control with adaptively scheduled gains. The engineered practicability of adaptive technology is verified. A series of flight tests verify the practicability of the designed adaptive technology. The results of these tests demonstrate the enhanced performance of the proposed control scheme with nonlinear parameter estimations under insufficient stability margins and validate its robustness in the event of actuator failures.     

Solar Magnetism: The State of Our Knowledge and Ignorance

We review some longstanding scientific mysteries related to solar magnetism, with final attention to the mystery of the “turbulent diffusion” essential for the theoretical α ω-dynamo that is believed to be the source of the magnetic fields of the Sun. Fundamental difficulties with the concept of turbulent diffusion of magnetic fields suggest that the solar dynamo problem needs to be reformulated. An alternative dynamo model is proposed, but it remains to be shown that the model can provide the quantitative aspects of the cyclic magnetic fields of the Sun.     

基于特征模型的Halo轨道维持控制

平动点，尤其是共线平动点轨道在未来深空探测活动中具有重要的应用价值，但由于共线平动点轨道不稳定，运行在其上的航天器在无控情况下将很快偏离标称轨道，因此在实际任务中，轨道维持必不可少。针对地月系L₂点附近的Halo轨道维持问题，首先在圆型限制性三体模型下，利用Richardson三阶近似解析解、微分修正以及打靶法获得了用于维持控制的标称轨道；然后设计了基于特征模型理论的黄金分割控制器用于速度跟踪以及PD控制器用于位置跟踪；最后分别在圆型限制性三体模型和双圆限制性四体模型下进行了仿真分析。结果表明：在两种模型下，位置和速度的跟踪精度分别优于100 m和0.003 m/s，但双圆限制性四体模型下所需总的速度增量比圆型限制性三体模型下所需总的速度增量高一个数量级。     

太阳辐射对ER3BP拉格朗日点的影响

为了研究Halo轨道中继通信以及行星际轨道转移等问题,需要首先对拉格朗日点(平动点)的理论进行分析。使用考虑太阳辐射的椭圆型限制性三体建模,代替原有的圆型限制性三体建模,提高了建模的精度。提出了在L1,L2,L3点的更为精确的新解析解,经过Matlab仿真,其相对于精确数值解的误差,分别为Ammar的解析解的13.13%,29.51%,0.46%。     

Stratosphere NO y Species Measured by MIPAS and GOMOS Onboard ENVISAT During 2002–2010: Influence of Plasma Processes onto the Observed Distribution and Variability

We present observations of stratosphere NO _y species from 2002 to 2010 taken by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS) and by the Global Ozone Monitoring by Occultation of Stars (GOMOS) instruments on board ENVISAT. We jointly used observations of MIPAS NO₂, HNO₃, N₂O₅, ClONO₂ and N₂O, and GOMOS NO₂ and NO₃. MIPAS results are part of the MIPAS2D database retrieved adopting a full 2D tomographic approach. We describe the mean distribution and variability of NO _y species in the stratosphere, identifying changes induced by plasma processes. Beside enhancements due to sporadic solar proton events, we show that winter polar NO₂ has an almost linear relationship with the geomagnetic activity index Ap down to about 10?hPa. This indicates a dominant role of energetic precipitating particles in the production of upper atmosphere NO _y . The correlation has clear signatures extending to mid latitudes. Partitioning of the NO _y reservoir species are also traced, with HNO₃ and N₂O₅ showing a correlation with Ap extending to lower altitude within the polar regions. We found no large signatures of an impact of thunderstorm-induced plasma processes onto monthly means of NO _y species in the stratosphere.     

Assessment and improvement of k-ω-γ model for separation-induced transition prediction

The purpose of this work is to improve the k-ω-γ transition model for separation-induced transition prediction. The fundamental cause of the excessively small separation bubble predicted by k-ω-γ model is scrutinized from the perspective of model construction. On the basis, three rectifications are conducted to improve the k-ω-γ model for separation-induced transition. Firstly, a damping function is established via comparing the molecular diffusion timescale with the rapid pressure-strain timescale. The damping function is applied to prevent the effective length scale from incorrect distribution near the leading edge of the separation bubble. Secondly, the pressure gradient parameter λ_ζ, is proposed as an indicator for local susceptibility to the separation instability. Additionally, λ_ζ,-based separation intermittency γ_sep is constructed to accelerate the substantial growth of turbulent kinetic energy after flow separation. The improved model appropriate for both low- and high-speed flow has been calibrated against a variety of diverse and challenging experiments, including the subsonic T3L plate, Aerospatial A airfoil, transonic NLR-7301 airfoil and deformed hypersonic inflatable aerodynamic decelerator aeroshell. The improved model is strictly based on local variables and Galilean invariance. Besides, the proposed improvement for k-ω-γ model can be fairly convenient to incorporate into other existing intermittency-based transition models.     

一种大椭圆轨道卫星星历预报方法

针对大偏心率轨道长时间星历预报,稠密输出星历的低效率问题,提出了一种新的星历预报方法。该方法通过建立卫星摄动力模型,对卫星运动方程进行数值积分来获取卫星星历。该方法的轨道积分器以标准的8阶Adams-Cowell多步法为基础,针对大偏心率轨道的特性,在一个轨道周期的不同时段内采用不同的积分步长,同时引入了用于生成小间隔等间距卫星星历的插值公式。该方法降低了卫星运动方程右函数的计算次数,尤其适用于需要稠密输出卫星星历的情况。以STK(Satellite Tool Kit,卫星工具包)的HPOP(High Precision Orbit Propagator,高精度轨道预报)模型为验证基准,通过数学仿真校验了该方法的有效性和精度。算例表明,该方法在预报精度损失很小的情况下,计算时间远小于标准的多步法和HPOP模型。