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

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

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

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

鲁棒颤振裕度法的应用

介绍了一种新的确定颤振边界的方法——鲁棒颤振裕度法。利用结构奇异值理论将颤振理论模型和试飞数据有机结合起来,进行了颤振边界预测。仿真结果表明,该方法是可行的。     

小展弦比压气机转子叶片颤振研究

针对小展弦比压气机转子叶片特点,在能量法计算基础上提出了两种颤振评判方法,即首先利用结构动力学分析结果对叶片频率分离度、振型和共振进行分析,判断某压气机小展弦比转子叶片存在高阶耦合颤振的可能;然后将结构动力学得到的振动位移插值到气动网格上,进行三维非定常黏性流场计算,得到非定常气动功,判断叶片是否存在失速颤振的可能.最后综合利用两部分分析结果和试验及断口分析情况,综合给出该小展弦比压气机叶片因失速颤振和耦合颤振共同作用,导致叶片在短时间内出现了裂纹.      

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

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

模态跟踪技术在颤振计算中的应用

以p-k法为基础，在特征值选取过程中，采用了多点预测方法及优化选择方法进行模态跟踪，有效解决了模态之间的窜支问题。介绍了基于模态跟踪技术的防窜支颤振计算方法；编制了防窜支p-k法颤振计算程序；对某工程实例，分别使用该计算程序与MSC．Nastran软件进行计算，经结果对比分析，证明该方法能有效解决窜支问题。     

鲁棒稳定性分析及其在飞行颤振试验中的应用

结合国内外的研究现状和工程应用需要,基于控制理论利用结构奇异值μ分析法将设计理论模型数据和实测试飞数据结合起来,可以在飞行试验过程中进行颤振边界预测,直接给出保守的颤振临界动压或者颤振临界速度,从而能够增加试验的安全性、缩短试飞周期、减少试验点并节省试验经费,是目前颤振试飞方法的有益补充.     

高超声速全动翼面全时域耦合分析方法及应用

在流场-结构温度场同步计算方法的基础上,建立了多物理场全时域耦合分析方法,将方法应用于沿轨道运动的高超声速全动翼面热气动弹性稳定性分析。采用基于有限体积模型的CFD同步计算方法求解高超声速流场和结构温度场,建立映射关系实现结构有限元模型气动载荷加载和温度场赋值。采用移动坐标系和动网格相结合的方式描述变速度飞行和翼面偏转过程。通过坐标系变换将翼面偏转过程和振动过程的网格变形量叠加,考虑翼面振动和偏转的耦合非定常效应。针对沿轨道运动的高超声速飞行器,建立了同步计算方法与全时域耦合分析方法相结合的热气动弹性稳定性分析流程。研究发现,与同步计算方法相比,全时域耦合分析方法能够模拟结构振动对流场和结构温度场的影响,计算得到的监测点热流密度波动幅值占热流峰值的10%左右,而温度变化并不明显,相比于刚体模型,监测点温度只下降了0.3%左右。全时域耦合热气动弹性分析方法得到的颤振临界点在4-5号状态点之间,颤振形式为铰链扭转模态与一阶弯曲模态的耦合颤振,与"冻结"模态的热颤振方法结果一致。     

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

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

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

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

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.     

不同气流偏角下的壁板热颤振分析及多目标优化设计

 研究了考虑热效应的不同气流偏角下的壁板颤振问题及其多目标优化设计。采用考虑气流偏角影响的一阶活塞气动力、Von-Karmon大变形理论和准定常热应力理论建立了复合材料壁板热颤振方程。利用模拟退火算法,对不同温度场下的偏航壁板颤振速度进行计算。以偏航壁板热颤振速度和壁板重量为多目标函数,在不发生热屈曲的条件下进行优化设计。结果显示：温升使偏航壁板颤振发生“跳跃”现象,对应的气流偏角发生变化;当壁板热颤振模态不变时,偏航壁板颤振速度随温升呈下降趋势,两者呈线性关系;而当热颤振模态发生变化,即偏航壁板颤振发生“跳跃”现象时,偏航壁板颤振速度随温升先升高而后降低,两者呈非线性关系;Pareto解对应的多目标函数之间呈线性关系。     

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

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

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

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

结构非线性颤振半主动控制的理论及实验研究

对结构非线性颤振的半主动控制即“颤振驯化”方案进行了理论及实验研究。以二无非线性颤振系统为对象,分析了半主动控制的原理及带有半主动控制环节后系统的响应特性。设计了相应的颤振半主动控制的风洞实验模型,系统响应信号的监测处理、刚度调节机构的运转控制由以8031芯片为主体的单片机来完成。风洞实验取得了良好的控制效果。     

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

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

高速颤振模型设计中颤振主要模态的判断

颤振模型设计时难以实现完全的动力学相似,需要对颤振主要模态进行合理选择。采用数值分析方法,对颤振模型设计时主要模态的选取问题进行研究。通过各阶模态振型下气动刚度系数的比较、指定运动形式下广义非定常气动力的计算和不同模态截断下颤振结果的收敛特性分析,研究了颤振分析时不同模态运动之间的相互影响,对模态运动引起的气动力和颤振特性变化进行评估。以高超翼面模型为研究对象的数值算例结果表明,几种分析方法所判断的颤振主要模态基本一致。其中基于振型的广义气动刚度系数参数,避免了非定常气动力的计算,可作为颤振模型设计或颤振分析时主要模态选取的快速判断方法。