优化设计 有约束优化 无约束优化

《优化设计 有约束优化 无约束优化》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、目录1.多维有约束优化-3-1.1题目-3-1.2已知条件-3-1.3建立优化模型-3-1.3.1问题分析及设计变量的确定-3-1.3.2目标函数的确定-4-1.3.3约束条件的建立-4-1.4优化方法的选择-5-1.5数学模型的求解-5-1.5.1确定数学优化模型-5-1.5.2运用Matlab优化工具箱对数学模型求解-6-1.5.3最优解以及结果分析-7-2.多维无约束优化-8-2.1题目-8-2.2确定优化设计模型-8-2.3运用Matlab优化工具箱对数学模型求解-9-2.3.1编写目标函数-9-2.3.2绘制该函数的平面和空间等值线-9-2.3.3利用matla

2、b工具箱fminunc函数对该模型进行求解-11-2.3.3求解结果-11--11-1.多维有约束优化1.1题目对一对单级圆柱齿轮减速器，以体积最小为目标进行多维有约束优化设计。1.2已知条件已知数输入功p=58kw，输入转速n1=1000r/min，齿数比u=5，齿轮的许用应力[]H=550Mpa，许用弯曲应力[]F=400Mpa。1.3建立优化模型1.3.1问题分析及设计变量的确定由已知条件得求在满足零件刚度和强度条件下，使减速器体积最小的各项设计参数。由于齿轮和轴的尺寸（即壳体内的零件）是决定减速器体积的依据，故可按它们的体积之和最小的原则建立目标函数。单机圆柱齿轮

3、减速器的齿轮和轴的体积可近似的表示为：式中符号意义由结构图给出，其计算公式为由上式知，齿数比给定之后，体积取决于b、z1、m、l、dz1和dz2六个参数，则设计变量可取为-11-1.3.2目标函数的确定根据以上分析，可知，该齿轮减速器以体积最小的目标函数为：1.3.3约束条件的建立(1)为避免发生根切，应有，得(2)齿宽应满足，和为齿宽系数的最大值和最小值，一般取=0.9，=1.4，得:(3)动力传递的齿轮模数应大于2mm，得(4)为了限制大齿轮的直径不至过大，小齿轮的直径不能大于，得(5)齿轮轴直径的范围：得(6)轴的支撑距离按结构关系，应满足条件：（可取=20），得(

4、7)齿轮的接触应力和弯曲应力应不大于许用值，得-11-(8)齿轮轴的最大挠度不大于许用值，得(9)齿轮轴的弯曲应力不大于许用值，得1.4优化方法的选择由于该问题有6个设计变量，16个约束条件的优化设计问题，采用传统的优化设计方法比较繁琐，比较复杂，所以选用Matlab优化工具箱中的fmincon函数来求解此非线性优化问题，避免了较为繁重的计算过程。1.5数学模型的求解1.5.1确定数学优化模型将已知及数据代入上式，该优化设计的数学优化模型表示为：（1）求变量：（2）目标函数：（3）约束条件：-11-1.5.2运用Matlab优化工具箱对数学模型求解（1）首先在Matlab

5、优化工具箱中编写目标函数的M文件myfun.m，返回x处的函数值f：functionf=myfun(x)f=0.785398(4.75x(1)x(2)^2x(3)^2+85x(1)x(2)x(3)^2-85x(1)x(3)^2+0.92x(1)x(6)^2-x(1)x(5)^2+0.8x(1)x(2)x(3)x(6)-1.6x(1)x(3)x(6)+x(4)x(5)^2+x(4)x(6)^2+28x(5)^2+32x(6)^2)（2）由于约束条件中有非线性约束，故需要编写一个描述非线性约束条件的M文件mycon.m：function[c,ceq]=myobj(x)c=[1

6、7-x(2);0.9-x(1)/(x(2)x(3));x(1)/(x(2)x(3))-1.4;2-x(3);x(2)x(3)-300;100-x(5);x(5)-150;130-x(6);x(6)-200;x(1)+0.5x(6)-x(4)-40;1486250/(x(2)x(3)sqrt(x(1)))-550;7098/(x(1)x(2)x(3)^2(0.169+0.006666x(2)-0.0000854x(2)^2))-400;7098/(x(1)x(2)x(3)^2(0.2824+0.00177x(2)-0.0000394x(2)^2))-400;117.04x(

7、4)^4/(x(2)x(3)x(5)^4)-0.003x(4);(1/(x(5)^3))sqrt((2850000x(4)/(x(2)x(3)))^2+2.410^12)-5.5;(1/(x(6)^3))sqrt((2850000x(4)/(x(2)x(3)))^2+610^13)-5.5];ceq=[];-11-（3）最后求解，调用目标函数和约束条件，用matlab软件中工具箱里的fmincon函数，求解有约束的优化，在commandwindow里输入：x0=[230;21;8;420;120;160];%给定初始值[x,f