聚甲基丙烯酸甲酯与碳纳米管纳米复合材料玻璃化转变及其非线性力学行为的分子动力学模拟*

2021-11-19黄多辉万明杰杨俊升2

物理学报 2021年21期

黄多辉万明杰杨俊升2)†

1) (宜宾学院理学部,宜宾 644007)

2) (华南理工大学,华南软物质科学与技术高等研究院,分子科学与工程学院,广州 510640)

短纤维结构对聚合物材料玻璃化转变温度及非线性力学具有非常重要的影响.本文利用粗粒化分子动力学方法研究了碳纳米管(CNTs)含量对聚甲基丙烯酸甲酯(PMMA)玻璃化转变、扩散系数及非线性力学特性的影响.分子动力学模拟结果显示:短CNTs 纤维的加入的确会改变PMMA 体系的玻璃化转变温度,模拟结果与实验测量结果一致,而且随着CNTs 含量的增加其对应的玻璃化转变温度也会随着增加.进一步分析扩散特性发现,CNTs 加入PMMA 对于体系扩散特性的改变主要发生在玻璃化温度以上,玻璃化温度以下结构对应的扩散系数差异非常的小.聚合物材料在服役过程中难免要遭受应力-应变的作用,而且其结构对应的模量和韧性成反比.基于此,本文通过非平衡分子动力学探究了短CNTs 纤维添加PMMA 复合材料的非线性力学特性.模拟结果显示:随着CNTs 纤维的含量增加,其对应的屈服模量也在不断的增加,而且含有短CNTs纤维的体系还能够保持原来的韧性.因此,分子层面的理论研究策略可以为进一步的实验和加工提供理论指导.

1 引言

通过短纳米纤维增强聚合物材料的力学、热学及电学性能的研究最近几年引起了研究者广泛的关注[1-4].首先,短纤维增强聚合物复合材料不仅仅在力学性能方面得到了很大的提升,而且短纳米纤维材料可以使聚合物力学模量提升的同时其对应的韧性特征受到很小的损失.其次,短纳米纤维增强的聚合物纳米复合物材料相比于连续性纤维材料更容易加工[5],成本更低,更容易被市场接受.故而,在成本和性能之间的折衷使得短纤维聚合物复合材料成为材料领域优秀的替代品.经过前期调研,我们发现在不同类型的纤维中,随着碳纳米管(CNTs)制备技术的成熟以及其优异的力学及电学性能可以为高性能聚合物基复合材料的制备及产业化提供可能[4,6,7].

根据现有的实验,不难发现CNTs 纤维与聚合物分子复合形成的聚合物纳米材料的确可以提高其对应的力学特性[8,9],这其中CNTs 的长度、取向和含量等参数都起着非常重要的作用.然而,由于当前实验表征技术在分子层面表征的困难,实验上还不能很好地描述CNTs 和聚合物之间微观相互作用机理.如何解释其中的微观相互作用机理弥补实验上的不足,对于CNT 与聚合物复合物材料的制备至关重要.相比于实验技术,分子动力学(MD)模拟可以在纳米尺度揭示实验一些现象的详细信息[10-13],如链构象、结构堆积信息及非平衡条件下的结构与力学性质的关系等.MD 模拟也有助于解释一些实验现象并且为纳米复合材料的设计开辟一条新的策略.因此,MD 模拟是理解纳米粒子增强聚合物纳米复合材料力学行为非常有利的手段.

尽管MD 模拟已广泛应用于增强聚合物纳米复合材料的建模,但对于一个全原子模型,其模拟所需要的计算工作量完全限制了其在有限时间尺度上对分子系统的适用性.很多时候只能通过单根CNTs 与聚合物的界面相互作用来探究其对应的微观机理[14-16],但是对于实际聚合物纳米复合物体系,纳米管是无规则分布在聚合物熔体中,在某种意义上来说小的全原子体系并不能体现出材料的宏观特性.这些缺陷也阻碍了MD 模拟研究纤维含量、尺寸和取向对增强聚合物纳米复合材料力学行为的影响.为了克服这些困难和限制,粗粒化(CG)模型无疑是最优的选择.当然,基于这些问题,文献[17-19]中针对石墨和CNTs 已经发展了一系列粗粒化MD 方法,并且这些方法的可靠性也得到了进一步的验证.目前,基于CG 模型大部分的研究主要还是在探究其结构对应的断裂行为,或者模型可靠性方面的验证性工作,对于聚合物纳米复合材料的热力学转变及结构和非线性力学关系的研究还非常的缺乏.

针对以上问题,本文选取短CNTs 纤维的CG结构模型和PMMA 聚合物分子CG 模型构成的复合材料模型为研究对象.基于CG 的分子力场来探究随机分布的短CNTs 纤维增强聚合物基体的玻璃化转变行为及非线性力学行为,并且给出不同CNTs 纤维含量下PMMA 玻璃化转变、扩散系数及相关拉伸和取向结构与非线性力学的关系.

2 模型与计算方法

本文使用了前人最近开发的CG 模型.基于CG模型,构建了PMMA/CNTs 聚合物纳米复合材料的体系.在该模型中,将每个甲基丙烯酸甲酯(C5O2H8)单体作为如图1(a)所示的原子质量为100.12 amu 的珠子,选择珠的中心作为单体的质量中心.CG 模型也用于建模(5,5) CNT,将5 个原子环映射成原子质量为600.55 amu 的CG 珠子(见图1(b)).基于以上的单体模型,本文构建了3 种体系,分别命名为PMMA,PMMA/10CNTs和PMMA/20CNTs.其中,PMMA 体系结构每根链含有200 个珠子,体系中包含200 根链(40000个珠子).PMMA/10CNTs 和PMMA/20CNTs 体系分别包含10 根和20 根CNTs 纤维,每根CNT纤维包含30 个珠子,并且CNT 在最开始的时候是无规分布在PMMA 熔体结构中,图2 给出了PMMA 及其CNT 在PMMA 中的分布图.

图1(a) 1 个PMMA 聚合物链的3 个单体及其由3 个珠子组成的CG 模型;(b) 1 个(5,5)带有10 个碳原子环的CNT 及其由3 个珠子组成的CG 模型Fig.1.(a) Three monomers of a PMMA polymer chain and its CG model made of three beads;(b) a (5,5) CNT with 10 rings of carbon atoms and its CG model made of three beads.

图2 PMMA 及PMMA/CNTs 结构示意图Fig.2.Schematic diagram of PMMA and PMMA/CNTs.

描述体系的粗粒化力场主要由非键相互作用(EvdW)、键合(Eb)及键角(Ea)相互作用构成,具体形式如下所示[20]:

其中,D0和r0分别为与平衡井深和平衡距离相关的Lennard-Jones 参数,在计算范德瓦耳斯(vdW)相互作用时,截断半径为1.25 nm.kd和d0分别为键长弹簧常数和平衡键距;kθ和θ0分别为键角的弹簧常数和平衡键角;力场参数如表1 所列.本文中所有的模拟都是基于Lammps 开源并行MD 软件,在宜宾学院曙光超级计算机上进行.为了保证体系的合理性,将所研究的体系在高温T=650 K下进行20 ns 的NPT 系综时间弛豫,压强保持1 atm (1 atm=1.01×105Pa).动力学弛豫和非平衡动力学过程中,本文选取的时间步长都是10 fs.

表1 PMMA 及PMMA-CNT 粗粒化粒子之间相互作用力场参数Table 1.Parameters of the CG force field for PMMA and CNT beads.

3 结果与讨论

3.1 PMMA 及PMMA/CNTs 聚合物纳米复合材料玻璃化转变特性

玻璃化转变温度对于聚合物纳米复合材料体系的应用起着至关重要的作用.图3 中通过体系的比容变化给出了PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系对应的玻璃化转变温度.为了保证模拟的合理性,降温过程中,本文选取的降温速度为20 K/10 ns.在玻璃化转变温度(Tg)以上或者以下,比容基本都是趋于线性变化,通过交叉点可以确定体系对应的玻璃化转变温度.PMMA体系对应的玻璃化转变温度点约为434 K,与最近的实验及全原子模拟值基本相符[16,21,22],进一步证明本文中选取模型及方法的可靠性.通过图3 不难发现,随着CNTs 纤维含量的增加,玻璃化转变温度随之增加,两个体系对应的玻璃化转变温度为454 K 和475 K.玻璃化温度的提升可以进一步拓展PMMA 结构体系在极端环境的使用.根据自由体系理论,CNTs 的加入,限制了PMMA 体系结构的运动能力,这一点会在接下来作进一步阐述.

图3 PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系对应的比容随着温度的变化Fig.3.Specific volumes of PMMA,PMMA/10CNTs and PMMA/20CNTs under the different temperatures.

CNTs 纤维加入PMMA 之后,它们之间的相互作用会直接影响PMMA 的运动能力,进而会影响其玻璃化温度的转变.最近的实验数据也指出,如果分子自由空间未被占据,对应的分子的运动能力也比较强.例如,在玻璃化温度以下,由于分子没有足够的自由空间,分子就无法运动和转动.但是目前,针对这一点实验上很难从分子层面解释其对应的构效关系和微观相互作用机理.为了进一步在分子层面解释清楚玻璃化温度提升的原因,分子层面上测定PMMA 及PMMA/CNTs 体系的运动及松弛动力学过程就尤为重要.基于平衡的MD方法,给出了所建3 个体系中PMMA 分子链整体的均方位移(MSD),根据MSD 的变化,计算了其对应的扩散系数.根据爱因斯坦方程,扩散系数与MSD 的关系如下[13]:

其中,r(t0)是粒子的初始坐标,r(t+t0)是粒子在时刻t的坐标.通过MSD (图4(a),PMMA/10CNTs与PMMA/20CNTs 对应的MSD 图已略去)与扩散系数满足的关系(4)式,图4(b)给出了PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系中PMMA 分子链整体运动对应的扩散系数随着温度的变化.从图4(b)不难发现,很明显CNTs 纤维加入之后PMMA 的扩散系数在不同的温度下会发生明显的降低.而且各体系都是以玻璃化温度为转变点发生突然的增长.这一点也印证了实验上的自由体积假说理论.

图4(a) PMMA 体系对应的MSD;(b) PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系不同温度下对应的扩散系数Fig.4.(a) Evolution of mean square displacement (MSD) of PMMA system;(b) the diffusion coefficient of PMMA,PMMA/10CNTs and PMMA/20CNTs under the different temperatures.

3.2 PMMA/CNTs 聚合物纳米复合材料非线性力学特性

聚合物材料在服役过程中难免要遭受外力场的作用,深刻认识外力作用下材料的结构和性能关系不仅具有重要的学术价值,也是决定材料能否长期有效服役的关键因素.聚合物纳米复合材料凝聚态结构是由大量的高分子依靠分子内和分子间的范德瓦耳斯相互作用凝聚而成,而且材料最后的拉伸性能是这些微观结构相互作用的结果.因此,非常有必要探究聚合物纳米复合材料结构与其非线性力学关系.图5(a)中,我们施加了一个恒定的拉伸速率1×109s—1来探究PMMA,PMMA/10CNTs及PMMA/20CNTs 结构体系的应力-应变过程,温度为室温(300 K).随着应变的增加,应力会经历典型的弹性形变、应力屈服、应力软化、应变硬化及断裂区.进一步观察,不难发现随着CNTs 的含量不同,PMMA 结构对应的力学性能会有很大的差别.随着CNT 的含量增加,其对应的屈服模量也在不断的增加.但是3 个体系的断裂伸长率并没有发生太大的改变,PMMA/10CNTs 相比于PMMA 体系稍微有所增加,但是随着含量的增加PMMA/20CNTs 相比于PMMA 就有所减小,总体而言断裂应变的变化都不大.

图5 PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系在300 K 下对应的应力-应变(a)、拉伸比(b)和取向参数(c)变化曲线Fig.5.Evolutions of stress-strain curves (a),stretch ratio (b) and orientational parameter (c) of PMMA,PMMA/10CNTs and PMMA/20CNTs under the temperature of 300 K.

聚合物纳米复合材料的力学特质主要由链构象对应的链结构微观拉伸比和链沿着拉伸方向的取向参数决定.图5(b)和图5(c)分别给出了PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系中PMMA 分子链对应的链微观拉伸比和末端矢量取向参数.采用λ=(L-L0)/L0计算链结构在宏观拉伸过程对应的微观拉伸比,其中,L0和L分别表示链原始末端距和受到拉伸以后的末端距,采用公式计算取向参数,其中,θ代表末端矢量与拉伸方向的夹角.不难发现,最开始阶段CNTs含量越多,PMMA 分子拉伸比增加得越快,这也是为什么在弹性形变区CNTs 含量越高,其结构对应的模量越高.但在3 个应变之后,PMMA/10CNTs体系对应更高的拉伸比和断裂应变.进一步分析取向参数,不难发现同样是PMMA/10CNTs 体系中PMMA 分子取向参数最大.根据拉伸比和取向参数不难得出,应力-应变曲线主要是分子链拉伸和取向协同作用的结果.

除此之外,还探究了不同温度对于PMMA,PMMA/10CNTs 及PMMA/20CNTs 结构体系的应力-应变过程的影响.图6(a)—(c)给出了PMMA,PMMA/10CNTs 及PMMA/20CNTs 体系在不同温度下对应的应力-应变曲线.不难发现,玻璃化温度以下3 个结构体系相比于玻璃化温度以上都发生了明显的应力软化现象及应变硬化现象.而且,随着CNTs 纤维含量的增加,硬化和屈服现象也越发得明显.为了进一步探究相同过冷度下,图6(d)给出了3 个体系在过冷度(ΔT=Tg– T=135 K)下的应力-应变曲线.不难发现,相同过冷度下随着CNTs 含量的增加,屈服强度明显的增加.但是PMMA/20CNTs 体系的断裂应变也会随着CNTs含量的增加而减小.

图6 PMMA (a),PMMA/10CNTs (b)及PMMA/20CNTs (c)结构体系在不同温度下对应的应力-应变曲线及在相同过冷度(ΔT= 135 K)下3 个体系结构对应的应力-应变曲线(d)Fig.6.E of stress-strain curves of PMMA (a),PMMA/10CNTs (b) and PMMA/20CNTs (c) under the different temperatures and the same cooling depth (ΔT=135 K)(d).

为了进一步说明PMMA,PMMA/10CNTs及PMMA/20CNTs 结构体系在应力-应变过程的结构失效机制,图7 和图8 给出了PMMA 及PMMA/10CNTs 结构的分子构象演化及其对应的原子维诺体积的变化图.从局部结构维诺体积变化发现,PMMA 及PMMA/10CNTs 体系在屈服、软化过程及应变硬化过程并没有发生结构的断裂,断裂发生在2 次屈服软化点附近.