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

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

NUMERICAL INVESTIGATIONS OF THE AERODYNAMICS AND THE AEROELASTIC STABILITY OF OSCILLATING ANNULAR WINGS

The nacelles of modern aeroengines are constantly increasing in size. Thus, engine air-loads are becoming more powerful and their importance for the aeroelastic stability is becoming more significant. The principal goal of this study is to answer the question of how unsteady airloads vary while shifting to transonic Mach numbers. The investigations are carried out by applying a finite volume Euler method to a harmonically oscillating annular wing. The results show that transonic effects in the case of an annular wing are essentially weaker than in the case of an airfoil. The order of magnitude of the variations is around 10%. Possible consequences for the aeroelastic stability are examined with the example of an elastically mounted annular wing in transonic flow. The shifts of the stability curves also remain within a range of 10%. In addition, an actuator disk method, which is frequently used for the simulation of the fan jet, is expanded in such a way that unsteady flows can be treated. Some unsteady air-loads are strongly dependent on the pressure jump across the fan.     

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

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

可渗透壁对翼型声涡相互作用下流场及声场的影响

由仿生学原理构建的可渗透翼型对湍流气动噪声抑制作用已展现良好的应用前景。对NACA 0012可渗透翼型和实体翼型进行了数值计算,得到了声涡相互作用下气动噪声声场和流场,分析了可渗透壁对翼型流场和声场的影响。研究表明,相对实体翼型,可渗透壁通过减小声源强度降低了主纯音噪声声压级幅值和远场总声压级,消除了高阶离散纯音,但对噪声的指向性没有较大改变。进一步的流场分析表明,可渗透壁对翼型气动性能影响不大的情况下能够降低边界层扰动和翼型后缘大尺度涡旋强度,并推迟分离泡转捩和再附位置。     

壁板颤振的分析模型、数值求解方法和研究进展

研究壁板颤振问题需要计及大挠度变形下结构的几何非线性效应,不仅涉及气动弹性稳定性,而且关心结构的非线性颤振响应.该文回顾了飞行器壁板颤振问题的国内外研究情况,评述了在壁板颤振研究中采用的分析模型、数值求解方法以及在理论分析和试验方面的研究成果,并提出了今后壁板颤振问题的4个研究方向.      

大长细比导弹的气动弹性降阶模型

对于大长细比导弹,需要在设计阶段准确计算气动弹性/气动伺服弹性,但其复杂的气动力给计算带来困难,因此气动力降阶模型是突破大长细比导弹跨音速气动弹性分析与控制瓶颈的关键技术.虽然气动力模型降阶方法已在预测二维机翼结构的气动弹性方面取得重要进展,但几乎未见关于全机模型的气动力降阶模型研究报道.本文基于递归Wiener模型的气动力降阶方法,利用CFD计算的气动力作为模型辨识数据,用鲁棒子空间和Levenberg-Marquardt算法辨识降阶模型参数,建立了大长细比导弹气动力降阶模型.在此基础上与大长细比导弹有限元模型相结合,构造出气动弹性降阶模型,并在数值仿真中测试气动弹性降阶模型在不同马赫数下的适用性.数值仿真结果表明,该气动弹性降阶模型能够精确预测导弹模型在不同飞行条件下的非定常气动力和导弹模型的气动弹性频率响应特性.     

Fluid-structure-acoustic coupling for a flat plate

The present work deals with the aeroacoustic sound radiated by a forward–backward facing step in combination with a flexible wall behind the step. A numerical flow computation with coupled aeroacoustic and vibroacoustic simulation was carried out. The structural deformations of the oscillating plate like structure in the wake of the forward–backward facing step were considered to be small and therefore not affecting the flow field. The presented approach enables a separate consideration for the aeroacoustic as well as the structural borne noise. The influence of the interactions of the acoustic medium with the flexible structure on the vibroacoustic sound radiation is investigated. One-sided and two-sided coupling approaches for the vibroacoustic analysis are introduced. The two-sided vibroacoustic computation allows for considering the damping influence of the ambient fluid on the flexible plate vibration and therefore on the sound radiation. Additional to the simulations, aeroacoustic measurements in an acoustic wind tunnel were performed for validation purposes.     

Aeroelastic effect on modal interaction and dynamic behavior of acoustically excited metallic panels

This paper details the study of the aeroelastic effect on modal interaction and dynamic behavior of acoustically excited square metallic panels with fully clamped edges using finite element method. The first-order shear deformation plate theory and von Karman nonlinear strain–displacement relationships are employed to consider the structural geometric nonlinearity caused by large vibration deflections. Piston aerodynamic theory and Gaussian white noise are used to simulate the aerodynamic load and the acoustic load, respectively. Motion equations are derived by the principle of virtual work in the physical coordinates and then transformed into the truncated modal coordinates with reduced orders. Runge–Kutta method is employed to obtain the system response, and the modal interaction mechanism is quantitatively valued by the modal participation distribution. Results show that in the pre-/near-flutter regions, in addition to the dominant fundamental resonant mode, the first twin companion antisymmetric modes can be largely excited by the aeroelastic coupling mechanism; thus, aeroelastic modal participation distribution and the spectrum response can be altered, while the dynamic behavior still exhibits linear random vibrations. In the post-flutter region, the dominant flutter motion can be enriched by highly ordered odd order super-harmonic motion occurs due to 1:1 internal resonances. Correspondingly, the panel dynamic behavior changes from random vibration to highly ordered motions in the fashion of diffused limit-cycle oscillations (LCOs). However, this LCOs motion can be affected by the intensifying acoustic excitation through changing the aeroelastic modal interaction mechanism. Accompanied with these changes, the panel can experience various stochastic bifurcations.     

航天器噪声试验中结构振动响应预示方法研究

航天器在随运载火箭发射过程中要承受严酷的噪声环境,需通过噪声试验来检验航天器承受噪声环境并能正常工作的能力.航天器噪声试验中结构振动的响应特性是结构强度设计应该考虑的因素之一,更是制定器上组件随机振动试验条件的重要依据,因此有必要在航天器研制初期对噪声载荷作用下的结构振动进行响应预示.文章应用商用有限元分析软件MSC.Patran和MSC.Nastran建立了某型号航天器结构舱板的有限元模型,将噪声载荷声压谱转换为脉动压力功率谱密度,进而采用模态法分析结构在噪声载荷作用下的随机振动响应,并将仿真预示结果与试验结果进行对比研究,在仿真分析中考虑阻尼参数模型和流场附加质量效应等因素的影响;通过研究表明:采用阻尼比随频率提高而减小的经验阻尼参数模型可以较好地反映中高频响应特性、得到较为准确的总均方根响应分析结果,进一步采用虚拟质量法考虑流场附加质量效应可以得到较为准确的功率谱密度响应分析结果.文章提出的仿真分析方法建模简便、计算成本低,适用于在航天器研制初期对航天器噪声试验中的结构振动进行响应预示.      

LAUNCH VEHICLE BUFFETING WITH AEROELASTIC COUPLING EFFECTS

Launch vehicle structural responses can couple with transonic flow state transitions at the nose of payload fairings. This self-sustained coupling yields a nonlinear equation of motion that can be analyzed using the force–response relationship and the periodicity condition. The traditional analysis approach for this phenomenon, however, linearizes the equation of motion by converting the alternating flow forces into an aerodynamic damping term and defines a stability criterion as the response amplitude that yields zero net system damping. This work clarifies the relationship between the present and traditional methods, and compares results and conclusions. The feasibility of modifying a launch vehicle response analysis of buffeting (random pressure fluctuations caused by turbulent flow) to include aeroelastic coupling effects is also explored. The aerodynamic stiffness and damping terms formulated herein are consistent with trends observed in wind-tunnel test data. It is shown, however, that the modified buffet analysis can be inaccurate, particularly when the aeroelastic coupling contribution does not dominate the system response.     

不同路面激励下车内噪声预测与板件声学贡献量分析

利用有限元模拟计算的方法对不同路面激励下车内噪声进行了预测, 并分析了板件对声 学的贡献量. 使用MATLAB软件生成了B级和C级路面谱. 由NASTRAN软件计算 出车身结构在不同路面谱激励下1$\sim$200\,Hz的频率响应结果, 在LMS Virtual. Lab声学软件中, 根据声传递向量技术预测驾驶员和乘客耳旁声压级响应. 还探讨了 声腔结构表面阻抗特性对车内噪声的影响. 并针对C级路面谱激励, 通过车身各板件对 驾驶员右耳声压的面板贡献量分析, 找出了结构优化的重点区域, 为降低车内声压检测点噪 声进行相关的声学处理提供参考依据.     

Frequency-dependent random fatigue of panel-type structures made of ceramic matrix composites

The panel-type structures used in aerospace engineering can be subjected to severe highfrequency acoustic loadings in service. This paper evaluates the frequency-dependent random fatigue of panel-type structures made of ceramic matrix composites(CMCs) under acoustic loadings. Firstly, the high-frequency random responses from the broadband random excitation will result in more stress cycles in a deinite period of time. The probability density distributions of stress amplitudes will be different in different frequency bandwidths, though the peak stress estimations are identical. Secondly, the fatigue properties of CMCs can be highly frequency-dependent. The fatigue evaluation method for the random vibration case is adopted to evaluate the fatigue damage of a representative stiffened panel structure. The frequency effect through S-N curves on random fatigue damage is numerically veriied. Finally, a parameter is demonstrated to characterize the mean vibration frequency of a random process, and hence this parameter can further be considered as a reasonable loading frequency in the fatigue tests of CMCs to obtain more reliable S-N curves.Therefore, the inluence of vibration frequency can be incorporated in the random fatigue model from the two perspectives.     

IDENTIFICATION AND PREDICTION OF UNSTEADY TRANSONIC AERODYNAMIC LOADS BY MULTI-LAYER FUNCTIONALS

Nonlinear unsteady aerodynamic effects present major modelling difficulties in the analysis and control of aeroelastic response. A rigorous mathematical framework, that can account for the complex nonlinearities and time-history effects of the unsteady aerodynamic response, is provided by the use of functional representations. A recent development, based on functional approximation theory, has achieved a new functional form; namely, multi-layer functionals. The development of a multi-layer functional for discrete-time, finite memory, causal systems has been shown to be realizable via finite impulse response neural networks. Identification of an appropriate temporal neural network model of the nonlinear transonic aerodynamic response is facilitated via a supervised training process using multiple input–output sets, with data obtained by an Euler CFD code. The training process is based on a genetic algorithm to optimize the network architecture, combined with a random search algorithm to update weight and bias values. The approach is examined for two different multiple aerodynamic input–output data sets, and in both cases, the prediction properties of the network model establish the multi-layer functional as a suitable representation of unsteady aerodynamic response.     

Development of numerical techniques for near‐field aeroacoustic computations

This work is concerned with the development of a numerical scheme capable of producing accurate simulations of sound propagation in the presence of a mean flow field. The method is based on the concept of variable decomposition, which leads to two separate sets of equations. These equations are the linearised Euler equations and the Reynolds‐averaged Navier–Stokes equations. This paper concentrates on the development of numerical schemes for the linearised Euler equations that leads to a computational aeroacoustics (CAA) code. The resulting CAA code is a non‐diffusive, time‐ and space‐staggered finite volume code for the acoustic perturbation, and it is validated against analytic results for pure 1D sound propagation and 2D benchmark problems involving sound scattering from a cylindrical obstacle. Predictions are also given for the case of prescribed source sound propagation in a laminar boundary layer as an illustration of the effects of mean convection. Copyright © 1999 John Wiley & Sons, Ltd.     

Nonlinear dynamics of an aeroelastic airfoil with free-play in transonic flow

Nonlinear dynamic behaviors of an aeroelastic airfoil with free-play in transonic air flow are studied. The aeroelastic response is obtained by using time-marching approach with computational fluid dynamics (CFD) and reduced order model (ROM) techniques. Several standardized tests of transonic flutter are presented to validate numerical approaches. It is found that in time-marching approach with CFD technique, the time-step size has a significant effect on the calculated aeroelastic response, especially for cases considering both structural and aerodynamic nonlinearities. The nonlinear dynamic behavior for the present model in transonic air flow is greatly different from that in subsonic regime where only simple harmonic oscillations are observed. Major features of the responses in transonic air flow at different flow speeds can be summarized as follows. The aeroelastic responses with the amplitude near the free-play are dominated by single degree of freedom flutter mechanism, and snap-though phenomenon can be observed when the air speed is low. The bifurcation diagram can be captured by using ROM technique, and it is observed that the route to chaos for the present model is via period-doubling, which is essentially caused by the free-play nonlinearity. When the flow speed approaches the linear flutter speed, the aeroelastic system vibrates with large amplitude, which is dominated by the aerodynamic nonlinearity. Effects of boundary layer and airfoil profile on the nonlinear responses of the aeroelastic system are also discussed.     

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

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

High‐order numerical simulations of flow‐induced noise

In this paper, the flow/acoustics splitting method for predicting flow‐generated noise is further developed by introducing high‐order finite difference schemes. The splitting method consists of dividing the acoustic problem into a viscous incompressible flow part and an inviscid acoustic part. The incompressible flow equations are solved by a second‐order finite volume code EllipSys2D/3D. The acoustic field is obtained by solving a set of acoustic perturbation equations forced by flow quantities. The incompressible pressure and velocity form the input to the acoustic equations. The present work is an extension of our acoustics solver, with the introduction of high‐order schemes for spatial discretization and a Runge–Kutta scheme for time integration. To achieve low dissipation and dispersion errors, either Dispersion‐Relation‐Preserving (DRP) schemes or optimized compact finite difference schemes are used for the spatial discretizations. Applications and validations of the new acoustics solver are presented for benchmark aeroacoustic problems and for flow over an NACA 0012 airfoil. Copyright © 2010 John Wiley & Sons, Ltd.     

Cavity Tones by Computational Aeroacoustics

The flow past open cavities is a problem that is encountered in many engineering applications and can result in intense acoustic tones. The flow physics and acoustics of cavity configurations are complex and computational simulation techniques provide an opportunity to gain further understanding and provide a tool to predict not only cavity tone frequencies but their amplitude. In this paper, we describe the available techniques for performing computational aeroacoustic simulations of cavity flows, and review recent applications for the prediction and control of cavity tones in subsonic, transonic and supersonic regimes.     

Direct numerical simulation of flow over a forward-facing step – Flow structure and aeroacoustic source regions

A direct numerical simulation of the flow over a forward-facing step at a Reynolds number of 8000 based on the step height is presented. Calculations were performed using second-order finite volume discretisation in space on co-located meshes. A hybrid calculation approach based on Lighthill’s acoustic analogy is explained. Results of the simulation are intended to be used as a database for the validation of different discretisation schemes for the flow computation and simulation approaches for the calculation of sound radiation using a hybrid approach. Turbulent statistics are presented along with aeroacoustic source regions. Strong and weak forms of the aeroacoustic source term are presented and compared. For visualization purposes, the strong form is more suitable, whereas for the calculation of sound radiation both forms can be used. From the visualization of the aeroacoustic sources, it can be seen that they mainly concentrate on the region of the leading edge of the step and the shear layer close to the step.     

点、线接触真实粗糙表面的弹流润滑研究

本文给出了点、线接触的真实粗糙表面的微弹流数值解。在给定随机粗糙表面样本后,求解大小不同的载荷和粗糙的弹流问题。从计算结果可以看出,由于Reynolds方程中速度项的作用,在表面对应粗糙的位置处引起了压力变化,从而因其产生的弹性变形使粗糙变得平滑。对此光滑表面解可以看出,因粗糙引起的压力和膜厚的变化在光滑解附近波动。载荷较大时,压力分节接近固体接触情况。