[go: up one dir, main page]

CN112233729A - General form and solving method of polymer insulation pyrolysis reaction mechanism function - Google Patents

General form and solving method of polymer insulation pyrolysis reaction mechanism function Download PDF

Info

Publication number
CN112233729A
CN112233729A CN202011067826.1A CN202011067826A CN112233729A CN 112233729 A CN112233729 A CN 112233729A CN 202011067826 A CN202011067826 A CN 202011067826A CN 112233729 A CN112233729 A CN 112233729A
Authority
CN
China
Prior art keywords
reaction mechanism
function
formula
polymer insulation
pyrolysis reaction
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011067826.1A
Other languages
Chinese (zh)
Inventor
李庆民
刘红磊
任鹏
彭鹏
刘秋实
金虎
李锐海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202011067826.1A priority Critical patent/CN112233729A/en
Publication of CN112233729A publication Critical patent/CN112233729A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/20Identification of molecular entities, parts thereof or of chemical compositions
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16CCOMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/04Ageing analysis or optimisation against ageing
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation
    • GPHYSICS
    • G06COMPUTING OR CALCULATING; COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Computing Systems (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • General Health & Medical Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Addition Polymer Or Copolymer, Post-Treatments, Or Chemical Modifications (AREA)

Abstract

本发明涉及一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法。包括如下步骤:在反应机理函数微分表达式的基础上,构造一种4参数的反应机理函数积分表达式模型;基于聚合物绝缘材料在非等温条件下的动力学方程,给出实验函数的构造形式;基于实验函数求解聚合物绝缘热解反应机理函数积分表达式。本发明可表征各种已有的反应机理函数,具有更好的普适性和应用范围,可有效反映聚合物绝缘的老化与劣化过程,为聚合物绝缘材料的老化机理研究和剩余寿命预测提供基础方法。

Figure 202011067826

The invention relates to a general characterization form and a solution method of a polymer insulation pyrolysis reaction mechanism function. It includes the following steps: on the basis of the differential expression of the reaction mechanism function, a 4-parameter integral expression model of the reaction mechanism function is constructed; based on the kinetic equation of the polymer insulating material under non-isothermal conditions, the construction of the experimental function is given Form; based on the experimental function to solve the integral expression of the polymer insulation pyrolysis reaction mechanism function. The invention can characterize various existing reaction mechanism functions, has better universality and application range, can effectively reflect the aging and deterioration process of polymer insulation, and provides research on the aging mechanism of polymer insulation materials and remaining life prediction. basic method.

Figure 202011067826

Description

一种聚合物绝缘热解反应机理函数的通用形式与求解方法A general form and solution method of a polymer insulation pyrolysis reaction mechanism function

技术领域technical field

本发明涉及绝缘材料老化评估领域,尤其涉及一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法。The invention relates to the field of insulation material aging assessment, in particular to a general characterization form and a solution method of a polymer insulation pyrolysis reaction mechanism function.

背景技术Background technique

聚合物绝缘材料因其优异的绝缘性能、耐热性能和易于加工的特点,被广泛应用于电力设备的绝缘系统中。然而,由于电、磁、热等多种应力作用,聚合物绝缘材料会出现老化与劣化现象,严重时可引发安全事故。许多电力设备投运至今已有十多年甚至几十年的时间,亟需制定科学的状态维修策略。因此,对聚合物绝缘材料的老化机理和寿命预测方法研究具有重要的工程应用价值。Polymer insulating materials are widely used in electrical equipment insulation systems due to their excellent insulating properties, heat resistance and easy processing. However, due to the effects of electrical, magnetic, thermal and other stresses, the polymer insulating material will experience aging and deterioration, which may lead to safety accidents in severe cases. Many power equipment has been put into operation for more than ten years or even decades, and it is urgent to formulate a scientific condition maintenance strategy. Therefore, the research on the aging mechanism and life prediction method of polymer insulating materials has important engineering application value.

材料在老化与劣化过程中,必定伴随着化学变化,存在物质的产出或溢出,在物质内部产生不同的形变和内应力。这些变化表现在固相反应中,会在反应物界面形成一个活性区域或活性点,这就是晶核。晶核在多种应力作用下,不断地生长、长大、扩散,最终导致材料失效。反应机理函数作为对晶核生长行为的数学描述,被用来表征石油与煤的高温裂解、高聚物的聚合固化和无机物的脱水分解等诸多演化过程,揭示其化学反应的发生路径和分解特性,以确定材料的稳定性和寿命。In the process of aging and deterioration, materials must be accompanied by chemical changes, and there will be the output or overflow of substances, resulting in different deformations and internal stresses inside the substances. These changes are manifested in solid-phase reactions, and an active region or active site is formed at the reactant interface, which is the nuclei. Under the action of various stresses, the crystal nucleus continuously grows, grows, and diffuses, eventually leading to material failure. As a mathematical description of the growth behavior of crystal nuclei, the reaction mechanism function is used to characterize many evolutionary processes such as high-temperature cracking of petroleum and coal, polymerization and solidification of polymers, and dehydration and decomposition of inorganic substances, revealing the occurrence path and decomposition of chemical reactions. properties to determine material stability and longevity.

反应机理函数研究在晶核概念提出后得到了迅速发展,在过去的近一个世纪中,许多热动力学研究者提出并完善了不同类型的固态动力学模型。尽管这些动力学模式函数对许多固态物质的反应过程给出了基本描述,但由于非均相反应本身的复杂性,加上实际样品颗粒几何形状的非规整性和堆积的非规则性,以及反应物质物理化学性质的多变性,常常会出现实际的热分解曲线与预估机理不相符合的情况,求解的反应机理函数必然失真。长期以来,因缺少一个较通用的反应机理函数形式,造成材料热动力分析的差异性和复杂性。原则上,可通过增加动力学模型函数的个数来克服这一局限,这些模型虽可解释某些材料反应机理偏离理想模型的情况,但缺少实际物理意义,且模型数量也不能无限增加。另一种解决方案是寻求一个更加普适的数学模型,使其涵盖尽可能多的反应机理,具有更好的推广性。The study of reaction mechanism function has been developed rapidly after the concept of crystal nucleus was proposed. In the past century, many thermodynamic researchers have proposed and perfected different types of solid-state kinetic models. Although these kinetic mode functions give a basic description of the reaction process of many solid-state substances, due to the complexity of the heterogeneous reaction itself, coupled with the irregularity of the particle geometry and packing of the actual sample, and the reaction Due to the variability of physical and chemical properties of substances, the actual thermal decomposition curve often does not conform to the predicted mechanism, and the solved reaction mechanism function must be distorted. For a long time, the lack of a more general functional form of reaction mechanism has resulted in the difference and complexity of material thermodynamic analysis. In principle, this limitation can be overcome by increasing the number of kinetic model functions. Although these models can explain the deviation of some material reaction mechanisms from the ideal model, they lack practical physical meaning, and the number of models cannot be increased indefinitely. Another solution is to seek a more general mathematical model that covers as many reaction mechanisms as possible and has better generalizability.

发明内容SUMMARY OF THE INVENTION

针对现有技术的不足,本发明提供一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法,包括以下步骤。Aiming at the deficiencies of the prior art, the present invention provides a general characterization form and a solution method of a polymer insulation pyrolysis reaction mechanism function, including the following steps.

1. 一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法,包括以下步骤:1. A general characterization form and solution method of a polymer insulation pyrolysis reaction mechanism function, comprising the following steps:

(1)在反应机理函数微分表达式的基础上,构造一种4参数的反应机理函数积分表达式模型,具体表达式为:(1) On the basis of the differential expression of the reaction mechanism function, a 4-parameter integral expression model of the reaction mechanism function is constructed. The specific expression is:

Figure DEST_PATH_IMAGE001
(1-1)
Figure DEST_PATH_IMAGE001
(1-1)

对于具体的材料热解反应,式中qmnp为常数;For a specific material pyrolysis reaction, q , m , n , and p are constants in the formula;

(2)基于聚合物绝缘材料在非等温条件下的动力学方程,给出实验函数的构造形式;(2) Based on the kinetic equation of the polymer insulating material under non-isothermal conditions, the structural form of the experimental function is given;

(3)基于实验函数求解聚合物绝缘热解反应机理函数积分表达式。(3) Based on the experimental function, the integral expression of the polymer insulation pyrolysis reaction mechanism function is solved.

2. 可选的,反应机理函数的微分表达式为:2. Optionally, the differential expression of the reaction mechanism function is:

Figure DEST_PATH_IMAGE002
(1-2)
Figure DEST_PATH_IMAGE002
(1-2)

其中,mnp代表不同数值的指数因子,在应用时其中的一个指数因子总是为0。Among them, m , n , p represent exponential factors of different values, and one of the exponential factors is always 0 when applied.

3. 可选的,实验函数的构造具体包括以下步骤:3. Optionally, the construction of the experimental function specifically includes the following steps:

(1)聚合物绝缘材料在非等温条件下的动力学方程为:(1) The kinetic equation of the polymer insulating material under non-isothermal conditions is:

Figure DEST_PATH_IMAGE003
(1-3)
Figure DEST_PATH_IMAGE003
(1-3)

(2)对式(1-3)移项并且两端同时积分可得:(2) Shift the term of equation (1-3) and integrate both ends simultaneously to obtain:

Figure DEST_PATH_IMAGE004
(1-4)
Figure DEST_PATH_IMAGE004
(1-4)

其中u=E/RTP(u)称为温度积分,如下式(1-5)所示:Where u=E/RT , P(u) is called the temperature integral, as shown in the following formula (1-5):

Figure DEST_PATH_IMAGE005
(1-5)
Figure DEST_PATH_IMAGE005
(1-5)

(3)温度积分P(u)用Senum-Yang温度积分近似表达,如下式(1-6)所示:(3) The temperature integral P(u) is approximately expressed by the Senum-Yang temperature integral, as shown in the following formula (1-6):

Figure DEST_PATH_IMAGE006
(1-6)
Figure DEST_PATH_IMAGE006
(1-6)

(4)以α=0.5为参考点,由式(1-4)可得:(4) Taking α = 0.5 as the reference point, it can be obtained from formula (1-4):

Figure DEST_PATH_IMAGE007
(1-7)
Figure DEST_PATH_IMAGE007
(1-7)

式中u 0.5α=0.5时,所对应的u值;where u 0.5 is the corresponding u value when α = 0.5;

(5)式(1-4)除以式(1-7)可得:(5) Equation (1-4) divided by Equation (1-7) can be obtained:

Figure DEST_PATH_IMAGE008
(1-8)
Figure DEST_PATH_IMAGE008
(1-8)

(6)将式(1-1)和(1-6)代入(1-8)可得:(6) Substitute equations (1-1) and (1-6) into (1-8) to obtain:

Figure DEST_PATH_IMAGE009
(1-9)
Figure DEST_PATH_IMAGE009
(1-9)

(7)令:(7) Order:

Figure DEST_PATH_IMAGE010
(1-10)
Figure DEST_PATH_IMAGE010
(1-10)

(8)对聚合物绝缘材料进行热失重分析,根据热失重曲线求取不同转化率α下材料的活化能;(8) The thermal weight loss analysis of the polymer insulating material is carried out, and the activation energy of the material at different conversion rates α is obtained according to the thermal weight loss curve;

(9)将在聚合物绝缘材料转化率α范围内得到的活化能E和温度T代入式(1-10),构成实验函数曲线,如果热解反应可用单一的反应机理函数描述,则不论升温速率是多少,各条实验曲线应互相重合。(9) Substitute the activation energy E and temperature T obtained in the range of the conversion rate α of the polymer insulating material into the formula (1-10) to form the experimental function curve. If the pyrolysis reaction can be described by a single reaction mechanism function, regardless of the temperature rise What is the rate, and each experimental curve should coincide with each other.

4. 可选的,聚合物绝缘热解反应机理函数积分表达式的求解具体包括以下步骤:4. Optionally, the solution of the integral expression of the polymer insulation pyrolysis reaction mechanism function specifically includes the following steps:

(1)令:(1) Order:

Figure DEST_PATH_IMAGE011
(1-11)
Figure DEST_PATH_IMAGE011
(1-11)

(2)将不同转化率α代入式(1-11),由最小二乘法拟合式(1-10)得到的实验曲线即可求得所对应参数mnp(2) Substitute different conversion rates α into formula (1-11), and the corresponding parameters m , n , and p can be obtained from the experimental curve obtained by fitting formula (1-10) by the least square method;

(3)将式(1-1)、式(1-6)代入式(1-4)可得:(3) Substitute equations (1-1) and (1-6) into equations (1-4) to obtain:

Figure DEST_PATH_IMAGE012
(1-12)
Figure DEST_PATH_IMAGE012
(1-12)

(4)式(1-12)中参数mnpu、E、βR已知,未知量只有qA值,任取两个转化率α所对应的参数代入式(1-12),求解方程组即可得q值;(4) In formula (1-12), the parameters m , n , p , u, E, β , and R are known, and the unknowns are only the values of q and A , and the parameters corresponding to the two conversion rates α are arbitrarily selected and substituted into formula (1) -12), the q value can be obtained by solving the equation system;

(5)将求得的参数qmnp代入式(1-1)即得所求材料热解反应过程所遵循的反应机理函数。(5) Substitute the obtained parameters q , m , n , and p into formula (1-1) to obtain the reaction mechanism function followed by the pyrolysis reaction process of the obtained material.

本发明提供一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法,该模型可以表征各种已有的反应机理函数,具有更好的普适性和应用范围,可更好地描述一些多元反应,有效表征绝缘材料的老化与劣化过程,为绝缘材料的老化机理研究和剩余寿命预测提供基础方法。The present invention provides a general characterization form and a solution method of a polymer insulation pyrolysis reaction mechanism function, the model can characterize various existing reaction mechanism functions, has better universality and application range, and can better describe Some multivariate reactions can effectively characterize the aging and deterioration process of insulating materials, and provide a basic method for studying the aging mechanism of insulating materials and predicting the remaining life.

附图说明Description of drawings

图1为本发明的实验流程示意图。FIG. 1 is a schematic diagram of the experimental flow of the present invention.

图2为交联聚乙烯的热失重曲线。Figure 2 is the thermogravimetric curve of cross-linked polyethylene.

图3为改进Flynn-Wall-Ozawa法拟合直线。Figure 3 is a straight line fitted by the improved Flynn-Wall-Ozawa method.

图4为交联聚乙烯不同升温速率下的实验曲线。Figure 4 shows the experimental curves of cross-linked polyethylene at different heating rates.

具体实施方式Detailed ways

下面结合具体实施例和附图对本发明作进一步说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。The present invention will be further described below with reference to specific embodiments and accompanying drawings. It should be understood that the specific embodiments described herein are only used to explain the present invention, but not to limit the present invention.

如图1所示,一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法,以交联聚乙烯为例,主要包括以下步骤。As shown in Figure 1, a general characterization form and solution method of a polymer insulation pyrolysis reaction mechanism function, taking cross-linked polyethylene as an example, mainly includes the following steps.

(1)反应机理函数的微分表达式为:(1) The differential expression of the reaction mechanism function is:

Figure DEST_PATH_IMAGE013
(1-1)
Figure DEST_PATH_IMAGE013
(1-1)

其中,mnp代表不同数值的指数因子,在应用时其中的一个指数因子总是为0;Among them, m , n , p represent exponential factors of different values, and one of the exponential factors is always 0 when applied;

一般情况下,材料剩余寿命计算和反应机理函数的积分表达形式G(α)密切相关,因此研究材料热解反应遵循的积分式表达模型具有重要意义。在反应机理函数微分表达式的基础上,构造一种4参数的反应机理函数积分表达式模型,具体表达式为:In general, the calculation of the residual life of a material is closely related to the integral expression form G ( α ) of the reaction mechanism function, so it is of great significance to study the integral expression model followed by the pyrolysis reaction of the material. On the basis of the differential expression of the reaction mechanism function, a 4-parameter integral expression model of the reaction mechanism function is constructed. The specific expression is:

Figure 596463DEST_PATH_IMAGE001
(1-2)
Figure 596463DEST_PATH_IMAGE001
(1-2)

对于具体的材料热解反应,式中qmnp为常数。For a specific material pyrolysis reaction, q , m , n , and p are constants in the formula.

(2)基于非等温条件下的动力学方程,给出实验函数的构造形式:(2) Based on the kinetic equation under non-isothermal conditions, the structural form of the experimental function is given:

1)非等温条件下的动力学方程为:1) The kinetic equation under non-isothermal conditions is:

Figure 368110DEST_PATH_IMAGE003
(1-3)
Figure 368110DEST_PATH_IMAGE003
(1-3)

2)对式(1-3)移项并且两端同时积分可得:2) Shift the term of Equation (1-3) and integrate both ends at the same time to obtain:

Figure 621543DEST_PATH_IMAGE004
(1-4)
Figure 621543DEST_PATH_IMAGE004
(1-4)

其中u=E/RTP(u)称为温度积分,如下式(1-5)所示:Where u=E/RT , P(u) is called the temperature integral, as shown in the following formula (1-5):

Figure 504049DEST_PATH_IMAGE005
(1-5)
Figure 504049DEST_PATH_IMAGE005
(1-5)

3)温度积分P(u)用Senum-Yang温度积分近似表达,如下式(1-6)所示:3) The temperature integral P(u) is approximately expressed by the Senum-Yang temperature integral, as shown in the following formula (1-6):

Figure 856533DEST_PATH_IMAGE006
(1-6)
Figure 856533DEST_PATH_IMAGE006
(1-6)

4)以α=0.5为参考点,由式(1-4)可得:4) Taking α = 0.5 as the reference point, it can be obtained from formula (1-4):

Figure 861398DEST_PATH_IMAGE007
(1-7)
Figure 861398DEST_PATH_IMAGE007
(1-7)

式中u 0.5α=0.5时,所对应的u值;where u 0.5 is the corresponding u value when α = 0.5;

5)式(1-4)除以式(1-7)可得:5) Equation (1-4) divided by Equation (1-7) can be obtained:

Figure 112382DEST_PATH_IMAGE008
(1-8)
Figure 112382DEST_PATH_IMAGE008
(1-8)

6)将式(1-2)和(1-6)代入(1-8)可得:6) Substitute equations (1-2) and (1-6) into (1-8) to get:

Figure 782266DEST_PATH_IMAGE009
(1-9)
Figure 782266DEST_PATH_IMAGE009
(1-9)

7)令:7) Order:

Figure 254836DEST_PATH_IMAGE010
(1-10)
Figure 254836DEST_PATH_IMAGE010
(1-10)

8)采用瑞士梅特勒公司的TGA-DSC3+同步热分析联用仪,对电缆交联聚乙烯样本进行热重测试,每次取样10mg,温度上限为600℃,在50ml/min的氮气条件下进行实验,选用氧化铝坩埚,升温速率分别设置为5、10、15、20和25°C/min。测得交联聚乙烯的热失重曲线(TG)如图2所示;8) Using the TGA-DSC3+ synchronous thermal analysis analyzer from Mettler, Switzerland, the cable cross-linked polyethylene sample was tested by thermogravimetry, 10 mg per sample, the upper temperature limit was 600 °C, and under the condition of 50 ml/min nitrogen Carry out the experiment, select the alumina crucible for use, and the heating rate is set to 5, 10, 15, 20 and 25 ℃/min respectively. The thermogravimetric curve (TG) of the measured cross-linked polyethylene is shown in Figure 2;

Flynn-Wall-Ozawa法不需设定反应机理函数就可以求解活化能,传统的Flynn-Wall-Ozawa法选择Doyle近似方法进行表征,具有较大的计算误差。本发明选择改进的Flynn-Wall-Ozawa法,其具体表达形式为 The Flynn-Wall-Ozawa method can solve the activation energy without setting the reaction mechanism function. The traditional Flynn-Wall-Ozawa method selects the Doyle approximation method for characterization, which has a large calculation error. The present invention selects the improved Flynn-Wall-Ozawa method, and its concrete expression form is

Figure 100002_DEST_PATH_IMAGE014
(1-11)
Figure 100002_DEST_PATH_IMAGE014
(1-11)

针对不同升温速率下的β选择相同的转化率α,对于同一热解反应G(α)为一恒定值,则ln(β/T 1.89466100)与1/T满足线性关系,通过拟合斜率可求得活化能EThe same conversion rate α is selected for β under different heating rates, and G ( α ) is a constant value for the same pyrolysis reaction, then ln( β / T 1.89466100 ) and 1 / T satisfy a linear relationship, which can be obtained by fitting the slope. get activation energy E ;

不同转化率α下对应的温度值如表1所示,针对不同转化率α的拟合曲线如图3所示,由拟合曲线斜率计算求得交联聚乙烯在不同转化率α下的活化能,如表2所示;The corresponding temperature values at different conversion rates α are shown in Table 1, and the fitting curves for different conversion rates α are shown in Figure 3. The activation of cross-linked polyethylene at different conversion rates α is calculated from the slope of the fitting curve. can, as shown in Table 2;

表1 不同转化率对应的温度值 Table 1 Temperature values corresponding to different conversion ratios

<i>α</i><i>α</i> 5°C/min5°C/min 10°C/min10°C/min 15°C/min15°C/min 20°C/min20°C/min 25°C/min25°C/min 0.20.2 451.381451.381 461.798461.798 469.017469.017 474.94474.94 480.326480.326 0.30.3 457.378457.378 467.909467.909 475.196475.196 480.912480.912 486.474486.474 0.40.4 461.761461.761 472.349472.349 479.644479.644 485.553485.553 490.879490.879 0.50.5 465.219465.219 475.977475.977 483.347483.347 489.393489.393 494.685494.685 0.60.6 468.373468.373 479.288479.288 486.811486.811 493.183493.183 497.583497.583 0.70.7 471.267471.267 482.273482.273 489.836489.836 496.112496.112 501.012501.012 0.80.8 474.353474.353 485.591485.591 493.275493.275 499.345499.345 504.619504.619

表2 不同转化率下的活化能Table 2 Activation energies at different conversion rates

转化率<i>α</i>(%)Conversion rate <i>α</i>(%) 0.20.2 0.30.3 0.40.4 0.50.5 0.60.6 0.70.7 0.80.8 活化能(kJ/mol)Activation energy (kJ/mol) 241.43241.43 244.84244.84 246.9246.9 246.05246.05 247.68247.68 246.47246.47 245.33245.33

9)将在交联聚乙烯转化率α范围内得到的活化能E和温度T代入式(1-10),构成实验函数曲线,图4为不同升温速率下拟合得到的实验函数曲线,可以看出,5种升温速率下的实验曲线几乎全部重叠,表明交联聚乙烯热解反应遵循单一的反应机理函数。9) Substitute the activation energy E and temperature T obtained in the range of the cross-linked polyethylene conversion rate α into formula (1-10) to form the experimental function curve. Figure 4 shows the experimental function curve obtained by fitting under different heating rates, which can be It can be seen that the experimental curves under the five heating rates almost all overlap, indicating that the cross-linked polyethylene pyrolysis reaction follows a single reaction mechanism function.

(3)基于交联聚乙烯的实验函数求解热解反应机理函数积分表达式:(3) Based on the experimental function of cross-linked polyethylene, the integral expression of the pyrolysis reaction mechanism function is solved:

1)令: 1) Order:

Figure DEST_PATH_IMAGE015
(1-12)
Figure DEST_PATH_IMAGE015
(1-12)

2)将不同转化率α代入式(1-12),由拟合的实验曲线可求得m=0.6689、n=-0.1420、p=0.0902;2) Substitute different conversion rates α into formula (1-12), and from the fitted experimental curve, m = 0.6689, n = -0.1420, p = 0.0902;

3)将式(1-2)、式(1-6)代入式(1-4)可得:3) Substitute formula (1-2) and formula (1-6) into formula (1-4) to obtain:

Figure 745117DEST_PATH_IMAGE012
(1-13)
Figure 745117DEST_PATH_IMAGE012
(1-13)

4)式(1-13)中参数mnpu、E、βR已知,未知量只有qA值,任取两个转化率α=0.2和α =0.3所对应的参数代入式(1-13),求解方程组即可得q=2.1150;4) In formula (1-13), the parameters m , n , p , u, E, β , and R are known, and the unknowns are only the values of q and A , and any two conversion rates corresponding to α = 0.2 and α = 0.3 are taken. Substitute the parameters into equation (1-13), and solve the equation system to get q = 2.1150;

5)将求得的参数qmnp代入式(1-2)可得交联聚乙烯热解反应过程所遵循的反应机理函数为:5) Substituting the obtained parameters q , m , n , and p into formula (1-2), the reaction mechanism function of the cross-linked polyethylene pyrolysis reaction process can be obtained as follows:

Figure 100002_DEST_PATH_IMAGE016
(1-14)
Figure 100002_DEST_PATH_IMAGE016
(1-14)

综上所述,本发明提供一种聚合物绝缘热解反应机理函数的通用表征形式与求解方法,该模型可以表征各种已有的反应机理函数,具有更好的普适性和应用范围,为绝缘材料的老化机理研究和剩余寿命预测提供基础方法。To sum up, the present invention provides a general characterization form and solution method for the reaction mechanism function of polymer insulation pyrolysis, the model can characterize various existing reaction mechanism functions, and has better universality and application range, It provides a basic method for studying the aging mechanism of insulating materials and predicting the remaining life.

上述实施例为本发明较佳的实施方式,但本发明的实施方式并不受上述实施例的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。The above-mentioned embodiments are preferred embodiments of the present invention, but the embodiments of the present invention are not limited by the above-mentioned embodiments, and any other changes, modifications, substitutions, combinations, The simplification should be equivalent replacement manners, which are all included in the protection scope of the present invention.

Claims (4)

1. A general representation form and solving method of a polymer insulation pyrolysis reaction mechanism function is characterized by comprising the following steps:
(1) on the basis of a differential expression of a reaction mechanism function, a 4-parameter reaction mechanism function integral expression model is constructed, and the specific expression is as follows:
Figure DEST_PATH_IMAGE002A
(1-1)
for a particular material pyrolysis reaction, in the formulaqmnpIs a constant;
(2) based on a kinetic equation of the polymer insulating material under a non-isothermal condition, giving a construction form of an experimental function;
(3) and solving a function integral expression of the polymer insulation pyrolysis reaction mechanism based on the experimental function.
2. The general characterization form and solution method for polymer insulation pyrolysis reaction mechanism function according to claim 1, wherein the differential expression of the reaction mechanism function is:
Figure DEST_PATH_IMAGE004A
(1-2)
wherein,mnpindex factors representing different values, one of which is always 0 when applied.
3. The general characterization form and solution method for polymer insulation pyrolysis reaction mechanism function according to claim 1, wherein the configuration of the experimental function specifically comprises the following steps:
(1) the kinetic equation of the polymer insulation material under non-isothermal conditions is as follows:
Figure DEST_PATH_IMAGE006A
(1-3)
(2) shifting the term of equation (1-3) and integrating both ends simultaneously can yield:
Figure DEST_PATH_IMAGE008A
(1-4)
whereinu=E/RTP(u)Referred to as temperature integral, as shown in the following formulas (1-5):
Figure DEST_PATH_IMAGE010A
(1-5)
(3) integral of temperatureP(u)Expressed approximately by Senum-Yang temperature integration, as shown in the following formulas (1-6):
Figure DEST_PATH_IMAGE012A
(1-6)
(4) to be provided withα=0.5 as a reference point, which can be obtained from the formula (1-4):
Figure DEST_PATH_IMAGE014
(1-7)
in the formulau 0.5Is composed ofαWhen =0.5, the correspondinguA value;
(5) the formula (1-4) is divided by the formula (1-7):
Figure DEST_PATH_IMAGE016
(1-8)
(6) substituting formulae (1-1) and (1-6) into (1-8) can give:
Figure DEST_PATH_IMAGE018
(1-9)
(7) order:
Figure DEST_PATH_IMAGE020
(1-10)
(8) carrying out thermal weight loss analysis on the polymer insulating material, and obtaining different conversion rates according to a thermal weight loss curveαThe activation energy of the lower material;
(9) will be in the conversion of polymer insulation materialαActivation energy obtained within the rangeEAnd temperatureTThe equations (1-10) are substituted to form experimental function curves, and if the pyrolysis reaction can be described by a single reaction mechanism function, the experimental curves should coincide with each other regardless of the temperature rise rate.
4. The general characterization form and solution method for polymer insulation pyrolysis reaction mechanism function according to claim 1, wherein the solution of the integral expression of the polymer insulation pyrolysis reaction mechanism function specifically comprises the following steps:
(1) order:
Figure DEST_PATH_IMAGE022
(1-11)
(2) will vary in conversion rateαSubstituting the formula (1-11), and obtaining the corresponding parameters by fitting the formula (1-10) with the least square methodmnp
(3) The compound represented by the formula (1-1) or (1-6) may be substituted with the compound represented by the formula (1-4):
Figure DEST_PATH_IMAGE024
(1-12)
(4) parameters of the formulae (1-12)mnpu、E、βRThe unknown quantity is known onlyqAndAvalue, arbitrarily two conversion ratesαThe corresponding parameters are substituted into the formula (1-12), and the equation set is solved to obtainqA value;
(5) the parameters to be obtainedqmnpSubstituting the formula (1-1) to obtain a reaction mechanism function followed by the pyrolysis reaction process of the required material.
CN202011067826.1A 2020-10-07 2020-10-07 General form and solving method of polymer insulation pyrolysis reaction mechanism function Pending CN112233729A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011067826.1A CN112233729A (en) 2020-10-07 2020-10-07 General form and solving method of polymer insulation pyrolysis reaction mechanism function

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011067826.1A CN112233729A (en) 2020-10-07 2020-10-07 General form and solving method of polymer insulation pyrolysis reaction mechanism function

Publications (1)

Publication Number Publication Date
CN112233729A true CN112233729A (en) 2021-01-15

Family

ID=74120462

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011067826.1A Pending CN112233729A (en) 2020-10-07 2020-10-07 General form and solving method of polymer insulation pyrolysis reaction mechanism function

Country Status (1)

Country Link
CN (1) CN112233729A (en)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2798036A1 (en) * 2012-12-05 2014-06-05 Nova Chemicals Corporation Reduction of fouling in high pressure reactors
US20140284198A1 (en) * 2011-10-10 2014-09-25 Lepez Conseils Finance Innovations-Lcfi Process and installation for pyrolysis of a product in the form of divided solids, in particular polymer waste
CN108918985A (en) * 2018-06-29 2018-11-30 华北电力大学 A kind of improved method for restoring distribution of space charge based on pulse electroacoustic method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140284198A1 (en) * 2011-10-10 2014-09-25 Lepez Conseils Finance Innovations-Lcfi Process and installation for pyrolysis of a product in the form of divided solids, in particular polymer waste
CA2798036A1 (en) * 2012-12-05 2014-06-05 Nova Chemicals Corporation Reduction of fouling in high pressure reactors
CN108918985A (en) * 2018-06-29 2018-11-30 华北电力大学 A kind of improved method for restoring distribution of space charge based on pulse electroacoustic method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
REN, P等: "Thermal Decomposition Mechanism of GIS Basin Insulator and Kinetic Parameters-Based Lifetime Prediction Methodology", 《POLYMERS》, pages 653 *
UGUR, M等: "A study on the ageing process for polyester resin using improved Weibull statistics", 《ELECTRICAL ENGINEERING》, pages 283 - 288 *
刘红磊: "变压器油纸绝缘老化的动力学特性与寿命评估方法", 《中国优秀硕士学位论文全文数据库》, pages 042 - 1294 *
李登淑等: "基于集对分析动态加权法的XLPE绝缘老化状态评估", 《水电能源科学 》, pages 194 - 197 *
林智勇等: "回复电压微分谱线特性的变压器绝缘老化研究", 《仪器仪表学报》, pages 1954 - 1960 *

Similar Documents

Publication Publication Date Title
Ye et al. Temperature control of battery modules through composite phase change materials with dual operating temperature regions
Ran et al. Spiral‐structured dielectric polymers exhibiting ultrahigh energy density and charge–discharge efficiency at high temperatures
CN104793111B (en) Based on reason, change, the insulated cable residual life comprehensive estimation methods of electrical characteristics
Huang et al. Plasticized hybrid network solid polymer electrolytes for lithium‐metal batteries
Newsome et al. High-temperature oxidation of SiC-based composite: rate constant calculation from ReaxFF MD simulations, Part II
Sang et al. Improving the thermal properties of ternary carbonates for concentrating solar power through simple chemical modifications by adding sodium hydroxide and nitrate
CN102331377B (en) Method for evaluating creep performance of T/P92 steel
CN102183539A (en) Method for quickly detecting crosslinked polyethylene cable insulating material structure
McGuire et al. A simple method of extrapolating the coexistence curve and predicting the melting point depression curve from cloud point data for polymer-diluent systems
Du et al. Structure and lithium ion diffusion in lithium silicate glasses and at their interfaces with lithium lanthanum titanate crystals
Zhao et al. Demonstration of Mg (NO3) 2· 6H2O-based composite phase change material for practical-scale medium-low temperature thermal energy storage
CN112233729A (en) General form and solving method of polymer insulation pyrolysis reaction mechanism function
Yuan et al. Effects of local thermal accumulation conditions on the thermal characteristics of lithium-ion batteries under high-rate charging
CN110940801A (en) Method for obtaining activation energy of dry insulation equipment based on equal conversion rate method
Brix et al. Modeling char conversion under suspension fired conditions in O2/N2 and O2/CO2 atmospheres
CN105970333B (en) polytetrafluoroethylene fibre
Li et al. Development of a pyrolysis model for poly (vinylidene fluoride-co-hexafluoropropylene) and its application in predicting combustion behaviors
Han et al. Research on thermal management and waste heat utilization of building distributed energy storage battery
Mensah et al. Kinetic parameter estimation from thermogravimetry and microscale combustion calorimetry
Harmathy Design of fire test furnaces
US12024667B2 (en) Cooling liquid composition
Tao et al. Pyrolysis of waste plastics: Effect of heating rate on product yields and oil properties
CN113258155B (en) Method for judging infiltration state of lithium ion battery
CN114062189A (en) Method for measuring rubber content in silicone rubber
CN115406932A (en) Device and method for measuring heat value of hydrocarbon high-temperature pyrolysis gas on line

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20210115

WD01 Invention patent application deemed withdrawn after publication