现代设计方法技术总结.ppt

MATLAB的实现过程: 1、编一个目标函数的M文件jxyh5_6fun(x).m,返回x处的函数值f。 function y=jxyh5_6fun(x) y=0.785398*(4.75*x(1)*x(2)^2*x(3)^2+85*x(1)*x(2)*x(3)^285*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文件 jxyh5_6con(x).m: function [c,ceq]=jxyh5_6con(x) c(1)=17-x(2); c(2)=0.9-x(1)/(x(2)*x(3)); c(3)=x(1)/(x(2)*x(3))-1.4; c(4)=2-x(3); c(5)=x(2)*x(3)-300; c(6)=100-x(5); c(7)=x(5)-150; c(8)=130-x(6); c(9)=x(6)-200; c(10)=x(1)+0.5*x(6)-x(4)+40; c(11)=1486250/(x(2)*x(3)*x(1)^0.5)-550; c(12)=9064860*2.76*2.176/(x(1)*x(2)*x(3)^2)-400; c(13)=9064860*2.76*2.176/(x(1)*x(2)*x(3)^2)*1.56*1.794/(2.76*2.176)-400; c(14)=1.11704*x(4)^4/(x(2)*x(3)*x(5)^4)-0.003*x(4); c(15)=1/x(5)^3*((2.85*10^6*x(4)/(x(2)*x(3)))^2+2.4*10^13)^0.5-5.5; c(16)=1/x(6)^3*((2.85*10^6*x(4)/(x(2)*x(3)))^2+6*10^13)^0.5-5.5; ceq=[]; 3、给定初值,给定变量的下限约束,并调用优化过程。 x0=[230,21,8,420,120,160]; options=optimset(Display,iter,LargeScale,off); [x,fval,exitflag,output,hess]=fmincon(@jxyh5_6fun,x0,[],[],[],[],[],[],@jxyh5_6con,options) 计算结果为: x = 180.8483 25.0000 8.0377 285.8483 104.9563 130.0000 fval = 5.6271e+007 exitflag = 1 output = iterations: 5 funcCount: 35 lssteplength: 1 stepsize: 2.7552e-009 algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorderopt: 0.0016 constrviolation: 8.8854e-009 message: [1x817 char] hess = lower: [6x1 double] upper: [6x1 double] eqlin: [0x1 double] eqnonlin: [0x1 double] ineqlin: [0x1 double] ineqnonlin: [17x1 double] 例5-5 求两级圆柱齿轮减速器优化设计． 已知某两级斜齿轮圆柱齿轮减速器（图8-6）的参数：高速级输入功率P1=6.2kw,转速 n1=1450 r/min．总传动比i=31.5，齿轮宽度系数为 齿轮材料和热处理：大齿轮45号钢 正火187～207．小齿轮45号钢调质228 ～ 255HBS，工作寿命年以上．要求按照总中心矩a∑最小来 确定传动方案． 解 1) 建立优化设计的数学模型． 将影响齿轮传动总中心矩a∑的六个独立参数作为设计变量： x=(mn1, mn2, z1, z3, i, )=(x1, x2, x3, x4, x5, x6)． Mn1,mn2分别为高速级和低速级齿轮副的模数； z1, z3　　分别为高速级和低速级小齿轮齿数；