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

    隨機時變下帶時間窗的取送貨車輛路徑問題優(yōu)化研究

    2022-03-21 06:49:10張歆悅合肥工業(yè)大學管理學院安徽合肥230009
    物流科技 2022年3期
    關鍵詞:車場模擬退火時變

    靳 鵬,張歆悅 (合肥工業(yè)大學 管理學院,安徽 合肥 230009)

    0 引 言

    車輛路徑問題(Vehicle Routing Problem,VRP) 自提出以來,一直是物流組織優(yōu)化領域中的研究難點。在實際物流系統(tǒng)中,為滿足取貨或配送及服務時間窗限制的客戶需求,帶時間窗的取送貨車輛路徑問題(Pickup and Delivery with Time Windows,PDPTW) 逐漸受到眾多學者的關注。然而,近年來城市交通日趨擁堵,受交通管制、交通事故等不確定因素的影響,城市交通路網存在一定的時變性和隨機性,車輛的行駛速度是時變的。在此背景下,研究隨機時變下帶時間窗的取送貨車輛路徑問題(Stochastic Time-dependent Pickup and Delivery with Time Windows, STDPDPTW) 具有重要意義。

    帶時間窗的取送貨車輛路徑問題(PDPTW) 是經典車輛路徑問題(Vehicle Routing Problem, VRP) 的擴展,屬于NP-hard問題,最早由Wilson提出。Bent 等人首次提出兩階段混合算法,第一階段選用模擬退火算法優(yōu)化路徑數(shù)量,第二階段選用大鄰域搜索算法降低路徑成本,有效求解PDPTW 問題。Al Chami 等人使用貪婪隨機的自適應搜索方法生成初始解,后使用遺傳算法優(yōu)化解決方案,有效提高求解PDPTW 問題的效率。目前用于求解帶時間窗的取送貨車輛路徑問題的啟發(fā)式算法主要包括遺傳算法、模擬退火算法、自適應大鄰域搜索算法、蟻群算法等。

    然而,城市化進程逐步加快,道路交通超負荷運轉,加之交通管制和交通事故等因素的影響,在物流運輸過程中,考慮車輛行駛時間的時變性和隨機性更貼近實際。隨機時變下車輛路徑問題(Stochastic Time-dependent Vehicle Routing Problem,STDVRP) 最早由隨機時變最優(yōu)路徑問題(Stochastic Time-dependent Optimal Path Problem, STDOPP) 發(fā)展而來,Hang 等以行駛時間負效用函數(shù)計算隨機時變路徑行駛時間,隨機路網中的最佳路徑受隨機依賴性影響。張春苗指出以車輛行駛時間的概率分布為基礎,網絡拓撲結構復雜,算法求解時間復雜度隨路網規(guī)模增大呈指數(shù)增長。由Sim提出魯棒優(yōu)化模型,將隨機時變路網轉化為確定型時變,曹慧等將該方法應用于隨機時變路網的最優(yōu)路徑問題,結果證明滿足隨機一致性條件,車輛行駛時間可依據確定型時變路網計算。Duan也指出魯棒優(yōu)化時間模型適用于實際的道路網絡。段征宇等人建立STDVRP 魯棒優(yōu)化模型,考慮客戶的時間窗約束,增大了求解的復雜度,并設計改進蟻群算法進行有效求解。

    綜上,多數(shù)文獻獨立研究STDVRP 和PDPTW,同時考慮隨機時變、時間窗約束和取送貨的車輛路徑問題的研究較少,本文建立了隨機時變下帶時間窗的取送貨車輛路徑問題(STDPDPTW) 模型,基于魯棒優(yōu)化方法計算隨機時變路網下車輛行駛時間,為快速且按時滿足客戶需求,以車輛總行駛時間最小化為目標函數(shù)建立混合整數(shù)規(guī)劃模型,并設計一種混合遺傳模擬退火算法進行求解。

    1 問題描述與數(shù)學模型

    1.1 問題描述。STDPDPTW 研究一定數(shù)量的車輛從地理位置不同的車場出發(fā)在客戶能接受的時間范圍內為客戶提供取貨或配送服務,受道路類型和行駛時段的影響,同一路段可能具有不同的行駛速度,制定總行駛時間最小的路徑規(guī)劃方案。對研究問題做如下描述:(1) 配送網絡中有多個車場,每個車場有一定數(shù)量相同類型的運輸車輛,每輛車的起止點為同一車場;(2)每個客戶僅由一輛車進行服務;(3) 客戶點的需求量不大于車輛載重量;(4) 每個客戶點有服務時間窗的要求,車輛可早于最早開始服務時間到達客戶點,等待直至最早開始服務時間進行服務,但不能晚于最遲開始時間到達。

    1.3 問題建模。以服務完所有客戶所需的總行駛時間最小化為目標函數(shù),則隨機時變下帶時間窗的取送貨車輛路徑問題的數(shù)學模型建立如下:

    約束條件如下:

    式(1) 表示最小化車輛總行駛時間;式(2) 表示每個客戶點僅由一輛車進行訪問;式(3) 表示每個客戶點的流量平衡約束;式(4) 表示每輛車僅使用一次;式(5)、式(6) 表示車輛到達客戶點的時間和客戶點的開始服務時間之間的關系;式(7) 表示一個訂單包含的取貨和送貨請求均由同一輛車提供服務;式(8) 表示車輛從車場出發(fā)完成任務后并返回該車場;式(9) 表示完成客戶點i 的取貨時間早于完成客戶點n+m+i 的送貨時間;式(10) 表示車輛從車場出發(fā)以及返回車場承載貨物量均為0;式(11) 表示車輛在執(zhí)行任務過程中承載貨物容量不能超過車輛容量限制;式(12)、式(13) 表示保證車輛離開上一節(jié)點到下一節(jié)點之間容量的一致性,車輛在客戶點之間承載容量的變化;式(14) 表示決策變量取值約束。

    1.4 時變路網下車輛行駛時間計算方法。在實際交通中,車輛在不同時間段內道路的擁堵情況不同,因而車輛的行駛速度存在差異,考慮實際生活中城市路段早晚高峰的情況,車輛的行駛速度不是階躍變化,而是更近似平滑變化。將一天的工作時間等分為η 個分段,改進Figliozzi的時間依賴函數(shù),使用平滑變化的速度時間函數(shù),如圖1 表示車輛的行駛速度時間函數(shù)。

    圖1 速度時間函數(shù)

    輸入:車輛k 離開節(jié)點i 的出發(fā)時刻λ,時間分段τ,時段范圍(ts,te],節(jié)點i 與節(jié)點j 的距離d。

    1.5 隨機時變車輛行駛時間的魯棒優(yōu)化。隨機時變路網的車輛行駛時間可以通過路段概率分布函數(shù)以及期望估算模型進行優(yōu)化,然而車輛必須在客戶規(guī)定的時間窗內進行服務,完成所有客戶的需求,而不是以一定概率進行滿足。因此,選用最大最小準則作為車輛行駛時間的魯棒優(yōu)化方法,對于車輛k 以時刻λ從節(jié)點i 到節(jié)點j 的最壞情況發(fā)生,從最壞的情況中找到最優(yōu)路徑的行駛時間為:

    2 遺傳模擬退火算法設計

    遺傳算法通過適應度函數(shù)對搜索空間中的多個解進行評估,能夠快速搜索到新解,適用于本文組合優(yōu)化問題的求解,但后期搜索能力減弱致使該算法易過早收斂。本文針對問題特點和遺傳算法的不足,結合模擬退火算法能有效跳出局部循環(huán)的特點,提出了一種高效求解STDPDPTW 模型的混合遺傳模擬退火算法(Hybrid Genetic Simulated Annealing Algorithm, HGSA)。

    2.1 染色體編碼方式。本文設計了一種三行染色體編碼方式,染色體的第一行為車輛訪問的客戶點,第二行為與第一行對應的車輛編號,第三行為客戶點的需求類型,1 為取貨點,-1 為送貨點,0 為車場。每條染色體的編碼方式如圖2 所示,以8 個客戶點為例,車輛1 的路徑為D-8-3-1-5-D,車輛2 的路徑為D-2-4-7-6-D。

    圖2 染色體編碼示意圖

    2.2 初始種群的生成。根據客戶點的需求類型和車輛的容量約束,按照以下5 個步驟生成初始種群:

    Step1:從客戶點集合P 中隨機選擇一個取貨點和一個送貨點進行兩兩組合形成待插入客戶點序列P;

    Step2:根據車輛容量約束,按待插入客戶點序列P順序為車輛劃分客戶點,得到車輛所訪問的客戶點序列H,形成染色體的第一行編碼;

    Step3:在每輛車所訪問的客戶點序列H的最前面和最后面加入該車輛對應的車場編號,用來表示車輛出發(fā)的起點和返回的終點,得到每輛車輛的初始路徑訪問的客戶點序列P,形成染色體的第一行編碼;

    Step4:在每輛車的初始路徑的編碼的第二行添加對應的車輛編號,在第三行添加客戶點對應的任務類型編號,形成一條完整的染色體,如圖3 所示;

    圖3 種群初始化過程示意圖

    Step5:根據預設的種群規(guī)模N重復Step1~Step4,得到初始種群。

    2.3 多段多點交叉算子。多點交叉操作不僅增加染色體的多樣性,而且能夠有效避免個體過早收斂,根據染色體編碼方式的特點,染色體中的每一段代表一輛車訪問的客戶點序列,本文設計了多段多點交叉算子,算法步驟如下:

    Step1:采用輪盤賭方式從種群中選擇兩條染色體作為父代染色體F和F;

    Step3:分別從父代染色體F和F中相同車輛編號的分段F和F中隨機選擇連續(xù)的多個基因進行交換;

    圖4 染色體交叉操作示意圖

    2.4 修復算子。使用多段多點交叉方式后的染色體可能存在部分客戶點缺失和重復的問題,則需要將未訪問的客戶點插入到路徑規(guī)劃方案中,且為縮短車輛的行駛時間,加快全局尋優(yōu)過程,設計了修復算子,算法步驟如下:

    Step1:從交叉操作后的染色體中標記出所有重復訪問的客戶點,利用公式(16) 依次計算染色體中重復點位上一位客戶點n到達該客戶點n的行駛時間與等待時間之和st;

    Step2:重復訪問客戶點中相同編號的客戶點兩兩比較st的大小,保留st較小的客戶點,并去除標記;

    Step3:利用公式(17) 從未訪問客戶點集合N中選擇客戶點依次替換標記出的重復訪問客戶點,生成修復后的新子代個體。

    式(16) 表示染色體中重復點位上一位客戶點n到達該客戶點n的行駛時間與等待時間之和。式(17) 表示從未訪問客戶點集合N中選擇待插入客戶點n,以依次待替換的染色體中重復點位的上一位客戶點n到達未插入的客戶點n之間的行駛時間和等待時間最小化為替換客戶點的選擇準則。

    2.5 變異算子。變異操作有助于提升該算法的全局搜索能力,避免陷入局部最優(yōu)的困境。本文選擇兩點交換變異方式,隨機從染色體中選擇兩個客戶點進行位置交換,得到新的子代染色體,具體操作如圖5 所示。

    圖5 染色體變異操作示意圖

    2.6 種群更新。通過選擇部分父代種群中的精英染色體替代新產生的子代種群中相同數(shù)量的染色體來完成種群更新操作,選擇替代染色體的數(shù)量為:N=N× 1-Ga( )p ,其中N為種群規(guī)模,Gap 為代溝。具體操作如下:一個新的種群由父代種群中適應度值前N位的染色體和子代種群中前(N-N)位的染色體組成,作為遺傳算法下一輪迭代的新種群。

    2.7 鄰域結構。鄰域搜索操作可以擴大搜索范圍增加解的多樣性,避免算法陷入局部最優(yōu),增大算法找到全局最優(yōu)解的可能性。本文采用了四種鄰域搜索算子,包括逆序搜索算子、單點交換算子、合并搜索算子和兩點交換算子,鄰域算子操作如下:(1) 逆序搜索算子:隨機選擇兩個客戶點進行交換,并將兩個客戶點中的客戶點反轉逆序排列,如圖6(a) 所示。(2) 單點交換算子:從原路徑中隨機選擇客戶點移除,插入到其它路徑中,如圖6(b) 所示。(3) 合并搜索算子:隨機選擇一條路徑,將該路徑上的客戶點插入到其它路徑中,如圖6(c) 所示。(4) 兩點交換算子:從兩條不同的路徑中分別各選擇一個客戶點,并交換他們的位置,如圖6(d) 所示。

    圖6 領域算子示意圖

    2.8 混合遺傳模擬退火算法求解步驟?;旌线z傳模擬退火算法將遺傳算法與模擬退火算法結合,通過遺傳算法得到較優(yōu)解,以此作為模擬退火算法的初始解,再通過鄰域搜索最終得到近似最優(yōu)路徑方案。具體算法步驟如下:

    Step1:初始化混合遺傳模擬退火算法的參數(shù),種群規(guī)模N、代溝Gap、交叉概率P、變異概率P、初始溫度T、終止溫度T、步長L、冷卻速率R;

    Step2:種群初始化,獲得車輛初始路徑序列方案集合;

    Step3:通過式(1) 計算染色體的適應度;

    Step4:根據交叉概率P對選擇的兩條父代染色體進行交叉操作生成交叉后的子代染色體;

    Step5:對交叉后的染色體使用修復算子進行修復,補全因交叉過程缺失的客戶點;

    Step6:根據變異概率P對交叉后的子代染色體進行變異操作生成變異后的子代染色體;

    Step7:種群更新操作;

    Step8:判斷種群中最優(yōu)解的適應度值連續(xù)10 代是否下降,若沒有下降,則輸出遺傳算法的當前最優(yōu)解S,轉Step9,否則,轉Step4;

    Step9:初始化當前溫度T=T,當前溫度下迭代次數(shù)L=0,將Step8 輸出的當前最優(yōu)解S 作為當前解S;

    Step10:如果溫度T>T,轉Step11,否則,輸出算法終止輸出最優(yōu)解;

    Step11:如果迭代次數(shù)L<L,則L=L+1,轉12,否則,L=0,T=T×R,轉Step10;

    3 實驗分析

    本文所有實驗使用Java 語言進行編寫,eclipse2018-12 的編程環(huán)境,Windows 10 操作系統(tǒng),Intel Core i7-7700 CPU@3.60 GHz 16.0 GB 的運行環(huán)境。

    3.1 參數(shù)設置。通過合理的參數(shù)設置能夠有效提高算法的求解質量,考慮到實際應用需求,選取多個不同算例,在合理的時間范圍內,采用控制變量法,最終確定一個最佳的參數(shù)組合。首先根據經驗確定參數(shù)的測試區(qū)間,再通過更改參數(shù)組合進行多組數(shù)值實驗,記錄實驗結果和CPU 運行時間,最終確定HGSA 參數(shù)設置最佳組合如表1 所示。

    表1 算法的最佳參數(shù)值

    3.2 PDPTW 算例。目前關于STDPDPTW 研究所進行的實驗并沒有公認的數(shù)據集,且STDPDPTW 是在經典的PDPTW 問題中加部分約束形成的衍生問題。因此,選取Li & Lim PDPTW 標準數(shù)據集pdp_100的56 個算例進行數(shù)值實驗,驗證本文算法有效性。

    本實驗選用Li & Lim PDPTW 標準數(shù)據集中的56 個算例實驗,算例包括LC 數(shù)據集(地點規(guī)則分布)、LR 數(shù)據集(地隨機分布) 和LRC 數(shù)據集(地點混合分布)。所有車輛行駛道路類型相同,且以速度為1 的恒定速度行駛,以車輛總行駛距為優(yōu)化目標。表2 記錄了PDPTW 的已知最好解R、HGSA 算法的總距離S及其運算時間,Gap 表示R與S之間的差百分比,Gap=(R-S)/S。增的 點離值

    表2 PDPTW 數(shù)據集實驗結果

    由表2 中數(shù)據可知,HGSA 算法更新了5 個PDPTW 算例分別為lc109、lc203、lr108、lrc106 和lrc204 的已知最好解。在其它的51 個PDPTW 算例中,HGSA 算法的求解結果劣于算例已知最好解,但與已知最好解的差距均值在1%之內,實驗數(shù)據說明了HGSA 算法求解PDPTW 的有效性。

    3.3 STDPDPTW 算例。基于Li & Lim 基準測試實例,根據實例規(guī)模將單一車場擴充為5 個車場,每個車場位置隨機生成,且每個車場分配5 輛車用于完成任務,并增加路段中車輛行駛速度的波動范圍,構建STDPDPTW 測試實例。依據城市道路通行情況,將道路分為五種類型,總服務時域分為9 個等分的時間分段,依據圖1 將車輛的行駛速度時間函數(shù)設置為連續(xù)函數(shù),不同的時段車輛的最高行駛速度和最低行駛速度不同,路段的道路類型由式(18) 判斷:

    式(18) 中,mod 為取余計算,例如從7 號節(jié)點到9 號節(jié)點的道路類型為2。

    表3 表示不同道路類型在9 個等分的時間分段中車輛行駛速度的變化范圍。對各道路類型和時間分段設置高、中、低三種波動范圍,對于時間分段τ、τ和τ中行駛速度增加±15%的波動范圍,對于時間分段τ、τ和τ中行駛速度增加±20%的波動范圍,對于時間分段τ、τ和τ中行駛速度增加±10%的波動范圍。根據行駛速度波動范圍通過隨機時變車輛行駛時間的魯棒優(yōu)化方法將隨機時變路網轉化為確定型路網。

    表3 各道路類型車輛行駛速度變化范圍

    表4 記錄了隨機選取的10 個STDPDPTW 算例的實驗結果,在相同實驗條件下每個算例運行20 次。其中,20 次運行結果的最大值記錄在C列中,平均值記錄在C列中,最小值記錄在C列中,20 次運行的平均運行時間記錄在CPU列中,Gap表示C與C之間的差值百分比,Gap表示C與C之間的差值百分比。其中,Gap=(C-C)/C,Gap=(C-C)/C。

    表4 STDPDPTW 數(shù)據集實驗結果

    由表4 數(shù)據可知,考慮車輛行駛時間的時變性和隨機性,在車輛數(shù)量有限的情況下,通過優(yōu)化車輛的出發(fā)時間和出發(fā)的車場位置可以縮短完成配送任務的時間。10 個STPDPTW 算例均在較短的時間內獲得近似最優(yōu)解,將引入多段多點交叉算子和修復算子的遺傳算法與模擬退火算法結合,有效提高了局部尋優(yōu)能力和收斂速度。另外,表4 的C與C之間的差值百分比的平均值僅相差2.01%,C與C之間的差值百分比的平均值僅相差1.68%,表明HGSA 具有較強的穩(wěn)定性和尋優(yōu)能力。

    4 結 論

    本文考慮城市交通擁堵、交通管制等不確定因素的影響,建立了隨機時變下帶時間窗的取送貨車輛路徑問題模型,提出考慮車輛加速度的行駛時間計算方式,運用隨機時變車輛行駛時間的魯棒優(yōu)化方法將隨機性時變轉化為確定型時變,根據STDPDPTW 模型特點,設計了兩階段的混合遺傳模擬退火算法,引入多段多點交叉算子和修復算子,保證了前期迭代的收斂速度后使用多領域搜索算子,提高了算法的搜索效率。通過標準數(shù)據集進行驗證,HGSA 算法更新了56 個PDPTW 算例中5 個已知最好解,其余算例與已知最好解的差距均值在1%之內,同時STDPDPTW 算例表明,HGSA 算法求解STDPDPTW 具有較強的穩(wěn)定性,能有效縮短車輛的總行駛時間,降低物流配送成本。

    猜你喜歡
    車場模擬退火時變
    結合模擬退火和多分配策略的密度峰值聚類算法
    時變干擾下AUV三維軌跡跟蹤反步滑??刂?/a>
    城市軌道交通車場乘降所信號設計方案研究
    模擬退火遺傳算法在機械臂路徑規(guī)劃中的應用
    測控技術(2018年3期)2018-11-25 09:45:08
    基于神經網絡的高速鐵路動車存車場火災識別算法研究
    電子測試(2018年11期)2018-06-26 05:56:10
    鐵路客車存車場火災自動報警系統(tǒng)設計
    基于時變Copula的股票市場相關性分析
    智富時代(2017年4期)2017-04-27 17:08:47
    基于模糊自適應模擬退火遺傳算法的配電網故障定位
    鈾礦山井底車場巷道內氡及其子體濃度分布規(guī)律研究
    SOA結合模擬退火算法優(yōu)化電容器配置研究
    電源技術(2015年5期)2015-08-22 11:18:24
    联通29元200g的流量卡| 丰满少妇做爰视频| 精品99又大又爽又粗少妇毛片| 真实男女啪啪啪动态图| 欧美区成人在线视频| 最新中文字幕久久久久| 久久精品久久精品一区二区三区| 成人亚洲精品av一区二区| 国产伦在线观看视频一区| 天天躁夜夜躁狠狠久久av| 欧美精品国产亚洲| 噜噜噜噜噜久久久久久91| 青春草亚洲视频在线观看| 国产精品一及| 久久精品影院6| 久久久精品94久久精品| 国产视频首页在线观看| 长腿黑丝高跟| 男女边吃奶边做爰视频| 国产国拍精品亚洲av在线观看| 国产高清视频在线观看网站| 亚洲国产欧美人成| 青春草亚洲视频在线观看| 亚洲国产精品成人综合色| 久久6这里有精品| 国产成人a∨麻豆精品| 日韩亚洲欧美综合| 秋霞在线观看毛片| 别揉我奶头 嗯啊视频| 亚洲三级黄色毛片| 特级一级黄色大片| 看免费成人av毛片| 国产一级毛片七仙女欲春2| 婷婷六月久久综合丁香| 好男人视频免费观看在线| 女的被弄到高潮叫床怎么办| 熟妇人妻久久中文字幕3abv| 免费观看a级毛片全部| 欧美一区二区亚洲| 久久精品国产亚洲网站| 高清视频免费观看一区二区 | 夜夜看夜夜爽夜夜摸| 亚洲美女视频黄频| 毛片一级片免费看久久久久| 久久99热6这里只有精品| 特大巨黑吊av在线直播| 乱码一卡2卡4卡精品| 亚洲性久久影院| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品久久久久久婷婷小说 | 国产老妇伦熟女老妇高清| 啦啦啦啦在线视频资源| 亚洲av电影不卡..在线观看| 亚洲国产精品久久男人天堂| 十八禁国产超污无遮挡网站| 国产成人精品一,二区| 尾随美女入室| 少妇高潮的动态图| 精品人妻偷拍中文字幕| 美女国产视频在线观看| 免费看美女性在线毛片视频| 亚洲av一区综合| 哪个播放器可以免费观看大片| av免费在线看不卡| 视频中文字幕在线观看| 久久综合国产亚洲精品| 国产视频内射| 久久久欧美国产精品| 国产爱豆传媒在线观看| 国产欧美另类精品又又久久亚洲欧美| av在线播放精品| 能在线免费观看的黄片| 欧美一级a爱片免费观看看| 国内精品一区二区在线观看| 99久国产av精品| 亚洲天堂国产精品一区在线| 欧美激情国产日韩精品一区| 久久精品人妻少妇| 国产伦精品一区二区三区视频9| 高清av免费在线| 毛片一级片免费看久久久久| 老司机福利观看| 99久久精品一区二区三区| 国产成人freesex在线| 日韩人妻高清精品专区| 五月玫瑰六月丁香| 丰满少妇做爰视频| av线在线观看网站| 免费无遮挡裸体视频| 国产精品无大码| 成人漫画全彩无遮挡| 国产老妇伦熟女老妇高清| 午夜福利成人在线免费观看| 内地一区二区视频在线| 亚洲国产日韩欧美精品在线观看| 亚洲第一区二区三区不卡| av.在线天堂| 亚洲人成网站在线播| 精品久久久久久久人妻蜜臀av| 国内少妇人妻偷人精品xxx网站| 亚洲精品色激情综合| 丰满少妇做爰视频| 国产色爽女视频免费观看| 国产av码专区亚洲av| 26uuu在线亚洲综合色| 99九九线精品视频在线观看视频| 好男人在线观看高清免费视频| 精品人妻一区二区三区麻豆| 亚洲av二区三区四区| 日韩精品有码人妻一区| 视频中文字幕在线观看| 韩国av在线不卡| 赤兔流量卡办理| 日韩av不卡免费在线播放| av又黄又爽大尺度在线免费看 | 女人被狂操c到高潮| 免费一级毛片在线播放高清视频| 国产精品一区二区三区四区免费观看| 啦啦啦韩国在线观看视频| 亚洲丝袜综合中文字幕| 小蜜桃在线观看免费完整版高清| 黄色一级大片看看| 成人亚洲欧美一区二区av| 边亲边吃奶的免费视频| 最近最新中文字幕免费大全7| 美女xxoo啪啪120秒动态图| 亚洲国产欧美在线一区| 久久这里只有精品中国| 一区二区三区四区激情视频| 亚洲国产成人一精品久久久| 特级一级黄色大片| 少妇丰满av| 亚洲av熟女| 美女高潮的动态| 中国美白少妇内射xxxbb| 国产极品精品免费视频能看的| 精品国产三级普通话版| 久久99热这里只频精品6学生 | 国产视频内射| 国产 一区 欧美 日韩| 国产亚洲精品av在线| 婷婷色麻豆天堂久久 | 大又大粗又爽又黄少妇毛片口| 精品国产一区二区三区久久久樱花 | 一级毛片aaaaaa免费看小| 欧美日本视频| 免费观看在线日韩| 久久这里只有精品中国| 99国产精品一区二区蜜桃av| 在现免费观看毛片| 欧美另类亚洲清纯唯美| 亚洲丝袜综合中文字幕| 网址你懂的国产日韩在线| 1000部很黄的大片| 亚洲精品成人久久久久久| 狂野欧美激情性xxxx在线观看| 最近中文字幕高清免费大全6| 成人性生交大片免费视频hd| 国产综合懂色| 成人高潮视频无遮挡免费网站| 神马国产精品三级电影在线观看| 最近最新中文字幕免费大全7| 日本av手机在线免费观看| 国产精品一区www在线观看| 观看免费一级毛片| 在线天堂最新版资源| 久久99热这里只有精品18| 欧美色视频一区免费| 久久久久免费精品人妻一区二区| 99热网站在线观看| 免费看av在线观看网站| 麻豆久久精品国产亚洲av| 最近视频中文字幕2019在线8| 内地一区二区视频在线| 汤姆久久久久久久影院中文字幕 | 2021少妇久久久久久久久久久| 午夜视频国产福利| 亚洲一级一片aⅴ在线观看| 51国产日韩欧美| 日韩,欧美,国产一区二区三区 | 亚洲欧美成人精品一区二区| 免费电影在线观看免费观看| 久久精品国产亚洲网站| 亚洲av日韩在线播放| 综合色丁香网| 啦啦啦韩国在线观看视频| 亚洲精品日韩av片在线观看| 国产亚洲av嫩草精品影院| 国产午夜精品一二区理论片| 国产成人91sexporn| 午夜激情福利司机影院| 国产午夜福利久久久久久| 国产精品女同一区二区软件| 久久久久九九精品影院| 国产淫语在线视频| 少妇被粗大猛烈的视频| 国产精品一区www在线观看| 九草在线视频观看| 变态另类丝袜制服| 亚洲国产精品国产精品| 国产 一区精品| 国产精品嫩草影院av在线观看| 搞女人的毛片| 亚洲欧洲日产国产| 午夜久久久久精精品| 午夜日本视频在线| 国产亚洲精品久久久com| 日日摸夜夜添夜夜添av毛片| www.色视频.com| 欧美xxxx性猛交bbbb| 日本熟妇午夜| 国产av不卡久久| 男人狂女人下面高潮的视频| 国产精品99久久久久久久久| 看非洲黑人一级黄片| 亚洲av熟女| 一卡2卡三卡四卡精品乱码亚洲| www.av在线官网国产| 少妇人妻一区二区三区视频| 狠狠狠狠99中文字幕| 久久久久性生活片| 啦啦啦观看免费观看视频高清| 国产精品久久久久久久电影| 亚洲国产精品sss在线观看| 国产av不卡久久| 丰满少妇做爰视频| 精品99又大又爽又粗少妇毛片| 亚洲精华国产精华液的使用体验| 蜜臀久久99精品久久宅男| 欧美性猛交黑人性爽| 日本一二三区视频观看| 国产精品麻豆人妻色哟哟久久 | 身体一侧抽搐| 亚洲内射少妇av| 内地一区二区视频在线| 国产极品天堂在线| 亚洲国产色片| 亚洲成色77777| 综合色av麻豆| 亚洲三级黄色毛片| 亚洲怡红院男人天堂| 亚洲人与动物交配视频| 天堂中文最新版在线下载 | 最新中文字幕久久久久| 99热6这里只有精品| 在线a可以看的网站| 深爱激情五月婷婷| 午夜免费男女啪啪视频观看| 床上黄色一级片| 亚洲精品乱码久久久v下载方式| 在线观看av片永久免费下载| 日韩亚洲欧美综合| 一级黄色大片毛片| 岛国在线免费视频观看| 欧美激情国产日韩精品一区| 青春草亚洲视频在线观看| 亚洲在线观看片| 日韩中字成人| 亚洲精品456在线播放app| 国产高潮美女av| 亚洲国产最新在线播放| 国产白丝娇喘喷水9色精品| 亚洲最大成人av| 最近最新中文字幕大全电影3| 亚洲欧美日韩无卡精品| 美女内射精品一级片tv| 精品一区二区三区视频在线| 三级国产精品片| 人妻制服诱惑在线中文字幕| 亚洲美女视频黄频| 日本一二三区视频观看| 天堂影院成人在线观看| 久久久久久久久久成人| 国产真实伦视频高清在线观看| 91久久精品国产一区二区三区| 久99久视频精品免费| 啦啦啦韩国在线观看视频| 自拍偷自拍亚洲精品老妇| 一个人免费在线观看电影| 久久草成人影院| 国产免费男女视频| 神马国产精品三级电影在线观看| 欧美一区二区亚洲| 国产91av在线免费观看| 一个人看视频在线观看www免费| 日韩中字成人| 欧美成人一区二区免费高清观看| 又黄又爽又刺激的免费视频.| 91午夜精品亚洲一区二区三区| 免费观看性生交大片5| 十八禁国产超污无遮挡网站| 亚洲真实伦在线观看| 91精品伊人久久大香线蕉| 国产视频首页在线观看| 成人二区视频| 日韩,欧美,国产一区二区三区 | 亚洲va在线va天堂va国产| 97超碰精品成人国产| 我的老师免费观看完整版| 国产视频首页在线观看| 久久久久久久久久久丰满| 寂寞人妻少妇视频99o| 少妇裸体淫交视频免费看高清| 亚洲三级黄色毛片| 蜜臀久久99精品久久宅男| 国语对白做爰xxxⅹ性视频网站| 中国美白少妇内射xxxbb| 国内揄拍国产精品人妻在线| 免费观看的影片在线观看| 一级爰片在线观看| 日韩欧美精品v在线| 1024手机看黄色片| 搞女人的毛片| 精品午夜福利在线看| 中文字幕免费在线视频6| 免费观看在线日韩| 精品免费久久久久久久清纯| 日本一本二区三区精品| 亚洲成人精品中文字幕电影| 国产免费视频播放在线视频 | 免费一级毛片在线播放高清视频| 日韩中字成人| 熟女人妻精品中文字幕| 国产精品国产高清国产av| 久久精品影院6| 国产精品永久免费网站| 亚洲国产精品sss在线观看| 亚洲五月天丁香| 在线播放国产精品三级| 亚洲欧美精品专区久久| 两个人视频免费观看高清| 精品久久久久久成人av| 日产精品乱码卡一卡2卡三| 1024手机看黄色片| 国产老妇女一区| 舔av片在线| 国产日韩欧美在线精品| 男女下面进入的视频免费午夜| 最近2019中文字幕mv第一页| 亚洲成色77777| 久久久久久久午夜电影| 男人狂女人下面高潮的视频| 国产免费福利视频在线观看| 国产免费视频播放在线视频 | 久久久欧美国产精品| 99久久中文字幕三级久久日本| 麻豆乱淫一区二区| 夜夜看夜夜爽夜夜摸| 色播亚洲综合网| 午夜免费男女啪啪视频观看| 我的女老师完整版在线观看| 嘟嘟电影网在线观看| 日本免费一区二区三区高清不卡| 亚洲高清免费不卡视频| 赤兔流量卡办理| 亚洲av成人av| 欧美成人精品欧美一级黄| ponron亚洲| 嫩草影院入口| 亚洲欧美日韩东京热| 高清av免费在线| 最近最新中文字幕大全电影3| 国产一区二区在线观看日韩| 麻豆精品久久久久久蜜桃| 中文天堂在线官网| 99久国产av精品国产电影| 国产毛片a区久久久久| 国产精品日韩av在线免费观看| 十八禁国产超污无遮挡网站| 亚洲欧美成人综合另类久久久 | 亚洲精品乱码久久久久久按摩| 国产高清有码在线观看视频| 亚洲成av人片在线播放无| 天堂√8在线中文| 十八禁国产超污无遮挡网站| 最新中文字幕久久久久| 色视频www国产| 22中文网久久字幕| 亚洲精品国产av成人精品| 女人久久www免费人成看片 | 特级一级黄色大片| 色综合站精品国产| 九九爱精品视频在线观看| 在线观看av片永久免费下载| 日韩欧美精品免费久久| 黄色一级大片看看| 亚洲av电影在线观看一区二区三区 | 中文欧美无线码| 精品欧美国产一区二区三| 国产高清不卡午夜福利| 久久6这里有精品| 中文资源天堂在线| 淫秽高清视频在线观看| 午夜久久久久精精品| 五月伊人婷婷丁香| 非洲黑人性xxxx精品又粗又长| 欧美区成人在线视频| av在线播放精品| 精品久久久噜噜| eeuss影院久久| 久久久国产成人免费| 欧美人与善性xxx| 亚洲av.av天堂| 久久久国产成人精品二区| 黄色日韩在线| 男女国产视频网站| av黄色大香蕉| 午夜福利在线在线| 欧美色视频一区免费| 一个人看视频在线观看www免费| 日韩大片免费观看网站 | 两性午夜刺激爽爽歪歪视频在线观看| 在线免费观看的www视频| 国产人妻一区二区三区在| 国产高清不卡午夜福利| 久久99热这里只频精品6学生 | 成人三级黄色视频| 亚洲一区高清亚洲精品| 亚洲精品久久久久久婷婷小说 | 搞女人的毛片| 好男人视频免费观看在线| 日本黄色视频三级网站网址| 熟女人妻精品中文字幕| 亚洲欧美日韩无卡精品| 亚洲精品一区蜜桃| 长腿黑丝高跟| av在线亚洲专区| 欧美zozozo另类| 插逼视频在线观看| 亚洲人与动物交配视频| 亚洲精品日韩av片在线观看| 又爽又黄a免费视频| 亚洲国产欧美在线一区| 亚洲欧美成人精品一区二区| 91久久精品国产一区二区成人| 网址你懂的国产日韩在线| 亚洲综合精品二区| 欧美精品一区二区大全| a级毛色黄片| 国产一区二区在线观看日韩| av在线老鸭窝| 久久久久久大精品| 国产精品乱码一区二三区的特点| 国语对白做爰xxxⅹ性视频网站| 成人毛片60女人毛片免费| 秋霞在线观看毛片| 免费播放大片免费观看视频在线观看 | 成年版毛片免费区| 亚洲最大成人手机在线| 欧美精品一区二区大全| 男人狂女人下面高潮的视频| 亚洲av男天堂| 麻豆一二三区av精品| 免费大片18禁| 国产高清不卡午夜福利| 熟女人妻精品中文字幕| 精华霜和精华液先用哪个| 舔av片在线| 内射极品少妇av片p| 亚洲婷婷狠狠爱综合网| 1000部很黄的大片| 免费看av在线观看网站| 国产成人免费观看mmmm| 可以在线观看毛片的网站| 国产精品女同一区二区软件| 国产高清不卡午夜福利| 熟女人妻精品中文字幕| 欧美日韩精品成人综合77777| 久久久久久久久久久丰满| 长腿黑丝高跟| 国产伦在线观看视频一区| 国产免费视频播放在线视频 | 国产高清视频在线观看网站| 中国美白少妇内射xxxbb| 久久韩国三级中文字幕| 麻豆久久精品国产亚洲av| 麻豆精品久久久久久蜜桃| 热99在线观看视频| 中文字幕久久专区| 天美传媒精品一区二区| 精品国产一区二区三区久久久樱花 | 可以在线观看毛片的网站| 一级毛片久久久久久久久女| 国产高清不卡午夜福利| 三级毛片av免费| 亚洲成人中文字幕在线播放| 一级爰片在线观看| 亚洲av一区综合| 在线a可以看的网站| 亚洲成人久久爱视频| 亚洲国产精品sss在线观看| 级片在线观看| 亚洲欧洲日产国产| 色综合色国产| 2021少妇久久久久久久久久久| 亚洲精品aⅴ在线观看| 国产美女午夜福利| 色噜噜av男人的天堂激情| 中文字幕亚洲精品专区| 日韩中字成人| 欧美另类亚洲清纯唯美| 草草在线视频免费看| 免费黄网站久久成人精品| 久久精品国产自在天天线| 高清毛片免费看| 99久久精品热视频| 国产久久久一区二区三区| 成人三级黄色视频| 干丝袜人妻中文字幕| 欧美日本亚洲视频在线播放| 女人被狂操c到高潮| 精品人妻视频免费看| 成人鲁丝片一二三区免费| 亚洲自拍偷在线| 久久精品国产99精品国产亚洲性色| 亚洲精品,欧美精品| 欧美丝袜亚洲另类| 国产爱豆传媒在线观看| 亚洲成人久久爱视频| 久久久亚洲精品成人影院| 51国产日韩欧美| 国产精品国产三级专区第一集| 一级黄片播放器| 黄片无遮挡物在线观看| 久久草成人影院| 久久精品国产自在天天线| 嘟嘟电影网在线观看| 欧美日本亚洲视频在线播放| 建设人人有责人人尽责人人享有的 | 别揉我奶头 嗯啊视频| 久久精品影院6| 欧美激情国产日韩精品一区| 最近手机中文字幕大全| 午夜日本视频在线| 别揉我奶头 嗯啊视频| 精品久久久久久久久av| 熟妇人妻久久中文字幕3abv| 男人的好看免费观看在线视频| 日韩人妻高清精品专区| 成年av动漫网址| 天堂√8在线中文| 不卡视频在线观看欧美| 欧美xxxx黑人xx丫x性爽| 欧美一区二区精品小视频在线| 国产精品伦人一区二区| 国产私拍福利视频在线观看| 在线播放无遮挡| 性插视频无遮挡在线免费观看| 我的老师免费观看完整版| 一级av片app| 国产精品,欧美在线| 人人妻人人看人人澡| 嫩草影院入口| 日韩av在线免费看完整版不卡| 中文字幕制服av| 又粗又硬又长又爽又黄的视频| 欧美人与善性xxx| 久久国内精品自在自线图片| 黄片无遮挡物在线观看| 国产色爽女视频免费观看| 一级二级三级毛片免费看| 精品国产露脸久久av麻豆 | 最近2019中文字幕mv第一页| 欧美bdsm另类| 国产精品电影一区二区三区| 欧美激情在线99| 白带黄色成豆腐渣| 亚洲精品国产av成人精品| 国产午夜福利久久久久久| 国产片特级美女逼逼视频| 国产午夜福利久久久久久| 亚洲欧美日韩高清专用| 免费观看在线日韩| 亚洲最大成人手机在线| 两个人的视频大全免费| 精品人妻熟女av久视频| 日日啪夜夜撸| 干丝袜人妻中文字幕| eeuss影院久久| 亚洲国产精品国产精品| 国产欧美另类精品又又久久亚洲欧美| 大又大粗又爽又黄少妇毛片口| 欧美一区二区国产精品久久精品| 在线播放无遮挡| 欧美激情在线99| 欧美区成人在线视频| 国产老妇伦熟女老妇高清| 91午夜精品亚洲一区二区三区| 国产精品一二三区在线看| 深夜a级毛片| 亚洲国产精品久久男人天堂| 在线免费观看的www视频| 大又大粗又爽又黄少妇毛片口| 男女边吃奶边做爰视频| 18+在线观看网站| 日本免费a在线| 黄色一级大片看看| 久久久久九九精品影院| 91久久精品电影网| 一级黄片播放器| 久久精品国产自在天天线| 别揉我奶头 嗯啊视频| 国产成年人精品一区二区| 秋霞伦理黄片| 我要搜黄色片| 69av精品久久久久久| 一个人看视频在线观看www免费| 一边摸一边抽搐一进一小说| 欧美一区二区精品小视频在线| 久久久久性生活片| www日本黄色视频网| 日韩欧美国产在线观看| 亚洲精品色激情综合| 国产成人a区在线观看| 国产一区二区亚洲精品在线观看|