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

    基于不同動力引力輔助模型的木星轉(zhuǎn)移軌道設(shè)計

    2019-07-05 02:06:54楊洪偉
    上海航天 2019年3期
    關(guān)鍵詞:木星引力機動

    楊 彬,楊洪偉,李 爽,尤 偉

    (1. 南京航空航天大學(xué) 航天學(xué)院,江蘇 南京 210016; 2. 南京航空航天大學(xué) 航天新技術(shù)實驗室,江蘇 南京 210016; 3. 上海衛(wèi)星工程研究所,上海 201109)

    0 引言

    隨著空間技術(shù)的快速發(fā)展,我國的深空探測能力逐漸提高。以木星為代表的外太陽系天體成為我國下一步深空探測的目標(biāo)[1]。對于這類距離地球遙遠(yuǎn)的目標(biāo),傳統(tǒng)的直接轉(zhuǎn)移方式因燃料消耗巨大而不再適用。引力輔助技術(shù)作為重要的節(jié)能技術(shù),對于深空探測軌道設(shè)計具有重要意義。引力輔助技術(shù)的應(yīng)用能有效減少任務(wù)的燃料消耗,使得原本無法完成的任務(wù)成為可能。而且,除了完成對既定目標(biāo)的探測外,探測器在行星引力輔助過程中還能完成對引力輔助行星的飛越探測,增加任務(wù)的科學(xué)回報。國外包括旅行者號[2]、卡西尼號[3]在內(nèi)的一系列深空探測器均采用了引力輔助技術(shù),極大減少了任務(wù)的燃料消耗。

    引力輔助技術(shù),即在探測器靠近某一星球時,借助該星球的引力作用來改變探測器的軌道與速度,是一種無消耗的軌道機動方式。引力輔助技術(shù)的應(yīng)用使得深空探測軌道設(shè)計變得更加復(fù)雜多樣。為獲得較好的引力輔助效果,針對復(fù)雜引力輔助軌跡的優(yōu)化設(shè)計一直是軌道動力學(xué)的熱點話題。CESARONE[4]推導(dǎo)了引力輔助過程中的參數(shù)表達(dá)式,并分析了近心點高度等引力輔助參數(shù)對引力輔助軌跡的影響。盡管引力輔助技術(shù)在節(jié)省燃料、提高效率等方面具有顯著優(yōu)勢,但其對所涉及天體的相位關(guān)系要求十分嚴(yán)格,借力前后的速度必須嚴(yán)格匹配[5]。這使得引力輔助轉(zhuǎn)移軌跡的發(fā)射窗口變得十分狹窄,造成搜索過程收斂困難。為提高發(fā)射窗口搜索算法的穩(wěn)定性,LONGUSKI等[6]提出了特征能量C3匹配法,基于引力輔助前后探測器的特征能量進行引力輔助轉(zhuǎn)移軌跡匹配。張旭輝等[7]采用蒂塞朗準(zhǔn)則對C3匹配方法進行了驗證。喬棟等[8]提出了軟匹配策略,用以進行軌道拼接。軟匹配策略的應(yīng)用有效解決了傳統(tǒng)方法對滿足約束條件轉(zhuǎn)移軌道方案的遺漏問題。雖然采用了各種方法來提高匹配的穩(wěn)定性,但等待相關(guān)行星到達(dá)理想相位關(guān)系,往往需要數(shù)年時間。而且,針對多天體引力輔助情況,等待時間將更長。

    為弱化引力輔助技術(shù)對能量匹配的約束,一些學(xué)者提出在引力輔助過程中,通過施加主動脈沖來改善引力輔助效果,增加任務(wù)窗口,減少相位等待時間。PRADO[9]對帶有脈沖機動的引力輔助變軌進行了描述,并推導(dǎo)了探測器速度、能量變化量等參數(shù),給出了特定條件下的最優(yōu)脈沖機動方向。引力輔助過程中脈沖施加的具體位置目前仍是軌道動力學(xué)研究的熱點問題[10-12]。李小玉等[13]研究了附加深空機動的引力輔助模型,推導(dǎo)了引力輔助后離開剩余速度的解析表達(dá)式。賈建華等[14]分析了平面引力輔助過程中施加脈沖的最優(yōu)位置,給出了平面情況下的最優(yōu)脈沖施加位置確定方法。CASALINO等[15]研究了在近拱點施加脈沖機動的引力輔助,但其所提方法只適用于脈沖幅值較小的情況。賈建華等[16]分析了三維情況下近拱點施加脈沖的大小和方向?qū)μ綔y器軌道的影響。除了在近拱點施加脈沖機動外,另一種方案是在引力輔助完成后施加脈沖機動。侯艷偉等[17]采用這種方式設(shè)計了基于金星引力輔助的火星轉(zhuǎn)移軌道。

    本文針對有動力引力輔助技術(shù)和木星探測轉(zhuǎn)移軌跡進行了研究,主要對近心點脈沖機動和甩擺后脈沖機動2種引力輔助模型進行了對比分析,并完成了木星引力輔助轉(zhuǎn)移軌跡的初始設(shè)計。首先,分別構(gòu)建近心點機動引力輔助模型和甩擺后機動引力輔助模型。然后,提出基于引力輔助模型進行深空探測的轉(zhuǎn)移軌道設(shè)計方法。隨后,以木星探測任務(wù)為工程應(yīng)用背景,基于不同引力輔助行星方案測試2種動力引力輔助模型的性能。數(shù)值仿真結(jié)果表明:近心點施加脈沖機動的引力輔助方案具有更好的節(jié)能效果。最后,應(yīng)用近心點機動引力輔助模型,采用地球-金星-金星引力輔助序列,完成木星探測任務(wù)轉(zhuǎn)移軌跡的設(shè)計。

    1 動力引力輔助模型

    引力輔助作為一種高效的節(jié)能技術(shù),多次為深空探測任務(wù)所采用。但對于軌道設(shè)計而言,包含引力輔助的轉(zhuǎn)移軌跡頗為復(fù)雜。因此,各種簡化模型被采用,以在保證一定精度的前提下盡量提高計算效率。其中,圓錐曲線拼接法所采用的影響球模型是目前初始設(shè)計最常用的簡化模型。根據(jù)圓錐曲線拼接法的原理,引力輔助過程被假設(shè)為一個瞬時速度突變。其中,速度突變量的大小根據(jù)影響球模型估算求解。無動力引力輔助方式受到安全飛越高度等因素限制,在大多數(shù)情況下不能滿足任務(wù)要求。因此,在引力輔助過程中或引力輔助結(jié)束后通過施加脈沖機動來匹配引力輔助前后的轉(zhuǎn)移軌跡是一種有效方式。這種將主動變軌和引力輔助相結(jié)合的方式被稱為有動力引力輔助技術(shù)。本文將針對近心點機動和引力輔助后機動2種動力引力輔助方式進行建模分析。

    1.1 近心點機動引力輔助模型

    近心點機動引力輔助過程中只考慮引力輔助行星的引力作用,即探測器進入借力行星影響球后沿雙曲線軌道接近借力行星。當(dāng)其飛至近心點時,沿速度方向施加脈沖機動,以修正探測器引力輔助完成后的狀態(tài)。隨后,探測器沿調(diào)整后的逃逸軌道飛出借力行星的影響球,完成引力輔助。

    相比于漫長的星際轉(zhuǎn)移過程,引力輔助過程持續(xù)時間很短。因此,假設(shè)引力輔助瞬時完成,則引力輔助前后探測器的位置不變且與引力輔助行星位置重合。探測器引力輔助前后的速度v+和v-分別通過求解引力輔助前后2段轉(zhuǎn)移軌跡的Lambert問題得到。因此,探測器進入借力行星影響球的剩余速度和理想的離開剩余速度分別為

    (1)

    式中:vP為借力行星在引力輔助時刻的速度。若探測器以理想剩余速度進入和離開引力輔助行星的影響球,則進入和離開影響球的真近點角等于雙曲線軌道漸近線的真近點角。進入和離開的真近點角的數(shù)學(xué)表達(dá)式分別為

    (2)

    根據(jù)雙曲線的性質(zhì),進入和離開真近點角和雙曲線剩余速度偏轉(zhuǎn)角δ滿足幾何約束,f+-f--δ-π=0。將式(2)帶入幾何約束方程,得到關(guān)于近心點高度的非線性方程

    (3)

    為求解非線性方程(3),對其求一階導(dǎo)數(shù)和二階導(dǎo)數(shù),即

    (4)

    (5)

    (6)

    1.2 甩擺后機動引力輔助模型

    動力引力輔助除在近心點添加機動外,另一種常見的方案是在引力輔助完成后立即施加速度脈沖進行調(diào)整,使得探測器的狀態(tài)滿足后續(xù)飛行任務(wù)的要求。施加的速度增量滿足

    v+=v-+ΔvGA+ΔvI

    (7)

    式中:ΔvGA為引力輔助行星提供的速度增量;ΔvI為引力輔助完成后探測器施加的速度增量。在式(7)中,當(dāng)任務(wù)時間確定后,引力輔助前后探測器的速度是固定的。因此,為節(jié)省燃料,只能調(diào)節(jié)引力輔助提供的速度增量,使得探測器提供的速度增量最小。

    為方便描述,本文采用引力輔助坐標(biāo)系[12]描述引力輔助過程的參數(shù),如圖1所示。

    圖1 引力輔助坐標(biāo)系Fig.1 Gravity assist reference frame

    圖中:坐標(biāo)原點O位于引力輔助行星的中心,ξ軸沿雙曲線進入剩余速度方向,ζ軸垂直于引力輔助前探測器的速度與引力輔助行星速度所決定的平面,η軸與ζ軸和ξ軸構(gòu)成右手坐標(biāo)系。將坐標(biāo)軸方向的單位矢量分別記為i,j和k,用剩余速度矢量和行星狀態(tài)表示,如式(8)所示,具體推導(dǎo)過程詳見文獻[12],此處不再贅述。

    (8)

    在引力輔助坐標(biāo)系中,探測器進入和離開雙曲線剩余速度表示為

    (9)

    (10)

    在引力輔助坐標(biāo)系中,引力輔助行星提供的速度增量可表示為

    (sinδsinψ)j+(cosδ-1)i]

    (11)

    此外,將引力輔助前后的速度差投影至引力輔助坐標(biāo)系,即

    v+-v-=v∞(κ1i+κ2j+κ3k)

    (12)

    將式(11),(12)帶入式(7),可得

    ΔvI=v∞[(κ1-cosδ+1)2+(κ2-sinδsinψ)2+

    (13)

    因此,只要確定剩余速度偏轉(zhuǎn)角δ和甩擺角度參數(shù)ψ,就能求得引力輔助完成后所需施加的脈沖機動速度增量。

    2 轉(zhuǎn)移軌跡設(shè)計

    轉(zhuǎn)移軌跡設(shè)計是空間任務(wù)規(guī)劃的基礎(chǔ)。圓錐曲線拼接法是經(jīng)典的轉(zhuǎn)移軌跡初始設(shè)計方法。采用該方法能快速完成轉(zhuǎn)移軌跡初始設(shè)計。本文基于上述引力輔助簡化模型,采用圓錐曲線拼接法進行復(fù)雜的引力輔助轉(zhuǎn)移軌跡設(shè)計。嵌入引力輔助簡化模型的圓錐曲線拼接法的執(zhí)行過程如圖2所示。

    圖2 嵌入動力引力輔助模型的圓錐曲線拼接法流程圖Fig.2 Flowchart of patched conic method with powered gravity-assist model

    首先,在已知探測器地球出發(fā)時刻、行星借力時刻和目標(biāo)到達(dá)時刻的條件下,通過星歷可求得相關(guān)天體在相應(yīng)時刻的狀態(tài)。然后,基于圓錐曲線拼接法的思想,求解Lambert問題,得到相鄰天體之間的轉(zhuǎn)移軌跡,如圖2中所示的軌跡D和軌跡1。最后,利用引力輔助模型將這些轉(zhuǎn)移軌跡段按順序拼接,并計算引力輔助過程中所需施加的脈沖速度增量,進而得到任務(wù)所需的總速度增量。

    因此,轉(zhuǎn)移軌跡設(shè)計問題被轉(zhuǎn)化為一個多變量優(yōu)化問題,即通過挑選最優(yōu)的任務(wù)時間節(jié)點和引力輔助參數(shù)(近心點高度、甩擺角度等),使整個任務(wù)的燃料消耗最少。對于這種多變量全局優(yōu)化問題,遺傳算法具有明顯優(yōu)勢。所以,本文后續(xù)的軌道設(shè)計將采用遺傳算法進行優(yōu)化,其性能指標(biāo)為探測器提供的總速度增量。

    3 仿真試驗與結(jié)果分析

    3.1 引力輔助方式比較

    以基于單次引力輔助的木星探測任務(wù)為例,對比分析2種簡化引力輔助方案。為方便對比,除引力輔助模型外,其他任務(wù)條件完全一致。探測器于2032年11月8日從地球出發(fā),于2033年4月27日完成金星引力輔助后飛往木星,并最終于2036年9月1日抵達(dá)木星。根據(jù)第2節(jié)描述的轉(zhuǎn)移軌跡設(shè)計方法,分別帶入近心點機動引力輔助模型和甩擺后機動引力輔助模型,計算最終的速度增量。仿真結(jié)果見表1。

    表 1 不同地球引力輔助方案的對比

    由表可知:金星引力輔助過程中,近心點機動所需的速度增量要遠(yuǎn)遠(yuǎn)小于甩擺后機動所需的速度增量,且前者在引力輔助過程中的近心點半徑也較小。這表明:在引力輔助過程中,探測器能從金星獲取更多能量,從而減小自身所需提供的速度增量。由此可見:相比于甩擺后的機動方案,近心點機動方案能節(jié)省更多能量。為驗證這一結(jié)論,針對地球和火星引力輔助任務(wù)進行仿真,結(jié)果見表2。

    表 2 不同情況下的引力輔助模型對比

    相比于表1,表2中引力輔助過程中需要施加的速度增量大大減小,且不同動力引力輔助方案的近心點半徑差異并不明顯。但是,相對于甩擺后機動的引力輔助方案,近心點機動方式的燃料消耗仍然較小,只是其優(yōu)勢并沒有表1中那么明顯。綜上所述,近心點機動的引力輔助方式能使探測器以更少的燃料消耗完成引力輔助前后的軌跡拼接。因此,后續(xù)木星探測任務(wù)轉(zhuǎn)移軌跡的設(shè)計將采用近心點機動引力輔助模型。

    3.2 木星轉(zhuǎn)移軌跡設(shè)計

    木星是太陽系最大的行星,對于研究太陽系的起源具有重要意義。迄今為止,包括NASA,ESA在內(nèi)的多家航天機構(gòu)已發(fā)射了伽利略號、朱諾號等木星探測器。此外,木星探測也已納入了我國未來的深空探測規(guī)劃。

    木星距離地球十分遙遠(yuǎn),現(xiàn)有探測器均采用行星引力輔助技術(shù)轉(zhuǎn)移至木星,且為進一步節(jié)省燃料,往往在抵達(dá)木星前實施多次行星引力輔助。金星-地球-地球引力輔助序列是常見的木星探測任務(wù)方案之一。本文將基于此序列完成木星轉(zhuǎn)移軌跡的初始設(shè)計。結(jié)合我國的木星探測計劃,將發(fā)射窗口搜索區(qū)間設(shè)定為2030年1月1日至2036年12月31日。性能指標(biāo)為探測器提供的總速度增量,即

    J=Δvtotal=‖ΔvD‖+‖ΔvA‖+

    ∑‖ΔvGAk‖

    (14)

    式中:ΔvD為探測器離開地球所需的速度增量;ΔvGAk為第k次引力輔助所需施加的速度增量,其值可通過1.1節(jié)中所述模型求解;ΔvA為探測器到達(dá)木星后完成木星捕獲所需的速度增量。考慮到木星高強度的輻射,將木星目標(biāo)捕獲軌道的近木點半徑設(shè)定為不得小于4倍的木星半徑。另外,為減小木星捕獲所需的速度增量,將捕獲軌道的遠(yuǎn)木點半徑設(shè)定為100倍的木星半徑。因此,探測器的目標(biāo)捕獲軌道為100RJ×4RJ的大橢圓環(huán)木軌道。木星捕獲所需的速度增量為

    (15)

    式中:vA∞為探測器到達(dá)木星時的剩余速度;μJ為木星引力常數(shù);ra和rp分別為目標(biāo)捕獲軌道的近木點半徑和遠(yuǎn)木點半徑。

    為確保探測器在引力輔助過程中的安全,對不同引力輔助行星的引力輔助軌跡的近心點半徑添加必要的約束,約束條件見表3。

    表3 不同行星引力輔助近心點半徑下限

    采用遺傳算法進行尋優(yōu),其進化過程如圖3所示。經(jīng)過25代遺傳搜索的最優(yōu)速度增量為1.170 km/s,最優(yōu)解的詳細(xì)信息見表4。

    圖3 遺傳算法搜索進化圖Fig.3 Search process of genetic algorithm

    探測器于2033年1月11日從地球發(fā)射,經(jīng)過159 d的星際轉(zhuǎn)移,于2033年7月3日到達(dá)金星,完成金星引力輔助后返回地球。2034年6月2日,探測器首次返回地球,完成第一次地球引力輔助。探測器于334 d后再次返回地球,在完成第二次地球引力輔助后,進入前往木星的轉(zhuǎn)移軌道,并于2039年6月12日到達(dá)木星。整個轉(zhuǎn)移過程持續(xù)2 343 d,需要探測器提供的總速度增量為1.17 km/s。任務(wù)轉(zhuǎn)移軌跡如圖4所示。

    表4 遺傳算法搜索結(jié)果詳細(xì)信息表

    圖4 木星探測轉(zhuǎn)移軌跡示意圖Fig.4 Transfer trajectory for Jupiter mission

    4 結(jié)論

    本文對比分析了近心點機動和甩擺后機動2種引力輔助方式。仿真結(jié)果表明:采用行星引力輔助時,在近心點施加脈沖機動比引力輔助完成后施加機動更加節(jié)能。基于此結(jié)論,本文完成了木星多天體引力輔助轉(zhuǎn)移軌跡的初始設(shè)計,其中包含1次金星引力輔助和2次地球引力輔助。本文研究內(nèi)容和結(jié)論對我國后續(xù)木星探測等深空探測任務(wù)具有一定的參考意義。但是,本文提出的模型只適用于脈沖機動方式。隨著電推進等連續(xù)小推力推進技術(shù)的成熟和在工程中的應(yīng)用,針對小推力和引力輔助相結(jié)合的軌道設(shè)計與優(yōu)化將是未來軌道設(shè)計的熱點。因此,后續(xù)將針對小推力引力輔助技術(shù)進行深入研究。

    猜你喜歡
    木星引力機動
    飛向木星
    軍事文摘(2023年12期)2023-06-12 07:50:26
    五彩斑斕的木星
    熱木星是什么星?
    裝載機動臂的疲勞壽命計算
    12萬畝機動地不再“流浪”
    機動三輪車的昨天、今天和明天
    木星
    引力
    初中生(2017年3期)2017-02-21 09:17:40
    感受引力
    A dew drop
    日本熟妇午夜| 色吧在线观看| 男女视频在线观看网站免费| 国产97色在线日韩免费| 中文资源天堂在线| 久久久久久人人人人人| 亚洲激情在线av| 露出奶头的视频| 日韩欧美精品免费久久 | 欧美乱色亚洲激情| 亚洲精品一区av在线观看| 99精品欧美一区二区三区四区| 69av精品久久久久久| 精品熟女少妇八av免费久了| 国产精品电影一区二区三区| 亚洲 欧美 日韩 在线 免费| 国产精品一及| 99riav亚洲国产免费| 国产精品亚洲av一区麻豆| 丁香六月欧美| 嫩草影视91久久| 在线免费观看的www视频| 亚洲av熟女| 亚洲欧美激情综合另类| 欧美高清成人免费视频www| 免费在线观看影片大全网站| 免费高清视频大片| 99在线视频只有这里精品首页| 欧洲精品卡2卡3卡4卡5卡区| 变态另类成人亚洲欧美熟女| 国产免费av片在线观看野外av| 国产精品女同一区二区软件 | 久久亚洲真实| 在线播放国产精品三级| 在线观看66精品国产| 午夜福利欧美成人| 国产精品 国内视频| 久久精品国产99精品国产亚洲性色| 国产欧美日韩精品一区二区| 日韩人妻高清精品专区| 国产精品日韩av在线免费观看| 一本一本综合久久| 免费人成在线观看视频色| 久久国产精品影院| 亚洲av免费在线观看| 亚洲国产精品成人综合色| 在线免费观看不下载黄p国产 | 婷婷丁香在线五月| 99久久九九国产精品国产免费| 日韩欧美精品v在线| 亚洲精品粉嫩美女一区| 午夜激情欧美在线| 亚洲无线观看免费| 亚洲自拍偷在线| 色综合站精品国产| av片东京热男人的天堂| 特级一级黄色大片| 男人舔女人下体高潮全视频| 日本 欧美在线| 精品99又大又爽又粗少妇毛片 | 国产精品野战在线观看| 国产av在哪里看| 一二三四社区在线视频社区8| www.www免费av| 国产精品久久电影中文字幕| 精品不卡国产一区二区三区| 国产精品,欧美在线| 久久精品91无色码中文字幕| 国产在视频线在精品| 婷婷亚洲欧美| www.色视频.com| 狠狠狠狠99中文字幕| 国产亚洲欧美98| 韩国av一区二区三区四区| 国产三级中文精品| 变态另类丝袜制服| 女警被强在线播放| 少妇的逼水好多| 老熟妇仑乱视频hdxx| 十八禁网站免费在线| 毛片女人毛片| 国产精品电影一区二区三区| 亚洲真实伦在线观看| 丰满乱子伦码专区| 老司机午夜十八禁免费视频| 特级一级黄色大片| 男女做爰动态图高潮gif福利片| 久久6这里有精品| 又爽又黄无遮挡网站| 国产成+人综合+亚洲专区| 精品人妻1区二区| 亚洲精品粉嫩美女一区| 最新在线观看一区二区三区| 亚洲欧美日韩高清在线视频| 熟女电影av网| 亚洲熟妇中文字幕五十中出| 观看美女的网站| 国产69精品久久久久777片| 欧美乱色亚洲激情| 久久欧美精品欧美久久欧美| а√天堂www在线а√下载| 蜜桃亚洲精品一区二区三区| 三级男女做爰猛烈吃奶摸视频| 国产精品久久久久久精品电影| 午夜福利在线观看吧| 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 免费在线观看影片大全网站| 国产精品久久久人人做人人爽| 岛国在线免费视频观看| 国产又黄又爽又无遮挡在线| 日本免费一区二区三区高清不卡| 在线观看美女被高潮喷水网站 | 亚洲第一电影网av| 午夜福利高清视频| 18禁美女被吸乳视频| 成熟少妇高潮喷水视频| 身体一侧抽搐| 成年女人永久免费观看视频| 精品无人区乱码1区二区| 熟妇人妻久久中文字幕3abv| 黄片小视频在线播放| 1024手机看黄色片| 又黄又粗又硬又大视频| 国产精品一及| 国内精品久久久久久久电影| 日韩精品中文字幕看吧| 久久久久久久久久黄片| 国产美女午夜福利| 夜夜躁狠狠躁天天躁| 日韩国内少妇激情av| 观看美女的网站| 免费人成视频x8x8入口观看| 亚洲成av人片免费观看| 久99久视频精品免费| 午夜精品在线福利| 欧美日韩亚洲国产一区二区在线观看| 精品国内亚洲2022精品成人| 国产精品98久久久久久宅男小说| 999久久久精品免费观看国产| 悠悠久久av| 国产真人三级小视频在线观看| 久久这里只有精品中国| 人人妻人人澡欧美一区二区| 高清日韩中文字幕在线| 淫秽高清视频在线观看| 一个人免费在线观看电影| 欧美成人a在线观看| 亚洲五月婷婷丁香| 亚洲人成电影免费在线| 免费看日本二区| 99久久成人亚洲精品观看| 亚洲在线自拍视频| 夜夜夜夜夜久久久久| 午夜福利成人在线免费观看| 夜夜躁狠狠躁天天躁| 一本综合久久免费| 男女午夜视频在线观看| 免费大片18禁| 国产在视频线在精品| 亚洲中文字幕日韩| 国产精品三级大全| 亚洲一区高清亚洲精品| 成年版毛片免费区| 男女床上黄色一级片免费看| 两人在一起打扑克的视频| 午夜免费男女啪啪视频观看 | 精品国产三级普通话版| 欧美乱妇无乱码| 久久国产精品人妻蜜桃| 一区二区三区激情视频| 五月玫瑰六月丁香| 桃红色精品国产亚洲av| 亚洲欧美日韩高清在线视频| 悠悠久久av| 国产三级黄色录像| 亚洲精品456在线播放app | 成人鲁丝片一二三区免费| 午夜免费男女啪啪视频观看 | 国产麻豆成人av免费视频| av专区在线播放| 国产黄片美女视频| 99精品欧美一区二区三区四区| 国产精品嫩草影院av在线观看 | 国产精品久久久久久人妻精品电影| 特大巨黑吊av在线直播| 色哟哟哟哟哟哟| 国产黄a三级三级三级人| 日韩亚洲欧美综合| 中文资源天堂在线| 久久亚洲精品不卡| 国内精品久久久久久久电影| 欧美不卡视频在线免费观看| 十八禁网站免费在线| 精品一区二区三区视频在线 | 国产精品精品国产色婷婷| 亚洲一区高清亚洲精品| 国产精品1区2区在线观看.| 精品国产美女av久久久久小说| 亚洲精华国产精华精| 国产激情偷乱视频一区二区| 美女大奶头视频| 在线观看美女被高潮喷水网站 | 高清在线国产一区| 99国产极品粉嫩在线观看| 欧美日韩乱码在线| 美女高潮的动态| 亚洲电影在线观看av| 少妇的逼好多水| 精品一区二区三区视频在线 | 观看美女的网站| 亚洲欧美精品综合久久99| 又黄又粗又硬又大视频| 国产主播在线观看一区二区| 毛片女人毛片| 欧美最黄视频在线播放免费| 亚洲国产精品999在线| 久久精品综合一区二区三区| 国产日本99.免费观看| 在线国产一区二区在线| 99热只有精品国产| 91麻豆av在线| 亚洲精品亚洲一区二区| 日本 欧美在线| 久久久精品大字幕| 成人三级黄色视频| 欧美乱妇无乱码| 国产亚洲精品综合一区在线观看| av欧美777| 夜夜爽天天搞| 日本在线视频免费播放| 日本黄色视频三级网站网址| 最近视频中文字幕2019在线8| 欧美丝袜亚洲另类 | 国产精华一区二区三区| 国产伦精品一区二区三区四那| 亚洲成人免费电影在线观看| 99国产极品粉嫩在线观看| 国产精品久久久久久久久免 | 久久久国产精品麻豆| av视频在线观看入口| 综合色av麻豆| 午夜免费激情av| 欧美成人免费av一区二区三区| 亚洲成av人片在线播放无| 欧美日韩黄片免| 岛国视频午夜一区免费看| 亚洲av免费高清在线观看| 尤物成人国产欧美一区二区三区| 国产精品国产高清国产av| 成人av一区二区三区在线看| 色在线成人网| 国产精品久久久久久精品电影| 亚洲欧美日韩卡通动漫| 在线看三级毛片| 国产色婷婷99| av天堂中文字幕网| 国内精品美女久久久久久| 久久人人精品亚洲av| 午夜免费成人在线视频| 中文字幕人妻丝袜一区二区| 中国美女看黄片| 亚洲avbb在线观看| 中文资源天堂在线| 国产一区二区三区视频了| 欧美在线黄色| 舔av片在线| 精品熟女少妇八av免费久了| av黄色大香蕉| 国产精品乱码一区二三区的特点| 中文字幕av成人在线电影| 波多野结衣高清作品| or卡值多少钱| 99国产精品一区二区蜜桃av| 夜夜夜夜夜久久久久| 小蜜桃在线观看免费完整版高清| 麻豆久久精品国产亚洲av| 久久久久久久久大av| 99riav亚洲国产免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 19禁男女啪啪无遮挡网站| 欧美在线一区亚洲| 国产成人啪精品午夜网站| 婷婷精品国产亚洲av| 国产不卡一卡二| 亚洲av一区综合| 非洲黑人性xxxx精品又粗又长| 中文在线观看免费www的网站| 一个人看的www免费观看视频| 日韩成人在线观看一区二区三区| 90打野战视频偷拍视频| 在线免费观看不下载黄p国产 | 亚洲一区二区三区不卡视频| 国产精品综合久久久久久久免费| 他把我摸到了高潮在线观看| 日本免费a在线| 最近最新免费中文字幕在线| 性欧美人与动物交配| 成人一区二区视频在线观看| 午夜福利在线观看免费完整高清在 | 两个人的视频大全免费| 日韩成人在线观看一区二区三区| 香蕉丝袜av| 欧美丝袜亚洲另类 | 一区二区三区免费毛片| 久久久久久久久中文| 国产高清激情床上av| 欧美日韩福利视频一区二区| 国产精品99久久久久久久久| 丰满人妻一区二区三区视频av | 日本成人三级电影网站| 舔av片在线| 精品久久久久久久久久久久久| 嫩草影视91久久| 美女cb高潮喷水在线观看| 亚洲人成网站高清观看| 中文字幕av成人在线电影| 久久6这里有精品| 亚洲不卡免费看| 日韩有码中文字幕| 国产欧美日韩一区二区精品| 日本在线视频免费播放| 99视频精品全部免费 在线| 色噜噜av男人的天堂激情| 亚洲av日韩精品久久久久久密| 三级毛片av免费| 制服人妻中文乱码| 精品国产超薄肉色丝袜足j| 亚洲av美国av| 91久久精品国产一区二区成人 | 国产精品美女特级片免费视频播放器| 国产伦精品一区二区三区四那| 精品电影一区二区在线| 美女高潮喷水抽搐中文字幕| 一级作爱视频免费观看| 亚洲av电影不卡..在线观看| 久久久久久久久中文| 国产三级中文精品| 中文字幕熟女人妻在线| 国产亚洲欧美98| 好看av亚洲va欧美ⅴa在| 此物有八面人人有两片| 最近最新中文字幕大全免费视频| 午夜福利成人在线免费观看| 免费无遮挡裸体视频| 亚洲精品日韩av片在线观看 | 看黄色毛片网站| 日本成人三级电影网站| 国产精品自产拍在线观看55亚洲| 国产单亲对白刺激| 波野结衣二区三区在线 | 性色av乱码一区二区三区2| 免费在线观看日本一区| 中文字幕人成人乱码亚洲影| 可以在线观看的亚洲视频| 国产精品99久久99久久久不卡| 天天添夜夜摸| 看免费av毛片| 国产精品,欧美在线| 久久久久精品国产欧美久久久| 在线视频色国产色| 熟女人妻精品中文字幕| 91久久精品电影网| 亚洲人成电影免费在线| 麻豆久久精品国产亚洲av| 毛片女人毛片| 国产精品久久久久久久电影 | 少妇的逼水好多| 国产私拍福利视频在线观看| 国产精品久久电影中文字幕| 国产精品av视频在线免费观看| 亚洲欧美激情综合另类| 国产男靠女视频免费网站| 真人一进一出gif抽搐免费| 变态另类成人亚洲欧美熟女| а√天堂www在线а√下载| 欧美日本视频| 亚洲无线在线观看| 国产在线精品亚洲第一网站| 五月玫瑰六月丁香| 色尼玛亚洲综合影院| 亚洲乱码一区二区免费版| 国产日本99.免费观看| 色哟哟哟哟哟哟| 91在线精品国自产拍蜜月 | 黄片小视频在线播放| 桃色一区二区三区在线观看| 久久国产乱子伦精品免费另类| 亚洲人与动物交配视频| 免费在线观看影片大全网站| 精品一区二区三区视频在线 | 国产亚洲av嫩草精品影院| 九九在线视频观看精品| 亚洲aⅴ乱码一区二区在线播放| 国内精品一区二区在线观看| 亚洲狠狠婷婷综合久久图片| 91在线观看av| 亚洲av日韩精品久久久久久密| 国产亚洲精品综合一区在线观看| 琪琪午夜伦伦电影理论片6080| svipshipincom国产片| 久久国产乱子伦精品免费另类| 51国产日韩欧美| 亚洲精品亚洲一区二区| 精品熟女少妇八av免费久了| 免费无遮挡裸体视频| 中亚洲国语对白在线视频| 欧美日韩瑟瑟在线播放| 亚洲五月天丁香| 国产精品国产高清国产av| 久久天躁狠狠躁夜夜2o2o| 亚洲va日本ⅴa欧美va伊人久久| 最近视频中文字幕2019在线8| 人妻丰满熟妇av一区二区三区| 日韩欧美国产一区二区入口| 婷婷亚洲欧美| 久久精品国产亚洲av香蕉五月| 久久久国产成人免费| 看黄色毛片网站| 中文字幕人成人乱码亚洲影| 久久久久国内视频| 久久精品国产亚洲av涩爱 | 亚洲精品久久国产高清桃花| 亚洲av成人精品一区久久| 人人妻,人人澡人人爽秒播| 69av精品久久久久久| 免费在线观看日本一区| 在线视频色国产色| 久久性视频一级片| 内射极品少妇av片p| 黄色女人牲交| 欧美丝袜亚洲另类 | 蜜桃亚洲精品一区二区三区| 亚洲国产精品999在线| 最近在线观看免费完整版| 小说图片视频综合网站| 日韩亚洲欧美综合| 精品电影一区二区在线| 久久精品国产99精品国产亚洲性色| 国产三级在线视频| 欧美色欧美亚洲另类二区| 国产乱人视频| 午夜激情福利司机影院| 免费看a级黄色片| 国产日本99.免费观看| 日韩欧美精品v在线| 欧美3d第一页| 91av网一区二区| 欧美一级毛片孕妇| 久久精品国产清高在天天线| av黄色大香蕉| 欧美黑人欧美精品刺激| 国产在视频线在精品| 变态另类丝袜制服| 中文字幕高清在线视频| 国产精品爽爽va在线观看网站| 久久久久久久亚洲中文字幕 | 90打野战视频偷拍视频| 两个人看的免费小视频| 成人国产一区最新在线观看| 久久亚洲真实| 91在线观看av| 超碰av人人做人人爽久久 | 2021天堂中文幕一二区在线观| 很黄的视频免费| 97碰自拍视频| 亚洲久久久久久中文字幕| 国产高清激情床上av| 日韩欧美在线二视频| 搡女人真爽免费视频火全软件 | 婷婷亚洲欧美| 99国产综合亚洲精品| 国产亚洲精品一区二区www| 偷拍熟女少妇极品色| 99久久精品一区二区三区| 国产精品日韩av在线免费观看| 无限看片的www在线观看| 51午夜福利影视在线观看| 亚洲va日本ⅴa欧美va伊人久久| 免费高清视频大片| 搞女人的毛片| 免费人成视频x8x8入口观看| 亚洲精品一卡2卡三卡4卡5卡| 国产精品 国内视频| 两个人的视频大全免费| 精品熟女少妇八av免费久了| 亚洲人成网站在线播放欧美日韩| 俺也久久电影网| 国产精品一区二区三区四区免费观看 | 在线观看日韩欧美| ponron亚洲| 99精品久久久久人妻精品| 欧美一区二区精品小视频在线| 国产成人aa在线观看| 国产爱豆传媒在线观看| 51国产日韩欧美| 午夜福利在线观看免费完整高清在 | 一进一出抽搐动态| 国产野战对白在线观看| 最新在线观看一区二区三区| 97人妻精品一区二区三区麻豆| 精品久久久久久成人av| 亚洲av成人av| 看片在线看免费视频| 亚洲狠狠婷婷综合久久图片| 又黄又粗又硬又大视频| 黄片大片在线免费观看| 久久久久久人人人人人| 国产视频内射| 一级a爱片免费观看的视频| 国产高清三级在线| 亚洲av成人av| 亚洲成av人片在线播放无| 欧美黑人巨大hd| 日韩中文字幕欧美一区二区| 在线观看免费午夜福利视频| 少妇人妻一区二区三区视频| 亚洲久久久久久中文字幕| 午夜免费观看网址| 亚洲人成网站高清观看| 日韩欧美免费精品| 精品福利观看| 欧美成狂野欧美在线观看| 无遮挡黄片免费观看| 91麻豆av在线| 欧美午夜高清在线| 观看免费一级毛片| 偷拍熟女少妇极品色| 岛国在线观看网站| 一边摸一边抽搐一进一小说| 国产三级中文精品| 亚洲第一欧美日韩一区二区三区| 丰满的人妻完整版| 欧美乱码精品一区二区三区| 欧美zozozo另类| 三级男女做爰猛烈吃奶摸视频| 久久久久九九精品影院| 日韩欧美精品v在线| 在线免费观看不下载黄p国产 | 天天躁日日操中文字幕| 国产午夜精品论理片| 综合色av麻豆| 国产激情偷乱视频一区二区| 一进一出好大好爽视频| 首页视频小说图片口味搜索| 亚洲欧美日韩无卡精品| 国产色婷婷99| 亚洲性夜色夜夜综合| 又黄又粗又硬又大视频| 日韩欧美精品免费久久 | 在线观看66精品国产| 国产亚洲av嫩草精品影院| 日韩人妻高清精品专区| 国产精品永久免费网站| 天堂动漫精品| 9191精品国产免费久久| 午夜免费激情av| 国产三级中文精品| 欧美一级毛片孕妇| 久久久久久国产a免费观看| 日韩国内少妇激情av| 午夜老司机福利剧场| 人人妻人人澡欧美一区二区| 日本三级黄在线观看| 国产精品一区二区三区四区久久| 18+在线观看网站| 嫁个100分男人电影在线观看| 久久久久久九九精品二区国产| 又爽又黄无遮挡网站| 国产精品免费一区二区三区在线| 一区二区三区国产精品乱码| 少妇裸体淫交视频免费看高清| 国产精品乱码一区二三区的特点| 欧美日韩福利视频一区二区| 2021天堂中文幕一二区在线观| 狂野欧美白嫩少妇大欣赏| 国产99白浆流出| 亚洲精华国产精华精| 成人高潮视频无遮挡免费网站| 国产一区二区三区在线臀色熟女| 久久久久久久亚洲中文字幕 | 欧美乱码精品一区二区三区| 在线天堂最新版资源| 日本 欧美在线| 国产精品国产高清国产av| 狠狠狠狠99中文字幕| 欧美绝顶高潮抽搐喷水| 免费看美女性在线毛片视频| 亚洲国产欧美网| 国产一区二区激情短视频| 精品国内亚洲2022精品成人| 五月玫瑰六月丁香| 精品久久久久久久久久免费视频| 高清日韩中文字幕在线| 午夜久久久久精精品| 他把我摸到了高潮在线观看| 久久伊人香网站| 日本 av在线| 成人亚洲精品av一区二区| 色噜噜av男人的天堂激情| 亚洲av电影不卡..在线观看| 国产一区二区在线观看日韩 | 日日摸夜夜添夜夜添小说| 俄罗斯特黄特色一大片| 色哟哟哟哟哟哟| 女人被狂操c到高潮| 18禁国产床啪视频网站| 又黄又爽又免费观看的视频| 免费观看的影片在线观看| 亚洲乱码一区二区免费版| 可以在线观看毛片的网站| 成人午夜高清在线视频| 国产免费一级a男人的天堂| 高清在线国产一区| 亚洲av日韩精品久久久久久密|