首页 / 专利分类库 / 飞行器;航空;宇宙航行 / 宇宙航行;及其所用的飞行器或设备 / (转入B64G 99/00) / 应用借力机制选择星际探测目标的探测器发射方法

应用借机制选择星际探测目标的探测器发射方法

申请号 CN200610010008.1 申请日 2006-04-30 公开(公告)号 CN100393585C 公开(公告)日 2008-06-11
申请人 哈尔滨工业大学; 发明人 崔平远; 乔栋; 崔祜涛; 栾恩杰;
摘要 应用借 力 机制选择星际探测目标的探测器发射方法,属于深空探测转移轨道技术领域。为了解决轨道半长轴较大或者偏心率较大的探测目标的可接近性评价这一难题,本 发明 的探测器发射方法包括:最优两脉冲转移轨道求取、借力飞行轨道类型选择和借力飞行轨道拼接三个部分。先求取从目标星到发射星体的最优两脉冲转移轨道,把发射星体作为借力星体,两脉冲转移在发射星体处的参数作为借力飞行时的匹配参数,然后采用周期约为发射星体公转周期整数倍的日心大椭圆轨道和远日点处的深空机动,搜索满足匹配条件的发射参数。本发明通过引入借力机制实现对原有两脉冲转移的扩展,提高半长轴较大或偏心率较大目标的可接近性。
权利要求

1、应用借机制选择星际探测目标的探测器发射方法,其特征在于所述方法为: 一、采用异面非共轴椭圆轨道的最优两脉冲计算方法,通过公式(1)~(4)和迭代方法绘制给定真近点的最优总的速度增量等高线图,确定全局最优解两脉冲转移可能出现的区域,采用SQP方法得到从目标星体到发射星体的全局最优的两脉冲转移轨道;

说明书全文

应用借机制选择星际探测目标的探测器发射方法

技术领域

发明属于深空探测转移轨道技术领域,涉及一种应用借力机制选择星际 探测目标的探测器,方法。 背景脉
二十世纪九十年代初,美国、欧空局及日本等国家与机构的空间研究单位, 再次掀起了旨在加深对太阳系的了解,进一步探索深空奥秘,从而揭示宇宙的 演变规律的深空探测活动的高潮。此时,探测目标的重点开始转向太阳系内为 数众多的小行星彗星,即那些体积较小、质量较轻、形状不规则、旋,性 不固定、引力场较弱且分布不均匀的非合作型小天体目标,以及八大行星及其
对于星p示探测任务而言,科学探测目标的选择无疑是任务设计与规划的第 一步,也是重要的一步。在选择探测目标时,首先要考虑它的可接近性。因此 开展星际探测目标的可接近性评价方法研究是非常重要和必要的。目标可接近 性的优劣通常是通过实现与目标星的交会任务所需的总的速度增量的大小来评 价的。
在己有的星际探测目标可接近性评价方法中,在先方法[l](参见Hulkower, N.D" Lau, CO" and Bender, D.F., Optimum Two-Impulse Transfer for Preliminary In鄉lanetaiy Trajectory Design, Journal of Guidance, Control and Pynamic, 1984, 7:458461; Helin, E.F" Hulko窗,N.D., Bender, D.F" The Discovery of脂DB, the Most Accessible Asteroid Known, Icarus, 1984, 57:42~47; C.O.Lau, N.D. Hulkower, "Accessibility of Near-Earth Asteroids", Journal of Guidance, Control and Dynamic, May-June 1987,10(3):225-232)基于异面非共轴椭圆轨道的最优两脉巧
计算方法,给出了通过绘制给定真近点的最优总的速度增量等高线图与圆ii 曲线拼接法相结合来确定目标可接近性的方法。这种方法不但可以求得全局最 优的两脉冲转移轨道的解,估算出从地球直接转移至小天体所需的最小两脉冲 转移的速度增量,而且这种方法是'时间自由"的。该方法是进行目标选择与可
^a性分析的经典方法。
在先方法[2](参见Ettore Perozzi, Alessandro Rossi, Giovanni B. Valsecchi, Basic targeting strategies for rendezvous and flyby missions to the near-Earth asteroids, Planetary and Space Science, 2001,49:3-22; Apostolos A.Christou^ The statistics of flight opporftinities to accessible near-Earth asteroids, Planetary and Space Science, 2003, 51:221-231; Richard P. Binzel, Ettore Perozzi, Andrew S.Rivkin, Alessandro Rossi, Alan W.harris, Schelte J. Bus, Giovanni B.Valsecchi, Stephen M.Slivan, Dynamical and compositioml assessment of near-Earth object mission targets, Meteoritics & Planetmy Science 39, No.3, 351-366,2004)通过采用经典的霍曼转移 方法计算交会任务所需的总的速度增量来评价探测目标的可接近性。该方法计 算简便,适合于大规模的计算与搜索。
随着技术的进步和深空探测任务的进一步开展,对探测目标的选择M趋 向于要求探测目标兼备科学价值和技术可实现性。就技术上的可实现性而言,
主要考虑它的可接近性。评价探测目标可接近性的在先方法主要有两种: 一种 是通过采用全局最优两脉冲转移方式计算交会任务所需的总的速度增量,另一 种是通过采用经典的霍曼转移策略计算交会任务所需的总的速度增量。采用这 些经典的两脉冲转移策略,评价的结果是那些轨道与地球相近的小天体的可接 近性很好(实现交会任务所需的总的速度增量较小),而那些具有较大科学价值, 但轨道半长轴较大或者偏心率较大的小天体的可接近性很差(实现交会任务所 需的总的速度增量和发射能量均很大,已经超出了目前人类技术所能实现的程 度)。这就使得任务的设计者很有可能把科学家们认为极具科学价值的目标排除 在可选的探测目标之外。 发明内容
本发明要解决的技术问题是克服上述已有方法的技术困难,提供一种设计 方式简单,技术上可行的应用借力机制选择星际探测目标的探测器,方法,
从而可以解决轨道半长轴较大或者偏心率较大的探测目标的可接近性评价i^; 难题。
本发明的应用借力机制选择星际探测目标的探测器MT方纟跑括:最优两 脉冲转移mit求取、借力飞行轨道类型选择和借力飞行轨道皿三个部分。先
求取从目标星到,星体的最优两脉冲转移车iiji, teMf星体作为借力星体, 两脉冲转移在皿星体处的参数作为借力飞行时的匹配参数,然后采用周期为 ,星体公转周期»倍的日心大椭圆«和远日点处的深空机动,^i足 匹配劍的,参数,即反向递推设计策略。具,术:^法为:
i、应用借力机帝跑«际探测目标的探测器皿方法,,征在于ff^方
法为:
一、采用异面非共轴椭圆轨道的最优两脉冲计算方法,iiii公式(i)〜(4) 和牛顿迭代方法織U给定真近点角的最优总的ffi增歸高线图,确定全局最 體两脉冲转移可能出现的区域,采用sqp方法得到从目标星体到MT星体的 全局最优的两脉冲转移,
formula see original document page 6 (i)
formula see original document page 6 (2)
formula see original document page 6 (3)
formula see original document page 6 (4)
式中:^和A分别标鄉星体和目标星体^a的单位方向矢量,^为引力常 数,p,为转移车爐的半交弦,《和《分别为初始和糊的健矢量,e为角度改
二、选择两脉冲转移的皿星体作为借力星体,探测器imf星体^M进
入飞行周期为,星体公转周期MW咅,近日点,为,星体平均公转,
的日心大椭圆勒道,根据目标星体的参数,M公式(5)〜(7)确定日心大
椭圆转移IWt的周期P,、远日点速度v。、远日点,r。:
formula see original document page 6 (6)
formula see original document page 6(7)
式中:^是太阳的引力常数,〜为近日点判5;
三、以两脉冲转移在微星体处的参数作为借力飞行时的匹配参数,采用 周期为发射星体公转周期整数倍的日心大椭圆轨道和调整远日点处的深空机 动,在远日点雖加一个深空机动吣,使得近日点^g:小于織星体的鹏半 径,探测器的鹏和鄉星体的轨道相切,M采用C,即双曲^i的平方的 方法将不同的執道段,起来,利用公式(8)〜(12)确定探领!t^深空机动后
轨道的长半轴。r,轨道周期g,偏心率e。探测器返回iW的远日点速度v、从 ,到与皿星体交会时的飞行时间7;:
formula see original document page 7
式中:、为探测器从交会点至顿始l^日点的时间。
本发明的基本原理:,经典的两脉冲转移方法计算交会任务所需的« 增量包括两锦分:魏时的驗增量和交会时的驗增量。对于半长轴较大 或者偏心報大的目标,通常会出现交会时相对速度增量较小,而鄉时所需 的速度增量很大,从而导致总的舰增量和織能翻艮大。正是基于这一点, 劍考虑将借力飞行机制弓l入到星际探测目标的可舰性i愤方法中。借力飞 行可以有效降低星际探测任务所需的鎌倉糧,进而M^、总的纖增量。同时, 为了减少动力学的要求和避^f时间的^W,这里OT皿星体借力和深空机 动相结合的策略设计转移M,取代两脉冲转移方法。
本发明与深空探测转移fiii有关,特别是深空探测目标的选择与可,性
i愤分析。该方法将用在深空探测(特别是小天條测)的目織择、可舰
性iffc等方面,也可为深空探测任务的设计与规划ii^一种im与分析的ifii径。本发明与在先方法[l]、在先方銜2]的不同之处在于,本发明M过引入借 力机制实5舰原有两脉冲转移的扩展,提高半长轴较大或偏心率较大目标的可 接近性。其优点在于-
1) 减少半长轴较大或偏心M大的目标星体实现交会任务所需的速度增量 和发射能量,为兼备科学价值和技术可实现性探测目标的选取,提供重要的参 考;
2) 选择两脉冲转移的发射星体作为借力星体,减小了设计方法的旨性, 同时也避免了将时间约束弓l入刑介方法;
3) 反向递推的设计策略,有效地保存了原有最优两脉冲转移轨道的«信 息,避免了不必要的重复计算;
4) 这种弓i力机制的星际探测目标可接近性刑介方法也是对间自由"的,可 用于交会机会的预测。
附图说明
图1为两脉冲转移的示意图,图2为交会4015小行星交会的最优两脉冲总 的速度增量等高线图,图3为发射星体借力飞行示意图,图4为,星体借力 飞行转移的性能关系图,图5为猜测的似,,。和匹配误差5(:之间的关系图,图6 为采用2年周期地球借力转移方案与4015小行星交会的飞行轨迹图。
具体实施方式一:本实施方式结合图1~6对本发明进行详细介绍:
l)鄉星体到目标星体的最优两脉冲转移fiii该段轨道的示意图,如图i所示。
这里定义:〖《,(发射时的双曲线M), /2《2 (交会时的双曲线,)。 假设探测器与鄉星体分离的速度为^、与目标星体交会时的職为^,则为:
R, =±("z《)=7,(1) ^ =±(》—z《)=(^ —72 (2)
这里^,和6分别是Mt星体和目标星体绕太阳公转的速度矢量,以上两式 中+号是措'短程"(即鹏小于;rM的角度改变量0实现的f爐转移),-号是 措'长程"(即通过大于;r弧度而小于2;r弧度的改变量实现的轨道转移)。《和仏
分别表示发射星体和目标星体位置的单位方向矢量,^和z可由如下方程(3)和formula see original document page 9这里采用牛顿迭代的方法可以求得满足最优解的A 。对于发射星体和目标 星,所在的公转轨道的每一对位置,都可以求得满足最优解的p,,而一个p,可
以计算得到两个J (探测驗日心转移轨道所需的总的驗增量)即"长程"情况 下的J—和"短程"情况下的L ,这里仅保留其中较小的一个。给定不同鄉星
体的平近点角和交会星体的平近点角,通过(5)〜(8)式和牛顿迭代方法就可 以得到最优两脉冲转移等高线图。这里以4015小行星为目标星体,地球为皿 星体,给出其最优两脉冲转移flit的等高线图,如图2所示。
由图2可以得到两脉冲转移所需总的速度增量最小的区域,从而猜测出较 好的初值,通过采用SQP (序列二次规划)方法求解出精确的^tt两脉冲转移 轨道,具体的参数如表1所示。全局最优两脉冲转移方法求解出来的地球在发
射时的平近点角可作为借力方案中地球在借力时的平近点角。 2)借力飞行轨道类型选择
借力飞行可以有效的斷氐星际探测任务所需的发射能量和总的速度增量, 这里采用借力飞行技术去扩展经典的两脉冲转移策略。为了减少动力学的要求
和避免时间的,,采用mt星体借力的策略。这种借力飞行策略可以描述为: 一个探测器从皿星体,进入飞行周期比发射星体公转周期MW咅稍大,近 日点半径为,星体平均公转半径的日心大椭圆轨道。在远曰点处施加一, 度脉冲,使得近日点判圣小于发射星体的轨道半径,探测器的轨道和发射星体 的车爐相切,这样就可以利用鄉天体的甩摆作用,减小探澳!器的鄉能量和 总的速度增量。该借力飞行策略如图3所示。
假设发射星体公转轨道为圆轨道,从发射星体发射的探测器速度、比标称
轨i!3I度稍大,方向平行于发射星体的速度K,.,这种情况下,飞行器在日心轨道
的近日点开,动,则近日点的速度^,为:
v,K,.十K, (9) 探测器从发射星体发射时相对于日心的位置为。,,则有探测器日心位置等 于魏星体的日心位置,艮P-
(10)
己知。,和v,,,即可以确定探测器远日点半径r。、速度v。和探测l^刀始轨道 周期S为:
这里是太阳的引力常数。探测器初始轨道周期比标称的轨道周期稍 长,如果将此轨道递推到近日点,发射星体必定不在原来的位置。若假设在远 日点处加一个深空机动气,加深空机动的目的是修正远日点处的速度v。,使 其能与地球交会,因此返回地球的远日点速度v为:
-化, (】2) 探测器在远日点处的半径不变,但在近日点处的半径减少,这样就可以与
織星体公转轨道相交,实现借力飞行。远日点速度^与远日点半径r。已知, 就可以确定从远日点返回与发射星体交会段轨道的一些其它的性质。探测器深 空机动后轨道的长半轴。,,轨道周期e,偏心率e.即为:
table see original document page 10
(13)
7>(», (14)
探测行器从,到与,星体交会时的飞行时间为-
这里、为探测器从交会点到初始轨,日点的时间,该时间可以通过幵普 勒方程求解。为了选择借力轨道的类型,这里以地球作为发射星体,给出该转 移轨道的性能。借力飞行后的远日点半经与总的速度增量的关系如图4所示。
由图4可以看出,对于借力S行后远日点半经在1.6至5.5AU (1AU =
1.4959787xWkm)之间的借力类型中,2年周期的地球(皿天体)借力所需的
总的速度增量是最低的。这是一个重要的性质,因为小行星和彗星多数都在这 个距离范围内。因此,选择2年周期的地球借力转移轨道,该转移方案包括一 次深空机动和一次地球借力,从地球发射和到地球借力探测器的飞行时间大约 为2年。
3)借力飞行轨道,
探测器从发射天体,进入"个以发射天体公转周期整数倍为周期,近日 点半径为发射天体平均公转轨道半径的日心大椭圆轨道。在远日点处,施加一 ^NI度脉冲,使f雜测器鹏与鄉天体轨道相交。若探测驗魏时对应发 射天体的平近点角为化,,mt天体在借力飞行时对应的平近点角为A,则从深 空机动点到借力飞行点的飞行时间/,为:
(T/2;r) ;r < Af9 - A/,.0 < 2;r
(/W、. 一 W,.。)7T
i80
(仏—M》
+ 2;r
(15)
180
这里r为发射天体公转的轨道周期。通过求解从深空机动点到借力— 兰伯特问题,可以得到飞入借力天体时的双曲线超速^ 。飞出借力天体时的双
曲^3I^可以从两脉冲转移轨道的参数中得到。通过采用c,匹配(c,即双曲 线超速的平方〉方法将不同的轨道段拼接起来。c,匹配,即探测器飞入借力天
体的双曲线超速和飞出借力天体的双曲线超速的大小相等。假设匹配误差c^
为:
W=|W-W (.16) 若假设iW为皿天体,4015小行星为目标星体,采用2年周期的地球借力,
则地球在g时的平近点角和匹配误差sc',的关系如图5所示。
如图5所示,有两次借力,因为皿曲线经过两次零点。在零点附近取初 值,然后采用牛顿法等可以精确的求得匹配的参数。2年周期的it^H昔力转移方
案的总的速度增量包括:发射时的速度增量W,、远日点处的深空机动气和交
会时的iIJt脉冲AF。。与4015小行星交会的最优两脉冲转移和2年周期的iW 借力转移方案的参数如表1所示。2年周期的地球借力(2)飞行方案的飞行轨 迹如图6所示。
表l在先方法与本发明方法的比较
汲)'
修州凍 平贴龙M'(度)
AK 一)
7}沃)
最优两脉冲(在先方法) 257.175 3U.274 6.799 67.U9 1 328.9
2年周期的itel借力(l〉(本发明方法) 210.871 311274 6.409 26.497 2106.4
2年周期的ife^借力(2)(本发B月方法) 301289 311274 5.632 25.630 2014.7 从表l可以看出,与在先方法相比,采用本发明方法最少可以使总的速度
增量AK和皿能量c,分别减少1.168 km/s和41.489 km2/s2。
从借力转移轨道的设计过程可以看出:探测器的l行轨迹可M发射时的 平近点角A^,借力飞行时的平近点角仏,交会时的平近点角M。来描述,并且
这种转移方案对时间是自由的或者对星历是自由的。这种方法比较适合评价星 际目标星体的可接近性。本发明方法的总的速度增量由发射时的速度增量,远 日点处的深空机动和交会时的速度增量组成。这种方法能够m^皿能量和总 的速度增量是以增加一个深空机动和飞行时间为代价的。这种方法适合于半长 轴较大或者偏心率较大的小天体。本发明方法的评价结果与在先方法[l]和在先
方法[2]评价的结果比较,如表2所示。
表2本发明方法与在先方法对^im大半长轴a偏心率目标zM亍星可,性的比较
本次工作
在先方法 在先方法
小行星 偏心率 本发明方'汰 在先方法 [2] 卩〗
(AU) △ K,。,。/ c、 Af/',,,。/ c3 A K'of 。/
(km/s) (km2" (knVs) (km/s) (km/s)
(7341)1991 VK 2.776 (X506 4.968 25.34 5.776 56.80 7.6 —
(4179)Toutatis 4.122 0.635 5.092 25.15 6.159 63.45 83 —
(3288)Selewus 2.962 0.457 5.180 26.48 5.887 51,24 8.0 5.912
(3908) Nyx 2.812 0.459 5.193 25.41 5.595 38.62 6.9
(8034) 1992 LR 2.579 0.409 5339 26.49 5.578 34.73 6.6
(1627) lvar 2.603 0397 5.389 2636 6.!08 50.59 82 6.117
(3551)Verenia 3.113 0.488 5.398 25.65 6.476 63.92 8.9 6.483
(6489)Golevka 4.009 0.605 5馬 26.48 6.499 63.94 83 —
(433) Eras 1.783 0223 5.547 25.40 5.935 38.38 7.6 5.947
(3352)McAuliffe 2.572 0.369 5.552 25.26 5.891 36.48 7.4 5,900
(4015) WilsonHarrington 4.285 0.623 5.632 25.63 6.799 67.12 8.6 一
(887)Alinda 3.885 0.563 5.635 25.38 6.812 69.24 9.4 6.852
(13651)丽BR 1.744 0.306 5.668 26.20 8. !55 128.7 9.5 一
(31345) 1998 PG 雇 0392 5.948 2537 6.437 42.75 82 —
(3】02)Krok 3.116 0.449 6.039 26.02 6.845 53.54 8.8 6.855
(1685) Toro 1.963 0.436 6.051 25.47 7.675 71.52 7.6 7.675
(1620)Geographos 1.663 0.335 6.148 25.67 88.50 8.2 8.537
(2】00)Ra-Shalom 1.195 0.437 6.353 25.62 7.939 85.52 9.7 7.949
(6178) 1986 DA 4.457 0.587 6.355 26.42 7.622 74.16 9.1 —
(2063) Bacchus 1.455 (X349 6J86 25.30 7.094 49.36 6.8 7.105
(7753) 1卿XB 2.174 0.482 6.421 25.65 6.750 36.61 6.8 一
(B3671) Dionysus 3.389 0.542 6.631 26.48 7.795 68.75 9.8 一
(35396)簡XFll 2.141 0.484 6.701 26.05 7.174 41.07 7.0 一
卿l) 1994AH2 4.330 0.709 6.899 26.3i 8.111 67.69 9.9 —
(22(H) Oljato 3.721 0.713 6.911 26.40 謝0 68.34 9.4 8,037
在表2中,e和e分别表示小行星的远日点距离和偏心率。列在表2中的小
行星都有很高的科学价值。如表2所示,通过采用本发明方法,小行星交会任
务所需的总的速度增量AK和发射能量C,都有明显较少。例如可旨跑源于彗星的 4015小行星(e = 0.623 ,结果如表1所示);据雷达观测,有特殊自旋状态的小 行星4179 (e-0.635),与在先方法相比,采用本发明方樹吏得总的速度增量和
皿能量分别减少17.32%和60.36%;可能起源于4号巨小行星Vesta表面残片 的6489小行星(e = o.6o5),采用本发明方法,使得总的速度增量和,肖讀分 别减少1.0530km/s和37.46kmW等等。从表2中,还可以看出一些目标(例如 4179, 1627, 3551, 6489, 4015, 887, 13651等)都表现出了良好的可接近性。
QQ群二维码
意见反馈