杜耀珂，楊盛慶，完備，王文妍，陳筠力

1. 上海航天控制技術(shù)研究所，上海 201109 2. 上海市空間智能自主控制技術(shù)重點(diǎn)實(shí)驗(yàn)室，上海 201109 3. 上海航天技術(shù)研究院，上海 201109

近地衛(wèi)星的軌道保持控制的目的是修正空間攝動(dòng)力對(duì)衛(wèi)星軌道參數(shù)的長(zhǎng)期累積影響，使衛(wèi)星保持在所需的目標(biāo)軌道附近。常規(guī)的近地軌道控制任務(wù)主要包括：衛(wèi)星太陽同步軌道的地面軌跡維持控制、大氣阻力對(duì)衛(wèi)星高度造成的衰減修正控制等。保持控制要求的精度較低，控制的時(shí)間間隔一般在幾個(gè)月甚至更長(zhǎng)的時(shí)間。因此，主要采用軌道參數(shù)超調(diào)的控制策略[1-2]。由于受到測(cè)量精度限制，以及保持控制殘差在幾個(gè)月甚至是幾年的長(zhǎng)時(shí)間累積影響，往往會(huì)造成衛(wèi)星實(shí)際的運(yùn)動(dòng)軌跡與期望的目標(biāo)軌道存在較大偏差。一般近地衛(wèi)星的軌道保持控制對(duì)回歸的精度要求并不太高，存在幾千米到幾十千米的星下點(diǎn)軌跡偏差，并不影響到衛(wèi)星的飛行任務(wù)。

隨著衛(wèi)星載荷對(duì)衛(wèi)星軌道精度要求的提高，以往的軌道設(shè)計(jì)和保持控制已經(jīng)不能滿足高精度對(duì)地觀測(cè)衛(wèi)星的使用需求。例如對(duì)于長(zhǎng)期地球形變或者地質(zhì)災(zāi)害監(jiān)測(cè)的重軌測(cè)量任務(wù)，就需要在整個(gè)壽命期間將衛(wèi)星軌道控制在一條預(yù)先設(shè)定的嚴(yán)格回歸軌道附近運(yùn)動(dòng)，偏差最大約為幾百米，這樣的軌道控制就稱為嚴(yán)格回歸軌道控制[3-4]。文獻(xiàn)[5-8]給出了嚴(yán)格回歸軌道的定義及設(shè)計(jì)方法：在獲取了嚴(yán)格回歸軌道起始點(diǎn)參數(shù)之后，以地球引力場(chǎng)模型為基礎(chǔ)進(jìn)行高精度軌道遞推，就能夠獲取目標(biāo)參考軌跡。文獻(xiàn)[9-11]描述了美國(guó)國(guó)家航空航天局(NASA)的LightSAR衛(wèi)星和日本的ALOS-2在自主高精度回歸軌道保持控制方面的工作，結(jié)果能夠有效提升近地衛(wèi)星軌道保持的精度。

本文在綜合了參考文獻(xiàn)研究的基礎(chǔ)上，將衛(wèi)星編隊(duì)的運(yùn)動(dòng)學(xué)理論[12-13]引入嚴(yán)格回歸軌道的分析和控制過程，提出了虛擬編隊(duì)構(gòu)形的概念。在此基礎(chǔ)上，分析實(shí)際衛(wèi)星相對(duì)于空間參考軌跡的偏差情況，并采取最優(yōu)的長(zhǎng)期軌道保持控制策略。本文從衛(wèi)星相對(duì)于參考軌跡的相對(duì)運(yùn)動(dòng)關(guān)系和軌道參數(shù)長(zhǎng)期累積變化兩個(gè)方面進(jìn)行分析和仿真。一方面相對(duì)運(yùn)動(dòng)關(guān)系反映了實(shí)際衛(wèi)星偏離空間參考軌跡的偏差狀況，能夠確定衛(wèi)星是否在軌道管道內(nèi)，為星上是否需要進(jìn)行自主保持控制提供依據(jù)；而另一方面，累積的空間攝動(dòng)影響通過虛擬編隊(duì)構(gòu)形參數(shù)變化分析出來，可以轉(zhuǎn)換為精確的衛(wèi)星軌道平根數(shù)偏差，從而計(jì)算出嚴(yán)格回歸軌道保持控制的控制量。

1 嚴(yán)格回歸軌道保持控制需求

嚴(yán)格回歸軌道是針對(duì)近年來的衛(wèi)星重軌對(duì)地觀測(cè)任務(wù)提出的一種具有新特點(diǎn)的衛(wèi)星運(yùn)行參考軌道。嚴(yán)格回歸軌道在設(shè)計(jì)上的特點(diǎn)是始末狀態(tài)的閉合性。經(jīng)過迭代設(shè)計(jì)得到的嚴(yán)格回歸軌道在空間的運(yùn)行軌跡如圖1中曲線所示。該軌跡是在地固坐標(biāo)系下一組完全確定的參考點(diǎn)的集合。嚴(yán)格回歸軌道的保持控制的目標(biāo)是衛(wèi)星始終運(yùn)行在以參考軌跡為中心，管道半徑為rt的管道內(nèi)，如圖2所示，因此也可以將嚴(yán)格回歸軌道的保持控制稱為“管道保持控制”。

圖1 衛(wèi)星嚴(yán)格回歸參考軌道設(shè)計(jì)結(jié)果Fig.1 Design results of satellite strictly-regressive orbit

圖2 管道保持控制需求Fig.2 Requirement of tube maintenance control

1.1 定義與實(shí)現(xiàn)

衛(wèi)星嚴(yán)格回歸軌道是指在僅考慮地球重力場(chǎng)模型情況下，在一個(gè)回歸周期的時(shí)間段內(nèi)，地固坐標(biāo)系下空間軌跡的初始點(diǎn)和末端點(diǎn)的位置和速度“完全”相同。

在上述定義中，地球重力場(chǎng)模型一般是基于特定的重力場(chǎng)模型文件，如：EGM2008?！巴耆毕嗤?，是指初末狀態(tài)的偏差限定在一個(gè)很小的范圍內(nèi)(通常位置誤差小于1 m)。根據(jù)以上條件，通過迭代計(jì)算，就可以得到一組WGS84坐標(biāo)系下時(shí)間、位置和速度的序列，該序列就稱為空間參考軌跡。在給定的時(shí)刻點(diǎn)情況下，對(duì)于參考軌跡中的一組位置和速度參數(shù)，經(jīng)過坐標(biāo)變換能夠得到嚴(yán)格回歸軌道的軌道參數(shù)。

1.2 高階軌道動(dòng)力學(xué)與模擬

為了能夠盡量真實(shí)地模擬地球重力場(chǎng)對(duì)衛(wèi)星軌道的影響，應(yīng)采用地球引力場(chǎng)位函數(shù)的表示方法，便于將地球重力場(chǎng)文件進(jìn)行衛(wèi)星軌道積分。在地固系中，一個(gè)質(zhì)點(diǎn)所承受的地球引力加速度為

(1)

式中：R=R(X,Y,Z)為質(zhì)點(diǎn)的地心位置矢量，R為R的矢量模值；V(R)為地球引力場(chǎng)位函數(shù)；μ為地球引力常數(shù)。

參考文獻(xiàn)[14]，可以得到如下以勒讓德多項(xiàng)式的表示形式

(2)

嚴(yán)格回歸參考軌道設(shè)計(jì)過程僅考慮地球重力場(chǎng)影響，而衛(wèi)星實(shí)際軌道不僅受地球重力場(chǎng)影響，還包括大氣阻力、日月三體引力、太陽光壓、固體潮、相對(duì)論效應(yīng)等，因此需要將上述攝動(dòng)力進(jìn)行建模，并仿真得到對(duì)應(yīng)時(shí)刻J2000坐標(biāo)系下的加速度，從而引入仿真模型。具體過程可以參考文獻(xiàn)[15-16]。

1.3 虛擬衛(wèi)星編隊(duì)及其坐標(biāo)系定義

本文的研究和分析過程采用的坐標(biāo)系為OIRITINI和OERETENE，其坐標(biāo)系定義如圖3所示，其中ωe為地球自轉(zhuǎn)角速度，v為參考衛(wèi)星的飛行速度。兩種坐標(biāo)的差異主要是地球自轉(zhuǎn)對(duì)坐標(biāo)系原點(diǎn)參數(shù)影響和衛(wèi)星飛行方向在兩個(gè)坐標(biāo)系投影的不同。OERETENE坐標(biāo)系用于分析相對(duì)參考軌跡的狀態(tài)，OIRITINI坐標(biāo)系用于分析衛(wèi)星相對(duì)軌道根數(shù)的變化。

圖中：取參考軌跡點(diǎn)對(duì)應(yīng)的慣性系位置坐標(biāo)為原點(diǎn)OI，OIRI為從地心指向OI，OINI指向軌道面法線方向，OITI根據(jù)右手法則確定。該坐標(biāo)系原點(diǎn)為嚴(yán)格回歸軌道的參考點(diǎn)，并非衛(wèi)星的真實(shí)坐標(biāo)，結(jié)合文獻(xiàn)[12]中編隊(duì)坐標(biāo)系的定義，稱該坐標(biāo)系為“虛擬編隊(duì)坐標(biāo)系”。

OERETENE為基于地固系坐標(biāo)的直角坐標(biāo)系，稱為“軌跡坐標(biāo)系”。取參考軌跡點(diǎn)對(duì)應(yīng)的地固系位置坐標(biāo)為原點(diǎn)OE，OERE為從地心指向OE，OETE垂直于OERE指向后續(xù)的參考軌跡點(diǎn)，OENE根據(jù)右手法則確定。

由慣性系到虛擬編隊(duì)坐標(biāo)系的轉(zhuǎn)換關(guān)系為

(3)

式中：rI、vI分別為慣性系位置和速度矢量。

圖3 OIRITINI與OERETENE的坐標(biāo)系定義Fig.3 Definition of OIRITINI and OERETENE coordination

由地固系到軌跡坐標(biāo)系的轉(zhuǎn)換關(guān)系為

(4)

式中：rE和vE分別為地固系位置和速度矢量。

衛(wèi)星的嚴(yán)格回歸軌道控制是指對(duì)于給定的參考軌跡上的任意一個(gè)檢驗(yàn)點(diǎn)k，其空間誤差Ek定義為：參考軌跡的檢驗(yàn)點(diǎn)k到真實(shí)軌跡上一點(diǎn)Pa的矢量在檢驗(yàn)點(diǎn)k的OT向(切向)分量為零，檢驗(yàn)點(diǎn)k到真實(shí)軌跡上一點(diǎn)Pa的距離定義為Ek。若對(duì)于參考軌跡上的所有Ek均小于一個(gè)正常數(shù)Rt，則稱真實(shí)軌跡位于參考軌跡附近半徑為Rt的管道內(nèi)。以上即是衛(wèi)星嚴(yán)格回歸軌道保持控制的基本任務(wù)，稱為“管道保持控制”。

2 衛(wèi)星軌道偏離參考軌跡分析

基于上述嚴(yán)格回歸軌道和坐標(biāo)系定義，取軌道參數(shù)初值如表1所示。其中參考軌道和衛(wèi)星軌道初值均采用表中所示參數(shù)，但在軌道積分模型中，空間參考軌跡僅采用70階地球重力場(chǎng)模型文件，而真實(shí)衛(wèi)星除相同的地球重力場(chǎng)模型文件之外，還加入了其他空間攝動(dòng)力模型。

軌道參數(shù)用衛(wèi)星相對(duì)于參考軌跡點(diǎn)的軌道參數(shù)差表示，選取軌道半長(zhǎng)軸偏差Δa，升交點(diǎn)赤經(jīng)差ΔΩ，軌道傾角差Δi和緯度輻角差Δu來表示，為了便于比較，將ΔΩ、Δi和Δu與軌道半長(zhǎng)軸a相乘， 轉(zhuǎn)換為長(zhǎng)度單位進(jìn)行表示。通過衛(wèi)星軌道在OIRITINI和OERETENE坐標(biāo)系內(nèi)的投影，確定衛(wèi)星相對(duì)于參考軌跡的累積狀態(tài)，以確定是否超出以參考軌道為中心的“管徑”范圍。

表1 初始軌道平根數(shù)[17]Table 1 Mean elements of initial orbit [17]

2.1 大氣阻力影響

使用太陽活動(dòng)高年的參數(shù)，衛(wèi)星面質(zhì)比取為0.01，經(jīng)過仿真，大氣阻力引起衛(wèi)星軌道半長(zhǎng)軸每天衰減約15 m，仿真時(shí)間長(zhǎng)度t為3天。經(jīng)仿真，在此期間，Δa的變化很明顯，衛(wèi)星緯度輻角u的累積變化a·Δu將達(dá)到約10 km，a·ΔΩ和a·Δi變化很小。以上衛(wèi)星軌道參數(shù)的變化轉(zhuǎn)換到軌跡坐標(biāo)系下，相對(duì)于參考軌道中心的偏差達(dá)到約600 m。Δa的模值越大，Δu累積越快，隨著時(shí)間的累積，造成的管徑偏差還將迅速增加。衛(wèi)星軌道受大氣阻力影響偏離參考軌道的情況如圖4 所示，對(duì)應(yīng)偏離嚴(yán)格回歸軌道在OIRINI平面的投影如圖5(a)所示，在此坐標(biāo)系下可以顯示衛(wèi)星虛擬編隊(duì)構(gòu)形參數(shù)變化。在OERENE平面的投影如圖5(b)所示，在此坐標(biāo)系下可以顯示相對(duì)于空間參考軌跡的管道半徑偏差，簡(jiǎn)稱管徑偏差。

2.2 日月引力攝動(dòng)影響

日月引力為保守力，量級(jí)較小，經(jīng)過模擬衛(wèi)星1年的運(yùn)行情況，ΔΩ、Δi存在一定的周期性波動(dòng)，且a·ΔΩ的波動(dòng)范圍達(dá)到了約2 800 m，a·Δi的波動(dòng)既包含了約14天為周期的短期波動(dòng)項(xiàng)，也包括了以約半年為周期的長(zhǎng)期波動(dòng)項(xiàng)，波動(dòng)范圍到達(dá)了約300 m。但半長(zhǎng)軸無持續(xù)衰減情況，僅在小范圍內(nèi)波動(dòng)，因此緯度輻角的累積變化Δu較小。相對(duì)于空間參考軌道中心的偏差達(dá)到約2 300 m。衛(wèi)星軌道受日月攝動(dòng)影響偏離參考軌道的情況如圖6所示，對(duì)應(yīng)的在OIRINI平面的投影如圖7(a)所示，在OERENE平面的投影如圖7(b)所示。

圖4 衛(wèi)星受大氣阻力影響偏離參考軌道的情況Fig.4 Deviation of satellite’s orbit to reference orbit due to atmospheric drag

圖5 衛(wèi)星軌跡受大氣阻力影響下在OIRINI與OERENE平面的投影(3天)Fig.5 Projection of satellite’s track in OIRINI plane and OERENE plane due to atmospheric drag (3 days)

圖6 衛(wèi)星受日月攝動(dòng)影響偏離參考軌道的情況Fig.6 Deviation of satellite’s orbit relative to reference orbit due to sun and moon’s perturbation

圖7 衛(wèi)星軌跡在受日月攝動(dòng)影響在OIRINI與OERENE平面上的投影(365天)Fig.7 Projection of satellite’s orbit in OIRINI plane and OERENE plane due to sun and moon’s perturbation (365 days)

2.3 太陽光壓、固體潮、相對(duì)論效應(yīng)等攝動(dòng)影響

太陽光壓、固體潮和相對(duì)論效應(yīng)等空間攝動(dòng)雖然包含了太陽光壓這種非保守力的長(zhǎng)期影響，但由于攝動(dòng)的量級(jí)總體較小。ΔΩ、Δi存在一定的周期性波動(dòng)，且a·ΔΩ的波動(dòng)范圍達(dá)到了約500 m，a·Δi的波動(dòng)既包含了約14天為周期短期的波動(dòng)項(xiàng)，也包括了以約半年為周期的長(zhǎng)期波動(dòng)項(xiàng)，波動(dòng)范圍到達(dá)了約50 m。半長(zhǎng)軸僅在小范圍內(nèi)波動(dòng)，量級(jí)為米級(jí)，從模擬衛(wèi)星運(yùn)行1年的結(jié)果來看，衛(wèi)星偏離嚴(yán)格回歸軌道中心較小，最大約為750 m。衛(wèi)星軌道受太陽光壓、固體潮和相對(duì)論效應(yīng)攝動(dòng)影響的變化如圖8所示，對(duì)應(yīng)的在OIRINI平面的投影如圖9(a)所示，在OERENE平面的投影如圖9(b)所示。

2.4 攝動(dòng)影響小結(jié)

將本節(jié)攝動(dòng)仿真結(jié)果梳理如表2所示。

3 管道保持控制設(shè)計(jì)

3.1 管道偏差來源分析

衛(wèi)星偏離嚴(yán)格回歸軌道中心的管徑偏差包括法向(OENE)分量EN和徑向(OERE)分量ER兩個(gè)部分。EN的影響因素是衛(wèi)星軌道與嚴(yán)格回歸軌道半長(zhǎng)軸偏差Δa引起的累積影響，即隨著時(shí)間累積Δu偏差逐漸增加，隨著地球自轉(zhuǎn)，該偏差會(huì)轉(zhuǎn)移到軌跡坐標(biāo)系的法線方向。根據(jù)衛(wèi)星軌道運(yùn)動(dòng)學(xué)，EN與Δu的關(guān)系可以表示為

圖8 衛(wèi)星受太陽光壓、固體潮和相對(duì)論效應(yīng)攝動(dòng)影響偏離參考軌道的情況Fig.8 Deviation of satellite’s orbit relative to reference orbit due to sunray pressure, solid tide and relativistic effect

圖9 衛(wèi)星軌跡受太陽光壓、固體潮和相對(duì)論效應(yīng)攝動(dòng)影響在OIRINI與OERENE平面上的投影(365天)Fig.9 Projection of satellite’s track in OIRINI plane and OERENE plane due to sunray pressure, solid tide and relativistic effect (365 days)

表2 空間攝動(dòng)對(duì)衛(wèi)星軌道相對(duì)參考軌跡的影響匯總Table 2 Summery of deviation of satellite’s orbit relative to reference track under space perturbations

PerturbationsRelative orbit elements effectsTube errorCharacteristicsAtmosphere dragΔaaccumulated by time, max:-15 m/day;Δu accumulated by time andΔa, Δu increase faster as Δa growsMax: 600m in 3 daysAccumulated quickly, periodi-cally maintenance control nee-dedThree body gravityΔi,ΔΩ change periodically and constantly; a·Δi: 300 m, with short-term range 150 m, long-term range 100 m; a·ΔΩ:2 800 mMax: 2 300 m/aOnly the constantly change nee-ded control when it exceed the range demandsSun light pressure, Solid tide, Relativity effectΔaaccumulated by time, very small;Δu accumulated slow-ly;Δi,ΔΩ change periodically;a·Δi:about 50 m;a·ΔΩ:about 500 mMax: 750 m/aVery small,could be controlled when the front two perturba-tions is corrected

(5)

式中：a為參考軌跡對(duì)應(yīng)的軌道半長(zhǎng)軸；ωe為地球自轉(zhuǎn)角速度；n為參考軌跡對(duì)應(yīng)的軌道角速度；u為衛(wèi)星緯度輻角。

軌道根數(shù)控制殘差和受攝運(yùn)動(dòng)累積造成衛(wèi)星軌道與嚴(yán)格回歸軌道的軌道傾角差Δi，升交點(diǎn)赤經(jīng)差ΔΩ，也會(huì)造成EN的增加。其影響關(guān)系可以參照編隊(duì)衛(wèi)星相對(duì)運(yùn)動(dòng)的關(guān)系來進(jìn)行分析，影響大小可以表示為[18]

EN2(u)=a(Δisinu-ΔΩsinicosu)

(6)

式中：i為嚴(yán)格回歸軌道的軌道傾角，其他參數(shù)定義同前。

徑向偏差分量ER的影響因素包括實(shí)際衛(wèi)星軌道與參考軌道半長(zhǎng)軸偏差Δa以及衛(wèi)星偏心率矢量受攝運(yùn)動(dòng)與嚴(yán)格回歸軌道偏心率矢量的差異引起[18]。

ER(u)=-Δexcosu-Δeysinu

(7)

(8)

式中：e為參考軌跡對(duì)應(yīng)的軌道偏心率；ω為參考軌跡對(duì)應(yīng)的近地點(diǎn)輻角；es為實(shí)際衛(wèi)星的軌道偏心率；ωs為實(shí)際衛(wèi)星的近地點(diǎn)輻角。

由于衛(wèi)星實(shí)際軌道具備準(zhǔn)凍結(jié)特性，對(duì)主要攝動(dòng)源的考慮與參考軌跡設(shè)計(jì)考慮地球非球形影響比較一致，偏心率矢量演化是相對(duì)緩慢的。因此，引起管道半徑超出設(shè)定的管道邊界閾值的因素主要為管徑法向分量EN的變化。

3.2 嚴(yán)格回歸特性保持策略設(shè)計(jì)

根據(jù)以上運(yùn)動(dòng)學(xué)分析，由Δu造成的衛(wèi)星軌道偏離參考軌跡狀態(tài)在不同的緯度帶上的影響不同，在地球兩極附近影響最小，地球赤道附近影響最大。平面外軌道根數(shù)偏差造成的偏離參考軌跡狀態(tài)在低緯度主要由ΔΩ決定，在高緯度主要有Δi決定。因此，通過合理偏置衛(wèi)星軌道的緯度輻角ΔuP，就能夠能夠抵消ΔΩ偏差造成的偏離參考軌跡狀態(tài)。

考慮控制殘差以及大氣阻力、日月三體引力等攝動(dòng)影響，從節(jié)省燃料消耗以及盡量減少控制頻次需求出發(fā)，管道控制保持方案采用“逐次超調(diào)+ 閾值觸發(fā)”相結(jié)合的控制方式。

1) 超調(diào)與偏置

綜合以上分析，為了抵消大氣阻力對(duì)衛(wèi)星軌道半長(zhǎng)軸衰減的影響，應(yīng)采用半長(zhǎng)軸超調(diào)控制策略，這樣可以盡可能利用自然攝動(dòng)，減少軌道保持控制頻次。但是半長(zhǎng)軸超調(diào)量的確定還需要考慮到緯度輻角u的偏置需求。通過將u偏置到uP，能夠在低緯度地區(qū)抵消升交點(diǎn)赤經(jīng)偏差造成的管徑偏差。取衛(wèi)星偏離嚴(yán)格回歸軌道中心的偏差包括法向分量邊界值為ENmax，則由式(5)和式(6)可得緯度輻角的偏置量ΔuP可以表示為

(9)

式中：第1項(xiàng)為抵消升交點(diǎn)赤經(jīng)偏差影響的偏置量；第2項(xiàng)為考慮最大管徑偏差之后還能夠利用半長(zhǎng)軸衰減可以使用偏置量。

由于衛(wèi)星實(shí)際半長(zhǎng)軸相對(duì)于參考軌跡會(huì)不斷降低，因此為了減小控制頻次，在軌只進(jìn)行半長(zhǎng)軸正偏置。取當(dāng)前實(shí)際緯度輻角偏差為Δu

若Δu-ΔuP>0，則半長(zhǎng)軸超調(diào)量ΔaP應(yīng)采用式(10)計(jì)算：

(10)

否則，可以利用軌道參數(shù)自然衰減，實(shí)現(xiàn)偏離參考軌跡中心偏差的減小，不需要進(jìn)行主動(dòng)保持控制。

衛(wèi)星實(shí)際在軌運(yùn)行期間，半長(zhǎng)軸偏置量由實(shí)際半長(zhǎng)軸衰減速率、緯度輻角的可用偏置量、以及平面外參數(shù)軌道參數(shù)控制殘差與攝動(dòng)變化等因素綜合考慮后得到。

此外，由偏心率矢量攝動(dòng)造成的偏離參考軌跡偏差相對(duì)于軌道高度變化引起的偏差為一個(gè)小量，通過式(5)可以確定出偏心率Δe的控制量。在對(duì)半長(zhǎng)軸修正控制過程中，選擇合適的控制點(diǎn)達(dá)到半長(zhǎng)軸和偏心率聯(lián)合修正的目的。

2) 閾值觸發(fā)

由于控制殘差的存在以及日、月引力的影響，衛(wèi)星的軌道傾角i在一段時(shí)間存在單向變化規(guī)律以及極限邊界，需要進(jìn)行修正控制。由相對(duì)運(yùn)動(dòng)分析，Δi的單向變化逐漸壓縮逐次進(jìn)行半長(zhǎng)軸超調(diào)與緯度輻角偏置的作用范圍，削弱逐次超調(diào)與偏置的管道保持控制效果。

當(dāng)Δi超出特定閾值時(shí)，必須進(jìn)行平面外觸發(fā)控制。將軌道傾角的偏差閾值設(shè)定為大于250 m的參數(shù)值，就能夠達(dá)到忽略周期波動(dòng)項(xiàng)，僅通過逐次超調(diào)控制，保證衛(wèi)星運(yùn)行軌跡不超出管徑偏差的要求。

ΔΩ的受攝運(yùn)動(dòng)引起的管徑偏差變化范圍約為Δi影響的10倍，但是經(jīng)過前面的運(yùn)動(dòng)學(xué)分析，可以通過緯度輻角u的偏置來抵消，因此不需要主動(dòng)控制。

4 數(shù)值試驗(yàn)

取參考衛(wèi)星軌道高度約為743 km，大氣阻力引起衛(wèi)星軌道半長(zhǎng)軸每天衰減約15 m。參考軌跡僅用70階地球重力場(chǎng)模型進(jìn)行遞推仿真，時(shí)間衛(wèi)星起始時(shí)刻與參考軌跡起始點(diǎn)狀態(tài)一致，遞推仿真除考慮階地球重力場(chǎng)之外，還包括了大氣阻力，日月三體引力，太陽光壓，固體潮，相對(duì)論效應(yīng)等攝動(dòng)力影響。

控制策略采用脈沖軌道半長(zhǎng)軸和相位偏置控制方式，根據(jù)偏離參考軌跡的狀態(tài)，采用超差觸發(fā)方式，觸發(fā)控制策略的邊界設(shè)置為450 m，軌道傾角誤差的啟控邊界設(shè)置為400 m。仿真時(shí)間為1年， 在保持控制情況下，衛(wèi)星相對(duì)軌道根數(shù)偏差、相對(duì)參考軌跡的管徑偏差分布情況、投影曲線及統(tǒng)計(jì)結(jié)果如圖10～圖13所示。

根據(jù)嚴(yán)格回歸軌道保持控制的仿真結(jié)果，1年 之內(nèi)，共進(jìn)行了72次嚴(yán)格回歸軌道保持控制，其中軌道平面內(nèi)保持控制為70次，軌道平面外保持為2次，平均約5天進(jìn)行一次保持控制，需要的總速度增量為3.88 m/s。通過保持控制，衛(wèi)星在全部時(shí)間內(nèi)，處于參考軌跡為中心偏差500 m 范圍內(nèi)的概率優(yōu)于98.95%，而且經(jīng)過軌道相位偏置能夠避免對(duì)衛(wèi)星升交點(diǎn)赤經(jīng)偏差的修正控制，極大的節(jié)省了衛(wèi)星所需的推進(jìn)劑消耗(如果修正ΔΩ，所需增加的控制速度增量將約為當(dāng)前策略速度增量的5倍)。

圖10 相對(duì)軌道根數(shù)變化(366天)Fig.10 Variation of relative orbit elements (366 days)

圖11 衛(wèi)星管徑偏差分布Fig.11 Distribution of satellite tube radius error

圖12 衛(wèi)星管徑偏差在OENERE平面上的投影(366天)Fig.12 Projection of satellite tube radius error in OENERE plane (366 days)

圖13 衛(wèi)星滿足管道控制范圍統(tǒng)計(jì)Fig.13 Statistic of satellite meeting requirements of tube control

5 結(jié) 論

1) 本文描述了基于地固坐標(biāo)系參考軌跡的衛(wèi)星嚴(yán)格回歸保持軌道控制所涉及的動(dòng)力學(xué)與運(yùn)動(dòng)學(xué)分析，以及據(jù)此設(shè)計(jì)的管道保持控制問題的一種解決方法。

2) 管道保持控制通過半長(zhǎng)軸超調(diào)和相位偏置的設(shè)計(jì)思想，一方面盡可能利用自然攝動(dòng)減少對(duì)衛(wèi)星軌道的修正次數(shù)，另一方面避免了對(duì)衛(wèi)星升交點(diǎn)赤經(jīng)的頻繁修正控制，能夠有效減少對(duì)星載推擠劑的消耗。

3) 對(duì)于長(zhǎng)期累積的傾角攝動(dòng)，設(shè)計(jì)了閾值觸發(fā)控制方式，避免了該項(xiàng)周期性波動(dòng)可能造成的頻繁控制。

4) 長(zhǎng)時(shí)間仿真結(jié)果表明，該方法能夠確保衛(wèi)星運(yùn)行在設(shè)定的參考軌跡偏差管道內(nèi)。