代理模型在飞行器多学科优化中的应用研究_蒋鲁佳.pdf

第 43 卷第 3 期2023 年 6 月弹箭与制导学报Journal of Projectiles，ockets，Missiles and GuidanceVol.43 No.3Jun.2023DOI:10 15892/j cnki djzdxb 2023 03 008收稿日期:2022 08 10作者简介:蒋鲁佳(1979)，男，高级工程师，博士，研究方向:飞行器设计。引用本文:蒋鲁佳，辛万青，张鸣，等 代理模型在飞行器多学科优化中的应用研究 J 弹箭与制导学报，2023，43(3):56-61JIANG Lujia，XIN Wanqing，ZHANG Ming，et al Application of surrogate model in aircraft multidisciplinary design optimization J Journal of Projectiles，ockets，Missiles and Guidance，2023，43(3):56-61代理模型在飞行器多学科优化中的应用研究蒋鲁佳，辛万青，张鸣，王威，陈红波(北京宇航系统工程研究所，北京100076)摘要:飞行器多学科优化能够提高飞行器总体设计的质量和综合性能，然而，在优化设计过程中，庞大的数据量与计算量、复杂的信息交互、非线性设计空间等会耗费大量的优化计算时间，从而降低优化效率及优化结果的准确性。文中针对飞行器多学科设计模型的特点，从工程应用的角度出发，通过构造代理模型，简化了优化迭代过程中需要进行的复杂的学科分析，有效提升了飞行器多学科设计优化的效率。关键词:飞行器总体设计;多学科设计优化;代理模型;总体设计中图分类号:V19文献标志码:A文章编码:1673-9728(2023)03-0056-06Application of Surrogate Model in Aircraft MultidisciplinaryDesign OptimizationJIANG Lujia，XIN Wanqing，ZHANG Ming，WANG Wei，CHEN Hongbo(Beijing Institute of Astronautical Systems Engineering，Beijing 100076，China)Abstract:Aircraft system multidisciplinary design optimization can improve design quality and integrated performance How-ever，huge amounts of calculation，complex data exchanges，and complicated nonlinearities require a lot of calculation time，which may reduce optimization efficiency and accuracy of optimization An engineering way to simplify complex subjects analy-sis during optimization iteration with surrogate model is presented，according to the characteristics of system multidisciplinarydesign model of aircraft esults show that the method can improve aircraft system multidisciplinary design optimization efficiencyKeywords:aircraft system design;MDO;surrogate model;system design0引言多学科设计优化是为解决复杂系统设计困难提出的概念，通过实现多个学科的自动设计和优化能够缩短产品研制周期、提高设计质量、降低研制成本，具有重要的研究意义和应用前景1，尽管在很多领域得到了应用，并不断向规范化、商业化发展，然而，在多学科设计优化过程中，大量的计算、学科之间复杂的数据交互以及学科分析中非线性响应都会消耗大量的计算时间，从而降低优化算法的寻优效率以及优化结果的准确性。代理模型技术通过构造近似函数可以将复杂的学科分析简化，从而避免了在优化过程中进行的极其耗时的多学科分析，从而降低了设计空间搜索的工作量。目前该技术已应用于风力涡轮机、巡航导弹头部外形、运输机机翼等优化设计中，并取得了较好的设计结果2-4，此外，还有人提出了一些基于代理模型技术的优化算法，并验证该方法能够有效提升优化效率5-6。文中以飞行器多学科设计优化为例，探寻代理模型技术在其中应用的思路和方法。1代理模型原理代理模型主要有多项式响应面模型、Kriging 模型、径向基函数模型以及人工神经网络模型等。1)多项式响应面模型多项式响应面的前几项数学表达式为:f(x)=0+mi=1ixi+mi=1mjiijxixj(1)式中:xi为 m 维自变量 x 的第 i 个分量;0，i，ij为待定参数。多项式响应面模型具有良好的可微性和连续性，易于实现优化，且可以通过系数判断各项对系统影响的大小7，但是，在处理复杂程度较高的非线性问题时，多项式响应面模型预测和拟合的结果相对较差。第 3 期蒋鲁佳等:代理模型在飞行器多学科优化中的应用研究此外，当建立模型的阶数较高时，过拟合的情况也有可能发生8。2)Kriging 模型Kriging 模型一般表示为:f(x)=g0(x)+z(x)(2)式中:g0(x)视作模型中确定的部分，可用多项式来表示;z(x)为涨落，是一个均值为 0、方差为 2的随机函数。Kriging 模型处理复杂程度较高的非线性问题时，往往能够取得较好的结果，但是构造模型所消耗的时间较长，且由于算法复杂，较难通过软件实现，因此在实际使用中会受到一定限制7，9-11。3)径向基函数模型径向基函数的基本形式为:f(x)=ni=1wi(r)=wT(3)式中:w(x)=w1w2wn;=(r1)(r2)(rn)，wi为权系数;(r)为径向函数;ri=x xi为样本点 xi与待测点 x 之间的欧几里得度量12。径向基函数模型的基本思路是首先确定一组样本点，然后再以这些样本点为中心，以径向函数作为基函数，通过把这些基函数线性叠加后计算待测点处的响应值13，径向基函数模型会根据选取径向函数的变化而变化14。径向基函数模型是一个结构相对简单且计算量较小的代理模型15-16。4)神经网络模型神经元的结构模型如图 1 所示。图 1神经元结构模型Fig 1Neuron structural model其中，wi为神经元与输入 xi间的连接权重;i为神经元的阈值;xi为神经元的输入信号;f 为传递函数17。神经网络模型使用输入和输出样本训练网络，通过调节各层神经元之间的权系数，实现模型的构建。神经网络模型建立的过程较其他代理模型更为简单，且可以通过实现并行计算来提高效率18-20。2飞行器多学科设计建模文中以飞行器为背景，采用由最小起飞质量和最大射程组成的多目标函数作为飞行器总体多学科设计优化的目标函数，并将影响该目标函数以及设计方案可行性的几个专业，包括发动机结构及内弹道设计、飞行器外形设计、飞行气动设计、飞行器结构设计、质量参数分析、飞行轨迹设计、姿控设计以及飞行载荷设计作为实现飞行器总体多学科设计优化集成的学科21。各学科使用方案论证阶段的设计方法或能够满足该阶段要求的工程计算方法建立各学科设计模型。通过封装每个学科模型，使各学科成为具备独立分析和设计能力的模块。对建立的各学科设计模型的输入、输出数据进行分析，根据分析结果建立的飞行器总体多学科设计关系如图 2所示。图 2飞行器总体多学科设计关系Fig 2elationship of aircraft system design subjects从图 2 中可以看出多学科设计模型之间数据传递关系复杂，且存在多个反馈回路，而反馈回路会为飞行器总体多学科设计系统分析和优化带来较大的困难，为此，采用多学科设计优化耦合关系处理方法22 进行解耦，解耦后的飞行器总体多学科设计关系消除了原系统中存在的反馈回路，如图 3 所示，系75弹 箭 与 制 导 学 报第 43 卷统在进行多学科分析时，只需由设计结构矩阵左上角的学科开始，依次进行分析即可23。图 3调整后的飞行器总体设计 DSMFig 3DSM of aircraft system design after adjustment使用 iSIGHT-FD 软件平台对飞行器总体多学科设计进行集成，集成的流程图和数据关系分别如图 4和图 5 所示24。图 4飞行器总体多学科设计集成的流程图Fig 4Flow chart of integrated aircraft systemmultidisciplinary design图 5飞行器总体多学科设计集成数据关系Fig 5elationship of integrated aircraft systemmultidisciplinary design3飞行器多学科设计优化采用非支配遗传算法对建立的飞行器总体多学科设计优化模型的设计空间进行搜索，设置的种群数量为 50，迭代次数为 50，交叉概率为 0 9，模型运行共耗时 91 h9 min51 s，平均每次迭代用时约 2 min11 s，得到的优化结果如图 6 所示，其中，纵坐标代表飞行器航程，横坐标代表飞行器起飞质量，图中隐去各坐标轴的度量(下同)，图中蓝色的点即为Pareto解集，红色的点为满足各种约束条件的可行解集。图 6优化结果Fig 6esult of optimization由以上优化结果可以看出，Pareto 解集大致形成了一个 Pareto 面，由于 Pareto 解集中的点是由优化迭代范围内的非劣解集组成，因此，所有可行解包括优化的初始点都在这个 Pareto 面的右下方。在纵轴或横轴的任意位置画一条与横坐标轴或纵坐标轴平行的直线，该直线与 Pareto 面的交点就是在射程一定的条件下起飞质量最小的单目标优化问题或在起飞质量一定的条件下射程最大的单目标优化问题的最优解，Pareto 解集也正是由这些“起飞质量一定的最大射程”或“射程一定的最小起飞质量”单目标优化问题的最优解组成。通过对优化结果分析发现，飞行器总体多学科设计优化 Pareto 解集中两个端点设计变量的变化非常小，特别是一些外形参数变化都是在 1 5 mm 左右，这对于工程设计来说，只相当于在实际生产制造过程中产生的误差，最终得到的起飞质量和射程目标函数85第 3 期蒋鲁佳等:代理模型在飞行器多学科优化中的应用研究也非常的接近。经进一步分析可知，虽然耗费了大量的计算时间，进行了 2 500 次迭代计算，也得到了Pareto解集，但这仅是在 2 500 次迭代结果中得到的Pareto 解集，迭代次数与设计空间上的可行解集相比还是非常有限的，尽管选择的非支配排序遗传算法在设计空间上具有全局搜索的能力，即理论上只要有充足的时间一定可以得到最优的 Pareto 解集，但是平均每次迭代需要耗费 2 min11 s 的模型，要进行全局搜索消耗的计算时间是非常巨大的。对于建立的飞行器总体多学科设计优化模型存在的分析时间长、计算量大的问题，在当前计算机发展水平下，采用直接优化的方法获得全局最优解还存在一定的困难。为此，提出使用代理模型对复杂的多学科设计模型进行近似，并基于此开展优化设计的方法。4基于代理模型的飞行器多学科设计优化将飞行器总体多学科设计优化模型中各学科分析模型作为一个整体建立代理模型，即在不考虑其内部各学科设计模型的输入输出以及各学科之间数据传递关系的条件下，建立飞行器总体多学科设计模型的代理模型。首先使用了典型样本点分别采用响应面代理模型和神经网络代理模型对飞行器总体多学科设计模型进行近似。构造的四阶响应面代理模型和神经网络代理模型误差分析部分结果分别如图 7、图 8 所示，其中，左侧“Total Error”列为构造代理模型所有对应的样本点平均误差归一化后的数值，右侧响应图则清晰的表达了代理模型响应变量偏离实际值的情况。图 7四阶响应面代理模型误差分析Fig 7Error analysis of quartic SM图 8神经网络代理模型误差分析Fig 8Error analysis of BF model由以上误差分析对比可以看出，由于建立的飞行器总体多学科设计模型较为复杂，采用响应面法构造的代理模型误差较大，已经不能很好的反映原模型的真实情况;而采用神经网络构造的代理模型，虽然模型文件规模较大，但能够对样本点实现较好的近似，为此，选择神经网络代理模型实现对原模型的近似。为了使样本点集具有良好的散布均匀性，并在取值范围内能够以较少的组合方式获得设计空间较全面的认识，使用最优拉丁超方设计方法进行了 DOE(design of experiment)实验25，并将该实验结果作为构造代理模型的样本点，并通过训练完成飞行器总体多学科设计神经网络代理模型的构造。然后，再对构造的代理模型进行优化，共耗时约 2 h，进行了 40 000次迭代，得到的 Pareto 解集如图 9 所示。从图中可以明显的看出起飞质量与射程的变化趋势，即随着装药量的增加，获得的射程也越大。其中，起飞质量在变化的过程中存在一个跳变，这可能是由于在构造代理模型样本点不够充分引起的。图 9代理模型优化结果Fig 9esult of surrogate model optimization95弹 箭 与 制 导 学 报第 43 卷在 Pareto 解集中选取了某起飞质量下的设计结果，并将该结果重新在飞行器总体多学科设计模型中进行分析计算，得到该起飞质量下的最大射程设计结果。通过比对发现，基于代理模型的飞行器多学科设计优化结果较原来的设计结果将最大射程提高了约6%，较好的说明基于代理模型的飞行器多学科设计优化在大幅提升优化效率的同时，还能够充分利用各学科模型相互作用产生出来的协同效应，进一步挖掘设计潜力26-27。5结论文中针对飞行器总体多学科、多目标、多约束优化设计带来的计算困难问题，首先对复杂多学科优化问题计算困难进行了分析;其次提出一套使用代理模型技术实现飞行器总体多学科设计优化的思路和方法;然后提出使用代理模型将飞行器总体多学科设计模型作为一个整体进行近似的方法;最后提出了使用最优拉丁超方设计方法对设计空间进行 DOE 实验，并将实验结果作为建立代理模型的样本点的方法。优化结果表明，该方法可以有效降低多学科设计优化问题的复杂性，提高优化效率。构造代理模型进行优化能够缩短计算需要的时间，虽然优化结果与原优化模型存在一定差距，但是随着构造代理模型的样本点增加，理论上能够在一定精度范围内反映设计变量与目标函数之间的大致关系，为设计人员权衡折衷、初步选择初始设计参数提供参考。此外，将代理模型的优化结果作为精确模型的初始设计点进行优化，也能够加快整个优化过程的收敛速度。参考文献:1 王振国，陈小前，罗文采，等 飞行器多学科设计优化理论与应用研究M 北京:国防工业出版社，2006:1 25WANG Z G，CHEN X Q，LUO W C，et al esearch onthe theory and application of multidisciplinary design opti-mization for aircraft M Beijing:National Defense Indus-try Press，2006:1 25 2 EFOT G，BACKSTOM T W V，VENTE G Numeri-cal optimisation of a small-scale wind turbine through theuse of surrogate modelling J Journal of Energy in South-ern Africa，2017，28:79 91 3 GUO S Z，ZHENG X M，ANG H S，et al Cruise missilehead shape optimisation using an adaptive sampling surro-gate modelJ The Aeronautical Journal，2018，122:1145 1162 4 李育超，齐婵颖，高通锋 基于 Kriging 代理模型的运输机机翼多学科优化设计 J 航空科学技术，2018，29(3):20 24LI Y C，QI C Y，GAO T F Multi-disciplinary optimizationof transport aircraft wing based on kriging surrogate model J Aeronautical Science and Technology，2018，29(3):20 24 5 MULLE J SOCEMO:surrogate optimization of computa-tionally expensive multiobjective problemsJ InformsJournal on Computing，2017，29:581 596 6 LEE Y H，COMAN E，EWOLDT H，et al A mul-tiobjective adaptive surrogate modeling-based optimization(Mo-Asmo)framework using efficient sampling strategies C/ASME Proceedings of the 2017 International DesignEngineering Technical Conferences Computers Infor-mation in Engineering ConferenceNew York:ASME，2017:1127 1139 7 李伟华 基于遗传算法的翼型性能优化 D 武汉:华中科技大学，2008LI W H Optimization of aerofoil performance based on thegenetic algorithmD Wuhan:Huazhong University ofScience and Technology，2008 8 余德海，宋保维，李家旺，等 鱼雷多学科设计优化中的代理模型研究 J 机械科学与技术，2007，27(2):184 187YU D H，SONG B W，LI J W，et al Surrogate models forthe multi-disciplinary design optimization of a torpedoJ Mechanical Science and Technology for Aerospace Engi-neering，2007，27(2):184 187 9 LIU X Y，YUAN Q Q，ZHAO M，et al Multiple objectivemultidisciplinary design optimization of heavier-than-waterunderwater vehicle using CFD and approximation model J Journal of Marine Science and Technology，2017，22(1):135 148 10 高轶男 基于组合代理模型的汽车碰撞安全性多目标优化研究 D 重庆:重庆大学，2016GAO Y N esearch on multi-objective optimization of ve-hicle crash safety based on combined surrogate model D Chongqing:Chongqing University，2016 11 张磊 基于高斯过程代理模型的粒子群算法研究D 杭州:浙江科技学院，2019ZHANG L esearch on particle swarm optimization basedon gaussian process surrogate model D Hangzhou:Zhe-jiang University of Science and Technology，2019 12 吴梦阳 某型装甲车底部防护组件结构材料一体化设计技术研究 D 南京:南京理工大学，2019WU M Yesearch on integrated design technology ofstructural materials for bottom protection components of acertain armored vehicleD Nanjing:Nanjing Universityof Science and Technology，201906第 3 期蒋鲁佳等:代理模型在飞行器多学科优化中的应用研究 13 刘彬，肖守讷，朱涛，等 自由外翻式组合吸能装置的设计及优化 J 中国机械工程，2018，29(17):2052 2059LIU B，XIAO S N，ZHU T，et al Design and optimizationof combined energy absorption devices based on free inver-sion J China Mechanical Engineering，2018，29(17):2052 2059 14 赵涛 基于代理模型验证的混合不确定性优化方法及其在空间绳系设计中的应用D 成都:电子科技大学，2015ZHAO T Hybrid uncertainty optimization method based onsurrogate model verification and its application in spacerope system design D Chengdu:University of ElectronicScience and Technology of China，2015 15 KOCH P N，WUJEK B Facilitating probabilistic multidis-ciplinary design optimization using kriging approximationmodels:AIAA 2002-5415 eston:AIAA，2002 16 杨路春，胡迟，高阳，等 基于多目标优化和多属性决策的数据采集装置总体优化设计 J 舰船科学技术，2021，43(1):164 169YANG L C，HU C，GAO Y，et al Data acquisition unitoptimal design based on multi-objective optimization andmulti-attribute decision makingJ Ship Science andTechnology，2021，43(1):164 169 17 杨华 机翼气动力双向代理模型设计 D 南京:南京航空航天大学，2007YANG H Design of a two-way surrogate model aerodynam-ic model for wingsD Nanjing:Nanjing University ofAeronautics and Astronautics，2007 18 王洪元，史国栋 人工神经网络技术及其应用 M 北京:中国石化出版社，2002:21 58WANG H Y，SHI G D Artificial neural network technologyand its applicationsM Beijing:China PetrochemicalPress，2002:21 58 19 梁晗 复合地层小半径曲线隧道盾构掘进关键技术研究 D 北京:北京交通大学，2020LIANG H Study on the key technology of shield tunnelingin small radius curve tunnel in composite stratum D Bei-jing:Beijing Jiaotong University，2020 20 葛学伟 气体辅助注射成型穿透机理及工艺优化方法研究 D 上海:上海交通大学，2010GE X W Mechanism study and process optimization on gaspenetration behavior for gas-assisted injection molding D Shanghai:Shanghai Jiao Tong University，2010 21 周军，林鹏，周敏，等 降低气动耦合特性的飞行器总体随控优化方法:201310060969 3 P 2013 07 10ZHOU J，LIN P，ZHOU M，et al Optimization method foroverall control of aircraft to reduce aerodynamic couplingcharacteristics:201310060969 3 P 2013 07 10 22 蒋鲁佳，辛万青，布向伟 导弹总体多学科设计优化耦合关系处理方法 J 导弹与航天运载技术，2009(6):5 7JIANG L J，XIN W Q，BU X W Processing methods ofcoupled relationship of multidisciplinary design optimizationfor missile J Missiles and Space Vehicles，2009(6):5 7 23 刘峰 载人潜器总体设计优化方法研究D 哈尔滨:哈尔滨工程大学，2014LIU F esearch on general design optimization method ofhuman occupied vehicle D Harbin:Harbin EngineeringUniversity，2014 24 蒋鲁佳 多学科设计优化在飞行器总体设计中的应用研究 J 弹箭与制导学报，2012，32(3):39 43JIANG L J The application of multidisciplinary design opti-mization in aircraft system designJ Journal of Projec-tiles，ockets，Missiles and Guidance，2012，32(3):39 43 25 赵选民 试验设计方法 M 北京:科学出版社，2006:123 169ZHAO X M Experimental design methodM Beijing:Science Press，2006:123 169 26 施雨阳 远程空空反辐射导弹总体优化设计技术研究 D 西安:西北工业大学，2015SHI Y Y esearch on long range air-to-air anti radiationmissile general design and optimization D Xi an:North-western Polytechnical University，2015 27 马燕芳 显式模型预测控制的多胞体近似方法及其应用 D 杭州:浙江工业大学，2014MA Y F Polytopic approximation of explicit model predic-tive control and its applicationsD Hangzhou:ZhejiangUniversity of Technology，201416