2019, Vol. 37
20世纪80年代，美国B2隐身轰炸机的问世标志着飞翼布局发展到了成熟阶段。近年来，飞翼布局以其特有的气动、结构和隐身优势，受到了无人机设计领域的高度重视[1]。美国是研究飞翼布局较早的国家之一，相继推出了X-45、X-47等系列低可探测性无人作战飞机研制计划[2-3]。与常规布局相比飞翼布局隐身性能好，巡航升阻比高，结构重量轻[4-5]；但是，飞翼布局飞机没有常规的平尾、垂尾，航向稳定性和操纵性较差[6-7]，俯仰操纵力臂较短，要获得与常规布局等效的俯仰操纵效果，需付出较大的升力损失[8-9]，飞翼布局飞机在起飞着陆阶段增升装置下偏增升的同时，将会产生较大的附加低头力矩，不满足起降操纵要求[10]。
为了提高飞翼布局无人机的起降气动性能，本文引入了一种安装在无人机下表面的腹襟翼装置[11]，打开腹襟翼可以增加飞机起降阶段的实用升力系数和俯仰力矩系数，能够改善飞机的起降性能。
近年来，国内外对腹襟翼开展了初步研究。Yann D. Staelens和Ron F. Blackwelder等人分别从试验和计算的角度简单研究了腹襟翼对飞翼布局飞机起降性能的影响[12-13]。Daniel Dougherty、Adam Vore和Hamid Hefazi对腹襟翼安装位置和尺寸的影响进行了试验研究[14]。康辉论述了腹襟翼对飞翼布局无人机的影响[11]。北京航空航天大学的张曙光等人通过飞翼布局小型验证机的飞行试验，验证了腹襟翼能够提高飞翼布局无人机的着陆安全性[15]。西北工业大学的沈冬等人对腹襟翼进行了二维的相关研究工作[16]。本文采用基于有限体积法的RANS方程，研究分析了腹襟翼对飞翼布局无人机起降气动性能的影响，并对腹襟翼的安装位置、偏角和板面实度进行了详细的气动优化设计研究。研究结果表明，打开腹襟翼能够增加起飞离地迎角的升力系数和俯仰力矩系数，从而提高了飞翼布局无人机的起降气动特性，通过对腹襟翼关键参数的气动优化设计研究，得到了最优的参数组合。
1 腹襟翼的工作原理如图 1所示，腹襟翼是安装在飞机腹部的一种增升和俯仰控制装置。图 2给出了腹襟翼绕流的流场示意图，由图可见由于腹襟翼的存在，当气流流经飞机下表面时，使得腹襟翼前区域气流速度减小，因此该区域承受的压力增大，腹襟翼后区域流动出现分离，因此该区域承受的压力减小。
|
图 1 腹襟翼示意图 Fig.1 Diagram of the belly-flap |
|
图 2 腹襟翼绕流示意图 Fig.2 Diagram of flow around the belly-flap |
图 3给出了压力变化对气动力影响的示意图，飞机下表面腹襟翼前后压力的改变对飞机产生了附加的抬头力矩，这种附加的抬头力矩对飞机起降阶段的俯仰力矩配平以及俯仰控制是有利的。
|
图 3 腹襟翼工作原理示意图 Fig.3 Diagram of working principles of the belly-flap |
本文数值模拟研究采用的是基于有限体积法的三维定常RANS方程求解器，湍流模型选择的是SA模式。
三维积分形式的雷诺平均N-S方程：
| $ \frac{\partial }{{\partial t}}\int {\int {\int {\boldsymbol{Q}dV + \int {\int {\boldsymbol{f} \cdot \boldsymbol{n}{\rm{d}}} } } } } S = 0 $ |
式中：V为控制体体积，S为控制体表面积，Q为守恒量，f为通过表面S的无黏通量和黏性通量之和，n为控制体表面S的外法向单位矢量。
为了验证计算方法的可靠性，采用该方法对某民机翼身组合体风洞试验标模的气动特性进行了数值模拟研究，并将气动力特性的计算结果与实验数据进行了对比。模型的几何参数如表 1所示，模型的俯视图如图 4所示。
| 表 1 模型的几何参数 Table 1 Geometry parameters of the model |
|
|
图 4 某民机标模的俯视图 Fig.4 Top view of the civil aircraft model |
计算网格为采用ANSYS-ICEM软件生成的结构化网格，网格量约为1200万，模型表面网格和对称面部分网格如图 5所示，表 2为计算网格的分布情况。
|
图 5 模型的网格 Fig.5 Grid of the model |
| 表 2 网格的分布 Table 2 Distribution of grid |
|
|
低速计算模拟状态与风洞试验状态一致，计算迎角为-5°到13°。图 6和图 7分别为模型的升力系数和俯仰力矩系数计算结果与实验数据的对比，图中实线为本文RANS方程的计算结果，红色矩形符号为风洞试验结果。由气动力特性的对比结果可以看出，计算结果与试验数据吻合良好，验证了本文计算方法的可靠性。
|
图 6 升力系数计算值与试验数据的对比 Fig.6 Comparison of lift coefficient between calculated value and experimental data |
|
图 7 俯仰力矩系数计算值与试验数据的对比 Fig.7 Comparison of moment coefficient between calculated value and experimental data |
为了研究腹襟翼对飞翼布局飞机起降气动性能的影响，本文对如图 8所示三种计算模型进行了CFD数值模拟研究。其中，a构型(Model a)为干净构型；b构型(Model b)为起降构型；c构型(Model c)为腹襟翼打开的起降构型，c构型腹襟翼位于无人机下表面重心后距重心15%MAC(平均气动弦长)位置处，腹襟翼打开角度为90°。
|
图 8 本文用到的三种计算模型 Fig.8 Three models used in this paper |
计算网格同样为采用ANSYS-ICEM软件生成的结构化网格，以c构型为例，网格如图 9所示，网格量约为1000万，网格分布如表 3所示。
|
图 9 c构型的网格 Fig.9 Grid of the model c |
| 表 3 网格的分布 Table 3 Distribution of grid |
|
|
采用RANS方程求解器对如图 8所示的三种构型进行了低速数值模拟研究，计算马赫数为0.2，计算高度为标准海平面高度。图 10为随迎角变化，b构型和c构型相对a构型升力系数增量的变化曲线，图 11为随迎角变化，c构型相对b构型升力系数和俯仰力矩系数增量的变化曲线。
|
图 10 b构型和c构型升力系数增量的变化曲线 Fig.10 Curve of CL increment of Model b and Model c |
|
图 11 c构型升力系数增量和俯仰力矩系数增量的变化曲线 Fig.11 Curve of CL increment and Cm increment of Model c |
如图 10所示，与a构型相比，由于打开增升装置，b构型和c构型的升力系数明显增大；如图 11所示，与b构型相比，迎角小于10°时，c构型升力系数有所增大，俯仰力矩系数有所增加。在起飞离地迎角(8°)处，c构型升力系数比b构型大0.0244，c构型俯仰力矩系数比b构型大0.0254。因此，与起降构型相比，打开腹襟翼能够提高飞机起飞离地迎角的升力系数和俯仰力矩系数，对飞机的起飞着陆是有利的。
截取如图 12所示距对称面2 m的剖面，图 13和图 14分别为c构型和b构型在该剖面处的压力云图和流线图谱。与图 14相比，图 13腹襟翼前气动压力明显增强，腹襟翼后气动压力明显减弱，因此与b构型相比，c构型抬头力矩增大。图 15为c构型和b构型在该剖面处压力分布的对比，如图 15所示，同样可以得到c构型腹襟翼前后气动载荷的变化，且与蓝线相比，黑线腹襟翼前面积的增加量比腹襟翼后面积的减小量要大，因此，c构型比b构型升力系数大。
|
图 12 剖面位置示意图 Fig.12 Diagram of profile position |
|
图 13 c构型剖面压力云图和流线 Fig.13 Pressure clouds and streamlines on the profile of Model c |
|
图 14 b构型剖面压力云图和流线 Fig.14 Pressure clouds and streamlines on the profile of Model b |
|
图 15 剖面压力分布的对比 Fig.15 Comparison of profile pressure distribution |
为了寻求适用于本文飞翼布局无人机最优的腹襟翼，本文对腹襟翼的弦向位置、偏角和板面实度参数进行了气动优化设计研究，为保证计算结果的可对比性，计算网格与上文c模型的数量和分布情况完全一致。
4.1 弦向位置气动优化以起降构型为基准，如图 16所示分别将腹襟翼安装在无人机起降构型重心后15%MAC、20%MAC、25%MAC、30%MAC和35%MAC位置处，研究腹襟翼的安装位置对无人机低速气动性能的影响。图 17和图 18分别为随迎角变化，在不同位置安装上腹襟翼的起降构型相对原起降构型升力系数增量和俯仰力矩系数增量的变化曲线，图 19和图 20为8°迎角时相对起降构型，随腹襟翼位置的改变升力系数增量和俯仰力矩系数增量的变化曲线。如图 19和图 20所示随安装位置不断后移，当地迎角的升力系数增量不断增大，但是俯仰力矩系数增量不断减小。
|
图 16 腹襟翼安装位置示意图 Fig.16 Diagram of belly-flap locations |
|
图 17 升力系数增量的变化曲线 Fig.17 Curve of CL increment |
|
图 18 俯仰力矩系数增量的变化曲线 Fig.18 Curve of Cm increment |
|
图 19 升力系数增量随腹襟翼位置变化的关系曲线 Fig.19 Curve of CL increment with belly-flap locations |
|
图 20 俯仰力矩系数增量随腹襟翼位置变化的关系曲线 Fig.20 Curve of Cm increment with belly-flap locations |
由此可见, 腹襟翼安装在重心后15%MAC和20%MAC位置，升力增量有限；安装在30%MAC和35%MAC位置，俯仰力矩系数增量有限。因此，权衡升力系数和俯仰力矩系数的增量，重心后25%MAC处为本文飞翼布局无人机腹襟翼的最佳安装位置。
4.2 偏角气动优化选取重心后25%MAC处为腹襟翼的安装位置，研究腹襟翼的偏角对低速气动性能的影响，腹襟翼分别偏转30°、45°、60°和90°。图 21和图 22分别为随迎角变化，安装不同偏角的腹襟翼的起降构型相对
|
图 21 升力系数增量的变化曲线 Fig.21 Curve of CL increment |
|
图 22 俯仰力矩系数增量的变化曲线 Fig.22 Curve of Cm increment |
原起降构型升力系数增量和俯仰力矩系数增量的变化曲线，图 23和图 24为8°迎角时相对起降构型，随腹襟翼偏角的改变升力系数增量和俯仰力矩系数增量的变化曲线。由图 23和图 24得到，在8°迎角时，随着腹襟翼偏角不断增大，当地迎角的升力系数增量和俯仰力矩系数增量都不断增大。因此，90°为本文飞翼布局无人机腹襟翼的最佳偏角。
|
图 23 升力系数增量随腹襟翼偏角变化的关系曲线 Fig.23 Curve of CL increment with deflection angles |
|
图 24 俯仰力矩系数增量随腹襟翼偏角变化的关系曲线 Fig.24 Curve of Cm increment with deflection angles |
如图 25在腹襟翼半模板面上设计了两个缝、四个缝和八个缝三种不同板面实度的腹襟翼模型，并分别安装在起降构型重心后25%MAC处，其中每个缝的形状和面积是相同的，保证每一种腹襟翼打开角度都为90°，研究腹襟翼板面实度对无人机低速气动性能的影响。
|
图 25 半模腹襟翼板面实度设计模型 Fig.25 Design models of solidity of belly-flap half models |
图 26和图 27分别为随迎角变化，安装不同板面实度腹襟翼的起降构型相对原起降构型升力系数增量和俯仰力矩系数增量的变化曲线，图中solid线为腹襟翼板面没有缝的计算结果；图 28和图 29为8°迎角时相对起降构型，随板面实度的改变当地迎角的升力系数增量和俯仰力矩系数增量的变化曲线，其中横坐标为腹襟翼半模板面缝的个数。如图 28和图 29所示，与其他模型相比，半模两个缝升力系数增量最大为0.0579；未开缝模型俯仰力矩系数增量最大为0.01619，因此，本文飞翼布局无人机腹襟翼的最优板面实度为半模开两个缝或者不开缝。
|
图 26 升力系数增量的变化曲线 Fig.26 Curve of CL increment |
|
图 27 俯仰力矩系数增量的变化曲线 Fig.27 Curve of Cm increment |
|
图 28 升力系数增量随腹襟翼板面实度变化的关系曲线 Fig.28 Curve of CL increment with the solidities |
|
图 29 俯仰力矩系数增量随腹襟翼板面实度变化的关系曲线 Fig.29 Curve of Cm increment with the solidities |
通过上述腹襟翼的气动优化设计研究，得到了本文飞翼布局无人机最优的腹襟翼参数组合，即腹襟翼半模板面开两个缝或者不开缝，偏角为90°，且安装位置为无人机重心后25%MAC处。图 30为腹襟翼最优的设计结果(板面不开缝)相对起降构型升力系数增量和俯仰力矩系数增量的变化曲线，在起飞离地迎角(8°)处，最优结果比起降构型升力系数大0.057，俯仰力矩系数大0.0162。
|
图 30 最优结果相对起降构型气动力特性增量的变化曲线 Fig.30 Curve of aerodynamic characteristics increment of the optimal result versus the takeoff and landing configuration |
本文采用基于有限体积法的三维定常RANS方程求解器，研究了腹襟翼对飞翼布局无人机起降气动性能的影响，并对腹襟翼的关键参数进行了详细地优化设计研究，有以下结论：
1) 通过对比研究飞翼布局无人机干净构型、起降构型和腹襟翼打开的起降构型三种模型的低速气动特性得到，打开腹襟翼能够提高飞机起飞离地迎角的升力和俯仰力矩，有利于飞机的起降；
2) 对于本文研究的飞翼布局无人机，腹襟翼安装在重心后25%MAC位置处，偏角为90°，且板面实度较大时对飞机的起飞性能最好。在起飞离地迎角(8°)处，最优结果比起降构型升力系数大0.057，俯仰力矩系数大0.0162。
| [1] |
李永红, 黄勇, 王义庆. 翼身相对厚度对小展弦比飞翼布局跨声速气动特性及流动机理的影响研究[J]. 空气动力学学报, 2015, 33(3): 302-306. LI Y H, HUANG Y, WANG Y Q. Wing-body thickness effects on aerodynamic and vortex flow characteristics of common low-aspect-ratio flying-wing configuration at transonic flow[J]. Acta Aerodynamica Sinica, 2015, 33(3): 302-306. DOI:10.7638/kqdlxxb-2015.0045 (in Chinese) |
| [2] |
沈林成, 朱华勇, 牛轶峰. 从X-47B看美国无人作战飞机发展[J]. 国防科技, 2013, 34(5): 28-36. SHEN L C, ZHU H Y, NIU Y F. A survey of unmanned combat aircraft system development from X-47B[J]. National Defense Science & Technology, 2013, 34(5): 28-36. DOI:10.3969/j.issn.1671-4547.2013.05.006 (in Chinese) |
| [3] |
胡添元.飞翼布局飞机总体总体多学科设计优化研究[D].南京航空航天大学. HU T Y. Multidisciplinary design optimization for flying wing aircraft preliminary design[D]. Nanjing University of Aeronautics and Astronautics. (in Chinese) |
| [4] |
Daniel Dougherty, Adam Vore, Hamid Hefazi. High lift flap system for a tailless transport aircraft[R]. AIAA 2013-0965.
|
| [5] |
王彦雄, 祝小平, 周洲. 穿越微下冲气流的飞翼布局无人机控制方法[J]. 航空学报, 2015, 36(5): 1673-1683. WANG Y X, ZHU X P, ZHOU Z. A control method of flying wing UAV for penetration of microburst[J]. Acta Aeronautica et Astronautica Sinica, 2015, 36(5): 1673-1683. (in Chinese) |
| [6] |
单继祥, 黄勇, 苏继川. 小展弦比飞翼布局新型嵌入面航向控制特性研究[J]. 空气动力学学报, 2015, 33(3): 296-301. SHAN J X, HUANG Y, SU J C. Effect of the novel embedded control surfaces on direction control characteristic of low-aspect-ratio flying-wing configuration[J]. Acta Aerodynamica Sinica, 2015, 33(3): 296-301. DOI:10.7638/kqdlxxb-2015.0048 (in Chinese) |
| [7] |
吴军飞, 秦永明, 黄湛. 小展弦比飞翼标模纵航向气动特性低速实验研究[J]. 空气动力学学报, 2016, 34(1): 125-130. WU J F, QIN Y M, HUANG Z. Low speed experiment on longitudinal and lateral aerodynamic characteristics of the low aspect ratio flying wing calibration model[J]. Acta Aerodynamica Sinica, 2016, 34(1): 125-130. (in Chinese) |
| [8] |
孔轶男, 何开锋, 王立新, 等. 增加喷流控制的飞翼飞机机动仿真[J]. 空气动力学学报, 2013, 31(2): 186-191. KONG Y N, HE K F, WANG L X. Simulation of flying wing aircraft maneuver with leading edge jet control[J]. Acta Aerodynamica Sinica, 2013, 31(3): 186-191. (in Chinese) |
| [9] |
刘尚民, 孙健, 刘朝君. 大展弦比飞翼布局模型验证机飞行控制技术研究[J]. 飞行力学, 2013, 31(6): 558-560. LIU S M, SUN J, LIU Z J. Research on flight control technology of the high-aspect-ratio flying wing model aircraft[J]. Flight dynamics, 2013, 31(6): 558-560. DOI:10.3969/j.issn.1002-0853.2013.06.018 (in Chinese) |
| [10] |
马超, 王立新. 飞翼布局作战飞机起降特性分析[J]. 北京航空航天大学学报, 2009, 35(4): 429-433. MA C, WANG L X. Take-off and landing features of flying-wing configuration fighter[J]. Journal of Beijing University of Aeronautics and Astronautics, 2009, 35(4): 429-433. (in Chinese) |
| [11] |
康辉. 提高飞翼飞行器起降性能的腹部襟翼[J]. 国际航空杂志, 2008, 2: 47. KANG H. Belly flap innovation for blended-wing body[J]. International Aviation, 2008, 2: 47. (in Chinese) |
| [12] |
JAMES CHUNG, ERIC HALLBERG. Landing pitch control analysis for a blended wing body UCAV[R]. AIAA 2010-1035.
|
| [13] |
YANN D STAELENS, RON F BLACKWELDER. Novel pitch control effectors for a blended wing body airplane in takeoff and landing configuration[R]. AIAA 2007-68.
|
| [14] |
DANIEL DOUGHERTY, ADAM VORE, HAMID HEFAZI. High lift flap system for a tailless transport aircraft[R]. AIAA 2013-0965.
|
| [15] |
LU Y, ZHANG S G, LI X Q. A hazard analysis-based approach to improve the landing safety of a BWB remotely piloted vehicle[J]. Chinese Journal of Aeronautics, 2012, 25: 846-853. DOI:10.1016/S1000-9361(11)60454-8 |
| [16] |
沈冬, 张彬乾. 翼型俯仰控制的腹部扰流板研究[J]. 机械科学与技术, 2012, 31(2): 224-228. SHEN D, ZHANG B Q. Computational investigation on belly-flap for aerofoil[J]. Mechanical Science and Technology for Aerospace Engineering, 2012, 31(2): 224-228. (in Chinese) |