文章快速检索     高级检索
  空气动力学学报  2021, Vol. 39 Issue (6): 136-148  DOI: 10.7638/kqdlxxb-2021.0315

引用本文  

张卫国, 孙俊峰, 招启军, 等. 旋翼翼型气动设计与验证方法[J]. 空气动力学学报, 2021, 39(6): 136-148.
ZHANG W G, SUN J F, ZHAO Q J, et al. Aerodynamic design and verification methods of rotor airfoils[J]. Acta Aerodynamica Sinica, 2021, 39(6): 136-148.

基金项目

装备预先研究项目(41406040202,30103010304);基础和前沿技术研究基金(PJD20190002);风雷青年创新基金(PJD20190003)

作者简介

张卫国(1975-),男,山东兖州人,研究员,研究方向:直升机空气动力学. E-mail:zwglxy@163.com

文章历史

收稿日期:2021-09-18
修订日期:2021-11-02
优先出版时间:2021-12-25
旋翼翼型气动设计与验证方法
张卫国1,2 , 孙俊峰3 , 招启军4 , 武杰2 , 李国强2,5 , 马帅2 , 吴霖鑫2     
1. 西北工业大学 航空学院,西安 710000;
2. 中国空气动力研究与发展中心 低速空气动力研究所,绵阳 621000;
3. 中国空气动力研究与发展中心 计算空气动力研究所,绵阳 621000;
4. 南京航空航天大学 直升机旋翼动力学国家级重点实验室,南京 210016;
5. 国防科技大学 空天科学学院,长沙 410073
摘要:为满足我国直升机研发对自主旋翼翼型的迫切需求,开展了旋翼翼型气动优化设计与验证方法研究。发展了旋翼翼型指标分析与给定方法,给出了我国直升机旋翼翼型谱系规划设想,以进化多目标算法为基础结合PCA算法建立了旋翼翼型多点/多目标优化设计方法,突破了旋翼翼型气动特性精准测量风洞试验技术。利用所建立的方法对典型厚度翼型进行了优化设计并开展了风洞试验验证,计算评估与试验验证均表明自主设计翼型综合性能较国外参考翼型有一定提升。在此基础上,进一步构建了旋翼性能理论计算与试验验证综合评估方法,两种方式得到的旋翼主要性能数据偏差小于5%,表明该方法具有较高的可靠性。评估结果显示,基于设计翼型的旋翼模型气动性能较基于参考翼型的旋翼模型提升了3%。
关键词旋翼    翼型    指标给定    气动优化设计    风洞试验    洞壁干扰修正    
Aerodynamic design and verification methods of rotor airfoils
ZHANG Weiguo1,2 , SUN Junfeng3 , ZHAO Qijun4 , WU Jie2 , LI Guoqiang2,5 , MA Shuai2 , WU Linxin2     
1. School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China;
2. Low Speed Aerodynamics Institute of China Aerodynamics Research and Development Center, Mianyang 621000, China;
3. Computational Aerodynamics Institute of China Aerodynamics Research and Development Center, Mianyang 621000, China;
4. National Key Laboratory of Rotorcraft Aeromechanics, Nanjing University of Aeronautics & Astronautics, Nanjing 210016, China;
5. College of Aerospace Science and Engineering, National University of Defense Technology, Changsha 410073, China
Abstract: In order to meet the urgent needs of next-generation helicopters for independent-designed high-performance rotor airfoils, a study on the aerodynamic optimization design and verification method of rotor airfoils has been carried out. The analyze and setting methods of rotor airfoil performance indicators are developed, and the helicopter rotor airfoil pedigree planning is proposed as well. Based on the evolutionary multi-objective algorithm combined with the principal component analysis (PCA) algorithm, a multi-point/multi-objective optimization design method for rotor airfoils is established. The method is then used to optimize the design of airfoils with typical thicknesses. Aerodynamic performance of the optimized airfoils are verified by wind-tunnel experiments in which breakthrough in the accurate measurement of rotor airfoil aerodynamic characteristics has been achieved. Both computational and experimental results show that the comprehensive performance of the designed airfoil has an improvement to some extent compared with the reference one. On this basis, a comprehensive evaluation method for theoretical calculation and experimental verification of rotor performance is constructed. The deviation of the main performance indicators of the rotor obtained by the two methods is less than 5%, indicating that the method is highly reliable. The evaluation results show that the aerodynamic performance of the rotor model based on the designed airfoil is improved by 3% compared with the one based on the reference airfoil.
Keywords: rotor    airfoil    indicators setting    aerodynamic optimization design    wind-tunnel test    wall interference correction    
0 引 言

翼型是旋翼的基本构成要素,其气动特性优劣对旋翼性能有关键性影响,进而影响直升机的载重、速度、航程、噪声和振动水平等[1-2]

旋翼翼型设计技术是直升机设计中的一项核心技术[3-6]。20世纪70年代起,设计者们逐渐认识到[7],旋翼翼型宽马赫数、变迎角、非定常等特殊的气动运行环境决定了其设计较固定翼翼型设计更为复杂,同时需要综合考虑前飞、悬停、机动等多种飞行条件,设计重点为提高阻力发散马赫数、最大升力、减小力矩等。基于以上因素,国外研究者发展并建立了旋翼翼型高精度数值计算方法、多目标多约束设计方法、精细的风洞测试验证技术等,建立了从二维翼型设计[8-12]至三维旋翼评估[13-15]再至全机飞行验证的较为完善的研发体系[16-17],催生了以美国VR和SC系列、法国OA系列、俄罗斯TsAGI系列等为典型代表的高性能旋翼翼型成果,这些翼型在第三代和第四代直升机上得到了广泛应用[18-20]

旋翼翼型系列发展方向主要分为两种:一种是追求更高的升力和阻力发散马赫数,谱系化方向发展(见图1);另一种是根据任务需求定制,专用化方向发展;其中谱系化是专用化的基础。当前,国外传统布局旋翼翼型谱系已经发展到第五代,翼型性能的整体提升空间愈发受限,翼型专用化方向发展更占优势。进入21世纪,共轴刚性、倾转旋翼等高速直升机面临更为复杂的气动环境,必须需要先进翼型作为支撑,从而推动了新的旋翼翼型设计技术研发[20-22]。例如,X-2在桨叶内侧配置双钝头翼型,在桨叶中段使用尖后缘翼型,在桨尖区域采用薄翼型以提升旋翼气动效率。倾转旋翼机在设计旋翼桨叶时,桨根采用厚翼型而桨尖采用薄翼型,使桨叶具有较大的升力系数和较高的阻力发散马赫数。


图 1 典型旋翼翼型气动特性比较 Fig.1 Comparison of aerodynamic characteristics between typical rotor airfoils

国内旋翼翼型研发工作从20世纪后期才开始启动。“八五”、“九五”期间,西北工业大学等单位对法国OA212翼型进行了改进,形成了OA212MK翼型、OA212MKT翼型,但改进翼型存在力矩系数偏大的问题。此外,掌握的翼型数据对国外翼型开展了反设计工作,设计了HA翼型族HA12、HA09和HA07等,但是未形成完整的旋翼翼型设计、验证方法,也没有形成具有独立知识产权的旋翼专用翼型。“十一五”至“十三五”期间,中国空气动力研究与发展中心(下称气动中心)、中国直升机设计研究所、西北工业大学、南京航空航天大学等单位[23-25]持续开展旋翼翼型自主设计与试验验证工作,进行了旋翼翼型气动力特性和优化设计方法研究,在翼型谱系规划、指标给定、设计方法、风洞试验和计算校核、性能评估等方面开展了大量的基础性研究工作,设计、计算和风洞试验验证的能力大幅进步,研发了旋翼翼型气动设计与评估软件平台,旋翼翼型的自主设计研发能力得到了质的提升[26-36]

但与国外相比,我国仍未建立自主旋翼翼型谱系,专用化翼型更为欠缺,从而严重制约了自主翼型在大型直升机上的工程应用。究其根源在于我国旋翼翼型的研究多限于单点的技术攻关,而旋翼翼型研发是一项精细化的系统工程,需要开展深入的系统研究才能构建我国自主的旋翼翼型气动设计与验证技术体系。体系示意图见图2,主要包括6个环节、3个层次。近15年来,气动中心重点对前5个环节开展了相关研究工作,初步将翼型设计与验证技术体系拓展至第三层次。

本文给出了我国旋翼翼型谱系规划设想,基于多点/多目标降维优化方法,设计了整体性能优于国外参考系列翼型的自主旋翼翼型,建立了旋翼性能理论评估与试验综合验证方法,开展了基于设计翼型的旋翼模型验证,推动了我国旋翼翼型的自主研发及技术体系的建设。


图 2 旋翼翼型气动设计与验证技术体系 Fig.2 Aerodynamic design and verification system for rotor airfoils
1 旋翼翼型设计指标给定 1.1 我国直升机旋翼翼型谱系规划

直升机旋翼桨叶通常由多种厚度翼型构成(见图3),以适应旋翼桨叶复杂的气动环境,这是旋翼翼型具有多个系列的主要根源,由此推动了国外旋翼翼型谱系化及专用化发展。旋翼翼型的升级换代极大地促进了直升机的换代升级。


图 3 典型直升机型号旋翼桨叶翼型构成 Fig.3 Rotor blade airfoil configurations of typical helicopters

目前,我国直升机旋翼仍然以使用国外第三代旋翼翼型为主,所带来的问题是常规构型旋翼改型升级空间受限,高速构型无专用翼型可用,需要规划并逐步建立具有自主知识产权的旋翼翼型谱系。

我国直升机旋翼翼型谱系的发展首先需要适应我国直升机研制发展需要。在常规构型直升机方面,需要在提高承载性能和速度性能等方面发展,对应于高性能型和高速度型旋翼翼型;在高速构型直升机方面,侧重发展复合式高速直升机旋翼翼型,以及倾转旋翼机的翼型;在无人直升机方面,需要侧重其总体性能特点和新概念布局形式加以发展,如涵道式、电推进式、变形旋翼式、微型直升机等。

表1给出了我国旋翼翼型谱系规划设想,以CRA(China Rotor Airfoil首字母缩写)命名,旋翼翼型谱系规划的基础立足于国外第三代水平进行提升(即以3XX系列起步),以6%~25%厚度范围内旋翼翼型为主体,但并不局限于此范围。

表 1 我国旋翼翼型谱系规划 Table 1 Pedigree planning of China rotor airfoils
1.2 旋翼翼型设计指标给定方法

旋翼翼型设计指标的给定尤其重要和困难,是旋翼翼型设计的输入和约束,直接决定了设计翼型是否符合使用要求。通常的给定方法如下:首先,按照直升机设计总体任务要求,确定旋翼主要参数;其次,利用旋翼气动性能及气动环境评估方法研究直升机在不同飞行状态下旋翼升阻比随总距、周期变距等操纵的变化规律,得到旋翼在不同飞行状态下的最优化操纵策略;然后,在旋翼主要参数、参考或设计翼型、桨叶气动外形的约束下,以旋翼最优操纵策略为前提,进行不同飞行状态下的旋翼气动特性计算和翼型气动环境评估;最后,依据计算和评估结果判定旋翼性能是否满足要求,进而给定旋翼翼型指标。在此过程中,首先使用参考翼型进行旋翼评估,如果不满足则提出新翼型设计指标,其次对应用新翼型的旋翼进行再评估,如果不满足则修改翼型设计,直到满足旋翼的使用要求。其给定方法研究过程如图4所示。


图 4 旋翼翼型指标给定方法研究过程 Fig.4 Process of setting performance indicatorsfor rotor airfoils

通过旋翼翼型指标给定方法的计算及分析,旋翼翼型设计指标给定实例如表2所示。

表 2 旋翼翼型设计指标 Table 2 Design indicators for rotor airfoils
2 旋翼翼型气动优化设计 2.1 旋翼翼型气动特性数值模拟技术

数值模拟技术是旋翼翼型设计的基础,用于评估翼型气动特性优劣,目前主要采用基于RANS的计算方法,但单纯使用此方法存在两点不足:一是旋翼翼型气动特性(如阻力系数、阻力发散马赫数等)计算精度较低;二是旋翼翼型在大迎角流动下对翼型外形、迎角、雷诺数非常敏感,后缘分离和前缘吸力峰强烈耦合,RANS计算方法的可靠性不能完全满足要求。

为增加RANS方法的适用性,采用转捩模型来提高阻力计算精度,主要基于Menterhe和Langtry发展的γ-Reθ[36]模型进行转捩模型的改进。模型共包括四个微分输运方程,其本质是:通过求解雷诺动量厚度输运方程来判断湍流转捩位置,而后采用间歇因子输运方程来模拟湍流转捩长度,最后,在完全发展湍流阶段回归到SST湍流模式。间歇因子、转捩雷诺动量厚度的输运方程表达式为:

$ \frac{{\partial (\rho \gamma )}}{{\partial t}} + \frac{{\partial (\rho \gamma {u_j})}}{{\partial {x_j}}} = \frac{1}{{ R e}}\frac{\partial }{{\partial {x_j}}}[(\mu + {\sigma _\gamma }{\mu _i})\frac{{\partial \gamma }}{{\partial {x_j}}}] + {P_\gamma } - {E_\gamma } $ (1)
$\begin{split}& \frac{{\partial (\rho \mathop {Re}\limits^{\sim} {\kern 1pt} {_{\theta t}})}}{{\partial t}} + \frac{{\partial (\rho \mathop R\limits^{\sim} {\kern 1pt} {e_{\theta t}}{u_j})}}{{\partial {x_j}}} = \\&\qquad \frac{1}{{{Re}{\kern 1pt} }}\frac{\partial }{{\partial {x_j}}}[{\sigma _{\theta t}}(\mu + {\mu _t})\frac{{\partial \mathop {Re}\limits^{\sim} {\kern 1pt} {_{\theta t}}}}{{\partial {x_j}}}] + {P_{\theta t}} \end{split} $ (2)

湍流模型采用SA及SST模型,同时使用多重网格和预处理等技术缩短计算周期,提高计算精度。

对于大迎角气流分离计算,考虑到LES能够高保真度地模拟分离流动,发展了RANS/LES混合高精度数值模拟技术,开展了DES类方法构建研究(图5)。


图 5 DES 得到的瞬时Q等值面 Fig.5 Q isosurface obtained by DES

表3给出了大迎角情况下9%厚度旋翼翼型多种计算方法与试验数据的对比,综合来看,DES计算方法的数值模拟结果与试验数据吻合度最高,较其他计算方法,精度更高,可靠性也更高。

表 3 9%厚度旋翼翼型多种计算方法比较 Table 3 Comparison between computational and experimental results for a 9%-thickness rotor airfoil 
2.2 高维多目标降维优化设计方法

旋翼运行过程中面临复杂的非定常环境,旋翼翼型的设计需要在多种飞行下满足旋翼气动性能的要求,具有明显的高维多目标、强约束的特点。综合采用全局多目标进化算法、翼型参数化技术、代理模型技术及高精度CFD方法构建了旋翼翼型的多目标优化设计方法。基于方差分析和Pareto前沿技术相结合,实现了设计空间的动态扩展,提高了设计空间的完备性;采用鲁棒优化设计技术,提高了设计点与非设计点的综合性能;针对旋翼翼型高维优化设计的特点,采用基于PCA(主成分分析)方法的多目标降维技术,降低了问题的复杂性,提高了优化效率。图6给出了基于PCA分析的优化设计流程。


图 6 基于PCA分析的多目标优化设计流程图 Fig.6 Flow chart of the multi-objective optimization design based on PCA analysis

以12%厚度翼型作为参考翼型,开展多目标/多点优化设计,利用PCA方法实现多目标降维。问题可描述为:

(1)提高前飞性能:

$ MDD\;\;{{\rm{subject}}\;\;\;\;{\rm{to}}}:\left|{C}_{MDD}\right|\leqslant \left|{C}_{MDD0}\right|$

(2)提高机动性能:

$\begin{split}& Ma = 0.3,{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0},{K}_{\mathrm{max}}\geqslant {K}_{\mathrm{max},0}\\& Ma = 0.4,{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0},{K}_{\mathrm{max}}\geqslant {K}_{\mathrm{max},0}\\& Ma = 0.5,{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0},{K}_{\mathrm{max}}\geqslant {K}_{\mathrm{max},0}\\& Ma = 0.6,{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0},{K}_{\mathrm{max}}\geqslant {K}_{\mathrm{max},0}\end{split}$

(3)提高悬停性能:

$\begin{split}& K\geqslant {K}_{0},{\;\;{\rm{subject}}\;\;{\rm{to}}:\;\;}Ma = 0.5,C_{L}{ = }0.{6}\\& K\geqslant {K}_{0},\;\;{{\rm{subject}}\;\;{\rm{to}}:\;\;}Ma = 0.6,C_{L}{ = }0.{6}\end{split}$

(4)保持厚度不减:

$ {{\rm{subject}}\;\;{\rm{to}}:\;\;}{t}_{\mathrm{max}}\geqslant {t}_{\mathrm{max},0} \qquad\qquad\qquad$

对原始优化问题进行PCA分析降维后,得到简化优化问题:

$ \begin{split} & MDD\;\;{{\rm{subject}}\;\;{\rm{to}}}:\left|{C}_{MDD}\right|\leqslant \left|{C}_{MDD0}\right|\\& Ma = 0.3,\;\;{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0}\;\;\\& Ma = 0.4,\;\;{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0}\;\;\\& Ma = 0.5,\;\;{K}_{\mathrm{max}}\geqslant {K}_{\mathrm{max},0}\\& Ma = 0.6,\;\;{C}_{L,\mathrm{max}}\geqslant {C}_{L,\mathrm{max},0}\\& K\geqslant {K}_{0},{\;\;{\rm{subject}}\;\;{\rm{to}}:\;\;}Ma = 0.6,C_{L}\text{ = }0.\text{6}\\& {{\rm{subject}}\;\;{\rm{to}}:\;\;}{t}_{\mathrm{max}}\geqslant {t}_{\mathrm{max},0}\end{split} $

经过PCA分析,剩余了7个设计目标,剔除了原始优化问题中的5个冗余目标,降低了原始设计问题的复杂性,优化目标更加明确。

图7给出了优化前后翼型外形的比较,表4给出了优化前后翼型气动特性的比较。从结果可以看出,采用PCA降维后得到的优化设计翼型的综合性能较参考翼型有明显提高,优化设计方法具有较高的可靠性。


图 7 基准翼型和优化翼型外形比较 Fig.7 Comparison between the baseline and the optimized airfoil

表 4 12%厚度基准翼型和优化翼型性能比较 Table 4 Performance comparison between a 12%-thickness reference airfoil and an optimized airfoil
2.3 旋翼翼型气动优化设计

根据旋翼翼型指标要求,结合不同位置剖面翼型的性能需求,利用构建的旋翼翼型多点/多目标优化设计方法,开展自主旋翼翼型的优化设计,选择内段、中段以及桨尖3个典型位置的翼型进行设计,反复迭代后给出设计翼型的外形数据,其研究过程如图8


图 8 旋翼翼型气动设计方法优化研究过程 Fig.8 Optimization design process for rotor airfoils

选取典型厚度自主设计的CRA309、CRA312翼型进行重点分析,翼型外形曲线见图9图10给出了9%厚度翼型与参考翼型的性能计算对比曲线,自主设计的翼型性能皆有不同程度的提升。


图 9 CRA3XX旋翼翼型系列外形曲线 Fig.9 Geometries of CRA3XX rotor airfoil series


图 10 9%厚度翼型参考翼型和优化翼型不同马赫数下性能曲线比较 Fig.10 Performance comparisons between a 9%-thickness airfoil and an optimized airfoil under conditions with different Mach numbers
3 旋翼翼型气动特性风洞试验验证 3.1 旋翼翼型升阻特性精准测量技术

高速风洞中通过增压实现旋翼翼型高雷诺数试验,此种工况下翼型风洞试验的阻力测量精准度除与模型制造水平、来流速压(或马赫数)控制、洞内温度等方面密切相关外,更与尾流测量技术有关,这是获取高精准度阻力数据的关键。

一般通过提高设备性能、改进试验方法和优化数据处理程序等几个方面来提高升阻特性测量的精准度。基于移测式尾流测量技术、阻力精确积分方法、洞壁干扰修正方法[35]等关键技术来获得高可信度的风洞试验验证数据。将经过验证的高速风洞试验技术应用到典型厚度设计翼型模型高速风洞试验验证中,验证设计翼型的风洞试验结果是否满足翼型设计指标。其研究过程如图11所示。


图 11 旋翼翼型高速风洞试验验证研究过程 Fig.11 Verification process of rotor airfoils in high-speed wind tunnel test verification
3.2 高精准度洞壁干扰修正技术

本节对洞壁干扰修正方法进行阐述。以法国S3MA风洞的壁压信息法为基础,建立了具有数值解的上、下壁干扰效应修正方法。

旋翼翼型高速风洞洞壁干扰壁压信息法修正公式采用单参数线化壁压信息法。扰动速势 $ \phi $ 在风洞壁附近的区域中满足亚声速小扰动线化速势方程:

$ {\beta ^2}\frac{{{\partial ^2}\phi }}{{\partial {x^2}}} + \frac{{{\partial ^2}\phi }}{{\partial {y^2}}} = 0 $ (5)
$ \beta = \sqrt {1 - {Ma^2}} $ (6)
$ \phi = {\phi _m} + {\phi _b} $ (7)

其中, $ {\phi _m} $ 为模型在洞壁附近的远场扰动速势, $ {\phi _b} $ 为待求的洞壁扰动速势。

通过适当的工程计算方法可以得到修正的气动力系数。例如,翼型试验可以使用以下公式:

$ {C_{Lc}} = {C_{Lu}} + \frac{{\partial {C_L}}}{{\partial \alpha }}\Delta \alpha + \frac{{\partial {C_L}}}{{\partial Ma}}\Delta Ma $ (8)
$ {C_{mc}} = {C_{mu}} + \frac{{\partial {C_m}}}{{\partial \alpha }}\Delta \alpha + \frac{{\partial {C_m}}}{{\partial Ma}}\Delta Ma + {d}_{C_m} $ (9)
$ {C_{Dc}} = {C_{Du}} + \frac{{\partial {C_D}}}{{\partial \alpha }}\Delta \alpha + \frac{{\partial {C_D}}}{{\partial Ma}}\Delta Ma + {d}_{C_{Df}} + {d}_{C_{Di}} $ (10)

式中: ${d}_{C_m} $ 为由于模型区中△Ma、△α不均匀而造成的附加俯仰力矩, ${d}_{C_{Df}} $ 是模型区中△Ma不均匀引起的浮阻力修正量, ${d}_{C_{D}} $ 是△α引起的诱导阻力修正量。

3.3 高雷诺数下设计翼型模型验证分析

验证试验在气动中心FL-20连续式跨声速风洞中开展(如图12所示)。对设计的CRA309旋翼翼型缩比模型进行考核验证。图13给出了参考翼型OA309和设计翼型CRA309在Ma = 0.4、0.6时的气动系数对比图。由图可见,设计翼型CRA309在失速后升力系数仍保持在较高水平,且在常用迎角0°~10°之间的阻力系数明显小于参考翼型OA309,最大升阻比明显优于参考翼型。


图 12 旋翼翼型高速风洞验证试验 Fig.12 Verification test of rotor airfoils in a high-speed wind tunnel


图 13 OA309翼型与CRA309翼型气动力系数曲线对比图 Fig.13 Comparison of aerodynamic coefficients between OA309 and CRA309 airfoils
4 基于设计翼型的旋翼理论和试验评估 4.1 旋翼气动特性理论评估方法

翼型优化设计的最终目的是提升旋翼气动性能,因此,设计翼型是否符合要求,需要在旋翼气动环境下进行考核评估。本文开展了两种评估方法研究,包括旋翼气动特性快速评估方法、旋翼气动特性数值评估方法。

4.1.1 旋翼气动特性快速评估方法

基于动量理论、叶素理论、涡流理论以及旋翼桨叶的挥舞运动建模等,开展旋翼气动特性快速评估方法研究。根据动量理论和涡流理论分别计算诱导速度,然后进行诱导速度、气动力以及挥舞运动之间的配平迭代,最终给出旋翼的气动性能数据。通过该方法可以实现旋翼气动特性的快速评估,从而满足翼型/旋翼的快速设计迭代的需要。流程图见图14


图 14 旋翼气动特性快速评估方法 Fig.14 Rapid method for evaluating rotor aerodynamic characteristics

基于参考翼型(OA3系列)和自主设计翼型(CRA3系列)配置两幅旋翼桨叶,旋翼桨叶其他参数一致,仅翼型不同。采用旋翼气动特性快速评估方法进行悬停性能评估。图15给出了悬停状态两种模型桨叶的扭矩-拉力系数曲线以及悬停效率-拉力系数曲线的比较。计算结果显示,自主设计翼型构成的桨叶悬停性能优于参考翼型的,大拉力系数下优化翼型性能的优化提升效果要优于小拉力系数下的,最大可实现3%的提升。


图 15 旋翼气动特性快速评估方法理论评估结果 Fig.15 Theoretical evaluations of rotor performance
4.1.2 旋翼气动特性数值评估方法

为进一步真实模拟旋翼气动特性,立足于国内现有的数值计算方法,发展了结构/笛卡尔网格重叠、结构网格变形、笛卡尔网格自适应、桨叶动力学建模、大规模并行计算和CFD/CSD耦合求解等先进方法,并结合低耗散格式和先进湍流模型方法、旋翼配平方法等,开展了旋翼气动特性数值评估方法研究。

基于建立的旋翼气动特性数值评估方法,对由基准翼型和设计翼型配置的两副旋翼分别开展了气动特性计算,对比分析两副旋翼的悬停和前飞性能以及流场特性。

图16图17给出了100%转速即桨尖马赫数为0.65时,基准旋翼和设计旋翼的悬停性能对比曲线。随着拉力系数的增加,两副旋翼的扭矩系数逐渐增加,悬停效率增加至最大值后开始减小。旋翼二(设计旋翼)的气动性能略好于旋翼一(基准旋翼);旋翼一的最大悬停效率可达到0.765,旋翼二的最大悬停效率可达到0.774。优化翼型性能在大拉力系数下的优化提升效果要优于小拉力系数下的,最大可实现1%幅度的提升。值得注意的是,工程快速评估方法和数值评估方法得到的结果之间还存在一定的差异。建立既准确又快速的旋翼气动特性理论评估方法是下一步研究的重要方向。


图 16 旋翼扭矩系数随拉力系数变化曲线对比(100%转速) Fig.16 Variations of the rotor torque coefficient as a function of the tension coefficient (100% rotating speed)


图 17 旋翼悬停效率随拉力系数变化对比(100%转速) Fig.17 Variations of the hover coefficient as a function of the pulling force coefficient (100% rotating speed)
4.2 基于设计翼型的旋翼模型试验验证

在气动特性评估的基础上,对基于参考翼型(OA3系列)和自主设计翼型(CRA3系列)配置的两幅旋翼桨叶模型,开展悬停和前飞工况的试验考核。

验证试验在气动中心FL-13低速风洞中基于直升机全域试验台开展,如图18所示。对基于参考翼型和设计翼型的两幅旋翼的气动性能进行了考核分析,验证提升效果。


图 18 旋翼模型风洞试验 Fig.18 Wind-tunnel test of a rotor model

旋翼桨叶模型直径4 m,弦长0.125 m,桨尖形状为抛物线,图19给出了桨叶翼型配置示意图。


图 19 旋翼桨叶翼型配置 Fig.19 Configuration of a rotor blade airfoil

桨叶模型采用复合材料制成,共2幅,分别为基于OA3系列参考翼型的旋翼(简称“参考翼型旋翼”,见图20)和基于CRA系列设计翼型的旋翼(简称“设计翼型旋翼”,见图21)。


图 20 参考翼型旋翼模型实物图 Fig.20 Photo of reference airfoils


图 21 设计翼型旋翼模型实物图 Fig.21 Photo of optimized airfoils

图22给出了参考翼型旋翼和设计翼型旋翼悬停性能结果对比曲线。同等拉力系数条件下,设计翼型旋翼的悬停效率较参考翼型旋翼有一定提升。参考翼型旋翼的悬停效率最大为0.735,设计翼型旋翼的悬停效率最大达到了0.760。在常用拉力系数范围内(0.010~0.020),设计翼型旋翼的悬停效率较参考翼型旋翼提升了约3%,与旋翼气动特性快速评估获得的结果较为一致。


图 22 悬停状态旋翼气动性能对比 Fig.22 Comparisons of aerodynamicperformance of a hovering rotor

图23给出了在保持桨尖马赫数相似的前提下,垂向力系数Cw = 0.0119时前飞工况的功率系数和升阻比变化曲线。如图所示,随前进比μ的增加,功率系数和升阻比逐渐增大,μ = 0.25时设计翼型旋翼当量升阻比达到了7.86;相同工况下,与参考翼型旋翼相比,设计翼型旋翼的功率系数减小了约3%~5%,升阻比增大了约3%~5%。


图 23 前飞状态旋翼性能对比(Cw = 0.0119) Fig.23 Comparisons of rotor performance at the forward flight state
4.3 旋翼气动特性理论评估与试验数据综合分析

基于建立的旋翼气动特性理论评估方法和试验验证方法,开展设计翼型旋翼性能综合对比分析。图24给出了额定转速(n = 1050 r/min)下设计翼型旋翼悬停状态气动性能理论计算与风洞试验结果对比。可以看出,两种方法获得的曲线趋势吻合度较好,扭矩系数和旋翼效率的误差均在5%以内。但前飞工况下,理论评估和风洞试验结果的差异较大,仅典型工况下二者偏差可达到10%之内。


图 24 设计翼型旋翼悬停状态理论评估与试验结果对比曲线 Fig.24 Comparisons of aerodynamic performance of a hovering rotor airfoil obtained by computation and experiment
5 结论与展望

本文开展了旋翼翼型气动设计和验证方法研究,建立了旋翼翼型气动优化设计方法,开展了自主翼型设计。翼型验证试验表明,设计翼型相较于参考翼型,性能获得了一定提升。进一步,开展了从二维翼型设计到三维旋翼评估的探索性工作。理论评估和试验验证表明,应用自主设计翼型可获得旋翼气动性能的提升,可为我国旋翼翼型自主研发体系的建立和自主翼型的工程应用提供重要基础支撑。

1) 制定了我国直升机旋翼翼型谱系规划设想,开展了旋翼翼型多点/多目标强约束优化设计研究,设计的CRA系列翼型相对于参考翼型(OA系列翼型),在俯仰特性无显著变化情况下,最大升力系数、最大升阻比、阻力发散迎角都有显著提升。

2) 自主设计翼型旋翼的悬停效率较参考翼型旋翼提升了约3%;配平前飞状态下,与参考翼型相比,设计翼型旋翼气动性能提升了约3%~5%。

3) 设计翼型旋翼的悬停状态气动性能的理论评估与风洞试验结果对比,吻合较好,扭矩系数和旋翼效率的误差均小于5%。

为了真正建立自主研发体系和研发自主翼型谱系,今后还需要在以下方面持续开展工作:

1)需要针对研发体系的六个环节系统性地开展研究工作,加强和完善已有的技术环节,补齐全尺寸飞行空白环节,重视各环节之间的支撑和衔接,使六个环节形成一个有机的整体;

2)推动我国旋翼翼型谱系化和专用化研发,形成涵盖不同构型的旋翼翼型谱系化数据库,推动自主旋翼翼型的工程应用。同时开展特定要求的新概念旋翼飞行器旋翼翼型设计探索性工作,以此支撑概念布局创新研发;

3)努力寻求旋翼翼型设计思想和技术上的新突破。通过翼型和旋翼两种设计技术的结合,充分考虑翼型配置、变弦长、扭转规律、桨尖形状等旋翼复杂布局设计,以达到最大化旋翼性能提升的目的。

参考文献
[1]
王适存, 徐国华. 直升机旋翼空气动力学的发展[J]. 南京航空航天大学学报, 2001, 33(3): 203-211.
WANG S C, XU G H. Progress of helicopter rotor aerodynamics[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2001, 33(3): 203-211. DOI:10.3969/j.issn.1005-2615.2001.03.001 (in Chinese)
[2]
吴希明. 高速直升机发展现状、趋势与对策[J]. 南京航空航天大学学报, 2015, 47(2): 173-179.
WU X M. Current status, development trend and countermeasure for high-speed rotorcraft[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2015, 47(2): 173-179. (in Chinese)
[3]
PADULA S, LI W. Options for robust airfoil optimization under uncertainty[C]//Proc of the 9th AIAA/ISSMO Symposium on Multidisciplinary Analysis and Optimization, Atlanta, Georgia. Reston, Virigina: AIAA, 2002. doi: 10.2514/6.2002-5602
[4]
LI W, KRIST S, CAMPBELL R. Transonic airfoil shape optimization in preliminary design environment[J]. Journal of Aircraft, 2006, 43(3): 639-651. DOI:10.2514/1.13921
[5]
HICKS R M, HENNE P A. Wing design by numerical optimization[J]. Journal of Aircraft, 1978, 15(7): 407-412. DOI:10.2514/3.58379
[6]
吴希明, 牟晓伟. 直升机关键技术及未来发展与设想[J]. 空气动力学学报, 2021, 39(3): 1-10.
WU X M, MU X W. A perspective of the future development of key helicopter technologies[J]. Acta Aerodynamica Sinica, 2021, 39(3): 1-10. DOI:10.7638/kqdlxxb-2021.0042 (in Chinese)
[7]
李萍, 庄开莲, 李静. 国外直升机旋翼翼型研究综述[J]. 直升机技术, 2007(3): 103-109. DOI:10.3969/j.issn.1673-1220.2007.03.025
[8]
WILBY P G. The development of rotor airfoil testing in the UK[J]. Journal of the American Helicopter Society, 2001, 46(3): 210. DOI:10.4050/jahs.46.210
[9]
DUMLUPINAR E, MURTHY V. CFD analysis of unsteady aerodynamics of airfoils for rotor applications[C]//36th AIAA Fluid Dynamics Conference and Exhibit, San Francisco, California. Reston, Virginia: AIAA, 2006. doi: 10.2514/6.2006-3375
[10]
CARADONNA F. Developments and challenges in rotorcraft aerodynamics[C]//Proc of the 38th Aerospace Sciences Meeting and Exhibit, Reno, NV, USA. Reston, Virigina: AIAA, 2000. doi: 10.2514/6.2000-109
[11]
VIANA F A C, SIMPSON T W, BALABANOV V, et al. Special section on multidisciplinary design optimization: metamodeling in multidisciplinary design optimization: how far have we really come?[J]. AIAA Journal, 2014, 52(4): 670-690. DOI:10.2514/1.j052375
[12]
MARTIN J D, SIMPSON T W. Use of kriging models to approximate deterministic computer models[J]. AIAA Journal, 2005, 43(4): 853-863. DOI:10.2514/1.8650
[13]
THOMPSON J F. Grid generation techniques in computational fluid dynamics[J]. AIAA Journal, 1984, 22(11): 1505-1523. DOI:10.2514/3.8812
[14]
CARADONNA F X, TUNG C. Experimental and analytical studies of a model helicopter rotor in hover[R]. NASA-TM-81232, 1981. https://ntrs.nasa.gov/api/citations/19820004169/downloads/19820004169.pdf
[15]
CEZE M, HAYASHI M, VOLPE E. A study of the CST parameterization characteristics[C]//27th AIAA Applied Aerodynamics Conference, San Antonio, Texas. Reston, Virginia: AIAA, 2009. doi: 10.2514/6.2009-3767
[16]
LANE K, MARSHALL D. Inverse airfoil design utilizing CST parameterization[C]//48th AIAA Aerospace Sciences Meeting Including the New Horizons Forum and Aerospace Exposition, Orlando, Florida. Reston, Virigina: AIAA, 2010. doi: 10.2514/6.2010-1228
[17]
DADONE L U. A review of design objectives for advanced helicopter rotor airfoil[C]//National Symposium on Helicopter Aerodynamic Efficiency, 1975.
[18]
PERRY F J. Aerodynamic of the helicopter world speed record[C]//43rd Annual National Form of the American Helicopter Society, 1987.
[19]
FELKER F F, BETZINA M D, SIGNOR D B. Performance and loads data form a hover test of a full-scale XV-15 rotor[R]. NASA-TM-86833, 1985. https://ntrs.nasa.gov/api/citations/19860005773/downloads/19860005773.pdf
[20]
FELKER F F, SIGNOR D B, YOUNG L A, et al. Performance and load data from a hover test of a 0.658-scale XV-15 rotor and wing [R]. NASA-TM-89419, 1987. https://ntrs.nasa.gov/api/citations/19870015863/downloads/19870015863.pdf
[21]
MCVEIGH M A. The V-22 tilt-rotor large-scale rotor performance/wing download test and comparison with theory [C]//11th European Rotorcraft Forum, 1985: 18-25. https://dspace-erf.nlr.nl/xmlui/bitstream/handle/20.500.11881/3002/ERF%201985-97.pdf?sequence = 1&isAllowed = y
[22]
孙俊峰, 刘刚, 江雄, 等. 基于Kriging模型的旋翼翼型优化设计研究[J]. 空气动力学学报, 2013, 31(4): 437-441.
SUN J F, LIU G, JIANG X, et al. Research of rotor airfoil design optimization based on the Kriging model[J]. Acta Aerodynamica Sinica, 2013, 31(4): 437-441. (in Chinese)
[23]
肖中云, 江雄, 陈作斌, 等. 低速预处理方法在悬停旋翼模拟中的应用[J]. 航空学报, 2008, 29(2): 321-326.
XIAO Z Y, JIANG X, CHEN Z B, et al. Application of low speed preconditioning to hovering rotor simulation[J]. Acta Aeronautica et Astronautica Sinica, 2008, 29(2): 321-326. DOI:10.3321/j.issn:1000-6893.2008.02.011 (in Chinese)
[24]
袁红刚, 杨永东, 杨炯, 等. 典型直升机旋翼翼型气动特性试验研究[J]. 实验流体力学, 2013, 27(1): 20-24.
YUAN H G, YANG Y D, YANG J, et al. Experimental investigation on aerodynamic performance of one typical helicopter rotor airfoil[J]. Journal of Experiments in Fluid Mechanics, 2013, 27(1): 20-24. DOI:10.3969/j.issn.1672-9897.2013.01.004 (in Chinese)
[25]
王清, 招启军. 基于遗传算法的旋翼翼型综合气动优化设计[J]. 航空动力学报, 2016, 31(6): 1486-1495.
WANG Q, ZHAO Q J. Synthetical optimization design of rotor airfoil by genetic algorithm[J]. Journal of Aerospace Power, 2016, 31(6): 1486-1495. (in Chinese)
[26]
招启军, 王清, 赵国庆. 旋翼翼型定常-非定常特性综合优化设计新方法[J]. 南京航空航天大学学报, 2014, 46(3): 355-363.
ZHAO Q J, WANG Q, ZHAO G Q. New optimization design method for rotor airfoil considering steady-unsteady characteristics[J]. Journal of Nanjing University of Aeronautics & Astronautics, 2014, 46(3): 355-363. DOI:10.3969/j.issn.1005-2615.2014.03.004 (in Chinese)
[27]
WANG Q, ZHAO Q J. Rotor airfoil profile optimization for alleviating dynamic stall characteristics[J]. Aerospace Science and Technology, 2018, 72: 502-515. DOI:10.1016/j.ast.2017.11.033
[28]
夏露, 高正红, 苏伟. Pareto遗传算法在气动外形优化中的应用[J]. 空气动力学学报, 2007, 25(2): 194-198.
XIA L, GAO Z H, SU W. The application of Pareto GAs to aerodynamic configuration design[J]. Acta Aerodynamica Sinica, 2007, 25(2): 194-198. DOI:10.3969/j.issn.0258-1825.2007.02.010 (in Chinese)
[29]
卜月鹏, 宋文萍, 韩忠华, 等. 基于CST参数化方法的翼型气动优化设计[J]. 西北工业大学学报, 2013, 31(5): 829-836.
BU Y P, SONG W P, HAN Z H, et al. Aerodynamic optimization design of airfoil based on CST parameterization method[J]. Journal of Northwestern Polytechnical University, 2013, 31(5): 829-836. DOI:10.3969/j.issn.1000-2758.2013.05.028 (in Chinese)
[30]
许瑞飞, 宋文萍, 韩忠华. 改进Kriging模型在翼型气动优化设计中的应用研究[J]. 西北工业大学学报, 2010, 28(4): 503-510.
XU R F, SONG W P, HAN Z H. Application of improved kriging-model-based optimization method in airfoil aerodynamic design[J]. Journal of Northwestern Polytechnical University, 2010, 28(4): 503-510 (in Chinese).
[31]
HAN Z H, GÖRTZ S. Hierarchical kriging model for variable-fidelity surrogate modeling[J]. AIAA Journal, 2012, 50(9): 1885-1896. DOI:10.2514/1.J051354
[32]
许和勇, 叶正寅, 王刚, 等. 基于非结构嵌套网格的旋翼前飞流场计算[J]. 西北工业大学学报, 2006, 24(6): 763-767.
XU H Y, YE Z Y, WANG G, et al. Improving numerical simulation of rotor forward flight flow field with unstructured dynamic overset grids[J]. Journal of Northwestern Polytechnical University, 2006, 24(6): 763-767. DOI:10.3969/j.issn.1000-2758.2006.06.020 (in Chinese)
[33]
陈笑天, 吴裕平, 田旭. 旋翼翼型中高速综合气动优化设计方法研究[J]. 航空科学技术, 2019, 30(9): 19-24.
CHEN X T, WU Y P, TIAN X. Research on comprehensive aerodynamic optimum design method of rotor airfoil at medium and high speed[J]. Aeronautical Science & Technology, 2019, 30(9): 19-24. (in Chinese)
[34]
ZHANG W G, SUN J F, WANG L Q, et al. Rotor airfoil aerodynamic design method and wind tunnel test verification[J]. Chinese Journal of Aeronautics, 2020, 33(8): 2123-2132. DOI:10.1016/j.cja.2020.03.013
[35]
肖中云, 郭永恒, 张露, 等. 直升机CFD仿真现状与发展趋势分析[J]. 空气动力学学报, 2021, 39(4): 14-25.
XIAO Z Y, GUO Y H, ZHANG L, et al. An analysis of current status and prospects of CFD based simulation of rotorcrafts[J]. Acta Aerodynamica Sinica, 2021, 39(4): 14-25. DOI:10.7638/kqdlxxb-2020.0067 (in Chinese)
[36]
LANGTRY R B, MENTER F R. Correlation-based transition modeling for unstructured parallelized computational fluid dynamics codes[J]. AIAA Journal, 2009, 47(12): 2894-2906. DOI:10.2514/1.42362