一种高度交联环氧模塑化合物和铜基片之间拉伸强度的分子动力学研究资料.docVIP

一种高度交联环氧模塑化合物和铜基片之间 拉伸强度的分子动力学研究 Shaorui Yang a, Feng Gao b, Jianmin Qu a,b,* a Department of Mechanical Engineering, Northwestern University, 2145 Sheridan Rd, Evanston, IL 60208, USA b Department of Civil and Environmental Engineering, Northwestern University, 2145 Sheridan Rd, Evanston, IL 60208, USA 摘要：本文是基于经典的分子动力学模拟，以了解纯拉伸环氧/铜双材料垂直于界面的变形和破坏行为数值模拟研究。考虑到环氧树脂是一种高度交联的环氧苯酚酚醛清漆，以及铜基板是标准面心立方单晶111)面作为环氧/铜界面。应力位移/应变曲线，了解双材料行为和预测环氧树脂/铜界面的拉伸强度。，界面破坏是环氧原子从铜基底脱离脆性引起的，且所述界面抗拉强度破坏强度到达之前卸载和重新加载几乎不受影响。温度环氧交联密度环氧官能团的影响也进行了研究。这项研究的结果环氧树脂/铜双材料界面的变形和破坏行为的机制提供了显著见解。 1 引言 3.53×3.53 nm2在x和y方向有数千个原子），在他们的许多研究中，环氧网络除了原子非常接近的界面之外都假定为刚性。正如在我们的研究中看到的那样，由于环氧网络中原子长时间的相互作用，这些假设可能不适合用于获得界面的强度或韧性。 为了调查高度交联的环氧树脂/铜界面的变形和破坏行为，我们对环氧/铜双材料进行了大规模的分子动力学模拟。环氧/铜双材料是由附着在铜衬底上的高度交联的环氧树脂成型化合物。环氧树脂的网络结构是通过在其中共价键的形成和该系统的MD平衡被迭代地执行一个模拟的交联方法构成。交联的分子网络涂布在Cu基板上形成环氧/铜双材料。通过施加速度到衬底上的分子动力学模拟中，得到正应力-位移/应变关系。对模拟单元尺寸进行研究，以确保结果是独立的模拟单元尺寸。除了应力-应变关系，我们还研究了循环载荷、温度、转化率和环氧官能单体对应力-应变关系的影响。本文的其余部分安排如下：在第二节中，介绍原子模型和加载设置。在第三节中讨论仿真结果和各种结构并加载参数的影响。最后，一些结论的言论总结在第四节。 这一节中，我们注意到，此处获得的环氧化合物建模的建模方法并不局限于当前的工作，而是可以扩展到其他各种环氧体系。因此，从这个研究结果提供指引，设计聚合物/金属界面的电子封装产业。 2 模拟方法 2.1 材料体系和MD模拟单元 图1 （a） 图2 环氧/铜双材料示意图 2.2 原子间相互作用 (1) 该PCFF既包括化合价和非键关系。化合价关系包括键、角、二面角，以及它们之间的耦合作用。非键关系库仑力包括部分带电原子与Lennard-Jones（LJ）9-6之间的范德华力。在本文所提出的示例中，LJ势和库仑能量的现实空间的贡献进行计算的截止距离为12.5?。库仑长期的互惠空间的贡献计算是通过粒子颗粒与颗粒目方法。 由于在铜基板的变形可以忽略不计，我们假设Cu基体是刚性的。此外，它是合理的假设[37]，该环氧树脂是通过范德华力（无共价键）键合到铜基板。因此，该聚合物的原子和Cu原子之间的相互作用可以通过描述LJ9-6潜在的铜和环氧树脂原子参数得到，如下第六电力混合法则[38]，其中，下标E和W分别代表环氧树脂和铜基板。 (2) 2.3 聚合物的 为了模拟聚合物的交联，边界条件被规定，使得细胞周期中的环氧基/铜界面在平面内模拟，并且牵引-自由的方向垂直于界面。在这些边界条件下，单体的混合物动态交联网络使用Material Studio软件脚本运行。在该交联过程中，化学键被一对相对的反应性的原子，即在树脂的结束环氧碳原子和固化剂中的羟基氧之间形成，只要它们是在规定的截止距离之内。毕竟这些活性原子键合，该系统达到平衡由最小化1000步的能量的和10000步NVT MD模拟在500K下释放所引起的新形成的键的内部应力，并把剩余的反应性原子对接近。此交联的平衡过程反复进行，直到得到所希望的转化率。为了加快这一进程，该键的形成截止距离由3增加到10?。在每个截止，五交联平衡循环进行。在达到所需的转化率后，将未反应的反应性位点中和氢原子和部分电荷，使整个系统的电荷