CN103991559B - 一种洛伦兹航天器悬停控制方法 - Google Patents
一种洛伦兹航天器悬停控制方法 Download PDFInfo
- Publication number
- CN103991559B CN103991559B CN201410230654.3A CN201410230654A CN103991559B CN 103991559 B CN103991559 B CN 103991559B CN 201410230654 A CN201410230654 A CN 201410230654A CN 103991559 B CN103991559 B CN 103991559B
- Authority
- CN
- China
- Prior art keywords
- lorentz
- hovering
- omega
- centerdot
- spacecraft
- 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 27
- 230000001133 acceleration Effects 0.000 claims abstract description 120
- 239000000446 fuel Substances 0.000 claims abstract description 10
- 230000005358 geomagnetic field Effects 0.000 claims abstract description 10
- 238000004364 calculation method Methods 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 239000000725 suspension Substances 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 2
- 238000009795 derivation Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 239000011159 matrix material Substances 0.000 claims description 2
- 230000017105 transposition Effects 0.000 claims description 2
- 101100042630 Caenorhabditis elegans sin-3 gene Proteins 0.000 claims 1
- 239000000126 substance Substances 0.000 abstract description 3
- 230000005669 field effect Effects 0.000 abstract 1
- 230000001141 propulsive effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 5
- 230000007774 longterm Effects 0.000 description 4
- 230000010354 integration Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000005389 magnetism Effects 0.000 description 1
- 238000005312 nonlinear dynamic Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
Landscapes
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
传统航天器轨道控制一般采用化学推进方式,因而任务周期受所携推进工质质量制约,尤其是对能量需求较高的低轨悬停控制,所需能量之多限制了采用传统推进方式对低轨空间目标长期悬停的可实现性。本发明提出了一种采用带电洛伦兹航天器在地磁场中运动时受到的洛伦兹力与传统推力器推力联合控制的燃料最优悬停控制方法,创新性地解决了这一问题,将悬停可实现范围延伸至低轨。本发明基于建立的地球非球形摄动作用下的洛伦兹航天器相对运动模型,给出了燃料最优条件下实现悬停所需的洛伦兹航天器最优荷质比与推力器提供的最优控制加速度。
Description
技术领域
本发明涉及航天器悬停控制方法,特别涉及一种采用洛伦兹力与传统推力联合进行悬停控制的方法。
背景技术
所谓航天器悬停,是指通过力和力矩控制,保持航天器相对于某空间目标的位置始终保持不变,航天器相对于空间目标仿佛是静止“悬停”于某个固定点上。悬停轨道在空间任务中的应用广泛。若将航天器悬停视为一种特殊形式的空间编队飞行,则最简单的悬停构型为跟飞构型,该构型中,追踪器与目标器位于同一圆轨道,仅提前或落后目标一固定相位角。除该特殊悬停构型,为实现悬停,需对航天器持续施加控制加速度以保持其对目标的相对位置不变。因而,悬停轨道多为非开普勒轨道。传统航天器轨道控制一般采用化学推进方式。因而,空间任务周期往往受航天器所携推进工质质量制约,尤其是对能量需求较高的低轨悬停控制。例如,对500km赤道圆轨道的目标航天器在其径向上方5km悬停一天所需的速度增量为1.586km/s。所需能量之多限制了采用传统推进方式对低轨空间目标长期悬停的可实现性,成为制约低轨航天器执行悬停任务的技术难点。
空间带电洛伦兹航天器在地磁场中运动时通常受到洛伦兹力作用,该洛伦兹力可用以进行航天器绝对轨道与相对轨道控制,从而实现无推进工质消耗的轨道机动。洛伦兹力由当地磁场以及航天器与当地磁场相对速度决定。由于低轨(轨道高度:100km至2000km)磁场强度明显强于高轨(轨道高度:20000km以上)且相对速度更快,因而,洛伦兹航天器在低轨应用效率更高。采用洛伦兹力辅助悬停将有效减少燃料消耗,使得长期悬停任务的可实现范围延伸至低轨。
发明内容
本发明针对洛伦兹力辅助悬停控制问题,假设地磁场为一随地球自转的倾斜磁偶极子且洛伦兹航天器可视为带电质点,采用拉格朗日动力学方法建立了地球非球形摄动二阶带谐系数(以下简称J2)摄动作用下的洛伦兹力辅助悬停动力学模型,基于所建立的悬停动力学模型,分析给出了实现悬停所需的开环总控制加速度,由于洛伦兹力的方向始终与当地磁场方向以及航天器相对当地磁场速度方向垂直,这一特性决定了洛伦兹力方向并不一定与悬停所需开环总控制加速度方向相同。因此,需要采用推力器推力与洛伦兹力联合控制,其中洛伦兹力作为辅助推力,减少推进工质的消耗。如何确定洛伦兹力与推力器推力的最优配比从而使得悬停任务周期内实现悬停所需速度增量消耗最少,是本发明设计中的技术难点。本发明基于欧拉-拉格朗日方程提出了一种洛伦兹力与推力器推力联合控制的燃料最优混合动力悬停控制方法。
本发明的技术方案如下:
一种航天器悬停控制方法,按以下步骤进行:
步骤一:给定初始时刻目标轨道参数与地磁场参数,给定悬停构型;
步骤二:各轴悬停距离计算:计算追踪器与目标器的相对位置矢量在相对运动坐标系中的各轴向分量;
步骤三:悬停所需总控制加速度计算:推导J2摄动作用下洛伦兹力辅助悬停动力学方程,根据悬停定义,计算实现悬停所需的各轴向总控制加速度;
步骤四:悬停所需洛伦兹航天器最优荷质比计算:以任务历时内实现悬停所需的速度增量消耗最优为指标,求解欧拉-拉格朗日方程得到该指标下洛伦兹航天器最优荷质比轨迹;
步骤五:悬停所需最优洛伦兹加速度及推力器提供的最优控制加速度计算:将洛伦兹航天器最优荷质比代入洛伦兹加速度定义式得到最优洛伦兹加速度,通过力的分解得到实现悬停所需的由推力器提供的最优控制加速度。
由于在地球静止或地球同步轨道(轨道高度约为35786km、轨道周期等于地球自转周期的圆轨道,以下简称GEOs),洛伦兹航天器与地磁场不存在相对运动速度,因而无法产生洛伦兹力,也就无法对此类特殊轨道目标进行悬停,除去该类特殊轨道外,洛伦兹航天器均可适用;考虑到洛伦兹航天器在低轨应用效率高,且对低轨目标悬停所需速度增量大,优选的,本发明应用于对低轨空间目标的悬停,即轨道高度在100km至2000km;
在步骤一中所述的悬停构型Γ=[ρ,E,A],其中ρ为悬停距离,E与A分别为悬停高低角与悬停方位角,其定义方式如图1所示,图中,OEXIYIZI为地心惯性坐标系,OE为地心,相对运动在相对运动坐标系OTxyz中描述,该坐标系原点位于目标器质心OT,其中x轴沿目标器径向,z轴与目标器轨道面法向重合,y轴与x、z轴构成笛卡尔右手直角坐标系,OL为洛伦兹航天器质心,悬停高低角E的范围为E∈[-π/2,π/2],规定追踪器悬停于目标器上方时为正,反之为负;悬停方位角A的范围为A∈[-π,π],从飞行方向(+y)看,规定追踪器悬停于目标器左方时为正,反之为负,根据上述定义,图中的悬停高低角与方位角均为正;
在步骤二中所述的计算追踪器与目标器的相对位置矢量在相对运动坐标系中的各轴向分量,其计算方法为:
式中,ρ=[xyz]T为追踪器与目标器的相对位置矢量在相对运动坐标系中的描述,上标T为矩阵转置符号;
在步骤三中所述的实现悬停所需的各轴向总控制加速度,其计算方法如下:
(1)建立J2摄动作用下的洛伦兹航天器相对运动动力学模型:
基于拉格朗日动力学方法推导得到的J2摄动作用下的洛伦兹力辅助悬停动力学方程:
其中,
式中,洛伦兹航天器地心距为RL=[(RT+x)2+y2+z2]1/2;
ω=[ωxωyωz]T与ε=[εxεyεz]T分别为相对运动坐标系的角速度与角加速度,即
式中,RLz=(RT+x)siniTsinuT+ysiniTcosuT+zcosiT为洛伦兹航天器地心距在地惯系中ZI方向的投影,iT与uT分别为目标器的轨道倾角与纬度幅角,系数其中地球非球形摄动二阶带谐系数为J2=1.0826×10-3,μ=3.986005×1014m3/s2为地球引力常数,RE为地球半径;
aC=[aRaSaW]T与aL=λl=λ[lxlylz]T分别为作用于洛伦兹航天器的控制加速度与洛伦兹加速度,其中λ=q/m为洛伦兹航天器的荷质比,且
式中,B0=8.0×1015T·m3为地磁场地磁矩值。Vrel=[VxVyVz]T为洛伦兹航天器与当地磁场相对速度在相对运动坐标系中的描述:
式中,β=ΩM-ΩT,其中ΩM=ωEt+Ω0为地磁轴在地惯系赤道平面内的相位角,ωE为地球自转角速度,Ω0为初始时刻地磁轴在地惯系赤道平面内的相位角,ΩT为目标器轨道升交点赤经,倾斜地磁轴与地球自转轴的夹角为α,各角度定义及其与地惯系、相对运动坐标系位置关系如图2所示;
六个轨道参数(RT,Vr,hT,iT,uT,ΩT)分别为目标器的轨道半径、径向速度、轨道角动量、轨道倾角、近地点幅角以及升交点赤经,其满足动力学方程
式(2)至式(10)即J2摄动作用下的洛伦兹航天器非线性相对运动方程;
(2)计算悬停所需的开环控制总加速度:
根据悬停定义,在相对运动坐标系中,追踪器与目标器的相对位置矢量保持不变,因而相对位置各轴分量对时间的各阶导数为零,即
将式(11)代入式(2)中得,
aC=h-λl(12)
式中,h为实现悬停所需的总控制加速度,即
在步骤四中所述的实现悬停所需的最优洛伦兹航天器荷质比,其计算方法如下:
由于洛伦兹加速度方向与当地磁场方向及洛伦兹航天器相对磁场速度方向垂直,当洛伦兹航天器轨道位置确定,洛伦兹加速度方向也随之确定;若洛伦兹加速度方向与悬停所需加速度方向重合,即l∥h,则若洛伦兹航天器荷质比大小为λ=||h||/||l||,产生的洛伦兹加速度即可完全抵偿悬停所需加速度;然而,由于洛伦兹力作用方向的局限性,并不能保证洛伦兹加速度方向与悬停所需加速度方向重合,因此,为实现悬停,需施加额外的控制加速度,本发明中假设该额外控制加速度由推力器提供,为最大限度地节省燃料,选取拉格朗日形式的燃料最优目标函数,即:
式中,tf为悬停任务历时;
求解欧拉-拉格朗日方程:
得到实现燃料最省的洛伦兹航天器荷质比最优轨迹
当磁场B方向与相对速度Vrel方向平行时,无洛伦兹加速度产生,即||l||=0,因而,悬停所需加速度将完全由推力器提供,该条件下设定洛伦兹航天器荷质比为0,如式(16)所示;
在步骤五中所述的实现悬停所需最优洛伦兹加速度及推力器提供的最优控制加速度,其计算方法如下:
将洛伦兹航天器最优荷质比λ*代入洛伦兹加速度定义式,得到实现悬停所需的最优洛伦兹加速度为
将式(17)代入式(12),得到最优控制加速度
因此,对于带电洛伦兹航天器,实现悬停所需的速度增量为
而对于非带电航天器,实现悬停所需的控制加速度均由推力器提供,因而所需速度增量为
与现有技术相比,本发明具有以下优点:
1)本发明直接基于洛伦兹航天器相对运动非线性动力学模型设计,考虑了航天器相对运动的非线性特性,克服了传统相对运动模型由线性化引起的误差。基此模型求解得到的悬停控制加速度更为精确。
2)本发明考虑到J2摄动是低地球轨道的主要摄动力且悬停构型保持对摄动力较为敏感,建立了J2摄动力与洛伦兹力共同作用下的洛伦兹航天器相对运动模型,更符合低地球轨道空间环境特性,因而求解得到的悬停控制加速度更为精确。
3)本发明采用洛伦兹力与推力器推力联合进行航天器悬停控制,创新性地解决了采用传统化学推进方式由于推进工质质量制约而无法对低轨道空间目标长期悬停的问题,利用洛伦兹力作为辅助推力,节省了速度增量消耗,延长了悬停任务周期,将长期悬停任务可实现范围延伸至低地球轨道。
控制工程师在应用过程中可以根据实际空间目标给定任意悬停构型,并将由该方法得到的控制量传输至执行机构实现对空间目标的悬停。
附图说明
图1为本发明坐标系与悬停构型定义图;
图2为本发明倾斜偶极子空间指向与相关角度定义图;
图3为本发明所述控制方法流程图;
图4为本发明实现悬停一天所需总控制加速度图;
图5为本发明实现悬停一个周期所需总控制加速度图;
图6为本发明实现悬停一天所需洛伦兹航天器最优荷质比图;
图7为本发明实现悬停一个周期所需洛伦兹航天器最优荷质比图;
图8为本发明实现悬停一天所需最优洛伦兹加速度图;
图9为本发明实现悬停一个周期所需最优洛伦兹加速度图;
图10为本发明实现悬停一天所需最优控制加速度对比图;
图11为本发明实现悬停一个周期所需最优控制加速度对比图;
图中及本发明中符号说明如下:
A悬停方位角,deg或rad
aC控制加速度矢量,m/s2
最优控制加速度矢量,m/s2
||aC||控制加速度大小,m/s2
aL洛伦兹加速度矢量,m/s2
最优洛伦兹加速度矢量,m/s2
||aL||洛伦兹加速度大小,m/s2
aR控制加速度径向分量,m/s2
aS控制加速度迹向分量,m/s2
aW控制加速度法向分量,m/s2
ax洛伦兹加速度径向分量,m/s2
ay洛伦兹加速度迹向分量,m/s2
az洛伦兹加速度法向分量,m/s2
B磁场强度矢量,T
Bx磁场强度径向分量,T
By磁场强度迹向分量,T
Bz磁场强度法向分量,T
B0地磁场地磁矩值,T·m3
E悬停高低角,deg或rad
h悬停所需总加速度矢量,m/s2
||h||悬停所需总加速度大小,m/s2
hT目标航天器轨道角动量,m2/s
hx悬停所需总加速度径向分量,m/s2
hy悬停所需总加速度迹向分量,m/s2
hz悬停所需总加速度法向分量,m/s2
iT目标航天器轨道倾角,deg或rad
J指标函数,m/s
J2地球非球形摄动二阶带谐系数
L拉格朗日函数
l洛伦兹加速度方向矢量(非单位矢量),(m·T)/s
lx洛伦兹加速度方向矢量径向分量,(m·T)/s
ly洛伦兹加速度方向矢量迹向分量,(m·T)/s
lz洛伦兹加速度方向矢量法向分量,(m·T)/s
m洛伦兹航天器质量,kg
单位偶极子矢量
nx单位偶极子矢量径向分量
ny单位偶极子矢量迹向分量
nz单位偶极子矢量法向分量
OEXIYIZI地心惯性坐标系
OL洛伦兹航天器质心
OT目标航天器质心
OTxyz相对运动坐标系
q洛伦兹航天器带电量,C
q/m洛伦兹航天器荷质比,C/kg
RE地球半径,m
RL洛伦兹航天器地心距,m
RLz洛伦兹航天器地心距在地心惯性坐标系OEXIYIZI中OEZI轴上的投影,m
RT目标航天器地心距,m
t时间,s
tf悬停任务历时,s
uT目标航天器纬度幅角,deg或rad
Vr目标航天器径向速度,m/s
Vrel洛伦兹航天器与当地磁场相对速度矢量,m/s
Vx洛伦兹航天器与当地磁场相对速度径向分量,m/s
Vy洛伦兹航天器与当地磁场相对速度迹向分量,m/s
Vz洛伦兹航天器与当地磁场相对速度法向分量,m/s
x追踪器与目标器径向相对距离,m
y追踪器与目标器迹向相对距离,m
z追踪器与目标器法向相对距离,m
α地磁轴倾角,deg或rad
β地惯系赤道面内地磁轴相位角与目标器轨道升交点赤经之差,deg或rad
ΔV实现悬停所需速度增量,m/s
ε相对运动坐标系转动角加速度矢量,rad/s2
εx相对运动坐标系转动角加速度矢量径向分量,rad/s2
εy相对运动坐标系转动角加速度矢量迹向分量,rad/s2
εz相对运动坐标系转动角加速度矢量法向分量,rad/s2
η系数
ηL系数
Γ悬停构型
λ洛伦兹航天器荷质比,C/kg
λ*最优洛伦兹航天器荷质比,C/kg
μ地球引力常数,m3/s2
ρ追踪器与目标器相对位置矢量,m
ρ追踪器与目标器相对距离,m
ΩM地磁轴在地惯系赤道面内相位角,deg或rad
ΩT目标航天器轨道升交点赤经,deg或rad
Ω0初始时刻地磁轴在地惯系赤道面内相位角,deg或rad
ω相对运动坐标系转动角速度矢量,rad/s
ωE地球自转角速度,rad/s
ωx相对运动坐标系转动角速度矢量径向分量,rad/s
ωy相对运动坐标系转动角速度矢量迹向分量,rad/s
ωz相对运动坐标系转动角速度矢量法向分量,rad/s
ξ系数
ξL系数
具体实施方案
下面结合附图,对本发明中的设计方法作进一步的说明:
本发明“一种航天器悬停控制方法”,其具体步骤如下:
步骤一:给定初始时刻目标轨道参数、地磁场参数与悬停构型:
如图3中所示,给定初始时刻目标航天器轨道参数如表1所示,目标轨道周期约为1.6h,属于低轨近圆轨道,初始时刻地磁轴在地惯系赤道面内的相位角Ω0=-60°,悬停构型为Γ=[4km,π/2,0],其中悬停距离ρ为4km,悬停高低角E为90°,悬停方位角A为0°;
表1初始时刻目标航天器轨道根数
步骤二:各轴悬停距离计算
如图3中所示,计算悬停距离在相对运动坐标系各轴向分量。将悬停构型代入式(1),得到各轴向的悬停距离分别为:
步骤三:悬停所需总控制加速度计算
如图3中所示,计算实现悬停所需的各轴向总控制加速度。将步骤一与步骤二中的参数代入式(13),得到悬停任务过程中任一时刻所需的总控制加速度h,其中,目标航天器的六个轨道参数(RT,Vr,hT,iT,uT,ΩT)为对式(10)进行数值积分的结果,一天内实现该悬停构型所需的各轴向总控制加速度如图4所示,一个轨道轨道周期内实现该悬停构型所需的各轴向总控制加速度如图5所示;
步骤四:悬停所需洛伦兹航天器最优荷质比计算
如图3中所示,计算燃料最优指标条件下实现悬停所需的洛伦兹航天器最优荷质比。将步骤一与步骤二中的参数代入式(7)中,得到悬停任务过程总任一时刻洛伦兹加速度的方向矢量。将该方向矢量l与悬停所需总加速度矢量h代入式(16),即得到悬停任务过程中任一时刻实现悬停所需的洛伦兹航天器最优荷质比,如图6与图7所示所示,其中,图6为一天内的最优洛伦兹航天器荷质比,图7为一个周期内的最优洛伦兹航天器荷质比;
步骤五:悬停所需最优洛伦兹加速度及推力器提供的最优控制加速度计算
如图3中所示,计算燃料最优指标条件下实现悬停所需的最优洛伦兹加速度以及由推力器提供的最优控制加速度;将步骤四中得到的洛伦兹航天器最优荷质比λ*与洛伦兹加速度方向矢量l代入式(17),即得到实现悬停所需的最优洛伦兹加速度,如图8和图9所示,其中,图8为一天内的最优洛伦兹加速度,图9为一个周期内的最优洛伦兹加速度,同理,将步骤三中得到的悬停所需总加速度矢量h与步骤四中得到的洛伦兹加速度方向矢量l代入式(18),即得到采用洛伦兹航天器实现悬停所需的由推力器提供的最优控制加速度如图10和图11所示,图10中给出了一天内采用传统非带电航天器与带电洛伦兹航天器实现悬停所需的由推力器提供的控制加速度对比图,同理,图11为一个周期内的对比图。
将实现悬停所需的总控制加速度h轨迹代入式(20)进行数值积分,得到非带电航天器悬停所需的速度增量,悬停一个周期所需的增量为82.71m/s,悬停一天(按24h计)所需的速度增量为1.24×103m/s,将洛伦兹航天器实现悬停所需的控制加速度轨迹代入式(19)进行数值积分,得到采用洛伦兹力辅助悬停所需的速度增量,采用洛伦兹力辅助悬停,一个周期内所需速度增量为4.54m/s,节省了约94.51%速度增量,一天内所需速度增量为365.85m/s,节省了约70.51%速度增量,对比发现,采用洛伦兹力辅助悬停,可有效减少速度增量消耗,延长悬停任务历时,将长期悬停可实现范围延伸至低地球轨道。
Claims (2)
1.一种航天器悬停控制方法,应用于非地球静止或地球同步轨道空间目标的悬停,按以下步骤进行:
步骤一:给定初始时刻目标轨道参数与地磁场参数,给定悬停构型;
所述的悬停构型Γ=[ρ,E,A],其中ρ为悬停距离,E与A分别为悬停高低角与悬停方位角,OEXIYIZI为地心惯性坐标系,OE为地心,相对运动在相对运动坐标系OTxyz中描述,该坐标系原点位于目标器质心OT,其中x轴沿目标器径向,z轴与目标器轨道面法向重合,y轴与x、z轴构成笛卡尔右手直角坐标系,OL为洛伦兹航天器质心,悬停高低角E的范围为E∈[-π/2,π/2],规定追踪器悬停于目标器上方时为正,反之为负;悬停方位角A的范围为A∈[-π,π],从飞行方向(+y)看,规定追踪器悬停于目标器左方时为正,反之为负;
步骤二:各轴悬停距离计算:计算追踪器与目标器的相对位置矢量在相对运动坐标系中的各轴向分量;
所述的计算追踪器与目标器的相对位置矢量在相对运动坐标系中的各轴向分量,其计算方法为:
式中,ρ=[xyz]T为追踪器与目标器的相对位置矢量在相对运动坐标系中的描述,上标T为矩阵转置符号;
步骤三:悬停所需总控制加速度计算:推导J2摄动作用下洛伦兹力辅助悬停动力学方程,根据悬停定义,计算实现悬停所需的各轴向总控制加速度,其中J2为地球非球形摄动二阶带谐系数;
所述的实现悬停所需的各轴向总控制加速度,其计算方法如下:
(1)建立J2摄动作用下的洛伦兹航天器相对运动动力学模型:
基于拉格朗日动力学方法推导得到的J2摄动作用下的洛伦兹力辅助悬停动力学方程:
其中,
式中,洛伦兹航天器地心距为RL=[(RT+x)2+y2+z2]1/2;
ω=[ωxωyωz]T与ε=[εxεyεz]T分别为相对运动坐标系的角速度与角加速度,即
式中,RLz=(RT+x)siniTsinuT+ysiniTcosuT+zcosiT为洛伦兹航天器地心距在地惯系中ZI方向的投影,iT与uT分别为目标器的轨道倾角与纬度幅角,系数其中地球非球形摄动二阶带谐系数为J2=1.0826×10-3,μ=3.986005×1014m3/s2为地球引力常数,RE为地球半径;
aC=[aRaSaW]T与aL=λl=λ[lxlylz]T分别为作用于洛伦兹航天器的控制加速度与洛伦兹加速度,其中λ=q/m为洛伦兹航天器的荷质比,且
式中,B0=8.0×1015T·m3为地磁场地磁矩值,Vrel=[VxVyVz]T为洛伦兹航天器与当地磁场相对速度在相对运动坐标系中的描述:
式中,β=ΩM-ΩT,其中ΩM=ωEt+Ω0为地磁轴在地惯系赤道平面内的相位角,ωE为地球自转角速度,Ω0为初始时刻地磁轴在地惯系赤道平面内的相位角,ΩT为目标器轨道升交点赤经,倾斜地磁轴与地球自转轴的夹角为α;
六个轨道参数(RT,Vr,hT,iT,uT,ΩT)分别为目标器的轨道半径、径向速度、轨道角动量、轨道倾角、近地点幅角以及升交点赤经,其满足动力学方程
式(2)至式(10)即J2摄动作用下的洛伦兹航天器非线性相对运动方程;
(2)计算悬停所需的开环控制总加速度:
根据悬停定义,在相对运动坐标系中,追踪器与目标器的相对位置矢量保持不变,因而相对位置各轴分量对时间的各阶导数为零,即
将式(11)代入式(2)中得,
aC=h-λl(12)
式中,h为实现悬停所需的总控制加速度,即
步骤四:悬停所需洛伦兹航天器最优荷质比计算:以任务历时内实现悬停所需的速度增量消耗最优为指标,求解欧拉-拉格朗日方程得到该指标下洛伦兹航天器最优荷质比轨迹;
在步骤四中所述的实现悬停所需的最优洛伦兹航天器荷质比,其计算方法如下:
由于洛伦兹加速度方向与当地磁场方向及洛伦兹航天器相对磁场速度方向垂直,当洛伦兹航天器轨道位置确定,洛伦兹加速度方向也随之确定;若洛伦兹加速度方向与悬停所需加速度方向重合,即l∥h,则若洛伦兹航天器荷质比大小为λ=||h||/||l||,产生的洛伦兹加速度即可完全抵偿悬停所需加速度;然而,由于洛伦兹力作用方向的局限性,并不能保证洛伦兹加速度方向与悬停所需加速度方向重合,因此,为实现悬停,需施加额外的控制加速度,该额外控制加速度由推力器提供,为最大限度地节省燃料,选取拉格朗日形式的燃料最优目标函数,即:
式中,tf为悬停任务历时;
求解欧拉-拉格朗日方程:
得到实现燃料最省的洛伦兹航天器荷质比最优轨迹
当磁场B方向与相对速度Vrel方向平行时,无洛伦兹加速度产生,即||l||=0,因而,悬停所需加速度将完全由推力器提供,该条件下设定洛伦兹航天器荷质比为0,如式(16)所示;
步骤五:悬停所需最优洛伦兹加速度及推力器提供的最优控制加速度计算:将洛伦兹航天器最优荷质比代入洛伦兹加速度定义式得到最优洛伦兹加速度,通过力的分解得到实现悬停所需的由推力器提供的最优控制加速度;
在步骤五中所述的实现悬停所需最优洛伦兹加速度及推力器提供的最优控制加速度,其计算方法如下:
将洛伦兹航天器最优荷质比λ*代入洛伦兹加速度定义式,得到实现悬停所需的最优洛伦兹加速度为
将式(17)代入式(12),得到最优控制加速度
因此,对于带电洛伦兹航天器,实现悬停所需的速度增量为
而对于非带电航天器,实现悬停所需的控制加速度均由推力器提供,因而所需速度增量为
2.权利要求1所述一种航天器悬停控制方法,其特征在于:该方法应用于对低轨空间目标的悬停,即轨道高度在100km至2000km。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410230654.3A CN103991559B (zh) | 2014-05-28 | 2014-05-28 | 一种洛伦兹航天器悬停控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410230654.3A CN103991559B (zh) | 2014-05-28 | 2014-05-28 | 一种洛伦兹航天器悬停控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103991559A CN103991559A (zh) | 2014-08-20 |
CN103991559B true CN103991559B (zh) | 2016-02-17 |
Family
ID=51305961
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410230654.3A Active CN103991559B (zh) | 2014-05-28 | 2014-05-28 | 一种洛伦兹航天器悬停控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103991559B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104554825B (zh) * | 2014-12-08 | 2016-04-27 | 南京航空航天大学 | 航天器磁推进装置 |
CN104787345A (zh) * | 2015-04-29 | 2015-07-22 | 田鸿波 | 一种电磁动力航天飞机 |
CN104932520B (zh) * | 2015-05-25 | 2018-01-09 | 北京航空航天大学 | 一种存在测量误差下的航天器雨滴形状悬停构型控制方法 |
CN105527974B (zh) * | 2016-01-31 | 2018-07-06 | 中国人民解放军国防科学技术大学 | 一种缺失径向控制的欠驱动航天器悬停渐近控制方法 |
CN107367942B (zh) * | 2017-08-31 | 2020-04-21 | 北京理工大学 | 一种小天体定点悬停常推力控制方法 |
CN108181925A (zh) * | 2017-12-22 | 2018-06-19 | 西北工业大学 | 一种基于星间洛伦兹力的多卫星编队构型设计方法 |
CN108280258A (zh) * | 2017-12-22 | 2018-07-13 | 西北工业大学 | 一种基于洛伦兹力的伴飞轨道设计方法 |
CN108052008B (zh) * | 2018-01-10 | 2020-04-07 | 大连理工大学 | 一种倾转旋翼飞行器过渡态切换过程的几何最优控制器设计方法 |
CN109240340B (zh) * | 2018-11-27 | 2020-04-24 | 北京理工大学 | 一种基于拟周期轨道的洛伦兹力多星编队构型方法 |
CN111446888B (zh) * | 2020-04-27 | 2021-07-30 | 北京理工大学 | 一种基于人工磁场的洛伦兹力悬浮方法 |
CN111552308B (zh) * | 2020-05-15 | 2021-07-02 | 清华大学 | 一种控制方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2815730A1 (fr) * | 2000-10-19 | 2002-04-26 | Astrium Sas | Procede de commande d'attitude et de stabilisation d'un satellite en orbite basse, par couplage avec le champ magnetique terrestre |
EP1308813A1 (fr) * | 2001-10-19 | 2003-05-07 | Astrium SAS | Procédé de commande d'attitude d'un satellite en orbite terrestre basse |
CN102759928A (zh) * | 2012-04-28 | 2012-10-31 | 中国人民解放军国防科学技术大学 | 一种平流层飞艇航迹控制方法 |
-
2014
- 2014-05-28 CN CN201410230654.3A patent/CN103991559B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2815730A1 (fr) * | 2000-10-19 | 2002-04-26 | Astrium Sas | Procede de commande d'attitude et de stabilisation d'un satellite en orbite basse, par couplage avec le champ magnetique terrestre |
EP1308813A1 (fr) * | 2001-10-19 | 2003-05-07 | Astrium SAS | Procédé de commande d'attitude d'un satellite en orbite terrestre basse |
CN102759928A (zh) * | 2012-04-28 | 2012-10-31 | 中国人民解放军国防科学技术大学 | 一种平流层飞艇航迹控制方法 |
Non-Patent Citations (2)
Title |
---|
卫星相对空间目标任意位置悬停的方法研究;闫野;《中国空间科学技术》;20090225(第1期);第1-3页 * |
黄煦等.赤道圆轨道洛伦兹航天器交会最优控制.《第三十二届中国控制会议论文集(B卷)》.2013, * |
Also Published As
Publication number | Publication date |
---|---|
CN103991559A (zh) | 2014-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103991559B (zh) | 一种洛伦兹航天器悬停控制方法 | |
Tsiotras et al. | Satellite attitude control and power tracking with energy/momentum wheels | |
CN104527994B (zh) | 异面交叉快变轨道固定时间稳定姿态指向跟踪控制方法 | |
CN101381004B (zh) | 基于大气阻力的微小卫星编队飞行控制方法及控制装置 | |
Hoffmann et al. | Precision flight control for a multi-vehicle quadrotor helicopter testbed | |
CN104035335B (zh) | 基于高精度纵、横程解析预测方法的平稳滑翔再入制导方法 | |
Takada et al. | Control moment gyro singularity-avoidance steering control based on singular-surface cost function | |
CN103076807B (zh) | 一种欠驱动柔性航天器姿态稳定的控制方法 | |
CN103092208A (zh) | 基于sgcmg和rw的航天器高精度快速姿态机动方法 | |
CN104309822B (zh) | 一种基于参数优化的航天器单脉冲水滴形绕飞轨迹悬停控制方法 | |
CN106096204B (zh) | 一种基于太阳帆推进技术的航天器日心椭圆悬浮轨道设计方法 | |
CN107487458A (zh) | 一种全电推进卫星平台姿轨控执行机构的*** | |
Chen et al. | Path planning for UAVs formation reconfiguration based on Dubins trajectory | |
CN110450991B (zh) | 微纳卫星集群捕获空间非合作目标的方法 | |
CN104249816A (zh) | 非合作目标绕飞悬停的姿轨协同控制方法 | |
CN106483466B (zh) | 一种卫星入轨阶段太阳电池阵输出电流的估算方法 | |
CN112572835B (zh) | 一种具有姿态切换的卫星在轨角动量管理及控制方法 | |
Zappulla et al. | Experiments on autonomous spacecraft rendezvous and docking using an adaptive artificial potential field approach | |
CN105912819A (zh) | 一种地月l1拉格朗日点转移轨道的快速设计方法 | |
US8262028B2 (en) | Spacecraft and method for operating the spacecraft | |
CN106595673A (zh) | 面对地球静止轨道目标操作的空间多机器人自主导航方法 | |
Huang et al. | Guidance, navigation, and control system design for tripropeller vertical-take-off-and-landing unmanned air vehicle | |
CN110334439A (zh) | 一种利用气动力辅助的超轨道速度飞行技术方案 | |
CN110647163B (zh) | 对geo空间目标持续可见光探测的绕飞轨道设计方法 | |
CN107918400B (zh) | 一种空天飞行器在轨操作相对位置姿态联合控制方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |