《数学建模嫦娥三号运行轨迹和着陆点分析》由会员分享，可在线阅读，更多相关《数学建模嫦娥三号运行轨迹和着陆点分析（24页珍藏版）》请在金锄头文库上搜索。
1、2021高教社杯全国大学生数学建模竞赛编 号 专 用 页赛区评阅编号由赛区组委会评阅前进行编号：赛区评阅记录可供赛区评阅时使用：评阅人评分备注全国统一编号由赛区组委会送交全国前编号：全国评阅编号由全国组委会评阅前进行编号：嫦娥三号软着陆轨道设计与控制策略摘要：根据题目附录和文献4中提供的嫦娥三号的运行参数，利用角动量守恒及向量几何的方法，分别确定了近日点、远日点的位置向量和速度向量。与文献4的真实数据比拟发现吻合良好。本文重点关注优化减速控制与着陆点避障两方面：前者燃耗最大，后者决定着陆成败。首先，在多重坐标变换根底上，建立了飞行器制动的动力学方程。并以燃耗为最优化性能指标、近月点状态为初始条
2、件、着陆点状态为终端条件，利用极值原理求解飞行器的着陆轨迹，及其最优控制参数。其次，对避障阶段采集的高程图采取水平剖分、比拟高程方差的方法，解出最优降落点。关键词：软着陆；最优轨道；避障241、 问题重述嫦娥三号于2021年12月2日1时30分成功发射，12月6日抵达月球轨道，于北京时间12月14号在月球外表实施软着陆。嫦娥三号在着陆准备轨道上的运行质量为2.4t，安装在其下部的主减速发动机能够产生1500N到7500N的可调节推力，其比冲即单位质量的推进剂产生的推力为2940m/s，可以满足调整速度的控制要求。嫦娥三号四周安装了姿态调整的发动机，在给定主减速发动机的推力方向后，能够自动通过多
3、个发动机的脉冲组合实现各种姿态的调整控制。嫦娥三号的预定着陆点为19.51W，44.12N，海拔为-2641m。嫦娥三号在高速飞行的情况下，为了保证嫦娥三号能准确地在月球预定区域内实现软着陆，关键的问题是着陆轨道与控制策略的设计。其着陆轨道设计的根本要求如下：着陆准备轨道为近月点15km，远月点100km的椭圆形轨道；着陆轨道为从近月点至着陆点，其软着陆过程共分为6个阶段，要求满足每个阶段在关键点所处的状态；尽量减少软着陆过程的燃料消耗。根据上述的根本要求，建立数学模型解决下面的问题：1计算其着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。2确定嫦娥三号的着陆轨道和在6个阶
4、段的最优控制策略。3对于设计的着陆轨道和控制策略进行相应的误差分析和敏感性分析。2、问题分析月球软着陆的方式有两种，分别是：1从月球转移轨道直至月球轨道，从而实现软着陆；2从月球的停泊轨道变轨到近月点，然后实现软着陆。2 嫦娥三号要求准确的在月球预定区域内实现软着陆，它降落的过程为如下几个阶段：嫦娥三号先进行霍曼变轨，由高度约100km的停泊轨道进入到近月点高度约15km的椭圆轨道；到达近月点时，探测器的制动发动机点火，进入动力下降阶段，即上图中的DA段；在距离月面约2km时，水平速度减为0，即上图中的AF段，调整姿态后，探测器以自由落体的方式降落到月面，即上图中的FL段。过程的主要考虑因素着
5、陆期间主要考虑两大因素：一、主减速段是用时最长、推进剂消耗最多，因此最小化推进剂消耗是该段制导律的主要设计目标。二、避障阶段决定飞行器能否顺利着陆，因此统计分析高程数据并采取相应横向推力最关重要。如上图示，以月心为坐标原点，月球自传轴为OZ轴，OX轴过零经度线，利用右手准那么确定OY轴，建立空间直角坐标系。嫦娥三号探测器绕月运动中，受到指向月心的引力F引力，以及相背于速度方向的由发动机产生的推力F推力两个力同时作用。引力是在月心坐标下考虑。F引= GMmr2rr其参数为：l G为万有引力常量且G = 6.6710-11l M为月球的质量， M=7.3491022kg。l m为航天器的质量l r
6、=（X,Y,Z）T是航天器质心的坐标，r就是航天器距月心的距离。发动机推力是在飞行器轨道坐标系下考虑的：F推=VedmdtFl Ve为发动机的冲力，是由发动机构造、性质、燃料共同决定。飞行过程中是不变的。Ve=2940ms=2.94m/sl m 为探测器中燃烧燃料的质量，l dmdt是秒耗量。表示飞行器燃料每秒燃烧消耗的质量。l F是推力方向，由飞行器的飞行姿态决定。F=1l 因此运动的加速度为：dr2dt2 =-GMmr2rr+F推m直角坐标系的旋转我们假设月球是一个不旋球体，不考虑侧向运动，我们建立的月球探测器的数学模型，利用燃耗最优原那么，依据庞特里亚金最大值原理设计了最优轨道。首先，我
7、们建立了以OXYZ为原点的月心惯性坐标系， OX轴是月球的0经度线； OZ是月球的自转方向；OY是根据右手坐标来确立的。要在月心惯性坐标下计算推力，需进行坐标旋转变换。a=a引+a推根据欧拉姿态角的变化过程，OXYZ经过两次有顺序的坐标轴旋转，得到两个矩阵：先将轨道坐标OX1Y1Z1沿OZ方向旋转角度，旋转变换的矩阵为：M=cos-sin0 sincos0001再将之沿绕OX方向旋转角度-，旋转变换的矩阵为：M-=0cossin0-sincos100因此得到总的旋转矩阵：T=M-M=0cossin0-sincos100cos-sin0 sincos0001=sincoscoscossin-si
8、nsin-coscoscoscos-sin0 AX1Y1Z1为原点在探测器的轨道坐标系, AZ1指向从月心到着陆器的延伸线方向，AY1垂直指向运动方向, AX1按右手坐标系确定。制动发动机推力F的方向与探测器纵轴重合，为F与AZ1轴正向所成夹角，为F在AX1Y1平面上的投影与AY1轴负向所成夹角。在月心惯性坐标中，推力为：F推=VedmdtFTl F是推力方向，如上图所示。l AZ1是以月心指向飞行器质心的向量。l 推力方向与AZ1的夹角。l 是推力在AX1Y1投影于AX1的夹角。 和随飞行器飞行状态而变，那么：F=sinsinsincoscos同理，在月心惯性坐标下，总的加速度为：a=Vem
9、dmdtTF-GMr2rr如上图所示，以月心为坐标原点，月球自转轴为竖轴建立空间直角坐标系OXYZ，假设近月点的空间直角坐标为(x,y,z)，用和分别表示近月点位置的经度和纬度，可由以下经纬度转换成为直角坐标系坐标的公式：Z=rsinY=rcoscos()X=rcossin1) 不考虑月球的自转。因为飞行器的速度相比月球自转的速度快、着陆时间快得多。2) 不考虑除月球以外其他天体引力对飞行器的扰动。3) 不考虑着陆区域附近因地形不同而造成的质量不同，进而影响引力场。即假设引力场在同样高度是均匀的。4) 假设月球引力场是从月心点发出的，引力指向月心点。5) 假设飞行器燃料足够使用。4、符号及字母
10、描述 5.1 近月点的参数根据文献资料数据 4，由北京工程研究所和空间智能控制技术国家级重点实验室的研究员们在这篇文献中对嫦娥三号的发射与着陆的设计做出了详细的说明,参与嫦娥三号设计与制造的他们在这篇文献中对嫦娥三号在登月过程中的各个阶段的各项参数都做出了详细而精确的推到与证明，通过这篇文献我们可以得到探测器近月点的位置信息为： 经度：19.0464W纬度：28.9989N高度：15km俯视姿态角：85资料中所给的近月点的速度计算近月点直角坐标由中构建的方程组，设点N为探测器在近月点处的位置，且该点所处在0度经线的位置上；近月点的位置为19.0464W，28.9989N，高度为月球半径+15k
11、m，着陆点的位置为19.51W，44.12N可得到近月点的直角坐标：XYZ=1719.5-182.6281，单位：km 5.2 远月点的参数用近月点计算远月点位置如上图所示，由于近月点、月心、远月点在同一直线上，近月点的位置为经度：19.0464W，纬度：28.9989N。运用投影的方法可得到远月点的位置，计算方法如下： 经度 = 180- 近月点经度 = 160.9536W 纬度 = 28.9989S远月点的直角坐标由中的方程组结合远月点的经纬度数据可以得到远月点的直角坐标为X1Y1Z1=1135.41413.8294.7,单位km5.3.1推导、计算月球的第一宇宙速率根据万有引力定律： G
12、Mmr2=mv2r v=GMr (月球质量 M=7.349*1022kg，G = 6.67*10-11106m )由此可以得到月球的第一宇宙速度为：V1。由于远月点处于嫦娥三号探测卫星的变轨点也就是变轨后的轨道与环月运动轨道的交点处，因此，远月点的速度即为探测器绕月运行的速度，即月球的第一宇宙速度，既得远月点的速度为V1。首先由角动量的定义可知：L=rp=m(rv)其中L为角动量，r为运动轨迹向量，p为动量，v为速度向量，m为探测器的重量。然后，刚根据角动量守恒定律：rv=常量。，远日点初次变轨时的位置的路径长度： r1=100km+月球半径r。由远日点与近日点角动量守恒得：r1v1=r2v2
13、 v2=r1v1r2代入数据可得 v2。与资料中所给的 v2=1.6957km/s相比拟，计算所得的数据与资料所给数据产生的误差3.8% 5%，因此计算所得的近日点速度是吻合且正确的。由5.1.1中图所示，设点N为探测器在近月点处的位置，且该点所处在0度经线的位置上；点L为探测器着陆点的位置，由此构建方向向量OL，ON以及探测器的速度方向向量V，近月点的位置为19.0464W，28.9989N，高度为月球半径r+15km，着陆点的位置为19.51W，44.12N根据向量的原理整理总结可得如下方程组：ONV=0ONON+V=ONOL 将ON，OL数据代入结合5.1.1中所计算出的近月点N的空进直
14、角坐标，运用matlab进行运算后整理可得如下线性方程组：1719.5-182.62810-281-182.62810-1719.5XYZ=-128210-254050619390综上所述及计算我们可以得到近月点的速度方向向量为：V=XYZ=145.31616.1-492.66、主减速阶段的轨道方程根据2.3中牛顿力学分析；经过2.4.1坐标方程变换后，飞行器的动力学方程如下：axayaz=TQVemsinsinQVemsincosQVemcos-gxgygz(其中Q=dmdt表示秒耗量)整理6.1中的方程，将其翻开引入简记符号后，系统状态方程可表示为：dXdt=VXdYdt=VYdZdt=VZdVXdt=aQVem+gXdVYdt=bQVem+gYdVZdt=cQVem+gZdmdt=-Q其中，a=sincossinsin+coscossincos+sincos b=-sinsinsinsin-coscossincos+coscos c=cossinsin-sinsincos
