资源描述
利用MATLAB求解机械设计优化问题的分析周婷婷(能源与动力学院,油气0701)摘要:MATLAB是目前国际上最流行的科学与工程计算的软件工具,它具有强大的数值分析、矩阵运算、信号处理、图形显示、模拟仿真和最优化设计等功能。本文浅谈MATLABS机械设计优化问题的几点应用。关键词:MATLAB约束条件机械设计优化引言:在线性规划和非线性规划等领域经常遇到求函数极值等最优化问题,当函数或约束条件复杂到一定程度时就无法求解,而只能求助于极值分析算法,如果借助计算器进行手工计算的话,计算量会很大,如果要求遇到求解极值问题的每个人都去用BASIC,C和FORTRAN之类的高级语言编写一套程序的话,那是非一朝一日可以解决的,但如用MATLAB语言实现极值问题的数值解算,就可以避免计算量过大和编程难的两大难题,可以轻松高效地得到极值问题的数值解,而且可以达到足够的精度。1无约束条件的极值问题的解算方法设有Rosenbrock函数如下:f(X1,X2)=100(X2-X1*X1)2+(1-X1)2求向量X取何值时,F(x)的值最小及最小值是多少?先用MATLAB语言的编辑器编写求解该问题的程序如下:%把函数写成MATLAB语言表达式fun=100*(X(2)-X(1)*X(1)2+(1-X(1)2%猜自变量的初值X0=-12;%所有选项取默认值options=;%调用最优化函数进行计算。%函数最小值存放在数组元素options(8)中%与极值点对应的自变量值存放在向量X里%计算步数存放在数组元素options(10)中X,options=fmins(fun,X0,options);%显示与极值点对应的自变向量X的值。%显示函数最小值options(8)%显示函数计算步数options(10)把上面这段程序保存为m文件,然后用“Tools”菜单中的“Run”命令行这段程序,就可以轻松的得到如下结果:X=9.999908938395383e-0019.99982742178110e-001ans=1.706171071794760e-001ans=195显然,计算结果与理论结果的误差小到e-10级,这里调用了MATLAB的最优化函数fmins(),它采用Nelder-Mead的单纯形算法,就是因为这个函数的采用,使最小值问题的解算变得非常简单。2.带约束条件的极值问题的解法设目标函数和约束条件如下:f(x)=-3X1+X2+X3-X1+2X2-X3=-114X1-X2-2X3=0,X2=0,X3=0;求X向量取何值时函数取极小值?对条件极值问题通常的做法都是将约束条件标准化(即把等式约束条件写成等号为0的形式,把不等式写成=0的形式)。然后把条件极值问题转换为非条件极值问题,MATLAB也采用同样的做法。面是求解该问题的MATLAB语言程序。funf=f=-3*X(1)+X(2)+X(3);%写出目标函数表达式。fung=g=2*X(1)-X(3)+1;X(1)-2*X(2)+X(3)-11;4*X(1)-X(2)-2*X(3)+3;%把约束条件标准化,写成向量函数。注意等式约束条件要放在前面。fun=funf,fung%把目标函数表达式和约束条件表达式合成一个向量函数。X0=101%猜初值options=options13=1%在options(13)中指定等式约束条件的个数。vlb=zeros(1,3);%指定向量X的下界为0。vub=;%对向量的上界不做要求。X,options=constr(fun,X0,options,Vlb,Vlb);%调constr()函数显示与极值点对应的X向量值。options(8)%显示极小值options(10)%显示计算次数g=2*X(1)-X(3)+1;X(1)-2*X(2)+X(3)-11;4*X(1)-X(2)-2*X(3)+3%显示约束条件表达式的取值。运行这段MATLAB程序得到如下结果:X=4.000000000000000e+0001.00000000000000le+0009.000000000000005e+000ans=-1.99999999999999le+000ans=17g=-6.21724522790877e-0151.77635683940025le-015-1.421085471520200e-014显然,计算结果是非常精确的,这里调用了MATLAB的最优化函数constr(),它是专门用来解算条件极值问题的3.机械优化设计应用实例机械优化设计把数学规划理论与数值方法应用于设计中,用计算机从大量可行方案中找出最优化设计方案,从而大大提高设计质量和设计效率。MATLAB具有解决线性规划和非线性规划、约束优化和无约束优化问题的内部函数,因而可以完成这一功能。现举一例:螺栓组联结的优化设计如图4所示的压力容器螺栓组联接中,已知Di=400mm,D2=250mm,缸内工作压力为p=1.5MPa,螺栓材料为35号钢,cs=320Mpa,安全系数S=3,取残余预紧力Qp=1.6F,采用铜皮石棉密封垫片。现从安全、可靠、经济的角度来选择螺栓的个数n和螺栓的直径d。3.1设计问题分析若从经济性考虑,螺栓数量尽量少些、尺寸小些,但这会使降低联结的强度和密封性,不能保证安全可靠的工作;若从安全、可靠度考虑,螺栓数量应多一些、尺寸大一些为好,显然经济性差,甚至造成安装扳手空间过小,操作困难。为此,该问题的设计思想是:在追求螺栓组联结经济成本最小化的同时,还要保证联结工作安全、可靠。3.2设计变量目标函数约束条件3.2.1设计变量选取螺栓的个数n和直径d(mm)为设计变量:X=ndT=X1X2t3.22目标函数追求螺栓组联结经济成本G最小为目标。而当螺栓的长度、材料和加工条件一定时,螺栓的总成本与nd值成正比,所以本问题优化设计的目标函数为minF(X)=Cn=nd=X1X2强度约束条件为了保证安全可靠地工作,螺栓组联结必须满足强度条件cacr320其中十加3.06Mpa;Q=QpF=1.6FF=2.6F=2.6屯4np=2.61.5-25060937二.n4nN;对于粗牙普通螺纹:由文献3推荐,小径d1=0.85d所以,强度约束条件为: 密圭寸约束条件为了保证密圭寸安全,螺栓间距应小于10d,所以,密封约束条件为:g2(X)二7:Di-10dn4二-10X20X1x-xfgi(X-10-10-100nd1nd安装扳手空间约束条件为了保证足够的扳手空间,螺栓间距应大于5d,二Di所以,安装约束条件为:g3(X)=5d-n5x?-叫0X1边界约束条件g4(X)二一为乞033.3建立数学模型综上所述,本问题的数学模型可表达为:设计变量:X=X1X2】T目标函数:minf(x)=X1X2约束条件:s.t.gj(X)岂0(i现运用MATLAB的优化函数进行求解先编写M文件functionc,ceq=mynas(x)c(1)=146192/(x(1)*x(2)9)-106;c(2)=400*pi/x(1)-10*x(2);c(3)=-400*pi/x(1)+5*x(2);ceq=;在MATLAB命令窗口输入:fun=x(1)*x(2);x0=4,6;A=-1,0;0,-1;b=0;0;Aeq=;beq=;lb=;ub=;=1,2,3,4,5,)%非线性不等式约束%非线性等式约束%目标函数%设计变量初始值%线性不等式约束矩阵%线性等式约束矩阵%边界约束矩阵图5数学模型要素图x,fval=fmincon(fun,xO,A,b,Aeq,beq,lb,ub,mynIsub)%调用有约束优化函数运行结果如下:X=11.449910.9751fval=125.6637所以,该问题优化结果为:n=11.4499,d=10.9751,目标函数最小值:F(X)=125.6637。根据实际问题的意义取整、标准化:n=12,d=12由此例可以看出,与其它编程语言相比,MATLAB语言可以简化编程。图5是调用MATLAB绘图函数自动对上例绘制的数学模型要素图(标注数字的曲线为目标函数的等值线),为此在MATLAB命令窗口输入:x1=0.1:20;y1=146192./(106.*x1.A2);y2=400.*pi./(10.*x1);y3=400.*pi./(5.*x1);Plot(y1,x1,y2,x1,y3,x1,x(1),x(2),o)y4=0.1:0.1:20;y4,x1=meshgrid(y4,x1);Q=y4.*x1;holdon;c,h=contour(y4,x1,Q);holdon;clabel(c,h);4结束语从上述实例可以看出,利用求解最优化问题具有编程简单,精度很高,速度很快,各种工形式的最优化问题都适用等优点,巧妙各种利用MATLAB语言可以取得事半功倍的效果。MATLAB具有科学计算的强大能力,不管处理什么样的对象算法、图形、图像、报告或者算法仿真MATLAB都能够帮助大家不断提高工作效率。MATLAB环境能够完成算法开发、数据分析和可视化、高性能数据分析等工作,相对于传统的C、C+或者FORTRAN语言,在MATLAB中完成这些工作所消耗的时间仅仅是传统手段的极小一部分。开放性的MATLAB软件和开发语言M语言,能够快速实现用户好的想法和概念。MATLAB本身包含了600余个用于数学计算、统计和工程处理的函数,这样,就可以迅速完成科学计算任务而不必进行额外的开发。业内领先的工具箱算法极大的扩展了MATLAB的应用领域,例如信号处理,数字图像处理、数据分析和统计以及算法建模和仿真等。MATLAB的工具箱都是业内的专家、工程师结合多年来的经验和专业知识,专门开发的用户处理特殊数学计算、分析和图形可视化的功能函数集合。利用这些工具箱,用户可以通过对比应用一系列不同的算法,而避免了应用程序的开发。使用MATLAB软件进行科学计算,能够极大加快科研人员进行研究开发的进度,减少在编写程序和开发算法方面所消耗的时间和有限的经费,从而获得最大的效能。参考文献1薛定宇控制系统计算机辅助设计-MATLAB语言及应用,清华大学出版社,19962濮良贵机械零件北京.高等教育出版社。3蒲俊,吉家锋MATLAB数学手册上海.浦东电子出版社,2002。4黄华梁彭文生机械设计基础高等教育出版社利用MATLAB求解机械设计优化问题的分析学院:能源与动力工程学院班级:油气0701姓名:周婷婷学号:0120705690137
展开阅读全文