搜索

x

留言板

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

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

石墨烯/柔性基底复合结构双向界面切应力传递问题的理论研究

白家豪 郭建刚

引用本文:
Citation:

石墨烯/柔性基底复合结构双向界面切应力传递问题的理论研究

白家豪, 郭建刚

Theoretical studies on bidirectional interfacial shear stress transfer of graphene/flexible substrate composite structure

Bai Jia-Hao, Guo Jian-Gang
PDF
HTML
导出引用
  • 界面力学性能是影响石墨烯/柔性基底复合结构整体力学性能的关键因素, 因此对该结构界面切应力传递机理的研究十分必要. 考虑了石墨烯和基底泊松效应的影响, 本文提出了二维非线性剪滞模型. 对于基底泊松比相比石墨烯较大的情况, 利用该模型理论研究了受单轴拉伸石墨烯/柔性基底结构的双向界面切应力传递问题. 在弹性粘结阶段, 导出了石墨烯双向正应变和双向界面切应力的半解析表达式, 分析了不同位置处石墨烯正应变和界面切应力的分布规律. 导出了石墨烯/柔性基底结构发生界面滑移的临界应变, 结果表明该临界应变低于利用经典一维非线性剪滞模型得到的滑移临界应变, 并且明显受到石墨烯宽度尺寸以及基底泊松比大小的影响. 基于二维非线性剪滞模型建立有限元模型 (FEM), 研究了界面滑移阶段石墨烯双向正应变和双向界面切应力的分布规律. 与一维非线性剪滞模型的结果对比表明, 当石墨烯宽度较大时, 二维模型和一维模型对石墨烯正应变、界面切应力以及滑移临界应变的计算结果均存在较大差别, 但石墨烯宽度很小时, 二维模型可近似被一维模型代替. 最后, 通过与拉曼实验结果的对比, 验证了二维非线性剪滞模型的可靠性, 并得到了石墨烯/聚对苯二甲酸乙二醇酯 (PET) 基底结构的界面刚度 (100 TPa/m) 和界面剪切强度 (0.295 MPa).
    Interfacial mechanical properties have a great influence on the overall mechanical performance of graphene/flexible substrate composite structure. Therefore, it is necessary to study interfacial shear stress transfer between graphene and flexible substrate. In this paper, a two-dimensional nonlinear shear-lag model (2D model) is presented. Taking the effects of Poisson’s ratio of the graphene and substrate into consideration, the bidirectional interfacial shear stress transfer between graphene and flexible substrate subjected to uniaxial tension is investigated by the 2D model when the Poisson’s ratio of substrate is larger than that of graphene. In the elastic bonding stage, the semi-analytical solutions of the bidirectional normal strains of the graphene and bidirectional interfacial shear stresses are derived, respectively, and their distributions at different positions are illustrated. The critical strain for interfacial sliding is derived by the 2D model, and the results show that the critical strain has a micron-scaled characteristic width. The width size of graphene has a significant influence on the critical strain when it is less than the characteristic width, but the size effect can be ignored when the width of graphene is larger than the characteristic width. In addition, the Poisson’s ratio of substrate can also affect the critical strain. Based on the 2D model, the finite element simulations are made to investigate the distribution of graphene's normal strains and interfacial shear stresses in the interfacial sliding stage. Furthermore, compared with the results obtained via one-dimensional nonlinear shear-lag model (1D model), the distributions of graphene’s normal strains and interfacial shear stresses calculated by 2D model show obvious bidimensional effects both in the elastic bonding stage and in the interfacial sliding stage when the width of graphene is large. There exists a compression strain in the graphene and a transverse (perpendicular to the tensile direction) shear stress in the interface, which are neglected in the 1D model. And the distributions of graphene’s tensile strain and longitudinal (along the tensile direction) interfacial shear stress are not uniform along the width, which are also significantly different from the results of 1D model. Moreover, the critical strain for interfacial sliding derived by the 2D model is lower than that obtained by the 1D model. However, when the width of graphene is small enough, the 2D model can be approximately replaced by the 1D model. Finally, by fitting the Raman experimental results, the reliability of the 2D model is verified, and the interfacial stiffness (100 TPa/m) and shear strength (0.295 MPa) between graphene and polyethylene terephthalate (PET) substrate are calculated.
      通信作者: 郭建刚, guojg@tju.edu.cn
    • 基金项目: 国家级-国家自然科学基金(11872268)
      Corresponding author: Guo Jian-Gang, guojg@tju.edu.cn
    [1]

    Geim A K 2009 Science 324 1530Google Scholar

    [2]

    Bolotin K I, Sikes K J, Jiang Z, Klima M, Fudenberg G, Hone J, Kim P, Stormer H L 2008 Solid State Commun. 146 351Google Scholar

    [3]

    Lee C G, Wei X D, Kysar J W, Hone J 2008 Science 321 385Google Scholar

    [4]

    Li X, Zhang R J, Yu W J, Wang K L, Wei J Q, Wu D H, Cao A Y, Li Z H, Cheng Y, Zheng Q S, Ruoff R S, Zhu H W 2012 Sci. Rep. 2 870Google Scholar

    [5]

    Young R J, Kinloch I A, Gong L, Novoselov K S 2012 Compos. Sci. Technol. 72 1459Google Scholar

    [6]

    Gong L, Kinloch I A, Young R J, Riaz I, Jalil R, Novoselov K S 2010 Adv. Mater. 22 2694Google Scholar

    [7]

    Jiang T, Huang R, Zhu Y 2014 Adv. Funct. Mater. 24 396Google Scholar

    [8]

    Xu C C, Xue T, Guo J G, Qin Q H, Wu S, Song H B, Xie H M 2015 J. Appl. Phys. 117 164301Google Scholar

    [9]

    Xu C C, Xue T, Guo J G, Kang Y L, Qiu W, Song H B, Xie H M 2015 Mater. Lett. 161 755Google Scholar

    [10]

    仇巍, 张启鹏, 李秋, 许超宸, 郭建刚 2017 物理学报 66 166801Google Scholar

    Qiu W, Zhang Q P, Li Q, Xu C C, Guo J G 2017 Acta Phys. Sin. 66 166801Google Scholar

    [11]

    Cox H L 1952 Br. J. Appl. Phys. 3 72Google Scholar

    [12]

    Guo G D, Zhu Y 2015 J. Appl. Mech. 82 031005Google Scholar

    [13]

    Cui Z, Guo J G 2016 AIP Adv. 6 125110Google Scholar

    [14]

    Zhang S L, Li J C M 2004 J. Polym. Sci., Part B: Polym. Phys. 42 260Google Scholar

    [15]

    Kurennov S S 2014 Mech. Compos. Mater. 50 105Google Scholar

    [16]

    Mathias J D, Grédiac M, Balandraud X 2006 Int. J. Solids Struct. 43 6921Google Scholar

    [17]

    Randrianalisoa J, Dendievel R, Bréchet Y 2011 Compos. Part B: Eng. 42 2055Google Scholar

    [18]

    Park K, Paulino G H 2011 Appl. Mech. Rev. 64 060802Google Scholar

    [19]

    Dourado N, Silva F G A, de Moura M F S F 2018 Constr. Build. Mater. 176 14Google Scholar

    [20]

    Högberg J L 2006 Int. J. Fract. 141 549Google Scholar

    [21]

    Camanho P P, Davila C G, de Moura M F 2003 J. Compos. Mater. 37 1415Google Scholar

    [22]

    Faccio R, Denis P A, Pardo H, Goyenola C, Mombrú A W 2009 J. Phys. Condens. Matter 21 285304Google Scholar

    [23]

    许超宸 2019 博士学位论文 (天津: 天津大学)

    Xu C C 2019 Ph. D. Dissertation (Tianjin: Tianjin University) (in Chinese)

    [24]

    Mohiuddin T M G, Lombardo A, Nair R R, Bonetti A, Savini G, Jalil R, Bonini N, Basko D M, Galiotis C, Marzari N, Novoselov K S, Geim A K, Ferrari A C 2009 Phys. Rev. B 79 205433Google Scholar

    [25]

    Sakata H, Dresselhaus G, Dresselhaus M S, Endo M 1988 J. Appl. Phys. 63 2769Google Scholar

    [26]

    Ni Z H, Yu T, Lu Y H, Wang Y Y, Feng Y P, Shen Z X 2008 ACS Nano 2 2301Google Scholar

    [27]

    Yu T, Ni Z H, Du C L, You Y M, Wang Y Y, Shen Z X 2008 J. Phys. Chem. C 112 12602

    [28]

    Koukaras E N, Androulidakis C, Anagnostopoulos G, Papagelis K, Galiotis C 2016 Extreme Mech. Lett. 8 191Google Scholar

  • 图 1  受单轴拉伸载荷的石墨烯/基底结构示意图

    Fig. 1.  Schematic diagram of the graphene/substrate structure under uniaxial tension.

    图 2  石墨烯单元的应力状态示意图

    Fig. 2.  The force balance of an element of graphene.

    图 3  局部界面上一点的切应力分析

    Fig. 3.  Analysis of interfacial shear stresses at local interface.

    图 4  (a) 二维非线性剪滞模型; (b) 双线性内聚力模型(Ⅱ + Ⅲ型混合模式)

    Fig. 4.  (a) Two-dimensional nonlinear shear-lag model; (b) bilinear cohesive shear-mode (Ⅱ + Ⅲ) law.

    图 5  弹性粘结阶段 (εsx = 0.2%) 时石墨烯正应变 (a) εx和(b) εy以及界面切应力 (c) τzx和(d) τzy的分布

    Fig. 5.  Distributions of graphene’s normal strains (a) εx and (b) εy; distributions of interfacial shear stresses (c) τzx and (d) τzy at the elastic bonding stage (εsx = 0.2%).

    图 6  典型线处石墨烯正应变 (a) εx和(b) εy以及界面切应力 (c) τzx和(d) τzy的分布

    Fig. 6.  Distributions of graphene’s strains (a) εx and (b) εy; distributions of interfacial shear stresses (c) τzx and (d) τzy along several representative lines.

    图 7  不同基底泊松比情况下滑移临界应变εsxc随石墨烯宽度W的变化(线为理论值, 散点为有限元值)

    Fig. 7.  Variation of the critical strain for sliding with the width of graphene at different Poisson's ratio of substrate (the lines are the theoretical results, and the scatter points are the FEM results).

    图 8  界面滑移阶段示意图

    Fig. 8.  Schematic diagram of interfacial sliding stage.

    图 9  界面滑移阶段 (εsx = 1%) 时石墨烯正应变 (a) εx和(b) εy以及界面切应力 (c) τzx和(d) τzy.的分布

    Fig. 9.  Distributions of graphene’s normal strains (a) εx and (b) εy; distributions of interfacial shear stresses (c) τzx and (d) τzy at the interfacial sliding stage (εsx = 1%).

    图 10  不同基底应变下石墨烯中心C点处压缩应变εyC随石墨烯宽度的变化

    Fig. 10.  Variation of compressive strain εyC at the center point C of graphene with its width when the strain of substrate is different.

    图 11  二维模型与一维模型结果的比较 (W = 21.8 μm) (a) εx和(b) τzx在弹性粘结阶段(εsx = 0.2%); (c) εx和(d) τzx在界面滑移阶段 (εsx = 1%)

    Fig. 11.  Comparisons of the results obtained via one-dimensional and two-dimensional models (W = 21.8 μm): (a) εx and (b) τzx at the elastic bonding stage (εsx = 0.2%); (c) εx and (d) τzx at the interfacial sliding stage (εsx = 1%).

    图 12  二维模型与一维模型计算结果的比较 (W = 1 μm) (a) εx和(b) τzx在界面滑移阶段 (εsx = 1%)

    Fig. 12.  Comparisons of the results obtained via one-dimensional and two-dimensional models (W = 1 μm): (a) εx and (b) τzx at the interfacial sliding stage (εsx = 1%).

    图 13  利用二维模型与实验数据拟合 (a) 基底拉伸应变εsx = 0.25%时εm沿石墨烯中心线(y = W/2)的分布; (b) 不同基底载荷作用下石墨烯中心C点处$ \varepsilon_{\rm m}^{C} $的大小

    Fig. 13.  Fitting results of experimental data by using 2D model: (a) εm along the centerline (y = W/2) when the tensile strain εsx = 0.25%; (b) $ \varepsilon_{\rm m}^{C} $at the center point C under different tensile loads.

  • [1]

    Geim A K 2009 Science 324 1530Google Scholar

    [2]

    Bolotin K I, Sikes K J, Jiang Z, Klima M, Fudenberg G, Hone J, Kim P, Stormer H L 2008 Solid State Commun. 146 351Google Scholar

    [3]

    Lee C G, Wei X D, Kysar J W, Hone J 2008 Science 321 385Google Scholar

    [4]

    Li X, Zhang R J, Yu W J, Wang K L, Wei J Q, Wu D H, Cao A Y, Li Z H, Cheng Y, Zheng Q S, Ruoff R S, Zhu H W 2012 Sci. Rep. 2 870Google Scholar

    [5]

    Young R J, Kinloch I A, Gong L, Novoselov K S 2012 Compos. Sci. Technol. 72 1459Google Scholar

    [6]

    Gong L, Kinloch I A, Young R J, Riaz I, Jalil R, Novoselov K S 2010 Adv. Mater. 22 2694Google Scholar

    [7]

    Jiang T, Huang R, Zhu Y 2014 Adv. Funct. Mater. 24 396Google Scholar

    [8]

    Xu C C, Xue T, Guo J G, Qin Q H, Wu S, Song H B, Xie H M 2015 J. Appl. Phys. 117 164301Google Scholar

    [9]

    Xu C C, Xue T, Guo J G, Kang Y L, Qiu W, Song H B, Xie H M 2015 Mater. Lett. 161 755Google Scholar

    [10]

    仇巍, 张启鹏, 李秋, 许超宸, 郭建刚 2017 物理学报 66 166801Google Scholar

    Qiu W, Zhang Q P, Li Q, Xu C C, Guo J G 2017 Acta Phys. Sin. 66 166801Google Scholar

    [11]

    Cox H L 1952 Br. J. Appl. Phys. 3 72Google Scholar

    [12]

    Guo G D, Zhu Y 2015 J. Appl. Mech. 82 031005Google Scholar

    [13]

    Cui Z, Guo J G 2016 AIP Adv. 6 125110Google Scholar

    [14]

    Zhang S L, Li J C M 2004 J. Polym. Sci., Part B: Polym. Phys. 42 260Google Scholar

    [15]

    Kurennov S S 2014 Mech. Compos. Mater. 50 105Google Scholar

    [16]

    Mathias J D, Grédiac M, Balandraud X 2006 Int. J. Solids Struct. 43 6921Google Scholar

    [17]

    Randrianalisoa J, Dendievel R, Bréchet Y 2011 Compos. Part B: Eng. 42 2055Google Scholar

    [18]

    Park K, Paulino G H 2011 Appl. Mech. Rev. 64 060802Google Scholar

    [19]

    Dourado N, Silva F G A, de Moura M F S F 2018 Constr. Build. Mater. 176 14Google Scholar

    [20]

    Högberg J L 2006 Int. J. Fract. 141 549Google Scholar

    [21]

    Camanho P P, Davila C G, de Moura M F 2003 J. Compos. Mater. 37 1415Google Scholar

    [22]

    Faccio R, Denis P A, Pardo H, Goyenola C, Mombrú A W 2009 J. Phys. Condens. Matter 21 285304Google Scholar

    [23]

    许超宸 2019 博士学位论文 (天津: 天津大学)

    Xu C C 2019 Ph. D. Dissertation (Tianjin: Tianjin University) (in Chinese)

    [24]

    Mohiuddin T M G, Lombardo A, Nair R R, Bonetti A, Savini G, Jalil R, Bonini N, Basko D M, Galiotis C, Marzari N, Novoselov K S, Geim A K, Ferrari A C 2009 Phys. Rev. B 79 205433Google Scholar

    [25]

    Sakata H, Dresselhaus G, Dresselhaus M S, Endo M 1988 J. Appl. Phys. 63 2769Google Scholar

    [26]

    Ni Z H, Yu T, Lu Y H, Wang Y Y, Feng Y P, Shen Z X 2008 ACS Nano 2 2301Google Scholar

    [27]

    Yu T, Ni Z H, Du C L, You Y M, Wang Y Y, Shen Z X 2008 J. Phys. Chem. C 112 12602

    [28]

    Koukaras E N, Androulidakis C, Anagnostopoulos G, Papagelis K, Galiotis C 2016 Extreme Mech. Lett. 8 191Google Scholar

  • [1] 刘东静, 周福, 胡志亮, 黄家强. 石墨烯/GaN异质结构界面热输运性质的分子动力学研究. 物理学报, 2024, 73(13): 137901. doi: 10.7498/aps.73.20240021
    [2] 王伟华. 二维有限元方法研究石墨烯环中磁等离激元. 物理学报, 2023, 72(8): 087301. doi: 10.7498/aps.72.20222467
    [3] 刘瑛, 郭斯琳, 张勇, 杨鹏, 吕克洪, 邱静, 刘冠军. 1/f噪声及其在二维材料石墨烯中的研究进展. 物理学报, 2023, 72(1): 017302. doi: 10.7498/aps.72.20221253
    [4] 孙颖慧, 穆丛艳, 蒋文贵, 周亮, 王荣明. 金属纳米颗粒与二维材料异质结构的界面调控和物理性质. 物理学报, 2022, 71(6): 066801. doi: 10.7498/aps.71.20211902
    [5] 徐翔, 张莹, 闫庆, 刘晶晶, 王骏, 徐新龙, 华灯鑫. 不同堆垛结构二硫化铼/石墨烯异质结的光电化学特性. 物理学报, 2021, 70(9): 098203. doi: 10.7498/aps.70.20201904
    [6] 郭丽娟, 胡吉松, 马新国, 项炬. 二硫化钨/石墨烯异质结的界面相互作用及其肖特基调控的理论研究. 物理学报, 2019, 68(9): 097101. doi: 10.7498/aps.68.20190020
    [7] 周愈之. 过渡金属硫族化合物柔性基底体系的模型与应用. 物理学报, 2018, 67(21): 218102. doi: 10.7498/aps.67.20181571
    [8] 白清顺, 沈荣琦, 何欣, 刘顺, 张飞虎, 郭永博. 纳米微结构表面与石墨烯薄膜的界面黏附特性研究. 物理学报, 2018, 67(3): 030201. doi: 10.7498/aps.67.20172153
    [9] 施佳妤, 蓝尤钊. 类石墨烯结构二维层状碳化硅的非线性二次谐波特性的第一性原理研究. 物理学报, 2018, 67(21): 217803. doi: 10.7498/aps.67.20181337
    [10] 仇巍, 张启鹏, 李秋, 许超宸, 郭建刚. 单层单晶石墨烯与柔性基底界面性能的实验研究. 物理学报, 2017, 66(16): 166801. doi: 10.7498/aps.66.166801
    [11] 蒋钊, 陈学康. 界面合金化控制柔性Al/PI薄膜应力的研究. 物理学报, 2015, 64(21): 216802. doi: 10.7498/aps.64.216802
    [12] 娄利飞, 潘青彪, 吴志华. 基于石墨烯用于微弱能量获取的柔性微结构研究. 物理学报, 2014, 63(15): 158501. doi: 10.7498/aps.63.158501
    [13] 黄坤, 殷雅俊, 吴继业. 单层石墨烯片的非线性板模型. 物理学报, 2014, 63(15): 156201. doi: 10.7498/aps.63.156201
    [14] 韩文鹏, 史衍猛, 李晓莉, 罗师强, 鲁妍, 谭平恒. 石墨烯等二维原子晶体薄片样品的光学衬度计算及其层数表征. 物理学报, 2013, 62(11): 110702. doi: 10.7498/aps.62.110702
    [15] 贾汝娟, 王苍龙, 杨阳, 苟学强, 陈建敏, 段文山. 二维Frenkel-Kontorova模型中六角对称结构的摩擦现象. 物理学报, 2013, 62(6): 068104. doi: 10.7498/aps.62.068104
    [16] 林力, 李云, 顾兆林, 刘兆杰, 程光旭. 计算二维声腔传递矩阵的正方形线声源模型. 物理学报, 2009, 58(8): 5484-5490. doi: 10.7498/aps.58.5484
    [17] 倪向贵, 殷建伟. 拉伸条件下双壁碳纳米管弹性性能的原子模拟. 物理学报, 2006, 55(12): 6522-6525. doi: 10.7498/aps.55.6522
    [18] 倪培根, 马博琴, 程丙英, 张道中. 二维LiNbO3非线性光子晶体. 物理学报, 2003, 52(8): 1925-1928. doi: 10.7498/aps.52.1925
    [19] 方前锋. 低应力振幅下非线性滞弹性内耗峰(P′1峰)的数值分析. 物理学报, 1997, 46(3): 536-543. doi: 10.7498/aps.46.536
    [20] 任尚元. 面心立方晶体的泊松比. 物理学报, 1983, 32(5): 664-669. doi: 10.7498/aps.32.664
计量
  • 文章访问数:  11581
  • PDF下载量:  259
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-11-11
  • 修回日期:  2019-12-11
  • 刊出日期:  2020-03-05

/

返回文章
返回