gjian.m

过程检测习题 过程检测习题

qv=22000; p=0.5*(10^6); t=100; d20=0.35; yD=17*10^(-6);yd=17*10^(-6); p1=3.48; u1=4.98e-5; k=1.4; pai=3.1415;
qvcom=0.8*qv
qvmin=1/3*qv
d=d20*(1+yd*(t-20))
disp('不锈钢管内壁粗糙度K比较少，能满足粗糙度的要求');
remin=4*qvmin/3600/(pai*d*u1)
recom=4*qvcom/3600/(pai*d*u1)
disp('确定差压计上限');
pmax=(4*qv/3600*sqrt(1-0.5^4)/(pai*0.6*0.5^2*d^2))^2/(2*p1)
if(pmax>6000)
    pcom=0.64*10e004
else
    pcom=0.64*6000
end
a2=4*qvcom/3600/(pai*d^2*sqrt(2*p1*pcom)) 
if(remin>200000)
    b0=(1+(0.6/a2)^2)^(-0.25)
else
    b0=(1+(0.6/a2+0.06)^2)^(-0.25)
end
 b(1,1)=b0;
for i=1:10
    e(1,i)=1-(0.41+0.35*b(1,i)^4)*pcom/(k*p1);
    c(1,i)=0.5959+0.0321*b(1,i)^2.1-0.1840*b(1,i)^8+0.0029*b(1,i)^2.5*(10^6/recom)^0.75;
    b(1,i+1)=(1+((e(1,i)*c(1,i)/a2)^2))^(-0.25);
    if (abs(b(1,i+1)-b(1,i))<0.0001)
       bn=b(1,i+1)
       cn=c(1,i)
       en=e(1,i)
        break
    end
end
d1=d*bn;
disp(['节流孔的直径=' num2str(d1)]);
disp('验算流量:');
q1vcom=cn*bn*4*pai*d1^2*sqrt(2*p1*pcom)/sqrt(1-bn^4)
if (q1vcom-qvcom)/qvcom<0.2e-2
    disp('合格');
else
    disp('不合格,请重新计算');
end
    d20=d/(1+yd*(t-20));
    opd20=0.0005*d20;
    disp(['opd20=+_' num2str(opd20)]);
    dp=(sqrt(1-bn^4)-cn*bn^2)/(sqrt(1-bn^4)-cn*bn^2)*1e4;
    disp(['压力损失=' num2str(dp)]);