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

    基于自適應(yīng)升溫模擬退火算法的農(nóng)業(yè)機(jī)器人全區(qū)域覆蓋策略

    2021-11-15 17:24:58張彥斐宮金良蘭玉彬
    關(guān)鍵詞:模擬退火柵格障礙物

    王 偉,張彥斐,宮金良,蘭玉彬,3

    (1 山東理工大學(xué) 機(jī)械工程學(xué)院,山東 淄博 255000; 2 山東理工大學(xué) 農(nóng)業(yè)工程與食品科學(xué)學(xué)院,山東 淄博 255000; 3 華南農(nóng)業(yè)大學(xué) 電子工程學(xué)院/人工智能學(xué)院,廣東 廣州 510642)

    智慧化無(wú)人農(nóng)業(yè)是將傳感數(shù)據(jù)與智能決策作為管控全局的要素[1],并結(jié)合多傳感信息融合技術(shù)、人機(jī)交互技術(shù)與云邊協(xié)同計(jì)算技術(shù)對(duì)農(nóng)業(yè)機(jī)器人進(jìn)行信息化管理的現(xiàn)代農(nóng)業(yè)[2]。解決農(nóng)業(yè)機(jī)器人在智慧化無(wú)人農(nóng)場(chǎng)中的全區(qū)域覆蓋問(wèn)題可減少機(jī)器人工作區(qū)域的重復(fù)率,達(dá)到對(duì)農(nóng)田工作區(qū)域的全覆蓋,是當(dāng)下急需解決的關(guān)鍵問(wèn)題。

    針對(duì)機(jī)器人全區(qū)域覆蓋的環(huán)境不同,全區(qū)域覆蓋問(wèn)題可分為靜態(tài)已知環(huán)境的覆蓋與動(dòng)態(tài)未知環(huán)境的覆蓋[3-4]。農(nóng)業(yè)機(jī)器人所面對(duì)的大田作業(yè)一般為靜態(tài)已知環(huán)境下的離線預(yù)規(guī)劃地圖覆蓋問(wèn)題,目前已有的研究多通過(guò)柵格法與單元分解法簡(jiǎn)化機(jī)器人全區(qū)域覆蓋的復(fù)雜工作環(huán)境[5]。通過(guò)對(duì)復(fù)雜工作環(huán)境進(jìn)行柵格化處理,完成全區(qū)域覆蓋,可保證機(jī)器人遍歷的精度,如賀利樂(lè)等[6]通過(guò)機(jī)器人本體上的傳感器構(gòu)建基于動(dòng)態(tài)柵格法的工作環(huán)境,從而保證機(jī)器人實(shí)現(xiàn)路徑全覆蓋,但系統(tǒng)開(kāi)展機(jī)器人全區(qū)域覆蓋的計(jì)算量會(huì)隨區(qū)域地圖的增大而呈指數(shù)式增加[7-8];單元分解法則根據(jù)工作區(qū)域中現(xiàn)有障礙物的形狀、大小、位置與區(qū)域分解規(guī)則將整體工作區(qū)域劃分成多個(gè)虛擬子區(qū)域,以完成對(duì)復(fù)雜工作區(qū)域的簡(jiǎn)化[9],如胡詩(shī)宇[10]利用柵格法對(duì)環(huán)境地圖建模,通過(guò)矩形分解的思想對(duì)地圖單元分解,并通過(guò)深度優(yōu)先搜索算法求解子區(qū)域之間的遍歷順序,但深度優(yōu)先搜索算法只適用于求解可行解,不能找到區(qū)域間最優(yōu)遍歷順序。

    本文通過(guò)柵格法與單元分解法相結(jié)合,簡(jiǎn)化農(nóng)業(yè)機(jī)器人的復(fù)雜工作環(huán)境;通過(guò)改進(jìn)的模擬退火算法,求解機(jī)器人在分區(qū)間的最優(yōu)遍歷順序;通過(guò)A*算法與八鄰域搜索法相結(jié)合的方法,進(jìn)行機(jī)器人的跨區(qū)域銜接路徑規(guī)劃,以此完成農(nóng)業(yè)機(jī)器人對(duì)工作區(qū)域的全區(qū)域覆蓋。

    1 柵格化農(nóng)田與矩形分區(qū)

    本文從山東理工大學(xué)蘭玉彬團(tuán)隊(duì)與淄博禾豐種業(yè)科技股份有限公司共建的智慧化無(wú)人農(nóng)場(chǎng)的實(shí)際農(nóng)業(yè)生產(chǎn)環(huán)境出發(fā),定義農(nóng)業(yè)機(jī)器人的復(fù)雜工作環(huán)境,即智慧化無(wú)人農(nóng)場(chǎng)中存在的由縱橫交錯(cuò)的田間道路自然分割的離散分布的數(shù)片農(nóng)田,以及農(nóng)田中存在的一些分散障礙物如農(nóng)田中的風(fēng)力發(fā)電機(jī)、電線桿、固定安裝的傳感器等構(gòu)成了農(nóng)業(yè)機(jī)器人復(fù)雜的工作環(huán)境[11]。

    1.1 柵格化農(nóng)田與障礙物膨脹處理

    根據(jù)本文建立的復(fù)雜農(nóng)田工作環(huán)境概念,定義工作環(huán)境地圖并參考文獻(xiàn)[12]方法對(duì)此區(qū)域進(jìn)行柵格化處理,將農(nóng)業(yè)機(jī)器人的工作范圍設(shè)置為柵格圖的單位長(zhǎng)度。農(nóng)田環(huán)境與農(nóng)田柵格化結(jié)果如圖1a所示。為避免農(nóng)業(yè)機(jī)器人在全遍歷過(guò)程中陷入死角,本文針對(duì)邊界與柵格線不重合的障礙物進(jìn)行二值化膨脹處理[13],處理后的效果如圖1b所示。

    1.2 柵格分區(qū)與分區(qū)合并

    為簡(jiǎn)化農(nóng)業(yè)機(jī)器人進(jìn)行大田作業(yè)時(shí)的復(fù)雜工作環(huán)境,本文將由田間道路自然分割形成的單片農(nóng)田作為一級(jí)工作區(qū)域分區(qū),并在一級(jí)工作區(qū)域分區(qū)內(nèi)根據(jù)農(nóng)田區(qū)域中障礙物的形狀、大小與位置劃分農(nóng)業(yè)機(jī)器人可自由遍歷的二級(jí)工作區(qū)域分區(qū),通過(guò)農(nóng)業(yè)機(jī)器人在一級(jí)分區(qū)與二級(jí)分區(qū)內(nèi)的遍歷完成對(duì)農(nóng)田工作區(qū)域的全區(qū)域覆蓋。

    劃分二級(jí)柵格分區(qū)具體方案為:分別過(guò)障礙物膨脹處理后的左下角頂點(diǎn)與右上角頂點(diǎn)作垂直于X軸的分區(qū)線,遇到其他障礙物或田間道路即停止;分別過(guò)障礙物膨脹處理后的左上角頂點(diǎn)與右下角頂點(diǎn)作垂直于Y軸的分區(qū)線,遇到其他障礙物或田間道路即停止。在圖1b的基礎(chǔ)上進(jìn)行柵格分區(qū)的結(jié)果如圖2a所示。

    圖1 柵格化農(nóng)田建模 (a)與障礙物膨脹處理 (b)Fig.1 Rasterized farmland modeling (a) and obstacle expansion treatment (b)

    圖2 柵格分區(qū) (a)與柵格分區(qū)合并 (b)Fig.2 Grid partition (a) and grid partition merging (b)

    工作區(qū)域中二級(jí)分區(qū)的數(shù)量與機(jī)器人遍歷面積重復(fù)率成正比,為減少機(jī)器人的路徑重復(fù)率,本文在柵格分區(qū)的基礎(chǔ)上進(jìn)行分區(qū)合并操作。一個(gè)柵格分區(qū)與其鄰近分區(qū)共同邊長(zhǎng)度相等則可以進(jìn)行分區(qū)合并[14],分區(qū)時(shí)先合并可縱向合并的分區(qū),再合并可橫向合并的分區(qū)。分區(qū)合并結(jié)果如圖2b所示。

    2 基于模擬退火算法求解旅行商問(wèn)題

    優(yōu)化農(nóng)業(yè)機(jī)器人在分區(qū)間的遍歷順序可減少農(nóng)業(yè)機(jī)器人在分區(qū)間的銜接路徑,進(jìn)而減少其遍歷路徑重復(fù)率、提高工作效率。本文在優(yōu)化各一級(jí)分區(qū)遍歷順序之后再優(yōu)化各一級(jí)分區(qū)內(nèi)二級(jí)分區(qū)的遍歷順序,通過(guò)農(nóng)業(yè)機(jī)器人在一級(jí)分區(qū)與二級(jí)分區(qū)內(nèi)的遍歷,實(shí)現(xiàn)整體工作區(qū)域的全覆蓋。尋找農(nóng)業(yè)機(jī)器人在分區(qū)間的最佳遍歷順序本質(zhì)上是一個(gè)旅行商問(wèn)題[15]。

    2.1 旅行商問(wèn)題

    將圖2b中的各分區(qū)抽象化為一個(gè)個(gè)由分區(qū)重心點(diǎn)代表的質(zhì)點(diǎn),則旅行商問(wèn)題可以描述為:旅行商從某個(gè)質(zhì)點(diǎn)出發(fā)訪問(wèn),尋找一條遍歷全部r個(gè)質(zhì)點(diǎn)且每個(gè)質(zhì)點(diǎn)僅遍歷1次的最短路徑[16],即搜索自然子集N={1,2,…,r}中各元素(質(zhì)點(diǎn)編號(hào))的排序,使得遍歷各質(zhì)點(diǎn)的路徑長(zhǎng)度(R)最?。?/p>

    2.2 模擬退火算法

    模擬退火算法因其隨機(jī)跳變接受新解的特性而具有較好的收斂速度與尋優(yōu)能力[18],故本文通過(guò)模擬退火算法求解旅行商問(wèn)題。算法求解旅行商問(wèn)題時(shí)首先設(shè)定初始溫度、降溫速率與鏈長(zhǎng)等參數(shù),并建立模擬退火算法所得可行解i的適應(yīng)度值為:

    算法在不斷降溫的過(guò)程中通過(guò)多種變換法生成新解,并根據(jù)Metropolis準(zhǔn)則以概率的形式選擇是否保留此解,則算法溫度為T(mén)時(shí)新解i被接受的概率P為:

    3 基于自適應(yīng)升溫的模擬退火算法

    本文提出基于貪婪機(jī)制的優(yōu)質(zhì)可行解生成方法與基于自適應(yīng)升溫的模擬退火算法改進(jìn)方法,以提高模擬退火算法求解旅行商問(wèn)題的尋優(yōu)能力與收斂性能。

    3.1 基于貪婪機(jī)制的優(yōu)質(zhì)可行解生成方法

    傳統(tǒng)模擬退火算法通常采用2變換法與3變換法生成新解,本文在以上2種方法的基礎(chǔ)上借鑒遺傳算法變異操作的思想,提出關(guān)于模擬退火算法的第3種新解生成方法即遺傳變異法。模擬退火算法新解生成方法如圖3所示。

    圖3 模擬退火算法新解生成方法Fig.3 New solution generation methods of simulated annealing algorithm

    2變換法任意選擇2個(gè)原解中的位置,將2個(gè)位置中間的排列順序進(jìn)行倒置操作以生成新解,將此種方法得到的新解適應(yīng)度值記為swap2(f);3變換法任意選擇3個(gè)位置,將前兩個(gè)位置間的數(shù)字置于第3個(gè)數(shù)字后以生成新解,將此種方法得到的新解適應(yīng)度值記為swap3(f);遺傳變異法通過(guò)交換隨機(jī)選擇的2個(gè)數(shù)字的位置以生成新解,將此種方法得到的新解適應(yīng)度值記為swap4(f)。

    本文計(jì)算以上3種方法生成的新解的適應(yīng)度值,并借鑒貪婪機(jī)制逐步尋找局部最優(yōu)以達(dá)到全局最優(yōu)的思想,只保留適應(yīng)度值高的新解對(duì)應(yīng)的狀態(tài)作為算法的實(shí)時(shí)新?tīng)顟B(tài)(New state),即:

    3.2 基于自適應(yīng)升溫的模擬退火算法

    在模擬退火算法陷入局部最優(yōu)解時(shí),人為提高算法運(yùn)行的溫度,有助于提高算法對(duì)于較差解的接受概率,進(jìn)而有更大的可能跳出局部最優(yōu)解。若算法運(yùn)行中產(chǎn)生的可行解的鄰域內(nèi)無(wú)比該解更優(yōu)的可行解,即:則將其定義為算法已陷入局部最優(yōu),式中k表示算法陷入局部時(shí)在當(dāng)前溫度下已迭代的次數(shù),L表示算法鏈長(zhǎng)。

    進(jìn)行局部升溫操作時(shí)若升溫幅度過(guò)小,則較差解的接受概率提升效果有限,達(dá)不到跳出局部最優(yōu)解的效果;若升溫幅度過(guò)大,則算法對(duì)較差解的接受概率接近100%[19],算法尋優(yōu)又會(huì)進(jìn)入漫長(zhǎng)的全局搜索,嚴(yán)重降低算法的收斂速度。故本文建立解集多樣性的概念以設(shè)計(jì)自適應(yīng)升溫策略。算法運(yùn)行過(guò)程中產(chǎn)生的各可行解組成一個(gè)解集,根據(jù)算法所記錄的解集實(shí)時(shí)最優(yōu)適應(yīng)度值與解集實(shí)時(shí)平均適應(yīng)度值的關(guān)系定義當(dāng)前解集多樣性

    4 基于 A*算法的分區(qū)間銜接路徑

    農(nóng)業(yè)機(jī)器人在一個(gè)分區(qū)內(nèi)完成遍歷工作時(shí),需根據(jù)改進(jìn)模擬退火算法規(guī)劃的分區(qū)間最優(yōu)遍歷順序轉(zhuǎn)至下個(gè)分區(qū)起點(diǎn)繼續(xù)遍歷。由于農(nóng)業(yè)機(jī)器人工作環(huán)境的復(fù)雜性,上個(gè)分區(qū)終點(diǎn)與下個(gè)分區(qū)起點(diǎn)可能并不連續(xù),所以,為了提高農(nóng)業(yè)機(jī)器人的工作效率,降低農(nóng)業(yè)機(jī)器人的遍歷路徑重復(fù)率,本文以路徑重復(fù)率為優(yōu)化目標(biāo),通過(guò)A*算法規(guī)劃?rùn)C(jī)器人分區(qū)間的銜接路徑。

    A*算法通過(guò)建立從起始點(diǎn)到目標(biāo)點(diǎn)的評(píng)價(jià)函數(shù)規(guī)劃兩點(diǎn)間的最優(yōu)路徑,評(píng)價(jià)函數(shù)計(jì)算式為:

    啟發(fā)代價(jià)函數(shù)表示一種預(yù)估距離,此函數(shù)的存在能夠保證機(jī)器人始終向目標(biāo)點(diǎn)移動(dòng)[20],本文中啟發(fā)代價(jià)函數(shù)通過(guò)曼哈頓距離實(shí)現(xiàn):

    通過(guò)八鄰域搜索法搜索當(dāng)前點(diǎn)(m,n)的周?chē)鷸鸥?,?(m-1,n+1)、(m,n+1)、(m+1,n+1)、(m+1,n)、(m+1,n-1)、(m,n-1)、(m-1,n-1)、(m-1,n)。將以上柵格作為當(dāng)前點(diǎn)代入評(píng)價(jià)函數(shù)計(jì)算其路徑代價(jià),在路徑代價(jià)最小的柵格基礎(chǔ)上繼續(xù)使用A*算法與八鄰域搜索法重復(fù)以上過(guò)程探索路徑直至目標(biāo)點(diǎn)。

    5 算法仿真與試驗(yàn)驗(yàn)證

    為驗(yàn)證改進(jìn)模擬退火算法的收斂能力與尋優(yōu)能力,以及農(nóng)業(yè)機(jī)器人全區(qū)域覆蓋策略對(duì)路徑重復(fù)率的優(yōu)化效果。本文分別通過(guò)Matlab 2014軟件對(duì)改進(jìn)模擬退火算法與農(nóng)業(yè)機(jī)器人全區(qū)域覆蓋策略進(jìn)行仿真分析。

    5.1 改進(jìn)模擬退火算法仿真

    定義3組分別包含20、30、40個(gè)分區(qū)的重心點(diǎn)坐標(biāo),對(duì)比傳統(tǒng)遺傳算法、模擬退火算法與本文改進(jìn)模擬退火算法規(guī)劃各組最優(yōu)路徑所得到的路徑長(zhǎng)度、算法收斂時(shí)的迭代次數(shù),結(jié)果如表1所示。3個(gè)算法規(guī)劃40個(gè)分區(qū)的最優(yōu)遍歷路徑如圖4所示。

    表1 3 種算法對(duì)不同分區(qū)規(guī)模的路徑規(guī)劃結(jié)果Table 1 Planning results of three algorithms for different partitioning sizes

    由表1可知,在分區(qū)數(shù)量為20、30時(shí),3種算法均能得到相近的路徑長(zhǎng)度,但在分區(qū)數(shù)量為20時(shí),改進(jìn)模擬退火算法收斂時(shí)的迭代次數(shù)較模擬退火算法和傳統(tǒng)遺傳算法分別減少了63.3%和66.7%;在分區(qū)數(shù)量為30時(shí),改進(jìn)模擬退火算法收斂時(shí)的迭代次數(shù)較模擬退火算法和傳統(tǒng)遺傳算法分別減少了35.5%和72.0%。由表1和圖4可知,在分區(qū)數(shù)量為40時(shí),改進(jìn)的模擬退火算法所獲得的路徑長(zhǎng)度分別比模擬退火算法和傳統(tǒng)遺傳算法減少了14.7%和10.1%,收斂時(shí)的迭代次數(shù)分別減少了9.8%和59.1%。

    圖4 3 種算法對(duì) 40 個(gè)分區(qū)的最優(yōu)遍歷路徑規(guī)劃圖Fig.4 The optimal traversal path planning diagram of 40 partitions by three algorithms

    5.2 路徑遍歷仿真

    為減少機(jī)器人的轉(zhuǎn)彎次數(shù),機(jī)器人在分區(qū)內(nèi)遍歷時(shí)沿矩形分區(qū)的長(zhǎng)邊作往返式遍歷。選擇下個(gè)分區(qū)中最靠近上個(gè)分區(qū)終點(diǎn)的單元格作為下個(gè)分區(qū)的起點(diǎn)。農(nóng)業(yè)機(jī)器人在圖2b基礎(chǔ)上進(jìn)行的路徑遍歷仿真如圖5所示。

    圖2b中農(nóng)業(yè)機(jī)器人可自由遍歷的工作面積為189 m2,圖5 中機(jī)器人移動(dòng)總面積為 222 m2,機(jī)器人遍歷路徑重復(fù)率為14.86%,工作區(qū)域覆蓋率接近100%。

    圖5 遍歷路徑規(guī)劃圖Fig.5 Traversal path planning diagram

    5.3 試驗(yàn)驗(yàn)證

    為驗(yàn)證本文全區(qū)域覆蓋策略對(duì)農(nóng)業(yè)機(jī)器人遍歷面積重復(fù)率與覆蓋率的優(yōu)化作用,以山東理工大學(xué)蘭玉彬團(tuán)隊(duì)研發(fā)的高地隙噴藥機(jī)器人為對(duì)象進(jìn)行全區(qū)域覆蓋遍歷試驗(yàn)。

    通過(guò)智慧化無(wú)人農(nóng)場(chǎng)云平臺(tái)人機(jī)交互界面將整體工作區(qū)域劃分為由虛擬道路分割成的4個(gè)子區(qū)域,在整體工作區(qū)域中設(shè)置分散的11個(gè)異形障礙物,并以此為環(huán)境地圖進(jìn)行機(jī)器人路徑規(guī)劃。

    本試驗(yàn)中高地隙噴藥機(jī)器人雙臂噴藥桿展開(kāi)后的工作范圍為10 m,行駛速度為1.7 m/s,試驗(yàn)工作區(qū)域面積為0.8 hm2。高地隙噴藥機(jī)器人、部分異形障礙物放大圖與虛擬道路分割線如圖6所示。

    圖6 農(nóng)業(yè)機(jī)器人試驗(yàn)現(xiàn)場(chǎng)Fig.6 Experimental site of agricultural robot

    根據(jù)高地隙噴藥機(jī)器人遍歷過(guò)程中實(shí)時(shí)回傳至農(nóng)場(chǎng)云平臺(tái)的工作數(shù)據(jù)可知,高地隙噴藥機(jī)器人使用9.1 min完成整體工作區(qū)域的遍歷,實(shí)際遍歷面積為0.927 hm2,遍歷面積重復(fù)率為15.83%,工作區(qū)域覆蓋率接近100%。

    6 結(jié)論

    本文引入遺傳算法變異操作的思想,建立基于貪婪機(jī)制的模擬退火算法優(yōu)質(zhì)可行解生成方法;設(shè)定模擬退火算法陷入局部最優(yōu)解的判斷依據(jù)并建立算法解集多樣性的概念,設(shè)計(jì)基于自適應(yīng)升溫的模擬退火算法改進(jìn)方案,以提高算法的尋優(yōu)能力與收斂速度。仿真試驗(yàn)表明,在分區(qū)規(guī)模為40時(shí),改進(jìn)的模擬退火算法所獲得的路徑長(zhǎng)度分別比模擬退火算法和傳統(tǒng)遺傳算法減少了14.7%和10.1%,收斂時(shí)的迭代次數(shù)分別減少了9.8%和59.1%。

    本文根據(jù)農(nóng)場(chǎng)實(shí)際工作環(huán)境建立一級(jí)分區(qū)的概念,在柵格化環(huán)境建模與障礙物膨脹處理的基礎(chǔ)上,在一級(jí)分區(qū)內(nèi)部建立二級(jí)分區(qū)的柵格分區(qū)和分區(qū)合并規(guī)則。通過(guò)改進(jìn)的模擬退火算法優(yōu)化機(jī)器人在分區(qū)間的遍歷順序,通過(guò)A*算法與八鄰域搜索法相結(jié)合進(jìn)行機(jī)器人跨區(qū)域銜接路徑規(guī)劃。仿真試驗(yàn)表明,機(jī)器人遍歷路徑重復(fù)率為14.86%。機(jī)器人現(xiàn)場(chǎng)遍歷試驗(yàn)表明,高地隙噴藥機(jī)器人遍歷面積重復(fù)率為15.83%,工作區(qū)域覆蓋率接近100%。

    猜你喜歡
    模擬退火柵格障礙物
    基于鄰域柵格篩選的點(diǎn)云邊緣點(diǎn)提取方法*
    高低翻越
    SelTrac?CBTC系統(tǒng)中非通信障礙物的設(shè)計(jì)和處理
    模擬退火遺傳算法在機(jī)械臂路徑規(guī)劃中的應(yīng)用
    基于模糊自適應(yīng)模擬退火遺傳算法的配電網(wǎng)故障定位
    SOA結(jié)合模擬退火算法優(yōu)化電容器配置研究
    不同剖面形狀的柵格壁對(duì)柵格翼氣動(dòng)特性的影響
    基于遺傳-模擬退火算法的城市軌道交通快慢車(chē)停站方案
    基于CVT排布的非周期柵格密度加權(quán)陣設(shè)計(jì)
    土釘墻在近障礙物的地下車(chē)行通道工程中的應(yīng)用
    免费在线观看亚洲国产| 成人三级黄色视频| 最近在线观看免费完整版| 宅男免费午夜| 成人性生交大片免费视频hd| 欧美绝顶高潮抽搐喷水| 久久久久久久精品吃奶| 国产精品久久久久久久电影 | 久久久色成人| 精品久久久久久,| 麻豆一二三区av精品| 长腿黑丝高跟| 久久精品亚洲精品国产色婷小说| 国产亚洲av嫩草精品影院| 97超视频在线观看视频| 91九色精品人成在线观看| 18美女黄网站色大片免费观看| 99国产极品粉嫩在线观看| 国产久久久一区二区三区| 久久久精品欧美日韩精品| 白带黄色成豆腐渣| 久久久国产成人精品二区| 午夜福利免费观看在线| 日日夜夜操网爽| 桃红色精品国产亚洲av| 老司机在亚洲福利影院| 欧美3d第一页| 国产精品免费一区二区三区在线| 国产精品1区2区在线观看.| 国产精品久久电影中文字幕| 搡老岳熟女国产| 国产99白浆流出| 日日摸夜夜添夜夜添小说| 特大巨黑吊av在线直播| 国产美女午夜福利| 国产高清有码在线观看视频| 日韩三级视频一区二区三区| 国产精品香港三级国产av潘金莲| 午夜免费成人在线视频| 精品电影一区二区在线| 黄色 视频免费看| 国产精品98久久久久久宅男小说| 日本在线视频免费播放| 免费在线观看成人毛片| 老鸭窝网址在线观看| 成在线人永久免费视频| 999精品在线视频| 黄色日韩在线| 丁香六月欧美| xxx96com| 精品久久久久久成人av| 禁无遮挡网站| 免费看美女性在线毛片视频| 国产精品日韩av在线免费观看| 亚洲片人在线观看| 好男人电影高清在线观看| 国产精品久久久久久久电影 | 性欧美人与动物交配| 亚洲精品久久国产高清桃花| 日本免费a在线| 色综合站精品国产| 亚洲,欧美精品.| avwww免费| 婷婷精品国产亚洲av在线| 身体一侧抽搐| 真实男女啪啪啪动态图| 久久久久国产一级毛片高清牌| 欧美日韩福利视频一区二区| 久久久久九九精品影院| 欧美性猛交╳xxx乱大交人| 欧美黑人欧美精品刺激| 少妇丰满av| 午夜福利欧美成人| 日韩成人在线观看一区二区三区| 欧美又色又爽又黄视频| 欧美日韩综合久久久久久 | 久久香蕉国产精品| 国产一区二区在线观看日韩 | 中文资源天堂在线| 九色国产91popny在线| 久久久久精品国产欧美久久久| 亚洲国产中文字幕在线视频| 97超级碰碰碰精品色视频在线观看| 欧美日韩福利视频一区二区| 日韩有码中文字幕| 亚洲五月婷婷丁香| 国产成人影院久久av| 狠狠狠狠99中文字幕| 国内毛片毛片毛片毛片毛片| 国产一区二区激情短视频| 他把我摸到了高潮在线观看| 精品国产美女av久久久久小说| 久久国产乱子伦精品免费另类| 国产野战对白在线观看| 国产精品影院久久| 午夜福利高清视频| 免费在线观看影片大全网站| 欧美日韩精品网址| 又黄又粗又硬又大视频| 成人永久免费在线观看视频| 国产精品美女特级片免费视频播放器 | 99热6这里只有精品| 欧洲精品卡2卡3卡4卡5卡区| 免费av不卡在线播放| 免费一级毛片在线播放高清视频| 日韩欧美国产在线观看| 精品99又大又爽又粗少妇毛片 | 国产成人精品久久二区二区91| 观看美女的网站| 欧美黄色片欧美黄色片| 丁香欧美五月| 波多野结衣高清无吗| 精品国产超薄肉色丝袜足j| 成人特级黄色片久久久久久久| 久久国产精品影院| 日本黄色视频三级网站网址| 亚洲人成网站在线播放欧美日韩| 日本成人三级电影网站| 亚洲午夜精品一区,二区,三区| 日韩精品中文字幕看吧| 网址你懂的国产日韩在线| 国产黄片美女视频| 精品一区二区三区视频在线观看免费| 久久国产乱子伦精品免费另类| 18禁观看日本| 国产亚洲精品久久久久久毛片| 99视频精品全部免费 在线 | 狠狠狠狠99中文字幕| 欧美性猛交╳xxx乱大交人| 91av网站免费观看| 性色av乱码一区二区三区2| 日本a在线网址| 18美女黄网站色大片免费观看| www.www免费av| 日本一二三区视频观看| 亚洲五月婷婷丁香| 国产欧美日韩一区二区精品| 别揉我奶头~嗯~啊~动态视频| 欧美日韩一级在线毛片| 精品一区二区三区四区五区乱码| 91麻豆av在线| 亚洲五月天丁香| 久久久久久国产a免费观看| 巨乳人妻的诱惑在线观看| 国产免费av片在线观看野外av| 亚洲av美国av| 成年女人永久免费观看视频| 亚洲中文av在线| 国产一区二区三区在线臀色熟女| 一个人免费在线观看电影 | 久久热在线av| 最近最新中文字幕大全电影3| 日韩精品中文字幕看吧| 亚洲午夜理论影院| 国产伦在线观看视频一区| 亚洲午夜精品一区,二区,三区| 国产乱人伦免费视频| 黄色成人免费大全| 身体一侧抽搐| 大型黄色视频在线免费观看| 一进一出抽搐gif免费好疼| 午夜视频精品福利| 久久精品人妻少妇| 亚洲精品在线观看二区| 日韩 欧美 亚洲 中文字幕| 精品国产美女av久久久久小说| 麻豆国产97在线/欧美| 欧美日韩乱码在线| 在线观看午夜福利视频| 1024手机看黄色片| 亚洲国产色片| 成人鲁丝片一二三区免费| 成人18禁在线播放| 午夜日韩欧美国产| 757午夜福利合集在线观看| 不卡一级毛片| av中文乱码字幕在线| 男人舔女人的私密视频| 国产精品久久久av美女十八| 99久久精品热视频| 日韩三级视频一区二区三区| a级毛片a级免费在线| 国产精品av久久久久免费| 天天躁日日操中文字幕| 18禁裸乳无遮挡免费网站照片| 国产亚洲精品久久久com| 香蕉国产在线看| 欧美中文日本在线观看视频| 真实男女啪啪啪动态图| 99国产精品99久久久久| 成年版毛片免费区| 我的老师免费观看完整版| 亚洲av成人不卡在线观看播放网| 亚洲真实伦在线观看| 亚洲国产看品久久| 国产一区二区在线av高清观看| 欧美性猛交╳xxx乱大交人| 91久久精品国产一区二区成人 | 18禁国产床啪视频网站| 老汉色av国产亚洲站长工具| 老汉色av国产亚洲站长工具| 欧美激情在线99| 在线观看免费午夜福利视频| 成人18禁在线播放| 国产亚洲精品av在线| 三级男女做爰猛烈吃奶摸视频| av在线蜜桃| 亚洲精品一区av在线观看| 一进一出好大好爽视频| 欧美黑人欧美精品刺激| а√天堂www在线а√下载| 午夜免费激情av| 首页视频小说图片口味搜索| 九色国产91popny在线| 欧美日本视频| 99视频精品全部免费 在线 | 好男人在线观看高清免费视频| 成年免费大片在线观看| 伦理电影免费视频| 可以在线观看毛片的网站| 成人永久免费在线观看视频| 精品国产超薄肉色丝袜足j| 一区福利在线观看| 琪琪午夜伦伦电影理论片6080| 一夜夜www| 国产精品美女特级片免费视频播放器 | 在线观看免费午夜福利视频| 亚洲欧美精品综合久久99| 日韩大尺度精品在线看网址| 一本综合久久免费| 亚洲专区国产一区二区| 人人妻人人澡欧美一区二区| 18禁国产床啪视频网站| 99国产精品一区二区蜜桃av| 久9热在线精品视频| 特大巨黑吊av在线直播| 国产精品野战在线观看| 久久国产精品影院| 欧美av亚洲av综合av国产av| 亚洲国产欧洲综合997久久,| 99久久无色码亚洲精品果冻| 男人的好看免费观看在线视频| 少妇的丰满在线观看| 国产精品电影一区二区三区| 久久国产精品影院| 丁香欧美五月| 国产成人福利小说| 亚洲aⅴ乱码一区二区在线播放| 久久这里只有精品19| 伦理电影免费视频| 又黄又粗又硬又大视频| 成人av一区二区三区在线看| avwww免费| 露出奶头的视频| 少妇裸体淫交视频免费看高清| 国产欧美日韩一区二区精品| 成人亚洲精品av一区二区| 免费在线观看日本一区| 国产精品亚洲美女久久久| av天堂中文字幕网| 亚洲色图 男人天堂 中文字幕| 老司机福利观看| 老司机深夜福利视频在线观看| 男人舔女人的私密视频| 一区福利在线观看| 偷拍熟女少妇极品色| 在线观看午夜福利视频| 两人在一起打扑克的视频| 亚洲在线观看片| 亚洲国产精品合色在线| 免费看a级黄色片| 亚洲av免费在线观看| 久久香蕉国产精品| 欧美xxxx黑人xx丫x性爽| av片东京热男人的天堂| 国产精品电影一区二区三区| 亚洲电影在线观看av| 1024香蕉在线观看| 舔av片在线| 国产麻豆成人av免费视频| 人妻丰满熟妇av一区二区三区| 欧美不卡视频在线免费观看| 午夜免费成人在线视频| 成年女人永久免费观看视频| 国产成人一区二区三区免费视频网站| 怎么达到女性高潮| 一二三四在线观看免费中文在| 欧美日本亚洲视频在线播放| 亚洲熟妇中文字幕五十中出| 啪啪无遮挡十八禁网站| 精品无人区乱码1区二区| 日韩欧美国产在线观看| 亚洲av五月六月丁香网| 亚洲国产精品999在线| 国产伦人伦偷精品视频| 中文字幕精品亚洲无线码一区| 成人18禁在线播放| 少妇熟女aⅴ在线视频| 伊人久久大香线蕉亚洲五| 亚洲,欧美精品.| 精品久久久久久久末码| 久久这里只有精品19| 国内精品久久久久精免费| 在线a可以看的网站| 中文字幕人成人乱码亚洲影| 精品免费久久久久久久清纯| 99热6这里只有精品| 国产激情欧美一区二区| 一二三四在线观看免费中文在| 亚洲成人精品中文字幕电影| 琪琪午夜伦伦电影理论片6080| 禁无遮挡网站| 亚洲欧洲精品一区二区精品久久久| 免费无遮挡裸体视频| 午夜影院日韩av| 偷拍熟女少妇极品色| aaaaa片日本免费| 国产av麻豆久久久久久久| 麻豆久久精品国产亚洲av| 99精品欧美一区二区三区四区| 国产成人精品久久二区二区91| 色视频www国产| 91在线观看av| 国产三级中文精品| 噜噜噜噜噜久久久久久91| 成人av在线播放网站| 国产探花在线观看一区二区| 老汉色∧v一级毛片| 九色成人免费人妻av| 给我免费播放毛片高清在线观看| 精品久久久久久久久久久久久| 精品福利观看| av中文乱码字幕在线| 国产一区二区三区在线臀色熟女| 成人鲁丝片一二三区免费| 国产免费男女视频| 亚洲自拍偷在线| 天天添夜夜摸| 美女 人体艺术 gogo| 午夜日韩欧美国产| 亚洲乱码一区二区免费版| 天天添夜夜摸| 久久久久国产一级毛片高清牌| 亚洲自偷自拍图片 自拍| 一a级毛片在线观看| 国产成人精品无人区| 国产精品 欧美亚洲| 麻豆成人av在线观看| 99精品欧美一区二区三区四区| 精品熟女少妇八av免费久了| 亚洲在线自拍视频| 熟女电影av网| 丰满人妻一区二区三区视频av | 亚洲片人在线观看| 又大又爽又粗| 精品一区二区三区av网在线观看| 亚洲精品在线美女| 国产黄a三级三级三级人| 欧美日韩福利视频一区二区| xxx96com| 欧美一区二区精品小视频在线| 视频区欧美日本亚洲| 可以在线观看的亚洲视频| 精品国产乱子伦一区二区三区| 午夜视频精品福利| 成人特级av手机在线观看| 精品无人区乱码1区二区| 超碰成人久久| 久久亚洲精品不卡| 丝袜人妻中文字幕| 男女那种视频在线观看| 一区二区三区激情视频| 五月玫瑰六月丁香| 老司机午夜十八禁免费视频| 亚洲一区二区三区不卡视频| 国产一区二区在线观看日韩 | 精品人妻1区二区| 九九热线精品视视频播放| 国产视频内射| 曰老女人黄片| 国内揄拍国产精品人妻在线| 麻豆av在线久日| 久久热在线av| 一进一出抽搐gif免费好疼| 女生性感内裤真人,穿戴方法视频| 成人亚洲精品av一区二区| 香蕉久久夜色| 国产一区二区在线观看日韩 | 国产主播在线观看一区二区| 最近最新免费中文字幕在线| 最新美女视频免费是黄的| 午夜两性在线视频| 亚洲av免费在线观看| 久久这里只有精品19| 亚洲自偷自拍图片 自拍| 中文字幕人成人乱码亚洲影| 成人无遮挡网站| 中文字幕久久专区| 日本五十路高清| 18禁美女被吸乳视频| 高潮久久久久久久久久久不卡| 一进一出抽搐gif免费好疼| 国产精品女同一区二区软件 | 特大巨黑吊av在线直播| 亚洲专区国产一区二区| 午夜免费观看网址| 久久亚洲真实| 黄色丝袜av网址大全| 日韩免费av在线播放| 亚洲精品粉嫩美女一区| 成年免费大片在线观看| 国产精品一区二区精品视频观看| 美女黄网站色视频| 国产精品美女特级片免费视频播放器 | 999久久久精品免费观看国产| or卡值多少钱| 国产精品 欧美亚洲| 人人妻人人澡欧美一区二区| 亚洲色图 男人天堂 中文字幕| 黄色日韩在线| or卡值多少钱| 1024手机看黄色片| 非洲黑人性xxxx精品又粗又长| 成年免费大片在线观看| 国产高清视频在线播放一区| 国产精品久久视频播放| 中文资源天堂在线| 琪琪午夜伦伦电影理论片6080| 午夜精品一区二区三区免费看| 免费电影在线观看免费观看| 熟女少妇亚洲综合色aaa.| 美女大奶头视频| 99热6这里只有精品| 日韩精品青青久久久久久| 在线观看免费午夜福利视频| 两个人的视频大全免费| 香蕉国产在线看| 国产亚洲精品久久久久久毛片| 91久久精品国产一区二区成人 | 久久久国产成人免费| 国产午夜福利久久久久久| 一边摸一边抽搐一进一小说| 免费在线观看日本一区| 十八禁网站免费在线| 淫妇啪啪啪对白视频| 搡老岳熟女国产| 老司机午夜十八禁免费视频| 日本黄大片高清| 色播亚洲综合网| 精品久久久久久久毛片微露脸| 非洲黑人性xxxx精品又粗又长| 哪里可以看免费的av片| 久久久久久久午夜电影| 亚洲精品乱码久久久v下载方式 | 又粗又爽又猛毛片免费看| netflix在线观看网站| av视频在线观看入口| 午夜精品久久久久久毛片777| netflix在线观看网站| 床上黄色一级片| 免费看日本二区| av国产免费在线观看| 中文亚洲av片在线观看爽| 免费电影在线观看免费观看| 国产视频内射| 免费在线观看影片大全网站| 麻豆av在线久日| 日韩精品中文字幕看吧| 成人特级av手机在线观看| 日韩欧美国产在线观看| 此物有八面人人有两片| 男女视频在线观看网站免费| 黑人操中国人逼视频| 床上黄色一级片| 国内精品久久久久精免费| www日本黄色视频网| 国产麻豆成人av免费视频| 欧美日韩黄片免| 老鸭窝网址在线观看| 国产成+人综合+亚洲专区| 九九热线精品视视频播放| 观看免费一级毛片| 12—13女人毛片做爰片一| 天天躁日日操中文字幕| 国产精品电影一区二区三区| 色吧在线观看| 欧美极品一区二区三区四区| 亚洲精品在线美女| 欧美成人性av电影在线观看| 色综合婷婷激情| 99在线人妻在线中文字幕| 精品久久久久久久久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 一本一本综合久久| 免费无遮挡裸体视频| 亚洲九九香蕉| 99精品在免费线老司机午夜| 九九久久精品国产亚洲av麻豆 | 国产高清三级在线| 丁香欧美五月| 久久久久性生活片| 午夜视频精品福利| 高清毛片免费观看视频网站| 久久久久国产精品人妻aⅴ院| 99久久99久久久精品蜜桃| 免费看a级黄色片| 精品熟女少妇八av免费久了| 身体一侧抽搐| 18美女黄网站色大片免费观看| 国产精品乱码一区二三区的特点| 最近视频中文字幕2019在线8| 天堂av国产一区二区熟女人妻| 天堂动漫精品| 日韩欧美国产在线观看| 日本 av在线| 久久精品国产99精品国产亚洲性色| 国产极品精品免费视频能看的| 麻豆av在线久日| 亚洲av五月六月丁香网| 成人午夜高清在线视频| 国产熟女xx| 琪琪午夜伦伦电影理论片6080| 亚洲人成网站高清观看| 黄色女人牲交| 又爽又黄无遮挡网站| 国产野战对白在线观看| 老司机午夜十八禁免费视频| 国产av在哪里看| 99在线视频只有这里精品首页| 精品国产亚洲在线| cao死你这个sao货| 亚洲精品一卡2卡三卡4卡5卡| 老鸭窝网址在线观看| 特级一级黄色大片| 露出奶头的视频| 欧美乱码精品一区二区三区| 国产亚洲欧美在线一区二区| 99热精品在线国产| 久久亚洲精品不卡| av黄色大香蕉| 国产欧美日韩一区二区三| 熟女少妇亚洲综合色aaa.| 美女午夜性视频免费| 久久久久国产一级毛片高清牌| 美女免费视频网站| 国产v大片淫在线免费观看| 搡老岳熟女国产| 色综合欧美亚洲国产小说| 女人高潮潮喷娇喘18禁视频| 国内精品久久久久久久电影| 97碰自拍视频| 黄色片一级片一级黄色片| 日韩免费av在线播放| 午夜精品久久久久久毛片777| 每晚都被弄得嗷嗷叫到高潮| 麻豆国产av国片精品| 又黄又爽又免费观看的视频| 岛国在线免费视频观看| 成人国产综合亚洲| 亚洲国产精品成人综合色| 亚洲欧美精品综合一区二区三区| 嫩草影院入口| 免费av不卡在线播放| 欧美xxxx黑人xx丫x性爽| 亚洲欧美日韩东京热| 亚洲七黄色美女视频| 国产成人精品久久二区二区免费| or卡值多少钱| 国产伦一二天堂av在线观看| 日本免费a在线| 亚洲国产欧美网| www.999成人在线观看| 免费观看精品视频网站| 88av欧美| 国产精品电影一区二区三区| 亚洲无线在线观看| 99国产精品99久久久久| 九九热线精品视视频播放| 中文资源天堂在线| 我要搜黄色片| 欧美日韩精品网址| 狠狠狠狠99中文字幕| 欧美日韩瑟瑟在线播放| 久久久久久久精品吃奶| 99国产精品一区二区三区| 国产精品自产拍在线观看55亚洲| 香蕉久久夜色| 男女之事视频高清在线观看| www日本在线高清视频| 欧美成人性av电影在线观看| 手机成人av网站| 99久久无色码亚洲精品果冻| 在线观看66精品国产| 亚洲国产欧美一区二区综合| 国内精品一区二区在线观看| 亚洲黑人精品在线| 一个人看的www免费观看视频| 丝袜人妻中文字幕| 成年女人看的毛片在线观看| 久久久国产成人免费| 国产乱人视频| 夜夜夜夜夜久久久久| 国产一区二区在线观看日韩 | 国产成人影院久久av| 1000部很黄的大片| 99在线视频只有这里精品首页| 亚洲专区中文字幕在线| 美女高潮的动态| 男女下面进入的视频免费午夜| 亚洲第一电影网av| 亚洲熟妇中文字幕五十中出| 国产欧美日韩精品一区二区| 久久久久久国产a免费观看| 亚洲国产欧美网| 免费大片18禁|