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

    基于高維多目標候鳥優(yōu)化算法的柔性作業(yè)車間調度

    2022-11-21 10:50:20王秋蓮段星皓
    中國機械工程 2022年21期
    關鍵詞:鳥群參考點飛鳥

    王秋蓮 段星皓

    南昌大學經濟管理學院,南昌,330031

    0 引言

    柔性作業(yè)車間調度問題(flexible job shop scheduling problem,F(xiàn)JSP)是指在一定的約束條件下,同時確定工件加工機器和加工順序,以達到理想的調度目標[1]。FJSP是NP-hard問題,大規(guī)模的FJSP很難通過數(shù)學規(guī)劃方法、拉格朗日松弛法等精確方法求得理論解,因此更多的研究采用啟發(fā)式算法對FJSP進行求解。如何改進啟發(fā)式算法使得其收斂性、穩(wěn)定性達到更優(yōu)是當前研究的熱點。

    FJSP根據(jù)調度目標的數(shù)量分為單目標FJSP和多目標FJSP(multi-objective FJSP,MOFJSP)。其中,MOFJSP可進一步細分為低維(調度目標不超過3個)MOFJSP和高維MOFJSP。目前,大部分學者主要對低維MOFJSP進行研究,如孫麗珍等[2]以最小化最大完工時間、機器負載的均方差為目標,設計了一種改進解碼方法和初始化規(guī)則的遺傳算法;孟冠軍等[3]提出一種混合人工蜂群算法求解以最大完工時間、最大機器負荷和機器總負荷為目標的低維MOFJSP。實際生產中,制造企業(yè)不僅需要考慮上述目標,還需要考慮能源消耗。李進宇等[4]研究了基于遞歸分析和圖像處理技術的混流車間能耗。李明等[5]提出一種新型帝國競爭算法來求解以最小化完工時間、最大拖期、最大機器負荷和總能耗為目標的高維MOFJSP,但該方法對高維多目標解集個體的選擇壓力不足,可能導致算法的過早收斂。黃夏寶等[6]提出一種改進的NSGAⅢ算法,引入基于參考點的小生境選擇機制來提高算法的多樣性,求解以最小化最大完工時間、機器總負荷、最大機器負荷和機器能耗為目標的高維MOFJSP,但該算法的局部搜索能力還有待提高。由此可見,考慮能耗目標的高維MOFJSP研究是當前車間調度問題研究的熱點及難點,但處于起步階段[7]。

    候鳥優(yōu)化(migrating birds optimization,MBO)算法[8]可用于求解二次分配問題,且所得的最終結果優(yōu)于遺傳算法、模擬退火算法等常用算法。QUAN等[9]將MBO算法用于求解混合流水車間的調度問題;任彩樂等[10]提出兩階段解碼方法及4種鄰域結構的改進MBO算法,求解混合流水車間調度問題。此外,有學者通過MBO算法來求解FJSP,如姚妮[11]以最小化最大完工時間為目標,結合MBO算法和變鄰域搜索策略來提高算法局部搜索能力;劉雪紅等[12]設計一種改進的MBO算法,引入精英分批策略和可行鄰域結構策略來解決以最小化最大完工時間和最大批次數(shù)目為目標的低維MOFJSP。

    盡管這些學者證明MBO算法對FJSP能獲得高質量解集,但用MBO算法進行高維MOFJSP的研究較少,因此本文建立一種以最小化完工時間、機器總負荷、延期時間、生產總能耗為優(yōu)化目標的MOFJSP模型,在改進MBO算法的基礎上,引入Pareto思想和基于參考點的選擇算子,設計了一種高維多目標候鳥優(yōu)化(multi-objective migrating birds optimization,MOMBO)算法,通過2種選擇算子給予鳥群選擇壓力,以提高鳥群的多樣性和收斂性,并基于屬性層次模型和灰色關聯(lián)度分析法的組合權重決策選擇最優(yōu)解。最后通過算例和實驗驗證算法的有效性。

    1 高維多目標柔性作業(yè)車間調度模型

    FJSP可以描述為:工件集J={J1,J2,…,Jn}中的n個工件在機器集M={M1,M2,…,Mm}中的m臺機器上加工,其中,第i個工件Ji有ei個工序。MOFJSP主要需要解決兩類子問題:機器選擇和工序排序,其中,機器選擇是指從可選機器集Mij(i∈[1,n],j∈[1,ei])中選擇一個機器完成工序Oij,工序排序是確定每個機器的加工順序。

    工件加工過程中,有以下假設:①所有機器在0時刻均可開始加工工件,所有工件在0時刻均被釋放;②每臺機器一次只能完成一道工序,每道工序在同一時刻只能由一臺機器完成;③工件的加工需嚴格按照工序依次完成;④加工一旦開始就不能中止,直到完成;⑤機器完成最后一道加工就停止運行;⑥加工過程中,工件的轉移時間忽略不計。

    本文提出的多目標柔性作業(yè)車間調度模型的調度目標的具體描述如下:

    (1)最大完工時間。最大完工時間是指所有工件加工完成的時間,即最后一個工件加工的完工時間,其計算公式為

    f1=max(Ci|i=1,2,…,n)

    (1)

    式中,Ci為工件i最后一道工序的完工時間。

    (2)總拖期。客戶和企業(yè)會約定產品的交貨期,交貨時間超過交貨期,企業(yè)會損失信譽并被罰款??偼掀诳梢苑从称髽I(yè)的信譽情況,其計算公式為

    (2)

    式中,Di為工件i的交貨期。

    (3)機器總負荷。機器總負荷指所有加工機器的負荷總和,可以反映企業(yè)資源的利用效率,其計算公式為

    (3)

    式中,tijk為工序Oij在機器k上的加工時間;xijk為0-1變量,xijk=1表示工序Oij在機器k上加工,xijk=0表示Oij不在機器k上加工。

    (4)總能耗??偰芎闹饕ㄜ囬g固定能耗、機器空載能耗、機器加工能耗,其計算公式為

    (4)

    (5)

    (6)

    Eg=Pgmax(Ci|i=1,2,…,n)

    (7)

    2 高維多目標候鳥優(yōu)化算法

    候鳥在飛行過程中會自然排成V字形隊列,這是因為領飛鳥通過翅膀扇動形成一股氣流,氣流向后方流動從而加大后方壓力,跟飛鳥利用壓力差飛行得更加輕松。MBO算法模仿候鳥的遷徙過程,將表現(xiàn)最優(yōu)的可行解作為領飛鳥,將其他可行解作為跟飛鳥。MBO算法主要分為鳥群初始化、領飛鳥進化、跟飛鳥進化和替換領飛鳥四個階段。本文在此基礎上,引入基于Pareto支配和基于參考點的選擇算子,對鳥群初始化、飛鳥進化階段進行改進,提出一種高維多目標候鳥優(yōu)化算法,該算法流程如圖1所示。

    圖1 高維多目標候鳥優(yōu)化算法流程圖

    2.1 基于Pareto支配和參考點的選擇規(guī)則

    本文引入基于Pareto支配的選擇算子和基于參考點的選擇算子來增強對Pareto前沿面的選擇壓力,進而對解進行篩選和替換。

    2.1.1基于Pareto支配的選擇算子

    Pareto支配的定義為:對于解p和q,當且僅當解p的所有目標值都不大于解q的目標值且至少存在一個目標值小于解q的目標值,則稱p支配q(pq)。

    定義si為鳥群中被個體i支配的個體集合,zi為鳥群中支配個體i的個體數(shù)量。

    基于上述定義,提出一種快速非支配排序算法,具體步驟如下:

    (1)將種群中不被其他個體支配的個體標記為第1級非支配前端F1;

    (2)F1中的個體p支配的個體集合為sp,對于sp中的個體q,將支配q的個體數(shù)zq減1即zq←zq-1后,若zq=0即不存在可以支配q的其他個體,則將q保存至集合Q中。對F1的每個支配個體集合中的全部個體執(zhí)行上述操作后,將最終得到的集合Q作為鳥群的第2級非支配前端F2;

    (3)遍歷F2中的個體,參照步驟(2)中的操作,找到第3級非支配前端F3。以此類推,直到種群中的每一個個體都分配到對應的非支配前端等級。

    通過基于Pareto思想的快速非支配排序算法將鳥群分為不同等級,等級越低說明該個體越優(yōu)。

    2.1.2基于參考點的選擇算子

    基于參考點的選擇算子首先在目標空間內采用單層參考點生成法[13]得到均勻分布的多個參考點,即將A個目標均勻劃分成H份,然后根據(jù)參考點依附的個體數(shù)量選擇較優(yōu)個體,具體流程如下:

    (8)

    然后通過下式歸一化鳥群個體的目標值:

    (9)

    歸一化鳥群個體的目標值后,此時理想點變?yōu)樵c,然后將理想點與各參考點的連線作為參考線,計算歸一化后的目標點與所有參考線的距離,選擇最近的參考線作為個體依附對象;

    (10)

    最后依次遍歷鳥群個體,統(tǒng)計每個參考點的依附個體數(shù)量即該參考點的小生境數(shù)。在MOMBO算法鳥群進化操作時,對于同一級非支配前端中的個體集合,優(yōu)先選擇小生境數(shù)更少的個體,保證算法的多樣性。

    2.2 多目標候鳥優(yōu)化算法的主要步驟

    2.2.1編碼

    運用MOMBO算法求解MOFJSP時首先要解決編碼問題,即將FJSP的潛在解轉換至可被MOMBO算法搜索的空間中。本算法采用二段式編碼規(guī)則[14],將解分為機器選擇(machine selection,MS)和工序排序(operation sequence,OS)兩部分。MS按工件的工序順序依次排列,每個位置的數(shù)值為該道工序選擇的加工機器號;OS中,每個位置的數(shù)值為工件號,從左往右遍歷OS,同一工件號出現(xiàn)第j次數(shù)表示對應工件的第j道工序。MS和OS如圖2所示。

    圖2 二段式編碼示意圖

    2.2.2解碼

    解碼則是將解空間轉換為問題空間。本算法采用貪婪插入式解碼法[15],以保證每道工序在所選加工機器上盡可能早地開始,對于前插的工序,它在OS上的位置也要相應調整。

    2.2.3鳥群初始化

    初始化是車間調度問題的一個重要步驟,它生成解集的好壞對算法的收斂性和多樣性有很大的影響。為保證初始解的質量和多樣性,本文分別對機器選擇和工序排序采用多種初始化規(guī)則。

    2.2.3.1 MS初始化規(guī)則

    (1)改進的全局選擇規(guī)則[16]。設定一個初始值為0的機器負荷數(shù)組,從工件集合中隨機選擇一個工件,從該工件的工序中隨機選擇一個工序,將該工序的可選機器加工時間與機器負荷數(shù)組的已有負荷相加,得到臨時機器負荷數(shù)組。從臨時機器負荷數(shù)組中選擇臨時負荷最小的機器作為該工序的加工機器;若臨時負荷最小的機器不只一個,則選擇加工時間最短的機器作為該工序的加工機器;更新機器負荷數(shù)組。當前工件的所有工序選擇完加工機器后,隨機選擇下一個工件,依次類推,直到所有工件的工序都選擇完機器。

    (2)改進的局部選擇規(guī)則[16]。設定一個初始值為0的機器負荷數(shù)組,從工件集合中依次選擇一個工件,從該工件的工序集合中隨機選擇一個工序,將該工序的可選機器加工時間與機器負荷數(shù)組的已有負荷進行相加,得到臨時機器負荷數(shù)組。從臨時機器負荷數(shù)組中選擇臨時負荷最小的機器作為該工序的加工機器;若臨時負荷最小的機器不只一個,則選擇加工時間最短的機器作為該工序的加工機器;更新機器負荷數(shù)組。當前工件的所有工序選擇完加工機器后,將機器負荷數(shù)組重新置0,選擇下一個工件,依次類推,直到所有工件的工序都選擇完機器。

    (3)隨機選擇。在加工工序的可選機器中隨機選擇一臺機器進行加工。

    2.2.3.2 OS初始化規(guī)則

    (1)剩余負荷最大規(guī)則[17]。將工件按剩余工序加工時間總和進行降序排列,優(yōu)先加工剩余加工時間長的工件。

    (2)最短加工時間規(guī)則[18]。將工序按照加工時間升序排列,優(yōu)先完成加工時間短的工序。

    (3)隨機選擇。隨機選擇工序進行加工。

    參考文獻[16,19],首先按照鳥總數(shù)(即個體總數(shù))4∶4∶2的比例,分別通過改進的全局選擇規(guī)則、改進的局部選擇規(guī)則和隨機選擇規(guī)則生成初始鳥群的MS,再按照上述比例生成初始鳥群的OS。

    初始化之后,按如下步驟將鳥群分為領飛鳥和左右跟飛鳥群:首先根據(jù)基于Pareto支配的選擇算子進行非支配排序,找到初始種群中的非支配解集;然后通過基于參考點的選擇算子,選擇小生境數(shù)最小的非支配解作為領飛鳥,若小生境數(shù)最小的非支配解有多個,則隨機選擇一個作為領飛鳥,并將剩下個體按非支配等級均勻分為左右跟飛鳥群。

    2.2.4鳥群進化

    首先定義三種鄰域結構:

    N1:隨機產生i(i=2,3,4)個變異工序。將i個工序順序重新排列組合,然后計算不同排列組合的目標值,通過上文中的兩種選擇算子選擇其中最好的新解作為鄰域解。

    N2:基于工序變異,交換兩個不同工件的任意工序位置。

    N3:基于機器變異,隨機選擇兩個變異工件,并在這兩個變異工件上分別隨機選擇一道工序,在不考慮原加工機器的情況下,重新選擇這兩道工序加工時間最短的機器。

    基于上述3種鄰域結構,增加鳥群個體的局部尋優(yōu)能力,鳥群進化具體步驟如下。

    2.2.4.1 領飛鳥進化

    按上述3種鄰域結構產生領飛鳥的3個鄰域解后,將領飛鳥與3個鄰域解合并作為進化解集,接著從該進化解集中選擇1個個體作為新的領飛鳥,具體的選擇規(guī)則如下:①若鄰域解支配領飛鳥,則替換領飛鳥;②若領飛鳥支配鄰域解,則不替換;③若領飛鳥與鄰域解互不支配,則進行領飛鳥選擇,具體的領飛鳥選擇操作如下。

    首先將原鳥群中的進化個體剔除,計算剔除進化個體后的鳥群的小生境數(shù);然后找出具有最小小生境數(shù)的參考點集J,遍歷J中的每個參考點。假設Ja是J中的一個參考點,判斷進化解集(包括進化個體、進化個體的鄰域解)的非支配解集中是否有個體依附于參考點Ja,若進化解集的非支配解集中沒有個體依附于參考點Ja,則將參考點Ja的小生境數(shù)ρa設為無窮大。若進化解集的非支配解集中有個體依附于參考點Ja,則進一步判斷ρa是否為0,若ρa為0,則從進化解集的非支配解集中依附Ja的個體中將與參考點Ja對應參考線距離最短的個體作為新的領飛鳥加入鳥群;若ρa不為0,則在進化解集的非支配解集中隨機選擇1個依附該參考點的個體作為新的領飛鳥加入鳥群。進化過程中選擇領飛鳥的流程如圖3所示。

    圖3 進化過程中基于參考點的選擇算子流程圖

    更新完領飛鳥之后,選擇一個次優(yōu)解作為第一個跟飛鳥的共享解,以增強跟飛鳥群的尋優(yōu)能力。次優(yōu)解的選擇過程如下:判斷進化解集中剩余個體的非支配等級,選擇非支配等級最低的個體作為次優(yōu)解;若非支配等級最低的個體數(shù)大于1,則比較其依附參考點在鳥群中的小生境數(shù),優(yōu)先選擇小生境數(shù)更小的個體,若小生境數(shù)相同,則隨機選擇一個體作為次優(yōu)解。

    2.2.4.2 跟飛鳥進化

    跟飛鳥進化包括對左右跟飛鳥群的進化,即首先將領飛鳥進化解集中的次優(yōu)解作為左右跟飛鳥群中第一只跟飛鳥的共享解,將共享解與生成的3個鄰域解組成該跟飛鳥的進化解集。然后,選擇進化解集中的最優(yōu)解作為新的跟飛鳥,同時將進化解集中的次優(yōu)解作為下一只跟飛鳥的共享解,替換跟飛鳥和選擇次優(yōu)解規(guī)則與上文中領飛鳥的進化過程相同。按上述步驟依次遍歷所有跟飛鳥,完成跟飛鳥的進化。

    2.2.5跟飛鳥交叉

    每次對左右跟飛鳥群執(zhí)行完進化操作后,從左右跟飛鳥群中隨機選擇Cr(Cr為交叉?zhèn)€數(shù))對個體,并對每對個體依次進行交叉操作。

    MS、OS的交叉操作分別采用均勻交叉策略和POX(precedence operation crossover)交叉策略[15]。MS的均勻交叉策略是:隨機產生一個長度為工序總數(shù)的0、1數(shù)組,交換父代數(shù)字1對應位的機器號,產生2個新的子代。OS的POX交叉策略是:首先將工件隨機地分為子集合J1和J2;然后將父代1中J1部分的工件號復制到子代1對應位置,將父代2中J2部分的工件號復制到子代2中對應位置;最后將父代2中J2部分的工件號按順序依次填入子代1中的剩余位,將父代1中J1部分的工件號按順序依次填入子代2中的剩余位。MS交叉如圖4所示,OS交叉如圖5所示。

    圖4 MS均勻交叉圖

    圖5 OS POX交叉圖

    2.2.6篩選個體產生新鳥群

    個體交叉生成子代鳥群Q后,將子代加入至原鳥群Y(數(shù)量為n)中,將Q和Y合并為R,對R重新進行非支配排序,在非支配解集中隨機選擇小生境數(shù)最小的個體作為新鳥群S的領飛鳥,同時更新S的小生境數(shù);然后按照非支配層等級由低到高將每層個體加入至S中,直至S的數(shù)量等于n;若S的飛鳥數(shù)量大于n,則需要重新在最后一層非支配層中通過基于參考點的選擇算子選擇一定數(shù)量個體保留在S中,未被選擇個體則剔除出S,使S的數(shù)量等于n。新鳥群生成如圖6所示。

    圖6 新鳥群生成圖

    生成新鳥群之后,若未達到最大迭代次數(shù),則返回至鳥群進化階段進行下一次迭代;若達到最大迭代次數(shù),則將新鳥群的非支配解集作為最終優(yōu)化結果輸出。

    2.2.7基于AHM和灰色關聯(lián)度分析法的組合權重決策

    實際生產需要選擇一個確定的方案。MOMBO最終的非支配解集可能存在多個解,本節(jié)采用組合賦權法進行方案選擇。組合賦權是主觀賦權和客觀賦權的線性加權,其中,主觀賦權采用屬性層次模型,客觀賦權采用灰色關聯(lián)度分析法。

    (1)屬性層次模型(attribute hierarchical model,AHM)先將各個指標劃分為相互關聯(lián)的有序結構,再結合專家的主觀意見和分析者的判斷來確定指標權重?;贏HM 的權重計算方法如下:

    ①建立判斷矩陣A=[aij],比較同一等級指標之間的重要性,其中,aij=1/aji(i≠j),aii=1。

    ②設μ=[μij]為測度矩陣,μij是由aij換算而成的測度指標,令

    (11)

    k=1,2,…

    這里設置β=1。

    ③計算指標權重。指標主觀權重向量Wa=(w1,w2,…,wm)按下式計算:

    (12)

    其中,w1+w2+…+wm=1且0≤wi≤1。

    (2)灰色關聯(lián)度分析(gray relation analysis,GRA)法的思想是:先根據(jù)某個問題的實際情況確定理想的最優(yōu)序列;然后,通過方案的序列曲線和幾何形狀與理想最優(yōu)序列的曲線和幾何形狀的相似程度來判斷方案的序列與理想最優(yōu)序列之間的關聯(lián)程度,曲線和幾何形狀越接近說明關聯(lián)度越大,方案越接近理想最優(yōu),則該方案的權重越大。GRA的主要步驟如下:

    ①數(shù)據(jù)規(guī)范化處理。假設有n個方案,每個方案有m個目標值,首先對各個指標進行規(guī)范化處理,規(guī)范化處理計算式為

    (13)

    ②由下式計算灰色關聯(lián)系數(shù):

    (14)

    灰色關聯(lián)系數(shù)反映yij與理想值之間的關聯(lián)程度。

    ③權重計算。根據(jù)灰色關聯(lián)系數(shù)矩陣計算得到客觀權重向量Wb,其中,第j個目標的權重為

    (15)

    (16)

    (3)綜合權重確定。將AHM得出的權重Wa與灰色關聯(lián)度分析法的權重Wb結合,求得綜合權重:

    w*=θWa+(1-θ)Wb

    (17)

    其中,θ為主觀偏重系數(shù),θ∈[0,1],由決策者依據(jù)偏好給出。組合賦權法既能充分利用專家長期積累的生產加工經驗,使結果更貼近實際生產,又能在一定程度上從數(shù)據(jù)本身特性來確定其重要性。

    3 性能測試與案例研究

    為驗證本算法的有效性,首先將MOMBO分別與三個變體、NSGAⅡ及NSGAⅢ的最終解集進行對比,然后通過仿真對MOMBO求得的非支配解集進行分析。本節(jié)所有算法均由MATLAB R2014實現(xiàn),MATLAB在Intel i5-8300H、內存8.00 GB的計算機上運行。

    3.1 算法性能測試

    本小節(jié)通過10個BRdata測試問題[20]來驗證MOMBO的性能,將測試問題的加工時間作為本調度問題中各個工序的基本加工時間,算例的其他數(shù)據(jù)如交貨期、機器加工功率、空載功率等數(shù)據(jù)如表1所示。

    表1 算例數(shù)據(jù)

    多目標優(yōu)化算法的解集有多種評價指標,本文采用的評價指標是反向世代距離(inverted generational distance,IGD)[21]和超體積(hypervolume,HV)[22]。

    IGD能綜合評價解集的收斂性和多樣性,其計算公式為

    (18)

    其中,P*為一組理想解集,本文設定為所有算法多次運算后得到的全部解中的非支配解集;|P*|為解集P*中解的個數(shù);x為解集P*中的某個解;解集P是某算法求出的該前沿面的一個近似解集,y是解集P中的某個解;d(x,y)為解x和y在目標空間中的歐氏距離。IGD值越小,算法性能越好。

    超體積能綜合反映解集的收斂性和多樣性,其計算公式為

    (19)

    其中,解集P是Pareto前沿面的近似解集;fm為解集P中某個非支配解f=(f1,f2,…,fm)T的第m個目標值;r=(r1,r2,…,rm)T為目標空間中的一個參考向量,它被解集P中所有目標向量支配;volume(*)表示近似Pareto前沿面與參考向量所圍成的超立方體的超體積。那么關于參考向量r的超體積指的是被解集P所支配且以參考向量r為邊界的目標空間的體積。

    將超體積定義為被解集P所支配且以參考點r為邊界的目標空間的體積,HV(P,r)越大,算法性能越好。對于4個目標的柔性作業(yè)車間調度問題,本文采用BADER等[23]提出的蒙特卡羅模擬方法計算HV(P,r)。

    3.1.1MOMBO與MOMBO變體的比較

    為驗證所做改進對MOMBO算法性能的影響,構建以下幾種變體類型:

    MOMBO1:采用半主動解碼。

    MBMBO2:種群初始化階段采用隨機初始化方式。

    MOMBO3:飛鳥進化過程僅采用N1鄰域結構產生鄰域解。

    MOMBO算法參數(shù)設置如表2所示。多次運行MOMBO及其3種變體,求得多個近似解集,然后計算IGD值和HV值,最終取平均值進行對比,結果如表3所示,其中,最優(yōu)的結果用加粗數(shù)字表示。由表3可知,除了算例MK04和MK05,MOMBO算法的IGD值、HV值都優(yōu)于3種變體算法,由此可知基于MBO算法提出的貪婪插入式解碼法、改進的多種初始化規(guī)則,以及改進的鄰域搜索策略在大部分情況下更適用于高維多目標柔性作業(yè)車間調度模型。

    表2 MOMBO參數(shù)設置

    表3 MOMBO及其變體的IGD、HV

    3.1.2MOMBO與NSGAⅡ、NSGAⅢ的比較

    NSGAⅡ、NSGAⅢ的參數(shù)設置如表4所示。多次運行MOMBO、NSGAⅡ和NSGAⅢ,求得多個近似解集,然后計算IGD值和HV值,最終取平均值進行對比,結果如表5所示,其中,最優(yōu)的結果用加粗數(shù)字表示。由表5可知,MOMBO的IGD值和HV值都顯著優(yōu)于NSGAⅡ和NSGAⅢ,驗證了MOMBO求解該問題的有效性。

    表4 NSGAⅡ、NSGAⅢ的參數(shù)設置

    表5 MOMBO、NSGAⅡ和NSGAⅢ的IGD、HV

    3.1.3算法收斂性分析

    為驗證算法運行過程中的收斂性,繪制MOMBO、3種變體、NSGAⅡ和NSGAⅢ的HV值隨迭代次數(shù)變化的進化軌跡圖。由圖7可以看出,6種算法的HV值都隨著迭代次數(shù)的增大而增大,這意味著它們的收斂性能穩(wěn)定;MOMBO每代的HV值在MK08、MK09算例上遠大于其他算法的HV值,說明MOMBO的進化效果最好;MOMBO2、NSGAⅡ、NSGAⅢ的HV值明顯小于MOMBO、MOMBO1、MOMBO的HV值,這說明本文提出的多種初始化方式是有效的。接下來考慮非支配解集占比情況,MOMBO求解MK01的初始解集中,非支配等級共有5層,其中,非支配解集占所有解集的比例為28.6%;MOMBO求解MK01的最終解集的非支配等級為4層,其中,非支配解集占所有解集的68.6%。這說明所提出的兩種選擇算子能給予算法有效的進化壓力,能有效促進算法的收斂。

    (a)MK04 (b)MK08 (c)MK09

    3.2 實例研究

    為驗證本文算法在實際生產加工中的效果,在某工廠進行加工實驗。實驗內容是,5臺加工機器上加工4個工件,通過Precision Power Analyzer LMG671收集機床加工工件過程中的功率。對得到的數(shù)據(jù)進行預處理,獲得每臺機器的加工功率和空載功率曲線圖。為方便計算,通過取均值求得各機器的加工功率和空載功率。車間固有功率為2kW,每個工件的交貨期和加工時間如表6所示,其中,“-”表示該工序不能在此機器上完成。

    表6 實驗數(shù)據(jù)

    3.2.1MOMBO與單目標MBO的比較

    通過MOMBO算法求得Pareto前沿解集,如表7所示。通過MBO算法分別求解出以最大完工時間f1、總拖期f2、機器最大負荷f3、總能耗f4為目標的單目標最優(yōu)解,如表8所示。觀察表8可以發(fā)現(xiàn),雖然單目標最優(yōu)解在某個目標值上最優(yōu),但其他3個目標值較差,單目標最優(yōu)解與通過MOMBO算法求出的最優(yōu)解集中的某個解的目標值之差即Δf1~Δf4如表9所示。由表9可見,在MOMBO求出的最優(yōu)解集中,總能找出優(yōu)于單目標最優(yōu)解的解,因此通過MOMBO求出的最優(yōu)解集相對于單目標最優(yōu)解具有更好的質量。

    表7 實例Pareto前沿解集

    表8 實例單目標最優(yōu)解

    表9 單目標最優(yōu)解與多目標最優(yōu)解的目標值之差

    3.2.2MOMBO與NSGAⅡ、NSGAⅢ的比較

    NSGAⅡ、NSGAⅢ求得的最優(yōu)解集分別如表10、表11所示。對比表7可以看出,MOMBO求出的最優(yōu)解大部分都優(yōu)于NSGAⅡ、NSGAⅢ求得的解。為進一步對比MOMBO與NSGAⅡ、NSGAⅢ的收斂性和多樣性,通過式(18)、式(19)求得3種算法的IGD值和HV值,如表12所示。在求解過程中,本文設定P*為三種算法多次運算后得到的非支配解集。如表12所示,對于本實例,MOMBO取得了優(yōu)于另外2個算法的結果。

    表10 NSGAⅡ的實例Pareto前沿解集

    表11 NSGAⅢ的實例Pareto前沿解集

    表12 MOMBO、NSGAⅡ和NSGAⅢ求解實例的結果

    圖8、圖9所示為3種算法求解實例的HV值和IGD值隨迭代次數(shù)的進化曲線。3種算法的HV值都隨著迭代次數(shù)的增大而增大,IGD值均隨著迭代次數(shù)的增大而減小,這說明它們的收斂性能穩(wěn)定。本文提出的MOMBO的HV值大于另外兩種算法的HV值,IGD值小于另外兩種算法的IGD值,因此MOMBO的進化效果最好。

    圖8 3種算法求解實例所得HV值隨迭代次數(shù)的進化軌跡

    圖9 3種算法求解實例所得IGD值隨迭代次數(shù)的進化軌跡

    3.2.3基于組合權重法的最優(yōu)方案決策

    首先根據(jù)專家對不同目標的打分情況建立AHM需要的判斷矩陣,然后通過計算得到各個目標的主觀權重向量Wa=(0.264,0.361,0.264,0.111);接著根據(jù)表7中的Pareto前沿解集計算出灰色關聯(lián)度分析法的各個數(shù)據(jù):

    Wb=(0.261,0.189,0.281,0.269)

    進一步取θ=0.5,求得綜合權重w*=(0.263,0.275,0.272,0.190)。由于本文是逆指標,故需要對y采用減法一致法處理得到y(tǒng)′,即利用每一列的最大值Mj依次減去該列的原始數(shù)據(jù),然后將y′×w*,得到最終得分矩陣:

    由R可知第二組方案的得分最高,其加工時間為53.42 s,拖期為0 s,機器負載為178.50 s,總能耗為364 681 J。圖10為按第二個調度方案加工的甘特圖。

    圖10 解2的甘特圖

    4 結論

    本文提出以最小化最大完工時間、總拖期、機器總負荷、總能耗為目標的高維多目標柔性作業(yè)車間調度模型,并運用一種MOMBO進行求解。算法性能測試表明MOMBO相對于其3種變體以及NSGAⅡ、NSGAⅢ來說,其最優(yōu)解集在IGD和HV兩種性能指標上表現(xiàn)更優(yōu)。實例研究證明本文提出的MOMBO算法相對于單目標MBO算法、NSGAⅡ和NSGAⅢ有著更好的解集質量,在實際生產加工中能夠給予決策者更好的選擇。

    猜你喜歡
    鳥群參考點飛鳥
    飛翔的鳥群
    飛鳥與少年
    心聲歌刊(2020年5期)2020-11-27 01:56:10
    FANUC數(shù)控系統(tǒng)機床一鍵回參考點的方法
    為什么鳥要成群飛翔?
    為什么鳥群飛行時不會彼此沖撞?
    參考點對WiFi位置指紋算法的影響
    測控技術(2018年5期)2018-12-09 09:04:24
    數(shù)控機床返回參考點故障維修
    飛鳥
    島與飛鳥
    螢火(2016年7期)2016-07-19 08:18:54
    FANUC數(shù)控機床回參考點故障分析與排除
    日韩熟女老妇一区二区性免费视频| 亚洲精品国产av蜜桃| 久久青草综合色| 国产老妇伦熟女老妇高清| 飞空精品影院首页| 尾随美女入室| 亚洲精品,欧美精品| 国产精品人妻久久久影院| 卡戴珊不雅视频在线播放| a 毛片基地| 亚洲欧美色中文字幕在线| 久久精品人人爽人人爽视色| 人人妻人人澡人人看| 人人澡人人妻人| 一级a爱视频在线免费观看| 丝袜喷水一区| 久久久久久久久久久免费av| 又大又黄又爽视频免费| 黄色怎么调成土黄色| 丝袜脚勾引网站| av在线老鸭窝| 亚洲七黄色美女视频| 日本vs欧美在线观看视频| av在线播放精品| 亚洲三区欧美一区| a级片在线免费高清观看视频| 午夜福利网站1000一区二区三区| 日韩av免费高清视频| 亚洲欧美清纯卡通| 国产亚洲av片在线观看秒播厂| 五月开心婷婷网| 中文欧美无线码| 久久久久精品国产欧美久久久 | 亚洲国产av影院在线观看| 国产av国产精品国产| 精品国产超薄肉色丝袜足j| 80岁老熟妇乱子伦牲交| 超色免费av| 国产精品 欧美亚洲| 欧美国产精品va在线观看不卡| 97在线人人人人妻| 久久精品国产亚洲av高清一级| 免费不卡黄色视频| 日韩欧美精品免费久久| 国产一区有黄有色的免费视频| 天天添夜夜摸| 国产精品亚洲av一区麻豆 | 中文天堂在线官网| 制服人妻中文乱码| 另类精品久久| 精品人妻熟女毛片av久久网站| 久久99精品国语久久久| 校园人妻丝袜中文字幕| 国产欧美日韩一区二区三区在线| 免费高清在线观看日韩| 极品人妻少妇av视频| 美女福利国产在线| 日日啪夜夜爽| 在线天堂最新版资源| 亚洲欧美成人综合另类久久久| 最近中文字幕高清免费大全6| 天堂俺去俺来也www色官网| 精品人妻在线不人妻| 天天操日日干夜夜撸| 我要看黄色一级片免费的| av女优亚洲男人天堂| 两性夫妻黄色片| 超色免费av| 久久久精品94久久精品| 午夜福利视频精品| 亚洲三区欧美一区| 亚洲国产精品国产精品| 久久亚洲国产成人精品v| 国产伦人伦偷精品视频| 免费在线观看黄色视频的| 精品一区在线观看国产| 我要看黄色一级片免费的| 九九爱精品视频在线观看| 深夜精品福利| 欧美xxⅹ黑人| 久久精品亚洲av国产电影网| 视频区图区小说| 日韩大片免费观看网站| 校园人妻丝袜中文字幕| 美女福利国产在线| 好男人视频免费观看在线| 色婷婷久久久亚洲欧美| 日韩中文字幕欧美一区二区 | 亚洲欧美精品综合一区二区三区| 热re99久久精品国产66热6| 悠悠久久av| 久久国产精品男人的天堂亚洲| 亚洲欧美成人精品一区二区| 欧美最新免费一区二区三区| 在线看a的网站| 欧美人与性动交α欧美精品济南到| 久久99热这里只频精品6学生| 欧美亚洲日本最大视频资源| 99精国产麻豆久久婷婷| 夫妻午夜视频| 超碰97精品在线观看| 国产探花极品一区二区| 精品一区二区免费观看| 亚洲精品日本国产第一区| 欧美97在线视频| 天天操日日干夜夜撸| 免费av中文字幕在线| 精品少妇一区二区三区视频日本电影 | 日韩欧美一区视频在线观看| 亚洲人成电影观看| 美女大奶头黄色视频| √禁漫天堂资源中文www| 在线免费观看不下载黄p国产| 蜜桃在线观看..| 在线观看免费日韩欧美大片| 丰满少妇做爰视频| 精品一区二区三卡| 亚洲美女黄色视频免费看| 亚洲精品av麻豆狂野| 丝瓜视频免费看黄片| 天天添夜夜摸| 成人影院久久| 亚洲伊人色综图| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产中文字幕在线视频| 久久久久精品性色| 蜜桃国产av成人99| 欧美人与性动交α欧美软件| 一本—道久久a久久精品蜜桃钙片| 桃花免费在线播放| 久久精品久久久久久噜噜老黄| av网站在线播放免费| 久久久久久久精品精品| 日本欧美国产在线视频| 自线自在国产av| 18禁国产床啪视频网站| 欧美激情高清一区二区三区 | 久久韩国三级中文字幕| 最近的中文字幕免费完整| 久久久国产欧美日韩av| 人人妻人人爽人人添夜夜欢视频| 母亲3免费完整高清在线观看| 操美女的视频在线观看| 日韩 亚洲 欧美在线| 免费在线观看完整版高清| 99香蕉大伊视频| 一本久久精品| 国产毛片在线视频| 18禁裸乳无遮挡动漫免费视频| 九色亚洲精品在线播放| 亚洲,欧美,日韩| 国产国语露脸激情在线看| 日韩免费高清中文字幕av| 精品一区二区三卡| 国产精品偷伦视频观看了| 久久国产亚洲av麻豆专区| videos熟女内射| 天天操日日干夜夜撸| 黑人欧美特级aaaaaa片| 午夜精品国产一区二区电影| 交换朋友夫妻互换小说| 国产成人a∨麻豆精品| 亚洲国产毛片av蜜桃av| 日韩av不卡免费在线播放| 人人澡人人妻人| 极品少妇高潮喷水抽搐| 制服诱惑二区| www日本在线高清视频| 国产又色又爽无遮挡免| 99久久精品国产亚洲精品| 美女国产高潮福利片在线看| 午夜福利,免费看| 亚洲av成人精品一二三区| 女性被躁到高潮视频| 18禁国产床啪视频网站| 亚洲美女视频黄频| 久久久久国产精品人妻一区二区| 90打野战视频偷拍视频| 日本猛色少妇xxxxx猛交久久| 一区二区三区激情视频| 久久精品国产亚洲av高清一级| 精品第一国产精品| 巨乳人妻的诱惑在线观看| 国产极品天堂在线| 欧美人与性动交α欧美软件| www.av在线官网国产| 日韩中文字幕视频在线看片| 女的被弄到高潮叫床怎么办| 日韩人妻精品一区2区三区| 青青草视频在线视频观看| www.自偷自拍.com| 天天操日日干夜夜撸| 久久精品久久精品一区二区三区| 一区福利在线观看| 亚洲国产av新网站| 亚洲精品在线美女| 男女免费视频国产| 毛片一级片免费看久久久久| 一本一本久久a久久精品综合妖精| 美女主播在线视频| 一级a爱视频在线免费观看| 精品酒店卫生间| 亚洲国产欧美在线一区| 男人舔女人的私密视频| 亚洲精品久久久久久婷婷小说| 午夜日韩欧美国产| 久久av网站| 中国三级夫妇交换| 18禁动态无遮挡网站| 国产免费福利视频在线观看| 久久国产亚洲av麻豆专区| 在线观看免费视频网站a站| av福利片在线| 欧美日韩精品网址| 最近中文字幕2019免费版| 搡老乐熟女国产| 女人高潮潮喷娇喘18禁视频| 大香蕉久久成人网| 超碰成人久久| 亚洲人成77777在线视频| www.精华液| 自线自在国产av| 视频在线观看一区二区三区| 国产精品熟女久久久久浪| 国产伦理片在线播放av一区| 一边亲一边摸免费视频| 飞空精品影院首页| 麻豆av在线久日| 亚洲婷婷狠狠爱综合网| 日韩,欧美,国产一区二区三区| 大码成人一级视频| 亚洲图色成人| 国产 一区精品| 亚洲国产精品成人久久小说| 欧美日韩综合久久久久久| 精品一区在线观看国产| 美女中出高潮动态图| 观看av在线不卡| 久久久久精品性色| 天天添夜夜摸| 国产免费视频播放在线视频| 午夜福利,免费看| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲国产欧美日韩在线播放| 又大又爽又粗| a级毛片黄视频| 欧美最新免费一区二区三区| 中文字幕av电影在线播放| 精品一品国产午夜福利视频| 日韩大片免费观看网站| 亚洲av成人不卡在线观看播放网 | 人体艺术视频欧美日本| 高清不卡的av网站| 99国产精品免费福利视频| 精品午夜福利在线看| 成年av动漫网址| 视频在线观看一区二区三区| 亚洲欧美中文字幕日韩二区| 精品国产乱码久久久久久男人| 国产男人的电影天堂91| 永久免费av网站大全| 成人影院久久| 亚洲精品久久成人aⅴ小说| 成人国产麻豆网| 少妇 在线观看| 亚洲成色77777| av又黄又爽大尺度在线免费看| 狂野欧美激情性bbbbbb| 国产一卡二卡三卡精品 | 国产成人91sexporn| 国产片内射在线| 国产精品免费视频内射| 黄片小视频在线播放| 成人手机av| 黄片播放在线免费| 热re99久久精品国产66热6| 在线精品无人区一区二区三| 欧美日韩福利视频一区二区| 美女福利国产在线| 亚洲国产精品国产精品| 免费观看av网站的网址| 国产一区亚洲一区在线观看| 中文字幕精品免费在线观看视频| 国产爽快片一区二区三区| 黑丝袜美女国产一区| 成人三级做爰电影| 欧美精品高潮呻吟av久久| 别揉我奶头~嗯~啊~动态视频 | 好男人视频免费观看在线| 在线精品无人区一区二区三| 黄片播放在线免费| 免费在线观看完整版高清| 国产男女内射视频| 日韩电影二区| 国产男人的电影天堂91| 午夜影院在线不卡| 最新在线观看一区二区三区 | 亚洲欧美色中文字幕在线| 亚洲欧美成人综合另类久久久| 另类精品久久| 亚洲美女黄色视频免费看| 国产黄色视频一区二区在线观看| 大话2 男鬼变身卡| 老司机深夜福利视频在线观看 | 亚洲男人天堂网一区| 欧美中文综合在线视频| 亚洲精品久久成人aⅴ小说| 亚洲精品乱久久久久久| 天天躁狠狠躁夜夜躁狠狠躁| 一本色道久久久久久精品综合| 自线自在国产av| 午夜av观看不卡| 最近2019中文字幕mv第一页| 天天操日日干夜夜撸| 欧美人与性动交α欧美软件| 国产伦人伦偷精品视频| 免费在线观看视频国产中文字幕亚洲 | 国产无遮挡羞羞视频在线观看| 美女大奶头黄色视频| 欧美 亚洲 国产 日韩一| 看十八女毛片水多多多| 亚洲,一卡二卡三卡| 亚洲精品一二三| 夫妻午夜视频| 另类亚洲欧美激情| 国产成人啪精品午夜网站| 亚洲国产毛片av蜜桃av| av又黄又爽大尺度在线免费看| 日日啪夜夜爽| 色吧在线观看| 欧美精品高潮呻吟av久久| 免费不卡黄色视频| 在线观看人妻少妇| 国产一区二区在线观看av| 女的被弄到高潮叫床怎么办| 久久久久精品久久久久真实原创| 免费黄网站久久成人精品| 99香蕉大伊视频| 亚洲精品国产区一区二| 69精品国产乱码久久久| 狠狠精品人妻久久久久久综合| 只有这里有精品99| 咕卡用的链子| avwww免费| 日韩不卡一区二区三区视频在线| 国产精品熟女久久久久浪| 免费在线观看完整版高清| 久久久久久久久久久久大奶| 老熟女久久久| 久久久欧美国产精品| 亚洲欧洲国产日韩| 久久久久网色| 欧美老熟妇乱子伦牲交| 十分钟在线观看高清视频www| 久久久精品区二区三区| 国产一区二区激情短视频 | 国产野战对白在线观看| 国产精品国产三级专区第一集| 我要看黄色一级片免费的| www.av在线官网国产| 成人18禁高潮啪啪吃奶动态图| 国产深夜福利视频在线观看| 国产精品久久久久久人妻精品电影 | 精品国产一区二区久久| av女优亚洲男人天堂| 人人妻人人澡人人爽人人夜夜| 我的亚洲天堂| 国产一区有黄有色的免费视频| 国产极品天堂在线| 欧美乱码精品一区二区三区| 国产在线视频一区二区| 欧美少妇被猛烈插入视频| 国产深夜福利视频在线观看| 丝袜人妻中文字幕| 99久久99久久久精品蜜桃| 精品人妻在线不人妻| 美女主播在线视频| 水蜜桃什么品种好| 亚洲中文av在线| 亚洲精品国产一区二区精华液| 亚洲精品第二区| 亚洲精品一二三| 日本wwww免费看| 国产精品99久久99久久久不卡 | 亚洲国产精品一区三区| 一个人免费看片子| 欧美日韩视频精品一区| 下体分泌物呈黄色| 国产精品一区二区在线观看99| 久久久精品94久久精品| 十八禁高潮呻吟视频| 欧美xxⅹ黑人| 国产精品 国内视频| 一二三四中文在线观看免费高清| 亚洲男人天堂网一区| 少妇人妻久久综合中文| 另类精品久久| 欧美精品一区二区大全| 久久精品国产综合久久久| a级片在线免费高清观看视频| 久久久久国产一级毛片高清牌| 国产成人欧美| 电影成人av| 久久青草综合色| 美女中出高潮动态图| 免费在线观看黄色视频的| 最近的中文字幕免费完整| 国产av精品麻豆| 校园人妻丝袜中文字幕| 精品人妻在线不人妻| 满18在线观看网站| 成人18禁高潮啪啪吃奶动态图| 一个人免费看片子| 99热国产这里只有精品6| 99热全是精品| 最近最新中文字幕大全免费视频 | 男女边吃奶边做爰视频| 日韩一区二区视频免费看| av一本久久久久| 久久亚洲国产成人精品v| 综合色丁香网| 久久精品人人爽人人爽视色| 免费在线观看视频国产中文字幕亚洲 | 熟妇人妻不卡中文字幕| 一本色道久久久久久精品综合| 男人添女人高潮全过程视频| 久久影院123| 国产在线视频一区二区| 欧美乱码精品一区二区三区| 交换朋友夫妻互换小说| 精品国产超薄肉色丝袜足j| 婷婷色综合大香蕉| 高清在线视频一区二区三区| 日韩av免费高清视频| 蜜桃在线观看..| 最近中文字幕高清免费大全6| 成人亚洲精品一区在线观看| 在现免费观看毛片| 国产极品天堂在线| 18禁观看日本| xxx大片免费视频| 午夜福利一区二区在线看| 久久久久久久久免费视频了| 久久女婷五月综合色啪小说| 男人爽女人下面视频在线观看| 日韩人妻精品一区2区三区| 激情视频va一区二区三区| 久久精品久久久久久噜噜老黄| 精品人妻一区二区三区麻豆| 国产精品三级大全| 国产精品二区激情视频| 欧美日韩国产mv在线观看视频| 国产伦人伦偷精品视频| 秋霞在线观看毛片| 国产国语露脸激情在线看| 天天躁狠狠躁夜夜躁狠狠躁| 国产亚洲最大av| 十八禁网站网址无遮挡| 精品国产国语对白av| 久久99一区二区三区| 久久久久国产精品人妻一区二区| 99精品久久久久人妻精品| 国产亚洲精品第一综合不卡| 欧美少妇被猛烈插入视频| 国产探花极品一区二区| 日韩一卡2卡3卡4卡2021年| 涩涩av久久男人的天堂| 黑人猛操日本美女一级片| 国产毛片在线视频| netflix在线观看网站| 永久免费av网站大全| 国产精品一二三区在线看| 日韩伦理黄色片| 丁香六月天网| 成人影院久久| 国产成人精品在线电影| 老司机在亚洲福利影院| 观看av在线不卡| 亚洲成人免费av在线播放| 免费久久久久久久精品成人欧美视频| 99国产精品免费福利视频| 亚洲欧美一区二区三区久久| 日韩av不卡免费在线播放| 欧美成人精品欧美一级黄| 色吧在线观看| 亚洲欧美成人精品一区二区| 色播在线永久视频| a级毛片在线看网站| 晚上一个人看的免费电影| 一本—道久久a久久精品蜜桃钙片| 午夜福利一区二区在线看| 国精品久久久久久国模美| 国产伦人伦偷精品视频| 丰满饥渴人妻一区二区三| 亚洲av成人精品一二三区| 久久狼人影院| 啦啦啦中文免费视频观看日本| 一级毛片黄色毛片免费观看视频| 亚洲欧美激情在线| 香蕉国产在线看| 丝袜在线中文字幕| 纵有疾风起免费观看全集完整版| 最近中文字幕2019免费版| 满18在线观看网站| 亚洲成人手机| 国产在线免费精品| 亚洲,一卡二卡三卡| 捣出白浆h1v1| 黄色怎么调成土黄色| 欧美激情 高清一区二区三区| 久久久久人妻精品一区果冻| kizo精华| 久久热在线av| 一本一本久久a久久精品综合妖精| 国产精品女同一区二区软件| 国产成人一区二区在线| 黄片播放在线免费| 中文字幕人妻丝袜一区二区 | 国产日韩欧美视频二区| 国产精品国产三级专区第一集| 亚洲精品成人av观看孕妇| 亚洲人成电影观看| 丰满饥渴人妻一区二区三| 两性夫妻黄色片| 悠悠久久av| 久久这里只有精品19| 精品视频人人做人人爽| 一区二区三区四区激情视频| 亚洲精品aⅴ在线观看| 亚洲精品中文字幕在线视频| xxx大片免费视频| 精品少妇黑人巨大在线播放| 在现免费观看毛片| av福利片在线| 久久久精品免费免费高清| 男男h啪啪无遮挡| 色精品久久人妻99蜜桃| www.自偷自拍.com| 亚洲精品国产色婷婷电影| 夫妻性生交免费视频一级片| 国产亚洲av高清不卡| 性少妇av在线| 1024香蕉在线观看| 天天躁日日躁夜夜躁夜夜| 超色免费av| 无遮挡黄片免费观看| 人人澡人人妻人| 制服丝袜香蕉在线| 国产成人欧美| 国产男女内射视频| 人人妻人人澡人人爽人人夜夜| 少妇 在线观看| 性少妇av在线| 一本一本久久a久久精品综合妖精| 欧美人与性动交α欧美软件| 日本色播在线视频| 亚洲精品久久成人aⅴ小说| 午夜福利视频精品| 又大又爽又粗| tube8黄色片| 最近中文字幕高清免费大全6| 男女下面插进去视频免费观看| 狠狠精品人妻久久久久久综合| 免费av中文字幕在线| 国产男女超爽视频在线观看| 一区二区三区激情视频| 大码成人一级视频| www.av在线官网国产| 亚洲欧美精品综合一区二区三区| 国产成人91sexporn| 久久av网站| 在现免费观看毛片| 久久久久精品久久久久真实原创| 热99久久久久精品小说推荐| 亚洲国产最新在线播放| 亚洲美女黄色视频免费看| 男人添女人高潮全过程视频| 国产成人a∨麻豆精品| 十八禁人妻一区二区| 波多野结衣av一区二区av| 亚洲精品久久午夜乱码| 纵有疾风起免费观看全集完整版| 丝瓜视频免费看黄片| 亚洲国产精品一区三区| 亚洲av日韩在线播放| 国产欧美日韩一区二区三区在线| 亚洲国产av新网站| 日本欧美视频一区| 国产精品亚洲av一区麻豆 | 伊人久久国产一区二区| 国产97色在线日韩免费| 午夜福利视频在线观看免费| 成人国语在线视频| 波多野结衣av一区二区av| 黑人猛操日本美女一级片| 久久久久久久国产电影| 十分钟在线观看高清视频www| 久久精品aⅴ一区二区三区四区| 黑人欧美特级aaaaaa片| 日韩熟女老妇一区二区性免费视频| 久久精品国产综合久久久| 男女之事视频高清在线观看 | 两个人看的免费小视频| 久久久国产精品麻豆| 最近手机中文字幕大全| 久久免费观看电影| 如何舔出高潮| 日韩中文字幕视频在线看片| 国产一卡二卡三卡精品 | 97在线人人人人妻| 欧美中文综合在线视频| 99re6热这里在线精品视频| 日韩欧美精品免费久久| 亚洲精品aⅴ在线观看| 亚洲精品中文字幕在线视频|