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

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

车载颤振试验可行性研究

用数值仿真和风洞试验的方法对车载颤振试验进行了可行性验证研究。制作了一个单梁式半机翼弹性模型作为试验模型,建立了其结构与气动力有限元模型。非定常气动力的计算基于偶极子格网法,并由PK法估算其颤振临界速度。在颤振风洞试验和车载颤振试验中,除了用目测颤振现象记录颤振临界速度外,还采集和分析翼尖加速度的亚临界响应信号以推算出颤振临界速度。对比数值仿真、风洞试验和车载试验得到的颤振临界速度可知,三者的结果一致,从而验证了车载颤振试验的可行性。     

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

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

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

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

大攻角机翼的气动弹性计算方法

 本文利用大攻角带分离流机翼的完全非定常非线性气动力计算方法,通过与机翼运动方程的同时求解,在时间域内实现了大攻角机翼非线性气动弹性的数值模拟,根据弹性机翼各种状态下的运动过程,可以得到大攻角机翼的颤振速度等重要参数以及亚临界、超临界等飞行状态的运动规律。算例结果表明,大攻角下机翼的气动弹性问题需引起足够的重视。     

考虑机翼几何非线性的气动弹性建模与分析

大展弦比柔性机翼结构重量轻、气动效率高,广泛应用于高空长航时无人机(UAVs)。飞行过程中,这类机翼在气动力作用下发生大变形,线性结构模型不再适用,需要建立考虑几何大变形的结构模型。采用牛顿力学方法推导了考虑结构几何非线性的机翼结构动力学模型,该方法推导过程简洁、物理意义明确,可以与Hodges基于哈密顿原理的推导方法相互补充,相互验证。为了能够更准确地求解大展弦比柔性机翼的非定常气动力,建立了能够考虑机翼三维效应且适用于机翼空间大变形的非定常气动力模型。基于建立的非线性结构模型和非定常气动力模型,采用松耦合方法建立了非线性气动弹性模型,并通过算例验证了气弹模型的准确性。研究结果表明,大展弦比柔性机翼颤振速度对来流迎角和机翼的展长均较为敏感;当来流速度大于颤振速度时,由于几何非线性,机翼振动并未发散而是形成稳定的极限环振荡(LCO);随着来流速度进一步增加,机翼再次穿过临界稳定点,由不稳定系统变为稳定系统,直到随着速度的增加系统再次达到临界稳定状态。     

热颤振地面模拟试验技术

地面颤振模拟试验是一项以真实飞行器结构作为试验对象，并利用激振器模拟非定常气动力的颤振验证试验技术。本文通过在地面颤振模拟试验的基础上引入热环境模拟设备，进一步研究热颤振地面模拟试验技术。建立了综合考虑多工况的气动插值点优化方法，然后利用Kriging代理模型构建了适用于时变温度场中结构的非定常气动力降阶模型，同时设计了气动加热环境地面模拟及热环境下结构的激励与响应测试方案，最终基于钛合金机翼模型搭建了热颤振地面模拟试验系统，并对时变颤振边界进行跟踪测试。试验结果表明，在激振力控制器的设计控制频带内试验结果与仿真结果吻合较好，但鲁棒控制器较窄的控制带宽限制了热颤振地面模拟试验的适用范围。     

考虑翼吊发动机推力的机翼颤振分析

分析了带翼吊发动机的大展弦比机翼颤振特性,其中考虑了机翼大变形和发动机推力的作用。采用梁单元对双梁式机翼和发动机吊舱进行了结构建模,采用修正的片条理论计算曲面定常气动力,修正的Theodorsen方法计算曲面非定常气动力,发动机推力始终沿发动机轴线。研究机翼在不同攻角时的定常气动力及发动机推力联合作用下固有振动和颤振特性的变化情况。结果表明,发动机推力对机翼静变形和固有频率改变不大,但对颤振临界速度的变化量不容忽视,有可能成为大型运输机机翼设计的限制因素。     

高超声速飞行器受热壁板的气动弹性声振分析

高超声速飞行器壁板在非定常气动力、热载荷和噪声载荷构成的多物理场联合作用下,将表现出复杂的非线性气动弹性声振响应,特别是在颤振临界动压附近,受热载荷以及声载荷作用,壁板表现出复杂的跳变运动。基于von Karman大变形板理论,建立了热-声载荷和气动力共同作用下的壁板运动方程,分析了超声速气流中受热壁板的屈曲变形及热屈曲稳定性,借助势阱概念初步分析了壁板跳变运动产生的机理。通过定义“穿零频次”给出了跳变运动定量的分类方法,并计算得到不同温升和动压情况下,壁板发生跳变运动所对应的临界声压级。结果表明：在颤振临界动压之前,随着动压的增加,受热壁板势阱的深度先增大后减小,且受热壁板的势阱深度随着温升的增加而增大。     

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

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

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

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

基于CFD方法的机翼阵风响应研究

通过引入＂网格速度＂方法模拟阵风条件,求解非定常欧拉方程实现了不同展弦比0012平直机翼阵风响应的数值模拟。首先采用该方法对NACA0006翼型迎角阶跃型阵风的气动力响应进行了计算,计算结果与文献结果、理论值吻合良好。进一步对展弦比分别为5、10的0012平直机翼在迎角阶跃型、One-Minus-Cosine型阵风作用下的气动力响应过程进行了模拟分析。研究结果表明,当展弦比增大时,阵风作用下机翼的升力系数响应会增加,机翼翼根部位气动响应幅值大于翼尖部位响应特性。     

机翼-副翼-调整片的非定常气动力和颤振计算

本文应用偶极子格网法对一个机翼-副翼-调整片颤振模型进行了三元非定常气动力和颤振计算,还研究了网格数目对结果收敛性的影响。计算结果表明本方法具有较好的收敛性,和试验结果相比较表明本方法有一定的工程精度,在飞机设计中能用于操纵面-调整片构型的非定常气动力和颤振计算。     

基于修正偶极子格网法的跨音速颤振分析

针对工程上翼面的跨音速颤振计算存在的问题，阐述了修正偶极子格网法的跨音速颤振计算方法。采用较为精确的定常气动力数据，对工程上常用的偶极子格网法（DLM）的气动力影响系数矩阵（AIC）进行修正，考虑了翼型、厚度及粘性等影响非定常气动力的因素，并将改进后的非定常气动力用于某超临界机翼的跨音速域的颤振计算，结果较好地预测到了跨音速“凹坑”现象。     

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

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

亚跨风洞中舵面亚临界颤振试验

设计了舵面颤振试验装置,在亚跨超风洞中对展弦比2．0的NACA0012矩形舵面开展了颤振试验研究。试验马赫数范围为0．3～0．75。试验采用直接观测法获得舵面在不同质量特性条件下的亚声速和接近跨声速的颤振特性。同时还采用亚临界数据分析方法对试验的扭转应变信号进行了离线分析,即通过采用ARMA方法识别扭转应变信号的阻尼和频率,并通过阻尼外插得到颤振临界动压值。研究结果表明：该试验装置可以用于在现有亚跨超风洞中开展舵面颤振问题研究。当采集的亚临界信号为典型指数衰减信号时,以ARMA方法为基础的亚临界颤振试验技术可以稳定地识别出信号阻尼和频率,并较为准确地获得舵面的颤振临界动压、颤振频率等颤振参数。     

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

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

弹性飞机跨声速机动载荷计算方法

基于跨声速非定常气动力求解的重叠场源法,开展了弹性飞机跨声速机动载荷计算方法研究,为现代飞机结构强度设计提供更加可靠和精确的临界载荷计算方法。首先通过采用重叠场源法求解关于62%根弦俯仰振荡的LANN机翼在马赫数为0.822的非定常气动力并与实验结果进行对比,验证了重叠场源法对跨声速激波效应预测和非定常气动力计算的能力;其次应用频域气动力有理近似技术拟合场源法计算得到的广义气动力系数矩阵,建立了机动载荷分析的状态空间模型;然后完成了某型民用飞机俯仰机动载荷分析,研究了俯仰机动飞行情况下飞机机体状态量及飞机部件载荷响应规律。计算结果表明：考虑机体弹性变形后,机翼和平尾气动载荷响应最大值分别减小了5.1%和10.6%,升降舵气动载荷响应最大值增大了16.2%,在飞机结构强度设计中必须考虑机体弹性效应对飞机部件载荷的影响。     

基于降阶模型的不同厚度翼型颤振边界预测

为了解决CFD/CSD计算效率低的问题,基于CFD技术,构造降阶的非定常气动力模型,并耦合结构运动方程,建立频域/时域气动弹性系统ROM,采用线性自回归滑动平均模型的系统辨识方法,分析了气动弹性系统的标准模型Isogai二维翼型的颤振边界。结果表明,在翼型最大厚度所在位置保持不变时,计算不同翼型厚度下对应的颤振边界得出,随着翼型厚度增加,跨声速凹坑逐渐左移。因此,当翼型最大厚度所在位置保持不变时,为了达到提高颤振速度的目标,通过采用该方法的计算结果来调整机翼翼型厚度,提高机翼对飞行环境的适应能力。     

大展弦比复合材料机翼的非线性颤振分析

大展弦比机翼在气动力作用下产生较大变形,颤振速度和颤振频率随之发生明显变化,线性理论难以获得比较合理的解答。综合考虑了结构几何非线性、气动非线性和材料各向异性对机翼运动状态的影响,将复合材料机翼建模为非线性薄壁单闭室截面梁,建立机翼的运动方程,并使用小扰动分析的方法得到机翼在平衡位置附近的振动方程。采用Theodorsen非定常气动理论构建气动模型,获得机翼在平衡位置附近的非线性颤振方程,并利用v-g法判定机翼颤振稳定性。通过算例演示了一些非线性颤振的特点,讨论了铺层角、展弦比、机翼线密度等参数对颤振速度的影响,并与线性理论得到的结果进行对比。