5一级过热器计算.c

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

#include <stdio.h>
#include <math.h>
#define Bj 6.1115
#define D 36.1
#define T1 817.1
#define I1 8021.3
#define t1 257
#define i1 2809
#define ibh 1109.5
main()
{
	float T2,I2,t2,i2,XD,ig,Qq,Tpj,wy,ad,a0,tpj,v,wq,a0q,a2,kr;
	float ku,kps,a,th,af,aff,a1,K,td,tx,Xt,Q,Xe;
	float h1,h2,v1,v2;
	ig=3016.6;
	XD=2.4;            /*假定*/
	i2=(ig*D-XD*ibh)/(D-XD);
	h1=3116.3-0.1*(3116.3-3093.7);
	h2=3140.8-0.1*(3140.8-3119.6);
	t2=10*(i2-h1)/(h2-h1)+360;
	Qq=(D-XD)*(i2-i1)/Bj;
	I2=I1-Qq/0.9936+0.015*193;
    T2=600+100*(I2-5787)/(6840-5787);
	Tpj=(T1+T2)/2;
	wy=Bj*6.377*(Tpj+273)/(273*15.63);
	a0=76;             /*临时计算*/
	ad=a0*1*1*0.96;
	tpj=(t1+t2)/2;
    v1=0.058821-0.1*(0.058821-0.045301);
	v2=0.060422-0.1*(0.060422-0.046723);
	v=v1+(tpj-300)*0.1*(v2-v1);
	wq=(D-XD)*v/0.095;
	a0q=1415;           /*临时计算*/
	a2=a0q*0.99;
    kr=(1.0-0.37*(Tpj+273)*0.001)*6.19/(1.0-0.37*0.961);
    ku=3.06*(pow(961,0.6667))/(pow((Tpj+273),0.6667));
	kps=(kr+ku)*0.0199;
	a=1.0-(exp(-kps));
	th=tpj+Qq*1000*Bj*(0.0048+1/a2)/470.7;
	af=118.6*a;
	aff=af*(1+0.4*pow((Tpj+273)*0.001,0.25)*pow((0.955/1.164),0.07));
    a1=ad+aff;
	K=0.65*a1*a2/(a1+a2);
	td=T1-t2;
	tx=T2-t1;
	Xt=(td-tx)/log(td/tx);
	Q=K*Xt*470.7/(1000*Bj);
	Xe=100*(Qq-Q)/Qq;
	printf("Qq=%f,i2=%f,t2=%f,a0q=%f\n",Qq,i2,t2,a0q);
	printf("I2=%f,T2=%f,Tpj=%f,ad=%f,tpj=%f,wy=%f,wq=%f\n",I2,T2,Tpj,ad,tpj,wy,wq);
	printf("ad=%f,v=%f,wq=%f,a2=%f,kr=%f,ku=%f,kps=%f,a=%f\n",ad,v,wq,a2,kr,ku,kps,a);
	printf("th=%f,af=%f,aff=%f,a1=%f,K=%f,Xt=%f,Q=%f,Xe=%f\n",th,af,aff,a1,K,Xt,Q,Xe);
}