365天读文献，Day3，二维超材料的力学各向异性

超材料的力学性能不仅取决于结构的组成，其材料也起着重要的作用。本研究采用经典的分子动力学模拟方法和通用的粗粒化模型，研究了二维正弦弯曲材料的变形机制以及力学各向异性对材料力学性能的影响，特别是负泊松比的影响。结果表明：在轴向拉伸载荷作用下，材料的变形过程可分为弯曲和拉伸两个阶段，其中弯曲阶段接头的旋转导致了材料的拉胀行为。此外，拉伸行为在很大程度上取决于正弦韧带的振幅/波长比和轴向与横向之间的刚度比。更有趣的是，随着刚度比增加，负泊松比也在增加。研究结果揭示了拉胀超材料的变形机制，为设计具有可调负泊松比的二维拉胀晶格结构提供了有益的指导。

Background（背景）

超材料是合理设计的多尺度复合材料系统，由于其广泛的应用潜力，其表现出大多数传统材料所没有的材料特性。虽然超材料通常是由他们的结构为主，其组成的影响也不能忽视。然而，由二维材料制成的超材料的微观结构-性能关系和成分-性能关系仍然是一个悬而未决的问题。二维材料的力学各向异性将显著影响相关超材料的稳定性和性能可调性。

Methods and models（方法和模型）

1、具有弯曲正弦梁的二维模型

2、2D材料的通用粗粒分子动力学（CG-MD）模型

相互作用可以分为三类，即键的相互作用，控制面内的杨氏模量，角的相互作用，定义面内剪切模量，和二面角的相互作用，影响面外弯曲刚度。这种通用CG-MD模型的优点之一是可以通过改变相互作用来调整力学性能。

3、分子动力学模拟

软件：LAMMPS

环境：1fs步长；1K；NPT系综；应变速率1x107s-1

Results and discussions （结果与讨论）

1、拉伸应力应变响应

如图2(a)所示，变形历史可大致分为两个阶段，即弯曲为主阶段和拉伸为主阶段。两个阶段的力学响应本质上是不同的。

在弯曲主导阶段，轴向组由于其初始弯曲形状而经历拉伸和弯曲，而横向组经历由轴向组弯曲引起的纯弯曲，导致非常规横向膨胀，其显示负泊松比。在以弯曲为主的阶段，局部应变远低于整个试样的名义应变，变形主要归因于“刚体转动”，应力缓慢增加。同时，沿y方向的横向膨胀沿着与沿x方向的轴向应变沿着几乎呈线性关系。请注意，在此阶段，可忽略面外屈曲，如图2（b）所示。然而，当轴向组进一步“拉直”超过0.2的应变时，横向组倾向于屈曲以释放其累积的弯曲能量。应变低于0.2时，整个结构的面外位移很小。相应地，横向群收缩，并且整个样品经历从负泊松比到正泊松比的急剧转变。在图S1和S2中，可以看出，当应变大于0.2时，整个结构的平面外位移迅速增加，其中整个结构的形状从加载方向为正弦曲线。

此外，随着轴向逐渐变直，而不是弯曲，伸展开始占主导地位。因此，如图2（a）所示，名义应力急剧增加。

通过图3(a)的简图中的应力分布图，进一步证实了在两个不同阶段中看到的变形模式。由于独特的结构，横向组不能承受沿着x方向的载荷，而轴向组承受沿沿着x方向的最大载荷。此外，沿着x方向的力不仅在轴向组内引起张力，而且还产生弯曲，这引起接合部的旋转，这也有助于拉直横向组。。当轴向应变超过0.2时，由于轴向组的拉直形状，沿x方向沿着的载荷增量只能产生较小的弯矩。同时，轴向组在沿x方向沿着的载荷下倾向于被拉伸而不是被拉直。

2、能量分析

图4中的能量分析结果显示了总能量变化，分为键能、角能和二面角能。从图4（a）可以看出，键能和角能的变化是总势能变化的主要贡献者，而二面角能的变化至少比键能和角能的变化小一个数量级，表明面内变形占相对主导地位。相比之下，如图4（b）和（c）所示，在拉伸主导阶段（>0.2应变），所有三项的势能都急剧增加，并且大部分能量变化来自轴向组。如图4（b-d）所示，这两个阶段之间的转变反映在键能和角能图中势能的急剧下降，以及横群二面角能的急剧增加。

3、振幅/波长参数的影响

从图5（a）中可以看出，曲线可以分为两个阶段，这与图2（a）中观察到的结果一致。有趣的是，随着振幅/波长比 A 的减小，从负泊松比到正泊松比的临界应变也相应减小，并且这一变化与早期观察到的面外屈曲现象密切相关。当 A=0.1 时，应变-应变关系表现出较强的非线性特征。图5（b）展示了低应变条件下泊松比与 A 之间的关系，其中所给出的泊松比是通过对应变小于0.05的数据点进行线性拟合得到的。可以观察到，随着 A 的增加，泊松比的大小逐渐增大，并趋近于0.63，这与文献中给出的值0.65非常接近。

波长 l和带宽 w在所提出的二维晶格结构的力学性能中也起着重要作用，但这两个几何参数对泊松比的影响较为有限。关于波长 l和带宽 w的影响的相关结果和讨论，详见图 S3。

4、力学各向异性效应

此部分非我研究重点，不过多阐述

Conclusions（结论）

本文系统研究了所提出的二维晶格结构的二维材料，通过几何因素和调谐机械各向异性对拉胀结构的影响，采用了CG-MD模拟方法。研究结果表明，所提出的二维晶格结构的变形过程可以分为两个阶段：弯曲主导阶段和拉伸主导阶段。在弯曲主导阶段，泊松比恒定且为负值；而在拉伸主导阶段，泊松比变为正值。拉胀性现象可以归因于轴向组和横向组的协同变形，轴向组上的张力产生的弯矩引起轴向/横向连接处的旋转，进而拉直横向组。在弯曲主导阶段和拉伸主导阶段之间的过渡伴随面外屈曲现象，导致泊松比从负值转变为正值。

此外，本研究还系统探讨了几何因素对泊松比的影响。研究发现，随着振幅/波长比A的增加，泊松比也随之增大，A值越大，表明结构的曲率越大；当A达到0.225时，泊松比趋于饱和。

最后，研究进一步分析了形成晶格的底层二维材料的机械各向异性对结构拉胀性的影响。结果显示，泊松比的大小随着加载方向的变化以及横向杨氏模量刚度比的增大而增大。通过磷烯、石墨烯和二硫化钼的全原子MD模拟，验证了刚度比与泊松比之间的关系，表明这些二维材料的规律具有广泛的适用性。这一规律可解释为横向组在较低刚度下更容易拉直，与轴向组的协同变形共同作用的结果。

总之，这些研究成果为由各向异性二维材料构成的拉胀二维晶格结构的设计提供了宝贵的指导，有助于实现可调节的负泊松比。