热管式换热器设计计算程序源代码

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);
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);
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},
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,
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;
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);
loop:j=flag+1;
cp1=insertcal(t1p[flag],t1p[j],cp1p[flag],cp1p[j],t1);
Q1=density10*V10*cp1*dett1;
Q=0.95*Q1;
t2ou=80;
do{
for(flag=0,i=0;i<33;i++){
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;
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);
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};
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;
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,
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));
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);
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;
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++){
dett2=t2ou-t2in;
Q2=density20*cp2*V20*dett2;
if(Q<Q2){
step=-0.1;
}
else{
step=+0.1;}
t2ou=t2ou+step;
}while(abs(Q-Q2)>0得空气出口温度\t%.2f('c)\n",t2ou);
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,
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,
#include <stdio.h>
#include <math.h>
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,
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);
printf("烟气入口温度('c)\n");
scanf("%f",&t1in);
printf("烟气出口温度('c)\n");
scanf("%f",&t1ou);dett1=t1in-t1ou;t1=(t1ou+t1in)/2;