火星轨道交会自主导航与制导方法

对火星采样返回任务中的火星轨道交会自主导航和制导技术进行了研究。采用光学自主导航敏感器测量的火星中心方向和视半径，相对敏感器测量的相对位置等观测量，设计了导航滤波器同时估计轨返组合体和上升器的轨道。在导航滤波器设计中，针对光学自主导航敏感器更新频率远低于滤波解算频率的问题，设计了一种连续观测量构造算法，确保每个滤波周期均可进行测量更新，以提高导航精度。基于导航滤波器估计结果，采用T-H制导设计了4脉冲共椭圆交会策略实施轨道控制，从而构成近程交会自主导航和制导方案用于完成火星轨道交会任务。通过数学仿真校验了所提出方法的有效性。     

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

同步双小行星系统共振轨道设计

研究了同步双小行星系统中共振轨道的设计方法及演化规律。首先，基于双椭球模型建立探测器运动方程，并给出共振轨道初值选取方法。然后，利用改进并行打靶法，提出一种双小行星系统平面共振轨道两步修正方法。同时结合稳定性理论及分岔理论，给出双小行星系统三维共振轨道生成和延拓方法；最后，以双小行星系统1999KW4为例，设计了共振比为1∶1,1∶2,1∶3,1∶4,2∶3的平面和空间共振轨道族，并分析了共振轨道的特性及轨道周期和轨道能量的变化规律。给出的双小行星系统中共振轨道的设计方法具有普适性，对未来双小行星系统探测任务中的轨道设计具有一定的参考意义与借鉴价值。     

基于TLE的弹道系数计算方法及应用分析

针对仅使用两行要素（Two Line Element，TLE）作为数据源的应用需求，研究了基于TLE轨道衰减的弹道系数计算方法。介绍了一种常用的基于两组TLE的直接计算法，分析TLE选取间隔对结果精度的影响；提出了一种基于多组TLE的迭代计算方法，以降低异常TLE对计算结果的影响；从弹道系数计算效果、在再入预报中的应用等方面对这两种方法进行比较分析。结果表明，两种方法各有优劣，基于多组TLE的迭代计算法稳定性更高、受TLE精度的影响更小；由于数据区间更短，基于两组TLE的计算结果对短期轨道衰减特性反应得更准确，用于临近再入时的预报效果更好。     

地球静止轨道微波辐射计技术

静止轨道微波辐射计可实现全天时、全天候、高频次云雨大气观测,在台风、流域性强降水探测和预报方面发挥重要作用。但在波束扫描、多频复用、系统定标等方面的运用存在较大困难,世界上尚无在轨应用先例。现从探测需求出发,对静止轨道微波辐射计进行了系统方案设计。提出了卫星平台与辐射计部件快慢结合扫描的方法,以解决波束扫描和系统定标难题;提出了低频选入射角层叠式布局准光学频段分离方法,以解决多频段复用难题。研制了微波辐射计原理样机,通过实验室和外场试验验证了系统设计,为静止轨道微波探测卫星研制提供了一定指导。     

Orbit determination and prediction of GEO satellite of BeiDou during repositioning maneuver

In order to establish a continuous GEO satellite orbit during repositioning maneuvers, a suitable maneuver force model has been established associated with an optimal orbit determination method and strategy. A continuous increasing acceleration is established by constructing a constant force that is equivalent to the pulse force, with the mass of the satellite decreasing throughout maneuver. This acceleration can be added to other accelerations, such as solar radiation, to obtain the continuous acceleration of the satellite. The orbit determination method and strategy are illuminated, with subsequent assessment of the orbit being determined and predicted accordingly. The orbit of the GEO satellite during repositioning maneuver can be determined and predicted by using C-Band pseudo-range observations of the BeiDou GEO satellite with COSPAR ID 2010-001A in 2011 and 2012. The results indicate that observations before maneuver do affect orbit determination and prediction, and should therefore be selected appropriately. A more precise orbit and prediction can be obtained compared to common short arc methods when observations starting 1 day prior the maneuver and 2 h after the maneuver are adopted in POD (Precise Orbit Determination). The achieved URE (User Range Error) under non-consideration of satellite clock errors is better than 2 m within the first 2 h after maneuver, and less than 3 m for further 2 h of orbit prediction.     

Effects of in-track maneuver on Sun illumination conditions of near-circular low Earth orbits

By developing approximate analytical models considering the J₂ perturbation, the effects of an in-track maneuver on the orbital Sun illumination conditions of near-circular low Earth orbits are analyzed. First, two approximate models for the variations in orbital sunshine angles are developed, one for variations at a given time and the other for variations at a given argument of latitude. Next, two approximate models for variations in orbital arc in Earth shadow are developed, one considers the small eccentricity and the other uses the zero eccentricity. Finally, the developed approximate models are applied to analyzing the Sun illumination conditions of a typical in-track maneuver mission on a near-circular low Earth orbit. From the results obtained, three major conclusions can be drawn. First, the variations in orbital sunshine angles at a given time may reach tens of degrees when the drifting time reaches hundreds of orbital periods, and the approximate model for that situation cannot effectively approach the numerical results. Second, the variations in orbital sunshine angles for any given argument of latitude are only a couple of degrees even when the drifting time reaches 500 orbital periods, and the approximation model developed can effectively approach the numerical results. Third, for variations in orbital arc in Earth shadow, the approximate model considering the small eccentricity has simple expressions and can effectively approach the numerical results; in contrast, the approximate model using the zero eccentricity has relatively worse precision.     

用两颗GPS卫星确定低轨卫星轨道的初值方法

研究了用两颗全球定位系统（ＧＰＳ）卫星粗略确定低轨卫星轨道的问题,并着重讨论了其中的初值问题,给出了一种有效的初值方法。实验结果证明：这种方法能给出一个确定的、实用的初值,大大减少了初轨计算的计算量,并使得算法收敛更加有保障。在工程实用方面比传统的猜测方法具有较大的优势。     

原子氧对航天器表面材料作用的数值模拟

对氏地球轨道环境对航天器表面材料的影响和航天器表面材料的低地球轨道环境寿命评定方法进行了介绍;并对低地球轨道环境和地面试验环境下,有无保护涂层的业胺所受冲蚀作用进行了成功的数值模拟,获得了具有工程应用价值的数值计算结果。该工作对航天器的设计具有指导意义。     

速率控制下的空间系绳的伸展

以地球同步轨道高度上的空间系绳的伸展为例,阐述了长度速率控制律中的参数对伸展过程的影响,控制律中参数选择的方法,最大系绳拉力的估算以及伸展过程中系绳应力的计算。