基于改进复合形法的平流层飞艇外形优化设计.pdf

61第40 卷第6 期2023年6 月真机计仿算文章编号：10 0 6-9 348(2 0 2 3)0 6-0 0 6 1-0 5基于改进复合形法的平流层飞艇外形优化设计曹胜鸿1-2,杨燕初1.2,张天晓12,张航悦1.2（1.中国科学院空天信息创新研究院，北京10 0 0 9 4；2.中国科学院大学，北京10 0 19 0）摘要：针对传统飞艇外形设计过程中仅将阻力作为优化目标的问题，先参照多学科优化设计理念，在原有阻力模型的基础上，引进体积与表面积对飞艇外形优化设计的影响，建立飞艇优化的数学模型，构建优化框架，再结合最优化设计中的复合形法，依据双向并行探索理念，提出一种改进复合形法，并将其运用到飞艇外形优化中，结合优化值给出不同需求下的优化设计方案，最后结合实际算例仿真分析得出改进复合形法对于飞艇外形的优化是可行的结论。关键词：平流层飞艇；多学科优化；气动阻力；体面比；复合形法中图分类号：TP391.9文献标识码：BOptimization Design of Stratospheric Airship ShapeBased on Improved Composite MethodCAO Sheng-hong,YANG Yan-chu-2,ZHANG Tian-xiao-2,ZHANG Hang-yue1,2(1.Institute of Aerospace Information Innovation,Chinese Academy of Sciences,Beijing 100094,China;2.University of Chinese Academy of Sciences,Beijing 100190,China)ABSTRACT:In view of the problem of only considering resistance as the optimization objective in the traditional air-ship shape design process,this paper first refers to the concept of multidisciplinary optimization design,and introducesthe influence of volume and surface area on the shape optimization design of airship based on the original resistancemodel.Then,this paper establishes the target function of airship optimization,and constructs the airship shape optimi-zation process framework.In addition,this paper proposes an improved composite form method by combining the com-posite form method and the concept of bidirectional parallel exploration in optimization design,and it is applied to theoptimization of the airship.Moreover,this algorithm can provide targeted optimization design scheme according to dif-ferent needs.Finally,it is concluded that the optimization algorithm based on the improved composite form method isfeasible for the optimization of airship shape.KEYWORDS:Stratospheric airship;Multidisciplinary optimization;Aerodynamic drag;Decent ratio;Complex method1引言临近空间（NearSpace)通常是指2 0 10 0 km的空域，在地球大气中处于平流层和中层，其在环境监测、区域通信、实时监视、预警与导弹防御等领域有很大的应用价值。飞艇是一种轻于空气的航空器，具有垂直起降、长航时、载重大、低能耗、噪声及排放小等特点2 ，在临近空间的探测方面具有广泛的应用价值。平流层飞艇是临近空间飞行器的一种，近几十年来，随着航空技术的发展，采用新材料和新技术的大型飞艇得到飞速发展。收稿日期：2 0 2 1-0 9-15修回日期：2 0 2 1-11-0 9因此，有针对性地对平流层飞艇进行优化设计，即多方面考虑飞艇的设计因素就显得尤为重要。平流层飞艇总体设计与优化技术作为方案设计的关键，将会直接决定飞艇系统的整体性能和操控3。国内多家研究机构开展了大量的飞艇设计优化工作，例如，中航通飞研究所的秦何军针对飞艇的气动外形进行优化设计4，中国特种飞行研究所的张超针对重在飞艇低阻囊体进行外形优化5。但上述优化设计仅以阻力系数作为优化对象,缺少对飞艇结构因素的考虑。通过对比分析发现，在改变长细比时，体面比与长细比成正相关，阻力系数与长细比呈现负相关，因此，本文保留阻力、体面比两个飞艇设计因素，建立数学模型，结合实际需求分配两者权重，结合改进的复合形法，构建优化设计框架，通62过Matlab软件运算得出艇体长度与艇体最大截面直径，依照实际需求，提供三种合理的飞艇外形设计方案，并结合实际算例仿真分析，验证优化艇形和优化算法的可行性。2算法介绍2.1飞艇优化程序框架本文在进行飞艇优化设计过程中，结合最优化设计中的复合形法进行改进，进行目标函数的最优值求解，首先输入原始艇形数据，结合多学科优化设计理念，建立飞艇优化目标函数，并根据实际需求，合理分配各项函数模型的权重，利用Matlab软件结合改进复合形法探索优化值，依照实际需求，有侧重的针对算法所得飞艇长度、飞艇最大截面直径、飞艇长细比进行定向飞艇优化，并提供后期仿真分析的艇形数据，其飞艇优化流程框架如图1所示。需要指出的是，本文所提出的改进复合形算法也可与其它最优值求解方法相结合使用，推广至其它最优值求解问题。开始输人初始艇形相关参数构建飞艇优化目标函数F=Wfi+W2f2+w.J,分配目标函数权重设置相关约束条件基于改进复合形法的飞艇优化算法F是否满足收敛条件min(F)T截面最大直径长细比飞艇长度优化飞艇外形参数结束图1飞艇优化流程框架2.2改进复合形法在约束问题的最优化求解过程中，往往会出现导数和梯度较难求得的情况，尤其是多元函数，甚至会没有明确的工程函数6 。因此常常采用直接求解的方式，复合形法的基本求解思想是基于无约束最优化问题直接求解方式的单纯形法，仅需在每次扩张、收缩之后验证新得到的点是否在可行域内即可。单纯形法最早由Spendley、H e x t 和Himsworth提出7 ，后由Nelder和Mead改进成为变型的单纯形法8 ,其基本思想是：三点构造单纯形，沿最差点至剩余点所构中心点的方向寻找最差点的替代点，从而构造新的单纯形，其探索过程如图2 所示。X1X图2单纯形法探索过程取三个点并计算其函数值，分析三个点的函数值进行排列得到最差点X，次差点X。最好点X,，取次差点X。与最好点X,的形心X，连接最差点X，与形心点X。做第一次等步长的扩张得到反射点X，若反射点的函数值优于最差点Xh则再进一步扩张，否则，将步长减半，如果新得到的点函数值依旧差于最差点函数值，则直接改变探索方向，其中单纯形法的探索过程，还需要将新得到的点与次差点进行比较,从而考虑新得到点的删减。经典单纯形法需要计算剩余点所构中心点，并将新得到的点与最差点、次差点进行比较，且随着探索迭代次数的增加，在最优值附近的探索速度较慢，探索步骤缓慢、繁琐，极大地降低了约束问题的求解效率。单纯形法求解最优化工程函数问题的关键在于单纯形的构建，而在单纯形构建的过程中，选点将会直接影响探索方向的准确性，经典单纯性法求解时需要计算次差点与最好点的形心，并与最差点相连，构成新的探索方向，杨晋提出可以直接连接最差点与最好点构成新的探索方向9 ,李庆高在非线性约束问题同样提出直接在最好点与最差点的方向探索得到最差点的替代点10 本文在杨晋单纯形法的基础上进一步改进，提出双向探索的求解思想，即沿着最差点和最好点的方向探索最差点的替代点的同时，沿着次差点和最好点的方向探索次差点的替代点，在两个收敛方向同时求得第二次单纯形的构造点，精简单纯形构造的计算步骤，提高约束问题的求解效率，其探索步骤如图3所示。XXh2XXh3XhX92图3改进算法探索过程示意图63改进单纯形法的计算过程如下：1)取三个点并计算其函数值，分析三个点的函数值并进行排列，从而得到最差点X，次差点X。，最好点X，连接次差点X，与最好点X,做第一次等步长的扩张，同时连接最差点X,与最好点X，做第一次等步长的扩张，如式(1)所示(Xh=X,+(X,-X,)(1)(Xg,=X,+(X,-X,)：反射系数，单纯形法中通常取为1,复合形法中常取1.3;Xh：最差点X,相对于最好点X,的反射点；X：次差点X相对于最好点X,的反射点；2)分别比较两个反射点与最差点、次差点的函数值，如果反射点函数值小于最差点、次差点的函数值，则说明探索正确，进一步扩张,如式(2)所示(X4=X,+(X,-Xh)(2)(Xg4=X,+(X,-Xgr)：扩张系数，一般取为2；Xh4：反射点Xh相对于最好点X，的扩张点；Xg4:反射点X相对于最好点X,的扩张点；3)如果反射点函数值比最差点、次差点函数值要大，则减少扩张步长，即收缩。如式(3）所示(Xh2=X,+(X,-Xhr)(3)(X2=X,+(X,-Xg,):收缩系数，一般取为0.5；Xh2：反射点Xh相对于最好点X，的收缩点；X2：反射点X，相对于最好点X,的收缩点；4)如若收缩点的函数值仍旧大于最差点、次差点的函数值,则说明探索方向有误，改变探索方向，即缩小单纯形，如式(4)所示(X3=0.5(X,+X,)(4)(Xr3=0.5(X,+X,)Xh3：最差点X,相对于最好点X,的压缩点；Xg3：次差点X，相对于最好点X,的压缩点；5)探索所得点与第一次单纯形的最好点构成第二次探索单纯形的三个点，图3所示蓝色线即为第二次单纯形的可能情况。单纯形构造完成后，逐次迭代，直到最好点与形心点几乎为一个点，终止迭代的判据如式（5）所示1n+12(x()-f(X(2)J)28(5)（n+X）：最好点；Xt2:形心点。针对收敛性判据，除上述判断形心点与最好点是否几近重合之外，也可以判断最后一次所构单纯形是否几近为一点，其收敛性判断的主体思想都是一样的，即判断最终构造点几乎为一点。其改进复合形法的算法流程图如图4所示开始计算初始单纯形点的函数值/(x),j=1,2,3排序最好点X，次差点X2，最差点X3计算次差点、最差点相对于最好点的反射点所得反射点函数值是否F计算次差点、最差点相对小于对应的次差点、于反射点的收缩点最差点函数值TY计算次差点、最差点相对于反射点的扩张点所得收缩点函数值是否小于对应的次差点、T最差点函数值所得点即为次差点计算次差点、最差点相对最差点的替代点于最好点的压缩点F是否满足收敛性判据T结束图4改进算法程序框图上述单纯形法的改进，省略了经典单纯形法中对于形心点的计算，同时对杨晋算法9 单一方向的探索进行改进，在次差点与最好点、最差点与最好点两个方向进行计算和探索，以此求得新的单纯形探索构建点，省略了反射点与次差点的比较过程，算法简单，探索速度快。本文借用一个案例函数对改进算法和经典算法进行比较，验证改进算法的可行性。案例目标函数如式(6)所示,收敛性要求0.0 1。f=X+2X2-4X,-2X,X2(6)两种方法单次单纯形最优值与送代次数的关系对比图，如图(5)所示。经典算法改进算法的送代效率对比6050经典算法40一改进算法3020100-100510152025X-送代次数图5经典算法与改进算法单次最优值对比通过图5的对比分析，针对式（6)所示案例目标函数，改进算法在第3次迭代过程中，出现大幅度改变，从第6 次迭代开始进入函数收敛区域，而经典算法在第13次送代才有大幅度的转变，从第2 0 次迭代开始进人函数收敛区域。改进算法的迭代速度相较于经典算法提高了58.33%，由此可64(10）得，改进算法的探索速度高于经典单纯形法，极大地提高了函数最优值求解的效率。3飞艇优化目标函数以往的飞艇优化设计着重于阻力系数的优化，忽略其它设计因素，本文在进行多次设计优化分析后，针对飞艇初始优化阶段，保留阻力系数和体面比两个重要的设计因素。其中，阻力系数直接决定飞艇推进功率，进而影响能源供给；体面比是艇囊体积与面积的比值，决定着系统浮力的有效利用率1,其目标函数数学模型如式(7)所示FcD入F=W1+W2(7)FCDO入。W1,W2为阻力和体面比的权重，Fcp为阻力，Fcpo为初始阻力，入为体面比，入。为初始体面比值。3.1阻力模型Hoerner研究了流线型物体的阻力，对于多数飞艇,其艇体的气流是流12 ,因此阻力系数关系如式(8)所示CDV111=4(+6(）1.2+24(2.7(8)C其中1为艇体长度，d为艇体最大截面直径。Hoerner提出对于实际水准的表面光洁度和R。510，表皮摩擦系数C,如式(9)所示0.043C(9)R综合式（8）式（9），给出阻力系数公式，如式（10)所示1.230.172+0.252+1.032CDVR1其中雷诺数计算参见式（11）PaulR.=(11)其中为粘性系数，p。为气体密度，为风速。故阻力模型函数如式（12）所示FDDVP.2.Vi2=0.5(12)其中k，为整艇与裸艇之间的阻力系数比,由飞艇的实际构型决定，V。为囊体体积。3.2体面比模型针对飞艇初始化设计阶段，王全保2 仅考虑表面积带来的影响，忽略优化过程中囊艇体积的变化，故此，本文引人体积模型，建立体面比影响模型，体面比模型可参照式（13）所示，式（14）-（15）为优化前期体积和表面积估算公式，其精确计算公式参照式(16)-（17)入=(13)S。其中囊体体积估算公式参照式（14）0.4653(14)d对于飞艇优化初始阶段，艇体表面积估算公式，可参照式(15)S。=(2.14+(l-6 0)0.08-(d-15)0.028 2T(V.)(15)其中体积精确值的计算公式参照式（16）=T(16)一0对于飞艇表面积精确值的计算公式，可参照式（17）2S。=2 元（1+dx(17)dx4算例优化仿真分析本文选取英国物理实验室所提出的一种低阻力飞艇形状,即后期Skyship飞艇形状的基础,NPL低阻力系数飞艇形状12 ，作为飞艇前期优化的外形母线函数,其数学表达式如式(18)所示210 xa+a6222+510,以及实际飞艇设计案例,取初始飞艇长度6 0 m，最大截面直径2 0 m，长细比3。取飞艇高度为2 0 km高空，大气压力547 4.8 Pa，风速2 5m/s，空气密度取0.0 8 8 kg/m，根据实际设计经验，以一定的最佳长细比范围为约束，此外考虑到飞艇的最大浮力，故以固定的囊体体积为约束,，如式(19)所示652.54.1d(19)(v。=1117 2.6 m3运用改进复合形法进行目标函数的优化探索，借用Matlab软件生成优化后的飞艇外形，如图6 至图8 所示，其中图6 对应阻力权重为1,体面比权重为0,以此求得以阻力为优化重点时的优化艇形，图7 对应阻力权重为0，体面比权重为1，以此求得以体面比为优化重点时的优化艇形，图8 对应阻力权重为0.5,体面比权重0.5,以此求得阻力与体面比平衡时的优化艇形，需要指出的是，在实际设计优化工作过程中，设计者可以根据实际需求更改调配、本文只给出常用的几组权重分配比例，此外，也可以针对目标函数增添其它设计影响因素，本文保留两项随长细比变化呈现相反变化的阻力和体面比因素项。此外，根据所得优化值点，本文给出三个不同侧重点的优化结果，分别对应每组图中的除原始艇形外的三条优化曲线，若设计考虑目标函数值最优，则选择以优化点的最大直径为恒定值的优化结果，若设计考虑长细比的范围，则选择以优化点的长细比为恒定值的优化结果，若设计考虑飞艇的微小外形优化，则以优化点的艇长为恒定值的优化结果。权重分配为阻力为1，体面比为0 时的优化结果如图7所示。权重为1:0 时NPL阻力型飞艇优化外形初始艇形最大截面直径优化艇形一长细比优化艇形艇长优化艇形105026-10X-飞艇长度2034图7权重为1：0 时飞艇优化外形从图7 可以看出，当以阻力作为优化目标时，优化后的艇形会变细长，即整体艇体长度变长，最大截面直径变小，抗阻能力增强,其详细相关对比数值参见表1。表1权重为1：0 时优化艇形与原始艇形对比最大截面直长细比优艇长优初始艇形径优化艇形化艇形化艇形阻力(N)436.8340384.5486398.6215415.4314阻力优化率11.97%8.75%4.90%体面比3.53683.11063.34593.4544体面比优化率-12.05%-5.40%-2.33%从表1中可以看出，当权重分配为1：0 时，从三个侧重点优化所得艇形的阻力得到有效地改善，体面比则会相应的减小。权重分配为阻力为0，体面比为1时的优化结果如图8所示。权重为0:1时NPL阻力型飞艇优化外形10初始艇形一最大截面直径优化艇形长细比优化艇形一艇长优化艇形20-1020X-飞艇长度3图8权重为0：1时飞艇优化外形从图8 可以看出，在以体面比为优化目标时，优化后的艇形会变短胖，即整体艇体长度变短，最大截面直径变大，更有利于艇体浮力的利用，其详细相关对比数值参见表2。表2权重为0：1时优化艇形与原始艇形对比初始最大截面直长细比优艇长优艇形径优化艇形化艇形化艇形阻力（N)436.8340439.7354506.1164557.3620阻力优化率-0.66%-15.86%-27.59%体面比3.53683.54553.66623.7099体面比优化率0.25%3.66%4.89%从表2 可以看出，当权重因子全部分配给体面比时，所得优化艇形的体面比得到显著提升，相应的截面变大，艇长变短，阻力会增大。通过图7 与图8 的优化结果可以看出，体面比和阻力是优化过程中两个呈现相反增量的目标量，所以，在实际的优化过程中，还需要依照实际需求分配阻力与体面比的权重。本文同时给出权重分配为1：1时的优化艇形结果，如图9,表3所示。(下转第113页）表3权重为1：1时优化艇形与原始艇形对比初始最大截面直长细比优艇长优艇形径优化艇形化艇形化艇形阻力(N)436.8340384.3684394.3462416.9757阻力优化率12.01%9.73%4.55%体面比3.53682.95873.30403.4618体面比优化率-16.35%-6.58%-2.12%权重为1:1时NPL阻力型飞艇优化外形初始艇形最大截面直径优化艇形长细比优化艇形艇长优化艇形44020-10102080X-飞艇长度图9权重为1：1时飞艇优化外形113上接第6 5页）Meeting.IEEE,2015:1-5.IEEE 1159-2009.7公茂法，侯林源，梁龙金，等，王宁.基于小波包和PNN的电能质量扰动定位与分类J.工矿自动化，2 0 16,4 2（5）：4 0-4 4.8覃星福，龚仁喜.基于广义S变换与PSO-PNN的电能质量扰动识别J.电力系统保护与控制，2 0 16,4 4(15）：10-17.9Faramarzi A,Heidarinejad M,Mirjalili S,et al.Marine predators al-gorithm;a nature-inspired MetaheuristicJ.Expert Systems withApplications,2020,152(152).10Sahani M,Dash P K.Automatic power quality events recognitionbased on Hilbert Huang transform and weighted bidirectional ex-treme learning machine J.IEEE Transactions on Industrial In-formatics,2018,14(9):3849-3858.11Jaekel B W.Electromagnetic environments Phenomena,classifica-tion,compatibility and immunity levels C.Eurocon,Eurocon 09从上述表1、2、3的仿真分析可以看出，优化过程中体面比与阻力呈现相反优化增量关系，故在实际应用中应根据实际需求分配权重因子，此外，本文考虑实际优化设计的三个关键变量，最大截面直径、长细比和飞艇长度，分别有侧重地提出优化方案，从优化结果来看，改进算法在飞艇优化中是可行的。5乡结论本文结合传统飞艇优化设计过程中单一考虑阻力因素的问题，首先引人体面比，构建优化函数数学模型，构建优化流程框架，其次针对传统最优值求解的复合形法进行改进，依照双向同步探索思维，提出最优值求解的改进算法，并将其应用至飞艇外形优化设计中，运用Matlab进行函数的计算和仿真，通过仿真分析验证了算法的可行性。参考文献：1 Anthony Colozza,James LDolce.High-Altitude,Long-EnduranceAirships for Coastal SurveillanceR.NASA PTM-2005-21347,2005.2王全保，陈吉安，段登平，付功义.平流层飞艇外形的设计优化J.计算机仿真，2 0 10,2 7(9)：4 4-4 7.3杨燕初，王生，顾逸东，李英堂。基于遗传算法的临近空间飞艇多学科优化设计J.计算机仿真，2 0 12,2 9（4）：4 9-53.4秦何军，程志航，曾友兵飞艇气动外形的工程优化设计方法C.沈阳：第六届中国航空学会青年科技论坛，2 0 14：10 9-113.5张超，席懿.重载飞艇低阻囊体外形优化J.科学技术创新，2016(15):49.IEEE.2009.作者简介陈一镖（19 9 2-），男（汉族），浙江省温州人，博士研究生，讲师，硕士研究生导师，主要研究领域为电能质量扰动检测技术。倪陈义（19 9 6-），男（汉族），浙江省温州人，硕士研究生，主要研究领域为电能质量扰动检测技术。陈浩（19 9 8-），男（汉族），浙江省温州人，硕士研究生，主要研究领域为群体智能优化算法。李彬彬（19 8 9-），男（汉族），浙江省温州人，工程师，主要研究领域为光伏逆变器控制、有源滤波器产品控制算法实现。6刘惟信.机械最优化设计M.北京.清华大学出版社.19 9 4.9.7Nelder J A,Mead R.A simplex method for function minimizationJ.Computer J.1965,(7):308-313.8马昌凤.最优化方法及Matlab程序设计M北京.科学出版社.2 0 10.9杨晋.无约束优化问题的再改进单纯形法J.太原重型机械学院学报.19 9 0(11):2 6-31.10李庆高.非线性约束极值问题的改进单纯形法J.经济数学.1990:58-62.11杨燕初，王生，顾逸东，李英堂，苗景刚.临近空间飞艇外形优化设计与仿真J.中国空间科学技术,2 0 11（5）：7 6-8 3.12Khoury G A,Gillett J D.Airship Technology M.Cambridge U-niversity Press,1999.13安伟刚，李为吉，王海峰.某型飞艇外形多目标优化设计及其决策J.西北工业大学学报，2 0 0 7,2 5（6）：7 8 9-7 9 3.作者简介曹胜鸿（19 9 8-），男（汉族），河北省三河市人，硕士研究生，主要研究领域为飞行器结构设计，CAD/CAE。杨燕初（19 8 1-），男（汉族），河南省郑州市人，研究员，博士研究生导师，主要研究领域为浮空飞行器总体与结构。张天晓（19 9 7-），男（汉族），河北省石家庄市人，硕士研究生，主要研究领域为飞行器结构设计，浮空器装备设计与制造。张航悦（19 9 4-），男（汉族），四川南充人，博士研究生，主要研究领域为飞行动力学仿真。