劉世興，史東華，陳菊，易中貴，王本亮，龔自正，郭永新*

(1.遼寧大學(xué)物理學(xué)院，沈陽 110036；2.遼寧大學(xué)空間科學(xué)與技術(shù)研究院，沈陽 110036；3.北京理工大學(xué)數(shù)學(xué)與統(tǒng)計學(xué)院，北京 100081；4.北京理工大學(xué)宇航學(xué)院，北京 100081；5.北京衛(wèi)星環(huán)境工程研究所，北京 100094)

1 引言

“杞人憂天”出自《列子·天瑞》：杞國有人憂天地崩墜，身亡所寄，廢寢食者。雖然這里的“杞人”是被人嘲諷的對象，但實際上，小行星撞擊地球的事件在歷史上時有發(fā)生。夜空劃過的美麗流星，每年11月份出現(xiàn)的獅子座流星雨等都是流星顆粒進入地球大氣層和空氣摩擦而形成的。我們將圍繞太陽運行、尺寸在1m～800km內(nèi)，且不容易釋放出氣體和塵埃的天體，稱作小行星。太陽系中存在大量的小行星，主要分布在火星與木星軌道之間的小行星帶和海王星外的柯伊伯帶，如圖1所示。在天文學(xué)上，將距離地球最小距離在0.3AU(1AU=1.496×108km)范圍內(nèi)的小行星稱為近地小行星。地球曾發(fā)生過22次不同程度的生物滅絕事件，其中至少10次以上是由小行星撞擊地球所引起[1]，如6500萬年前的恐龍滅絕事件，1908年的通古斯空爆事件，2013年的俄羅斯車?yán)镅刨e斯克地區(qū)空爆事件等[2-5]。僅2021年內(nèi)全球已發(fā)生近地小行星飛掠地球事件1605次，觀測到29顆小行星進入地球大氣層發(fā)生火流星事件。我國境內(nèi)也經(jīng)常發(fā)生火流星事件，比如，2014年11月5日在內(nèi)蒙古錫林郭勒，2017年10月4日在云南香格里拉，2018年在云南西雙版納，2019年10月11日在吉林松原地區(qū)，2020年在青海玉樹地區(qū)，2021年在河南駐馬店等發(fā)生了火流星事件。在2022年3月11日，編號為2022 EB5，直徑大約3m的小行星在冰島北部上空進入了地球大氣層。雖然上述事件均未造成災(zāi)難，但已威脅到了人們的生產(chǎn)和生活安全。由此可見，近地小行星撞擊風(fēng)險概率雖小，但危害極大，但近地小行星撞擊地球的危害可測可防。防范近地小行星撞擊風(fēng)險，具有非常迫切的現(xiàn)實需求和深遠的戰(zhàn)略意義[6-10]。

圖1 太陽系中的八大行星和近地小行星帶分布圖

應(yīng)對近地小行星撞擊風(fēng)險，給國際天文界、航天界等多個領(lǐng)域帶來重大科學(xué)和技術(shù)挑戰(zhàn)。21世紀(jì)以來，國際社會高度重視小行星撞擊風(fēng)險防御工作，聯(lián)合國成立多個機構(gòu)專門負責(zé)協(xié)調(diào)應(yīng)對近地小行星撞擊風(fēng)險，各主要航天國家也積極致力于相關(guān)工作，并有力牽引、推動了空間科學(xué)研究、空間技術(shù)探索和空間資源利用。我國近地小行星撞擊風(fēng)險應(yīng)對工作起步較晚、基礎(chǔ)薄弱、國際貢獻率低、國際話語權(quán)小，這與我國航天大國地位不符[8]。近幾年我國已經(jīng)開始開展近地小行星防御問題研究，受到國家高度重視，并已上升為國家戰(zhàn)略[7]，且在小行星的軌道預(yù)測、領(lǐng)航監(jiān)測器的軌道設(shè)計與控制、小行星的在軌處置和監(jiān)測預(yù)警等方面取得了一定的研究成果[11-24]。應(yīng)對小行星撞擊風(fēng)險主要涉及兩個方面：(1)監(jiān)測預(yù)警。搜索發(fā)現(xiàn)、跟蹤定軌和數(shù)據(jù)更新、物性測量、計算撞擊概率、預(yù)估撞擊風(fēng)險走廊和落點，預(yù)報撞擊風(fēng)險。即發(fā)現(xiàn)小行星和確定撞擊風(fēng)險。(2)在軌處置。依據(jù)監(jiān)測預(yù)警信息，發(fā)射防御航天器，改變小行星軌道避免其撞擊地球(即偏轉(zhuǎn)，Deflection)，或?qū)⑿⌒行欠至殉伤槠苊饣蚪档推鋵Φ厍虻奈：?即摧毀，Disruption)。這里無論是小行星軌道精確預(yù)報、在大天體攝動影響下的軌道不確定性，還是在人為外力作用下(如動能撞擊等瞬時作用，或拖曳等長期作用下)的小行星軌道偏移、防御航天器自身的深空自主導(dǎo)航與控制技術(shù)等，都離不開對小行星運行軌道的精確計算測定，離不開精確的軌道設(shè)計與動力學(xué)研究。小行星的軌道動力學(xué)問題是防御的基礎(chǔ)、前提和核心問題。精確的軌道設(shè)計與計算需要采用基于約束系統(tǒng)的幾何動力學(xué)方法進行建模，同時采用高精度的幾何數(shù)值積分方法進行數(shù)值模擬。而軌道動力學(xué)一直是分析力學(xué)發(fā)展進程中的核心問題之一，但它針對航天工程提出的更復(fù)雜的軌道動力學(xué)問題研究的不夠深入，特別是航天工程中需要考慮到小行星的非規(guī)則形狀和自轉(zhuǎn)、引力源非球?qū)ΨQ性對牛頓引力場的修正、太陽引力場的相對論修正、小行星運行環(huán)境的隨機性等問題，給小行星軌道動力學(xué)建模、軌道高精度計算以及在軌處置航天器(群)的控制提出了更高要求，也是分析力學(xué)面臨的新課題。

2 幾何力學(xué)與幾何數(shù)值積分的發(fā)展現(xiàn)狀

現(xiàn)代分析力學(xué)最具標(biāo)志性的成果是幾何力學(xué)的發(fā)展，以及在此基礎(chǔ)上構(gòu)造的保結(jié)構(gòu)算法——幾何數(shù)值積分方法，將二者與控制結(jié)合起來所興起的計算幾何力學(xué)與控制使得若干動力學(xué)系統(tǒng)實現(xiàn)了長時間、高精度的數(shù)值建模分析，也將為小行星運行軌道的高精度數(shù)值模擬與測量、航天器(群)的精確軌道設(shè)計、動力學(xué)建模與控制提供一種理論基礎(chǔ)。同時，小行星軌道動力學(xué)建模和高精度計算也必將進一步完善分析力學(xué)的高精度算法，特別是Hamel場變分積分子算法。

幾何力學(xué)是基于流形、纖維叢、Lie群與Lie代數(shù)等現(xiàn)代微分幾何工具，將經(jīng)典分析力學(xué)升級成無坐標(biāo)語言的整體表述形式，形式更加簡潔優(yōu)美，為多體系統(tǒng)、流體、彈性力學(xué)、場論、量子系統(tǒng)及幾何控制理論等提供了統(tǒng)一的框架，它在空間探索的軌道設(shè)計、機器人運動控制、等離子體、超流、洋流動力學(xué)和氣候建模、計算解剖學(xué)等諸多領(lǐng)域已經(jīng)得到成功的應(yīng)用[25,26]。利用現(xiàn)代微分幾何語言表述Lagrange力學(xué)和Hamilton力學(xué)有以下幾點優(yōu)勢：首先，適用于位形流形上的一般力學(xué)系統(tǒng)，以及由纖維叢上聯(lián)絡(luò)所刻畫的約束力學(xué)系統(tǒng)[27,28]。聯(lián)絡(luò)可以自然地描述由非完整約束所決定的不可積分布，通過它將系統(tǒng)的基底變量和纖維變量聯(lián)系起來，并具有鮮明的物理意義，進而采用規(guī)范理論的觀點統(tǒng)一且有力地處理非完整力學(xué)和非線性控制理論中的諸多問題，包括力學(xué)系統(tǒng)的幾何表示和幾何結(jié)構(gòu)、穩(wěn)定性和鎮(zhèn)定、最優(yōu)控制等，在車輛、機器人等工程領(lǐng)域有重要的應(yīng)用；其次，方便用于描述系統(tǒng)的大范圍運動，也便于進行幾何類比而將有限維連續(xù)系統(tǒng)的結(jié)構(gòu)推廣至無窮維和離散力學(xué)系統(tǒng)[29,30]。特別地，對于位形空間是Lie群情況，采用對稱性約化理論[31]可實現(xiàn)代數(shù)空間上的線性化，還可以得到判別穩(wěn)定性的能量—動量方法。最后，基于離散Hamilton原理的變分積分子構(gòu)造保持力學(xué)系統(tǒng)幾何結(jié)構(gòu)的數(shù)值格式——幾何數(shù)值積分算法，為獲得長時間高精度的動力系統(tǒng)模擬提供有效的方法[32-36]。將計算幾何力學(xué)同幾何控制理論結(jié)合，得到能長時間保持高精度的數(shù)值控制算法，這一研究將為精確、高效、實時地控制航天器、潛水器、各類機器人等現(xiàn)代工程系統(tǒng)提供合理的設(shè)計原理。

目前，在無窮維系統(tǒng)中的幾何力學(xué)和幾何數(shù)值積分算法是國內(nèi)外相關(guān)領(lǐng)域的研究熱點，這既是分析力學(xué)研究的邏輯延伸，也是若干工程科學(xué)的實際需要，這其中Hamel形式及Hamel場變分積分子在動力學(xué)系統(tǒng)建模和高精度計算中發(fā)揮核心作用。采用無坐標(biāo)的幾何語言，可以通過類比，將具有對稱性的有限維力學(xué)系統(tǒng)的成熟方法推廣到無窮維系統(tǒng)，特別是借助于幾何類比，可將分析力學(xué)中不依賴于坐標(biāo)的活動標(biāo)架法推廣到無窮維系統(tǒng)和場論情況，通過直接對非物質(zhì)速度進行變分的方式構(gòu)造Hamel場變分積分子[36,37]，這個算法的能量保持顯著優(yōu)于Shake和Rattle算法的數(shù)值特性[38]；對基于SE(3)的幾何精確梁進行動力學(xué)分析，驗證了該算法能夠長時間保持幾何精確梁的能量與動量[39,40]；從協(xié)變場論觀點與時空離散角度研究Hamel場變分積分子算法的保動量性質(zhì)，在證明該算法保幾何結(jié)構(gòu)性質(zhì)的同時，間接說明該方法具有長時間數(shù)值穩(wěn)定的性質(zhì)[41]。此外，Hamel場變分積分子在計算無窮維非完整力學(xué)系統(tǒng)時具有明顯優(yōu)勢[42]，該方法表述Chaplygin-Timoshenko雪橇模型時表述簡單，便于分析該動力學(xué)方程解的適定性和解釋其特有的非線性行為。Hamel場變分積分子已用在工業(yè)機械臂、球面上移動機器人等的最優(yōu)控制，展現(xiàn)出了較好的數(shù)值特性[39,43,44]。

盡管Hamel場變分積分子算法在無窮維系統(tǒng)的動力學(xué)建模和高精度計算中得到了驗證，并且有理由相信不依賴于坐標(biāo)選取的Hamel形式適于描述不規(guī)則形狀和自轉(zhuǎn)的小行星軌道動力學(xué)和姿-軌耦合問題，但針對小行星的精確軌道計算和預(yù)測，還需要考慮各種實際非理想因素，如小行星不規(guī)則形狀、后牛頓效應(yīng)的攝動等對構(gòu)造Hamel變分積分子的影響。此外，小行星的運行環(huán)境是動態(tài)變化的，甚至是隨機的，盡管在有關(guān)隨機微分方程的處理上，幾何力學(xué)和幾何數(shù)值積分方法具有明顯的優(yōu)勢[45]，但是Hamel場變分積分子對隨機問題和基于學(xué)習(xí)的離散場論保結(jié)構(gòu)算法的拓展還缺乏研究。

3 近地小行星軌道動力學(xué)高精度建模研究進展

小行星在太陽系中的運動非常復(fù)雜，存在各種攝動因素的影響，如各大行星和小行星帶上質(zhì)量較大小行星的質(zhì)點引力攝動、太陽和大行星扁率產(chǎn)生的引力攝動、太陽引力場的相對論效應(yīng)、小行星自轉(zhuǎn)和形狀不規(guī)則以及表面熱物理環(huán)境的差異導(dǎo)致的亞爾科夫斯基(Yarkovsky)效應(yīng)所產(chǎn)生的非引力攝動項等。因此，小行星的運動不是簡單的Kepler問題，要精確計算和預(yù)測小行星運行軌道，一方面需要考慮諸多影響小行星運行的攝動因素，包括小行星在中心天體(太陽)和第三天體(大行星)的不規(guī)則形狀造成的引力攝動[46-49]。為精確計算小行星軌道，并考慮大天體(八大行星和月球)與主帶67顆質(zhì)量較大的小行星的質(zhì)點引力攝動、地球扁率項的引力攝動、太陽質(zhì)心的后牛頓效應(yīng)、太陽輻射壓對小行星軌道的攝動，通過引入測量隨機誤差、利用運動方程的積分計算、觀測量的理論歸算等求解小天體運行軌道[12, 14]。文獻[14]考慮各大行星的引力攝動、后牛頓效應(yīng)等因素，采用改進的顯式辛算法和嵌套的RKF7(8)積分器對43顆已命名的近地小行星軌道演化進行數(shù)值研究，在一定程度上彌補了對于大偏心率問題采用定步長的不足。除此之外，由于小行星自轉(zhuǎn)、形狀不規(guī)則以及表面熱物理環(huán)境的差異導(dǎo)致了Yarkovsky效應(yīng)也要產(chǎn)生非引力攝動項。另一方面在計算方法上，構(gòu)造包含非保守力項的高階Galerkin變分積分子能夠處理小天體姿-軌耦合、太陽輻射力等非保守力項的情況[50]，而使用Chebyshev控制點與Lagrange插值方法構(gòu)造的高階譜變分積分子比常規(guī)的Galerkin變分積分子具有更好的穩(wěn)定性[51]。采用變步長的數(shù)值算法，構(gòu)造精細的力學(xué)模型并減少計算機截斷誤差，以達到小行星軌道計算精度和效率要求[11, 13]。雖然以上成果對影響小行星運行軌道的各種攝動因素有所研究，但是針對不同的小行星，對這些引力攝動和非引力效應(yīng)如何對比觀測數(shù)據(jù)進行系統(tǒng)考慮，以提高計算精度和效率，還缺乏深入研究。現(xiàn)有的小行星運行軌道的計算方法難以在動態(tài)環(huán)境下實現(xiàn)長時間高精度的效果，缺乏對復(fù)雜攝動因素的系統(tǒng)分析。此外，小行星運行環(huán)境瞬息萬變，各種攝動要素的影響交替變化，具有一定的隨機性。為此需要采用分析力學(xué)的幾何數(shù)值積分方法，特別是Hamel場變分積分子算法，確保在引力和非引力攝動效應(yīng)下的長時間高精度計算需要加強幾何數(shù)值積分方法與隨機問題的結(jié)合研究。

4 近地小行星探測/在軌處置航天器(群)軌道動力學(xué)與控制研究進展

近地小行星具有弱引力場、形狀不規(guī)則的特點，從而使其外部引力場環(huán)境極其復(fù)雜，從而導(dǎo)致其運動形態(tài)復(fù)雜多樣。這些復(fù)雜性導(dǎo)致軌道優(yōu)化算法對初始值極其敏感，為探測/在軌處置航天器(群)的軌道設(shè)計和控制帶來了極大困難和挑戰(zhàn)[52-54]。

4.1 小行星處置航天器軌道設(shè)計與最優(yōu)控研究進展

目前，近地小行星處置航天器的軌道設(shè)計通常利用Lagrange點L1和L2 周圍的周期軌道的穩(wěn)定和不穩(wěn)定流形，揭示相空間中的管道網(wǎng)絡(luò)，這解釋了圓型限制性三體問題(CR3BP)中的遷移現(xiàn)象[55-59]。從不變流形理論出發(fā)，利用行星際高速公路來設(shè)計返回的轉(zhuǎn)移軌道[15]，基于三體問題的不變流形設(shè)計低成本登月軌道[16,17]。最近發(fā)現(xiàn)了太陽系中的一個從小行星帶延伸到天王星及更遠地方的流形管道結(jié)構(gòu)[60]；利用Lagrange擬序結(jié)構(gòu)(LCS)理論[61]計算橢圓型限制性三體問題(ER3BP)中的遷移勢壘，證明了ER3BP相空間中存在周期性脈動形式的Lagrange擬序結(jié)構(gòu)，同時限制性四體問題中也存在依賴時間的不變流形[19-20]。將不變流形理論和離散力學(xué)與最優(yōu)控制方法(DMOC)結(jié)合，成功計算出從地球到月球的具有脈沖推力的優(yōu)化軌跡[62]。文獻[18]綜述了利用不變流形理論和弱穩(wěn)定邊界理論研究三體系統(tǒng)中低能量轉(zhuǎn)移與捕獲軌道設(shè)計的研究進展。在航天器群重構(gòu)制導(dǎo)方面，基于笛卡爾坐標(biāo)的重構(gòu)技術(shù)，采用隨機樹搜索、遺傳算法、非線性優(yōu)化技術(shù)、序列二次規(guī)劃、相關(guān)狀態(tài)Riccati方程(SDRE)技術(shù)和腦風(fēng)暴優(yōu)化算法等，計算了跟隨航天器從給定初態(tài)到給定末態(tài)的最省能量控制軌跡[63-68]。

4.2 伴飛航天器軌道動力學(xué)的綜合建模研究進展

考慮伴飛航天器受太陽光壓、其它星體引力攝動、隨機擾動等綜合影響的高精度軌道動力學(xué)模型具有重要意義。已有工作主要考慮類地J2攝動而忽略其它的引力項對小行星弱引力場中的攝動軌道建立引力攝動模型[69]，也有僅考慮太陽光壓效應(yīng)而不考慮其它因素，建立繞小行星的軌道攝動模型，用于小行星航天器群編隊控制開發(fā)[70]。為實現(xiàn)從密切狀態(tài)到平均狀態(tài)的轉(zhuǎn)換，現(xiàn)有的轉(zhuǎn)換公式往往也是需要類地的關(guān)于引力場J2項的假設(shè), 利用一個基于J2冪展開的Hamilton量來轉(zhuǎn)換密切根數(shù)和平均根數(shù)，需假設(shè)J2比余緯向勢大幾個數(shù)量級，有的工作甚至忽略了除J2之外的任何攝動項[71]；對于中心天體的旋轉(zhuǎn)速率等于或大于軌道平均運動，使用頻率分解從相關(guān)信號中剔除密切(短周期)頻率的做法計算代價大，且需頻繁測量來檢測短周期頻率[72]。另一種思路是設(shè)計消除密切效應(yīng)的濾波器，利用Kalman濾波器，考慮二階引力勢從平均絕對根數(shù)中重構(gòu)密切絕對根數(shù)[73]。

4.3 探測/在軌處置航天器群編隊制導(dǎo)與控制研究進展

近地小行星防御需要航天器群在深空軌道高精度編隊飛行，與繞地軌道的編隊場景(星鏈(圖2)、GPS(圖3))相比，具有更大的挑戰(zhàn)性，描述其相對運動的方程更為復(fù)雜，必須考慮攝動影響，編隊高精度控制的連續(xù)小推力往往控制在納牛到微牛的量級[21]。要完成航天器群小行星防御任務(wù)需要解決兩個基本問題：編隊制導(dǎo)和編隊控制，前者指集群軌道設(shè)計問題，后者為集群在軌保持或重構(gòu)隊形。在編隊制導(dǎo)方面，常見的方法通常使用線性規(guī)劃、迭代學(xué)習(xí)控制、線性二次調(diào)節(jié)器(LQR)、相關(guān)狀態(tài)Riccati方程(SDRE)、Lagrange乘子、模型預(yù)測控制和感知框架等技術(shù)，通過跟蹤預(yù)定義的參考點實現(xiàn)[74-81]。另外，基于相對軌道根數(shù)的動力學(xué)模型是將編隊制導(dǎo)描述為跟蹤特定相對軌道根數(shù)的形式[82-84]。在編隊控制方面，可利用人工勢場來實現(xiàn)重構(gòu)[85,86]，也有使用多種生物優(yōu)化技術(shù)來確定最優(yōu)的重構(gòu)方案[87,88]，還有提出一種基于相對軌道根數(shù)的解析式燃料最優(yōu)脈沖編隊重構(gòu)策略[89]，以及基于相對軌道根數(shù)的ΔV最優(yōu)方案和相關(guān)軌道動力學(xué)的精確線性化方法[90,91]。繞地球運行的集群編隊控制方法，需要考慮大氣阻力和J2對相對運動的影響[82-92]。為避免集群中個體的相對漂移、內(nèi)部的碰撞，并減少整個任務(wù)中使用的推進劑，可以采用幾種新的開環(huán)控制方法[91]。另外，編隊保持方法也可針對特定軌道根數(shù)來設(shè)置[93,94]。

圖2 美國星鏈計劃[95]

圖3 GPS星座[96]

5 問題與結(jié)論

小行星防御涉及數(shù)學(xué)、物理、化學(xué)、力學(xué)和天體力學(xué)等多個學(xué)科領(lǐng)域，是一個復(fù)雜的系統(tǒng)工程，但其防御離不開四步：一來、二去、三繞、四控。小行星防御過程如圖4所示。

圖4 小行星防御過程示意圖

首先，所謂的“一來”，即如何對近地小行星軌道進行精確計算與預(yù)測問題。盡管國內(nèi)外對影響小行星運行軌道的各種攝動因素有所研究，包括各大行星和小行星帶上質(zhì)量較大小行星的質(zhì)點引力攝動、太陽和大行星扁率產(chǎn)生的引力攝動、太陽引力場的相對論效應(yīng)、小行星自轉(zhuǎn)和形狀不規(guī)則以及表面熱物理環(huán)境的差異導(dǎo)致的Yarkovsky效應(yīng)所產(chǎn)生的非引力攝動項等，但是針對不同的小行星，對這些引力攝動和非引力效應(yīng)如何對比觀測數(shù)據(jù)進行系統(tǒng)考慮，以提高計算精度和效率，還缺乏深入研究。同時，現(xiàn)有的小行星運行軌道的計算方法難以在動態(tài)環(huán)境下實現(xiàn)長時間高精度的效果，影響了復(fù)雜攝動因素的系統(tǒng)分析，為此需要采用分析力學(xué)的幾何數(shù)值積分方法，特別是Hamel場變分積分子算法，確保在引力和非引力攝動效應(yīng)下的長時間高精度計算。此外，小行星運行環(huán)境瞬息萬變，各種攝動要素的影響交替變化，具有一定的隨機性，需要加強幾何數(shù)值積分方法與隨機問題的結(jié)合研究。

其次，所謂的“二去”，即如何設(shè)計發(fā)射近地小行星處置航天器的軌道與實現(xiàn)最優(yōu)控制問題。根據(jù)時變系統(tǒng)的Lagrange擬序結(jié)構(gòu)理論，通過有限時間Lyapunov指數(shù)(FTLE)可以找到時變系統(tǒng)中的分界線，以此可以區(qū)分流場中有不同動力學(xué)特征的區(qū)域。然而，關(guān)于小行星防御航天器的軌道設(shè)計方面，多采用不變流形理論來設(shè)計軌道，缺乏采用Lagrange擬序結(jié)構(gòu)為最優(yōu)控制算法提供初始值的研究。另一方面，航天器群制導(dǎo)未利用領(lǐng)航器的優(yōu)化軌道，且基于笛卡爾坐標(biāo)的重構(gòu)技術(shù)計算量很大，此外，采用DMOC方法計算編隊飛行衛(wèi)星群重構(gòu)的最優(yōu)開環(huán)控制律，利用Halo軌道，限制了對小行星軌道的應(yīng)用。

再次，所謂的“三繞”，即如何對近地小行星伴飛航天器進行綜合建模問題。在構(gòu)建繞小行星的受擾Kepler軌道模型時，主要考慮類地球小行星的J2攝動建立引力攝動模型，無法達到任意小行星所需精度水平，也缺乏考慮相對運動中太陽光壓、其它第三星體引力攝動、隨機擾動等綜合影響，因而這樣的軌道動力學(xué)模型精度不高。此外，天體力學(xué)N體問題中存在一種規(guī)范對稱性，忽略其規(guī)范自由度可能造成軌道演化動力學(xué)在數(shù)值上的不穩(wěn)定性[72]。由于現(xiàn)有模型主要考慮類地球系統(tǒng)的引力場J2項，密切狀態(tài)到平均狀態(tài)轉(zhuǎn)換方法計算效率和計算精度較低，不適用于任意非地球小行星和其它行星體。有的剔除密切(短周期)頻率的頻率分解法具有計算代價大且需頻繁檢測的缺點，也不適于旋轉(zhuǎn)速度慢的小行星。即使采用Kalman濾波器，僅考慮低階引力場并不能準(zhǔn)確得到小行星中心軌道的密切或平均軌道根數(shù)。

最后，所謂的“四控”，即如何對探測/處置航天器群實現(xiàn)有效的編隊制導(dǎo)與控制，從而對具有潛在威脅的小行星實行干預(yù)或清除問題。目前，關(guān)于小行星防御航天器群隊形重構(gòu)和保持研究的多數(shù)技術(shù)是簡單地推廣近地編隊飛行的方法和成果，不能滿足集群任務(wù)關(guān)鍵要求。而且，航天器群制導(dǎo)未利用領(lǐng)航器的優(yōu)化軌道，且基于笛卡爾坐標(biāo)的重構(gòu)技術(shù)計算量很大，缺乏計算量少、能有效利用ΔV、基于相對軌道根數(shù)的大規(guī)模編隊重構(gòu)控制方案。當(dāng)前的一些方法只能考慮一個軌道周期的短時間跨度，動力學(xué)的線性化只在極短的時間內(nèi)是準(zhǔn)確的，而頻繁的線性化又會增加計算成本?；谙鄬壍栏鶖?shù)的編隊重構(gòu)控制方案通常局限于兩個航天器，并不能解決大群體的避碰問題。另外，缺乏解決在小行星特有的動態(tài)環(huán)境下航天器群全姿態(tài)最優(yōu)控制方法，忽略航天器姿態(tài)可能會影響任務(wù)的執(zhí)行。

綜上所述，目前針對小行星防御問題的軌道預(yù)測、偏移與設(shè)計、航天器群的編隊最優(yōu)制導(dǎo)與控制的研究還不完善，尚缺乏綜合考慮各種攝動效應(yīng)的精確的數(shù)值計算方法和合理的小行星動力學(xué)模型，以及帶隨機擾動的在軌航天器(群)全位姿動力學(xué)模型；同時，圍繞地球以外小行星特有動態(tài)環(huán)境下的航天器群編隊最優(yōu)制導(dǎo)與控制方法仍需進一步發(fā)展。上述問題的解決可為小行星的軌道偏轉(zhuǎn)、天基監(jiān)測系統(tǒng)的建立、引力拖車等編隊任務(wù)提供理論指導(dǎo)和技術(shù)支持。鑒于此，采用基于幾何力學(xué)和幾何數(shù)值積分建立動態(tài)環(huán)境下的高精度Hamel變分積分子，通過高精度積分子研究長期作用力下小行星軌道的演化規(guī)律，對小行星的軌道偏轉(zhuǎn)問題進行動力學(xué)建模，實現(xiàn)移步易景，實時計算小行星的運動狀態(tài)，做到對小行星偏轉(zhuǎn)軌道的精準(zhǔn)預(yù)測，給小行星的災(zāi)害防御提供精準(zhǔn)的動態(tài)數(shù)據(jù)。在此基礎(chǔ)上，構(gòu)建處置航天器基于Lagrange擬序結(jié)構(gòu)的最優(yōu)控制算法；建立考慮各種攝動的小行星在軌航天器全位姿動力學(xué)模型；針對航天器群編隊的構(gòu)建、保持和重構(gòu)等任務(wù)，研究其場論框架下柔性集群編隊控制理論和算法，為小行星防御的在軌處置提供理論支撐。