1热平衡计算.c

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

#include <stdio.h>
#include <math.h>
#define Qnet 18359.0
#define tlk 30
#define ilk 193
#define q3 0.5
#define q4 1.5
#define q5 0.66
#define igr 3332.4
#define ibh 1109.5
#define igs 485.93
main()
{
	float q2,ggl,Qgl,B,Bj,Qky,qky,y;
	float Tpy,Ipy,trk,irk;
    Tpy=155;                        /*假定*/
	Ipy=967.0+0.01*(Tpy-100.0)*991;
	q2=(Ipy-1.33*ilk)*(100.0-q4)/Qnet;
	ggl=100.0-(q4+q3+q2+q5);
	Qgl=36.1*(igr-igs)+0.02*36.1*(ibh-igs);
	B=100.0*Qgl/(ggl*Qnet);
	Bj=B*(1.0-q4*0.01);
	trk=335.0;
	irk=1958.0+(2634.0-1958.0)*0.01*(trk-300.0);
	Qky=1.06*(irk-ilk);
	qky=Qky*100.0/Qnet;
	y=1.0/(1.0+q5/(ggl+qky));
	printf("Tpy=%4.1f,Ipy=%4.1f,q2=%2.2f\n",Tpy,Ipy,q2);
	printf("ggl=%2.2f,Qgl=%6.1f,B=%1.3f\n",ggl,Qgl,B);
	printf("Bj=%1.4f,Qky=%4.1f,qky=%2.2f,y=%1.4f\n",Bj,Qky,qky,y);
}