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

高级检索

留言板

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

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

火星探测器制动捕获策略研究

王靓玥 郭延宁 马广富

王靓玥, 郭延宁, 马广富. 火星探测器制动捕获策略研究[J]. 深空探测学报(中英文), 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
引用本文: 王靓玥, 郭延宁, 马广富. 火星探测器制动捕获策略研究[J]. 深空探测学报(中英文), 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
WANG Liangyue, GUO Yanning, MA Guangfu. Research About Strategies of Mars Probe Capture[J]. Journal of Deep Space Exploration, 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
Citation: WANG Liangyue, GUO Yanning, MA Guangfu. Research About Strategies of Mars Probe Capture[J]. Journal of Deep Space Exploration, 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001

火星探测器制动捕获策略研究

doi: 10.15982/j.issn.2095-7777.2020.20171123001
基金项目: 国家自然科学基金资助项目(61973100,61673135,61876050)
详细信息
    作者简介:

    王靓玥(1996–),女,博士,主要研究方向:航天器姿轨一体化制导与控制。通讯地址:哈尔滨市西大直街92号哈尔滨工业大学327信箱主楼421室(150001)E-mail:wlyue96@163.com

    马广富(1963–),男 教授,主要研究方向:航天器制导与控制、最优控制。通讯地址:哈尔滨工业大学327信箱(150001)电话:(0451)86402726 E-mail:magf@hit.edu.cn

    通讯作者:

    郭延宁(1985–),男,博士,副教授,主要研究方向:航天器控制,深空探测导航、制导与控制。本文通讯作者。通讯地址:哈尔滨市西大直街92号哈尔滨工业大学327信箱主楼421室(150001)电话:(0451)86413411-8606 E-mail:guoyn@hit.edu.cn

  • ● Three typical capture models using thrusters are systematically compared and the merits and demerits aresummarized. ● The optimal ignition point is discovered by the numerical simulation. ● A closed-loop guidance strategy based on the Zero Effort Miss/Zero Effort Velocity is proposed to implement specific transition track.
  • 中图分类号: V47

Research About Strategies of Mars Probe Capture

  • 摘要: 针对火星探测器制动捕获问题进行研究,通过对火星探测器捕获问题的动力学建模、典型工况数值仿真,对3种典型的直接捕获策略:定向捕获、匀速捕获、变速捕获进行比较、分析,系统地总结了3种捕获策略的优缺点,得到变速捕获推力器燃料消耗低,但需要实时速度信息等结论。进一步地,通过数值分析点火时刻对制动捕获效果的影响,发现探测器最优点火时刻为到达近火点前制动耗时一半左右所对应的时刻。为实现特定捕获过渡轨道,设计了基于ZEM/ZEV闭环制导的捕获策略。通过与变速捕获对比表明:所设计捕获策略的制动效果与变速捕获相当,且能满足特定的末端位置/速度约束,可为实际任务提供一种备选方案。
    Highlights
    ● Three typical capture models using thrusters are systematically compared and the merits and demerits aresummarized. ● The optimal ignition point is discovered by the numerical simulation. ● A closed-loop guidance strategy based on the Zero Effort Miss/Zero Effort Velocity is proposed to implement specific transition track.
  • 图  1  探测器运行轨迹图

    Fig.  1  Trajectory of probe

    图  2  探测器运行的轨道平面图

    Fig.  2  Orbital plane of probe running

    图  3  不同点火时刻耗时、耗燃料

    Fig.  3  Time and fuel consumption of ignition time

    图  4  不同点火时刻距火星最远、近距离

    Fig.  4  The furthest and nearest distance from Mars of ignition time

    表  1  仿真条件

    Table  1  Simulation conditions

    项目参数
    探测器质量/kg4.948 × 103
    发动机推力/N 3 × 103
    比冲/s 314
    点火时刻速度矢量/(km·s–1[–1.980;4.416;–2.544]
    点火时刻位置矢量/km[3 330.9;1 981.4;829.8]
    下载: 导出CSV

    表  2  仿真结果及优缺点

    Table  2  Simulation result, merits and demerits

    项目定向捕获匀速捕获变速捕获
    捕获时间/s1 446.81 360.01 358.8
    消耗燃料/kg1 410.01 325.11 323.9
    距火星表面
    最近距离/km
    635.425291.875247.065
    参数个数565
    优点实现简单实现简单,推力器利用率高推力器利用率高
    缺点推力器利用率低制动效果严重依赖于角速度的选取推力器实现复杂
    下载: 导出CSV

    表  3  不同工况仿真条件及结果

    Table  3  Simulation results of different conditions

    发动机推力/
    N
    比冲/
    s
    最优点火刻/
    s
    制动耗时/
    s
    燃料消耗/
    kg
    3 000314–6881 268.41 235.9
    3 500314–5721 082.41 230.4
    3 000250–6691 246.81 363.1
     注:探测器到达近火点时刻为0 s。
    下载: 导出CSV

    表  4  ZEV、ZEM/ZEV仿真条件

    Table  4  Simulation conditions of ZEV and ZEM/ZEV

    时间速度矢量/(km·s–1位移矢量/(km)
    制动捕获开始 [–1.980;4.416;–2.544] [3330.9;1981.4;829.8]
    制动捕获结束[–2.558;1.726;–1.828] [–189.9;6010;–2292]
    下载: 导出CSV

    表  5  新型捕获策略仿真结果

    Table  5  Simulation results of the new capture strategy

    捕获
    策略
    捕获
    时间/s
    消耗
    燃料/kg
    终点位置/km终点速度/(km·s–1
    变速
    捕获
    1 358.81 324[–189.9;6 010;
    –2 292]
    [–2.558;1.726;
    –1.828]
    ZEV1 3591 324[–140.2;6 012;
    –2 270]
    [–2.558;1.728;
    –1.828]
    ZEM/
    ZEV
    1 3591 323[–189.6;6 010;
    –2 292]
    [–2.559;1.726;
    –1.828]
    下载: 导出CSV
  • [1] 于登云,孙泽洲,孟林智,等. 火星探测发展历程与未来展望[J]. 深空探测学报,2016,3(2):108-113.

    YU D Y,SUN Z Z,MENG L,et al. The development process and prospects for Marsexploration[J]. Journal of Deep Space Exploration,2016,3(2):108-113.
    [2] 方宝东, 张建刚, 张文巧. 火星探测发展趋势及关键技术浅析[C]//中国宇航学会深空探测技术专业委员会第三届学术会议. [S.l.]: 中国宇航学会深空探测技术专业委员会, 2006.
    [3] 李军锋,龚胜平. 有限推力模型火星探测捕获策略分析[J]. 中国科学:物理学,力学,天文,2013(6):781-786.

    LI J F,GONG S P. Analysis of capture strategies for Mars explorer with finite-thrust[J]. Scienti Sinica(Physic,Mechanic&Astronomic),2013(6):781-786.
    [4] 闵学龙,潘腾,郭海林. 火星探测器使命轨道捕获策略研究[J]. 航天器工程,2008,17(6):39-43. doi:  10.3969/j.issn.1673-8748.2008.06.004

    MIN X L,PAN T,GUO H L. Analysis of orbit capture method for Mars vehicle[J]. Spacecraft Engineering,2008,17(6):39-43. doi:  10.3969/j.issn.1673-8748.2008.06.004
    [5] 吕静,张明明,龚胜平. 旋转大气下火星探测器轨道捕获[J]. 北京航空航天大学学报,2013,39(3):348-354.

    LV J,ZHANG M M,GONG S P. Aerocaptureperiod under rotating atmospheric environment for Mars vehicle[J]. Journal of Beijing University of Aeronautics and Astronautics,2013,39(3):348-354.
    [6] 罗续盛,麻娜,荆武兴,等. 采用有限推力的火星制动捕获策略[J]. 西北工业大学学报,2017,35(2):348-354. doi:  10.3969/j.issn.1000-2758.2017.02.028

    LUO X S,MA N,XING J W,et al. Capture and brake strategy for Mars with explorer with finite-thrust[J]. Journal of Northwestern Polytechnical University,2017,35(2):348-354. doi:  10.3969/j.issn.1000-2758.2017.02.028
    [7] 郭延宁, 于彦博, 马广富. 基于路径点优化的火星着陆器动力下降段反馈制导律研究[J]. 中国宇航学会深空探测技术专业委员会第九届学术年会. 杭州: 中国宇航学会深空探测技术专业委员会, 2010.
    [8] GUO Y N,HAWKINS M,WIE B. Applications of generalized Zero-Effort-Miss/Zero-Effort-Velocity feedback guidance algorithm[J]. Journal of Guidance,Control,and Dynamics,2013,36(3):810-820.
    [9] 肖冰. 带有输入饱和的挠性航天器姿态跟踪鲁棒控制研究[D]. 哈尔滨: 哈尔滨工业大学, 2010.

    XIAO B. Research on robust attitude tracking control for flexible spacecraft with input saturation[D]. Harbin: Harbin Institute of Technology, 2010
    [10] 盛伟强. 火星近火点捕获制动姿轨一体化控制[D]. 哈尔滨: 哈尔滨工业大学, 2015.

    SHENG W Q. Study on coupling control of attitude and orbit for Mars orbit insertion at periareon[D]. Harbin: Harbin Institute of Technology, 2015.
    [11] HECHLER M,YÁ EZ A. Mars express orbit design[J]. Acta Astronautica,2003,53(4):497-507.
    [12] 方宝东,吴美平,张伟. 火星引力捕获动力学与动态误差分析[J]. 力学学报,2015,47(1):15-23. doi:  10.6052/0459-1879-14-327

    FANG B D,WU M P,ZHANG W. Mars gravity capture dynamic model and error analysis[J]. Chinese Journal of Theoretical and Applied Mechanics,2015,47(1):15-23. doi:  10.6052/0459-1879-14-327
  • [1] 张宝明, 朱岩, 王连国, 杨建峰, 周斌, 徐卫明, 孙树全, 蔡治国, 徐欣锋, 杜庆国.  中国首次火星探测任务火星车有效载荷定标试验 . 深空探测学报(中英文), 2020, 7(5): 481-488. doi: 10.15982/j.issn.2096-9287.2020.20200043
    [2] 滕锐, 韩宏伟, 乔栋.  火星探测最优离轨制导方法研究 . 深空探测学报(中英文), 2020, 7(2): 184-190. doi: 10.15982/j.issn.2095-7777.2020.20190315001
    [3] 薛彬, 刘生润, 杨建峰.  用于火星表面生命信息探测的激光拉曼技术进展 . 深空探测学报(中英文), 2019, 6(5): 503-512. doi: 10.15982/j.issn.2095-7777.2019.05.012
    [4] 徐侃彦, 马玲玲, 印红, 张轶男.  火星无人探测与行星保护 . 深空探测学报(中英文), 2019, 6(1): 9-15. doi: 10.15982/j.issn.2095-7777.2019.01.002
    [5] 孟庆宇, 付中梁, 董吉洪, 王栋.  火星探测高分辨率可见光相机光学系统设计 . 深空探测学报(中英文), 2018, 5(5): 458-464. doi: 10.15982/j.issn.2095-7777.2018.05.008
    [6] 李春来, 刘建军, 耿言, 曹晋滨, 张铁龙, 方广有, 杨建峰, 舒嵘, 邹永廖, 林杨挺, 欧阳自远.  中国首次火星探测任务科学目标与有效载荷配置 . 深空探测学报(中英文), 2018, 5(5): 406-413. doi: 10.15982/j.issn.2095-7777.2018.05.002
    [7] 杨甲森, 刘明洁, 陈托, 智佳, 张华伟, 王炜, 陈志敏.  中国首次火星探测任务有效载荷地面综合测试系统设计 . 深空探测学报(中英文), 2018, 5(5): 442-449. doi: 10.15982/j.issn.2095-7777.2018.05.006
    [8] 申智春, 林小艳, 程坤, 王海鹏.  火星探测器器箭分离冲击响应影响分析与评价 . 深空探测学报(中英文), 2018, 5(5): 483-487. doi: 10.15982/j.issn.2095-7777.2018.05.012
    [9] 刘建军, 苏彦, 左维, 任鑫, 孔德庆, 温卫斌, 张洪波, 李春来.  中国首次火星探测任务地面应用系统 . 深空探测学报(中英文), 2018, 5(5): 414-425. doi: 10.15982/j.issn.2095-7777.2018.05.003
    [10] 耿言, 周继时, 李莎, 付中梁, 孟林智, 刘建军, 王海鹏.  我国首次火星探测任务 . 深空探测学报(中英文), 2018, 5(5): 399-405. doi: 10.15982/j.issn.2095-7777.2018.05.001
    [11] 刘庆会.  火星探测VLBI测定轨技术 . 深空探测学报(中英文), 2018, 5(5): 435-441. doi: 10.15982/j.issn.2095-7777.2018.05.005
    [12] 叶斌龙, 赵健楠, 黄俊.  美国2020火星车着陆区遴选进展及对2020中国火星任务着陆探测部分的一些思考 . 深空探测学报(中英文), 2017, 4(4): 310-324. doi: 10.15982/j.issn.2095-7777.2017.04.002
    [13] 朱岩, 白云飞, 王连国, 沈卫华, 张宝明, 王蔚, 周盛雨, 杜庆国, 陈春红.  中国首次火星探测工程有效载荷总体设计 . 深空探测学报(中英文), 2017, 4(6): 510-514,534. doi: 10.15982/j.issn.2095-7777.2017.06.002
    [14] 李贺, 全齐全, 王鑫剑, 姜生元, 邓宗全.  一种基于压电驱动的火星岩石钻探器的研制 . 深空探测学报(中英文), 2016, 3(2): 156-161. doi: 10.15982/j.issn.2095-7777.2016.02.010
    [15] 陈晓, 尤伟, 黄庆龙.  火星探测巡航段天文自主导航方法研究 . 深空探测学报(中英文), 2016, 3(3): 214-218. doi: 10.15982/j.issn.2095-7777.2016.03.003
    [16] 于登云, 孙泽洲, 孟林智, 石东.  火星探测发展历程与未来展望 . 深空探测学报(中英文), 2016, 3(2): 108-113. doi: 10.15982/j.issn.2095-7777.2016.02.002
    [17] 高朝辉, 童科伟, 时剑波, 申麟.  载人火星和小行星探测任务初步分析 . 深空探测学报(中英文), 2015, 2(1): 10-19. doi: 10.15982/j.issn.2095-7777.2015.01.002
    [18] 傅惠民, 娄泰山, 肖强.  火星进入段探测器自校准状态估计 . 深空探测学报(中英文), 2015, 2(3): 224-228. doi: 10.15982/j.issn.2095-7777.2015.03.006
    [19] 侯建文, 周杰.  “火星科学实验室”巡航段导航、制导与控制 . 深空探测学报(中英文), 2014, 1(2): 110-116.
    [20] 陈颖, 周璐, 王立.  一种火星多模式组合探测任务设想 . 深空探测学报(中英文), 2014, 1(2): 156-160.
  • 加载中
图(4) / 表 (5)
计量
  • 文章访问数:  217
  • HTML全文浏览量:  65
  • PDF下载量:  29
  • 被引次数: 0
出版历程
  • 收稿日期:  2017-11-23
  • 修回日期:  2018-01-18
  • 网络出版日期:  2020-08-10
  • 刊出日期:  2020-08-10

火星探测器制动捕获策略研究

doi: 10.15982/j.issn.2095-7777.2020.20171123001
    基金项目:  国家自然科学基金资助项目(61973100,61673135,61876050)
    作者简介:

    王靓玥(1996–),女,博士,主要研究方向:航天器姿轨一体化制导与控制。通讯地址:哈尔滨市西大直街92号哈尔滨工业大学327信箱主楼421室(150001)E-mail:wlyue96@163.com

    马广富(1963–),男 教授,主要研究方向:航天器制导与控制、最优控制。通讯地址:哈尔滨工业大学327信箱(150001)电话:(0451)86402726 E-mail:magf@hit.edu.cn

    通讯作者: 郭延宁(1985–),男,博士,副教授,主要研究方向:航天器控制,深空探测导航、制导与控制。本文通讯作者。通讯地址:哈尔滨市西大直街92号哈尔滨工业大学327信箱主楼421室(150001)电话:(0451)86413411-8606 E-mail:guoyn@hit.edu.cn
  • ● Three typical capture models using thrusters are systematically compared and the merits and demerits aresummarized. ● The optimal ignition point is discovered by the numerical simulation. ● A closed-loop guidance strategy based on the Zero Effort Miss/Zero Effort Velocity is proposed to implement specific transition track.
  • 中图分类号: V47

摘要: 针对火星探测器制动捕获问题进行研究,通过对火星探测器捕获问题的动力学建模、典型工况数值仿真,对3种典型的直接捕获策略:定向捕获、匀速捕获、变速捕获进行比较、分析,系统地总结了3种捕获策略的优缺点,得到变速捕获推力器燃料消耗低,但需要实时速度信息等结论。进一步地,通过数值分析点火时刻对制动捕获效果的影响,发现探测器最优点火时刻为到达近火点前制动耗时一半左右所对应的时刻。为实现特定捕获过渡轨道,设计了基于ZEM/ZEV闭环制导的捕获策略。通过与变速捕获对比表明:所设计捕获策略的制动效果与变速捕获相当,且能满足特定的末端位置/速度约束,可为实际任务提供一种备选方案。

注释:
1)  ● Three typical capture models using thrusters are systematically compared and the merits and demerits aresummarized. ● The optimal ignition point is discovered by the numerical simulation. ● A closed-loop guidance strategy based on the Zero Effort Miss/Zero Effort Velocity is proposed to implement specific transition track.

English Abstract

王靓玥, 郭延宁, 马广富. 火星探测器制动捕获策略研究[J]. 深空探测学报(中英文), 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
引用本文: 王靓玥, 郭延宁, 马广富. 火星探测器制动捕获策略研究[J]. 深空探测学报(中英文), 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
WANG Liangyue, GUO Yanning, MA Guangfu. Research About Strategies of Mars Probe Capture[J]. Journal of Deep Space Exploration, 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
Citation: WANG Liangyue, GUO Yanning, MA Guangfu. Research About Strategies of Mars Probe Capture[J]. Journal of Deep Space Exploration, 2020, 7(2): 178-183. doi: 10.15982/j.issn.2095-7777.2020.20171123001
    • 火星是距离地球最近的类地行星,自1960年起,人类就开始了火星探测的步伐,截止目前一共进行了43次火星探测任务[1]。2016年1月,中国火星探测项目正式立项。然而,已发射的火星探测器失败率高达50%,较大一部分探测器是由于未成功进行捕获制动而任务失败,如俄罗斯1973年7月发射的“火星5号”,美国1999年12月发射的“火星气候探测者”等[2]。火星制动捕获阶段器地时延大,且机会唯一,是探测器飞行全过程中最为重要的阶段。火星制动捕获是指在探测器接近火星时对探测器制动减速,使探测器能够被火星引力吸引进入火星环绕轨道,而不飞掠或撞击火星。捕获策略的选择对探测器的总体设计及探测任务的规划有直接影响[3],高性能的捕获策略是火星环绕及着陆任务开展的重要前提条件。因此分析研究火星制动捕获策略对火星探测任务的成功开展至关重要。

      火星同地球一样,具有大气,因此捕获制动可分为两类。一类为大气捕获制动,即探测器不经制动直接进入火星大气,利用大气的摩擦阻力进行减速[4]。吕静[5]等针对大气制动捕获,分别建立了在旋转大气和静止大气环境下的火星捕获模型,对任务耗时、探测器承受热量等进行比较分析。另一类为直接捕获制动,即利用探测器的推力器制动实现轨道捕获,可分为直接捕获、匀速捕获、变速捕获3类。李军锋[3]等针对有限推力模型,利用粒子群优化算法和最优控制对轨道半长轴和偏心率进行约束,对燃料最优捕获、姿态匀速转动和姿态惯性定向捕获策略即本文中的变速捕获、匀速捕获、直接捕获策略进行仿真比较,得到燃料最优捕获制动效果最好,但工程上不易实现等结论,但对3种捕获策略的比较不够全面系统。罗续盛[6]等利用多目标粒子群算法优化控制参数,以消耗燃料和入轨近火点的高度误差为双目标,对3种制动捕获策略进行分析,得到匀速捕获策略可为实际任务优先选用,且双目标提供更全面的决策信息。闵学龙[4]等研究了这两类火星制动捕获策略,并对任务耗时及速度增量进行了仿真分析比较,得到大气捕获策略速度需求小,制动耗时长,过载大,直接捕获策略速度需求大,制动耗时短等结论。由于火星大气成分复杂多变,大气捕获对探测器隔热性能要求高,制动过程耗时长,直接捕获更具有工程应用价值。

      影响火星制动捕获性能的因素较多,比如制动捕获方案的选取、推力发动机点火时刻的选取、推力发动机关机条件设计、推力发动机最大推力约束、目标过渡轨道设计等[12]。对于直接捕获制动情况,目前文献中提到的已有的制动捕获方案包括定向捕获、匀速捕获、变速捕获等3种,而目前公开发表的文献中关于这3种捕获方案的系统对比分析较少,此外,关于推力发动机的点火时刻、关机条件或过渡轨道设计等因素对制动捕获效果的影响研究也较少。

      基于此,本文建立了火星探测器制动捕获过程中简化的动力学二体模型和姿态动力学模型,针对3种典型的制动捕获方案分别给出了详细的期望姿态定义,通过设定特定远火点高度作为过渡轨道条件,在只进行一次点火制动的情况下,定量地对3种策略进行分析比较。并以变速捕获制动为例,对不同点火时刻的制动效果进行仿真比较,以探寻燃料最优的点火时刻。鉴于零控速度偏差(ZEV)和零控位移偏差/零控速度偏差(ZEM/ZEV)的制导方法已成功应用于弹道导弹拦截、小行星拦截等实际任务,并已取得成功案例[7-8],本文探索了ZEV和ZEM/ZEV制导策略在火星制导捕获中的应用。

    • 理论上,天体引力场可以达到空间无穷远处,但在一定范围内,某个天体产生的引力起主导作用。火星探测器在接近火星至被火星捕获的过程中,火星探测器受到太阳、火星、火卫及太阳系其他大质量天体的引力作用,属于多体动力学问题。由于距离等多因素的影响,其它天体对火星探测器产生的引力大小相对于火星对探测器产生的引力可以忽略,同时,考虑火星为匀质标准球体,即忽略火星的非球形引力摄动项,可将火星制动捕获问题简化为二体模型。

      探测器的轨道动力学方程如式(1)所示

      $$\left\{ {\begin{array}{*{20}{l}} {{{\dot r}} = {{v}}} \\ {{{\dot v}} = \dfrac{{{F}}}{m} - {\mu _{\rm m}}\dfrac{{ {{r}} }}{{{{{||r||}}^{{3}}}}}} \\ {\dot m = \dfrac{{{||F||}}}{{{I_{\rm sp}}g}}} \end{array}} \right.$$ (1)

      其中:${{r}}$为探测器的在火心惯性坐标系(MCI)的位置矢量;${{v}}$为探测器的速度矢量;${\mu _{\rm m}}$为火星的引力常数;${I_{\rm sp}}$为探测器推力发动机真空比冲;$m$为探测器的质量;${{F}}$为推力器产生的制动力。

    • 在火星制动捕获过程中,由于仅配置一个大推力发动机用于产生控制加速度,在施加不同方向的推力时,必须配合姿态控制系统完成。姿态跟踪控制的数学模型可由式(2)表示[9]

      $$ \left\{ {\begin{array}{*{20}{l}} {{{J}}{{{{\dot \omega }}}_{\bf{e}}} = - {{{\omega }}^ \times }{\bf{J\omega }} + {{J}}({{{\omega }}_{\bf{e}}}^ \times {{R}}{{{\omega }}_{\bf{d}}} - {{R}}{{{{\dot \omega }}}_{\bf{d}}}) + {{u}} + {{d}}} \\ {{{{{\dot q}}}_{\rm e}} = \dfrac{1}{2}\Xi ({{{q}}_{\rm e}}){\omega _{\rm e}}} \end{array}} \right. $$ (2)

      其中:$\Xi ({{{q}}_{\rm e}}) = \left[ {\begin{array}{*{20}{c}} { - {{{q}}_{\rm ev}}^{\rm T}} \\ {{{E}}({{{q}}_{\rm e}})} \end{array}} \right]$${{E}}({{{q}}_{\rm e}}) = {{{q}}_{\rm e0}}{{{I}}_3} + {{{q}}_{\rm ev}}^ \times $${{{q}}_{\rm e}} = {{{q}}_{\rm d}}^{ - 1} \otimes {{q}}$${{d}}$为干扰力矩;${{u}}$为控制力矩;${{\omega }}$为本体坐标系相对于惯性坐标系的角速度;${{{\omega }}_{\bf{d}}}$为期望坐标系相对于惯性坐标系的角速度;${{{\omega }}_{\bf{e}}}$为误差角速度;${{{q}}_{\rm d}}$为期望坐标系相对于惯性坐标系的姿态四元数;${{q}}$为探测器本体坐标系相对于惯性坐标系的姿态四元数;${{{q}}_{\rm e}}$${{q}} $${{q}}_{\rm d} $之间的误差四元数;${{R}} \in {\Re ^{3 \times {\rm{3}}}}$为本体坐标系相对于期望坐标系的转换矩阵,可由${{q}}_{\rm d} $求得。

    • 设期望坐标系三轴的单位矢量在MCI的坐标为$\left( {{{i}},{{j}},{{k}}} \right)$,探测器捕获点火时刻的初始速度矢量为${{{v}}_{\rm{0}}}$,位置矢量为${{{r}}_{\rm{0}}}$

      火心惯性坐标系到期望坐标系的转换矩阵,可以如式(3)表示

      $${{R}}_{\rm MCI} = {\left[ {\begin{array}{*{20}{c}} {{i}}&{{j}}&{{k}} \end{array}} \right]^{\rm T}}$$ (3)
    • 定向捕获方式是指探测器的指向在MCI中始终保持固定,推力发动机向一个固定的方向输出推力。

      定向捕获方式以探测器点火时刻的姿态为整个制动捕获时刻的期望姿态指向,即始终保持探测器姿态不变,即推力方向始终为制动捕获点火时刻的反方向。当探测器轨道参数中的远地点高度下降到一定阈值以后即可停止捕获。

      $\left( {{{{i}}_{\bf{d}}},{{{j}}_{\bf{d}}},{{{k}}_{\bf {d}}}} \right)$的确定方式如式(4)

      $$ {{{i}}_{\rm d}} = \dfrac{{ - {{{v}}_{\rm{0}}}}}{{\left\| {{{{v}}_{\rm{0}}}} \right\|}},\;\;{{{j}}_{\rm d}} = \dfrac{{{{{r}}_{\rm{0}}} \times {{{v}}_{\rm{0}}}}}{{\left\| {{{{r}}_{\rm{0}}} \times {{{v}}_{\rm{0}}}} \right\|}},\;\;{{{k}}_{\rm d}} = {{{i}}_{\rm d}} \times {{{j}}_{\rm d}} $$ (4)

      其中:${{{i}}_{\bf{d}}}$方向为速度反方向;${{{j}}_{\bf{d}}}$方向为轨道法向;${{{k}}_{\bf{d}}}$方向与${{{i}}_{\bf{d}}}$${{{j}}_{\bf{d}}}$成右手系。

    • 由于定向捕获方式推力方向与速度方向不始终在同一直线上,只有沿速度矢量方向的推力分量才会对探测器产生阻尼制动作用,推力利用率低。为了提高推力的利用率,将探测器的捕获轨道近似为圆轨道,使探测器在轨道平面内匀速旋转,旋转角速度的选取如文献[10]所示。停止捕获条件同定向捕获。

      $\left( {{{{i}}_{\bf{c}}},{{{j}}_{\bf{c}}},{{{k}}_{\bf{c}}}} \right)$的确定方式如式(5)

      $$\left[ {\begin{array}{*{20}{c}} {{{{i}}_{\rm c}}}&{{{{j}}_{\rm c}}}&{{{{k}}_{\rm c}}} \end{array}} \right] = {{{R}}_2}({{{\omega }}_{\rm Mean}}\Delta t)\left[ {\begin{array}{*{20}{c}} {{{{i}}_{\rm d}}}&{{{{j}}_{\rm d}}}&{{{{k}}_{\rm d}}} \end{array}} \right]$$ (5)

      其中:${{{\omega }}_{\rm Mean}}$为探测器匀速变化的角速度;$\Delta t$为制动捕获过程开始后经历的时间;${{{R}}_2}$为期望坐标系匀速旋转的转换矩阵。

    • 为了最大限度地利用推力,主推力器在空间中输出推力方向应始终沿着速度矢量的反方向。变速捕获方式需要实时测量探测器的速度矢量,制动捕获过程的期望姿态应随着速度矢量的改变而改变。停止捕获条件同定向捕获。

      $\left( {{{{i}}_{\bf{t}}},{{{j}}_{\bf{t}}},{{{k}}_{\bf{t}}}} \right)$的求法如下

      $$ {{{i}}_{\rm t}} = \dfrac{{ - {{v}}}}{{\left\| {{v}} \right\|}}\;\;\;\;{{{j}}_{\rm t}} = \dfrac{{{{r}} \times {{v}}}}{{\left\| {{{r}} \times {{v}}} \right\|}}\;\;\;\;{{{k}}_{\rm t}} = {{{i}}_{\rm t}} \times {{{j}}_{\rm t}} $$ (6)
    • 取制动捕获点火时刻为探测器距离火星位置最近的时刻,其速度矢量、位置矢量由进入火星影响球半径的速度矢量、位置矢量自由运动求得,仿真条件参照罗续盛[6]等的论文,如表1所示。

      表 1  仿真条件

      Table 1.  Simulation conditions

      项目参数
      探测器质量/kg4.948 × 103
      发动机推力/N 3 × 103
      比冲/s 314
      点火时刻速度矢量/(km·s–1[–1.980;4.416;–2.544]
      点火时刻位置矢量/km[3 330.9;1 981.4;829.8]

      通过Matlab软件对公式(4)~(6)对应的3种捕获模式进行数值仿真,3种捕获模式下探测器的运行轨迹如图1所示,轨道平面如图2所示。3种捕获模式的仿真结果及优缺点如表2所示。

      图  1  探测器运行轨迹图

      Figure 1.  Trajectory of probe

      图1可以看出,定向捕获探测器的运动轨迹最长,间接反映了该捕获模式耗时最长。主推力器开启后,燃料消耗与捕获耗时成正比,故定向捕获燃料消耗最多。从表2图2可得,当匀速捕获角速度选取较为合适时,匀速捕获与变速捕获的制动效果相当。

      图  2  探测器运行的轨道平面图

      Figure 2.  Orbital plane of probe running

      表 2  仿真结果及优缺点

      Table 2.  Simulation result, merits and demerits

      项目定向捕获匀速捕获变速捕获
      捕获时间/s1 446.81 360.01 358.8
      消耗燃料/kg1 410.01 325.11 323.9
      距火星表面
      最近距离/km
      635.425291.875247.065
      参数个数565
      优点实现简单实现简单,推力器利用率高推力器利用率高
      缺点推力器利用率低制动效果严重依赖于角速度的选取推力器实现复杂

      定向捕获推力发动机在空间中向一个固定方向输出推力,捕获过程实现简单,但探测器的速度矢量实时变化,其推力器利用率低,制动捕获过程耗时长,燃料消耗多,一般在实际任务中较少使用。

      匀速捕获在角速度选取合理时,耗时短,消耗燃料少,制动效果与变速捕获相当,但其制动效果的优劣严重依赖于角速度的选取。一般在实际任务中常作为备选方案使用。

      变速捕获的制动效果最优,推力发动机产生制动力方向始终沿着速度矢量的反方向,推力器燃料消耗低,但需要实时速度信息反馈,制动效果受反馈信息精度影响。在实际任务中一般选择变速捕获方式,如欧洲航天局(European Space Agency,ESA)发射的“火星快车号”(Mars Express)[11]

      3种捕获策略的近火点的高度均在200 km以上,有足够的余量保证探测器因误差原因而受到火星大气摩擦,导致不可控影响或因火星表面不平撞到火星。

    • 以2.4节选择变速捕获策略为例,探测器的点火时刻为近火点时刻。为了定量的比较不同点火时刻的制动捕获效果,其中考虑近火点时刻为0 s,仿真可得不同点火时刻对应的捕获过程耗时和耗燃料曲线如图3所示,不同点火时刻对应的捕获轨道中距火星表面的最近和最远距离如图4所示。

      图  3  不同点火时刻耗时、耗燃料

      Figure 3.  Time and fuel consumption of ignition time

      图  4  不同点火时刻距火星最远、近距离

      Figure 4.  The furthest and nearest distance from Mars of ignition time

      图3~4仿真数据可得,点火的最优时刻在探测器到达近火点前688 s,此时的制动捕获时间为1 268.4 s,消耗燃料为1 235.9 kg。燃料消耗较近火点时刻点火减少6.7%左右,制动耗时减少6.7%左右,捕获后轨道近火点最高。探测器的点火时刻不宜在探测器到达近火点后,其制动捕获的耗时过长,消耗燃料过多,且制动捕获后轨道的近火点高度过低。

      为推广得到更为一般性的结论,改变2.4节中的仿真条件,对各个工况下的制动捕获进行打靶仿真,仿真条件及结果如表3所示,不同点火时刻对应的捕获过程耗时和耗燃料曲及不同点火时刻对应的捕获轨道中的最近和最远距离曲线由于篇幅原因不详细给出。

      表 3  不同工况仿真条件及结果

      Table 3.  Simulation results of different conditions

      发动机推力/
      N
      比冲/
      s
      最优点火刻/
      s
      制动耗时/
      s
      燃料消耗/
      kg
      3 000314–6881 268.41 235.9
      3 500314–5721 082.41 230.4
      3 000250–6691 246.81 363.1
       注:探测器到达近火点时刻为0 s。
    • 由仿真示例推广可得,探测器燃料最优点火时刻为探测器到达近火点前制动耗时一半左右对应的时刻。点火时刻不宜选择探测器到达近火点后,其制动捕获的耗时过长,消耗燃料过多,不利于探测器的设计与控制,且制动捕获后轨道的近火点高度过低,探测器有撞到火星表面的危险。

    • 零控速度偏差(ZEV)、零控位移偏差(ZEM)最早应用于制导技术,在弹道导弹拦截,小行星拦截等工程应用中已取得一定成功。本文将ZEV、ZEM/ZEV制导策略应用于火星捕获任务,设计基于ZEV、ZEM/ZEV闭环制导的制导捕获策略,以实现特定捕获过渡轨道。

    • 给定制动捕获任务初末时刻${t_0}$${t_f}$,以及制动捕获点速度矢量${{{r}}_0}$,位置矢量${{{v}}_0}$,期望时刻的速度矢量${{{v}}_f}$,考虑基于加速度平方积分形式的性能指标如式(7)所示。

      $$J = \dfrac{1}{2}\int_{{t_0}}^{{t_f}} {{{{a}}^{\rm T}}{{a}}} {\rm d}t$$ (7)

      ${{ZEV}}$指在不控情况下,探测器的期望速度矢量与由当前时刻自由运行到期望时刻的速度矢量之差。${{ZEV}}$表征当前状态同期望状态的速度差。可由式(8)表示

      $${{ZEV}} = {{{v}}_f} - {{{\tilde v}}_f}$$ (8)

      其中:${{{\tilde v}}_f}$为期望时刻的位置矢量,利用数值积分预测在${t_f}$时刻的${{{\tilde v}}_f}$

      利用性能指标及约束条件得到最优加速度如式(9)

      $${{a}} = \dfrac{{{{ZEV}}}}{{{t_{\rm go}}}}$$ (9)

      其中:${t_{\rm go}} = {t_f} - t$为剩余飞行时间。

    • 给定制动捕获任务初末时刻${t_0}$${t_f}$,以及制动捕获点火时刻速度矢量${{{r}}_0}$,位置矢量${{{v}}_0}$,期望时刻的速度矢量${{{v}}_f}$,位置矢量${{{r}}_f}$,考虑基于加速度平方积分形式的性能指标如式(7)所示。

      ${{{{ZEM}}} / {{{ZEV}}}}$指在不控情况下,探测器的期望速度矢量、位移矢量与由当前时刻自由运行到期望时刻的速度、位移矢量之差。${{ZEV}}$由式(8)所示,${{ZEM}}$可由式(10)表示

      $${{ZEM}} = {{{r}}_f} - {{{\tilde r}}_f}$$ (10)

      其中:${{{\tilde r}}_f}$为期望时刻的位置矢量,利用数值积分预测在${t_f}$时刻的${{{\tilde r}}_f}$

      由文献[7~8]可得,最优控制律表达式如式(11)所示

      $$ {{a}} = \dfrac{6}{{{t_{\rm go}}^2}}{{ZEV}} - \dfrac{2}{{{t_{\rm go}}}}{{ZEM}} $$ (11)
    • 以变速捕获制动结束时的速度、位置矢量作为期望时刻的速度、位置矢量,取2.4节的仿真结果作为本次仿真条件,即制动捕获时间为1 358.8 s,初始及期望速度、位置矢量如表4所示。仿真结果如表5所示,仿真结果图由于篇幅原因不详细给出。

      表 4  ZEV、ZEM/ZEV仿真条件

      Table 4.  Simulation conditions of ZEV and ZEM/ZEV

      时间速度矢量/(km·s–1位移矢量/(km)
      制动捕获开始 [–1.980;4.416;–2.544] [3330.9;1981.4;829.8]
      制动捕获结束[–2.558;1.726;–1.828] [–189.9;6010;–2292]

      表 5  新型捕获策略仿真结果

      Table 5.  Simulation results of the new capture strategy

      捕获
      策略
      捕获
      时间/s
      消耗
      燃料/kg
      终点位置/km终点速度/(km·s–1
      变速
      捕获
      1 358.81 324[–189.9;6 010;
      –2 292]
      [–2.558;1.726;
      –1.828]
      ZEV1 3591 324[–140.2;6 012;
      –2 270]
      [–2.558;1.728;
      –1.828]
      ZEM/
      ZEV
      1 3591 323[–189.6;6 010;
      –2 292]
      [–2.559;1.726;
      –1.828]

      由仿真结果可以看出ZEV和ZEM/ZEV制动捕获策略可以获得和变速捕获相当的制动效果,并且能满足特定的末端位置/速度约束。

    • 本文首先对火星探测器捕获问题进行动力学建模,并阐述了3种典型捕获策略的原理,通过典型数值仿真进行分析比较,系统地总结了3种捕获策略的优缺点,得到定向捕获策略控制简单、实现容易但燃料消耗高;匀速捕获燃料消耗低但制动捕获效果严重依赖角速度的选取;变速捕获燃料消耗低但需要实时速度信息的结论。进一步地,通过数值分析不同的点火时刻,得到结论:探测器自点火点飞行至近火点消耗时间为制动捕获总时间一半时对应的点火时刻为最优点火时刻。同时,当期望时刻位置矢量、速度矢量已知时,设计的基于ZEM/ZEV的闭环制导捕获策略不仅能够获得与变速捕获相当的捕获效果,而且能够满足特定的末端位置/速度约束,可作为实际任务的一种备选方案。此外,未来可以分析最佳点火距离或设计新型制动捕获策略以期为实际探火任务提供理论支撑。

参考文献 (12)

目录

    /

    返回文章
    返回