传热问题设计与编程.doc
#include<stdio.h>#include<math.h>main()float ms2=7500.00/3600,/*丁二烯的质量流量*/r=373000, /*冷凝潜热*/a2=950,a1, /*冷凝给热系数,管内给热系数*/T2=25,T1=15, /*冷却水出，进口温度*/t2=40,t1=40, /*冷凝温度*/d1=0.02,d2=0.025, /*内径，外径*/u=1.5, /*水的经济流速*/P=998.2, /*水的密度*/i=0.001005, /*水的粘度*/Pr, /*普朗克常数*/Cp=4183, /*水的比热容*/x=0.599, /*管内摩擦因数*/Nu, /*努赛尔常数*/ b=0.0025, /*壁厚*/X=45, /*钢管导热率*/k,K, /*假设的总传热系数，题目中的理论传热系数*/N,n, /*管程数与单程管数*/ms1=18.577, /*水的质量流量*/t,p,Q,tm,dm,A,a,l=6,L,Re, /*辅助参数*/Rs1=5.8e-4,Rs2=1.76e-4;/*水侧与丁二烯侧的热阻*/ /*变量定义*/Q=ms2*r;printf("n在本问题中：钢管的导热率取45，水的经济流速取1.5m/s。nn理论需要值得计算：nn总热量：Q=%.3fJn",Q);tm=19.58;printf("n平均温差：tm=%.3f摄氏度n",tm);Re=d1*u*P/i;printf("n雷诺数：Re=%.3fn",Re);Pr=Cp*i/x;printf("n普朗特常数：Pr=%.3fn",Pr);Nu=0.023*pow(Re,0.8)*pow(Pr,0.4);printf("n努塞尔数：Nu=%.3fn",Nu);a1=x*Nu/d1;printf("n管内给热系数：a1=%.3fn",a1);dm=(d1+d2)/2;printf("n管子的平均半径：dm=%.3f米n",dm);K=1/(d2/(a1*d1)+Rs1*d2/d1+b*d2/(X*dm)+Rs2+1/a2);printf("n总的给热系数：K=%.3fn",K);A=Q/(tm*K);printf("n换热面积：A=%.3f平方米nnn",A);int y=1;for(k=390;(K/k<=1.3&&K/k>=1.1);k+)a=Q/(tm*k);/*其中a是设出总的传热系数后k后计算出的面积*/n=4.0*ms1/(P*0.02*0.02*3.14*u);/*单程管数的计算*/L=a/(3.14*d2*n);N=L/l;p=(a-A)/A;if(p<=0.3)&&(p>=0.1)printf("nnn实际设计值%d：nn当1.1<K/k<1.3,面积余量在0.10.3时根据题目设计出的参数：nn换热面积： A=%.3f平方米nn单程管数： n=%.0fnn管程数： N=%.0fnn总的给热系数：K=%.3fnnn",y+,a+0.5,n,N+0.5,k);
