基于降阶模型的翼型颤振主动抑制研究

为了解决CFD/CSD计算效率低下的问题,基于气动力降阶模型建立开环气动弹性分析系统,耦合操纵面模态对应的气动状态空间方程,发展了一套伺服气动弹性分析系统,设计主动控制律从而进行二维翼型颤振抑制。结果表明,基于ROM计算的无量纲颤振速度准确性及精度高;使用开、闭环气弹系统时基于ROM和CFD分别仿真时域仿真结果吻合较好;基于ROM闭环气弹系统的NACA0012翼型气弹颤振边界比开环气弹系统提高了约20%。研究结果为有效抑制翼型颤振提供一种快速计算和设计方法。     

基于ARMA/ROM的颤振边界模拟和快速预测

为实现对颤振边界进行快速预测,采用基于系统辨识的降阶气动力方法,将ARMA/ROM训练数据用于构建闭环气动弹性系统的连续时间状态空间模型,并将所得到的状态空间模型用来确定颤振边界;对AGARD445.6机翼进行了颤振边界和颤振频率的模拟,将ARMA/ROM方法的计算结果与全阶耦合计算结果、实验结果进行了对比。研究结果表明,该方法是高效、精确的,可以用于飞行器颤振边界的模拟和快速预测。     

操纵面对跨声速机翼气动弹性特性的影响

 运用基于非定常CFD的气动力辨识技术，得到跨声速非定常气动力降阶模型。耦合结构动力学方程，建立了基于状态空间的跨声速气动弹性分析模型。分析了典型三自由度二元机翼的颤振边界，分析结果与CFD/CSD直接耦合方法吻合。然后研究了操纵面结构参数（固有频率和重心位置）对跨声速气动弹性特性的影响。研究发现,一些传统的结构设计准则和颤振排除技术未必适用于跨声速状态;操纵面偏转模态常常成为诱发跨声速颤振的主要模态;经典的质量平衡技术可能会降低跨声速气动弹性系统的稳定性。     

基于模态耦合方法的机翼颤振数值模拟

针对机翼气动弹性问题,发展了一种流固耦合模拟技术,将固体结构的广义位移表示为结构固有模态的线性组合,求解结构的运动,再与流体进行数据耦合、迭代计算.同时为了确定机翼颤振边界,发展了一种线性插值确定颤振边界方法.采用标准气弹机翼AGARD445.6算例,验证所采用的模态耦合技术和颤振边界确定方法的有效性及工程实用意义.     

非定常气动跨声速机翼颤振分析

颤振是气动弹性分析中影响最大的问题,在飞行器设计中得到了极大的重视。目前,在飞行器气动弹性分析中普遍采用基于偶极子网格法空气动力的机翼颤振计算方法。本文采用p-k法,以AGARD445.6机翼标准模型为算例,计算其在不同高度与飞行马赫数下的颤振速度和频率特性曲线,发现机翼在跨声速区域存在明显的"凹坑"物理特性,并且这些计算结果与试验值较为吻合。     

地效对二维气动弹性特性的影响

为研究地面效应对飞行器的气动弹性特性的影响,基于计算流体力学(CFD),耦合求解结构动力学方程,建立了气动弹性求解系统。通过涡诱导的圆柱自由振荡算例和Isogai的二元气动弹性算例,分别研究了地面效应对气动弹性动态响应和颤振边界的影响。研究结果表明:在地效作用下,涡街的对称性及发展受到影响,圆柱的动态响应幅值减小,频率增加,越靠近地面影响越显著;地效作用下,二维翼型在亚声速阶段的颤振速度下降,跨声速阶段,颤振边界凹坑范围扩大,颤振速度与激波位置密切相关。     

基于多重网格方法的跨声速颤振数值模拟研究

使用多重网格的方法对二维翼型和三维机翼的跨声速颤振进行了数值求解.流场控制方程为N-S方程,湍流模型采用SA模型网格.计算网格采用结构网格,空间离散采用中心格式,使用双时间法进行时间推进.流场与结构之间的数据通过径向基函数(RBF)插值方法进行交换.通过耦合求解流场控制方程和结构运动方程得到Isogai(case A)机翼模型和Agard445.6机翼的跨声速颤振边界.通过与使用单重网格时得到的结果相比较得出,多重网格方法能够节省计算时间,提高颤振分析计算效率.     

翼梢装置对某机翼气动弹性行为影响研究

以大客某方案机翼为基本翼,通过数值模拟的方法研究了翼梢装置对机翼气动弹性特性影响,包括静气动弹性及颤振特性。其中通过CFD/CSD弱耦合求解的方法研究其静气动弹性响应,气动力计算采用面元法,结构响应计算采用结构有限元法,通过插值实现翼面气动力与有限元节点力之间的传递,以及有限元模型与气动网格之间的变形传递。对基本翼及带翼梢装置机翼静力学有限元模型局部修改得到动力学模型,应用MSC NASTRAN进行颤振特性分析。研究发现翼梢装置使得机翼的气动弹性特性不同程度均有降低,而不同翼梢装置对其影响又有所不同,可见,翼梢装置的设计在追求气动特性改善的同时必须关注其带来的结构特性的损失。     

变弯度机翼参数化气动弹性建模与颤振特性分析

变体飞行器在变体过程中结构质量、刚度和阻尼特性会发生显著变化，导致其气动弹性效应十分复杂。如何高效、准确预测变体过程颤振边界是变体机翼结构动力学设计的难点问题之一。现有的非参数化气动弹性建模方法仅能针对单一变体构型进行颤振分析，对变构型的颤振需重复建模，效率低下且可能存在颤振边界丢失问题。针对后缘连续变弯度的变体机翼颤振分析问题，提出了一种参数化气动弹性建模新方法，并综合非定常气动力、流-固耦合插值和气动弹性建模，对变弯度机翼的参变颤振特性进行系统性分析。为验证该参数化建模方法在预测参变颤振特性的准确性，在变弯度机翼的参变模态特征、气动力计算和颤振预测等方面，进行了数值计算与对比研究。仿真结果表明，该方法可高效准确地预测全参数空间内变体机翼的参变颤振特性。     

状态空间气动弹性模型及其在颤振计算中的应用

本文评述了已有的状态空间气动弹性模型,并对它们提出了一些改进意见。本文还提出了一种新的状态空间气动弹性模型。用两种机翼为例进行了颤振计算。算例表明:这种新的状态空间气动弹性模型是一种精度较高而阶数较低的模型。本文还提出了用状态空间气动弹性模型进行颤振计算时自动识别模态及自动确定颤振点的方法。     

变弯度后缘与常规舵面机翼的颤振主动抑制对比

后缘变弯度机翼的气动弹性建模与稳定性分析日益受到关注。为了探究变弯度后缘相比常规偏转舵面机翼颤振主动抑制的方法与特点，以一个小展弦比后缘变弯度机翼为对象，首先建立结构有限元模型，并引入变弯度后缘变形模态和常规舵面偏转模态，采用亚声速偶极子格网法计算非定常气动力；然后使用基于最小状态法的有理函数拟合进行频域到时域模型的转换，建立两种构型机翼的气动弹性模型，并在建模时考虑了变弯度后缘与常规舵面控制带宽的差异；最后利用线性高斯二次型方法设计控制律进行颤振主动抑制，分析对比两种控制方式的特性差异。结果表明：采用变弯度后缘的闭环系统能够将颤振临界速度提高22%，其提升效果优于常规舵面，所需舵面偏转峰值更小。     

基于CFD/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

二阶时间精度的CFD/CSD耦合算法研究

非线性气动弹性分析中,涉及到非线性流体动力学(CFD)和非线性结构动力学(CSD)的耦合计算问题。本文分析比较了目前几种耦合算法:全耦合、松耦合和紧耦合,并从耦合边界能量传递守恒上就松耦合及紧耦合方法进行了时间精度的分析,得出传统松耦合中即使流体和结构子系统达到高阶时间精度,耦合算法的时间精度仍仅为一阶,紧耦合方法虽然可以达到二阶时间精度,但没有明显提高计算效率。随后,本文基于松耦合流程改进了耦合格式,分析表明其具有二阶时间精度,通过AGARD445.6机翼颤振模型的算例验证,说明该方法可以在保证计算精度的基础上明显提高计算效率,并保持了传统松耦合方法模块化的优点,在模拟非线性气动弹性问题时具有很高的优越性。     

基于CFD技术的高效叶栅耦合颤振分析方法

为了考虑叶轮机叶片结构与流体之间的耦合效应,同时提高叶轮机颤振数值研究的效率,发展了一种基于非定常气动力降阶模型(ROM)的叶栅耦合颤振分析方法。该方法运用时域计算流体力学(CFD)技术计算少数几个叶片的非定常气动力,通过系统辨识及一些假设构建整个叶栅振动的非定常气动力降阶模型,并在状态空间耦合叶栅结构动力学方程建立叶栅气动弹性方程,采用特征值和时域仿真分析该系统稳定性。运用该降阶耦合方法对STCF4(Standard Test Configuration 4)以及NASA Rotor67叶栅系统的稳定性进行了计算。通过与直接计算流体力学/计算结构动力学(CFD/CSD)耦合方法和非耦合方法计算结果的比较验证了该方法的准确性,且该降价耦合方法的计算效率相对于直接CFD/CSD耦合方法提高了1~2个量级,为叶轮机气动弹性参数研究、失谐研究以及多模态耦合计算等提供了便利。     

基于径向基函数的网格变形及非线性气动弹性时域仿真研究

为开展非线性气动弹性研究,基于非线性结构有限元软件NASTRAN和自主研制的多块结构化计算流体力学(CFD)求解器,开发了一套基于计算流体力学/计算结构动力学(CFD/CSD)耦合求解方法的气动弹性时域仿真程序.该程序采用径向基函数(RBF)交换两套求解器之间的数据并进行网格变形.为提高RBF方法的效率,构造了基于多次插值的空间待插值点精简算法.在多次插值过程中,每次插值的对象为上次插值的误差,并同时限制插值区域,以此实现了空间待插值网格数的精简.数个网格变形的算例表明该方法可支持大变形运动,并且具有较高的计算效率.采用此程序开展了AGARD 445.6机翼颤振计算、大展弦比机翼的静气动弹性计算与切尖三角翼极限环振荡(LCO)现象的动气动弹性仿真,结果揭示了当机翼展弦比较大或者响应幅值较大时,结构非线性对于气动弹性有显著影响.     

复合材料大展弦比机翼动力学建模与颤振分析

新一代航空结构广泛采用复合材料，对复合材料机翼的气动弹性工程化建模和分析是飞机设计的重要任务。应用气动弹性分析理论和方法，对复合材料大展弦比机翼进行了结构有限元建模、模型修正、固有振动特性计算、部件发散与颤振工程分析。本文使用MSC／NASTRAN软件，在复合材料大展弦比机翼的初步静力分析模型基础上，依据结构图纸、相关试验结果反复修改得到合理的机翼结构动力学有限元模型，固有振动计算中采用动力减缩方法消除局部模态并提高计算精度，采用亚音速偶极子格网法求解非定常气动力，并对单独机翼进行了发散和颤振计算分析。     

一种基于亚临界响应的颤振稳定性边界预测新方法

基于地面模态试验事先获取的结构模态参数,通过弹性机翼在气流中的亚临界响应,解算出作用在机翼上的非定常模态气动力系数。以弹性机翼的模态位移作为系统输入,模态气动力系数作为系统输出,通过辨识方法获得弹性机翼振动的气动力模型。在时域内耦合结构运动方程和气动力模型,建立基于试验数据的气动弹性稳定性分析模型。通过分析系统稳定性随动压的变化规律,获得弹性机翼的颤振稳定性特性。与经典颤振边界外推方法的主要区别在于该方法实质上只需要一次亚临界响应试验即可预测颤振临界点,可极大降低颤振试验的风险和成本。该方法既可用于颤振风洞试验,也可用于颤振试飞。     

考虑随机型不确定性的浸入式颤振求解方法

气动弹性模型中的参数不确定性一般具有一定的分布规律,为了定量分析随机型参数不确定性对颤振的影响特性,考虑气动弹性系统中广义刚度的随机型不确定性,基于浸入式随机多项式展开(PCE)方法,在传统的颤振求解方法——p-k法的基础上,提出了针对不确定性气动弹性系统稳定性分析的增广p-k法——PCEPK(Polynomial Chaos Expansion with p-k)法,并将该方法应用到某机翼的颤振分析中,研究了均匀分布下的广义刚度不确定性对颤振边界的影响,并同基于结构奇异值μ理论的鲁棒颤振分析的结果和计算效率进行了对比。最后,采用标准的蒙特卡罗模拟(Monte Carlo Simulation,MCS)方法验证了结果的正确性。研究结果表明,PCEPK法计算的颤振边界范围是该分布下的"确定"结果,不因随机样本数而改变,克服了随机方法依赖样本数的缺点。同时,与基于结构奇异值理论的鲁棒颤振分析方法相比,它能够考虑不确定性参数分布对颤振特性的影响,具有更广泛的适用范围。     

基于非定常气动力辨识技术的气动弹性数值模拟

选择离散型输入输出差分模型,运用最小二乘方法进行非定常气动力建模,并将辨识得到的降阶模型用于气动弹性的数值模拟。1个马赫数下的颤振临界点的计算仅需调用一次非定常流场求解器。计算精度保持与非定常欧拉方程计算方法相当的同时计算效率提高了1~2个量级。计算了跨声速具有S型颤振边界的气动弹性标准算例－Isogaiwing和三维气动弹性标模算例AGARD445.6,辨识模型计算边界与非定常Euler方程计算结果吻合。证明非定常气动力辨识技术可以提供高效的高精度的气动弹性分析。