高效非结构网格变形与流场插值方法

非结构动网格在含动边界的非定常流动模拟中应用十分广泛。为了提高网格变形的计算效率和对复杂外形的适用性,发展了高效的K近邻-径向基函数（KNN-RBF）动网格方法。并采用高效的交界面节点一致性方法实现了网格变形的并行化。由于边界的运动和变形形式多种多样,在自主开发的计算流体力学（CFD）程序中实现了动网格的用户自定义函数（UDF）功能模块,可以灵活地设置此类问题中的边界运动。在模拟大变形问题时,采用了全局网格重构的策略对网格质量较差的区域进行处理。提出了快速查找新旧网格最近点的两级KD（K-Dimensional）树方法,并根据最近点在新旧网格间进行流场插值。从而建立了高效健壮的并行流场插值方法。通过对若干算例进行模拟,验证了所发展方法的正确性和健壮性。     

动基座近舰面流场数值模拟

航母尾迹流场对舰载机的着舰有较大影响,所以需要对其流场特点进行研究,分析不同状态下舰载机气动特性的变化。采用嵌套网格技术对航母处于垂荡状态下无人机的着舰进行了模拟。首先,利用SFS2舰船进行数值计算,验证了舰船流场的数值模拟方法。然后,对比了单独无人机定常与非定常计算结果,表明所建立的嵌套网格适用于无人机流场的模拟。接着,对航母单相流和两相流的流场结果进行了分析,结果显示甲板下方的流动对甲板上方流场没有大的影响。因此,忽略了水的影响只对航母在空气流场中的特性进行研究,结果表明航母尾迹非定常特性明显,静止航母下滑轨迹上的速度均处于周期性波动状态,且波动幅值随着远离航母而逐渐衰减;而在垂荡情况下航母尾迹变得更加紊乱,水平方向速度波动的周期性减弱,但垂向速度的波动幅值进一步增大。对于静止航母,无人机在不同时刻着舰气动特性的变化也存在差异;当航母处于垂荡状态时,无人机的升力和俯仰力矩在短时间内会有更大的波动。     

转子叶栅非同步振荡发声特性研究

轴流压气机转子叶片排振荡疲劳失效是常见的气动弹性失稳问题。当转子叶栅处于非同步振荡状态下时,压气机管道内部将伴随着异常噪声的产生,这类噪声的频率既不是分布在叶片通过频率及其谐波上,也不是分布在转子轴频率及其谐波上,同时也不满足简单多普勒效应。为了解释这种异常噪声现象,以三维升力面理论为基础,讨论叶片对其附近流体施加非定常载荷的发声问题,给出了转子叶栅振荡异常发声问题的物理解释,建立了声场频率及模态特性与转子声源特性的直接关系式,并给出了频率特性不同于叶片通过频率且不符合简单多普勒效应的完整解释。在此基础上,通过机理性实验研究证实了该模型的正确性。实验结果表明,理论预测声场的频率和周向模态特性与实验结果完全一致。     

控制力矩陀螺磁悬浮转子的指数趋近积分滑模控制研究

针对微重力、动框架等多源扰动下磁悬浮控制力矩陀螺磁轴承转子系统的强非线性、参数摄动及未建模动态等问题,开展控制力矩陀螺的磁轴承转子鲁棒控制研究。建立了磁轴承转子轴向通道动力学模型与电磁力工作点线性化模型,分析了系统参数摄动及未建模动态,并在此基础上提出了一种非线性系统指数趋近积分滑模控制方法,设计了含有电流及位移的积分滑模平面,采用指数趋近律与饱和函数抑制抖振实现到达阶段滑模控制律。实验结果表明：采用该方法可有效抑制滑模的抖振现象,参数摄动时系统具有较佳动态性能,同时在转子转动、框架工作时保证系统良好鲁棒性。     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

悬停状态下小型共轴直升机操纵响应特性分析

从在研共轴双旋翼直升机的工程实际出发,建立了悬停状态下旋翼非定常气动特性的计算模型.引入Leishman-Beddoes指数函数的半经验公式,建立了二维翼型非定常气动模型;引入干扰因子到动态入流模型,建立了反映共轴双旋翼直升机上下旋翼气动干扰的诱导速度模型;从桨叶的挥舞动力学模型出发,利用四阶Runge-Kutta算法求解桨叶刚性挥舞角的数值解.通过计算分析,得到了悬停状态,总距突增时上下旋翼升力和扭矩的动态响应特性,并分析比较了半差动和全差动航向控制方式的操纵响应差别,为共轴式直升机机动特性的研究作了必要的准备.      

基于动态逆的空空导弹奇异摄动中制导律

结合奇异摄动和动态逆方法,以飞行时间和末端能量的线性组合为性能指标,推导了一种近最优中制导律.导弹状态方程根据状态变量时间常数写成三时间尺度系统.降阶系统包括纵程、横程和比能动力学方程.用最优控制理论求解出降阶系统最优解,即原系统的外解.动态逆应用于边界层修正,得到一个闭环控制器,其性能可根据爬升、高空等高飞行和下滑阶段或其他情况的不同要求进行调整.仿真结果表明:该制导律满足了对弹道平滑度、飞行时间和末端能量的要求.中段末端航向误差小,为中、末制导的平稳交接班提供了保障.      

计入畸变修正的旋翼尾迹前飞状态稳定性分析

现有的旋翼尾迹稳定性分析方法在涡线过近处会得出不现实的过大发散率结果,因而存在一定的局限性。为了消除这种前飞状态下尾迹高度畸变对其稳定性分析所造成的不合理干扰,提出了一个旋翼尾迹在正弦扰动下,结合涡线畸变修正的线性化稳定性分析方法。在该方法中,将尾迹涡线离散为一系列涡元,并对其端点施加沿尾迹涡线分布的正弦扰动;考虑了桨尖涡的自诱导、互诱导以及与桨叶的干扰作用,并创新地引入在伪涡核对稳定性分析中,因尾迹高度畸变产生的不合理干扰而进行的修正。以H-34型直升机旋翼模型为例,在前飞状态下,重点计算分析了尾迹向后飘移、涡核大小和前飞速度等对旋翼尾迹稳定性的影响。研究结果表明:引入伪涡核模型可以有效地修正在涡线过近处所产生的不现实稳定性分析结果;前飞速度、涡核大小等因素会在一定条件下有限地影响尾迹的不稳定性,但都不改变前飞状态旋翼尾迹的本质不稳定性。     

高超飞行器内流道激波振荡问题的数值研究及试验验证

采用数值模拟和风洞试验的方法,对吸气式高超声速飞行器盲腔状态的流动特性进行了研究.采用“双时间步”方法进行了内外流一体化的非定常数值模拟,利用彩色纹影系统对高超声速飞行器前体流场进行显示,并采用动态压力传感器测得了飞行器内流道的壁面压力随时间的变化.结果表明:盲腔构型在高超声速飞行中会出现周期性的激波振荡现象.数值模拟所得流场变化特征、内流道壁面压力振荡周期和壁面压力变化趋势与试验结果吻合良好.     

双三角翼俯仰振荡运动的流场特性数值模拟

 建立了适用于双三角翼大迎角非定常分离流场模拟的数值方法,研究双三角翼俯仰振荡时的动态流场特性,给出动态流场结构和气动力性能随迎角的变化规律,重点考察了减缩频率、转轴位置、平均迎角和振幅等参数对动态流场迟滞效应和气动力曲线迟滞环的影响。研究结果表明:俯仰振荡到相同大迎角时上仰和下俯的流场存在明显差异;减缩频率对气动力迟滞效应的影响相对大于转轴位置;平均迎角的变化导致双三角翼背风区流场结构呈现不同流态,而振幅的大小决定这些流态的数目,事实上俯仰运动时如果跨越的流态数目越多则流场结构的动态响应滞后现象就越显著。通过数值分析,有利于提高对双三角翼在俯仰振荡运动条件下的非定常特性和流场滞后效应等非线性现象的认识。