非稳定动态过程非定常气动力建模

现有的大迎角非定常气动力建模方法，通常是以一个或多个频率的稳定振动试验数据来预测稳定滞环。然而，飞机快速机动如过失速机动的过程，不可能是持续的稳定振动，而是一个非稳定的动态过程。因此，这个过程中的气动力不会达到稳定滞环，而是始终处于进入滞环的初始非稳定过程中。基于振动理论分析得出，非定常气动力的动态响应过程存在非稳定和稳定两个阶段，传统建模方法着眼于稳定阶段，而飞机的真实机动过程在非稳定阶段。设计了一种适于非线性系统辨识的激励输入，并以最小二乘支持向量机（LS-SVM）方法为例，实现了在大迎角区幅值和频率范围内任意运动的非定常气动力建模。模型训练完成后，用来预测某机翼在不同基准状态下大迎角范围内做俯仰运动时的升力系数、阻力系数和俯仰力矩系数。结果表明，不仅稳定滞环实现了准确预测，进入滞环的初始非稳定过程也得到了准确预测；此外，基准状态对气动力在初始非稳定过程中的特性存在明显影响。进一步的验证还表明，基于稳定滞环数据只能预测到稳定滞环，无法预测进入滞环的非稳定过程。     

三角翼俯仰振荡中若干参数影响的实验研究

本文研究一个 6 0°三角翼以不同振荡频率 ,不同振幅从不同起始迎角进行俯仰振荡的纵向气动力特性 ,研究发现跨越具有不同时间尺度的流态的振荡会产生明显的气动力迟滞现象 ,这一迟滞现象随减缩频率的增大有一个极限值 ;在同一种流态下振荡 ,气动力迟滞现象不明显 ,这一结果对充分利用非定常气动力有指导意义。     

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

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

基于多小波的Volterra级数非定常气动力建模方法

非定常气动力建模除了要准确描述气动力的非定常特性,还要反映其非线性特性。Volterra级数因为对系统非线性具有很强的描述能力正日益受到重视。一阶Volterra核只能表达线性特性,要建立反映非线性特性的非定常气动力模型,需要引入二阶核甚至更高阶核的影响。高阶Volterra核辨识的主要困难在于待辨识参数的数量随着核阶次的增加而呈指数增加导致计算难度急剧加大,即出现所谓维数灾难问题。以一种分段二次多小波为基函数将Volterra核展开,求解一个高维病态方程组来计算展开系数,利用小波的多分辨分析在时间和频率两个维度的分解特性将方程降维,最终将问题转化为求解一个低维方程组得到稳定解。通过NACA0012翼型在马赫数0.8下作沉浮运动时,升力系数、阻力系数和俯仰力矩系数的二阶核和三阶核的辨识构建非定常气动力模型,然后由此计算不同减缩频率下的气动力并与CFD结果进行比较,验证了Volterra级数对非线性非定常气动力的描述能力和多小波处理方法的有效性。     

基于CFD/CSD紧耦合的旋翼翼型气弹分析

采用CFD/CSD（计算流体力学/计算结构力学）紧耦合的方法,以Fluent软件作为主控平台,通过UDF（用户自定义函数）及I/O（输入/输出）文件读写的方式实现结构响应和气动载荷的数据交换,耦合求解了旋翼桨叶剖面的气动力和振动响应.在此基础上研究旋翼桨叶剖面在变距、沉降（挥舞）和周期交变来流条件下的气动特性和振动响应特性.结果表明:桨叶剖面在轻失速情况下,气动载荷周期性比较好,表现出光滑的迟滞环曲线,结构沉降响应也表现出光滑的周期性现象,扭转响应出现局部轻微振荡.深失速情况下,气动载荷及结构响应都表现出强烈的非线性振荡,高频成分较为明显.      

翼型动态特性数值模拟及其影响因素分析

采用非结构动网格技术模拟了NACA 0012翼型绕1/4弦长位置俯仰的强迫振荡过程,验证了翼型的周期性动态气动特性。分析了不同的湍流模型、单周期内步数和内迭代次数等因素对计算结果的影响;研究了马赫数和减缩频率对动态气动迟滞特性的影响。结果表明, Transtion SST模型得到的结果更准确;单周期内计算步数大于200时,计算结果收敛性更好;随着马赫数变大,迟滞效应先增大后减小;在亚音速区,迟滞环沿逆时针方向,动导数小于零,模型动稳定;而在跨音速区,迟滞环沿顺时针方向,动导数大于零,模型动不稳定;减缩频率越大,翼型动态迟滞效应越明显。     

飞行器自激振荡的流动物理与动力学机制

全局亚迭代耦合求解流体动力学方程和刚体动力学方程,研究轴对称飞行器单自由度俯仰失稳运动的非定常特征.高超音速锥-柱-裙飞行器的单自由度俯仰失稳运动发展为极限环形式,周期性运动伴随着波系结构的非定常变化;自激振荡源于静稳定构型的振动回复机制和尾裙激波流动迟滞效应所构成的阻尼机制.基于第二拉格朗日方程和虚功原理,导出能够描述迟滞现象的参数化非线性动力学模型.多尺度近似分析获得参数化运动特征:自激振动过程是拟简谐运动;静稳定构型才能出现自激振荡周期运动;平衡点阻尼是决定运动稳定特性的分叉参数;振幅特性与阻尼非线性相关,频率特性与刚度非线性相关.基于数值结果的参数辨识,非线性模型的理论分析、重构都与数值结果高度一致,从而有效地佐证了自激振荡建模研究的合理性.     

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.     

Unsteady aerodynamic modeling at high angles of attack using support vector machines

Accurate aerodynamic models are the basis of flight simulation and control law design.Mathematically modeling unsteady aerodynamics at high angles of attack bears great difficulties in model structure determination and parameter estimation due to little understanding of the flow mechanism.Support vector machines(SVMs)based on statistical learning theory provide a novel tool for nonlinear system modeling.The work presented here examines the feasibility of applying SVMs to high angle-of-attack unsteady aerodynamic modeling field.Mainly,after a review of SVMs,several issues associated with unsteady aerodynamic modeling by use of SVMs are discussed in detail,such as selection of input variables,selection of output variables and determination of SVM parameters.The least squares SVM(LS-SVM)models are set up from certain dynamic wind tunnel test data of a delta wing and an aircraft configuration,and then used to predict the aerodynamic responses in other tests.The predictions are in good agreement with the test data,which indicates the satisfying learning and generalization performance of LS-SVMs.     

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

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

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

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

适用于气动弹性的小波非定常气动力降阶方法

发展非定常气动力模型降阶技术旨在缩减计算耗费并且使得计算流体力学信息能够应用于气动伺服弹性以及设计优化当中。采用小波方法建立基于Volterra级数的非定常气动力降阶模型。在模型识别过程中,激励是光滑连续的且以零阶保持方式进行离散,采样频率选择二的整数幂次。近似的一阶Volterra核基于Haar尺度函数族展开,鉴于Volterra核在系统响应中的衰减特性,可在合适的有限时长截断一阶核。为了获得一阶核的各项展开系数,需要求解由输入/输出数据组成的超定方程,其中涉及到奇异值分解算法。为了验证小波方法的有效性,算例选取了二维的NACA64a010翼型。数值仿真结果表明该方法能够比较准确地预测结构小扰动引起的非定常气动力响应且能描述一定的非线性现象。     

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

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

非定常气动力非线性微分方程建模方法

为准确描述过失速机动中非定常气动力特性,研究了以非线性微分方程为基本结构的非定常气动建模方法.基于动力学系统建模思想,分析揭示该模型的物理机理,并发展和改进了基于风洞强迫振荡试验的模型参数辨识方法:基于小振幅试验数据,采用线性回归参数辨识方法辨识确定气动模型中特征时间常数等线性参数;基于大振幅试验数据,采用遗传算法全局...     

民机极限飞行状态的动态气动力试验与建模

飞行失控是造成民机灾难性航空事故的重要因素,飞行失控中飞机难以避免超出正常飞行包线范围,进入具有复杂非线性和非定常动态气动特性的极限飞行状态。本文开展典型民机布局飞机极限飞行状态的动导数、大振幅试验,对大迎角动态气动力的参数影响规律以及非线性、非定常特性进行分析和建模。结果表明,在飞机失速到过失速区域,飞行姿态快速变化过程中动态气动力的非线性和非定常特征显著;在动导数试验和建模中,考虑运动角速率的影响,可以预示气动力非线性的迎角范围,并捕捉到关于飞机动稳定性演化的关键特征;利用Goman-Khrabrov状态空间模型结合大振幅试验,可以确定模型中表征非定常特征的关键时间常数,获得特定极限飞行状态运动中的非定常动态气动力特性。研究方法和结果为开展民机极限飞行状态的动态气动力风洞试验设计与建模提供了一个可行途径,能改进飞机飞行失控预防、极限状态改出、飞行模拟训练和飞行事故分析等。     

基于分层思路的动态非线性气动力建模方法

很多非线性气动力模型难以精确预测系统的小扰动线性特征。针对这一局限,提出了一种非线性分层模型,用于辨识跨声速非线性非定常气动力。分层建模需要同时提供微幅振荡和大幅振荡两套训练样本,首先通过线性模型(如带外输入的自回归(ARX)模型)对微幅振荡样本进行建模,而后采用非线性模型(如径向基函数神经网络(RBFNN))对大幅振荡的样本与线性模型的差量进行建模,进而把线性模型和非线性模型的输出相叠加,得到分层非线性动力学模型。算例表明建立的分层气动力模型与单一自回归径向基函数(AR-RBF)神经网络模型相比不仅具有更高的数值精度,可以精确预测大幅运动中的非线性特征,而且在小扰动状态下自动退化为线性模型,能够精确刻画结构微幅振荡下的线性动力学特性。     

非定常气动力的结构自适应神经网络建模方法

讨论了非线性非定常气动力的结构自适应神经网络模型建模方法,该方法具有同时进行结构辨识与参数辨识的优点;利用纵向大振幅强迫振荡风洞试验数据,验证了建模方法及所建模型的有效性,结果表明:结构自适应神经网络模型对非定常气动力有很好的逼近能力,由于采用飞行参数的时间离散数据作为输入量,模型可直接应用于飞行仿真研究。     

Reduced order model for unsteady aerodynamic performance of compressor cascade based on recursive RBF

Based on Recursive Radial Basis Function (RRBF) neural network, the Reduced Order Model (ROM) of compressor cascade was established to meet the urgent demand of highly efficient prediction of unsteady aerodynamics performance of turbomachinery. One novel ROM called ASA-RRBF model based on Adaptive Simulated Annealing (ASA) algorithm was developed to enhance the generalization ability of the unsteady ROM. The ROM was verified by predicting the unsteady aerodynamics performance of a highly-loaded compressor cascade. The results show that the RRBF model has higher accuracy in identification of the dimensionless total pressure and dimensionless static pressure of compressor cascade under nonlinear and unsteady conditions, and the model behaves higher stability and computational efficiency. However, for the strong nonlinear characteristics of aerodynamic parameters, the RRBF model presents lower accuracy. Additionally, the RRBF model predicts with a large error in the identification of aerodynamic parameters under linear and unsteady conditions. For ASA-RRBF, by introducing a small-amplitude and high-frequency sinusoidal signal as validation sample, the width of the basis function of the RRBF model is optimized to improve the generalization ability of the ROM under linear unsteady conditions. Besides, this model improves the predicting accuracy of dimensionless static pressure which has strong nonlinear characteristics. The ASA-RRBF model has higher prediction accuracy than RRBF model without significantly increasing the total time consumption. This novel model can predict the linear hysteresis of dimensionless static pressure happened in the harmonic condition, but it cannot accurately predict the beat frequency of dimensionless total pressure.     

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

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

翼型舵面偏转非定常流动数值模拟

针对NACA0012翼型舵面偏转问题,数值模拟了不同参数对翼型气动特性的影响。基于非结构动网格技术,采用ALE有限体积描述下的二维可压缩非定常N-S方程,计算通量采用Vanleer格式、时空二阶格式,利用Venkatakrishnan限制器抑制数值振荡。非定常计算结果表明,NACA0012翼型绕1/4弦点作周期性俯仰振动的升力系数和俯仰力矩系数结果与实验数据吻合良好,验证了数值方法的准确性;在翼型舵面表面有分离区产生,升力系数和俯仰力矩系数形成滞回环,在亚声速情况下,滞回环幅值较小,进入超声速阶段以后,幅值增大,随着翼型间缝隙宽度逐渐增加,翼型升力系数和俯仰力矩系数与无缝翼型相比逐渐降低。