- 1、本文档共4页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
C语言编程换热器模拟编程.doc
# include stdio.h
# include math.h
main()
{double Do=0.025,Di=0.021,L=6,Wh=29.5; Th1=280,Wc=37.5, Tc1=160,Ro=0.0005,Ri=0.0001,D1,D2,D2O1=0.85,D2O2=0.919,Tc2,Th2,
MD2O1,MD2O2,K=12.5,Cp1,Cp2,ramda1,ramda2,niu1,niu2,
yita1,yita2,a1,b1,a2,b2,Tc20,Th20,Tmc,Tmh,C,Qc,Qh,
Ho,Hi,Reo,Rei,Pro,Pri,Si,rou1,ui,Hi0,Twi0,
Ai,Twi,niuwi,yitawi,Dwi,Xwi,rouwi,Ao,
Smax,B=0.23,D=0.8,t=0.032,De,uo,rou2,
Ho0,Two0,Two,niuwo,yitawo,Dwo,Xwo,rouwo,
Ko,NTU,E,F,Cmin,Cmax,Xh,Xc,ld=44,b=0.002,Dm=0.0023,e=2.71828;
int n=0;
Tc2=180;
Th2=260;
do{ n++;
printf(n=%d\n,n);
Th20=Th2;
Tmh=(Th1+Th20)/2;
do{Tc20=Tc2;
Tmc=(Tc1+Tc20)/2;
Cp1=(0.7072+(0.00147-0.00051*D2O1)*Tmc-0.318*D2O1)*(0.055*K+0.35)*4.18*1000;
Cp2=(0.7072+(0.00147-0.00051*D2O2)*Tmh-0.318*D2O2)*(0.055*K+0.35)*4.18*1000;
C=Wh*Cp2/(Wc*Cp1);
Tc2=Tc1+C*(Th1-Th20);
}while(fabs(Tc2-Tc20)=0.1);
Qh=Wh*Cp2*(Th1-Th2);
Qc=Wc*Cp1*(Tc2-Tc1);
Ai=3.14*Di*L*324;
Si=0.25*3.14*Di*Di*324/2;
Xc=1+Tmc/100.0;
MD2O1=pow(D2O1,2);
D1=0.942+0.248*Xc+0.174*MD2O1+0.0841/(Xc*D2O1)-0.312*Xc/D2O1-0.556*exp(-Xc);
rou1=1000*D1;
ui=Wc/(rou1*Si);
b1=log((log(90.0+1.22)/log(13+1.22)))/(log((50.0+273)/(100.0+273)));
a1=log(log(90+1.22))-b1*log(50.0+273);
niu1=exp(exp(a1+b1*log(Tmc+273)))-1.22;
yita1=niu1*rou1/1000000;
Rei=Di*ui*rou1/yita1;
ramda1=0.4213*(1-0.00054*Tmc)/D2O1/3.6;
Pri=Cp1*yita1/ramda1;
Hi0=0.027*pow(Rei,0.8)*pow(Pri,0.33)*ramda1/Di;
Twi0=Tmc+Qc/(Hi0*Ai);
do
{ niuwi=pow(e,pow(e,a1+b1*log(Twi0+273)))-1.22;
Xwi=1+Twi0/100.0;
Dwi=0.942+0.248*Xwi+0.174*MD2O1+0.0841/(Xwi*D2O1)-0.312*Xwi/D2O1-0.556*exp(-Xwi);
rouwi=1000*Dwi;
yitawi=niuwi*rouwi/1000000;
Hi=Hi0*(pow((yita1/yitawi),0.14));
Twi=Tmc+Qc/(Hi*Ai);
Twi0=Twi;
}while(fabs(Twi-Twi0)=0.5);
Ao=3.14*Do*L*324;
Smax=B*D*(1-Do/t);
De=4*(t*t-0.25*3.14*Do*Do)/(3.14*Do);
Xh=1+Tmh/100.0;
MD2O2=pow(D2O2,2);
D2=0.
文档评论(0)