📄 declare.pas

📁 《SESyetem节能模拟系统(版本:1.01a Final)》为自由软件,对热风炉进行热平衡计算、传热模拟、燃烧优化、操作优化等设计。本软件用Dephi6.0编写
💻 PAS
📖 第 1 页 / 共 5 页
字号:
  end;

Function AirExpendCoefficient(O2y:Real;COy:Real;H2y:Real;CH4y:Real;
                              N2y:Real;N2:Real;RO2y:Real;CO2:Real;
                              CO:Real;CH4:Real;m:Integer;CmHn:Real;
                              H2S:Real):Real;
var x,y:Real;
   begin
   x:=O2y-0.5*COy-0.5*H2y-2*CH4y;
   y:=N2y-N2*(RO2y+COy+CH4y)/(CO2+CO+CH4+m*CmHn+H2S);
   Result:=21/(21-79*x/y);
   end;

Function IncompleteComparecomplete(a:Real;COy:Real;H2y:Real;
                                   CH4y:Real;O2y:Real):Real;
   begin
   if a>=1 then Result:=100/(100-0.5*COy-0.5*H2y) else
   Result:=100/(100+1.88*COy+1.88*H2y+9.52*CH4y-4.762*O2y);
   end;

Function TheoryDryAirQ(H2:Real;CO:Real;CH4:Real;m:Integer;
                       n:Integer;CmHn:Real;H2S:Real;O2:Real):Real;
   begin    
   Result:=0.0238*(H2+CO)+0.095*CH4+0.0476*(m+n/4)*CmHn+0.0714*H2S-0.0476*O2;
   end;

Function TheoryHumidityFumeQ(CO:Real;CH4:Real;m:Integer;n:Integer;
                             CmHn:Real;CO2:Real;H2:Real;H2S:Real;
                             N2:Real;H2O:Real;Log:Real):Real;
   begin
   Result:=0.01*(CO+3*CH4+(m+n/2)*CmHn+CO2+H2+2*H2S+N2+H2O)+0.79*Log;
   end;

Function FactHumidityFumeQ(V0:Real;a:Real;gk:Real;Log:Real):Real;
   begin    
   Result:=V0+(a*(1+0.00124*gk)-1)*Log;
   end;

//某一温度下的比热计算公式    单位：J/(kg.℃)
   //采用沈维道编的《工程热力学》P499的公式
Function CpO2(T:Real):Real;
   begin
   Result:=1e+3*(25.48+15.2*1e-3*power(T+273,1)+5.062*1e-6*power(T+273,2)+1.312*1e-9*power(T+273,3))/32;
   end;
Function CpCO2(T:Real):Real;
   begin
   Result:=1e+3*(22.26+59.811*1e-3*power(T+273,1)-35.01*1e-6*power(T+273,2)+7.47*1e-9*power(T+273,3))/44.09;
   end;
Function CpN2(T:Real):Real;
   begin
   Result:=1e+3*(28.184+0.846*1e-3*power(T+273,1)+5.661*1e-6*power(T+273,2)-2.145*1e-9*power(T+273,3))/28.016;
   end;
Function CpH2O(T:Real):Real;
   begin
   Result:=1e+3*(32.24+19.24*1e-3*power(T+273,1)+10.56*1e-6*power(T+273,2)-3.595*1e-9*power(T+273,3))/18.016;
   end;
Function CpH2Oq(T:Real):Real;
   begin
   Result:=1e+3*(30.01+10.7*1e-3*power(T+273,1)+0.33*1e+5*power(T+273,-2))/18;
   end;
Function Cpair(T:Real):Real;
   begin          //空气比热计算公式
   Result:=1e+3*(28.15+1.967*1e-3*power(T+273,1)+4.801*1e-6*power(T+273,2)-1.966*1e-9*power(T+273,3))/28.97;
   //Result:=6.761+1.0158*1.0e-3*(T+273.15)-0.084*1.0e5/power((T+273.15),2)+0;
   end;
Function CpCO(T:Real):Real;
   begin
   Result:=1e+3*(28.16+1.675*1e-3*power(T+273,1)+5.372*1e-6*power(T+273,2)-2.22*1e-9*power(T+273,3))/28.01;
   end;
Function CpH2(T:Real):Real;
   begin
   Result:=1e+3*(29.21-1.916*1e-3*power(T+273,1)-4.004*1e-6*power(T+273,2)-0.8705*1e-9*power(T+273,3))/2.016;
   end;
Function CpCH4(T:Real):Real;
   begin
   Result:=1e+3*(19.89+50.24*1e-3*power(T+273,1)+12.69*1e-6*power(T+273,2)-11.01*1e-9*power(T+273,3))/16.04;
   end;
Function CpC2H4(T:Real):Real;
   begin
   Result:=1e+3*(4.026+155.0*1e-3*power(T+273,1)-81.56*1e-6*power(T+273,2)+16.98*1e-9*power(T+273,3))/32;
   end;
Function QdwEnergy(CO:Real;H2:Real;CH4:Real;C2H6:Real;C2H4:Real;H2S:Real):Real;
   begin      //低位发热量 Qdws  单位 kcal/m^3
   Result:=30.2*CO+25.8*H2+85.7*CH4+152*C2H6+143*C2H4+56*H2S;
   end;
Function CoalNeedQ(Vm:Real;Tr:Real;Vf:Real;Tf:Real):Real;
   begin         //煤气用量,单位：m^3煤气/m^3热风
   Result:=Vm*Tr/Vf/Tf;
   end;
Function FuelPhysicsEnergy(B:Real;Cm:Real;tm:Real;Cme:Real;te:Real):Real;
   begin       //燃料的物理热量
   Result:=B*(Cm*tm-Cme*te);
   end;
Function TheoryHumidityAirQ(Log:Real;gk:Real):Real;
   begin      //理论湿空气量
   Result:=Log*(1+0.00124*gk);
   end;
Function HelpAirPhysicsEnergy(B:Real;a:Real;Los:Real;Ck:Real;tk:Real;Cke:Real;te:Real):Real;
   begin   //助燃空气的物理热
   Result:=B*a*Los*(Ck*tk-Cke*te);
   end;
Function CoolWindEnergy(Cf1:Real;tf1:Real;Cfe:Real;te:Real):Real;
   begin    //冷风带入的热量
   Result:=Cf1*tf1-Cfe*te;
   end;

Function HotStoveEnergy(Cf2:Real;tf2:Real;Cfe:Real;te:Real):Real;
   begin      //热风炉带出的热量
   Result:=Cf2*tf2-Cfe*te;
   end;
Function FumePhysicsEnergy(B:Real;bw:Real;Vns:Real;Cy2:Real;ty2:Real;Cye:Real;te:Real):Real;
   begin   //烟气带出的物理热量
   Result:=B*bw*Vns*(Cy2*ty2-Cye*te);
   end;
Function ChemistryIncompleteEnergy(B:Real;bw:Real;Vns:Real;Qdwsh:Real):Real;
   begin    //化学不完全燃烧热损失的热量
   Result:=B*bw*Vns*Qdwsh;
   end;
Function CoalMachineWaterEnergy(B:Real;gmj:Real;tm:Real;Cq:Real;ty2:Real):Real;
   begin    //煤气机械水的吸热量
   Result:=B*gmj*(100-tm+539+1.244*(Cq*ty2-36))*1e-3;
   end;
Function CoolWaterEnergy(Gs:Real;T:real;Vf:Real;Tf:Real;ts2:Real;ts1:Real):Real;
   begin    //冷却水的吸热量
   Result:=Gs*T*(ts2-ts1)/Vf/Tf;
   end;
Function CoolPipeEnergy(T:Real;Vf:Real;Tf:Real;qi:Real;Ai:Real):Real;
   begin     //冷风管道表面散热量
   Result:=T*qi*Ai/Vf/Tf;
   end;
Function RadiateEnergy(e:Real;tb:Real;te:Real;ad:Real):Real;
   begin
   Result:=4.88*e*(power((273+tb)/100,4)-power((273+te)/100,4))+ad*(tb-te);
   end;
Function ConvectionEnergy(A:Real;tb:Real;te:Real):Real;
   begin
   Result:=A*power(tb-te,0.25);
   end;
Function Cpy(CO2y:Real;N2y:Real;O2y:Real;CH4y:Real;COy:Real;H2Oy:Real;t:Real):Real;
   begin
     Result:=0.92*1e-5*(CO2y*CpCO2(t)*44/22.4+N2y*CpN2(t)*28/22.4+O2y*CpO2(t)*32/22.4+CH4y*CpCH4(t)*16/22.4+COy*CpCO(t)*28/22.4+H2Oy*CpH2Oq(t)*18/22.4);
   end;
Function space(n:Integer):String;
var i:integer;
var s:String;
begin
s:='';
for i:=1 to n do s:=s+' ';
Result:=s;
end;

//模拟
Function AverageRadLen(V:Real;F:Real):Real;
  begin
   Result:=3.6*V/F;
  end;

Function EmissivityCO2(PCO2:Real;T:Real;l:Real):Real;
   begin
     Result:=4.07*Power(PCO2*l,1/3)/(5.67*Power((T+273)/100,0.5));
   end;
Function EmissivityH2O(PH2O:Real;T:Real;l:Real):Real;
   begin
     Result:=40.7*Power(PH2O,0.8)*Power(l,0.6)/(5.67*Power((T+273)/100,1));
   end;
Function AbsorptivityCO2(eCO2:Real;T:Real;ts:Real):Real;
  begin
    Result:=eCO2*Power((T+273)/(ts+273),0.65);
  end;
Function AbsorptivityH2O(eH2O:Real;T:Real;ts:Real):Real;
  begin
    Result:=eH2O*Power((T+273)/(ts+273),0.45);
  end;
Function RadiationEnergy(es:Real;e:Real;T:Real;A:Real;ts:Real):Real;
   begin       //气体与固体间的辐射换热量
     Result:=(es+1)/2*5.67*(e*Power((T+273)/100,4)-A*Power((ts+273)/100,4));
   end;
Function ConvectionCoefficient(de:Real;r:Real;Pr:Real;w:Real;v:Real):Real;
   begin  //对流传热系数  单位：W/(m^2.C)
     Result:=0.03*r*Power(Pr,0.4)*Power(w,0.8)/(Power(de,0.17)*Power(v,0.83));
   end;
Function RadiationCoefficient(qr:Real;T:Real;Ts:Real):Real;
   begin    //辐射传热系数  单位：W/(m^2.C)
    if T=Ts then Result:=0 else  Result:=qr/(T-Ts);
   end;
Function FumeP(T:Real;p:Real):Real;
  begin    //实际烟气密度
  Result:=273*p/(T+273);
  end;
Function FumeW(Lk:Real;Af:Real):Real;
   begin                  //烟气实际速度
   Result:=Lk/(3600*Af); //计算速度 ，m/s
   end;
Function FumeCp(T:Real):Real;
var x,y:Real;   //烟气比热
  begin
  x:=SESystemElement.o2_fume[2]*32*CpO2(T)+SESystemElement.co2_fume[2]*44*CpCO2(T)+28*SESystemElement.n2_fume[2]*CpN2(T)+18*SESystemElement.h2o_fume*CpH2Oq(T);
  y:=SESystemElement.o2_fume[2]*32+SESystemElement.co2_fume[2]*44+28*SESystemElement.n2_fume[2]+18*SESystemElement.h2o_fume;
  Result:=x/y;
  end;
Function uO2(T:Real):Real;
   begin
    Result:=19.2*1e-6*(273+125)*power(T/273,1.5)/(T+125);
   end;
Function uCO2(T:Real):Real;
   begin
    Result:=13.8*1e-6*(273+254)*power(T/273,1.5)/(T+254);
   end;
Function uN2(T:Real):Real;
   begin
    Result:=16.6*1e-6*(273+104)*power(T/273,1.5)/(T+104);
   end;
Function uH2O(T:Real):Real;
   begin
    Result:=8.93*1e-6*(273+961)*power(T/273,1.5)/(T+961);
   end;
Function uair(T:Real):Real;
   begin
    Result:=17.09*1e-6*(273+111)*power(T/273,1.5)/(T+111);
   end;
Function pair(T:Real):Real;
  begin   //实际空气密度
  Result:=273*1.293/(T+273);
  end;
Function rO2(T:Real):Real;
   begin
   Result:=0.0245*power((T+273.16)/273.16,0.87);
   end;
Function rCO2(T:Real):Real;
   begin
   Result:=0.0145*power((T+273.16)/273.16,1.23);
   end;
Function rN2(T:Real):Real;
   begin
   Result:=0.0242*power((T+273.16)/273.16,0.8);
   end;
Function rH2O(T:Real):Real;
   begin
   Result:=0.015*power((T+273.16)/273.16,1.45);
   end;
Function rair(T:Real):Real;
   begin
   Result:=0.0236*power((T+273.16)/273.16,0.82);
   end;
Function FumeU(T:Real):Real;
var x,y:Real;
   begin
   x:=SESystemElement.o2_fume[2]*32*uO2(T)+SESystemElement.co2_fume[2]*44*uCO2(T)+28*SESystemElement.n2_fume[2]*uN2(T)+18*SESystemElement.h2o_fume*uH2O(T);
   y:=SESystemElement.o2_fume[2]*32+SESystemElement.co2_fume[2]*44+28*SESystemElement.n2_fume[2]+18*SESystemElement.h2o_fume;
   Result:= x/y;
   end;
Function FumeR(T:Real):Real;
var x,y:Real;
  begin
  x:=SESystemElement.o2_fume[2]*32*rO2(T)+SESystemElement.co2_fume[2]*44*rCO2(T)+28*SESystemElement.n2_fume[2]*rN2(T)+18*SESystemElement.h2o_fume*rH2O(T);
  y:=SESystemElement.o2_fume[2]*32+SESystemElement.co2_fume[2]*44+28*SESystemElement.n2_fume[2]+18*SESystemElement.h2o_fume;
  Result:=x/y;
  end;
Function FumePr(Cp:Real;u:Real;r:Real):Real;
   begin
   Result:=Cp*u/r;
   end;

Function SolidX(h:Real;ts:Real):Real;
var x:Real;   // 固体导热系数
  begin
  x:=(1.68-0.23*1e-3*ts);
  if h<12 then
    begin
    x:=(1.68-0.23*1e-3*ts);
    end else
    if (h<24) and (h>=12)  then
    begin
    x:=1.52-0.186*1e-3*ts;
    end else
  if(h>=24) and (h<=36) then x:=0.84+0.58*1e-3*ts else
  if h>36 then
    begin
    x:=0.84+0.58*1e-3*ts;
    end;
    Result:=x;
  end;
Function TotalTransferCoefficient(a:Real;hd:Real;r:Real):Real;
  begin    //
   Result:=1/(1/a+hd/6/r);
  end;
Function C1(az:Real;A:Real;W:Real;C:Real;L:Real;yc:Real;u:Real):Real;
  begin
   Result:=az*A*yc/(2*u*W*C*L);
  end;
Function C2(az:Real;A:Real;M:Real;Cs:Real;Tc:Real):Real;
  begin
   Result:=az*A*Tc/(2*M*Cs);
  end;
Function RegeneratorCp(h:Real;ts:Real):Real;
var x:Real; //耐火材料的比热  kJ/(kg.c)
  begin
    x:=1e+3*(0.837+0.234*1e-3*ts);
   if h<12 then
    begin
    x:=1e+3*(0.837+0.234*1e-3*ts);
    end else
    if (h<24) and (h>=12)  then
    begin
    x:=1e+3*(0.837+0.234*1e-3*ts);
    end else
    if (h>=24) then x:=1e+3*(0.837+0.264*1e-3*ts);
   Result:=x;
end;
Function LossEnergy(Aws:Real;qw:Real;W:Real;C:Real;yc:Real):Real;
  begin
   Result:=Aws*qw*yc/(W*C);
  end;

Function AirW(Lk:Real;Af:Real):Real;
   begin                  //烟气实际速度
   Result:=Lk/(3600*Af); //计算速度 ，m/s
   end;
Function AirPr(Cp:Real;u:Real;r:Real):Real;
   begin
   Result:=Cp*u/r;
   end;
Function A1(C1:Real):Real;
  begin
   Result:=(1-C1)/(1+C1);
  end;
Function A2(C1:Real):Real;
  begin
   Result:=C1/(1+C1);
  end;
⌨️ 快捷键说明

复制代码 Ctrl + C
搜索代码 Ctrl + F
全屏模式 F11
切换主题 Ctrl + Shift + D
显示快捷键 ?
增大字号 Ctrl + =
减小字号 Ctrl + -