基于Volterra级数降阶模型的气动弹性分析

采用Volterra级数进行非定常气动力辨识,建立了基于状态空间的非定常气动力降阶 模型,耦合结构方程,建立了降阶的气动弹性系统和气动伺服弹性系统,开展了颤振分析和 颤振主动抑制的初步研究,计算效率提高了2～3个数量级。通过气动弹性标模AGARD445.6和 气动伺服弹性标模BACT算例验证,证明降阶方法正确,可以提供高效、高精度的气动弹性 分析。
     

前缘钝度对高速翼型颤振边界的影响

高速飞行器控制面前缘不同的钝度对其气动特性具有影响,同时也影响了其气动弹性特性。采用了CFD/CSD耦合计算方法,研究3种不同钝度的前缘对2维翼型在Mach数为5时的颤振边界的影响。计算结果显示,随着翼型的钝度增加,颤振边界不断提高,即前缘钝度增加了气动弹性稳定性。     

一种在时域内计算三维翼颤振的方法

使用模态法构造结构模型,研究了在时域内耦合结构模型和气动模型预测三维翼颤振的方法;求解N-S方程计算非定常气动力.Newmark算法简化颤振方程.采用无限平面插值方法耦合结构网格和气动网格;使用弹性系数算法和局部网格重构法重构非结构气动网格,以确保结构方程和气动方程进行耦合求解.AGARD 445.6翼算例检验了计算方法的可行性.     

基于POD降阶模型的气动弹性快速预测方法研究

CFD/CSD耦合数值模拟是解决复杂气动弹性问题精度最高的方法,但同时也是计算效率最低的方法。研究了气动弹性系统的时域POD降阶模型方法,并引入平衡截断技术进一步降低时域POD/ROM的阶数,从而有效克服了时域POD/ROM阶数过高的缺点。以AGARD445.6机翼为例,说明了时域POD/ROM建模的各个细节,并将其用于气动弹性动响应及颤振边界的预测。计算结果表明,POD/ROM具有接近CFD/CSD耦合计算的精度,同时又大大提高了计算效率约1到2个量级。     

一种基于当地流活塞理论的超音速导弹气动伺服弹性分析方法

运用当地流活塞理论计算导弹的非定常气动力，在状态空间内实现了气动弹性系统与控制系统的耦合，并进行了时域分析。对某典型导弹的气动伺服弹性问题算例，计算了不同马赫数、不同迎角以及不同的传感器安装位置对导弹气动伺服弹性特性的影响，同时考察了控制系统的引入对不同刚度弹体气动伺服弹性特性的影响。算例结果与理论分析相吻合，表明该方法能简便的对全机、全弹复杂外形的气动伺服弹性正问题进行高效准确的分析，在较大范围的马赫数和迎角内都能得到可靠的结果。     

一种新的界面映射推进方法及其在气动弹性力学中的应用

非线性气动弹性研究中,涉及到非线性的流体动力学(CFD)和非线性的结构动力学(CSD )耦合问题,在耦合界面上不仅要进行两场之间信息的传递,而且要有匹配的时间推进格式 。针对计算流体动力学（CFD）和计算结构动力学（CSD）的耦合计算方法,基于边界能 量守恒,发展了一种新的界面处理方法,包括：① 界面映射方法,该方法基于耦合边界上 局部的网格信息,将两场之间载荷信息和位移信息的转换放在同一个映射矩阵中来处理,克 服了占用大量CPU时间和内存的需求;② 界面推进方法,该法基于松耦合计算流程,通过引 入半步交错推进,达到了二阶时间精度,提高了计算精度和效率。最后将该界面处理方法应 用于柔性大展弦比机翼的气动弹性计算和AGARD445.6机翼的动响应分析中。结果表明该方法 能够高效、高精度地处理不同物理场之间的数据交换和时间推进,并具有处理复杂几何体非 线性气动弹性问题的能力。     

用于快速仿真优化的改进差分进化算法及其应用

提出一种改进的差分进化算法,采用一种“位变”方法计算收缩因子,该方法首先根 据适应值对种群排序,然后根据各个体的排列位置确定收缩因子;采用正态分布函数对算法 参数进行随机扰动来维持种群的多样性;该算法还提出一种新的变异算子,并将其与基本的 差分变异算子结合使用以提高算法的寻优精度。经过对多个Benchmark函数的测试、分析和 比较,结果表明该算法具有较高的收敛精度和较快的收敛速度。最后将该算法用于火箭发动 机涡轮气动优化,以较小的计算成本将涡轮气动效率提高了2.5％。应用结果表明该算法 适用于快速仿真优化问题,能有效地节约计算成本。
     

吸气式高超声速飞行器鲁棒反演控制器设计

针对吸气式高超声速飞行器气动／推进／结构弹性耦合控制问题,提出了鲁棒反演控制器设计方法。采用反演和动态逆方法设计虚拟控制量和实际控制量,通过引入一阶低通滤波器来获取虚拟控制量的导数,解决了虚拟控制量求导复杂问题;为了增强控制器的鲁棒性,采用充分光滑投影算子对模型非匹配不确定项进行估计和补偿,同时避免了可能出现的参数漂移问题。仿真结果表明,该控制器对模型气动参数拟合误差、攻角和升降舵偏角摄动、气动弹性影响具有鲁棒性,对速度指令和高度指令具有很好的跟踪效果。     

大攻角翼面超声速热颤振分析

分析了大攻角超声速翼面气动加热条件下的热颤振。根据模态叠加法建立颤振运动方程,用当地流活塞及小扰动线性化理论分别计算超声速区和亚声速区的非定常气动力,状态空间法与Runge-Kutta法结合仿真求解结构动响应。给出的两个工程实例计算结果表明,该方法计算精度可满足工程要求。另外还讨论了气动加热效应对颤振的影响。     

基于滑模微分器的吸气式高超声速飞行器鲁棒反演控制

针对吸气式高超声速飞行器气动\推进\结构弹性耦合运动学模型,提出了一种基于滑模微分器的鲁棒反演控制器设计方法。分别采用反演和动态逆方法设计虚拟控制量和实际控制量。引入滑模微分器来精确估计虚拟控制量的导数,避免了传统反演控制方法"微分项膨胀"问题。基于滑模微分器,设计了一种非线性干扰观测器,可对模型不确定项进行精确估计和补偿,从而增强了控制器的鲁棒性。仿真结果表明,该控制器对模型不确定性和气动弹性影响具有较强的鲁棒性,且能实现对速度指令和高度指令很好的跟踪效果。     

统一分区加强解算法解物体穿越激波非定常流场

使用统一分区加强解算法求解非定常可压缩Navier-Stokes方程，模拟了物体分别以超音速和亚声速穿越与之同向运动激波的全过程，并且通过和单区计算结果的详细对比，说明统一分区加强解算法能够简便、有效的处理复杂几何外形和提高网格分辨率。     

用基于混合网格的多重网格法数值模拟绕物体的紊流流动

本文在混合网格基础上用多重网格法求解了紊流N-S方程。混合网格是由在物面附近采用的柱状网格与其他区域则采用的非结构网格组成。方程的求解采用Jamson的有限体积法，紊流模型采用两层Baldwin-Lomax代数紊流模型。用多重网格法来加速解的收敛。数值算例表明，用混合网格及多重网格法来求解紊流Navier-Stokes是非常有效的。     

非结构网格技术应用于固体火箭发动机内流场数值模拟

从二维欧拉方程出发，采用非结构网格技术对固体火箭发动机内流场进行数值模拟，研究了点火初期、燃烧过程中和燃烧结束时固体火箭发动机燃烧室不规则物理区域内形成的复杂流动，分析了流场结构及特性，并对在这一领域应用非结构网格技术存在的问题进行分析。     

用非结构直角网格和欧拉方程计算运载火箭绕流

采用八叉树结构，物面附近用投影法，生成运载火箭绕流计算的非结构直角网格。并将Jameson的有限体积法推广用于这种网格的欧拉方程计算。算例表明，网格生成的时机化费很少，质量好；欧拉方程解算的收敛质量也好。     

基于多块对接网格的多侧喷流瞬态干扰特性研究

基于无限插值理论和多块网格技术,生成了带多姿控发动机的拦截弹结构体网格。利用非定常N-S方程数值求解技术,计算了大气层内高超音速拦截弹带多个横向喷流的瞬态干扰流场,为拦截弹瞬态多侧向喷流的精确控制力建模提供了理论参考依据。研究结果表明,与单喷瞬态干扰相比,瞬态多侧向喷流干扰延迟效应更明显;多侧喷瞬态喷流尾流对尾翼存在更强干扰现象。     

压制性干扰条件下的双基地雷达探测能力研究

从双基地雷达方程出发,提出了双基地雷达探测目标的判定式。在雷达干扰方程的基础上,建立了双基地雷达在自卫式压制性干扰和远距离支援式压制性干扰条件下的干扰方程,并确立了干扰暴露区应满足的条件。根据双基地雷达探测目标判定式和干扰方程,在基于网格剖分算法的基础上,对双基地雷达探测区和干扰暴露区进行了仿真,得出了相应的探测区和暴露区面积。仿真结果表明：基于网格剖分算法的双基地雷达探测能力仿真更贴近实际,有效可行。     

喷管内横流流场的数值模拟

采用显式的原始ＭａｃＣｏｒｍａｃｋ时间分裂格式和改进的ＴＶＤ ＭａｃＣｏｒｍａｃｋ格式分别求解了二维轴对称喷管内横流流场的Ｎ－Ｓ方，得到了与理论分析结果相近似的流场结果图谱，并对以上两方法得到的结果进行了比较，结果表明：ＴＶＤ的方法具有更高的精度。文中还介绍了一种基于ＴＴＭ方法的改进的网格生成方法，用于生成贴壁性和正交性都比较好的网格，实现了完全由边界网点分布来确定内部网格点的分布。     

潜入式摆动喷管两相内流场数值模拟

气相采用三维薄层近似N-S方程，固相采用拉格朗日坐标系下颗粒轨道模型，计算方法采用显式ENO差分格式和分区多块网格技术，建立了固体火箭发动机内流场模拟软件。对轴对称JPL喷管进行模拟，计算结果与文献进行比较，符合较好。推广到实际三维装药型面时，发现颗粒数目增加导致计算效率降低。采用软件中的气相计算模块对潜入式摆动喷管内流场进行数值模拟，比较了喷管摆动角在发动机喉道附近引起流动变化，研究了非对称流动对喷管内表面压强的影响。     

基于连通检测和区域扩展的多重聚类方法研究

针对现有聚类算法在分选过程中存在辐射源定位和脉冲划分的困难，提出了基于连通检测和区域扩展的多重聚类方法。使用网格统计的方法对脉冲列分布密度进行量化，再借鉴图像分割算法，将量化的网格矩阵划分连通区域。基于峰值检测和梯度下降的原则扩展连通区域，通过边缘检测得到最终聚类结果。仿真结果表明，该算法能保证更多的聚类脉冲列参与参数分选，提高辐射源识别率。