Aeroelastic effect on aerothermoacoustic response of metallic panels in supersonic flow

A finite element formulation is presented for the analysis of the aeroelastic effect on the aerothermoacoustic response of metallic panels in supersonic flow. The first-order shear deforma-tion theory (FSDT) and the von Karman nonlinear strain-displacement relationships are employed to consider the geometric nonlinearity induced by large deflections. The piston theory and the Gaus-sian white noise are used to simulate the mean flow aerodynamics and the turbulence from the boundary layer. The thermal loading is assumed to be steady and uniformly distributed, and the material properties are assumed to be temperature independent. The governing equations of motion are firstly formulated in structural node degrees of freedom by using the principle of virtual work, and then transformed and reduced to a set of coupled nonlinear Duffing oscillators in modal coor-dinates. The dynamic response of a panel is obtained by the Runge-Kutta integration method. The results indicate that the increasing aeroelastic effect can lead the panel vibration from a random motion to a highly ordered motion in the fashion of diffused limit cycle oscillations (LCOs), and remarkably alter the stochastic bifurcation and the spectrum of the aerothermoacoustic response. On the other hand there exists a counterbalance mechanism between the external random loading and the aeroelastic effect, which mainly functions through the nonlinear frequency-amplitude response. It is surmised that the aeroelastic effect must be considered in sonic fatigue analysis for panel structures in supersonic flow.     

Response regimes of narrow-band stochastic excited linear oscillator coupled to nonlinear energy sink

This paper draws attention to the issue of the vibration absorption of nonlinear mechanical system coupled to nonlinear energy sink(NES) under the impact of the narrow band stochastic excitation. Firstly, based on the complex-averaging method and frequency detuning methodology,response regimes of oscillators have been researched under the linear impact of coupling a nonlinear attachment with less relativistic mass and an external sinusoidal forcing, of which results turn out that the quasi-periodicity response regime of system which occurs when the external excitation amplitude exceeds the critical values will be the precondition of the targeted energy transfer.Secondly, basing on the path integration method, vibration suppression of NES has been researched when it is affected by a main oscillator with a narrow band stochastic force in the form of trigonometric functions, of which results show that response regimes are affected by the amplitude of stochastic excitation and the disturbance strength. Finally, all these conclusions have been approved by the numerical verification and coincided with the theoretical analysis; meanwhile, after the comparing analysis with the optimal linear absorber, it turns out that the NES which is affected by the narrow band stochastic force could also suppress the vibration of system with a better effect.     

Supersonic flutter control and optimization of metamaterial plate

A metamaterial plate is designed by embedding a periodic array of local nonlinear resonators for its supersonic flutter control. Based on the von Karman large deformation theory and supersonic piston aerodynamic theory, the nonlinear aeroelastic equations of the metamaterial plate are obtained by using the Hamilton principle. The comparisons for aeroelastic behaviors of the metamaterial plate and pure plate show that the proposed metamaterial plate can lead to an enlarged flutter boundary and lower vibration amplitude. Furthermore, a parametric optimization strategy for local nonlinear resonators is proposed to improve the nonlinear flutter behaviors of the metamaterial plate, and a significant enhancement of passive control performance can be achieved through optimization design. The present study demonstrates that the design of the metamaterial plate can provide an effective approach and potential application for nonlinear flutter suppression of supersonic plate.     

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

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

两自由度非线性吸振器在简谐激励下的振动抑制

 为了分析两自由度(2-DOF)非线性吸振器(NES)在主共振附近产生周期及非周期响应时的振动抑制效果,首先应用复变量-平均法,推导了线性主结构连接一个2-DOF NES的系统在主共振附近的慢变力学模型,并对系统平衡点的稳定性进行了分析,进一步对系统响应类型进行了预测。然后,应用增量谐波平衡(IHB)法,对系统的慢变模型进行了数值验证。之后,对系统在稳态响应段及强调制响应(SMR)段的振动抑制效果进行了比较。最后研究了该系统在激励频率小幅变化时的振动抑制效果,并同单自由度(SDOF)NES进行了比较。仿真结果表明,2-DOF NES在强调制响应段的振动抑制效果更好,并且得到了比同质量优化后的S-DOF NES更好的振动抑制效果。     

Recent advance in nonlinear aeroelastic analysis and control of the aircraft

A review on the recent advance in nonlinear aeroelasticity of the aircraft is presented in this paper. The nonlinear aeroelastic problems are divided into three types based on different research objects, namely the two dimensional airfoil, the wing, and the full aircraft. Different non- linearities encountered in aeroelastic systems are discussed firstly, where the emphases is placed on new nonlinear model to describe tested nonlinear relationship. Research techniques, especially new theoretical methods and aeroelastic flutter control methods are investigated in detail. The route to chaos and the cause of chaotic motion of two-dimensional aeroelastic system are summarized. Var- ious structural modeling methods for the high-aspect-ratio wing with geometric nonlinearity are dis- cussed. Accordingly, aerodynamic modeling approaches have been developed for the aeroelastic modeling of nonlinear high-aspect-ratio wings. Nonlinear aeroelasticity about high-altitude long- endurance （HALE） and fight aircrafts are studied separately. Finally, conclusions and the chal- lenges of the development in nonlinear aeroelasticity are concluded. Nonlinear aeroelastic problems of morphing wing, energy harvesting, and flapping aircrafts are proposed as new directions in the future.     

Nonlinear thermomechanical deformation behaviour of P-FGM shallow spherical shell panel

In the present article, the linear and the nonlinear deformation behaviour of functionally graded (FG) spherical shell panel are examined under thermomechanical load. The temperature-dependent effective material properties of FG shell panel are evaluated using Voigt’s micro-mechanical rule in conjunction with power-law distribution. The nonlinear mathematical model of the FG shell panel is developed based on higher-order shear deformation theory and Green-Lagrange type geometrical nonlinearity. The desired nonlinear governing equation of the FG shell panel is computed using the variational principle. The model is discretised through suitable nonlinear finite element steps and solved using direct iterative method. The convergence and the val-idation behaviour of the present numerical model are performed to show the efficacy of the model. The effect of different parameters on the nonlinear deformation behaviour of FG spherical shell panel is highlighted by solving numerous examples.     

NES对二维机翼气弹不稳定性的抑制作用

采用数值方法研究了加装非线性能量汇（NES）的二维机翼在不同速度来流下的振动响应机制,着重探索了NES对系统振动的抑制以及系统内的靶向能量传递（TET）特性。首先,建立了加装在机翼前缘及后缘的NES与二维机翼的耦合系统模型,该模型考虑了机翼的沉浮与扭转振动。然后,从非线性振动响应和能量传递等几个方面研究了前NES与后NES对机翼振动的抑制效果和机制。进一步,应用频谱分析发现了此非线性耦合系统振动中存在共振捕获（resonance captures）特性,同时研究了前、后NES与机翼振动模式（沉浮与俯仰）间的靶向能量传递现象与机翼不同的极限环运动之间的对应关系。结果表明,采用前、后都加装NES的方法能够拓宽NES与机翼振动模式间发生靶向能量传递与共振捕获的频率范围,从而提升NES对机翼振动进行有效抑制的临界来流速度。     

大展弦比机翼的非线性气弹响应分析

高空长航时飞行器将在未来军事与民用领域具有广阔的应用前景,这类飞行器具有大展弦比、大柔度机翼,在对该类型结构的气弹响应分析中,必须要考虑机翼的非线性因素.目前,对大型的工程结构MSC.NASTRAN提供了线性颤振、线性气弹的计算功能,但无法进行非线性结构的气弹响应分析.本文在非线性运动方程中引入ONERA气动力的线性项,利用DMAP语言在NASTRAN非线性瞬态响应分析程序中加入相关气动力模块,实现了大型有限元机翼的非线性气弹响应分析.文中给出了算例,结果验证了该法分析大展弦比机翼结构气弹响应的有效性,具有重要的工程应用价值.     

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.     

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.     

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

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

LPV modeling and controller design for body freedom flutter suppression subject to actuator saturation

In recent years, the Active Flutter Suppression (AFS) employing Linear Parameter-Varying (LPV) framework has become a hot spot in the research field. Nevertheless, the flutter suppression technique is facing two severe challenges. On the one hand, due to the fatal risk of flight test near critical airspeed, it is hard to obtain the accurate mathematical model of the aeroelastic system from the testing data. On the other hand, saturation of the actuator may degrade the closed-loop performance, which was often neglected in the past work. To tackle these two problems, a new active controller design procedure is proposed to suppress flutter in this paper. Firstly, with the aid of LPV model order reduction method and State-space Model Interpolation of Local Estimates (SMILE) technique, a set of high-fidelity Linear Time-Invariant (LTI) models which are usually derived from flight tests at different subcritical airspeeds are reduced and interpolated to construct an LPV model of an aeroelastic system. And then, the unstable aeroelastic dynamics beyond critical airspeed are ‘predicted’ by extrapolating the resulting LPV model. Secondly, based on the control-oriented LPV model, an AFS controller in LPV framework which is composed of a nominal LPV controller and an LPV anti-windup compensator is designed to suppress the aeroelastic vibration and overcome the performance degradation caused by actuator saturation. Although the nominal LPV controller may have superior performance in linear simulation in which the saturation effect is ignored, the results of the numerical simulations show that the nominal LPV controller fails to suppress the Body Freedom Flutter (BFF) when encountering the actuator saturation. However, the LPV anti-windup compensator not only enhances the nominal controller’s performance but also helps the nominal controller to stabilize the unstable aeroelastic system when encountering serious actuator saturation.     

Reliability analysis of aero-engine blades considering nonlinear strength degeneration

To comprehensively consider the effects of strength degeneration and failure correlation, an improved stress-strength interference (SSI) model is proposed to analyze the reliability of aeroengine blades with the fatigue failure mode. Two types of TC4 alloy experiments are conducted for the study on the damage accumulation law. All the parameters in the nonlinear damage model are obtained by the tension-compression fatigue tests, and the accuracy of the nonlinear damage model is verified by the damage tests. The strength degeneration model is put forward on the basis of the Chaboche nonlinear damage theory and the Griffith fracture criterion, and determined by measuring the fatigue toughness during the tests. From the comparison of two kinds of degeneration models based on the Miner’s linear law and the nonlinear damage model respectively, the nonlinear model has a significant advantage on prediction accuracy especially in the later period of life. A time-dependent SSI reliability model is established. By computing the stress distribution using the finite element (FE) technique, the reliability of a single blade during the whole service life is obtained. Considering the failure correlation of components, a modified reliability model of aero-engine blades with common cause failure (CCF) is presented. It shows a closer and more reasonable process with the actual working condition. The improved reliability model is illustrated to be applied to aero-engine blades well, and the approach purposed in this paper is suitable for any actual machinery component of aero-engine rotor systems.     

一种CFD/CSD耦合计算方法

针对柔性大展弦比机翼气动弹性分析和主动弹性机翼(AAW)设计发展了一种计算流体动力学(CFD)和计算结构动力学(CSD)的耦合计算方法。其主要思想是采用在同一物理时间弱耦合求解CFD/CSD技术。气动力采用非定常N-S方程的双时间有限体积求解技术,结构响应则采用有限元数值求解技术。CFD和CSD耦合计算的边界信息(气动力和网格)由所设计的界面程序传输。网格信息传输采用守恒体积转换(CVT)方法将CSD计算结构响应位移插值到CFD网格点上。变形已有的CFD网格技术用以确定CFD的变形网格。以位移或载荷的迭代误差为判断耦合计算的收敛标准。最后得到了机翼在Ma=0.8395,α=5.06°时CFD/CSD耦合计算的收敛值。针对计算结果分析了机翼受静气动弹性过程中结构响应和气动特性随时间变化的效应。初步研究结果表明:这种弱耦合方法求解非线性气动弹性问题是可行的。     

非线性壁板颤振分析

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

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

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

Stochastic analysis of two dimensional nonlinear panels with structural damping under random excitation

Stochastic behavior of panels in supersonic flow is investigated to assess the significance of including the damping caused by the strains resulting from axial extension of the panel. The governing equations of motion are based on the Von Karman's large deflection equation and are considered with Kelvin's model of structural damping. The panel under study is two dimensional and simply supported for which the first order piston theory is used to account for the unsteady aerodynamic loading. Transformation of the governing partial differential equation to a set of ordinary differential equations is performed through the Galerkin averaging technique. The statistical response moment equations are generated for two modes using the Fokker–Planck equation with a Gaussian closure scheme. The response moments history and their steady state values are considered. Comparison with previous studies shows that the new nonlinear damping term developed through modeling has a significant effect on the response moments. This is verified through analysis of the effects of external pressure and in-plane force spectral density, air-to-structure mass ratio and structural damping ratio on the mean square value of modal amplitudes.     

气体静压止推轴承静动态性能及振动抑制

为探讨抑制气体轴承自激振动的新手段,研究了圆盘形单个小孔节流的气体静压止推轴承振动问题。通过两平行圆盘间的气体流动方程,推导建立了气体轴承静动态特性分析模型。通过数值分析的方法,研究并总结了供气压力和气膜厚度对气体轴承静态承载力、静态气体质量流量、静动态刚度以及动态阻尼的影响,计算得到了一系列气体轴承静动态特性变化曲线。然后针对系统出现负阻尼情况,通过引入非线性能量阱(NES)来抑制自激振动的问题。研究结果表明,提高供气压力可以有效地提高气体轴承的静态性能;气体轴承的自激振动主要源于气膜挤压的负阻尼特性;同时,当NES的阻尼超过临界阻尼时,可以使系统的振动得到有效的抑制。     

Hypersonic flutter and flutter suppression system of a wind tunnel model

The aeroelastic behavior of a thin flat rudder model was numerically simulated and experimentally investigated in a hypersonic wind tunnel. In particular, a flutter suppression system taking advantage of collision within small gaps was proposed and a novel system for the flutter simulation of the whole nonlinear aeroelastic system including the flutter suppression system was developed. First, the critical flutter dynamic pressure of the rudder without the flutter suppression system was calculated with different methods. Then, the whole nonlinear aeroelastic system, including theflutter suppression system, was simulated to design the gap size. Finally, the flutter suppression system was experimentally validated in a hypersonic wind tunnel operating at Mach number 5. The typical phenomenon of Limit Cycle Oscillation(LCO) was observed, avoiding the structural failure of the model and the consistency between numerical and experimental results was demonstrated.The proposed suppression system can improve the design and reusability of test models of hypersonic flutter experiments.