CN111540413A - 一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法 - Google Patents
一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法 Download PDFInfo
- Publication number
- CN111540413A CN111540413A CN202010298185.4A CN202010298185A CN111540413A CN 111540413 A CN111540413 A CN 111540413A CN 202010298185 A CN202010298185 A CN 202010298185A CN 111540413 A CN111540413 A CN 111540413A
- Authority
- CN
- China
- Prior art keywords
- model
- epoxy resin
- sio
- thermal resistance
- interface
- 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
Links
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N Silicium dioxide Chemical class O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 title claims abstract description 102
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000003822 epoxy resin Substances 0.000 claims abstract description 55
- 239000000377 silicon dioxide Substances 0.000 claims abstract description 55
- 229920000647 polyepoxide Polymers 0.000 claims abstract description 54
- 239000002131 composite material Substances 0.000 claims abstract description 36
- 238000004132 cross linking Methods 0.000 claims abstract description 35
- 229910052681 coesite Inorganic materials 0.000 claims abstract description 29
- 229910052906 cristobalite Inorganic materials 0.000 claims abstract description 29
- 235000012239 silicon dioxide Nutrition 0.000 claims abstract description 29
- 229910052682 stishovite Inorganic materials 0.000 claims abstract description 29
- 229910052905 tridymite Inorganic materials 0.000 claims abstract description 29
- 238000004088 simulation Methods 0.000 claims abstract description 28
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000013215 result calculation Methods 0.000 claims abstract 2
- 238000005457 optimization Methods 0.000 claims description 40
- VILCJCGEZXAXTO-UHFFFAOYSA-N 2,2,2-tetramine Chemical compound NCCNCCNCCN VILCJCGEZXAXTO-UHFFFAOYSA-N 0.000 claims description 31
- 238000004422 calculation algorithm Methods 0.000 claims description 20
- 239000003795 chemical substances by application Substances 0.000 claims description 20
- 239000002245 particle Substances 0.000 claims description 19
- 229960001124 trientine Drugs 0.000 claims description 18
- 230000000737 periodic effect Effects 0.000 claims description 16
- 239000013078 crystal Substances 0.000 claims description 15
- 239000004593 Epoxy Substances 0.000 claims description 14
- 239000004841 bisphenol A epoxy resin Substances 0.000 claims description 13
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical group [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims description 10
- IISBACLAFKSPIT-UHFFFAOYSA-N bisphenol A Chemical compound C=1C=C(O)C=CC=1C(C)(C)C1=CC=C(O)C=C1 IISBACLAFKSPIT-UHFFFAOYSA-N 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 125000004430 oxygen atom Chemical group O* 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 8
- 238000005094 computer simulation Methods 0.000 claims description 5
- 238000007306 functionalization reaction Methods 0.000 claims description 5
- 125000004435 hydrogen atom Chemical group [H]* 0.000 claims description 5
- 125000002887 hydroxy group Chemical group [H]O* 0.000 claims description 5
- 230000003647 oxidation Effects 0.000 claims description 5
- 238000007254 oxidation reaction Methods 0.000 claims description 5
- 229910052710 silicon Inorganic materials 0.000 claims description 5
- 239000013598 vector Substances 0.000 claims description 5
- 230000033444 hydroxylation Effects 0.000 claims description 3
- 238000005805 hydroxylation reaction Methods 0.000 claims description 3
- 239000000178 monomer Substances 0.000 claims description 3
- 210000004027 cell Anatomy 0.000 claims 6
- 210000002858 crystal cell Anatomy 0.000 claims 3
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims 1
- 238000005520 cutting process Methods 0.000 claims 1
- 239000001301 oxygen Substances 0.000 claims 1
- 229910052760 oxygen Inorganic materials 0.000 claims 1
- 238000010008 shearing Methods 0.000 claims 1
- 239000000463 material Substances 0.000 abstract description 10
- 238000002474 experimental method Methods 0.000 abstract description 4
- 238000012546 transfer Methods 0.000 abstract description 4
- 230000007246 mechanism Effects 0.000 abstract description 3
- 238000004364 calculation method Methods 0.000 description 30
- 229910004298 SiO 2 Inorganic materials 0.000 description 26
- 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 13
- 230000004907 flux Effects 0.000 description 7
- 229920000642 polymer Polymers 0.000 description 7
- 125000004429 atom Chemical group 0.000 description 6
- 239000011231 conductive filler Substances 0.000 description 5
- 238000000329 molecular dynamics simulation Methods 0.000 description 5
- 239000010703 silicon Substances 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000000137 annealing Methods 0.000 description 2
- 238000005311 autocorrelation function Methods 0.000 description 2
- 229910052799 carbon Inorganic materials 0.000 description 2
- 125000004432 carbon atom Chemical group C* 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 239000000945 filler Substances 0.000 description 2
- 230000003993 interaction Effects 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 125000004433 nitrogen atom Chemical group N* 0.000 description 2
- 239000002861 polymer material Substances 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000010292 electrical insulation Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 125000000524 functional group Chemical group 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000010438 heat treatment Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 229920005989 resin Polymers 0.000 description 1
- 239000011347 resin Substances 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- 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)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，包括以下步骤：1)建立SiO2超晶胞模型；2)构建环氧树脂晶胞模型；3)构建复合材料层模型；4)层模型的交联过程；5)动力学平衡；6)模拟结果计算。本发明属于热界面复合材料界面热阻预测技术领域，具体是一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，有助于解决目前通过理论和实验方法不能准确分析热界面材料传热机理的问题。
Description
技术领域
本发明属于热界面复合材料界面热阻预测技术领域，具体是指一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法。
背景技术
热界面材料可以分为很多种，但其中的聚合物基热界面材料由于具有电绝缘性能优异，质轻，易加工成型和低成本等优点越来越受到重视。然而由于大多数聚合物材料的热导率仅为0.1～0.3W/(m·K)，导热能力较差，仅仅依赖聚合物自身的导热性能，难以满足电子元器件的散热要求，因此为了提高聚合物材料的导热性能通常是在聚合物基体材料中添加导热填料。而导热填料与高分子聚合物的界面相容性较差，导热填料粒子的表面不易被聚合物基体润湿，加之两者表面张力的差异，使得导热填料颗粒在高分子聚合物中难以充分有效地分散。因此，造成导热填料颗粒与聚合物基体的界面间有一定的空气间隙(空气热导率只有0.023W/(m·K)，从而使得复合材料的界面热阻较大。目前，主要是通过界面改性的方法减小热阻。热界面材料的界面热阻是微尺度传热领域的热点问题，由于它反映的是界面声子的传输性质，因此直接影响材料的热传导性能，从而对微/纳米器件的设计和热优化产生影响。
目前，通过微尺度传热理论分析和实验方法尚未对热界面材料界面热阻的内部机制和变化规律建立完善的理论体系和研究方法，从实验研究的角度而言，界面热阻难以直接测量。近几年发展起来的分子动力学模拟是理论和实验有效的补充手段，它通过原子间相互作用势，求出每一个原子所受的力，在选定的时间步长、边界条件、初始位置和初始速度下，对有限数目的分子(原子)建立其牛顿动力学方程组，用数值方法求解得到这些原子经典运动轨迹和速度，然后对足够长时间的结果求统计平均，从而得到所需要的宏观物理量和力学量。
据文献报道，计算界面热阻的非平衡态分子动力学方法主要有以下两种模型：一种是周期性的反向非平衡分子动力学(RNEMD)模型，模型一端和模型中部的粒子通过交换速度实现外加热流，周期性边界条件应用于x，y，z三个方向，模型中产生的局部温度具有对称性；另一种是绝热壁形式的非平衡分子动力学(NEMD)模型，对模拟体系外加热流或温度梯度，周期性边界条件只应用于模型的y，z方向上，分别在x方向两端设置绝热壁，模型中产生的局部温度不具有对称性。周期性的RNEMD模型由于采用的是对称体系，所以该方法的计算精度高，但计算量也大，而绝热壁的NEMD模型计算量小，通过固定边界原子位置来实现绝热。本发明中的模型如图2所示，不同之处在于：通过在z方向上添加真空层来实现z方向上的绝热条件，通过交换体系两端的粒子速度来实现外加热流，计算量可进一步减小，这样，既能研究填料与基体界面热阻的大小，又能研究填料与基体单侧的界面相互作用能，计算得到的界面热阻通过与界面能量的对比分析，可以保证仿真结果的可靠性。
发明内容
为解决上述现有难题，本发明为了解决目前通过理论和实验方法不能准确分析热界面材料传热机理的问题，提供一种计算热界面材料界面热阻的新方法。
本发明采用的技术方案如下：一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，包括以下步骤：
1)建立SiO2超晶胞模型：根据α-SiO2的空间群和晶格参数，建立SiO2的晶体模型，并切割晶体构建超晶胞模型，然后对超胞模型进行几何优化；
2)构建环氧树脂晶胞模型：构建双酚A型环氧树脂和三乙烯四胺固化剂的单体分子模型，建立45个双酚A型环氧树脂分子和15个三乙烯四胺固化剂分子组成的初始密度为0.6g/cm3的3维周期性晶胞，对晶胞进行几何优化；
3)构建复合材料层模型：建立SiO2和环氧树脂的层模型，选择优化算法对层模型进行几何优化；
4)层模型的交联过程：通过运行xink.pl交联脚本，对优化后层模型中的双酚A型环氧树脂和三乙烯四胺固化剂分子进行交联反应，得到环氧树脂分子链；
5)动力学平衡：对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟，使模型达到平衡；
6)模拟结果计算：通过运行TC.pl脚本，得到SiO2/环氧树脂层状模型的热流密度和温度分布，计算得到界面热阻。
进一步地，步骤1)根据α-SiO2的空间群P3221和晶格参数：
a＝b＝4.910，c＝5.402；
在分数坐标下，硅原子和氧原子的位置参数分别为：
Si：a＝0.480781，b＝0.480781，c＝0；
O：a＝0.150179，b＝0.414589，c＝0.116499；
得到SiO2的晶体模型，剪切(0 0 
)晶面，表面向量U(1 1 0)、V(1 -1 0)，建立7×4超胞，模型尺寸为
的SiO2超晶胞模型；对SiO2表面羟基化处理，选择表面硅原子上的断裂键并将其与羟基结合、将表面氧原子上的断裂键与氢原子结合以模拟真实的氧化过程；选择Smart优化算法对SiO2超晶胞模型进行几何优化。
进一步地，步骤2)构建双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子模型，然后构建一个由45个DGEBA分子和15个TETA分子组成的初始密度为0.6g/cm3的3维周期性晶胞，选定COMPASSII力场，模型精度设置为Medium；选择Smart优化算法对晶胞模型进行几何优化。
进一步地，步骤3)将上面得到的SiO2超晶胞模型设置为Layer 1，环氧树脂晶胞模型设置为Layer 2，并在Layer 2下设置
的真空层，其他设置默认，建立二氧化硅/环氧树脂复合材料的层模型；选择Smart优化算法对层模型进行几何优化。
进一步地，步骤4)在分子动力学软件中打开xink.pl交联脚本文件，定义反应位点，设置好交联参数，运行脚本文件进行交联模拟，得到环氧树脂分子链，然后将研究体系移到盒子中间。
进一步地，步骤5)对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟。
进一步地，步骤6)运行TC.pl脚本，对平衡后的模型进行界面热阻的计算；采用交换模型两端粒子速度实现外加热流的NEMD方法，周期性边界条件只应用于模型的x，y方向上，在模型的z方向上引入30埃的真空层作为边界条件，得到层状模型的热流密度和温度分布，基于传统的傅里叶定律计算不同功能化和不同接枝率的二氧化硅与环氧树脂之间的界面热阻。
采用上述方案本发明取得有益效果如下：本发明一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，可研究不同功能化和不同接枝率的二氧化硅与环氧树脂之间的界面热阻，同时计算了界面能量和振动能谱，对界面问题的研究更加深入细致，且模型原子数目减少，计算量减少，但有界面能量和振动能谱印证界面热阻的模拟结果，以保证模拟结果的准确性，这有效节省了计算成本，并通过仿真计算和理论研究复合材料的界面热阻，不仅可以减少实验成本，对于设计具备优异性能的热界面材料具有重要的研究意义和应用价值。
附图说明
图1为本发明二氧化硅/环氧树脂复合材料界面热阻的预测新方法的流程图；
图2为二氧化硅/环氧树脂复合材料界面热阻的预测新方法的仿真模型示意图；
图3为二氧化硅/环氧树脂复合材料界面热阻的预测新方法的温度分布图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
本发明采用的技术方案如下：一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，包括以下步骤：
1)建立SiO2超晶胞模型：根据α-SiO2的空间群和晶格参数，建立SiO2的晶体模型，并切割晶体构建超晶胞模型，然后对超胞模型进行几何优化；
2)构建环氧树脂晶胞模型：构建双酚A型环氧树脂和三乙烯四胺固化剂的单体分子模型，建立45个双酚A型环氧树脂分子和15个三乙烯四胺固化剂分子组成的初始密度为0.6g/cm3的3维周期性晶胞，对晶胞进行几何优化；
3)构建复合材料层模型：建立SiO2和环氧树脂的层模型，选择优化算法对层模型进行几何优化；
4)层模型的交联过程：通过运行xink.pl交联脚本，对优化后层模型中的双酚A型环氧树脂和三乙烯四胺固化剂分子进行交联反应，得到环氧树脂分子链；
5)动力学平衡：对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟，使模型达到平衡；
6)模拟结果计算：通过运行TC.pl脚本，得到SiO2/环氧树脂层状模型的热流密度和温度分布，计算得到界面热阻。
进一步地，步骤1)根据α-SiO2的空间群P3221和晶格参数：
a＝b＝4.910，c＝5.402；
在分数坐标下，硅原子和氧原子的位置参数分别为：
Si：a＝0.480781，b＝0.480781，c＝0；
O：a＝0.150179，b＝0.414589，c＝0.116499；
得到SiO2的晶体模型，剪切(0 0 
)晶面，表面向量U(1 1 0)、V(1 -1 0)，建立7×4超胞，模型尺寸为
的SiO2超晶胞模型；对SiO2表面羟基化处理，选择表面硅原子上的断裂键并将其与羟基结合、将表面氧原子上的断裂键与氢原子结合以模拟真实的氧化过程；选择Smart优化算法对SiO2超晶胞模型进行几何优化。。
步骤2)构建双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子模型，然后构建一个由45个DGEBA分子和15个TETA分子组成的初始密度为0.6g/cm3的3维周期性晶胞，选定COMPASSII力场，模型精度设置为Medium；选择Smart优化算法对晶胞模型进行几何优化。
步骤3)将上面得到的SiO2超晶胞模型设置为Layer 1，环氧树脂晶胞模型设置为Layer 2，并在Layer 2下设置
的真空层，其他设置默认，建立二氧化硅/环氧树脂复合材料的层模型；选择Smart优化算法对层模型进行几何优化。
步骤4)在分子动力学软件中打开xink.pl交联脚本文件，定义反应位点，设置好交联参数，运行脚本文件进行交联模拟，得到环氧树脂分子链，然后将研究体系移到盒子中间。
步骤5)对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟。
步骤6)运行TC.pl脚本，对平衡后的模型进行界面热阻的计算；采用交换模型两端粒子速度实现外加热流的NEMD方法，周期性边界条件只应用于模型的x，y方向上，在模型的z方向上引入30埃的真空层作为边界条件，得到层状模型的热流密度和温度分布，基于传统的傅里叶定律计算不同功能化和不同接枝率的二氧化硅与环氧树脂之间的界面热阻。
具体使用时，建立模型：采用分子动力学软件建立不同功能化和不同接枝率的二氧化硅与环氧树脂复合材料模型，运行xink.pl交联脚本对双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子进行交联反应，得到环氧树脂分子链，然后进行动力学平衡，最后通过运行TC.pl脚本计算界面热阻，计算二氧化硅/环氧树脂复合材料界面热阻的仿真模型示意图如图2；仿真计算：对平衡后的模型进行界面热阻的计算；采用交换模型两端粒子速度实现外加热流的NEMD方法，但边界条件不同于原来的边界条件，周期性边界条件只应用于模型的x、y方向，而在模型的z方向引入一个
的真空层作为边界条件；为了获得温度分布，首先将整个复合材料体系，而不是晶格，划分成20个等份，从K1到K20标记，沿着热流方向(模拟中的z轴方向)，如图2所示；整个体系末端的K1段称为“热”端，另一端的K20段称为“冷”端；然后使冷端动能最大的粒子和热端动能最小的粒子之间的能量和动量交换，就像是两个选定的粒子之间会有一个假设的弹性碰撞一样；这样的交换不仅可以保持整个系统的总能量和总动量不变，而且交换后，热端温度升高，冷端温度降低；
动量守恒公式如下：
能量守恒公式如下：
从公式(1)和(2)，可以得到：
值得注意的是，当mc＝mh时，冷、热端粒子的新速度分别表示为：
其中下标c和h分别表示冷端和热端中的粒子，mc和mh分别表示在冷端和热端中被选择的粒子的质量，
和
分别表示交换前后冷端粒子的速度，
和
分别表示交换前后热端粒子的速度；
通过计算单位时间单位面积内的能量交换，可以得到热流密度，公式如下：
其中〈Jz(t)〉是z方向的热流密度，A是横截面积，t是粒子速度进行交换的时间间隔；
根据理想气体气动理论的温度关系公式：
得到整个仿真体系的局部温度计算公式：
其中，Tk表示第k层的温度，nk是第k层的粒子数，kB是玻尔兹曼常数；界面热阻计算：根据上述公式(7)和(9)分别得到的热流密度和局部温度可计算界面热阻：
其中Rint表示界面热阻，ΔT是不同位置的温度差。
实施列1
未经修饰的二氧化硅/环氧树脂复合材料界面热阻的计算过程，具体步骤如下：
建立SiO2超晶胞模型：根据α-SiO2的空间群P3221(第154号)和晶格参数：
a＝b＝4.910，c＝5.402；
在分数坐标下，硅原子和氧原子的位置参数分别为：
Si：a＝0.480781，b＝0.480781，c＝0；
O：a＝0.150179，b＝0.414589，c＝0.116499；
得到SiO2的晶体模型，剪切(0 0 
)晶面，表面向量U(1 1 0)、V(1 -1 0)，建立7×4超胞，模型尺寸为
的SiO2超晶胞模型，然后对SiO2表面羟基化处理，即选择表面硅原子上的断裂键并将其与羟基结合、将表面氧原子上的断裂键与氢原子结合以模拟真实的氧化过程，选择Smart优化算法对SiO2超晶胞模型进行几何优化，在分子动力学软件中选择Geometry Optimization，模型精度设置为Medium，优化的算法选择Smart，分子动力学的力场选择COMPASSII，Summation method选择Ewald，点击Run进行优化计算。随着计算步长的增加总的能量值趋于稳定时，得到收敛曲线，即模型达到平衡稳定状态；
构建环氧树脂晶胞模型：构建双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子模型，然后构建一个由45个DGEBA分子和15个TETA分子组成的初始密度为0.6g/cm3的3维周期性晶胞，其中COMPASSII力场被选中，模型精度设置为Medium；选择Smart优化算法对晶胞模型进行几何优化，随着计算步长的增加,总的能量值趋于稳定时，得到收敛曲线，即模型达到平衡稳定状态；
构建复合材料层模型：将上面得到的SiO2超晶胞模型设置为Layer 1，环氧树脂晶胞模型设置为Layer 2，并在Layer 2下设置
的真空层，其他设置默认，建立二氧化硅/环氧树脂复合材料的层模型；然后选择Smart优化算法对层模型进行几何优化，随着计算步长的增加,总的能量值趋于稳定时，得到收敛曲线，即模型达到平衡稳定状态；
层模型的交联过程：在分子动力学软件中打开xink.pl交联脚本文件，在优化后层模型中的双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子上定义反应位点，即对DGEBA分子上的末端碳原子和TETA分子上的氮原子进行重新命名，分别命名为R1和R2；设置的初始截断距离为
最大截断距离为
步长为
目标交联度为85％，反应温度为300K；然后运行脚本文件进行交联模拟，得到环氧树脂分子链，最后将研究体系移到盒子中间；
动力学平衡：对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟，在分子动力学软件中分别进行Anneal和Dynamics计算；NVT退火温度范围为300K～600K，时间步长为1fs，总模拟时长为200ps，力场为COMPASSII，控温的方法为Andersen；然后选择正则系综(NVT系综)对体系进行1次平衡，随后选择微正则系综(NVE系综)再次进行平衡，2次模拟均设置温度为298K，步长分别为1fs和0.1fs，模拟时间长分别为2000ps和10ps，力场为COMPASSII，控温的方法为Andersen。
模拟结果计算：Dynamics计算完成后选择Interaction_Energy.pl脚本计算界面能量；选择速度自相关函数收集原子轨迹，并对收集的轨迹数据进行分析，研究界面声子的匹配性；运行TC.pl脚本，对平衡后的模型进行界面热阻的计算；
未经修饰的二氧化硅/环氧树脂复合材料界面模型，计算得到的界面能量值为-284.81kcal/mol，计算振动能谱发现界面声子匹配性较差，界面热阻值约为0.625×10-8m2KW-1。
实施例2
KH550修饰的二氧化硅/环氧树脂复合材料界面热阻的计算过程，具体步骤如下：
建立SiO2超晶胞模型：根据α-SiO2的空间群P3221(第154号)和晶格参数：
a＝b＝4.910，c＝5.402。
在分数坐标下，硅原子和氧原子的位置参数分别为：
Si：a＝0.480781，b＝0.480781，c＝0；
O：a＝0.150179，b＝0.414589，c＝0.116499；
得到SiO2的晶体模型，剪切(0 0 
)晶面，表面向量U(1 1 0)、V(1 -1 0)，建立7×4超胞，模型尺寸为
的SiO2超晶胞模型；然后对SiO2表面羟基化处理，即选择表面硅原子上的断裂键并将其与羟基结合、将表面氧原子上的断裂键与氢原子结合以模拟真实的氧化过程，设定接枝率5.4％，手动接枝KH550官能团到SiO2表面；选择Smart优化算法对SiO2超晶胞模型进行几何优化，在分子动力学软件中选择GeometryOptimization，模型精度设置为Medium，优化的算法选择Smart，分子动力学的力场选择COMPASSII，Summation method选择Ewald，点击Run进行优化计算。随着计算步长的增加，总的能量值趋于稳定时，得到收敛曲线，即模型达到平衡稳定状态；
构建环氧树脂晶胞模型：构建双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子模型，然后构建一个由45个DGEBA分子和15个TETA分子组成的初始密度为0.6g/cm3的3维周期性晶胞，其中COMPASSII力场被选中，模型精度设置为Medium；选择Smart优化算法对晶胞模型进行几何优化，随着计算步长的增加，总的能量值趋于稳定时，得到收敛曲线，即模型达到平衡稳定状态；
构建复合材料层模型：将上面得到的SiO2超晶胞模型设置为Layer 1，环氧树脂晶胞模型设置为Layer 2，并在Layer 2下设置
的真空层，其他设置默认，建立二氧化硅/环氧树脂复合材料的层模型；然后选择Smart优化算法对层模型进行几何优化，随着计算步长的增加，总的能量值趋于稳定时，得到收敛曲线，即模型达到平衡稳定状态；
层模型的交联过程：在分子动力学软件中打开xink.pl交联脚本文件，在优化后层模型中的双酚A型环氧树脂(DGEBA)和三乙烯四胺(TETA)固化剂分子上定义反应位点，即对DGEBA分子上的末端碳原子和TETA分子上的氮原子进行重新命名，分别命名为R1和R2；设置的初始截断距离为
最大截断距离为
步长为
目标交联度为85％，反应温度为300K；然后运行脚本文件进行交联模拟，得到环氧树脂分子链，最后将研究体系移到盒子中间；
动力学平衡：对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟，在分子动力学软件中分别进行Anneal和Dynamics计算；NVT退火温度范围为300K～600K，时间步长为1fs，模拟时间长为200ps，力场为COMPASSII，控温的方法为Andersen；然后选择正则系综(NVT系综)对体系进行1次平衡，随后选择微正则系综(NVE系综)再次进行平衡，2次模拟均设置温度为298K，步长分别为1fs和0.1fs，模拟时间长分别为2000ps和10ps，力场为COMPASSII，控温的方法为Andersen；
模拟结果计算：Dynamics计算完成后选择Interaction_Energy.pl脚本计算界面能量；选择速度自相关函数收集原子轨迹，并对收集的轨迹数据进行分析，研究界面声子的匹配性；运行TC.pl脚本，对平衡后的模型进行界面热阻的计算；
KH550修饰的二氧化硅/环氧树脂复合材料界面模型，计算得到的界面能量值为-353.03kcal/mol，计算振动能谱发现界面声子匹配性较好，界面热阻值约为0.505×10-8m2KW-1。
以上所述仅为本发明的实施例，并非因此限制本发明的专利范围，凡是利用本发明说明书内容所作的等效结构或等效流程变换，或直接或间接运用在其它相关的技术领域，均同理包括在本发明的专利保护范围内。
Claims (7)
1.一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，包括以下步骤：
1)建立SiO2超晶胞模型：根据α-SiO2的空间群和晶格参数，建立SiO2的晶体模型，并切割晶体构建超晶胞模型，然后对超胞模型进行几何优化；
2)构建环氧树脂晶胞模型：构建双酚A型环氧树脂和三乙烯四胺固化剂的单体分子模型，建立45个双酚A型环氧树脂分子和15个三乙烯四胺固化剂分子组成的初始密度为0.6g/cm3的3维周期性晶胞，对晶胞进行几何优化；
3)构建复合材料层模型：建立SiO2和环氧树脂的层模型，选择优化算法对层模型进行几何优化；
4)层模型的交联过程：通过运行xink.pl交联脚本，对优化后层模型中的双酚A型环氧树脂和三乙烯四胺固化剂分子进行交联反应，得到环氧树脂分子链；
5)动力学平衡：对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟，使模型达到平衡；
6)模拟结果计算：通过运行TC.pl脚本，得到SiO2/环氧树脂层状模型的热流密度和温度分布，计算得到界面热阻。
2.根据权利要求1所述的一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，步骤1)根据α-SiO2的空间群P3221和晶格参数：
a＝b＝4.910，c＝5.402；
在分数坐标下，硅原子和氧原子的位置参数分别为：
Si：a＝0.480781，b＝0.480781，c＝0；
O：a＝0.150179，b＝0.414589，c＝0.116499；
得到SiO2的晶体模型，剪切
晶面，表面向量U(1 1 0)、V(1 -1 0)，建立7×4超胞，模型尺寸为
的SiO2超晶胞模型；对SiO2表面羟基化处理，选择表面硅原子上的断裂键并将其与羟基结合、将表面氧原子上的断裂键与氢原子结合以模拟真实的氧化过程；选择Smart优化算法对SiO2超晶胞模型进行几何优化。
3.根据权利要求1所述的一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，步骤2)构建双酚A型环氧树脂和三乙烯四胺固化剂分子模型，然后构建一个由45个双酚A型环氧树脂和15个三乙烯四胺固化剂分子组成的初始密度为0.6g/cm3的3维周期性晶胞，选定COMPASSII力场，模型精度设置为Medium；选择Smart优化算法对晶胞模型进行几何优化。
4.根据权利要求1所述的一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，步骤3)将上面得到的SiO2超晶胞模型设置为Layer1，环氧树脂晶胞模型设置为Layer 2，并在Layer 2下设置
的真空层，其他设置默认，建立二氧化硅/环氧树脂复合材料的层模型；选择Smart优化算法对层模型进行几何优化。
5.根据权利要求1所述的一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，步骤4)在分子动力学软件中打开xink.pl交联脚本文件，定义反应位点，设置好交联参数，运行脚本文件进行交联模拟，得到环氧树脂分子链，然后将研究体系移到盒子中间。
6.根据权利要求1所述的一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，步骤5)对交联完成后的SiO2/环氧树脂层状模型进行动力学模拟。
7.根据权利要求1所述的一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法，其特征在于，步骤6)运行TC.pl脚本，对平衡后的模型进行界面热阻的计算；采用交换模型两端粒子速度实现外加热流的NEMD方法，周期性边界条件只应用于模型的x，y方向上，在模型的z方向上引入30埃的真空层作为边界条件，得到层状模型的热流密度和温度分布，基于传统的傅里叶定律计算不同功能化和不同接枝率的二氧化硅与环氧树脂之间的界面热阻。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010298185.4A CN111540413A (zh) | 2020-04-16 | 2020-04-16 | 一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010298185.4A CN111540413A (zh) | 2020-04-16 | 2020-04-16 | 一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN111540413A true CN111540413A (zh) | 2020-08-14 |
Family
ID=71974999
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010298185.4A Pending CN111540413A (zh) | 2020-04-16 | 2020-04-16 | 一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111540413A (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112164425A (zh) * | 2020-09-27 | 2021-01-01 | 齐鲁工业大学 | 一种基于Materials Studio的二氧化硅包覆氟化钙的模拟方法 |
| CN112908422A (zh) * | 2021-01-07 | 2021-06-04 | 国网内蒙古东部电力有限公司电力科学研究院 | 复合绝缘子芯棒环氧树脂材料的交联模型构建方法和老化仿真方法 |
| CN113030175A (zh) * | 2021-04-13 | 2021-06-25 | 中南大学 | 一种计算沥青混合料界面热阻和界面热导系数的方法 |
| CN113392554A (zh) * | 2021-06-15 | 2021-09-14 | 南京航空航天大学 | 一种热致相变膜调控刀具涂层应力的评估方法 |
| CN113470756A (zh) * | 2021-07-05 | 2021-10-01 | 中国科学院化学研究所 | 一种构建芳香环化合物与脂肪族化合物交联模型的方法及模型 |
| CN113990413A (zh) * | 2021-10-26 | 2022-01-28 | 西安交通大学 | 建立聚乙烯链结构与挤出加工性能关系的仿真方法及系统 |
| CN114974468A (zh) * | 2022-03-22 | 2022-08-30 | 武汉大学 | 金刚石基复合材料界面热阻的优化方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100038544A1 (en) * | 2008-08-14 | 2010-02-18 | The Boeing Company | Thermal Effect Measurement with Mid-Infrared Spectroscopy |
| CN107967403A (zh) * | 2017-12-28 | 2018-04-27 | 南京理工大学 | 石墨烯/环氧树脂复合热界面材料热导率的预测方法 |
-
2020
- 2020-04-16 CN CN202010298185.4A patent/CN111540413A/zh active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100038544A1 (en) * | 2008-08-14 | 2010-02-18 | The Boeing Company | Thermal Effect Measurement with Mid-Infrared Spectroscopy |
| CN107967403A (zh) * | 2017-12-28 | 2018-04-27 | 南京理工大学 | 石墨烯/环氧树脂复合热界面材料热导率的预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 张晓星，等: "改性SiO2 增强环氧树脂热力学性能的分子动力学模拟", 《高压电技术》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112164425A (zh) * | 2020-09-27 | 2021-01-01 | 齐鲁工业大学 | 一种基于Materials Studio的二氧化硅包覆氟化钙的模拟方法 |
| CN112164425B (zh) * | 2020-09-27 | 2023-03-28 | 齐鲁工业大学 | 一种基于Materials Studio的二氧化硅包覆氟化钙的模拟方法 |
| CN112908422B (zh) * | 2021-01-07 | 2022-07-08 | 国网内蒙古东部电力有限公司电力科学研究院 | 复合绝缘子芯棒环氧树脂材料的交联模型构建方法和老化仿真方法 |
| CN112908422A (zh) * | 2021-01-07 | 2021-06-04 | 国网内蒙古东部电力有限公司电力科学研究院 | 复合绝缘子芯棒环氧树脂材料的交联模型构建方法和老化仿真方法 |
| CN113030175A (zh) * | 2021-04-13 | 2021-06-25 | 中南大学 | 一种计算沥青混合料界面热阻和界面热导系数的方法 |
| CN113392554A (zh) * | 2021-06-15 | 2021-09-14 | 南京航空航天大学 | 一种热致相变膜调控刀具涂层应力的评估方法 |
| CN113392554B (zh) * | 2021-06-15 | 2024-08-09 | 南京航空航天大学 | 一种热致相变膜调控刀具涂层应力的评估方法 |
| CN113470756A (zh) * | 2021-07-05 | 2021-10-01 | 中国科学院化学研究所 | 一种构建芳香环化合物与脂肪族化合物交联模型的方法及模型 |
| CN113470756B (zh) * | 2021-07-05 | 2023-11-28 | 中国科学院化学研究所 | 芳香环化合物与脂肪族化合物交联模型及其构建方法 |
| CN113990413A (zh) * | 2021-10-26 | 2022-01-28 | 西安交通大学 | 建立聚乙烯链结构与挤出加工性能关系的仿真方法及系统 |
| CN113990413B (zh) * | 2021-10-26 | 2023-09-12 | 西安交通大学 | 建立聚乙烯链结构与挤出加工性能关系的仿真方法及系统 |
| CN114974468A (zh) * | 2022-03-22 | 2022-08-30 | 武汉大学 | 金刚石基复合材料界面热阻的优化方法 |
| CN114974468B (zh) * | 2022-03-22 | 2024-04-16 | 武汉大学 | 金刚石基复合材料界面热阻的优化方法 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111540413A (zh) | 一种二氧化硅/环氧树脂复合材料界面热阻的预测新方法 | |
| Artrith et al. | Constructing first-principles phase diagrams of amorphous LixSi using machine-learning-assisted sampling with an evolutionary algorithm | |
| Zhai et al. | Effective thermal conductivity of polymer composites: Theoretical models and simulation models | |
| Yuan et al. | Effects of surface functionalization on thermal and mechanical properties of graphene/polyethylene glycol composite phase change materials | |
| Yu et al. | Hard carbon nanotube sponges for highly efficient cooling via moisture absorption–desorption process | |
| Aslfattahi et al. | Efficiency enhancement of a solar dish collector operating with a novel soybean oil-based-MXene nanofluid and different cavity receivers | |
| CN110442887B (zh) | 一种研究重复聚酰亚胺体系热力学性质的模型构建方法 | |
| Stournara et al. | From ab initio calculations to multiscale design of Si/C core–shell particles for Li-ion anodes | |
| CN108314777B (zh) | 一种环氧/酸酐固化物交联分子模型构建及其玻璃化温度提取方法 | |
| Su et al. | Optimization of effective thermal conductivity of thermal interface materials based on the genetic algorithm-driven random thermal network model | |
| Peng et al. | Enhancement on thermal and mechanical properties of insulating paper cellulose modified by silane coupling agent grafted hBN | |
| CN107967403A (zh) | 石墨烯/环氧树脂复合热界面材料热导率的预测方法 | |
| Li et al. | Multi-particle molecular dynamics simulation: shell thickness effects on sintering process of Cu-Ag core-shell nanoparticles | |
| Liu et al. | Interfacial thermal conductance across hexagonal boron nitride & paraffin based thermal energy storage materials | |
| Shen et al. | Investigation of physical properties of epoxy-functionalized graphene nanoplatelets composite coatings on DC-GIL insulators by molecular dynamics simulation | |
| Wu et al. | Effect analysis on thermal management of power batteries utilizing a form-stable silicone grease/composite phase change material | |
| Pan et al. | Numerical Simulation of Manifold Microchannel Heat Sinks for Thermal Management in a Li‐Ion Battery | |
| Zhu et al. | Molecular dynamics simulation study of the solid polymer electrolyte that PEO grafted POSS | |
| Kasim et al. | Effect of thermal stratification on MHD free convection with heat and mass transfer over an unsteady stretching surface with heat source, Hall current and chemical reaction | |
| Wang et al. | Performance simulation of L-shaped heat pipe and air coupled cooling process for ternary lithium battery module | |
| Chen et al. | Unified deep learning network for enhanced accuracy in predicting thermal conductivity of bilayer graphene, hexagonal boron nitride, and their heterostructures | |
| Yuan et al. | Multi-objective optimization of thermal modules in high heat flux laptops | |
| Kumar et al. | Numerical and analytical modelling of effective thermal conductivity of multi-walled carbon nanotubes polymer nanocomposites including the effect of nanotube orientation and interfacial thermal resistance | |
| CN113257366A (zh) | 一种支化度对超支化聚酯改性的纳米复合环氧树脂导热性影响的分子动力学模拟方法 | |
| Zhao et al. | Reducing interfacial thermal resistance between polyethylene oxide-based solid-state polymer electrolyte and lithium anode by using IVA group two-dimensional materials: A molecular dynamics study |
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 | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20200814 |
|
| RJ01 | Rejection of invention patent application after publication |