天体力学：第六章 受摄运动方程.pptVIP

2、数值法 数值法是解微分方程的数值积分法。可以直接按照直角坐标形式的受摄运动方程(6.3)，进行数值积分，以已知的某一初始时刻天体的坐标和速度向量作为初值，逐步求得在任一时刻天体的坐标和速度向量。也可以按照拉格朗日(6.52)或牛顿受摄运动方程(6.96)进行数值积分，推求任一时刻卫星的轨道根数。但直接按(6.3)解算，计算公式简单，可以比较严密地顾及各种摄动力，因此解算精度较高。 数值法的缺点 不能给出问题的形式解，而且解算只能按选定的步长，逐步推算，运算时间较长，受累积误差影响较大。但是由于现在计算机技术的发展，不论在计算速度还是精度上都不是主要问题，因此数值法在实际工作中得到广泛应用。 常用方法： 龙格-库塔(Runge-Kutta)法 亚当姆斯-考威尔(Adams-Cowell)法 龙格－库塔法 亚当姆斯－考威尔法 第八节 常见的摄动力 在天体运动过程中，若以二体问题轨道为无摄运动，则会受到很多种摄动力。以人造地球卫星为例： 日月行星摄动：日月行星对卫星的引力产生的摄动 阻尼摄动：大气阻力引起的摄动 形状摄动：地球非球形对称引起的摄动 光压摄动：太阳光压产生的摄动 ……等等 1、其他天体的引力摄动 行星P在中心力作用下的运动方程： 行星Pi产生的摄动加速度（引潮加速度）： 多体引力摄动： 2、阻尼摄动 3、形状摄动 太阳光压引起的轨道摄动 太阳光压：卫星受太阳照射而产生的推力 摄动加速度： ：与卫星表面反射性能有关的系数 ：卫星在垂直于太阳方向上的截面积 ：与太阳相距一个地球平均半径处太阳光压的压强 ：太阳中心的地心向径 ：卫星的质量 在高出地面1000km以上时，太阳光压的影响大于大气阻力的影响，但在高度较低时，大气阻力的影响要大得多； 当太阳黑子最大时，在900km左右的高度两种摄动相等； 当太阳黑子最小时，在500km左右的高度两种摄动相等； 间接太阳光压引起的摄动一般小于直接太阳光压影响的15%； 王正涛.卫星跟踪卫星测量确定地球重力场的理论与方法[博士论文].武汉大学，2005 小结 摄动的概念 正则方程组、正则根数 拉格朗日行星运动方程 牛顿方程 泊松括号与拉格朗日括号 摄动导数 常见的摄动力 [Am, Bn] 不管Am是哪一个轨道根数，只有Bn＝p时，[Am,Bn]才可能不等于零 [Bm, Bn] 以上计算出了关于轨道根数（Ω，i，ω，p，e，τ）所有组合的拉格朗日括号。 第六节 牛顿方程 1、摄动导数 定理一： 摄动导数： 定理二： 摄动导数的作用 例一：平近点角M的摄动导数 例二：偏近点角E的摄动导数 例三：真近点角 f 的摄动导数 求下列函数的摄动导数 2、用摄动加速度分量表示的受摄方程 针对非保守摄动力的情况，牛顿提出了用摄动加速度分量表示天体受摄运动方程的方案，故称用摄动加速度分量表示的受摄运动方程为牛顿方程，牛顿方程实际上就是建立了轨道根数变化率同受摄加速度分量之间的关系，这一关系可利用摄动导数来推导。 牛顿方程 牛顿方程与拉格朗日行星运动方程 拉格朗日行星运动方程只是牛顿方程的一种特殊形式，即存在摄动函数的情况 第七节 受摄运动方程的解法 受摄运动方程的形式 直角坐标形式（6.3）式 轨道根数形式：拉格朗日运动方程（6.52）式 牛顿方程（6.96）式 求解方法： 分析法 数值法 1、分析法 分析法是利用摄动力与地球中心引力之比小于10-3这一特性，用级数展开和迭代的方法解拉格朗日或牛顿受摄运动方程。 解算时把轨道根数的摄动分为长期摄动、长周期摄动和短周期摄动。 地球引力场引起轨道摄动的特性 轨道摄动：长期摄动，周期性摄动 长期摄动：呈系统性变化，与时间成线性关系，由轨道根数a、e 和 i 引起 周期性摄动：短周期摄动（周期小于1天）和长周期摄动（周期为几十至几百天），其中前者是真近点角 f 或格林尼治恒星时SG的周期函数，后者是近地点角距ω或升交点赤经Ω的周期函数 长期摄动、长周期摄动和短周期摄动按数值的大小可分为一阶摄动、二阶摄动、三阶摄动等，以量级为地球引力场C2,0的项称为一阶，量级为C22,0的项称为二阶，等等。 分析法的解按精度分为一阶解、二阶解等，一阶解就是求出轨道根数的全部一阶摄动(包括一阶长期摄动、一阶长周期摄动和一阶短周期摄动)和二阶长期摄动。一阶解的精度达到10-6。二阶解要求出全部二阶摄动和三阶长期摄动，精度达到10-8。 分析法的缺点： 只能对各种摄动力单独求解，然后相加，不能顾及各种摄动力之间的相互影响，而且难以推导出完整的二阶