2022年MATLAB软件在机械优化设计中的应用 .pdf
《2022年MATLAB软件在机械优化设计中的应用 .pdf》由会员分享,可在线阅读,更多相关《2022年MATLAB软件在机械优化设计中的应用 .pdf(3页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、第 21 卷第 7 期2 0 0 4 年 7 月机械设计JOURNAL OF MACHIN E DESIGNVol.21No.7Jul.2004MAT LAB 软件在机械优化设计中的应用王春香 , 冯慧忠( 内蒙古科技大学机械工程学院, 内蒙古包头014010)摘要 :应用 MAT LA B 软件 , 对已有的优化实例进行求解, 并加以研究和分析, 着重说明如何利用MATL AB 优化工具箱中的函数fmincon, 来解决单目标多变量约束非线性优化问题的思路和方法 , 并给出了示例目标函数、 约 束函数图形 显示,以及求优的 MAT LAB 程序 。关键词 :MATL AB; 优化 ;机械设计
2、 ; 软件中图分类号 :TH122文献标识码 : A文章编号 : 1001- 2354( 2004)07- 0052- 03MATLAB 是美国 MathWorks 公司推出的一套功能强大的工程计算软件 , 它将计算与可视化集成到一个灵 活的计算环 境中,并提供了大量的内置函数, 在解决广泛 的工程问题 时, 可 以直接利用这 些函数获得 数值 解, 故被广 泛地 应用于 自动 控制、数理统计、数值分析、流体力学和机械设计等许多工程领域。鉴于目前国内关于探讨MAT LA B软件在工 程机械优化 问题和机械优化设计教学领域中应用的例子并不多见 , 加之出 于教学和研究的需要, 近年来对运用MA
3、TLAB 软件 , 如何解决 实际工程优化问题和优化设计教 学中的 问题 进行 了初步 的探 索和研 究。 限 于 篇 幅, 仅 列 举 一 简 单 实 例 , 重 点 说 明 如 何 用MA TLAB, 以及它的优化工具箱提供的功能来解决某些机械 优化问题。关于MATLAB 及其优 化工具 箱功 能的详 细情 况, 请参见文献 1 4 。1采用 M atlab 软件进行优化分析1. 1优化问题示例有一圆形 的等截 面销轴 , 一端 固定在 机架上 , 另 一端受 集中载荷 P = 10 kN 和转矩 M= 100 N. m 的 作用 , 其简化模型 如图 1 所示。由于结构的需要, 轴的长度
4、 l8 cm, 已知销轴材 料的许用弯曲应力 w= 120 MPa;许用切应 力 = 80 MPa; 允许挠度 f = 0. 01 cm; 密度= 7. 8 t/ m3; 弹性模量 E = 2105M Pa。现要求设计这根销轴, 在满 足使 用要 求的 前 提下 , 使 其质量最轻 5。图 1等截面销轴受力简图1. 1. 1数学模型的建立优化设计追求的指标为销轴质量Q 最轻 , Q 的 计算式为 :Q ( X) = Q( x1, x2) =14x21x2=14d2l= 0. 006 13d2l式中 : 设计变量X = x1, x2T= d, l T。显然 , d、 l 越小 , Q 值越小。但
5、二者的取值要受强 度、 刚 度及结构尺寸等各种使用要求的限制,将 以上各限 制条件整理 后可得 :( 1) 抗弯强度条件 :max=Pl0. 1 d3wd3-8. 33l0( 2) 抗扭强度条件 :max=M0. 2 d3d3-6. 250( 3) 刚度条件 :fmax=P l33 EJ=64Pl33E d4 f d4-0.34 l30( 4) 结构尺寸限制 :llminl-80据此 , 可写出优化设计的数学模型为:min Q( X ) = 0. 006 13d2lX = x1, x2 = d , l Ts. td3-8. 33 l0d3-6. 250d4-0. 34 l30l-80这是一个单
6、目标非线性二维约束优化问题。1. 1. 2求优结果为了便于对采用不同手段得到的优化结果进行比 较,以 及确认 MA TLA B分析结果的 可信 度,特 意选 择了 已经采 用其 它方法得到最优结果的例子(引自文献 5 , 并假定 其结果正确 ),最优方案为 : d*=4. 309 cm , l*=8 cm , Q*=0. 911 kg。1. 2采用 MATL AB 软件对优化问题示例的分析1. 2. 1约束函数的性态分析应用 MAT LA B语言编程 : x, y = meshgrid( linspace(0, 17, 19) , linspace( 0,11, 19) ) ;% 根 据函数的
7、定义划分网格区域Q= 0. 006 13*y.2* x ;mesh( x, y, Q) ;% 通过三维网格模拟目标函数图形xlabel( l ) ;ylabel( d ) ;zlabel( Q ) ;% 标注三坐标轴title(目标函数的图像)% 标注图形名称收稿日期 :2003-09-25; 修订日期 : 2004-01-30作者简介 :王春香 ( 1962-) , 女, 内蒙古包头人, 副教授 , 硕士 , 主要研究方向: 机械优化设计方法及CAD/CAM/ CAI 技术。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师
8、精心整理 - - - - - - - 第 1 页,共 3 页 - - - - - - - - - 运行程序可画出目标函数的三维图形, 如 图 2 所示 。图 2函数 Q= 0. 006 13d2l 的图像同理 ,可作如下编程 :ezplot( x3- 8. 33*y , 0, 10, 0, 20 )% 弯曲强度条件的图形hold ony= 0 0. 01 20;x= 6. 25(1/ 3) ;plot( x, y,k-)% 扭转强度条件的图形hold onezplot( x4- 0. 34*y3 , 0, 10, 0, 20 )% 刚度条件的图形hold onx= 0 0. 01 10;y=
9、8;plot( x, y,k-)% 长度边界条件的图形title(各约束函数的图像)text( 6. 5, 11, 可行域)hold off执行程序绘出各约束函数的图形, 如图 3 所示。从图 2可以看出 :目标函数的图像规则, 即性态好 ,对于多数优化方法均适用 ; 但由图 3 可行域的构成分析可知: d3-6.250( 即抗扭强度条件 ) 不是起作用约束 , 设计时可不予考虑。图 3设计变量的可行域1. 2. 2应用 MATL AB 软件的 优化工 具箱 对优 化问题 示例 求解a.取设计变量的初值为: x0= 410t首先 ,编写目标函数 的 m 文件 : Objfun. m, 返回 x
10、 处的 函数值 f 。functi on f= Objfun( x)f =0. 006 13* x ( 1)2*x( 2) ;b. 因设 计约束含 3 个非线 性约束 , 故需编写一个描述非线性约束的 m 文件 :N onLinConstr. mfuncti on c, ceq=NonLinConstr ( x)c( 1) =-x( 1) 3+ 8. 33*x( 2) ;c( 2) =-x( 1) 3+ 6. 25;c( 3) =-x( 1) 4+ 0. 34*x( 2)3;ceq= ;c. 设置线性约束的系数:A = 0-100 ;b= -80 ;d.给定 变量的初值 ,并且调用优化函数:x
11、0= 410 ;A = 0-100 ;b= -80 ;1b= zeros(2, 1) ;options= optimset( Display , iter ,LargeScale , off ) ; x, fval, exitflag, output = fmincon(Objfun , x0, A, b, , , lb, ,NonLinConstr, options)e.计算结果 :x=% 最后的优化结果:4. 05438. 000 0% d*= 4. 054 3, l*= 8 cmfval=% 优化后的最小质量0. 806 1% Qmin= 0. 806 1 kgexitflag=% 算法
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 2022年MATLAB软件在机械优化设计中的应用 2022 MATLAB 软件 机械 优化 设计 中的 应用
限制150内