基于分段设计的动力翼伞轨迹规划作者：朱二琳高海涛赵俊杰
来源：《江苏理工学院学报》2020年第04期
摘要：針对动力翼伞的飞行特点，对飞行及回收过程的轨迹规划问题进行研究。根据实际飞行工况，在传统翼伞轨迹规划中，加入了动力翼伞的任务执行阶段。根据各阶段的几何关系，将轨迹规划问题转换为几何参数的寻优问题。以动力翼伞的回收精度最优作为优化目标，采用量子遗传算法对目标函数进行寻优计算，得出各段轨迹的几何参数。仿真分析验证了该设计方法的可行性和有效性。
关键词：动力翼伞;轨迹规划;分段设计;量子遗传算法
中图分类号：V455.2 文献标识码：A 文章编号：2095-7394（2020）04-0028-08
动力翼伞作为一类新型的柔性飞行器，在传统翼伞系统的基础上加入了动力推进装置，与传统翼伞相比，除了能够完成滑翔和转弯运动，它还能在空中维持高度、甚至爬升[1]。凭借这些优良的飞行特性，目前动力翼伞在军事侦察、农业播撒、广告宣传等领域得到了广泛的应用。
传统翼伞的归航方法一般分为三种[2-6]：简单归航法（包括径向归航和带盲区的锥形归航）、最优控制归航法、分段归航法。20世纪90年代以来，分段归航法以其方法简单、易于实现的特点被广泛地应用于实际工程。根据翼伞自身的飞行特性，将翼伞从投放到着陆的整个过程分为几个阶段，再根据各个阶段的特点，进行分段轨迹的设计，这就是工程中普遍应用的翼伞分段归航法。熊菁对传统翼伞的分段归航做了详细介绍，采用改进的遗传算法对分段轨迹进行优化[7]。张兴会考虑了能量最优，提出了基于能量约束的分段轨迹优化方法，既保证了着陆精度又节省了能量消耗[8]。然而，目前对于动力翼伞的轨迹规划问题还没见到相关的研究文献。
近年来，随着航空航天科技的发展，对翼伞自主归航的要求越来越高，而飞行前的轨迹规划有利于能耗的降低及回收，特别是对于动力翼伞，其所执行的任务并不仅限于物资的空投和回收，所以有效的轨迹规划显得尤为重要。本文采用量子遗传算法研究了基于分段设计的动力翼伞的轨迹规划问题，并通过仿真对所提算法的可行性进行了验证。
1 动力翼伞质点模型
对于飞行器的轨迹规划问题，复杂的动力学模型将增加运算的复杂度，会影响轨迹的准确性[9]。本文采用较为简单的3自由度质点模型，在建模之前，先做如下几点基本假设：
（1）在飞行过程中，将动力翼伞看作一个质点，不考虑姿态的变化，水平速度保持不变;
（2）在任务执行阶段，动力翼伞等高飞行，垂直速度为零，在其他阶段，垂直速度为定值;
（3）只考虑水平面的风场，且风场信息已知。
如图1所示为动力翼伞水平面的飞行示意图，大地坐标系原点为着陆的目标点，X轴的方向同风场方向保持一致，Z轴方向垂直向下，Y轴同X轴、Z轴共同构成右手坐标系。
根据飞行示意图，可将动力翼伞模型简化成如下形式：
其中：（[x y z]）代表动力翼伞在大地坐标系下的位置;[Vs]和[Vz]分别代表动力翼伞的水平速度和垂直速度（在任务执行阶段[Vz=0]）;[Ψ]代表水平飞行方向同大地坐标系X轴之间的夹角;[Ψ]代表转弯角速度;[u]代表伞绳控制量（同转弯角速度之间存在对应关系）;[W]代表水平风场信息。
对于动力翼伞而言，任务执行段的轨迹应是固定不变的，不能随风产生漂移，所以动力翼伞的轨迹规划问题不能在风坐标系下进行讨论。本文的研究基于上述的大地坐标系，在轨迹规划过程中，不考虑风的因素。
2 动力翼伞分段轨迹设计
2.1 分段策略
传统翼伞的归航轨迹可大致分为三个阶段[7]：向心段，从翼伞初始投放点到着陆区域;能量控制段，在着陆区上空盘旋削高直到满足一定高度条件;着陆段，寻找逆风方向、滑翔，最后进行雀降。对于动力翼伞的轨迹规划问题，是在传统翼伞的归航过程中加入特定的任务执行段，当任务执行完毕后再进行归航。动力翼伞目前主要应用于特定的工况，通常将动力翼伞的任务执行阶段设定为拥有固定圆心和半径的圆。
动力翼伞水平方向运动只能通过单拉伞绳和双拉伞绳来进行控制，所以动力翼伞水平面的轨迹可以看作是直线和圆弧的组合。图2为动力翼伞分段轨迹的结构示意图。
如图2所示，A点是动力翼伞的初始点，O点是着陆的目标点，BC段和DE段是动力翼伞的向心段，分别对应任务执行区域和着陆区域。以[O3]点为原点、[RTP]为半径的圆形轨迹为动力翼伞的任务执行阶段，此阶段动力翼伞高度保持不变。FG段为能量控制段，此阶段动力翼伞盘旋削高。HO段为着陆阶段。AB、CD、EF、GH段为各阶段的过渡阶段，动力翼伞以最小的转弯半径转弯飞行。[βii=1，…，5]代表最小半径圆弧对应的圆心角弧度。根据各段轨迹的几何关系，动力翼伞轨迹规划的关键就是确定任务执行段的进入点D和能量消耗段的进入点F。R代表圆弧的半径，[θ]代表进入点和对应阶段圆的原点连线同大地坐标系X轴正向的夹角。下标TP代表任务执行段进入点，EP则代表能量控制段的进入点。[RTP]根据实际工况可事先设定，所以待优化的轨迹参数为[θTP，REP，θEP]。
2.2 各段轨迹的计算
[n1，n2=1]代表动力翼伞逆时针转弯，[n1，n2=-1]代表动力翼伞顺时针转弯。[x0y0z0]代表初始位置，[α0]代表初始飞行方向，[Rmin]代表过渡段最小转弯半径。
2.3 目标函数
动力翼伞归航同传统翼伞类似，需要保证以下两点：（1）着陆点距目标点误差最小;（2）着陆方向为风向的反方向。在分段轨迹设计中已经考虑了逆风对准，所以动力翼伞轨迹规划的目标函数选择为从初始点到目标点飞行的水平面轨迹同定滑翔比条件下初始高度对应的水平距离之间偏差的绝对值。由于动力翼伞在任务执行段飞行高度保持不变，因此，在目标函数中这一阶段的水平面轨迹可忽略。
目标函数如下式所示：[J=Rmin∙β1+β2+β3+β5+REP∙β4+BC+DE+REP-2Rmin∙REP-f∙z0] 。（11）
其中：[f]代表滑翔比，[f=Vs/Vz]，[z0]代表初始高度。分段轨迹的参数变量为[θTP，REP，θEP]，以下通过参数寻优的方法对轨迹规划问题进行求解。
3 量子遗传算法
量子遗传算法（Quantum genetic algorithm，QGA）[10-12]是量子计算和遗传算法相结合的产物。染色体在编码过程当中，引入量子的概率幅表示方法，从而使每个染色体能够表示多个状态的叠加。通过量子旋转门策略对染色体进行更新，可达到对目标函数的最优求解，因此，只要是遗传算法能够求解的优化问题，都可以用量子遗传算法来解决。
对于上文描述的动力翼伞分段轨迹优化问题，目标函数如式（11）所示，这是一个标准的多变量参数寻优问题，很适合用量子遗传算法来进行求解计算。根据上述问题描述，优化的参数变量为[θTP，REP，θEP]。转弯方向根据初始点位置的不同可有不同的组合方式，在优化之前进行确定，针对不同的[n1，n2]值分别对目标函数进行量子遗传算法的求解，并选取最优的解。
量子遗传算法对分段轨迹参数寻优的具体步骤如下：
（1）采用小生境协同进化策略，对种群[Q（g0）]进行初始化;
（2）测量种群[Q（g0）]中每个染色体，可得二进制編码;
（3）计算目标函数，对种群[Q（g0）]中所有个体逐个进行适应度的评价;
