CN106802667B - 一种基于双停泊轨道的Walker星座部署方法 - Google Patents

一种基于双停泊轨道的Walker星座部署方法 Download PDF

Info

Publication number
CN106802667B
CN106802667B CN201611263948.1A CN201611263948A CN106802667B CN 106802667 B CN106802667 B CN 106802667B CN 201611263948 A CN201611263948 A CN 201611263948A CN 106802667 B CN106802667 B CN 106802667B
Authority
CN
China
Prior art keywords
orbit
satellite
working
parking
deployment
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
Application number
CN201611263948.1A
Other languages
English (en)
Other versions
CN106802667A (zh
Inventor
李成
魏世隆
何民
郭碧波
廖祥
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen Aerospace Dongfanghong Satellite Co.,Ltd.
Original Assignee
Aerospace Dongfanghong Development Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Aerospace Dongfanghong Development Ltd filed Critical Aerospace Dongfanghong Development Ltd
Priority to CN201611263948.1A priority Critical patent/CN106802667B/zh
Publication of CN106802667A publication Critical patent/CN106802667A/zh
Application granted granted Critical
Publication of CN106802667B publication Critical patent/CN106802667B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
    • G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
    • G05D1/10—Simultaneous control of position or course in three dimensions
    • G05D1/101—Simultaneous control of position or course in three dimensions specially adapted for aircraft
    • G05D1/104—Simultaneous control of position or course in three dimensions specially adapted for aircraft involving a plurality of aircrafts, e.g. formation flying

Landscapes

  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明提供一种基于双停泊轨道的Walker星座部署方法，包括以下步骤：步骤S1，将用于部署Walker星座的卫星分别依次送入第一停泊轨道和第二停泊轨道；步骤S2，第一颗卫星通过轨道机动进入工作轨道，其余卫星停留在第一停泊轨道和第二停泊轨道；步骤S3，等待当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，下一颗卫星通过轨道机动进入工作轨道；依此直到用于部署Walker星座的所有卫星依次进入工作轨道。本发明可以解决一箭多星部署单轨道面的局限，同时解决一箭多星单停泊轨道异面部署周期较长的缺点；在此基础上，还优化了星座部署时间和卫星相对位置的部署精度角度。

Description

一种基于双停泊轨道的Walker星座部署方法
技术领域
本发明涉及一种星座部署方法，尤其涉及一种基于双停泊轨道的Walker星座部署方法。
背景技术
轨道通信系统(Orbcomm系统)属于专用数据通信卫星系统，是美国轨道通信公司(Orbcomm LLC)于世纪90年开始建设运营的短数据通信星座系统，Orbcomm系统卫星星座由47颗卫星组成，包括7个轨道面，采用轨道高度低于1000km的LEO轨道，Orbcomm星座部署采用一箭多星部署一个轨道面，7个轨道面需要发射7枚火箭。
铱星通信星座、全球星通信星座、GPS导航卫星星座以及北斗导航卫星星座等卫星星座系统均是利用一箭多星部署一个轨道面的卫星，在卫星质量轻的情况下，如微小卫星，这种部署方法需采用多枚火箭进行异面部署，系统发射部署成本高。
COSMIC星座为美国和台湾地区联合研制的掩星大气探测星座，由6颗卫星组成，分布于6个轨道面，即每个轨道部署1颗卫星；该COSMIC星座于2006年一箭6星发射部署，采用一箭多星发射单停泊轨道异面部署的方法，6颗卫星首先进入500km左右的低停泊轨道，此后每颗卫星到达预定轨道面时，依靠自身轨道机动能力进入轨道高度为800km左右的任务轨道。
以上说的这些现有卫星星座系统或相关专利文献提出的星座部署方法，存在以下几项缺点中的一点或几点：第一、一箭多星只部署一个轨道面，多个轨道面需多次发射，发射部署准备周期长、成本高，对于单轨道面数量少、卫星质量轻的星座，发射部署效费比低；如Orbcomm星座、铱星星座、全球星星座、GPS导航卫星星座和北斗卫星星座采用这种发射部署方式，但上述星座一个轨道面卫星数量较多或卫星质量重，一箭部署一个轨道面有其合理性。第二、一箭多星发射，利用单停泊轨道实现异面部署，可以解决一箭多星只部署一个轨道所带来的局限性，COSMIC星座部署采用这种方法；但在卫星数量较多且轨道面分布范围较大的情况下，利用一箭多星单停泊轨道实现异面部署，周期较长。
发明内容
本发明所要解决的技术问题是需要提供一种能够解决一箭多星部署单轨道面的局限，并同时解决一箭多星单停泊轨道异面部署周期较长的缺点的Walker星座部署方法。
对此，本发明提供一种基于双停泊轨道的Walker星座部署方法，包括以下步骤：
步骤S1，将用于部署Walker星座的卫星分别依次送入第一停泊轨道和第二停泊轨道；
步骤S2，第一颗卫星通过轨道机动进入工作轨道，其余卫星停留在第一停泊轨道和第二停泊轨道；
步骤S3，等待当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，下一颗卫星通过轨道机动进入工作轨道；依此直到用于部署Walker星座的所有卫星依次进入工作轨道。
本发明的进一步改进在于，所述步骤S2中，第一颗卫星通过轨道机动自第一停泊轨道进入工作轨道；
所述步骤S3包括以下子步骤：
步骤S301，等待第二停泊轨道的轨道面相对于相邻的第一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，第二颗卫星通过轨道机动自第二停泊轨道进入工作轨道；
步骤S302，等待第一停泊轨道的轨道面相对于第一停泊轨道中在先一颗卫星所在工作轨道的轨道面进动，直到其升交点赤经相差预设角度时，第三颗卫星通过轨道机动自第一停泊轨道进入工作轨道；等待第二停泊轨道的轨道面相对于第二停泊轨道中在先一颗卫星所在工作轨道的轨道面进动，直到其升交点赤经相差预设角度时，第四颗卫星通过轨道机动自第二停泊轨道进入工作轨道；依此循环运行所述步骤S302，直到用于部署Walker星座的所有卫星依次进入工作轨道。
本发明的进一步改进在于，所述第一停泊轨道的轨道高度高于工作轨道的轨道高度，所述第二停泊轨道的轨道高度低于工作轨道的轨道高度。
本发明的进一步改进在于，所述预设角度为360°/用于部署Walker星座的卫星数量。
本发明的进一步改进在于，不同高度的当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面的进动速率差计算公式为
Figure BDA0001200308600000031
其中，a0为卫星工作轨道标称半长轴，e0为卫星工作轨道标称偏心率，i0为卫星工作轨道标称倾角，J2为地球的非球形引力，μ为地球引力常数，Re为地球参考半径，a为当前的停泊轨道半长轴；并根据轨道面升交点经度差及部署时间需求公式
Figure BDA0001200308600000032
求得当前的停泊轨道半长轴a，其中，ΔΩ为当前的停泊轨道的轨道面与工作轨道的轨道面之间的升交点赤经，t为卫星在停泊轨道上的驻留时间。
本发明的进一步改进在于，用于部署Walker星座的卫星数量为六颗，所述步骤S1中三颗卫星送入第一停泊轨道，另外三个卫星送入第二停泊轨道；所述预设角度为60°。
本发明的进一步改进在于，设
Figure BDA0001200308600000033
为工作轨道的升交点赤经进动速率，
Figure BDA0001200308600000034
为第二停泊轨道的升交点赤经进动速率，
Figure BDA0001200308600000035
为第一停泊轨道的升交点赤经进动速率，则
Figure BDA0001200308600000036
其中，
Figure BDA0001200308600000037
和
Figure BDA0001200308600000038
这三个进动速率可以通过公式
Figure BDA0001200308600000039
得到。
本发明的进一步改进在于，用于部署Walker星座的六颗卫星分别为A卫星、B卫星、C卫星、D卫星、E卫星和F卫星，所述步骤S1中，所述A卫星、B卫星和C卫星分别依次送入第二停泊轨道，所述D卫星、E卫星和F卫星分别依次送入第一停泊轨道；则所述Walker星座的部署过程如下：
步骤A，自入轨tD天后，D卫星首先通过降轨进入工作轨道，tD为D卫星部署完成所需时间；
步骤B，自入轨tA天后，第二停泊轨道的轨道面相对D卫星工作轨道的轨道面升交点赤经差异达到60°，A卫星升轨至工作轨道，tA为A卫星部署完成所需时间；
步骤C，自D卫星降轨tE天后，第一停泊轨道的轨道面相对D卫星工作轨道面升交点赤经差异达到60°，E卫星降轨至工作轨道，tE为E卫星部署完成所需时间；
步骤D，自A卫星升轨tB天后，第二停泊轨道的轨道面相对A卫星工作轨道面升交点赤经差异达到60°，B卫星升轨至工作轨道，tB为B卫星部署完成所需时间；
步骤E，自E卫星降轨tF天后，第一停泊轨道的轨道面相对E卫星工作轨道面升交点赤经差异达到60°，F卫星降轨至工作轨道，tF为F卫星部署完成所需时间；
步骤F，自B卫星升轨tC天后，第二停泊轨道的轨道面相对B卫星工作轨道面升交点赤经差异达到60°，C卫星升轨至工作轨道，tC为C卫星部署完成所需时间。
本发明的进一步改进在于，A卫星部署完成所需时间为
Figure BDA0001200308600000041
D卫星部署完成所需时间为
Figure BDA0001200308600000042
其中，Δt为预先设置好的隔离时间，该隔离时间用于避免出现两颗卫星在同一时期进行轨道机动。
本发明的进一步改进在于，E卫星部署完成所需时间公式为
Figure BDA0001200308600000043
其中，
Figure BDA0001200308600000044
为E卫星所在停泊轨道升交点赤经进动速率；B卫星部署完成所需时间公式为
Figure BDA0001200308600000045
其中，
Figure BDA0001200308600000046
为B卫星所在停泊轨道升交点赤经进动速率；F卫星部署完成所需时间公式为
Figure BDA0001200308600000051
其中，
Figure BDA0001200308600000052
为F卫星所在停泊轨道升交点赤经进动速率；C卫星部署完成所需时间公式为
Figure BDA0001200308600000053
其中，
Figure BDA0001200308600000054
为C卫星所在停泊轨道升交点赤经进动速率。E卫星所在停泊轨道、B卫星所在停泊轨道以及C卫星所在停泊轨道指的就是第一停泊轨道和第二停泊轨道。
与现有技术相比，本发明的有益效果在于：采用异面卫星星座部署方法实现一箭多星发射，并利用第一停泊轨道和第二停泊轨道的双停泊轨道实现异面部署，可以解决一箭多星部署单轨道面的局限，同时解决一箭多星单停泊轨道异面部署周期较长的缺点；在此基础上，本发明优化了星座部署时间和卫星与卫星之间相对位置的部署精度角度，提出将卫星从第一停泊轨道和第二停泊轨道依次送至工作轨道的机动策略。
附图说明
图1是本发明一种实施例的工作流程结构示意图；
图2是本发明一种实施例的星座部署原理示意图；
图3是本发明一种实施例的星座部署时序原理示意图。
具体实施方式
下面结合附图，对本发明的较优的实施例作进一步的详细说明：
如图1所示，本例提供一种基于双停泊轨道的Walker星座部署方法，包括以下步骤：
步骤S1，将用于部署Walker星座的卫星分别依次送入第一停泊轨道和第二停泊轨道；
步骤S2，第一颗卫星通过轨道机动进入工作轨道，其余卫星停留在第一停泊轨道和第二停泊轨道；
步骤S3，等待当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，下一颗卫星通过轨道机动进入工作轨道；依此直到用于部署Walker星座的所有卫星依次进入工作轨道。
本例所述第一停泊轨道的轨道高度优选高于工作轨道的轨道高度，所述第二停泊轨道的轨道高度优选低于工作轨道的轨道高度。所述预设角度一般设置为360°/用于部署Walker星座的卫星数量。
所述第一停泊轨道和第二停泊轨道统称为停泊轨道，即停泊轨道包括第一停泊轨道和第二停泊轨道；停泊轨道的用途类似于辅助车道，以此为跳板最终实现卫星进入设计工作轨道。本例所述双停泊轨道包括一个高于工作轨道的第一停泊轨道，另一个低于工作轨道的第二停泊轨道，这样设计的好处是使得“运行在高于工作轨道的第一停泊轨道上的卫星”与“运行在低于工作轨道的第二停泊轨道上的卫星”轨道升交点朝着相反方向进动，以此来实现轨道面之间的设计相位。
本例所述的高于工作轨道和低于工作轨道的限制范围相对宽松，其主要受到以下几个方面的约束：一、卫星自身速度增量的约束，变轨消耗的燃料不能过多，要保证轨道进动；二、部署时间约束，低地球轨道卫星设计寿命一般较短，过长的部署时间将增加星座构建风险以及星座系统的效能和实用性；三、轨道摄动的约束，第二停泊轨道高度不能过低，否则大气阻力影响多大，卫星需消耗多余的燃料用于维持轨道高度，对于需要长时间运行于停泊轨道的卫星尤其如此。此外，摄动力引起的轨道倾角差异和星间相对相位差异等也是需要考虑的问题。
本例的设计思路是：利用地球J2项非球形引力引起的轨道面进动实现星座的异面部署；对于地球卫星轨道，J2项非球形引力造成的升交点赤经进动速率为
Figure BDA0001200308600000061
该公式中，
Figure BDA0001200308600000062
为卫星所在轨道面升交点赤经进动速率，单位为rad/s；J2为地球J2项的非球形引力，无量纲，常取：J2＝0.00108263；μ为地球引力常数，单位为km3/s2，常取μ＝398600.5km3/s2；a为卫星轨道半长轴，单位为km；Re为地球参考半径，单位为km，常取：Re＝6378.14km；e为卫星轨道偏心率，无量纲；i为卫星轨道倾角，单位为deg。
由上式可知，在轨道偏心率e、轨道倾角i一定的情况下，轨道面升交点赤经进动速率随半长轴的变化而变化，即不同卫星之间轨道高度的差异将引起卫星轨道面升交点赤经漂移速率的差异，进而实现星座的异面部署。
如图2所示，所述Walker星座基本部署过程如下：(1)通过运载火箭等将卫星依次送入轨道高度低于或高于工作轨道的第二停泊轨道或第一停泊轨道；(2)第一颗卫星首先通过轨道机动进入工作轨道，其余卫星停留在停泊轨道；(3)由于轨道面进动速率差异，停泊轨道面相对首颗卫星的工作轨道面进动；(4)当停泊轨道面进动至与相邻卫星的工作轨道面升交点赤经相差约60°时，第二颗卫星通过轨道机动进入工作轨道，其余卫星通过这种方式依次进入预期工作轨道。
也就是说，所述步骤S2中，第一颗卫星通过轨道机动自第一停泊轨道进入工作轨道；
所述步骤S3包括以下子步骤：
步骤S301，等待第二停泊轨道的轨道面相对于相邻的第一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，第二颗卫星通过轨道机动自第二停泊轨道进入工作轨道；
步骤S302，等待第一停泊轨道的轨道面相对于第一停泊轨道中在先一颗卫星所在工作轨道的轨道面进动，直到其升交点赤经相差预设角度时，第三颗卫星通过轨道机动自第一停泊轨道进入工作轨道；等待第二停泊轨道的轨道面相对于第二停泊轨道中在先一颗卫星所在工作轨道的轨道面进动，直到其升交点赤经相差预设角度时，第四颗卫星通过轨道机动自第二停泊轨道进入工作轨道；依此循环运行所述步骤S302，直到用于部署Walker星座的所有卫星依次进入工作轨道。
根据近地轨道进动规律，不同高度的当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面的进动速率差计算公式为
Figure BDA0001200308600000071
其中，a0为卫星工作轨道标称半长轴，e0为卫星工作轨道标称偏心率，i0为卫星工作轨道标称倾角，J2为地球的非球形引力，μ为地球引力常数，Re为地球参考半径，a为当前的停泊轨道半长轴；并根据轨道面升交点经度差及部署时间需求公式
Figure BDA0001200308600000072
求得当前的停泊轨道半长轴a，其中，ΔΩ为当前的停泊轨道的轨道面与工作轨道的轨道面之间的相对升交点赤经，t为卫星在停泊轨道上的驻留时间。
本例以用于部署Walker星座的卫星数量为六颗为优选实施例，所述步骤S1中三颗卫星送入第一停泊轨道，另外三个卫星送入第二停泊轨道；所述预设角度为60°，由于设计星座相邻轨道面之间的升交点赤经间隔为60°，停泊轨道(包括第一停泊轨道和第二停泊轨道)的轨道面进动至与相邻工作轨道的轨道面升交点赤经相差约60°时意味着运行在当前的停泊轨道(包括第一停泊轨道和第二停泊轨道)上的卫星完成了轨道调整。此时，当前的停泊轨道(包括第一停泊轨道和第二停泊轨道)上的卫星需变轨进入工作轨道。
在轨道偏心率e和轨道倾角i一定的情况下，轨道面升交点进动速率随半长轴的变化而变化，即不同卫星之间轨道高度的差异将引起卫星轨道面升交点进动速率的差异，可以此来实现星座的异面部署。假设
Figure BDA0001200308600000081
为工作轨道的升交点赤经进动速率，单位为rad/s；
Figure BDA0001200308600000082
为第二停泊轨道的升交点赤经进动速率，单位为rad/s；
Figure BDA0001200308600000083
为第一停泊轨道的升交点赤经进动速率，单位为rad/s；则
Figure BDA0001200308600000084
因而，随着时间的积累，运行在第二停泊轨道上的卫星就与运行在第一停泊轨道上的卫星产生了明显的升交点赤经差异，当这个差异达到60°后部署一颗位于第一停泊轨道或第二停泊轨道上的卫星进入工作轨道；当这个差异达到120°时部署下一颗位于第二停泊轨道或第一停泊轨道上的卫星进入工作轨道；当这个差异达到180°时部署下一颗位于第一停泊轨道或第二停泊轨道上的卫星进入工作轨道；当这个差异达到240°时部署下一颗位于第二停泊轨道或第一停泊轨道上的卫星进入工作轨道；当这个差异达到300°时部署下一颗位于第一停泊轨道或第二停泊轨道上的卫星进入工作轨道。由于第一颗卫星作为基准星进入工作轨道，故升交点赤经差异达到300°后星座即可部署完成。也就是说，等待当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，下一颗卫星通过轨道机动自另一个停泊轨道进入工作轨道；依此直到用于部署Walker星座的所有卫星依次进入工作轨道。
本例所述Walker星座优选共6颗卫星，根据任务要求，需将每颗卫星部署到6个不同的轨道面上。这6个轨道面均匀分布在360°升交点赤经上，相邻两个轨道面之间的升交点赤经间隔为60°。由于6颗卫星采用一箭6星的方式发射入轨，入轨后6星处于升交点赤经相近的轨道上。
本例用于部署Walker星座的六颗卫星分别为A卫星、B卫星、C卫星、D卫星、E卫星和F卫星，所述步骤S1中，所述A卫星、B卫星和C卫星分别依次送入第二停泊轨道，所述D卫星、E卫星和F卫星分别依次送入第一停泊轨道；如图3所示，则所述Walker星座的部署过程如下：
步骤A，自入轨tD天后，D卫星首先通过降轨进入工作轨道，tD为D卫星部署完成所需时间，其余卫星停留在各自的停泊轨道(包括第一停泊轨道和第二停泊轨道)；
步骤B，自入轨tA天后，第二停泊轨道的轨道面相对D卫星工作轨道的轨道面升交点赤经差异达到60°，A卫星升轨至工作轨道，tA为A卫星部署完成所需时间，届时A卫星和D卫星之间的轨道面升交点赤经相差60°，A卫星相位滞后D卫星240°；
步骤C，自D卫星降轨tE天后，第一停泊轨道的轨道面相对D卫星工作轨道面升交点赤经差异达到60°，E卫星降轨至工作轨道，tE为E卫星部署完成所需时间，届时D卫星和E卫星之间的轨道面升交点赤经相差60°，E卫星相位超前D卫星240°；
步骤D，自A卫星升轨tB天后，第二停泊轨道的轨道面相对A卫星工作轨道面升交点赤经差异达到60°，B卫星升轨至工作轨道，tB为B卫星部署完成所需时间，届时B卫星与A卫星之间的轨道面升交点赤经相差60°，B卫星相位滞后A卫星240°；
步骤E，自E卫星降轨tF天后，第一停泊轨道的轨道面相对E卫星工作轨道面升交点赤经差异达到60°，F卫星降轨至工作轨道，tF为F卫星部署完成所需时间，届时F卫星与E卫星之间的轨道面升交点赤经相差60°，F卫星相位超前E卫星240°；
步骤F，自B卫星升轨tC天后，第二停泊轨道的轨道面相对B卫星工作轨道面升交点赤经差异达到60°，C卫星升轨至工作轨道，tC为C卫星部署完成所需时间，届时C卫星与B卫星之间的轨道面升交点赤经相差60°，C卫星相位滞后B卫星240°。
经过以上描述，如图3所示，6颗卫星部署序列为D卫星→A卫星→E卫星→B卫星→F卫星→C卫星，轨道面自西向东排序为C卫星→B卫星→A卫星→D卫星→E卫星→F卫星，整个星座部署时序示意图如下图3所示。
本例所述Walker星座的星座部署时序示意图如图3所示，所述A卫星部署完成所需时间为
Figure BDA0001200308600000101
D卫星部署完成所需时间为
Figure BDA0001200308600000102
其中，Δt为预先设置好的隔离时间，该隔离时间用于避免出现两颗卫星在同一时期进行轨道机动。
所述E卫星部署完成所需时间为
Figure BDA0001200308600000103
其中，
Figure BDA0001200308600000104
为E卫星所在停泊轨道升交点赤经进动速率；B卫星部署完成所需时间为
Figure BDA0001200308600000105
其中，
Figure BDA0001200308600000106
为B卫星所在停泊轨道升交点赤经进动速率；F卫星部署完成所需时间为
Figure BDA0001200308600000107
其中，
Figure BDA0001200308600000108
为F卫星所在停泊轨道升交点赤经进动速率；C卫星部署完成所需时间为
Figure BDA0001200308600000109
其中，
Figure BDA00012003086000001010
为C卫星所在停泊轨道升交点赤经进动速率。
若不考虑初始入轨时间差异及升降轨所需时间，根据图3所示的星座部署时序示意图，第二停泊轨道的全部卫星完成升轨部署所需时间为tL＝tA+tB+tC；第一停泊轨道的全部卫星完成降轨部署所需时间为tH＝tD+tE+tF+Δt。卫星所在停泊轨道运行时间与相对升交点赤经存在以下关系式
Figure BDA00012003086000001011
该公式中，ΔΩAD为初始时刻A卫星相对D卫星的升交点赤经之间的差值。为了使部署时间T取得最小值的tA，其为
Figure BDA00012003086000001012
将tA代入前式求得tD为
Figure BDA0001200308600000111
在各卫星的停泊轨道参数为标称值的情况下，tA求解式简化为
Figure BDA0001200308600000112
tD求解式简化为
Figure BDA0001200308600000113
根据霍曼变轨规律，从停泊轨道至工作轨道可以通过两次脉冲速度增量实现，单次速度增量及总速度增量分别为
Figure BDA0001200308600000114
该公式中，ΔV1为霍曼转移第一次脉冲所需速度增量，单位为km/s；ΔV2为霍曼转移第二次脉冲所需速度增量，单位为km/s；ΔV为霍曼转移所需速度增量总和，单位为km/s；a为停泊轨道半长轴，单位为km；a0为工作轨道半长轴，单位为km。
基准星部署轨道机动规划以星座总部署时间最优为目标，同时综合考虑部署机动时间、推进和测控等条件约束，权衡确定基准星降轨时机及降轨高度。星座部署包括基准星和非基准星的轨道机动部署。基准星的部署以星座部署所需时间为优化目标，基准星最佳降轨时机为
Figure BDA0001200308600000115
考虑工程实施约束，基准星的降轨过程尽量保证在测控圈内，离tDopt最近的时间点作为次优的升降轨时间点。非基准星部署轨道机动规划以相对相位最优且相对升交点经度次优为目标，确定部署机动日期及变轨量。升交点赤经偏差最小的首次部署机动日期为
Figure BDA0001200308600000121
上式中，floor()为向下取整函数，如floor(5.6)＝5，floor(5.3)＝5；t0为设定起始时刻；tm-t0之后第N天卫星实施首次部署机动，满足测控可见性及变轨点纬度幅角要求的时刻；Δu0为t0初始时刻相对目标部署相位的偏差，单位为rad；
Figure BDA0001200308600000122
为非基准星相对基准星的轨道角速率，单位为rad/s。Δumr为升降轨期间累积的相对相位，其符号与k相同。k为整数，算式如下：
Figure BDA0001200308600000123
ΔΩ0为t0初始时刻非基准星相对目标升交点赤经的偏差。p的算式如下：
Figure BDA0001200308600000124
该公式中，J2为地球非球形引力J2带谐项系数。Re为地球参考半径。ao为基准星工作轨道半长轴，单位为km。io为基准星工作轨道倾角，单位为deg。c＝ap/ao，ap为非基准星停泊轨道半长轴，无量纲。
在确定首次部署机动日期的情况下，根据测控及纬度幅角条件选定tm由于实际变轨过程总是存在误差，当天的变轨量需对之前的累积相位偏差进行修正，以保证后续相位部署精度。设第d天机动之前每天相对基准星的实际累积相对相位为Δλi(i＝1...d-1,d＜6)，重新规划第d天相对基准星需累积的相位为
Figure BDA0001200308600000125
单位为rad/day，其中
Figure BDA0001200308600000126
第d天轨道半长轴机动目标量为
Figure BDA0001200308600000127
该公式中，tin为当天霍曼转移第一次变轨时刻，根据卫星测控可见性具体选定；tout为后一天霍曼转移第一次变轨时刻，根据卫星测控可见性具体选定；这两个时刻根据卫星测控可见性具体选定。no为基准星轨道角速率，rad-1，按下式计算：
Figure BDA0001200308600000131
本例采用异面卫星星座部署方法实现一箭多星发射，并利用第一停泊轨道和第二停泊轨道的双停泊轨道实现异面部署，可以解决一箭多星部署单轨道面的局限，同时解决一箭多星单停泊轨道异面部署周期较长的缺点；在此基础上，本例优化了星座部署时间和卫星与卫星之间相对位置的部署精度角度，提出将卫星从第一停泊轨道和第二停泊轨道依次送至工作轨道的机动策略。
以上内容是结合具体的优选实施方式对本发明所作的进一步详细说明，不能认定本发明的具体实施只局限于这些说明。对于本发明所属技术领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干简单推演或替换，都应当视为属于本发明的保护范围。

Claims (9)

1.一种基于双停泊轨道的Walker星座部署方法，其特征在于，包括以下步骤：
步骤S1，将用于部署Walker星座的卫星分别依次送入第一停泊轨道和第二停泊轨道；
步骤S2，第一颗卫星通过轨道机动进入工作轨道，其余卫星停留在第一停泊轨道和第二停泊轨道；
步骤S3，等待当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，下一颗卫星通过轨道机动进入工作轨道；依此直到用于部署Walker星座的所有卫星依次进入工作轨道；所述步骤S2中，第一颗卫星通过轨道机动自第一停泊轨道进入工作轨道；
所述步骤S3包括以下子步骤：
步骤S301，等待第二停泊轨道的轨道面相对于相邻的第一颗卫星所在工作轨道的轨道面进动直到其升交点赤经相差预设角度时，第二颗卫星通过轨道机动自第二停泊轨道进入工作轨道；
步骤S302，等待第一停泊轨道的轨道面相对于第一停泊轨道中在先一颗卫星所在工作轨道的轨道面进动，直到其升交点赤经相差预设角度时，第三颗卫星通过轨道机动自第一停泊轨道进入工作轨道；等待第二停泊轨道的轨道面相对于第二停泊轨道中在先一颗卫星所在工作轨道的轨道面进动，直到其升交点赤经相差预设角度时，第四颗卫星通过轨道机动自第二停泊轨道进入工作轨道；依此循环运行所述步骤S302，直到用于部署Walker星座的所有卫星依次进入工作轨道。
2.根据权利要求1所述的基于双停泊轨道的Walker星座部署方法，其特征在于，所述第一停泊轨道的轨道高度高于工作轨道的轨道高度，所述第二停泊轨道的轨道高度低于工作轨道的轨道高度。
3.根据权利要求1所述的基于双停泊轨道的Walker星座部署方法，其特征在于，所述预设角度为360°/用于部署Walker星座的卫星数量。
4.根据权利要求1至3任意一项所述的基于双停泊轨道的Walker星座部署方法，其特征在于，不同高度的当前停泊轨道的轨道面相对于相邻一颗卫星所在工作轨道的轨道面的进动速率差计算公式为
Figure FDA0002249520990000021
其中，a0为卫星工作轨道标称半长轴，e0为卫星工作轨道标称偏心率，i0为卫星工作轨道标称倾角，J2为地球的非球形引力，μ为地球引力常数，Re为地球参考半径，a为当前的停泊轨道半长轴；并根据轨道面升交点经度差及部署时间需求公式
Figure FDA0002249520990000022
求得当前的停泊轨道半长轴a，其中，ΔΩ为当前的停泊轨道的轨道面与工作轨道的轨道面之间的升交点赤经，t为卫星在停泊轨道上的驻留时间。
5.根据权利要求1至3任意一项所述的基于双停泊轨道的Walker星座部署方法，其特征在于，用于部署Walker星座的卫星数量为六颗，所述步骤S1中三颗卫星送入第一停泊轨道，另外三个卫星送入第二停泊轨道；所述预设角度为60°。
6.根据权利要求5所述的基于双停泊轨道的Walker星座部署方法，其特征在于，设
Figure FDA0002249520990000023
为工作轨道的升交点赤经进动速率，
Figure FDA0002249520990000024
为第二停泊轨道的升交点赤经进动速率，
Figure FDA0002249520990000025
为第一停泊轨道的升交点赤经进动速率，则
Figure FDA0002249520990000026
7.根据权利要求6所述的基于双停泊轨道的Walker星座部署方法，其特征在于，用于部署Walker星座的六颗卫星分别为A卫星、B卫星、C卫星、D卫星、E卫星和F卫星，所述步骤S1中，所述A卫星、B卫星和C卫星分别依次送入第二停泊轨道，所述D卫星、E卫星和F卫星分别依次送入第一停泊轨道；则所述Walker星座的部署过程如下：
步骤A，自入轨tD天后，D卫星首先通过降轨进入工作轨道，tD为D卫星部署完成所需时间；
步骤B，自入轨tA天后，第二停泊轨道的轨道面相对D卫星工作轨道的轨道面升交点赤经差异达到60°，A卫星升轨至工作轨道，tA为A卫星部署完成所需时间；
步骤C，自D卫星降轨tE天后，第一停泊轨道的轨道面相对D卫星工作轨道面升交点赤经差异达到60°，E卫星降轨至工作轨道，tE为E卫星部署完成所需时间；
步骤D，自A卫星升轨tB天后，第二停泊轨道的轨道面相对A卫星工作轨道面升交点赤经差异达到60°，B卫星升轨至工作轨道，tB为B卫星部署完成所需时间；
步骤E，自E卫星降轨tF天后，第一停泊轨道的轨道面相对E卫星工作轨道面升交点赤经差异达到60°，F卫星降轨至工作轨道，tF为F卫星部署完成所需时间；
步骤F，自B卫星升轨tC天后，第二停泊轨道的轨道面相对B卫星工作轨道面升交点赤经差异达到60°，C卫星升轨至工作轨道，tC为C卫星部署完成所需时间。
8.根据权利要求7所述的基于双停泊轨道的Walker星座部署方法，其特征在于，A卫星部署完成所需时间为
Figure FDA0002249520990000031
D卫星部署完成所需时间为
Figure FDA0002249520990000032
其中，Δt为预先设置好的隔离时间，该隔离时间用于避免出现两颗卫星在同一时期进行轨道机动。
9.根据权利要求7所述的基于双停泊轨道的Walker星座部署方法，其特征在于，E卫星部署完成所需时间公式为
Figure FDA0002249520990000033
其中，
Figure FDA0002249520990000034
为E卫星所在停泊轨道升交点赤经进动速率；B卫星部署完成所需时间公式为
Figure FDA0002249520990000035
其中，
Figure FDA0002249520990000036
为B卫星所在停泊轨道升交点赤经进动速率；F卫星部署完成所需时间公式为
Figure FDA0002249520990000037
其中，
Figure FDA0002249520990000038
为F卫星所在停泊轨道升交点赤经进动速率；C卫星部署完成所需时间公式为
Figure FDA0002249520990000041
其中，
Figure FDA0002249520990000042
为C卫星所在停泊轨道升交点赤经进动速率。
CN201611263948.1A 2016-12-30 2016-12-30 一种基于双停泊轨道的Walker星座部署方法 Active CN106802667B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611263948.1A CN106802667B (zh) 2016-12-30 2016-12-30 一种基于双停泊轨道的Walker星座部署方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611263948.1A CN106802667B (zh) 2016-12-30 2016-12-30 一种基于双停泊轨道的Walker星座部署方法

Publications (2)

Publication Number Publication Date
CN106802667A CN106802667A (zh) 2017-06-06
CN106802667B true CN106802667B (zh) 2020-03-10

Family

ID=58985347

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611263948.1A Active CN106802667B (zh) 2016-12-30 2016-12-30 一种基于双停泊轨道的Walker星座部署方法

Country Status (1)

Country Link
CN (1) CN106802667B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108614575A (zh) * 2018-06-20 2018-10-02 北京电子工程总体研究所 一种卫星静止轨道定点位置调整方法
CN109765926B (zh) * 2018-12-27 2021-12-07 宁波天擎航天科技有限公司 一种升交点赤经精确控制方法、系统及装置
CN110198184B (zh) * 2019-05-22 2021-09-07 航天科工空间工程发展有限公司 一种低轨星座系统间频谱共存星座设计方法
CN111177948A (zh) * 2020-01-15 2020-05-19 东方红卫星移动通信有限公司 一种利用轨道要素进行描述的典型Walker星座数学模型
CN113525719B (zh) * 2021-07-19 2023-03-14 中国人民解放军战略支援部队航天工程大学 一种大规模低轨卫星星座的部署方法、装置、设备及存储介质
CN114006646B (zh) * 2021-09-27 2023-09-29 中国人民解放军战略支援部队航天工程大学 一种Walker星座构型维持的轨道控制频次分析方法及装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102197605A (zh) * 2008-08-28 2011-09-21 法国国家太空研究中心 用于接收并存储卫星数据的地面站网格
CN104965982A (zh) * 2015-06-17 2015-10-07 中国人民解放军63920部队 一种地月平动点星座的布设方法
CN105335541A (zh) * 2014-08-12 2016-02-17 中国人民解放军装备学院 导航卫星星座的工程设计方法
CN105353384A (zh) * 2014-08-18 2016-02-24 中国人民解放军装备学院 一种导航星座备份方法
CN106156417A (zh) * 2016-07-04 2016-11-23 中国人民解放军国防科学技术大学 等时间间隔快速重访卫星星座构型优化设计方法
CN106249253A (zh) * 2016-07-20 2016-12-21 中国人民解放军装备学院 低轨通信和导航增强混合星座的优化设计方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102197605A (zh) * 2008-08-28 2011-09-21 法国国家太空研究中心 用于接收并存储卫星数据的地面站网格
CN105335541A (zh) * 2014-08-12 2016-02-17 中国人民解放军装备学院 导航卫星星座的工程设计方法
CN105353384A (zh) * 2014-08-18 2016-02-24 中国人民解放军装备学院 一种导航星座备份方法
CN104965982A (zh) * 2015-06-17 2015-10-07 中国人民解放军63920部队 一种地月平动点星座的布设方法
CN106156417A (zh) * 2016-07-04 2016-11-23 中国人民解放军国防科学技术大学 等时间间隔快速重访卫星星座构型优化设计方法
CN106249253A (zh) * 2016-07-20 2016-12-21 中国人民解放军装备学院 低轨通信和导航增强混合星座的优化设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
卫星星座构形控制与设计研究;项军华;《工程科技II辑》;20090715;第C031-43页 *
星座与星座轨道控制技术;盘科炎等;《航天控制》;20020331;第51-57页 *
木卫停泊轨道间低耗能小推力转移轨道设计方法研究;路毅等;《力学与实践》;20161031;第38卷(第5期);第501-510页 *

Also Published As

Publication number Publication date
CN106802667A (zh) 2017-06-06

Similar Documents

Publication Publication Date Title
CN106802667B (zh) 一种基于双停泊轨道的Walker星座部署方法
Woffinden et al. Navigating the road to autonomous orbital rendezvous
WO2020158505A1 (ja) 衛星コンステレーション形成システム、衛星コンステレーション形成方法、デブリ除去方式、衛星コンステレーション構築方式、地上設備、宇宙交通管理システム、宇宙物体管理部、および軌道降下時衝突回避運用方法
CN102880184A (zh) 一种静止轨道卫星自主轨道控制方法
US10752383B1 (en) Systems and methods for satellite constellation launch using air-launched vehicles
CN111367167B (zh) 一种有限时间约束的高轨目标抵近观测方法
CN113525719B (zh) 一种大规模低轨卫星星座的部署方法、装置、设备及存储介质
WO2015057102A1 (ru) Космическая обслуживающая система и способ ее построения
CN112179207B (zh) 航天器快速交会的运载火箭任务规划方法
Finley et al. Techniques for leo constellation deployment and phasing utilizing differential aerodynamic drag
CN103885456A (zh) 一箭多星异轨面分布转移方法
Kominato et al. Optical hybrid navigation and station keeping around Itokawa
CN113255109B (zh) 针对高轨非零倾角卫星在轨健康监测轨道设计方法及系统
Bérend et al. ALTAIR innovative air-launch system-consolidated design, lessons learned and perspective
Akim et al. Mechanics and Motion Control of a Space Vehicle in the Project of Relict Substance Delivery on Earth (the Project" Phobos_Grunt")
Martin et al. Flight dynamics analysis for philae landing site selection
De Pasquale et al. ATV Jules Verne reentry observation: Mission design and trajectory analysis
Wermuth et al. Operational Concept of a Picosatellite Release from a LEO Satellite
Levi et al. Strategy for Mitigating Collision Between Landsat-5 and the Afternoon Constellation
Finger et al. State Spaceports-If You Build it, Will They Come?
RU2771191C1 (ru) Космическая обслуживающая система и способ ее функционирования
RU2774915C1 (ru) Способ обслуживания искусственных спутников на средних и высоких околоземных орбитах
Bonnal et al. Optimization of Tethered De-Orbitation of Spent Upper Stages: the``MAILMAN''Process
Ockels et al. EuroMoon 2000 a plan for a European lunar south Pole expedition
Clark Operations with Tiangong 2

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
CP03 Change of name, title or address

Address after: 518000 whole building of satellite building, 61 Gaoxin South Jiudao, Yuehai street, Nanshan District, Shenzhen City, Guangdong Province

Patentee after: Shenzhen Aerospace Dongfanghong Satellite Co.,Ltd.

Address before: 518057 1-6 / F, Weixing building, 2002 Keyuan Road, Nanshan District, Shenzhen City, Guangdong Province

Patentee before: AEROSPACE DONGFANGHONG DEVELOPMENT Ltd.

CP03 Change of name, title or address