CN106094854B - 电磁编队卫星姿态和轨道相对控制方法 - Google Patents
电磁编队卫星姿态和轨道相对控制方法 Download PDFInfo
- Publication number
- CN106094854B CN106094854B CN201610561867.3A CN201610561867A CN106094854B CN 106094854 B CN106094854 B CN 106094854B CN 201610561867 A CN201610561867 A CN 201610561867A CN 106094854 B CN106094854 B CN 106094854B
- Authority
- CN
- China
- Prior art keywords
- satellite
- attitude
- magnetic
- dipole
- control
- 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
-
- 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/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
-
- 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
Abstract
本发明提供一种电磁编队卫星姿态和轨道相对控制方法,步骤1,包括以下步骤:建立姿态动力学方程;根据步骤1建立的四元数描述的线性形式的姿态动力学方程、体坐标系下的姿态动力学方程和刚体卫星的动力学方程,设计得到卫星姿控系统反馈控制律。步骤3:根据步骤2提出的反馈控制律,对电磁编队卫星姿态和轨道进行相对控制。优点为:建立了姿态动力学方程、相对运动动力学方程及电磁力和电磁力矩的模型,由此设计的卫星姿控系统反馈控制律、飞轮角动量磁卸载算法、大角度姿态机动方法、电磁编队飞行卫星运动控制器和磁偶极子的控制分配方法具有较好的工作性能,能够实现高效的电磁编队卫星姿态和轨道相对控制。
Description
技术领域
本发明属于电磁编队卫星姿态和轨道相对控制技术领域,具体涉及一种电磁编队卫星姿态和轨道相对控制方法。
背景技术
卫星电磁编队飞行(EMFF)具有无需推进剂,而且没有光学传感器污染等多种优点,在近十年得到了广泛关注。电磁编队飞行这一概念能够应用于多种未来的任务,例如分布式航天器,天文的干涉量度分析法,以及相机装配等。
Miller教授的团队首先提出了电磁编队飞行这一概念,随后,MIT的空间系统实验室开始对其进行研究,主要研究了磁偶极子模型并进行了动态分析、控制器的设计以及地面模拟试验,这些研究为电磁编队飞行的研究打下了坚实的基础。随后,Sedwick教授的团队提出了谐振感应近场发电系统,这是电磁编队飞行在微重力环境下的首次验证。近年来,在中国也进行了很多关于电磁编队飞行的研究,这些研究主要是关于电磁编队飞行的相对平衡状态、动力学性能和对其控制方面。研究提出了多种控制器,例如Zeng et al.于2012年提出的有限时间控制器,Cai et al.于2013年提出的误差反馈控制器,Zhang et al.于2014年提出的自适应控制器和Huang et al.于2015年提出的基于线性矩阵不等式的控制器等。
电磁编队卫星的姿态能够通过给分布式卫星配备三个正交的反作用飞轮来实现,分布式卫星的相对位置可以通过电磁相互作用来控制。以往的研究主要是假定能够对其姿态进行很好的控制,来研究运动动力学模型和对电磁编队飞行的控制。针对卫星姿态控制系统,主要都是用反作用飞轮来实现,反作用飞轮通过低速摩擦补偿和鲁棒控制器的设计来进行控制。通过以上设计方法得到的卫星姿态和轨道相对控制系统,仍然具有较大的局限性,难以更为高效的有针对性的对卫星姿态和轨道进行相对控制。
发明内容
针对现有技术存在的缺陷,本发明提供一种电磁编队卫星姿态和轨道相对控制方法,可有效提高电磁编队卫星姿态和轨道相对控制精度。
本发明采用的技术方案如下:
本发明提供一种电磁编队卫星姿态和轨道相对控制方法,包括以下步骤:步骤1:建立姿态动力学方程;姿态动力学方程包括:
(一)用四元数描述的线性形式的姿态动力学方程:
其中:
—描述航天器姿态四元数的一阶导数;
ωxb、ωyb、ωzb分别代表体坐标系三个轴的转动角速度;
q1=cos(φ/2);
q2=exsin(φ/2);
q3=eysin(φ/2);
q4=ezsin(φ/2);
φ—欧拉轴的旋转角;
[ex ey ez]—欧拉轴的单位矢量;
(二)体坐标系下的姿态动力学方程:
其中:
Ixb、Iyb、Izb分别为体坐标系下三个轴的转动惯量;
Lxb、Lyb、Lzb分别为体坐标系下三个轴的力矩;
(三)刚体卫星的动力学方程:
其中:
IT—卫星及飞轮整体的转动惯量;
—卫星的角加速度;
—飞轮的角加速度;
Iw—反作用飞轮的转动惯量;
ωB—卫星的角速度;
ωw—体坐标系下飞轮的角速度;
LB—体坐标系下的力矩;
步骤2:根据步骤1建立的四元数描述的线性形式的姿态动力学方程、体坐标系下的姿态动力学方程和刚体卫星的动力学方程,设计得到卫星姿控系统反馈控制律;
本步骤具体包括:
步骤2.1,假定当前四元数为q,期望四元数为qc,则误差四元数qe的定义为:
假定当前角速度为ω,期望角速度为ωc;则误差角速度ωe的定义为:
ωe=ω-ωc (2)
步骤2.2,基于误差四元数和误差角速度的姿态动力学模型可被描述为
其中:
I—整星的惯量矩阵;
—误差角加速度;
—期望角加速度;
H—整星的角动量,H=IBωB+Iwωw;
IB—卫星的转动惯量;
u—反作用飞轮的控制力矩,u=Iwωw;
Md—外部干扰力矩;
ME—地磁场引起的干扰力矩;
—误差四元数的一阶导数;
步骤2.3,将地球磁场近似为一个以地球为中心的偶极子;非线性控制系统的非线性部分能够利用反馈线性化方法被抵消,从而提出反馈控制律:
其中:
P,D—正定增益矩阵,通常被设计为对角矩阵;
qe123—误差四元数中的三个分量qe1、qe2、qe3组成的列向量,qe123=[qe1 qe2 qe3]T;
步骤3:根据步骤2提出的反馈控制律,对电磁编队卫星姿态和轨道进行相对控制。
优选的,步骤3具体为:
在根据步骤2提出的反馈控制律对电磁编队卫星姿态和轨道进行相对控制时,在控制过程中,实时判断是否发生磁饱和现象,如果发生,采用飞轮角动量磁卸载算法卸载磁力矩。
优选的,在控制过程中,实时判断是否发生磁饱和现象,如果发生,采用飞轮角动量磁卸载算法卸载磁力矩,具体为:
步骤3.1,假设反作用飞轮的冗余转动惯量为ΔH;在控制过程中,实时测量飞轮所处磁场B,判断是否满足以下的磁卸载条件:
45°<(B,ΔH)<135° (30)
如果满足,进一步判断磁场B是否垂直于ΔH,如果不垂直,则按式(37)计算得到磁卸载力矩值Tunloading,并将该磁卸载力矩值Tunloading施加于飞轮中,使其反向转动,进行磁卸载:
|B|—飞轮所处磁场B的大小;
K—增益系数;
如果B垂直于ΔH,则按下式(38)计算得到磁卸载力矩值Tunloading,并将该磁卸载力矩值Tunloading施加于飞轮中,使其反向转动,进行磁卸载:
Tunloading=-KΔH(K>0) (38)
优选的,还包括:
步骤4,设计大角度姿态机动方法,包括:
步骤4.1,假设机动前,即初始时刻t0时刻的理论姿态为q(t0),理论角速度为ω(t0)=[0 0 0]T;机动后,即tf时刻目标姿态为q(tf),目标角速度为ω(tf)=[0 0 0]T;
步骤4.2,根据下式计算得到控制力矩Lc:
其中:
IT—卫星及飞轮整体的转动惯量;
q123—q1、q2、q3组成的列向量,q123=[q1 q2 q3]T;
Kd—大角度姿态机动控制律的微分系数;
Kp—大角度姿态机动控制律的比例系数;
步骤4.3,将控制力矩Lc作用于飞轮,使其以最快速度按下式规划的路径进行姿态机动:
其中:
ωc(t)—期望角速度关于时间的函数;
—期望角加速度关于时间的函数;
qc(t)—期望四元数关于时间的函数;
—瞬时欧拉轴;
—绕欧拉轴旋转角加速度关于时间的函数;
—绕欧拉轴旋转角速度关于时间的函数;
θ(t)—绕欧拉轴旋转角度大小关于时间的函数。
优选的,还包括:
步骤5,建立相对运动动力学方程;
相对运动动力学方程为:
其中:
ρ—环绕星相对参考星的相对位置矢量;
—环绕星相对参考星的相对速度矢量;
—环绕星相对参考星的相对加速度矢量;
ωn—参考卫星的轨道角速度;
r—参考星在地惯系中的位置矢量;
um—电磁控制相对加速度;
R—干扰项和未建模的动力学部分;
C(ωn)∈R3×3为
R3×3—三维实数集;
非线性矩阵N(·)∈R3×3定义为如下等式:
μE—地球引力常数;
x、y、z—环绕星与参考星的相对位置矢量在Hill坐标系下的分量;
步骤6,根据步骤5建立的相对运动动力学方程,设计得到电磁编队飞行卫星运动控制器;
所述电磁编队飞行卫星运动控制器是一种终端滑模控制器,用于实现轨道跟踪误差保持在零附近。
优选的,轨道跟踪误差可被定义为
其中:
e—轨道跟踪位置误差;
—轨道跟踪速度误差;
ρd—目标相对位置矢量;
—目标相对速度矢量;
因此,设计的终端滑模控制器描述为:
其中:β>0,ki>0,i=1,2,3,0<γ<1;
us—终端滑模控制相对加速度;
—目标相对加速度矢量;
β—控制系数β>0;
γ—控制系数0<γ<1;
k—控制系数ki>0,i=1,2,3;
sat(·)—饱和函数,
φ0—一严格正常数,表示所设计的边界层的厚度;
S—滑模面,滑模面可描述为
|S|—滑模面的模。
优选的,为了消除振动,用连续的饱和函数代替了公式(51)中的饱和函数,连续的饱和函数可描述为
其中,设计系数τ满足0<τ<1;公式(51)可被写为
由此得到终端滑模控制器的另一种描述。
优选的,还包括:
步骤7,建立电磁力和电磁力矩的模型,包括:
步骤7.1,卫星j对卫星i的电磁力Fij和电磁力矩Tij可被描述为
其中:
μ0—真空磁导率,μ0=4π×10-7H/m;
rij—由卫星i到卫星j的相对位置矢量;
|rij|—由卫星i到卫星j的相对位置矢量的模;
μi—卫星i的磁偶极子的强度,μi=NiIiSi;
μj—卫星j的磁偶极子的强度,μj=NjIjSj;
Ni、Nj分别为偶极子A、B的绕线匝数;
Ii、Ij分别为偶极子A、B的线圈电流;
Si、Sj分别为偶极子A、B的线圈面积;
由公式(55)和公式(56),可得到如下三自由度的远场电磁模型:
其中:
—偶极子A受到的电磁力在三个轴方向的分量;
—偶极子A受到的电磁力矩在三个轴方向的分量;
μA—卫星A的磁偶极子强度;
μB—卫星B的磁偶极子强度;
α、β分别为偶极子A、B与ox轴之间的夹角;
χ、δ分别为偶极子A、B绕ox轴的转角;
步骤8,在步骤6设计的电磁编队飞行卫星运动控制器的基础上,基于步骤7的电磁力和电磁力矩的模型,为获取卫星间相对位置的控制过程中每颗卫星的磁偶极子,还设计两种计算磁偶极子的方法:
第一种是根据牛顿第三定律,将电磁编队飞行卫星的某一偶极子固定,计算可得其他偶极子的电磁力;
第二种是将偶极子的计算问题转换为一个优化问题,利用一种粒子群算法和序列二次规划方法混合算法以实现自由偶极子的最优化,最终实现电磁干扰力矩MT和地磁场干扰力矩ME最小的优化目标。
优选的,步骤7中,第一种计算磁偶极子的方法包括:
根据牛顿第三定律,电磁编队飞行卫星运动控制器所对应的公式(54)代表针对n个卫星的3(n-1)个标量公式,因此,系统在此约束下方程会有无穷多解,如果将电磁编队飞行卫星的某一偶极子固定,其他偶极子就可以通过如下公式进行计算
其中:
μ0—真空磁导率,μ0=4π×10-7H/m;
rij—由卫星i到卫星j的相对位置矢量;
|rij|—由卫星i到卫星j的相对位置矢量的模;
μi—卫星i的磁偶极子强度;
μj—卫星j的磁偶极子强度;
步骤7中,第二种计算磁偶极子的方法为:
建立以下目标函数:
μoptimal—最优磁矩;
arg(·)—幅值;
n—卫星个数;
MT—电磁干扰力矩;
ME—地磁场干扰力矩;
采用粒子群算法和序列二次规划方法混合算法以实现自由偶极子的最优化,其中,首先采用粒子群算法设定初始参数值,然后采用粒子群算法进行渐进估算,当收敛值满足预先设定的临界值时,将收敛值提供给序列二次规划方法,采用序列二次规划方法进行非线性约束下问题的非线性解算。
本发明提供的电磁编队卫星姿态和轨道相对控制具有以下优点:
本发明提供的电磁编队卫星姿态和轨道相对控制方法,充分考虑了电磁编队卫星姿态和轨道相对控制方法设计每个阶段需要考虑的因素,建立了姿态动力学方程、相对运动动力学方程及电磁力和电磁力矩的模型,由此设计的卫星姿控系统反馈控制律、飞轮角动量磁卸载算法、大角度姿态机动方法、电磁编队飞行卫星运动控制器和磁偶极子的控制分配方法具有较好的工作性能,能够实现高效的电磁编队卫星姿态和轨道相对控制。
附图说明
图1为本发明提供的电磁编队卫星姿态和轨道相对控制方法的流程示意图;
图2为远场电磁模型;
图3为俯仰控制回路框架图;
图4为角加速度随时间变化的曲线;
图5为姿态机动流程图。
具体实施方式
为了使本发明所解决的技术问题、技术方案及有益效果更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
结合图1,本发明提供一种电磁编队卫星姿态和轨道相对控制方法,在建立姿态动力学方程、建立相对运动动力学方程和建立电磁力和电磁力矩的模型的基础上,提出以下创新点:1、设计一种卫星姿控系统反馈控制律;2、设计了一种飞轮角动量磁卸载算法;3、设计了一种大角度姿态机动方法;4、设计一种电磁编队飞行卫星运动控制器;5、设计一种磁偶极子的控制分配。按照本发明设计方法,能够根据电磁编队卫星的姿态控制方法,实现电磁编队卫星的对日定向;同时还能够实现磁卸载,电磁编队卫星相位位置保持的非线性控制,以及磁矩分配。下面对上述五大部分分别详细介绍:
(一)卫星姿控系统反馈控制律
步骤1:建立姿态动力学方程;姿态动力学方程包括:用四元数描述的线性形式的姿态动力学方程、体坐标系下的姿态动力学方程和刚体卫星的动力学方程,其建立过程包括以下步骤:
S1.1,用四元数描述航天器的姿态,将四元数描述的姿态动力学方程描述为线性形式;
用四元数描述航天器的姿态为:
q=(q1 q2 q3 q4)T (1)
其中:
q—航天器姿态四元数;
q1=cos(φ/2);
q2=exsin(φ/2);
q3=eysin(φ/2);
q4=ezsin(φ/2);
φ—欧拉轴的旋转角;
[ex ey ez]—欧拉轴的单位矢量;
T—矩阵的转置;
四元数满足如下约束:
因此,四元数描述下的姿态动力学可以描述为公式(6)所示的线性形式:
其中:
—描述航天器姿态的四元数的一阶导数;
ωxb、ωyb、ωzb分别代表体坐标系的三个轴上的转动角速度;
S1.2,根据轨道坐标系到体坐标系的转换矩阵,得到体坐标系下的卫星角动量;
由轨道坐标系到体坐标系的转换矩阵BG为:
其中:
—航天器的滚转角;
θ—航天器的俯仰角;
ψ—航天器的偏航角;
且三个角均为欧拉角;
由于欧拉角很小,因此上式可以近似为:
惯性坐标系下的卫星角动量为HI=II·ωI,
II—惯性坐标系下的转动惯量;
ωI—惯性坐标系下的角速度;
下标I—惯性坐标系下的向量;
S1.3,根据角动量定理,得到体坐标系下滚转轴的动力学方程;
根据刚体的角动量定理,得到如下等式:
LI—惯性坐标系下的力矩;
定义惯性坐标系到体坐标系的转换矩阵为BI,那么体坐标系下的转动惯量IB、体坐标系下的角速度ωB和体坐标系下的对应力矩LB分别为:
IB=BIIIBI T,ωB=BI·ωI,LB=BI·LI (7)
因此,体坐标系下的卫星角动量HB为:
HB=IB·ωB=BIIIBI TBI·ωI=BI·HI (8)
根据角动量定理,能够得到:
因此:
Δθ、Δψ—航天器的滚转角、俯仰角、偏航角的角度变化量;则体坐标系下滚转轴的动力学方程为等式(14),其中变量都定义在体坐标系下:
其中:
体坐标系下的卫星角动量
体坐标系下的力矩
体坐标系下的卫星转动角速度
其中:Ixb、Iyb、Izb分别为体坐标系下三个轴的转动惯量;
Lxb、Lyb、Lzb分别为体坐标系下三个轴的力矩;
Hxb、Hyb、HZb分别为体坐标系下三个轴的角动量;
因此,在体坐标系下,姿态动力学方程可被写作:
其中:
IT—卫星及飞轮整体的转动惯量;
—卫星的角加速度;
—飞轮的角加速度;
Iw—反作用飞轮的转动惯量;
ωB—卫星的角速度;
ωw—体坐标系下飞轮的角速度;
LB—体坐标系下的力矩;
S1.4,由于电磁编队飞行卫星的姿态能够通过三个正交的反作用飞轮控制,最后推导了刚性卫星的动力学公式;
电磁编队飞行卫星的姿态能够通过三个正交的反作用飞轮控制,这三个反作用飞轮可以沿体坐标系的三个坐标轴安装;设反作用飞轮的转动惯量为Iw,三个轴的转动惯量分别为Iwx,Iwy和Iwz,体坐标系下飞轮的角速度为ωw,卫星的角速度为ωB,卫星及飞轮整体的转动惯量为IT;因此,
HB=ITωB+Iwωw (13)
那么,刚性卫星的动力学方程可被写作:
其中:
IT—卫星及飞轮整体的转动惯量;
—卫星的角加速度;
—飞轮的角加速度;
Iw—反作用飞轮的转动惯量;
ωB—卫星的角速度;
ωw—体坐标系下飞轮的角速度;
LB—体坐标系下的力矩;
步骤2:根据步骤1建立的四元数描述的线性形式的姿态动力学方程、体坐标系下的姿态动力学方程和刚体卫星的动力学方程,设计得到卫星姿控系统反馈控制律;
本步骤具体包括:
步骤2.1:假定当前四元数为q,期望四元数为qc。那么,误差四元数qe的定义为:
假定当前角速度为ω,期望角速度为ωc;那么,误差角速度ωe的定义为:
ωe=ω-ωc (2)
步骤2.2,基于误差四元数和误差角速度的姿态动力学模型可被描述为
I—整星的惯量矩阵;
—误差角加速度;
—期望角加速度;
H—整星的角动量,H=IBωB+Iwωw;
IB—卫星的转动惯量;
u—反作用飞轮的控制力矩,u=Iwωw;
Md—外部干扰力矩;
ME—地磁场引起的干扰力矩;
—误差四元数的一阶导数;
为便于计算ME,将地球磁场近似为一个以地球为中心的偶极子;非线性控制系统的非线性部分能够利用反馈线性方法被抵消,从而提出反馈控制律:
其中:
P,D—正定增益矩阵,通常被设计为对角矩阵;
qe123—误差四元数中的三个分量qe1、qe2、qe3组成的列向量,qe123=[qe1 qe2 qe3]T;
姿态控制系统中,在反馈线性化后的闭环控制是一个典型的二阶线性系统,因此能够确定增益矩阵P和D;卫星姿态控制回路中包括动力学模型、感应器、控制器和反作用飞轮模型。以俯仰控制回路为例,该闭环回路的框架图如图3所示;
反作用飞轮和陀螺的模型可以视为传递函数,反作用飞轮的传递函数为惯性环节:
ωwheel—飞轮的振荡角频率;
s—拉氏算子;
陀螺的传递函数Wgyro为振荡环节:
ωg—陀螺的振荡角频率;
ζg—陀螺的阻尼系数;
在图3中,姿态控制系统开环传递函数由动力学模型的传递函数、陀螺的传递函数、反作用飞轮的传递函数和PD控制环节组成,系统的开环传递函数Wopen可表示为
p—PD控制的环节比例系数;
d—PD控制的环节微分系数;
Iy—干扰力矩产生的转动惯量;
系统的闭环传递函数Wclose为
外部干扰力矩Ty(s)到姿态输出θ(s)的闭环传递函数Wdisturb为
从而能够得到输入θc与外部干扰力矩Ty、输出θout的简化关系式
如公式(24)所示,p的增加能够抑制干扰对姿态的影响;与之相反,p的增加导致高频段的增益增大,这会影响弹性振动的稳定性;因此,参数p需要合理设计。
(二)飞轮角动量磁卸载算法
步骤3:根据步骤2提出的反馈控制律,对电磁编队卫星姿态和轨道进行相对控制。具体为:在根据步骤2提出的反馈控制律对电磁编队卫星姿态和轨道进行相对控制时,在控制过程中,实时判断是否发生磁饱和现象,如果发生,采用飞轮角动量磁卸载算法卸载磁力矩。
电磁力矩是由卫星磁偶极子的相互作用和地球磁场产生的,对于低轨小卫星,可以利用反作用飞轮对磁力矩器进行磁卸载;假定反作用飞轮的冗余转动惯量为ΔH,则其磁动量为
K—增益系数;
B—飞轮所处磁场;
|B|—飞轮所处磁场B的大小;
飞轮的磁力矩Tw为:
如果B垂直于ΔH,那么
Tw=-KΔH(K>0) (27)
如果B不垂直于ΔH,B可以被分解为沿着ΔH的方向和垂直于ΔH的方向两部分,分别记为BΔH和BN。那么对卸载有效的力矩可写为
对卸载无效的力矩为
由等式(29)可得,如果B和ΔH之间的角度很小,那么不适合磁卸载;磁卸载通常在|b·Δh|<0.707时进行,其中b和Δh分别是B和ΔH的单位向量;因此,得到磁卸载条件:
45°<(B,ΔH)<135° (30)
如果将比例反馈控制律用于磁卸载,那么会得到陀螺速度ωb和体坐标系下飞轮的角速度ωw;因此,整体的角动量可被表示为
H=I·ωb+Iw·ωw (31)
轨道的角速度为ω0,则理论角动量可被写为
H0=I·ω0 T (32)
因此,多余的角动量为
ΔH=H-H0≈Iw·ωw (33)
通过如下比例反馈控制律,可以得到磁卸载的指令:
Tw=-KΔH (34)
由于Tw=M×B,因此
B×T=B×(M×B)=M(B·B)-B(B·M) (35)
如果M与B平行,那么不能够产生力矩,即B·M=0;因此,
M=B×Tw/B2 (36)
磁卸载的施加的力矩为
如果B垂直于ΔH,则按下式(38)计算得到磁卸载力矩值Tunloading,并将该磁卸载力矩值Tunloading施加于飞轮中,使其反向转动,进行磁卸载:
Tunloading=-KΔH(K>0) (38)
(三)设计大角度姿态机动方法
为使太阳能帆板产生足够电流,调整电磁编队飞行卫星的姿态使其保持对日定向是最好的方法;传统的过程是调整姿态,使其由三轴稳定到对日定向;反作用飞轮用于控制姿态,因此,最大力矩和最大转动速度会导致大角度机动产生新问题;在飞轮的最大力矩和转动速度的限制下,基于欧拉旋转理论,能够设计出最理想的机动算法;此外,根据误差四元数和角速度的反馈控制律能够实现跟踪控制;
步骤4.1,假设机动前,即初始时刻t0时刻的理论姿态为q(t0),理论角速度为ω(t0)=[0 0 0]T;机动后,即tf时刻目标姿态为q(tf),目标角速度为ω(tf)=[0 0 0]T;
步骤4.2,根据欧拉定理,通过将瞬时欧拉轴旋转θ角,能够实现由q(t0)到q(tf)的姿态机动,且机动路径是最短的;为了满足最大力矩和转动速度的约束,机动路径可以规划为
Lc—飞轮的转动力矩;
—飞轮的转动速度;
—飞轮的转动速度;
—绕欧拉轴旋转角加速度;
—瞬时欧拉轴;
假定飞轮的最大输出力矩为Lcmax,那么绕欧拉轴的最大角加速度为
为以最快速度实现姿态机动,机动路径可以规划为
其中:
ωc(t)—期望角速度关于时间的函数;
—期望角加速度关于时间的函数;
qc(t)—期望四元数关于时间的函数;
—瞬时欧拉轴;
—绕欧拉轴旋转角加速度关于时间的函数;
—绕欧拉轴旋转角速度关于时间的函数;
θ(t)—绕欧拉轴旋转角度大小关于时间的函数。
为满足飞轮的磁饱和约束,机动过程要满足如图4所示变化律;
如图4所示,确定t1和tf以确定机动路径;通过忽略机动期间的干扰,可以保持角动量;因此,
θf—机动的角度;
为了确保飞轮在机动过程中不会达到磁饱和,应确定机动时间tf;假定最大角速度为那么
t1—从初始时刻开始以最大角加速度机动的结束时刻;
t2—第二次以最大角加速度机动的结束时刻,机动角加速度与t0~t1相反;其中,t1和tf的计算过程流程图如图5所示;
基于动量守恒定律能够得到ωwt的计算方程
ωwt—t时刻飞轮的角速度;
ωBt—体坐标系下卫星t时刻的角速度;
ωw0—初始时刻飞轮的角速度;
ωB0—体坐标系下卫星初始时刻的角速度;
θt—t时刻机动的角度;
—t时刻机动的角速度;
那么,控制律可写为
其中:
IT—卫星及飞轮整体的转动惯量;
q123—q1、q2、q3组成的列向量,q123=[q1 q2 q3]T;
Kd—大角度姿态机动控制律的微分系数;
Kp—大角度姿态机动控制律的比例系数;
(四)电磁编队飞行卫星运动控制器
步骤5,建立相对运动动力学方程;
相对运动动力学方程为:
其中:
ρ—环绕星相对参考星的相对位置矢量;
—环绕卫星相对参考卫星的相对速度矢量;
—环绕卫星相对参考卫星的相对加速度矢量;
ωn—参考卫星的轨道角速度;
r—参考星在地惯系中的位置矢量;
um—电磁控制相对加速度;
R—干扰项和未建模的动力学部分;
C(ωn)∈R3×3为
R3×3—三维实数集;
非线性矩阵N(·)∈R3×3定义为如下等式:
μE—地球引力常数;
x、y、z—环绕星与参考星的相对位置矢量在Hill坐标系下的分量;
步骤6,根据步骤5建立的相对运动动力学方程,设计得到电磁编队飞行卫星运动控制器;
所述电磁编队飞行卫星运动控制器是一种终端滑模控制器,用于实现轨道跟踪误差保持在零附近。轨道跟踪误差可被定义为
其中:
e—轨道跟踪位置误差;
—轨道跟踪速度误差;
ρd—目标相对位置矢量;
—目标相对速度矢量;
因此,设计的终端滑模控制器描述为:
其中β>0,ki>0,i=1,2,3,0<γ<1;
us—终端滑模控制相对加速度;
—目标相对加速度矢量;
β—控制系数β>0;
γ—控制系数0<γ<1;
k—控制系数ki>0,i=1,2,3;
sat(·)—饱和函数,
φ0—一严格正常数,表示所设计的边界层的厚度;
S—滑模面,滑模面可描述为
|S|—滑模面的模。
注:为了消除振动,用连续的饱和函数代替了公式(51)中的饱和函数;
其中,设计系数τ满足0<τ<1;公式(51)可被写为
由此得到终端滑模控制器的另一种描述。
(五)磁偶极子的控制分配
步骤7,建立电磁力和电磁力矩的模型;
本发明提出了电磁力和电磁力矩的远场模型,如图5所示。远场模型在卫星距离大于线圈半径的6.76倍时是有效的,此时远场模型的误差小于10%;
卫星j对卫星i的电磁力Fij和电磁力矩Tij可被描述为
μ0—真空磁导率,μ0=4π×10-7H/m;
rij—由卫星i到卫星j的相对位置矢量;
|rij|—由卫星i到卫星j的相对位置矢量的模;
μi—卫星i的磁偶极子的强度,μi=NiIiSi;
μj—卫星j的磁偶极子的强度,μj=NjIjSj;
Ni、Nj分别为偶极子A、B的绕线匝数;
Ii、Ij分别为偶极子A、B的线圈电流;
Si、Sj分别为偶极子A、B的线圈面积;
由公式(55)和公式(56),可得到如下三自由度的远场电磁模型:
—偶极子A受到的电磁力在三个轴方向的分量;
—偶极子A受到的电磁力矩在三个轴方向的分量;
μA—卫星A的磁偶极子强度;
μB—卫星B的磁偶极子强度;
α、β分别为偶极子A、B与ox轴之间的夹角;
χ、δ分别为偶极子A、B绕ox轴的转角;
在公式(57)中,电磁力和电磁力矩由两颗卫星的磁偶极子乘积决定,并且随两卫星的相对距离和角度变化;电磁力与相对距离的四次方成反比,而电磁力矩与相对距离的三次方成反比;因此,电磁力和电磁力矩是强非线性化和强耦合的。
步骤8,在步骤6设计的电磁编队飞行卫星运动控制器的基础上,基于步骤7的电磁力和电磁力矩的模型,为获取卫星间相对位置的控制过程中每颗卫星的磁偶极子,还设计两种计算磁偶极子的方法:
第一种是根据牛顿第三定律,将电磁编队飞行卫星的某一偶极子固定在,计算可得其他偶极子的电磁力;
根据牛顿第三定律,电磁编队飞行卫星运动控制器所对应的公式(54)代表针对n个卫星的3(n-1)个标量公式,因此,系统在此约束下方程会有无穷多解,如果将电磁编队飞行卫星的某一偶极子固定,其他偶极子就可以通过如下公式进行计算
其中,μ0—真空磁导率,μ0=4π×10-7H/m;
rij—由卫星i到卫星j的相对位置矢量;
|rij|—由卫星i到卫星j的相对位置矢量的模;
μi—卫星i的磁偶极子强度;
μj—卫星j的磁偶极子强度;
第二种是将偶极子的计算问题转换为一个优化问题,利用一种粒子群算法和序列二次规划方法(PSO-SQP)混合算法以实现自由偶极子的最优化,最终实现电磁干扰力矩MT和地磁场干扰力矩ME最小的优化目标。
建立以下目标函数:
μoptimal—最优磁矩;
arg(·)—幅值;
n—卫星个数;
MT—电磁干扰力矩;
ME—地磁场干扰力矩;
本发明提出的粒子群算法和序列二次规划方法(PSO-SQP)混合算法以实现该自由偶极子的最优化,其中粒子群算法是一种渐进估算算法,而序列二次规划方法是一种能够高度适用于非线性约束下问题的非线性方法;这一混合算法利用粒子群算法设定初始参数值,将其提供给序列二次规划方法;其优点在于序列二次规划方法能够减少计算时间并且提高收敛性能;当收敛值满足预先设定的临界值时,应用序列二次规划方法。
本发明提供的电磁编队卫星姿态和轨道相对控制方法,充分考虑了电磁编队卫星姿态和轨道相对控制方法设计每个阶段需要考虑的因素,建立了姿态动力学方程、相对运动动力学方程及电磁力和电磁力矩的模型,由此设计的卫星姿控系统反馈控制律、飞轮角动量磁卸载算法、大角度姿态机动方法、电磁编队飞行卫星运动控制器和磁偶极子的控制分配方法具有较好的工作性能,能够实现高效的电磁编队卫星姿态和轨道相对控制。按照本发明提出的电磁编队卫星姿态和轨道相对控制方法,能够实现电磁编队卫星的对日定向;同时还能够实现磁卸载,电磁编队卫星相位位置保持的非线性控制,以及磁矩分配。具体具有以下优点:
(1)在飞轮的最大力矩和转动速度的限制下,基于欧拉旋转理论,能够设计出最理想的大角度姿态。此外,根据误差四元数和角速度的反馈控制律能够实现跟踪控制。将反作用飞轮用于控制姿态,能够实现最短机动路径的姿态机动,并满足飞轮的磁饱和约束,最终确定机动过程的控制律,实现姿态控制对日定向和大角度机动。
(2)提出了一种终端滑模控制器来实现轨道跟踪误差保持在零附近。设计的终端滑模控制器能够使给出的相对运动系统是大范围渐进稳定的,具有鲁棒性。
(3)提出两种计算磁偶极子的方法:一种是根据牛顿第三定律,将某一偶极子固定在电磁编队飞行卫星上,计算可得其他偶极子的电磁力;另一种是将偶极子的计算问题转换为一个优化问题,利用一种粒子群算法和序列二次规划方法(PSO-SQP)混合算法以实现该自由偶极子的最优化,最终实现飞轮力矩和地球磁场干扰力矩最小化的优化目标;这两种方法均能够实现磁矩的最优控制分配。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视本发明的保护范围。
Claims (9)
1.一种电磁编队卫星姿态和轨道相对控制方法,其特征在于,包括以下步骤:
步骤1:建立姿态动力学方程;姿态动力学方程包括:
(一)用四元数描述的线性形式的姿态动力学方程:
其中:
—描述航天器姿态四元数的一阶导数;
ωxb、ωyb、ωzb分别代表体坐标系三个轴的转动角速度;
q1=cos(φ/2);
q2=exsin(φ/2);
q3=eysin(φ/2);
q4=ezsin(φ/2);
φ—欧拉轴的旋转角;
[ex ey ez]—欧拉轴的单位矢量;
(二)体坐标系下的姿态动力学方程:
其中:
Ixb、Iyb、Izb分别为体坐标系下三个轴的转动惯量;
Lxb、Lyb、Lzb分别为体坐标系下三个轴的力矩;
(三)刚体卫星的动力学方程:
其中:
IT—卫星及飞轮整体的转动惯量;
—卫星的角加速度;
—飞轮的角加速度;
Iw—反作用飞轮的转动惯量;
ωB—卫星的角速度;
ωw—体坐标系下飞轮的角速度;
LB—体坐标系下的力矩;
步骤2:根据步骤1建立的四元数描述的线性形式的姿态动力学方程、体坐标系下的姿态动力学方程和刚体卫星的动力学方程,设计得到卫星姿控系统反馈控制律;
本步骤具体包括:
步骤2.1,假定当前四元数为q,期望四元数为qc,则误差四元数qe的定义为:
假定当前角速度为ω,期望角速度为ωc;则误差角速度ωe的定义为:
ωe=ω-ωc (2)
步骤2.2,基于误差四元数和误差角速度的姿态动力学模型可被描述为
其中:
I—整星的惯量矩阵;
—误差角加速度;
—期望角加速度;
H—整星的角动量,H=IBωB+Iwωw;
IB—卫星的转动惯量;
u—反作用飞轮的控制力矩,u=Iwωw;
Md—外部干扰力矩;
ME—地磁场引起的干扰力矩;
—误差四元数的一阶导数;
步骤2.3,将地球磁场近似为一个以地球为中心的偶极子;非线性控制系统的非线性部分能够利用反馈线性化方法被抵消,从而提出反馈控制律:
其中:
P,D—正定增益矩阵,通常被设计为对角矩阵;
qe123—误差四元数中的三个分量qe1、qe2、qe3组成的列向量,qe123=[qe1 qe2 qe3]T;
步骤3:根据步骤2提出的反馈控制律,对电磁编队卫星姿态和轨道进行相对控制。
2.根据权利要求1所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,步骤3具体为:
在根据步骤2提出的反馈控制律对电磁编队卫星姿态和轨道进行相对控制时,在控制过程中,实时判断是否发生磁饱和现象,如果发生,采用飞轮角动量磁卸载算法卸载磁力矩。
3.根据权利要求2所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,在控制过程中,实时判断是否发生磁饱和现象,如果发生,采用飞轮角动量磁卸载算法卸载磁力矩,具体为:
步骤3.1,假设反作用飞轮的冗余转动惯量为ΔH;在控制过程中,实时测量飞轮所处磁场B,判断是否满足以下的磁卸载条件:
45°<(B,ΔH)<135° (30)
如果满足,进一步判断磁场B是否垂直于ΔH,如果不垂直,则按式(37)计算得到磁卸载力矩值Tunloading,并将该磁卸载力矩值Tunloading施加于飞轮中,使其反向转动,进行磁卸载:
|B|—飞轮所处磁场B的大小;
K—增益系数;
如果B垂直于ΔH,则按下式(38)计算得到磁卸载力矩值Tunloading,并将该磁卸载力矩值Tunloading施加于飞轮中,使其反向转动,进行磁卸载:
Tunloading=-KΔH(K>0) (38) 。
4.根据权利要求1所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,还包括:
步骤4,设计大角度姿态机动方法,包括:
步骤4.1,假设机动前,即初始时刻t0时刻的理论姿态为q(t0),理论角速度为ω(t0)=[0 0 0]T;机动后,即tf时刻目标姿态为q(tf),目标角速度为ω(tf)=[0 0 0]T;
步骤4.2,根据下式计算得到控制力矩Lc:
其中:
IT—卫星及飞轮整体的转动惯量;
q123—q1、q2、q3组成的列向量,q123=[q1 q2 q3]T;
Kd—大角度姿态机动控制律的微分系数;
Kp—大角度姿态机动控制律的比例系数;
步骤4.3,将控制力矩Lc作用于飞轮,使其以最快速度按下式规划的路径进行姿态机动:
其中:
ωc(t)—期望角速度关于时间的函数;
—期望角加速度关于时间的函数;
qc(t)—期望四元数关于时间的函数;
—瞬时欧拉轴;
—绕欧拉轴旋转角加速度关于时间的函数;
—绕欧拉轴旋转角速度关于时间的函数;
θ(t)—绕欧拉轴旋转角度大小关于时间的函数。
5.根据权利要求1所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,还包括:
步骤5,建立相对运动动力学方程;
相对运动动力学方程为:
其中:
ρ—环绕星相对参考星的相对位置矢量;
—环绕星相对参考星的相对速度矢量;
—环绕星相对参考星的相对加速度矢量;
ωn—参考卫星的轨道角速度;
r—参考星在地惯系中的位置矢量;
um—电磁控制相对加速度;
R—干扰项和未建模的动力学部分;
C(ωn)∈R3×3为
R3×3—三维实数集;
非线性矩阵N(·)∈R3×3定义为如下等式:
μE—地球引力常数;
x、y、z—环绕星与参考星的相对位置矢量在Hill坐标系下的分量;
步骤6,根据步骤5建立的相对运动动力学方程,设计得到电磁编队飞行卫星运动控制器;
所述电磁编队飞行卫星运动控制器是一种终端滑模控制器,用于实现轨道跟踪误差保持在零附近。
6.根据权利要求5所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,轨道跟踪误差可被定义为
其中:
e—轨道跟踪位置误差;
—轨道跟踪速度误差;
ρd—目标相对位置矢量;
—目标相对速度矢量;
因此,设计的终端滑模控制器描述为:
其中:β>0,ki>0,i=1,2,3,0<γ<1;
us—终端滑模控制相对加速度;
—目标相对加速度矢量;
β—控制系数β>0;
γ—控制系数0<γ<1;
k—控制系数ki>0,i=1,2,3;
sat(·)—饱和函数,
φ0—一严格正常数,表示所设计的边界层的厚度;
S—滑模面,滑模面可描述为
|S|—滑模面的模。
7.根据权利要求6所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,为了消除振动,用连续的饱和函数代替了公式(51)中的饱和函数,连续的饱和函数可描述为
(53)
其中,设计系数τ满足0<τ<1;公式(51)可被写为
由此得到终端滑模控制器的另一种描述。
8.根据权利要求6所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,还包括:
步骤7,建立电磁力和电磁力矩的模型,包括:
步骤7.1,卫星j对卫星i的电磁力Fij和电磁力矩Tij可被描述为
其中:
μ0—真空磁导率,μ0=4π×10-7H/m;
rij—由卫星i到卫星j的相对位置矢量;
|rij|—由卫星i到卫星j的相对位置矢量的模;
μi—卫星i的磁偶极子的强度,μi=NiIiSi;
μj—卫星j的磁偶极子的强度,μj=NjIjSj;
Ni、Nj分别为偶极子A、B的绕线匝数;
Ii、Ij分别为偶极子A、B的线圈电流;
Si、Sj分别为偶极子A、B的线圈面积;
由公式(55)和公式(56),可得到如下三自由度的远场电磁模型:
其中:
—偶极子A受到的电磁力在三个轴方向的分量;
—偶极子A受到的电磁力矩在三个轴方向的分量;
μA—卫星A的磁偶极子强度;
μB—卫星B的磁偶极子强度;
α、β分别为偶极子A、B与ox轴之间的夹角;
χ、δ分别为偶极子A、B绕ox轴的转角;
步骤8,在步骤6设计的电磁编队飞行卫星运动控制器的基础上,基于步骤7的电磁力和电磁力矩的模型,为获取卫星间相对位置的控制过程中每颗卫星的磁偶极子,还设计两种计算磁偶极子的方法:
第一种是根据牛顿第三定律,将电磁编队飞行卫星的某一偶极子固定,计算可得其他偶极子的电磁力;
第二种是将偶极子的计算问题转换为一个优化问题,利用一种粒子群算法和序列二次规划方法混合算法以实现自由偶极子的最优化,最终实现电磁干扰力矩MT和地磁场干扰力矩ME最小的优化目标。
9.根据权利要求8所述的电磁编队卫星姿态和轨道相对控制方法,其特征在于,步骤7中,第一种计算磁偶极子的方法包括:
根据牛顿第三定律,电磁编队飞行卫星运动控制器所对应的公式(54)代表针对n个卫星的3(n-1)个标量公式,因此,系统在此约束下方程会有无穷多解,如果将电磁编队飞行卫星的某一偶极子固定,其他偶极子就可以通过如下公式进行计算
其中:
μ0—真空磁导率,μ0=4π×10-7H/m;
rij—由卫星i到卫星j的相对位置矢量;
|rij|—由卫星i到卫星j的相对位置矢量的模;
μi—卫星i的磁偶极子强度;
μj—卫星j的磁偶极子强度;
步骤7中,第二种计算磁偶极子的方法为:
建立以下目标函数:
μoptimal—最优磁矩;
arg(·)—幅值;
n—卫星个数;
MT—电磁干扰力矩;
ME—地磁场干扰力矩;
采用粒子群算法和序列二次规划方法混合算法以实现自由偶极子的最优化,其中,首先采用粒子群算法设定初始参数值,然后采用粒子群算法进行渐进估算,当收敛值满足预先设定的临界值时,将收敛值提供给序列二次规划方法,采用序列二次规划方法进行非线性约束下问题的非线性解算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610561867.3A CN106094854B (zh) | 2016-07-15 | 2016-07-15 | 电磁编队卫星姿态和轨道相对控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610561867.3A CN106094854B (zh) | 2016-07-15 | 2016-07-15 | 电磁编队卫星姿态和轨道相对控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106094854A CN106094854A (zh) | 2016-11-09 |
CN106094854B true CN106094854B (zh) | 2019-03-01 |
Family
ID=57221427
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610561867.3A Active CN106094854B (zh) | 2016-07-15 | 2016-07-15 | 电磁编队卫星姿态和轨道相对控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106094854B (zh) |
Families Citing this family (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107422741B (zh) * | 2017-04-28 | 2019-12-27 | 西北工业大学 | 基于学习的保预设性能集群飞行分布式姿态追踪控制方法 |
CN107065565B (zh) * | 2017-05-18 | 2019-03-01 | 中国人民解放军战略支援部队航天工程大学 | 一种用于集群航天器电磁拖拽的自抗扰控制方法 |
CN107239036B (zh) * | 2017-06-27 | 2019-10-01 | 上海航天控制技术研究所 | 一种近地卫星冗余飞轮角动量自主管理方法 |
CN107168372B (zh) * | 2017-06-30 | 2019-05-10 | 清华大学 | 基于在轨参数辨识和偏置的卫星跟飞长期摄动补偿方法 |
CN108181925A (zh) * | 2017-12-22 | 2018-06-19 | 西北工业大学 | 一种基于星间洛伦兹力的多卫星编队构型设计方法 |
CN108303990A (zh) * | 2017-12-26 | 2018-07-20 | 中国人民解放军战略支援部队航天工程大学 | 一种基于电磁力的集群航天器控制力/力矩共享方法 |
CN108279699B (zh) * | 2018-01-02 | 2019-11-08 | 东南大学 | 一种时空可变流场下飞行器的球面轨道编队跟踪控制方法 |
CN108363840B (zh) * | 2018-01-23 | 2022-03-01 | 中国人民解放军战略支援部队航天工程大学 | 一种基于电磁力的集群航天器磁矩最优分配方法 |
CN108762282B (zh) * | 2018-05-02 | 2023-07-04 | 中国人民解放军战略支援部队航天工程大学 | 一种基于电磁力的多集群航天器构形重构控制方法 |
CN108628345B (zh) * | 2018-06-28 | 2020-03-27 | 北京航空航天大学 | 一种电磁航天器编队悬停协同控制方法及系统 |
CN108958277A (zh) * | 2018-08-02 | 2018-12-07 | 西北工业大学 | 一种mems固体微推进器的微纳卫星姿轨一体化控制方法 |
JP2020032874A (ja) * | 2018-08-30 | 2020-03-05 | 国立研究開発法人宇宙航空研究開発機構 | 人工衛星の位置・姿勢制御システム及び人工衛星の位置・姿勢制御方式 |
CN110083170A (zh) * | 2019-04-11 | 2019-08-02 | 北京航空航天大学 | 一种使用固体微推力器进行轨道保持的优化控制方法 |
CN110231828B (zh) * | 2019-05-31 | 2020-07-21 | 燕山大学 | 基于非奇异快速终端滑模的四旋翼无人机视觉伺服控制方法 |
CN110316402A (zh) * | 2019-06-03 | 2019-10-11 | 上海航天控制技术研究所 | 一种编队控制模式下的卫星姿态控制方法 |
CN110231830B (zh) * | 2019-06-26 | 2022-05-24 | 北京控制工程研究所 | 一种有利于热平衡的航天器对日定向目标坐标系确定方法 |
CN110389592B (zh) * | 2019-07-18 | 2020-12-01 | 南京航空航天大学 | 一种基于分布式智能飞轮的航天器姿态控制方法 |
CN112319857B (zh) * | 2020-10-12 | 2022-09-13 | 中山大学 | 一种用于远距离分布式卫星的组合姿态控制方法及系统 |
CN112455726B (zh) * | 2020-12-10 | 2022-04-19 | 长光卫星技术股份有限公司 | 一种低轨光学遥感卫星多点成像快速机动控制方法 |
CN113071713B (zh) * | 2021-03-11 | 2022-11-22 | 中国空间技术研究院 | 卫星磁矩的分配方法和装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008047493A1 (en) * | 2006-10-18 | 2008-04-24 | National University Corporation NARA Institute of Science and Technology | Remote control system |
CN101794154A (zh) * | 2009-11-25 | 2010-08-04 | 哈尔滨工业大学 | 一种编队卫星相对轨道与姿态解耦控制方法 |
CN103092208A (zh) * | 2013-01-09 | 2013-05-08 | 哈尔滨工业大学 | 基于sgcmg和rw的航天器高精度快速姿态机动方法 |
CN103412571A (zh) * | 2013-07-22 | 2013-11-27 | 西北工业大学 | 一种基于反馈线性化技术的航天器相对姿态控制方法 |
CN104391446A (zh) * | 2014-11-05 | 2015-03-04 | 北京航空航天大学 | 一种平流层卫星六自由度动力学建模与级联控制方法 |
CN104406598A (zh) * | 2014-12-11 | 2015-03-11 | 南京航空航天大学 | 一种基于虚拟滑模控制的非合作航天器姿态估计方法 |
CN104635495A (zh) * | 2015-02-09 | 2015-05-20 | 北京航空航天大学 | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 |
CN104881036A (zh) * | 2015-06-19 | 2015-09-02 | 哈尔滨工业大学 | 基于代数Lyapunov方程的控制受限小卫星三轴磁力矩姿态控制方法 |
CN105468007A (zh) * | 2015-12-02 | 2016-04-06 | 哈尔滨工业大学 | 一种基于干扰观测器的挠性卫星轨迹线性化姿态控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080004758A1 (en) * | 2006-06-29 | 2008-01-03 | Northrop Grumman Corporation | Apparatus and method for tracking an orbital body relative to a planetary body using a single sensor |
US8000381B2 (en) * | 2007-02-27 | 2011-08-16 | Hemisphere Gps Llc | Unbiased code phase discriminator |
-
2016
- 2016-07-15 CN CN201610561867.3A patent/CN106094854B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2008047493A1 (en) * | 2006-10-18 | 2008-04-24 | National University Corporation NARA Institute of Science and Technology | Remote control system |
CN101794154A (zh) * | 2009-11-25 | 2010-08-04 | 哈尔滨工业大学 | 一种编队卫星相对轨道与姿态解耦控制方法 |
CN103092208A (zh) * | 2013-01-09 | 2013-05-08 | 哈尔滨工业大学 | 基于sgcmg和rw的航天器高精度快速姿态机动方法 |
CN103412571A (zh) * | 2013-07-22 | 2013-11-27 | 西北工业大学 | 一种基于反馈线性化技术的航天器相对姿态控制方法 |
CN104391446A (zh) * | 2014-11-05 | 2015-03-04 | 北京航空航天大学 | 一种平流层卫星六自由度动力学建模与级联控制方法 |
CN104406598A (zh) * | 2014-12-11 | 2015-03-11 | 南京航空航天大学 | 一种基于虚拟滑模控制的非合作航天器姿态估计方法 |
CN104635495A (zh) * | 2015-02-09 | 2015-05-20 | 北京航空航天大学 | 一种含模型参数不确定性平流层卫星的经向偏移分级控制方法 |
CN104881036A (zh) * | 2015-06-19 | 2015-09-02 | 哈尔滨工业大学 | 基于代数Lyapunov方程的控制受限小卫星三轴磁力矩姿态控制方法 |
CN105468007A (zh) * | 2015-12-02 | 2016-04-06 | 哈尔滨工业大学 | 一种基于干扰观测器的挠性卫星轨迹线性化姿态控制方法 |
Non-Patent Citations (2)
Title |
---|
基于Simulink/State的小卫星姿态控制模块设计与仿真研究;吕波 等;《科学技术与工程》;20110831;第11卷(第23期);第5609-5611页 |
基于处理器在回路的卫星编队仿真验证系统;胡敏 等;《系统仿真学报》;20101130;第22卷(第11期);第2618-2620页 |
Also Published As
Publication number | Publication date |
---|---|
CN106094854A (zh) | 2016-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106094854B (zh) | 电磁编队卫星姿态和轨道相对控制方法 | |
Bayat | Model predictive sliding control for finite-time three-axis spacecraft attitude tracking | |
US8131409B2 (en) | Gyroless transfer orbit sun acquisition using only wing current measurement feedback | |
JP4783366B2 (ja) | Cmgアレイにおけるトルクを最適化するための方法及びシステム | |
US20080315039A1 (en) | System and methods for space vehicle torque balancing | |
Hu et al. | Attitude control and vibration suppression for flexible spacecraft using control moment gyroscopes | |
Hu et al. | Dynamics and vibration suppression of space structures with control moment gyroscopes | |
CN103488081B (zh) | 惯性稳定平台控制方法 | |
Kök | Comparison and analysis of attitude control systems of a satellite using reaction wheel actuators | |
Hallaj et al. | Sliding mode control of electromagnetic tethered satellite formation | |
Liu et al. | Dynamics and control of a flexible solar sail | |
Ye et al. | Fuzzy sliding mode control of nonparallel-ground-track imaging satellite with high precision | |
Kojima et al. | Steering control law for double-gimbal scissored-pair CMG | |
Guo et al. | Integrated vibration isolation and attitude control for spacecraft with uncertain or unknown payload inertia parameters | |
Blonk | Modeling and control of a ball-balancing robot | |
Fan et al. | Rotational and relative translational control for satellite electromagnetic formation flying in low earth orbit | |
Xu et al. | PD-type magnetic-based sun-pointing attitude control scheme | |
Kuchwa-Dube et al. | Altitude and attitude tracking of a quadrotor-based aerial manipulator using super twisting sliding mode control | |
Bayat et al. | A heuristic design method for attitude stabilization of magnetic actuated satellites | |
Zhang et al. | Fixed-time sliding mode attitude control of a flexible spacecraft with rotating appendages connected by magnetic bearing | |
Li et al. | Modeling and Analysis of a 3-DOF spherical momentum exchange actuator for spacecraft attitude maneuver | |
Sun et al. | Optimal spacecraft rendezvous using the combination of aerodynamic force and Lorentz force | |
Reijneveld et al. | Attitude control system of the Delfi-n3Xt satellite | |
Lavezzi et al. | Hybrid Open/Closed-Loop Attitude Control Method for Imaging Satellites | |
Carabellese et al. | Magnetorquer-only nonlinear attitude control for cubesats |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: No. 1, Bayi Road, Huairou District, Beijing City, Beijing Applicant after: China People's Liberation Army Strategic Support Unit Space Engineering University Address before: No. 1, Bayi Road, Huairou District, Beijing City, Beijing Applicant before: PLA 'S Equipment College |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |