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

    一種用于PFSP節(jié)能優(yōu)化的混合禁忌搜索算法

    2021-01-12 13:33:04張雨晨熊福力
    計(jì)算機(jī)測(cè)量與控制 2020年12期
    關(guān)鍵詞:交貨期搜索算法凈利潤(rùn)

    張雨晨, 熊福力

    (西安建筑科技大學(xué) 信息與控制工程學(xué)院,西安 710055)

    0 引言

    隨著綠色制造的全球化與生產(chǎn)力的提高,能源需求量日益增加,能源消耗造成的環(huán)境問(wèn)題日益突出[1]。通過(guò)節(jié)能與環(huán)境保護(hù)的手段可有效地推進(jìn)可持續(xù)發(fā)展政策,具有顯著的經(jīng)濟(jì)效益和社會(huì)效益。為了滿足綠色生產(chǎn)和客戶滿意度的要求,企業(yè)必須同時(shí)考慮訂單接受和生產(chǎn)計(jì)劃。訂單接受與調(diào)度(OAS,order acceptance and scheduling)問(wèn)題決定接受哪些訂單并決策接受訂單的加工順序,根據(jù)文獻(xiàn)[2-3]的研究,它有著與大多數(shù)車間調(diào)度問(wèn)題相同的特點(diǎn)。置換流水車間調(diào)度問(wèn)題(PFSP,permutation flow shop problem)是流水車間調(diào)度問(wèn)題中一類典型的子問(wèn)題,也是調(diào)度問(wèn)題的研究熱點(diǎn)[4]。它通常側(cè)重于優(yōu)化生產(chǎn)效率指標(biāo),如總完成時(shí)間、拖期、加權(quán)完工時(shí)間等[5-6]。而能源效率優(yōu)化是企業(yè)實(shí)現(xiàn)可持續(xù)發(fā)展的關(guān)鍵因素之一,因此,進(jìn)行訂單選擇與生產(chǎn)調(diào)度的統(tǒng)一決策,使企業(yè)能夠獲得更加穩(wěn)定有效的決策方案。車間能效優(yōu)化相關(guān)研究已有杰出成果,文獻(xiàn)[7-8]提出同時(shí)優(yōu)化能源消耗與總完工時(shí)間的多目標(biāo)數(shù)學(xué)規(guī)劃模型。劉向[9]提出一種混合遺傳算法求解面向節(jié)能降耗的調(diào)度模型。Shabtay[10-11]提出了集成提前/拖期懲罰和交貨期分配的數(shù)學(xué)模型。優(yōu)化PFSP的方法有很多,如精確求解算法、構(gòu)造啟發(fā)式算法與智能優(yōu)化算法。其中智能優(yōu)化算法更適合優(yōu)化PFSP,尤其是禁忌搜索[12]、遺傳算法、迭代局部搜索、粒子群算法等。

    本文突破了傳統(tǒng)PFSP的角度,建立了置換流水車間訂單接受與調(diào)度(PFSOAS,permutation flow shop order acceptance and scheduling)問(wèn)題,旨在最大化企業(yè)生產(chǎn)總凈利潤(rùn)(TNR,total net revenue)。本文于傳統(tǒng)禁忌搜索的基礎(chǔ)上,提出了一種節(jié)能混合禁忌搜索算法(EHTS,energy-saving hybrid tabu search),它整合了一種能耗調(diào)整策略與交貨期配置策略,實(shí)現(xiàn)對(duì)能耗的二次降低,加入了一種訂單接受與拒絕(OAR,order acceptance and rejection)的編碼方式,并設(shè)計(jì)了NEH構(gòu)造啟發(fā)式算法產(chǎn)生禁忌搜索算法的初始解,以改善算法優(yōu)化質(zhì)量。通過(guò)上述改進(jìn),提高了算法解決本文所提出的特殊的PFSOASP的尋優(yōu)能力。通過(guò)大量實(shí)際算例的統(tǒng)計(jì)實(shí)驗(yàn),分析了該算法的優(yōu)異性能,接著描述了所提出的PFSOAS模型,并且通過(guò)實(shí)驗(yàn)分析了EHTS算法求解PFSOAS問(wèn)題的有效性與優(yōu)越性。

    1 問(wèn)題描述

    由客戶提供一批訂單集合I,與訂單i相關(guān)的指標(biāo)分別是固定收入Ri,相同的截止日期D與交貨期d,拖期懲罰權(quán)重wi,客戶期望的交貨期區(qū)間。PFSOAS問(wèn)題通過(guò)決定接受那些訂單,決策已接受訂單的加工順序與確定訂單開(kāi)始加工時(shí)間使總利潤(rùn)最大化。

    1.1 相關(guān)假設(shè)與參數(shù)說(shuō)明

    相關(guān)假設(shè):每臺(tái)機(jī)器上每個(gè)訂單的處理時(shí)間為已知且恒定的;每個(gè)階段在同一時(shí)刻只能處理一個(gè)訂單;生產(chǎn)過(guò)程中無(wú)搶占,即某訂單只有在上一個(gè)訂單完工后才可被加工;任何接受的訂單只有當(dāng)前階段的完成處理后才能在下一階段處理,準(zhǔn)備時(shí)間為0;工作臺(tái)之間緩沖區(qū)無(wú)限大;每臺(tái)機(jī)器在調(diào)度開(kāi)始時(shí)處于空閑狀態(tài),并且一次最多可處理一個(gè)接受的訂單;機(jī)器上訂單的建立時(shí)間和釋放時(shí)間忽略不計(jì)。

    表1 模型參數(shù)說(shuō)明

    1.2 PFSOAS模型

    目標(biāo)函數(shù):

    maxTNR

    (1)

    (2)

    其中:TEC為對(duì)接受訂單進(jìn)行加工的總能源成本,ECi為加工訂單i時(shí)所需的機(jī)器能耗。

    (3)

    (4)

    s.t.:

    (5)

    2≤i≤n,2≤k≤m

    (6)

    (7)

    si,k=max(ci,(k-1),c(i-1),k), 2≤i≤n,2≤k≤m

    (8)

    ci,k=si,k+pi,k·xi

    ?i∈I,k∈K

    (9)

    Ci=ci,m,?i∈Io

    (10)

    (11)

    (12)

    ci,k≥ci,(k-1)+xipi,k,?i∈I,k∈K

    (13)

    ci,(k-1)≤si,k,?i∈Io,k∈K

    (14)

    c(i-1),k≤si,k,?i∈Io,k∈K

    (15)

    (16)

    (17)

    βi,k,t+γi,k,t≤1,?i∈I,k∈K,t∈T

    (18)

    yi,k,t=γi,k,t,?i∈I,k∈K,t∈T

    (19)

    式(5)計(jì)算拖期懲罰權(quán)重。式(6)與式(7)計(jì)算機(jī)器空閑時(shí)間。式(8)與式(9)分別計(jì)算開(kāi)始時(shí)間與完工時(shí)間。式(10)表示訂單i的總完工時(shí)間。式(11)表示生產(chǎn)效率的約束。式(12)表示交貨日期的約束。式(13)表示只有在第k-1階段處理了接受的訂單后,才能在第k階段處理接受的訂單。式(14)和式(15)表示任務(wù)不可中斷和機(jī)器不可被搶占,式(16)表示對(duì)訂單集中每個(gè)訂單進(jìn)行一次接受與拒絕,式(17)確保每個(gè)階段的所有訂單僅處理一次,式(18)確保每臺(tái)機(jī)器一次只能處理一個(gè)訂單,式(19)表示若γi,k,t=0,則第k個(gè)機(jī)器在t時(shí)刻關(guān)閉;否則γi,k,t=1。

    2 節(jié)能混合禁忌搜索算法

    本文提出的EHTS算法以NEH構(gòu)造啟發(fā)式算法產(chǎn)生初始解,提出了一種OAR編碼方式,同時(shí)設(shè)計(jì)了一種能耗調(diào)整策略,在保持加工所有接受訂單的最大完工時(shí)間不變下,使算法在求解的過(guò)程中利用分時(shí)電價(jià)的特點(diǎn),實(shí)現(xiàn)對(duì)能耗的二次降低,與一種交貨期配置策略相結(jié)合,實(shí)現(xiàn)TNR的提高。

    2.1 EHTS算法流程

    根據(jù)以上描述,EHTS算法總體為三階段優(yōu)化過(guò)程,第一階段用于優(yōu)化總凈利潤(rùn),第二階段作為節(jié)能調(diào)整階段,降低能耗成本,第三階段配置最佳交貨期,使總凈利潤(rùn)得到上升的空間。EHTS算法流程如圖1所示。

    圖1 EHTS算法流程圖

    針對(duì)相同交貨期下FPFSOAS問(wèn)題的EHTS算法流程描述如下:

    1)以總凈利潤(rùn)作為目標(biāo)值,采用NEH啟發(fā)式算法生成初始解,初始化禁忌表;

    2)對(duì)步驟1)后產(chǎn)生的初始解進(jìn)行禁忌搜索,每次迭代計(jì)算目標(biāo)值時(shí)使用訂單接受與拒絕方法重新編碼,更新局部最優(yōu)解和禁忌表;

    3)檢驗(yàn)是否滿足終止準(zhǔn)則,若滿足則終止算法迭代,解碼得到相應(yīng)解決方案,否則返回步驟2)進(jìn)行下一次迭代;

    4)解碼并保存最大總凈利潤(rùn)近優(yōu)解,并計(jì)算該解決方案中每個(gè)訂單每道工序的完工時(shí)間;

    5)根據(jù)能耗調(diào)整策略,對(duì)高峰期兩邊的工序加工時(shí)間進(jìn)行二次調(diào)整。

    6)在交貨期有效區(qū)間內(nèi)進(jìn)行交貨期配置。

    2.2 節(jié)能調(diào)整

    第一階段提出的基于OAR編碼方式與NEH產(chǎn)生初始解的改進(jìn)禁忌搜索方法,能夠以較低的耗電成本生成最大總凈利潤(rùn)的調(diào)度,但某些訂單仍會(huì)產(chǎn)生較高的電能消耗。使用能耗調(diào)整策略調(diào)整了第一階段優(yōu)化生成的解決方案,從而避免某些訂單的生產(chǎn)仍處于用電高峰期。

    調(diào)整規(guī)則:從最后一個(gè)接受的訂單的最后一道工序開(kāi)始,并以逆序的方式從最后一個(gè)訂單調(diào)整到從用電高峰期開(kāi)始生產(chǎn)的第一個(gè)訂單,使用這種調(diào)整策略對(duì)非高峰期的所有訂單的加工順序重新調(diào)度,并將這些工序從高峰時(shí)段移開(kāi),直到非高峰期的所有訂單完成調(diào)度。節(jié)能調(diào)整流程如圖2所示。

    圖2 節(jié)能調(diào)整流程圖

    2.3 交貨期配置策略

    本文以最大化考慮能耗成本與拖期懲罰的總凈利潤(rùn)目標(biāo),組成總凈利潤(rùn)的兩個(gè)指標(biāo);能耗成本與拖期懲罰均與時(shí)間有關(guān),企業(yè)可以根據(jù)總凈利潤(rùn)的優(yōu)化情況,向客戶提交不超過(guò)最晚交貨期的交付時(shí)間,因此交貨期配置策略會(huì)對(duì)總凈利潤(rùn)的再次提高具有一定很大的可能性。交貨期配置策略從上一階段已經(jīng)調(diào)整好的非高峰期生產(chǎn)的第一道工序開(kāi)始,保持非高峰期時(shí)段的工序不變,整體在交貨期區(qū)間內(nèi)移動(dòng),根據(jù)設(shè)置好的移動(dòng)時(shí)間長(zhǎng)度,尋找到使得總凈利潤(rùn)增長(zhǎng)的交貨期。交貨期配置策略的流程如圖3所示。

    圖3 交貨期配置策略

    2.4 編碼方式

    訂單接受與拒絕的編碼方式為π=(πOAS,πREJ),即拒絕訂單的編碼拼接于接受訂單的編碼之后,其中πOAS,πREJ∈{1,2, ,n}。在每一次迭代中檢驗(yàn)每條解上各訂單的完工時(shí)間是否超過(guò)交貨期上限,最初I為根據(jù)針對(duì)最大完工時(shí)間和能耗成本的編碼方式產(chǎn)生的一條解,πOAS為一個(gè)空集合,πREJ為一個(gè)空拒絕訂單集,調(diào)度I中所有工件并檢驗(yàn)每一個(gè)訂單的總完工時(shí)間;若訂單i在交貨期上限之前完工,則從I刪除訂單i,并將訂單i的編碼從πOAS集合中的第一個(gè)有效位置開(kāi)始添加;否則將訂單i添加到πREJ;重復(fù)該過(guò)程,直到檢測(cè)到I編碼集合變?yōu)榭占稀?/p>

    例如6訂單2工序的情況,且每個(gè)訂單具有相同交貨期,判斷每一個(gè)訂單的最大完工時(shí)間是否超過(guò)交貨期上限后,拒絕最后2個(gè)訂單,則該解表示為π={6,5,2,4,3,1},其中接受訂單集合為{6,5,2,4},拒絕訂單集合為{3,1}。

    2.5 初始解的產(chǎn)生

    對(duì)于置換流水線調(diào)度問(wèn)題,NEH 算法是目前最有效的啟發(fā)式算法之一[13],它最早由文獻(xiàn)[14]提出。本文使用NEH算法產(chǎn)生初始解。初始訂單集的訂單按目標(biāo)值的降序排序,調(diào)度排序前兩個(gè)的訂單以獲得部分解決方案,依次將其余訂單插入當(dāng)前排序中所有空位,保留目標(biāo)值最優(yōu)的當(dāng)前排序,直至未排序的訂單數(shù)為零為止。NEH算法流程如圖4所示。

    圖4 NEH算法流程

    2.6 鄰域移動(dòng)方式

    假設(shè)π為當(dāng)前解,通過(guò)兩點(diǎn)交換或插入方式獲得一定大小的當(dāng)前解鄰域N(π)。本文中選擇兩點(diǎn)交換和插入操作的概率相等。

    兩點(diǎn)交換:隨機(jī)生成兩個(gè)位置,直接交換π中兩個(gè)位置上的個(gè)體以獲得N(π);

    前插/后插:將隨機(jī)位置上的個(gè)體插入到另一位置上的個(gè)體之前或之后。

    2.7 禁忌表

    禁忌表作為禁忌搜索的存儲(chǔ)結(jié)構(gòu),可防止搜索過(guò)程中出現(xiàn)循環(huán)和避免陷入局部最優(yōu)。禁忌表的主要指標(biāo)為禁忌對(duì)象和禁忌長(zhǎng)度。為了避免算法迭代中出現(xiàn)重復(fù)移動(dòng),將每次迭代的局部最優(yōu)解的鄰域操作放入禁忌表中,這些元素在下次搜索時(shí)將不會(huì)被考慮;禁忌表長(zhǎng)度為禁忌表中禁忌對(duì)象的最大值,其大小會(huì)影響算法效果,通過(guò)Taguchi方法調(diào)整禁忌表長(zhǎng)度。

    2.8 藐視準(zhǔn)則與終止準(zhǔn)則

    在本文中,當(dāng)禁忌搜索算法中出現(xiàn)候選解全部被禁忌,或者某個(gè)禁忌候選解的適應(yīng)度值優(yōu)于當(dāng)前最好解的適應(yīng)度值,則解禁某個(gè)禁忌候選解作為新的當(dāng)前最好解。

    本文將最大迭代次數(shù)設(shè)置為終止準(zhǔn)則。

    3 實(shí)驗(yàn)結(jié)果與分析

    所有實(shí)驗(yàn)在環(huán)境Intel? CoreTMi5-4300U,主頻為2.2 GHz,內(nèi)存為16 GB的個(gè)人電腦上運(yùn)行,編程環(huán)境為Matlab R2017a。實(shí)驗(yàn)設(shè)定訂單數(shù)量n為10、20、30、40或50,工序數(shù)m為6、8、10,組成15種例如10/6,20/8,30/10等形如n/m的訂單規(guī)模,15種規(guī)模對(duì)應(yīng)150個(gè)隨機(jī)算例,每個(gè)算例運(yùn)行30次。根據(jù)文獻(xiàn)[15],訂單加工時(shí)間pi,k在[0.3, 2,5]中隨機(jī)生成的,機(jī)器的不同運(yùn)行方式的能耗在[2,20]內(nèi)隨機(jī)產(chǎn)生。工廠為24小時(shí)輪班式工作制,分時(shí)電價(jià)如表2所示。每種規(guī)模里所有訂單的相同固定收入Ri在[1 000,2 000]內(nèi)隨機(jī)產(chǎn)生。本文用于測(cè)試同一算例的所有算法的終止條件為最大迭代次數(shù)Max_Iter=200。

    表1 電鈴故障模式

    交貨期根據(jù)所有訂單的總加工時(shí)間平均值產(chǎn)生,如式(20)所示:

    (20)

    3.1 算法性能對(duì)比

    大量研究表明,遺傳算法(GA,genetic algorithm)、禁忌搜索算法(TS,tabu search)、迭代局部搜索算法(ILS,iterated local search)等用于求解傳統(tǒng)PFSP是有效的。通過(guò)對(duì)比本文所提出的EHTS算法、僅加入OAR編碼方式的改進(jìn)禁忌搜索(ITS-OAR,improved tabu search-OAR)與傳統(tǒng)TS算法在優(yōu)化PFSOAS模型下的算法性能,分析EHTS算法的優(yōu)勢(shì);并且通過(guò)EHTS優(yōu)化隨機(jī)實(shí)例,分析節(jié)能調(diào)整與交貨期配置策略的有效性。

    由于元啟發(fā)式算法的性能對(duì)參數(shù)敏感,因此使用Taguchi正交實(shí)驗(yàn)設(shè)計(jì)法對(duì)EHTS、ITS-OAR、TS進(jìn)行算法調(diào)參。調(diào)整后參數(shù)如表3所示,其中n為訂單總數(shù)量,max{}為取最大值函數(shù),round{}為取整函數(shù),sqrt{}為開(kāi)平方函數(shù)。

    表3 算法調(diào)整后的參數(shù)

    表4給出了15種規(guī)模實(shí)例問(wèn)題在傳統(tǒng)TS、ITS-OAR、EHTS的優(yōu)化性能對(duì)比,指標(biāo)分別為總凈利潤(rùn)平均值、拒絕訂單數(shù)、能耗成本平均值以及相應(yīng)的總凈利潤(rùn)增長(zhǎng)率、能耗成本降低率。TNRbest為EHTS優(yōu)化后的訂單總凈利潤(rùn)值,TECbest為加工EHTS優(yōu)化解決方案所耗費(fèi)的電費(fèi)成本??們衾麧?rùn)增長(zhǎng)率(IRTNR)與能耗成本降低率(DRTEC)分別如式(21)和式(22)。其中,TNR與TEC分別為TS、ITS-OAR優(yōu)化后的總凈利潤(rùn)值與能耗成本值??紤]到算例總數(shù)和規(guī)模總數(shù)較大,對(duì)150個(gè)算例運(yùn)行30次后的優(yōu)化結(jié)果取平均值進(jìn)行統(tǒng)計(jì)實(shí)驗(yàn)。根據(jù)表4所繪制的不同工序下三種算法性能對(duì)比如圖5所示。

    (21)

    (22)

    從圖5可以看出,在優(yōu)化PFSOAS問(wèn)題方面,EHTS算法的優(yōu)化性能更優(yōu)于傳統(tǒng)TS算法和ITS-OAR算法。根據(jù)表4,從訂單總利潤(rùn)的角度來(lái)看,EHTS算法對(duì)比傳統(tǒng)TS算法,其總凈利潤(rùn)值可以提高5%~11%左右,平均提高大約8%;EHTS算法比ITS-OAR算法優(yōu)化得到的總凈利潤(rùn)值高1%~7%左右,平均提高大約4%;說(shuō)明通過(guò)能耗調(diào)整和交貨期配置策略,可以使得企業(yè)訂單總利潤(rùn)值顯著提高。由于該目標(biāo)中帶有拖期懲罰項(xiàng),會(huì)與總凈利潤(rùn)增長(zhǎng)相互制約,隨著訂單規(guī)模的增長(zhǎng),某些訂單的生產(chǎn)仍處于高峰時(shí)段,總凈利潤(rùn)增長(zhǎng)幅度一直保持在5%左右。從能源消耗成本的角度來(lái)看,EHTS算法比前兩種算法具有顯著優(yōu)勢(shì),能耗成本減少大約9%。拒絕訂單數(shù)也相應(yīng)減少。因此EHTS算法有著較好的有效性和穩(wěn)定性。

    圖5 3種算法的性能比較

    表4 不同規(guī)模下不同算法的優(yōu)化結(jié)果

    4 結(jié)束語(yǔ)

    本文在傳統(tǒng)能效優(yōu)化的置換流水車間調(diào)度能耗模型角度上加入了訂單接受與調(diào)度,針對(duì)問(wèn)題的特殊性,提出了一種新型的節(jié)能禁忌搜索算法,提高傳統(tǒng)算法對(duì)企業(yè)總凈利潤(rùn)和能耗成本的優(yōu)化效果,訂單總凈利潤(rùn)平均提高了13%,拒絕訂單數(shù)與能耗成本也相應(yīng)減少。通過(guò)NEH構(gòu)造啟發(fā)式產(chǎn)生初始解、節(jié)能調(diào)整策略的改進(jìn),可提高算法搜索效率,改善初始解質(zhì)量,通過(guò)多種算法的性能對(duì)比,較單一算法相比,該EHTS算法具有較強(qiáng)全局搜索能力的優(yōu)勢(shì),具有一定的實(shí)用價(jià)值。本文的研究?jī)?nèi)容和結(jié)論將推動(dòng)綠色制造領(lǐng)域通過(guò)科學(xué)可行的調(diào)度方法,實(shí)現(xiàn)節(jié)能降耗、高效益的生產(chǎn)目標(biāo),具有一定的實(shí)踐意義。

    猜你喜歡
    交貨期搜索算法凈利潤(rùn)
    改進(jìn)的和聲搜索算法求解凸二次規(guī)劃及線性規(guī)劃
    帶有安裝時(shí)間與維修活動(dòng)的單機(jī)排序問(wèn)題
    三季報(bào)凈利潤(rùn)大幅增長(zhǎng)股(滬市主板)
    成本結(jié)構(gòu)離散的兩屬性電子逆向拍賣機(jī)制設(shè)計(jì)
    2014年凈利潤(rùn)增(降)幅最大的50家公司
    基于汽車接力的潮流轉(zhuǎn)移快速搜索算法
    基于逐維改進(jìn)的自適應(yīng)步長(zhǎng)布谷鳥(niǎo)搜索算法
    前三季度凈利潤(rùn)同比增(減)幅前50家公司
    復(fù)雜環(huán)境下上海WT企業(yè)交貨期優(yōu)化研究
    帶有退化效應(yīng)的多個(gè)交貨期窗口單機(jī)排序問(wèn)題
    国产一卡二卡三卡精品| 免费观看a级毛片全部| 两人在一起打扑克的视频| 亚洲欧洲日产国产| av国产精品久久久久影院| 国语对白做爰xxxⅹ性视频网站| 侵犯人妻中文字幕一二三四区| 欧美 日韩 精品 国产| 波多野结衣av一区二区av| 日韩欧美一区视频在线观看| 人体艺术视频欧美日本| 一区二区三区精品91| 国产欧美日韩综合在线一区二区| 国产精品一区二区在线不卡| 欧美精品啪啪一区二区三区 | 欧美另类一区| 人人澡人人妻人| 日韩免费高清中文字幕av| 亚洲av综合色区一区| av在线老鸭窝| 人人妻人人澡人人爽人人夜夜| 亚洲精品自拍成人| 欧美人与性动交α欧美精品济南到| 精品久久蜜臀av无| 一级毛片我不卡| 一区二区三区激情视频| kizo精华| 亚洲一码二码三码区别大吗| 极品人妻少妇av视频| 久久午夜综合久久蜜桃| 一边摸一边做爽爽视频免费| 在线观看免费视频网站a站| 90打野战视频偷拍视频| 国产熟女午夜一区二区三区| 久久人人爽av亚洲精品天堂| 色婷婷久久久亚洲欧美| 一个人免费看片子| 色94色欧美一区二区| 嫁个100分男人电影在线观看 | 97在线人人人人妻| 日本黄色日本黄色录像| 日韩一本色道免费dvd| e午夜精品久久久久久久| 久久热在线av| 激情视频va一区二区三区| 男女无遮挡免费网站观看| 午夜福利在线免费观看网站| 色婷婷av一区二区三区视频| 国产精品三级大全| 成人国产av品久久久| 丝袜脚勾引网站| 不卡av一区二区三区| 久久久久久久精品精品| www.熟女人妻精品国产| 精品久久蜜臀av无| 亚洲国产av新网站| 欧美成狂野欧美在线观看| 中文欧美无线码| 免费在线观看黄色视频的| 欧美精品一区二区免费开放| 亚洲av日韩精品久久久久久密 | 午夜免费成人在线视频| 丝袜在线中文字幕| 精品人妻1区二区| 国产野战对白在线观看| 亚洲欧美激情在线| 在线观看免费午夜福利视频| 欧美在线黄色| 国产av一区二区精品久久| 久久天堂一区二区三区四区| 91字幕亚洲| 久久精品国产亚洲av涩爱| 日韩 欧美 亚洲 中文字幕| 日本欧美视频一区| 国产深夜福利视频在线观看| 精品人妻1区二区| 精品视频人人做人人爽| 蜜桃在线观看..| www.999成人在线观看| 一本一本久久a久久精品综合妖精| 久久久久精品人妻al黑| 成人亚洲欧美一区二区av| 满18在线观看网站| 国产精品久久久久久人妻精品电影 | 欧美人与善性xxx| 精品国产一区二区三区久久久樱花| 成人亚洲精品一区在线观看| 亚洲欧美日韩高清在线视频 | 国产精品一区二区精品视频观看| 中文字幕精品免费在线观看视频| 一本色道久久久久久精品综合| 两人在一起打扑克的视频| 精品一区二区三区四区五区乱码 | 日本猛色少妇xxxxx猛交久久| 欧美少妇被猛烈插入视频| 亚洲av电影在线观看一区二区三区| 久久热在线av| 国产av国产精品国产| 精品少妇内射三级| 天天躁日日躁夜夜躁夜夜| 亚洲午夜精品一区,二区,三区| 亚洲精品成人av观看孕妇| 老司机靠b影院| 丁香六月天网| 久久av网站| 大陆偷拍与自拍| 成人午夜精彩视频在线观看| 男女午夜视频在线观看| 韩国高清视频一区二区三区| 亚洲情色 制服丝袜| 18禁黄网站禁片午夜丰满| 日韩一本色道免费dvd| 欧美日韩成人在线一区二区| 后天国语完整版免费观看| 视频区图区小说| 国产97色在线日韩免费| 国产精品免费视频内射| 亚洲国产最新在线播放| 亚洲一区中文字幕在线| 欧美日韩黄片免| 搡老岳熟女国产| 久久精品国产亚洲av高清一级| xxx大片免费视频| 色婷婷久久久亚洲欧美| 男女床上黄色一级片免费看| 69精品国产乱码久久久| 一区二区三区乱码不卡18| 国产亚洲av片在线观看秒播厂| 国产熟女欧美一区二区| 国产成人欧美在线观看 | 国产有黄有色有爽视频| 国产片内射在线| 韩国高清视频一区二区三区| 免费高清在线观看日韩| 2018国产大陆天天弄谢| 啦啦啦啦在线视频资源| 成人免费观看视频高清| 各种免费的搞黄视频| 欧美 亚洲 国产 日韩一| 日本欧美视频一区| 国产精品一国产av| 国语对白做爰xxxⅹ性视频网站| 精品亚洲乱码少妇综合久久| 久久99一区二区三区| 久久久精品94久久精品| 嫩草影视91久久| 午夜福利在线免费观看网站| 欧美黄色淫秽网站| 久久久久网色| 久久精品国产亚洲av涩爱| 女警被强在线播放| 曰老女人黄片| 人人妻人人爽人人添夜夜欢视频| 欧美xxⅹ黑人| 亚洲图色成人| 久久热在线av| 国产一区二区 视频在线| 日日摸夜夜添夜夜爱| 国产精品一区二区精品视频观看| 91成人精品电影| 搡老岳熟女国产| 亚洲av在线观看美女高潮| 国产爽快片一区二区三区| 久久久久久免费高清国产稀缺| 国产精品免费视频内射| 老司机深夜福利视频在线观看 | 国产精品 国内视频| 人体艺术视频欧美日本| 18禁黄网站禁片午夜丰满| 精品一区在线观看国产| 9191精品国产免费久久| 国产爽快片一区二区三区| 午夜福利一区二区在线看| av有码第一页| 国产精品一国产av| 久久综合国产亚洲精品| 国产日韩欧美视频二区| 欧美 日韩 精品 国产| 中文字幕精品免费在线观看视频| 韩国高清视频一区二区三区| 中文字幕最新亚洲高清| 日本色播在线视频| 91精品三级在线观看| 国产91精品成人一区二区三区 | 国产成人91sexporn| 亚洲色图综合在线观看| 美国免费a级毛片| 国产精品99久久99久久久不卡| 不卡av一区二区三区| 在线观看www视频免费| 午夜福利乱码中文字幕| 美女脱内裤让男人舔精品视频| 亚洲欧洲国产日韩| 校园人妻丝袜中文字幕| 亚洲熟女毛片儿| 午夜福利一区二区在线看| 国产成人精品无人区| 狂野欧美激情性xxxx| 欧美日韩国产mv在线观看视频| 午夜两性在线视频| 考比视频在线观看| 纵有疾风起免费观看全集完整版| 色综合欧美亚洲国产小说| 超碰成人久久| 国产精品久久久av美女十八| 国产高清不卡午夜福利| 欧美成狂野欧美在线观看| 大码成人一级视频| 日韩制服骚丝袜av| 国产成人啪精品午夜网站| 欧美另类一区| 人妻人人澡人人爽人人| 亚洲国产欧美日韩在线播放| av欧美777| 免费女性裸体啪啪无遮挡网站| 国产成人一区二区三区免费视频网站 | 最新在线观看一区二区三区 | 男人舔女人的私密视频| 一本一本久久a久久精品综合妖精| 久久免费观看电影| 高清av免费在线| 婷婷色综合www| av在线app专区| 老司机影院毛片| 七月丁香在线播放| 欧美日韩黄片免| 中文字幕人妻丝袜制服| 亚洲精品av麻豆狂野| 午夜日韩欧美国产| 精品国产一区二区三区四区第35| 午夜av观看不卡| 一区二区三区乱码不卡18| 少妇被粗大的猛进出69影院| 在线av久久热| 久久久久久久国产电影| 亚洲av日韩在线播放| xxxhd国产人妻xxx| 另类亚洲欧美激情| 国产成人av教育| 男人爽女人下面视频在线观看| a级片在线免费高清观看视频| 黄色片一级片一级黄色片| 下体分泌物呈黄色| 久久国产精品大桥未久av| 欧美中文综合在线视频| 伊人久久大香线蕉亚洲五| 久久人妻熟女aⅴ| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产毛片av蜜桃av| 日本午夜av视频| 国产在视频线精品| 久久天堂一区二区三区四区| 欧美成狂野欧美在线观看| av天堂久久9| 久热爱精品视频在线9| 成年动漫av网址| 国产亚洲精品久久久久5区| 亚洲伊人久久精品综合| 最近手机中文字幕大全| 男女午夜视频在线观看| 久久女婷五月综合色啪小说| 日韩一卡2卡3卡4卡2021年| 热re99久久国产66热| 亚洲 国产 在线| 19禁男女啪啪无遮挡网站| 午夜av观看不卡| 50天的宝宝边吃奶边哭怎么回事| www.自偷自拍.com| 一级毛片黄色毛片免费观看视频| 婷婷色综合大香蕉| 久久人人爽av亚洲精品天堂| 亚洲av成人精品一二三区| 久久这里只有精品19| 国产av国产精品国产| 免费看av在线观看网站| 亚洲伊人色综图| 午夜福利免费观看在线| 一级毛片 在线播放| 欧美人与性动交α欧美软件| av在线老鸭窝| 天天操日日干夜夜撸| www.精华液| 韩国高清视频一区二区三区| 久久鲁丝午夜福利片| 18禁黄网站禁片午夜丰满| 欧美日韩国产mv在线观看视频| 涩涩av久久男人的天堂| 欧美日韩亚洲综合一区二区三区_| 最新在线观看一区二区三区 | 在线观看www视频免费| 99国产综合亚洲精品| www.精华液| 亚洲三区欧美一区| 一边摸一边做爽爽视频免费| 精品久久久精品久久久| 美女视频免费永久观看网站| 亚洲国产看品久久| 亚洲欧美一区二区三区黑人| 成人国产av品久久久| 激情视频va一区二区三区| 国产不卡av网站在线观看| 成人手机av| 丁香六月欧美| 亚洲精品一区蜜桃| 你懂的网址亚洲精品在线观看| 中文精品一卡2卡3卡4更新| 欧美日韩av久久| 国产精品国产av在线观看| 久久性视频一级片| 久久 成人 亚洲| 免费黄频网站在线观看国产| 天天操日日干夜夜撸| 成人18禁高潮啪啪吃奶动态图| 亚洲成人免费av在线播放| 在线观看免费午夜福利视频| 欧美 日韩 精品 国产| 丝袜人妻中文字幕| 国产老妇伦熟女老妇高清| 天天添夜夜摸| 精品高清国产在线一区| 亚洲国产日韩一区二区| 人人澡人人妻人| 国产精品久久久人人做人人爽| 十分钟在线观看高清视频www| 悠悠久久av| 麻豆av在线久日| 国产老妇伦熟女老妇高清| 美女脱内裤让男人舔精品视频| 亚洲精品av麻豆狂野| 国产伦理片在线播放av一区| 国产高清videossex| 搡老岳熟女国产| 99热网站在线观看| 久久精品国产a三级三级三级| 香蕉国产在线看| 国产精品免费视频内射| 国产精品 国内视频| 国产老妇伦熟女老妇高清| 一级毛片黄色毛片免费观看视频| 午夜福利影视在线免费观看| 久久人妻福利社区极品人妻图片 | 久久中文字幕一级| 美国免费a级毛片| 90打野战视频偷拍视频| 男男h啪啪无遮挡| 亚洲成av片中文字幕在线观看| 亚洲欧洲精品一区二区精品久久久| 国产在线观看jvid| 亚洲精品乱久久久久久| 国产在线观看jvid| 中文精品一卡2卡3卡4更新| 午夜福利视频在线观看免费| 菩萨蛮人人尽说江南好唐韦庄| 黄色 视频免费看| 久久精品亚洲av国产电影网| 亚洲熟女毛片儿| 亚洲精品日韩在线中文字幕| 亚洲欧美精品综合一区二区三区| 亚洲精品自拍成人| 99热全是精品| 久久性视频一级片| 黄网站色视频无遮挡免费观看| 久久久久久久久久久久大奶| 亚洲专区国产一区二区| 18在线观看网站| 一级a爱视频在线免费观看| 精品福利永久在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 国产精品国产三级专区第一集| 国产精品99久久99久久久不卡| cao死你这个sao货| 亚洲欧美精品综合一区二区三区| 又大又爽又粗| 99热全是精品| 成年人黄色毛片网站| 久久精品久久久久久久性| 免费高清在线观看视频在线观看| 国产日韩欧美视频二区| 两性夫妻黄色片| videos熟女内射| 国产精品秋霞免费鲁丝片| 在线观看www视频免费| 国产精品成人在线| 欧美 日韩 精品 国产| 两性夫妻黄色片| 久久久久久免费高清国产稀缺| 久久精品亚洲av国产电影网| 悠悠久久av| 欧美日韩一级在线毛片| 亚洲av在线观看美女高潮| 亚洲国产av新网站| 免费黄频网站在线观看国产| 午夜激情久久久久久久| 国产淫语在线视频| 亚洲欧洲日产国产| 宅男免费午夜| 一级毛片女人18水好多 | 天天躁日日躁夜夜躁夜夜| 成在线人永久免费视频| av电影中文网址| 视频区欧美日本亚洲| 国产免费视频播放在线视频| 欧美日韩综合久久久久久| 欧美乱码精品一区二区三区| 亚洲七黄色美女视频| 亚洲自偷自拍图片 自拍| 亚洲欧美一区二区三区黑人| 国产精品99久久99久久久不卡| 久久免费观看电影| 亚洲熟女毛片儿| 久久九九热精品免费| 美女视频免费永久观看网站| 亚洲欧洲国产日韩| 亚洲少妇的诱惑av| 777米奇影视久久| 精品视频人人做人人爽| 搡老乐熟女国产| 91麻豆精品激情在线观看国产 | 国产精品一区二区精品视频观看| 超色免费av| 另类精品久久| 精品国产乱码久久久久久小说| 久久久久久久大尺度免费视频| 国产国语露脸激情在线看| 亚洲精品一区蜜桃| 精品欧美一区二区三区在线| 国产精品久久久人人做人人爽| 久久午夜综合久久蜜桃| 欧美日本中文国产一区发布| 国产成人a∨麻豆精品| 日韩熟女老妇一区二区性免费视频| 王馨瑶露胸无遮挡在线观看| 久久精品国产亚洲av高清一级| 下体分泌物呈黄色| 亚洲欧美一区二区三区黑人| 亚洲成人手机| www.999成人在线观看| 2018国产大陆天天弄谢| 国产极品粉嫩免费观看在线| 亚洲精品国产av成人精品| 十八禁网站网址无遮挡| 亚洲av男天堂| 国产精品国产三级国产专区5o| 18禁观看日本| videos熟女内射| 免费高清在线观看视频在线观看| 国产成人a∨麻豆精品| 国产在线观看jvid| 欧美激情高清一区二区三区| 欧美国产精品一级二级三级| 久久久久久人人人人人| 亚洲欧美一区二区三区黑人| 国产福利在线免费观看视频| 精品第一国产精品| 一区二区日韩欧美中文字幕| 国产97色在线日韩免费| 欧美性长视频在线观看| 在线 av 中文字幕| 亚洲视频免费观看视频| bbb黄色大片| 亚洲欧美日韩另类电影网站| 在线看a的网站| 美国免费a级毛片| 美女高潮到喷水免费观看| 中文字幕最新亚洲高清| 国产精品香港三级国产av潘金莲 | 欧美乱码精品一区二区三区| 国产精品九九99| 人人妻人人添人人爽欧美一区卜| 亚洲精品久久成人aⅴ小说| 国产97色在线日韩免费| 亚洲图色成人| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品免费视频内射| 高清黄色对白视频在线免费看| 欧美日韩国产mv在线观看视频| 久久精品国产亚洲av涩爱| 激情五月婷婷亚洲| 日韩电影二区| 色视频在线一区二区三区| 波野结衣二区三区在线| 手机成人av网站| 亚洲人成电影观看| av线在线观看网站| av网站在线播放免费| 国产伦理片在线播放av一区| 亚洲国产欧美网| 大型av网站在线播放| 久久人人爽人人片av| 操出白浆在线播放| tube8黄色片| 欧美另类一区| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲 | 黑丝袜美女国产一区| 欧美 亚洲 国产 日韩一| 亚洲av国产av综合av卡| 中文字幕亚洲精品专区| 亚洲成国产人片在线观看| 午夜免费观看性视频| 精品少妇久久久久久888优播| 久久女婷五月综合色啪小说| 午夜激情av网站| 中文欧美无线码| 97在线人人人人妻| 久久久久精品国产欧美久久久 | 国产伦人伦偷精品视频| 国产日韩欧美视频二区| 菩萨蛮人人尽说江南好唐韦庄| 精品一品国产午夜福利视频| av国产精品久久久久影院| 肉色欧美久久久久久久蜜桃| 亚洲成人免费电影在线观看 | 国产真人三级小视频在线观看| 亚洲精品国产色婷婷电影| 国产亚洲欧美在线一区二区| 国产片特级美女逼逼视频| 在线观看免费日韩欧美大片| av线在线观看网站| 丰满人妻熟妇乱又伦精品不卡| 欧美激情极品国产一区二区三区| 不卡av一区二区三区| 国产黄频视频在线观看| 亚洲精品一区蜜桃| 亚洲精品国产色婷婷电影| 亚洲av日韩精品久久久久久密 | 成在线人永久免费视频| 亚洲欧美日韩另类电影网站| 国产麻豆69| 两个人免费观看高清视频| 免费一级毛片在线播放高清视频 | 岛国毛片在线播放| 亚洲午夜精品一区,二区,三区| 日本色播在线视频| 日韩视频在线欧美| 欧美日韩视频高清一区二区三区二| 天天添夜夜摸| 一边摸一边做爽爽视频免费| 操美女的视频在线观看| 熟女av电影| 女人爽到高潮嗷嗷叫在线视频| 天天躁夜夜躁狠狠久久av| 欧美日韩亚洲高清精品| xxx大片免费视频| 99国产精品免费福利视频| 欧美 亚洲 国产 日韩一| 两性夫妻黄色片| 国产精品久久久av美女十八| 午夜免费鲁丝| 男女免费视频国产| 大陆偷拍与自拍| 一级a爱视频在线免费观看| 人人妻人人澡人人看| 国产淫语在线视频| 国产1区2区3区精品| 久久99一区二区三区| 国产高清视频在线播放一区 | 日本欧美国产在线视频| 成人黄色视频免费在线看| 少妇精品久久久久久久| 一级毛片女人18水好多 | 超碰97精品在线观看| 少妇的丰满在线观看| 99国产精品免费福利视频| 欧美精品啪啪一区二区三区 | 一本色道久久久久久精品综合| 午夜福利一区二区在线看| 国产免费一区二区三区四区乱码| 又大又爽又粗| 久久久国产欧美日韩av| 高潮久久久久久久久久久不卡| 国产亚洲精品久久久久5区| 国产无遮挡羞羞视频在线观看| 一本—道久久a久久精品蜜桃钙片| 高清视频免费观看一区二区| 中文字幕精品免费在线观看视频| 成人黄色视频免费在线看| 久久狼人影院| 丁香六月欧美| 亚洲精品日韩在线中文字幕| 亚洲av在线观看美女高潮| 久久午夜综合久久蜜桃| 免费女性裸体啪啪无遮挡网站| 丝袜在线中文字幕| 亚洲欧洲日产国产| 欧美日韩综合久久久久久| 妹子高潮喷水视频| 国产精品人妻久久久影院| 欧美激情极品国产一区二区三区| 一二三四在线观看免费中文在| 精品国产一区二区久久| 亚洲人成网站在线观看播放| 一本综合久久免费| 欧美人与性动交α欧美精品济南到| 午夜福利乱码中文字幕| 午夜91福利影院| 制服诱惑二区| 精品国产一区二区三区四区第35| 99国产精品99久久久久| 免费高清在线观看视频在线观看| 国产在线观看jvid| 国产免费视频播放在线视频| 国产精品国产三级专区第一集| 亚洲欧美色中文字幕在线| 啦啦啦视频在线资源免费观看| 亚洲一区二区三区欧美精品| tube8黄色片| 亚洲精品av麻豆狂野| 免费看十八禁软件| 国产精品久久久av美女十八| 啦啦啦视频在线资源免费观看| 亚洲精品日韩在线中文字幕| 欧美精品高潮呻吟av久久| 免费一级毛片在线播放高清视频 | av有码第一页|