2016年第37卷第2期 (总第166期) 中北大学学报(自然科学版) Vo1.37 No.2 2016 (Sum No.166) JOURNALOFNORTHUNIVERSITYOFCHINA(NATURAL SCIENCEEDmON)文章编号:1673—3193(2016)02—0137—05 基于自适应遗传算法的 空射运载火箭弹道优化方法 胡 钰 ,王华 (1.中北大学机电工程学院,山西太原030051;2.北京宇航系统工程研究院,北京100076) 摘要: 在传统的运载火箭能力优化设计基础上,针对空中发射运载火箭,将运载能力优化转换为最大速 度优化,将复杂环境转换为多种弹道约束,提出一种在满足多约束条件下优化终端速度的弹道设计方 法.将约束条件以罚函数形式引入目标函数,在传统遗传算法基础上,通过自适应调节交叉概率和变异概 率,实现了一级飞行段攻角的优化.仿真结果验证了自适应遗传算法在弹道优化上的可行性,解决了传统运 载能力优化方法设计复杂、难以满足各项过程约束的问题,为解决运载能力优化问题提供了一种新的解决 思路,对各类弹道优化的工程实践具有一定的参考价值. 关键词:空中发射;运载火箭;弹道设计;遗传算法; 自适应 中图分类号: V448.235 文献标识码:A doi:10.3969/j.issn.1673—3193.2016.02.008 Traj ectory Optimization Method of Air-Launched Vehicle Based on Adaptive Genetic Algorithm HU Yu 。WANG Hua (1.School of Mechatronic Engineering,North University of China,Taiyuan 030051,China; 2.Beijing Institute of Astronautical System Engineering,Beijing 100076,China) Abstract:Based on optimized design of the traditionallaunching vehicle capability,this paper converts the capacity optimization to maximum speed optimization,and converts complex environmental restric— tion to trajectories constraint,proposing tr ̄ectory design method of optimize terminal velocity meeting muti—constrained conditions air-launched vehicle.The constraints in the form of penalty function intro- duced in the obj ective function,optimization of primary fight angle of attack is realized by adaptive con— trol crossover probability and mutation probabiting based on traditional genetic algorithm.Simulation results testify the effectively and feasibility of the adaptive genetic algorithm on multi·。constraints trajec‘- tory optimization,thus providing a new idea of solving the problem of carrying capacity optimization and providing advantages for engineering applications. Key words:air launch;launch vehicle;traj ectory optimization;genetic algorithm;adaptive control 相比于陆地发射,空射运载火箭由于具有运 载能力大、机动范围广、突防能力强、对地面基础 轨道,为运载火箭的空射技术打开了成功之i'q c . 2003年美国国防部启动了“猎鹰”计划,其首要目 标就是采用一种低成本的小型运载火箭将卫星快 设施依赖性小、发射费用低等特点备受各国青 睐 .1990年4月5日,美国采用B-52成功发射 了“飞马座”固体运载火箭,将两颗小型卫星送入 速送到地球低轨道[引. 空射运载火箭通过载机将其运送到预定区域 收稿日期:2015—08—26 作者简介:胡钰(1987一),女,博士生,主要从事兵器科学与技术研究 138 中北大学学报(自然科学版) 2016年第2期 并达到预定高度、速度和航线后,基于运载器自 身重力进行水平投放,投放后为保证飞机的飞行 安全,需使用空气舵进行姿态控制__4],自由下落 3~5 S后一级发动机点火.然而,由于空射运载 火箭具有初始速度并且初始弹道倾角为零,为了 快速冲出大气层减小动压需进行大攻角拉起,利 用发动机推力和气动升力爬高,从而增加弹道倾 角及飞行高度.与陆地垂直发射的运载火箭相比, 一级飞行段飞行环境更为恶劣,姿控、分离、载 荷、热很容易超出其包络范围,为弹道设计提出 了更多要求. 遗传算法以自然选择和遗传理论为基础,通 过模拟物种进化过程中的复制、交叉和变异等现 象,对目标空间进行随机优化搜索,找到最适合 环境的个体L5],但是普通遗传算法存在早熟的现 象,同时其收敛速度较慢,为了解决这个问题,文 献[6]提出自适应遗传算法,实现交叉概率和变异 概率随群体的适应度自动改变.本文采用遗传算 法对弹道进行优化,在满足多学科约束的条件下 寻找到使关机点速度最大的程序攻角. 1 弹道建模 1.1火箭模型 在发射坐标系中建立火箭质心运动方程 ] , (1) ]一 ]+巨]一匡三]一匡兰], 2 * , (3) 式中:发射坐标系原点与发射点0固连,OX轴在 发射点水平面内,指向发射瞄准方向,o 轴垂直 于发射点水平面,指向上方,OZ轴与xo.y面相垂 直并构成右手坐标系.3.7,Y, 为火箭在发射坐标 系中的位置三分量, , , 为火箭速度在发 射坐标系上的三分量, ,W ,W 和g ,g ,g 分别为火箭视速度和引力加速度在发射坐标系中 的位置三分量,n , , 和以h,a 口 分别为牵 连加速度和哥氏加速度在发射坐标系上的三分量, m为火箭质量, ,P P 和R ,R ,R 分别 为发动机推力和气动力在箭体坐标系上的三分量, GB为箭体坐标系到发射坐标系的转移矩阵. 1.2总体思路 ~ 对于运载火箭弹道设计而言,最重要的指标 为运载能力的大小.根据齐奥尔科夫斯基公式, 为讨论方便,假设各级真空比冲相等均为 ,子级 结构系数相等均为 ,火箭的级问比相等均为e=== ,E为火箭载荷比,即火箭起飞总质量与有效 质量之比.对于级数为 的火箭,理想速度可由 式(4)估算 Vr—niln—————l(1一正■ )+点■ _——— . (4) 由式(4)可以看出,对于以固定速度入轨的弹 道,运载能力增大等效为固定运载能力终端速度 增大.对于空中发射运载火箭,二级飞行段以上 由于约束条件较少,采用零攻角飞行可使推力方 向与速度方向相同,进而使能量最大限度地用于 终端速度地增加上,而一级飞行段由于环境恶劣, 约束条件较多,因此需要根据不同的约束要求对 弹道进行调整,运载能力的优化问题转换为一级 飞行段弹道的优化问题. 一级飞行段恶劣的环境及诸多约束条件为小 回路的闭合增加了难度.传统上,在弹道专业完 成弹道设计后,将弹道结果传输给气动、载荷、姿 控等专业进行复算,若不满足要求,需重新调整 弹道.反复的迭代不仅效率低,同时难以设计出 满足各专业要求的最优弹道. 本文转换思路,通过寻找气动、载荷、姿控等 与弹道设计的关系,将环境约束转变成过程约束, 实现一级飞行弹道设计. 1.3弹道设计方案 火箭在稠密大气层中的飞行任务是爬高和加 速,因此轨道设计一方面要充分利用火箭发动机 推力和气动升力进行爬高,另一方面要力求减少 气动阻力的影响__5]. 首先,考虑到投放后飞机的飞行安全,在投 放初始3 S采用O。攻角进行飞行;其次,发动机点 火后由于飞行速度较小,可采用较大攻角进行迅 速爬升;之后,火箭获得一定高度和速度后进行 两次下压,通过减小攻角实现动压和载荷的降低; 最后,为姿控和分离创造条件,在发动机耗尽前 将攻角归零.其程序攻角示意图如图1所示,其 (总第166期) 基于自适应遗传算法的空射运载火箭弹道优化方法(胡钰等) 139 中横坐标t为时间,纵坐标a为程序攻角, , 固定不变,因此存在早熟和收敛速度较慢的问题, 自适应遗传算法根据群体自适应度设计交叉概率 P 和变异概率P ,使得适应度高于群体平均适 应度的个体能够IJ ̄.N遗传至下一代,而适应度低 Otcon2,a。 为程序攻角设计的3次定轴值. 0 叫. 于群体平均适应度的个体尽量被淘汰.式(5)给出 ‰2 0 蛳 0 图1程序攻角示意图 Fig.1 Structure of advanced guidance system 各段攻角上抬(下压)开始时刻及结束时刻由 各段攻角定轴值、火箭姿态控制能力及高度决定, 因此一级飞行段程序攻角由3次定轴值O/con、, Otcon2,a 决定. 1.4优化变量及约束条件 由上文可知,空射火箭一级飞行段优化变量 为三次程序攻角定轴值O/con,,Otoc ̄,O/con3,优化目标 为一级飞行段关机点速度 最大.在一级飞行 过程中,主要受到气动热和载荷的,其中气 动热主要受飞行过程中马赫数及攻角的影响Is3, 根据载荷设计原理,火箭结构由平移和转动加速 度产生的惯性力与气动力平衡,使火箭产生沿箭 体纵轴的分布剪力和弯矩为 n F1 一 QSaf(C ̄,zT,zy,△z卜1)一 =1 QSaf(Cu ̄, r,,z ,△z卜1), n M1 一 F1,卜1△z卜1+QSaf(C ̄,oTT, ,△z卜1), 一1 式中: 一1,2,…,N为火箭离散点数;Q为飞行 动压;S为火箭气动参考面积;f( ,z ,37.,, H)是总法向气动系数导数、火箭质心、火箭气 动压心和分站站点的函数. 可以发现,对于弹道设计,具体表现为最大 动压值Qr 及动压与攻角乘积的最大值Qa ,同时考虑到二级姿控系统控制能力,还要 对一级关机点动压Q 一z进行约束. 2遗传算法的实现 2.1自适应遗传算法 传统遗传算法中交叉概率P。和变异概率P 了P 和P 的求解公式[ . P 一 f 一厂a ≥ ,’ l 足z 厂< , Pm =:::== j ,一x一  ̄avg f ̄f.vg, (5) 【 是 厂< , 式中:厂衄 为群体中个体适应度最大值; 为群 体的平均适应度; 为要交叉的两个个体中较大 的适应度;f为要变异个体的适应度;k ,kz,ks,k 为介于(O,1)之间的常值. 2.2算法步骤 由上文可知,弹道优化的目标是一级飞行段 关机点速度最大,同时Q眦 , 一和Q 一 不大于 约束值,设置遗传算法的目标函数为 一一 nd+(max(0, 一Q一 )) + 0.1(max(0,Q 一Qrllaxc)) + (max(0,Q1—2一Q1-2c)) , (6) 式中:Q , 和Q 分别为一级飞行过程中 的动压最大约束值、动压与攻角乘积的最大约束 值和关机点动压最大约束值. 3仿真结果分析 采用本文设计的优化方法,对一级飞行段攻 角进行优化,假设Q口 ===1 200 000 Pa·deg, QI 一120 000。,Q1_2c一80 000。,同时由于风洞试 验的局限及姿态控制能力的,优化变量O/con1, a ,a ,均不大于25。.图2为采用遗传算法得到 的攻角随时间变化曲线,图3和图4分别为动压 和动压与攻角乘积随时间变化曲线. 从仿真结果可以看出按照图1的攻角设计思 路,采用带罚函数的自适应遗传算法可以在满足 各项过程约束的同时,实现关机点速度最大,很 好地解决了空中发射运载火箭一级飞行段飞行环 境恶劣、姿控、分离、载荷、热很容易超出其包络 范围的问题,为小回路的设计提供了方便. 140 中北大学学报(自然科学版) 2016年第2期 一一。一南≈d)I\ t/S 图2攻角随时间变化曲线 Fig.2 Curve of angle of attack changes with time 图3动压随时间变化曲线 Fig.3 Curve of dynamic pressure changes with time t/s 圈4动压与攻角乘积随时间变化曲线 igF.4 Curve of the product of dynamic pressure with angle of attack changes ith tiwme 4结 论 本文将空射运载火箭运载能力优化转换为一 级飞行段关机点速度的最大,将一级飞行过程恶 劣的飞行环境转换为弹道设计的过程约束.以某 一特定空射火箭为背景,将约束条件以罚函数形 (总第166期) 基于自适应遗传算法的空射运载火箭弹道优化方法(胡钰等) 141 式引入目标函数,通过自适应遗传算法优化求解 得到满足约束条件的最大终端速度.仿真结果表 明,采用此方法可以在工程上解决此问题,同时, 可以将此运用到其他弹道设计中,对工程实践具 有一定的研究价值. 参考文献: Eli翔天.空射弹道导弹/运载火箭技术深度解析[J].现 [4]茹家欣.空中发射运载火箭发射轨道设计[J].现代防 御技术,2004,32(2):28—30. Ru Jiaxin.Trajectory design for air launch vehicle[J]. Modern Defense Technology,2004,32(2):28—30.(in Chinese) [5] 玄光男,陈润伟.遗传算法与工程优化[M].北京:清 华大学出版社,2004. nivas M,Patnaik L M Adaptive probabilities of [6] Sricrossover and mutation in genetic algorithm[J].IEEE 代兵器,2010,27(4):25—36. Trans on Systems,Man&Cybernetics,1994(24): Xiang Tian.Air-launched ballistic missile/rocket tech— 656. nology depth analysis[-J].Modern Weaponry,2010, [7]贾沛然.远程火箭弹道学[M].湖南:国防科技大学 27(4):25—36.(in Chinese) 出版社,2004. [2]一九九零年世界航天十大新闻[J].世界导弹与航天, [8]杨恺,高效伟.高超声速气动热环境工程算法[J].导 1991,14(1):4. 弹与航天运载技术,2010,39(4):19—28. In 1 990 top ten news of the world space[J ̄.Missiles Yang Kai,Gao Xiaowei.Engineering algorithm for &Spacecraft,1991,14(1):4.(in Chinese) aeroheating environment of hypersonic aircrafts[J]. [3]钱锟.空射小型运载火箭飞行试验项目[J].国际航空 Missiles and Space Vehicles,2010,39(4):19—28.(in 杂志,2008,53(8):21—23. Chinese) Qian Kun.Airdrop small launch vehicle flight test [9]韩瑞锋.遗传算法原理与应用实例[M].天津:兵器 program[J].International Aviation,2008,53(8):21— 工业出版社,2009. 23.(in Chinese) (上接第136页) [2]柯铧,柯科.模糊综合评判在步作战效能评估中的 系统中的应用[J].系统仿真学报,2007,19(4):772— 应用[J].数学的时间与认识,2011,41(5):180—185. 775. Ke Hua,Ke Ke.Applications in evaluation of the rifle Hu Liangming,Xu Cheng,Flang Jun Case-based rea— fighting efficiency by fuzzy synthetical evalution[J-1. soning and its application in scheme design of small Mathematics in Practice and Theory,2011,41(5): arms expert system[J].Journal of System Simulation, 180—185.(in Chinese) 2007,19(4):772—775.(in Chinese) [3]顾丽娟,李彪,孙慧静.基于模糊聚类分析的步作 [73楼俐,徐诚,张飞猛.灰色系统理论及量化分析法在 战效能评估[J].海军航空工程学院学报,2014,29 械方案设计评价中的应用[J].南京理工大学学报, (3):280—284. 2010,34(5):581—585. Gu Lijuan,Li Biao,Sun Huijing.Evaluation of the ri— Lou Li,Xu Cheng,Zhang Feimeng.Application of fie fighting efficiency based on fuzzy clustering analysis grey system theory and quantitative analysis method to [J].Journal of Naval Aeronautical and Astronautical scheme design evaluation of small arms[J].Nanjing University,2014,29(3):280—284.(in Chinese) University of Science Technology,2010,34(5): [43腾臣.步作战效能评估方法研究[D].南京:南京理 581—585.(in Chinese) 工大学,2014. [8]朴尚哲,超木日力格,于剑.模糊C均值算法的聚类 [53郭凯,徐诚.基于武器射击效率的步系统作战效能 有效性评价[J].模式识别与人工智能,2015,28(5): 评估方法[刀.兵工学报,2007,28(2):148-152. 452—461. Guo Kai,Xu Cheng.Evaluation method of fighting ef— Piao Shangzhe,Chaomurilige,YU Jian.Cluster valid— fectiveness of small arlTls system based on shooting ef— ity indexes for FCM clustering algorithm.PR&AI[J]. ficiency[J].ActaArmamentarll,2007,28(2):148— Pattern Recognition and Artificial Intelligence,2015, 152.(in Chinese) 28(5):452—461.(in Chinese) [6]胡良明,徐诚,方峻.CBR技术在械方案设计专家