大气阻力引起卫星轨道衰减的数值模拟

摘要

本文以德国低轨道卫星CHAMP为例,联合考虑地球扁率和大气阻力摄动的影响,对相应摄动方程进行数值积分,计算轨道根数变化,并进而计算得到卫星空间位置,由此模拟考察大气阻力引起的轨道高度衰减.模拟中使用综合考虑了太阳辐射和磁暴等多种因素影响的最新国际大气标准JB2008模式来计算热层大气密度.选取CHAMP卫星轨道高度自然衰减(无点火提升卫星高度操作)的2005全年进行模拟;为了考察不同年份阻力系数的可能变化,对2002年1-3月处在较大高度的轨道也进行了模拟.考虑到CHAMP卫星的特殊几何构形及飞行高度的热层温度条件,取阻力系数大于2.8,并在一定范围内变化,以求得模拟与实际轨道衰减符合较好.结果表明,对于2005年,阻力系数为2.91时模拟得到的轨道高度的衰减与实际轨道衰减符合得最好,模拟与实际轨道半长轴全年的标准偏差为81 m;在卫星高度稍高的2002年,模拟的最佳阻力系数为3.0;模拟所得最佳阻力系数值比传统使用的值2.2大30％以上.由于在模拟中忽略了高阶保守力分量,所得近/远地点高度没有出现实际轨道所显示的周期性起伏.