机械优化设计实例及matlab优化工具课件

1、机械优化设计实例 及matlab优化工具 机械优化设计实例 机械优化设计的一般过程 建立数学模型的基本原则 机械优化设计实例 机械优化设计全过程一般可分为：1)建立优化设计的数学模型。2)选择适当的优化方法。3)编写计算机程序。4)准备必要的初始数据并上机计算。5)对计算机求得的结果进行必要的分析。 机械优化设计的一般过程1）设计变量的选择： 尽量减少设计变量数目 设计变量应当相互独立2）目标函数的确定： 选择最重要指标作为设计追求目标3）约束条件的确定： 性能约束和边界约束建立数学模型的基本原则试设计一重量最轻的空心传动轴。空心传动轴的D、d分别为轴的外径和内径。轴的长度不得小于3m。轴的材

2、料为45钢，密度为7.810-6/，弹性模量E=2105MPa，许用切应力=60MPa。轴所受扭矩为M=1.5106Nmm。 设计实例1:分析设计变量：外径D、内径d、长度l设计要求：满足强度，稳定性和结构尺寸要求外，还应达到重量最轻目的。所设计的空心传动轴应满足以下条件：（1）扭转强度 空心传动轴的扭转切应力不得超过许用值，即 空心传动轴的扭转切应力: 设计实例1:空心传动轴的扭切应力: 设计实例1:经整理得: （2）抗皱稳定性 扭转切应力不得超过扭转稳定得临界切应力:空心传动轴的扭转稳定的临界切应力为: 设计实例1:整理得: 设计实例1:（3）结构尺寸 设： 则数学模型为： 设计实例1:

3、设计实例1:平面连杆机构优化设计一曲柄摇杆机构，M为连秆BC上一点，mm为预期的运动轨迹，要求设计该曲柄摇杆机构的有关参数，使连杆上点M在曲柄转动一周中，其运动轨迹(即连杆曲线)MM最佳地逼近预期轨迹mm。 设计实例2: 设计一再现预期轨迹mm的曲柄摇杆机构。已知xA67mm，yA10mm，等分数s12，对应的轨迹mm上12个点的坐标值见表，许用传动角300。设计实例2:一、建立优化设计的数学模型 点M的坐标: 设计实例2:点M的坐标: 设计实例2:该问题有8个设计变量，记为： 设计实例2:2）确定目标函数将曲柄一周转角分为s等分，要求连秆曲线最佳地逼近预期轨迹mm，具体可由连杆曲线上的s个点

4、M最佳地逼近预期轨迹上的s个点m予以实现。由此可按点距和最小的原则建立如下目标函数：设计实例2:3）确定约束条件(1)由曲柄存在条件，可得:(2)由杆长必须大于零及曲柄1为最短杆，可得： 设计实例2:(3)由满足传动角条件，可得： 设计实例2:优化设计工具优化设计工具第1部分 MATLAB基础第2部分 优化计算工具第1部分 MATLAB基础1.1 MATLAB环境简介1.2 数据表示1.3 数组1.4 源文件（M文件）1.1 MATLAB窗口 启动MATLAB其窗口如右1、Command Window (命令窗口)2、Workspace (工作区)1.2 数据表示 1、变量 变量用标识符表示（

5、字母打头、字母、数字、下划线组成，长度19）。可以合法出现而定义。 区分大小写字母，以当前值定义其类型。2、函数名 函数名用标识符表示。1.3 数组 行向量、列向量、矩阵1.3.1 创建数组的常用方法 1）直接列表定义数组例如： x=1 2 3 4 5 6 7 8 9 0 y=1,2,3,4,5,6,7,8,9,0 z=1,2,3,4,5;2,3,4,5,6,71.3 数组 2）域表定义数组 变量初值：增量：终值初值：终值 变量（初值：增量：终值）常数 例如： x=0:0.02:10 y=1:80 1.3 数组 1.3.2、 数组的访问（一维） 数组名 表示全体元素数组名(k) 表示第k元素数

6、组名(k1:k2) 表示第k1到k2元素1.3.3 数组运算1). 纯量与数组的算术运算 a c1 或 c1 a 其中可为、 结果为a1c1 a2c1 anc1 或c1a1 c1a2 c1an2). 数组加（减） 使两数组的对应各元素相加(减)3). 数组点乘两数组的对应元素相乘a.*b 结果：a1b1 a2b2anbn （a与b的维数必须相同）1.3.3 数组运算4). 数组点正除（右除）使两数组的对应元素正除 a./b 结果为: （a、b维数必须相同）1.3.3 数组运算1.4 源文件（M文件）分为两类： 函数文件和非函数文件 都用扩展名.M 1.4.1 函数文件（相当于子程序） 格式 f

7、unction 输出表函数名(输入表) 函数体1.4.2 非函数文件 无函数头的M文件，由若干命令和注释构成。相当于主程序如： %Filename is a sine.m x = 0:0.1:2*pi; y = sin(x); plot(x,y) 1.4.3 M-文件的操作 1.4.3 M-文件的操作 1.4.3 M-文件的操作 1.4.3 M-文件的操作 第2部分 优化计算工具2.1 线性规划优化函数2.2 无约束非线性优化函数2.3 约束优化函数MATLAB解决的线性规划问题的标准形式为：2.1 线性规划优化函数函数 linprog格式 x = linprog(f,A,b,Aeq,beq)

8、 x = linprog(f,A,b,Aeq,beq,lb,ub,x0)f是优化参数x的系数矩阵；A线性不等式约束系数阵b线性不等式约束常数向量Aeq线性等式约束系数阵Beq线性等式约束常数向量 2.1 线性规划及其优化函数应用举例 求使函数 取最小值的x值， 且满足约束条件： 2.1线性规划及其优化函数代码 f = -5; -4; -6; A = 1 -1 1;3 2 4;3 2 0; b = 20; 42; 30; lb = zeros(3,1); x,fval = linprog(f,A,b,lb)结果 x = 0.0000 15.0000 3.0000 fval = -78.00002

9、.2 无约束非线性优化函数函数fminunc格式 x = fminunc(fun,x0) x = fminunc(fun,x0,options) x,fval = fminunc() 2.2 无约束非线性优化函数应用举例 求 的最小值代码%首先编写目标的.m文件 function f=myfun(x) f=3*x(1)2+2*x(1)*x(2)+x(2)2%然后调用函数 fminunc x0=1,1; x,fval=fminunc(myfun,x0)2.2 无约束非线性优化函数结果x = 1.0e-008 * -0.7512 0.2479fval = 1.3818e-0162.3 约束优化函数

10、函数 fmincon格式 x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)x = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon,options)x,fval = fmincon()2.3 约束优化函数盖板优化实例： 目标函数：约束： 盖板优化实例盖板优化实例盖板优化实例运行结果：x = 0.6332 25.3264fval = 101.3056盖板优化实例前面空心轴的问题： clear allx0=23,19,4;options=optimset(largescale,off,display,iter,tolx,1e

11、-6);x,fval,exitflag,output=fmincon(myfun1,x0,confun1,options)function f=myfun1(x) f=6.12*(x(1)2-x(2)2)*x(3)*10e-6x = 33.7505 12.8830 3.0000fval = 0.1787exitflag = 4output = iterations: 7 funcCount: 39 stepsize: 1 algorithm: medium-scale: SQP, Quasi-Newton, line-search前面空心轴的问题：装载机优化设计:任何一个行业，一个市场，都是先

12、来的有肉吃，后来的汤都没的喝。7月-227月-22Wednesday, July 27, 2022年轻人欠缺经验，但请不要忘记：年轻是你最大的本钱。不要怕犯错，也不要畏惧挑战，你应该坚持到底，在出人头地的过程中努力再努力。00:11:5100:11:5100:117/27/2022 12:11:51 AM君子心胸开朗，思想上坦率洁净，外貌动作也显得十分舒畅安定。小人心里欲念太多，心理负担很重，就常忧虑、担心，外貌、动作也显得忐忑不安，常是坐不定，站不稳的样子。7月-2200:11:5100:11Jul-2227-Jul-22开始努力吧！在这个过程中你必须放弃很多东西，但你要明白它们都不是你最终

13、想要的，你要相信在你成功以后，总有一天它们会再回来，而且比现在更美好！。00:11:5100:11:5100:11Wednesday, July 27, 2022想赢个三回两回，三五，有点智商就行;想做个百老店，想一辈子赢，没有德商绝对不行。我们这个世界，从不会给一个伤心的落伍者颁发奖牌。7月-227月-2200:11:5100:11:51July 27, 2022没有承受困难的能力，就没有希望了。在漫长的人生旅途中，有时要苦苦撑持暗无天日的境遇；有时却风光绝项，无人能比。2022年7月27日12:11 上午7月-227月-22我知道什么是劳动:劳动是世界上一切欢乐和一切美好事情的源泉。企业不景气，问题不在员工，而在老板的管理方法不当。27 七月 202212:11:51 上午00:11:517月-22一个公司要发展迅速得力于聘用好的人才，尤其是需要聪明的人才。七月 2212:11 上午7月-2200:11July 27, 2022防民之口，甚于防川。国语。2022/7/27 0:11:5100:11:5127 July 2022与人交往的时候，多听少说。这就是，上帝为什么给我们一个嘴巴两个耳朵的原因。12:11:51 上午12:11 上午0