MATLAB优化设计

1、MATLAB优化设计学院:专业:班级:学号:姓名:指导老师：机电学院机械设计制造及其自动化072&&&-*20131*大禹陈祯2015年10月25日题目11求解如下最优化问题minf(x)=-2x1-6x2x2-2x1x22x2subjectto-x12x2工2x1x2_2x1，x2-0步骤一：对已有的数学模型matlab编程1 .编写.m文件并保存：%实对称矩阵%列向量%对应维数矩阵%列向量h=2,-2;-2,4;f=-2;-6;a=1,1;-1,2;b=2;2;lb=zeros(2,1);x,value=quadprog(h,f,a,b,lb)2 .运行.m文件结果

2、如图1.0所示：Voluc：IMinIVIrk1.1.-1刁-124】22-6-22P-2,-2,4-24O；OJO-7.2000-Z.2OOQ-7.3300CiSOOD:1.200010.8000-1.2OOCeeHunabfhlbva蕈xm.isea-tian.七ohhizLgk：七已d-£=O-30001.2000wwlIue=-7.2000>>图1.0题目一文件运行结果步骤二：matlab运行结果分析阶段由图1.0知，当x1=0.8,x2=1.2时，minf(x)=-7.2。题目22、某农场拟修建一批半球壳顶的圆筒形谷仓，计划每座谷仓容积为300立方米，圆筒半径不

3、得超过3米，高度不得超过10米。半球壳顶的建筑造价为每平方米150元，圆筒仓壁的造价为每平方米120元，地坪造价为每平方米50元，求造价最小的谷仓尺寸为多少？步骤一：题目分析阶段设：圆筒的半径为R,圆筒的高度为Ho谷仓的容积为300立方米，可得：223二R2H-二R3=3003圆筒高度不得超过10米，可得：0-H-10圆筒半径不得超过3米，可得：0-R-3当造价最小时：minf(R,H)=1502二R21202二R2H50二R2步骤二：数学模型建立阶段minf(R,H)=1502二R21202二R2H50二R2_22_3二RH二R=300s.t.30,H_100_R_3步骤三：matlab编程

4、阶段1 .编写myfun3.m文件并保存：functionf=myfun2(x)f=350*3.14*x(1)A2+240*3.14*x(1)*x(2);%目标函数2 .编写mycon3.m文件并保存：functiong,ceq=mycon2(x)q=;ceq=2*3.14*x(1)A3/3+3.14*x(1)A2*x(2)-300;%约束等式3 .编写文件final.m保存：X=fmincon(myfun2,3;3,0;0,3;10,mycon2)4 .运行final.m结果如图2.0所示：闻.a之qValueMiihMg2.9369s+042,93693+042,9369e+G4田k3田.

5、615733,6157ActiveinequalitieslowerupperJ(towithinoptions,TolCon=le006):meqlinineqnonlin3.00008.6157fval=2.9369e+004图2.0题目二文件运行结果步骤四:matlab运行结果分析阶段由图2.0知，当圆筒的半径R为3m,高度H为8.6157m时,满足体积要求的谷仓最小造价为29369元。题目33、已知轴一端作用载荷P=1000N/cm,扭矩M=100Nm,轴长不小于8cm,材料的许用弯曲应力为120MPa,许用扭剪应力为80MPa,许用挠度为0.01cm,密度为7.8t/m,弹性模量为2

6、00000MPa,设计该轴，使得满足上述条件，且重量最轻。P图3.0示意图步骤一：受力分析阶段。根据题意轴端受力F=pd轴受最大的弯矩M弯=Fs=pdl=100000dld4轴的横截面对中性轴的惯性矩Iz为：Iz=64几轴的最大弯曲正应力(7须满足MyM弯d32000001xIz=2Iz=d，<(r=120pa轴受最大扭切应力须满足M工亘1600宜ax=Wp=16=71d3<T=MCpa轴的最大挠度WB须满足_333Fl64pdl64134 _,3wb=3EI=3Edn=630000000d7r<f=0.0001m轴的质量m为d2一冗3M=pV=pA1=41汉.8X03步骤二

7、：数学模型建立阶段轴径记为x1米，轴长1为x2米，建立数学模型：x12冗3minf(x1,x2)=4"2>7.8M03200000x25 .t.g1(x1,x2)=.x12n-120000000<01600g2(x1,x2)=以13-8000000000-364x27zg3(x1,x2)=630000000x1冗-0.0001<0g4(x1,x2)x1三0,g5(x1,x2)-x20.08<0步骤三：matlab编程阶段1 .编写myfun.m文件并保存：functionf=myfun(x)f=x(1)A2*x(2)*3.14*7800/4;%目标函数2 .编

8、写mycon.m文件并保存：functiong,ceq=mycon(x)g(1)=3200000*x(2)/(x(1)A2*3.14)-120000000;%约束条件g(2)=1600/(x(1)A3*3.14)-80000000;%约束条件g(3)=64*x(2)A3/(x(1)A3*3.14159*630000000)-0.0001;%约束条件g(4)=-x(1);%约束条件g(5)=-x(2)+0.08;%约束条件ceq=;3 .编写文件final.m保存:x0=0.002;0.009;lb=0.001;0.008;ub=1000;1000;x,fval=fmincon(myfun,x0

9、,lb,ub,mycon);x1=x(1);x2=x(2);display(x(1);display(x(2);display(fval);4 .运行final.m结果如图3.1所示：ValueMinMax.0.3328O33280.332B1,0JD0e-O3:0.QOS01.DO.0.00001000;1000IODOWOOO2fi1;0OSDQ)O02610800CXJ20；00090)00020OCW0.026100261,02610.08000.08000.0800Actqua.1iltxes11a'wi±ltirioptions-TdLCon=1e-006)zlo

10、werupperxneQlinineqjionlin5OiThS=602ei3X15=o-osoafval=CL332E图3.1题目三文件运行结果步骤四：matlab运行结果分析阶段由图3.1知，当轴径d为26.1mm，轴长l为80mm时，轴的质量最轻为0.3328千克。总结优化设计是近年发展起来的一门新学科，它为工程领域提供了一种新的重要的科学设计方法。当然它在机械领域也是运用非常广泛。优化设计课程主要以Matlab软件作为运算处理工具，充分利用了Matlab交互式程序设计及强大的数值分析能力。学校给优化设计这门课程仅仅安排了五周的课时，我相信要深入学习这些课时是远远不够的，我们只有自己动手

11、去写程序，自己去摸索公式，才能将优化设计这种方法运用在机械设计之中。五周的时间虽然很短，但在陈祯老师的悉心指导下，我们走进这门学科的大门，所谓师傅领进门修行靠个人，我还得不断地去摸索，我相信在以后从事机械行业后会不断运用到优化设计的思想，做出最好用的设计产品。之前在上学期通过机械原理课程设计的学习我就开始接触到Matlab软件，但是程度非常有限。这学期通过优化设计的学习，再次让我感到Matlab软件的强大计算功能，同样在这个过程中我的编程能力也得到提高，也能在马上就要学的单片机课程中大显身手。在解答本次的大作业的时候，刚开始看到这三个工程实例我觉得好难，这下子该咋办，根本不知如何入手，于是我又重新去钻研优化设计课本，但是课本关于编程的知识很少，我去图书馆借了好几本关于Matlab编程的书本，苦心钻研，然后结合老师讲课的PPT和老师给的资料，发现好多东西都是相通的，马上找到方法，对症下药，很快也就得出结果了，终于满头的雾水烟消云散。短短五周的优化设计课程结束了，但我相信在以后从事机械行业后会不断运用到优化设计的思想和方法。感谢陈祯悉心指导把我们带进了“优化设计”的大门，让我们能够通