碳纳米管表面缺陷对聚甲基丙烯酸甲酯基体弹性的影响

2020-06-04李云龙杨斌唐黎明李皓王照智刘冬

中南大学学报（自然科学版） 2020年4期

李云龙，杨斌，唐黎明，李皓，王照智，刘冬

(1. 汕头大学工学院，广东汕头，515063；2. 沈阳工业大学机械工程学院，辽宁沈阳，110870；3. 沈阳化工大学能源与动力工程学院，辽宁沈阳，110142)

碳纳米管(CNTs)因其良好的力学性能、热力学性能和表面吸附活性等优势，被视为理想的纳米补强剂，广泛地应用在高分子复合材料领域[1-3]。近年来，许多学者研究了CNTs增强高分子复合基体的力学性能。DEPLANCKE 等[4]探索CNTs 的预设分布对超高分子量聚乙烯(UHMWPE)纳米复合材料的力学性能的影响，发现CNTs/UHMWPE 复合材料的屈服应力和应变硬化特性可在混入较高的CNTs 浓度后得到改善；TARFAOUI 等[5]分别制备了不同体积分数的CNTs纳米织物复合材料，探索了不同体积分数下的CNTs对织物基复合材料基体弹性性能的影响。实验观察到加入的CNTs可有效地控制复合材料试件的损伤萌生和裂纹扩展。同时，织物基复合材料的弹性模量可提升2%以上。为深入掌握CNTs/高分子纳米复合材料的结构和力学特性，特别是在探索CNTs与高分子基体间的界面作用方面，基于连续力学的理论探索方法得到了长足发展。ZALAMEA等[6]利用连续介质力学中的2 种方法计算了多壁碳纳米管(MWCNT)结构中多个壳层之间的界面相互作用，并认为MWCNT 能有效地将应力从基质传递到其内部；GAO 等[7]提出了由CNTs 增强的聚合物复合材料的剪切滞后模型，并使用多尺度方法研究从聚合物基体到CNTs最内层的载荷转移以及影响载荷传递的因素；SPANOS 等[8]应用细观力学混合有限元方法研究了单壁碳纳米管(SWCNT)增强复合材料的应力传递，并详细分析了SWCNT体积分数、界面刚度和基体弹性模量对应力传递的影响。值得注意的是，基于宏观连续力学的理论模型均采取对CNTs 的无缺陷结构的连续体假设，但在实际的CNTs/高分子纳米复合材料的制备过程中，通常需要采取如温度、压力和化学处理等复杂的工艺流程处理。由于CNTs 自身的sp2杂化结构具有较好反应活性，目前CNTs表面的缺陷现象是无法避免的[9]。此时，传统的连续力学理论模型无法捕捉到CNTs 自身微观结构的改变，导致很难通过实验以及传统连续力学等手段预测与分析CNTs/高分子纳米复合基体准确的力学性能。分子动力学(MD)是一种有效的计算机模拟方法，用来探究不同环境下分子构象与形态以及分子之间相互作用过程。MD 模拟可以有效地捕捉CNTs 表面缺陷的微观信息，并同时计算其与高分子之间的相互作用以及高分子纳米复合材料的宏观力学性能。LU 等[10]利用了MD 理论，模拟研究了Stone-Wales 缺陷对单壁纳米管(SCNT)的力学性能的影响，发现SCNT的断裂总是从缺陷的位置开始发生，而对于无缺陷管，裂纹位置随机产生。对于给定的缺陷数，与扶手式SCNT管相比，锯齿形结构SCNT的平均强度和极限应变相对较小。目前，部分研究利用MD 模拟来探索CNTs 缺陷对其弹性性能的影响，但含缺陷CNTs/高分子纳米复合材料的MD模拟工作却鲜有报道，其内在原子层面的微观机理仍有待探索与揭示。为此，本文建立了含不同缺陷种类和缺陷数量的CNTs增强高分子复合基体分子模型，利用MD理论模拟了高分子复合材料的弹性模量，通过提取相关的微观参数，重点对比与分析CNTs 的表面缺陷种类与数量对高分子复合材料弹性模量的内在微观影响机理，为工程实际提供微观理论支撑。

1 分子建模与动力学模拟

为探索表面含缺陷的CNTs对高分子纳米复合基体弹性模量的微观影响机理，建立含不同缺陷种类和数量的CNTs增强高分子复合基体分子无定型模型。

第1 步，构建不同种类CNTs 分子结构，直径为10.85 nm，长度为59.03 nm 的非缺陷(8，8)CNTs 为碳管的基础构型。图1 所示为表面含缺陷的CNTs 模型，选定与构建thrower-stone-wales(TSW)和single vacancy(SV)这2 种实验中常见的缺陷类型下的CNTs 分子结构[11-12]，分别如图1(a)和(b)所示。其中，TSW 缺陷的构建原理是将C—C键旋转90°后重新组合共价键，原有的4 个相邻的六角形环转换成2个相邻的五角形环和2个分离的七角型环。此时，原有CNTs上的sp2结构并未受到任何影响。SV 缺陷的成型原理则是将表面1 个碳原子删除，由于碳原子的高活性，用氢原子填充余下的不饱和共价键。为探索CNTs表面不同数量缺陷对高分子基体弹性模量的微观影响机理，针对SV和TSW这2种缺陷，分别随机构建了排布有7，14，21和28个缺陷的CNTs，如图1(c)所示。

第2 步，选取分子式为CH2=C[CH3]CO2CH3的聚甲基丙烯酸甲酯(PMMA)为高分子材料。采用Materials Studio(MS)分子模拟软件的Visualizer 建模模块建立聚合度为5的分子链模型，如图2(a)所示。为构建CNTs/PMMA纳米复合基体，建立带有周期性边界条件且长×宽×高为6.5 nm×5.0 nm×5.0 nm 的单元模拟盒子，将CNTs 放置盒子中心位置(CNTs长度方向与盒子长度方向一致)，如图2(b)所示。然后利用Amorphous Cell (AC)模块，基于蒙特卡罗“随机数”的计算方法，将建立的分子链随机填充到含有CNTs的周期性单元盒子，直至填充密度达到1.10 g/cm3[13]。依照此方法，分别拓展建立全部种类CNTs/PMMA的无定型高分子复合材料模型，如图2(c)所示。

图1 表面含缺陷的CNTs模型Fig.1 Molecular models of CNTs with different surface defects

图2 CNTs/PMMA 复合材料模型Fig.2 Molecular modeling process of CNTs/PMMA composites

第3步，为得到局部最优分子构型，对构建的CNTs/PMMA 模型进行几何结构优化处理。其中，优化方法选取共轭梯度法[14]，其能量收敛精度设置为0.418 6 J/mol。为寻找高分子复合基体的全局能量最低构象，对几何优化后的所有模型继续进行NPT 系统下的动力学平衡处理，温度设置为298 K，压力设置为1 个标准大气压，总模拟时间为2×10-9s，模拟步长为1×10-15s。在所有建模、优化与动力学模拟过程中，选取COMPASS[15-18]力场。其中，静电作用能的求解方法为Ewald[19]，截断半径设置为1.25 nm。选取Andersen[20]算法控制温度，选取Berendsen算法[21]控制压力。

2 不同缺陷种类CNTs 增强PMMA基体拉伸计算

为探索CNTs表面不同缺陷种类对高分子纳米复合材料弹性模量的影响，首先，构建表面分别含有14 个TSW 和SV 缺陷的CNTs/PMMA 和无缺陷CNTs/PMMA 的分子无定型复合基体模型；随后，采用恒定应变法(strain-constant method)[22]，对3个分子模型进行MD拉伸模拟。拉伸模拟方法如下：以1.0×10-4为拉伸应变率，在盒子的长度方向均匀扩展晶格尺寸，如图3所示。

图3 CNTs/PMMA基体MD拉伸方法Fig.3 Method of tensile process of CNTs/PMMA composites by MD simulations

其中，每拉伸1步后，分子体系需在NPT系综下(298 K 温度；y和z轴方向施加0.1 MPa 压力)进行1 ps的动力学平衡，使盒子内所有原子重新适应拉伸后尺寸，再根据式(1)计算模型在每步拉伸中的最后一帧的拉伸方向应力，并记录。所有MD拉伸模拟均通过LAMMPS软件实现[23]。

维里应力σ在模型的每个增量方向的计算公式如下：

式中：mα为原子α的质量；分别为原子α在i和j方向的速度分量；为原子β施加在原子α上的i方向的力；为α和β原子在j方向的距离；Vcell为模型体积。

最后，3 个分子模型在线性阶段的应力-应变曲线如图4 所示。经计算，无缺陷CNTs/PMMA，TSW-CNTs/PMMA 和SV-CNTs/PMMA 的弹性模量分别为3.131，2.913 和2.696 GPa。由此可知，当CNTs 表面存在TSW 和SV 缺陷时，其复合基体的弹性模量均呈现降低趋势，分别降低了6.96%和13.89%。其中，相对于TSW 缺陷，含SV 缺陷的CNTs/PMMA复合基体的弹性模量更小。

图4 无缺陷CNTs/PMMA，TSW-CNTs/PMMA，SVCNTs/PMMA的应力-应变曲线Fig.4 Strain-stress curves of non-defects CNTs/PMMA,TSW-CNTs/PMMA and SV-CNTs/PMMA composites

3 CNTs不同缺陷种类的影响机理

根据ARASH 等[24]的研究可知：纳米复合材料基体的弹性模量由补强剂的弹性模量和补强剂与高分子基体之间的界面作用共同决定。因此，为探究CNTs的表面缺陷影响复合基体拉伸性能的微观机理，着重对这2个因素进行探索与分析。

3.1 不同缺陷种类CNTs 与PMMA 基体界面作用能计算与分析

首先根据式(2)，分别提取3种复合材料模型在拉伸过程中CNTs与PMMA基体之间的界面能量变化曲线，如图5 所示。由图5 可见：在拉伸过程中，表面含有无缺陷、TSW 与SV 缺陷的CNTs 与PMMA 基体的界面作用能量的平均值分别为34.95，45.42和40.48 kJ/mol。

式中：UInter为CNTs 与PMMA 基体之间的界面能量；UTotal为CNTs/PMMA 复合基体的总势能；UPMMA为PMMA基体势能；Unano为CNTs的势能。

图5 拉伸过程中不同缺陷种类CNTs与高分子基体界面作用能量变化及CNTs自身结构示意图Fig.5 Variations of interfacial energy between CNTs and polymer matrix during tensile process and molecular models of different CNTs

由图5 可知：当CNTs 表面存在TSW 和SV 缺陷时，在拉伸过程中与PMMA 基体的界面相互作用力得到了一定程度提升。这是因为当SV缺陷存在时，由于苯环缺少1 个碳原子，其CNTs 表面的SP2杂化数量减少，原则上表面吸附力应降低，但由于CNTs未破坏的表面仍有较好吸附活性，同时缺陷处引入的3 个氢原子降低了CNTs 表面的平滑性，CNTs 结构受到一定程度的形变，进而出现了波浪构型与褶皱现象(如图5(c)所示)，因此，在CNTs 结构变化和氢原子的综合作用下，SV-CNTs表面对PMMA分子链的吸附力比无缺陷CNTs模型的吸附力提高了15.8%。

而当CNTs表面存在TSW缺陷时，其缺陷并没有破坏CNTs 原本的SP2的数量与杂化结构，CNTs表面也因结构的改变出现明显的波浪与褶皱现象(如图5(b)所示)，进而改变了CNTs 原有的线性光滑表面特性，同时TSW-CNTs 仍保持原有CNTs 表面的较强静电吸附作用力，导致PMMA 分子链更易于缠绕和吸附在CNTs表面。当复合材料基体被拉伸时，这些表面波浪构型和褶皱现象可阻碍高分子链向拉伸方向移动，从而相比无缺陷CNTs和SV-CNTs，可呈现更高的界面作用能量，可分别提升29.94%和12.20%。

值得注意的是，虽然含TSW 与SV 缺陷的CNTs 与PMMA 基体的界面作用能量更高，但TSW和SV-CNTs/PMMA复合基体的弹性模量相对于无缺陷CNTs/PMMA基体的弹性模量却分别降低了7.03%和13.74%。由此可推测TSW 和SV-CNT与PMMA 之间的界面作用对整体复合材料的拉伸性能影响较小，TSW 和SV-CNTs 自身的弹性模量对复合材料的弹性模量起到主要作用。

3.2 不同缺陷种类CNTs拉伸计算与分析

为证实以上观点，再对无缺陷CNTs，TSWCNTs 与SV-CNTs 这3 种碳管单独模拟应变率为1.0×10-4的分子动力学拉伸，计算并提取其应力应变曲线，如图6 所示。根据图6 可计算出3 种碳纳米管的弹性模量分别为616.6，543.1和469.4 GPa。由此可知：无缺陷CNTs的弹性模量最大，含TSW和含SV缺陷的CNT分别降低了11.92%和23.87%。对于TSW 缺陷，其表面结构相对SV 缺陷比较完整，但其微观构型仍比无缺陷CNTs 的稳定性低；在SV 缺陷结构中，因缺少了1 个碳原子，原本CNTs 的连续苯环结构被打破，缺陷处易诱发应力集中现象，因此，其应力降低，弹性模量最低。

图6 无缺陷CNTs，TSW-CNTs与SV-CNTs的应力-应变曲线Fig.6 Strain-stress curves of non-defects CNTs,TSWCNTs and SV-CNTs

通过对CNTs与PMMA之间的界面能量与单独CNTs弹性模量的综合分析和比对可知，对于TSW和SV-CNTs/PMMA 复合基体，虽然含缺陷CNTs与PMMA基体的界面作用能要比无缺陷CNTs模型的界面作用略高，但TSW 和SV-CNTs 的弹性模量相对于无缺陷CNT 的弹性模量分别降低了11.92%和23.87%。显而易见，此时TSW和SV-CNTs自身的弹性模量在影响复合基体的弹性模量中占据主导作用，导致高分子基体的弹性模量最终呈降低趋势。同时，由于TSW-CNTs 的弹性模量要高于SV-CNTs 的弹性模量15.7%，因此，TSW-CNTs/PMMA 复合基体的弹性模量高于SV-CNTs/PMMA复合基体的弹性模量。

4 不同缺陷数量CNTs 增强PMMA基体拉伸计算

为进一步探索CNTs表面缺陷数量对高分子复合基体拉伸性能的影响，利用MD模拟方法，分别计算表面含有7，14，21和28个TSW和SV缺陷的CNTs/PMMA复合基体分子模型的拉伸行为，其应力-应变曲线如图7 所示。由图7 可见：表面含有7，14，21 和28 个TSW 复合材料基体的弹性模量分别为3.017，2.913，2.622 和2.533 GPa，表面含有7，14，21 和28 个SV-CNTs/PMMA 复合材料基体的弹性模量分别为2.781，2.696，2.412 和2.387 GPa。含有2 种缺陷类型的CNTs/高分子基体均呈现出随着缺陷数量增加而弹性模量减少的趋势。

5 CNTs不同缺陷数量的影响机理

5.1 不同缺陷数量CNTs 与PMMA 基体界面作用能计算与分析

首先，根据式(2)，计算并提取了MD 拉伸模拟过程中，所有模型内TSW和SV-CNTs与PMMA之间的界面作用能量，如图8所示。

图7 不同缺陷种类与缺陷数量下复合材料基体的应力-应变曲线Fig.7 Strain-stress curves of PMMA matrix reinforced by different defected CNTs and different defects numbers

图8 含不同缺陷数量的CNTs与PMMA界面作用能Fig.8 Variations of interfacial energy between PMMA matrix and CNTs with different defects number

由图8 可见：随着缺陷数量增加，TSW-CNT与PMMA 界面间的作用能平均值在缺陷数量为7，14，21 和28 时分别为41.990，45.422，46.758 和47.783 kJ/mol，SV-CNT与PMMA界面间作用能平均值在缺陷数量为7，14，21 和28 时分别为38.670，40.483，41.542 和44.770 kJ/mol。根据此前分析的结论，随着缺陷数量增加，CNTs 本身的构型越发不规则，表面的褶皱和波浪特征逐渐增加，导致在拉伸过程中，CNTs 表面吸附周围PMMA 分子链，提升了界面间的作用力。在相同缺陷数量的情况下，TSW-CNTs 比SV-CNTs 与PMMA 基体之间的界面能相对要高。这是由于引入SV 缺陷，虽然可引入更多的氢原子，但CNTs表面的sp2杂化结构受到损坏，杂化结构的数量随着缺陷增加逐渐减少，而引入TSW 缺陷并没有改变sp2杂化结构和数量，保持着CNTs本身较强的吸附力，同时，其表面的褶皱和波浪特性也随着缺陷数量增加而增强，因此，TSW-CNTs 与PMMA基体之间的界面作用逐渐增加。

5.2 不同缺陷数量CNTs拉伸计算与分析

含有不同缺陷数量下单独TSW 和SV-CNTs 的应力-应变曲线如图9 所示。从图9 可见：随着缺陷数量由7增加到28，TSW-CNTs的弹性模量逐渐由578.73 GPa减少到461.60 GPa；SV-CNTs的弹性模量逐渐由524.56 GPa 减少到320.21 GPa。对于2种缺陷下的CNTs，其弹性模量均随缺陷数量增加而减少，这是因为随着缺陷数量增加，其结构越发不稳定，应力集中现象增多，进而碳管本身承受拉伸能力逐渐降低。

图9 不同缺陷数量下CNTs的应力-应变曲线Fig.9 Strain-stress curves of TSW-CNTs and SV-CNTs with different defects number

由以上计算结果可知：CNTs 自身的弹性模量在影响复合材料基体的弹性模量方面占据主动作用，这与此前的结论相一致。同时，由于随缺陷数量增加，CNTs的弹性模量逐渐降低，导致CNTs增强复合材料基体的弹性模量也随着CNTs表面缺陷数量降低而降低。

在含有相同缺陷数量的情况下，SV-CNTs 的弹性模量要比TSW-CNTs 的小。该现象可解释为：对于SV 缺陷，由于碳管表面碳六边形结构被破坏，结构不连续特性相对更明显，SV-CNTs 结构承载外力能力比TSW-CNTs 更低，因此，由于CNTs 的主导作用，相对于SV-CNTs，在相同缺陷数量下TSW-CNTs 增强高分子复合基体的弹性模量更大。

由文献[13]可知，当纳米补强材料有具有较高的界面作用力时，可有效提升高分子基体的韧性和抗断裂性能。因此，在实际应用中，当在工程需求中强调高分子复合材料的韧性和抗断裂时，可通过技术手段适当增加CNTs 表面的缺陷数量，尤其是TSW 缺陷类型的数量。但当要求确保高分子复合材料的强度时，可尽量减少CNTs表面上缺陷的发生。