机械优化大作业-MATLAB复合型法.doc

__________________________________________________ __________________________________________________ __________________________________________________ __________________________________________________ __________________________________________________ __________________________________________________ 机床主轴机构优化设计 机床主轴是机床中重要的零件之一，一般多为支撑空心阶梯轴。为了便于使用材料力学进行结构分析，常常将阶梯轴简化成以当量直径表示的等截面轴。 如图0.1所示是一个已经简化处理的专用机床双主轴的力学模型。从机床主轴制造成本较低和加工精度较高的要求出发，需要考虑主轴的自重和外身段挠度这样两个重要因素。对于专用机床来说，并不追求过高的加工精度。因此选取零件自身重量轻为设计目标，将主轴的刚度作为约束条件。 设计变量和目标函数 当主轴的材料选定之后，与主轴重量设计方案有关设计变量主轴的外径D、孔径d、两支承跨度l、和外伸段长度a，如图0.1所示。由于机床主轴的孔径主要取决于待加工棒料的直径，不能作为设计变量处理。因此，设计变量为 机床主轴重量最轻优化设计目标函数为 约束条件 刚度性条件 机床的加工质量在很大程度上取决于主轴的刚度，主轴刚度是一个很重要的性能指标。因此，要求主轴悬臂端挠度不超过给定的变形量y0 。 根据材料力学可知，主轴悬臂端挠度为 式中，是空心主轴的惯性矩；E=2.1x105MPa是主轴的弹性模量；F是作用主轴外伸端的力。 整理得到主轴刚度的约束条件为 2.设计边界条件 三个设计变量的边界约束条件为 3.具体算例 已知某机床主轴悬臂端受到的切削力F=15000N，主轴内径d=30mm，悬臂端许用挠度y0=0.05mm。要求主轴两支承跨度，外径，悬臂长度。 建立优化的数学模型 4.基于matlab程序的复合形法实现优化 （1）算法原理 复合形法来源于无约束优化问题的单纯形法，通过构造复合形来求得最优解，新的复合形通过替换旧的复合形中的坏点（目标函数最大或次大的点）得到，替换方式仍是单纯性法中的反应、压缩、扩展这几个基本方法。 （2）算法步骤 其程序框图如下： 是否否否否是是是是否求反射点计算各顶点的目标函数值f（xj）（j=1,2，……，k）一次坏点xG代替最坏点xH结 束xR可行？xC可行？？ 是 否 否 否 否 是 是 是 是 否 求反射点 计算各顶点的目标函数值 f（xj）（j=1,2，……，k） 一次坏点xG代替最坏点xH 结 束 xR可行？ xC可行？ ？ 计算出去xH后的各顶点中心将各定点的目标函数值和坐标按目标函数值的大小排序形成初始复合形的k个顶点xj（j=1,2，…… 计算出去xH后的各顶点中心 将各定点的目标函数值和坐 标按目标函数值的大小排序 形成初始复合形的k个顶点 xj（j=1,2，……，k） 输入 n，k，ε 开 始 （3）算法的MATLAB实现 在MATLAB中利用编程语句编写实现复合形法的函数，命名：jiangfun.m 功能：用复合形法求解约束多为函数的极值 调用格式：[x,minf]= jiangfun(f,g,X,alpha,sita,gama,beta,var,eps) 其中符号表示如下： f g X alpha gama beta 目标函数 约束函数 初始复合形 反射系数 扩展系数 收缩系数 var eps x minf 自变量向量 自变量精度 目标变量最小值时的自变量值 目标函数的最小值 由目标函数级约束条件，取初始复合形 取α（反射系数）=1.2，β（收缩系数）=0.3，θ（紧缩系数）=0.5，γ（扩展系数）=2。 在MATLAB中输入下列命令 syms s t r f=0.7854*(s+r)*(t^2-900); g=[1-9.7009*r^2*(s+r)/(t^4-30^4);s/300-1;1-s/650;t/60-1;1-t/140;r/90-1;1-r/150]; X=[325 380 500 600;90 120 100 135;100 95 120 130]; [x,minf]=jiangfun(f,g,X,1.2,0.5,8,0.3,[s t r]) 所得结果为： x =300.0357 75.2442 90.0013 Minf = 1458663