编程模拟飞船加速变轨过程-物理基础篇（3）Kepler轨道及其描述（上）

编程模拟飞船加速变轨过程物理基础篇（3）

Kepler轨道及其描述（上）

开普勒轨道及其描述是用若干要素表示出飞船在地球中心引力场中绕地球飞行具体形状的一种方法，我们希望在已知这些要素的前提下精确的描述运动轨迹，进而方便我们下一步在知道某时刻的速度与位矢时，求出这些轨道要素。

开普勒轨道·正文开始

对于地球与飞船组成的系统，由二体问题的知识可以推导出：

外力f可能包括稀薄大气阻力，光压力等。此式为非线性方程，求解困难，因此先考虑f=0的情况（实际上，要模拟的过程也是在忽略这些的前提下进行的）即：

因为地球质量远大于卫星，系统质心与地心可视为重合，就可以看作是中心引力运动。为了方便，下面都以中心引力运动来讲解开普勒轨道。

速度分布

结合第一篇，速度可表示为横向、径向两个分量：

横向速度，径向速度

代入 以及，将径向速度中的替掉，有：

所以速度大小

也即

注：a指半长轴。这也说明了离引力中心越远，速度越小

几个具有特殊意义的速度如下：

1、第一宇宙速度，在圆形轨道上运行的速度，轨道半径=天体半径，

2、第二宇宙速度，对于抛物线，a=∞，此时，当飞船位于近地点时速度最大。若近地点刚好在天体表面上，则称该点速度为第二宇宙速度，其意义是刚好不再做椭圆运动的速度临界点，所以也叫逃逸速度

3、若在近地点处速度比逃逸速度大，那么飞船运动轨迹为双曲线，此时

可知它即使是在无穷远处也有一个速度，成为双曲线的过速度。

 

通常所讲的开普勒轨道指椭圆轨道，后面提到Kepler轨道时也特指椭圆轨道

 

接下来介绍一下椭圆的偏心率e，半通径p，长半轴a，短半轴b之间的关系，以便后面轨道能量、轨道周期、轨道要素关系的表达。

，

因椭圆上r>0，且r有限，所以

上式表示角速度方向不变，即在椭圆上沿一个方向运动。称此椭圆为轨道，运动为椭圆轨道运动。

 

轨道能量

因中心引力为有势力，则单位质量的质点在轨道上能量可表示为

因能量守恒，代入e,a,p,b的关系消去v和r：

可见能量大小取决于半长轴

轨道周期

需先引入“面积公式”概念

取极限，因为θ趋近于0时，sinθ~θ，再次代入就会有：

说明：单位时间扫过面积（面积速度）是常数，且为动量矩的一半

若周期用P表示，那么运动一周时间内，径矢扫过面积为πab，算定积分可得：

可以将H替掉：（这就是开普勒第三定理）

轨道要素

Kepler轨道是飞船相对地球的运动轨道，该轨道位于一个过地球球心、方向不变的动平面内，该平面称为轨道面

Kepler轨道是周期性椭圆运动，地球固定在椭圆的一个焦点上。称椭圆的近-远地点连线为拱线，拱线在轨道面上的方向不变。另外，引入升降交点概念：由南半球穿过赤道，轨道与赤道平面的交点叫升交点；由北半球穿过赤道与赤道的交点叫降交点

因此，想确定飞船的运动轨迹，需要知道：轨道面的方向，椭圆轨道的大小形状，拱线所在轨道平面内的方向，卫星在轨道上的位置。

 

用于确定飞船轨道的要素有6个，换句话说，已知这6个要素，便可确定一个轨道。这些要素分别是：

升交点赤经Ω，轨道倾角i，近地点幅角ω，半长轴a，偏心率e，卫星过近地点的时刻t_p  实际上，本程序开发时无需考虑Ω和i，然后由于篇幅有限，相关定义还请自行百度。

 

轨道的时间历程

由可看出，只要知道了真近点角f与时间的关系，就能知道r随时间的变化规律。向中代入和以引入t，消去H，就可推出，进而消去r得

该方程没有解析解，因此f和t的关系不明显。其实，若改用另一个几何参数，偏近点角E来描述，则更容易给出位置和时间的关系。

 

图中几何关系：acosE-ae=rcosf

由轨道方程可导出，带入到几何关系中，有：

。由此得到轨道方程的另一形式：

r=a(1-ecosE)

将该式与比较，得到

以及

*注：该微分方程很难从上面直接推导出来，现仅结合结论和几何关系证明该式的正确性。

，因此只需证明：

由上图几何关系可看出

将该关系和r--cosf，r--cosE关系式带入，整理得到：

该式显然成立，证明完毕

 

将f--E的两关系式带到f与t的微分关系式中，得到

积分，左边从t_p到t，右边从0到E，结果为：

定义平近角点M，，显然，M和t有线性关系。上面三式结合，可推出：

M=E-esinE

此式成为Kepler方程。

顺便讲一下迭代关系。

由以上的分析可知：

1.若飞船位置（用E表示）已知，过近地点的时刻t_p已知，那么可以求得对应该位置的时刻t

2.同理已知t和t_p，也可通过M算出E

迭代关系可取：

，或

可以取E(0)=M

另外，f和E还有一个更通用的公式：

相关坐标系的定义

原教材中还有对坐标系的介绍（第一轨道坐标系，第二轨道坐标系）和其间的变换矩阵。但由于本程序模拟的内容将原本三维的运动简化为二维，因此无需了解太多，想了解的读者可自行百度。