中文核心期刊

中国科技核心期刊

中国科学引文数据库(CSCD)来源期刊

中国高校优秀科技期刊

中国宇航学会深空探测技术专业委员会会刊

高级检索

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

与地共轨小行星附近测绘轨道的性能分析与评估

石玉 舒磊正 张皓

石玉, 舒磊正, 张皓. 与地共轨小行星附近测绘轨道的性能分析与评估[J]. 深空探测学报(中英文), 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
引用本文: 石玉, 舒磊正, 张皓. 与地共轨小行星附近测绘轨道的性能分析与评估[J]. 深空探测学报(中英文), 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
SHI Yu, SHU Leizheng, ZHANG Hao. Analysis and Evaluation of Mapping Orbits in the Vicinity of Co-Orbital Asteroids[J]. Journal of Deep Space Exploration, 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
Citation: SHI Yu, SHU Leizheng, ZHANG Hao. Analysis and Evaluation of Mapping Orbits in the Vicinity of Co-Orbital Asteroids[J]. Journal of Deep Space Exploration, 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044

与地共轨小行星附近测绘轨道的性能分析与评估

doi: 10.15982/j.issn.2096-9287.2022.20220044
基金项目: 国家自然科学基金青年基金资助项目
详细信息
    作者简介:

    石玉(1993– ),男,博士,助理研究员,主要研究方向:航天器轨道动力学与控制。通讯地址:中国科学院北京新技术基地空间应用工程与技术中心,北京市海淀区邓庄南路9号(100094)电话:18811591132E-mail:shiyu@csu.ac.cn

    张皓(1985– ),男,博士,研究员,主要研究方向 :分布式航天器系统设计与仿真;航天器编队先进控制方法。通讯地址:中国科学院北京新技术基地空间应用工程与技术中心,北京市海淀区邓庄南路9号(100094)E-mail: hao.zhang.zhr@csu.ac.cn

  • ● The stable region of fly-around orbits varies a lot for different asteroids and different orbit radii. ● A method to calculate the coverage rate of the orbits was proposed considering the sunlit condition and topographical shading. ● The performances of mapping orbits were analyzed in the perspective of coverage rate, mapping time, fuel cost and robustness.
  • 中图分类号: V412.4

Analysis and Evaluation of Mapping Orbits in the Vicinity of Co-Orbital Asteroids

  • 摘要: 在小行星探测任务中,航天器抵达小行星附近后需对小行星的各项基本参数测量,为后续任务提供支持。分析了小行星附近环绕型轨道和飞掠型轨道在小行星测绘方面的性能,采用数值方法分析了不同距离下小行星附近的稳定环绕轨道和双曲线飞掠轨道,提出了考虑光照条件和小行星不规则形状遮挡的覆盖率计算方法,从覆盖率、测绘时间、燃料消耗,以及鲁棒性等方面评价测绘轨道。以近球形小行星Bennu 101955和细长型小行星Eros 433为例,比较分析了其附近绕飞轨道和飞掠轨道的测绘性能,为未来小行星探测中的测绘任务规划提供参考。
    Highlights
    ● The stable region of fly-around orbits varies a lot for different asteroids and different orbit radii. ● A method to calculate the coverage rate of the orbits was proposed considering the sunlit condition and topographical shading. ● The performances of mapping orbits were analyzed in the perspective of coverage rate, mapping time, fuel cost and robustness.
  • 图  1  小行星附近航天器所受各个作用力大小

    Fig.  1  Magnitudes of forces of spacecraft in the vicinity of asteroids

    图  2  Bennu附近稳定绕飞的区域

    Fig.  2  Stable region for fly-around orbits near asteroid Bennu

    图  3  Eros附近稳定绕飞的区域

    Fig.  3  Stable region for fly-around orbits near asteroid Eros

    图  4  Bennu和Eros附近的飞掠轨道

    Fig.  4  Flyby orbits around asteroid Bennu and Eros

    图  5  小行星表面观测模型

    Fig.  5  Observation model of the asteroid’s surface

    图  6  小行星表面的遮挡关系

    Fig.  6  Shading on the asteroid’s surface

    图  7  小行星表面遮挡指定三角形的区域

    Fig.  7  Regions of the asteroid’s surface blocking the selected triangle

    图  8  小行星可能被遮挡的三角形区域

    Fig.  8  Triangles that may be blocked by the asteroid’s surface

    图  9  环绕型轨道一个轨道周期的覆盖率

    Fig.  9  Coverage rate of fly-around orbits in 1 revolution

    图  10  Bennu与Eros附近飞掠型轨道的覆盖率

    Fig.  10  Coverage rate of flyby orbits near Bennu and Eros

    图  11  Bennu与Eros附近飞掠型轨道的覆盖区域

    Fig.  11  Covered region of flyby orbits near Bennu and Eros

    图  12  30天环绕轨道时间–覆盖率曲线

    Fig.  12  Time-coverage curve of fly-around orbits for 30 days

    图  13  飞掠轨道的两种方式

    Fig.  13  Two types of flyby orbits

    表  1  小行星参数信息

    Table  1  Physical parameters of asteroids

    参数BennuEros
    引力常数/(km3·s−25.2×10−94.5×10−4
    平均半径/km0.248.42
    自转周期/h5.274.30
    引力场4阶球谐函数[13]4阶球谐函数[14]
    下载: 导出CSV

    表  2  Bennu/Eros附近飞掠轨道的时间、燃料和覆盖率

    Table  2  Time, fuel cost and coverage rate for flyby orbits near Bennu/Eros

    飞掠次数时间/d燃料消耗/(m·s−1覆盖率/%
    南北飞掠54.5/3.20.7/33.374.1/54.2
    98.2/5.81.4/66.990.5/71.2
    1311.8/8.32.1/100.496.6/77.7
    东西飞掠54.5/3.20.7/33.275.1/52.9
    98.2/5.81.4/67.995.0/80.8
    1311.8/8.32.1/100.896.4/87.9
    下载: 导出CSV

    表  3  Bennu/Eros附近环绕轨道的鲁棒性

    Table  3  Robustness of fly-around orbits near Bennu/Eros

    误差大小/%最近点距离/km最远点距离/km覆盖率/%
    −151.05/41.71.83/61.299.0/84.1
    −100.99/41.51.70/53.798.9/83.7
    −50.93/40.41.57/48.998.3/83.8
    00.89/38.51.49/46.297.8/83.7
    50.85/33.71.42/49.897.3/82.7
    100.81/32.51.37/47.496.8/81.8
    150.78/26.61.32/389.594.2/80.8
    下载: 导出CSV

    表  4  Bennu/Eros附近飞掠轨道的鲁棒性

    Table  4  Robustness of flyby orbits around Bennu/Eros

    误差大小/%近星点距离/km覆盖率/%
    −151.27/44.234.0/21.2
    −101.25/43.533.7/19.6
    −51.22/42.832.9/18.5
    01.20/42.132.7/17.5
    51.18/41.431.5/16.8
    101.15/40.730.6/15.8
    151.13/40.129.5/15.1
    下载: 导出CSV
  • [1] TSUDA Y,YOSHIKAWA M,ABE M,et al. System design of the Hayabusa 2—asteroid sample return mission to 1999 JU3[J]. Acta Astronautica,2013,91:356-362. doi:  10.1016/j.actaastro.2013.06.028
    [2] WILLIAMS B,ANTREASIAN P,CARRANZA E,et al. OSIRIS-REx flight dynamics and navigation design[J]. Space Science Reviews,2018,214(4):69. doi:  10.1007/s11214-018-0501-x
    [3] 李翔宇,乔栋,黄江川,等. 小行星探测近轨操作的轨道动力学与控制[J]. 中国科学:物理学 力学 天文学,2019,49(8):69-80.

    LI X Y,QIAO D,HUANG J C,et al. Dynamics and control of proximity operations for asteroid exploration mission[J]. Scientia Sinica Physica,Mechanica & Astronomica,2019,49(8):69-80.
    [4] 林翰文,曾祥远. 小行星433 Eros本体悬停轨道的稳定控制[J]. 动力学与控制学报,2017,15(5):476-480. doi:  10.6052/1672-6553-2017-017

    LIN H W,ZENG X Y. Stability control of body-fixed hovering over asteroid 433 Eros[J]. Journal of Dynamics and Control,2017,15(5):476-480. doi:  10.6052/1672-6553-2017-017
    [5] BROSCHART S B,SCHEERES D J. Control of hovering spacecraft near small bodies:application to asteroid 25143 Itokawa[J]. Journal of Guidance,Control,and Dynamics,2005,28(2):343-354.
    [6] 倪彦硕,宝音贺西,李俊峰. 考虑太阳摄动的小行星附近轨道动力学[J]. 深空探测学报(中英文),2014,1(1):67-74.

    NI Y S,BAOYIN H X,LI J F. Orbit dynamics in the vicinity of asteroids with Solar perturbation[J]. Journal of Deep Space Exploration,2014,1(1):67-74.
    [7] SCHEERES D J,SUTTER B M,ROSENGREN A J. Design,dynamics and stability of the OSIRIS-REx Sun-terminator orbits[J]. Advances in the Astronautical Sciences,2013,148:3263-3282.
    [8] GIANCOTTI M,CAMPAGNOLA S,TSUDA Y,et al. Families of periodic orbits in Hill’s problem with solar radiation pressure:application to Hayabusa 2[J]. Celestial Mechanics and Dynamical Astronomy,2014,120(3):269-286. doi:  10.1007/s10569-014-9564-5
    [9] 胡维多,SCHEERES D J,向开恒. 飞行器近小行星轨道动力学的特点及研究意义[J]. 天文学进展,2009,27(2):152-166. doi:  10.3969/j.issn.1000-8349.2009.02.004

    HU W D,SCHEERES D J,XIANG K H. The characteristics of near asteroid orbital dynamics and its implication to mission analysis[J]. Progress in Astronomy,2009,27(2):152-166. doi:  10.3969/j.issn.1000-8349.2009.02.004
    [10] 姜宇,宝音贺西. 强不规则天体引力场中的动力学研究进展[J]. 深空探测学报(中英文),2014,1(4):250-261.

    JIANG Y,BAOYIN H X. Research trend of dynamics in the gravitational field of irregular celestial body[J]. Journal of Deep Space Exploration,2014,1(4):250-261.
    [11] RAPPAPORT N J,GIAMPIERI G,ANDERSON J D. Perturbations of a spacecraft orbit during a hyperbolic flyby[J]. Icarus,2001,150(1):168-180. doi:  10.1006/icar.2000.6559
    [12] 崔祜涛,史雪岩,崔平远,等. 航天器环绕小行星Ivar的运动分析[J]. 宇航学报,2004(3):251-255. doi:  10.3321/j.issn:1000-1328.2004.03.002

    CUI H T,SHI X Y,CUI P Y,et al. Motion analysis of spacecraft around the asteroid Ivar[J]. Journal of Astronautics,2004(3):251-255. doi:  10.3321/j.issn:1000-1328.2004.03.002
    [13] MCMAHON J W,SCHEERES D J,HESAR S G,et al. The OSIRIS-REx radio science experiment at Bennu[J]. Space Science Reviews,2018,214(1):1-41. doi:  10.1007/s11214-017-0435-8
    [14] SCHEERES D J,MILLER J K,YEOMANS D K. The orbital dynamics environment of 433 Eros:a case study for future asteroid missions[J]. Interplanetary Network Progress Report,2003,42(152):1-26.
    [15] HESAR S G,SCHEERES D J,MCMAHON J W. Sensitivity analysis of the OSIRIS-REx terminator orbits to maneuver errors[J]. Journal of Guidance,Control,and Dynamics,2016,40(1):81-95.
    [16] ANTREASIAN P, HELFRICH C, MILLER J, et al. Preliminary considerations for NEAR's low-altitude passes and landing operations at 433 Eros[C]//AIAA/AAS Astrodynamics Specialist Conference and Exhibit. Boston, MA, USA: AIAA, 1998.
    [17] HU W D,SACHEERES D J. Numerical determination of stability regions for orbital motion in uniformly rotating second degree and order gravity fields[J]. Planetary and Space Science,2004,52(8):685-692.
  • [1] 崔书豪, 王悦, 张瑞康.  双小行星系统不规则引力场中的共振轨道动力学研究 . 深空探测学报(中英文), 2022, 9(4): 382-390. doi: 10.15982/j.issn.2096-9287.2022.20220024
    [2] 郝志鑫, 郑建华, 李明涛.  黄道面内多目标小行星飞越探测任务轨道优化设计 . 深空探测学报(中英文), 2022, 9(4): 373-381. doi: 10.15982/j.issn.2096-9287.2022.20210143
    [3] 丹尼尔J.谢尔斯.  小行星远距离抵近轨道 . 深空探测学报(中英文), 2019, 6(5): 448-455. doi: 10.15982/j.issn.2095-7777.2019.05.005
    [4] 杨洪伟, 宝音贺西.  小行星附近制导与控制研究综述 . 深空探测学报(中英文), 2019, 6(2): 179-188. doi: 10.15982/j.issn.2095-7777.2019.02.010
    [5] 李春来, 刘建军, 严韦, 封剑青, 任鑫, 刘斌.  小行星探测科学目标进展与展望 . 深空探测学报(中英文), 2019, 6(5): 424-436. doi: 10.15982/j.issn.2095-7777.2019.05.003
    [6] 张荣桥, 黄江川, 赫荣伟, 耿言, 孟林智.  小行星探测发展综述 . 深空探测学报(中英文), 2019, 6(5): 417-423,455. doi: 10.15982/j.issn.2095-7777.2019.05.002
    [7] 李宗良, 高俊, 刘国西, 周成, 汤章阳, 邹达人.  小行星探测电推进系统方案研究 . 深空探测学报(中英文), 2018, 5(4): 347-353. doi: 10.15982/j.issn.2095-7777.2018.04.004
    [8] 张韵, 刘岩, 李俊峰.  小行星防御动能撞击效果评估 . 深空探测学报(中英文), 2017, 4(1): 51-57. doi: 10.15982/j.issn.2095-7777.2017.01.008
    [9] 刘磊, 曹建峰, 胡松杰, 唐歌实.  地月L2点周期轨道的月球背面覆盖分析 . 深空探测学报(中英文), 2017, 4(4): 361-366. doi: 10.15982/j.issn.2095-7777.2017.04.008
    [10] 徐青, 王栋, 邢帅, 蓝朝桢.  小行星形貌测绘与表征技术 . 深空探测学报(中英文), 2016, 3(4): 356-362. doi: 10.15982/j.issn.2095-7777.2016.04.007
    [11] 郑茂繁, 耿海, 梁凯, 唐福俊, 黄永杰, 柯于俊.  用于小行星探测的离子推力器技术研究 . 深空探测学报(中英文), 2015, 2(3): 236-240. doi: 10.15982/j.issn.2095-7777.2015.03.008
    [12] 曾祥远, 龚胜平, 李俊峰, 蒋方华, 宝音贺西.  应用太阳帆悬停探测哑铃形小行星 . 深空探测学报(中英文), 2015, 2(1): 48-52. doi: 10.15982/j.issn.2095-7777.2015.01.007
    [13] 张泽旭, 郑博, 周浩, 崔祜涛.  载人小行星探测任务总体方案研究 . 深空探测学报(中英文), 2015, 2(3): 229-235. doi: 10.15982/j.issn.2095-7777.2015.03.007
    [14] 高朝辉, 童科伟, 时剑波, 申麟.  载人火星和小行星探测任务初步分析 . 深空探测学报(中英文), 2015, 2(1): 10-19. doi: 10.15982/j.issn.2095-7777.2015.01.002
    [15] 杨福全, 赵以德, 李娟, 耿海, 张天平, 周海燕.  主带小行星采样返回任务中的离子电推进应用方案 . 深空探测学报(中英文), 2015, 2(2): 168-173. doi: 10.15982/j.issn.2095-7777.2015.02.011
    [16] 周必磊, 陆希, 尤伟.  载人小行星探测的总体方案设想 . 深空探测学报(中英文), 2015, 2(1): 43-47. doi: 10.15982/j.issn.2095-7777.2015.01.006
    [17] 王峰, 杨波, 胡存明, 吴昊, 费晓星.  小行星探测用双谱段相机设计 . 深空探测学报(中英文), 2015, 2(2): 174-179. doi: 10.15982/j.issn.2095-7777.2015.02.012
    [18] 倪彦硕, 宝音贺西, 李俊峰.  考虑太阳摄动的小行星附近轨道动力学 . 深空探测学报(中英文), 2014, 1(1): 67-74.
    [19] 姜宇, 宝音贺西.  强不规则天体引力场中的动力学研究进展 . 深空探测学报(中英文), 2014, 1(4): 250-261. doi: 10.15982/j.issn.2095-7777.2014.04.002
    [20] 尚海滨, 崔平远, 熊旭, 武小宇.  载人小行星探测目标选择与轨道优化设计 . 深空探测学报(中英文), 2014, 1(1): 36-43.
  • 加载中
图(13) / 表 (4)
计量
  • 文章访问数:  37
  • HTML全文浏览量:  13
  • PDF下载量:  3
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-05-10
  • 修回日期:  2022-06-01
  • 网络出版日期:  2022-08-27
  • 刊出日期:  2022-09-30

与地共轨小行星附近测绘轨道的性能分析与评估

doi: 10.15982/j.issn.2096-9287.2022.20220044
    基金项目:  国家自然科学基金青年基金资助项目
    作者简介:

    石玉(1993– ),男,博士,助理研究员,主要研究方向:航天器轨道动力学与控制。通讯地址:中国科学院北京新技术基地空间应用工程与技术中心,北京市海淀区邓庄南路9号(100094)电话:18811591132E-mail:shiyu@csu.ac.cn

    张皓(1985– ),男,博士,研究员,主要研究方向 :分布式航天器系统设计与仿真;航天器编队先进控制方法。通讯地址:中国科学院北京新技术基地空间应用工程与技术中心,北京市海淀区邓庄南路9号(100094)E-mail: hao.zhang.zhr@csu.ac.cn

  • ● The stable region of fly-around orbits varies a lot for different asteroids and different orbit radii. ● A method to calculate the coverage rate of the orbits was proposed considering the sunlit condition and topographical shading. ● The performances of mapping orbits were analyzed in the perspective of coverage rate, mapping time, fuel cost and robustness.
  • 中图分类号: V412.4

摘要: 在小行星探测任务中,航天器抵达小行星附近后需对小行星的各项基本参数测量,为后续任务提供支持。分析了小行星附近环绕型轨道和飞掠型轨道在小行星测绘方面的性能,采用数值方法分析了不同距离下小行星附近的稳定环绕轨道和双曲线飞掠轨道,提出了考虑光照条件和小行星不规则形状遮挡的覆盖率计算方法,从覆盖率、测绘时间、燃料消耗,以及鲁棒性等方面评价测绘轨道。以近球形小行星Bennu 101955和细长型小行星Eros 433为例,比较分析了其附近绕飞轨道和飞掠轨道的测绘性能,为未来小行星探测中的测绘任务规划提供参考。

注释:
1)  ● The stable region of fly-around orbits varies a lot for different asteroids and different orbit radii. ● A method to calculate the coverage rate of the orbits was proposed considering the sunlit condition and topographical shading. ● The performances of mapping orbits were analyzed in the perspective of coverage rate, mapping time, fuel cost and robustness.

English Abstract

石玉, 舒磊正, 张皓. 与地共轨小行星附近测绘轨道的性能分析与评估[J]. 深空探测学报(中英文), 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
引用本文: 石玉, 舒磊正, 张皓. 与地共轨小行星附近测绘轨道的性能分析与评估[J]. 深空探测学报(中英文), 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
SHI Yu, SHU Leizheng, ZHANG Hao. Analysis and Evaluation of Mapping Orbits in the Vicinity of Co-Orbital Asteroids[J]. Journal of Deep Space Exploration, 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
Citation: SHI Yu, SHU Leizheng, ZHANG Hao. Analysis and Evaluation of Mapping Orbits in the Vicinity of Co-Orbital Asteroids[J]. Journal of Deep Space Exploration, 2022, 9(4): 391-399. doi: 10.15982/j.issn.2096-9287.2022.20220044
    • 在小行星探测任务中,对小行星表面测绘是航天器抵近小行星后需完成的最基本且至关重要的任务,既能测量小行星形状、大小、表面地形和引力场等基本信息,又能为后续着陆、采样或撞击任务提供支持。

      在小行星测绘中,对于不同的小行星探测任务,选取的测绘轨道都不尽相同。日本的“隼鸟2号”(Hayabusa 2)探测器对小行星“龙宫”(Ryugu 162173 )的探测选择了悬停与飞掠的测绘方案[1],美国的“欧西里斯号”(OSIRIS-REX)探测器则选择了飞掠与绕飞的方案[2]。因此在小行星探测任务设计中,有必要对不同的轨道在小行星测绘方面的性能进行评估。

      根据航天器与小行星的相对运动关系,测绘轨道主要分为悬停、飞掠和环绕3种方式,且已得到了较为充分的研究,小行星附近航天器的悬停控制研究[3-5],小行星附近的周期轨道的计算与延拓[6-8]、小行星不规则引力场对环绕和飞掠轨道的影响等[9-12]。目前对于小行星测绘轨道的研究主要集中在动力学与控制领域,即根据小行星附近的复杂动力学环境,分析轨道的动力学特性,并在此基础上设计相应的控制器结构。而在轨道的测绘性能方面的研究非常少,即从实际的测绘任务指标出发,分析轨道的特性。

      本文分析了小行星附近的动力学环境,计算了小行星附近的环绕型和飞掠型轨道,并提出了考虑小行星表面遮挡的观测覆盖模型与计算方法。论文以近球形小行星Bennu 101955和细长型小行星Eros 433为例,计算分析了其附近环绕型轨道和飞掠轨道的覆盖性、测绘时间、燃料消耗、鲁棒性等性能指标,从而衡量不同类型测绘轨道对小行星的观测效果,为未来小行星探测任务中的测绘方案规划提供支持。

    • 在小行星附近,小行星的中心引力为航天器所受到的主导作用力,航天器所受到的主要摄动力为小行星非中心引力、太阳光压和太阳引力摄动。

      假设小行星运行在半径1 AU的圆轨道上,其自转轴与公转轴平行,则定义如下坐标系。太阳–小行星旋转坐标系:原点在小行星质心;x轴方向沿太阳指向小行星方向;z轴沿轨道平面法向。定义小行星惯性坐标系与初始时刻的太阳–小行星重合。定义小行星本体坐标系,其中z轴沿自转轴方向,初始时刻本体系与惯性系重合。

      设航天器在小行星惯性系中的位置为${\boldsymbol r}$,航天器运动的动力学方程可表示为

      $$ {\ddot {\boldsymbol{r}}} = {{{\boldsymbol{a}}}_{\text{g}}} + {{{\boldsymbol{a}}}_{\text{s}}} + {{{\boldsymbol{a}}}_{{\text{srp}}}} $$ (1)

      其中:$ {{{\boldsymbol{a}}}_{\text{g}}} $为小行星引力; $ {{{\boldsymbol{a}}}_{\text{s}}} $太阳引力摄动; $ {{{\boldsymbol{a}}}_{{\text{srp}}}} $为太阳光压。

      采用球谐函数模型描述引力场,在本体坐标系中有

      $$ {{\boldsymbol g}_{\text b}} = \nabla \frac{\mu }{r}{\sum\limits_{l = 0}^\infty {\sum\limits_{m = 0}^l {\left( {\frac{{{r_{\text e}}}}{{{r_{\text b}}}}} \right)} } ^l}{P_{lm}}({\rm{sin}}\varphi )({C_{lm}}{\rm{cos}} m\lambda + {{\text{S}}_{lm}}{\rm{sin}}m\lambda ) $$ (2)

      其中:$\mu $为小行星引力常数;${r_{\text e}}$为小行星参考半径;${r_{\text b}}$$ \varphi $$ \lambda $分别为航天器在本体系中的距离、纬度和经度,则惯性系中的小行星引力为

      $$ {{\boldsymbol a}_{\text g}} = {\boldsymbol T} \cdot {{\boldsymbol g}_{\text b}} $$ (3)

      其中:$ {\boldsymbol T} $为本体系到惯性系的转换矩阵。

      航天器所受太阳光压作用为

      $$ {{\boldsymbol a}_{{\text{srp}}}} = \beta {{\hat {\boldsymbol{d}}}_{\text s}} $$ (4)
      $$ \beta = \frac{{{P_0}(1 + \eta )}}{{{d^2}(m/A)}} $$ (5)

      其中:${{\hat {\boldsymbol{d}}}_{\text s}}$为从太阳指向小行星的单位矢量;$ \eta $为航天器反照率;$ P $为太阳常数;$ m/A $为航天器质量与面积之比。

      太阳引力摄动可以表示为

      $$ {{\boldsymbol a}_{\text{s}}} = {\mu _{\text{s}}}\left( {\frac{{{{\boldsymbol r}_{\text{s}}} - {\boldsymbol r}}}{{{{\left| {{{\boldsymbol r}_{\text{s}}} - {\boldsymbol r}} \right|}^3}}} - \frac{{{{\boldsymbol r}_{\text{s}}}}}{{{{\left| {{{\boldsymbol r}_{\text{s}}}} \right|}^3}}}} \right) $$ (6)

      其中:${\mu _{\text{s}}}$为太阳引力常数;${{\boldsymbol r}_{\text{s}}}$为太阳相对小行星的位置矢量。

      本文选取尺寸较小,近似球形的小行星Bennu 101955和尺寸较大,细长型小行星Eros 433 的物理参数作为算例,分析小行星附近的动力学环境,其基本参数如表1所示。

      表 1  小行星参数信息

      Table 1.  Physical parameters of asteroids

      参数BennuEros
      引力常数/(km3·s−25.2×10−94.5×10−4
      平均半径/km0.248.42
      自转周期/h5.274.30
      引力场4阶球谐函数[13]4阶球谐函数[14]

      取航天器的反照率$ \eta $为0.4,质量与面积之比62,航天器在距离小行星不同距离下所受各个作用力的大小如图1所示,横坐标为航天器距离小行星距离与小行星平均半径的比值。可以发现,对于尺寸较小的小行星,其引力场较弱,太阳光压的影响非常明显。当航天器距离小行星Bennu约3倍平均半径时,太阳光压成为最主要的摄动力,且数量级逐渐逼近中心引力,因而分析附近轨道的动力学特性时,太阳光压是最重要的影响因素。当小行星尺寸较大时,引力场较强,此时非中心引力是最主要的摄动力,太阳光压远小于中心引力。此外,航天器在小行星附近时,所受到的太阳引力摄动非常小,其影响最弱。

      图  1  小行星附近航天器所受各个作用力大小

      Figure 1.  Magnitudes of forces of spacecraft in the vicinity of asteroids

    • 环绕型轨道是一种经典的测绘轨道,可以实现对小行星表面的周期性观测。而由于小行星附近复杂的动力学环境,环绕型轨道会受到极强的扰动,从而影响轨道的稳定性,增加撞击和逃逸的风险。因此,在设计环绕型的轨道时,一般需要考虑相对稳定的轨道作为测绘轨道。

      本节采用数值方法分析小行星附近环绕型轨道的稳定性。取飞行时间为60 d,如果航天器到小行星的距离始终维持在给定的范围内,则认为该轨道适合作为小行星的环绕轨道。在数值仿真中,以Eros和Bennu为例,选取轨道的初始偏心率为0,计算不同半长轴和升交点赤经的轨道,选择轨道到小行星的距离的允许偏差范围为30%。在Bennu附近,分别计算初始距离为3、5、10倍小行星平均半径下的轨道演化情况如图2所示,其中标记颜色的区域为飞行60 d后仍处于允许距离范围内的轨道。结果表明,轨道的倾角为90°附近,升交点赤经为90°或270°附近时,轨道能够稳定飞行。其余状态下,轨道均无法实现稳定绕飞。这主要是因为Bennu附近的主要摄动力为太阳光压,绝大多数轨道在太阳光压的影响下偏心率逐渐增大,使得近星点高度逐渐降低,从而无法实现稳定绕飞。事实上,这类能够实现稳定绕飞的轨道称为晨昏线轨道[15],其轨道平面始终近似垂直于太阳矢量方向。

      图  2  Bennu附近稳定绕飞的区域

      Figure 2.  Stable region for fly-around orbits near asteroid Bennu

      Eros附近,初始距离为3、4、5倍和10倍小行星平均半径下的轨道演化情况如图3所示。其中标记颜色的区域为飞行60 d后仍处于允许距离范围内的轨道。可以看到,与Bennu不同,Eros附近的轨道演化情况与轨道半径密切相关,这主要是由于Eros附近的主要摄动力为非中心引力,其作用会随着轨道高度的增加逐渐减弱。当航天器处于3倍平均半径时,顺行轨道均不能稳定绕飞,只有大倾角的逆行轨道才能实现相对稳定的绕飞。这种现象是由C22项摄动引起的,由于C22项摄动与小行星自旋相关,从而对顺行轨道与逆行轨道有不同的影响,与已有文献结论相吻合。即小行星附近的C22项摄动会使得顺行轨道的稳定性发生变化,而对逆行轨道的影响则由于小行星的相对快速自旋而平均化,具体分析见文献[16~17]。而随着航天器距离增加,C22项摄动的影响逐渐减弱,稳定的区域逐渐增加,直到当航天器处于5倍平均半径时,绝大多数轨道均能实现稳定绕飞。这是由于航天器所受的非中心引力作用降低,轨道的稳定范围增加所致。当航天器处于10倍平均半径时,非中心引力作用继续降低,太阳光压的摄动作用逐渐明显,稳定范围减小。在特定的升交点赤经下,轨道无法实现稳定绕飞。因此,考虑到测绘轨道一般需要较大的轨道倾角增加对表面的覆盖区域,在近距离时,环绕型的轨道只适用于对低纬度地区的测绘。对高纬度地区的观测需要增加环绕轨道的高度。

      图  3  Eros附近稳定绕飞的区域

      Figure 3.  Stable region for fly-around orbits near asteroid Eros

    • 小行星引力场较弱,轨道机动可以显著改变飞行轨迹,因而可以通过双曲线轨道飞掠小行星的方式实现对小行星的观测。在设计飞掠型轨道时,不需要考虑轨道的长期稳定性,而更多需要关注飞掠轨道近星点附近的状态。

      采用二体模型中的双曲线轨道设计飞掠轨道,其中近星点的位置由近星点的赤经和赤纬$ (\lambda ,\phi ) $和近星点距离${r_p}$描述,而近星点的速度大小由双曲线的转向角$\delta $确定,速度方向与当地子午面的夹角记为$\gamma $,则轨道偏心率为

      $$ e = \frac{1}{{\sin (\delta /2)}} $$ (7)

      半通径为

      $$ p = {r_p}(1 + e) $$ (8)

      值得注意的是,轨道只有在飞掠近星点附近的较短时段内受到的非中心引力场较为显著,其余轨道段受到的非中心引力场较弱,且会随着小行星的自旋而被平均化,因而近星点时刻对轨道的影响相对较弱。

      为了修正摄动力对轨道的影响,设航天器初状态的时间和位置为${t_0}$${{\boldsymbol x}_0}$,末状态的时间和位置为${t_{\text {f}}}$${{\boldsymbol x}_{\text f}}$,构造两点间转移的Lambert问题,即求解初始速度${{v}_0}$满足

      $$ F({{\boldsymbol v}_0}) - {{\boldsymbol x}_{\text f}} = 0 $$ (9)

      其中:$F({\boldsymbol v})$${t_0}$时刻,初始状态为$ [{{\boldsymbol x}_0},{{\boldsymbol v}_0}] $运动至${t_{\text f}}$时刻的位置。

      以二体模型下的双曲线轨道为初始猜测,数值迭代求解方程零点,从而得到修正后的飞掠轨道。

      近星点距离为5倍平均半径时,Bennu 和Eros附近从不同经度飞掠赤道面的双曲线轨道如图4所示,其中黑色轨迹为二体下的轨道,蓝色轨迹为受到扰动后的轨道,绿色轨迹为修正后的轨道。可以看到,在Bennu附近轨道飞掠轨道受到的扰动较强,使得近星点的位置产生变化,这与Bennu附近的主要摄动力是太阳光压相关,经过数值修正后,轨道可以实现对指定区域上空的飞掠。Eros附近光压影响相对较弱,因而飞掠轨道受到的影响较小。

      图  4  Bennu和Eros附近的飞掠轨道

      Figure 4.  Flyby orbits around asteroid Bennu and Eros

      总的来看,飞掠轨道大体沿双曲线飞行,虽然扰动会使得轨道产生偏差,但从覆盖性的角度看,由于小行星快速自旋,轨道受到的扰动对覆盖性的影响较小。

    • 小行星表面区域的可见性取决于航天器、小行星和太阳的相对几何关系。如果太阳的入射光线照射到小行星表面,经过反射后的光线能够被航天器捕捉到,则该区域视为可见,如图5所示。对小行星表面上的区域,设其法向量为${\boldsymbol n}$、光线的入射角为$\phi $、反射角为$\varepsilon $、相对于航天器的视角为$\theta $、航天器的位置为${{\boldsymbol r}_{\text sc}}$,则可观测条件需满足

      图  5  小行星表面观测模型

      Figure 5.  Observation model of the asteroid’s surface

      $$\left\{ \begin{gathered} 0 \leqslant \phi \leqslant {\phi _{\text{m}}} \\ 0 \leqslant \varepsilon \leqslant {\varepsilon _{\text{m}}} \\ 0 \leqslant \theta \leqslant {\theta _{\text{m}}} \\ \end{gathered} \right.$$ (10)

      其中:${\phi _{\text{m}}}$${\varepsilon _{\text{m}}}$为限定的最大入射角与反射角,即当地太阳高度角和卫星高度角;${\theta _{\text{m}}}$为航天器的视场角。除此之外,航天器需要在特定的高度范围内对地面进行观测,即高度约束满足

      $$ {h_{\min }} \leqslant \left| {{{\boldsymbol r}_{{\text{sc}}}}} \right| \leqslant {h_{\max }} $$ (11)

      影响上述几何参数范围的最主要因素是航天器上的载荷性能,如相机的分辨率、视场角等极大的影响轨道对小行星的覆盖性。然而,这类参数的影响趋势是显而易见的。例如,增大视场角,增大允许的最大入射角和反射角,均可以极大的提高观测效果。本文主要研究测绘轨道性能分析,因而在后续的仿真计算中,不再分析这类几何参数对观测性能的影响,直接选定允许的最大入射角与反射角均为75°,允许的观测高度不超过标称高度的130%,航天器的视场角设为2.5°,且航天器始终指向小行星中心。

    • 小行星的形状不规则,表面可能存在互相遮挡的情况,还需要增加入射光线和反射光线遮挡的判据。本文采用多面体描述小行星的几何形状,其表面的每个面是三角形,如果光线同时穿过多面体上两个三角形,则距离较远的三角形被遮挡。因此,遮挡问题转化为求解空间中线段与三角形的相交问题。

      设反射光线的起点为小行星表面一点$ P $,终点为航天器$ Q $$ O $为天体质心,三角形$ ABC $(顶点逆时针排列)的法向量为${{\boldsymbol n}_t}$,方向指向小行星外侧,如图6所示,为了判断三角形$ ABC $是否会遮挡光线,首先计算三角形面法向量与光线的点积,若

      图  6  小行星表面的遮挡关系

      Figure 6.  Shading on the asteroid’s surface

      $$ {\boldsymbol {PQ}} \cdot {{\boldsymbol n}_t} \geqslant 0 $$ (12)

      则光线与三角形平行或者从三角形背面进入,不构成遮挡关系。

      其次,计算光线与三角形所在平面的交点。设三角形三个顶点的坐标为$ ({x}_{A},{y}_{A},{z}_{A}) $$ ({x}_{B},{y}_{B},{z}_{B}) $$ ({x}_{C},{y}_{C},{z}_{C}) $,则三角形所在平面方程为

      $$ ax + by + cz = 1 $$ (13)

      其中

      $$ \left[ {\begin{array}{*{20}{c}} a \\ b \\ c \end{array}} \right] = {\left[ {\begin{array}{*{20}{c}} {{x_A}}&{{y_A}}&{{z_A}} \\ {{x_B}}&{{y_B}}&{{z_B}} \\ {{x_C}}&{{y_{\text{C}}}}&{{z_C}} \end{array}} \right]^{ - 1}}\left[ {\begin{array}{*{20}{c}} 1 \\ 1 \\ 1 \end{array}} \right] $$ (14)

      光线所在直线方程为

      $$ \frac{{x - {x_P}}}{{{x_Q} - {x_P}}} = \frac{{y - {y_P}}}{{{y_Q} - {y_P}}} = \frac{{z - {z_P}}}{{{z_Q} - {z_P}}} = t $$ (15)

      联立方程求解交点,可得

      $$ t=\frac{1-(a{x}_{p}+b{y}_{P}+c{z}_{p})}{a({x}_{Q}-{x}_{P})+b({y}_{Q}-{y}_{P})+c({z}_{Q}-{z}_{P})} $$ (16)

      $ t\notin (0,1) $,则交点在不在线段$ {\boldsymbol {PQ}} $内,三角形不构成遮挡关系。

      如果交点在线段上,则需要判断交点是否在三角形上,若交点$ M $满足

      $$ \begin{gathered} {{\boldsymbol n}_t} \cdot ({\boldsymbol {AB}} \times {\boldsymbol {AM}} ) \geqslant 0 \\ {{\boldsymbol n}_t} \cdot ({\boldsymbol {BC}} \times {\boldsymbol {BM}} ) \geqslant 0 \\ {{\boldsymbol n}_t} \cdot ({\boldsymbol {CA}} \times {\boldsymbol {CM}} ) \geqslant 0 \\ \end{gathered} $$ (17)

      则交点在三角形$ ABC $上。

      本文在计算可见性时,以三角形形心表示该区域,如果多面体表面三角形的形心可见,则视为该三角形可见。表面覆盖率定义为已观测三角形面积的百分比。多面体表面的三角形数量众多,在计算被观测区域时,需要对其余每个三角面进行遮挡判断,会极大地增加计算量。因此可以根据小行星的多面体模型,计算得到多面体每个三角形所在地平线以上的其余三角形,并预先存储。在计算时,只需计算对应地平线上的三角形的遮挡关系即可。小行星表面遮挡指定三角形的区域如图7所示,其中绿色区域为蓝色三角形当地地平线以上,且仰角高于某一最小仰角(本文取15°)的部分,对蓝色三角形的观测构成潜在的遮挡关系。Eros表面所有可能被遮挡的三角形如图8所示,可以看到,Eros表面可能被遮挡的区域主要集中在小行星中部凹陷处,对相应的观测带来更多约束。

      图  7  小行星表面遮挡指定三角形的区域

      Figure 7.  Regions of the asteroid’s surface blocking the selected triangle

      图  8  小行星可能被遮挡的三角形区域

      Figure 8.  Triangles that may be blocked by the asteroid’s surface

    • 对于环绕轨道而言,为提高覆盖率,通常会选择高轨道倾角的环绕型轨道,以实现对小行星高纬度地区的覆盖,因此,根据第2节计算得到的环绕型轨道,选定环绕型轨道的倾角为90°,轨道半径为5倍平均半径,计算不同升交点赤经的轨道在一个轨道周期内对小行星表面的覆盖率,其中观测模型的参数在第3节给出,计算结果如图9所示。结果表明,对Bennu而言,由于太阳在180°方向附近,因此当升交点赤经在180°附近时,覆盖率最高,在90°和270°时,覆盖率较低。注意到覆盖率随升交点赤经的变化曲线并非严格对称,这主要是因为轨道受到的扰动非常大,从而改变的轨道的高度,影响了覆盖率。Eros附近环绕轨道的覆盖率受升交点赤经的影响不大,可见Eros的不规则形状极大地影响了轨道对表面的覆盖率。

      图  9  环绕型轨道一个轨道周期的覆盖率

      Figure 9.  Coverage rate of fly-around orbits in 1 revolution

      对飞掠型轨道而言,由于经度的影响在环绕型轨道中已经体现,这里讨论近星点纬度和飞掠方向对覆盖性的影响。取近星点距离为5倍平均半径,转向角为60°,纬度分别取南北纬80°,南北纬40°和0°,近星点速度与当地子午面的夹角为0°~360°,计算轨道的覆盖率,结果如图10所示。可以发现,Eros的不规则形状导致其覆盖率低于Bennu。随着纬度的升高,覆盖率显著下降。这是因为赤道面地区的面积大,光照条件好,所以覆盖率高。对于单次飞掠而言,覆盖率最大的方式为从赤道附近,以较低的轨道倾角,逆着小行星的自旋方向的飞掠轨道,其覆盖的面积如图11所示。

      图  10  Bennu与Eros附近飞掠型轨道的覆盖率

      Figure 10.  Coverage rate of flyby orbits near Bennu and Eros

      图  11  Bennu与Eros附近飞掠型轨道的覆盖区域

      Figure 11.  Covered region of flyby orbits near Bennu and Eros

    • 无论是环绕轨道还是飞掠轨道,无法通过短时间绕飞或单次飞掠实现对小行星表面绝大多数地区的覆盖,因而需要讨论轨道实现大范围覆盖所需的时间。对于环绕型轨道而言,长时间的环绕需要选择稳定的环绕轨道,根据第2节对轨道稳定性的分析结果,选择Bennu附近的初始倾角90°、升交点赤经80°的轨道,以及Eros附近初始倾角90°、升交点赤经180°的轨道,计算轨道飞行20 d的覆盖性曲线,如图12所示。其中Bennu附近的环绕轨道在20 d内能达到90%以上的覆盖率,并趋于稳定,Eros附近的环绕型轨道在10 d内达到80%以上的覆盖率并趋于稳定。Eros附近的环绕轨道的光照条件明暗交替,其覆盖曲线呈现明显的阶跃特征,而Bennu附近的环绕轨道位于晨昏线附近,阶跃特征不明显,但随着摄动作用的影响,轨道面产生小扰动,阶跃特征逐渐呈现。此外,由于环绕型轨道可以实现较长时间的稳定绕飞,在理想情况下环绕型轨道在测绘过程中燃料消耗为0。

      图  12  30天环绕轨道时间–覆盖率曲线

      Figure 12.  Time-coverage curve of fly-around orbits for 30 days

      对于飞掠轨道而言,一般需要在飞掠过程中施加机动脉冲将多个单次飞掠轨道拼接成多次飞掠轨道,从而提高覆盖率。本文考虑南北向和东西向的两种飞掠方式,近星点距离为5倍平均半径,最远点距离为15倍平均半径,轨道转向角60°,其中南北向轨道的沿赤经为180°的双曲线轨道多次往返飞掠,东西向轨道的近星点纬度从南纬90°到北纬90°均匀分布,结果如图13所示。

      图  13  飞掠轨道的两种方式

      Figure 13.  Two types of flyby orbits

      Bennu和Eros附近飞掠轨道的时间、覆盖率和燃料消耗$ (\Delta v) $随飞掠次数变化的统计结果如表2所示。可以看到,随着飞掠次数的增加,两种飞掠方式的时间,燃料消耗和覆盖率均增加。与环绕轨道相比,两种飞掠方式的轨道均能通过提高飞掠次数达到与环绕轨道相当的覆盖率,并且所需时间减少了一半以上。然而在燃料消耗方面,Bennu附近飞掠轨道的消耗远小于Eros附近的轨道,这主要是由于Bennu的引力场较弱,轨道机动作用更明显。因此,对于本文所考虑的测绘问题,多次飞掠模式不适合作为引力场较强的小行星探测任务中的测绘轨道。

      表 2  Bennu/Eros附近飞掠轨道的时间、燃料和覆盖率

      Table 2.  Time, fuel cost and coverage rate for flyby orbits near Bennu/Eros

      飞掠次数时间/d燃料消耗/(m·s−1覆盖率/%
      南北飞掠54.5/3.20.7/33.374.1/54.2
      98.2/5.81.4/66.990.5/71.2
      1311.8/8.32.1/100.496.6/77.7
      东西飞掠54.5/3.20.7/33.275.1/52.9
      98.2/5.81.4/67.995.0/80.8
      1311.8/8.32.1/100.896.4/87.9
    • 由于各种扰动的存在,航天器总会偏离预定的轨道。本节考虑模型误差对测绘轨道的影响(主要由小行星的引力场数的不确定性引起),即讨论当小行星所受到的作用力存在误差的情况下,设计的环绕与飞掠型轨道是否能够实现目标。这里将航天器所受的实际作用力按比例增大和缩小,分析轨道对于模型误差的鲁棒性。

      取4.2节的环绕轨道为例,仿真时间为30 d,其最近点距离、最远点距离和覆盖率如表3所示。结果表明,在10%的模型误差范围内,Bennu和Eros的模型误差对环绕轨道的影响在允许范围内,当误差达到15%时,Eros附近的环绕轨道无法保持稳定。在覆盖性方面,模型误差对覆盖性的影响相对较小,且随着航天器所受作用力的增大,覆盖率逐渐降低。

      表 3  Bennu/Eros附近环绕轨道的鲁棒性

      Table 3.  Robustness of fly-around orbits near Bennu/Eros

      误差大小/%最近点距离/km最远点距离/km覆盖率/%
      −151.05/41.71.83/61.299.0/84.1
      −100.99/41.51.70/53.798.9/83.7
      −50.93/40.41.57/48.998.3/83.8
      00.89/38.51.49/46.297.8/83.7
      50.85/33.71.42/49.897.3/82.7
      100.81/32.51.37/47.496.8/81.8
      150.78/26.61.32/389.594.2/80.8

      对于飞掠轨道而言,选取4.2节中的南北向飞掠轨道,仿真结果如表4所示。结果表明模型误差对飞掠轨道的影响相对较小,可以看到,在15%的误差范围内,Bennu和Eros附近的飞掠轨道依然能够完成预定的观测,其近星点的距离和覆盖率的变化相对较小。随着航天器所受作用力增大,近星点高度降低,相应的覆盖面积减小。因此飞掠型的轨道相比于环绕型轨道更适合模型不确定性较大的小行星。

      表 4  Bennu/Eros附近飞掠轨道的鲁棒性

      Table 4.  Robustness of flyby orbits around Bennu/Eros

      误差大小/%近星点距离/km覆盖率/%
      −151.27/44.234.0/21.2
      −101.25/43.533.7/19.6
      −51.22/42.832.9/18.5
      01.20/42.132.7/17.5
      51.18/41.431.5/16.8
      101.15/40.730.6/15.8
      151.13/40.129.5/15.1
    • 本文针对小行星任务中的探测侦察问题,分析比较了小行星附近的环绕型与飞掠型测绘轨道的特性。总的来看,环绕型轨道和飞掠型轨道均能实现对小行星的测绘。具体而言,环绕型轨道一般观测周期较长,燃料消耗低,但强烈依赖小行星附近的动力学环境,对于引力场较弱的小行星,多数环绕型轨道无法实现稳定绕飞。飞掠型轨道能够快速实现对小行星的观测,且对小行星的模型不确定具有较强的鲁棒性,但是需要消耗燃料进行机动实现对目标的多次飞掠,因而对于引力场较强的小行星,其燃料消耗较大。本文以近球形小行星Bennu 101955 和细长型小行星Eros 433为例,考虑小行星附近的复杂动力学环境与小行星不规则形状的相互遮挡关系,计算了环绕型轨道和飞掠型轨道的稳定性、覆盖率、测绘时间、燃料消耗以及鲁棒性等指标,并分析了环绕型轨道和飞掠型轨道的测绘特性与适用范围,为未来小行星探测任务中的测绘方案规划提供支持。

参考文献 (17)

目录

    /

    返回文章
    返回