管壳式程序c语言,换热器模型原程序（C语言）.doc

换热器模型原程序(C语言)

#include

#include

float density(float t,float d)

{float x,z;

x=1+(t/100);

z=0.942+0.248*x+0.174*d*d+0.0841/(x*d)-0.312*x/d-0.556*exp(-x);

z=z*1000;

return(z);

}

/*以下为比热容的计算：*/

float hcap(float t,float d)

{float cp;

cp=(0.7072+(0.00147-0.000551*d)*t-0.318*d)*(0.055*12.5+0.35);

cp=cp*4186;

return(cp);

}

/*以下为热导率的计算：*/

float hcon(float t,float d)

{float cn;

cn=0.1008*(1-0.00054*t)/d;

cn=cn*4186/3600;

return(cn);

}

float nianc(float t,float d)

{float ni;

ni=exp(exp(22.81142-3.68738*log(t+273)))-1.22;

ni=ni*density(t,d)*0.000001;

return(ni);

}

float nianh(float t,float d)

{float ni;

ni=exp(exp(18.9173-2.92782*log(t+273)))-1.22;

ni=ni*density(t,d)*0.000001;

return(ni);

}

main()

{ int n,N;

float tc1,tc2,th1,th2,k,q,e,h,hi,ho,ntu,rei,reo,pri,pro,wc,tao,

wh,tmc,tmh,ai,ao,twi,two,cpc,cph,rdc,rdh,ndc,ndh,di,dou,dk,

d1,d2,s1,s2,cmin,cmax,de,p,l,v1,v2,rs,rt,thx,tx,ty,thy,tcx,tcy;

tc1=160; th1=280; wh=29.5; N=324; wc=37.5; di=0.021; dou=0.025; dk=0.8;

l=6.0; p=0.032; n=0; rs=0.0005; rt=0.0001;

thx=th1;

do{

tmh=(th1+thx)/2;

cph=hcap(tmh,0.919);

q=wh*cph*(th1-thx);

tcx=tc1;

do{

tmc=(tc1+tcx)/2;

tc2=tc1+q/(wc*hcap(tmc,0.850));

tcy=tcx;

tcx=tc2;

}

while(fabs(tc2-tcy)>0.1);

tmc=(tc1+tc2)/2;

d1=density(tmc,0.850);

d2=density(tmh,0.919);

cpc=hcap(tmc,0.850);

rdc=hcon(tmc,0.850);

rdh=hcon(tmh,0.919);

ndc=nianc(tmc,0.850);

ndh=nianh(tmh,0.919);

s1=3.1415926*0.25*di*di*N*0.5;

v1=wc/d1/s1;

ai=3.1415926*di*N*l;

rei=di*v1*d1/ndc;

pri=cpc*ndc/rdc;

s2=0.230*dk*(1-dou/p);

v2=wh/d2/s2;

de=4*(p*p-3.1415926*0.25*dou*dou)/(3.1415926*dou);

reo=de*v2*d2/ndh;

pro=cph*ndh/rdh;

ao=3.1415926*dou*N*l;

h=0.023*pow(rei,0.8)*pow(pri,0.33)*rdc/di;

hi=h;

twi=tmc+q/(h*ai);

do{h