- 1、本文档共11页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
优化设计有约束优化无约束优化
目 录
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利用matlab工具箱fminunc函数对该模型进行求解 - 11 -
2.3.3求解结果 - 11 -
1.多维有约束优化
1.1 题目
对一对单级圆柱齿轮减速器,以体积最小为目标进行多维有约束优化设计。
1.2 已知条件
已知数输入功p=58kw,输入转速n1=1000r/min,齿数比u=5,齿轮的许用应力[]H=550Mpa,许用弯曲应力[]F=400Mpa。
1.3 建立优化模型
1.3.1问题分析及设计变量的确定
由已知条件得求在满足零件刚度和强度条件下,使减速器体积最小的各项设计参数。由于齿轮和轴的尺寸(即壳体内的零件)是决定减速器体积的依据,故可按它们的体积之和最小的原则建立目标函数。
单机圆柱齿轮减速器的齿轮和轴的体积可近似的表示为:
式中符号意义由结构图给出,其计算公式为
由上式知,齿数比给定之后,体积取决于b、z1 、m、l、dz1 和dz2 六个参数,则设计变量可取为
1.3.2目标函数的确定
根据以上分析,可知,该齿轮减速器以体积最小的目标函数为:
1.3.3 约束条件的建立
(1)为避免发生根切,应有,得
(2)齿宽应满足,和为齿宽系数的最大值和最小值,一般取=0.9,=1.4,得:
(3)动力传递的齿轮模数应大于2mm,得
(4)为了限制大齿轮的直径不至过大,小齿轮的直径不能大于,得
(5)齿轮轴直径的范围:得
(6)轴的支撑距离按结构关系,应满足条件:(可取=20),得
(7)齿轮的接触应力和弯曲应力应不大于许用值,得
(8)齿轮轴的最大挠度不大于许用值,得
(9)齿轮轴的弯曲应力不大于许用值,得
1.4 优化方法的选择
由于该问题有6个设计变量,16个约束条件的优化设计问题,采用传统的优化设计方法比较繁琐,比较复杂,所以选用Matlab优化工具箱中的fmincon函数来求解此非线性优化问题,避免了较为繁重的计算过程。
1.5 数学模型的求解
1.5.1 确定数学优化模型
将已知及数据代入上式,该优化设计的数学优化模型表示为:
(1)求变量:
(2)目标函数:
(3)约束条件:
1.5.2运用Matlab优化工具箱对数学模型求解
(1)首先在Matlab优化工具箱中编写目标函数的M文件 myfun.m,返回x处的函数值f:
function f = myfun(x)
f=0.785398*(4.75*x(1)*x(2)^2*x(3)^2+85*x(1)*x(2)*x(3)^2-85*x(1)*x(3)^2+0.92*x(1)*x(6)^2-x(1)*x(5)^2+0.8*x(1)*x(2)*x(3)*x(6)-1.6*x(1)*x(3)*x(6)+x(4)*x(5)^2+x(4)*x(6)^2+28*x(5)^2+32*x(6)^2)
(2)由于约束条件中有非线性约束,故需要编写一个描述非线性约束条件的M文件mycon.m:
function[c,ceq]=myobj(x)
c=[17-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.5*x(6)-x(4)-40;1486250/(x(2)*x(3)*sqrt(x(1)))-550; 7098/(x(1)*x(2)*x(3)^2*(0.169+0.006666*x(2)-0.0000854*x(2)^2))-400;7098/(x(1)*x(2)*x(3)^2*(0.2824+0.00177*x(2)-0.0000394*x(2)^2))-400;117.04*x(4)^4/(x(2)*x(3)*x(5)^4)-0.003*x(4);(1/(x(5)^3))*sqrt((2850000*x(4)/(x(2)*x(3)))
文档评论(0)