卫星编队轨道参数的确定方法.docxVIP

卫星编队轨道参数的确定方法

1卫星轨道参数取值的动力学分析方法

随着空间技术的发展，世界上越来越多的空间飞机正在发射。大卫星的研制成本高,研制周期长。与单颗大卫星相比,编队小卫星以其高效,灵活,低成本,高可靠性等诸多优势越来越受到广泛的关注,多颗小卫星协同作业,实现单颗大卫星的功能,应用前景十分广泛。美国国家航空航天局(NationalAeronauticsandSpaceAdministration,NASA)已成功发射OE-1和Landsat-7进行编队实验。

如何进行编队构型设计是当前研究的热点。目前有运动学方法,如文献提出了一种通过求解C-W方程的解析解和任务要求的初始条件求得卫星轨道参数的方法。C-W方法要求解二阶微分方程组的解析解,并且要根据要求确定初始6个坐标,再根据坐标确定卫星轨道参数。计算方法相对复杂。也有轨道动力学方法,如文献提出了利用相对轨道参数求解绕飞卫星轨道参数的方法,并且给出了在水平面进行构型设计的例子。

本文基于动力学方法,从轨道参数描述的相对运动方程出发,得出卫星的轨道参数取值。由于轨道参数是相对运动方程解析解的积分常数,在二体的条件下保持为常数［1］［3］,更加便于队形的设计。

2基于卫星轨道动力学理论的基础

2.1u3000坐标系

首先定义两个坐标系。

地球惯性坐标系。以地心为坐标原点O,地心和春分点连线为X轴,赤道面指向北极点为Z轴,由右手定则确定Y轴的坐标系为地球惯性坐标系。

卫星轨道坐标系。坐标系原点与参考卫星质心重合,Y轴与速度方向重合,X轴从地心指向参考卫星质心的方向,Z轴与X轴构成右手系。在描述编队相对运动时,是以参考星的轨道坐标系作为动力学方程的基准坐标系,定义相对运动状态［2］。

2.2卫星轨道参数的选择

通过六个经典的轨道参数来描述卫星在宇宙空间的位置,分别是:轨道半长轴a,偏心率e,轨道倾角i,近地点幅角ω,升交点赤经Ω,近地点时刻tp。

卫星的轨道参数分为三组,第一组轨道半长轴和偏心率决定轨道的大小和形状,第二组轨道的倾角,升交点赤经,近地点幅角决定了卫星轨道的位置,第三组参数近地点时刻决定了卫星在轨道上的瞬时位置［4］。

M为卫星的平近点角,有时也可以用参数M来代替近地点时刻描述卫星的瞬时位置。其中n表示卫星的平均角速度。

由上述参数构成的系统可以精确描述卫星的运动。

2.3卫星和绕飞卫星轨道倾角和偏心率分布

构成编队的卫星应该满足以下条件［3］:

1)轨道长半轴相同(即轨道周期相同),确保绕飞轨道闭合;

2)轨道偏心率是小量,即参考轨道是近圆轨道,参考卫星和绕飞卫星的轨道倾角和偏心率有微小差异,确保编队飞行的各个卫星适当分离;

3)平纬度幅角(λ=ω+M),或者纬度幅角(u=ω+f)基本相同,确保航天器编队的伴随运动。

用Δ表示两颗卫星参数的差,即Δσ=σh-σc,σ表示卫星的各个轨道参数,以下标c和h分别表示参考卫星和环绕卫星,于是,上述条件可以描述为:

2)ec,eh,Δe,Δi,ΔΩ,Δu,Δλ均为小量。

2.4卫星计划的基本模型

2.4.1绕坐标轴运动模型

假设地球惯性坐标系为坐标系E,参考卫星的轨道坐标系为坐标系C,跟踪卫星的坐标系为坐标系H,则根据坐标系的定义可以求得从坐标系E到坐标系C和坐标系H的转移矩阵。

Mj(θ)表示绕坐标轴j旋转角度θ的初等转移矩阵。在参考卫星的轨道坐标系中建立环绕卫星的运动方程,得到相对运动模型［5］:

其中Rc=ac(1-eccosEc),Rh=ah(1-ehcosEh),E表示卫星轨道的偏近点角。式(2)就是精确的卫星相对运动方程。

2.4.2卫星编制的描述

由数学知识可知:当Δθ为小量时,有cos(θ+Δθ)≈cosθ-sinθ·Δθ,sin(θ+Δθ)≈sinθ+cosθ·Δθ。因为ec,eh,Δe,Δi,ΔΩ,Δu,Δλ均为小量,将其代入相对运动方程并略去高阶小量,有:

为了直观地描述卫星编队的任务要求,可以将卫星编队描述为如下五个参数:1绕飞椭圆的半短轴:p=aeA,该参数决定了绕飞椭圆的大小。2垂直方向的振幅:s=aeB,该参数决定了绕飞轨迹在z方向运动的振幅。3初始相位差:α=θ-ωc-η,该参数和p、s一起决定了绕飞轨迹平面在惯性系下的空间指向。4绕飞椭圆的初始相位:ψ=Mc+θ,该参数决定了初始时刻环绕卫星在绕飞椭圆上的位置。5绕飞椭圆