80°后掠三角翼强迫俯仰、自由滚转双自由度耦合运动特性数值研究

耦合求解NS方程和刚体动力学方程数值研究了80°后掠三角翼强迫俯仰、自由滚转双自由度运动特性。内容包括:对比双自由度俯仰/滚转运动和单自由度滚转及俯仰运动特性的差异,研究动态俯仰角对三角翼滚转运动振幅和频率的影响以及三角翼双自由度运动条件下的气动力特性。结果表明:在动态俯仰角的影响下,细长三角翼不能形成极限环形式的振荡;俯仰运动的迟滞效应造成滚转振幅的变化落后于俯仰角的变化;俯仰耦合滚转运动将会导致飞行器升力降低,横向稳定性变差。     

三角翼俯仰滚转耦合运动气动特性研究

介绍一套用于３ｍ 低速风洞的俯仰滚转两自由度大振幅非定常实验系统,并利用该系统对—三角翼单独俯仰和滚转及俯仰滚转耦合运动时的非定常气动特性进行了研究。结果表明,飞行器俯仰滚转耦合运动时的气动特性比单独俯仰和滚转时的气动特性复杂得多。     

某典型飞行器模型俯仰/滚转两自由度耦合动态气动特性

先进飞行器在大迎角高机动飞行过程中,气动力复杂演化与运动姿态剧烈变化的相互作用,容易诱发俯仰拉起的滚转或侧偏等耦合运动,呈现出十分强烈的多自由度动态气动力特性。为此,选取某典型四代机修型设计模型,利用FL-26风洞专用的俯仰/滚转两自由度动态试验装置,开展了单自由度和俯仰/滚转两自由度耦合运动的动态气动力特性研究。结果表明,研究模型在迎角30°左右会出现滚转单自由度的侧偏和小振幅摇滚运动现象;俯仰单自由度振动条件下,模型法向力和俯仰力矩会出现明显的气动迟滞特性;在强迫俯仰/自由滚转（或强迫滚转）两自由度耦合运动情况下,模型的侧偏运动幅值得到了一定的抑制,但纵横向气动力都出现了迟滞现象,甚至是多个气动迟滞环结构,容易引发耦合失稳发散。     

2.4m跨声速风洞大振幅动态试验技术

现代先进飞行器在高机动飞行过程中,容易诱发非指令的多自由度耦合运动,呈现出复杂的动态气动力及耦合运动现象。为满足先进飞行器多自由度动态气动特性研究与试验评估的需求,中国空气动力研究与发展中心高速空气动力研究所在2.4m×2.4m跨声速风洞建立了可以开展大迎角静态、单自由度俯仰振动、快速拉起、俯仰/滚转双自由度耦合运动等试验的模拟技术。通过典型的70°三角翼验证模型的试验研究,结果表明试验技术获取的试验数据合理可靠,变化规律正确,能够准确反映模型的动态气动力迟滞特性,实现了俯仰/滚转两自由度耦合大振幅运动的纵横向动态气动力测量,可以为飞行器的试验鉴定评估提供技术支撑。     

三角翼受迫俯仰滚转耦合运动的气动特性研究

基于刚性动网格的技术,选用B-L湍流模型,利用有限控制体积法对N-S方程进行数值离散,对76°大后掠三角翼的受迫俯仰滚转耦合运动进行了数值模拟,在此基础上,对俯仰滚转耦合运动的气动力特性和流场结构进行了分析。计算结果表明:俯仰滚转耦合运动时,三角翼上表面的涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩,滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大,但对法向力影响不大。     

细长三角翼滚转/侧滑耦合运动效应分析

通过耦合求解NS方程组和刚体动力学方程组,数值研究了80°后掠三角翼自由滚转/自由侧滑的耦合运动特性。通过与单自由度翼摇滚特性比较,发现滚转/侧滑耦合条件下三角翼的摇滚现象更为剧烈:在相同的滚转角下,耦合运动的滚转角速度更大,具有相同角速度时,耦合运动所达到的滚转角更大。研究表明:滚转/侧滑耦合条件下,涡流的非对称迟滞振荡仍是维持三角翼周期性等幅摇滚振荡的流动机理;摇滚过程中法向力的侧向分力是产生侧滑运动的主要原因;左滚右侧滑、右滚左侧滑是三角翼滚转/侧滑双自由度耦合运动的作动机制;侧滑运动的引入,导致迎角和侧滑角随滚转角迟滞变化,强化了三角翼摇滚过程中的迟滞效应,加剧流动结构的非对称特性。     

三角翼受迫俯仰滚转运动气动特性的数值模拟研究

基于动网格的技术，利用有限体积法对N—S方程进行数值离散，选用湍流模型，对65°大后掠三角翼的受迫俯仰滚转运动进行了数值模拟，在此基础上，对俯仰滚转耦合运动的气动力特性和流场结构进行了分析。计算结果表明：俯仰滚转耦合运动时，三角翼上表面涡分布的非对称性将产生横侧方向的偏航力矩和滚转力矩，滚转力矩和偏航力矩随着滚转振幅角和滚转缩减频率的增大而增大，但对法向力影响不大。     

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

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

俯仰—滚摆耦合复杂流场试验研究

给出 8 0°三角翼模型在低速风洞中实施俯仰 滚摆耦合大幅度振动的复杂流场测量结果 ,包括静态和动态、单自由度俯仰和滚摆、俯仰 滚摆耦合运动条件下的流动显示结果和强迫俯仰、自由滚摆条件下的滚转角时间历程测量结果。分析了耦合运动对前缘分离涡和气动特性的影响 ,以及不同相位角和不同频率对耦合气动特性的影响。结果表明 ,双自由度的耦合运动比单自由度运动具有更为明显的动态迟滞特性 ,进一步推迟分离涡的破碎 ,推迟开始滚摆的攻角αwr。     

细长三角翼滚转/侧滑耦合运动的数值研究

 采用亚迭代耦合求解流体动力学方程和刚体动力学方程数值研究方法，建立无时间滞后的耦合数值手段来分析研究飞行器非定常的耦合运动特征，研究细长三角翼多自由度非定常的耦合关联运动。研究表明：在流动结构失稳后机体运动逐步形成极限环振荡的自维持运动，滚转力矩系数滞回曲线呈现典型的“双8”稳定形态；在滚转和侧滑两自由度下，三角翼进行“落叶式飘荡”耦合运动，自激滚转振荡更为剧烈。     

一种近空间高超声速飞行器滚转稳定性研究

研究近空间高升阻比构型高超声速飞行器面临的复杂稳定性问题,对其机动飞行安全和操控具有重要意义。为此,通过静、动态风洞试验和数值模拟,开展了类HTV-2飞行器布局典型马赫数和迎角下,滚转方向静稳定性、单自由度和两自由度俯仰/滚转耦合运动动稳定性研究。结果表明：飞行器滚转方向为静稳定;自由振动的滚转单自由度运动为衰减振动曲线,动导数为负值,表现为动稳定;强迫俯仰/自由滚转两自由度耦合运动情况下,由于惯性力矩和非对称滚转力矩共同作用,导致滚转出现振幅周期性变化的极限环运动,甚至失稳发散。     

吸气式高超声速飞行器俯仰/滚转耦合运动特性

针对一种类似SR-72构型的吸气式高超声速飞机开展了进气道通流状态下俯仰/滚转耦合运动相关研究。通过数值模拟获得了滚转单自由度静稳定性、动稳定性以及强迫俯仰/自由滚转运动下的两自由度耦合动稳定性，研究了飞行器转动惯量以及俯仰运动频率对耦合运动的影响，简要分析了耦合运动的机理。研究发现虽然此飞行器具有滚转静稳定性和动稳定性，但是在强迫俯仰/自由滚转运动过程中，滚转通道却出现了小幅度振荡与大振幅振荡交替出现的情况，最大滚转角超过70°。小幅度振荡出现在正弦俯仰振荡的上半周期，其振荡频率随轴向转动惯量增加而降低，幅值随俯仰振荡频率增加而增大；大振幅振荡出现在下半周期，其幅值基本不变，而振荡频率与俯仰振荡一致。这种现象基本不受惯性耦合作用影响，可以认为是由气动力主导的。     

双三角翼及其翼身组合的滚转运动特性比较研究

采用Navier-Stokes方程与滚转运动方程耦合计算方法,比较研究了不同后掠角的双三角翼和翼身组合体的滚转运动特性,分析了机翼前缘后掠角及细长机身对非定常滚转力矩时滞环、动态流场结构和物面瞬时压力分布的影响。研究结果表明：主翼迎风面上的融合涡能量在80°/60°双三角翼上耗散较小,而在76°/40°双三角翼上耗散严重,这是造成两模型滚转力矩稳定性与时滞特性差异的主要因素;机身对气流的扰动作用,大幅增强了滚转力矩的线性分量;机身对气流的上洗作用,增强了边条涡与融合涡吸力及其时滞性,同时加剧了主翼背风面的两涡干扰;大滚转角时机身对横流流动的干扰,使得主翼背风面压力分布的时滞差异显著增加。该研究结果有助于认识后掠角与细长机身影响双三角翼滚转运动特性的物理机理。     

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

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

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

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

带副翼偏转的三角翼自由滚转运动数值模拟

 通过耦合求解非定常Euler/Navier-Stokes方程和单自由度滚转运动方程,对带副翼偏转的65°后掠角尖前缘三角翼WI1-SLE自由滚转运动进行了研究,Navier-Stokes方程的求解采用基于Spalart-Allmaras湍流模型的脱体涡模拟(DES)。在多块结构网格上,应用基于弧长的无限插值理论(TFI)生成变形网格,实现副翼偏转,而三角翼的滚转运动则通过网格的整体旋转实现。结果表明:Euler方程和DES方法均准确地模拟出了三角翼在滚转运动过程中存在的3个平衡位置。出现平衡位置的原因分别是:①流动对称性;②机翼左侧发生涡破裂的分离涡与右侧分离涡相互平衡使得滚转力矩为0,并且平衡位置仅与三角翼两侧涡强的差有关;③副翼偏转和左右机翼不对称分离涡涡强差产生的滚转力矩相互平衡。此外,滚转运动对副翼偏角幅值很敏感,幅值的微小改变会影响最终的平衡位置和向平衡位置运动的路径。     

细长三角翼摇滚运动数值研究

采用亚迭代耦合求解流体动力学方程和刚体动力学方程数值研究方法,建立流体与刚体运动耦合的无时间滞后的数值手段来分析研究飞行器非定常的耦合运动特征,预测失稳运动的发生与发展。研究表明:在流动结构失稳后机体运动逐步形成极限环振荡的周期性自维持运动,滚转力矩系数滞回曲线呈现典型的”双8”稳定形态,前缘涡的非对称飘起与再附是典型三角翼摇滚的驱动形式;三角翼构型的自由滚振幅随攻角变化呈现阶跃式的变化;自激振荡的幅度及频率与滚动惯量的大小存在对应的规律性关系。非定常自激振动的翼摇滚具有多控制参数、非稳定、非线性的运动特征。     

80°／65°双三角翼滚转稳定特性预测研究

在FL-23风洞中开展了80°／65°双三角翼大迎角下的滚转特性研究．包括静态测力试验．动导数试验和自由滚转试验,通过静态测力试验及动导数试验获得了双三角翼模型在大迎角条件下的滚转力矩特性以及动导数特性,从而对双三角翼大迎角条件下的滚转运动特性进行了预测．最后通过自由摇滚试验对预测结果进行了验证。研究结果表明随着模型迎角的增加,双三角翼呈现不同的滚转运动形态,包括静态稳定、舣周期震荡、准极限环摇滚．通过静态气动力及动导数可以较准确地对模型的运动形态及对应的迎角范围进行预测。     

静态滚转角对双三角翼涡流特性影响的数值研究

采用微可压缩模型（SCM）做为双三角翼低速流动控制方程,采用二阶Roe格式和LU-SGS方法离散控制方程。研究20°根弦名义攻角下,76°/40°双三角翼具有不同静态滚转角下的定姿态涡流特性。研究给出了不同滚转角下的空间涡流及表面极限流线,对滚转引起的攻角、侧滑角效应以及后掠角效应进行了分析,并对这些效应引起的流场结构、表面压力分布及气动力/力矩特性的变化进行了分析。结果表明：SCM能够较好地模拟低速双三角翼背风区的复杂涡流场,方法稳定性、收敛性较好 静态滚转导致双三角翼背风区涡流呈现复杂的非对称现象和涡破裂特征,并导致升力下降、横向稳定性变差。     

三角翼受迫滚转运动气动特性的数值模拟研究

基于动网格的技术,利用有限体积法对N—S方程进行数值离散,选用湍流模型,对65°大后掠三角翼的受迫滚转运动进行了数值模拟,在此基础上．对滚转运动的气动力特性和流场结构进行了分析。计算结果表明：滚转运动时．三角翼上表面涡分布的非对称性产生横侧方向的偏航力矩和滚转力矩,偏航力矩、滚转力矩系数滞回曲线呈现典型的“单环”形态,俯仰力矩、法向力滞回曲线呈现典型的“单8”形态。