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

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

飞行器大迎角下俯仰静、动导数的数值计算

 采用一种适用于非定常流动计算的逐次超松弛（SOR）时间推进方法耦合求解非定常Navier-Stokes方程和强迫运动方程,数值模拟了NACA0012翼型、弹道外形(HBS)和有翼导弹等飞行器在不同初始迎角下做强迫俯仰振荡的黏性动态流场,给出动态气动力及力矩系数的时间历程;采用积分法获得飞行器的静导数和动导数,将计算结果与实验结果及半经验理论预测方法进行比较并分析飞行器在大迎角下动态特性的非线性特征;计算结果与实验值相吻合,表明本文采用的非定常数值模拟方法在大迎角下能够较准确地模拟飞行器的动态特性,具有工程应用价值。     

扑翼飞行器平面形状对气动力的影响研究

为了了解扑翼飞行器在扑动过程中的气动力,采用非定常涡格法模拟扑动过程中的气动力计算,得到了刚性扑翼在扑动周期内气动力的变化及尾涡的形态。在此基础上,研究了扑翼的周期平均气动力随扑动频率、幅度、迎角及来流速度的变化关系,并进一步研究了不同的机翼平面形状对气动力的影响。计算结果表明,迎角和扑动频率的增大能够增加扑翼的升力,不同翼面形状的扑翼在升力特性及推力特性方面具有不同的优势。     

基于非定常CFD方法跨声速民用飞机标模俯仰气动导数计算研究

采用非定常计算流体力学(CFD)方法和动网格技术,计算了跨声速民用飞机全机构型TCR气动标准模型在0°~30°迎角范围内气动力系数的时均值、俯仰静导数和俯仰动导数;通过与试验结果对比,研究了非定常CFD方法在较大迎角范围内对包含鸭翼、垂尾等部件的复杂外形气动导数的计算精度。与试验结果对比表明,在16°~24°迎角范围内,受鸭翼和机翼上流动分离再附以及内外翼分离涡相互作用等复杂流动结构的影响,俯仰力矩系数的动导数与试验数据有相对较大的差异。     

有地面干扰的动导数计算方法

提出了一种有地面干扰情况的飞行器动导数计算方法,其基本思路是从非定常升力面气动力理论出发,利用飞行器的低频运动的气动力,通过振动频率的级数展开,选用适当的气动力部分进行积分,得到所需的气动力导数,为阿了面的影响,将传统的镜像法结合到非定常用功力理论理论方法中,计算结果表明,地面干扰随距地面的高度变化而迅速变化,而且某些动导数的符号都发生了变化,说明地面干扰城飞行动力学的分析是应该加以考虑的。     

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

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

环量控制翼型气动力迟滞效应研究

为研究环量控制翼型在非定常运动中的气动力迟滞效应,通过求解环量控制翼型的强迫俯仰振动数值模拟数据,研究了不同减缩频率、振幅和射流动量系数对环量控制翼型气动特性的影响规律,分析了不同动量系数对非定常流场变化的作用机理。研究结果表明:减缩频率和迎角振幅越大,气动力迟滞效应越明显;小迎角下动量系数越大,迟滞效应越微弱,而大迎角下动量系数越大,迟滞效应越明显,并表现出提前动态失速的特征;气动力迟滞效应主要受翼型上表面压力分布影响,与后缘涡和上表面分离涡关系密切。所得结论对于环量控制技术在无人飞行器控制和操纵方面的应用具有一定的实际意义。     

飞机大迎角滚转运动频率域建模与稳定性分析

本文运用频率域建模法,分别建立了飞机大迎角滚转运动的Ｆｏｕｒｉｅｒ变换模型及非线性代数模型,并从数学模型中提取动导数。基于非线性非定常数学模型与局部线化的动导数数学模型,计算了飞机滚转运动的稳定性,并运用李雅普诺夫稳定性分析理论,进行飞机滚转运动稳定性分析。结果表明,大振幅非定常滚转运动实验可以用于飞机绕体轴滚转运动的稳定性分析。从大振幅实验数据提取的动导数,基本反映了飞机的非定常气动力特性。通过     

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

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

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

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

平头增阻再入体俯仰动态特性计算与流动机理分析

采用自由振荡法数值模拟了一种平头增阻再入体的高超声速俯仰动态特性,沿用稳定性导数模型作为动态特性的指标,通过把稳定性导数视为攻角的函数来考虑非线性影响,并应用奇异分解线性最小二乘法辨识稳定性导数。对于本文研究对象,计算结果表明,平头外形在较高马赫数和小攻角范围内存在动不稳定现象,随着再入时马赫数的降低,动不稳定的攻角范围不断缩小,直至演化为动稳定的。文中对平头外形出现这种动态特性现象的空气动力学原因进行了分析,有必要作进一步的确证与研究。     

超,高超声速机翼俯仰导数的欧拉方程摄动解法

应用片条理论，将文献（１）的尖前缘任意翼型的俯仰稳定性导数的计算方法作、推广、用于计算任意平面形状的三维机翼。在各种特殊的情况下，其计算结果与现有的理论都符合得很好。在激波附体的条件下，本文的计算方法可用于计算任意迎角。任意平面形状机翼的超、高速声速俯仰稳定性导数。     

迎角动态变化对二元高超声速进气道气动特性的影响

针对高超声速飞行器受到扰动后迎角可能会发生瞬时大幅度改变或振荡的问题,对来流马赫数为4.03、迎角动态变化的二元高超声速进气道流场进行了非定常数值模拟,分别在进气道等速上仰和按正弦波振荡条件下研究了迎角动态变化对进气道气动特性的影响。计算结果表明:迎角动态变化时,进气道流场特征和性能参数的变化规律和稳态时情况基本保持一致,但是存在明显的迟滞现象;二元高超声速进气道动态上仰的迎角速率越大,进气道发生不起动时的迎角值越大;在进气道下壁面附近的低速区,非定常效应的影响显著,受进气道固壁面运动的影响,低速区的产生、发展及消失影响进气道气动特性。     

小迎角高超声速非旋转钝锥非线性运动分析

采用分段拟合技术对小迎角高超声速(马赫数Ma=6)非旋转钝锥双平面拍摄风洞自由飞试验结果进行了分析,获得了非旋转钝锥在小迎角高超声速下气动导数非线性的具体形式,分析了非旋转钝锥在小迎角下的非线性运动特点。研究发现在小迎角范围内钝锥的动导数系数均呈现明显的非线性,而静导数系数非线性较弱,可近似为线性。各组试验的静、动导数系数的非线性形式表明,除模型Ⅱ外的其余模型均为两方向振幅不同的极限圆锥运动,模型Ⅱ为极限平面运动。每组试验两个平面的静导数系数在小迎角范围内保持基本相等。而无论出现极限圆锥运动的四组试验还是出现极限平面运动的一组试验,模型在俯仰和偏航两个方向的动导数系数均存在较大差异。     

超声速和高超声速机翼俯仰导数—当地流活塞理论解法

本文给出计算超声速和高超声速尖前缘三维机翼俯仰导数的一个近似解析方法。该方法是根据当地流活塞理论和片条理论导得的。在激波附着于机翼前缘情况下,本方法可用于计算任意迎角,颇为一般的平面形状、尖前缘剖面机翼的俯仰导数。中小平均迎角的数值结果与有限的试验数据和其它理论作了比较,示例表明本方法的结果令人满意。     

高超声速飞行器的滑模观测器设计

对一个高超声速飞行器设计了滑动模态观测器来估计攻角和航迹倾角。通过实时的求解一个极点配置问题,使观测器的切换增益能够根据观测器状态的变化而改变,从而保证了观测误差的动态特性在整个飞行包线内的一致性,给出了单输出和多输出非线性系统的滑模观测器设计过程。设计和仿真算例验证了这种方法的有效性。     

超声速高超声速球锥绕流的边界层稳定性特点初探

本文对一个球锥组合体,在超声速来流及零迎角条件下,对其边界层的稳定性特点进行了分析并讨论了转捩预测问题。结果表明转捩总是最先出现在来流马赫数不太大,或来流雷诺数较高并且具有较大的周向速度或周向波数的情况下。     

逃逸飞行器模型高超声速风洞自由翻滚动导数试验技术

高超声速风洞模型自由翻滚动导数试验技术是为满足航空航天飞行器 0°～ 360°全攻角范围内的动导数测量及产生极限环振动现象的研究之急需而研制。简要介绍了试验技术研制的自由翻滚试验装置、液体轴承、角度测试系统与系统建模的大攻角非线性数据处理技术。该项技术已成功地为逃逸飞行器模型提供了满意的试验数据     

冲压发动机进气道不起动边界分析

不起动边界是超声速进气道研究的重要内容,它是进气道保护控制的基础和前提。针对这一问题,对冲压发动机进气道进行了不同边界条件下的数值模拟,对稳态流场结果进行了分析。基于仿真数据,利用量纲分析工具对进气道最大抗反压能力进行了分析,讨论并给出了两进气道流动的相似条件。分析结果表明:在不考虑攻角变化的条件下,进气道不起动边界主要与来流马赫数有关。     

高超声速滑翔飞行器摆动式机动突防弹道设计

 为提高大升阻比高超声速滑翔飞行器机动突防能力,提出了一种侧向摆动式机动弹道的设计方法。基于动态逆的思想,建立了侧向摆动式机动弹道的弹道形式和所需倾侧角间的关系模型;在平衡滑翔假设下,将包括动压、过载和热流在内的飞行约束转化为迎角约束,从而确定了迎角-速度飞行走廊;在此基础上,设计了平衡滑翔情况下满足侧向摆动式机动及飞行约束所需的迎角变化规律。根据所设计的迎角和倾侧角,即可实现平衡滑翔情况下预定机动模式的侧向机动。以HTV(Hypersonic Technology Vehicle)为例进行仿真分析,结果表明该方法能够快速设计出预定机动幅度和机动频率的侧向摆动式机动突防弹道。