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

    鐵路集裝箱空箱調(diào)運(yùn)問題的遺傳算法*

    2011-08-08 04:00:04李引珍劉永莉陳志忠
    關(guān)鍵詞:空車空箱調(diào)運(yùn)

    段 剛,張 慧,陳 莉,李引珍,劉永莉,陳志忠

    (1.蘭州交通大學(xué)交通運(yùn)輸學(xué)院,甘肅 蘭州 730070;2.蘭州鐵路局貨運(yùn)處,甘肅 蘭州 730070;3.蘭州城市學(xué)院數(shù)學(xué)學(xué)院,甘肅 蘭州 730070)

    由于我國地區(qū)經(jīng)濟(jì)發(fā)展不平衡、各區(qū)域貨物流量不均衡以及集裝箱辦理站之間集裝箱到發(fā)量的不平衡,導(dǎo)致鐵路適箱貨源和箱源分布也不平衡。即使在同一地區(qū),也會(huì)因季節(jié)和時(shí)間的變化而出現(xiàn)某一段時(shí)間集裝箱空箱積壓,另一段時(shí)間短缺的現(xiàn)象,從而導(dǎo)致空箱調(diào)運(yùn)。頻繁地進(jìn)行空箱調(diào)運(yùn)不但浪費(fèi)鐵路運(yùn)能,而且空箱在運(yùn)輸途中支出的各種費(fèi)用也給鐵路運(yùn)輸企業(yè)增加了負(fù)擔(dān)。集裝箱的日常管理費(fèi)用通常可以分為:空箱調(diào)運(yùn)費(fèi)、堆存費(fèi)、修理費(fèi)、集裝箱租賃費(fèi)、吊裝費(fèi)和其他費(fèi)用,而集裝箱空箱調(diào)運(yùn)費(fèi)用占總費(fèi)用的1/4。因此,優(yōu)化空箱調(diào)運(yùn)方案并設(shè)計(jì)出高效的算法,對于減少空箱調(diào)運(yùn)成本,提高鐵路運(yùn)能和集裝箱使用效率,具有十分現(xiàn)實(shí)的意義。

    1 研究現(xiàn)狀

    鐵路空箱調(diào)運(yùn)與空車調(diào)配問題無論從模型構(gòu)造上還是求解方法上都非常相似,因此可以借鑒。熊紅云等[1]和張得志等[2]采用矩陣表示染色體的遺傳算法分別求解鐵路空車和空箱調(diào)運(yùn)問題,他們設(shè)計(jì)的交叉算子均對2個(gè)染色體的均值運(yùn)算和取余運(yùn)算的結(jié)果進(jìn)行代數(shù)和計(jì)算,從而得到新的染色體,變異算子也都是對染色體子矩陣中的元素進(jìn)行調(diào)整,區(qū)別是一個(gè)采用整數(shù)編碼,一個(gè)采用十進(jìn)制編碼。果鵬文等[3]采用振蕩法解決大規(guī)模路網(wǎng)上的空車調(diào)配問題,在事先不知道某個(gè)區(qū)段空車排空方向的前提下,預(yù)先人為指定該區(qū)段的空車排空方向,作為初始方案,使得區(qū)段上中間站的空車流,能夠按照一定的原則歸并到前方技術(shù)站;然后對初始方案進(jìn)行計(jì)算,對區(qū)段空車方向不斷進(jìn)行調(diào)整,反復(fù)計(jì)算,直到指定的空車方向與計(jì)算結(jié)果相符合時(shí)為止。雷中林等[4]針對運(yùn)輸需求等不確定性建立了空車調(diào)配的隨機(jī)機(jī)會(huì)約束模型,并設(shè)計(jì)了遺傳算法求解。林柏梁等[5]設(shè)計(jì)了分步優(yōu)化迭代算法求解線路能力約束下的鐵路空車調(diào)配問題。還有應(yīng)用神經(jīng)網(wǎng)絡(luò)法[6]、物體重心法[7]和區(qū)段中心優(yōu)化法[8]求解問題。在此,作者從鐵路運(yùn)輸實(shí)際出發(fā),針對空箱調(diào)運(yùn)模型,設(shè)計(jì)了一種新的遺傳算法,可以高效求解空箱調(diào)運(yùn)問題。

    2 空箱調(diào)運(yùn)問題的數(shù)學(xué)模型

    2.1 參數(shù)與變量

    S為空箱供應(yīng)站集合,S={1,2,…,m};D為空箱需求站集合,D={1,2,…,n};ai為 i站可供應(yīng)的空箱數(shù),ai∈Z+,i∈S;bj為 j站需求的空箱數(shù),bj∈Z+,j∈D;dij為 i站與 j站間的距離,i∈S,j∈D;xij為i站向j站調(diào)運(yùn)的空箱數(shù),決策變量。

    2.2 空箱調(diào)運(yùn)模型

    空箱調(diào)運(yùn)模型如下:

    (1)式為目標(biāo)函數(shù),表示空箱調(diào)運(yùn)的總距離最小;(2)式與(3)式為空箱供應(yīng)與需求的平衡約束;(4)式為變量的非負(fù)和整數(shù)約束。

    3 空箱調(diào)運(yùn)問題的遺傳算法

    自從1975年Holland提出遺傳算法后,因?yàn)槠浜唵?、魯棒性好和并行?jì)算等特點(diǎn),從而在工程技術(shù)、自動(dòng)控制、圖像處理、人工智能和管理科學(xué)等眾多領(lǐng)域得到了廣泛應(yīng)用。本文提出一種新的遺傳算法,通過交叉與變異操作,在保證解的可行性基礎(chǔ)上,能夠快速高效求解空箱調(diào)運(yùn)問題。

    3.1 染色體的構(gòu)成及種群的初始化

    采用整數(shù)編碼,以矩陣V表示一個(gè)染色體,

    其中,xij∈Z+∪{0},i∈S,j∈D 。顯然,基因 xij表示i站向j站調(diào)運(yùn)的空箱量,這樣,染色體的基因型與其表現(xiàn)型完全一致,無需解碼,因此以下對染色體和解不加區(qū)分。但要使得染色體可行,必須滿足供需平衡的條件。以下種群的初始化,交叉算子和變異算子的操作,均圍繞染色體的可行性來進(jìn)行。

    初始種群的產(chǎn)生采用下面的算法[1-2]:

    步驟1 隨機(jī)產(chǎn)生下標(biāo) i和 j,其中i∈ S,j∈D;

    步驟2 計(jì)算 xij,xij=min{ai,bj};

    步驟3 更新 ai和 bj,ai= ai- xij,bj= bjxij;

    步驟4 重復(fù)步驟1—步驟3,直到所有ai和bj都為0。

    由上述算法產(chǎn)生的染色體必然對應(yīng)一組可行解,其分配思想就是隨機(jī)化的西北角法。

    3.2 交叉算子

    按交叉概率pc隨機(jī)選擇染色體V1和V2作為父個(gè)體,其中。為了使得交叉操作不破壞染色體的可行性,采用下面的方法產(chǎn)生子代染色體V3=(yij)和V4=(zij)。首先,令

    其中,i∈ S,j∈ D,i和 j不同時(shí)為1,0 < α ≤0.5。當(dāng)i=1,j≥2時(shí),(9)式和(11)式的最小值皆取第1項(xiàng),而當(dāng)j=1,i≥2時(shí),兩式的最小值皆取第2項(xiàng)。α為交叉參數(shù),是本遺傳算法中最重要的參數(shù)。實(shí)際上,文獻(xiàn)[1]和[2]中的交叉方法可以看作是本交叉算子的特例,即α取值為0.5,也就是對2個(gè)父代染色體的對應(yīng)基因求均值。而本文中,α取值更廣泛,這樣就能夠擴(kuò)大種群的多樣性。

    (9)式和(11)式可以保證子代染色體每行的調(diào)運(yùn)量之和不會(huì)超過供應(yīng)量,每列的調(diào)運(yùn)量之和不會(huì)超過需求量,(10)式和(12)式能夠保證調(diào)運(yùn)量都非負(fù),但卻不能保證等式成立,故需要對調(diào)運(yùn)量進(jìn)行調(diào)整以滿足平衡條件。

    更新 yi1,j1,令 yi1,j1= φ +yi1,j1;遍歷 i,j,直至所有行列都滿足平衡。

    3.3 變異算子

    求解運(yùn)輸問題的表上作業(yè)法,其對基變量的調(diào)整采用閉合回路法。該方法雖然能保證解的可行性,但尋找閉合回路卻相當(dāng)耗時(shí)。本文借鑒該思想,采用基于矩形閉合回路的方法設(shè)計(jì)變異算子,可以大大減少計(jì)算量。

    按變異概率pm選擇染色體V5=(wij),隨機(jī)選擇 i0,j0,i1,j1,其中 i0∈S,j0∈D,i1∈S/{i0},j1∈ D/{j0}。當(dāng) wi0,j0wi1,j1> 0 ,即 wi0,j0和 wi1,j1同時(shí)為正數(shù),取δ=min{wi0j0,wi1j1}。然后在矩形閉合回路 wi0j0,wi0j1,wi1j1,wi1j0中進(jìn)行調(diào)整,令

    其余基因wij取值不變,組成新的染色體V6。

    3.4 評價(jià)函數(shù)和選擇操作

    由于染色體的可行性以及目標(biāo)函數(shù)為線性函數(shù),所以取目標(biāo)函數(shù)為評價(jià)函數(shù),即:

    其中,V=(xij)。

    采用基于輪盤賭的選擇策略,詳見文獻(xiàn)[9]。

    3.5 算法步驟

    步驟1 初始化種群;

    步驟2 通過交叉操作和變異操作產(chǎn)生子代染色體;

    步驟3 計(jì)算每個(gè)染色體的評價(jià)函數(shù),并記錄最優(yōu)染色體;

    步驟4 通過輪盤賭選擇染色體;

    步驟5 重復(fù)步驟2—步驟4,直到滿足停止條件。

    其中,停止條件為可以定義為進(jìn)化的最大代數(shù),也可以通過比較前后兩代最優(yōu)染色體的變化情況確定停止條件,本文采用前者。

    4 實(shí)驗(yàn)仿真

    以蘭州鐵路局為例,對該局內(nèi)集裝箱辦理站間的空箱調(diào)運(yùn)進(jìn)行優(yōu)化。圖1所示為蘭州局的18個(gè)集裝箱辦理站。這里取2009年某月的集裝箱運(yùn)輸數(shù)據(jù),用以驗(yàn)證算法的有效性。各個(gè)集裝箱辦理站的空箱需求量、供應(yīng)量以及兩站之間的距離如表1所示,其中,距離為兩站之間的最短路徑長度。由于綠化、中堡、張家祠和鴛鴦鎮(zhèn)4個(gè)集裝箱辦理站到發(fā)集裝箱數(shù)量平衡,故沒有出現(xiàn)在表1中。因?yàn)榭傂枨罅繛?54箱,而總供應(yīng)量只有344箱,所以增加虛設(shè)供應(yīng)站S6,其供應(yīng)量為210箱。

    圖1 蘭州鐵路局集裝箱辦理站示意圖Fig.1 Lanzhou Railway Bureau container handling stations

    為了尋求交叉概率、變異概率以及交叉參數(shù)α的最佳取值,利用 VC6.0編程,在 Core 2 Duo&2.2GHz、1G DDR2的惠普筆記本上運(yùn)行通過,在每組概率組合下都做了10次測試,其中群體規(guī)模為20,遺傳代數(shù)為600。首先令 α =0.1,Pm=0.1,交叉概率Pc在(0,1)區(qū)間每次變化0.05,評價(jià)函數(shù)的平均值和標(biāo)準(zhǔn)差的變化見圖2和圖3。顯然 Pc的最優(yōu)取值為0.65,在此基礎(chǔ)上,變異概率Pm在(0,1)內(nèi)每次變化0.05,同樣可以得到評價(jià)函數(shù)的均值和標(biāo)準(zhǔn)差變化情況,見圖4和圖5。從圖中可以看出,當(dāng)Pm大于0.65時(shí),評價(jià)函數(shù)的均值幾乎相同,其標(biāo)準(zhǔn)差也比較接近,且都比較小,因此,選取 Pm=0.65。

    表1 集裝箱辦理站間距離及供需表Table 1 The distance between container handling stations and supply and demand quantity

    圖2 不同交叉概率下評價(jià)函數(shù)平均值Fig.2 The average value of evaluation function under different crossover probability

    圖3 不同交叉概率下評價(jià)函數(shù)標(biāo)準(zhǔn)差Fig.3 The standard deviation value of evaluation function under different crossover probability

    圖4 不同變異概率下評價(jià)函數(shù)平均值Fig.4 The average value of evaluation function under different mutation probability

    圖5 不同變異概率下評價(jià)函數(shù)標(biāo)準(zhǔn)差Fig.5 The standard deviation value of evaluation function under different mutation probability

    在最優(yōu)交叉概率和變異概率組合下,令α在區(qū)間[0.05,0.5]內(nèi)每次變化0.05,得到評價(jià)函數(shù)的均值和標(biāo)準(zhǔn)差,如圖6和圖7所示。當(dāng)α取值為0.1,0.35和0.5時(shí),評價(jià)函數(shù)平均值都比較小,標(biāo)準(zhǔn)差也相差不多,所以可以任選一個(gè),本文取α=0.1。

    圖6 不同交叉參數(shù)值下評價(jià)函數(shù)均值Fig.6 The average value of evaluation function under different crossover parameter

    圖7 不同交叉參數(shù)值下評價(jià)函數(shù)標(biāo)準(zhǔn)差Fig.7 The standard deviation value of evaluation function under different crossover probability

    在以上各參數(shù)最優(yōu)取值下,通過本文設(shè)計(jì)的遺傳算法,當(dāng)計(jì)算到520代時(shí)達(dá)到最優(yōu),評價(jià)函數(shù)為143689,收斂過程如圖8所示(其中,從390代到510代,評價(jià)函數(shù)都為143720,與最優(yōu)評價(jià)函數(shù)只相差31,圖中顯示不明顯)。采用LINGO 9.0軟件計(jì)算得到的最優(yōu)解如表2所示,計(jì)算時(shí)間不到1 s(LINGO 9.0最小計(jì)時(shí)單位為s)。本遺傳算法的計(jì)算時(shí)間僅為0.14 s,不但能得到上述最優(yōu)解,而且還能得到另外一組最優(yōu)解,如表3所示。不難看出,在表2中由金昌,蘭州北,嘉峪關(guān)和白銀市組成的閉合回路里,對調(diào)運(yùn)量進(jìn)行了1個(gè)單位的調(diào)整,即得到表3,而這正是本文中的變異操作。所以使用該遺傳算法,不僅能在不同的計(jì)算中得到這2組最優(yōu)解,還能在一次計(jì)算中同時(shí)得到這2組最優(yōu)解,而這是LINGO軟件做不到的。

    圖8 遺傳算法收斂圖Fig.8 Genetic algorithm convergence figure

    表2 應(yīng)用LINGO計(jì)算得到的最優(yōu)調(diào)運(yùn)量Table 2 The best solution solved by LINGO software container

    表3 另外一組最優(yōu)調(diào)運(yùn)量Table 3 The other best solution solved by genetic algorithm container

    5 結(jié)論

    鐵路集裝箱空箱調(diào)運(yùn)優(yōu)化模型與空車調(diào)運(yùn)模型完全相同,因此,設(shè)計(jì)高效的算法就成為解決這類問題的關(guān)鍵。本文提出的遺傳算法,交叉操作是對父代染色體進(jìn)行線性組合并取整后經(jīng)過適當(dāng)調(diào)整得到子代染色體,變異操作則采用矩形閉合回路法。通過對蘭州鐵路局集裝箱辦理站的實(shí)例計(jì)算,表明該算法不但效率非常高,而且還能得到不同的最優(yōu)解。關(guān)于交叉參數(shù)α的取值問題,本文采用固定值,通過測試找到其最佳取值。如果在計(jì)算過程中,對α進(jìn)行調(diào)整,采用變參數(shù)法,或者直接對其進(jìn)行編碼并參與交叉和變異操作進(jìn)行自適應(yīng)調(diào)整,這樣會(huì)進(jìn)一步提高該算法的效率,以利于求解更大規(guī)模的問題,這將是我們下一步要研究的內(nèi)容。

    [1]熊紅云,魯五一,溫紅艷.鐵路空車調(diào)配問題的遺傳啟發(fā)算法[J].中國鐵道科學(xué),2002,23(4):118-121.XIONG Hong-yun,LU Wu-yi,WEN Hong-yan.Hereditary heuristic algorithm for empty car distribution in railway transportation[J].China Railway Science,2002,23(4):118-121.

    [2]張得志,謝如鶴,黃孝章.鐵路集裝箱空箱調(diào)度模型及求解算法[J].中國鐵道科學(xué),2003,24(3):125-129.ZHANG De-zhi,XIE Ru-he,HUANG Xiao-zhang.Dispatch model and solution algorithm of railway empty container[J].China Railway Science,2003,24(3):125 -129.

    [3]果鵬文,褚 江,林柏梁.用振蕩法解大規(guī)模路網(wǎng)上的空車調(diào)配問題[J].中國鐵道科學(xué),2002,23(4):111-117.GUO Peng-wen,CHU Jiang,LIN Bo-liang.Reiterative adjusting method of empty wagon on large-scale railway network[J].China Railway Science,2002,23(4):111-117.

    [4]雷中林,何世偉,宋 瑞,等.鐵路空車調(diào)配問題的隨機(jī)機(jī)會(huì)約束模型及遺傳算法[J].鐵道學(xué)報(bào),2005,27(5):1-5.LEI Zhong-lin,HE Shi-wei,SONG Rui,et al.Stochastic chance-constrained model and genetic algorithm for empty car distribution in railway transportation[J].Journal ofthe China Railway Society,2005,27(5):1-5.

    [5]林柏梁,喬國會(huì).基于線路能力約束下的鐵路空車調(diào)配迭代算法[J].中國鐵道科學(xué),2008,29(1):93-96.LIN Bo-liang,QIAO Guo-hui.Iterative algorithm of railway network empty cars distribution based on restriction of route capacity[J].China Railway Science,2008,29(1):93-96.

    [6]黨建武.神經(jīng)網(wǎng)絡(luò)在鐵路空車調(diào)度問題中的應(yīng)用[J].蘭州鐵道學(xué)院學(xué)報(bào),1999,18(1):77-85.DANG Jian-wu.Application of neural networks to empty carriage scheduling in railway transportation[J].Journal of Lanzhou Railway Institute,1999,18(1):77-85.

    [7]紀(jì)嘉倫,林柏梁,李福志,等.用重心優(yōu)化方法求解鐵路網(wǎng)上空車調(diào)配問題[J].鐵道學(xué)報(bào),2001,23(3):109-113.JI Jia-lun,LIN Bo-linag,LI Fu-zhi,et al.Distribution of empty wagons on railway network by using the center-of-gravity- optimization method[J].Journal of the China Railway Society,2001,23(3):109-113.

    [8]果鵬文,林柏梁,余 洋.大規(guī)模路網(wǎng)上空車調(diào)配的區(qū)段中心優(yōu)化法[J].中國鐵道科學(xué),2001,22(2):122-128.GUO Peng-wen,LIN Bo-liang,YU Yang.The section center optimization method of empty car adjustment on large - scale railway network[J].China Railway Science,2001,22(2):122-128.

    [9]Liu B.Stackelberg-Nash equilibrium for multilevel programming with multiple followers using genetic algorithms[J].Computers and Mathematics with Applications,1998,36(7):79-89.

    猜你喜歡
    空車空箱調(diào)運(yùn)
    基于人員分配的艦載機(jī)出動(dòng)調(diào)運(yùn)指揮模型
    以翻車機(jī)空車線為例對自動(dòng)防溜系統(tǒng)的分析和思考
    農(nóng)業(yè)部:鼓勵(lì)規(guī)模養(yǎng)殖,集中屠宰,限制畜禽調(diào)運(yùn)
    火車翻車機(jī)空車調(diào)車系統(tǒng)的優(yōu)化改進(jìn)
    山東冶金(2017年2期)2017-05-10 08:20:50
    基于ANSYS空箱扶壁式高大翼墻動(dòng)力分析
    集裝箱碼頭殘損空箱規(guī)范化管理措施
    集裝箱化(2016年12期)2017-03-20 02:30:53
    新型集裝箱設(shè)計(jì)將減少空箱運(yùn)輸量
    基于時(shí)間窗的鐵路重載運(yùn)輸空車回送優(yōu)化
    基于時(shí)間約束的鐵路空車調(diào)配系統(tǒng)可靠性分析
    調(diào)運(yùn)肉牛應(yīng)激反應(yīng)繼發(fā)癥的診斷和治療
    亚洲七黄色美女视频| 好看av亚洲va欧美ⅴa在| 亚洲国产欧美一区二区综合| 欧美乱码精品一区二区三区| 两性夫妻黄色片| 女人高潮潮喷娇喘18禁视频| 人妻久久中文字幕网| 国产伦人伦偷精品视频| www.自偷自拍.com| 免费无遮挡裸体视频| 国产又黄又爽又无遮挡在线| 国产真实乱freesex| 国内毛片毛片毛片毛片毛片| 久久99热这里只有精品18| 婷婷精品国产亚洲av在线| 色老头精品视频在线观看| 免费在线观看成人毛片| 成人午夜高清在线视频 | 天天躁夜夜躁狠狠躁躁| av电影中文网址| 亚洲国产精品久久男人天堂| 最近最新中文字幕大全免费视频| 一二三四在线观看免费中文在| 99国产综合亚洲精品| 熟女少妇亚洲综合色aaa.| 午夜福利一区二区在线看| 日韩免费av在线播放| 一二三四在线观看免费中文在| 两性午夜刺激爽爽歪歪视频在线观看 | 不卡一级毛片| 中文字幕高清在线视频| 国内久久婷婷六月综合欲色啪| 欧美乱妇无乱码| 成人免费观看视频高清| 色综合亚洲欧美另类图片| 欧美av亚洲av综合av国产av| 动漫黄色视频在线观看| 亚洲中文日韩欧美视频| 50天的宝宝边吃奶边哭怎么回事| 欧美另类亚洲清纯唯美| 中出人妻视频一区二区| 99在线视频只有这里精品首页| 久久久精品欧美日韩精品| 国产在线观看jvid| 久久精品人妻少妇| 精品不卡国产一区二区三区| 超碰成人久久| xxx96com| 日本在线视频免费播放| aaaaa片日本免费| 国产成人影院久久av| 91av网站免费观看| 久久香蕉国产精品| 99精品欧美一区二区三区四区| 波多野结衣av一区二区av| 黄色 视频免费看| 国产精品亚洲一级av第二区| 亚洲精品在线美女| 精品一区二区三区av网在线观看| x7x7x7水蜜桃| 搡老妇女老女人老熟妇| 啦啦啦 在线观看视频| 在线天堂中文资源库| 精品国产一区二区三区四区第35| 亚洲真实伦在线观看| 久久久久久久久中文| 黑人操中国人逼视频| 一级毛片高清免费大全| 91大片在线观看| 欧美一级a爱片免费观看看 | 这个男人来自地球电影免费观看| 在线永久观看黄色视频| 久久精品国产亚洲av香蕉五月| 亚洲精品在线观看二区| 99国产极品粉嫩在线观看| 亚洲成人久久性| 亚洲中文字幕日韩| 免费在线观看视频国产中文字幕亚洲| 美女扒开内裤让男人捅视频| 成人手机av| 亚洲狠狠婷婷综合久久图片| 国产av不卡久久| 波多野结衣高清作品| 岛国视频午夜一区免费看| 一级a爱片免费观看的视频| 午夜免费成人在线视频| 国产精品久久久久久亚洲av鲁大| 嫩草影视91久久| 亚洲成人久久性| 日韩欧美 国产精品| 国产一区二区三区视频了| 久久久久久九九精品二区国产 | 国产激情偷乱视频一区二区| 中文字幕av电影在线播放| 人人澡人人妻人| 中文字幕最新亚洲高清| 人妻久久中文字幕网| 亚洲avbb在线观看| 精品午夜福利视频在线观看一区| 桃色一区二区三区在线观看| 亚洲欧美激情综合另类| 久久久久久亚洲精品国产蜜桃av| 亚洲成a人片在线一区二区| 91麻豆av在线| 国内久久婷婷六月综合欲色啪| 一本精品99久久精品77| 熟女少妇亚洲综合色aaa.| 亚洲avbb在线观看| 欧美另类亚洲清纯唯美| 观看免费一级毛片| 丰满人妻熟妇乱又伦精品不卡| 国产aⅴ精品一区二区三区波| 人妻丰满熟妇av一区二区三区| 精品国产乱子伦一区二区三区| 亚洲在线自拍视频| 成人一区二区视频在线观看| 一级毛片女人18水好多| 国产在线精品亚洲第一网站| 亚洲天堂国产精品一区在线| 亚洲国产精品久久男人天堂| 美女扒开内裤让男人捅视频| 首页视频小说图片口味搜索| 两人在一起打扑克的视频| 亚洲成a人片在线一区二区| 久久草成人影院| 老司机午夜福利在线观看视频| 国产极品粉嫩免费观看在线| 午夜福利在线观看吧| 亚洲狠狠婷婷综合久久图片| 国产免费av片在线观看野外av| 久久久国产欧美日韩av| 亚洲精品一区av在线观看| 亚洲精品一区av在线观看| 国产一区二区三区视频了| 在线观看66精品国产| 欧美日韩精品网址| 欧美精品亚洲一区二区| 91九色精品人成在线观看| 高潮久久久久久久久久久不卡| 中文字幕久久专区| 亚洲在线自拍视频| 黄色丝袜av网址大全| 国产人伦9x9x在线观看| 国内精品久久久久精免费| 久久 成人 亚洲| 又黄又爽又免费观看的视频| 中文亚洲av片在线观看爽| 亚洲专区中文字幕在线| 美女午夜性视频免费| 日韩欧美三级三区| 国产主播在线观看一区二区| 美女高潮喷水抽搐中文字幕| 国产一区二区三区在线臀色熟女| 91av网站免费观看| а√天堂www在线а√下载| 在线天堂中文资源库| 国产亚洲欧美在线一区二区| 亚洲精品中文字幕一二三四区| 母亲3免费完整高清在线观看| 男人的好看免费观看在线视频 | 亚洲熟妇熟女久久| 国产三级黄色录像| 久久中文字幕人妻熟女| 午夜精品在线福利| 国产高清有码在线观看视频 | 又黄又爽又免费观看的视频| 亚洲五月色婷婷综合| 男人舔奶头视频| 18禁黄网站禁片免费观看直播| 丰满的人妻完整版| 午夜福利免费观看在线| 国产一区二区三区在线臀色熟女| 日日干狠狠操夜夜爽| 欧美色视频一区免费| xxxwww97欧美| 午夜久久久久精精品| 91九色精品人成在线观看| 亚洲熟女毛片儿| 成年版毛片免费区| 国产乱人伦免费视频| 国产亚洲av高清不卡| 亚洲 国产 在线| 欧美av亚洲av综合av国产av| 两个人免费观看高清视频| 变态另类丝袜制服| 亚洲七黄色美女视频| 欧美成狂野欧美在线观看| 夜夜夜夜夜久久久久| 欧美大码av| 脱女人内裤的视频| 国产免费男女视频| 两性午夜刺激爽爽歪歪视频在线观看 | 久久精品国产亚洲av高清一级| 国产精品自产拍在线观看55亚洲| 黄频高清免费视频| 久久精品成人免费网站| 国产熟女午夜一区二区三区| 免费在线观看成人毛片| 久久午夜综合久久蜜桃| 日日夜夜操网爽| 女性被躁到高潮视频| 成人特级黄色片久久久久久久| 女人高潮潮喷娇喘18禁视频| 亚洲精品美女久久av网站| 亚洲精品中文字幕一二三四区| 在线观看日韩欧美| 精品久久蜜臀av无| 精品国内亚洲2022精品成人| 两性午夜刺激爽爽歪歪视频在线观看 | 成熟少妇高潮喷水视频| www.熟女人妻精品国产| 中文字幕精品免费在线观看视频| 亚洲中文字幕日韩| 美女免费视频网站| 国产欧美日韩一区二区三| 视频区欧美日本亚洲| 久久天躁狠狠躁夜夜2o2o| 免费在线观看黄色视频的| 桃色一区二区三区在线观看| 欧美国产日韩亚洲一区| 中文字幕高清在线视频| 黄网站色视频无遮挡免费观看| 一级a爱片免费观看的视频| 99精品欧美一区二区三区四区| 亚洲av第一区精品v没综合| 99久久国产精品久久久| svipshipincom国产片| 中文字幕精品免费在线观看视频| 国产高清视频在线播放一区| 日韩成人在线观看一区二区三区| 99热只有精品国产| 亚洲精华国产精华精| 婷婷丁香在线五月| 一级毛片高清免费大全| 欧美三级亚洲精品| 久久精品夜夜夜夜夜久久蜜豆 | 12—13女人毛片做爰片一| 欧美黄色片欧美黄色片| 欧美性猛交╳xxx乱大交人| 18禁黄网站禁片午夜丰满| 一本精品99久久精品77| 午夜激情av网站| 国产野战对白在线观看| 长腿黑丝高跟| 99riav亚洲国产免费| 亚洲av电影不卡..在线观看| 日韩国内少妇激情av| 午夜影院日韩av| 极品教师在线免费播放| 午夜老司机福利片| 黄色毛片三级朝国网站| 欧美日韩黄片免| 亚洲成人国产一区在线观看| 午夜福利在线观看吧| 国产精品永久免费网站| 久久久久国产精品人妻aⅴ院| 日韩成人在线观看一区二区三区| 久久久久久国产a免费观看| 欧美色欧美亚洲另类二区| 免费一级毛片在线播放高清视频| 欧美性长视频在线观看| 日本五十路高清| 一二三四社区在线视频社区8| 青草久久国产| 高清在线国产一区| 国产亚洲精品综合一区在线观看 | 欧美精品啪啪一区二区三区| 91字幕亚洲| 变态另类丝袜制服| 欧美人与性动交α欧美精品济南到| 老司机福利观看| 日本熟妇午夜| √禁漫天堂资源中文www| 国产亚洲精品综合一区在线观看 | 丝袜在线中文字幕| 亚洲激情在线av| 又大又爽又粗| 丰满人妻熟妇乱又伦精品不卡| 国产精品99久久99久久久不卡| 黄色女人牲交| 欧美黑人巨大hd| 色尼玛亚洲综合影院| 精品不卡国产一区二区三区| 一区二区日韩欧美中文字幕| 满18在线观看网站| www国产在线视频色| 亚洲va日本ⅴa欧美va伊人久久| 国产又爽黄色视频| 性色av乱码一区二区三区2| 中文字幕精品亚洲无线码一区 | 女性被躁到高潮视频| 亚洲自拍偷在线| 性色av乱码一区二区三区2| 国产日本99.免费观看| 精品国产国语对白av| 午夜福利高清视频| 哪里可以看免费的av片| 精品熟女少妇八av免费久了| 国产aⅴ精品一区二区三区波| 亚洲av五月六月丁香网| 久热这里只有精品99| 最近最新免费中文字幕在线| 国产精品久久电影中文字幕| 国产精品免费视频内射| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区在线臀色熟女| 真人一进一出gif抽搐免费| 欧美大码av| 亚洲专区国产一区二区| 在线永久观看黄色视频| 男女视频在线观看网站免费 | 国产欧美日韩一区二区精品| 日本在线视频免费播放| 日本成人三级电影网站| av欧美777| 日韩av在线大香蕉| 日韩欧美在线二视频| 老司机深夜福利视频在线观看| 国产精品综合久久久久久久免费| 黄片小视频在线播放| av视频在线观看入口| 成年女人毛片免费观看观看9| 国产成+人综合+亚洲专区| 少妇粗大呻吟视频| 精品乱码久久久久久99久播| 国产av又大| 啦啦啦韩国在线观看视频| 999精品在线视频| 女生性感内裤真人,穿戴方法视频| 久9热在线精品视频| 成人欧美大片| 日韩欧美国产一区二区入口| 成人亚洲精品一区在线观看| aaaaa片日本免费| 免费在线观看视频国产中文字幕亚洲| 久久久国产欧美日韩av| 免费高清在线观看日韩| www国产在线视频色| 国产免费男女视频| 国产在线观看jvid| 搡老妇女老女人老熟妇| 午夜久久久在线观看| 亚洲自偷自拍图片 自拍| 美女大奶头视频| 老鸭窝网址在线观看| 中文在线观看免费www的网站 | 欧美精品啪啪一区二区三区| 久久精品人妻少妇| 国产伦人伦偷精品视频| 可以免费在线观看a视频的电影网站| 欧美日本视频| 欧美又色又爽又黄视频| 日本免费a在线| 老司机福利观看| 妹子高潮喷水视频| 亚洲,欧美精品.| 精品熟女少妇八av免费久了| netflix在线观看网站| 97碰自拍视频| 黄色视频,在线免费观看| 国产区一区二久久| 岛国在线观看网站| 久久精品国产清高在天天线| 后天国语完整版免费观看| 久久青草综合色| 午夜视频精品福利| 久久人妻福利社区极品人妻图片| 国产亚洲欧美精品永久| 一级黄色大片毛片| 日韩三级视频一区二区三区| 日本撒尿小便嘘嘘汇集6| 手机成人av网站| 日本 av在线| 香蕉国产在线看| 人妻丰满熟妇av一区二区三区| 久久欧美精品欧美久久欧美| 国产黄色小视频在线观看| 中文在线观看免费www的网站 | 中文亚洲av片在线观看爽| 国产av不卡久久| av片东京热男人的天堂| 欧美激情久久久久久爽电影| 午夜免费观看网址| 日韩有码中文字幕| 两性午夜刺激爽爽歪歪视频在线观看 | 国产成人啪精品午夜网站| 久久香蕉精品热| 亚洲人成77777在线视频| 欧美另类亚洲清纯唯美| 日韩中文字幕欧美一区二区| 国产精品爽爽va在线观看网站 | 91在线观看av| 国产私拍福利视频在线观看| 两性夫妻黄色片| 黑人操中国人逼视频| 精品国产一区二区三区四区第35| 香蕉国产在线看| 人妻丰满熟妇av一区二区三区| 亚洲精品在线观看二区| svipshipincom国产片| 色综合站精品国产| 怎么达到女性高潮| 叶爱在线成人免费视频播放| 欧美日韩福利视频一区二区| 一区二区三区国产精品乱码| 在线视频色国产色| 桃红色精品国产亚洲av| 91九色精品人成在线观看| 精品国产一区二区三区四区第35| 制服诱惑二区| 两性午夜刺激爽爽歪歪视频在线观看 | 精品久久久久久久人妻蜜臀av| 精品国产一区二区三区四区第35| 香蕉国产在线看| 欧美精品啪啪一区二区三区| 久久精品影院6| 亚洲成国产人片在线观看| 欧美精品亚洲一区二区| 国产精品二区激情视频| av片东京热男人的天堂| 国产精品亚洲一级av第二区| 欧美日韩中文字幕国产精品一区二区三区| svipshipincom国产片| 一区二区三区高清视频在线| 51午夜福利影视在线观看| 久久热在线av| 亚洲国产精品sss在线观看| 好男人在线观看高清免费视频 | 色播在线永久视频| 久久精品亚洲精品国产色婷小说| 一进一出好大好爽视频| 国产精品美女特级片免费视频播放器 | 国产高清激情床上av| 视频区欧美日本亚洲| 精品欧美一区二区三区在线| 国产精品综合久久久久久久免费| 久久中文字幕一级| 午夜久久久久精精品| 人人妻人人澡欧美一区二区| 18美女黄网站色大片免费观看| 免费观看精品视频网站| 草草在线视频免费看| 国产国语露脸激情在线看| 国产精品亚洲一级av第二区| 中文字幕另类日韩欧美亚洲嫩草| 日本一本二区三区精品| 午夜影院日韩av| 男人操女人黄网站| 成人永久免费在线观看视频| 欧美激情 高清一区二区三区| 两人在一起打扑克的视频| 亚洲人成网站在线播放欧美日韩| 琪琪午夜伦伦电影理论片6080| 国产高清videossex| 日韩成人在线观看一区二区三区| 青草久久国产| 岛国在线观看网站| 在线天堂中文资源库| 中出人妻视频一区二区| 成人av一区二区三区在线看| cao死你这个sao货| 在线国产一区二区在线| 国产一卡二卡三卡精品| 亚洲成a人片在线一区二区| 又紧又爽又黄一区二区| 91av网站免费观看| 免费在线观看日本一区| 成人三级黄色视频| 黄片大片在线免费观看| 中文字幕精品免费在线观看视频| 18禁国产床啪视频网站| 97超级碰碰碰精品色视频在线观看| 国产午夜福利久久久久久| 一个人观看的视频www高清免费观看 | 成人三级做爰电影| 亚洲精品在线观看二区| 啪啪无遮挡十八禁网站| 色播在线永久视频| 1024香蕉在线观看| 亚洲国产中文字幕在线视频| av天堂在线播放| 久久午夜亚洲精品久久| 黄色成人免费大全| 成人亚洲精品一区在线观看| 中文在线观看免费www的网站 | 波多野结衣av一区二区av| 99精品久久久久人妻精品| 午夜激情av网站| 男女之事视频高清在线观看| 男男h啪啪无遮挡| 88av欧美| 欧美一级毛片孕妇| www.自偷自拍.com| 午夜免费激情av| 白带黄色成豆腐渣| 俄罗斯特黄特色一大片| 窝窝影院91人妻| 欧美激情 高清一区二区三区| √禁漫天堂资源中文www| 精华霜和精华液先用哪个| 亚洲欧美激情综合另类| 99久久久亚洲精品蜜臀av| a在线观看视频网站| 国产亚洲av高清不卡| 日韩av在线大香蕉| 国产成+人综合+亚洲专区| 国产成年人精品一区二区| 国产成人欧美在线观看| 嫩草影视91久久| 亚洲人成电影免费在线| 动漫黄色视频在线观看| 一级片免费观看大全| 亚洲专区字幕在线| 亚洲狠狠婷婷综合久久图片| 啦啦啦韩国在线观看视频| 欧美黄色片欧美黄色片| 99热这里只有精品一区 | av视频在线观看入口| 国产成人精品久久二区二区免费| 99在线人妻在线中文字幕| 久久久久久久久免费视频了| 国产免费av片在线观看野外av| 欧美黑人精品巨大| 国产日本99.免费观看| 久久久久亚洲av毛片大全| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲一区二区三区色噜噜| 一区二区三区精品91| 激情在线观看视频在线高清| 国产成人精品久久二区二区91| 精品高清国产在线一区| 女生性感内裤真人,穿戴方法视频| 男女那种视频在线观看| 久久天堂一区二区三区四区| 国内久久婷婷六月综合欲色啪| 日韩精品免费视频一区二区三区| 国产视频内射| 久久人人精品亚洲av| 桃色一区二区三区在线观看| 人人妻人人澡欧美一区二区| aaaaa片日本免费| www.自偷自拍.com| 三级毛片av免费| 国产爱豆传媒在线观看 | 久9热在线精品视频| 妹子高潮喷水视频| 欧美乱妇无乱码| 麻豆久久精品国产亚洲av| 真人一进一出gif抽搐免费| 每晚都被弄得嗷嗷叫到高潮| 国产av又大| 黄色 视频免费看| 最近在线观看免费完整版| 久久久久久免费高清国产稀缺| 禁无遮挡网站| 黄色 视频免费看| 欧美色视频一区免费| 国产高清视频在线播放一区| 亚洲一卡2卡3卡4卡5卡精品中文| 国产高清视频在线播放一区| 欧美日韩福利视频一区二区| 国产高清视频在线播放一区| 无人区码免费观看不卡| 久久久久亚洲av毛片大全| 97超级碰碰碰精品色视频在线观看| 韩国精品一区二区三区| 亚洲人成伊人成综合网2020| 亚洲精品一区av在线观看| 亚洲免费av在线视频| 美女免费视频网站| 亚洲男人的天堂狠狠| 久久狼人影院| 在线播放国产精品三级| 国产蜜桃级精品一区二区三区| 国内久久婷婷六月综合欲色啪| 国产精品免费一区二区三区在线| 一本精品99久久精品77| 岛国视频午夜一区免费看| 啦啦啦观看免费观看视频高清| 国产一区二区在线av高清观看| 精品卡一卡二卡四卡免费| 欧美激情高清一区二区三区| 国产精品亚洲美女久久久| 一夜夜www| 亚洲av五月六月丁香网| 久久久久国内视频| 亚洲精品国产区一区二| 非洲黑人性xxxx精品又粗又长| 嫩草影视91久久| 国产亚洲欧美98| 国产主播在线观看一区二区| 精品久久久久久久久久久久久 | 少妇裸体淫交视频免费看高清 | 亚洲国产精品久久男人天堂| 午夜免费观看网址| 精品国产乱子伦一区二区三区| 久久久久久久午夜电影| 香蕉久久夜色| 一个人观看的视频www高清免费观看 | 妹子高潮喷水视频| 亚洲成人久久爱视频| 久久人人精品亚洲av| 午夜久久久在线观看| 99热只有精品国产| 久久国产亚洲av麻豆专区| 制服丝袜大香蕉在线| 欧美不卡视频在线免费观看 | 日韩高清综合在线| 麻豆久久精品国产亚洲av| 99久久精品国产亚洲精品| 18禁黄网站禁片免费观看直播| 午夜成年电影在线免费观看| 88av欧美| 韩国av一区二区三区四区| 国产成人精品久久二区二区91|