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

    不確定環(huán)境下車輛配送路徑魯棒優(yōu)化及求解算法*

    2019-12-26 06:11:04夏緒輝曹建華
    組合機床與自動化加工技術 2019年12期
    關鍵詞:魯棒遺傳算法抗體

    趙 瀟,夏緒輝,王 蕾,曹建華

    (1.武漢科技大學 機械傳動與制造工程湖北省重點實驗室,武漢 430081;2.湖北文理學院 機械工程學院,湖北 襄陽 441053)

    0 引言

    車輛配送路徑規(guī)劃是供應鏈及物流、生產制造運作管理領域的基本問題,它對于提高物料流轉效率起著至關重要的作用。因此,自車輛路徑問題提出以來便受到國內外學者的廣泛關注,時至今日仍然是該領域的熱點問題之一[1-3]。其中,難點主要體現(xiàn)在兩個方面。一是由于客觀環(huán)境的不確定性,需求以及交通運輸相關情況往往多變,由此可能導致當前方案的可行性降低,進而造成總體成本的增加或效率的降低;二是車輛路徑規(guī)劃問題一般屬于NP-Hard問題,其求解效率隨著問題規(guī)模的增加呈指數(shù)式驟然下降。配送路徑優(yōu)化是決定整個系統(tǒng)運作效率與成本的關鍵之一[4]。

    目前,針對該類問題的研究已有大量研究成果[5]。但是,在不確定環(huán)境下,多數(shù)研究成果均是在隨機規(guī)劃與動態(tài)規(guī)劃框架下完成[6]。在隨機規(guī)劃下,車輛路徑問題中的隨機因素,如客戶需求或道路通行時間,通常被看作是服從某已知概率分布的隨機事件。但是在現(xiàn)實條件下,不確定因素往往不會嚴格服務某一分布,而且,要準確獲得或求出該概率分布信息基本難以實現(xiàn)。在動態(tài)規(guī)劃框架下,按照配送階段或路徑路段可將運輸問題劃分為多個階段,基于此可構建面向車輛路徑問題的動態(tài)規(guī)劃模型。該方面雖然直觀易懂,但是,動態(tài)規(guī)劃算法由于其天生存在的維數(shù)災難問題使得其難以處理大規(guī)模的優(yōu)化問題。因此,動態(tài)規(guī)劃下的車輛路徑問題其有效性也具有一定局限性。

    針對上述兩個難題,本文將研究在不確定環(huán)境下,特別是當需求信息或運輸時間等不確定因素相關數(shù)據無法完全獲知時,如何合理規(guī)劃車輛配送路徑,以實現(xiàn)運輸成本的最小化這一為基本目標,同時控制配送或回收貨物的延遲率,進而提高配送方案在不確定環(huán)境下的運作效率。

    1 問題描述與模型構建

    1.1 問題描述

    為了模型構建的科學性和有效性,本文做出如下假設:

    (1) 每一個需求點只能被一輛配送車輛所服務;

    (2) 上游具備的庫存總量可以滿足下游所有客戶的需求;

    (3) 不考慮物資在配送中心以及運輸過程中的裝卸貨時間;

    (4) 物資運輸?shù)膯挝毁M用和配送中心的存儲費用都是已知的;

    (5) 產品的生產量都是已知的;

    (6) 從供應商到配送中心以及配送中心到需求點的距離都是已知的。

    假設當前在某一網絡包含{1,…,n}需求點,其組成集合V={1,…,n}。假設任意兩點之間的單位配送費用為cij,通行時間為tij,每個點的配送需求為di,且當前共有K輛具有相同運輸能力C的運輸車輛。在此,對于每一個車輛的配送任務而言,其可允許的總的配送時間上限假設為L,即如果假設每一輛車在初始時刻出發(fā),那么它應在L時間到達最后一個目的地。通過L可反映不同配送需求點在時間上的要求。當車輛到達時間超出L時,還將接受βk的單位懲罰成本。

    在上述背景下,要求完成所有配送需求的總成本最低。為此,涉及到的決策變量如下所示:

    uki:車輛k服務完i點后剩余的物資數(shù)量或搭載能力,

    wki:車輛到達最后需求點的時間。

    1.2 魯棒優(yōu)化模型

    現(xiàn)實條件下,由于外部環(huán)境的復雜性,獲得關鍵參數(shù)(如需求或行程時間)的精確取值或者準確的概率分布往往是極其困難甚至是不可能的。由此可能導致基于傳統(tǒng)模型得到的路徑規(guī)劃方案在現(xiàn)實中的可行性較低,甚至失效,即上述模型得到的方案在不確定環(huán)境下的魯棒性較低。為此,本部分將研究對應的魯棒優(yōu)化模型,以提高在不能得到需求或行程時間準確值時的方案可行性。魯棒優(yōu)化是當前實現(xiàn)上述目標的主要方法之一,基于當前國內外相關研究成果[7-8],本文將構建面向配送車輛路徑規(guī)劃的魯棒優(yōu)化模型。

    假設各個配送點的需求是相互獨立的,當不確定因素未知時,可假定其屬于某一不確定集合。具體地,不確定因素有各配送點需求di和兩點之間行程時間tij(不失一般性,假設任意車輛在任意相同兩點之間的行程時間相同),那么定義不確定集合U={Ud,Ut},以及

    (1)

    (2)

    Y3={y∈Rs|yTQy≤1};

    特別地,當Z滿足以下條件時,不確定集合Ut分別可構成凸包集合、箱型集合及橢圓形集合等。

    Z3={z∈Rs|zTQz≤1}。

    基于上述定義,便可得到在不確定需求或行程時間參數(shù)下的配送車輛路徑規(guī)劃魯棒優(yōu)化模型(Robust Delivery Vehicle Routing Problem,RDVRP)。

    (3)

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    ukj-uki+C(1-xkij)≥dj,?k∈K,i,j∈V,d∈Ud

    (10)

    di≤uki≤C,?i∈V,d∈Ud

    (11)

    xkij∈{0,1},?(i,j)∈V,k∈K,

    (12)

    wk,uki≥0,?k∈K

    (13)

    (14)

    (15)

    目標函數(shù)1旨在最小化總成本,包含所有車輛的運輸配送成本以及懲罰成本,其中(wk-L)+={wk-L,0}。約束條件4和5保證每個需求點都被服務一次。約束條件6保證運輸網絡內的流量平衡,即在任意需求點,車輛進出守恒。約束條件7表示每一輛車都從車輛或配送中心出發(fā)。約束條件8表示任意車輛的服務路徑不存在子回路。約束條件9表示當車輛到達時間超出L時,還將接受βk的單位懲罰成本,其中同時時間信息無法準確獲取,但屬于由約束15定義的不確定集合。約束條件10和11保證車輛剩余的貨物或搭載能力必須滿足該需求點的需要,且該需求信息無法準確獲取,但屬于由約束14定義的不確定集合。約束條件12到13是對決策變量的限制。

    基于現(xiàn)有研究成果,模型RDVRP都為0~1混合整數(shù)線性規(guī)劃問題。為保證路徑規(guī)劃方案的魯棒性,可取模型RDVRP的魯棒對應模型,即考慮在最壞情形下的不確定因素取值,此時,約束9以及10和11可改寫為:

    (16)

    (17)

    (18)

    基于此,得到的優(yōu)化模型便可計算在最大配送需求以及最長行程時間下的規(guī)劃方案,因此該方案在不確定環(huán)境下具備較高的魯棒性。

    2 面向大規(guī)模配送問題的免疫遺傳算法

    配送車輛路徑優(yōu)化問題屬于NP-Hard問題,雖然模型RDVRP可通過現(xiàn)有一些軟件直接求解,如Cplex等,但是面對大規(guī)模決策問題時,這些軟件求解效率往往較低[5]。為此,本部分將研究具有較高求解效率的啟發(fā)式算法,以提高魯棒優(yōu)化模型的可行性。

    免疫遺傳算法(Immune genetic algorithm, IGA)是免疫算法和遺傳算法的結合體[9-10]。免疫遺傳算法主要包含7個模塊:抗原識別,初始抗體產生,抗體適應度評估,記憶細胞分化,抗體促進和抑制,抗體生產和種群更新。目標函數(shù)被看著抗原,而候選解被看做抗體。可行解與最優(yōu)解間的近似程度是用抗原和抗體間的相關關系來描述的,抗體差異化的計算是保持種群多樣化的一種基本手段。

    (1)初始抗體編碼

    在基本遺傳算法和免疫遺傳算法中,抗體都需要通過編碼來表示基因型,編碼方法通常包括二進制編碼、實數(shù)編碼等。一般來說,二進制編碼比實數(shù)編碼具有更強的搜索能力,但實數(shù)編碼在變異操作中可以保持更好的種群多樣性,并且不需要解碼,可以有效地提高算法的精度和運算速度。因此,在本文配送車輛路徑優(yōu)化問題研究中,實數(shù)編碼方法將被用來對抗體進行編碼。設置初始變化區(qū)間[m(i),n(i)],φ(i)表示第i個決策變量在區(qū)間[0,1]中對應的實數(shù),g(i)表示基因。初始變化區(qū)間和實數(shù)區(qū)間滿足線性變化關系,

    φ(i)=m(i)+g(i)(n(i)-m(i))。

    通過匹配初始區(qū)間和實數(shù)區(qū)間,可以獲得一系列基因,并且這些基因按順序同問題解決方案(g(1),g(2),g(3),...,g(p))的編碼相關聯(lián)。

    (2)多樣性評估

    為了有效地克服算法過早收斂的弊端,基于近似向量距離的個體選擇概率方法將被運用到本文。抗原、細胞B和抗體分別表示目標函數(shù)、最優(yōu)解ki和候選解。第M個抗體由非空免疫系統(tǒng)集K和抗體在集合K上的距離組成,其中抗體的向量距離:

    (19)

    依據公式(19),抗體密度為:

    (20)

    因此,個體選擇概率可以表達為:

    (21)

    集合K中相似抗體數(shù)量越多,抗體i被選中的可能性越低。反之,集合K中同抗體i基因相似的抗體數(shù)越少,抗體i被選中的可能性越高。這種選擇概率使得具有有效進化基因的低適應性個體能夠獲得生殖(再生)的機會。然而,這種個體選擇的概率只同每一個抗體的適應度相關,沒有考慮個體抗體之間的編碼相似性。因此,歐幾里得距離將被引用以應對這一問題??贵wφ1,φ2,φ3,...,φn和φ1,φ2,φ3,...,φn間的歐幾里得距離可以定義為:

    (22)

    歐幾里得距離越大,2個抗體之間的相似度越低;d=0表示2個抗體是相同的。在引入歐幾里得距離后,抗體密度的選擇概率可以表示為如下:

    (23)

    (3)免疫操作

    在選擇操作過程中,以P值作為選擇概率進行抗體群體的選擇,并且輪盤賭法被用來確保好的個體可以繼承到下一代。

    變異操作是指在一個或多個位點上以一定的概率改變基因的值。變異本身是一種局部隨機搜索機制,與選擇算子相結合的它可以保證免疫遺傳算法的有效性。本文算法將隨機選擇一組個體進行變異操作以產生更好的個體:

    Mi,G=Xi,G+H(Xr1,G-Xr2,G)

    (24)

    其中,Xi,G是從優(yōu)秀個體集里隨機選擇的個體。Xr1,G和Xr2,G是從當前種群中隨機選擇的個體。Mi,G是變異操作后生成的染色體。H是調整因子,其迭代公式如下:

    (25)

    randk,k={1,2},服從[0,1]均勻分布。θ1是一個固定參數(shù),表示調整后的控制參數(shù)的概率。Hu,G和Hl,G都是固定參數(shù),分別代表控制參數(shù)H的上下界。Hi,G和Hi,G+1表示第i個個體變異過程中的調整因子。

    綜上所述,免疫遺傳算法的步驟如圖1所示。

    圖1 算法流程圖

    3 案例分析

    3.1 案例背景

    企業(yè)A是本地的大型連鎖超市,其主要業(yè)務是向本地區(qū)分銷配送某些商品。然而,隨著配送業(yè)務量的擴大,配送成本逐漸增加,延遲配送現(xiàn)象時有發(fā)生。因此,有必要研究配送車輛的配送網絡問題。在假設已知環(huán)境下,配送中心與需求點之間的需求、配送時間、最大配送時間等相關數(shù)據分別如表1、表2、表3所示。該問題中共包含40個分銷點配送需求點,其中有5個需求點,如圖2所示。目前擁有5輛型號一致的配送車輛,每輛車的配送上限為5t。

    圖2 供應鏈空間分布情況

    表1 配送中心/超市到分銷點的平均配送需求量

    表2 配送中心/超市到分銷點的平均運輸時間

    表3 配送中心/超市到分銷點的最大配送時間

    3.2 算法效果分析

    基于上述數(shù)據,本文將通過免疫遺傳算法對選擇模型進行求解。所有的數(shù)值實驗均在Windows 10.0系統(tǒng)MATLAB中完成。種群的最大數(shù)值設置為40,最大迭代次數(shù)設置為500,個體交叉的概率為0.75。

    由圖3、圖4可以看出,GA效率快于IGA。但是,IGA在解的最優(yōu)性上顯然好于GA。如圖5所示,Cplex所使用的精確算法在在迭代約400次后達到收斂狀態(tài),即相對于IGA效率較低。為進一步驗證IGA解的可靠性,如表4所示為當問題規(guī)模為50,90,130時各算法性能差異。通過對比可以發(fā)現(xiàn),IGA在求解時間上的優(yōu)勢非常明顯,解的質量雖然隨著問題規(guī)模的增加下降明顯,但相對于GA來說依然保持著較大優(yōu)勢。

    圖3 免疫遺傳算法的平均及最優(yōu)適應度曲線

    圖4 免疫遺傳算法與遺傳算法目標函數(shù)收斂對比

    圖5 Cplex收斂曲線

    此外,通過改變個體交叉概率(0.5、0.7、0.9),研究其對免疫遺傳算法的影響。圖6表明,在種群規(guī)模為20時,隨著交叉概率的增加,其對應的適應度值和最佳解的質量亦增加。最后,為了驗證不同種群大小對算法表現(xiàn)的影響,本文分別計算了種群數(shù)量為20、30、40、50、60下算法的適應度值,具體見圖7。從圖7可以發(fā)現(xiàn),當種群規(guī)模是40時,效果是最好的。

    圖6 不同交叉概率下免疫遺傳算法的適應度曲線

    圖7 不同種群大小下免疫遺傳算法的適應度曲線

    3.3 方案效果分析

    為驗證魯棒路徑優(yōu)化方案效果,本部分將通過對比一般隨機優(yōu)化模型下得到的方案加以分析。

    如圖8、圖9、表4、表5所示,傳統(tǒng)配送路徑方案和魯棒配送路徑方案在局部存在一定差別。特別是如圖中橢圓標注區(qū)域中距離配送中心較遠地帶。在魯棒優(yōu)化中,這些需求點在保證配送成本盡可能低的前提下,相對均為地分配給了配送中心,以保證配送的及時性。而傳統(tǒng)模型對配送及時性未做考慮,因此導致橢圓區(qū)域內的需求點以配送總距離最短分布,此時有些配送中心服務對象較多,有些較少。如圖10、圖11所示,通過對比配送延誤比發(fā)現(xiàn),魯棒優(yōu)化模型在各個問題規(guī)模下相對傳統(tǒng)模型均較低,且低幅隨著問題規(guī)模的增加而增大。這表明在本問題中,魯棒優(yōu)化模型可以較好的保證各個需求點配送的及時性,能更好的滿足各個點的配送需求。

    表4 算法效果對比

    表5 模型結果對比

    圖8 魯棒優(yōu)化模型所得方案

    圖9 傳統(tǒng)風險中性模型所得方案

    圖10 不同規(guī)模下配送成本對比

    圖11 不同規(guī)模下配送延誤對比

    4 結論

    本文提出了不確定環(huán)境下隨機因素數(shù)據不完全時的車輛配送路徑魯棒優(yōu)化模型以及求解算法。構建了基于魯棒優(yōu)化的車輛配送路徑規(guī)劃模型以及確定型的魯棒對等式。同時,為了提高該模型在大規(guī)模問題下的可行性,提出了基于免疫遺傳算法的求解算法。通過算例分析發(fā)現(xiàn),魯棒優(yōu)化模型可顯著提升在不確定環(huán)境下的配送及時率,相比一般隨機優(yōu)化方法,配送延誤可降低28.6%。提出的免疫遺傳算法在算法效率上顯著優(yōu)于一般遺傳算法,且對比CPLEX等精確求解方法具有明顯優(yōu)勢。

    然而,路徑規(guī)劃問題屬于典型NP-Hard問題,隨著問題規(guī)模的增加,其求解難度急劇增大,因此,研究面向大規(guī)模問題的魯棒路徑優(yōu)化問題將是下一步的主要工作。

    DOI:10.1016/j.ejor.2018.07.039.

    DOI:10.1016/j.cor.2018.02.006.

    猜你喜歡
    魯棒遺傳算法抗體
    基于學習的魯棒自適應評判控制研究進展
    自動化學報(2019年6期)2019-07-23 01:18:18
    基于自適應遺傳算法的CSAMT一維反演
    目標魯棒識別的抗旋轉HDO 局部特征描述
    自動化學報(2017年4期)2017-06-15 20:28:54
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應用
    基于遺傳算法和LS-SVM的財務危機預測
    抗BP5-KLH多克隆抗體的制備及鑒定
    基于改進的遺傳算法的模糊聚類算法
    基于Cauchy魯棒函數(shù)的UKF改進算法
    目標軌跡更新的點到點魯棒迭代學習控制
    乙肝抗體從哪兒來
    肝博士(2015年2期)2015-02-27 10:49:44
    少妇人妻精品综合一区二区 | a在线观看视频网站| 国产亚洲欧美在线一区二区| 欧美午夜高清在线| 亚洲专区中文字幕在线| 亚洲国产高清在线一区二区三| 欧美激情久久久久久爽电影| 日韩国内少妇激情av| 黄色丝袜av网址大全| 国产成人啪精品午夜网站| 成年女人毛片免费观看观看9| 色播亚洲综合网| 亚洲国产欧美人成| 久久香蕉国产精品| 亚洲精品色激情综合| 亚洲真实伦在线观看| 12—13女人毛片做爰片一| 亚洲精品一卡2卡三卡4卡5卡| 美女被艹到高潮喷水动态| 亚洲不卡免费看| 99热6这里只有精品| av女优亚洲男人天堂| 99在线视频只有这里精品首页| 亚洲精品影视一区二区三区av| 淫秽高清视频在线观看| 亚洲激情在线av| 搡老岳熟女国产| 欧美又色又爽又黄视频| 亚洲 国产 在线| 亚洲av免费在线观看| 成人av一区二区三区在线看| 国产精品一区二区三区四区免费观看 | 国产精品野战在线观看| 国产v大片淫在线免费观看| 亚洲欧美日韩高清专用| 一个人观看的视频www高清免费观看| 岛国视频午夜一区免费看| 亚洲精品在线观看二区| 精品免费久久久久久久清纯| 国产精品av视频在线免费观看| 97碰自拍视频| 两个人的视频大全免费| 国产精品99久久久久久久久| 波多野结衣巨乳人妻| 制服人妻中文乱码| 亚洲不卡免费看| 国产精品 欧美亚洲| 成人18禁在线播放| 亚洲欧美精品综合久久99| 给我免费播放毛片高清在线观看| 亚洲欧美日韩卡通动漫| 国产野战对白在线观看| 1024手机看黄色片| 国产aⅴ精品一区二区三区波| 桃红色精品国产亚洲av| av国产免费在线观看| 精品一区二区三区人妻视频| 在线国产一区二区在线| 午夜久久久久精精品| 51国产日韩欧美| 精品免费久久久久久久清纯| 亚洲精品粉嫩美女一区| 丁香欧美五月| 亚洲一区二区三区色噜噜| 久久久久国产精品人妻aⅴ院| 午夜视频国产福利| 日韩欧美在线二视频| 老司机午夜十八禁免费视频| 色播亚洲综合网| 亚洲性夜色夜夜综合| www.999成人在线观看| eeuss影院久久| 亚洲精品在线观看二区| 国内精品美女久久久久久| 免费无遮挡裸体视频| 亚洲午夜理论影院| 淫妇啪啪啪对白视频| 99久国产av精品| 久久草成人影院| 无人区码免费观看不卡| 亚洲欧美日韩高清专用| 成人无遮挡网站| 国产精品久久久久久人妻精品电影| 欧美性感艳星| 男女床上黄色一级片免费看| 3wmmmm亚洲av在线观看| 中文字幕人成人乱码亚洲影| 午夜影院日韩av| 日韩欧美精品v在线| 搡老岳熟女国产| 男女下面进入的视频免费午夜| 在线视频色国产色| 精华霜和精华液先用哪个| 精华霜和精华液先用哪个| 午夜a级毛片| 久久亚洲精品不卡| 亚洲国产精品999在线| 亚洲成人久久性| 国产亚洲精品综合一区在线观看| 精品人妻一区二区三区麻豆 | 久久精品国产自在天天线| 亚洲国产精品sss在线观看| 男女午夜视频在线观看| 亚洲自拍偷在线| 国产精品影院久久| 国产精品影院久久| 国产三级在线视频| 男人舔奶头视频| 国内少妇人妻偷人精品xxx网站| 欧美一区二区亚洲| aaaaa片日本免费| 国产亚洲av嫩草精品影院| 国产淫片久久久久久久久 | ponron亚洲| 亚洲精品一区av在线观看| 日韩 欧美 亚洲 中文字幕| eeuss影院久久| 日韩 欧美 亚洲 中文字幕| av中文乱码字幕在线| 久久精品人妻少妇| 在线观看av片永久免费下载| 亚洲中文字幕一区二区三区有码在线看| 亚洲av二区三区四区| 欧美日韩瑟瑟在线播放| 一本综合久久免费| 亚洲午夜理论影院| 在线观看舔阴道视频| 天天一区二区日本电影三级| 在线观看免费视频日本深夜| 好男人电影高清在线观看| 身体一侧抽搐| 999久久久精品免费观看国产| 亚洲中文字幕日韩| 黄色成人免费大全| 婷婷亚洲欧美| 少妇高潮的动态图| 亚洲精品一卡2卡三卡4卡5卡| 成年版毛片免费区| 欧美乱色亚洲激情| 一本综合久久免费| 内地一区二区视频在线| 麻豆成人午夜福利视频| 久久久久国内视频| 亚洲av免费高清在线观看| 一本久久中文字幕| 美女cb高潮喷水在线观看| 成人无遮挡网站| 嫁个100分男人电影在线观看| 最后的刺客免费高清国语| 一边摸一边抽搐一进一小说| 夜夜爽天天搞| 一个人观看的视频www高清免费观看| 亚洲aⅴ乱码一区二区在线播放| 两个人视频免费观看高清| av视频在线观看入口| 亚洲七黄色美女视频| 69av精品久久久久久| 久久天躁狠狠躁夜夜2o2o| 国产精华一区二区三区| 性色avwww在线观看| 桃红色精品国产亚洲av| 国产精华一区二区三区| 午夜福利视频1000在线观看| 黄片大片在线免费观看| 99精品久久久久人妻精品| 亚洲不卡免费看| 午夜福利成人在线免费观看| www.熟女人妻精品国产| 啦啦啦观看免费观看视频高清| 在线国产一区二区在线| 伊人久久精品亚洲午夜| 黄色成人免费大全| 三级国产精品欧美在线观看| 久久欧美精品欧美久久欧美| 欧美av亚洲av综合av国产av| 欧美bdsm另类| 国产精品av视频在线免费观看| 热99在线观看视频| 日韩欧美 国产精品| 免费av不卡在线播放| 色综合婷婷激情| 久久久久国产精品人妻aⅴ院| 精品福利观看| 国产精品99久久99久久久不卡| 一个人看的www免费观看视频| 三级国产精品欧美在线观看| 久久久久亚洲av毛片大全| 黄色女人牲交| 色综合站精品国产| 男女之事视频高清在线观看| 亚洲精品影视一区二区三区av| 亚洲av免费高清在线观看| 91麻豆av在线| 日日摸夜夜添夜夜添小说| 亚洲最大成人中文| 国产乱人视频| 九九久久精品国产亚洲av麻豆| 一进一出好大好爽视频| АⅤ资源中文在线天堂| 91在线精品国自产拍蜜月 | 性欧美人与动物交配| 级片在线观看| 内射极品少妇av片p| 久久香蕉精品热| 日韩欧美 国产精品| 少妇人妻一区二区三区视频| 五月玫瑰六月丁香| 欧美性猛交╳xxx乱大交人| 国产不卡一卡二| 国产91精品成人一区二区三区| 国产欧美日韩精品亚洲av| 少妇的逼水好多| 国产精品99久久99久久久不卡| 人妻夜夜爽99麻豆av| 在线视频色国产色| 亚洲av免费在线观看| 日日干狠狠操夜夜爽| 麻豆成人午夜福利视频| 精品99又大又爽又粗少妇毛片 | 欧美区成人在线视频| 亚洲片人在线观看| 一夜夜www| 禁无遮挡网站| www国产在线视频色| 精品久久久久久久人妻蜜臀av| 岛国在线观看网站| 一区二区三区激情视频| 悠悠久久av| 欧美bdsm另类| 欧美性猛交黑人性爽| 青草久久国产| 欧美中文综合在线视频| 18+在线观看网站| 黄色成人免费大全| 搡老岳熟女国产| 变态另类丝袜制服| 日本 欧美在线| 欧美丝袜亚洲另类 | av黄色大香蕉| 欧美高清成人免费视频www| 三级男女做爰猛烈吃奶摸视频| 久久欧美精品欧美久久欧美| 久久香蕉国产精品| 毛片女人毛片| 亚洲天堂国产精品一区在线| 亚洲欧美日韩东京热| 成人国产一区最新在线观看| 欧美黄色片欧美黄色片| 在线天堂最新版资源| 伊人久久大香线蕉亚洲五| 精品国产三级普通话版| 亚洲精品国产精品久久久不卡| 亚洲精华国产精华精| 午夜精品久久久久久毛片777| 91麻豆av在线| 真人做人爱边吃奶动态| 成年版毛片免费区| 嫩草影视91久久| 国内揄拍国产精品人妻在线| 国产精品日韩av在线免费观看| 99国产精品一区二区蜜桃av| 久久性视频一级片| 蜜桃亚洲精品一区二区三区| 脱女人内裤的视频| 国产精品一区二区三区四区免费观看 | 三级毛片av免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲精华国产精华精| 午夜激情福利司机影院| 亚洲天堂国产精品一区在线| 午夜日韩欧美国产| 国产精品久久久久久久电影 | 欧美日韩福利视频一区二区| av片东京热男人的天堂| 一本精品99久久精品77| 老熟妇乱子伦视频在线观看| 午夜免费成人在线视频| 久久香蕉国产精品| 午夜免费激情av| 国产成人啪精品午夜网站| 国内精品一区二区在线观看| 久久久久亚洲av毛片大全| 亚洲精品色激情综合| 日韩欧美精品免费久久 | 精品一区二区三区视频在线 | 欧美在线黄色| 国产一区二区亚洲精品在线观看| 亚洲黑人精品在线| aaaaa片日本免费| 日韩大尺度精品在线看网址| 手机成人av网站| 亚洲美女黄片视频| 久久香蕉国产精品| 99riav亚洲国产免费| 内地一区二区视频在线| 内射极品少妇av片p| 亚洲中文字幕一区二区三区有码在线看| 亚洲精品成人久久久久久| 天天添夜夜摸| 精品欧美国产一区二区三| 观看免费一级毛片| 香蕉久久夜色| 最好的美女福利视频网| 午夜免费观看网址| 亚洲精品久久国产高清桃花| 精品无人区乱码1区二区| 看片在线看免费视频| 一区二区三区激情视频| 九九久久精品国产亚洲av麻豆| 国产伦人伦偷精品视频| 一级黄片播放器| 99热精品在线国产| 久久精品91蜜桃| 日本 av在线| 亚洲av熟女| 欧美黄色淫秽网站| 看片在线看免费视频| 在线观看免费视频日本深夜| 91av网一区二区| 久久久久久大精品| 99在线视频只有这里精品首页| 丰满人妻熟妇乱又伦精品不卡| 亚洲一区二区三区不卡视频| 在线观看av片永久免费下载| 97超视频在线观看视频| 成人午夜高清在线视频| x7x7x7水蜜桃| 免费观看精品视频网站| 亚洲精品一卡2卡三卡4卡5卡| 91麻豆精品激情在线观看国产| 我要搜黄色片| 国产在视频线在精品| 欧美绝顶高潮抽搐喷水| 2021天堂中文幕一二区在线观| 久久香蕉国产精品| 久9热在线精品视频| 中文字幕人妻熟人妻熟丝袜美 | 国产精品一区二区三区四区久久| 一边摸一边抽搐一进一小说| 国产精品 欧美亚洲| 757午夜福利合集在线观看| 欧美av亚洲av综合av国产av| 淫秽高清视频在线观看| 可以在线观看的亚洲视频| 欧美一级a爱片免费观看看| 全区人妻精品视频| 午夜免费观看网址| 日韩欧美精品免费久久 | 免费观看的影片在线观看| 99久久成人亚洲精品观看| 国产精品98久久久久久宅男小说| 男女下面进入的视频免费午夜| а√天堂www在线а√下载| 女人被狂操c到高潮| 伊人久久精品亚洲午夜| 美女免费视频网站| 最近最新中文字幕大全免费视频| 精品99又大又爽又粗少妇毛片 | 欧美黑人巨大hd| 国产在视频线在精品| 国产一级毛片七仙女欲春2| 欧美性猛交黑人性爽| 熟女少妇亚洲综合色aaa.| 老司机在亚洲福利影院| 欧美黄色片欧美黄色片| 一级作爱视频免费观看| 精品久久久久久久人妻蜜臀av| 日本与韩国留学比较| 国产精品久久久久久人妻精品电影| 成人18禁在线播放| 在线天堂最新版资源| 九色成人免费人妻av| 无人区码免费观看不卡| 国产伦一二天堂av在线观看| 精华霜和精华液先用哪个| 我要搜黄色片| 老司机福利观看| 欧美最新免费一区二区三区 | 又黄又爽又免费观看的视频| av国产免费在线观看| 毛片女人毛片| 亚洲午夜理论影院| 午夜福利欧美成人| 欧美xxxx黑人xx丫x性爽| 国产成人av激情在线播放| 国产午夜精品久久久久久一区二区三区 | 国产色婷婷99| 精品电影一区二区在线| 亚洲性夜色夜夜综合| 亚洲欧美精品综合久久99| 国产毛片a区久久久久| 中文字幕人成人乱码亚洲影| 久久午夜亚洲精品久久| 窝窝影院91人妻| 国产精品久久久久久亚洲av鲁大| 欧美日韩乱码在线| 在线播放国产精品三级| 身体一侧抽搐| 噜噜噜噜噜久久久久久91| 99久久综合精品五月天人人| 成年女人看的毛片在线观看| 天堂√8在线中文| 久久久久性生活片| 免费av不卡在线播放| 国产av不卡久久| 99久久精品国产亚洲精品| 日韩精品中文字幕看吧| 精品国产美女av久久久久小说| 成年人黄色毛片网站| 国产精品野战在线观看| 村上凉子中文字幕在线| 国产黄色小视频在线观看| 日韩欧美国产在线观看| 日韩免费av在线播放| 免费看美女性在线毛片视频| 亚洲第一电影网av| 欧美日韩黄片免| 欧美性猛交╳xxx乱大交人| 天堂网av新在线| 好男人电影高清在线观看| 69av精品久久久久久| 久久久精品大字幕| 人人妻人人澡欧美一区二区| 最近最新免费中文字幕在线| 欧美高清成人免费视频www| 男女做爰动态图高潮gif福利片| 男人和女人高潮做爰伦理| 亚洲国产日韩欧美精品在线观看 | www.www免费av| 波多野结衣巨乳人妻| 国产主播在线观看一区二区| 久久久色成人| 欧美绝顶高潮抽搐喷水| 亚洲精品色激情综合| x7x7x7水蜜桃| 网址你懂的国产日韩在线| 十八禁网站免费在线| 欧美日韩亚洲国产一区二区在线观看| 中文亚洲av片在线观看爽| 九九热线精品视视频播放| 欧美极品一区二区三区四区| 全区人妻精品视频| a级毛片a级免费在线| 亚洲久久久久久中文字幕| 又粗又爽又猛毛片免费看| 亚洲av第一区精品v没综合| 中文在线观看免费www的网站| 国内精品美女久久久久久| 一夜夜www| 久久国产精品人妻蜜桃| 日日干狠狠操夜夜爽| 午夜精品久久久久久毛片777| 香蕉久久夜色| 色吧在线观看| 三级毛片av免费| 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 91麻豆av在线| 日韩欧美 国产精品| 欧美黑人巨大hd| 亚洲无线观看免费| 19禁男女啪啪无遮挡网站| 黄色丝袜av网址大全| 最好的美女福利视频网| 青草久久国产| 久久久久国产精品人妻aⅴ院| 女生性感内裤真人,穿戴方法视频| 首页视频小说图片口味搜索| 成人特级黄色片久久久久久久| 观看美女的网站| 国产又黄又爽又无遮挡在线| 午夜福利在线在线| 最后的刺客免费高清国语| 韩国av一区二区三区四区| 国产精品久久久久久久久免 | 热99re8久久精品国产| 黄色日韩在线| 淫秽高清视频在线观看| 久久久久精品国产欧美久久久| 国产高潮美女av| 免费搜索国产男女视频| 国产激情偷乱视频一区二区| 午夜a级毛片| 欧美日韩福利视频一区二区| 少妇高潮的动态图| 亚洲精品粉嫩美女一区| 亚洲人成网站在线播| 久久精品人妻少妇| a在线观看视频网站| 欧美绝顶高潮抽搐喷水| 丰满乱子伦码专区| 人妻丰满熟妇av一区二区三区| 国产精品久久久久久精品电影| 国内精品久久久久久久电影| 国产精品亚洲一级av第二区| 日本a在线网址| 9191精品国产免费久久| 久久久久久久久大av| 熟女少妇亚洲综合色aaa.| 精品人妻1区二区| tocl精华| 小说图片视频综合网站| 性色av乱码一区二区三区2| 真人做人爱边吃奶动态| 亚洲av免费在线观看| 婷婷丁香在线五月| 久久久久亚洲av毛片大全| 亚洲五月天丁香| 国产伦在线观看视频一区| АⅤ资源中文在线天堂| 国产精品永久免费网站| 村上凉子中文字幕在线| 美女 人体艺术 gogo| 1000部很黄的大片| 欧美xxxx黑人xx丫x性爽| 亚洲色图av天堂| 欧美性猛交黑人性爽| 日本一本二区三区精品| 久久精品夜夜夜夜夜久久蜜豆| 桃色一区二区三区在线观看| 国产免费一级a男人的天堂| 99riav亚洲国产免费| 一级黄色大片毛片| 欧美成人性av电影在线观看| 中文字幕av在线有码专区| www.999成人在线观看| 欧美国产日韩亚洲一区| 午夜免费成人在线视频| 亚洲人成电影免费在线| 国产精品久久久久久久电影 | 淫妇啪啪啪对白视频| 日韩 欧美 亚洲 中文字幕| 国产成人av教育| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久精品吃奶| 亚洲成av人片在线播放无| 男女床上黄色一级片免费看| 亚洲不卡免费看| 亚洲国产精品999在线| 午夜激情欧美在线| 久久久久久九九精品二区国产| 嫁个100分男人电影在线观看| 三级男女做爰猛烈吃奶摸视频| 国产伦人伦偷精品视频| 美女被艹到高潮喷水动态| 波多野结衣高清作品| 精品久久久久久久久久久久久| 日韩 欧美 亚洲 中文字幕| 亚洲第一欧美日韩一区二区三区| 日韩欧美国产一区二区入口| 日韩亚洲欧美综合| 亚洲人成网站在线播放欧美日韩| 日本免费a在线| 国产亚洲精品综合一区在线观看| 波多野结衣高清作品| 欧美国产日韩亚洲一区| 国产精品香港三级国产av潘金莲| 精品日产1卡2卡| 波多野结衣高清无吗| 一级黄片播放器| 在线观看av片永久免费下载| 青草久久国产| 亚洲国产高清在线一区二区三| 免费看十八禁软件| 波野结衣二区三区在线 | 国产中年淑女户外野战色| 99久久久亚洲精品蜜臀av| 十八禁人妻一区二区| 又紧又爽又黄一区二区| 99国产精品一区二区蜜桃av| 国产精品乱码一区二三区的特点| 国内毛片毛片毛片毛片毛片| 国产亚洲精品综合一区在线观看| 久久久久久久午夜电影| 窝窝影院91人妻| 女生性感内裤真人,穿戴方法视频| 人妻丰满熟妇av一区二区三区| 国产成年人精品一区二区| 久久国产精品人妻蜜桃| 香蕉av资源在线| 丰满人妻熟妇乱又伦精品不卡| 男人舔女人下体高潮全视频| 欧美国产日韩亚洲一区| 国产91精品成人一区二区三区| xxx96com| 亚洲av美国av| 啦啦啦韩国在线观看视频| 精品久久久久久久毛片微露脸| 亚洲中文日韩欧美视频| 综合色av麻豆| 91在线观看av| 我要搜黄色片| 免费av观看视频| or卡值多少钱| 在线播放无遮挡| 国产一区二区三区在线臀色熟女| 欧美乱码精品一区二区三区| 香蕉丝袜av| 免费看光身美女| 最后的刺客免费高清国语| 国产极品精品免费视频能看的| 最后的刺客免费高清国语| 国产视频一区二区在线看| 黄片大片在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 成人av在线播放网站| 婷婷精品国产亚洲av在线| 久久精品国产99精品国产亚洲性色| 亚洲精品成人久久久久久| 亚洲男人的天堂狠狠| 国产蜜桃级精品一区二区三区| 亚洲天堂国产精品一区在线| 日韩欧美三级三区| 成人一区二区视频在线观看| 丰满乱子伦码专区|