月球软着陆控制系统综合仿真简析.doc

月球软着陆控制系统综合仿真及分析 姓名：冯思松 学号：2120120040 学院：宇航学院 专业：飞行器总体设计 引言： 在月球探测带来巨大利益的驱使下，世界各国纷纷出台了自己的探月计划，再一次掀起了新一轮探月高潮。在月球上着陆分为两种，一种称为硬着陆，顾名思义，就是探测器在接近月球时不利用制动发动机减速而直接撞击月球。另一种称为软着陆，这种着陆方式要求探测器在距月面一定高度时开启制动系统，把探测器的速度抵消至零，然后利用小推力发动机把探测器对月速度控制在很小的范围内，从而使其在着陆时的速度具有几米每秒的数量级。显然，对于科学研究，对探测器实施月球软着陆的科学价值要大于硬着陆。 研究方法： 本实验要求对一种控制方式下的月球软着陆过程进行仿真，并分析相关参数变动对软着陆过程的影响。需要利用Matlab/Simulink建立月球软着陆制导控制系统的仿真模型。需要变动的相关参数包括：发动机推力、比冲、初始速度方向。 理论根据： 月球软着陆的动力学方程与两个坐标系有关。月心惯性坐标系定义为 ：原点 选在月心，轴指向动力下降起始点，轴垂直于轴指向着陆点方向，轴按右手法则确定。表示月心到探测器的距离，和表示月球经度和纬度。探测器轨道坐标系定义为：原点为探测器质心，轴与从月心到探测器质心的矢径方向重合，背离月心方向为正，轴垂直于轴指向运动方向为正， 按右手法则确定。制动推力F 的方向与探测器本体轴重合，和为在轨道坐标系中表示的推力方向角。假设制动发动机为常推力液体发动机，忽略月球自转，则月球软着陆动力学方程可表示为： 式中：，为发动机比冲，为地表重力加速度常数。 月球软着陆器的制导控制方案，采用一种燃耗次优制导率。其作用是求解推力方向角和。可以将公式表示为 式中：，，分别为，，的终值。 通过制导率求解出的推力方位角带入动力学方程可求出软着陆器的运动轨迹。 系统仿真： 推力器的参数为：，，； 月球常数：，月球半径； 初始参数：，，，，，； 终端参数：，，，。 对于初始参数，取0会导致公式分母为0，故在实际计算中定义。 通过Simulink搭建模型如下： 图 41 Simulink模型 图 42 动力学Dynamic子系统 图 43 制导控制率Control子系统 仿真分析： 通过系统仿真，可得到速度、位置、推力方位角等参数随时间的变化曲线。如图 51所示为着陆器到月心距离随时间变化曲线。着陆器下降到具月球表面2km高度用时534s，推进剂燃烧后飞船质量下面将分别针对发动机推力、比冲和初始速度方向的偏差，分析其对着陆器飞行过程的影响。 图 51 到月心距离随时间变化曲线 发动机推力偏差对软着陆过程的影响 发动机推力偏差为±10%，标准推力为1500N，最小推力为1350N，最大推力为1500N。由于制导控制率不变，着陆器仍然能下降到具月球表面2km的高度，但时间会缩短。由图 52可见推力的增大可以缩短着陆器下降的时间。 由于比冲不变，推力变化会引起燃料消耗速度的变化，如图 53所示。然而，由于飞行时间的减小，大推力下总的燃料消耗量会减小。可见采用大推力发动机可以减小燃料的使用量。但实际情况下还要考虑大推力发动机是否会增加额外的重量，因为推力增大对燃料的节约很有效。如这个系统中，增大10%的推力只节约了1%的燃料。 图 52 到月心距离虽时间变化曲线 图 53 着陆器质量变化曲线 飞行时间的不同会导致最终下降的目的地不同，如图 54和图 55所示。但这只会导致纬度的不同，而经度不受影响。而且在整个下降过程中着陆器所处的经度都没有变化。 图 54 纬度变化曲线 图 55 经度变化曲线 推力的变化会影响减速后的末速度。虽然任务要求下降段结束后的末速度为零，但实际情况仍然具有很低的末速度。由图 56和图 58可知，大推力会使减速后的末速度减小，但缺点是会增大下降时的过载。从图 58可以明显看出下降过程中垂直方向分速度的变化趋势。在减速初期，着陆器会加速下降，末期下降速度会逐渐减小。但无论哪个过程，大推力下加速度都会更大，这也说明了大推力工作时下降过程时间更短的原因。由于y轴的分速度很小，所以不做考虑。 图 56 y轴方向分速度变化曲线 图 57 z轴方向分速度变化曲线 图 58 x轴方向分速度变化曲线 比冲偏差对软着陆过程的影响 发动机比冲偏差为±10%，标准比冲为300，最小比冲为270，最大比冲为330。比冲偏差对着陆过程的影响要小于推力偏差的影响。在相同的推力下，推进剂的比冲越大会延长下降段的时间。这会需要发动机工作更长时间。但由于当推力一定时，大比冲的推进剂单位时间的消耗量更