机械优化设计_经典实例PPT课件

1、机械优化设计_经典实例机械优化设计实例机械优化设计实例 机械优化设计的一般过程机械优化设计的一般过程 建立数学模型的基本原则建立数学模型的基本原则 机械优化设计实例机械优化设计实例 机械优化设计全过程一般可分为：机械优化设计全过程一般可分为：1)1)建立优化设计的数学模型。建立优化设计的数学模型。2)2)选择适当的优化方法。选择适当的优化方法。3)3)编写计算机程序。编写计算机程序。4)4)准备必要的初始数据并上机计算。准备必要的初始数据并上机计算。5)5)对计算机求得的结果进行必要的分析。对计算机求得的结果进行必要的分析。 机械优化设计的一般过程机械优化设计的一般过程1）设计变量设计变量的选

2、择：的选择： 尽量减少设计变量数目尽量减少设计变量数目 设计变量应当相互独立设计变量应当相互独立2）目标函数目标函数的确定：的确定： 选择最重要指标作为设计追求目标选择最重要指标作为设计追求目标3）约束条件约束条件的确定：的确定： 性能约束和边界约束性能约束和边界约束建立数学模型的基本原则建立数学模型的基本原则试设计一重量最轻的空心传动轴。空心传动轴试设计一重量最轻的空心传动轴。空心传动轴的的D D、d d分别为轴的外径和内径。轴的长度不得分别为轴的外径和内径。轴的长度不得小于小于3m3m。轴的材料为。轴的材料为4545钢，密度为钢，密度为7.87.81010-6-6/ /，弹性模量，弹性模量

3、E=2E=210105 5MPaMPa，许用切应力，许用切应力 =60MPa=60MPa。轴所受扭矩为。轴所受扭矩为M=1.5M=1.510106 6N Nmmmm。 设计实例设计实例1: 1:分析分析设计变量：外径设计变量：外径D D、内径、内径d d、长度、长度l l设计要求：满足强度，稳定性和结构尺寸要设计要求：满足强度，稳定性和结构尺寸要求外，还应达到重量最轻目的。求外，还应达到重量最轻目的。所设计的空心传动轴应满足以下条件：所设计的空心传动轴应满足以下条件：（1 1）扭转强度）扭转强度 空心传动轴的扭转切应力不得超过许用值，即空心传动轴的扭转切应力不得超过许用值，即 空心传动轴的扭转

4、切应力空心传动轴的扭转切应力: : )(1644dDD设计实例设计实例1: 1:空心传动轴的扭切应力空心传动轴的扭切应力: : )(1644dDD01027. 1544DDd设计实例设计实例1: 1:经整理得经整理得: : 2/3)2(7 . 0DdDE（2 2）抗皱稳定性）抗皱稳定性 扭转切应力不得超过扭转稳定得临界切应力扭转切应力不得超过扭转稳定得临界切应力: :空心传动轴的扭转稳定的临界切应力为空心传动轴的扭转稳定的临界切应力为: : 设计实例设计实例1: 1:034.1542/344DdDdDD整理得整理得: : 设计实例设计实例1: 1:（3 3）结构尺寸）结构尺寸 minll 0d

5、0dD设：设： ldDxxx321则数学模型为：则数学模型为： 62210)(12. 6)(minldDf63222110)(12. 6xxx 设计实例设计实例1: 1:01027. 11027. 1)(541425441xxDDdXg034.15434.154)(2/3121424112/3442xxxxxxDdDdDDg033)(33xlXg0)(24xdXg0)(215xxdDXg 设计实例设计实例1: 1:平面连杆机构优化设计平面连杆机构优化设计一曲柄摇杆机构，一曲柄摇杆机构，M为连秆为连秆BC上一点，上一点，mm为预期的运动为预期的运动轨迹，要求设计该轨迹，要求设计该曲柄摇杆机构的有

6、曲柄摇杆机构的有关参数，使连杆上关参数，使连杆上点点M在曲柄转动一在曲柄转动一周中，其运动轨迹周中，其运动轨迹(即连杆曲线即连杆曲线)MM最最佳地逼近预期轨迹佳地逼近预期轨迹mm。 设计实例设计实例2:2: 设计一再现预期轨迹设计一再现预期轨迹mm的曲柄摇杆机构。已知的曲柄摇杆机构。已知xA67mm，yA10mm，等分数，等分数s12，对应的轨迹，对应的轨迹mm上上12个点的坐标值见表，许用传动角个点的坐标值见表，许用传动角300。设计实例设计实例2:2:一、建立优化设计的数学模型一、建立优化设计的数学模型 点点M的坐标的坐标: )sin()sin()cos()cos(5151llyyllxx

7、AMAMcossincos22cos2arccos)(14141242124124232221lllarctgl lllll lllll设计实例设计实例2:2:点点M的坐标的坐标: )sin()sin()cos()cos(5151llyyllxxAMAMcossincos22cos2arccos)(14141242124124232221lllarctgl lllll lllll出。将由设计的已知条件给其中 0设计实例设计实例2:2:该问题有该问题有8个设计变量，记为：个设计变量，记为： Txxxx821,设计实例设计实例2:2:Tlllll054321,2）确定目标函数）确定目标函数将曲柄一

8、周转角分为将曲柄一周转角分为s等分，要求连秆曲线最佳地逼近等分，要求连秆曲线最佳地逼近预期轨迹预期轨迹mm，具体可由连杆曲线上的，具体可由连杆曲线上的s个点个点M最佳地逼最佳地逼近预期轨迹上的近预期轨迹上的s个点个点m予以实现。由此可按点距和最小予以实现。由此可按点距和最小的原则建立如下目标函数：的原则建立如下目标函数：simiMimiMiyyxxxf122)()()(设计实例设计实例2:2:3）确定约束条件）确定约束条件(1)由曲柄存在条件，可得由曲柄存在条件，可得:(2)由杆长必须大于零及曲柄由杆长必须大于零及曲柄1为最短杆，可得：为最短杆，可得： 0)(0)(0)(32413423124

9、3211llllxgllllxgllllxg0)(14lexg设计实例设计实例2:2:(3)由满足传动角条件由满足传动角条件，可得：，可得： 02)(arccos180)(02)(arccos)(32214232263221423225llllllxgllllllxg。设计实例设计实例2:2:优化设计工具优化设计工具优化设计工具优化设计工具r第第1部分部分 MATLAB基础基础r第第2部分部分 优化计算工具优化计算工具第第1 1部分部分 MATLABMATLAB基础基础1.1 MATLAB环境简介环境简介1.2 数据表示数据表示1.3 数组数组1.4 源文件（源文件（M文件）文件）1.1 MA

10、TLAB1.1 MATLAB窗口窗口 启动启动MATLAB其窗口如右其窗口如右n1、Command Window (命令窗口命令窗口)n2、Workspace (工作区工作区)1.2 1.2 数据表示数据表示 1、变量 变量用标识符表示（字母打头、字母、数字、下划线组成，长度19）。可以合法出现而定义。 区分大小写字母，以当前值定义其类型。2、函数名 函数名用标识符表示。1.3 1.3 数组数组 行向量、列向量、矩阵行向量、列向量、矩阵1.3.1 1.3.1 创建数组的常用方法创建数组的常用方法 1 1）直接列表定义数组）直接列表定义数组例如：例如： x=1 2 3 4 5 6 7 8 9 0

11、 y=1,2,3,4,5,6,7,8,9,0 z=1,2,3,4,5;2,3,4,5,6,71.3 1.3 数组数组 2 2）域表定义数组）域表定义数组 变量初值：增量：终值初值：终值变量初值：增量：终值初值：终值 变量（初值：增量：终值）常数变量（初值：增量：终值）常数 例如：例如： x=0:0.02:10 x=0:0.02:10 y=1:80 y=1:80 1.3 1.3 数组数组 1.3.21.3.2、 数组的访问（一维）数组的访问（一维） 数组名数组名 表示全体元素表示全体元素数组名数组名(k) (k) 表示第表示第k k元素元素数组名数组名(k(k1 1:k:k2 2) ) 表示第表

12、示第k k1 1到到k k2 2元素元素1.3.3 1.3.3 数组运算数组运算1). 1). 纯量与数组的算术运算纯量与数组的算术运算 a ca c1 1 或或 c c1 1 a a 其中其中可为、可为、 结果为结果为 a1c1 a2c1 anc1 或或 c1a1 c1a2 c1an 2). 数组加（减）数组加（减） 使两数组的对应各元素相加使两数组的对应各元素相加(减减)3). 3). 数组点乘数组点乘两数组的对应元素相乘两数组的对应元素相乘a.*b 结果：结果： a1b1 a2b2anbn （a a与与b b的维数必须相同）的维数必须相同）1.3.3 1.3.3 数组运算数组运算4).

13、数组点正除（右除）数组点正除（右除）使两数组的对应元素正除使两数组的对应元素正除 a./b 结果为结果为: （a a、b b维数必须相同）维数必须相同）nnbababa11122111.3.3 1.3.3 数组运算数组运算1.4 1.4 源文件（源文件（MM文件）文件）分为两类：分为两类： 函数文件和非函数文件函数文件和非函数文件 都用扩展名都用扩展名.M .M 1.4.1 1.4.1 函数文件（相当于子程序）函数文件（相当于子程序） 格式格式 function 输出表输出表函数名函数名(输入表输入表) 函数体函数体1.4.2 1.4.2 非函数文件非函数文件 无函数头的无函数头的M文件，由若

14、干命令和注释构成。文件，由若干命令和注释构成。相当于主程序相当于主程序如：如： %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-文件的操作文件的操作 )5 . 75 . 7, 5 . 75 . 7()sin(2222yxyxyxz1.4.3 M-1.4.3 M-文件的操作文件的操作 1.4.3 M-1.4.3 M-文件的操作文件的操作 1.4.3 M-1.4.3 M-文件的操作文件的操作 第第2 2部分部分 优化计算工具优化计算工具2.1 2.1 线性规划优化函数线性规划优化函数2.2 2.

15、2 无约束非线性优化函数无约束非线性优化函数2.3 2.3 约束优化函数约束优化函数MATLAB解决的线性规划问题的标准形式为：min. .tsxcT0,xbAxTnnmijxxxxxaA),.,(,)(32, 10,),.,(,),.,(2121bbbbbccccTmTn且2.1 2.1 线性规划优化函数线性规划优化函数 函数函数 linprog linprog 格式格式 x = linprog(f,A,b,Aeq,beq) x = linprog(f,A,b,Aeq,beq) x = linprog(f,A,b,Aeq,beq,lb,ub,x0) x = linprog(f,A,b,Aeq

16、,beq,lb,ub,x0)f是优化参数是优化参数x的系数矩阵；的系数矩阵；A线性不等式约束系数阵线性不等式约束系数阵b线性不等式约束常数向量线性不等式约束常数向量Aeq线性等式约束系数阵线性等式约束系数阵Beq线性等式约束常数向量线性等式约束常数向量 2.1 2.1 线性规划及其优化函数线性规划及其优化函数 应用举例应用举例 求使函数求使函数 取最小值的取最小值的x x值，值， 且满足约束条件：且满足约束条件： 321645)(xxxxf20321xxx42423321xxx302321 xx0, 0, 0321xxx2.12.1线性规划及其优化函数线性规划及其优化函数代码 f = -5;

17、-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.2 2.2 无约束非线性优化函数无约束非线性优化函数 函数函数fminuncfminunc 格式格式 x = fminunc(fun,x0) x = fminunc(fun,x0) x = fminunc(fun,x0,options) x = fminunc(fun,x0,options) x,fval = fminunc

18、( x,fval = fminunc() ) 2.2 2.2 无约束非线性优化函数无约束非线性优化函数 应用举例应用举例 求求 的最小值的最小值 代码代码%首先编写目标的首先编写目标的.m.m文件文件 function f=myfun(x)function f=myfun(x) f=3 f=3* *x(1)2+2x(1)2+2* *x(1)x(1)* *x(2)+x(2)2x(2)+x(2)222212123)(xxxxxf% %然后调用函数然后调用函数 fminunc fminunc x0=1,1; x0=1,1; x,fval=fminunc(myfun,x0)x,fval=fminunc

19、(myfun,x0)2.2 2.2 无约束非线性优化函数无约束非线性优化函数 结果结果 x =x = 1.0e-008 1.0e-008 * * -0.7512 0.2479 -0.7512 0.2479fval =fval = 1.3818e-016 1.3818e-0162.3 2.3 约束优化函数约束优化函数函数 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 2.3 约束优化函数约束优化函数盖板优化实例：盖板优化实例： 211205 . 02602)(xxhtxf0)(0)(01321115 . 1)(014571)(014571)(01411)(2615221max4231max32