一种无人机旋翼设计方法

下面是小编为大家整理的一种无人机旋翼设计方法,供大家参考。

　

　一种无人机旋翼设计方法

　摘要：为了提高多旋翼无人机的航时，本文提出了一种无人机旋翼设计方法。该方法根据给定的爬升速度、旋翼转速、拉力、旋翼直径、桨叶数、翼型，能够计算出最大力效旋翼的几何特性，包括弦长分布、桨叶角分布。利用该方法设计了某型无人机的旋翼，得到了桨叶角和弦长沿径向的分布。制作了旋翼模型，并且开展了旋翼台架试验。在不同转速下，测量出拉力、扭矩、功率、力效等参数，根据计算出的旋翼数据建模，采用多重参考系模型进行气动性能计算。经验证，试验结果与计算吻合，并且旋翼具有较高力效。

　关键词：无人机；旋翼；设计方法；拉力；力效 中 图 分 类 号 ：

　V211.44 文 献 标 识 码 ：

　ADOI ：10.19452/j.issn1007-5453.2022.04.001 基金项目：辽宁省教育厅重点攻关和服务地方项目（JYT19004，JYT19022，JYT2022161，JYT2022162）

　0 20 世纪 0 30 年代，英国研制出第一架可复用无人机 “ 蜂王 ” ，开启了无人机研制的先河。经过近 90 年的发展，无人机已经成为一类专门的飞行器。无人机具有成本低、机动性好、使用方便等优点。在军事领域，无人机成为新型的武器平台，可以执行武装巡逻、电子干扰、侦察监视、对地攻击、通信中继、目标定位和攻击评估等任务。在民用领域，无人机具有航空拍摄、农业植保、航空物流、交通监控、地貌测量和

　救灾防火等用途。旋翼无人机采用旋翼产生拉力，为了提高无人机的航时，必须设计出高效率的旋翼。

　旋翼桨叶气动设计采用的基本理论包括动量理论、叶素理论和涡流理论。国内外很多学者开展了旋翼的气动设计与性能分析研究。J.M.Grasmeyer 等[1]以诱导损失最小为目标设计了“黑寡妇”微型飞行器的旋翼，并且测试了旋翼的拉力和效率。H.I.Kwon 等[2]采用可变信度气动分析和多级优化方法设计了某型旋翼，并开展了计算流体力学（CFD）分析和风洞试验。P.Hajela 等[3]采用遗传算法对旋翼桨叶进行多学科优化设计。J.Ahmad 等[4]基于移动重叠网格进行了非定常旋翼叶片的气动计算。王立群等[5]采用有限体积法在桨叶固接坐标系直接求解三维欧拉方程，不附加任何尾迹模型，计算了桨叶表面压力系数和沿展向升力系数，将计算结果与试验数据进行了对比。杨爱明等[6]采用运动嵌套网格方法模拟了旋翼前飞非定常流场，计算了旋翼悬停流场，旋翼桨叶表面压力分布的计算值与试验值吻合很好。招启军等[7]充分考虑旋翼尾迹对流场的影响和减少尾迹的数值耗散，建立了一个基于 Navier-Stokes 方程/自由尾迹分析/全位势方程的旋翼流场求解的混合方法，以两叶的 Caradorma&Tung模型旋翼和四叶的 UH-60A 直升机旋翼为算例，计算给出了旋翼桨叶表面的压强分布以及桨尖涡的位置。陈平剑等[8]研究了桨叶气动外形设计方法，设计参数主要包括桨叶的弦

　长、扭转角、翼型选择和配置以及桨尖形状，以此方法设计了某型直升机旋翼桨叶。谢辉等[9]针对低速中小型无人机，设计了一种矩形薄型直桨叶两叶螺旋桨，该螺旋桨在静拉力、功率、效率、拉力系数、功率系数等关键性指标上均表现出优异特性。臧士新[10]搭建了旋翼桨叶快速设计平台，建立了气动设计模块，辅以结构设计模块，编制了旋翼桨叶设计软件，并用算例验证了平台设计功能的有效性。招启军等[11]基于运动嵌套网格技术建立了适用于共轴刚性双旋翼悬停流场模拟的 CFD 方法，研究了桨叶平面外形参数对共轴刚性旋翼悬停性能的影响。叶靓等[12]将非结构嵌套网格方法和网格自适应技术相结合，发展了一套适合于共轴式双旋翼流场数值计算的求解器，应用该求解器计算了共轴式双旋翼的桨尖涡轨迹和拉力分布特性。唐正飞等[13]全面比较了共轴式双旋翼与单旋翼尾迹流场中轴向、径向和周向速度的特点与差异，研究结果表明，双旋翼流场与单旋翼明显不同，而这些不同沿各向异性，轴向最大，径向次之，周向最小。童自力等[14]通过数值求解带有动量源项的三维不可压缩 N-S方程，模拟了共轴式双旋翼的流动，探讨了两旋翼间相互干扰的特性。童自力等[15]将 N-S 方程中加入动量源项进行共轴式双旋翼的气动力计算，该方法能够较好地预测双旋翼的气动特性。许和勇等[16]对悬停共轴双旋翼复杂流场进行了数值模拟研究，研究结果表明，双旋翼之间的干扰使得两者

　的拉力性能下降，上旋翼对下旋翼的干扰程度远大于下旋翼对上旋翼的影响，双旋翼之间距离增加，则上旋翼拉力系数增加，下旋翼拉力系数减小，总拉力系数减小。张睿等[17]通过不同的网格模型用以优化机翼外形，最终得出了理想的机翼。李春华等[18]阐述了共軸刚性旋翼的概念以及主要需要的技术，并且对技术的发展提出了建议。

　本文提出了一种无人机旋翼设计方法，该方法使用弦长分布函数拟合旋翼弦长，利用高阶面元法计算截面的最大升阻比迎角，最终获得最高力效旋翼的弦长分布和桨叶角分布。加工了两叶旋翼，并且进行了性能测试，将性能试验结果与性能计算结果进行了对比。由该方法设计出的旋翼在能满足要求的同时，也能具有较高的力效，应用价值较高。

　1 旋翼设计方法 在 旋翼拉力、无人机爬升速度、旋翼转速给定的情况下，最高力效旋翼的能量损失是最小的。当桨叶所有截面都处于最大升阻比迎角时，旋翼的力效是最高的。旋翼桨叶的某一截面如图 1 所示。

　2 某型无人机旋翼的设计 径 某型无人机采用两叶旋翼，旋翼半径 R=0.279m ，桨毂半径Rh=0.045m ，采用 6 RAF6 翼型，将桨叶划分成 6 6 段（7 7 个截面），无人机爬升速度 1m/s ，旋翼最大转速 5000r/min ，需用拉力T=4kgf （ 1kgf≈0.98N ）。根据上节所述的旋翼设计方法，对

　某型两叶旋翼进行了设计，得到的桨叶角分布如图 2 所示，弦长分布如图 3 所示。

　從图 2 2 可以看出，桨根桨叶角最大，桨尖桨叶角最小，从桨， 根到桨尖，桨叶角逐渐减小， 75% 半径处的桨叶角为 6.46° 。

　从图 3 3 可以看出，从桨根到桨尖，弦长先增大然后减小，最大弦长位于第二个截面，为 0.053m ；桨尖弦长最小，为 0.011m 。根据桨叶角分布及弦长分布，建立了旋翼三维数模，如图 4所示。根据旋翼三维数模，在辽宁飓风科技有限公司加工了两叶旋翼，材料为德国榉木，旋翼如图 5 所示。

　3 旋翼气动性能计算与试验 本文采用 F MRF 方法对旋翼气动性能进行计算，其主要思想是在空间建立两个流场域：一个是旋转域， 另一个是静止域。旋转域用来模拟旋转运动，静止域用来模拟空间流场，旋转域和静止域通过交界面建立流场间的联系，旋转域网格如图6 所示。

　计算采用 k k- -T ωSST 湍流模型，k k- -T ωSST 湍流模型具有较好的普适性，在工程上得到广泛应用。网格数量为 593 万。边界条件设置为：速度入口、压力出口。将入口定义为速度入口，添加来流速度变量，出口定义为环境压力 101325Pa，流场域如图 7 所示。

　当各项计算数值均收敛后，进行数据的后处理。取转速为1895r/min 时的流场特性进行分析，旋翼的压力云图以及流

　线图如图 8、图 9 所示。由图 8 可知，旋翼压力最小值出现在旋翼的尖端，从尖端到根部压力逐渐增加，在 r=0.65R～R上压力变化趋势较大，而后趋于平缓。

　由图 9 9 可知，靠近旋翼的部分流速较大，随着与旋翼的距离逐渐增大，流速逐渐变小，且整个区域的流线没有出现分离现象，证明此旋翼具有较高的力效。

　试验台测试旋翼的最大直径为 106.7cm ，拉力传感器量程为0 0 ～ 12kgf ，转速传感器量程为 1500 ～ 300000r/min ，扭矩传感器量程为 0 0 ～ 8N·m 。本文将拉力、扭矩以及旋翼吸收功率的试验结果与 CFD 计算结果进行比较，以此验证计算的精度，对比结果如图 10～图 12 所示。

　图 由图 0 10 可知，随着转速增加，旋翼拉力的试验值与计算值差距逐渐增加，差距维持在 10% 以内。

　图 由图 1 11 可知，旋翼扭矩的试验值与计算值始终相差不超过10% 。由此可见，计算结果与实际结果相差不大，较为吻合。由图 12 可知，旋翼吸收功率试验值与计算值十分接近，CFD计算得出的结果与实际相符，精度较高。

　4 结论 本文提出了一种无人机旋翼设计方法，该方法计算出桨叶所有截面最大升阻比迎角，使得旋翼具有较高力效。利用该方法设计了某型无人机的旋翼，采用MRF方法对旋翼进行了CFD计算，并进行了试验台测试对比验证，得出以下结论：（1）

　本文方法设计的旋翼具有高的力效，能够满足某型无人机的拉力和力效需求；（2）CFD 计算结果与试验结果较为吻合，证明 MRF 方法可以为此类工程问题提供较为准确的预估与计算；（3）本文提出的无人机旋翼设计方法具有较高的工程应用价值。

　参考文献 [1]GrasmeyerJM ，KeennonMT.Developmentoftheblackwidowmicroairvehicle[J].ProgressinAstronauticsandAeronautics，2022，195：519-535. [2]KwonHI ， YiS ， ChoiS ，etal.Designofefficientpropellersusingvariable-fidelityaerodynamicanalysisandmultileveloptimization[J].JournalofPropulsion&Power，2022，31（4）：1-16. [3]HajelaP ，LeeJ.Geneticalgorithmsinmultidisciplinaryrotorbladedesign[C]//StructuralDynamicsandMaterialsConference ，1995. [4]AhmadJ ，DuqueEP.Helicopterrotorbladecomputationinunsteadyflowsusingmovingoversetgrids[J].AIAAJournalofAircraft ，1996，33（1）：54-60.

　[5]王立群，乔志德，张茹.旋翼定常流场的欧拉方程计算[J].计算物理，1997，14（6）：841-845.WangLiqun，QiaoZhide，ZhangRu.TheEulerequationscalculationofthesteadyflowfieldoftherotor[J].ChineseJournalofComputationalPhysics，1997，14（6）：841-845.（inChinses）

　[6]杨爱明，乔志德.基于运动嵌套网格的前飞旋翼绕流 N-S方 程 数 值 计 算 [J]. 航 空 学 报 ， 2000 ， 22 （ 5 ）：434-436.YangAiming ，QiaoZhide.N-Sequationsofflowaroundaforward-flyingrotorbasedonmovingnestedgrids[J].ActaAeronauticaetAstronauticaSinica，2000，22（5）：434-436.（inChinese）

　[7]招启军，徐国华.基于 Navier-Stokes 方程/自由尾迹/全位势方程的旋翼流场模拟混合方法[J].空气动力学学报，2022 ， 27 （ 1 ）

　：

　15-21.ZhaoQijun ，XuGuohua.HybridmethodforrotorflowfieldsimulationbasedonNavier-Stokesequation/freewake/fullpotentialequation[J].ActaAerodynamicaSinica，2022，27（1）：15-21.（inChinese）

　[8]陈平剑，徐玉貌.旋翼桨叶气动外形设计[J].直升机技术，2022 （ 3 ）

　：

　16-19.ChenPingjian ，XuYumao.Aerodynamicconfigurationdesignofrotorblade[J].HelicopterTechnique，2022（3）：16-19.（inChinese）

　[9]谢辉，王力，张琳.一种适用于中小型无人机的新型螺旋桨设计[J].航空工程进展，2022，6（1）：71-76.XieHui，WangLi ，ZhangLin.Anewpropellerdesignforsmallandmedium-sizedUAV[J].AdvancesinAeronauticalScienceandEngineering ，2022，6（1）：71-76.（inChinese）

　[10]臧士新.直升機旋翼桨叶气动外形快速设计技术研究[D]. 南 京 ：

　南 京 航 空 航 天 大 学 ，2022.ZangShixin.Researchonrapidaerodynamicdesigntechnologyofhelicopterrotorblade[D].Nanjing ：NanjingUniversityofAeronautics&Astronautics ， 2022.（inChinese）

　[11]招启军，朱正，原昕.桨叶外形对共轴刚性旋翼悬停性能影响的 CFD 分析[J].南京航空航天大学学报，2022，79 （5）：653-661.ZhaoQijun ， ZhuZheng ，YuanXin.CFDanalysisoftheinfluenceofbladeshapeonthehoverperformanceofcoaxialrigidrotors[J].JournalofNanjingUniversityofAeronautics&Astronautics，2022，79（5）：653-661.（inChinese）

　[12]叶靓，徐国华.共轴式双旋翼悬停流场和气动力的 CFD计算[J].空气动力学学报，2022（4）：437-442.YeLiang，XuGuohua.CFDcalculationofflowfieldandaerodynamicfor

　ceofcoaxialdoublerotorsinHover[J].ActaAerodynamicaSinica，2022（4）：437-442.（inChinese）

　[13]唐正飞，高正.共轴双旋翼与单旋翼悬停流场试验测量值的对比[J].南京航空航天大学学报，1997，29（6）：627-632.TangZhengfei ，GaoZheng.Comparisonofmeasuredvaluesbetweencoaxialrotorandsinglerotorinhoverflowfield[J].JournalofNanjingUniversityofAeronautics&Astronautics，1997，29（6）：627-632.（inChinese）

　[14]童自力，孙茂.共轴式双旋翼流动的 N-S 方程模拟[J].航 空 学 报 ， 1998 ， 19 （ 1 ）：

　1-5.TongZili ，SunMao.N-Sequationssimulationofcoaxialrotorflow[J].ActaAeronauticaetAstronauticaSinica，1998，19（1）：1-5.（inChinese）

　[15]童自力，孙茂.共轴式双旋翼气动力特性的计算研究[J].航 空 学 报 ， 1999 ， 20 （ 4 ）：

　348-350.TongZili ，SunMao.Computationalstudyontheaerodynamiccharacteristicsofcoaxialdoublerotors[J].ActaAeronauticaetAstronauticaSinica，1999，20（4）：348-350.（inChinese）

　[16]许和勇，叶正寅.悬停共轴双旋翼干扰流动数值模拟[J].航 空 动力 学 报， 2022， 26（2 ）：453-457.XuHeyong ，YeZhengyin.Numericalsimulationofcoaxialrotorinterac

　tionflowinhover[J].JournalofAerospacePower，2022，26（2）：453-457.（inChinese）

　[17]张睿，崔德刚，裴志，等.基于不同网格模型的机翼型架外形分布优化设计[J].航空科学技术，2022，32（12）：18-26.ZhangRui ， CuiDegang ， PeiZhi ，etal.Shapedistributionoptimizationdesignofairfoilframebasedondifferentmeshmodels[J].AeronauticalScience&Technology，2022，32（12）：18-26.（inChinese）

　[18]李春華，樊枫，徐明.共轴刚性旋翼构型高速直升机发展研究[J].航空科学技术，2022，32（1）：47-52.LiChunhua，FanFeng ，XuMing.Developmentofhigh-speedhelicopterwithcoaxialrigidrotorconfiguration[J].AeronauticalScience&T...