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

    地月低能軌道轉(zhuǎn)移的混沌控制方法

    2018-08-10 02:54:08泮斌峰
    宇航學(xué)報(bào) 2018年7期
    關(guān)鍵詞:龐加萊航天器月球

    鄭 越,泮斌峰,唐 碩

    (1. 西北工業(yè)大學(xué)航天學(xué)院,西安 710072;2. 陜西省空天飛行器設(shè)計(jì)重點(diǎn)實(shí)驗(yàn)室,西安 710072)

    0 引 言

    三體模型下航天器軌道常常呈現(xiàn)混沌運(yùn)動(dòng)狀態(tài)。利用混沌的遍歷特性,能夠使航天器在不耗能的情況下飛行到想要的合理位置。但由于混沌運(yùn)動(dòng)本身具有內(nèi)在隨機(jī)性和長(zhǎng)期趨勢(shì)不可預(yù)見性,導(dǎo)致航天器在混沌區(qū)域的滑行時(shí)間過長(zhǎng)。因此,在利用三體模型下混沌運(yùn)動(dòng)低能量特性的同時(shí),如何有效減少軌道轉(zhuǎn)移時(shí)間是實(shí)現(xiàn)大規(guī)模、低成本探月活動(dòng)的迫切需求,對(duì)探月工程的發(fā)展具有非常重要的意義。

    目前采用混沌控制理論設(shè)計(jì)軌道的研究尚處于起步階段。1995年,Bollt等[1]基于龐加萊截面對(duì)流形進(jìn)行多次截?cái)鄟韺?shí)現(xiàn)地月低能轉(zhuǎn)移,有效地減少航天器在混沌區(qū)域的滑行時(shí)間,轉(zhuǎn)移過程中施加了8次脈沖。隨后,Schroer等[2]通過對(duì)一組不穩(wěn)定周期軌道進(jìn)行打靶,利用比Bollt的方法略高的能量縮短了地月轉(zhuǎn)移時(shí)間,轉(zhuǎn)移過程中需要至少4次脈沖。Macau[3]將從地球和月球附近出發(fā)的軌道進(jìn)行拼接,通過3次脈沖,進(jìn)一步用能量的增加換取了地月轉(zhuǎn)移時(shí)間的縮短。近期,Salazar等[4]發(fā)展了Bollt方法[1],通過增大龐加萊截面流形截?cái)嗉s束值,利用更大的能量和8次脈沖,以更短的時(shí)間到達(dá)月球;此外,借助于月球引力輔助變軌方式,也將混沌控制拓展至到地月三體系統(tǒng)L4,L5拉格朗日點(diǎn)的低能轉(zhuǎn)移任務(wù)。

    為了減小地月轉(zhuǎn)移過程中所需要的能量,現(xiàn)有的混沌控制方法[1-4]都利用了軌道在龐加萊截面圖上的投影點(diǎn)的穩(wěn)定方向與不穩(wěn)定方向的交點(diǎn),分別采用了同一軌道截?cái)郲1,4]和不同軌道拼接[2-3]的方法。對(duì)于同一軌道截?cái)嗟姆椒?,需要軌道本身存在位置接近的投影點(diǎn)。而每次截?cái)嘈枰谲壍擂D(zhuǎn)移過程中施加兩個(gè)脈沖,導(dǎo)致軌道轉(zhuǎn)移過程所需要的脈沖過多。即使減小龐加萊截面流形截?cái)嗉s束值以保證低能耗轉(zhuǎn)移,軌道轉(zhuǎn)移所需要的時(shí)間也較長(zhǎng)。對(duì)于不同軌道拼接的方法,通常需要利用周期軌道作為中間軌道[2]而實(shí)現(xiàn),具有依賴性。

    針對(duì)現(xiàn)有混沌控制方法存在的需要脈沖次數(shù)多和必須通過固定周期軌道等缺點(diǎn),本文提出一種新的低能地月軌道轉(zhuǎn)移方法。該方法針對(duì)地月三體系統(tǒng)中混沌運(yùn)動(dòng)的內(nèi)在隨機(jī)性,將混沌區(qū)域在龐加萊截面圖上進(jìn)行分層。進(jìn)一步分析航天器在混沌區(qū)域各層之間的運(yùn)動(dòng)趨勢(shì),找到了在盡可能減少飛行時(shí)間的前提下完成地月軌道轉(zhuǎn)移的途徑。通過粒子群算法確定地月轉(zhuǎn)移軌道通過L1瓶頸區(qū)的狀態(tài)點(diǎn),利用向前-向后法在盡可能減少飛行時(shí)間的前提下通過2~4次脈沖實(shí)現(xiàn)地月低能軌道轉(zhuǎn)移。

    1 地月三體系統(tǒng)

    1.1 動(dòng)力學(xué)模型

    考慮地月圓型限制性三體問題(Circular restricted three-body problem,CRTBP),質(zhì)量可以忽略的航天器在地球和月球兩個(gè)主天體的引力作用下的運(yùn)動(dòng),其中主天體繞其公共質(zhì)心做勻速圓周運(yùn)動(dòng)[5-6]。歸一化地球質(zhì)量m1和月球質(zhì)量m2之和為單位質(zhì)量,則:

    (1)

    (2)

    其中,Ω為旋轉(zhuǎn)坐標(biāo)系下的等效勢(shì)能:

    (3)

    設(shè)r1,r2為航天器到兩個(gè)主天體的距離,有:

    (4)

    CRTBP系統(tǒng)中,存在雅可比積分如下:

    (5)

    在地月三體系統(tǒng)中,單位長(zhǎng)度為地月之間的距離,L=3.844×105km,單位時(shí)間為T=104 h,單位速度為V=1024 m/s。

    1.2 龐加萊截面圖及分層

    圖1為地月三體系統(tǒng)中雅可比能量C=3.17948的龐加萊截面圖,E,M分別表示地球和月球。以平動(dòng)點(diǎn)L1[11]為分界線,將龐加萊截面區(qū)域分為地球引力區(qū)和月球引力區(qū)。在地球引力區(qū)中存在四組穩(wěn)定的KAM環(huán),分別用Ⅰ,Ⅱ,Ⅲ,Ⅳ表示。本文利用KAM環(huán)對(duì)龐加萊截面地區(qū)引力區(qū)中的混沌區(qū)域進(jìn)行分層。首先繪制第Ⅱ, Ⅲ, Ⅳ組KAM環(huán)的包絡(luò)線,形成的封閉區(qū)域分別為2個(gè),3個(gè)和1個(gè)。然后將第Ⅲ組KAM環(huán)靠近地球一端的包絡(luò)線相連,靠近月球一端的包絡(luò)線與第Ⅳ組包絡(luò)線靠近月球一端的邊緣連接,則可以定義第Ⅲ組KAM環(huán)和第Ⅳ組KAM環(huán)的包絡(luò)線之間的區(qū)域?yàn)镃層。最后將第Ⅱ組KAM環(huán)靠近地球一端的包絡(luò)線相連,靠近月球一端的包絡(luò)線以其邊緣的曲度延伸并將龐加萊截面地球引力區(qū)分為兩個(gè)部分,定義靠近地球部分的混沌區(qū)域?yàn)锳層,而靠近月球部分C層以外的混沌區(qū)域?yàn)锽層。

    2 混沌控制在低能地月軌道轉(zhuǎn)移中的應(yīng)用

    混沌系統(tǒng)具有內(nèi)在的敏感性,表現(xiàn)為對(duì)系統(tǒng)施加小的擾動(dòng),會(huì)使系統(tǒng)偏離固有軌跡而運(yùn)動(dòng)[12-13]。這個(gè)擾動(dòng)可以非常小,不改變系統(tǒng)動(dòng)力學(xué)特征,而驅(qū)使固有的系統(tǒng)動(dòng)力學(xué)到最終想要的狀態(tài)。這個(gè)混沌控制的過程也被稱為目標(biāo)命中[14]。由于在混沌系統(tǒng)中,微小的擾動(dòng)能夠?qū)е戮薮蟮慕K端狀態(tài)偏差,混沌控制的目的就是確定何時(shí)以及如何選擇微小的擾動(dòng)來達(dá)到預(yù)期的終端狀態(tài)。在地月三體系統(tǒng)中,混沌控制的目的就是通過微小擾動(dòng)產(chǎn)生的能耗改變航天器軌跡以盡可能縮短飛行時(shí)間。

    本文根據(jù)地月系統(tǒng)龐加萊截面圖上軌道投影點(diǎn)的分布情況,進(jìn)一步提出了地月軌道轉(zhuǎn)移時(shí)間優(yōu)化原則。利用向前-向后法,在節(jié)省時(shí)間的前提下,拼接出地月低能轉(zhuǎn)移軌道。

    2.1 地月軌道轉(zhuǎn)移時(shí)間優(yōu)化原則

    航天器在地月三體系統(tǒng)下不耗費(fèi)任何能量的軌道轉(zhuǎn)移可以認(rèn)為是在地球引力區(qū)域的滑行過程中受到月球小引力的影響,與地月之間連線的相交位置逐漸向月球靠近,最終滑行至月球引力區(qū)而到達(dá)月球附近的過程。由于在航天器的滑行過程中會(huì)經(jīng)過混沌區(qū)域的很多層,而過多的混沌層和在每一層內(nèi)過長(zhǎng)的游蕩時(shí)間,是地月軌道轉(zhuǎn)移時(shí)間過長(zhǎng)的主要原因。因此要減少地月軌道轉(zhuǎn)移時(shí)間就要讓航天器只經(jīng)過必經(jīng)的混沌層并且盡可能減少在每一層的游蕩時(shí)間。

    圖1為地月三體系統(tǒng)中雅可比能量C=3.17948的龐加萊截面圖。由地月CRTBP三體動(dòng)力學(xué)模型的公式(2)可以得出,在每個(gè)主天體的引力區(qū),離主天體距離越遠(yuǎn),航天器可以達(dá)到的最大速度越小。航天器在地月系統(tǒng)的地球引力區(qū)的滑行過程是一個(gè)在動(dòng)能作用下遠(yuǎn)離地球,而在地球引力作用下速度逐漸變小而再向地球附近滑行的一個(gè)往復(fù)過程。一旦航天器動(dòng)能可以克服地球引力,航天器將穿過L1瓶頸區(qū)域并受月球引力作用從而到達(dá)月球附近。而龐加萊截面圖的包絡(luò)線可以認(rèn)為是航天器在離主天體相應(yīng)距離下可以達(dá)到的最大速度,可以直接影響航天器滑行時(shí)的軌跡形狀。

    第1.2節(jié)中分析了地月軌道初始位置的選擇原則,即應(yīng)該從軸上第Ⅰ組KAM環(huán)的右側(cè)開始向月球方向移動(dòng)。圖2為初始位置從第Ⅰ組KAM環(huán)的右側(cè)向月球方向過渡的一組周期軌道(圓圈表示周期軌道的初始位置);圖3為地月系統(tǒng)的龐加萊截面圖。圖2中的周期軌道圖2(a),2 (b),2 (c), 2(d),2 (e), 2(f) ,2 (g),2 (h)在圖3中的投影點(diǎn)分別對(duì)應(yīng)于周期軌道a, b, c, d, e, f, g, h。圖3中,投影點(diǎn)位于KAM環(huán)圓心的周期軌道f為穩(wěn)定周期軌道,對(duì)其周圍軌跡流形具有吸引作用;而投影點(diǎn)位于混沌區(qū)域的周期軌道a, b, c, d, e, g, h為不穩(wěn)定周期軌道,可以作為三體系統(tǒng)下軌跡流形狀態(tài)隨著初始點(diǎn)變化的分析依據(jù)。將圖3中的周期軌道與圖1中混沌區(qū)域的分區(qū)圖相對(duì)應(yīng),周期軌道a, b在A層中,通過A層在龐加萊截面圖上所處的位置可知軌道流形在A層內(nèi)距離地球近且速度較大,流形遠(yuǎn)離地球再回到地球附近的一個(gè)過程中軌跡流形如軌道a, b形成交叉點(diǎn)。同時(shí),A層中的周期軌道隨著初始位置的改變,會(huì)經(jīng)過一系列如軌道b的過渡而進(jìn)入B層,演變?yōu)閳D2(c)中周期軌道的形狀(圖2(c)中周期軌道為文獻(xiàn)[2]中規(guī)定的中間目標(biāo))。在B層中,航天器遠(yuǎn)離和靠近地球的一個(gè)往復(fù)過程與地球的距離顯著增大,軌道形狀也如圖2(d), 2(e)發(fā)生了很大的變化。而在B層內(nèi)航天器已經(jīng)可以脫離地球引力的束縛而進(jìn)入月球引力區(qū),如圖2(e)所示。當(dāng)軌道在龐加萊截面圖中處于以周期軌道f為中心的KAM組的內(nèi)圈即C層時(shí),航天器遠(yuǎn)離和靠近地球的一個(gè)往復(fù)過程與地球的距離進(jìn)一步增大,且在C層內(nèi)航天器非常容易滑行至月球引力區(qū)。圖2(h)為軌道從C區(qū)向月球附近轉(zhuǎn)移的流形圖。

    通過圖2(a), 2(b), 2(c), 2(d),2(e), 2(f), 2(g),2(h)可以演變出航天器從地球附近向月球附近軌道轉(zhuǎn)移的過程。轉(zhuǎn)移軌道在混沌區(qū)域的過渡過程可以通過A→B→C或者A→B。為了節(jié)省轉(zhuǎn)移時(shí)間,軌道應(yīng)該通過A→B轉(zhuǎn)移至月球附近,所以軌道經(jīng)過C層是軌道轉(zhuǎn)移耗時(shí)過長(zhǎng)的一個(gè)重要原因。軌道轉(zhuǎn)移耗時(shí)過長(zhǎng)的另一個(gè)原因是軌道流形會(huì)在相鄰兩層之間來回滑行。如圖2(g)所示,其軌跡在龐加萊截面圖3上的投影g在B, C層之間滑行會(huì)消耗過多的飛行時(shí)間。所以,為了節(jié)省地月軌道轉(zhuǎn)移時(shí)間,就是要使航天器盡快從A層滑行至B層然后再滑行至月球附近,且避免軌道流形在兩層之間來回滑行。

    2.2 地月軌道拼接方法

    地月系統(tǒng)的混沌區(qū)域內(nèi),在初始端施加小的擾動(dòng)可以使軌跡終端的狀態(tài)發(fā)生很大的變化,而混沌本身具有遍歷性,可以在初始端速度改變很小的前提下使航天器軌跡遍布廣泛區(qū)域,通過兩條軌跡的拼接來實(shí)現(xiàn)軌道轉(zhuǎn)移。

    軌道間的拼接是一個(gè)目標(biāo)命中的過程,可以用Shinbrot等提出的向前-向后法(Forward-backward method)來實(shí)現(xiàn)[15]。算法的思想如圖4所示,為了找到從初始狀態(tài)點(diǎn)進(jìn)入以εt為半徑的目標(biāo)域的一段軌跡,對(duì)以(p0,q0)為中心長(zhǎng)度為δ的控制段(p,q)施加控制(p,q分別為位置和速度)并向前迭代,而小的目標(biāo)域向后迭代,如果存在一條同時(shí)通過初始位置附近和目標(biāo)區(qū)域的軌跡,則控制段和目標(biāo)區(qū)域的延展方向?qū)⒆罱K相交,從交點(diǎn)演變到控制段的點(diǎn)pint出發(fā)可以直接到達(dá)目標(biāo)區(qū)域。

    由于在軌道轉(zhuǎn)移過程中施加脈沖只能改變航天器的速度,不能改變其位置,本文將龐加萊截面圖上地球附近的初始點(diǎn)P為中心的區(qū)域作為目標(biāo)域,而控制段上狀態(tài)點(diǎn)位置坐標(biāo)位于L1瓶頸處,橫向速度的大小由長(zhǎng)度為δ的線段確定,且控制段上所有狀態(tài)點(diǎn)速度大小相同。目標(biāo)域不穩(wěn)定方向上的點(diǎn)正向時(shí)間積分會(huì)逐漸遠(yuǎn)離以P為初始點(diǎn)的軌跡,控制段穩(wěn)定方向上的點(diǎn)反向時(shí)間積分會(huì)逐漸遠(yuǎn)離以控制段中心點(diǎn)為初始點(diǎn)的軌跡。目標(biāo)域和控制段經(jīng)過數(shù)次迭代后,其延展方向在龐加萊截面圖上相交,通過交點(diǎn)反向時(shí)間和正向時(shí)間積分,可以得到一條從P點(diǎn)附近出發(fā)到達(dá)L1瓶頸的一段軌跡。由于盡可能縮短航天器在混沌區(qū)域的運(yùn)送時(shí)間是地月軌道轉(zhuǎn)移的重要目標(biāo),應(yīng)該盡量減少控制段和目標(biāo)域的迭代次數(shù)。當(dāng)控制段與目標(biāo)域迭代過程中位置很近,延展方向又幾乎平行不易存在交點(diǎn)時(shí),可以通過一個(gè)脈沖來拼接軌[16]。

    2.3 控制段的確定方法

    在PSO算法中,有m個(gè)n維的向量的粒子組成一個(gè)種群,每個(gè)粒子對(duì)應(yīng)一個(gè)位置向量和一個(gè)速度向量。每個(gè)粒子根據(jù)自己的飛行經(jīng)歷和種群中其他粒子的飛行經(jīng)歷來調(diào)整自己的飛行,通過跟蹤個(gè)體最優(yōu)粒子和全局最優(yōu)粒子來控制自己的運(yùn)動(dòng)。在粒子尋優(yōu)過程中,第i個(gè)粒子第s維向量的運(yùn)動(dòng)軌跡描述如下:

    (6)

    式中:c1和c2為學(xué)習(xí)因子,是非負(fù)常數(shù),r1s和r2s服從[0,1]上的均勻分布的隨機(jī)數(shù),vis是粒子的速度,pis為局部最優(yōu)位置,是第i個(gè)粒子迄今為止搜索到的最優(yōu)位置,pgs為全局最優(yōu)位置,ω是慣性權(quán)值。

    本文算法中,控制段中心點(diǎn)的候選狀態(tài)點(diǎn)向量與每個(gè)粒子的二維位置向量(y,θ)相對(duì)應(yīng),y表示狀態(tài)點(diǎn)在縱坐標(biāo)上的位置,在固定雅可比積分下,該點(diǎn)的速度可以通過式(5)得到:

    (7)

    θ表示狀態(tài)點(diǎn)速度的方向,則:

    (8)

    由第2.2節(jié)可知,軌道可以通過目標(biāo)域和控制段的延展方向在龐加萊截面圖上的交點(diǎn)完成拼接。將目標(biāo)域數(shù)次迭代后形成的區(qū)域定義為拼接區(qū)域F,本文根據(jù)軌跡與F區(qū)的最短縱坐標(biāo)距離計(jì)算粒子適應(yīng)度值,方法如圖5所示。軌跡從每個(gè)狀態(tài)點(diǎn)出發(fā),如果在預(yù)定時(shí)間TM內(nèi),正向時(shí)間積分能夠到達(dá)月球附近,在預(yù)定時(shí)間TF內(nèi),反向時(shí)間積分不經(jīng)過C層,則將反向時(shí)間積分的軌跡在龐加萊截面圖上的投影點(diǎn)與區(qū)域F最短的縱坐標(biāo)距離作為這個(gè)狀態(tài)點(diǎn)所對(duì)應(yīng)粒子的適應(yīng)度值;否則將該粒子的適應(yīng)度值設(shè)定為最大值。計(jì)算每個(gè)粒子的最小適應(yīng)度值找出局部最優(yōu)位置和全局最優(yōu)位置,最后得出的全局最優(yōu)位置即為控制段中心點(diǎn)對(duì)應(yīng)的狀態(tài)向量。

    2.4 地月軌道轉(zhuǎn)移設(shè)計(jì)方法

    通過地月三體系統(tǒng)中混沌區(qū)域分層和軌道在各層間轉(zhuǎn)移的規(guī)律能夠得到盡可能快的實(shí)現(xiàn)地月軌道轉(zhuǎn)移的途徑;向前-向后法可以盡可能地減小地月軌道轉(zhuǎn)移過程中所需要的能量。將這兩種優(yōu)點(diǎn)相結(jié)合,能夠在節(jié)省時(shí)間的前提下,有效實(shí)現(xiàn)低能地月軌道轉(zhuǎn)移。

    為了盡可能減少軌道轉(zhuǎn)移所需要的時(shí)間,一方面,目標(biāo)域正向時(shí)間積分并記錄每次迭代在龐加萊截面圖上的投影點(diǎn),當(dāng)?shù)鶱次時(shí)如果有部分投影點(diǎn)能夠離開A層,則將這些投影點(diǎn)組成的區(qū)域設(shè)定為拼接區(qū)域F;另一方面,從控制段出發(fā)的參考軌道正向時(shí)間積分可以接近月球附近的KAM 環(huán)Q,反向時(shí)間積分不經(jīng)過C層且盡可能快地到達(dá)拼接區(qū)域F。通過向前-向后法將目標(biāo)域和控制段的延展方向在區(qū)域F進(jìn)行拼接。目標(biāo)域內(nèi)地月軌道轉(zhuǎn)移的起始點(diǎn)通過拼接點(diǎn)反向時(shí)間積分得到,而控制段上的轉(zhuǎn)移點(diǎn)通過拼接點(diǎn)正向時(shí)間積分得到。地月軌道轉(zhuǎn)移過程如圖6所示。

    在整個(gè)地月轉(zhuǎn)移過程中,本文方法所需施加的脈沖為2~4次。必要的兩次脈沖為航天器脫離地球附近周期區(qū)域和進(jìn)入月球附近周期區(qū)域所施加的。在控制段和拼接區(qū)域F的脈沖可以根據(jù)不同情況而改變。在控制段,如果從轉(zhuǎn)移點(diǎn)正向時(shí)間積分能夠達(dá)到月球附近,則地月軌道轉(zhuǎn)移完成,否則在轉(zhuǎn)移點(diǎn)施加沖量轉(zhuǎn)變?yōu)榭刂贫沃行狞c(diǎn),再完成地月軌道轉(zhuǎn)移。在拼接區(qū)域,如果目標(biāo)域和控制段的延展方向在區(qū)域F相交,則這段軌跡的拼接不需要施加沖量,否則在龐加萊截面圖上選擇控制段與區(qū)域F橫坐標(biāo)相同,縱坐標(biāo)距離最近的點(diǎn)進(jìn)行拼接,施加沖量的大小由縱坐標(biāo)間的距離決定。

    3 數(shù)值仿真校驗(yàn)

    地月軌道轉(zhuǎn)移就是對(duì)地球附近的周期軌道上的航天器施加脈沖,使其飛行至月球附近,然后再對(duì)其施加脈沖使其繞月飛行。在地月三體系統(tǒng)中,地月轉(zhuǎn)移軌道的初始點(diǎn)選擇與文獻(xiàn)[1-4]相同的P點(diǎn)。該點(diǎn)位于半徑為59669 km的圍繞地球旋轉(zhuǎn)的圓軌道上,在圖3所示的龐加萊截面圖中用‘*’表示。在雅可比能量C=7.17218時(shí)航天器通過P點(diǎn)圍繞著地球做周期運(yùn)動(dòng)(由于參數(shù)精度選擇有差異,本文所得到的周期軌道雅可比能量與文獻(xiàn)[1-4]中的C=7.1738不一致),施加一個(gè)平行的脈沖使航天器進(jìn)入混沌區(qū)域,這時(shí)系統(tǒng)的雅可比能量為C=3.17948(實(shí)現(xiàn)地月軌道轉(zhuǎn)移的必要條件是系統(tǒng)雅可比能量必須小于CL=3.1883)。這要求速度的改變?chǔ)=748.06 m/s(文獻(xiàn)[1-4]中ΔV=744.4 m/s),而如圖3所示,點(diǎn)P在C=3.17948時(shí)包圍在混沌區(qū)域內(nèi),當(dāng)軌道轉(zhuǎn)移到月球附近并距離繞月周期軌道為中心的KAM環(huán)Q距離很近時(shí),便可以通過一個(gè)很小的脈沖使軌道流形進(jìn)入KAM環(huán)的邊緣而繞月運(yùn)動(dòng)。

    當(dāng)系統(tǒng)雅可比能量C=3.17948,航天器在不施加外力的情況下從龐加萊截面圖上的P點(diǎn)滑行到月球附近需要27年多的時(shí)間,而且航天器并沒有長(zhǎng)時(shí)間停留在月球附近,而是在月球附近飛行數(shù)圈后便很快遠(yuǎn)離。圖7(a)為航天器在系統(tǒng)雅可比能量C=3.17948下從P點(diǎn)自由滑行的軌跡,圖7(b)為其對(duì)應(yīng)的龐加萊截面圖,對(duì)應(yīng)圖1中龐加萊截面圖的分區(qū)??梢钥闯觯教炱髟邶嫾尤R截面圖中的A層游蕩很長(zhǎng)時(shí)間后,其軌跡轉(zhuǎn)變?yōu)樵贏層和B層之間飛行;然后經(jīng)過相當(dāng)長(zhǎng)的時(shí)間,轉(zhuǎn)變?yōu)樵贐層和C層之間飛行;最后通過C層軌跡轉(zhuǎn)移到月球附近。

    本文選擇P點(diǎn)為中心的區(qū)域作為目標(biāo)域,用本文方法確定控制段,規(guī)定預(yù)定時(shí)間TM和TF,進(jìn)而得到地月低能轉(zhuǎn)移軌道。令TM=5T,圖8(a),圖8(b)分別為預(yù)定時(shí)間取TF=50T和TF=25T得到的兩種轉(zhuǎn)移軌道。圖8(a)中,航天器從距離P點(diǎn)非常接近的P1點(diǎn)出發(fā),施加ΔVP1=736.3 m/s的脈沖使航天器進(jìn)入混沌區(qū)域。在控制段需要ΔVB1=0.002 m/s,使從B1出發(fā)的流形與從P1出發(fā)的流形拼接形成地月轉(zhuǎn)移軌道。在離繞月準(zhǔn)周期軌道KAM最近的點(diǎn)A1施加一個(gè)速度增量使軌道流形進(jìn)入KAM環(huán)內(nèi)需要ΔVA1=4.6 m/s。將這3個(gè)脈沖相加,可得到在混沌軌道的控制中所需要的總能量為740.9 m/s,飛行時(shí)間為303天。圖8(b)中,航天器從距離P點(diǎn)非常接近的P2點(diǎn)出發(fā),施加ΔVP2=752.4 m/s的脈沖使航天器進(jìn)入混沌區(qū)域。在控制段需要ΔVB2=0.2 m/s,使從B2出發(fā)的流形與從P2出發(fā)的流形拼接形成地月轉(zhuǎn)移軌道。在拼接點(diǎn)C2需要ΔVC2=141.6 m/s,在離繞月準(zhǔn)周期軌道KAM最近的點(diǎn)A2施加一個(gè)速度增量使軌道流形進(jìn)入KAM環(huán)內(nèi)需要ΔVA2=2.9 m/s。將這4個(gè)脈沖相加,可得到在混沌軌道的控制中所需要的總能量為

    897.1 m/s,飛行時(shí)間為184天。與現(xiàn)有文獻(xiàn)相比,通過本文方法實(shí)現(xiàn)地月軌道轉(zhuǎn)移所需的脈沖次數(shù)較少,具體的軌跡可以根據(jù)轉(zhuǎn)移過程中需要的總能量和飛行時(shí)間的要求來選擇。表1為文獻(xiàn)[1-4]與圖8中軌道所需總能量、轉(zhuǎn)移時(shí)間和脈沖次數(shù)的對(duì)比。與文獻(xiàn)[1-4]相比,圖8(a)中所設(shè)計(jì)的軌道使用最小的能耗得到了飛行時(shí)間較短的轉(zhuǎn)移軌道,而圖8(b)中所設(shè)計(jì)的軌道通過稍大的能耗顯著縮短了飛行時(shí)間。

    表1 地月軌道轉(zhuǎn)移所需的總能量、飛行時(shí)間和脈沖次數(shù)Table 1 Required energy, time of flight and impulse times in Earth-Moon orbit transfer

    4 結(jié) 論

    本文通過將地月三體系統(tǒng)的龐加萊截面圖分為不同的層,得到了地月系統(tǒng)在實(shí)現(xiàn)地月轉(zhuǎn)移軌道時(shí)耗時(shí)過長(zhǎng)的原因,提出地月軌道轉(zhuǎn)移時(shí)間優(yōu)化原則。結(jié)合向前-向后法,將從地球附近目標(biāo)域正向時(shí)間積分的軌跡和L1瓶頸處的控制段軌跡進(jìn)行拼接進(jìn)而完成低能地月軌道轉(zhuǎn)移。該方法僅需施加2~4次脈沖,在盡可能節(jié)省能量的情況下,有效地減少了航天器在混沌區(qū)域的滑行時(shí)間。而轉(zhuǎn)移所需的能耗和時(shí)間的折衷標(biāo)準(zhǔn)是進(jìn)一步需要完成工作。

    猜你喜歡
    龐加萊航天器月球
    2022 年第二季度航天器發(fā)射統(tǒng)計(jì)
    到月球上“飆車”
    軍事文摘(2022年12期)2022-07-13 03:11:44
    龐加萊偵察術(shù)
    中外文摘(2022年7期)2022-05-17 09:36:42
    陪我去月球
    飛碟探索(2022年1期)2022-03-12 14:11:20
    龐加萊偵察術(shù)
    青年文摘(2021年20期)2021-12-11 18:45:12
    龐加萊偵查術(shù)
    月球上的另一個(gè)我
    軍事文摘(2021年18期)2021-12-02 01:28:54
    2019 年第二季度航天器發(fā)射統(tǒng)計(jì)
    2018 年第三季度航天器發(fā)射統(tǒng)計(jì)
    2018年第二季度航天器發(fā)射統(tǒng)計(jì)
    日韩成人在线观看一区二区三区| 麻豆av在线久日| www.自偷自拍.com| 国产精品亚洲一级av第二区| 午夜免费鲁丝| 国产精品久久电影中文字幕 | 性色av乱码一区二区三区2| tube8黄色片| 亚洲成人手机| 亚洲免费av在线视频| 久久精品熟女亚洲av麻豆精品| 国产精品熟女久久久久浪| 少妇的丰满在线观看| 国内毛片毛片毛片毛片毛片| 乱人伦中国视频| 9191精品国产免费久久| 人人妻人人澡人人爽人人夜夜| 侵犯人妻中文字幕一二三四区| 亚洲综合色网址| 国产精品久久久久久人妻精品电影 | 99国产精品免费福利视频| 丰满人妻熟妇乱又伦精品不卡| 欧美激情极品国产一区二区三区| 国产aⅴ精品一区二区三区波| 国产在线免费精品| 国产亚洲精品一区二区www | 一区在线观看完整版| 亚洲欧美日韩另类电影网站| 黄频高清免费视频| 97人妻天天添夜夜摸| 国产成人精品无人区| 精品国产国语对白av| 国产成人欧美在线观看 | 法律面前人人平等表现在哪些方面| 少妇裸体淫交视频免费看高清 | 80岁老熟妇乱子伦牲交| 亚洲欧美激情在线| 91老司机精品| bbb黄色大片| 51午夜福利影视在线观看| 亚洲精品一二三| 人妻一区二区av| 一进一出好大好爽视频| h视频一区二区三区| 老熟妇乱子伦视频在线观看| 五月天丁香电影| 精品高清国产在线一区| 性少妇av在线| 黄频高清免费视频| 亚洲精品国产精品久久久不卡| 嫩草影视91久久| 亚洲专区字幕在线| 国产高清国产精品国产三级| 夫妻午夜视频| 国产成人免费观看mmmm| 亚洲精品自拍成人| 国产男女超爽视频在线观看| tube8黄色片| 欧美+亚洲+日韩+国产| 亚洲午夜精品一区,二区,三区| 亚洲人成伊人成综合网2020| 欧美精品一区二区大全| 91麻豆精品激情在线观看国产 | 精品高清国产在线一区| 亚洲国产av新网站| 亚洲一区中文字幕在线| 久久午夜亚洲精品久久| 亚洲av美国av| 久久精品国产a三级三级三级| 女人爽到高潮嗷嗷叫在线视频| 成人三级做爰电影| 久久久国产精品麻豆| 亚洲精品一卡2卡三卡4卡5卡| av超薄肉色丝袜交足视频| 女性被躁到高潮视频| 69av精品久久久久久 | a级毛片黄视频| 黄频高清免费视频| 午夜福利视频精品| 亚洲七黄色美女视频| 少妇 在线观看| 99香蕉大伊视频| 三上悠亚av全集在线观看| 色尼玛亚洲综合影院| 成年人黄色毛片网站| 久久久国产精品麻豆| 日韩精品免费视频一区二区三区| 日韩中文字幕欧美一区二区| 成人精品一区二区免费| 夫妻午夜视频| 国产精品秋霞免费鲁丝片| 人人妻人人澡人人爽人人夜夜| 丝袜喷水一区| 亚洲专区中文字幕在线| 久久久精品区二区三区| 我要看黄色一级片免费的| 精品国内亚洲2022精品成人 | 2018国产大陆天天弄谢| 亚洲欧美一区二区三区黑人| 露出奶头的视频| 蜜桃在线观看..| 国产不卡一卡二| 美女视频免费永久观看网站| 黑人操中国人逼视频| 高清在线国产一区| 视频在线观看一区二区三区| 精品国产乱码久久久久久男人| 十八禁网站免费在线| 一区福利在线观看| 亚洲国产欧美日韩在线播放| 午夜91福利影院| 99九九在线精品视频| 精品久久久精品久久久| 日日爽夜夜爽网站| 中文字幕色久视频| 精品国产一区二区久久| 男女边摸边吃奶| 国产在线免费精品| 亚洲久久久国产精品| 精品一区二区三卡| 丰满饥渴人妻一区二区三| 啦啦啦中文免费视频观看日本| aaaaa片日本免费| 汤姆久久久久久久影院中文字幕| 99精国产麻豆久久婷婷| 国产成人精品在线电影| 另类亚洲欧美激情| av线在线观看网站| 999久久久国产精品视频| 久久久精品区二区三区| 又黄又粗又硬又大视频| 国产成人av激情在线播放| 亚洲 欧美一区二区三区| 黑人欧美特级aaaaaa片| www.熟女人妻精品国产| av国产精品久久久久影院| 成人永久免费在线观看视频 | 国产精品久久电影中文字幕 | 9191精品国产免费久久| xxxhd国产人妻xxx| 男女高潮啪啪啪动态图| 精品欧美一区二区三区在线| 欧美乱妇无乱码| 一边摸一边抽搐一进一小说 | 女性被躁到高潮视频| av又黄又爽大尺度在线免费看| 亚洲人成电影免费在线| 亚洲免费av在线视频| 亚洲熟女精品中文字幕| 一进一出抽搐动态| 精品福利永久在线观看| 大型av网站在线播放| 成年人黄色毛片网站| 日本vs欧美在线观看视频| 多毛熟女@视频| 看免费av毛片| 久久精品亚洲精品国产色婷小说| 在线观看66精品国产| 嫩草影视91久久| 国产又色又爽无遮挡免费看| 亚洲久久久国产精品| 欧美日韩亚洲综合一区二区三区_| 日本a在线网址| 日本av免费视频播放| 久久香蕉激情| 午夜激情久久久久久久| 国产高清国产精品国产三级| 亚洲国产欧美网| 久久人人爽av亚洲精品天堂| 国产亚洲精品一区二区www | 亚洲色图av天堂| 巨乳人妻的诱惑在线观看| 国产成人精品久久二区二区91| 国产激情久久老熟女| 国产精品一区二区精品视频观看| 午夜福利,免费看| 久久精品亚洲av国产电影网| 亚洲 欧美一区二区三区| 午夜视频精品福利| 成在线人永久免费视频| 2018国产大陆天天弄谢| 久热爱精品视频在线9| 久热这里只有精品99| 精品国产乱码久久久久久小说| 老司机亚洲免费影院| 高清视频免费观看一区二区| 国产男女内射视频| 好男人电影高清在线观看| 久久久久久久精品吃奶| 国产成人精品久久二区二区免费| 亚洲男人天堂网一区| 免费不卡黄色视频| 国产精品久久久久久精品古装| 久久中文字幕一级| xxxhd国产人妻xxx| av视频免费观看在线观看| 久久香蕉激情| 一边摸一边抽搐一进一出视频| 精品人妻熟女毛片av久久网站| 国产成人欧美在线观看 | 女同久久另类99精品国产91| 成年人免费黄色播放视频| 精品久久久久久电影网| 色婷婷久久久亚洲欧美| 少妇被粗大的猛进出69影院| 久久性视频一级片| 精品高清国产在线一区| 国产在线观看jvid| 午夜福利在线观看吧| 青青草视频在线视频观看| 中文字幕制服av| 国产一区二区激情短视频| 精品久久久久久电影网| 在线观看免费高清a一片| 日韩欧美国产一区二区入口| 天堂中文最新版在线下载| 91麻豆av在线| 天堂俺去俺来也www色官网| 日韩三级视频一区二区三区| 国产av又大| 热re99久久国产66热| 波多野结衣av一区二区av| 交换朋友夫妻互换小说| 国产99久久九九免费精品| 97在线人人人人妻| 精品国产一区二区三区四区第35| 电影成人av| 精品国产乱子伦一区二区三区| www.自偷自拍.com| 久久精品国产亚洲av高清一级| 亚洲人成伊人成综合网2020| 国产精品国产高清国产av | 亚洲第一青青草原| 久9热在线精品视频| 久久精品国产综合久久久| 亚洲第一av免费看| 在线观看免费日韩欧美大片| 国产在视频线精品| 国产精品麻豆人妻色哟哟久久| 亚洲七黄色美女视频| 久久午夜综合久久蜜桃| 女人高潮潮喷娇喘18禁视频| 欧美日韩成人在线一区二区| 18在线观看网站| 国产精品 国内视频| 18禁国产床啪视频网站| 亚洲av国产av综合av卡| 国产精品 欧美亚洲| 男男h啪啪无遮挡| 婷婷成人精品国产| 一进一出好大好爽视频| 日韩大码丰满熟妇| 性少妇av在线| 丰满迷人的少妇在线观看| 啦啦啦在线免费观看视频4| 搡老熟女国产l中国老女人| 亚洲欧美激情在线| 久久中文看片网| 精品久久久久久久毛片微露脸| 中文字幕人妻熟女乱码| 一区在线观看完整版| 丰满人妻熟妇乱又伦精品不卡| 在线十欧美十亚洲十日本专区| 乱人伦中国视频| 国产一区二区激情短视频| 日韩人妻精品一区2区三区| 国产男女内射视频| 国产老妇伦熟女老妇高清| 亚洲欧美色中文字幕在线| 国内毛片毛片毛片毛片毛片| 日本撒尿小便嘘嘘汇集6| 日韩中文字幕视频在线看片| 日韩免费高清中文字幕av| 纯流量卡能插随身wifi吗| 欧美黑人精品巨大| 纵有疾风起免费观看全集完整版| 亚洲国产欧美在线一区| 超碰成人久久| 亚洲美女黄片视频| 69av精品久久久久久 | 又紧又爽又黄一区二区| 欧美黄色片欧美黄色片| 欧美人与性动交α欧美精品济南到| 精品卡一卡二卡四卡免费| 久热爱精品视频在线9| 久久精品aⅴ一区二区三区四区| 国产福利在线免费观看视频| 操出白浆在线播放| 国产精品久久久久成人av| av欧美777| 欧美成狂野欧美在线观看| 真人做人爱边吃奶动态| 国产亚洲精品第一综合不卡| av超薄肉色丝袜交足视频| 精品少妇久久久久久888优播| 精品人妻1区二区| 交换朋友夫妻互换小说| 91老司机精品| 99久久人妻综合| 大型av网站在线播放| 黄片播放在线免费| 十八禁高潮呻吟视频| 国产精品久久久久久人妻精品电影 | 婷婷丁香在线五月| 久热爱精品视频在线9| 亚洲全国av大片| 久久国产精品大桥未久av| 亚洲伊人久久精品综合| 精品少妇内射三级| 丁香六月天网| 黄频高清免费视频| 法律面前人人平等表现在哪些方面| 黑人巨大精品欧美一区二区mp4| 精品国产一区二区三区四区第35| 男女免费视频国产| 国产单亲对白刺激| 欧美日韩av久久| 少妇被粗大的猛进出69影院| 亚洲人成电影观看| av天堂久久9| 成人特级黄色片久久久久久久 | 久久天堂一区二区三区四区| 久久久久国内视频| 久久热在线av| 麻豆成人av在线观看| 在线永久观看黄色视频| 久久久久国内视频| 国产男女内射视频| 国产老妇伦熟女老妇高清| 精品一品国产午夜福利视频| 亚洲七黄色美女视频| 首页视频小说图片口味搜索| 看免费av毛片| 天天躁夜夜躁狠狠躁躁| 欧美日本中文国产一区发布| 午夜精品久久久久久毛片777| 91老司机精品| 嫩草影视91久久| h视频一区二区三区| 精品亚洲乱码少妇综合久久| 免费高清在线观看日韩| 久久久久视频综合| 曰老女人黄片| 巨乳人妻的诱惑在线观看| 国产伦人伦偷精品视频| 老司机午夜福利在线观看视频 | 国产1区2区3区精品| 老熟妇仑乱视频hdxx| 国产精品九九99| 日本vs欧美在线观看视频| 国产成人免费无遮挡视频| 亚洲精品久久成人aⅴ小说| 亚洲国产av新网站| 久久久久久久国产电影| 欧美精品一区二区免费开放| 欧美日韩av久久| 黑丝袜美女国产一区| 天堂8中文在线网| 国产在线观看jvid| 日本撒尿小便嘘嘘汇集6| 超碰97精品在线观看| 露出奶头的视频| 我的亚洲天堂| 蜜桃在线观看..| 成人国产av品久久久| 国产精品免费一区二区三区在线 | 日韩欧美免费精品| 中文字幕av电影在线播放| 欧美变态另类bdsm刘玥| 成人三级做爰电影| 午夜福利视频在线观看免费| 国产熟女午夜一区二区三区| 欧美成人午夜精品| 亚洲欧美一区二区三区黑人| 大片免费播放器 马上看| 99精品欧美一区二区三区四区| 色综合欧美亚洲国产小说| av线在线观看网站| 欧美成狂野欧美在线观看| 黄色视频在线播放观看不卡| 久久免费观看电影| 激情视频va一区二区三区| 欧美乱妇无乱码| 久久ye,这里只有精品| 757午夜福利合集在线观看| 天天躁夜夜躁狠狠躁躁| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 欧美日韩一级在线毛片| 国产欧美日韩一区二区三区在线| 无人区码免费观看不卡 | 操美女的视频在线观看| 国产xxxxx性猛交| 午夜激情久久久久久久| 亚洲精品久久成人aⅴ小说| 丰满人妻熟妇乱又伦精品不卡| 一进一出好大好爽视频| 黄色怎么调成土黄色| 建设人人有责人人尽责人人享有的| 亚洲精品美女久久av网站| 国产日韩一区二区三区精品不卡| www.999成人在线观看| 久久热在线av| 五月开心婷婷网| 9191精品国产免费久久| 亚洲七黄色美女视频| 亚洲色图av天堂| 亚洲av国产av综合av卡| 国产精品久久久久久精品古装| 最黄视频免费看| 99久久人妻综合| 成人国语在线视频| av欧美777| 国产人伦9x9x在线观看| 亚洲熟妇熟女久久| 搡老熟女国产l中国老女人| 久久人妻熟女aⅴ| 久久久久网色| 免费人妻精品一区二区三区视频| 蜜桃国产av成人99| 菩萨蛮人人尽说江南好唐韦庄| 色在线成人网| 熟女少妇亚洲综合色aaa.| 国产精品久久久人人做人人爽| 丰满迷人的少妇在线观看| 天堂动漫精品| 最新美女视频免费是黄的| 国产成人系列免费观看| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产亚洲av香蕉五月 | 脱女人内裤的视频| 极品教师在线免费播放| 欧美日韩精品网址| 制服诱惑二区| 乱人伦中国视频| 久久久久国产一级毛片高清牌| 一区二区日韩欧美中文字幕| 成人18禁高潮啪啪吃奶动态图| 捣出白浆h1v1| 亚洲av片天天在线观看| 国产福利在线免费观看视频| 亚洲中文字幕日韩| 一进一出抽搐动态| 巨乳人妻的诱惑在线观看| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人欧美精品刺激| 精品国产国语对白av| 最近最新中文字幕大全免费视频| 亚洲三区欧美一区| 精品国产一区二区久久| 国产又色又爽无遮挡免费看| 黑人欧美特级aaaaaa片| 国产男女超爽视频在线观看| 成年版毛片免费区| 欧美成狂野欧美在线观看| 亚洲精品在线观看二区| 高清毛片免费观看视频网站 | 久久中文看片网| 精品欧美一区二区三区在线| 蜜桃国产av成人99| 一边摸一边抽搐一进一小说 | 国产精品久久久久成人av| 五月开心婷婷网| 精品国内亚洲2022精品成人 | 亚洲美女黄片视频| bbb黄色大片| 女人被躁到高潮嗷嗷叫费观| 又大又爽又粗| 免费日韩欧美在线观看| 俄罗斯特黄特色一大片| 男男h啪啪无遮挡| 欧美激情极品国产一区二区三区| 欧美大码av| 久久狼人影院| 国产成人欧美在线观看 | 一二三四社区在线视频社区8| 真人做人爱边吃奶动态| 久久精品国产亚洲av高清一级| 日韩一区二区三区影片| 免费在线观看黄色视频的| videos熟女内射| 黑人猛操日本美女一级片| 亚洲精品美女久久久久99蜜臀| 国产成人精品久久二区二区免费| 亚洲av欧美aⅴ国产| 亚洲精品国产区一区二| 亚洲精品乱久久久久久| 亚洲,欧美精品.| 国产成人精品久久二区二区免费| 亚洲黑人精品在线| 纯流量卡能插随身wifi吗| 无人区码免费观看不卡 | 亚洲欧美日韩高清在线视频 | 国产精品99久久99久久久不卡| 一区在线观看完整版| 午夜福利视频在线观看免费| 成人18禁高潮啪啪吃奶动态图| 人人妻人人澡人人爽人人夜夜| 亚洲午夜精品一区,二区,三区| tocl精华| 无遮挡黄片免费观看| 欧美大码av| 免费高清在线观看日韩| 侵犯人妻中文字幕一二三四区| bbb黄色大片| 国产又色又爽无遮挡免费看| av国产精品久久久久影院| 一夜夜www| 热re99久久精品国产66热6| 桃花免费在线播放| 他把我摸到了高潮在线观看 | 母亲3免费完整高清在线观看| 悠悠久久av| 91麻豆av在线| 久久久水蜜桃国产精品网| 热99国产精品久久久久久7| 欧美日韩成人在线一区二区| 91成人精品电影| 国产aⅴ精品一区二区三区波| 亚洲av日韩精品久久久久久密| 亚洲午夜精品一区,二区,三区| 国产熟女午夜一区二区三区| 99re在线观看精品视频| av国产精品久久久久影院| 狠狠狠狠99中文字幕| 999精品在线视频| 国产成人欧美在线观看 | 久久久国产成人免费| 亚洲第一av免费看| 啦啦啦中文免费视频观看日本| 亚洲国产精品一区二区三区在线| 黄色 视频免费看| 欧美成人免费av一区二区三区 | 一夜夜www| 国产精品久久久久久精品古装| 狂野欧美激情性xxxx| 深夜精品福利| 欧美 亚洲 国产 日韩一| 高清毛片免费观看视频网站 | 欧美国产精品va在线观看不卡| 亚洲色图 男人天堂 中文字幕| 国产日韩欧美在线精品| 国产成人av激情在线播放| 51午夜福利影视在线观看| 丝袜在线中文字幕| 成人三级做爰电影| 亚洲精品国产精品久久久不卡| www日本在线高清视频| 亚洲色图 男人天堂 中文字幕| 日本欧美视频一区| 一级毛片精品| 亚洲国产欧美一区二区综合| 久久精品国产a三级三级三级| 亚洲欧洲精品一区二区精品久久久| 99精国产麻豆久久婷婷| 热re99久久精品国产66热6| 精品欧美一区二区三区在线| 黄色视频不卡| 欧美精品亚洲一区二区| 日本vs欧美在线观看视频| 精品乱码久久久久久99久播| 亚洲自偷自拍图片 自拍| 国产欧美日韩一区二区精品| 成年动漫av网址| 精品少妇内射三级| 嫁个100分男人电影在线观看| 少妇精品久久久久久久| 高潮久久久久久久久久久不卡| 国产成人精品无人区| 亚洲 国产 在线| 中文字幕最新亚洲高清| 国内毛片毛片毛片毛片毛片| 久久婷婷成人综合色麻豆| 国产激情久久老熟女| 亚洲精品粉嫩美女一区| 国产精品1区2区在线观看. | 91大片在线观看| 天堂动漫精品| 97人妻天天添夜夜摸| 99久久精品国产亚洲精品| 美女高潮喷水抽搐中文字幕| 国产老妇伦熟女老妇高清| 午夜老司机福利片| 日韩欧美三级三区| 国产老妇伦熟女老妇高清| 三级毛片av免费| 亚洲第一青青草原| 久久久久网色| 菩萨蛮人人尽说江南好唐韦庄| 人妻一区二区av| 我要看黄色一级片免费的| 欧美激情久久久久久爽电影 | 国产主播在线观看一区二区| 久久精品国产综合久久久| 香蕉国产在线看| 日韩成人在线观看一区二区三区| 国产av国产精品国产| 亚洲av国产av综合av卡| 在线观看一区二区三区激情| 91精品国产国语对白视频| 麻豆av在线久日| 国产在线观看jvid| 黄色a级毛片大全视频| 久久人妻福利社区极品人妻图片| 亚洲欧美精品综合一区二区三区| 国产成人精品久久二区二区91| 一本色道久久久久久精品综合| 久久精品国产亚洲av香蕉五月 | 亚洲欧美一区二区三区久久| 久久久久久久久久久久大奶| 久久久久久久大尺度免费视频| av又黄又爽大尺度在线免费看| 99国产精品免费福利视频| 两性午夜刺激爽爽歪歪视频在线观看 | 我的亚洲天堂|