第六章航天器轨道摄动

第六章航天器轨道摄动主讲教师：杏建军2020年5月20日2授课内容1.近地空间轨道摄动力2.轨道摄动理论3.具体应用31.1什么是轨道摄动卫星在空间绕地球运行，除了受地球中心引力作用外，还将受到地球非球形引力、太阳、月亮和其它天体引力的影响，以及太阳光压、大气阻力和地球潮汐力等因素的影响，使航天器的运动偏离二体问题的轨道。41.1什么是轨道摄动轨道摄动是相对于二体问题来说的，二体问题的轨道称为无摄轨道或正常轨道，考虑其它摄动力影响的轨道为摄动轨道。其他力相对地球中心引力而言都是小量，因此统称为摄动力，这些力作用下的航天器运动称为摄动运动。51.2近地轨道空间摄动力1、地球引力势函数M地球质量R地球参考半径coscoscossinsinECFECFxRRyRzRnmP勒让德函数61.2近地轨道空间摄动力勒让德函数勒让德多项式71.2近地轨道空间摄动力nmCnmS地球引力系数01000mmm归一化的系数81.2近地轨道空间摄动力带谐项，与经度无关EGM2008JGM3WGS84常用的地球引力模型3220(10)JCO最大的一项EGM2008J2=0.0010826261738520000101110111000CSCCSS91.2近地轨道空间摄动力地球引力（地心固定坐标系ECF）ECIECIECFECFRrr101.2近地轨道空间摄动力地球J2项摄动力220202(sin)GMGMRUPCrrr22231(3sin1)2GMGMRUJrr2222222232311(31)(3sinsin1)22JRRZUJJiurrr111.2近地轨道空间摄动力地球J2项摄动力2222222232311(31)(3sinsin1)22JRRZUJJiurrr2222TJJJECFJUUUFUxyz22222522235551132TECIRZZZFJxyzrrrr121.2近地轨道空间摄动力地球J2项摄动力2222222232311(31)(3sinsin1)22JRRZUJJiurrr222211sinTJJJRICJUUUFUrrurui222224313sinsin2sinsincos2sincoscos2TRICRFJiuiuuiiur131.2近地轨道空间摄动力讨论题：仅考虑地球中心引力和J2项摄动力时，卫星的动力学方程是什么？作业：以x0=[-5292392.072;-4862.201380;3111662.355;-4136.781314;3101114.660;-4147.028008]为人造地球轨道卫星的初始条件，仅考虑地球中心引力和J2项摄动力，计算卫星轨道根数，并于二体结果进行比较。141.2近地轨道空间摄动力半长轴a在地球J2项摄动力作用下的运动151.2近地轨道空间摄动力偏心率e在地球J2项摄动力作用下的运动161.2近地轨道空间摄动力轨道倾角i在地球J2项摄动力作用下的运动171.2近地轨道空间摄动力升交点赤经在地球J2项摄动力作用下的运动181.2近地轨道空间摄动力近地点角距在地球J2项摄动力作用下的运动191.2近地轨道空间摄动力平近点角在地球J2项摄动力作用下的运动201.2近地轨道空间摄动力2、其他天体引力3GMsrrsrmMrs质点M对m的引力加速度讨论：这个引力加速度是其他天体对卫星最终引力摄动吗？211.2近地轨道空间摄动力2、其他天体引力3GMssr质点M对地球的引力加速度综合这两项，得到其他天体对卫星的引力摄动33GMssrsrsr331niiiiiGMssrsrsr多个天体221.2近地轨道空间摄动力2、其他天体引力33GMssrsrsr33222232321212rrsrsrssssrsreeeesrsr假设的合理性：地月距离约38万公里，地日距离约1.5亿公里33113rsssreesr33GMrssrsrreeee231.2近地轨道空间摄动力2、其他天体引力33GMrssrsrreeeesree当32GMrsrresree当3GMrsrre讨论：估计地面物体受到的月球引力摄动（地球质量约是月球的81倍）241.2近地轨道空间摄动力2、其他天体引力33GMssrsrsr计算其他天体引力摄动，需要知道其他天体相对于地球的位置矢量s。s矢量的计算有两种方法：简易方法和精确方法（JPL星历DE450，包含了1600年到2170年的星历数据）251.2近地轨道空间摄动力2、其他天体引力JPL星历261.2近地轨道空间摄动力3、大气阻力大气阻力是影响近地轨道卫星的最大非引力摄动力，是一种耗散力，对卫星轨道的寿命有很大影响，也是天然的太空清洁师。大气阻力的精确建模非常困难：高层大气密度变化复杂，很难精确测量；与卫星的姿态密切相关，对于外形复杂的卫星迎风面积较难预测；与卫星的外表面特性密切相关。271.2近地轨道空间摄动力3、大气阻力212DrrACvmreDC阻力系数，无单位量，描述了大气与卫星表面材料的相互作用情况，一般取2.0~2.3。rvvωr卫星相对于大气的速度矢量高层大气密度281.2近地轨道空间摄动力3、大气阻力高层大气模型Harris-priester模型Jacchia1971模型Jacchia1977模型MSIS86模型291.2近地轨道空间摄动力4、其他摄动力太阳光压摄动地球辐射压力摄动地球潮汐摄动相对论效应30授课内容1.近地空间轨道摄动力2.轨道摄动理论3.具体应用312.1摄动方程322.1摄动方程332.1摄动方程342.2一般摄动法常数变易法352.2一般摄动法常数变易法二体运动，α为常数(,)(,)ttrrαvvα受摄运动，α为时间变量增加一个密切轨道条件362.2一般摄动法pddtvαaα完整的摄动方程61610iiiipiiddtddtrva拉格朗日和高斯分别根据这组方程得到一组计算行星摄动的方程。372.3拉格朗日行星摄动方程022220222222222022111cos1sin1cos1sin11sin12daRdtnaMdeeReRdtnaeMnaediRRidtnaeideRiRdtnaeeinaeidRdtinaeidMeRRdtnaeenaa0MMnt382.4高斯摄动方程222222222[sin(1cos)]11[sin(coscos)]cos1sin1sin1[cos(1)sin]cos1[(cos2)(1)sin]RIRICCRIRIaeffeffneeefffEfnaruifnaerufnaeierffffinaeperrMnfefffnaepp392.5轨道摄动相关的轨道根数402222222232311(31)(3sinsin1)22JRRZUJJiurrr2.6平均轨道根数地球J2项摄动的势函数地球J2项摄动的平均势函数22012avgJRUdM222222223/20131sin122(1)23avgJnriRUdfnRJhedMndt2hdfdtr412222222220003cos2352sin22331sin122EEEaeiJRnipJRnipJRMnniepJ2项摄动对轨道根数的长期影响2.6平均轨道根数422.6平均轨道根数太阳和月球的引力引起所有轨道根数周期性变化，但只有升交点赤经、近地点角距和平近点角有长期变化0.00338cos/Moonin0.00154cos/Sunin20.00169(45sin)/Moonin20.00077(45sin)/Moonin3na432.6平均轨道根数大气阻力：主要为轨道半长轴与偏心率的长期衰减，即轨道高度不断降低，偏心率也减小，因而轨道向圆轨道演化。442.6平均轨道根数212DDCAaVm大气阻力220DrevCAaaem每圈变化量卫星寿命revHLa45授课内容1.近地空间轨道摄动力2.轨道摄动理论3.具体应用463.1太阳同步轨道在地球J2项摄动影响下，卫星轨道面的升交点赤经的长期变化率与平太阳赤经变化率相同。223cos2EJRnip360==0.9856deg/365.25day太阳同步轨道与轨道高度、偏心率和轨道倾角相关；太阳同步轨道的轨道倾角要大于90°。473.1太阳同步轨道200400600800100012001400160096979899100101102103轨道倾角°轨道高度km483.1太阳同步轨道493.2冻结轨道（CriticalInclination）在地球J2项摄动影响下，卫星轨道的拱线不发生长期移动的轨道。冻结轨道只与轨道倾角相关；冻结轨道卫星与地面的高度始终不变。222352sin022EJRnip0063.43116.57ior503.2冻结轨道（CriticalInclination）主要用途：大椭圆轨道卫星Molniya轨道e=0.73，轨道周期12小时51本章小结近地空间轨道摄动力：地球非球形摄动三体引力摄动大气阻力摄动太阳光压摄动………….331niiiiiGMssrsrsr212DrrACvmre52本章小结轨道方程的解法：特殊摄动法：数值积分一般摄动法：求解析解61610iiiipiiddtddtrva53本章小结022220222222222022111cos1sin1cos1sin11sin12daRdtnaMdeeReRdtnaeMnaediRRidtnaeideRiRdtnaeeinaeidRdtinaeidMeRRdtnaeenaa54本章小结222222222[sin(1cos)]11[sin(coscos)]cos1sin1sin1[cos(1)sin]cos1[(cos2)(1)sin]RIRICCRIRIaeffeffneeefffEfnaruifnaerufnaeierffffinaeperrMnfefffnaepp55J2项摄动分析J2项摄动的平均势函数222223/231sin12(1)23avgiRnRJe本章小结2222222220003cos2352sin22331sin122EEEaeiJRnipJRnipJRMnniep