一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法

摘要

本发明涉及一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法，包括：S1.基于分子动力学的方法构建环氧树脂分子和固化剂分子模型，将其置于晶胞中，并完成进行交联过程，随后对模型进行几何优化和分子动力学平衡，得到纯环氧树脂模型；S2.基于分子动力学的方法构建了纳米SiO2球形颗粒和不同支化度的超支化聚酯，并在SiO2表面接枝不同支化度的超支化聚酯；S3.将所述的SiO2置于晶胞中心，添加环氧树脂分子和固化剂分子，完成交联过程，并对模型进行几何优化和分子动力学平衡，得到不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型；S4.分别计算出纯环氧树脂模型和不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型的导热率，定量表征材料的导热性能。此发明公开了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法，从微观层面上解释了导热性能变化机制。

说明书

技术领域

本发明涉及绝缘材料改性技术领域，尤其涉及一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法。

背景技术

环氧树脂因其优良的电气性能和机械性能广泛应用于绝缘套管和干式变压器等电力设备中。随着我国电网电压等级的不断提高和电力设备的体积不断缩小，电力设备的发热问题尤为突出，严重时会造成电力设备的停运，引发大规模停电，这就对环氧树脂导热性能提出了更为严苛的要求，但本征环氧树脂的导热率仅为0.17～0.21W/(m·K)-

向环氧树脂中填充高导热率的纳米颗粒可以提升环氧树脂的导热率，但纳米颗粒表面能高，与聚合物的相容性差，在掺杂过程中易出现团聚现象，往往达不到预期效果。超支化聚酯作为一种高度支化的大分子表面改性剂，用其接枝纳米颗粒表面，能有效加强纳米颗粒与聚合物的界面结合作用，提升分散性，导热率提升效果明显。

将超支化聚酯接枝到纳米颗粒表面主要有“grafting from”法和“grafting to”法，其中“grafting from”法利用纳米颗粒表面的活性基团作为引发点来引发超支化单体进行聚合反应形成高度支化的聚酯，在接枝过程中受空间位阻的影响较小因而接枝率高。但由于聚合反应程度无法控制，因而纳米颗粒表面的聚酯的支化度存在差异，而支化度的不同对纳米复合材料的导热率影响很大，故有必要研究支化度对复合纳米环氧树脂的导热率影响。传统材料制备实验无法准确控制超支化聚酯的支化度，同时缺乏从微观层面上揭示材料导热性能变化的手段。

近年来，随着计算机技术的发展，分子模拟广泛应用于高分子材料设计领域。分子模拟可以在原子尺度构建材料，大大减少实验成本和开发周期，同时结合材料的微观结构与宏观性能，揭示材料改性的机制。

发明内容

本发明的目的在于提供了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法。

本发明的目的是这样实现的：

一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法，包括如下步骤：

S1：基于分子动力学的方法构建环氧树脂分子和固化剂分子模型，将其放置于晶胞中，完成交联过程，随后对模型进行几何优化和分子动力学平衡，得到纯环氧树脂模型；

S2：基于分子动力学的方法构建纳米SiO

S3：将所述的SiO

S4：分别计算出纯环氧树脂模型和不同支化度下的超支化聚酯改性SiO

进一步的，在所构建的模型中预先设置了环氧树脂分子与固化剂分子的比例k，k由公式

进一步的，所述步骤S1或S3中环氧树脂分子和固化剂分子的交联过程具体为：

(1)标记反应原子，所述反应原子为环氧树脂中环氧基团的碳原子和固化剂中氨基团的氮原子，分别标记为R1和R2；

(2)设置交联温度、截断半径距离R、初始截断半径R

(3)计算被标记为R1的原子和被标记为R2的原子之间的距离D。若R大于R

(4)R

进一步的，所述步骤S2中设置了纳米SiO

进一步的，所述步骤S2中构建的超支化聚酯单体为AB

羧基的转化率x由表达式

比例系数r比例系数r表征线性聚合反应和支化聚合反应的竞争关系，由表达式

接枝方法采用“Grafting form”模式，具体过程为：

(1)在SiO

(2)硅烷偶联剂末端氨基与超支化聚酯单体中羧基进行反应，接入单体，随后完成接枝单体中氨基与未接枝的单体中羧基反应，进而形成超支化聚酯。保持羧基的转化率x不变，改变比例系数r，接枝不同支化度的超支化聚酯。

进一步的，所述步骤S3中预先设置了SiO

进一步的，所述步骤S1或S3中，几何优化采用共轭梯度法，分子动力学模拟先后在NVT和NPT系综下进行，控温方法选择Nose，控压方法选择Berendsen，力场选择COSPASSII，分别选用Ewald和Atom based法来描述分子间静电力和范德华力。

进一步的，所述步骤S4中，导热率的计算模型需利用Symmetry模块中的Spuercell功能对原有模型进行扩晶胞处理，构建1×1×3的模型，导热率的计算方法为反向扰动非平衡动力学方法，通过交换模型高温层和低温层中的不同速度粒子实现能量不断从两侧传递到中间区域；

导热率λ的计算公式为：

式中其中J为模型z方向上的能量通量，

J的计算公式为

式中N是该区域中动能交换的原子数，m和v分别是原子的质量和速度，下标1和2表示速度互换的最热和最冷的原子，t是在模拟中施加热通量的时间间隔，L

温度T的计算公式为

式中k

与现有技术相比，本发明的有益效果是：

本发明提出了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法，针对实际试验过程纳米颗粒表面采用“grafting form”法接枝超支化聚酯时支化度难以控制的特点，通过分子动力学的方法，通过改变支化度定量研究了支化度对材料导热率的影响，节约了试验成本和开发周期，更好地应用于生产实践。

附图说明

图1为本发明的操作流程图。

图2为环氧树脂分子和固化剂分子交联原理图。

图3为纳米颗粒接枝超支化聚酯原理图。

图4为不同支化度的超支化聚酯修饰的纳米SiO

图5为不同支化度下的超支化聚酯改性SiO

图6为导热率计算示意图。

图7为能量通量随模拟时间的变化。

图8为模型z方向上温度梯度分布图。

具体实施方式

下面结合附图和实施例，对本发明的具体实施方式作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案，而不能以此来限制发明的保护范围。

本发明的目的是针对现有技术的缺陷，提供了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法。

实施例

本实施例提供了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法，如图1所示，包括步骤：

S1：基于分子动力学的方法构建了环氧树脂分子和固化剂分子模型，环氧树脂分子选用双酚A二缩水甘油醚(DGEBA)，DGEBA中环氧基团为2个，固化剂为4,4’-二氨基二苯砜(44DDS)，44DDS中有2个氨基团且含有4个氢原子，故两者比例k＝2，初始晶胞中DGEBA分子为50条，44DDS分子为25条，初始密度设置为0.6g/cm

(1)标记反应原子，所述反应原子为环氧树脂中环氧基团的碳原子和固化剂中氨基团的氮原子，分别标记为R1和R2；

(2)设置交联温度为580K，截断半径距离R为

(3)计算被标记为R1的原子和被标记为R2的原子之间的距离D。若D大于R

(4)R

图5(a)为纯环氧树脂模型；

紧接着对模型进行几何结构优化和分子动力学模拟，具体为：

(1)几何结构优化：基于交联后的模型，使用共轭梯度法进行几何优化，获得能量最低构象。能量和力的最大容差分别设置为0.001kcal/mol和

(2)分子动力学模拟：基于能量最低构象，首先在定容定温(NVT)系综中，在600K下进行200ps的分子动力学模拟，然后在300K下再进行500ps的分子动力学模拟；之后在定压定温(NPT)系综中，温度设置为300K，压力设置为1atm，进行500ps的分子动力学模拟。模拟过程中力场选择COMPASSII，控温方法选择Nose，控压方法选择Berendsen，分别选用Ewald和Atom based法来描述分子间静电力和范德华力。

S2：基于分子动力学的方法构建了半径为

(1)在SiO

(2)超支化聚酯单体选用3,5-二氨基苯甲酸(DABA)，DABA同时含有2个氨基和1个羧基，硅烷偶联剂末端氨基与超支化聚酯单体中羧基进行反应，接入单体，随后完成接枝单体中氨基与未接枝的单体中羧基反应，进而形成超支化聚酯。

图4分别接枝了超支化聚酯的SiO

S3：将所述的SiO

S4：采用基于NEMD改进的反向扰动非平衡动力学方法(RNEMD)来计算不同支化度下的超支化聚酯的改性SiO

以纯环氧树脂为例，基于分子动力学，模型导热率λ的计算公式为：

式中其中J为模型z方向上的能量通量，

式中N是该区域中动能交换的原子数，m和v分别是原子的质量和速度，下标1和2表示速度互换的最热和最冷的原子，t是在模拟中施加热通量的时间间隔，L

温度T的计算公式为

式中k

在上述实施例中，仅对本发明进行示范性描述，但是本领域技术人员在阅读本专利申请后可以在不脱离本发明的精神和范围的情况3下对本发明进行各种修改。