边条翼和近距鸭翼布局模型动态气动特性分析

针对边条翼与近距鸭翼这两类典型战斗机布局模型,在中国空气动力研究与发展中心FL-24风洞进行了大振幅俯仰动态试验与模型自由摇滚试验,并对比分析了边条翼与近距鸭翼布局模型高速大迎角的动态气动特性。结果表明:边条翼模型纵向动态特性明显优于近距鸭翼模型,尤其是俯仰力矩迟滞效应更强;近距鸭翼模型在攻角26°~45°区间出现了较大的滚转力矩,容易诱发摇滚运动;最后,通过自由摇滚试验验证了俯仰动态试验分析结论,即近距鸭翼模型在迎角大于30°后出现了极限环摇滚现象。     

典型超声速/高超声速动导数计算方法研究

以当地流活塞理论为基础完成超声速/高超声速动导数计算公式推导,发展了一种高效的动导数计算方法.选取1个超声速、2个高超声速动导数算例标模,计算了超声速/高超声速下飞行器动导数变化规律.研究结果表明:BFM超声速流动算例的俯仰组合动导数绝对值随马赫数的增大而减小,计算重心位置后移亦会导致俯仰组合动导数绝对值减小,纵向动稳定性降低;在0° ～20°迎角范围内,高超声速流动中的尖锥和钝头旋成体俯仰组合动导数绝对值随迎角的增大而增大,纵向动稳定性增强;所提方法计算时间约为双时间动导数计算方法的1/27.     

稳定性导数计算方法研究

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

鸭翼边条翼布局俯仰非线性抑制研究

采用CFD方法,研究鸭翼抑制边条翼布局俯仰非线性上仰的可行性。研究结果表明,边条翼布局具有复杂的多涡系旋涡流动结构,而导致俯仰非线性的主要因素是主翼涡的破裂。鸭翼对边条翼翼身组合体俯仰非线性上仰具有很强的抑制作用;鸭翼负偏度可提高鸭翼对俯仰非线性上仰的抑制能力,在更大迎角范围内对俯仰非线性具有明显的减缓作用。     

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

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

超音速、高超音速飞行器动导数的高效计算方法

将超音速、高超音速气动力工程计算方法综合起来,推广到非定常气动力的计算,并用于求解任意外形飞行器作强迫运动时的非定常气动力.在此基础上进一步发展了一种适合于超音速、高超音速大迎角状态的飞行器动导数计算方法和程序.计算结果与国际标准实验数据作了比较,证明了该程序的正确性和精度.研究了减缩频率和基本迎角对动导数的影响,结果表明基本迎角对动稳定性导数影响较大,而减缩频率在小迎角范围内对动导数影响很小.     

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

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

近耦合鸭式布局鸭翼展向吹气涡控技术数值模拟研究

采用鸭翼展向吹气间接涡控技术,对后掠角为50°的主翼和鸭翼气动布局进行数值模拟,给出不同吹气动量系数下的数值模拟结果,建立了鸭翼吹气动量系数与布局气动力系数之间的关系。并针对该布局模型,将风洞测力、测压以及水洞流动显示试验结果与数值模拟结果进行了详细的分析比较,结果表明,对鸭翼实施展向吹气技术,确实可以延迟和控制主翼涡破裂、增大升力的效果,在大迎角下把鸭翼作为涡发生器对主翼进行控制是可行的,计算结果与试验结果定性上是吻合的,是可以模拟这种复杂流场的。     

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

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

考虑进排气的飞翼静、动态气动特性数值模拟

采用计算流体力学（CFD）数值模拟技术,以二维模型为例,建立了基于小幅度非定常运动的动导数计算方法,构建了考虑喷流引射效应的进排气模型、考虑管道效应的通气模型以及常规的保形模型,分别对其静态以及动态气动特性进行了数值模拟计算,并辨识了其纵向俯仰力矩系数组合动导数。研究表明:非定常动导数辨识方法较为准确可靠,相比保形模型,进排气以及通气模型的静态失速攻角（AOAs）增大,升力系数、阻力系数及俯仰力矩系数也增大,但力矩系数斜率基本不变,说明静态稳定性差异不大。而相同攻角下,进排气影响下的俯仰力矩系数组合动导数绝对值最大,表明了进排气模型具有最大的动态阻尼,而通气模型次之,保形模型的动导数绝对值最小。      

飞翼布局飞行器大迎角非定常气动建模

为进一步研究复杂流场气动力模型和辨识方法,通过常规飞翼布局飞行器动态单独俯仰和单独滚转试验数据获得"同相"和"反相"导数。在采用状态空间法建立非定常数学模型的基础上,利用两步线性回归辨识出模型参数,给出带有误差限的特征时间常数分布。数学模型辨识结果与大迎角大振幅运动的试验数据对比表明,该模型建模结果基本上与试验气动数据变化一致,能够反映模型大迎角时的非线性和非定常特性,具有可用性。     

飞行器双自由度耦合振荡的数值模拟方法

针对模拟双自由度(2DOF)运动的"PQR"强迫振荡试验装置,通过坐标变换,将机构运动转化为空间欧拉角进行表示,实现对机构俯仰/偏航、俯仰/滚转和偏航/滚转等3种耦合运动的统一描述,采用预定运动轨迹的动网格计算技术,构造了基于非定常雷诺平均Navier-Stokes(URANS)方程的飞行器双自由度振荡数值模拟方法。围绕某复杂构型飞行器低速流动,在模型迎角10°、滚转角±40°、偏航角±40°振荡预设条件下分别对单自由度偏航、滚转振荡和双自由度偏航/滚转振荡进行了模拟,得到的气动系数迟滞回线在形态和量值上与风洞试验结果吻合,证明当前方法可以作为飞行器复杂耦合运动的一种有效研究手段。针对偏航/滚转耦合振荡,从运动形式和气动特性上与单独偏航、滚转运动进行了对比分析,结果表明,耦合运动在气流角、绕体轴的旋转角速度等方面与单自由度叠加效果不同,力矩系数迟滞回线存在曲线交叉的现象,表现出与单自由度振荡不同的阻尼特性。在当前模拟状态下,偏航/滚转振荡流场以横向分离涡运动为主要特征,可以推测大迎角情况下纵横向分离涡流动结构更加复杂,耦合作用更强,需要进一步在分离流模拟方法上开展研究。     

气动导数非定常数值模拟研究

非定常数值模拟飞行器的强迫振动过程,获得任意时刻的气动参数,可以解决复杂外形以及多种飞行条件的气动导数问题。通过AVICFD-X软件计算0?~80?初始攻角下的BFM带翼导弹的俯仰动导数,验证该方法的可行性。计算结果与实验值具有很好的一致性,反映了大攻角流场非线性特征。计算中充分考虑了工程实际需求,尽可能地减小计算代价,获得了理想结果,表明该方法适用于复杂外形气动导数计算。     

微型扑翼飞行器非定常运动对平尾的影响

以西北工业大学自行研制的微型扑翼飞行器ASN211为研究对象,利用其简化的二维扑翼及平尾串列翼模型进行了非定常数值模拟,分析了扑翼俯仰运动及沉浮运动对平尾气动性能的影响。在数值模拟模块中,模型的俯仰运动及沉浮运动由动网格技术实现。通过计算流体力学(CFD)软件Fluent对此非定常流场进行数值计算,重点研究了扑翼非定常运动尾流对平尾气动效率的影响。定常状态与非定常时均条件下平尾升力曲线的对比分析表明,扑翼的非定常运动能够增大平尾的失速迎角及最大升力系数,因而使平尾的失速特性得到改善。     

钝前缘三角翼无人机气动特性研究

三角翼布局因其优良的气动特性在军用飞机和无人机上获得了广泛应用.为了研究钝前缘三角翼无人机的气动特性,首先采用求解雷诺平均N-S方程的方法对NASA钝前缘三角翼标模进行对比计算,以验证计算方法的可靠度;然后对无人机四个升降舵偏角的气动力和流场特性进行分析研究.结果表明:三角翼无人机在升力系数较小时具有较高的升阻比,当迎角小于1 5°时,钝前缘三角翼前缘气流附体、吸力较高,翼面的横向流动不明显,使飞机的升阻比提高;当迎角大于15°后,涡流特征起主导作用,使得飞机在直到40°迎角范围内没有出现大面积气流分离,具有良好的俯仰稳定性,升降舵效率较高.钝前缘三角翼气动布局在翼展受限、翼载较小的条件下具有一定的气动特性优势.     

考虑多种物理效应的钝锥俯仰稳定性参数影响分析

首先采用二阶ROE格式计算了非定常NS方程下高超声速钝锥体的俯仰阻尼导数,通过与文献中内伏牛顿法和实验结果对比验证了本文所采用非定常计算方法的正确性.同时采用二阶ROE格式及双时间步长、刚性旋转动网格技术等,数值分析了NS方程计算条件下高空高马赫数钝锥体模型的强迫振动非定常流场特征,并运用最小二乘法辨识出其俯仰静稳定性导数与俯仰阻尼导数.通过计算结果对比分析了不同高度(含稀薄气体效应)、Maxwell滑移边界条件及五组元化学非平衡模型对钝锥体模型非定常流动特征及静、动导数的影响.结果分析表明,本文所计算的钝锥模型在所假设强迫振动下的俯仰稳定性参数随高度变化剧烈,且考虑滑移边界条件与化学反应会使得俯仰力矩迟滞曲线形态发生较大程度改变.准确计算高空、高马赫数条件下飞行器俯仰阻尼导数应考虑化学非平衡效应和稀薄气体效应的影响,是否考虑这些物理效应会关系到飞行器静、动稳定的程度.     

旋转流场下的振荡动导数试验技术研究

为研究飞机在旋转流场下的非定常气动特性,中国空气动力研究与发展中心低速所在Φ5m 立式风洞开展了旋转流场下的振荡动导数试验技术研究。本文推导了在旋转流场下识别组合动导数的方法,介绍了试验设备,获得了在旋转的同时,由振荡产生的3个组合动导数,并对试验结果进行了分析与讨论。将单自由度动导数结果与Φ3．2m 风洞试验结果进行了对比,旋转/振荡耦合试验结果表明：旋转运动使得俯仰组合动导数变得不稳定,而对于横向组合动导数,大转速则会显著增大非线性。该试验技术能够为研究旋转流场下的非定常气动特性提供一个有效的试验平台。     

操纵面作动对无尾布局无人机纵向气动特性的影响

通过风洞测力实验,研究了不同操纵面作动对某无尾布局无人机纵向气动特性的影响。实验结果表明：升降副翼以及襟副翼正向偏转都会使全机升力系数、阻力系数以及低头力矩增加。升降副翼作动引起的增量要高于襟副翼,并且舵偏角度越大增量越大。全动翼尖作动对全机纵向气动特性基本没有影响。在线性段,鸭翼作动对升力系数和阻力系数影响不大;线性段之外,鸭翼作动使得升力系数和阻力系数减小。迎角α〈16°以及α〉38°时,鸭翼正向作动使得低头力矩减小,负向作动使得低头力矩增加。操纵面作动对低头力矩的控制效率由高到低依次为：升降副翼、襟副翼、鸭翼和全动翼尖。进一步分析表明不同操纵面的控制效率与舵容量系数具有较大关系。     

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

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

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。