机床主轴结构优化设计指导老师： 姓名： 学号：机床主轴结构优化设计一、 机械优化设计的一般过程建立优化设计的数学模型选择适当的优化方法编写计算机程序准备必要的初始数据并上机计算对计算机求得的结果进行必要的分析其中，建立优化设计的数学模型是首要的和关键的一步，其基本原则有：1、 设计变量的选择在充分了解设计要求的基础上，应根据各设计参数对目标函数的影响程度认真分析其主次，尽量减少设计变量的数目，以简化优化设计问题。另外，还应注意设计变量应当相互独立，否则会使目标函数出现“山脊”或“沟谷”，给优化带来困难。2、 目标函数的确定常取其中最主要的指标作为目标函数，而其余的指标列为约束条件。3、 约束条件的确定在选取约束条件时应当避免出现相互矛盾的约束。因为相互矛盾的约束必然导致可行域为一空集，使问题的解不存在。另外应当尽量减少不必要的约束。不必要的约束不仅增加优化设计的计算量，而且可能使可行域缩小，影响优化结果。二、 优化实例机床主轴是机床中重要零件之一，一般为多支承空心阶梯轴。为了便于使用材料力学公式进行结构分析，常将阶梯轴简化成以当量直径表示的等截面轴。下面以两支承主轴为例，说明优化设计的全过程。右图所示的是一个已经简化的机床主轴。已知主轴内径d=30mm，外力F=15000N，许用挠度y0=0.05mm。主轴材料是铸钢。密度=7.810-6Kg/ mm3，弹性模量E=210GPa设计变量数n=3，约束函数个数m=5，设计变量的初值、上下限列于表8-1中。表8-1 初始数据设计变量X1X2X3初始值480100120下限值3006090上限值650140150三、 实例分析-数学模型的建立在设计这根主轴时，有两个重要因素需要考虑。一是主轴的自重；一是主轴伸出端的挠度。对于普通机床，并不追求过高的加工精度，对机床主轴的优化设计，以选取主轴的自重最轻为目标，外伸端的挠度是约束条件。（1） 设计变量的确定当主轴的材料选定时，其设计方案由四个设计变量决定。即孔径d、外径D、跨距l及外伸端长度a。由于机床主轴内孔常用于通过待加工的棒料，其大小由机床型号决定，不能作为设计变量。故设计变量 取为x=x1 x2 x3T=l D aT（2） 目标函数的确定机床主轴优化设计的目标函数则为f(x)=14x1+x3(x22-d2)式中，-材料的密度。（3） 约束条件的确定主轴刚度是一个重要性能指标，其外伸端的挠度y不得超过规定值y0，据此建立性能约束g(x)=y-y00在外力F给定的情况下，y是设计变量x的函数，其值按下式计算y=Fa2(l+a)3EI式中，I=64(D4-d4)则 g(x)=64Fx32(x1+x3)3E(x24-d4）-y00此外，通常还应考虑主轴内最大应力不得超过许用应力。由于机床主轴对刚度要求比较高，当满足要求时，强度尚有相当富裕，因此应力约束条件可不考虑。边界约束条件为设计变量的取值范围，即lminllmaxDminDDmaxaminaamax综上所述，将所有约束函数规格化，主轴优化设计的数学模型可表示为minf(x)=14x1+x3(x22-d2) g1(x)=64Fx32(x1+x3)3E(x24-d4）/y0-10g2(x)=1-x1/lmin0 g3(x)=1-x2/Dmin0 g4(x)=x2/Dmax-10g5(x)=1-x3/amin0这里未考虑两个边界约束，x1lmax和x3amax，这是因为无论从减小伸出端挠度上看，都要求主轴跨距x1、伸出端长度x3往小处变化，所以对其上限可以不作限制。这样可以减少一些不必要的约束，有利于优化计算。四、 优化方法的选择该实例中，设计变量数n=3，约束函数个数m=5，且有一个约束函数是非线性的，故属于非线性多变量约束优化问题。可以选择MATLAB中的fmincon()函数。因为fmincon函数是实现多变量约束优化， 解决非线性多变量约束问题的 一种优化方法。当然也可以选择遗传算法ga()。五、 程序及编译结果1.编写一个M文件返回目标函数在点x处的值： function f=my1(x)p=7.8*10(-6);d=30;f=0.25*pi*p*(x(1)+x(3)*(x(2)2-d2);2.将非线性约束编进一个M文件内function c,ceq=nonlcon1(x)F=15000;d=30;E=2.1*105;c=(64*F*x(3)2*(x(1)+x(3)/(3*E*pi*(x(2)4-d4)*0.05)-1;ceq=;3.（1）调用fmincon（）函数： lb=300 60 90; %设计变量的下限ub=650 140 150;%设计变量的上限x0=480 100 120;%设计变量的初始值x,f,exitflag,output=fmincon(my1,x0,lb,ub,nonlcon1) %调用fmincon（），返回exitflag参数，描述函数计算的退出条件；返回output机构数组，其中包含了优化信息Warning: Trust-region-reflective algorithm does not solve this type of problem, using active-set algorithm.You could also try the interior-point or sqp algorithms: set the Algorithm option to interior-point orsqp and rerun. For more help, see Choosing the Algorithm in the documentation. In fmincon at 472Local minimum found that satisfies the constraints.Optimization completed because the objective function is non-decreasing in feasible directions, to within the default value of the function tolerance,and constraints were satisfied to within the default value of the constraint tolerance.Active inequalities (to within options.TolCon = 1e-006): lower upper ineqlin ineqnonlin 1 1 3 x = 300.0000 74.8898 90.0000f = 11.2494exitflag = 1output = iterations: 11 funcCount: 44 lssteplength: 1 stepsize: 4.8867e-007 algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorderopt: 1.7487e-007 constrviolation: 2.6791e-008 message: 1x788 char(2)调用ga（）函数： lb=300 60 90;%设计变量的下限ub=650 140 150;%设计变量的上限x,f,exitflag,output=ga(my1,3,lb,ub,nonlcon1)%调用ga（），返回exitflag参数，描述函数计算的退出条件；返回output机构数组，其中包含了优化信息Optimization terminated: average change in the fitness value less than options.TolFun and constraint violation is less than options.TolCon.x = 300.0002 74.8904 90.0006f = 11.2497exitflag = 1output = problemtype: nonlinearconstr rngstate: 1x1 struct generations: 4 funccount: 5180 message: 1x140 charmaxconstraint: 0结果分析：调用fmincon()函数和遗传算法函数ga()都得到了最优解，结果近似，本次优化设计可认为达到了目的。六、进一步的考虑上述主轴优化设计中是把阶梯轴简化成当量直径的等截面轴进行结构分析的，这只是一种近似分析方法，而其近似程度往往不能令人满意。尤其是对一些受力、形状和支承都比较复杂的轴，不可能做出那样的简化，况且机床主轴的设计还对其动力学性能提出一定的要求。因此，将主轴简化后用材料力学公式进行分析的方法也不能满足工程设计的需要。图8-3所示的机床主轴为三支承系统，受有力和力矩的作用。对其进行重量最轻结构优化设计时，不仅对伸出端点的挠度有要求，而且对主轴系统的第一阶段自振频率也有要求。对于这样复杂的系统，材料力学分析方法已显得无能为力了。这时常使用有限元法来计算系统的应力、变形、自振频率等。研究怎样以尽可能少的有限元分析次数而获得优化结果，是结构优化设计研究中的一个重要课题。七、 心得通过本次试验，我有了优化设计的初步概念，并进行了尝试，掌握了建立数学模型的方法，并能根据不同的数学模型合理选择优化方法，最终能找到最优解，达到优化目的。同时，也基本掌握了MATLAB编程。这一过程也锻炼了我解决问题、分析问题的能力，受益匪浅。