侧向随动力作用下大展弦比柔性机翼的稳定性

 随动力能够诱发弹性结构发生颤振失稳。以侧向随动力和集中质量分别模拟发动机推力和外挂质量,考虑机翼垂直弯曲-扭转刚度比、集中质量大小、侧向随动力和集中质量的位置以及机翼后掠角和上反角的影响,研究了受侧向随动力作用的大展弦比柔性机翼的气动弹性稳定性。数值模拟所采用的大展弦比柔性机翼非线性气动弹性模型耦合了几何精确完全本征运动梁模型和ONERA动失速非定常气动力模型,该模型考虑了几何非线性、动失速和材料各向异性。模拟结果表明,侧向随动力对机翼颤振可以具有稳定作用,其具体表现依赖于若干变参数的影响,如：减小机翼垂直弯曲-扭转刚度比;发动机吊舱靠近翼根布置;使发动机推力作用点在法向上与机翼弹性轴靠近;单纯的集中质量避免布置在柔性机翼中部,且布置在机翼弹性轴之前或下方,这些设计或布置均有利于提高带发动机吊舱/有效载荷外挂的柔性机翼的气动弹性稳定性。     

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

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

柔性飞机非线性气动弹性与飞行动力学耦合静、动态特性

高空长航时(HALE)飞机结构细长、具有柔性,在常规飞行条件下可发生结构大变形、气动失速以及结构低频振动与刚体运动耦合,这些现象显著影响其静、动态特性.基于几何精确完全本征运动梁模型、ONERA动失速气动力模型和六自由度刚体运动模型,建立了考虑几何非线性、动失速和材料各向异性等因素的大展弦比柔性飞机非线性气动弹性与飞行...     

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.     

基于CFD/CSD的非线性气动弹性分析方法

提出了一种基于计算流体力学/计算结构动力学(CFD/CSD)的非线性气动弹性分析方法,并应用于切尖三角翼的非线性颤振和极限环振荡(LCO)研究。该方法将非线性有限元(FEM)和CFD计算相结合,并辅以高精度的界面插值,能够分析结构和气动非线性共存的气动弹性问题。结构部分以四边形平板壳元为基础,采用更新的拉格朗日(UL)方法分析结构大变形引起的几何非线性问题。气动部分以Navier-Stokes(N-S)方程作为控制方程,采用CFD方法计算跨声速气动力。机翼的非线性颤振计算表明了方法的有效性。最后应用该方法研究了切尖三角翼的LCO现象,其计算精度明显优于已有结果。     

基于非线性试验气动力的飞机静气动弹性响应分析

分别使用线性和非线性静气动弹性分析方法,对某飞机纵向静气动弹性响应特性随飞行动压、攻角、平尾偏度和纵向过载变化的趋势进行了分析,并将分析结果和飞行试验进行了比较。线性方法的气动力计算使用亚音速偶极子格网法。非线性方法由风洞试验提供刚体气动力,并使用线性气动力影响系数矩阵对其进行弹性化处理。通过对两种方法的计算结果进行对比分析可以看出:①使用刚体试验气动力的非线性分析方法能够获得和飞行试验比较一致的结果;②在线性方法所提供的结果中,翼面的弯曲变形及剪力、弯矩和扭矩等部分参数能为初步设计提供大致参考,但气动力系数的弹性增量、翼面的扭转变形、翼面的剪力、弯矩和扭矩的弹性部分可能不能为型号设计提供正确的指导;③不论刚体气动力为线性还是非线性,气动力系数随动压变化的趋势均呈线性。     

低速流场中柔性悬臂板的后颤振响应

建立了一个新的非线性气动弹性模型,对低速流场中柔性悬臂板的后颤振响应特性进行了分析。建模中考虑了结构几何非线性、气动力非线性以及两者之间的强耦合效应。通过实验数据对所建立的气动弹性模型进行了验证。发现在低速流场中柔性悬臂板可能会以周期加倍的方式进入混沌运动。结构几何非线性效应和翼尖涡引起的非定常气动力效应对柔性悬臂板的结构响应有显著影响,而尾涡变形引起的非定常气动力对结构运动的影响较小。还研究了不同耦合算法的差异,给出了小展弦比大柔性结构非线性气动弹性数值仿真时耦合策略的选择依据。     

直升机的气动弹性问题

直升机的气动弹性问题与固定翼飞机不同,不仅要考虑单片桨叶,更要将旋翼视为一个整体,考虑其动态入流、尾迹影响以及旋翼与机身之间的相互耦合等。就单片桨叶而言,在结构动力学上,需要考虑离心力场、几何非线性以及桨叶的非线性挥舞-摆振-扭转耦合;在气动力上,需要考虑动态入流以及桨尖处可能的失速效应,本质上属于非线性气动弹性力学范畴。由于旋翼气动力通常是以周期形式通过旋翼轴传给机身,并引起机身振动,而机身运动又通过改变桨叶根部形态反过来影响旋翼的气动弹性特性,这种旋翼/机身耦合问题,也是近年来直升机气动弹性问题研究中的重要方向和热点之一。此外,随着旋翼流场数值分析方法的日趋成熟,采用动态重叠网格或滑移网格方法来实现桨叶运动,并通过动网格技术来实现桨叶的弹性变形,从而实现弹性旋翼流场的数值模拟,目前正呈现出勃勃生机,成为直升机气动弹性研究的又一重要方向和热点。随着各种新构型直升机的相继出现,如倾转旋翼机、前行桨叶概念旋翼(ABC)直升机和复合式直升机等,也带来了新的气动弹性问题。不断发现问题、解决问题,推动本学科持续发展,永远是气动弹性工作者终身奋斗的目标。     

Flutter and gust response analysis of a wing model including geometric nonlinearities based on a modified structural ROM

In this paper, a framework is established for nonlinear flutter and gust response analyses based on an efficient Reduced Order Model (ROM). The proposed method can be used to solve the aeroelastic response problems of wings containing geometric nonlinearities. A structural modeling approach presented herein describes the stiffness nonlinearities with a modal formulation. Two orthogonal spanwise modes describe the foreshortening effects of the wing. Dynamic linearization of the ROM under nonlinear equilibrium states is applied to a nonlinear flutter analysis, and the fully nonlinear ROM coupled with the non-planar Unsteady Vortex Lattice Method (UVLM) is applied to gust response analysis. Furthermore, extended Precise Integration Method (PIM) ensures accuracy of the dynamic equation solutions. To demonstrate applicability and accuracy of the method presented, a wind tunnel test is conducted and good agreements between theoretical and test results of nonlinear flutter speed and gust response deflection are reached. The method described in this paper is suitable for predicting the nonlinear flutter speed and calculating the gust responses of a large-aspect-ratio wing in time domain. Meanwhile, the results derived highlight the effects of geometric nonlinearities obviously.     

Active Aerothermoelastic Control of Hypersonic Double-wedge Lifting Surface

Designing re-entry space vehicles and high-speed aircraft requires special attention to the nonlinear thermoelastic and aerodynamic instability of their structural components. The thermal effects are important since temperature environment brings dramatic influences on the static and dynamic behaviors of flight structures in supersonic/hypersonic regimes and is likely to cause instability, catastrophic failure and oscillations resulting in structural failure due to fatigue. In order to understand the dynamic behaviors of these ＂hot＂ structures, a double-wedge lifting surface with combining freeplay and cubic structural nonlinearities in both plunging and pitching degrees-of-freedom operating in supersonic/hypersonic flight speed regimes has been analyzed. A third order piston theory aerodynamic is used to estimate the applied nonlinear unsteady aerodynamic loads. Also considered is the loss of torsional stiffness that may be incurred by lifting surfaces subject to axial stresses induced by aerodynamic heating. The aerodynamic heating effects are estimated based on the adiabatic wall temperature due to high speed airstreams. As a recently emerging technology, the active aerothermoelastic control is aimed at providing solutions to a large number of problems involving the aeronautical/aerospace flight vehicle structures. To prevent such damaging phenomena from occurring, an application of linear and nonlinear active control methods on both flutter boundary and post-flutter behavior has been fulfilled. In this paper, modeling issues as well as numerical simulation have been presented and pertinent conclusions outlined. It is evidenced that a serious loss of torsional stiffness may induce the dynamic instability; however active control can be used to expand the flutter boundary and convert unstable limit cycle oscillations （LCO） into the stable LCO and/or to shift the transition between these two states toward higher flight Mach numbers.     

弹性飞行器飞行载荷与动态特性分析

在综合考虑了刚体运动自由度和弹性自由度的弹性飞行器运动方程的基础上 ,采用了线性气动力的影响系数形式 ,提出了一套飞行载荷与动态特性分析的方程和方法 ,用于解决飞行力学和气动弹性力学中的静气动弹性发散、配平与结构变形、气动载荷分布、颤振与飞行动态特性等问题。算例表明 ,该方法有效可行 ,且对于大柔性飞行器 ,刚体运动模态与弹性模态之间的耦合是显著的 ,应予以重视。     

Static aeroelastic analysis of very flexible wings based on non-planar vortex lattice method

A rapid and efficient method for static aeroelastic analysis of a flexible slender wing when considering the structural geometric nonlinearity has been developed in this paper. A non-planar vortex lattice method herein is used to compute the non-planar aerodynamics of flexible wings with large deformation. The finite element method is introduced for structural nonlinear statics analysis. The surface spline method is used for structure/aerodynamics coupling. The static aeroelastic characteristics of the wind tunnel model of a flexible wing are studied by the nonlinear method presented, and the nonlinear method is also evaluated by comparing the results with those obtained from two other methods and the wind tunnel test. The results indicate that the traditional linear method of static aeroelastic analysis is not applicable for cases with large deformation because it produces results that are not realistic. However, the nonlinear methodology, which involves combining the structure finite element method with the non-planar vortex lattice method, could be used to solve the aeroelastic deformation with considerable accuracy, which is in fair agreement with the test results. Moreover, the nonlinear finite element method could consider complex structures. The non-planar vortex lattice method has advantages in both the computational accuracy and efficiency. Consequently, the nonlinear method presented is suitable for the rapid and efficient analysis requirements of engineering practice. It could be used in the preliminary stage and also in the detailed stage of aircraft design.     

Airfoil Aeroelastic Flutter Analysis Based on Modified Leishman-Beddoes Model at Low Mach Number

Based on modified Leishman-Beddoes (L-B) state space model at low Mach number (lower than 0.3), the airfoil aeroelastic system is presented in this paper. The main modifications for L-B model include a new dynamic stall criterion and revisions of normal force and pitching moment coefficient. The bifurcation diagrams, the limit cycle oscillation (LCO) phase plane plots and the time domain response figures are applied to investigating the stall flutter bifurcation behavior of airfoil aeroelastic systems with symmetry or asymmetry. It is shown that the symmetric periodical oscillation happens after subcritical bifurcation caused by dynamic stall, and the asymmetric periodical oscillation, which is caused by the interaction of dynamic stall and static divergence, only happens in the airfoil aeroelastic system with asymmetry. Validations of the modified L-B model and the airfoil aeroelastic system are presented with the experimental airload data of NACA0012 and OA207 and experimental stall flutter data of NACA0012 respectively. Results demonstrate that the airfoil aeroelastic system presented in this paper is effective and accurate, which can be applied to the investigation of airfoil stall flutter at low Mach number.     

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.     

Gust response analysis and wind tunnel test for a high-aspect ratio wing

A theoretical nonlinear aeroelastic response analysis for a flexible high-aspect ratio wing excited by harmonic gust load is presented along with a companion wind tunnel test. A multidisci-plinary coupled numerical calculation is developed to simulate the flexible model wing undergoing gust load in the time domain via discrete nonlinear finite element structural dynamic analysis and nonplanar unsteady vortex lattice aerodynamic computation. A dynamic perturbation analysis about a nonlinear static equilibrium is also used to determine the small perturbation flutter bound-ary. A novel noncontact 3-D camera measurement analysis system is firstly used in the wind tunnel test to obtain the spatial large deformation and responses. The responses of the flexible wing under different static equilibrium states and frequency gust loads are discussed. The fair to good quanti-tative agreements between the theoretical and experimental results demonstrate that the presented analysis method is an acceptable way to predict the geometrically nonlinear gust response for flex-ible wings.     

大展弦比复合材料机翼结构建模与气动弹性分析

建立了具有任意截面形状的大展弦比复合材料翼的结构模型。研究了机翼的两种铺设方式：即周向均匀刚度配置（CUS）和周向反对称刚度配置（CAS）。对于复合箱梁情形,目前的结构建模方法正确性通过ANSYS有限单元软件得到了验证。为了研究具有NACA0012翼型的大展弦比复合材料机翼的气动弹性问题,利用线性ONERA空气动力模型来描述小攻角情形下的非定常空气动力载荷。最后,利用U—g法预示了机翼在各种复合层铺设角度下的颤振速度。     

Numerical studies of a non-linear aeroelastic system with plunging and pitching freeplays in supersonic/hypersonic regimes

The flutter and post flutter of a two-dimensional double-wedge lifting surface with combined freeplay and cubic stiffness nonlinearities in both plunging and pitching degrees-of-freedom operating in supersonic/hypersonic flight speed regimes have been analyzed. In addition to the structural nonlinearities, the third-order piston theory aerodynamics is used to evaluate the unsteady non-linear aerodynamic force and moment. Such model accounts for stiffness and damping contributions produced by the aerodynamic loads. Responses involving limit cycle oscillation and chaotic motion are observed over a large number of parameters that characterizes the aeroelastic system. Results of the present study show that the freeplay in the pitching degree-of-freedom and soft/hard cubic stiffness in the pitching and plunging degrees-of-freedom have significant effects on the LCOs exhibited by the aeroelastic system in the supersonic/hypersonic flight speed regimes. The simulations also show that the aeroelastic system behavior is greatly affected by physical structural parameters, such as the radius of gyration and the frequency ratio, especially in post-flutter regimes, when accounting for all system nonlinearities. It has been shown that at high Mach numbers the non-linear aerodynamic stiffness yields detrimental effects from the aeroelastic point of view, while the damping one do not.     

带有气动及结构非线性的二元机翼颤振分析

研究了翼型在低马赫数条件下的非定常气动特性,从翼型表面气流运动的角度对Leishman-Beddoes(L-B)模型进行了修正,并在此基础上建立了适合低马赫数颤振研究且带有气动及结构非线性的二元机翼气弹系统分析模型.对比低马赫数翼型气动载荷试验结果表明对L-B模型的修正是有效的,且机翼颤振试验结果亦验证了二元机翼气弹分析模型.研究结果表明:二元机翼气弹系统的失速颤振与初始变距角和来流速度密切相关,且耦合的三次非线性变距和浮沉刚度是造成系统呈现准周期运动的主要原因.      

新型倾转旋翼机气动弹性稳定性分析模型

针对新型半铰接式、摆振柔软的倾转旋翼机,建立了旋翼/动力舱/机翼气动弹性稳定性分析模型,并导出了系统的动力学方程。分析模型采用了当量铰、弹性约束的刚性桨叶和弹性机翼假设,采用Pitt-Peters动力入流模型计入了非定常气动力的影响。所建模型为进行新型倾转旋翼机大速度飞行时气弹稳定性分析奠定了基础。     

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.