Курсовая работа: Математическое моделирование тепловой работы вращающейся печи

IF TF=0 THEN CT:=2.522+0.0005815*TF ELSE

CT:=(-317.9+2.522*TF+0.0005815*sqr(TF)+86840/(TF+273))/TF;

{Объемная теплоемкость воздуха, кДж/м3*К (н.у.)}

CB:=1.287+0.0001201*TF;

{Объёмная теплоёмкость углекислого газа, кДж/(м3*К) (н.у.)}

IF TF=0 THEN C1:=2.081+0.0002017*TF ELSE

C1:=(-139.5+2.081*TF+0.0002017*sqr(TF)+38110/(TF+273))/TF;

{Объёмная теплоёмкость водяного пара, кДж/(м3*К) (н.у.)}

C2:=1.49+0.0002303*TF;

{Объёмная теплоёмкость азота, кДж/(м3*К) (н.у.)}

C3:=1.28+0.0001103*TF;

{Объёмная теплоёмкость продуктов горения, кДж/м3*К (н.у.)}

CG:=C1*P1+C2*P2+C3*P3;

{Объёмная теплоёмкость газов в факеле, кДж/(м3*К) (н.у.)}

CV:=CT*PT+CB*PB+CG*PG;

{Массовая теплоёмкость газов в факеле, кДж/кг*К}

CF:=CV/ROV;

End;

{============================================================================}

Function tFurnace.fHdis:real;

{============================================================================}

{Расчеттеплотыдиссоциациифакела}

var PCO2, PH2O, KCO2, KH2O, ACO2, AH2O, PS, LR, K1, K2: real;

Begin

fHdis:=0;

if TF>1500 then begin

{Объёмные доли трехатомных газов}

PCO2:=P1*PG; {…углекислого газа}

PH2O:=P2*PG; {…водяного пара}