CN108490963B - 全电推进卫星电推力器故障模式下的位置保持方法及系统 - Google Patents
全电推进卫星电推力器故障模式下的位置保持方法及系统 Download PDFInfo
- Publication number
- CN108490963B CN108490963B CN201810128311.4A CN201810128311A CN108490963B CN 108490963 B CN108490963 B CN 108490963B CN 201810128311 A CN201810128311 A CN 201810128311A CN 108490963 B CN108490963 B CN 108490963B
- Authority
- CN
- China
- Prior art keywords
- electric thruster
- ignition
- control
- electric
- thruster
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 60
- 238000012423 maintenance Methods 0.000 claims abstract description 13
- 230000008878 coupling Effects 0.000 claims abstract description 9
- 238000010168 coupling process Methods 0.000 claims abstract description 9
- 238000005859 coupling reaction Methods 0.000 claims abstract description 9
- 238000005259 measurement Methods 0.000 claims abstract description 5
- 230000001133 acceleration Effects 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 238000010304 firing Methods 0.000 claims description 3
- 238000012804 iterative process Methods 0.000 claims description 3
- 230000001174 ascending effect Effects 0.000 claims 2
- 238000004364 calculation method Methods 0.000 abstract description 4
- 239000003380 propellant Substances 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 5
- 238000009434 installation Methods 0.000 description 3
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000003203 everyday effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course, altitude or attitude of land, water, air or space vehicles, e.g. using automatic pilots
- G05D1/08—Control of attitude, i.e. control of roll, pitch, or yaw
- G05D1/0808—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft
- G05D1/0816—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability
- G05D1/0825—Control of attitude, i.e. control of roll, pitch, or yaw specially adapted for aircraft to ensure stability using mathematical models
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Algebra (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Automation & Control Theory (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开了一种全电推进卫星电推力器故障模式下的位置保持方法及系统,该方法包括如下步骤:当全电推进卫星某台电推力器发生故障时,该电推力器所在的分支不再使用,则使用另一分支的两个电推力器进行位置保持控制;其中,所述位置保持控制包括如下步骤:步骤一:根据测轨数据计算各轨道要素所需控制量;步骤二:根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;步骤三:计算电推力器点火位置偏角;步骤四:计算电推力器点火速度增量以及点火弧段中点赤经;步骤五:计算电推力器点火时刻及点火时长。本发明可适用于星上自主计算,解决电推力器故障模式下的位置保持问题,并实现较高的推进剂利用效率。
Description
技术领域
本发明属于卫星轨道姿态动力学与控制技术领域,尤其涉及一种全电推进卫星电推力器故障模式下的位置保持方法及系统。
背景技术
全电推进卫星(All-Electric Propulsion Satellite)采用高比冲的电推进系统(如离子电推进或霍尔电推进等),取代复杂的双组元统一化学推进系统,来实现星箭分离后的卫星转移轨道变轨、卫星入轨后的在轨位置保持、角动量卸载及离轨等任务。全电推进卫星由于电推力器配置数量、安装位置以及在轨供电功率的限制,其东西、南北位置保持控制相互耦合,且电推力器故障模式下控制参数减少,位置保持策略设计非常复杂。
全电推进卫星一般配置4台电推力器,呈矩形安装于卫星背地板或南北板,包含2个分支,每个分支包含对角线2台电推力器,2个分支互为备份。正常模式下,4台电推力器均能参与位置保持任务,每天(或每2天)每台电推力器点火1次。位置保持控制方式可通过地面定期向星上注入电推力器点火计划,由星上自主执行;或完全由星上根据位置保持策略自主计算电推力器点火参数,自主执行点火过程。当星上某台电推力器发生故障时,则需要制定新的位置保持策略,以尽量少的推进剂消耗实现位置保持控制。
北京理工大学CN104015938A号专利提供了一种电推进静止轨道卫星的位置保持方法,该方法针对配置4台矩形构型电推力器的卫星,实现了东西、南北位置保持的联合控制,达到了较高的位置保持精度和较高的推进剂利用效率。但是该方法没有给出电推力器故障模式下的位置保持策略,对于全电推进卫星电推力器发生故障情况下的位置保持,该方法并不适用。
发明内容
本发明解决的技术问题是:克服现有技术的不足,提供了一种全电推进卫星电推力器故障模式下的位置保持方法及系统,可适用于星上自主计算,解决电推力器故障模式下的位置保持问题,并实现较高的推进剂利用效率。
本发明目的通过以下技术方案予以实现:根据本发明的一个方面,提供了一种全电推进卫星电推力器故障模式下的位置保持方法,所述方法包括如下步骤:当全电推进卫星某台电推力器发生故障时,该电推力器所在的分支不再使用,则使用另一分支的两个电推力器进行位置保持控制;其中,所述位置保持控制包括如下步骤:步骤一:根据测轨数据计算各轨道要素所需控制量;步骤二:根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;步骤三:计算电推力器点火位置偏角;步骤四:计算电推力器点火速度增量以及点火弧段中点赤经;步骤五:根据点火弧段中点赤经得到电推力器点火弧段中点时刻,根据电推力器点火速度增量得到电推力器点火时长,根据电推力器点火弧段中点时刻和电推力器点火时长得到电推力器点火起始时刻。
上述全电推进卫星电推力器故障模式下的位置保持方法中,在步骤一中,各轨道要素所需控制量包括倾角矢量控制量(Δix,Δiy)、偏心率矢量控制量(Δex,Δey)和平经度漂移率控制量ΔD;其中,
倾角控制量(ΔixC,ΔiyC)按照如下控制律计算:
其中,ns为太阳视运动角速率,ix0、iy0为倾角矢量测量值,T为控制周期,λs=nst为太阳矢量与春分点方向的夹角,t为相对春分时刻的积秒,Ax=0.0235°、Ay=0.0215°分别为倾角矢量半年周期摄动项的幅值,Ωms为白道升d交点黄经,以18.6年为周期在0-360°范围内变化,其变化可记为Ωms=12.11°-0.052954T,其中T为相对儒略日,相对于1950年1月1日。
上述全电推进卫星电推力器故障模式下的位置保持方法中,偏心率矢量控制量(Δex,Δey)按照如下控制律计算:
其中,(ex0,ey0)为偏心率矢量,is=23.45°为黄赤交角,Re=0.011CrS/m为偏心率摄动圆半径,Cr为光压系数,S/m为卫星面积质量比。
上述全电推进卫星电推力器故障模式下的位置保持方法中,平经度控制需经过2个控制周期才能实现,平经度漂移率控制量ΔD包括2个控制周期的平经度漂移率控制量,2个控制周期的平经度漂移率控制量分别为ΔD1、ΔD2,平经度漂移率控制量ΔD1、ΔD2按照如下步骤计算:
首先计算平经度与漂移率的控制目标(λ*,D*):
采用优化算法进行求解,根据平经度与漂移率摄动变化方程:
首先给定ΔD1、ΔD2一组初值,固定ΔD1求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD2,再固定ΔD2求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD1,最终使实际控制结果(λf,Df)与控制目标(λ*,D*)的误差满足一定要求,即得到平经度所需的平经度漂移率控制量ΔD1、ΔD2,其中,D0为平经度漂移率,λ0为平经度。
上述全电推进卫星电推力器故障模式下的位置保持方法中,在步骤二中,总偏心率矢量控制量的公式如下:
上述全电推进卫星电推力器故障模式下的位置保持方法中,在步骤三中,电推力器点火位置偏角的公式如下:
σ=σe-σb
其中,σe为偏心率矢量偏角:
σb为径向速度增量所产生的偏心率矢量修正量偏角:
上式中包含σ角,即需要迭代过程,δ的迭代初值可选为:
上述全电推进卫星电推力器故障模式下的位置保持方法中,在步骤四中,电推力器点火速度增量的公式如下:
点火弧段中点赤经的公式如下:
其中a、a’、b、b’分别为:
其中,Rs为静止轨道半径,Vs为静止轨道卫星飞行速度,ΔV1为第一个电推力器第1次点火产生的速度增量,ΔV1'为第一个电推力器第2次点火产生的速度增量,ΔV4为第一个电推力器所在对角线上的另一电推力器第1次点火产生的速度增量,ΔV4'为第一个电推力器所在对角线上的另一电推力器第2次点火产生的速度增量,l1为第一个电推力器第1次点火过程的弧段中点赤经,l1'为第一个电推力器第2次点火过程的弧段中点赤经,l4为第一个电推力器所在对角线上的另一电推力器第1次点火过程的弧段中点赤经,l4'为第一个电推力器所在对角线上的另一电推力器第2次点火过程的弧段中点赤经,a、a’、b、b’均为转换量。
上述全电推进卫星电推力器故障模式下的位置保持方法中,在步骤五中,电推力器点火弧段中点时刻的公式如下:
其中,t0为过春分点时刻,ωe为地球自转角速率,li为点火弧段中点赤经。
上述全电推进卫星电推力器故障模式下的位置保持方法中,在步骤五中,电推力器点火时长的公式如下:
其中mc为航天器重量,Fp为电推力器推力,ΔVi为电推力器点火速度增量。
根据本发明的另一方面,还提供了一种全电推进卫星电推力器故障模式下的位置保持系统,所述系统包括:第一模块,用于根据测轨数据计算各轨道要素所需控制量;第二模块,用于根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;第三模块,用于计算电推力器点火位置偏角;第四模块,用于计算电推力器点火速度增量以及点火弧段中点赤经;第五模块,用于根据点火弧段中点赤经得到电推力器点火弧段中点时刻,根据电推力器点火速度增量得到电推力器点火时长,根据电推力器点火弧段中点时刻和电推力器点火时长得到电推力器点火起始时刻。
本发明与现有技术相比具有如下有益效果:
(1)本发明提供了一种简单易行的全电推进卫星电推力器故障模式下的位置保持方法,计算过程简单,电推力器点火参数明确,适合星上自主执行;
(2)本发明所提供的全电推进卫星电推力器故障模式下的位置保持方法,能够以较高的推进剂利用效率实现卫星位置保持。
附图说明
通过阅读下文优选实施方式的详细描述,各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的,而并不认为是对本发明的限制。而且在整个附图中,用相同的参考符号表示相同的部件。在附图中:
图1是本发明实施例提供的全电推进卫星电推力器故障模式下的位置保持方法的流程示意图;
图2是本发明实施例提供的星体电推力器布局示意图;
图3是本发明实施例提供的电推力器点火位置示意图;
图4是本发明实施例提供的一年内东西方向的位置保持控制结果的示意图;
图5是本发明实施例提供的一年内南北方向的位置保持控制结果的示意图。
具体实施方式
下面将参照附图更详细地描述本公开的示例性实施例。虽然附图中显示了本公开的示例性实施例,然而应当理解,可以以各种形式实现本公开而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了能够更透彻地理解本公开,并且能够将本公开的范围完整的传达给本领域的技术人员。需要说明的是,在不冲突的情况下,本发明中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本发明。
图2是本发明的实施例所采用的星体电推力器布局示意图,四台电推力器安装在卫星背地板,按安装位置区分为NW(1)、NE(2)、SW(3)、SE(4)。图3是本发明的实施例中电推力器点火位置示意图,图中圆周表示一个轨道周期,X轴表示春分点方向,Δt1、Δt2表示NW(1)、NE(2)两台电推力器的点火弧段,均位于约赤经90°附近,Δt3、Δt4表示SW(3)、SE(4)两台电推力器的点火弧段,均位于约赤经270°附近。图4为实施例中一年内东西方向的位置保持控制结果,图中横坐标为相对初始历元时刻的天数,纵坐标为卫星地理经度。图5为实施例中一年内南北方向的位置保持控制结果,图中横坐标为相对初始历元时刻的天数,纵坐标为卫星地理纬度。
图1是本发明实施例提供的全电推进卫星电推力器故障模式下的位置保持方法的流程示意图。如图1所示,该方法包括如下步骤:当全电推进卫星某台电推力器发生故障时,该电推力器所在的分支不再使用,则使用另一分支的两个电推力器进行位置保持控制;其中,所述位置保持控制包括如下步骤:
步骤一:根据测轨数据计算各轨道要素所需控制量;
步骤二:根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;
步骤三:计算电推力器点火位置偏角;
步骤四:计算电推力器点火速度增量以及点火弧段中点赤经;
步骤五:根据点火弧段中点赤经得到电推力器点火弧段中点时刻,根据电推力器点火速度增量得到电推力器点火时长,根据电推力器点火弧段中点时刻和电推力器点火时长得到电推力器点火起始时刻。
(1)步骤一、根据测轨数据计算各轨道要素所需控制量。各轨道要素测量值分别为倾角矢量(ix0,iy0)、偏心率矢量(ex0,ey0)、平经度λ0、平经度漂移率D0,再根据下述轨道要素控制律得到控制量。
倾角矢量控制量:
其中T为控制周期,λs=nst为太阳矢量与春分点方向的夹角,Ax=0.0235°,Ay=0.0215°分别为倾角矢量半年周期摄动项的幅值。Ωms为白道升d交点黄经,以18.6年为周期在0~360°范围内变化,其变化可记为Ωms=12.11°-0.052954T,其中T为相对儒略日,相对于1950年1月1日。
偏心率矢量控制量:
其中is=23.45°为黄赤交角。
平经度漂移率控制量ΔD:首先计算平经度与漂移率的控制目标(λ*,D*):
其中,为卫星定点位置等效的平经度漂移加速度;tday为时间,为86400s;平经度控制需经过2个控制周期才能实现,平经度漂移率控制量ΔD包括2个控制周期的控制量,即2个控制周期的控制量分别为ΔD1、ΔD2,可采用优化算法进行求解,根据平经度与漂移率摄动变化方程:
首先给定ΔD1、ΔD2一组初值,固定ΔD1求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD2,再固定ΔD2求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD1,最终使实际控制结果(λf,Df)与控制目标(λ*,D*)的误差满足一定要求,即得到平经度所需的控制量ΔD1、ΔD2。
(2)步骤二、根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量:
(3)步骤三、计算电推力器点火位置偏角:
σ=σe-σb
其中σe为偏心率矢量偏角:
其中σb为径向速度增量所产生的偏心率矢量修正量偏角:
上式中包含σ角,即需要迭代过程,δ的迭代初值可选为:
上述迭代过程只需执行1~3步即可。
(4)步骤四、计算电推力器点火速度增量以及点火弧段中点赤经:
其中a、a’、b、b’分别为:
其中,Rs为静止轨道半径,Vs为静止轨道卫星飞行速度,ΔV1为第一个电推力器第1次点火产生的速度增量,ΔV1'为第一个电推力器第2次点火产生的速度增量,ΔV4为第一个电推力器所在对角线上的另一电推力器第1次点火产生的速度增量,ΔV4'为第一个电推力器所在对角线上的另一电推力器第2次点火产生的速度增量,l1为第一个电推力器第1次点火过程的弧段中点赤经,l1'为第一个电推力器第2次点火过程的弧段中点赤经,l4为第一个电推力器所在对角线上的另一电推力器第1次点火过程的弧段中点赤经,l4'为第一个电推力器所在对角线上的另一电推力器第2次点火过程的弧段中点赤经,a、a’、b、b’均为转换量。
(5)步骤五、计算电推力器点火时刻及点火时长。电推力器点火弧段中点赤经可以转换为相应的时刻:
其中t0为过春分点时刻,ωe为地球自转角速率。再根据各电推力器所需产生的速度增量求得实际的电推力器点火时长:
其中mc为航天器重量,Fp为电推力器推力。与电推力器点火弧段中点时刻作差可得到电推力器点火起始时刻:
重复执行以上步骤,即可实现全电推进卫星电推力器故障模式下的位置保持。
具体的,卫星电推力器布局如附图2所示,电推力器初始安装角度为θ=50°,α=20°,电推力器沿轨道切向、径向、法向的投影系数分别为:KT=0.2620、KR=0.7198、KN=0.6428。卫星初始重量3000kg,有效面积75m2,光压系数1.5。
卫星运行过程中电推力器NE(2)故障,其所在分支不再工作,仅通过NW(1)与SE(4)两台电推力器进行位置保持。其中每2天为1个控制周期,每台电推力器点火2次,点火顺序分别为NW(1)点火、SE(4)点火、NW(1)第二次点火、SE(4)第二次点火。
(1)步骤一、根据测轨数据计算各轨道要素所需控制量。测轨数据:(ix0=0,iy0)=(0,0.01°),偏心率矢量(ex0,ey0)=(0,0)、平经度λ0=-0.021°、平经度漂移率D0=1.81×10-10,由倾角控制律得到每个控制周期的倾角控制量:
ΔixC=0.0002°
ΔiyC=-0.0067°
偏心率控制量为:
Δex=-0.2769×10-5
Δey=-0.0242×10-5
平经度漂移率控制量为:
ΔD1=-8.4591×10-11
ΔD2=-1.7534×10-9
其中ΔD1是当前控制周期的平经度漂移率控制量、ΔD2是下一个控制周期的平经度漂移率控制量。
(2)步骤二、根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量:
(3)步骤三、计算电推力器点火位置偏角:
σ=σe-σb=2.7725°
以及电推力器点火位置半夹角:
上式为迭代3次之后的计算结果。
(4)步骤四、计算电推力器点火速度增量以及点火弧段中点赤经:
(5)步骤五、计算电推力器点火时刻及点火时长。
测轨数据历元时刻为2020年1月1日00:00:00,过春分点时刻为2020年1月1日09:17:34,相对历元秒长为35454s,将电推力器点火弧段中点赤经转换为相对历元秒长:
t1=65648s
t1'=136654s
t4=108735s
t4'=179731s
电推力器各次点火时长:
Δt1=8787.6s
Δt1'=648.5s
Δt4=1461.4s
Δt4'=8052.7s
再求得电推力器点火起始时刻(以第1个控制周期为例):
tstart1=52830s
tstart1'=138997s
tstart4=95765s
tstart4'=181903s
重复以上步骤,在电推力器控制下卫星将保持在定点位置附近东西、南北±0.05°的范围内。
图4-图5为一年内卫星位置保持与角动量卸载联合控制结果。由图4可知采用本方法,可在单台电推力器故障情况下完成东西方向位置保持,满足±0.05°的经度控制范围要求。由图5可知采用本方法,可在单台电推力器故障情况下完成南北方向位置保持,满足±0.05°的纬度控制范围要求。
本实施例还提供了一种全电推进卫星电推力器故障模式下的位置保持系统,该系统包括:第一模块、第二模块、第三模块、第四模块和第五模块。其中,
第一模块,用于根据测轨数据计算各轨道要素所需控制量;第二模块,用于根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;第三模块,用于计算电推力器点火位置偏角;第四模块,用于计算电推力器点火速度增量以及点火弧段中点赤经;第五模块,用于根据点火弧段中点赤经得到电推力器点火弧段中点时刻,根据电推力器点火速度增量得到电推力器点火时长,根据电推力器点火弧段中点时刻和电推力器点火时长得到电推力器点火起始时刻。
以上所述的实施例只是本发明较优选的具体实施方式,本领域的技术人员在本发明技术方案范围内进行的通常变化和替换都应包含在本发明的保护范围内。
Claims (7)
1.一种全电推进卫星电推力器故障模式下的位置保持方法,其特征在于,所述方法包括如下步骤:
当全电推进卫星某台电推力器发生故障时,该电推力器所在的分支不再使用,则使用另一分支的两个电推力器进行位置保持控制;其中,
所述位置保持控制包括如下步骤:
步骤一:根据测轨数据计算各轨道要素所需控制量;
步骤二:根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;
步骤三:计算电推力器点火位置偏角;
步骤四:计算电推力器点火速度增量以及点火弧段中点赤经;
步骤五:根据点火弧段中点赤经得到电推力器点火弧段中点时刻,根据电推力器点火速度增量得到电推力器点火时长,根据电推力器点火弧段中点时刻和电推力器点火时长得到电推力器点火起始时刻;
在步骤一中,各轨道要素所需控制量包括倾角矢量控制量(ΔixC,ΔiyC)、偏心率矢量控制量(Δex,Δey)和平经度漂移率控制量ΔD;其中,
倾角矢量控制量(ΔixC,ΔiyC)按照如下控制律计算:
其中,ns为太阳视运动角速率,ix0、iy0为倾角矢量测量值,T为控制周期,λs=nst为太阳矢量与春分点方向的夹角,t为相对春分时刻的积秒,Ax=0.0235°、Ay=0.0215°分别为倾角矢量半年周期摄动项的幅值,Ωms为白道升交点黄经,以18.6年为周期在0-360°范围内变化,其变化可记为Ωms=12.11°-0.052954T1,其中T1为相对儒略日,相对于1950年1月1日;
偏心率矢量控制量(Δex,Δey)按照如下控制律计算:
其中,(ex0,ey0)为偏心率矢量,is=23.45°为黄赤交角,Re=0.011CrS/m为偏心率摄动圆半径,Cr为光压系数,S/m为卫星面积质量比;
平经度控制需经过2个控制周期才能实现,平经度漂移率控制量ΔD包括2个控制周期的平经度漂移率控制量,2个控制周期的平经度漂移率控制量分别为ΔD1、ΔD2,平经度漂移率控制量ΔD1、ΔD2按照如下步骤计算:
首先计算平经度与漂移率的控制目标(λ*,D*):
采用优化算法进行求解,根据平经度与漂移率摄动变化方程:
首先给定ΔD1、ΔD2一组初值,固定ΔD1求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD2,再固定ΔD2求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD1,最终使实际控制结果(λf,Df)与控制目标(λ*,D*)的误差满足一定要求,即得到平经度所需的平经度漂移率控制量ΔD1、ΔD2,其中,D0为平经度漂移率,λ0为平经度;
在步骤二中,总偏心率矢量控制量的公式如下:
3.根据权利要求2所述的全电推进卫星电推力器故障模式下的位置保持方法,其特征在于:在步骤四中,电推力器点火速度增量的公式如下:
点火弧段中点赤经的公式如下:
其中a、a’、b、b’分别为:
其中,Rs为静止轨道半径,Vs为静止轨道卫星飞行速度,ΔV1为第一个电推力器第1次点火产生的速度增量,ΔV1'为第一个电推力器第2次点火产生的速度增量,ΔV4为第一个电推力器所在对角线上的另一电推力器第1次点火产生的速度增量,ΔV4'为第一个电推力器所在对角线上的另一电推力器第2次点火产生的速度增量,l1为第一个电推力器第1次点火过程的弧段中点赤经,l1'为第一个电推力器第2次点火过程的弧段中点赤经,l4为第一个电推力器所在对角线上的另一电推力器第1次点火过程的弧段中点赤经,l4'为第一个电推力器所在对角线上的另一电推力器第2次点火过程的弧段中点赤经,a、a’、b、b’均为转换量。
7.一种全电推进卫星电推力器故障模式下的位置保持系统,其特征在于,所述系统包括:
第一模块,用于根据测轨数据计算各轨道要素所需控制量;
第二模块,用于根据倾角控制量以及南北、东西位置保持耦合关系计算总偏心率矢量控制量;
第三模块,用于计算电推力器点火位置偏角;
第四模块,用于计算电推力器点火速度增量以及点火弧段中点赤经;
第五模块,用于根据点火弧段中点赤经得到电推力器点火弧段中点时刻,根据电推力器点火速度增量得到电推力器点火时长,根据电推力器点火弧段中点时刻和电推力器点火时长得到电推力器点火起始时刻;
各轨道要素所需控制量包括倾角矢量控制量(ΔixC,ΔiyC)、偏心率矢量控制量(Δex,Δey)和平经度漂移率控制量ΔD;其中,
倾角矢量控制量(ΔixC,ΔiyC)按照如下控制律计算:
其中,ns为太阳视运动角速率,ix0、iy0为倾角矢量测量值,T为控制周期,λs=nst为太阳矢量与春分点方向的夹角,t为相对春分时刻的积秒,Ax=0.0235°、Ay=0.0215°分别为倾角矢量半年周期摄动项的幅值,Ωms为白道升交点黄经,以18.6年为周期在0-360°范围内变化,其变化可记为Ωms=12.11°-0.052954T1,其中T1为相对儒略日,相对于1950年1月1日;
偏心率矢量控制量(Δex,Δey)按照如下控制律计算:
其中,(ex0,ey0)为偏心率矢量,is=23.45°为黄赤交角,Re=0.011CrS/m为偏心率摄动圆半径,Cr为光压系数,S/m为卫星面积质量比;
平经度控制需经过2个控制周期才能实现,平经度漂移率控制量ΔD包括2个控制周期的平经度漂移率控制量,2个控制周期的平经度漂移率控制量分别为ΔD1、ΔD2,平经度漂移率控制量ΔD1、ΔD2按照如下步骤计算:
首先计算平经度与漂移率的控制目标(λ*,D*):
采用优化算法进行求解,根据平经度与漂移率摄动变化方程:
首先给定ΔD1、ΔD2一组初值,固定ΔD1求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD2,再固定ΔD2求得实际控制结果(λf,Df)与控制目标(λ*,D*)最近的ΔD1,最终使实际控制结果(λf,Df)与控制目标(λ*,D*)的误差满足一定要求,即得到平经度所需的平经度漂移率控制量ΔD1、ΔD2,其中,D0为平经度漂移率,λ0为平经度;
总偏心率矢量控制量的公式如下:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810128311.4A CN108490963B (zh) | 2018-02-08 | 2018-02-08 | 全电推进卫星电推力器故障模式下的位置保持方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810128311.4A CN108490963B (zh) | 2018-02-08 | 2018-02-08 | 全电推进卫星电推力器故障模式下的位置保持方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108490963A CN108490963A (zh) | 2018-09-04 |
CN108490963B true CN108490963B (zh) | 2021-03-26 |
Family
ID=63339930
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810128311.4A Active CN108490963B (zh) | 2018-02-08 | 2018-02-08 | 全电推进卫星电推力器故障模式下的位置保持方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108490963B (zh) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109063380B (zh) * | 2018-09-12 | 2023-02-17 | 北京理工大学 | 一种静止轨道电推进卫星故障检测方法及位置保持方法 |
CN109765926B (zh) * | 2018-12-27 | 2021-12-07 | 宁波天擎航天科技有限公司 | 一种升交点赤经精确控制方法、系统及装置 |
CN111268177B (zh) * | 2020-03-10 | 2021-09-03 | 北京控制工程研究所 | 一种静止轨道卫星分散式闭环自主位置保持控制方法 |
CN111487065B (zh) * | 2020-04-24 | 2022-03-04 | 中国空间技术研究院 | 一种卫星控制系统及电推进系统闭环联试方法 |
CN112278330B (zh) * | 2020-09-27 | 2022-02-01 | 北京控制工程研究所 | 一种基于星时驱动的电推进位置保持方法 |
CN113581494B (zh) * | 2021-07-22 | 2023-06-06 | 中国长城工业集团有限公司 | 一种geo卫星化电混合推进变轨方法 |
CN113734469B (zh) * | 2021-08-30 | 2023-02-03 | 北京控制工程研究所 | 一种电推力器位置保持推力分配方法和系统 |
CN115682919B (zh) * | 2022-10-28 | 2023-08-11 | 哈尔滨工业大学 | 一种基于曲面积分的霍尔推力器推力矢量偏心计算方法 |
CN116946392B (zh) * | 2023-07-12 | 2024-03-05 | 中国西安卫星测控中心 | 基于多维姿态偏置的地球同步卫星电推进倾角控制方法 |
CN117406251B (zh) * | 2023-09-14 | 2024-04-30 | 中国西安卫星测控中心 | 基于相对位置信息解算相对轨道要素的方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941749A (zh) * | 2014-04-30 | 2014-07-23 | 北京控制工程研究所 | 东西位置保持与角动量卸载联合控制方法 |
CN105373133A (zh) * | 2015-11-26 | 2016-03-02 | 中国空间技术研究院 | 一种同步轨道电推进位置保持与角动量卸载联合控制方法 |
CN105511490A (zh) * | 2015-12-15 | 2016-04-20 | 北京理工大学 | 一种静止轨道卫星位置保持-角动量卸载联合控制方法 |
CN105883008A (zh) * | 2014-12-15 | 2016-08-24 | 中国空间技术研究院 | 卫星推力器布局方法 |
CN106628264A (zh) * | 2016-11-23 | 2017-05-10 | 中国空间技术研究院 | 一种针对全电推进卫星的推力器布局方法 |
CN106697332A (zh) * | 2016-11-23 | 2017-05-24 | 中国空间技术研究院 | 一种全电推进卫星平台位置保持电推力器冗余配置方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9487309B2 (en) * | 2013-05-20 | 2016-11-08 | Kratos Integral Holdings, Llc | Eccentricity control for geosynchronous satellites |
-
2018
- 2018-02-08 CN CN201810128311.4A patent/CN108490963B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103941749A (zh) * | 2014-04-30 | 2014-07-23 | 北京控制工程研究所 | 东西位置保持与角动量卸载联合控制方法 |
CN105883008A (zh) * | 2014-12-15 | 2016-08-24 | 中国空间技术研究院 | 卫星推力器布局方法 |
CN105373133A (zh) * | 2015-11-26 | 2016-03-02 | 中国空间技术研究院 | 一种同步轨道电推进位置保持与角动量卸载联合控制方法 |
CN105511490A (zh) * | 2015-12-15 | 2016-04-20 | 北京理工大学 | 一种静止轨道卫星位置保持-角动量卸载联合控制方法 |
CN106628264A (zh) * | 2016-11-23 | 2017-05-10 | 中国空间技术研究院 | 一种针对全电推进卫星的推力器布局方法 |
CN106697332A (zh) * | 2016-11-23 | 2017-05-24 | 中国空间技术研究院 | 一种全电推进卫星平台位置保持电推力器冗余配置方法 |
Non-Patent Citations (1)
Title |
---|
GEO卫星基于电推进系统的倾角与偏心率联合控制方法;李强 等;《中国空间科学技术》;20160625;第36卷(第3期);第77-84页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108490963A (zh) | 2018-09-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108490963B (zh) | 全电推进卫星电推力器故障模式下的位置保持方法及系统 | |
CN102880184B (zh) | 一种静止轨道卫星自主轨道控制方法 | |
CN105373133B (zh) | 一种同步轨道电推进位置保持与角动量卸载联合控制方法 | |
US8439312B2 (en) | System and methods for simultaneous momentum dumping and orbit control | |
US7918420B2 (en) | System and methods for simultaneous momentum dumping and orbit control | |
CN109063380B (zh) | 一种静止轨道电推进卫星故障检测方法及位置保持方法 | |
JP6271043B2 (ja) | 軌道制御装置および衛星 | |
EP3112273B1 (en) | Efficient stationkeeping design for mixed fuel systems in response to a failure of an electric thruster | |
RU2721036C2 (ru) | Эффективная схема удержания на орбите для систем на смешанном топливе | |
CN111268177B (zh) | 一种静止轨道卫星分散式闭环自主位置保持控制方法 | |
CN114394260B (zh) | 地球同步卫星小推力南北保持控制解析计算方法 | |
Zhao et al. | Minimum-fuel station-change for geostationary satellites using low-thrust considering perturbations | |
CN112061424B (zh) | 一种基于融合目标姿态的机动过程能源角动态跟踪方法 | |
CN112357122B (zh) | 一种带有充气天线的卫星 | |
CN113447043B (zh) | 一种基于gnss的卫星天文导航系统误差自主标定方法及系统 | |
Varma | Control of satellites using environmental forces: aerodynamic drag/solar radiation pressure | |
RU2535353C2 (ru) | Способ удержания космического аппарата на геосинхронной 24-часовой орбите | |
CN114063645A (zh) | 基于偏心率倾角矢量的倾斜绕飞保持控制效果评估方法 | |
CN109507909B (zh) | 基于混合动力轨控技术的卫星共位分析建模方法及系统 | |
CN108490966B (zh) | 基于微分代数的静止轨道摄动相对轨迹高阶制导方法 | |
Somov et al. | Guidance and control of a space robot at additional launching and approaching a communication geostationary satellite | |
Chao | Semiautonomous stationkeeping of geosynchronous satellites | |
Buckingham et al. | Orbit position control using solar pressure. | |
Pitchaimani et al. | Strategies for Enhancing Power Generation in Satellites with Large Local Time Angles | |
Hintz et al. | Error analyses for the delivery of a spinning probe to Jupiter |
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 |