基于动态近似边界条件的气动弹性数值模拟

发展了一种利用欧拉方程计算非定常气动力的数值方法,通过在固定物面边界上满足动态近似边界条件计算出非定常气动力,避免了在每个时间步重新生成网格或需用动网格技术进行网格变形处理过程,提高了计算效率。运用这种方法计算了一系列非定常气动力算例,并与非结构动网格准确边界条件下的欧拉方程解和实验数据进行了比较,进一步分析了翼型俯仰角和马赫数对非定常气动力相对误差的影响。将气动力解算器与结构方程耦合进行气动弹性数值模拟,计算了跨音速具有S型颤振边界的二元气动弹性标准算例-Isogaiwing。算例结果表明,利用动态近似边界条件的欧拉方程具有简便、高效的特点,并能在小振幅情况下得到与精确边界条件精度相当的非定常流场解,还可以用于气动弹性分析。     

一种考虑攻角影响的超音速翼面颤振分析方法

采用活塞理论计算非定常气动力，用新的方法考虑了攻角的影响，和颤振运动方程耦合求解三维翼面的颤振临界条件。分别计算了四种翼型(平板翼，圆弧翼，六边形翼和四边形翼)的颤振临界速度，计算结果和风洞试验数据比较，二者相当一致。     

桥梁颤振气动导数识别的迭代法

本文提出一种利用自由振动响应通过年顿－拉夫逊迭代同时识别出桥梁全部8个颤振气动导的方法。此法具有对迭代初值要求不高，识别结果稳定的特点，并具有一定的抗噪声干扰能力。数字仿真与实物试验结果表明本文方法有效，可行。     

A CFD/CSD Model for Transonic Flutter

In this paper, a rapid deforming technique is developed to generate dynamic, three-dimensional, multi-block, mesh. The second-order Runge-Kutta time-marching method is used to solve the structural equations of motion. A dual-time method and finite volume discretization are applied for the unsteady Euler/Navier-Stokes equations to calculate the aerodynamic forces, in which the physical time step is synchronous with the structural equations of motion. The Spalart-Allmaras turbulence model is adopted for a turbulent flow. Due to mass dissimilarity, exiting in flutter calculations for a compressible flow, methods of variable mass and variable stiffness are developed to calculate the dynamic pressure of flutter at the point of mass similarity, and the flutter characteristics are then obtained in accordance with similarity rule. For completeness, the calculated transonic flutter characteristic results are presented and discussed for a double-wing and an aircraft model.     

前缘激波脱体的大迎角翼面颤振工程计算方法

本文利用二阶线化理论对当地流活塞理论作了修正 ,利用激波脱体的临界折角确定翼型上声速线位置 ,并用插值法估计脱体激波到声速线之间局部亚音速区的当地气流参数 ,声速线后的当地气流参数仍由膨胀波公式计算 ,从而形成了适合于激波脱体流态的小展弦比大迎角翼面颤振工程估算的一整套方法。如考虑翼前缘圆角的存在 ,全部Mn >1的小展弦翼面都可用本文方法统一进行颤振分析。经Ma =1.1～ 2 .5的带迎角舵模型风洞颤振试验验证 ,表明本文方法具有一定的精度 ,能满足工程设计要求。     

非结构运动网格下的三维机翼颤振数值分析

利用非结构运动网格技术,考虑三维机翼弹性变形求解三维跨音速非定常流场,同时耦合多自由度结构运动方程进行跨音速标模算例AGARD445.6机翼的颤振数值分析。采用的流动控制方程为三维非定常Euler方程。用中心有限体积法进行空间离散,用双时间隐式方法进行时间推进。数值分析获得的颤振临界速度与马赫数关系曲线明显反映了跨音速颤振临界速度随马赫数变化出现“凹坑”的物理特征。同时,数值分析结果与风洞实验结果一致性较好验证了所使用的方法。     

大展弦比机翼气动颤振的有限元分析

摘要：在Theodorsen二元气动力的基础上,建立非定常气动力时域内积分形式的表达式或者等价的频域表达式,利用粘弹性结构振动分析中对积分方程的等价变换将其写成与结构动力学方程一致的二阶常微分方程,将气动力的影响作为对结构有限元模型质量阵、刚度阵和阻尼阵的补充,保留了结构原有的所有动力学特性,并且能够直接用计算结构动力学的通用有限元软件进行空气-结构耦合的整体动力学分析,适合应用于具有复杂结构的气弹问题。气动力模型的建立可以利用各种试验及数值方法得到的气动力数据,适用性强。算例给出了大展弦比机翼的颤振边界计算结果。     

基于Euler方程的B型和C型嗡鸣特性数值研究

运用非定常Eu ler方程求解气动力,耦合结构运动方程,采用改进的四步龙格-库塔推进方法,在时间域内研究了音速附近发生的单自由颤振问题。其中气动力模块采用非结构网格技术进行空间离散,运用线性弹簧方法实现非结构网格的运动。空间上采用Jam son中心格式的有限体积法,时间上采用双时间推进方法。从无量纲化的结构方程出发,把M数、音速、密度、舵面刚度、转动惯量、弦长等一系列气动、结构参数对嗡鸣特性的影响归结为三个无量纲参数的影响。通过研究发现,运用非定常Eu ler方程可以模拟B型和C型嗡鸣,但不能模拟因激波附面层干扰而产生的A型嗡鸣。C型嗡鸣比B型嗡鸣发展缓慢,提高系统阻尼对于排除C型嗡鸣更有效。     

带有结构刚度非线性的跨音速翼型颤振特性研究

以非定常纳维尔－斯托克斯方程为主管方程,计算翼型振荡的瞬态非定常气动力,并与带有结构刚度非线性的颤振方程耦合求解,用时间推进的方法,计算了带有结构刚度非线性的结构响应特性,及带有三次型刚度和间隙型刚度非线性的跨音速翼型库振特性。计算研究表明,岂于同时具有结构和气动非线性,导致了振荡极限环极为复杂的特性。     

Application of proper orthogonal decomposition to the discrete Euler equations

The response of a fluid moving above a panel to localized oscillation of the panel is predicted using reduced‐order modelling (ROM) with the proper orthogonal decomposition technique. The flow is assumed to be inviscid and is modelled with the Euler equations. These non‐linear equations are discretized with a total‐variation diminishing algorithm and are projected onto an energy‐optimal subspace defined by an energy‐threshold criterion applied to a modal representation of time series data. Results are obtained for a bump oscillating in a Mach 1.2 flow. ROM is found to reduce the degrees of freedom necessary to simulate the flowfield by three orders of magnitude while preserving solution accuracy. Other observed benefits of ROM include increased allowable time step and robustness to variation of oscillation amplitude. Published in 2002 by John Wiley & Sons, Ltd.     

Control of the aerodynamic characteristics of wing airfoils by nonstationary energy supply in transonic flow

The possibility of controlling the aerodynamic characteristics of wing airfoils in transonic regimes of flight using one-sided pulse-periodic energy supply has been studied. The flow structure near the symmetric airfoil at different angles of attack and its aerodynamic characteristics as functions of the value of energy in its nonsymmetric (about the airfoil) supply have been determined by numerical solution of two-dimensional nonstationary gasdynamic equations. A comparison of the obtained results and the data of calculation of flow past the airfoil at different angles of attack without energy supply has been made. It has been established that a prescribed lift can be obtained, using energy supply, with a much higher fineness ratio of the airfoil than that in the case of flow past it at an angle of attack. Translated from Inzhenerno-Fizicheskii Zhurnal, Vol. 82, No. 1, pp. 18–22, January–February, 2009.