張則強(qiáng) ，蔣 晉 ，尹 濤 ，許培玉

（1. 西南交通大學(xué)機(jī)械工程學(xué)院，四川 成都 610031；2. 西南交通大學(xué)軌道交通運(yùn)維技術(shù)與裝備四川省重點(diǎn)實(shí)驗(yàn)室，四川 成都 610031）

拆卸線是實(shí)現(xiàn)廢舊機(jī)電產(chǎn)品規(guī)?；妥詣?dòng)化生產(chǎn)的重要組織方式，而如何提高拆卸效率和生產(chǎn)線平衡率，引起學(xué)者們的廣泛關(guān)注，由此提出拆卸線平衡問題（disassembly line balancing problem，DLBP)[1].

關(guān)于DLBP問題的研究多集中于直線型DLBP問題優(yōu)化[1-4]. 在實(shí)際的生產(chǎn)作業(yè)中，相較于直線型布局方式，最大的不同之處在于U型拆卸線在設(shè)計(jì)過程中可從前后雙向搜索可分配作業(yè)任務(wù)，并可將首尾兩個(gè)作業(yè)任務(wù)分配到同一個(gè)工作站. 此外，U型布局具有生產(chǎn)柔性強(qiáng)、占地面積小、效率高等特點(diǎn)[5].肖欽心等[6]考慮U型布局下拆卸過程存在的多種約束限制，以最小化拆卸節(jié)拍時(shí)間和工作站空閑時(shí)間指標(biāo)為優(yōu)化目標(biāo)，采用改進(jìn)的并行鄰域搜索算法進(jìn)行求解. 張則強(qiáng)等[7]提出一種Pareto蟻群遺傳算法，對多目標(biāo)U型拆卸線平衡問題進(jìn)行協(xié)同優(yōu)化，以彌補(bǔ)傳統(tǒng)求解方法的不足. 拆卸作業(yè)過程中，由于部分零件存在需求性或危害性，需要后續(xù)相關(guān)處理操作，因此必須拆卸，其余零部件則可選擇性進(jìn)行拆卸，由此能節(jié)約成本，提高拆卸效率. Ren等[3]建立了以利潤為導(dǎo)向的不完全拆卸線問題模型，并通過引力搜索算法求解. Bentaha等[4]利用AND/OR優(yōu)先關(guān)系圖，通過基于拉格朗日松弛和蒙特卡洛采樣技術(shù)的求解方法，解決以利潤為最高要求的不完全拆卸線平衡問題. 因此，本文綜合考慮U型布局的柔性生產(chǎn)特點(diǎn)和不完全拆卸的高效性，提出U型不完全拆卸線平衡問題（U-shaped partial disassembly line balance problem，UPDLBP）.

對于DLBP問題的求解方法研究，Güng?r等[8]首次證明DLBP是NP-hard問題，其解空間隨著問題規(guī)模的增長呈指數(shù)增長，因此，數(shù)學(xué)規(guī)劃方法[9-10]不適用于求解大規(guī)模DLBP. 而啟發(fā)式算法[11-12]依賴啟發(fā)式規(guī)則，采用先決策后優(yōu)化的形式，將多目標(biāo)轉(zhuǎn)化為單目標(biāo)進(jìn)行計(jì)算，只能獲得一個(gè)容易受決策者主觀因素影響的解. 目前，群智能算法例如人工蜂群算法[13]、變鄰域搜索算法[14]、蟻群算法[15]等由于全局尋優(yōu)能力強(qiáng)和適用范圍廣等特點(diǎn)，在DLBP中得到廣泛應(yīng)用.

與其他算法相比，狼群算法（wolfpack algorithm，WPA)[16]受參數(shù)設(shè)置影響性較低，在求解速度和求解質(zhì)量上更具優(yōu)勢. 本文提出一種自適應(yīng)反向?qū)W習(xí)多目標(biāo)狼群算法（adaptive opposition-based learning multiobjective wolfpack algorithm，AOBL-MWPA），通過對算法進(jìn)行離散化操作，利用輪盤賭操作劃分狼群；采用自適應(yīng)游走行為，在迭代前期增強(qiáng)全局尋優(yōu)能力，后期增強(qiáng)局部搜索能力，加強(qiáng)探狼間信息交互，其余人工狼通過奔襲行為和召喚行為向當(dāng)前迭代次數(shù)的全局最優(yōu)解靠攏；引入反向?qū)W習(xí)策略（oppositionbased learning，OBL)[17]避免算法陷入局部最優(yōu)解；通過精英保留策略記錄全局最優(yōu)解，最后引入Pareto解集思想[18]和非支配排序遺傳算法Ⅱ（non-dominated sorting genetic algorithm Ⅱ，NSGA-Ⅱ）擁擠距離機(jī)制[19]獲得分布均勻的非劣解集.

1 U型不完全DLBP問題

1.1 問題描述

本文從實(shí)際角度出發(fā)，考慮采用U型布局和不完全拆卸方式，如圖1所示，工人同時(shí)對U型入口側(cè)和出口側(cè)零部件進(jìn)行拆卸，非必須拆卸的零部件直接進(jìn)行粉碎處理，降低拆卸成本，提高拆卸效率.對UPDLBP展開研究，并作出以下假定：1） 待拆卸產(chǎn)品單一，且能保證長期供應(yīng)不中斷；2） 零部件完整，且拆卸作業(yè)時(shí)間為標(biāo)準(zhǔn)定值；3） 忽略任何突發(fā)情況；4） 為簡化問題，不考慮工人行走、物料運(yùn)輸?shù)绕渌蛩赜绊?

圖1 U型不完全拆卸線Fig. 1 U-shaped partial disassembly line

1.2 符號說明

i,j：任務(wù)編號，i,j= 1,2,···,N；

k：工作站編號，k= 1,2,···,K；

TC：拆卸節(jié)拍時(shí)間；

Tk：第k個(gè)工作站的實(shí)際工作時(shí)間；

ti：任務(wù)i作業(yè)時(shí)間；

hi：若任務(wù)i具有危害性，hi=1 ，否則hi=0；

di：若任務(wù)i具有需求性，di=1 ，否則di=0；

wk：若工作站k的拆卸任務(wù)集中存在危害性零部件，wk=1 ，否則wk=0；

Pij：若任務(wù)i為任務(wù)j的緊前任務(wù)，Pij=1，否則Pij=0. 優(yōu)先關(guān)系矩陣P= (Pij)N×N；

Ci：任務(wù)i單位時(shí)間作業(yè)成本；

Ch：工作站單位時(shí)間附加無害化處理成本；

Cs：工作站單位時(shí)間待機(jī)成本；

xik：決策變量，若任務(wù)i分配到工作站k的入口側(cè)，xik=1，否則xik= 0 ；

yik：決策變量，若任務(wù)i分配到工作站k的出口側(cè)，yik=1，否則yik= 0；

Sk：決策變量，若開啟工作站k，Sk=1，否則Sk=0 .

1.3 優(yōu)化目標(biāo)

實(shí)際DLBP問題中，決策者希望同時(shí)優(yōu)化多個(gè)目標(biāo)，既能保證企業(yè)的經(jīng)濟(jì)效益，降低企業(yè)成本，還能保證拆卸線的平衡性. 因此，本文考慮協(xié)同優(yōu)化以下4個(gè)目標(biāo)：最小化工作站數(shù)量f1、最小化空閑時(shí)間均衡指標(biāo)f2、最小化拆卸深度f3及最小化拆卸成本f4，分別如式（1） ～ （4）所示，其中，拆卸成本由拆卸作業(yè)成本、工作站待機(jī)空閑成本及附加無害化處理成本三部分組成.

1.4 約束條件

約束條件包括節(jié)拍時(shí)間約束，任務(wù)分配約束，優(yōu)先關(guān)系約束，需求、危害指標(biāo)約束和工作約束.

節(jié)拍時(shí)間約束為

由于采用不完全拆卸方式，允許存在部分零部件不拆卸，但拆卸任務(wù)有且僅能分配到一個(gè)工作站.任務(wù)分配約束如式（6）所示.

采用U型布局形式，待拆卸零部件可在入口側(cè)或出口側(cè)進(jìn)行拆卸，分別加以優(yōu)先關(guān)系約束. 優(yōu)先關(guān)系約束如式（7）所示.

具有危害屬性和經(jīng)濟(jì)效益的零部件必須拆卸.需求、危害指標(biāo)約束如式（8）所示.

需要開啟的工作站存在理論上下閾值約束，且不允許開啟空的工作站內(nèi)無任務(wù)分配. 工位約束如式（9）、（10）所示.

2 AOBL-MWPA算法設(shè)計(jì)

本文結(jié)合DLBP實(shí)際問題特點(diǎn)，提出一種自適應(yīng)反向?qū)W習(xí)多目標(biāo)狼群算法進(jìn)行求解計(jì)算，由于DLBP問題中需要優(yōu)化多個(gè)存在一定制衡性的目標(biāo)，很難同時(shí)達(dá)到最優(yōu)值，因此考慮去除頭狼作用，保留其他算法機(jī)制. 采用實(shí)數(shù)編碼方式生成可行拆卸序列，通過對算法操作進(jìn)行離散化操作，并引入反向?qū)W習(xí)策略跳出局部最優(yōu)，通過Pareto解集和NSGA-Ⅱ擁擠距離機(jī)制評價(jià)非劣解集，引入精英保留策略加快算法向全局最優(yōu)解靠攏.

2.1 編碼及解碼設(shè)計(jì)

在拆卸線設(shè)計(jì)過程中通過人工識別確定零部件的組成結(jié)構(gòu)、拆卸作業(yè)時(shí)間、危害和需求屬性等拆卸信息，并整合形成優(yōu)先關(guān)系矩陣P和拆卸信息矩陣B，以便作為算法求解計(jì)算所能識別處理的問題參數(shù)輸入，在滿足特定約束條件下通過編碼操作生成可行解，通過解碼操作將拆卸任務(wù)合理分配到各工作站內(nèi)并計(jì)算目標(biāo)適應(yīng)度值，為實(shí)現(xiàn)廢舊產(chǎn)品的大規(guī)模拆卸奠定數(shù)據(jù)基礎(chǔ).

采用基于實(shí)數(shù)編碼方式生成可行拆卸序列，將優(yōu)先關(guān)系矩陣P中無緊前約束任務(wù)放入可選任務(wù)集中，隨機(jī)挑選其中的一個(gè)任務(wù)i放入當(dāng)前拆卸序列位置m中，將P中任務(wù)i所在行全部置為0，以釋放對其他任務(wù)的約束關(guān)系，同時(shí)將任務(wù)i所在列全部置為1，以加強(qiáng)自身約束，保證不會重復(fù)選取到該任務(wù)，然后m=m+1，重復(fù)以上步驟，直至所有任務(wù)均分配完成，由此得到的實(shí)數(shù)序列即為可行拆卸序列.

由于在U型不完全拆卸過程中，并不是所有任務(wù)都需要進(jìn)行拆卸，只要零部件存在需求性或危害性這兩種特性中的其中一個(gè)，就必須要拆卸該零部件，但也會存在某一個(gè)零部件同時(shí)存在需求性與危害性的情況，在解碼時(shí)僅需確保必須拆卸零部件總數(shù)滿足要求即可，不會影響算法性能. 因此，將解碼過程分為兩個(gè)階段，第一階段是得到實(shí)際需要拆卸的任務(wù)序列，第二階段是將所得的拆卸序列在滿足所有約束情況下分配到工作站中，具體步驟如下：

步驟1輸入可行拆卸序列X，計(jì)算拆卸產(chǎn)品中必須拆除的零部件總數(shù)Nnum，設(shè)置開啟工作站序號ks=1，工作站空閑時(shí)間TR=TC，危害性或需求性零部件計(jì)數(shù)iex=1；

步驟2從解序列位置編號m= 1開始遍歷X，若某位置m上的任務(wù)i∈{i|di+hi≥1}，則iex=iex+1；若iex=Nnum，則將對應(yīng)位置編號m及之前的序列單獨(dú)取出，作為新的不完全拆卸序列Yp輸出，序列長度為l，解碼第一階段完成；

步驟3初始化序列位置編號，入口位置p=1，出口位置q=l；

步驟4判斷位置p、q上的任務(wù)i、j的作業(yè)時(shí)間是否小于當(dāng)前工作站空閑時(shí)間TR，若條件成立，則將滿足條件的任務(wù)放入可分配任務(wù)集合S中；

步驟5確定可分配任務(wù)集合大小NS,num，若NS,num=0 ，則開啟新的工作站ks=ks+1 ，重置空閑時(shí)間TR=TC；否則將拆卸作業(yè)時(shí)間較長的零部件作為待分配任務(wù)；

步驟6確定待分配任務(wù)的編號，若為i，則分配到入口側(cè)，p=p+1 ；否則，分配到出口側(cè)，q=q?1； 修改工作站剩余空閑時(shí)間TR=TR?ti；

步驟7重復(fù)步驟4～6，直至序列Yp中任務(wù)分配完畢.

2.2 “輪盤賭”選擇操作

本文采用“輪盤賭”法進(jìn)行種群劃分，根據(jù)個(gè)體適應(yīng)度值計(jì)算個(gè)體選擇概率及累加概率選取As匹人工狼組成探狼群，具體計(jì)算如下：

式中：NW為種群規(guī)模；f(Xi)、p(Xi)分別為第i個(gè)個(gè)體Xi的適應(yīng)度值和選擇概率；P(Xi)為累加概率，其值越大，個(gè)體選中機(jī)率越大，由此選擇探狼個(gè)體.

2.3 自適應(yīng)游走行為

探狼在解空間內(nèi)朝R個(gè)方向分別進(jìn)行游走搜索獵物，根據(jù)獵物氣味濃度（可認(rèn)定為目標(biāo)適應(yīng)度值）自主決策前進(jìn)方向，具體表示為

式中：xie,r為探狼i在第e（e= 1,2,···,E）維空間朝第r（r= 1,2,···,R）個(gè)方向前進(jìn)后的位置；xie為探狼i在第e維空間的位置；sa為游走步長.

為了使算法在迭代前期擴(kuò)大全局搜索尋優(yōu)能力，同時(shí)在迭代后期有更大概率朝全局最優(yōu)解靠攏，加強(qiáng)穩(wěn)定性. 考慮隨著迭代次數(shù)的增加，減少游走方向數(shù)R，加快算法收斂速度. 計(jì)算方法如下：

式中：Rmax、Rmin分別為游走方向數(shù)的最大、最小值；g為當(dāng)前迭代次數(shù)；Mgen為算法最大迭代次數(shù).

結(jié)合DLBP編碼特點(diǎn)，設(shè)計(jì)一種基于隨機(jī)擾動(dòng)的自適應(yīng)游走行為，產(chǎn)生若干鄰域解，具體步驟如下：

步驟1輸入拆卸序列X，設(shè)置游走次數(shù)kt= 1，游走方向計(jì)數(shù)r= 1，游走次數(shù)閾值為kt,max，根據(jù)式（13）計(jì)算游走方向數(shù)R；

步驟2生成表示序列位置的隨機(jī)整數(shù)u（u∈{1,2,···,U})，確定u上的任務(wù)編號iask；

步驟3利用優(yōu)先關(guān)系矩陣確定iask相鄰的緊前任務(wù)Ptask和緊后任務(wù)Stask，將緊前緊后任務(wù)間的序列作為可操作序列片段Yo；

步驟4生成與序列Y相同長度的單調(diào)遞增隨機(jī)數(shù)組，將iask對應(yīng)的隨機(jī)數(shù)rnum按式（12）進(jìn)行計(jì)算，并將隨機(jī)數(shù)組按遞增順序重排，對應(yīng)的iask位置發(fā)生變動(dòng)，生成新的序列片段Z；

步驟5用序列片段Z代替Yo，放入序列X中，生成新的鄰域解XNew并保存，更新游走方向計(jì)數(shù)r=r+1；

步驟6若r

步驟7若kt

2.4 召喚行為

結(jié)合DLBP實(shí)際特點(diǎn)，利用基于交叉操作接收召喚信息，并通過生成隨機(jī)整數(shù)確定變異點(diǎn)，基于變異操作執(zhí)行猛狼向頭狼奔襲行為. 具體步驟如下：

步驟1輸入探狼序列X，猛狼序列Y，生成兩個(gè)表示序列位置的隨機(jī)整數(shù)p、q；

步驟2確定序列Y中位置p、q間的序列片段Z，并遍歷序列X，若某位置上的任務(wù)編號在序列Z中存在，則依次存放到新的序列片段ZNew中；

步驟3用ZNew替換Y中原拆卸序列片段Z，生成新的拆卸序列YNew1，完成召喚信息傳遞；

步驟4生成表示拆卸序列位置的隨機(jī)整數(shù)m，確定YNew1中該位置的任務(wù)編號iask；

步驟5根據(jù)優(yōu)先關(guān)系矩陣確定iask在YNew1中的相鄰最近的緊前緊后任務(wù)位置a、b，將位置a和b之前的序列除去iask作為可操作序列，在其中隨機(jī)插入iask，由此生成新的拆卸序列YNew2；

步驟6比較拆卸序列Y、YNew1及YNew2的適應(yīng)度值，更新猛狼位置.

2.5 目標(biāo)驅(qū)動(dòng)圍攻行為

為快速收斂，本文考慮將單個(gè)目標(biāo)最優(yōu)的人工狼作為驅(qū)動(dòng)解按式（14）執(zhí)行圍攻行為，為算法每次的迭代尋優(yōu)提供單目標(biāo)參考值.

式中：λ為 [?1,1]間的隨機(jī)數(shù)；sc為圍攻步長；和分別為第k次迭代時(shí)人工狼與獵物在第e維空間所處位置.

步驟1計(jì)算交換序列對數(shù)ENum，生成ENum個(gè)不大于拆卸序列長度N的互不相等的隨機(jī)整數(shù)，用以確定交換對位置集EPos；

步驟2根據(jù)EPos逐個(gè)確定G、xi對應(yīng)位置上的任務(wù)交換對{GTask,xi,Task}；

步驟3基于DLBP實(shí)數(shù)編碼的特點(diǎn)，同一任務(wù)有且僅能分配一次，因此確定GTask在拆卸序列xi中的位置，并交換序列對GTask與xi,Task的位置；

步驟4判斷交換序列對位置后是否滿足優(yōu)先關(guān)系約束，條件成立的交換對如圖2中實(shí)線表示的交換對{5,4}、{3,2}，將其予以保留，將不滿足條件的如圖2中虛線表示的交換對{10,7}舍棄.

圖2 目標(biāo)驅(qū)動(dòng)圍攻行為Fig. 2 Goal-driven besieging behavior

2.6 反向?qū)W習(xí)策略

反向?qū)W習(xí)策略[17]基本思想是為擴(kuò)大種群搜索范圍，基于當(dāng)前個(gè)體在可行域范圍內(nèi)生成一個(gè)對應(yīng)的反向解個(gè)體，比較兩個(gè)個(gè)體的適應(yīng)度值，選擇較優(yōu)的個(gè)體進(jìn)入下一次迭代優(yōu)化. 具體定義如下：設(shè)解Xi={xi1,xi2,···,xie}為e維空間中的一個(gè)可行解，xie為第e維空間的變量，其定義域?yàn)?[ae,be]，

式中：rn為 [0,1] 間的隨機(jī)數(shù).

在滿足DLBP優(yōu)先關(guān)系約束情況下，生成人工狼反向種群，從中挑選適應(yīng)度值較優(yōu)的人工狼替換原種群中較差的個(gè)體，保證種群的優(yōu)越性.

2.7 多目標(biāo)處理方法

由于DLBP問題涉及多個(gè)量綱不同的目標(biāo)優(yōu)化，不能簡單的將其轉(zhuǎn)化為單目標(biāo)問題求解，同時(shí)缺少決策者的個(gè)人偏好等關(guān)鍵信息，需要采取先優(yōu)化再決策的方式，因此考慮引入Pareto解集理論獲得多個(gè)具有不同側(cè)重點(diǎn)的非劣解.

對于最小化多目標(biāo)優(yōu)化問題，給定兩個(gè)可行解X1、X2，其第w個(gè)目標(biāo)適應(yīng)度值分別為fw(X1)，fw(X2)，若滿足式（16）則稱解X1Pareto支配解X2，記為X1?X2.

本文考慮采用基于外部檔案的精英保留策略，將每次迭代尋優(yōu)獲得的Pareto解存放于外部檔案中，并不斷更新外部檔案，使算法加速朝全局最優(yōu)解收斂. 當(dāng)非劣解個(gè)數(shù)超過設(shè)置的外部檔案容量時(shí)，通過NSGA-Ⅱ擁擠距離機(jī)制[19]篩選去除部分非劣解.定義邊界個(gè)體的擁擠距離為∞，其余解X對應(yīng)的單目標(biāo)值fw(X)分別按升序排列后，計(jì)算如式（17）所示.

式中：L(X)為解X的擁擠距離；Xg?1、Xg+1分別為X排序前、后的解.

2.8 算法流程

AOBL-MWPA具體算法流程如圖3所示.

圖3 BL-MWPA算法流程Fig. 3 AOBL-MWPA flowchart

3 算法性能驗(yàn)證和實(shí)例應(yīng)用

在硬件配置為Intel(R) Core(TM) i5-9500 CPU，3.00 GHz主頻，8.00 GB內(nèi)存的計(jì)算機(jī)上，通過Win10系統(tǒng)上的MATLAB 2018開發(fā)算法程序. 為驗(yàn)證所以算法的性能，將其應(yīng)用于不同規(guī)模的算例中求解計(jì)算，并與現(xiàn)有文獻(xiàn)中的算法結(jié)果進(jìn)行對比.

3.1 基準(zhǔn)算例驗(yàn)證

通過文獻(xiàn)[20]的19個(gè)基準(zhǔn)算例來驗(yàn)證AOBLMWPA算法的性能，定義任務(wù)拆卸數(shù)量N為首項(xiàng)為8，公差為4，末項(xiàng)為80的一組等比數(shù)列，不考慮零部件間的優(yōu)先關(guān)系約束，任務(wù)拆卸作業(yè)時(shí)間ti∈{3,5,7,11} ，節(jié)拍時(shí)間TC= 26 s，定義最后一個(gè)作業(yè)時(shí)間為11 s的任務(wù)為危害性零部件，最后一個(gè)作業(yè)時(shí)間為7 s的任務(wù)為需求零部件. 優(yōu)化目標(biāo)有：最小化工作站數(shù)量f1、最小化空閑時(shí)間均衡指標(biāo)f2、最小化拆卸深度f3及最小化拆卸成本f4，同時(shí)附加計(jì)算求解時(shí)間進(jìn)行對比說明. 現(xiàn)有文獻(xiàn)中的求解算法有蟻群算法（ant colony optimization，ACO）[20]、改進(jìn)蟻 群算法（improved ant colony optimization，IACO）[15]、多目標(biāo)免疫機(jī)制協(xié)作遺傳算法（multi-objective immune mechanism cooperative genetic algorithm，MIGA）[21]及變鄰域搜索算法（variable neighborhood search，VNS）[14]. 通過正交試驗(yàn)設(shè)計(jì)確定算法參數(shù)，因篇幅問題不展開說明，從求解時(shí)間和求解質(zhì)量兩方面綜合考慮，設(shè)置算法參數(shù)為NW= 60，Mgen=120，sa=1.5，kt,max= 8，計(jì)算對比結(jié)果如圖4所示.

圖4 19個(gè)基準(zhǔn)算例對比結(jié)果Fig. 4 Comparison results of 19 benchmark instances

由基準(zhǔn)算例的構(gòu)造方法可得，真實(shí)Pareto前沿解最優(yōu)結(jié)果為{N/4,0,1,2}或{N/4,0,2,1}. 通過對所提算法和現(xiàn)有文獻(xiàn)求解結(jié)果的均值進(jìn)行對比，在工作站數(shù)量f1和空閑時(shí)間均衡指標(biāo)f2上，所提AOBLMWPA算法和MIGA算法在小規(guī)模（N≤20）上均能求解到最優(yōu)解f1=N/4，f2= 0，但在中規(guī)模（20

3.2 實(shí)例應(yīng)用

現(xiàn)以某型報(bào)廢汽車為對象，分析AOBL-MWPA算法在U型不完全拆卸線方案規(guī)劃中的應(yīng)用情況.現(xiàn)獲得的拆卸信息包括零部件拆卸作業(yè)時(shí)間t（s），危害屬性h，需求屬性d，單位時(shí)間作業(yè)成本C（元/s），具體如表1所示，優(yōu)先關(guān)系如圖5所示.

圖5 某汽車拆卸優(yōu)先關(guān)系Fig. 5 Disassembly precedence relations of a car

表1 某汽車拆卸信息Tab. 1 Disassembly information of a car

根據(jù)實(shí)際生產(chǎn)情況，設(shè)置其生產(chǎn)節(jié)拍TC= 640 s.經(jīng)過正交試驗(yàn)測試計(jì)算，綜合考慮算法求解精度和求解時(shí)間，確定算法參數(shù)為NW= 60，Mgen= 120，sa=1.5，Tmax= 8，外部檔案大小NQ= 10. 算法運(yùn)行求解10次，取其中較好的一次結(jié)果如表2所示，拆卸方案中帶負(fù)號的編號零件表示該零件在U型出口側(cè)進(jìn)行拆卸作業(yè)，反之，則在U型入口側(cè)進(jìn)行拆卸作業(yè). 通用工作站需要對入口側(cè)與出口側(cè)的產(chǎn)品進(jìn)行拆卸作業(yè)，專用工作站只需對流水線某一側(cè)待拆產(chǎn)品進(jìn)行作業(yè).

從表2結(jié)果看出，根據(jù)Pareto解集思想，所提算法能求解得到10個(gè)非劣解，開啟的工作站數(shù)量均為9個(gè)，拆卸零部件數(shù)量為38個(gè)，空閑時(shí)間均衡指標(biāo)變化范圍為28.142 5～33.015 1 s，拆卸成本變化范圍為381.946 4～387.954 9 元/s，因此決策者可根據(jù)實(shí)際情況選擇合適的拆卸方案，若決策者注重空閑時(shí)間均衡指標(biāo)可選用方案2，該優(yōu)化目標(biāo)為28.142 5 s；若決策者注重企業(yè)的拆卸成本，則可選擇方案1，最小成本為381.946 4元/s.

表2 U型汽車拆卸線任務(wù)分配方案Tab. 2 Task allocation plan for U-shaped car disassembly line

4 結(jié) 論

1）本文結(jié)合企業(yè)實(shí)際情況，對U型不完全拆卸線平衡問題展開研究，構(gòu)建了以最小化工作站數(shù)量，空閑時(shí)間均衡指標(biāo)，拆卸深度和拆卸成本為優(yōu)化目標(biāo)的數(shù)學(xué)模型.

2）針對問題特點(diǎn)，設(shè)計(jì)一種自適應(yīng)反向?qū)W習(xí)多目標(biāo)狼群算法，為兼顧算法前期的全局搜索能力和后期的穩(wěn)定性，提出改進(jìn)的自適應(yīng)游走行為，同時(shí)引入反向精英學(xué)習(xí)策略有效避免算法陷入局部最優(yōu)解，利用Pareto解集思想和NSGA-Ⅱ擁擠距離篩選評價(jià)機(jī)制，實(shí)現(xiàn)不同側(cè)重點(diǎn)的精英解集的保留，指導(dǎo)算法搜索尋優(yōu)方向，同時(shí)加快算法收斂速度.

3）通過將所提算法應(yīng)用于不同規(guī)模問題的19個(gè)基準(zhǔn)算例，并與現(xiàn)有文獻(xiàn)中的求解算法進(jìn)行對比，結(jié)果表明所提算法在不同評價(jià)指標(biāo)上的求解結(jié)果優(yōu)于其他算法，驗(yàn)證了其求解性能更優(yōu). 最后，將所提算法應(yīng)用于采用U型布局不完全拆卸方式的具有40項(xiàng)任務(wù)的某汽車拆卸實(shí)例中，求解獲得10個(gè)Pareto非劣解，為決策者提供了廣泛的決策空間，驗(yàn)證了所提算法和模型的正確性和有效性.

致謝：中車“十四五”科技重大專項(xiàng)課題（2021CHZ010-3）.