某改型飞机时差动导数对尾旋时间历程估算的影响

提供了利用欧拉方法计算的某改型飞机尾旋时间历程结果，文中比较了使用和没使用洗流时差动导数数据对计算结果的影响，说明洗流时差动导数数据在这种模拟计算中是非常重要的，它使计算结果更合理。     

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

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

某型飞机时差导数对预测尾旋特性的影响

本文提供了利用欧拉方法计算的某型飞机尾旋时间历程结果，文中比较了是否使用洗流时差动导数数据对计算结果的影响，说明洗流时差动导数数据在这种模拟计算中是非常重要的，它使计算结果更合理。     

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

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

基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

绕低频振动旋成体的跨音速非定常流核函数方法

本文对无粘跨音速流中作俯仰振动的细长旋成体,给出了一种计算表面非定常压力分布和动稳定性导数的近似方法,即用时间线化近似,在引入局部线化假设后,将跨音速小扰动方程化为一个积分方程,并用精确物面边界条件求数值解。本方法算例与其它理论和实验结果作了比较,表明用它可以得出较合理的结果。     

机动飞行数值模拟关键技术及初步验证

建立飞机机动飞行仿真及飞行性能评估能力是未来航空CFD发展的重要目标之一,为此,对机动飞行数值模拟的方法、原理、关键技术及潜在的应用方向进行了探讨。采用惯性系下任意拉格朗日-欧拉(ALE)形式描述的Navier-Stokes方程与机体运动学方程相耦合的方法,结合动网格技术建立了通用的机动飞行数值模拟应用框架。从微分方程形式的运动学方程出发,建立了与流动控制方程进行紧耦合和松耦合计算的两套计算方案,并提出了利用预报-校正技术提高耦合计算精度的方法。计算结果表明,所提方案合理、可行,为单自由度/多自由度机动飞行模拟、动导数分析等计算提供了可靠的新工具。     

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

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

准静弹性飞机动态特性的计算研究

本文在小扰动和准静弹性假设下对弹性飞机纵向动态特性进行了计算和分析。对弹性变形对动态特性的影响进行了初步讨论。本文采用格林函数法求得飞机的气动力,用正则模态法计入弹性变形影响并引入准静弹性假设,得到了研究纵向动态特性所需的全部刚体飞机气动导数和弹性飞机气动导数。利用这些导数,对准静弹性飞机动态特性进行计算研究。并将计算结果与刚化模型的同一飞机进行比较和分析。     

滚转和锥进运动对弹箭动导数求解的影响

现有的弹箭动稳定性导数计算通常忽略了耦合运动效应的影响,而弹箭在实际飞行中,通常同时存在俯仰、滚转以及锥进运动。为此,本文提出了基于欧拉转动定理和滑移网格技术的复杂角运动模拟方法,该方法利用著名的罗德里格斯转换矩阵插值求得弹箭在每个时间步的角速度修正值,并指定给球形滑移网格区,从而模拟弹箭俯仰滚转锥进相耦合的三自由度耦合角运动,在此基础上,通过对非定常气动参数进行求解辨识,分析了滚转频率和锥进频率对其俯仰组合动导数和升力系数动导数的影响规律。结果表明:本文提出的复杂角运动模拟方法,可有效消除姿态角计算的累积误差,实现对弹箭任意给定角运动的准确模拟;在耦合角运动状态下,弹箭气动力系数迟滞环发生明显的振荡和偏移,俯仰组合动导数和升力系数动导数均随滚转频率和锥进频率增加而显著变化,且变化规律具有非线性特性。     

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

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

飞机/外挂物组合体绕流的数值计算

采用嵌合体埋入方法计算带外挂翼-身组合体复杂外形的流动,将飞机/外挂物外形的流动区域分为机翼-机身和外挂物两个子区域,在机翼-机身子域上求解欧拉方程,在外挂物子域上求解全速势方程,通过两个流场的耦合迭代计算全流场的解。对C-5A飞机及某战斗机的计算表明,数值模拟与实验吻合。实现了不同子域流场使用不同方程的计算方法。     

飞机机翼表面霜冰的三维数值模拟

基于欧拉两相流理论对三维情况下飞机机翼表面的霜冰进行了数值模拟.根据水滴拟流体模型建立三维水滴控制方程;提出一套水滴控制方程的数值求解方法;由三维水滴流场的求解结果计算机翼表面的水滴收集特性,提出一种三维积冰外形的生成方法,完成了对飞机机翼表面霜冰的三维数值模拟.对ONERA M6机翼在不同迎角下霜冰的积冰情况进行了数值预测,并分析了结冰条件对积冰的影响.      

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

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

折叠机翼变体飞机纵向操纵性与稳定性研究

针对一种飞翼布局折叠机翼变体飞机方案,建立了相应的研究模型,使用涡格估算方法计算出机翼折叠角度对全机纵向静稳定性的影响.利用飞机纵向小扰动运动方程,分别得到了机翼展开和折叠状态的长、短周期模态,并对其动稳定性进行了讨论.结合工程估算方法和风洞试验方法,计算出内段机翼折叠过程全机力矩系数、升降舵操纵导数以及升降舵配平偏角...     

F12全机动态特性数值模拟

采用非定常正弦振荡和准定常旋转坐标系两种数值模拟方法,对F12全机模型三个攻角进行动态特性数值模拟,得到飞行器的纵向组合动导数和阻尼动导数,计算结果与试验数据对比具有良好的一致性,同时用无粘和粘性两种计算方程进行准定常旋转坐标系模拟,计算结果一致,表明两种数值模拟方法准确性高,计算速度快,可为设计提供快速可信的数据支持。     

超音速、高超音速有攻角的三维薄翼非定常二次理论及其应用

 动导数对于飞行器的动态品质分析和控制系统的设计是必不可少的参数,它直接影响飞行器的动态特性。因此,国外自七十年代以来,无论在风洞实验、飞行试验,还是分析、计算等方面都作了很多工作,取得了较大的进展。 动导数计算说到底主要是非定常气动力计算问题,已有的线化理论不能满足需要,对考虑非线性效应的动导数计算方法的需求则日显迫切。二次理论是计算非线性效应的一种方法,它是由第二次叠代或物面偏角参数的二次展开法构成,主要非线性效应由二次方程算出。     

稳定性导数计算方法研究

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

外挂物投放技术的数值模拟

应用Icemcfd和Cart3d软件对飞机投放外挂物的流场进行了数值模拟：对欧拉方程用有限体积法进行离散，运用自动及人工干预网格细化、多重网格技术的加速收敛、自动位置和网格更新的六自由度投掷分析，得到外挂物的分离轨迹。计算结果与AEDC的试验结果作对比，两者吻合较好。     

计算高亚声速情况下三维机翼动稳定性导数的面元法

本文应用我们提出的跨-亚音速三维谐振机翼的局部线化面元法,同时采用了与之配偶的亚临界定常流绕三维厚度机翼的面元法迭代解,对几种不同平面形状的机翼,进行了俯仰和滚转阻尼导数的计算和分析,并和现有的实验数据以及线化理论的计算结果进行了比较,从而表明,本方法可以得到满意的动导数计算结果;还表明,当接近临界马赫数时,动导数计算必须计及跨音速非线性效应。