雙二體幾何切面法及天梯地月轉(zhuǎn)移軌道分析*

童科偉，程炳琳，汪小衛(wèi)

(1. 中國運載火箭技術(shù)研究院 研究發(fā)展部， 北京 100076; 2. 中國運載火箭技術(shù)研究院， 北京 100076)

雙二體軌道拼接設(shè)計方法在月球探測軌道設(shè)計中具有基礎(chǔ)性作用，國內(nèi)外學者以雙二體軌道拼接法為基礎(chǔ)對月球探測軌道設(shè)計方法進行了大量的研究[1-8]。

在基于雙二體問題求解地月轉(zhuǎn)移軌道時，Battin、郗曉寧、黃文德、賀波勇等都不約而同地以入口點月心經(jīng)緯度作為中間變量來拼接地球影響球內(nèi)飛行軌道和月球影響球內(nèi)飛行軌道[1-2,4-5,8]。這種基于月球影響球球坐標來描述入口點位置的方法比較直觀，但是沒有充分利用轉(zhuǎn)移軌道的物理意義。

對于著陸地點已知的奔月任務(wù)，賀波勇推導了一種類解析的改進雙二體模型，該方法與傳統(tǒng)的從地球計算到月球不同，在計算窗口時針對已知的月心軌道六根數(shù)逆向計算地月轉(zhuǎn)移軌道六根數(shù)，計算量極小[9]。

對于工程約束不明、具體探測地點不定的地月轉(zhuǎn)移軌道窗口計算問題，還需要按照傳統(tǒng)的計算量較大的搜索方法，進而得到一些運動特性。與常規(guī)的以入口點月心經(jīng)緯度作為中間變量的雙二體拼接方法不同，本文提出一種直接基于軌道特征的地月轉(zhuǎn)移飛行軌道的幾何方法，物理意義明確，求解地月軌道更加直觀，通過兩次降維解耦操作，降低問題求解維數(shù)，進一步提升計算效率。

考慮到目前大多數(shù)計算機都是多核中央處理器(central processing unit, CPU)，最簡單的提高算力的途徑是采用并行計算，充分利用各核的計算能力，并行計算技術(shù)在國內(nèi)外受到高度重視。并行程序設(shè)計編程語言通?；谙鬟f接口標準(message passing interface, MPI)、OpenMP等，二者甚至已經(jīng)分別形成了集群并行和單機多核并行領(lǐng)域的事實標準[10-13]。

本文以Java相關(guān)的多核并行編程技術(shù)為基礎(chǔ)[13]，設(shè)計了一種地月轉(zhuǎn)移軌道的多核并行計算方法，并應(yīng)用于天梯地月轉(zhuǎn)移軌道計算與特性分析。按照國際天梯協(xié)會(international space elevator consortium, ISEC)的定義，天梯是一種把有效載荷從地球表面提升到空間的系統(tǒng)，該系統(tǒng)是一個長10萬千米的繩系，質(zhì)心在地球靜止軌道(geostationary orbit, GEO)，并固定于地表的某個錨點。隨著爬升高度的變化，天梯攀爬器在電能作用下抬升高度并獲得地球自轉(zhuǎn)帶來的圓速度，從而使有效載荷獲得勢能和動能[14]。由于天梯位于赤道上空，軌道傾角受限，與常規(guī)月球轉(zhuǎn)移軌道設(shè)計有一定區(qū)別。本文以一種幾何切面法雙二體模型來分析天梯發(fā)射航天器進入月球轉(zhuǎn)移軌道的特性。

1 雙二體幾何切面法

雙二體假設(shè)下的地月轉(zhuǎn)移軌道計算也稱為圓錐曲線拼接法。本質(zhì)上是以月球影響球為邊界, 將軌道分成多段來拼接：進入影響球前, 飛行軌跡為地心圓錐曲線(一般為橢圓)；在影響球內(nèi)，飛行軌跡為月心圓錐曲線(一般為雙曲線)。對于自由返回軌道等任務(wù)，還需考慮飛出影響球后的飛行軌跡，其為地心圓錐曲線(一般也為橢圓)。在進入和離開影響球的兩個邊界點(一般稱為入口點和出口點)，將這幾段圓錐曲線拼接成完整的運動軌跡。

圓錐曲線拼接法的設(shè)計參數(shù)通常以月球影響球入口點/出口點在月心白道坐標系下的經(jīng)緯度作為獨立變量，通過該經(jīng)緯度來描述入口點/出口點的空間位置[1-2,4-6,8]。

采用傳統(tǒng)的圓錐曲線拼接法的設(shè)計參數(shù)計算轉(zhuǎn)移窗口時問題待求解參數(shù)的維度大，計算轉(zhuǎn)移窗口時的計算時間過長。對于在天梯系統(tǒng)出發(fā)這一類地月轉(zhuǎn)移軌道設(shè)計問題，由于出發(fā)軌道的傾角已知、升交點赤經(jīng)(right ascension of ascending node, RAAN)也可求出，本文借助空間幾何關(guān)系把原始問題降維，從而大大降低求解問題的規(guī)模。其基本思想是：基于轉(zhuǎn)移軌道面和月球影響球的幾何關(guān)系來描述空間位置關(guān)系，地月/月地轉(zhuǎn)移軌道面與月球影響球必須相交(否則一定不可能被月球捕獲)。通過轉(zhuǎn)移軌道面切割月球影響球形成的幾何關(guān)系就可以描述地月/月地轉(zhuǎn)移軌道，結(jié)合Lambert原理就可以求解出所需的轉(zhuǎn)移軌道參數(shù)。這樣就把原來的三維球面搜索算法降維成了二維平面內(nèi)切割月球影響球形成的圓上的搜索算法。另外本文的描述方式把地月轉(zhuǎn)移軌道求解問題解耦成平面內(nèi)轉(zhuǎn)移軌道形狀計算問題以及轉(zhuǎn)移軌道空間定向問題，從而大大降低了問題求解維度。

1.1 轉(zhuǎn)移三角形形狀計算問題

對于從天梯出發(fā)的奔月任務(wù)來說，為了擴大任務(wù)窗口，可以適當允許出發(fā)轉(zhuǎn)移軌道有小傾角，以盡可能節(jié)約推進劑，比如要求i∈[-0.5°,0.5°]。由于天梯質(zhì)心位于地球定點位置，相應(yīng)的升交點赤經(jīng)為：

Ω=Ω0+ωEt

其中，Ω0為初始點降交點赤經(jīng)，ωE為地球自轉(zhuǎn)角速度，t為據(jù)初始點時刻。

如圖1所示，定義地心為E，月心為L，軌道面與月球影響球切割形成圓的圓心為O。轉(zhuǎn)移軌道在地球近地點出發(fā)，該點定義為P，到達月球影響球的入口點A。

圖1 軌道面與月球影響球的切割示意圖Fig.1 Demonstration of section of orbital plane and the sphere of influence of moon

轉(zhuǎn)移軌道面與月球影響球切割的交線必然形成一個圓，該圓上各點到月心的距離等于影響球半徑rs，rs=66 200 km。轉(zhuǎn)移軌道面法方向或軌道動量矩方向nOL為：

nOL=[sinisinΩ,-sinicosΩ,cosi]T

(1)

cos∠OLE|nOL·rL|/rL

(2)

其中，rL為地月距離。

由圖(1)所示的幾何關(guān)系可知，順行軌道(i,Ω)和逆行軌道(π-i,Ω+π)都滿足以上幾何關(guān)系(二者軌道面法向相差180°)，相應(yīng)的式(2)取絕對值，可根據(jù)工程實際約束來選取。

l=rLcos∠OLE

(3)

地月轉(zhuǎn)移軌道面與月球影響球必須相交，否則轉(zhuǎn)移軌道與月球影響球沒有交點，就無法形成環(huán)月軌道。對應(yīng)的數(shù)學約束條件就是：

l≤rs

(4)

發(fā)射窗口搜索時間區(qū)間出現(xiàn)在滿足式(4)的時間段，該約束條件可以極大地節(jié)約轉(zhuǎn)移窗口的計算量。

軌道面與影響球切割而成的圓半徑OA的大小r0為：

(5)

(6)

k-r0≤rEA≤k+r0

(7)

圖2 軌道面與月球影響球的截面平面示意圖Fig.2 Section plane of orbital plane and the sphere of influence of moon

由地月幾何關(guān)系，rEA取值范圍可擴大為：

rL-rs≤k-r0≤rEA≤k+r0≤rL+rs

(8)

已知地月平均距離為384 400 km，因此地月距離rL考慮一定余量就可以得到rEA的取值范圍，如取余量Δ=0.1可得：

(rL-rs)(1-Δ)≤rEA≤(rL+rs)(1+Δ)

即

286 380 km≤rEA≤495 660 km

(9)

在給定時間tPA內(nèi)，從近地出發(fā)點P運行到入口點A的地月轉(zhuǎn)移軌道構(gòu)成了一個Lambert問題。根據(jù)Lambert定理，地月轉(zhuǎn)移出發(fā)點P到入口點A的飛行時間tPA僅與軌道半長軸、兩點與地心距離之和rPE+rEA以及弦長rPA有關(guān)。又由余弦定理，給定rPE和rEA，則弦長rPA和真近點角θPA一一對應(yīng)。因此地月轉(zhuǎn)移軌道問題轉(zhuǎn)化為已知tPA求θPA的問題。E、P、A確定了一個與轉(zhuǎn)移軌道相關(guān)的三角形，稱為轉(zhuǎn)移三角形。

賀波勇推導了tPA和θPA的導數(shù)關(guān)系。將tPA作為設(shè)計變量，用牛頓迭代法求解θPA[9]。本文的算例測試發(fā)現(xiàn)，該方法對初值較敏感，收斂性一般，計算量較大，且計算精度也不太高，為此改用Brent法求解。

Brent法在20世紀60年代，由阿姆斯特丹數(shù)學中心的Van Wijngaarden和Dekker等研究成功，并在1973年由Brent進行了改進。該方法特別適合于求解一維非線性方程的根，而無須提供函數(shù)的導數(shù)。該方法具有超線性收斂特性，又能夠保證二分法的收斂確定性。具體算法參見文獻[15]。

按照Brent算法，選擇tPA作為設(shè)計變量，求解θPA，二者的求解關(guān)系如下：

地月轉(zhuǎn)移軌道由近地點出發(fā)，存在如下關(guān)系：

rEA=p/(1+ecosθPA)

p=rEP(1+e)

可得到：

(10)

(11)

其中：p為半通徑；e為偏心率；a為半長軸。

至此由θPA、a、e不難求得轉(zhuǎn)移軌道飛行時間tPA，可以將求得的tPA與給定的轉(zhuǎn)移軌道飛行時間tPA0之差作為待求的一維非線性目標函數(shù)，對應(yīng)的未知量θPA可由Brent法求解。

轉(zhuǎn)移三角形形狀求解問題在整個搜索過程中重復(fù)出現(xiàn)，但由本文模型的描述方式可知，該三角形與軌道傾角、升交點赤經(jīng)和近地點幅角無關(guān)，因此該求解問題與整體窗口搜索算法解耦，可以把轉(zhuǎn)移三角形所有可能的軌道組合提前求解出來，在接下來的窗口搜索中就只需要匹配這些已知的轉(zhuǎn)移三角形，從而大大降低問題的搜索維度，節(jié)省計算量。此外轉(zhuǎn)移三角形形狀求解問題也從原來的三維影響球的球面經(jīng)緯度搜索問題降維成為平面內(nèi)圓上點的搜索問題，進一步節(jié)省了計算量。

1.2 轉(zhuǎn)移三角形的空間定向問題

1)對于如近地軌道空間站出發(fā)的地月轉(zhuǎn)移軌道，i和初始Ω0一般已知，則當時的Ω可由到達入口點時刻和飛行時間來決定，未知量為近地點幅角ω；

2)對于天梯奔月任務(wù)來說，天梯位于赤道上方，ω等于0°(順行軌道)或者180°(逆行軌道)，未知量為i。

以上兩種情況均可由Brent算法求解。

由地心慣性系到軌道系的坐標旋轉(zhuǎn)矩陣為：

其中：Lx(?)表示繞x軸旋轉(zhuǎn)?角度的坐標旋轉(zhuǎn)矩陣；Lz(?)表示繞z軸旋轉(zhuǎn)?角度的坐標旋轉(zhuǎn)矩陣；u=ω+θ為緯度幅角，其中θ為真近點角。

地心慣性系下的入口點位置矢量rEA為：

(12)

由圖1可知，rEA與月球位置矢量rL(由星歷計算)及入口點到月心的矢量rLA滿足如下幾何關(guān)系：

rLA=rEA-rL

而rLA的大小等于月球影響球半徑rs，則：

|rEA-rL|-rs=0

(13)

因此根據(jù)未知的ω或者i，由Brent算法預(yù)先判斷所求區(qū)間是否有解，再求解式(13)對應(yīng)的一維非線性方程的根，就可以得到滿足條件的未知的ω或者i[15]。即：

1)對于i給定的轉(zhuǎn)移軌道，ω∈[0°,360°]，由式(12)和式(13)，根據(jù)已知的五個軌道根數(shù)，由Brent算法求解ω；

2)對于天梯任務(wù)，i∈[-0.5°,0°]和i∈[0°,0.5°]對應(yīng)的逆行和順行軌道，同樣由式(12)和式(13)，根據(jù)已知的其他五個軌道根數(shù)，由Brent算法求解i。

然后根據(jù)式(1)～(6)計算轉(zhuǎn)移三角形數(shù)據(jù)庫中待評估的轉(zhuǎn)移三角形是否滿足式(7)。

以上求解軌道面參數(shù)的過程稱為轉(zhuǎn)移三角形的空間定向問題。

以上求解出了地月轉(zhuǎn)移段軌道六個軌道根數(shù)，還需要進一步判斷該軌道進入月球影響球被月球捕獲以后形成的環(huán)月軌道是否滿足給定條件。此時需把地心慣性系下的入口點處的矢量rLA和vLA轉(zhuǎn)換到月心白道系下[2]：

M=Lz(uL+π)Lx(iL)Lz(ΩL)

其中：M為地心慣性系到月心瞬時白道系的旋轉(zhuǎn)矩陣。上式中的上標L表示月球相關(guān)的參數(shù)。

基于上式求得的月心系下的位置和速度，就可以求得月心段繞月軌道參數(shù)，之后可以判斷是否滿足需要的環(huán)月軌道條件，比如環(huán)月軌道高度約束等。

以上是普通地月轉(zhuǎn)移軌道的求解流程。對于載人航天任務(wù)常用的自由返回軌道，可基于對稱性原理，進一步求得月地轉(zhuǎn)移段軌道參數(shù)，再判斷是否滿足載人航天任務(wù)常用的自由返回軌道約束條件，如返回地球時的近地點高度、再入角等。本文只計算一般的地月轉(zhuǎn)移軌道，自由返回軌道求解的具體流程可參考Battin[1]、黃文德[4-5]、賀波勇[8]等的報道。

2 天梯奔月解法

2.1 CPU并行計算技術(shù)簡介

當前CPU主頻的提升已經(jīng)明顯遇到了瓶頸，計算機硬件界主要關(guān)注如何發(fā)展多核CPU。目前民用CPU已經(jīng)擁有十多個核心。針對CPU核心數(shù)越來越多的現(xiàn)狀，科學計算必須考慮并行/并發(fā)程序設(shè)計，避免浪費計算資源。

并行計算一般可分為計算密集型、數(shù)據(jù)密集型、網(wǎng)絡(luò)密集型[10]，本文主要關(guān)注計算密集型計算。1994年誕生的MPI已經(jīng)成為集群并行程序設(shè)計事實上的標準。1997年誕生的OpenMP已經(jīng)成為單機多核并行事實上的標準[10]。本文的并行程序主要針對單機多核并行架構(gòu)，而單機多核圖形處理器(graphics processing unit, GPU)加速并行架構(gòu)和多機多核并行架構(gòu)不在本文考慮。常用的并行分解技術(shù)有兩種：遞歸分解和數(shù)據(jù)分解[11-12]。本文著重關(guān)注單機CPU并行程序設(shè)計，并主要采用數(shù)據(jù)分解的方式來設(shè)計并行程序。

本文綜合考慮到可移植性、編程方便性、程序穩(wěn)健性等，避免采用復(fù)雜的MPI，而是采用Java語言，并行計算程序設(shè)計基于Parallel Java 2 Library(PJ2庫)[13]。

PJ2庫是RIT(Rochester Institute of Technology)的Kaminsky教授開發(fā)的Java語言的并行計算調(diào)度庫，并應(yīng)用于RIT計算機學院的超級計算機集群上。其主要設(shè)計思想是簡化并行計算程序設(shè)計的難度，特別是避免使用復(fù)雜的MPI，能夠極好地適用于單機CPU、GPU并行以及分布式集群異構(gòu)并行計算。PJ2庫有五種CPU核心并行調(diào)度機制：Fixed、Leapfrog、Dynamic、Proportional、Guided。后四種調(diào)度機制代表了當前常見的先進并行調(diào)度機制，具體介紹參見文獻[13]。本文選擇Dynamic調(diào)度機制作為示例。

2.2 天梯地月轉(zhuǎn)移軌道窗口計算并行設(shè)計

天梯系統(tǒng)的基本組成如圖3所示，主要包括天梯繩索、地球表面節(jié)點、GEO節(jié)點、頂點錨、繩索攀爬器、地球表面運輸系統(tǒng)和運營中心[14]。

圖3 天梯系統(tǒng)組成示意圖[15]Fig.3 Composition of space elevator system[15]

天梯繩索是天梯最重要的組成部分，其主要功能為向攀爬器提供依附的途徑，使攀爬器能夠沿著纜繩從地表節(jié)點進入空間，同時其還要平衡整個天梯的重力與系統(tǒng)繞地球轉(zhuǎn)動形成的離心力；攀爬器為運輸有效載荷進入空間的運載器，其沿著繩索爬升，進入軌道預(yù)定位置后將有效載荷釋放；地球表面節(jié)點和頂點錨分別位于繩索的兩端，起到固定繩索位置、調(diào)整繩索姿態(tài)的功能。整個天梯的質(zhì)心設(shè)置在GEO節(jié)點位置，這樣整個天梯系統(tǒng)就以地球的自轉(zhuǎn)速度圍繞地球旋轉(zhuǎn)，從而保證天梯與地面相對靜止[14,16]。

由于天梯從地表一直延伸到GEO更遠端，在執(zhí)行行星際發(fā)射任務(wù)時，在天梯不同的高度釋放可以獲得不同的能量。比如，在天梯遠端頂點錨處(距離地心100 000 km)釋放甚至無須施加速度增量就可以飛往火星[16]。本文設(shè)計的從天梯飛往月球的轉(zhuǎn)移軌道算法，重點解決在天梯哪一高度釋放才能以最小的代價飛往月球。

本文假定天梯初始位于升交點赤經(jīng)60°赤道上空。轉(zhuǎn)移軌道約束地月轉(zhuǎn)移軌道真近點角差θPA<π，近月點軌道高度為50 km≤rH≤1 000 km?？紤]到2025年3月白道面與赤道面相對夾角達到最大值，約為28.72°，本文窗口計算的時間段取為2025年3月。

本文采用基于數(shù)據(jù)分解的并行方式，主要設(shè)計兩級并行算法：轉(zhuǎn)移三角形形狀并行計算問題和轉(zhuǎn)移三角形的空間定向并行計算問題。具體流程如圖4所示。

圖4 軌道設(shè)計流程Fig.4 Flow chart of trajectory design

第一個問題是轉(zhuǎn)移三角形形狀計算問題并行化：在初估轉(zhuǎn)移軌道特性時，rEA根據(jù)式(9)均勻選擇若干個點作為并行化設(shè)計參數(shù)，從近地點出發(fā)飛行到入口點的飛行時間tPA取30～170 h，基于Brent算法求解轉(zhuǎn)移三角形的真近點角，求得的結(jié)果存儲下來用于下一個子問題的搜索。

第二個問題是轉(zhuǎn)移三角形的空間定向問題并行化：天梯地月轉(zhuǎn)移軌道窗口計算時，以到達入口點時刻tA作為劃分并行顆粒度的設(shè)計變量，計算1個月的時間周期。以到達入口點時刻作為并行顆粒度劃分的主要好處是可以在并行計算之前計算好tA時刻的月球星歷。在計算空間定向時，軌道傾角i∈[-0.5°,0.5°]。在前一問題中求得的滿足約束的轉(zhuǎn)移三角形數(shù)據(jù)庫中，求解滿足式(13)的Brent子問題，從而求解出對應(yīng)的轉(zhuǎn)移軌道空間定向參數(shù)軌道傾角。然后判斷所求得的軌道是否滿足地月轉(zhuǎn)移軌道條件，如果滿足則找到一條可行的轉(zhuǎn)移軌道。

2.3 天梯奔月窗口計算結(jié)果

下面分別計算GEO點釋放和天梯遠端頂點錨處釋放奔月的計算結(jié)果，分別如圖5和圖6所示。本文的算例中，采用一臺Windows 10系統(tǒng)的電腦(AMD 3700X八核十六線程CPU，主頻3.6 GHz，內(nèi)存8G)，rEA取500個點，tPA間隔取0.2 h，tA間隔取10 min，并行計算在天梯GEO點釋放時1個月的奔月窗口需要437.57 s，串行計算需要3 040.32 s，并行加速比為6.95。

由圖5和圖6的計算結(jié)果表明，從天梯奔月每月內(nèi)存在兩次窗口，每次窗口持續(xù)4～6天。對于奔月任務(wù)，天梯頂點錨處釋放所需的速度增量反而比在GEO節(jié)點直接奔月所需的速度增量大，這主要是由于在遠端頂點錨處釋放的地月轉(zhuǎn)移軌道獲得了過多的能量，此能量甚至能飛行到火星，但不適合用于月球軌道轉(zhuǎn)移。為了詳細計算窗口，還需要把rEA、tPA和tA的間隔取得更密集一些，但這里用于了解轉(zhuǎn)移軌道的特性已經(jīng)足夠。

(a) 飛行時間結(jié)果(a) Results of flight time

(a) 飛行時間結(jié)果(a) Results of flight time

為尋找到底在何處釋放可以以較小代價飛往月球，計算了天梯不同軌道高度釋放的奔月窗口，計算結(jié)果表明：在天梯距地心51 000 km處釋放進入地月轉(zhuǎn)移軌道，所需要的速度增量較小，所需速度增量如圖7所示。在此處釋放飛行器獲得的相對地心慣性系的線速度為3.719 km/s。

(a) 飛行時間結(jié)果(a) Results of flight time

為詳細研究從地心距51 000 km處釋放的地月轉(zhuǎn)移軌道特性，rEA取2 000個點，tPA間隔取0.1 h，tA間隔取1 min，并行計算天梯在GEO點釋放時1個月的奔月窗口需要30 220.12 s(合503.7 min)。

如圖8所示，在距地心51 000 km處的天梯奔月環(huán)月軌道傾角分布在10°～170°之間。

圖8 月球軌道傾角與升交點赤經(jīng)Fig.8 Inclination and RAAN of Moon orbit

天梯地月轉(zhuǎn)移軌道的月球影響球入口點如圖9所示，在月球影響球上西半球，在南、北緯度均有分布。

圖9 入口點經(jīng)緯度分布Fig.9 Latitude and longitude of entry point

以尋找到的速度增量最小的轉(zhuǎn)移軌道為例，輸入STK進行驗證，計算得到的轉(zhuǎn)移軌道根數(shù)如表1所示，到入口點飛行時間111.9 h，入口點對應(yīng)歷書時儒略日(對應(yīng)STK的JED時間)為2 460 747.973 611 11，求得的最小速度增量為0.855 m/s，近月點高度為517.040 km。

表1 地月轉(zhuǎn)移軌道入口點根數(shù)

基于商業(yè)軟件STK的Astrogator模塊進行驗證，STK/Astrogator在地球段和月球段軌道預(yù)報器分別選地球點質(zhì)量和月球點質(zhì)量模型，STK計算得到的近月點高度為517.645 km。本文計算結(jié)果與之接近，顯示了本文算法的正確性。設(shè)計結(jié)果如圖10所示。

圖10 天梯地月轉(zhuǎn)移軌道Fig.9 Earth-Moon transfer orbit from space elevator

3 結(jié)論

本文提出一種基于飛行軌道面參數(shù)描述地月轉(zhuǎn)移軌道的雙二體模型表達方式，避免采用傳統(tǒng)的基于月球影響球入口點經(jīng)緯度的描述方式，該表達方式由軌道定義出發(fā)，通過地月幾何關(guān)系來表述，更加明確直觀。通過表述方式的轉(zhuǎn)換把原來的三維球面搜索算法降維成為二維平面內(nèi)圓上的搜索算法，結(jié)合一維非線性方程求根算法Brent算法和Lambert原理，可以高效地求解地月轉(zhuǎn)移軌道的形狀參數(shù)。

通過合理的解耦分解，將原始的軌道窗口多變量搜索問題進一步分解為兩個子問題：一個為轉(zhuǎn)移三角形形狀計算問題，通過轉(zhuǎn)移軌道形成的三角形幾何關(guān)系，結(jié)合Lambert原理求解出轉(zhuǎn)移軌道的形狀參數(shù)(半長軸、偏心率及真近點角)；另一個為轉(zhuǎn)移三角形的空間定向問題，基于已求得的轉(zhuǎn)移軌道形狀數(shù)據(jù)匹配搜索可行的轉(zhuǎn)移軌道，確定剩余的三個軌道面相關(guān)參數(shù)。最后為了加速窗口搜索過程，充分發(fā)揮多核計算機算力，形成一種轉(zhuǎn)移三角形形狀計算問題和空間定向問題的兩級并行算法。

基于本文的并行圓錐曲線幾何切面法，計算了天梯不同高度進行奔月任務(wù)的軌道窗口，結(jié)果表明：

1)在天梯遠端釋放執(zhí)行奔月任務(wù)效果甚至不如在GEO點直接釋放的效果。

2)在距地心51 000 km左右釋放執(zhí)行奔月任務(wù)所需的速度增量極小。

3)天梯出發(fā)的奔月環(huán)月軌道軌道傾角分布在10°～170°之間。入口點在月球影響球上西半球，在南、北緯度均有分布。

4)計算天梯奔月1個月的發(fā)射窗口，本文的兩級并行算法加速比達6.95。