异常分离情况下组合体GEO变轨策略生成方法及系统与流程

异常分离情况下组合体geo变轨策略生成方法及系统
技术领域
1.本发明涉及宇航飞行器轨道动力学，具体地，涉及一种异常分离情况下组合体geo变轨策略生成方法及系统。


背景技术：

2.地球静止轨道航天器(geo)具有对地球相对“静止”特点，被广泛用于通信、导航、中继等领域。相较于传统单体航天器而言，组合体航天器由于每个舱段都携带推进系统，一方面在应对主动段与多次轨道转移段异常情况时具有较好的鲁棒性，另一方面也就要求在组合体异常分离情况下需要快速对组合体或组合体每个舱段设计轨道转移策略。
3.国内外目前已经应用的转移轨道策略设计方案较多，但多是基于有限推力的复杂积分寻优或基于混合推进系统的分段变轨方法。刘俊尧、赵建伟、宗岩等在中国专利“一种脉冲轨道转移向有限推力轨道转移转换的方法”(专利文献：cn112455725a)中首先确定脉冲点火点位置，然后在此位置进行有限推力积分寻优，并对点火方向进行修正，对于单体航天器多次点火轨道转移问题，该方法计算量大，迭代寻优收敛速度慢，对于组合航天器多舱段该缺陷尤为显著。林宝军、蒋国伟、范苑等在中国专利“小型地球静止卫星轨道转移方法及系统”(专利文献：cn111891396a)中提出了一种首先利用化学推进系统进行航天器近地点高度提升，其次利用电推进系统进行航天器近地点高度、倾角、偏心率调整，最后利用化学推进系统进行航天器定点捕获，该方法一方面需要单体航天器同时具备化学与电推进两套系统，局限性较强；一方面将近地点高度、倾角与偏心率分阶段调整，燃料消耗较大，影响航天器在轨寿命。
4.综上所述，需要针对异常分离条件的组合体航天器geo轨道转移策略快速生成需求，开展变轨策略方法设计优化。


技术实现要素：

5.针对现有技术中的缺陷，本发明的目的是提供一种异常分离情况下组合体geo变轨策略生成方法及系统。
6.根据本发明提供的一种异常分离情况下组合体geo变轨策略生成方法，包括：
7.步骤a：根据组合体与运载器或组合体之间异常分离时刻状态，确定组合体分离点轨道参数及与组合体轨道转移策略生成相关的平台参数；
8.步骤b：根据平台参数确定轨道转移策略的设计约束条件；
9.步骤c：根据所述约束条件以轨道转移所需速度增量最小为原则，进行最优变轨次数估计；
10.步骤d：根据估计结果，建立描述单次变轨状态变化及多次变轨间联系的数学模型；
11.步骤e：利用所述数学模型以geo轨道为目标，通过迭代寻优生成满足设计约束的轨道转移策略。
12.优选地，所述步骤a包括：
13.步骤s2.1：组合体包含多个舱段，每个舱段都携带推进系统，异常分离包括运载器未能将组合体送至预定轨道即提前分离的发射段异常与组合体中某舱段未能将组合体送至准地球同步轨道即提前分离的转移段异常；首先明确异常分离时刻的轨道状态，采用geo轨道通用形式进行描述，半长轴a0、倾角i0、偏心率e0、地理经度λ0、地理纬度η0；
14.步骤s2.2：明确组合体航天器、主星或主星与推进舱的相关平台参数，包括分离时刻重量m0、可用燃料剩余量m'0、发动机推力f0、发动机比冲isp0。
15.优选地，所述步骤b中：
16.结合实际工程归纳出三方面设计约束：地面测控条件对点火段地理经度跨度及测控时长约束、推力器特性、弧段损失及安全保护需求对单次最长点火时间约束、变轨对象可用燃料约束。
17.优选地，所述步骤c中，为了确定轨道转移优化参数维数需要依照如下步骤以转移所需速度增量最小为原则，进行最优变轨次数估计：
18.步骤s4.1：依照如下公式计算得到航天器从异常分离时刻状态转移至geo目标轨道的最小速度增量δv：
[0019][0020]
δv＝v
n
‑
v0[0021]
其中，r0为转移轨道设计对象异常分离时刻的地心距，由其轨道根数计算得到；a0为轨道半长轴；μ为地球引力常数；v0为转移轨道设计对象异常分离时刻速度；v
n
为目标轨道速度；a
n
为目标轨道半长轴；
[0022]
步骤s4.2：根据火箭公式，计算得到轨道转移最小速度增量对应的燃料消耗δm：
[0023][0024]
其中，g＝9.80665m/s2为地球引力加速度；
[0025]
步骤s4.3：根据发动机比冲计算推进剂秒流量dm，并结合燃料消耗计算总共的变轨时长t：
[0026][0027][0028]
步骤s4.4：结合推力器单次最长点火时间t
max
计算得到最优变轨次数估计值n：
[0029][0030]
其中[x]为屋檐函数，表示取≥x的最小整数。
[0031]
优选地，所述步骤d包括：
[0032]
步骤s5.1：根据每次点火待优化变量：变轨目标半长轴a
k
、变轨前偏移圈次q
k
，k＝1,
…
,n，计算点火前参数：
[0033]
轨道周期t
k
：
[0034][0035]
经度漂移率
[0036][0037]
其中，ω
e
＝7.2921
×
10
‑5rad/s，π为圆周率。
[0038]
升交点地理经度λ
k
：
[0039][0040]
步骤s5.2：根据每次变轨目标半长轴a
k
，k＝1，
…
,n，计算得到单次点火点及目标轨道参数：
[0041]
由活力公式，计算目标轨道速度v
k
：
[0042][0043]
其中，r
k
为转移轨道设计对象点火时刻的地心距，由其轨道根数计算得到；
[0044]
变轨速度增量δv
k
：
[0045][0046]
其中，α与β
k
依照如下公式计算
[0047][0048]
β
k
＝π
‑
i
k
‑1‑
α
[0049]
变轨后倾角i
k
：
[0050][0051]
变轨燃料消耗δm
k
：
[0052][0053]
点火时长t
k
：
[0054][0055]
优选地，所述步骤e包括：
[0056]
步骤s6.1：将geo轨道转移策略求解问题转化为多变量、多目标、多约束寻优问题，优化模型描述如下：
[0057]
目标函数：
[0058][0059]
其中，表示轨道目标罚函数；表示第j个轨道参数变轨终值与目标轨道偏差；δorb＝(δa,δi,δe,δλ,δη)，表示变轨结果与目标轨道偏差；γ表示权重系数；
[0060]
约束条件：
[0061][0062]
步骤s6.2：针对上述优化模型选择多变量、多目标、多约束优化算法进行迭代寻优；
[0063]
步骤s6.3：输出寻优结果，确定异常分离条件的组合体航天器geo轨道转移策略。
[0064]
根据本发明提供的一种异常分离情况下组合体geo变轨策略生成系统，包括：
[0065]
模块a：根据组合体与运载器或组合体之间异常分离时刻状态，确定组合体分离点轨道参数及与组合体轨道转移策略生成相关的平台参数；
[0066]
模块b：根据平台参数确定轨道转移策略的设计约束条件；
[0067]
模块c：根据所述约束条件以轨道转移所需速度增量最小为原则，进行最优变轨次数估计；
[0068]
模块d：根据估计结果，建立描述单次变轨状态变化及多次变轨间联系的数学模型；
[0069]
模块e：利用所述数学模型以geo轨道为目标，通过迭代寻优生成满足设计约束的轨道转移策略。
[0070]
优选地，所述模块a包括：
[0071]
模块s2.1：组合体包含多个舱段，每个舱段都携带推进系统，异常分离包括运载器未能将组合体送至预定轨道即提前分离的发射段异常与组合体中某舱段未能将组合体送至准地球同步轨道即提前分离的转移段异常；首先明确异常分离时刻的轨道状态，采用geo轨道通用形式进行描述，半长轴a0、倾角i0、偏心率e0、地理经度λ0、地理纬度η0；
[0072]
模块s2.2：明确组合体航天器、主星或主星与推进舱的相关平台参数，包括分离时刻重量m0、可用燃料剩余量m'0、发动机推力f0、发动机比冲isp0。
[0073]
所述模块b中：
[0074]
结合实际工程归纳出三方面设计约束：地面测控条件对点火段地理经度跨度及测控时长约束、推力器特性、弧段损失及安全保护需求对单次最长点火时间约束、变轨对象可用燃料约束。
[0075]
优选地，所述模块c中，为了确定轨道转移优化参数维数需要依照如下模块以转移所需速度增量最小为原则，进行最优变轨次数估计：
[0076]
模块s4.1：依照如下公式计算得到航天器从异常分离时刻状态转移至geo目标轨道的最小速度增量δv：
[0077][0078]
δv＝v
n
‑
v0[0079]
其中，r0为转移轨道设计对象异常分离时刻的地心距，由其轨道根数计算得到；a0为轨道半长轴；μ为地球引力常数；v0为转移轨道设计对象异常分离时刻速度；v
n
为目标轨道速度；a
n
为目标轨道半长轴；
[0080]
模块s4.2：根据火箭公式，计算得到轨道转移最小速度增量对应的燃料消耗δm：
[0081][0082]
其中，g＝9.80665m/s2为地球引力加速度；
[0083]
模块s4.3：根据发动机比冲计算推进剂秒流量dm，并结合燃料消耗计算总共的变轨时长t：
[0084][0085][0086]
模块s4.4：结合推力器单次最长点火时间t
max
计算得到最优变轨次数估计值n：
[0087][0088]
其中[x]为屋檐函数，表示取≥x的最小整数。
[0089]
优选地，所述模块d包括：
[0090]
模块s5.1：根据每次点火待优化变量：变轨目标半长轴a
k
、变轨前偏移圈次q
k
，k＝1,
…
,n，计算点火前参数：
[0091]
轨道周期t
k
：
[0092][0093]
经度漂移率
[0094][0095]
其中，ω
e
＝7.2921
×
10
‑5rad/s，π为圆周率。
[0096]
升交点地理经度λ
k
：
[0097][0098]
模块s5.2：根据每次变轨目标半长轴a
k
，k＝1,
…
,n，计算得到单次点火点及目标轨道参数：
[0099]
由活力公式，计算目标轨道速度v
k
：
[0100][0101]
其中，r
k
为转移轨道设计对象点火时刻的地心距，由其轨道根数计算得到；
[0102]
变轨速度增量δv
k
：
[0103][0104]
其中，α与β
k
依照如下公式计算
[0105][0106]
β
k
＝π
‑
i
k
‑1‑
α
[0107]
变轨后倾角i
k
：
[0108][0109]
变轨燃料消耗δm
k
：
[0110][0111]
点火时长t
k
：
[0112][0113]
优选地，所述模块e包括：
[0114]
模块s6.1：将geo轨道转移策略求解问题转化为多变量、多目标、多约束寻优问题，优化模型描述如下：
[0115]
目标函数：
[0116][0117]
其中，表示轨道目标罚函数；表示第j个轨道参数变轨终值与目标轨道偏差；δorb＝(δa,δi,δe,δλ,δη)，表示变轨结果与目标轨道偏差；γ表示权重系数；
[0118]
约束条件：
[0119][0120]
模块s6.2：针对上述优化模型选择多变量、多目标、多约束优化算法进行迭代寻优；
[0121]
模块s6.3：输出寻优结果，确定异常分离条件的组合体航天器geo轨道转移策略。
[0122]
与现有技术相比，本发明具有如下的有益效果：
[0123]
本发明解决了组合体航天器异常分离条件下，传统有限推力复杂设计方法计算量大、计算速度慢等不足，具有一定的工程实用性。
附图说明
[0124]
通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、目的和优点将会变得更明显：
[0125]
图1异常分离情况下组合体geo变轨策略快速生成方法原理框图；
[0126]
图2组合体航天器(两舱段)构型示意图；
[0127]
图3组合体航天器与运载异常分离(倾角超差)轨道示意图；
[0128]
图4基于异常分离条件的组合体航天器geo轨道转移策略仿真图。
具体实施方式
[0129]
下面结合具体实施例对本发明进行详细说明。以下实施例将有助于本领域的技术人员进一步理解本发明，但不以任何形式限制本发明。应当指出的是，对本领域的普通技术人员来说，在不脱离本发明构思的前提下，还可以做出若干变化和改进。这些都属于本发明的保护范围。
[0130]
如附图1所示，本发明提供了一种异常分离情况下组合体geo变轨策略快速生成方法。具体包含如下步骤：
[0131]
步骤a：根据组合体与运载器或组合体之间异常分离时刻状态，确定组合体分离点轨道参数及与组合体轨道转移策略生成相关的平台参数。
[0132]
步骤b：根据地面测控需求、单次点火能力等因素，确定轨道转移策略设计约束条件。
[0133]
步骤c：以轨道转移所需速度增量最小为原则，进行最优变轨次数估计。
[0134]
步骤d：如图4所示，建立描述单次变轨状态变化及多次变轨间联系的数学模型。
[0135]
步骤e：以geo轨道为目标，通过迭代寻优快速生成满足设计约束的轨道转移策略。
[0136]
步骤a包括：
[0137]
步骤s2.1：如图2所示，组合体航天器通常包含多个舱段，并且每个舱段都携带推进系统。如图3所示，其异常分离包括运载器未能将组合体航天器送至预定轨道即提前分离的发射段异常与组合体航天器中某舱段未能将组合体航天器送至准地球同步轨道即提前
分离的转移段异常，为了进行航天器轨道转移策略设计，首先需要明确异常分离时刻的轨道状态。不失一般性，采用geo轨道通用形式进行描述，即半长轴a0、倾角i0、偏心率e0、地理经度λ0、地理纬度η0；
[0138]
步骤s2.2：此外，还需要明确转移策略设计对象(组合体航天器、主星或主星与推进舱)的相关平台参数，一般包括分离时刻重量m0、可用燃料剩余量m'0、发动机推力f0、发动机比冲isp0；
[0139]
步骤b中，结合实际工程需求归纳出3方面设计约束：地面测控条件对点火段地理经度跨度及测控时长约束、推力器特性、弧段损失及安全保护需求对单次最长点火时间约束、变轨对象可用燃料约束；
[0140]
步骤c中，为了确定轨道转移优化参数维数需要依照如下步骤以转移所需速度增量最小为原则，进行最优变轨次数估计：
[0141]
步骤s4.1：依照如下公式计算得到航天器从异常分离时刻状态转移至geo目标轨道的最小速度增量δv：
[0142][0143]
δv＝v
n
‑
v0[0144]
其中，r0为转移轨道设计对象异常分离时刻的地心距，可由其轨道根数计算得到；v0为转移轨道设计对象异常分离时刻速度；v
n
为目标轨道速度；a
n
为目标轨道半长轴，通常geo轨道取值为42164km；
[0145]
步骤s4.2：根据火箭公式，计算得到轨道转移最小速度增量对应的燃料消耗δm：
[0146][0147]
其中，g＝9.80665m/s2为地球引力加速度；
[0148]
步骤s4.3：根据发动机比冲计算推进剂秒流量dm，并结合燃料消耗计算总共的变轨时长t：
[0149][0150][0151]
步骤s4.4：结合推力器单次最长点火时间t
max
计算得到最优变轨次数估计值n：
[0152][0153]
其中[x]为屋檐函数，表示取≥x的最小整数。
[0154]
步骤d包括：
[0155]
步骤s5.1：根据每次点火待优化变量，即变轨目标半长轴a
k
、变轨前偏移圈次q
k
，k＝1,
…
,n计算点火前参数：
[0156]
轨道周期t
k
：
[0157][0158]
经度漂移率
[0159][0160]
其中，ω
e
＝7.2921
×
10
‑5rad/s
[0161]
升交点地理经度λ
k
：
[0162][0163]
步骤s5.2：根据每次变轨目标半长轴a
k
，k＝1,
…
,n，计算得到单次点火点及目标轨道参数：
[0164]
由活力公式，计算目标轨道速度v
k
：
[0165][0166]
其中，r
k
为转移轨道设计对象点火时刻的地心距，可由其轨道根数计算得到；
[0167]
变轨速度增量δv
k
：
[0168][0169]
其中，α与β
k
依照如下公式计算
[0170][0171]
β
k
＝π
‑
i
k
‑1‑
α
[0172]
变轨后倾角i
k
：
[0173][0174]
变轨燃料消耗δm
k
：
[0175][0176]
点火时长t
k
：
[0177][0178]
步骤e包括：
[0179]
步骤s6.1：将geo轨道转移策略求解问题转化为多变量、多目标、多约束寻优问题，优化模型描述如下：
[0180]
目标函数：
[0181]
[0182]
其中，表示轨道目标罚函数；表示第j个轨道参数变轨终值与目标轨道偏差；δorb＝(δa,δi,δe,δλ,δη)，表示变轨结果与目标轨道偏差；γ表示权重系数，取较大正整数，用于对不满足精度指标参数进行大权重惩罚；
[0183]
约束条件：
[0184][0185]
步骤s6.2：针对上述优化模型选择现有多变量、多目标、多约束优化算法进行迭代寻优；
[0186]
步骤s6.3：输出寻优结果，确定异常分离条件的组合体航天器geo轨道转移策略。
[0187]
本实施例中，假设异常发生在运载器与组合体航天器分离时刻，由于运载器故障导至分离时刻轨道倾角发生偏差，依照步骤a确定分离轨道参数：
[0188]
orb0＝(a0,i0,e0,λ0,η0)＝(24473.64km,30
°
,0.731215,
‑
9.21
°
,30
°
)
[0189]
确定组合体与轨道转移相关的平台参数为：分离时刻重量m0＝3863kg、可用燃料剩余量m'0＝1820kg、发动机推力f0＝350n、发动机比冲isp0＝315s；
[0190]
依照步骤b确定轨道转移策略设计约束：最短测控时长t
fire
≥30min、地面测控地理经度范围[30
°
e，170
°
e]、发动机单次点火最长时间t
engine
≤4050s、轨道转移最多可用燃料
[0191]
依照步骤c，结合geo目标轨道参数orb
n
＝(42164km,30
°
,0,60
°
,0
°
)对变轨次数进行估计：
[0192][0193]
最后依照步骤e所建立的优化模型，结合步骤d进行单次变轨状态变化及多次变轨间参数更新，通过迭代寻优生成如表1所示的满足设计约束的轨道转移策略。
[0194]
表1轨道转移策略迭代寻优解
[0195]
[0196][0197]
本发明还提供一种异常分离情况下组合体geo变轨策略生成系统，包括：
[0198]
模块a：根据组合体与运载器或组合体之间异常分离时刻状态，确定组合体分离点轨道参数及与组合体轨道转移策略生成相关的平台参数。
[0199]
模块b：根据平台参数确定轨道转移策略的设计约束条件。
[0200]
模块c：根据所述约束条件以轨道转移所需速度增量最小为原则，进行最优变轨次数估计。
[0201]
模块d：根据估计结果，建立描述单次变轨状态变化及多次变轨间联系的数学模型。
[0202]
模块e：利用所述数学模型以geo轨道为目标，通过迭代寻优生成满足设计约束的轨道转移策略。
[0203]
本领域技术人员知道，除了以纯计算机可读程序代码方式实现本发明提供的系统及其各个装置、模块、单元以外，完全可以通过将方法步骤进行逻辑编程来使得本发明提供的系统及其各个装置、模块、单元以逻辑门、开关、专用集成电路、可编程逻辑控制器以及嵌入式微控制器等的形式来实现相同功能。所以，本发明提供的系统及其各项装置、模块、单元可以被认为是一种硬件部件，而对其内包括的用于实现各种功能的装置、模块、单元也可以视为硬件部件内的结构；也可以将用于实现各种功能的装置、模块、单元视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
[0204]
以上对本发明的具体实施例进行了描述。需要理解的是，本发明并不局限于上述特定实施方式，本领域技术人员可以在权利要求的范围内做出各种变化或修改，这并不影响本发明的实质内容。在不冲突的情况下，本技术的实施例和实施例中的特征可以任意相互组合。