Designand Research设计与研究 基于模块间约束的机床布局优化方法 李爱平何琪刘雪梅 (同济大学机械与能源工程学院,上海201804) 摘要:机床的占地面积是生产线布局的一个重要约束,为获得占地面积最小的机床布局方案,提出了基于 模块间约束的机床布局优化方法。首先将待优化机床三维模型主要模块映射为二维矩形布局图元, 并在此基础上构建了机床二维布局数学模型,然后根据不干涉算法设计了相邻算法,用来更精确地 表达机床模块间的不干涉和相邻位置约束,最后应用改进的遗传算法通过Matlab计算得到机床布 局优化后的方案。通过案例研究,得到较好的机床布局优化结果,验证了该优化方法的有效性。 关键词:机床布局 矩形图元中图分类号:TH122 相邻算法遗传算法 文献标识码:A A layout optimization method of machine tool with constraint between modules LI Aiping,HE Qi,LIU Xuemei (College of Mechanical and Energy Engineering,Ton ̄i University,Shanghai 20 1 804,CHN) Abstract:The land occupation of machine tool is one important constraint of production line layout.In order to ob— tain the least layout area of machine tool,a layout optimization method of machine tool with constraint between modules is proposed.Firstly,establish two-dimensional rectangular layout unit according to the main modules of machine tool and build the two—dimensional layout model of the machine too1.Then, design neighboring algorithm on the base of non——interference algorithm to express the non--interference constraint and neighboring constraint more accurately.At last,an improved genetic algorithms is used to optimize the layout of the machine tool by Maflab calculation.A detailed case is presented in this pa- per,and the result shows that the method is effective. Keywords:Machine Tool Layout;Rectangle Graph Unit;Neighboring Algorithm;Genetic Algorithm 机床布局是生产线布局设计中急需解决的问题。 性能最优与空间最小是机械系统布局的两个主要目 局中的图元作了精确的数学表示,并构建了数学模型, 给出了一些理论上的算法,并将改进的遗传算法应用 到卫星舱的布局问题中。徐义春 在研究了矩形布 局优化问题时设计了一种构造式算法——定位法,即 将一个矩形围绕另一个已经确定位置的矩形作为参照 标。在以满足性能要求为目的的机械系统布局方面, 国外Moon…和Kota 等人提出了分解运动变换矩阵 以实现机械结构自动重组。国内的张广鹏 以同样 的思路研究机床布局,根据刀具位姿矩阵导出了运动 进行部署,结合遗传算法寻优,得到了较好的计算结 果。Jacquenot 等在2009年的ASME设计工程国际 级联矩阵,进而生成运动功能方案,最终实现机床的结 构形态方案。韩泽光 研究了机械产品的运动方案 设计与结构方案设计,对这两个设计进行了精确的数 会议上提出了多目标自由布局问题的方法,对机械产 品的二维和三维布局都有一定的借鉴意义。 国内外研究学者在考虑性能或空间两者之一单个 因素的机床布局研究都已经较为成熟,但是对于如何 在机床空间优化布局的同时保证机床的性能这一问题 还有待研究。本文在对机床的布局优化中充分考虑了 性能与空间两方面的因素。在保证机床加工性能与操 学描述与表达,为机床的结构与布局方案设计提供了 一种理论设计方法。王德伦 提出运动构型、尺度与 结构方案相融合的机械产品创新设计状态空间方法体 系,并将其应用到机床的方案设计中。在以空间最小 为目的进行的机械系统布局方面,冯恩民等 对布 国家高档数控机床与基础制造装备科技重大专项(项目编号:2011ZX04015—022);上海市基础性研究重点项目(11JC1413200);中央高校基本 科研业务费专项资金 9・ 设计与研究 gnand Re,search 作性能的基础上,对机床的空间布局进行优化。具体 研究的机床为其与地面接触的各个模块(床身、电气 柜、排屑机构、刀库)的底面形状为矩形或可简化为矩 形的这一类机床,总体思路是将三维的机床划分成多 个二维矩形块作为矩形布局图元,以机床所占矩形包 络面积最小为优化目标,用数学方法精确描述机床二 维矩形块之间的约束。应用改进的遗传算法针对图元 进行布局优化计算,得到机床布局的优化结果。 根据以上的定义,可以通过 、0 、a 唯一确定矩 形图元和其在布局空间上的位置。即 R =R ( ,0 ,a。)={ +Al U +A 2 l A^∈l—n , 0 ], =1,2} n个二维矩形布局图元表示为R (i∈, )。m个二 维矩形功能图元不参与布局,它们的位置参数由二维 矩形布局图元所确定。表示为R (i∈lm),其中, = {1,2,…,m}。 1.3布局方案 l机床布局优化模型 1.1问题描述 对于布局参数a (i∈, )确定的布局图元,其布局 位置c 可由 、Oi确定,记c =( ,0 )∈R ,称 c=(c ,c ,…,C )∈R 为一个布局方案。 记P矶和Pm为矩形布局图元R 中的顶点坐标P ( ∈L)的横坐标和纵坐标;称IX(C)=(P 一P…) (Pyma-P…)为C的矩形包络,其中P…=max{Pm}, P =rain{P },P =max{P协},P =rain{Pm},i∈ , , ∈/4。 机床作为生产线上的主要加工设备,它的占地面 积影响着整个生产线的布局。机床所占的矩形包络面 积越小,生产线可以布局得越灵活、紧凑、高效。如何 使机床在满足功能要求的同时,合理布局各个功能模 块,使机床的占地面积达到最小便是本文研究的问题。 1.2机床模块数学模型构建 机床各个模块的相对位置的确定就是机床的布 局。在机床的布局过程中,首先要将机床划分为若干 个相对独立的模块。由于本文研究的机床与地面接触 的各个模块底面形状为矩形,可将三维机床的各个模 1.4约束描述 (1)不干涉约束 机床在布局的过程中,各个模块之间不能干涉,操 作工人占据的操作空问与机床的模块不能干涉。用数 学方法如下所述。 块简化为二维矩形布局图元(床身、立柱、电气柜等)。 操作工人在操作机床的时候会占据一定的操作空间, 机床各个模块在布局过程中为了不与这些操作空问相 干涉,本文将这些操作空间分隔出来,定义为二维矩形 功能图元(控制观察区,工件装夹区等)。经过以上步 骤,可得到n个二维矩形布局图元和m个二维矩形功 图元R 除去边界后的那部分称为 的内部,记为 int ,有int R = ( ,0 ,a )={ +AlH。+A 2l, IA ∈(一 。扪口 ),Ii2=1,2}cR ,机床模块间、机床功能区的不 干涉约束可表示为int尺 n int R,=0(i≠ ,i, E,…), In+ ={1,2,…,n+m}。 能图元,然后对矩形图元进行精确的数学表达。具体 如下所述。 (2)相邻约束 为保证机床的性能,要求机床特定模块之间具有 相邻的关系。如机床床身与立柱具有相邻关系,刀库 与主轴运动区中特定换刀的位置相邻,电气柜与控制 观察区相邻等。这种相邻关系的数学表达如下所述。 图元R 的边界记为bou R ,bou R。=尺 -int R 。任 以机床所占地的矩形区 域 为布局的对象,首先以 待布局区域的矩形 的型 心为坐标原点,建立直角坐 标系f 一。一Z (f,轴与矩形的 长平行,z:轴与矩形的宽平 行)。在矩形 上布n个二 意点JP与矩形尺 相邻约束表为P∈R ;任意两矩形相 邻表为bou R nbou R ≠0,且满足不干涉约束。根据 具体机床模块的约束,记点与矩形图元之间所有图元 约束为dot(PR ), ∈,…;记二维矩形图元之间所有相 邻约束为st(R。尺 ),其中i, ∈,…,i 。 1.5机床二维布局数学模型 维矩形布局图元R ,i∈, = 图1矩形图元Ri的示意图 {1,2,…, },记尺 的型心为 =( ,y )∈R ;R 的4 个顶点按逆时针方向分别记为P ∈R , ∈,4={1,2, 3,4},其中 P 与 轴的正方向所成角度记为0 ∈ (一丌/2,Ir/2);记平行P P 的单位矢量为H =(cosO , sinO ),根据右手定则确定与ll 正交的单位向量l, ,', = (一sinO ,eosO );记 的长宽分别为2口 20 2(0 2≤ 口 1);记a =(0 l,0 2),女口图1所示。 ・ 在满足不干涉约束和相邻约束的条件下,使矩形 包络值最小的机床布局问题的二维数学模型如下: 目标函数: (C)=(P…一P )(P 一P ) (1) 110 ・ ‘u10等 0卅】 Design0nd Research设计与研究 不干涉约束:int R nint置= (i#j,i, ∈,…) 其中,n+ ={1,2,…,n+m} 相邻约束:d0t( ),i∈In+ ,st(R ) 其中i j∈I….i≠j g( )=s(P .)s(P ,)。同理,可以计算得到矩形R 对 矩形 的相邻度记为g( 两个矩形尺 、 ,之间的 相邻度记为g =min{g( 叫),g( )}。 将以上算法进行运算,得到的结果与用Matlab画 图进行比较,验证了不干涉算法和相邻算法的正确性。 2.3改进的遗传算法 2约束表达算法与优化算法 2.1不干涉算法 根据王秀梅等 6]研究布局优化模型关于不干涉 本文针对机床带约束的布局优化问题,采用遗传 算法来搜寻全局最优解。遗传算法中对布局变量进行 算法的证明得到,机床布局问题的不干涉约束int n int R = (i≠ ,i, ∈,…)成立的充要条件为 k ∈,4 或z ∈,4,使得max{min{(P P ) /'it l z∈厶},min {Pik-eft ,n∥l k∈,4}}/>0,其中P (k∈厶)为尺 的4 个顶点,,l 为P Pm+I1的单位外法向量,规定P 川= PI。。 记 =R nR ,计算 的面积S( )来表示R 与 足两者之间的干涉程度,具体计算方法参见冯恩民 等l7 关于卫星舱布局研究的论文。 2.2相邻算法 (1)点与矩形的相邻度计算 记点P。与矩形尺 的相邻度为s(Po),则S(P。): min{lP。P f},P ∈bot R 。计算点与矩形相邻度的步 骤如下: 步骤一,初始化,分别输人点Po、R 的型心 、 , 旋转角 相邻度s(Po)=0; 步骤二,判断点 与矩形R 是否干涉,判断方法 与判断两个矩形是否干涉的方法相同(将点 看作 矩形,即Po=P扪k∈,4),若干涉,则s(P。)=0,并终 止,若不干涉,跳步骤三; 步骤三,计算 的各个顶点坐标P 和外法向量 ; 步骤四,计算点 到 各个顶点之间的向量 PoP 确定k ,使得I PoP I=min{IPoP I}; 步骤五,判断P0P ・,lm ≥0是否成立,若成 立,则计算P0到 , 、所在直线的距离s i (点到 直线的计算公式),s(P。)=s 。,终止,若不成立,转步 骤六: 步骤六,判断P。P ・,l叭l+1 ≥0是否成立,若成 立,则计算P。到P…川 ,P ,所在直线的距离s (点到 直线的计算公式),s(P。)=s i ,终止,若不成立,则s (P。)=IPoP I,终止。 (2)两个矩形之间的相邻度计算 矩形R 对矩形R 的相邻度记为g 的4个 顶点与矩形R 的相邻度分别为s(P ),k∈,4,在4个 相邻度中选取最小的两个相邻度s( ,),S(P ,),则 划、《 ≯ -?5 十 ‘ 十进制浮点数(实数)编码,采用单点杂交和正态变 异,终止准则为达到最大的进化代数或者迭代20代, 种群的质量都没有改进为止。在该遗传算法中,为了 使重新布局的机床满足性能要求,设计了罚函数,与原 目标函数结合作为整体的适应值函数,用来淘汰不满 足约束的布局方案。建立适应值函数,应用遗传算法 的主要步骤如下所述。 (1)适应值函数 在建立适应值函数过程中,考虑了机床占地面积、 图元之间的干涉量、点与矩形之间的相邻度和两个矩 形图元之间的相邻度四方面。适应值函数包括一个目 标项和3个约束项(罚函数),它们之间统一了量纲, 分配了相应的权重(罚因子),即 P 口 f(c )= (c )+∞ s( )+∑OJ2us (P )+∑C03vg “=o =o (2) 式中:5( )为干涉度,s(P, )为第u对点与矩形之间的 相邻度,g, 为第 对矩形之间的相邻度,P为需求相邻 度的点与矩形对的总数,q为需求相邻度的矩形对的 总数。09。、 、 为各项的权重,它们的值可采用层 次分析法根据实际机床的情况通过经验确定。 (2)染色体编码、单点杂交与正态变异 染色体编码:[ ,Y,0],其中 ,Y,0分别表示矩形 图元型心的 坐标、Y坐标和图元相对于z 轴转过的角 度: 单点杂交:对个体c.,c,,产生一个[0,1]的随机 r…若r <P ,则按均匀分布选取一个基因位,将c。、C 的前后两部分位串进行交换。 正态变异:对于个体C =c , =( , ,…, ), ~N(o, ),i∈f ,且 , ∈, 相互独立。 (3)改进的遗传算法的主要步骤 步骤一,初始化,设定最大的进化代数 ,种群规 模Ⅳ,矩形布局图元的个数n,矩形功能图元m,长宽 o ,需求相邻度的点与矩形对的总数P,需求相邻度的 矩形对的总数q,权重系数 。、 : 、09, ,令进化代数 为0,根据直角坐标,从DCR孤 中选取Ⅳ个个体点,作 得到 (P )。 为保证机床的加工性能,尺 (床身和工作台)中的 加工区与R (立柱、主轴和床身)中的主轴运动区相 邻。采用矩形之间的相邻度算法得到g 。 为便于对机床的操作观察,尺 (控制柜)与 (操 作观察区)相邻。采用矩形之间的相邻度算法得到 g5"7。 由机床以上约束知:P=1,q=2,结合式(2),建立 的适应值函数为 c )= (c )+ 1S( )+∞21 s(ed) + 31gl2+ 32g57 (3) 3.3确定权重 式(3)中60 、 :。、 、∞ 为各个约束的罚因子,也 可以看作是各个目标的权重。应用层次分析法,对该 机床的各个性能约束进行两两比较建立成对的比较矩 阵,如表2所示,经过归一化处理,令项 (c )前的系 数确定为1,可得到各个约束的权重分别为: =15. 137,∞2l:3.277, 3l=5.907,∞32=5.907。 表2布局约束的比较矩阵 性能 (C ) S(力) S(Pd) g12 g57 (C ) 1 1/9 l/5 l/7 1/7 S(力) 9 1 6 3 3 S(Pd) 5 1/6 1 l/2 1/2 g12 7 l/3 2 1 l g57 7 1/3 2 l 1 3.4优化计算 将机床的原布局方案Cn作为一个初始个体,应用 不干涉算法得到初始干涉量S( )=0,采用相邻度算 法得到初始s(P )=0,g =0,g ,:0,将干涉量和相邻 度值代人式(3)得初始适应值f(C )=2.64X10 mm 。 根据经验改变HMC63h机床各模块的布局位置,取可 行的9种布局方案c ,c ,…,Cq,将它们作为9个初 始个体,加上原布局方案 ,建立初始种群 ,c , …,cq,其规模Ⅳ=10。数学表达第后代种群,即 P =(Cl ,C2 ,…,Cl0 ) =( 1 l,Y1 1,0l l,…,CP1 ,Y1h,0l ,…, 1o l,Y10 1, 010 l,…, l0 ,Y1o ,0lo ) 其中 Ⅲ,YⅢ,Oi kj为一条染色体上的3个布局参 数,分别表示第k代群体中第i个个体点的第 个图元 的型心的 坐标、Y坐标和图元相对于Z.轴转过的角 度;当k=0时,P 表示初始种群,这时C =C ,i=1, 2,3,…,10。 输入得到的初始值N、n、m、a 、P、q、 1、 2l、0331、 ‘u 00等 平 I删 Desjgnand Research设计与研究 ,令 =20,应用Matlab按遗传算法步骤编程实现 杂交、变异和筛选的过程,计算收敛得到优化结果。由 于遗传算法变异与交叉的随机性,每次得到布局图元 的位置参数都会不同,为克服这一困难,将结果中R 的布局位置变换到原布局位置上,其他图元作相同的 变换,变换后可得到统一的布局结果。从对称性上对 布局结果略加调整,得到的机床布局结果如表3所示。 表3各图元布局优化后的位置、形状参数 Xil 0 口 Rl —l 300.5 —315.5 0 (3 632,1 556) 2 1 054 —315 O (1 077,2 080) 262.5 946.7 0 (2 647,277.559) 兄 -2 666.5 779.9  ̄r/2 (634,900) R5 -2 505 一l 673.1 (2 411,600) -3 566 —3l5.1 O (900,1 556) -572.5 一l 531.9 0 (1 456,1 029) 图6为最终的二维布局图,由计算得到最终的矩 形包络的值为1.63 X10 mm ,二维占地面积减少了 38.26%。对照图3和图6可以发现,位置变动的模块 有尺 (排屑机构)、R。(刀库)、R (电气柜),减少的大 部分占地面积是由于尺 (排屑机构)的位置变动。 图6最终的机床二维矩形布局图 图7优化后机床三维俯视图 根据得到的机床二维布局图,布置相应的机床三 设计与研究DesIgnand 维模块,得到图7机床三维俯视图和图8机床三维立 体图。 HMC63h机床经过二维矩形布局图元转换,模块 参考文献 [1]Moon Y M,Kota S.Automated synthesis of mechanism using dual—vee- tor algebra[J].Mechanism and Machine Theory,2002,37(2):146— 166. 问约束添加,优化计算3大步骤,得到优化的布局方 案。通过优化前后布局的对比,发现机床占地面积优 化很大程度来自机床辅助模块的位置调整。 [2]Moon YM,Kota S.Design of reconfigurable machine tools[J].Journal of Manufacturing Science and Engineering—Transactions of the ASME, 2002,124(2):480-483. [3]张广鹏,史文浩,黄玉美,等.数控车床结构布局形态方案创成研究 [J].中国机械工程,2003,14(21):1805—1808. 『4]韩泽光.机械系统运动与结构方案集成设计理论与方法[D].大连: 大连理工大学.2008. [5]王德伦,汪伟,刘彪,等.机构创新与应用[c].中国机构与机器科学 应用国际会议(2011 CCMMS)暨中国轻工机械协会科技研讨会论 文集,2011:10—12. [6]王秀梅,冯恩民,藤弘飞.旋转舱布局优化模型的主要性质及不干涉 算法[J].大连理工大学学报,1995,35(2):125—131. 『7]冯恩民,官召华,刘重阳,等.带性能约束的卫星舱布局问题改进遗 传算法[J].大连理工大学学报,2005,45(3):459—463. 图8优化后机床三维立体图 [8]徐义春,董方敏,刘勇,等.带平衡约束矩形布局优化问题的遗传算 法[J].模式识别与人工智能,2010,23(6):794—801. 4 结语 本文提出了相邻算法,引用了不干涉算法。给出了 改进的遗传算法。以机床所占的矩形包络面积最小为 [9]Jacquenot G,Bennis F,Maisonneuve JJ,et a1.2D—multi—objective placement algorithm for free—from components[EB/OL].[2010一O1— 29]. 第一作者:李爱平,女,1951年生,博士,教授,博 士生导师,研究方向为数字化设计与制造、制造系统集 优化目标,在不干涉约束和相邻约束条件下,给出了对 机床与地接触的底面模块可简化为矩形的这一类机床 布局优化的一般数学建模的方法。最后引用实例,对 HMC63h型号的机床进行了优化布局,得到了布局优 化后的机床,其占地面积减少38.26%,验证了该方法 对机床二维布局优化的有效性。 (上接第108页) 成及自动化等,已发表论文60余篇,出版专/编著3 鄙。 (编辑文章编号:130135 谭弘颖) (收稿日期.2012—07—30) 如果您想发表对本文的看法,请将文章编号填入读者意见调查表中的相应位鼍。 使工作台绕 轴倾斜10。,切割木板和塑料板做实验, 结果很理想,可以作为以后高速数控切割机床设计的 一缝的作用,便不会出现工件滑落摔坏的现象。 4.4实例验证 用于切割木材,当转速选n=9 000 r/min档位时, 个参考。为使切屑下滑路径短,可以如图1所示使 参1995,10(1):9—13. 工作台绕 坐标轴转动在l,一z平面内倾斜。 考文献 [1]赵俊伟,武良臣.超高速数控铣床的设计[J].湖北工学院学报, 工作台面振动为136~150 Hzl3],取f_-145 Hz,取工作 面振动的单振幅A=0.0l mm,动摩擦因数近似取木在 木上的静摩擦因数,则f0=0.3, 。=arctanO.3= 16.7O。,代入式(13)及式(16)可得台面及工作面倾斜 角范围为2.63。<or <32.21。。 [2]闻邦椿,李以农,张义民,等.振动利用工程[M].北京:科学出版社, 2005. [3]李志仁,戴大力,黄文波.木工机床振动质量的评价[J].林业机械与 木工设备,2001,29(6):35—37. 5 结语 w2sin/xo f Eh当台面倾斜角满足t.to-arcsin 4通讯作者:赵庆志,男,1962年生,博士,教授,硕 —<o/O< 士生导师,研究方向为数控技术、机电一体化技术,已 0 发表论文70篇。 (编辑余捷) (收稿日期:2011—08—10) 文章编号:130134 如果您想发表对本文的看法。请将文章编号填入读者意见调查表中的相应位置。 。。。 g 时,切屑顺利下滑且切割工件不会出现 摔坏现象。切屑薄层是指脆性材料的切屑,根据这个 倾斜角度范围,把图1所示机床一侧3个支脚垫起来, t 0 I^ 20 13 ̄1