桥梁颤振气动导数在均匀流与紊流作用下的统一辨识理论研究

桥梁断面颤振气动导数的准确识别是桥梁空气动力稳定性研究的基础,一直是桥梁风工程研究领域的前沿课题。基于工况模态分析理论,本文首先导出了紊流条件下桥梁颤振气动导数辨识的完整方法;其次,将上述数据处理方式运用到由均匀流条件引起的自由衰减振动中,得出与紊流情况类似的结果。本文提出的方法能够用同样的辨识原理、同一个数据处理程序统一地对紊流和均匀流下的桥梁气动导数进行辨识;最后,利用数值仿真算例验证上述理论。     

CHAOTIC MOTIONS AND LIMIT CYCLE FLUTTER OF TWO-DIMENSIONAL WING IN SUPERSONIC FLOW

Based on the piston theory of supersonic flow and the energy method, the flutter motion equations of a two-dimensional wing with cubic stiffness in the pitching direction are established. The aeroelastic system contains both structural and aerodynamic nonlinearities. Hopf bifurcation theory is used to analyze the flutter speed of the system. The effects of system parameters on the flutter speed are studied. The 4th order Runge-Kutta method is used to calculate the stable limit cycle responses and chaotic motions of the aeroelastic system. Results show that the number and the stability of equilibrium points of the system vary with the increase of flow speed. Besides the simple limit cycle response of period 1, there are also period-doubling responses and chaotic motions in the flutter system. The route leading to chaos in the aeroelastic model used here is the period-doubling bifurcation. The chaotic motions in the system occur only when the flow speed is higher than the linear divergent speed and the initial condition is very small. Moreover, the flow speed regions in which the system behaves chaos axe very narrow.     

A general modal frequency-domain vortex lattice method for aeroelastic analyses

The generalized aerodynamic force (GAF) matrix is derived for the Unsteady Vortex Lattice Method (UVLM) without the assumption of out-of-plane dynamics. As a result, the approach naturally includes in-plane motion and forces unlike the doublet lattice method (DLM). The derived UVLM GAF is therefore applicable to industry-standard techniques for aeroelastic stability analyses, such as the p–k method. In this work, the fluid–structure interpolation is performed with radial basis functions for surface interpolation. The generalized aerodynamic forces computed with the UVLM are verified against the DLM from NASTRAN on a simple flat plate configuration. The ability of the UVLM to include steady loads is verified with a T-tail flutter case and the results confirm the importance of including steady loads for T-tail flutter analysis. The modal frequency domain VLM therefore provides the same level of efficiency and accuracy than the DLM, but without the restrictions and with the ability to handle complex geometries. It is therefore a viable replacement to the DLM.     

用“参数多项式方法”确定风力透平叶片的失速振荡流

用“参数多项式方法”确定风力透平叶片的失速振荡流陈佐一,孙永忠,杨玲（清华大学热能工程系北京１０００８４）关键词：参数多项式方法,风力透平,颤振在现代各类流动诱导振动问题的研究中,流体激振的安全性分析与诊断,要求在较短的时间内,对大量的复杂流动工况进...     

A theoretical computational model of a plate in hypersonic flow

A theoretical model of an elastic panel in hypersonic flow is derived to be used for design and analysis. The nonlinear von Kármán plate equations are coupled with 1st order Piston Theory and linearized at the nonlinear steady-state deformation due to static pressure differential and thermal loads. Eigenvalue analysis is applied to determine the system’s stability, natural frequencies and mode shapes. Numerically time marching the equations provides transient response prediction which can be used to estimate limit cycle oscillation amplitude, frequency and time to onset. The model’s predictive capability is assessed by comparison to an experiment conducted at a free stream flow of Mach 6. Good agreement is shown between the theoretical and experimental natural frequencies and mode shapes of the fluid–structure system. Stability analysis is performed using linear and nonlinear methods to plot stability, flutter and buckling zones on a free stream static pressure vs temperature differential plane.     

Numerical bifurcation analysis of static stall of airfoil and dynamic stall under unsteady perturbation

By the finite element method combined with Arbitrary-Lagrangian-Eulerian (ALE) frame and explicit Characteristic Based Split Scheme (CBS), the complex flows around stationary and sinusoidal pitching airfoil are studied numerically. In particular, the static and dynamic stalls are analyzed in detail, and the natures of the static stall of NACA0012 airfoil are given from viewpoint of bifurcations. Following the bifurcation in Map, the static stall is proved to be the result from saddle-node bifurcation which involves both the hysteresis and jumping phenomena, by introducing a Map and its Floquet multiplier, which is constructed in the numerical simulation of flow field and related to the lift of the airfoil. Further, because the saddle-node bifurcation is sensitive to imperfection or perturbation, the airfoil is then subjected to a perturbation which is a kind of sinusoidal pitching oscillation, and the flow structure and aerodynamic performance are studied numerically. The results show that the large-scale flow separation at the static stall on the airfoil surface can be removed or delayed feasibly, and the ensuing lift could be enhanced significantly and also the stalling incidence could be delayed effectively. As a conclusion, it can be drawn that the proper external excitation can be considered as a powerful control strategy for the stall. As an unsteady aerodynamic behavior of high angle of attack, the dynamic stall can be investigated from viewpoint of nonlinear dynamics, and there exists a rich variety of nonlinear phenomena, which are related to the lift enhancement and drag reduction.     

机翼大攻角下失速颤振的气动弹性研究

本文采用流固耦合的数值方法研究了机翼在0°～50°攻角下的颤振。计算结果表明,随着来流攻角α0的增大,机翼的固有频率对颤振的影响越来越大,颤振由线性的强迫振动逐渐发展成为非线性的自激振动,而当α0增加到一定程度以后,大尺度分离流交替地从机翼头部和尾部产生,机翼和流场会发生共振,引起机翼的失速颤振。     

Control of a turbulent flow separated at mid-chord along an airfoil with DBD plasma actuators

The goal of the present experimental study is to investigate the ability of surface DBD plasma actuators to delay flow separation along the suction side of a NACA0015 airfoil. Three single surface DBD actuators that can operate separately are mounted on the suction side of the profile, at 18%, 27% and 37% of the chord length. The boundary layer is transitioned by a tripper to be sure that the flow control is not due to the laminar-to-turbulent transition. The angle of attack is equal to 11.5° and the free-stream velocity to U₀ = 40 m/s, resulting in a chord-based Reynolds number of Re_c = 1.33 × 10⁶. The flow is studied with a high-resolution PIV system. In such conditions, the baseline flow separation occurs at 50% of chord. Then, the different single DBD have been switched on separately, in order to investigate the actuator location effect. One highlights that the DBD located at x_c/c = 18% is more effective than the two others ones, with a separation delay up to 64% of chord. When the three DBDs operate simultaneously, the separation point moves progressively toward the trailing edge when the high voltage is increased, up to 76% of chord at 20 kV. Finally, the effect of the actuation frequency on the control authority has been investigated, by varying the value of the operating frequency and by burst-modulation. For frequencies equal to 50 Hz and 500 Hz (reduced frequency F⁺ = 0.31 and 3.1), the separation has been delayed at 76 and 80% of chord, respectively.     

Numerical simulation of flow over an airfoil in heavy rain via a two-way coupled Eulerian–Lagrangian approach

Airfoil performance degradation in heavy rain has attracted many aeronautical researchers’ eyes. In this work, a two-way momentum coupled Eulerian–Lagrangian approach is developed to study the aerodynamic performance of a NACA 0012 airfoil in heavy rain environment. Scaling laws are implemented for raindrop particles. A random walk dispersion approach is adopted to simulate raindrop dispersion due to turbulence in the airflow. Raindrop impacts, splashback and formed water film are modeled with the use of a thin liquid film model. The steady-state incompressible air flow field and the raindrop trajectory are calculated alternately through a curvilinear body-fitted grid surrounding the airfoil by incorporating an interphase momentum coupling term. Our simulation results of aerodynamic force coefficients agree well with the experimental results and show significant aerodynamic penalties at low angles of attack for the airfoil in heavy rain. An about 3° rain-induced increase in stall angle of attack is predicted. The loss of boundary momentum by raindrop splashback and the effective roughening of the airfoil surface due to an uneven water film are testified to account for the degradation of airfoil aerodynamic efficiency in heavy rain environment.