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

    基于改進(jìn)遺傳算法的浮法玻璃優(yōu)化切割

    2018-08-20 03:44:00王聰聰楊明順高新勤
    關(guān)鍵詞:原片成品規(guī)格

    王聰聰,楊明順,原 欣,高新勤

    WANG Congcong,YANG Mingshun,YUAN Xin,GAO Xinqin

    西安理工大學(xué) 機(jī)械與精密儀器工程學(xué)院,西安 710048

    School of Mechanical and Precision Instrument Engineering,Xi’an University of Technology,Xi’an 710048,China

    1 引言

    浮法玻璃是平板玻璃的一種,具有表面平整光滑、透視性佳、厚度均勻等優(yōu)點(diǎn),廣泛應(yīng)用于制鏡玻璃、汽車(chē)擋風(fēng)玻璃、建筑門(mén)窗等多個(gè)方面[1]。和普通平板玻璃不同,浮法玻璃生產(chǎn)過(guò)程中的許多因素都會(huì)導(dǎo)致玻璃表面產(chǎn)生各種缺陷,而以目前的生產(chǎn)水平,浮法玻璃廠家很難杜絕缺陷的出現(xiàn)[2],這些缺陷的種類(lèi)、數(shù)量和大小會(huì)不同程度地影響玻璃質(zhì)量,進(jìn)而影響玻璃價(jià)格。因此,在浮法玻璃切割過(guò)程中,僅關(guān)注提高玻璃原板利用率的傳統(tǒng)切割方法對(duì)浮法玻璃不再適用,必須根據(jù)浮法玻璃的生產(chǎn)特性,提出新的優(yōu)化切割方案。

    玻璃切割問(wèn)題是一個(gè)典型的組合優(yōu)化問(wèn)題,類(lèi)似的問(wèn)題還有二維下料問(wèn)題、矩形件排樣問(wèn)題、二維裝箱問(wèn)題等。近年來(lái),學(xué)者們對(duì)該類(lèi)問(wèn)題進(jìn)行了大量研究。文獻(xiàn)[3]針對(duì)玻璃優(yōu)化切割問(wèn)題,以玻璃原板利用率最大為優(yōu)化目標(biāo),提出了一種啟發(fā)式優(yōu)化排樣算法,并通過(guò)實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[4]通過(guò)對(duì)玻璃切割問(wèn)題的研究,以玻璃原板利用率最大為優(yōu)化目標(biāo),提出了一種融合量子粒子群優(yōu)化和蟻群優(yōu)化的混合算法,并通過(guò)實(shí)驗(yàn)對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[5]針對(duì)有約束的兩階段二維切割問(wèn)題,以板材利用率最大為優(yōu)化目標(biāo),提出了一種并行算法,并以實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[6]針對(duì)矩形件二維下料問(wèn)題,以板材利用率最大為優(yōu)化目標(biāo),提出了一種單毛坯條帶四塊排樣方法,并通過(guò)實(shí)驗(yàn)對(duì)該方法進(jìn)行了驗(yàn)證;文獻(xiàn)[7]針對(duì)多規(guī)格的兩階段二維一刀切下料問(wèn)題,以板材利用率最大為優(yōu)化目標(biāo),提出了一種列生成啟發(fā)式算法,并以實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[8]針對(duì)一刀切矩形件排樣問(wèn)題,以板材利用率最大為優(yōu)化目標(biāo),提出了一種將啟發(fā)式遞歸與免疫克隆算法相結(jié)合的混合優(yōu)化算法,并以兩個(gè)實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[9]針對(duì)二維矩形件優(yōu)化排樣問(wèn)題,以板材利用率最大為優(yōu)化目標(biāo),提出了一種確定性啟發(fā)式最佳匹配算法,并以兩個(gè)實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[10]針對(duì)二維裝箱問(wèn)題,以填充率最大為優(yōu)化目標(biāo),提出了一種混合進(jìn)化算法,并通過(guò)實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[11]針對(duì)離線定向和不定向二維裝箱問(wèn)題,以填充率最大為優(yōu)化目標(biāo),提出了一種超啟發(fā)式算法,并通過(guò)實(shí)驗(yàn)對(duì)該算法進(jìn)行了驗(yàn)證。

    圖1 浮法玻璃優(yōu)化切割過(guò)程圖示

    上述算法都不同程度達(dá)到了減少材料浪費(fèi),提高原板利用率(填充率)的目的,但僅考慮提高玻璃原板利用率,并不能完全滿足浮法玻璃優(yōu)化切割的基本要求。浮法玻璃的價(jià)格是由切割獲得的玻璃成品的面積和玻璃質(zhì)量共同決定的,因此,即使玻璃原板利用率最大,也并不能保證玻璃成品總價(jià)格最高。

    迄今為止,僅有少數(shù)學(xué)者在解決以上組合優(yōu)化問(wèn)題時(shí),將缺陷對(duì)成品質(zhì)量的影響考慮在內(nèi)。文獻(xiàn)[12]針對(duì)考慮質(zhì)量要求的一維下料問(wèn)題,提出了一種二階段序列啟發(fā)式算法,并以實(shí)例對(duì)該算法進(jìn)行了驗(yàn)證;文獻(xiàn)[13]針對(duì)考慮缺陷的二維切割問(wèn)題,以切割獲得的矩形塊的價(jià)值最大為優(yōu)化目標(biāo),提出了一種基于動(dòng)態(tài)規(guī)劃的啟發(fā)式算法,并通過(guò)實(shí)例驗(yàn)證了算法的有效性;文獻(xiàn)[14]針對(duì)紡織工業(yè)中考慮缺陷的布料切割問(wèn)題,以總利潤(rùn)最大為優(yōu)化目標(biāo)提出了一種改進(jìn)模擬退火算法,并通過(guò)實(shí)例證明了算法的有效性;文獻(xiàn)[15]針對(duì)考慮缺陷的玻璃優(yōu)化切割問(wèn)題,以最小化切割獲得的廢品數(shù)量為優(yōu)化目標(biāo),提出了一種基于動(dòng)態(tài)規(guī)劃的算法,并通過(guò)實(shí)例證明了算法的有效性。

    本文在上述研究的基礎(chǔ)上,結(jié)合浮法玻璃的生產(chǎn)實(shí)際,考慮相應(yīng)的缺陷分級(jí)和玻璃分級(jí)標(biāo)準(zhǔn),同時(shí)考慮提高玻璃成品質(zhì)量和玻璃原板利用率,以玻璃成品總價(jià)格最大為優(yōu)化目標(biāo)建立數(shù)學(xué)模型,提出了一種符合浮法玻璃切割問(wèn)題特性的遺傳算法。

    2 浮法玻璃優(yōu)化切割問(wèn)題簡(jiǎn)介

    浮法玻璃的切割過(guò)程主要包括:缺陷檢測(cè)、生成切割方案、縱向切割和橫向切割4個(gè)步驟,如圖1所示。在實(shí)際生產(chǎn)過(guò)程中,玻璃帶被輸送輥道緩慢勻速向前輸送,在玻璃帶運(yùn)行過(guò)程中,首先,通過(guò)缺陷檢測(cè)裝置,獲得玻璃帶上各缺陷的位置、種類(lèi)和大小等信息;其次,在已知玻璃缺陷信息的基礎(chǔ)上,根據(jù)切割訂單要求,在某一固定長(zhǎng)度的玻璃帶上生成有效的切割方案;最后,按照確定的切割方案實(shí)施縱向切割和橫向切割,獲得玻璃成品。

    與普通平板玻璃切割相比,浮法玻璃優(yōu)化切割的特性主要體現(xiàn)在以下三方面:

    (1)浮法玻璃的切割工藝要求。縱切必須切整個(gè)板寬,中間不能抬刀,橫切可以改變玻璃規(guī)格。據(jù)此可知,在規(guī)劃切割方案時(shí),只有長(zhǎng)度相等的玻璃原片才能排放在玻璃原板的同一列。以表1中的切割訂單為例:訂單編號(hào)為1的玻璃原片可以和自身或編號(hào)為2的玻璃原片規(guī)劃在同一列。

    表1 切割訂單

    (2)浮法玻璃原片質(zhì)量等級(jí)的判定。按照浮法玻璃國(guó)家標(biāo)準(zhǔn)或企業(yè)自行制定的標(biāo)準(zhǔn),可將切割獲得的玻璃成品劃分為若干個(gè)質(zhì)量等級(jí),不同等級(jí)的玻璃允許含有的缺陷種類(lèi)、數(shù)量和等級(jí)不同,玻璃質(zhì)量等級(jí)越高,其每平方米價(jià)格越高,因此應(yīng)盡可能切出質(zhì)量等級(jí)較高的玻璃成品。

    (3)玻璃原片的質(zhì)量等級(jí)與其排放位置息息相關(guān)。如圖2所示,假設(shè)排放在第1列第1個(gè)位置和第2個(gè)位置的訂單編號(hào)分別為i和j(i≠j),若交換這兩塊玻璃的排放位置,并不會(huì)改變玻璃原板的利用率,但由于玻璃原片的位置變動(dòng),其上包含的缺陷信息也隨之改變。因此,當(dāng)某塊玻璃原片的位置發(fā)生變動(dòng),其質(zhì)量等級(jí)需要根據(jù)新的缺陷信息進(jìn)行重新判定。

    圖2 排放位置對(duì)玻璃原片質(zhì)量的影響

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

    浮法玻璃優(yōu)化切割問(wèn)題是指:在一塊已知長(zhǎng)度、寬度以及缺陷信息的浮法玻璃板面上,根據(jù)訂單要求切割一定規(guī)格的玻璃原片時(shí),應(yīng)優(yōu)先規(guī)劃質(zhì)量較好的玻璃原片,并保證玻璃原板的利用率較大,以使得切割所得玻璃成品的銷(xiāo)售總價(jià)最高。

    該問(wèn)題中各變量及其符號(hào)表示如下:

    Length為玻璃原板的長(zhǎng)度;Width為玻璃原板的寬度;n為玻璃原板上的缺陷點(diǎn)總數(shù);D為缺陷集,其中即每個(gè)缺陷點(diǎn)均含有3個(gè)基本屬性:缺陷位置pi,缺陷種類(lèi)ti,缺陷大小si,其中pi=(xi,yi),xi、yi分別為缺陷點(diǎn)中心到切割原點(diǎn)的橫向和縱向距離,ti∈{氣泡,夾雜,光畸變,錫灰,玻筋,線道,其他}。

    m為切割訂單包含的訂單規(guī)格總數(shù);CO表示切割訂單,其中Oj=(lj,wj),lj和wj分別表示切割訂單中各訂單規(guī)格對(duì)應(yīng)的長(zhǎng)度和寬度,l和w分別表示切割訂單中所有訂單規(guī)格對(duì)應(yīng)的長(zhǎng)度和寬度的集合,

    DC為缺陷分級(jí)參數(shù),用于判定缺陷的等級(jí),DC=根據(jù)缺陷分級(jí)參數(shù)可將缺陷劃分為H+1個(gè)等級(jí),等級(jí)為H+1的缺陷也稱(chēng)為“不允許缺陷”。

    QC為玻璃分級(jí)參數(shù),用于判定玻璃原片的質(zhì)量等級(jí),其中O≤M1≤M2≤…≤MT,且矩陣Mt互不相等(“≤”表示前一個(gè)矩陣的所有元素均不超過(guò)后一個(gè)矩陣對(duì)應(yīng)位置的元素)。Mt表示一個(gè)N×(H+1)的矩陣,O表示一個(gè)和矩陣Mt同型的零矩陣,N為玻璃原板上包含的缺陷種類(lèi)數(shù),矩陣Mt中第i行第j列的元素Mt(i,j)表示等級(jí)為t的玻璃每平方米允許含有種類(lèi)為i等級(jí)為j的缺陷的數(shù)量,根據(jù)玻璃分級(jí)參數(shù)可將玻璃劃分為T(mén)+1個(gè)質(zhì)量等級(jí),等級(jí)為T(mén)+1的玻璃也稱(chēng)為“廢品級(jí)玻璃”。

    SP為玻璃單位面積價(jià)格的集合,SP={St|t=1,其中St表示等級(jí)為t的玻璃每平方米價(jià)格為St。

    xc,r為決策變量,其值等于排放在玻璃原板第c列的第r塊玻璃原片的訂單規(guī)格編號(hào);colMax為玻璃原板上排放玻璃原片的最大列數(shù);glaNumMax為每列最多可以排放的玻璃原片塊數(shù);colReal為玻璃原板上實(shí)際排放的玻璃原片的列數(shù);glaNum(c)為第c列實(shí)際排放的玻璃原片的塊數(shù);集合cset為列數(shù)集合cset={1,2,…,colMax};rset為塊數(shù)集合,rset={1,2,…,glaNum(c)};lrnd等于集合l中的任意一個(gè)值,即glassArr表示一個(gè)可行的切割方案;g(c,r)表示該切割方案中排放在玻璃原板第c列的第r塊玻璃原片的信息集:

    其中,xst(xc,r)和yst(xc,r)分別為該玻璃原片左下角距離切割原點(diǎn)的橫向和縱向距離;lr(xc,r)和wr(xc,r)分別表示該玻璃原片的長(zhǎng)度和寬度;Dr(xc,r)為該玻璃原片上包含的缺陷集,Dr(xc,r)∈D,nr為該玻璃原片上的缺陷總數(shù),dpk、dtk、dsk分別表示該玻璃原片上各缺陷的位置、種類(lèi)和大小,為該玻璃原片對(duì)應(yīng)的缺陷計(jì)數(shù)矩陣,用來(lái)存儲(chǔ)該玻璃原片實(shí)際每平方米包含的不同種類(lèi)、不同等級(jí)缺陷的數(shù)量;glag(xc,r)為該玻璃原片的質(zhì)量等級(jí);pr(xc,r)為該塊玻璃原片每平方米的價(jià)格。圖3為一個(gè)可行切割方案示例。

    該問(wèn)題的目標(biāo)函數(shù)及約束條件如下:

    圖3 可行切割方案示例

    式(1)為目標(biāo)函數(shù)的計(jì)算公式,其值為玻璃成品的總價(jià)格;式(2)、(3)為玻璃原板最多可以排放玻璃原片的列數(shù)colMax,及每列最多可以排放的玻璃原片的塊數(shù)glaNumMax的計(jì)算公式;式(4)、(5)表示玻璃原板實(shí)際排放的玻璃原片的列數(shù)不超過(guò)colMax,每列實(shí)際排放的玻璃原片的塊數(shù)不超過(guò)glaNumMax;式(6)、(7)分別為玻璃原板的長(zhǎng)度約束和寬度約束,表示排放的所有玻璃原片均不能超過(guò)玻璃原板的邊界;式(8)為排放在同一列的玻璃原片的等長(zhǎng)約束,即排放在玻璃原板同一列的玻璃原片的長(zhǎng)必須相等;式(9)為決策變量的取值范圍,表示排放的玻璃原片的規(guī)格必須來(lái)自切割訂單;式(10)、(11)為該玻璃原片上缺陷的位置約束,表示同時(shí)滿足式(10)、(11)約束的缺陷點(diǎn)在該玻璃原片上;式(12)為缺陷等級(jí)的判定公式;式(13)為玻璃質(zhì)量等級(jí)的判定公式,其含義為已知某塊玻璃的缺陷計(jì)數(shù)矩陣CM,玻璃分級(jí)參數(shù)QC,首先判斷CM≤M1是否成立(CM的所有元素均不超過(guò)玻璃分級(jí)參數(shù)矩陣M1對(duì)應(yīng)位置的元素),若成立,則該玻璃等級(jí)為1,否則降級(jí),繼續(xù)判斷CM≤M2是否成立,若成立,則該玻璃的等級(jí)為2,否則降級(jí),繼續(xù)判斷,以此類(lèi)推,若直至CM≤MT仍不成立,則該玻璃等級(jí)為T(mén)+1(即廢品級(jí));式(14)為玻璃原片單位面積價(jià)格的判定公式,用于根據(jù)玻璃原片的質(zhì)量等級(jí)進(jìn)一步判定其每平方米的價(jià)格。

    4 算法設(shè)計(jì)

    浮法玻璃優(yōu)化切割屬于NP難類(lèi)組合優(yōu)化問(wèn)題,隨著訂單規(guī)格的增加及缺陷信息的多樣化,尋優(yōu)過(guò)程將變得十分復(fù)雜,此時(shí)很難在有限時(shí)間內(nèi)遍歷所有解空間,從中找出最優(yōu)解,因此精確求解算法不再適用。遺傳算法是基于進(jìn)化理論原理發(fā)展起來(lái)的一種高效的隨機(jī)搜索算法,在解決組合優(yōu)化問(wèn)題時(shí)具有較好的性能。本文根據(jù)浮法玻璃優(yōu)化切割問(wèn)題的特點(diǎn),提出了一種基于玻璃原片排放位置編碼的遺傳算法。以下對(duì)該算法的基本步驟進(jìn)行說(shuō)明。

    4.1 編碼

    根據(jù)浮法玻璃切割“縱切必須切整個(gè)板寬”的要求(即排在同一列的訂單規(guī)格的長(zhǎng)度必須相等),以及玻璃原片排放位置對(duì)其質(zhì)量等級(jí)的影響,提出了一種基于玻璃原片排放位置的序列編碼方式。每條染色體含有g(shù)n=colMax×glaNumMax個(gè)基因,每個(gè)基因gi的取值為切割訂單中某一個(gè)訂單規(guī)格的編號(hào),即區(qū)間 [glaNumMax×(k-1)+1,k×glaNumMax]上的基因值對(duì)應(yīng)的玻璃規(guī)格的長(zhǎng)度相等,

    以表1中的訂單規(guī)格為例,假設(shè)玻璃原板的長(zhǎng)Length為12200,寬Width為4000,根據(jù)式(2)、(3)可知該玻璃原板最多可以排放6列,每列最多可以排放2塊玻璃原片,染色體基因數(shù)gn=6×2=12。經(jīng)過(guò)編碼得到一個(gè)可能的染色體如圖4所示,該染色體表示玻璃原板上每個(gè)位置擬排放的訂單規(guī)格。如:第一個(gè)基因位上的元素7表示排放在玻璃原板第一列的第一塊玻璃原片的訂單規(guī)格編號(hào)為7;而編號(hào)為7和編號(hào)為6的訂單規(guī)格的長(zhǎng)度相等,可以排在同一列。

    圖4 染色體編碼示例

    4.2 解碼

    為得到具體規(guī)劃方案,需要對(duì)染色體進(jìn)行解碼。根據(jù)玻璃原板的寬度約束和長(zhǎng)度約束從左向右對(duì)染色體依次進(jìn)行解碼,獲得每列實(shí)際允許排放的玻璃規(guī)格。圖4所示染色體的解碼結(jié)果為glassArr={(7,6),(5),(1,2),(4,3),(2,1)},其含義為:玻璃原板第一列實(shí)際排放的訂單規(guī)格編號(hào)依次為7和6,第二列排放的訂單規(guī)格編號(hào)為5,以此類(lèi)推。

    4.3 種群初始化

    為保證初始種群的多樣性,采用一種啟發(fā)式搜索算法獲得符合浮法玻璃切割要求的初始種群。一個(gè)可行初始解的生成過(guò)程如下:

    步驟1輸入玻璃原板上可排放玻璃原片的最大列數(shù)colMax,每列可排放的玻璃原片的最大塊數(shù)glaNumMax。

    步驟2令c=1,r=2。

    步驟3從切割訂單中隨機(jī)選擇一個(gè)訂單規(guī)格編號(hào)firstp,排放在第c列的第1個(gè)位置,即染色體的第glaNumMax×(c-1)+1個(gè)位置。

    步驟4在切割訂單中尋找和firstp等長(zhǎng)的所有玻璃規(guī)格的編號(hào),存放在集合eql中。

    步驟5在等長(zhǎng)訂單規(guī)格集合eql中隨機(jī)選擇一個(gè)訂單規(guī)格編號(hào),排放在第c列第r個(gè)位置,即染色體的第glaNumMax×(c-1)+r個(gè)位置。

    步驟6判斷r=glaNumMax是否成立,若成立,繼續(xù)進(jìn)行步驟7,否則令r=r+1,轉(zhuǎn)步驟5。

    步驟7判斷c=colMax是否成立,若成立,結(jié)束程序,否則,令r=2,c=c+1,轉(zhuǎn)步驟3。

    將以上過(guò)程循環(huán)N(p)次,即可獲得該問(wèn)題的初始種群。

    4.4 計(jì)算適應(yīng)度函數(shù)值

    采用該問(wèn)題的目標(biāo)函數(shù)值作為該算法的適應(yīng)度函數(shù)值,第i個(gè)染色體的適應(yīng)度函數(shù)值的計(jì)算公式為:

    某塊玻璃原片的價(jià)格是由其面積lr×wr及其單位面積的價(jià)格pr共同決定的,當(dāng)訂單規(guī)格已知時(shí),該玻璃原片的面積是確定的,但該玻璃原片單位面積的價(jià)格需要由其包含的缺陷信息判定。

    以下給出根據(jù)某塊玻璃原片上包含的缺陷信息,判定其質(zhì)量等級(jí),進(jìn)而判定其單位面積價(jià)格的步驟:

    步驟1輸入缺陷分級(jí)參數(shù)DC,玻璃分級(jí)參數(shù)QC,玻璃單位面積價(jià)格的集合SP,玻璃原片上包含的缺陷集合Dr。

    步驟2根據(jù)缺陷分級(jí)參數(shù)DC判定該玻璃原片上各缺陷的等級(jí),見(jiàn)式(12),整理缺陷的種類(lèi)和等級(jí)信息,獲得該玻璃原片的缺陷計(jì)數(shù)矩陣CM。

    步驟3根據(jù)缺陷計(jì)數(shù)矩陣CM判定該玻璃原片等級(jí),見(jiàn)式(13)。

    步驟4根據(jù)玻璃原片等級(jí)判定其每平方米價(jià)格,見(jiàn)式(14)。

    根據(jù)以上過(guò)程,依次判斷glassArr中每塊玻璃原片的單位面積價(jià)格,再根據(jù)式(15)計(jì)算該染色體對(duì)應(yīng)的適應(yīng)度函數(shù)值。

    4.5 選擇算子

    采用隨機(jī)按比例的方式選擇個(gè)體進(jìn)入交配庫(kù),使得適應(yīng)度函數(shù)值較大的個(gè)體被選擇概率較大。

    步驟1將所有個(gè)體的適應(yīng)度函數(shù)值按照升序排列,得到集合Fsort。

    步驟3將步驟2循環(huán)N(p)次,直至選滿N(p)個(gè)個(gè)體為止。

    4.6 交叉算子

    為保證交叉后染色體的可行性,本文針對(duì)浮法玻璃切割問(wèn)題的特點(diǎn),提出了一種兩點(diǎn)交叉法。

    步驟1從選擇操作產(chǎn)生的交配庫(kù)中隨機(jī)選取兩個(gè)染色體F1、F2,作為交叉操作的父代。

    步驟2隨機(jī)產(chǎn)生兩個(gè)交叉點(diǎn)P1、P2,為保證兩個(gè)父代經(jīng)過(guò)交叉后不破壞染色體的可行性,P1,P2∈{k×glaNumMax},其中k=1,2,…,colMax-1,且P1≠P2。

    步驟3將父代染色體F1中第[1,P1]、[P2+1,gn]之間的基因復(fù)制到子代染色體O1的對(duì)應(yīng)位置,其中g(shù)n為染色體包含的基因總數(shù)。

    步驟4將父代染色體F2中第[P1+1,P2]之間的基因復(fù)制到染色體O1的對(duì)應(yīng)位置,獲得子代染色體O1。

    步驟5將F1和F2的位置交換,同理,可以得到后代O2。

    如圖5所示,對(duì)父代F1、F2執(zhí)行交叉操作,選擇交叉點(diǎn)P1=4,P2=8,按照以上方式執(zhí)行交叉操作,得到子代O1、O2。

    圖5 染色體交叉示例

    4.7 變異算子

    為避免搜索過(guò)程過(guò)早地陷入局部最優(yōu),變異算子隨機(jī)選擇部分染色體的部分基因進(jìn)行變異,具體步驟如下:

    步驟1依據(jù)變異概率Pm和變異臨時(shí)隨機(jī)數(shù)Pmt的大小關(guān)系判斷是否對(duì)該個(gè)體進(jìn)行變異,若Pmt≤Pm(Pmt,Pm∈[0,1])成立,則對(duì)該個(gè)體進(jìn)行變異操作,否則,不進(jìn)行變異操作。

    步驟2隨機(jī)選擇排放在玻璃原板的第k(k=1,2,…,colMax)列的所有玻璃原片進(jìn)行變異,即選擇染色體在區(qū)間[glaNumMax×(k-1)+1,k×glaNumMax]之間的基因進(jìn)行變異。

    步驟3隨機(jī)生成一組可行的列訂單組合替換排放在第k列的訂單規(guī)格,得到變異后的染色體。一個(gè)可行列訂單組合的生成步驟為:從切割訂單中隨機(jī)選擇一個(gè)玻璃規(guī)格編號(hào)firstp,排放在第1個(gè)位置;在切割訂單中尋找和firstp等長(zhǎng)的所有玻璃規(guī)格的編號(hào),存放在集合eql中;在等長(zhǎng)訂單規(guī)格集合eql中隨機(jī)選擇訂單規(guī)格編號(hào),排放在第[2,glaNumMax]個(gè)位置。

    為避免種群收縮于局部最優(yōu)解,又使種群進(jìn)化保持一定的連續(xù)性,在進(jìn)化過(guò)程中采用動(dòng)態(tài)變異率。在進(jìn)化的初級(jí)階段,采用較高的變異率,使種群具有足夠的突變性,避免收斂到局部最優(yōu)解,而隨著進(jìn)化代數(shù)的提高,降低種群的變異率,避免種群無(wú)序化,使其在較優(yōu)情況下向著最優(yōu)解靠攏。取第t代個(gè)體的變異率Pm,t=N/(N+t)×Pm,其中Pm為最初設(shè)定的變異率,N為開(kāi)始階段設(shè)定的最大迭代次數(shù),t為當(dāng)前迭代次數(shù)。

    4.8 精英保留策略

    為加快算法的收斂速度,防止進(jìn)化結(jié)果停滯不前或發(fā)生隨機(jī)漫游現(xiàn)象,采用最優(yōu)保留策略,即從本代染色體中選擇適應(yīng)度值最高的染色體中的一個(gè),來(lái)取代子代染色體中適應(yīng)度值最低的染色體中的一個(gè)。

    4.9 算法整體流程

    設(shè)P(t)為第t代種群,算法整體流程如下:

    步驟1輸入相關(guān)參數(shù)。輸入浮法玻璃切割相關(guān)參數(shù),即種群規(guī)模N(p),交叉概率Pc,變異概率Pm,最大迭代次數(shù)N。

    步驟2初始化種群。令t=0,產(chǎn)生個(gè)體數(shù)為N(p)的初始種群P(0)。

    步驟3計(jì)算適應(yīng)度函數(shù)值。根據(jù)缺陷分級(jí)參數(shù)、玻璃分級(jí)參數(shù)判定每塊玻璃原片的質(zhì)量等級(jí),判定其單位面積的價(jià)格,進(jìn)而計(jì)算第t代種群P(t)中每個(gè)個(gè)體的適應(yīng)度函數(shù)值。

    步驟4選擇操作。采用隨機(jī)按比例的方式,從P(t)中選擇N(p)個(gè)個(gè)體復(fù)制到P(t+1)中。

    步驟5交叉操作。依據(jù)交叉臨時(shí)隨機(jī)數(shù)Pct和交叉概率Pc的大小關(guān)系,Pct、Pc∈[0,1],從選擇操作所得種群中隨機(jī)選擇一對(duì)臨時(shí)父代進(jìn)行交叉,共進(jìn)行N(p)/2次,產(chǎn)生N(p)個(gè)后代形成的臨時(shí)種群。

    步驟6變異操作。對(duì)于交叉操作產(chǎn)生的臨時(shí)種群,依據(jù)變異臨時(shí)隨機(jī)數(shù)Pmt和當(dāng)代變異概率Pm,t的大小關(guān)系,對(duì)個(gè)體進(jìn)行變異操作,產(chǎn)生N(p)個(gè)后代形成新的后代種群。

    步驟7最優(yōu)保留策略。用P(t)中適應(yīng)度函數(shù)值最大的一個(gè)個(gè)體取代P(t+1)中適應(yīng)度函數(shù)值最小的一個(gè)個(gè)體。

    步驟8令t=t+1,判斷是否達(dá)到最大迭代次數(shù),即若t=N,輸出計(jì)算結(jié)果,結(jié)束程序;否則,轉(zhuǎn)步驟3。

    5 實(shí)例分析

    5.1 實(shí)例說(shuō)明

    現(xiàn)有一塊長(zhǎng)12200 mm、寬4000 mm的浮法玻璃原板,經(jīng)檢測(cè),該玻璃板面上共含有40個(gè)缺陷點(diǎn),表2為各缺陷點(diǎn)的詳細(xì)信息。已知缺陷等級(jí)參數(shù)DC={0.5,1.0,1.5,3,5},玻璃分級(jí)參數(shù)QC={M1,M2,M3},各等級(jí)玻璃單位面積的價(jià)格依次為11.2 元/m2,10.5 元/m2,9.6 元/m2,0元/m2,在該玻璃原板上切割表1所示訂單規(guī)格,要求盡可能使獲得的玻璃成品總價(jià)格最高。

    表2 缺陷信息

    表3 實(shí)驗(yàn)結(jié)果統(tǒng)計(jì)

    玻璃分級(jí)參數(shù)QC中M1、M2、M3的取值如下:

    5.2 求解結(jié)果與分析

    采用本文提出的改進(jìn)遺傳算法解決上述算例。令種群規(guī)模N(p)=10,迭代次數(shù)N=50,染色體交叉概率Pc=0.8,變異概率Pm=0.2。對(duì)遺傳算法而言,其優(yōu)化結(jié)果并不是唯一、確定的,因此針對(duì)以上浮法玻璃切割問(wèn)題,將該遺傳算法運(yùn)行10次,計(jì)算結(jié)果見(jiàn)表3。

    根據(jù)表3所示實(shí)驗(yàn)結(jié)果,可得出以下結(jié)論:

    (1)對(duì)浮法玻璃優(yōu)化切割問(wèn)題而言,切割獲得的玻璃成品塊數(shù)越多,玻璃成品總價(jià)不一定最大。如:第1次和第2次實(shí)驗(yàn)均獲得12塊玻璃成品,但其玻璃成品的總價(jià)格均不是最大。

    (2)切割獲得的玻璃成品總面積最大,玻璃成品總價(jià)不一定最大。如:第6次實(shí)驗(yàn)獲得的玻璃成品總面積最大,但其玻璃成品總價(jià)并不是最大。

    (3)玻璃成品的總價(jià)是由玻璃成品面積和玻璃質(zhì)量等級(jí)共同決定的。如:第7次實(shí)驗(yàn)獲得10塊玻璃成品,玻璃成品總面積為39.7788 m2,但其所得玻璃成品總價(jià)最大,是上述實(shí)驗(yàn)中的最佳方案。

    以上10次實(shí)驗(yàn)中,第7次實(shí)驗(yàn)獲得的玻璃成品總價(jià)最大,將該次實(shí)驗(yàn)對(duì)應(yīng)的切割方案作為上述實(shí)例的最終切割方案,圖6為最終切割方案的示意圖。

    圖6 玻璃原板缺陷信息及切割方案圖示

    5.3 算法性能比較

    針對(duì)浮法玻璃優(yōu)化切割問(wèn)題,以玻璃成品總價(jià)格最大為優(yōu)化目標(biāo),本文還提出了一種基于局部遍歷的啟發(fā)式搜索算法,該算法在保證排放在每個(gè)位置上的玻璃原片的質(zhì)量等級(jí)最高的基礎(chǔ)上,使其面積盡可能大,以保證排放每個(gè)位置的玻璃原片的價(jià)格最大。針對(duì)表4中6組不同的切割訂單,將其優(yōu)化結(jié)果與改進(jìn)遺傳算法進(jìn)行比較。

    由表5可知,和隨機(jī)啟發(fā)式算法相比,本文提出的改進(jìn)遺傳算法具有更好的性能,且針對(duì)一個(gè)確定切割訂單和玻璃原板,隨機(jī)啟發(fā)式算法只能找到一組解決方案,而遺傳算法每次運(yùn)算可以獲得最優(yōu)解對(duì)應(yīng)的多個(gè)解決方案,可以為規(guī)劃人員提供多種選擇。

    表4 切割訂單

    表5 算法優(yōu)化性能對(duì)比

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

    本文針對(duì)浮法玻璃切割問(wèn)題,根據(jù)缺陷分級(jí)參數(shù)和玻璃分級(jí)參數(shù)對(duì)玻璃原片的質(zhì)量等級(jí)進(jìn)行判定,進(jìn)而獲得玻璃原片的價(jià)格,以玻璃成品總價(jià)最大為優(yōu)化目標(biāo),建立該問(wèn)題的數(shù)學(xué)模型,提出了一種基于玻璃原片排放位置編碼的遺傳算法。對(duì)算法的各個(gè)模塊進(jìn)行了設(shè)計(jì),并以實(shí)例驗(yàn)證了所建模型的正確性及算法的有效性。實(shí)驗(yàn)結(jié)果表明:浮法玻璃的成品總價(jià)格是由玻璃成品面積和玻璃質(zhì)量等級(jí)共同決定的,僅追求玻璃利用率最高或玻璃質(zhì)量等級(jí)最高,均不能滿足生產(chǎn)利益最大化的要求。本文提出的遺傳算法,在一定程度上能夠提高玻璃成品的銷(xiāo)售價(jià)格,具有較好的性能。

    猜你喜歡
    原片成品規(guī)格
    轉(zhuǎn)換性使用在短視頻解說(shuō)著作權(quán)中的應(yīng)用研究
    法制博覽(2023年35期)2023-02-13 14:07:46
    近3成苗企難以維持!規(guī)格越大越虧,2022如何讓泥鰍賺錢(qián)?
    閉月羞花
    寶藏(2021年6期)2021-07-20 06:12:22
    一種太陽(yáng)能光伏玻璃深加工連線工藝探討
    建材世界(2020年5期)2020-11-05 04:57:04
    汽車(chē)前風(fēng)窗玻璃光畸變問(wèn)題影響因素及改善方案
    山鋼日鋼1.2mm極限薄規(guī)格帶鋼的生產(chǎn)實(shí)踐
    山東冶金(2019年1期)2019-03-30 01:35:30
    2017年1—4月熱帶作物及其成品進(jìn)出口情況
    2017年1—3月熱帶作物及其成品進(jìn)出口情況
    2017年2月熱帶作物及其成品進(jìn)出口情況(續(xù))
    2017年1—2月熱帶作物及其成品進(jìn)出口情況(續(xù))
    啦啦啦在线免费观看视频4| 男女无遮挡免费网站观看| 岛国毛片在线播放| 香蕉国产在线看| 欧美少妇被猛烈插入视频| 天天躁日日躁夜夜躁夜夜| 午夜视频精品福利| 亚洲成av片中文字幕在线观看| cao死你这个sao货| 午夜成年电影在线免费观看| 美女国产高潮福利片在线看| 亚洲精品久久成人aⅴ小说| 天堂中文最新版在线下载| 99久久国产精品久久久| 亚洲国产欧美网| 国产欧美日韩一区二区三区在线| 女人久久www免费人成看片| 日本av免费视频播放| 国产人伦9x9x在线观看| 国产亚洲一区二区精品| 国产一区二区三区综合在线观看| 欧美日韩福利视频一区二区| 中文字幕高清在线视频| 好男人电影高清在线观看| 欧美老熟妇乱子伦牲交| 国产一区二区三区在线臀色熟女 | 丁香六月天网| 99精品欧美一区二区三区四区| 男女免费视频国产| 黄片大片在线免费观看| 免费久久久久久久精品成人欧美视频| 丁香六月欧美| 久久国产精品影院| 在线观看人妻少妇| 国产伦理片在线播放av一区| 黄色视频不卡| 高清黄色对白视频在线免费看| 免费高清在线观看视频在线观看| 中文字幕人妻熟女乱码| av片东京热男人的天堂| 搡老乐熟女国产| 日日夜夜操网爽| 岛国在线观看网站| 婷婷色av中文字幕| av天堂在线播放| 免费少妇av软件| 亚洲专区中文字幕在线| 曰老女人黄片| 婷婷色av中文字幕| 91成人精品电影| 国产伦理片在线播放av一区| 老司机影院毛片| 亚洲成人免费电影在线观看| 一本久久精品| 久久久久网色| 在线观看www视频免费| 丝瓜视频免费看黄片| 日韩 欧美 亚洲 中文字幕| 91麻豆精品激情在线观看国产 | 高清黄色对白视频在线免费看| 麻豆国产av国片精品| 中文字幕色久视频| 亚洲av电影在线观看一区二区三区| 热re99久久国产66热| 日韩 亚洲 欧美在线| 国产欧美日韩一区二区三区在线| 国产精品av久久久久免费| 日本撒尿小便嘘嘘汇集6| 99九九在线精品视频| 国产在线一区二区三区精| 最近中文字幕2019免费版| 亚洲av男天堂| 9热在线视频观看99| 成年动漫av网址| 久久久久久久久久久久大奶| 久久久久久久久免费视频了| 99久久人妻综合| 黑丝袜美女国产一区| 日本vs欧美在线观看视频| 精品国产国语对白av| 久久精品国产a三级三级三级| 免费高清在线观看日韩| 亚洲国产中文字幕在线视频| 国产亚洲av高清不卡| 亚洲国产精品999| 国产色视频综合| 日本av手机在线免费观看| 成人av一区二区三区在线看 | 一本大道久久a久久精品| 日本欧美视频一区| 十八禁人妻一区二区| 国产一区二区 视频在线| 蜜桃国产av成人99| 亚洲性夜色夜夜综合| 人人妻人人澡人人看| 黑人操中国人逼视频| 国产日韩欧美亚洲二区| 国产成人精品无人区| 午夜成年电影在线免费观看| 女人被躁到高潮嗷嗷叫费观| 国产在线免费精品| 亚洲综合色网址| 精品久久久精品久久久| 国产精品久久久人人做人人爽| 久久天躁狠狠躁夜夜2o2o| 久久影院123| 建设人人有责人人尽责人人享有的| 午夜视频精品福利| av有码第一页| 久久久久久亚洲精品国产蜜桃av| 狂野欧美激情性bbbbbb| 国产三级黄色录像| 久久久国产一区二区| 久久久久久亚洲精品国产蜜桃av| 国产精品一区二区精品视频观看| 在线天堂中文资源库| 久久久精品国产亚洲av高清涩受| 久久久久网色| 欧美日本中文国产一区发布| 精品一区二区三区四区五区乱码| 青青草视频在线视频观看| 三上悠亚av全集在线观看| 国产精品一区二区在线观看99| 日韩大片免费观看网站| 国产精品一二三区在线看| 成年人免费黄色播放视频| 一区二区av电影网| av电影中文网址| 欧美黑人精品巨大| 欧美日韩福利视频一区二区| 美女国产高潮福利片在线看| 黄色视频不卡| 免费久久久久久久精品成人欧美视频| 美女视频免费永久观看网站| 亚洲欧美一区二区三区黑人| √禁漫天堂资源中文www| 亚洲成人免费电影在线观看| 丝袜美腿诱惑在线| 多毛熟女@视频| 啦啦啦啦在线视频资源| 亚洲七黄色美女视频| 成人影院久久| 新久久久久国产一级毛片| 最近最新中文字幕大全免费视频| 亚洲色图综合在线观看| 精品人妻1区二区| 亚洲七黄色美女视频| 正在播放国产对白刺激| 涩涩av久久男人的天堂| 亚洲五月婷婷丁香| 精品高清国产在线一区| 欧美午夜高清在线| 国产成人免费观看mmmm| 国产精品1区2区在线观看. | 精品亚洲乱码少妇综合久久| 视频在线观看一区二区三区| 亚洲欧美一区二区三区久久| 亚洲成人国产一区在线观看| 成年女人毛片免费观看观看9 | 精品第一国产精品| tube8黄色片| 日韩人妻精品一区2区三区| 50天的宝宝边吃奶边哭怎么回事| 欧美97在线视频| 在线观看免费日韩欧美大片| 中文字幕精品免费在线观看视频| 国产在线一区二区三区精| 日本黄色日本黄色录像| 性色av乱码一区二区三区2| 麻豆av在线久日| 精品久久蜜臀av无| 亚洲欧美一区二区三区久久| 国产精品1区2区在线观看. | 亚洲成人免费av在线播放| 亚洲av日韩精品久久久久久密| 后天国语完整版免费观看| 国产高清videossex| 男人爽女人下面视频在线观看| av在线老鸭窝| 国产有黄有色有爽视频| 国产麻豆69| 国产免费视频播放在线视频| 亚洲全国av大片| 女人久久www免费人成看片| 精品欧美一区二区三区在线| 久久久久久久大尺度免费视频| 香蕉丝袜av| 欧美精品av麻豆av| 欧美av亚洲av综合av国产av| 午夜精品久久久久久毛片777| 午夜福利影视在线免费观看| 精品少妇黑人巨大在线播放| 久久久久久亚洲精品国产蜜桃av| 亚洲国产欧美日韩在线播放| 国产伦人伦偷精品视频| 国精品久久久久久国模美| 天天影视国产精品| 一进一出抽搐动态| 动漫黄色视频在线观看| 满18在线观看网站| 在线观看舔阴道视频| 在线永久观看黄色视频| 精品一区二区三区四区五区乱码| 天堂俺去俺来也www色官网| 性高湖久久久久久久久免费观看| 中文字幕最新亚洲高清| 搡老岳熟女国产| 亚洲熟女精品中文字幕| 九色亚洲精品在线播放| 精品国产一区二区三区久久久樱花| 国产成人啪精品午夜网站| 久久久久久人人人人人| 男女边摸边吃奶| 99国产精品99久久久久| 久久精品aⅴ一区二区三区四区| 日韩有码中文字幕| 国产精品久久久久久精品电影小说| 十八禁网站网址无遮挡| 欧美精品av麻豆av| 日韩电影二区| 精品国产一区二区久久| 久久国产精品影院| 最近最新中文字幕大全免费视频| netflix在线观看网站| 日韩欧美一区二区三区在线观看 | 美女脱内裤让男人舔精品视频| 久久影院123| 少妇精品久久久久久久| 视频区图区小说| 97精品久久久久久久久久精品| 久久精品人人爽人人爽视色| 久久人人爽av亚洲精品天堂| 精品人妻在线不人妻| 国产成人av教育| 黑人巨大精品欧美一区二区mp4| 大陆偷拍与自拍| 黑人猛操日本美女一级片| 日本五十路高清| 桃红色精品国产亚洲av| 国产精品99久久99久久久不卡| av有码第一页| 一级毛片女人18水好多| a级片在线免费高清观看视频| 国产成人免费观看mmmm| 老司机靠b影院| 啪啪无遮挡十八禁网站| 亚洲国产欧美网| 免费在线观看完整版高清| 国产成人欧美在线观看 | 日韩 亚洲 欧美在线| 国产一区二区三区在线臀色熟女 | 亚洲九九香蕉| 久久99热这里只频精品6学生| 青春草视频在线免费观看| 日本一区二区免费在线视频| 丝袜在线中文字幕| 在线观看免费日韩欧美大片| 老司机亚洲免费影院| 岛国毛片在线播放| xxxhd国产人妻xxx| 久久av网站| 久久久久国产一级毛片高清牌| 中文字幕制服av| 亚洲精品日韩在线中文字幕| 99国产精品一区二区三区| 国产在线免费精品| 国产高清国产精品国产三级| 精品人妻在线不人妻| 十分钟在线观看高清视频www| 蜜桃国产av成人99| 两性夫妻黄色片| www日本在线高清视频| 90打野战视频偷拍视频| 99国产精品99久久久久| 国产精品九九99| 久久午夜综合久久蜜桃| 久久久欧美国产精品| 国产成人欧美在线观看 | 999精品在线视频| 欧美在线一区亚洲| 热99国产精品久久久久久7| 国产成人精品在线电影| 亚洲欧美日韩高清在线视频 | 欧美亚洲 丝袜 人妻 在线| 超碰成人久久| 成在线人永久免费视频| 王馨瑶露胸无遮挡在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 99国产综合亚洲精品| 脱女人内裤的视频| 高清视频免费观看一区二区| 国产精品一区二区在线不卡| 日韩视频在线欧美| 可以免费在线观看a视频的电影网站| 久久久久久久国产电影| 国产欧美日韩一区二区精品| 另类精品久久| 久久精品亚洲av国产电影网| 国产日韩欧美视频二区| 最新的欧美精品一区二区| av不卡在线播放| 国产欧美日韩精品亚洲av| 伦理电影免费视频| 男女无遮挡免费网站观看| 久久久久久亚洲精品国产蜜桃av| 国产深夜福利视频在线观看| 女警被强在线播放| 好男人电影高清在线观看| 制服诱惑二区| 另类精品久久| 亚洲熟女精品中文字幕| 老汉色∧v一级毛片| 超色免费av| 国产精品免费视频内射| 亚洲国产欧美一区二区综合| 老司机影院成人| 五月开心婷婷网| 欧美一级毛片孕妇| 动漫黄色视频在线观看| 视频区欧美日本亚洲| 桃花免费在线播放| 国产精品久久久久久精品电影小说| 一本一本久久a久久精品综合妖精| 在线观看免费午夜福利视频| 丰满少妇做爰视频| 午夜福利乱码中文字幕| 69精品国产乱码久久久| 精品久久久精品久久久| 少妇被粗大的猛进出69影院| 国产av又大| 欧美人与性动交α欧美软件| 国产91精品成人一区二区三区 | 成人三级做爰电影| 蜜桃国产av成人99| 国产伦人伦偷精品视频| 制服人妻中文乱码| 人人妻,人人澡人人爽秒播| 久久久久国内视频| 99久久综合免费| 亚洲精品美女久久久久99蜜臀| 男女高潮啪啪啪动态图| 91字幕亚洲| 精品福利永久在线观看| 男男h啪啪无遮挡| 国产主播在线观看一区二区| 搡老岳熟女国产| 涩涩av久久男人的天堂| 操出白浆在线播放| 中文字幕精品免费在线观看视频| 黄片小视频在线播放| 日本wwww免费看| 亚洲精品中文字幕在线视频| 69av精品久久久久久 | 精品第一国产精品| 午夜福利视频在线观看免费| 视频在线观看一区二区三区| 午夜日韩欧美国产| 成年人免费黄色播放视频| 欧美另类一区| 午夜视频精品福利| 午夜福利在线观看吧| 亚洲av成人一区二区三| 黑人巨大精品欧美一区二区mp4| 90打野战视频偷拍视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲五月色婷婷综合| 女性生殖器流出的白浆| svipshipincom国产片| 亚洲精品国产av蜜桃| 99国产精品一区二区蜜桃av | 电影成人av| 99久久综合免费| 黄片小视频在线播放| 爱豆传媒免费全集在线观看| 另类精品久久| 国产免费福利视频在线观看| 50天的宝宝边吃奶边哭怎么回事| 国产又爽黄色视频| 欧美人与性动交α欧美精品济南到| 欧美激情久久久久久爽电影 | 成年动漫av网址| 午夜久久久在线观看| 在线观看免费视频网站a站| 国产成人啪精品午夜网站| 80岁老熟妇乱子伦牲交| 国产一区二区 视频在线| 男女床上黄色一级片免费看| 一级毛片精品| 俄罗斯特黄特色一大片| 色婷婷久久久亚洲欧美| 久久天躁狠狠躁夜夜2o2o| 又紧又爽又黄一区二区| 一区二区三区精品91| 亚洲精品在线美女| 一进一出抽搐动态| 无遮挡黄片免费观看| av又黄又爽大尺度在线免费看| 欧美精品人与动牲交sv欧美| 正在播放国产对白刺激| 人人妻人人爽人人添夜夜欢视频| 亚洲色图综合在线观看| 日韩视频一区二区在线观看| 我要看黄色一级片免费的| 性色av乱码一区二区三区2| 欧美国产精品va在线观看不卡| 日本a在线网址| 亚洲一区中文字幕在线| 精品福利永久在线观看| 国产精品一二三区在线看| 久久久久国内视频| 12—13女人毛片做爰片一| 操出白浆在线播放| 国产亚洲av片在线观看秒播厂| 最新在线观看一区二区三区| 欧美97在线视频| 一级毛片电影观看| 亚洲男人天堂网一区| 50天的宝宝边吃奶边哭怎么回事| 成人亚洲精品一区在线观看| 99久久99久久久精品蜜桃| 精品国产乱码久久久久久男人| 热99国产精品久久久久久7| 美女大奶头黄色视频| 国产日韩欧美视频二区| 亚洲人成电影免费在线| 精品亚洲乱码少妇综合久久| 久久久国产一区二区| 国产精品国产av在线观看| 青青草视频在线视频观看| 亚洲精品国产精品久久久不卡| 国产免费一区二区三区四区乱码| 满18在线观看网站| 亚洲熟女毛片儿| 美女主播在线视频| 在线观看人妻少妇| 久久久精品94久久精品| 老汉色av国产亚洲站长工具| av天堂久久9| 国产不卡av网站在线观看| 淫妇啪啪啪对白视频 | 亚洲精品国产一区二区精华液| 啦啦啦免费观看视频1| e午夜精品久久久久久久| 悠悠久久av| 久久久久国产精品人妻一区二区| 人人妻人人澡人人爽人人夜夜| 国产精品一区二区免费欧美 | 丝袜人妻中文字幕| 在线观看免费高清a一片| 久久久欧美国产精品| 久久久久视频综合| 好男人电影高清在线观看| 桃红色精品国产亚洲av| 色婷婷久久久亚洲欧美| 啦啦啦 在线观看视频| 欧美久久黑人一区二区| 日本欧美视频一区| 国产一区二区在线观看av| 国产一级毛片在线| 欧美日韩成人在线一区二区| 曰老女人黄片| 免费女性裸体啪啪无遮挡网站| 国产av国产精品国产| 久久久欧美国产精品| 免费久久久久久久精品成人欧美视频| tube8黄色片| 女人精品久久久久毛片| 老熟女久久久| 亚洲国产av新网站| 成年动漫av网址| 色视频在线一区二区三区| 亚洲人成77777在线视频| 色婷婷久久久亚洲欧美| 两性夫妻黄色片| 精品亚洲乱码少妇综合久久| 欧美国产精品一级二级三级| 日韩三级视频一区二区三区| 亚洲va日本ⅴa欧美va伊人久久 | a级毛片黄视频| 久久九九热精品免费| 自线自在国产av| 成人三级做爰电影| 91精品三级在线观看| 一区二区三区乱码不卡18| 中亚洲国语对白在线视频| 黄色 视频免费看| 最新在线观看一区二区三区| 国产不卡av网站在线观看| 在线观看免费高清a一片| 满18在线观看网站| 欧美亚洲日本最大视频资源| 国产成人欧美在线观看 | 午夜福利乱码中文字幕| 啦啦啦中文免费视频观看日本| 最近中文字幕2019免费版| 国产精品99久久99久久久不卡| 久久香蕉激情| 国产一区有黄有色的免费视频| 欧美大码av| 亚洲av日韩精品久久久久久密| 久久久久网色| 别揉我奶头~嗯~啊~动态视频 | 久久久久久久久免费视频了| 叶爱在线成人免费视频播放| 麻豆乱淫一区二区| 操美女的视频在线观看| 热99久久久久精品小说推荐| 亚洲国产欧美网| 最近最新免费中文字幕在线| 自线自在国产av| 精品久久蜜臀av无| 国产淫语在线视频| www.精华液| xxxhd国产人妻xxx| 欧美97在线视频| 午夜久久久在线观看| 宅男免费午夜| 啦啦啦在线免费观看视频4| 人人妻人人澡人人爽人人夜夜| 午夜日韩欧美国产| 亚洲中文av在线| 国产成人av激情在线播放| 欧美日本中文国产一区发布| 日本av手机在线免费观看| 人人妻人人爽人人添夜夜欢视频| 国产国语露脸激情在线看| 中文字幕制服av| 日本一区二区免费在线视频| 999久久久国产精品视频| 国产欧美日韩综合在线一区二区| 国产成人欧美| 久久香蕉激情| 久久人人97超碰香蕉20202| 99久久国产精品久久久| 日本五十路高清| 日韩中文字幕视频在线看片| 精品一区二区三卡| 国产成人啪精品午夜网站| 另类亚洲欧美激情| 久久午夜综合久久蜜桃| 日本精品一区二区三区蜜桃| 亚洲精品久久午夜乱码| 国产在线免费精品| 王馨瑶露胸无遮挡在线观看| 免费观看av网站的网址| 欧美乱码精品一区二区三区| 黄频高清免费视频| 在线 av 中文字幕| 免费看十八禁软件| 国产一级毛片在线| 热99re8久久精品国产| av在线app专区| 可以免费在线观看a视频的电影网站| 日韩精品免费视频一区二区三区| av一本久久久久| 91精品三级在线观看| 午夜福利乱码中文字幕| 麻豆av在线久日| 日本av免费视频播放| 不卡一级毛片| 少妇人妻久久综合中文| 欧美97在线视频| 亚洲国产精品一区三区| 大香蕉久久网| 国产亚洲精品一区二区www | 国产欧美日韩一区二区精品| 亚洲全国av大片| 亚洲第一欧美日韩一区二区三区 | 丝袜在线中文字幕| 一区二区三区精品91| 国产高清视频在线播放一区 | 亚洲全国av大片| 18禁裸乳无遮挡动漫免费视频| 久久狼人影院| 青春草亚洲视频在线观看| 19禁男女啪啪无遮挡网站| 中文字幕高清在线视频| 欧美激情高清一区二区三区| 欧美 亚洲 国产 日韩一| 亚洲av欧美aⅴ国产| 一级a爱视频在线免费观看| 又大又爽又粗| 黄色 视频免费看| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 日韩大片免费观看网站| 精品视频人人做人人爽| 91大片在线观看| 久久亚洲国产成人精品v| 久久国产亚洲av麻豆专区| 99国产综合亚洲精品| 黄色视频在线播放观看不卡| 久久性视频一级片| 亚洲全国av大片| 亚洲国产欧美在线一区| 欧美av亚洲av综合av国产av| 欧美另类一区| 亚洲精品国产精品久久久不卡| 一本大道久久a久久精品| 日韩大码丰满熟妇| 久久国产精品男人的天堂亚洲| 亚洲精品av麻豆狂野| 视频在线观看一区二区三区| 亚洲欧美色中文字幕在线| 免费一级毛片在线播放高清视频 | 啦啦啦中文免费视频观看日本| 最近最新免费中文字幕在线| 免费黄频网站在线观看国产| 丁香六月欧美| 国产精品偷伦视频观看了| 日韩精品免费视频一区二区三区| 丰满少妇做爰视频| 最黄视频免费看| 久久亚洲国产成人精品v| 亚洲国产精品成人久久小说|