法律状态公告日
法律状态信息
法律状态
2018-10-23
未缴年费专利权终止 IPC(主分类):G06F17/50 授权公告日:20170714 终止日期:20171029 申请日:20141029
专利权的终止
2017-07-14
授权
授权
2015-04-29
实质审查的生效 IPC(主分类):G06F17/50 申请日:20141029
实质审查的生效
2015-04-01
公开
公开
技术领域:
本发明涉及飞行器轨道设计技术领域,具体涉及一种飞船返回预定落点回归轨道设计方 法。
背景技术
为了保证载人飞船每天能够为航天员提供返回着陆国内预定着陆场的机会,我国载人航 天工程载人飞船轨道设计一直遵循回归轨道的设计原则,一期工程载人飞行和二期工程交会 对接的载人飞船采用的是2天回归轨道。三期空间实验室和空间站工程中,载人飞船运行轨 道高度和轨道倾角与一、二期工程相比发生了变化。三期工程飞船拟采用3天回归轨道,不 再通过轨道机动控制对空间实验室和空间站相对飞船的初始相位进行调整。在轨道倾角和轨 道高度在一定范围内变化、轨道初始相位不确定的条件下,飞船轨道的回归特性和返回特性 将发生较大的改变。为了保证三期工程载人飞船仍能返回预定着陆场,需进行基于准确返回 预定着陆场的三天回归轨道设计,为三期工程载人飞船返回轨道设计和控制实施提供依据。 因此,飞船返回预定落点回归轨道设计是空间实验室和空间站工程轨道设计的技术要点之一, 目前尚没有相关的轨道设计方法。
发明内容:
本发明需解决技术问题在于提供一种飞船返回预定落点回归轨道设计方法。
为解决上述技术问题,本发明提供的一种飞船返回预定落点回归轨道设计方法包括如下 步骤:
1)设置初始轨道参数,包括轨道半长轴初值a0、轨道偏心率初值e0、轨道倾角初值i0、 升交点经度初值D0、轨道近地点幅角初值ω0、轨道平近点角初值M0、飞船质量和迎风面积、 轨道力模型,轨道力模型包括:地球中心引力、地球形状引力摄动;
2)设置轨道设计条件,包括:标称落点经度L0、纬度B0;初始返回圈号qr;轨道回归 周期天数i;轨道回归周期圈数ki;落点经度收敛门限εL;回归圈升交点经度收敛门限ελn;
3)数值积分计算各圈纬度为B0的星下点经度L,计算各圈升交点经度λn;
4)计算返回圈星下点经度差ΔL=L0-L(qr),L(qr)为初始返回圈号qr星下点纬度为B0的 经度;计算回归圈升交点经度差Δλn=λ(qr+ki)-λ(qr),λ(qr)为初始返回圈号qr的升交点经 度,λ(qr+ki)为初始返回圈号qr+ki的升交点经度;
5)根据下列公式由返回圈星下点经度差ΔL计算轨道平近点角偏差ΔM;
其中
为航天器轨道平面在惯性空间进动的平均角速度,μ为地球引力常数,Re为地球赤道 参考半径,a为轨道半长轴,i为轨道倾角,e为轨道偏心率,n为轨道平均角速度,ωe为 地球自转角速度;
6)修正轨道平近点角,M=M0+ΔM;
7)数值积分计算各圈纬度为B0的星下点经度L,计算各圈升交点经度λn;
8)计算返回圈星下点经度差ΔL=L0-L(qr),计算回归圈升交点经度差 Δλn=λ(qr+ki)-λ(qr);
9)判断返回圈星下点经度差是否满足ΔL<εL;若不满足,则返回到5),重复5)-9); 若满足,则转到10);
10)根据下列公式由回归圈升交点经度差Δλn计算轨道半长轴修正量Δa;
其中,
ki为回归周期圈数,i=1,2,3,k1=k;k2=2k+1;k3=3k+1和3k+2,ωe为地球自转角速度, 为航天器轨道平面在惯性空间进动的平均角速度,J2为地球引力摄动项,ω为轨道近地点 幅角;
11)修正轨道半长轴,a=a0+Δa;
12)数值积分计算各圈纬度为B0的星下点经度L,计算各圈升交点经度λn;
13)计算返回圈星下点经度差ΔL=L0-L(qr),计算回归圈升交点经度差 Δλn=λ(qr+ki)-λ(qr);
14)同时判断返回圈星下点经度差是否满足ΔL<εL和回归圈升交点经度差是否满足 Δλn<ελn;若不满足,则返回到5),重复5)-14);若满足,则转到15);
15)计算结束,得到满足返回预定落点和回归特性的设计轨道参数轨道半长轴a和轨道 平近点角M,即满足返回预定落点和回归特性的设计轨道根数为a,e0,i0,D0,ω0,M。
本发明提出了一种飞船返回预定落点回归轨道设计方法,通过选择轨道平近点角实现星 下点准确经过预定落点,通过调整轨道半长轴使轨道满足回归特性,建立轨道半长轴和轨道 平近点角2个参数双层迭代求解流程,得到轨道倾角、轨道半长轴和轨道平近点角相互匹配 的设计参数,保证飞船星下点轨迹每回归周期准确经过预定落点。
附图说明:
图1是本发明飞船返回预定落点回归轨道设计的计算流程图。
图2是实例3天回归轨道初始返回圈为15圈的设计轨道参数关系图。
图3是实例3天回归轨道初始返回圈为30圈的设计轨道参数关系图。
图4是实例3天回归轨道初始返回圈为46圈的设计轨道参数关系图。
图5是实例2天回归轨道初始返回圈为15圈的设计轨道参数关系图。
图6是实例2天回归轨道初始返回圈为31圈的设计轨道参数关系图。
具体实施方式:
一种飞船返回预定落点回归轨道设计方法,是通过选择轨道平近点角实现星下点准确经 过预定落点,通过调整轨道半长轴使轨道满足回归特性,建立轨道半长轴和轨道平近点角2 个参数双层迭代求解流程,得到轨道倾角、轨道半长轴和轨道平近点角相互匹配的设计参数, 保证飞船星下点轨迹每回归周期准确经过预定落点。
该方法包括如下步骤:
步骤一、建立星下点经过预定落点的算法
考虑地球引力摄动J2项的影响,航天器轨道平面在惯性空间进动的平均角速度为:
式中,μ为地球引力常数,Re为地球赤道参考半径,a为轨道半长轴,i为轨道倾角,e为轨 道偏心率;
航天器星下点轨迹可用下述方程描述:
Δλ=arctan(cos i tan u) (3)
式中,为地心纬度,L为经度,u为轨道纬度幅角,λ0为初始升交点经度,Δλ为轨道纬 度幅角u相距升交点的经度差,ωe为地球自转角速度,t为从初始升交点起算的时间;
式(4)两端求经度L对时间t的差分:
已知轨道平近点角M与轨道平均角速度n的关系为:
M=nt (6)
式(6)两端求轨道平近点角M对时间t的差分:
ΔM=n·Δt (7)
则式(5)可改写为:
因此得到:
式(9)表明,星下点经度偏差ΔL可通过调整轨道平近点角ΔM消除;
设轨道平近点角初值为M0,则消除星下点经度偏差ΔL的轨道平近点角M为:
M=M0+ΔM (10)
因此,通过选择轨道平近点角可实现星下点准确经过预定落点;当进行星下点轨迹设计 时,ΔL即为按照预定落点纬度为基准的星下点与预定落点的经度差,通过调整轨道平近点角 ΔM,使得ΔL=0,即得到了星下点轨迹经过预定落点的轨道平近点角M;
步骤二、建立回归轨道的算法
轨道每运行一圈升交点经度的西退量λn为:
式中,Tn为地球引力摄动J2项影响下的轨道交点周期;
式中,T为轨道周期,ω为轨道近地点幅角;
轨道周期与轨道半长轴的关系为:
式中,μ为地球引力常数;
定义:
则有:
Tn=CTn·T (15)
轨道运行1天后升交点经度再次回到初始升交点附近,相邻圈升交点的最宽间隔即为轨 道每运行一圈升交点经度的西退量;
设k为1天内星下点轨迹的整数圈次,k=int[2π/λn],1天回归轨道升交点经度间隔的 条件为:kλn=2π,k为1天回归轨道的回归周期圈数;同理,2天回归轨道升交点经度间隔 的条件为:(2k+1)λn=4π,2k+1为2天回归轨道的回归周期圈数;3天回归轨道升交点经 度间隔的条件为:(3k+1)λn=6π和(3k+2)λn=6π,3k+1和3k+2为3天回归轨道的回归周 期圈数;
设回归周期圈数为ki,统一1天-3天回归轨道升交点经度间隔的条件为:
当轨道交点周期存在偏差ΔTn时,则有:
式中,Δλn为相邻回归圈的升交点经度偏差;
整理式(17),得到:
式(15)两端求轨道交点周期偏差ΔTn对轨道周期T的差分:
ΔTn=CTn·ΔT (19)
则有:
式(13)两端求轨道半长轴a对轨道周期T的差分:
因此有:
式(22)表明,通过调整轨道半长轴Δa可消除相邻回归圈的升交点经度差Δλn;
设轨道半长轴初值为a0,则消除相邻回归圈的升交点经度差Δλn的轨道半长轴a为:
a=a0+Δa (23)
当进行回归轨道设计时,Δλn为回归轨道相邻回归圈的升交点经度差,通过调整轨道半 长轴Δa,使得Δλn=0,即得到了满足回归轨道特性的轨道半长轴a;
步骤三、建立轨道半长轴和轨道平近点角2个参数双层迭代求解流程
基于步骤一和步骤二,建立轨道半长轴和轨道平近点角2个参数的双层迭代求解流程, 精确求解飞船满足返回预定落点和回归特性的轨道参数,见图1。
1)初始轨道参数设置。初始轨道参数包括:轨道半长轴初值a0、轨道偏心率初值e0、轨 道倾角初值i0、升交点经度初值D0、轨道近地点幅角初值ω0、轨道平近点角初值M0、飞船 质量和迎风面积、轨道力模型等。轨道力模型包括:地球中心引力、地球形状引力摄动;
2)设置轨道设计条件。轨道设计条件包括:标称落点经度L0、纬度B0;初始返回圈号qr; 轨道回归周期天数i;轨道回归周期圈数ki;落点经度收敛门限εL;回归圈升交点经度收敛 门限ελn;
3)数值积分计算各圈纬度为B0的星下点经度L,计算各圈升交点经度λn;
4)计算返回圈星下点经度差ΔL=L0-L(qr),L(qr)为初始返回圈号qr星下点纬度为B0的 经度;计算回归圈升交点经度差Δλn=λ(qr+ki)-λ(qr),λ(qr)为初始返回圈号qr的升交点经 度,λ(qr+ki)为初始返回圈号qr+ki的升交点经度;
5)根据公式(9),由返回圈星下点经度差ΔL计算轨道平近点角修正量ΔM;
6)修正轨道平近点角,M=M0+ΔM;
7)数值积分计算各圈纬度为B0的星下点经度L,计算各圈升交点经度λn;
8)计算返回圈星下点经度差ΔL=L0-L(qr),计算回归圈升交点经度差 Δλn=λ(qr+ki)-λ(qr);
9)判断返回圈星下点经度差是否满足ΔL<εL;若不满足,则返回到5),重复5)-9); 若满足,则转到10);
10)根据公式(22),由回归圈升交点经度差Δλn计算轨道半长轴修正量Δa;
11)修正轨道半长轴,a=a0+Δa;
12)数值积分计算各圈纬度为B0的星下点经度L,计算各圈升交点经度λn;
13)计算返回圈星下点经度差ΔL=L0-L(qr),计算回归圈升交点经度差 Δλn=λ(qr+ki)-λ(qr);
14)同时判断返回圈星下点经度差是否满足ΔL<εL和回归圈升交点经度差是否满足 Δλn<ελn;若不满足,则返回到5),重复5)-14);若满足,则转到15);
15)计算结束,得到满足返回预定落点和回归特性的设计轨道参数轨道半长轴a和轨道 平近点角M,即满足返回预定落点和回归特性的设计轨道根数为a,e0,i0,D0,ω0,M。
为验证本发明所述方法的有效性,本具体实施方式给出了2个实例,分别进行了返回预 定落点的3天回归轨道参数设计和2天回归轨道参数设计。
载人飞船计算用轨道根数设定为:轨道偏心率e=0.0005,轨道升交点经度ΩG=350°,轨 道近地点幅角ω=130°,轨道倾角变化为:42.0°-43.0°;轨道摄动模型考虑地球非球形引力项 32×32阶次。
载人飞船预定落点设定为:经度L=110.0°,大地纬度B=42.0°。
对载人飞船轨道星下点轨迹经过预定落点进行3天回归和2天回归的轨道参数设计,求 解不同轨道倾角对应的回归轨道的半长轴和平近点角。
(1)3天回归轨道参数设计实例
3天回归轨道的回归周期为46圈,经过落点的初始返回圈可选为第1天的15圈、第2 天的30圈、第3天的46圈。
表1、表2、表3分别为初始返回圈为15圈、30圈、46圈的3天回归轨道设计结果。轨 道倾角i从42.0°变化至43.0°,对应每一个轨道倾角都有一组轨道半长轴和轨道平近点角的 设计参数与之对应,H为轨道半长轴对应的轨道高度,H=a-6371km,N为设计回归轨道的返 回圈升交点经度。
图2、图3、图4分别为初始返回圈为15圈、30圈、46圈的3天回归轨道设计参数关系 图。横坐标为轨道倾角,从42.0°到43.0°;左侧纵坐标为轨道半长轴,右侧纵坐标为轨道平 近点角。对应每一个轨道倾角都有一组轨道半长轴和轨道平近点角的设计参数与之对应。
表1 3天回归轨道初始返回圈为15圈的设计轨道参数
表2 3天回归轨道初始返回圈为30圈的设计轨道参数
表3 3天回归轨道初始返回圈为46圈的设计轨道参数
3天回归轨道参数设计实例结果表明:
按照不同天返回圈经过预定落点,设计的回归轨道参数不同。轨道倾角在42.0°-43.0°变 化时,初始返回圈为15圈,轨道半长轴变化近似为下凹的抛物线,最小值为6760.3km,最 大值为6768.3km,对应轨道高度389.3km到397.3km,变化范围约8km;轨道平近点角变化 近似为线性增大,从68.3°到212.9°,变化范围约145°。初始返回圈为30圈,轨道半长轴变 化近似为上凸的抛物线,最小值为6760.0km,最大值为6769.2km,对应轨道高度389.0km到 398.2km,变化范围约9km;轨道平近点角变化近似为线性增大,从318.2°到94.7°,变化范 围约137°。初始返回圈为46圈,轨道半长轴变化近似为先上凸后下凹的2组抛物线,最小 值为6760.5km,最大值为6768.8km,对应轨道高度389.5km到397.8km,变化范围约8km; 轨道平近点角变化近似为线性增大,从191.7°到337.0°,变化范围约145°。
轨道倾角在42.0°-43.0°变化时,3天回归轨道15圈返回的升交点经度范围为31.9°-41.5 °,30圈返回的升交点经度范围为32.6°-41.5°,46圈返回的升交点经度范围为32.3°-41.7°。 3天回归轨道不同返回圈的升交点经度范围大致一致,在31.9°-41.7°之间。
因此,返回预定落点的3天回归轨道,按照1、2、3不同天返回圈设计时,轨道半长轴 参数范围基本一致,在6760km到6769km约9km范围内变化,轨道高度在389km到398km 之间,但随着轨道倾角的变化趋势不同。对于1、2、3不同天返回圈,轨道平近点角参数将 360°近似平分为3个区间,区间边界相互重叠。
(2)2天回归轨道参数设计实例
2天回归轨道的回归周期为31圈,经过落点的返回圈可选为第1天的15圈、第2天的 31圈。
表4、表5分别为初始返回圈为15圈、31圈的2天回归轨道设计结果。轨道倾角i从42.0 °变化至43.0°,对应每一个轨道倾角都有一组轨道半长轴和轨道平近点角的设计参数与之对 应,H为轨道半长轴对应的轨道高度,H=a-6371km,N为设计回归轨道的返回圈升交点经度。
图5、图6分别为初始返回圈为15圈、31圈的2天回归轨道设计参数关系图。横坐标为 轨道倾角,从42.0°到43.0°;左侧纵坐标为轨道半长轴,右侧纵坐标为轨道平近点角。对应 每一个轨道倾角都有一组轨道半长轴和轨道平近点角的设计参数与之对应。
表4 2天回归轨道初始返回圈为15圈的设计轨道参数
表5 2天回归轨道初始返回圈为31圈的设计轨道参数
2天回归轨道参数设计实例结果表明:
轨道倾角在42.0°-43.0°变化时,初始返回圈为15圈和31圈,轨道半长轴变化规律基本 一致,轨道半长轴变化近似为先上凸后下凹的2组抛物线,最小值为6710.0km,最大值为 6718.6km,对应轨道高度339.0km到347.6km,变化范围约9km;初始返回圈为15圈和31 圈,轨道平近点角变化近似为线性增大。15圈为初始返回圈,轨道平近点角从18.1°到164.7 °,范围约147°;31圈为初始返回圈,轨道平近点角从198.2°到344.7°,范围约147°。2组 轨道平近点角区间没有覆盖360°,有2段约33°的空缺。
轨道倾角在42.0°-43.0°变化时,2天回归轨道15圈返回和31圈返回的升交点经度范围 基本一致,在32.2°-41.6°之间。
因此,返回预定落点的2天回归轨道,按照1、2不同天返回圈设计时,轨道半长轴参数 范围基本一致,在6710km到6719km约9km范围内变化,轨道高度在340km到348km之间, 而且随着轨道倾角的变化趋势相同。对于1、2不同天返回圈,轨道平近点角参数将360°近 似平分为2个区间,但区间边界有2段空缺。
3天回归轨道和2天回归轨道的参数设计实例结果表明,本发明提出的返回预定落点回 归轨道设计方法,通过选择轨道平近点角实现星下点轨迹准确经过预定落点、调整轨道半长 轴实现轨道回归特性,设计了轨道倾角、轨道半长轴和轨道平近点角相互匹配的回归轨道, 保证了飞船星下点轨迹每回归周期准确经过预定落点。
按照不同天返回,3天回归轨道存在1天返回、2天返回和3天返回的3组回归轨道,2 天回归轨道存在1天返回和2天返回的2组回归轨道。
不同天返回的回归轨道参数不同,但对于一个轨道倾角,轨道半长轴和轨道平近点角是 唯一匹配的。对于轨道倾角42.0°-43.0°,3天回归轨道的轨道高度在389km到398km之间; 对于1、2、3不同天返回圈,轨道平近点角范围略大于120°。2天回归轨道的轨道高度在330km 到348km之间;对于1、2不同天返回圈,轨道平近点角范围略小于180°。
3天回归轨道和2天回归轨道的参数设计实例验证了飞船返回预定落点回归轨道设计方 法,具有较好的工程应用价值。
以上所述的具体描述,对本发明的目的、技术方案和有益效果进行了详细说明,所应理 解的是,以上所述仅为本发明的具体实例而已,并不用于限定本发明的保护范围。
机译: 改变跟随飞机的道路的方法,飞机最初跟随预定的轨道,将飞机辅助返回到预定的轨道的方法。
机译: 轨道空间飞船和返回系统
机译: 一种设备,用于根据定义的事件并根据预定定律确保已遭受干扰的简单周期或复合,同步或频率在预定的频率和相位中返回