删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

直接再入大气的月地转移轨道设计*

本站小编 Free考研考试/2021-12-25

自中国探月工程开展以来,在地月转移轨道设计、中途修正等方面取得了丰富的研究成果,相关技术已成功应用于工程实践[1-3]。根据探月工程总体规划,中国将继续开展月球无人采样返回任务,轨道设计及轨道控制面临着新的挑战和难点。月地转移并不是地月转移简单的逆过程,相关约束条件更加复杂。根据月球采样返回任务的需求,从月球停泊轨道出发直接再入大气的月地转移轨道是月球采样返回任务的首选方案。
月地转移轨道设计不同于地月转移轨道。首先通过调整发射窗口,可容易获取共面的地月转移轨道,但由于月球停泊轨道面的固定,月地转移的初始状态不一定满足共面转移的条件,则需要进行轨道面的调整,并且探测器再入大气返回地球时,受回收场的约束,需要匹配落点的地理位置,避免轨道面的调整,保证任务可靠实施。
目前,国内外针对月地转移轨道设计展开了深入的研究。Ocampo和Saudemont[4]研究了利用单脉冲和三脉冲机动方式从环月轨道出发进入月地转移轨道的问题,但未考虑再入点约束。高玉东等[5]基于双二体模型假设,利用二分法搜索从月球停泊轨道到地球停泊轨道的共面月地转移轨道,但未涉及直接再入大气的约束。张磊等[6-7]以近月点速度、月球停泊轨道升交点赤经及近月点幅角为设计参数,以再入点地心距、再入角和再入倾角为终端参数,通过两级修正设计满足约束的月地转移轨道,在此基础上通过分别调整再入倾角和飞行时间完成对落点地理位置匹配,求解过程复杂。郑爱武等[8-10]以探测器出月球影响球的时刻,位置和速度为中间变量,将月地转移轨道分为地心段和月心段,通过调整影响球边界处的位置和速度获得一条连续的月地转移轨道,在此基础上给出了精确月地转移轨道设计方法和月地返回窗口搜索策略,但初值选择依赖于设计者的经验且未涉及月球停泊轨道面调整。汪中生等[11]分析了多种月地转移方案,并给出了相应的基本设计方法,但其求解效率尚有改进和优化的空间。
由于三体问题不存在解析解,月地转移轨道设计必须采用数值方法。一般的思路是:根据边值约束,采用简化模型设计初值,之后考虑精确摄动模型,采用微分修正等数值方法求解精确转移轨道。在目前的研究成果中,伪状态理论[12]是一种计算三体问题的近似解析形式,由于考虑地月球中心引力共同作用,预报精度相对较高,广泛应用于三体轨道的初值设计。
针对以上研究现状,本文基于伪状态理论, 以月球停泊轨道为出发点,考虑再入点的边界约束条件,通过简单迭代设计单脉冲月地转移轨道的初值,在此基础上利用微分修正方法求解精确转移轨道,并对月地返回窗口进行搜索和分析。
1 动力学模型 在地球惯性系下,建立探测器动力学模型,模型中考虑了地球非球形项、月球及太阳的影响,其中,地球非球形引力取8×8阶,月球和太阳星历均采用DE405。
(1)

式中:rv分别为探测器的位置和速度矢量;μeμmμh分别为地球、月球和太阳的引力常数;rhrm分别为太阳和月球的位置;rhdrmd分别为探测器到太阳和月球的位置;R为地球引力摄动位函数。
2 设计约束 由于探测器运行在月球停泊轨道上,机动时刻探测器的位置无法改变,只能将速度和转移时间作为自由变量,故初始轨道参数有速度和转移时间。则初始设计参数选取为
(2)

对于直接再入大气的月地转移轨道,通常对再入时刻的地心距、再入角、再入倾角有严格要求,为匹配落点的地理位置,再入轨道面与回收场必须共面,可采用回收场位置与再入轨道角动量间夹角来描述,则目标终端参数选取为
(3)

式中:Re为再入时刻的地心距;γe为再入时刻的再入角;ie为再入时刻的再入倾角;φe为再入时刻回收场位置与再入轨道角动量的夹角。
3 伪状态理论 伪状态理论是Wilson[12]提出的一种考虑地月中心引力作用下轨道近似计算方法。
考虑地球、月球及探测器三体系统,认为探测器对地球和月球运动没有影响。给定探测器tI时刻相对地心的位置速度ΣI=(RI,VI),利用伪状态理论可以计算tk时刻相对地心的位置速度Σk=(Rk, Vk)和相对月心的位置速度σk=(rk,vk)。定义以tk时刻月球位置为中心的伪状态转换球(Pseudostate Transformation Sphere,PTS),半径为RPTS。假设RI在PTS外部,Rk在PTS内部,如图 1所示。
图 1 伪状态理论示意图 Fig. 1 Schemetic of pseudostate theory
图选项




具体计算步骤如下:
1) 地心二体段。根据tI时刻ΣI,按地心二体轨道求解tk时刻相对地心的位置速度ΣIX=(RIX,VIX)。
2) 线性回溯段。求解tk时刻相对月心的位置速度σkx=(rkx, vkx),由式(4)计算tc时刻相对月心的位置速度σcx=(rcx,vcx)。
(4)

式中:tc为探测器到达PTS边界的时刻,满足关系式tc < tk,且由约束来确定。
3) 月心二体段。根据tc时刻相对月心的位置速度σcx,按月心二体轨道求解tk时刻相对月心的位置速度σk,即可得到相对地心的位置速度Σk
4 初值设计 本节给出了月地转移轨道的初值设计方法。若机动时刻探测器的位置在近月点前,选择近月点对应的伪状态位置作为迭代变量,转移轨道分为两段处理,转移轨道在近月点前按月球二体轨道处理。若机动时刻探测器的位置在近月点后,选择机动时刻探测器对应的伪状态位置作为迭代变量。
定义机动时刻t0探测器在月球惯性系下位置为r0,近月点前转移时间Δt1,近月点后转移时间Δt2,伪状态位置为Rs
步骤1 计算t0时刻月球位置速度RmVm,令Rs= RmΔt1=0。
步骤2 定义回收场的经度为λl,纬度为?l,位置为Rlhe为地心段角动量,ke为地轴的单位方向。根据余弦定理,可得
(5)

式中:α的符号对应了升降轨方式,β的符号对应了轨道面的方位,决定了再入航程大小(见图 2)。
图 2 月地转移轨道的几何约束 Fig. 2 Geometric constraints of Moon-to-Earth transfer orbit
图选项




定义t0时刻Rs在地球固连系下的经度为Gs,再入时刻Rl在地球固连系下的经度为Gl=Gs+α+β,记G=G1λ1,则地心段的飞行时间为
(6)

式中:ωe为地球自转角速度;N为指定的天数。
Rs为返回轨道地球段的起始点,探测器到达再入点时地心距为Re,再入角为γe,采用Luo等[13]提出的飞行路径角约束的Lambert算法可以求解该类边值问题,这里不再赘述。
, 地球段转移角为A1,则再入点位置为
(7)

再入点速度为
(8)

Rs处的飞行角为γs,速度大小为VsVs由式(9)计算:
(9)

步骤3 令探测器在PTS边界处相对月球的速度为vs,位置为rs,其大小为PTS半径,在月球无穷远处的速度大小为
(10)

(11)

轨道半长轴为
(12)

定义Δf为初始位置r0v间的夹角,由几何关系可知:
(13)

式中:f为月球无穷远处真近点角;f0为机动时刻的真近点角。
em为偏心率,存在以下关系:
(14)

,则
(15)

根据式(13)~式(15),可以求解f0em。令κ=-am/r0,则
(16)

(17)

Δ=2arcsin(1/em),,则探测器在近月点的位置和速度方向分别为
(18)

为完成月地转移,机动时刻探测器速度应为
(19)

步骤4 令,则PTS边界处
(20)

f0≥0,根据开普勒飞行时间公式计算r0rs的转移时间ts,否则转向步骤5。
步骤5 此时机动时刻探测器位置在近月点前,根据开普勒飞行时间公式分别计算r0rp的转移时间Δt1rprs的转移时间ts,同时更新Rm为近月点时刻的月球位置。
步骤6 更新Rs。根据伪状态理论得到
(21)

如果大于给定阈值,则令Rs= Rs,返回步骤2。否则迭代结束,将速度v 0作为初值。
5 精确解求解 初始设计参数P和目标终端参数Q存在确定的函数关系为
(22)

QP的误差之间的关系可线性近似为
(23)

为提高求解算法的速度,将轨道积分停止条件设置为到达给定再入角γe,则初始设计参数P中剔除了转移时间Δt,目标终端参数Q中解除了再入角γe的约束,则月地转移轨道设计降价为一个三维的两点边值问题。
X =[r, v]T,则式(23)可以重写为
(24)

当积分停止条件设置为到达给定再入角,则存在以下关系:
(25)

(26)

式(27)给出了消除时间变量后偏差传递矩阵形式:
(27)

式中:除矩阵外均有确定的形式,需要注意的是,对目标终端参数cos φe偏导数计算时,需要考虑由于再入时间变化引起的返回场惯性系下位置变化。由于三体轨道动力学特性导致矩阵没有解析表达式,本文采用文献[2-3, 14-15]提出的方法求解。
得到误差传递矩阵后,可利用各种迭代算法搜索精确解,本文采用微分修正方法求解精确解。
6 算例与分析 本节首先给出了一个月地转移轨道设计算例。相关算法均采用C++编程,程序运行环境:CPU为Intel Core 2.53G。
假设探测器在2022年1月1日进行单脉冲机动,此时探测器在月球惯性系下位置为[1 937.4, 0, 0]km,转移时间为3~4d,再入点地心高度为120km,再入倾角为45°,再入角为-6°,返回场的经度为110°,纬度为40°。
根据第4节的初值设计方法,存在4条月地转移轨道满足给定约束条件,下面以其中一条为例,表 1给出了伪状态位置误差的3个分量迭代过程。在初值设计时,迭代6次,耗时1 ms,在精确值求解时,迭代4次,耗时3s, 即可满足设定精度要求,表 2给出了再入点状态偏差的迭代过程。比较初值[-337.339, 2218.060, 914.856]m/s与精确解[-346.252, 2219.854, 914.926]m/s,两者模值的差为2.925m/s,两者矢量差的模值为9.092m/s,为初始速度大小的0.37%,在任务分析阶段,可用初值代替精确解,为任务分析提供参考。
表 1 初值设计时伪状态位置误差迭代过程 Table 1 Iteration process of pseudostate position error during initial solution design
迭代次数 Δx/m Δy/m Δz/m
1 7101140 -8735400 -3641240
2 453967 101658 104456
3 1622.07 279.576 393.42
4 102.868 13.296 21.020
5 5.799 0.804 1.249
6 0.350 0.0425 0.069


表选项






表 2 精确解求解时再入点状态偏差迭代过程 Table 2 Iteration process of reentry point state error during exact solution search
迭代次数 ΔRe/m Δcos ie Δcos φe
1 413183 0.0112 0.0007
2 6385.410 0.0019 0.0015
3 321.208 0.00012 0
4 2.334 0 0


表选项






表 3给出满足给定约束条件的4条月地转移轨道的关键特征参数,图 3图 4给出相应的3D和2D示意图。分析结果可知,4条月地转移轨道通过调整转移时间满足了再入轨道与返回场的共面要求,再入点速度大小基本相同,再入航程因轨道方位不同而出现差别。
图 3 月地转移轨道3D示意图 Fig. 3 Three-dimensional illustration of Moon-to- Earth transfer orbit
图选项




图 4 月地转移轨道2D示意图 Fig. 4 Two-dimensional illustration of Moon-to- Earth transfer orbit
图选项




表 3 月地转移轨道特征参数 Table 3 Characteristic parameters of Moon-to-Earth transfer orbit
编号 再入
方式
转移
时间/d
再入点速度/
(m·s-1)
再入
航程/km
1 升轨 3.689 10.985 4969
2 升轨 3.871 10.984 10378
3 降轨 3.334 10.987 32748
4 降轨 3.516 10.985 38124


表选项






下面计算2022年1月1日到2022年1月6日从月球停泊轨道出发的月地返回窗口,其中,月球停泊轨道在2022年1月1日的具体参数:高度为200km,偏心率为0,倾角为22.4°,升交点赤经为0°,纬度幅角为0°,目标轨道参数与上述算例相同。其中2022年1月1日满足共面转移条件,这是月地转移能量最省方式之一,以此为基准分析速度增量变化趋势。由于月球停泊轨道周期较小,角速度较快,速度方向变化较大,在每个轨道周期内变轨速度增量变化剧烈,如图 5所示,在每个轨道周期内,存在一纬度幅角使变轨速度增量最小。
图 5 一个轨道周期内变轨速度增量变化趋势 Fig. 5 Track change velocity increment versus time in one orbit period
图选项




在窗口搜索时选取每个轨道周期内最小变轨速度增量,分析其变化趋势。由图 6可知,结果有明显的周期性。在每一天内变轨速度增量逐渐增大,变轨速度增量差别最大为40m/s。随着探测器出发时刻延迟,月球停泊轨道与月地转移轨道面夹角增大,如图 7所示。若实现月地转移,需要调整轨道平面,变轨速度增量逐渐增大,再入航程同时逐渐增大,若月球出发时间延迟6 d,变轨速度增量增加130m/s。
图 6 变轨速度增量变化趋势 Fig. 6 Track change velocity increment versus time
图选项




图 7 轨道平面改变量变化趋势 Fig. 7 Change of orbit plane versus time
图选项




通过对月地返回窗口搜索可知,单脉冲转移方式通过调整轨道面,增加了月地返回机会,保证任务顺利进行。但随着对月球探测的继续深入,月球探测任务更加复杂,对探测器中止任务能力的需求相应增加,显然单脉冲转移方式在轨道平面调整量较大时,需要消耗更多能量,无法满足任务的要求,因此对多脉冲月地转移轨道的研究是后续工作的重点。
7 结论 本文给出了从月球停泊轨道出发直接再入大气的月地转移轨道设计方法。
1) 初值设计算法精度高,耗时少,在任务分析阶段可作为代替精确解,为任务分析提供参考依据。精确解求解算法求解效率高,可以作为设计月地转移轨道及搜索月地返回窗口的一个有力工具。
2) 给出了满足给定约束的4条月地转移轨道,并对比了关键特征参数,为月地转移轨道的选择提供依据。
3) 通过月地返回窗口分析可知,由于匹配返回场,每天内变轨速度增量逐渐增大。当月球停泊轨道与月地转移轨道面夹角逐渐增大时,使用单脉冲机动,则轨道面调整量增大,导致燃料消耗增大,限制了月地转移机会。

参考文献
[1] 杨维廉, 周文艳. 嫦娥一号月球探测卫星轨道设计[J]. 航天器工程, 2007, 16(6): 16-24.
YANG W L, ZHOU W Y. Orbit design for lunar exploration satellite CE1[J]. Spacecraft Engineering, 2007, 16(6): 16-24. (in Chinese)
[2] 李立涛, 杨涤, 崔祜涛. 奔月轨道的一种求解方法[J]. 宇航学报, 2003, 24(2): 150-155.
LI L T, YANG D, CUI H T. An algorithm for the solution of trans-lunar trajectory[J]. Journal of Astronautics, 2003, 24(2): 150-155. (in Chinese)
[3] 周文艳, 杨维廉. 月球探测器转移轨道的中途修正[J]. 宇航学报, 2004, 25(1): 89-92.
ZHOU W Y, YANG W L. Mid-correction of trans-lunar trajectory[J]. Journal of Astronautics, 2004, 25(1): 89-92. (in Chinese)
[4] OCAMPO C, SAUDEMONT R. Initial trajectory model for amulti-maneuver Moon-to-Earth abort sequence[J]. Journal of Guidance, Control, and Dynamics, 2010, 33(4): 1184-1194.
[5] 高玉东, 郗晓宁, 白玉铸, 等. 月球探测器返回轨道快速搜索设计[J]. 宇航学报, 2008, 29(3): 765-771.
GAO Y D, XI X N, BAI Y Z, et al. Fast searching design method for return transfer trajectory of lunar probe[J]. Journal of Astronautics, 2008, 29(3): 765-771. (in Chinese)
[6] 张磊, 于登云, 张熇. 直接再入大气的月球返回轨道设计研究[J]. 航天器工程, 2010, 19(5): 50-55.
ZHANG L, YU D Y, ZHANG H. Design of moon return trajectory with direct atmospheric reentry[J]. Spacecraft Engineering, 2010, 19(5): 50-55. (in Chinese)
[7] 张磊, 于登云, 张熇. 月地转移轨道快速设计与特性分析[J]. 中国空间科学技术, 2011, 31(3): 62-70.
ZHANG L, YU D Y, ZHANG H. Preliminary design and characteristic analysis of Moon-to-Earth transfer trajectories[J]. Chinese Space Science and Technology, 2011, 31(3): 62-70. (in Chinese)
[8] 郑爱武, 周建平, 刘勇. 月地转移轨道快速设计方法[J]. 北京航空航天大学学报, 2014, 40(3): 344-349.
ZHENG A W, ZHOU J P, LIU Y. Fast design method of Moon-to-Earth transfer trajectory[J]. Journal of Beijing University of Aeronautics and Astronautics, 2014, 40(3): 344-349. (in Chinese)
[9] 郑爱武, 周建平. 月地转移轨道精确轨道设计[J]. 飞行器测控学报, 2014, 33(1): 52-58.
ZHENG A W, ZHOU J P. Accurate design of Moon-to-Earth transfer orbit[J]. Journal of Spacecraft TT&C Technology, 2014, 33(1): 52-58. (in Chinese)
[10] 郑爱武, 周建平. 直接再入大气的月地返回窗口搜索策略[J]. 航空学报, 2014, 35(8): 2243-2250.
ZHENG A W, ZHOU J P. A search strategy of back windows for Moon-to-Earth trajectories directly returning to the earth[J]. Acta Aeronautica et Astronautica Sinica, 2014, 35(8): 2243-2250. (in Chinese)
[11] 汪中生, 孟占峰, 高珊, 等. 探测器月地转移入射变轨策略优化设计[J]. 航天器工程, 2019, 28(1): 10-18.
WANG Z S, MENG Z F, GAO S, et al. Optimization design of spacecraft trans-earth insertion strategy[J]. Spacecraft Engineering, 2019, 28(1): 10-18. (in Chinese)
[12] WILSON S W.A pseudostate theory for the approximation of three-body trajectories: AIAA-70-1061[R].Reston: AIAA, 1970. https://www.researchgate.net/publication/23621131_A_pseudostate_theory_for_the_approximation_of_three-body_trajectories
[13] LUO Q, MENG Z, HAN C. Solution algorithm to a quasi-Lambert's problem with fixed flight-direction angle constraint[J]. Celestial Mechanics and Dynamical Astronomy, 2011, 108(4): 409-427.
[14] ZIMMER S, OCAMPO C. Analytical gradients for gravity assist trajectories using constant specific impulse engines[J]. Journal of Guidance, Control, and Dynamics, 2005, 28(4): 753-760.
[15] OCAMPO C, MUNOZ J P. Variational model for the optimization of constrained finite-burn escape sequences[J]. Journal of Guidance, Control, and Dynamics, 2010, 33(5): 1615-1621.


相关话题/设计 地球 计算 大气 过程

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • 基于矩独立重要度的电路系统容错设计方法*
    随着电路系统在工业控制、国防军工、日常生活等领域的应用越发广泛,各领域对其可靠性的要求越来越高,因此电路系统的高可靠性设计也越发重要[1]。容错设计作为提高电路系统可靠性的一种重要方法[1],在航天器电源电路[1-2]、核电站控制电路[1,3]等要求高可靠性的领域已经发挥了重要作用。目前,电路系统容 ...
    本站小编 Free考研考试 2021-12-25
  • 近程动态范围激光雷达测距系统设计及误差分析*
    激光雷达作为一种新兴的主动式遥感技术,能够快速、准确、实时地获取地面目标的三维空间信息,近年来得到了极大的发展和应用。一个完整的激光雷达系统由激光测距仪、动态差分GPS接收机、惯性导航系统和成像装置组成。而激光测距仪是系统的核心部分,其采用的测距方法也有所不同。常用的方法有:脉冲法、相位法、三角法、 ...
    本站小编 Free考研考试 2021-12-25
  • 几何不确定性区间分析及鲁棒气动优化设计*
    航空航天工业中,飞行器设计的各参数往往被看作确定的量,但实际上,工程中是无法避免不确定性存在的。例如飞行器结构设计由于制造水平、测量误差等会造成参数的不确定性,飞行器飞行环境及载荷作用也会存在更强的不确定性等[1]。不确定性因素严重影响了与其相对的确定性设计的精度和可靠度,因此,想要获得更可靠更稳定 ...
    本站小编 Free考研考试 2021-12-25
  • 发动机短舱泄压过程瞬态仿真*
    中国民用航空规章第25部《运输类飞机适航标准》CCAR25.1103(d)条款规定[1]:对于涡轮发动机和辅助动力装置的引气导管系统,如在空气导管引气口与飞机用气装置之间任意部位上的导管破损,不得造成危害。发动机引气管道发生破裂或泄漏时,迅速上升的短舱内部压力可能会破坏短舱结构甚至损坏发动机,因此需 ...
    本站小编 Free考研考试 2021-12-25
  • 高峰时段下离港航空器绿色滑行策略设计与评价*
    近年来,大型枢纽机场航班数量迅速增加,特别是高峰时段,机场容量受限,场面拥堵严重,导致航空器“时走时停”或是在跑道端长时间排队等待,造成大量额外燃油消耗与资源浪费的同时,也产生了大量不必要的污染物排放,极不利于绿色民航的发展。航空运输业产生的排放污染物主要包括碳氢(HC)、一氧化碳(CO)和氮氧化物 ...
    本站小编 Free考研考试 2021-12-25
  • 液体火箭发动机故障诊断器设计及其HIL验证*
    液体火箭发动机是在极端物理条件下运行的复杂热力学系统,其故障的发生和发展具有极端的快速性和极大的破坏性。液体火箭对发动机的安全性和可靠性提出了更高的要求,发展可靠有效的发动机故障诊断系统,可以对发动机工作过程中出现的故障予以预警和判断,并能及时采取有效措施,保护液体火箭和载荷安全,有效避免由于发动机 ...
    本站小编 Free考研考试 2021-12-25
  • 一种区域参考大气密度的建模与应用方法*
    参考大气模型在空天试验中有着重要的应用[1-3]。它表征了大气的理想气体性质,包括压力、密度、温度、各元素分子量等随时间、空间的变化情况。模型的输入可以是高度、纬度、或者在一年中所处的天数。不同的参考大气模型有不同的输入和输出。例如美国标准大气模型(USSA—1962)以及其后续型号(USSA—19 ...
    本站小编 Free考研考试 2021-12-25
  • 机动发射条件下空间飞行器上升段弹道设计*
    在现代战争条件下,精确制导武器发展迅速,可以实现“点对点”打击,导致提高武器系统的生存能力变得十分重要。与传统的固定阵地发射方式相比,机动发射在提高生存能力的同时带来了打击的突然性,可以有效提高武器系统作战效能。同样,在应急状态下,能够实现卫星的机动发射,可以迅速对报废卫星进行补充,有效支持作战行动 ...
    本站小编 Free考研考试 2021-12-25
  • 基于深度学习的无人机数据链信噪比估计算法*
    地-空数据链作为无人机系统的重要组成部分,发挥着发送上行遥控指令和回传遥测侦察信息等重要作用[1]。信噪比(Signal-to-NoiseRatio,SNR)是评价无人机通信系统信道环境和通信质量的重要指标,精确的信噪比估计既可以为无人机数据链提供功率控制、信道分配所需要的信息,又可以促使数据链系统 ...
    本站小编 Free考研考试 2021-12-25
  • 空间光学遥感器真空热试验工装模块化设计*
    空间光学遥感器是搭载在卫星或航天飞行器上利用光学技术获取目标属性和相关信息的精密科学设备。为了适应光学遥感技术的快速发展,对光学遥感器的要求逐渐向高空间分辨率、高光谱分辨率、高辐射分辨率性能及轻量化的方向发展。而光学遥感器作为卫星有效载荷,对环境温度均匀性的要求将变得更加苛刻,均匀的温度场可为其提供 ...
    本站小编 Free考研考试 2021-12-25