直升机动态失速研究

本文简要介绍了旋翼和翼型的动态失速特性、动态失速主动和被动控制、旋翼综合法预测使用的一些动态失速模拟方法.UH-60A直升机飞行试验11029飞行是一个典型的旋翼动态失速实例,表演出旋翼后行桨叶上出现高攻角引起的动态失速,而前行桨叶上则出现因为激波诱导前缘分离引起的动态失速.翼型的动态失速特性方面介绍了典型的动态失速迟滞回线(攻角正弦变化和锯齿变化)、缩减频率和雷诺数对翼型动态失速特性的影响等.动态失速主动控制介绍可动前缘翼型的情况,被动控制则是使用前缘涡发生器的情况.     

Numerical investigation of dynamic stall suppression of rotor airfoil via improved co-flow jet

The decrease in aerodynamic performance caused by the shock-induced dynamic stall of an advancing blade and the dynamic stall of a retreating blade at low speed and high angles of attack limits the flight speed of a helicopter. However, little research has been carried on the flow control methods employed to suppress both the dynamic stall induced by a shock wave and the dynamic stall occurring at high angles of attack. The dynamic stall suppression of a rotor airfoil by Co-Flow Jet(CFJ) is nume...     

改进型CLOR桨尖旋翼气动特性试验研究及数值分析

 通过风洞试验及数值模拟对具有改进型CLOR(CLOR-Ⅱ)桨尖的旋翼悬停和前飞状态气动特性开展研究。在CLOR桨尖旋翼试验及数值分析的基础上,考虑旋翼非定常流场特点,兼顾旋翼悬停和前飞气动性能,对旋翼桨叶的气动外形进行了改进,主要包括采用多种翼型优化配置以综合改善旋翼前行侧压缩性及后行侧桨叶失速特性,并考虑旋翼前飞状态对其桨叶动力学特性的需求,重新设计了桨尖前后缘的外形。在风洞中分别对3种旋翼进行多种状态条件下的试验研究,为从流动细节上获得不同桨尖旋翼的气动特性差别,采用计算流体力学(CFD)方法对试验状态进行了数值模拟对比。对更高转速状态进行模拟,结果表明相对于其他两种旋翼,CLOR-Ⅱ桨尖旋翼在改善跨声速特性和提高失速迎角等方面具有明显优势,而且综合提高了旋翼悬停和前飞气动性能。     

旋翼非定常动态失速的三维效应分析

采用建立的高精度计算流体动力学(CFD)方法,针对旋翼非定常动态失速的三维（3D）效应特性进行研究。以Helishape 7AD旋翼为基准,开展三维效应对旋翼非定常气动特性的影响分析。研究了来流马赫数对旋翼翼型动态失速特性影响。在此基础上,针对三维情形旋翼动态失速非定常涡流动特性及诱导分离特征开展了数值分析,通过与二维情形对比表明:受旋翼旋转、轴向诱导速度等三维效应影响,旋翼桨叶剖面动态失速涡的产生、对流和脱落明显滞后于二维翼型情形,并且涡强度也更弱。越靠近桨叶内段,桨叶剖面非定常动态失速特性与二维旋翼翼型情形的差距越为明显。      

多片后缘小翼对直升机旋翼桨叶动态失速及桨毂振动载荷的控制

 减缓直升机后行桨叶动态失速发生、降低直升机桨毂振动载荷是提高直升机飞行速度、改进直升机飞行性能的重要途径。本文研究了直升机在高速高载情况下利用多片受控的桨叶后缘小翼对直升机的后行桨叶动态失速和桨毂振动载荷同时进行控制的有效方法。建立了弹性桨叶和后缘刚性小翼的结构动力学模型。桨叶剖面气动载荷采用Leishman-Beddoes 二维非定常动态失速模型计算,后缘小翼剖面气动载荷采用Hariharan-Leishman二维亚声速非定常气动模型计算。采用伽辽金和数值积分相结合的方法求解旋翼系统的气弹响应。建立了有效的多片后缘小翼控制策略和控制方法,分析了3片后缘小翼的运动规律及对后行桨叶动态失速和桨毂振动载荷的控制效果,结果表明利用多片小翼的运动是控制桨叶动态失速和桨毂振动载荷的有效方法。     

高前进比旋翼气动特性分析方法研究

提出一种前进比高达0.8的旋翼气动特性分析方法,该方法针对高前进比旋翼前行桨叶压缩性、后行桨叶失速效应严重以及桨叶偏流作用和反流区大的特点,建立了旋翼气动力模型以及与之相适应的旋翼诱导速度时变非均布模型与桨叶非定常挥舞运动模型,然后根据高前进比旋翼气动力、旋翼诱导速度和桨叶挥舞运动三者之间的内在耦合关系提出了高前进比旋翼气动特性的动态响应计算方法,最后以H-34旋翼为例计算了该旋翼高前进比状态的气动特性,并将计算结果与风洞试验数据进行对比验证,结果合理。     

旋翼桨叶非定常挥舞运动的分析计算方法

从旋翼桨叶的实际气动环境出发导出了桨叶的非定常挥舞运动方程，其中考虑了影响旋翼桨叶非定常挥舞运动的各种因素，包括旋翼入流的时滞效应及前行桨叶的压缩性和后行桨叶的失速特性。给出了计算桨叶非定常挥舞运动的递推公式，并进而描述了旋的整体挥舞运动。最后以某型直升机为例，计算了直升机悬停及前飞中旋加变距操纵后旋翼桨叶挥舞运动响应的过渡过程，结果表明：用本文所述方法能合理地计算旋翼桨叶挥舞运动的非定常过渡过程     

柔性自适应壁风洞的翼型实验技术

给出一种新的二维跨音速柔性自适应壁风洞实验迭代方案,计算方法和实验验证结果。根据实验时实验段上下壁和模型上下表面实测压强分布对风洞内、外流场进行非线性数值模拟;计算流线化壁面的形状,进行自适应实验。用该迭代方案,在堵塞比ε=8％,实验段高与翼型弦长比值H/c=1.5情况下,对NACA-0012翼型进行了高亚音速验证性实验。实验结果与国外大风洞无干扰实验结果吻合很好。实验时迭代次数仅需1～2次。实验结果展示了自适应壁风洞实验技术用于翼型跨音速实验的前景。     

Parametric analyses for synthetic jet control on separation and stall over rotor airfoil

Numerical simulations are performed to investigate the effects of synthetic jet control on separation and stall over rotor airfoils. The preconditioned and unsteady Reynolds-averaged Navier–Stokes equations coupled with a k x shear stream transport turbulence model are employed to accomplish the flowfield simulation of rotor airfoils under jet control. Additionally,a velocity boundary condition modeled by a sinusoidal function is developed to fulfill the perturbation effect of periodic jets. The validity of the present CFD procedure is evaluated by the simulated results of an isolated synthetic jet and the jet control case for airfoil NACA0015. Then, parametric analyses are conducted specifically for an OA213 rotor airfoil to investigate the effects of jet parameters(forcing frequency, jet location and momentum coefficient, jet direction, and distribution of jet arrays) on the control effect of the aerodynamic characteristics of a rotor airfoil. Preliminary results indicate that the efficiency of jet control can be improved with specific frequencies(the best lift-drag ratio at F+= 2.0) and jet angles(40 or 75) when the jets are located near the separation point of the rotor airfoil. Furthermore, as a result of a suitable combination of jet arrays, the lift coefficient of the airfoil can be improved by nearly 100%, and the corresponding drag coefficient decreased by26.5% in comparison with the single point control case.     

Flapping and flexible wings for biological and micro air vehicles

Micro air vehicles (MAVs) with wing spans of 15 cm or less, and flight speed of 30–60 kph are of interest for military and civilian applications. There are two prominent features of MAV flight: (i) low Reynolds number (10⁴–10⁵), resulting in unfavorable aerodynamic conditions to support controlled flight, and (ii) small physical dimensions, resulting in certain favorable scaling characteristics including structural strength, reduced stall speed, and low inertia. Based on observations of biological flight vehicles, it appears that wing motion and flexible airfoils are two key attributes for flight at low Reynolds number. The small size of MAVs corresponds in nature to small birds, which do not glide like large birds, but instead flap with considerable change of wing shape during a single flapping cycle. With flapping and flexible wings, birds overcome the deteriorating aerodynamic performance under steady flow conditions by employing unsteady mechanisms. In this article, we review both biological and aeronautical literatures to present salient features relevant to MAVs. We first summarize scaling laws of biological and micro air vehicles involving wing span, wing loading, vehicle mass, cruising speed, flapping frequency, and power. Next we discuss kinematics of flapping wings and aerodynamic models for analyzing lift, drag and power. Then we present issues related to low Reynolds number flows and airfoil shape selection. Recent work on flexible structures capable of adjusting the airfoil shape in response to freestream variations is also discussed.     

Parametric analyses on dynamic stall control of rotor airfoil via synthetic jet

The effects of synthetic jet control on unsteady dynamic stall over rotor airfoil are investigated numerically. A moving-embedded grid method and an Unsteady Reynolds Averaged Navier-Stokes(URANS) solver coupled with k-x Shear Stress Transport(SST) turbulence model are established for predicting the complex flowfields of oscillatory airfoil under jet control. Additionally, a velocity boundary condition modeled by sinusoidal function has been developed to fulfill the perturbation effect of periodic jet. The validity of present CFD method is evaluated by comparisons of the calculated results of baseline dynamic stall case for rotor airfoil and jet control case for VR-7 B airfoil with experimental data. Then, parametric analyses are conducted emphatically for an OA212 rotor airfoil to investigate the effects of jet control parameters(jet location, dimensionless frequency, momentum coefficient, jet angle, jet type and dual-jet) on dynamic stall characteristics of rotor airfoil. It is demonstrated by the calculated results that efficiency of jet control could be improved with specific momentum coefficient and jet angle when the jet is located near separation point of rotor airfoil. Furthermore, the dual-jet could improve control efficiency more obviously on dynamic stall of rotor airfoil with respect to the unique jet, and the influence laws of dual-jet's angles and momentum coefficients on control effects are similar to those of the unique jet. Finally,unsteady aerodynamic characteristics of rotor via synthetic jet which is located on the upper surface of rotor blade in forward flight are calculated, and as a result, the aerodynamic characteristics of rotor are improved compared with the baseline. The results indicate that synthetic jet has the capability in improving aerodynamic characteristics of rotor.     

具有两个测量控制面的适应壁风洞理论与实验

本文研究了具有两个测量控制面的适应壁风洞。在超临界马赫条件下,做了二维NACA-0012翼型实验。在近洞壁处的两个控制面上测取了静压。基于上述测量,叙述了在适应壁风洞中得到无干扰流场的迭代方法。给出了为外场计算这两个控制面上静压的函数关系。评述了所选取迭代方法的收敛性。得到了一步收敛公式,并且在数字模拟风洞中得到证实。结果指出,应用一步迭代公式,所选取的迭代方法可以加速使流场收敛到无界条件。     

Experimental investigation of shock-buffet criteria on a pitching airfoil

An experimental investigation of the shock-buffet phenomenon subject to unsteady pitching supercritical airfoil around its quarter chord has been conducted in a transonic wind tunnel. The model was equipped with pressure taps connected to the fast response pressuretransducers. Measurements were conducted at different free-stream Mach number from 0.61 to0.76. The principle goal of this investigation was to experimentally discuss the shock-buffet criterion over a SC(2)-0410 supercritical pitching ...     

高亚声速翼型柔壁自适应壁风洞实验技术

本文通过解跨声速小振动方程求解假想外场，给出进行自适应壁风洞实验的迭代方案。结果表明：二者具有较好的一致性。本文研究表明：两种方案均有很好的应用前景。     

A methodology for simulating 2D shock-induced dynamic stall at flight test-based fluctuating freestream

A comprehensive methodology for simulating 2 D dynamic stall at fluctuating freestream is proposed in this paper.2 D CFD simulation of a SC1095 airfoil exposed to a fluctuating freestream of Mach number 0.537 ± 0.205 and Reynolds number 6.1 × 10~6(based on the mean Mach number) and undergoing a 10° ± 10° pitch oscillation with a frequency of 4.25 Hz was conducted.These conditions were selected to be representative of the flow experienced by a helicopter rotor airfoil section in a real-life fast forward flight.Both constant freestream dynamic stall as well as fluctuating freestream dynamic stall simulations were conducted and compared.The methodology was carefully validated with experimental data for both transonic flow and dynamic stall under fluctuating freestream.Overall, the results suggest that the fluctuating freestream alters the dynamic stall mechanism documented for constant freestream in a major way, emphasizing that inclusion of this effect in the prediction of dynamic stall related rotor loads is imperative for rotor performance analysis and blades design.     

直升机机动飞行的逆模拟

本文给出了一种直升机机动飞行的逆模拟方法以计算跟随预定飞行轨迹的驾驶员操纵,根据这一方法可以确定为完成直升机机动飞行所需的驾驶员操纵输入及直长机的飞行速度、角速度和的变化历程。直同飞行动力学模型没有作任何线化假设,其中考虑了旋翼入流的时滞效应、前行桨叶的压缩性物后行桨叶的失速特性及旋翼桨叶的非定常挥舞运动,引入了旋翼尾迹对直升机机身、尾翼和尾桨的气动干扰。最后以黑鹰直升机为例计算了鱼跃越障机动飞行     

Modern helicopter rotor aerodynamics

The helicopter rotor wake is among the most complex fluid dynamic structures being three dimensional and in many cases unsteady. The wake begins at the blade(s) where the flow can be transonic near the blade tip and undergo compressible dynamic stall. Farther down in the wake, the flow is essentially incompressible. Moreover, the rotor blades undergo complex unsteady motions because of the necessity to balance moments; they are elastic as well. In this paper, the fundamental aeromechanics of the wake and the flow on the blade is discussed and the primary methods of analysis, computation, and experiment employed to uncover the physics of the rotor wake are described.     

翼型前缘变形对动态失速效应影响的数值计算

翼型或机翼的动态失速效应所引起的低头力矩和正气动阻尼限制了飞行器气动性能的提高,甚至可能诱导发生不稳定运动。应用于小尺寸机翼的前缘动态变形（DDLE）技术,通过实时改变前缘形状,能够改善翼型前缘区域的速度梯度,进而抑制动态失速效应。采用转捩剪切应力输运（SST）黏性模型结合分区混合动态网格技术,研究了这种前缘变形对机翼俯仰运动所引起的非定常流动的影响,得到通过小幅度前缘变形抑制和延迟动态失速的方法,从而提高翼型的气动性能。翼型NAC A0012的数值模拟结果与动态失速风洞试验结果比较表明：所使用的数值计算方法能够较为准确地模拟翼型在动态失速过程中升力系数与俯仰力矩系数的变化情况,可用于研究前缘变形对翼型俯仰运动所引起的非定常流动的影响。前缘动态变形翼型俯仰运动过程的非定常流场的数值模拟表明：在大迎角下不同幅度的前缘下垂运动能够抑制流动分离的发生,从而抑制动态失速,但在大迎角下小幅度高频率的前缘下垂变形能更高效地抑制动态失速;前缘变形幅度以及变形沿中弧线的分布对升力系数和俯仰力矩系数的影响并不明显。     

现代超临界翼型设计及其风洞试验

开展了现代超临界翼型的设计研究,对现役飞机的压力分布形态进行了分析,针对现役飞机在巡航状态和阻力发散点的压力分布进行对比,提取了现役飞机超临界剖面设计的要点。采用类函数/型函数变换(CST)参数化方法、基于二阶震荡及自然选择的随机权重混合粒子群算法(RwSecSelPSO)、雷诺平均Navier-Stokes(RANS)方程、Kriging代理模型结合定期望值型的目标函数建立了优化设计系统。针对提高阻力发散马赫数和降低巡航低头力矩的设计指标,利用优化设计系统通过调整目标期望值设计了一系列满足设计指标但阻力发散马赫数不同的超临界翼型,并选择了其中具有典型特性的翼型进行了对比分析,验证了提高阻力发散马赫数和低速失速特性的设计方法,指出了在阻力发散点形成平顶形压力分布的超临界翼型具有较好的综合性能。对设计的超临界翼型进行了高、低速风洞试验验证,试验结果表明:设计结果达到了设计指标要求,提出的低速改进方案有效,层流对超临界翼型失速特性影响较大。     

Aerodynamic shape optimization for alleviating dynamic stall characteristics of helicopter rotor airfoil

In order to alleviate the dynamic stall effects in helicopter rotor, the sequential quadratic programming(SQP) method is employed to optimize the characteristics of airfoil under dynamic stall conditions based on the SC1095 airfoil. The geometry of airfoil is parameterized by the class-shape-transformation(CST) method, and the C-topology body-fitted mesh is then automatically generated around the airfoil by solving the Poisson equations. Based on the grid generation technology, the unsteady Reynolds-averaged Navier-Stokes(RANS) equations are chosen as the governing equations for predicting airfoil flow field and the highly-efficient implicit scheme of lower–upper symmetric Gauss–Seidel(LU-SGS) is adopted for temporal discretization. To capture the dynamic stall phenomenon of the rotor more accurately, the Spalart–Allmaras turbulence model is employed to close the RANS equations. The optimized airfoil with a larger leading edge radius and camber is obtained. The leading edge vortex and trailing edge separation of the optimized airfoil under unsteady conditions are obviously weakened, and the dynamic stall characteristics of optimized airfoil at different Mach numbers, reduced frequencies and angles of attack are also obviously improved compared with the baseline SC1095 airfoil. It is demonstrated that the optimized method is effective and the optimized airfoil is suitable as the helicopter rotor airfoil.