punish1.m

使用惩罚函数法对机械设计中的问题进行优化处理

function y1=punish1(x0,x,r)
syms d D n;
d1=x0(1);
D1=x0(2);
n1=x0(3);
if(d1-D1+16<=0)
    g1=0;
else
    g1=d-D+16;
end
if(2.5-d1<=0)
    g2=0;
else
    g2=2.5-d;
end
if(2-n1<=0)
    g3=0;
else
    g3=2-n;
end
if((n1+1)*d1-50<=0)
    g4=0;
else
    g4=(n+1)*d-50;
end
C1=D1/d1;
K1=0.615/C1+(4*C1-1)/(4*C1-4);
C=D/d;
K=0.615/C+(4*C-1)/(4*C-4);
if(K1*8*700*D1/(pi*d1^3)-444<=0)
    g5=0;
else
    g5=K*8*700*D/(pi*d^3)-444;
end
if(10+n1*0.15*d1+(n1+1)*d1-2.6*D1<=0)
    g6=0;
else
    g6=(10+n*0.15*d+(n+1)*d)/D-2.6;
end
if(8*700*D1^3*n1-10*8.1*10^4*d1^4<=0)
    g7=0;
else
    g7=8*700*D^3*n/(8.1*10^4*d^4)-10;
end
if(d1-6<=0)
    g8=0;
else
    g8=d-6;
end
if(10-D1<=0)
    g9=0;
else
    g9=10-D;
end
if(D1-60<=0)
    g10=0;
else
    g10=D-60
end
m=0.0078*pi*d^2/4*(n+1.5)*pi*D/cos(pi/36);
y=m+r*(g1^2+g2^2+g3^2+g4^2+g5^2+g6^2+g7^2+g8^2+g9^2+g10^2);
y1=[diff(y,d);diff(y,D);diff(y,n)];
d=x(1);
D=x(2);
n=x(3);
y1=subs(y1);