1、04 热分析动力学热分析动力学热分析动力学概述热分析动力学概述五十年代科学技术的迅速发展特别是航天技术的兴起，迫切需五十年代科学技术的迅速发展特别是航天技术的兴起，迫切需要耐高温的高分子材料。研究高分子材料的热稳定性和使用寿要耐高温的高分子材料。研究高分子材料的热稳定性和使用寿命促进了热重法用于反应动力学的研究。日前，热重法已广泛命促进了热重法用于反应动力学的研究。日前，热重法已广泛用于无机物的脱水、绝食物的热分解、石油高温裂解和煤的热用于无机物的脱水、绝食物的热分解、石油高温裂解和煤的热裂解等的反应动力学研究。裂解等的反应动力学研究。虽然热分析研究反应动力学有许多优点如快速、试样用量少、虽然
2、热分析研究反应动力学有许多优点如快速、试样用量少、不需要分析反应物和产物等，但是由于热分析方法的影响因素不需要分析反应物和产物等，但是由于热分析方法的影响因素多、重复性差和误差较大等缺点，因此在利用热分析法研究反多、重复性差和误差较大等缺点，因此在利用热分析法研究反应动力学时要谨慎，并不是所有反应都适用。应动力学时要谨慎，并不是所有反应都适用。热分析动力学概述热分析动力学概述近几十年来，热重法在测定动力学参数方面，不仅应用领域近几十年来，热重法在测定动力学参数方面，不仅应用领域宽，而且研究的反应类型比较多如热分解反应、脱水反应、结宽，而且研究的反应类型比较多如热分解反应、脱水反应、结品反应等等
3、，由此而积累了大量有关动力学方面的研究资料。品反应等等，由此而积累了大量有关动力学方面的研究资料。在实验方法、数据处理和理论上也有较大的发展，达些都为热在实验方法、数据处理和理论上也有较大的发展，达些都为热重法研究反应动力学打下了牢固的基础。重法研究反应动力学打下了牢固的基础。热分析动力学特点热分析动力学特点1.热分析动力学方法的信息来源是体系变化过程中的物理性质的变热分析动力学方法的信息来源是体系变化过程中的物理性质的变化，因而它对体系所测物理性质以外的其它性质没有任何限制条件，化，因而它对体系所测物理性质以外的其它性质没有任何限制条件，即具有即具有非特异性非特异性 的特点。但这种非特异性是
4、相对的，即热分析方法的特点。但这种非特异性是相对的，即热分析方法只对其测定的物理性质的变化有响应。只对其测定的物理性质的变化有响应。2.现代热分析仪器灵敏度高，热分析动力学方法具有现代热分析仪器灵敏度高，热分析动力学方法具有响应速度快，响应速度快，样品用量少，分析时间短样品用量少，分析时间短 等优点。等优点。3.热分析动力学方法直接检测的是体系的某一物理性质的变化，可热分析动力学方法直接检测的是体系的某一物理性质的变化，可以以同时得到反应过程中相应物理性质变化的静态信息和动态动力学同时得到反应过程中相应物理性质变化的静态信息和动态动力学信息信息。热分析动力学特点热分析动力学特点4.热分析动力学
5、方法可以热分析动力学方法可以 原位、在线、不干扰地连续检测原位、在线、不干扰地连续检测 一个反应，一个反应，从而具有以下优点：从而具有以下优点：（1）可以得到整个过程完整的动力学信息；）可以得到整个过程完整的动力学信息；（2）动力学测量结果比非原位的采样方法更为准确；）动力学测量结果比非原位的采样方法更为准确；（3）测量过程中无需在体系中添加任何试剂，反应后的体系可以很）测量过程中无需在体系中添加任何试剂，反应后的体系可以很方便地进行后续研究与分析；方便地进行后续研究与分析；（4）操作比较简便，不需要在特定的时间点进行采样分析。）操作比较简便，不需要在特定的时间点进行采样分析。热分析动力学特点
6、热分析动力学特点5.热分析方法的影响因素很多，往往重复性较差，实验误差较大，热分析方法的影响因素很多，往往重复性较差，实验误差较大，而且不是所有的化学反应都可以用热分析动力学研究。而且不是所有的化学反应都可以用热分析动力学研究。热分析方法常常用于研究凝聚态特别是固态反应，不同的热分析方热分析方法常常用于研究凝聚态特别是固态反应，不同的热分析方法只适用于相应的反应体系。法只适用于相应的反应体系。例如，气例如，气-气反应不宜用热分析方法，高聚物晶型转变动力学不宜采气反应不宜用热分析方法，高聚物晶型转变动力学不宜采用用TG法进行研究，而法进行研究，而 DTA或或DSC是研究高聚物晶型转变动力学的最是
7、研究高聚物晶型转变动力学的最佳方法。佳方法。热分析动力学的基本原理热分析动力学的基本原理当全自动的热分析仪诞生后，研究者在热分析的动力学研究领域进当全自动的热分析仪诞生后，研究者在热分析的动力学研究领域进行了开创性的工作。行了开创性的工作。在上世纪在上世纪 50年代，年代，Borchardt等提出了最广泛采用的动力学方法，并等提出了最广泛采用的动力学方法，并采用采用DTA技术研究了氯化重氮苯的热分解反应动力学。技术研究了氯化重氮苯的热分解反应动力学。Freeman等采用等采用TG进行了早期的热分解动力学研究。进行了早期的热分解动力学研究。Kissinger提出了一个从提出了一个从 DTA曲线的
8、峰尖温度求算反应活化能的常用曲线的峰尖温度求算反应活化能的常用方法。方法。早期的热分析动力学研究方法是建立在假定反应机理是简单级数反早期的热分析动力学研究方法是建立在假定反应机理是简单级数反应的基础上。然而，许多反应，特别是一些固态反应、高聚物的降应的基础上。然而，许多反应，特别是一些固态反应、高聚物的降解反应，反应机理非常复杂，常常用一个通式解反应，反应机理非常复杂，常常用一个通式f(a)来代表反应机理。来代表反应机理。热分析动力学的基本原理热分析动力学的基本原理热分析动力学是建立在化学热力学、化学动力学及热分析技术基础热分析动力学是建立在化学热力学、化学动力学及热分析技术基础上的一门分支学
9、科。它的基本思想是，用化学动力学的知识，研究上的一门分支学科。它的基本思想是，用化学动力学的知识，研究用热分析方法测定得到的物理量（如质量、温度、热量、模量和尺用热分析方法测定得到的物理量（如质量、温度、热量、模量和尺寸等）的变化速率与温度之间的关系。寸等）的变化速率与温度之间的关系。热分析动力学方法从根本上说是基于这样一个基本原理：在程序控热分析动力学方法从根本上说是基于这样一个基本原理：在程序控制温度下，用物理方法（如制温度下，用物理方法（如 TG法、法、DTA法或法或DSC法等）监测研究体法等）监测研究体系在反应过程中物理性质（如质量、样品与参比物之间的温度差、系在反应过程中物理性质（如
10、质量、样品与参比物之间的温度差、热流差或功率差等）随反应时间或温度的变化，并且监测的物理性热流差或功率差等）随反应时间或温度的变化，并且监测的物理性质的变化正比于反应进度或反应速率。质的变化正比于反应进度或反应速率。非等温法研究动力学过程的优点非等温法研究动力学过程的优点（1）能在反应开始到结束的整个温度范围内连续计算动力学参数；）能在反应开始到结束的整个温度范围内连续计算动力学参数；（2）在等温法过程中，样品必须升到一定温度并有明显的反应才可）在等温法过程中，样品必须升到一定温度并有明显的反应才可测定，很难严格控制反应的始末态，这样的结果往往令人怀疑，而测定，很难严格控制反应的始末态，这样的
11、结果往往令人怀疑，而非等温法无此问题；非等温法无此问题；（3）一条非等温热分析曲线相当于无数条等温热分析曲线，实验样）一条非等温热分析曲线相当于无数条等温热分析曲线，实验样品用量少；品用量少；（4）对于反应进度的分析简单快速，节省时间。）对于反应进度的分析简单快速，节省时间。因此，非等温动力学逐渐成为热分析动力学因此，非等温动力学逐渐成为热分析动力学（Thermal AnalysisKinetics，TAK）的核心。近半个世纪以来在各个方面均有很大发展。）的核心。近半个世纪以来在各个方面均有很大发展。非等温法研究动力学过程的特点非等温法研究动力学过程的特点非等温法研究非均相体系的非等温法研究非
12、均相体系的 TAK过程中，基本上沿用了等温、均相过程中，基本上沿用了等温、均相体系的动力学理论和动力学方程，并作了相应的调整以适应非等温体系的动力学理论和动力学方程，并作了相应的调整以适应非等温非均相体系的需要非均相体系的需要。1.均相体系的浓度（均相体系的浓度（c）的概念在非均相体系中不再适用，用反应转）的概念在非均相体系中不再适用，用反应转化百分率（化百分率（a）来表示非均相体系中的反应进度。考虑到非均相反应）来表示非均相体系中的反应进度。考虑到非均相反应的复杂性，除了均相反应中的简单级数反应动力学方程外，从的复杂性，除了均相反应中的简单级数反应动力学方程外，从20世世纪纪30年代以来建立
13、了许多不同的动力学模型函数年代以来建立了许多不同的动力学模型函数f(a)来描述非均相反来描述非均相反应的动力学过程。应的动力学过程。2.早期的动力学研究工作都是在等温条件下进行的，后来在线性升早期的动力学研究工作都是在等温条件下进行的，后来在线性升温条件下进行动力学研究，通常升温速率为温条件下进行动力学研究，通常升温速率为?，动力学方程作了如下，动力学方程作了如下变形：变形：dt=dT/?。非等温法研究动力学过程的特点非等温法研究动力学过程的特点3.在非等温非均相体系中继续沿用在等温均相反应体系中的动力学在非等温非均相体系中继续沿用在等温均相反应体系中的动力学方程。在绝大多数场合使用的是方程。
14、在绝大多数场合使用的是Arhenius公式来描述反应速率常数公式来描述反应速率常数k(T)与热力学温度与热力学温度T关系：关系：k?Ae?E RT其中其中A为指前因子，为指前因子，E为活化能，为活化能，R为普适气体常数。由此，在升为普适气体常数。由此，在升温速率为温速率为b时，非等温非均相反应的动力学方程就有如下形式：时，非等温非均相反应的动力学方程就有如下形式：d?A?E RT?ef?dT?动力学研究的目的就是求算能描述某反应的动力学研究的目的就是求算能描述某反应的“动力学三因子动力学三因子”(Kinetic Triplet)，即指前因子，即指前因子A、活化能、活化能E和动力学模型函数和动力
15、学模型函数f(a)。微分法微分法在热分析实验过程中，仪器直接记录的信息曲线是在热分析实验过程中，仪器直接记录的信息曲线是a-t的曲线（或的曲线（或a-T的曲线）。热分析仪附带微分单元，或配上计算机进行图形转换处的曲线）。热分析仪附带微分单元，或配上计算机进行图形转换处理，得到理，得到da/dt-T曲线（或曲线（或da/dT-T曲线）采用上式即可进行动力学处曲线）采用上式即可进行动力学处理。由于采用的是理。由于采用的是a对对t（或（或a对对T）一阶微分数据，这种方法常常叫）一阶微分数据，这种方法常常叫微分法，微分法，f(a)又称为微分形式的动力学模型函数。又称为微分形式的动力学模型函数。积分法积
16、分法d?A?E RT?ef?dT?上式进行移项并两端同时积分得到上式进行移项并两端同时积分得到G?A?exp?E?T0TRT?dT?A?exp?E?00TRT?dT?AE?R?P?u?式中，积分下限式中，积分下限 T0的积分值趋近于的积分值趋近于 0，积分下限可由，积分下限可由 0代替。代替。P(u)称称为温度积分（为温度积分（Temperature Integral），其形式如下），其形式如下P?u?e?uu2?du式中式中u=E/RT?由于由于P(u)在数学上得不到有限的精确解，常常由一个近似公式代替。在数学上得不到有限的精确解，常常由一个近似公式代替。直接将直接将a-T数据引入上式，同样
17、可以进行动力学处理。这种数据处数据引入上式，同样可以进行动力学处理。这种数据处理方法常常叫积分法，理方法常常叫积分法，G(a)又称为积分形式的动力学模型函数。又称为积分形式的动力学模型函数。u非等温法研究动力学过程的局限性非等温法研究动力学过程的局限性1.从理论上非等温法的结果与等温法的结果能保持一致。由于等温从理论上非等温法的结果与等温法的结果能保持一致。由于等温反应动力学至少在方法论上比较成熟，其结果的可靠性更高。非等反应动力学至少在方法论上比较成熟，其结果的可靠性更高。非等温法的结果常常拿来与等温法的结果进行比较，来证明非等温法结温法的结果常常拿来与等温法的结果进行比较，来证明非等温法结
18、果的可靠性。大量的事实表明，在很多反应体系中这两种结果很难果的可靠性。大量的事实表明，在很多反应体系中这两种结果很难保持一致。保持一致。2.非均相反应实际上包含多个基元反应平行、连续进行。其转化百非均相反应实际上包含多个基元反应平行、连续进行。其转化百分率是多个基元反应综合的结果，需要对非均相反应的复杂本质进分率是多个基元反应综合的结果，需要对非均相反应的复杂本质进行进一步认识。行进一步认识。非等温法研究动力学过程的局限性非等温法研究动力学过程的局限性3.采用采用Arhenius公式描述热分解反应速率常数与热力学温度公式描述热分解反应速率常数与热力学温度 T关系时，关系时，首先遇到的问题是首先
19、遇到的问题是 Arhenius公式能否适用于非等温非均相体系，寻公式能否适用于非等温非均相体系，寻找更合适的关系式一直是关注的焦点。其次是怎样解释找更合适的关系式一直是关注的焦点。其次是怎样解释Arhenius公公式中两个参数指前因子式中两个参数指前因子 A和活化能和活化能 E的物理含义，求算得到的活化能的物理含义，求算得到的活化能E的数值随转化率发生变化也是一个不容回避的事实。的数值随转化率发生变化也是一个不容回避的事实。非等温法研究动力学过程的局限性非等温法研究动力学过程的局限性4.尽管提出了许多动力学模型函数来描述非均相反应的动力学过程，尽管提出了许多动力学模型函数来描述非均相反应的动力
20、学过程，但是非均相反应本身非常复杂，样品的几何形状的非规整性及反应但是非均相反应本身非常复杂，样品的几何形状的非规整性及反应的理化性质的多变性常常导致实际动力学过程与理想过程推导出来的理化性质的多变性常常导致实际动力学过程与理想过程推导出来的机理不相符合。另外，推导出来的动力学模型函数如此之多，在的机理不相符合。另外，推导出来的动力学模型函数如此之多，在应用这些模型时往往难以入手。如何用尽可能精炼的动力学模型函应用这些模型时往往难以入手。如何用尽可能精炼的动力学模型函数来描述多变的实际动力学过程也不容忽视。数来描述多变的实际动力学过程也不容忽视。等温法等温法为了得到物质有关热现象的动力学数据，
21、样品在指定条件下恒温受为了得到物质有关热现象的动力学数据，样品在指定条件下恒温受热，获得转化百分率热，获得转化百分率a对时间对时间t的曲线，然后根据等温法的动力学方程的曲线，然后根据等温法的动力学方程G?Aexp?E RT?dt?kt0t在在a-t曲线上选取一组带入可能的动力学模型函数曲线上选取一组带入可能的动力学模型函数G(a)式中，如果式中，如果G(a)-t图为一直线，斜率为图为一直线，斜率为k，选取线性相关系数最大的，选取线性相关系数最大的G(a)为最可能为最可能的机理函数。的机理函数。采用同样步骤在不同温度下一系列获得转化百分率采用同样步骤在不同温度下一系列获得转化百分率a对时间对时间
22、t的曲线，的曲线，从而求算得到一组从而求算得到一组k值。值。由式由式 lnk=-E/RT+lnA可知，作可知，作lnA-1/T图得到一条直线，由斜率和图得到一条直线，由斜率和截距可分别得到指前因子截距可分别得到指前因子A和活化能和活化能E的数值。的数值。单升温速率法单升温速率法(非等温法非等温法)单升温速率法是通过在一个升温速率下，对反应测定得到的一条热单升温速率法是通过在一个升温速率下，对反应测定得到的一条热分析曲线上的数据进行动力学分析的方法。将实验得到的分析曲线上的数据进行动力学分析的方法。将实验得到的da/dT-T数数据或据或a-T数据分别引入微分式或积分，尝试将所有可能的动力学模型数
23、据分别引入微分式或积分，尝试将所有可能的动力学模型函数函数f(a)或或G(a)分别带入两式，通过移项两边取对数将方程线性化，分别带入两式，通过移项两边取对数将方程线性化，当温度积分采用当温度积分采用MKN近似公式时，得到近似公式时，得到?ln?1 f?d?dT?ln A?E RT?1.921503?ln?G?/T?ln?AE?R?3.772050?1.921503lnE?1.921503E RT?从回归直线的斜率和截距可以求算动力学参数（指前因子从回归直线的斜率和截距可以求算动力学参数（指前因子A和活化和活化能能E），并根据线性相关性的好坏来判定反应最可能遵循的动力学），并根据线性相关性的好坏
24、来判定反应最可能遵循的动力学模型模型f(a)或或G(a)。单升温速率法单升温速率法(非等温法非等温法)通常根据所选方程是源于微分式还是源于积分式将单升温速率法分通常根据所选方程是源于微分式还是源于积分式将单升温速率法分为微分法和积分法两大类。为微分法和积分法两大类。两类方法各有利弊：两类方法各有利弊：微分法不涉及难解的温度积分，形式简单，但要用到精确的转化率微分法不涉及难解的温度积分，形式简单，但要用到精确的转化率对反应时间或温度的一阶微商数据；对反应时间或温度的一阶微商数据；积分法可以直接用转化率对反应时间或温度的数据，但不能回避温积分法可以直接用转化率对反应时间或温度的数据，但不能回避温度
25、积分问题及由此产生的近似方法的误差。度积分问题及由此产生的近似方法的误差。单升温速率法局限单升温速率法局限在单升温速率法中，由于在单升温速率法中，由于k(T)和和f(a)或或G(a)不能分离，因此在求算动不能分离，因此在求算动力学参数时只能同时得到动力学三因子。力学参数时只能同时得到动力学三因子。这样产生的后果是：良好的线性关系不能保证所选机理模型函数的这样产生的后果是：良好的线性关系不能保证所选机理模型函数的合理性，往往一组实验数据有多个机理模型函数与之相匹配。合理性，往往一组实验数据有多个机理模型函数与之相匹配。在实际操作过程中为了选择合理的机理模型函数，常常采用多种方在实际操作过程中为了
26、选择合理的机理模型函数，常常采用多种方法并用，如非等温法与等温法相结合，微分法与积分法相结合等。法并用，如非等温法与等温法相结合，微分法与积分法相结合等。由于这一方法的科学性正在遭到怀疑，近年来由于这一方法的科学性正在遭到怀疑，近年来ICTAC已不再推荐用已不再推荐用单升温速率法来进行动力学分析。单升温速率法来进行动力学分析。多升温速率法（等转化率法）多升温速率法（等转化率法）多升温速率法是指用不同升温速率所测得的几条热分析曲线来进行多升温速率法是指用不同升温速率所测得的几条热分析曲线来进行动力学分析的一类方法。动力学分析的一类方法。由于大多数多升温速率法常用到几条热分析曲线上同一转化率由于大
27、多数多升温速率法常用到几条热分析曲线上同一转化率a处的处的数据进行动力学处理，所以多升温速率法常常叫做等转化率法。数据进行动力学处理，所以多升温速率法常常叫做等转化率法。这类方法的特点是能将这类方法的特点是能将k(T)和和f(a)或或G(a)分离，在相同转化率分离，在相同转化率a下下f(a)或或G(a)的值不随升温速率的不同发生改变，从而在不引入动力学模型的值不随升温速率的不同发生改变，从而在不引入动力学模型函数的前提条件下得到比较可靠的动力学参数活化能函数的前提条件下得到比较可靠的动力学参数活化能E的数值，因此的数值，因此多升温速率法又称为多升温速率法又称为“Model-free Metho
28、d”。采用多升温速率法得到的活化能采用多升温速率法得到的活化能E的数值常常用来验证单升温速率法的数值常常用来验证单升温速率法结果的可靠性。结果的可靠性。Friedman法法ln?d?dT?ln?Af?E RT在在TA曲线上截取不同升温速率曲线上截取不同升温速率b下相同转化率下相同转化率a时时da/dt-T的值，由的值，由ln(da/dT)b对对1/T作图，用最小二乘法进行线性回归，由斜率可求得作图，用最小二乘法进行线性回归，由斜率可求得在该转化率在该转化率a时活化能时活化能E的数值。的数值。Kissinger法法2?ln?Tp?ln?ARE?f?p?E RTp?当认为当认为f(ap)与与?无关
29、，对于所有的动力学模型函数，其值近似等于无关，对于所有的动力学模型函数，其值近似等于 1，因此在不同升温速率因此在不同升温速率?下由对作图，可得一条直线，由直线斜率和截下由对作图，可得一条直线，由直线斜率和截距可分别求算得到活化能距可分别求算得到活化能E和指前因子和指前因子A的数值。的数值。Flynn-Wall-Ozawa(FWO)法法把温度积分的把温度积分的 Doyle近似式代入，得近似式代入，得ln?ln?AE RG?5.3308?1.0516E RT在在TA曲线上截取不同升温速率曲线上截取不同升温速率?下相同转化率下相同转化率a时时T的值，由的值，由ln?对对1/T作图，用最小二乘法进行
30、线性回归，由斜率可求得在该转化率作图，用最小二乘法进行线性回归，由斜率可求得在该转化率a时活时活化能化能E的数值。的数值。KAS法法把温度积分的把温度积分的 C-R近似式代入，得近似式代入，得2?ln?T?ln AR EG?E RT?在在TA曲线上截取不同升温速率曲线上截取不同升温速率?下相同转化率下相同转化率a时时T的值，由的值，由ln?/T2对对1/T作图，用最小二乘法进行线性回归，由斜率可求得在该转化率作图，用最小二乘法进行线性回归，由斜率可求得在该转化率a时活化能时活化能E的数值。的数值。温度积分近似式温度积分近似式温度积分概念的引入温度积分概念的引入在绝大多数热分析实验过程中，反应体
31、系按照程序线性升高温度。在绝大多数热分析实验过程中，反应体系按照程序线性升高温度。同样，和微分法相比，积分法处理线性升温过程的动力学数据更有同样，和微分法相比，积分法处理线性升温过程的动力学数据更有优势。因此，积分法在现代热分析动力学的研究上得到了广泛的应优势。因此，积分法在现代热分析动力学的研究上得到了广泛的应用。然而，积分法不可避免地面临用。然而，积分法不可避免地面临“温度积分温度积分”这一难题。这一难题。温度积分，又叫温度积分，又叫Arrhenius积分，在热分析动力学发展的历史上起积分，在热分析动力学发展的历史上起着极其重要的作用。着极其重要的作用。积分法克服了微分法的一些缺点，但温度
32、积分又引发新的问题，其积分法克服了微分法的一些缺点，但温度积分又引发新的问题，其来源于不能用一个简单的有限的表达式来精确逼近温度积分。来源于不能用一个简单的有限的表达式来精确逼近温度积分。温度积分概念的引入温度积分概念的引入G?A?exp?E?T0TRT?dT?A?exp?E?00TRT?dT?AE?R?P?u?式中，积分下限式中，积分下限T0的积分值趋近于的积分值趋近于 0，积分下限可由，积分下限可由 0代替。代替。P(u)称称为温度积分（为温度积分（Temperature Integral），其形式如下），其形式如下P?u?e?u?uu?du式中式中u=E/RT2由于由于P(u)在数学上得
33、不到有限的精确解，常常由一个近似公式代替。在数学上得不到有限的精确解，常常由一个近似公式代替。直接将直接将a-T数据引入上式，同样可以进行动力学处理。数据引入上式，同样可以进行动力学处理。推导温度积分近似公式推导温度积分近似公式推导方法有三类：推导方法有三类：（1）级数展开公式；）级数展开公式；（2）复杂的近似公式；）复杂的近似公式；（3）简单的近似公式。）简单的近似公式。温度积分被转化为各种近似的有理函数或有限级数，统称为温度积温度积分被转化为各种近似的有理函数或有限级数，统称为温度积分近似公式。分近似公式。温度积分的分步积分表达式温度积分的分步积分表达式P?u?e?u?u?uu?du?A?
34、2?uu?1 u de?e?e?e?e?e?u22?e?uu2u?e du?u2u?u?2?e?u?uu?e?u?2u?u?2?u?3du?B?u?2u de?3?u?eu?u?2eu?e?u?6?u?4du?C?3?uu?2eu?2eu?2eu?2e2222?uu?6u?4de?uu?6eu?6eu?6e333?u3?u?uu4u?4?6e?udu?4?u?u?u?u?uu?24u?5de?uu?24e4?u?u?u?uu5u?24e du?u?u?5e?u?2u?2!3!4!?1?2?3?L?uuu?D?常用的温度积分近似公式常用的温度积分近似公式1Coats-Redfern近似近似截取温
35、度积分的分步积分表达式（式截取温度积分的分步积分表达式（式D）括号内的第一项和第）括号内的第一项和第二项，而忽略其它项，并带入表达式二项，而忽略其它项，并带入表达式uE/RT，得到温度积分的，得到温度积分的Coats-Redfern近似公式：近似公式：?T02RT?E RT?E RTedT?1?2RT E?eECoats-Redfern 方程是较早推导出来著名近似公式之一，一直被广方程是较早推导出来著名近似公式之一，一直被广泛应用。泛应用。Coats-Redfern 方程通常被进一步简化，通过忽略方程通常被进一步简化，通过忽略2E/RT项而得到所谓的修正的项而得到所谓的修正的 Coats-Re
36、dfern方程：方程：Coats-Redfern方程方程Coats-Redfern 方程是较早推导出来著名近似公式之一，一直被广方程是较早推导出来著名近似公式之一，一直被广泛应用。泛应用。Coats-Redfern 方程通常被进一步简化，通过忽略方程通常被进一步简化，通过忽略2E/RT项而得到所谓的修正的项而得到所谓的修正的 Coats-Redfern方程：方程：?T0e?E RTRT?E RTdT?eE2由于易于实现动力学方程的线性化，这个方程也得到了广泛应用，由于易于实现动力学方程的线性化，这个方程也得到了广泛应用，特别是应用在特别是应用在 KAS等转化率法中。这个方程的精度是相当不够的，
37、等转化率法中。这个方程的精度是相当不够的，因此必须慎重应用。因此必须慎重应用。Doyle方程方程取温度积分的分步积分表达式（式取温度积分的分步积分表达式（式D）括号内的第一项和第二项，）括号内的第一项和第二项，而忽略其它项，并考虑到而忽略其它项，并考虑到 u的取值区间范围的取值区间范围 20u60，经过适当的，经过适当的数学处理及近似，得到温度积分的数学处理及近似，得到温度积分的 Doyle近似式：近似式：P?u?0.00484e?1.0516uDoyle方程也是较早推导出来著名近似公式之一，并且非常容易实方程也是较早推导出来著名近似公式之一，并且非常容易实现动力学方程的线性化，因此得到广泛应
38、用。现动力学方程的线性化，因此得到广泛应用。Doyle方程方程Doyle近似公式在整个近似公式在整个u的取值区间范围内精确度不高。的取值区间范围内精确度不高。Gao等为等为了进一步提高精度，把了进一步提高精度，把u的取值区间范围的取值区间范围 10u70以间隔为以间隔为 u5分解成若干个小的范围，每个小的范围经过各自数学处理及近似，分解成若干个小的范围，每个小的范围经过各自数学处理及近似，得到温度积分在每一小段得到温度积分在每一小段 u的取值区间的的取值区间的 Doyle近似公式，从而大近似公式，从而大大提高大提高Doyle近似公式在整个近似公式在整个u的取值区间范围的精确度。的取值区间范围的
39、精确度。Starink方程方程Starink仔细研究了仔细研究了 Doyle公式的形式，提出用一个通式可以代表公式的形式，提出用一个通式可以代表这样一类近似公式，即这样一类近似公式，即P?u?exp?Au?B?uk式中，式中，A、B和和k为常量。经过在为常量。经过在u的取值区间范围的取值区间范围 20u60计算，计算，指出指出A1时，时，k1.95，B0.235，温度积分近似公式的精度最，温度积分近似公式的精度最高。高。在最近的工作中，在最近的工作中，Starink 认为如果认为如果A的值不等于的值不等于1，将会得到更高，将会得到更高精度的温度积分近似公式，其中精度的温度积分近似公式，其中A1
40、.0008，k1.92，B0.312。Gorbachev-Lee-Beck方程方程式中的方程式中的方程A和和B移项并化简得到移项并化简得到考虑到考虑到u的取值区间范围的取值区间范围 20u60，(12/u)的值趋近于的值趋近于 1，可以认为是一个，可以认为是一个e?u1P?u?2 (A)u 1?2 u常量，因此可将（常量，因此可将（12/u）移出积分）移出积分?ue1?2 u号外，从而得到号外，从而得到Gorbachev-Lee-Beck?(B)22u 1?4 u近似公式近似公式:?e?ue?u?u2?1?2u?du?u2u带入表达式带入表达式uE/RT，得到，得到?T0e?E RTRT?E
41、RTdT?eE?2RT2Li Chung-Hsung方程方程采用完全相同的步骤，式中的方程（采用完全相同的步骤，式中的方程（A）和（）和（C）移项并化简得到）移项并化简得到?e?u2u?ue?1?6 u?du?u2?1?2 u?2?u同样考虑到同样考虑到u的取值区间范围的取值区间范围 20u60，（，（16/u2）的值趋近于）的值趋近于 1，并认为是一个常量，因此可将（并认为是一个常量，因此可将（16/u2）移出积分号外，从而得）移出积分号外，从而得到到Li Chung-Hsung近似公式近似公式:e1?2 uP?u?22u 1?6 u?uAgrawal和冉全印叶素方程和冉全印叶素方程为了进一
42、步提高温度积分近似式的精度，为了进一步提高温度积分近似式的精度，Agrawal，冉全印等经过，冉全印等经过研究研究Gorbachev-Lee-Beck 近似公式和近似公式和Li Chung-Hsung近似式的形式近似式的形式以及它们分别与温度积分数值解的偏差与以及它们分别与温度积分数值解的偏差与u的关系，发现通过调节的关系，发现通过调节式右端分母中式右端分母中u2的系数，可以大大提高温度积分近似式的精度。的系数，可以大大提高温度积分近似式的精度。e1?2 uP?u?22u 1?5 u?ue1?2 uP?u?22u 1?4.6 u?ue1?2 uP?u?22u 1?5.2 u?uSenum-Ya
43、ng方程方程在数学分析中在数学分析中 函数定义为函数定义为?x?e t dt,x?0?t x?10?由由函数定义可知，温度积分可表示为补余不完全函数定义可知，温度积分可表示为补余不完全 函数函数?1,u?P?u?(-1,u)函数可以用连分数表示式进行计算，温度积分可用连分数函数可以用连分数表示式进行计算，温度积分可用连分数表示式展开为表示式展开为eP?u?uu?1?u?u12131?24u?1?LSenum-Yang方程方程对上述连分数的分母截取到第一级、第二级、第三级和第四级，则对上述连分数的分母截取到第一级、第二级、第三级和第四级，则分别得到分别得到Senum-Yang的一级、二级、三级和
44、四级有理近似表达式。的一级、二级、三级和四级有理近似表达式。对连分数的分母截取的级数越高，可以得到更高级数的有理近似表对连分数的分母截取的级数越高，可以得到更高级数的有理近似表达式。常见的温度积分达式。常见的温度积分 Senum-Yang第四级有理近似表达式如下第四级有理近似表达式如下?u32eu?18u?86u?96P?u?432u u?20u?120u?240u?120(-1,u)函数的分母截取的级数越高，温度积分的展开近似表达式函数的分母截取的级数越高，温度积分的展开近似表达式就越复杂，所得到的近似表达式的精确度越大。就越复杂，所得到的近似表达式的精确度越大。Madhusudanan-K
45、rishnan-Ninan方程方程Madhusudanan 等通过对温度积分展开式取二级近似，得到等通过对温度积分展开式取二级近似，得到e1?uP?u?2u 3?u并由该公式出发，采用数学方法与数值计算相结合推导出了温度积并由该公式出发，采用数学方法与数值计算相结合推导出了温度积分的分的MKN近似表达式近似表达式。MKN近似表达式近似表达式的对数形式为的对数形式为?u?lnP(u)?0.297580?1.921503lnu?1.000953u后来又通过对温度积分的展开式取三级近似和系列近似，采用与推后来又通过对温度积分的展开式取三级近似和系列近似，采用与推导上式完全相同的方法推导出来的导上式完
46、全相同的方法推导出来的MKN近似表达式近似表达式 和和。由于。由于MKN近似表达式的对数形式中，变量活化能近似表达式的对数形式中，变量活化能E和温度和温度T通过取对数而通过取对数而分离，非常容易实现热分析动力学方程的线性化，给动力学处理带分离，非常容易实现热分析动力学方程的线性化，给动力学处理带来极大的方便，来极大的方便，评价温度积分评价温度积分P(u)的标准的标准（1）由于发表时间的限制，文献中的陈旧数据不是经过计算机计算由于发表时间的限制，文献中的陈旧数据不是经过计算机计算出来的，这意味着这些数据含有错误。新的出来的，这意味着这些数据含有错误。新的P(u)必须是计算机计算的必须是计算机计算
47、的结果。结果。（2）文献的文献的P(u)数据表包含活化能变量。新的数据表包含活化能变量。新的P(u)数据表必须直接数据表必须直接计算并展示结果。计算并展示结果。（3）不同的文献的不同的文献的P(u)数据表数据表u的取值范围变动很大。新给出的的取值范围变动很大。新给出的P(u)数据表的使用范围最好在数据表的使用范围最好在 0.5u100。（4）Doyle文献中文献中P(u)数据表的有效数字为数据表的有效数字为 4个，个，Zsako 给出的数据给出的数据的有效数字为的有效数字为 5个。新的个。新的P(u)数据表的有效数字必须为数据表的有效数字必须为 8-10个。个。（5）必须在计算机中计算）必须在
48、计算机中计算P(u)数据，在保证计算精度的前提下计算数据，在保证计算精度的前提下计算程序运行必须足够快。程序运行必须足够快。（6）计算的方法必须保证在整个）计算的方法必须保证在整个u的取值范围有足够的精度。的取值范围有足够的精度。温度积分近似表达式温度积分近似表达式的导出的导出在求算动力学参数时，常常采用最小二乘法和迭代法来进行热分析在求算动力学参数时，常常采用最小二乘法和迭代法来进行热分析数据处理。采用形式过于复杂的温度积分近似表达式，必然增加计数据处理。采用形式过于复杂的温度积分近似表达式，必然增加计算的难度，占用大量的浮点计算时间，降低程序运行的效率。算的难度，占用大量的浮点计算时间，降
49、低程序运行的效率。如果考虑到测量误差，也完全没有必要过分追求计算精度。在迭代如果考虑到测量误差，也完全没有必要过分追求计算精度。在迭代计算中，精度越高，迭代越容易出现发散现象。实际上，文献上采计算中，精度越高，迭代越容易出现发散现象。实际上，文献上采用更多的温度积分近似表达式是用更多的温度积分近似表达式是 Coats-Redfern近似公式和近似公式和Doyle近似近似公式。然而这两式的精确度比较低，导致得到的动力学参数的结果公式。然而这两式的精确度比较低，导致得到的动力学参数的结果不可靠。不可靠。温度积分近似表达式温度积分近似表达式的导出的导出目的目的:推导出一种简单可靠（可与推导出一种简单
50、可靠（可与 Coats-Redfern 近似式和近似式和 Doyle近似近似式相比）、精度高（可与一些较复杂的式相比）、精度高（可与一些较复杂的Agrawal近似式等相比）的近似式等相比）的温度积分近似表达式。温度积分近似表达式。另外，在固态反应中，另外，在固态反应中，uE/RT10通常毫无意义，因此也没有刻通常毫无意义，因此也没有刻意追求在非常低的意追求在非常低的u的取值区间赋予新的温度积分近似表达式以特的取值区间赋予新的温度积分近似表达式以特别高的精度。别高的精度。温度积分近似表达式温度积分近似表达式的导出的导出通过考察通过考察 Gorbachev-Lee-Beck 方程的推导过程，发现在
