2008航天器动力学10-摄动理论答案.ppt

上述方程组原则上可以求解。但由于严重非线性，只可能进行数值计算，没有解析解。 另一方面，由于摄动的影响与万有引力相比及其微小，因此轨道要素的变化是缓慢的，我们关心的是轨道的长期摄动。 的变化缓慢。 的变化快。 当方程中的变量变化快慢明显不同时，理论分析中可考虑采用平均法对方程进行简化。 4、长期摄动影响 在数值计算中，如果变量变化缓慢，计算的步长可以大些；如果变量变化剧烈，计算的步长就要小些。如果两种情况同时存在（特征值相差大） ，方程称为“刚性方程” ，也称为“病态方程” ，需要特殊的算法。 平均法的思路是： 用 统一表示各轨道要素 在卫星运行的一个周期内，近似认为轨道要素没有变化，将方程在一个周期内积分，得到每运行一周轨道要素的微小增量 由于?i不是时间的显函数，而是真点角的显函数，因此… 虽然轨道根数的变化很复杂，但是其对 f 的积分有可能有解析解。经过积分，得 clear all syms mu e p theta f omega J2 Re inc disp [－－－－－符号运算－－积分－－－－－－ ] ; theta f+omega; r p/ 1+e*cos f ; dt_df r^2/sqrt mu*p ; %转换关系 fr -3*mu*J2*Re^2/2/r^4* 1-3* sin inc ^2* sin theta ^2 ; ft -3*mu*J2*Re^2/2/r^4* sin inc ^2*sin 2*theta ; fn -3*mu*J2*Re^2/2/r^4*sin 2*inc *sin theta ; dp_dt dt_df*2*r*sqrt p/mu *ft; dOMEGA_dt dt_df*r*sin theta *fn/sqrt mu*p /sin inc ; di_dt dt_df*r*cos theta *fn/sqrt mu*p ; p simple int dp_dt,f,0,2*pi inc simple int di_dt,f,0,2*pi OMEGA simple int dOMEGA_dt,f,0,2*pi 符号运算的程序 利用Matlab的积分功能 很容易得到结果！ 上述平均化过程消除了轨道要素在每一个运行周期内的波动，但能准确反映轨道要素的长期变化规律，即轨道的长期摄动。轨道要素在每个运行周期内的波动称为周期摄动。 地球扁平率不影响轨道的几何形状，只引起轨道平面的转动。 机械能守恒，动量矩不守恒 为什么？ 前面，在太阳同步轨道中，引用过公式 问题 但现在得到的公式是 利用 两者有什么关系？ 物理含义？ §3.4 大气阻力引起的摄动 地球周围的大气层厚度为1000－1500公里。 由于稀薄气体的持续作用，航天器的能量会不断耗散。 即航天器的速度与高度不断下降。 1、大气阻力对轨道高度和速度的影响 ? × 注意 由于阻力的影响，航天器的高度不断下降，但速度却是在增加的！ 与一些人的直观结论不同！ 但是，我们根据右图，又可以直观地证明。 因此，你能否根据有关的公式，推出这一结论？ V小 V大 2、大气阻力引起的摄动加速度 大气密度随高度变化，近似表示为： 为大气密度， r0为轨道上的某已知点，一般选近地点。 H为常数，与r0 有关。 为r0处的大气密度， 由于 r 随真近点角 f（半长轴与矢径的夹角） 变化，因此有 * Page * 第三章 航天器轨道运动的摄动 §3.1 轨道摄动的原因及处理思路 §3.2 轨道根数摄动方程 §3.3 非中心引力场引起的摄动 §3.4 大气阻力引起的摄动 §3.1 轨道摄动的原因及处理思路 前面介绍了航天器轨道运动的理想模型－－开普勒轨道。由于各种干扰因素的存在，航天器的轨道不再遵循开普勒运动规律。 这些干扰因素包括 ? 非均匀的地球引力场 ? 大气阻力 ? 太阳光压 ? 太阳、月球的引力 这些干扰因素将使航天器偏离其预期的轨道。但由于这些干扰因素与地球引力相比是微量，这种微小扰动通常称为摄动。 轨道摄动的处理方法 1、考威尔方法（Cowell,1910年） 基本思想：直接考虑干扰因素和中心万有引力。 二体系统的动力学方程 附加摄动项 在惯性坐标系中，可以写成 优点：该方法简单直观，适用范围广，对摄动项没有什么限制，可同时处理不同的扰动。 缺点：由于摄动项远小于地球引力项，为了在计算中反映出摄动的影响，要求计算精度高，积分步长小，运算量大，误差积累可能严重。 目前由于计算水平提高，该方法重新得到重视。 2、恩克方法（Encke,1857年） 基本思想：以开普勒轨道为基准轨道， 列写航天器偏离此轨道的动力学方程。 建立以? r为变量的动力学方程： 优点：由于轨道偏差变化缓慢，因此方程积分步