A Study of the Reflection of Radio Waves from the Plasma Region on the Track behind an Axisymmetric Body Moving in the Earth’s Atmosphere

The results of numerical calculation of the dependences of the electron density, the eigenfrequency and the dielectric plasma permeability on the geometric parameters and the altitude of body motion in the near and far wake behind a thin conical body with a spherical nose blunting have been presented. The electron density maximum has been shown to be located in the region of the neck of the near wake behind the body, which determines the type of this region (supercritical or subcritical). This in turn affects the propagation of radio waves through this plasma region. A comparative analysis was performed for two different bodies with the same ballistic coefficient values. No characteristic distinctions were revealed in the values of electron density or the plasma eigenfrequency in the near and far wake behind these bodies. However, it has been shown that there are differences in the values of the distance from the bottom cross section to the neck of the near wake behind these bodies.     

Design of Interplanetary Transfers with Passive Gravity Assists and Deep Space Maneuvers

In the paper, the problem of designing interplanetary trajectories with several swing-bys and deep-space maneuvers is solved using the method of virtual trajectories developed by the authors. The algorithms for the calculation of both heliocentric and planetocentric trajectory arcs are presented, including the case of resonant trajectories. The results of applying the method of virtual trajectories to the problem of designing an interplanetary transfer to Jupiter are given and compared with the baseline trajectories for the Juno, Europa Clipper, and Laplace missions.     

Detection of Molecular Scattering Field from a Polarization Analysis of the Sky Background during Transitive Twilight and Temperature Measurements near the Stratopause

The simplest version of the method of detecting the single molecular scattering field based on the polarization measurements of the twilight sky background by all-sky cameras has been considered. The method can be used during transitive twilight (with solar zenith angles of 94°–98°), when effective single scattering occurs in the upper stratosphere and lower mesosphere. The long-term measurements conducted using this method in the Moscow region and Apatity make it possible to determine the temperature of these atmospheric layers and estimate the contribution and properties of multiple scattering during the transitive twilight.     

Gravity-oriented satellite dynamics subject to gravitational and active damping torques

The dynamics of the rotational motion of a satellite moving in the central Newtonian field of force over a circular orbit under the effect of gravitational and active damping torques, which depend on the satellite angular velocity projections, has been investigated. The paper proposes a method of determining all equilibrium positions (equilibrium orientations) of a satellite in the orbital coordinate system for specified values of damping coefficients and principal central moments of inertia. The conditions of their existence have been obtained. For a zero equilibrium position where the axes of the satellite-centered coordinate system coincide with the axes of the orbital coordinate system, the necessary and sufficient conditions for asymptotic stability are obtained using the Routh–Hurwitz criterion. A detailed analysis of the regions where the conditions of the asymptotic stability of a zero equilibrium position are fulfilled have been obtained depending on three dimensionless parameters of the problem, and the numerical study of the process of attenuation of satellite’s spatial oscillations for various damping coefficients has been carried out. It has been shown that there is a wide range of damping parameters from which, by choosing the necessary values, one can provide the asymptotic stability of satellite’s zero equilibrium position in the orbital coordinate system.     

天基强银河背景下的空间碎片探测试验及分析

本文设计并完成了银河背景下碎片探测试验,数据分析表明,强的银河背景会降低探测器对暗弱碎片的探测能力,但对碎片的天文定位精度影响不大。本文验证了采用基于探测概率的极限探测星等方法表征探测器极限探测能力的合理性,以及基于仿真的方法增加数据完备度的可行性。本文的工作将为后续天基探测器的设计及探测试验的开展提供依据。     

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

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

基于深度学习的组合体航天器模型预测控制

利用模型预测算法先预测控制结果后控制的类人行为特点,借助深度学习在多参数寻优上的优势,提出了一种基于卷积神经网络的模型预测控制算法,满足航天工程低硬件需求,实现组合航天器多场景下姿态控制律的重构。该算法首先利用模型预测控制将组合航天器从初始状态控制到预期状态,然后将控制过程中状态量用于3层3核卷积神经网络的训练,训练完成后,用该卷积神经网络代替模型预测对组合航天器进行控制,从而降低计算资源需求。仿真校验表明：该算法可预测5个控制周期内的控制参数,相比传统模型预测算法所需硬件计算时间降低约5倍,在一般硬件环境下30 s内即可完成各场景下的组合航天器姿态控制,控制精度在10 -4 量级。     

椭圆轨道卫星对地全球覆盖电推进控制

针对椭圆轨道卫星近/远地点的星下点对全球或特定纬度区域的访问问题,提出一种连续小推力下的对地覆盖控制策略。首先,推导了自然摄动对卫星拱线变化的影响,并探讨了进行小推力覆盖控制的必要性。然后,针对燃料消耗的优化问题,将控制方程展开成含傅里叶级数的形式,用以获得便于星上计算的解析形式的次优解,同时探讨了截取阶数与优化程度的关系。在进行拱线控制的同时,通过合理设置约束,对椭圆轨道的近地点高度进行保护,确保卫星安全运行。仿真结果表明,提出的方法能够以适当的燃料消耗代价实现椭圆轨道的近/远地点的全球覆盖控制或特定纬度区域的反复推扫,且控制力在可接受的范围内。     

分数阶控制器实现方法研究

研究了分数阶控制器的模拟和数字实现方法,其中数字实现是重点。数字实现是通过构造一个新的IIR型数字滤波器来对函数的分数阶微积分计算进行近似求解,然后对该滤波器进行连分式展开来完成的。以PID控制为例进行的数值仿真结果表明,用此方法设计的离散分数阶PID控制器比传统离散整数阶PID控制器具有更优的控制性能,同时也验证了该数字实现方法的正确性与可行性。     

旋转弹气动力建模与飞行轨迹仿真

运用一种自回归滑动平均(ARMA)的时域气动力建模方法,以计算流体力学与刚体动力学(CFD/RBD)耦合仿真的输出结果为样本,对旋转弹的非定常气动力进行建模。利用建立的气动力模型与刚体动力学方程求解模块耦合,实现了旋转弹轨迹的快速仿真,并讨论了不同的建模方式对仿真精度的影响。算例结果表明:采用气动力模型与刚体动力学方程耦合仿真技术可以在不同初始发射条件下进行旋转弹飞行姿态与运动轨迹预测,且与CFD/RBD仿真结果吻合较好,证明ARMA气动力建模方法可以在保证旋转弹轨迹预测精度的同时大幅缩短仿真时间,节省计算资源。