回转表面上非测地线缠绕方程及其解

本文利用微分几何推导了回转表面上的非测地线缠绕方程,并用Runge-Kutta-Fehlberg法求解。此法与文献的分段圆锥法对比,两者计算结果很一致,说明在工程中都是可用的。     

非线性飞行稳定性研究的新综合方法

本文提出“分歧和突变理论-时间历程-龙格库塔灵敏度分析”的综合方法,用其研究作为非线性系统的飞机稳定性。该方法利用分歧理论来确定飞机的全局稳定性,并由平衡面找出应予检验的平衡点,对这些平衡点,计算状态参数的时间历程,由此确定飞机运动的性质。最后,利用龙格-库塔灵敏度分析方法来研究动稳定性随时间的变化,由此可以解析地求得任意瞬间影响飞行稳定性的空气动力项,结构项及组合项。 利用上述方法所得到的有关稳定性的信息是全面的,它对飞机的操稳设计有重要参考价值。     

低温气动阀动态响应特性的非线性分析方法

为了在液体火箭发动机试验过程中缩短低温气动阀的响应时间,减少非线性求解阀门动态响应的计算量,提升试验过程中的安全性与可靠性,根据低温气动阀的工作原理,构建阀门的动态响应非线性传递模型,并在此基础上,通过引入变步长系数与定义遗忘因子,在四阶龙格-库塔法基础上,构建了一种动态递推龙格-库塔法的非线性动力学求解方法,并提出响应速度的改进方案.最后通过实例的验证与仿真,为提升低温气动阀响应速度提供了一种方法指导,为提高液体火箭发动机试验技术提供了支持.     

高分辨率格式在非定常无粘可压缩流动中的应用研究

本文发展了一种计算非定常无粘可压缩流动的高分辨率格式。通量计算采用AUSMPW+Rie-mann求解器,其中使用了基于单调性的压力加权函数f,并与Roe Riemann求解器进行比较。为提高精度,使用五阶WENO格式进行左右状态插值,且时间推进采用三阶TVD Runge-Kutta方法。用该方法对移动接触间断问题、Sod激波管问题、二维激波反射问题和双马赫反射问题进行了数值计算,数值结果表明基于五阶WENO插值的AUSMPW+格式有很高的激波及接触间断分辨能力,并比基于五阶WENO插值的Roe格式有更高的计算效率。     

基于非结构网格流场计算的网格重排序

提出基于非结构网格流场计算时的网格重排序,优化了网格信息在内存中的存贮位置,使得计算过程中调用数据的处理速度加快,有效提高了计算效率。以M6机翼及DLR-F4翼身组合体的跨音速无粘流场为算例,分别通过显式四步龙格-库塔推进格式和隐式LU-SGS推进格式来验证网格重排序后对计算效率的提高作用。原始网格均用Delaunay方法生长,并分别使用了两种重排序方法进行比较。计算结果表明:重排序后的计算效率可以提高到未排序时的3到5倍,很具实用价值。     

旋翼桨叶绕流Euler方程数值模拟

本研究应用格心格式有限体积法求解欧拉方程,模拟绕旋翼浆叶的流动。旋翼尾流的作用通过自由尾流分析来求解局部的诱导下洗速度以修正翼型攻角。选择O—H型式生成弦向和展向网格。通过算例计算为进一步寻求有效的Euler或NS方程数值解积累了经验。     

弹射救生数值仿真及不利姿态下救生性能分析

 国内弹射座椅救生性能预测分析主要依赖地面弹射试验。由于国内相关试验设备的缺乏,不利姿态下弹射座椅的救生性能很难通过试验获得,因此国内不利姿态下弹射座椅救生性能的分析研究比较薄弱。采用数值仿真的方法对弹射座椅的救生性能进行预测和分析,针对弹射座椅弹射过程的4个阶段分别建立了数学模型,其中自由飞阶段采用四元数法代替欧拉速率方程,以解决传统的六自由度方程出现的奇异性问题。采用四阶龙格-库塔法对数学模型进行求解,计算结果与地面弹射试验结果吻合较好。在此基础上,对某新型弹射座椅的弹射姿态轨迹进行计算,分析研究了弹射座椅在不利姿态下的救生性能,得到不利姿态下弹射座椅救生性能不佳的主要原因是座椅高速稳定性差以及飞机向下的牵连速度和飞机滚转角的影响,并根据分析结果提出了改善弹射座椅不利姿态下救生性能的基本方法。     

辛算法的纠飘研究

辛算法较RK(Runge-Kutta)方法,保持辛结构不变或保持哈密顿系统规律性不变是突出的优点,但点态数值精度并不理想.推导出了三阶、四阶辛算法的漂移量计算公式,通过补偿漂移量就能提高点态数值精度,既保辛结构又保证点态数值高精度,适合于工程应用.建立了分数步对称辛算法(即FSJS算法)的纠漂公式,制定了漂移的约束标准.相关算例的数值结果表明:三阶FSJS算法漂移量最小,点态数值精度更高.