二维翼型结冰数值计算

翼型结冰计算是随着结冰时间步长变化的准定常过程,因此在每个时间步长后都要更新各个模块的计算工作,如此反复直到达总结冰计算时间.在各模块计算中,分别采用了泊松方程法生成网格、求解雷诺时均方程计算流场、拉格朗日追踪法求解水滴轨迹、基于Messinger模型计算冰形、采用变距滤点法重构边界.引入动态链接库技术将各模块集成为结冰计算平台,并将三种冰形(楔形冰、槽状冰及混合冰)计算结果与国外软件及试验数据进行了对比.      

基于直接配点法的滑翔轨迹快速优化设计

介绍了基于五次Gauss-Lobatto多项式的直接配点法在再入飞行器三维轨迹最优化问题中的应用。首先给出了再入飞行器轨迹优化问题模型,其中运动方程为三自由度模型,性能指标选为到达指定地点飞行时间最短,控制变量则为无量纲升力系数和倾侧角。再入飞行过程中受到加热率、过载和动压约束,终端状态受到目标位置约束。然后,应用直接配点法将最优控制问题离散化为非线性规划问题,将动态优化问题转化为静态参数最优化问题。选取各节点和配点上的状态量和控制量作为优化参数。最后应用SNOPT软件包对参数最优化问题进行求解。仿真结果表明直接配点法对于再入飞行器轨迹初始参数取值不敏感,且求解过程具有一定的实时性。     

考虑转捩影响的风力机翼型气动特性计算研究

基于二维雷诺平均Navier-Stokes(RANS)方程,耦合eN数据库方法进行流动转捩判断,采用SA湍流模型计算了风力机专用翼型DU91-W2-250的气动特性,得到了雷诺数Re=1.0×106时,攻角在-11°~12°范围内该翼型的升力和阻力特性曲线及压力分布。通过对全湍流计算结果、耦合转捩判断的计算结果及实验数据的对比,表明考虑转捩因素后的计算结果与实验结果吻合更好,说明要准确模拟风力机翼型的绕流,必须考虑转捩影响。     

无人机对地自动攻击占位轨迹生成研究

基于无人机对地自动攻击过程,应用最小值原理求解了无人机对地自动攻击占位段的应飞轨迹,该方法将威胁模型、飞机运动方程和对地攻击任务有机结合,使生成的占位轨迹既能有效规避各种威胁又能为武器投放创造条件。对占位轨迹计算中的两点边值问题采用了共轭梯度法求解。通过实例仿真表明该方法是正确有效的,所得占位轨迹波动平缓没有剧烈折线现象,对占位机动研究具有一定的参考价值。     

高阶间断有限元法的并行计算研究

根据间断有限元法的数据结构特点,基于METIS网格分区技术,设计并行计算策略,在非结构网格上实现了并行高阶间断有限元法。控制方程的数值通量项使用Local Lax-Friedrichs(LLF)格式计算。设计了并行的牛顿-块高斯赛德尔法(Newton-Block GS)来加速收敛,提高迭代效率。并行性能分析表明,所设计的并行算法能够得到较好的加速比和并行效率,有效地节省计算时间,合理分配内存。这使得采用高阶间断有限元法计算更为复杂的问题成为可能。     

应用多块对接结构网格方法的直升机涵道尾桨气动特性分析

基于湍流Navier-Stokes方程,建立了一套直升机涵道尾桨流场及气动特性的CFD分析方法.在该方法中,针对直升机涵道尾桨流场的特点以及动量源方法对网格系统的要求,提出了一种多块对接网格生成方法,分块贴体网格采用求解Poisson方程的方法生成.通过结合涵道尾桨的运动方式、几何特征及气动特征,建立一个包含动量源项的N-S方程的涵道尾桨流场计算方法和迭代流程;为较好地捕捉涵道壁附近存在的旋涡及分离现象,采用了一方程S-A湍流模型.通过对ONERA M6机翼的绕流特性的计算分析,以及对涵道尾桨在悬停状态下的拉力、诱导速度、载荷分布等的计算分析,验证了该CFD方法的有效性.在此基础上,采用该方法进行直升机悬停、侧飞和前飞状态下涵道尾桨流场与气动特性的数值分析,得到了关于涵道尾桨流场和气动特性的一些有意义的结论.     

基于轨道方程线性化的中段飞行轨道修正方法

对于深空探测、大范围远程交会、抵近观察等中段飞行时间较长的航天任务,中段轨道修正是完成任务的重要保证.基于摄动制导的思路,通过对轨道方程进行线性化,研究了二体条件下的中段轨道修正方法.在近地航天任务中,各种摄动力特别是J2项对终端位置精度的影响不可忽略,因此研究了考虑J2摄动影响时线性化方程的补偿方法.最后通过数值仿真...     

悬停状态直升机桨叶扭转分布的优化数值计算

建立了一套基于高精度计算流体力学(CFD)技术和代理模型优化算法的旋翼气动外形设计方法。在该方法中,旋翼流场气动性能的计算采用了基于Navier-Stokes/Euler方程的CFD方法,并根据流场特点、精度和效率的要求采用Baldwin-Lomax(B-L)湍流模型,通量计算采用Roe-MUSCL格式进行。为提高网格生成质量和便于流场控制方程的求解,将流场分成两个区域,即围绕旋翼的黏性区和无黏的背景网格区。其中,桨叶网格使用了基于二维翼型网格的参数化方法生成,数值计算结果表明该方法有效地提高了网格生成质量及效率。在参考旋翼流场及桨叶细节流动分析的基础上给出设计变量及范围,有效减小了优化问题的规模;为满足优化和机理分析的需要,将基于置换遗传算法优化的拉丁超立方(PermGA LHS)方法和径向基函数(RBF)的代理模型优化方法引入到桨叶外形的优化设计中。首先以Helishape 7A旋翼为算例,检验了数值模拟方法的准确性。然后,应用所建立的优化方法针对旋翼负扭转分布进行了优化计算,结果表明优化后的旋翼悬停气动性能比优化前有了明显提高。     

基于速度增量的空间绳系机器人中距离逼近过程最优轨迹规划

针对空间绳系机器人中距离逼近过程最优轨迹规划问题,提出了基于速度增量的多目标逼近轨迹优化方法,优化指标为总速度增量及逼近时间。首先建立逼近过程相对动力学模型及最优逼近轨迹优化模型,然后利用改进型非劣分类遗传算法得到相对逼近距离1.5 km内逼近轨迹的Pareto最优解。仿真结果表明,该方法可以揭示空间绳系机器人逼近距离1.5 km内逼近时间、燃料消耗、相对目标的面内视界角及速度增量次数之间的相互关系,能满足针对不同任务需求提供相应最优轨迹的要求。     

Footprint Problem with Angle of Attack Optimization for High Lifting Reentry Vehicle

A formal analysis to footprint problem with effects of angle of attack (AOA) is presented. First a flexible and rapid standardized method for footprint generation is developed. Zero bank angle control strategy and the maximum crossrange method are used to obtain virtual target set; afterward, closed-loop bank angle guidance law is used to find footprint by solving closest approach problem for each element in virtual target set. Then based on quasi-equilibrium glide condition, the typical inequality reentry trajectory constraints are converted to angle of attack lower boundary constraint. Constrained by the lower boundary, an original and practical angle of attack parametric method is proposed. By using parametric angle of attack profile, optimization algorithm for angle of attack is designed and the impact of angle of attack to footprint is discussed. Simulations with different angle of attack profiles are presented to demonstrate the performance of the proposed footprint solution method and validity of optimal algorithm.