搜索

x

留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

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

黄多辉 万明杰 杨俊升

引用本文:
Citation:

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

黄多辉, 万明杰, 杨俊升

Mmolecular dynamics study of glass transition and nonlinear mechanical behavior of poly(methyl methacrylate)/carbon nanotubes nanocomposites

Huang Duo-Hui, Wan Ming-Jie, Yang Jun-Sheng
PDF
HTML
导出引用
  • 短纤维结构对聚合物材料玻璃化转变温度及非线性力学具有非常重要的影响. 本文利用粗粒化分子动力学方法研究了碳纳米管(CNTs)含量对聚甲基丙烯酸甲酯(PMMA)玻璃化转变、扩散系数及非线性力学特性的影响. 分子动力学模拟结果显示: 短CNTs纤维的加入的确会改变PMMA体系的玻璃化转变温度, 模拟结果与实验测量结果一致, 而且随着CNTs含量的增加其对应的玻璃化转变温度也会随着增加. 进一步分析扩散特性发现, CNTs加入PMMA对于体系扩散特性的改变主要发生在玻璃化温度以上, 玻璃化温度以下结构对应的扩散系数差异非常的小. 聚合物材料在服役过程中难免要遭受应力-应变的作用, 而且其结构对应的模量和韧性成反比. 基于此, 本文通过非平衡分子动力学探究了短CNTs纤维添加PMMA复合材料的非线性力学特性. 模拟结果显示: 随着CNTs纤维的含量增加, 其对应的屈服模量也在不断的增加, 而且含有短CNTs纤维的体系还能够保持原来的韧性. 因此, 分子层面的理论研究策略可以为进一步的实验和加工提供理论指导.
    The glass transition temperature and nonlinear mechanics of polymer nanocomposites are strongly influenced by the short fibers. In this paper, coarse-grained molecular dynamics simulations are used to study the effects of single-walled carbon nanotube (CNT) content on the glass transition, diffusion coefficient, viscosity and nonlinear mechanical properties of poly(methyl methacrylate) (PMMA)/CNT nanocomposites. The glass transition temperature Tg is very important for the application of the materials. The Tg is related to the specific volume of the system. Generally, the location of the discontinuity on the curve of specific volume vs. temperature is the position of Tg. Our simulation results show that the Tg of PMMA/CNT composite increases with CNT content, and the result is consistent with the experimental value (434 K). This increase of Tg is evidently due to the presence of CNTs, which imposes a limit on the mobility of the molecules of PMMA. For the free volume in the liquid state, recent experiments pointed out that the molecular mutation is relatively easy to occur because the unoccupied volume is large. Further analysis of the diffusion coefficient of the PMMA/CNT indicates that the difference in diffusion characteristic occurs above the glass transition temperature, and the diffusion coefficient of PMMA system and PMMA/CNT system are the same below the glass transition temperature. Polymer materials in the service process will inevitably suffer the deformation, and the modulus and toughness of material are inversely proportional. Based on this problem, the nonlinear mechanical properties of short CNTs added PMMA composite are studied by nonequilibrium molecular dynamics. Our results show that the yield modulus increases with the CNT content increasing. However, the toughness is almost unchanged. In order to further understand the origin of stress of PMMA/CNT nanocomposites, the stretch ratio and orientation parameters of MPPA chains are also investigated in the present work. According to the stretch ratio and orientation parameters, it is not difficult to conclude that the stress-strain curve is mainly the result of the synergistic effect of molecular chain stretching and orientation. This work provides a theoretical guidance for further experiments and processing at the atomic and molecular level.
      通信作者: 杨俊升, yangjunsheng2005@163.com
    • 基金项目: 宜宾学院预研项目 (批准号: 2019YY06) 和宜宾学院计算物理四川省高等学校重点实验室开放课题基金(批准号: YBXYJSWL-ZD-2020-003)资助的课题.
      Corresponding author: Yang Jun-Sheng, yangjunsheng2005@163.com
    • Funds: Project supported by the Pre-Research Project of Yibin University, China (Grant No. 2019YY06) and the Open Research Fund of Computational Physics Key Laboratory of Sichuan Province, Yibin University, China (Grant No. YBXYJSWL-ZD-2020-003)
    [1]

    Mackay M E, Dao T T, Tuteja A, Ho D L, van Horn B, Kim H C, Hawker C J 2003 Nat. Mater. 2 762Google Scholar

    [2]

    Polizos G, Tuncer E, Agapov A L, Stevens D, Sokolov A P, Kidder M K, Jacobs J D, Koerner H, Vaia R A, More K L, Sauers I 2012 Polymer 53 595Google Scholar

    [3]

    Mammeri F, Teyssandier J, Connan C, Le Bourhis E, Chehimi M M 2012 RSC Adv. 2 2462Google Scholar

    [4]

    Luo J T, Wen H C, Wu W F, Chou C P 2008 Polym. Compos. 29 1285Google Scholar

    [5]

    De S K 1996 Short Fibre-polymer Composites (Woodhead: Woodhead Publishing) p257

    [6]

    Salami-Kalajahi M, Haddadi-Asl V, Behboodi-Sadabad F, Rahimi-Razin S, Roghani-Mamaqani H 2012 Polym. Compos. 33 215Google Scholar

    [7]

    Van Loock F, Fleck N A 2018 Polymer 148 259Google Scholar

    [8]

    Haggenmueller R, Gommans H H, Rinzler A G, Fischer J E, Winey K I 2000 Chem. Phys. Lett. 330 219Google Scholar

    [9]

    Du F, Fischer J E, Winey K I 2003 J. Polym. Sci. , Part B:Polym. Phys. 41 3333Google Scholar

    [10]

    Wang J F, Yang J P, Tam L h, Zhang W 2021 Mech. Syst. Singal PR 153 107530Google Scholar

    [11]

    杨俊升, 朱子亮, 曹启龙 2020 物理学报 69 038101

    Yang J S, Zhu Z L, Cao Q L 2020 Acta Phys. Sin. 69 038101

    [12]

    陈超, 段芳莉 2020 物理学报 69 1 93102Google Scholar

    Chen C, Duan F L 2020 Acta Phys. Sin. 69 1 93102Google Scholar

    [13]

    Yang J S, Yang C L, Wang M S, Chen B D, Ma X G 2011 Phys. Chem. Chem. Phys. 13 15476Google Scholar

    [14]

    Skountzos E N, Mermigkis P G, Mavrantzas V G 2018 J. Phys. Chem. B 122 9007Google Scholar

    [15]

    Jiang Q, Tallury S S, Qiu Y, Pasquinelli M A 2020 Nanotechnol. Rev. 9 136Google Scholar

    [16]

    Mohammadi M, fazli H, karevan M, Davoodi J 2017 Eur. Polym. J. 91 121Google Scholar

    [17]

    Zhao J, Jiang J W, Wang L, Guo W, Rabczuk T 2014 J. Mech. Phys. Solids 71 197Google Scholar

    [18]

    Arash B, Park H S, Rabczuk T 2015 Compos. Struct. 134 981Google Scholar

    [19]

    Meng Z, Soler-Crespo R A, Xia W, Gao W, Ruiz L, Espinosa H D, Keten S 2017 Carbon 117 476Google Scholar

    [20]

    Mousavi A A, Arash B, Zhuang X, Rabczuk T 2016 Compos. B. Eng. 95 404Google Scholar

    [21]

    Ash B J, Siegel R W, Schadler L S 2004 J. Polym. Sci. , Part B:Polym. Phys. 42 4371Google Scholar

    [22]

    Mohammadi M, Davoodi J, Javanbakht M, Rezaei H 2018 Mater. Res. Express 6 035309Google Scholar

  • 图 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.

    图 3  PMMA, PMMA/10CNTs及PMMA/20CNTs体系对应的比容随着温度的变化

    Fig. 3.  Specific volumes of PMMA, PMMA/10CNTs and PMMA/20CNTs under the different temperatures.

    图 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.

    图 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.

    图 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).

    图 7  PMMA体系维诺体积随着应变的变化过程

    Fig. 7.  Voronoi volume of PMMA system under the different stain.

    图 8  PMMA/10CNTs体系维诺体积随着应变的变化过程

    Fig. 8.  Voronoi volume of PMMA/10CNTs system under the different stain.

    表 1  PMMA及PMMA-CNT粗粒化粒子之间相互作用力场参数

    Table 1.  Parameters of the CG force field for PMMA and CNT beads.

    Type of interactionParametersPMMACNTPMMA
    /CNT
    Bondkd/(kcal·mol–1·Å–2)97.3805.15
    d04.059.45
    Anglekθ/(kcal·mol–1·Å–2)40032140
    θ0/(o)84.6180
    vdWD0/(kcal·mol–1)0.565.341.4
    r06.539.457.7
    下载: 导出CSV
  • [1]

    Mackay M E, Dao T T, Tuteja A, Ho D L, van Horn B, Kim H C, Hawker C J 2003 Nat. Mater. 2 762Google Scholar

    [2]

    Polizos G, Tuncer E, Agapov A L, Stevens D, Sokolov A P, Kidder M K, Jacobs J D, Koerner H, Vaia R A, More K L, Sauers I 2012 Polymer 53 595Google Scholar

    [3]

    Mammeri F, Teyssandier J, Connan C, Le Bourhis E, Chehimi M M 2012 RSC Adv. 2 2462Google Scholar

    [4]

    Luo J T, Wen H C, Wu W F, Chou C P 2008 Polym. Compos. 29 1285Google Scholar

    [5]

    De S K 1996 Short Fibre-polymer Composites (Woodhead: Woodhead Publishing) p257

    [6]

    Salami-Kalajahi M, Haddadi-Asl V, Behboodi-Sadabad F, Rahimi-Razin S, Roghani-Mamaqani H 2012 Polym. Compos. 33 215Google Scholar

    [7]

    Van Loock F, Fleck N A 2018 Polymer 148 259Google Scholar

    [8]

    Haggenmueller R, Gommans H H, Rinzler A G, Fischer J E, Winey K I 2000 Chem. Phys. Lett. 330 219Google Scholar

    [9]

    Du F, Fischer J E, Winey K I 2003 J. Polym. Sci. , Part B:Polym. Phys. 41 3333Google Scholar

    [10]

    Wang J F, Yang J P, Tam L h, Zhang W 2021 Mech. Syst. Singal PR 153 107530Google Scholar

    [11]

    杨俊升, 朱子亮, 曹启龙 2020 物理学报 69 038101

    Yang J S, Zhu Z L, Cao Q L 2020 Acta Phys. Sin. 69 038101

    [12]

    陈超, 段芳莉 2020 物理学报 69 1 93102Google Scholar

    Chen C, Duan F L 2020 Acta Phys. Sin. 69 1 93102Google Scholar

    [13]

    Yang J S, Yang C L, Wang M S, Chen B D, Ma X G 2011 Phys. Chem. Chem. Phys. 13 15476Google Scholar

    [14]

    Skountzos E N, Mermigkis P G, Mavrantzas V G 2018 J. Phys. Chem. B 122 9007Google Scholar

    [15]

    Jiang Q, Tallury S S, Qiu Y, Pasquinelli M A 2020 Nanotechnol. Rev. 9 136Google Scholar

    [16]

    Mohammadi M, fazli H, karevan M, Davoodi J 2017 Eur. Polym. J. 91 121Google Scholar

    [17]

    Zhao J, Jiang J W, Wang L, Guo W, Rabczuk T 2014 J. Mech. Phys. Solids 71 197Google Scholar

    [18]

    Arash B, Park H S, Rabczuk T 2015 Compos. Struct. 134 981Google Scholar

    [19]

    Meng Z, Soler-Crespo R A, Xia W, Gao W, Ruiz L, Espinosa H D, Keten S 2017 Carbon 117 476Google Scholar

    [20]

    Mousavi A A, Arash B, Zhuang X, Rabczuk T 2016 Compos. B. Eng. 95 404Google Scholar

    [21]

    Ash B J, Siegel R W, Schadler L S 2004 J. Polym. Sci. , Part B:Polym. Phys. 42 4371Google Scholar

    [22]

    Mohammadi M, Davoodi J, Javanbakht M, Rezaei H 2018 Mater. Res. Express 6 035309Google Scholar

  • [1] 邓永和, 张宇文, 谭恒博, 文大东, 高明, 吴安如. NiCu双金属纳米粒子的表面偏析、结构特征与扩散. 物理学报, 2021, 70(17): 177601. doi: 10.7498/aps.70.20210336
    [2] 廖晶晶, 蔺福军. 混合手征活性粒子在时间延迟反馈下的扩散和分离. 物理学报, 2020, 69(22): 220501. doi: 10.7498/aps.69.20200505
    [3] 张恒, 黄燕, 石旺舟, 周孝好, 陈效双. Al原子在Si表面扩散动力学的第一性原理研究. 物理学报, 2019, 68(20): 207302. doi: 10.7498/aps.68.20190783
    [4] 李亚雄, 刘先贵, 胡志明, 高树生, 端祥刚, 常进. 页岩气滑脱、扩散传输机理耦合新方法. 物理学报, 2017, 66(11): 114702. doi: 10.7498/aps.66.114702
    [5] 杨文龙, 韩浚生, 王宇, 林家齐, 何国强, 孙洪国. 聚酰亚胺/功能化石墨烯复合材料力学性能及玻璃化转变温度的分子动力学模拟. 物理学报, 2017, 66(22): 227101. doi: 10.7498/aps.66.227101
    [6] 林家齐, 李晓康, 杨文龙, 孙洪国, 谢志滨, 修翰江, 雷清泉. 聚酰亚胺/钽铌酸钾纳米颗粒复合材料结构与机械性能分子动力学模拟. 物理学报, 2015, 64(12): 126202. doi: 10.7498/aps.64.126202
    [7] 刘式达, 付遵涛, 刘式适. 间歇湍流的分数阶动力学. 物理学报, 2014, 63(7): 074701. doi: 10.7498/aps.63.074701
    [8] 向辉, 刘大欢, 阳庆元, 密建国, 仲崇立. 骨架柔性对短链烷烃分子在金属-有机骨架材料中扩散的影响. 物理学报, 2011, 60(9): 093602. doi: 10.7498/aps.60.093602
    [9] 陈德彝, 王忠龙. 白交叉关联色噪声驱动的线性振子的扩散. 物理学报, 2010, 59(1): 111-115. doi: 10.7498/aps.59.111
    [10] 杨春, 冯玉芳, 余毅. 动力学研究AlN/α-Al2O3(0001)薄膜生长初期的吸附与扩散. 物理学报, 2009, 58(5): 3553-3559. doi: 10.7498/aps.58.3553
    [11] 李美丽, 张 迪, 孙宏宁, 付兴烨, 姚秀伟, 李 丛, 段永平, 闫 元, 牟洪臣, 孙民华. 二元Lennard-Jones液体的相分离过程及其扩散性质的分子动力学研究. 物理学报, 2008, 57(11): 7157-7163. doi: 10.7498/aps.57.7157
    [12] 刘 磊, 任晓敏, 周 静, 王 琦, 熊德平, 黄 辉, 黄永清. 横向外延过生长磷化铟材料的生长速率模型. 物理学报, 2007, 56(6): 3570-3576. doi: 10.7498/aps.56.3570
    [13] 曹 博, 包良满, 李公平, 何山虎. Cu/SiO2/Si(111)体系中Cu和Si的扩散及界面反应. 物理学报, 2006, 55(12): 6550-6555. doi: 10.7498/aps.55.6550
    [14] 乔永红, 王绍青. 硅晶体中点缺陷结合过程的分子动力学研究. 物理学报, 2005, 54(10): 4827-4835. doi: 10.7498/aps.54.4827
    [15] 常福宣, 陈 进, 黄 薇. 反常扩散与分数阶对流-扩散方程. 物理学报, 2005, 54(3): 1113-1117. doi: 10.7498/aps.54.1113
    [16] 唐 鑫, 张 超, 张庆瑜. Cu(111)三维表面岛对表面原子扩散影响的分子动力学研究. 物理学报, 2005, 54(12): 5797-5803. doi: 10.7498/aps.54.5797
    [17] 杨 春, 余 毅, 李言荣, 刘永华. 温度对ZnO/Al2O3(0001)界面的吸附、扩散及生长初期模式的影响. 物理学报, 2005, 54(12): 5907-5913. doi: 10.7498/aps.54.5907
    [18] 王秀英, 孙力玲, 刘日平, 姚玉书, 张 君, 王文魁. 高压下Co在Zr46.75Ti8.25Cu7.5Ni10Be27.5大块金属玻璃过冷液相区中的扩散. 物理学报, 2004, 53(11): 3845-3848. doi: 10.7498/aps.53.3845
    [19] 胡晓君, 戴永兵, 何贤昶, 沈荷生, 李荣斌. 空位在金刚石近(001)表面扩散的分子动力学模拟. 物理学报, 2002, 51(6): 1388-1392. doi: 10.7498/aps.51.1388
    [20] 金进生, 叶高翔, 钱昌吉, 翟国庆, 叶全林, 焦正宽. 金原子在熔融玻璃表面的凝聚特性. 物理学报, 2001, 50(3): 544-549. doi: 10.7498/aps.50.544
计量
  • 文章访问数:  1082
  • PDF下载量:  40
  • 被引次数: 0
出版历程
  • 收稿日期:  2021-04-20
  • 修回日期:  2021-06-22
  • 上网日期:  2021-08-15
  • 刊出日期:  2021-11-05

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

  • 1. 宜宾学院理学部, 宜宾 644007
  • 2. 华南理工大学, 华南软物质科学与技术高等研究院, 分子科学与工程学院, 广州 510640
  • 通信作者: 杨俊升, yangjunsheng2005@163.com
    基金项目: 宜宾学院预研项目 (批准号: 2019YY06) 和宜宾学院计算物理四川省高等学校重点实验室开放课题基金(批准号: YBXYJSWL-ZD-2020-003)资助的课题.

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

English Abstract

    • 通过短纳米纤维增强聚合物材料的力学、热学及电学性能的研究最近几年引起了研究者广泛的关注[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玻璃化转变、扩散系数及相关拉伸和取向结构与非线性力学的关系.

    • 本文使用了前人最近开发的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模型

      Figure 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结构示意图

      Figure 2.  Schematic diagram of PMMA and PMMA/CNTs.

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

      $ {E}_{\mathrm{v}\mathrm{d}\mathrm{W}}\left(r\right)={D}_{0}\left[{\left(\frac{{r}_{0}}{r}\right)}^{12}-2{\left(\frac{{r}_{0}}{r}\right)}^{6}\right], $

      $ {E}_{\rm b}\left(d\right)=\frac{{k}_{d}}{2}{\left(d-{d}_{0}\right)}^{2}, $

      $ {E}_{\rm a}\left(\theta \right)=\frac{{k}_{\theta }}{2}{\left(\theta -{\theta }_{0}\right)}^{2}, $

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

      Type of interactionParametersPMMACNTPMMA
      /CNT
      Bondkd/(kcal·mol–1·Å–2)97.3805.15
      d04.059.45
      Anglekθ/(kcal·mol–1·Å–2)40032140
      θ0/(o)84.6180
      vdWD0/(kcal·mol–1)0.565.341.4
      r06.539.457.7

      表 1  PMMA及PMMA-CNT粗粒化粒子之间相互作用力场参数

      Table 1.  Parameters of the CG force field for PMMA and CNT beads.

    • 玻璃化转变温度对于聚合物纳米复合材料体系的应用起着至关重要的作用. 图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体系对应的比容随着温度的变化

      Figure 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]:

      $ D=\underset{t\to \infty }{\mathrm{lim}}\frac{1}{6t}\langle {\left[r\left(t+{t}_{0}\right)-r\left({t}_{0}\right)\right]}^{2} \rangle , $

      其中, 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体系不同温度下对应的扩散系数

      Figure 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.

    • 聚合物材料在服役过程中难免要遭受外力场的作用, 深刻认识外力作用下材料的结构和性能关系不仅具有重要的学术价值, 也是决定材料能否长期有效服役的关键因素. 聚合物纳米复合材料凝聚态结构是由大量的高分子依靠分子内和分子间的范德瓦耳斯相互作用凝聚而成, 而且材料最后的拉伸性能是这些微观结构相互作用的结果. 因此, 非常有必要探究聚合物纳米复合材料结构与其非线性力学关系. 图5(a)中, 我们施加了一个恒定的拉伸速率1×109 s–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)变化曲线

      Figure 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分子链对应的链微观拉伸比和末端矢量取向参数. 采用$\lambda =({L-{L}_{0}})/{{L}_{0}}$计算链结构在宏观拉伸过程对应的微观拉伸比, 其中, L0L分别表示链原始末端距和受到拉伸以后的末端距, 采用公式$P=\dfrac{3}{2}\cos$θ-$\dfrac{1}{2}$计算取向参数, 其中, θ代表末端矢量与拉伸方向的夹角. 不难发现, 最开始阶段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)

      Figure 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次屈服软化点附近.

      图  7  PMMA体系维诺体积随着应变的变化过程

      Figure 7.  Voronoi volume of PMMA system under the different stain.

      图  8  PMMA/10CNTs体系维诺体积随着应变的变化过程

      Figure 8.  Voronoi volume of PMMA/10CNTs system under the different stain.

    • 本文通过粗粒化的MD方法对PMMA及PMMA/CNTs聚合物复合体系的玻璃化转变、扩散、黏性及非线性力学特性进行了系统研究, 具体结论如下:

      1) 通过比容变化, 给出了不同温度下PMMA及PMMA/CNTs体系的特征体积, 随着CNT含量的增加引起特征体积的进一步减小(密度增加). 根据特征体系拟合出了PMMA及PMMA/CNTs体系的玻璃化转变温度, 玻璃化转变温度会随着CNTs纤维含量的增加而增加, 这些结果与实验一致.

      2) 通过分析PMMA及PMMA/CNTs体系扩散特性发现, CNTs纤维加入PMMA对于体系扩散特性的改变主要发生在玻璃化温度以上. 而且PMMA, PMMA/10CNTs及PMMA/20CNTs体系扩散系数在玻璃化温度以上成指数增长.

      3) 通过CNTs纤维添加PMMA复合材料的非线性力学特性分析发现, 随着CNTs纤维的含量增加, 其对应的屈服模量也在不断的增加, 但是含有CNTs的体系还能够保持原来的韧性, 这一点可以为进一步的实验和加工提供理论指导.

参考文献 (22)

目录

    /

    返回文章
    返回