航天器最优受控绕飞轨迹推力幅值延拓设计方法

针对航天器燃料最优可变周期绕飞轨迹的求解问题, 提出了一种以推力幅值为延拓参数的延拓方法. 问题求解从最为简单的双脉冲绕飞轨迹出发, 首先利用有限推力替代脉冲推力, 设定推力序列为“开— 关— 开”,然后逐步减小推力幅值, 最终得到最小推力绕飞轨迹; 此后, 再逐步增加推力幅值, 结合主矢量曲线判断最优推力开关序列, 将最小推力解延拓至有限推力以及脉冲推力燃料最优解. 该方法通过对推力幅值的延拓, 实现了有限推力bang-bang 控制与脉冲推力燃料最优绕飞轨迹优化问题的一并求解, 同时避免了最优控制问题中协态变量的随机猜测. 慢速与快速绕飞算例的优化结果验证了所提出方法的有效性.      

基于虚拟中心引力场方法的航天器转移轨道设计

空间机动技术是实现空间操作任务的基础,具有重要的研究价值. 研究了连续推力作用下航天器转移轨道设计问题,提出了一种基于虚拟中心引力场的轨道设计方法. 该方法有两大特点:(1) 能够将机动轨道设计问题转化为虚拟中心引力场参数的优化问题,简化了设计过程;(2) 对轨道形状或推力方向、大小不做任何假定,能够应用于一般情况下的机动轨道设计. 将该方法应用于航天器二维和三维的转移轨道设计,并和形状方法进行了对比分析. 仿真结果分析表明,采用该方法简化了轨道设计过程,为航天器快速轨道设计提供了新思路.      

基于星历匹配法的载人小行星探测轨迹优化问题求解

对基于星历匹配法的载人小行星探测轨迹优化问题解法进行了研究. 提出星历匹配法的概念,利用其求出最优发射窗口初值与候选最优探测序列,然后利用遗传算法对特定探测序列的时间节点等参数进行优化,最后利用基于庞德里亚金极大值原理的同伦法进行小推力转化以求得最终探测轨道及其推力控制律. 结果显示星历匹配法可以快速准确地求出最优发射窗口初值与候选探测序列,这大大提高了载人小行星探测等多目标交会探测问题的求解效率.     

On optimization of many-revolution low-thrust orbit transfers: Part 1

The paper proposes a numerical approach to the problem of optimal control of low-thrust spacecraft in a strong central gravity field. The approach employs the solution of the averaged equations of optimal motion. For an optimal variable-thrust maneuver, it is shown that quasioptimal trajectories are close to the averaged ones. It is established that in the case of constant thrust, the averaged solutions can be used as a satisfactory first approximation for finding a quasioptimal solution by minimizing the discrepancy at the right end of the trajectory __________ Translated from Prikladnaya Mekhanika, Vol. 44, No. 7, pp. 110–119, July 2008.     

Optimal design of near-Earth asteroid sample-return trajectories in the Sun–Earth–Moon system

In the 6th edition of the Chinese Space Trajectory Design Competition held in 2014, a near-Earth asteroid sample-return trajectory design problem was released, in which the motion of the spacecraft is modeled in multi-body dynamics, considering the gravitational forces of the Sun,Earth, and Moon. It is proposed that an electric-propulsion spacecraft initially parking in a circular 200-km-altitude low Earth orbit is expected to rendezvous with an asteroid and carry as much sample as possible back to the Earth in a10-year time frame. The team from the Technology and Engineering Center for Space Utilization, Chinese Academy of Sciences has reported a solution with an asteroid sample mass of 328 tons, which is ranked first in the competition.In this article, we will present our design and optimization methods, primarily including overall analysis, target selection, escape from and capture by the Earth–Moon system,and optimization of impulsive and low-thrust trajectories that are modeled in multi-body dynamics. The orbital resonance concept and lunar gravity assists are considered key techniques employed for trajectory design. The reported solution, preliminarily revealing the feasibility of returning a hundreds-of-tons asteroid or asteroid sample, envisions future space missions relating to near-Earth asteroid exploration.     

Low-thrust trajectory optimization for multiple target bodies tour mission

Spacecraft science missions to planets or asteroids have historically visited only one or several celestial bodies per mission. The research goal of this paper is to create a trajectory design algorithm that generates trajectory allowing a spacecraft to visit a significant number of asteroids during a single mission. For the problem of global trajectory optimization, even with recent advances in low-thrust trajectory optimization, a full enumeration of this problem is not possible. This work presents an algorithm to traverse the searching space in a practical fashion and generate solutions. The flight sequence is determined in ballistic scenario, and a differential evolution method is used with constructing a three-impulse transfer problem, then the local optimization is implemented with low-thrust propulsion on the basis of the solutions of impulsive trajectories. The proposed method enables trajectory design for multiple asteroids tour by using available low thrust propulsion technology within fuel and time duration constraints.     

THE EVOLUTION CHARACTERIZATION OF LIBRATIN POINT ORBITS AND APPLICATION RESEARCH

三体问题中, 轨道的受力和运动规律非常复杂. 对于特定的任务, 如何选择轨道的初始解是一大难题.针对平面三体问题, 利 用近拱点庞加莱映射, 对平动点顺行轨道和逆行轨道的长期和短期演化性质进行分析.根据轨道的初始状态将其分为逃逸轨道和捕获轨道.对于逃逸轨道, 给出了同宿轨道和异宿轨道的设计方法, 并利用两级微分修正法消除了拼接点处的位置不连续问题.对于捕获轨道, 得到了几类典型的周期和准周期轨道.对逆行轨道的演化性质进行分析时发现, 逆行轨道通常为准周期轨道, 比顺行轨道更加稳定.利用近拱点庞加莱映射可以快速确定不同类型轨道对应的初始状态, 为特定任务需求下的轨道设计提供了一种快速而有效的选择方案.     

勒让德伪谱法求解三维刚体摆姿态运动最优控制问题

研究伪谱法求解三维刚体摆姿态运动最优控制问题. 针对三维刚体摆这类含有约束的力学模型,提出了基于勒让德伪谱法的三维刚体摆姿态最优控制方法. 利用插值逼近设计了三维刚体摆姿态运动最优控制算法,得到了三维刚体摆的姿态最优控制轨迹,并结合松弛参数来控制插值点的取值,寻找满足的可行解. 仿真结果表明,基于勒让德伪谱法的最优控制算法使得三维刚体摆能以较小的误差运动到期望的末端姿态,且计算速度快,能够获得精度较高的控制输入量.     

多体系统轨迹跟踪的瞬时最优控制保辛方法

随着近年来机器人在各行业领域的广泛应用,对机器人的动力学与控制性能不断提出新的要求,特别是对设计越来越复杂、操作越来越灵巧的智能机器人,要求其能够对目标轨迹实现高精度跟踪以满足实际工作需求. 因此,针对机器人多体系统对目标轨迹跟踪的任务需求,基于微分代数方程提出瞬时最优控制保辛方法. 首先,采用多体动力学绝对坐标建模方法建立机器人系统的普适动力学方程,即微分代数方程;然后,采用保辛方法将连续时间域内的微分代数方程进行离散化,进而得到以当前位置、速度和拉式乘子为未知量的非线性代数方程组;其次,通过引入对目标轨迹跟踪以及对控制加权的瞬时最优性能指标,根据瞬时最优控制理论获得当前最优控制输入;最后,通过离散时间步的更新完成对目标轨迹的跟踪任务. 为了验证本文方法的有效性,以双摆轨迹跟踪控制为例进行了数值仿真,结果表明:针对机器人轨迹跟踪任务所提出的瞬时最优控制保辛方法能够实现对目标轨迹的高精度跟踪,且瞬时最优控制由受控微分代数方程推导获得,更具一般性,能够适应其他复杂多体系统的轨迹跟踪控制问题.      

HYBRID OPTIMIZATION STRATEGY FOR MOTION PLANNING OF SPACE ROBOT SYSTEM WITH PRISMATIC JOINT

研究了自由漂浮带滑移铰空间机器人非完整运动规划的最优控制问题,提出一种由高斯伪谱法求解可行解与直接打靶法求解最优解相结合的混合优化策略.首先,根据多体系统动力学理论建立空间机器人的动力学模型,给定系统的初始和目标位形,将空间机器人运动规划问题描述成博尔察（Bolza）型最优控制问题;然后,利用高斯伪谱法将最优控制问题离散为非线性规划问题,求解在较少勒让德-高斯（Legendre-Gauss,LG）点时状态变量和控制变量对应的可行解;最后,在LG点处离散控制变量,作为直接打靶法的初值,利用序列二次规划算法求解空间机器人系统的优化运动轨迹和最优控制输入.通过数值仿真,系统优化运动轨迹光滑平稳,最优控制输入也能很好地满足各种约束条件,仿真结果验证了该混合优化策略的鲁棒性和有效性.     

基于hp自适应伪谱法的组合动力可重复使用运载器轨迹优化

针对传统化学火箭难以重复使用、发射成本高的问题,提出了一种水平起飞/降落的新一代可重复使用运载器飞行方案,并对其动力模式设计和上升段轨迹优化方法等关键技术进行深入研究。设计了一种涡轮\冲压发动机结合火箭\冲压发动机的组合动力模式,建立了发动机推力与高度、马赫数等变量间的耦合模型,根据动力形式将上升段轨迹分为两段并采用全局搜索法确定动力切换的最佳时机。根据分段结果,分别以燃料最省和终端速度最大为指标,利用hp自适应伪谱法对两段进行轨迹优化设计。该算法基于双层策略求解最优控制问题,兼备伪谱法和有限元法的优点,与打靶法、伪谱法和间接法相比,初值更易选取,收敛速度更快。     

Multi-swingby optimization of mission to Saturn using global optimization algorithms

Based on the trajectory design of a mission to Saturn, this paper discusses four different trajectories in various swingby cases. We assume a single impulse to be applied in each case when the spacecraft approaches a celestial body. Some optimal trajectories ofEJS, EMS, EVEJS and EVVEJS flying sequences are obtained using five global optimization algorithms： DE, PSO, DP, the hybrid algorithm PSODE and another hybrid algorithm, DPDE. DE is proved to be supe- rior to other non-hybrid algorithms in the trajectory optimi- zation problem. The hybrid algorithm of PSO and DE can improve the optimization performance of DE, which is vali- dated by the mission to Saturn with given swingby sequences. Finally, the optimization results of four different swingby sequences are compared with those of the ACT of ESA.     

柔性双尾鳍效率的测量和优化

建立了仿生拍动推进推力和效率的测量平台,通过电机的转动扭矩、转动角速度及推进力计算得到仿生推进系统的推进效率.对本课题组研制的组合仿生无人水下航行器中柔性双尾鳍模型进行了测量与分析,结果表明柔性尾鳍推进时小摆幅摆动运动效率要高于大摆幅运动;双尾鳍推进方式在小摆幅模式下无论在推进力还是在效率方面都要优于单尾鳍.最后根据模型实验数据,优化了双尾鳍推进的运动参数,并用于优化UUV性能.     

推力电磁轴承的电磁场分析

利用有限元法对推力电磁轴承的磁场分布进行了分析计算,并通过对某一实际电磁轴承的推力轴承磁场分布及漏磁的分析计算,得出了推力轴承和推力盘之间的气隙及推力轴承和转子之间的气隙之比的最佳设计值,所得结果可作为推力电磁轴承实际设计的准则     

改进型扩展比例边界有限元法

结合了扩展有限元法(extended finite elementmethods,XFEM)和比例边界有限元法(scaled boundary finite elementmethods,SBFEM)的主要优点,提出了一种改进型扩展比例边界有限元法(improvedextended scaled boundary finite elementmethods,$i$XSBFEM),为断裂问题模拟提供了一条新的途径.类似XFEM,采用两个正交的水平集函数表征材料内部裂纹面,并基于水平集函数判断单元切割类型;将被裂纹切割的单元作为SBFE的子域处理,采用SBFEM求解单元刚度矩阵,从而避免了XFEM中求解不连续单元刚度矩阵需要进一步进行单元子划分的缺陷;同时,借助XFEM的主要思想,将裂纹与单元边界交点的真实位移作为单元结点的附加自由度考虑,赋予了单元结点附加自由度明确的物理意义,可以直接根据位移求解结果得出裂纹与单元边界交点的位移;对于含有裂尖的单元,选取围绕裂尖单元一圈的若干层单元作为超级单元,并将此超级单元作为SBFE的一个子域求解刚度矩阵,超级单元内部的结点位移可通过SBFE的位移模式求解得到,应力强度因子可基于裂尖处的奇异位移(应力)直接获得,无需借助其他的数值方法.最后,通过若干数值算例验证了建议的$i$XSBFEM的有效性,相比于常规XFEM,$i$XSBFEM的基于位移范数的相对误差收敛性较好;采用$i$XSBFEM通过应力法和位移法直接计算得到的裂尖应力强度因子均与解析解吻合\较好.      

Low-thrust trajectory optimization in a full ephemeris model

The low-thrust trajectory optimization with complicated constraints must be considered in practical engineering. In most literature, this problem is simplified into a two-body model in which the spacecraft is subject to the gravitational force at the center of mass and the spacecraft＇s own electric propulsion only, and the gravity assist （GA） is modeled as an instantaneous velocity increment. This paper presents a method to solve the fuel-optimal problem of low-thrust trajectory with complicated constraints in a full ephemeris model, which is closer to practical engineering conditions. First, it introduces various perturbations, including a third body＇s gravity, the nonspherical perturbation and the solar radiation pressure in a dynamic equation. Second, it builds two types of equivalent inner constraints to describe the GA. At the same time, the present paper applies a series of techniques, such as a homotopic approach, to enhance the possibility of convergence of the global optimal solution.     

TRAJECTORY TRACKING CONTROL AND VIBRATION SUPPRESSION OF SPACE-BASED ROBOT WITH ELASTIC FOUNDATION

讨论了基座存在弹性情况下,载体位置无控、姿态受控的漂浮基空间机械臂惯性空间轨迹跟踪控制及基座弹性振动主动抑制问题。由系统位置几何关系及动量守恒关系,建立了系统运动Jacobi 关系;之后利用拉格朗日方法并结合系统动量守恒关系建立了系统动力学方程。基于奇异摄动理论的两种时间尺度假设,将该方程分解为描写系统刚性运动的慢变子系统与描写系统弹性振动的快变子系统。对慢变子系统设计了基于计算力矩法的轨迹跟踪控制器;对于快变子系统则设计了线性二次最优控制方案。数值仿真证实了提出的控制方法的有效性。     

具有内域的双层加筋圆柱壳动响应特性简

采用外域双渐近法模拟外部流体域,内域双渐近法模拟内部流体域,采用非线性有限元软件ABAQUS 模拟结构,建立流固耦合数值模型并验证了其有效性. 利用建立的数值模型研究内部流体对双层加筋圆柱壳动响应的影响. 研究结果表明：一阶双渐近的解在中后期振荡周期变短、幅值变小,而二阶双渐近的解与解析解吻合良好;内部流体的存在减弱了双层加筋圆柱壳外壳的塑性变形,对其内壳塑性变形的影响较小;内部流体的存在减弱了双层加筋圆柱壳的速度响应.     

两串列扑翼的相位差对平均推力影响机理的实验研究

在一个低雷诺数的循环水洞中,实验研究了前后翅翼之间的相位差对两串列扑翼平均推力的影响.利用一个三分量的Kistler 压力传感器来测量扑翼的瞬时力;利用一个数字粒子测速仪系统(TSI DPIV) 来测量扑翼的前缘涡以及其周围的流场. 当相位差从0° 增加到360°,前翅的平均推力随着相位差正弦变化;前翅平均推力的增加是由于后翅的前缘涡和滞止区域增加了前翅的有效攻角. 后翅平均推力曲线有一个明显的V 字形低谷.低谷处较小的平均推力是由于前翅的脱落涡抑制了后翅前缘涡的形成并且减小了其有效攻角.当间距为0.5倍弦长相位差约为290°时,前后翅翼平均推力系数的合值能达到最大值0.667,明显大于两倍的单翼平均推力系数(2×0.255).      

BILU implicit multiblock Euler/Navier–Stokes simulation for rotor tip vortex and wake convection

In this paper, a block incomplete lower–upper (BILU) decomposition method is incorporated with a multiblock three‐dimensional Euler/Navier–Stokes solver for simulation of hovering rotor tip vortices and rotor wake convection. Results of both Euler and Navier–Stokes simulations are obtained and compared with experimental observations. The comparisons include surface pressure distributions and tip vortex trajectories. The comparisons suggest that resolution of the boundary layer is important for the accurate evaluation of the blade surface loading, but is less so for the correct prediction of the vortex trajectory. Numerical tests show that, using Courant–Friedrichs–Lewy (CFL) number of 10 or 30 with the developed BILU implicit scheme can be 6–7 times faster than an explicit scheme. The importance of solution acceleration schemes that increase the permitted time‐step is illustrated by comparing the evolving wake structures at different stages of the calculation. In contrast to fixed wing simulations, the extent of the wake structures is shown to require resolution of large physical time. This observation explains the poor performance that is obtained when employing convergence acceleration strategies originally intended for solution of equilibrium problems, such as the multigrid methods. Copyright © 2007 John Wiley & Sons, Ltd.