姜 宇，李恒年 ，寶音賀西

(1.宇航動力學國家重點實驗室,西安710043; 2.西安衛(wèi)星測控中心,西安710043;3.清華大學航天航空學院,北京100084)

全球星座可用于導航、通信、環(huán)境監(jiān)測、數(shù)據(jù)中繼、目標識別等諸多領(lǐng)域.全球星座一般采用Walker星座或包含Walker星座的復合星座[1].Walker星座[2]由具有相同軌道半長軸和軌道傾角的多個衛(wèi)星組成，分為若干個軌道面，每個軌道面內(nèi)的衛(wèi)星數(shù)目相同，不同軌道面內(nèi)的衛(wèi)星之間的相對相位保持特定的關(guān)系.Walker星座一般以N/P/F表示，N代表星座中的衛(wèi)星總數(shù)，P為星座軌道面數(shù)，F(xiàn)為相位因子.典型的全球星座有美國的GPS、俄羅斯的Glonass、歐盟的Galileo等.Glonass采用24/3/1的Walker星座，Galileo采用27/3/1的Walker星座，早期GPS采用24/3/2的Walker星座.星座衛(wèi)星之間由于入軌偏差及所受攝動力的不同，而發(fā)生星座衛(wèi)星軌道逐漸漂移，進而引起星座整體構(gòu)形的漂移，從而影響星座的導航、覆蓋等性能.為了保持星座的整體構(gòu)形，需要對星座衛(wèi)星的軌道進行控制.

星座軌道控制主要包括星座衛(wèi)星入軌后轉(zhuǎn)移軌道段的軌道控制和長期保持控制[1,3].星座軌道控制的目的是一方面保持星座中每個衛(wèi)星的軌道偏差在允許的偏差范圍內(nèi)，另一方面保持星座的整體空間構(gòu)形漂移較小，從而滿足應(yīng)用對星座構(gòu)形的需求.

本文主要針對運行在中高軌道上的全球Walker星座的構(gòu)形長期保持控制問題，通過研究星座衛(wèi)星的軌道攝動及構(gòu)形相對漂移，研究了星座構(gòu)形長期保持的控制方法及控制策略.對于星座構(gòu)形的保持，研究了星座構(gòu)形攝動補償法和星座構(gòu)形數(shù)值微分修正法.最后給出了全球星座構(gòu)形長期保持控制策略并進行了仿真驗證.

1 星座衛(wèi)星軌道攝動運動分析

對于中高軌道的星座衛(wèi)星來說，其主要的攝動力為地球扁率J2項、日月引力以及太陽光壓.地球扁率J2項對衛(wèi)星軌道半長軸、偏心率、傾角的長期影響為零，對衛(wèi)星的升交點赤經(jīng)、近地點幅角和平近點角均有長期影響.日月引力對衛(wèi)星軌道傾角和升交點赤經(jīng)都有長期影響.衛(wèi)星軌道升交點赤經(jīng)、緯度幅角的長期攝動主要由J2項引起，軌道傾角的長期攝動主要由日月引力引起，偏心率矢量的長期攝動主要由太陽光壓引起.

1.1 地球扁率攝動分析

考慮地球扁率J2項攝動，則星座衛(wèi)星平均軌道根數(shù)的長期變化率[4]為

(1)

可見地球扁率J2項攝動對衛(wèi)星軌道半長軸、偏心率和傾角沒有長期影響，對升交點赤經(jīng)、近地點幅角和平近點角均有長期影響.在J2項攝動的作用下，升交點赤經(jīng)的變化率、近地點幅角的變化率以及平近點角的變化率均與衛(wèi)星的軌道半長軸、偏心率和傾角有關(guān).

1.2 日月引力攝動分析

太陽引力攝動引起的衛(wèi)星軌道傾角長期變化率為

(cosβssinisinΩ-sinβscosissinicosΩ+sinβssiniscosi)

(2)

式中，t為時間，βs為太陽視運動的黃經(jīng)，is為黃道傾角，ns為地球繞太陽公轉(zhuǎn)的角速率.

如果不計βs的周期項，則化為

(3)

對太陽同步軌道，則式(2)化簡為熟知的結(jié)論作為其特例

(4)

月球引力攝動引起的衛(wèi)星軌道傾角長期變化率為

sin2imcosi)-2cos2imsinicos(Ω+Ωm)]

(5)

可見日月引力攝動對衛(wèi)星軌道傾角長期變化率的影響與衛(wèi)星的軌道半長軸、傾角和升交點赤經(jīng)有關(guān)，且軌道半長軸越大，日月引力攝動引起的衛(wèi)星軌道傾角長期變化率的大小越大.

1.3 太陽光壓攝動分析

(6)

其中

(7)

(8)

A為航天器垂直于光線方向的截面積，m為衛(wèi)星質(zhì)量，K為光壓系數(shù)，全吸收時K=1，ps為光壓強度.

2 星座構(gòu)形漂移運動分析

(9)

可見，Walker星座的衛(wèi)星間的平根數(shù)半長軸差、傾角差、偏心率差均可引起星座衛(wèi)星軌道面法向、軌道沿跡方向的相對演化.即使將Walker星座的平均軌道根數(shù)的半長軸、偏心率和傾角均完全按照設(shè)計結(jié)果控制至相等，此時地球扁率J2項攝動不會引起星座構(gòu)形的演化；但日月引力引起Walker星座不同軌道面上的衛(wèi)星傾角攝動不同，太陽光壓引起Walker星座不同軌道面上的衛(wèi)星偏心率矢量攝動不同，進而引起星座衛(wèi)星之間偏心率差δe、傾角差δi的絕對值逐漸增大，由式(9)知存在半長軸差、傾角差或偏心率差的Walker星座衛(wèi)星間的J2項攝動不同，故而引起星座構(gòu)形的漂移.

3 控制策略計算方法

Walker星座構(gòu)形的控制需采用相對控制的方法，由于在入軌偏差、控制執(zhí)行偏差、J2項攝動力、日月引力等的作用下，星座構(gòu)形整體漂移，控制方法只需保證星座整體幾何構(gòu)形的穩(wěn)定和單星軌道滿足特定要求即可.Walker星座整體構(gòu)形的控制可采取攝動補償法和數(shù)值微分修正法兩種方法.

3.1 星座構(gòu)形攝動補償法

星座構(gòu)形攝動補償法[5-8]通過控制星座衛(wèi)星間的半長軸差和傾角差來達到間接控制升交點赤經(jīng)差和相位差的目的，記

(10)

令

λ1=ω1+M1

(11)

(12)

Φ為J2項作用下，星座衛(wèi)星軌道半長軸偏差和傾角偏差對升交點赤經(jīng)以及相位角的影響矩陣.Δt時間內(nèi)在J2項攝動力的長期作用下，星座衛(wèi)星間升交點赤經(jīng)差和相位差的漂移分別為ΔΩ和Δλ，則要消除J2項攝動力對星座衛(wèi)星間升交點赤經(jīng)差和相位差的漂移影響，所需的星座衛(wèi)星間半長軸和傾角的初始調(diào)整量為Δa和Δi，則

(13)

其中

(14)

如此，通過控制Walker星座衛(wèi)星之間的半長軸差和傾角差，達到了對升交點赤經(jīng)差和相位差聯(lián)合控制的目的.

Walker星座軌道的初始設(shè)計是平均軌道根數(shù)，星座各衛(wèi)星的半長軸、傾角均相等.設(shè)預定的星座控制周期為Δt，如果星座設(shè)計壽命為10年，且10年內(nèi)不做控制，則Δt=86400×365×10s.如果星座設(shè)計壽命為10年，但控制周期為3年，則Δt=86400×365×3s.

星座構(gòu)形攝動補償法的計算過程為：考慮地球形狀攝動力、日月引力、太陽光壓，首先根據(jù)初始設(shè)計的軌道平根數(shù)做平根數(shù)軌道預報，得出Δt時間后的ΔΩ和Δλ，由公式(13)計算所需的星座衛(wèi)星間半長軸和傾角的初始調(diào)整量Δa和Δi，更新星座衛(wèi)星初始軌道根數(shù)，再做平根數(shù)軌道預報，看Δt時間后的ΔΩ和Δλ是否滿足要求，如此迭代循環(huán).該過程一般執(zhí)行1至3次就可獲得滿意的結(jié)果.

3.2 星座構(gòu)形數(shù)值微分修正法

同星座構(gòu)形攝動補償法不同，星座構(gòu)形數(shù)值微分修正方法是基于瞬根數(shù)的.根據(jù)是否做傾角修正，構(gòu)形數(shù)值微分修正方法又有兩種情形，其數(shù)值計算步驟為：

方法1.僅做半長軸修正情形

如果衛(wèi)星攜帶推進劑剩余不多，則希望在盡量節(jié)約推進劑的情況下進行星座構(gòu)形的保持，此時，不進行軌道面外的控制，僅采用半長軸修正來最大限度地保證星座構(gòu)形的穩(wěn)定.根據(jù)入軌后的軌道根數(shù)，進行高精度軌道預報，根據(jù)Δt時間后的廣義相位差增量Δφ來計算半長軸的修正量，廣義相位差增量等于相位差的增量與升交點赤經(jīng)差的增量在衛(wèi)星軌道面內(nèi)的投影的和

Δφ=ΔΩcosi+Δλ

(15)

微分修正法的半長軸修正步長為

(16)

方法2.半長軸和傾角都做修正情形

根據(jù)入軌后的軌道根數(shù)，進行高精度軌道預報，根據(jù)Δt時間后的相位差增量Δλ、升交點赤經(jīng)差增量ΔΩ，微分修正法的半長軸和傾角修正步長為

(17)

星座構(gòu)形數(shù)值微分修正方法的計算過程為：首先根據(jù)初始設(shè)計的軌道根數(shù)，加上高精度的動力學模型，做軌道預報，得出Δt時間后的ΔΩ和Δλ，由式(16)或(17)計算所需的星座衛(wèi)星間軌道根數(shù)調(diào)整量，更新星座衛(wèi)星初始軌道根數(shù)，再做軌道預報，看Δt時間后的ΔΩ和Δλ是否滿足要求，如此迭代循環(huán).該過程一般執(zhí)行4至8次就可獲得滿意的結(jié)果.

同半長軸和傾角都做修正相比，僅做半長軸修正能更多地節(jié)省推進劑的消耗，而從第4節(jié)的算例可見，星座構(gòu)形的漂移主要由相位差的漂移引起，而僅作半長軸修正能較好地對星座衛(wèi)星間相位差的漂移進行控制，兼顧了構(gòu)形保持與節(jié)省推進劑消耗量之間的平衡.

4 算 例

4.1 星座構(gòu)形攝動補償法算例

以構(gòu)形為Walker24/3/1全球星座為例進行星座控制策略的仿真計算，取軌道高度為21528km，傾角為55°，回歸周期為7天，一個回歸周期內(nèi)回歸13圈.取星座不同軌道面上的兩個衛(wèi)星軌道進行試算，初始軌道平根數(shù)如下所示.歷元時刻為2015年9月21日10點30分0.0秒.設(shè)副星和主星的升交點赤經(jīng)差-120°，相位差45°.

如果不做星座的保持控制，則星座10年后的構(gòu)形漂移如表1所示.

表1 不做控制的星座構(gòu)形漂移值

可見如果不做星座的保持控制，則10年后星座兩星的升交點赤經(jīng)差漂移量達到3.9°，相位差漂移達到9.34°.按照星座構(gòu)形攝動補償法計算的副星控制策略如表2所示.按照構(gòu)形攝動補償法計算出補償量后，進行平均軌道根數(shù)的預報，考慮地球非球形攝動項J2、J3、J4，日月引力，計算補償后的Δi、ΔΩ、Δλ.

表2 副星控制策略

可見，進行構(gòu)形攝動補償后，10年后兩星的升交點赤經(jīng)差漂移量為3.48°，相位差漂移達到3.58°，對相位差的漂移補償較好.對應(yīng)該攝動補償控制策略的星座構(gòu)形漂移如表3所示.

表3 采用構(gòu)形攝動補償控制策略后的星座構(gòu)形漂移值

由表3可見，進行構(gòu)形攝動補償控制后，升交點赤經(jīng)差的漂移方向不變，而相位差的漂移方向會反向，即相位差先逐漸減小，再逐漸增大，5年后相位差漂移達到-7.16°.

4.2 星座構(gòu)形數(shù)值微分修正法算例

仍以構(gòu)形為Walker24/3/1的全球星座為例，計算使用星座構(gòu)形數(shù)值微分修正方法的控制策略.為了對比，我們?nèi)允褂蒙弦恍」?jié)的軌道根數(shù)，將其看作瞬根.歷元時刻為2015年9月21日10點30分0.0秒.初始軌道根數(shù)如下所示：

如果不做星座的保持控制，考慮10×10階地球引力場、日月引力、太陽光壓.衛(wèi)星質(zhì)量均取為1900kg，截面積50m2，則星座的構(gòu)形漂移如表4所示.

表4 不做控制的星座構(gòu)形漂移值

由表4可知，星座衛(wèi)星主星和副星的相位差漂移速率較大.按照星座構(gòu)形數(shù)值微分修正法計算控制策略，如表5所示，由于星座構(gòu)形數(shù)值微分修正法直接使用瞬根進行計算，所以避免了瞬平轉(zhuǎn)換的誤差，因此，我們將控制周期取為5年進行試算.計算結(jié)果如表5和表6所示.

表5 副星控制策略

表6 采用構(gòu)形數(shù)值微分修正控制策略后的星座構(gòu)形漂移值

由表5和表6可知，星座構(gòu)形數(shù)值微分修正法對升交點赤經(jīng)差和相位差的控制效果極好.升交點赤經(jīng)差的漂移方向不變，而相位差的漂移方向會反向，即相位差先逐漸減小，再逐漸增大，5年后升交點赤經(jīng)差漂移1.76°，相位差漂移3.469°.

5 結(jié) 論

通過研究星座衛(wèi)星的軌道攝動及構(gòu)形相對漂移，研究了星座構(gòu)形長期保持的控制方法及控制策略.對于星座構(gòu)形的保持，研究了星座構(gòu)形攝動補償法，提出了星座構(gòu)形數(shù)值微分修正法.星座構(gòu)形攝動補償法是基于平均軌道根數(shù)的，對相位差漂移的補償效果較好.星座構(gòu)形數(shù)值微分修正法是基于瞬根數(shù)的，在實際應(yīng)用中，避免了瞬平轉(zhuǎn)換的誤差，對相位差的漂移控制效果較好.對比兩種方法，從實際應(yīng)用的角度和計算精確程度來講，星座構(gòu)形數(shù)值微分修正法略好于星座構(gòu)形攝動補償法.

參 考 文 獻

[1]潘科炎,王旭東,李果.星座與星座軌道控制技術(shù)[J].航天控制,2002,20(3)： 51-57

Pan K Y, Wang X D, Li G. Constellations and their orbit control techniques[J].Aerospace Control,2002,20(3):51-57

[2]Walker J G. Satellite constellations [J]. Journal of the British Interplanetary Society, 1984, 37: 559-572

[3]胡松杰,陳力,劉林.衛(wèi)星星座的結(jié)構(gòu)演化[J] .天文學報，2003,44(1): 46-53

Hu S J, Chen L, Liu L. The structure evolution of satellite constellation [J]. Acta Astronomica Sinica, 2003, 44(1): 46-53

[4]劉林.航天器軌道理論[M] .北京：國防工業(yè)出版社,2000

[5]項軍華,張育林.地球非球形對衛(wèi)星軌道的長期影響及補償研究[J].飛行力學，2007,25(2)：85-88

Xiang J H, Zhang Y L. Study on secular effect caused by earth non-spherical and compensation for its effect [J]. Flight Dynamics, 2007, 25(2)：85-88

[6]項軍華,張育林.區(qū)域覆蓋星座構(gòu)型優(yōu)化及協(xié)同控制策略研究[J].飛行力學，2007,25(3)：87-91

Xiang J H, Zhang Y L. Study on configuration optimization and coordination control strategy for the regional coverage satellite constellation [J]. Flight Dynamics, 2007, 25(3)：87-91

[7]項軍華,范麗,張育林.衛(wèi)星星座結(jié)構(gòu)自穩(wěn)定設(shè)計研究[J] .飛行力學，2007,25(4)：81-85

Xiang J H, Fan L, Zhang Y L. Study on design of the structure self-stabilization for satellite constellation [J]. Flight Dynamics, 2007, 25(4)：81-85

[8]李恒年,李濟生,焦文海.全球星攝動運動及攝動補償運控策略研究[J].宇航學報， 2010，31(7)：1756-1761

Li H N, Li J S, Jiao W H. Analyzing perturbation motion and studying configuration maintenance strategy for compass-m navigation constellation [J]. Journa1 of Astronautics, 2010，31(7)：1756-1761

[9]王瑞,向開恒,馬興瑞.平均軌道要素及其在衛(wèi)星星座設(shè)計中的應(yīng)用[J] .中國空間科學技術(shù),2002(5)： 14-20

Wang R, Xiang K H, Ma X R. The mean orbit elements and its application in the satellite constellation design [J]. Chinese Space Science and Technology，2002(5)： 14-20

[10]蒙波,葉立軍,韓潮.衛(wèi)星星座組網(wǎng)過程的策略規(guī)劃[J] .宇航學報，2009,30(1)：150-154

Meng B, Ye L J, Han C. Programming of satellite constellation establishment strategy [J]. Journa1 of Astronautics, 2009, 30(1)：150-154