浅谈小卫星虚拟卫星的发展.docxVIP

浅谈小卫星虚拟卫星的发展

1小卫星绕飞轨道动力学和控制

几个小卫星的飞行组成了一定的飞行轨迹，在分布方式上形成了“虚拟卫星”或网络卫星，这些卫星构成的板块称为网络板块。

从应用入手,研究卫星编队飞行完全是一个创新和开拓性任务,非常富有挑战性。文献研究了小卫星绕飞轨道动力学和控制,这是卫星编队飞行的一个特例,也是最简单的编队飞行,仅在同一轨道平面内编队,而且研究动力学模型只停留在线性化的Hill方程。国内外也有学者对卫星编队飞行进行过研究,并获得一些很好结果,但是大部份仅从简化模型或者仅从控制策略某一方面研究。本文从应用入手,根据实际应用需求,研究编队飞行三种动力学模型,同时进行定位卫星和区域导航卫星两个实例数学仿真,并且比较这些模型差别和适合的应用场合。

2大量燃料难以实现

根据需要,小卫星编队飞行在理论上可以设计任意队形,但在实际中由于需要消耗大量燃料所以很难实现。本文着重研究多颗卫星相对轨道动力学并且尽量利用它们的动力学特性,低燃耗或不消耗燃料(运行时间较短)来设计编队飞行的轨迹,从而满足工作应用需要。

2.1动态模型

根据上述要求,作者按不同应用要求研究了三种编队飞行动力学模型。

(1)构建w方程的一般方程

Hill方程的坐标系定义为X轴指向飞行方向,Y轴背向地心,Z轴垂直于轨道平面,构成右手坐标系。Hill方程表示如下

??x-2n˙y=Fx/m=Ax??y+2n˙x-3n2y=Fy/m=Ay??z+n2z=Fz/m=Az}(1)x???2ny˙=Fx/m=Axy??+2nx˙?3n2y=Fy/m=Ayz??+n2z=Fz/m=Az???????(1)

式中n为轨道角速度,对圆轨道来说n=常数。

(2)飞行周边卫星到地心距离

当对重力不进行一次近似,即不进行线性化,则方程式(1)将变为下列方程:

??x-2n˙y-(W2-n2)x=Fx/m??y+2n˙x-(W2-n2)(y+R0)=Fy/m??z+W2z=Fz/m}(2)x???2ny˙?(W2?n2)x=Fx/my??+2nx˙?(W2?n2)(y+R0)=Fy/mz??+W2z=Fz/m???????(2)

式中W2=μ/R3,R为编队飞行周边卫星到地心的距离;n2=μ/R3030,R0为编队飞行中心卫星到地心的距离;R=√x2+(y+R0)2+z2R=x2+(y+R0)2+z2????????????????√。

(3)飞行动力学方程

描述编队飞行以两个卫星相对运动为基础,即一个卫星为中心,另一个卫星相对它运动。精确模型是以地心坐标系来描述相对运动,这是以两个单独存在的卫星建立相对运动方程(即开普勒方程),然后把相应变量相减,得到周边卫星相对中心卫星的运动方程。这就是开普勒相对运动,它是两个椭圆运动的合成。

中心卫星以地心赤道坐标系为基础的运动方程:

??x1=-μx1r31+A1x??y1=-μy1r31+A1y??z1=-μz1r31+A1z}(3)x??1=?μx1r31+A1xy??1=?μy1r31+A1yz??1=?μz1r31+A1z?????????????(3)

周边卫星以地心赤道坐标系为基础的运动方程:

??x2=-μx2r32+A2x??y2=-μy2r32+A2y??z2=-μz2r32+A2z}(4)x??2=?μx2r32+A2xy??2=?μy2r32+A2yz??2=?μz2r32+A2z?????????????(4)

式中r1=√x21+y21+z21r1=x21+y21+z21??????????√;r2=√x22+y22+z22r2=x22+y22+z22??????????√;Ax,Ay,Az为三轴方向的轨道摄动加速度。利用式(3)、(4)得到两个卫星地心坐标系的绝对位置,然后相减,得到相对位置Δx,Δy,Δz;再根据编队飞行轨道坐标系和地心坐标系之间坐标转换,最后得到所需编队飞行动力学方程:

[xyz]=B[ΔxΔyΔz](5)

式中B为坐标转换矩阵。

当中心卫星初始条件(x10,y10,z10和˙x10?˙y10?˙z10)与编队卫星初始条件(x20,y20,z20和˙x20?˙y20?˙z20)和摄动加速度A都已知,则方程式(3)、(4)可以解算。

根据地心坐标系参数x,y,z和˙x,˙y,˙z与轨道六要素(a,e,i,Ω,ω,f)两者可以互换,则方程式(3)、(4)当A=0时,其解可以写成下列形式:

[x1y1z1]=a1(1-e21)1+e1cosf1[cosΩ1cos(ω1+f1)-sinΩ1sin(ω1+f1)cosi1sinΩ1cos(ω1+f1)+cosΩ1sin(ω1+f1)cosi1sin(ω1+f1)sini1](6)和[x2y2z2]=a2(1-e22)1+e2cosf2[