基于ROM技术的阵风响应分析方法

阵风响应分析是大型飞机设计过程中必不可少的环节. 现有的阵风响应分析主要采用基于线化升力面理论的气动力模型,不能考虑到各种非线性效应,不适合于跨音速气动弹性的分析. 基于CFD技术,采用系统辨识方法,在状态空间内建立了降阶的非定常气动力模型(reduced order model, ROM). 耦合结构运动方程、非定常气动力模型(结构运动)、外激阵风的气动力模型,建立了基于CFD技术的阵风响应分析模型.算例研究了某一典型机翼在方波激励下的阵风响应问题,对比了各阶模态位移的响应以及翼根弯矩的响应. 基于ROM技术的计算结果与CFD/CSD直接耦合仿真结果吻合,证明了该方法的正确性和精度.      

基于特征正交分解的非定常气动力建模技术

采用特征正交分解(proper orthogonal decomposition, POD)方法, 建立了基于状态空间的非定常气动力降阶模型, 并耦合结构方程, 建立了降阶的气动弹性系统, 开展了颤振分析的初步研究, 计算效率提高了2~3个数量级. 具体过程是:首先获取全阶系统的频域快照构成关联矩阵, 通过对关联矩阵进行奇异值分解提取流场模态(或流场基), 对低能量模态截断形成降阶子空间, 并将其映射到全阶系统, 从而形成基于状态空间的降阶非定常气动力模型. 对气动弹性标模AGARD445.6进行算例验证, 证明了降阶方法正确, 可以提供高效、高精度的气动弹性分析.      

基于POD和代理模型的静气动弹性分析方法

静气动弹性问题考虑弹性结构与定常气动力间的相互耦合作用,对飞行器的性能和安全具有显著的影响.在现代飞行器设计阶段,计算流体力学(CFD)/计算结构力学(CSD)直接耦合方法是精确考察静气动弹性影响的重要手段.然而,基于CFD技术的气动力仿真手段在耦合过程中计算量大且耗时长,难以满足设计阶段的需求.因此,为了兼顾计算精度与效率,文章采用本征正交分解(POD)和Kriging代理模型相结合的模型降阶方法,替代CFD求解过程并耦合有限元分析(FEA)方法,建立了高效、准确的静气动弹性分析框架.相较于传统的以模态法为主的静气动弹性分析方法,该方法能够解决更为复杂的静气动弹性问题以及提供静气动弹性变形过程中的气动分布载荷.针对典型三维跨声速HIRENASD机翼模型开展的马赫数、迎角变化的算例验证表明:由建立的静气动弹性分析方法与CFD/CSD直接耦合方法计算得到机翼翼梢处的静变形量间的相对误差在5%以内;同时该方法预测静平衡位置处的气动分布载荷的误差在5%以内,静气动弹性分析的计算效率至少提升了6倍.     

基于SCI/ERA方法的高效气动力降阶模型

由正交Walsh函数构造Walsh-单信号-复合-输入,对其作用下的计算流体力学响应采用单信号-复合-输入/特征系统实现算法SCI/ERA(Single-Composite-Input/Eigensystem Realization Algorithm)辨识得到离散时间非定常气动力状态空间降阶模型。通过对Isogai机翼剖面气动弹性算例的计算证明该方法具有和非定常计算流体力学方法相当的精度同时模型维数降低2个数量级;在模型构造时间上,SCI/ERA方法比脉冲/ERA方法计算效率提高24%,同时内存占用减小34%;由理论分析可知当耦合结构模态数目增加时,SCI/ERA方法所需的计算开销增幅远小于脉冲/ERA方法;采用频域平衡特征正交分解BPOD(Balanced Proper Orthogonal Decomposition)方法可以准确地从降阶模型中提取出一个低频二次降阶模型,同时保持与原模型相当的精度。二次降阶后模型维数进一步减小88%。     

基于流形切空间插值的折叠翼参数化气动弹性建模

变体飞行器的气动弹性力学建模是当前先进飞行器设计的研究热点和难点. 然而传统的气动弹性动力学建模方法对于具有结构参变特性的变体飞行器气动弹性力学研究存在建模效率低、计算复杂等问题. 本研究提出了一种基于流形切空间插值的可折叠式变体机翼参数化气动弹性建模方法. 首先, 该方法建立若干个典型折叠角下的折叠翼结构有限元模型, 通过流形切空间插值方法建立折叠翼参数化结构动力学模型. 其次, 采用偶极子网格法得到参数化非定常气动力模型, 进而建立气动和结构相互耦合的折叠翼参数化气动弹性模型. 为了验证该参数化建模方法在折叠翼气动弹性分析中的准确性, 本文以一小展弦比折叠翼为研究对象, 从折叠翼自由振动时的参变模态特性、颤振边界预测两方面进行了算例验证, 并与直接计算方法进行了对比, 进一步验证了参数化气动弹性模型的有效性. 研究结果表明, 该参数化气动弹性模型对上述两类问题的计算精度与直接计算方法一致, 并且有着计算效率更高的优势.      

基于CFD的气动力建模及其在气动弹性中的应用

CFD技术为带有气动力非线性的气动弹性分析提供了一种研究途径,但是基于CFD的气动弹性直接数值模拟方法的计算量很大,不便开展定性分析和参数设计.基于CFD的非定常气动力模型的降阶技术为缓解计算效率与计算精度之间的矛盾以及系统的复杂性与易分析、易设计性之间的矛盾提供了行之有效的技术途径.综述了近年来发展的两类基于CFD技术的非定常气动力降阶技术和一种非线性气动力的谐波平衡方法,以及这些方法在非线性气动弹性研究中的运用.对比了各种方法的优越性并作了进一步的展望.      

Parametric aeroelastic modeling based on component modal synthesis and stability analysis for horizontally folding wing with hinge joints

To investigate the aeroelastic stability of a folding wing effectively, a parametric aeroelastic analysis approach is proposed. First, the fixed interface component modal synthesis is used to derive the structural dynamic equation for a folding wing, in which the elastic connection is considered. The unsteady aerodynamic model is established by the doublet lattice method (DLM), and the aeroelastic model is achieved from integration of the DLM with the component modal analysis. The generalized aerodynamic influence coefficient matrix is established by modes kept and constraint modes of each component. The aeroelastic stability of a folding wing is investigated based on the Gram matrix in control theory. The effectiveness of the proposed method is verified via comparison with traditional flutter eigenvalue analysis for both extended and folded configurations. The proposed method identifies coupled modes and improves computational efficiency when compared to classical aeroelastic stability analysis methods, such as the p–k method.     

跨音速极限环型颤振的高效数值分析方法

事先建立一个低阶的非线性、非定常气动力模型是开展非线性流场中气动弹性问题研究的一个捷径. 基于CFD方法, 通过计算结构在流场中自激振动的响应来获得系统的训练数据. 采用带输出反馈的循环RBF神经网络, 建立时域非线性气动力降阶模型.耦合结构运动方程和非线性气动力降阶模型, 采用杂交的线性多步方法计算结构在不同速度(动压)下的响应历程, 从而获得模型极限环随速度(动压)变化的特性. 两个典型的跨音速极限环型颤振算例表明, 基于气动力降阶模型方法的计算结果与直接CFD仿真结果吻合很好, 与后者相比其将计算效率提高了1~2个数量级.      

Parametric active aeroelastic control of a morphing wing using the receptance method

This study proposes the parametric active aeroelastic control of a folding wing, which is a promising concept of morphing wings, by integrating the parameterized aeroservoelastic model and the receptance-based control strategy. It starts with establishing the parameterized aeroservoelastic model of the folding wing with respect to its folding angle and air speed. The transfer functions between the embedded sensors and the actuators of the folding wing are then efficiently obtained using the parameterized aeroservoelastic model. Finally, control gains for varying folding angle and air speed of the wing are synthesized using the receptance method. The results of numerical analyses show that smooth evolutions of the control gains can be obtained over a wide range of folding angles and air speeds, although the dynamic behavior of the folding wing is sensitive to these parameters. In addition, the active aeroelastic control can effectively suppress aeroelastic vibrations and expand flutter boundaries of the folding wing. The benefit of the parametric active aeroelastic control for the folding wing, i.e., smoothly switching among the control laws without exciting undesirable vibrations is highlighted.     

尾流激励的叶片气动力降阶模型

气动力降阶模型是研究叶片气动弹性振动快速高效的新方法。现有气动力降阶模型的研究主要集中在叶片颤振方面,没有涉及更为常见的上游尾流激励的叶片振动问题。本文提出基于Volterra级数的尾流激励叶片气动力降阶模型,为尾流激励下叶片振动和动静叶干涉振动研究提供了新的思路。采用行波法简化尾流的参数个数,用阶跃信号法识别降阶模型的核函数。二维叶片的算例结果表明,本文方法可以较准确地描述尾流激励引起的叶片气动力振荡,而且计算效率极高。     

Unsteady aerodynamic reduced-order modeling of an aeroelastic wing using arbitrary mode shapes

Computational fluid dynamics (CFD) based unsteady aerodynamic reduced-order model (ROM) can offer significant improvements to the efficiency of transonic aeroelastic analysis. To construct a ROM based on mode shapes, one run of CFD solver is needed to compute aerodynamic responses corresponding to mode excitations. When mode shapes change with structure, another run of the CFD solver is required to construct the new ROM. The typically large computational cost associated with repeated runs of the CFD solver impedes the application of existing unsteady aerodynamic reduced-order modeling methods to transonic aeroelastic design optimization and aeroelastic uncertainty analysis. This paper demonstrates a method that can replace the CFD solver used in the process of existing unsteady aerodynamic reduced-order modeling. It can produce aerodynamic responses corresponding to mode excitations for arbitrary mode shapes within a few seconds. Computational cost can be reduced by two orders of magnitude using the mode excitations and the corresponding aerodynamic responses computed by the method to construct the ROMs used for flutter analyses in aeroelastic design optimization or aeroelastic uncertainty analysis in transonic regime compared with the existing unsteady aerodynamic reduced-order modeling methods. Results show that the method can accurately produce the aerodynamic responses corresponding to the mode excitations and predict the flutter characteristics of AGARD 445.6 wings root-attached in three different ways.     

An efficient implementation of aeroelastic tailoring based on efficient computational fluid dynamics-based reduced order model

Current and future trends in the aerospace industry leverage on the potential benefits provided by lightweight materials that can be tailored to realize desired mechanical characteristics when loaded. For aircraft design, the deployment of aeroelastic tailoring is hindered by the need to re-compute, for any possible modification of the structure, the dependence of the aerodynamic field on the underlying structural properties. To make progress in this direction, the work presents a rapid computational fluid dynamics based aeroelastic tool which is built around a reduced order model for the aerodynamics that is updated for any modification of the structure by using the structural dynamics reanalysis method. The aeroelastic tailoring tool is demonstrated in transonic flow for the AGARD 445.6 wing, suitably modified with composite materials. It was found that the proposed method provides accurate engineering predictions for the aeroelastic response and stability when the structure is modified from the baseline model.     

基于气动力降阶模型的跨音速气动弹性稳定性分析

基于离散型输入输出差分模型,运用非定常CFD方法训练信号,然后运用最小二乘方法进行参数辨识,得到降阶的非定常气动力模型,再将该离散差分模型转换为连续时间域内的状态方程。耦合气动状态方程和结构状态方程,得到耦合系统的气动弹性状态方程。求解不同动压下状态矩阵的特征值,根据根轨迹图分析系统的稳定性特性。分析结果与直接耦合CFD/CSD方法结果相吻合,可以计算跨音速非线性气动弹性问题。其计算效率比直接耦合CFD/CSD方法提高1~2个数量级。针对Isogai wing在跨音速出现的S型颤振边界进行了较为细致的分析,阐述了该现象是由于系统诱发颤振的分支随着速度(来流动压)的提高而发生转移所导致的。     

Suppression of the wing-rock phenomenon using nonlinear controllers

This paper deals with the design of nonlinear controllers for the wing-rock phenomenon of a delta wing aircraft. A fifth order dynamic model is used to describe this phenomenon. A state transformation is introduced such that the transformed dynamic model is in a form which is suitable for a variety of control designs. A feedback linearization control scheme and a sliding-mode control (SMC) scheme are then proposed to suppress wing rock oscillations. It is shown that the two controllers successfully suppress the undesired oscillations and guarantee the asymptotic convergence of all system trajectories to their desired values. The effectiveness of the proposed controllers is verified through simulation studies.

     

Evaluation of linear,inviscid, viscous,and reduced-order modelling aeroelastic solutions of the AGARD 445.6 wing using root locus analysis

Reduced-order modelling (ROM) methods are applied to the Computational Fluid Dynamics (CFD)-based aeroelastic analysis of the AGARD 445.6 wing in order to gain insight regarding well-known discrepancies between the aeroelastic analyses and the experimental results. The results presented include aeroelastic solutions using the inviscid Computational Aeroelasticity Programme–Transonic Small Disturbance (CAP-TSD) code and the FUN3D code (Euler and Navier–Stokes). Full CFD aeroelastic solutions and ROM aeroelastic solutions, computed at several Mach numbers, are presented in the form of root locus plots in order to better reveal the aeroelastic root migrations with increasing dynamic pressure. Important conclusions are drawn from these results including the ability of the linear CAP-TSD code to accurately predict the entire experimental flutter boundary (repeat of analyses performed in the 1980s), that the Euler solutions at supersonic conditions indicate that the third mode is always unstable, and that the FUN3D Navier–Stokes solutions stabilize the unstable third mode seen in the Euler solutions.     

变后掠翼的参变气动弹性建模与分析

参变气动弹性建模是可变后掠翼气动弹性研究中的难点之一. 当地建模技术是一种构造线性参变 (LPV) 模型的实用方法, 但一直缺乏有效途径来解决当地气动弹性模型的不一致问题. 气动弹性模型的不一致性体现在当地结构动力学模型和非定常气动力模型随参数变化的不连续性. 本文提出了一种自下而上的方法, 对变后掠翼不一致的当地气动弹性模型进行了一致性处理. 首先, 采用匈牙利算法跟踪结构模态并按模态分支进行排序, 使得匹配后的模态能够保证结构动力学模型的一致性; 其次, 对有理函数拟合表达式中的系数矩阵进行缩放处理, 解决了空气动力系数矩阵的不一致问题. 采取上述两项措施后, 最终生成了一致的当地状态空间气动弹性模型. 这样, 对一致的状态空间模型进行插值, 就可快速生成任意后掠角下的气动弹性模型, 使得系统的稳定性分析和慢参变响应计算得以高效进行. 数值仿真结果验证了一致性处理的必要性:未经一致性处理的原始模型经插值后得到的模型会出现严重的模型误差. 本文为变后掠翼的参变气动弹性系统建模提供了一种实用、准确和高效的建模方法.      

Reduced order nonlinear aeroelasticity of swept composite wings using compressible indicial unsteady aerodynamics

Nonlinear dynamic aeroelasticity of composite wings in compressible flows is investigated. To provide a reasonable model for the problem, the composite wing is modeled as a thin walled beam (TWB) with circumferentially asymmetric stiffness layup configuration. The structural model considers nonlinear strain displacement relations and a number of non-classical effects, such as transverse shear and warping inhibition. Geometrically nonlinear terms of up to third order are retained in the formulation. Unsteady aerodynamic loads are calculated according to a compressible model, described by indicial function approximations in the time domain. The aeroelastic system of equations is augmented by the differential equations governing the aerodynamics lag states to derive the final explicit form of the coupled fluid-structure equations of motion. The final nonlinear governing aeroelastic system of equations is solved using the eigenvectors of the linear structural equations of motion to approximate the spatial variation of the corresponding degrees of freedom in the Ritz solution method. Direct time integrations of the nonlinear equations of motion representing the full aeroelastic system are conducted using the well-known Runge–Kutta method. A comprehensive insight is provided over the effect of parameters such as the lamination fiber angle and the sweep angle on the stability margins and the limit cycle oscillation behavior of the system. Integration of the interpolation method employed for the evaluation of compressible indicial functions at any Mach number in the subsonic compressible range to the derivation process of the third order nonlinear aeroelastic system of equations based on TWB theory is done for the first time. Results show that flutter speeds obtained by the incompressible unsteady aerodynamics are not conservative and as the backward sweep angle of the wing is increased, post-flutter aeroelastic response of the wing becomes more well-behaved.     

AEROELASTIC OPTIMIZATION DESIGN OF COMPOSITE WING WITH CONSIDERATION OF STRUCTURAL INSTABILITY

从工程数学求解和有限元分析角度对复合材料结构的稳定性分析方法进行研究,基于这两个方面分别建立了同时考虑壁板稳定性约束和气动弹性约 束的气动弹性优化技术,并以大展弦比复合材料机翼为对象,进行气动弹性综合优化设计。研究表明,机翼气动弹性优化中若不考虑稳定性约束条件,虽然可以获得较小结构重量,但往往不满足稳定性要求;相比从有限元角度考虑结构失稳特征的气动弹性综合优化设计方法,通过工程数学方法对机翼结构进行分区失稳分析优化可以更加精准地控制变量,在满足各项性能指标,特别是稳定性约束的同时,进一步减轻了结构重量,提高了结构失稳因子。     

An sequential optimization and aeroelastic constraint transformation method for strength-aeroelastic comprehensive design

With the development of modern aircraft technology, aeroelasticity plays a more and more crucial role in aircraft structural design. However, low efficiency of present aeroelastic analysis and optimization methods makes it difficult to apply in engineering practice. This paper presents a sequential optimization and aeroelastic constraint transformation method (SOACTM) for comprehensive design of airplane wings with strength and aeroelastic constraints. Optimization with structural strength constraint and aeroelastic constraint is transformed into a serial of cycles of decoupled structural strength sub-optimizations and aeroelastic sub-optimizations based on sequential optimization strategy. In structural strength sub-optimization, structural strength constraint is translated along its normal direction to make optimal design point satisfying aeroelastic constraint. And the goal of aeroelastic sub-optimization is to find the translational distance of structural strength constraint. Aeroelastic constraint is transformed to equivalent structural strength constraint via above approach. In this way, number of aeroelastic analyses in SOACTM is less than that in traditional optimization method and total computational time decreases. SOACTM is verified based on two examples. Traditional optimization method is applied for the sake of validation. The results demonstrated the accuracy and efficiency of SOACTM for wing comprehensive optimization considering both structural strength and aeroelastic constraints.     

Efficient unsteady aerodynamic loads prediction based on nonlinear system identification and proper orthogonal decomposition

In the present work, an efficient surrogate-based framework is developed for the prediction of motion-induced surface pressure fluctuations and integral force and moment coefficients. The model construction is realized by performing forced-motion computational fluid dynamics (CFD) simulations, while the result is processed via the proper orthogonal decomposition (POD) to obtain the predominant flow modes. Subsequently, a nonlinear system identification is carried out with respect to the applied excitation and the resulting POD coefficients. For the input/output model identification task, a recurrent local linear neuro-fuzzy approach is employed in order to capture the linear and nonlinear characteristics of the dynamic system. Once the reduced-order model (ROM) is trained, it can substitute the flow solver within unsteady aerodynamic or aeroelastic simulation frameworks for a given configuration at fixed freestream conditions. For demonstration purposes, the ROM approach is applied to the LANN wing in high subsonic and transonic flow. Due to the characteristic lambda-shock system, the unsteady aerodynamic surface pressure distribution is dominated by nonlinear effects. Numerical investigations show a good correlation between the results obtained by the ROM methodology in comparison to the full-order CFD solution. In addition, the surrogate approach yields a significant speed-up regarding unsteady aerodynamic calculations, which is beneficial for multidisciplinary computations.