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

  1. 1、本文档共11页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 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)

xy88118 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档