耦合运动非定常气动模型对飞机飞行特性仿真的影响

 在大振幅偏航滚转单自由度强迫振荡和耦合强迫振荡运动风洞试验的基础上,分别运用准定常气动力建模方法和非定常气动力建模方法,获取飞机的气动力模型;运用两种不同的气动力模型,对飞机的纵向正弦振荡机动飞行进行了仿真。结果表明,准定常模型在小迎角下的仿真结果与非定常模型的结果基本一致;在较大迎角时,准定常模型仿真结果趋于稳定,而非定常模型仿真结果使飞机进入横航向振荡飞行状态,这一现象与F-16XL的试飞结果类似,其主要原因是耦合运动的非定常模型与准定常模型所产生的阻尼特性不一致。这也表明,分析大迎角飞行特性时应考虑耦合运动下的非定常气动特性。     

跨声速机翼非定常气动力的全位势粘位迭代计算

采用Ｃ－Ｈ型网格,守恒型非定常全位势方程的时间精确近似因式分解差分地计算二维,三维的跨声速非定常位势流,用准定常,准二维方法计算边界层位移厚度,通过粘位迭代得到的跨声速翼型,机翼的非定常气动力,所得结果与实验数据吻合很好。     

F12全机动态特性数值模拟

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

基于涡尾迹方法的风力机载荷与响应计算

采用涡尾迹方法计算准定常气动性能,结合Leishman-Beddoes动态失速模型计算叶片的非定常气动力。应用有限元分析软件,分析风力机在不同转速下的风轮模态,塔架的前后、左右振动模态。基于非定常气动力和结构振动模态,针对风力机风轮,建立其结构运动方程。求解出的叶片振动速度加入到非定常气动力的计算中,从而建立了考虑叶片振动的风力机载荷、响应计算模型。分析了Phase VI叶片和某1.5MW风力机各个叶片截面的气动力、振动位移和振动速度随时间的响应曲线,表明该方法能够较好地计算出风力机的气动性能、载荷和响应。     

旋转总压畸变对压气机稳定性影响的三维非定常模型

基于三维非定常激盘模型,发展了一个可用于预测压气机动态失速特性及旋转总压畸变对压气机气动稳定性影响的三维非定常计算模型,并利用该模型计算分析了旋转总压畸变的幅值和旋转频率变化对一台双级低速轴流压气机性能和气动稳定性的影响.计算结果表明,旋转总压畸变的幅值和旋转频率对压气机的性能、动态失速过程特性和失速边界有强烈的影响;...     

直升机旋翼/机身非定常气动干扰数值分析

基于非定常面元/时间步进全展自由尾迹建立了旋翼/机身非定常气动干扰分析方法。方法中耦合了非定常面元法和时间步进自由尾迹,以准确模拟旋翼非定常气动力、旋翼尾迹及桨叶对机身的非定常干扰效应。为模拟前飞状态下具有升力的机身,将机身离散为汇/偶极子面元,并采用涡线镜像法模拟旋翼尾迹靠近机身表面产生的加速效应。通过计算前飞状态的Maryland、ROBIN(Rotor Body INteraction)旋翼/机身干扰下的非定常压力分布,并与可得到的实验值、CFD计算结果对比,验证方法的准确性。随后分析前飞速度、旋翼与机身距离对旋翼/机身非定气动干扰的影响。计算结果表明机身头部和中部非定常压力主要受桨叶的通过性影响,而机身尾梁主要受尾迹/机身干扰影响,机身非定常气动力频率为桨叶片数的倍频。随前飞速度的增加,机身非定常压力幅值增加,尾梁压力幅值先增加后减小;增加旋翼与机身距离将减小机身和尾梁非定常压力幅值。     

非结构动网格技术及其在超声速飞行器头罩分离模拟中的应用

针对激波诱导物体运动的二维动力学模型,比较了采用流动定常、结合非定常弹道方程的所谓"准定常"计算方法和非结构动网格计算方法,通过流场和运动特性分析,指出发展有相对运动的非定常流动模拟方法的重要性;另外采用修正弹簧近似和局部重构的非结构动网格技术,耦合求解Euler方程及6DOF弹道方程,数值模拟高度H=25km、飞行马赫数Ma=6.0、迎角α=0°状态下火箭抛罩的动力学过程.研究表明,稠密大气层内头罩受到的气动力比惯性力大几个数量级,导致头罩的运动特性完全由气动力控制.分离过程中,头罩和火箭头部之间的气动干扰在时间和空间上都表现出强烈的非线性特征,头罩姿态变化剧烈.本文还对分离过程中部分位置上的流场进行定常模拟,比较了定常和非定常的气动力,所得结论对于分析地面风洞可控轨迹系统(CTS)试验具有指导意义.     

螺旋桨滑流与机翼气动干扰数值模拟研究

分别采用非定常方法和激励盘理论研究螺旋桨滑流对机翼的干扰作用,并对计算结果进行对比分析。非定常方法基于动态面搭接结构化网格,运用有限体积法求解非定常RANS方程。对于激励盘理论,首先采用点对点结构化网格,随单独螺旋桨旋转的方式进行非定常模拟得到螺旋桨桨叶压力分布,然后将桨叶压强分布转化为激励盘载荷分布,最后将载荷分布作为边界条件对干扰流场进行定常模拟。数值结果表明:除了在较大来流攻角时有所差异,两种方式得到的机翼气动力吻合良好;非定常方法可以得到螺旋桨与机翼的瞬时干扰流场,而激励盘理论则可以大大缩短数值模拟的时间并降低计算要求。     

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

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

可变形翼型超声速非定常气动力的研究

由小扰动线化方程出发,利用基于格林公式的边界元方法求解二维超声速流动速度场,进而对变形翼型气动特性进行理论分析和数值计算。计算结果显示虚拟质量力导致的非定常附加升力对厚度变化远比弦长或弯度变化敏感;弦长或弯度变化诱导的非定常效应影响很小,此时非定常气动升力和准定常气动升力近似相等,和不可压缩或亚声速时相类似,变形翼型的非定常气动升力近似等于准定常计算结果叠加上非定常附加升力,与具体的变形历史过程无关。同时还分析了超声速时翼型往复变形时的气动特性,和亚声速情况完全不同,随着马赫数升高,完成往复变形所需要外界输入的功逐渐减小,趋近一个常值。另外翼型完成往复变形所需做的功将不依赖于来流攻角,且与翼型参数的变化加速度a存在简单的数学关系:A ∝a～(1/2)。     

利用转动前缘的动边界效应控制翼面分离的数值模拟

本文使用数值模拟方法研究了活动翼面对翼型背风面分离的控制作用。数值实验的模型为一头部装有一可旋转圆柱的对称Joukowsky翼型。使用的数值方法为在贴体网格上用有限差分法求解可压流的NS方程。数值实验表明,当头部圆柱以一定角速度旋转时,其升力比不旋转时有明显的增加,分离也受到一定程度的控制。     

变展长战术导弹非定常气动特性研究

以某轴对称基准弹为模型,通过数值模拟对其进行了准定常和非定常气动特性研究,探索可变翼飞行器在变展长过程中的非定常气动特性。结果表明:在弹翼连续变展长的过程中,非定常气动特性曲线在对应的准定常气动特性曲线周围形成滞回环,若来流马赫数增大,则滞回环变小。当来流为超声速时,轴对称基准弹飞行迎角和变形周期对其非定常气动特性影响较小,而非定常气动特性并不显著,接近于准定常气动特性,在可变翼飞行器设计过程中可以忽略不计;当来流为亚声速时,展长变化对非定常气动特性有一定的影响,不能忽略。     

Mechanism analysis of hysteretic aerodynamic characteristics on variable-sweep wings

Variable-sweep wings have large shape-changing capabilities and wide flight envelops, which are considered as one of the most promising directions for intelligent morphing UAVs. Aerodynamic investigations always focus on several static states in the varying sweep process, which ignore the unsteady aerodynamic characteristics. However, deviations to static aerodynamic forces are inevitably caused by dynamic sweep motion. In this work, first, unsteady aerodynamic characteristics on a typical variable-sweep UAV with large aspect ratio were analyzed. Then, deep mechanism of unsteady aerodynamic characteristics in the varying sweep process was studied. Finally, numerical simulation method integrated with structured moving overset grids was applied to solve the unsteady fluid of varying sweep process. The simulation results of a sweep forward-backward circle show a distinct dynamic hysteresis loop surrounding the static data for the aerodynamic forces. Compared with the static lift coefficients , at the same sweep angles, dynamic  lift coefficient in sweep forward process are all smaller, while dynamic sweep backward  lift coefficient are all larger. In addition, dynamic deviations to static  lift coefficient are positively related with the varying sweep speeds. Mechanism study on the unsteady aerodynamic characteristics indicates that three key factors lead to the dynamic hysteresis loop in varying sweep process. They are the effects of additional velocity caused by varying sweep motion, the effects of flow hysteresis and viscosity. The additional velocity induced by sweep motion affects the transversal flow direction along the wing and the effective angle of attack at the airfoil profile. The physical properties of flow, the hysteresis and viscosity affect the unsteady aerodynamic characteristics by flow separation and induced vortexes.     

Aerodynamic modelling of insect-like flapping flight for micro air vehicles

Insect-like flapping flight offers a power-efficient and highly manoeuvrable basis for a micro air vehicle capable of indoor flight. The development of such a vehicle requires a careful wing aerodynamic design. This is particularly true since the flapping wings will be responsible for lift, propulsion and manoeuvres, all at the same time. It sets the requirement for an aerodynamic tool that will enable study of the parametric design space and converge on one (or more) preferred configurations. In this respect, aerodynamic modelling approaches are the most attractive, principally due to their ability to iterate rapidly through various design configurations. In this article, we review the main approaches found in the literature, categorising them into steady-state, quasi-steady, semi-empirical and fully unsteady methods. The unsteady aerodynamic model of Ansari et al. seems to be the most satisfactory to date and is considered in some detail. Finally, avenues for further research in this field are suggested.     

滑动蒙皮变后掠无人机非定常气动特性研究

为了探索飞机变形过程中的非定常气动特性及其机理,基于滑动蒙皮变后掠无人机这一可变形飞行器大尺度全局变形研究平台,通过风洞实验对其进行了非定常气动特性研究。结果表明:无人机变后掠过程中非定常气动特性曲线在相应的准定常曲线周围形成滞回环;无人机变后掠速率越快,滞回效应越显著;无人机以约7.5°/s的速率变后掠,会使非定常气动特性数值偏离相应的准定常数值5%以上;变后掠无人机非定常气动特性产生的主要原因可能在于流场结构迟滞和机翼附加速度。     

双旋翼微型飞行器的非定常气动特性

发展了一种基于计算流体力学(CFD)的非定常气动特性预测方法,计算方法包括了动量源模型、预处理方法、非结构嵌套网格和Spalart-Allmaras(S-A)湍流模型等技术.通过计算悬停Caradonna-Tung算例和俯仰振荡NACA0012算例,验证了计算方法模拟双旋翼微型飞行器动态流场的有效性.数值模拟了双旋翼微型飞行器动态流场,给出了非定常气动系数的迟滞曲线,分析了缩减频率、前飞速度和螺旋桨转速对非定常气动特性的影响.计算结果表明:力矩系数迟滞效应随缩减频率、前飞速度和螺旋桨转速增大而增大,升力系数迟滞效应随缩减频率和前飞速度增大而增大,但基本不随螺旋桨转速变化而变化.      

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

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

三维旋转效应对桨叶气动特性影响的计算

在旋转的作用下,流动分离向下游延迟,结果大大改变了风力机桨叶的气动特性,这种由旋转所带来的影响被称为三维旋转效应。本文以桨叶表面边界层理论分析的结果为基础,给出了三维旋转效应对气流分离影响的解析关系,将其与非定常气动和动态失速模型相结合,并适当考虑旋转条件下桨叶翼型载荷的分布特性,使其能计算三维旋转效应下风力机桨叶的非定常气动特性.通过算例表明,当这一三维旋转效应模型与非定常动态失速模型耦合后再计算桨叶非定常气动特性时,能大大改善计算结果。     

冲压翼伞流场与气动操纵特性的数值模拟

采用有限体积法求解shear stress transport(SST) k-ω二方程湍流模型下的Navier-Stokes(N-S)控制方程,对冲压翼伞的气动特性进行数值模拟,分析翼伞的流场机理和气动操纵特性.模拟得到的升阻特性与试验数据较吻合,在此基础上分析前缘切口、弧形下反以及稳定幅对升阻特性的影响.通过软件Fluent实现了非定常流动模拟,研究了翼伞的非定常升力特性,其升力系数的脉动受脱体涡的非定常过程影响,当迎角为16°时,翼伞升力变化周期为0.36s.最后分析了翼伞稳定滑翔阶段的纵向静稳定性,相比于单边后缘下拉方式,通过闭合翼伞一侧进气口实现航向操纵更稳定有效.