三维弹道有理Bezier曲线造型与优化方法

为有效克服传统弹道优化方法解的收敛性和最优性受搜索算法、初始猜测等影响大的问题,提出了基于有理Bezier曲线的三维弹道造型与优化计算方法。根据边界条件用光滑且只含少量造型变量的有理Bezier曲线形成待优化三维造型弹道,采用逆动力学方法计算造型弹道倾角、弹道角速度和法向加速度等,结合弹道其它动力学方程和约束计算造型弹道性能指标,通过对弹道造型变量寻优得到最优弹道。这种方法将无限维弹道优化问题转换为对很少造型变量的参数优化问题, 而且得到的最优弹道完全光滑可飞。与间接法相比,无需求解两点边值问题使用方便、鲁棒性强;与直接法相比,不需要对动力学方程离散化、寻优变量少、求解收敛性好,而且解的全局最优性和光滑性更高。仿真结果及与自适应伪谱法的比较校验了方法的实用性和有效性。     

多目标融合的冗余空间机械臂碰前轨迹优化

面向在轨对接与装配任务,提出了多目标融合的冗余度空间机械臂碰前轨迹优化方法。通过分析任务特点,阐述了空间机械臂碰前的多重任务目标。针对由空间机械臂连续轨迹跟踪、碰撞姿态控制及碰撞冲量幅值优化三个目标组成的第一类任务,通过设置任务优先级融合多目标并采用基于主任务零空间的方法实现对碰前轨迹的优化;而对由空间机械臂点到点规划、碰撞姿态控制及碰撞冲量幅值优化三个目标组成的第二类任务,通过设置任务权重将多目标融合并采用遗传算法对碰前轨迹进行优化。以典型七自由度机械臂为例开展仿真实验,验证了提出的轨迹优化方法的有效性。     

太阳帆日心定点悬浮转移轨道设计

研究了太阳帆航天器日心定点悬浮轨道(HFDO)的转移轨道设计问题,以球坐标形式建立了太阳帆的动力学模型,基于该模型给出在日心悬浮轨道基础上实现定点悬浮的条件,提出了一种实现日心定点悬浮的转移轨道设计方法。首先,确定定点悬浮的位置;然后,设计经过该位置的绕日极轨轨道;最后,实施轨道减速实现定点悬浮,并给出了解析形式的轨道控制律。结合太阳极地观测任务,设计了定点悬浮在太阳北极1AU处的太阳帆转移轨道。仿真结果表明:该轨道转移方案总耗时3.5年,太阳帆定点到黄北极距日心1AU处,此后只要保持太阳光垂直照射帆面,即可维持稳定的悬浮状态。     

Rapid design and optimization of low-thrust rendezvous/interception trajectory for asteroid deflection missions

Asteroid deflection techniques are essential in order to protect the Earth from catastrophic impacts by hazardous asteroids. Rapid design and optimization of low-thrust rendezvous/interception trajectories is considered as one of the key technologies to successfully deflect potentially hazardous asteroids. In this paper, we address a general framework for the rapid design and optimization of low-thrust rendezvous/interception trajectories for future asteroid deflection missions. The design and optimization process includes three closely associated steps. Firstly, shape-based approaches and genetic algorithm (GA) are adopted to perform preliminary design, which provides a reasonable initial guess for subsequent accurate optimization. Secondly, Radau pseudospectral method is utilized to transcribe the low-thrust trajectory optimization problem into a discrete nonlinear programming (NLP) problem. Finally, sequential quadratic programming (SQP) is used to efficiently solve the nonlinear programming problem and obtain the optimal low-thrust rendezvous/interception trajectories. The rapid design and optimization algorithms developed in this paper are validated by three simulation cases with different performance indexes and boundary constraints.     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

蜂窝夹层板超高速撞击极限方程分析

为分析撞击极限方程对蜂窝夹层结构的穿透特性的预测能力,调研得到了3类撞击极限方程的表达形式和等效方法,以及131个采用碳纤维复合材料(CFRP)面板的蜂窝夹层板结构的试验数据,并对撞击极限方程的预测能力进行了比较计算。结果发现,MET方程对他源数据的四种(未失效、失效、总体以及安全)预测率均大于80%,进行在轨航天器结构的失效分析时可优先选用;SRL方程对本源数据的安全预测率达到了100%,在他源数据上的安全预测率也很高,适用于航天器防护结构设计。探讨了撞击极限方程中的系数、速度分界值的优化思路,以提高撞击极限方程的预测率。     

吸气式高超声速推进助推段内流振荡及其抑制

为了掌握高超声速飞行器盲腔流场的激波振荡机理,采用"双时间步"方法对轴对称吸气式高超声速飞行器内流道盲腔流场进行了非定常数值模拟。结果表明,飞行器内流道出现了非定常激波振荡现象,内流道壁面压力出现了低频高幅的周期性变化。为了消除高超声速飞行器内流道的激波振荡现象,设计了级间段泄流方案,对飞行器进行了内外流一体化定常和非定常的数值模拟。结果表明,级间段泄流方案可以有效地消除内流道的激波振荡现象。     

三阶幻方问题的代数解法

介绍了中国古代数学家甄鸾和杨辉关于三阶幻方的一个构造方法。用线性代数的方法构造三阶幻方约束方程组时,当自由未知量分别取1,2,3,…,9时,对应的方程组共有56组整数解,为了找出幻方的全部解,需要确定自由未知量取值的限制条件,由此对方程组中的自由未知量的取值进行了分类探讨,进而求出三阶幻方的全部解。利用线性代数的方法去探求三阶幻方的解法种数显然不是最简捷的,但却有完备的理论做保证。此外,提供了一个如何利用线性代数知识来解决实际问题的实例,并简介了三阶幻方所具有的奇妙的特性。     

各向异性Chaboche黏塑性本构方程隐式应力积分算法

针对各向异性Chaboche黏塑性本构方程,利用有效参数控制法推导了该方程的隐式积分算法,该算法把所有未知量都用同一个控制参数表示,并形成控制方程,从而把本构方程应力积分的问题转化成一个非线性控制方程的求解问题.该积分算法程序,与通用有限元软件Marc结合,通过单个单元对单向拉伸、循环、松弛的响应的实验数据进行了数值模拟,通过与显式算法对比,该隐式算法的稳定性和收敛性有所提高.构件的对比分析也表明隐式算法比显式算法具有较高的效率.      

隐式耦合求解燃烧室内喷雾两相流场

在任意非结构网格下,采用开源程序对某型涡喷发动机环形燃烧室内的喷雾燃烧过程进行数值模拟.首先选用k-ε双方程模型、EBU(eddy break up)-Arrhenius紊流燃烧模型、颗粒随机轨道模型等模拟流场中的物理过程.然后用隐式耦合算法求解气相流动、拉格朗日方法跟踪液滴相在流动区域的运动,并考虑气相与液滴相的完全双向耦合作用.仿真实验得到的 燃烧室内温度分布、马赫数分布、速度分布、压力分布均较合理,温度分布与文献中的基本一致,这表明隐式耦合方法能够较好地模拟燃烧室内流场,反映其内部流动特点.