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

单层MoS<sub>2</sub>的热弹耦合非线性板模型

本站小编 Free考研考试/2021-12-29
摘要:单层二硫化钼(${\rm{Mo}}{{\rm{S}}_2}$)是有广泛应用前景的二维纳米材料, 但其力学性质还没有被深入研究, 特别是其热弹耦合力学行为迄今还没有被关注到. 本文首次提出了考虑热应力影响的单层${\rm{Mo}}{{\rm{S}}_2}$的非线性板理论, 并对比研究了其与石墨烯的热弹耦合力学性质. 对于不可移动边界, 结果显示: 1)有限温度产生的热应力降低了${\rm{Mo}}{{\rm{S}}_2}$的刚度, 但提高了石墨烯的刚度; 2)在相同几何尺寸和温度条件下, 变形较小时${\rm{Mo}}{{\rm{S}}_2}$的刚度大于石墨烯, 但伴随变形的增大, ${\rm{Mo}}{{\rm{S}}_2}$的刚度将小于石墨烯. 研究结果表明, 边界预加轴向外力和环境温度可以调节单层二维纳米结构力学性质. 本文建立的热弹耦合板模型, 可以推广至其他单层二维纳米结构.
关键词: 单层MoS2/
石墨烯/
热应力/
非线性板理论

English Abstract


全文HTML

--> --> -->
1.引　言
自机械剥离获得单层石墨烯以来, 多种单层二维纳米结构不断被发现[1-4]. 比如, 与石墨烯一样具有单原子厚度的六方氮化硼(h-BN)[5]; 原子不在同一面内的准二维纳米单层, 例如硅烯[6]; 以及不同原子构成的“三片单层结构”(three-sheet single layer), 例如二硫化钼(${\rm{Mo}}{{\rm{S}}_2}$)[7,8]等. 最近还发现了“五片单层结构”的${\rm{MoS}}{{\rm{i}}_2}{{\rm{N}}_4}$[9]. 石墨烯因优异的物理化学性质, 具有广泛的潜在应用价值. 其力学性质引起了广泛的兴趣, 并得了大量的研究成果[10-12]. 但针对${\rm{Mo}}{{\rm{S}}_2}$单层的研究还主要集中在电学、热学和摩擦性质上, 对其力学性质的研究还较少[13]. 特别是, ${\rm{Mo}}{{\rm{S}}_2}$面外抗弯能力的研究还很少被关注[13-15]. 此外对于二维纳米结构, 实际应用中常常需要考虑环境温度变化的影响. 虽然${\rm{Mo}}{{\rm{S}}_2}$的热学性质已经引起了研究者的关注[16-18]. 但迄今为止, 还没有针对单层${\rm{Mo}}{{\rm{S}}_2}$热弹耦合(热-力耦合)的系统研究. 任何材料在工程应用前, 都需要对其力学性质有基本的认识, ${\rm{Mo}}{{\rm{S}}_2}$也不会例外. 纳米结构通常对温度变化比较敏感, 因此在力学模型中考虑温度的影响是必要的.
需要指出的是, 迄今还没有统一的理论来描述单层二维结构的力学性质[19-20]. 这是因为, 如何建立化学键理论和二维材料力学性质的关系, 是具有高度挑战性的难题. 考虑到以量子化学为基础, 来建立“多片”单层二维结构力学理论的巨大困难, 研究者通常退而求其次, 使用唯像的宏观连续介质理论来描述其力学性质[19]. 连续介质理论中的力学参数, 则使用量子计算和实验来拟合. 其中密度泛函紧束缚(FDTB)计算和分子动力学(MD)是两种被广泛使用的方法[19-24]. 迄今的文献中, 多数${\rm{Mo}}{{\rm{S}}_2}$的力学参数是在绝对零度下得到的, 忽略了有限温度的影响[13-15]. 但温度对二维材料力学性质的影响, 在理论和实践中都是需要考虑的. 对于石墨烯, 温度对弹性模量, Poisson比、以及面外弯曲和扭转刚度的影响不大[21]. 对于${\rm{Mo}}{{\rm{S}}_2}$, 弹性参数和温度的关系还缺乏系统的研究. 近期的分子动力学计算(MD)显示, ${\rm{Mo}}{{\rm{S}}_2}$在温度较低情况下可以忽略温度对弹性参数的影响[18]. 虽然温度对${\rm{Mo}}{{\rm{S}}_2}$的弹性参数影响较小, 但温度会产生显著的热膨胀效应. 一个有趣的现象是, 在有限温度下(1000 K以下), ${\rm{Mo}}{{\rm{S}}_2}$热胀冷缩, 而同为单层的石墨烯却是热缩冷胀[16,17]. 造成这种截然相反结果的原因, 是石墨烯和${\rm{Mo}}{{\rm{S}}_2}$有着不同的热变形机制[22]: 由于面外弯曲刚度很小, 石墨烯在温度场中波浪式起伏, 并造成负的热膨胀系数; 但${\rm{Mo}}{{\rm{S}}_2}$的面外弯曲刚度远远大于石墨烯, 可有效抑制温度造成的起伏; 这使得${\rm{Mo}}{{\rm{S}}_2}$在高温下的体积膨胀, 导致正的热膨胀系数. 但无论是正或负的热膨胀系数, 在边界受限时都将产生热应力. 二维材料的抗弯能力通常较弱, 热应力可能导致热屈曲, 并使结构丧失功能. 因此澄清二维材料的热弹耦合力学性质, 在工程应用中具有重要的意义. 本文将基于唯像的宏观F?ppl-von Karman 板理论, 来建立具有独立面内和面外力学参数的单层${\rm{Mo}}{{\rm{S}}_2}$的热弹耦合力学理论. 该理论可方便的推广到其他单层二维纳米材料.
2.单层${\rm{Mo}}{{\rm{S}}_2}$的热弹耦合非线性板理论
对于二维材料, 被广泛采用的宏观连续介质力学模型是F?ppl-von Karman板理论[19,22-24].本文将建立 图1所示单层${\rm{Mo}}{{\rm{S}}_2}$的热弹耦合非线性板理论. 虽然还没有通过量子力学(或化学键)理论证明, 三原子层的二维结构的变形能(自由能)可写为经典的板理论形式. 但该唯像理论可以得到和量子计算一致的结构变形 [14-15]. 经典的F?ppl-von Karman板壳理论中, 绝度零度时的变形能包含面内的张拉和剪切变形能, 以及面外的弯曲和扭转变形能[25], 如方程(1)所示. 出人意料的是, 该变形能很好地描述了单层石墨烯的力学性质[25-27]:
图 1 单层${\rm{Mo}}{{\rm{S}}_2}$计算简图:　(a) 顶视图; (b)侧视图; (c)等效板立体图; (d)边界载荷
Figure1. Computational model of single-layer ${\rm{Mo}}{{\rm{S}}_2}$: (a) Top view of the structure; (b) Side view of the structure; (c) Stereo plate model of the structure; (d) Applied edge loads

$\begin{split}\;&U = {U_{\rm{e}}} + {U_{\rm{s}}} =\\ & \int\nolimits_S \left[ \frac{1}{2}{k_{\rm{B}}}{{\left( {2H} \right)}^2} + {k_{\rm{G}}}K + \frac{1}{2}{k_{\rm{b}}}{{\left( {2J} \right)}^2} - {k_{\rm{g}}}Q \right]{\rm{d}}S .\end{split}$
方程(1)中, $H$, $K$为平均曲率和Gauss曲率, kB和kG分别为弯曲刚度和扭转刚度; $Q = {\rm{tr}}\left( {{\varepsilon _0}} \right)$, $J = \det \left( {{\varepsilon _0}} \right)$是板中面的二维应变张量${\varepsilon _0}$的两个不变量, ${k_{\rm{b}}}$和${k_{\rm{g}}}$是面内刚度. 在宏观F?ppl-von Karman板理论中, 面内和面外刚度间有如下关系: ${k_{\rm{B}}}/{k_{\rm{b}}} = - {k_{\rm{G}}}/{k_{\rm{g}}} = {h^2}/12$, $h$为板厚. 和经典的板理论不同, 二维材料的面内和面外刚度间没有此关系[19]. 即不能通过面内力学参数, 直接通过二维材料的厚度得到面外的弯曲和扭转刚度. 这种不一致性称为Yakobson悖论(Yakobson paradox)[19]. 迄今为止, 对Yakobson悖论还没有得到一致认可的物理解释. 文献[23]指出, 石墨烯的在有限变形时的各向异性导致了Yakobson悖论. 使用板的等效厚度概念, 可以方便地使用连续介质力学理论, 为设计提供方便[24]. 本文将放弃引入等效厚度, 直接把面内刚度和面外刚度作为独立的力学参数. 这样处理既回避了Yakobson悖论产生的不确定性, 在理论上也没有逻辑矛盾. 对比方程(1)和经典板壳的变形能, 有:
${k_{\rm{b}}} = \frac{Y}{{1 - {\nu ^2}}},\quad {k_{\rm{g}}} = \frac{Y}{{1 + \nu }}.$
其中, $Y$, $\nu $分别为二维弹性模量和Poisson比. kB和kG 需要通过原子计算或实验确定. 在直角坐标系中, 令$u, v, w$分别为$x, y, z$方向的位移, 如图1所示. 使用von Karman非线性应变, 则面内应变张量的分量为[25,27]
$\begin{split} &\varepsilon _{xx}^0 = \frac{{\partial u}}{{\partial x}} + \frac{1}{2}{\left( {\frac{{\partial w}}{{\partial x}}} \right)^2}, \\ & \varepsilon _{yy}^0 = \frac{{\partial v}}{{\partial y}} + \frac{1}{2}{\left( {\frac{{\partial w}}{{\partial y}}} \right)^2}, \\ &\varepsilon _{xy}^0 = \frac{1}{2}\left( {\frac{{\partial u}}{{\partial y}} + \frac{{\partial v}}{{\partial x}} + \frac{{\partial w}}{{\partial x}}\frac{{\partial w}}{{\partial y}}} \right). \end{split} $
面内应变张量的两个不变量为[25]
$\begin{split} & 2J = \frac{{\partial u}}{{\partial x}} + \frac{{\partial v}}{{\partial y}} + \frac{1}{2}\left[ {{{\left( {\frac{{\partial w}}{{\partial x}}} \right)}^2} + {{\left( {\frac{{\partial w}}{{\partial y}}} \right)}^2}} \right], \\& Q = \left[ {\frac{{\partial u}}{{\partial x}} + \frac{1}{2}{{\left( {\frac{{\partial w}}{{\partial x}}} \right)}^2}} \right]\left[ {\frac{{\partial v}}{{\partial y}} + \frac{1}{2}{{\left( {\frac{{\partial w}}{{\partial y}}} \right)}^2}} \right]\\&\;\;\;\;\;\;- \frac{1}{4}{\left( {\frac{{\partial u}}{{\partial y}} + \frac{{\partial v}}{{\partial x}} + \frac{{\partial w}}{{\partial x}}\frac{{\partial w}}{{\partial y}}} \right)^2}. \end{split} $
平均曲率和Gauss曲率为[28]
$\begin{split}&H \approx \frac{1}{2}\left( {\frac{{{\partial ^2}w}}{{\partial {x^2}}} + \frac{{{\partial ^2}w}}{{\partial {y^2}}}} \right),\quad \\&K \approx \frac{{{\partial ^2}w}}{{\partial {x^2}}}\frac{{{\partial ^2}w}}{{\partial {y^2}}} - {\left( {\frac{{{\partial ^2}w}}{{\partial x\partial y}}} \right)^2}.\end{split}$
在此需要指出的是, 方程(5)是经典F?ppl-von Karman板理论中对曲率的近似表达. 对于宏观板而言, 在小变形和小截面转动条件下, 可以忽略近似的曲率表达和精确表达的差别[25]. 若定义应力Airy函数为$F$, 可仿照经典板壳理论, 定义面内的二维应力为[25,29]
${N_{xx}} = \frac{{{\partial ^2}F}}{{\partial {y^2}}},\quad {N_{xy}} = - \frac{{{\partial ^2}F}}{{\partial x\partial y}},\quad {N_{yy}} = \frac{{{\partial ^2}F}}{{\partial {x^2}}}.$
从变形能(1), 可得:
$\begin{split} &{N_{xx}} = {k_{\rm{b}}}\varepsilon _{xx}^0 + \left( {{k_{\rm{b}}} - {k_{\rm{g}}}} \right)\varepsilon _{yy}^0,\quad {N_{xy}} = {k_{\rm{g}}}\varepsilon _{xy}^0, \\& {N_{yy}} = \left( {{k_{\rm{b}}} - {k_{\rm{g}}}} \right)\varepsilon _{xx}^0 + {k_{\rm{b}}}\varepsilon _{yy}^0, \end{split} $
进一步可把面内应变表示为应力的形式:
$\begin{split}&\varepsilon _{xx}^0 = \frac{1}{\chi }{N_{xx}} - \frac{1}{\lambda }{N_{yy}},\quad \varepsilon _{yy}^0 = \frac{1}{\lambda }{N_{xx}} - \frac{1}{\chi }{N_{yy}},\\&\varepsilon _{xy}^0 = \frac{{{N_{xy}}}}{{{k_g}}}.\end{split}$
其中$\chi = {k_{\rm{g}}}\left( {2{k_{\rm{b}}} - {k_{\rm{g}}}} \right)k_{\rm{b}}^{{\rm{ - 1}}}$, $\lambda = {k_{\rm{g}}}\left( {2{k_{\rm{b}}} - {k_{\rm{g}}}} \right)( {k_{\rm{b}}} - $$ {k_{\rm{g}}} )^{ - 1}$. 此时面内变形能密度可写为
$\begin{split}&\frac{1}{2}{k_{\rm{b}}}{\left( {2J} \right)^2} - {k_{\rm{g}}}Q = \frac{1}{{2\chi }}{\left( {{N_{xx}} + {N_{yy}}} \right)^2} \\&- \frac{1}{{{k_{\rm{g}}}}}\left[ {{N_{xx}}{N_{yy}} - {{\left( {{N_{xy}}} \right)}^2}} \right].\end{split}$
把方程(6)代入 (9) 式, 有:
$\begin{split}{U_s} =\;& \int\nolimits_S \left[ \frac{1}{{2\chi }}{{\left( {\frac{{{\partial ^2}F}}{{\partial {y^2}}} + \frac{{{\partial ^2}F}}{{\partial {x^2}}}} \right)}^2}\right.\\&\left.- \frac{1}{{{k_{\rm{g}}}}}\left( {\frac{{{\partial ^2}F}}{{\partial {y^2}}}\frac{{{\partial ^2}F}}{{\partial {x^2}}} - {{\left( {\frac{{{\partial ^2}F}}{{\partial x\partial y}}} \right)}^2}} \right) \right]{\rm{d}}S .\end{split}$
为保证位移场的连续单值性, 应变场满足如下的完备性条件[25]: $2\dfrac{{{\partial ^2}\varepsilon _{xy}^0}}{{\partial x\partial y}} - \dfrac{{{\partial ^2}\varepsilon _{xx}^0}}{{\partial {y^2}}} - \dfrac{{{\partial ^2}\varepsilon _{yy}^0}}{{\partial {x^2}}} = K$, 并可等价的表示为应力函数的形式, ${\Delta ^2}F = - \lambda K$. 应力函数引入时, 需要在势能泛函中引入Lagrange乘子$l\left( {x, y} \right)$, 则(1)式可写为
$\begin{split} U =\;& \int\nolimits_S {\Bigg\{ {\frac{1}{2}{k_{\rm{B}}}{{\left( {2H} \right)}^2} + {k_{\rm{G}}}K + \frac{1}{{2\chi }}{{\left( {\frac{{{\partial ^2}F}}{{\partial {x^2}}} + \frac{{{\partial ^2}F}}{{\partial {y^2}}}} \right)}^2}} } \\& - \frac{1}{{{k_{\rm{g}}}}}\left[ {\frac{{{\partial ^2}F}}{{\partial {x^2}}}\frac{{{\partial ^2}F}}{{\partial {y^2}}}\left. { - {{\left( {\frac{{{\partial ^2}F}}{{\partial x\partial y}}} \right)}^2}} \right]} \right.\\&+ l\left( {x,y} \right)\left( {{\Delta ^2}F + \chi K} \right) \Bigg\}{\rm{d}}S.\\[-15pt] \end{split} $
对(11)式进行复杂但直接的计算, 并识别Lagrange乘子后(具体计算方法可参考文献[25]), (11)式可化为
$\begin{split} U = \;&\int\nolimits_S {\Bigg\{ {\frac{1}{2}{k_{\rm{B}}}{{\left( {2H} \right)}^2} + {k_{\rm{G}}}K} } + \frac{1}{{2\chi }}\left( {\frac{{{\partial ^2}F}}{{\partial {x^2}}} + \frac{{{\partial ^2}F}}{{\partial {y^2}}}} \right)^2 \\ &- \frac{1}{{{k_{\rm{g}}}}}\left[ {\frac{{{\partial ^2}F}}{{\partial {x^2}}}\frac{{{\partial ^2}F}}{{\partial {y^2}}}} - {{\left( {\frac{{{\partial ^2}F}}{{\partial x\partial y}}} \right)}^2} \right] \\ &+ \frac{1}{2}\left[ \frac{{{\partial ^2}F}}{{\partial {y^2}}}{{\left( {\frac{{\partial w}}{{\partial x}}} \right)}^2} + \frac{{{\partial ^2}F}}{{\partial {x^2}}}{{\left( {\frac{{\partial w}}{{\partial y}}} \right)}^2}\right.\\ &\left.- 2\frac{{{\partial ^2}F}}{{\partial x\partial y}}\frac{{\partial w}}{{\partial x}}\frac{{\partial w}}{{\partial y}} \right] \Bigg\}{\rm{d}}S.\\[-15pt] \end{split} $
为了考虑温度对结构的影响, 假设边界轴向外力和热应力的势能为
$\begin{split} W =\;& \int\nolimits_S {q\left( {x,y,t} \right){\rm{d}}S} + \frac{1}{2}\int\nolimits_S {\left[ {\left( {N_{xx}^0 - N_{xx}^{\rm{T}}} \right){{\left( {\frac{{\partial w}}{{\partial x}}} \right)}^2}} \right.} \\ & \left. { + 2N_{xy}^0\frac{{\partial w}}{{\partial x}}\frac{{\partial w}}{{\partial y}} + \left( {N_{yy}^0 - N_{yy}^{\rm{T}}} \right){{\left( {\frac{{\partial w}}{{\partial y}}} \right)}^2}} \right]{\rm{d}}S. \end{split} $
在此, $q\left( {x, y, t} \right)$为面内载荷; $N_{xx}^0, N_{yy}^0$为边界$x, y$轴方向的初始二维应力, $N_{xx}^{\rm{T}}, N_{yy}^{\rm{T}}$为边界$x, y$轴方向的二维热应力. 边界上的轴向应力, 压为负, 拉为正. 当不考虑温度对弹性参数的影响时, 均匀温度场下的热应力为[29]
$N_{xx}^{\rm{T}} = {k_b}\varepsilon _{xx}^{\rm{T}} = {k_{\rm{b}}}\alpha T,\quad N_{yy}^{\rm{T}} = {k_b}\varepsilon _{yy}^{\rm{T}} = {k_b}\alpha T.$
在此, $\varepsilon _{xx}^{\rm{T}}$, $\varepsilon _{yy}^{\rm{T}}$为热应变, $\alpha $为热膨胀系数. 把$w, F$作为独立变量, 构造Lagrange函数, $L = {U_{\rm{e}}} + {U_{\rm{s}}} - $$ W$, 并对其变分, 令$\delta L = 0$[30], 可得考虑均匀温度场产生的热应力影响的受力平衡方程为
$\begin{split} &{k_{\rm{B}}}{\nabla ^4}w = - \left( {\frac{{{\partial ^2}F}}{{\partial {y^2}}} + N_{xx}^0 - N_{xx}^{\rm{T}}} \right)\frac{{{\partial ^2}w}}{{\partial {x^2}}} - \bigg( {\frac{{{\partial ^2}F}}{{\partial {x^2}}}} \\& { + N_{yy}^0 - N_{yy}^{\rm{T}}} \bigg)\frac{{{\partial ^2}w}}{{\partial {y^2}}} + 2\left( {\frac{{{\partial ^2}F}}{{\partial x\partial y}} + N_{xy}^0} \right)\frac{{{\partial ^2}w}}{{\partial x\partial y}} + q, \\ &{\nabla ^4}F = - \chi \left[ {\frac{{{\partial ^2}w}}{{\partial {x^2}}}\frac{{{\partial ^2}w}}{{\partial {y^2}}} - {{\left( {\frac{{{\partial ^2}w}}{{\partial x\partial y}}} \right)}^2}} \right]. \\[-15pt] \end{split} $
由于变形能(12)和经典板的势能形式相似, $F, w$的边界条件和经典板理论中的相同. 对于四边铰支的板, 有:
$\begin{split} &w\left( {x = 0,a} \right) = 0,\quad \frac{{{\partial ^2}w}}{{\partial {x^2}}}\left( {x = 0,a} \right) = 0, \\ & w\left( {y = 0,b} \right) = 0,\quad \frac{{{\partial ^2}w}}{{\partial {y^2}}}\left( {x = 0,b} \right) = 0, \\ & \int\nolimits_0^a {\left[ {\frac{1}{\chi }\frac{{{\partial ^2}F}}{{\partial {y^2}}} - \frac{1}{\lambda }\frac{{{\partial ^2}F}}{{\partial {x^2}}} - \frac{1}{2}{{\left( {\frac{{\partial w}}{{\partial x}}} \right)}^2}} \right]} {\rm{d}}x = 0, \\ &\int\nolimits_0^a {\left[ {\frac{1}{\lambda }\frac{{{\partial ^2}F}}{{\partial {y^2}}} - \frac{1}{\chi }\frac{{{\partial ^2}F}}{{\partial {x^2}}} - \frac{1}{2}{{\left( {\frac{{\partial w}}{{\partial y}}} \right)}^2}} \right]} {\rm{d}}x = 0. \end{split}$
在此需要注意的是, 当前的理论具有4个独立的力学参数: 面内和面外的力学参数是独立的！该理论没有人为的引入二维结构的厚度, 很好的解决Yakobson悖论引发的不确定性. 这是和已有的单层二维纳米结构力学模型最明显的区别, 例如文献[10, 11, 26]中的板壳理论. 事实上, 本文的二维结构热弹耦合非线性板理论具有4个独立的力学参数和一个热学参数(热膨胀系数). 对于不同的单层二维结构, 只要通过原子计算(例如FDTB或MD), 得到5个参数, 就可以得到相应的力学模型. 因此, 该理论可以用来描述其他单层原子结构的力学性质, 例如石墨烯、h-BN等. 下节将通过平衡方程(15)来对比研究石墨烯和${\rm{Mo}}{{\rm{S}}_2}$的力学性质.
3.算例及讨论
本节通过方程组(15), 来讨论单层${\rm{Mo}}{{\rm{S}}_2}$和石墨烯的热弹耦合力学性质. 假设有四边铰支, 边长为分别为$a$和$b$的单层片(如图1所示), 此时的边界条件为(16)式. 由于方程组(15)是非线性的, 除了个别的情况, 精确的解析解是很难得到的[25,27,29]. 故本文使用Galerkin方法[30]近似求解方程(15). 方程(15)和(16)的形式和经典板壳理论相似(但力学参数和板厚无关), 可以引入和F?ppl-von Karman板相似的近似解析解[29]:
$\begin{split}& w = \sum\limits_{m,n = 1}^\infty {{\eta _{mn}}\sin \frac{{n{\rm{\pi }}x}}{a}} \sin \frac{{m{\rm{\pi }}y}}{b}, \\& F = \sum\limits_{j,k = 1}^\infty {{\xi _{jk}}\sin \frac{{j{\rm{\pi }}x}}{a}} \sin \frac{{k{\rm{\pi }}y}}{b}. \end{split} $
为了简化讨论, 仅取一项, 并令${\eta _{11}} = \eta, \quad {\xi _{11}} = \xi $. 把(17)式代入(15)式, 等号两边同时乘以$\sin \left( {{\text{π}}{a^{ - 1}}x} \right) $$ \sin \left( {{\text{π}}{b^{ - 1}}y} \right)$, 并在矩形期间积分(Galerkin积分[30]), 整理得到:
${k_1}\eta - {k_2}\eta \xi - {k_0}q = 0,\quad \xi = - {k_3}{\eta ^2},$
在此
$\begin{split} & {k_1} = \frac{{{k_{\rm{B}}}{{\rm{\pi }}^4}{{\left( {{b^2} + {a^2}} \right)}^2}}}{{{a^4}{b^4}}} + \frac{{{{\rm{\pi }}^2}\left( {N_{xx}^0 - {k_{\rm{b}}}\alpha T} \right)}}{{{a^2}}} \\&\qquad+ \frac{{{{\rm{\pi }}^2}\left( {N_{yy}^0 - {k_{\rm{b}}}\alpha T} \right)}}{{{b^2}}}, \\ & k_2= \frac{{32{{\rm{\pi }}^2}}}{{3{a^2}{b^2}}},\quad {k_3} = \frac{{16{a^2}{b^2}\chi }}{{3{{\rm{\pi }}^2}{{\left( {{a^2} + {b^2}} \right)}^2}}},\quad {k_0} = \frac{{16}}{{{{\rm{\pi }}^2}}}. \end{split}$
方程组(18)可简化为一个方程:
${k_1}\eta + {k_2}{k_3}{\eta ^3} = {k_0}{\kern 1pt} q.$
(20)式即确定结构中心点的变形幅值和载荷关系的非线性代数方程.
这里需要指出的是, 对于${\rm{Mo}}{{\rm{S}}_2}$的面内力学参数, 二维杨氏模量和Poisson比, 已经有比较系统的研究[13,31-33]. 面外的弯曲刚度${k_{\rm{B}}}$, 近期已经受到了关注[14-15]. 但是, 扭转刚度(Gauss 刚度)${k_{\rm{G}}}$, 迄今还没有相关报到. 如何确定Gauss刚度是一个需要进一步研究的问题. 但根据微分几何中的经典理论, Gauss-Bonnet积分公式[28], 有$\displaystyle \iint\nolimits_S {K{\rm{d}}S }= $$ 2{\text{π}} - \displaystyle \sum\nolimits_{i = 1}^n {{a_i}} - \displaystyle \oint\nolimits_C {{k}{\rm{d}}s}$. 其中${k}$为边界曲线的测地曲率, ${a_i}$为角点的外角. 若石墨烯片的边界在变形前后处于相同的平面内, 则$\displaystyle \oint\nolimits_C {{k}{\rm{d}}s} = 0$, 此时矩形的石墨烯片有$2{\text{π}} - \displaystyle \sum\nolimits_{i = 1}^n {{a_i}} = 0$, 故$\displaystyle \iint\nolimits_S {K{\rm{d}}S = }0$. 因此Gauss曲率不出现在平衡方程和边界条件中, 不影响本文的结果. 但当结构出现自由边界时, Gauss曲率刚度将出现在边界条件中, 并影响结构的受力变形, 我们将在后续研究中处理该问题. 在此取${\rm{Mo}}{{\rm{S}}_2}$的力学参数为[14-15,33]: $Y = 120\, {\rm{N}}/{\rm{m}}$, $\nu = $$ 0.23$, ${k_{\rm{{\rm B}}}} = 9.61\, {\rm{eV}}$. 为直观了解${\rm{Mo}}{{\rm{S}}_2}$的力学性质, 用石墨烯来对比研究. 石墨烯的力学参数为[19]: $Y = 340\, {\rm{N}}/{\rm{m}}$, $\nu = 0.165$, ${k_{\rm{{\rm B}}}} = 1.6\, {\rm{eV}}$. 热膨胀系数[17]: ${\rm{Mo}}{{\rm{S}}_2}$为$\alpha = 6.49 \times 1{0^{ - 5}}~{{\rm K}^{ - 1}}$, 石墨烯为$\alpha = $$ - 2.14 \times 1{0^{ - 5}}\, {{\rm K}^{ - 1}}$. 根据方程(2), ${\rm{Mo}}{{\rm{S}}_2}$有, ${k_{\rm{b}}} = $$ 792\, {\rm{eV}}/{\rm{n}}{{\rm{m}}^2}$, ${k_{\rm{g}}} = 610\, {\rm{eV}}/{\rm{n}}{{\rm{m}}^2}$; 石墨烯有, ${k_{\rm{b}}} = $$ 2814\, {\rm{eV}}/{\rm{n}}{{\rm{m}}^2}$, ${k_{\rm{g}}} = 1822\, {\rm{eV}}/{\rm{n}}{{\rm{m}}^2}$. 为便于讨论, 几何尺寸取为正方形, 即$a = b$. 根据方程(20)可得${\rm{Mo}}{{\rm{S}}_2}$和石墨烯, 在不同情况下的中心点变形和参数的关系, 如图2—图6所示.
图 2 $a = b = 6\;{\rm{nm}}$, $N_{xx}^0 = N_{yy}^0 = 0$时, 结构在两个不同温度下的载荷变形幅值曲线
Figure2. Loads-response curves with two temperatures for $a = b = 6\;{\rm{nm}}$ and $N_{xx}^0 = N_{yy}^0 = 0$.

图 6 给定边界轴向力和温度条件下的载荷、几何尺寸及变形幅值曲面
Figure6. Loads- dimensions-response surfaces with the given stretching stresses and the temperature.

先讨论温度对结构的影响. 令$N_{xx}^0 = N_{yy}^0 = 0$, 从方程(20)可得石墨烯和${\rm{Mo}}{{\rm{S}}_2}$在相同载荷下的变形幅值, 如图2所示. 从图2可以看出: 对于0 K条件下的正方形片, 载荷较小时石墨烯的变形大于${\rm{Mo}}{{\rm{S}}_2}$; 但伴随载荷的增大, 石墨烯的变形将小于${\rm{Mo}}{{\rm{S}}_2}$的变形. 这说明, 小变形时石墨烯的刚度小于${\rm{Mo}}{{\rm{S}}_2}$, 而大变形时石墨烯的刚度大于${\rm{Mo}}{{\rm{S}}_2}$. 出现此反常现象的原因是: 1)${\rm{Mo}}{{\rm{S}}_2}$的弯曲刚度远远大于石墨烯, 而小变形时面外的抗弯刚度是抵抗变形的主要来源, 因此${\rm{Mo}}{{\rm{S}}_2}$的变形小于石墨烯; 2)伴随变形的增大, 面内的刚度逐渐成为抵抗变形的主要因素, 而石墨烯的面内弹性模量远远大于${\rm{Mo}}{{\rm{S}}_2}$, 这使得大变形时石墨烯的变形小于${\rm{Mo}}{{\rm{S}}_2}$. 但是, 当温度出现小幅升高时, 例如仅升高$1\;\rm K$(如图2所示), 石墨烯的变形就一致地小于${\rm{Mo}}{{\rm{S}}_2}$. 这表明两种材料的力学性质对温度很敏感.
来看轴向力的影响, 此时令$T = 0$. 从图3可知, 虽然轴向的拉力可显著强化石墨烯和${\rm{Mo}}{{\rm{S}}_2}$的刚度, 但变形趋势和图2一致. 相对于轴向力的影响, 温度对石墨烯和${\rm{Mo}}{{\rm{S}}_2}$力学性质的影响更复杂和有趣. 在有限温度条件下(T < 1000 K), 石墨烯的热膨胀系数$\alpha < 0$, 是热缩冷胀材料; 而${\rm{Mo}}{{\rm{S}}_2}$的热膨胀系数$\alpha > 0$, 是热胀冷缩材料[16-18]. 这使得温度强化了石墨烯的刚度, 但降低了${\rm{Mo}}{{\rm{S}}_2}$的刚度. 图2和图4显示, 相同条件下轻微的温度升高就使得石墨烯的变形幅值小于${\rm{Mo}}{{\rm{S}}_2}$. 但是, 当同时考虑温度和轴向力的影响时, 情况会变得复杂. 如图5所示, 温度较低时, 在$N_{xx}^0 = N_{yy}^0 = 7.0\;{\rm{nN}}/{\rm{nm}}$的轴向拉力作用下${\rm{Mo}}{{\rm{S}}_2}$的变形, 小于在$N_{xx}^0 = N_{yy}^0 = $$ 0.1\;{\rm{nN}}/{\rm{nm}}$的轴向拉力作用下石墨烯的变形; 但当温度超过临界值时, 变形幅值大小的顺序将颠倒. 从方程(20)的系数表达式(19)还可以看出, 变形幅值和结构的几何尺寸紧密相关, 如图6所示. 以上结果表明, 温度和边界轴向力对两种二维材料的力学性质均有显著影响. 在此需要指出的是, 本文仅考虑轴向拉力对两种材料力学性质的影响. 当边界出现轴向压力时, 结构可能屈曲. 此时温度的提高将增加石墨烯的抗屈曲能力, 但会降低${\rm{Mo}}{{\rm{S}}_2}$的抗屈曲能力. 屈曲将使得结构的力学行为, 特别是动力学行为变得异常复杂[34-35]. 我们将另文讨论该问题. 在此还需要指出的是, 本文根据唯像的连续介质力学建立了单层${\rm{Mo}}{{\rm{S}}_2}$的热-弹耦合模型, 并依据原子计算结果取定力学和热膨胀系数. 尽管新的模型建立在严谨的逻辑推理之上, 但该理论的正确性最终需要实验和量子计算的检验. 相对于被广泛深入研究的石墨烯而言, ${\rm{Mo}}{{\rm{S}}_2}$还缺乏系统量子计算和实验研究[1,22,23]. 特别是多场耦合条件下的力学问题还需要深入研究.
图 3 $a = b = 6\;{\rm{nm}}$, $T = 0\;\rm K$时, 在两个不同边界拉力下的载荷变形幅值曲线
Figure3. Loads-response curves with two edge stretching stresses for $a = b = 6\;{\rm{nm}}$ and $T = 0\; \rm K$.

图 4 $a = b = 6\;{\rm{nm}}$时, 在两个不同温度和边界载荷下的载荷变形幅值曲线
Figure4. Loads-response curves with two edge stresses and two temperatures for $a = b = 6\;{\rm{nm}}$.

图 5 $a = b = 6\;{\rm{nm}}$时, 给定边界轴向力条件下的载荷、温度及变形幅值曲面
Figure5. Loads-temperatures-response surfaces with the given stretching stresses for $a = b = 6\;{\rm{nm}}$.

4.结　论
假设面内张拉和剪切变形能, 与面外弯曲和扭转变形能相互独立, 本文建立了单层${\rm{Mo}}{{\rm{S}}_2}$的热弹耦合非线性板理论, 并通过该理论对比研究了单层${\rm{Mo}}{{\rm{S}}_2}$和石墨烯的力学性质. 结果显示, 在边界不移动条件下, 边界的预加轴向力和有限温度对结构的力学性质具有显著影响. 因此可以通过边界预加轴向载荷或环境温度来改变二维纳米结构的力学性质. 例如, 通过提高环境温度来弱化${\rm{Mo}}{{\rm{S}}_2}$的刚度, 或强化石墨烯的刚度. 此外, 本文提出的热弹耦合非线性板理论不依赖单层二维结构的厚度, 可以方便的应用于描述其他单层二维纳米结构的力学性质.
相关话题/力学 结构 计算 材料 纳米
  • 领限时大额优惠券,享本站正版考研考试资料!
    大额优惠券
    优惠券领取后72小时内有效，10万种最新考研考试考证类电子打印资料任你选。涵盖全国500余所院校考研专业课、200多种职业资格考试、1100多种经典教材，产品类型包含电子书、题库、全套资料以及视频，无论您是考研复习、考证刷题，还是考前冲刺等，不同类型的产品可满足您学习上的不同需求。 ...
    考试优惠券 本站小编 Free壹佰分学习网 2022-09-19
  • In<sub>1+</sub><i><sub>x</sub></i>Te化合物的结构及热电性能研究
    摘要:InTe化合物中In+的孤对电子引发的晶格非谐性振动使得其具有本征极低的热导率,因此被认为是一种具有潜力的中温热电材料.然而,较低的电输运性能使得InTe的热电性能不高.在本工作中,采用熔融、退火结合放电等离子活化烧结工艺制备了一系列In1+xTe(x=0,0.001,0.003,0.005, ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 第一性原理研究Mg掺杂对LiCoO<sub>2</sub>正极材料结构稳定性及其电子结构的影响
    摘要:LiCoO2作为商业化最早的锂离子电池正极材料,至今仍受到许多研究人员的广泛关注.高电压下LiCoO2面临着严重的容量衰减和性能下降等问题,实验上通常采用体相元素掺杂以稳定LiCoO2在高电压下的晶体结构,从而提高其电化学性能.Mg元素掺杂被认为是一种能够提高LiCoO2高电压循环稳定性的有效 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 溶剂热制备铬掺杂硫化锌和硫化纳米结构和磁性能
    摘要:以乙醇胺和乙二胺为混合溶剂,通过简单溶剂热法,用S,ZnO和CdO为源,成功制备了不同量Cr掺杂ZnS和CdS半导体纳米结构.X-射线衍射测试表明,纳米结构ZnS和CdS具有纤锌矿结构.扫描电子显微镜给出了不同铬含量的ZnS和CdS的形貌.用电子能量散射谱观察到产物的成分为Cr,Zn,Cd和S ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 大面积二维磁性材料的制备及居里温度调控
    摘要:当前,尽管集成电路制造工艺水平不断提高,但受到量子效应的限制,器件尺寸的缩小使业界遇到了可靠性低、功耗大等瓶颈,微电子行业延续了近50年的“摩尔定律”将难以持续.2004年二维材料—石墨烯的问世,为突破集成电路的功耗瓶颈带来了新的机遇.由于低维特性,二维材料在一层或者几层原子厚度中表现出丰富多 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 二维材料中贝里曲率诱导的磁性响应
    摘要:二维材料中由贝里曲率诱导的新型磁学响应是近年来的新兴领域.这些二维材料所表现出的磁学特性及量子输运与贝里曲率直接相关,而贝里曲率又与晶体的对称性、电子的轨道磁性、自旋轨道耦合以及磁电效应等息息相关.研究这些新型磁性响应一方面有益于研究不同量子效应间的耦合作用,另一方面可探索量子效应在电子与信息 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 层状磁性拓扑材料中的物理问题与实验进展
    摘要:层状磁性材料与拓扑材料的交汇点同时结合了二者的优势,形成了在最小二维单元下同时具有磁序和拓扑性的材料体系,即层状磁性拓扑材料.这类材料的电子结构中可能存在狄拉克点、外尔点、节线等具有螺旋性或手性的拓扑电子态,同时涵盖了绝缘体、半金属和金属等的材料分类,导致新物性、新现象成为可能,因此引起了广泛 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 二维磁性材料及多场调控研究进展
    摘要:二维磁性材料是二维材料家族的新成员,其在单原胞层厚度依然保持长程磁序且易受外场调控,这为二维极限下的磁性以及其他新奇物理效应的研究提供了理想的平台,又为低功耗自旋电子学/磁存储器件的研制开辟了新的途径,成为国际上备受关注的前沿热点.本综述首先系统介绍了近年来发现的各类本征二维磁性材料的晶体结构 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 基于二维材料的自旋-轨道矩研究进展
    摘要:在具有自旋-轨道耦合效应的材料中,电荷流能够诱导产生垂直于电流方向的纯自旋流,当其注入近邻的磁性层时,会对其磁矩产生自旋-轨道矩.自旋-轨道矩能够快速、高效地翻转磁矩,为开发高性能的自旋电子器件提供了一种极佳的磁矩操控方式.二维材料由于具有很多的优点,如种类丰富、具有多样化的晶体结构和对称性、 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 二维磁性材料的物性研究及性能调控
    摘要:以石墨烯和二硫化钼为代表的二维材料,由于具有良好的电学、热学、光学以及力学性质,近年来成为了科学界一大研究热点.而作为二维材料的分支,二维磁性材料由于具有磁各向异性、单层磁有序等特殊性质,特别是磁性还可借助多种物理场进行调控,使其具有丰富的物理特性和潜在的应用价值,逐渐受到研究者的普遍关注.本 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29
  • 轴子拓扑绝缘体候选材料层状<inline-formula><tex-math id="M1">\begin{document}${\bf{Eu}}_{
    摘要:二维磁性材料的研究推动了现代纳米电子器件的发展.寻找本征的具有磁性的层状材料,为探索研究新的二维磁性材料、制备二维电子器件提供了重要的材料基础.近来,本征二维反铁磁拓扑材料的发现引起了人们的广泛关注和兴趣.${m{EuIn}}_{2}{m{As}}_{2}$被预言是一种轴子拓扑绝缘材料 ...
    中科院物理研究所 本站小编 Free考研考试 2021-12-29

Free考研考试FreeKaoYan.Com
欢迎来到Free考研考试，"为实现人生的Free而奋斗"

© 2020 FreeKaoYan! . All rights reserved.