煤热解模型预测

 %test 计算PE（0）
Fb0=0.4412;
P0=0.897;
T0=650; %初始温度
beta=2000; %升温速率
%不稳桥在全部链接中的概率
PB=Fb0*P0;
%计算反应物中不稳定桥的概率
w=[P0Fb0]; %w（1）---p（0）--完整链接在全煤中的分数   w（2）--Fb（0）--不稳定桥在完整链接中所占比例
Prb0=PRb（w）;
ww=[Prb0P0];
BR_0=BR（ww）;
%计算末端含有官能团的概率
%x=[275101Fb0*P01.8560.9110.8387] %x（1）--MWmon-单体单元分子量   x（2）--MWA--核的分子量  x（3）--pb--不稳定桥占完整链接比例  x（4）--MWB/MWA   x（5）--p    x（6）--MWC/MWA%
%Pe0=pe（x）
Pe0=1;
%计算反应物中的外围结构
y=[Pe0P0]; %y（1）--pe（0）  y（2）--p（0）
SR0=SR（y）; 
%计算后形质体的摩尔数

m5_0=m5（P0）;
m4_0=m4（P0）;
m3_0=m3（P0）;
m2_0=m2（P0）;
m1_0=（1-P0）^2;

%计算后形质体的不稳定桥

z=[Fb0P0];
BM5_0=BM5（z）;
BM4_0=BM4（z）;
BM3_0=BM3（z）;
BM2_0=BM2（z）;
BM1_0=0;

%计算SM
SM0=SM（y）;

%计算中间体
n10_0=n10（P0）;
n9_0=n9（P0）;
n8_0=n8（P0）;
n7_0=n7（P0）;
n6_0=n6（P0）;

%计算中间体的不稳定桥
BN10_0=BN10（z）;
BN9_0=BN9（z）;
BN8_0=BN8（z）;
BN7_0=BN7（z）;
BN6_0=BN6（z）;
%计算中间体的外围结构
SN_0=SN（y）;

%初始非冷凝气体 表III  fuel and energy 55-1991663
YG_0=0;
%计算KB不稳定桥断裂速率常数的表达式及值 在T温度的值
t=[0:0.014:0.275];
%T=t*beta+T0;
%kB=kB_splinet（T）;
%外围结构脱除反应速率
%kG=1d15*exp（（-230000）*power（8.314*T-1））; %kG--气体逸出速率常数表达式