一种基于欧拉方程的动导数简化计算方法

动导数是飞机飞行品质分析的关键参数,在飞机概念设计阶段快速可靠地预测动导数可缩短设计周期并减少成本。基于欧拉方程提出一种飞机组合动导数的简化计算方法,在模拟小幅强迫振动过程中忽略网格变形,通过对物面处流体添加法向扰动速度项模拟物面运动对流体造成的扰动。采用欧洲动导数研究标准模型DLR-F12对该方法进行了验证计算,计算结果与欧拉方法计算结果及试验数据具有较好的一致性。     

基于CFD数值仿真技术的飞行器动导数计算

为了研究表征飞行器动态稳定特性的动导数,结合计算流体力学方法,基于滑移网格技术,采用小幅度强迫振动法和差分法两种非定常方法,推导了俯仰组合动导数以及滚转阻尼导数的计算表达式,并结合多参考系模型采用准定常方法计算滚转阻尼导数。建立了基于计算流体力学技术的飞行器动导数计算方法。对国际动导数标模Finner导弹进行了计算验证。两种非定常方法对俯仰组合动导数的计算误差分别为0.65%、6.13%,对滚转阻尼导数的计算误差分别为2.5%、0.06%。求解滚转阻尼的准定常方法对于滚转阻尼导数的计算误差为2.67%。三种方法的计算结果与风洞试验结果吻合的很好。非定常方法计算精度高,准定常方法可以快捷迅速得到结果。本文的数值方法不仅适用于超声速,同样可以用于亚、跨声速时的动导数计算。     

某型导弹动导数风洞试验和数值模拟

动导数的精确预测对于弹箭类武器的成功研制至关重要。以某型导弹动态特性的研究为对象,在FD06风洞中采用自由振动法完成动态特性试验,获取该型导弹俯仰和滚转两个方向的典型振幅衰减曲线和动导数结果;针对该型导弹,基于双时间推进法数值求解非定常N-S方程和强迫运动方程,获取该型导弹俯仰动导数和滚转动导数。结果表明:CFD计算结果与风洞试验结果有较好的一致性;在给定马赫数和攻角下CFD数值模拟结果可以为后续工程型号研制提供参考。     

基于刚性动网格技术的动导数数值模拟

基于刚性动网格技术,建立了俯仰动导数的非定常数值计算方法.首先使用小幅度强迫俯仰振荡方法求解俯仰组合动导数,然后利用小幅度强迫升沉振荡方法求解洗流时差导数,通过两者相减即可得到俯仰阻尼导数.利用国际动导数标准模型Finner导弹进行算例验证,计算得到的俯仰组合动导数与试验值误差为2.76%,洗流时差导数值约为俯仰阻尼导数的11.5%,与文献的结果一致.结论表明:动导数单独模拟方法具有较好的工程实用价值,且可以推广到横向以及航向的动导数数值模拟.      

基于CFD的动导数计算与减缩频率影响分析

以CFD流场与气动力计算为基础,分别利用小幅度强迫振动模式和差分两种方法推导了飞行器的纵向组合动导数计算表达式;建立了基于CFD技术的飞行器动导数计算方法;以国际标模Finner导弹为例,利用所推导的两种方法进行该导弹的动导数计算,重点分析了减缩频率对动导数计算的影响,提出了利用CFD方法开展动导数计算时减缩频率的选择...     

基于结构动网格的无人机地面效应研究

为了更好地分析地面效应对无人机气动特性的影响,采用一种新的结构动网格生成技术,运用CFD手段对某无人机的地面效应进行了研究。在一套静态网格CFD计算基础上,运用结构动网格生成技术,获得非定常的网格移动,得到一系列的不同攻角、地面效应高度的CFD网格及初始流场。在此基础上进行定常流场求解,计算得到地面效应的定常CFD结果。最后对某双尾撑布局无人机的地面效应进行了CFD计算研究,对相应的全机气动特性的影响进行了分析。计算结果表明：地面效应对平尾的影响大于对机翼的影响并使平尾的失速迎角提前;该方法是可行的,具有一定的工程应用价值。     

飞艇动导数与附加质量相互融合的方法

飞艇的动导数和附加质量分别表征其在有黏流和无旋无环流中所受非定常气动力/力矩。为了解决在飞艇建模中这两种气动系数的融合问题，分别研究了动导数和附加质量的成分划分问题以及同成分气动系数的融合方法。通过介绍能同时兼容有黏流和无旋无环流的气动力和力矩分析理论，得出同成分气动力/力矩融合时应当取有黏流中的结果并摒弃对应的无旋无环流结果。通过研究气动力/力矩与运动体当前运动参数的关联性，建立依据当前运动参数划分气动力/力矩或气动系数的方法。为了使两种流场中的气动系数分类方法相同，对飞艇的当前运动参数进行重构，使得气动系数在两种流场中均可按重构后运动参数明确划分和计算。根据研究结果，建立了一种动导数与附加质量的新融合方法，并讨论了它与现有文献方法的差异。通过算例分析不同融合方法对飞艇纵向扰动运动特性的影响，说明采用新融合方法的必要性。     

基于阶跃响应的飞行器单独静、动导数数值求解方法

基于阶跃响应方法,结合刚性动网格技术,对飞行器的单独静、动导数的精细化数值计算进行了研究.以纵向为例,通过给物面施加恒定附加攻角,求解得到阶跃响应运动过程的非定常气动力,求导得到静导数.同样给物面施加恒定的俯仰角速度,并同时强迫物面平动以抵消俯仰转动产生的附加攻角影响,可由非定常气动力求导得到动导数值.分别利用NACA0012翼型和三维SACCON飞翼无人机进行了计算验证,各攻角下的静、动导数值与文献、试验结果吻合得很好,最大误差不超过5%.结论表明:基于阶跃响应的单独静、动导数直接模拟方法计算耗时仅为传统强迫振动方法的21%,效率相对较高,且可推广到横航向的动导数计算,为飞行器的稳定性研究提供参考.      

基于准定常的飞行器动导数的高效计算方法

结合工程设计经验和工程计算方法,针对飞行器动导数计算方法展开研究。基于求解准定常欧拉方程发展了一种适合于亚跨超声速状态下,任意气动外形气动导数的高效计算方法,并将计算结果与一般非定常方法的求解结果以及国际标准实验数据进行了对比,研究结果表明该方法结果可靠、效率较高。     

升力浮力复合型飞艇动导数分析

计算了球体和椭球体的平移、转动加速度导数,数值模拟了NACA0015翼型的俯仰振荡,验证了计算流体动力学(CFD)方法计算动导数的有效性.模拟升力浮力复合型飞艇的强迫振荡运动,得到了全艇以及机翼、艇身、平尾和垂尾的动导数,并与常规飞艇进行了比较分析.机翼对速度导数的贡献导致复合型飞艇升沉速度导数、滚转角速度导数较常规飞艇显著增大.常规飞艇加速度导数以惯性力成分为主,但黏性力成分以及机翼对艇身和尾翼的干扰对复合型飞艇升沉、滚转和俯仰3个方向的加速度导数有重要影响,不能将加速度导数简化为附加质量系数.      

稳定性导数计算方法研究

稳定性导数的准确计算对于飞行器的操稳特性具有重要意义.应用计算流体力学(CFD)方法中的非结构化动网格技术建立能够模拟飞行器做周期性俯仰运动的强迫振荡法,以国际动导数标模Finner导弹为验证算例,获得不同马赫数下俯仰力矩系数的迟滞环曲线,进而计算Finner导弹在不同马赫数下的静稳定性导数和动稳定性导数.结果表明:本文计算得到的俯仰静稳定性导数与试验数据非常接近;在亚音速和超音速范围内,动稳定性导数计算结果与试验值很接近,但在跨音速范围内,本文计算结果与试验曲线的规律性一致,但误差较大.     

超/高超声速飞行器动态稳定性导数极快速预测方法

飞行器设计早期阶段需要预测大量工况下的动导数。本文发展了一种面向超/高超声速飞行器的动导数极快速预测方法：首先基于当地流活塞理论，将飞行器进行小幅非定常运动所受到的气动力分为受自由来流引起的无附加扰动项以及受物面变形或运动引起的附加扰动项；通过当地表面斜度法、激波后等熵关系求解物面当地流动参数，进而结合非定常运动规律求出飞行器所受非定常气动力；再采用待定系数法对非定常气动力进行提取、辨识，最终得到超/高超声速飞行器动导数。该方法克服了传统方法对CFD流场参数的依赖和耦合，具有极高的计算效率；同时典型算例验证表明，该方法在超声速、高超声速工况下都能够很好预测动导数变化趋势。将该方法应用于复杂外形飞行器动导数预测，并讨论了与CFD方法的误差来源。本文方法可作为高速飞行器总体设计阶段布局选型的工具。     

飞艇续航时间论证方法研究

由于飞艇续航时间受飞行速度、重力浮力配平方式、载油量变化、风场条件等多种因素的影响,飞机续航时间计算方法并不适用于飞艇,因此提出了一种新的基于试验设计和统计分析方法的飞艇续航时间论证方法。分析了工作高度的风场速度分布规律,分别给出完全靠发动机推力矢量平衡和完全靠动升力平衡两种不同情况下飞艇续航时间计算方法,以某大型对流层飞艇为研究对象,计算了初始方案在不同平均巡航速度、配平重量和配平方式下的平均燃油消耗率和最大续航时间,对经过减阻、减重、减少耗油率等技术改进后优化方案能达到的续航时间和概率进行了分析和论证。研究表明,某飞艇初始方案在使用区域风场条件下基本能满足留空时间72h的指标要求,经技术改进后,指标仍有较大的提升的空间。方法考虑了飞艇续航时间的多种影响因素、取值变化和交互影响,比较适宜在飞艇这类涉及重浮力配平和任务耗油率变化的浮空器上使用。     

飞行器动态稳定性参数计算方法研究进展

动态稳定性参数（简称动导数）是飞行器控制系统设计、飞行器动不稳定发生边界分析及相应动态稳定性判据研究的关键气动参数。在对飞行稳定性问题进行概述的基础上,介绍飞行器动态稳定性参数数值模拟的国内外研究进展。并按照理论方法、工程近似方法及计算流体力学（CFD）模拟方法的动导数发展方向对近年来主要的动导数计算方法进行了综述,评价了各种动导数预测方法的优缺点,指出了动导数数值模拟在理论基础、非定常气动力建模、预测方法精度和效率等方面存在的问题。最后对动导数数值模拟的发展趋势进行了展望。     

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

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

飞机平尾偏转对大迎角动态气动特性的影响

给出了 BJ-1背景机有不同平尾偏角时静态和大振幅俯仰振荡情况下的气动特性。根据风洞实验结果,采用基于三角剖分数据库建模的方法,计算了任意平尾偏角和任意振动频率时的背景机动态气动特性,与实验结果比较表明,提出的有平尾偏角的建模方法能提供任意舵偏时的全机气动特性。在非定常建模的基础上,采用动导数仿真,分析了不同平尾偏角时的动导数,结果表明,水平尾翼偏角主要对零动态阻尼所对应迎角有较大影响,随着舵偏角的增大,该迎角减小。     

非结构混合网格自适应并行技术

计算流体力学（CFD）模拟实际工程问题所采用的网格规模可达千万量级，并行技术是减少计算时间的有效方法。耦合流场信息的网格自适应技术能有效动态优化计算网格，被NASA视为一项亟待发展的CFD关键技术。混合网格自适应系统包含网格分布优化、表面网格投影和空间网格匹配等关键技术。针对以上3项关键技术分别建立了高效的并行算法。首先，提出了"先唯一后同一"的两步法策略实现了网格单元分布优化过程的并行相容性；其次，基于局部曲面拟合思想，实现了曲面重构和新增物理网格点投影的完全并行；再次，提出了空间网格匹配技术的半并行算法，快速解决了网格单元交错问题。为了提高后续流场计算的并行效率，发展了基于并行重分区-网格数据迁移方法的动态负载平衡技术，并采用圆柱激波流场自适应模拟对动态负载平衡技术进行初步验证。最后，采用三角翼自适应加密测试了自适应系统的并行效率。结果表明，建立的混合网格自适应系统并行效率较高，且相比流场求解耗费总时间的比例低于1%。     

CFD法中的动态失速模拟

本文介绍了CFD动态失速模拟所需要捕捉的关键物理,CFD法一般情况和较先进的CFD法,最后重点介绍了CFD法动态失速模拟的当前水平.总的来说在动态失速预测方面,CFD法优于综合法,CFD/CSD耦合结果与试验结果之间的吻合度还是可以接受的.俯仰力矩曲线的相位、大小和形状与试验数据吻合得较好,耦合解与试验数据之间的主要偏差在于法向力预测上.