原创力文档- 1、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
- 2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
- 4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
- 5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传，每下载1次， 网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
管壳式换热器数学模型模拟与计算
第 PAGE 1 页
换热器模型原程序（C语言）
#includemath.h
#includestdio.h
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);
文档评论（0）