学文网1月地转移轨道设计约束条件分析
月心段返回轨道相对月心为双曲线轨道，但为了节省燃料，地心段轨道一般采用地心椭圆轨道。如***1所示。本文所使用的边界条件有：1）目标矢量rT。rT是从地球质心到月球影响球（SphereofInfluence）上一点的矢量，代表航天器出现在月球影响球上的位置。因此，该矢量与航天器穿过月球影响球的历元TSoI相关。2）再入点矢量ri。ri是地心到再入点的矢量。忽略地球半径在两极和赤道之间大约20km的差别及其他差异，可近似认为再入高度为一常数，本文将再入高度取为120km。3）返回轨道相对地球赤道面的倾角iR。从rT到ri的轨道有2条，θ＜π或θ＞π，对应***1中的轨道1和轨道2，倾角分别为iR和180°－iR。由于逆行返回轨道相对大气的速度会增大，而顺行返回轨道则相反，因此逆行返回轨道不适合工程应用。4）再入面的飞行速度倾角γi。γi（即再入角）总是为负值，一般取绝对值。设计中要求其在一个范围（即再入走廊）内，γimin＜γi＜γimax，其中γimin为允许的最小再入角，γimax为允许的最大再入角。5）出月球影响球方向。必须选择是降轨出月球影响球还是升轨出月球影响球，否则将有2个不同的轨道面包含目标矢量。6）飞行圈数。从目标矢量到再入矢量可以单圈到达，也可以飞行N圈（N＞0）后到达。飞行圈数大于0的情况在航天器交会对接中经常使用。在月地转移轨道中，为了尽快返回，只考虑单圈返回的情况，即取N＝0。7）从月球影响球到再入的飞行时间TFR。为了尽快返回地球，脉冲转移方式下，月地转移时间一般为3～5天。如果忽略月球重力场，并将月地转移轨道近似为半长轴等于地月平均距离（385000km）的椭圆轨道，那么该轨道半个周期约为5天（～120h），该值可以看作是月地转移飞行时间上限的近似值［11］。研究结果表明，当月球离地球最近时，从月球影响球开始的最大飞行时间大约为3．5天；而当月球离地球最远时，从月球影响球开始的最大飞行时间大约为4天。抛物线飞行时间可以作为TFR的下限。可以用Barkers方程得到2个矢径（rT＜ri）间抛物线轨道飞行时间的表达式，见式（1）。该方程假设rT和ri定义的点之间没有近拱点。根据经验，对月球转移轨道而言，当月球离地球最近时TFR大约为1．7天，最远时大约为2．1天。
2摄动项分析和轨道动力学模型
2．1摄动力及其量级对低轨（高度为100～300km）月球卫星来说，所受的摄动源包含下列10类，对应的摄动量级见表1［12］。对高轨地球卫星，比如GEO（GeostationaryEarthOrbit，地球静止轨道）卫星而言，所受的摄动力及其量级见表2。从上表可以看出，对于一般的轨道分析，只要考虑非球形引力摄动和第三体引力摄动，这2类摄动因素的影响可体现航天器轨道变化的主要特征。2．2轨道动力学方程根据摄动力及其量级分析结果和月地转移轨道的精度要求可以看出，在实际工作中，对于一个具体问题，并不需要考虑所有的摄动力因素，而是可以对具体的摄动力进行取舍，采用合适的有限精度模型。下面以月球影响球为界，采用2个受摄二体问题建立月地转移轨道的轨道力学运动方程。当月球探测器运动于月球影响球中时，月球引力作为主要力源，地球引力和太阳引力作为摄动力，根据摄动力的大小，月心段轨道运动考虑的摄动因素主要为地球引力摄动和太阳引力摄动，运动方程为¨r＝－μmr3r－μe（redr3ed＋rer3e）－μs（rsdr3sd＋rsr3s）其中μe，μm和μs分别是地球、月球和太阳的引力常数；re，rs和r分别为地球、太阳和月球探测器在月心天球坐标系中的位置矢量；red表示探测器到地球的矢径；rsd表示探测器到太阳的矢径；地球和太阳的星历从DE405中获取。当月球探测器出月球影响球后，地球引力成为主要力源，月球和太阳的引力为摄动力。根据摄动力的大小，地心段轨道考虑的摄动因素有地球非球形引力摄动、月球引力摄动和太阳引力摄动，建立如下运动方程¨r＝－μer3r－μm（rmdr3md＋rmr3m）－μs（rsdr3sd＋rsr3s）＋ae其中rm，rs和r分别为月球、太阳和卫星相对于地心的矢径；rmd表示探测器到月球的矢径；ae为地球非球形摄动加速度；月球和太阳的星历从DE405中获取。
3月地转移轨道精确设计方法
3．1微分改正法以基于双二体模型得到的出月球影响球时刻和相应的位置、速度为初值，通过迭代修正出月球影响球的速度，使从月球影响球出发积分到再入时刻的位置与指定再入点的位置偏差满足精度要求。该迭代过程中，初始参数为出月球影响球的速度矢量，终端参数为再入点的位置矢量。终端参数可以表示成初始参数的函数，记为q＝f（p）其中q为终端参数，q＝［qxqyqz］T为再入点的位置；p为初始参数，p＝V＝［vxvyvz］T为出月球影响球的速度。按线性摄动法，初始参数的小偏差Δp传播至再入时刻带来的终端参数的变化Δq可近似表示为Δq＝fpΔp（2）记误差传递矩阵为K，K＝fp为3×3的矩阵，也可以表示为3个偏导数矩阵的乘积K＝fp＝qx1x1x0x0p其中x0和x1分别为探测器的初始和终端状态，即位置和速度。利用以上关系，一个初始参数可以得到对应的唯一一个终端参数。但对转移轨道来说，根本无法找出K的解析表达式，而我们的目的主要是得到它的具体数值。因此本文采用有限差分法进行求解，其基本原理就是在初始参数的每个分量上增加一个小量，计算终端参数qn。计算偏差Δpn＝qn－qtarget其中qtarget为要求的目标终端参数；n为迭代次数。解式（2），得到出月球影响球速度的修正量ΔVΔp＝ΔV＝K－1Δq然后更新初始参数pn＋1＝pn＋Δp重复该过程，直到ΔV为零，得到最终解。3．2精确轨道设计的计算流程本文采用了一种双向嵌套循环搜索算法，求解同时满足两端约束条件的精确月地转移轨道。以出月球影响球的时刻和位置、速度为中间变量，一方面采用前向数值积分和微分改正法搜索满足地球再入端的轨道，另一方面采用后向数值积分并进行倾角和近月距修正得到满足月球端的轨道。由于月地转移轨道与月球影响球交点的位置初值是在无摄动的情况下得到的，这样，一旦用受摄模型反向积分到近月点时，得到的近月距和倾角将不是原先设定的月球停泊轨道近月距和倾角，所以必须对中间变量进行倾角和近月距修正。一般修正2～3次就可满足容差要求。通过这种双向嵌套循环搜索算法，可以求解同时满足两端约束条件的精确轨道，并使得两段轨道在月球影响球边界处的位置和速度连续，从而获得一条完整的满足两端约束条件的月地转移精确轨道。月地转移轨道精确设计的计算流程如***2所示。
4设计算例
4．1约束条件和参数设置选取2017年1月26日出月球影响球作为返回窗口，下面给出该窗口下月地转移轨道精确设计的一个算例。表3为该算例对应的约束条件设置。根据表3给出的约束条件，以快速轨道设计结果作为计算初值，采用精确轨道动力学模型和数值积分法，通过3．2节的设计流程，最后可以得到4条满足两端约束条件的月地转移轨道，分别为：降轨出月球影响球顺行轨道、降轨出月球影响球逆行轨道、升轨出月球影响球顺行轨道和升轨出月球影响球逆行轨道。由于逆行轨道返回所需的速度增量远大于顺行轨道，速度增量差可以达到113m／s，因此在工程实际中不考虑逆行返回。4．2设计结果顺行的月地转移轨道入轨所需的速度增量小，适合工程使用，这里仅给出其中一条升轨出月球影响球的顺行精确轨道的最后计算结果，见表4。在STK中采用同样的模型对该轨道进行仿真，结果见***3～***6。***3给出了从100m高的环月停泊轨道经一次制动进入月地转移轨道的惯性轨迹，***4为该轨道在3D下的全景，***5和***6分别给出了地面和月面上的星下点轨迹。可以看出，STK的仿真结果与本文程序计算结果完全一致。5结束语本文以我国后续月球探测任务的实际工程需求为背景，以基于Lambert算法的快速轨道设计结果为初值，开展了精确轨道设计研究。通过对整个返回飞行阶段的摄动项和量级分析，建立了适合于月地转移轨道精确设计的轨道动力学方程。精确轨道设计以出月球影响球的时刻和位置、速度为中间变量，将轨道分为地心段和月心段分别进行计算。采用微分改正法和前向、后向数值积分的双向嵌套循环搜索算法，求解同时满足两端约束条件的精确轨道。通过这种双向嵌套循环搜索算法，获得分别满足月球端约束条件的月心段轨道和满足再入端约束条件的地心段轨道，并使得两段轨道在月球影响球边界处的位置和速度连续，从而获得一条完整的满足两端约束条件的月地转移精确轨道。最后以2017年1月26日出月球影响球作为返回窗口，给出了具体的设计算例，并将程序计算结果在STK中进行了仿真验证，进一步证明了其正确性。
作者：郑爱武 周建平 单位：航天飞行动力学技术重点实验室 北京航天飞行控制中心 中国载人航天工程办公室