《热分析动力学幻灯片.ppt》由会员分享，可在线阅读，更多相关《热分析动力学幻灯片.ppt（86页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、热分析动力学第1页，共86页，编辑于2022年，星期日l “What do bread and chocolate,hair and finger-nail clippings,coal and rubber,oint-ments and suppositories,explosives,kidney stones and ancient Egyptian papyri have in common?Many interesting answers could probably be suggested,but the connection wanted in this context is
2、that they all under-go interesting and practically important changes on heating”M.E.BrownIntroduction to Thermal Analysis:Techniques and Applications 第2页，共86页，编辑于2022年，星期日引引 言言 What?Why?When?Where?第3页，共86页，编辑于2022年，星期日What?定义和结果定义和结果 什么是热分析动力学什么是热分析动力学(KCE)？用热分析技术研究某种物理变化或化学用热分析技术研究某种物理变化或化学 反应反应(以下
3、统称反应以下统称反应)的动力学的动力学 热分析动力学获得的信息是什么？热分析动力学获得的信息是什么？判断反应遵循的机理、得到反应的动力判断反应遵循的机理、得到反应的动力学速率参数学速率参数(活化能活化能E E和指前因子和指前因子A等等)。即动力学即动力学 “三联体三联体”（kinetic triplet）第4页，共86页，编辑于2022年，星期日Why?条件和目的条件和目的l 为什么热分析能进行动力学为什么热分析能进行动力学 研究？研究？为什么要做动力学分析？为什么要做动力学分析？第5页，共86页，编辑于2022年，星期日物理性质物理性质（质量、能量等）（质量、能量等）温度（温度（T）过程进度
4、过程进度()时间（时间（t）=HT/H程序控温程序控温T=To+t 动力学关系动力学关系热分析：在程序控温下，测量物质的物理性质与热分析：在程序控温下，测量物质的物理性质与热分析：在程序控温下，测量物质的物理性质与热分析：在程序控温下，测量物质的物理性质与 温度的关系的一类技术温度的关系的一类技术温度的关系的一类技术温度的关系的一类技术 (5 5thth ICTA ICTA)条条 件件第6页，共86页，编辑于2022年，星期日目的目的v理论上：理论上：探讨物理变化或化学反应的机理（尤其探讨物理变化或化学反应的机理（尤其是非均相、不等温）是非均相、不等温）v生产上：生产上：提供反应器设计参数提供
5、反应器设计参数l应用上：应用上：建立过程进度、时间和温度之间的关系，建立过程进度、时间和温度之间的关系，可用于预测材料的使用寿命和产品的保质稳定期，可用于预测材料的使用寿命和产品的保质稳定期，评估含能材料的危险性，从而提供储存条件。此外评估含能材料的危险性，从而提供储存条件。此外可估计造成环境污染物质的分解情况可估计造成环境污染物质的分解情况第7页，共86页，编辑于2022年，星期日When?历史历史化学动力学化学动力学 源于源于1919世纪末世纪末2020世纪初世纪初热分析动力学热分析动力学 始于始于2020世纪世纪3030年代、盛于年代、盛于5050年代年代（主要应评估高分子材料在航空航天
6、应用（主要应评估高分子材料在航空航天应用 中的稳中的稳定性和使用寿命研究的需要）定性和使用寿命研究的需要）第8页，共86页，编辑于2022年，星期日Where?理论基础理论基础 等温、等温、均相均相 dT/dt=c 不等温、非均相不等温、非均相 Arrhenius 常数：常数：k(T)=Aexp(-E/RT)第9页，共86页，编辑于2022年，星期日1.1.回顾篇回顾篇 How?Idealized and Empirical Kinetic Models for Heterogeneous ReactionsMethodology of Kinetic Analysis第10页，共86页，编辑
7、于2022年，星期日 How?动力学模式动力学模式(机理机理)函数函数 均相反应：均相反应：f(c）=(1 c）n 非均相反应：非均相反应：根据控制反应速率的根据控制反应速率的“瓶颈瓶颈”v 气体扩散气体扩散 v 相界面反应相界面反应v 成核和生长成核和生长 第11页，共86页，编辑于2022年，星期日 均相反应均相反应(液相液相/气相气相)浓度浓度C C表示进程表示进程,级数反应级数反应非均相反应非均相反应(固体或固气反应固体或固气反应)转化率转化率表示进程表示进程引入相界面引入相界面与体积之比与体积之比速引率速引率入步控入步控骤制骤制气体扩散气体扩散相界面推进相界面推进反应物界面收缩反应物
8、界面收缩引引入入收收缩缩维维数数一维一维二维二维三维三维成核和生长成核和生长一维一维二维二维三维三维瞬瞬间间成成核核引入成核速率引入成核速率引引入入维维数数第12页，共86页，编辑于2022年，星期日 常见固态反应的机理函数（理想化）常见固态反应的机理函数（理想化）1.Acceleratory(The shape of a T curve)Symbol f(a)g(a)Pn n(a)1-1/n a 1/n E1 a lna2.Sigmoid Am m(1-a)-ln(1-a)1-1/m -ln(1-a)1/m B1 a(1-a)lna/(1-a)B2 (1/2)(1-a)-ln(1-a)-1
9、-ln(1-a)2 B3 (1/3)(1-a)-ln(1-a)-2 -ln(1-a)3 B4 (1/4)(1-a)-ln(1-a)-3 -ln(1-a)4第13页，共86页，编辑于2022年，星期日3.Deceleratory R2 2(1-a)2(1-a)1/21/2 1-(1-a)1-(1-a)1/21/2 R3 3(1-a)3(1-a)2/32/3 1-(1-a)1-(1-a)1/31/3 D1 1/2a a 1/2a a2 2 D2 1-ln(1-a)1-ln(1-a)-1-1 (1-a)ln(1-a)+a (1-a)ln(1-a)+a D3 (3/2)(1-a)(3/2)(1-a)2
10、/32/31-(1-a)1-(1-a)2/32/3-1-1 1-(1-a)1-(1-a)1/31/3 2 2 D4 (3/2)(1-a)(3/2)(1-a)-1/3-1/3-1-1-1-1 1-2a/3-(1-a)1-2a/3-(1-a)2/32/3 D5 (-3/2)(1-a)(-3/2)(1-a)2/32/3(1-a)(1-a)1/31/3-1-1-1-1 (1-a)(1-a)1/31/3-1-12 2 D6 (3/2)(1-a)(3/2)(1-a)4/34/3(1-a)(1-a)-1/3-1/3-1-1-1 -1 (1-a)(1-a)-1/3-1/3-1-12 2 第14页，共86页，编
11、辑于2022年，星期日F1*1-a -1n(1-a)F2 (1a)2 1/(1a)F3 (1a)3/2 (1/1-a)2 F(3/2)2(1-a)3/2 (1-a)-1/2 F(5/2)(2/3)(1-a)5/2 (1-a)-3/2*F1 is the same as A1 Sestak-Berggren empirical function(1971)f(a)=aa)=a m(1-a)1-a)n 第15页，共86页，编辑于2022年，星期日How?方法方法微分式：微分式：积分式积分式:第16页，共86页，编辑于2022年，星期日1.实验数据的准备实验数据的准备 TG：TWW0WWT第17页，
12、共86页，编辑于2022年，星期日DSC:HT/HTdH/dtHTH第18页，共86页，编辑于2022年，星期日 2.热分析方法热分析方法 等温等温(isothermal)法法不等温不等温(non-isothermal)法法 l 按动力学方程形式：按动力学方程形式：微商法微商法 积分法积分法v 按加热速率方式：按加热速率方式：单个扫描速率法单个扫描速率法 (single scanning method)多重扫描速率法多重扫描速率法 (multiple scanning method)(等转化率法，等转化率法，iso-conversional)第19页，共86页，编辑于2022年，星期日2-1
13、等温法等温法:2-1-1 模式适配法模式适配法(model-fitting method)a)测定几种不同测定几种不同T 下下(在该温度范围内反应能发生在该温度范围内反应能发生)的等温的等温 t 曲线。曲线。b)作作 t/t0.5 或或 t/t0.9 的约化时间图的约化时间图(reduced time plot,t0.5、t0.9 分别为分别为 0.5或或0.9 的时间的时间)，与文献报道，与文献报道的标准图的标准图(master plot)比较比较,判定最可几机理函数。判定最可几机理函数。c)根据上式计算在该温度下的根据上式计算在该温度下的k 值，如此重复可得一组值，如此重复可得一组 k 1
14、，T1；k 2，T2；k i，Ti；代入代入 由线性方程斜率由线性方程斜率 E;截矩截矩 A 第20页，共86页，编辑于2022年，星期日 2-1-2 等温等转化率法等温等转化率法 （isothermal isoconversional method)无需预先获得最可几机理函数（无需预先获得最可几机理函数（model-free)求取活化求取活化能能E值，且可得到活化能随着反应进程的关系（值，且可得到活化能随着反应进程的关系（E）选定某选定某值，则可从不同温度值，则可从不同温度T 的等温的等温 t 曲线中得曲线中得到对应于该到对应于该值的一组值的一组 t、T 数据，代入经两边取对数、数据，代入经
15、两边取对数、重排后得到的重排后得到的 因因在定值时，等式右边前两项为常数，则由斜率可在定值时，等式右边前两项为常数，则由斜率可求求 E。此式亦为一旦动力学三联体都获得后，建立时间此式亦为一旦动力学三联体都获得后，建立时间t、温度温度T 和分解百分数和分解百分数之间关系的基础之间关系的基础第21页，共86页，编辑于2022年，星期日 2-2 不等温法不等温法2-2-1 2-2-1 微商法微商法:l Kissinger-Akahira-Sunose(K-A-S)(1956)l Freeman-Carroll(1958)l Newkirk(1960)l Friedman(1964)l Achar-B
16、rindly-Sharp(A-B-S)(1966)第22页，共86页，编辑于2022年，星期日 Kissinger-Akahira-Sunose equation Anal.Chem.,29(1957)1702作多重加热速率作多重加热速率下的测定，选择下的测定，选择TA曲线峰值对应的温度曲线峰值对应的温度Tp 由线性方程斜率由线性方程斜率E，然后由截矩，然后由截矩A 注：注：1.Kissinger(1956):在最大速率处，适于在最大速率处，适于n级反应级反应 2.Akahira-Sunose(1969):指定指定处亦可处亦可 3.Ozawa:不限于不限于n级反应级反应第23页，共86页，编辑
17、于2022年，星期日 Freeman-Carroll equation J.Phys.Chem.,62(1958)394设设动力学方程微分式取对数，再用差值表示，则有：动力学方程微分式取对数，再用差值表示，则有：作图，由斜率作图，由斜率E;截矩截矩nAnderson-Freeman equation(1961)若取若取 为等值，则上式可简化为：为等值，则上式可简化为：第24页，共86页，编辑于2022年，星期日 Newkirk equation Anal.Chem.,32(12)(1960)1558若若 且且 n=1，则有：，则有：取两个实验点取两个实验点 T1和和T2,则有：则有：可求得可求
18、得 E第25页，共86页，编辑于2022年，星期日 Friedman equation J.Polym.Sci.Part C,6(1964)183作多重加热速率作多重加热速率下的测定，选择等下的测定，选择等处处 斜率斜率 E；截矩；截矩 若若则：则：斜率斜率E;截矩截矩 A第26页，共86页，编辑于2022年，星期日 Achar-Bridly-Sharp equationProc.Int.Clay Conf.Jerusalem,1(1966)67 Anal.Chem.41(1969)2060尝试不同的尝试不同的 f(),由线性方程由线性方程 斜率斜率 E 截矩截矩 A能获得最佳线性的能获得最佳
19、线性的 f()为最可几机理函数为最可几机理函数 第27页，共86页，编辑于2022年，星期日 2-2-2 2-2-2 积分法积分法:动力学方程积分式动力学方程积分式 TT0,=0 温度积分温度积分（Temperature Integral)为一非收敛级数，无精确解为一非收敛级数，无精确解 第28页，共86页，编辑于2022年，星期日 温度积分的近似表达式温度积分的近似表达式Doyle近似式近似式(J.Appl.Polym.Sci.,6(1962）639)Schlomlich展开级数（展开级数（Doyle,Nature,207(1965)290)经验公式（经验公式（Zsaco,J.Thermal
20、 Anal.8(1975)593)第29页，共86页，编辑于2022年，星期日Senum-Yang近似近似（J.Thermal Anal.11(1977)445)渐近展开级数渐近展开级数（Zsaco,Thermal Analysis,p167,1984)第30页，共86页，编辑于2022年，星期日 常用积分法常用积分法l Horowitz-Metzger equation(1963)l Coats-Redfern equation(1964)l Flynn-Wall-Ozawa equation(1965)l Zsaco equation(1968)l Maccallum-Tanner equ
21、ation(1968)l Satava-Sestak equation(1971)第31页，共86页，编辑于2022年，星期日 Horowitz-Metzger equation Anal.Chem.,35(10)(1963)1464 Tr 能使能使 1=1/e 的参考温度（一般为峰温）的参考温度（一般为峰温）T Tr 作作 ln g(）图，斜率图，斜率 E Dharwadkar-Karkhanavala 修正修正（J.Them.Anal.,2(1969)1049)：Ti，Tf 反应始、终温反应始、终温第32页，共86页，编辑于2022年，星期日Coats-Redfern equation N
22、ature，201（1964）68设设 则有：则有：（n=1)(n1)第33页，共86页，编辑于2022年，星期日 Flynn-Wall-Ozawa equation Bull.Chem.Soc.Jpn.,38(1965)1881 取不同取不同下曲线的等下曲线的等处之温度处之温度 T 作作 lg1/T图，由斜率图，由斜率 E 注：注：Ozawa(1965):在最大转化速率处在最大转化速率处 Flynn-Wall(1966):指定指定处亦可处亦可第34页，共86页，编辑于2022年，星期日 Zsaco equation J.Phys.Chem.,72(7)(1968)2406令令 P(x)用用
23、Zsaco近似，近似，B对于某反应和定对于某反应和定下为常数，假设下为常数，假设g(）的形式，并选定的形式，并选定E(如由如由F-W-O法获得），则：法获得），则：Bi 的平均值，的平均值，r 实验数据个数；实验数据个数；标准偏差标准偏差能使能使最小者，其最小者，其g(）和和E为最可能者为最可能者由由 求得求得 A 第35页，共86页，编辑于2022年，星期日 Maccallum-Tanner equation J.Eur.Polym.,4(1968)333 作作 lg g()1/T 图，由斜率图，由斜率 E第36页，共86页，编辑于2022年，星期日 Satava-Sestak equati
24、on Thermochim Acta,2(1971)423选择选择 g(），作，作 lg g(）1/T 图，最佳线性者为正确图，最佳线性者为正确之之 g(）由斜率由斜率 E,然后然后,由截矩由截矩 A第37页，共86页，编辑于2022年，星期日2.2.思考篇思考篇Uncertainties of kinetic resultskinetic compensation effect(KCE)Simple or multi-step reaction?第38页，共86页，编辑于2022年，星期日Problem 动力学结果的不确定性动力学结果的不确定性不等温法的优点不等温法的优点（Coats-Red
25、fern,Analyst,88(1963)906)1.1.要求的实验数据较之等温法少，快速方便要求的实验数据较之等温法少，快速方便;一条非一条非等温等温TATA曲线相当于多条等温曲线相当于多条等温TATA曲线。曲线。2.2.可以用一种连续的方式研究在整个温度范围内的可以用一种连续的方式研究在整个温度范围内的动力学过程。动力学过程。3.3.严格的等温法实际上的不可行性。严格的等温法实际上的不可行性。4.4.只用一个样品，可避免样品间的误差（包括样品只用一个样品，可避免样品间的误差（包括样品用量、粒度大小和形状、堆积方式等）用量、粒度大小和形状、堆积方式等）第39页，共86页，编辑于2022年，星
26、期日问题：问题：1.同一种样品的同一反应，有不同的动力学结果，而同一种样品的同一反应，有不同的动力学结果，而 且超出且超出了实验允许的误差分析范围之外。了实验允许的误差分析范围之外。如：如：CaCO3 的分解、的分解、PS 的降解等，的降解等，E 值值 2.当选择不同机理函数当选择不同机理函数 f()或或 g()代入各种动力学方法（微代入各种动力学方法（微商法和积分法），所得到的线性都很好，造成比较上的困商法和积分法），所得到的线性都很好，造成比较上的困难。难。3.即使不同的机理函数代入后线性略有不同，但是由于不即使不同的机理函数代入后线性略有不同，但是由于不同的模式函数所对应的动力学参数结果
27、相差很大，因此同的模式函数所对应的动力学参数结果相差很大，因此依赖线性判断机理函数，而后获得动力学参数的过程就依赖线性判断机理函数，而后获得动力学参数的过程就使人有使人有“失之毫厘，差之千里失之毫厘，差之千里”的担心。的担心。第40页，共86页，编辑于2022年，星期日80 100120140160abcdfgijklmeh12中不同方法处理中不同方法处理PS热分解热分解TG数据的数据的E值置信范围值置信范围F.Carrasco,Thermochim.Acta,213(1993)115EkJ/M第41页，共86页，编辑于2022年，星期日 ANDAND的不等温热分解反应动力学参数的不等温热分解
28、反应动力学参数模式模式 E/kJ mol-1-1 lnA/min-1-1 P4 24.5 3.9 0.9783P3 35.1 6.9 0.9813P2 56.2 12.7 0.9837P2/3 182.9 46.2 0.9862D1 246.2 62.8 0.9865F1 139.4 35.7 0.9928A4 29.5 5.3 0.9903A3 41.7 9.0 0.9913A2 66.1 15.9 0.9921D3 269.1 67.4 0.9928R3 131.0 32.0 0.9924R2 127.6 31.3 0.9910第42页，共86页，编辑于2022年，星期日 错误动力学三联体
29、错误动力学三联体 Vs 正确动力学三联体正确动力学三联体 F(),Eapp,Aapp f(),E,A Koga,et al,Thermochim.Acta,188(1991)333 p 在在 DSC/DTG 曲线峰温处的曲线峰温处的第43页，共86页，编辑于2022年，星期日Analysis 分析分析 1.动力学方程的缺陷：均相动力学方程的缺陷：均相非均相，等温非均相，等温不等不等温。（温。（20世纪世纪7080年代）年代）2.动力学补偿效应（动力学补偿效应（Kinetic Compensation Effect,KCE)的存在。（的存在。（20世纪世纪90年代）年代）3.固态反应（非均相反应
30、）本身的复杂性固态反应（非均相反应）本身的复杂性 (21世纪）世纪）第44页，共86页，编辑于2022年，星期日Kinetic Equation?Kinetic Compensation Effect?Complexity of Solid Reaction?第45页，共86页，编辑于2022年，星期日 动力学补偿效应（动力学补偿效应（KCE)It was shown that the kinetic exponent in the kinetic model function,the activation energy and the preexponential factor are mu
31、tually correlated.As a consequence of this correlation any TA curve can be by an apparent kinetic model instead of the appropriate one of a certain value of apparent Arrhenius parameters N.Koga et al Netsu Sokutei,20(1993)210N.Koga et al Netsu Sokutei,20(1993)210第46页，共86页，编辑于2022年，星期日 KCE的分类的分类 1.1.
32、存在于同一系列的物质（如分子母体结构存在于同一系列的物质（如分子母体结构相同的系列衍生物）的反应动力学结果相同的系列衍生物）的反应动力学结果中（第一类中（第一类 KCE)2.2.存在于同一物质，但不同实验条件下（如样存在于同一物质，但不同实验条件下（如样品用量、升温速率等）的反应动力学结果中品用量、升温速率等）的反应动力学结果中（第二类（第二类 KCE)3.3.存在于同一组数据，用不同数据处理方法存在于同一组数据，用不同数据处理方法的动力学结果中（第三类的动力学结果中（第三类 KCE?)?)第47页，共86页，编辑于2022年，星期日 From this point of view,it se
33、ems that the method of kinetic analysis aiming to ascertain all kinetic parameters from only one experimental TA curve are somewhat problematic.We have to realize that this problem cannot be solved even using the most sophisticated non-linear regression algorithms unless the kinetic models or at lea
34、st one kinetic para-meters is a priori known.J.Sestak&J.Malek,Solid State Ionics,63-65(1993)245 第48页，共86页，编辑于2022年，星期日The ICTAC kinetic project 在在第第1111届届国国际际热热分分析析和和量量热热学学会会议议(ICTAC，1996，美美国国费费城城)上上，由由动动力力学学工工作作组组（The Kinetic Workshop）发发起起，给给各各国国志志愿愿参参加加的的TA工工作作者者下下列列数数据据，然然后后由由参参与与者者自自选选数数据据处处理理方方
35、法法进进行行动动力力学分析：学分析：l 在真空和氮气气氛中，在不同扫描速率下，分别测得的在真空和氮气气氛中，在不同扫描速率下，分别测得的 CaCO3、NH4ClO4的热分解的热分解T曲线和数据曲线和数据l 在真空和氮气气氛中，在不同温度的等温条件下，分别在真空和氮气气氛中，在不同温度的等温条件下，分别 测得的测得的CaCO3、NH4ClO4的热分解的热分解t曲线和数据曲线和数据l 几种不同温度的等温模拟数据和几种不同加热速率的不几种不同温度的等温模拟数据和几种不同加热速率的不 等温模拟数据（反应机理设定为等权平行的一级反应）等温模拟数据（反应机理设定为等权平行的一级反应）第49页，共86页，编
36、辑于2022年，星期日结果与讨论结果与讨论(Thermochim.Acta,355(2000)125-180)(Thermochim.Acta,355(2000)125-180)1.1.同一反应在不同实验条件同一反应在不同实验条件(如气氛如气氛)下，甚至同一下，甚至同一套数据用不同处理方法所得到的动力学结果也不套数据用不同处理方法所得到的动力学结果也不同。同。2.2.单个加热速率法应该尽量避免使用，由于单个加热速率法应该尽量避免使用，由于KCE的存的存在，光凭线性优劣判断在，光凭线性优劣判断(model-fit)可以导致结果可以导致结果的不一致性。的不一致性。3.3.无模式无模式(model-
37、free)方法对于获得准确、一致的方法对于获得准确、一致的活化能和活化能和E(或或 E T)关系极为有效。关系极为有效。4.4.许多表面上简单许多表面上简单(single),),实际上却是包含了多步实际上却是包含了多步(multi-step)的复杂反应，判断这一点的初步方法的复杂反应，判断这一点的初步方法是考察其是考察其 E(或或 E T)关系。关系。第50页，共86页，编辑于2022年，星期日 Solution 解决办法解决办法 1.尝试引进机理函数分数指数的形式和兼容的经验机理尝试引进机理函数分数指数的形式和兼容的经验机理函数代替理想化的非均相反应机理函数。函数代替理想化的非均相反应机理函
38、数。2.多重扫描速率法（如多重扫描速率法（如F-W-O,K-A-S或或Fried-man公式）可在不公式）可在不知道机理函数的情况下知道机理函数的情况下(model-free)求取活化能求取活化能 E,摆脱摆脱KCE的影响，因此推荐使用。的影响，因此推荐使用。3.用等转化率法研究活化能用等转化率法研究活化能E随着反应进度随着反应进度的变化，如基本的变化，如基本保持常数，则为单步（保持常数，则为单步（single-step)简单反应，否则有可能简单反应，否则有可能是多步是多步(multi-step)复杂反应。复杂反应。第51页，共86页，编辑于2022年，星期日3.3.进展篇进展篇v New I
39、deasv New Methodv New Techniques第52页，共86页，编辑于2022年，星期日新概念新概念（new ideas)1.1.只有真正的简单反应才能仅用一套（而且必只有真正的简单反应才能仅用一套（而且必 须是整套）须是整套）“动力学三联体动力学三联体”来表征。来表征。2.2.许多非均相反应其实往往包含了多步反应，故其许多非均相反应其实往往包含了多步反应，故其动力学三联体往往是多步反应共同贡献的表观结动力学三联体往往是多步反应共同贡献的表观结果，并会随着反应进度果，并会随着反应进度或温度或温度T T的变化而变化。的变化而变化。3.3.等温和非等温方法处理的结果不一定相同，
40、因为其温等温和非等温方法处理的结果不一定相同，因为其温度范围不同（前者常常会掩盖反应的复杂本质）；度范围不同（前者常常会掩盖反应的复杂本质）；反之，两者的一致也不能作为获得正确动力学结果反之，两者的一致也不能作为获得正确动力学结果的充要判断条件。的充要判断条件。第53页，共86页，编辑于2022年，星期日新方法新方法（new methods)q 无模式法无模式法（model-free methods)多重加热速率法多重加热速率法 等转化率法等转化率法 等温等转化率法等温等转化率法(popescu method)非参数法非参数法(NPK)l 非非 Arrhenius 温度关系法温度关系法l 其它
41、其它第54页，共86页，编辑于2022年，星期日 等转化率法（等转化率法（isoconversional method)50506060年代：年代：Kissinger-Akahira-Sunose equation(1956,1969)Friedman equation(1964)Flynn-Wall-Ozawa equation(1965)90 90年代年代：Starink equation(1996)Vyazovkin equation(1997,2001)Li-Tang equation(1999)Budrugeac equation(2002)第55页，共86页，编辑于2022年，星期
42、日Starink approach(Thermochim.Acta,288(1996)97)温度积分温度积分 K-A-S法(50 x20)F-W-O法(60 x20)通式通式：（C-取决于取决于 s 的常数）的常数）S=2;B=1 K-A-S 法法 S=0;B=1.0518 F-W-O法法 S=1;B=1 Boswell法法 S=1.8;B=1.0071.2105 E1.004 Starink 法法 第56页，共86页，编辑于2022年，星期日 NL-INT&MNL-INT approach S.Vyazovkin,J.Comput.Chem.18(1997)393,22(2001)178 l
43、NL-INL(non-linear integral)等转化率法：等转化率法：测一组测一组 n 个不同个不同下的下的 TA 曲线，能使下式之值为最曲线，能使下式之值为最 小的小的E 值，即为活化能值值，即为活化能值:MNL-INL(modified non-linear integral)等转化率法：等转化率法：减小积分区间，使减小积分区间，使 E 为常数为常数 第57页，共86页，编辑于2022年，星期日 NL-DIF approach P.Budrugeac,J.Therm.anal.Calor.,68(2002)131NL-DIF(non-linear differential)等转化率
44、法等转化率法 实验方法同前，固定某个实验方法同前，固定某个值，能使下式最小值之值，能使下式最小值之 E为活化能为活化能值：值：第58页，共86页，编辑于2022年，星期日 Chao-Rui Li&Tong B.Tang approach Thermochim.Acta,325(1999)43 微分式：微分式：积分式：积分式：综合式：综合式：式中：式中：（在定（在定处为常数）处为常数）则由一组则由一组测得的测得的TA曲线上取等曲线上取等处数据代入综合式，处数据代入综合式，由该直线由该直线方程之斜率可以求得方程之斜率可以求得E 第59页，共86页，编辑于2022年，星期日 Popescu appr
45、oach(Thermochim.Acta,285(1996)309)“变异的变异的”（Variant)F-W-O法法 定义：定义：第60页，共86页，编辑于2022年，星期日l mn5W T(K)TmTnn4n3n2nn1m1m2m3m4m5Tm1Tm2Tm3Tm4Tm5Tn5Tn4Tn3Tn2Tn1m与与15曲线交点的温度坐标为曲线交点的温度坐标为Tm1 Tm5n与与15曲线交点的温度坐标为曲线交点的温度坐标为Tn1 Tn5Tm与与15曲线交点的质量坐标曲线交点的质量坐标 为为m1 m5Tn与与15曲线交点的质量坐曲线交点的质量坐标为标为n1 n5第61页，共86页，编辑于2022年，星期日
46、1.1.判断机理函数判断机理函数 f()用各种用各种f()进行尝试，计算进行尝试，计算Fmn1 Fmn5，由于，由于Imn 是是常数，故常数，故Fmn1/应为一截矩为零的直线，符合此应为一截矩为零的直线，符合此条件者，为最可机机理函数。条件者，为最可机机理函数。2.2.计算计算 E 和和 A 令令 第62页，共86页，编辑于2022年，星期日 m,n (Tm 1,Tn1)(Tn5,Tn5)作作 1/Tc 直线直线 斜率斜率-E/R；截矩截矩ln A/Fmn 优点：优点：f()与与 k(T)无关无关第63页，共86页，编辑于2022年，星期日 NPK(non-parametric kinetic
47、s)method R.Serra,R.Nomen,J.Sempere,Thermal.Anal.52(1998)933 f(T)不一定是不一定是 k(T)在几个不同在几个不同下测到的下测到的d/d t 的数值组成一个矩阵的数值组成一个矩阵(nm);矩阵的行对应于转化率矩阵的行对应于转化率1 n,列对应于列对应于不同温度不同温度T1 Tm，矩阵，矩阵A的元素为的元素为 于是函数于是函数 f()和和f(T)可以分别表示为可以分别表示为a 和和b 两个列两个列矢量。矢量。第64页，共86页，编辑于2022年，星期日 于是反应速率可以表为一个矩阵积的形式：于是反应速率可以表为一个矩阵积的形式：NPK法
48、应用法应用singular value decomposition(SVD)算法将矩阵算法将矩阵 A 分解为两个矢量分解为两个矢量a 和和 b，这些矢量可进，这些矢量可进一步结合考核实验所得的一步结合考核实验所得的d/d t和和d/d tT关系进关系进一步分析以分别确定模式函数和其它动力学参数。一步分析以分别确定模式函数和其它动力学参数。注注:Refer to J.Sempere et al,Thermochim.Acta,388(2002)407-414 for detail第65页，共86页，编辑于2022年，星期日 Dollimore approach(Thermochim.Acta,2
49、82(1996)13）Arrhenius(Z.Phys.,4(1889)226)Harcourt-Esson(Phil.Trans.R.Soc.London Ser A,212(1913)187)第66页，共86页，编辑于2022年，星期日1.1.判断机理函数判断机理函数 根据以下根据以下TA曲线的特征曲线的特征:a)TG曲线上反应起始和终止温度处的转折情曲线上反应起始和终止温度处的转折情况况(模糊模糊-diffuse,还是清晰还是清晰-sharp)b)DTG曲线峰的低曲线峰的低,高温侧半峰宽之比值高温侧半峰宽之比值c)c)峰温处的峰温处的值值(max)第67页，共86页，编辑于2022年，星
50、期日TWTd/dTlohi组别组别 机理函数机理函数 Ti,Tf 特征特征 lo/hi A A2,A3,A4 Ti sharp,Tf sharp 1 B R2,R3,D1,D2,D3,D4 Ti diffuse,Tf sharp 1 C F1,F2,F3 Ti diffuse,Tf diffuse 1Ti diffuseTf sharp第68页，共86页，编辑于2022年，星期日TG/DTG CurveD2D1R2,D4F1,A2,A3,A4,D3,R3B1F2,F3P1,E1,D1R2D4F1A2,A3,A4D3,R3A2A3A4D3R3max1.0-0.90.9-0.80.8-0.70.7