l.c

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

#include <stdio.h>
#include <math.h>
#define T1 385.5
#define I1 3742.5
#define igs 485.93
#define tgs 115.0
#define Qw 100.85
#define i1 586.8
#define t1 138.9
#define Xa 0.02
#define H 697.7
#define Bj 6.0863
#define a 0.9936
#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,ku,ksp,kps,e,ad,a1,ea1;
float af,aff,K;
    t2=186.5;
	i2=(794.22+789.77)/2.0;
	D=36.1;
	Qs=1.02*D*(i2-i1)/Bj;
    I2=I1-Qs/a+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);
	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));
    ksp=(kr+ku)*0.0133;
    kps=exp(-ksp);
	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;
    ad=78.5;
    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("i2=%6.0f,I2=%6.1f,T2=%3.1f,K=%4.0f\n",i2,I2,T2);
	printf("/ntd=%3.1f,tx=%3.1f,tlog=%2.6f\n",td,tx,tlog);
	printf("/nt1=%4.1f,Qs=%5.2f,Xt=%4.0f,Q=%5.2f,Xe=%3.1f",t1,Qs,Xt,Q,Xe);

}