资源描述：
《机械优化设计综述.doc》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、关于MATLAB的机械优化设计综述摘要：机械优化设计是以数学规划理论为基础，以计算机为工具，一种自动寻优、先进的、现代的设计方法。MATLAB作为一种集计算功能、符号运算功能和图形处理功能于一身的科学计算语言，其中的优化工具箱，可为线性、非线性最小化、非线性最小二乘、二次规划、方程求解、多目标优化等问题提供求解方法。运用MATLAB解决机械设计的优化问题不仅方便实用，且编程简单，效率高。关键词：机械优化，MATLAB正文：线性规划在机械设计中应用非常少，主要是因其要求目标函数、约束条件都是线性的，而机械设计问题一般都是非线性的。目前，线性规划
2、在一维下料问题中应用比较多。求解线性规划问题常用的方法有单纯形法、大M法等。在MATLAB屮由linprog函数来求解线性规划问题。二次规划在机械设计中的应用比较少，在MATLAB中求解二次规划问题可用quadprog函数。无约束优化有很多种算法，如最速卜•降法＞Newton法、拟Newton法、共轨梯度法、信赖域法等。其中，拟Newton法利用了Newton法的二次收敛性使得可靠性更高，同时这种算法收敛速度也很快，所以在目前应用得相当广泛。而拟Newton法中应用比较多的主要有BFGS算法和DFP变尺度法。在MATLAB优化工具箱中求解无约
3、束优化问题釆用fminbnd函数、fminunc函数和fminsearch函数，其中默认的算法就是BFGS算法ofminunc函数和fminbnd函数要求目标函数必须连续，而函数fminsearch常用来处理不连续的函数□fminbnd函数可以求解区间：*1,*2］内单变量非线性函数的最小值。在机械优化设计中大多数问题是有约束优化问题。为了保证设计的机械零件安全可靠，设计时必须确定相应的设计准则，而这些设计准则就成了目标函数的约束条件。这些设计准则包括强度准则、刚度准则、寿命准则、振动稳定性准则和可靠性准则，大多数是非线性的。对于非线性约束优
4、化问题有很多解法，早期都是通过构造惩罚函数将有约束优化传换为无约束优化问题求解，而现在已经被更为有效的基于K-T方程解的方法所取代。在MATLAB优化工具箱中采用序列二次规划法（SQP）求解约束优化问题，采用的函数为fmincom函数。对于大量的工程设计方案要评价其优劣，往往需要同时考虑多个目标。如在机械加T中，往往要求机械加丁-成本最低、牛•产率最高和刀具寿命最长。显然，这样的问题就属于多目标优化问题。多目标优化问题的求解方法很多，其中一种常用的方法就是将多冃标优化问题新构造一个函数，即评价函数，从而将多目标优化问题转化为求评价函数的单冃标
5、优化问题。这类方法主要有：主要冃标法、线性加权和法、理想点法和极大极小法等。另一•种是将多目标优化问题转化为一系列单目标优化问题来求解，这一类主耍有分层序列法。在MATLAB中常用fminimax函数求解最大最小化问题，用fgoalat—tain函数求解多目标达到问题。实例：通过课本上人字架的优化设计例子，参考一些论文，运用MATLAB软件对这个二维非线性优化问题进行求解：选择设计变量：以人字架高度h和钢管平均直径D作为设计变量；确定冃标函数：使人字架结构的质量m最小，确定冃标函数；确定约束条件：通过强度条件和稳定约束条件建立；建立数学模型：
6、用MATLAB进行优化求解：用MATLAB求解步骤如下：(1)冃标函数分析：[x,y]=meshgrid(linspace(0,17,19)Jinspace(0,11/19))；B=76;T=0.25;rou=7.8e-3;m=2*pi*rou*T*x*sqrt(BA2+yA2)；mesh(x,y,m);holdon；xlabel(,D(cm),);ylabel(,h(cm),)；zlabel(,m(kg)')；titled标函数图像分析J目标函数图像分析••161.5丫・:.•二・・・・••・161•••••・■:1...t叫..・]・・•
7、•h(cm)159.5—•]D(cm)目标函数是人字架的高h和钢管平均直径D的函数，并且目标函数图像规则。(2)约束函数分析：»ezplot(,((1.5e4*sqrt(76A2+yA2))/(pi*0.25*x*y))-4200,/[0/12/0,120])»holdon;»y=0:0.01:120;x=2;»plot(x,y/k・')»holdon»y=0:0.01:120;x=10;»holdon»ezplotC((1.5e4*sqrt(76A2+yA2))/(pi*0.25*x*y))-((piA2*2.1e6*(0.25A2+xA
8、2))/(8*(76A2+yA2)))[0,12,0,120])»holdon»x=0:0.01:15;y=30;»plot(x,y；k-')»holdon»x=
