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

    基于博弈論的虛擬制造網(wǎng)絡(luò)車間調(diào)度優(yōu)化方法

    2019-07-11 11:32:40張國輝王小剛白躍偉
    中國機(jī)械工程 2019年12期
    關(guān)鍵詞:道工序納什實(shí)例

    聶 黎 張國輝 王小剛 白躍偉

    1.上海第二工業(yè)大學(xué)智能制造與控制工程學(xué)院,上海,201209 2.鄭州航空工業(yè)管理學(xué)院管理科學(xué)與工程學(xué)院,鄭州,450015

    0 引言

    隨著全球供應(yīng)鏈的迅速發(fā)展和廣泛應(yīng)用,越來越多的中小企業(yè)聯(lián)合構(gòu)建虛擬制造網(wǎng)絡(luò)(virtual manufacturing network, VMN)來生產(chǎn)單個(gè)企業(yè)無法完成的特定復(fù)雜產(chǎn)品[1]。構(gòu)建及運(yùn)行VMN的關(guān)鍵問題之一是對生產(chǎn)過程進(jìn)行科學(xué)管理,以較高的柔性和效率滿足客戶個(gè)性化需求。在作業(yè)調(diào)度決策過程中,需要考慮顧客和制造商的利益:來自不同顧客的訂單都應(yīng)盡快得以滿足;VMN應(yīng)以最高的效率和最低的成本完成所有加工任務(wù)。然而顧客和制造商雙方利益有時(shí)互相促進(jìn),有時(shí)互相沖突,因此有必要開發(fā)一種具有協(xié)調(diào)機(jī)制的作業(yè)調(diào)度方法以均衡雙方利益。

    博弈論是平衡多方利益的有效工具之一,近年來許多學(xué)者將之用于作業(yè)調(diào)度領(lǐng)域。TAYEB等[2]將博弈論原理用于帶機(jī)器故障維修計(jì)劃的車間作業(yè)調(diào)度問題。LI等[3]引入納什均衡方法來解決集成工藝規(guī)劃與調(diào)度問題,并設(shè)計(jì)了混合遺傳算法。ZHOU等[4]構(gòu)建了網(wǎng)絡(luò)化制造環(huán)境下作業(yè)調(diào)度問題的博弈論數(shù)學(xué)模型,并開發(fā)了基于遺傳算法的求解算法來求解該數(shù)學(xué)模型。周光輝等[5-6]針對服務(wù)型制造車間關(guān)鍵任務(wù)調(diào)度問題和服務(wù)型制造系統(tǒng)外包任務(wù)分配決策問題,開發(fā)了基于博弈論的Stackelberg模型。SUN等[7]、ZHANG等[8]應(yīng)用非合作博弈理論,建立了動(dòng)態(tài)生產(chǎn)環(huán)境下柔性作業(yè)車間調(diào)度問題的調(diào)度博弈模型。但這些研究都忽略了顧客和制造商雙方利益的均衡。本文針對虛擬制造網(wǎng)絡(luò)的車間調(diào)度問題,利用博弈論的納什均衡策略使各顧客和制造商VMN之間的利益達(dá)到均衡。

    1 虛擬制造網(wǎng)絡(luò)車間調(diào)度問題描述

    VMN車間調(diào)度問題可以表述如下:一個(gè)VMN由多個(gè)企業(yè)組成,每個(gè)企業(yè)提供不少于一臺(tái)的加工設(shè)備用于共享。這些共享的設(shè)備記作Mj(j= 0,1,…,m-1),其中,m為設(shè)備的臺(tái)數(shù)。該VMN需要加工一批來自不同客戶訂單的工件,記作Ji(i=0,1,…,n-1),其中,n為工件的個(gè)數(shù)。工件Ji包含hi道工序,記作Oi,l(l=0,1,…,hi-1)。每個(gè)工件的工序必須按照事先約定的順序進(jìn)行操作。工序Oi,l有ki,l個(gè)加工設(shè)備可供選擇,工序Oi,l的候選加工設(shè)備集合記作Mi,l。工序Oi,l在其第q個(gè)候選加工設(shè)備上的加工時(shí)間記作pi,l,q(q=0,1,…,ki,l-1)。工序Oi,l的最小加工時(shí)間記作pi,l,pi,l=min(pi,l,q)。本文假設(shè)VMN中所有加工設(shè)備都已準(zhǔn)備就緒且不發(fā)生故障。

    2 虛擬制造網(wǎng)絡(luò)車間調(diào)度博弈模型

    本文采用博弈論思想,將上述VMN車間調(diào)度問題建模為一個(gè)包含n+1個(gè)參與者的非合作博弈模型。該模型表述為三元組G=(n+1,S,U),其中,S為所有參與者的行動(dòng)策略集,S={s0,s1,…,sn},st表示第t個(gè)參與者的策略,t=0,1,…,n;U為所有參與者的收益函數(shù)集,U={u0,u1,…,un},ur為第r個(gè)參與者的收益函數(shù),r=0,1,…,n。

    2.1 參與者

    在該博弈模型中,VMN及需要加工的n個(gè)工件都視為博弈的參與者。VMN決定工件的加工順序,每個(gè)工件選擇自己的加工路徑。

    2.2 行動(dòng)策略

    該模型包含行動(dòng)策略不同的兩類參與者:工件及整個(gè)制造系統(tǒng)VMN。加工工件Ji的策略si=(mi,1,mi,2,…,mi,hi),其中,mi,l∈Mi,l。假設(shè)某工件有3道工序,且每道工序有3臺(tái)可供選擇的機(jī)器,那么該工件可能的行動(dòng)策略有27個(gè)。VMN的策略sn= (l1,l2,…,ln)表示n個(gè)加工工件之間的相對加工順序關(guān)系。

    2.3 收益函數(shù)

    博弈過程中,每個(gè)參與者根據(jù)自己的策略采取行動(dòng),從而獲得收益,并期望最大化自身利益。加工工件的目標(biāo)是希望以最快的速度完成加工,VMN的目標(biāo)是以最短的時(shí)間完成所有工件。

    加工工件的流程時(shí)間越短,收益越大,故加工工件收益函數(shù)設(shè)計(jì)為加工工件流程時(shí)間的遞減函數(shù),即加工工件Ji的收益函數(shù)為

    ui=ui(S)=1/fi

    (1)

    式中,fi為工件Ji的流程時(shí)間。

    VMN完成所有工件所耗時(shí)間makespan越短,收益就越大,故其收益函數(shù)設(shè)計(jì)為makespan的遞減函數(shù),即VMN的收益函數(shù)為

    un=un(S)=1/(ms)

    (2)

    每個(gè)參與者的收益不僅與自身的行動(dòng)策略有關(guān),還受其他參與者行動(dòng)策略的影響。有時(shí)它們之間的相互影響是積極的,若所有工件都能以最短的流程時(shí)間完成加工,則VMN的makespan最短。有時(shí)它們之間的相互影響是消極的,某個(gè)工件在追求自身的最短流程時(shí)間時(shí),可能要延遲其他工件的加工,從而引起其他工件的流程時(shí)間延長,有時(shí)甚至導(dǎo)致整個(gè)VMN的makespan延長。所以,需要找到所有參與者利益的平衡點(diǎn)。

    2.4 納什均衡

    (3)

    本文的n+1非合作博弈模型不一定存在純納什均衡。為了解決此問題,本文定義使得每個(gè)參與者都獲得最大收益的解為理想納什均衡。很明顯,如果博弈模型滿足

    fi=fi

    (4)

    則該博弈模型存在理想納什均衡,反之亦然。

    式(4)是個(gè)極強(qiáng)條件,大多數(shù)調(diào)度問題都不具備該性質(zhì)。本文針對不存在理想納什均衡的博弈模型定義一個(gè)解與理想納什均衡之間的距離:

    fi)+(ms-m)

    (5)

    可見,使得D越小的解越優(yōu)。本文把具有最小D的解稱為D最小納什均衡。如果一個(gè)解是理想納什均衡,那么它對應(yīng)的D為0。

    由此,本文將調(diào)度問題轉(zhuǎn)化為尋找D最小納什均衡問題。下一小節(jié)將具體闡述利用遺傳算法求解D最小納什均衡。在這此前,通過一個(gè)例子來理解本文的博弈模型。

    表1給出的調(diào)度實(shí)例考察了一個(gè)包含3個(gè)加工機(jī)器的VMN(假設(shè)時(shí)間單位為min)?,F(xiàn)有3個(gè)工件需要安排加工,工件J0只要完成工序O0,0和O0,1就可以交付。由于工序O0,0有2臺(tái)可選加工機(jī)器,工序O0,1有3臺(tái)可選加工機(jī)器,所以工件J0有6種不同的加工路徑策略。同理,工件J1有12種加工路徑策略,工件J2有6種加工路徑策略。3個(gè)工件有3!種不同加工順序關(guān)系,所以一共有6×12×6×6=2 592種可能的調(diào)度方案。該調(diào)度實(shí)例對應(yīng)的博弈模型有4個(gè)參與者,一共有2 592種不同的行動(dòng)策略組合。在不同的行動(dòng)策略組合下,4個(gè)參與者可能得到不同的收益。表2比較了2個(gè)不同行動(dòng)策略組合下的各參與者收益情況。由表2發(fā)現(xiàn),策略S1和S2的區(qū)別在于:根據(jù)策略S1,工序O2,0在機(jī)器M2上加工,而根據(jù)策略S2,工序O2,0在機(jī)器M0上加工;另外,3個(gè)工件加工順序也不同。不難發(fā)現(xiàn),策略S2是一個(gè)納什均衡,而S1不是。

    表1 調(diào)度實(shí)例的相關(guān)數(shù)據(jù)

    表2 兩種行動(dòng)策略下參與者收益比較

    注:表中的收益是按照式(1)或式(2)計(jì)算結(jié)果的倒數(shù)

    3 基于遺傳算法的D最小納什均衡求解算法

    本文設(shè)計(jì)了基于遺傳算法的優(yōu)化算法來尋找博弈模型的D最小納什均衡。算法流程如下:

    (1)設(shè)定算法的參數(shù)。如果式(4)成立,則該實(shí)例存在理想納什均衡,輸出理想納什均衡,算法終止;否則該實(shí)例不存在理想納什均衡,進(jìn)入下一步驟。

    (2)隨機(jī)生成初始種群,評(píng)估每個(gè)個(gè)體的適合度。

    (3)通過各種遺傳操作,如選擇、變異和交叉,由當(dāng)前種群演變成新一代種群。

    (4)判斷本次迭代是否滿足終止條件即是否達(dá)到一定的迭代次數(shù)。如果滿足條件,則算法終止,最佳個(gè)體即為所求;否則,轉(zhuǎn)到步驟(3),進(jìn)入下一次迭代。

    由于篇幅有限,僅側(cè)重于編碼與解碼機(jī)制以及適應(yīng)度函數(shù)的討論。

    3.1 編碼與解碼機(jī)制

    編碼是將調(diào)度問題可行解描述為遺傳算法“染色體”的過程?!叭旧w”包含的調(diào)度問題可行解的信息就是博弈模型中每個(gè)參與者的行動(dòng)策略。編碼設(shè)計(jì)中,一方面需要考慮如何在染色體中完整包含這些信息,另一方面還需要盡量以緊湊的方式表達(dá)這些信息,以減少存儲(chǔ)空間開銷,同時(shí)還要方便遺傳算法在進(jìn)化過程中對染色體進(jìn)行的各種遺傳算子操作。

    假設(shè)調(diào)度問題包括n個(gè)工件,則每個(gè)染色體被設(shè)計(jì)成n+1個(gè)部分。前n個(gè)部分對應(yīng)n個(gè)工件的行動(dòng)策略,即工件的各道工序所選擇的加工機(jī)器。最后一個(gè)部分對應(yīng)VMN,該部分染色體決定了所有工件的加工順序。以表1中的實(shí)例為例,((1,0), (1,2,0), (2,0), (0,2,1))就是一個(gè)染色體。該染色體由4個(gè)部分組成。第1部分(1,0)是第1個(gè)工件J0的行動(dòng)策略,它表達(dá)的意思是該工件的第1道工序O0,0在其候選加工設(shè)備集合中的機(jī)器M2上加工;該工件的第2道工序O0,1在其候選加工設(shè)備集合中的機(jī)器M0上加工。第2部分(1,2,0)是第2個(gè)工件J1的行動(dòng)策略,它表達(dá)的意思是該工件的第1道工序O1,0在其候選加工設(shè)備集合中的機(jī)器M1上加工;該工件的第2道工序O1,1在其候選加工設(shè)備集合中的機(jī)器M2上加工;該工件的第3道工序O1,2在其候選加工設(shè)備集合中的機(jī)器M1上加工。第3部分(2,0)是第3個(gè)工件J2的行動(dòng)策略。它表達(dá)的意思是該工件的第1道工序O2,0在其候選加工設(shè)備集合中的機(jī)器M2上加工;該工件的第2道工序O2,1在其候選加工設(shè)備集合中的機(jī)器M0上加工。最后一部分(0,2,1)是VMN給各個(gè)工件排序的策略,它表達(dá)的意思是第1個(gè)工件最先加工,其次是第3個(gè)工件,最后加工第2個(gè)工件。該編碼機(jī)制的優(yōu)點(diǎn)在于保持了染色體的線性表達(dá),也保證了各種基本變異和交叉操作能很容易地在染色體上實(shí)現(xiàn)而不會(huì)產(chǎn)生非法解。

    解碼是編碼的逆過程,把按照規(guī)則編寫的染色體翻譯成為可行的調(diào)度方案。以表1的實(shí)例來說明如何將染色體((1,0), (1,2,0), (2,0), (0,2,1))解碼為活動(dòng)調(diào)度方案。第1個(gè)工件J0的第1道工序O0,0在M2上加工并且在加工順序中處于第1位,所以工序O0,0在0~4 min內(nèi)在M2上加工。然后,J0的第2道工序O0,1在4~6 min內(nèi)在M0上加工。第2個(gè)工件J1的第1道工序O1,0選擇在M1上加工,由于此時(shí)M1上沒有其他工件搶占機(jī)器,所以工序O1,0在0~4 min內(nèi)在機(jī)器M1上加工。然后,J1的第2道工序O1,1選擇在M2上加工。然而,第3個(gè)工件J2的第1道工序O2,0此時(shí)也選擇在M2上加工。由于工件J2加工順序相對于工件J1較前,所以工序O2,0在4~6 min內(nèi)在機(jī)器M2上加工,而工序O1,1在6~7 min內(nèi)在機(jī)器M2上加工。剩余的解碼過程不再贅述。最終該染色體可解碼成的活動(dòng)調(diào)度方案如圖1所示。

    圖1 染色體對應(yīng)的甘特圖Fig.1 Gantt chart corresponding the chromosome

    3.2 適應(yīng)度函數(shù)

    遺傳算法評(píng)估每個(gè)個(gè)體的適應(yīng)度時(shí),適應(yīng)度函數(shù)起到非常重要的作用。本文設(shè)計(jì)的適應(yīng)度函數(shù)為

    m)]-1

    (6)

    前期仿真實(shí)驗(yàn)得到的合適遺傳算法參數(shù)如下:種群的大小為100,最大迭代次數(shù)為100,交叉和變異操作的概率分別為0.8和0.2。

    4 仿真實(shí)驗(yàn)及結(jié)果分析

    為了驗(yàn)證本文優(yōu)化調(diào)度方法的有效性,從相關(guān)文獻(xiàn)中選擇了若干具有代表性的benchmark問題進(jìn)行測試。首先,本文以文獻(xiàn)[4]中實(shí)例為例進(jìn)行了測試,該實(shí)例規(guī)模為6×6,不存在理想納什均衡。圖2為用傳統(tǒng)方法和用本文方法得到的優(yōu)化調(diào)度方案甘特圖。表3給出了兩種方法的結(jié)果,其中,S1表示傳統(tǒng)方法得到的最優(yōu)解所對應(yīng)的各個(gè)個(gè)體行動(dòng)策略,S2表示本文方法得到的D最小納什均衡所對應(yīng)的各個(gè)參與者行動(dòng)策略。當(dāng)參與者是工件時(shí),“加工機(jī)器/加工順序”欄是工件對應(yīng)的加工機(jī)器;當(dāng)參與者是V(表示VMN)時(shí),“加工機(jī)器/加工順序”欄是所有工件的加工順序。

    由表3可以看出,兩種方法都可以得到最優(yōu)的makespan=36。然而,傳統(tǒng)方法中并沒有考慮各個(gè)顧客的利益,即不能在滿足整個(gè)VMN生產(chǎn)成本最低的情況下,保證每個(gè)顧客訂單盡可能早地完工交付,而本文方法在保證不延長makespan的基礎(chǔ)上,使得工件J1、J3和J5比傳統(tǒng)方法分別提前2 min、2 min和1 min單位完工。

    圖2 在文獻(xiàn)[4]實(shí)例上的調(diào)度甘特圖Fig.2 Gantt chart for the instance in [4]

    行動(dòng)策略S1行動(dòng)策略S2參與者加工機(jī)器/加工順序收益加工機(jī)器/加工順序收益J0M0,M3,M0,M1,M3,M132M0,M3,M0,M2,M5,M132J1M1,M0,M5,M4,M2,M028M1,M0,M5,M4,M1,M026J2M4,M0,M5,M0,M4,M136M4,M0,M3,M0,M4,M136J3M3,M1,M2,M3,M1,M235M3,M1,M2,M3,M1,M233J4M2,M4,M3,M4,M2,M536M2,M1,M5,M4,M3,M536J5M5,M0,M1,M4,M0,M336M5,M2,M1,M2,M0,M335VJ0,J1,J2,J3,J4,J536J0,J1,J2,J3,J4,J536

    注:表中的收益是按照式(1)或式(2)計(jì)算結(jié)果的倒數(shù)

    另外,還選擇了文獻(xiàn)[9]中的10個(gè)實(shí)例進(jìn)行了測試。該10個(gè)實(shí)例規(guī)模從2×2到4×5,其中,實(shí)例SFJS07存在理想納什均衡,其他實(shí)例不存在理想納什均衡。表4列出了各個(gè)實(shí)例的結(jié)果。當(dāng)參與者是工件時(shí),“加工機(jī)器/加工順序”欄是工件對應(yīng)的加工機(jī)器;當(dāng)參與者是V(表示VMN)時(shí),“加工機(jī)器/加工順序”欄是所有工件的加工順序。從表4中結(jié)果可以看到,實(shí)例SFJS04和SFJS08不可能同時(shí)使VMN的makespan和各個(gè)工件的流程時(shí)間都達(dá)到最小,但本文方法可以適當(dāng)?shù)鼐釼MN和各個(gè)工件之間的利益。實(shí)例SFJS08不考慮各個(gè)工件的利益時(shí),VMN的最優(yōu)makespan是253 min;均衡各個(gè)工件利益時(shí),VMN的最優(yōu)makespan是256 min。雖然比傳統(tǒng)方法得到的makespan多3 min,但保證了每個(gè)工件都獲得盡可能小的流程時(shí)間,這意味著VMN以較小的成本代價(jià)使得所有顧客的要求盡早得以滿足。圖3、圖4為2種方法在實(shí)例SFJS04和SFJS08的甘特圖。

    表4 D最小納什均衡對應(yīng)的策略和收益

    注:表中的收益是按照式(1)或式(2)計(jì)算結(jié)果的倒數(shù)

    圖3 在實(shí)例SFJS04上的調(diào)度甘特圖Fig.3 Gantt chart for the instance of SFJS04

    圖4 在實(shí)例SFJS08上的調(diào)度甘特圖Fig.4 Gantt chart for the instance of SFJS08

    5 結(jié)語

    本文將虛擬制造網(wǎng)絡(luò)車間調(diào)度問題視為一場博弈,采用博弈論的概念和建模方法,提出了n+1非合作博弈調(diào)度優(yōu)化模型,并設(shè)計(jì)了基于遺傳算法的D最小納什均衡求解算法。對多個(gè)實(shí)例的測試驗(yàn)證了所提調(diào)度優(yōu)化方法是有效的,該方法能夠均衡多個(gè)顧客和VMN之間的利益,保證制造商生產(chǎn)效率和顧客滿意度的雙贏。進(jìn)一步的研究將使用博弈論的Stackelberg策略來開發(fā)虛擬制造網(wǎng)絡(luò)環(huán)境下車間調(diào)度問題的優(yōu)化方法,讓虛擬制造網(wǎng)絡(luò)中的某個(gè)企業(yè)占主導(dǎo)地位,其他企業(yè)作為跟隨者,為支撐生產(chǎn)經(jīng)營活動(dòng)良好運(yùn)作提供優(yōu)化的車間調(diào)度解決方案。

    猜你喜歡
    道工序納什實(shí)例
    “瓷中君子”誕生記
    例析求解排列組合問題的四個(gè)途徑
    THE ROLE OF L1 IN L2 LEARNING IN CHINESE MIDDLE SCHOOLS
    THE ROLE OF L1 IN L2 LEARNING IN CHINESE MIDDLE SCHOOLS
    修鐵鏈
    水泥各工序單位產(chǎn)品綜合電耗正確計(jì)算的實(shí)證研究
    四川水泥(2014年2期)2014-09-10 07:53:30
    完形填空Ⅱ
    完形填空Ⅰ
    師傅領(lǐng)進(jìn)門,修行靠個(gè)人
    愛比死亡更倔強(qiáng)
    幸福·悅讀(2007年8期)2007-08-14 04:36:58
    中文字幕人妻熟人妻熟丝袜美| 久9热在线精品视频| 国产一区二区三区在线臀色熟女| 女生性感内裤真人,穿戴方法视频| 深夜精品福利| 宅男免费午夜| 美女免费视频网站| 少妇人妻一区二区三区视频| 国产高潮美女av| 国产视频一区二区在线看| 国产免费男女视频| 国产成人av教育| 国产成人福利小说| 久久精品国产清高在天天线| 男人舔女人下体高潮全视频| 大型黄色视频在线免费观看| av天堂在线播放| 免费看a级黄色片| 亚洲美女黄片视频| 色在线成人网| 最好的美女福利视频网| 亚洲精品亚洲一区二区| 99国产极品粉嫩在线观看| 在线a可以看的网站| 麻豆成人av在线观看| 国产白丝娇喘喷水9色精品| 亚洲色图av天堂| 中文字幕av在线有码专区| 国内精品久久久久精免费| 高潮久久久久久久久久久不卡| 夜夜躁狠狠躁天天躁| 精品久久久久久久久久免费视频| 观看免费一级毛片| 淫妇啪啪啪对白视频| 成人美女网站在线观看视频| 一区二区三区免费毛片| 91麻豆精品激情在线观看国产| 国内揄拍国产精品人妻在线| 搡女人真爽免费视频火全软件 | 波多野结衣高清无吗| 欧美+亚洲+日韩+国产| 亚洲av成人不卡在线观看播放网| 老熟妇仑乱视频hdxx| 国产高潮美女av| 两个人的视频大全免费| 91久久精品电影网| 老熟妇仑乱视频hdxx| 国产精品综合久久久久久久免费| 欧美日韩黄片免| 亚洲成人久久性| 欧美一区二区精品小视频在线| 三级国产精品欧美在线观看| 亚洲一区二区三区色噜噜| 精品不卡国产一区二区三区| 久久国产精品影院| 国产乱人伦免费视频| 午夜福利在线在线| 国内揄拍国产精品人妻在线| 欧美性感艳星| 老熟妇乱子伦视频在线观看| 少妇的逼水好多| 麻豆一二三区av精品| 久久久久久久久久成人| 网址你懂的国产日韩在线| 一区二区三区激情视频| 免费在线观看影片大全网站| 此物有八面人人有两片| 最好的美女福利视频网| 最新中文字幕久久久久| 精品久久国产蜜桃| 露出奶头的视频| 夜夜爽天天搞| 国产亚洲av嫩草精品影院| 亚洲欧美日韩高清专用| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 在线观看舔阴道视频| 特级一级黄色大片| 成年版毛片免费区| 亚洲成人免费电影在线观看| 午夜日韩欧美国产| 亚洲av一区综合| 在线观看美女被高潮喷水网站 | 日韩亚洲欧美综合| 别揉我奶头~嗯~啊~动态视频| 国产中年淑女户外野战色| 极品教师在线视频| 国产又黄又爽又无遮挡在线| 日韩大尺度精品在线看网址| 色综合婷婷激情| 一个人看视频在线观看www免费| 国产白丝娇喘喷水9色精品| 午夜两性在线视频| 国产高潮美女av| 日韩成人在线观看一区二区三区| 国产主播在线观看一区二区| 色吧在线观看| 天天躁日日操中文字幕| 国产一区二区亚洲精品在线观看| 国产黄色小视频在线观看| 久久午夜福利片| 国产欧美日韩一区二区精品| 无人区码免费观看不卡| 老司机福利观看| 天堂影院成人在线观看| 麻豆一二三区av精品| 国产精品亚洲美女久久久| 成年版毛片免费区| 精品国产亚洲在线| 日韩欧美精品v在线| 国产一级毛片七仙女欲春2| 日本免费一区二区三区高清不卡| 久久香蕉精品热| 久久精品人妻少妇| 91av网一区二区| 欧美+亚洲+日韩+国产| 亚洲国产精品sss在线观看| 亚洲美女视频黄频| 久久精品影院6| 国产精品爽爽va在线观看网站| 一个人免费在线观看电影| 久久午夜福利片| 欧美色视频一区免费| 国产精品美女特级片免费视频播放器| 亚洲人成伊人成综合网2020| 亚洲 欧美 日韩 在线 免费| 中文亚洲av片在线观看爽| 不卡一级毛片| 村上凉子中文字幕在线| 久久这里只有精品中国| 久久99热这里只有精品18| 深夜精品福利| 国产成人福利小说| 黄色视频,在线免费观看| 我要看日韩黄色一级片| 免费人成视频x8x8入口观看| 亚洲精品久久国产高清桃花| 久久午夜亚洲精品久久| 99热6这里只有精品| 天堂av国产一区二区熟女人妻| 中文在线观看免费www的网站| 美女被艹到高潮喷水动态| 两性午夜刺激爽爽歪歪视频在线观看| 99久久九九国产精品国产免费| 国产亚洲精品久久久久久毛片| 一区二区三区高清视频在线| 一进一出抽搐动态| 亚洲成a人片在线一区二区| 成人鲁丝片一二三区免费| 日韩欧美在线二视频| 99精品在免费线老司机午夜| 91狼人影院| 亚洲精品一卡2卡三卡4卡5卡| 久久久久久九九精品二区国产| 国产一区二区在线观看日韩| 天堂√8在线中文| 亚洲av熟女| 日本 av在线| 男人舔女人下体高潮全视频| 国产三级黄色录像| 观看美女的网站| 亚洲精品久久国产高清桃花| 真人一进一出gif抽搐免费| 狂野欧美白嫩少妇大欣赏| 精品久久久久久成人av| 在线播放国产精品三级| 成人欧美大片| 免费看光身美女| 精品不卡国产一区二区三区| 中文字幕av成人在线电影| 五月伊人婷婷丁香| 国产成年人精品一区二区| 亚洲av一区综合| 一边摸一边抽搐一进一小说| 亚洲欧美激情综合另类| 99久久成人亚洲精品观看| 亚洲一区二区三区不卡视频| 精品乱码久久久久久99久播| 少妇的逼好多水| 两个人视频免费观看高清| 久久性视频一级片| 桃色一区二区三区在线观看| 亚洲一区高清亚洲精品| 久久精品国产亚洲av香蕉五月| av国产免费在线观看| 久久精品夜夜夜夜夜久久蜜豆| 变态另类丝袜制服| 欧美一区二区亚洲| 亚洲最大成人手机在线| 老司机午夜福利在线观看视频| av在线老鸭窝| 国产精品亚洲美女久久久| 欧美日韩黄片免| 18禁黄网站禁片免费观看直播| 亚洲国产欧洲综合997久久,| 国内毛片毛片毛片毛片毛片| 国产精品1区2区在线观看.| 亚洲精品一卡2卡三卡4卡5卡| 国产爱豆传媒在线观看| 国产主播在线观看一区二区| 亚洲性夜色夜夜综合| 国产极品精品免费视频能看的| 亚洲一区二区三区色噜噜| 国产午夜精品论理片| 国产熟女xx| 亚洲精品在线观看二区| 欧美黄色片欧美黄色片| 我要看日韩黄色一级片| 嫩草影院入口| 午夜激情福利司机影院| av在线老鸭窝| 人妻制服诱惑在线中文字幕| 在线观看av片永久免费下载| 亚洲精品成人久久久久久| 床上黄色一级片| 午夜免费激情av| 亚洲av电影不卡..在线观看| 啦啦啦观看免费观看视频高清| 免费看光身美女| 1000部很黄的大片| 国产私拍福利视频在线观看| 我要搜黄色片| 精品人妻熟女av久视频| 一本精品99久久精品77| 99国产精品一区二区三区| av在线蜜桃| 亚洲美女视频黄频| 18禁在线播放成人免费| 国产精品久久久久久人妻精品电影| 精品免费久久久久久久清纯| 国产精华一区二区三区| 男人和女人高潮做爰伦理| 欧美另类亚洲清纯唯美| 欧美日韩国产亚洲二区| 12—13女人毛片做爰片一| 成年女人毛片免费观看观看9| 老司机福利观看| 国产欧美日韩精品一区二区| 麻豆成人av在线观看| 少妇人妻一区二区三区视频| 国产欧美日韩精品亚洲av| 69av精品久久久久久| 亚洲人成电影免费在线| 国产精品久久久久久久久免 | 两性午夜刺激爽爽歪歪视频在线观看| 一区二区三区激情视频| 欧美bdsm另类| 日本a在线网址| 国产中年淑女户外野战色| 国产精品一区二区免费欧美| 又黄又爽又免费观看的视频| 麻豆久久精品国产亚洲av| 国产91精品成人一区二区三区| 免费大片18禁| 久9热在线精品视频| 精品人妻一区二区三区麻豆 | 日韩成人在线观看一区二区三区| 亚洲av日韩精品久久久久久密| 久久性视频一级片| 亚洲不卡免费看| 欧美日韩瑟瑟在线播放| 深夜a级毛片| 黄色视频,在线免费观看| 精品熟女少妇八av免费久了| 日韩av在线大香蕉| 亚洲国产精品成人综合色| 亚洲欧美激情综合另类| 怎么达到女性高潮| 全区人妻精品视频| 久久这里只有精品中国| 国产黄色小视频在线观看| 亚洲欧美激情综合另类| 国产欧美日韩精品一区二区| 精品一区二区三区人妻视频| 99热6这里只有精品| 亚洲美女搞黄在线观看 | 欧美一区二区亚洲| 狂野欧美白嫩少妇大欣赏| 丁香欧美五月| 丰满乱子伦码专区| 国产成人av教育| 人人妻人人看人人澡| 久久久久久九九精品二区国产| 波多野结衣高清作品| 欧美成狂野欧美在线观看| 国产又黄又爽又无遮挡在线| 日韩中文字幕欧美一区二区| 日韩欧美一区二区三区在线观看| 性插视频无遮挡在线免费观看| 女人被狂操c到高潮| 麻豆一二三区av精品| 一进一出抽搐gif免费好疼| 久久精品影院6| 久99久视频精品免费| 赤兔流量卡办理| 日本五十路高清| or卡值多少钱| 日本撒尿小便嘘嘘汇集6| 午夜日韩欧美国产| 亚洲美女黄片视频| 美女xxoo啪啪120秒动态图 | 91字幕亚洲| 国产成人福利小说| 别揉我奶头~嗯~啊~动态视频| 国产欧美日韩精品一区二区| 悠悠久久av| 久久精品影院6| 好男人在线观看高清免费视频| 怎么达到女性高潮| 日韩欧美国产一区二区入口| 午夜亚洲福利在线播放| 天堂网av新在线| 搡老熟女国产l中国老女人| 亚洲成人久久性| 国产亚洲欧美在线一区二区| 欧美色欧美亚洲另类二区| 国产精品电影一区二区三区| 一个人看的www免费观看视频| 国产av在哪里看| 成人特级av手机在线观看| 51国产日韩欧美| 亚洲国产日韩欧美精品在线观看| 亚洲真实伦在线观看| 成人特级av手机在线观看| 成人av在线播放网站| 99久久精品一区二区三区| 神马国产精品三级电影在线观看| 国产黄色小视频在线观看| 色综合欧美亚洲国产小说| 黄色日韩在线| 99久久精品热视频| 97热精品久久久久久| 亚洲在线自拍视频| 国产免费一级a男人的天堂| 国产一区二区激情短视频| 国产日本99.免费观看| 欧美日韩国产亚洲二区| av国产免费在线观看| 最近最新免费中文字幕在线| 国内久久婷婷六月综合欲色啪| 亚洲电影在线观看av| 哪里可以看免费的av片| 久久精品人妻少妇| 亚洲内射少妇av| 久久久国产成人免费| 国产精品野战在线观看| 国产精品久久久久久亚洲av鲁大| 琪琪午夜伦伦电影理论片6080| 久久草成人影院| 国语自产精品视频在线第100页| 亚洲美女视频黄频| 国产亚洲av嫩草精品影院| 搡老熟女国产l中国老女人| 日韩欧美精品v在线| 男女视频在线观看网站免费| 国产精品永久免费网站| 亚洲一区二区三区色噜噜| 变态另类丝袜制服| 最新在线观看一区二区三区| 色综合亚洲欧美另类图片| 十八禁人妻一区二区| 可以在线观看毛片的网站| 99国产极品粉嫩在线观看| 性色av乱码一区二区三区2| 色综合欧美亚洲国产小说| 久久国产精品影院| 亚洲自偷自拍三级| 精品熟女少妇八av免费久了| 国产熟女xx| 免费电影在线观看免费观看| 国产成人aa在线观看| 亚洲国产精品合色在线| 亚洲av.av天堂| 能在线免费观看的黄片| 国产亚洲精品av在线| 欧美bdsm另类| 午夜福利视频1000在线观看| 午夜精品一区二区三区免费看| 亚洲专区中文字幕在线| 免费大片18禁| 九色国产91popny在线| 97碰自拍视频| 69人妻影院| 中文在线观看免费www的网站| 国产爱豆传媒在线观看| 日本a在线网址| 男女那种视频在线观看| 神马国产精品三级电影在线观看| 成人午夜高清在线视频| 高潮久久久久久久久久久不卡| 亚洲乱码一区二区免费版| 日韩av在线大香蕉| 免费观看精品视频网站| 无遮挡黄片免费观看| 色精品久久人妻99蜜桃| 免费在线观看日本一区| 久99久视频精品免费| 内射极品少妇av片p| 欧美3d第一页| 直男gayav资源| 精品人妻偷拍中文字幕| a级毛片a级免费在线| 亚洲国产精品久久男人天堂| 成人毛片a级毛片在线播放| 天美传媒精品一区二区| 尤物成人国产欧美一区二区三区| 一进一出抽搐gif免费好疼| 欧美性感艳星| 国产aⅴ精品一区二区三区波| 国产免费一级a男人的天堂| 深爱激情五月婷婷| 天堂影院成人在线观看| 午夜福利高清视频| 日本a在线网址| 啦啦啦观看免费观看视频高清| 婷婷色综合大香蕉| 午夜福利视频1000在线观看| 日日摸夜夜添夜夜添av毛片 | 又粗又爽又猛毛片免费看| 日韩av在线大香蕉| 国产精品久久电影中文字幕| 又粗又爽又猛毛片免费看| 嫩草影院入口| 久久久久久久久中文| 久久久久久久久大av| 免费在线观看亚洲国产| 在线天堂最新版资源| 听说在线观看完整版免费高清| 亚洲天堂国产精品一区在线| 免费看a级黄色片| 嫁个100分男人电影在线观看| 在线观看舔阴道视频| 久久久久久久久大av| 国产人妻一区二区三区在| 免费无遮挡裸体视频| 国产精品日韩av在线免费观看| 国内久久婷婷六月综合欲色啪| 欧美日韩综合久久久久久 | 婷婷精品国产亚洲av| 琪琪午夜伦伦电影理论片6080| 国产精品永久免费网站| 日韩人妻高清精品专区| 国产黄a三级三级三级人| 天美传媒精品一区二区| 99国产精品一区二区蜜桃av| 中国美女看黄片| .国产精品久久| 精品不卡国产一区二区三区| 九九热线精品视视频播放| 国产 一区 欧美 日韩| 五月伊人婷婷丁香| 精品熟女少妇八av免费久了| 老熟妇乱子伦视频在线观看| 床上黄色一级片| 别揉我奶头 嗯啊视频| 国产野战对白在线观看| 午夜精品在线福利| 亚洲va日本ⅴa欧美va伊人久久| 亚洲内射少妇av| www.熟女人妻精品国产| 亚洲自拍偷在线| 亚洲人成电影免费在线| 精品不卡国产一区二区三区| 一夜夜www| 日日摸夜夜添夜夜添小说| 国产精品一及| 成熟少妇高潮喷水视频| 我的老师免费观看完整版| 岛国在线免费视频观看| 精品国内亚洲2022精品成人| 欧美bdsm另类| 久久99热6这里只有精品| 一区二区三区激情视频| 欧美成人性av电影在线观看| 人人妻人人看人人澡| 精品乱码久久久久久99久播| 久久99热6这里只有精品| 日韩欧美 国产精品| 91九色精品人成在线观看| 亚洲精品一区av在线观看| 真人一进一出gif抽搐免费| 亚洲精品一卡2卡三卡4卡5卡| 黄色视频,在线免费观看| 可以在线观看毛片的网站| 欧美午夜高清在线| 成熟少妇高潮喷水视频| 村上凉子中文字幕在线| 成年人黄色毛片网站| 黄色一级大片看看| 偷拍熟女少妇极品色| 男女那种视频在线观看| 黄色女人牲交| 色视频www国产| 91午夜精品亚洲一区二区三区 | 波野结衣二区三区在线| 精品免费久久久久久久清纯| a级一级毛片免费在线观看| 亚洲激情在线av| 桃色一区二区三区在线观看| 国产精品一区二区性色av| 高清毛片免费观看视频网站| 日本熟妇午夜| 变态另类丝袜制服| 国产色爽女视频免费观看| 久久欧美精品欧美久久欧美| 偷拍熟女少妇极品色| 中文字幕久久专区| 国产三级中文精品| 欧美在线一区亚洲| 亚洲无线观看免费| 哪里可以看免费的av片| 精品一区二区三区av网在线观看| 国内揄拍国产精品人妻在线| 久久草成人影院| 亚洲成人精品中文字幕电影| 一进一出抽搐动态| 国产高潮美女av| 久久性视频一级片| 国产精品久久久久久久电影| 亚洲aⅴ乱码一区二区在线播放| 日韩国内少妇激情av| 一个人免费在线观看的高清视频| 窝窝影院91人妻| 97超视频在线观看视频| 美女高潮的动态| 亚洲国产高清在线一区二区三| 色播亚洲综合网| 欧美绝顶高潮抽搐喷水| 欧美性感艳星| 在线国产一区二区在线| 国产欧美日韩精品亚洲av| 最近最新中文字幕大全电影3| 国产伦精品一区二区三区视频9| 亚洲自拍偷在线| 在现免费观看毛片| 国产亚洲精品久久久久久毛片| 欧美黑人巨大hd| 黄色配什么色好看| 国产高清有码在线观看视频| 内射极品少妇av片p| 国产伦人伦偷精品视频| 国产精品亚洲美女久久久| 90打野战视频偷拍视频| 日日夜夜操网爽| 国产精品一区二区三区四区久久| 精品一区二区三区av网在线观看| 国产精品久久久久久精品电影| 精品久久久久久久久久久久久| 久久精品国产亚洲av香蕉五月| 国产成人啪精品午夜网站| 热99在线观看视频| 内地一区二区视频在线| 久久这里只有精品中国| 蜜桃久久精品国产亚洲av| 久久99热这里只有精品18| 日韩欧美在线二视频| 日韩欧美一区二区三区在线观看| 欧美高清性xxxxhd video| 男人舔女人下体高潮全视频| 欧美另类亚洲清纯唯美| 色噜噜av男人的天堂激情| 精品日产1卡2卡| 男人舔女人下体高潮全视频| 校园春色视频在线观看| 国产一区二区在线观看日韩| 免费看日本二区| 亚洲成人精品中文字幕电影| 最新在线观看一区二区三区| 国产av在哪里看| 成年免费大片在线观看| www日本黄色视频网| 如何舔出高潮| 日韩中字成人| 亚洲精品乱码久久久v下载方式| 九九久久精品国产亚洲av麻豆| 高清在线国产一区| 国产伦精品一区二区三区视频9| 色尼玛亚洲综合影院| 国产精品免费一区二区三区在线| 亚洲三级黄色毛片| 丰满人妻一区二区三区视频av| 久久精品国产99精品国产亚洲性色| 深夜a级毛片| 蜜桃久久精品国产亚洲av| avwww免费| 色尼玛亚洲综合影院| 国产精品一及| 免费在线观看日本一区| 国产一区二区亚洲精品在线观看| 赤兔流量卡办理| 俄罗斯特黄特色一大片| 欧美绝顶高潮抽搐喷水| 两个人视频免费观看高清| 激情在线观看视频在线高清| 日本 欧美在线| 国产老妇女一区| 嫁个100分男人电影在线观看| 久久精品91蜜桃| 能在线免费观看的黄片| 国产爱豆传媒在线观看| 18禁黄网站禁片午夜丰满| 国产三级中文精品| 欧美3d第一页| 国产一区二区在线观看日韩| 精品人妻一区二区三区麻豆 | 男人的好看免费观看在线视频| 亚洲国产高清在线一区二区三| 真人一进一出gif抽搐免费| 每晚都被弄得嗷嗷叫到高潮| 免费在线观看成人毛片| av中文乱码字幕在线| 特大巨黑吊av在线直播| 男女下面进入的视频免费午夜| 精品人妻一区二区三区麻豆 | 一本一本综合久久|