CN113257366B - 一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 - Google Patents
一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 Download PDFInfo
- Publication number
- CN113257366B CN113257366B CN202110557686.4A CN202110557686A CN113257366B CN 113257366 B CN113257366 B CN 113257366B CN 202110557686 A CN202110557686 A CN 202110557686A CN 113257366 B CN113257366 B CN 113257366B
- Authority
- CN
- China
- Prior art keywords
- epoxy resin
- model
- thermal conductivity
- molecular dynamics
- hyperbranched polyester
- 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.)
- Active
Links
- 239000003822 epoxy resin Substances 0.000 title claims abstract description 65
- 229920000647 polyepoxide Polymers 0.000 title claims abstract description 65
- 238000000034 method Methods 0.000 title claims abstract description 61
- 229920006150 hyperbranched polyester Polymers 0.000 title claims abstract description 57
- 238000000329 molecular dynamics simulation Methods 0.000 title claims abstract description 43
- 239000002114 nanocomposite Substances 0.000 title claims abstract description 21
- 238000004132 cross linking Methods 0.000 claims abstract description 25
- 239000003795 chemical substances by application Substances 0.000 claims abstract description 19
- 229910004298 SiO 2 Inorganic materials 0.000 claims abstract description 18
- 230000008569 process Effects 0.000 claims abstract description 16
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 239000002131 composite material Substances 0.000 claims abstract description 10
- 239000000463 material Substances 0.000 claims abstract description 9
- 239000012798 spherical particle Substances 0.000 claims abstract description 5
- 229920000728 polyester Polymers 0.000 claims abstract description 4
- 125000004429 atom Chemical group 0.000 claims description 30
- 238000006243 chemical reaction Methods 0.000 claims description 24
- 239000000178 monomer Substances 0.000 claims description 21
- 125000003178 carboxy group Chemical group [H]OC(*)=O 0.000 claims description 20
- 238000004364 calculation method Methods 0.000 claims description 15
- 125000003277 amino group Chemical group 0.000 claims description 13
- 125000003700 epoxy group Chemical group 0.000 claims description 12
- 229910052799 carbon Inorganic materials 0.000 claims description 9
- 125000004432 carbon atom Chemical group C* 0.000 claims description 9
- 125000002887 hydroxy group Chemical group [H]O* 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 7
- 239000006087 Silane Coupling Agent Substances 0.000 claims description 6
- 229910052757 nitrogen Inorganic materials 0.000 claims description 6
- 125000004433 nitrogen atom Chemical group N* 0.000 claims description 6
- 238000006116 polymerization reaction Methods 0.000 claims description 6
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 5
- 125000004430 oxygen atom Chemical group O* 0.000 claims description 5
- 238000004088 simulation Methods 0.000 claims description 5
- GOJUJUVQIVIZAV-UHFFFAOYSA-N 2-amino-4,6-dichloropyrimidine-5-carbaldehyde Chemical group NC1=NC(Cl)=C(C=O)C(Cl)=N1 GOJUJUVQIVIZAV-UHFFFAOYSA-N 0.000 claims description 4
- 239000002245 particle Substances 0.000 claims description 4
- 238000005411 Van der Waals force Methods 0.000 claims description 3
- 230000010261 cell growth Effects 0.000 claims description 3
- 238000002939 conjugate gradient method Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 239000004848 polyfunctional curative Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 125000000467 secondary amino group Chemical class [H]N([*:1])[*:2] 0.000 claims description 3
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical group [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 2
- 125000002915 carbonyl group Chemical group [*:2]C([*:1])=O 0.000 claims description 2
- 238000005984 hydrogenation reaction Methods 0.000 claims description 2
- 230000000737 periodic effect Effects 0.000 claims description 2
- 125000002924 primary amino group Chemical class [H]N([H])* 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 claims 1
- 230000008859 change Effects 0.000 abstract description 3
- 230000007246 mechanism Effects 0.000 abstract description 2
- 239000002105 nanoparticle Substances 0.000 description 9
- LCFVJGUPQDGYKZ-UHFFFAOYSA-N Bisphenol A diglycidyl ether Chemical compound C=1C=C(OCC2OC2)C=CC=1C(C)(C)C(C=C1)=CC=C1OCC1CO1 LCFVJGUPQDGYKZ-UHFFFAOYSA-N 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 239000004593 Epoxy Substances 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000009191 jumping Effects 0.000 description 2
- 238000012900 molecular simulation Methods 0.000 description 2
- 229920000642 polymer Polymers 0.000 description 2
- 150000003141 primary amines Chemical class 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- UENRXLSRMCSUSN-UHFFFAOYSA-N 3,5-diaminobenzoic acid Chemical group NC1=CC(N)=CC(C(O)=O)=C1 UENRXLSRMCSUSN-UHFFFAOYSA-N 0.000 description 1
- MQJKPEGWNLWLTK-UHFFFAOYSA-N Dapsone Chemical compound C1=CC(N)=CC=C1S(=O)(=O)C1=CC=C(N)C=C1 MQJKPEGWNLWLTK-UHFFFAOYSA-N 0.000 description 1
- 235000010627 Phaseolus vulgaris Nutrition 0.000 description 1
- 244000046052 Phaseolus vulgaris Species 0.000 description 1
- 238000005054 agglomeration Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003631 expected effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- LNEPOXFFQSENCJ-UHFFFAOYSA-N haloperidol Chemical compound C1CC(O)(C=2C=CC(Cl)=CC=2)CCN1CCCC(=O)C1=CC=C(F)C=C1 LNEPOXFFQSENCJ-UHFFFAOYSA-N 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000977 initiatory effect Effects 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000003607 modifier Substances 0.000 description 1
- 239000002861 polymer material Substances 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 229910052710 silicon Inorganic materials 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational 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
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Epoxy Resins (AREA)
Abstract
本发明涉及一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,包括:S1.基于分子动力学的方法构建环氧树脂分子和固化剂分子模型,将其置于晶胞中,并完成进行交联过程,随后对模型进行几何优化和分子动力学平衡,得到纯环氧树脂模型;S2.基于分子动力学的方法构建了纳米SiO2球形颗粒和不同支化度的超支化聚酯,并在SiO2表面接枝不同支化度的超支化聚酯;S3.将所述的SiO2置于晶胞中心,添加环氧树脂分子和固化剂分子,完成交联过程,并对模型进行几何优化和分子动力学平衡,得到不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型;S4.分别计算出纯环氧树脂模型和不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型的导热率,定量表征材料的导热性能。此发明公开了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,从微观层面上解释了导热性能变化机制。
Description
技术领域
本发明涉及绝缘材料改性技术领域,尤其涉及一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法。
背景技术
环氧树脂因其优良的电气性能和机械性能广泛应用于绝缘套管和干式变压器等电力设备中。随着我国电网电压等级的不断提高和电力设备的体积不断缩小,电力设备的发热问题尤为突出,严重时会造成电力设备的停运,引发大规模停电,这就对环氧树脂导热性能提出了更为严苛的要求,但本征环氧树脂的导热率仅为0.17~0.21W/(m·K)-1,无法满足电力设备的散热需求。
向环氧树脂中填充高导热率的纳米颗粒可以提升环氧树脂的导热率,但纳米颗粒表面能高,与聚合物的相容性差,在掺杂过程中易出现团聚现象,往往达不到预期效果。超支化聚酯作为一种高度支化的大分子表面改性剂,用其接枝纳米颗粒表面,能有效加强纳米颗粒与聚合物的界面结合作用,提升分散性,导热率提升效果明显。
将超支化聚酯接枝到纳米颗粒表面主要有“grafting from”法和“grafting to”法,其中“grafting from”法利用纳米颗粒表面的活性基团作为引发点来引发超支化单体进行聚合反应形成高度支化的聚酯,在接枝过程中受空间位阻的影响较小因而接枝率高。但由于聚合反应程度无法控制,因而纳米颗粒表面的聚酯的支化度存在差异,而支化度的不同对纳米复合材料的导热率影响很大,故有必要研究支化度对复合纳米环氧树脂的导热率影响。传统材料制备实验无法准确控制超支化聚酯的支化度,同时缺乏从微观层面上揭示材料导热性能变化的手段。
近年来,随着计算机技术的发展,分子模拟广泛应用于高分子材料设计领域。分子模拟可以在原子尺度构建材料,大大减少实验成本和开发周期,同时结合材料的微观结构与宏观性能,揭示材料改性的机制。
发明内容
本发明的目的在于提供了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法。
本发明的目的是这样实现的:
一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,包括如下步骤:
S1:基于分子动力学的方法构建环氧树脂分子和固化剂分子模型,将其放置于晶胞中,完成交联过程,随后对模型进行几何优化和分子动力学平衡,得到纯环氧树脂模型;
S2:基于分子动力学的方法构建纳米SiO2球形颗粒和超支化聚酯单体,并在SiO2表面接枝不同支化度的超支化聚酯;
S3:将所述的SiO2置于晶胞中心,添加环氧树脂分子和固化剂分子,完成交联过程,并对模型进行几何优化和分子动力学平衡,得到不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型;
S4:分别计算出纯环氧树脂模型和不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型的导热率,定量表征材料的导热性能。
进一步的,在所构建的模型中预先设置了环氧树脂分子与固化剂分子的比例k,k由公式确定,Nh为固化剂分子中氨基团中氢原子的数量,其中Nep为环氧树脂中环氧基团的数量,预设模型的初始密度并对模型施加周期性边界。
进一步的,所述步骤S1或S3中环氧树脂分子和固化剂分子的交联过程具体为:
(1)标记反应原子,所述反应原子为环氧树脂中环氧基团的碳原子和固化剂中氨基团的氮原子,分别标记为R1和R2;
(2)设置交联温度、截断半径距离R、初始截断半径R0和最大截断半径Rmax和目标交联度S0;
(3)计算被标记为R1的原子和被标记为R2的原子之间的距离D。若R大于R0,则跳到步骤(4),反之,R1原子和R2原子之间形成共价键,即首先是环氧基团中碳原子与伯胺中氮原子反应,其次是环氧基团中的碳原子与第一步反应之后的产物仲胺反应,反应完成后统计当前的交联度S;
(4)R0=R0+R,重复步骤(3),直到R0大于Rmax或者当前交联度S大于S0,交联结束。
进一步的,所述步骤S2中设置了纳米SiO2的半径,并对SiO2进行氢化处理,即其表面上的硅原子上的断裂键并与羟基结合,氧原子上的断裂键与氢原子结合。
进一步的,所述步骤S2中构建的超支化聚酯单体为AB2型单体,含有2个氨基和1个羧基,支化度Db由公式Db=xe-r,其中x为羧基的转化率,r为比例系数,x的限定区间为[0.6,0.9],r的限定区间为[0.2,0.8];
羧基的转化率x由表达式确定,其中NCO为超支化聚酯中羰基的数量,NCOOH为超支化聚酯中羧基的数量;
比例系数r比例系数r表征线性聚合反应和支化聚合反应的竞争关系,由表达式确定,其中kL表示线性单元反应速率,kT表示支化单元反应速率;
接枝方法采用“Grafting form”模式,具体过程为:
(1)在SiO2表面的羟基中氧原子上接枝硅烷偶联剂,引入氨基;
(2)硅烷偶联剂末端氨基与超支化聚酯单体中羧基进行反应,接入单体,随后完成接枝单体中氨基与未接枝的单体中羧基反应,进而形成超支化聚酯。保持羧基的转化率x不变,改变比例系数r,接枝不同支化度的超支化聚酯。
进一步的,所述步骤S3中预先设置了SiO2的体积分数Dv,Dv由表达式其中Vpar为SiO2的体积,V0为纳米复合环氧树脂的体积,接枝密度ρ由表达式/>确定,其中NOH1为SiO2接枝超支化聚酯的羟基的数量,NOH为SiO2上羟基的总数目。
进一步的,所述步骤S1或S3中,几何优化采用共轭梯度法,分子动力学模拟先后在NVT和NPT系综下进行,控温方法选择Nose,控压方法选择Berendsen,力场选择COSPASSII,分别选用Ewald和Atom based法来描述分子间静电力和范德华力。
进一步的,所述步骤S4中,导热率的计算模型需利用Symmetry模块中的Spuercell功能对原有模型进行扩晶胞处理,构建1×1×3的模型,导热率的计算方法为反向扰动非平衡动力学方法,通过交换模型高温层和低温层中的不同速度粒子实现能量不断从两侧传递到中间区域;
导热率λ的计算公式为:
式中其中J为模型z方向上的能量通量,为模型z方向上温度梯度;
J的计算公式为
式中N是该区域中动能交换的原子数,m和v分别是原子的质量和速度,下标1和2表示速度互换的最热和最冷的原子,Δt是在模拟中施加热通量的时间间隔,Lx和Ly是x轴和y轴的长度;
温度T的计算公式为
式中kB是玻尔兹曼常数,温度梯度由模型z方向上的温度分布线性拟合求得。
与现有技术相比,本发明的有益效果是:
本发明提出了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,针对实际试验过程纳米颗粒表面采用“grafting form”法接枝超支化聚酯时支化度难以控制的特点,通过分子动力学的方法,通过改变支化度定量研究了支化度对材料导热率的影响,节约了试验成本和开发周期,更好地应用于生产实践。
附图说明
图1为本发明的操作流程图。
图2为环氧树脂分子和固化剂分子交联原理图。
图3为纳米颗粒接枝超支化聚酯原理图。
图4为不同支化度的超支化聚酯修饰的纳米SiO2颗粒。
图5为不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型。
图6为导热率计算示意图。
图7为能量通量随模拟时间的变化。
图8为模型z方向上温度梯度分布图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步描述。以下实施例仅用于更加清楚地说明本发明的技术方案,而不能以此来限制发明的保护范围。
本发明的目的是针对现有技术的缺陷,提供了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法。
实施例
本实施例提供了一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,如图1所示,包括步骤:
S1:基于分子动力学的方法构建了环氧树脂分子和固化剂分子模型,环氧树脂分子选用双酚A二缩水甘油醚(DGEBA),DGEBA中环氧基团为2个,固化剂为4,4’-二氨基二苯砜(44DDS),44DDS中有2个氨基团且含有4个氢原子,故两者比例k=2,初始晶胞中DGEBA分子为50条,44DDS分子为25条,初始密度设置为0.6g/cm3,环氧树脂分子和固化剂分子的交联原理如图2所示,具体为:
(1)标记反应原子,所述反应原子为环氧树脂中环氧基团的碳原子和固化剂中氨基团的氮原子,分别标记为R1和R2;
(2)设置交联温度为580K,截断半径距离R为初始截断半径为R0为/>最大截断半径Rmax为/>预设目标交联度S0=90%;
(3)计算被标记为R1的原子和被标记为R2的原子之间的距离D,若D大于R0,则跳到步骤(d),反之,R1原子和R2原子之间形成共价键,即第一步是环氧基团中碳原子与伯胺中氮原子反应,第二步是环氧基团中的碳原子与第一步反应之后的产物仲胺反应,反应完成后计算目前的交联度S;
(4)R0=R0+R,重复步骤(c),直到R0>Rmax或者S>S0,交联过程结束;
图5(a)为纯环氧树脂模型;
紧接着对模型进行几何结构优化和分子动力学模拟,具体为:
(1)几何结构优化:基于交联后的模型,使用共轭梯度法进行几何优化,获得能量最低构象。能量和力的最大容差分别设置为0.001kcal/mol和最大迭代次数为10000。
(2)分子动力学模拟:基于能量最低构象,首先在定容定温(NVT)系综中,在600K下进行200ps的分子动力学模拟,然后在300K下再进行500ps的分子动力学模拟;之后在定压定温(NPT)系综中,温度设置为300K,压力设置为1atm,进行500ps的分子动力学模拟。模拟过程中力场选择COMPASSII,控温方法选择Nose,控压方法选择Berendsen,分别选用Ewald和Atom based法来描述分子间静电力和范德华力。
S2:基于分子动力学的方法构建了半径为的纳米SiO2球形颗粒和超支化聚酯单体分子模型,羧基的转化率x=0.75,比例系数r为0.25,0.5和0.75,对应的支化度分别58.4%、45.5%和35.4%,之后在SiO2表面接枝不同支化度的超支化聚酯,接枝原理如图3所示,具体为:
(1)在SiO2表面的羟基中氧原子上接枝硅烷偶联剂KH550,引入氨基;
(2)超支化聚酯单体选用3,5-二氨基苯甲酸(DABA),DABA同时含有2个氨基和1个羧基,硅烷偶联剂末端氨基与超支化聚酯单体中羧基进行反应,接入单体,随后完成接枝单体中氨基与未接枝的单体中羧基反应,进而形成超支化聚酯。
图4分别接枝了超支化聚酯的SiO2球形颗粒,超支化聚酯的支化度分别为58.4%、45.5%和35.4%;
S3:将所述的SiO2置于晶胞中心,添加环氧树脂单体和固化剂,纳米SiO2的体积分数为7.2%,进行如步骤S1相同的交联过程、几何优化和分子动力学平衡,得到支化度分别为58.4%、45.5%和35.4%的超支化聚酯改性的SiO2/环氧树脂复合模型,如图5(b)、(c)、(d)所示;
S4:采用基于NEMD改进的反向扰动非平衡动力学方法(RNEMD)来计算不同支化度下的超支化聚酯的改性SiO2/环氧树脂复合模型的导热率,RNEMD的计算模型示意图如图6所示,将模型沿着z轴方向分为40层,两侧为高温层,中间为低温层,通过交换高温层和低温层之间的粒子实现能量不断从两侧传递到中间区域,最终形成稳定的热流和温度梯度。为了保证计算结果的准确性,利用Symmetry模块中的Spuercell功能对原有模型进行扩晶胞处理,构建1×1×3的模型。
以纯环氧树脂为例,基于分子动力学,模型导热率λ的计算公式为:
式中其中J为模型z方向上的能量通量,为模型z方向上温度梯度;J的计算公式为:
式中N是该区域中动能交换的原子数,m和v分别是原子的质量和速度,下标1和2表示速度互换的最热和最冷的原子,Δt是在模拟中施加热通量的时间间隔,Lx和Ly是x轴和y轴的长度,图7为分子动力学中J随时间的变化情况;
温度T的计算公式为
式中kB是玻尔兹曼常数,温度梯度由模型z方向上的温度分布线性拟合求得,如图8所示。
在上述实施例中,仅对本发明进行示范性描述,但是本领域技术人员在阅读本专利申请后可以在不脱离本发明的精神和范围的情况3下对本发明进行各种修改。
Claims (8)
1.一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于,包括以下步骤:
S1:基于分子动力学的方法构建环氧树脂分子和固化剂分子模型,将其放置于晶胞中,完成交联过程,随后对模型进行几何优化和分子动力学平衡,得到纯环氧树脂模型;
S2:基于分子动力学的方法构建纳米SiO2球形颗粒和超支化聚酯单体,并在SiO2表面接枝不同支化度的超支化聚酯;
S3:将所述的SiO2置于晶胞中心,添加环氧树脂分子和固化剂分子,完成交联过程,并对模型进行几何优化和分子动力学平衡,得到不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型;
S4:分别计算出纯环氧树脂模型和不同支化度下的超支化聚酯改性SiO2/环氧树脂复合模型的导热率,定量表征材料的导热性能。
2.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:在所构建的模型中预先设置了环氧树脂分子与固化剂分子的比例k,k由公式确定,其中Nh为固化剂分子中氨基团中氢原子的数量,Nep为环氧树脂分子中环氧基团的数量,预设模型的初始密度并对模型施加周期性边界。
3.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:所述步骤S1或S3中环氧树脂分子和固化剂分子的交联过程具体为:
(1)标记反应原子,所述反应原子为环氧树脂分子中环氧基团的碳原子和固化剂分子中氨基团的氮原子,分别标记为R1和R2;
(2)设置交联温度、截断半径距离R、初始截断半径R0和最大截断半径Rmax和目标交联度S0;
(3)计算被标记为R1的原子和被标记为R2的原子之间的距离D,若D大于R0,则跳到步骤(4),反之,R1原子和R2原子之间形成共价键,即首先是环氧基团中碳原子与伯胺中氮原子反应,其次是环氧基团中的碳原子与第一步反应之后的产物仲胺反应,反应完成后统计当前的交联度S;
(4)R0=R0+R,重复步骤(3),直到R0大于Rmax或者当前交联度S大于S0,交联结束。
4.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:所述步骤S2中设置了纳米SiO2的半径,并对SiO2进行氢化处理,即其表面上的硅原子上的断裂键与羟基结合,氧原子上的断裂键与氢原子结合。
5.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:所述步骤S2中构建的超支化聚酯单体为AB2型单体,含有2个氨基和1个羧基,支化度Db由公式Db=xe-r确定,其中x为羧基的转化率,r为比例系数,x的限定区间为[0.6,0.9],r的限定区间为[0.2,0.8];
羧基的转化率x由表达式确定,其中NCO为超支化聚酯中羰基的数量,NCOOH为超支化聚酯中羧基的数量;
比例系数r表征线性聚合反应和支化聚合反应的竞争关系,由表达式确定,其中kL表示线性单元反应速率,kT表示支化单元反应速率;
接枝方法采用“Grafting form”模式,具体过程为:
(1)在SiO2表面的羟基中氧原子上接枝硅烷偶联剂,引入氨基;
(2)硅烷偶联剂末端氨基与超支化聚酯单体中羧基进行反应,接入单体,随后完成接枝单体中氨基与未接枝的单体中羧基聚合反应,形成超支化聚酯,保持羧基的转化率x不变,改变比例系数r,接枝不同支化度的超支化聚酯。
6.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:所述步骤S3中预先设置了SiO2的体积分数Dv,Dv由表达式其中Vpar为SiO2的体积,V0为纳米复合环氧树脂的体积,接枝密度ρ由表达式/>确定,其中NOH1为SiO2接枝超支化聚酯的羟基的数量,NOH为SiO2上羟基的总数目。
7.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:所述步骤S1或S3中,几何优化采用共轭梯度法,分子动力学模拟先后在NVT和NPT系综下进行,控温方法选择Nose,控压方法选择Berendsen,力场选择COSPASSII,分别选用Ewald和Atom based法来描述分子间静电力和范德华力。
8.根据权利要求1所述的一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法,其特征在于:所述步骤S4中,导热率的计算模型需利用Symmetry模块中的Spuercell功能对原有模型进行扩晶胞处理,构建1×1×3的模型,导热率的计算方法为反向扰动非平衡动力学方法,通过交换模型高温层和低温层中的不同速度粒子实现能量不断从两侧传递到中间区域;
导热率λ的计算公式为:
式中其中J为模型z方向上的能量通量,为模型z方向上温度梯度;
J的计算公式为:
式中N是该区域中动能交换的原子数,m和v分别是原子的质量和速度,下标1和2表示速度互换的最热和最冷的原子,Δt是在模拟中施加热通量的时间间隔,Lx和Ly是x轴和y轴的长度;
温度T的计算公式为:
式中kB是玻尔兹曼常数,温度梯度由模型z方向上的温度分布线性拟合求得。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110557686.4A CN113257366B (zh) | 2021-05-21 | 2021-05-21 | 一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110557686.4A CN113257366B (zh) | 2021-05-21 | 2021-05-21 | 一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113257366A CN113257366A (zh) | 2021-08-13 |
CN113257366B true CN113257366B (zh) | 2024-01-26 |
Family
ID=77183622
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110557686.4A Active CN113257366B (zh) | 2021-05-21 | 2021-05-21 | 一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113257366B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114882958B (zh) * | 2022-05-16 | 2024-07-09 | 天津大学 | 全面反映交联的酸酐固化环氧树脂分子动力学模拟方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201413703Y (zh) * | 2009-05-27 | 2010-02-24 | 苏州巨峰金属线缆有限公司 | 高导热涤纶双玻璃丝铜扁线 |
CN111505064A (zh) * | 2020-04-28 | 2020-08-07 | 中国矿业大学 | 一种催化燃烧式甲烷传感器服役状态评价方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8293831B2 (en) * | 2008-10-30 | 2012-10-23 | E I Du Pont De Nemours And Company | Thermoplastic composition including thermally conductive filler and hyperbranched polyesteramide |
WO2012066520A1 (en) * | 2010-11-19 | 2012-05-24 | Ecole Polytechnique Federale De Lausanne (Epfl) | Dual cure compositions, related hybrid nanocomposite materials and dual cure process for producing same |
-
2021
- 2021-05-21 CN CN202110557686.4A patent/CN113257366B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201413703Y (zh) * | 2009-05-27 | 2010-02-24 | 苏州巨峰金属线缆有限公司 | 高导热涤纶双玻璃丝铜扁线 |
CN111505064A (zh) * | 2020-04-28 | 2020-08-07 | 中国矿业大学 | 一种催化燃烧式甲烷传感器服役状态评价方法 |
Non-Patent Citations (7)
Title |
---|
Enhanced mechanical and thermal properties of mono-component high performance epoxy resin by blending with hydroxyl terminal polyethersulfone;Minqiang Jiang等;《Polymer Testing》;302-309 * |
Thermomechanical properties of silica-epoxy nanocomposite modified by hyperbranched polyester: A molecular dynamics simulation;Jianwen Zhang等;《High Performance Polymers》;1-12 * |
环氧树脂/超支化聚酯/纳米SiO2复合材料的制备及性能;蒋玉梅等;《高分子材料科学与工程》;134-137 * |
端基对超支化聚酯改性纳米SiO2/环氧树脂绝缘性能的影响研究;卜梦晨;《中国优秀硕士学位论文全文数据库工程科技II辑》;C042-317 * |
纳米ZnO与SiO2改性环氧树脂的绝缘特性及导热性能;杨国清等;《大电机技术》;38-43 * |
超支化聚酯改性SiO2/环氧树脂复合材料热力学性能的分子动力学模拟;王冬伟;《中国优秀硕士学位论文全文数据库》;B020-285 * |
面向虚拟仿真实验的匝间短路故障异步电动机模型设计与实现;王攀攀;邓先明;徐瑞东;张建文;;实验室研究与探索(第11期);111-117 * |
Also Published As
Publication number | Publication date |
---|---|
CN113257366A (zh) | 2021-08-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Progress in low dielectric polyimide film–A review | |
Zhang et al. | Strengthened epoxy resin with hyperbranched polyamine-ester anchored graphene oxide via novel phase transfer approach | |
CN113257366B (zh) | 一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 | |
Zhao et al. | A comparative study on 3D printed silicone-epoxy/acrylate hybrid polymers via pure photopolymerization and dual-curing mechanisms | |
CN108314777B (zh) | 一种环氧/酸酐固化物交联分子模型构建及其玻璃化温度提取方法 | |
CN109265922B (zh) | 一种高韧性自催化环氧树脂及制备方法 | |
CN103642050A (zh) | 具有poss基的树枝状聚合物及其制备方法 | |
Yu et al. | Highly populated and nearly monodispersed nanosilica particles in an organic medium and their epoxy nanocomposites | |
JP5026244B2 (ja) | バイポーラプレート製造用成形組成物 | |
JPH0329259B2 (zh) | ||
CN102775755A (zh) | 一种聚芳醚腈和羰基铁粉复合磁性材料及其制备方法 | |
CN111499837B (zh) | 一种自修复原位荧光示踪型耐热耐候环氧树脂及其制备方法 | |
CN113088033A (zh) | 基于原位聚合的导电/导热石墨烯/形状记忆聚合物复合材料、制备方法及应用 | |
Mei et al. | 4D printing of polyamide 1212 based shape memory thermoplastic polyamide elastomers by selective laser sintering | |
Li et al. | Vat polymerization-based 3D printing of nanocomposites: A mini review | |
Guo et al. | Advances and outlook in modified graphene oxide (GO)/epoxy composites for mechanical applications | |
Zhou et al. | UV-assisted direct ink writing 4D printing of benzoxazine/epoxy thermosets | |
Guan et al. | Self-assembled proton conduction networks consisting of SPEEK, NH2-POSS, and IL with enhanced proton conduction and decreased IL loss | |
Santiago et al. | Enhancement of epoxy thermosets with hyperbranched and multiarm star polymers: a review | |
Li et al. | Hybrid filler with nanoparticles grown in situ on the surface for the modification of thermal conductive and insulating silicone rubber | |
Chaudhary et al. | Study of the cross-linking density effect on the mechanical properties of h-BNNS reinforced epoxy nanocomposite part-1: a molecular dynamics simulation | |
Wang et al. | Optimizing curing process of graphene oxide/waterborne epoxy blends by curing kinetics simulation considering the coupling of heat conduction and curing reaction | |
CN113045715A (zh) | 一种可水洗耐高温导电的3d打印光固化材料及其打印工艺 | |
Bekeshev et al. | The Influence of Pristine and Aminoacetic Acid-Treated Aluminum Nitride on the Structure, Curing Processes, and Properties of Epoxy Nanocomposites | |
Kim et al. | Effect of the surface modification of silica nanoparticles on the viscosity and mechanical properties of silica/Epoxy nanocomposites |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |