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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

一种高效的组合动导数分量数值模拟算法

基于传统网格速度法,发展了在流场网格不变的情况下,数值模拟沉浮运动、纯俯仰运动和俯仰振荡运动的算法,提出了一套直接求解旋转阻尼导数和洗流时差导数的方法,该方法不需要实时更新网格,减少了计算时间和所需内存,避免了负体积的出现。首先通过NACA 0006的阵风响应和NACA 0012翼型的强迫俯仰振荡,验证了程序和方法的正确性;然后数值模拟了BFM(basic finner missile)标准模型的沉浮运动、纯俯仰运动和俯仰振荡运动,并计算了其旋转阻尼导数、洗流时差导数和组合动导数,所得结果与动网格方法和风洞试验值吻合,最大误差不超过4%,进一步验证了该方法的正确性。     

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

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

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

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

考虑地面效应的翼型动态特性数值模拟

为了研究考虑地面干扰的飞行器非定常气动特性,基于滑移网格技术,通过求解二维非定常Euler方程,对不同近地高度下的NACA0012翼型的动态振动特性进行了数值模拟。分析了存在地面干扰时的升力系数以及俯仰力矩系数迟滞环变化,并结合强迫振动法,计算了翼型各个近地高度下的纵向组合动导数。结果表明,地面效应不仅对定常流场产生影响,更显著地影响了非定常气动力及力矩,近地高度越小,这种效应越明显,使得升力系数的迟滞环面积变小,幅值减小,而俯仰力矩系数的迟滞环变化不规律,力矩系数的动导数随着近地高度减小而增大(绝对值减小),纵向的动态稳定性产生损失。地面效应干扰在飞行器非定常气动研究中应该引起重视。     

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

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

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

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

稳定性导数计算方法研究

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

窄条翼导弹俯仰机动中滚转失稳及其控制过程

窄条翼布局导弹通常具有复杂的横向气动特性,在大迎角飞行及快速机动中很容易诱发出现滚转非指令偏离和连续振荡,可能导致飞行失控,影响落点精度。为了研究窄条翼导弹俯仰快速机动对滚转失稳的诱发过程及滚转失稳对俯仰机动控制效果的影响,并验证三通道解耦控制方法的有效性,针对典型俯仰机动过程,分别利用2.4 m跨声速风洞虚拟飞行试验平台和耦合气动/运动/控制的一体化数值计算方法开展了相关研究。结果表明,风洞试验和数值模拟均成功预测了俯仰拉起和保持过程中的滚转自激失稳运动及其引起的纵、横向耦合运动,针对该机动过程,三通道解耦控制方法能够有效抑制滚转运动,保持姿态稳定。     

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

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

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

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

非定常超音速侧滑机翼的升力、力矩系数及动导数的数值计算

 本文利用分布在机翼平面上的非定常超音速点源的特征线网格法给出非定常超音速侧滑机翼的数值解。本文给出侧滑机翼作沉浮、俯仰、滚转谐振时的振型以及计算升力、力矩和动导数（包括交叉导数）的计算方法。按本文用非定常理论数值解对典型的无侧滑机翼计算出的动导数与常用的按准定常理论解析公式计算结果进行比较和鉴定,结果表明本文方法是令人满意的。     

不同攻角条件下偏转头弹箭流场与气动特性

基于DES（detached eddy simulation,分离涡模拟）方法及AUSM（advection upstream splitting method,对流迎风分裂方法）格式,对不同偏角的偏转头弹箭在非零攻角下的超声速流场和气动特性进行了数值模拟.结果表明:在非零攻角情况下,与常规弹箭流场相比,偏转头弹箭头部迎风面激波更强;而当有效入射角相等时,即使弹箭的攻角不一样,其头部流场结构仍相同,且尾翼流场基本不受头部偏角影响.另外,由气动力系数变化曲线可知,偏转头弹箭具有良好的纵向静稳定性,且升力系数随攻角的增长率大于阻力系数,即升阻比大于常规弹箭,因此适当增大攻角可以提高偏转头弹箭的气动特性.