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

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

高速飞行器壁板颤振分析的研究进展

壁板颤振是壁板结构在高速气流中发生的一种自激振动现象，特别是超音速和高超音速飞行器上特别容易发生这种现象。壁板颤振引发的非线性振动将对高速飞行器结构的疲劳强度、飞行性能和飞行安全带来不利的影响。随着高速飞行器研发工作的开展，壁板颤振问题将得到越来越多的重视。根据目前国内外高速飞行器壁板颤振的研究现状，介绍了壁板颤振的六种分析模型并从结构理论和气动力理论出发详述了这种分类的依据。阐述了温度、气流偏角、壁板几何尺寸及边界条件对壁板颤振的影响规律。并介绍了目前用于分析壁板颤振问题的频域和时域方法并总结了各种分析方法的优缺点。最后归纳了目前对高速飞行器壁板颤振研究得出的几个重要结论，提出了今后在高速飞行器壁板颤振研究中需要解决的若干问题。     

超声速流动中功能梯度曲壁板的热气动弹性颤振机理

对高超声速环境中功能梯度曲壁板的热气动弹性颤振机理及分岔特性进行了研究。分别采用活塞理论和Eckert参考焓方法模拟气动力以及气动加热效应,在求解板内二维热传导方程以及考虑温升对材料物性影响的基础上,建立了一个气动加热-气动弹性双向耦合的功能梯度曲壁板的热气动弹性颤振模型。采用有限元方法对曲壁板控制方程进行了数值模拟,着重分析了不同拱高下曲壁板的分岔行为,探讨了拱高变化对曲壁板分岔图的影响,发现了曲壁板颤振三种典型的颤振行为,即:热屈曲、混沌以及规则振动。对初始拱高板厚比为1时,曲壁板的两种规则振动行为进行对比发现,随着马赫数的增大,气动加热效应所引起的热内力会使曲壁板的规则振动更加复杂,同时振动的主振型及频率均会发生变化。     

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

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

超音速气流中受热壁板的非线性颤振特性

随着高速飞行器技术的快速发展,壁板热颤振成为国内外研究人员的关注热点.壁板热颤振对飞行器性能有重大影响,甚至影响飞行安全.以超音速气流下的无限展长二维壁板结构作为研究对象,计入热效应的影响,根据Kirchhoff平板理论和Von Karman大变形几何非线性壁板理论建立系统的运动微分方程.并以壁板在x =0.25处为算例绘制了其前三阶弯曲构型的静态分岔图,对不同轴向载荷条件下壁板的屈曲构型进行了分析.考虑到温差△T是影响热应力的重要因素,对比了不同温差条件下,壁板的静气动弹性变形图,结果表明,温差越大,壁板偏离静平衡位置的位移越大.     

基于流固耦合的激波/边界层干扰作用下壁板颤振特性

为了解激波/边界层干扰作用下壁板气动弹性及其对流动分离的影响，采用自主开发的双向流固耦合求解器，对不同激波冲击位置下壁板的振动响应和流动特性进行了数值模拟研究。壁板几何非线性运动方程采用有限差分法求解，基于有限体积法求解Navier-Stokes方程组，对流通量采用MUSCL和AUSMPW+格式离散，双向流固耦合采用交错迭代算法。研究结果表明：激波/边界层干扰作用下壁板振动位移先增大后减小，经若干振荡周期后达到稳定颤振状态，呈现二阶振动模态，壁板变形相对于激波冲击位置呈现非对称性，壁板前部分的振幅始终小于壁板后部分；激波冲击位置可显著改变壁板的颤振振幅、频率及分离区长度，当激波冲击位置靠近壁板两端时，壁板振动最终收敛达到静稳定状态；壁板振动响应与流场特征不随激波冲击位置的改变而单调变化，对于激波冲击位置x/a=0.35工况，壁板颤振可有效抑制激波/边界层干扰流动分离。     

壁板的非线性热稳定分析

<正>壁板热颤振~([1])一般指在超声速气流中,由于惯性力、弹性力、流经壁板表面气流引发的气动力和摩擦产生的热应力耦合作用下壁板产生的一种自激振动现象。在实际研究中发现,当受热壁板发生颤振时,壁板横向振动的振幅大小与壁板的厚度数量级相当,此时属于大变形,壁板的面内张力与壁板的横向位移出现了一种非线性的关系。过去几十年,对于受热板的非线性振动分析受到了很     

非线性壁板颤振分析

利用一种改进的计算流体力学与计算结构动力学(CFD/CSD)耦合方法研究了由气动和结构几何非线性引起的壁板颤振问题。在非定常气动力计算中,考虑了通量分裂格式、隐式时间推进方法和几何守恒律;二维和三维壁板的结构几何非线性建模则采用了有限元的协同旋转理论,并利用一种近似能量守恒算法求解结构的非线性响应。流场和结构求解器采用二阶松耦合方法联立求解,并将其应用于壁板在超声速、跨声速和亚声速的颤振计算中。当考虑结构几何非线性和气动非线性时,出现了典型的极限环振荡现象,并对颤振边界和极限环振荡幅度进行了对比分析。     

高温环境下二维正交编织C/SiC复合材料壁板振动模态演化

为了研究高温环境下二维正交编织C/SiC复合材料壁板的固有振动特性随温度的变化规律,进一步揭示模态跃迁和丢失现象。首先建立了二维正交编织复合材料的细观非均匀有限元模型,基于细观模型采用体积平均法计算得到了均匀化后材料的模量、热膨胀系数和热传导系数等宏观性能参数。在此基础上,分别研究了均匀温度和线性非均匀温度载荷下,四边简支复合材料板在屈曲前后固有频率及模态振型随温度的变化规律。研究结果表明:在均匀温度场下,如果仅考虑热载荷对材料模量的影响时,随着温度的升高,各阶固有频率逐渐降低,但降低的幅度不大;仅考虑热应力对结构刚度影响时,在临界屈曲温度后固有频率反而上升,并且随着温度的升高出现了频率交错,其对应振型发生了跃迁现象;当同时考虑热物性和热应力影响时,各阶频率值变化趋势与仅考虑热应力时类似,只是对应温度点的频率值有所下降。而在线性非均匀温度场下,通过对比不同工况下的模态振型发现,随着温度的升高某些模态振型还会出现丢失,并且这种丢失是从低阶到高阶依次发生的。     

考虑气动热的机翼颤振特性分析

为研究气动热对高超音速飞行器气动弹性特性的影响,首先对小展弦比机翼施加不同温度场,通过热应力分析理论获得不同温度载荷下机翼的热模态,然后应用二阶活塞理论计算非定常气动力,最后采用p-k法进行颤振计算。计算结果表明,气动加热后结构的模态特性和颤振特性均发生变化,由于温度效应降低了各阶固有频率,并且改变了它们之间的差距,从而导致机翼颤振速度降低。     

高速流场中变刚度复合材料层合板颤振分析

变刚度复合材料层合板在高速流场中的颤振行为是设计中需要考虑的问题。本文研究了高速流场中的曲线纤维变刚度层合复合材料壁板非线性颤振响应,分析了边界条件和纤维方向对颤振特性的影响。利用von-Karman大变形应变-位移关系,采用气动力活塞理论,根据虚功原理和有限元法建立变刚度复合材料壁板颤振的气动弹性力学模型,采用Newmark法对壁板的颤振方程求解。给出了不同边界条件和纤维方向条件下层合复合材料壁板的颤振特性。计算结果表明：随着纤维在板中心处或在边界±a/2处与x方向夹角（T₀或T₁）的增大,颤振临界动压减小;相同动压下,随着T₀或T₁的增大,极限环振幅增大。研究表明采用曲线纤维进一步提高了复合材料层合板的可设计性,通过调整曲线纤维路径可以改变复合材料壁板的颤振特性。     

传热时间迟滞影响的薄板热气动弹性耦合 振荡模型

为了对高超声速下飞行器薄壁结构的振动行为进行研究,建立了一种考虑传热时间迟滞影响的二维无限长薄板的热气动弹性耦合振荡模型.同时,给出了只考虑有限历史时间影响的板内温度分布表达式,并对提出的温度表达式的精确性进行了数学证明,从而将以往被忽略掉的热应力产生的板内力矩引入薄板振动方程.数值模拟结果表明:气动力是薄板振动的主要驱动力,板内力的作用相当于外部预紧力,而板内力矩的变化会驱使薄板进行小振幅振动.综合以上三种因素,对薄板在这三种因素耦合情况下的振动行为进行了研究,发现在来流马赫数较大时,薄板会进入十分复杂的振动状态.最后,通过非线性动力学理论对薄板的振动行为进行了分析,发现了薄板振动中的分岔、混沌等现象,以及通向该类现象的途径.      

机翼对螺旋桨发动机旋转颤振的影响研究

为了研究机翼对螺旋桨发动机旋转颤振的影响规律,揭示其影响机理,通过片条理论计算螺旋桨气动载荷,在MSC.NASTRAN软件平台上进行二次开发,对不同结构参数的吊舱-螺旋桨系统和机翼-吊舱-螺旋桨系统进行了旋转颤振分析。研究表明,机翼结构的弹性会大幅降低发动机俯仰模态频率,而小幅降低发动机偏航模态频率,从而改变两模态频率之差,影响模态耦合的程度,进而改变旋转颤振速度。另外,当发动机的运动与机翼翼面的运动耦合紧密时,机翼翼面的气动载荷能够显著提高发动机的旋转颤振速度。     

Nonlinear aeroelastic analysis of the folding fin with freeplay under thermal environment

The nonlinear aeroelastic behavior of a folding fin in supersonic flow is investigated in this paper. The finite element model of the fin is established and the deployable hinges are represented by three torsion springs with the freeplay nonlinearity. The aerodynamic grid point is assumed to be at the center of each aerodynamic box for simplicity. The aerodynamic governing equation is given by using the infinite plate spline method and the modified linear piston theory. An improved fixed-interface modal synthesis method, which can reduce the rigid connections at the interface, is developed to save the problem size and computation time. The uniform temperature field is applied to create the thermal environment. For the linear flutter analyses, the flutter speed increases first and then decreases with the rise of the hinge stiffness due to the change of the flutter coupling mechanism. For the nonlinear analyses, a larger freeplay angle results in a higher vibration divergent speed. Two different types of limit cycle oscillations and a multiperiodic motion are observed in the wide range of airspeed under the linear flutter boundary. The linear flutter speed shows a slight descend in the thermal environment, but the effect of the temperature on the vibration divergent speed is different under different hinge stiffnesses when there exists freeplay.     

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

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

Aeroelastic stability analysis of heated flexible panel subjected to an oblique shock

The critical conditions for aeroelastic stability and the stability boundaries of a flexible two-dimensional heated panel subjected to an impinging oblique shock are considered using theoretical analysis and numerical computations, respectively. The von-Karman large deflection theory of isotropic flat plates is used to account for the geometrical nonlinearity of the heated panel, and local first-order piston theory is employed in the region before and after shock waves to estimate the aerodynamic pressure. The coupled partial differential governing equations, according to the Hamilton principle, are established with thermal effect based on quasi-steady thermal stress theory. The Galerkin discrete method is employed to truncate the partial differential equations into a set of ordinary differential equations, which are then solved by the fourth-order Runge-Kutta numerical integration method. Lyapunov indirect method is applied to evaluate the stability of the heated panel. The results show that a new aeroelastic instability (distinct from regular panel flutter) arises from the complex interaction of the incident and reflected wave system with the panel flexural modes and thermal loads. What’s more, stability of the panel is reduced in the presence of the oblique shock. In other words, the heated panel becomes aeroelastically unstable at relatively small flight aerodynamic pressure.