基于后缘小翼的旋翼翼型动态失速控制分析

针对后缘小翼(TEF)的典型运动参数对旋翼翼型动态失速特性的控制进行了研究。发展了一套适用于带有后缘小翼控制的旋翼翼型非定常流动特性模拟的高效、高精度CFD方法。通过求解Poisson方程生成围绕旋翼翼型的黏性贴体和正交网格,为保证后缘小翼附近的网格生成质量,建立了基于翼型点重构的方法来描述后缘小翼的偏转运动;为克服变形网格方法可能导致网格畸变的不足,发展了一套适用于带有后缘小翼控制的旋翼翼型运动嵌套网格方法。基于非定常雷诺平均Navier-Stokes(URANS)方程、双时间法、Spalart-Allmaras(S-A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,发展了旋翼翼型非定常气动特性分析的高精度数值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隐式时间推进方法及并行技术提高计算效率。以有试验结果验证的HH-02翼型和SC1095翼型为算例,精确捕捉了动态失速状态下的气动力迟滞效应,验证了本文方法的有效性。着重针对SC1095旋翼翼型的动态失速状态开展后缘小翼的控制分析,提出了可以体现翼型升力、阻力及力矩综合特性的关系式Po和Pc,揭示了后缘小翼振荡频率、相位差和偏转幅值对动态失速特性影响的规律。研究结果表明:当后缘小翼偏转的相对运动频率为1.0,且小翼运动规律与翼型振荡规律之间的相位差为0°时,后缘小翼能够更好地抑制翼型动态失速现象;在此状态下,当偏转幅值为10°时,SC1095翼型最大阻力系数和最大力矩系数可以分别降低19%和27%。     

旋翼非定常气动特性CFD模拟的通用运动嵌套网格方法

针对直升机旋翼非定常气动特性CFD模拟中的网格生成难题,提出了一套高效、通用的运动嵌套网格生成方法.首先,基于Poisson方程求解和翻折法生成旋翼桨叶的正交贴体网格.其次,针对旋翼桨叶的扭转分布及变距、挥舞等复杂运动,建立了一套通用的洞单元识别的扰动衍射法;为保证洞包络面的封闭性,完善了挖洞过程中网格加密策略;在洞边界确立基础上,提出了一种高效、鲁棒的最小距离法贡献单元搜索的改进方法.在此基础上,建立了基于RANS(Reynolds-averaged Navier-Stokes)方程的旋翼非定常流场CFD模拟方法.最后,采用所建立的方法分别对悬停和前飞状态下的C-T(Caradonna-Tung)旋翼和7A(Helishape 7A)旋翼的气动特性、桨尖涡的位置进行了计算,计算结果与试验值误差小于5%,验证了该运动嵌套网格生成方法在旋翼非定常气动特性CFD模拟中的有效性.      

基于Navier-Stokes方程的旋翼前飞状态气动特性数值模拟研究

建立了一个适用于旋翼前飞状态非定常流场和气动载荷计算的数值方法。在该方法中,控制方程为惯性坐标系下的三维非定常Navier-Stokes方程,空间方向上将三阶逆风格式(MUSCL)与通量差分分裂方法相结合。为较好地模拟流动分离,采用了一方程的Spalart-Allmaras湍流模型,时间方向则采用双时间法推进求解。为计入桨叶的旋转、周期性挥舞和变距运动,采用了运动嵌套网格方法。此外,为了更真实地反映旋翼桨叶的实际运动,发展了一个旋翼配平分析模型及求解方法。应用所建立的方法,首先计算了具有先进气动外形的四片桨叶的UH-60A直升机旋翼悬停状态的流场及气动性能,验证了该方法对悬停状态旋翼气动性能数值模拟的有效性。在此基础上,着重对Caradonna模型旋翼无升力前飞状态和AH-1G直升机旋翼桨-涡干扰前飞状态的气动特性(表面压强分布、扭矩、气动载荷)进行了计算,数值结果和试验值吻合良好,表明本文建立的CFD方法对旋翼前飞状态流场的模拟和气动载荷的计算是有效的。     

基于CFD方法的直升机旋翼/尾桨非定常气动干扰计算

采用CFD方法针对直升机旋翼/尾桨气动干扰问题进行了数值计算研究.建立了一个适用于旋翼/尾桨干扰气动力计算的方法.在该方法中,控制方程采用三维非定常雷诺平均Navier-Stokes(RANS)方程,网格系统使用嵌套网格技术以计及旋翼和尾桨的相对运动.基于所建立的方法,分别计算研究了悬停状态、低速前飞状态以及巡航平飞状态时干扰状态的旋翼和尾桨气动力,并深入分析了尾桨旋转方向对干扰特性的影响.计算结果表明:旋翼对尾桨的气动干扰在各计算状态都存在,且干扰影响不总是负面的,这取决于飞行状态;而尾桨对旋翼的干扰作用则在悬停状态时较为明显,但在前飞状态时可忽略不计;尾桨旋转方向对尾桨气动力的影响较大,合理地选择尾桨的旋转方向可以改善其气动特性.      

倾转四旋翼机气动特性分析及参数影响研究

采用机身/机翼/尾翼贴体网格与旋翼动量源相结合的CFD方法,开展了某倾转四旋翼机气动布局参数对气动特性的影响分析研究。首先,建立了一套适用于倾转四旋翼机全机干扰流场模拟的数值方法。然后,基于所建立方法开展了不同部件干扰情况下旋翼气动特性分析。最后,分析了前飞速度、旋翼旋向组合、旋翼横向距离对旋翼与机翼气动特性的影响规律。结果表明：后旋翼气动特性受全机干扰的影响较大,后旋翼拉力相对于孤立旋翼减小了41%。随着前飞速度的增加,前旋翼和机翼对后旋翼的气动干扰增强。从机头前方正视,机翼的左旋翼右旋和机翼的右旋翼左旋,有利于改善全机气动性能。前后旋翼的横向距离增大会提高后机翼和后旋翼的气动性能,横向距离为0.5 m时,后旋翼拉力比0 m时增大47%。     

基于CFD/CSD耦合方法的新型桨尖旋翼气动弹性载荷计算

为提高旋翼非定常气动弹性载荷的分析精度,在刚性旋翼计算流体力学(CFD)方法中引入计算结构动力学(CSD)方法,建立了一套适合于新型桨尖旋翼气动弹性载荷分析的CFD/CSD耦合方法。旋翼流场分析采用NavierStokes/Euler方程作为控制方程,围绕旋翼生成运动嵌套网格。在流场求解中,采用双时间法推进,通量计算采用Jameson中心格式,并采用B-L(Baldwin-Lomax)湍流模型。基于Hamilton变分原理和中等变形梁理论开展桨叶弹性运动变形分析,并发展了一套具有任意转角梁单元的新方法以提高新型桨尖旋翼的动力学分析精度。采用基于代数变换方法的网格变形策略,建立了一套CFD/CSD松耦合方法,桨叶运动变形和旋翼气动力信息通过流固交接面传递。首先分别对CSD和CFD模块进行了验证,然后计算了UH-60A旋翼在高速前飞状态下的气动弹性载荷,并与试验值进行了对比,最后重点对旋翼桨尖形状进行了参数分析。计算结果表明,相比于升力线理论和刚性旋翼CFD方法,CFD/CSD耦合方法可以显著提高旋翼非定常气动弹性载荷的分析精度,并能更准确地反映新型桨尖旋翼的气动弹性耦合效应;同时采用后掠桨尖在桨叶前行侧30°~90°方位角范围可以显著降低激波强度,有利于改善旋翼的气动特性。     

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

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

基于黏性涡粒子和非定常面元法耦合的倾转四旋翼飞行器气动干扰研究

为研究倾转四旋翼飞行器复杂的气动干扰问题，首先基于黏性涡粒子和非定常面元法建立了倾转四旋翼飞行器干扰流场的数值计算方法，并验证了该方法的有效性；其次，根据建立的方法对悬停和前飞状态下的气动特性和干扰情况进行了计算和分析；最后，针对干扰情况对倾转四旋翼的气动布局进行了部分优化设计。结果表明：悬停时的干扰主要来自旋翼下洗流对机翼作用产生的向下载荷，占旋翼拉力的25.3%;前飞时的干扰主要表现为旋翼尾流对机翼“增阻”作用，致使全机升阻比降低13.2%,增加前后旋翼的横向间距能使前飞时后旋翼的拉力和效率分别提升13.5%和4.2%,全机的升阻比提高3.1%;旋翼下方机翼后缘向下偏转的布局能使悬停时机翼受到的向下载荷减小40.7%。     

螺旋桨滑流对自转旋翼气动特性影响分析

自转旋翼机在前飞时螺旋桨滑流穿过桨盘平面,为研究螺旋桨滑流对自转旋翼的非定常气动干扰,基于RANS(雷诺平均Navier-Stokes)方程,采用运动嵌套网格方法建立了适用于自转旋翼-螺旋桨气动干扰流场的计算分析方法,并对模型进行模拟。分析低速情况下,孤立状态自转旋翼和组合状态自转旋翼非定常气动特性及流场特性,研究不同速度及螺旋桨位置对自转旋翼气动特性的影响。计算结果表明螺旋桨滑流会影响自转旋翼在各方位角的升阻力特性,并引起自转旋翼尾迹在0°方位角附近发生畸变;相同升力下,来流速度越大旋翼后倾角越小,螺旋桨滑流对自转旋翼影响越小;增大螺旋桨与自转旋翼间距可以减弱螺旋桨滑流对自转旋翼的气动干扰。     

基于主动后缘小翼控制的旋翼厚度噪声试验

基于声压相消原理,提出了一种基于主动后缘小翼控制的旋翼厚度噪声降噪方法。基于Ffowcs Williams-Hawkings(FW-H)方法给出了后缘小翼主动控制条件下旋翼气动噪声的计算公式;进行了试验总体方案设计,在全消声室开展了噪声主动抑制的开环试验(扫频、扫相和扫幅),并与数值计算结果进行了对比,验证了计算方法及试验的有效性;分析了后缘小翼偏转参数(幅度、频率、初始相位)对旋翼厚度噪声降噪的影响规律;分析了观测距离对厚度噪声降噪的影响。试验结果表明:不同频率的偏转规律均有一个最优的初始控制相位,合适的偏转规律相位能有效降低桨盘平面内观测点的厚度噪声2～3 dB。     

电控旋翼气动特性建模与风洞试验验证

首先建立了带襟翼翼型的非定常气动力模型,继而基于Peters-He广义动态尾迹理论,考虑襟翼偏转对电控旋翼叶素环境的影响,建立了电控旋翼有限状态尾迹模型;进一步基于Theodorsen理论推导出电控旋翼桨叶挥舞响应与桨叶变距和襟翼操纵量的关系,综合以上建立了电控旋翼气动特性分析模型.以改进型电控旋翼试验系统为平台进行了风洞试验,测量了不同风速、不同襟翼操纵条件下的电控旋翼气动力、桨距、襟翼偏角及旋翼挥舞角的变化情况.理论计算结果与试验数据符合情况良好,验证了所建立的分析模型的正确性,并得出以下结论:旋翼转速一定时,桨叶变距与襟翼操纵基本呈线性关系;旋翼拉力随襟翼总距的增加而逐渐减小,襟翼总距较大时,其实际气动效率略有下降;前飞状态时,襟翼总距操纵会引起桨叶的纵向周期变距.      

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.     

Panel/full-span free-wake coupled method for unsteady aerodynamics of helicopter rotor blade

A full-span free-wake method is coupled with an unsteady panel method to accurately predict the unsteady aerodynamics of helicopter rotor blades in hover and forward flight. The unsteady potential-based panel method is used to consider aerodynamics of finite thickness multi-bladed rotors, and the full-span free-wake method is applied to simulating dynamics of rotor wake. These methods are tightly coupled through trailing-edge Kutta condition and by converting doublet-wake panels to full-span vortex filaments. A velocity-field integration technique is also adopted to overcome singularity problem during the interaction between the rotor wake and blades. Helicopter rotors including Caradonna-Tung, UH-60A, and AH-1G rotors, are simulated in hover and forward flight to validate the accuracy of this approach. The predicted aerodynamic loads of rotor blades agree well with available measured data and computational fluid dynamics (CFD) results, and the unsteady dynamics of rotor wake is also well simulated. Compared to CFD, the present method obtains accurate results more efficiently and is suitable to rotorcraft aeroelastic analysis.     

GAW-1翼型前后缘变弯度气动性能研究

传统增升装置主要用于提高飞机起降气动性能。利用计算流体力学(CFD)的方法,引入了通用飞机翼型的前后缘变弯装置的概念,数值模拟了GAW-1翼型在爬升状态时,前缘变弯装置、后缘襟翼/副翼偏转以及前后缘装置综合偏转对翼型气动特性的影响。研究表明,前缘变弯装置可以有效地改善翼型的失速特性,失速迎角提高了3°左右,最大升力系数提高了4.56%;同时提高升阻比50%~120%;但在设计升力系数下,升力系数和阻力系数都略微减小。另一方面,后缘变弯装置可以改变最大升阻比所对应的迎角,以及在小迎角时,提高升力系数6%左右。翼型综合偏转可以在小迎角时增加升力系数,在大迎角时增加升阻比。     

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

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

旋翼气动力数值模拟准确度影响因素计算

针对旋翼气动力CFD计算时的一些常见影响因素的选择问题,进行了数值计算讨论.计算求解器选择自研代码,该代码采用基于共享内存(OMP)的并行化方法,在嵌套网格系统下求解Navier-Stokes方程,空间方向上使用2阶迎风格式,时间方向采用双时间方法.讨论的问题包括背景网格数目对悬停桨尖涡捕捉、对前飞状态激波捕捉和气动力的影响,远场动量理论修正作用对计算结果影响等.对于讨论的不同网格数目情况下,计算的悬停状态旋翼拉力差异在6%左右,扭矩差异2.6%左右,旋翼前飞阻力和扭矩最大差异为2.58%和3.28%.是否采用动量修正旋翼拉力相差3%~5%.通过计算,给出了针对以上因素计算结果差异的范围,可以为旋翼气动力计算工作提供一定的参考.      

倾转四旋翼飞行器垂直飞行状态气动特性

综合采用基于滑移网格技术的计算流体力学(CFD)方法与悬停状态气动干扰试验方法,对倾转四旋翼(QTR)飞行器垂直飞行状态的流场进行模拟与试验,研究飞行器垂直飞行状态气动特性以及部分参数对气动特性的影响.结果表明:倾转四旋翼飞行器在垂直飞行状态,前后旋翼之间干扰不明显,但旋翼与机翼的干扰明显;旋翼旋向对旋翼与机翼的干扰不...     

Robust Navier-Stokes method for predicting unsteady flowfield and aerodynamic characteristics of helicopter rotor

A robust unsteady rotor flowfield solver CLORNS code is established to predict the complex unsteady aerodynamic characteristics of rotor flowfield. In order to handle the difficult problem about grid generation around rotor with complex aerodynamic shape in this CFD code,a parameterized grid generated method is established, and the moving-embedded grids are constructed by several proposed universal methods. In this work, the unsteady Reynolds-Averaged Navier-Stokes(RANS) equations with Spalart-Allmaras are selected as the governing equations to predict the unsteady flowfield of helicopter rotor. The discretization of convective fluxes is accomplished by employing the second-order central difference scheme, third-order MUSCL-Roe scheme, and fifth-order WENO-Roe scheme. Aimed at simulating the unsteady aerodynamic characteristics of helicopter rotor, the dual-time scheme with implicit LU-SGS scheme is employed to accomplish the temporal discretization. In order to improve the computational efficiency of holecells and donor elements searching of the moving-embedded grid technology, the ‘‘disturbance diffraction method" and ‘‘minimum distance scheme of donor elements method" are established in this work. To improve the computational efficiency, Message Passing Interface(MPI) parallel method based on subdivision of grid, local preconditioning method and Full Approximation Storage(FAS) multi-grid method are combined in this code. By comparison of the numerical results simulated by CLORNS code with test data, it is illustrated that the present code could simulate the aerodynamic loads and aerodynamic noise characteristics of helicopter rotor accurately.