Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Transonic Aeroelastic Numerical Simulation in Aeronautical Engineering

A lower-upper symmetric Gauss-Seidel (LU-SGS) subiteration scheme is constructed for time-marching of the fluid equations. The Harten-Lax-van Leer-Einfeldt-Wada (HLLEW) scheme is used for the spatial discretization. The same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Through subiteration between the fluid and structural equations, a fully implicit aeroelastic solver is obtained for the numerical simulation of fluid/structure interaction. To improve the ability for application to complex configurations, a multiblock grid is used for the flow field calculation and transfinite interpolation (TFI) is employed for the adaptive moving grid deformation. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between the fluid and structure. The developed code was first validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. Then, the flutter character of a tail wing with control surface was analyzed. Finally, flutter boundaries of a complex aircraft configuration were predicted.

Numerical Analyses of Transonic Fluid/Structure interaction in Aircraft Engineering

Through the coupling between aerodynamic and structural governing equations, a fully implicit multiblock aeroelastic solver was developed for transonic fluid/stricture interaction. The Navier-Stokes fluid equations are solved based on LU-SGS (lower-upper symmetric Gauss-Seidel) Time-marching subiteration scheme and HLLEW (Harten-Lax-van Leer-Einfeldt-Wada) spacing discretization scheme and the same subiteration formulation is applied directly to the structural equations of motion in generalized coordinates. Transfinite interpolation (TFI) is used for the grid deformation of blocks neighboring the flexible surfaces. The infinite plate spline (IPS) and the principal of virtual work are utilized for the data transformation between fluid and structure. The developed code was fort validated through the comparison of experimental and computational results for the AGARD 445.6 standard aeroelastic wing. In the subsonic and transonic range, the calculated flutter speeds and frequencies agree well with experimental data, however, in the supersonic range, the present calculation overpredicts the experimental flutter points similar to other computations. Then the flutter character of a complete aircraft configuration is analyzed through the calculation of the change of structural stiffness. Finally, the phenomenon of aileron buzz is simulated for the weakened model of a supersonic transport wing/body model at Mach numbers of 0.98 and l.05. The calculated unsteady flow shows, on the upper surface, the shock wave becomes stronger as the aileron deflects downward, and the flow behaves just contrary on the lower surface of the wing. Corresponding to general theoretical analysis, the flow instability referred to as aileron buzz is induced by a stronger shock alternately moving on the upper and lower surfaces of wing. For the rigid structural model, the flow is stable at all calculated Mach numbers as observed in experiment

Improvement of Aircraft Rolling Power by Use of Piezoelectric Actuators

Piezoelectric actuators are distributed on both side of a rectangular wing model,and the possibility of improvement of aircraft rolling power is investigated. The difference between the model with aileron deflection and the model without aileron (fictitious control surface, FCS) is studied. The analytical results show that these two cases are substantial different. In aileron deflection case, the aeroelastic effect is disadvantageous, so the structural stiffness should be high until the electrical voltage is not necessary. But in the case of FCS,the aeroelastic effect is advantageous and it means that lower structural stiffness can lead to lower voltage. Compared with aileron project, the FCS project can save structure weight.

飞行器跨声速气动弹性数值分析

将流体和结构运动方程分别构造为含子迭代的计算格式，发展了一种紧耦合气动弹性分析方法．其中流体计算的空间离散采用改进的HLLEW（Harten-Lax-van Leer-Einfeldt-Wada）格式．TFI（transfinite interpolation）方法用于生成随结构变形的自适应多块动网格．利用所发展的方法，对一翼-身-尾气动外形，数值预测了马赫数在0．3-1．3范围内的气动颤振边界．并详细研究了时间步长、子迭代步数、初始流场、耦合方法、疏密网格对颤振计算结果的影响．

采用分布式压电驱动器升力面的颤振主动抑制

对采用分布式压电驱动器升力面的颤振主动抑制进行了理论与试验研究。应用LQG最优控制法设计了主动控制律，在控制律降阶时提出了平衡实现与LK法结合使用的新途径，在对不定常气动力进行有理函数拟合时对LS法进行了改进。试验中利用激光测速仪非接触测量模型的速度响应并在地面共振试验中用压电驱动器激振模型。颤振风洞试验结果表明，理论计算合理并与试验结果吻合良好。

静气动弹性计算方法研究

对基于结构网格的Euler方程及N-S方程求解器和基于非结构网格的Euler方程求解器,采用结构模态分析方法和柔度矩阵方法,对无人机大展弦比机翼在Ma=0.6,α=2°,飞行高度20 km的巡航状态下的静气动弹性特性进行了数值模拟.验证了两种求解器对静气动弹性模拟的准确性.同时,对模态分析方法和柔度矩阵方法进行了对比研究,发现柔度矩阵方法更适用于静气动弹性数值模拟.另外,对应用物面法向偏转方法替代网格变形技术模拟静气动弹性进行了研究,计算表明物面法向偏转方法可以大大提高静气动弹性计算效率和克服机翼结构变形过大时动网格技术无法处理的不足.

压电驱动器的气动弹性应用

随着压电智能材料与结构的发展,压电驱动器在气动弹性控制领域占据重要地位.使用压电驱动器控制翼面变形,利用而不是抵抗气动弹性效应可以控制升力、力矩以及它们的分布.采用基本相同的智能结构翼面控制系统,根据不同的控制目标需求,使用压电智能材料驱动器可以达到多种目的,包括静态的形状控制与动态的颤振抑制、抖振控制与阵风响应控制.静态控制方面例如改变翼面形状获得附加空气动力以增加升力、提供横滚力矩、改变升力分布以减小诱导阻力或减小翼根弯矩等;动态控制例如利用改变翼面形状产生的附加空气动力作为控制载荷,改变气动弹性系统的耦合程度,根据控制效果要求可作为气动阻尼、气动刚度或气动质量.这种控制方法可以减轻结构重量,提高操纵效率,扩大飞行包线,提高材料利用率,已成为可变形飞行器的重要研究内容.本文主要阐述压电驱动器气动弹性应用的动机与机理、发展与成就以及问题与展望.