CN119047147B 实现火星探测发射窗口任意时刻起飞诸元设计方法及系统 （深空探测实验室(天都实验室)）.docxVIP

(19)国家知识产权局

(12)发明专利

(10)授权公告号CN119047147B(45)授权公告日2025.07.08

(21)申请号202411033507.7

(22)申请日2024.07.30

(65)同一申请的已公布的文献号申请公布号CN119047147A

(43)申请公布日2024.11.29

(73)专利权人深空探测实验室(天都实验室)地址230088安徽省合肥市高新区创新产

业园三期A1栋10-13楼

(72)发明人刘继忠耿光有王建明罗泰茗尤志鹏汪中生康焱高朝阳李洋史露莎马骁孙宇

(74)专利代理机构北京高文律师事务所11359专利代理师马晓田

(51)Int.CI.

GO6F30/20(2020.01)

B64G1/24(2006.01)

GO6F30/15(2020.01)

G06F119/14(2020.01)

(56)对比文件

CN114818243A,2022.07.29审查员蒋悦

权利要求书4页说明书9页附图5页

(54)发明名称

实现火星探测发射窗口任意时刻起飞诸元设计方法及系统

(57)摘要

CN119047147B本发明公开了一种实现火星探测发射窗口任意时刻起飞的诸元设计方法，涉及火星探测发射轨道设计技术领域。本发明包括：根据发射时间计算探测器出发参数所需初值；根据探测器出发参数初值，设计运载火箭入轨参数，并结合运载火箭入轨参数中的计算量和设计量，通过迭代完成发射轨道设计；调整目标轨道倾角，计算生成随起飞时刻变化规律性强的发射轨道目标诸元。本发明能够使运载火箭实现在当天发射窗口内任意时刻发射起飞，满足探测器奔火星轨道逃逸飞行的精确出发参数要求，节省探测器分离后因起飞时刻偏差引起探测器的额外推进剂消耗，

CN119047147B

CN119047147B权利要求书1/4页

2

1.一种实现火星探测发射窗口任意时刻起飞的诸元设计方法，其特征在于，包括以下步骤：

根据发射时间计算探测器出发参数所需初值；

根据探测器出发参数初值，设计运载火箭入轨参数，并结合运载火箭入轨参数中的计算量和设计量，通过迭代完成发射轨道设计，将对应的起飞时刻作为基准起飞时刻；

调整目标轨道倾角，计算生成随起飞时刻变化规律性强的发射轨道目标诸元；

对生成的多条发射轨道目标诸元，采用最小二乘原理将目标诸元拟合为起飞时刻相对于基准起飞时刻的多项式，获取拟合系数；

将得到的多项式拟合系数装订运载火箭中；

按照实时发射窗口和各根数对应的多项式系数，由多项式函数计算得出实时零窗口发射诸元，实现实时零窗口发射；

其中，运载火箭计算量包括升交点经度、真近角，设计量包括近地点高度、逃逸C3、轨道倾角、近地点幅角；

将窗口内不同起飞时刻轨道倾角、近地点幅角和滑行时间拟合为起飞时刻相对于基准起飞时刻之差的多项式，获取拟合系数，计算方法如下：

c=(ATA)-1Ab(9)

式中c表示拟合系数，b表示样本点结果序列，A是由发射时刻组成的矩阵，对于n个样本点，n≥4,以通过三次曲线拟合轨道倾角为例，表达式如下：

iEAt=c?△t3+c?△t2+c?△t+c?(10)

式中，△t为在发射窗口内发射时刻相对于基准发射时刻的时间差，iE△t为与基准发射时刻相差△t时发射对应的轨道倾角，ck,k=0,1,2,3为拟合系数，已知不少于4组样本点时，即可通过式(9)求解出拟合系数，同理可计算近地点幅角、逃逸C3和滑行时间的拟合系数。

2.根据权利要求1所述的一种实现火星探测发射窗口任意时刻起飞的诸元设计方法，其特征在于：所述根据发射时间计算探测器出发参数所需初值，具体如下：

从探测器分离开始，至到达近火点结束，整个探测火星轨道在日心坐标系下完成积分，飞行动力学方程如下：

式中μs为太阳引力系数，μ;太阳系下扰动天体引力系数，对于火星探测转移轨道，需要考虑8大行星及月球的影响；F和F两项分别指地球和火星的非球形摄动力项，8⑧和δ是开关函数，当探测器在地球和火星引力影响球界内取1,其余取0,8V(t?,t?)表示t?和t?时刻施加和结束深空机动，FsRp表示太阳辐射压项，Fothers表示其它次要摄动影响项，r表示射惯性系下地心至运载火箭质心的矢径；

结合地球出发和火星到达要求