• <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
    久久久久免费精品人妻一区二区| 国产麻豆成人av免费视频| 国产精品久久久久久精品电影| 欧美乱色亚洲激情| 久久精品国产99精品国产亚洲性色| 一个人免费在线观看电影| 最好的美女福利视频网| 在线播放无遮挡| 亚洲av不卡在线观看| 最近最新免费中文字幕在线| 91麻豆精品激情在线观看国产| 午夜两性在线视频| 啦啦啦观看免费观看视频高清| 国产淫片久久久久久久久 | 九九久久精品国产亚洲av麻豆| 亚洲第一区二区三区不卡| 中文亚洲av片在线观看爽| 老熟妇乱子伦视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 男女床上黄色一级片免费看| 国内精品美女久久久久久| 国产激情偷乱视频一区二区| 国产成人欧美在线观看| 毛片一级片免费看久久久久 | 色尼玛亚洲综合影院| 青草久久国产| 丰满乱子伦码专区| 少妇被粗大猛烈的视频| 两个人的视频大全免费| 国产在线男女| 成年女人毛片免费观看观看9| 天美传媒精品一区二区| av欧美777| 欧美三级亚洲精品| 中文字幕久久专区| 老司机午夜福利在线观看视频| av在线天堂中文字幕| av国产免费在线观看| 午夜久久久久精精品| 午夜福利在线观看吧| 深夜a级毛片| 亚洲 欧美 日韩 在线 免费| 99精品在免费线老司机午夜| 国产亚洲欧美98| 亚洲av免费高清在线观看| 成人特级av手机在线观看| 国产不卡一卡二| 午夜两性在线视频| 中文资源天堂在线| 十八禁国产超污无遮挡网站| 亚洲人成伊人成综合网2020| 一本久久中文字幕| 欧美性感艳星| 九色成人免费人妻av| 久久精品国产自在天天线| 老熟妇仑乱视频hdxx| 婷婷亚洲欧美| 国产色婷婷99| 亚洲avbb在线观看| 99热6这里只有精品| 亚洲精品影视一区二区三区av| 免费av毛片视频| x7x7x7水蜜桃| 超碰av人人做人人爽久久| 九九久久精品国产亚洲av麻豆| 免费搜索国产男女视频| 老熟妇仑乱视频hdxx| 国产白丝娇喘喷水9色精品| 成人欧美大片| 国产精品久久电影中文字幕| 午夜精品久久久久久毛片777| av视频在线观看入口| 亚洲av日韩精品久久久久久密| 美女黄网站色视频| 神马国产精品三级电影在线观看| 成人av一区二区三区在线看| 国产精品一区二区三区四区免费观看 | 国内精品久久久久久久电影| 日日干狠狠操夜夜爽| 国产黄a三级三级三级人| 大型黄色视频在线免费观看| 久久精品国产亚洲av涩爱 | 很黄的视频免费| 中国美女看黄片| 精品免费久久久久久久清纯| 国产精品三级大全| 97超视频在线观看视频| av欧美777| 亚洲国产精品sss在线观看| 午夜福利高清视频| 黄色日韩在线| 波野结衣二区三区在线| 色综合婷婷激情| 国产又黄又爽又无遮挡在线| 国产v大片淫在线免费观看| 亚洲成人久久爱视频| 男人狂女人下面高潮的视频| 麻豆av噜噜一区二区三区| 俺也久久电影网| 91久久精品国产一区二区成人| 久久久久久久午夜电影| 香蕉av资源在线| 亚洲欧美日韩无卡精品| 亚洲 欧美 日韩 在线 免费| 国产黄片美女视频| 欧美日韩国产亚洲二区| 99久久99久久久精品蜜桃| 亚洲av成人精品一区久久| 亚洲国产精品久久男人天堂| 非洲黑人性xxxx精品又粗又长| 性色avwww在线观看| 精品免费久久久久久久清纯| 欧美日韩乱码在线| 国产精品自产拍在线观看55亚洲| 国产精品,欧美在线| 婷婷色综合大香蕉| 色哟哟哟哟哟哟| 日本黄色视频三级网站网址| 亚洲经典国产精华液单 | 亚州av有码| 久久精品夜夜夜夜夜久久蜜豆| 午夜免费成人在线视频| 悠悠久久av| 亚洲中文字幕一区二区三区有码在线看| 久久国产精品影院| 亚洲最大成人中文| 又黄又爽又刺激的免费视频.| 亚洲一区高清亚洲精品| 午夜精品在线福利| 久久久久久久精品吃奶| 国产一区二区激情短视频| 99热只有精品国产| 热99在线观看视频| 国产精品久久久久久精品电影| 国产一区二区在线av高清观看| 日韩高清综合在线| 欧美乱妇无乱码| 成年免费大片在线观看| 国产成人欧美在线观看| 欧美成人一区二区免费高清观看| 天美传媒精品一区二区| 少妇裸体淫交视频免费看高清| 久久精品人妻少妇| 男女那种视频在线观看| 婷婷丁香在线五月| 97超级碰碰碰精品色视频在线观看| 久久久色成人| 毛片女人毛片| av在线观看视频网站免费| 日韩人妻高清精品专区| 人妻夜夜爽99麻豆av| 人妻夜夜爽99麻豆av| 久久久久久九九精品二区国产| 此物有八面人人有两片| 免费看a级黄色片| 亚洲美女黄片视频| 日本五十路高清| 午夜福利在线在线| 少妇人妻一区二区三区视频| 亚洲精品日韩av片在线观看| 午夜福利18| 久久精品国产自在天天线| 国产精品av视频在线免费观看| 日本 av在线| 久久亚洲精品不卡| 丁香欧美五月| av欧美777| 日韩欧美一区二区三区在线观看| 久久久久久久久大av| 美女高潮喷水抽搐中文字幕| 91麻豆精品激情在线观看国产| 欧美精品国产亚洲| 日韩欧美三级三区| 国产国拍精品亚洲av在线观看| 淫秽高清视频在线观看| 日本免费a在线| 亚洲国产精品成人综合色| 亚洲五月天丁香| 日本 欧美在线| 免费看日本二区| 精品国内亚洲2022精品成人| 欧美日韩黄片免| 神马国产精品三级电影在线观看| 亚洲自拍偷在线| 精华霜和精华液先用哪个| 亚洲精品日韩av片在线观看| 日本 av在线| 丰满人妻一区二区三区视频av| 九色成人免费人妻av| 90打野战视频偷拍视频| 亚洲七黄色美女视频| 69人妻影院| 国产中年淑女户外野战色| 人妻久久中文字幕网| 精品一区二区三区av网在线观看| 国产激情偷乱视频一区二区| 欧美日韩瑟瑟在线播放| 成人国产综合亚洲| 成年女人看的毛片在线观看| 精品人妻偷拍中文字幕| 日韩精品青青久久久久久| 美女大奶头视频| 波多野结衣高清作品| 偷拍熟女少妇极品色| 少妇裸体淫交视频免费看高清| 女人十人毛片免费观看3o分钟| 男插女下体视频免费在线播放| 好男人在线观看高清免费视频| 十八禁人妻一区二区| 欧美色视频一区免费| 97超级碰碰碰精品色视频在线观看| 国产69精品久久久久777片| 国产精品久久电影中文字幕| 真实男女啪啪啪动态图| 淫妇啪啪啪对白视频| 精品一区二区三区av网在线观看| 成人精品一区二区免费| 特级一级黄色大片| 国产精品1区2区在线观看.| 12—13女人毛片做爰片一| 99国产极品粉嫩在线观看| 中文字幕人成人乱码亚洲影| 亚洲,欧美精品.| a级毛片免费高清观看在线播放| 深夜精品福利| 九色成人免费人妻av| 国产黄a三级三级三级人| 午夜亚洲福利在线播放| 色精品久久人妻99蜜桃| 精品熟女少妇八av免费久了| www.熟女人妻精品国产| 免费一级毛片在线播放高清视频| 色5月婷婷丁香| 婷婷精品国产亚洲av| 三级男女做爰猛烈吃奶摸视频| 在线天堂最新版资源| 国产乱人视频| 亚洲美女视频黄频| 国产精品日韩av在线免费观看| 国内毛片毛片毛片毛片毛片| 国产亚洲av嫩草精品影院| 在线看三级毛片| 亚洲久久久久久中文字幕| 午夜日韩欧美国产| 国产白丝娇喘喷水9色精品| 日韩精品青青久久久久久| 男女做爰动态图高潮gif福利片| 一个人免费在线观看的高清视频| 亚洲欧美日韩卡通动漫| a在线观看视频网站| x7x7x7水蜜桃| 亚洲在线自拍视频| 国产午夜精品论理片| 亚洲 欧美 日韩 在线 免费| 少妇人妻一区二区三区视频| 国内揄拍国产精品人妻在线| 精品午夜福利在线看| av在线天堂中文字幕| 一本久久中文字幕| 特大巨黑吊av在线直播| 麻豆一二三区av精品| 黄色一级大片看看| 亚洲成人久久性| 99热只有精品国产| 欧美激情在线99| 免费在线观看成人毛片| 国内久久婷婷六月综合欲色啪| 中文资源天堂在线| 悠悠久久av| 亚洲片人在线观看| 床上黄色一级片| 听说在线观看完整版免费高清| 国产精华一区二区三区| 国产成人欧美在线观看| 天堂av国产一区二区熟女人妻| 免费av不卡在线播放| 好男人电影高清在线观看| 久久精品国产清高在天天线| 嫩草影院新地址| 国产中年淑女户外野战色| 久久久精品大字幕| 久久精品国产亚洲av涩爱 | 九九在线视频观看精品| 嫩草影视91久久| 天堂动漫精品| 一个人免费在线观看电影| av女优亚洲男人天堂| 精品久久久久久久久久久久久| 内射极品少妇av片p| 国产午夜精品久久久久久一区二区三区 | 国产欧美日韩一区二区三| 精品免费久久久久久久清纯| 别揉我奶头~嗯~啊~动态视频| 久久国产精品人妻蜜桃| 一个人看视频在线观看www免费| 国产精品电影一区二区三区| 亚洲人与动物交配视频| 在线观看一区二区三区| 老司机福利观看| 免费看日本二区| 欧美日韩福利视频一区二区| 国产精品爽爽va在线观看网站| 一级毛片久久久久久久久女| 国产亚洲av嫩草精品影院| 成人亚洲精品av一区二区| 欧美乱色亚洲激情| 欧美又色又爽又黄视频| 长腿黑丝高跟| 欧美乱妇无乱码| av在线老鸭窝| 熟女电影av网| 日韩中字成人| а√天堂www在线а√下载| 精品无人区乱码1区二区| 亚洲在线自拍视频| 小说图片视频综合网站| 久久精品夜夜夜夜夜久久蜜豆| 在线观看美女被高潮喷水网站 | 国内久久婷婷六月综合欲色啪| 亚洲国产日韩欧美精品在线观看| 亚洲av五月六月丁香网| 欧美日本亚洲视频在线播放| 亚洲18禁久久av| 男插女下体视频免费在线播放| 欧美另类亚洲清纯唯美| 特级一级黄色大片| 欧美一区二区精品小视频在线| www.色视频.com| 一个人观看的视频www高清免费观看| 亚洲第一欧美日韩一区二区三区| 在线国产一区二区在线| 午夜福利在线观看免费完整高清在 | 欧美日韩瑟瑟在线播放| 婷婷色综合大香蕉| 淫秽高清视频在线观看| 日韩高清综合在线| 欧美另类亚洲清纯唯美| 午夜日韩欧美国产| 精品久久久久久久久亚洲 | 国产欧美日韩精品亚洲av| 九九热线精品视视频播放| 嫁个100分男人电影在线观看| 久久草成人影院| 天堂av国产一区二区熟女人妻| 国产精品久久久久久久电影| avwww免费| 午夜福利免费观看在线| 久久久久国内视频| 亚洲av不卡在线观看| 特级一级黄色大片| 精品人妻1区二区| 久久中文看片网| 中亚洲国语对白在线视频| 极品教师在线视频| 欧美绝顶高潮抽搐喷水| 午夜影院日韩av| 九九热线精品视视频播放| 国产精品美女特级片免费视频播放器| 淫妇啪啪啪对白视频| .国产精品久久| 午夜福利欧美成人| 最近视频中文字幕2019在线8| aaaaa片日本免费| 久久久久免费精品人妻一区二区| 中出人妻视频一区二区| 美女被艹到高潮喷水动态| 国产午夜福利久久久久久| 宅男免费午夜| 性插视频无遮挡在线免费观看| 天堂影院成人在线观看| 国产精品久久久久久久电影| 亚洲熟妇中文字幕五十中出| 一级作爱视频免费观看| a级一级毛片免费在线观看| 国产精品久久久久久亚洲av鲁大| 有码 亚洲区| 欧美xxxx性猛交bbbb| 18+在线观看网站| 日本与韩国留学比较| 舔av片在线| 成人鲁丝片一二三区免费| 日日摸夜夜添夜夜添小说| ponron亚洲| 国产欧美日韩精品亚洲av| 别揉我奶头~嗯~啊~动态视频| 可以在线观看毛片的网站| 97碰自拍视频| 国产伦人伦偷精品视频| 亚洲人成网站在线播| 国产黄片美女视频| 国产91精品成人一区二区三区| 俺也久久电影网| 成人一区二区视频在线观看| 欧美又色又爽又黄视频| 制服丝袜大香蕉在线| 久久草成人影院| 国产毛片a区久久久久| 最好的美女福利视频网| 丰满的人妻完整版| 性色av乱码一区二区三区2| 久久久色成人| 久久人妻av系列| 嫩草影视91久久| 搡老妇女老女人老熟妇| 亚洲精品456在线播放app | 韩国av一区二区三区四区| 国产精品影院久久| 亚洲七黄色美女视频| 高清在线国产一区| 久久精品91蜜桃| 精品午夜福利在线看| 精品无人区乱码1区二区| 国内精品一区二区在线观看| 亚洲片人在线观看| 最近最新免费中文字幕在线| 十八禁国产超污无遮挡网站| 日韩人妻高清精品专区| 搡老妇女老女人老熟妇| 精品99又大又爽又粗少妇毛片 | 一个人看视频在线观看www免费| 国产欧美日韩精品亚洲av| 成年女人毛片免费观看观看9| 又紧又爽又黄一区二区| 国产av不卡久久| 一级黄片播放器| 伦理电影大哥的女人| 久久婷婷人人爽人人干人人爱| 一本一本综合久久| 少妇的逼水好多| 久久国产乱子伦精品免费另类| 久久天躁狠狠躁夜夜2o2o| 窝窝影院91人妻| 久久国产乱子伦精品免费另类| 日韩欧美精品免费久久 | 久久精品国产亚洲av香蕉五月| 亚洲av免费高清在线观看| 18+在线观看网站| 在线观看一区二区三区| 国产欧美日韩一区二区精品| 欧美+日韩+精品| 国产欧美日韩一区二区三| 精品人妻熟女av久视频| 成人高潮视频无遮挡免费网站| 国产精品久久久久久精品电影| 亚洲国产色片| h日本视频在线播放| 亚洲综合色惰| 亚洲人成电影免费在线| 欧美激情国产日韩精品一区| 国产精品,欧美在线| 热99re8久久精品国产| 内射极品少妇av片p| 美女高潮喷水抽搐中文字幕| 两人在一起打扑克的视频| 精品不卡国产一区二区三区| 丰满乱子伦码专区| 国产久久久一区二区三区| 女同久久另类99精品国产91| 欧美不卡视频在线免费观看| 亚洲激情在线av| 国内精品久久久久久久电影| 亚洲专区国产一区二区| 亚洲国产精品999在线| 日日干狠狠操夜夜爽| 国产色婷婷99| 日韩中字成人| 中国美女看黄片| 一级av片app| 日本黄大片高清| 丰满人妻一区二区三区视频av| 一个人观看的视频www高清免费观看| 夜夜躁狠狠躁天天躁| 精品久久久久久久久亚洲 | 久久6这里有精品| 级片在线观看| 免费看a级黄色片| 午夜福利视频1000在线观看| 国产精品99久久久久久久久| 亚洲18禁久久av| 99久国产av精品| 国产爱豆传媒在线观看| 夜夜爽天天搞| 丁香六月欧美| 日本 av在线| 97超级碰碰碰精品色视频在线观看| 亚洲国产高清在线一区二区三| 深爱激情五月婷婷| 天堂网av新在线| 美女免费视频网站| 一级作爱视频免费观看| 日韩精品中文字幕看吧| 欧美一区二区国产精品久久精品| 日本三级黄在线观看| 在线国产一区二区在线| 露出奶头的视频| 在线观看舔阴道视频| 99国产综合亚洲精品| 国产伦在线观看视频一区| 18禁在线播放成人免费| 国产亚洲精品av在线| 乱人视频在线观看| 久久久久免费精品人妻一区二区| 国产精品精品国产色婷婷| 欧美日韩国产亚洲二区| 俄罗斯特黄特色一大片| 搞女人的毛片| а√天堂www在线а√下载| 又黄又爽又免费观看的视频| 69av精品久久久久久| 白带黄色成豆腐渣| 老鸭窝网址在线观看| 国内精品一区二区在线观看| 午夜福利免费观看在线| 国内精品久久久久精免费| 黄色女人牲交| 婷婷丁香在线五月| 久久久久亚洲av毛片大全| 亚洲欧美清纯卡通| 欧美不卡视频在线免费观看| 国产三级黄色录像| 国产在视频线在精品| 亚洲成人精品中文字幕电影| 禁无遮挡网站| 黄色日韩在线| 男女之事视频高清在线观看| 日日摸夜夜添夜夜添小说| av国产免费在线观看| 国内精品久久久久久久电影| 熟女电影av网| 色播亚洲综合网| netflix在线观看网站| 亚洲熟妇中文字幕五十中出| 午夜久久久久精精品| 男女做爰动态图高潮gif福利片| 欧美中文日本在线观看视频| av天堂在线播放| 国产美女午夜福利| 我的老师免费观看完整版| 嫁个100分男人电影在线观看| 激情在线观看视频在线高清| 制服丝袜大香蕉在线| 中文资源天堂在线| 综合色av麻豆| 窝窝影院91人妻| 亚洲国产欧洲综合997久久,| 十八禁国产超污无遮挡网站| 一进一出好大好爽视频| 亚洲avbb在线观看| 精品久久久久久久久亚洲 | 久久久久性生活片| 91字幕亚洲| 日韩有码中文字幕| 不卡一级毛片| 午夜精品在线福利| 精品国产亚洲在线| 久久九九热精品免费| 亚洲成人免费电影在线观看| 午夜福利在线观看吧| 黄色日韩在线| 最近最新免费中文字幕在线| 少妇高潮的动态图| 日本a在线网址| 国产精品99久久久久久久久| 黄色视频,在线免费观看| 在线国产一区二区在线| 国产一区二区亚洲精品在线观看| 中文在线观看免费www的网站| 精品一区二区免费观看| 亚洲美女搞黄在线观看 | 亚洲国产精品久久男人天堂| 色综合亚洲欧美另类图片| 久久久久精品国产欧美久久久| 欧美午夜高清在线| 国产大屁股一区二区在线视频| 久久久成人免费电影| 久久精品国产亚洲av天美| 久久婷婷人人爽人人干人人爱| 国产久久久一区二区三区| 久久精品夜夜夜夜夜久久蜜豆| 亚洲av不卡在线观看| avwww免费| 亚洲久久久久久中文字幕| 婷婷六月久久综合丁香| 亚洲久久久久久中文字幕| 亚洲美女视频黄频| 亚洲第一电影网av| 免费av毛片视频| 国产 一区 欧美 日韩| 亚州av有码| 国产成人a区在线观看| 我的女老师完整版在线观看| av在线老鸭窝| 久久人人爽人人爽人人片va | 亚洲中文字幕日韩| 毛片女人毛片| 12—13女人毛片做爰片一| 简卡轻食公司| 嫩草影视91久久| 老司机午夜十八禁免费视频| 国产野战对白在线观看| 国产精品伦人一区二区| 午夜福利欧美成人| 国产成人啪精品午夜网站| 国产精品久久久久久精品电影| 久久人妻av系列| 国产黄a三级三级三级人| 人人妻人人看人人澡| 国产精品女同一区二区软件 | 免费av观看视频| 少妇裸体淫交视频免费看高清| 国产精品乱码一区二三区的特点| 91久久精品国产一区二区成人| 欧美最新免费一区二区三区 | 亚洲精品乱码久久久v下载方式| 香蕉av资源在线|