换热器程序.docVIP

化工过程综合与分析换热器程序 源程序： /*化工过程分析与综合换热器*/ #includestdio.h #includemath.h /*Crude 原油壳程冷物料 Residue 渣油管程热物料*/ #define e 2.71828 /*自然数*/ #define D20_C 0.850 /*20度时 原油（冷 壳程) 的相对密度，可查取或实测*/ #define D20_R 0.919 /*20度时 渣油（热 管程） 的相对密度，可查取或实测*/ #define K_C 12.5 /*原油特性因数，可查取或实测*/ #define K_R 12.5 /*渣油特性因数，可查取或实测*/ #define er 0.5 /*壁温迭代收敛判据/0C */ #define ec 0.01 /*换热器流体出口温度迭代收敛判据/0C */ #define pi 3.1415926 double frd(double Tm, double D20) /*密度Relative density计算函数*/ { double D=0,X=0; X=1+Tm/100; D=0.942+0.248*X+0.174*pow(D20,2)+0.0841/(X*D20)-0.312*X/D20-0.556*pow(e,-X); return D*1000; } double fcp(double Tm,double D20,double K) /*比定压热容计算函数*/ { double a=0; a=4.1855*(0.6811-0.308*(0.99417*D20+0.009181)+(1.8*Tm+32)*(0.000815-0.000306* (0.99417*D20+0.009181)))*(0.055*K+0.35); return a*1000; } double ftc(double Tm,double D20) /*热导率Thermal conductivity计算函数*/ { double a=0; a=0.4213*(1-0.00054*Tm)/D20; return a*1000.0/3600; } double fdv_C(double Tm) /*原油动力粘度Kinematic viscosity计算函数*/ { double v=0,a=0,b=0, t1=50,t2=100,v1=90,v2=13; b=log(log(v1+1.22)/log(v2+1.22))/log((t1+273)/(t2+273)); a=log(log(v1+1.22))-b*log(t1+273); v=pow(e,pow(e,a+b*log(Tm+273)))-1.22; return v*frd(Tm,D20_C)*pow(10,-6); /* 动力黏度等于运动粘度乘以密度*/ } double fdv_R(double Tm) /*渣油动力粘度Kinematic viscosity计算函数*/ { double v=0,a=0,b=0, t1=50,t2=100,v1=1500,v2=120; b=log(log(v1+1.22)/log(v2+1.22))/log((t1+273)/(t2+273)); a=log(log(v1+1.22))-b*log(t1+273); v=pow(e,pow(e,a+b*log(Tm+273)))-1.22; return v*frd(Tm, D20_R)*pow(10,-6); /* 动力黏度等于运动粘度乘以密度*/ } void main() { int i=0,j=0,m=0,l=0,NB=0,no=0; double Lg=0,L=0,N=0,Np=0,P=0,Di=0,Do=0,Dm=0,De=0,Ds=0,A=0,AT=0,AS=0, Ai=0,Ao=0, Wc=0,Wh=0,rdc=0,rdh=0, cpc=0,cph=0,tcc=0,tch=0,kvc=0, kvh=0, dvc=0,dvh=0,dvwi=0,dvwo=0, Tc1=0,T