《数学建模获奖论文材料-嫦娥三号软着陆轨道设计与控制策略》由会员分享，可在线阅读，更多相关《数学建模获奖论文材料-嫦娥三号软着陆轨道设计与控制策略（32页珍藏版）》请在金锄头文库上搜索。
1、.2014高教社杯全国大学生数学建模竞赛承 诺 书我们仔细阅读了全国大学生数学建模竞赛章程和全国大学生数学建模竞赛参赛规则（以下简称为“竞赛章程和参赛规则”，可从全国大学生数学建模竞赛网站下载）。我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛章程和参赛规则的，如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我们郑重承诺，严格遵守竞赛章程和参赛规则，以保证竞赛的公正、公平性。如有违反竞赛章程和参
2、赛规则的行为，我们将受到严肃处理。我们授权全国大学生数学建模竞赛组委会，可将我们的论文以任何形式进行公开展示（包括进行网上公示，在书籍、期刊和其他媒体进行正式或非正式发表等）。我们参赛选择的题号是（从A/B/C/D中选择一项填写）： A 我们的报名参赛队号为（8位数字组成的编号）： 所属学校（请填写完整的全名）： 华南理工大学 参赛队员 (打印并签名) ：1. 2. 3. 指导教师或指导教师组负责人 (打印并签名)： （论文纸质版与电子版中的以上信息必须一致，只是电子版中无需签名。以上内容请仔细核对，提交后将不再允许做任何修改。如填写错误，论文可能被取消评奖资格。） 日期：2014 年 9 月
3、 15 日赛区评阅编号（由赛区组委会评阅前进行编号）：2014高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）：嫦娥三号软着陆轨道设计与控制策略摘 要嫦娥三号的核心任务是实施高可靠高安全的月面软着陆, 要求着陆器必需具备自主障碍识别与规避能力。本文将通过研究月球软着陆主减速段、快速调整段、避障段和缓速下降段的飞行动力学模型，在动力学方程的基础上对各阶段制导律进行优化设计后，结合工程实际确立嫦娥三号软着陆轨道的下降轨迹
4、和控制策略。针对第一问，本文基于开普勒定律求解着陆准备轨道的运动方程，得到近月点与远月点的速度。由于着陆点的位置是由近月点的位置和软着陆轨迹共同决定的。因此本文通过求解嫦娥三号在主减速段和粗避障段的经过的横向距离，利用着陆点的位置近似确定近月点的位置。实验结果表明，此近似求解的近月点位置误差在10%以内。同时，本文所建立的基于飞行器动力学模型的主减速段燃料最优制导律，通过蚁群算法的参数优化，也取得误差5%以内的良好结果。本文还将证明主减速段的燃料最优制导律等效于时间最优制导律。针对第二问，本文通过查阅相关资料，给出了软着陆各个阶段的动力学方程和制导律。主减速段是耗时最长，燃料消耗最多的阶段，因
5、此采用了燃料最优制导律。快速调整阶段是过渡阶段，拟采用推力线性变化制导律。粗避障段的主要目的是避开大的陨石坑，由于着陆器此时接近垂直下降，因而采用最优开关制导律。为实现粗避障，本文采用对光学图像进行分格，并且判断各格为障碍区（0）或非障碍区（1），最终通过在0-1图上寻找具有安全半径的连片非障碍区，找到合适的着陆区域，取得了良好的实验效果。精避障段则要构建三维高程图和平均坡面，计算平均坡度，估算出障碍高度，最终找出安全着陆区。缓速下降段主要考虑的是着陆的安全性，本文将利用变推力制导律和垂直软着陆模型进行分析。针对第三问，本文将在模型的分析与检验中给出相应的误差分析，实验结果表明，本文所提出的模
6、型与嫦娥三号实际结果相符合，具有实用性和准确性。关键词：制导律 蚁群算法 数字高程图 嫦娥三号 软着陆1.问题的重述1.1 背景实施在月球表面的软着陆是月球勘探计划的重要一步。所谓月球软着陆,是指着陆器在制动系统作用下以很小的速度准确降落到月面指定区域,以保证试验设备和宇航员的安全。由于月球表面没有大气，因此整个软着陆过程需要制导发动机的控制。正确地建立着陆器的飞行动力学模型，并在此基础对着陆器的下降轨迹和制导律进行优化设计，对于着陆器飞行程序的设计和燃料资源消耗等方面具有重要的指导意义。1.2 需要解决的问题根据题目信息和附件所示内容，分析嫦娥三号软着陆过程中的6个阶段，研究以下的问题：（1
7、）根据附件1与附件2研究嫦娥三号的着陆准备轨道，确定该椭圆轨道上近月点与远月点的位置，以及嫦娥三号相应的速度大小与方向。（2）建立合理的着陆器飞行动力学方程，根据软着陆过程中6个阶段不同的特点与要求，结合附件1与附件2的内容，确立嫦娥三号的下降轨迹和各个阶段的最优控制策略。（3）根据问题一和问题二建立的数学模型，对模型进行相应的误差分析和敏感性分析。2.模型假设（1）假设诸如月球引力非球项、日月引力摄动等影响因素均可忽略不计（2）月球表面附近没有大气，所以在飞行器的动力学模型中没有大气阻力项（3）嫦娥三号在软着陆全过程不发生意外碰撞和震荡3.符号说明符号符号说明嫦娥三号进行霍曼转移的初始时刻嫦
8、娥三号到达近月点的时刻嫦娥三号在主减速段的终止时刻着陆器的月心距月球半径远月点与月面的距离着陆器在近月点的月心距近月点与月面的月心距月球的引力常数在时刻的横向速度在时刻的径向速度在时刻的横向速度在时刻的径向速度横向速度径向速度着陆器的月心距着陆器的极角推力加速度发动机的推力大小着陆器的初始质量推力方向角的系数，待优化参量燃料消耗率燃料消耗质量比冲终止时刻合速度城市和的距离t时刻蚂蚁k由城市转移到城市的概率信息量表示蚂蚁k下一步允许走过的城市集合信息量对蚂蚁的影响力大小由城市i转移到城市j的期望程度终止时刻，待优化的时间参量的作用信息增量接近段的水平距离主减速段的水平距离近月点与着陆点的距离粗避
9、障段的初始时刻粗避障段的终止时刻月面高度下降速度与垂直方向的夹角，称为下降角月球的重力加速度制动力的开关量等效推力加速度的斜率近月点与着陆点的水平距离相对误差近月点速度相对误差主减速仰幅角的相对误差主减速段时间相对误差单格内障碍点数单格内障碍点阈值过暗灰度临界值过亮灰度临界值注：未列出的符号及重复的符号以出现处为准。4.问题一的建模与求解4.1 问题分析根据题目信息，着陆准备轨道为近月点15km，远月点100km的椭圆轨道，根据开普勒定律，可以建立椭圆轨道上的运动方程，求解出近月点与远月点的速度的大小与方向。着陆点的位置是由近月点在月心坐标系中的位置和软着陆形态共同决定的。因此可以通过求解软着
10、陆的下降轨迹，结合着陆点的位置，确定近月点在月心坐标系中的切向距离。此外，利用月球的半径参数，可以求得近月点的径向距离，进而确定近月点的位置。相应地，远月点的位置可以结合近月点的位置，利用椭圆的轨迹方程进行确定。在本问题中，软着陆的下降轨迹的求解是一个难点，考虑到着陆器在近月点处离月面距离较高，因此可以在不引起较大的误差的情况下，建立简单的模型进行求解。根据相关参考资料1有：在接近段中，着陆器采用接近直线下降方式逐步接近着陆区，利用附件1与2提供的接近段的高度信息，可以求解出对应的水平距离。而主减速段的轨迹可以建立飞行动力学模型进行求解。考虑到主减速段是软着陆过程用时最长、推进剂消耗最多的任务
11、段。该段的主要任务是消除较大的动力下降段初始水平速度(约 1.7 km/s), 因此推进剂消耗优化是该段制导律的主要设计目标。因此，求解基于动力学模型的燃料消耗最低制导律方程，可以确定主减速段的水平距离。综合上述求解出的距离信息，就可以在一定的误差允许范围内给出近月点的位置，进而得到远月点的位置。总体而言，问题一的流程图如下所示：图4-1问题一的流程图4.2 模型的建立4.2.1 着陆准备轨道的运动学方程假设嫦娥三号在时刻开始进行霍曼转移，在时刻到达近月点，则在惯性坐标系中，以月心为原点的极坐标形式的着陆准备轨道的运动学方程2为：在时刻：（4-1）其中的为着陆器的月心距，为月球半径，可取值为月
12、球的平均半径1737.013km，为远月点与月面的距离，取值为100km，为着陆器在近月点的月心距，为近月点与月面的距离，取值为15km，为月球的引力常数，其值为：，为初始横向速度，即着陆器从远月点开始运动的速度，为初始径向速度，其值为0.在时刻：（4-2）其中的为初始横向速度，即着陆器从近月点开始运动的速度，为初始径向速度，其值为0.4.2.2 主减速段的动力学方程在惯性坐标系中，以月心为原点的极坐标形式受控飞行器动力学方程为2：（4-3）其中为月球引力常数，为径向速度，为横向速度，为着陆器的月心距，为着陆器的极角。其中为推力加速度，其值可以由下式计算：（4-4）其中为发动机推力的大小，查阅
13、相关资料1可得：主减速段的主要任务是消除较大的动力下降段初始水平速度，所以发动机一直工作在最大推力模式即7500N。为着陆器的初始质量，其值为2400kg，为比冲（即单位质量的推进剂产生的推力），其值为2940m/s。为推力方向角，即推力方向与当地水平线的夹角，查阅相关资料3可得：（4-5）其中,为待优化的参量。4.2.3 基于蚁群算法3的主减速段燃料最优制导律4.2.3.1 燃料最优制导律设主减速段的初始时刻为，终止时刻为，对于推力幅值恒定的飞行器，燃料消耗最省的性能指标为：（4-6）其中为燃料消耗率，为了使消耗的燃料质量最小，由于初始时刻，发动机推力和比冲均为定值，因此只能使主减速段的终止时刻最小。所以，对于主减速段而言，燃料最优制导律等同于时间最优制导律。查阅相关资料1可得，主减速末期的俯仰角为，因此可列出如下等式：（4-7）其中为终止时刻的合速度，为57m/s.因此，只要确定了的函数
