第20卷第4期 哈尔滨理工大学学报 JOURNAL OF HARBIN UNIVERSITY OF SCIENCE AND TECHNOLOGY VoI.2O No.4 Aug.2015 2015年8月 遗传算法在机械优化设计中的应用研究 孙全颖, 王艺霖, 杜须韦 (哈尔滨理工大学机械动力工程学院,黑龙江哈尔滨150080) 摘要:为了寻求更优更高效的现代机械优化设计方法,依据遗传算法的基本理论,提出了一 种新型的智能优化方法.分析了遗传算法的基本原理、特点及运行步骤,并运用Matlab软件遗传算 法工具箱进行求解,结合蜗轮齿圈体积优化机械设计实例对遗传算法和传统优化方法中的序列二 次规划法进行对比,得出传统优化后体积是常规设计的75.28%,而采用遗传算法优化后体积仅为 常规设计的65.53%.由此可得出遗传算法的优化结果较传统优化算法的优化结果更优,过程更高 效的结论. 关键词:机械设计;遗传算法;Matlab优化工具箱;蜗轮蜗杆 DOI:10.15938/j.jhust.2015.04.009 中图分类号:TH122 文献标志码:A 文章编号:1007—2683(2015)04—0046-05 Application and Research of Genetic Algorithm in Mechanical Optimization SUN Quan—ying, WANG Yi—lin,DU Xu—wei (School of Mechanical Power Engineering,Harbin University of Science and Technology,Harbin 150080,China) Abstract:In order to search for a more better and more efficient modern mechanical optimization design meth— od,a new intelligent optimization method based on the basic theory of genetic algorithm is proposed.The basic principle,characteristics and operation steps of genetic algorithm are analyzed and using the genetic algorithm tool— box of MATLAB software tO solve,combined with the worm gear ring volume optimization mechanical design exam— pie of genetic algorithm and compared with the traditional optimization method of sequential quadratic programming method,it is concluded that the traditional optimized volume is 75.28%of the conventional design,and using ge— netic algorithm to optimize volume after only 65.53%of the conventional design.It is concluded that the optimiza— tion result of genetic algorithm is more obvious,the process is more efficient which is compared with traditional op— timization algorithm. Keywords:mechanical design;genetic algorithm;the Matlab optimization toolbox;WOITn gear and worm 合的设计方法,始于20世纪60年代,目的是提高工 0 引 言 优化设计是一种将优化原理与计算机技术相结 收稿日期:2014—08—27 程及其产品的设计质量和设计效率 .优化设计被 广泛的应用在机械设计中,在规定的工况及设计要 求条件下,满足机械产品的几何关系和性能约束的 作者简介:孙全颖(196O一),男,教授,硕士研究生导师,E—mail:mmjz1314@126.con 王艺霖(1989一),女,硕士研究生; 杜须韦(199O一),女,硕士研究生. 第4期 孙全颖等:遗传算法在机械优化设计中的应用研究 47 限制下,使机械设计的某项或某几项指标获得最 优值 . 机械优化设计的应用愈来愈广泛,同时也面临 着许多问题需要解决改进,其中寻求更优更高效的 设计方法便是其中之一l3 J. 近年来,随着计算机科学和智能理论的推广与 应用,出现了许多现代智能优化算法,到目前为止, 智能优化算法主要包括:模拟自然界遗传机制和生 物进化论的遗传算法,源于对鸟类捕食行为的研究 而发明的粒子群算法,基于模拟金属冷却过程的退 火法,模拟人类和动物记忆功能的禁忌搜索法,和以 人脑结构为参考模型的神经网络法等 J. 遗传算法由美国密执根大学的霍兰德教授提 出,现被广泛地应用与诸多领域,本文将在研究遗传 算法的运行机理基础上,结合具体的机械优化设计 问题,将遗传算法求解结果和特点进行分析,并与传 统优化算法求解进行比较,找出遗传算法的优 越性 |¨1. 1 遗传算法 遗传算法(genetic algorithm简称GA)是源于研 究生物系统过程中所进行的计算机模拟,由美国密 歇根大学的Holland教授在1975年首次提出.遗传 算法通过模拟达尔文生物进化论的自然选择和孟德 尔的遗传学机理的生物进化过程的计算模型,是一 种通过模拟自然进化过程搜索最优解的方法提出了 一种效率高、并行性强的全局性搜索方法,它能够自 动的搜索可行解空间的信息,通过全局搜索获得最 优解 . 遗传算法与传统算法的不同点在于遗传算法作 用对象不是参数本身,而是对变量运用编码技术得 到的基因个体,遗传算法可以直接对各种结构的对 象进行求解,遗传算法的适应度函数摆脱了函数连 续性和可微性的束缚,只需凭借适应度函数即可来 评价最优单位个体,并行性强,鲁棒性好,避免使计 算过早收敛于局部最优解 . 遗传算法主要有3种算子:选择、交叉、变异,其 数学模型如下: SGA=(C, ,Pn, , ,厂, , ). 其中:C表示个体的编码方法;E表示个体适应度评 价函数;尸n表示初始种群; 表示种群大小; 表示 选择算子;F表示交叉算子; 表示变异算子;T表 示遗传运算终止条件_1 驯.遗传算法基本流程图如 图1所示 图1遗传算法基本流程图 可行解的编码、遗传算子的选取、适应度函数的 确定是运用遗传算法求解优化问题的关键因素,因 为编码方式与遗传算子的不同配比形式对于整体问 题有不同的求解方法,故对于不同的优化问题需要 应用不同的求解方法. 遗传算法的步骤: 1)染色体编码与解码.编码对应关系为: 000000…0000=O— l 000000…0001:1_+ 1十 000000…0010=2一 1+26 111111…1111=2 一1一 6= . 解码公式: :U +( 一 ) . 2)个体适应度的检查评估.将目标函数转化为 适应度,通过适应度对函数进行评价,寻求最优值. 3)遗传算子.选择算子:设种群数为 ,个体i 的适应度 ,则个体i被选取的概率 P = /∑ , 1,2,…,n. 48 哈尔滨理工大学学报 第20卷 交叉算子:依交换概率P ,随机选择个体进行 交换操作.交叉体现了信息交换的思想.交叉操作如 图2所示. 单点 菜个 体[l Io1l J —00 —}r=兰—兰 0二缸1l l 11 1个体 图2单点交叉示意图 变异算子:根据突变率,对种群中个体进行变异 操作.如图3所示. ——11011]L 11.,t 图3变异操作示意图 2 遗传算法与传统算法求解分析 2.1机械优化设计实例 以某电梯曳引机传动装置的普通圆柱蜗杆减 速器为例,其输入功率P。:5.5 kW,转速 。= 1 400 r/min,传动比i=20,载荷系数K=1.1,蜗杆 选用低碳合金钢20CrMnTi,蜗轮选用铸造锡青铜 ZCuSnl0Pbl,蜗轮齿圈许用应力[or ]=220 MPa.优 化目的:在满足各方强度、刚度、稳定性的条件下,将 蜗轮齿圈体积最小作为设计目标,以减少贵重的青 铜材料成本. 2.1.1建立目标函数和选取设计变量 蜗轮齿圈体积为: (1) : . b≤ 6d。l, d。l=(q+2)m, d =d +tf,。m, (2) d0= 一2m, d。=(z2+2)m=(iz1+2)m, d,=(z2—2.4)m=(izl一2.4)m 其中:b是齿宽,d 是蜗杆齿顶圆直径,d 是蜗轮齿 顶圆直径, 是蜗轮齿根圆直径,。。是蜗杆头数,z 是蜗轮齿数.蜗轮齿圈结构尺寸如图4所示. 图4蜗轮齿圈结构尺寸 体积公式可转换成: : [(i 。+ +2):一(i 一4.4)z], (3) 因此,选取蜗杆头数z。,模数m,直径系数g为设计 变量,记作 X=( l, 2, 3) =( l,m,q) 根据式(3),目标函数改写成: ): [( 。 +2)z一( 。一44)z]. (4) 2.1.2确定约束条件 蜗轮齿面接触强度条件限制,见式(5),其中载 荷系数K=1~1.4,当精度等级大于7级,齿面滑动 的速度不高于3 m/s,并且载荷平稳时,K取最小值, [盯 ]是蜗轮齿圈材料的许用接触应力,其中蜗杆的 传动效率和蜗轮的扭矩见式(6)和式(7). m ≥ 3丽, (5) 田~-O.0354 ̄,(6) P, =i叼 1=9550ir/ , (7) A: %/ ̄丽tl+厂F2,l≤l。一 d,, (8) 2 2 “一d.一mq’ (9) 2 tan20。 2 tan20。 一 d, 一 。m ’ (10) E=2.1×10 MPa. (11) .J, = 簖= m (g一2.4)m Lg一・4)’ . (12) 第4期 孙全颖等:遗传算法在机械优化设计中的应用研究 49 模数约束:对于功率处于中小型的蜗杆传动, 的优化结果圆整,蜗杆头数z =3,模数m=4 mm,直 3≤m≤5; 蜗杆直径系数约束:对于功率处于中小型的蜗 杆传动,5≤g≤16; 蜗杆头数约束:对于功率处于中小型的蜗杆传 动,2≤z1≤4; 综上,优化实例约束条件如式(13)所示: ( 2-X 加√(x ̄x2]一2T2ta—n2—0 ̄)2・ 0.729i。 ;一157.51r 2( 2—2.4) ≤0, g3( )= l一5≤0, g4( )=3一 1≤0, g5( )= 2—16≤O, g6( ):5一 2≤O, g7( )= 3—4≤0, g8( )=2一 3≤0. (13) 2.2遗传算法求解分析 Matlab遗传算法工具箱集合了一些列函数,其 中ga.Ill是工具箱的核心函数,提供遗传算法工具 箱与外部的接口,在Matlab软件运行环境下,设定 参数,执行M文件,完成优化.工具箱中遗传算法的 主函数为:[x fva1] ga(@fitness,nvars,options)其 中X表示最终值达到的点;fval表示适应度函数在 最终点X的值.@myfitness表示适应度函数M文件 的句柄,nvars表示适应度函数变量个数,options是 遗传算法参数结构体. 遗传算法工具有一个图形用户GUI如图所示, 可以在命令窗口键入“gatool”直接调用,免去了在命 令窗口调用函数等复杂过程. 编写适应度函数M文件,运用该工具箱进行求 解,采用二进制编码,初始种群选取20,最大迭代次 数取100代,交叉概率0.8,变异概率取0.2,计算结 果如下:迭代次数止于66代时,目标函数最优值 X=( 1, 2, 3) =(zl,m,q) =(2.973,4.214, 13.144) 蜗轮齿圈体积 in=622 455.108 0 mm , 求解结果和最佳适应度与最佳个体分布如图5、图6 所示;迭代次数止于100代时,目标函数最优值X= (2.976,4.047,14.808) 。蜗轮齿圈体积 ; = 612 515.823 3 mm ,求解结果和最佳适应度与最佳 个体分布如图7、图8所示,对迭代次数止于100代 径系数g=15,V:602 696.498 4 mm . 图5迭代66次的求解结果 1O 越 s 略 20.40 60 80 l0o 迭代次数 培l5 r L,S 10 刍s 图7迭代100次的求解结果 ×10 0 2 嘣0.5 O 0 20 40 60 80 100 迭代次数 15 茎Ⅷ s 0 . 1 2 3 Number of variables 图8迭代100次的最佳适应度与最佳个体 50 哈尔滨理工大学学报 第20卷 2.3传统优化算法求解分析 结合Matlab软件,上述优化数学模型是属于约 束非线性规划问题,采用序列二次规划法(SQP),运 用求解工具箱中的fmincon函数,在极值点存在的 必要条件下来求极值点,在已知一系列二次规划问 题的解的前提下,逼近Kuhn—Tucker条件方程的 解,其调用格式为:[x,fva1]=fmincon(fun,x0,A,b, Aeq,beq,lb,ub,nonlcon,options,P1,P2…). 设初始值X =[2,5,18r,变量的下限lb=[2, 3,5],变量的上限ub=[3,5,16];采用传统优化算 法,调用fmincon函数得到最优解X=( 1, 2, 3) : (z ,m,q) =(3.000 0,5.000 0,7.727 7) ,V= 673 921.337 4 mm ,对优化结果按设计要求圆整,蜗 杆头数z.=3,模数m=5 mm,直径系数q=8,得出 蜗轮齿圈体积V=692 436.234 4 mm . 2.4求解结果分析 按照常规方法设计结果为蜗杆头数 =2,模数 m=5 mm,直径系数q=18,得出蜗轮齿圈体积 = 919 759.968 8 mm ,传统优化后体积是常规设计的 75.28%,而采用遗传算法优化后体积仅为常规设计 的65.53%,对比得出采用遗传算法优化蜗轮齿圈 体积的更为理想,优化结果较明显. 表1求解结果对比表 3 结 论 本文通过对机械优化设计领域现状的介绍和对 遗传算法基本理论的分析,结合具体曳引机中蜗轮 蜗杆优化实例,分别运用遗传算法和SQP序列二次 规划法进行优化,首先在求解方法上,通过利用 Matlab软件的遗传算法工具箱对比SQP法调用 fmincon函数,前者只需编写目标函数M文件,省去 了大量优化算法程序的编写,大大提高了设计效率; 其次在求解过程方面,遗传算法求解过程更简单化、 直接化,不受目标函数可微性的限制,全局搜索能力 高,避免了局部最优解等问题;最后将两种方法求解 结果与常规设计结果对比分析,得出遗传算法优化 结果更为明显.因此,本文提出运用遗传算来优化蜗 轮蜗杆传动装置,该方法为机械优化设计提供了一 种新的尝试,然而对于不同类型的优化问题,遗传算 法不可能优于其他所有算法,因此不断总结分析各 种优化问题的共同特征和求解遇到的共性问题,从 而来选择最适用的优化方法将成为未来机械优化研 究的主要方向之一,新型优化方法的研究具有一定 的经济价值和现实意义. 参考文献: [1]孙全颖.机械优化设计[M].哈尔滨:哈尔滨工业大学出版社, 2007:1—5. [2]裴飞飞.量子遗传算法在机械优化问题中的应用研究[D].武 汉:武汉科技大学,2010:1—55 [3] 贺春华.传统优化算法与竞选算法应用于机械优化设计的比 较[j].广东工业大学学报,2010,27(3):46—49. [4] 海丽切木阿不来提.浅谈几种智能优化算法[J].电脑知识与 技术,2011,7(19):4628—4630. [5] VOSE M.D.Modeling Simple Genetic Algorithms,Foundations of Genetic Algoritms[c]//Morgan Kaufmann Publishers,1993:63 73 [6] 葛继科.遗传算法研究综述[J].计算机应用研究,2008,25 (10):2911—2915. [7]IJl J,FENG Z P,CHANG J Z,et a1.Aerodynamic Optimum Design of Transonic Turbine Cascades Using Genetic Algorithms[J].Jour— hal ofThermal Science,1997,6(2):】11一l16. [8] 吕德文.MATLAB遗传算法工具箱的应用与研究[J].湖南农 机,2013,40(3):130—131. [9]李书全.遗传算法中的交叉算子的述评[J].计算机工程与应 用,2012,48(1):36—39. [10] 巩敦卫,孙晓燕.基于模式定理的遗传算法交叉和变异概率 上限[J].控制与决策,2004,19(5):0554—0556. [11] 曹俊,朱如鹏.遗传算法的发展及其存机械程领域中的应用 [J].机械技术史,2002,42(5):410—416. [12]雷英杰.MATLAB遗传算法工具箱及应用[M].西安,西安电 子科技大学出版社,2004:149—155. [13] 褚洪生,杜增吉,阎金华,等.MATLAB7.2优化设计例指导教 程[M].北京,机械工业出版社,2007:198. [14] 宋茂福,赵勇.基于MATLAB遗传算法工具箱的圆柱螺旋弹 簧迷糊可靠性优化设计[J].机械,2008,35(8):l一4. [15] 何大阔,王福利,贾明兴.改进的遗传算法在优化设计中的应 用[J].东北大学学报:自然科学版,2005,26(12):1123—1126. [16] 黄康,许志伟.改进型的遗传算法在机械优化设计中的应用 研究[J].现代制造工程,2005(10):80—81. [17] 余成树,吴作伟,高永芳,等.改进遗传算法在结构优化中的 应用[J].机械工程与自动化,2007,141(2):86—87. [18]林远艳,徐成福.基于改进遗传算法的行星轮系的优化设计 [J].桂林航天工业高等专科学校学报,2007,45(1):51—53. [19] 郭卫,赵栓峰,杨桂红.基于遗传算法的叠弹簧性联轴器的模 糊优化设计[J].机械传动,2004,28(4):13一l6. [2O] 宋光辉,阮米庆.基于混合遗传算法的两轴式机械变速器的 优化设计[J].机械工程师,2008(2):46—49 (编辑:关毅)