哈工大深空探测轨道设计课程作业..docxVIP

地球—火星转移轨道设计1.研究问题描述在未来3年内寻找发射机会，设计地球－火星转移轨道，以总能量最小为指标。假设地球停泊轨道为高度200km、轨道倾角28.5o的圆轨道，目标轨道为高度500km、倾角90 o的火星绕飞轨道。轨道动力学模型只考虑太阳、地球、火星的引力作用及地球J2项摄动。分别给出初步轨道设计参数和精确轨道设计参数。2.解决方案首先，基于pork-chop图寻找合适的发射窗口，然后采用圆锥曲线拼接法，进行轨道的初步设计，分段用二体模型设计各阶段轨道，最后在精确力学模型下，采用微分修正法，基于B平面参数进行精确轨道设计。2.1地火转移轨道的初步设计对于初步设计，是以圆锥曲线拼接法(Patched Conic Method)为基本原理，将探测轨道划分为几个分段过程，每个分段过程可以近似成二体问题处理，然后通过接口，将每一段拼接在一起。它是在二体模型假设下，通过求解Lambert问题确定发射窗口，利用圆锥曲线拼接法确定地心段、日心段、火星段的初始轨道参数。2.1.1日心轨道设计及发射窗口的搜索本文采用等高线图法搜索发射窗口，通过给出既定时间段内所有的发射和到达时间情况，获得初始和终了位置，进而解算Lambert问题，从而绘制“猪排”图，观察得到最优发射机会。基于等高线图的最优发射机会搜索算法的主要步骤如下：根据任务的需要确定出发射机会搜索的目标函数（性能指标）、发射时间的区间以及飞行时间的区间；选取一组出发时刻t0和达到时刻tf，根据行星历表计算地球的位置RE(t0)、速度VE(t0)和火星的位置RM(tf)、速度VM(tf)；利用转移时间（tf- t0）、RE(t0)、以及RM(tf)，通过求解Lambert问题，可以得到探测器在始末位置处的速度矢量V1(t0), V2(tf)；确定发射机会的目标函数，并绘制出等高线图；根据等高线图，找到目标函数取值较小的区域，确定出性能指标指标最优的发射时间；以上步骤中涉及到的目标函数（性能指标）通常是指双曲线超速v∞，发射能量C3，速度增量Δv，它们的具体定义为：双曲线在达到地球引力影响球边缘时速度有剩余，这个双曲线剩余速度v∞通常称为双曲线超速，计算公式为式中v1是飞行器的速度矢量，vLE是发射时刻地球绕太阳公转的速度矢量。发射能量是影响任务初始设计的关键参数，在飞行器质量一定的情况下，发射能量越大，所需运载火箭的运载能力越强，它的大小是发射时双曲线超速的平方当假设飞行器从停泊轨道开始转移，这在停泊轨道上施加的速度增量Δv1为类似地，可以定义飞行器到达目标天体的双曲线超速以及到达目标停泊轨道，交会时的速度增量式中v2是飞行器到达目标行星时的速度矢量，vAT是发射时刻目标行星绕太阳公转的速度矢量，rTp为飞行器在目标天体俘获是近心点的高度。在整个飞行过程中，所需的总的速度增量Δvtotal为2.1.2地心段参数的确定对地心段参数进行初步设计时，假设在地球影响球内探测器只受到地球的引力作用，从而按照二体轨道特性进行轨道参数的计算，逃逸速度增量假设为脉冲。假设探测器的发射过程为:探测器从地面发射后，首先进入圆形的停泊轨道;然后在经历一段时间的滑行后，在特定的时间，特定的位置，经过末级火箭加速后进入双曲线轨道进行逃逸，如图1所示。为了尽量避免不必要的燃料消耗，假定停泊轨道和双曲线轨道是共面的，且双曲线的近地点半径与圆轨道半径相同，速度相切。因此，地心段参数初步设计的内容主要包括两个方面:停泊轨道参数的特性分析以及停泊轨道参数选择。图1地心逃逸双曲线轨道面内参数的几何关系设停泊轨道半径(即双曲线近地点半径)为今1，根据双曲线轨道参数的相关公式容易算得双曲线轨道半长轴、偏心率、动量矩以及近地点速度:其中a1为轨道半长轴，e1为偏心率，为地心引力常量，v∞1为地心双曲线剩余速度大小，hl为双曲线轨道的动量矩大小，vp1为双曲线近地点速度大小。而根据剩余速度矢量v∞1可以计算出转移双曲线渐近线对应的赤经α∞1（RLA: Right ascension）和赤纬δ∞1（DLA: Declination）进一步可以确定B平面中的渐近线的单位方向矢量而B平面的其他两个坐标系可以表示为式中，B平面的夹角可以由渐近线的赤纬和轨道倾角确定双曲线的单位角动量矢量，可以表示为当速度趋于无穷时，真近点角的正弦和余弦可以确定为进一步可以确定转移双曲线轨道的近地点矢量从而可以确定航天器在双曲线转移轨道的近地点处的位置矢量为，而速度的单位矢量为而航天器在近地点的速度矢量可以表示为。根据近地点的位置矢量和速度矢量，就可以计算出双曲线轨道的轨道六根数。考虑到行星际距离比地球影响球大很多，故地心段双曲线轨道与地球影响球的交点位置对轨道几乎没有影响，而v∞1矢量的大小和方向完全决定了探测器在脱离地球作用范围之外的运动，因此控制