一种考虑目标天体影响球的借力飞行仿真方法技术

一种考虑目标天体影响球的借力飞行仿真方法，针对借力目标天体引力影响球大小与探测器飞行距离之比较大时，忽略借力目标天体影响球内探测器的运动会带来较大分析误差的不足，本方法将借力目标天体的影响球简化为球形，并将探测器在影响球内的运动简化为沿双曲线轨道的二体运动进行分析，该双曲线轨道与飞出影响球后探测器轨道采用圆锥曲线拼接法进行处理，完成对借力飞行前后的仿真分析工作，使仿真分析更接近探测器在空间的真实运动情况。

【技术实现步骤摘要】

本专利技术涉及一种空间探测器借力飞行的仿真方法，特别是借力目标天体的引力影响球大小与飞行距离之比较大，其影响球内飞行轨道不宜忽略时的天体借力飞行仿真方法。
技术介绍
随着人类对空间探测活动认识的逐渐加深，天体借力飞行技术得到了越来越多的 关注。在空间探测活动中天体借力飞行技术的意义主要有三一是可以节省探测器所带的 燃料；二是可以縮短星际航行时间；三是可使探测器进入非常特殊的轨道进行难得的探测 活动。 在国夕卜的文献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 =rinfluence target—center附target附V "central—body 其中，rtogrt 。mtCT为借力目标天体与惯性中心之间的平均距离，即所述的影响球到 惯性中心的距离； mtarget为借力目标天体的质量； mcentral body为惯性中心天体的质量。 计算借力目标天体影响球半径与其到惯性中心的距离之比 ratio = ~^^"target—center 需要说明的是，当ratio《0. 01时，忽略探测器在借力目标天体影响球内的运动， 飞行轨迹按照中心天体内两段轨道拼接进行仿真，简单的讲就是将中心天体内的两段轨道 连接在一起即可，至于所述的两段轨道的确定，不是本专利技术的创新点，同时也属于本领域技 术人员的公知常识，这里不再进行赘述本文档来自技高网...

【技术保护点】
一种考虑目标天体影响球的借力飞行仿真方法，其特征在于步骤如下：（１）将借力目标天体影响球假设为球形，确定该影响球的半径ｒ↓［ｉｎｆｌｕｅｎｃｅ］及其与借力目标天体与惯性中心之间的平均距离ｒ↓［ｔａｒｇｅｔ＿ｃｅｎｔｅｒ］之比ｒａｔｉｏ；当ｒａｔｉｏ＞０．０１时，考虑探测器在借力目标天体影响球内的运动，转步骤（２），否则，忽略探测器在借力目标天体影响球内的运动，将探测器在惯性空间中的借力飞行前后的轨迹在拼接点处进行拼接仿真；（２）根据具体的发射任务，确定探测器在飞向借力目标天体的转移轨道在惯性坐标系下的初始位置＊↓［０］及发射时间ｔ↓［０］；（３）假设探测器在借力目标天体影响球内的轨道为双曲线轨道，根据指定的飞越借力目标天体时的轨道高度ｈ↓［ｐ］，利用双曲线轨道特性，求解出探测器进入借力目标天体影响球时的惯性中心－借力目标天体中心连线与借力目标天体中心－探测器位置矢量之间的夹角λ↓［１］，进而求出此时探测器的惯性位置矢量大小ｒ↓［１］；根据夹角λ↓［１］和位置矢量大小ｒ↓［１］，求得此时的位置＊↓［１］；（４）根据发射能力或探测任务指定探测器在转移轨道上的耗时ｔ↓［１］，通过求解兰伯特问题，确定探测器发射及进入借力目标天体影响球时，在惯性坐标系下的速度矢量大小ｖ↓［０］、ｖ↓［１］；（５）利用借力目标天体的星历表、ｔ↓［０］及ｔ↓［１］，将探测器进入借力目标天体影响球时的惯性坐标系下的轨道根数，即步骤（４）中的速度矢量大小ｖ↓［０］、ｖ↓［１］转换为借力目标天体参考系下的轨道根数；（６）根据双曲线轨道运动特性，求解探测器在借力目标天体影响球内的轨道根数和运行的时间ｔ↓［２］；（７）利用探测器飞出目标天体影响球时的星历和坐标转换关系，得到探测器在飞出借力目标天体影响球时惯性坐标系下的轨道根数，以此为初始根数计算借力飞行后探测器的惯性轨道；（８）将步骤（６）求解得到的借力目标天体影响球内探测器的轨道根数转换为惯性参考坐标系下的轨道根数，得到惯性参考系下探测器在目标天体影响球内的运动轨道描述，并将其与借力飞行前后的惯性轨道进行圆锥曲线拼接，从而得到完整的借力飞行轨道，按照得到的完整的借力飞行轨道对探测器进行飞行仿真，模拟探测器在空间的运动情况。...

【技术特征摘要】

【专利技术属性】
技术研发人员：张燕，杨芳，刘胜利，
申请(专利权)人：航天东方红卫星有限公司，
类型：发明
国别省市：11[中国|北京]