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

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

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

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

具有随机不确定性的机翼颤振优化

针对不确定性参数在颤振优化分析中的影响,提出了一种含颤振约束的、具有随机不确定性的机翼气动弹性优化设计方法.在MSC.NASTRAN平台上建立了结构的动力学模型、采用ZAERO程序进行了非定常气动力计算与颤振分析,用响应面方法计算了隐式失效功能函数,以重量最小为设计目标,对包含颤振概率约束的复合材料机翼进行了优化设计....     

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

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

含不确定参数的复合材料壁板热颤振分析

研究了含有不确定结构参数的壁板颤振问题,利用vonKarman大变形应变-位移关系、气动力活塞理论和准定常热应力理论建立了复合材料壁板颤振的气动弹性力学模型,考虑在壁板颤振分析模型中存在的不确定参数,将其用区间向量定量化,基于区间扩张理论和Taylor级数展开,并结合有限元计算方法,提出了区间分析的方法来估计含有不确定参数的壁板结构颤振临界风速的区间,以及发生极限环振动时振幅的变化区间。通过数值算例,将本文提出的壁板颤振的区间有限元模型与随机有限元模型进行了比较,显示了本文方法的有效性和可行性。这种方法的优点是只需要知道不确定参数的所在范围界限,为解决含有不确定参数的壁板颤振这类复杂的气动弹性动力学问题提供了一个途径。     

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

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

超声速气流中受热壁板的二次失稳型颤振

研究了超声速气流中受热壁板的非线性气动弹性响应,发现了一种新的动态失稳现象——二次失稳型颤振。基于von Karman非线性应变-位移关系、Reissner-Mindlin板理论和一阶活塞理论建立超声速气流中三维壁板的有限元模型。通过数值算例,研究了超声速气流中受热壁板发生二次失稳型颤振的条件,并运用非线性振动理论分析了二次失稳型颤振的机理。研究表明,超声速气流中受热壁板在平衡态的稳定性未发生变化时,也会因系统参数的变化引起气动弹性响应性质的突变,导致壁板的二次失稳型颤振。二次失稳型颤振能否发生不仅受到气流速压和壁板温升的影响,而且还与初始扰动有关。当扰动引起壁板的初始变形较小时,不能激发出二次失稳型颤振,壁板的气动弹性响应最终收敛到屈曲平衡态。应用二次失稳型颤振理论和分析方法,确定了前人给出的一个金属壁板模型的热颤振边界的风洞试验结果,而且计算结果与试验结果符合良好,从而对这一壁板热颤振现象的风洞试验结果作出了较合理的理论解释。     

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

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

有限长壁面处理对叶片气动弹性稳定性的影响

应用壁面处理控制叶片颤振与其他方法相比有其特殊优势,但由于加入软壁面特别是有限长的壁面处理,给求解方法带来很多困难.本文提出一种新的模型,采用一种新方法--传递单元法,从而有效避免了界面匹配和特征值求解中引起的各种困难问题,并且通过该模型可以探讨研究者们更为关心的有限长壁面处理对叶片颤振特性的影响.数值计算结果表明,有限长壁面处理阻抗值的大小对气动弹性稳定性有显著影响.因此,本文提出的模型为实现应用壁面处理控制叶片颤振,进一步丰富了其理论依据.     

基于扫频激励的高安全性颤振边界预测方法

针对一种基于脉冲激励的高安全性颤振边界预测方法,提出改变其原有的脉冲激励形式,采用扫频激励以改善脉冲激励激振频带较窄和信噪比较低等问题。根据采样数据和气动弹性运动方程,获得亚临界速度下的模态位移响应和模态气动力系数。采用系统辨识构建相应的非定常气动力模型,耦合状态空间的结构运动方程和气动力方程得到气动弹性稳定性分析状态方程,进而预测颤振边界。基于两种激励形式预测的颤振边界的对比结果表明,在较低的测试速度下采用扫频激励具有更高的预测精度。     

一种建立非定常气动力频域模型的简单方法

本文给出的方法,线化了非定常气动力的karpel近似式的非线性部分,并吸取了Roger近似式中同一矩阵各元素独立确定以得到较高精度的特点。它不仅简化了拟合中的计算,改善了拟合精度,还保持了Karpel方法所引入的最小阶次。     

基于NN-PSM的航空发动机机载自适应稳态模型

为建立一种适用于大包线、变状态的高精度、高实时性航空发动机机载自适应稳态模型，提出一种基于神经网络和推进系统矩阵相融合（NN-PSM）的机载自适应稳态模型建模方法。该方法基于小偏差线性化方法对发动机进行线性化来提取推进系统矩阵，用于表征机载模型与发动机之间的输出偏差量。基于神经网络建立发动机基线模型，用于映射飞行条件与发动机输出量之间的关系，利用神经网络的强拟合能力提高机载模型的稳态精度；设计卡尔曼滤波器实时估计发动机健康参数，提高模型的自适应能力。在大包线、变状态的飞行条件下进行仿真验证，并与传统的复合推进系统模型（CPSM）进行对比，结果表明:NN-PSM模型的平均精度在0.66%以内，而CPSM的平均精度为2.07%以内，运行时间仅为CPSM的1/10，且具有数据存储量少的特点。      

某型飞机颤振试飞数据的模态参数识别

 地面模态试验表明某机型的机翼对称一阶弯曲与发动机对称俯仰频率比较接近。为了检查这两阶模态是否会在设计的飞行速度下引起颤振,在不同飞行速度下利用脉冲激励方式测量了频响函数。这组频响函数不仅在关心频段存在密集模态,而且受噪声干扰严重。采用了修正的整体正交多项式识别方法(引入额外的加权函数来降低估计的方差,分子和分母引入采用相同的正交基) , 对这组数据进行了模态分析,并利用一些先验信息和规则从拟合模态鉴别出真实模态。研究表明机翼对称一弯和发动机对称俯仰的两阶模态的频率随速度增加而趋近,两阶阻尼比均随速度增加而下降,但是在设计速度下,模态参数离颤振有较大裕量。     

用逐步扩阶双递推时域法识别线性振动系统复模态参数

本文给出了从线性振动系统的自由响应识别该系统复模态参数的方法,把自由响应的表达式变为一个自回归方程和一个多项式方程,利用最小二乘递推与逐步扩阶递推的双递推法来进行参数识别。文中讨论了几个应用中的问题,并给出了计算机的模拟计算结果以及越野汽车车架的参数识别试验结果。     

Dryden大气紊流频谱的等价确定性函数及其在飞机响应问题中的应用

本文扼要评述了B.Etkin等人最近发表的多输入随机过程的等价确定性输入方法（EDIT）,并且加以扩展,导出线性变换公式,分析了确定性输入函数的奇偶性。按照EDIT把目前常用的Dryden大气紊流频谱化成等价确定性时间函数,并拟合成简单的表达式。说明了在时域内计算飞机对大气紊流响应特性的方法和步骤。在算例中比较了用EDIT和常规的功率频谱法所得的结果。     

采用乘性速压摄动的频域颤振预测方法

 在鲁棒颤振分析方法基础上,提出一种乘性速压摄动表达式,使得颤振系统成为不确定性系统,并直接应用简谐气动力建立其μ框架,从而可以通过频域μ分析预测颤振临界点。导出了判定颤振临界点的μ极值条件,并给出其算法实现。该方法属于频域方法,在理论上具有良好的一致性;且无需求解颤振特征值,从根本上避免了经典的颤振分析方法中可能出现的“窜支”问题。计算结果表明,该方法与p-k法所得结果精度相当,适用于工程颤振分析。     

Flutter analysis of an airfoil with bounded random parameters in incompressible flow via Gegenbauer polynomial approximation

The effects of parameter uncertainty on the flutter characteristics of a two-dimensional airfoil in an incompressible flow were investigated through Gegenbauer polynomial approximation. The uncertain parameters, such as the linear and cubic pitch stiffness coefficients are modeled as bounded random variables with λ-PDFs (probability density functions). With the aid of Gegenbauer polynomial approximation, the two-dimensional stochastic airfoil system is transformed at first into its equivalent deterministic one. Then the Hopf-bifurcation point is determined through the equivalent deterministic system, and the onset of the flutter, together with the flutter frequency against the probability density distribution parameter and the intensity of the random variable is explored. In addition, the ranges of the peak pitch and plunge responses against the flight speed and the PDFs of the peak pitch response are obtained. Numerical results show how the probability density distribution parameters and the intensities of random parameters affect the flutter onset speed, flutter angular frequency, the upper and lower boundaries of peak responses, and the PDFs of peak responses.     

基于自然激励技术的颤振边界预测

为了预测紊流激励条件下机翼的颤振边界,基于自然激励技术提取紊流响应的自由衰减信号,采用矩阵束方法识别模态参数,最后通过Z-W(Zimmerman-Weissenburger)方法计算稳定性判据,拟合判据变化曲线并外推颤振边界.对平板机翼模型进行了数值仿真分析,对单独机翼模型风洞颤振试验数据进行了计算.结果表明:采用自然激励技术与矩阵束方法能够较准确地识别紊流激励响应的模态参数,频率识别误差小于6％,阻尼比识别误差小于30％,结合Z-W方法能够在较低风速较早地预测颤振边界,有助于提高试验的安全性.     

状态变量模型建立方法

为了克服小扰动方法精度不高、稳态终值响应法动态过程不一致以及传统拟合法的随着维数增加精度下降和拟合时间长等缺点,综合了3种算法（小扰动法、稳态终值响应法和传统拟合法）的优点,设计了基于变尺度法的混合求解方法建立了航空发动机状态变量线性模型。仿真结果表明：建立的航空发动机状态变量线性模型与非线性部件级模型在动态过程响应中吻合良好,而且具有较高的稳态精度,能够保证航空发动机最终稳定值的一致性。     

嫦娥三号动力落月段轨迹确定策略

针对CE-3（嫦娥三号）月球探测器动力下降弧段,特别是悬停避障段频繁机动的特点,提出了采用B样条函数逼近方法进行落月轨迹确定.仿真分析表明：在动力下降运动较平滑弧段,B样条逼近法计算结果略优于多项式拟合法;而在频繁机动弧段,B样条逼近法有明显优势.计算结果表明,加入VLBI（Very Long Baseline Interferometry,甚长基线干涉测量）数据后能有效提高落月轨迹确定精度,在没有系统误差的情况下联合定位后位置精度优于50 m.此外,还分析了三向测量系统差对定位的影响,可对CE 3任务提供参考.最后对CE 3实测数据进行处理,动力落月段末点位置和着陆器定位计算值相差不到200 m.