再入飞行器高超声速绕流的一步和二步有限元计算

本文针对再入飞行器的空气动力问题，采用一步和二步有限元格式求解可缩守恒型Ｅｕｌｅｒ方程，并结合不同的人工粘性模型来捕捉激波，给出了高超声速二维钝头绕流的结果。     

钝锥大攻角超声速分离流场的数值模拟

最近,张涵信等人在传统的Beam-Warming隐式、无迭代、空间推进技术的基础上,根据边界层方程的性质,设计了一种可用小步长推进求解抛物化NS(PNS)方程、而不会引起解的“漂移现象”发生的方法。这种方法对轴对称流动的计算是成功的。本文就是将这一思想推广应用于大攻角有周向分离的流场计算。求解的区域为具有薄亚声速层的有粘与无粘干扰的整个激波层内的流场。在对攻角α=0°和α=20°的球钝锥的计算中,关于壁面上的压力、热流率及流场的涡旋结构均得到了满意的结果。文中特别研究了钝锥大攻角绕流的流动分离图象。 为了增强块三对角矩阵的主对角优势,通常在差分方程的左端附加二阶增量项。本文以选取适当小的推进步长的方法来达到增强主对角优势的目的,不需再附加二阶增量项,从而提高了解的精度。     

用欧拉方程计算机翼与翼身组合体绕流

本文用保角转绘加剪切变换生成三维机翼与翼身组合体Ｃ－Ｈ型网格，并在这类网格拓扑上，用有限体积法研制出可供分析机翼与翼身组合体绕流的三维欧拉方程计算程序。本方法的特色是改进了机翼表面网格点的分布，使机翼后缘与网格线一致；采用了当量机身，使机岙表面到对称面光滑过渡；嵌入了边界层粘性修正。改善了计算结果。数值计算表明，对粘性影响较大的超临界绕流，边界层修正效果显著。     

仿真与响应曲面法的结合及其在综合生产计划中的应用

 论述响应曲面法与仿真相结合的优化结构,着重讨论了一阶和二阶响应曲面的仿真拟合原理和寻优过程,从而可在仿真环境下给出最优决策参数。利用所开发的GASPRSM仿真优化软件,对飞机工厂生产计划安排进行了仿真运行,证明可以取得显著的经济效果。     

人工释能法求解线性代数方程组

本文建立求解线性代数方程组的人工释能法迭代格式,并给出该迭代格式收敛的充分必要条件和最佳释能时机,同时给出实现人工释能法的数值计算方法及其数值稳定性条件。     

基于μ分析的高超声速飞行器再入轨迹评估

为解决高超声速飞行器再入过程中存在的不确定性问题,将μ分析理论应用到标称再入轨迹的评估中.采用线性稳定储备准则,在密度、升力系数和阻力系数存在不确定性时,对标称轨迹在制导律下的制导稳定性进行了分析;根据所获得的最坏参数组合,采用系统的开环尼克尔斯曲线验证了μ分析的结果;分析了系统所能容纳的最大不确定性参数,并研究了系统对不同不确定性参数的敏感性.结果表明:在密度、升力系数和阻力系数不确定性中,阻力系数对该标称轨迹的鲁棒稳定性影响最为明显.     

A coupling VWM/CFD/CSD method for rotor airload prediction

A coupling fluid-structure method with a combination of viscous wake model(VWM),computational fluid dynamics(CFD) and comprehensive structural dynamics(CSD) modules is developed in this paper for rotor unsteady airload prediction. The hybrid VWM/CFD solver is employed to model the nonlinear aerodynamic phenomena and complicated rotor wake dynamics;the moderate deflection beam theory is implemented to predict the blade structural deformation; the loose coupling strategy based on the ‘delt method' is used to couple the fluid and structure solvers.Several cases of Helishape 7A rotor are performed first to investigate the effect of elastic deformation on airloads. Then, two challenging forward flight conditions of UH-60 A helicopter rotor are investigated, and the simulated results of wake geometry, chordwise pressure distribution and sectional normal force show excellent agreement with available test data; a comparison with traditional CFD/CSD method is also presented to illustrate the efficiency of the developed method.     

基于分数阶微积分理论的最优导引律设计

随着技术的发展及战场环境的日益复杂化,拦截机动目标的需求越来越大。然而传统制导律在拦截机动目标时存在制导精度差、末端过载突变的问题,故提出了一种基于分数阶微积分理论的最优导引律。首先,介绍了分数阶微积分的定义、性质及其数字实现方法;然后,分析了弹目相对运动关系,通过分数阶变阶次建模和最优控制理论推导出了分数阶导引律;最后,仿真结果表明：与传统比例导引法相比,所设计的分数阶最优导引律能够保持比例导引法良好的追踪性能且拦截时间能够缩短2s,末端过载值趋近于0。该方法解决了传统比例导引法在末端由视线角速率发散而导致的末端过载突变问题。     

基于线化稳定性理论的后掠翼边界层内横流稳定性研究

采用Mack方法求解忽略曲率和压缩性的三维线化稳定性方程（Orr-Sommerfeld方程）,针对无限展长后掠翼,计算分析边界层内不同波长横流驻波沿弦向的放大率,确定最不稳定波的波长。计算结果与参考文献中的实验结论和计算结果符合较好,表明该方法能够准确的预测出最不稳定波的波长,可以用来指导横流驻波主导下的后掠翼流动稳定性问题研究。     

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

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