• <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)車間
    欧美日韩国产亚洲二区| 亚洲在线自拍视频| 精品福利观看| 日本 av在线| 最近视频中文字幕2019在线8| 国产精品精品国产色婷婷| 国产女主播在线喷水免费视频网站 | 亚洲精品粉嫩美女一区| 男女视频在线观看网站免费| 亚洲精品乱码久久久v下载方式| 麻豆一二三区av精品| 在线观看免费视频日本深夜| 黄色欧美视频在线观看| 一个人观看的视频www高清免费观看| 2021天堂中文幕一二区在线观| 丰满的人妻完整版| 免费观看的影片在线观看| 国产精品久久电影中文字幕| 亚洲av二区三区四区| 国产 一区精品| 亚洲av美国av| 精品不卡国产一区二区三区| 成人欧美大片| 能在线免费观看的黄片| 老女人水多毛片| 国产白丝娇喘喷水9色精品| 美女大奶头视频| 亚洲性夜色夜夜综合| 在线免费十八禁| 五月玫瑰六月丁香| 国产精品国产三级国产av玫瑰| 99久久无色码亚洲精品果冻| 国产精品永久免费网站| 在线天堂最新版资源| av国产免费在线观看| 久久久久久久久大av| 免费看a级黄色片| 一进一出抽搐动态| 99国产精品一区二区蜜桃av| 亚洲精品456在线播放app| 中文资源天堂在线| 久久久午夜欧美精品| 亚州av有码| 成人毛片a级毛片在线播放| 欧美绝顶高潮抽搐喷水| 啦啦啦韩国在线观看视频| 网址你懂的国产日韩在线| 97人妻精品一区二区三区麻豆| 亚洲内射少妇av| 在现免费观看毛片| 看片在线看免费视频| 99视频精品全部免费 在线| 欧美国产日韩亚洲一区| 一卡2卡三卡四卡精品乱码亚洲| 啦啦啦啦在线视频资源| 欧美3d第一页| 啦啦啦韩国在线观看视频| 一进一出抽搐gif免费好疼| 长腿黑丝高跟| 一级黄片播放器| 国产在线男女| 日韩国内少妇激情av| 免费一级毛片在线播放高清视频| 亚洲精品456在线播放app| 亚洲成人久久爱视频| 国产精品久久久久久久久免| 午夜视频国产福利| 日产精品乱码卡一卡2卡三| 淫秽高清视频在线观看| 色av中文字幕| 欧美高清成人免费视频www| 女人被狂操c到高潮| 亚洲欧美日韩高清专用| 69av精品久久久久久| 亚洲av五月六月丁香网| 麻豆av噜噜一区二区三区| 最近2019中文字幕mv第一页| 给我免费播放毛片高清在线观看| 国产一区二区激情短视频| 久久久久久久久大av| 国产免费一级a男人的天堂| 国产午夜精品论理片| 久久久精品大字幕| 午夜久久久久精精品| 国产精品国产高清国产av| 国产一区二区三区av在线 | 日本在线视频免费播放| 国产黄色视频一区二区在线观看 | 日韩高清综合在线| 欧美潮喷喷水| 国产精品福利在线免费观看| 亚洲七黄色美女视频| 午夜精品国产一区二区电影 | 22中文网久久字幕| 国产av不卡久久| 最近2019中文字幕mv第一页| 国产伦在线观看视频一区| 亚洲成人精品中文字幕电影| 亚洲婷婷狠狠爱综合网| 亚洲av中文字字幕乱码综合| 欧美三级亚洲精品| 成熟少妇高潮喷水视频| 嫩草影院新地址| 国产成人aa在线观看| 黄色日韩在线| av免费在线看不卡| 最近2019中文字幕mv第一页| 99久久九九国产精品国产免费| www日本黄色视频网| 麻豆久久精品国产亚洲av| 国产精品亚洲美女久久久| 日韩成人伦理影院| 亚洲欧美日韩无卡精品| 久久午夜福利片| 真实男女啪啪啪动态图| 神马国产精品三级电影在线观看| 亚洲av熟女| 日本三级黄在线观看| 少妇高潮的动态图| 一个人看视频在线观看www免费| 成年免费大片在线观看| 久久精品久久久久久噜噜老黄 | 欧美日本视频| 男女做爰动态图高潮gif福利片| 亚洲欧美清纯卡通| 级片在线观看| 欧美一级a爱片免费观看看| 一区二区三区免费毛片| 成年版毛片免费区| 国产真实伦视频高清在线观看| 嫩草影院新地址| 久久精品国产鲁丝片午夜精品| 日韩成人伦理影院| 国产综合懂色| 老师上课跳d突然被开到最大视频| 真实男女啪啪啪动态图| 夜夜看夜夜爽夜夜摸| 深夜精品福利| 国产v大片淫在线免费观看| 网址你懂的国产日韩在线| 久久亚洲精品不卡| 校园人妻丝袜中文字幕| 91久久精品国产一区二区成人| 悠悠久久av| 日日撸夜夜添| 俺也久久电影网| 欧美极品一区二区三区四区| 国产精品人妻久久久久久| 99久久久亚洲精品蜜臀av| 国产探花极品一区二区| 成人三级黄色视频| 免费看美女性在线毛片视频| 免费观看的影片在线观看| 午夜精品国产一区二区电影 | 亚洲人成网站高清观看| 校园春色视频在线观看| 免费看av在线观看网站| 大香蕉久久网| 日韩强制内射视频| 午夜久久久久精精品| 日韩人妻高清精品专区| 久久久久久大精品| 国产精品日韩av在线免费观看| 日本与韩国留学比较| 自拍偷自拍亚洲精品老妇| 成人特级av手机在线观看| 99热这里只有是精品在线观看| 网址你懂的国产日韩在线| 亚洲熟妇熟女久久| 亚洲精品乱码久久久v下载方式| 麻豆久久精品国产亚洲av| 成人漫画全彩无遮挡| 一级黄片播放器| 日韩亚洲欧美综合| 国产真实乱freesex| 久久久久国内视频| 中文资源天堂在线| 国产av不卡久久| 女人被狂操c到高潮| 波野结衣二区三区在线| 最近视频中文字幕2019在线8| 九九爱精品视频在线观看| 国内精品美女久久久久久| 国国产精品蜜臀av免费| 成人国产麻豆网| 色噜噜av男人的天堂激情| 亚洲最大成人av| 国产精品伦人一区二区| 日韩欧美 国产精品| 久久综合国产亚洲精品| 国产探花在线观看一区二区| 少妇被粗大猛烈的视频| 国产麻豆成人av免费视频| 久久鲁丝午夜福利片| 亚洲成人中文字幕在线播放| 日韩欧美 国产精品| 成人毛片a级毛片在线播放| 日日撸夜夜添| 老司机福利观看| 麻豆久久精品国产亚洲av| 日韩成人av中文字幕在线观看 | 在线免费观看的www视频| 非洲黑人性xxxx精品又粗又长| 看十八女毛片水多多多| 国产伦精品一区二区三区视频9| 人妻夜夜爽99麻豆av| 人人妻人人澡欧美一区二区| 亚洲最大成人av| 国产乱人偷精品视频| 日韩一本色道免费dvd| 日韩成人av中文字幕在线观看 | 精品久久久久久久久av| 国产精品福利在线免费观看| 欧美日本亚洲视频在线播放| 久久精品影院6| 韩国av在线不卡| 九九在线视频观看精品| 国产黄色小视频在线观看| av在线老鸭窝| 亚洲美女视频黄频| 美女高潮的动态| 成年女人看的毛片在线观看| 亚洲精品国产成人久久av| 高清日韩中文字幕在线| 搡女人真爽免费视频火全软件 | 一本精品99久久精品77| 日本色播在线视频| 亚洲国产精品成人久久小说 | 日韩中字成人| 韩国av在线不卡| 在线天堂最新版资源| 亚洲无线在线观看| www日本黄色视频网| 成人毛片a级毛片在线播放| 亚洲成人精品中文字幕电影| 亚洲va在线va天堂va国产| 91狼人影院| 久久久久久久久大av| 美女高潮的动态| 亚洲中文字幕一区二区三区有码在线看| 丰满的人妻完整版| 国产精品一及| 男女做爰动态图高潮gif福利片| 一进一出好大好爽视频| 精品无人区乱码1区二区| 99视频精品全部免费 在线| 午夜亚洲福利在线播放| 在线播放无遮挡| 综合色av麻豆| 亚洲精品久久国产高清桃花| 国产片特级美女逼逼视频| 亚洲婷婷狠狠爱综合网| 99热精品在线国产| 久久精品国产清高在天天线| 成人亚洲精品av一区二区| 国产高清视频在线播放一区| 一本久久中文字幕| 在线观看66精品国产| 中国国产av一级| 成人美女网站在线观看视频| 日韩制服骚丝袜av| 69人妻影院| 日本欧美国产在线视频| 国产精品伦人一区二区| 精品少妇黑人巨大在线播放 | 五月伊人婷婷丁香| 免费看av在线观看网站| av在线观看视频网站免费| 久久久色成人| 欧美日韩一区二区视频在线观看视频在线 | 免费高清视频大片| 色综合站精品国产| 舔av片在线| 久久午夜福利片| 亚洲国产日韩欧美精品在线观看| 国产成人福利小说| 亚洲在线自拍视频| 日韩制服骚丝袜av| 别揉我奶头 嗯啊视频| 亚洲乱码一区二区免费版| 在现免费观看毛片| 国产极品精品免费视频能看的| 哪里可以看免费的av片| 久久久精品欧美日韩精品| 国产亚洲精品久久久com| 高清毛片免费看| 日本-黄色视频高清免费观看| 91在线精品国自产拍蜜月| 日韩在线高清观看一区二区三区| 成人特级av手机在线观看| 精品久久久噜噜| 菩萨蛮人人尽说江南好唐韦庄 | 国产黄a三级三级三级人| 一级毛片电影观看 | 99热6这里只有精品| 天堂av国产一区二区熟女人妻| eeuss影院久久| 午夜激情欧美在线| 欧美+亚洲+日韩+国产| 日本欧美国产在线视频| 精品一区二区三区视频在线观看免费| 亚洲av第一区精品v没综合| 欧美中文日本在线观看视频| 亚洲欧美精品综合久久99| 国内精品久久久久精免费| 国产精品一区二区三区四区免费观看 | 日日干狠狠操夜夜爽| 国产精品人妻久久久久久| 国产国拍精品亚洲av在线观看| 国产大屁股一区二区在线视频| 在线观看一区二区三区| 中文字幕免费在线视频6| 午夜福利在线在线| www日本黄色视频网| 看非洲黑人一级黄片| 精品午夜福利在线看| 免费搜索国产男女视频| 内射极品少妇av片p| 亚洲最大成人av| 国产精品野战在线观看| 国国产精品蜜臀av免费| 赤兔流量卡办理| 国产黄a三级三级三级人| 精品一区二区三区av网在线观看| 免费av毛片视频| 好男人在线观看高清免费视频| 欧美3d第一页| 99久久成人亚洲精品观看| 成人精品一区二区免费| 国产一区二区在线观看日韩| 色av中文字幕| 在线观看一区二区三区| 欧美高清性xxxxhd video| 午夜福利18| 欧美色欧美亚洲另类二区| 亚洲精品亚洲一区二区| 免费在线观看成人毛片| av天堂在线播放| 狠狠狠狠99中文字幕| 最后的刺客免费高清国语| 久久久久久久久久久丰满| 午夜精品一区二区三区免费看| 国产大屁股一区二区在线视频| 成人毛片a级毛片在线播放| 午夜精品在线福利| 不卡一级毛片| 欧美色视频一区免费| 永久网站在线| 亚洲av免费高清在线观看| 精品熟女少妇av免费看| 麻豆国产av国片精品| 日韩,欧美,国产一区二区三区 | 性色avwww在线观看| 亚洲成a人片在线一区二区| 国语自产精品视频在线第100页| 成人午夜高清在线视频| 少妇高潮的动态图| 国产精品久久视频播放| 小说图片视频综合网站| 插阴视频在线观看视频| 在线天堂最新版资源| 性色avwww在线观看| 女同久久另类99精品国产91| 国产一区二区三区在线臀色熟女| 欧美极品一区二区三区四区| 国产激情偷乱视频一区二区| 三级经典国产精品| 看非洲黑人一级黄片| 午夜日韩欧美国产| 两性午夜刺激爽爽歪歪视频在线观看| 菩萨蛮人人尽说江南好唐韦庄 | 一本久久中文字幕| 成熟少妇高潮喷水视频| 真人做人爱边吃奶动态| 小说图片视频综合网站| 亚洲内射少妇av| 干丝袜人妻中文字幕| 真人做人爱边吃奶动态| 欧美又色又爽又黄视频| 美女内射精品一级片tv| 欧美激情久久久久久爽电影| 国产一区二区三区av在线 | 免费不卡的大黄色大毛片视频在线观看 | 欧美性感艳星| 两性午夜刺激爽爽歪歪视频在线观看| 日韩av在线大香蕉| 九九在线视频观看精品| 精品一区二区三区视频在线| 亚洲av中文av极速乱| 三级毛片av免费| 又粗又爽又猛毛片免费看| 在现免费观看毛片| 亚洲天堂国产精品一区在线| 国产成人一区二区在线| 婷婷六月久久综合丁香| 国产成人aa在线观看| 午夜福利高清视频| 国产视频内射| 午夜福利18| 精品午夜福利视频在线观看一区| 午夜a级毛片| 99久久成人亚洲精品观看| 在线免费观看的www视频| 国产色爽女视频免费观看| 高清日韩中文字幕在线| 三级毛片av免费| av在线亚洲专区| 看非洲黑人一级黄片| 国产精华一区二区三区| 大又大粗又爽又黄少妇毛片口| 日本与韩国留学比较| 老司机午夜福利在线观看视频| 在线观看一区二区三区| 真实男女啪啪啪动态图| 久久久久久久午夜电影| 欧美激情国产日韩精品一区| 亚洲五月天丁香| 在线观看免费视频日本深夜| 久久午夜亚洲精品久久| 色哟哟哟哟哟哟| 一区二区三区四区激情视频 | 麻豆乱淫一区二区| 亚洲精品一卡2卡三卡4卡5卡| 欧美一区二区亚洲| 三级经典国产精品| 亚洲av免费高清在线观看| 欧美日韩一区二区视频在线观看视频在线 | 国产精品综合久久久久久久免费| 国产精品爽爽va在线观看网站| 国产精品野战在线观看| 99热只有精品国产| 亚洲精品乱码久久久v下载方式| 一级a爱片免费观看的视频| 国产精品美女特级片免费视频播放器| 身体一侧抽搐| 男插女下体视频免费在线播放| 中文亚洲av片在线观看爽| 一级毛片我不卡| 婷婷精品国产亚洲av| 99九九线精品视频在线观看视频| 午夜日韩欧美国产| 99九九线精品视频在线观看视频| 欧美日本视频| 欧美一区二区亚洲| 色5月婷婷丁香| 免费看美女性在线毛片视频| 欧美+日韩+精品| 精品久久久久久久久av| av专区在线播放| 国产av不卡久久| 欧美最黄视频在线播放免费| 麻豆国产97在线/欧美| 天堂网av新在线| 亚洲av免费在线观看| 又爽又黄a免费视频| 亚洲自拍偷在线| 99热这里只有是精品在线观看| 男女啪啪激烈高潮av片| 精品人妻熟女av久视频| 国产白丝娇喘喷水9色精品| 我要搜黄色片| 尾随美女入室| or卡值多少钱| 亚洲在线自拍视频| 国产一区二区在线观看日韩| 国产精品久久久久久久久免| 久久九九热精品免费| 日本黄色视频三级网站网址| 久久精品国产亚洲网站| 搡女人真爽免费视频火全软件 | 国产一区二区在线观看日韩| 久久久久久久亚洲中文字幕| 91久久精品国产一区二区三区| 久久鲁丝午夜福利片| 18禁裸乳无遮挡免费网站照片| 精品人妻偷拍中文字幕| 午夜福利成人在线免费观看| 国产精品国产三级国产av玫瑰| 亚洲中文日韩欧美视频| 亚洲aⅴ乱码一区二区在线播放| 免费观看在线日韩| 国产亚洲欧美98| 久久婷婷人人爽人人干人人爱| 大又大粗又爽又黄少妇毛片口| 亚洲性夜色夜夜综合| 亚洲经典国产精华液单| 可以在线观看的亚洲视频| 91av网一区二区| 丝袜喷水一区| av在线天堂中文字幕| .国产精品久久| 乱系列少妇在线播放| 精华霜和精华液先用哪个| 真实男女啪啪啪动态图| 深爱激情五月婷婷| 全区人妻精品视频| 成年免费大片在线观看| 国产精品,欧美在线| 欧美高清成人免费视频www| 男女视频在线观看网站免费| 午夜福利高清视频| 亚洲人成网站在线播| 国产国拍精品亚洲av在线观看| 国产免费男女视频| 久久久久九九精品影院| 人人妻人人澡欧美一区二区| 一本精品99久久精品77| 听说在线观看完整版免费高清| 欧美最新免费一区二区三区| 三级毛片av免费| 久久亚洲精品不卡| 99精品在免费线老司机午夜| 亚洲欧美日韩卡通动漫| 成人毛片a级毛片在线播放| 亚洲欧美精品综合久久99| 精品午夜福利在线看| 国产精品99久久久久久久久| 亚洲va在线va天堂va国产| 亚洲精品亚洲一区二区| 久久精品国产鲁丝片午夜精品| 久久久久久久午夜电影| 人人妻人人澡人人爽人人夜夜 | 亚洲一区二区三区色噜噜| 久久精品国产99精品国产亚洲性色| 丰满人妻一区二区三区视频av| 中文资源天堂在线| 国产精华一区二区三区| 2021天堂中文幕一二区在线观| 中文在线观看免费www的网站| 久久久精品欧美日韩精品| 小说图片视频综合网站| 99热这里只有是精品在线观看| 少妇人妻一区二区三区视频| 如何舔出高潮| 婷婷精品国产亚洲av在线| 少妇人妻一区二区三区视频| 免费在线观看成人毛片| 久久精品人妻少妇| 少妇熟女欧美另类| 嫩草影视91久久| 嫩草影院入口| 亚洲人成网站在线播放欧美日韩| 成人性生交大片免费视频hd| 午夜福利高清视频| 美女高潮的动态| 2021天堂中文幕一二区在线观| 99久国产av精品| 亚洲欧美日韩无卡精品| 国产日本99.免费观看| 国内少妇人妻偷人精品xxx网站| 在现免费观看毛片| 欧美区成人在线视频| 久久午夜福利片| 男人狂女人下面高潮的视频| 神马国产精品三级电影在线观看| 可以在线观看的亚洲视频| 日韩欧美 国产精品| 国产欧美日韩一区二区精品| 国产亚洲精品久久久久久毛片| 1024手机看黄色片| 国产午夜精品久久久久久一区二区三区 | 亚洲欧美日韩高清专用| 天天躁夜夜躁狠狠久久av| 成人鲁丝片一二三区免费| 国产精品av视频在线免费观看| 干丝袜人妻中文字幕| 黄片wwwwww| 国产乱人偷精品视频| 免费不卡的大黄色大毛片视频在线观看 | 美女高潮的动态| 日本与韩国留学比较| 女的被弄到高潮叫床怎么办| 亚洲熟妇中文字幕五十中出| 白带黄色成豆腐渣| 久久精品夜色国产| 丰满的人妻完整版| 欧美激情国产日韩精品一区| 亚洲国产精品成人久久小说 | 俺也久久电影网| 美女大奶头视频| 国产高清三级在线| 国产精品女同一区二区软件| 少妇高潮的动态图| 男插女下体视频免费在线播放| 亚洲中文字幕日韩| 麻豆成人午夜福利视频| 亚洲中文字幕日韩| 麻豆国产av国片精品| 一进一出抽搐gif免费好疼| 久久人人爽人人片av| 亚洲精品成人久久久久久| 免费黄网站久久成人精品| 高清日韩中文字幕在线| 午夜免费男女啪啪视频观看 | 级片在线观看| 天美传媒精品一区二区| 久久热精品热| 一本精品99久久精品77| 日本a在线网址| 在线a可以看的网站| 午夜爱爱视频在线播放| 国产精品福利在线免费观看| 亚洲欧美精品综合久久99| 九九久久精品国产亚洲av麻豆| 99久久成人亚洲精品观看| 十八禁网站免费在线| 久久久欧美国产精品| 免费观看人在逋| 黄色日韩在线| 亚洲成a人片在线一区二区| 国产av不卡久久| 女同久久另类99精品国产91| 欧美成人精品欧美一级黄| 国产一区二区亚洲精品在线观看| 高清午夜精品一区二区三区 |