基于边界元方法的气动弹性降阶模型及仿真

传统气动弹性的时域计算耗费了大量时间,为了提高计算效率,本文发展了基于边界元方法的降阶模型技术。首先基于边界元方法建立非定常流场的求解模型,结合特征值分析技术建立了非定常气动力的低阶模型;然后,利用边界元方法建立了气动网格和结构网格之间的信息转换矩阵;最后将非定常气动力降阶模型和结构动力学方程联合,建立了气动弹性系统的低阶状态空间模型。将所发展的降阶模型方法应用于NACA0012翼型的非定常气动力求解中,结果表明降阶模型可以在保证原系统计算精度的同时提高了计算效率;将降阶模型技术应用到三维机翼的气动弹性响应计算中,在系统阶数仅为12阶的情况下可以得到与原系统一致的极限环响应,说明降阶模型技术在求解气动弹性问题中的巨大优势。     

Method of calculating unsteady aerodynamic characteristics of a lifting surface at high subsonic speed

A method is presented for calculating the unsteady aerodynamic characteristics of harmonically oscillating thin wings traveling at high subsonic speed. The medium is assumed ideal. The aerodynamic coefficients are expressed in terms of the rotational derivatives, which are determined for a Strouhal number of zero. The calculation of the rotational derivatives of the aerodynamic coefficients in a compressible medium reduces to the conversion of the corresponding characteristics of a transformed wing, determined in an incompressible medium for altered boundary conditions. To calculate the aerodynamic characteristics of the transformed wing in the incompressible medium we use a technique based on replacement of the lifting surface by a system of discrete unsteady vortices. The problem is solved in general form, and together with the new relations for the rotational derivatives with dots we derive the known formulas for the rotational derivatives without dots.     

应用当地流活塞理论的大攻角升力面颤振气动力表达式

本文应用当地流活塞理论,给出了弹性的夺动翼面的非定常压分分布以及用模态坐标的广义气动力系数表达式。同时提供了配套的当地流参数计算公式和采用样条函数计算广义气动力系数数值积分的表达式。     

基于三路舵机驱动的扑翼飞行模型气动特性研究

扑动而形成非定常气动现象是扑翼飞行过程中产生高升力的主要原因。本文以Ellington实验的鹰蛾翅膀为原形,设计扑翼实验及数值计算模型。通过压差传感器对翅膀模型上翼面固定位置进行测压,分析前缘涡的产生及脱落情况(考虑动压效应)。测量上下翼面固定位置处的压差,揭示扑翼飞行中产生高升力的主要原因。利用烟风洞观察扑翼模型周围流场结构及特殊涡产生变化情况。另外,根据Ellington提供的升力关系式估算了扑翼模型在一个周期内的平均升力。最后,基于三维欧拉方程对扑翼飞行气动特性进行数值模拟,计算结果与实验吻合良好。     

Calculation of nonlinear aerodynamic characteristics of a wing using a 3‐D panel method

The nonlinear aerodynamic characteristic of a wing is investigated using the frequency‐domain panel method. To calculate the nonlinear aerodynamic characteristics of a three‐dimensional wing, the iterative decambering approach is introduced into the frequency‐domain panel method. The decambering approach uses the known nonlinear aerodynamic characteristic of airfoil and calculates two‐variable decambering function to take into consideration the boundary‐layer separation effects for the each section of the wing. The multidimensional Newton iteration is used to account for the coupling between the different sections of wing. The nonlinear aerodynamic analyses for a rectangular wing, a tapered wing, and a wing with the control surface are performed. Present results are given with experiments and other numerical results. Computed results are in good agreement with other data. This method can be used for any wing having different nonlinear aerodynamic characteristics of airfoil. The present method will contribute to the analysis of aircraft in the conceptual design because the present method can predict the nonlinear aerodynamic characteristics of a wing with a few computing resources and significant time. Copyright © 2007 John Wiley & Sons, Ltd.     

Mechanism of unsteady aerodynamic heating with sudden change in surface temperature

The characteristics and mechanism of unsteady aerodynamic heating of a transient hypersonic boundary layer caused by a sudden change in surface temperature are studied. The complete time history of wall heat flux is presented with both analytical and numerical approaches. With the analytical method, the unsteady compressible boundary layer equation is solved. In the neighborhood of the initial and final steady states, the transient responses can be expressed with a steady-state solution plus a perturbation series. By combining these two solutions, a complete solution in the entire time domain is achieved. In the region in which the analytical approach is applicable, numerical results are in good agreement with the analytical results, showing reliability of the methods. The result shows two distinct features of the unsteady response. In a short period just after a sudden increase in the wall temperature, the direction of the wall heat flux is reverted, and a new inflexion near the wall occurs in the profile of the thermal boundary layer. This is a typical unsteady characteristic. However, these unsteady responses only exist in a very short period in hypersonic flows, meaning that, in a long-term aerodynamic heating process considering only unsteady surface temperature, the unsteady characteristics of the flow can be ignored, and the traditional quasi-steady aerodynamic heating prediction methods are still valid.     

Calculation of unsteady supersonic flow past a two-dimensional cascade of plates ween vorticity inhomogeneities of the flow act on it

In the framework of the linear theory of small perturbations the problem of unsteady subsonic flow past a two-dimensional cascade of plates has been considered in a number of papers. Thus, the unsteady aerodynamic characteristics of a cascade of vibrating plates were calculated in [1] by the method of integral equations, while the same method was used in [2, 3] to calculate the sound fields that are excited when sound waves Coming from outside or vorticity inhomogeneities of the oncoming flow act on the cascade. The problem of a two-dimensional cascade of vibrating plates in a supersonic flow was solved in [4, 5]. In [4] the solution was constructed on the basis of the well-known solution of the problem of vibrations of a single plate, while in [5] a variant of the method of integral equations was used which differed slightly from the usual formulation of this method [1–3]. The approach proposed in [5] is used below to calculate the unsteady flow past a two-dimensional cascade of plates in the case when vorticity inhomogeneities of a supersonic oncoming flow act on it. Equations are obtained for the strength of the unsteady pressure jumps arising in such a flow and the vortex wakes shed from the trailing edges of the plates. Examples of the calculations illustrating the accuracy of the method and its possibilities are given.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp, 152–160, May–June, 1986.     

模型昆虫翼作非定常i运动时的气动力特性

基于Navier-Stokes方程的数值解,研究了一模型昆虫翼在小雷诺数(Re=100)下作非定常运动时的气动力特性.这些运动包括翼启动后的常速转动,快速加、减速转动,常速转动中快速上仰(模拟昆虫翼的上挥或下拍、翻转等运动).有如下结果在小雷诺数下,模型昆虫翼以大攻角(α=35°)作常速转动运动时,由于失速涡不脱落,可产生较大的升力系数.其机理是翼转动时,翼尖附近(该处线速度大)上翼面压强比翼根附近(该处线速度小)的小得多,因而存在展向压强梯度,同时存在着沿展向的离心力,此展向压强梯度和离心力导致的展向流动在失速涡的轴向方向,其可避免失速涡脱落.模型昆虫翼在快速加、减速转动和快速上仰运动中,虽然雷诺数小,但由于在短时间内产生了大涡量,也可产生十分大的气动力,例如在快速上仰运动中,升力系数可大于10.     

Predicting the wind noise from the pantograph cover of a train

Finite element and boundary element calculations are combined to predict the flow noise radiated from a 1/10th-scale model of an aerodynamic cover used around the pantograph on a train at 250 km h⁻¹. The solutions of the unsteady air flow over the cover and the resulting sound propagation are divided into two parts in order to keep the problem tractable. First the unsteady fluid flow is solved using large-eddy simulation (LES). The pressure histories on the cover are then used to predict the radiated sound, using a boundary element method to solve the Helmholtz equation. The result thus leans heavily on assumptions about the coupling of the two solutions, the propagation of sound in a disturbed medium and the efficacy of LES. The predicted sound pressure levels are compared with experimental measurements made in an anechoic wind tunnel. © 1997 John Wiley & Sons, Ltd.     

局部粗糙机翼表面气动特性影响研究

采用柱型粗糙元,以尖劈和机翼外形为基础,利用CFD数值模拟方法,研究机翼表面局部粗糙区域对其周围气动特性的影响。研究相同粗糙元高度下,粗糙元位于尖劈表面不同位置时局部边界层和气动特性的变化情况;基于流动分区理论,采用空气动力学理论分析与数值模拟结合的方法,分析F16机翼可接受的局部粗糙元高度;根据分析结果,在改进平板外形基础上,验证不同粗糙元高度对改进目标块区域气动特性的影响程度,并给出流经局部粗糙区域的流体发展状况。为了验证数值结果的准确性,采用S-A与SST湍流模型进行对比求解。本文工作对复杂大气环境引起的飞行器局部粗糙表面区域气动特性变化的研究具有指导意义。     

Analytical formulation of 2-D aeroelastic model in weak ground effect

This paper deals with the aeroelastic modeling and analysis of a 2-D oscillating airfoil in ground effect, elastically constrained by linear and torsional springs and immersed in an incompressible potential flow (typical section) at a finite distance from the ground. This work aims to extend Theodorsen theory, valid in an unbounded flow domain, to the case of weak ground effect, i.e., for clearances above half the airfoil chord. The key point is the determination of the aerodynamic loads, first in the frequency domain and then in the time domain, accounting for their dependence on the ground distance. The method of images is exploited in order to comply with the impermeability condition on the ground. The new integral equation in the unknown vortex distribution along the chord and the wake is solved using asymptotic expansions in the perturbation parameter defined as the inverse of the non-dimensional ground clearance of the airfoil. The mathematical model describing the aeroelastic system is transformed from the frequency domain into the time domain and then in a pure differential form using a finite-state aerodynamic approximation (augmented states). The typical section, which the developed theory is applied to, is obtained as a reduced model of a wing box finite element representation, thus allowing comparison with the corresponding aeroelastic analysis carried out by a commercial solver based on a 3-D lifting surface aerodynamic model. Stability (flutter margins) and response of the airfoil both in frequency and time domains are then investigated. In particular, within the developed theory, the solution of the Wagner problem can be directly achieved confirming an asymptotic trend of the aerodynamic coefficients toward the steady-state conditions different from that relative to the unbounded domain case. The dependence of flutter speed and the frequency response functions on ground clearance is highlighted, showing the usefulness of this approach in efficiently and robustly accounting for the presence of the ground when unsteady analysis of elastic lifting surfaces in weak ground effect is required.     

An analytic approach to theoretical modeling of highly unsteady viscous flow excited by wing flapping in small insects

Numerous studies on the aerodynamics of insect wing flapping were carried out on different approaches of flight investigations, model experiments, and numerical simulations, but the theoretical modeling remains to be explored. In the present paper, an analytic approach is presented to model the flow interactions of wing flapping in air for small insects with the surrounding flow fields being highly unsteady and highly viscous. The model of wing flapping is a 2-D flat plate, which makes plunging and pitching oscillations as well as quick rotations reversing its positions of leading and trailing edges, respectively, during stroke reversals. It contains three simplified aerodynamic assumptions: (i) unsteady potential flow; (ii) discrete vortices shed from both leading and trailing edges of the wing; (iii) Kutta conditions applied at both edges. Then the problem is reduced to the solution of the unsteady Laplace equation, by using distributed singularities, i.e., sources/sinks, and vortices in the field. To validate the present physical model and analytic method proposed via benchmark examples, two elemental motions in wing flapping and a case of whole flapping cycles are analyzed, and the predicted results agree well with available experimental and numerical data. This verifies that the present analytical approach may give qualitatively correct and quantitatively reasonable results. Furthermore, the total fluid-dynamic force in the present method can be decomposed into three parts: one due to the added inertial (or mass) effect, the other and the third due to the induction of vortices shed from the leading-and the trailing-edge and their images respectively, and this helps to reveal the flow control mechanisms in insect wing flapping. The project supported by the National Natural Science Foundation of China (10072066) and the Chinese Academy of Sciences (KJCX-SW-LO4, KJCX2-SW-L2)     

Effect of thrust on the aeroelastic instability of a composite swept wing with two engines in subsonic compressible flow

This paper aims to investigate aeroelastic stability boundary of subsonic wings under the effect of thrust of two engines. The wing structure is modeled as a tapered composite box-beam. Moreover, an indicial function based model is used to calculate the unsteady lift and moment distribution along the wing span in subsonic compressible flow. The two jet engines mounted on the wing are modeled as concentrated masses and the effect of thrust of each engine is applied as a follower force. Using Hamilton's principle along with Galerkin's method, the governing equations of motion are derived, then the obtained equations are solved in frequency domain using the K-method and the aeroelastic instability conditions are determined. The flutter analysis results of four example wings are compared with the experimental and analytical results in the literature and good agreements are achieved which validate the present model. Furthermore, based on several case studies on a reference wing, some attempts are performed to analyze the effect of thrust on the stability margin of the wing and some conclusions are outlined.     

基于当地流活塞理论的气动弹性计算方法研究

发展了一种高效、高精度的超音速、高超音速非定常气动力计算 方法------基于定常CFD技术的当地流活塞理论. 运用当地流活塞理论计算非定常 气动力，耦合结构运动方程，实现超音速、高超音速气动弹性的时域模拟. 运用这 种方法计算了一系列非定常气动力算例和颤振算例，并和原始活塞理论、非定 常Euler方程结果作了比较. 由于局部地使用活塞理论假设，这种方法大大地克服 了原始活塞理论对飞行马赫数、翼型厚度和飞行迎角的 限制. 与非定常Euler方程方法相比，当地流活塞理论的效率很高.     

低压涡轮内部流动及其气动设计研究进展

随着高空无人飞行器研究的不断升温, 高空低雷诺数条件下动力装置的研究越来越受到人们的重视.结合近年来国内外相关领域的研究工作, 对低雷诺数低压涡轮内部复杂流动机理的研究进展进行了介绍, 包括低雷诺数情况下低压涡轮内部非定常流动的特点, 叶片边界层分离及转捩现象机理, 上游周期性尾迹与下游叶片边界层相互作用机理等. 在此基础上给出了适合低雷诺数条件的低压涡轮气动设计方法:尾迹通过与边界层的相互作用, 能够抑制分离, 进而减小叶型损失, 在气动设计中有效引入非定常效应可以大幅度提高低压涡轮的气动负荷或降低气动损失, 最终达到提高性能的目的;数值及实验结果验证了这种设计方法的有效性.      

Aeroelastic analysis of helicopter rotor blade in hover using an efficient reduced-order aerodynamic model

This paper presents a coupled flap–lag–torsion aeroelastic stability analysis and response of a hingeless helicopter blade in the hovering flight condition. The boundary element method based on the wake eigenvalues is used for the prediction of unsteady airloads of the rotor blade. The aeroelastic equations of motion of the rotor blade are derived by Galerkin's method. To obtain the aeroelastic stability and response, the governing nonlinear equations of motion are linearized about the nonlinear steady equilibrium positions using small perturbation theory. The equilibrium deflections are calculated through the iterative Newton–Raphson method. Numerical results comprising steady equilibrium state deflections, aeroelastic eigenvalues and time history response about these states for a two-bladed rotor are presented, and some of them are compared with those obtained from a two-dimensional quasi-steady strip aerodynamic theory. Also, the effect of the number of aerodynamic eigenmodes is investigated. The results show that the three-dimensional aerodynamic formulation has considerable impact on the determination of both the equilibrium condition and lead-lag instability.     

基于RBF插值技术的CFD/CSD非线性耦合分析方法研究

基于非结构混合网格的N-S方程求解器和结构柔度影响系数法,发展了一种考虑气动、结构非线性的基于RBF插值技术CFD/CSD耦合分析方法,适用于解决现代大展弦比飞机的非线性静气动弹性问题。该方法采用时间相关法（即求解非定常方程组,用长时间的渐近解趋于定常状态）求解静气弹分析时的定常流动。考虑大展弦比飞机结构变形问题为大变形小应力问题,在利用柔度系数法求解结构方程时,假设每次求解结构方程时应力与应变为线性关系,整体静气弹分析过程为非线性关系,因此每次求解结构方程时要更新柔度影响系数矩阵。在非定常N-S方程每求解一个时间步耦合一次结构有限元分析,由于结构有限元分析的时间相对于气动分析时间是很短的,所以这种方法实际上近似使用了一次求解非定常气动力的时间完成了整个静气动弹性分析的过程。对于气动网格与结构有限元网格不一致性,本文采用径向基函数（RBF）插值方法中的TPS方法进行结构弹性变形和气动载荷插值,采用虚功原理完成气动载荷数据交换。为了节省气弹分析时间,采用动网格方法对气动网格进行更新,本文基于RBF插值方法发展一种适用于混合网格（四面体、三棱柱、金字塔和六面体）变形的动网格方法,可以保证附面层网格的质量与分布从而准确模拟其流动。利用该方法对M6机翼、DLR-F6翼身组合体和某大型客机机翼进行了静气动弹性特性分析,结果验证了本文开发的非线性CFD/CSD耦合分析方法的可行性、精确性和高效性。     

Modeling three-dimensional dynamic stall

The dynamic stall process in three-dimensional (3D) cases on a rectangular wing undergoing a constant rate ramp-up motion is introduced to provide a qualitative analysis about the onset and development of the stall phenomenon. Based on the enhanced understanding of the mechanism of dynamic stalls, a 3D dynamic stall model is constructed with the emphasis of the onset, the growth, and the convection of the dynamic stall vortex on the 3D wing surface. The results show that this engineering dynamic stall model can simulate the 3D unsteady aerodynamic performance appropriately.     

Separated Flow and Buffeting Control

In transonic flow conditions, the shock wave/turbulent boundary layer interaction and the flow separations on the upper wing surfaces of civil aircraft induce flow instabilities, ‘buffet’ and then structural vibrations, ‘buffeting’. Buffeting can greatly affect aerodynamic behavior. The buffeting phenomenon appears when the aircraft's Machnumber or angle of attack increases. This phenomenon limits the aircraft's flight envelope. The objectives of this study are to cancel out or decrease the aerodynamic instabilities (unsteady separation, movement of the shock position) due to this type of flow by using control systems. The following actuators can be used: ‘Vortex Generators’ situated upstream of the shock position, a ‘Bump’ located at the shock position, and a new moving part designed by ONERA, situated on the trailing edge of the wing, the ‘Trailing Edge Deflector’ or TED. It looks like an adjustable ‘Divergent Trailing Edge’. It is an active actuator and can take different deflections or be driven by dynamic movements up to 250 Hz. Tests were performed in transonic 2D flow with models well equipped with unsteady pressure transducers. For high lift coefficients, a selected static position of the ‘Trailing Edge Deflector’ increases the wing's aerodynamic performances and delays the onset of buffet. Furthermore, in 2D flow buffet conditions, the ‘Trailing Edge Deflector’, driven by a closed-loop active control using the measurements of the unsteady wall static pressures, can greatly reduce buffet. The aerodynamic performances are not improved to the same extent by the bump actuator. From our experience, there is no effect on buffet or separated flow because of the incorrect positioning of the bump. All that can be observed is a local improvement on the intensity of the shock wave when the bump is very precisely situated at the shock position. Vortex generators have a great impact on the separated flow. The separated flow instabilities are greatly reduced and buffet is totally controlled even for strong instabilities. The aerodynamic performances of the airfoil are also greatly improved.     

Aeroelastic stability of a cantilevered plate in yawed subsonic flow

The aeroelastic stability of cantilevered plates with their clamped edge oriented both parallel and normal to subsonic flow is a classical fluid–structure interaction problem. When the clamped edge is parallel to the flow the system loses stability in a coupled bending and torsion motion known as wing flutter. When the clamped edge is normal to the flow the instability is exclusively bending and is referred to as flapping flag flutter. This paper explores the stability of plates during the transition between these classic aeroelastic configurations. The aeroelastic model couples a classical beam structural model to a three-dimensional vortex lattice aerodynamic model. The aeroelastic stability is evaluated in the frequency domain and the flutter boundary is presented as the plate is rotated from the flapping flag to the wing configuration. The transition between the flag-like and wing-like instability is often abrupt and the yaw angle of the flow for the transition is dependent on the relative spacing of the first torsion and second bending natural frequencies. This paper also includes ground vibration and aeroelastic experiments carried out in the Duke University Wind Tunnel that confirm the theoretical predictions.