动力学计算软件使用说明书2.docVIP

动力学计算软件使用说明书 一 原理 本动力学计算软件主要根据胡荣祖、史启祯主编的《热分析动力学》一书提供的计算模型使用C++语言编写而成，主要采用迭代和最小二乘算法进行数据处理，根据所得直线的线性程度、偏差的大小以及普适积分法和微分方程法所求得的活化能及指前因子的一致性进行综合分析判断，并且利用多重速率扫描Kissinger法所求得的活化能进一步进行核实，最后确定机理函数，求得动力学各参数及动力学方程。 二 计算模型 1 普适积分法： 由 为了用积分方程逻辑选择出最概然机理函数，并求出相应的动力学参数E、A，为此把方程（4）改写成如下对数形式： 由一条TG或DSC曲线可得到原始数据：Ti、αi（i=1，2，…m）、β和T0，利用这些数据和线性最小二乘法处理方程，由斜率求E，截距求A。其中G(α)选取表1中所列的各种不同形式。 2 微分方程法 对方程（4）两边微分，整理得 将方程（7）代入方程（6），且两边取对数，则得计算非定温反应动力学的微分方程模型（8）  函数序号 积分形式机理函数G(α) 微分形式机理函数f(α) 1 2 3 4,5 (n=2,1/2) 6 7 8 9 10～16 （n=2/3,1/2.1/3,4,1/4,2,3） n=(2/3,1/2,1/3,4,1/4,2,3) 17～22 (n=1/2,3,2,4,1/3,1/4) (n=1/2,3,2,4,1/3,1/4) 23～27 （n=1,3/2,1/2,1/3,1/4） （n=1,3/2,1/2,1/3,1/4） 28 29 30 方程（8）左边与1/T呈线性关系，对于表1中给出的每一个微分形式的机理函数f(α)，方程（8）可用迭代法求解。 给定任意一个大于零的E值，利用这个值对每个数据点都可以计算出左边表达式相应的值，然后利用线性最小二乘法从斜率得到新的E值，而从截距得到新的A值。把E的这个修正值作为初值，再次迭代，可得到 另一个修正值，这样经几次迭代后，就会得到较合理的E和A值。对表1中的每一个机理函数f(α)都做上述处理，最后找出相关系数大，偏差小、且逻辑上合理的机理函数f(α)。 具体计算公式如下： 令： 则方程（5）变为 其中，a = E/R，b = ln(A/β)，L为数据点的个数。 用最小二乘法解方程组（6）得 这里，， 计算出a和b后就可得到E = aR，A = βeb。相关系数r和剩余方差Q按下式计算 由一条TG或DSC曲线可得到原始数据：Ti、αi（i=1，2，…m）、T0、β和dα/dT，利用这些数据并选取表1中所列的f(α)各种不同形式，即可求出动力学各参数。 3 通用Kissinger法 由第一类动力学方程 出发，当T = Tp 时，有 将 代入（10）式，整理，取对数，得 作－关系图，由直线斜率求E. 三 最概然机理函数的推断 由张同来－胡荣祖－杨正权－李富平法提出的选择合理动力学参数及最概然机理函数的五条判据是： 1 用普适积分法和微分方程法求得的动力学参数E和A值应在材料热分解反应动力学参数值的正常范围内，即活化能E值在80～250kJ/mol之间，指前因子的对数（lgA/s-1）值在7～30之间； 2 用微分法和积分法计算结果的线性相关系数要大于0.98； 3 用微分法和积分法计算结果的标准偏差应小于0.3； 4 根据上述原则选择的机理函数f(α)应与研究对象的状态相符； 5 与两点法、Kissinger法、Ozawa法和其它微分法求得的动力学参数值应尽量一直。 四 使用方法 本软件支持文本文件输入数据，输入时其排列格式依次为（dα/dT，T/K，αi）；升温速率β采用K·s-1为单位，温度采用国际单位，活化能采用焦耳为单位，指前因子的单位为s-1（或者可以不标注）。结果输出为积分法和微分法所对应的每种机理函数所求得的活化能、指前因子、相关系数、标准偏差和两种方法活化能之差，以及Kissinger法所求得的活化能。，结果分别可以输出到临时文本下（当前目录）下，和当前运行程序下，通过连续击打Enter键，不断输出每个函数对应求得的结果。 Huazhong Univers