9下级省煤器.c

这个源码是从其它网站上转载来的

#include <stdio.h>
#include <math.h>
#define T1 375.6
#define I1 3643.1
#define igs 485.93
#define tgs 115.0
#define Xa 0.02
#define H 697.7
#define Bj 6.1115
#define Ilk 193.0
#define Vy 6.709
#define Fy 11.0
main()
{
    float t2,i2,Qs,I2,T2,td,D,tx,Xt,Q,Xe,h1,h2,tlog,Tpj;
	float Tpjj,Wy,tpj,tb,kr,Qw,ku,kps,e,ad,a1,ea1;
float af,aff,K,XD,i1,t1,a0,a;
    D=36.1;
    XD=2.4;
	Qw=(XD*1688)/(1.02*D);
	i1=igs+Qw;
    t1=142+(i1-597.8)/(602.1-597.8);
    t2=186;              /*假定*/
	i2=794.22;  /*临时计算*/
	Qs=1.02*D*(i2-i1)/Bj;
    I2=I1-Qs/0.9936+Xa*Ilk;
	h1=1919.0;
	h2=2915.0;
    T2=100.0*(I2-h1)/(h2-h1)+200.0;
	Tpj=(T1+T2)/2;
	Tpjj=Tpj+273.0;
    Wy=Bj*Vy*Tpjj/(Fy*273.0);
    a0=81;                        /*临时计算*/
	ad=a0*1*1*1;
	tpj=(t1+t2)/2;
    tb=tpj+25.0;
    kr=8.96*(1-370.0*Tpjj)/(1.0-370.0*590.2);
    ku=4.0*(pow(Tpjj,-0.6667))/(pow(590.2,-0.6667));
    kps=(kr+ku)*0.0133;
    a=1.0-(exp(-kps));
	af=33.7*(1.0-exp(-kps));
	aff=af*(1+0.4*pow(Tpjj/1000,0.25)*pow(1/1.8,0.07));
    e=0.0027;
    a1=ad+aff;
	ea1=1/(1+e*a1);
	K=a1*ea1;
	td=T1-t2;
	tx=T2-t1;
    tlog=log(td/tx);
	Xt=(td-tx)/tlog;
	Q=K*H*Xt/(1000.0*Bj);
	Xe=100.0*(Qs-Q)/Qs;
	printf("a=%f,af=%f,aff=%f,a1=%f,t1=%f,t2=%f\n",a,af,aff,a1,t1,t2);
	printf("Qw=%f,i1=%f,i2=%6.1f,I2=%f,T2=%f,K=%4.1f\n",Qw,i1,i2,I2,T2,K);
	printf("Wy=%f,Qs=%f,Tpj=%f,tpj=%f,tb=%f,kr=%f,ku=%f,kps=%f\n",Wy,Qs,Tpj,tpj,tb,kr,ku,kps);
	printf("ad=%f,td=%3.1f,tx=%3.1f,tlog=%2.6f\n",ad,td,tx,tlog);
	printf("/nt1=%4.1f,Qs=%5.2f,Xt=%4.1f,Q=%5.2f,Xe=%3.2f",t1,Qs,Xt,Q,Xe);

}