CN113311854A - 一种月球采样返回任务中的定点着陆轨道设计方法 - Google Patents
一种月球采样返回任务中的定点着陆轨道设计方法 Download PDFInfo
- Publication number
- CN113311854A CN113311854A CN202110545890.4A CN202110545890A CN113311854A CN 113311854 A CN113311854 A CN 113311854A CN 202110545890 A CN202110545890 A CN 202110545890A CN 113311854 A CN113311854 A CN 113311854A
- Authority
- CN
- China
- Prior art keywords
- moon
- lunar
- point
- orbit
- braking
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005070 sampling Methods 0.000 title claims abstract description 24
- 238000013461 design Methods 0.000 title abstract description 12
- 238000012546 transfer Methods 0.000 claims abstract description 69
- 238000012937 correction Methods 0.000 claims abstract description 51
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims abstract description 50
- 238000004364 calculation method Methods 0.000 claims description 50
- 230000009467 reduction Effects 0.000 claims description 6
- 238000005457 optimization Methods 0.000 abstract description 8
- 238000003032 molecular docking Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 12
- 238000004422 calculation algorithm Methods 0.000 description 8
- 238000004590 computer program Methods 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 238000012986 modification Methods 0.000 description 5
- 230000004048 modification Effects 0.000 description 5
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000001174 ascending effect Effects 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 230000008685 targeting Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 210000001503 joint Anatomy 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0833—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using limited authority control
Landscapes
- Engineering & Computer Science (AREA)
- Computer Security & Cryptography (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种月球采样返回任务中的定点着陆轨道设计方法,该方法包括:结合月球采样返回任务地月转移段近月点倾角优化,通过在近月制动机动中引入轨道面法向速度增量分量来瞄准预定月面落点经度,通过环月降轨变轨纬度幅角和轨道面内速度增量分量来瞄准预定月面落点纬度和动力下降开始点高度,并使得月面落点在预定月面上升时刻位于目标轨道面内;另外,通过引入环月降轨速度增量法向分量来修正落点经度和调整动力下降时刻来消除纵向航程误差。本发明解决了现有定点着陆技术中没有考虑轨道设计对后续月球轨道交会对接的影响,导致较大的月球轨道交会对接轨道面修正速度增量的问题。
Description
技术领域
本发明涉及航天技术领域,尤其涉及一种月球采样返回任务中的定点着陆轨道设计方法。
背景技术
在月球采样返回任务或载人登月任务中,都可能涉及探测器从环月轨道开始下降,通过变轨机动进入月面动力下降弹道,并在预定月面落点实现软着陆。月面定点着陆轨道设计是我国实现月球采样返回和载人登月必须突破的一项关键轨道设计技术。
月球采样返回任务轨道全飞行过程包括运载发射、地月转移、近月制动、环月飞行、月面动力下降、月面工作、月面动力上升、交会对接、月地转移、再入回收等飞行阶段。参考图1和图2所示,探测器在完成地月转移和近月制动后,最初在一个圆形轨道上飞行,稍后轨道器和着陆器上升器组合体实现分离。着陆器上升器组合体以下简称着上组合体。着上组合体在轨道上运行几圈后,在一个合适的纬度幅角实施切向变轨机动,进入近月点高度约为15km的椭圆轨道。着上组合体将在该椭圆轨道上运行几圈,当着上组合体到达动力下降开始点并取得预定的高度15km,开始沿动力下降轨迹实现在月球表面着陆。这里动力下降开始点为近月点。
国内外诸多学者已对月面定点着陆轨道设计方法进行了深入的研究,包括通过落月圈次的调整以及环月轨道面的调整等方法来实现对预定落点经度的瞄准,还包括通过环月降轨平面内速度增量和变轨位置来瞄准预定月面落点纬度和动力下降点高度。其中,环月轨道面的调整可以通过调整升交点赤经和修正轨道倾角来实现。
在实际工程应用特别是在月球采样返回任务工程设计中,遇到定点着陆理论研究没有回答的新问题。在采样返回任务中,着上组合体落月后在一定工作时间内,例如在1~2天内进行采样工作,然后上升器必须在预定月面上升时间起飞,使得其入轨后的环月轨道面与轨道器的轨道面基本一致,以减小后续二者实现月球轨道交会对接所需的平面修正速度增量。为实现这一目标,需要在地月转移段对近月点的目标环月轨道倾角进行调整优化,使得月面起飞时上升器位于轨道器轨道面内。但倾角的调整涉及轨道面的变化,会影响月面落点位置。另一方面,现有定点着陆技术没有考虑倾角优化的影响,不能保证月面上升时月面落点位于轨道器轨道面内,所以可能导致较大的月球轨道交会对接轨道面修正速度增量。而且,现有方法没有考虑两次经度瞄准,仅采用环月降轨一次控制来瞄准预定落点经度,不能实现大的落点经度调整。所以,采样返回任务中实现月面定点着陆时不能简单套用现有技术,其中如何保证月面起飞时上升器位于轨道器轨道面内是一大难题。
发明内容
本发明的技术解决问题是:针对月球采样返回任务中要求月面起飞时上升器位于轨道器轨道面内导致的环月轨道倾角调整对月面落点的影响,对环月轨道倾角优化和月面定点着陆变轨方法进行联合设计,从而同时满足月面定点着落和月面共面起飞的设计要求,并尽量减少实际飞行中误差因素对月面落点散布的影响,保障月球采样返回任务目标的顺利完成。
第一方面,本申请实施例提供的一种月球采样返回任务中的定点着陆轨道设计方法,该方法包括:
步骤1.在每一地月转移中途修正中,采用变轨速度增量三分量瞄准预定近月点倾角、近月点高度和真近点角,其中,所述真近点角的取值为0度;
步骤2.对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度和轨道器轨道参数,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角;其中,在进行近月制动计算时,引入近月制动速度增量法向分量,其初值为0或根据公式计算,对应的有限推力计算中为引入推力偏航角,偏航角初值为0;
步骤3.调整地月转移近月点倾角的目标值和近月制动速度增量法向分量,重复步骤1和步骤2,并应用2对2微分修正计算来瞄准预定月面落点经度、所述月面上升时月面落点月心位置矢量和所述轨道器轨道法向的夹角,所述夹角的取值为90度;
步骤4.完成地月转移后,实施近月制动时,通过1对1微分修正来调整近月制动速度增量法向分量,瞄准预定月面落点经度;
步骤5.完成近月制动后,实施环月降轨时,增加速度增量法向分量,对应的有限推力计算中为引入推力偏航角,并联合变轨纬度幅角和变轨速度增量横向分量,或联合变轨速度增量横向和径向分量,通过3对3微分修正来瞄准预定落点经度、落点纬度和动力下降起始点高度;
步骤6.完成环月降轨后,在实施动力下降前,在预报的近月点时刻±1分钟内对动力下降开始点时刻进行线性调整,对每个调整后的动力下降开始点时刻数值,进行动力下降试算,直到计算所得月面落点离目标月面落点月面距离取得最小值,并根据最后一次调整后的动力下降开始点时刻执行动力下降,实现月面软着陆。
可选地,对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度和轨道器轨道参数,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角,包括:
若设计近月制动次数为两次,则两次近月制动的法向速度增量或推力方向偏航角取值相同;在进行环月降轨计算时,采用变轨位置纬度幅角和变轨速度增量横向分量或速度反方向的变轨速度增量大小,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定的落点纬度和动力下降起始点高度,其中,动力下降起始点高度的取值为15km。
可选地,完成地月转移后,实施近月制动时,通过1对1微分修正来调整近月制动速度增量法向分量,瞄准预定月面落点经度,包括:
对于近月制动速度增量法向分量的每次迭代取值,依次计算近月制动、环月降轨、动力下降轨道,得到月面落点经度;其中,在进行环月降轨计算时,采用变轨纬度幅角和变轨速度增量横向分量,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定落点纬度和动力下降起始点高度。
可选地,完成近月制动后,实施环月降轨时,增加速度增量法向分量,对应的有限推力计算中为引入推力偏航角,包括:
如果环月降轨包含两次变轨,则两次变轨的推力方向偏航角取值相同。
与现有技术相比,本发明具有如下有益效果:
(1)本发明首次解决了采样返回任务中实现月面定点着陆时如何保证月面起飞时上升器位于轨道器轨道面的难题,通过设计密切结合月面采样任务要求的近月点倾角优化的月面定点着陆轨道,以便实现高精度月面定点着陆,同时采用倾角优化来保证月面上升时月面落点位于轨道器轨道面内,从而能够最小化月球轨道交会对接时的轨道面修正速度增量,保障任务目标的实现。
(2)本发明方法中创新性地引入了内外两层微分修正嵌套计算结构。其中,在进行地月转移中途修正规划时,外层采用2对2的微分修正:采用近月点倾角和近月制动速度增量法向分量瞄准预定月面落点经度、月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角90度,内层采用2对2的微分修正:采用变轨纬度幅角和变轨速度增量横向分量瞄准落点纬度和动力下降起始点高度。数值仿真表明,该算法计算效率高,收敛速度快;实际飞控中已成功应用该算法,在短时间内完成了变轨策略生成,保障了任务的顺利开展。
(3)本发明方法充分考虑了实际飞行控制中误差因素的影响,创新性地提出了分两步来瞄准预定月面落点经度,即步骤1~4首先通过近月制动时引入速度增量法向分量瞄准预定月面落点经度,简称为“粗瞄准”,步骤5再在环月降轨时引入速度增量法向分量修正近月制动控制误差,进一步实现对预定月面落点经度的“精瞄准”。这种通过两次控制瞄准预定落点经度的设计有助于尽可能减小误差因素造成的落点经度方向散布。
附图说明
图1为本发明提供的一种地月转移、近月制动和环月降轨轨道方案示意图;
图2为本发明提供的一种环月轨道在月面投影和预定着陆点的示意图;
图3为本发明提供的一种月球采样返回任务中的定点着陆轨道设计方法的流程示意图;
图4为本发明提供的一种采样返回任务全飞行过程以及月面定点着陆和环月轨道倾角优化联合设计算法的示意图。
具体实施方式
本申请实施例提供的方案中,所描述的实施例仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本申请保护的范围。
以下结合说明书附图对本申请实施例所提供的一种月球采样返回任务中的定点着陆轨道设计方法做进一步详细的说明,该方法具体实现方式可以包括以下步骤(方法流程如图3所示):
步骤1.在每一地月转移中途修正中,采用变轨速度增量三分量瞄准预定近月点倾角、近月点高度和真近点角,其中,所述真近点角的取值为0度;
步骤2.对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度和轨道器轨道参数,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角;其中,在进行近月制动计算时,引入近月制动速度增量法向分量,其初值为0或根据公式计算,对应的有限推力计算中为引入推力偏航角,偏航角初值为0;
步骤3.调整地月转移近月点倾角的目标值和近月制动速度增量法向分量,重复步骤1和步骤2,并应用2对2微分修正计算来瞄准预定月面落点经度、所述月面上升时月面落点月心位置矢量和所述轨道器轨道法向的夹角,所述夹角的取值为90度;
步骤4.完成地月转移后,实施近月制动时,通过1对1微分修正来调整近月制动速度增量法向分量,瞄准预定月面落点经度;
步骤5.完成近月制动后,实施环月降轨时,增加速度增量法向分量,对应的有限推力计算中为引入推力偏航角,并联合变轨纬度幅角和变轨速度增量横向分量,或联合变轨速度增量横向和径向分量,通过3对3微分修正来瞄准预定落点经度、落点纬度和动力下降起始点高度;
步骤6.完成环月降轨后,在实施动力下降前,在预报的近月点时刻±1分钟内对动力下降开始点时刻进行线性调整,对每个调整后的动力下降开始点时刻数值,进行动力下降试算,直到计算所得月面落点离目标月面落点月面距离取得最小值,并根据最后一次调整后的动力下降开始点时刻执行动力下降,实现月面软着陆。
在一种可能实现的方式中,对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度和轨道器轨道参数,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角,包括:
若设计近月制动次数为两次,则两次近月制动的法向速度增量或推力方向偏航角取值相同;在进行环月降轨计算时,采用变轨位置纬度幅角和变轨速度增量横向分量或速度反方向的变轨速度增量大小,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定的落点纬度和动力下降起始点高度,其中,动力下降起始点高度的取值为15km。
在一种可能实现的方式中,完成地月转移后,实施近月制动时,通过1对1微分修正来调整近月制动速度增量法向分量,瞄准预定月面落点经度,包括:
对于近月制动速度增量法向分量的每次迭代取值,依次计算近月制动、环月降轨、动力下降轨道,得到月面落点经度;其中,在进行环月降轨计算时,采用变轨纬度幅角和变轨速度增量横向分量,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定落点纬度和动力下降起始点高度。
在一种可能实现的方式中,完成近月制动后,实施环月降轨时,增加速度增量法向分量,对应的有限推力计算中为引入推力偏航角,包括:
如果环月降轨包含两次变轨,则两次变轨的推力方向偏航角取值相同。
具体地,如图4所示,本发明的月球采样返回任务中的定点着陆轨道设计方法,包括内外两层微分修正的嵌套计算结构,具体算法和步骤如下所述:
1)、首先,在地月转移中途修正中,进行环月轨道倾角优化。
式中,ωm是月球的自转角速率。根据上式计算可得倾角初值。
对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角。
其中,在进行近月制动计算时,引入速度增量法向分量Δvn,其初值为0或根据下述公式计算,对应的有限推力计算中为引入推力偏航角,偏航角初值取为0。具体的,推力偏航角计算公式如下所示:
式中,Δλ是需要进行的月面落点经度调整量,r为变轨时的月心距,n为轨道角速率,a为轨道半长轴,e为轨道偏心率,u为轨道纬度幅角,Ω为升交点赤经。
另外,在进行环月降轨计算时,采用变轨位置纬度幅角和变轨速度增量横向分量或速度反方向的速度增量大小,或采用速度增量的横向和径向分量,通过微分修正计算来瞄准落点纬度和动力下降起始点高度15km。
重复上述步骤来调整地月转移近月点倾角的目标值和近月制动速度增量法向分量,即应用2对2微分修正计算来瞄准预定月面落点经度、月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角90度。
上述计算过程为一个内外两层微分修正嵌套计算结构,外层采用2对2的微分修正,即采用近月点倾角和近月制动法向分量瞄准预定月面落点经度、月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角90度,内层采用2对2的微分修正,即采用变轨纬度幅角和变轨速度增量横向分量瞄准预定落点纬度和动力下降起始点高度。
2)、其次,在近月制动实施中,调整速度增量法向分量,实现预定月面落点经度的“粗瞄准”。
具体地,完成地月转移后,实施近月制动时,通过1对1微分修正来调整速度增量法向分量,瞄准预定月面落点经度。该计算过程中,对于近月制动速度增量法向分量的每次迭代取值,依次计算近月制动、环月降轨、动力下降轨道,得到月面落点经纬度。
在进行环月降轨计算时,采用变轨纬度幅角和变轨速度增量横向分量,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定落点纬度和动力下降起始点高度。
上述计算过程为一个内外两层微分修正嵌套计算结构,外层采用1对1的微分修正,即采用近月制动速度增量法向分量瞄准预定月面落点经度,内层采用2对2的微分修正,即采用变轨纬度幅角和变轨速度增量横向分量瞄准预定落点纬度和动力下降起始点高度。
3)、在实施环月降轨时,瞄准月面落点纬度,同时进行月面落点经度的“精瞄准”。
具体来说,完成近月制动后,实施环月降轨时,增加速度增量法向分量,联合变轨纬度幅角和变轨速度增量横向分量,或联合速度增量横向和径向分量,通过3对3微分修正来瞄准预定落点经度、落点纬度和动力下降起始点高度。
其中,速度增量法向分量的初值取为0,或采用前面月面落点经度调整公式估算。
4)、最后,进行动力下降开始点时刻调整。
完成环月降轨后,在实施动力下降前,在预报的近月点时刻±1分钟内对动力下降开始点时刻进行线性调整,对调整后的每个动力下降开始点时刻的数值,进行动力下降试算,直到计算所得月面落点离目标月面落点的月面距离取得最小值。
这里所说的动力下降开始点时刻线性调整,可以采用拉网格的方法进行搜索。
月面距离的计算可采用球面三角公式:
需要补充说明的是:(1)在图3所示的算法示意图中,月面工作时间是通过指定月面上升绝对时刻,并取月面工作前月心对国内深空站可见后首次到达近月点为动力下降点来保证的;如果要求定时着陆,则通过调整近月制动后环月轨道的半长轴或轨道周期的目标值来实现;(2)在图3所示的算法示意图中,采用的是外层2对2微分修正、内层2对2微分修正的嵌套计算结构,或者采用其它变通形式的嵌套计算结构,例如:外层1对1微分修正进行倾角优化、中层1对1微分修正来调整近月制动速度增量法向分量以便瞄准预定月面落点经度、内层2对2微分修正为采用环月降轨纬度幅角和变轨速度增量横向分量瞄准预定落点纬度和动力下降起始点高度,或者外层1对1微分修正进行倾角优化、内层3对3微分修正为采用近月制动速度增量法向分量、环月降轨纬度幅角和变轨速度增量横向分量瞄准预定月面落点经度、纬度和动力下降起始点高度;(3)在实施近月制动时,外层采用环月降轨速度增量法向分量,和近月制动速度增量法向分量一起瞄准预定月面落点经度、月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角90度。
本发明说明书中未作详细描述的内容属本领域技术人员的公知技术。例如:地月转移、近月制动和环月降轨的过程,近月制动轨道面内速度增量一般用于瞄准控后环月轨道的半长轴或轨道周期,轨道器调相采用多脉冲变轨瞄准月球轨道交会对接交班点的纬度幅角、半长轴和0偏心率,以及1对1、2对2、3对3的微分修正算法等。
综上所述,一方面本发明首次解决了采样返回任务中实现月面定点着陆时如何保证月面起飞时上升器位于轨道器轨道面的难题,通过设计密切结合月面采样任务要求的近月点倾角优化的月面定点着陆轨道,为实现高精度月面定点着陆创造了必要条件,同时倾角优化可以保证月面上升时月面落点位于轨道器轨道面内,从而能够最小化月球轨道交会对接时的轨道面修正速度增量,对保障任务目标的实现具有重要意义。
另一方面,本发明方法中创新性地引入了内外两层微分修正嵌套计算结构。其中,在进行地月转移中途修正规划时,外层采用2对2的微分修正:采用近月点倾角和近月制动速度增量法向分量瞄准预定月面落点经度、月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角90度,内层采用2对2的微分修正:采用变轨纬度幅角和变轨速度增量横向分量瞄准预定落点纬度和动力下降起始点高度。数值仿真表明,该算法计算效率高,收敛速度快;实际飞控中已成功应用该算法,在短时间内完成了变轨策略生成,保障了任务的顺利开展。
第三方面,本发明方法充分考虑了实际飞行控制中误差因素的影响,创新性地提出了分两步来瞄准预定月面落点经度,即步骤1~4首先通过近月制动时引入速度增量法向分量瞄准预定月面落点经度,简称“粗瞄准”,步骤5再在环月降轨时引入速度增量法向分量修正近月制动控制误差,进一步实现对预定月面落点经度的“精瞄准”。这种通过两次控制瞄准预定落点经度的设计有助于尽可能减小误差因素造成的落点经度方向散布。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器和光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (4)
1.一种月球采样返回任务中的定点着陆轨道设计方法,其特征在于,包括:
步骤1.在每一地月转移中途修正中,采用变轨速度增量三分量瞄准预定近月点倾角、近月点高度和真近点角,其中,所述真近点角的取值为0度;
步骤2.对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度和轨道器轨道参数,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角;其中,在进行近月制动计算时,引入近月制动速度增量法向分量,其初值为0或根据公式计算,对应的有限推力计算中为引入推力偏航角,偏航角初值为0;
步骤3.调整地月转移近月点倾角的目标值和近月制动速度增量法向分量,重复步骤1和步骤2,并应用2对2微分修正计算来瞄准预定月面落点经度、所述月面上升时月面落点月心位置矢量和所述轨道器轨道法向的夹角,所述夹角的取值为90度;
步骤4.完成地月转移后,实施近月制动时,通过1对1微分修正来调整近月制动速度增量法向分量,瞄准预定月面落点经度;
步骤5.完成近月制动后,实施环月降轨时,增加速度增量法向分量,对应的有限推力计算中为引入推力偏航角,并联合变轨纬度幅角和变轨速度增量横向分量,或联合变轨速度增量横向和径向分量,通过3对3微分修正来瞄准预定落点经度、落点纬度和动力下降起始点高度;
步骤6.完成环月降轨后,在实施动力下降前,在预报的近月点时刻±1分钟内对动力下降开始点时刻进行线性调整,对每个调整后的动力下降开始点时刻数值,进行动力下降试算,直到计算所得月面落点离目标月面落点月面距离取得最小值,并根据最后一次调整后的动力下降开始点时刻执行动力下降,实现月面软着陆。
2.如权利要求1所述的方法,其特征在于,对于一个具体的近月点倾角数值,依次计算地月转移、近月制动、轨道器调相、环月降轨、动力下降轨道,得到月面落点经纬度和轨道器轨道参数,并计算在预定时间实施月面上升时月面落点月心位置矢量和轨道器轨道法向的夹角,包括:
若设计近月制动次数为两次,则两次近月制动的法向速度增量或推力方向偏航角取值相同;在进行环月降轨计算时,采用变轨位置纬度幅角和变轨速度增量横向分量或速度反方向的变轨速度增量大小,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定的落点纬度和动力下降起始点高度,其中,动力下降起始点高度的取值为15km。
3.如权利要求1所述的方法,其特征在于,完成地月转移后,实施近月制动时,通过1对1微分修正来调整近月制动速度增量法向分量,瞄准预定月面落点经度,包括:
对于近月制动速度增量法向分量的每次迭代取值,依次计算近月制动、环月降轨、动力下降轨道,得到月面落点经度;其中,在进行环月降轨计算时,采用变轨纬度幅角和变轨速度增量横向分量,或采用速度增量横向和径向分量,通过微分修正计算来瞄准预定落点纬度和动力下降起始点高度。
4.如权利要求1所述的方法,其特征在于,完成近月制动后,实施环月降轨时,增加速度增量法向分量,对应的有限推力计算中为引入推力偏航角,包括:
如果环月降轨包含两次变轨,则两次变轨的推力方向偏航角取值相同。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110545890.4A CN113311854A (zh) | 2021-05-19 | 2021-05-19 | 一种月球采样返回任务中的定点着陆轨道设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110545890.4A CN113311854A (zh) | 2021-05-19 | 2021-05-19 | 一种月球采样返回任务中的定点着陆轨道设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113311854A true CN113311854A (zh) | 2021-08-27 |
Family
ID=77373878
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110545890.4A Pending CN113311854A (zh) | 2021-05-19 | 2021-05-19 | 一种月球采样返回任务中的定点着陆轨道设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113311854A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114115330A (zh) * | 2021-10-29 | 2022-03-01 | 北京空间飞行器总体设计部 | 一种兼顾火星环绕、进入、着陆探测的轨道设计方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109344449A (zh) * | 2018-09-07 | 2019-02-15 | 北京空间技术研制试验中心 | 航天器月地转移轨道逆向设计方法 |
US10289761B1 (en) * | 2013-06-12 | 2019-05-14 | The United States Of America, As Represented By The Secretary Of The Navy | Method for modeling dynamic trajectories of guided, self-propelled moving bodies |
CN109992823A (zh) * | 2019-02-15 | 2019-07-09 | 北京空间飞行器总体设计部 | 一种月面高精度定时定点着陆轨道控制方法 |
CN110765504A (zh) * | 2019-10-29 | 2020-02-07 | 北京空间技术研制试验中心 | 航天器环月轨道交会对接的轨道设计方法 |
-
2021
- 2021-05-19 CN CN202110545890.4A patent/CN113311854A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10289761B1 (en) * | 2013-06-12 | 2019-05-14 | The United States Of America, As Represented By The Secretary Of The Navy | Method for modeling dynamic trajectories of guided, self-propelled moving bodies |
CN109344449A (zh) * | 2018-09-07 | 2019-02-15 | 北京空间技术研制试验中心 | 航天器月地转移轨道逆向设计方法 |
CN109992823A (zh) * | 2019-02-15 | 2019-07-09 | 北京空间飞行器总体设计部 | 一种月面高精度定时定点着陆轨道控制方法 |
CN110765504A (zh) * | 2019-10-29 | 2020-02-07 | 北京空间技术研制试验中心 | 航天器环月轨道交会对接的轨道设计方法 |
Non-Patent Citations (2)
Title |
---|
李革非 等: "定时定点月面着陆全程轨道控制设计", 《宇航学报》 * |
汪中生 等: "月面定点着陆变轨策略研究", 《航天器工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114115330A (zh) * | 2021-10-29 | 2022-03-01 | 北京空间飞行器总体设计部 | 一种兼顾火星环绕、进入、着陆探测的轨道设计方法 |
CN114115330B (zh) * | 2021-10-29 | 2022-07-29 | 北京空间飞行器总体设计部 | 一种兼顾火星环绕、进入、着陆探测的轨道设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109398762B (zh) | 一种固体火箭基于椭圆转移轨道的入轨式弹道设计方法 | |
CN111591469B (zh) | 一种低轨星座系统相位保持方法、系统、设备和存储介质 | |
CN109625323B (zh) | 一种卫星化学推进变轨方法及系统 | |
CN109911249B (zh) | 低推重比飞行器的星际转移有限推力入轨迭代制导方法 | |
CN110104219B (zh) | 一种控制探测器着陆地外天体的方法及装置 | |
CN110765504B (zh) | 航天器环月轨道交会对接的轨道设计方法 | |
CN109080854B (zh) | 航天器返回预定落点的大椭圆轨道变轨规划方法 | |
CN106114910B (zh) | 一种航天器飞行轨道滚动时域控制方法 | |
CN109539903A (zh) | 一种固体运载火箭椭圆转移轨道迭代制导控制方法 | |
CN109269504B (zh) | 一种具有末端约束的姿态机动路径规划方法 | |
RU2424954C1 (ru) | Способ управления движением разгонного блока на участке доразгона | |
CN112607065B (zh) | 一种基于电推进系统的高精度相位控制方法 | |
CN113602535B (zh) | 一种微纳卫星在轨自主交会控制的方法及计算机设备 | |
CN112061424B (zh) | 一种基于融合目标姿态的机动过程能源角动态跟踪方法 | |
CN102878997A (zh) | 一种大偏心率轨道的星上快速高精度外推方法 | |
CN107506505B (zh) | 高精度地月自由返回轨道设计方法 | |
CN113311854A (zh) | 一种月球采样返回任务中的定点着陆轨道设计方法 | |
CN112486196A (zh) | 一种满足严格时间位置约束的飞行器快速轨迹优化方法 | |
CN109992823B (zh) | 一种月面高精度定时定点着陆轨道控制方法 | |
Kos et al. | Altair descent and ascent reference trajectory design and initial dispersion analyses | |
CN108082538B (zh) | 一种考虑始末约束的多体系统低能量捕获轨道方法 | |
CN113569391A (zh) | 地月转移轨道的参数确定方法、装置、设备及介质 | |
CN113310496A (zh) | 一种确定月地转移轨道的方法及装置 | |
CN115343960B (zh) | 一种地月系统中航天器光照阴影规避控制方法 | |
CN113741197B (zh) | 高轨目标的单次抵近三脉冲控制方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210827 |
|
RJ01 | Rejection of invention patent application after publication |