发明名称 一种燃料最省的火星复杂地形区安全着陆轨迹生成方法
摘要 本发明涉及一种燃料最省的火星复杂地形区安全着陆轨迹生成方法,属于行星着陆技术领域。本发明方法首先建立火星动力下降动力学模型,然后根据目标着陆地形建立了导航函数,通过导航函数求解的避障控制力能够有效避免着陆器与障碍发生碰撞,并安全着陆到目标着陆点。本发明将求解的避障控制力引入到火星动力下降动力学方程中,控制力有一部分用于实现障碍规避,再对改造的动力学方程进行着陆轨迹优化,则能够实现障碍规避的同时节省动力下降轨迹的燃料消耗。本发明给出的火星复杂地形区的着陆轨迹优化方法既能考虑到目标着陆区的地形,也能有效降低燃料消耗,避免了传统最优轨迹和障碍规避中分别存在无法对障碍进行规避和燃料消耗多的缺点。
申请公布号 CN104267734B 申请公布日期 2016.11.02
申请号 CN201410376959.5 申请日期 2014.08.01
申请人 北京理工大学 发明人 崔平远;胡海静;朱圣英;高艾;徐瑞
分类号 G05D1/10(2006.01)I 主分类号 G05D1/10(2006.01)I
代理机构 代理人
主权项 一种燃料最省的火星复杂地形区安全着陆轨迹生成方法,其特征在于:具体包括如下步骤:步骤一,建立火星动力下降动力学模型;为了描述着陆器的运动,首先建立着陆点固连坐标系<img file="FDA0000986876760000011.GIF" wi="205" he="70" />其中原点为预定的着陆点,<img file="FDA0000986876760000012.GIF" wi="30" he="46" />和<img file="FDA0000986876760000013.GIF" wi="35" he="62" />轴在水平面,<img file="FDA0000986876760000014.GIF" wi="36" he="61" />轴竖直向上;x、y、h为着陆器位置,v<sub>x</sub>,v<sub>y</sub>和v<sub>h</sub>分别为三轴方向的速度,m为着陆器的质量,为关于时间的函数,随时间变化;将着陆器视为质点,对时间求导,建立动力下降动力学模型为:<img file="FDA0000986876760000015.GIF" wi="117" he="54" /><img file="FDA0000986876760000016.GIF" wi="126" he="62" /><img file="FDA0000986876760000017.GIF" wi="117" he="69" /><img file="FDA0000986876760000018.GIF" wi="406" he="122" /><img file="FDA0000986876760000019.GIF" wi="413" he="127" /><img file="FDA00009868767600000110.GIF" wi="510" he="119" /><img file="FDA00009868767600000111.GIF" wi="269" he="135" />其中,g和g<sub>0</sub>分别为火星和地球表面的重力加速度,n为推力发动机的数量,I<sub>sp</sub>为推力发动机的脉冲,T为推力发动机的最大推力,φ为推力发动机与合力方向的夹角,u=[u<sub>x</sub> u<sub>y</sub> u<sub>h</sub>]为三轴方向的控制向量;通过调节控制向量能够控制着陆器的运动轨迹;由于推力发动机打开后,则在整个着陆过程持续工作,因而着陆过程中控制向量需要满足如下关系:<img file="FDA00009868767600000112.GIF" wi="746" he="91" />其中u<sub>min</sub>和u<sub>max</sub>分别为控制向量的最小值和最大值;其中,u=0表示推力发动机处于关闭状态,u=1表示推力发动机以最大推力运行;步骤二,建立导航函数;根据已获得的火星地形数据,通过最小二乘拟合方式获得目标着陆区的障碍信息,构建如下形式的导航函数:<img file="FDA00009868767600000113.GIF" wi="574" he="183" />其中,q=[x y]<sup>T</sup>为着陆器的平面位置向量,q<sub>d</sub>为目标着陆点的位置向量,β(q)为障碍函数,定义为如下形式:<img file="FDA0000986876760000021.GIF" wi="326" he="135" /><img file="FDA0000986876760000022.GIF" wi="686" he="95" />其中,M为障碍的数量,q<sub>j</sub>为第j个障碍的位置,ρ<sub>j</sub>为第j个障碍的大小;利用导航函数的梯度信息生成避障力F能够阻止着陆器接近障碍:<img file="FDA0000986876760000023.GIF" wi="766" he="151" />其中,F<sub>x</sub>,F<sub>y</sub>,F<sub>z</sub>分别为三个轴方向产生的障碍规避力,K为增益系数,用于调节生成规避力的大小;步骤三,获得最小燃料的障碍规避轨迹;通过将步骤二利用导航函数生成的避障力F引入步骤一的动力学模型中,则推力发动机的推力被分为避障力F和控制力两部分,形成的改造后动力学方程如下:<img file="FDA0000986876760000024.GIF" wi="126" he="54" /><img file="FDA0000986876760000025.GIF" wi="126" he="62" /><img file="FDA0000986876760000026.GIF" wi="117" he="70" /><img file="FDA0000986876760000027.GIF" wi="499" he="118" /><img file="FDA0000986876760000028.GIF" wi="509" he="126" /><img file="FDA0000986876760000029.GIF" wi="603" he="119" /><img file="FDA00009868767600000210.GIF" wi="404" he="144" />令u<sub>F</sub>=[u<sub>Fx</sub> u<sub>Fy</sub> u<sub>Fz</sub>]<sup>T</sup>,u<sub>F</sub>为u<sub>F</sub>的大小,<img file="dest_path_FDA00005478182600000211.GIF" wi="835" he="146" /><img file="FDA00009868767600000212.GIF" wi="382" he="127" />则合力的控制向量为u<sub>net</sub>=u+u<sub>F</sub>;在实际着陆任务中,利用导航函数生成的避障力实际上也是由推力发动机产生,因而,合力的控制向量受如下约束:0<u<sub>min</sub>≤u<sub>net</sub>≤u<sub>max</sub>≤1其中,u<sub>net</sub>是合力u<sub>net</sub>的大小;采用高斯伪谱方法对改造的动力学方程进行优化,优化初始值为着陆器开始着陆时刻的初始状态,目标为着陆器目标着陆状态,优化性能指标为燃料消耗最小,优化结果即为燃料最省的着陆轨迹。
地址 100081 北京市海淀区中关村南大街5号北京理工大学