第l5卷第1期 工业工程 Vo1.15 No.1 2012年2月 Industrial Engineering Journal February 2012 基于改进粒子群算法的制造单元设施布局问题研究 郑永前,丁奎学 (同济大学机械工程学院,上海201804) 摘要:为避免单元系统布局和单元内设施布局分开孤立研究所导致的问题解空间损失,利用并行工程的思想对单元 布局的两个环节集成考虑,对单元系统布局、单元内设施布置、设施摆放方向进行同时描述,并建立多目标集成优化 模型。针对模型的复杂性,设计了改进粒子群算法,算法吸收了遗传算法中的交叉操作算子,具有跳出局部最优解 的能力。最后通过求解单元设施布置实例,验证模型和算法的有效性。 关键词:制造单元;设施布局;粒子群算法;多目标优化 中图分类号:TH165 文献标志码:A 文章编号:1007-7375(2012)01.0125.06 Layout Design in Cellular Manufacturing Based on Improved Particle Swarm Optimization Zhang Yong—qian,Ding Kui-xue (School of Mechanical Engineering,Tongji University,Shanghai 201804,China) Abstract:Traditionally,the intra—cell and inter—cell layout problems in layout design of cellular manufac— turing systems are solved separately,which may result in a local optimal solution.In order to avoid the lo- cal optimal solution,these two problems are solved concurrently in this paper.This problem is formulated as an integer programming with multiple objectives.In the model,the orientation of the facilities,the in- tra-cell and inter—cell layout are described simultaneously.Due to the complexity of the model,an im- proved particle swarm optimization(PSO)algorithm is proposed.To improve its performance,the PSO al— gorithm is modiifed by adopting the crossover operator used in genetic algorithm.A simulation experiment veriifes the validity of the proposed method. Key words:cellular manufacturing;facility layout;particle swarln optimization(PSO);multi—objective optimization 单元生产方式是适应于目前多品种、少批量的 成,解决的是单元内设备布置问题;Wangl6 假定单 市场需求环境所产生的先进生产方式,它融合了 元内和单元间布局为u形,考虑解决单元内设施布 Job—shop高柔性与Flow-shop高生产率的特点,解决 局和单元系统布局的方法;Shahram等 7j假定各设 了传统生产方式生产率与生产柔性的矛盾 J。为构 施面积相等且各单元面积相等,对于设备和单元面 造单元生产系统,需研究以下几方面:1)单元构建; 积相差较大的实际问题,Shahram所提的方法已不再 2)单元内设备布局;3)单元系统布局。大量学者 适用。 对单元构建问题进行了研究:Yin等_2 考虑了多因素 本文利用并行工程的思想,研究了单元制造系 单元构建问题;Kaebernick等 考虑了单元数量的 统中的单元内设施布局和单元系统布局问题,提出 约束;王爱民等 考虑了加工设备的生产能力和任 了对单元设施布局的两个阶段进行集成优化的方 务交货期的约束。然而,人们对单元设施布局方面 法,建立多目标优化模型,并利用改进的粒子群算法 的研究较少。Attahiru等 假定单元系统布局已完 对模型进行求解。 收稿日期:201I一03.15 基金项目:上海市自然科学基金资助项目(10ZR1431700) 作者简介:郑永前(1965一),男,河南省人,副教授,博士,主要研究方向为生产系统工程与管理、先进制造技术、企业信息化技术 与应用. 工业工程 第15卷 1 制造单元设施布局问题 关信息和要求确定一个折中解(compromise),也称 妥协解。为了保证两个优化目标统一的量纲,本文 设施布局是指在一定的生产环境下,制造系统 采用归一化因子 加]。 设计人员根据生产目标确定系统中设施的布置形式 定义归一化因子叩, : 和位置 引。假定制造单元和设施均为矩形块状结 M M .p 构,其长和宽已知;各单元内的设施统一采用线性布 叼=1/∑∑∑ cI 1 J 1 P 1 eil(v/ti)ti, (5) 局,且设施有横向和纵向两种摆放方式,单元间采用 多行布局。 =1/∑Si, (6) I=l M M P 本文要解决的问题可描述为:生产单元构建阶段 已完成,现对单元内设施布局和单元系统布局进行优 arin z= ∑∑∑nL=1 =I P=I  ̄ceil( / )d +卢 s。 化,使得总物料搬运成本最小,车间面积利用率最高。 (7) 2 问题建模 式中,Ot, 分别为搬运距离和布局面积的权重, +/3=1。 为了描述方便,引入以下基本符号:P(P=1,2, 2.4约束分析 …,P)为零件编号;i,.『(i, =1,2,…, )为设施编 2。4.1单元内设施布局约束 号;( ,Y )为设施i的中心点坐标; 、 分别为设 单元内设施布局模型引入以下变量:h 、 分别 施 的长度和宽度;c(c=1,2,…,c)为单元编号。 表示设备i的横向和纵向长度,当设备横向摆放时, 2.1物料搬运距离计算 h :f , = i;否则,h = , =f ,d 为单元内相邻 离散型制造企业多采用托盘、容器等物料装载 设备间的安全距离,如图1所示。 工具,为了更符合实际应用情况,物料搬运距离可描 y 述为 M 村P ∑∑∑n ̄eeil( / )・d 。 (1) 式中,n 为加工零件P时,该零件从设施 到设 施 的转移次数; 为零件p的需求量; 为零件P 的搬运批量;ceil()表示对数值进行向上取整,d 为 x 设施 √之间的距离。工厂中物料一般是沿着通道 图1单元内布局参数示意图 流动的,本文取折线距离作为设施间距。 Fig.1 Parameters of intra-cell layout d =I 一 ,I+lY 一),,I。 (2) l — ,I≥(h + f)/2+d , (8) 2.2设施面积利用率计算 Y =乃。 (9) 设施布置应考虑紧凑性原则,即应提高设施的 其中,式(8)保证设备i和 保持一定的安全距 面积利用率,设施面积利用率可描述为 离,式(9)保证同一单元内的设备布置在同一行中。 max∑Si=1 /S。 (3) 2.4.2卑无 局约慕 根据物流顺畅的要求,各制造单元间要保持一 式中,Si表示设施i的矩形包络面积,S表示整 定间距,本文将单元间距d加入单元包络矩形,形成 体布局的矩形包络面积。由于各设施的矩形包络面 单元拓扑结构,如图2所示。 积为定值,所以该优化目标可等价于minS。 S=[max(x ̄)一rain( ∞)]・[max(), )一 min(Y )], i,J=1,2,…,m。 (4) 其中, 、 分别为设施i包络矩形的左下角的 、 坐标值, 分别为设施 包络矩形的右上角的 、Y坐标值。 2.3目标函数 本文要解决的单元设备布置问题有两个目标, x 属于多目标优化问题。求解多目标优化问题就是求 图2制造单元拓扑结构 出问题的Pareto最优解集,然后再由决策者根据相 Fig.2 The extented structure of manufacturing cell 第1期 郑永前,丁奎学:基于改进粒子群算法的制造单元设施布局问题研究 定义决策变量 ,。 f1,单元c在第r行 一问题解的一个集合开始。但与GA相比,PSO只有最 优粒子提供信息给其他粒子,信息的流动方向是单 向的,整个搜索过程为跟随当前最优解的过程,所有 粒子很有可能更快地收敛于最优解,但也相对容易 陷入局部最优解。文献[10]利用遗传算法对单元布 10,其他 (1o) (11) 生产单元应满足约束条件: r=l ∑S =1,c=1,2,…c; C 局问题求解,文献[1 1]求解单向环形设备布局问题。 本文在文献[10]、[11]的基础上,吸收遗传算法中 ∑S ≤C,r=1,2,…R; [ cb一( .b+ .b)]・[( 出+ cb)一 。,b]I>0, S ,=S ,,; (12) [Y b一(Y |h+ |h)]・[(Y b+ b)一Y .b]≥0, S ,≠S , ; (13) b+ ≤max ; (14) Y b+ ≤max Y。 (15) 式中,约束条件(10)表示各生产单元所在的行 是唯一的;式(11)表示每一行单元数小于等于单元 总数;式(12)表示若单元C和单元c 在同一行,则单 元的拓扑形状在 方向上不重叠;式(13)表示若单 元C和单元c 不在同一行,则单元的拓扑形状在Y方 向上不重叠,式(14)、(15)表示单元c必须满足布局 空间约束,不能超出所要求的布局空间。 3 问题求解与改进粒子群算法设计 3.1基本粒子群算法 粒子群算法兼有模拟进化算法和群智能的特 点,与其他进化算法相类似,PSO算法也是通过个体 间的协作与竞争,实现复杂空间中最优解的搜索。 在每一次迭代中,粒子通过跟踪两个极值来更新自 己:第一个极值是粒子本身找到的历史最优解,称为 个体极值点( );另一个极值是整个种群目前找到 的最好解,称为全局极值点(g )。在找到这两个最 好解后,粒子根据式(16)~(18)来更新自己的速度 和位置。 =n, +z1r1(Pb。。t— )+12r2(gb t一 ), (16) = + +1。 (17) 其中,Z ,Z 为学习参数,r。,r 为[0,1]区间上的 随机数, 为第 个粒子在第k次迭代中的位置,P 为第 个粒子自身的历史最优位置,g 为所有粒子 的历史最优位置,∞为设定的约束系数,定义60为 = 一( ax一 i )k/T, (18) 式中, ∞ i 为已知常量, 是迭代的总次数, k为已迭代次数。 3.2杂交粒子群算法 PSO与GA都具有隐含并行性,搜索过程都是从 的交叉操作,对粒子群算法进行改进。 粒子群中的粒子被赋予一个杂交概率P,该概 率由用户给出,不依粒子的适应度改变,在迭代中依 据杂交的概率选取指定数量的粒子,将选定的粒子 随机地两两杂交,产生相同数目的子代取代父代群 体,让a和b表示被选择的2个粒子,计算公式为 k = : +(1-A) , (19) = +(1-A) :, (20) k = : +(1一 ) , (21) k“= :“+(1-A) :。 (22) 以上各式中, 为0—1之间的随机数。 3.3粒子的结构化编码/解码 单元设施布局问题求解有2个难点:1)单元设 施布置问题,属于离散空间的非数值优化问题,而 PSO属于连续空间域的优化算法;2)一般的编码方 法只用于解决一个方面的问题,而对单元设施布局 问题要求设计的编码方式能同时描述单元系统布 局、单元内设施布置、设施摆放方向。本文设计了一 种粒子的多维实数编码方法,实现从粒子到设施布 局形式的映射,使PSO成功应用到离散的设备布局 问题求解中。 本文采用结构化编码方法,如表1所示,粒子编 码被分成3部分。 表1粒子位置编码示例 Tab.1 Example of the encoding of a particle 各部分编码/解码解释如下。 1)单元位置编码。该部分为长度为C的正实 数编码串,每一个编码位对应一个生产单元,解码时 对编码位上的实数大小进行比较,该实数的排列序 工业工程 第15卷 号即为对应的设备编号。本文采用自动换行策略, 若一行中各单元横向长度和单元间距离之和超出设 备布局范围,则该行的最后一个单元在下一行中进 行布置。 2)设备位置编码。该部分是长度为 的正实 数编码串,解码时对属于同一单元的设备所对应的 实数编码进行排序,实数的排列顺序即为该单元中 设备的编号的排列顺序。 3)设备摆放方向编码。该部分是长度为 的 正整数编码,若该整数为奇数,这单元横向放置;若 为偶数,则单元纵向放置。 该示例有7台设备( , =1,2,…,7),现已通过 单元构建得到单元划分结果为:C :{1,4,5},c = {3,6},C ={2,7}。单元位置编码为(1O.5,3.2, 11.6),解码后的单元排列次序为2—1—3;设备位置 编码为(2.1,1.3,13.6,0.3,8.2,3.5,1.6),解码后 的设备排列顺序为6—3 l5—1—4I7—2;设备摆放方 向编码为(2,7,3,4,15,19,10),解码后为纵一横一 横一纵一横一横一纵。 3.4问题求解和改进PSO算法流程设计 采用改进PSO求解单元设施布置问题,算法流 程设计如下。 步骤1 确定粒子种群和数量,初始化粒子位置 和速度,令g 。 未优化的次数n=0。 步骤2 判断算法是否达到最大迭代次数,若 是,则停止迭代,输出结果;否则,转步骤3。 步骤3 判断gbest是否已经连续Ⅳ代没有被优 化,若是,则停止迭代,输出结果;否则,转步骤4。 步骤4按式(16)一(18)对粒子的速度和位置 进行更新。 步骤5比较粒子的适应值和自身最优值P 。 如果当前值比Pbest更优,则置P 为当前值。 步骤6 比较粒子的适应值与种群最优值g 。 如果当前值比g 。 更优,则重置g 。 步骤7 选取各子群中的非最优粒子,按式 (19)~(22)进行交叉操作,更新各子群,转步骤2。 4实验与仿真分析 某车间接到一批生产任务,需要加工工件的种 类、批量、工艺路径如表2所示,各设备的面积如表3 所示。公司考虑到对现有的车间布局不能满足生产 要求,决定对车间设备重组,并进行重新布局(车间 面积为20 m x 18 m)。 表2零件生产工艺路径 Tab.2 The process paths ofthe parts 善需求量 工路艺径 表3设施面积 Tab.3 The area of the facilities m 编号 面积(长X宽) 编号 面积(长x宽) l 2.8×1.8 9 3.2×2.4 2 3.6×1.6 10 2.0×1.6 3 3.0×1.6 11 3.6×2.0 4 3.0×1.8 12 3.0×2.2 5 4.8×2.4 13 2.4×2.0 6 3.2×2.8 14 2.2×1.8 7 2.6×1.6 l5 4.2×2.2 8 2.8×1.8 现已知制造单元的划分结果:C ={m ,m ,m , m9,m13},c2={m3,m4,m7,m8,in15},c3={m6, m 。,m m。 ,m }。同一单元内,相邻设施间距离 为1.2 m,相邻单元的间距为3.0m。 文献[6]中假定单元构建已完成,以总物流费用 最小为目标,用模拟退火算法对单元内设施布局和 单元系统布局问题集成求解,不考虑车间面积利用 率问题。 本文算法设定粒子最大迭代次数gen=1 000, 粒子数量psosize=40,g 连续没有进化的代数N= 200,学习参数f =Z =2,最大惯性系数权重W = 0.9,最小惯性系数W i =0.4,设施位置编码数值范 围为[0,l0],单元位置编码数值范围为[0,1],设备 摆放方向编码数值范围为[0,100],目标函数中物料 搬运距离和布局面积所占权重分别为0.7、0.3。改 进粒子群算法将所有粒子随机地划分为4个种群, 杂交概率P=0.4, 取0—1之间的随机数。 分别用文献[6]中的方法、基本粒子群算法和改 进粒子群算法求解本文实例,3种算法各运行3O次 并取最优值,得到单元系统布局和单元内设施布局 结果,以及单元内各设施的摆放方向(1表示设施横 向摆放,0表示设施纵向摆放),结果如表4所示。 第1期 郑永前,丁奎学:基于改进粒子群算法的制造单元设施布局问题研究 129 可以看出,文献[6]的总搬运距离最短,但占地 Wang Ai-min,Ding Guo-zhi,Ning Ru—xin.Rapid desing 面积也是最大的,由本文模型得到的物料搬运距离 technology for manufactruing cells[J].Transactions of Bei— 略大,但占地面积减少很多;在给定的搬运距离和布 jing Institute of Technology,2006,26(10):850-854. 局面积的权重下,改进粒子群算法得到的结果是最 [5]Attahiru S A,Mingyuan C,Sunderesh S.Integrating the grouping and layout problems in cellulra manufacturing sys。 优的。 tems[J].Computers and Industirla Engineering,1992,23 5结束语 (4):55-58. [6]Wang Tai—yue,Lin Her-chang,Wu Kudi—bin.An improved 本文研究解决了制造单元的设施布局问题,对 simulated annealing for facility layout problems in cellular 单元系统布局和单元内设施布局集成考虑并进行多 manufacturing systems[J].Computers and Industiral Engi- 目标优化。单元布局问题属于NP完全问题,文中提 neering,1998,34(2):309-319. 出的基于结构化编码的改进PSO算法能进行有效的 [7]Shahram A,Napsiah I.An improved lagorithm for layout de。 求解。吸收了交叉算子的改进PSO算法既具有快速 sing in cellulra manufacturing systems[J].Journal of Manu— 搜索、鲁棒性好的特点,又具有使算法跳出局部最优 facturing Systems,2009,28(4):132—139. 解的能力。通过仿真实验可以看出,模型在兼顾设 [8]Russell D,Kai—yin G.The facility layout problem:Recent 施布局的面积利用率方面更具优势。 nad emerging trends and perspectives[J].Journal of Manu・ 本文所提出的单元布局方法适用于解决所生产 facturing Systems,1996,15(5):351—363. [9]锁小红,刘战强.制造系统设备布局的建模理论与求解方 产品的工艺路径已确定且产品的需求量比较稳定的 法[J].计算机集成制造系统,2007,13(10):1941—1951. 单元布局问题,对于多工艺生产路径和动态产品需 Suo Xiao—hong,Liu Zhan—qiang.Modeling and solution for 求的单元布局方法是今后研究的重点。 facility layout of manufacturing systems[J].Computer Inte— 参考文献: rgated Manufacturing Systems,2007,13(10):1941—1951. [10]王定益,王丽亚.一种单元制生产方式设备布局问题的 [1]Ah kioon S,Bulgak A A,Bektas T.Integrated cellular man— 研究方法[J].工业工程与管理,2006,16(1):10—13. ufacturing systems design with production plnaning and dy— Wang Ding-yi,Wang Li-ya.The facility layout study based namic system reconfiguration[J].European Journal of Oper- on cellulra manufacturing[J].Industirla Engineering and ational Research,2009,192(2):414-428. Management,2006,16(1):10—13. [2]Yin Yong,Yasuda K.Manufacturing ceNs’desing in con— [11]曾议,竺长安.基于群智能算法的设备布局离散优化研 sideration of various production factors[J].International 究[J].计算机集成制造系统,2007,13(3):541-552. Journal of Production Research,2002,40(4):885-906. Zeng Yi,Zhu Chang-an.Discrete optimization problem of [3]Kaebernick H,Bazargan—Lari M.An integrated approach to machine layout based on swarln intelligence lagorithm[J]. the design of cellulra manufacturing[J].CIRP Annals— Computer Integrated Manufacturing Systems,2007,13(3): Manufacturing Technology,1996,45(1):421-425. 541-552. [4]王爱民,丁国智,宁汝新.制造单元快速构建技术研究 [12]Kennedy J,Eberhart R.Particle Swarm Optimiztaion[C]. [J].北京理工大学学报,2006,26(10):850-854. Proceedings of IEEE International Conference on Neural 130 Networks,Australia:Perth,1995:1942—1948. 工业工程 第15卷 rithm for facility layout problems under variable demand in [13]靳雁霞,韩燮,周汉昌.改进的粒子群优化算法[J].计 算机工程与设计,2009,30(17):4074—4076. Jin Yan—xia,Han Xie,Zhou Han—chang.Improved particle cellular manufacturing systems[J].Computers in Indus— try,2001,46(2):181—188. [15]Ho Y C,Mooddie C L,A hybrid approach for concurrent swarl/1 optimization algorithm.Computer engineering and decision[J].2009,30(17):4074.4076. [14]Wang T Y,Wu K B,Liu Y W.A simulated annealing lago— (上接第124页) [12]张沙清,陈新度,陈庆新,等.不确定环境下模具制造项 目群随机调度[J].计算机集成制造统,2009,15(7): 1389.1396. Zhang Sha-qing,Chen Xin—du,Chen Qing—xin,et a1.Sto- chastic scheduling for multiple mould and die manufacturing projects under uncertainty[J].Computer Integrated Manu— facturing Systems,2009,15(7):1389-1396. [13]张沙清,陈新度,陈庆新,等.基于多步Q学习的模具制 造项目群随机调度算法[J].中国机械工程,2009,20 (12):1439-1445. Zhang Sha—qing,Chen Xin—du,Chen Qing-xin,et a1.Sto— chastic scheduling algorithm for multiple mould and die manufacturing projects based on multi—step Q learning[J]. China Mechanical Engineeirng,2009,20(12):1439—145. [14]李英杰,陈庆新,陈新度,等.多属性虚拟企业部分并行 协商项目规划[J].计算机集成制造系统,2005,11(6): layout design of cell and their flow paths in a tree configura— tion[J].International Journal of Production Research, 2000,38(4):895-928. 81O一817. Li Ying-jie,Chen Qing—xin,Chen Xin-du,et a1.Partilaly parallel multi-・issue negotiation ofr project scheduling in vir・ tual enterprise[J].Computer Integrated Manufacturing Sys— terns,2005,11(6):810-817,850. [15]Wu L H,Chen X,Chen X D.The research on proactive— reactive scheduling framework based on real・-time manufac・・ turing information[J].Materials Science Forum,2009(626— 627):789-794. [16]袁健.几种新的混合遗传算法研究[D].长沙:湖南大 学,2007. Yuan Jian.The Research of Several New Kinds of Hybrid Genetic Algorithm[D].Changsha:Hunan University, 2007. [17]王凌.车间作业调度及其遗传算法[M].北京:清华大学 出版社,2003,138—144.