刘 毅,赵晓霞,欧阳绍修
(中航飞机股份有限公司研发中心,陕西汉中 723000)
螺旋桨飞机升力失速特性研究
刘 毅*,赵晓霞,欧阳绍修
(中航飞机股份有限公司研发中心,陕西汉中 723000)
飞机的最大升力系数CLmax直接影响翼载的选取,进而影响飞机的重量和经济性。螺旋桨飞机在带动力条件下CLmax随滑流强度增加而提高,但按照常规理论采用无动力CLmax数据选取的翼载偏小偏保守,没有充分发掘飞机的性能潜力。结合适航规定以及某四发螺旋桨飞机飞行实际情况,提出了一种基于发动机慢车状态确定CLmax的概念,并将飞机带动力CLmax分解为无动力CLmax、螺旋桨拉力、螺旋桨法向力、滑流增升效应等4部分贡献,通过无动力和多天平带动力风洞实验完成了上述分量的模拟、测量和修正。计算表明某四发螺旋桨飞机在发动机慢车失速试飞条件下的滑流强度约为0.1,螺旋桨系统的动力增升作用使不同襟翼构型的CLmax增加8% ~9%。该方法获得的CLmax与飞机试飞结果较为吻合,充分挖掘了飞机的低速性能潜力,并为同类螺旋桨飞机设计提供了一定的参考。
:螺旋桨飞机;滑流;最大升力系数;法向力;风洞实验
螺旋桨发动机是载人航空初期的主要动力系统,直至二战结束,主要的固定翼飞机均为螺旋桨飞机。其后由于飞行速度的提高,军/民用航空进入了涡喷和涡扇动力时代,但在飞行M数低于0.6的低速范围内,螺旋桨动力由于推力大、油耗低等优点仍得到了广泛运用,其范例包括经典的C-130系列、国产运7、运8、运12以及最新的A-400M、V-22倾转旋翼机等。
螺旋桨气动特性的理论研究起源于19世纪,经历了动量理论、叶素理论、涡流理论和计算流体力学(CFD)等主要阶段[1-2],螺旋桨与机翼之间的相互干扰研究历程与之相似。单独螺旋桨流场的主要特点是使通过桨盘的气流加速、旋转,当螺旋桨前进速度为0时滑流直径将收缩为桨盘直径的0.816至0.92倍,而安装到飞机后在正常飞行速度时滑流直径与桨盘直径相当,在几个桨盘直径后方滑流大致以15度的角度扩张。螺旋桨安装至机翼后二者之间存在相互干扰作用,一方面螺旋桨滑流使机翼局部速压增加并附加有旋转速度分量,另一方面机翼的上洗效应使螺旋桨的有效迎角增加,法向力随之增加,总的结果是二者组合的升力特性大于单独部件升力之和[3]。线性段的带滑流效应升力特性可由动量理论进行估算,文献[4]指出有滑流状态下的最大升力系数需要考虑滑流导致的失速迎角推迟效应,还将进一步增加。
螺旋桨与飞机的相互干扰可以通过计算和实验两类手段来进行研究。滑流影响的计算分析方法包括涡格法[5-6]和求解Euler/N-S方程方法两大类,其中后者对螺旋桨效应的模拟主要有激励盘(等效盘)[7-10]和多参考系滑移网格[11-14]两种,上述研究结果表明计算方法可以对滑流影响的趋势和流场基本特征提供有意义的参考,但Veldhuis的对比表明CFD获得的气动力与风洞实验数据的差异相对涡格法并无缩小[2]。带动力风洞实验则具有相对较高的精度以及对真实流态的模拟能力,赵学训、李征初和Catalano分别采用测压和流场测量等方法对飞机滑流影响区的压力分布特性、空间流场和附面层特性进行了研究[15-17];Russell、Gentry、Petrov和唐克兵对带动力的气动力和力矩特性进行了实验研究[18-21];李尚斌研究了螺旋桨与飞机分别独立支撑的滑流影响实验,获得了纯滑流的气动力影响量[22];欧阳绍修研究了单体及安装于飞机上的螺旋桨法向力特性[23]。
上述动力影响风洞实验对于螺旋桨飞机的设计实践和具体应用还存在一定盲区:第一是多数没有对螺旋桨直接力(如拉力、扭矩、法向力等)和滑流影响进行独立测量,其中滑流影响量一般按照间接模拟法原理模拟滑流效应,但直接力和力矩与真实飞机一般不相似,导致对飞机带动力气动特性评估的误差;第二是滑流带来的CLmax收益与适航规范和实际飞行状态脱节,难以直接应用于飞机设计,例如AN-10飞机带动力CLmax可达5.4,AN-70的CLmax甚至达到7[24],但各种军民用飞机规范对失速条件下的动力系统状态均有严格规定,失速试飞的动力系统状态与实验条件并不对应。
本文采用了一种多天平测力和法向力模拟风洞实验方法,实现了螺旋桨滑流影响与直接力影响的独立测量,经修正后二者均满足相似条件。同时根据民机适航规范要求,确定了失速试飞条件下的滑流强度进而获得了可应用于飞机设计使用的带动力CLmax数据,从而挖掘了螺旋桨飞机的低速性能潜力。
1.1 失速点的动力状态
螺旋桨滑流效应的强弱可以由滑流强度B来衡量,其计算公式为:
其中,T为拉力,D为螺旋桨桨盘直径,q为飞行速压。
螺旋桨飞机失速试飞时往往要求发动机在小功率状态,但由于此时飞行速度也较低,拉力相对较大,最终仍可能产生一定的拉力和滑流强度。以下结合民机适航规定对失速试飞时的滑流强度进行分析。
CCAR-25-R4的25.103条规定失速速度试飞条件是“发动机慢车,或者如果产生的推力导致失速速度明显下降,在此失速速度时不超过零推力”,25.201条规定失速演示应在“无动力”状态进行。CCAR-23的规定与CCAR-25类似[25-26]。规范中“慢车”、“零推力”和“无动力”的规定意图是避免过大的动力增升效应带来失速速度显著降低的冒进结果,保证飞行安全。应用于某四发涡桨飞机时从安全角度考虑不可能进行无动力失速试飞,在设计的各个飞行阶段临界发动机停车后至少须保证2至3台发动机工作,因此应当采用发动机最小功率稳定工作的慢车条件作为失速特性的考核点。根据某型机发动机数据计算得到失速试飞状态的滑流强度见表1,可见发动机慢车时仍产生了0.08~0.1的滑流强度。
表1 失速试飞动力状态的滑流强度Table 1 Slipstream intensity at stall flight test
1.2 螺旋桨飞机的CLmax分解
根据前文分析,螺旋桨飞机在失速试飞时是保持(多发)螺旋桨以较小功率工作的,因此实际升力应包含4部分贡献,即全机无动力升力、滑流增量(动力间接影响)、螺旋桨法向力和拉力(动力直接影响),示意图见图1。图中无动力升力定义为不考虑动力影响的升力,滑流产生的升力是指螺旋桨后气流动能增加带来的飞机升力增加量,二者均为风轴系;螺旋桨拉力定义为体轴中向前的驱动力,法向力与拉力垂直,拉力和法向力均有升力方向的分量。上述4部分作用在升力方向上的和即为飞机的带动力升力,因此CLmax也应考虑这些影响。
图1 螺旋桨飞机升力的构成Fig.1 Breakdown of the lift of propeller driven aircraft
全机无动力升力可由实验或计算方法获得,螺旋桨拉力在升力方向的分量可由发动机数据获取,量值一般较小。螺旋桨滑流增量与法向力可通过带动力风洞实验来获取,其主要难点是滑流影响量通常采用间接模拟法来模拟滑流效果,此时螺旋桨的法向力与真实飞机一般是不相似的,若不将滑流与直接力单独测量则会在结果中引入系统误差。
1.3 滑流与法向力的分离测量
飞机滑流影响风洞实验采用了固定拉力系数的间接模拟法,通过选择4-5个涵盖实际飞机拉力系数范围的状态点进行动力影响实验,模拟相似参数兼顾了拉力系数、扭矩系数和前进比与实际状态的一致性以满足螺旋桨绕流特性的相似。为解决滑流与直接力单独测量的难题,采用了多天平带动力实验技术,即每套电机与螺旋桨组成的动力系统与飞机模型不直接接触,该系统的力和力矩通过电机天平再传递至飞机模型,最后在全机气动力读数中扣除所有电机天平的读数,可以得到“纯滑流”影响下的飞机气动特性,电机模型示意图见图2。螺旋桨直接力主要包括拉力和法向力两部份,其中拉力由发动机厂家直接提供,法向力则采用直接模拟桨叶角、前进比并考虑与飞机之间相互干扰影响的实验方法获得,其详细原理可参见相关文献[23]。
图2 风洞实验的电机螺旋桨Fig.2 Electrical propeller used in the wind tunnel test
本文的带动力风洞实验方法在中国空气动力研究与发展中心低速所的8 m×6 m风洞完成了验证,实验模型见图3。8 m×6 m风洞是一座直流式、闭口、串列双实验段的大型风洞,本文使用的实验段宽8 m、高6 m、长15 m、有效截面积47.4 m2,常用风速20 m/s至85 m/s。螺旋桨动力系统主要由60 kW无刷电机、无刷直流电动机调速装置、操作控制台和冷却水系统组成,电机额定功率为60 kW。
从图片可见带动力飞机实验模型将电机内置于发房内,相应的驱动、散热、控制等管线全部布置在机翼内部,最大限度地避免了对气动外形的破坏,保证了实验效果。
图3 带动力风洞实验模型Fig.3 Model of the powered wind tunnel test
气动力测试采用了5台盒式六分量应变天平分别测量全机和四台电机的气动力,主天平、小天平法向力设计载荷分别为25 000 N、1 600 N,静态校准天平性能参数的测量不确定度不超过0.05%,准确度误差不超过0.2%;轴向力设计载荷分别为8 000 N、1500 N,测量不确定度不超过0.08%,准确度误差不超过0.2%,满足GJB 2244A-2011指标和动力影响实验的精度要求。风洞实验的Re数为1.6×106,动力系统实现的滑流强度最大为2.2。
某型机巡航、起飞、着陆构型在纯滑流影响下的升力曲线见图4~图6。从升力曲线可见滑流使飞机升力线斜率、CLmax均增加,增量正比于滑流强度。巡航构型时不同滑流强度的零升迎角基本相等,可见对单段翼型剖面的机翼而言滑流增升的机理主要是增加了滑流影响区域内的速压,区域内机翼翼段的净升力增加,除以无动力参考面积后得到了增加的升力系数。某型机由于采用了4台涡桨发动机,滑流吹洗面积达到毛机翼面积的50%,因此在较小的滑流强度下也有可观的升力增量。滑流对非线性段升力和CLmax的贡献进一步增加,失速迎角推迟约1°~2°,此现象主要与较大迎角下机翼上翼面受到加速气流吹洗,附面层内气流能量增加,气流分离推迟有关。
放襟翼条件下滑流的增升效应更为显著,其主要原因是襟翼几乎都位于滑流影响区,对滑流的响应也更为敏感。在起飞大拉力条件下滑流增升带来的CLmax增量可达到与无动力CLmax相当的量级,从而具有大幅降低离地速度的潜力。着陆构型由于襟翼偏度和后退量最大,在同等滑流强度下增升最为显著,在前文计算得到的发动机慢车时B=0.1的条件下CLmax也有明显的增加。
图4 巡航构型的带动力升力特性(扣除螺旋桨直接力)Fig.4 Lift curves of cruise configuration with powered propeller(direct force of propeller subtracted)
图5 起飞构型的带动力升力特性(扣除螺旋桨直接力)Fig.5 Lift curves of takeoff configuration with powered propeller(direct force of propeller subtracted)
图6 着陆构型的带动力升力特性(扣除螺旋桨直接力)Fig.6 Lift curves of landing configuration with powered propeller(direct force of propeller subtracted)
不同襟翼下滑流效应产生的CLmax增量见图7,可见ΔCLmax~B变化曲线呈单调递增的趋势,随着B的增加增速变缓。在图中标示出了起飞、巡航和着陆发动机慢车失速试飞点,可见失速试飞时滑流带来的ΔCLmax远小于动力增升的能力边界,但仍具有较可观的量值,在飞机设计中考虑这一因素可以更加准确地评估飞机的失速特性。
图7 滑流产生的CLmax增量Fig.7 CLmaxincrement due to slipstream
图8 不同构型的CLmax分解Fig.8 CLmaxbreakdown of different configuration
考虑全机无动力、滑流、拉力和法向力全部作用后的CLmax总量见图8。拉力分量的贡献是直接将发动机慢车状态的拉力无量纲化后得到的。法向力根据文献[23]的原理和方法按照相似准则测量修正得到。无动力数据根据风洞实验结果修正后得到。从各组成部分来看飞机无动力状态固有的CLmax仍是带动力CLmax的主要贡献因素,各构型条件下拉力和法向力贡献量基本相当,滑流贡献量随襟翼偏度增加而增加。总的动力影响使各襟翼CLmax均有8% ~9%的增长,考虑此影响后在飞机设计初期选择翼载时可提高相应比例,从而起到降低机翼面积,降低飞机重量并提升效能的作用。
本文根据适航规定和某型机飞行时动力系统的实际工作状态,在发动机慢车状态获取了更加符合真实情况的飞机CLmax。通过将动力影响量分解为滑流、法向力和拉力分量并采用改进的风洞实验方法进行相似模拟,定量测定表明滑流可使各襟翼失速点的CLmax增加约8% ~9%。通过对比某4发涡桨飞机试飞反推的CLmax与风洞实验修正值,表明通过合理的相关性修正,二者的差异可控制在±0.04以内,相对比例不超过±2%,吻合良好。本文的研究方法有效挖掘了螺旋桨飞机的低速性能潜力,提升了飞机效能,可以为螺旋桨飞机总体气动布局设计提供重要参考。
[1]Li Shangbin,Jiao Yuqin.The investigate development of propeller slipstream’s effect[J].Science Technology and Engineering,2012,12(8):1867-1873.(in Chinese)李尚斌,焦予秦.螺旋桨滑流影响的研究进展[J].科学技术与工程,2012,12(8):1867-1873.
[2]Veldhuis L L M.Review of propeller-wing aerodynamic interference[C].24th international congress of the aeronautical science,2004.
[3]Hoerner S F,Borst H V.Fluid-dynamic lift 2ed[M].Published by Mrs.Liselotte A.Hoerner:1985,(12-1)-(12-24).
[4]Committee of Aircraft Design Manual.Aircraft design manual(Volume 6)aerodynamic design[M].Beijing:National Defense Industry Press,2002:522-531.(in Chinese)《飞机设计手册》总编委会.飞机设计手册(第6册)气动设计[M].北京:国防工业出版社,2002:522-531.
[5]Bronswijk N.The effects of propeller power on the stability and control of a tractor-propeller powered single-engine low-wing monoplane[D].University of Sydney,2001.
[6]Wang Chuanbin.Investigation of the engineering calculation method of the aerodynamic interference between propeller slipstream and aircraft components[D].Nanjing:Nanjing University of Aeronautics and Astronautics,2006.(in Chinese)王传斌.螺旋桨滑流与飞机部件气动干扰的工程估算方法研究[D].南京:南京航空航天大学,2006.
[7]Duan Yiqian,Shi Ai’ming.A new and effective actuator disk model approach for the simulation of propeller slipstream[J].Journal of Northwestern Polytechnical University,2012,30(6):841-846.(in Chinese)段义乾,史爱明.一种新型的螺旋桨滑流激励盘模型的研究方法[J].西北工业大学学报,2012,30(6):841-846.
[8]Li Bo,Liang Dewang,Huang Guoping.Propeller slipstream effects on aerodynamic performance of turbo-prop airplane based on equivalent actuator disk model[J].Acta Aeronautica et Astronautica Sinica,2008,29(4):845-852.(in Chinese)李博,梁德旺,黄国平.基于等效盘模型的滑流对涡桨飞机气动性能的影响[J].航空学报,2008,29(4):845-852.
[9]He Hui.Investigation of the extreme low-altitude flow field numerical simulation method of propeller aircraft[D].Xi'an:Northwestern Polytechnical University,2007.(in Chinese)何辉.螺旋桨飞机超低空流场数值模拟方法研究[D].西安:西北工业大学,2007.
[10]Zuo Suihan.Investigation of the interference of propeller slipstream to wing/aircraft using numerical method[D].Northwestern Polytechnical University,2007.(in Chinese)左岁寒.螺旋桨滑流对机翼/全机干扰的数值模拟研究[D].西北工业大学,2007.
[11]Ren Xiaofeng,Yang Shipu,Duan Zhuoyi,et al.Investigation of the propeller slipstream effects on wing aerodynamic characteristics based on multiple frame of reference[C]//Proceedings of 14thnational computational fluid dynamics conference,2008:582-585.(in Chinese)任晓峰,杨士普,段卓毅,等.基于多参考坐标系的螺旋桨滑流对机翼气动特性影响分析[C]//第14届全国计算流体力学会议论文集,2008:582-585.
[12]Cheng Xiaoliang,Li Jie.Unsteady computational method for the propeller/wing interaction[J].Science Technology and Engineering,2011,11(14):3229-3235.程晓亮,李杰.螺旋桨滑流对机翼气动特性影响的方法研究[J].科学技术与工程,2011,11(14):3229-3235.
[13]Zhang Liu,Bai Junqiang,Li Huaxing,et al.Research on aerodynamic interference for propeller slipstream over the wing[J].Aeronautical Computing Technique,2012,42(2):87-92.(in Chinese)张刘,白俊强,李华星,等.螺旋桨滑流与机翼之间气动干扰影响研究[J].航空计算技术,2012,42(2):87-92.
[14]Zhang Xin,Yang Yong.Application of dynamic patched-grid on unsteady simulation of propeller flows[J].Science Technology and Engineering,2012,12(7):1564-1567.(in Chinese)张鑫,杨永.应用动态面搭接网格数值模拟双发螺旋桨流场[J].科学技术与工程,2012,12(7):1564-1567.
[15]Zhao Xuexun.Experiment research of the airflow which surrounds aircraft under the influence of propeller slipstream[J].Aerodynamic Experiment and Measurement&Control,1995,9(4):48-52.(in Chinese)赵学训.螺旋桨滑流对飞机绕流影响的实验研究[J].气动实验与测量控制,1995,9(4):48-52.
[16]Li Zhengchu,Wang Xunnian,Chen Hong,et al.Experimental research of influence of propeller slipstream on wing flow field[J].Measurements in Fluid Mechanics.2000,14:44-48.(in Chinese)李征初,王勋年,陈洪,等.螺旋桨滑流对飞机机翼流场影响实验研究[J].流体力学实验与测量,2000,14:44-48.
[17]Catalano F M.On the effects of an installed propeller slipstream on wing aerodynamic characteristics[J].Acta Polytechnica.2004,44 (3):8-14.
[18]Russell J S,McCoy H M,Propeller characteristics and slipstream effects on a high wing monoplane from wind tunnel tests[D].California Institute of Technology,1935.
[19]Gentry G L,Takallu M A,Applin Z T,Aerodynamic characteristics of a propeller-powered high-lift semispan wing[R].NASA technical memorandum 4541.1994.
[20]Petrov A V,Stepanov Y G,Shmakov M V.Development of a technique and method of testing aircraft models with turboprop enginesimulators in a small scale wind tunnel-results of tests[J].Acta Polytechnica.2004,44(2):27-31.
[21]Tang Kebin,Xia Shenglin.Effect investigation of slipstream to propeller-driven aircraft’s aerodynamic characteristics[C].Proceedings of the 1stcontemporary experimental aerodynamic conference 2007:405-408.(in Chinese)唐克兵,夏生林.螺旋桨滑流对飞机气动特性影响研究[C].第一届近代实验空气动力学会议文集,2007:405-408.
[22]Li Shangbin,Jiao Yuqin.The experiment investigation of the effect of propeller slipstream based on wing-tip support and independent propeller support[J].Engineering Mechanics,2013,30(7):288-293.(in Chinese)李尚斌,焦予秦.基于翼尖支撑和螺旋桨独立支撑的螺旋桨滑流影响实验研究[J].工程力学.2013,30(7):88-293.
[23]Ouyang Shaoxiu,Zhao Xiaoxia,Jiang Zonghui,et al.The research on wind tunnel test method of propeller’s radial force[J].Journal of Experiments in Fluid Mechanics.2012,26(3):91-94.(in Chinese)欧阳绍修,赵晓霞,江宗辉,等.螺旋桨径向力风洞实验方法研究[J].实验流体力学.2012,26(3):91-94.
[24]Petrov A V.Aerodynamics of STOL airplanes with powered high-lift systems[C].28th International congress of the aeronautical science.2012.
[25]CAAC.CCAR-25-R4 China civil aviation regulations part 25 Airworthiness standards:transport category airplanes[S].2011.(in Chinese)中国民用航空局.CCAR-25-R4中国民用航空规章第25部:运输类飞机适航标准[S].2011.
[26]CAAC.CCAR-23-R3 China civil aviation regulations part 23 Airworthiness standards:normal,utility,acrobatic,and commuter category airplanes[S].2004.(in Chinese)中国民用航空局.CCAR-23-R3中国民用航空规章第23部:正常类、实用类、特技类和通勤类飞机适航规定[S].2004.
Investigation on lift stall characteristics of propeller aircraft
Liu Yi*,Zhao Xiaoxia,Ouyang Shaoxiu
(Research and Development Center,AVIC Aircraft Co.,LTD,Hanzhong 723000,China)
The maximum lift coefficient CLmaxis directly connected to the selection of wing loading,which further affects aircraft’s mass and economy.CLmaxof propeller aircraft is improved with increased slipstream intensity at powered condition,however the conventional theory trends to some conservative and smaller power-off value,hence does not utilize the full potential performance of aircraft.Combining the airworthiness regulations and the actual flight condition of an aircraft with four turbo-prop engines,a new concept is proposed to find the optimal CLmaxbased on engine standby condition.The power-on CLmaxis separated into four parts:power-off CLmax,propeller thrust component,propeller normal force and lift due to slipstream,which are simulated,measured and corrected by power-off and multi-balance power-on wind tunnel tests.The calculation shows that the propeller aircraft with four engines has slipstream intensity in the order of 0.1 at stall flight test condition with engines standby.The slipstream effect is the major factor of lift increment even at small slipstream intensity,and the increment increases as flap angle or slipstream intensity increase.The four engines at standby condition also generate considerable thrust and normal force due to lift.The power effect of turbo-prop engine increases the CLmaxof different flap configurations up to 8% ~9%.The CLmaxobtained by the method agrees well with flight test,and the low speed performance of the aircraft is well utilized.The research can be an important reference for propeller aircraft design.
propeller aircraft;slipstream;maximum lift coefficient;normal force;wind tunnel test
V211.71
:Adoi:10.7638/kqdlxxb-2014.0024
0258-1825(2015)05-0655-06
2014-04-17;
:2014-08-06
刘毅*(1982-),男,硕士,工程师,研究方向:飞机气动设计.E-mail:evanliuyi@hotmail.com
刘毅,赵晓霞,欧阳绍修.螺旋桨飞机升力失速特性研究[J].空气动力学学报,2015,33(5):655-660.
10.7638/kqdlxxb-2014.0024 Liu Y,Zhao X X,Ouyang S X.Investigation on lift stall characteristics of propeller aircraft[J].Acta Aerodynamica Sinica,2015,33(5):655-660.