[发明专利]一种基于分子动力学模拟的纳米材料塑性模量计算方法

说明书

技术领域

本发明涉及纳米器件设计制造领域，具体是一种基于分子动力学模拟的纳米材料塑性模量计算方法。

背景技术

纳米材料广泛应用于纳米器件设计制造领域，例如逻辑门、量子器械中的导线、力反馈传感器。作为制造微纳机电系统中微/纳器件的材料，纳米材料所具有的表面效应和尺寸效应使其表现出与宏观材料截然不同的力学特性，比如塑性模量。纳米尺度下的塑性模量表现出了与宏观弹性模量性质相反的特性，在纳力反馈传感器的设计制造中有重要应用前景。因此，对纳米材料的塑性模量进行分析研究势在必行。

塑性模量的计算方法有很多，而基于扫描电子显微镜和透射电子显微镜等设备的物理实验方法是最可靠，也是最有效的方法。然而，倘若直接利用扫描电子显微镜等设备对纳米材料和器件进行力学性能实验，往往很难实现，甚至是困难重重。其原因具体表现为：一、需要制造出适合拉伸实验和夹具装夹的微小试件；二、需要精度极高的测量装置；三、需要制造出适用于拉伸实验的试验台；四、研发周期的长和实验数据的可重复性；五、高昂的试验费用。因此，物理实验方法令研究学者们望而却步。

分子动力学模拟方法被公认为本世纪以来除理论分析和实验观察之外的第三种科学研究手段，称之为“计算机实验”手段。采用分子动力学模拟方法不仅可以节约实验成本，有效、快捷地建立实验模拟环境，而且还可以观测瞬态原子运动状态。随着计算机速度、存储技术、硬件性能和并行处理技术的飞速发展，分子动力学模拟方法被广泛地应用于材料科学、生物物理和药品设计等科学研究中。目前，分子动力学模拟方法因其自身的突出特点已经成为研究纳米线力学性能机制的重要方法，并且正处在学术界和工业界全面重视的阶段。因此，采用分子动力学模拟方法计算纳米材料塑性模量是可行的。

发明内容：

本发明的目的在于提供一种提高了计算方法的实用性和有效性的基于分子动力学模拟的纳米材料塑性模量计算方法，以解决上述背景技术中提出的问题。

为实现上述目的，本发明提供如下技术方案：

一种基于分子动力学模拟的纳米材料塑性模量计算方法，包括以下步骤:

步骤一、建立纳米线分子动力学仿真模型；

步骤二、设置仿真参数；

步骤三、采用分子动力学模拟方法计算并输出纳米材料的应力应变数据；

步骤四、将应力应变数据绘制出应力应变曲线；

步骤五、找到应力应变曲线上的第一峰值点和第一下屈服点；

步骤六、按照塑性模量的定义求出塑性模量值。

作为本发明进一步的方案：步骤一中，建立纳米线分子动力学仿真模型是指基于LAMMPS软件建立三维尺寸为30a*6a*6a的长方体或者长为30a横截面直径为6a的圆柱体纳米线三维分子动力学仿真模型，a为所研究材料的晶格常数。

作为本发明进一步的方案：步骤二中，仿真参数是在LAMMPS软件中确定分子动力学模拟计算所必须的相关参数，包括应变率、温度、压强、时间步、势函数、系综选择、边界条件确定、外部载荷类型、系统控温和控压方法。

作为本发明进一步的方案：步骤五中，第一峰值点是指弹性极限点，第一下屈服点是连接后续非周期性波动的第一个最低临界点，也是第一次应力骤降幅度最大区段与下一个应力波动周期的拐点。

作为本发明进一步的方案：步骤六，塑性模量的定义是指应力应变曲线上在第一峰值点之后到第一下屈服点之间的一段呈线性函数关系的应力与应变的比值，塑性模量的公式为（1）-（3）；

式中：m是指在应力应变曲线中第一峰值点之后到第一下屈服点之间的一段呈现线性函数关系的的曲线段的数据数目，第一下屈服点下的应力，为第一峰值点下的应力，第一下屈服点下的应变，为第一峰值点下的应变。

作为本发明进一步的方案：由于计算机模拟输出的数据数量及其庞大并且数据提取和计算容易受到人为主观因素的干扰，实际计算求解塑性模量的值时是使用线性拟合的数据处理方法对塑性模量曲线部分进行线性拟合，拟合后获得的斜率K作为塑性模量的值，斜率K的计算公式如公式（4）所示；

（4）。

与现有技术相比，本发明的有益效果是：

本发明通过分子动力学模拟方法对纳米尺度下的塑性模量进行求解，既降低了计算费用，而且也提高了计算的精确度，同时相对于物理实验方法也提高了计算效率。

附图说明

图1为本发明所述基于分子动力学模拟的纳米材料塑性模量计算方法方法的流程图。

图2为纳米线分子动力学仿真模型。