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

    基于人工蜂群算法的柔性工藝與車間調(diào)度集成優(yōu)化

    2017-04-20 05:38:02宋栓軍楊佩莉石雯麗
    計(jì)算機(jī)應(yīng)用 2017年2期
    關(guān)鍵詞:車間工序柔性

    宋栓軍,楊佩莉,石雯麗

    (西安工程大學(xué) 機(jī)電工程學(xué)院,西安 710048)

    (*通信作者電子郵箱447090664@qq.com)

    基于人工蜂群算法的柔性工藝與車間調(diào)度集成優(yōu)化

    宋栓軍,楊佩莉*,石雯麗

    (西安工程大學(xué) 機(jī)電工程學(xué)院,西安 710048)

    (*通信作者電子郵箱447090664@qq.com)

    為實(shí)現(xiàn)柔性工藝與車間調(diào)度集成優(yōu)化,在考慮工件特征的加工工藝、次序及加工機(jī)器的柔性基礎(chǔ)上,以最小化最大完工時(shí)間為優(yōu)化目標(biāo),提出一種基于交叉變異的人工蜂群算法。該算法針對(duì)柔性工藝與車間調(diào)度集成問題的離散性特征,對(duì)工藝路線進(jìn)行序列編碼,工件調(diào)度采用基于工序的編碼方式。通過工藝種群與調(diào)度種群的交叉變異操作,分別使采蜜蜂及觀察蜂進(jìn)行局部尋優(yōu),偵查蜂進(jìn)行全局尋優(yōu),以此提高算法性能。在此基礎(chǔ)上用兩部分測(cè)試實(shí)例分別驗(yàn)證了集成研究的必要性及改進(jìn)算法的有效性。

    柔性工藝規(guī)劃;車間調(diào)度;人工蜂群算法

    0 引言

    工藝規(guī)劃與車間調(diào)度集成(Integrated Process Planning and Scheduling,IPPS)是制造系統(tǒng)中急需解決的關(guān)鍵問題之一。在傳統(tǒng)的企業(yè)生產(chǎn)過程中,工藝規(guī)劃和車間調(diào)度往往被當(dāng)作獨(dú)立的系統(tǒng)進(jìn)行研究,導(dǎo)致了工藝與調(diào)度的脫節(jié)[1]。如果將工藝規(guī)劃和車間調(diào)度系統(tǒng)進(jìn)行集成,則可避免制造資源的浪費(fèi),優(yōu)化生產(chǎn)節(jié)拍,提高制造系統(tǒng)的工作效率[2-3]。

    Chryssolouris等[4]首先提出了工藝規(guī)劃與車間調(diào)度集成的構(gòu)想;之后Beckendorff等[5]在二者集成時(shí)考慮了工藝路線的柔性問題,首次在IPPS問題中考慮到工藝的柔性特點(diǎn)。上述研究使學(xué)者們開始關(guān)注工藝與調(diào)度集成問題的研究。劉倩雯[6]利用人工蜂群算法(Artificial Bee Colony, ABC)求解作業(yè)車間調(diào)度問題,通過仿真實(shí)驗(yàn)驗(yàn)證了算法的有效性;瞿璨[7]將人工蜂群與模擬退火算法相結(jié)合求解面向綠色制造的工藝路線優(yōu)化,通過算例測(cè)試獲得了高效穩(wěn)定的解。這兩個(gè)研究都是將工藝與調(diào)度系統(tǒng)進(jìn)行獨(dú)立研究。呂盛坪等[8-9]對(duì)工藝規(guī)劃與車間調(diào)度集成現(xiàn)狀和趨勢(shì)以及集成模型作了分析,為IPPS問題的后續(xù)研究奠定基礎(chǔ);崔曉康等[10]提出三階段工藝與調(diào)度的集成方法,并通過實(shí)例驗(yàn)證了方法的有效性;王進(jìn)峰等[11]對(duì)蟻群算法中信息素?fù)]發(fā)速度以及螞蟻轉(zhuǎn)移概率進(jìn)行改進(jìn),將其應(yīng)用于IPPS優(yōu)化中,最后通過仿真實(shí)例表明了算法的有效性。這四個(gè)研究都是針對(duì)單目標(biāo)IPPS問題,并且沒有考慮工藝的柔性特征。Liu等[12]提出一種多目標(biāo)微型人工蜂群算法,對(duì)有約束條件的多目標(biāo)彈性作業(yè)車間調(diào)度問題求解,優(yōu)化目標(biāo)是盡量減少機(jī)器的工作量以及生產(chǎn)周期;文笑雨[13]采用蜜蜂繁殖優(yōu)化算法,分別對(duì)多目標(biāo)及多目標(biāo)不確定性IPPS問題進(jìn)行研究,之后研究了多目標(biāo)IPPS決策方法,最后通過實(shí)例驗(yàn)證算法的有效性。這兩篇文章都是針對(duì)多目標(biāo)IPPS問題進(jìn)行研究。從上述研究看,雖然針對(duì)IPPS的研究較多,但考慮到工藝柔性的較少;人工蜂群算法本身具備全局收斂優(yōu)勢(shì)[14],通過三種蜂的轉(zhuǎn)換,借助啟發(fā)式搜索策略,具有很好的局部搜索及全局尋優(yōu)的能力,因此本文在研究以上方法的基礎(chǔ)上,考慮到工件的柔性特征,提出了一種基于交叉變異的人工蜂群算法,對(duì)柔性工藝與車間調(diào)度集成(Integrated Flexible Process Planning and Scheduling,IFPPS)問題進(jìn)行優(yōu)化求解,并通過工藝的序列編碼及基于工序的調(diào)度編碼使算法能夠適應(yīng)IFPPS問題的離散性特征,提出工藝種群與調(diào)度種群的交叉變異操作,不僅使產(chǎn)生的新解能夠滿足工序的順序約束,還能夠通過此算法最終獲得更優(yōu)的加工路線和調(diào)度方案。

    1 IFPPS優(yōu)化數(shù)學(xué)模型

    本文以單目標(biāo)IFPPS問題為例,以加工任務(wù)的最大完工時(shí)間(makespan)最小為目標(biāo)進(jìn)行優(yōu)化研究,車間調(diào)度問題為作業(yè)車間調(diào)度。在建立優(yōu)化數(shù)學(xué)模型前作如下假設(shè)[15-16]:

    1)初始時(shí)刻,制造車間的所有機(jī)器均可正常使用;

    2)在進(jìn)行排產(chǎn)優(yōu)化時(shí),不同工件相互之間沒有優(yōu)先級(jí),不同工件的加工工序之間也沒有優(yōu)先級(jí),同一工件的各個(gè)工序之間存在優(yōu)先級(jí);

    3)在同一時(shí)間,每臺(tái)機(jī)器僅可加工一道工序;

    4)所有工序的加工時(shí)間是指其實(shí)際加工與其準(zhǔn)備加工時(shí)間之和。

    基于以上假設(shè),本文研究的IFPPS數(shù)學(xué)模型如下所示,其中:A表示一個(gè)非常大的正數(shù);n表示需要加工的工件總數(shù);m表示機(jī)器數(shù);gi表示工件i可選工藝路線數(shù);pil表示工件i的第l條工藝路線中的工序數(shù);oijl表示工件i在第l條工藝路線下的第j道工序;k表示工序oijl的可選加工機(jī)器;bijlk表示工序oijl在機(jī)器k上的開工時(shí)間;tijlk表示工序oijl在機(jī)器k上的加工時(shí)間;cijlk表示工序oijl在機(jī)器k上的完工時(shí)間;ci(j-1)lq表示工序oijl的上一道加工工序oi(j-1)l在其加工機(jī)器q上的完工時(shí)間;cxylk表示機(jī)器k上排在工序oijl緊前加工的工序oxyl的完工時(shí)間;zi(j-1,j)l(q,k)表示工序oi(j-1)l的加工機(jī)器q轉(zhuǎn)向工序oijl的加工機(jī)器k的機(jī)器轉(zhuǎn)換時(shí)間;ci表示工件i的完工時(shí)間;vijlk表示工序oijl在機(jī)器k上的加工成本;

    Yijlpqsk=

    目標(biāo)函數(shù)如式(1)所示:

    f1=min{makespan=max(ci)}

    (1)

    約束條件如式(2)~(10)所示:

    1)工件i第l條工藝路線的首道加工工序的最早完工時(shí)間大于或等于其加工時(shí)間:

    (ci1lk×Zi1lk×Xil)+A(1-Xil)≥(ti1lk×Zi1lk×Xil);

    ?i∈[1,n], ?l∈[1,gi], ?k∈[1,m]

    (2)

    2)工件i第l條工藝路線的末道加工工序的最早完工時(shí)間小于或等于工件i的最大完工時(shí)間:

    (cipillk×Zipillk×Xil)-A(1-Xil)≤makespan;

    ?i∈[1,n], ?l∈[1,gi], ?k∈[1,m]

    (3)

    3)一個(gè)工件同一時(shí)刻僅可有一道工序被加工:

    (cijlk×Zijlk×Xil)-(ci(j-1)lk1×Zi(j-1)lk1×Xil)+

    A(1-Xil)≥(tijlk×Zijlk×Xil); ?i∈[1,n],

    ?j∈[1,pil], ?l∈[1,gi], ?k,k1∈[1,m]

    (4)

    4)每臺(tái)機(jī)器同一時(shí)刻僅可加工一道工序:

    (5)

    5)工件i僅可選取一條工藝路線輸入到調(diào)度系統(tǒng):

    (6)

    6)一道工序僅可選取一臺(tái)加工機(jī)器:

    (7)

    7)所有工序的開工時(shí)間大于或等于0:

    bijlk×Zijlk×Xil≥0; ?i∈[1,n], ?j∈[1,pil],

    ?l∈[1,gi], ?k∈[1,m]

    (8)

    8)工序oijl在機(jī)器k上的完工時(shí)間等于其開工時(shí)間與加工時(shí)間之和:

    cijlk×Zijlk×Xil=(bijlk+tijlk×Zijlk×Xil; ?i∈[1,n],

    ?j∈[1,pil], ?l∈[1,gi], ?k∈[1,m]

    (9)

    9)工序oijl在機(jī)器k上的開工時(shí)間等于工序oi(j-1)l在其加工機(jī)器q上的完工時(shí)間與工序oi(j-1)l的加工機(jī)器q轉(zhuǎn)向工序oijl的加工機(jī)器k的機(jī)器轉(zhuǎn)換時(shí)間之和以及之前最后一個(gè)在加工機(jī)器k上進(jìn)行加工的工序完工時(shí)間中較大的數(shù):

    bijlk=max(oi(j-1)lq+zi(j-1, j),l(q,k),cxylk); ?i∈[1,n],

    ?j∈[2,pil], ?l∈[1,gi], ?k,q∈[1,m]

    (10)

    2 IFPPS優(yōu)化算法設(shè)計(jì)

    2.1 柔性工藝調(diào)度的編碼及解碼

    1)柔性工藝規(guī)劃編碼及解碼。

    根據(jù)工藝的柔性特征,本文創(chuàng)新性地提出一種序列編碼方式對(duì)工藝規(guī)劃部分的各個(gè)序列進(jìn)行編碼,表1為工件1的加工工藝信息。對(duì)工件1加工工藝信息進(jìn)行序列編碼,編碼結(jié)果如圖1所示。

    表1 工件1加工工藝信息表

    圖1 工件1加工工藝信息序列編碼方案

    特征序列中的第n個(gè)數(shù)字m表示特征m的加工順序是n;工序序列中的第a個(gè)數(shù)字b表示工件第a個(gè)特征選擇第b種可選工藝加工;機(jī)器序列中的第i個(gè)數(shù)字j表示第i個(gè)工藝用j號(hào)機(jī)器進(jìn)行加工。用序列編碼方式很容易實(shí)現(xiàn)解碼,首先根據(jù)特征序列,可以獲得工件1各個(gè)特征的加工順序?yàn)镕4—F2—F1—F6—F3—F5;然后根據(jù)工序序列獲得各特征所選的加工工藝,再按照特征順序進(jìn)行排列,可以得到工件1具體的加工工序序列為O5—O2—O1—O7—O4—O6;最后根據(jù)機(jī)器序列,結(jié)合加工工藝得到各個(gè)工藝所選擇的加工機(jī)器,從而確定出該工件的一條可行工藝路線O5(M3)—O2(M4)—O1(M1)—O7(M1)—O4(M1)—O6(M4)。

    2)車間調(diào)度編碼及解碼。

    在車間調(diào)度編碼的相關(guān)研究中,用得最多的是基于工序的編碼方法,本文也利用此方式對(duì)個(gè)體進(jìn)行編碼[17]。設(shè)一個(gè)企業(yè)制造車間中的待加工工件數(shù)為4,在初始工藝路線中按照適應(yīng)度函數(shù)值挑選出各工件較優(yōu)的路線后,就可以知道要加工的工序總數(shù)。假設(shè)工件1共包含6道工序,工件2包含5道工序,工件3包含3道工序,工件4包含5道工序, 那么[1 1 2 3 4 4 2 4 1 1 2 4 2 1 1 4 2 3 3]就為一種可行的調(diào)度編碼方案。在此方案中,工件序號(hào)出現(xiàn)的次數(shù)與其所包含的工序數(shù)相同,例如工件1包含6道工序,那么1就在調(diào)度編碼方案中一共出現(xiàn)6次。此外方案中第6個(gè)位置上的4是從左到右數(shù)數(shù)字4出現(xiàn)的第2次,所以第6個(gè)位置代表了工件4的工藝路線中第2道要加工的工序。

    由于本文選用的優(yōu)化目標(biāo)為最小化最大完工時(shí)間,屬于正規(guī)調(diào)度指標(biāo),并且基于工序的編碼方式產(chǎn)生的調(diào)度方案均為可行的調(diào)度方案,因此選擇使用文獻(xiàn)[18]中的貪婪解碼方法將編碼方案解碼為活動(dòng)調(diào)度類型。貪婪解碼方法的具體步驟可由文獻(xiàn)[18]獲得,此處不再贅述。

    2.2 基于交叉變異的人工蜂群算法

    使用人工蜂群算法求解IFPPS優(yōu)化問題時(shí),每個(gè)蜜源代表一個(gè)可行的調(diào)度方案,每個(gè)蜜源的參數(shù)從開始到結(jié)束代表某個(gè)任務(wù)中各工序可能的工作順序,因此可由改變蜜源參數(shù)來改變各工序的加工順序,并計(jì)算其完工時(shí)間,尋找最優(yōu)方案[14]。為了保證種群的多樣性,本文對(duì)基本人工蜂群算法進(jìn)行改進(jìn),提出了一種基于交叉變異的人工蜂群算法,并設(shè)計(jì)了算法求解的主要內(nèi)容。

    1)蜂群初始化。

    由于工藝及工件順序的表示需要是一個(gè)大于或等于1的整數(shù),而基本人工蜂群算法中,種群初始化方法是產(chǎn)生連續(xù)的隨機(jī)數(shù),會(huì)出現(xiàn)非整數(shù)的情況,因此本文對(duì)人工蜂群算法的初始化方式進(jìn)行了改進(jìn)。通過Matlab程序直接產(chǎn)生符合要求的可行性隨機(jī)解,在工藝規(guī)劃階段,由于工序有次序約束,因此隨機(jī)生成的特征序列可能不符合要求,采用文獻(xiàn)[18]中的約束調(diào)整方法對(duì)特征序列進(jìn)行調(diào)整,使生成的各工件的初始工藝路線可行。由于基于工序的編碼方式不會(huì)生成不可行解,因此本文對(duì)車間調(diào)度部分的初始化蜂群進(jìn)行隨機(jī)生成,然后將編碼方案解碼后,計(jì)算蜂群中每個(gè)蜜蜂所代表的IFPPS方案的完工時(shí)間。

    2)工藝種群交叉變異操作。

    為了使工藝規(guī)劃系統(tǒng)可以源源不斷地為調(diào)度系統(tǒng)輸入各個(gè)工件不同的工藝路線方案,在進(jìn)行柔性工藝規(guī)劃與車間調(diào)度集成優(yōu)化時(shí),集成優(yōu)化總流程每迭代一次,就需要對(duì)工藝種群更新一次。針對(duì)工藝的序列編碼方式,本文提出工藝種群交叉變異操作對(duì)其進(jìn)行種群更新,保證種群的多樣性。

    在工藝種群交叉變異操作中,使采蜜蜂和觀察蜂利用變異操作進(jìn)行局部搜索。具體步驟如下:

    步驟1 將特征序列中兩個(gè)隨機(jī)位置的數(shù)值進(jìn)行交換,產(chǎn)生新的特征序列,并通過約束調(diào)整方法對(duì)其進(jìn)行調(diào)整,使其符合特征之間的次序約束條件;

    步驟2 將工藝序列中一個(gè)隨機(jī)位置所選擇的加工工藝轉(zhuǎn)變成另一種可選加工工藝;

    步驟3 將機(jī)器序列中一個(gè)隨機(jī)位置所選擇的加工機(jī)器轉(zhuǎn)變成另一種可選加工機(jī)器。

    在工藝種群交叉變異操作中,使偵查蜂利用交叉操作進(jìn)行全局搜索。當(dāng)偵查蜂進(jìn)行全局搜索隨機(jī)產(chǎn)生新蜜源時(shí),由于算法初始化已經(jīng)產(chǎn)生了符合工件特征順序約束的序列,優(yōu)化過程就是適當(dāng)調(diào)整這個(gè)序列中參數(shù)的順序,以得到一個(gè)能夠使某個(gè)目標(biāo)值達(dá)到最優(yōu)的順序,而不需要重新生成任何參數(shù)或信息,因此采用兩點(diǎn)交叉操作來對(duì)其進(jìn)行調(diào)整,能夠使產(chǎn)生的新序列滿足工件的特征次序約束。在偵查蜂進(jìn)行交叉操作時(shí),分別對(duì)工件的三個(gè)序列進(jìn)行兩點(diǎn)交叉操作,產(chǎn)生新的可行工藝路線。在文獻(xiàn)[2]的基礎(chǔ)上,對(duì)其特征序列交叉操作進(jìn)行改進(jìn),將原始的復(fù)制交叉點(diǎn)中斷的參數(shù)信息到新的特征序列相應(yīng)位置改為復(fù)制交叉點(diǎn)兩端的參數(shù)信息到新的特征序列的相應(yīng)位置,此方式可以確保生成的新的可行解滿足工序的次序約束條件,具體步驟如下:

    步驟1 隨機(jī)挑選兩個(gè)蜜源(蜜源1、蜜源2),對(duì)其三個(gè)序列分別進(jìn)行交叉操作。

    步驟2 在蜜源的序列中隨機(jī)選擇兩個(gè)交叉點(diǎn),特征序列用步驟3交叉操作,工藝序列和機(jī)器序列用步驟4進(jìn)行交叉操作。

    步驟3 將特征序列1交叉點(diǎn)兩端參數(shù)信息復(fù)制給新特征序列的相應(yīng)位置,刪掉特征序列2中新特征序列中已有的特征數(shù)值,將特征序列2中剩余的數(shù)值按順序依次填入到新特征序列的空位置中,這樣就生成了新的特征序列。

    步驟4 將工藝序列1交叉點(diǎn)中斷參數(shù)信息復(fù)制到新工藝序列的相應(yīng)位置,將工藝序列2的交叉點(diǎn)兩端的參數(shù)信息按順序填入新工藝序列的空位置中,這樣就生成了新的工藝序列,機(jī)器序列同工藝序列。

    3)調(diào)度種群交叉變異操作。

    在調(diào)度系統(tǒng)的優(yōu)化流程迭代中,需要不斷尋找新的調(diào)度方案,對(duì)調(diào)度種群進(jìn)行更新。針對(duì)基于工序的調(diào)度編碼方式,本文提出一種調(diào)度種群交叉變異操作進(jìn)行種群更新。

    在調(diào)度種群交叉變異操作中,使采蜜蜂及觀察蜂利用變異操作進(jìn)行局部搜索。即在原來的調(diào)度方案中,隨機(jī)挑選兩個(gè)不同的參數(shù),將其位置進(jìn)行互換,即可得到變異后的調(diào)度方案。

    在調(diào)度種群交叉變異操作中,使偵查蜂利用交叉操作進(jìn)行全局搜索。本文采用基于工序先后順序的交叉操作[2],具體步驟如下:

    步驟1 在調(diào)度種群中隨機(jī)挑選兩個(gè)蜜源(可行解);

    步驟2 將工件集{1,2,…,n}隨機(jī)劃分為兩個(gè)非空子集JobSet1和JobSet2;

    步驟3 將蜜源1、蜜源2中與JobSet1相同的工件分別復(fù)制到新蜜源1、新蜜源2的對(duì)應(yīng)位置,并保持順序不變;

    步驟4 將蜜源1、蜜源2中與JobSet2相同的工件分別復(fù)制到新蜜源2、新蜜源1中,并保持順序不變。

    以圖1的編碼方案為例,對(duì)上述蜜源進(jìn)行交叉操作,結(jié)果如圖2所示。

    圖2 基于工序先后順序的交叉操作

    3 IFPPS優(yōu)化流程設(shè)計(jì)

    本文采用集成式優(yōu)化策略,以最大完工時(shí)間最小化為優(yōu)化目標(biāo),由于人工蜂群算法中適應(yīng)度值越大代表蜜源越好,因此使適應(yīng)度函數(shù)等于目標(biāo)函數(shù)的倒數(shù),如式(11)所示,f是目標(biāo)函數(shù)值,F(xiàn)itness是適應(yīng)度函數(shù)值。

    Fitness=1/f

    (11)

    根據(jù)上述基于交叉變異的人工蜂群算法,對(duì)工藝規(guī)劃系統(tǒng)和車間調(diào)度系統(tǒng)進(jìn)行集成優(yōu)化,最終確定各個(gè)工件的最優(yōu)工藝路線和車間調(diào)度方案,具體流程如圖3所示。

    圖3 優(yōu)化策略流程

    求解步驟如下:

    步驟1 設(shè)加工工件數(shù)為n,根據(jù)各工件的工藝信息,初始化n個(gè)工藝種群。

    步驟2 計(jì)算每個(gè)工件初始路線的適應(yīng)度函數(shù)值,選擇值最大的一條工藝路線。

    步驟3 把各工件所選擇的工藝路線輸入車間調(diào)度系統(tǒng)。

    步驟4 通過基于交叉變異的人工蜂群算法對(duì)車間調(diào)度系統(tǒng)進(jìn)行優(yōu)化,主要包含以下步驟:

    1)根據(jù)工藝系統(tǒng)輸送的工藝路線,初始化調(diào)度種群。

    2)計(jì)算種群中所有個(gè)體的適應(yīng)度值,并將值最大的調(diào)度方案記為最優(yōu)調(diào)度方案。

    3)判斷是否達(dá)到車間調(diào)度系統(tǒng)優(yōu)化的終止條件(達(dá)到調(diào)度系統(tǒng)最大迭代次數(shù)時(shí)算法結(jié)束):如果是,就輸出最優(yōu)調(diào)度方案和各工件對(duì)應(yīng)的工藝路線,并將其記為最優(yōu)解,之后進(jìn)行步驟5;不是就進(jìn)行步驟4中的4)。

    4)根據(jù)本文2.2節(jié),利用調(diào)度種群的交叉變異操作產(chǎn)生新的調(diào)度種群。

    5)轉(zhuǎn)步驟4中的2),迭代次數(shù)加1。

    步驟5 根據(jù)基本ABC算法的貪婪準(zhǔn)則更新最優(yōu)解。

    步驟6 判斷是否滿足集成系統(tǒng)的終止條件(達(dá)到集成系統(tǒng)最大運(yùn)行次數(shù)時(shí)算法結(jié)束),如果是,就輸出最優(yōu)解,不是就轉(zhuǎn)到步驟7。

    步驟7 根據(jù)本文2.2節(jié)工藝種群交叉變異操作方法,更新所有工藝種群。

    步驟8 轉(zhuǎn)到步驟2,迭代次數(shù)加1。

    4 實(shí)例驗(yàn)證

    本文進(jìn)行了兩組實(shí)例測(cè)試:第一組是為了驗(yàn)證柔性工藝規(guī)劃與車間調(diào)度集成問題的必要性,對(duì)基于交叉變異的人工蜂群算法用于求解柔性工藝規(guī)劃與車間調(diào)度集成和非集成的柔性工藝規(guī)劃與車間調(diào)度(NonIFPPS)的結(jié)果進(jìn)行比較;第二組是為了評(píng)估本文基于交叉變異的人工蜂群算法求解該集成問題的有效性,將本文算法求解結(jié)果與目前較好的其他算法進(jìn)行比較。兩部分實(shí)例算法參數(shù)設(shè)置如下:初始化蜜蜂總數(shù)NP=100(工藝系統(tǒng)與調(diào)度系統(tǒng)相同);采蜜蜂數(shù)量FoodNumber=NP/2(工藝系統(tǒng)與調(diào)度系統(tǒng)相同);偵查蜂數(shù)量SearchNumber=5(工藝系統(tǒng)與調(diào)度系統(tǒng)相同);集成最大運(yùn)行次數(shù)RunTime=200;調(diào)度系統(tǒng)最大迭代次數(shù)maxCycle=50;蜜源最大開采次數(shù)Limit=5(工藝系統(tǒng)與調(diào)度系統(tǒng)相同)。

    4.1 集成必要性驗(yàn)證

    為了驗(yàn)證研究柔性工藝規(guī)劃與車間調(diào)度集成問題的必要性,以六種工件在五臺(tái)機(jī)器上加工的IFPPS問題為例進(jìn)行分析,工件的參數(shù)信息如表1~3所示。通過Matlab對(duì)算法進(jìn)行編程,柔性工藝規(guī)劃與車間調(diào)度集成優(yōu)化與非集成優(yōu)化具體的測(cè)試結(jié)果對(duì)比如表4所示,集成式最優(yōu)調(diào)度方案甘特圖如圖4所示。

    從對(duì)比結(jié)果可以看出,單就工藝路線優(yōu)化而言,非集成式中的工件1、工件2、工件4及工件5的最優(yōu)工藝路線對(duì)應(yīng)的加工時(shí)間小于集成式中對(duì)應(yīng)的加工時(shí)間,而集成式柔性工藝規(guī)劃與車間調(diào)度的最優(yōu)調(diào)度方案的完工時(shí)間小于非集成式柔性工藝規(guī)劃與車間調(diào)度。這是因?yàn)樵趥鹘y(tǒng)的車間生產(chǎn)系統(tǒng)中,車間調(diào)度方案往往是在工藝規(guī)劃之后進(jìn)行的,在執(zhí)行調(diào)度時(shí),工藝規(guī)劃必須限定在之前所產(chǎn)生的工藝路線之中,這樣車間調(diào)度不可避免會(huì)受到工藝路線的約束,從而就會(huì)限制調(diào)度系統(tǒng)產(chǎn)生優(yōu)化方案的效果。對(duì)柔性工藝規(guī)劃與車間調(diào)度進(jìn)行集成優(yōu)化研究可以同時(shí)考慮工件的工藝路線優(yōu)化和車間調(diào)度優(yōu)化,不僅能夠增加車間生產(chǎn)系統(tǒng)的柔性,同時(shí)能夠減少工藝路線對(duì)調(diào)度優(yōu)化的限制,因此即使單個(gè)工件的工藝路線不是最優(yōu),但是整體的調(diào)度方案卻是最優(yōu)的。由此說明研究柔性工藝規(guī)劃與車間調(diào)度集成優(yōu)化問題是十分必要的。

    4.2 算法有效性驗(yàn)證

    由于現(xiàn)有單目標(biāo)工藝規(guī)劃與車間調(diào)度集成優(yōu)化研究中,缺乏同時(shí)考慮工藝規(guī)劃多種柔性特征的實(shí)例,因此,選用Morad等[19]文章中的一個(gè)實(shí)例,該實(shí)例是5個(gè)工件在5臺(tái)機(jī)器上進(jìn)行加工,每個(gè)工件不存在加工工藝柔性及特征次序柔性,只存在加工機(jī)器柔性,并且工件在不同機(jī)器之間的轉(zhuǎn)換時(shí)間忽略不計(jì)。具體工件的工藝信息如表5所示。

    通過Matlab對(duì)實(shí)例進(jìn)行編程,本文基于交叉變異的人工蜂群算法求得的優(yōu)化結(jié)果為makespan=31,與其他優(yōu)化算法求解結(jié)果對(duì)比如表6所示。從表6中可以看出,本文人工蜂群算法的優(yōu)化結(jié)果優(yōu)于其他算法,主要原因在于本文針對(duì)柔性工藝規(guī)劃與車間調(diào)度集成優(yōu)化問題設(shè)計(jì)了較好的編碼方法與搜索策略,使得算法在求解效率及搜索能力上得到增強(qiáng),因此驗(yàn)證了本文人工蜂群算法的有效性及優(yōu)越性。本文算法求得的最優(yōu)調(diào)度方案甘特圖如圖5所示。

    表2 工件的加工工藝信息表(測(cè)試1)

    表3 工件在不同機(jī)器之間的轉(zhuǎn)換時(shí)間(測(cè)試1)

    表4 集成必要性驗(yàn)證測(cè)試結(jié)果(測(cè)試1)

    圖4 集成式最優(yōu)調(diào)度方案甘特圖(測(cè)試1)

    表5 工件的工藝信息[19](測(cè)試2)

    Tab.5Processinformationofworkpieces(test2)

    工件加工特征加工工藝可選加工機(jī)器對(duì)應(yīng)加工時(shí)間工件加工特征加工工藝可選加工機(jī)器對(duì)應(yīng)加工時(shí)間1F1O1M1,M25,3F2O2M27F3O3M26F4O4M4,M53,44F1O1M2,M32,6F2O2M38F3O3M3,M43,8F4O4M4,M57,42F1O1M17F2O2M2,M34,6F3O3M3,M47,7F4O4M5105F1O1M1,M33,5F2O2M37F3O3M4,M59,6F4O4M533F1O1M1,M2,M34,5,8F2O2M45F3O3M2,M36,5F4O4M54

    圖5 車間調(diào)度甘特圖(makespan=31) (測(cè)試2)

    表6 算法有效性驗(yàn)證測(cè)試結(jié)果(測(cè)試2)

    Tab.6 Algorithm verification results (test 2)

    方法makespan方法makespanHeuristic[19]38IGA[1]33SA[19]33ABC31Agent?based[20]33

    5 結(jié)語

    為了使制造車間能夠及時(shí)響應(yīng)訂單的個(gè)性化需求,消除缺乏柔性的缺陷,提高其生產(chǎn)效率,本文提出一種基于交叉變異的人工蜂群算法對(duì)柔性工藝規(guī)劃與車間調(diào)度集成進(jìn)行優(yōu)化求解,重點(diǎn)對(duì)算法的具體操作及工藝與調(diào)度集成優(yōu)化流程進(jìn)行了設(shè)計(jì),最后選取最小化最大完工時(shí)間為優(yōu)化目標(biāo),通過兩部分實(shí)例對(duì)集成研究的必要性以及算法求解的有效性進(jìn)行了驗(yàn)證。

    References)

    [1] 高亮,李新宇.工藝規(guī)劃與車間調(diào)度集成研究現(xiàn)狀及進(jìn)展[J].中國機(jī)械工程,2011,22(8): 1001-1007.(GAO L, LI X Y.Current research on integrated process planning and scheduling [J].China Mechanical Engineering, 2011, 22(8): 1001-1007.)

    [2] 宋栓軍,楊佩莉,石雯麗.考慮多目標(biāo)的柔性工藝與調(diào)度集成優(yōu)化算法[J/OL].計(jì)算機(jī)應(yīng)用研究,2017 [2016- 08- 05].http://www.arocmag.com/getarticle?aid=912497cbe6174f12.(SONG S J, YANG P L, SHI W L.Integrated optimization algorithm of flexible process planning and shop scheduling with consideration of multi-objectives [J/OL].Application Research of Computers, 2017 [2016- 08- 05].http://www.arocmag.com/getarticle?aid=912497cbe6174f12.

    [3] 巴黎,李言,楊明順,等.考慮不確定加工時(shí)間的工藝規(guī)劃與調(diào)度集成問題研究[J].中國機(jī)械工程,2015,26(24):3348-3355.(BA L, LI Y, YANG M S, et al.Research on integrated process planning and scheduling problem with consideration of uncertain processing time [J].China Mechanical Engineering, 2015, 26(24): 3348-3355.)

    [4] CHRYSSOLOURIS G, CHAN S, SUH N P.An integrated approach to process planning and scheduling [J].CIRP Annals — Manufacturing Technology, 1985, 34(1): 413- 417.

    [5] BECKENDORFF U, KREUTZFELDT J.Reactive workshop scheduling based on alternative routings [C]// Proceedings of a Conference on Factory Automation and Information Management.Bocaraton, Florida: CRC Press, 1991: 875-885.

    [6] 劉倩雯.人工蜂群算法及其在調(diào)度問題中的應(yīng)用研究[D].北京: 北京交通大學(xué),2014: 5-61.(LIU Q W.Research on artificial bee colony and its application in scheduling problems [D].Beijing: Beijing Jiaotong University, 2014: 5-61.)

    [7] 瞿璨.面向綠色制造的工藝規(guī)劃決策研究[D].武漢: 華中科技大學(xué),2013:2-54.(QU C.Study on the process planning for green manufacturing [D].Wuhan: Huazhong University of Science and Technology, 2013: 2-54.)

    [8] 呂盛坪,喬立紅.工藝規(guī)劃與車間調(diào)度及兩者集成的研究現(xiàn)狀和發(fā)展趨勢(shì)[J].計(jì)算機(jī)集成制造系統(tǒng),2014,20(2):290-300.(LYU S P, QIAO L H.Current status and developing trend of process planning and job shop scheduling [J].Computer Integrated Manufacturing Systems, 2014, 20(2): 290-300.)

    [9] 呂盛坪,喬立紅.工藝規(guī)劃與車間調(diào)度復(fù)合式集成模型[J].計(jì)算機(jī)集成制造系統(tǒng),2014,20(1):111-120.(LYU S P, QIAO L H.Hybrid integration model for process planning and scheduling [J].Computer Integrated Manufacturing Systems, 2014, 20(1): 111-120.)

    [10] 崔曉康,肖艷秋,王磊,等.面向工藝規(guī)劃與車間調(diào)度的集成方法研究[J].機(jī)床與液壓,2015,43(3): 134-137.(CUI X K, XIAO Y Q, WANG L, et al.Research on integrated method of process planning and job-shop scheduling [J].Machine Tool & Hydraulics, 2015, 43(3): 134-137.)

    [11] 王進(jìn)峰,陰國富,雷前召,等.蟻群算法在工藝規(guī)劃與車間調(diào)度集成優(yōu)化中的應(yīng)用[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2012,42(S1): 173-177.(WANG J F, YIN G F, LEI Q Z, et al.Integrated process planning and scheduling based on an ant colony algorithm [J].Journal of Southeast Universit (Natural Science Edition), 2012, 42(S1): 173-177.)

    [12] LIU Z, MA S, SHI Y, et al.Solving multi-objective flexible job shop scheduling with transportation constraints using a micro artificial bee colony algorithm [C]// Proceedings of the 2013 IEEE 17th International Conference on Computer Supported Cooperative Work in Design.Piscataway, NJ: IEEE, 2013: 427-432.

    [13] 文笑雨.多目標(biāo)集成式工藝規(guī)劃與車間調(diào)度問題的求解方法研究[D].武漢:華中科技大學(xué),2014: 5-61.(WEN X Y.Research on the solution methods for multi-objective integrated process planning and scheduling problem [D].Wuhan: Huazhong University of Science and Technology, 2014: 5-61.)

    [14] 江銘炎,袁東風(fēng).人工蜂群算法及其應(yīng)用[M].北京:科學(xué)出版,2014:16-264.(JIANG M Y, YUAN D F.Artificial Bee Colony Algorithm and Its Application [M].Beijing: Science Press, 2014: 16-264.)

    [15] 吳正佳,羅月勝,周玉瓊,等.一種求解典型JSP的改進(jìn)離散粒子群優(yōu)化算法[J].計(jì)算機(jī)應(yīng)用研究,2013,30(8):2405-2409.(WU Z J, LUO Y S, ZHOU Y Q, et al.Improved discrete particle swarm optimization algorithm for typical Job-Shop scheduling problem [J].Application Research of Computers, 2013, 30(8): 2405-2409.)

    [16] 張國輝,張凌杰,吳立輝,等.精英進(jìn)化策略求解柔性作業(yè)車間調(diào)度問題[J].計(jì)算機(jī)應(yīng)用研究,2016,33(12): 1-5.(ZHANG G H, ZHANG L J, WU L H, et al.Solving flexible Job-Shop scheduling with elite evolution strategy [J].Application Research of Computers, 2016, 33(12): 1-5.)

    [17] 李新宇.工藝規(guī)劃與車間調(diào)度集成問題的求解方法研究[D].武漢:華中科技大學(xué),2009:58.(LI X Y.Research on the solution methods of integrated process planning and scheduling [D].Wuhan: Huazhong University of Science and Technology, 2009: 58.)[18] FANG H-L, ROSS P, CORNE D.A promising genetic algorithm approach to job-shop scheduling, re-scheduling and open-shop scheduling problems [C]// Proceedings of the 5th International Conference on Genetic Algorithms.San Mateo, CA: Morgan Kaufmann Publishers, 1993: 375-382.

    [19] MORAD N, ZALZALA A M S.Genetic algorithms in integrated process planning and scheduling [J].Journal of Intelligent Manufacturing, 1999, 10(2): 169-179.

    [20] FUJI N, INOUE R, UEDA K.Integration of process planning and scheduling using multi-Agent learning [C]// Proceedings of the 41st CIRP Conference on Manufacturing Systems.Berlin: Springer-Verlag, 2008: 297-300.

    This work is partially supported by the Scientific Research Fund by Shaanxi Provincial Education Department (15JK1311), Research Start-up Fund for Doctor of Xi’an Polytechnic University (BS1301), Innovation Fund for Graduate Students of Xi’an Polytechnic University (CX201628).

    SONG Shuanjun, born in 1974, Ph.D., associate professor.His research interests include optimization of production system, supply chain management.

    YANG Peili, born in 1992, M.S.candidate.Her research interests include optimization of production system.

    SHI Wenli, born in 1991, M.S.candidate.Her research interests include supply chain management.

    Optimization of integrated flexible process planning and job shop scheduling based on artificial bee colony

    SONG Shuanjun, YANG Peili*, SHI Wenli

    (CollegeofMechanicalandElectricalEngineering,Xi’anPolytechnicUniversity,Xi’anShaanxi710048,China)

    To achieve the optimization of integrated flexible process planning and job shop scheduling, taking the flexibility of manufacturing process and order and manufacturing machine of the workpieces into account, for minimizing the maximum completion time of the product processing task, an artificial bee colony algorithm based on crossover and mutation was proposed.Aiming at the discrete characteristics of integrated flexible process and job shop scheduling, the process route was coded in sequence, and the job scheduling was based on the working procedure.To improve the performance of the algorithm, by means of crossover and mutation operation of process population and scheduling population, the employed foragers and onlookers bees seeked local optimality, and the scouts seeked global optimality.On this basis, the necessity of the integration research and the effectiveness of the improved algorithm were verified by two test cases.

    flexible process planning; job shop scheduling; Artificial Bee Colony (ABC)

    2016- 07- 18;

    2016- 08- 24。 基金項(xiàng)目:陜西省教育廳科研基金資助項(xiàng)目(15JK1311);西安工程大學(xué)博士科研啟動(dòng)基金資助項(xiàng)目(BS1301);西安工程大學(xué)研究生創(chuàng)新基金資助項(xiàng)目(CX201628)。

    宋栓軍(1974—),男,陜西西安人,副教授,博士,主要研究方向:生產(chǎn)系統(tǒng)優(yōu)化、供應(yīng)鏈管理; 楊佩莉(1992—),女,陜西興平人,碩士研究生,主要研究方向:生產(chǎn)系統(tǒng)優(yōu)化; 石雯麗(1991—),女,甘肅隴南人,碩士研究生,主要研究方向:供應(yīng)鏈管理。

    1001- 9081(2017)02- 0523- 07

    10.11772/j.issn.1001- 9081.2017.02.0523

    TH166;TP

    A

    猜你喜歡
    車間工序柔性
    一種柔性拋光打磨頭設(shè)計(jì)
    120t轉(zhuǎn)爐降低工序能耗生產(chǎn)實(shí)踐
    昆鋼科技(2022年2期)2022-07-08 06:36:14
    灌注式半柔性路面研究進(jìn)展(1)——半柔性混合料組成設(shè)計(jì)
    石油瀝青(2021年5期)2021-12-02 03:21:18
    100MW光伏車間自動(dòng)化改造方案設(shè)計(jì)
    智能制造(2021年4期)2021-11-04 08:54:28
    高校學(xué)生管理工作中柔性管理模式應(yīng)用探索
    大理石大板生產(chǎn)修補(bǔ)工序詳解(二)
    石材(2020年4期)2020-05-25 07:08:50
    土建工程中關(guān)鍵工序的技術(shù)質(zhì)量控制
    招工啦
    “扶貧車間”拔窮根
    把農(nóng)業(yè)搬進(jìn)車間
    亚洲综合色惰| 亚洲人成网站高清观看| 国产黄片美女视频| 美女大奶头视频| 国产一区二区三区在线臀色熟女| 国产真实乱freesex| 亚洲国产精品久久男人天堂| 亚洲欧美日韩东京热| 老熟妇仑乱视频hdxx| 99久久中文字幕三级久久日本| 99国产精品一区二区蜜桃av| 99久久精品国产国产毛片| 国产精品久久视频播放| 亚洲av二区三区四区| 欧美激情久久久久久爽电影| 一级黄色大片毛片| 精品人妻偷拍中文字幕| 日韩强制内射视频| 亚洲av免费在线观看| 免费大片18禁| bbb黄色大片| 亚洲第一电影网av| 国产精品人妻久久久久久| 亚洲国产高清在线一区二区三| 亚洲欧美日韩无卡精品| 波多野结衣高清作品| 搡老妇女老女人老熟妇| 少妇猛男粗大的猛烈进出视频 | 少妇裸体淫交视频免费看高清| 日韩强制内射视频| 亚洲av第一区精品v没综合| 97超级碰碰碰精品色视频在线观看| 啦啦啦韩国在线观看视频| 亚洲第一区二区三区不卡| 亚洲av.av天堂| 一区二区三区免费毛片| 在线播放国产精品三级| 麻豆成人av在线观看| 国产主播在线观看一区二区| 亚洲国产精品成人综合色| 国产v大片淫在线免费观看| av在线亚洲专区| 人人妻,人人澡人人爽秒播| 综合色av麻豆| 日本撒尿小便嘘嘘汇集6| 午夜视频国产福利| 嫩草影视91久久| 国国产精品蜜臀av免费| 国产成人aa在线观看| 可以在线观看毛片的网站| 日本黄大片高清| 国产高清视频在线观看网站| 国产精品嫩草影院av在线观看 | 国产午夜精品久久久久久一区二区三区 | 1000部很黄的大片| 国产av不卡久久| 麻豆国产av国片精品| 国产亚洲精品久久久com| 久久人妻av系列| 国产亚洲精品久久久久久毛片| 国产精品久久久久久精品电影| 亚洲五月天丁香| 免费黄网站久久成人精品| 伦理电影大哥的女人| 搡女人真爽免费视频火全软件 | 麻豆国产av国片精品| 99热精品在线国产| 嫩草影视91久久| 噜噜噜噜噜久久久久久91| 在线免费观看不下载黄p国产 | 国产国拍精品亚洲av在线观看| 变态另类丝袜制服| 两性午夜刺激爽爽歪歪视频在线观看| 又黄又爽又免费观看的视频| 婷婷精品国产亚洲av| 九九在线视频观看精品| 动漫黄色视频在线观看| 日韩精品青青久久久久久| 欧美bdsm另类| 亚洲国产欧洲综合997久久,| 丰满乱子伦码专区| 一个人观看的视频www高清免费观看| 午夜老司机福利剧场| 熟女人妻精品中文字幕| 亚洲av不卡在线观看| 久久久久久久亚洲中文字幕| 在线观看舔阴道视频| 在线a可以看的网站| av福利片在线观看| 九九在线视频观看精品| 国产av在哪里看| 国内少妇人妻偷人精品xxx网站| 日韩中字成人| 亚洲国产欧美人成| 很黄的视频免费| 欧美日韩国产亚洲二区| 久久精品久久久久久噜噜老黄 | 级片在线观看| 日日夜夜操网爽| 欧美不卡视频在线免费观看| 国产爱豆传媒在线观看| 最好的美女福利视频网| 国产免费av片在线观看野外av| 成年女人永久免费观看视频| 欧美在线一区亚洲| 少妇人妻精品综合一区二区 | 国产真实伦视频高清在线观看 | 欧美高清性xxxxhd video| 十八禁国产超污无遮挡网站| 色5月婷婷丁香| 精品一区二区三区视频在线观看免费| 日日啪夜夜撸| 亚洲综合色惰| 99九九线精品视频在线观看视频| 亚洲五月天丁香| bbb黄色大片| 又爽又黄a免费视频| 人人妻人人看人人澡| 国产精品亚洲一级av第二区| 久久亚洲精品不卡| 岛国在线免费视频观看| 黄色丝袜av网址大全| 午夜福利在线观看免费完整高清在 | 欧美区成人在线视频| 国产欧美日韩精品亚洲av| 人妻少妇偷人精品九色| 熟妇人妻久久中文字幕3abv| 99久久精品热视频| 最近在线观看免费完整版| 人人妻人人看人人澡| 91av网一区二区| 女人十人毛片免费观看3o分钟| 久久久久久久久久久丰满 | 免费搜索国产男女视频| 国产色爽女视频免费观看| 精品人妻偷拍中文字幕| 国产高清有码在线观看视频| 国产精品永久免费网站| 精品一区二区三区视频在线观看免费| 久久久久久大精品| 国产精品久久久久久久久免| 波野结衣二区三区在线| 长腿黑丝高跟| 亚洲一区高清亚洲精品| 波野结衣二区三区在线| 少妇的逼好多水| 亚洲欧美激情综合另类| 日韩国内少妇激情av| 午夜a级毛片| 国产亚洲av嫩草精品影院| 国产欧美日韩精品一区二区| 色在线成人网| 免费搜索国产男女视频| 亚洲在线观看片| 99国产极品粉嫩在线观看| 丝袜美腿在线中文| 国产精品电影一区二区三区| 国产精品福利在线免费观看| 免费搜索国产男女视频| 男女那种视频在线观看| 波多野结衣高清无吗| 亚洲av熟女| 久久人妻av系列| 此物有八面人人有两片| 欧美在线一区亚洲| 精品一区二区三区av网在线观看| 日韩欧美精品v在线| 欧美xxxx黑人xx丫x性爽| 直男gayav资源| 老熟妇仑乱视频hdxx| 男女之事视频高清在线观看| 又黄又爽又刺激的免费视频.| 大型黄色视频在线免费观看| 一区二区三区四区激情视频 | 日本欧美国产在线视频| 亚洲,欧美,日韩| 看黄色毛片网站| 久久久久久九九精品二区国产| 18禁在线播放成人免费| 国产色爽女视频免费观看| 精品一区二区三区视频在线观看免费| 国产伦精品一区二区三区视频9| 听说在线观看完整版免费高清| 亚洲avbb在线观看| 日韩一本色道免费dvd| 狠狠狠狠99中文字幕| 最近最新中文字幕大全电影3| 日韩欧美 国产精品| or卡值多少钱| 日本撒尿小便嘘嘘汇集6| 中文字幕高清在线视频| 一个人看的www免费观看视频| 精品人妻1区二区| 精品久久国产蜜桃| 国产淫片久久久久久久久| 尤物成人国产欧美一区二区三区| 日韩精品中文字幕看吧| 亚洲精品一区av在线观看| 久久久久久国产a免费观看| 男人的好看免费观看在线视频| 精品久久久久久久久久久久久| 国产高清不卡午夜福利| 国产蜜桃级精品一区二区三区| www.www免费av| 国产av麻豆久久久久久久| 此物有八面人人有两片| 久久久久国内视频| 成人永久免费在线观看视频| 国产成年人精品一区二区| 久久99热6这里只有精品| 熟女电影av网| 少妇裸体淫交视频免费看高清| 久久人妻av系列| 国产精品无大码| 男人和女人高潮做爰伦理| 亚洲最大成人中文| 国产精品日韩av在线免费观看| 九九爱精品视频在线观看| 在线观看av片永久免费下载| 精品久久久久久久末码| 国产伦一二天堂av在线观看| 99热只有精品国产| 国产精品伦人一区二区| 日本免费一区二区三区高清不卡| 最近最新免费中文字幕在线| 91在线观看av| 深夜精品福利| 亚洲av电影不卡..在线观看| 97碰自拍视频| 嫁个100分男人电影在线观看| 欧美+日韩+精品| 亚洲精品色激情综合| 老司机午夜福利在线观看视频| 国产一级毛片七仙女欲春2| 亚洲av第一区精品v没综合| 精品国内亚洲2022精品成人| 一区福利在线观看| 久久久久性生活片| 亚洲一区二区三区色噜噜| av视频在线观看入口| 可以在线观看毛片的网站| 久久婷婷人人爽人人干人人爱| 99热6这里只有精品| 精品不卡国产一区二区三区| 亚洲av第一区精品v没综合| 国产精品电影一区二区三区| 99久久中文字幕三级久久日本| 免费看av在线观看网站| 两个人的视频大全免费| 日韩欧美免费精品| 欧美另类亚洲清纯唯美| 可以在线观看毛片的网站| 老熟妇仑乱视频hdxx| 国产黄色小视频在线观看| 久久久精品大字幕| 九九久久精品国产亚洲av麻豆| 日本撒尿小便嘘嘘汇集6| 永久网站在线| 中亚洲国语对白在线视频| 亚洲成av人片在线播放无| 国产精品国产三级国产av玫瑰| 女的被弄到高潮叫床怎么办 | 国产精品伦人一区二区| 此物有八面人人有两片| 日韩亚洲欧美综合| 变态另类成人亚洲欧美熟女| 欧美一区二区亚洲| 看片在线看免费视频| 91麻豆精品激情在线观看国产| 成人二区视频| av天堂在线播放| 在现免费观看毛片| 日本成人三级电影网站| 69av精品久久久久久| 亚洲中文日韩欧美视频| 国产男人的电影天堂91| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| 一级黄片播放器| 身体一侧抽搐| 春色校园在线视频观看| 日本熟妇午夜| 免费看美女性在线毛片视频| 看片在线看免费视频| 亚洲精品色激情综合| 内地一区二区视频在线| 18禁黄网站禁片免费观看直播| 欧美xxxx性猛交bbbb| 夜夜看夜夜爽夜夜摸| 免费搜索国产男女视频| 亚洲国产精品合色在线| 淫秽高清视频在线观看| 国产伦精品一区二区三区四那| 在线免费观看的www视频| 蜜桃久久精品国产亚洲av| 日本与韩国留学比较| 麻豆成人午夜福利视频| 一区二区三区激情视频| 日韩高清综合在线| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av一区综合| 美女被艹到高潮喷水动态| 免费不卡的大黄色大毛片视频在线观看 | 美女 人体艺术 gogo| 又黄又爽又免费观看的视频| 亚洲综合色惰| x7x7x7水蜜桃| 国产伦一二天堂av在线观看| 日韩欧美精品v在线| 少妇猛男粗大的猛烈进出视频 | 婷婷精品国产亚洲av在线| 国产黄色小视频在线观看| 色尼玛亚洲综合影院| 国产一区二区三区视频了| 九色成人免费人妻av| 日本撒尿小便嘘嘘汇集6| 成人国产综合亚洲| 国产亚洲91精品色在线| 色视频www国产| 亚洲经典国产精华液单| 国产一区二区三区在线臀色熟女| 波多野结衣巨乳人妻| 国产av不卡久久| 久久久久久久亚洲中文字幕| 亚洲av美国av| 十八禁国产超污无遮挡网站| 免费不卡的大黄色大毛片视频在线观看 | 日韩中字成人| 性欧美人与动物交配| 亚洲精品在线观看二区| 很黄的视频免费| 久久精品人妻少妇| 非洲黑人性xxxx精品又粗又长| 成人毛片a级毛片在线播放| 日本欧美国产在线视频| 三级国产精品欧美在线观看| 亚洲欧美日韩高清专用| 12—13女人毛片做爰片一| 嫁个100分男人电影在线观看| 在线观看66精品国产| 久久亚洲真实| 国产高清不卡午夜福利| 国产高清有码在线观看视频| 又黄又爽又免费观看的视频| 欧美成人性av电影在线观看| 99久久九九国产精品国产免费| 校园人妻丝袜中文字幕| 亚洲国产日韩欧美精品在线观看| 18禁在线播放成人免费| 国产精品99久久久久久久久| 91在线精品国自产拍蜜月| 老司机深夜福利视频在线观看| 18禁在线播放成人免费| 老司机深夜福利视频在线观看| 免费看美女性在线毛片视频| 日韩国内少妇激情av| 天美传媒精品一区二区| 黄色一级大片看看| 欧美黑人巨大hd| 国语自产精品视频在线第100页| 久久久久久久精品吃奶| 少妇高潮的动态图| 毛片一级片免费看久久久久 | 一级毛片久久久久久久久女| 欧美人与善性xxx| 国产精品女同一区二区软件 | 午夜福利在线观看吧| 国产白丝娇喘喷水9色精品| 亚洲国产精品成人综合色| 伦理电影大哥的女人| 免费搜索国产男女视频| 欧美日韩综合久久久久久 | 看免费成人av毛片| 欧美在线一区亚洲| 日本成人三级电影网站| 一本一本综合久久| 男女啪啪激烈高潮av片| 国产69精品久久久久777片| 精品久久久久久,| 久久草成人影院| 日本色播在线视频| 亚洲精品456在线播放app | 淫妇啪啪啪对白视频| 亚洲自偷自拍三级| a级一级毛片免费在线观看| 欧美日韩中文字幕国产精品一区二区三区| 欧美日本视频| 国产伦人伦偷精品视频| 国语自产精品视频在线第100页| 亚洲精华国产精华精| 日韩精品中文字幕看吧| 亚洲最大成人av| 久久久久精品国产欧美久久久| 午夜影院日韩av| 性色avwww在线观看| 在线天堂最新版资源| 毛片女人毛片| 国产爱豆传媒在线观看| 国产亚洲欧美98| 国产精品无大码| 亚洲欧美精品综合久久99| 最新中文字幕久久久久| 久久人妻av系列| 亚洲熟妇中文字幕五十中出| 看黄色毛片网站| 美女xxoo啪啪120秒动态图| 欧美人与善性xxx| 狠狠狠狠99中文字幕| 白带黄色成豆腐渣| 一本久久中文字幕| 在线观看av片永久免费下载| 亚洲一区高清亚洲精品| 亚洲美女视频黄频| 成人美女网站在线观看视频| 久久久久久大精品| 亚洲va日本ⅴa欧美va伊人久久| 成人午夜高清在线视频| 久久香蕉精品热| 日韩亚洲欧美综合| 99国产极品粉嫩在线观看| 欧美3d第一页| 久久亚洲精品不卡| 国产色爽女视频免费观看| 国产av不卡久久| 人人妻人人看人人澡| 亚洲av一区综合| 伦理电影大哥的女人| www.www免费av| av在线蜜桃| 热99re8久久精品国产| 精品人妻1区二区| 女人被狂操c到高潮| av在线观看视频网站免费| 日韩 亚洲 欧美在线| 九色国产91popny在线| 免费无遮挡裸体视频| 嫩草影院新地址| 一个人看视频在线观看www免费| 免费看a级黄色片| 国产精品美女特级片免费视频播放器| 国产亚洲精品久久久久久毛片| 亚洲av免费在线观看| 99在线视频只有这里精品首页| 黄色配什么色好看| 国产精品女同一区二区软件 | 久久精品国产亚洲av涩爱 | 91狼人影院| 欧美一区二区国产精品久久精品| 精品久久久久久久久久免费视频| 无遮挡黄片免费观看| 美女 人体艺术 gogo| 伊人久久精品亚洲午夜| or卡值多少钱| 久久久国产成人精品二区| 日韩精品青青久久久久久| 三级毛片av免费| 又粗又爽又猛毛片免费看| videossex国产| 久久久久久久久久成人| 天堂√8在线中文| 精品一区二区三区av网在线观看| 亚州av有码| 精品人妻一区二区三区麻豆 | 国内少妇人妻偷人精品xxx网站| 免费观看的影片在线观看| 欧美日韩综合久久久久久 | 亚洲不卡免费看| 少妇人妻一区二区三区视频| x7x7x7水蜜桃| 成人国产一区最新在线观看| 日韩 亚洲 欧美在线| 不卡一级毛片| 免费观看的影片在线观看| 亚洲成人免费电影在线观看| 日本a在线网址| 日韩在线高清观看一区二区三区 | 成人性生交大片免费视频hd| 两个人的视频大全免费| 高清毛片免费观看视频网站| 国产精品98久久久久久宅男小说| 免费搜索国产男女视频| av在线蜜桃| 美女高潮的动态| www.色视频.com| 美女被艹到高潮喷水动态| 天美传媒精品一区二区| 日本黄色视频三级网站网址| xxxwww97欧美| 久久精品国产亚洲av香蕉五月| 国产探花在线观看一区二区| 欧美日韩乱码在线| 日本黄色视频三级网站网址| 简卡轻食公司| 91久久精品电影网| 别揉我奶头~嗯~啊~动态视频| 三级国产精品欧美在线观看| 男插女下体视频免费在线播放| 国产爱豆传媒在线观看| 色播亚洲综合网| 欧美三级亚洲精品| 国产日本99.免费观看| 午夜a级毛片| 联通29元200g的流量卡| 国产黄色小视频在线观看| 国产亚洲精品久久久久久毛片| 午夜激情欧美在线| 久久精品国产亚洲av香蕉五月| 国产欧美日韩精品一区二区| 欧美激情在线99| 国产午夜福利久久久久久| 色视频www国产| 91麻豆精品激情在线观看国产| 久久久久久久久中文| 国产精品久久久久久av不卡| 不卡一级毛片| 中出人妻视频一区二区| 十八禁网站免费在线| 国产伦一二天堂av在线观看| 国产精华一区二区三区| av在线观看视频网站免费| 久久久久精品国产欧美久久久| 欧美性猛交黑人性爽| 国产黄片美女视频| 国产成人av教育| 搡女人真爽免费视频火全软件 | 国产伦在线观看视频一区| 18禁裸乳无遮挡免费网站照片| 欧美又色又爽又黄视频| 午夜久久久久精精品| 最近视频中文字幕2019在线8| 一级黄色大片毛片| 欧美色欧美亚洲另类二区| 久久精品夜夜夜夜夜久久蜜豆| 日本一本二区三区精品| 校园人妻丝袜中文字幕| 三级国产精品欧美在线观看| 精品日产1卡2卡| 99在线人妻在线中文字幕| 欧美日韩瑟瑟在线播放| 在线免费十八禁| 男女下面进入的视频免费午夜| 国产av麻豆久久久久久久| 少妇猛男粗大的猛烈进出视频 | 亚洲欧美日韩东京热| 男人和女人高潮做爰伦理| 美女免费视频网站| 国产探花在线观看一区二区| 91久久精品国产一区二区成人| 午夜福利高清视频| 国内久久婷婷六月综合欲色啪| a级毛片a级免费在线| 在线观看免费视频日本深夜| 亚洲成a人片在线一区二区| 亚洲一级一片aⅴ在线观看| 国产成人一区二区在线| 久久久色成人| 精品久久久久久久久亚洲 | 亚洲欧美日韩卡通动漫| 日韩欧美 国产精品| 内地一区二区视频在线| 狂野欧美激情性xxxx在线观看| 欧美区成人在线视频| 美女高潮的动态| 波多野结衣高清作品| 国产精品国产高清国产av| 久99久视频精品免费| 黄色女人牲交| 国产精品不卡视频一区二区| 午夜免费激情av| 97人妻精品一区二区三区麻豆| 亚洲三级黄色毛片| 一个人免费在线观看电影| 日本黄大片高清| 精品99又大又爽又粗少妇毛片 | h日本视频在线播放| 韩国av在线不卡| 亚洲精品色激情综合| 久久久久国产精品人妻aⅴ院| 18禁黄网站禁片午夜丰满| 男女做爰动态图高潮gif福利片| netflix在线观看网站| 午夜免费激情av| www日本黄色视频网| 欧美潮喷喷水| 毛片一级片免费看久久久久 | 午夜精品久久久久久毛片777| 国产主播在线观看一区二区| 91在线观看av| 精品欧美国产一区二区三| 黄片wwwwww| 看黄色毛片网站| 99久久精品国产国产毛片| 久久精品国产自在天天线| 国产精品人妻久久久久久| 给我免费播放毛片高清在线观看| 一a级毛片在线观看| 亚洲午夜理论影院| 久久久色成人| 婷婷六月久久综合丁香| 久久久久免费精品人妻一区二区| 亚洲欧美日韩高清专用| 国产伦精品一区二区三区视频9| 中国美白少妇内射xxxbb| 久9热在线精品视频| 久久久成人免费电影| 国产免费男女视频| 国产人妻一区二区三区在| 午夜福利欧美成人| 精品国内亚洲2022精品成人| 午夜福利在线观看吧| 成人三级黄色视频| 日韩一区二区视频免费看| 亚洲va在线va天堂va国产| 久久精品国产自在天天线| 国产精品不卡视频一区二区|