机械优化设计答案.doc

长江大学 机械工程学院 机械优化设计大作业 班级 姓名 序号 2013年5月 19-24 1、 2、 3、 4、 5、求函数的极小点。 6、某工厂有一张边长为的正方形的铁板，欲制成一个方形无盖水槽，问在该铁板的4个角处剪去多大的相等的正方形才能使水槽的容积最大？ 7、某车间生产甲（如轴）、乙（如齿轮）两种产品。生产甲种产品每件需要用材料9㎏，3个工时、4kw电，可获利60元；生产乙种产品每件需要用材料4㎏、10个工时， 5kw电,可获利120元。若每天能供应材料360㎏，有300个工时，能供电200kw电，问每天生产甲、乙两种产品各多少件，才能够获得最大的利润。 8、已知：轴上作用均布载荷 q=100N/ cm，扭矩 M=100N·m；轴长不得小于8cm；材料的许用弯曲应力 [σw]=120MPa，许用扭剪应力 [τ]= 80MPa，许用挠度 [f] = 0.01cm；密度[ρ] = 7.8t /m，弹性模量E=2×105MPa。要求：设计销轴，在满足上述条件的同时，轴的质量应为最轻。 1、 解：1) 将目标函数写成二次函数的形式f(x)=xTHx+CTx,其中： x=(x1;x2), H=(2,-4;-4,4), C=(1;-1) 线性不等式约束函数的向量矩阵和常数向量为 A=(2,1;-1,2), b=(2;2) 2)编制求解二次规划的M文件 H=[2,-4;-4,4]; C=[1;-1]; A=[2,1;-1,2]; b=[2;2]; Aeq=[]; beq=[]; lb=[0;0]; [x,f,exitflag]=quadprog(H,C,A,b,Aeq,beq,lb) M文件运行结果 Optimization terminated. x = 0 0.2500 f = -0.1250 exitflag = 1 说明 QUADPROG converged with a solution X. 则当x1=0,x2=0.25时，f有最小值-0.125. 2. 解：将 maxf 改写为 minf=-x1-2x2-x3 ,将不等式约束中的“大于等于”改写为“小于等于”， 即 2x1-x2+5x3≤6 ,则有： f=(-1,-2,-1); A=(2,1,-1;2,-1,5;4,1,1); b=（2；6；6） matlab程序为： f=[-1,-2,-2]; A=[2,1,-1;2,-1,5;4,1,1]; b=[2;6;6] lb=[0;0;0] [x,fval,exitflag]=linprog(f,A,b,[],[],lb) 运行结果为： Optimization terminated. x = 0.0000 4.0000 2.0000 fval = -12.0000 exitflag = 1 故当x1=0,x2=4,x3=2时，f有最大值12. 3、 解：1）先建立M—文件fun0.m定义目标函数，即 function f=fun0(x); f=(x(1)-x(4))^2+(x(2)-x(5))^2+(x(3)-x(6))^2 2）在建立M文件ljw0.m定义非线性约束，即 function[c,ceq]=ljw0(x) c=[x(1)^2+x(2)^2+x(3)^2-5;(x(4)-3)^2+x(5)^2-1]; ceq=[]; 3)主程序： x0=[2;0;1;4;0;1]; A=[0,0,0,0,0,1;0,0,0,-1,0,0]; b=[8;-4] [x,fval,exitflag,output]=fmincon(fun0,x0,A,b,[],[],ljw0) 运行结果为： x = 108.0495 106.0875 119.0000 108.0496 106.0874 8.0000 fval = 1.2321e+004 exitflag = 5 output = iterations: 4 funcCount: 28 lssteplength: 1 stepsize: 2.8654e-004 algorithm: medium-scale: SQP, Quasi-Newton, line-search firstorderopt: 1.7737e-004 constrviolati