机械优化设计案例11. 题目对一对单级圆柱齿轮减速器，以体积最小为目标进行优化设计。

2.已知条件已知数输入功p=58kw，输入转速n=1000r/min，齿数比u=5，1??]=400Mpa。

，许用弯曲应力[ 齿轮的许用应力[=550Mpa]FH3.建立优化模型3.1问题分析及设计变量的确定由已知条件得求在满足零件刚度和强度条件下，使减速器体积最小的各项设计参数。

由于齿轮和轴的尺寸（即壳体内的零件）是决定减速器体积的依据，故可按它们的体积之和最小的原则建立目标函数。

单机圆柱齿轮减速器的齿轮和轴的体积可近似的表示为：222222??)?0.25(b?c)(D?0.25?b(d?ddv?0.25(bd)?d?)22ggzz121222222????d?)?7c?0.258l(dd?dd22z101zzz22222222??)10mb(mzub?d?0?.25b[m?z0b?db?m.z8u1z2z1112222]3228dd?.6d)?d?l.205bd10?0.05b(mzu? m?12z1zz222z1z式中符号意义由结构图给出，其计算公式为d?mz,d?mz2211D?umz?10m1g2d?1.6d,d?0.25(umz?10m?1.6d)2z22g10z c?0.2b由上式知，齿数比给定之后，体积取决于b、z、m、l、d 和d z21 z1六个参数，则设计变量可取为TT]ddl[bzx[xxxxx]m?x?21z2314z5163.2目标函数为222222f(x)?0.785398(4.75xxx?85xxx?85xx?0.92xx?xx?5631323111122222220.8xxxx?1.6xxx?xx?xx?28x ?32x)?min61236541364563.3约束条件的建立z?z?17，得1 ）为避免发生根切，应有min10x??17?g(x)21b???????maxmin d的最大值，为齿宽系数齿宽应满足和2 )dmaxmin??，，得和最小值，一般取=1.4=0.9maxmin0(xx)?xg()?0.9?x322101.4??(x)?x(xx)g 3123 3）动力传递的齿轮模数应大于2mm，得02?x?g(x)?34）为了限制大齿轮的直径不至过大，小齿轮的直径不能大于4d，得axm10??300g(x)?xx325d?d?d）齿轮轴直径的范围：5得maxminzzz0gx?)?100?(x560?x??150g(x)570?130?xg(x)?680?)?x?200g(x69d??0.5?b?2ll按结构关系，应满足条件：）轴的支撑距离62minz? =20），得（可取min0??400?x?.5x?x(gx)41106）齿轮的接触应力和弯曲应力应不大于许用值，得7g(x)?1468250(xxx)?550?0123117098?400??0g(x)4222??12)0.85410x?..xxx(0169?06666?10?x221237 098?)(gx?400?04222?13)10177028240(xxx.?.?x10?.0?394x221322??][，得不大于许用值8）齿轮轴的最大挠度axm440??0.003x117x)?.04x(xxx)g(4324145??][不大于许用值，得9）齿轮轴的弯曲应力ww6x1012.85?12240510??5.x)?()?2.4?g(153xxx3526x102.85?112240??5.5)?()?6?10g(x163xxx362 4.优化方法的选择个约束条件的优化设计问题，16由于该问题有6个设计变量，Matlab采用传统的优化设计方法比较繁琐，比较复杂，所以选用函数来求解此非线性优化问题，避免了fmincon优化工具箱中的较为繁重的计算过程。

数学模型的求解5. 该优化设计的数学优化模型表示为：5.1.1将已知及数据代入上式，2222?85xxxx?85xx?0f(x)?.785398(4.75xxmin33211132222x?x.6xxxxxx?0.8xxx?1.092xx?5533461121661222 )xx?28?32?xx6645Subject to:g(x)?17?x?021g(x)?0.9?x(xx)?03212g(x)?x(xx)?1.4?03321g(x)?2?x?034g(x)?xx?300?0352g(x)?100? x?056g(x)?x?150?057g(x)?130?x?06830?xg?200(x)?690??40?x?0.5x?xg(x)410610550?x)?1468250(xx)?g(x1311270980?g(x)?400?4?222?12)?0.6666?10x.x?0854?10(xxx0.1692213270980xg()??400?42?2213).?0394?10177(xxx0.2824?0.?10x x22312440?003xx04x(xx)?0.)g(x?117.43451426x1021.85?12240?10??5.5(g(x)?)?2.4153xxx3256x10852.?1122405.?)?()6?10??5(gx163xxx362优化工具箱对数学模型进行程序求解5.1.2运用Matlabmyfun.m,Matlab首先在优化工具箱中编写目标函数的M文件返回x：处的函数值ffunction f = myfun(x)f=0.785398*(4.75*x(1)*x(2)^2*x(3)^2+85*x(1)*x(2)*x(3)^2-85*x(1)*x(3)^2+0.92* x(1)*x(6)^2-x(1)*x(5)^2+0.8*x(1)*x(2)*x(3)*x(6)-1.6*x(1)*x(3)*x(6)+x(4)*x(5)^2+ x(4)*x(6)^2+28*x(5)^2+32*x(6)^2)由于约束条件中有非线性约束，故需要编写一个描述非线性：文件约束条件的Mmycon.mfunction[c,ceq]=myobj(x)c=[17-x(2);0.9-x(1)/(x(2)*x(3));x(1)/(x(2)*x(3))-1.4;2-x(3);x(2)*x(3)-300;100-x(5);x( 5)-150;130-x(6);x(6)-200;x(1)+0.5*x(6)-x(4)-40;1486250/(x(2)*x(3)*sqrt(x(1)))-550; 7098/(x(1)*x(2)*x(3)^2*(0.169+0.006666*x(2)-0.0000854*x(2)^2))-400;7098/(x(1)* x(2)*x(3)^2*(0.2824+0.00177*x(2)-0.0000394*x(2)^2))-400;117.04*x(4)^4/(x(2)*x( 3)*x(5)^4)-0.003*x(4);(1/(x(5)^3))*sqrt((2850000*x(4)/(x(2)*x(3)))^2+2.4*10^12)-5 .5;(1/(x(6)^3))*sqrt((2850000*x(4)/(x(2)*x(3)))^2+6*10^13)-5.5];ceq=[];里输入：最后在command window给定初始值x0=[230;21;8;420;120;160];%4[x,fval,exitflag,output]=fmincon(@myfun,x0,[],[],[],[],[],[],@myobj,output)%调用优化过程5.1.3最优解以及结果分析运行结果如下图所示：由图可知，优化后的最终结果为x=[123.3565 99.8517 1.7561 147.3157 150.4904129.5096]7 f(x)=2.36e*10由于齿轮模数应为标准值，齿数必须为整数，其它参数也要5 进行圆整,所以最优解不能直接采用，按设计规范，经标准化和圆整后：x=[124 100 2 148 150 130]7 f(x)=6.16 *106.结果对比分析73，而优化后×10mm若按初始值减速器的体积V大约为6.3273,优化结果比初始值体积减少为：mm则为6.16×10 的体积V77)×100%10＝2.5%＝1－(6.16×10/6.32×Δν所以优化后的体积比未优化前减少了2.5%，说明优化结果相对比较成功。

7.学习心得体会学习机械优化设计课程的心得体会通过将近一学期的学习，对这门课有了初步的了解和认识，学期伊始，浏览全书，发现全是纯理论知识，觉得这门课会很枯燥，但是又回过头来想想，作为21世纪的大学生，要使自己适应社会需求，首先在做任何事之前都应该有正确的态度看待问题，把这些想法作为促使自己进步的动力，再去学习课本知识，效果应该很不一样，有了想法就付诸行动，随着对课本内容的学习跟老师的讲解，发现并不是像自己在学期初想的那样困难，特别是在老师介绍了一些与机械优化设计相关的计算机语言和计算机软件后，真正体会到科学优化设计的强大跟简洁明了，与传统优化设计方法相比较，大大提高了设计效率和质量。

传统设计方法常在调查分析的基础上，参照同类产品通过估算，经验类比或试验来确定初始设计方案，如不能满足指标要求，则进行反复分析计算—性能检验—参数修改，到满足设计指标要求为止。

整个传统设计过程就是人工凑试和定性分析比较的过程，是被动地重复分析产品性能，不是主动设计产品参数。

按照传统设计方法做出的设计方案，有改进余地，但不是最佳设计方案。

而现代化设计工作是借助电子计算机，，应用一些精确度较高的力学数值分析方法，优化软件进行分析计算，找最优设计方案，实现理论设计代替经验设计，用精确计算代替近似计算，用优化设计代替一般的安全寿命可行性设计。

在进行程序求解的过程中，因为是初学Matlab软件，对很多问题的关键点不能够掌握，非线性约束如何书写，上、下限如何选择，函数格式如何书写，变量未定义等等或大或小的问题，但是在一步步排除错误、重新编写程序的过程中，渐渐的对Mtalab熟悉起来，懂得了一些优化方法的简单计算过程和原理，省去了6繁琐复杂的优化计算过程在学完课程之后，反思自己在学习过程中的得失，深深体会到，不论在人生的哪个阶段，都要对自己负责，做任何事都要耐心，细致，“千里之行，始于足下”，学会在物欲横流的社会大潮中，坚持踏踏实实走好人生的每一步。

8.参考文献[1] 孙靖民,梁迎春. 机械优化设计. 北京:机械工业出版社,2006.[2] 濮良贵,纪名刚. 机械设计. 8版. 北京:高等教育出版社,2006.[3] 孙桓,陈作模,葛文杰. 机械原理. 7版. 北京:高等教育出版社,2006.[4]李涛,贺勇军,刘志俭. MATLAB工具箱应用指南—应用数学篇[M].北京:电子工业出版社,2000.7机械优化设计案例2复杂刀具优化设计数学模型的建立及算法改进摘要: 目的建立复杂刀具优化的数学模型,提高优化算法速度.方法采用优化设计与CAD相结合的方法. 结果与结论解决了传统刀具设计的缺点,改进后的算法速度大幅度提高.关键词: 数学模型；优化；算法在传统的刀具设计中,通过查表和经验公式来确定各种结构参数和几何参数,然后,反复计算来得到相对较优的刀具参数.这种方法使设计过程复杂费时,且得不到最优化的参数,设计出的刀具成本高,加工效率低.因而刀具的计算机辅助设计应采用优化设计与CAD相结合的方法,欲进行优化设计,必需首先建立刀具优化设计的数学模型,由于复杂刀具的种类繁多,结构变[1]，,因而需分门别类地建立模型此篇仅以轮切式拉刀化多样,优化目标不同.为例1 拉刀优化设计的数学模型在拉刀参数设计过程中需要选择的主要参数有拉削余量A,齿升量a,f齿距t,容屑槽形状和深度h,容屑系数k,同时工作齿数等,这些参数可分为两类,一类是独立参数,如拉削余量和容屑槽形状等,这些参数基本不受其他参数的影响.另一类参数是非独立参数,如齿升量、齿距、容屑槽深度、容屑系数等,这些参数既相互限制又相互依赖,第一类参数的选择比较容易.可以用经验公式和数据库来解决.第二类参数比较复杂,只有通过优化的方法才能得到较好的结果.粗切齿升量的选择是一个比较复杂的问题.增大a可使齿数减少,拉刀f长度变短,但同时又要求容屑槽深度增加.另外齿升量的增加又会引起拉削力的增大,受到拉床和拉刀拉应力的限制.齿距是决定拉刀长度的一个重要因素,t越大,拉刀越长,同时工作齿数越少.这样会在拉削过程中引起振动,生产效率低,降低刀具的使用寿命;t过小,又会使容屑空间变小,从而限制了齿升量的增大.其他参数如同时工作齿数z,容屑槽深度h,容屑系数k都是a和t的fi8函数,只有当a和t选择后才能确定.从上述参数分析可知,a和t是拉刀ff设计的关键,在a和t 之间应有一最佳组合值,使得a在拉床的额定应力和ff拉刀的许用应力范围内达到最高,即使拉刀的长度最小.1.1 目标函数的建立确定以af和t为优化的自变量,A为切削余量.拉刀长度是与拉削生产率、成本及其工艺性能有关的参数,拉刀越短对使用和制造越有利,因而取粗切齿[4]部分长度L作为优化目标F= minL(a,t) =tA/(2a). ff约束条件的建立(1)1.2h- 1.13(kaL )≥0. 2/1制槽1)容屑空间的限fw(2)式中 h是与t有关的参数;k为容屑系数,是与t和af有关的参数;Lw为拉削长度.F-pπD≥0. zcwzi/拉额定力的限制拉2)床e(3)式中 Fe为拉床额定拉力;Dw为拉削后孔直径;p为单位切削力;zi为同时工作齿数,zi=INT(Lw/t)+1;zc为组齿数.[σ] –2pD≥0. 用3)拉刀许拉应力制限的inwzi/zcdm(4)式中 [σ]为拉刀许用拉应力;dmin为拉刀最小直径.11 -z≥0. 时大4)最同工制限数齿作的i(5)9z - 3≥0. 制5)最小同时的限作齿数工i(6)25 -t≥0. 制的6)最大齿距限(7)t- 4≥0. 制限最小齿距的7)(8)8)弧形槽能保证稳定的分屑要求的最大齿升量h-a≥0;f(D,n,z) -a≥0. cffz(9)t- Int(2t)/2 = 0. 倍9)齿距整0.5的数应为(10)1.3 优化模型:a,t;自变量f:F=minL(a,t)=tA/2a; 目标函数ff:g(1)=h-1.13(kafLw)≥0;g(2)=F-pπD≥1/2约束方程wzi/zc e0;g(3)=[σ]-2pD≥0;g(4)=11-z≥0;g(5)=z-3≥wzi/zcdmin ii0;g(6)=25-t≥0;g(7)=t-4≥0;g(8)=h-a≥f0;g(9)=f(D,n,z)-a≥0;g(10)=t-Int(2t)/2=0. cfz2 优化算法102.1 标准算法复合形法是一种采用直接搜索方式求解非线性规划问题的数值计算方法,这个方法可以在N维非线性约束的空间中自动选择并改进设计点,该方法的:]3[2.一般步骤为1)m>n+1{x}(i=1, 2,…,n,n+1,…,m)构在可行域内生成个点i①初始顶点的形成,:可以人工选定,成初始复合形,这里需要注意两个问题②需要检验初始顶点是否满足约束条件,;即检验其可行性;也可随机产生2),f(x)≤将其由小到大的顺序重新编号计算各顶点的目标函数值,1f(x)≤…≤f(x); 2m3)xm-1,即心点点个的中坏点后复合形中其余确定除去最mi?11?;xc=xi1m?1m?4)xxx=x+α(x-x),α为映对中心点确定最坏点的映射点cccamm1.3;～,0.9一般取射系数5)x按则件,约某个束条:的可行性如果违背了射检验映点a(x+x)/2x,xa向中心移动一半距离,反复直至映射点把映射点a ac x;是可行点a6)f(x),x,完成一次计算新的可行点的函数值用它代替最坏点m a;回到第二步,迭代7)f(x)-f(x)<ε直到满足.,则终止以上称为复合形,重复以上过程1m法的“标准算法”,由于该算法的概念简单、容易实现,且能有效灵活地处理不等式约束问题,所以在结构化设计中得到广泛的应用.112.2 存在的问题把上述标准算法应用于工程实际时,就会发现它还存在以下几个问题:1)过多的可行性检验限制了其在优化设计中的有效应用.初始顶点生成和映射点的确定,都要进行可行性检验,在结构优化设计中,可行性检验其实质上就是结构分析过程,其计算量通常要占总工作量的80%以上,因此结构分析次数过多,必然会导致因计算时间过长而降低算法的效率.2)迭代过程中向极值点逼近的速度问题.开始若干次迭代(一次迭代是对于选取一个既满足约束条件又使目标函数值有所改善的新点所需的计算),目标函数值下降得很快,各顶点迅速接近极值点,一般来说,最初的(5～10)次迭代函数值下降得最快.随着迭代次数的增加,函数值的变化却越来越缓慢,也就是说,这时要使目标函数值有微小的改善,都要付出宝贵的计算时间.3)局部最优点问题.用上述算法得到的最优点有可能是局部最优点,虽然可通过多取几个初始点,经计算后得到几个最优点,然后比较得到全局最优点,但这样必然会导致计算工作量的成倍增加.2.3 分层复合形法针对标准算法中存在的问题,采用“分层复合形法”,它是对标准复合形法的改进,其基本思想是:充分利用复合形法开始时目标函数值急剧下降的特点,以迭代次数为控制参数,进行两层优化计算,为避免产生局部最优点,在第一层迭代中,选取多组复合形分别地进行计算,经过若干次有效地迭代,各顶点迅速地逼近最优点,分布在最优点附近.分层复合形法的基本步骤如下:(n )g, (i=1,}…,m)n,n1){x}…, {x(1)个初组初始顶点构成选择ggii n=Int[n/2]+1,n为设计变量数.始复合形,这里只要设计变量所取的值g不太小且相互间离得远些,就可不对初始顶点作可行性检验.2)2)~ 6)步的计算是第一层迭代对各初始复合形标准算法第,取映射率为αn;迭代次数为,t113)n次以后,第一层迭代结束当各复合形都迭代,取两个最好设计点组成t1α(α<α),迭代次数为,取映射率为新的复合形进入第二层迭代221n;.第二层迭代得到的最优点可被认为全局最优点t212分层复合形法有以下几个优点:①迭代次数大大减少;②以迭代次数为停止准则,可根据需要人工控制计算工作量;③第二层迭代能有效地产生全局最优点.3 结论依据本文所述方法,已开发出具有高效率优化CAD系统,证明对传统算法的改进是有效的.参考文献:[1] 唐锡荣. CAD/CAM技术[M].北京:北京航空航天大学出版社, 1994.18-36.[2] 蔡锁章.计算方法[M].北京:中国科学技术出版社, 1993. 54-60.[3] 徐灏.机械设计手册.第二卷[M].北京:机械工业出版社, 1991. 40-41.[4] 吴伏家,刘兆华.圆孔拉刀CAD系统研制[J].华北工学院院报, 1996,(增刊): 74-78.13机械优化设计案例3直齿圆柱齿轮传动的优化设计摘要：一、问题描述：现有一单级渐开线直齿圆柱齿轮减速器，其输入功率N=280kW，输入转速n=980r/min，传动比i=5。