CN111460615B - 一种地月l2点任务发射窗口设计方法 - Google Patents
一种地月l2点任务发射窗口设计方法 Download PDFInfo
- Publication number
- CN111460615B CN111460615B CN202010141056.4A CN202010141056A CN111460615B CN 111460615 B CN111460615 B CN 111460615B CN 202010141056 A CN202010141056 A CN 202010141056A CN 111460615 B CN111460615 B CN 111460615B
- Authority
- CN
- China
- Prior art keywords
- moon
- orbit
- window
- earth
- task
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C21/00—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
- G01C21/24—Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 specially adapted for cosmonautical navigation
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Aviation & Aerospace Engineering (AREA)
- Astronomy & Astrophysics (AREA)
- Automation & Control Theory (AREA)
- General Physics & Mathematics (AREA)
- Navigation (AREA)
Abstract
本发明公开了一种地月L2点任务发射窗口设计方法。本发明根据工程约束、任务需求和轨道设计情况,获取窗口设计约束条件、轨道和任务重要参数等;然后同时对转移段和使命段轨道进行约束条件相关的特性分析,分别筛选出满足任务要求的发射窗口日期集合;最后求解两者并集,生成相应日期下完整的标称飞行轨道进行确认,得到最终的发射窗口设计结果。该方法可满足工程任务对发射窗口的多约束需求。
Description
技术领域
本发明涉及深空探测轨道设计技术领域,具体涉及地月系统L2平动点Halo轨道任务的发射窗口设计方法。
背景技术
平动点是第三体在受两个大天体的万有引力作用时,引力合力与第三体自身惯性力的动平衡点。地月系统L2平动点位于地月连线的延长线上,随月球以相同公转角速率围绕地球旋转,在该点附近运行的航天器能够始终运行在月球背面上空。由于这一特殊的空间几何特性,地月L2点在月球背面着陆探测任务、空间科学观测等任务中具有重要的应用价值。目前,仅有美国和中国开展过地月L2点探测任务,其中我国的嫦娥四号任务通过发射鹊桥号中继星至地月L2点的Halo轨道,建立地球与月球背面探测器的中继通信链路,有效支持了月背着陆及巡视探测任务。
月球探测任务的发射窗口设计一般需要考虑诸如运载发射要求、光照和阴影条件等工程需求约束。对于环月轨道的探测任务而言,环月轨道的光照和阴影为周期性交替,除月食外阴影一般是受月球遮挡所致,光照阴影时长可根据轨道高度、β角等参数利用解析公式快速求解。而地月L2点附近的轨道由于特殊的空间几何位置和轨道构型,光照和阴影的规律与环月轨道差异很大,其阴影是受到地球或月球遮挡而造成的,与初始轨道相位有一定相关性,但既无可用的解析计算公式,也不具有显著的周期性,且持续时间也远长于环月轨道,甚至可达数十小时,严重威胁航天器的在轨安全性。而在鹊桥号中继星之前已开展的地月L2点探测任务中,由于均为拓展试验任务,不涉及运载发射约束,同时由于在地月L2点运行时间不长,也未开展长期运行的阴影规律研究,因此地月平动点任务发射窗口设计方法并无先例可参照。
发明内容
有鉴于此,本发明提供了一种地月L2点任务发射窗口设计方法,通过对地月转移轨道和Halo轨道分别进行特性分析和参数统计,筛选满足工程任务需求的发射窗口。
本发明的地月L2点任务发射窗口设计方法,包括如下步骤:
步骤1,根据工程约束和任务需求,获取发射窗口分析范围、运载火箭对地月转移轨道初始状态的要求、航天器平台对最长阴影时长的约束;
步骤2,根据轨道设计结果,获取地月转移轨道的终端参数、Halo轨道振幅参数、运行时长;
步骤3,利用地月转移轨道的初始和终端参数、转移时间,在发射窗口分析范围内构造一组不同出发时间的地月转移轨道,获得其出发时刻的轨道根数;
步骤4,结合运载火箭对近地点幅角的要求,筛选出满足要求的地月转移轨道,获得相应的发射日期集合;
步骤5,根据发射窗口分析范围和转移时间,计算得到初步的到达地月L2点轨道的时间范围;
步骤6,根据Halo轨道的振幅参数、运行时长和到达时间范围,生成一组不同初始相位的Halo轨道星历;
步骤7,分析不同初始相位的Halo轨道在任务期间的光照阴影条件相关参数,结合航天器平台对最长阴影时长的约束,筛选满足要求的Halo轨道到达时间,并根据转移时间反推出相应的发射日期集合;
步骤8,根据步骤4和步骤7的结果,选择两者并集内的发射日期,生成相应发射窗口的标称轨道星历,确认相关参数均可满足任务要求,最终获取满足工程约束的发射窗口设计结果。
有益效果:
本发明将Halo轨道特性分析与发射窗口筛选相结合,通过对转移轨道和使命轨道解耦并分别开展与约束条件相关的特性分析并需求规律,进而实现多约束条件下的窗口筛选工作。
附图说明
图1为本发明的算法流程图。
图2为基于运载火箭约束的发射窗口设计结果。
图3为基于光照条件约束的发射窗口设计结果。
图4为轨道最长有效阴影统计。
图5为轨道阴影次数统计。
具体实施方式
下面结合附图并举实施例,对本发明进行详细描述。
本发明提供了一种地月L2点任务发射窗口设计方法,用于地月L2点Halo轨道任务的发射窗口分析和设计工作。
为了尽量减少探测器在地月L2点轨道所遭遇的阴影时长和次数,在地月L2点任务的发射窗口设计中不仅需要开展转移轨道对运载约束的适应性分析,还必须开展Halo轨道分析,获取其光照阴影的规律特性,最终筛选出满足任务需求的发射窗口。
本发明通过对转移段和使命段轨道解耦并分别进行特性分析,在满足运载发射约束的同时,尽量减少航天器在Halo轨道上长期运行遭遇的阴影次数和时长。
首先,根据工程约束、任务需求和轨道设计情况,获取窗口设计约束条件、轨道和任务重要参数等;然后同时对转移段和使命段轨道进行约束条件相关的特性分析,分别筛选出满足任务要求的发射窗口日期集合;最后求解两者并集,生成相应日期下完整的标称飞行轨道进行确认,得到最终的发射窗口设计结果。该方法可满足工程任务对发射窗口的多约束需求。
(1)根据工程约束和任务需求,获取发射窗口分析范围[t0,tf]、运载火箭对近地点幅角的要求[ω0,ωf]、地月转移轨道的初始参数,包括近地点高度h0、出发轨道倾角i0和出发轨道类型、航天器平台对最长阴影时长的约束TEclipse。
(2)根据轨道设计结果,获取地月转移轨道的终端参数,包括近月点高度hf、倾角if、到达轨道类型;获取转移轨道飞行时长Δt;获取Halo轨道振幅参数Az和运行时长。
(3)根据地月转移轨道的初始参数、终端参数和转移时间,建立如下两点边值问题:
利用二者线性化误差关系:
给定合适的初值,即可以通过数次迭代就求出转移轨道的数值解。
按上述求解方法,在发射窗口分析范围[t0,tf]内,以天为步长,生成一组满足初始和终端条件的地月转移轨道参数,获取其初始轨道根数。
(4)结合运载火箭对近地点幅角的要求[ω0,ωf],在步骤(3)得到的地月转移轨道中筛选出满足近地点幅角要求的轨道,获得相应的发射日期集合。设计结果如图2所示。
(5)根据发射窗口分析范围[t0,tf]和转移时间Δt,计算得到初步的到达地月L2点轨道的时间范围[t0+Δt,tf+Δt],到达时刻范围即是Halo轨道的初始时刻范围。
(6)根据Halo轨道的振幅参数Az,利用Richardson三阶解析公式,生成南向Halo轨道远月端在地月L2点会合坐标系下的位置速度以此作为轨道初值进行数值积分。由于地月L2点的不稳定性,长时间数值积分后轨道将会逐渐发散,因此为了获取满足运行时长要求的Halo轨道,需要在轨道每圈穿越地月L2点会合坐标系XZ面时施加修正速度脉冲,通过微分修正使其控后第三次穿越XZ面时X向速度为零。按上述步骤依次积分,直至轨道飞行时间满足运行时长要求,即可获得一条Halo轨道星历。
(7)对步骤(6)中获得的Halo轨道组进行遍历,统计Halo轨道的阴影时长、次数、最低光照强度以及间隔时长等参数随初始相位的变化规律;结合航天器平台对最长阴影时长的约束TEclipse,筛选出最长阴影小于TEclipse的Halo轨道,并根据Halo轨道的到达时间和转移时间反推出相应的发射日期,如图3所示;Halo轨道在不同初始相位下的最长有效阴影统计示例见图4和图5。
(8)根据步骤(4)和步骤(7)获取的发射日期结果,选择两者并集内的日期,生成相应的地月转移至Halo轨道的标称轨道星历,并对照运载发射约束和航天器光照条件约束筛选,最终获取满足工程约束的发射窗口设计结果。
需要说明的是,本发明说明书中未作详细描述的内容属本领域专业技术人员的公知技术。
综上所述,以上仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种地月L2点任务发射窗口设计方法,其特征在于,包括如下步骤:
步骤1,根据工程约束和任务需求,获取发射窗口分析范围、运载火箭对地月转移轨道初始状态的要求、航天器平台对最长阴影时长的约束;
步骤2,根据轨道设计结果,获取地月转移轨道的终端参数、Halo轨道振幅参数和运行时长;
步骤3,利用地月转移轨道的初始和终端参数、转移时间,在发射窗口分析范围内构造一组不同出发时间的地月转移轨道,获得其出发时刻的轨道根数;
步骤4,结合运载火箭对近地点幅角的要求,筛选出满足要求的地月转移轨道,获得相应的发射日期集合;
步骤5,根据发射窗口分析范围和转移时间,计算得到初步的到达地月L2点轨道的时间范围;
步骤6,根据Halo轨道的振幅参数、运行时长和到达时间范围,生成一组不同初始相位的Halo轨道星历;
步骤7,分析不同初始相位的Halo轨道在任务期间的光照阴影条件相关参数,结合航天器平台对最长阴影时长的约束,筛选满足要求的Halo轨道到达时间,并根据转移时间反推出相应的发射日期集合;
步骤8,根据步骤4和步骤7的结果,选择两者并集内的发射日期,生成相应发射窗口的标称轨道星历,确认相关参数均可满足任务要求,最终获取满足工程约束的发射窗口设计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010141056.4A CN111460615B (zh) | 2020-03-03 | 2020-03-03 | 一种地月l2点任务发射窗口设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010141056.4A CN111460615B (zh) | 2020-03-03 | 2020-03-03 | 一种地月l2点任务发射窗口设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111460615A CN111460615A (zh) | 2020-07-28 |
CN111460615B true CN111460615B (zh) | 2020-12-11 |
Family
ID=71680089
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010141056.4A Active CN111460615B (zh) | 2020-03-03 | 2020-03-03 | 一种地月l2点任务发射窗口设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111460615B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112525001B (zh) * | 2021-02-10 | 2021-04-16 | 星河动力(北京)空间科技有限公司 | 运载火箭发射窗口的确定方法、装置、设备及存储介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106672266A (zh) * | 2017-01-04 | 2017-05-17 | 北京理工大学 | 考虑时间约束的平衡点Halo轨道调相轨道转移方法 |
CN106682274A (zh) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
CN108460496A (zh) * | 2018-03-28 | 2018-08-28 | 上海航天控制技术研究所 | 一种深空探测全局发射窗口搜索方法 |
CN110077626A (zh) * | 2019-03-14 | 2019-08-02 | 北京空间飞行器总体设计部 | 一种三脉冲地月L2点Halo轨道的捕获方法 |
CN110489779A (zh) * | 2019-07-03 | 2019-11-22 | 上海卫星工程研究所 | 一种木星探测借力飞行轨道优化设计方法 |
-
2020
- 2020-03-03 CN CN202010141056.4A patent/CN111460615B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106682274A (zh) * | 2016-12-05 | 2017-05-17 | 北京理工大学 | 考虑振幅约束的一种Halo轨道在轨保持方法 |
CN106672266A (zh) * | 2017-01-04 | 2017-05-17 | 北京理工大学 | 考虑时间约束的平衡点Halo轨道调相轨道转移方法 |
CN108460496A (zh) * | 2018-03-28 | 2018-08-28 | 上海航天控制技术研究所 | 一种深空探测全局发射窗口搜索方法 |
CN110077626A (zh) * | 2019-03-14 | 2019-08-02 | 北京空间飞行器总体设计部 | 一种三脉冲地月L2点Halo轨道的捕获方法 |
CN110489779A (zh) * | 2019-07-03 | 2019-11-22 | 上海卫星工程研究所 | 一种木星探测借力飞行轨道优化设计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111460615A (zh) | 2020-07-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Goodman | History of space shuttle rendezvous and proximity operations | |
Zimpfer et al. | Autonomous rendezvous, capture and in-space assembly: past, present and future | |
CN110909461B (zh) | 基于可达集概念的地月/月地直接转移轨道设计方法 | |
CN106595673B (zh) | 面对地球静止轨道目标操作的空间多机器人自主导航方法 | |
Bodin et al. | PRISMA: An in-orbit test bed for guidance, navigation, and control experiments | |
Montenbruck | An epoch state filter for use with analytical orbit models of low earth satellites | |
CN111460615B (zh) | 一种地月l2点任务发射窗口设计方法 | |
Kam et al. | LADEE flight dynamics: overview of mission design and operations | |
CN113093246A (zh) | 地面多目标点成像快速判定及任务参数计算方法 | |
Bosanac et al. | A low-thrust-enabled SmallSat heliophysics mission to Sun-Earth L5 | |
Conte et al. | Targeting the Martian moons via direct insertion into Mars’ orbit | |
Sorensen et al. | The Clementine mission—A 10-year perspective | |
Mesarch et al. | Maneuver operations results from the lunar reconnaissance orbiter (LRO) mission | |
Liu et al. | Direction/distance/velocity measurements deeply integrated navigation for venus capture period | |
Guinn et al. | The Deep Space Positioning System (DPS) Navigator Concept for the Lunar Gateway | |
Grebow et al. | Design and execution of Dawn HAMO to LAMO transfer at Ceres | |
Pelletier | Cassini orbit determination performance (July 2008-December 2011) | |
Tolson et al. | Atmospheric modeling using accelerometer data during Mars Atmosphere and Volatile Evolution (MAVEN) flight operations | |
Scala et al. | GNSS-based navigation for a remote sensing three-satellite formation flying | |
Chen et al. | Preliminary mission design and analysis of a lunar far-side positioning CubeSat mission | |
You et al. | Navigating Mars reconnaissance orbiter: Launch through primary science orbit | |
Parker et al. | Navigating a crewed lunar vehicle using liaison | |
CN114894199B (zh) | 一种地月空间航天器的天基测定轨方法 | |
Genova et al. | Transfers from TLI to Lunar Frozen Orbits with Applications to NASA's CLPS & ARTEMIS Programs | |
Ocampo et al. | Multiple-Spacecraft Orbit-Transfer Problem: The No-Booster Case |
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 |