地球卫星动力学定轨中摄动模型的选取

地球卫星动力学定轨中需要依据卫星轨道高度考虑不同的摄动项,以提高定轨精度。针对这一需求,分析了不同轨道高度下卫星所受的地球非球形、日月引力、太阳光压和大气阻力4种摄动项对二体模型的校正量,数值仿真表明实验结果与摄动项的力学分析相吻合,另外,给出了地球卫星动力学定轨中不同轨道下卫星摄动项的选取原则。     

高阶次月球重力场模型截断下的低轨环月卫星轨道受摄分析

研究低轨月球卫星在月球非球形摄动和地球第三体引力摄动作用下轨道高度变化问题.首先依据Kaula准则比较分析目前国际上公认的最精确的两个重力场模型GLGM-2和LP165P,提出了在一定阶次截断重力场模型的问题,然后通过仿真不同阶次重力场模型作用下轨道高度为50km的圆形极轨道环月卫星轨道特征的变化,验证了 50km以上高度卫星非球形摄动分析时可以将重力场模型截断至一定阶次的结论,并利用截断至70阶次的重力场模型仿真得到了50km和200km圆轨道卫星无控条件下正常运行的时间.最后在仿真地球引力对200km圆轨道卫星高度影响的基础E仿真其在月球非球形和地球引力摄动作用下轨道要素变化,对低轨环月卫星轨道保持控制提供依据.     

Sagnac效应对星间时间同步的影响

文章介绍了星间双向时间比对原理,分析了Sagnac效应对双向时间同步的影响,推出了地球非旋转坐标系同轨和异轨卫星星间时间同步的Sagnac效应公式,并以静止同步轨道为例,进行了分析计算。结果表明,对于同轨卫星,星间时间同步的Sagnac影响与轨道半径的平方根以及卫星与地心夹角的正弦成正比。     

多约束条件下月亮对IGSO/MEO卫星干扰预报方法

针对IGSO/MEO卫星倾斜轨道(55 °)和偏航主动姿态控制的特殊性,研究了适用倾斜中高轨卫星月球干扰的亮度模型、地球遮挡模型、几何满足模型,提出了多元素约束条件下的月球干扰预报方法,解决了传统预报算法虚警概率较大的问题.实际控制结果表明,该方法实现了对在轨IGSO/MEO卫星地球敏感器月球干扰的准确预报,取得了很好的应用效果.     

J2摄动作用下近地轨道卫星编队构形长期演化机理分析

研究近地轨道卫星编队构形在引力摄动作用下的长期演化机理。把用转移矩阵法得出的相对运动分析解表达成一般轨道根数形式，为研究摄动影响而进一步表达成无奇点轨道根数形式；根据各轨道根数的J2摄动解分别研究编队构形在轨道坐标系三个方向上的振幅与相位的长期变化情况，揭示了构形引力摄动的长期演化机理；给出编队构形的摄动表达式，并以仿真实例验证其可信性；最后给出构形长期引力摄动的相关结论，并对编队构形设计提出建议，这将为研究编队构形的规划技术带来方便。     

GPS卫星光压摄动建模发展与启示

对中高轨道的导航卫星而言,光压摄动模型误差是目前导航卫星动力学模型最大的误差源,因此,需要对光压摄动进行分析并精确建模。文章按照光压涉动建模原理的不同,研究分析了GPS卫星光压模型随工程需求的发展过程、建模方法的改进以及最新的建模方法,比较了不同建模方法下光压摄动模型的工程应用特点,提出了建立适用于北斗导航卫星姿态控制规律的高精度光压摄动模型有待研究的问题,以上内容可为北斗卫星精密定轨与轨道预报工作提供借鉴。     

基于星历拟合的短弧运动学定轨

当导航卫星在姿轨控和轨道恢复期间,传统的统计定轨理论难以实现精密定轨。首次提出 了基于10参数星历拟合的短弧运动学定轨方法,建立和推导了相应的理论模型和定轨解算方 法。其优点在于不仅能够反映卫星运动的物理学特征,提高了速度和轨道预报精度,而且不 需要累积数据,实现近实时快速计算,克服了动力学法定轨发散和单点定位无法获得速度信 息的不足。对COMPASS M-01导航卫星实测数据的处理表明,10分钟短弧运动学定轨的位置精 度优于10 m,速度精度为2 cm／s,预报5分钟轨道精度为15.02 m,满足了短弧跟踪条件下R DSS对轨道精度的要求,实现了卫星精密定轨。
     

月球卫星轨道摄动及冻结轨道研究

利用理论分析、数值仿真与相图分析,论述了月球卫星冻结轨道与地球卫星冻结轨道的区别,分析结果表明,月球重力场存在较大异常,会引起月球卫星轨道发生较大漂移。月球冻结轨道在田谐项影响下,还存在中等周期的漂移。仅简单考虑带谐项系数,无法求得完美的月球冻结系数。月球重力场异常对绕月卫星的影响与地球相比存在很大区别。月球轨道卫星的长期运行与控制策略的设计,不能按照地球轨道卫星的传统方法。目前使用的月球引力模型精度较差,尽管基于这些不可靠的引力模型,可以得出很多有用结论,但对未来高精度的月球探测任务来说,还存在不足,需要在将来的月球探测任务中,探测高精度的月球重力场,以利于未来月球探测航天系统的任务分析与设计。     

小推力地球卫星圆轨道同轨调相设计方法研究

在研究小推力地球卫星圆轨道同轨调相任务设计问题的基础上,提出了一种半解析调相参数分析方法,并据此发展了一种高效的精确设计方法.首先,根据推力方向对卫星相位角的影响规律,采用推力方向假设和轨道平均技术推导了调相轨道参数满足的函数关系,通过微分修正可快速获得调相轨道的关键参数;然后,基于摄动轨道动力学模型,通过对控制量进行离散化建立了复杂约束条件下的燃料最省调相轨道设计模型,并以初始分析结果为初值,采用序列二次规划算法进行求解;最后,以地球同步轨道调相任务为例对所提方法进行了数值验证.数值仿真结果表明:提出的初始分析方法可为调相轨道设计提供合理的初值猜测,发展的鲁棒设计方法可有效用于摄动模型复杂约束下的小推力调相轨道设计.     

基于J2项和大气阻力摄动的卫星编队保持方法研究

研究J2摄动和大气阻力对低轨编队卫星相对位置的影响,在此基础上给出一种编队保 持方案。文中定量分析了J2摄动对编队卫星三轴相对位置的影响,给出了大气阻力对编队 卫 星相对轨道要素影响表达式。在同时考虑两种摄动力前提下,推导给出了x方向受摄动 的漂移量Δx与编队卫星轨道长半轴之差Δa的周期变化量之间的解析关系式,基于 该关系式,设计了单边极限环形式的卫星长期编队保持控制方案。最后通过数学仿真验证了 该方案的可靠性,仿真结果与理论分析相符。该控制方案在计算控制量时只需知道编队卫星 的轨道长半轴之差,容易实现,为工程实践提供依据。
     

CE5T拓展试验轨道精度分析

针对嫦娥5T服务舱(CE5T)拓展试验中的绕地大椭圆轨道,分析了轨道动力学演化趋势,通过测轨数据类型组合策略分析了统一S频段测量(USB)和甚长基线干涉测量(VLBI)在定轨中的贡献,得到了百米级的精密定轨精度;针对地月第二平动点(L2)绕飞轨道,分析了地心和月心积分的轨道动力学差异,制定了精密定轨的参数求解策略,得到了百米级的精密定轨精度;针对月球交会对接轨道的特点,选取三种不同的重力场模型定轨,比较了三者在轨道预报和数据拟合的差异,并与嫦娥3号(CE3)环月轨道的定轨精度进行比对,验证了不同重力场的适用范围,从计算精度和效率两方面制定了优化的定轨策略。     

月球卫星的自主轨道确定

分别利用紫外三轴姿态敏感器测出的月心方向和测距仪测出的月心距作为观测量,采用推广卡尔曼滤波方法实时确定出月球卫星环月期间的轨道,考虑到月球卫星在环月期间太阳对月球光照条件的限制,提出分段虑波策略,并给出了仿真分析.     

系统参数对电动力绳系动力学的影响

将电动力绳系(EDT)的主星质量、子星质量、绳系质量以及绳系中的电流视为系统参数,研究这些参数对系统的摆动动力学和轨道动力学的影响。哑铃模型下的电动力绳系摆动动力学方程存在不稳定的周期解,通过Floquet理论来衡量周期解的不稳定程度,从而研究各系统参数对摆动动力学的影响。建立了用春分点轨道元素的形式描述的电动力绳系轨道动力学方程,并以降轨时间来衡量电动力绳系的降轨效率,从而研究系统参数对轨道动力学的影响。运用算例对周期解迁移矩阵的特征值、降轨时间随各系统参数的变化关系进行了仿真,分别得出了各系统参数对系统摆动动力学和轨道动力学的影响。综合本文的仿真结果,并考虑实际发射及空间运行中的其它因素,对电动力绳系的设计和降轨策略提出了建议。     

解析计算在月球中继卫星Halo轨道设计中的应用

面向月球中继卫星工程轨道设计需求,研究解析计算方法在地月系L2点halo轨道设计中的应用问题。在讨论圆型限制性三体问题三阶解析近似计算方法的基础上,分析了解析计算与数值计算的差异,给出了解析近似计算在工程约束下的适用范围,进而提出了基于解析计算的轨道设计和特征筛选方法。分别采用解析初值和数值初值进行halo轨道外推,比对验证采用解析计算设计轨道的可行性。研究结果表明,解析计算方法适用于月球中继卫星轨道的初步设计、特征分析和构型筛选。     

大行星、月球和小天体环绕型探测器的轨道问题

针对环绕型探测器的轨道,以具有不同物理特征的大天体(金星、火星和月球等)为例,根据它们的质量大小(决定其引力作用范围的大小)、密度分布和形状特征(决定其非球形引力场的特征),具体阐述了它们各自环绕型探测器轨道的可能形式和不同的变化特征,另外介绍了小天体探测时的两种伴飞形式。关于环绕型探测器的轨道特征,如月球无大气,对其低轨探测器而言没有能量耗散影响,却同样有轨道寿命问题,而且由于其质量分布不均匀,这种现象还与轨道倾角有密切关系;而对火星的环绕型探测器而言,尽管有类似现象,但与轨道倾角的关系却大不相同。这些重要的轨道变化特征,可为深空探测中目标轨道的选择和设计提供依据。     

Development of a high-precision selenodetic coordinate system for the physical surface of the Moon based on LED beacons on its surface

The paper presents a mathematical algorithm for processing an array of angular measurements of light beacons on images of the lunar surface onboard a polar artificial lunar satellite (PALS) during the Luna–Glob mission and coordinate–time referencing of the PALS for the development of reference selenocentric coordinate systems. The algorithm makes it possible to obtain angular positions of point light beacons located on the surface of the Moon in selenocentric celestial coordinates. The operation of measurement systems that determine the position and orientation of the PALS during its active existence have been numerically simulated. Recommendations have been made for the optimal use of different types of measurements, including ground radio trajectory measurements, navigational star sensors based on the onboard star catalog, gyroscopic orientation systems, and space videos of the lunar surface.     

嫦娥五号平动点拓展任务轨道方案研究

中国嫦娥五号探测器成功实现月球采样返回任务,为最大限度利用任务资源,研究了利用嫦娥五号轨道器的平动点拓展任务轨道方案,设计了平动点轨道及其转移轨道.首先,给出了任务轨道设计的轨道动力学模型,包括圆型限制性三体问题模型和精确力模型.其次,基于嫦娥二号和嫦娥5T1平动点拓展任务设计经验,介绍了平动点轨道直接转移与入轨等轨道...     

Theory of physical libration of the Moon with the liquid core: Forced librations

For a two-layer model of the Moon that consists of a solid nonspherical mantle and an ellipsoidal homogeneous liquid core, a theory of forced librations under the effect of gravitational Earth’s moments has been developed. The motion of the Moon over its orbit has been described by the high-accuracy theory of DE/LE-4 orbital motion. Tables have been constructed that present forced librations of the Moon caused by the second harmonic of its force function, in the neighborhood of its motion according to the generalized Cassini laws. Disturbances of the first-order with respect to dynamic compressions of the Moon and its core are obtained in analytical form for Andoyer variables and Poincare variables and for the projection of the angular velocity vector of Moon’s mantle rotation and the Poincare coordinate system (relative to which core’s liquid accomplishes simple motion) on its major central axes of inertia, as well as for the classical variables in the Moon libration theory, etc. Constructed tables of the forced librations theory give the amplitudes and periods of librations and combinations of arguments of the orbital motion theory that correspond to libration parameters. The interpretation of basic variations has been given and a comparison with the previous theories has been carried out, in particular with the modern empirical theory constructed based on the laser observation data.     

Optimal low-thrust transfers between close near-circular coplanar orbits

Four types of optimal solutions are demonstrated to exist for transfers (time of flight is not fixed) between close near-circular coplanar orbits. One solution is realized with the help of fixed orientation of the propulsion system (PS) along a transversal in the orbital coordinate system. Another is reached at fixed orientation of the PS in the inertial coordinate system. The third and fourth types of solutions change the PS orientation in the process of executing the maneuver. Regions of existence are established for all types of solutions, and algorithms for determination of parameters of these maneuvers are suggested. The algorithms were used to calculate parameters of the maneuvers of transfer from a launching orbit to a working Sun-synchronous orbit, and to calculate the maneuvers of supporting the parameters of such an orbit in a specified range.     

Manifold dynamics in the Earth–Moon system via isomorphic mapping with application to spacecraft end-of-life strategies

Recently, manifold dynamics has assumed an increasing relevance for analysis and design of low-energy missions, both in the Earth–Moon system and in alternative multibody environments. With regard to lunar missions, exterior and interior transfers, based on the transit through the regions where the collinear libration points L₁ and L₂ are located, have been studied for a long time and some space missions have already taken advantage of the results of these studies. This paper is focused on the definition and use of a special isomorphic mapping for low-energy mission analysis. A convenient set of cylindrical coordinates is employed to describe the spacecraft dynamics (i.e. position and velocity), in the context of the circular restricted three-body problem, used to model the spacecraft motion in the Earth–Moon system. This isomorphic mapping of trajectories allows the identification and intuitive representation of periodic orbits and of the related invariant manifolds, which correspond to tubes that emanate from the curve associated with the periodic orbit. Heteroclinic connections, i.e. the trajectories that belong to both the stable and the unstable manifolds of two distinct periodic orbits, can be easily detected by means of this representation. This paper illustrates the use of isomorphic mapping for finding (a) periodic orbits, (b) heteroclinic connections between trajectories emanating from two Lyapunov orbits, the first at L₁, and the second at L₂, and (c) heteroclinic connections between trajectories emanating from the Lyapunov orbit at L₁ and from a particular unstable lunar orbit. Heteroclinic trajectories are asymptotic trajectories that travels at zero-propellant cost. In practical situations, a modest delta-v budget is required to perform transfers along the manifolds. This circumstance implies the possibility of performing complex missions, by combining different types of trajectory arcs belonging to the manifolds. This work studies also the possible application of manifold dynamics to defining suitable, convenient end-of-life strategies for spacecraft orbiting the Earth. Seven distinct options are identified, and lead to placing the spacecraft into the final disposal orbit, which is either (a) a lunar capture orbit, (b) a lunar impact trajectory, (c) a stable lunar periodic orbit, or (d) an outer orbit, never approaching the Earth or the Moon. Two remarkable properties that relate the velocity variations with the spacecraft energy are employed for the purpose of identifying the optimal locations, magnitudes, and directions of the velocity impulses needed to perform the seven transfer trajectories. The overall performance of each end-of-life strategy is evaluated in terms of time of flight and propellant budget.