2019数学建模a题

问题一分析。

对于问题一，需确定着陆准备轨道近月点和远月点的位置，以及嫦娥三号相应速度的大小与方向。文章中给定了着陆点坐标，该段中，将主要减速阶段（即从15km到3km）作为制动段。制动终点到达着陆点。

采用逆推的方法，通过着陆点坐标和建立的动力学方程，逆推出近月点的坐标以及速度的大小方向。然后通过角动量守恒定理，算出远月点的速度。着陆器距离月面相对较高，且着陆器走过的月面距离比较长，将月球视为平面建立模型会带来较大的偏差。

因此，制动段有必要将月球视为球体来建立均匀球体下的三维软着陆模型。制动段推进系统采用常值推力方式，通过姿态控制来完成制动力方向的改变。

3.1.2均匀球体三维动力学模型。

首先定义几个坐标系：

参考惯性坐标系oxryrzr原点o位于月球中心，zr轴由月心指向初始软着陆点，xr轴位于环月轨道平面内且指向前进方向，yr轴与构成直角坐标系。该坐标系仅用于软着陆下降轨迹和制导律设计中；

2)下降轨道参考坐标系ox0y0z0.原点o位于着陆器质心，zo轴由月心指向着陆器质心为正，xo轴位于当地水平面内且指向着陆器前进方向，yo轴与xo和zo轴构成直角坐标系；

3)着陆器oxbybzb体坐标系。原点o位于着陆器质心，xb轴在制动推力矢量延长线上，沿推力方向为正，yb,zb轴分别根据着陆器上仪器设备的安装而定，并与xb轴构成直角坐标系。坐标系示意图及着陆器位置与推力矢量关系如图2

所示：图1软着陆坐标系定义与推力矢量空间关系

图。2(a)

给出了上面各坐标系的示意和着陆器在坐标系中的位置，图。

2(b)给出了f在下降轨道参考坐标系中的位置。

其中，α为在xryr平面内的横向月心角；β为下降轨道平面内的纵向月心角；推力f与坐标系ox0y0z0之间的2个推力方向角分别为推力方位角ψ和推力仰角θ,他们定义为：

推力方位角绕正zo轴旋转为正，推力仰角绕负yo轴旋转为正。

分别用u,v,w表示着陆器下降速度在坐标系ox0y0z0三轴上的分量，于是有。

若不考虑摄动影响且忽略月球自转，同时引入质量方程，可利用球坐标系与直角坐标系的关系最终得到下降轨道参考坐标系下的软着陆动力学模型

下降轨道上分析可得α是个定值，则v应始终为0，则软着力动力方程可简化为。

运用matlab

软件对上述方程组（3）进行求解，结果如下，m=2400-(f*t)/(isp*g)

u=(pi*isp*t)/900-(pi*isp)/2+(isp^2*m*pi*log(f*t-isp*m))/900*f)-(isp^2*m*pi*log(450*f-isp*m))/900*f) w=isp/2-(isp*t)/450+(225*um)/ryue^2-(t*um)/ryue^2+(isp^2*m*log(225*f-isp*m))/450*f)-(isp^2*m*log(f*t-isp*m))/450*f)+57/2r= int(((pi*isp)/2-(pi*isp*t)/900-(isp^2*m*pi*log(f*t-isp*m))/900*f)+(isp^2*m*pi*log(450*f-isp*m))/900*f))/isp*t^2)/900-t*((isp^2*m*log(-isp*m))/450*f)+1700)+(t^2*um)/(2*ryue^2)-(isp^3*m^2*log(f*t-isp*m))/450*f^2)-(isp^2*m*t)/(450*f)+(isp^3*m^2*log(-isp*m))/450*2)+(isp^2*m*t*log(f*t-isp*m))/450*f)-15000),t)

注：：dbeta为β，deta为β）

程序见附件程序。

通过以上求解，我们发现完全求解软着陆的动力学方程组。

有一定难度，并且解得的解析。

解不利于数值的求取和问题分析，主要表现在两方面。

求出的解无法进行积分，导致后面。

不能求出某一时刻的速度和位置。

积分求出的解，带入不同的初始值，解的差异很大，且。

都不符合事实依据。但是以上求解的切向分速度。

u的结果比较符合事实和推理，且数值准。

确，容易运算分析，其图像如下图（

所示，因此我们将动力学方法进行简化改进得方程。组（

进而求解。附件程序1clear

isp=2940;

m=dsolve('dm=-f/isp/g','m(0)=2400')

u=dsolve('du=-f/(m-(f*t)/isp)*(pi/900*t)',u(450)=0')

w=dsolve('dw=f/450*t*isp/(m*isp-f*t)-um/ryue^2','w(225)=28.5')

r=dsolve('d2r=f/450*t*isp/(m*isp-f*t)-um/ryue^2','r(0)=15000,dr(0)=17000')

b=u/rbeta=int(b,t)

2019数学建模a题

2019数学建模A题

2019数学建模A题

2019数学建模B题

其他用户还读了