一种模拟闪络发展中高温对硅橡胶材料裂解的方法技术

本发明专利技术涉及一种模拟闪络发展中高温对硅橡胶材料裂解的方法，包括：S1.基于分子动力学模拟构建了硅氧烷分子和硫化剂分子模型，置于晶胞中进行交联，并对其几何优化和分子动力学平衡，得到硅橡胶交联模型；S2.基于ReaxFF力场模拟闪络工况，构建一定电场强度下不同温度的模拟环境来模拟整个裂解过程；S3.基于分子动力学模拟对硅橡胶在闪络工况下不同温度进行仿真模拟，对其电热裂解产物的种类、数量及速率进行统计；S4.基于裂解温升、裂解产物、裂解速率等参数定量表征材料裂解的程度。分析闪络工况下温升对材料裂解的影响。此发明专利技术公开了一种模拟闪络发展中高温对硅橡胶材料裂解的方法，从微观层面上解释了在一定电场下温度对材料电热裂解的作用机制。料电热裂解的作用机制。料电热裂解的作用机制。

【技术实现步骤摘要】
一种模拟闪络发展中高温对硅橡胶材料裂解的方法


[0001]本专利技术涉及绝缘材料性能劣化研究领域，尤其涉及一种模拟闪络发展中高温对硅橡胶材料裂解的方法。

技术介绍

[0002]复合绝缘子由于具有质量轻、憎水性强、耐污性能好等优点，逐步取代普通的瓷、玻璃绝缘子，已广泛应用于特高压输电线路中。复合绝缘子的绝缘材料主要为有机硅橡胶，硅橡胶材料在闪络放电等工况下受到热、电、化学等因素的作用会迅速发生裂解，导致绝缘性能降低，给输电线路的安全运行带来严重威胁。
[0003]温度是复合绝缘子在运行过程中最重要的影响因素，特别是在闪络工况下，闪络电弧温度会在短时间内发生陡升，最高温度甚至超过10000K，在急剧变化的温度作用下，硅橡胶材料发生裂解，绝缘性能会迅速劣化甚至失效。
[0004]目前国内外学者通过各种宏观试验手段研究复合绝缘子在闪络工况下材料裂解和性能下降情况。不同宏观试验测试从多个角度与因素对材料性能变化做出了合理论断，但由于闪络过程十分短暂，且短时间内温度变化很大，宏观试验很难分析闪络工况下硅橡胶材料的动态裂解过程，更无法在微观层面解释高温对硅橡胶材料电热裂解作用机制。
[0005]近年来，随着计算机技术的发展，反应分子动力学模拟广泛应用于高分子材料设计领域。反应分子动力学模拟可以在原子尺度构建材料，可研究宏观实验方法尚无法考察的现象，从微观角度模拟化学键的自由断裂与生成，揭示原子层面的化学反应过程，同时能够大大减少实验成本和开发周期。

技术实现思路

[0006]本专利技术的目的在于提供一种模拟闪络发展中高温对硅橡胶材料裂解的方法。
[0007]本专利技术的目的是这样实现的：
[0008]一种模拟闪络发展中高温对硅橡胶材料裂解的方法,包括如下步骤：
[0009]S1.基于分子动力学的方法构建硅橡胶分子和硫化剂分子模型，将其置于晶胞中，并完成进行交联过程，随后对模型进行几何优化和分子动力学平衡，得到接近实际构型且稳定的硅橡胶交联模型。
[0010]S2.基于ReaxFF反应力场模拟闪络工况，构建出在一定电场强度下，不同温度的的反应分子动力学模拟环境，进而模拟整个化学反应的过程。
[0011]S3.基于分子动力学模拟方法对硅橡胶在闪络工况下不同高温情况进行仿真模拟，并对其电热裂解产物的种类，以及各种裂解产物出现的数量、速率进行统计。
[0012]S4.基于初始裂解温度变化、裂解产物、裂解速率等参数定量表征材料的电热裂解的程度。分析闪络工况下温度变化对材料裂解的影响。
[0013]进一步的，在所构建的模型中预先设置了硅橡胶分子与硫化剂分子的比例b，b由公式确定，其中N
S
为硫化剂分子中原子的数量，N
EP
为硅橡胶分子中乙烯基团的数量，
预设模型的初始密度并对模型施加周期性边界。
[0014]进一步的，所属步骤S1中硅橡胶分子和硫化剂分子的交联过程具体为：
[0015](1)标记反应原子，所述反应原子为硅橡胶分子中乙烯基团的两条链上的碳自由基，分别标记为R1和R2；
[0016](2)设置交联温度、选择合适的交联力场、截断半径距离R、初始截断半径R0和最大截断半径R
max
和目标交联度S0；
[0017](3)计算被标记为R1的原子和被标记为R2的原子之间的距离D。若D大于R0，则跳到步骤(4)，反之，R1原子和R2原子之间形成共价键，完成两条硅橡胶链的交联，反应完成后统计当前的交联度S；
[0018](4)R0＝R0+R，重复步骤(3)，直到R0大于R
max
或者当前交联度S大于S0，交联结束。
[0019]进一步的，所属步骤S1中，为使模型更接近实际构象，对交联后的硅橡胶模型进行稳定性优化，首先利用共轭梯度法对模型进行能量最小化计算，得到能量相对较低的分子模型，随后在NPT等温等压系综下优化模型密度，使得模型分子排布更接近实际构型。在NVT正则系综下进一步降低体系能量，使结构更稳定。优化后的模型总能量大幅下降且趋于平衡，处于较为稳定的状态。
[0020]进一步的，步骤S2中，基于ReaxFF反应力场构建的分子动力学模拟环境文件中包含基本参数、原子参数、键参数、键角参数、二面角参数、非对角作用项、氢键参数共7个部分，涵括12项能量值，共100余个参数，并在已有的力场文件中添加了13个参数用于表征范德华作用、键角能量，并将此作为初始力场文件。
[0021]进一步的，步骤S3中，对于优化后的模型，本文选用Lammps软件，基于NVT系综进行反应分子动力学模拟。
[0022]进一步的步骤S3中，为保证电热裂解反应的精确性，反应体系单位统一为real，原子类型设置为charge，结合量子化学计算结果，选用的反应力场为CHONSSi.ff。反应体系选用Nose
‑
Hoover法进行温度与压力控制，选用高斯分布算法设置体系初始速度，采用EEM法计算体系电荷分布。
[0023]进一步的，在步骤S3中，在复合绝缘子表面闪络电弧发展阶段，电弧温度出现了一个大陡度的提升，为进一步研究温度对硅橡胶材料裂解的影响，分别设置了温度为4000K、5000K、6000K、7000K、8000K、9000K、10000K的，沿x轴施加0.5kV/cm的电场。为保证计算精度，反应步长设为0.1fs，每隔100步输出当前体系的反应状态。
[0024]进一步的，在步骤S4中，裂解产物中占比较大主产物为氢气H2，其中作为过渡产物的甲烷CH4随着裂解的不断发展会进一步分解为H2，占比极小副产物为水H2O、乙烯C2H4、乙炔C2H2等，产物类型与参考文献中对甲基乙烯基硅氧烷裂解实验中测得的结果相一致。设定不同温度下单位甲基乙烯基硅氧烷裂解产出H2分子数N
H2
，裂解速度V
H2
。副产物分子数N
FCW
，裂解速度V
FCW
。
[0025]裂解产物数量需进行归一化处理，归一化处理公式为N＝n/n
a
，式中N为归一化处理结果，n为裂解过程中的不同产物数，n
a
为体系原子总数。
[0026]系数γ
k
定量表征材料在特定温度下的电热裂解的程度，k定义为特定高温情况。
[0027]γk的计算公式为
[0028][0029]式中α为裂解主产物H2的权重系数，取0.758；β为裂解副产物的权重系数，取0.242；N
tH
为裂解初始甲基乙烯基硅氧烷主链中H原子个数，N
tC
为裂解初始甲基乙烯基硅氧烷主链中C原子个数，N
tO
为裂解初始甲基乙烯基硅氧烷主链中O原子个数；n裂解初始体系中总分子个数；ΔU为材料x轴所施加的电场强度。
[0030]与现有技术相比，本专利技术的有益效果是：
[0031]本专利技术的通过提供一种模拟闪络发展中高温对本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种模拟闪络发展中高温对硅橡胶材料裂解的方法，其特征在于，包括以下步骤：S1.基于分子动力学的方法构建硅橡胶分子和硫化剂分子模型，将其置于晶胞中，并完成进行交联过程，随后对模型进行几何优化和分子动力学平衡，得到接近实际构型且稳定的硅橡胶交联模型；S2.基于ReaxFF反应力场模拟闪络工况，构建出在一定电场强度下，不同温度的的反应分子动力学模拟环境，进而模拟整个化学反应的过程；S3.基于分子动力学模拟方法对硅橡胶在闪络工况下不同高温情况进行仿真模拟，并对其电热裂解产物的种类，以及各种裂解产物出现的数量、速率进行统计；S4.基于初始裂解温度变化、裂解产物、裂解速率等参数定量表征材料的电热裂解的程度。分析闪络工况下温度变化对材料裂解的影响。2.根据权利要求1所述的一种模拟闪络发展中高温对硅橡胶材料裂解的方法，其特征在于：在所构建的模型中预先设置了硅橡胶分子与硫化剂分子的比例b，b由公式确定，其中N
S
为硫化剂分子中原子的数量，N
EP
为硅橡胶分子中乙烯基团的数量，预设模型的初始密度并对模型施加周期性边界。3.根据权利要求1所述的一种模拟闪络发展中高温对硅橡胶材料裂解的方法，其特征在于：所属步骤S1中硅橡胶分子和硫化剂分子的交联过程具体为：(1)标记反应原子，所述反应原子为硅橡胶分子中乙烯基团的两条链上的碳自由基，分别标记为R1和R2；(2)设置交联温度、选择合适的交联力场、截断半径距离R、初始截断半径R0和最大截断半径R
max
和目标交联度S0；(3)计算被标记为R1的原子和被标记为R2的原子之间的距离D。若D大于R0，则跳到步骤(4)，反之，R1原子和R2原子之间形成共价键，完成两条硅橡胶链的交联，反应完成后统计当前的交联度S；(4)R0＝R0+R，重复步骤(3)，直到R0大于R
max
或者当前交联度S大于S0，交联结束。4.根据权利要求1所述的一种模拟闪络发展中高温对硅橡胶材料裂解的方法，其特征在于：所属步骤S1中，为使模型更接近实际构象，对交联后的硅橡胶模型进行稳定性优化，首先利用共轭梯度法对模型进行能量最小化计算，得到能量相对较低的分子模型，随后在NPT等温等压系综下优化模型密度，使得模型分子排布更接近实际构型。在NVT正则系综下进一步降低体系能量，使结构更稳定。优化后的模型总能量大幅下降且趋于平衡，处于较为稳定的状态。5.根据权利要求1所述的一种模拟闪络发展中高温对硅橡胶材料裂解的方法，其特征在于：步骤S2中，基于ReaxFF反应力场构建的分子动力学模拟环境文件中包含基本参数、原子参数、键参数、键角参数、二面角参数、非对角作用项、氢键参数共7个部分，涵括12项能量值，共100余个参数，并在已有的力场文件中添加了13个参数用于表征范德华作用、键角能量，并将此作为初始力场文件。6.根据...

【专利技术属性】
技术研发人员：王路伽，刘嘉伟，包丁煜，马晓华，张建文，
申请(专利权)人：中国矿业大学，
类型：发明
国别省市：