第二章衛星運動理論.docVIP

第三章　衛星運動理論 　　人造衛星運動軌道理論是描述人造衛星太空中運動的一門學科，是衛星大地測量的重要理論基礎。人造衛星進入軌道自動飛行階段後，也和自然天體一樣在萬有引力及其它力的作用下，遵循牛頓運動定律，因而可以利用天體力學研究其運動規律。但和自然天體相比，人造衛星的運動又有著特殊性，主要表現為（李征航等，1998）： 人造衛星離地球較近，因不能將人造衛星如同行星視為質點，必須考慮地球引力對衛星運動的影響。 人造衛星受到作用力較自然星體複雜，除地球引力場的作用外，還受到其它力的影響，如海潮、固體潮、大氣阻力、太陽幅射壓…等(Tapley,1989) 人造衛星的運動周期較一般自然天體的運動周期短。 自然天體只存在定軌問題，而人造衛星還存在軌道設計、衛星回收等問題（Montenbruck and Gill，2001）確算出。 根據牛頓萬有引力定律，可知宇宙中任意兩質點皆相互吸引；而從牛頓第二運動定律得知物體受外力作用時所產生加速度的大小與外力大小成正比，與該物體的質量成反北，加速度的方向與外力方向相同牛運動定律只成立於慣性坐標系統，故假設地球之質量為M，衛星質量為m，衛星受地球中心引力之影響而產生之加速度可表示為（Montenbruck and Gill，2001）： (3-1) 其中頓萬有引力常數為衛星在地心慣性座標系之位置向量由於衛星質量遠小於地球質量，可忽略不計，則(3-1)式可寫成： (3-2) 　　由式(3-2)可知二體運動方程為二階常微分方程，對3-2)式作一次積分，確定三個積分常數便得，再積分一次，確定另三個積分常數則，如下所示：　 (3-3)式中，C1~C6為六個積分常數，ｔ為時間函數。因此，只要能決定出六個積分常數，此微分方程的解與便可確定，依據能量、動量、角動量守恆積分等原則，可找出幾何或物理意義彼此獨立的六個軌道元素，如刻卜勒（Keplerian）元素(Escobal,1975)、Hill元素(Cui, 1992)、Equinoctial元素(Velez,1972)等。在此僅介紹刻卜勒軌道元素（曾清涼及儲慶美，1999），如圖3-1（a）所示： 決定刻卜勒橢圓的形狀和大小 a：軌道長半徑(semi-major axis) e：軌道離心率(eccentricity) 決定衛星軌道平面與地球體之間的相對位置 Ω：昇交點赤經(right ascension of the ascending node)，為地球赤道平面上，昇交點與春分點之間的地心夾角。 i：軌道面傾角 (orbital inclination)，為地球赤道面與衛星軌道面夾角。 決定刻卜勒橢圓在軌道平面上的方向 ω：近地點變角 (argument of perigee)，為軌道平面上昇交點與近地點之地心夾角。 決定衛星在軌道上的瞬時位置 f：真異常角 (true anomaly)，軌道平面上衛星與近地點之夾角。在幾何圖中，無法標示平異mean anomaly)，可由真異常角與平異常角之間的關係，透過離心異常角Ｅ(eccentric anomaly)求得，真異常角與離心異常角的關係可圖3-1：(a)為刻卜衛星軌道元素　(b)為衛星軌道真異常角與離心異常角關係 據幾何學中定義出離心異常角與真異常角之關係為： (3-4) 又由卜勒方程式可得E與M之關係式： (3-5) 在二體問題中，由於只考慮到地球的引力，此力為保守力，故能達到能量守恆，而對於一個能量守恆之衛星運動系統，a,e,i,,皆為時間常數，只有異常角會隨時間變化(蕭泰中，1999)。另外，在二體問題中，軌道平面是固定的，不會擺動，地心則是橢道橢圓的一個焦點。 (3-6) 其中，：地球引力位引起之擾動加速度 　　　：多體引起之擾動加速度 　　　：地球固體潮引起之擾動加速度 　　　：海潮引起之擾動加速度 　　　：空氣阻力引起之擾動加速度 　　　：太陽幅射壓引起之擾動加速度 　　　：地球幅射壓引起之擾動加速度 　　　：相對論效應引起之擾動加速度 表3-1提供了各種擾動量之間數量級的比較，以下各節則將深入介紹上述8種擾動力模式，而擾動力對於軌道的影響，將詳述於第五章 10-3 6*10-3 2*10-3 5*10-5 非球體引力場其餘項總和 1*10-4 9*10-5 5*10-6 3*10-7 多體擾動 1*10-6 1-2*10-6 2*10-6 5*10-6 固體潮擾動 2*10-7 1-2*10-7 3*10-8 1*10-9 海潮擾動 3*10-8 2*10-8 2*10-9 1*10-10 大氣阻力擾動 1-2*10-10 1-2*10-10 3*10-12 0 太陽幅射壓擾動 5*10-9 5*10-8 4*10-9 1*10-7 地