航空学报 > 2016, Vol. 37 Issue (8): 2539-2551   doi: 10.7527/S1000-6893.2015.0335

基于CFD和混合配平算法的直升机旋翼地面效应模拟

朱明勇, 招启军, 王博   

  1. 南京航空航天大学 直升机旋翼动力学国家级重点实验室, 南京 210016
  • 收稿日期:2015-09-03 修回日期:2015-12-04 出版日期:2016-08-15 发布日期:2012-08-11
  • 通讯作者: 招启军,Tel.:025-84893753。E-mail:zhaoqijun@nuaa.edu.cn E-mail:zhaoqijun@nuaa.edu.cn
  • 作者简介:朱明勇,男,硕士研究生。主要研究方向:直升机空气动力学和直升机计算流体力学。E-mail:zhumingyong@nuaa.edu.cn;招启军,男,博士,教授,博士生导师。主要研究方向:直升机计算流体力学、直升机空气动力学及流动控制、气动噪声、总体设计。Tel:025-84893753。E-mail:zhaoqijun@nuaa.edu.cn;王博,男,博士,讲师。主要研究方向:直升机计算流体力学、直升机空气动力学、新概念旋翼飞行器气动设计等。Tel:025-84893753。E-mail:wangbo@nuaa.edu.cn
  • 基金资助:

    国家自然科学基金(11272150)

Simulation of helicopter rotor in ground effect based on CFD method and hybrid trim algorithm

ZHU Mingyong, ZHAO Qijun, WANG Bo   

  1. National Key Laboratory of Science and Technology on Rotorcraft Aeromechanics, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
  • Received:2015-09-03 Revised:2015-12-04 Online:2016-08-15 Published:2012-08-11
  • Supported by:

    National Natural Science Foundation of China (11272150)

摘要:

建立了一套基于非结构网格技术和动量源模型的直升机旋翼计算流体力学(CFD)方法,用来模拟贴地飞行时直升机旋翼非定常气动特性。其中,控制方程采用三维Navier-Stokes方程,空间方向上采用Jameson格式,时间方向上采用五步Runge-Kutta迭代法,选用Spalart-Allmaras湍流模型。旋翼对流场的作用采用动量源项模拟,为更真实地模拟地面效应(IGE)的作用,采用了“移动地面”的物面边界来代替常规的“固定地面”边界,并对旋翼附近及旋翼与地面之间的网格进行加密处理,以提高地面涡的捕捉精度。考虑实际飞行环境下旋翼的运动和操纵,在流场计算时考虑旋翼配平特性。其中,配平方程的旋翼气动力通过CFD方法和动量-叶素组合理论模型的耦合计算给出,为了提高配平方法的鲁棒性和效率,提出并建立了基于遗传算法/拟牛顿法的高效混合迭代算法。运用所建立的方法,首先,选用有试验结果可供对比的算例计算了地面效应作用下的旋翼拉力增益、功率变化,验证了计算方法的有效性,解决了涡流理论方法较难模拟的“小速度前飞旋翼需用功率突增”问题。然后,着重研究了UH-60A直升机旋翼在不同离地高度、不同前进比状态,旋翼需用功率、诱导速度、地面涡及旋翼操纵的变化规律。计算结果表明:地面涡出现在较小的前进比范围内,随前进比的增大,地面涡在纵向平面将顺来流方向移动,在轴向方位靠近地面方向移动,直至最后不断减弱消失。

关键词: 旋翼, 地面效应, 非定常气动特性, 地面涡, 配平分析, Navier-Stokes方程

Abstract:

A computational fluid dynamics (CFD) method based on the unstructured grid technique and momentum source method is developed aimming at simulating the unsteady aerodynamic characteristics of a rotor in ground effect (IGE). In this method, the three-dimensional Navier-Stokes equations have been taken as governing equations, the discretization of convective fluxes and the time marching are completed by Jameson scheme and the five-step Runge-Kutta iteration method respectively, and one equation Spalart-Allmaras turbulence model has been employed. The rotor is modeled as a distribution of momentum source to simulate the ground effect more realistically, the boundary of "moving ground" is used instead of the conventional "stationary ground", and the grids near the rotor plane and between the rotor and the ground are refined to capture the ground vortex more accurate. Considering the motion and control of the rotor in the actual flight environment, the rotor trim is taken into account in the simulation of the rotor flowfield. The aerodynamic forces of trim equations are given by coupling of CFD method and momentum-blade element theory model. In order to improve the robustness and efficiency of the trim method, the genetic algorithm/quasi Newton hybrid iterative algorithm is proposed and established. Firstly, the method is used to calculate the rotor thrust increment and rotor power required in ground effect. The calculated results are compared with the experimental results aimming at verifying the validity of the method, the problem of ramp increment of rotor power required at a low speed flight has been solved which is difficult for wake analysis methods. Then, the flowfields of UH-60A helicopter rotor in different flight heights different advance ratios have been calculated, meanwhile the difference of the rotor power required, induced velocity, ground vortex and rotor control are investigated. The calculated results show that the ground vortex appears at small advance ratio; with increase of the advance ratio, it moves along inflow direction in the longitudinal plane and moves close to the ground in axial direction, then, it becomes weak continuously and finally disappears.

Key words: rotor, in ground effect, unsteady aerodynamic characteristics, ground vortex, trim analyses, Navier-Stokes equations

中图分类号: