Numerical investigation of the effects of structural geometric and material nonlinearities on limit-cycle oscillation of a cropped delta wing

This article presents numerical simulations of the limit-cycle oscillation (LCO) of a cropped delta wing in order to investigate the effects of structural geometric and material nonlinearities on aeroelastic behavior. In the computational model, the structural part included both the geometric nonlinearity that arises from large deflections, and the material nonlinearity that originates from plasticity. The Euler equations were employed in the fluid part to describe the transonic aerodynamics. Moreover, the load transfer was conducted using a 3-D interpolating procedure, and the interfaces between the structural and aerodynamic domains were constructed in the form of an exact match. The flutter and LCO behaviors of the cropped delta wing were simulated using the coupling model, and the results were compared with existing experimental measurements. For lower dynamic pressures, the geometric nonlinearity provided the proper mechanism for the development of the LCO, and the numerical results correlated with the experimental values. For higher dynamic pressures, the material nonlinearity led to a rapid rise in the LCO amplitude, and the simulated varying trend was consistent with the experimental observation. This study demonstrated that the LCO of the cropped delta wing was not only closely related to geometric nonlinearity, but was also remarkably affected by material nonlinearity.     

EFFECTS OF STEADY ANGLE OF ATTACK ON NONLINEAR GUST RESPONSE OF A DELTA WING MODEL

The effects of a steady angle of attack on the nonlinear aeroelastic response of a delta wing model to a periodic gust have been studied. For the theoretical analysis, a three-dimensional time-domain vortex lattice aerodynamic model and a reduced order aerodynamic technique were used and the structure was modelled using von Karman plate theory that allows for geometric strain–displacement nonlinearities in the delta wing structure. Also, an experimental investigation has been carried out in the Duke wind tunnel using a rotating slotted cylinder gust generator and an Ometron VPI 4000 Scanning Laser Vibrometer measurement system to measure deflections (velocities) of a delta wing test model. The fair to good quantitative agreement between theory and experiment verifies that the present analytical approach has reasonable accuracy and good computational efficiency for nonlinear gust response analysis in the time-domain. The results also contribute to a better physical understanding of the nonlinear aeroelastic response of a delta wing model to gust loads when the steady angle of attack is varied.     

适用于几何非线性气动弹性分析的结构动力学降阶模型方法研究

几何非线性是壁板颤振和大展弦比机翼气动弹性等问题的一个主要特征,在进行数值仿真分析时往往需要采用商业非线性有限元求解器,存在计算量大和耦合迭代策略不易控制等问题。本文发展了一种适用于几何非线性的结构动力学降阶模型(CSD-ROM),利用广义坐标的非线性多项式表征非线性内力,采用参数识别方法获取多项式系数,并通过增加额外的线性模态来改善模型预测精度。基于此方法,分别针对壁板颤振、切尖三角翼的CFD/CSD-ROM非线性颤振问题开展了时域响应分析。计算结果表明,通过CSD-ROM计算出的壁板颤振速度为590 m/s,颤振频率为174 Hz,与有限元结果误差分别为0.8%和1.7%。马赫数0.879时切尖三角翼的颤振动压预测结果为2.25 psi,与非线性有限元相比的误差为3.8%。本文采用的非线性和线性模态基底组合方法,在保证计算精度的基础上可有效降低训练样本数量,一定程度上可替代非线性有限元开展气动弹性分析。     

Gust response: Simulation of an aeroelastic experiment by a fluid–structure interaction method

In this paper fluid–structure interaction simulations regarding a gust generator experiment are presented, which has been conducted in 2010 in the Transonic Wind Tunnel in Göttingen (DNW-TWG), Germany. The main objective of the experiment was the investigation of the dynamic response problem of an elastic wing model concerning an encountering generic gust induced by a gust generator. Fluid–structure simulations, using a finite element structural model and a computational fluid dynamics model based on time-accurate, Reynolds-averaged Navier–Stokes equations, are compared to the experiment to validate the numerical methodology. Comparisons include steady and unsteady deflections of the elastic wing and pressure distributions. Finally, the results of simulated transfer functions of the gust generator to the elastic wing are presented in comparison to the test data.     

BEM–FEM coupling for the analysis of flexible propellers in non-uniform flows and validation with full-scale measurements

The first part of the paper presents a partitioned fluid–structure interaction (FSI) coupling for the non-uniform flow hydro-elastic analysis of highly flexible propellers in cavitating and non-cavitating conditions. The chosen fluid model is a potential flow solved with a boundary element method (BEM). The structural sub-problem has been modelled with a finite element method (FEM). In the present method, the fully partitioned framework allows one to use another flow or structural solver. An important feature of the present method is the time periodic way of solving the FSI problem. In a time periodic coupling, the coupling iterations are not performed per time step but on a periodic level, which is necessary for the present BEM–FEM coupling, but can also offer an improved convergence rate compared to a time step coupled method. Thus, it allows to solve the structural problem in the frequency domain, meaning that any transients, which slow down the convergence process, are not computed. As proposed in the method, the structural equations of motion can be solved in modal space, which allows for a model reduction by involving only a limited number of mode shapes.The second part of the paper includes a validation study on full-scale. For the full-scale validation study a purposely designed composite propeller with a diameter of 1 m has been manufactured. Also an underwater measurement set-up including a stereo camera system, remote control of the optics and illumination system has been developed. The propeller design and the underwater measurement set-up are described in the paper. During sea trials blade deflections have been measured in three different positions. A comparison between measured and calculated torque shows that the measured torque is much larger than computed. This is attributed to the differences between effective and nominal wakefields, where the latter one has been used for the calculations. To correct for the differences between measured and computed torque the calculated pressures have been amplified accordingly. In that way the deformations which have been computed with the BEM–FEM coupling for non-uniform flows became very similar to the measured results.     

Navier–Stokes computation of transonic vortices over a round leading edge delta wing

A 3D Navier–Stokes solver has been developed to simulate laminar compressible flow over quadrilateral wings. The finite volume technique is employed for spatial discretization with a novel variant for the viscous fluxes. An explicit three-stage Runge–Kutta scheme is used for time integration, taking local time steps according to the linear stability condition derived for application to the Navier–Stokes equations. The code is applied to compute primary and secondary separation vortices at transonic speeds over a 65° swept delta wing with round leading edges and cropped tips. The results are compared with experimental data and Euler solutions, and Reynolds number effects are investigated.     

Numerical investigation of structural geometric nonlinearity effect in high-aspect-ratio wing using CFD/CSD coupled approach

An efficient and robust fluid–structure coupled methodology has been developed to investigate the linear and non-linear static aeroelastic behavior of flexible high-aspect-ratio wing. A three-dimensional open source finite element solver has been loosely coupled with an in-house Reynolds-averaged Navier–Stokes solver, designed for hybrid-unstructured meshes, to perform aero-structural coupled simulations. For volume mesh deformation and two-way data interpolation over non-matching grids interface, a radial basis function methodology combined with a data reduction algorithm has been used. This technique is efficient in handling large deflections and provides high-quality deformed meshes. Structural geometric nonlinearity has been considered to predict the deformations in the vertical and torsional directions caused by gravitational and aerodynamic loading. A multi-material finite element model has been generated to match the experimental configuration. Computational aeroelastic simulations were performed on an experimental high-aspect-ratio aeroelastic wing model with a slender body at the tip to get non-linear static deflections, twist and structure natural frequencies. The effect of the geometric nonlinearity is significant for large deformation analysis and has been highlighted in the predicted maximum tip deflection and twist. Good qualitative and quantitative agreement has been achieved between the predicted results and the available experimental data.     

DEVELOPMENT OF A THREE-DIMENSIONAL VISCOUS AEROELASTIC SOLVER FOR NONLINEAR PANEL FLUTTER

A new three-dimensional (3-D) viscous aeroelastic solver for nonlinear panel flutter is developed in this paper. A well-validated full Navier–Stokes code is coupled with a finite-difference procedure for the von Karman plate equations. A subiteration strategy is employed to eliminate lagging errors between the fluid and structural solvers. This approach eliminates the need for the development of a specialized, tightly coupled algorithm for the fluid/structure interaction problem. The new computational scheme is applied to the solution of inviscid two-dimensional panel flutter problems for subsonic and supersonic Mach numbers. Supersonic results are shown to be consistent with the work of previous researchers. Multiple solutions at subsonic Mach numbers are discussed. Viscous effects are shown to raise the flutter dynamic pressure for the supersonic case. For the subsonic viscous case, a different type of flutter behavior occurs for the downward deflected solution with oscillations occurring about a mean deflected position of the panel. This flutter phenomenon results from a true fluid/structure interaction between the flexible panel and the viscous flow above the surface. Initial computations have also been performed for inviscid, 3-D panel flutter for both supersonic and subsonic Mach numbers.     

Limit cycle oscillation behavior of transonic control surface buzz considering free-play nonlinearity

The limit cycle oscillation (LCO) behaviors of control surface buzz in transonic flow are studied. Euler equations are employed to obtain the unsteady aerodynamic forces for Type B and Type C buzz analyses, and an all-movable control surface model, a wing/control surface model and a three-dimensional wing with a full-span control surface are adopted in the study. Aerodynamic and structural describing functions are used to deal with aerodynamic and structural nonlinearities, respectively. Then the buzz speed and buzz frequency are obtained by V-g method. The LCO behavior of the transonic control surface buzz system with linear structure exhibits subcritical or supercritical bifurcation at different Mach numbers. For nonlinear structural model with a free-play nonlinearity in the control surface deflection stiffness, the double LCO phenomenon is observed in certain range of flutter speed. The free-play nonlinearity changes the stability of LCOs at small amplitudes and turns the unstable LCO into a stable one. The LCO behavior is dominated by the aerodynamic nonlinearity for the case with large control surface oscillation amplitude but by the structural nonlinearity for the case with small amplitude. Good agreements between LCO behaviors obtained by the present method and available experimental data show that our study may help to explain the experimental observation in wind tunnel tests and to understand the physical mechanism of transonic control surface buzz.     

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

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

Aeroheating CFD analysis of missile slots

This paper develops a hypersonic aerothermal simulation method for missile slot flow. The finite volume method of structure grid solver is developed for solving Euler and Navier-Stokes equations. The solver includes Park's two temperature model and the air multi-species reaction model. The second-order accuracy TVD numerical method was deduced to compute the hypersonic aeroheating which improves the computational efficiency. Computational results are given to show the high accuracy comparing to the existing experimental data.     

Non-linear flexural oscillations of a partially tapered annular plate

The free finite amplitude axisymmetric oscillations of an isotropic annular plate with partially tapered thickness are investigated. The time variable is eliminated by a Ritz-Kantorovich averaging method. The von Karman plate equations are then reduced to two non-linear ordinary differential equations, which form a non-linear eigenvalue problem. Solutions to the problem are obtained by utilizing a direct computational method. The results reveal the effects of large amplitude upon the dynamic responses. Also, an annulus of constant thickness, which has the same boundary conditions and the same volume as the partially tapered one, is investigated. Their results, which may shed light on the optimal design of annular plates, are compared.     

Application of the SAS turbulence model to predict the unsteady flow field behaviour in a forward centrifugal fan

In the current study, the unsteady flow in a centrifugal fan is carried out using Computational Fluid Dynamics calculation based on the Scale Adaptive Simulation (SAS) approach to model the turbulence phenomenon. The SAS concept is based on the introduction of the von Karman length scale into the turbulence scale equation. The information provided by the von Karman length scale allows SAS models to dynamically adjust to resolved structures in an Unsteady Reynolds-Averaged Navier–Stokes (URANS) simulation, which results in a Large Eddy Simulation-like behaviour in unsteady regions of the flow field. At the same time, the model provides standard RANS capabilities in stable flow regions. The introduction of the von Karman length scale is based on the reformulation of Rottas's equation for the integral length scale. To validate the numerical results, the overall performances of the fan and the wall pressure fluctuations computed upon the volute casing surface are compared with the unsteady measured data.     

Flow-induced vibration of curved cylinder arrays subject to loose support

In this paper, an improved model is developed aimed at analyzing the fluidelastic vibration of a single flexible curved tube which is surrounded by rigid cylinders and subjected to cross-flow and loose support. Based on the previous model, the axial extension of the curved tube described by von Karman nonlinearity has been accounted for in the current research. Simulations are performed to explore the effect of quasi-steady fluid force model and velocity-limited friction model on the post-instability behavior. Numerical results show that the out-of-plane response is confined to a limit cycle by the von Karman nonlinearity and the in-plane vibration is induced by the out-of-plane vibration through the nonlinear coupling. When the loose support comes into play, the nonlinear impact forces become dominate. The results are presented; comparisons are made to analyze the parameters influencing the fretting-wear damage, such as normal work rate, contact ratio and impact force level.     

Numerical simulation of transonic limit cycle oscillations using high-order low-diffusion schemes

This paper simulates the NLR7301 airfoil limit cycle oscillation (LCO) caused by fluid–structure interaction (FSI) using Reynolds averaged Navier–Stokes equations (RANS) coupled with Spalart–Allmaras (S–A) one-equation turbulence model. A low diffusion E-CUSP (LDE) scheme with 5th order weighted essentially nonoscillatory scheme (WENO) is employed to calculate the inviscid fluxes. A fully conservative 4th order central differencing is used for the viscous terms. A fully coupled fluid–structural interaction model is employed. For the case computed in this paper, the predicted LCO frequency, amplitudes, averaged lift and moment, all agree excellently with the experiment performed by Schewe et al. The solutions appear to have bifurcation and are dependent on the initial fields or initial perturbation. The developed computational fluid dynamics (CFD)/computational structure dynamics (CSD) simulation is able to capture the LCO with very small amplitudes measured in the experiment. This is attributed to the high order low diffusion schemes, fully coupled FSI model, and the turbulence model used. This research appears to be the first time that a numerical simulation of LCO matches the experiment. The simulation confirms several observations of the experiment.     

In-plane and out-of-plane dynamics of a curved pipe conveying pulsating fluid

This paper investigates the in-plane and out-of-plane dynamics of a curved pipe conveying fluid. Considering the extensibility, von Karman nonlinearity, and pulsating flow, the governing equations are derived by the Newtonian method. First, according to the modified inextensible theory, only the out-of-plane vibration is investigated based on a Galerkin method for discretizing the partial differential equations. The instability regions of combination parametric resonance and principal parametric resonance are determined by using the method of multiple scales (MMS). Parametric studies are also performed. Then the differential quadrature method (DQM) is adopted to discretize the complete pipe model and the nonlinear dynamic equations are carried out numerically with a fourth-order Runge–Kutta technique. The nonlinear dynamic responses are presented to validate the out-of-plane instability analysis and to demonstrate the influence of von Karman geometric nonlinearity. Further, some numerical results obtained in this work are compared with previous experimental results, showing the validity of the theoretical model developed in this paper.     

超音速气流中受热曲壁板的非线性颤振特性

基于von Karman 大变形理论及带有曲率修正的一阶活塞理论, 用Galerkin方法建立了超音速气流中受热二维曲壁板的非线性气动弹性运动方程; 采用牛顿迭代法计算得到由静气动载荷和热载荷引起的静气动弹性变形; 根据李雅谱诺夫间接法分析了壁板初始曲率与温升对颤振边界的影响; 对二维曲壁板的非线性气动弹性方程组进行数值积分求解,分析了动压参数对受热二维曲壁板分岔特性的影响, 给出了典型状态下曲壁板非线性颤振响应的时程图与相图. 分析结果表明对小初始曲率的曲壁板, 温升对其静气动弹性变形影响较大, 且随着温升的增加其颤振临界动压急剧减小; 对具有较大初始曲率的曲壁板, 温升对其静气动弹性变形的影响较弱, 且随着温升的增加颤振临界动压基本保持不变. 初始几何曲率与气动热效应使得曲壁板具有复杂的动力学特性, 不再像平壁板一样, 经过倍周期分岔进入混沌, 而会出现由静变形状态直接进入混沌运动的现象, 且在混沌运动区域中还会出现静态稳定点或谐波运动, 在大曲率情况下, 曲壁板不会产生混沌运动, 而是幅值在一定范围内的极限带振荡.      

A theoretical and experimental investigation of the effects of a steady angle of attack on the nonlinear flutter of a delta wing plate model

Limit cycle oscillations (LCO) of wings on certain modern high performance aircraft have been observed in flight and in wind tunnel experiments. Whether the physical mechanism that gives rise to this behavior is a fluid or structural nonlinearity or both is still uncertain. It has been shown that an aeroelastic theoretical model with only a structural nonlinearity can predict accurately the limit cycle behavior at low subsonic flow for a plate-like wing at zero angle of attack. Changes in the limit cycle and flutter behavior as the angle of attack is varied have also been observed in flight. It has been suggested that this sensitivity to angle of attack is due to a fluid nonlinearity. In this investigation, we study the flutter and limit cycle behavior of a wing in low subsonic flow at small steady angles of attack. Experimental results are compared to those predicted using an aeroelastic theoretical model with only a structural nonlinearity. Results from both experiment and theory show a change in flutter speed as the steady angle of attack is varied. Also the LCO magnitude increased at a given velocity as the angle of attack was increased for both the experiment and theory. While not proving that the observed sensitivity to angle of attack of LCO in aircraft is due to a structural nonlinearity, the results do show that a change in the aeroelastic behavior at angles of attack can be caused by a structural nonlinearity as well as a fluid nonlinearity. In this paper, only structural nonlinearities are considered, but an extension to include aerodynamic nonlinearities would be very worthwhile.     

基于聚合多重网格方法的跨音速颤振的数值模拟

采用流固耦合方法对跨音速颤振进行了数值模拟。流体方面在非结构网格上用有限体积方法求解了Euler方程;结构方面则求解了后掠机翼典型剖面的结构模态方程。时间推进采用双时间步长：对每一真实时间步,都通过基于聚合多重网格方法的伪时间步推进,对流体和结构方程交替迭代．得到一个稳态的流固耦合的解。文章最后给出了NACA64A010翼型剖面的跨音速颤振边界．与相关文献的计算结果符合良好。