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

    求解最優(yōu)月球軟著陸軌道的隱式打靶法

    2019-08-15 03:01:18彭坤彭睿黃震張柏楠
    航空學(xué)報(bào) 2019年7期
    關(guān)鍵詞:狀態(tài)變量著陸器邊值問(wèn)題

    彭坤,彭睿,黃震,張柏楠

    1.中國(guó)空間技術(shù)研究院 載人航天總體部,北京 100094 2.北京航空航天大學(xué),北京 100083

    近年來(lái),各航天大國(guó)紛紛將月球探測(cè)與開(kāi)發(fā)作為研究重點(diǎn)[1-2]。而實(shí)施月球軟著陸是月球資源開(kāi)發(fā)、原位資源利用以及載人登月任務(wù)的基礎(chǔ)。月球軟著陸過(guò)程一般分為霍曼轉(zhuǎn)移段和動(dòng)力下降段[3-4],其中動(dòng)力下降段會(huì)消耗大量燃料,因此需設(shè)計(jì)燃料消耗最優(yōu)動(dòng)力下降段著陸軌道。

    月球軟著陸軌道設(shè)計(jì)是一種終端時(shí)刻自由的最優(yōu)控制問(wèn)題,其主要求解方法分為3大類[5-6]:間接法、直接法和混合法。間接法是最早用來(lái)求解月球軟著陸軌道的方法,其利用龐特亞金極大值原理將最優(yōu)控制問(wèn)題轉(zhuǎn)化為兩點(diǎn)邊值問(wèn)題進(jìn)行求解[7]。該方法的不足是兩點(diǎn)邊值問(wèn)題的求解會(huì)因初值敏感而很難收斂。更多的學(xué)者應(yīng)用直接法求解月球軟著陸軌道,利用配點(diǎn)法[8-9]、函數(shù)逼近方法[10-11]以及偽譜法[12-14]將軌道狀態(tài)變量和控制變量離散化并作為優(yōu)化變量,再采用序列二次規(guī)劃或智能尋優(yōu)算法進(jìn)行求解。直接法雖然收斂性好,但其求解時(shí)間隨計(jì)算精度的提高而增加;同時(shí)其是一種次優(yōu)解法,燃料消耗非最優(yōu)?;旌戏▌t利用間接法的最優(yōu)控制律,將伴隨變量初值作為優(yōu)化變量,采用非線性規(guī)劃[15]或智能優(yōu)化算法[16-17]進(jìn)行求解。其求解結(jié)果相比于直接法更接近最優(yōu)解,同時(shí)計(jì)算量小;相比于間接法降低了初值敏感性。但混合法的尋優(yōu)時(shí)間仍遠(yuǎn)高于間接法的迭代求解時(shí)間,同時(shí)其舍棄了部分最優(yōu)性條件,求解結(jié)果不能嚴(yán)格保證最優(yōu)性,還需要與間接法結(jié)果進(jìn)行比對(duì)。因此,只有間接法能求得理論最優(yōu)解,同時(shí)其一旦能收斂,求解速度最快,計(jì)算精度最高。

    為此,學(xué)者們開(kāi)始重新研究通過(guò)間接法來(lái)求解最優(yōu)月球軟著陸軌道問(wèn)題,研究重點(diǎn)是對(duì)終端時(shí)刻的修正和收斂性強(qiáng)的兩點(diǎn)邊值問(wèn)題求解方法。趙吉松等[18]采用線性攝動(dòng)法[19]同時(shí)對(duì)伴隨變量初值和終端時(shí)刻進(jìn)行修正。肖堯等[20]引入一種時(shí)間尺度變換方法[21]來(lái)解決終端時(shí)刻自由問(wèn)題,借用MATLAB的bvp4c函數(shù)求解月球軟著陸兩點(diǎn)邊值問(wèn)題。

    為提高求解兩點(diǎn)邊值問(wèn)題的收斂性,本文建立了月球軟著陸的歸一化系統(tǒng)模型,減少狀態(tài)變量之間和伴隨變量之間的量級(jí)差距。同時(shí),采用時(shí)間尺度變換方法將終端時(shí)刻自由問(wèn)題轉(zhuǎn)化為時(shí)間區(qū)間為[0,1]的終端時(shí)刻固定問(wèn)題。進(jìn)而把終端時(shí)刻看作狀態(tài)變量并引入哈密爾頓函數(shù)終端值作為隱式終端條件,將原顯式終端條件的兩點(diǎn)邊值問(wèn)題轉(zhuǎn)化為隱式終端條件的兩點(diǎn)邊值問(wèn)題。最后,采用變分法對(duì)隱式終端條件進(jìn)行處理,提出一種隱式打靶法并對(duì)該問(wèn)題進(jìn)行求解,得到燃料消耗最優(yōu)的月球軟著陸軌道,并分析了月球軟著陸過(guò)程中燃料消耗與著陸器發(fā)動(dòng)機(jī)推重比之間的關(guān)系。

    1 系統(tǒng)模型

    忽略動(dòng)力下降段的橫向偏差,采用二維極坐標(biāo)系描述月球軟著陸過(guò)程,并以初始時(shí)刻月面著陸器的月心距和質(zhì)量為參考量,可建立月球軟著陸過(guò)程的歸一化狀態(tài)方程為[17]

    (1)

    (2)

    (3)

    (4)

    設(shè)月球軟著陸過(guò)程發(fā)動(dòng)機(jī)一直開(kāi)機(jī),則月面著陸器的燃料消耗與飛行時(shí)間成正比,故可設(shè)性能指標(biāo)為

    (5)

    2 隱式終端條件兩點(diǎn)邊值問(wèn)題

    2.1 最優(yōu)控制問(wèn)題向兩點(diǎn)邊值問(wèn)題的轉(zhuǎn)化

    利用龐特亞金極大值原理可得哈密爾頓函數(shù)H為

    (6)

    1) 最優(yōu)推力方向角u控制律為

    tanu=(-λvr)/(-λvθ)

    (7)

    2) 伴隨變量方程為

    (8)

    3) 由橫截條件[22]可推導(dǎo)出終端條件為

    (9)

    (10)

    兩點(diǎn)邊值問(wèn)題的邊界條件為

    (11)

    2.2 顯式終端條件向隱式終端條件的轉(zhuǎn)化

    (12)

    (13)

    (14)

    該邊界條件的最后一個(gè)終端約束與其他終端約束不同,它不是某個(gè)狀態(tài)變量的終端值,而是包含某些狀態(tài)變量的一個(gè)等式約束,屬于隱式終端條件。

    3 隱式打靶法

    3.1 隱式終端條件處理

    (15)

    隱式終端條件可簡(jiǎn)單表示為

    (16)

    定義g(1)的變分[22]為

    δg(1)=g*(1)-g(1)=-g(1)

    (17)

    式中:g*(1)為終端約束的真值,即g*(1)=0;g(1)為迭代計(jì)算值,存在誤差。另一方面,通過(guò)變分法則可得

    (18)

    式中:

    (19)

    (20)

    (21)

    (22)

    A5x0=-g(1)

    (23)

    式中:

    (24)

    通過(guò)共軛函數(shù)法還可得到4個(gè)顯式終端條件的關(guān)系式,再加入式(23)可得

    (25)

    式中:

    (26)

    (27)

    3.2 算法流程

    隱式打靶法求解月球軟著陸軌道兩點(diǎn)邊值問(wèn)題的步驟如下:

    3) 判斷終止條件

    (28)

    是否滿足,若滿足,則停止迭代;若滿足第1分式,則求得最優(yōu)解Y*(τ)=Y(k)(τ);若滿足第2分式,則迭代失敗,輸出失敗信息。否則,轉(zhuǎn)第4)步。

    (29)

    4 仿真驗(yàn)證

    月球軟著陸初始條件[10]為r0=1 753 km,θ0=0°,vr0=0 m/s,vθ0=(9.65×10-4·r0)m/s。終端約束為rf=1 738 km,vrf=0 m/s,vθf(wàn)=0 m/s。月球引力常數(shù)為μL=4.902 75×1012m3/s2。以載人月面著陸器為例,設(shè)著陸器質(zhì)量為m0=18 t,發(fā)動(dòng)機(jī)推力為Ft=50 000 N,排氣速度為w=310×9.8 m/s。

    由表1可知,隱式打靶法可以搜索到最優(yōu)月球軟著陸軌道,收斂精度高,歸一化的終端約束精度可以達(dá)到10-10量級(jí),且收斂時(shí)間短,僅為1 s左右。最終求得的載人月面著陸器的月球軟著陸過(guò)程總時(shí)間為486.5 s,燃料消耗為8 007.719 kg。圖1 給出了歸一化系統(tǒng)模型中4個(gè)狀態(tài)變量的變化曲線,圖2給出了對(duì)應(yīng)的4個(gè)伴隨變量的變化曲線,圖3給出了軟著陸過(guò)程中推力方向角的變化曲線,圖4給出了間接法構(gòu)造的哈密爾頓函數(shù)的變化曲線。由圖1可知,所求得的軌道嚴(yán)格滿足軟著陸的月心距、徑向速度、橫向速度約束。由圖2和圖4可知,所求得的軌道嚴(yán)格滿足間接法的式(9)的橫截條件要求,證明了其最優(yōu)性。同時(shí)由圖2可知,由于系統(tǒng)模型的歸一化處理,伴隨變量的量級(jí)相差較小,更有利于問(wèn)題求解過(guò)程的收斂。由圖3可知,推力方向角初始值在180°左右,說(shuō)明λvr(0)<0和λvθ(0)>0,且|λvθ(0)|≥|λvr(0)|,驗(yàn)證了初值猜測(cè)的合理性。

    將隱式打靶法與直接法和混合法進(jìn)行對(duì)比分析。直接法采用文獻(xiàn)[10]的插值函數(shù)逼近法將控制函數(shù)參數(shù)化,將燃料消耗作為性能指標(biāo),終端目標(biāo)約束作為懲罰項(xiàng),采用人工免疫算法進(jìn)行求解?;旌戏ú捎梦墨I(xiàn)[17]的方法,將伴隨變量初值和飛行時(shí)間作為控制變量,評(píng)價(jià)函數(shù)與直接法相同,也采用人工免疫算法進(jìn)行求解。3種方法的求解結(jié)果如表2所示。由表可看出,隱式打靶法的終端約束誤差比直接法和混合法小4~6個(gè)數(shù)量級(jí),收斂時(shí)間也比直接法和混合法低2~3個(gè)數(shù)量級(jí)。直接法求解的燃料消耗大于隱式打靶法和混合法?;旌戏ㄓ捎诳刂坡膳c隱式打靶法相同,優(yōu)化結(jié)果接近隱式打靶法結(jié)果?;旌戏ㄇ蠼鈺r(shí)將終端約束設(shè)為評(píng)價(jià)函數(shù)的懲罰項(xiàng)予以考慮,用隨機(jī)搜索算法人工免疫算法進(jìn)行搜索,在前期搜索中能快速收斂到最優(yōu)解附近,但在后期搜索中局部收斂過(guò)程比較慢,無(wú)法將終端約束誤差控制得特別低。而隱式打靶法直接通過(guò)終端約束誤差修正控制變量,可以得到較低誤差的終端約束,故收斂精度高于混合法。

    表1 隱式打靶法搜索結(jié)果Table 1 Search results of implicit shooting method

    圖1 歸一化狀態(tài)變量隨時(shí)間變化曲線Fig.1 Temporal variation curves of normalization state variables

    圖2 伴隨變量隨時(shí)間變化曲線Fig.2 Temporal variation curves of adjoint variables

    圖3 推力方向角隨時(shí)間變化曲線Fig.3 Temporal variation curves of thrust direction angle

    圖4 哈密爾頓函數(shù)隨時(shí)間變化曲線Fig.4 Temporal variation curves of Hamiltonian function

    表2 隱式打靶法與直接法和混合法搜索結(jié)果對(duì)比Table 2 Comparisons of search results of implicit shooting method with direct method and hybrid method

    為進(jìn)一步驗(yàn)證隱式打靶法的收斂性,對(duì)η∈[0.1, 1.0](η=Ft/(m0g)為發(fā)動(dòng)機(jī)推重比,g為

    表3 隱式打靶法與顯式打靶法搜索結(jié)果對(duì)比

    Table 3 Comparisons of search results of implicit shooting method with explicit shooting method

    方法隱式打靶法顯式打靶法燃料消耗Δm/kg8007.7198007.985月心距r(tf)/m1737999.9991737737.201徑向速度vr(tf)/(m·s-1)-1.437×10-8-6.887×10-1橫向速度vθ(tf)/(m·s-1)-4.936×10-9-1.035×10-11收斂時(shí)間tc/min0.020.01

    圖5 燃料消耗隨推重比變化曲線Fig.5 Variation curve of fuel consumption with thrust-weight ratio

    5 結(jié) 論

    1) 時(shí)間尺度變換法和本文推導(dǎo)的隱式打靶法的組合能夠快速而精確地求解出終端時(shí)刻自由的最優(yōu)月球軟著陸軌道,與直接法和混合法相比,隱式打靶法優(yōu)化精度高,收斂速度快。

    2) 系統(tǒng)模型的歸一化處理減小了伴隨變量的量級(jí)差,進(jìn)而縮小了伴隨變量的猜測(cè)范圍,同時(shí)結(jié)合隱式打靶法中伴隨變量初值和終端時(shí)刻的聯(lián)合修正,大大降低了迭代過(guò)程的敏感性,從而提高了間接法的收斂性。

    3) 本文提出的隱式打靶法收斂性強(qiáng),可適應(yīng)參數(shù)值大范圍變化的月球軟著陸問(wèn)題求解,也可應(yīng)用于其他終端時(shí)刻自由的軌道優(yōu)化問(wèn)題求解。

    4) 對(duì)于最優(yōu)月球軟著陸軌道,推重比η∈[0.1,1.0]時(shí),燃料消耗隨推重比增加先急劇減小后緩慢增大,推重比為0.6時(shí)燃料消耗最小,該變化趨勢(shì)可為載人月面著陸器下降級(jí)發(fā)動(dòng)機(jī)選型提供參考。

    猜你喜歡
    狀態(tài)變量著陸器邊值問(wèn)題
    一階動(dòng)態(tài)電路零狀態(tài)響應(yīng)公式的通用拓展
    非線性n 階m 點(diǎn)邊值問(wèn)題正解的存在性
    基于TwinCAT3控制系統(tǒng)的YB518型小盒透明紙包裝機(jī)運(yùn)行速度的控制分析
    基于嵌套思路的飽和孔隙-裂隙介質(zhì)本構(gòu)理論
    帶有積分邊界條件的奇異攝動(dòng)邊值問(wèn)題的漸近解
    嫦娥四號(hào)巡視器拍攝的著陸器圖像
    航天器工程(2019年4期)2019-11-11 03:41:00
    中國(guó)完成首次火星探測(cè)任務(wù)著陸器懸停避障試驗(yàn)
    伙伴(2019年12期)2019-01-17 04:38:56
    基于多工況的新型著陸器軟著陸性能優(yōu)化
    Mining the Moon Becomes a Serious Prospect
    Recent Development and Emerged Technologies of High-Tc Superconducting Coated Conductors
    国产亚洲精品第一综合不卡 | 18禁动态无遮挡网站| 亚洲精品日韩av片在线观看| 日韩欧美一区视频在线观看| 国产亚洲精品第一综合不卡 | 久久久久久久大尺度免费视频| 九九久久精品国产亚洲av麻豆| 女的被弄到高潮叫床怎么办| 最近中文字幕高清免费大全6| 26uuu在线亚洲综合色| 精品久久蜜臀av无| 欧美bdsm另类| 人人妻人人澡人人爽人人夜夜| 久久99热这里只频精品6学生| 亚洲国产日韩一区二区| 纯流量卡能插随身wifi吗| 精品国产国语对白av| 色94色欧美一区二区| 亚洲精品中文字幕在线视频| 久久久久网色| 丝袜脚勾引网站| 久久ye,这里只有精品| 中文字幕人妻丝袜制服| 精品99又大又爽又粗少妇毛片| 99九九线精品视频在线观看视频| 菩萨蛮人人尽说江南好唐韦庄| 久久热精品热| 精品一区二区三卡| 一本一本综合久久| 青春草亚洲视频在线观看| 精品久久蜜臀av无| 久久精品国产自在天天线| 国产淫语在线视频| 大香蕉97超碰在线| a级毛片免费高清观看在线播放| 久久久久视频综合| 少妇的逼好多水| 成人毛片60女人毛片免费| 午夜免费观看性视频| 欧美精品亚洲一区二区| 午夜福利视频精品| 一区二区三区四区激情视频| 国产精品欧美亚洲77777| 人妻夜夜爽99麻豆av| 青春草视频在线免费观看| 久久久久久久国产电影| 欧美亚洲日本最大视频资源| 18禁动态无遮挡网站| 少妇 在线观看| 国产视频首页在线观看| 日本与韩国留学比较| 国产欧美亚洲国产| 精品少妇黑人巨大在线播放| 国产淫语在线视频| 国产精品人妻久久久久久| 色94色欧美一区二区| av.在线天堂| 少妇被粗大的猛进出69影院 | 国产精品一区二区在线观看99| 亚洲精品乱码久久久v下载方式| 成人二区视频| 十八禁网站网址无遮挡| 国产伦精品一区二区三区视频9| 人人妻人人添人人爽欧美一区卜| 狂野欧美激情性bbbbbb| 午夜精品国产一区二区电影| 香蕉精品网在线| 性色avwww在线观看| 亚洲欧洲日产国产| 国产熟女欧美一区二区| 日日啪夜夜爽| 亚洲伊人久久精品综合| 久久99蜜桃精品久久| 国产熟女欧美一区二区| 欧美3d第一页| 插阴视频在线观看视频| 国产 精品1| 天堂俺去俺来也www色官网| 美女国产高潮福利片在线看| 成人18禁高潮啪啪吃奶动态图 | 欧美日韩成人在线一区二区| 欧美xxⅹ黑人| 久久99蜜桃精品久久| av一本久久久久| 王馨瑶露胸无遮挡在线观看| 亚洲精品久久午夜乱码| 免费久久久久久久精品成人欧美视频 | 男女边摸边吃奶| 最新的欧美精品一区二区| 久久热精品热| 大香蕉久久网| 国产精品人妻久久久久久| av专区在线播放| av黄色大香蕉| 欧美激情国产日韩精品一区| 亚洲av国产av综合av卡| 精品一区二区免费观看| 久久99一区二区三区| 国产又色又爽无遮挡免| av在线app专区| 日本91视频免费播放| kizo精华| 国产男女超爽视频在线观看| 国国产精品蜜臀av免费| 日韩 亚洲 欧美在线| 成人午夜精彩视频在线观看| 69精品国产乱码久久久| 能在线免费看毛片的网站| 少妇精品久久久久久久| 久久99一区二区三区| 最新中文字幕久久久久| 日本色播在线视频| 成人亚洲精品一区在线观看| 久久av网站| 黄色一级大片看看| 久久久久国产网址| 夜夜看夜夜爽夜夜摸| 欧美日韩在线观看h| 男女国产视频网站| 在线观看免费高清a一片| 久久人人爽av亚洲精品天堂| 久久影院123| 高清av免费在线| 国产亚洲最大av| 精品一区在线观看国产| 日日爽夜夜爽网站| 黄色欧美视频在线观看| 亚洲欧美一区二区三区黑人 | 丰满少妇做爰视频| 日日爽夜夜爽网站| 久久精品久久精品一区二区三区| 久久人人爽av亚洲精品天堂| av不卡在线播放| 亚洲人成网站在线播| 日韩在线高清观看一区二区三区| 人妻制服诱惑在线中文字幕| 人体艺术视频欧美日本| 亚洲精品aⅴ在线观看| 免费观看在线日韩| 色网站视频免费| 亚洲av日韩在线播放| 国产精品一国产av| 日韩大片免费观看网站| 欧美日韩视频高清一区二区三区二| 国产精品人妻久久久影院| 少妇被粗大猛烈的视频| 亚洲精品日韩在线中文字幕| 亚洲一区二区三区欧美精品| 日本av免费视频播放| 亚州av有码| 日本wwww免费看| 视频区图区小说| 亚洲伊人久久精品综合| 激情五月婷婷亚洲| av国产精品久久久久影院| 好男人视频免费观看在线| 夜夜骑夜夜射夜夜干| 久久久久久久精品精品| 亚洲,一卡二卡三卡| 亚洲国产色片| 国语对白做爰xxxⅹ性视频网站| 免费日韩欧美在线观看| 亚洲精华国产精华液的使用体验| 精品少妇黑人巨大在线播放| 久久精品国产亚洲av涩爱| av有码第一页| 黄色一级大片看看| 国产欧美亚洲国产| 22中文网久久字幕| 亚洲av成人精品一二三区| 日韩制服骚丝袜av| 97超视频在线观看视频| 极品人妻少妇av视频| 国产精品免费大片| 黄色怎么调成土黄色| 亚洲欧美日韩另类电影网站| 亚洲精品456在线播放app| 三级国产精品片| 97精品久久久久久久久久精品| 亚洲五月色婷婷综合| 在线亚洲精品国产二区图片欧美 | 天堂中文最新版在线下载| 夫妻性生交免费视频一级片| 满18在线观看网站| 免费观看在线日韩| 免费观看av网站的网址| 国产一区二区三区av在线| 国产男女内射视频| 街头女战士在线观看网站| 又大又黄又爽视频免费| 精品国产一区二区久久| 中文字幕最新亚洲高清| 国产在视频线精品| 在线观看www视频免费| 在线观看一区二区三区激情| 国产深夜福利视频在线观看| 乱人伦中国视频| 精品视频人人做人人爽| 亚洲欧美成人精品一区二区| 日本免费在线观看一区| 中国三级夫妇交换| 国产精品久久久久久精品古装| 久久亚洲国产成人精品v| 色婷婷av一区二区三区视频| 欧美亚洲日本最大视频资源| 母亲3免费完整高清在线观看 | av在线老鸭窝| 国产高清三级在线| 国产在线一区二区三区精| 80岁老熟妇乱子伦牲交| 80岁老熟妇乱子伦牲交| 欧美激情国产日韩精品一区| 成人毛片60女人毛片免费| 男女高潮啪啪啪动态图| 国产乱来视频区| 久久久久久久久久久丰满| 亚洲精品乱码久久久v下载方式| 超碰97精品在线观看| 午夜福利视频在线观看免费| 久久精品熟女亚洲av麻豆精品| 高清视频免费观看一区二区| 亚洲综合色惰| 99热这里只有是精品在线观看| 三级国产精品片| 亚洲经典国产精华液单| 国产日韩一区二区三区精品不卡 | 国产成人精品一,二区| 99久久精品一区二区三区| 国产色婷婷99| 人人妻人人爽人人添夜夜欢视频| 国产探花极品一区二区| 久久99一区二区三区| 亚洲国产av新网站| 久久精品夜色国产| 性高湖久久久久久久久免费观看| 免费高清在线观看视频在线观看| 男人操女人黄网站| 国产欧美日韩综合在线一区二区| 欧美日韩精品成人综合77777| 亚洲怡红院男人天堂| 免费大片18禁| 国产乱人偷精品视频| 中文天堂在线官网| 欧美另类一区| 国产成人aa在线观看| 久久国产精品男人的天堂亚洲 | 午夜激情av网站| 一级毛片我不卡| 欧美亚洲日本最大视频资源| a级片在线免费高清观看视频| 黑丝袜美女国产一区| 七月丁香在线播放| 免费av中文字幕在线| 日日爽夜夜爽网站| 久久99热6这里只有精品| 黑人巨大精品欧美一区二区蜜桃 | 久久久久久久大尺度免费视频| 80岁老熟妇乱子伦牲交| 满18在线观看网站| 又大又黄又爽视频免费| 午夜免费鲁丝| 欧美精品亚洲一区二区| 久久久精品免费免费高清| 三级国产精品欧美在线观看| 最黄视频免费看| 一区二区三区四区激情视频| 狂野欧美白嫩少妇大欣赏| 欧美成人精品欧美一级黄| 精品人妻在线不人妻| 免费高清在线观看日韩| 亚洲欧洲国产日韩| 国产高清国产精品国产三级| 熟女人妻精品中文字幕| 乱人伦中国视频| 国产午夜精品一二区理论片| 中文字幕人妻丝袜制服| 日韩中字成人| 制服人妻中文乱码| 最后的刺客免费高清国语| 久久久久久久大尺度免费视频| 亚洲久久久国产精品| 国产精品国产三级专区第一集| 日韩中字成人| 国产男女内射视频| 亚洲精品久久午夜乱码| 欧美 亚洲 国产 日韩一| 久久久国产欧美日韩av| 少妇被粗大猛烈的视频| 在现免费观看毛片| 精品一品国产午夜福利视频| 亚洲精品视频女| xxxhd国产人妻xxx| 国产欧美日韩一区二区三区在线 | 亚洲第一区二区三区不卡| 欧美日韩在线观看h| 在线观看三级黄色| 亚洲国产精品成人久久小说| 亚洲少妇的诱惑av| av.在线天堂| 国产一区二区在线观看av| 综合色丁香网| 精品一区二区三区视频在线| 国产亚洲一区二区精品| 亚洲人成网站在线观看播放| 亚洲av中文av极速乱| 丰满乱子伦码专区| 热99国产精品久久久久久7| 免费少妇av软件| 成人亚洲精品一区在线观看| 蜜臀久久99精品久久宅男| 少妇的逼好多水| videossex国产| 99热这里只有精品一区| 欧美激情国产日韩精品一区| 91aial.com中文字幕在线观看| 一边亲一边摸免费视频| 黄色毛片三级朝国网站| 如何舔出高潮| 一级毛片aaaaaa免费看小| 国产成人精品久久久久久| 国产男女内射视频| 狂野欧美白嫩少妇大欣赏| 最近中文字幕2019免费版| 久久久精品免费免费高清| 日本-黄色视频高清免费观看| 高清黄色对白视频在线免费看| 九九久久精品国产亚洲av麻豆| 久久人人爽av亚洲精品天堂| 韩国高清视频一区二区三区| kizo精华| 精品国产一区二区久久| 国产精品久久久久久精品电影小说| 国产白丝娇喘喷水9色精品| 欧美性感艳星| 午夜免费男女啪啪视频观看| 中文字幕最新亚洲高清| 97超碰精品成人国产| 欧美国产精品一级二级三级| 午夜免费男女啪啪视频观看| 高清av免费在线| 我要看黄色一级片免费的| 国产高清三级在线| 久久久国产精品麻豆| 日本色播在线视频| 亚洲av福利一区| 国产精品一区www在线观看| 男女无遮挡免费网站观看| 这个男人来自地球电影免费观看 | 99热全是精品| 看十八女毛片水多多多| 在线观看免费日韩欧美大片 | 高清毛片免费看| 国产免费现黄频在线看| 最近最新中文字幕免费大全7| 国语对白做爰xxxⅹ性视频网站| 97在线人人人人妻| 边亲边吃奶的免费视频| 国产av国产精品国产| 99久久精品一区二区三区| 欧美日韩在线观看h| 久久国产亚洲av麻豆专区| 美女中出高潮动态图| 精品久久久久久久久亚洲| 五月天丁香电影| 日韩熟女老妇一区二区性免费视频| 制服诱惑二区| 国产黄片视频在线免费观看| 香蕉精品网在线| 国产免费现黄频在线看| 欧美日韩在线观看h| 亚洲情色 制服丝袜| 如何舔出高潮| 人人妻人人澡人人爽人人夜夜| 日韩三级伦理在线观看| 亚洲精品乱码久久久v下载方式| av在线观看视频网站免费| 亚洲精品av麻豆狂野| 99久久人妻综合| 亚洲成人av在线免费| 桃花免费在线播放| 亚洲国产精品成人久久小说| 国产精品国产av在线观看| 国产一区二区在线观看日韩| 亚洲国产精品一区三区| 国产精品国产三级国产专区5o| 久久av网站| 一本色道久久久久久精品综合| 建设人人有责人人尽责人人享有的| 精品人妻在线不人妻| 欧美精品一区二区大全| 欧美少妇被猛烈插入视频| 狂野欧美白嫩少妇大欣赏| 免费看光身美女| 亚洲精品日韩av片在线观看| 日本av免费视频播放| 少妇 在线观看| 亚洲av欧美aⅴ国产| 一级毛片我不卡| 久久精品国产亚洲av涩爱| 一级毛片 在线播放| 久久久久网色| 国产精品国产av在线观看| 国产黄频视频在线观看| 制服人妻中文乱码| 成年av动漫网址| 久久人妻熟女aⅴ| av天堂久久9| 亚洲国产毛片av蜜桃av| 精品午夜福利在线看| 国产精品久久久久久精品古装| 高清午夜精品一区二区三区| 777米奇影视久久| 狂野欧美激情性bbbbbb| 日日摸夜夜添夜夜添av毛片| 最黄视频免费看| 国产 精品1| 成年人免费黄色播放视频| 各种免费的搞黄视频| 国产69精品久久久久777片| 国产白丝娇喘喷水9色精品| 亚洲欧美日韩卡通动漫| 波野结衣二区三区在线| 成人无遮挡网站| 美女脱内裤让男人舔精品视频| 一区二区三区免费毛片| 免费观看av网站的网址| 狂野欧美白嫩少妇大欣赏| 在线观看免费日韩欧美大片 | 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 十分钟在线观看高清视频www| 99久久综合免费| 国产av一区二区精品久久| 亚洲天堂av无毛| 亚洲婷婷狠狠爱综合网| 在线播放无遮挡| 观看av在线不卡| 亚洲色图综合在线观看| 亚洲图色成人| av.在线天堂| 午夜激情久久久久久久| 国产69精品久久久久777片| 国产成人精品福利久久| 午夜激情久久久久久久| 欧美成人精品欧美一级黄| 日本黄大片高清| 欧美国产精品一级二级三级| 建设人人有责人人尽责人人享有的| 欧美人与善性xxx| 欧美日韩精品成人综合77777| 夜夜骑夜夜射夜夜干| 午夜免费男女啪啪视频观看| 草草在线视频免费看| 久久精品国产亚洲网站| 欧美+日韩+精品| 日韩av在线免费看完整版不卡| 美女xxoo啪啪120秒动态图| 欧美日韩亚洲高清精品| 肉色欧美久久久久久久蜜桃| av黄色大香蕉| 国产精品99久久99久久久不卡 | 人妻 亚洲 视频| 亚洲国产成人一精品久久久| 边亲边吃奶的免费视频| 99九九在线精品视频| 日韩,欧美,国产一区二区三区| 成人国产av品久久久| 亚洲精华国产精华液的使用体验| 久久久久精品久久久久真实原创| a级毛色黄片| 男的添女的下面高潮视频| 最后的刺客免费高清国语| 99国产精品免费福利视频| 国产综合精华液| 国产不卡av网站在线观看| 夫妻午夜视频| 少妇精品久久久久久久| 女人精品久久久久毛片| av线在线观看网站| 亚洲美女视频黄频| 纯流量卡能插随身wifi吗| 成人黄色视频免费在线看| 日韩一本色道免费dvd| 国产亚洲午夜精品一区二区久久| 亚洲美女搞黄在线观看| 哪个播放器可以免费观看大片| 欧美精品亚洲一区二区| 最近中文字幕2019免费版| 亚洲高清免费不卡视频| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久久a久久爽久久v久久| 最新的欧美精品一区二区| 国产精品久久久久成人av| 国产在线免费精品| 99九九线精品视频在线观看视频| tube8黄色片| 3wmmmm亚洲av在线观看| 色婷婷久久久亚洲欧美| 亚洲av二区三区四区| 男女边吃奶边做爰视频| 精品亚洲成a人片在线观看| 亚洲欧洲国产日韩| 亚洲精品国产av蜜桃| 欧美日韩精品成人综合77777| 欧美精品亚洲一区二区| 日本黄大片高清| 国产成人a∨麻豆精品| 日日摸夜夜添夜夜添av毛片| 啦啦啦在线观看免费高清www| 99热这里只有是精品在线观看| 性高湖久久久久久久久免费观看| 一级毛片电影观看| 夜夜爽夜夜爽视频| 日韩不卡一区二区三区视频在线| 在线观看www视频免费| 久久免费观看电影| 午夜激情av网站| 国产高清不卡午夜福利| 激情五月婷婷亚洲| 国产成人精品一,二区| 五月伊人婷婷丁香| 国产 精品1| 永久免费av网站大全| 久久国内精品自在自线图片| av视频免费观看在线观看| 亚洲内射少妇av| videos熟女内射| 男女高潮啪啪啪动态图| 人人妻人人澡人人爽人人夜夜| 亚洲婷婷狠狠爱综合网| 欧美日韩综合久久久久久| 国产免费一区二区三区四区乱码| 国产成人午夜福利电影在线观看| 国产精品秋霞免费鲁丝片| 久热久热在线精品观看| 日产精品乱码卡一卡2卡三| 2021少妇久久久久久久久久久| 十八禁网站网址无遮挡| 黑人猛操日本美女一级片| 秋霞伦理黄片| 狂野欧美激情性xxxx在线观看| h视频一区二区三区| 日本爱情动作片www.在线观看| 中国美白少妇内射xxxbb| 国产成人精品婷婷| 又黄又爽又刺激的免费视频.| 久久99蜜桃精品久久| 亚洲国产色片| 99视频精品全部免费 在线| 99九九线精品视频在线观看视频| 大又大粗又爽又黄少妇毛片口| av免费在线看不卡| 久久精品国产亚洲av天美| 3wmmmm亚洲av在线观看| 999精品在线视频| 国产精品人妻久久久影院| 一个人免费看片子| 精品午夜福利在线看| 国产亚洲av片在线观看秒播厂| 日日撸夜夜添| 国产在线一区二区三区精| 草草在线视频免费看| 在线亚洲精品国产二区图片欧美 | 欧美日韩视频高清一区二区三区二| 成年女人在线观看亚洲视频| 有码 亚洲区| av在线老鸭窝| 色婷婷av一区二区三区视频| av黄色大香蕉| 最近最新中文字幕免费大全7| 久久久久久久亚洲中文字幕| 日韩中文字幕视频在线看片| 欧美一级a爱片免费观看看| 男人操女人黄网站| 国产69精品久久久久777片| 国产永久视频网站| 丝袜脚勾引网站| 欧美精品人与动牲交sv欧美| 国产精品久久久久久精品电影小说| 精品熟女少妇av免费看| 精品卡一卡二卡四卡免费| 亚洲成人av在线免费| 一本大道久久a久久精品| 男女免费视频国产| 999精品在线视频| 色哟哟·www| 亚洲欧美成人精品一区二区| 色视频在线一区二区三区| 国产欧美亚洲国产| 在线观看美女被高潮喷水网站| 两个人免费观看高清视频| 2018国产大陆天天弄谢| 日韩欧美一区视频在线观看| 欧美国产精品一级二级三级| 99久久精品一区二区三区| 午夜激情久久久久久久| 成人综合一区亚洲| 亚洲精品色激情综合| 一区二区日韩欧美中文字幕 | 国产在线免费精品| 高清不卡的av网站| 波野结衣二区三区在线| av国产精品久久久久影院| 黄色一级大片看看| 亚洲精品久久午夜乱码| av有码第一页| 精品卡一卡二卡四卡免费| 一级毛片黄色毛片免费观看视频| 精品国产国语对白av| 国产免费福利视频在线观看| 热re99久久精品国产66热6| 国产国拍精品亚洲av在线观看| 成年人免费黄色播放视频| 99热6这里只有精品| 日韩不卡一区二区三区视频在线| 久久精品国产鲁丝片午夜精品|