删除或更新信息,请邮件至freekaoyan#163.com(#换成@)

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

本站小编 Free考研考试/2021-12-29

摘要:界面力学性能是影响石墨烯/柔性基底复合结构整体力学性能的关键因素, 因此对该结构界面切应力传递机理的研究十分必要. 考虑了石墨烯和基底泊松效应的影响, 本文提出了二维非线性剪滞模型. 对于基底泊松比相比石墨烯较大的情况, 利用该模型理论研究了受单轴拉伸石墨烯/柔性基底结构的双向界面切应力传递问题. 在弹性粘结阶段, 导出了石墨烯双向正应变和双向界面切应力的半解析表达式, 分析了不同位置处石墨烯正应变和界面切应力的分布规律. 导出了石墨烯/柔性基底结构发生界面滑移的临界应变, 结果表明该临界应变低于利用经典一维非线性剪滞模型得到的滑移临界应变, 并且明显受到石墨烯宽度尺寸以及基底泊松比大小的影响. 基于二维非线性剪滞模型建立有限元模型 (FEM), 研究了界面滑移阶段石墨烯双向正应变和双向界面切应力的分布规律. 与一维非线性剪滞模型的结果对比表明, 当石墨烯宽度较大时, 二维模型和一维模型对石墨烯正应变、界面切应力以及滑移临界应变的计算结果均存在较大差别, 但石墨烯宽度很小时, 二维模型可近似被一维模型代替. 最后, 通过与拉曼实验结果的对比, 验证了二维非线性剪滞模型的可靠性, 并得到了石墨烯/聚对苯二甲酸乙二醇酯 (PET) 基底结构的界面刚度 (100 TPa/m) 和界面剪切强度 (0.295 MPa).
关键词: 石墨烯/柔性基底结构/
泊松比/
界面切应力传递/
二维非线性剪滞模型

English Abstract


--> --> -->
石墨烯是一种具有单原子层厚度的二维纳米材料[1], 由于具有独特的电学、热学和力学性能[2,3], 被认为在柔性电子器件和纳米复合材料等领域具有广阔的应用前景[4,5]. 在应用于柔性电子器件和纳米复合材料中时, 石墨烯通常被置于柔性基底上, 或嵌入基质中. 与碳纳米管/基体结构类似, 决定石墨烯/柔性基底复合结构的应用性能和可靠性的最关键问题之一在于其界面力学性能, 尤其是石墨烯与基底之间的剪切相互作用. 因此, 深入研究石墨烯与基底间界面切应力传递问题具有十分重要的意义.
近年来, 已有众多****对石墨烯与柔性基底间的界面切应力传递问题进行了研究. 在实验方面, Gong等[6]用拉曼实验技术对SU-8光刻胶/石墨烯/聚甲基丙烯酸甲酯 (PMMA)“三明治”结构界面的切应力传递进行了监测, 并且观察到与受较大应变载荷时碳纤维和基体间界面的断裂行为不同, 石墨烯与聚合物基底间的应力传递仍然通过界面摩擦进行. 这种界面纳米摩擦行为在Jiang等[7]的实验中也被观察到, Jiang等利用原位拉曼光谱技术对置于单轴拉伸基底上石墨烯的应变进行了测量, 发现当基底受较大拉伸载荷时, 界面通过纳米摩擦作用进行切应力传递. 在Jiang等的研究中将这种界面摩擦行为称为界面滑移(下文沿用这一称谓), 并且计算得到基底所受拉伸应变载荷为0.3%时, 界面边缘开始出现滑移区域. 此外, Xu等[8,9]和仇巍等[10]在更大载荷范围对石墨烯/聚对苯二甲酸乙二醇酯 (PET)基底界面应力传递进行了研究, 同样在基底受到较大拉伸载荷情况下观察到了界面滑移行为, 并且发现石墨烯与PET基底间界面应力传递存在尺寸效应问题.
在理论研究方面, Cox的剪切-滞后理论[11](剪滞理论)已经被广泛应用到受单轴载荷的石墨烯/柔性基底结构界面切应力传递问题的研究中. Gong等[6]利用一维线性剪滞模型与实验结果进行拟合对比, 计算得出SU-8光刻胶/石墨烯/PMMA“三明治”结构的界面剪切强度约为2.3 MPa. Jiang等[7]在线性剪滞模型基础上, 发展出了考虑界面滑移的一维非线性剪滞模型. 模型假设界面在弹性粘结阶段之后, 进入滑移阶段, 界面切应力大小随界面相对位移的增加保持不变. 利用该理论模型拟合实验数据, 获得石墨烯与PET基底间界面剪切强度为0.46—0.69 MPa. 另外, Guo 和Zhu[12]以及Cui和Guo[13]在理论研究中, 将受单轴载荷的石墨烯/基底结构沿加载方向上的单向界面切应力传递视为Ⅱ型裂纹沿界面切向的扩展, 并利用双线性内聚力模型近似代替一维非线性剪滞模型对该问题进行了研究, 计算了复合结构界面滑移的临界应变, 结果显示滑移临界应变明显受到石墨烯长度尺度的影响. 但是, 常用来作为石墨烯复合结构基底的柔性材料(PET, PMMA等)往往具有比石墨烯更大的泊松比[6-10]. 然而一维剪滞模型只考虑了受单轴载荷的石墨烯/基底结构沿加载方向的界面切应力传递, 忽略了石墨烯和基底泊松效应引起的界面横向切应力传递作用, 可能会对石墨烯应变分布、界面切应力分布以及滑移临界应变的计算结果造成一定的影响.
本文在考虑沿加载方向界面切应力传递基础上, 考虑了石墨烯和基底泊松效应引起的界面横向切应力传递作用, 提出了可用于研究石墨烯/柔性基底结构双向界面切应力传递问题的二维非线性剪滞模型. 对于基底泊松比相比石墨烯较大的情况, 利用该模型理论研究了受单轴拉伸石墨烯/柔性基底复合结构的双向界面切应力传递问题. 推导了弹性粘结阶段石墨烯双向正应变、双向界面切应力的半解析表达式以及滑移临界应变的解析表达式. 基于二维非线性剪滞模型建立有限元模型, 计算了界面滑移阶段石墨烯双向正应变和双向界面切应力的分布. 然后, 与一维非线性剪滞模型进行比较, 分析两种计算模型的区别. 最后采用二维非线性剪滞模型与实验数据进行拟合, 得到石墨烯与PET基底间界面参数, 验证模型的可靠性.
2
2.1.问题描述
-->一块长L、宽W的石墨烯薄膜完美地粘附在受单轴拉伸载荷的柔性基底上, 如图1所示, 以O为原点, 建立坐标系. 沿x轴方向, 由于界面切应力作用, 基底的拉伸应变可以引起石墨烯薄膜的拉伸. 沿y轴方向, 由于泊松效应, 石墨烯与基底产生横向压缩应变, 当基底泊松比相比石墨烯较大时(例如, PET泊松比为0.30—0.35[14]), 基底压缩应变同样会通过界面横向切应力(垂直于拉伸方向)传递到石墨烯薄膜上. 所以, 石墨烯与受单轴拉伸载荷的基底间界面上实际存在双向切应力传递作用, 石墨烯实际处于二维应力状态下.
图 1 受单轴拉伸载荷的石墨烯/基底结构示意图
Figure1. Schematic diagram of the graphene/substrate structure under uniaxial tension.

模型中, 石墨烯和基底均被视为各向同性线弹性连续体. 设石墨烯杨氏模量为E, 泊松比为υ, 厚度为h. 基底杨氏模量为Es, 泊松比为υs. 石墨烯的本构方程为:
$\left\{ \begin{aligned}&\frac{{\partial u}}{{\partial x}} = {\varepsilon _x} = \frac{1}{E}\left( {{\sigma _x} - \upsilon {\sigma _y}} \right),\\& \frac{{\partial v}}{{\partial y}} = {\varepsilon _y} = \frac{1}{E}\left( {{\sigma _y} - \upsilon {\sigma _x}} \right),\end{aligned} \right.$
式中u, v, σx, σy, εxεy分别代表了石墨烯沿xy方向上的位移、正应力和正应变. 同时, 基底的位移与正应变有如下几何关系:
$\left\{ \begin{aligned}&\frac{{\partial {u_{\rm{s}}}}}{{\partial x}} = {\varepsilon _{{\rm{s}}x}},\\&\frac{{\partial {v_{\rm{s}}}}}{{\partial y}} = {\varepsilon _{{\rm{s}}y}},\end{aligned} \right.$
式中us, vs, εsxεsy分别代表了基底沿xy方向上的位移和正应变, 其中${\varepsilon _{{\rm{s}}y}} = - {\upsilon _{\rm{s}}} \times {\varepsilon _{{\rm{s}}x}}$. 石墨烯与基底沿界面切向的相对位移为:
$\left\{ \begin{aligned}& {\delta _x} = {u_{\rm{s}}} - u,\\& {\delta _y} = {v_{\rm{s}}} - v,\end{aligned} \right.$
式中δxδy分别代表了石墨烯与基底沿xy两个切向方向的相对位移.
对石墨烯单元进行受力分析, 如图2所示. 仿照与石墨烯/基底结构相类似的矩形胶结补片结构与胶结接头结构求解双向应力解时的经典假设[15-17], 忽略石墨烯薄膜面内切应力, 从而得到平衡方程:
图 2 石墨烯单元的应力状态示意图
Figure2. The force balance of an element of graphene.

$\left\{ \begin{aligned}& h\frac{{\partial {\sigma _x}}}{{\partial x}} = - {\tau _{zx}},\\& h\frac{{\partial {\sigma _y}}}{{\partial y}} = - {\tau _{zy}},\end{aligned} \right.$
其中τzx, τzy分别代表界面沿xy方向的切应力.
2
2.2.二维非线性剪滞模型
-->一维非线性剪滞模型[7]对石墨烯/基底结构单向界面切应力传递问题进行了描述. 界面局部处于弹性粘结阶段时, 切应力τziK为界面刚度随基底与石墨烯界面相对位移δi的增加而线性增加, 当切应力达到界面剪切强度τc后, 界面进入滑移损伤阶段, 切应力大小保持不变. 本构关系为:
${\tau _{z{\rm{i}}}} = \left\{ {\begin{aligned}&K{\delta _{\rm{i}}}, \quad\quad\quad\quad\quad\quad\; {\; 0 \leqslant {\delta _{\rm{i}}} \leqslant {\delta _{\rm{c}}},}\\&\left( {1 - d} \right)K{\delta _{\rm{i}}} = {\tau _{\rm{c}}}, \quad {{\delta _{\rm{i}}} \geqslant {\delta _{\rm{c}}}\; ,\; }\end{aligned}} \right.$
其中损伤变量$d = 1 - {{{\delta _{\rm{c}}}} / {{\delta _{\rm{i}}}}}$, δc为滑移起始时界面的相对位移.
值得注意的是, 内聚力模型与非线性剪滞模型描述的物理概念一致, 均在描述牵引力和界面相对位移的关系, 非线性剪滞模型本质上相当于具有无限大断裂能的双线性内聚力模型. 在复合结构界面分层等问题的研究中, 众多****在描述单一断裂模式(Ⅰ, Ⅱ, Ⅲ型裂纹)的内聚力模型基础上, 发展出了可用于研究混合裂纹扩展的二维与三维内聚力模型[18-20]. Guo和Zhu[12]以及Cui和Guo[13]将石墨烯/基底单向界面切应力传递问题视为Ⅱ型裂纹扩展. 本文将石墨烯/基底双向界面切应力传递问题视为Ⅱ + Ⅲ型混合裂纹沿界面切向的扩展, 并仿照内聚力模型的做法, 将一维非线性剪滞模型(以下简称一维模型)扩展为可用于研究双向界面切应力传递的二维非线性剪滞模型(以下简称二维模型).
石墨烯/基底结构局部界面上一点的切应力分析如图3所示, 假设界面沿切向各向同性. 界面切应力合力τm
图 3 局部界面上一点的切应力分析
Figure3. Analysis of interfacial shear stresses at local interface.

${\tau _{\rm{m}}} = \sqrt {{\tau _{zx}}^2 + {\tau _{zy}}^2},$
界面实际相对位移δm
${\delta _{\rm{m}}} = \sqrt {{\delta _x}^2 + {\delta _y}^2} .$
二维模型本构关系如图4(a)所示, 当界面处于弹性粘结阶段时, 界面切应力合力形式本构关系为${\tau _{\rm{m}}} = K{\delta _{\rm{m}}}$, 分量形式本构关系为:
图 4 (a) 二维非线性剪滞模型; (b) 双线性内聚力模型(Ⅱ + Ⅲ型混合模式)
Figure4. (a) Two-dimensional nonlinear shear-lag model; (b) bilinear cohesive shear-mode (Ⅱ + Ⅲ) law.

$ {\tau _{zx}} = K{\delta _x}, ~~ {\tau _{zy}} = K{\delta _y},$
界面剪切强度为τc, 界面滑移起始时, 切应力合力τm
${\tau _{\rm{m}}} = \sqrt {{\tau _{zx}}^2 + {\tau _{zy}}^2} = {\tau _{\rm{c}}}.$
结合(7)式— (9)式, 得到滑移起始时界面实际相对位移δm, c:
$\delta _{{\rm{m,c}}}^{} = \sqrt {{\delta _x}^2 + {\delta _y}^2} = \frac{{{\tau _{\rm{c}}}}}{K}.$
进入滑移阶段时, 界面切应力合力τm大小保持不变, 方向与实际相对位移δm方向相同. 切应力合力形式本构关系为${\tau _{\rm{m}}} = \left( {1 - D} \right)K{\delta _{\rm{m}}} = {\tau _{\rm{c}}}$, 分量形式本构关系为:
$\left\{ \begin{aligned}& {\tau _{zx}} = \left( {1 - D} \right)K{\delta _x},\\& {\tau _{zy}} = \left( {1 - D} \right)K{\delta _y},\end{aligned} \right.$
其中损伤变量$D = 1 - {{{\delta _{{\rm{m, c}}}}} / {{\delta _{\rm{m}}}}}$.
Guo和Zhu[12]以及Cui和Guo[13]利用研究Ⅱ型裂纹的双线性内聚力模型代替一维模型对石墨烯/基底结构单向界面切应力传递问题进行了研究. 而对于研究Ⅱ + Ⅲ型混合裂纹的双线性内聚力模型[21](本构关系如图4(b)), 当采用Benzeggagh-Kenane准则(B-K准则)判断断裂起始时, 断裂能Gc
${G_{\rm{c}}} = {G_{{\rm{shear}}}} = \frac{1}{2}{\tau _{\rm{c}}}{\delta _{{\rm{m,f}}}},$
其中δm, f表示混合裂纹断裂起始时界面的相对位移, Gshear代表Ⅱ + Ⅲ型混合裂纹起裂的能量释放率. 当${G_{{\rm{shear}}}} \to + \infty $时, ${\delta _{{\rm{m, f}}}} \to + \infty $. 此时双线性内聚力模型的损伤变量Dcoh[21]
${D_{{\rm{coh}}}} \!=\! \mathop {\lim }\limits_{{\delta _{{\rm{m,f}} \to + \infty }}} \frac{{{\delta _{{\rm{m,f}}}}\left( {{\delta _m} - {\delta _{{\rm{m,c}}}}} \right)}}{{{\delta _{\rm{m}}}\left( {{\delta _{{\rm{m,f}}}} - {\delta _{{\rm{m,c}}}}} \right)}} \approx D \!=\! 1 \!-\! \frac{{{\delta _{{\rm{m,c}}}}}}{{{\delta _{\rm{m}}}}}.$
${D_{{\rm{coh}}}} \approx D$时, 两模型本构关系一致, 双线性内聚力模型软化阶段相当于二维模型的界面滑移阶段. 所以可以在ABAQUS中, 利用研究Ⅱ + Ⅲ型混合裂纹的双线性三维内聚力单元(沿垂直界面方向相对位移${\delta _z} = 0$)代替二维模型对石墨烯/基底结构双向界面切应力传递问题进行仿真模拟研究.
2
3.1.弹性粘结阶段
-->当基底所受拉伸载荷较小时, 石墨烯/基底界面处于弹性粘结阶段, 此时采用二维模型弹性粘结阶段本构关系描述界面力学行为. 将(8)式代入(4)式, 结合(1)式— (3)式, 得到如下偏微分方程组:
$\left\{ \begin{aligned}& \frac{{{\partial ^2}{\sigma _x}}}{{\partial {x^2}}} = A{\sigma _x} + B{\sigma _y} + C,\\& \frac{{{\partial ^2}{\sigma _y}}}{{\partial {y^2}}} = B{\sigma _x} + A{\sigma _y} + D,\end{aligned} \right.$
方程组中$A = \dfrac{K}{{Eh}}$, $B = - \dfrac{{K\upsilon }}{{Eh}}$, $C = - \dfrac{{K{\varepsilon _{{\rm{s}}x}}}}{h}$, $D = - \dfrac{{K{\varepsilon _{{\rm{s}}y}}}}{h}$. 相应边界条件如下, 沿石墨烯边缘$x = 0$, $x = L$上有${\sigma _x} = 0$
$\frac{{{\partial ^2}{\sigma _y}}}{{\partial {y^2}}} = A{\sigma _y} + D.$
对(15)式进行求解, 得到σy$x = 0$, $x = L$上边界条件为
${\sigma _y} = {\varepsilon _{{\rm{s}}y}}E\left( {1 - \cosh (\alpha y) + \frac{{\cosh (\alpha W) - 1}}{{\sinh (\alpha W)}}\sinh (\alpha y)} \right).\tag{16a}$
同理可得沿$y = 0$, $y = W$上边界条件为${\sigma _y} = 0$
${\sigma _x} = {\varepsilon _{{\rm{s}}x}}E\left( {1 - \cosh (\alpha x) + \frac{{\cosh (\alpha L) - 1}}{{\sinh (\alpha L)}}\sinh (\alpha x)} \right),\tag{16b}$
其中$\alpha = \sqrt {{K / {(Eh)}}} $. 上述耦合偏微分方程组(14)式已在文献[16]中得到了解决. 当基底受双轴拉伸载荷时(即${\varepsilon _{{\rm{s}}x}} > 0$, ${\varepsilon _{{\rm{s}}y}} > 0$), 可以直接利用傅里叶级数法计算得到方程组中正应力σxσy的半解析解, 将计算结果代入(4)式即可求得界面切应力τzxτzy的半解析解. 其结果具体形式如下:
$\begin{split}{\sigma _x} =\; & \sum\limits_{n = 1,3,5,\cdots}^\infty \left\{ \sum\limits_{m = 1,3,5,.\cdots}^\infty {{f_{mn,x}}\sin \left( {\frac{{m{\text{π}}x}}{L}} \right)}\right. \\ & \times \left.\sin \left[ {\frac{{n{\text{π}}\left( {{\varPhi _{m,x}} + y} \right)}}{{\left( {2{\varPhi _{m,x}} + W} \right)}}} \right] \right\},\end{split}\tag{17a}$
$\begin{split}{\sigma _y} =\; & \sum\limits_{n = 1,3,5,\cdots}^\infty \left\{ \sum\limits_{m = 1,3,5,\cdots}^\infty {{f_{mn,y}}\sin \left( {\frac{{m{\text{π}}y}}{W}} \right)}\right.\\ & \times \left.\sin \left[ {\frac{{n{\text{π}}\left( {{\varPhi _{m,y}} + x} \right)}}{{\left( {2{\varPhi _{m,y}} + L} \right)}}} \right] \right\},\end{split}\tag{17b}$
$\begin{split}{\tau _{zx}} = \; & - h\sum\limits_{n = 1,3,5,\cdots}^\infty \left\{\sum\limits_{m = 1,3,5,\cdots}^\infty \frac{{m{\text{π}}}}{L}{f_{mn,x}}\right.\\ & \times \left.\cos \left( {\frac{{m{\text{π}}x}}{L}} \right) \sin \left[ {\frac{{n{\text{π}}\left( {{\varPhi _{m,x}} + y} \right)}}{{\left( {2{\varPhi _{m,x}} + W} \right)}}} \right] \right\} ,\end{split}\tag{18a}$
$\begin{split}{\tau _{zy}} = \; & - h\sum\limits_{n = 1,3,5,\cdots}^\infty \left\{\sum\limits_{m = 1,3,5,\cdots}^\infty \frac{{m{\text{π}}}}{W}{f_{mn,y}}\right.\\ & \times \left.\cos \left( {\frac{{m{\text{π}}y}}{W}} \right) \sin \left[ {\frac{{n{\text{π}}\left( {{\varPhi _{m,y}} + x} \right)}}{{\left( {2{\varPhi _{m,y}} + L} \right)}}} \right] \right\},\end{split}\tag{18b}$
式中fmn,x, fmn,y, Φm,xΦm,y均为与材料尺寸和力学性质有关的参数, 表达式详见附录.
对于本文研究的基底受单轴拉伸的情况(即${\varepsilon _{{\rm{s}}x}} > 0$, ${\varepsilon _{{\rm{s}}y}} = - {\upsilon _{\rm{s}}} \times {\varepsilon _{{\rm{s}}x}} < 0$), 需使用叠加法计算. 取石墨烯杨氏模量$E = 1\;{\rm{TPa}}$, 泊松比$\upsilon = 0.165$[22], 厚度$h = 0.335\; \;{\rm{nm}}$, $L = W = 21.8\;{{\text{μ}}\rm{m}}$, 基底 (PET) 泊松比${\upsilon _{\rm{s}}} = 0.33$[14]. 根据文献[7]选取石墨烯/基底 (PET) 界面刚度$K = 74\;{\rm{TPa/m}}$, 界面剪切强度${\tau _{\rm{c}}} = 0.5\;{\rm{MPa}}$. 当基底所受拉伸载荷${\varepsilon _{{\rm{s}}x}} = 0.2 $%时, 基底横向压缩应变${\varepsilon _{{\rm{s}}y}} = - {\upsilon _{\rm{s}}} \times 0.2 $%, 复合结构处于弹性粘结阶段. 假设两种双轴拉伸加载情况. Case 1: ${\varepsilon _{{\rm{s}}x}} = {\varepsilon _{{\rm{s}}y}} = \xi$%; case 2: ${\varepsilon _{{\rm{s}}x}} = \left( {0.2 + \xi } \right)$%, ${\varepsilon _{{\rm{s}}y}} = \left( {\xi - {\upsilon _{\rm{s}}} \times 0.2} \right)$% (ξ可取任意大于${\upsilon _{\rm{s}}} \times 0.2$的值). 将case 1与case 2分别代入(17a)式、(17b)式和(18a)式、(18b)式进行计算, case 2– case 1即为基底受单轴拉伸载荷${\varepsilon _{{\rm{s}}x}} = 0.2$%的结果(本文取$\xi = 0.3$), 石墨烯正应力为:
$\left\{ \begin{aligned}& {\sigma _x}\left| {_{{\varepsilon _{{\rm{s}}x}} = 0.2\% }} \right. = {\sigma _x}\left| {_{\rm case2}} \right. - {\sigma _x}\left| {_{{\rm{case1}}},} \right.\\& {\sigma _y}\left| {_{{\varepsilon _{{\rm{s}}x}} = 0.2\% }} \right. = {\sigma _y}\left| {_{{\rm{case2}}}} \right. - {\sigma _y}\left| {_{{\rm{case1}}},} \right.\end{aligned} \right.$
将(19)式代入(20)式, 求得石墨烯正应变为:
$\left\{ \begin{aligned}& {\varepsilon _x}\left| {_{{\varepsilon _{{\rm{s}}x}} = 0.2\%}} \right. = \frac{1}{E}\left( {{\sigma _x} - \upsilon {\sigma _y}} \right),\\& {\varepsilon _y}\left| {_{{\varepsilon _{{\rm{s}}x}} = 0.2\%}} \right. = \frac{1}{E}\left( {{\sigma _y} - \upsilon {\sigma _x}} \right).\end{aligned} \right.$
界面切应力为:
$\left\{ \begin{aligned}& {\tau _{zx}}\left| {_{{\varepsilon _{{\rm{s}}x}} = 0.2\% }} \right. = {\tau _{zx}}\left| {_{{\rm{case}}2}} \right. - {\tau _{zx}}\left| {_{{\rm{case}}1}} \right.,\\& {\tau _{zy}}\left| {_{{\varepsilon _{{\rm{s}}x}} = 0.2\% }} \right. = {\tau _{zy}}\left| {_{{\rm{case}}2}} \right. - {\tau _{zy}}\left| {_{{\rm{case1}}}.} \right.\end{aligned} \right.$
计算过程详见文献[16], 由于篇幅原因此处不做记述.
将上文所述材料参数代入(20)式和(21)式, 石墨烯正应变εxεy分布如图5(a)图5(b)所示, 界面切应力τzxτzy分布如图5(c)图5(d)所示. 从图5(a)图5(b)中可以看出, εxεy均在石墨烯薄膜中心区域达到最大, 且εy最大值大约是εx最大值的33%. 此外, 由于石墨烯泊松比引起的xy方向的耦合效应, εx在靠近边缘($x = 0$$x = L$)区域内沿宽度(y方向)的分布与εy在靠近边缘($y = 0$$y = W$)区域内沿长度(x方向)的分布皆具有一定的不均匀性; 图5(c)图5(d)结果显示, τzxτzy在界面边角位置达到最大值, 且τzy最大值大约是τzx最大值的30%, 横向切应力τzy的传递作用影响较大. 另外, 与正应变分布相似, τzx沿宽度方向与τzy沿长度方向的分布同样具有一定的不均匀性. 综上结果可知, 二维模型计算的石墨烯正应变分布与界面切应力分布表现出明显的二维应力分布状态, 这是一维模型所不能解释的.
图 5 弹性粘结阶段 (εsx = 0.2%) 时石墨烯正应变 (a) εx和(b) εy以及界面切应力 (c) τzx和(d) τzy的分布
Figure5. 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), 有限元模型建立具体过程详见3.3节. 典型线分别选取沿石墨烯两个边缘($x = 0$$y = 0$)与两条中心线(x = $ {L / 2}$$y = {W / 2}$). 从图6中可以清楚地看出εxτzx沿宽度方向分布以及εyτzy沿长度方向分布的不均匀性, 而且理论模型结果与有限元结果符合较好.
图 6 典型线处石墨烯正应变 (a) εx和(b) εy以及界面切应力 (c) τzx和(d) τzy的分布
Figure6. Distributions of graphene’s strains (a) εx and (b) εy; distributions of interfacial shear stresses (c) τzx and (d) τzy along several representative lines.

2
3.2.滑移临界应变
-->将界面开始出现滑移区域时基底所受拉伸应变载荷定义为滑移临界应变εsxc. 根据图5(c)图5(d)可知, 弹性粘结阶段时最大切应力出现在界面的四个边角(图1中点O, A, BD), 所以界面边角位置最先发生滑移损伤. 当边角处切应力合力τm达到剪切强度τc时, 界面开始出现滑移区域. 由于结构对称性, 选取O点处界面切应力进行计算, 根据(9)式有
${\tau _{\rm{m}}}\left( {0,0} \right) = \sqrt {\tau _{zx}^{\rm{2}}\left( {{\rm{0,0}}} \right) + \tau _{zy}^2\left( {0,0} \right)} = {\tau _{\rm{c}}}.$
将(16a)式和(16b)式代入(4)式中可得:
$\left\{ \begin{aligned}& {\tau _{zx}}\left( {0,0} \right) = \alpha {\varepsilon _{{\rm{s}}x}}Eh\frac{{\left( {1 - \cosh \alpha L} \right)}}{{\sinh \alpha L}},\\& {\tau _{zy}}\left( {0,0} \right) = \alpha {\varepsilon _{{\rm{s}}y}}Eh\frac{{\left( {1 - \cosh \alpha W} \right)}}{{\sinh \alpha W}},\end{aligned} \right.$
式中${\varepsilon _{{\rm{s}}y}} = - {\upsilon _{\rm{s}}} \times {\varepsilon _{{\rm{s}}x}}$. 由(22)式和(23)式计算得到滑移临界应变εsxc:
$\begin{split}& {\varepsilon _{{\rm{s}}x{\rm{c}}}} = \frac{1}{{\sqrt {KEh} }} \times \\& \frac{{{\tau _{\rm{c}}}}}{{\sqrt {{{\left( {\dfrac{{\cosh (\alpha L) - 1}}{{\sinh (\alpha L)}}} \right)}^2} + \upsilon _{\rm{s}}^2{{\left( {\dfrac{{\cosh (\alpha W) - 1}}{{\sinh (\alpha W)}}} \right)}^2}} }}.\end{split}$
Cui和Guo[13]利用一维模型研究了石墨烯长度尺寸对滑移临界应变的影响. 而依据二维模型推导 (24)式可知, 石墨烯宽度尺寸以及基底泊松比大小同样会对滑移临界应变造成影响. 将3.1节所述材料参数代入(24)式, 图7给出了不同基底泊松比情况下滑移临界应变εsxc随石墨烯宽度的变化关系. 从图7可以看出, 滑移临界应变存在一个特征宽度, 当石墨烯宽度小于该特征量时, 滑移临界应变εsxc表现出明显的宽度尺寸效应, 石墨烯宽度越大, εsxc越小, 界面边缘越易发生滑移损伤; 当石墨烯宽度大于该特征量时, 宽度尺寸对εsxc的影响可被忽略. 对于本文选取的石墨烯与界面力学参数, εsxc的特征宽度大约为8 μm. 而通过分析 (24)式中参数量α可知, α为界面刚度K与石墨烯面内刚度Eh比值的二分之一次方, 1/α在微米量级, 所以滑移临界应变的特征宽度出现在微米量级是由界面力学性能与石墨烯自身力学性质所共同决定的. 另外, 对比图7中三条曲线结果可知, 石墨烯宽度相同时, 基底泊松比越大, εsxc 越小. 此外, 图7还给出了一维模型与二维模型对滑移临界应变计算结果的比较. 结果显示, 一维模型过高估计了滑移临界应变εsxc. 当石墨烯宽度大于8 μm, 基底泊松比为0.49时, 一维模型与二维模型计算结果的偏差可达11.3%, 但当石墨烯宽度小于0.5 μm时, 两种模型对εsxc计算结果的差别低于0.5%, 可被忽略.
图 7 不同基底泊松比情况下滑移临界应变εsxc随石墨烯宽度W的变化(线为理论值, 散点为有限元值)
Figure7. 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).

2
3.3.界面滑移阶段
-->当基底所受拉伸载荷较大时, 界面开始出现滑移区域, 石墨烯/基底结构进入界面滑移阶段, 如图8所示. 在ABAQUS中建立石墨烯/基底结构三维有限元模型, 计算界面滑移阶段时石墨烯正应变以及界面切应力的分布. 分别采用M3D4R单元与C3D8R单元对石墨烯和基底进行建模, 材料参数按照3.1节所述进行设置, 基底 (PET) 杨氏模量${E_{\rm{s}}} = 3\;{\rm{GPa}}$. 根据2.2节所述方法, 采用具有双线性内聚力法则的三维COH3 D8单元代替二维模型对界面进行建模. 在切向两方向上, 分别设置界面刚度$K = 74\;{{{\rm{TPa}}} / {\rm{m}}}$, 界面剪切强度${\tau _{\rm{c}}} = 0.5\;{\rm{MPa}}$[7]. 在垂直方向上, 界面刚度和界面强度被视为惩罚参数, 设置为较大的数值, 防止垂直界面方向相对位移的产生. 当界面切应力合力达到剪切强度τc时, 界面进入软化阶段, 采用B-K准则判断断裂起始, 将界面切向断裂能设置为一个相当大的数值, 此时双线性内聚力模型软化阶段可近似替代二维模型界面滑移阶段.
图 8 界面滑移阶段示意图
Figure8. Schematic diagram of interfacial sliding stage.

基底所受拉伸应变${\varepsilon _{{\rm{s}}x}} = $1%时, 复合结构进入界面滑移阶段. 石墨烯正应变εxεy分布数值结果如图9(a)图9(b)所示, 界面切应力τzxτzy分布数值结果如图9(c)图9(d)所示. 图9结果显示, 界面滑移阶段时, 石墨烯正应变与界面切应力的分布同样表现出明显的二维应力分布状态. εxεy仍在石墨烯中心区域达到最大, 而界面切应力则由于界面边缘发生滑移损伤, 最大值位置发生变化. τzx在界面中心线 ($y = {W / 2}$) 上的滑移区域内达到最大, 而τzy则在边缘线($y = 0$$y = W$)上滑移区域与弹性粘结区域交界点处达到最大值, 且τzy最大值大约是τzx最大值的60%, 横向切应力τzy的传递作用进一步增大. 此外, 界面滑移区域的出现也使得εx, τzx沿宽度方向分布与εy, τzy沿长度方向分布的不均匀性进一步增大.
图 9 界面滑移阶段 (εsx = 1%) 时石墨烯正应变 (a) εx和(b) εy以及界面切应力 (c) τzx和(d) τzy.的分布
Figure9. 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%).

2
4.1.与一维模型计算结果的对比
-->根据3.1节相关参数的取值, 将二维模型的计算结果与一维模型结果进行细致对比. 首先, 一维模型只考虑了沿拉伸方向的界面切应力传递, 忽略了界面横向切应力传递作用, 无法考虑石墨烯宽度以及基底泊松比大小对滑移临界应变εsxc的影响, 也无法对石墨烯压缩应变εy以及界面横向切应力τzy的分布进行计算. 3.2节中已经将两模型对于滑移临界应变εsxc的计算结果进行了比较(如图7所示), 并发现一维模型过高估计了复合结构的滑移临界应变, 但当石墨烯宽度小于0.5 μm时, 两种模型对εsxc计算结果的差别可以忽略. 而对于石墨烯压缩应变εy, 由于没有考虑界面横向切应力传递作用, 可认为一维模型对其的计算结果为
${\varepsilon _y} = - \upsilon \times {\varepsilon _x}.$
图10给出了基底受不同拉伸载荷作用下, 二维模型对石墨烯中心C点处(如图1所示)压缩应变εyC的计算结果随其宽度的变化关系. 结果显示, 无论复合结构处于弹性粘结阶段 (${\varepsilon _{{\rm{s}}x}} = 0.2\%$) 还是界面滑移阶段 (${\varepsilon _{{\rm{s}}x}} = $0.6%, 1%), 石墨烯宽度大于20 μm时, 由于界面横向切应力传递区域足够大, 基底横向压缩应变可以完全传递到石墨烯上, 此时$\varepsilon _y^C \approx - {\upsilon _s} \times {\varepsilon _{sx}}$; 当宽度小于20 μm, 大于0.5 μm时, 虽然基底横向压缩应变并未完全传递到石墨烯上, 但界面横向切应力传递对石墨烯中心C点处压缩应变εyC的增大作用已经不可忽略, 并且随宽度增加, 该作用进一步增大; 当宽度小于0.5 μm时, 则有$\varepsilon _y^{C} \approx - {\upsilon _{}} \times \varepsilon _x^{C}$, 与(25)式计算结果相同, 此时由于传递区域过小, 界面横向切应力τzy的传递作用可被忽略.
图 10 不同基底应变下石墨烯中心C点处压缩应变εyC随石墨烯宽度的变化
Figure10. Variation of compressive strain εyC at the center point C of graphene with its width when the strain of substrate is different.

选取两处典型位置, 即沿边缘 ($y = 0$) 和中心线 ($y = {W / 2}$) 处, 将二维模型对石墨烯拉伸应变εx和界面切应力τzx的计算结果与一维模型的结果分别进行了比较(如图11所示). 图11(a)图11(b)给出了弹性粘结阶段 (${\varepsilon _{{\rm{s}}x}} = $0.2%)的比较结果. 由图11结果显示, 二维模型对εxτzx沿边缘线分布的计算结果与一维模型重合, 沿中心线分布的计算结果与一维模型具有一定差别. 图11(c)图11(d)给出了界面滑移阶段 (${\varepsilon _{{\rm{s}}x}} = $1%)的比较结果. 结果显示, 界面滑移阶段时, 二维模型对εxτzx沿两条典型线上分布的计算结果与一维模型结果均具有较大差别. 特别对于τzx, 由于二维模型考虑了界面横向切应力τzy, 界面滑移行为由τzxτzy共同决定. 在界面滑移阶段, 沿界面边缘线 (y = 0) 上的滑移区域内存在较大的横向切应力τzy, 而在沿界面中心线 (y = W/2) 上横向切应力τzy为零(如图9(d)), 一维模型则没有考虑横向切应力τzy. 由此导致二维模型对τzx沿界面边缘分布的计算结果与对沿中心线分布的计算结果以及一维模型结果均具有较大差异, 特别对于沿边缘线上x小于5.6 μm的界面滑移部分, 差异更加明显.
图 11 二维模型与一维模型结果的比较 (W = 21.8 μm) (a) εx和(b) τzx在弹性粘结阶段(εsx = 0.2%); (c) εx和(d) τzx在界面滑移阶段 (εsx = 1%)
Figure11. 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%).

但经计算发现, 当石墨烯宽度足够小时, 二维模型对εxτzx沿典型线分布的计算结果与一维模型结果的差别可被忽略, 且εxτzx沿宽度方向的分布基本均匀. 图12(a)图12(b)给出了石墨烯宽度为1 μm时界面滑移阶段的比较结果, 两模型计算结果基本重合. 弹性粘结阶段时比较结果与之相似, 此处不做赘述.
图 12 二维模型与一维模型计算结果的比较 (W = 1 μm) (a) εx和(b) τzx在界面滑移阶段 (εsx = 1%)
Figure12. 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%).

综上结果可知, 当石墨烯宽度较大时, 一维模型会高估复合结构的滑移临界应变, 对石墨烯拉伸应变εx与界面切应力τzx的计算结果与二维模型有一定误差, 尤其在界面滑移阶段时, 差别较大. 且通过二维模型结果可知, 石墨烯薄膜存在较大压缩应变εy, 界面横向切应力τzy的传递作用不可忽略. 但当石墨烯宽度小于0.5 μm时, 针对本文所代入的界面力学参数, 一维模型与二维模型对滑移临界应变计算结果的差别低于0.5%, 对εxτzx的计算结果基本重合, 且二维模型对εy的计算结果约等于石墨烯自身泊松比所引起的横向压缩应变, 与一维模型根据(25)式对其的计算结果相同, 此时界面横向切应力τzy的传递作用可被忽略, 二维模型可近似被一维模型代替.
2
4.2.与实验结果的对比
-->许超宸[23]利用原位拉曼光谱技术测量了放置在受单轴拉伸PET基底上的长方形石墨烯沿拉伸方向中心线上应变εx的分布. 实验结果显示, 小尺寸石墨烯与基底完美贴合, 符合本文所建立模型的界面完美粘结假设. 本节利用二维模型与许超宸[23]测得的小尺寸石墨烯 ($L = 20 $μm, $W = 3\;{\rm{mm}}$)的实验数据进行拟合对比, 确定界面参数, 验证模型可靠性.
与实验数据进行拟合, 就需对拉曼技术测量石墨烯应变的实验原理与细节进行了解. 许超宸[23]通过检测拉曼谱线中2D峰峰位ω2D对石墨烯应变分布进行测量. 2D峰频移变化Δω2D与石墨烯面内正应变分量之和${\varepsilon _x} + {\varepsilon _y}$有如下对应关系[24]:
$\Delta {\omega _{{\rm{2D}}}} = - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}}\left( {{\varepsilon _x} + {\varepsilon _y}} \right),$
其中$\Delta {\omega _{2{\rm{D}}}} = {\omega _{2{\rm{D}}}} - \omega _{2{\rm{D}}}^0$, $ \omega _{{\rm{2D}}}^{\rm{0}} $为被测位置2D峰初始峰位, γ2D为Grüneisen系数, 定义石墨烯面内正应变分量之和${\varepsilon _{\rm{m}}} = {\varepsilon _x} + {\varepsilon _y}$. 依据(26)式测量石墨烯应变, 首先要对$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项进行标定. Mohiuddin等[24]的标定实验中考虑了泊松效应的影响, 选择石墨烯中心区域为观测点, 当基底所受拉伸载荷较小时, 2D峰峰位ω2D随基底拉伸应变εsx的增大线性变化, 此时认为石墨烯中心区域应变与基底应变协同一致, 即有
${\varepsilon _x} = {\varepsilon _{{\rm{s}}x}},\; {\varepsilon _y} = - {\upsilon _{\rm{s}}}{\varepsilon _{{\rm{s}}x}}.$
将(27)式代入(26)式得到
$- \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} = \frac{{\Delta {\omega _{{\rm{2D}}}}}}{{\left( {1 - {\upsilon _{\rm{s}}}} \right){\varepsilon _{{\rm{s}}x}}}}.$
利用2D峰峰位ω2D随基底拉伸应变εsx变化的线性段的斜率结合式(28)式, 完成对$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项的标定. 但在许超宸[23]的实验中, 采用了一维模型假设, 忽略石墨烯与基底的泊松比(假设${\upsilon _{\rm{s}}} = \upsilon = 0$). 利用2D峰线性变化段对$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项进行标定时认为
$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} \approx \frac{{\Delta {\omega _{{\rm{2D}}}}}}{{{\varepsilon _{{\rm{s}}x}}}}.$
在完成标定后, 利用标定好的$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项建立2D峰峰位ω2D与石墨烯应变的对应关系时, 假设${\varepsilon _y} = 0$, 并根据(26)式得到如下应变测定公式:
${\varepsilon _x} \approx - \frac{{\Delta {\omega _{{\rm{2D}}}}}}{{\omega _{{\rm{2D}}}^0{\gamma _{{\rm{2D}}}}}}.$
但由(26)式可知, 实际与Δω2D成正比的是石墨烯面内正应变分量之和εm. 并且根据 (28)式与(29)式对比可知, 许超宸忽略石墨烯与基底泊松比作用导致实验中标定的$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项实际为$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $×(1 – υs), 所以其依据(30)式测出的应变εx与实际应变值可能存在一定的误差. 本文建立的二维模型考虑了石墨烯与基底的泊松比作用, 并解耦计算得到石墨烯拉伸应变分量εx和压缩应变分量εy. 所以与许超宸[23]实验数据拟合时, 需要考虑压缩应变分量εy, 将许超宸依据 (30)式所测的实验数据乘以$1 - {\upsilon _{\rm{s}}}$项, 将其还原回实验精确测量值石墨烯面内正应变之和εm. 此外在标定实验中, 许超宸只在石墨烯中心区域对$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项进行了标定, 并将之用于沿石墨烯整个长度上应变分布的测量中, 而$ - \omega _{{\rm{2D}}}^{\rm{0}}{\gamma _{{\rm{2D}}}} $项较为敏感[25-27], 可能会随标定位置的不同而改变, 所以其实验中相对准确的测量区域为石墨烯中心位置.
根据3.3节所述方法, 在ABAQUS中建立石墨烯/基底三维有限元模型. 图13(a)给出了基底拉伸载荷${\varepsilon _{{\rm{s}}x}} =$0.25%时, 石墨烯面内正应变分量之和εm沿中心线 ($y = {W / 2}$) 分布的模拟值与实验测量值的拟合结果. 图中结果显示, 模拟值与实验值拟合良好, 从而确定界面刚度K为100 TPa/m. 图13(b)则给出了不同基底载荷作用下石墨烯中心C点处$ \varepsilon_{\rm m}^{\rm C} $的模拟值与实验值的拟合结果, 由此确定界面剪切强度τc为0.295 MPa. 而许超宸[23]通过一维模型拟合实验数据得到同尺寸石墨烯与PET界面剪切强度为0.314 MPa, Guo和Zhu[12]通过一维模型拟合得到石墨烯与PET界面刚度为74 TPa/m.
图 13 利用二维模型与实验数据拟合 (a) 基底拉伸应变εsx = 0.25%时εm沿石墨烯中心线(y = W/2)的分布; (b) 不同基底载荷作用下石墨烯中心C点处$ \varepsilon_{\rm m}^{C} $的大小
Figure13. 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.

由于二维模型只考虑了石墨烯/基底沿界面切向的应力传递, 未考虑石墨烯在较大压缩应变下的面外屈曲行为. 实验研究显示[28], 放置在基底上的单层石墨烯, 当基底所受压缩应变载荷大于0.3%时, 石墨烯发生面外屈曲, 所以本文只在基底受较小载荷情况下, 用模拟值与实验值进行了拟合.
考虑了石墨烯和基底泊松效应的影响, 本文提出了二维非线性剪滞模型. 对于基底泊松比相比石墨烯较大的情况, 利用该模型理论研究了受单轴拉伸石墨烯/柔性基底复合结构的双向界面切应力传递问题. 在弹性粘结阶段, 导出了石墨烯双向正应变和双向界面切应力的半解析表达式, 分析了不同位置处石墨烯正应变和界面切应力的分布规律. 导出了石墨烯/基底复合结构发生界面滑移的临界应变, 分析了石墨烯宽度以及基底泊松比大小对滑移临界应变的影响. 基于二维非线性剪滞模型建立有限元模型, 研究了界面滑移阶段石墨烯正应变和界面切应力的分布规律. 将二维非线性剪滞模型计算结果与一维非线性剪滞模型的结果进行了比较. 研究结果表明: 当石墨烯宽度较大时, 无论在弹性粘结阶段还是界面滑移阶段, 石墨烯正应变和界面切应力都呈现出二维应力分布状态, 石墨烯内存在较大的压缩应变分量εy, 界面上存在较大的横向切应力分量τzy, 这在一维模型中是被忽略的. 而且εxτzx沿宽度方向的分布呈现出不均匀性, 与一维模型对其的计算结果具有显著差别. 滑移区域首先出现在界面的边角位置, 并且滑移临界应变εsxc存在一个微米量级的特征宽度, 当石墨烯宽度小于该特征量时, 宽度越大, 临界应变εsxc越小, 界面边缘更易发生滑移损伤. 但当石墨烯宽度大于该特征量时, 宽度尺寸变化对εsxc的影响可被忽略. 另外, 基底泊松比同样会对滑移临界应变εsxc造成影响, 基底泊松比越大, εsxc越小. 一维模型则会高估复合结构的滑移临界应变. 当石墨烯宽度很小时, 二维模型可以退化为一维模型, 界面横向切应力传递作用可被忽略.
最后, 通过与拉曼实验结果的对比, 验证了二维非线性剪滞模型的可靠性, 并拟合得到了石墨烯/PET基底结构的界面刚度 (100 TPa/m) 和界面剪切强度 (0.295 MPa), 而根据一维模型拟合实验数据所得界面剪切强度值为0.314 MPa.
3.1节中(17a)式与(18a)式中参数具体表达式如下:
$\small{f_{mn,x}} = \frac{{ - 16F}}{{mn{{\text{π}}^2}\left[ {{{\left( {\dfrac{{m{\text{π}}}}{L}} \right)}^2}{{\left( {\dfrac{{n{\text{π}}}}{{2{\varPhi _{m,x}} + W}}} \right)}^2} + A{{\left( {\dfrac{{n{\text{π}}}}{{2{\varPhi _{m,x}} + W}}} \right)}^2} + A{{\left( {\dfrac{{m{\text{π}}}}{L}} \right)}^2} + E} \right]}} ,\tag{A1}$
$\small{\varPhi _{m,x}} = - \frac{1}{{\sqrt d }}\ln \left( {1 - \frac{{Cd}}{F}} \right),\tag{A2}$
其中$E = {A^2} - {B^2}$, $F = A \times C - B \times D$, $ d = \dfrac{{A{{\left( {{{m{\text{π}}}}/{L}} \right)}^2} + E}}{{{{\left( {{{m{\text{π}}}}/{L}} \right)}^2} + A}}. $
(17b)式与(18b)式中参数具体表达式如下:
${f_{mn,y}} = \dfrac{{ - 16G}}{{mn{{\text{π}}^2}\left[ {{{\left( {\dfrac{{m{\text{π}}}}{W}} \right)}^2}{{\left( {\dfrac{{n{\text{π}}}}{{2{\varPhi _{m,y}} + L}}} \right)}^2} + A{{\left( {\dfrac{{n{\text{π}}}}{{2{\varPhi _{m,y}} + L}}} \right)}^2} + A{{\left( {\dfrac{{m{\text{π}}}}{W}} \right)}^2} + E} \right]}},\tag{A3}$
${\varPhi _{m,y}} = - \frac{1}{{\sqrt b }}\ln \left( {1 - \frac{{Db}}{G}} \right),\tag{A4}$
其中$G = A \times D - B \times C$, $ d = \dfrac{{A{{\left( {{{m{\text{π}}}}/{L}} \right)}^2} + E}}{{{{\left( {{{m{\text{π}}}}/{L}} \right)}^2} + A}}. $
相关话题/计算 结构 实验 石墨 基底

  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效,10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材,产品类型包含电子书、题库、全套资料以及视频,无论您是考研复习、考证刷题,还是考前冲刺等,不同类型的产品可满足您学习上的不同需求。 ...
    本站小编 Free壹佰分学习网 2022-09-19
  • H, F修饰单层GeTe的电子结构与光催化性质
    摘要:采用基于密度泛函理论的第一性原理平面波赝势方法,计算了单层GeTe、表面氢化及氟化单层GeTe的晶体结构、稳定性、电子结构和光学性质.计算结果表明,经过修饰后,GeTe的晶格常数、键角、键长增大,且均具有较好的稳定性.电子结构分析表明,单层GeTe为间接带隙半导体,全氢化修饰、全氟化修饰以及氢 ...
    本站小编 Free考研考试 2021-12-29
  • 考虑界面接触热阻的一维复合结构的热整流机理
    摘要:建立了考虑变截面、变热导率及界面接触热阻效应的组合热整流结构的温度场及热整流系数的理论模型和有限元解.数值算例证明了本文模型及算法的可靠性,进而通过参数影响研究确定了若干几何及材料参数对结构热整流系数的影响规律,揭示界面接触热阻对热整流效果的影响机理.研究结果表明长度比、截面半径变化率、热导率 ...
    本站小编 Free考研考试 2021-12-29
  • 肖特基钙钛矿太阳电池结构设计与优化
    摘要:有机-无机杂化钙钛矿材料有高吸收系数、低廉的制作成本以及较为简单的制备工艺,在近年来表现出良好的发展前景.本文采用wx-AMPS模拟软件对平面结构钙钛矿太阳电池和肖特基钙钛矿太阳电池进行建模仿真对比,从理论上分析无载流子传输层的肖特基钙钛矿太阳电池的优势.结果显示,器件两侧电极功函数和吸收层的 ...
    本站小编 Free考研考试 2021-12-29
  • 三维金字塔立体复合基底表面增强拉曼散射特性
    摘要:表面增强拉曼散射(SERS)是一种广泛应用于化学反应检测、医学诊断和食品分析等领域的高灵敏度光谱技术.基底结构的构建对提高探针分子的SERS信号有非常重要的影响.本文利用聚甲基丙烯酸甲酯(PMMA)包裹银纳米颗粒制备了一种三维金字塔立体复合SERS基底,实现了对罗丹明6G(R6G)分子的高灵敏 ...
    本站小编 Free考研考试 2021-12-29
  • Verwey相变处Fe<sub>3</sub>O<sub>4</sub>的结构、磁性和电输运特性
    摘要:作为典型的金属–绝缘体转变,Fe3O4的Verwey相变蕴涵的丰富物理现象与微观机制,因而受到了人们的广泛关注.在Verwey相变处,Fe3O4的晶体结构、电子结构以及磁各向异性等均发生转变,但其磁基态并未发生改变.与其他强关联体系相比,Fe3O4的Verwey相变不需要考虑磁交换耦合作用的变 ...
    本站小编 Free考研考试 2021-12-29
  • 高效外腔倍频产生426 nm激光的实验研究
    摘要:利用与铯原子吸收线对应的852nm半导体激光作为基频光,泵浦基于周期极化磷酸钛氧钾(PPKTP)晶体的环形腔,进行高效外腔谐振倍频并产生426nm激光.在理论分析小角度环形腔内的热透镜效应基础上,发现晶体中等效热透镜中心位置并非在晶体的几何中心.在理论分析的基础上,实验上通过精密平移台精细调节 ...
    本站小编 Free考研考试 2021-12-29
  • 高熵合金短程有序现象的预测及其对结构的电子、磁性、力学性质的影响
    摘要:如何有效预测高熵合金的稳态结构,是开展研究其物理及化学等性能的基础.以FeCuCrMnMo合金为例,在有限晶胞尺寸内,采用蒙特卡洛结合密度泛函理论杂化计算方法(MonteCarlo/densityfunctionaltheory,MC/DFT)预测高熵合金的平衡态结构.与准随机近似方法(spe ...
    本站小编 Free考研考试 2021-12-29
  • 石墨烯纳米网电导特性的能带机理:第一原理计算
    摘要:通过第一原理电子结构计算来研究有序多孔纳米网的电导特性变化的能带机理.能带结构分析结果表明:石墨烯纳米网超晶格(3m,3n)(m和n为整数)的电子本征态在布里渊区中心点发生四重简并;碳空位孔洞规则排列形成的石墨烯纳米网具有由简并态分裂形成的宽度可调带隙,无论石墨烯的两个子晶格是否对等.在具有磁 ...
    本站小编 Free考研考试 2021-12-29
  • 腔光子-自旋波量子耦合系统中各向异性奇异点的实验研究
    摘要:通过耦合三维微波腔中光子和腔内钇铁石榴石单晶小球中的自旋波量子形成腔-自旋波量子的耦合系统,并通过精确调节系统参数在该实验系统中观测到各向异性奇异点.奇异点对应于非厄米系统中一种特殊状态,在奇异点处,耦合系统的本征值和本征矢均简并,并且往往伴随着非平庸的物理性质.以往大量研究主要集中在各向同性 ...
    本站小编 Free考研考试 2021-12-29
  • Tl<sub>0.33</sub>WO<sub>3</sub>电子结构和太阳辐射屏蔽性能第一性原理研究
    摘要:节能减排已成为当今社会发展的主题,对节约能源、提高太阳能的高效综合利用的新型窗用透明隔热材料的理论设计和研究尤其重要.本文采用基于密度泛函理论的计算方法,研究了六方相三氧化钨Tl掺杂前、后的晶格参数、电子能带结构、形成能和光学性质.研究结果表明,Tl掺杂后晶格体积增大,系统能量降为负值,体系具 ...
    本站小编 Free考研考试 2021-12-29