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

    混合工作日歷下批量生產(chǎn)柔性作業(yè)車間多目標調度方法

    2022-01-20 07:35:22吳錦錦王培露李明金
    工業(yè)工程 2021年6期
    關鍵詞:道工序空閑時刻

    曾 強,吳錦錦,王培露,李明金

    (河南理工大學 1.能源科學與工程學院;2.圖書館,河南 焦作 454000)

    當前,批量生產(chǎn)已成為我國最主要的生產(chǎn)方式之一。批量生產(chǎn)作業(yè)車間調度對于提高生產(chǎn)率,降低生產(chǎn)成本,保證均衡生產(chǎn)等具有重要意義。隨著市場競爭的日益激烈,生產(chǎn)調度過程中需要綜合考慮多個目標。在同一個企業(yè)或車間里,不同型號、新舊不一而功能相同的設備共同存在,同一道工序往往能在多臺設備上加工。調查發(fā)現(xiàn),在一些企業(yè)或車間里不同設備有可能按不同的工作日歷運行。因此,受到多目標、設備柔性、混合工作日歷因素的綜合影響,調度人員往往無從尋優(yōu),主要依靠自身經(jīng)驗進行生產(chǎn)調度,導致調度效果較差,給企業(yè)帶來較大損失。因此,混合工作日歷下批量生產(chǎn)柔性作業(yè)車間調度成為一個重要研究方向。

    現(xiàn)有作業(yè)車間調度對工作日歷的處理方式主要有3種:完全不考慮工作日歷、調度中不考慮調度后考慮工作日歷、調度中考慮工作日歷。第1種處理方式在調度前后均不考慮工作日歷[1-3]。該方式計算量最小,但調度方案與實踐脫節(jié)。第2種處理方式在調度前不考慮工作日歷,調度后進行時間映射[4]。該方式計算量較小,當設備采用相同工作日歷運行時可采用這種方式。第3種處理方式在調度過程中即根據(jù)工作日歷準確推算設備調整和工序加工起止時刻。該方式雖計算量相對較大,但當設備采用混合工作日歷運行時,采用該方式才能得到有效調度方案。工作日歷是工作制和工作時段(或工作班次)的合稱。目前,關于調度中考慮工作日歷的作業(yè)車間調度研究成果很少,這些文獻對工作日歷的處理不夠全面,文獻[5-6]未考慮工作班次(工作時段),文獻[7]未考慮工作制。

    柔性作業(yè)車間多目標調度常用的求解方法有“間接法”和“直接法”兩種。間接法將多目標轉化為單目標求解,如極大極小法、線性加權和法[8]、主要目標法[9]、功效系數(shù)法[10]、分層排序法[11]等。該方法一次只能得到一個“最優(yōu)解”而不能得到最優(yōu)解集。直接法在多目標空間中尋找最優(yōu)解集,其中運用最多的是基于Pareto尋優(yōu)的方法,包括多目標遺傳算法[12]、強化Pareto遺傳算法[13]、小生境Pareto遺傳算法[14]、非支配排序遺傳算法[15]等。文獻[14, 16]通過研究已經(jīng)得出NSGA優(yōu)于其他算法的結論。隨著研究的不斷深入,第一代非支配排序遺傳算法(non-dominated sorting genetic algorithm,NSGA)在多目標優(yōu)化算法中的缺點不斷地暴露。為了更好地解決多目標優(yōu)化問題,學者提出了帶精英策略的非支配排序的遺傳算法(elitist non-dominated sorting genetic algorithm,NSGA-Ⅱ)。NSGA-Ⅱ降低了非劣排序遺傳算法的復雜性,具有運行速度快、解集的收斂性好的優(yōu)點。

    在批量生產(chǎn)方式下,生產(chǎn)周期的縮短,主要依賴于精細化調度。常用的措施有3種。1) 提前調整:將工序時間細分為設備調整時間和加工時間,若同一加工子批的前后2道工序所用設備不同,可使其下道工序提前進行設備調整。2) 交叉作業(yè):將同一個加工子批劃分為多個移動子批,采用平順移動方式移動,實現(xiàn)同一個加工子批前后工序的交叉作業(yè)。3) 擠壓安排:采用“間隙擠壓”法安排工序[17],最大程度地減少設備空閑時間。

    由于混合工作日歷下批量生產(chǎn)柔性作業(yè)車間多目標調度在柔性作業(yè)車間調度的基礎上進一步考慮了混合工作日歷、多目標優(yōu)化、批量生產(chǎn)等因素,因此是一個高度復雜的NP難問題。基于此,本文提出了一種混合工作日歷下批量生產(chǎn)柔性作業(yè)車間多目標調度方法以輔助生產(chǎn)調度。

    1 混合工作日歷下批量生產(chǎn)柔性作業(yè)車間多目標調度模型

    1.1 問題描述

    車間需在m臺設備上安排n批量產(chǎn)品(n個加工子批)的加工任務。假設如下:1) 各加工子批有多道工序,一道工序的可選加工設備可能不唯一;2) 各加工子批的工藝流程、設備調整時間、加工時間已知;3) 同一加工子批允許劃分為移動子批,采用平順移動方式移動,移動時間可忽略;4) 一個加工子批的加工不允許被其他加工子批打斷;5) 各加工子批優(yōu)先級相同;6) 在同一個調度周期內,各設備按設定的工作日歷(工作制和工作時段的合稱)運行,不允許改變;7) 當設備按其工作日歷啟停時,該設備上對應的工作同步啟停;8) 設備初始時間軸連續(xù);9) 調度從指定的調度起始時刻往后進行。在以上假設條件下,為調度人員提供一個Pareto調度集供其決策,優(yōu)化目標是生產(chǎn)周期最短,制造成本最低。

    1.2 符號說明

    表1給出了本模型涉及的符號。

    表1 符號定義Table1 Symbol definition

    1.3 模型構建

    1) 目標函數(shù)。

    2) 約束條件。

    其中,式(1)表示生產(chǎn)周期最小;式(2)表示制造成本最低;式(3)表示一道工序只能由m臺設備中的其中一臺完成;式(4)表示若Ji的工序j和工序j?1在同一臺設備上加工,則Ji的第一移動子批必須在工序j?1全部加工之后才能開始調整;式(5)表示若Ji的工序j和工序j?1不在同一臺設備上加工,則Ji的任意移動子批s的工序j可以提前進行調整,但必須在工序j?1完工后才可以開始加工;式(6)表示一臺設備不能同時加工兩道工序;式(7)表示任意工序的完工時刻與開始調整時刻之差不能小于其所需加工時間和調整時間之和;式(8)表示任意工序的設備加工開始時刻與加工結束時刻之間的有效時間與實際加工時間相等;式(9)表示任意工序的設備調整開始時刻與調整結束時刻之間的有效時間與實際調整時間相等;式(10)表示任意工序的設備調整開始時刻和結束時刻、任意移動子批加工開始時刻和結束時刻都在設備的工作時段內;式(11)表示任意工序的調整開始時刻不能早于調度起始時刻。

    2 NSGA II算法設計

    上述問題是一個高度復雜的NP難問題,一般的數(shù)學方法無法求解。鑒于NSGA-Ⅱ的優(yōu)點,以Excel VBA為平臺,設計NSGA-Ⅱ算法對其求解。

    2.1 定義對象、變量及數(shù)組

    根據(jù)算法需要,定義了圖1所示的自定義類型、表2所示的變量和數(shù)組。圖1中,Mb為調度前設備時間狀態(tài)數(shù)組,Mb(i).Ts為設備i的調度前時間狀態(tài)數(shù)組;J為加工子批數(shù)組,J(i)為加工子批i;Chr為種群數(shù)組,Chr(i)為個體i。

    圖1 自定義類型Figure1 Custom type

    2.2 算法流程

    本文算法流程[18]如圖2所示。

    圖2 算法流程Figure2 Flow of the algorithm

    2.3 讀取參數(shù)

    從工作表“設備”讀取參數(shù)賦給數(shù)組Mi;從工作表“加工子批”和“工藝流程”讀取參數(shù)賦給對象數(shù)組J;從工作表“工作制”讀取參數(shù)賦給數(shù)組Ws;按照假設條件8),令i從1~m變化,為Mb(i).Ts的每行賦2個元素,bt和tln;從工作表“算法參數(shù)”讀取參數(shù)賦給表2中的變量。

    表2 變量和數(shù)組Table2 Variables and arrays

    2.4 編碼方式

    采用“基于工序和設備的分段編碼”方式分別對工序和設備進行編碼。式(12)中,Chr(i).R為tpn×13的數(shù)組,第2、4列用于編碼,第1列為任務號,第3列為各加工子批的工序號,第5~13列用于解碼。第2列的自然數(shù)表示加工子批號,取值范圍為1~n,各自然數(shù)出現(xiàn)的次數(shù)與對應加工子批的工序數(shù)相等;第4列為各工序選擇的設備號,來自于該工序的可行加工設備號向量。式(12)表示加工子批1的第1道工序在設備4上加工(行2),第2道工序在設備1上加工(行4),加工子批2的第1道工序在設備2上加工(行1),第2道工序在設備3上加工(行3),依次類推。

    2.5 種群初始化

    對于個體Chr(i)的編碼獲取步驟如下。

    第1步,為Chr(i).R的第1列賦以從1~tpn的任務號;第2步,令j從1~n變化,在Chr(i).R第2列依次隨機尋找J(j).pn個空元素,賦為j,同時將對應工序號賦給Chr(i).R的同行第3列;第3步,令j從1~n變化,對每個j,產(chǎn)生1~J(Chr(i).R(j, 2)).Pr(Chr(i).R(j, 3)).an的隨機自然數(shù)k,令J(Chr(i).R(j, 2)).Pr(Chr(i).R(j, 4)=J(Chr(i).R(j, 2)).Pr(Chr(i).R(j, 3)).Mn(k)。易知,通過上述方式獲得的個體是可行的。

    令i從1~ps變化,分別對Chr(i).R賦值,完成種群初始化。

    2.6 遺傳操作

    遺傳操作包括交叉操作和變異操作。采用“基于工序和設備的分段”交叉和變異方式分別進行交叉和變異,并采用“遺傳算子改進策略”以保證交叉后子代個體可行[19-20]。

    2.7 解碼操作

    相對于傳統(tǒng)的柔性作業(yè)車間調度,本算法要同時考慮混合工作日歷、間隙擠壓、平順移動等因素,因此大大增加了解碼的難度。解碼操作的任務是對于個體Ch(i),計算Ch(i).R的第5~13列值,它們依次對應工序的移動子批對應的加工時間(第5列)、調整時間(第6列)、加工時間(第7列)、調整成本(第8列)、加工成本(第9列)、調整開始時刻(第10列)、調整結束時刻(第11列)、加工開始時刻(第12列)和加工結束時刻(第13列)。

    解碼操作的具體步驟如下。

    1) 令R=ch(i).R, Ma=Mb。

    2) 給第5~9列賦值。先獲取R(i, 4)在J(R(i, 2)).Pr(R(i, 3)).Mn數(shù)組中的索引號in, 再令R(i, 5)=J(R(i, 2)).Pr(R(i, 3)).Ct1(in);令R(i, 6)=J(R(i, 2)).Pr(R(i, 3)).St (in);令R(i, 7)=J(R(i, 2)).Pr(R(i, 3)).Ct (in);R(i, 8)=J(R(i, 2)).Pr(R(i, 3)).St(in)×J(R(i, 2)).Pr(R(i, 3)).Ps(in);R(i, 9)=J(R(i, 2)).Pr(R(i, 3)).Ct(in)×J(R(i, 2)).Pr(R(i, 3)).Pc (in)。

    3) 給第10~13列賦值。令i從1~n變化,對于任務i,操作如下。

    (1) 確定最早可能的調整開始時刻g。分3種情況處理。

    ① 若當前工序是加工子批R(i, 2)的第1道工序,令g=bt。

    ② 找到加工子批R(i, 2)的上一道工序所在的行號x,若R(x, 4)=R(i, 4)(圖3(a)),則 令g=max (t1,bt),表示若當前工序與上道工序所用設備相同,則須上道工序加工結束后才能開始當前工序的調整工作。

    ③找到加工子批R(i, 2)的上一道工序所在的行號x,若R(x, 4)≠R(i, 4)(圖3(b)),則令g=max (t2,t3,bt),t2取R(x, 13)逆 向 推 算[20]min (R(x, 5),R(i, 5))×(J(R(i, 2)).vn?1)+R(i, 6) 所得的時刻,t3取R(x, 12)正向推算[19](R(x, 5)所得的時刻。表示若當前工序與上道工序所用設備不相同,則按“平順移動原理”得到下一道工序最早可能的調整開始時刻。由于設備R(k, 4)與R(i, 4)所用的工作日歷有可能不同,導致t2有小于t3的可能,但上道工序的第一個移動子批未加工結束時,不可能開始下一道工序,故要取g=max (t2,t3, bt)。

    圖3 理論最早可調整開始時刻Figure3 Earliest start time of adjustment in theory

    (2) 從左向右,采用“間隙擠壓法”,根據(jù)Ma(R(i, 4)).Ts數(shù)組,依次尋找設備R(i, 4)的可行空閑時間段。如圖4所示,對于空閑時段k,取tb=max (g, Ma(R(i, 4)).Ts(2×k?1)),然后從tb正向推算[19]R(i, 6)+R(i, 7) 后得到加工結束時刻te,若te≤Ma(R(i,4)).Ts(2×k?1),表明當前工序能安排在設備R(i, 4)的第k個空閑時間段,否則判斷設備的下一個空閑時間段,直到找到可行空閑時間段為止。Ma(R(i, 4)).Ts(2×k?1))、Ma(R(i, 4)).Ts(2×k))分別是設備R(i, 4)的第k個空閑時間段的下限和上限。

    圖4 間隙擠壓法Figure4 Clearance extrusion method

    (3)安排當前任務。設找到設備R(i, 4)的第k個空閑時間段為可行時間段,令R(i, 10)=tb,R(i, 11)、R(i, 12)取tb正向推算R(i, 6)所得時刻,R(i, 13)取te。

    (4)更新設備R(i, 4)的時間狀態(tài)向量Ma(R(i, 4)).Ts。如圖5所示,調度開始時,Ma(R(i, 4)).Ts由bt、tln兩個元素組成,表示空閑時段為bt~tln(空閑時間段數(shù)frn=1)。若工序a被安排給設備R(i, 4),需將 tsba(設備調整開始時刻)和 t cea(工序加工結束時刻)插入到bt~tln之 間,此時空閑時間段分別為b t~tsba和 tcea~tln(空閑時間段數(shù)frn=2)。進一步,若工序b也被安排給設備R(i, 4)(排在工序a之后),此時,設備空閑時間段分別為 bt~tsba 、 tcea~tsbb 和 tceb~tln(空閑時間段數(shù)frn=3)?;诖?,Ma(R(i, 4)).Ts更新步驟如下。將Ma(R(i, 4)).Ts的長度增加2位,從設備R(i,4)的k空閑時段對應的第2個數(shù)據(jù)起,將其后的數(shù)據(jù)后移2位,空出2個空位,在這2個空位處分別填入待插入工序的設備調整開始時刻和加工結束時刻。

    圖5 更新Ma(R(i,4)).TsFigure5 Updating Ma(R(i,4)).Ts

    當所有任務安排完畢,令Ch(i).R=R,解碼結束。

    2.8 計算目標值

    在得到Ch(i).R的第5~13列之后,求出Ch(i).R第13列最大值與第10列最小值的差值賦給Chr(i).O(1),將第8、9列之和賦給Chr(i).O(2)。

    3 案例分析

    某機械加工車間擁有18臺設備,某次生產(chǎn)任務涉及到產(chǎn)品1~4,其工藝流程如表3所示。為縮短生產(chǎn)周期,調度前將4種產(chǎn)品進行批量劃分,形成8個加工子批,各加工子批的參數(shù)如表4所示。設備所用工作制共有X、Y、Z 3種,其設置如表5所示。表中每2列代表一個工作制。以工作制X為例,A列是非周末但不運行的日期,B列是周末但要運行的日期。各設備的工作日歷如表6所示?,F(xiàn)要合理安排8個加工子批的生產(chǎn)作業(yè)計劃,達到生產(chǎn)周期和生產(chǎn)成本最小化。

    表3 工藝流程Table3 Process flow

    表4 加工子批Table4 Processing sub-batch

    表5 工作制設置Table5 Settings of worksheet“work systems”

    表6 設備工作日歷Table6 Work calendars of machines

    取調度起始時刻bt =“2020/3/2 8:00”,ps=50,mg=100,cr =0.7,mr =0.3,mbs=2,tln=100 000(對應日期時間格式值4 637/11/26 0:00),移動批量均取5,多次獨立運算,得到的Pareto解集接近且均勻。圖6是Pareto解集之一,其中,Pareto解A對應的調度表如表7所示,加工子批甘特圖如圖7所示,設備甘特圖如圖8所示。

    圖6 Pareto解集Figure6 Pareto solution set

    圖7中每一行代表一個加工子批的生產(chǎn)作業(yè)計劃。圖8中每一行代表一個設備的生產(chǎn)排程,黑色框代表設備調整起止時刻,白色框代表工序加工起止時刻。從圖7和表7可見,同一加工子批的部分相鄰兩道工序之間存在時間上的交叉現(xiàn)象,這正是平順移動的結果。平順移動有利于下道工序盡早開工,從而有效縮短生產(chǎn)周期。

    圖7 Pareto解A對應的加工子批甘特圖Figure7 Gantt chart of processing sub-batch of Pareto solution A

    圖8 Pareto解A對應的設備甘特圖Figure8 Gantt chart of the machines of Pareto solution A

    由表7第10行可知,工序J4.4在設備1(200T)上加工,加工開始時刻為2020/3/6 14:30,加工結束時刻為2020/3/10 14:30,易知2020/3/6 14:30~2020/3/10 14:30的日歷時間為96 h,遠大于加工時間20 h。從表4可知,設備1采用“Z工作制”,每天運行10 h。從表5可知,按照“Z工作制”2020/3/7和2020/3/8為非運行日期,易知2020/3/6 14:30 ~ 2020/3/10 14:30的實際運行時間只有20 h,正好與工序的加工時間相等。同樣,由表7第16行知,工序J6.2在設備2(150T)上加工,其調整開始時刻為2020/3/7 21:18,結束時刻為2020/3/9 9:06,易知2020/3/7 21:18 ~ 2020/3/9 9:06的日歷時間為35.8 h,也遠遠大于其調整時間1.8 h。同理可得2020/3/7 21:18 ~ 2020/3/9 9:06的實際運行時間只有1.8 h,正好與設備的調整時間相等。綜上,算法在進行時間推算的時候考慮了設備的工作日歷,保證了所推算的設備調整起止時刻、工序加工起止時刻與實際情況吻合。

    為了比較移動批量對生產(chǎn)效率的影響,保持其他參數(shù)不變,分別取如下兩種移動批量進行進化計算。1) 30-30-20-20-20-20-20-20;2) 移動批量10-10-10-10-10-10-10-10。將得到的Pareto解集與移動批量5-5-5-5-5-5-5-5的Pareto解集進行對比,如圖9所示。

    圖9 移動批量變化對Pareto解集的影響Figure9 Impact to Pareto solution set by change of shift batch

    由圖9可見,在一定范圍內,移動批量越小,Pareto解集整體左移,即生產(chǎn)周期越短。這是因為平順移動下,移動批量的減小可以使單個移動批次的加工時間縮短。按照平順移動原則,下道工序可以提前更多的時間開始進行設備調整,最終使得生產(chǎn)周期縮短。然而,隨著移動批量的減少,移動次數(shù)有可能會增加。

    4 結論

    針對調度設備采用混合工作日歷的現(xiàn)象,提出了一種混合工作日歷下批量生產(chǎn)柔性作業(yè)車間多目標調度方法,主要結論如下。

    1) 算法能得到分布較均勻,收斂性較好的Pareto解集,能輔助調度人員快速找到滿意的調度方案。

    2) 算法在安排工序的過程中,采用了時間推算技術,根據(jù)各設備的工作日歷推算設備調整起止時刻和工序加工起止時刻,保證了調度方案的可行性。

    3) 加工批量不變,減小移動批量,通過平順移動生產(chǎn)周期能得到縮短,提高了調度方案的質量。

    猜你喜歡
    道工序空閑時刻
    恩賜
    詩選刊(2023年7期)2023-07-21 07:03:38
    冬“傲”時刻
    “瓷中君子”誕生記
    科學之友(2021年12期)2021-12-23 04:19:09
    例析求解排列組合問題的四個途徑
    捕獵時刻
    修鐵鏈
    “鳥”字謎
    小讀者之友(2019年9期)2019-09-10 07:22:44
    彪悍的“寵”生,不需要解釋
    WLAN和LTE交通規(guī)則
    CHIP新電腦(2016年3期)2016-03-10 14:09:48
    街拍的歡樂時刻到來了
    亚洲全国av大片| 禁无遮挡网站| 免费在线观看视频国产中文字幕亚洲| 国产人伦9x9x在线观看| 成年人黄色毛片网站| 精品国产亚洲在线| 青草久久国产| 日韩中文字幕欧美一区二区| 欧美国产日韩亚洲一区| 亚洲一卡2卡3卡4卡5卡精品中文| 久9热在线精品视频| 日本 欧美在线| 一级a爱视频在线免费观看| 波多野结衣高清作品| 欧美日韩亚洲综合一区二区三区_| 欧美激情高清一区二区三区| 国产麻豆成人av免费视频| 一夜夜www| 欧美精品啪啪一区二区三区| 午夜成年电影在线免费观看| 欧美日韩乱码在线| 夜夜躁狠狠躁天天躁| 母亲3免费完整高清在线观看| 搡老熟女国产l中国老女人| 国内毛片毛片毛片毛片毛片| 视频在线观看一区二区三区| 日日爽夜夜爽网站| 18美女黄网站色大片免费观看| 女同久久另类99精品国产91| 两个人视频免费观看高清| 色哟哟哟哟哟哟| 首页视频小说图片口味搜索| 日韩视频一区二区在线观看| 午夜福利18| 丰满人妻熟妇乱又伦精品不卡| 91av网站免费观看| 午夜福利18| 欧美黄色淫秽网站| 免费在线观看黄色视频的| 日本成人三级电影网站| 久久精品91蜜桃| 午夜福利视频1000在线观看| 亚洲精品在线美女| 国产精品野战在线观看| 国产成人一区二区三区免费视频网站| 色哟哟哟哟哟哟| 国产精品香港三级国产av潘金莲| 日本一区二区免费在线视频| 大型av网站在线播放| 色在线成人网| 久久久久久久久免费视频了| www.精华液| 老司机午夜福利在线观看视频| 色播亚洲综合网| 国产亚洲精品第一综合不卡| 精品不卡国产一区二区三区| 久99久视频精品免费| 国产av一区二区精品久久| 欧美精品亚洲一区二区| 1024手机看黄色片| 午夜激情av网站| 亚洲aⅴ乱码一区二区在线播放 | 男女视频在线观看网站免费 | 久久草成人影院| 制服丝袜大香蕉在线| 国产成人系列免费观看| 成年免费大片在线观看| 老鸭窝网址在线观看| 99国产极品粉嫩在线观看| 99久久99久久久精品蜜桃| 久久久久久久久久黄片| 国产精品电影一区二区三区| 91国产中文字幕| 啦啦啦免费观看视频1| 亚洲精品中文字幕在线视频| 亚洲色图av天堂| 一级毛片女人18水好多| 免费在线观看成人毛片| 亚洲人成伊人成综合网2020| 国产精品久久久av美女十八| 欧美乱色亚洲激情| 久久久久国产一级毛片高清牌| 制服丝袜大香蕉在线| 香蕉丝袜av| 一级a爱片免费观看的视频| 看黄色毛片网站| 男人的好看免费观看在线视频 | 999精品在线视频| 色av中文字幕| 精品熟女少妇八av免费久了| 日韩大码丰满熟妇| 88av欧美| 日韩精品青青久久久久久| 日韩国内少妇激情av| 人人妻人人澡人人看| 国产精品久久久久久人妻精品电影| 美女高潮到喷水免费观看| 麻豆成人av在线观看| 伦理电影免费视频| a级毛片a级免费在线| 亚洲 欧美 日韩 在线 免费| 男人操女人黄网站| 免费电影在线观看免费观看| av免费在线观看网站| 亚洲狠狠婷婷综合久久图片| 欧美色欧美亚洲另类二区| 久久久久久久精品吃奶| 少妇被粗大的猛进出69影院| 国产又色又爽无遮挡免费看| 老汉色∧v一级毛片| 亚洲五月色婷婷综合| 精品福利观看| 久久精品国产综合久久久| 国产精品永久免费网站| 国产精品久久久人人做人人爽| 国产精品99久久99久久久不卡| 免费看日本二区| 1024视频免费在线观看| 草草在线视频免费看| 成人国产一区最新在线观看| 成人亚洲精品av一区二区| 国产蜜桃级精品一区二区三区| 久久精品国产99精品国产亚洲性色| www.自偷自拍.com| 国产亚洲精品av在线| 成年女人毛片免费观看观看9| 亚洲性夜色夜夜综合| 国产午夜福利久久久久久| 亚洲一区中文字幕在线| 国产精品自产拍在线观看55亚洲| 午夜免费鲁丝| 亚洲av熟女| 美女大奶头视频| 午夜免费鲁丝| 一级a爱片免费观看的视频| 91老司机精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲国产毛片av蜜桃av| 成年人黄色毛片网站| 他把我摸到了高潮在线观看| 香蕉久久夜色| 欧美亚洲日本最大视频资源| 久久这里只有精品19| 久久久国产欧美日韩av| www国产在线视频色| 精品国产一区二区三区四区第35| 午夜激情福利司机影院| 草草在线视频免费看| 黑人巨大精品欧美一区二区mp4| 久久久水蜜桃国产精品网| 国产精品综合久久久久久久免费| 在线十欧美十亚洲十日本专区| 精品第一国产精品| 女人被狂操c到高潮| 丝袜美腿诱惑在线| 一a级毛片在线观看| 又紧又爽又黄一区二区| 精品久久久久久久久久免费视频| 免费在线观看影片大全网站| 亚洲av电影不卡..在线观看| 可以在线观看的亚洲视频| 国产一卡二卡三卡精品| 欧美日韩黄片免| 亚洲午夜精品一区,二区,三区| 两个人看的免费小视频| 国产精品av久久久久免费| 午夜福利在线在线| 不卡av一区二区三区| www.熟女人妻精品国产| 午夜免费激情av| 欧美性猛交黑人性爽| 国产不卡一卡二| 曰老女人黄片| svipshipincom国产片| 视频在线观看一区二区三区| 日本三级黄在线观看| 在线观看一区二区三区| 精品少妇一区二区三区视频日本电影| 无限看片的www在线观看| 久久午夜综合久久蜜桃| 欧美性猛交╳xxx乱大交人| 麻豆成人av在线观看| 国产黄色小视频在线观看| 国产精品,欧美在线| 老汉色∧v一级毛片| 亚洲专区字幕在线| 国产成+人综合+亚洲专区| 成年人黄色毛片网站| 久久99热这里只有精品18| netflix在线观看网站| 国产色视频综合| 国产三级在线视频| 欧美三级亚洲精品| 久久久久久久久免费视频了| 后天国语完整版免费观看| 日日干狠狠操夜夜爽| 亚洲中文av在线| 一区二区日韩欧美中文字幕| 村上凉子中文字幕在线| 看片在线看免费视频| 最近最新中文字幕大全电影3 | 免费在线观看日本一区| 精品国产乱子伦一区二区三区| 免费看日本二区| 一区二区日韩欧美中文字幕| 欧美日本亚洲视频在线播放| 91九色精品人成在线观看| 欧洲精品卡2卡3卡4卡5卡区| 黄片大片在线免费观看| 88av欧美| 亚洲中文日韩欧美视频| 亚洲天堂国产精品一区在线| 亚洲avbb在线观看| 午夜激情福利司机影院| 午夜精品在线福利| 国产一区二区激情短视频| 美女大奶头视频| 男人舔奶头视频| 天天躁夜夜躁狠狠躁躁| 99精品在免费线老司机午夜| 在线播放国产精品三级| 天堂影院成人在线观看| 国产久久久一区二区三区| 亚洲成av人片免费观看| 俺也久久电影网| 免费看美女性在线毛片视频| 91成人精品电影| 自线自在国产av| 精品一区二区三区视频在线观看免费| 亚洲欧美精品综合一区二区三区| 国产伦在线观看视频一区| 国产激情欧美一区二区| 淫妇啪啪啪对白视频| 亚洲欧洲精品一区二区精品久久久| 最近最新中文字幕大全电影3 | 亚洲自拍偷在线| 欧美亚洲日本最大视频资源| 午夜福利高清视频| 久久婷婷成人综合色麻豆| 久久久久久久久免费视频了| 中文字幕人妻熟女乱码| 欧美中文日本在线观看视频| 成年免费大片在线观看| 精品一区二区三区av网在线观看| 在线观看免费午夜福利视频| 国产99久久九九免费精品| 久久香蕉精品热| 看黄色毛片网站| 精品一区二区三区av网在线观看| 国产精品久久视频播放| 色婷婷久久久亚洲欧美| 狂野欧美激情性xxxx| 亚洲色图av天堂| 欧洲精品卡2卡3卡4卡5卡区| 成年人黄色毛片网站| 天堂√8在线中文| 亚洲国产高清在线一区二区三 | 啦啦啦观看免费观看视频高清| 国产极品粉嫩免费观看在线| 日韩精品免费视频一区二区三区| 欧美国产精品va在线观看不卡| 国产av一区在线观看免费| 在线观看午夜福利视频| 国内揄拍国产精品人妻在线 | 无遮挡黄片免费观看| 一边摸一边做爽爽视频免费| 日韩欧美国产在线观看| cao死你这个sao货| 国产在线观看jvid| 日韩一卡2卡3卡4卡2021年| 亚洲熟妇中文字幕五十中出| 久久久久亚洲av毛片大全| 日韩免费av在线播放| 国产精品香港三级国产av潘金莲| 国产91精品成人一区二区三区| 日韩 欧美 亚洲 中文字幕| 亚洲成人精品中文字幕电影| 韩国精品一区二区三区| 亚洲av五月六月丁香网| 99久久99久久久精品蜜桃| 99精品在免费线老司机午夜| 亚洲av中文字字幕乱码综合 | 亚洲精品美女久久久久99蜜臀| 精品国产亚洲在线| 亚洲专区国产一区二区| 亚洲电影在线观看av| 国产午夜福利久久久久久| 日韩免费av在线播放| 熟女电影av网| 免费无遮挡裸体视频| 一边摸一边做爽爽视频免费| 久久欧美精品欧美久久欧美| 久久久水蜜桃国产精品网| 黄色片一级片一级黄色片| 在线av久久热| 日韩国内少妇激情av| 国内揄拍国产精品人妻在线 | 欧美色欧美亚洲另类二区| 成人永久免费在线观看视频| 亚洲国产精品999在线| 色精品久久人妻99蜜桃| 国产99久久九九免费精品| 亚洲 国产 在线| 成在线人永久免费视频| 国内久久婷婷六月综合欲色啪| 天天躁狠狠躁夜夜躁狠狠躁| 夜夜爽天天搞| 怎么达到女性高潮| 精品国产一区二区三区四区第35| 狠狠狠狠99中文字幕| 久久国产精品人妻蜜桃| 97碰自拍视频| 亚洲黑人精品在线| 亚洲午夜精品一区,二区,三区| 他把我摸到了高潮在线观看| 中文字幕精品亚洲无线码一区 | 亚洲午夜精品一区,二区,三区| 美女高潮喷水抽搐中文字幕| 狂野欧美激情性xxxx| 欧美日韩精品网址| 亚洲av成人av| 午夜免费鲁丝| 18禁黄网站禁片免费观看直播| 久久精品国产亚洲av高清一级| 亚洲天堂国产精品一区在线| 日韩欧美一区二区三区在线观看| 99在线视频只有这里精品首页| 国产亚洲精品综合一区在线观看 | 欧美精品亚洲一区二区| 女性生殖器流出的白浆| 欧美色欧美亚洲另类二区| 日本一本二区三区精品| 久久婷婷成人综合色麻豆| 波多野结衣高清作品| 丝袜人妻中文字幕| 亚洲第一av免费看| 欧美日韩中文字幕国产精品一区二区三区| 男女那种视频在线观看| 人人妻人人澡欧美一区二区| 大型av网站在线播放| 国产伦在线观看视频一区| 美女午夜性视频免费| 大型av网站在线播放| 人人澡人人妻人| 亚洲专区中文字幕在线| 后天国语完整版免费观看| 波多野结衣巨乳人妻| 亚洲精品一区av在线观看| 成在线人永久免费视频| 亚洲专区字幕在线| 婷婷六月久久综合丁香| 18禁裸乳无遮挡免费网站照片 | 欧美激情久久久久久爽电影| 国产免费av片在线观看野外av| 成熟少妇高潮喷水视频| 国产成人欧美| 满18在线观看网站| 欧美一级a爱片免费观看看 | 久久久精品欧美日韩精品| 可以免费在线观看a视频的电影网站| 国产精品亚洲一级av第二区| 啦啦啦 在线观看视频| 久9热在线精品视频| 亚洲 欧美 日韩 在线 免费| 99国产综合亚洲精品| 首页视频小说图片口味搜索| 69av精品久久久久久| 少妇的丰满在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产欧洲综合997久久, | 国产成人系列免费观看| 女同久久另类99精品国产91| 美女大奶头视频| 欧美日韩中文字幕国产精品一区二区三区| 久久天躁狠狠躁夜夜2o2o| x7x7x7水蜜桃| 中文字幕av电影在线播放| 老司机在亚洲福利影院| 日本五十路高清| 日韩三级视频一区二区三区| 99riav亚洲国产免费| 国产97色在线日韩免费| 真人做人爱边吃奶动态| av欧美777| 亚洲黑人精品在线| 久久久久精品国产欧美久久久| 国产精品98久久久久久宅男小说| 哪里可以看免费的av片| 国产精品久久电影中文字幕| 亚洲自拍偷在线| av片东京热男人的天堂| 国语自产精品视频在线第100页| 欧美激情 高清一区二区三区| 日韩欧美免费精品| 热99re8久久精品国产| 波多野结衣巨乳人妻| 天天躁夜夜躁狠狠躁躁| 美女 人体艺术 gogo| 国产片内射在线| 悠悠久久av| a级毛片a级免费在线| 制服人妻中文乱码| 精品久久久久久,| 久久精品国产亚洲av高清一级| 手机成人av网站| 国产午夜精品久久久久久| 国产高清videossex| 18禁国产床啪视频网站| 可以免费在线观看a视频的电影网站| 亚洲精品在线美女| 极品教师在线免费播放| 国产久久久一区二区三区| 悠悠久久av| 一边摸一边抽搐一进一小说| 亚洲专区国产一区二区| 欧美zozozo另类| 看免费av毛片| 亚洲av成人不卡在线观看播放网| 亚洲欧美激情综合另类| 亚洲一区二区三区色噜噜| 午夜福利免费观看在线| 精品国内亚洲2022精品成人| av在线天堂中文字幕| 久久久久国产一级毛片高清牌| 亚洲午夜精品一区,二区,三区| 成人免费观看视频高清| 两性午夜刺激爽爽歪歪视频在线观看 | 国产极品粉嫩免费观看在线| 色婷婷久久久亚洲欧美| 国产精品自产拍在线观看55亚洲| 欧美中文综合在线视频| 黄色a级毛片大全视频| 国产主播在线观看一区二区| 桃色一区二区三区在线观看| 亚洲片人在线观看| 欧美黄色淫秽网站| 久久久久久人人人人人| x7x7x7水蜜桃| 老熟妇仑乱视频hdxx| 两个人免费观看高清视频| 国产国语露脸激情在线看| 精品少妇一区二区三区视频日本电影| 少妇裸体淫交视频免费看高清 | 可以在线观看的亚洲视频| 一本综合久久免费| 亚洲成人久久性| 亚洲国产高清在线一区二区三 | 久久香蕉激情| 性色av乱码一区二区三区2| 最近在线观看免费完整版| 校园春色视频在线观看| 夜夜看夜夜爽夜夜摸| 免费看十八禁软件| 精品久久蜜臀av无| 欧美国产精品va在线观看不卡| 欧美av亚洲av综合av国产av| 99久久无色码亚洲精品果冻| 法律面前人人平等表现在哪些方面| 欧美绝顶高潮抽搐喷水| 成人国产一区最新在线观看| 国产91精品成人一区二区三区| av在线天堂中文字幕| 国产成人av教育| 色av中文字幕| 国产精品99久久99久久久不卡| 婷婷六月久久综合丁香| 亚洲五月婷婷丁香| 亚洲av电影不卡..在线观看| 国产黄片美女视频| 国产精品免费一区二区三区在线| 桃红色精品国产亚洲av| 夜夜躁狠狠躁天天躁| 国产爱豆传媒在线观看 | 欧美性长视频在线观看| 天堂动漫精品| svipshipincom国产片| 亚洲 欧美 日韩 在线 免费| 三级毛片av免费| 亚洲成人免费电影在线观看| 国产av一区二区精品久久| 18禁黄网站禁片免费观看直播| 黄网站色视频无遮挡免费观看| 男女做爰动态图高潮gif福利片| 欧美黑人巨大hd| 在线十欧美十亚洲十日本专区| 亚洲av熟女| 一本大道久久a久久精品| av有码第一页| 日日夜夜操网爽| АⅤ资源中文在线天堂| xxx96com| 黄片小视频在线播放| 一级毛片高清免费大全| 99在线人妻在线中文字幕| xxx96com| 久久精品夜夜夜夜夜久久蜜豆 | 香蕉丝袜av| 精品久久久久久久人妻蜜臀av| 美女高潮喷水抽搐中文字幕| 精品久久久久久久久久免费视频| 天天添夜夜摸| 国产亚洲精品综合一区在线观看 | 欧美激情极品国产一区二区三区| 精品久久久久久久久久免费视频| 久久精品夜夜夜夜夜久久蜜豆 | 麻豆国产av国片精品| 国产视频一区二区在线看| 久久久久久久精品吃奶| 亚洲精品中文字幕在线视频| 男人的好看免费观看在线视频 | 午夜两性在线视频| 啪啪无遮挡十八禁网站| 国产精品一区二区精品视频观看| 国产真实乱freesex| 精品国产超薄肉色丝袜足j| 中亚洲国语对白在线视频| 婷婷亚洲欧美| 久久九九热精品免费| 亚洲 欧美 日韩 在线 免费| 久久午夜亚洲精品久久| 国产精品1区2区在线观看.| 国产黄色小视频在线观看| 91在线观看av| 久久精品国产亚洲av香蕉五月| 精品国产乱子伦一区二区三区| 亚洲精品色激情综合| 黄片播放在线免费| 久久人妻福利社区极品人妻图片| 看黄色毛片网站| 久久久国产精品麻豆| 高潮久久久久久久久久久不卡| 久久久国产精品麻豆| 国产1区2区3区精品| 午夜福利一区二区在线看| 国产日本99.免费观看| 久久精品国产99精品国产亚洲性色| 欧美+亚洲+日韩+国产| 精品熟女少妇八av免费久了| 麻豆一二三区av精品| 国产成年人精品一区二区| 婷婷丁香在线五月| 亚洲成人精品中文字幕电影| 欧美一级a爱片免费观看看 | 97人妻精品一区二区三区麻豆 | 国产单亲对白刺激| 亚洲精品色激情综合| 一本精品99久久精品77| 午夜福利欧美成人| 国产激情偷乱视频一区二区| 国产欧美日韩精品亚洲av| 黄色a级毛片大全视频| 日日干狠狠操夜夜爽| 女性被躁到高潮视频| 男人操女人黄网站| 午夜视频精品福利| 婷婷丁香在线五月| 久久精品国产清高在天天线| 亚洲专区字幕在线| 亚洲成av人片免费观看| 国产成人欧美| 国产成人啪精品午夜网站| 日韩欧美一区二区三区在线观看| 国产激情偷乱视频一区二区| 老鸭窝网址在线观看| 国产亚洲精品久久久久久毛片| 亚洲色图 男人天堂 中文字幕| 国产精品98久久久久久宅男小说| 日韩一卡2卡3卡4卡2021年| 久久精品国产清高在天天线| 国产高清videossex| 嫩草影院精品99| 欧美日韩瑟瑟在线播放| 99久久国产精品久久久| 午夜免费成人在线视频| 亚洲在线自拍视频| 视频区欧美日本亚洲| 麻豆成人av在线观看| 久久午夜亚洲精品久久| 亚洲片人在线观看| 最近最新中文字幕大全电影3 | 1024手机看黄色片| 18禁黄网站禁片免费观看直播| 国产精品野战在线观看| 日韩高清综合在线| 成人三级黄色视频| 亚洲免费av在线视频| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 黄频高清免费视频| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩一级在线毛片| 久久精品aⅴ一区二区三区四区| 久久午夜综合久久蜜桃| 国产蜜桃级精品一区二区三区| 国产在线精品亚洲第一网站| 亚洲中文av在线| 一本综合久久免费| 亚洲人成电影免费在线| 亚洲美女黄片视频| 成人欧美大片| 日韩欧美免费精品| 亚洲专区中文字幕在线| 国产麻豆成人av免费视频| 91成人精品电影| 免费在线观看日本一区| 亚洲九九香蕉| 国产一区二区三区视频了| 亚洲av熟女| 国产av在哪里看| 淫妇啪啪啪对白视频| 欧美日韩亚洲国产一区二区在线观看| 欧美日韩中文字幕国产精品一区二区三区| 日韩中文字幕欧美一区二区| 国产v大片淫在线免费观看| 99热只有精品国产|