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

    考慮運(yùn)輸時(shí)間的分布式柔性作業(yè)車(chē)間綠色調(diào)度

    2022-11-21 10:50:48張洪亮徐公杰潘瑞林
    中國(guó)機(jī)械工程 2022年21期
    關(guān)鍵詞:工序工件機(jī)器

    張洪亮 徐公杰 鮑 薔 潘瑞林

    安徽工業(yè)大學(xué)管理科學(xué)與工程學(xué)院,馬鞍山,243032

    0 引言

    在中國(guó),制造企業(yè)消耗了全國(guó)50%以上的電能,并產(chǎn)生了至少26%的二氧化碳排放量[1]。通過(guò)研發(fā)節(jié)能設(shè)備或采用新的加工技術(shù)減小節(jié)能減排壓力的方法通常需要大量的投入,而綠色調(diào)度能在不增加企業(yè)成本的情況下有效減少碳排放并提高能源效率[2]。

    隨著信息技術(shù)以及全球化的飛速發(fā)展,許多制造企業(yè)從傳統(tǒng)的單工廠模式逐漸向能降低人工成本、提高生產(chǎn)效率[3]的分布式多工廠模式轉(zhuǎn)變。分布式柔性作業(yè)車(chē)間調(diào)度問(wèn)題(distributed flexible job shop scheduling problem,DFJSP)作為分布式車(chē)間調(diào)度問(wèn)題的重要形式之一,重點(diǎn)關(guān)注多個(gè)柔性作業(yè)車(chē)間類(lèi)型的工廠協(xié)同生產(chǎn)。現(xiàn)在對(duì)該問(wèn)題的研究主要分為單目標(biāo)DFJSP和多目標(biāo)DFJSP。單目標(biāo)DFJSP主要以完工時(shí)間為優(yōu)化目標(biāo),對(duì)應(yīng)研究中,GIOVANNI等[4]設(shè)計(jì)了一種包含局部搜索策略的遺傳算法;LIU等[5]將概率融入實(shí)參數(shù)編碼方法,設(shè)計(jì)了一種改進(jìn)的遺傳算法;吳銳等[6]提出一種改進(jìn)的人工蜂群算法;MENG等[7]提出了4種混合整數(shù)線性規(guī)劃模型和1種約束規(guī)劃模型。對(duì)于多目標(biāo)DFJSP,吳秀麗等[8]設(shè)計(jì)了一種改進(jìn)的差分進(jìn)化算法來(lái)優(yōu)化總成本和提前/延期懲罰;LI等[9]設(shè)計(jì)了一種基于Pareto的混合禁忌搜索算法來(lái)優(yōu)化完工時(shí)間、機(jī)器負(fù)荷和提前/延期懲罰。

    上述研究存在以下不足:①要求同一工件的所有工序必須在同一工廠內(nèi)完成,沒(méi)有考慮工件在工廠之間的轉(zhuǎn)移。實(shí)際上,分布式生產(chǎn)系統(tǒng)中,一些工件可能在不同的工廠加工。②要求工件在前道工序加工完成后立即執(zhí)行下一道工序,忽略了工序之間的運(yùn)輸時(shí)間。實(shí)際上,工件在工序之間流轉(zhuǎn)需要運(yùn)輸時(shí)間,且加工時(shí)間和運(yùn)輸時(shí)間之間有較強(qiáng)的耦合關(guān)系。③主要關(guān)注時(shí)間相關(guān)的目標(biāo)且以單目標(biāo)為主,沒(méi)有考慮能耗相關(guān)的目標(biāo)。隨著綠色制造的推進(jìn),在生產(chǎn)調(diào)度決策中考慮與能耗相關(guān)的目標(biāo)變得至關(guān)重要。

    綜上所述,考慮運(yùn)輸約束的DFJSP研究具有重要意義,盡管一些學(xué)者在單個(gè)工廠的車(chē)間調(diào)度問(wèn)題中考慮了運(yùn)輸時(shí)間[1,10-13],但考慮運(yùn)輸時(shí)間的DFJSP研究還很少。此外,尚未有學(xué)者在DFJSP的研究中考慮能耗相關(guān)目標(biāo)。因此,本文以考慮運(yùn)輸時(shí)間的分布式柔性作業(yè)車(chē)間綠色調(diào)度問(wèn)題(distributed flexible job shop green scheduling problem with transportation time,DFJGSPT)為研究對(duì)象,建立最小化完工時(shí)間和總能耗的混合整數(shù)規(guī)劃模型,并設(shè)計(jì)一種改進(jìn)的非支配排序遺傳算法(improved non-dominated sorting genetic algorithm Ⅱ,INSGA-Ⅱ)。通過(guò)測(cè)試45個(gè)算例,以及與常用的多目標(biāo)進(jìn)化算法的對(duì)比來(lái)驗(yàn)證本文提出的INSGA-Ⅱ的有效性。

    1 問(wèn)題描述及數(shù)學(xué)模型

    1.1 問(wèn)題描述

    DFJGSPT可以描述為:N個(gè)工件需要在F個(gè)類(lèi)型為柔性作業(yè)車(chē)間的分布式工廠內(nèi)進(jìn)行加工。每個(gè)工件有ni道工序,每個(gè)工廠有Mf臺(tái)機(jī)器。工件需要通過(guò)運(yùn)輸工具在工廠以及機(jī)器之間運(yùn)輸。工件的加工時(shí)間,工件在機(jī)器、工廠間的運(yùn)輸時(shí)間,以及機(jī)器的相關(guān)能耗信息已知。調(diào)度目標(biāo)為最小化最大完工時(shí)間和總能耗。假設(shè):①所有工件和機(jī)器0時(shí)刻可用;②每臺(tái)機(jī)器同一時(shí)刻只能加工一道工序;③每個(gè)工件同一時(shí)刻只能在一臺(tái)機(jī)器上加工;④不考慮中斷情況;⑤同一工件不同工序之間有順序約束;⑥有足夠的運(yùn)輸工具完成工件的轉(zhuǎn)移;⑦不考慮裝卸載時(shí)間。

    1.2 數(shù)學(xué)模型

    基于以上描述可得DFJGSPT的數(shù)學(xué)模型:

    minCmax=min(max(Ci))

    (1)

    minTEC=min(PE+TE+IE)

    (2)

    (3)

    (4)

    (5)

    αilk+βifu≤1

    (6)

    cijkf=sijkf+pijkfxijkf

    (7)

    cijkf≥cghkf+pijkf-M(1-yijghkf)

    (8)

    (9)

    (10)

    (11)

    式中,Cmax為最大完工時(shí)間;Ci為工件i的完工時(shí)間;TEC為總能耗;PE為所有機(jī)器的加工能耗;TE為所有運(yùn)輸任務(wù)的運(yùn)輸能耗;IE為所有機(jī)器的空閑能耗;xijkf為0-1決策變量,如果Oij在工廠f的機(jī)器k上加工,則xijkf為1,否則為0;F為工廠數(shù)量;Mf為工廠f中機(jī)器的數(shù)量;yijghkf為0-1決策變量,工廠f的機(jī)器k加工完Ogh后加工Oij,則yijghkf為1,否則為0;n為工件數(shù)量;ni為工件i的工序數(shù)量;αilk為0-1決策變量,如果工件i從機(jī)器l運(yùn)送到機(jī)器k,則αilk為1,否則為0;βifu為0-1決策變量,如果工件i從工廠f運(yùn)送到工廠u,則βifu為1,否則為0;cijkf為Oij在工廠f的機(jī)器k上的加工結(jié)束時(shí)間;sijkf為Oij在工廠f的機(jī)器k上的加工開(kāi)始時(shí)間;pijkf為Oij在工廠f的機(jī)器k上的加工時(shí)長(zhǎng);M為一個(gè)足夠大的正數(shù);pekf為工廠f中機(jī)器k的單位加工能耗;iekf為工廠f中機(jī)器k的單位空閑能耗;fte為工廠之間的單位運(yùn)輸能耗;TTfu為工件從工廠f到工廠u的運(yùn)輸時(shí)間;mte為機(jī)器之間的單位運(yùn)輸能耗;ttlk為工件從機(jī)器l到機(jī)器k的運(yùn)輸時(shí)間;i、g為工件索引;j、h為工序索引;l、k為機(jī)器索引;f、u為工廠索引。

    式(1)、式(2)為優(yōu)化目標(biāo)——最小化最大完工時(shí)間和最小化總能耗;式(3)保證每道工序只能在一個(gè)工廠的一臺(tái)機(jī)器上加工;式(4)、式(5)表示每道工序的緊后或者緊前工序最多只有一道;式(6)確保每個(gè)工件不能同時(shí)在機(jī)器和車(chē)間之間進(jìn)行運(yùn)輸;式(7)表示工序的完工時(shí)間等于開(kāi)始時(shí)間加上加工時(shí)長(zhǎng);式(8)確保機(jī)器不能同時(shí)加工多個(gè)工件;式(9)~式(11)分別表示機(jī)器加工能耗、機(jī)器空閑能耗和運(yùn)輸能耗。

    2 改進(jìn)NSGA-Ⅱ

    NSGA-Ⅱ是一種有效求解多目標(biāo)優(yōu)化問(wèn)題的進(jìn)化算法,但在求解多目標(biāo)柔性作業(yè)車(chē)間調(diào)度問(wèn)題時(shí),存在早熟和易陷入局部最優(yōu)的不足。本文針對(duì)DFJGSPT的特性,對(duì)NSGA-Ⅱ進(jìn)行以下改進(jìn):①設(shè)計(jì)了同時(shí)考慮加工時(shí)間和能耗的初始化方法,以提高初始種群的質(zhì)量;②設(shè)計(jì)了考慮運(yùn)輸時(shí)間的貪婪插入解碼方法,將染色體轉(zhuǎn)換為可行有效的調(diào)度方案;③采用多父代交叉和兩點(diǎn)插入及隨機(jī)選擇的變異方式進(jìn)行種群更新;④設(shè)計(jì)了一種變鄰域搜索策略來(lái)提高Pareto前沿的質(zhì)量。

    2.1 基于工序和機(jī)器的雙層編碼

    本文采用基于工序和機(jī)器的雙層編碼方式對(duì)染色體編碼,工廠的分配則通過(guò)解碼策略確定。以表1所示的數(shù)據(jù)為例,編碼方案示例見(jiàn)圖1。

    表1 MODFJGSPT示例

    圖1 編碼示例

    2.2 考慮加工時(shí)間和能耗的種群初始化

    本文設(shè)計(jì)了一種考慮加工時(shí)間和能耗的種群初始化方法:工序?qū)拥木幋a通過(guò)隨機(jī)的方式產(chǎn)生;機(jī)器層編碼通過(guò)設(shè)計(jì)的全局選擇(globalselection,GS)、局部選擇(local selection,LS)、隨機(jī)選擇(random select,RS)3種策略產(chǎn)生。參考文獻(xiàn)[14],通過(guò)全局選擇、局部選擇和隨機(jī)選擇產(chǎn)生的初始解數(shù)量比例設(shè)為6∶3∶1。種群初始化的偽代碼如下。

    算法1:種群初始化 ∥ Gr、Lr和Rr分別為通過(guò)全局選擇、局部選擇和隨機(jī)選擇所產(chǎn)生的初始解數(shù)量在種群占比,其值分別為0.6、0.3和0.11:采用隨機(jī)的方式產(chǎn)生工序?qū)拥木幋aOS2:創(chuàng)建兩個(gè)值為0的數(shù)組分別記錄所有機(jī)器的加工時(shí)間PTk和加工能耗PEk(1≤k≤m)3:If n<=Nind*(1-Rr)4:隨機(jī)選擇一個(gè)工件i ∥每個(gè)工件只被選擇一次5: For j=1 to ni6: TempPTk=pijk +PTk;TempPEk=pijk×pek + PEk(k∈K_use) ∥ K_use:Oij的可用機(jī)器集7: 根據(jù)TempPTk和TempPEk計(jì)算Pareto支配關(guān)系8: 選擇具有非支配TempPTk和TempPEk的機(jī)器k的索引放入機(jī)器層編碼MA中9: 更新PTk和PEk10: End For11: If n>Nind*Gr&&n <=Nind*(Gr+Lr)12: 將數(shù)組PTk和PEk的值設(shè)為013: End If14: 轉(zhuǎn)到步驟5,直到所有作業(yè)都被選中一次15:Else16:采用隨機(jī)的方式產(chǎn)生機(jī)器層編碼MA17:End If

    2.3 考慮運(yùn)輸時(shí)間的貪婪插入解碼方法

    通過(guò)編碼建立了工件和機(jī)器的映射關(guān)系,為進(jìn)一步確定工件在對(duì)應(yīng)機(jī)器上的開(kāi)工時(shí)間,本文設(shè)計(jì)了一種考慮運(yùn)輸時(shí)間的貪婪插入解碼方法。首先通過(guò)機(jī)器編碼確定工件所在的工廠,然后考慮運(yùn)輸時(shí)間并引入貪婪的思想,將工件盡可能地插入機(jī)器的空閑時(shí)段內(nèi),以減少機(jī)器的空閑時(shí)間和能耗,提高生產(chǎn)效率。解碼方法的偽代碼如下。

    算法2:解碼方法 ∥ TO:工序總數(shù);ts/tf:第一個(gè)滿(mǎn)足插入條件的空閑時(shí)間段的開(kāi)始/結(jié)束時(shí)間1:NOMk=?(k=1,2,…,m) %機(jī)器k上已經(jīng)安排的工序數(shù)量2:For h=1 to TO3: Oij=OS(h)4: 獲取Oij的加工機(jī)器k,pijk和運(yùn)輸?shù)竭_(dá)機(jī)器k的時(shí)間ta5: If j==1 &&NOMk==?6: sijk=0;cijk=sijk + pijk7: Else If j ~=1 &&NOMk==?8: sijk=ta;cijk=sijk+pijk9: Else If j==1 &&NOMk~=?10: Ins=find(IDk>=pijk)11: If Ins==?12: sijk=max{TMk,ta};cijk=sijk+pijk13: Else14: sijk=ts;cijk=sijk+pijk15: End If16: Else If j ~=1 &&NOMk~=?17: Ins=find(IDk>=pijk)18: NIns=find(tf-max{ts,ta} >=pijk)19: If NIns==?20: sijk=max{TMk,ta};cijk=sijk+pijk21: Else22: sijk=max{ts,ta};cijk=sijk+pijk23: End If24: End If25: 更新機(jī)器k的空閑時(shí)間時(shí)段IDk、NOMk和TMk26:End For

    2.4 多父代交叉操作

    為了融合多父代的信息,獲得更高質(zhì)量的子代,并增加種群的多樣性,本文采用多父代交叉方式[15]更新種群。工序?qū)雍蜋C(jī)器層的交叉方式分別如圖2、圖3所示。

    圖2 工序?qū)咏徊?/p>

    圖3 機(jī)器層交叉

    工序?qū)咏徊娴牟襟E如下:

    (1)隨機(jī)選擇父代個(gè)體P1、P2和P3;

    (2)將工件索引集合Job={1,2,…,n}隨機(jī)劃分為2個(gè)互不包含的集合Job1和Job2;

    (3)將P1中包含集合Job2的元素復(fù)制到子代C1,將P3中包含集合Job1的元素復(fù)制到子代C1;

    (3)將P2中包含集合Job1的元素復(fù)制到子代C2,將P3中包含集合Job2的元素復(fù)制到子代C2。

    機(jī)器層交叉的步驟如下:

    (1)隨機(jī)選擇父代個(gè)體P1、P2和P3;

    (2)隨機(jī)產(chǎn)生一個(gè)長(zhǎng)度與工序總數(shù)一致且由 0、1 組成的集合R;

    (3)在P2和P3中隨機(jī)選出與R中的1位置對(duì)應(yīng)的編碼,復(fù)制到子代C1中的相應(yīng)位置上;

    (4)在P1和P3中隨機(jī)選出與R中的1位置對(duì)應(yīng)的編碼,復(fù)制到子代C2中的相應(yīng)位置;

    (5)將P1和P2中其他的編碼分別保留到子代C1和C2中。

    2.5 變異操作

    對(duì)工序?qū)雍蜋C(jī)器層編碼分別采用不同的變異方式:①工序?qū)?,采用兩點(diǎn)插入變異方式,隨機(jī)選擇2個(gè)互不相鄰的位置l1和l2(l1

    圖4 工序?qū)幼儺?/p>

    2.6 變鄰域搜索

    為進(jìn)一步提高個(gè)體的質(zhì)量,本文設(shè)計(jì)了一種用于Pareto前沿的變鄰域搜索策略。通過(guò)以下兩種鄰域結(jié)構(gòu)獲得新解:①隨機(jī)選擇工序?qū)泳幋a方案的2個(gè)位置,將這2個(gè)位置之間的編碼進(jìn)行翻轉(zhuǎn);②獲取運(yùn)輸時(shí)間最長(zhǎng)的2道相鄰工序,將后一道工序換到另一臺(tái)運(yùn)輸時(shí)間短的機(jī)器上進(jìn)行加工。如果新解更優(yōu),則更新。

    2.7 算法整體步驟

    INSGA-Ⅱ的主要步驟如下:①初始化參數(shù),產(chǎn)生初始種群;②交叉變異產(chǎn)生子代;③解碼并對(duì)種群進(jìn)行非支配排序;④對(duì)Pareto前沿進(jìn)行變鄰域搜索;⑤通過(guò)對(duì)非支配排序和擁擠距離計(jì)算,獲得下一代種群;⑥如果滿(mǎn)足終止條件則輸出Pareto前沿,否則轉(zhuǎn)到②。

    3 實(shí)驗(yàn)分析

    本文算法基于MATLAB語(yǔ)言實(shí)現(xiàn),并在Intel CPU E3-1240 V5、8 GB RAM的電腦上運(yùn)行。

    3.1 算例生成及評(píng)價(jià)指標(biāo)

    鑒于目前尚沒(méi)有DFJGSPT測(cè)試算例,本文對(duì)文獻(xiàn)[16]中的基準(zhǔn)算例進(jìn)行拓展,生成了45組測(cè)試算例。具體方式如下:對(duì)每個(gè)基準(zhǔn)算例分別考慮2、4和6個(gè)工廠的場(chǎng)景,工廠內(nèi)機(jī)器間的運(yùn)輸時(shí)間在[1,5] min內(nèi)隨機(jī)產(chǎn)生,工廠間的運(yùn)輸時(shí)間在[10,15] min內(nèi)隨機(jī)產(chǎn)生。機(jī)器單位加工能耗pek∈[10,15]kW,機(jī)器單位空閑能耗iek∈[1,3] kW,機(jī)器間的單位運(yùn)輸能耗mte=2 kW,工廠間的單位運(yùn)輸能耗fte=5 kW。采用多目標(biāo)優(yōu)化問(wèn)題中廣泛應(yīng)用的覆蓋率C[3]和反世代距離IGD[17]評(píng)估算法性能,具體含義和公式如下:

    (12)

    (13)

    式中,A、B分別為需要對(duì)比的兩種算法的Pareto前沿;|B|為集合B的大??;|PF*|為第一前沿的非支配解集PF*的大?。籨(x,y)為x點(diǎn)與y點(diǎn)之間的歐氏距離。

    C(A,B)越大,A越好;IGD(A,PF*)越小,A越好。在本文研究中,各實(shí)例的PF*由每一種算法運(yùn)行10次得到的非支配解集綜合形成。

    3.2 算法參數(shù)

    本文提出的INSGA-Ⅱ有4個(gè)參數(shù):種群規(guī)模(Nind)、交叉率(Pc)、突變率(Pm)和最大迭代次數(shù)(Maxit)。為保證參數(shù)設(shè)置的合理性,參考文獻(xiàn)[13],采用正交試驗(yàn)確定最佳的參數(shù)組合。每個(gè)參數(shù)選擇4個(gè)水平,各個(gè)水平的值如表2所示。根據(jù)參數(shù)和水平的個(gè)數(shù),選擇L16(45)正交表進(jìn)行試驗(yàn)。針對(duì)算例Mk2-01,算法在每種參數(shù)組合下運(yùn)行10次,以IGD作為各種參數(shù)組合的評(píng)價(jià)指標(biāo),計(jì)算結(jié)果如表3所示。

    表2 參數(shù)水平表

    表3 正交試驗(yàn)表

    根據(jù)表3的結(jié)果繪制出參數(shù)水平的趨勢(shì),見(jiàn)圖6??芍?dāng)Nind為第一水平、Pc為第四水

    (a)Nind的水平趨勢(shì)圖

    平、Pm為第二水平、Maxit為第三水平時(shí),IGD最小。因此,最終確定的參數(shù)組合為Nind=50,Pc=0.9,Pm=0.1,Maxit=300。

    3.3 改進(jìn)策略的有效性

    為驗(yàn)證本文設(shè)計(jì)的種群初始化策略和變鄰域搜索策略的有效性,本文將INSGA-Ⅱ(A1)與采用隨機(jī)初始化策略的INSGA-Ⅱ(A2)和不包含變鄰域搜索策略的INSGA-Ⅱ(A3)進(jìn)行對(duì)比。

    參考文獻(xiàn)[7],以Mk01、Mk04、Mk09、Mk12和Mk15在不同工廠數(shù)下對(duì)應(yīng)的15組算例為測(cè)試算例,對(duì)這3種算法進(jìn)行測(cè)試。每種算法運(yùn)行10次,采用C和IGD比較算法的性能。在15組算例測(cè)試結(jié)果中,C(A1,A2)=C(A1,A3)=1,C(A2,A1)=C(A3,A1)=0,說(shuō)明算法A1可以獲得更優(yōu)的Pareto解集。此外,IGD(A1)均為0,IGD(A2)、IGD(A2)如表4所示,可以看出IGD(A1)遠(yuǎn)小于IGD(A2)和IGD(A3),進(jìn)一步表明算法A1獲得的Pareto解集更優(yōu)。測(cè)試結(jié)果可以證明本文設(shè)計(jì)的種群初始化策略和變鄰域搜索策略的有效性。

    表4 改進(jìn)策略測(cè)試結(jié)果

    3.4 與其他算法對(duì)比分析

    為進(jìn)一步驗(yàn)證本文算法的有效性,將其與以被廣泛應(yīng)用的多目標(biāo)優(yōu)化算法NSGA-Ⅱ、SPEA2和MOEA/D進(jìn)行對(duì)比。對(duì)比算法的參數(shù)設(shè)置與INSGA-Ⅱ相同,其中,SPEA2和MOEA/D的外部存檔大小設(shè)置為種群大小。每個(gè)算法運(yùn)行10次,采用C和IGD比較算法的性能。所有算例下,C(INSGA-Ⅱ,NSGA-Ⅱ)=C(INSGA-Ⅱ,SPEA2)=C(INSGA-Ⅱ,MOEA/D)=1,C(NSGA-Ⅱ,INSGA-Ⅱ)=C(SPEA2,INSGA-Ⅱ)=C(MOEA/D,INSGA-Ⅱ)=0,這表明INSGA-Ⅱ得到的解優(yōu)于另外3種算法得到的解。IGD(INSGA-Ⅱ)均為0,另外3種算法的IGD如表5所示,可以看出,IGD(INSGA-Ⅱ)遠(yuǎn)小于IGD(NSGA-Ⅱ)、IGD(SPEA2)和IGD(MOEA/D),這與C值結(jié)果相互對(duì)應(yīng)。此外,隨著工廠數(shù)量的增加,IGD(NSGA-Ⅱ)、IGD(SPEA2)和IGD(MOEA/D)增大,而IGD(INSGA-Ⅱ)不變,這說(shuō)明INSGA-Ⅱ能解決不同規(guī)模的分布式柔性作業(yè)車(chē)間調(diào)度問(wèn)題。

    表5 平均IGD值

    本文設(shè)計(jì)的INSGA-Ⅱ之所以能夠獲得更優(yōu)的解,主要因?yàn)椋孩倏紤]加工時(shí)間和能耗的初始化方法能產(chǎn)生更高質(zhì)量的初始種群。該方法采用了新的全局和局部選擇策略來(lái)產(chǎn)生機(jī)器層編碼,均衡加工時(shí)間和能耗,同時(shí),為提高種群多樣性,應(yīng)用隨機(jī)選擇方式產(chǎn)生部分機(jī)器層編碼。②考慮運(yùn)輸時(shí)間的貪婪插入解碼方法通過(guò)將工件插入機(jī)器的空閑時(shí)間內(nèi),合理地安排工件在機(jī)器上的加工順序,能壓縮機(jī)器的空閑時(shí)間,進(jìn)而提高生產(chǎn)效率、降低能耗。③多父代交叉操作可以融合多父代的信息,獲得更高質(zhì)量的子代;新的變異操作不僅避免了算法陷入局部最優(yōu),并且增加了種群的多樣性。④變鄰域搜索策略能進(jìn)一步提高Pareto前沿的質(zhì)量,該策略的嵌入增強(qiáng)了算法的局部搜索能力。為了進(jìn)一步展現(xiàn)本文算法的性能,以算例Mk2-01、Mk2-04、Mk2-09、Mk2-12和Mk2-15為例,繪制4種算法一次求解得到的Pareto前沿分布。

    由圖7可以看出,本文算法Pareto前沿均在另外3種算法Pareto前沿的左下方,表明本文算法可以得到更優(yōu)的Pareto前沿,這與上文的C和IGD評(píng)價(jià)指標(biāo)計(jì)算結(jié)果相符。Pareto前沿分布圖進(jìn)一步表明INSGA-Ⅱ能有效解決考慮運(yùn)輸時(shí)間的分布式柔性作業(yè)車(chē)間綠色調(diào)度問(wèn)題。

    (a)Mk2-01 (b)Mk2-04

    圖8~圖15所示為MOEA/D、NSGA-Ⅱ、SPEA2和INSGA-Ⅱ?qū)λ憷齅k2-09求解得到的調(diào)度方案,可以看出,INSGA-Ⅱ調(diào)度方案的完工時(shí)間最短,并且機(jī)器的總空閑時(shí)間明顯也短于另外3種調(diào)度方案,這意味著INSGA-Ⅱ可以提高機(jī)器的利用率,進(jìn)而提高生產(chǎn)效率、減少能源消耗。

    圖8 Mk2-09的MOEA/D調(diào)度方案(工廠1)

    圖9 Mk2-09的MOEA/D調(diào)度方案(工廠2)

    圖10 Mk2-09的NSGA Ⅱ調(diào)度方案(工廠1)

    圖11 Mk2-09的NSGA Ⅱ調(diào)度方案(工廠2)

    圖12 Mk2-09的SPEA2調(diào)度方案(工廠1)

    圖13 Mk2-09的SPEA2調(diào)度方案(工廠2)

    圖14 Mk2-09的INSGA Ⅱ調(diào)度方案(工廠1)

    圖15 Mk2-09的INSGA Ⅱ調(diào)度方案(工廠2)

    4 結(jié)論

    (1)本文研究了考慮運(yùn)輸時(shí)間的分布式柔性作業(yè)車(chē)間綠色調(diào)度問(wèn)題,建立了最小化最長(zhǎng)完工時(shí)間和總能耗的混合整數(shù)規(guī)劃模型。

    (2)提出了一種改進(jìn)的快速非支配排序遺傳算法,通過(guò)基于工序和機(jī)器的雙層編碼方式、考慮運(yùn)輸時(shí)間的貪婪解碼方法、考慮加工時(shí)間和能耗的種群初始化方法、多父代交叉和變異操作及變鄰域搜索策略來(lái)提高算法的尋優(yōu)能力。

    (3)設(shè)計(jì)了多組考慮運(yùn)輸時(shí)間的分布式柔性作業(yè)車(chē)間綠色調(diào)度問(wèn)題的測(cè)試算例,實(shí)驗(yàn)結(jié)果表明本文設(shè)計(jì)的算法能夠有效解決該問(wèn)題。

    本文假設(shè)運(yùn)輸工具無(wú)限,沒(méi)有考慮運(yùn)輸資源的分配,因此,在未來(lái)的研究中將考慮運(yùn)輸資源的分配和機(jī)器故障等實(shí)際因素。此外,將進(jìn)一步提煉分布式車(chē)間調(diào)度的啟發(fā)式策略,以更好地提高算法性能。

    猜你喜歡
    工序工件機(jī)器
    機(jī)器狗
    120t轉(zhuǎn)爐降低工序能耗生產(chǎn)實(shí)踐
    昆鋼科技(2022年2期)2022-07-08 06:36:14
    機(jī)器狗
    大理石大板生產(chǎn)修補(bǔ)工序詳解(二)
    石材(2020年4期)2020-05-25 07:08:50
    土建工程中關(guān)鍵工序的技術(shù)質(zhì)量控制
    考慮非線性誤差的五軸工件安裝位置優(yōu)化
    未來(lái)機(jī)器城
    電影(2018年8期)2018-09-21 08:00:06
    三坐標(biāo)在工件測(cè)繪中的應(yīng)用技巧
    人機(jī)工程仿真技術(shù)在車(chē)門(mén)裝焊工序中的應(yīng)用
    焊接殘余形變?cè)诠ぜ苎b配中的仿真應(yīng)用研究
    焊接(2015年9期)2015-07-18 11:03:52
    亚洲精品国产色婷婷电影| 国产 一区精品| 黑人猛操日本美女一级片| 国模一区二区三区四区视频| 久久精品国产亚洲av天美| 亚洲国产日韩一区二区| 免费观看a级毛片全部| 亚洲图色成人| 午夜免费男女啪啪视频观看| 欧美精品人与动牲交sv欧美| 免费av中文字幕在线| 亚洲av免费高清在线观看| 国产精品蜜桃在线观看| 欧美日本中文国产一区发布| 少妇被粗大猛烈的视频| 观看美女的网站| 亚洲欧美日韩另类电影网站| 蜜桃在线观看..| 肉色欧美久久久久久久蜜桃| 日韩制服骚丝袜av| 婷婷色综合www| 乱人伦中国视频| 国产精品.久久久| 精品久久久久久久久亚洲| 制服人妻中文乱码| 999精品在线视频| 中国美白少妇内射xxxbb| 亚洲精品乱码久久久久久按摩| 亚洲国产欧美在线一区| 五月天丁香电影| 国产成人精品一,二区| 插逼视频在线观看| 国产精品人妻久久久影院| 成人无遮挡网站| 亚洲人成网站在线播| 看免费成人av毛片| 菩萨蛮人人尽说江南好唐韦庄| 久久精品久久精品一区二区三区| 精品99又大又爽又粗少妇毛片| 欧美人与善性xxx| 18在线观看网站| 亚洲欧洲日产国产| 建设人人有责人人尽责人人享有的| 在线观看免费视频网站a站| 国产成人精品无人区| 大码成人一级视频| 精品一区二区三区视频在线| 18+在线观看网站| 亚洲精品亚洲一区二区| 在线观看美女被高潮喷水网站| 日韩一区二区三区影片| 日韩,欧美,国产一区二区三区| 午夜av观看不卡| 天堂8中文在线网| 国产乱来视频区| 少妇被粗大的猛进出69影院 | 免费观看av网站的网址| 免费久久久久久久精品成人欧美视频 | 婷婷色麻豆天堂久久| 夫妻午夜视频| 婷婷成人精品国产| 亚洲人成77777在线视频| 少妇丰满av| 国产免费视频播放在线视频| 高清在线视频一区二区三区| 人妻夜夜爽99麻豆av| 啦啦啦在线观看免费高清www| 亚洲精品国产av成人精品| 久久久a久久爽久久v久久| 免费观看av网站的网址| 欧美少妇被猛烈插入视频| 秋霞在线观看毛片| 欧美激情极品国产一区二区三区 | 国产综合精华液| 日韩一区二区三区影片| 亚洲中文av在线| 又大又黄又爽视频免费| 夜夜骑夜夜射夜夜干| 久久久久久久亚洲中文字幕| 天天躁夜夜躁狠狠久久av| 免费av中文字幕在线| 亚洲精品成人av观看孕妇| 蜜桃在线观看..| 免费高清在线观看视频在线观看| 国精品久久久久久国模美| 狠狠精品人妻久久久久久综合| 精品酒店卫生间| 精品国产一区二区久久| 成年人午夜在线观看视频| 久久这里有精品视频免费| 母亲3免费完整高清在线观看 | 午夜福利网站1000一区二区三区| 精品久久国产蜜桃| 91精品国产九色| 五月开心婷婷网| 欧美日韩精品成人综合77777| 最近中文字幕高清免费大全6| 国产爽快片一区二区三区| 成人免费观看视频高清| 国产成人精品婷婷| 超碰97精品在线观看| 久久精品国产亚洲av涩爱| 日本色播在线视频| 高清欧美精品videossex| 一级爰片在线观看| 91国产中文字幕| xxxhd国产人妻xxx| 精品亚洲成国产av| 久久久久久久久久久久大奶| 91在线精品国自产拍蜜月| 蜜桃久久精品国产亚洲av| 美女脱内裤让男人舔精品视频| 中国三级夫妇交换| 少妇人妻精品综合一区二区| 尾随美女入室| 69精品国产乱码久久久| 一本一本综合久久| 中国美白少妇内射xxxbb| av免费观看日本| 美女内射精品一级片tv| 国产成人精品一,二区| 美女大奶头黄色视频| av天堂久久9| 久久久精品区二区三区| 全区人妻精品视频| 丰满乱子伦码专区| 91精品国产九色| 黑丝袜美女国产一区| 亚洲av福利一区| 少妇人妻精品综合一区二区| 亚洲经典国产精华液单| 国产精品欧美亚洲77777| 国产亚洲一区二区精品| 免费av不卡在线播放| 桃花免费在线播放| 99热这里只有是精品在线观看| 国产欧美亚洲国产| 国产欧美另类精品又又久久亚洲欧美| 国产视频内射| 免费av不卡在线播放| 精品久久蜜臀av无| 校园人妻丝袜中文字幕| 欧美精品高潮呻吟av久久| 日韩人妻高清精品专区| 国产综合精华液| 丝瓜视频免费看黄片| 狠狠婷婷综合久久久久久88av| 男女边吃奶边做爰视频| 高清毛片免费看| 大话2 男鬼变身卡| 国产精品一区二区在线不卡| 乱人伦中国视频| 寂寞人妻少妇视频99o| 精品久久蜜臀av无| 国产免费视频播放在线视频| 久久久国产欧美日韩av| 插逼视频在线观看| 777米奇影视久久| 免费黄色在线免费观看| videosex国产| 国产精品三级大全| 亚洲精品一区蜜桃| 亚洲精品aⅴ在线观看| 久久午夜综合久久蜜桃| 日日啪夜夜爽| 婷婷色麻豆天堂久久| 91在线精品国自产拍蜜月| 精品亚洲成国产av| 久久久午夜欧美精品| 一二三四中文在线观看免费高清| 大又大粗又爽又黄少妇毛片口| 亚洲精品久久久久久婷婷小说| 国产成人免费无遮挡视频| 涩涩av久久男人的天堂| freevideosex欧美| 性色avwww在线观看| 久久久精品免费免费高清| 69精品国产乱码久久久| 国产精品久久久久成人av| 日韩在线高清观看一区二区三区| 色吧在线观看| 亚洲精品第二区| 日韩熟女老妇一区二区性免费视频| 亚洲精品色激情综合| 久久人人爽人人片av| 国产白丝娇喘喷水9色精品| 啦啦啦在线观看免费高清www| 99热这里只有是精品在线观看| 久久女婷五月综合色啪小说| 成人漫画全彩无遮挡| 国产乱人偷精品视频| 久久免费观看电影| 一区二区三区乱码不卡18| 久久久欧美国产精品| 少妇丰满av| 曰老女人黄片| 亚洲国产日韩一区二区| 亚洲精品国产色婷婷电影| 亚洲美女搞黄在线观看| 91国产中文字幕| 国语对白做爰xxxⅹ性视频网站| 另类亚洲欧美激情| 桃花免费在线播放| 中文字幕av电影在线播放| 亚洲少妇的诱惑av| 日本色播在线视频| 亚洲欧美中文字幕日韩二区| 五月开心婷婷网| 九色成人免费人妻av| 建设人人有责人人尽责人人享有的| 亚洲av二区三区四区| 国产国拍精品亚洲av在线观看| 91精品国产九色| 天堂中文最新版在线下载| 国产日韩一区二区三区精品不卡 | 亚洲av国产av综合av卡| 99久久综合免费| 赤兔流量卡办理| 亚洲怡红院男人天堂| 久久久久国产精品人妻一区二区| 一级毛片黄色毛片免费观看视频| 最近中文字幕高清免费大全6| 久久国产精品男人的天堂亚洲 | 国产日韩一区二区三区精品不卡 | 3wmmmm亚洲av在线观看| 亚洲欧美一区二区三区国产| 免费观看在线日韩| 毛片一级片免费看久久久久| 欧美人与善性xxx| 免费不卡的大黄色大毛片视频在线观看| 少妇熟女欧美另类| 亚洲第一区二区三区不卡| 黄色怎么调成土黄色| 中文字幕久久专区| 国产成人91sexporn| 国产精品一二三区在线看| 免费观看在线日韩| 精品卡一卡二卡四卡免费| 亚洲精华国产精华液的使用体验| 啦啦啦在线观看免费高清www| 女的被弄到高潮叫床怎么办| 免费人妻精品一区二区三区视频| 亚洲四区av| 久久鲁丝午夜福利片| 亚洲欧洲日产国产| 涩涩av久久男人的天堂| 国产极品粉嫩免费观看在线 | 全区人妻精品视频| 高清不卡的av网站| 久久亚洲国产成人精品v| 一本色道久久久久久精品综合| 婷婷成人精品国产| 婷婷色麻豆天堂久久| 午夜福利影视在线免费观看| 九草在线视频观看| 人人澡人人妻人| 亚洲第一av免费看| 成人国产av品久久久| 一区二区三区精品91| 国产精品熟女久久久久浪| 国产精品一区二区三区四区免费观看| 亚洲精品国产色婷婷电影| av在线老鸭窝| 搡女人真爽免费视频火全软件| 国产 精品1| 国产高清不卡午夜福利| 91久久精品国产一区二区成人| 久久久欧美国产精品| 亚州av有码| 18禁观看日本| av免费观看日本| 51国产日韩欧美| 日本免费在线观看一区| av国产久精品久网站免费入址| 欧美+日韩+精品| 精品久久久久久久久av| 丝袜喷水一区| 九九久久精品国产亚洲av麻豆| 啦啦啦啦在线视频资源| 国产精品国产三级专区第一集| 99久国产av精品国产电影| 女性生殖器流出的白浆| 亚洲色图综合在线观看| a级片在线免费高清观看视频| 国产精品秋霞免费鲁丝片| 免费观看的影片在线观看| 曰老女人黄片| 热99国产精品久久久久久7| 久久国产精品男人的天堂亚洲 | 久久狼人影院| 这个男人来自地球电影免费观看 | 日本黄色日本黄色录像| 少妇丰满av| 中文字幕精品免费在线观看视频 | tube8黄色片| 交换朋友夫妻互换小说| 国产精品欧美亚洲77777| 亚洲成人一二三区av| 99re6热这里在线精品视频| 日韩精品有码人妻一区| 国模一区二区三区四区视频| 国产日韩欧美在线精品| 69精品国产乱码久久久| 韩国av在线不卡| 国产免费一级a男人的天堂| 欧美日韩国产mv在线观看视频| 成年人免费黄色播放视频| 亚洲精品国产色婷婷电影| 国产免费又黄又爽又色| 成人国产av品久久久| 亚洲精品自拍成人| 欧美精品亚洲一区二区| 欧美精品人与动牲交sv欧美| 久久国内精品自在自线图片| 免费高清在线观看日韩| 亚洲色图综合在线观看| 欧美日本中文国产一区发布| 亚洲国产精品一区三区| 考比视频在线观看| 午夜久久久在线观看| 国产一区有黄有色的免费视频| 九九久久精品国产亚洲av麻豆| 晚上一个人看的免费电影| 国产精品偷伦视频观看了| 欧美xxxx性猛交bbbb| 国产成人精品无人区| 美女cb高潮喷水在线观看| 国产成人精品在线电影| 王馨瑶露胸无遮挡在线观看| 国产精品.久久久| 久久青草综合色| 麻豆成人av视频| 久久午夜综合久久蜜桃| 亚洲人成网站在线观看播放| 伦理电影免费视频| 另类亚洲欧美激情| 制服人妻中文乱码| 99久国产av精品国产电影| 欧美日韩综合久久久久久| 午夜福利视频精品| .国产精品久久| 卡戴珊不雅视频在线播放| 国产成人一区二区在线| 青春草亚洲视频在线观看| 久久婷婷青草| 久久人人爽人人爽人人片va| 美女中出高潮动态图| 日韩三级伦理在线观看| 在线观看一区二区三区激情| 精品久久蜜臀av无| 夜夜爽夜夜爽视频| 大码成人一级视频| 男人添女人高潮全过程视频| 一级a做视频免费观看| 欧美97在线视频| 国产男人的电影天堂91| 一级毛片aaaaaa免费看小| 国产成人精品婷婷| 日韩一区二区视频免费看| 一级片'在线观看视频| 久久国产精品男人的天堂亚洲 | 51国产日韩欧美| 成人毛片a级毛片在线播放| 亚洲国产色片| 亚洲av中文av极速乱| 久久久久国产精品人妻一区二区| freevideosex欧美| 亚洲av综合色区一区| 亚洲色图 男人天堂 中文字幕 | 国产高清三级在线| 国产成人精品婷婷| 久久毛片免费看一区二区三区| av一本久久久久| 国产不卡av网站在线观看| 久久久久人妻精品一区果冻| 亚洲在久久综合| 最近中文字幕高清免费大全6| 国产在线视频一区二区| 国产精品久久久久久av不卡| 国产av码专区亚洲av| 午夜福利影视在线免费观看| 高清黄色对白视频在线免费看| a级毛片黄视频| 交换朋友夫妻互换小说| 韩国av在线不卡| av播播在线观看一区| 国产精品久久久久久精品古装| 成人黄色视频免费在线看| 女的被弄到高潮叫床怎么办| 精品一区二区免费观看| 免费少妇av软件| 国产亚洲欧美精品永久| kizo精华| 国产av国产精品国产| www.av在线官网国产| 卡戴珊不雅视频在线播放| 夫妻性生交免费视频一级片| 亚洲av免费高清在线观看| 两个人的视频大全免费| 国产日韩欧美视频二区| 欧美日韩精品成人综合77777| 精品99又大又爽又粗少妇毛片| 亚洲美女视频黄频| 又黄又爽又刺激的免费视频.| 国产成人精品婷婷| 日韩在线高清观看一区二区三区| 久久精品久久精品一区二区三区| 亚洲人成77777在线视频| 插阴视频在线观看视频| 欧美精品一区二区大全| 久久99热6这里只有精品| 一区在线观看完整版| 18禁观看日本| 免费不卡的大黄色大毛片视频在线观看| 国产免费又黄又爽又色| 久久久久人妻精品一区果冻| 国产欧美亚洲国产| av电影中文网址| 人妻夜夜爽99麻豆av| 最近的中文字幕免费完整| 国产黄频视频在线观看| 熟妇人妻不卡中文字幕| 国产欧美另类精品又又久久亚洲欧美| av视频免费观看在线观看| 99九九线精品视频在线观看视频| 亚洲中文av在线| 人妻人人澡人人爽人人| 母亲3免费完整高清在线观看 | 校园人妻丝袜中文字幕| 免费高清在线观看日韩| 国内精品宾馆在线| 男女边摸边吃奶| 热re99久久精品国产66热6| 少妇丰满av| videos熟女内射| 国产男女内射视频| av视频免费观看在线观看| 青青草视频在线视频观看| 国产精品.久久久| 亚洲精品国产av成人精品| av网站免费在线观看视频| 天堂俺去俺来也www色官网| 3wmmmm亚洲av在线观看| av在线老鸭窝| 国产高清国产精品国产三级| 精品酒店卫生间| 一本久久精品| 亚洲精品自拍成人| 国产欧美日韩一区二区三区在线 | 国产成人午夜福利电影在线观看| 青春草亚洲视频在线观看| 国产成人精品福利久久| 日本欧美国产在线视频| 免费人妻精品一区二区三区视频| 九九爱精品视频在线观看| 国产一区亚洲一区在线观看| 国产精品国产av在线观看| 最近中文字幕高清免费大全6| 国产男女超爽视频在线观看| 国产精品久久久久久久电影| 国产成人aa在线观看| 国产精品无大码| 国产在线一区二区三区精| 自拍欧美九色日韩亚洲蝌蚪91| 秋霞伦理黄片| 久久久精品区二区三区| 亚洲国产av新网站| 大又大粗又爽又黄少妇毛片口| 中国国产av一级| av有码第一页| 美女国产高潮福利片在线看| 日本av免费视频播放| 高清在线视频一区二区三区| 国国产精品蜜臀av免费| 国产在线免费精品| 欧美 亚洲 国产 日韩一| 精品人妻偷拍中文字幕| 中文字幕亚洲精品专区| 亚洲国产欧美在线一区| 亚洲,欧美,日韩| 成人无遮挡网站| 天天躁夜夜躁狠狠久久av| 91精品国产九色| 日韩大片免费观看网站| a级片在线免费高清观看视频| 久久韩国三级中文字幕| 91久久精品电影网| 蜜桃国产av成人99| 亚州av有码| 久久亚洲国产成人精品v| 丝袜美足系列| 日韩av不卡免费在线播放| 久久久久久久大尺度免费视频| 七月丁香在线播放| 在线播放无遮挡| 美女内射精品一级片tv| 一级,二级,三级黄色视频| 午夜精品国产一区二区电影| 又粗又硬又长又爽又黄的视频| 亚洲性久久影院| 最后的刺客免费高清国语| 亚洲精品自拍成人| 久久精品国产a三级三级三级| 日韩av免费高清视频| 国产在线视频一区二区| 免费日韩欧美在线观看| 成人黄色视频免费在线看| 国产色爽女视频免费观看| 狂野欧美激情性xxxx在线观看| 久久久国产精品麻豆| 丝袜脚勾引网站| 亚洲熟女精品中文字幕| 精品久久久噜噜| 午夜免费鲁丝| 寂寞人妻少妇视频99o| 妹子高潮喷水视频| 日韩三级伦理在线观看| 91精品三级在线观看| 中文字幕精品免费在线观看视频 | 午夜免费鲁丝| 国产免费福利视频在线观看| 嘟嘟电影网在线观看| 大码成人一级视频| 十八禁网站网址无遮挡| 日韩亚洲欧美综合| 欧美人与性动交α欧美精品济南到 | 亚洲综合色网址| 91精品三级在线观看| 大香蕉97超碰在线| 男人操女人黄网站| 中文字幕久久专区| 黄色毛片三级朝国网站| 成人毛片a级毛片在线播放| 亚洲av不卡在线观看| 日韩av免费高清视频| 欧美激情极品国产一区二区三区 | 人妻 亚洲 视频| 91精品伊人久久大香线蕉| 亚洲精品aⅴ在线观看| 美女视频免费永久观看网站| 大片电影免费在线观看免费| 精品酒店卫生间| 国产精品99久久99久久久不卡 | 久久精品熟女亚洲av麻豆精品| 亚洲色图综合在线观看| 男的添女的下面高潮视频| 久久久久久人妻| 久久久精品94久久精品| 色网站视频免费| 天堂俺去俺来也www色官网| 亚洲一级一片aⅴ在线观看| 伦精品一区二区三区| 在现免费观看毛片| 日韩伦理黄色片| 午夜视频国产福利| 日韩欧美精品免费久久| 欧美精品一区二区免费开放| 夜夜爽夜夜爽视频| 少妇猛男粗大的猛烈进出视频| 男女边摸边吃奶| av在线app专区| 青青草视频在线视频观看| 大话2 男鬼变身卡| 精品99又大又爽又粗少妇毛片| 丰满乱子伦码专区| 美女脱内裤让男人舔精品视频| 欧美丝袜亚洲另类| 午夜激情久久久久久久| videosex国产| 99热这里只有是精品在线观看| 国产精品人妻久久久影院| 国产欧美亚洲国产| 欧美xxxx性猛交bbbb| 久久久久久伊人网av| 国产精品99久久久久久久久| 亚洲国产精品成人久久小说| 日本wwww免费看| 色5月婷婷丁香| 久久久久网色| 亚洲怡红院男人天堂| 色5月婷婷丁香| 男人爽女人下面视频在线观看| 亚洲av成人精品一二三区| 熟女av电影| 老司机亚洲免费影院| 少妇人妻 视频| 国产伦理片在线播放av一区| 亚洲精品日本国产第一区| 日韩伦理黄色片| 97超碰精品成人国产| 热99国产精品久久久久久7| 久久精品久久久久久噜噜老黄| 国产精品人妻久久久久久| 插阴视频在线观看视频| 欧美日韩精品成人综合77777| 欧美变态另类bdsm刘玥| 夜夜骑夜夜射夜夜干| 少妇被粗大猛烈的视频| 国产精品人妻久久久久久| 日本欧美国产在线视频| 亚洲欧美成人综合另类久久久| 国产伦理片在线播放av一区| 国产片特级美女逼逼视频| 九九久久精品国产亚洲av麻豆| 欧美人与善性xxx| 亚洲欧美一区二区三区国产| 最近最新中文字幕免费大全7| 看免费成人av毛片| 少妇熟女欧美另类| 国产成人精品婷婷| 成人亚洲欧美一区二区av| 香蕉精品网在线| 狂野欧美激情性xxxx在线观看| 国产免费福利视频在线观看| 国产在线一区二区三区精|