文章快速检索  
  高级检索
基于后缘襟翼偏转的大型客机变弯度技术减阻收益
何萌1, 杨体浩1, 白俊强2, 杨一雄1     
1. 西北工业大学 航空学院, 西安 710072;
2. 西北工业大学 无人系统技术研究院, 西安 710072
摘要: 为了满足工程实际约束,针对宽体客机内外襟翼位置与偏角卡位进行了机翼后缘变弯度减阻收益研究。使用雷诺-平均Navier-Stokes(RANS)方程对襟翼不同后缘偏角采用遍历的方式进行了气动力评估,得到后缘襟翼最佳偏角;探究了变弯度技术在非设计点的减阻收益,以及变弯度技术对宽体客机阻力发散和抖振边界设计要求的拓展能力,进一步采用远场阻力分解方法探究了设计结果的减阻机理。结果表明,在变马赫数的非设计点,考虑俯仰力矩系数配平之前,阻力能获得一些收益,但考虑俯仰力矩系数配平后,变弯度后的阻力系数不减反增;当升力系数发生变化时,变弯度在考虑俯仰力矩配平的情况下,均能取得一定的收益;变弯度技术也减缓了抖振点激波诱导分离的趋势,对抖振特性有较好的改善作用。基于后缘襟翼偏转的变弯度减阻收益评估和机理分析,能为宽体客机机翼变弯度设计提供参考。
关键词: 大型客机设计    后缘变弯度    力矩配平    阻力发散    抖振    
Drag reduction benefits of variable camber technology of airliner based on trailing-edge flap deflection
HE Meng1, YANG Tihao1, BAI Junqiang2, YANG Yixiong1     
1. School of Aeronautics, Northwestern Polytechnical University, Xi'an 710072, China;
2. Unmanned System Research Institute, Northwestern Polytechnical University, Xi'an 710072, China
Abstract: To satisfy the practical constraints of engineering projects, this study examines the drag reduction caused by morphing trailing-edge wings based on the position of the inner and outer flaps of the wide-body aircraft and the deflection position. Aerodynamic evaluations based on the Reynolds-Averaged Navier-Stokes (RANS) equation are conducted on all the deflection angles of the flaps, obtaining the optimal deflection angle of the trailing edge flaps. Benefits of the variable camber technology are investigated, including the expansion of the design requirements for the buffet and drag-divergence boundary as well as the drag reduction at other non-design points. Furthermore, the far field drag decomposition method is used to explore the drag reduction mechanism of the design results. Results show that at the non-design points of the variable Mach number, the drag reduction can be obtained regardless of the trimming moment. However, the drag is not decreased while considering the trimming moment. With different cruise lift coefficients, the variable camber technology can invariably achieve certain drag reduction when considering the trimming moment. Moreover, the variable camber technology can also slow the trend of shock-induced separation and improve the buffeting characteristics of the wing. The drag reduction evaluation and mechanism analysis of the variable camber technology based on trailing edge flap deflection provide reference for the variable camber design of wide-body airliner wings.
Keywords: airliner design    trailing-edge variable camber    trimming moment    drag-divergence    buffet    

随着人们对绿色生活环境关注度的增加以及燃油价格的不断上涨,未来民航客机的经济性和环保性越来越重要[1]。机翼后缘变弯度的设计增强了机翼适应各种飞行条件的能力, 使工程师能在各工况下设计机翼不同的后缘偏角以改善其气动性能, 而不是一系列飞行条件下的气动性能之间的折衷,是一种有效的减小燃油消耗的技术。

Boeing和NASA共同开展了“可变弯度后缘襟翼”(Variable Camber Continuous Trailing Edge Flap,VCCTEF)项目研究[2-3],探索了变弯度技术用于未来质量更轻、机翼柔性更大的飞行器的性能优势,比如提高巡航效率、扩宽抖振边界、增加机动性等。Boeing757的原型机GTM(Generic Transport Model)构型上使用此项技术,有效减小了其巡航阻力[4-5]。Boeing787也采用了机翼后缘变弯度设计,所减少的巡航阻力相当于节省340~450 kg(750~1 000 lb)的重量[6]

变弯度的结构设计[7-9]目前具有较高的技术储备,比如FlexSys FlexFoil装置[10-11]改变襟翼偏角时,控制表面仍能呈现出平滑和连续的过度。国内外学者对机翼后缘变弯度技术在气动设计中的收益也已经开展了广泛的研究。Molinari等[12-13]和Lee等[14]使用低精度气动模型研究了后缘连续变弯的优点;Lyu和Martins[15-16]使用高精度的伴随求解器对CRM机翼后缘连续变弯度技术的收益进行研究,对407个可能的飞行状态进行后缘变弯度的优化,生成一个优化构型的数据库,分析了后缘变弯度技术对燃油消耗的减小量,但是没有考虑到配平阻力对减阻收益的影响。Burdette等[17-18]使用气动结构求解器对后缘变弯度的收益进行了研究,结果表明,可变弯的机翼有1.02%的燃油收益,考虑结构变形的影响,燃油可减小1.72%。梁煜[19]和陈钱[20]等以翼型为例,研究了翼型后缘变弯度对气动性能与压力分布的影响。郭同彪等[21-22]研究了后缘连续变弯度对跨声速翼型气动特性的影响,以及针对民用客机机翼-机身-平尾构型开展了后缘连续变弯度机翼气动优化设计。

之前的工作主要针对二维翼型以及机翼全翼展连续变弯对气动特性的影响研究,全翼展连续变弯是机翼成型各控制剖面后缘均可变弯,设计变量较多,且对机翼变形材料要求较高,在规律研究上具有指导意义。但是针对大型客机,连续变弯技术的技术成熟度和可靠性有待进一步提高,目前尚无法在工程实际中得到应用。波音公司为787飞机测试和开发了一种后缘可变弯度(TEVC)系统,在巡航时通过操纵后缘内外襟翼以0.5°为增量进行偏转[6]。本文考虑了工程实际约束的影响,只对内外襟翼进行变弯偏转(具有两个设计变量),对宽体客机进行了减阻收益评估。在考虑配平阻力的影响下探究了变弯度技术对阻力发散和抖振边界设计要求的拓展能力,以及变弯度技术在非设计点的减阻收益;内外襟翼偏角不同可能会使气动性能的改善效果更好,但会导致机翼后缘变弯机构更加复杂,因此对比了襟翼同偏和差偏在气动上的减阻收益;分析了在不同非设计点机翼后缘变弯度技术的压力分布变化规律、载荷分布变化规律,并进一步采用远场阻力分解方法探究了变弯度技术的减阻机理。

1 三维机翼后缘变弯度建模

变弯度时使用如图 1所示的襟翼平面形状,具有两个设计变量,分别为内襟翼偏角和外襟翼偏角。采用定轴偏转的方式偏转襟翼,偏转前后的效果对比如图 2所示,可以看到,变弯度后的翼面依然保持着连续、光滑的特征。同时由于机构的行程和精度限制,要求襟翼以0.5°的卡位运行,且最大偏角不超过2°,规定襟翼下偏为正,上偏为负。平尾也作为一个额外的自由度参与考虑,即通过平尾安定面的偏转来配平俯仰力矩, 使变弯度构型的俯仰力矩与基础构型一致。

图 1 襟翼平面形状 Fig. 1 Flap plane shape
图 2 变弯度偏转效果 Fig. 2 Variable camber deflection effect

由于挂架只对局部的流场产生影响,对整体气动特性影响不大[23-24],且为了减小网格量和减小流场的复杂性,建模时没有考虑挂架的影响[25]

2 数值求解方法

本文使用雷诺-平均Navier-Stokes (RANS)方程求解器进行定常求解。全湍流计算,使用的湍流模型为剪切应力输运(SST)湍流模型。选取CRM翼身组合体构型[26]对求解器精度进行校验,网格量为1 725万,计算状态为:马赫数Ma=0.85, 雷诺数Re=5×106, 升力系数CL=0.519。

4个控制剖面的压力系数Cp分布计算结果与试验数据的对比如图 3所示,求解器计算得到的压力分布形态与试验值吻合良好,尤其是下表面贴合较好。对比结果表明,本文采用的求解器具有可满足工程应用需求的计算精度和可靠性。宽体构型计算时采用多块结构网格,网格拓扑如图 4所示。

图 3 CRM构型计算结果与试验数据对比 Fig. 3 Comparison of CRM configuration calculation results with experimental data
图 4 计算网格拓扑 Fig. 4 Grid topology

为了明确网格量对计算结果的影响,本文计算了2 000万、3 000万和4 300万网格量下阻力系数的数值,计算结果如表 1所示。使用2 000万网格计算时,相较于3 000万的网格量,计算出来的阻力值增加7.5 counts(1 count=0.000 1),使用4 300万的网格相较于3 000万的网格,阻力值只增加了约2 counts。从工程应用研究角度,可认为3 000万网格与4 300万网格具有相当的计算精度水平。考虑到本文计算量大,为了兼顾计算时间和计算精度,最终采用3 000万的网格量用于变弯度技术减阻收益的评估。

表 1 基础构型不同网格量计算结果 Table 1 Computation results of different gridsizes for original configuration
网格编号 网格量/万 升力系数CL 阻力系数CD/count
L0 4 300 0.48 228.80
L1 3 000 0.48 230.76
L2 2 000 0.48 238.29
3 基础构型气动性能

表 2给出了未考虑襟翼变弯的基础构型气动特性,设计马赫数Ma=0.85,设计升力系数CL=0.48,雷诺数Re=4.7×107。巡航设计点力矩基本配平,巡航攻角为2.48°,阻力系数为230.98 counts,升阻比为20.78。马赫数0.87与马赫数0.85的阻力差量为17.53 counts,阻力发散特性较好。图 5为马赫数0.85机翼表面压力分布及流线图。

表 2 基础构型气动性能 Table 2 Aerodynamic performance of original configuration
马赫数Ma 升力系数CL 阻力系数CD/count 力矩系数Cmy 升阻比K
0.85 0.48 230.98 0.007 1 20.78
0.87 0.48 248.51 0.002 7 19.32
图 5 马赫数0.85机翼表面压力系数分布及流线图 Fig. 5 Pressure coefficient distribution and surface streamline of wing(Ma =0.85)
4 采用变弯度技术的状态选取

变弯度技术能够在巡航过程中通过改变后缘偏角,在非设计点状态下获得减阻收益,从而降低整个巡航段的平均阻力。因此,本文选取的变弯度技术减阻收益评估状态包括以下3类:

1) 随着飞行过程中燃油不断消耗,飞机的升力系数也会随时间发生变化[17, 27]。故选取了升力系数变化10%和20%的非设计点,即小升力系数点(CL=0.432和CL=0.384)和大升力系数点(CL=0.528和CL=0.576)。在此基础上,进一步分析整个升力系数变化范围内的减阻规律。

2) 受航空管制以及变高度巡航策略等因素的影响,民用客机难以在整个巡航阶段内都以一个固定的设计马赫数进行巡航,其巡航马赫数往往会在设计马赫数附近小幅变化[17]。故选取了马赫数偏小点(Ma=0.83,CL=0.48)和阻力发散点(Ma=0.87,CL=0.48)进行变弯度技术的减阻收益研究。

3) 研究变弯度技术对升力系数增加所引起的高速抖振特性的影响。对于大型民用客机,设计规范要求其具有1.3 g的抖振裕量,故重点研究了变弯度技术对抖振点(Ma=0.85,CL=0.624)气动特性的改善作用。

最终,选取的计算状态的分布如图 6所示。

图 6 变弯度设计的计算状态分布 Fig. 6 Calculation condition distribution of variable camber design
5 变升力系数状态后缘变弯技术减阻收益 5.1 小升力系数下的减阻收益

对内外襟翼的可能偏角进行遍历计算,表 3给出了CL=0.384状态下基础构型、不考虑力矩配平阻力的最佳襟翼偏角构型、考虑配平阻力的最佳襟翼同偏构型和最佳襟翼差动偏转构型结果。襟翼同偏表示内、外襟翼的偏角一致,差动偏转表示内、外襟翼的偏角不一致。

表 3 基础构型与最优襟翼偏角构型气动力系数(CL=0.384) Table 3 Aerodynamic coefficients of original configuration and optimal flap declination configuration(CL=0.384)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的ΔCD, modi/count
0 0 206.4 0.028 12 0 0 0
-1.5 -0.5 204.1 0.016 21 -2.3 -0.011 91 -1.6
-1 -1 204.5 0.031 78 -1.9 0.003 66 -2.1
-1 -1.5 204.8 0.040 45 -1.6 0.012 33 -2.3

图 7给出了CL=0.384时的计算结果,Inner表示内襟翼偏角,Outer表示外襟翼偏角,可以看到小升力系数时襟翼上偏有助于阻力减小,内襟翼上偏,外襟翼下偏,造成低头力矩增加,而襟翼同偏对力矩系数的影响不大。不考虑力矩配平阻力的最优襟翼偏角为内襟翼-1.5°,外襟翼-0.5°,减阻2.3 counts;考虑配平阻力后仅减阻1.6 counts;考虑配平阻力下最佳同偏襟翼偏角为-1°,减阻2.1 counts,最佳差动偏转襟翼偏角为内襟翼-1°,外襟翼-1.5°,减阻2.3 counts。

图 7 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(CL=0.384) Fig. 7 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(CL=0.384)

表 4图 8给出了CL=0.432的计算结果,不考虑力矩配平阻力的最佳襟翼偏角为内襟翼-1.0°,外襟翼-0.5°,减阻1.4 counts; 考虑配平阻力后仅减阻1.2 counts;考虑配平阻力下最佳襟翼偏角为内襟翼-0.5°,外襟翼-0.5°,减阻1.4 counts。可见,是否考虑力矩配平对最佳襟翼偏角和减阻量有较大影响,在进行变弯度减阻收益研究时有必要考虑配平阻力的影响。

表 4 基础构型与最优襟翼偏角构型气动力系数(CL=0.432) Table 4 Aerodynamic coefficients of original configuration and optimal flap declination configuration(CL=0.432)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的ΔCD, modi/count
0 0 218.6 0.018 35 0 0 0
-1.0 -0.5 217.2 0.013 45 -1.4 -0.004 90 -1.2
-0.5 -0.5 217.3 0.020 42 -1.3 0.002 06 -1.4
图 8 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(CL=0.432) Fig. 8 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(CL=0.432)

图 9图 10给出了变弯前后机翼上激波强度对比、环量分布对比及剖面压力系数Cp对比,Original表示基础构型,TEVC表示变弯构型(襟翼同偏最优结果)。对于小升力系数的计算工况,内外襟翼均上偏使得机翼后缘卸载,巡航攻角增大,引起激波位置后移,载荷前移。机翼载荷前移导致抬头力矩增加。同时机翼中后缘部分载荷减小,使压力恢复平缓,激波强度减小。通过机翼展向环量分布对比可以看出,机翼襟翼的同时上偏使载荷外移,更加贴近椭圆形环量分布。

图 9 变弯前后对比(内襟翼-1°,外襟翼-1°,CL=0.384) Fig. 9 Comparison of initial and flap deflection configurations(inner flap -1°, outer flap -1°, CL=0.384)
图 10 变弯前后对比(内襟翼-0.5°,外襟翼-0.5°,CL=0.432) Fig. 10 Comparison of initial and flap deflection configurations(Inner flap -0.5°, Outer flap -0.5°, CL=0.432)
5.2 大升力系数下减阻收益

表 5图 11给出了CL=0.528的计算结果,从图 11可以看出,大升力系数下,襟翼下偏有助于阻力减小,襟翼偏角变化对力矩系数的影响与小升力系数时基本一致。不考虑力矩配平阻力的最佳襟翼偏角为内襟翼0.5°,外襟翼1.0°,减阻5.5 counts; 考虑配平阻力后仅减阻4.2 counts;考虑配平阻力的最佳同偏襟翼偏角为1°,减阻5.0 counts,最佳差动偏转襟翼偏角为内襟翼1.5°,外襟翼0.5°,减阻5.7 counts。

表 5 基础构型与最优襟翼偏角构型气动力系数(CL=0.528) Table 5 Aerodynamic coefficients of original configuration and optimal flap declination configuration(CL=0.528)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的ΔCD, modi/count
0 0 258.2 -0.006 51 0 0 0
0.5 1 252.7 -0.015 12 -5.5 -0.008 61 -4.2
1 1 253.0 -0.007 99 -5.2 -0.001 48 -5.0
1.5 0.5 254.4 0.005 90 -3.9 0.012 41 -5.7
图 11 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(CL=0.528) Fig. 11 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(CL=0.528)

表 6图 12给出了CL=0.576的计算结果,不考虑力矩配平阻力的最佳襟翼偏角为内襟翼1.0°,外襟翼1.5°,减阻14.0 counts;考虑配平阻力后仅减阻12.2 counts;考虑配平阻力的最佳同偏襟翼偏角为1.5°,减阻12.9 counts,最佳差动偏转襟翼偏角为内襟翼2.0°,外襟翼1.0°,减阻13.5 counts。大升力系数下,变弯度技术具有较为明显的气动性能的改善。配平阻力依然对结果有较大影响,在进行变弯度减阻收益研究时有必要考虑配平阻力的影响。

表 6 基础构型与最优襟翼偏角构型气动力系数(CL=0.576) Table 6 Aerodynamic coefficients of original configuration and optimal flap declination configuration(CL=0.576)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的CD, modi/count
0 0 295.4 -0.020 26 0 0 0
1 1.5 281.4 -0.029 41 -14.0 -0.009 14 -12.2
1.5 1.5 282.1 -0.022 04 -13.3 -0.001 78 -12.9
2 1 284.2 -0.008 82 -11.2 0.011 44 -13.5
图 12 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(CL=0.576) Fig. 12 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(CL=0.576)

图 13图 14给出了变弯前后机翼上激波强度对比、环量分布对比及剖面压力系数Cp对比。大升力系数下,内外襟翼均下偏使得机翼后缘加载,上翼面的压力分布会变得饱满,压力恢复平缓,激波强度降低,波阻减小。通过机翼展向环量分布对比可以看出,机翼襟翼的同时下偏使载荷内移,环量分布更加偏离椭圆形环量分布。

图 13 襟翼同偏变弯前后对比(内襟翼1.0°,外襟翼1.0°,CL=0.528) Fig. 13 Comparison of initial and flap deflection configurations(Inner flap 1.0°, Outer flap 1.0°, CL=0.528)
图 14 襟翼同偏变弯前后对比(内襟翼1.5°,外襟翼1.5°,CL=0.576) Fig. 14 Comparison of initial and flap deflection configurations(Inner flap 1.5°, Outer flap 1.5°, CL=0.576)
6 变马赫数状态减阻收益 6.1 马赫数偏小(Ma=0.83)时的减阻收益

表 7图 15可看出在不考虑力矩配平阻力的影响下,内襟翼偏转-1°,外襟翼偏转0.5°,减阻0.5 counts,但此偏转条件下,低头力矩增加,考虑配平阻力后,阻力增加0.5 counts。

表 7 基础构型与最优襟翼偏角构型气动力系数(Ma=0.83) Table 7 Aerodynamic coefficients of original configuration and optimal flap declination configuration(Ma=0.83)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的ΔCD, modi/count
0 0 236.8 0.007 27 0 0 0
-1 0.5 236.3 -0.014 30 -0.5 -0.021 56 0.5
图 15 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(Ma=0.83) Fig. 15 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(Ma=0.83)

对于马赫数偏小的非设计点,首先由于激波位置在变弯区域之前且距离太远,变弯控制的效果较弱, 没有较为明显的减阻效果。其次最优偏角为内襟翼上偏,外襟翼下偏,导致机翼载荷外移,低头力矩增加,造成配平阻力增加。因此考虑力矩配平的影响之后,变弯度技术对马赫数偏小的非设计点没有带来气动特性的改善。

6.2 阻力发散点(Ma=0.87)的减阻收益

表 8图 16可看出在不考虑力矩配平阻力的影响下,内襟翼偏转-1°,外襟翼偏转0.5°,减阻1.4 counts,但此偏转条件下,低头力矩增加,考虑配平阻力后,阻力增加2.4 counts。

表 8 基础构型与最优襟翼偏角构型气动力系数(Ma=0.87) Table 8 Aerodynamic coefficients of original configuration and optimal flap declination configuration(Ma=0.87)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的ΔCD, modi/count
0 0 248.5 0.002 72 0 0 0
-1 0.5 247.1 -0.019 24 -1.4 -0.021 95 2.4
图 16 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(Ma=0.87) Fig. 16 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(Ma=0.87)

对于阻力发散点,首先由于激波位置靠后,以至于在变弯度过渡区之后,变弯控制的效果也很弱。其次考虑力矩的影响之后,配平阻力的增加使变弯度技术对阻力发散点也没有带来气动特性的改善。

在不同的马赫数下,不论是马赫数增大还是减小,内襟翼上偏、外襟翼下偏导致机翼载荷外移,低头力矩增加,考虑力矩系数后,变弯度的减阻收益均为0。

7 变弯度技术对抖振特性的改善作用

表 9图 17给出了CL=0.624抖振点的计算结果,考虑力矩配平阻力时,最佳同偏襟翼偏角为1.5°,减阻18.8 counts,最佳差动偏转襟翼偏角为内襟翼2.0°,外襟翼1.0°,减阻21.2 counts。

表 9 基础构型与最优襟翼偏角构型气动力系数(CL=0.624) Table 9 Aerodynamic coefficients of original configuration and optimal flap declination configuration(CL=0.624)
内襟翼偏角/(°) 外襟翼偏角/(°) 阻力系数CD/count 俯仰力矩系数Cmy ΔCD/count ΔCmy 考虑力矩后的ΔCD, modi/count
0 0 347.5 -0.028 89 0 0 0
1.5 1.5 329.0 -0.028 70 -18.5 0.000 19 -18.8
2 1 330.0 -0.017 28 -17.6 0.011 61 -21.2
图 17 不同内外襟翼偏角下的气动阻力以及俯仰力矩系数变化云图(CL=0.624) Fig. 17 Contour of aerodynamic drag and pitching moment coefficient at different inner and outer flaps(CL=0.624)

图 18给出了抖振点襟翼同偏变弯前后机翼剖面压力系数Cp对比。内外襟翼下偏,机翼后缘载荷增加。内翼段通过襟翼下偏,压力恢复平缓,呈无激波形态,中外翼段通过外襟翼下偏,激波位置后移,激波强度减小。

图 18 襟翼同偏变弯前后机翼剖面压力系数对比(内襟翼1.5°,外襟翼1.5°,CL=0.624) Fig. 18 Comparison of wing section pressure coefficient of original and flap deflection configurations(Inner flap 1.5°, Outer flap 1.5°, CL=0.624)

抖振是激波附面层干扰引起分离流产生的非定常气动力引起的振动现象[28]图 19为变弯前和变弯后的抖振特性对比,给出了一系列升力系数下的机翼上表面压力云图与极限流线。变弯前构型在CL=0.594和CL=0.604下的分离区大小分别接近于变弯后构型在CL=0.624和CL=0.634下的分离区大小,结果表明变弯度可改善抖振特性。变弯度通过调整载荷分布,明显降低中外翼激波强度,减缓了激波诱导分离的趋势。

图 19 变弯前和变弯后抖振特性对比 Fig. 19 Comparison of buffeting characteristics before and after flap deflection
8 变弯度技术非设计点减阻机理分析

图 20给出了不同升力系数下变弯度的减阻情况(考虑力矩系数配平后), 其中蓝线和绿线分别给出了同偏与差动偏转下的结果。可以看到,偏离设计点升力系数越远,变弯度获得的减阻收益越大,且大升力系数时的收益要大于小升力系数时的收益。对比同偏和差动偏转,两者减阻量随升力系数变化的变化趋势基本一致,差动偏转的收益略高于同偏的收益。

图 20 不同升力系数下变弯度的减阻量(考虑力矩系数配平后) Fig. 20 Drag reduction of variable camber under different lift coefficients(considering trimming moment coefficient)

进一步采用远场阻力分解方法[29-30]探究了变弯度技术的减阻来源,远场阻力分解方法是基于线性动量定理关系进行阻力计算,将阻力按照产生的机制分为黏性阻力、激波阻力、诱导阻力和伪阻力。

图 21给出了不同升力系数下变弯前后激波阻力值对比,小升力系数状态下波阻减小量不足1 count。大升力系数状态下,波阻有明显的减小。波阻的减小主要是因为后缘襟翼变弯度技术会调整载荷分布,改变激波位置和强度。

图 21 变弯前后激波阻力对比 Fig. 21 Comparison of shock drag of original and flap deflection configurations

图 22给出了不同升力系数下变弯前后诱导阻力值对比。在不同升力系数下,诱导阻力值没有明显改善。从之前的环量分布对比可以看出,在大升力系数时,载荷内移,更加偏离椭圆形环量分布,导致诱导阻力略微增加。

图 22 变弯前后诱导阻力对比 Fig. 22 Comparison of induced drag of original and flap deflection configurations
9 结 论

1) 变弯度引起全机俯仰力矩变化带来的配平阻力,对变弯度减阻收益影响明显,考虑配平阻力损失是进行变弯度机翼设计的关键。

2) 在整个巡航升力系数范围内的非设计点,变弯度均能取得一定减阻收益,大巡航升力系数下减阻收益更大。不同巡航升力系数下变弯度减阻机理存在明显差异。小巡航升力系数下诱导阻力以及激波阻力都得到不同程度的减小。大巡航升力系数下以适当增大诱导阻力的方式,显著削弱激波阻力。

3) 变弯度通过将机翼载荷内移,削弱激波强度,减小激波诱导的附面层分离,拓展抖振边界裕度。

4) 变弯度无法改善阻力发散点和阻力蠕增特性。

5) 对比差动偏转和同偏,两者的减阻收益随升力系数的变化趋势一致,但在某些情况下,差动偏转减阻收益更明显。

参考文献
[1] 陈迎春, 张美红, 张淼, 等. 大型客机气动设计综述[J]. 航空学报, 2019, 40(1): 522759.
CHEN Y C, ZHANG M H, ZHANG M, et al. Review of large civil aircraft aerodynamic design[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(1): 522759. (in Chinese)
Cited By in Cnki (5) | Click to display the text
[2] KAUL U K, NGUYEN N T. A 3-D drag optimization study of variable camber continuous trailing edge flap (VCCTEF) using OVERFLOW: AIAA-2015-2422[R]. Reston: AIAA, 2015.
Click to display the text
[3] URNES J, NGUYEN N, IPPOLITO C, et al. A mission adaptive variable camber flap control system to optimize high lift and cruise lift to drag ratios of future N+3 transport aircraft: AIAA-2013-0214[R]. Reston: AIAA, 2013.
Click to display the text
[4] IPPOLITO C, NGUYEN N, TOTAH J. Initial assessment of a variable-camber continuous trailing-edge flap system on a rigid wing for drag reduction in subsonic cruise: AIAA-2013-5143[R]. Reston: AIAA, 2013.
Click to display the text
[5] KAUL U K, NGUYEN N T. Drag optimization study of variable vamber continuous trailing edge flap (VCCTEF) using OVERFLOW: AIAA-2014-2444[R]. Reston: AIAA, 2014.
Click to display the text
[6] GUY N. Boeing unveils plans for trailing edge variable camber on 787 to reduce drag, save weight[EB/OL].[2019-09-10]. 2006, http://Tlightglobal.com/news/ar-ticles/.
[7] SOFLA A Y N, MEGUID S A, TAN K T, et al. Shape morphing of aircraft wing:Status and challenges[J]. Materials & Design, 2010, 31(3): 1284-1292.
Click to display the text
[8] 杨智春, 解江. 柔性后缘自适应机翼的概念设计[J]. 航空学报, 2009, 30(6): 1028-1034.
YANG Z C, XIE J. Concept design of adaptive wing with flexible trailing edge[J]. Acta Aeronautica et Astronautica Sinica, 2009, 30(6): 1028-1034. (in Chinese)
Cited By in Cnki (38) | Click to display the text
[9] 尹维龙, 石庆华, 田东奎. 变体后缘的索网传动机构设计与分析[J]. 航空学报, 2013, 34(8): 1824-1831.
YIN W L, SHI Q H, TIAN D K. Design and analysis of transmission mechanism with cable networks for morphing trailing-edge[J]. Acta Aeronautica et Astronautica Sinica, 2013, 34(8): 1824-1831. (in Chinese)
Cited By in Cnki (10) | Click to display the text
[10] KOTA S, HETRICK J, OSBORN R, et al. Design and application of compliant mechanisms for morphing aircraft structures[J]. Proceedings of SPIE, 2003, 5054: 24-33.
Click to display the text
[11] KOTA S, OSBORN R, ERVIN G, et al. Mission adaptive compliant wing-design, fabrication and flight tesT[C]//Proceedings of the Rto Applied Vehicle Technology Panel(AVT) Symposium. Brussels: The NATO Science and Technology Organization (STO), 2009.
Click to display the text
[12] MOLINARI G, QUACK M, DMITRIEV V, et al. Aero-structural optimization of morphing airfoils for adaptive wings[J]. Journal of Intelligent Material Systems and Structures, 2011, 22(10): 1075-1089.
Click to display the text
[13] MOLINARI G, QUACK M, ARRIETA A, et al. Design, realization and structural testing of a compliant adaptable wing[J]. Smart Materials and Structures, 2015, 24(10): 105027.
Click to display the text
[14] LEE D, GONZALEZ L F, PERIAUX J, et al. Robust aerodynamic design optimisation of morphing aerofoil/wing using distributed MOGA[C]//Proceedings of the 28th Congress of the International Council of the Aeronautical Sciences. Bonn: The International Council of the Aeronautical Sciences, 2012.
Click to display the text
[15] LYU Z. High-fidelity aerodynamic design optimization of aircraft configurations[D]. Ann Arbor: University of Michigan, 2014: 104-122.
Click to display the text
[16] LYU Z, MARTINS J R R A. Aerodynamic shape optimization of an adaptive morphing trailing edge wing[J]. Journal of Aircraft, 2015, 52(6): 1951-1970.
Click to display the text
[17] BURDETTE D A, KENWAY G K, MARTINS J. Aerostructural design optimization of a continuous morphing trailing edge aircraft for improved mission performance[C]//17th AIAA/ISSMO Multidisciplinary Analysis and Optimization Conference. Reston: AIAA, 2016.
Click to display the text
[18] BURDETTE D A. High-fidelity aerostructural design optimization of transport aircraft with continuous morphing trailing edge technology[D]. Ann Arbor: University of Michigan, 2017: 119-123.
Click to display the text
[19] 梁煜, 单肖文. 大型民机翼型变弯度气动特性分析与优化设计[J]. 航空学报, 2016, 37(3): 790-798.
LIANG Y, SHAN X W. Aerodynamic analysis and optimization design for variable camber airfoil of civil transport jet[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(3): 790-798. (in Chinese)
Cited By in Cnki (10) | Click to display the text
[20] 陈钱, 白鹏, 尹维龙, 等. 可连续光滑偏转后缘的变弯度翼型气动特性分析[J]. 空气动力学学报, 2010, 28(1): 46-53.
CHEN Q, BAI P, YIN W L, et al. Analysis on the aerodynamic characteristics of variable camber airfoil with continuous smooth morphing trailing edge[J]. Acta Aerodynamica Sinica, 2010, 28(1): 46-53. (in Chinese)
Cited By in Cnki (48) | Click to display the text
[21] 郭同彪, 白俊强, 杨体浩. 后缘连续变弯度对跨声速翼型气动特性的影响[J]. 航空学报, 2016, 37(2): 513-521.
GUO T B, BAI J Q, YANG T H. Influence of continuous trailing-edge variable camber on aerodynamic characteristics of transonic airfoils[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(2): 513-521. (in Chinese)
Cited By in Cnki (12) | Click to display the text
[22] 郭同彪, 白俊强, 李立, 等. 民用客机变弯度机翼优化设计[J]. 中国科学:技术科学, 2018, 48(1): 55-66.
GUO T B, BAI J Q, LI L, et al. The morphing trailing-edge wing optimization design of the civil aircraft[J]. Scientia Sinica Technologica, 2018, 48(1): 55-66. (in Chinese)
Cited By in Cnki (6) | Click to display the text
[23] 周洪升.翼下吊挂式短舱造型方法及其性能分析[D], 南京: 南京航空航天大学, 2009: 35-42.
ZHOU H S. Parameter model design and numerical simulation of under-the-wing nacellle for civil aircraft[D]. Nanjing: Nanjing University of Aeronautics and Astronautics, 2009: 35-42(in Chinese).
Cited By in Cnki | Click to display the text
[24] 党铁红. 翼吊布局民用飞机发动机安装设计[J]. 民用飞机设计与研究, 2008(2): 8-14.
DANG T H. Wing crane layout civil aircraft engine installation design[J]. Civil Aircraft Design and Research, 2008(2): 8-14. (in Chinese)
Cited By in Cnki (19) | Click to display the text
[25] 张冬云, 张美红, 王美黎, 等. 翼吊布局民机短舱位置气动影响[J]. 空气动力学学报, 2017, 35(6): 781-786.
ZHANG D Y, ZHANG M H, WANG M L, et al. Aerodynamic influence of nacelle position of a wing-mounted civil aircraft[J]. Acta Aerodynamica Sinica, 2017, 35(6): 781-786. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[26] BRODERSEN O, CRIPPA S. RANS-based aerodynamic drag and pitching moment predictions for the common research model[M]. Springer International Publishing, 2014: 485-493.
[27] LIEM R P, KENWAY G K W, MARTINS J R R A. Multimission aircraft fuel-burn minimization via multipoint aerostructural optimization[J]. AIAA Journal, 2015, 53(1): 104-122.
Click to display the text
[28] 高传强, 张伟伟. 机翼跨声速抖振数值模拟及模态分析[J]. 航空学报, 2019, 40(7): 122597.
GAO C Q, ZHANG W W. Numerical simulation and modal analysis of transonic buffet flow over wings[J]. Acta Aeronautica et Astronautica Sinica, 2019, 40(7): 122597. (in Chinese)
Cited By in Cnki (2) | Click to display the text
[29] VANDERVOOREN J, SLOOFF J W. CFD-based drag prediction: State-of-the-art, theory, prospects: T-1990-247[R]. Reston: AIAA, 1990.
[30] 何小龙.基于流场区域分解的远场阻力分解方法[C]//第一届中国空气动力学大会, 2018.
HE X L. Far field drag decomposition method based on flow field decomposition[C]//1st Chinese Conference of Aerodynamics, 2018(in Chinese).
http://dx.doi.org/10.7527/S1000-6893.2020.23462
中国航空学会和北京航空航天大学主办。
0

文章信息

何萌, 杨体浩, 白俊强, 杨一雄
HE Meng, YANG Tihao, BAI Junqiang, YANG Yixiong
基于后缘襟翼偏转的大型客机变弯度技术减阻收益
Drag reduction benefits of variable camber technology of airliner based on trailing-edge flap deflection
航空学报, 2020, 41(7): 123462.
Acta Aeronautica et Astronautica Sinica, 2020, 41(7): 123462.
http://dx.doi.org/10.7527/S1000-6893.2020.23462

文章历史

收稿日期: 2019-09-10
退修日期: 2020-02-19
录用日期: 2020-03-31
网络出版时间: 2020-04-13

相关文章

工作空间