楊 磊，徐曉鳴※，張海洋，曾偉鎮(zhèn)

（1.廣東海洋大學機械與動力工程學院，廣東湛江 524088；2.小天才科技有限公司，廣東東莞 523846）

0 引言

在需求多樣化市場下，客戶對產(chǎn)品的交付質(zhì)量和交付周期要求越來越高，如何將客戶需求與批量生產(chǎn)相結(jié)合來組織生產(chǎn)，是企業(yè)面臨的一個難題。國內(nèi)外企業(yè)紛紛引進精益生產(chǎn)方式，采取混流生產(chǎn)線，以期提高生產(chǎn)效率和降低節(jié)拍，提高企業(yè)自身的競爭力[1-3]。然而，進行混流生產(chǎn)線改造首先要解決其平衡問題。

針對混流生產(chǎn)平衡研究中，根據(jù)研究內(nèi)容不同，第一類是給定節(jié)拍時間，求最小工位數(shù)；第二類是給定生產(chǎn)線長度，求最小化節(jié)拍時間；第三類是優(yōu)化生產(chǎn)線工位負荷均衡。以往的研究大致有以下幾種研究方法：工業(yè)工程改善方法、啟發(fā)式算法、仿真方法和智能算法。吳智勇等[2]通過建立補充生產(chǎn)機制、模式排產(chǎn)、作業(yè)改善等精益改善方案，提升了混流生產(chǎn)的問題；許向川等[3]利用Flexsim 仿真的方法，建立仿真模型，通過精益改善，提升了生產(chǎn)線平衡率；周金平[4]利用Plant Simulation 仿真軟件研究混流生產(chǎn)的物料配送問題，優(yōu)化最大完工時間。包建軍等[5]運用啟發(fā)式算法對混流裝配線的平衡問題進行分析與優(yōu)化，并采用Flexsim 仿真驗證優(yōu)化方案的有效性。邱伊健等[6]利用Flexsim 對混流生產(chǎn)線的瓶頸問題，生產(chǎn)線平衡問題進行優(yōu)化，并進一步利用遺傳算法優(yōu)化投產(chǎn)序列；李偉等[7]針對混流裝配線第二類平衡問題，建立了基于作業(yè)空間和裝配工具約束等條件下的優(yōu)化模型，并利用遺傳算法進行求解；鄭晨鳴等[8]考慮了作業(yè)時間隨機情形下，最小工作站數(shù)和各產(chǎn)品負荷均勻為目標，設計改進的遺傳算法并驗證其有效性。

本文結(jié)合A 企業(yè)實際生產(chǎn)場景，通過合并工藝相似的產(chǎn)品生產(chǎn)線，在考慮工序之間優(yōu)先順序、產(chǎn)品需求信息和作業(yè)時間約束下，以生產(chǎn)節(jié)拍和負荷均衡為目標，建立混流生產(chǎn)線第二類平衡問題的多目標混合整數(shù)規(guī)劃模型，設計求解模型的有效算法，有效提高混流生產(chǎn)線平衡率和降低生產(chǎn)線均衡指數(shù)，從而提高企業(yè)的生產(chǎn)效率。研究內(nèi)容為企業(yè)改造混流生產(chǎn)線提供理論依據(jù)和技術支持，可以推廣到其他類型的混流生產(chǎn)線設計與改造上。

1 混流生產(chǎn)線平衡模型建立

在一個計劃期T內(nèi)，已知某混流生產(chǎn)線上有N個工作站完成M種產(chǎn)品的裝配作業(yè)，這些產(chǎn)品的工藝相近、結(jié)構(gòu)相似，屬于企業(yè)內(nèi)同一產(chǎn)品族的[9-10]。該混流生產(chǎn)線上所有類型產(chǎn)品所包含的工序數(shù)為S，每種產(chǎn)品需求為Dm，M種產(chǎn)品的總需求量為表示第m種產(chǎn)品在第i道工序上加工的時間，由于不同類型的產(chǎn)品工序上有所不同，所以如果第m種產(chǎn)品不含有第i道工序，那么tim= 0，qm為第m種產(chǎn)品占總需求的比例，其中qm=Dq∕D；Tmk為第m種產(chǎn)品在第k個工作站單件裝配的時間。

模型假設條件：（1）各個作業(yè)元素的操作時間不變；（2）混流生產(chǎn)線的工作站可以進行一種或者多種作業(yè)元素的操作，員工均為多技能工，可以完成任何任務；（3）一個作業(yè)元素僅能被分配到工作站，一個作業(yè)元素只能在一個工作站內(nèi)完成該作業(yè)元素的全部內(nèi)容；（4）產(chǎn)品作業(yè)元素的優(yōu)先關系已知，并且作業(yè)元素加工的優(yōu)先關系不可變動。

該問題數(shù)學模型如下：

式中：CT為混流生產(chǎn)線的實際生產(chǎn)節(jié)拍；C*為混流生產(chǎn)線的理論節(jié)拍，C*=T∕D；Pre(i)為依據(jù)作業(yè)元素的優(yōu)先關系，先于作業(yè)元素i完成的作業(yè)元素集合。

式（1）表示在相同權重下求取最小生產(chǎn)節(jié)拍和負荷指數(shù)；式（2）表示產(chǎn)品在生產(chǎn)線中最小生產(chǎn)節(jié)拍；式（3）表示生產(chǎn)線的負荷指數(shù)，式（4）表示每個作業(yè)元素都只能被分配到一個工作站進行加工；式（5）作業(yè)元素的分配必須遵守優(yōu)先關系集合；式（6）表示工作站的作業(yè)時間不得大于理論節(jié)拍時間；式（7）表示作業(yè)元素分配的決策變量，1表示作業(yè)元素i被分配到第k工作站到，0則表示作業(yè)元素i未被分配到k工作站。

2 雙種群遺傳算法

傳統(tǒng)的遺傳算法中，只有一個初始種群開始操作、交叉和編譯遺傳，也稱為單種群的遺傳算法。雖然這種算法簡單快捷，但是它很容易陷入局部最優(yōu)解而且搜索空間有限，因此需要在這種算法的基礎上擴大搜索空間，能夠跳出局部最優(yōu)解，所以需要一種改進的遺傳算法要在不止一個的種群之中來開展操作[11-12]。

雙種群遺傳算法思路是兩個子種群使其各自進化，再選擇合適的節(jié)點，使其按照設定的規(guī)則進行種群間的交流，這種方式使用起來并不算復雜，而且可以擴大了搜索的空間而避免了陷入局部最優(yōu)的缺點。其中這兩個子種群都有它各自的作用，種群1 中的交叉概率Pc和變異概率Pm的值較大，這個種群稱為探測子種群，主要是在進化過程不斷地超越平面，克服過早收斂，使算法可以保持較高的搜索效率；種群2 的交叉概率Pc和變異概率Pm的值較小，被稱為開發(fā)子種群，主要可以在局部范圍內(nèi)尋找優(yōu)秀個體，并把這些個體保持下來。

2.1 編碼設計

編碼直接決定了遺傳算法后續(xù)的操作，如選擇算子、交叉算子和變異算子的實現(xiàn)，遺傳算法的計算性能是受編碼直接影響的。應對不同的實際問題，應采取恰當?shù)木幋a方式，這樣才不會產(chǎn)生大量的非可行解，使得計算過程變得十分冗余。通過優(yōu)先關系圖來設計編碼方式，以保證種群解的可行性，編碼設計思路如圖1所示。

圖1 優(yōu)先作業(yè)關系編碼

2.2 譯碼設計

譯碼的思路如下：對未知的生產(chǎn)節(jié)拍，進行節(jié)拍預估，一般是預估一個理論節(jié)拍CT*，然后以一定的增量去進行試探計算。

步驟1：初始化節(jié)拍，計算理論節(jié)拍令CT*等于這個理論節(jié)拍。

步驟2：以CT*為暫定節(jié)拍，并將作業(yè)元素按照作業(yè)元素優(yōu)先關系圖的邏輯關系將m個作業(yè)元素分配到n個工作站中?？梢缘玫矫總€工作站的時間：T1，T2，…，Tk，…，Tn，如果這些工作站的時間小于CT*，那么CT*就是這個狀態(tài)下的最小節(jié)拍，搜索就停止，如果不是那么進行下一步計算。

步驟3：計算潛在增量Δ1,Δ2,Δ3,…,Δn，其中Δi是第i+1 個工作站的第一作業(yè)元素的時間(i= 1,2,3,…,k,…,n)。令目前的工作站中作業(yè)時間最大 的 作 為 實 際 節(jié) 拍CT，CT*= min{Tk+ Δi} ，如 果CT≤CT*那么CT就是這種分配方案下的最小節(jié)拍，搜索停止；如果不符合條件就返回步驟2繼續(xù)分配作業(yè)元素。

2.3 適應度函數(shù)設計

本生產(chǎn)線平衡建立模型第一目標是最小生產(chǎn)節(jié)拍，本文通過平衡率指標來映射生產(chǎn)節(jié)拍，第二目標函數(shù)則是負荷指數(shù)，所以構(gòu)建的適應度函數(shù)為：

2.4 選擇算子設計

采取輪盤賭方式進行選擇，將種群中的個體分別計算遺傳到下一代的概率，再計算每個個體的累計概率，當個體適應度越大被選中的概率也就越大。

2.5 交叉算子設計

交叉算子采取的是兩點三段交叉，首先在父代染色體選取兩個交叉點，將這交叉點兩段之間的基因截取，在父代2中尋找相同基因片段的排序，將父代2中排序的基因投放至父代1 的交叉點之間形成子代1。子代2 的形成，同理可得。如圖2所示。

圖2 交叉示例

2.6 變異算子設計

采用單點交叉，在父代染色體中獲取變異點，變異點之前的基因保持不變，變異點之后的基因根據(jù)刪除后的優(yōu)先關系矩陣重新編排。如圖3所示。

圖3 變異示例

3 案例應用

3.1 現(xiàn)有門體預裝生產(chǎn)線的平衡分析

目前A 公司的門體預裝生產(chǎn)線的生產(chǎn)模式是可變品種的批量生產(chǎn)，取該生產(chǎn)線生產(chǎn)量最大的3 種產(chǎn)品BD330、BD144、BD96，其日需求量分別為500、500、250，A 公司的日有效工作時間為10 h，各種型號的工序加工時間如表1所示。

表1 工序加工時間

BD330、BD144、BD96的理論節(jié)拍分別為72 s、72 s、144 s，生產(chǎn)平衡率分別為0.531、0.373 和0.509，現(xiàn)場測得BD330、BD144、BD96 的實際生產(chǎn)節(jié)拍分別為77.35 s、77.35 s和60.26 s，僅BD96的實際生產(chǎn)節(jié)拍在計劃節(jié)拍內(nèi)，在計劃期內(nèi)只有BD96 可以完成預期的計劃。剩下兩種產(chǎn)品都只能通過加班加點的方式來進行產(chǎn)品產(chǎn)量的補償。

3.2 門體的作業(yè)元素劃分

進行混流生產(chǎn)線平衡的研究要將在生產(chǎn)線上生產(chǎn)的產(chǎn)品最小作業(yè)單元劃分出來，根據(jù)A 公司門體預裝生產(chǎn)線的實際生產(chǎn)情況，按照作業(yè)元素劃分的原則，使每個作業(yè)工序劃分到不能進行下一步劃分的程度劃分作業(yè)元素時，遵循的原則如下。

（1）工序在被拆分之前可以進行標準時間測量。

（2）若遇到不能再進一步拆解的工序，那么該工序就是最小的作業(yè)單元。

（3）某些工序存在特殊的加工方法，在拆分工序后也必須保持加工方法。

考慮實際的生產(chǎn)情況和劃分作業(yè)元素遵循的原則，將作業(yè)元素劃分如表2所示。

表2 作業(yè)元素劃分表

3.3 繪制出綜合優(yōu)先關系圖

過各個產(chǎn)品的作業(yè)元素優(yōu)先加工順序表以及對應的加工時間，再通過各產(chǎn)品的作業(yè)元素優(yōu)先加工順序表和各產(chǎn)品之間的需求比例制作出混合作業(yè)元素的時間表，進而繪制出綜合優(yōu)先關系圖，如圖4所示。

圖4 綜合優(yōu)先關系

3.4 混流生產(chǎn)線改造

根據(jù)綜合優(yōu)先關系圖，按照關系圖中的作業(yè)元素的加工順序邏輯，一個作業(yè)元素是另外一個作業(yè)元素的緊前工序則標注1，若不是緊前工序則標注0，按照這樣的規(guī)則建立綜合優(yōu)先關系矩陣。

基于平衡模型和雙種群遺傳算法，matlab 編程求解，將現(xiàn)有的門體加工作業(yè)元素進行重新分配，將各作業(yè)元素分配至固定的工作站中。

作業(yè)元素分配的結(jié)果如下：

工作站1：作業(yè)元素1→2→4

工作站2：作業(yè)元素3→5

工作站3：作業(yè)元素6

工作站4：作業(yè)元素7

工作站5：作業(yè)元素8→10→13

工作站6：作業(yè)元素9→14→12→11

工作站7：作業(yè)元素15→16

工作站8：作業(yè)元素17

3.5 改造效果

通過混流改造后，計算可得混流生產(chǎn)線的平衡率為0.859、節(jié)拍為48.16 s，平滑指數(shù)為3.375。進行混流改造后實際生產(chǎn)節(jié)拍遠低于各個型號產(chǎn)品單獨生產(chǎn)時的實際生產(chǎn)節(jié)拍，平衡率也相比于之前各產(chǎn)品單獨生產(chǎn)時有了較大幅度的提升。

4 結(jié)束語

本文通過對A 公司生產(chǎn)線研究，對其生產(chǎn)線進行混流改造，通過建立生產(chǎn)線平衡問題，應用雙種群遺傳算法對其進行工序的重新編排，求解出較優(yōu)的工序安排解，改造后優(yōu)化了生產(chǎn)效率，提高了生產(chǎn)效率，

研究表明針對多樣化的產(chǎn)品，在工序以及結(jié)構(gòu)相似的情況下，進行相同工序的合并，相似工藝結(jié)構(gòu)的產(chǎn)品合并，進行混流改造，可以滿足客戶需求，有效地提升生產(chǎn)效率，降低企業(yè)生產(chǎn)成本。同時，利用雙種群遺傳算法對工序的重新編排，可以擴大搜索空間，防止陷入局部最優(yōu)，對混流生產(chǎn)線平衡問題優(yōu)化有重大意義，該方法和技術可以推廣到其他類型的混流生產(chǎn)線設計與改造上。