CN111361762B - 一种地月转移轨道发动机试喷方法 - Google Patents
一种地月转移轨道发动机试喷方法 Download PDFInfo
- Publication number
- CN111361762B CN111361762B CN202010141058.3A CN202010141058A CN111361762B CN 111361762 B CN111361762 B CN 111361762B CN 202010141058 A CN202010141058 A CN 202010141058A CN 111361762 B CN111361762 B CN 111361762B
- Authority
- CN
- China
- Prior art keywords
- correction
- midway
- engine
- midway correction
- time
- 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.)
- Active
Links
Images
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
Abstract
本发明公开了一种地月转移轨道发动机试喷方法,将中途修正策略与试喷需求进行联合设计,在确保转移轨道近月点状态的前提下,以推进剂最优的原则满足发动机试喷需求。将发动机试喷与地月转移中途修正结合,根据地月转移轨道某一时刻的位置速度数据,计算中途修正控制量并结合发动机试喷要求确定合适的试喷时机和策略,可以覆盖不同误差情况的推进剂最优试喷方法,可同时满足转移轨道终端目标要求和发动机试喷需求,具有良好的工程可操作性。
Description
技术领域
本发明涉及深空探测轨道设计技术领域,具体涉及一种地月转移轨道发动机试喷时机和策略的确定方法。
背景技术
在月球探测任务中,探测器由运载火箭直接送入地月转移轨道后,需要在转移飞行过程中对大推力轨控发动机进行一次试喷,用于验证大推力发动机的性能,以确保后续近月制动等重要变轨事件的可靠性和安全性。然而,轨控发动机试喷一般都有最短时长要求,在运载发射和定轨精度较高的情况下,这一最短时长所产生的速度增量一般较大,如果在地月转移过程中直接增加一次变轨用于发动机试喷,不但会改变近月点的目标状态,还会增大地月转移轨道控制的难度,导致飞行程序和测控资源的紧张。因此,需要对大推力轨控发动机试喷的时机和策略进行研究。
发明内容
有鉴于此,本发明提供了一种地月转移轨道发动机试喷方法,将中途修正策略与试喷需求进行联合设计,在确保转移轨道近月点状态的前提下,以推进剂最优的原则满足发动机试喷需求。
本发明的地月转移轨道发动机试喷方法,首先根据轨道状态计算第一次中途修正量,若发动机试喷产生速度增量小于或等于第一次中途修正量,则在第一次中途修正时进行发动机试喷;否则,计算在不执行第一次中途修正的基础上的第二次中途修正量,若发动机试喷产生速度增量小于或等于所述第二次中途修正量,则第一次中途修正不执行,在第二次中途修正时进行发动机试喷;若发动机试喷产生速度增量大于所述第二次中途修正量,则计算发动机试喷产生速度增量与所述第二次中途修正量的差值,将所述差值设定为第一次中途修正的反向修正量,执行第一次中途修正,在第二次中途修正时进行发动机试喷。
较优的,具体包括如下步骤:
步骤1,获取探测器的初始轨道参数,包括初始时刻状态、中途修正点时刻、转移轨道终端时刻、发动机试喷最小开机时长对应的速度增量ΔVE;
步骤2,将步骤1中的初始时刻位置速度分别预报至第一次中途修正点和转移轨道终端时刻,求得近月点目标状态与预定量的偏差;
步骤3,根据目标状态偏差计算第一次中途修正量;
步骤4,更新第一次中途修正点的速度,并预报至转移轨道终端时刻,若满足要求,则获得第一次中途修正量大小ΔV1;若不满足,返回步骤3;
步骤5,若ΔV1>ΔVE,则在第一次中途修正时实施发动机试喷;若ΔV1<ΔVE,则将步骤1中的初始时刻位置速度预报至第二次中途修正时刻;
步骤6,根据所得目标状态的偏差计算在不实施第一次中途修正前提下的第二次中途修正量;
步骤7,更新第二次中途修正点的速度,并预报至转移轨道终端时刻,若满足要求,则获得不实施第一次中途修正前提下的第二次中途修正量大小ΔV2;若不满足,返回步骤6;
步骤8,若ΔV2>ΔVE,则不实施第一次中途修正,在第二次中途修正时实施发动机试喷;若ΔV2<ΔVE,进入步骤9;
步骤9,对转移轨道终端时刻进行调整,直至第二次中途修正量ΔV2=ΔVE;若ΔV2始终小于ΔVE,进入步骤10;
步骤10,固定第二次修正量大小为ΔVE,按步骤2~4,调整第一次修正的速度分量瞄准近月点tf的目标状态;选择第二次中途修正的高度角和方位角作为优化变量,直到获得速度增量最优的第一次中途修正量;执行第一次中途修正,并在第二次中途修正时实施发动机试喷。
有益效果:
本发明巧妙地将发动机试喷与地月转移中途修正结合,根据地月转移轨道某一时刻的位置速度数据,计算中途修正控制量并结合发动机试喷要求确定合适的试喷时机和策略,能够覆盖不同误差情况的推进剂最优试喷方法,可同时满足转移轨道终端目标要求和发动机试喷需求,具有良好的工程可操作性。
附图说明
图1为本发明试喷方法流程图。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种地月转移轨道发动机试喷方法,考虑在地月转移中途修正时实施发动机试喷。地月转移中途修正用于修正运载入轨、定轨等误差导致的轨道偏差,以保证探测器能够以预期的近月点状态到达月球。一般地月转移段安排3次中途修正,实际飞行中会根据轨道参数和控制执行情况取消1-2次中途修正。但中途修正的变轨量通常都较小,不需采用大推力发动机实施。特别是在运载发射和定轨精度较高的情况下,轨控发动机试喷最短时长所产生的速度增量往往会超出转移轨道本身所需要的中途修正量。因此,需要综合考虑中途修正安排,合理选择试喷的时机和策略。
本发明首先根据轨道状态计算第一次中途修正量,若发动机试喷产生速度增量小于或等于第一次中途修正量,则在第一次中途修正时进行发动机试喷;否则,计算在不执行第一次中途修正的基础上的第二次中途修正量,若动机试喷产生速度增量小于或等于所述第二次中途修正量,则第一次中途修正不执行,在第二次中途修正时进行发动机试喷;若动机试喷产生速度增量大于所述第二次中途修正量,则计算动机试喷产生速度增量与所述第二次中途修正量的差值,将所述差值设定为第一次中途修正的反向修正量,执行第一次中途修正,在第二次中途修正时进行发动机试喷。
具体的,如图1所示,本发明试喷方法包括如下步骤:
(1)获取探测器的初始轨道参数,包括初始时刻状态、中途修正点时刻、转移轨道终端时刻,发动机试喷最小开机时长对应的速度增量ΔVE;
(2)将步骤(1)中的初始时刻位置速度分别预报至第一次中途修正点和转移轨道终端时刻,求得近月点目标状态与预定量的偏差;
(3)根据转移轨道终端目标状态偏差计算第一次中途修正量;
(4)更新第一次中途修正点的速度,并预报至转移轨道终端时刻,若目标状态满足要求,则获得第一次中途修正量大小ΔV1;若不满足,返回步骤(3);
(5)若ΔV1>ΔVE,则在第一次中途修正时实施发动机试喷,若ΔV1<ΔVE,则将步骤(1)中的初始时刻位置速度预报至第二次中途修正时刻;
(6)根据所得终端目标状态的偏差,在不实施第一次中途修正前提下,计算第二次中途修正量;
(7)更新第二次中途修正点的速度,并预报至转移轨道终端时刻,若满足要求,则获得不实施第一次中途修正前提下的第二次中途修正量大小ΔV2;若不满足,返回步骤(6);
(8)若ΔV2>ΔVE,则不实施第一次中途修正,在第二次中途修正时实施发动机试喷,若ΔV2<ΔVE,进入步骤(9);
(9)对转移轨道终端时刻进行调整,直至第二次中途修正量ΔV2=ΔVE;若ΔV2始终小于ΔVE,进入步骤(10);
(10)采用第一次和第二次中途修正联合控制策略,固定第二次修正量大小为ΔVE,按步骤(2)-(4),调整第一次修正的速度分量瞄准近月点tf的目标状态;选择第二次中途修正的高度角和方位角作为优化变量,直到获得速度增量最优的第一次中途修正量。
本文将发动机试喷需求与地月转移中途修正策略相结合进行设计,以保证在不同误差条件下,都能以中途修正速度增量最优的策略实施发动机试喷。
下面给出一个具体实例进行说明:
(1)获取探测器的初始轨道参数,包括初始时刻t0的位置速度(r0,v0)、第一次和第二次中途修正点时刻t1和t2、转移轨道终端时刻tf,发动机试喷最小开机时长对应的速度增量ΔVE;
(2)将步骤(1)中的初始时刻位置速度(r0,v0)分别预报至t1和转移轨道终端时刻tf。根据终端时刻的位置速度(rf,vf),求解终端目标状态。
工程上要求的地月转移终端目标变量一般为近月点高度Hm,近月点倾角和探测器相对月心的位置矢量和速度矢量垂直(即近月点),将终端状态用q表示:
由q2和q3的表达式可知:q2表示探测器在近月点时相对于月球的速度和月心矢径的夹角的余弦,q3表示探测器在近月点时相对于月球赤道面的瞬时轨道倾角。
求解终端目标状态与预定状态q*的偏差
Δq=q-q*
(3)探测器在第一次中途修正处的速度矢量记为
与终端目标状态的关系可描述为q=f(v1),则可根据终端目标状态偏差计算第一次中途修正速度增量
(4)更新第一次中途修正点的速度为v1=v1+Δv1,并预报至转移轨道终端时刻,若满足要求,则获得第一次中途修正速度增量Δv1;若不满足,返回步骤(3)继续迭代至满足目标要求;
(5)对步骤(4)获得的Δv1求模
ΔV1=|Δv1|
若ΔV1≥ΔVE,则在第一次中途修正时实施发动机试喷,若ΔV1<ΔVE,则将步骤(1)中的初始时刻位置速度预报至t2,获得第二次中途修正处的速度矢量
(6)参照步骤(3)的方法,将第二次中途修正处的速度矢量与终端目标状态的关系记为q=g(v2),则可根据所得终端目标状态的偏差计算在不实施第一次中途修正前提下的第二次中途修正速度增量
(7)更新第二次中途修正点的速度为v2=v2+Δv2,并预报至转移轨道终端时刻,若满足要求,则获得不实施第一次中途修正前提下的第二次中途修正速度增量Δv2;若不满足,返回步骤(6);
(8)对步骤(7)获得的Δv2求模
ΔV2=|Δv2|
若ΔV2≥ΔVE,则取消第一次中途修正,在第二次中途修正时实施发动机试喷,若ΔV2<ΔVE,进入步骤(9);
(9)在[tf-30min,tf+30min]的范围内对转移轨道终端时刻进行遍历,并参照步骤(6)(7)计算第二次中途修正量,直至第二次中途修正量ΔV2=ΔVE;若完成遍历后,ΔV2始终小于ΔVE,进入步骤(10);
(10)采用第一次和第二次中途修正联合控制策略,固定第二次中途修正量的大小为ΔVE,参照步骤(2)~(4),调整第一次修正的速度分量瞄准近月点tf的目标状态,求解Δv1;选择第二次中途修正的高度角α2和方位角β2作为优化变量,高度角α2、方位角β2与Δv2之间的转换关系如下
ΔV2=|Δv2| Δv2x=ΔV2sinα2sinβ2
α2=cos-1(ΔV2z/|Δv2|),Δv2y=ΔV2sinα2cosβ2
β2=tan-1(Δv2y/Δv2x) Δv2z=ΔV2cosα2
将高度角α2在[-90°,90°]范围内,方位角β2在[0°,360°]遍历,求解使得两次中途修正总速度增量最优的高度角α2和方位角β2。
需要说明的是,本发明说明书中未作详细描述的内容属本领域专业技术人员的公知技术。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种地月转移轨道发动机试喷方法,其特征在于,包括如下步骤:
步骤一,根据轨道状态计算第一次中途修正量:
S11,获取探测器的初始轨道参数,包括初始时刻状态、中途修正点时刻、转移轨道终端时刻、发动机试喷最小开机时长对应的速度增量ΔVE;
S12,将S11中的初始时刻位置速度分别预报至第一次中途修正点和转移轨道终端时刻,求得近月点目标状态与预定量的偏差;
S13,根据目标状态偏差计算第一次中途修正量;
S14,更新第一次中途修正点的速度,并预报至转移轨道终端时刻,若满足要求,则获得第一次中途修正量大小ΔV1;若不满足,返回S13;
步骤二,若发动机试喷最小开机时长对应的速度增量ΔVE小于或等于第一次中途修正量ΔV1,则在第一次中途修正时进行发动机试喷;否则,执行步骤三;
步骤三,计算在不执行第一次中途修正的基础上的第二次中途修正量:
S31,将S11中的初始时刻位置速度预报至第二次中途修正时刻;
S32,根据所得目标状态的偏差计算在不实施第一次中途修正前提下的第二次中途修正量;
S33,更新第二次中途修正点的速度,并预报至转移轨道终端时刻,若满足要求,则获得不实施第一次中途修正前提下的第二次中途修正量大小ΔV2;若不满足,返回S32;
步骤四,若发动机试喷最小开机时长对应的速度增量ΔVE小于或等于所述第二次中途修正量ΔV2,则第一次中途修正不执行,在第二次中途修正时进行发动机试喷;若发动机试喷最小开机时长对应的速度增量ΔVE大于所述第二次中途修正量ΔV2,则执行步骤五;
步骤五,对转移轨道终端时刻进行调整,直至第二次中途修正量ΔV2=ΔVE,则第一次中途修正不执行,在第二次中途修正时进行发动机试喷;若ΔV2始终小于ΔVE,进入步骤六;
步骤六,固定第二次修正量大小为ΔVE,按步骤一,调整第一次修正的速度分量瞄准近月点tf的目标状态;选择第二次中途修正的高度角和方位角作为优化变量,直到获得速度增量最优的第一次中途修正量;执行第一次中途修正,并在第二次中途修正时实施发动机试喷。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010141058.3A CN111361762B (zh) | 2020-03-04 | 2020-03-04 | 一种地月转移轨道发动机试喷方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010141058.3A CN111361762B (zh) | 2020-03-04 | 2020-03-04 | 一种地月转移轨道发动机试喷方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111361762A CN111361762A (zh) | 2020-07-03 |
CN111361762B true CN111361762B (zh) | 2020-12-11 |
Family
ID=71200257
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010141058.3A Active CN111361762B (zh) | 2020-03-04 | 2020-03-04 | 一种地月转移轨道发动机试喷方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111361762B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113071712B (zh) * | 2021-03-10 | 2022-07-29 | 北京空间飞行器总体设计部 | 一种月地转移入射变轨策略快速计算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1076005B1 (en) * | 1999-08-13 | 2007-01-03 | Hughes Electronics Corporation | Spacecraft orbit control using orbit position feedback |
CN100491204C (zh) * | 2007-12-26 | 2009-05-27 | 北京控制工程研究所 | 一种利用定轨数据标定加速度计的方法 |
CN102591209B (zh) * | 2012-02-06 | 2014-06-04 | 哈尔滨工业大学 | 卫星闭环测试喷气推力控制的计算方法 |
CN109606739B (zh) * | 2019-01-18 | 2020-10-27 | 哈尔滨工业大学 | 一种探测器地月转移轨道修正方法及装置 |
-
2020
- 2020-03-04 CN CN202010141058.3A patent/CN111361762B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN111361762A (zh) | 2020-07-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rayman et al. | Design of the first interplanetary solar electric propulsion mission | |
CN109398762B (zh) | 一种固体火箭基于椭圆转移轨道的入轨式弹道设计方法 | |
JP6271043B2 (ja) | 軌道制御装置および衛星 | |
US10046869B2 (en) | Inertial sensing augmentation for navigation of spacecraft | |
CN105353621A (zh) | 一种地球静止轨道卫星电推力器故障模式推力分配方法 | |
CN109625334A (zh) | 卫星动量轮在轨姿态无偏差起旋和消旋的控制方法 | |
CN111361762B (zh) | 一种地月转移轨道发动机试喷方法 | |
CN113602532A (zh) | 一种固体运载火箭入轨修正方法 | |
CN111460614B (zh) | 一种地月l2点转移轨道中途修正方法 | |
CN110304279A (zh) | 一种电推进卫星的质心在轨标定补偿方法 | |
CN115015983B (zh) | 光学遥感星座的长期在轨维持方法 | |
CN113525721A (zh) | 卫星轨道变换方法、装置、电子设备及存储介质 | |
Kos et al. | Altair descent and ascent reference trajectory design and initial dispersion analyses | |
CN112649006A (zh) | 一种太阳同步圆轨道的轨道规划方法 | |
CN115562001A (zh) | 一种电推进推力标定方法 | |
CN113343442B (zh) | 一种求解固定时间有限燃料多脉冲转移轨道的方法及系统 | |
Brunner et al. | Comparison of numerical predictor-corrector and Apollo skip entry guidance algorithms | |
US6317661B1 (en) | Argument of perigee correction with longitude control for inclined, eccentric, geosynchronous satellites | |
CN113447043A (zh) | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 | |
CN115574666B (zh) | 一种掠地巡航靶标定高方法 | |
CN111290433B (zh) | 一种长期自主编队联合管道保持方法 | |
CN115610705A (zh) | 一种定制化保持环约束下控制量自动计算方法 | |
Hawkins et al. | Guidance Modifications and Enhancements for Space Launch System Block-1 in Support of Artemis I and Beyond | |
CN113311854A (zh) | 一种月球采样返回任务中的定点着陆轨道设计方法 | |
CN107450580B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |