一类基于颤振行列式的颤振分析新方法

提出了一类基于颤振行列式的颤振分析新方法,分别采用颤振行列式值的模,或颤振矩阵的模最小特征值这两个指标的频域下界,来直接确定颤振临界点.根据这两个指标在颤振临界点附近取极小值的特点,编制了相应的算法.数值算例验证了这两个指标在确定颤振临界点时的等价性.该方法与传统颤振分析方法的结果精度相当,并且根据颤振临界点处的颤振特征向量还可以定量给出各模态参与颤振耦合的程度.      

基于矩阵奇异值理论的颤振分析新方法

根据颤振分析的基本概念,提出了一类基于矩阵奇异值理论的颤振分析新方法。该方法的特点是,以计算颤振矩阵最小奇异值或条件数的倒数来直接搜索颤振临界点。证明了这两个指标在颤振临界点处的等价性。根据指标在颤振临界点附近取极小值的特点,编制了相应的算法,在确定颤振临界点时无需计算颤振特征根,避免了“窜支”问题,从而减少了人工干预,提高了计算自动化程度。数值算例结果表明,采用该方法计算得到的颤振临界速度和颤振频率与p-k法计算结果的精度相当,且两个指标对应的计算结果一致,验证了其等价性。     

基于谐波平衡法和 V-g 法的高效颤振预测分析

提出了一种高效的颤振预测分析方法,将之应用于跨声速颤振边界分析及结构参数影响研究中。本方法基于频域颤振分析 V-g 方法,为气弹系统提供一定的人工阻尼使之保持简谐运动状态,从而将结构动力学方程转换到频域内。然后通过高效的谐波平衡法得到一系列简谐运动频率下的气动力描述函数矩阵,结合频域结构方程,将气弹系统的稳定性问题转化为广义特征值求解。结果表明：本方法计算得到的颤振边界与高精度的时间推进方法非常吻合,分析效率有了明显的提升,而且当结构参数发生变化后,只需进行若干次广义特征值求解即可得到新的颤振边界,无需像时间推进方法一样开展大量的气动/结构耦合数值模拟。     

基于功能原理的颤振模态参与度分析方法

颤振边界与颤振耦合机理对飞行器设计与颤振试验设计有着重要意义。模态坐标下的颤振计算通过各阶广义坐标的变化特性来确定颤振边界，并对颤振耦合机理进行分析。然而这常常依赖设计人员的工程经验，难以保证判别标准量化统一。基于功能原理提出了一种新的颤振模态参与度分析方法，通过广义坐标下广义气动力做功的能量累积，实现颤振模态参与度分析。分别采用频域方法与CFD/CSD方法对AGARD445.6标模进行颤振计算，验证了所提模态参与度分析方法的正确性。随后针对双体飞机颤振计算中发现的“漂移频率”现象，使用所提方法进行了解释，凸显了方法的优势。综合表明，所提模态参与度分析方法较好地反映了颤振耦合机理，具有指标正确可靠、结果归一化强、物理意义明确、适用于复杂结构复杂模态的特点。     

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

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

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

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

基于高阶谐波平衡的跨声速颤振高效预测方法

跨声速流场激波及其诱导的附面层分离等非线性因素导致跨声速颤振边界很难被准确预测,尤其是目前工程常用的偶极子格网法,在跨声速时该方法的预测精度大幅下降。在雷诺平均Navier-Stokes方程流场求解器的框架内,利用结构模态建立广义结构运动方程,利用径向基函数建立模态振型的插值方法,结合径向基函数和无限插值两种网格变形方法的优点实现高效高鲁棒性网格变形方法,从而实现颤振时间推进分析流程,利用国际颤振标模AGARD445.6机翼验证程序在跨声速颤振边界预测中的可靠性。然而,时域方法在气动/结构反复迭代,需要耗费大量的计算资源和时间。为了提高颤振预测效率,基于高阶谐波平衡（HOHB）方法快速获得广义力影响系数矩阵,利用该矩阵建立频域模态位移和气动力之间的关系,实现高效颤振频域分析方法。通过二维翼型和三维机翼算例进行验证,结果表明：在不对计算精度产生明显影响的前提下,HOHB方法能够提高颤振预测效率约6倍。     

不确定性颤振风险定量分析

 针对小风险颤振定量分析,提出一种概率鲁棒颤振初步分析方法。该方法结合鲁棒颤振分析和概率颤振分析,针对不同层次风险和不确定性水平提供颤振边界特性全面信息,给决策制定提供更好参考。在鲁棒颤振分析中,采用振型迭代的方法,通过结构奇异值( μ )分析和特征值求解得到零风险的鲁棒颤振临界速度。在概率颤振分析中,采用标准蒙特卡罗模拟(MCS)方法得到较大风险时的颤振速度分布和概率临界速度。采用概率比较和二分法估计概率鲁棒颤振不确定度半径和概率鲁棒颤振临界速度,同时采用样本重新利用的随机方法估计颤振风险和不确定性定量关系。通过一个考虑集中质量不确定性的大展弦比双梁式机翼的颤振计算分析表明,确定型μ方法只能用于零风险颤振分析,标准的MCS方法只能用于较大风险颤振速度分布问题。概率鲁棒颤振分析结果表明增加1%颤振风险可以使不确定度半径增大62%或颤振速度边界增大5%,从而放宽了设计要求。     

飞机颤振模态参数的频域子空间辨识

 研究了基于非参数噪声模型的频域子空间系统辨识法，并采用频域子空间辨识算法实现了基于多通道数据的颤振模态参数辨识，改变了传统的单通道颤振试验数据分析模式，试验结果表明子空间算法能有效提取多通道数据中包含的主要危险模态信息，且计算量小，适用于模态参数的在线分析。     

切削颤振在线监测的试验研究

 切削颤振发生过程最显著的频域特征是包含颤振频率的某个窄频段内振动能量的急剧上升。本文根据频段能量原理提出了一种可靠的判据,并利用微型计算机实现了切削颤振的在线监测。该方法具有快速、灵敏、准确和抗干扰能力强等优点。本文还探讨了实际应用中的某些具体措施。     

基于流固耦合的斜激波冲击作用下曲壁板气动弹性分析

采用自主发展的双向流固耦合求解器,研究了斜激波冲击作用下曲壁板的气动弹性响应特性。曲壁板的几何非线性大变形运动方程采用有限差分法求解,流体控制方程基于有限体积法求解,双向流固耦合采用交错迭代算法。计算结果表明:当动压小于临界颤振动压时,曲壁板表现出静平衡状态,且随着动压的增大,壁板变形的非对称性越明显。当动压大于临界颤振动压时,壁板振动位移先增大后减小,最终达到稳定颤振状态,且该极限环颤振并不关于初始位置正负对称的。同时,随着动压的增大,壁板颤振的正向峰谷值、负向峰谷值和振幅均逐渐增大,颤振频率则逐渐减小。壁板振动响应规律并不随着壁板弯曲高度的改变而单调递增或递减,较小的弯曲高度可以降低壁板颤振临界动压值,但是当弯曲高度进一步增大后,由于气动非线性特性增强,准周期无规则运动状态被激发了出来,临界颤振动压迅速升高。     

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

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

Match Point Solution for Robust Flutter Analysis in Constant-Mach Prediction

This paper presents a method for robust flutter computation which uses flight altitude as the perturbation variable in order to obtain a match point solution. The air density and sound speed of standard atmosphere model are approximated as the polynomial function of altitude, such that the flight altitude becomes the single perturbation variable that describes the aeroelastic system. The uncertainties of generalized stiffness and damping are considered and the uncertain aeroelastic system can be formulated as linear fractional transformation （LFT） representation which is suitable for/.t analysis framework. Finally, the match point solutions of robust flutter margins can be computed with structured singular value （SSV） theory. The robust flutter analysis method provided in this paper is suitable for constant-Mach flight flutter test and provides valuable reference for flight envelope expansion.     

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

活塞理论气动力静压修正方法及其在曲壁板颤振分析中的应用

提出了一种采用计算流体力学(CFD)计算的压力分布对活塞理论气动力进行静压修正的方法,将该方法应用到曲壁板的静气动弹性变形及颤振稳定性分析中,并与采用曲率修正活塞理论气动力的计算结果进行了对比.分析结果表明,采用本文提出的活塞理论气动力静压修正方法进行曲壁板的气动弹性分析,在圆柱曲壁板曲率较小的情况下,与采用曲率修正活塞理论气动力方法得到的静气动弹性变形、稳定性边界差别不大;而在曲率较大时,采用本文方法计算得到的曲壁板静气动弹性变形,其曲壁板靠近前缘部分被压的更低,而曲壁板的颤振稳定性边界更小,且这种差别随着圆柱曲壁板曲率的增加而不断增大.该方法突破了曲率修正活塞理论的小曲率限制,扩大了活塞理论气动力在曲壁板颤振分析中的适用范围.     

螺旋颤振分析

 对于装有螺旋桨发动机、涡轮螺桨发动机的飞机,存在着螺桨与发动机组合的螺旋颤振稳定性问题。本文介绍~种实用的分析方法,取沉浮、侧向平动、俯仰、偏航4个自由度,以准定常方法计算出螺旋桨气动力,并计入陀螺力矩,形成系统状态方程,然后提取特征值,以确定系统的稳定性。算例中还对影响稳定性的若干参数进行了分析。     

高空台飞行环境模拟腔μ综合控制设计

针对航空发动机在高空台(AGTF)试验中高空飞行环境模拟腔(FESV)控制问题以及建模过程中产生的未建模动态不确定性问题,提出了一种μ综合控制结构方案,同时对于截止频率相差近几百倍的温度回路和压力回路的大差异特性,提出了采用同一加权灵敏度函数的处理方法,获得了模拟腔控制的期望效果；采用D K迭代算法设计了μ综合控制器,为解决因系统矩阵不满秩带来的温度、压力同时控制的问题提供了一种设计方法和思路,不同的飞行轨迹仿真验证结果表明:μ综合控制器对于飞行环境模拟腔温度和压力的控制效果一致性良好,动态误差不大于09%,稳态误差不大于05%,因此设计的μ综合控制系统具有伺服跟踪性能和鲁棒性能。      

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

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