Optimization of low-energy resonant hopping transfers between planetary moons

In response to the scientific interest in Jupiter's Galilean moons, NASA and ESA have plans to send orbiting missions to Europa and Ganymede, respectively. The inter-moon transfers of the Jovian system offer obvious advantages in terms of scientific return, but are also challenging to design and optimize due in part to the large, often chaotic, sensitivities associated with repeated close encounters of the planetary moons. The approach outlined in this paper confronts this shortcoming by exploiting the multi-body dynamics with a patched three-body model to enable multiple “resonant-hopping” gravity assists. Initial conditions of unstable resonant orbits are pre-computed and provide starting points for the elusive initial guess associated with the highly nonlinear optimization problem. The core of the optimization algorithm relies on a fast and robust multiple-shooting technique to provide better controllability and reduce the sensitivities associated with the close approach trajectories. The complexity of the optimization problem is also reduced with the help of the Tisserand–Poincaré (T–P) graph that provides a simple way to target trajectories in the patched three-body problem. Preliminary numerical results of inter-moon transfers in the Jovian system are presented. For example, using only 59 m/s and 158 days, a spacecraft can transfer between a close resonant orbit of Ganymede and a close resonant orbit of Europa.     

Control of satellite imaging formations in multi-body regimes

Libration point orbits may be ideal locations for satellite imaging formations. Therefore, control of these arrays in multi-body regimes is critical. A continuous feedback control algorithm is developed that maintains a formation of satellites in motion that is bounded relative to a halo orbit. This algorithm is derived based on the dynamic characteristics of the phase space near periodic orbits in the circular restricted three-body problem (CR3BP). By adjusting parameters of the control algorithm appropriately, satellites in the formation follow trajectories that are particularly advantageous to imaging arrays. Image reconstruction and coverage of the (u, v) plane are simulated for interferometric satellite configurations, demonstrating potential applications of the algorithm and the resulting motion.     

Rapid generation of entry trajectories with waypoint and no-fly zone constraints

Waypoints are positions for multiple payload deployments or reconnaissance missions, and no-fly zones are exclusion zones that cannot be passed for threat avoidance or due to geopolitical restrictions. This paper proposes a rapid entry trajectory generation approach satisfying waypoint and no-fly zone constraints for entry vehicles with relatively high lift-to-drag ratio. A lateral planning algorithm based on a Newton iteration scheme is developed to simultaneously design both the magnitude and sign of the control variable according to waypoints and no-fly zones. The algorithm converts the highly constrained trajectory planning problem into a series of one-parameter search problems based on a reduced-order system. Then, the quasi-equilibrium glide phenomenon is employed to extract the remaining state variables corresponding to longitudinal motion. The algorithm is tested using the Common Aero Vehicle model, and the results demonstrate that the algorithm can generate flyable entry trajectories rapidly within allowable tolerances while satisfying all the flight constraints.     

基于Gauss伪谱法的高超声速飞行器多约束三维再入轨迹优化

高超声速飞行器再入飞行过程中,需要满足多种过程约束和终端状态约束,同时再入初始状态根据飞行任务不同会有较大变化,针对其特点的快速轨迹优化问题已成为当今热点。本文研究了一种基于"初值轨迹生成+Gauss伪谱法+SQP求解NLP"的方法,既利用了Gauss伪谱法收敛快、精度高的特点,又结合初值轨迹生成算法,弥补了Gauss伪谱法对初值敏感的不足。本文在仿真过程中选取再入总吸热量最小为性能指标,求解了满足多种约束的再入轨迹,并将优化的结果与数值积分的结果进行比较,验证了此算法有效性和可行性。     

Optimal trajectories for the aeroassisted flight experiment

This paper deals with the determination of optimal trajectories for the aeroassisted flight experiment (AFE). The intent of this experiment is to simulate a GEO-to-LEO transfer, where GEO denotes a geosynchronous Earth orbit and LEO denotes a low Earth orbit. Specifically, the AFE spacecraft is released from the Space Shuttle and is accelerated by means of a solid rocket motor toward Earth, so as to achieve atmospheric entry conditions identical with those of a spacecraft returning from GEO. During the atmospheric pass, the angle of attack is kept constant, and the angle of bank is controlled in such a way that the following conditions are satisfied: (a) the atmospheric velocity depletion is such that, after exiting, the AFE spacecraft first ascends to a specified apogee and then descends to a specified perigee; and (b) the exit orbital plane is identical with the entry orbital plane. The final maneuver, not analyzed here, includes the rendezvous with and the capture by the Space Shuttle. In this paper, the trajectories of an AFE spacecraft are analyzed in a 3D space, employing the full system of 6 ODEs describing the atmospheric pass. The atmospheric entry conditions are given, and the atmospheric exit conditions are adjusted in such a way that requirements (a) and (b) are met, while simultaneously minimizing the total characteristic velocity, hence the propellant consumption required for orbital transfer. Two possible transfers are considered: indirect ascent (IA) to a 178 NM perigee via a 197 NM apogee; and direct ascent (DA) to a 178 NM apogee. For both transfers, two cases are investigated: (i) the bank angle is continuously variable; and (ii) the trajectory is divided into segments along which the bank angle is constant. For case (ii), the following subcases are studied; 2, 3, 4 and 5 segments; because the time duration of each segment is optimized, the above subcases involve 4, 6, 8 and 10 parameters, respectively. It is shown that the optimal trajectories of cases (i) and (ii) coalesce into a single trajectory: a two-subarc trajectory, with the bank angle constant in each subarc (bang-bang control). Specifically, the bank angle is near 180° in the atmospheric entry phase (positive lift projection phase) and is near 0° in the atmospheric exit phase (negative lift projection phase). It is also shown that, during the atmospheric pass, the peak values of the changes of the orbital inclination and the longitude of the ascending node are nearly zero; hence, the peak value of the wedge angle (angle between the instantaneous orbital plane and the initial orbital plane) is nearly zero. This means that the motion of the spacecraft is nearly planar in an inertial space.     

舱体入水工况参数对冲击特性的影响分析

为分析大型舱体返回入水时入水垂直速度等工况参数对其冲击特性的影响,提出一种基于任意拉格朗日欧拉法的舱体入水过程模拟方法,应用球底结构入水问题理论计算方法验证该模拟方法的有效性。通过对不同入水工况的有限元模拟,分析入水垂直速度和入水角度对冲击特性的影响。结果表明本文建立的舱体-流体有限元模型能够有效模拟舱体入水过程,舱体冲击加速度峰值与入水垂直速度成正比,舱体以一定角度入水能够降低入水过程的冲击加速度。研究结果可以为新型舱体的结构设计和入水冲击试验提供指导,从而减少试验次数,缩短开发周期。     

多体卫星地面物理仿真的一致性研究

三轴气浮台不能直接用于多体卫星的地面物理仿真实验,这是由于部件转动引起台体 的质心变化,进而产生重力静不平衡力矩,使仿真过程无法进行,因此需要质心补偿系统对 气浮台的重力静不平衡力矩进行补偿。本文通过对多体卫星动力学和气浮台动力学的比较, 分析了在重力场中的地面物理仿真系统与真实卫星的动力学特性的不同点,给出地面物理仿 真系统能完全模拟真实卫星运动的一致性条件,并进行了数学仿真。当转动附件的质心在转 动轴或转动中心上时,气浮台系统可以完全模拟卫星系统的质量特性和力学耦合特性,否则 ,多体气浮台系统和多体卫星系统的动力学耦合特性不完全一致,需要采取一定的方法进行 补偿。
     

火星探测器进入阶段稳定性分析

文章针对＂火星探路者＂,对火星探测器在进入火星大气减速过程中的运动稳定性问题进行了研究,分析了质心位置、进入条件、大气参数等因素对火星探测器运动稳定性的影响。分析结果表明,在马赫数大于1.2时,探测器处于一个比较稳定的状态。     

大挠性多体卫星的自抗扰姿态控制系统设计

研究了最近提出的非线性自抗扰控制器 (ADRC)在大挠性多体卫星姿态控制中的应用问题。提出了一种双闭环自抗扰姿态控制器 (ADRAC) ,并从鲁棒性、干扰抑制、动静态性能指标和振动抑制等方面与使用于某挠性卫星的传统PID控制器进行了比较。在考虑了反作用轮和敏感器的实际非线性和敏感器噪声的影响情况下进行了仿真 ,结果表明 ,双闭环自抗扰姿态控制器各方面性能均优于传统的PID。本文提出的自抗扰姿态控制器 ,对实现挠性多体卫星的高精度高稳定度姿态控制 ,具有实际的应用价值     

Venus transfer design by combining invariant manifolds and low-thrust arcs

The design of interplanetary trajectories based on patched circular restricted three body models is gradually becoming a valuable alternative to the classical patched conic approach. The main advantage offered by such a model is the possibility to exploit the manifold dynamics to move naturally far from or toward a body. Generally, propulsive maneuvers are required to match these structures. Low-thrust arcs offer the possibility to have a significant propellant mass reduction when moving from manifold to manifold. The aim of this paper is to present a methodology to design low-thrust trajectories between two planetary orbits connecting the manifolds of two circular three body systems. The approach is based on a grid search on the main parameters governing the solution to identify those trajectories moving within the manifold images on given Poincarè sections. The value of the Jacoby constant of the target libration point periodic orbit is chosen as stop condition for the thrusting phases. Ballistic arcs follow up to the proper Poincarè section intersection. A grid search for an Earth to Venus transfer is presented as test case.     

复杂结构充液航天器晃动动力学与姿态稳定性

复杂充液航天器的晃动及其对控制系统的影响是当前空间高技术的重要问题。 本文首先由力学变分原理导出充液复杂系统的方程式,其中包括流体力学方程和相应的边界条件,以及弹性连续介质力学方程;特别是,考虑了强毛细作用力的影响。其次,考查了微重条件下带有隔板的贮箱内的液体晃动、带有网孔隔板腔体内的液体晃动和粘弹性隔板对抑制晃动的作用。另外,本文将分析带有多腔充液自旋卫星的姿态稳定性,以及晃动、液体涡旋、贮箱偏置,能量耗散和哥氏加速度对卫星运动的影响。最后,讨论了多体系统充液挠性空间飞行器的稳定性问题。 同时,本文给出了数值计算结果和某些试验结果,并与理论分析作了对比。     

Trajectories of in-plane periodic solutions of tethered satellite system projected on van der Pol planes

In this paper, in-plane periodic solutions for a dumbbell model in elliptic orbits are searched using bifurcation, and their trajectories are projected on the van der Pol planes. The trajectories projected on the van der Pol plane will be employed as a tool to predict when the control of delayed feedback control will need to act to maintain the periodic motions.     

基于地球引力变轨的日冕探测器轨道设计

研究了相对黄道面有一定倾角的探测器轨道设计的问题。以金星借力轨道设计为例,分析了轨道偏心率与轨道倾角增量之间的关系。根据C3匹配原理搜索了“地球-中间天体-地球”多天体交会的发射窗口。最后,设计了与地球轨道周期相等的三次地球借力轨道,该轨道倾角可以达到黄纬30°以上。理论分析及仿真结果表明:基于地球引力设计此类轨道时,应采用多天体交会方案,才能既保证地球逃逸能量低,又保证首次飞入地球影响球前轨道偏心率较大的双重指标;同时应采用多次地球借力方案,该方案具有每次借力后轨道偏心率逐渐减小的特点,当其减小到零时,再次借力后轨道倾角不会继续增加。     

The Motion of a Material Point in a Newtonian Field under the Action of Equiangular Thrust

The problems of investigation and optimization of the motion of spacecraft are extensively discussed in the literature. Nevertheless, in many cases a large variety of qualitative characteristics of their motion and of the form of their trajectories are still unclear. In this paper we consider a plane equiangular acceleration of a spacecraft both in a Newtonian field and in its absence (at a large distance from the center of attraction). The general equation of a trajectory of plane acceleration is presented with the introduction of a new variable, an index of an exponent, which allows one to obtain convenient solutions at different values of the time-independent angle of inclination of the vector of thrust to the spacecraft's radius vector (i.e., when equiangular acceleration takes place). Asymptotic solutions are constructed and an interesting fact is revealed. Namely, it is shown that when the center of attraction exists or is absent, for all initial conditions the trajectories appearing at the above equiangular acceleration of a material point tend to the standard logarithmic spirals at a large distance from the center. Specifically, when the value of transverse (perpendicular to the radius vector) thrust is constant, there appears a logarithmic spiral with an angle of inclination to the radius vector equal to 35.264°. Different forms of the trajectory of equiangular acceleration of spacecraft at a low thrust are also studied. The results obtained can be useful for the investigation and choice of optimum space trajectories.     

垂直回收火箭可伸缩着陆腿展开的动力学仿真分析

可伸缩着陆腿是实现运载火箭垂直回收的有效途径。相比于着陆过程，火箭着陆腿展开过程涉及的动力学问题更值得特别关注。根据多体动力学方法建立可伸缩着陆腿的虚拟样机模型，针对展开过程开展仿真分析，研究了支腿结构柔性、腿节滑移摩擦力、气动罩气动力及过载等因素对展开时长和展开到位角速度两个关键展开特性指标的影响。结果表明，结构柔性对展开特性影响可忽略，过载对展开特性的影响最大，其次是气动罩气动力和腿节摩擦力。相关研究结果可为用于火箭垂直回收的可伸缩式着陆腿设计提供参考。     

携带晃动燃料柔性航天器姿态机动中的同宿环分叉研究

本文研究了带液体晃动和柔性附件的耦合航天器系统在液体燃料耗散和柔性附件扭转振动的作用下,经历从最小惯量轴到最大惯量轴姿态机动中的混沌动力学行为.将液体晃动等效为球摆模型并由此建立了带柔性附件充液航天器多体耦合系统动力学模型.首先推导出耦合系统动力学方程并采用Melnikov积分预测受扰系统稳定与不稳定流形是否横截相交,得到了参数形式表达的混沌运动解析判据,这对航天器的设计有重要的指导意义.研究发现,混沌的发生依赖于刚体形状,阻尼比,充液比和扭转振动频率.此外,在经过被动再定向姿态机动后,由于液体晃动的本质非线性特性,充液航天器最终将进行大章动角的周期极限环运动而非绕着最大惯量轴自旋.     

商用微电子技术在卫星辐射环境中的应用

汇总了商用器件的辐射效应,讨论了这些器件的辐射敏感度,讨论了对其危害的评估方法和必要时的回避技术,以帮助设计师和项目管理人员在空间辐射环境中用好这些商用器件。     

Can solar sails be used to test fundamental physics?

The relative importance of certain general relativistic effects is enhanced by solar radiation pressure (SRP). The observation and study of the trajectories of a solar sail could potentially provide tests of various effects of general relativity. In particular, we study Keplerian and non-Keplerian orbits near the sun as well as escape trajectories for a solar sail, for which general relativistic effects and the solar radiation pressure are considered simultaneously. In contrast with the conventional solar mission, a solar sail allows for non-Keplerian orbits, for which the orbital plane lies above the sun. It is predicted that there is an analog of the Lense–Thirring effect for non-Keplerian orbits. Also the SRP increases the amount of precession per orbit due to the Lense–Thirring effect for polar heliocentric orbits. A solar sail would also enhance the relative importance of effects associated with a possible net charge on the sun and during many rotations this effect may be measurable.     

Trajectories for missions of low-thrust spacecraft aimed at delivery of soil samples from main belt asteroids and Phobos

Trajectories of spacecraft with electro-jet low-thrust engines are studied for missions planning to deliver samples of matter from small bodies of the Solar System: asteroids Vesta and Fortuna, and Martian moon Phobos. Flight trajectories are analyzed for the mission to Phobos, the limits of optimization of payload spacecraft mass delivered to it are determined, and an estimate is given to losses in the payload mass when a low-thrust engine with constant outflow velocity is used. The model of an engine with ideally regulated low thrust is demonstrated to be convenient for calculations and analysis of flight trajectories of a low-thrust spacecraft.     

On the evolution of weak discontinuities in radiative magnetogasdynamics

A method of wavefront analysis is used to analyze the formation of shock waves in a two-dimensional steady supersonic flow past plane and axisymmetric bodies in radiative magnetogasdynamics. The gas is assumed to be perfectly conducting and to be permeated by a magnetic field orthogonal to the trajectories of the gas particles. The medium is taken to be sufficiently hot for the effects of thermal radiation to be significant, which is treated by the optically thin approximation to the radiative transfer equation. Transport equations, which lead to the determination of the distance at which the first characteristic could intersect with a successive one and also to conditions, which insure that no shock will ever evolve on the wave front, are derived. The effect of upstream flow Mach number and the magnetic field strength on the behavior of the wavefront are examined.