热管式换热器设计计算程序源代码
#include
#include
float insertcal(float t1,float t2,float l1,float l2,float x){
return l1+(l2-l1)*(x-t1)/(t2-t1);
}
main(){
float t2p[33]={-50,-40,-30,-20,-10,0,10,20,30,40,50,60,70,80,90,
100,120,140,160,180,200,250,300,350,400,500,600,700,
800,900,1000,1100,1200},t1p[13]={0,100,200,300,400,500,600,
700,800,900,1000,1100,1200};
float cp2p[33]={1.013,1.013,1.013,1.009,1.009,1.005,
1.005,1.005,1.005,1.005,1.005,1.005,1.009,1.009,1.009,
1.009,1.009,1.013,1.017,1.022,1.026,1.038,1.047,1.059,
1.068,1.093,1.114,1.135,1.156,1.172,1.185,1.197,1.210},
cp1p[13]={1.042,1.068,1.097,1.122,1.151,1.185,1.214,1.239,
1.264,1.290,1.306,1.323,1.340};
float miu1p[13]={12.2,21.54,32.8,45.81,60.83,76.3,93.61,112.1,
131.8,152.5,174.3,197.1,221.0},miu2p[33]={9.23,10.04,10.80,
11.61,12.43,13.28,14.16,15.06,16.00,16.96,17.95,18.97,20.02,
21.09,22.10,23.13,25.45,27.8,30.09,32.49,34.85,40.61,48.33,
55.46,63.09,79.38,96.89,115.4,134.8,155.1,177.1,199.3,233.7};
float limuda1p[13]={2.28,3.13,4.01,4.84,5.7,6.56,7.42,8.27,9.15,10.0,
10.9,11.75,12.62},limuda2p[33]={2.04,2.12,2.20,2.28,2.36,2.44,2.51,
2.59,2.67,2.76,2.83,2.9,2.96,3.05,3.13,3.21,3.34,3.49,3.64,3.78,
3.93,4.27,4.6,4.91,5.21,5.74,6.22,6.71,7.18,7.63,8.07,8.5,9.15};
float Pr1p[13]={0.72,0.69,0.67,0.65,0.64,0.63,0.61,0.60,0.59,0.58,0.57,0.56},
Pr2p[33]={0.728,0.728,0.723,0.716,0.712,0.707,0.705,0.703,0.701,0.699, 0.698,0.696,0.694,0.692,0.690,0.688,0.686,0.684,0.682,0.681,0.680,0.677,
0.674,0.676,0.678,0.687,0.699,0.706,0.713,0.717,0.719,0.722,0.724};
float density1p[13]={1.295,0.950,0.748,0.617,0.525,0.457,0.405,0.363,0.330,0.301,
0.275,0.257,0.240},density2p[33]={1.584,1.515,1.453,1.395,1.342,1.293,
1.247,1.205,1.165,1.128,1.093,1.000,0.972,0.946,0.898,0.854,0.815,0.779,
0.746,0.674,0.615,0.566,0.524,0.456,0.404,0.362,0.329,0.301,0.277,0.257,
0.239};
int i,j,flag;
float Q1,V20,t2in,Q,density20,cp2,Q2,t2ou,step,dett2,t2;
float t1in,t1ou,t1,density10=1.295,V10,vv1,vv2,cp1,dett1;
float l1=1.25,l2=l1,u10=5,B,n,ST=0.065,zeita,kesai1,kesai2,u1,u2,u20,Re1,Re2,
d0=0.025,miu1,miu2,h1,h2,limuda1,limuda2;
float Sf=0.006,detf=0.001,S=Sf+detf,lf=0.0125,nf=floor(1/S),Pr1,Pr2;
float beita,df=d0+2*lf,pi=3.1415926,limudaw=45.3,m1,m2,yeta1,yeta2,e=2.718281828,R;
float di=0.020,detTm,q,m,N,NFV,SL=0.090,D1,D2,detP1,detP2,density1,density2,Nreal;
printf("请输入标准状态下烟气流量(m3/h)\n");
scanf("%f",&vv1);V10=vv1/3600;
printf("烟气入口温度('c)\n");
scanf("%f",&t1in);
printf("烟气出口温度('c)\n");
scanf("%f",&t1ou);dett1=t1in-t1ou;t1=(t1ou+t1in)/2;
printf("请输入标准状态下空气流量(m3/h)\n");
scanf("%f",&vv2);V20=vv2/3600;
printf("空气入口温度('c)\n");
scanf("%f",&t2in);
for(flag=0,i=0;i<33;i++){
if(t1p[i]<=t1&&t1<=t1p[i+1])
{flag=i;goto loop;}}
loop:j=flag+1;
cp1=insertcal(t1p[flag],t1p[j],cp1p[flag],cp1p[j],t1);
Q1=density10*V10*c
p1*dett1;
Q=0.95*Q1;
t2ou=80;
do{
for(flag=0,i=0;i<33;i++){
if(t2p[i]<=t2ou&&t2ou<=t2p[i+1])
{flag=i;goto loop1;}}
loop1:j=flag+1;
density20=1.293;
cp2=insertcal(t2p[flag],t2p[j],cp2p[flag],cp2p[j],t2ou);
dett2=t2ou-t2in;
Q2=density20*cp2*V20*dett2;
if(Q
}
else{
step=+0.1;}
t2ou=t2ou+step;
}while(abs(Q-Q2)>0.5);
printf("\n计算得空气出口温度\t%.2f('c)\n",t2ou);
/*printf("\nV10=%f m3/s\n",V10);
printf("\nV20=%f m3/s\n",V20);*/
t2=(t2in+t2ou)/2;
B=V10/(l1*u10);n=floor(B/ST);
zeita=(ST-d0-2*lf*detf*nf)/ST;
u1=u10/zeita;
u20=V20/(B*l2);u2=u20/zeita;
for(flag=0,i=0;i<33;i++){
if(t1p[i]<=t1&&t1<=t1p[i+1])
{flag=i;goto loop2;}}
loop2:j=flag+1;
miu1=insertcal(t1p[flag],t1p[j],miu1p[flag],miu1p[j],t1)/1000000;
limuda1=insertcal(t1p[flag],t1p[j],limuda1p[flag],limuda1p[j],t1)/100;
Pr1=insertcal(t1p[flag],t1p[j],Pr1p[flag],Pr1p[j],t1);
density1=insertcal(t1p[flag],t1p[j],density1p[flag],density1p[j],t1);
Re1=u1*d0/miu1;h1=0.137*(limuda1/d0)*pow(Re1,0.6338)*pow(Pr1,1/3);
for(flag=0,i=0;i<33;i++){
if(t2p[i]<=t2&&t2<=t2p[i+1])
{flag=i;goto loop3;}}
loop3:j=flag+1;
miu2=insertcal(t2p[flag],t2p[j],miu2p[flag],miu2p[j],t2)/1000000;
limuda2=insertcal(t2p[flag],t2p[j],limuda2p[flag],limuda2p[j],t2)/100;
Pr2=insertcal(t2p[flag],t2p[j],Pr2p[flag],Pr2p[j],t2);
density2=insertcal(t2p[flag],t2p[j],density2p[flag],density2p[j],t2);
Re2=u2*d0/miu2;h2=0.137*(limuda2/d0)*pow(Re2,0.6338)*pow(Pr2,1/3);
beita=((pi*(pow(df,2)-pow(d0,2))/2+pi*df*detf)*nf+pi*d0-pi*d0*detf*nf)/(pi*d0);
m1=sqrt(2*h1/limudaw/detf);m2=sqrt(2*h2/limudaw/detf);
yeta1=tanh(m1*lf)/(m1*lf);
yeta2=tanh(m2*lf)/(m2*lf);
R=1/(h1*pi*d0*beita*yeta1)+1/(pi*limudaw*l1)*log(d0/di)+1/(h2*pi*d0*beita*yeta2);
detTm=((t1in-t2ou)-(t1ou-t2in))/log((t1in-t2ou)/(t1ou-t2in));
q=detTm/R;
N=(Q*1000/q)*1.2;
m=ceil(N/n);
/*printf("\nRe1=%f\n",Re1);
printf("\nRe2=%f\n",Re2);
printf("\nbeita=%f\n",beita);*/
Nreal=m*n;
printf("\n计算所需热管数\t%.0f(根)\t%.0f(列)\t%.0f(排)\n",Nreal,n,m);
NFV=ST*SL-pi*d0*d0/4-pi*(df*df-d0*d0)*detf*nf/4;
D1=4*NFV/(B*d0*l1*beita);D2=4*NFV/(B*d0*l2*beita);
kesai1=1.92*pow(Re1,-0.145);
kesai2=1.92*pow(Re2,-0.145);
detP1=m*density1*u1*u1*kesai1/2;
detP2=m*density2*u2*u2*kesai2/2;
/*printf("\nu1=%f \n",u1);
printf("\nu2=%f \n",u2);
printf("\nzeita=%f \n",zeita);
printf("\nbeita=%f \n",beita);
printf("\nR=%f \n",R);
printf("\nyeta1=%f \n",yeta1);
printf("\nyeta2=%f \n",yeta2);
printf("\nq=%f \n",q);
printf("\nNFV=%f m3\n",NFV);
printf("\nD1=%f m\n",D1);
printf("\nD2=%f m\n",D2);
printf("\nkesai1=%f\n",kesai1);
printf("\nkesai2=%f\n",kesai2);*/
printf("\n烟气侧压力损失\t%.2f(Pa)\n\n空气侧压力损失\t%.2f(Pa)\n",detP1,detP2);
getch();
}