机械优化设计案例11. 题目对一对单级圆柱齿轮减速器,以体积最小为目标进行优化设计。
2.已知条件已知数输入功p=58kw,输入转速n=1000r/min,齿数比1?]=550Mpa,许用弯用应力[曲应力u=5,齿轮的许H?]=400Mpa。
[ F3.建立优化模型3.1问题分析及设计变量的确定由已知条件得求在满足零件刚度和强度条件下,使减速器体积最小的各项设计参数。
由于齿轮和轴的尺寸(即壳体内的零件)是决定减速器体积的依据,故可按它们的体积之和最小的原则建立目标函数。
单机圆柱齿轮减速器的齿轮和轴的体积可近似的表示为:222222??)?0.25(b?c)(.25Db(d?d?dv?0.25)b(d??d)?02gzz1g122222222????d?)?0.257l(d8?dddc?2112 zzzz022222222??)10m(mzu?d?b.25?[m0zb?d.b?m8zbub0?1112zz12222]3228dd6d)?d?l?05bd.?005 b(mzu?10m?1..2 2zz2zz2z121式中符号意义由结构图给出,其计算公式为d?mz,d?mz2112D?umz?10m12g d?1.6d,d?0.25(umz?10m?1.6d)2z2g210z c?0.2b由上式知,齿数比给定之后,体积取决于b、z、m、l、d 和z11d 六个参数,则设计变量可取为z2TT]ddbzmxxxxx]l?[xx?[23145z61z213.2目标函数为222222f(x)?0.785398(4.75xxx?85xxx?85xx?0.92xx?xx?5231116233112222220.8xxxx?1.6xxx?xx?xx?28x ?32x)?min6646213316545约束条件的建立3.3.zz?17?,得1)为避免发生根切,应有min0??17?xg(x)21b???????maxmin d的最大值为齿宽系数2 )齿宽应满足和,dmaxmin??,,得和最小值,一般取=1.4=0.9maxmin g(x)?0.9?x(xx)?03212g(x)?x(xx)?1.4?031323)动力传递的齿轮模数应大于2mm,得g(x)?2?x?0344)为了限制大齿轮的直径不至过大,小齿轮的直径不能大于d,得max1g(x)?xx?300?0352d?d?d5)齿轮轴直径的范围:得maxzminzz0?100?xxg()?560?x150?g(x)?570?x?g(x)?130680200?x)?x?g(69l按结构关系,应距离满足条件:撑6)轴的支?b?2??0.5d?l=20),得(可取2zminmin g(x)?x?0.5x?x?40?041610)齿轮的接触应力和弯曲应力应不大于许用值,得7.0550?xxx)?1468250g(x)?(1231170980?x)??400g(42??2212)x?0.854?10xxxx(0.169?0.6666?102223170 980g(x??400?)4?22213)x?10?0.?xxx(0.2824?0.17710394x23221??][ 8)齿轮轴的最大挠度,得不大于许用值max440?.003xxx(xx)?0g(x)?117.04 4521443??][ 9)齿轮轴的弯曲应力,得不大于许用值ww6x?102.8512124?5.5?2.4?100g(x)?()?153xxx3526x1085?12.2124?5.5?10?0?g(x)()?6163xxx3624.优化方法的选择由于该问题有6个设计变量,16个约束条件的优化设计问题,采用传统的优化设计方法比较繁琐,比较复杂,所以选用Matlab优化工具箱中的fmincon函数来求解此非线性优化问题,避免了较为繁重的计算过程。
5.数学模型的求解5.1.1将已知及数据代入上式,该优化设计的数学优化模型表示为:2222minf(x)?0.785398(4.75xxx?85xxx?85xx?311132322220.92xx?xx?0.8xxxx?1.6xxx?xx5631616154312222 ?xx?28x?32x)6465Subject to:0?17?x(gx)?210xx)?0.9?x(g(x)?32210.4?(xx)?1g(x)?x33120?x?g(x)?2340?300?x)?xxg(3520?x ?(x)?100g560150?x)?x?g(570x?x)?130?g(680??200(x)g?x690??405x?xg(x)?x?0.410160550?xx) ??g(x)1468250(x131127098?400)??0g(x4??22212)10x.854(0.169?0.6666?10?x?0xxx231227098?400? 0xg()?4?22213)x394?.2824?0177?1010x?0.xxx(0.2122344g(x)?117.04x(xxx)?0.003x?0 44514326x?1012.852124?5.5?0)?2.4?10(gx)?(153xxx3526x10.85?122124?5.10?5?0)g(x?()?6163xxx326 5.1.2运用Matlab优化工具箱对数学模型进行程序求解首先在Matlab优化工具箱中编写目标函数的M文件myfun.m,返回x处的函数值f:function 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)由于约束条件中有非线性约束,故需要编写一个描述非线性约束条件的M文件mycon.m:function[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(6)-200。
130-x(6)。
x(5)-150。
100-x(5)。
x(2)*x(3)-300。
x(3).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]。
%给定初始值[x,fval,exitflag,output]=fmincon(@myfun,x0,[],[],[],[],[],[],@myobj,output)%调用优化过程5.1.3最优解以及结果分析运行结果如下图所示:由图可知,优化后的最终结果为x=[123.3565 99.8517 1.7561 147.3157 150.4904 129.5096]7 f(x)=2.36e*10由于齿轮模数应为标准值,齿数必须为整数,其它参数也要进行圆整,所以最优解不能直接采用,按设计规范,经标准化和圆整后:x=[124 100 2 148 150 130]7f(x)=6.16 *106.结果对比分析73,而优化后mm大约为6.32×10V若按初始值减速器的体积73,优化结果比初始值体积减少为:10mm 的体积V则为6.16×77)×100%10/6.32×=101Δν=-(6.16×2.5%所以优化后的体积比未优化前减少了2.5%,说明优化结果相对比较成功。
7.学习心得体会学习机械优化设计课程的心得体会通过将近一学期的学习,对这门课有了初步的了解和认识,学期伊始,浏览全书,发现全是纯理论知识,觉得这门课会很枯燥,但是又回过头来想想,作为21世纪的大学生,要使自己适应社会需求,首先在做任何事之前都应该有正确的态度看待问题,把这些想法作为促使自己进步的动力,再去学习课本知识,效果应该很不一样,有了想法就付诸行动,随着对课本内容的学习跟老师的讲解,发现并不是像自己在学期初想的那样困难,特别是在老师介绍了一些与机械优化设计相关的计算机语言和计算机软件后,真正体会到科学优化设计的强大跟简洁明了,与传统优化设计方法相比较,大大提高了设计效率和质量。
传统设计方法常在调查分析的基础上,参照同类产品通过估算,经验类比或实验来确定初始设计方案,如不能满足指标要.求,则进行反复分析计算—性能检验—参数修改,到满足设计指标要求为止。
整个传统设计过程就是人工凑试和定性分析比较的过程,是被动地重复分析产品性能,不是主动设计产品参数。
按照传统设计方法做出的设计方案,有改进余地,但不是最佳设计方案。
而现代化设计工作是借助电子计算机,,应用一些精确度较高的力学数值分析方法,优化软件进行分析计算,找最优设计方案,实现理论设计代替经验设计,用精确计算代替近似计算,用优化设计代替一般的安全寿命可行性设计。
在进行程序求解的过程中,因为是初学Matlab软件,对很多问题的关键点不能够掌握,非线性约束如何书写,上、下限如何选择,函数格式如何书写,变量未定义等等或大或小的问题,但是在一步步排除错误、重新编写程序的过程中,渐渐的对Mtalab熟悉起来,懂得了一些优化方法的简单计算过程和原理,省去了繁琐复杂的优化计算过程在学完课程之后,反思自己在学习过程中的得失,深深体会到,不论在人生的哪个阶段,都要对自己负责,做任何事都要耐心,细致,“千里之行,始于足下”,学会在物欲横流的社会大潮中,坚持踏踏实实走好人生的每一步。
8.参考文献[1] 孙靖民,梁迎春. 机械优化设计. 北京:机械工业出版社,2006.[2] 濮良贵,纪名刚. 机械设计. 8版. 北京:高等教育出版社,2006.[3] 孙桓,陈作模,葛文杰. 机械原理. 7版. 北京:高等教育出版社,2006.[4]李涛,贺勇军,刘志俭. MATLAB工具箱应用指南—应用数学篇[M].北京:电子工业出版社,2000.机械优化设计案例2复杂刀具优化设计数学模型的建立及算法改进摘要: 目的建立复杂刀具优化的数学模型,提高优化算法速度.方法采用优化设计与CAD相结合的方法.结果与结论解决了传统刀具设计的缺点,改进后的算法速度大幅度提高.关键词: 数学模型;优化;算法在传统的刀具设计中,通过查表和经验公式来确定各种结构参数和几何参数,然后,反复计算来得到相对较优的刀具参数.这种方法使设计过程复杂费时,且得不到最优化的参数,设计出的刀具成本高,加工效率低.因而刀具的计算机辅助设计应采用优化设计与CAD相结合的方法,欲进行优化设计,必需首先建立刀具优化设计的数学模型,由于复杂刀具的种类繁多,结构变[1],此篇仅以轮切式拉因而需分门别类地建立模型化多样,优化目标不同,.刀为例1 拉刀优化设计的数学模型在拉刀参数设计过程中需要选择的主要参数有拉削余量A,齿升量a,f齿距t,容屑槽形状和深度h,容屑系数k,同时工作齿数等,这些参数可分为两类,一类是独立参数,如拉削余量和容屑槽形状等,这些参数基本不受其他如齿升量、齿距、容屑槽深度、容,另一类参数是非独立参数.参数的影响屑系数等,这些参数既相互限制又相互依赖,第一类参数的选择比较容易.可以用经验公式和数据库来解决.第二类参数比较复杂,只有通过优化的方法才能得到较好的结果.粗切齿升量的选择是一个比较复杂的问题.增大a可使齿数减少,拉刀f长度变短,但同时又要求容屑槽深度增加.另外齿升量的增加又会引起拉削力的增大,受到拉床和拉刀拉应力的限制.齿距是决定拉刀长度的一个重要因素,t越大,拉刀越长,同时工作齿数越少.这样会在拉削过程中引起振动,生产效率低,降低刀具的使用寿命。