油库上机程序2

#include
#include
main()
{
double D=15.25,H=11.74,h,tu=0,tqi=-3,pusai,tj,G,c,F,rou,tao,td,zh18=0;
double tz,ty,tbi,v,dao,B,gr,pr,x,yi,n,a1bi,a2bi,a3bi,Re,K,Kbi,GP=10,Q,Q1,Q2;
pusai=H/D;h=0.85*H;
tj=(tu+(pusai*4+1)*tqi)/(2+4*pusai);
rou=0.89-0.00066*(40-15);
G=rou*3.14*D*D*H*0.85*1000/4;
F=3.14*D*H+3.14*D*D/4+3.14*D*D*sqrt(2)/4;
tao=16*24;
td=(tqi+40)/2;
while(fabs(td-zh18)>0.001)
{
zh18=td;
c=(4.1868*(0.403+0.81*td/1000))*1000/sqrt(0.89)/3600;
tz=tj+(40-tj)*exp(-1.454*F*tao/G/c);
td=(tz+40)/2;
}
printf("储存16天后温度%.2f摄氏度\n\n",tz);
if((40-tj)/(tz-tj)<=2)
ty=(40+tz)/2;
else
ty=tj+(40-tz)/log((40-tj)/(tz-tj));
tbi=(ty+tj)/2;
while(GP>1)
{tbi=tbi+0.01;
td=(tbi+ty)/2;
rou=0.89-0.00066*(td-15);
v=8.921-3.439*log10(273.15+td);v=pow(10,v);v=pow(10,v)-0.6;v=v/1000000;
dao=117.5/890*(1-0.00054*td);
c=(4.1868*(0.403+0.81*td/1000))/sqrt(0.89);
B=((rou-0.88)*(0.718-0.739)/(0.89-0.88)+0.739)/1000;
gr=(pow(h,3)*B*(ty-tbi)*9.81)/pow(v,2);
pr=(v*c*rou*1000000)/dao;
x=pr*gr;
if(x>0.001&&x<500)
{yi=1.18;n=1.0/8;}
else if(x>500&&x<20000000)
{yi=0.54;n=1.0/4;}
else
{yi=0.135;n=1.0/3;}
a1bi=(yi*dao*pow(x,n))/h;
Re=3*D*1000000/13.03;
if(Re>5&&Re<30)
{yi=0.81;n=0.4;}
else if(Re>80&&Re<5000)
{yi=0.625;n=0.46;}
else if(Re>5000&&Re<500000)
{yi=0.197;n=0.6;}
else if(Re>500000)
{yi=0.023;n=0.8;}
a2bi=yi*2.418*pow(Re,n)/100/D;
a3bi=0.96*5.67*(pow((tbi+273)/100,4)-pow((tqi+273)/100,4))/(tbi-tqi);
Kbi=1/(1/a1bi+1/(a2bi+a3bi));
GP=fabs(tbi+Kbi*(ty-tqi)/a1bi-ty);}
printf("取罐壁温度为%.2f摄氏度时罐壁导热系数为%.2fW/m2.C\n",tbi,Kbi);
K=(Kbi*3.14*D*h+(3.14*D*(H-h)+3.14*D*D*sqrt(2)/4)*0.7+0.12*3.14*D*D/4)/F;
printf("取罐顶罐底传热系数为经验值求得总传热系数K为%.2fW/m2.C\n",K);
c=(4.1868*(0.403+0.81*ty/1000))/sqrt(0.89);
Q1=G*c*(40-tz);Q2=K*F*(ty-tj);
Q=Q1*1000/48/3600+Q2;
printf("单位时间所需总热量为%.2fW 加热器所用蒸汽量为%.2fkg/s\n",Q,Q/1000/2114.75);
tbi=110;GP=10;
while(GP>0.5)
{tbi=tbi+0.01;
td=(ty+tbi)/2;
v=8.921-3.439*log10(273.15+td);v=pow(10,v);v=pow(10,v)-0.6;v=v/1000000;
rou=0.89-0.00066*(td-15);
c=(4.1868*(0.403+0.81*td/1000))/sqrt(0.89);
dao=117.5/890*(1-0.00054*td);
B=((rou-0.88)*(0.718-0.739)/(0.89-0.88)+0.739)/1000;
gr=(pow(0.06,3)*B*(tbi-ty)*9.81)/pow(v,2);
pr=(v*c*rou*1000000)/dao;
x=pr*gr;
if(x>0.001&&x<500)
{yi=1.18;n=1.0/8;}
else if(x>500&&x<20000000)
{yi=0.54;n=1.0/4;}
else
{yi=0.135;n=1.0/3;}
a2bi=(yi*dao*pow(x,n))/0.06;
K=1/(1/a2bi+0.00172);
GP=fabs(ty+K*(151.11-ty)/a2bi-tbi);}
printf("加热器面积%.2fm2\n\n",Q/(151.11-ty)/K);
printf("如采用维温加热\n\n");
ty=40;tbi=25;GP=10;
while(GP>0.5)
{tbi=tbi+0.01;
td=(tbi+ty)/2;
rou=0.89-0.00066*(td-15);
v=8.921-3.439*log10(273.15+td);v=pow(10,v);v=pow(10,v)-0.6;v=v/1000000;
dao=117.5/890*(1-0.00054*td);
c=(4.1868*(0.403+

0.81*td/1000))/sqrt(0.89);
B=((rou-0.87)*(0.739-0.76)/(0.88-0.87)+0.76)/1000;
gr=(pow(h,3)*B*(ty-tbi)*9.81)/pow(v,2);
pr=(v*c*rou*1000000)/dao;
x=pr*gr;
if(x>0.001&&x<500)
{yi=1.18;n=1.0/8;}
else if(x>500&&x<20000000)
{yi=0.54;n=1.0/4;}
else
{yi=0.135;n=1.0/3;}
a1bi=(yi*dao*pow(x,n))/h;
Re=3*D*1000000/13.03;
if(Re>5&&Re<30)
{yi=0.81;n=0.4;}
else if(Re>80&&Re<5000)
{yi=0.625;n=0.46;}
else if(Re>5000&&Re<500000)
{yi=0.197;n=0.6;}
else if(Re>500000)
{yi=0.023;n=0.8;}
a2bi=yi*2.418*pow(Re,n)/100/D;
a3bi=0.96*5.67*(pow((tbi+273)/100,4)-pow((tqi+273)/100,4))/(tbi-tqi);
Kbi=1/(1/a1bi+1/(a2bi+a3bi));
GP=fabs(tbi+Kbi*(ty-tqi)/a1bi-ty);}
printf("取罐壁温度为%.2f摄氏度时罐壁导热系数为%.2fW/m2.C\n",tbi,Kbi);
K=(Kbi*3.14*D*h+(3.14*D*(H-h)+3.14*D*D*sqrt(2)/4)*0.7+0.12*3.14*D*D/4)/F;
printf("取罐顶罐底传热系数为经验值求得总传热系数K为%.2fW/m2.C\n",K);
c=(4.1868*(0.403+0.81*ty/1000))/sqrt(0.89);
Q=K*F*(ty-tj);
printf("单位时间所需总热量为%.2fW 加热器所用蒸汽量为%.2fkg/s\n",Q,Q/1000/2114.75);
tbi=110;GP=10;
while(GP>0.5)
{tbi=tbi+0.01;
td=(ty+tbi)/2;
v=8.921-3.439*log10(273.15+td);v=pow(10,v);v=pow(10,v)-0.6;v=v/1000000;
rou=0.89-0.00066*(td-15);
c=(4.1868*(0.403+0.81*td/1000))/sqrt(0.89);
dao=117.5/890*(1-0.00054*td);
B=((rou-0.87)*(0.739-0.76)/(0.88-0.87)+0.76)/1000;
gr=(pow(0.06,3)*B*(tbi-ty)*9.81)/pow(v,2);
pr=(v*c*rou*1000000)/dao;
x=pr*gr;
if(x>0.001&&x<500)
{yi=1.18;n=1.0/8;}
else if(x>500&&x<20000000)
{yi=0.54;n=1.0/4;}
else
{yi=0.135;n=1.0/3;}
a2bi=(yi*dao*pow(x,n))/0.06;
K=1/(1/a2bi+0.00172);
GP=fabs(ty+K*(151.11-ty)/a2bi-tbi);}
printf("加热器面积%.2fm2\n",Q/(151.11-ty)/K);
}