doi: 10.7638/kqdlxxb-2021.0282

基金项目: 翼型、叶栅国家重点实验室基金（6142201200106）详细信息

A multi-disciplinary global/local optimization method for flying-wing airfoils design

HTML全文

图 1 B-2轰炸机Figure 1. B-2 boomer

下载: 全尺寸图片 幻灯片20世纪90年代，波音公司提出了450座级的飞翼运输机方案—BWB-450飞翼运输机（图2），且在2004年与NASA合作设计了X-48系列的飞翼验证机（图3），以进行飞翼布局稳定性和操纵性研究[2-3]；同期欧盟也持续开展了多个飞翼布局运输机项目，包括MOB、VELA等[4-7]，2003年，英国剑桥大学联合美国麻省理工等团队，进行了静音客机项目（SAX）的研究，其中典型代表就是SAX-40[8-18]，如图4所示；俄罗斯重要空气研究中心（TsAGI）也提出了一种兼具飞翼气动效率和常规客机机身结构优势的750座级的混合型飞翼客机方案[19]；在国内，西北工业大学、北京航空航天大学、南京航空航天大学等也持续开展了飞翼客机的总体方案设计研究[20-24]。

图 2 波音BWB-450座飞翼运输机Figure 2. Boeing BWB-450 flying-wing transport aircraft

下载: 全尺寸图片 幻灯片

图 3 NASA X-48B验证机Figure 3. NASA X-48B testing aircraft

下载: 全尺寸图片 幻灯片

图 4 SAX-40 静音客机Figure 4. SAX-40 silent passenger plane下载: 全尺寸图片 幻灯片目前，无人侦查机、攻击机是一个新的发展潮流。美国已完成了多款飞翼布局无人攻击机的研制，如X-47B（图5）[25]、X-45C[26-27]、RQ180等。同时我国也成功试飞了“利剑”无人机。此外，俄罗斯的“电鳐”、法国主导的“神经元”、英国的“雷神”等无一例外地采取了飞翼布局。可见，飞翼布局无人机是未来无人机发展的一个重要方向。

图 5 X-47B无人机Figure 5. X-47B drone下载: 全尺寸图片 幻灯片当然无尾飞翼布局也存在着一些不足：1）飞翼布局存在先天性的纵向和航向的稳定性不足，操纵面设置复杂。同时其纵向尺度相对较短，使得提高纵向操作面效能困难较大，需要依赖多个操纵面和推力矢量等共同产生飞翼布局所需要的各种力和力矩，造成飞控系统操纵律设计困难。2）大型飞翼布局的巡航升阻特性不理想。由于飞翼布局的装载需求几乎完全依赖内翼区域，导致该区域的剖面厚度相对较大，限制了飞翼布局巡航速度的提高。当飞翼布局巡航速度超过其临界马赫数时，全机的激波阻力会显著提高，造成升阻特性损失严重。因此，提高飞翼布局飞机的临界马赫数尤为重要。飞翼翼型是提升飞翼布局性能的关键手段之一，为此关于飞翼布局翼型的研究层出不穷。根据设计模型的不同可将其分成三类：经典翼型设计、三维布局环境翼型剖面设计及二维/三维翼型综合设计。1）经典翼型设计。经典翼型是以机翼展向剖面为基础的二维问题。为此翼型设计首先需要以设计机翼飞行状态为基础，应用经验公式将三维设计状态映射为二维设计状态进行设计。Liebeck等[3]通过给出亚声速翼型的理想速度分布，然后进行了BWB-450内翼段翼型设计，得到了后缘反弯的LW109A的翼型；张彬乾等[28]进行了飞翼布局内翼段、外翼段翼型气动、隐身多学科设计，并分析了气动与隐身特性的关系；李权等[29]进行了飞翼布局超临界翼型的设计研究；黄江涛等[30]利用PCA目标降维方法进行了飞翼翼型多学科优化设计。经典翼型优化设计可以快速得到具有良好特性的翼型，但由于设计过程中并未考虑展向横流效应，设计结果在三维布局上难以得到理想的效果。此外，由于需要考虑气动、隐身和控制的多学科翼型设计对设计空间，包括维度和尺度均提出了更高的要求，大大增加了翼型设计难度，同时也降低了翼型设计结果对机翼性能的提升效果。2）三维布局环境翼型剖面设计。即直接在三维布局上进行翼型剖面优化设计。Pambagjo等[31]以NASA超临界翼型为基础，采用基于小扰动速势方程的余量修正法进行飞翼布局内翼段不同站位翼型剖面设计，设计翼型具有典型的超临界特征；Peigin等[32]以MOB项目的布局为基础，利用基于代理模型全局优化方法，开展了三维环境下翼型剖面外形的单点和多点设计；Hileman等[33]结合内翼段和外翼段翼型的设计要求，利用MSES软件开展了静音客机翼型的反设计；李沛峰等[34]采用直接曲率法（direct iterative surface curvature, DISC），在大型飞翼布局客机上进行外翼段翼型剖面外形设计，以削弱外翼段的激波阻力；Kuntawala等[35]利用基于样条曲线的FFD方法及伴随梯度方法，在飞翼布局上进行了翼型多剖面的气动设计；Zingg等[36]利用梯度算法进行了飞翼布局剖面外形和扭转角设计研究；周琳等[37]利用伴随方法建立了气动外形与隐身协同优化设计方法，并以类X47-B布局为例，对展向不同剖面翼型进行了气动、隐身协同设计；Kanazaki等[38]基于EGO方法，开展了小展弦比飞翼布局客机的翼型设计；Thomas等[39]对比了宽体型和窄体型飞翼布局客机翼型设计，并对比了二者的气动效率和燃油效率；Martins等[40]开展了考虑气动、配平和操纵性的飞翼布局客机设计，并对三个站位的翼型剖面进行优化设计。在飞翼布局三维环境进行翼型剖面设计，可以更好地依照当地流场条件，针对设计要求获得性能更优的翼型。但其最直接的问题就是三维状态物理模型计算量很大，涉及多学科优化问题，对计算资源的需求尤甚，更难实现在广阔设计空间的全局优化。梯度伴随方法的发展虽然可以大幅提高优化效率，但由于梯度算法先天的局部搜索特点，其设计结果对初始输入模型的性能依赖性较大，无法保证收敛到全局最优解。3）针对上述翼型设计存在的问题，Qin等[7]采用梯度优化方法，先利用经典翼型设计模型进行翼型设计，并以所得设计结果为初值，再进行飞翼布局三维环境下精细化剖面设计，实现了二维、三维条件结合的翼型设计，在一定程度上弥补了单纯经典翼型设计和三维环境剖面优化存在的问题。由此可见，飞翼布局翼型设计是结合了气动、隐身、控制等复杂多学科设计的问题。这些设计要求很难同时满足，给优化设计带来巨大的困难，进一步造成了经典翼型设计结果在三维布局上应用效果不理想。因此，本文基于飞翼布局展向气动特性分析，进行展向区域划分，分解翼型设计要求，构建分区经典翼型系列设计模型，并以之为基础，进一步在三维飞翼布局上进行多剖面翼型优化设计，提升机翼性能。1. 飞翼布局翼型设计翼型是飞行器的基础，翼型的特性对飞行器的特性有重要影响。常规布局翼型的主要目的是提高飞行器气动性能，而飞翼翼型设计则还需要考虑其对布局操稳的贡献。此外，对于考虑隐身影响的军用飞机而言，翼型对整机隐身特性影响尤为重要，在翼型设计中充分融入隐身影响设计，将对整机的RCS降低起重要作用。因此，飞翼布局翼型需考虑气动、控制、隐身等多学科综合设计要求。翼身融合的布局特性，造成飞翼翼型设计要求繁多且复杂，而飞翼布局展向不同区域的特性需求不一，因此对翼型的要求亦是迥异。因此，为了实现高效的飞翼翼型设计，本文基于飞翼布局研究和探索，以布局分区为基础，分解和简化多学科翼型设计要求，建立了分区翼型设计模型，并形成飞翼翼型系列。1.1 飞翼布局展向气动特性分析为了分析飞翼布局展向气动分布特点，本文选取典型双后掠飞翼布局（图6）为例。该布局外形类似X-47B飞机，内翼段后掠角接近60°，外翼段后掠角接近30°。由于飞翼布局内翼区承担了主要的装载任务，因此其当地厚度和弦长一般远大于外翼段。飞翼布局通常要求在保持纵向力矩配平的同时，具有高巡航升阻比以及良好的巡航效率。因此，本文选取典型跨声速巡航状态（[Math Processing Error]Ma=0.80,CL=0.30,Re=6.8×106）进行流场及气动特性分析。

图 6 飞翼布局展向分区示意图Figure 6. Spanwise partition of a flying-wing layout下载: 全尺寸图片 幻灯片图7展示了类X-47布局压力分布云图，图左侧是上表面压力分布，右侧是下表面压力分布。由图可知，在当前计算状态下，外翼段出现了更加明显的激波。图8和图9分别展示了相应的展向升力系数分布和环量分布。该布局内翼段当地升力系数较小，并沿着展长逐渐提高，在外翼区达到最大，随后逐渐降低。由于本文计算几何模型简化了翼尖区，导致翼尖区当地弦长过长，因此升力系数在翼尖区有个陡升。

图 7 布局压力分布云图Figure 7. Surface pressure on a flying wing下载: 全尺寸图片 幻灯片

图 8 展向升力系数分布Figure 8. Spanwise distribution of the lift coefficient下载: 全尺寸图片 幻灯片

图 9 展向环量分布Figure 9. Spanwise load distribution下载: 全尺寸图片 幻灯片由图8可知，飞翼布局内翼段虽然厚度较大，但由于升力系数较小，所以内翼段并未产生明显的激波，而外翼段高升力导致飞翼布局中间区域内侧开始产生强激波，横跨整个外翼区，直至延伸到翼尖区。而当地弦长的差距导致该布局的当地载荷从内翼段到外翼段逐渐降低（如图9所示），呈现近似三角形分布[7]，这种分布形态体现了激波阻力和诱导阻力权衡的结果。综上所述，对于双后掠飞翼无尾布局而言，中央体区和中间区承担着纵向力矩配平的任务，中间区至翼尖区则是气动减阻的重点关注对象，因此，本文将以双后掠飞翼布局，进行飞翼翼型设计模型和设计方法研究。1.2 飞翼布局翼型分区设计要求传统布局飞机，如运输机和无人机等，要求具有良好巡航气动效率和起降性能，因此其翼型一般要求如下：1）在飞机的整个使用范围内，翼型必须具有高巡航升阻比，良好的阻力发散特性以及较小的压缩性阻力增量等；2）具有足够的抖振边界，构成的机翼在超过飞机巡航马赫数一定范围内不出现抖振；3）后加载产生的低头力矩控制在一定范围内，以降低配平阻力和尾翼载荷；4）具有高的最大升力系数和失速特性，以简化增升装置的设计。与传统布局飞机相比，飞翼布局飞机对翼型气动设计进一步提出了更加严苛的要求，在传统的气动特性基础上，加强了对力矩特性、边界特性的要求等，具体来说有以下方面的设计要求：1）纵向自配平的力矩特性，由于尾翼的缺失和操纵面设置复杂等问题，飞翼翼型需要尽可能实现巡航设计点的力矩自配平，以缓解飞机操纵控制方面的压力，降低飞机控制系统的设计难度；2）更好的阻力发散特性，高跨声速乃至超声速巡航是未来飞机重要发展的趋势，而飞翼布局的装载要求导致其翼型厚度相对较大，因此需要提高翼型的阻力发散特性以满足飞机巡航效率、抖振边界等设计要求；3）与传统布局飞机不同，飞翼布局融合了机翼机身，因此内翼段机翼承担了装载要求，需要有大的容积和厚度；4）除气动设计要求外，良好特性的飞翼布局对其翼型还需其他学科要求的协同设计，如隐身、结构等。翼身高度融合的布局的空间外形、部件外形和外载荷分布等与传统布局不同，因此需要进行结构优化设计以满足飞行包线、飞行安全的设计要求等；而隐身特性直接关系到军用飞机的作战能力和战场生存能力，飞翼布局飞机作为天生的低RCS布局，更需要在翼型设计阶段融合RCS影响，以充分发挥布局优势。以上飞翼布局对翼型设计的整体要求，通常这些要求很难同时满足，而飞翼布局不同区域的性能需求不一，因此各处翼型的学科要求也不尽相同，需要结合三维布局分区的设计要求，提炼不同站位翼型的设计模型，分别进行优化设计。结合上述分析，本文基于经典的双后掠飞翼布局，建立了三维分区翼型设计理念，将布局分解为中央体区、中间区域和外翼区以及翼尖区，构建分区翼型设计模型，以满足整个飞机的设计要求。对于气动隐身特性需求较高的双后掠军用飞翼布局而言，其翼型通常需要满足以下要求：1） 中央体区，又称内翼区，该区域翼型在满足装载需求的同时，还需要具有良好的升阻特性和隐身特性。此外还需要一定抬头力矩，促进实现布局的纵向力矩配平；2） 中间区域和外翼区一般是激波集中产生的区域，因此该区域的翼型设计除保障装载需求以外，需要具有良好的巡航升阻特性以及阻力发散特性，且该区域距离重心位置较近，需要尽可能减少低头力矩乃至增加抬头力矩以满足纵向力矩的配平需求；3）翼尖区的翼型需要降低巡航阻力，提高升阻比，以获得高效的巡航效率，还要具备良好的低速升力特性和失速特性，推迟起降分离，为飞翼布局的起降阶段提供升力。此外，还需要良好的隐身特性。飞翼布局外翼区和内翼区翼型设计是两种截然不同的思路。气动设计方面，对于外翼段翼型，主要是减小巡航阻力同时提高隐身特性，以保证较高巡航效率和减小起飞着陆的失速速度，因此外翼段翼型大多采用超临界翼型；对于内翼段翼型，主要要求其具有较高的抬头力矩，以中和外翼段产生的低头力矩，实现巡航状态纵向力矩配平。1.3 飞翼布局翼型的典型特征以往飞翼翼型设计已经体现了分区的思路，主要是基于翼型力矩配平的设计策略。力矩系数和翼型的弯度线形状息息相关，为了满足纵向力矩配平的要求，BWB-450采用翼型后缘反弯[41]，即后缘反加载的策略设计内翼段翼型，以提供抬头力矩，而外翼段通过合理的翼型选配设计以及扭转分布，保证了配平状态下展向近似椭圆形状的载荷分布，减小了飞机的诱导阻力，使飞翼布局的气动特性优势得到了更加充分的发挥，如图10所示。但后缘反弯剖面的中央体设计，会导致巡航性能的牺牲，且需要更大的控制面和舵机实现飞机控制。

图 10 BWB-450飞翼布局运输机翼根翼型Figure 10. The entral body of a BWB-450 flying wing transport aircraft下载: 全尺寸图片 幻灯片针对这一问题，静音客机SAX-29、SAX-40采用了前缘下表面前弯的外形[33]。绕此类外形气动压心前移，因此这种翼型组成的中央体外形可实现无平尾的纵向平衡，且没有由尾翼平衡引起的升力损失。在巡航过程中燃油消耗引起的重心变化用增大矢量推力角来平衡，使静安定裕度保持为6％～9.5％，而采用后缘反弯的BWB布局，该裕度只能达到5％。图11展示了静音客机系列的翼型配置和压力分布，从中可以看出，内翼段采用了前加载翼型，当地升力系数较小，且主要集中于前缘部分，提供了抬头力矩，后缘升力系数基本为0。外翼段采用了典型的超临界翼型，当地升力系数较高，上表面存在弱激波，压力恢复较为和缓，且具有明显的后加载特征。此外，前加载翼型能够对气动特性和隐身特性进行合理折中，是一种理想的飞翼布局中央体翼型。如美国研制的B-2隐身轰炸机（图12）采用的就是前加载翼型。

图 11 静音客机翼型配置Figure 11. Airfoil configurations of silent aircrafts

下载: 全尺寸图片 幻灯片

图 12 B-2轰炸机中央体外形Figure 12. Central body of a B-2 boomer下载: 全尺寸图片 幻灯片综上所述，由于飞翼布局不同区域对翼型的设计需求不同，导致当地翼型的外形特征和压力分布不尽相同，因此飞翼翼型设计不可一概而论，应当结合不同区域的需求，建立分区设计模型，实现高效、便捷的翼型设计。2. 飞翼布局翼型设计方法2.1 传统修型设计早期的飞翼翼型设计主要依赖设计者的工程经验，结合设计对象的特性要求，对现有翼型外形进行修型。如波音公司的Liebeck通过给出翼型理想速度分布进行反弯翼型的修型设计[3]。2.2 反设计翼型反设计思想的提出可追溯到20世纪30年代NACA6系列翼型的设计。Jacobs提出翼型流场分析理论可逆向应用于翼型设计，即根据目标气动特性求解翼型几何[42-44]。常用的反设计方法包括余量修正法（Iterative Residual-correction Method）[15]、表面曲率法（Direct Iterative Surface Curvature，DISC）[16]、表面流线法（Constrained Direct Iterative Surface Curvature，CDISC）[20]以及Gappy POD方法[21]等。这些方法在飞翼翼型设计中得到广泛应用，虽然反设计方法效率很高，但需要预先给定压力分布，这一点依赖于设计人员的经验，而对于新的设计构型，压力分布很难给出，也不能保证存在对应的几何外形。其次，即使存在对应外形，也不能保证所给的目标流动分布就是最优的，这在很大程度上影响了设计的效果，反设计方法很难兼顾低速等其他设计要求[45]。2.3 基于数值优化设计方法基于数值优化的设计方法，主要设计思路是采用基于物理计算模型和现代数值优化算法结合的气动力设计方法对飞行器外形进行优化，其由参数化方法、优化算法以及物理计算模型构成。而根据优化算法的不同可将其分为两类：1）基于进化类算法的全局优化设计[46-47]；2）基于梯度算法的局部优化设计[48]。对于精细化优化设计问题，物理计算模型需要较大的内存和CPU资源，且进化类算法对物理计算模型的调用次数过多，复杂问题对计算资源的需求往往大到难以接受的地步，因此设计者们基于取样方法和代理模型技术，发展了基于代理模型的全局优化设计方法。优化设计方法也是应用最广泛的飞翼翼型设计方法。飞翼布局翼型设计是结合了气动、隐身和控制的复杂多学科设计问题。复杂的设计要求导致其与传统翼型设计迥异，需要在广阔的设计空间内进行全局优化，这对基于代理模型的优化设计造成了困难。针对代理模型优化的设计空间构造问题，本文采用了课题组开发的Bspline曲线的自适应方法[49]。该方法选择敏感区域加入新的变量，实现设计空间维度扩展，利用Bspline曲线的节点插入算法，对样本进行精确地高维重构。同时利用自适应设计空间扩展方法[50]对表现较好的有效样本的设计变量边界进行检测，进行尺度扩展，有效提高了样本的利用率，确保优化在充足的设计空间内高效地进行。结合分区翼型设计模型共同组成了二维翼型设计方法。针对经典翼型设计结果对飞翼布局性能提升不理想的问题，本文建立了三维环境下翼型修型设计。由于三维流动环境下翼型设计的设计变量数和计算代价较大，基于代理模型的全局优化方法在面临高维设计变量时会面临“维度灾难”问题，因此课题组发展了基于样条曲线的FFD（Free Form Deformation，FFD）方法以及梯度伴随的序列最小二乘算法（Sequential Least Programming，SLSQP）[51]的梯度优化设计方法，实现高效的三维环境翼型修型设计。FFD方法适用于任意网格形式和拓扑结构，可以简单、直接且高效地控制物体的变形，因此在飞行器三维气动外形优化设计中应用十分广泛。序列最小二乘算法是一种适用于非线性约束的梯度优化算法，可以处理等式和不等式约束。该算法可以逐次逼近目标函数的二阶近似方程和约束的一阶近似方程，能够有效处理拥有大规模设计变量的非线性多约束问题。2.4 气动隐身多学科设计为了充分发挥飞翼布局的气动隐身特性，需要对飞翼翼型进行多学科设计。气动隐身多学科设计包括气动、隐身计算模型以及设计方法。高效高精度的物理计算模型是优化设计的基础。在气动计算中，本文二维翼型/三维布局均采用高精度的基于雷诺平均Navier-Stokes（RANS）方程，实现高效的流场模拟。无黏项采用迎风格式，黏性项采用中心离散格式，采用近似因式分解法进行时间推进求解，湍流模型一律采用SST两方程模型。在隐身计算中，过去大多是将二维翼型沿展向拉成三维机翼，利用高频近似物理光学法（Physical Optics，PO）及其混合方法进行RCS计算，计算效率低。因此本文在二维翼型设计中，采用高效的二维矩量法[52]进行翼型隐身特性的计算。该方法是麦克斯韦方程的精确求解方法。文献[49]验证了二维矩量法的精度和可靠性。在三维布局优化设计中，采用物理光学法进行三维RCS特性求解。气动、隐身多学科设计是典型的多目标问题。对于多目标问题，优化解表现为一个解集，即Pareto前沿。设计人员根据Pareto前沿的分布进行多目标决策，获得最终的优化解。多目标遗传算法是求解多目标优化问题Pareto前沿的最有效方法之一。遗传算法（GA）借鉴生物进化的思想，通过计算机模拟物种繁衍父代遗传基因的重新组合与“优胜劣汰”自然选择机制的联合作用，解决科学与工程中的复杂问题。本文采用NSGA-Ⅲ（Non-dominated Sorting GeneticAlgorithm Ⅲ）[53]方法进行翼型多学科优化设计。2.5 翼型优化设计框架基于三维分区设计理念，构建了不同站位的二维设计模型，结合自适应参数化方法、自适应设计空间扩展方法、Kriging代理模型和NSGA-III算法，建立了二维翼型设计方法。并针对三维效应下二维设计结果难以充分发挥作用的问题，结合FFD方法和序列最小二乘法建立三维环境翼型修型设计方法。结合气动、隐身计算模型，建立了“全局+局部”的飞翼翼型多学科优化设计方法。优化流程如图13所示。

图 13 “全局+局部”飞翼翼型多学科优化设计方法流程图Figure 13. Flow chart of the global/local multidisciplinary optimization method for designing flying-wing airfoils下载: 全尺寸图片 幻灯片3. 飞翼布局翼型设计方法3.1 气动隐身多学科设计根据飞翼布局飞机翼型设计要求，以NACA 65016翼型[37]为基础翼型进行优化设计，设计状态为[Math Processing Error]Ma=[Math Processing Error]0.70、CL=0.25、Re=20×106。翼型设计模型如下：

图14和图15为设计翼型外形和压力分布对比，从中可以看出，比起初始翼型NACA 65016，设计翼型cm0.03前缘半径减小，后缘负弯度增加，具有明显的前缘正加载和后缘反加载特征，上表面外形趋于平坦，最大厚度位置后移。设计状态下初始翼型的强激波被削弱为两道弱激波，有利于提高阻力发散特性，如图16所示，阻力发散马赫数提高至0.74，同时翼型力矩特性明显提升，如图17所示。

图 14 翼型外形对比Figure 14. Comparison of airfoil shapes下载: 全尺寸图片 幻灯片

图 15 压力分布对比Figure 15. Comparison of pressure coefficients下载: 全尺寸图片 幻灯片

图 16 阻力发散特性对比Figure 16. Comparison of drag divergence下载: 全尺寸图片 幻灯片

图 17 力矩特性对比Figure 17. Comparison of momentum coefficients

下载: 全尺寸图片 幻灯片为了验证设计翼型的优化效果，将设计翼型的装配到类X47-B三维布局[37]。X47-B是典型的小展弦比无人机布局，要求有高效的巡航效率和隐身特性。本节首先研究设计翼型对其气动特性的改善。该布局几何外形和基本参数分别如图18和表1所示。其中，设计翼型cm0.03装配到内翼段，利用Xfoil改变设计翼型的相对厚度为13%和11%，其中13%厚度的cm0.03_13装配到中间区域，11%厚度的cm0.03_11装配到翼尖区域，形成设计方案OptModel_1，然后与对称翼型配置的布局方案BaseModel进行计算对比。计算状态如表2所示。

图 18 基准布局外形Figure 18. Shape of the base model

下载: 全尺寸图片 幻灯片表 1 基准布局几何参数Table 1. Geometric parameters of the base model

下载: 导出CSV | 显示表格表 2 气动计算状态Table 2. Conditions for aerodynamic computation

下载: 导出CSV | 显示表格图19为计算结果和压力分布云图对比，相较于BaseModel，OptModel_1阻力系数降低41.9 counts （1 counts = 0.0001），力矩系数提高了0.013。内翼段激波基本消除，外翼段前缘低压区向前扩张，后缘压力线分布更加均匀，压力恢复和缓。虽然由于布局参数和重心位置等原因，OptModel_1还未达到完全配平，但OptModel_1低头力矩有很明显的减弱，纵向力矩配平明显改善。但OptModel_1外翼段仍旧有明显的激波，且前缘低压区较小，升阻特性没有充分达到优化设计的效果，还有很大的改进余地。

图 19 三维布局气动计算结果对比 Figure 19. Comparison of surface pressure下载: 全尺寸图片 幻灯片

综上，对飞翼布局翼型设计而言，不同区域的翼型设计不可一概而言，应当结合当地区域的设计要求，进行分区翼型设计。3.2 分区翼型系列设计根据上述讨论，建立分区翼型设计模型，并形成基于区域划分的翼型系列。中央体区翼型设计状态为[Math Processing Error]Ma=0.65、CL=0.25。气动特性要求在保持抬头力矩同时提高升阻特性。

中间区翼型设计状态为典型跨声速设计状态[Math Processing Error]Ma=0.73、CL=0.55。基准翼型为相对厚度13%的对称翼型NACA 65013。Kink 站位翼型气动设计要求在保持抬头力矩同时提高升阻特性和阻力发散特性，设计模型如式（3）所示。

翼尖区翼型的设计要求为减阻和提升阻力发散特性，约束为设计翼型的力矩系数不小于初始翼型，同时保持最大厚度不减小。设计状态为[Math Processing Error]Ma=0.73、[Math Processing Error]CL=0.35。隐身设计要求与之前保持一致。基准翼型为NACA 65011。设计模型如式（4）所示。

表3为设计翼型的计算结果对比。相较于各初始翼型，设计翼型阻力特性均有明显的改善，力矩系数满足约束。图20为不同区域设计翼型外形和压力分布对比。可以看出，设计翼型前缘半径减小，呈现明显的前加载特征。由于抬头力矩约束，内翼段和中间区域翼型后缘均出现不同程度反加载特征。翼尖区域翼型呈现典型跨声速压力分布特征，跨声速设计状态的翼型最大厚度位置前移。在亚声速低载荷状态下，内翼段翼型设计前后的流场并未出现明显的激波。由于不同设计升力系数，中间区域设计翼型和翼尖区域翼型呈现了不同的跨声速特征。表 3 设计翼型计算结果对比Table 3. Comparison of computational results

下载: 导出CSV | 显示表格

图 20 不同区域设计翼型外形和压力分布对比Figure 20. Comparison of airfoil shapes and pressure distribution on different partitions

下载: 全尺寸图片 幻灯片从图21可以看出，装配分区设计翼型的布局阻力降低了79.3 counts，力矩系数提高0.02，改善明显。内翼段激波基本完全消除，外翼段前缘低压区扩张，后缘压力线分布更加均匀，压力恢复和缓。但对飞翼布局而言，机翼横向流动效应明显，对此在翼型设计模型中并未加以考虑，因此未能充分发挥翼型设计的优势，布局的激波依旧比较明显。鉴于此，本文提出了三维环境的翼型修型设计方法，进一步提升翼型性能。

图 21 三维布局气动计算结果对比Figure 21. Comparison of computational results for a 3D layout

下载: 全尺寸图片 幻灯片3.3 三维环境多剖面翼型设计通过二维翼型全局设计，三维布局的气动特性明显提升.为了进一步提升设计翼型性能。使之充分发挥设计效果，在二维设计完成后，利用基于气动伴随方程的梯度方法，进行三维环境翼型修型设计，即优化的初始外形为OptModel_2。设计目标为巡航状态下减阻和提高RCS特性，同时保持三个站位处翼型最大厚度不减小以及保持低头力矩不增加。计算模型如式（5）所示。

最终优化结果为OptModel_3。从图22可以看出，相较于OptModel_2，OptModel_3阻力系数减小35 counts，力矩系数保持不变，并且外翼段的激波明显减弱，翼尖区的激波被基本消除。图23为各站位翼型剖面和压力分布对比，红色为三维修型翼型和压力分布，黑色为初始翼型和压力分布。从图23中可以看出，内翼段的修型翼型没有明显变化；中间区的修型翼型前缘半径减小，中部正弯度增加，后缘负弯度增加；翼尖区的修型翼型前缘正弯度和后缘负弯度明显增加。

图 22 三维布局气动计算结果对比Figure 22. Comparison of computational results for a 3D layout下载: 全尺寸图片 幻灯片

图 23 三维优化设计翼型和压力分布对比Figure 23. Comparison of airfoil shapes and pressure distributions obtained by 3D optimization design

下载: 全尺寸图片 幻灯片由上文可知，后缘反加载的翼型特征对气动特性有不利影响，而前缘正加载翼型可以很好地协调翼型的气动特性。在二维翼型设计中，各站位设计翼型均呈现明显的前缘正加载和后缘反加载，共同配合以满足气动特性设计需求。而三维修型设计对二维设计结果的外形进行了微调，尤其是前缘正加载和后缘反加载的分布，即翼型弯度分布。通过调整翼型各剖面的弯度分布，改善了布局的载荷分布，使其气动特性得到进一步提升。3.4 气动隐身多学科设计由上文可知，军用飞机翼型除了良好的巡航升阻特性和纵向力矩配平特性以外，还需要良好的隐身特性。因此，本节利用多目标遗传算法进行各区翼型气动隐身设计。在此过程中，各区翼型气动设计要求不变，隐身特性要求在入射频率9 GHz下，前向[Math Processing Error]±30°均值降低。入射角示意图如图24所示。

图 24 入射角度示意图Figure 24. Sketch of the incident angle

下载: 全尺寸图片 幻灯片为了在不损失布局气动特性的基础上，进一步提升飞翼布局外翼段的RCS特性，翼尖区设计两组翼型，分别为气动特性占优的翼型和隐身特性占优的翼型。气动特性占优的翼型装配于三维布局的外翼段距离翼尖1/4处位置，隐身特性占优的翼型位于外翼段翼尖处。 表4为三维布局隐身计算状态。表 4 隐身计算状态Table 4. Conditions for the stealth computation

下载: 导出CSV | 显示表格随后利用气动隐身伴随方法进行进一步的三维环境翼型修型设计，设计目标为减阻和提高RCS特性，保持低头力矩系数不增加及各剖面翼型厚度不降低。设计模型如式（6）所示。

图25～图28为Aero/RCS_OptModel_Ⅰ与设计结果Aero/RCS_OptModel_Ⅱ气动隐身特性计算对比。Aero/RCS_OptModel_Ⅱ阻力系数降低35 counts，力矩特性保持。从压力分布云图对比可以看出，Aero/RCS_OptModel_Ⅱ中央体区、中间区域和翼尖区的激波基本消除，但外翼段中部还有较为明显的激波。在各区域设计翼型对比中，翼尖处翼型的外形改变较为明显，翼型前缘正弯度和后缘负弯度有明显的增加，其余剖面翼型外形变化不大。表5为前向RCS均值对比，可以看出 Aero/RCS_OptModel_Ⅱ的隐身特性与Aero/RCS_OptModel_Ⅰ基本持平，改善不大。

图 25 三维布局气动隐身计算结果对比Figure 25. Comparison of computational results for a 3D layout

下载: 全尺寸图片 幻灯片

图 28 俯仰RCS特性对比Figure 28. Comparison of pitch RCS

下载: 全尺寸图片 幻灯片

图 26 三维优化设计翼型和压力分布对比Figure 26. Comparisons of airfoil shapes and pressure distributions obtained by 3D optimization design

下载: 全尺寸图片 幻灯片

图 27 偏航RCS特性对比Figure 27. Comparison of yaw RCS下载: 全尺寸图片 幻灯片综上所述，在分区系列设计中，翼型气动隐身设计对三维飞翼布局的气动隐身特性有明显的改善。但由于二维翼型模型无法考虑横向流动效应，导致设计翼型装配到三维布局上无法充分发挥设计效果。三维环境的翼型修型设计可以进一步提升布局的气动特性。对于隐身特性而言，翼型分区设计效果显著，而三维环境的修型设计对隐身特性的改善不明显。表 5 前向RCS均值（单位：m2）Table 5. Average frontal RCS (unit: m2)

下载: 导出CSV | 显示表格4. 结论与展望为了充分发挥飞翼布局飞机气动、隐身特性的布局优势，需要对其翼型进行高效设计。本文针对飞翼布局翼型的气动、隐身设计问题，提出了一种“全局+局部”的翼型多学科优化设计方法，以类X47-B布局为例，进行了内翼区、中间区和翼尖区的翼型优化设计，得到结论如下：1） 飞翼布局各区翼型设计要求不尽相同。内翼区翼型需要抬头力矩以帮助布局达到纵向力矩配平，以及良好的隐身特性，中间区翼型要求具有良好的巡航升阻特性、阻力发散特性和一定的抬头力矩以配合实现纵向力矩配平，而翼尖区则面临气动减阻和隐身特性的设计要求。2） 由于抬头力矩约束和隐身设计要求，飞翼布局翼型均呈现明显的前缘正加载和后缘反加载的外形特征。翼型两种加载特征的合理配置，可以在保持纵向力矩配平的同时，实现飞翼布局气动、隐身特性的共赢。3） 基于飞翼布局各区翼型设计要求，建立了分区翼型设计模型，利用基于代理模型的全局优化方法进行优化设计，并形成了分区翼型系列。为了进一步提高翼型特性，在分区翼型系列的基础上，利用梯度伴随方法进一步进行三维环境翼型修型设计。建立的“全局+局部”的翼型多学科优化设计方法，可以提高优化设计效率，使设计结果更加满足工程实际的设计要求。4） 飞翼翼型因飞翼布局飞机的种类、需求导致其设计要求不可一概而论。本文主要针对飞翼布局无人军用机进行了翼型系列的设计，未来将在此基础上，开展考虑起降特性的大展弦比的飞翼运输机、考虑层流设计的高空长航时飞翼无人机等翼型。由于飞翼布局的三维流动效应很强，因此拟结合风洞试验，对设计翼型系列的气动性能进行验证，为飞翼翼型的选配提供进一步支撑。