• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

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

    2022-12-02 04:45:08童科偉程炳琳汪小衛(wèi)
    國防科技大學學報 2022年6期
    關(guā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。

    猜你喜歡
    天梯入口月球
    到月球上“飆車”
    軍事文摘(2022年12期)2022-07-13 03:11:44
    陪我去月球
    飛碟探索(2022年1期)2022-03-12 14:11:20
    百丈天梯上層云
    南風(2021年35期)2021-12-03 11:55:28
    月球上的另一個我
    軍事文摘(2021年18期)2021-12-02 01:28:54
    基于新一代稱重設(shè)備的入口治超勸返系統(tǒng)分析
    玩“火”男人蔡國強的還鄉(xiāng)“天梯”
    尼阿多天梯
    文苑(2020年10期)2020-11-07 03:15:26
    秘密入口
    作品三
    第九道 靈化閣入口保衛(wèi)戰(zhàn)
    国产黄片美女视频| 我的女老师完整版在线观看| 国产亚洲最大av| 人人澡人人妻人| 一本久久精品| 亚洲精品456在线播放app| av专区在线播放| 黑丝袜美女国产一区| 日本爱情动作片www.在线观看| 精品国产一区二区久久| 成人18禁高潮啪啪吃奶动态图 | 欧美变态另类bdsm刘玥| 日韩一区二区视频免费看| 麻豆成人午夜福利视频| 精品一区二区三区视频在线| 国产精品国产三级国产专区5o| 涩涩av久久男人的天堂| 亚洲精品一二三| 最近中文字幕2019免费版| 日本黄大片高清| 国产精品福利在线免费观看| 菩萨蛮人人尽说江南好唐韦庄| 日韩制服骚丝袜av| 夜夜看夜夜爽夜夜摸| av福利片在线观看| 免费av不卡在线播放| 人妻制服诱惑在线中文字幕| kizo精华| 久热这里只有精品99| 国产精品不卡视频一区二区| 不卡视频在线观看欧美| 国产白丝娇喘喷水9色精品| 高清黄色对白视频在线免费看 | 丝袜喷水一区| 夜夜看夜夜爽夜夜摸| 亚洲怡红院男人天堂| 国产日韩一区二区三区精品不卡 | 韩国av在线不卡| 国产亚洲av片在线观看秒播厂| 美女视频免费永久观看网站| 久久久久久久久久人人人人人人| 在线观看一区二区三区激情| 少妇 在线观看| 国产精品欧美亚洲77777| 黄片无遮挡物在线观看| 我的老师免费观看完整版| 欧美成人午夜免费资源| 亚洲欧美日韩卡通动漫| 交换朋友夫妻互换小说| 亚洲国产毛片av蜜桃av| 王馨瑶露胸无遮挡在线观看| www.色视频.com| 久久久亚洲精品成人影院| 亚洲国产成人一精品久久久| 国产精品一区二区在线不卡| 国产午夜精品久久久久久一区二区三区| 99国产精品免费福利视频| 久久久久久久久久成人| h日本视频在线播放| 成年美女黄网站色视频大全免费 | 亚洲精品,欧美精品| 男女边吃奶边做爰视频| 99热全是精品| 国产成人一区二区在线| 日韩av在线免费看完整版不卡| 中文精品一卡2卡3卡4更新| 一级毛片我不卡| 黑人巨大精品欧美一区二区蜜桃 | 黑人高潮一二区| 亚洲,一卡二卡三卡| 久久精品久久久久久久性| 老司机影院毛片| 亚洲熟女精品中文字幕| 日日爽夜夜爽网站| 内地一区二区视频在线| 国产av码专区亚洲av| 少妇丰满av| 久久精品久久精品一区二区三区| 国产欧美日韩一区二区三区在线 | 2022亚洲国产成人精品| 久久精品久久精品一区二区三区| 成年av动漫网址| 亚洲中文av在线| 最黄视频免费看| 亚洲精品久久午夜乱码| 国产深夜福利视频在线观看| 插逼视频在线观看| 国产亚洲欧美精品永久| 麻豆成人午夜福利视频| 中国美白少妇内射xxxbb| freevideosex欧美| 精品少妇久久久久久888优播| 免费看av在线观看网站| 成人免费观看视频高清| 国产欧美日韩综合在线一区二区 | 久久毛片免费看一区二区三区| 国产女主播在线喷水免费视频网站| 黑人高潮一二区| 热99国产精品久久久久久7| 欧美精品高潮呻吟av久久| 天堂中文最新版在线下载| 最近最新中文字幕免费大全7| 美女福利国产在线| 丝袜喷水一区| 亚洲精品日韩在线中文字幕| 91在线精品国自产拍蜜月| 欧美日韩亚洲高清精品| 国产有黄有色有爽视频| 一级a做视频免费观看| 简卡轻食公司| 99国产精品免费福利视频| av免费观看日本| 国产黄片视频在线免费观看| 又爽又黄a免费视频| 欧美日韩综合久久久久久| 亚洲国产精品专区欧美| 老司机影院成人| 日本欧美视频一区| 蜜臀久久99精品久久宅男| 99九九线精品视频在线观看视频| 午夜精品国产一区二区电影| 日本免费在线观看一区| av一本久久久久| 国产中年淑女户外野战色| 亚洲无线观看免费| 精品久久久久久电影网| 春色校园在线视频观看| 国产成人午夜福利电影在线观看| 国产熟女午夜一区二区三区 | 纯流量卡能插随身wifi吗| 九草在线视频观看| 各种免费的搞黄视频| 大香蕉97超碰在线| 免费看不卡的av| 国产精品国产三级专区第一集| 成人18禁高潮啪啪吃奶动态图 | 精品亚洲乱码少妇综合久久| 交换朋友夫妻互换小说| 日本爱情动作片www.在线观看| 亚洲人成网站在线播| 国产精品一区二区在线观看99| 亚洲怡红院男人天堂| 最近手机中文字幕大全| 久久久欧美国产精品| 久久99精品国语久久久| 最后的刺客免费高清国语| 国产精品麻豆人妻色哟哟久久| 亚洲国产精品一区二区三区在线| 99九九线精品视频在线观看视频| 青春草国产在线视频| 亚洲国产成人一精品久久久| 久久精品国产亚洲网站| 久久热精品热| 七月丁香在线播放| 日韩成人av中文字幕在线观看| 免费大片18禁| 久久久久网色| 亚洲精品456在线播放app| 国产色爽女视频免费观看| 免费看光身美女| 免费观看a级毛片全部| 我要看日韩黄色一级片| 卡戴珊不雅视频在线播放| 狂野欧美白嫩少妇大欣赏| 伊人久久国产一区二区| 91成人精品电影| 久久精品国产鲁丝片午夜精品| 成人午夜精彩视频在线观看| 人体艺术视频欧美日本| 精品国产乱码久久久久久小说| 亚洲精品国产成人久久av| 亚洲va在线va天堂va国产| 国产无遮挡羞羞视频在线观看| 久久狼人影院| 夫妻性生交免费视频一级片| 久久99蜜桃精品久久| 高清av免费在线| 欧美日韩国产mv在线观看视频| 国产精品久久久久久精品古装| 另类亚洲欧美激情| 极品人妻少妇av视频| 欧美 亚洲 国产 日韩一| 亚洲丝袜综合中文字幕| 久久精品国产自在天天线| 精品久久久精品久久久| 尾随美女入室| 一本—道久久a久久精品蜜桃钙片| 最近2019中文字幕mv第一页| 久久久久久久久久人人人人人人| 久久影院123| 特大巨黑吊av在线直播| 亚洲国产精品一区二区三区在线| 亚洲色图综合在线观看| 一级毛片aaaaaa免费看小| 久久99精品国语久久久| av不卡在线播放| 国产精品麻豆人妻色哟哟久久| 久久精品国产亚洲av涩爱| 日韩人妻高清精品专区| 大码成人一级视频| 国产免费一区二区三区四区乱码| 国产69精品久久久久777片| 99久久精品热视频| .国产精品久久| 亚洲精品,欧美精品| 久久精品久久精品一区二区三区| 水蜜桃什么品种好| 国产欧美日韩精品一区二区| 色94色欧美一区二区| 两个人的视频大全免费| 精品午夜福利在线看| 亚洲va在线va天堂va国产| 中文乱码字字幕精品一区二区三区| 欧美激情国产日韩精品一区| 伊人久久精品亚洲午夜| 男女免费视频国产| 99久久精品热视频| 美女福利国产在线| 亚洲成人一二三区av| 51国产日韩欧美| 一区二区av电影网| 在线观看免费高清a一片| 色视频在线一区二区三区| 中文精品一卡2卡3卡4更新| 久久国产乱子免费精品| 久久久亚洲精品成人影院| 伊人亚洲综合成人网| 日韩人妻高清精品专区| 男女边摸边吃奶| 国产精品秋霞免费鲁丝片| 亚洲怡红院男人天堂| 亚洲国产av新网站| 亚洲婷婷狠狠爱综合网| 十八禁高潮呻吟视频 | 99热国产这里只有精品6| 亚洲av电影在线观看一区二区三区| 男人爽女人下面视频在线观看| 中文字幕免费在线视频6| 高清视频免费观看一区二区| 久久99蜜桃精品久久| a级毛片免费高清观看在线播放| 大香蕉97超碰在线| 一级毛片aaaaaa免费看小| 久久人妻熟女aⅴ| 久久久久精品久久久久真实原创| 一区二区三区四区激情视频| 91精品国产九色| 中文资源天堂在线| 亚洲熟女精品中文字幕| 国精品久久久久久国模美| 老女人水多毛片| 亚洲精品,欧美精品| 午夜福利影视在线免费观看| 不卡视频在线观看欧美| 亚洲精品国产成人久久av| 亚洲av在线观看美女高潮| 在线观看av片永久免费下载| 亚洲,一卡二卡三卡| 男人添女人高潮全过程视频| 极品人妻少妇av视频| 99热全是精品| 欧美xxⅹ黑人| 黄色毛片三级朝国网站 | 美女xxoo啪啪120秒动态图| 高清在线视频一区二区三区| 少妇人妻一区二区三区视频| 国产精品99久久99久久久不卡 | 三级国产精品欧美在线观看| 王馨瑶露胸无遮挡在线观看| 色5月婷婷丁香| 十分钟在线观看高清视频www | 制服丝袜香蕉在线| 18禁在线无遮挡免费观看视频| 久久毛片免费看一区二区三区| 国产成人91sexporn| 丰满饥渴人妻一区二区三| 你懂的网址亚洲精品在线观看| 丝袜脚勾引网站| 久久精品国产亚洲av天美| 日本猛色少妇xxxxx猛交久久| 亚洲怡红院男人天堂| 啦啦啦在线观看免费高清www| 在线免费观看不下载黄p国产| 精品一区二区三卡| 一区二区三区精品91| 久久ye,这里只有精品| 香蕉精品网在线| 十八禁网站网址无遮挡 | 精品国产露脸久久av麻豆| 精品午夜福利在线看| 中文乱码字字幕精品一区二区三区| 国产淫语在线视频| 日本与韩国留学比较| 精品少妇久久久久久888优播| 国产精品一区www在线观看| 一本—道久久a久久精品蜜桃钙片| 插逼视频在线观看| 亚洲一级一片aⅴ在线观看| 国产精品嫩草影院av在线观看| 国产淫语在线视频| 少妇被粗大猛烈的视频| 午夜老司机福利剧场| 极品少妇高潮喷水抽搐| 少妇高潮的动态图| 美女主播在线视频| 我要看日韩黄色一级片| 黄色视频在线播放观看不卡| 两个人免费观看高清视频 | 日韩一区二区三区影片| 在线播放无遮挡| 国产真实伦视频高清在线观看| 精品少妇久久久久久888优播| 欧美日韩亚洲高清精品| 国产精品一区二区在线观看99| 桃花免费在线播放| 最后的刺客免费高清国语| 国产成人91sexporn| 天美传媒精品一区二区| 韩国av在线不卡| 亚洲国产毛片av蜜桃av| 国产白丝娇喘喷水9色精品| 国产伦精品一区二区三区视频9| 国产高清有码在线观看视频| 老熟女久久久| 亚洲电影在线观看av| √禁漫天堂资源中文www| 国产视频内射| 九草在线视频观看| 最近最新中文字幕免费大全7| 日韩亚洲欧美综合| 久久这里有精品视频免费| 亚洲精品,欧美精品| 国产成人免费无遮挡视频| 97在线人人人人妻| 午夜福利在线观看免费完整高清在| 午夜久久久在线观看| 欧美丝袜亚洲另类| av福利片在线观看| 精品少妇内射三级| 国产白丝娇喘喷水9色精品| 免费大片18禁| 人人妻人人爽人人添夜夜欢视频 | 国产精品久久久久久精品电影小说| 久久国产精品大桥未久av | 久久久久久久久久人人人人人人| 国产精品欧美亚洲77777| 日本黄色片子视频| 欧美区成人在线视频| 少妇的逼好多水| 国产69精品久久久久777片| 国产黄片美女视频| 国产亚洲午夜精品一区二区久久| 18禁动态无遮挡网站| 欧美最新免费一区二区三区| 伦理电影免费视频| 亚洲av不卡在线观看| 美女脱内裤让男人舔精品视频| 亚洲内射少妇av| 妹子高潮喷水视频| 日本免费在线观看一区| 天堂8中文在线网| av在线app专区| 成人毛片a级毛片在线播放| 在线观看美女被高潮喷水网站| 亚洲国产精品国产精品| 最新的欧美精品一区二区| 在线观看免费高清a一片| 国产成人午夜福利电影在线观看| 国产高清国产精品国产三级| 中国美白少妇内射xxxbb| 国产亚洲最大av| 国产免费又黄又爽又色| 国产成人精品婷婷| 久久这里有精品视频免费| 少妇的逼好多水| 国产色婷婷99| 亚洲欧洲日产国产| 少妇的逼水好多| 亚洲欧洲日产国产| 18+在线观看网站| 亚洲精品国产色婷婷电影| 国产av一区二区精品久久| 久久人人爽人人爽人人片va| 久久精品熟女亚洲av麻豆精品| 黑人高潮一二区| 少妇精品久久久久久久| 日韩中文字幕视频在线看片| 日韩电影二区| 亚洲精品国产成人久久av| 另类亚洲欧美激情| 美女xxoo啪啪120秒动态图| 熟女av电影| 成人漫画全彩无遮挡| 欧美高清成人免费视频www| 久久久精品94久久精品| 亚洲精品一二三| 精品一区在线观看国产| 亚洲欧美成人精品一区二区| 美女xxoo啪啪120秒动态图| 亚洲无线观看免费| 黄色怎么调成土黄色| 亚洲在久久综合| 久久午夜综合久久蜜桃| 老女人水多毛片| 黄色配什么色好看| 99国产精品免费福利视频| 欧美精品一区二区免费开放| 高清欧美精品videossex| 精品人妻熟女毛片av久久网站| 午夜久久久在线观看| 中文字幕免费在线视频6| 亚洲一区二区三区欧美精品| 18禁在线无遮挡免费观看视频| 久久毛片免费看一区二区三区| 免费不卡的大黄色大毛片视频在线观看| 99热这里只有精品一区| 在线观看免费日韩欧美大片 | 女人精品久久久久毛片| 一级黄片播放器| 五月开心婷婷网| 亚洲av福利一区| 国产精品一二三区在线看| 中文字幕人妻丝袜制服| 亚洲自偷自拍三级| 免费看日本二区| 久久国产精品大桥未久av | 久久久久久久亚洲中文字幕| 亚洲天堂av无毛| 免费观看的影片在线观看| 精品人妻熟女毛片av久久网站| 午夜激情久久久久久久| 精品99又大又爽又粗少妇毛片| 日韩av不卡免费在线播放| 亚洲国产最新在线播放| 国产毛片在线视频| 久久人人爽人人爽人人片va| av有码第一页| 亚洲精品日韩在线中文字幕| 亚洲欧美一区二区三区国产| 久久99一区二区三区| 午夜老司机福利剧场| 国产黄片美女视频| 免费观看无遮挡的男女| 国产午夜精品一二区理论片| 中文精品一卡2卡3卡4更新| 成人美女网站在线观看视频| 亚洲怡红院男人天堂| 国产一级毛片在线| .国产精品久久| 嫩草影院新地址| 欧美高清成人免费视频www| 一边亲一边摸免费视频| 欧美精品高潮呻吟av久久| 这个男人来自地球电影免费观看 | av在线观看视频网站免费| 国产老妇伦熟女老妇高清| 久久国产精品大桥未久av | 久久精品国产亚洲av涩爱| √禁漫天堂资源中文www| 久久久国产欧美日韩av| 最黄视频免费看| 国产视频内射| a级毛片在线看网站| 国语对白做爰xxxⅹ性视频网站| 国产伦理片在线播放av一区| 国产一区亚洲一区在线观看| av线在线观看网站| 最黄视频免费看| 天堂俺去俺来也www色官网| 91午夜精品亚洲一区二区三区| 亚洲va在线va天堂va国产| 亚洲国产精品一区二区三区在线| 免费在线观看成人毛片| 免费播放大片免费观看视频在线观看| 国产色爽女视频免费观看| 一级a做视频免费观看| 中文资源天堂在线| 99久久综合免费| 亚洲av日韩在线播放| 欧美xxxx性猛交bbbb| 久久国产乱子免费精品| 亚洲国产精品国产精品| 精品久久国产蜜桃| 免费看光身美女| 人妻人人澡人人爽人人| 久久久亚洲精品成人影院| 婷婷色av中文字幕| 久久久久久久久大av| 在线观看免费视频网站a站| 久久人人爽av亚洲精品天堂| 夫妻性生交免费视频一级片| 亚洲国产精品成人久久小说| 亚洲熟女精品中文字幕| 国产亚洲av片在线观看秒播厂| 色哟哟·www| 欧美人与善性xxx| 男男h啪啪无遮挡| 亚洲电影在线观看av| 免费黄色在线免费观看| 又大又黄又爽视频免费| 男人添女人高潮全过程视频| 天堂8中文在线网| 黑丝袜美女国产一区| 欧美日韩av久久| 国产真实伦视频高清在线观看| 久久久久久久久久久丰满| 国产国拍精品亚洲av在线观看| 日韩在线高清观看一区二区三区| 热re99久久国产66热| h日本视频在线播放| 91精品国产国语对白视频| 大香蕉久久网| 亚洲欧美清纯卡通| 国产亚洲5aaaaa淫片| 黑人巨大精品欧美一区二区蜜桃 | 成人亚洲欧美一区二区av| 日韩欧美精品免费久久| 日韩av在线免费看完整版不卡| 99热6这里只有精品| 亚洲激情五月婷婷啪啪| 一级黄片播放器| 搡女人真爽免费视频火全软件| 大又大粗又爽又黄少妇毛片口| 18禁在线播放成人免费| 国产精品一区二区在线观看99| 下体分泌物呈黄色| 最黄视频免费看| 成人美女网站在线观看视频| 青青草视频在线视频观看| 人体艺术视频欧美日本| 99热6这里只有精品| 日韩不卡一区二区三区视频在线| 一级毛片aaaaaa免费看小| 久久99热这里只频精品6学生| 夜夜看夜夜爽夜夜摸| 一级黄片播放器| 国产精品秋霞免费鲁丝片| 国产欧美日韩精品一区二区| 中文字幕av电影在线播放| 在线看a的网站| 国产精品久久久久久精品古装| 中国三级夫妇交换| 亚洲无线观看免费| 插阴视频在线观看视频| 亚洲真实伦在线观看| 欧美成人午夜免费资源| 国国产精品蜜臀av免费| 国产黄片美女视频| 欧美xxxx性猛交bbbb| 少妇人妻 视频| 美女中出高潮动态图| 久久99热这里只频精品6学生| 精品熟女少妇av免费看| 啦啦啦啦在线视频资源| 97超视频在线观看视频| 麻豆成人av视频| 高清欧美精品videossex| 日本与韩国留学比较| 日本91视频免费播放| 美女xxoo啪啪120秒动态图| 插逼视频在线观看| 在线观看www视频免费| 色吧在线观看| 欧美人与善性xxx| 99久国产av精品国产电影| 国产精品成人在线| 亚洲国产毛片av蜜桃av| 日韩成人av中文字幕在线观看| 永久网站在线| 秋霞在线观看毛片| 国产 一区精品| 久久久久久久国产电影| 18禁裸乳无遮挡动漫免费视频| 人人妻人人爽人人添夜夜欢视频 | 国产欧美亚洲国产| 国产老妇伦熟女老妇高清| 日韩中字成人| 亚洲丝袜综合中文字幕| 在线观看免费高清a一片| 日韩av不卡免费在线播放| 久久精品国产a三级三级三级| 国产亚洲最大av| 美女脱内裤让男人舔精品视频| 天堂俺去俺来也www色官网| 各种免费的搞黄视频| 精品人妻熟女毛片av久久网站| 亚洲精品日本国产第一区| 26uuu在线亚洲综合色| 久久这里有精品视频免费| 18禁在线无遮挡免费观看视频| 免费高清在线观看视频在线观看| 丝袜喷水一区| 久久青草综合色| 欧美精品人与动牲交sv欧美| 亚洲图色成人| 99九九线精品视频在线观看视频| 免费观看av网站的网址| 亚洲精品456在线播放app| 亚洲内射少妇av| .国产精品久久| 国产精品一二三区在线看| 亚洲自偷自拍三级| 少妇高潮的动态图| 秋霞在线观看毛片| 久久久久视频综合| 不卡视频在线观看欧美| 人妻少妇偷人精品九色| 极品人妻少妇av视频| 国产毛片在线视频| 国产老妇伦熟女老妇高清| 亚洲欧美成人综合另类久久久| 亚洲人成网站在线观看播放| 午夜91福利影院| 国产精品无大码| 亚洲丝袜综合中文字幕|