颤振边界高效精确预测分析

高效精确地确定多种飞机构型的颤振边界在飞机设计过程中具有重要意义。为了提高计算效率和计算结果的准确性,针对亚声速和跨声速两种马赫数区域,提出分别采用线性和非线性方法进行非定常气动力分析。非线性分析在引入精确的定常气动力的基础上,采用高效率跨声速小扰动方程进行求解;颤振求解统一采用g 法。对大型飞机的梁架—减缩刚度组合模型的空机及三种典型燃油构型进行涵盖飞行包线的全马赫数变高度颤振分析,结果表明：四种构型的颤振边界与颤振试飞边界一致,与其他分析方法相比,效率有明显提高,尤其是对多种飞机构型能够高效地获得准确的颤振边界,即说明本文采用的方法是目前适用于工程上的一种高效精确的预测大型飞机颤振边界的方法。     

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

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

Investigation of Vortex Interaction in Canard-FSW Configurations Based on the Numerical Wind Tunnel Method

The vortex interference mechanism on low Reynolds number between the canard and main wing of the canard-forward sweep wing (Canard-FSW) configurations is simulated numerically by employing the numerical wind tunnel method. The variations of aerodynamic characteristics of Canard-FSW configurations with different positions of the canard are investigated, finding that the aerodynamic interference and mutual coupling effect between the canard and main wing have made great contributions to the lift and stability characteristics of the whole aircraft. Canard can radically improve the surface flow pattern of the main wing. And its own vortex can have a favorable interference on the main wing and can effectively control the airflow boundary layer separation. At small angles of attack, the aerodynamic characteristics are sensitive to the positions of the canard and the main wing, but at high angles of attack, the aerodynamic performances of the configuration are not only related to the shape of the canard (forward or backward), but also with the size of control force as well as the features of the vortices generated above the main wing and the canard. The different configurations and vortices are illustrated using the velocity vector, streamlines and pressure contours.     

飞行器动态稳定性参数计算方法研究进展

动态稳定性参数（简称动导数）是飞行器控制系统设计、飞行器动不稳定发生边界分析及相应动态稳定性判据研究的关键气动参数。在对飞行稳定性问题进行概述的基础上,介绍飞行器动态稳定性参数数值模拟的国内外研究进展。并按照理论方法、工程近似方法及计算流体力学（CFD）模拟方法的动导数发展方向对近年来主要的动导数计算方法进行了综述,评价了各种动导数预测方法的优缺点,指出了动导数数值模拟在理论基础、非定常气动力建模、预测方法精度和效率等方面存在的问题。最后对动导数数值模拟的发展趋势进行了展望。     

面向复杂构型飞机的非定常气动力建模与辨识

不论是现代高机动隐身战斗机的设计需求还是常规布局飞机的飞行动力学分析,深入研究大迎角飞行时的非线性非定常气动力模型都极其重要。基于纵向运动小振幅及大振幅强迫振荡试验数据,分析了常规稳定导数模型的准确性,并从导数模型出发发展了简化涡流和分离流时间迟滞效应的非定常气动力线性模型和非线性模型,最后应用风洞典型机动历程模拟试验验证了模型的有效性。结果表明：对于复杂构型高机动飞机模型,发展并改进的非线性微分方程模型可以准确预测飞机不同机动下的非定常气动力特性,具有较强的工程可行性。     

旋转机翼飞机旋翼模式前飞状态干扰气动特性

和传统直升机相比,旋转机翼（CRW）飞机在旋翼模式前飞时各部件之间存在更为严重的气动干扰。为了获得旋转机翼/机身/鸭翼/平尾之间的非定常气动干扰规律,基于运动嵌套网格技术,通过求解三维非定常雷诺平均Navier-Stokes（URANS）方程,建立了旋翼前飞流场数值模拟方法。首先对传统直升机旋翼/机身干扰模型进行了计算,验证了方法的可靠性,然后对某旋转机翼飞机全机在旋翼模式前飞状态下的非定常流场进行了数值模拟,并对各个气动部件上的非定常气动力和力矩的变化进行了分析。结果表明：飞机在旋翼模式前飞时,机身部件对旋转机翼的干扰较弱,在经过机身上方时拉力峰值仅略有增加;旋转机翼对鸭翼和垂尾干扰较弱,对机身和平尾干扰较强,随着前飞速度增大,旋转机翼对平尾的干扰会产生较大的升力损失和抬头力矩,需要引起重视。计算结果为该类飞行器的总体综合设计提供了参考。     

一种提高风洞动态试验数据质量的模型姿态控制和测量技术

低速风洞动态试验装置通常采用“电机+减速器+偏心机构+线性传动机构”的方式,存在传递环节多、机械间隙大等问题,给模型运动的精确控制和模型姿态的精确测量带来较大困难。为满足大型客机等大飞机对低速风洞动态试验的要求,基于现有的静态试验通用支撑平台,采用“电子凸轮”技术,建立了一套迎角/侧滑角解耦、可进行飞机小振幅动导数和大振幅非定常气动特性研究的动态试验装置。利用该动态试验装置进行了某飞机大尺度动态试验模型动导数试验和大振幅振荡试验,获得了试验数据的重复性,并研究了动态试验数据和静态试验数据之间的相关性。试验结果表明：利用该装置获得的某飞机动态试验数据重复性较好、规律合理,能满足大型飞机动态试验要求。     

基于总体布局参数的飞机静稳定性研究

基于飞机初步设计阶段所得到的总体布局参数,建立飞机气动导数数值模型,对亚声速静稳定性进行估算。最后对所建立的数值模型进行算例分析,并与计算流体动力学(CFD)方法所求的气动导数、实际飞行试验数据进行对比,结果表明,基于总体布局参数计算的气动导数方法快速、合理。     

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

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

基于CFD动网格技术的飞艇动导数计算方法

由于飞艇和飞机这两种飞行器在飞行原理、气动外形等方面都存在较大的差别,因此现有的飞机动导数工程估算公式和经验数据并不适用于飞艇动导数求解。为了解决上述问题,提出了一种新的基于CFD动网格技术的飞艇动导数计算方法。通过非结构化动网格技术,建立飞艇在流场中的强迫振动模型,采用有限体积法对N—S方程求解,研究飞艇在强迫振动中的气动特性变化规律,并基于数值计算结果求解飞艇的动导数。研究发现,得到的计算结果与风洞试验结果基本一致,表明基于CFD动网格技术进行飞艇的动导数求解是有效可行的。     

战术弹纵横向气动力计算方法

本文介绍一种计算中等攻角、中等侧滑角情况下正常式或鸭式布局的战术弹纵、横向气动力特性的方法。它是一种数值计算与工程计算相结合的方法。位流气动力通过求解线性偏微分方程获得,粘性分离产生的非线性气动力通过对一系列脱体涡的计算获得。方法中考虑了翼面与尾面间近距干扰和远距干扰的影响。本方法可用于计算来流马赫数为0到3.5左右、攻角和侧滑角为0°到20°左右、正常式或鸭式布局的战术弹纵、横向气动力系数。计算结果与实验数据比较,有良好的一致性。     

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

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

串列翼飞行器机翼耦合特性研究

串列翼飞行器由于其前后翼以及机身之间的相互干扰,气动特性复杂且难以预测。针对一款串列翼飞行器,以前后翼之间的垂直距离为变量,设计了五种气动布局,并使用CFD方法进行了数值模拟计算。通过对五种布局升阻特性与俯仰特性的比较及分析,发现前后翼垂直方向距离会显著影响整机升阻比、俯仰稳定性、气动中心位置以及压力中心位置。两翼间垂直方向上的距离越大,飞行器升阻比越高,且气动中心更加靠后。而在两翼间距离相同的情况下,前翼在下的布局拥有更高的升阻比,而前翼在上的布局拥有更好的俯仰静稳定性。     

鸭式旋翼/机翼飞机悬停状态飞行动力学特性

针对鸭式旋翼/机翼(Canard Rotor/Wing,CRW)飞机独特的气动布局,常规的分析方法及经验公式很难准确地对CRW飞机进行飞行动力学研究,通过飞行辨识对CRW飞机悬停状态特性进行了研究。首先,设计了飞行试验并获得了高质量的飞行数据,基于频率响应对CRW飞机的状态空间模型进行了简化。然后,在频域内对飞机的动力学参数进行了拟合优化,获得了CRW飞机悬停状态的动力学模型,并用飞行数据对所建模型进行了验证。最后,用辨识所得参数与常规直升机悬停状态时的参数进行了对比。结果显示悬停时CRW飞机的操纵导数和阻尼导数均比常规直升机小,经分析,操纵导数的减小主要是CRW飞机独特的旋翼设计所致,阻尼导数减小的原因主要是旋翼气动影响以及鸭翼、平尾、垂尾的结构影响。动力学特性分析结果为CRW飞机旋翼模式总体设计的进一步优化提供了指引和参考,所建立的模型可用于控制系统设计。     

DLR-F6/FX2B翼身组合体构型高阶精度数值模拟

基于雷诺平均Navier-Stokes(RANS)方程和结构网格技术,采用五阶空间离散精度的加权紧致非线性格式(WCNS)和剪切应力输运(SST)两方程湍流模型,开展了DLR-F6和DLR-F6_FX2B 2种翼身组合体构型的高阶精度数值模拟,计算外形来自AIAA第三届阻力预测研讨会。主要目的是确认WCNS模拟跨声速典型运输机构型和预测局部构型变化引起的气动特性变化量的能力。在固定升力系数条件下,采用粗、中、细3套网格开展了网格收敛性研究,从气动力系数、压力系数分布、表面流态等方面研究了网格规模对DLR-F6和DLR-F6_FX2B翼身组合体数值模拟结果的影响;采用中等网格开展了来流迎角对2种翼身组合体气动特性的影响研究。通过与National Transonic Facility(NTF)的试验结果和CFL3D的计算结果对比,表明采用高阶精度计算方法得到了网格收敛的数值模拟结果,较好地模拟了DLR-F6翼身组合体局部修型引起的微小气动特性变化和翼身结合部流动特性的差异。     

微型涵道飞行器飞行力学模型研究

微型飞行器(MAV)非线性飞行力学模型研究是MAV设计中的一个重要环节。微型涵道飞行器由于大包线、尺寸小、速度低、气动布局特殊,其飞行力学特性具有显著的非线性特性。以低雷诺数风洞实验为基础,研究了微型涵道飞行器的空气动力学特性,并采用CFD方法计算了微型涵道飞行器的动导数。在此基础上建立了微型涵道飞行器的飞行力学模型,并计算了基本飞行性能和配平曲线。结果表明,微型涵道飞行器与常规飞行器相比有很大差异,可以完成悬停、垂直起降和前飞的大包线飞行。     

折叠翼变体飞行器非定常气动特性实验研究

折叠翼变体飞行器是一种可以在飞行中改变自身气动外形的新型飞行器。研制出了一种折叠翼变体飞行器的风洞实验模型,在风洞实验中测得了模型不同变体位置下的气动力以及进行变体运动时气动力的动态变化过程,并通过PIV实验手段获得模型周围的流场在变体运动过程中的变化情况。结果表明：在机翼变形过程中,折叠翼模型有明显的非定常气动现象产生,而且折叠变形的速度越大,非定常现象越明显。出现非定常现象的主要原因是变体运动对机翼前缘涡的影响。     

Numerical Simulation of Rotor-aerodynamic Surface Interaction in Hover Using Moving Chimera Grid

Three-dimensional unsteady Navier-Stokes equations are numerically solved to simulate the aerodynamic interaction of rotor, canard and horizontal tail in hover based on moving chimera grid. The variations of unsteady aerodynamic forces and moments of the canard and horizontal tail with respect to the rotor azimuth are analyzed with the deflection angle set at 0° and 50°, respectively. The pressure map of aerodynamic surfaces and velocity vector distribution of flow field are investigated to get better understanding of the unsteady aerodynamic interaction. The result shows that the canard and horizontal tail present different characteristics under the downwash of the rotor. The canard produces much vertical force loss with low amplitude fluctuation. Contrarily, the horizontal tail, which is within the flow field induced by the down wash of the rotor, produces only less vertical force loss, but the amplitudes of the lift and pitching moment are larger, implying that a potential deflection angle scheme in hover is 50° for the canard and 0° for the horizontal tail.     

飞机大攻角空间机动气动力建模研究

 准确的气动力模型是大攻角机动仿真与控制等飞行动力学问题研究的基础。针对纵向运动建立的非线性非定常气动力微分方程模型扩展应用于一般空间运动形式,在分析大攻角非定常流动机理和试验数据的基础上,发展了一种多自由度非线性非定常气动力数学模型,并探讨了模型辨识和参数估计问题。用全机构形的大攻角动态风洞试验数据对模型进行的验证表明所发展的气动力模型能够描述大攻角气动力主要的非线性和非定常特性。     

DLR-F4翼身组合体的阻力计算

为了考察自行研发的CFD软件的计算能力和阻力计算精度，本文采用LU—SGS方法、MUSCL差分格式和Baldwin—Lomax代数湍流模型，数值模拟了AIAA阻力计算工作室提供的DLR—F4翼身组合体的绕流流场，综合分析了easel和case2的气动力的计算结果，并与NASA Christopher L．Rumsey采用CFL3D6．0和AFRL/VAAC Don W．Kinsey采用Cobalt60提供的两组计算结果以及AGARD提供的两种不同风洞的测力试验结果作了比较。计算结果表明，本文计算精度与国外CFD软件相当。为了提高激波，边界层干扰的模拟精度，今后要重点加强湍流模型的应用研究。