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

    載人月面著陸器飛行策略?xún)?yōu)化研究

    2023-10-21 01:31:36孫興亮朱恩涌
    宇航學(xué)報(bào) 2023年9期
    關(guān)鍵詞:環(huán)月著陸器交會(huì)

    黃 震,孫興亮,余 歡,朱恩涌

    (北京空間飛行器總體設(shè)計(jì)部,北京 100094)

    0 引 言

    載人登月的飛行模式主要有兩種[1]:一種是一次發(fā)射、船器組合奔月模式,即載人月面著陸器與載人飛船共同飛行至環(huán)月軌道,如美國(guó)“阿波羅”計(jì)劃、蘇聯(lián)登月計(jì)劃[2];另一種是分別發(fā)射、船器獨(dú)立奔月模式,即載人月面著陸器提前飛行至月球軌道,在該軌道與載人飛船交會(huì)對(duì)接,實(shí)現(xiàn)落月前的人員和物資轉(zhuǎn)移,如美國(guó)“阿爾忒彌斯”計(jì)劃[3]、中國(guó)載人月球探測(cè)工程[4]。

    無(wú)論哪種飛行模式,載人登月飛行器的速度增量需求都非常大。其中,載人月面著陸器一般負(fù)責(zé)近月制動(dòng)、月面下降著陸和月面起飛,其總速度增量將超過(guò)5 000 m/s,約為近地載人飛行器的15倍,如阿波羅任務(wù)的載人月面著陸器推進(jìn)劑與發(fā)射質(zhì)量的占比達(dá)到了71.5%[5]。

    載人月面著陸器的設(shè)計(jì)難點(diǎn)在于重量控制[6],從圖1中不難看出,重量控制的重點(diǎn)在于控制推進(jìn)劑的使用量。根據(jù)飛行器自身特征,對(duì)任務(wù)進(jìn)行合理的規(guī)劃,選擇優(yōu)化的飛行策略是載人月面著陸器最為關(guān)鍵的設(shè)計(jì)因素。

    圖1 阿波羅著陸器質(zhì)量組成Fig.1 Mass allocation of the Apollo lunar lander

    1 載人月面著陸器及其飛行過(guò)程

    載人月面著陸器質(zhì)量約26 t[7],由登月艙和推進(jìn)艙組成,具備自主控制飛行能力,可將兩名航天員、月球車(chē)、相關(guān)載荷送達(dá)月球表面,采用輕量化設(shè)計(jì)提高了結(jié)構(gòu)效率,并具備動(dòng)力冗余功能來(lái)保證航天員的安全。載人月面著陸器是航天員登陸月球后的月面生活中心、能源中心及數(shù)據(jù)中心,支持開(kāi)展月面駐留和月面活動(dòng),如圖2所示[7]。

    圖2 載人月面著陸器構(gòu)型圖Fig.2 Configuration of the manned lunar lander

    載人登月的初步方案是[4]:采用兩枚運(yùn)載火箭分別將載人月面著陸器和載人飛船送至地月轉(zhuǎn)移軌道,載人飛船和載人月面著陸器在環(huán)月軌道交會(huì)對(duì)接,航天員從載人飛船進(jìn)入載人月面著陸器。之后,載人月面著陸器將下降著陸于月面預(yù)定區(qū)域,航天員登陸月球開(kāi)展科學(xué)考察與樣品采集。完成既定任務(wù)后,航天員乘坐載人月面著陸器上升至環(huán)月軌道與載人飛船交會(huì)對(duì)接,并攜帶樣品乘坐載人飛船返回地球。

    按照該載人登月方案,載人月面著陸器將主要經(jīng)歷以下飛行階段:

    1)地月轉(zhuǎn)移飛行

    載人月面著陸器從地球飛往月球的過(guò)程,一般為3~5天,期間需要進(jìn)行中途修正和近月制動(dòng),速度增量約為850 m/s[8-9]。

    2)環(huán)月飛行

    載人月面著陸器近月制動(dòng)后,環(huán)繞月球長(zhǎng)期飛行,等待載人飛船共同執(zhí)行登月任務(wù)。環(huán)月飛行期間,軌道面和軌道高度需要維持,總速度增量不宜超過(guò)50 m/s。

    3)動(dòng)力下降與著陸

    由于月球沒(méi)有大氣,航天員搭乘載人月面著陸器著陸到月球的過(guò)程主要依靠發(fā)動(dòng)機(jī)反推進(jìn)行減速和著陸,等效速度增量約2 000 m/s[9-10],為推進(jìn)劑消耗最多的飛行階段。

    4)起飛上升

    航天員搭乘載人月面著陸器上升入軌的過(guò)程主要依靠發(fā)動(dòng)機(jī)提供發(fā)射動(dòng)力,等效速度增量約1 850 m/s[8],推進(jìn)劑消耗量占比大。

    5)環(huán)月交會(huì)對(duì)接

    載人月面著陸器與載人飛船在月球軌道交會(huì)對(duì)接的過(guò)程中,速度增量約為80 m/s[11]。由于是任務(wù)末端,干重對(duì)推進(jìn)劑用量的影響非常顯著,可以達(dá)到1∶10左右。

    根據(jù)載人月面著陸器的兩艙功能分配,近月制動(dòng)和動(dòng)力下降主減速的推進(jìn)劑消耗由推進(jìn)艙提供[7],其余飛行階段的推進(jìn)劑消耗由登月艙提供。

    本文在討論各種優(yōu)化策略時(shí),給出了一些仿真算例,采用的主要輸入?yún)?shù)說(shuō)明如下。載人月面著陸器發(fā)射質(zhì)量取26 t,兩艙軌控發(fā)動(dòng)機(jī)比沖均取300 s,根據(jù)各飛行階段的速度增量需求,首次與載人飛船交會(huì)對(duì)接前的整器質(zhì)量取19.5 t,月面下降前的整器質(zhì)量取19.9 t,登月艙月面起飛前質(zhì)量取8 t。

    2 月面著陸器正常飛行策略?xún)?yōu)化

    本文根據(jù)不同飛行階段的特點(diǎn),提出飛行全過(guò)程的優(yōu)化思路,從任務(wù)規(guī)劃、軌道特性、實(shí)施方案等多個(gè)角度優(yōu)化載人月面著陸器推進(jìn)劑使用量。

    2.1 基于最低發(fā)射窗口的地月轉(zhuǎn)移設(shè)計(jì)

    飛行器地月轉(zhuǎn)移設(shè)計(jì)通常有一般轉(zhuǎn)移軌道、自由返回軌道和混合軌道3種。為適應(yīng)不同登月點(diǎn)和不同環(huán)月軌道傾角的任務(wù)需求,發(fā)射時(shí)為無(wú)人狀態(tài)的月面著陸器可采用一般轉(zhuǎn)移軌道。這種策略下,載人月面著陸器到達(dá)任意傾角環(huán)月軌道的近月制動(dòng)速度增量基本保持不變,實(shí)現(xiàn)全月面任意登月點(diǎn)可達(dá),任務(wù)適應(yīng)能力強(qiáng)、推進(jìn)劑消耗小。不同環(huán)月軌道傾角下的近月制動(dòng)速度增量算例見(jiàn)圖3。

    圖3 不同環(huán)月軌道傾角下近月制動(dòng)速度增量算例Fig.3 Velocity increments for lunar orbit insertion in lunar orbits with different inclinations

    近月制動(dòng)速度增量隨地月相對(duì)位置關(guān)系、地月轉(zhuǎn)移的升降軌方式變化,即:

    ΔvLOI(sC,tL)≈ΔvLOI(sC,tL+TM)

    (1)

    式中:ΔvLOI為近月制動(dòng)速度增量;sC為升軌或降軌出發(fā)標(biāo)志;tL為發(fā)射日期;TM為恒星月,約27.3天。對(duì)于同一種近地出發(fā)方式,ΔvLOI隨發(fā)射日期近似呈周期變化。

    根據(jù)近月制動(dòng)特性分析,綜合考慮推進(jìn)劑消耗代價(jià)、任務(wù)規(guī)劃靈活性設(shè)計(jì)發(fā)射窗口,具體方法建模如下:

    (2)

    式中:Δvmax為近月制動(dòng)速度增量上限約束;NT為一個(gè)恒星月內(nèi)滿(mǎn)足速度增量約束的發(fā)射窗口數(shù)量;Nmin為單月發(fā)射窗口數(shù)量下限約束,一般為2~3天。

    以2029年發(fā)射的5天一般轉(zhuǎn)移軌道為例,ΔvLOI的變化情況見(jiàn)圖4。取Δvmax為810 m/s,則NT不少于11天;降低Δvmax至800 m/s,則NT不少于8天;進(jìn)一步降低Δvmax至790 m/s,則NT不少于5天,不小于發(fā)射窗口數(shù)量下限。因此,在發(fā)射窗口數(shù)量滿(mǎn)足工程任務(wù)最低需求的前提下,可以?xún)?yōu)化一般轉(zhuǎn)移軌道的近月制動(dòng)速度增量分配,節(jié)約近月制動(dòng)推進(jìn)劑消耗量達(dá)140 kg。

    圖4 一般轉(zhuǎn)移軌道近月制動(dòng)速度增量變化情況Fig.4 Velocity increments for lunar orbit insertion in a general translunar orbit

    2.2 基于軌道攝動(dòng)的環(huán)月設(shè)計(jì)

    載人月面著陸器近月制動(dòng)后的初始環(huán)月軌道升交點(diǎn)經(jīng)度Ω0相對(duì)固定,落月時(shí)刻的軌道升交點(diǎn)經(jīng)度ΩL為

    ΩL=Ω0+(ωM+ωP)NL

    (3)

    式中:NL為近月制動(dòng)至落月的環(huán)月飛行時(shí)長(zhǎng),一般不進(jìn)行大幅調(diào)整,以免影響整個(gè)登月任務(wù)規(guī)劃;ωM為月球自轉(zhuǎn)引起的升交點(diǎn)經(jīng)度西退速率,約-13.2(°)/天;ωP為月球軌道攝動(dòng)引起的升交點(diǎn)赤經(jīng)東進(jìn)速率,對(duì)于相同的近月點(diǎn)位置和軌道傾角,ωP主要與遠(yuǎn)月點(diǎn)高度ha相關(guān)。

    載人登月任務(wù)中,落月當(dāng)圈的環(huán)月軌道升交點(diǎn)經(jīng)度Ωt根據(jù)登月點(diǎn)位置、月面駐留時(shí)長(zhǎng)確定,是由任務(wù)規(guī)劃安排唯一確定的。因此在ΩL偏離Ωt時(shí),需增加軌道機(jī)動(dòng)使ΩL=Ωt。

    為節(jié)約推進(jìn)劑消耗,基于環(huán)月軌道攝動(dòng)機(jī)理,本文提出兩次近月制動(dòng)策略,通過(guò)調(diào)整第1次近月制動(dòng)后的初始環(huán)月軌道遠(yuǎn)月點(diǎn)高度ha0改變?chǔ)豍,進(jìn)而改變?chǔ)?的變化速率,盡可能減小升交點(diǎn)經(jīng)度偏差。初始環(huán)月軌道設(shè)計(jì)可建模為如下優(yōu)化問(wèn)題:

    (4)

    式中:x為優(yōu)化設(shè)計(jì)變量;NL1為2次近月制動(dòng)間的環(huán)月飛行時(shí)長(zhǎng);hat為第2次近月制動(dòng)后進(jìn)入的目標(biāo)環(huán)月軌道遠(yuǎn)月點(diǎn)高度;NL2為第2次近月制動(dòng)至落月的環(huán)月飛行時(shí)長(zhǎng)。

    ΩL的調(diào)整范圍越大,越能夠更接近任務(wù)所需的Ωt,額外用來(lái)調(diào)整升交點(diǎn)經(jīng)度所需的推進(jìn)劑就越少。以目標(biāo)環(huán)月軌道傾角165°的登月任務(wù)為例,不同ha0范圍下的ΩL可達(dá)域?qū)挾纫?jiàn)圖5。按環(huán)月飛行60天考慮,若ha0在2 000~3 500 km范圍調(diào)整,對(duì)應(yīng)的ΩL的可達(dá)域?qū)挾葹?.4°;若ha0的可調(diào)范圍擴(kuò)充至2 000~7 000 km,則對(duì)應(yīng)的ΩL的可達(dá)域?qū)挾葹?6.5°。

    圖5 不同遠(yuǎn)月點(diǎn)高度下的升交點(diǎn)經(jīng)度可達(dá)域Fig.5 Attainable region of longitude of ascent node for lunar orbits with various apolune height

    采用遠(yuǎn)月點(diǎn)高度可調(diào)的大橢圓環(huán)月軌道設(shè)計(jì)策略,可降低環(huán)月飛行期間軌道修正的速度增量需求。若采用1次近月制動(dòng)策略,直接進(jìn)入遠(yuǎn)月點(diǎn)高度hat的目標(biāo)環(huán)月軌道,則需要為升交點(diǎn)經(jīng)度修正多預(yù)留一部分推進(jìn)劑。以5°的升交點(diǎn)經(jīng)度修正量為例,針對(duì)上述165°傾角的目標(biāo)環(huán)月軌道,需多消耗推進(jìn)劑約150 kg。

    2.3 基于兩艙均衡消耗的環(huán)月降軌設(shè)計(jì)

    載人月面著陸器登月艙和推進(jìn)艙均配置了推進(jìn)系統(tǒng),任何一個(gè)艙段的推進(jìn)劑使用量不僅與自身執(zhí)行任務(wù)的速度增量有關(guān),還由于整器質(zhì)量影響,與另一艙段推進(jìn)劑已使用情況有關(guān),兩艙推進(jìn)劑使用量耦合性很強(qiáng),如下式所示:

    (5)

    式中:Δmk為第k次軌控任務(wù)的推進(jìn)劑消耗量;m0為整器初始質(zhì)量;Δvi為第i次軌控任務(wù)的速度增量;Isp為發(fā)動(dòng)機(jī)比沖;g為地球重力加速度。

    假設(shè)第k次軌控任務(wù)由推進(jìn)艙執(zhí)行,第l次軌控任務(wù)(l≤k-1)由登月艙執(zhí)行,根據(jù)上式可得Δmk對(duì)Δvl的敏感度:

    (6)

    式中:mk為第k次軌控任務(wù)后的整器剩余質(zhì)量;Δvl為登月艙執(zhí)行第l次軌控任務(wù)的速度變化量。

    由上式可見(jiàn),Δmk與Δvl負(fù)相關(guān),即由登月艙執(zhí)行的第l次軌控任務(wù)Δvl越小,由推進(jìn)艙執(zhí)行的第k次軌控任務(wù)推進(jìn)劑消耗量Δmk就越大。因此,一般來(lái)說(shuō),登月艙推進(jìn)劑預(yù)算時(shí)需考慮Δvl取極大值的工況,推進(jìn)艙推進(jìn)劑預(yù)算時(shí)則要考慮Δvl取極小值的工況,即兩艙在推進(jìn)劑預(yù)算時(shí)都按照最大需求預(yù)留,造成了資源浪費(fèi)。

    例如,中途修正速度增量ΔvMCC由登月艙執(zhí)行,散布范圍較大,約為0~27 m/s,見(jiàn)圖6。根據(jù)兩艙推進(jìn)劑使用的耦合關(guān)系,該速度增量波動(dòng)使兩艙需要分別預(yù)留240 kg和130 kg以應(yīng)對(duì)中途修正不確定性。

    圖6 中途修正速度增量分布算例Fig.6 Distributions of velocity increments for mid-course correction

    本文提出月面下降兩艙均衡消耗方法,通過(guò)調(diào)整動(dòng)力下降初始高度,達(dá)到充分利用兩艙剩余推進(jìn)劑的目的。由于推進(jìn)艙推進(jìn)劑在動(dòng)力下降結(jié)束后將無(wú)法使用,應(yīng)盡可能實(shí)現(xiàn)推進(jìn)艙推進(jìn)劑耗盡。在ΔvMCC較小時(shí),登月艙剩余推進(jìn)劑較多,整器質(zhì)量大,通過(guò)降低動(dòng)力下降初始高度,減輕推進(jìn)艙動(dòng)力下降負(fù)擔(dān);反之,登月艙剩余推進(jìn)劑較少、整器質(zhì)量輕時(shí),通過(guò)提高動(dòng)力下降初始高度,可以利用推進(jìn)艙剩余推進(jìn)劑能力減輕登月艙環(huán)月降軌負(fù)擔(dān),從而實(shí)現(xiàn)兩艙推進(jìn)劑均衡消耗。

    月面下降兩艙均衡消耗方法可用如下優(yōu)化模型表示:

    (7)

    式中:hPDI為動(dòng)力下降初始高度,為設(shè)計(jì)變量;ΔvLO和ΔvPD1分別為環(huán)月降軌和動(dòng)力下降速度增量,分別由登月艙與推進(jìn)艙提供,與hPDI相關(guān);mD0為環(huán)月降軌前的整器質(zhì)量,與環(huán)月降軌前的整器推進(jìn)劑消耗量ΔmD0相關(guān);ΔmLO和ΔmPD1分別為環(huán)月降軌和動(dòng)力下降的推進(jìn)劑消耗量;mPR1為動(dòng)力下降前推進(jìn)艙推進(jìn)劑可用量,與前期推進(jìn)劑消耗情況相關(guān)。

    假設(shè)ΔvMCC取0 m/s時(shí),hPDI設(shè)計(jì)為15 km,則在ΔvMCC取27 m/s時(shí),推進(jìn)艙可從26 km處執(zhí)行動(dòng)力下降,在充分利用推進(jìn)艙剩余推進(jìn)劑的同時(shí),可節(jié)省登月艙推進(jìn)劑21 kg。

    2.4 基于中途分離的動(dòng)力下降設(shè)計(jì)

    載人月面著陸器的動(dòng)力下降過(guò)程如圖7所示。動(dòng)力下降過(guò)程一般可分為2個(gè)階段:第1個(gè)階段從動(dòng)力下降開(kāi)始至距月面較低高度,主要目標(biāo)是降低環(huán)月軌道高度和速度至目標(biāo)值;第2個(gè)階段從第1個(gè)階段結(jié)束至著陸月面,主要目標(biāo)是瞄準(zhǔn)一條穩(wěn)定、緩速的飛行軌跡,便于執(zhí)行避障所需的導(dǎo)航、地形識(shí)別和機(jī)動(dòng)控制,確保著陸安全性。

    圖7 動(dòng)力下降和動(dòng)力上升階段劃分Fig.7 Stage division for powered descent and powered ascent

    為減少推進(jìn)劑消耗,載人月面著陸器一般在上升前將推進(jìn)劑消耗殆盡的艙段分離,以較小的質(zhì)量規(guī)模上升入軌。分離點(diǎn)可選擇在月面,也可以選擇在下降中途,即動(dòng)力下降階段1和階段2的交班點(diǎn)。由于提前進(jìn)行艙段分離,登月艙執(zhí)行階段2的穩(wěn)定著陸任務(wù)負(fù)重小,任務(wù)效益更高。

    動(dòng)力下降全程的推進(jìn)劑消耗隨兩艙分離高度變化,以推進(jìn)劑消耗最少為目標(biāo)開(kāi)展中途分離高度的優(yōu)化設(shè)計(jì),優(yōu)化模型如下:

    (8)

    式中:hSEP為中途分離高度;x為設(shè)計(jì)變量;mPDI和mSEP分別為動(dòng)力下降前的整器質(zhì)量、中途分離后的登月艙質(zhì)量;mR1為中途分離時(shí)推進(jìn)艙剩余質(zhì)量,包括艙段干重和推進(jìn)劑不可用量;ΔvPD2和ΔmPD2分別為登月艙獨(dú)自下降著陸過(guò)程的速度增量和推進(jìn)劑消耗;ΔmPD為動(dòng)力下降的總推進(jìn)劑消耗;fs為設(shè)計(jì)約束,一般考慮安全性約束,如兩艙落點(diǎn)的安全距離、兩艙分離后低高度飛行的安全避障等。

    動(dòng)力下降的推進(jìn)劑總消耗一般隨中途分離高度升高而增加,主要因?yàn)閯?dòng)力下降第2階段并非沿推進(jìn)劑最優(yōu)軌跡飛行,該階段在整個(gè)動(dòng)力下降軌跡的占比越大,總速度增量越大。圖8給出了動(dòng)力下降的推進(jìn)劑總消耗量隨中途分離高度的變化算例,圖中縱坐標(biāo)表示的是不同分離高度的推進(jìn)劑消耗量相對(duì)于2 km分離高度的變化量。由結(jié)果可見(jiàn),若中途分離高度從3 km降至2 km,推進(jìn)劑總消耗量可節(jié)省135 kg。

    圖8 動(dòng)力下降推進(jìn)劑總消耗量隨分離高度變化Fig.8 Variation of total propellant consumption with separation height for powered descent

    2.5 基于最佳推重比的月面上升設(shè)計(jì)

    月面動(dòng)力上升的主要目標(biāo)是提高環(huán)月軌道高度和速度至目標(biāo)值,總速度增量達(dá)到1 850 m/s左右,是載人月面著陸器推進(jìn)劑消耗最多的飛行階段之一。該飛行階段受月球引力影響,推重比對(duì)速度增量的影響顯著,如圖9所示。

    圖9 上升段等效速度增量隨推重比變化情況Fig.9 Variation of equivalent velocity increment of powered ascent with the ratio of thrust to weight

    圖10 拋分質(zhì)量隨初始推重比的變化關(guān)系Fig.10 Variation of desert weight with the initial ratio of thrust to weight

    從圖中可以看出,推重比0.25與推重比0.5相比,需要多消耗速度增量約100 m/s。因此,在載人月面著陸器上升段的最佳推重比選擇為0.4~0.6。

    然而,載人月面著陸器為確保安全,需要在特定情況下采取動(dòng)力重構(gòu)的方式繼續(xù)上升[12],這就會(huì)使上升過(guò)程中推重比顯著減小,推進(jìn)劑消耗增大。因此提出在推力不足的情況下,通過(guò)減少起飛重量,以維持較好的推重比,達(dá)到節(jié)約推進(jìn)劑的目的。

    不同起飛推重比下的拋分質(zhì)量mED可由下式估計(jì):

    (9)

    式中:mPA0為月面起飛初始質(zhì)量;ΔmE-N為動(dòng)力重構(gòu)上升推進(jìn)劑使用量較正常飛行的增量設(shè)計(jì)目標(biāo);ΔvPAN和ΔvPAE分別為標(biāo)稱(chēng)和動(dòng)力重構(gòu)時(shí)動(dòng)力上升等效速度增量,與初始推重比ε相關(guān)。

    在初始推重比ε取0.5時(shí),若拋分質(zhì)量mED為297 kg,則動(dòng)力重構(gòu)后的上升推進(jìn)劑消耗與標(biāo)稱(chēng)模式一致。換言之,采用該方法能夠抵消動(dòng)力重構(gòu)造成的推重比損失,與不拋分方案相比,節(jié)約推進(jìn)劑約285 kg。

    2.6 基于快速自主的月軌交會(huì)對(duì)接設(shè)計(jì)

    載人月面著陸器返回環(huán)月軌道的質(zhì)量規(guī)模與交會(huì)對(duì)接的資源需求密切相關(guān),應(yīng)考慮快速交會(huì)對(duì)接方案,通過(guò)大幅壓縮飛行時(shí)長(zhǎng),最大限度地降低載人月面著陸器對(duì)能源、乘員消耗品等的需求。

    嫦娥探月任務(wù)中,環(huán)月交會(huì)對(duì)接模式飛行時(shí)長(zhǎng)約2天[13],快速交會(huì)對(duì)接模式可在環(huán)月軌道1~2圈內(nèi)完成,飛行時(shí)長(zhǎng)壓縮至4%~8%。

    基于霍曼變軌的環(huán)月快速交會(huì)對(duì)接方案見(jiàn)圖11。載人月面著陸器從月面上升后進(jìn)入環(huán)月橢圓軌道,入軌點(diǎn)在近月點(diǎn)附近,飛行半圈后在特征點(diǎn)(入軌軌道遠(yuǎn)月點(diǎn))執(zhí)行Δv1進(jìn)入調(diào)相軌道,在完成調(diào)相后的下一個(gè)特征點(diǎn)執(zhí)行霍曼轉(zhuǎn)移,在半圈內(nèi)通過(guò)兩次變軌Δv3, Δv4進(jìn)入交會(huì)目標(biāo)軌道,之后與載人飛船進(jìn)行近距離控制飛行,一般可再細(xì)分為尋的段、接近段、平移靠攏段等。實(shí)際任務(wù)過(guò)程中,根據(jù)入軌精度按需安排一次軌道面修正,如圖11的Δv2。

    圖11 基于霍曼變軌的環(huán)月快速交會(huì)對(duì)接方案Fig.11 Lunar rapid rendezvous strategy based on the Hohmann maneuver method

    霍曼交會(huì)策略的推進(jìn)劑消耗可用如下公式評(píng)估:

    (10)

    式中:ΔvR為總速度增量;i=1,3,4為各次霍曼轉(zhuǎn)移標(biāo)志;μM為月球引力常數(shù);k和l分別為霍曼轉(zhuǎn)移初始和末端的特征點(diǎn)標(biāo)志,若k代表近月點(diǎn),則l代表遠(yuǎn)月點(diǎn),反之亦然。

    采用快速自主交會(huì)對(duì)接能夠減少消耗品及相關(guān)設(shè)備攜帶量,如水、氣、食物以及廢物收集設(shè)備等,可為載人月面著陸器節(jié)省大量的質(zhì)量和空間。按照阿波羅任務(wù)2名乘員4小時(shí)快速交會(huì)對(duì)接設(shè)計(jì),可節(jié)省干重約43 kg,折合發(fā)射質(zhì)量超過(guò)300 kg。

    3 月面著陸器應(yīng)急飛行策略?xún)?yōu)化

    載人月面著陸器對(duì)航天員安全的要求高,需要在飛行全過(guò)程設(shè)計(jì)應(yīng)急飛行方案,確保在發(fā)生任何故障的情況下,均能將航天員送回在環(huán)月軌道等待的載人飛船。本文關(guān)注的是,通過(guò)飛行策略?xún)?yōu)化實(shí)現(xiàn)整器推進(jìn)劑使用量的合理分配,主要優(yōu)化思路包括:處置故障的推進(jìn)劑消耗量盡量與正常分配值相當(dāng),減少額外推進(jìn)劑配置;壓縮故障處置時(shí)間,減少應(yīng)對(duì)故障額外配置的專(zhuān)門(mén)設(shè)備和資源;綜合評(píng)估利用多重故障下的推進(jìn)劑,減少資源重復(fù)預(yù)留等,以下為幾個(gè)典型示例。

    3.1 基于物資拋分的異面上升設(shè)計(jì)

    環(huán)月軌道傾角由登月點(diǎn)緯度和月面工作時(shí)間共同確定。如圖12所示,根據(jù)軌道和登月點(diǎn)緯度幾何關(guān)系,可得:

    圖12 環(huán)月軌道傾角和著陸點(diǎn)緯度關(guān)系Fig.12 Relationship between the lunar orbit inclination and latitude of landing site

    (11)

    式中:Ds為月面工作時(shí)間;φ為登月點(diǎn)緯度;i為環(huán)月軌道傾角。

    應(yīng)急情況下,月面著陸器需要在圖12中A~B點(diǎn)(著陸點(diǎn)、起飛點(diǎn))之間緊急起飛,登月艙上升入軌后與載人飛船軌道不共面,需要消耗更多的推進(jìn)劑完成后續(xù)交會(huì)對(duì)接任務(wù)。以1°的軌道面偏差修正為例,需要30 m/s的速度增量,約42 kg推進(jìn)劑。在正常情況下,這部分為修正軌道面多攜帶的推進(jìn)劑將被帶回環(huán)月軌道,無(wú)法產(chǎn)生效益,造成資源浪費(fèi)。

    為節(jié)省推進(jìn)劑資源,盡量縮小應(yīng)急飛行與正常飛行的推進(jìn)劑消耗量差額,可采用物資拋分設(shè)計(jì),即在應(yīng)急上升前將部分可舍棄的設(shè)備和物資拋離,通過(guò)降低上升初始質(zhì)量減少應(yīng)急飛行推進(jìn)劑消耗。應(yīng)急上升后的交會(huì)對(duì)接推進(jìn)劑使用量較正常飛行的增量ΔmRD+與應(yīng)急拋分質(zhì)量mED的關(guān)系如下式所示:

    (12)

    式中:ΔmE-N為應(yīng)急飛行推進(jìn)劑使用量較正常飛行的增量設(shè)計(jì)目標(biāo),即通過(guò)應(yīng)急拋分設(shè)計(jì)約束應(yīng)急飛行的推進(jìn)劑使用量;ΔvPA為動(dòng)力上升的等效速度增量。

    圖13給出了應(yīng)急拋分質(zhì)量與軌道面偏差修正設(shè)計(jì)需求。若故障發(fā)生后,通過(guò)在月面拋分170 kg左右的物資,軌道面偏差修正能力可提升至2°,在不額外增加推進(jìn)劑的情況下,滿(mǎn)足大部分低緯度地區(qū)的應(yīng)急異面上升需求。

    圖13 應(yīng)急拋分質(zhì)量與應(yīng)急交會(huì)對(duì)接能力關(guān)系Fig.13 Variations of emergent desert mass with emergent rendezvous ability

    3.2 基于蘭伯特變軌的應(yīng)急交會(huì)策略

    航天員在月面活動(dòng)期間,一旦發(fā)生能源系統(tǒng)癱瘓、密封艙被微流星擊穿等緊急重大故障,需要載人月面著陸器立即起飛返回環(huán)月軌道,盡可能快地將航天員安全送回載人飛船。在2.6節(jié)采用霍曼交會(huì)策略實(shí)現(xiàn)2圈月軌快速交會(huì)的基礎(chǔ)上,進(jìn)一步考慮能源、壓力應(yīng)急等緊急故障下的月軌交會(huì)問(wèn)題,可以采用更加快速的蘭伯特交會(huì)策略,降低蓄電池容量、應(yīng)急氣體等配置需求。

    基于蘭伯特變軌的環(huán)月快速交會(huì)對(duì)接方案如圖14所示。載人月面著陸器入軌后在合適的相位執(zhí)行蘭伯特變軌的首脈沖ΔvLbt1,進(jìn)入橢圓轉(zhuǎn)移軌道,完成相位追及后執(zhí)行蘭伯特變軌的末脈沖ΔvLbt2,進(jìn)入目標(biāo)交會(huì)軌道。根據(jù)實(shí)際在軌執(zhí)行精度,ΔvLbt1和ΔvLbt2之間可視情安排幾次修正脈沖。若快速性要求有所降低,蘭伯特變軌策略也可結(jié)合特征點(diǎn)變軌執(zhí)行,即載人月面著陸器入軌后先通過(guò)特征點(diǎn)變軌進(jìn)入一個(gè)過(guò)渡軌道再執(zhí)行蘭伯特交會(huì),如阿波羅計(jì)劃的共橢圓交會(huì)策略[11]。

    圖14 基于蘭伯特變軌的環(huán)月快速交會(huì)對(duì)接方案Fig.14 Lunar rapid rendezvous strategy based on the Lambert maneuver method

    蘭伯特變軌的速度增量需求與交會(huì)對(duì)接時(shí)長(zhǎng)要求相關(guān),一般來(lái)說(shuō)初始相位偏差越大、交會(huì)時(shí)長(zhǎng)越短,速度增量需求越大,如下式所示:

    (13)

    式中:ΔtLbt為蘭伯特轉(zhuǎn)移段時(shí)長(zhǎng);aLbt為蘭伯特轉(zhuǎn)移軌道半長(zhǎng)軸;r1和r2分別為蘭伯特轉(zhuǎn)移的初始和末端月心距;c為初始和末端的月心矢量間弦長(zhǎng);v1和v2分別為蘭伯特轉(zhuǎn)移的初始和目標(biāo)速度;vL1和vL2分別為蘭伯特轉(zhuǎn)移軌道的初始和末端速度;φL1和φL2分別為初始和末端的速度矢量夾角。

    圖15給出了蘭伯特交會(huì)策略和霍曼交會(huì)策略在速度增量和交會(huì)時(shí)長(zhǎng)的對(duì)比算例。算例中,在6°的初始相位差下,蘭伯特交會(huì)策略通過(guò)增加63 m/s的速度增量可將交會(huì)時(shí)長(zhǎng)縮短93 min,為應(yīng)對(duì)該緊急模式的蓄電池、應(yīng)急氣瓶等配置可以減配,干重減少約30 kg,折合發(fā)射質(zhì)量超過(guò)250 kg。

    圖15 蘭伯特交會(huì)策略和霍曼交會(huì)策略對(duì)比Fig.15 Comparison between the Lambert and Hohmann rendezvous strategies

    3.3 基于PRA的多重故障應(yīng)對(duì)策略

    載人月面著陸器在飛行全任務(wù)階段具有多種故障模式。按照傳統(tǒng)的推進(jìn)劑分配方式,需要對(duì)每個(gè)飛行階段的故障模式進(jìn)行獨(dú)立分析,逐一獲取推進(jìn)劑消耗量,再將每個(gè)飛行階段中處置獨(dú)立故障的推進(jìn)劑消耗最大值按照飛行過(guò)程逐一加總,作為整個(gè)飛行器全任務(wù)階段處置故障的預(yù)留推進(jìn)劑。這種方法無(wú)疑會(huì)將不同飛行階段相關(guān)故障模式進(jìn)行重復(fù)累計(jì),造成推進(jìn)劑預(yù)留量過(guò)于保守。

    提出基于概率風(fēng)險(xiǎn)評(píng)估(PRA)的多重故障推進(jìn)劑預(yù)算方法,從飛行器的任務(wù)剖面出發(fā),將事件樹(shù)和故障樹(shù)相結(jié)合,對(duì)正常和故障飛行場(chǎng)景進(jìn)行建模,以準(zhǔn)確識(shí)別可能出現(xiàn)的多重故障應(yīng)急模式分支。

    圖16給出了基于事件樹(shù)的載人月面著陸器多重故障應(yīng)急模式分析示例。以“N_”標(biāo)注的是正常事件,以“E_”標(biāo)注的是應(yīng)急事件。在示例的事件樹(shù)中,共1個(gè)正常飛行模式和11個(gè)應(yīng)急飛行分支,其中包含1種應(yīng)急事件的分支有5條,包含2種應(yīng)急事件的分支有5條,包含3種應(yīng)急事件的分支有1條,不存在4種及以上應(yīng)急事件組合的分支。

    圖16 載人月面著陸器的事件樹(shù)建模示例Fig.16 An example of event tree modeling for a manned lunar lander

    以正常飛行模式的推進(jìn)劑使用量為基線(xiàn),若按傳統(tǒng)方法將處置應(yīng)急事件的推進(jìn)劑累加,需額外預(yù)留308 kg推進(jìn)劑。按照事件樹(shù)分析方法,11個(gè)應(yīng)急飛行分支中需要推進(jìn)劑最多的是包含3種應(yīng)急事件的分支:“入軌偏差大+月面提前起飛+動(dòng)力重構(gòu)上升”分支,處置故障需額外消耗推進(jìn)劑270 kg,比傳統(tǒng)方法有所減少。

    進(jìn)一步采用PRA方法對(duì)應(yīng)急飛行分支發(fā)生概率進(jìn)行定量分析,若“入軌偏差大+月面提前起飛+動(dòng)力重構(gòu)”和“入軌偏差大+動(dòng)力重構(gòu)”2條應(yīng)急分支的多重故障無(wú)相關(guān)性且發(fā)生概率極低,則可以作為系統(tǒng)殘余風(fēng)險(xiǎn),不額外預(yù)留推進(jìn)劑進(jìn)行處置。因此,在該示例中可以考慮將剩余的9條應(yīng)急分支的最大推進(jìn)劑消耗量作為系統(tǒng)處置故障的預(yù)留推進(jìn)劑,為190 kg,僅相當(dāng)于傳統(tǒng)方法的62%。

    4 結(jié) 論

    圍繞載人月面著陸器推進(jìn)劑使用量控制這一核心目標(biāo),在正常、應(yīng)急飛行策略中提出了一系列總體層次的優(yōu)化思路,可為月面著陸器飛行任務(wù)規(guī)劃和軌道設(shè)計(jì)提供借鑒和參考。在保證航天員安全的前提下,盡量實(shí)現(xiàn)月面運(yùn)輸能力的最大化,支持載人月球探測(cè)工程取得更多工程和科學(xué)成果。

    猜你喜歡
    環(huán)月著陸器交會(huì)
    “織交會(huì)”打造大朗樣板
    嫦娥五號(hào)探測(cè)器環(huán)月運(yùn)行
    上海航天(2020年6期)2021-01-18 04:50:12
    環(huán)月降軌實(shí)現(xiàn)月面著陸的控制策略
    嫦娥四號(hào)巡視器拍攝的著陸器圖像
    航天器工程(2019年4期)2019-11-11 03:41:00
    2019年,水交會(huì)來(lái)了!
    中國(guó)完成首次火星探測(cè)任務(wù)著陸器懸停避障試驗(yàn)
    伙伴(2019年12期)2019-01-17 04:38:56
    嫦娥四號(hào)順利進(jìn)入環(huán)月軌道
    基于多工況的新型著陸器軟著陸性能優(yōu)化
    環(huán)月遠(yuǎn)程快速交會(huì)任務(wù)規(guī)劃
    載人航天(2016年4期)2016-12-01 06:56:13
    立方體星交會(huì)對(duì)接和空間飛行演示
    太空探索(2016年9期)2016-07-12 09:59:53
    宅男免费午夜| 婷婷精品国产亚洲av在线| 午夜影院日韩av| 在线永久观看黄色视频| 香蕉丝袜av| 香蕉国产在线看| 国产单亲对白刺激| 亚洲一区二区三区不卡视频| 精品人妻1区二区| 亚洲国产看品久久| 另类亚洲欧美激情| 丁香欧美五月| 黄色 视频免费看| 亚洲午夜精品一区,二区,三区| 久久精品亚洲熟妇少妇任你| 久久久久国产精品人妻aⅴ院| 亚洲av五月六月丁香网| 手机成人av网站| 欧美日韩av久久| 黑人巨大精品欧美一区二区mp4| 国产亚洲av高清不卡| 日本黄色视频三级网站网址| 欧美成人午夜精品| 午夜精品国产一区二区电影| 亚洲精品在线观看二区| 欧美激情 高清一区二区三区| 天堂中文最新版在线下载| 在线观看免费午夜福利视频| xxxhd国产人妻xxx| 亚洲精品中文字幕在线视频| 欧美黑人欧美精品刺激| 中文欧美无线码| 可以免费在线观看a视频的电影网站| 久久 成人 亚洲| 午夜两性在线视频| 欧美午夜高清在线| 电影成人av| 性欧美人与动物交配| 啦啦啦在线免费观看视频4| 一边摸一边做爽爽视频免费| 国产精品国产av在线观看| 免费女性裸体啪啪无遮挡网站| 日韩有码中文字幕| 日日摸夜夜添夜夜添小说| 一级a爱视频在线免费观看| 波多野结衣av一区二区av| av有码第一页| 精品国产一区二区久久| 中文字幕高清在线视频| 999久久久国产精品视频| 天堂动漫精品| 69精品国产乱码久久久| 丁香六月欧美| 欧美+亚洲+日韩+国产| 国产精品二区激情视频| 后天国语完整版免费观看| 色尼玛亚洲综合影院| 亚洲va日本ⅴa欧美va伊人久久| 97碰自拍视频| 视频区欧美日本亚洲| 午夜精品久久久久久毛片777| 天堂影院成人在线观看| 国产亚洲精品久久久久久毛片| 男人舔女人下体高潮全视频| 久久亚洲真实| 久久久国产一区二区| 又黄又粗又硬又大视频| 国产精品av久久久久免费| 99热国产这里只有精品6| 自拍欧美九色日韩亚洲蝌蚪91| 多毛熟女@视频| 日韩高清综合在线| 国产精华一区二区三区| 最新在线观看一区二区三区| 性欧美人与动物交配| 一级作爱视频免费观看| 黄色毛片三级朝国网站| 最好的美女福利视频网| 国产精品av久久久久免费| av片东京热男人的天堂| 国产精品久久久久成人av| av天堂在线播放| 美女高潮到喷水免费观看| tocl精华| svipshipincom国产片| 免费日韩欧美在线观看| 午夜影院日韩av| 国产av精品麻豆| 99久久精品国产亚洲精品| 一级毛片精品| 久久人人精品亚洲av| 久久久久久大精品| 久久久久国产一级毛片高清牌| www日本在线高清视频| 国产亚洲精品第一综合不卡| 日本a在线网址| 可以免费在线观看a视频的电影网站| 欧美日韩国产mv在线观看视频| 成人国产一区最新在线观看| 亚洲九九香蕉| 美女扒开内裤让男人捅视频| 精品少妇一区二区三区视频日本电影| 日本精品一区二区三区蜜桃| 亚洲色图 男人天堂 中文字幕| 最新美女视频免费是黄的| 亚洲欧美精品综合一区二区三区| 国产精品二区激情视频| 国产国语露脸激情在线看| 99国产精品99久久久久| 好男人电影高清在线观看| 国产一区二区激情短视频| 国产精品成人在线| 亚洲av片天天在线观看| 免费在线观看视频国产中文字幕亚洲| 久久 成人 亚洲| 成人18禁高潮啪啪吃奶动态图| 国产片内射在线| 久久中文看片网| 黄片小视频在线播放| 午夜福利免费观看在线| 十分钟在线观看高清视频www| 欧美乱码精品一区二区三区| 欧美在线黄色| 夜夜爽天天搞| 欧美大码av| 亚洲精品成人av观看孕妇| 丰满的人妻完整版| 国产亚洲精品久久久久5区| 麻豆av在线久日| 亚洲午夜精品一区,二区,三区| 亚洲av成人一区二区三| 久久久久久人人人人人| 久久中文字幕人妻熟女| 两人在一起打扑克的视频| 亚洲精品一二三| 精品国产美女av久久久久小说| 嫩草影院精品99| 亚洲五月天丁香| 国产单亲对白刺激| 老司机福利观看| 欧美丝袜亚洲另类 | 最新美女视频免费是黄的| 欧美乱色亚洲激情| 丝袜在线中文字幕| av欧美777| 中文字幕另类日韩欧美亚洲嫩草| 国产成人影院久久av| av网站在线播放免费| 一区二区三区国产精品乱码| 国产亚洲精品久久久久久毛片| 巨乳人妻的诱惑在线观看| 欧美在线一区亚洲| 亚洲色图综合在线观看| 国产av在哪里看| 少妇裸体淫交视频免费看高清 | 午夜视频精品福利| 精品国产乱码久久久久久男人| 亚洲avbb在线观看| 亚洲成人久久性| 中文字幕人妻丝袜一区二区| 免费观看人在逋| 免费看十八禁软件| 久久草成人影院| 啦啦啦在线免费观看视频4| 国产伦一二天堂av在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品aⅴ一区二区三区四区| 国产精品偷伦视频观看了| 99精品在免费线老司机午夜| 久久人人精品亚洲av| 亚洲五月天丁香| 看免费av毛片| 一区福利在线观看| 在线观看免费午夜福利视频| 又黄又粗又硬又大视频| 亚洲精品美女久久久久99蜜臀| 欧美日韩福利视频一区二区| 在线十欧美十亚洲十日本专区| 男女做爰动态图高潮gif福利片 | www国产在线视频色| 欧美激情久久久久久爽电影 | 母亲3免费完整高清在线观看| 正在播放国产对白刺激| 少妇粗大呻吟视频| 国产成人精品无人区| 制服人妻中文乱码| 精品人妻1区二区| 18美女黄网站色大片免费观看| 黄色a级毛片大全视频| 成人国产一区最新在线观看| 亚洲五月天丁香| 国产精品电影一区二区三区| 亚洲自拍偷在线| 最近最新免费中文字幕在线| 亚洲av片天天在线观看| 国产精品国产av在线观看| 午夜精品久久久久久毛片777| 精品久久久精品久久久| 在线免费观看的www视频| 国产精品久久久av美女十八| 女同久久另类99精品国产91| 香蕉久久夜色| 日韩中文字幕欧美一区二区| 亚洲人成网站在线播放欧美日韩| a级毛片黄视频| 嫩草影视91久久| 长腿黑丝高跟| 十八禁人妻一区二区| 亚洲欧美激情综合另类| 成人18禁在线播放| 欧美 亚洲 国产 日韩一| 伊人久久大香线蕉亚洲五| 日韩免费av在线播放| 亚洲激情在线av| 视频区欧美日本亚洲| 欧美日韩福利视频一区二区| 我的亚洲天堂| 国产精品1区2区在线观看.| 国产亚洲精品第一综合不卡| 88av欧美| 欧美激情高清一区二区三区| 别揉我奶头~嗯~啊~动态视频| 操美女的视频在线观看| 一级毛片精品| 国产视频一区二区在线看| 亚洲少妇的诱惑av| 99热只有精品国产| 99国产精品一区二区三区| 欧美日韩黄片免| 午夜免费观看网址| 亚洲国产精品一区二区三区在线| 日本一区二区免费在线视频| 成人三级做爰电影| 亚洲色图综合在线观看| 久久精品91蜜桃| 日韩欧美三级三区| 大码成人一级视频| 老熟妇乱子伦视频在线观看| 日韩中文字幕欧美一区二区| 一进一出抽搐动态| 亚洲精品国产色婷婷电影| 国产精品98久久久久久宅男小说| 国产不卡一卡二| 亚洲精品国产精品久久久不卡| 久久精品成人免费网站| 91麻豆av在线| 日日夜夜操网爽| 999精品在线视频| 国产一卡二卡三卡精品| 日韩av在线大香蕉| 免费不卡黄色视频| 99精品在免费线老司机午夜| 妹子高潮喷水视频| 男人舔女人下体高潮全视频| 免费观看人在逋| 欧美丝袜亚洲另类 | 热99re8久久精品国产| 免费高清视频大片| 人成视频在线观看免费观看| 色婷婷av一区二区三区视频| 在线国产一区二区在线| 天天躁夜夜躁狠狠躁躁| 男女床上黄色一级片免费看| 久久午夜亚洲精品久久| 国内久久婷婷六月综合欲色啪| 国产亚洲精品一区二区www| 久久久久久免费高清国产稀缺| 欧美黑人精品巨大| 最近最新中文字幕大全免费视频| 亚洲成a人片在线一区二区| 久久影院123| 成人永久免费在线观看视频| 免费在线观看影片大全网站| 长腿黑丝高跟| 久久午夜亚洲精品久久| 50天的宝宝边吃奶边哭怎么回事| 曰老女人黄片| 午夜精品国产一区二区电影| 日本黄色日本黄色录像| 午夜精品在线福利| 免费在线观看日本一区| 777久久人妻少妇嫩草av网站| 97碰自拍视频| 午夜福利一区二区在线看| 交换朋友夫妻互换小说| 中文字幕人妻熟女乱码| 国产精品一区二区在线不卡| 久久国产乱子伦精品免费另类| 超色免费av| 超碰成人久久| 黑丝袜美女国产一区| 黑人猛操日本美女一级片| 成人国产一区最新在线观看| 久久久久久久久久久久大奶| 老司机深夜福利视频在线观看| 亚洲男人天堂网一区| 中文字幕色久视频| 神马国产精品三级电影在线观看 | 不卡av一区二区三区| 十八禁网站免费在线| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 日韩高清综合在线| 欧美日韩瑟瑟在线播放| 99国产综合亚洲精品| 中文字幕最新亚洲高清| 在线视频色国产色| 老司机午夜十八禁免费视频| 黄色 视频免费看| 亚洲欧美一区二区三区黑人| 女人爽到高潮嗷嗷叫在线视频| 在线免费观看的www视频| 色综合欧美亚洲国产小说| 三上悠亚av全集在线观看| 亚洲精华国产精华精| 午夜日韩欧美国产| 老汉色av国产亚洲站长工具| 伦理电影免费视频| 国产高清激情床上av| 免费不卡黄色视频| 黑人巨大精品欧美一区二区蜜桃| 国产极品粉嫩免费观看在线| 麻豆久久精品国产亚洲av | 日本免费一区二区三区高清不卡 | 国产精品一区二区在线不卡| 午夜a级毛片| 女人精品久久久久毛片| 国产男靠女视频免费网站| 中文欧美无线码| 国产精品99久久99久久久不卡| 丝袜美腿诱惑在线| 黑人操中国人逼视频| 亚洲精品美女久久av网站| 精品久久久久久久久久免费视频 | 一进一出抽搐动态| 日本免费一区二区三区高清不卡 | 精品一区二区三区av网在线观看| 男女做爰动态图高潮gif福利片 | 看片在线看免费视频| 男女床上黄色一级片免费看| 精品一区二区三区av网在线观看| 99久久人妻综合| 久久久国产精品麻豆| 18禁黄网站禁片午夜丰满| 两性午夜刺激爽爽歪歪视频在线观看 | 国产精品久久久人人做人人爽| 精品无人区乱码1区二区| 亚洲一区二区三区色噜噜 | 一本综合久久免费| 国产一区二区激情短视频| 精品第一国产精品| 久久中文看片网| 十八禁人妻一区二区| 可以在线观看毛片的网站| 一二三四社区在线视频社区8| 午夜成年电影在线免费观看| 精品国产一区二区三区四区第35| 国产成人av教育| 精品国产超薄肉色丝袜足j| www.www免费av| 亚洲色图av天堂| 叶爱在线成人免费视频播放| 欧美日韩国产mv在线观看视频| 久久精品人人爽人人爽视色| 99精品在免费线老司机午夜| 亚洲中文av在线| 国产精品电影一区二区三区| 琪琪午夜伦伦电影理论片6080| 国产精品电影一区二区三区| 黄色视频不卡| 欧洲精品卡2卡3卡4卡5卡区| 欧美日韩av久久| 亚洲欧洲精品一区二区精品久久久| 最新美女视频免费是黄的| 999久久久精品免费观看国产| 黑丝袜美女国产一区| 久久久久久久久中文| 五月开心婷婷网| 男男h啪啪无遮挡| 欧美日韩中文字幕国产精品一区二区三区 | 青草久久国产| 日本黄色视频三级网站网址| 欧美日韩一级在线毛片| 日韩欧美国产一区二区入口| 久久久精品国产亚洲av高清涩受| 久久欧美精品欧美久久欧美| 老司机午夜十八禁免费视频| 久久久国产欧美日韩av| www国产在线视频色| 亚洲国产毛片av蜜桃av| 十分钟在线观看高清视频www| 国产野战对白在线观看| 日本vs欧美在线观看视频| 久久 成人 亚洲| 啪啪无遮挡十八禁网站| 日韩一卡2卡3卡4卡2021年| 欧美乱码精品一区二区三区| 精品一区二区三卡| 欧美在线一区亚洲| 91大片在线观看| 国产男靠女视频免费网站| 免费人成视频x8x8入口观看| 黄色a级毛片大全视频| 99国产极品粉嫩在线观看| 中国美女看黄片| 怎么达到女性高潮| 亚洲第一av免费看| 麻豆国产av国片精品| 国产成+人综合+亚洲专区| 另类亚洲欧美激情| 欧美日韩一级在线毛片| 一本综合久久免费| 亚洲激情在线av| 99久久国产精品久久久| 女性被躁到高潮视频| 啦啦啦 在线观看视频| 欧洲精品卡2卡3卡4卡5卡区| 久久久久久免费高清国产稀缺| 国产亚洲欧美在线一区二区| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 久久亚洲真实| 99国产精品99久久久久| 他把我摸到了高潮在线观看| 两性夫妻黄色片| 国产精品国产高清国产av| 久久中文看片网| 亚洲第一欧美日韩一区二区三区| 天堂俺去俺来也www色官网| 极品教师在线免费播放| 国产精品永久免费网站| 自线自在国产av| 日本a在线网址| 亚洲午夜理论影院| 亚洲av第一区精品v没综合| 亚洲午夜精品一区,二区,三区| av超薄肉色丝袜交足视频| 国产精品一区二区精品视频观看| 最近最新免费中文字幕在线| 亚洲三区欧美一区| 国产欧美日韩综合在线一区二区| www.自偷自拍.com| 免费在线观看亚洲国产| 亚洲国产精品sss在线观看 | 午夜91福利影院| 国内毛片毛片毛片毛片毛片| 无限看片的www在线观看| 三上悠亚av全集在线观看| 69精品国产乱码久久久| 亚洲自拍偷在线| 亚洲国产中文字幕在线视频| 黄片小视频在线播放| 中文字幕av电影在线播放| 久久中文字幕一级| 一级黄色大片毛片| 真人一进一出gif抽搐免费| 乱人伦中国视频| 欧美成人免费av一区二区三区| 日本 av在线| 国产成人av激情在线播放| 琪琪午夜伦伦电影理论片6080| 欧美乱色亚洲激情| 99riav亚洲国产免费| 亚洲国产欧美日韩在线播放| e午夜精品久久久久久久| 啦啦啦免费观看视频1| 91大片在线观看| 国产精品日韩av在线免费观看 | 国产xxxxx性猛交| 成人影院久久| 国产免费现黄频在线看| 手机成人av网站| 久久久久九九精品影院| 久久久久亚洲av毛片大全| 欧美日韩亚洲国产一区二区在线观看| 亚洲成av片中文字幕在线观看| 99久久99久久久精品蜜桃| 亚洲欧美一区二区三区黑人| 在线观看www视频免费| x7x7x7水蜜桃| 久久精品国产亚洲av香蕉五月| 9色porny在线观看| 俄罗斯特黄特色一大片| 侵犯人妻中文字幕一二三四区| av福利片在线| 国产一区在线观看成人免费| 成人亚洲精品av一区二区 | 人人妻人人澡人人看| 国产成人欧美在线观看| 欧美日韩亚洲高清精品| 欧美日韩视频精品一区| 久热爱精品视频在线9| 午夜亚洲福利在线播放| 欧美激情久久久久久爽电影 | 51午夜福利影视在线观看| 国产黄色免费在线视频| 搡老乐熟女国产| 国产亚洲精品一区二区www| av天堂久久9| 亚洲国产精品合色在线| 国产一区二区在线av高清观看| 91麻豆av在线| 一区二区三区激情视频| 国产精品爽爽va在线观看网站 | 在线观看免费视频日本深夜| 1024香蕉在线观看| 在线观看日韩欧美| 亚洲国产精品合色在线| 久久精品国产清高在天天线| 亚洲欧美激情综合另类| 成人特级黄色片久久久久久久| 国产成人精品在线电影| 夫妻午夜视频| 欧美日韩亚洲高清精品| 日韩av在线大香蕉| 欧美大码av| 激情在线观看视频在线高清| 狠狠狠狠99中文字幕| 久久热在线av| 色哟哟哟哟哟哟| 女生性感内裤真人,穿戴方法视频| 精品第一国产精品| 自线自在国产av| 青草久久国产| 欧美日韩国产mv在线观看视频| 亚洲熟妇中文字幕五十中出 | 老司机深夜福利视频在线观看| 亚洲成av片中文字幕在线观看| 亚洲人成电影观看| 国产99久久九九免费精品| 国产激情久久老熟女| 妹子高潮喷水视频| 在线观看免费视频日本深夜| 免费少妇av软件| 中文字幕另类日韩欧美亚洲嫩草| 青草久久国产| 国产成人免费无遮挡视频| 亚洲 欧美 日韩 在线 免费| 在线观看www视频免费| 成人亚洲精品一区在线观看| e午夜精品久久久久久久| 国产免费男女视频| 日本免费一区二区三区高清不卡 | 欧美一级毛片孕妇| 好看av亚洲va欧美ⅴa在| 老熟妇仑乱视频hdxx| 女性生殖器流出的白浆| 男人操女人黄网站| 精品久久蜜臀av无| 色尼玛亚洲综合影院| 九色亚洲精品在线播放| 国产精品综合久久久久久久免费 | 亚洲色图 男人天堂 中文字幕| 国产在线精品亚洲第一网站| 搡老岳熟女国产| 亚洲五月婷婷丁香| aaaaa片日本免费| 不卡av一区二区三区| 欧美另类亚洲清纯唯美| 久久午夜综合久久蜜桃| 亚洲va日本ⅴa欧美va伊人久久| 女性被躁到高潮视频| 国产又色又爽无遮挡免费看| 中文字幕人妻丝袜制服| 成年女人毛片免费观看观看9| 丁香欧美五月| 欧美亚洲日本最大视频资源| 18禁美女被吸乳视频| 国产人伦9x9x在线观看| 搡老岳熟女国产| 国产精品美女特级片免费视频播放器 | 国产一区二区三区在线臀色熟女 | 国产片内射在线| 一二三四在线观看免费中文在| 久久狼人影院| 夜夜爽天天搞| aaaaa片日本免费| 高清欧美精品videossex| 久久国产精品影院| 免费少妇av软件| www.999成人在线观看| 久久亚洲真实| 怎么达到女性高潮| 国产精品影院久久| 在线天堂中文资源库| 18禁黄网站禁片午夜丰满| 侵犯人妻中文字幕一二三四区| av在线播放免费不卡| 亚洲国产精品一区二区三区在线| 精品国产国语对白av| 久久狼人影院| 成年人黄色毛片网站| 香蕉久久夜色| 国产精品免费一区二区三区在线| 视频在线观看一区二区三区| 欧美大码av| 超碰97精品在线观看| 久久久精品欧美日韩精品| 久久久久国产一级毛片高清牌| 精品少妇一区二区三区视频日本电影| 国产91精品成人一区二区三区| 免费在线观看亚洲国产| 国产精品国产高清国产av| 亚洲第一青青草原| 人人妻人人澡人人看| 亚洲五月色婷婷综合| 淫秽高清视频在线观看| 欧美中文综合在线视频| 亚洲国产精品999在线| 国产精品爽爽va在线观看网站 | 天天添夜夜摸| 丰满迷人的少妇在线观看| 久久欧美精品欧美久久欧美| 免费搜索国产男女视频| 怎么达到女性高潮|