专利名称：一种考虑目标天体影响球的借力飞行仿真方法
技术领域：
本发明涉及一种空间探测器借力飞行的仿真方法，特别是借力目标天体的引力影
响球大小与飞行距离之比较大，其影响球内飞行轨道不宜忽略时的天体借力飞行仿真方法。
背景技术：
随着人类对空间探测活动认识的逐渐加深，天体借力飞行技术得到了越来越多的 关注。在空间探测活动中天体借力飞行技术的意义主要有三一是可以节省探测器所带的 燃料；二是可以縮短星际航行时间；三是可使探测器进入非常特殊的轨道进行难得的探测 活动。 在国夕卜的文献C. Ocampo的"Trajectory Analysis for the Lunar FlybyRescue of AsiaSat_3/HGS_l，， 一 文中，Daniel Muhonen， Sylvia Davis等在"Alternative Gravity-Assist Sequences for the ISEE-3 Escape Trajectory，，一文中，Rollin W.禾口 Stanley Ross在"The Venus Swingby Mission Mode and itsRole in the Manned E邓loration of Mars"—文中，以及我国的高长生等在"黄道面外宇宙探测器轨道设计"一 文、马文臻等在"多天体交会借力飞行轨道设计的图解法"一文、卢征仁等在"借月球引力 节省火星探测器发射能量的方法"一文中对空间探测器的天体借力飞行技术进行了的仿真 与分析。虽然这些文献都对借力目标天体近旁沿双曲线轨道运动的借力飞行实质进行了描 述，但是，在分析过程中对整个探测器的轨道进行了简化，忽略了借力目标天体的引力影响 球内的探测器轨道段，即将探测器的整个轨道视为借力飞行前轨道段和借力飞行后轨道段 两部分，而将借力目标天体引力影响球内的轨道段简化为拼接点。这种分析方法在星际探 测活动中影响不大，因为探测器在太阳系内飞行时，主要受太阳引力影响，且借力天体的影 响球半径与其到太阳的距离之比小于0. 01。然而，对于飞行在与地_月系统这种由两颗距 离较近、质量相差不大的系统内的探测器而言，若要借助小质量天体——月球的引力飞行， 考虑到借力目标天体的引力影响球半径(如，月球影响球半径为66200km)与探测器的飞行 最大直线距离(如，地月平均距离384400km)之比较大，探测器在借力目标天体的影响球内 飞行时间较长，对探测器的轨道产生的影响较大，这时忽略借力目标天体影响球的影响会 带来较大的误差。
发明内容
本发明的技术解决问题是针对借力目标天体引力影响球大小与探测器飞行距离 之比较大时，忽略借力目标天体影响球内探测器的运动会带来较大分析误差的不足，本方 法将借力目标天体的影响球简化为球形，并将探测器在影响球内的运动简化为沿双曲线轨 道的二体运动进行分析，该双曲线轨道与飞出影响球后探测器轨道采用圆锥曲线拼接法进 行处理，完成对借力飞行前后的仿真分析工作，使仿真分析更接近探测器在空间的真实运 动情况。
5
本发明的技术解决方案是一种考虑目标天体影响球的借力飞行仿真方法，步骤如下 (1)将借力目标天体影响球假设为球形，确定该影响球的半径rif^及其与借力目标天体与惯性中心之间的平均距离rt，t ee一之比ratio ;当ratio > 0. 01时，考虑探测器在借力目标天体影响球内的运动，转步骤(2)，否则，忽略探测器在借力目标天体影响球内的运动，将探测器在惯性空间中的借力飞行前后的轨迹在拼接点处进行拼接仿真；
(2)根据具体的发射任务，确定探测器在飞向借力目标天体的转移轨道在惯性坐标系下的初始位置6及发射时间t。； (3)假设探测器在借力目标天体影响球内的轨道为双曲线轨道，根据指定的飞越借力目标天体时的轨道高度hp，利用双曲线轨道特性，求解出探测器进入借力目标天体影响球时的惯性中心-借力目标天体中心连线与借力目标天体中心-探测器位置矢量之间的夹角A"进而求出此时探测器的惯性位置矢量大小巧；根据夹角、和位置矢量大小巧，求
得此时的位置巧； (4)根据发射能力或探测任务指定探测器在转移轨道上的耗时、，通过求解兰伯特问题，确定探测器发射及进入借力目标天体影响球时，在惯性坐标系下的速度矢量大小v0、 v丄； (5)利用借力目标天体的星历表、t。及tp将探测器进入借力目标天体影响球时的惯性坐标系下的轨道根数，即步骤(4)中的速度矢量大小v。、^转换为借力目标天体参考系下的轨道根数； (6)根据双曲线轨道运动特性，求解探测器在借力目标天体影响球内的轨道根数和运行的时间t2 ; (7)利用探测器飞出目标天体影响球时的星历和坐标转换关系，得到探测器在飞出借力目标天体影响球时惯性坐标系下的轨道根数，以此为初始根数计算借力飞行后探测器的惯性轨道； (8)将步骤(6)求解得到的借力目标天体影响球内探测器的轨道根数转换为惯
性参考坐标系下的轨道根数，得到惯性参考系下探测器在目标天体影响球内的运动轨道描述，并将其与借力飞行前后的惯性轨道进行圆锥曲线拼接，从而得到完整的借力飞行轨道，
按照得到的完整的借力飞行轨道对探测器进行飞行仿真，模拟探测器在空间的运动情况。
本发明与现有技术相比有益效果为本发明考虑探测器在借力目标天体影响内的轨道运动，即将探测器在该影响球内的运动简化为沿双曲线轨道的二体运动，然后采用圆锥曲线拼接法将双曲线轨道与借力目标天体影响球外的轨道相接。通过这种与实际力学环境更近似的简化设计，利用增加的双曲线轨道段解决了借力飞行过程中目标天体引力影响球内探测器的运动问题，完成对探测器借力飞行前后的仿真分析工作。由于这种仿真方法更接近于探测器的实际运动规律，从而提高了仿真精度。
图1是从惯性中心轨道到借力目标天体影响球的转移轨道示意 图2是本仿真方法中三段轨道拼接的示意图； 图3是算例中三段轨道拼接得到的借力飞行轨道图，图中闭合曲线表示月球公转
6轨道，闭合曲线中心处的圆点表示地球位置。
具体实施例方式
—种考虑目标天体影响球的借力飞行仿真方法，包括下列步骤 (1)将借力目标天体影响球假设为球形，确定该影响球的半径rif^及其与影响
球到惯性中心的距离之比ratio ;当ratio > 0. 01时，考虑探测器在借力目标天体影响球
内的运动，转步骤(2)，否则，忽略探测器在借力目标天体影响球内的运动，飞行轨迹按照中
心天体内两段轨道拼接进行仿真；所述的中心天体在太阳坐标系中为太阳，在地球坐标系
中为地球。 影响球的半径确定公式如下 r =r
influence target—center
附target
附
V "central—body 其中，rtogrt 。mtCT为借力目标天体与惯性中心之间的平均距离，即所述的影响球到 惯性中心的距离； mtarget为借力目标天体的质量；
mcentral body为惯性中心天体的质量。 计算借力目标天体影响球半径与其到惯性中心的距离之比
ratio = ~^^
"target—center 需要说明的是，当ratio《0. 01时，忽略探测器在借力目标天体影响球内的运动， 飞行轨迹按照中心天体内两段轨道拼接进行仿真，简单的讲就是将中心天体内的两段轨道 连接在一起即可，至于所述的两段轨道的确定，不是本发明的创新点，同时也属于本领域技 术人员的公知常识，这里不再进行赘述。 (2)根据具体的发射任务，确定探测器在中心天体-借力目标转移轨道上的初始 位置6及发射时间t。； 至于初始位置^及发射时间t。与所选发射火箭的型号、借力目标天体的选择有关，
不是本发明的创新点，同时也属于本领域技术人员的公知常识，这里不再进行赘述。
(3)假设上述转移轨道为双曲线轨道，根据指定的飞越借力目标天体时的轨道高
度hp，利用双曲线轨道特性，求解出探测器进入借力目标天体影响球时的惯性中心-借力目
标天体中心连线与借力目标天体中心-探测器位置矢量之间的夹角A p进而求出此时探测
器的惯性位置矢量大小巧； 夹角A工的求解过程如下 首先，假设这时探测器沿双曲线轨道的渐近线进入月球影响球，那么，此时探测器
位置矢量与速度矢量间的夹角小为
=2、 t^"(e!sinhF-F)
/An 其中，COShi^二c。sy ， &为双曲线轨道的偏心率；ai为双曲线轨道的半长轴； ^ = 111 COShi^COSh2^-1为双曲线轨道的近点角；当Y G
时，F取正值；当
Y G [180° ，360° )时，F取负值。
函数coshF为F的双曲余弦值。 (7)利用探测器飞出目标天体影响球时的星历和借力目标天体参考系到惯性参考 系的坐标转换关系，得到探测器在飞出借力目标天体影响球时惯性坐标系下的轨道根数， 以此为初始根数计算借力飞行后探测器的惯性轨道； 上述坐标转换关系是坐标转换矩阵乘以借力目标天体参考系的矢量，转换矩阵的 求解方法与(5) —致。 (8)将借力目标天体影响球内探测器的轨道根数转换为惯性参考坐标系下的轨道 根数，得到惯性参考系下探测器在目标天体影响球内的运动轨道。并将其与借力飞行前后 的惯性轨道进行圆锥曲线拼接，从而得到完整的借力飞行轨道。 如图2所示，本发明将探测器的轨道划分为惯性中心体作用下的圆锥曲线轨道 段、借力目标天体影响球内双曲线轨道段和飞出借力目标天体影响球后的惯性体作用下的 轨道段三个阶段，并将三段圆锥曲线轨道拼接以进行仿真。
实施例 以在地月系统内运动的探测器借月球引力飞行为例详细介绍本发明，具体如下
(1)将月球简化为球形，利用下面的解析关系计算月球引力影响球半径和月球引 力影响球半径与其到惯性中心的距离之比。
计算月球引力影响球半径rm。。n—influence :
<formula>formula see original document page 9</formula> 其中，rem = 384400km为地心与月心之间的平均距离；
nim自为月球的质量； [OO"] mearth为地球的质量； 0.1722 ratio > 0. 01，故需要考虑探测器在月球影响球内的双曲线轨道运动。 (2)根据某次发射任务，可以确定探测器在地-月转移轨道上的轨道高度为
200km，所以初始位置r。 = 6578. 137km及发射时间t。 2014年1月13日正午12时； (3)根据某次借力飞行任务确定的近月距rpm = 200+1378 = 1578km。 假设这时探测器沿双曲线轨道的渐近线进入月球影响球，那么，此时探测器位置
矢量与速度矢量间的夹角小为
其中，A =1 + ，为双曲线轨道的偏心率，ai = 7000km为双曲线轨道的半长轴。
若记双曲线渐近线间的夹角为2u/，那么 + 工、tan^ = ±~^ 其中，= a/(e^-l) = 24582084km2。 所以，探测器进入借力目标天体影响球时地月连线与月心探测器位置矢量之间的 夹角入1: 、 = 180° -^ + ￥= 40.148° 根据余弦定理，由A 1确定进入月球影响球处的地心位置矢量大小 n = 丄-2v咖cc^ = 336515.69 km (4)指定进入月球影响球前探测器在转移轨道上需耗时^ = 2. 86天，通过求解兰 伯特问题，求出探测器的初始速度大小v。及进入月球影响球时的速度大小Vl :
v0 = 10. 91560657km/s
Vl = 0. 574426155km/s (5)利用月球星历表、t。及tp将探测器进入月球影响球时的地心赤道惯性坐标系 下的轨道根数转换为月心赤道参考系下的轨道根数； (6)以步骤(5)得到的轨道根数为初始值，通过求解双曲线轨道运动方程，得到 探测器飞出月球影响球时的月心赤道参考系下的轨道根数和在月球影响球内运行的时间 28. 65489998772743小时; (6)利用探测器飞出月球影响球时的月球星历，将探测器的轨道根数转换为地心 赤道参考系下的轨道根数，进一步求得探测器在地球引力作用下的借力后的轨道，如图3 所示。 本发明未详细说明部分属于本领域技术人员公知常识。
10
权利要求
一种考虑目标天体影响球的借力飞行仿真方法，其特征在于步骤如下(1)将借力目标天体影响球假设为球形，确定该影响球的半径rinfluence及其与借力目标天体与惯性中心之间的平均距离rtarget_center之比ratio；当ratio＞0.01时，考虑探测器在借力目标天体影响球内的运动，转步骤(2)，否则，忽略探测器在借力目标天体影响球内的运动，将探测器在惯性空间中的借力飞行前后的轨迹在拼接点处进行拼接仿真；(2)根据具体的发射任务，确定探测器在飞向借力目标天体的转移轨道在惯性坐标系下的初始位置及发射时间t0；(3)假设探测器在借力目标天体影响球内的轨道为双曲线轨道，根据指定的飞越借力目标天体时的轨道高度hp，利用双曲线轨道特性，求解出探测器进入借力目标天体影响球时的惯性中心-借力目标天体中心连线与借力目标天体中心-探测器位置矢量之间的夹角λ1，进而求出此时探测器的惯性位置矢量大小r1；根据夹角λ1和位置矢量大小r1，求得此时的位置(4)根据发射能力或探测任务指定探测器在转移轨道上的耗时t1，通过求解兰伯特问题，确定探测器发射及进入借力目标天体影响球时，在惯性坐标系下的速度矢量大小v0、v1；(5)利用借力目标天体的星历表、t0及t1，将探测器进入借力目标天体影响球时的惯性坐标系下的轨道根数，即步骤(4)中的速度矢量大小v0、v1转换为借力目标天体参考系下的轨道根数；(6)根据双曲线轨道运动特性，求解探测器在借力目标天体影响球内的轨道根数和运行的时间t2；(7)利用探测器飞出目标天体影响球时的星历和坐标转换关系，得到探测器在飞出借力目标天体影响球时惯性坐标系下的轨道根数，以此为初始根数计算借力飞行后探测器的惯性轨道；(8)将步骤(6)求解得到的借力目标天体影响球内探测器的轨道根数转换为惯性参考坐标系下的轨道根数，得到惯性参考系下探测器在目标天体影响球内的运动轨道描述，并将其与借力飞行前后的惯性轨道进行圆锥曲线拼接，从而得到完整的借力飞行轨道，按照得到的完整的借力飞行轨道对探测器进行飞行仿真，模拟探测器在空间的运动情况。FSA00000045632900011.tif,FSA00000045632900012.tif
2. 根据权利要求1所述的一种考虑目标天体影响球的借力飞行仿真方法，其特征在于 所述步骤(1)中影响球的半径确定公式如下<formula>formula see original document page 0</formula>附、'"central—body其中，mtargrt为借力目标天体的质量；mcentral—body为惯性中心天体的质量。
3.根据权利要求1所述的一种考虑目标天体影响球的借力飞行仿真方法，其特征在于 所述步骤(3)中夹角A工的求解过程如下首先，假设这时探测器沿双曲线轨道的渐近线进入月球影响球，那么，此时探测器位置矢量与速度矢量间的夹角小为<formula>formula see original document page 3</formula>其中<formula>formula see original document page 3</formula>为双曲线轨道的偏心率；<formula>formula see original document page 3</formula>，t为借力目标天体的半径；A为双曲线轨道的半长轴；然后，记双曲线渐近线间的夹角为2u/，那么<formula>formula see original document page 3</formula>其中，<formula>formula see original document page 3</formula>则探测器进入借力目标天体影响球时的夹角入i :<formula>formula see original document page 3</formula>
4.根据权利要求1所述的一种考虑目标天体影响球的借力飞行仿真方法，其特征在于所述步骤(3)中的探测器的惯性位置矢量大小巧<formula>formula see original document page 3</formula>
5.根据权利要求1所述的一种考虑目标天体影响球的借力飞行仿真方法，其特征在于所述步骤(4)中探测器在惯性坐标系下的速度矢量大小V。、 Vl公式如下<formula>formula see original document page 3</formula>其中，<formula>formula see original document page 3</formula>其中，P为向借力目标转移的轨道的半通径；a为向借力目标转移的轨道的半长轴；P为惯性中心天体的阴离常数；A v为发射时刻、进入借力目标天体影响球时刻的真近点角之差；AE为发射时刻、进入借力目标天体影响球时刻的偏近点角之差。
6.根据权利要求1所述的一种考虑目标天体影响球的借力飞行仿真方法，其特征在于所述步骤(6)中时间t2的确定过程如下首先，假设探测器进入/飞出借力目标天体影响球时的速度矢量平行于双曲线轨道的渐近线，那么，在目标天体影响球边界处探测器的真近点角Y :其中，小为探测器位置矢量与速度矢量间的夹角；V为双曲线渐近线间夹角的半角； 由此可知，探测器在月球影响球内运行的时间t2 :<formula>formula see original document page 4</formula>苴中，<formula>formula see original document page 4</formula>l二^e。sy ，61为双曲线轨道的偏心率；&1为双曲线轨道的半长轴; 为双曲线轨道的近点角；当Y e
时，F取正值；当[180° ，360° )时，F取负值c
全文摘要
一种考虑目标天体影响球的借力飞行仿真方法，针对借力目标天体引力影响球大小与探测器飞行距离之比较大时，忽略借力目标天体影响球内探测器的运动会带来较大分析误差的不足，本方法将借力目标天体的影响球简化为球形，并将探测器在影响球内的运动简化为沿双曲线轨道的二体运动进行分析，该双曲线轨道与飞出影响球后探测器轨道采用圆锥曲线拼接法进行处理，完成对借力飞行前后的仿真分析工作，使仿真分析更接近探测器在空间的真实运动情况。
文档编号G06F17/50GK101794336SQ201010120400
公开日2010年8月4日 申请日期2010年3月8日 优先权日2010年3月8日
发明者刘胜利, 张燕, 杨芳 申请人:航天东方红卫星有限公司