搜索

文章查询

x

留言板

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

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

预取向半晶态高分子片晶结构形成微观机理及其应力-应变响应特性的分子动力学模拟

杨俊升 朱子亮 曹启龙

预取向半晶态高分子片晶结构形成微观机理及其应力-应变响应特性的分子动力学模拟

杨俊升, 朱子亮, 曹启龙
PDF
HTML
导出引用
导出核心图
  • 通过分子动力学模拟方法对不同预取向聚乙烯醇熔体(polyvinyl alcohol, PVA)形成的半晶态高分子熔体形成核结晶及拉伸过程中的应力-应变响应特性进行了系统地研究. 模拟结果显示预取向度高的PVA熔体对应更快的成核动力学. 通过追踪全trans伸直链长度(d tt)、成核原子维诺体积(V )和中心对称参数(S)等序参量在不同取向度熔体下的等温成核与结晶演化过程, 给出了PVA熔体成核路径及形成半晶态的分子构象; 通过对形成的半晶态高分子结构进一步分析, 发现随着熔体取向度的增加, 晶体和无定型对应的取向度也会增加, 但是当应变剪切大于5时, 其对应的结晶度、晶体和无定型取向度不再发生变化; 通过对无定型区链结构的定义与分析, 可知取向度越高的熔体对应越高的Tie链数目, 随着熔体取向度的增加, Loop链的数目也会减小; 通过恒速拉伸应力测试可知, 所形成半晶态高分子结构力学响应会随着取向度及Tie链数增加而增加, 当取向及Tie链数目相同时, 应力-应变曲线形状大小也基本保持一致.
      通信作者: 杨俊升, yangjunsheng2005@163.com ; 曹启龙, qlcao@mail.ustc.edu.cn
    • 基金项目: 国家自然科学基金(批准号: 11704329)、四川省教育厅科研基金(批准号: 15ZB293)和宜宾学院计算物理四川省高等学校重点实验室开放课题基金(批准号: JSWL2014KFZ02)资助的课题
    [1]

    Li L, de Jeu W 2005 Adv. Polym. Sci. 181 75

    [2]

    Liu D, Cui K, Huang N, Wang Z, Li L 2015 Sci. China Chem. 58 1570

    [3]

    Stephanou P S, Tsimouri I C, Mavrantzas V G 2016 Macromolecules 49 3161

    [4]

    Luo C, Kröger M, Sommer J U 2016 Macromolecules 49 9017

    [5]

    Luo C, Kröger M, Sommer J U 2017 Polymer 109 71

    [6]

    Tang X, Yang J, Xu T, Tian F, Xie C, Li L 2017 Phys. Rev. Mate. 1 073401

    [7]

    Yang J, Tang X, Wang Z, Xu T, Tian F, Ji Y, Li L 2017 J. Chem. Phys. 146 014901

    [8]

    Tang X, Yang J, Tian F, Xu T, Xie C, Chen W, Li L 2018 J. Chem. Phys. 149 224901

    [9]

    Yamamoto T 2014 Macromolecules 47 3192

    [10]

    Baig C, Edwards B J 2010 Europhys. Lett. 89 36003

    [11]

    Yang J S, Yang C, Wang M, Chen B, Ma X 2011 Phys. Chem. Chem. Phys. 13 15476

    [12]

    Yang J S, Huang D H, Cao Q, Li Q, Wang L, Wang F 2013 Chin. Phys. B 22 098101

    [13]

    杨俊升, 黄多辉 2019 物理学报 68 138301

    Yang J S, Huang D H 2019 Acta Phys. Sin. 68 138301

    [14]

    杨文龙, 韩浚生, 王宇, 林家齐, 何国强, 孙洪国 2017 物理学报 66 227101

    Yang W L, Han J S, Wang Y, Lin J Q, He G Q, Sun H G 2017 Acta Phys. Sin. 66 227101

    [15]

    潘登, 刘长鑫, 张泽洋, 高玉金, 郝秀红 2019 物理学报 68 176801

    Pan D, Liu C X, Zhang Z Y, Gao Y J, Hao X H 2019 Acta Phys. Sin. 68 176801

    [16]

    Cui K, Ma Z, Wang Z, Ji Y, Liu D, Huang N, Chen L, Zhang W, Li L 2015 Macromolecules 48 5276

    [17]

    Cui K, Meng L, Ji Y, Li J, Zhu S, Li X, Tian N, Liu D, Li L 2014 Macromolecules 47 677

    [18]

    Luo C, Sommer J 2016 ACS Macro Lett. 5 30

    [19]

    Tang X, Chen W, Li L 2019 Macromolecules 52 3575

    [20]

    Meyer H, Müller-Plathe F 2001 J. Chem. Phys. 115 7807

    [21]

    Luo C, Sommer J 2009 Comp. Phys. Comm. 180 1382

    [22]

    Wang S, Wang Y, Cheng S, Li X, Zhu X, Sun H 2013 Macromolecules 46 3147

    [23]

    Kelchner C L, Plimpton S, Hamilton J 1998 Phys. Rev. B 58 11085

    [24]

    Wang Z, Ju J, Yang J, Ma Z, Liu D, Cui K, Yang H, Chang J, Huang N, Li L 2016 Sci. Rep. 6 32968

    [25]

    Wang Y, Jiang Z, Wu Z, Men Y 2013 Macromolecules 45 518

    [26]

    Siviour C R, Jordan J L 2016 J. Dyn. Behav. Mater. 2 15

    [27]

    Lin Y, Li X, Meng L, Chen X, Li L 2018 Macromolecules 51 2690

    [28]

    Lin Y, Li X, Meng L, Chen X, Lü F, Zhang Q, Li L 2018 Polymer 148 79

  • 图 1  (a)高分子熔体剪切示意图; (b)高分子熔体在剪切场温度为T = 1.0下的σxy, RetePall演化曲线

    Fig. 1.  (a) Schematic diagram of polymer melts under shear flow; (b) evolutions of σxy, Rete and Pall for polymer melts under shear with temperature of T = 1.0.

    图 2  (a)成核原子判定方法示意图; (b)伸直链判定方法示意图; (c) PVA半晶态高分子结构晶区与无定型区S值的分布

    Fig. 2.  Schematic diagram of (a) centro-symmetry parameter and (b) stretched chain segment for PVA crystal; (c) S distribution of crystal and melt of PVA semicrystalline polymers.

    图 3  不同应变下PVA熔体的${\phi _{\rm{c}}}$, d tt, VS随时间的演化过程 (a) γ = 0; (b) γ = 1; (c) γ = 1.5; (d) γ = 2; (e) γ = 4; (f) γ = 8

    Fig. 3.  Evolutions of ${\phi _{\rm{c}}}$, d tt, V and S of nucleation atoms under different stains: (a) γ = 0; (b) γ = 1; (c) γ = 1.5; (d) γ = 2; (e) γ = 4; (f) γ = 8.

    图 4  (a)半晶态高分子结构中无定型对应的链结构分类模型; (b), (c) PVA半晶态高分子结构中晶体和无定型的原子结构

    Fig. 4.  (a) Schematic diagram of crystalline and amorphous chain structure for the semicrystalline polymers; (b) and (c) atomic structure of crystal and melt of PVA semicrystalline polymers.

    图 5  (a)不同剪切应下PVA半晶态高分子对应的${\phi _{\rm{c}}}$、晶体取向参数Pc和无定型结构取向参数Pa; (b)不同剪切应变下PVA半晶态高分子无定型链结构数目的演化

    Fig. 5.  (a) ${\phi _{\rm{c}}}$, crystalline order parameter Pc, and amorphous order parameter Pa for PVA semicrystalline polymers with different shear strains; (b) the evolution of the numbers of amorphous chains for PVA semicrystalline polymers with different shear strains.

    图 6  不同预取向半晶态高分子结构在恒定速率为$1 \times {10^{{\rm{ - }}5}}{\tau ^{ - 1}}$拉伸场下对应的应力-应变曲线

    Fig. 6.  Stress-strain curves of PVA semicrystalline polymer with different orientation degree under the stretched rate of $1 \times {10^{{\rm{ - }}5}}{\tau ^{ - 1}}$.

    图 7  γ = 6时对应的单链构象随着应变ε的演化过程

    Fig. 7.  Structural evolution of single PVA chain as a function of strain ε when γ = 6.

  • [1]

    Li L, de Jeu W 2005 Adv. Polym. Sci. 181 75

    [2]

    Liu D, Cui K, Huang N, Wang Z, Li L 2015 Sci. China Chem. 58 1570

    [3]

    Stephanou P S, Tsimouri I C, Mavrantzas V G 2016 Macromolecules 49 3161

    [4]

    Luo C, Kröger M, Sommer J U 2016 Macromolecules 49 9017

    [5]

    Luo C, Kröger M, Sommer J U 2017 Polymer 109 71

    [6]

    Tang X, Yang J, Xu T, Tian F, Xie C, Li L 2017 Phys. Rev. Mate. 1 073401

    [7]

    Yang J, Tang X, Wang Z, Xu T, Tian F, Ji Y, Li L 2017 J. Chem. Phys. 146 014901

    [8]

    Tang X, Yang J, Tian F, Xu T, Xie C, Chen W, Li L 2018 J. Chem. Phys. 149 224901

    [9]

    Yamamoto T 2014 Macromolecules 47 3192

    [10]

    Baig C, Edwards B J 2010 Europhys. Lett. 89 36003

    [11]

    Yang J S, Yang C, Wang M, Chen B, Ma X 2011 Phys. Chem. Chem. Phys. 13 15476

    [12]

    Yang J S, Huang D H, Cao Q, Li Q, Wang L, Wang F 2013 Chin. Phys. B 22 098101

    [13]

    杨俊升, 黄多辉 2019 物理学报 68 138301

    Yang J S, Huang D H 2019 Acta Phys. Sin. 68 138301

    [14]

    杨文龙, 韩浚生, 王宇, 林家齐, 何国强, 孙洪国 2017 物理学报 66 227101

    Yang W L, Han J S, Wang Y, Lin J Q, He G Q, Sun H G 2017 Acta Phys. Sin. 66 227101

    [15]

    潘登, 刘长鑫, 张泽洋, 高玉金, 郝秀红 2019 物理学报 68 176801

    Pan D, Liu C X, Zhang Z Y, Gao Y J, Hao X H 2019 Acta Phys. Sin. 68 176801

    [16]

    Cui K, Ma Z, Wang Z, Ji Y, Liu D, Huang N, Chen L, Zhang W, Li L 2015 Macromolecules 48 5276

    [17]

    Cui K, Meng L, Ji Y, Li J, Zhu S, Li X, Tian N, Liu D, Li L 2014 Macromolecules 47 677

    [18]

    Luo C, Sommer J 2016 ACS Macro Lett. 5 30

    [19]

    Tang X, Chen W, Li L 2019 Macromolecules 52 3575

    [20]

    Meyer H, Müller-Plathe F 2001 J. Chem. Phys. 115 7807

    [21]

    Luo C, Sommer J 2009 Comp. Phys. Comm. 180 1382

    [22]

    Wang S, Wang Y, Cheng S, Li X, Zhu X, Sun H 2013 Macromolecules 46 3147

    [23]

    Kelchner C L, Plimpton S, Hamilton J 1998 Phys. Rev. B 58 11085

    [24]

    Wang Z, Ju J, Yang J, Ma Z, Liu D, Cui K, Yang H, Chang J, Huang N, Li L 2016 Sci. Rep. 6 32968

    [25]

    Wang Y, Jiang Z, Wu Z, Men Y 2013 Macromolecules 45 518

    [26]

    Siviour C R, Jordan J L 2016 J. Dyn. Behav. Mater. 2 15

    [27]

    Lin Y, Li X, Meng L, Chen X, Li L 2018 Macromolecules 51 2690

    [28]

    Lin Y, Li X, Meng L, Chen X, Lü F, Zhang Q, Li L 2018 Polymer 148 79

  • [1] 胡蕴成, 叶祥熙, 明辰, 宁西京. 体材料结晶能力的理论预测. 物理学报, 2009, 58(5): 3293-3301. doi: 10.7498/aps.58.3293
    [2] 徐树杰, 师春生, 赵乃勤, 刘恩佐. 热加工过程中动态再结晶现象的多相场研究. 物理学报, 2012, 61(11): 116101. doi: 10.7498/aps.61.116101
    [3] 杨俊升, 黄多辉. 环状聚合物及其对应的线性链熔体在启动剪切场下流变特性的分子动力学模拟研究. 物理学报, 2019, 68(13): 138301. doi: 10.7498/aps.68.20190403
    [4] 黄 文, 曾慧中, 张 鹰, 蒋书文, 魏贤华, 李言荣. 不同晶化工艺对非晶PZT纳米薄膜形核取向生长机理的影响. 物理学报, 2005, 54(3): 1334-1340. doi: 10.7498/aps.54.1334
    [5] 严大东, 张兴华. 聚合物结晶理论进展. 物理学报, 2016, 65(18): 188201. doi: 10.7498/aps.65.188201
    [6] 陶为俊, 浣石. 沿时间逐步求解应力的拉格朗日分析方法研究. 物理学报, 2012, 61(20): 200703. doi: 10.7498/aps.61.200703
    [7] 贝帮坤, 王华光, 张泽新. 有限尺寸胶体体系的二维结晶. 物理学报, 2019, 68(10): 106401. doi: 10.7498/aps.68.20190304
    [8] 李酽, 张琳彬, 李娇, 连晓雪, 朱俊武. 电场条件下氧化锌结晶特性及极化产物的拉曼光谱分析. 物理学报, 2019, 68(7): 070701. doi: 10.7498/aps.68.20181961
    [9] 王理林, 王志军, 林鑫, 王锦程, 黄卫东. 冷却速率对温敏聚N-异丙基丙烯酰胺胶体结晶过程的影响. 物理学报, 2016, 65(10): 106403. doi: 10.7498/aps.65.106403
    [10] 刘艳松, 陈 铠, 乔 峰, 黄信凡, 韩培高, 钱 波, 马忠元, 李 伟, 徐 骏, 陈坤基. 尺寸可控的纳米硅的生长模型和实验验证. 物理学报, 2006, 55(10): 5403-5408. doi: 10.7498/aps.55.5403
    [11] 王小霞, 罗积润, 赵青兰, 王其富. 颗粒状纳米碳酸钡锶钙的研制. 物理学报, 2010, 59(10): 7383-7389. doi: 10.7498/aps.59.7383
    [12] 王庆学. 异质结构的应变和应力分布模型研究. 物理学报, 2005, 54(8): 3757-3763. doi: 10.7498/aps.54.3757
    [13] 贾琳, 王理林, 申洁楠, 张忠明, 李俊杰, 王锦程, 王志军. 聚乙烯醇水溶液二维定向凝固的微观组织演化. 物理学报, 2017, 66(19): 196402. doi: 10.7498/aps.66.196402
    [14] 冯奇, 李梦凯, 唐海通, 王晓东, 高忠民, 孟繁玲. 石墨烯/聚乙烯醇/聚偏氟乙烯基纳米复合薄膜的介电性能. 物理学报, 2016, 65(18): 188101. doi: 10.7498/aps.65.188101
    [15] 陈梦谪, 柯俊, 李树尘. 冷轧取向硅钢二次再结晶过程的能量分析. 物理学报, 1990, 39(4): 672-676. doi: 10.7498/aps.39.672
    [16] 韩甫田, 刘平安, 唐振方, 施其宏, 郭立平. 半结晶聚酯(PET)的二相共存结构的表征. 物理学报, 2001, 50(6): 1132-1138. doi: 10.7498/aps.50.1132
    [17] 赵辉, 董宝中, 郭梅芳, 王良诗, 乔金梁. 小角x射线散射结晶聚合物结构的研究. 物理学报, 2002, 51(12): 2887-2891. doi: 10.7498/aps.51.2887
    [18] 田丽丽, 王楠, 彭银利, 姚文静. 定向结晶条件下聚乙二醇6000的强动力学效应. 物理学报, 2016, 65(9): 096401. doi: 10.7498/aps.65.096401
    [19] 宫本恭幸, 小田俊理, 姚 尧, 方忠慧, 周 江, 李 伟, 马忠元, 徐 骏, 黄信凡, 陈坤基. 激光干涉结晶法制备一维周期结构的纳米硅阵列. 物理学报, 2008, 57(8): 4960-4965. doi: 10.7498/aps.57.4960
    [20] 刘建军, 张万林, 张光寅. 掺镁铌酸锂晶体的缺陷结构及其结晶化学分析. 物理学报, 1996, 45(11): 1852-1858. doi: 10.7498/aps.45.1852
  • 引用本文:
    Citation:
计量
  • 文章访问数:  356
  • PDF下载量:  22
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-08-04
  • 修回日期:  2019-11-28
  • 刊出日期:  2020-02-01

预取向半晶态高分子片晶结构形成微观机理及其应力-应变响应特性的分子动力学模拟

  • 1. 宜宾学院, 计算物理四川省高等学校重点实验室, 宜宾 644000
  • 2. 潍坊科技学院, 寿光 262700
  • 3. 青岛大学材料科学与工程学院, 青岛 266003
  • 通信作者: 杨俊升, yangjunsheng2005@163.com ; 曹启龙, qlcao@mail.ustc.edu.cn
    基金项目: 国家自然科学基金(批准号: 11704329)、四川省教育厅科研基金(批准号: 15ZB293)和宜宾学院计算物理四川省高等学校重点实验室开放课题基金(批准号: JSWL2014KFZ02)资助的课题

摘要: 通过分子动力学模拟方法对不同预取向聚乙烯醇熔体(polyvinyl alcohol, PVA)形成的半晶态高分子熔体形成核结晶及拉伸过程中的应力-应变响应特性进行了系统地研究. 模拟结果显示预取向度高的PVA熔体对应更快的成核动力学. 通过追踪全trans伸直链长度(d tt)、成核原子维诺体积(V )和中心对称参数(S)等序参量在不同取向度熔体下的等温成核与结晶演化过程, 给出了PVA熔体成核路径及形成半晶态的分子构象; 通过对形成的半晶态高分子结构进一步分析, 发现随着熔体取向度的增加, 晶体和无定型对应的取向度也会增加, 但是当应变剪切大于5时, 其对应的结晶度、晶体和无定型取向度不再发生变化; 通过对无定型区链结构的定义与分析, 可知取向度越高的熔体对应越高的Tie链数目, 随着熔体取向度的增加, Loop链的数目也会减小; 通过恒速拉伸应力测试可知, 所形成半晶态高分子结构力学响应会随着取向度及Tie链数增加而增加, 当取向及Tie链数目相同时, 应力-应变曲线形状大小也基本保持一致.

English Abstract

    • 高分子材料成型加工过程中, 不可避免地要经历不同维度和强度流场的作用[1-3], 而这些流场的作用对于高分子材料最后的力学、光学和电学等物理特性都会产生一定的影响. 归根结底, 影响这些物理特性的本质原因是高分子熔体的结构发生了变化, 这就需要在分子层面对其结构变化的机理进行深入地研究. 在高分子制品加工初期对高分子熔体进行预取向是高分子加工的一般步骤, 如流延、挤出或者纺丝, 都是将高分子熔体进行预取向而进行初步加工的, 但是关于高分子预取向以后高分子熔体的成核与结晶特性的研究又非常有限, 因为实验在分子层面也很难解释其中的机理, 所以分子动力学模拟不失为一种有效可行的办法来解决相关的问题[4-15].

      实验上对于相关的研究已经发展了更高时间分辨率和更小空间的尺寸探测技术, 如宽角X射线散射(wide angle X-ray scattering, WAXS)、小角X射线散射(small angle X-ray scattering, SAXS)及小角中子散射(small angle neutron scattering, SANS)等探测技术[2,16,17], 其目的就是为了探索高分子材料在加工过程中结构转变的机理, 以及各结构参数之间的相互关系对加工后高分子材料最终物理特性的影响机制. 虽然目前实验技术已经有了很大的改进, 可以根据结构的原位X射线衍射花样来判断其对应的结构与性能的关系, 但是现在的实验技术还是很难从分子层面或者原子层面直接地检测到高分子成核与结晶的相关信息. 故而, 分子层面动力学模拟研究成核与结晶演化过程是十分必要的.

      目前理论模拟方面也有相关研究, 但是基本上都是基于各向同性的熔体进行的, 对于预取向熔体从无序到有序的成核与结晶机理的研究还十分有限. 关于聚合物各向同性的聚合物熔体的研究, 最近Luo和Sommer[18]的研究工作认为高分子的成核与结晶是由高分子最开始的缠结状态决定的. 关于预取向的聚合物熔体, Tang等[8]最近的模拟工作显示, 预取向会促进高密度区的形成, 从而促进晶体成核速率增加, 这个过程有别于各向同性聚合物熔体的结晶过程. Tang等[8]研究了温度带来的影响, 但是没有考虑不同取向度对无定型区和力学特性的影响, 这也是进一步研究的重点. 关于高分子成核与结晶的研究进展, 中国科学技术大学李良彬教授课题组[19]对于相关的工作也做了展望, 他们提出目前高分子独有的柔性和连接性是高分子有别于小分子的区别, 同时也指出了目前高分子结晶面临的问题: 1)柔性链段是如何变成刚性链段的?2)片晶间无定型是如何形成的?3)高分子链是如何到达晶体生长前端的?从而可以看出目前关于高分子的成核与结晶确实需要从分子和理论层面深入分析和解答加工过程中面对的这些问题.

      基于目前理论和实验方面的迫切需求, 本文将对各向同性高分子熔体的静态结晶及不同预取向的高分子熔体采用分子动力学的方法, 通过不同的序参量来表征高分子熔体在分子层面上对应的成核与结晶过程, 并比较分析给出二者之间的成核机理, 从而给出不同高分子熔体对应的成核结构及晶体生长过程对应的分子图像. 通过对形成的半晶态高分子结构进一步分析, 给出不同取向度半晶高分子结构的无定型结构特点和应力-应变响应的分子特性.

    • 本文主要以Meyer 和Müller-Plathe[20]发展的粗粒化聚乙烯醇熔体(coarse grain-polyvinyl alcohol, CG-PVA)模型为基础, 采用分子动力学方法研究CG-PVA的成核与结晶过程, 模拟软件采用平行的LAMMPS上的分子动力学进行. 对于所研究的体系, 粒子之间的非键相互作用采用Lennard-Jones-9-6 (LJ96)势来描述[21]:

      ${U_{{\rm{LJ}}}}(r) = \varepsilon \left[ {{{\left( {\dfrac{{{\sigma _0}}}{r}} \right)}^9} - {{\left( {\dfrac{{{\sigma _0}}}{r}} \right)}^6}} \right],$

      其中, r是两个粒子之间的距离, $\varepsilon {\rm{ = }}1.511{k_{\rm{B}}}T$, ${\sigma _0}{\rm{ = }}0.89\sigma $. 当${U_{{\rm{LJ}}}}(r)$最小时, 截断半径为1.02$\sigma $ (σε分别为无量纲的距离和能量常数). 两单体相邻之间的键长势由下式给出:

      ${U_{{\rm{bond}}}}(r) = \frac{1}{2}{k_{\rm{b}}}{(r - {b_0})^2},$

      其中, b0 = 0.5σ, kb = 2704kBT/σ2. CG-PVA模型的键角势为列表势, 具体参数是由Luo和Sommer[21]给出.

    • 对于研究体系, 本文选取的是具有缠结效应的N = 500的CG-PVA模型结构, 原因是将聚乙烯醇(polyvinyl alcohol, PVA)熔体从熔点以上降温到熔点以下时, 它具有更规整的结构, 这样的过程更有利于研究聚合物的结晶机理, 而且PVA可以代表类似PE的规整高分子结构的结晶过程. 为了避免边界效应的影响, 之后将1000根PVA分子采用无归行走的方式添加到周期性盒子里面, 盒子包含了50万个粒子, 其大小为61.2σ × 61.2σ × 61.2σ. 为了保证体系的合理性, 在高温T = 1.0下进行107步的时间弛豫. 高温弛豫之后, 将体系沿xy方向施加一个$1 \times {10^{{\rm{ - }}5}}{\tau ^{ - 1}}$的剪切流场(图1(a)). 在动力学弛豫和非平衡动力学过程, 本文选取的时间步长都是0.01$\tau $. 图1(b)给出了一个典型的高分子熔体在剪切场下的应力-应变曲线, 里面包含了应力上升、应力屈服、应力下降和应力平台, 结果与实验观察一致[22], 证明预剪切过程是合理的. 同时, 以剪切应变$\gamma $= 0—20对应的结构为等温结晶的研究对象, 探索不同预取向对于高分子等温成核与结晶的影响. 图1(b)还给出了剪切过程中, 链的末端距Rete及分子链沿着剪切方向的取向参数Pall的变化. 曲线参数主要根据链末端矢量与剪切方向夹角θij采用如下所示的经典取向参数定义:

      图  1  (a)高分子熔体剪切示意图; (b)高分子熔体在剪切场温度为T = 1.0下的σxy, RetePall演化曲线

      Figure 1.  (a) Schematic diagram of polymer melts under shear flow; (b) evolutions of σxy, Rete and Pall for polymer melts under shear with temperature of T = 1.0.

      $P =(1/2)\cdot{{3\left\langle {{{\cos }^2}\theta _{ij}^{} - 1} \right\rangle }} ,$

      其中θij角对应的是链与链末端矢量间的夹角大小. 不难发现, 链的末端距随着剪切应变的增加达到6个应变以后保持不变, 取向达到10个应变之后不再发生变化, 为进一步研究取向对于结晶的影响创造了条件. 在等温结晶过程中选取的结晶温度为T = 0.8, 处于结晶温度以下, 高分子熔体可以自发地成核与结晶在这个温度.

    • 有效地判断聚合物结晶过程中链内和链间的协同作用, 可以采用Luo和Sommer[21]的伸直链段判断晶体的方法, 如图2(b)所示(定义θ大于150°为trans构象, 如果trans连续单体数目大于13, 则其定义为有效伸直刚性链段), 但是并不是所有的伸直刚性链结构都是晶体结构, 故需要一个判定链间结构的方法来进一步判定. 本文在定义晶体与成核原子过程中引入了中心对称参数(centro-symmetry parameter, CSP)与伸直链判定相结合的方法来区分晶体[23], CSP的原理公式为

      图  2  (a)成核原子判定方法示意图; (b)伸直链判定方法示意图; (c) PVA半晶态高分子结构晶区与无定型区S值的分布

      Figure 2.  Schematic diagram of (a) centro-symmetry parameter and (b) stretched chain segment for PVA crystal; (c) S distribution of crystal and melt of PVA semicrystalline polymers.

      $S = \sum\limits_{i = 1}^{N/2} {\left| {{R_i} + {R_{{N}/{2} + i}}} \right|} .$

      原理图如图2(a)所示. 该方法已经被成熟地应用到了没有连接性的粒子体系中, 可以有效地判断这些粒子体系从无序到局部有序的过程. 将该方法引入到聚合物体系的成核与结晶过程, 可有效判断半晶态高分子晶区与无定型区的成核原子S值的两个峰值, 如图2(c)所示. 将这两相结合, 不难发现通过空间对称分布给出成核与结晶过程, 不只是包含了链内的构象变化, 还考虑了空间链间构象分布的影响.

    • 为了探索预取向对于结晶过程中链内构象、链间构象及密度的影响, 图3(a)(e)分别给出了不同预取向PVA熔体的结晶度${\phi _{\rm{c}}}$、伸直链段长度d tt、成核原子平均维诺体积V和平均对称参数S随时间的演化过程. 其中结晶度的计算方式与Luo和Sommer[21]定义晶体的方式是一致的, 将二面角大于150°时的结构定义为trans结构, 然后这样的全trans链段必须大于13个单体才能定义有效伸直链段, 将此判断方法与对称参数相结合, 如果S < 0.8, 将此类原子定义为成核原子.

      图  3  不同应变下PVA熔体的${\phi _{\rm{c}}}$, d tt, VS随时间的演化过程 (a) γ = 0; (b) γ = 1; (c) γ = 1.5; (d) γ = 2; (e) γ = 4; (f) γ = 8

      Figure 3.  Evolutions of ${\phi _{\rm{c}}}$, d tt, V and S of nucleation atoms under different stains: (a) γ = 0; (b) γ = 1; (c) γ = 1.5; (d) γ = 2; (e) γ = 4; (f) γ = 8.

      通过图3(a)(f)${\phi _{\rm{c}}}$的变化不难看出, 预取向的确可以促进高分子熔体的结晶过程, 而且随着剪切应变越大, 对应的成核诱导期也会越短. dtt代表链内构象的变化, 对于各向同性的聚合物熔体(图3(a)), 比较dtt${\phi _{\rm{c}}}$的变化不难发现, 伸直链长度开始明显增长的时间明显要比结晶度开始快速增长的时间要提前, 说明在均相成核与结晶过程中, 伸直链段的形成是成核前提, 而且随着伸直链段的形成, V先于伸直链段和中心对称参数降低, 密度降低之后, S也会进一步降低, 伸直链长度会发生一个突然的增加, 进一步说明在稳定核开始形成过程中, 成核链段需要通过成核原子密度的波动形成高密度区, 进而不断地调整自己的构象(S值的降低), 以促进稳定晶核的形成. 稳定晶核形成以后伴随的就是晶体稳定的生长过程(图3(a)). 对于预取向PVA熔体, 不难发现拉伸会产生更多的高密度区, 而这些高密度区会促进伸直链的进一步形成, 从而促进稳定核的形成. 与各向同性的熔体比较, 预取向主要是促进了链间的密度有序, 密度有序会进一步促进晶核的形成. 通过链内构象、链间密度及局部序参量, 可以构建出不同取向度高分子熔体的成核与结晶生长机理. 这个结论恰好验证了实验上通过SAXS或者WAXS研究结晶的一些猜想. Wang等[24]通过SAXS和WAXS对于高密度聚乙烯的研究表明, 预有序体可以促使晶核的形成, 而且预有序体属于高密度相, 这个结论与本文的模拟数据也是相符的.

    • 高分子在熔体阶段预取向的不同会导致结晶后力学性能和加工性能存在差异, 然而, 这些加工特性和物理特性的区别主要归因于产生的半晶态高分子晶区与无定型区的结构在预剪切过程中对于高分子成核的影响. 具体发生了那些改变?实验很难准确地检测到晶体和无定型结构的变化, 动力学结果的统计与分析恰好可以弥补实验上在这方面的缺陷和不足. 为了研究无定型结构的区别, 根据实验理解与可能存在的结构, 我们定义了三种类型的无定型结构(图4(a)). Tie 链是连接两个片晶之间的无定型链段, Loop链是折叠链结构, Tail链是链端结构. 从图4(c)可以看出, 本文判定方法可以准确有效地判定出晶体和无定型链的结构分布, 为进一步的研究提供了可靠的理论依据.

      图  4  (a)半晶态高分子结构中无定型对应的链结构分类模型; (b), (c) PVA半晶态高分子结构中晶体和无定型的原子结构

      Figure 4.  (a) Schematic diagram of crystalline and amorphous chain structure for the semicrystalline polymers; (b) and (c) atomic structure of crystal and melt of PVA semicrystalline polymers.

      根据图4所示的结构定义, 图5(a)给出了不同剪切应变下半晶态高分子结构晶区与无定型结构晶体取向和无定型取向的变化. 从图5(a)中不难发现, 应变越大的熔体形成的半晶态高分子结构对应的结晶度越大, 同时晶体和无定型的取向也越大. 分析不同的无定型链结构在半晶态高分子结构中的数量发现, 应变越大的熔体的Tie链数目也越多, 而Loop链的数目会随着应变的增加而减小. 在小应变下, Tie链数目越小, Loop链数目越多表明, 没有预取向的半晶态高分子晶区主要是以链内的折叠为主; 在高应变下, 主要是以链间的有序为主, 通过Tie链连接片晶与片晶. Tie链数目的改变及Loop链的数目进而会影响到力学性能的改变.

      图  5  (a)不同剪切应下PVA半晶态高分子对应的${\phi _{\rm{c}}}$、晶体取向参数Pc和无定型结构取向参数Pa; (b)不同剪切应变下PVA半晶态高分子无定型链结构数目的演化

      Figure 5.  (a) ${\phi _{\rm{c}}}$, crystalline order parameter Pc, and amorphous order parameter Pa for PVA semicrystalline polymers with different shear strains; (b) the evolution of the numbers of amorphous chains for PVA semicrystalline polymers with different shear strains.

    • 为了进一步探索不同取向度半晶态高分子结构的力学特性, 对形成的半晶态高分子结构施加了一个拉伸速率为$1 \times {10^{{\rm{ - }}5}}{\tau ^{ - 1}}$的恒应变拉伸场. 相同的拉伸速率下, 不难发现不同取向度半晶态高分子结构对应的力学性能会有很大的差别. 图6为不同初始剪切应变下体系的应力-应变曲线汇总. 拉伸的应力-应变曲线与最近实验的结果一致[25-28], 进一步从侧面反映了模拟结果的合理性, 但是实验目前只能从拉伸过程中X射线衍射花样来判断晶体的熔融重结晶过程. 相比于实验动力学的轨迹, 分子动力学可以直接给出结构的演化, 从而进一步验证实验的猜想和推断. 分析应力-应变曲线, 发现随着初始剪切应变的增加, 体系的屈服应变逐渐增加, 且初始模量也有所增长, 但是在初始应变达到一定数值后, 屈服应变和初始模量的变化都较小. 各体系的初始模量取决于各晶体的取向程度, 取向度越大, 体系初始模量越大; 通过进一步的构象分析可知, 在早期拉伸过程中, 晶体的旋转和破坏是重要的应力来源, 随着晶体与拉伸方向取向一致, 晶体破坏减少, 此时主要发生取向链段结晶和晶体间的无定形链拉伸, 因此在大应变初始预剪切体系中, 晶体取向较快, 力学曲线进而与无定形链末端距变化近乎一致; 从微观的原子受力角度来看, 晶区原子沿取向方向的受力大于非晶区原子, 整体上随着链段取向的增加, 原子沿取向方向受力也随之增加; 预剪切对缠结密度影响不大, 但是对于拉伸过程的缠结度变化率有较大影响, 整体上缠结度随着拉伸过程下降, 且初始预剪切应变越大体系最终缠结度越小.

      图  6  不同预取向半晶态高分子结构在恒定速率为$1 \times {10^{{\rm{ - }}5}}{\tau ^{ - 1}}$拉伸场下对应的应力-应变曲线

      Figure 6.  Stress-strain curves of PVA semicrystalline polymer with different orientation degree under the stretched rate of $1 \times {10^{{\rm{ - }}5}}{\tau ^{ - 1}}$.

      进一步分析单链构象的变化. 图7给出了应变等于6时形成的半静态高分子单链在不用拉伸应变下的结构. 力学曲线结合构象的演化可以得到以下结论. 1)预取向可以提高沿取向方向的力学强度, 大应变预剪切和小应变预剪切的力学曲线形状存在差异. 2) 拉伸前期, 应力来源于晶体的旋转和破坏, 之后是链段的拉伸和取向结晶. 3)链段沿拉伸方向力与取向相关性很强, 取向越大的链, 原子受力越大. 因为晶区有着更大的取向, 在拉伸方向力更大. 4)对整个拉伸过程的理解: 首先是键角尺度上的变化引起的弹性阶段, 然后是晶体旋转和破坏同时伴随无定形链拉伸取向结晶, 之后是晶体基本取向并且在拉伸下继续破坏, 随着链的取向进一步加强取向结晶数目增多, 反映为总结晶度开始上升.

      图  7  当γ = 6时对应的单链构象随着应变ε的演化过程

      Figure 7.  Structural evolution of single PVA chain as a function of strain ε when γ = 6.

    • 通过分子动力学模拟方法探究了各向同性与不同预取向PVA熔体的成核路径. 通过比较${\phi _{\rm{c}}}$, d tt, VS随着时间的演化过程, 发现取向度越高的PVA熔体对应越快的成核动力学. 并且得到以下更详细的成核动力学路径和力学响应特性: 1)对于各向同性的PVA熔体, 链段结构由柔性向刚性转变先于链间有序的发生, 进而是局部有序的发生, 最后才是成核或密度的转变发生, 这个分子图像被定义成链内有序到局部有序结构波动诱导的成核. 2)当PVA熔体在剪切场下取向后, 发现在成核结晶前发生了由于构象有序链段的聚集、共同取向而形成的高密度区, 及密度波动先于成核结晶的发生; 这样的图像被总结为流动下密度波动诱导的成核结晶. 3)预取向可以改变半晶态高分子结构的晶体和无定型的取向. 对于无定型区, 剪切应变不同, 对应的Tie链、Loop链的数目也会发生明显的变化, 进而影响其结晶加工过程的力学特性. 4)预取向可以提高沿取向方向的力学强度, 大应变预剪切和小应变预剪切的力学曲线形状存在差异. 拉伸前期, 应力来源有晶体的旋转和破坏, 之后是链段的拉伸和取向结晶. 链段沿拉伸方向力与取向相关性很强, 取向越大的链, 原子受力越大.

参考文献 (28)

目录

    /

    返回文章
    返回