一种调轨方案自动化设计方法技术

本发明专利技术提供了一种调轨方案自动化设计方法，包括：1)、设置计算条件；2)、基于扁率和大气摄动递推轨道至指定时刻，递推后对半长轴的衰减量和由此产生的相位漂移量进行补偿；3)、根据两星轨道参数生成调轨窗口；4)、根据两星轨道参数判定流程号，选择主流程；主流程分支1为顺行漂移分支，进行匹配时机预估，结合上述计算结果判定后续方案号并执行对应的调轨方案；主流程分支2为逆向漂移分支，执行逆向漂移方案，与分支1同时生成两套并行方案；5)：对生成的方案进行评估，用以判定流程结束或重新开始。本发明专利技术方法快速可靠，提高决策效率。

【技术实现步骤摘要】
一种调轨方案自动化设计方法
本专利技术涉及卫星轨道设计
，具体地，涉及一种近地圆轨道的两星联合使用相位差初始化设计方法。
技术介绍
两星或多星的联合使用是完成空间任务的一种非常有效的工作模式，可以达到单科卫星独立工作所难以达到的性能。合理设计联合使用方式可以提高时间分别率、加大覆盖范围，达到改善覆盖特性的目的。近地圆轨道的两星联合使用对两星之间的轨道关系提出了一定的要求，即两星构型应满足星座对地最大重访时间不超过要求值。为确保实现这一要求，两星轨道需要满足：(1)具有良好的星下点重复特性，这对两星相位差提出了要求；(2)为避免两星相位差呈现发散趋势，两星轨道应长期处在同一高度。实际工程中，往往是在一颗星(以下称“01星”)已在轨的情况下，发射第二颗星(以下称“02星”)。02星实际入轨后，两星高度不同是难以避免的，初始相位差一般不会正好等于标称设计值。因此，在02星入轨后，为达到相位匹配的轨道调整来对两星构型进行初始化是必要的。对于可随时进行速度切向变轨的卫星，两星构型初始化过程较为简便，在入轨后可得两星标称相位差，待漂移至标称值时(下文称“相位匹配”)调整两星轨道至同一高度即完成调整。但对于轨道机动能力受约束的卫星，其初始化过程的设计相对较复杂。本专利技术以变轨能力受β角约束为例进行说明。不同的入轨时间，轨道的β角不同，并按规律呈现变化，在一些β角区间内，卫星处于正飞状态，仅可提供正切向推力，即仅可向上调轨，在另一些区间内，卫星偏航机动180°，只可提供负切向，即仅可向下调轨，β角处于其它情况不可调轨。由此，产生调轨窗口的概念，两星可调轨的时段在轨道状态已知的情况下被分成若干区段。将调轨后两星轨道高度差变小的称为窗口A，高度差变大的称为窗口B，当预估的相位匹配时间不落于窗口A时，单次调轨不能直接完成初始化，则初始化方案必包含两次或两次以上的调轨。在这种情况下，两次调轨时机和调轨量的分配需要进行具体设计。02星发射日期的不确性定带来了相位差待漂移量的不确定性、入轨高度的不确定性带来了漂移速度甚至漂移方向的不确定性，加之调轨窗口的离散分布星使得问题呈现出多种可能，对于每一种可能，必须一一加以区分并给出不同的解决方案。
技术实现思路
本专利技术提供了一种针对两星联合使用构型初始化要求的调轨方案自动化设计方法，在调轨窗口受到显著约束的情况下，为了提高设计效率，适应02星入轨后快速设计快速决策的需求，免去02星入轨前因输入条件不确定性造成的大量繁琐设计，在两星任意初始轨道状态下均可提供构型初始化的解决方案。为达到上述目的，本专利技术所采用的技术方案如下：一种调轨方案自动化设计方法，针对两星联合使用构型初始化的要求，包括如下：1)、设置计算条件，包括，地球扁率参数、大气摄动造成的轨道高度下降速度，仿真时间；2)、基于扁率和大气摄动递推轨道至指定时刻，输入01星轨道状态并递推至历元时刻，采用映射J4与轨道半长轴线性衰减结合的方式，在J4递推后对半长轴的衰减量和由此产生的相位漂移量进行补偿；3)、根据两星轨道参数生成调轨窗口，β角对卫星姿态模式具有约束，继而对调轨窗口有约束，以历元时刻两星轨道参数出发，递推数十天，可得基于时间轴的β角数据，继而得到调轨窗口；4)、根据两星轨道参数判定流程号，选择主流程；主流程分支1为顺行漂移分支，进行匹配时机预估，结合上述计算结果判定后续方案号并执行对应的调轨方案；主流程分支2为逆向漂移分支，执行逆向漂移方案，与分支1同时生成两套并行方案；5)：对4)生成的方案进行评估，用以判定流程结束或重新开始。所述2)中，基于扁率和大气摄动递推轨道，在非球型摄动解析递推法的基础上加入了轨道高度变化的考虑，能够有效补偿在解析法递推中因轨道高度变化带来的显著相位误差，确保递推精度要求和快速性要求，所述补偿按下式进行：式中，t为递推时长，a0,e0,i0,Ω0,ω0,M0是初始时刻轨道根数；映射J4采用拟平根递推法；a2是半长轴二阶量，以轨道高度衰减率代入；M2是平近点角变化二阶量，即相位变化补偿量。所述3)中，调轨窗口由太阳面光照角进行判定，符合一定区间时，仅可向上调轨，符合另一区间时，仅可向下调轨，其它情况不可调轨，即不属于调轨窗口，将使相位差漂移率变慢的称为窗口A，反之称为窗口B。所述4)中，根据两星的轨道参数进行自动判定，将主流程结构分为两个层次：其一，根据两星轨道高度判定使用主流程1或是主流程2；其二，进入任一主流程中并行运行两个分支；顺行分支中进行匹配时机预估，并依判定结果执行预置方案；逆向漂移分支直接执行逆向漂移措施和预置方案；同时输出两个并行方案供后续环节评估；其中，标称相位差的计算由下式表达：式中，days为02星轨迹落后于01星轨迹的天数，ωe为地球自转速度，和分别为地球非球形摄动引起的相应轨道根数的变化率。匹配时机预估方法是：以当前两星相对相位漂移率预估其漂移至指定标称相位差为匹配的时机，其中涉及的标称相位差目标值依据两星相位匹配要求计算而得，计算所需的摄动参数和重访天数按需求灵活设置。所述调轨方案为适应在指定窗口区段达到匹配要求的超调设计值和漂移设计值，形成基于两次轨控的调轨方案，合理分配两次轨控调整量可达到前后双向调整匹配时机，并使匹配时机落于理想调轨窗口内。所述逆向漂移方案在主流程的分支2中使用，是两星相位朝与初始状态相反的方向漂移，与分支1并行使用并一同输出，供后续环节评估。本专利技术提供的一种针对两星联合使用构型初始化要求的调轨方案自动化设计方法，解决了调轨窗口受到显著约束情况下的两星相位匹配的两次以上的机动方案设计问题，按照输入计算过程量，逐级对问题进行分类讨论，通过设计合理的轨控时机和轨控量，使各次轨控均位于调轨窗口内并满足构型要求；能够实现自动化设计，解决了输入不确定性带来的复杂设计问题，便于在入轨后飞控组前方人员和资源有限的情况下，进行当场快速设计快速决策难的问题。附图说明图1为本专利技术方法流程框图；图2为本专利技术中调轨窗口示意图；图3为本专利技术中相位差标称值设计示意图。具体实施方式下面结合具体实施例对本专利技术进行详细说明。在本实施例中，如图1所示，本专利技术包括如下步骤：环节1：设置设计条件设置计算所需的条件包括，地球扁率参数、大气摄动造成的轨道高度下降速度，仿真时间。环节2：输入轨道状态并递推一般情况下，由于采集到的两星轨道历元不一致，这里以01星为例，递推其轨道至与02星相同的历元，便于后续环节的计算。输入01星轨道状态并递推至历元时刻，采用映射J4与轨道半长轴线性衰减结合的方式，在J4递推后对半长轴的衰减量和由此产生的相位漂移量按下式进行补偿。式中，t为递推时长，a0,e0,i0,Ω0,ω0,M0是初始时刻轨道根数；映射J4采用拟平根递推法，具有精度高、速度快的优点，此处不再详述；a2是半长轴二阶量，以轨道高度衰减率代入；M2是平近点角变化二阶量，即相位变化补偿量。环节3：生成窗口定义β角为轨道面法向量与太阳矢量夹角的余角，其对卫星姿态模式具有约束，继而对调轨窗口有约束，以初始情况下02星高于01星为例，两星调轨窗口定义如下：窗口A：01星向上或02星向下，02星轨道高于01星时减小相位漂移率；窗口B：02星向上或01星向下，02星轨道高于01星时增大相位漂移率。以历本文档来自技高网...

【技术保护点】
一种调轨方案自动化设计方法，针对01、02两星联合使用构型初始化的要求，其特征在于，包括如下：1)、设置计算条件，包括，地球扁率参数、大气摄动造成的轨道高度下降速度，仿真时间；2)、基于扁率和大气摄动递推轨道至指定时刻，输入01星轨道状态并递推至历元时刻，采用映射J4与轨道半长轴线性衰减结合的方式，在J4递推后对半长轴的衰减量和由此产生的相位漂移量进行补偿；3)、根据两星轨道参数生成调轨窗口，β角对卫星姿态模式具有约束，继而对调轨窗口有约束，以历元时刻两星轨道参数出发，递推数十天，可得基于时间轴的β角数据，继而得到调轨窗口；4)、根据两星轨道参数判定流程号，选择主流程；主流程分支1为顺行漂移分支，进行匹配时机预估，结合上述计算结果判定后续方案号并执行对应的调轨方案；主流程分支2为逆向漂移分支，执行逆向漂移方案，与分支1同时生成两套并行方案；5)：对4)生成的方案进行评估，用以判定流程结束或重新开始。

【技术特征摘要】
1.一种调轨方案自动化设计方法，针对01、02两星联合使用构型初始化的要求，其特征在于，包括如下：1)、设置计算条件，包括，地球扁率参数、大气摄动造成的轨道高度下降速度，仿真时间；2)、基于扁率和大气摄动递推轨道至指定时刻，输入01星轨道状态并递推至历元时刻，采用映射J4与轨道半长轴线性衰减结合的方式，在J4递推后对半长轴的衰减量和由此产生的相位漂移量进行补偿；3)、根据两星轨道参数生成调轨窗口，β角对卫星姿态模式具有约束，继而对调轨窗口有约束，以历元时刻两星轨道参数出发，递推数十天，可得基于时间轴的β角数据，继而得到调轨窗口；4)、根据两星轨道参数判定流程号，选择主流程；主流程分支1为顺行漂移分支，进行匹配时机预估，结合上述计算结果判定后续方案号并执行对应的调轨方案；主流程分支2为逆向漂移分支，执行逆向漂移方案，与分支1同时生成两套并行方案；5)：对4)生成的方案进行评估，用以判定流程结束或重新开始。2.根据权利要求1所述的调轨方案自动化设计方法，其特征在于，所述2)中，基于扁率和大气摄动递推轨道，在非球型摄动解析递推法的基础上加入了轨道高度变化的考虑，能够有效补偿在解析法递推中因轨道高度变化带来的显著相位误差，确保递推精度要求和快速性要求，所述补偿按下式进行：式中，t为递推时长，a0,e0,i0,Ω0,ω0,M0是初始时刻轨道根数；映射J4采用拟平根递推法；a2是半长轴二阶量，以轨道高度衰减率代入；M2是平近点角变化二阶量，即相位变化补偿量。3.根据权利要求1所述的调轨方案自动化设计方法，其特征在于，所述3)中，调轨窗口由太阳面光照角进行判定，符合一定区间时，仅可向上调轨，符合另一区间时，仅可向下调轨，其它情况不可调轨，即不属于调轨窗口，将使相位差漂移率变慢的称为窗口A，反之称为窗口B。4.根据权利要求1所述的调轨方案自动化设计方法，其特征在于，所述4...

【专利技术属性】
技术研发人员：许哲，沈庆丰，叶小舟，黄欣，李绿萍，李楠，陆启省，李鉴，邓武东，袁双，
申请(专利权)人：上海卫星工程研究所，
类型：发明
国别省市：上海,31