基于欧拉方程的机翼跨音速颤振的频域计算

机翼跨音速颤振的频域计算方法是以给定机翼模态分布下机翼上各点的模态值作为运动幅值，以三维非定常Euler方程为控制方程，采用有限体积法和双时间推进，求解三维机翼简谐运动下的非定常气动力。所求得的气动力作为已知值运用于颤振方程，利用v-g法进行求解。对得到的一系列的阻尼、速度和频率进行了线性插值，从而得到颤振速度和颤振频率。     

基于气动力分段表达的颤振求解方法与模态跟踪技术的应用

在频域颤振方程的求解中,PQI法采用分段气动力的表达形式将颤振方程进行分段直接求解,避免了迭代。由于广义气动力分段插值函数PQI的构建依赖折合频率的选取,将构造的分段插值函数嵌入颤振方程中求解存在模态窜支问题,对此发展了一种基于特征值摄动理论的预测跟踪法来确定各阶模态对应的特征值随参数变化的前后对应关系。文中介绍了PQI颤振分析的基本方法,并采用P-K法进行对比分析来验证算法的计算精度,最后将发展方法应用到某飞机全机轴对称状态颤振分析中来验证算法在实际工程颤振分析中的有效性,数值结果表明该方法计算精度高,模态跟踪效果良好。     

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

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

压电网络对叶盘结构气动阻尼的影响

首先通过引入气动影响系数,建立了振动相关气动力与系统振动位移成比例的线性模型。基于此导出了叶盘结构气弹耦合的线性振动方程。接着在叶盘结构中引入压电网络,建立了流 机 电耦合动力学方程；并采用集总参数等效模型研究了压电网络对叶盘结构气动阻尼的影响,对流 机 电耦合系统的特征值问题进行求解,获得了可作为系统气弹稳定性判据的模态气动阻尼比。研究表明,压电网络可以有效改善叶盘结构主要振动形式——节径型振动的气弹稳定性；通过对压电网络中电学参数进行优化设计,可以在较宽范围内提高叶盘结构的颤振边界；压电网络对弯曲振动模态气动阻尼的改善效果优于弯扭耦合振动模态。      

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

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

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

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

叶轮机气弹耦合方程求解的时间推进方法评估

为了提高叶轮机颤振预测时的计算效率,同时兼顾求解的稳定性,有必要选取最适合的时间推进方法。基于降阶结构 动力学方程的双向流固耦合方法常用于压气机气弹性能分析,针对气动弹性控制方程中结构动力学方程考查了多种时间推进方 法对计算结果的影响。通过单自由度弹簧系统的数值求解探讨了不同时间推进方法的特点;选取实际压气机算例NASA Rotor 67 进行颤振流固耦合分析,进行大量数值试验找到了各时间推进方法需要的最大时间步长,以此为依据对比出不同时间推进方法的 计算效率。结果表明：采用单自由度弹簧系统和压气机颤振耦合求解得出的结论基本一致,说明时间推进方法的性能更多地与方 法本身的数学性质有关,但是流固耦合求解时流场和结构的信息交互会带来一定影响。以实际颤振算例叶片瞬态响应振幅的对 数衰减率作为衡量精度的标准,经典龙格-库塔方法是颤振耦合计算达到相同精度耗时最少的方法,其次是4阶隐式Adams方法 和Newmark方法。     

跨声速颤振计算方法研究

采用改进的弹簧近似方法,以ALE(Arbitrary Lagrangian-Eulerian)描述的Euler方程为控制方程,耦合结构运动方程,计算了Isogai机翼模型的颤振边界。气动力求解采用双时间推进,外时间为物理时间,与结构运动方程同步;结构运动方程求解采用Rayleigh-Ritz方法,时间推进也采用双时间法,节省了大量计算时间。研究表明,计算结果与参考文献提供的结果吻合良好,证明了所提出的跨声速颤振计算方法的准确性。     

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

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

带结构刚度非线性的全动操纵面跨音速颤振特性研究

主要的研究工作是建立了以三维非定常N—S方程为主管方程的计算非定常气动力的高效计算方法，采用颤振方程和N—S方程耦合求解，时间推进的方法计算结构响应特性。发展了一种对既带有结构非线性又含有气动非线性的多重非线性的复杂颤振系统响应特性的分析方法。对全动操纵面转轴带有硬弹簧和间隙型非线性刚度的颤振系统进行了分析，计算结果表明这类系统具有复杂的极限环特性。     

考虑预应力和液体影响的圆柱壳颤振分析

 为了研究预应力和液体对圆柱壳颤振的影响,采用杂交元方法,建立了圆柱壳的气动弹性方程。由Sanders薄壳理论和经典有限元理论,从壳的精确解推导节点位移函数入手,并由一阶活塞理论得到结构的气动阻尼矩阵和刚度矩阵,最终推导出考虑预应力和内部液体影响的圆柱壳的混合有限元公式和气动弹性方程。通过特征值法验证了此种有限元方法的正确性,并重点研究了预应力和内部注有液体对圆柱壳稳定性的影响。计算结果表明:预应力和液体对圆柱壳的颤振特性有显著的影响。     

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

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

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

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

时域法在压气机转子气动弹性计算中的应用

发展了用于三维叶轮机气动弹性计算的软件系统AEAS(aeroelastie analysis system).该系统采用时间推进法求解转子叶片的瞬态位移及动应力,以模拟压气机转子中的气动弹性现象,AEAS主要应用于强迫响应分析和颤振分析,强迫响应分析中,在坎贝尔图上共振点对应的工况下,在转子叶片上施加周期性变化的非定常...     

压气机转子叶片颤振边界的预测方法

通过发展的基于能量法的颤振数值预测方法得到了压气机转子叶片的颤振边界.将计算结构动力学(CSD)分析得到的叶片表面节点位移插值到耦合面的流体网格点上,并将设计的多层动网格技术应用于计算流体动力学(CFD)方法,实现叶栅振荡作用下的非定常分析,得到叶片表面的非定常气动功以及模态气动阻尼比.以某第一级压气机转子叶片为例,对...     

突加高能载荷作用下航空发动机结构动态响应及安全性综述

航空发动机在服役期间可能遭受鸟撞、叶片丢失等突加高能载荷的作用,造成发动机整机/部件动力学特性恶化和关键构件的损伤,危及发动机的结构安全性。本文从突加高能载荷复现方法与传递规律、突加高能载荷作用下转子/整机结构响应研究、突加高能载荷作用下关键构件损伤机理三个方面综述了现有研究工作,并针对近年来发展的抗突加高能载荷的安全性设计方法进行了探讨,最后分析了突加高能载荷问题的科学本质及发展趋势,为突加高能载荷作用下航空发动机安全性设计提供了重要参考。     

Time-domain nonlinear aeroelastic analysis and wind tunnel test of a flexible wing using strain-based beam formulation

A theoretical formulation for time-domain nonlinear aeroelastic analysis of a flexible wing model is presented and validated by wind tunnel tests. A strain-based beam model for nonlinear structural analysis is combined with the Unsteady Vortex Lattice Method (UVLM) to form the complete framework for aeroelastic analysis. The nonlinear second-order differential equations are solved by an implicit time integration scheme that incorporates a Newton-Raphson sub-iteration technique. An advanced fiber optic sensing technique is firstly used in a wind tunnel for measuring large structural deformations. In the theoretical study, the nonlinear flutter boundary is determined by analyzing the transient response about the nonlinear static equilibrium with a series of flow velocities. The gust responses of the wing model at various gust frequencies are also studied. Comparisons of the theoretical and experimental results show that the proposed method is suitable for determining the nonlinear flutter boundary and simulating the gust response of flexible wings in the time domain.     

Structural Nonlinear Flutter Characteristics Analysis for an Actuator-fin System with Dynamic Stiffness

The flutter characteristics of an actuator-fin system are investigated with structural nonlinearity and dynamic stiffness of the electric motor. The component mode substitution method is used to establish the nonlinear governing equations in time domain and frequency domain based on the fundamental dynamic equations of the electric motor and decelerator. The existing describing function method and a proposed iterative method are used to obtain the flutter characteristics containing preload freeplay nonlinea...