防空导弹弹体弹性振荡的抑制

导弹在飞行过程中，除了质心运动和绕质心的角运动，还存在沿纵轴方向的结构弯曲，即弹性振荡。本文讨论了弹体的弹性振荡通过敏感元件耦合对弹上控制回路稳定性的影响，分析了两种抑制弹性振荡的方法，并通过实例仿真比较了两种方法各自的特点。     

基于欧拉方程的二维振荡机翼非定常气动设计反命题方法

建立了基于欧拉方程组的二维振荡机翼非定常气动设计反命题方法的数学模型.通过一系列变换,将物理时空中的求解域转换成映象坐标系中的规范区域;并导出了映象坐标下的欧拉方程,结合非定常反命题的边界条件,便可用有限差分方法求解.本方法引入映象坐标系解决反命题边界形状的不确定性,并能利用欧拉方程现有的各种差分格式于反命题求解.     

隔离段内正激波串受迫振荡特性研究

为了解隔离段在真实燃烧室反压作用下的工作特性,通过数值方法模拟燃烧室振荡燃烧引起的脉动反压,在来流Mai=2条件下,探讨了脉动反压引起的激波串受迫振荡特性。结果表明,在反压脉动条件下,激波串前缘位置周期振荡,其振荡频率与反压脉动频率一致,振幅与反压脉动频率负相关,反压脉动频率600Hz时,振幅仅为管道长2.97%;就抗反压能力而言,在反压脉动条件下,隔离段能承受比定常状态更大的反压峰值,且脉动频率越大,可承受瞬态反压峰值越大。激波串的受迫振荡包含两种形式:激波串整体位置的前后运动和激波节之间的相对运动。在一定程度内,反压脉动频率较低时(50Hz,100Hz),激波串整体位置的前后运动是主要的受迫振荡形式,激波节之间的相对运动很弱;随反压脉动频率升高(300Hz,600Hz),激波节之间的相对运动加剧,逐渐成为激波串受迫振荡的主要形式。另外,在激波串振荡过程中,激波串往上下游运动经过同一位置时,激波串形态会出现迟滞,这是导致总压恢复系数迟滞的原因。     

低速轴流压气机旋转失速边界的模型描述

从低速轴流压气机完整旋转失速过程所具有的突变性和迟滞性等典型动力学特征出发,通过对其拓扑性质的对比分析,证明压气机旋转失速与突变理论尖点突变模型在拓扑结构上具有相似性;以拓扑不变性定理为基础,通过发展有效的拓扑映射方法,建立了尖点突变模型的突变点集与压气机旋转失速突变点集之间的拓扑映射关系.在两台低速轴流压气机上的模型实际应用效果表明:尖点突变模型能够定性地描述具有突变型失速特点的低速压气机不同转速工况下完整旋转失速过程的突变性、迟滞性以及分岔特性,通过该发展的拓扑映射方法,可以将实验获取的失速点和恢复点完全映射到模型线上,这表明尖点突变模型能较为准确地描述旋转失速的失速边界.      

随机复合材料统计多尺度边界元算法研究

热传输行为是空天飞行器热防护系统服役过程中的重要问题,研究随机复合材料的热传输机理,可为热防护系统材料与结构的一体化设计提供理论依据。针对随机复合材料的热传导问题,发展一种统计多尺度边界元算法,首先建立等效材料参数的统计多尺度分析模型,然后给出统计意义下等效参数的多尺度边界元预测算法,并通过与理论结果的对比来验证算法的有效性,最后研究微结构分布状态对陶瓷多孔材料等效导热系数的影响。结果表明:采用统计多尺度模型及多尺度边界元算法预测随机复合材料的热传导性能是有效的。     

混合边界子结构的航空发动机机匣模型修正

提出了一种基于混合边界模态综合的复杂结构有限元模型修正方法。其主要步骤包括:①子结构划分,根据结构形式划分待修正区域,得到子结构和残余结构;②缩聚和装配,利用混合边界模态综合法,将子结构内部自由度集缩聚至混合边界自由度集,得到子结构缩聚矩阵,并与残余结构的系统矩阵进行装配;③修正,基于灵敏度分析方法,对装配后的残余结构进行参数修正。将该方法应用于航空发动机外机匣的精细化有限元建模及模型修正研究,针对局部连接结构参数修正,子结构模型修正方法的参数收敛后最大误差为064%,其计算效率提高了515倍。算例结果表明,该方法在保证修正精确度的同时,能提高大规模复杂结构有限元模型修正的计算效率。      

基于高阶WENO格式的旋翼非定常涡流场数值模拟

为提高直升机旋翼黏性涡流场计算流体力学（CFD）的计算精度及降低数值耗散,发展了一套基于五阶加权本质无振荡（WENO）格式的旋翼非定常涡流场数值计算方法。采用运动嵌套网格方法生成围绕前飞状态旋翼的网格系统,以Navier-Stokes方程作为主控方程,湍流模型采用一方程Spalart-Allmaras模型。为了提高旋翼流场中对涡的形成、演化等发展过程的模拟精度,采用高间断分辨率的Roe-WENO格式计算对流通量。对状态变量的插值,选取适当的加权因子,通过多个重构模板的凸组合来构造五阶WENO格式,在交界面附近获得具有五阶精度的状态变量。为了提高计算效率,采用高效的隐式LU-SGS（Lower-Upper Symmetric Gauss-Seidel）推进方法和并行计算策略。最后,运用该方法分别对悬停状态下C-T（Caradonna-Tung）旋翼和UH-60A旋翼以及前飞状态下C-T旋翼和SA349/2旋翼的涡流场及气动特性进行了数值模拟,并将数值结果与传统二阶精度格式的计算值进行了对比。结果表明：在计算时间增长10%～20%的前提下,五阶WENO格式能够捕捉到更精确的涡流场特性,反映了五阶WENO格式在计算旋翼非定常涡流场时具有更高的计算精度与更低的数值耗散特性。     

喷注方案对凹腔稳焰燃烧室低频振荡的影响研究

为了研究在入口来流马赫数2.52,总温1486K的超声速来流条件下,稳焰凹腔上游不同位置乙烯横向喷注对模型发动机燃烧室内低频燃烧振荡特性的影响,通过1kg/s直连式超燃试验平台,利用高频压力传感器、高速摄影相机等设备,对凹腔上游近距离、远距离喷注等方案的发动机内部压力与火焰动态特性进行了研究。试验结果表明:在当前当量比条件下,当稳焰凹腔上游近距离喷注燃料时,燃烧室存在较大范围亚声速区域,并出现由热声不稳定性激励的低频压力振荡,频率分布范围较宽(50～400Hz)且振幅较弱。对于燃料喷注位置到稳焰凹腔距离较远的情况,燃烧室内出现以火焰逆传和火焰吹脱为特征的周期性火焰振荡现象。分析认为较远喷注距离有利于燃料-空气充分混合并形成预混区,导致火焰快速逆传。火焰逆传与DDT (爆燃转爆震)中的火焰加速传播过程有关。周期性火焰逆传与火焰吹脱过程相耦合形成了具有特定主频(约85Hz)且振幅较大的低频压力振荡。     

飞行器双自由度耦合振荡的数值模拟方法

针对模拟双自由度(2DOF)运动的"PQR"强迫振荡试验装置,通过坐标变换,将机构运动转化为空间欧拉角进行表示,实现对机构俯仰/偏航、俯仰/滚转和偏航/滚转等3种耦合运动的统一描述,采用预定运动轨迹的动网格计算技术,构造了基于非定常雷诺平均Navier-Stokes(URANS)方程的飞行器双自由度振荡数值模拟方法。围绕某复杂构型飞行器低速流动,在模型迎角10°、滚转角±40°、偏航角±40°振荡预设条件下分别对单自由度偏航、滚转振荡和双自由度偏航/滚转振荡进行了模拟,得到的气动系数迟滞回线在形态和量值上与风洞试验结果吻合,证明当前方法可以作为飞行器复杂耦合运动的一种有效研究手段。针对偏航/滚转耦合振荡,从运动形式和气动特性上与单独偏航、滚转运动进行了对比分析,结果表明,耦合运动在气流角、绕体轴的旋转角速度等方面与单自由度叠加效果不同,力矩系数迟滞回线存在曲线交叉的现象,表现出与单自由度振荡不同的阻尼特性。在当前模拟状态下,偏航/滚转振荡流场以横向分离涡运动为主要特征,可以推测大迎角情况下纵横向分离涡流动结构更加复杂,耦合作用更强,需要进一步在分离流模拟方法上开展研究。     

飞机主起落架舱门变形分析与边界约束优化

基于有限元软件Patran和Nastran,对一种常见的短程民用飞机主起落架舱门在空中严重工况下的变形进行仿真分析,并通过试验验证有限元模型的准确性;最后,根据有限元仿真分析结果对舱门的边界约束进行优化,并将优化后的舱门变形情况与优化前进行对比分析;分析结果表明,优化边界约束后的舱门变形量显著减小,其结构应力水平明显降低。