CN105539881B - 一种仅使用一对斜对称推力器的位置保持优化方法 - Google Patents
一种仅使用一对斜对称推力器的位置保持优化方法 Download PDFInfo
- Publication number
- CN105539881B CN105539881B CN201510937585.4A CN201510937585A CN105539881B CN 105539881 B CN105539881 B CN 105539881B CN 201510937585 A CN201510937585 A CN 201510937585A CN 105539881 B CN105539881 B CN 105539881B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- msup
- delta
- mtr
- 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
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B64—AIRCRAFT; AVIATION; COSMONAUTICS
- B64G—COSMONAUTICS; VEHICLES OR EQUIPMENT THEREFOR
- B64G1/00—Cosmonautic vehicles
- B64G1/22—Parts of, or equipment specially adapted for fitting in or to, cosmonautic vehicles
- B64G1/24—Guiding or controlling apparatus, e.g. for attitude control
- B64G1/242—Orbits and trajectories
Abstract
本发明公开的一种仅使用一对斜对称推力器的位置保持优化方法,涉及一种静止轨道卫星的位置保持优化方法,属于卫星轨道控制技术领域。本发明可以应用于使用一对斜对称推力器的静止轨道卫星长期运行管理任务。本发明在一个位置保持周期内斜对称推力器各开机两次,同时控制轨道倾角、平经度飘移率和偏心率,共五个等式方程,将推力器的开机赤经和开机速度增量作为自由变量,相应的会有八个待优化自由自变量,通过优化求解相应的开机赤经和开机速度增量,完成若干次的位置保持周期,即可实现仅使用一对斜对称推力器的位置保持。本发明可以以燃料较优的方式仅利用一对斜对称推力器实现位置保持,即解决部分推力器失效的静止轨道卫星位置飘移的问题。
Description
技术领域
本发明涉及一种静止轨道卫星的位置保持优化方法,特别涉及一种仅使用一对斜对称安装推力器的位置保持优化方法,属于卫星轨道控制技术领域。
背景技术
由于地球非球形引力、日月第三体引力以及太阳光压等摄动的影响,静止轨道卫星相对标称点位置发生漂移。地球非球形引力摄动主要引起卫星经度漂移率的变化使卫星在经度方向上的漂移,日月引力摄动主要引起卫星倾角的变化使在纬度方向上的漂移,而太阳光压摄动主要引起卫星偏心率的变化使卫星在经度方向和径向方向震动。随着时间的推移,这些摄动力对卫星产生的影响会使卫星逐渐漂移出定点保持窗口。因此,静止轨道卫星需要在一定时间间隔内进行轨道机动来抵消摄动因素的影响使卫星能够在任务周期内始终保持在卫星的定点窗口内(李恒年.地球静止卫星轨道与共位控制技术[M].北京:国防工业出版社,2010.)。
电推进具有比冲高、质量轻、体积小,特别是消耗工质少等优点,用于航天器能显著提高有效载荷比,若增加推进剂携带量则可延长工作寿命。电推进发动机可以长时间持续工作,开关机基本不受限制,一般在几十到几百毫牛顿量级(高扬.电火箭星际航行:技术进展、轨道设计与综合优化[J].力学学报,2011,43(6):991~1019.)。两组斜对称安装在卫星的背地板上,分别为NW-SE和NE-SW,每组斜对称推力器包括一个北推力器和一个南推力器,两组斜对称的推力器构成四方形,可以实现在轨位置保持(Anzel B.,Stationkeepingthe Hughes HS 702Satellite with a Xenon Ion Propulsion System,IAF-98-A.1.09,49th International Astronautical Congress,Sept 28-Oct 2,1998,Melbourne,Australia.)。
但在某个推力器出现故障后,与其对称安装的另一个推力器也无法工作,仅能利用另外剩余的推力器对完成位置保持控制,所以仅有两对推力器中的一对可以工作:NW(1)-SE(4)或NE(2)-SW(3),这就是故障模式。如推力器NE(2)出现故障,则由推力器对NW(1)-SE(4)实现位置保持控制。在故障情况下,将推力器的开机的4个位置赤经作为变量,再加上4个推力器速度增量,相应的会有8个自变量,而被控目标轨道要素有5个变量,分别为轨道倾角矢量(2个变量)、偏心率矢量(2个变量)和经度漂移率(1个变量),相应的有5个等式方程,自变量个数大于等式方程约束的个数,通过优化求解8个变量,优化的目标是推力器总速度增量最小,即燃料最优。
发明内容
本发明公开的一种仅使用一对斜对称推力器的位置保持优化方法要解决的技术问题是,提供一种仅使用一对斜对称推力器即可实现静止轨道卫星位置保持优化方法,此外,所述的方法能够以燃料最优的方式解决部分推力器失效的静止轨道卫星东西和南北方向位置飘移的问题,即实现静止轨道卫星位置保持。
本发明的目的是通过下述技术方案实现的。
本发明公开的一种仅使用一对斜对称推力器的位置保持优化方法,在一个位置保持周期内斜对称推力器各开机两次,同时控制轨道倾角、平经度飘移率和偏心率,将推力器的开机赤经和开机速度增量作为自由变量,通过优化求解相应的开机赤经和开机速度增量,完成若干次的位置保持周期,即可实现仅使用一对斜对称推力器的位置保持。
本发明公开的一种仅使用一对斜对称推力器的位置保持优化方法,具体步骤如下:
步骤一、确定一个位置保持周期轨道保持所需的初始轨道要素以及摄动对轨道要素一个位置保持周期的改变量,作为轨道保持的计算输入。
根据地面站观测数据确定卫星的初始瞬时轨道要素,获取卫星一个位置保持周期内第一次经过开机处的平均轨道要素作为初始轨道要素。轨道保持所需的轨道要素包括轨道倾角矢量(ix,iy)、偏心率(ex,ey)、经度λ以及经度飘移率D,通过对轨道观测数据的处理可以获得初始平均轨道要素,分别为ix0、iy0、ex0、ey0、λ0和D0。通过对地球非球形引力、日月引力和太阳光压等摄动进行分析获得平均轨道要素一个位置保持周期内的改变量,分别为ΔixD、ΔiyD,ΔexD、ΔeyD和ΔDD。
步骤二、确定一个位置保持周期内所需的轨道倾角矢量改变量(Δix,Δiy)。
轨道倾角矢量(ix,iy)可以利用轨道法向的控制力直接改变,一个位置保持周期内轨道倾角矢量所需改变量(Δix,Δiy)为:
步骤三、确定一个位置保持周期内所需的偏心率矢量改变量(Δex,Δey)。
偏心率矢量可以利用轨道切向和径向的控制力直接改变,一个位置保持周期内偏心率矢量所需改变量为:
步骤四、确定一个位置保持周期内所需的经度飘移率改变量ΔD。
定义四个推力器分为两对斜对称推力器,推力器NW(1)和推力器NE(2)为北推力器,推力器SW(3)和推力器SE(4)为南推力器,推力器NW(1)和推力器SE(4)组成一对斜对称推力器,推力器NE(2)和推力器SW(3)组成一对斜对称推力器。以推力器对NW(1)-SE(4)为例说明仅使用一对斜对称推力器实现静止轨道卫星位置保持,推力器对NE(2)-SW(3)同理实现。
经度通过改变经度飘移率间接控制,而每个位置保持周期的经度飘移率的目标量由此位置保持周期的初始经度偏差ΔλO、一个位置保持周期内经度速度的改变量ΔDD和推力径向分量导致的经度改变量DR决定:
其中由初始经度偏差ΔλO决定。
位置保持过程中,经度变化范围相对于360°的经度区间较小,从而可以认为经度漂移加速度为常值。从而一个位置保持周期内漂移速度改变量ΔDD为:
ΔDD=ΓλTD (4)
其中,Γλ是地球非球形引力田谐项导致的经度变化的加速度,TD为一个位置保持周期的时间。
在一个位置保持周期由于速度增量径向分量引起的经度改变量DR为:
DR=2kR(ΔV1-+ΔV4-+ΔV1-'+ΔV4-')/(TDVs) (5)
其中,Vs是静止轨道平均速度,ΔV1-,ΔV4-,ΔV1-'和ΔV4-'为上一个位置保持周期四次开机导致的速度增量。
相应的经度飘移率的改变量ΔD为:
ΔD=Dt-D0 (6)
其中D0是初始经度漂移率。
步骤五、确定一个位置保持周期内所需的四个推力器的控制速度增量与经度漂移率的关系。推力器NW(1)在第一个和第二个象限工作的赤经分别为l1'和l1,导致的速度增量分别用ΔV1'和ΔV1代表;推力器SE(4)在第三个和第四个象限导致工作的赤经分别为l4'和l4,导致的速度增量分别用ΔV4'和ΔV4代表。令kT、kN以及kR代表推力器推力在切向、法向以及径向的推力分量大小。
推力器NW(1)工作导致的速度增量ΔV1和ΔV1'的切向分量符号为正,推力器SE(4)工作导致的速度增量ΔV4和ΔV4'的切向分量符号为负,四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与经度漂移率ΔD的关系为,
其中Rs是静止轨道标称半径。
步骤六、确定一个位置保持周期内所需的四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与轨道倾角矢量改变量(Δix,Δiy)的关系。
速度增量ΔV1在轨道法向的分量符号为负,速度增量ΔV4在轨道法向的分量符号为正。速度增量ΔV1的赤经为l1∈[90° 180°],所以sinl1=s1>0,速度增量ΔV4的赤经为l4∈[270° 360°],所以sinl4=-s4<0,速度增量ΔV1'在轨道法向的分量符号为负,速度增量ΔV4'在轨道法向的分量符号为正。速度增量ΔV1'的赤经为l1'∈[0° 90°],所以sinl1'=s1'>0,速度增量ΔV4'的赤经为l4'∈[180° 270°],所以sinl4'=-s4'<0,从而四次速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与轨道倾角矢量改变量(Δix,Δiy)的关系为
步骤七、确定一个位置保持周期内所需的四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与偏心率矢量改变量(Δex,Δey)的关系。
速度增量ΔV1和ΔV4在轨道径向的分量符号均为负。速度增量ΔV1的赤经为l1∈[90° 180°],所以sinl1=s1>0,速度增量ΔV4的赤经为l4∈[270°360°],所以sinl4=-s4<0,速度增量ΔV1'和ΔV4'在轨道径向的分量符号均为负。速度增量ΔV1'的赤经为l1'∈[0° 90°],所以sinl1'=s1'>0,速度增量ΔV4'的赤经为l4'∈[180° 270°],所以sinl4'=-s4'>0,从而四次速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与偏心率改变量(Δex,Δey)的关系为
步骤八、求解一个位置保持周期内所需的四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')和推力器工作赤经(l1,l1',l4和l4')。
以上方程组(7)、(8)和(9)组合共有5个等式方程,而自由变量有8个,分别为4个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')和4个推力器工作赤经(l1,l1',l4和l4'),自变量的个数多于被控量的个数,此处以总速度增量最小目标求解上述的方程组(7)、(8)和(9),即可实现燃料最优,其目标函数为,
J=ΔV1+ΔV4+ΔV1'+ΔV4' (10)
对于速度增量幅值存在如下不等式约束。
对于开机处赤经存在如下不等式约束。
上述的步骤二、步骤三和步骤四确定一个位置保持周期的轨道要素改变量(Δix,Δiy,Δex,Δey和ΔD),以方程式(10)作为目标函数,方程式(7)、(8)和(9)作为等式约束,方程式(11)和(12)作为不等式约束,利用常用数值求解软件便可求解一个位置保持周期内所需的四个控制速度增量(ΔV1,ΔV1',ΔV4和ΔV4')和推力器工作赤经(l1,l1',l4和l4')。
步骤九、确定一个位置保持周期内的推力器的开机时长以及开机时刻,将此作为控制系统的输入量,控制推力器开关机。
上述的步骤八在冲量假设下进行规划,获得了每次速度增量作用位置和大小。则推力器连续开机时长可以由Tsiolkovsky方程获得,从而可以利用开机时长对半分原则确定推力器开关机时刻。
步骤十、重复步骤一至步骤九,完成若干次的位置保持周期,即可实现仅使用一对斜对称推力器的位置保持,即解决部分推力器失效的静止轨道卫星东西和南北方向位置飘移的问题。
步骤四所述的由初始经度偏差ΔλO决定的计算公式优选公式(13):
所述的步骤八利用常用数值求解软件便可求解优选为利用MATLAB的自带函数fmincon便可求解。
有益效果
1、本发明的一种仅使用一对斜对称推力器的位置保持优化方法,在四个推力器中的部分推力器失效的情况下,可以仅使用剩余的一对斜对称推力器实现位置保持,即解决部分推力器失效的静止轨道卫星东西和南北方向位置飘移的问题。
2、本发明的一种仅使用一对斜对称推力器的位置保持优化方法,由于将四次开机的赤经和速度增量作为待优化变量联合求解,可以以燃料较优的方式利用一对斜对称推力器实现位置保持。
附图说明
图1为本发明的推力器布局方案;
图2为本发明中推力器对NW(1)-SE(4)工作的示意图;
图3为本发明中推力器对NE(2)-SW(3)工作的示意图。
具体实施方式
为了更好的说明本发明的目的和优点,下面结合附图和实例对发明内容做进一步说明。
实施例1:
以斜对称推力器对NW(1)-SE(4)工作为例介绍使用一对斜对称推力器的静止轨道卫星位置保持优化方法,斜对称推力器对NE(2)-SW(3)工作的原理与其类似。一个位置保持周期内推力器NW(1)和SE(4)各开机两次,为了以燃料最大效率实现轨道倾角修正,北推力器应在赤经90°附近工作,南推力器应在赤经270°附近工作。北推力器NW(1)分别在第一和第二象限工作;南推力器SE(4)分别在第三和第四象限工作。一个位置保持周期内,开机顺序为4→1’→1→4’,前两个的速度增量小于后两个的速度增量,施加前两次速度增量后不会造成较大的经度偏移。对于斜对称推力器对NE(2)-SW(3)工作,开机顺序为1天的2→3’→3→2’。
以一天为一个位置保持周期,一个位置保持周期内推力器NW(1)和SE(4)各开机两次,实现对轨道倾角矢量、偏心率矢量和经度漂移率的控制。
获取卫星第一次赤经270°处的平均轨道要素作为初始轨道要素。由于斜对称推力器NW(1)-SE(4)第一次开机的赤经l4接近270°,所以以卫星第一次经过赤经270°处的平均轨道要素作为初始轨道要素;对于斜对称推力器NE(2)-SW(3)第一次开机的赤经l2接近90°,所以以卫星第一次经过赤经90°处的平均轨道要素作为初始轨道要素。轨道保持所需的轨道要素包括轨道倾角矢量(ix,iy),偏心率(ex,ey),经度(λ)以及经度飘移率(D),通过对轨道观测数据的处理可以获得初始平均轨道要素,分别为ix0、iy0、ex0、ey0、λ0和D0。通过对地球非球形引力、日月引力和太阳光压等摄动进行分析获得平均轨道要素一天内的改变量,分别为ΔixD、ΔiyD,ΔexD、ΔeyD和ΔDD。
步骤二、确定一个位置保持周期内所需的轨道倾角矢量改变量。
轨道倾角矢量可以利用轨道法向的控制力直接改变,一个位置保持周期内轨道倾角矢量所需改变量为:
步骤三、确定一个位置保持周期内所需的偏心率矢量改变量。
偏心率矢量可以利用轨道切向和径向的控制力直接改变,一个位置保持周期内偏心率矢量所需改变量为:
步骤四、确定一个位置保持周期内所需的经度飘移率改变量。
经度通过改变经度飘移率间接控制,而每个位置保持周期的经度飘移率的目标量由此位置保持周期的初始经度偏差、一个位置保持周期内经度速度的改变量和推力径向分量导致的经度改变量决定。
其中由初始经度偏差ΔλO决定,本专利采用下式。
位置保持过程中,经度变化范围较小,从而可以认为经度漂移加速度为常值。从而一天内漂移速度改变量为:
ΔDD=ΓλTD (18)
其中,Γλ是地球非球形引力田谐项导致的经度变化的加速度,TD为一天的时间。
在一个位置保持周期由于速度增量径向分量引起的经度改变量为
DR=2kR(ΔV1-+ΔV4-+ΔV1-'+ΔV4-')/(TDVs) (19)
其中,Vs是静止轨道平均速度,ΔV1-,ΔV4-,ΔV1-'和ΔV4-'为上一个位置保持周期四次开机导致的速度增量。
相应的经度飘移率的改变量为:
ΔD=Dt-D0 (20)
其中D0是初始经度漂移率。
步骤五、确定一个位置保持周期内所需的四个推力器的控制速度增量与经度漂移率的关系。推力器NW(1)在第一个和第二个象限工作的赤经分别为l1'和l1,导致的速度增量分别用ΔV1'和ΔV1代表;推力器SE(4)在第三个和第四个象限导致工作的赤经分别为l4'和l4,导致的速度增量分别用ΔV4'和ΔV4代表。令kT、kN以及kR代表推力器推力在切向、法向以及径向的推力分量大小。
ΔV1和ΔV1'的切向分量符号为正,ΔV4和ΔV4'的切向分量符号为负,四个速度增量与经度漂移率的关系为:
其中Rs是静止轨道标称半径。
步骤六、确定一个位置保持周期内所需的四个推力器的控制速度增量与轨道倾角矢量改变量的关系。
ΔV1在轨道法向的分量符号为负,ΔV4在轨道法向的分量符号为正。ΔV1的赤经为l1∈[90° 180°],所以sinl1=s1>0,ΔV4的赤经为l4∈[270° 360°],所以sinl4=-s4<0,ΔV1'在轨道法向的分量符号为负,ΔV4'在轨道法向的分量符号为正。ΔV1'的赤经为l1'∈[0° 90°],所以sinl1'=s1'>0,ΔV4'的赤经为l4'∈[180° 270°],所以sinl4'=-s4'<0,从而四次速度增量对轨道倾角矢量的改变量为:
步骤七、确定一个位置保持周期内所需的四个推力器的控制速度增量与偏心率矢量改变量的关系。
ΔV1和ΔV4在轨道径向的分量符号均为负。ΔV1的赤经为l1∈[90° 180°],所以sinl1=s1>0,ΔV4的赤经为l4∈[270° 360°],所以sinl4=-s4<0,ΔV1'和ΔV4'在轨道径向的分量符号均为负。ΔV1'的赤经为l1'∈[0°90°],所以sinl1'=s1'>0,ΔV4'的赤经为l4'∈[180° 270°],所以sinl4'=-s4'>0,从而四次速度增量对偏心率的改变量为
步骤八、求解一个位置保持周期内所需的四个推力器的控制速度增量和推力器工作赤经。
以上方程组(21)、(22)和(23)组合共有5个等式方程,而自由变量有8个,分别为4个速度增量和4个推力器工作赤经,自变量的个数多于被控量的个数,此处将利用优化理论求解以上方程,其中目标函数为:
J=ΔV1+ΔV4+ΔV1'+ΔV4' (24)
对于速度增量幅值存在如下不等式约束。
对于开机处赤经存在如下不等式约束。
在轨道要素改变量确定的情况下,以式(24)作为目标函数,方程式(21)、(22)和(23)作为等式约束,方程式(25)和(26)作为不等式约束,利用MATLAB的自带函数fmincon便可求解。
步骤九、确定一个位置保持周期内的推力器的开机时长以及开机时刻,将此作为控制系统的输入量,控制推力器开关机。
步骤八在冲量假设下进行规划,获得了每次速度增量作用位置和大小。则推力器连续开机时长可以由Tsiolkovsky方程获得,从而可以利用开机时长对半分原则确定开机时刻。
步骤十、重复步骤一至步骤九,完成若干次的位置保持周期,即可实现使用一对斜对称推力器的位置保持。
以上所述的具体描述,对发明的目的、技术方案和有益效果进行进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (3)
1.一种仅使用一对斜对称推力器的位置保持优化方法,其特征在于:包括如下步骤:
步骤一、确定一个位置保持周期轨道保持所需的初始轨道要素以及摄动对轨道要素一个位置保持周期的改变量,作为轨道保持的计算输入;
根据地面站观测数据确定卫星的初始瞬时轨道要素,获取卫星一个位置保持周期内第一次经过开机处的平均轨道要素作为初始轨道要素;轨道保持所需的轨道要素包括轨道倾角矢量(ix,iy)、偏心率(ex,ey)、经度λ以及经度飘移率D,通过对轨道观测数据的处理可以获得初始平均轨道要素,分别为ix0、iy0、ex0、ey0、λ0和D0;通过对地球非球形引力、日月引力和太阳光压摄动进行分析获得平均轨道要素一个位置保持周期内的改变量,分别为ΔixD、ΔiyD,ΔexD、ΔeyD和ΔDD;
步骤二、确定一个位置保持周期内所需的轨道倾角矢量改变量(Δix,Δiy);
轨道倾角矢量(ix,iy)可以利用轨道法向的控制力直接改变,一个位置保持周期内轨道倾角矢量所需改变量(Δix,Δiy)为,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;i</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&Delta;i</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>x</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;i</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&Delta;i</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>i</mi>
<mrow>
<mi>y</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤三、确定一个位置保持周期内所需的偏心率矢量改变量(Δex,Δey);
偏心率矢量可以利用轨道切向和径向的控制力直接改变,一个位置保持周期内偏心率矢量所需改变量为,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;e</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&Delta;e</mi>
<mrow>
<mi>x</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>e</mi>
<mrow>
<mi>x</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;e</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>&Delta;e</mi>
<mrow>
<mi>y</mi>
<mi>D</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>e</mi>
<mrow>
<mi>y</mi>
<mn>0</mn>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤四、确定一个位置保持周期内所需的经度飘移率改变量ΔD;
定义四个推力器分为两对斜对称推力器,推力器NW(1)和推力器NE(2)为北推力器,推力器SW(3)和推力器SE(4)为南推力器,推力器NW(1)和推力器SE(4)组成一对斜对称推力器,推力器NE(2)和推力器SW(3)组成一对斜对称推力器;以推力器对NW(1)-SE(4)为例说明仅使用一对斜对称推力器实现静止轨道卫星位置保持,推力器对NE(2)-SW(3)同理实现;
经度通过改变经度飘移率间接控制,而每个位置保持周期的经度飘移率的目标量由此位置保持周期的初始经度偏差Δλ0、一个位置保持周期内经度速度的改变量ΔDD和推力径向分量导致的经度改变量DR决定;
<mrow>
<msub>
<mi>D</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<msub>
<mi>D</mi>
<mrow>
<msub>
<mi>&Delta;&lambda;</mi>
<mn>0</mn>
</msub>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;D</mi>
<mi>D</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>R</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中由初始经度偏差Δλ0决定;
位置保持过程中,经度变化范围相对于360°的经度区间较小,从而可以认为经度漂移加速度为常值;从而一个位置保持周期内漂移速度改变量ΔDD为,
ΔDD=ΓλTD (4)
其中,Γλ是地球非球形引力田谐项导致的经度变化的加速度,TD为一个位置保持周期的时间;
在一个位置保持周期由于速度增量径向分量引起的经度改变量DR为,
DR=2kR(ΔV1-+ΔV4-+ΔV1-'+ΔV4-')/(TDVs) (5)
其中,Vs是静止轨道平均速度,ΔV1-,ΔV4-,ΔV1-'和ΔV4-'为上一个位置保持周期四次开机导致的速度增量;
相应的经度飘移率的改变量ΔD为,
ΔD=Dt-D0 (6)
其中D0是初始经度漂移率;
步骤五、确定一个位置保持周期内所需的四个推力器的控制速度增量与经度漂移率的关系;推力器NW(1)在第一个和第二个象限工作的赤经分别为l1'和l1,导致的速度增量分别用ΔV1'和ΔV1代表;推力器SE(4)在第三个和第四个象限导致工作的赤经分别为l4'和l4,导致的速度增量分别用ΔV4'和ΔV4代表;令kT、kN以及kR代表推力器推力在切向、法向以及径向的推力分量大小;
推力器NW(1)工作导致的速度增量ΔV1和ΔV1'的切向分量符号为正,推力器SE(4)工作导致的速度增量ΔV4和ΔV4'的切向分量符号为负,四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与经度漂移率ΔD的关系为,
<mrow>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
<mo>=</mo>
<mo>-</mo>
<mfrac>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mn>3</mn>
</mfrac>
<mi>&Delta;</mi>
<mi>D</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Rs是静止轨道标称半径;
步骤六、确定一个位置保持周期内所需的四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与轨道倾角矢量改变量(Δix,Δiy)的关系;
速度增量ΔV1在轨道法向的分量符号为负,速度增量ΔV4在轨道法向的分量符号为正;速度增量ΔV1的赤经为l1∈[90° 180°],所以sinl1=s1>0,速度增量ΔV4的赤经为l4∈[270° 360°],所以sinl4=-s4<0,速度增量ΔV1'在轨道法向的分量符号为负,速度增量ΔV4'在轨道法向的分量符号为正;速度增量ΔV1'的赤经为l1'∈[0° 90°],所以sinl1'=s1'>0,速度增量ΔV4'的赤经为l4'∈[180° 270°],所以sinl4'=-s4'<0,从而四次速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与轨道倾角矢量改变量(Δix,Δiy)的关系为,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;i</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<msub>
<mi>k</mi>
<mi>N</mi>
</msub>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</msqrt>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</msqrt>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;i</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>k</mi>
<mi>N</mi>
</msub>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤七、确定一个位置保持周期内所需的四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与偏心率矢量改变量(Δex,Δey)的关系;
速度增量ΔV1和ΔV4在轨道径向的分量符号均为负;速度增量ΔV1的赤经为l1∈[90°180°],所以sinl1=s1>0,速度增量ΔV4的赤经为l4∈[270° 360°],所以sinl4=-s4<0,速度增量ΔV1'和ΔV4'在轨道径向的分量符号均为负;速度增量ΔV1'的赤经为l1'∈[0° 90°],所以sinl1'=s1'>0,速度增量ΔV4'的赤经为l4'∈[180° 270°],所以sinl4'=-s4'>0,从而四次速度增量(ΔV1,ΔV1',ΔV4和ΔV4')与偏心率改变量(Δex,Δey)的关系为,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;e</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mo>-</mo>
<mn>2</mn>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<mn>2</mn>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>+</mo>
<mn>2</mn>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</msqrt>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;e</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>V</mi>
<mi>s</mi>
</msub>
</mfrac>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<mn>2</mn>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<mn>2</mn>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>-</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>R</mi>
</msub>
<msqrt>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mrow>
<mo>&prime;</mo>
<mn>2</mn>
</mrow>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<msub>
<mi>k</mi>
<mi>T</mi>
</msub>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
步骤八、求解一个位置保持周期内所需的四个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')和推力器工作赤经(l1,l1',l4和l4');
以上方程组(7)、(8)和(9)组合共有5个等式方程,而自由变量有8个,分别为4个速度增量(ΔV1,ΔV1',ΔV4和ΔV4')和4个推力器工作赤经(l1,l1',l4和l4'),自变量的个数多于被控量的个数,此处以总速度增量最小目标求解上述的方程组(7)、(8)和(9),即可实现燃料最优,其目标函数J为,
J=ΔV1+ΔV4+ΔV1'+ΔV4' (10)
对于速度增量幅值存在如下不等式约束,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msup>
<msub>
<mi>&Delta;V</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>&GreaterEqual;</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
对于开机处赤经存在如下不等式约束,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msup>
<msub>
<mi>s</mi>
<mn>1</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>0</mn>
<mo>&le;</mo>
<msup>
<msub>
<mi>s</mi>
<mn>4</mn>
</msub>
<mo>&prime;</mo>
</msup>
<mo>&le;</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
上述的步骤二、步骤三和步骤四确定一个位置保持周期的轨道要素改变量(Δix,Δiy,Δex,Δey和ΔD),以方程式(10)作为目标函数,方程式(7)、(8)和(9)作为等式约束,方程式(11)和(12)作为不等式约束,利用常用数值求解软件便可求解一个位置保持周期内所需的四个控制速度增量(ΔV1,ΔV1',ΔV4和ΔV4')和推力器工作赤经(l1,l1',l4和l4');
步骤九、确定一个位置保持周期内的推力器的开机时长以及开机时刻,将此作为控制系统的输入量,控制推力器开关机;
上述的步骤八在冲量假设下进行规划,获得了每次速度增量作用位置和大小;则推力器连续开机时长可以由Tsiolkovsky方程获得,从而可以利用开机时长对半分原则确定推力器开关机时刻;
步骤十、重复步骤一至步骤九,完成若干次的位置保持周期,即可实现仅使用一对斜对称推力器的位置保持,即解决部分推力器失效的静止轨道卫星东西和南北方向位置飘移的问题。
2.如权利要求1所述的一种仅使用一对斜对称推力器的位置保持优化方法,其特征在于:步骤四所述的由初始经度偏差Δλ0决定的计算公式(13),
3.如权利要求1或2所述的一种仅使用一对斜对称推力器的位置保持优化方法,其特征在于:所述的步骤八利用MATLAB的自带函数fmincon便可求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937585.4A CN105539881B (zh) | 2015-12-15 | 2015-12-15 | 一种仅使用一对斜对称推力器的位置保持优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510937585.4A CN105539881B (zh) | 2015-12-15 | 2015-12-15 | 一种仅使用一对斜对称推力器的位置保持优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105539881A CN105539881A (zh) | 2016-05-04 |
CN105539881B true CN105539881B (zh) | 2018-02-02 |
Family
ID=55819537
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510937585.4A Active CN105539881B (zh) | 2015-12-15 | 2015-12-15 | 一种仅使用一对斜对称推力器的位置保持优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105539881B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106697332B (zh) * | 2016-11-23 | 2018-10-19 | 中国空间技术研究院 | 一种全电推进卫星平台位置保持电推力器冗余配置方法 |
CN107600462B (zh) * | 2017-08-22 | 2019-05-24 | 长光卫星技术有限公司 | 一种基于时分复用方式的小卫星轨道控制方法 |
CN109443813B (zh) * | 2018-09-28 | 2020-06-09 | 中国空间技术研究院 | 卫星电推进矢量调节机构转动测试方法 |
CN110254753B (zh) * | 2019-06-04 | 2021-02-12 | 北京理工大学 | 一种地球静止轨道卫星电推力器及其布局优化方法 |
CN112507461B (zh) * | 2020-12-15 | 2021-07-23 | 北京航天自动控制研究所 | 一种运载火箭动力软着陆段发动机开机方法 |
CN113184220B (zh) * | 2021-04-21 | 2021-11-19 | 中国人民解放军63923部队 | 一种地球同步轨道通信卫星的轨道控制方法及装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5020746A (en) * | 1989-09-29 | 1991-06-04 | Hughes Aircraft Company | Method for satellite station keeping |
US5443231A (en) * | 1993-11-17 | 1995-08-22 | Hughes Aircraft Company | Method and apparatus for a satellite station keeping |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6908064B2 (en) * | 2002-10-15 | 2005-06-21 | Lockheed Martin Corporation | High-efficiency REA optimized stationkeeping |
GB2460443A (en) * | 2008-05-29 | 2009-12-02 | Inmarsat Global Ltd | Propulsion system arrangement for decommissioning a satellite |
CN102649481B (zh) * | 2012-04-23 | 2015-08-26 | 上海卫星工程研究所 | 一种基于双力矩四棱锥的8推力器布局设计方法 |
US9296494B1 (en) * | 2012-08-10 | 2016-03-29 | Lockheed Martin Corporation | Thruster orbit control method and configuration |
-
2015
- 2015-12-15 CN CN201510937585.4A patent/CN105539881B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5020746A (en) * | 1989-09-29 | 1991-06-04 | Hughes Aircraft Company | Method for satellite station keeping |
US5443231A (en) * | 1993-11-17 | 1995-08-22 | Hughes Aircraft Company | Method and apparatus for a satellite station keeping |
Also Published As
Publication number | Publication date |
---|---|
CN105539881A (zh) | 2016-05-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105539881B (zh) | 一种仅使用一对斜对称推力器的位置保持优化方法 | |
Hasnain et al. | Capturing near-Earth asteroids around Earth | |
Macdonald et al. | Solar sail mission applications and future advancement | |
Janhunen et al. | Fast E-sail Uranus entry probe mission | |
Genta et al. | Optimal low-thrust trajectories for nuclear and solar electric propulsion | |
Macdonald et al. | Low-thrust-enabled highly-non-Keplerian orbits in support of future Mars exploration | |
CN103853047A (zh) | 一种基于状态量反馈的小推力跟踪制导方法 | |
Singh | Very low earth orbit propellant collection feasibility assessment | |
Prado | Using tethered gravity-assisted maneuvers for planetary capture | |
Gajeri et al. | A Titan mission using the direct fusion drive | |
Struck | The feasibility of shading the greenhouse with dust clouds at the stable lunar Lagrange points | |
Polzin et al. | Enabling deep space science missions with Nuclear Thermal Propulsion | |
Millis | Energy, incessant obsolescence, and the first interstellar missions | |
Murray | Continuous Earth-Moon payload exchange using motorised tethers with associated dynamics | |
Co | Operationally Responsive Spacecraft Using Electric Propulsion | |
Nan et al. | Global 4-D trajectory optimization for spacecraft | |
Poncy et al. | A preliminary assessment of an orbiter in the Haumean system: How quickly can a planetary orbiter reach such a distant target? | |
Hogan et al. | General high-altitude orbit corrections using electrostatic tugging with charge control | |
Yuan et al. | Three decades of progress in China’s space High-Tech Program empowered by modern astrodynamics | |
Simmonds et al. | Mini-CHT powered Formation Flying Mission for Magnetic Reconnection Research in Space | |
Santos et al. | Optimal low-thrust trajectories to reach the asteroid apophis | |
Nordley et al. | Mass beam propulsion, an overview | |
Malerba et al. | Book of Abstracts of The Sixth International Conference on Tethers in Space | |
Kezerashvili | Exploration of the solar system and beyond using a thermonuclear fusion drive | |
Grandl et al. | Near Earth Asteroids-Prospection, Orbit Modification and Mining |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |