叶轮机周期性非定常流高效降阶数值计算方法研究进展

简要介绍了叶轮机内非定常流对叶轮机性能的重要影响,总结了叶轮机内动静相干和叶片振动诱发的非定常流的时空周期性特点,并指出这些时空周期性特点可以加以利用, 发展分析这些非定常流的高效降阶数值计算方法。本文根据不同数值降阶方法的特点对其进行了归类,详细分析了各自的优缺点及在气动弹性力学领域的应用,并进一步展望了降阶数值方法在未来的几个关键技术及发展趋势。     

预估蒸汽轮机叶片颤振的三维方法

发展了一种可考虑叶片间相角变化的三维失速颤振预测方法。二维计算的最小气动阻尼比三维值低, 说明取特征截面的二维方法偏于保守, 且特征截面的取法具有任意性, 故二维方法仅能作为叶片颤振初步预估之用。为提高颤振预估的精度, 需用三维计算对二维计算结果进行校核。三维和二维计算均表明叶片间相角对蒸汽轮机叶片气弹稳定性影响很大, 叶片颤振预测必须考虑这一关键因素的影响。      

Aeroelastic scaling laws for gust load alleviation control system

Gust load alleviation (GLA) tests are widely conducted to study the effectiveness of the control laws and methods. The physical parameters of models in these tests are aeroelastic scaled, while the scaling of GLA control system is always unreached. This paper concentrates on studying the scaling laws of GLA control system. Through theoretical demonstration, the scaling criterion of a classical PID control system has been come up and a scaling methodology is provided and veri-fied. By adopting the scaling laws in this paper, gust response of the scaled model could be directly related to the full-scale aircraft theoretically under both open-loop and closed-loop conditions. Also, the influences of different scaling choices of an important non-dimensional parameter, the Froude number, have been studied in this paper. Furthermore for practical application, a compen-sating method is given when the theoretical scaled actuators or sensors cannot be obtained. Also, the scaling laws of some non-linear elements in control system such as the rate and amplitude sat-urations in actuator have been studied and examined by a numerical simulation.     

Variations of flutter mechanism of a span-morphing wing involving rigid-body motions

The present paper aims to reveal the significance of rigid-body motions for the flutter mechanism of a span-morphing wing model. The inclusion of rigid-body motions into aeroelastic formulation and flutter analysis is presented. A state-space aeroelastic equation combining the dynamics of stepped Euler-Bernoulli beam with unsteady strip aerodynamic theory is developed by quasi-static modeling. Using a numerical example, variations of flutter mechanism from the bending-torsional flutter to the body-freedom flutter are observed as the span increases. In addition, effects of some dimensionless parameters on the variations of flutter mechanism are investigated. The investigated parameters belonging to the fuselage have limited influences on the bending-torsional flutter but a significant impact on the body-freedom flutter.     

Wing aeroelasticity analysis based on an integral boundary-layer method coupled with Euler solver

An interactive boundary-layer method, which solves the unsteady flow, is developed for aeroelastic computation in the time domain. The coupled method combines the Euler solver with the integral boundary-layer solver(Euler/BL) in a ‘‘semi-inverse" manner to compute flows with the inviscid and viscous interaction. Unsteady boundary conditions on moving surfaces are taken into account by utilizing the approximate small-perturbation method without moving the computational grids. The steady and unsteady flow calculations for the LANN wing are presented. The wing tip displacement of high Reynolds number aero-structural dynamics(HIRENASD) Project is simulated under different angles of attack. The flutter-boundary predictions for the AGARD445.6 wing are provided. The results of the interactive boundary-layer method are compared with those of the Euler method and experimental data. The study shows that viscous effects are significant for these cases and the further data analysis confirms the validity and practicability of the coupled method.     

直升机旋翼桨叶气弹优化减振设计方法

从振源着手通过设计参数的最优选择设计直升机旋翼桨叶,使传递到机身的交变载荷最小达到降低振动水平的目的是旋翼桨叶设计思想的进步。本文在简单概述该方法的基础上,以某4桨摆振柔软的无铰复合材料桨叶为研究对象,将其大梁模拟为一单闭室复合材料盒形梁,研究了通过铺层角的优化选择,降低4次/转的桨毂交变力与力矩的情况。数值算例表明方法效果良好。      

迷宫式密封装置气弹力的计算

采用射流理论和BLASIUS圆管理论, 导出了密封中各个动力系数的表达式, 使密封中各个动力系数的计算更加实用、简单。      

Flutter and gust response analysis of a wing model including geometric nonlinearities based on a modified structural ROM

In this paper, a framework is established for nonlinear flutter and gust response analyses based on an efficient Reduced Order Model (ROM). The proposed method can be used to solve the aeroelastic response problems of wings containing geometric nonlinearities. A structural modeling approach presented herein describes the stiffness nonlinearities with a modal formulation. Two orthogonal spanwise modes describe the foreshortening effects of the wing. Dynamic linearization of the ROM under nonlinear equilibrium states is applied to a nonlinear flutter analysis, and the fully nonlinear ROM coupled with the non-planar Unsteady Vortex Lattice Method (UVLM) is applied to gust response analysis. Furthermore, extended Precise Integration Method (PIM) ensures accuracy of the dynamic equation solutions. To demonstrate applicability and accuracy of the method presented, a wind tunnel test is conducted and good agreements between theoretical and test results of nonlinear flutter speed and gust response deflection are reached. The method described in this paper is suitable for predicting the nonlinear flutter speed and calculating the gust responses of a large-aspect-ratio wing in time domain. Meanwhile, the results derived highlight the effects of geometric nonlinearities obviously.     

Aeroelastic simulation of the first 1.5-stage aeroengine fan at rotating stall

A massive parallel aeroelastic simulation platform has been built to investigate the first 1.5-stage fan of an aeroengine at rotating stall. The Computational Fluid Dynamics (CFD) solver and Computational Structural Dynamics (CSD) solver are coupled directly by non-matching mesh interfaces. The unsteady rotor/stator interaction is solved by the Sliding Mesh Interface method. The original rotor blades are shrouded by the midspan shrouds. An unshrouded fan is also created to investigate the effects of the midspan shrouds. Both the shrouded fan and unshrouded fan have stable aeroelasticity at the designed state. At rotating stall, the stalled region rotates at 30% of the rotor speed on the absolute reference frame. The energy spectrum of the rotating stalled flow is measured quantitatively. It shows that the first two order excitations are much stronger than the higher order excitations. In the flow of rotating stall, the fifth backward travelling wave mode of shrouded fan is resonated by the fifth excitation of the rotational stalled flow because the rotational speed of the stalled region coincides with the modal rotational speed, while for the unshrouded fan, the first bending mode is resonated by the second excitation of the rotational stalled flow, forming two waves in the circumference of the annulus blades. At rotating stall, the vibration of the shrouded blades is still under control but the vibration of the unshrouded blades is diverged and out of control. A novel tool, i.e., resonance map, is proposed to predict the resonance. It provides a perspective to explain the effects of midspan shrouds at a theoretical level, and it would also be helpful in the structural design of blades.     

Transonic aeroelastic simulation for instability searches and uncertainty analysis

In this paper the use of eigenvalue stability analysis of very large dimension aeroelastic numerical models arising from the exploitation of computational fluid dynamics is reviewed. A formulation based on a block reduction of the system Jacobian proves powerful to allow various numerical algorithms to be exploited, including frequency domain solvers, reconstruction of a term describing the fluid-structure interaction from the sparse data which incurs the main computational cost, and sampling to place the expensive samples where they are most needed. The stability formulation also allows non-deterministic analysis to be carried out very efficiently through the use of an approximate Newton solver. Finally, the system eigenvectors are exploited to produce nonlinear and parameterised reduced order models for computing limit cycle responses. The performance of the methods is illustrated with results from a number of academic and large dimension aircraft test cases.