CN108153956B - 基于分子动力学测定纳米材料压痕模量与硬度的计算方法 - Google Patents
基于分子动力学测定纳米材料压痕模量与硬度的计算方法 Download PDFInfo
- Publication number
- CN108153956B CN108153956B CN201711386231.0A CN201711386231A CN108153956B CN 108153956 B CN108153956 B CN 108153956B CN 201711386231 A CN201711386231 A CN 201711386231A CN 108153956 B CN108153956 B CN 108153956B
- Authority
- CN
- China
- Prior art keywords
- indentation
- hardness
- molecular dynamics
- calculating
- indentation modulus
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,包括:建立块体纳米材料和压头分子动力学模型;设置计算参数;采用分子动力学模拟的方法计算并输出纳米材料力与位移数据;将压入力与位移数据绘制出力位移曲线;取出卸载过程中的前1/3数据进行拟合,并得到拟合曲线的斜率;按照压头计算接触面积的方法得到接触面积;根据压痕模量和硬度的定义压痕模量和硬度值。本发明通过分子动力学模拟方法对纳米尺度下材料的压痕模量和硬度进行计算,能够克服已有实验方法的局限,并可以观察材料中的缺陷状况;不仅提高了计算精度和可重复性,而且降低了测试费用,节约成本。
Description
技术领域
本发明属于纳米材料测试方法技术领域,涉及基于分子动力学测定纳米材料压痕模量与硬度的方法。
背景技术
纳米材料被应用于制造微纳机电系统中的微/纳器件,其具有表面效应和尺寸效应因而表现出与宏观材料截然不同的力学性能,比如压痕模量和压痕硬度。纳米压痕作为最简单的力学测试方法已经被广泛应用,这种方式测得的力学性能在器件设计制造过程中有重要的应用前景。
目前纳米压痕技术在实验中已经被应用,并且这种方法是可靠而有效的。然而,这种实验的方法具有一定的局限性,具体表现在:(1)需要制造出适合纳米测试样品的压头和载物装置;(2)实验中纳米压痕仪可实现的最小载荷在微米级,无法达到纳牛;(3)实验的费用比较昂贵;(4)纳米尺度的材料在测试过程中受到环境的干扰大以致于实验的可重复性很差。
发明内容
为解决上述问题,本发明公开了一种基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,提高了计算方法的实用性和有效性。
为了达到上述目的,本发明提供如下技术方案:
基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,包括如下步骤:
步骤一,建立块体纳米材料和压头分子动力学模型;
步骤二,设置计算参数;
步骤三,采用分子动力学模拟的方法计算并输出纳米材料力与位移数据;
步骤四,将压入力与位移数据绘制出力位移曲线;
步骤五,取出卸载过程中的前1/3数据进行拟合,并得到拟合曲线的斜率;
步骤六,按照压头计算接触面积的方法得到接触面积;
步骤七、根据压痕模量和硬度的定义压痕模量和硬度值。
进一步的,所述步骤一中模型是采用LAMMPS软件建立的。
进一步的,所述步骤一中压头为球形。
进一步的,所述步骤二中计算参数包括:边界条件、势函数、系综选择、外部载荷类型、系统控温方式、加载和卸载速度、温度、时间步。
进一步的,所述步骤六中接触面积公式为:
上式中,A为接触面积,Ri为压头半径,hc为接触深度。
进一步的,所述步骤七中,
压痕模量公式为:
压痕硬度公式为:
与现有技术相比,本发明具有如下有点和有益效果:
本发明通过分子动力学模拟方法对纳米尺度下材料的压痕模量和硬度进行计算,能够克服已有实验方法的局限,并可以观察材料中的缺陷状况;不仅提高了计算精度和可重复性,而且降低了测试费用,节约成本。
附图说明
图1为本发明提供的基于分子动力学测定纳米材料压痕模量与硬度的计算方法流程图。
图2为本发明实施例建立的纳米压入分子动力学仿真模型。
图3为加载-保载-卸载过程的力位移曲线。
具体实施方式
以下将结合具体实施例对本发明提供的技术方案进行详细说明,应理解下述具体实施方式仅用于说明本发明而不用于限制本发明的范围。
本发明提供的基于分子动力学测定纳米材料压痕模量与硬度的计算方法,如图1所示,包括如下步骤:
第一步,建立纳米块体和球形压头分子动力学模型:
纳米块体由fcc结构的理想单晶Al组成,晶格常数a为0.405nm.通过LAMMPS软件建立了三维尺寸为70a*70a*60a的长方体分子动力学模型。图2为建立的模型,Al块体尺寸为70a*70a*60a,压头为刚性的金刚石压头,半径为4nm。本例中纳米块为长方体形,根据需要也可以建模其他形状的纳米块体,压头也可以采用圆锥形等其他形状。
第二步,设置仿真参数:
在LAMMPS软件中确定分子动力学模拟所必须的相关参数,包括温度为1K,边界条件是加载方向上为自由边界条件,其他方向上为周期边界条件,时间步为0.5飞秒,铝原子间的势函数为EAM势,压头与铝间的作用势采用Morse,系综在平衡阶段采用NVT,在压头工作阶段采用NVE,外部载荷类型为压入载荷、系统控温方式为速度标定法、加载和卸载速率均为50m/s,保载时间为10皮秒。
第三步,采用分子动力学模拟的方法计算并输出纳米材料力与位移数据:
采用lammps计算,输出加载-保载-卸载过程中压头受力和压头位移的数据。
第四步,将压入力与位移数据绘制出力位移曲线:
基于压入力与位移数据得到的力-位移曲线如图3所示。
第五步,取出卸载过程中的前1/3数据进行拟合,并得到拟合曲线的斜率:
卸载过程中的前1/3数据是指压头从最大压入深度到块体表面过程中前1/3的数据,对这1/3的数据进行二次函数y=a(x-b)c的拟合,拟合曲线的斜率指的是以过卸载前的一点(hmax,Pmax)作拟合曲线的斜率其中P为载荷,h为压入深度(即压头位移),Pmax为最大载荷,hmax为最大压头位移。
第六步,按照球形压头计算接触面积的方法得到接触面积:
球形压头计算接触面积的方法指的是通过已知的压头接触深度与压头半径来计算接触面积,接触面积公式为下式(1):
A=π(2×4×0.33-0.33×0.33)=7.95nm2。
第七步,根据压痕模量和硬度的定义压痕模量和硬度值:
压痕模量是斜率与接触面积的函数,其公式为下式(2):
压痕硬度是指用载荷除以面积,其公式为下式(3):
由此,得到纳米块体模型的压痕模量与压痕硬度。
本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也视为本发明的保护范围。
Claims (4)
1.基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,其特征在于,包括如下步骤:
步骤一,建立块体纳米材料和压头分子动力学模型;
步骤二,设置计算参数;
步骤三,采用分子动力学模拟的方法计算并输出纳米材料力与位移数据;
步骤四,将压入力与位移数据绘制出力位移曲线;
步骤五,取出卸载过程中的前1/3数据进行拟合,并得到拟合曲线的斜率;
步骤六,按照压头计算接触面积的方法得到接触面积;
接触面积公式为:
步骤七、根据压痕模量和硬度的定义压痕模量和硬度值;
压痕模量公式为:
压痕硬度公式为:
2.根据权利要求1所述的基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,其特征在于:所述步骤一中模型是采用LAMMPS软件建立的。
3.根据权利要求2所述的基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,其特征在于:所述步骤二中计算参数包括:边界条件、势函数、系综选择、外部载荷类型、系统控温方式、加载和卸载速度、温度、时间步。
4.根据权利要求3所述的基于分子动力学模拟的纳米材料压痕模量和压痕硬度的计算方法,其特征在于:所述步骤一中压头为球形。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711386231.0A CN108153956B (zh) | 2017-12-20 | 2017-12-20 | 基于分子动力学测定纳米材料压痕模量与硬度的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711386231.0A CN108153956B (zh) | 2017-12-20 | 2017-12-20 | 基于分子动力学测定纳米材料压痕模量与硬度的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108153956A CN108153956A (zh) | 2018-06-12 |
CN108153956B true CN108153956B (zh) | 2021-08-31 |
Family
ID=62464685
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711386231.0A Active CN108153956B (zh) | 2017-12-20 | 2017-12-20 | 基于分子动力学测定纳米材料压痕模量与硬度的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108153956B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110793983B (zh) * | 2019-11-21 | 2022-03-04 | 山东建筑大学 | 同步辐射原位测定金属玻璃有序化过程动力学曲线的方法 |
CN114323950A (zh) * | 2021-12-31 | 2022-04-12 | 华侨大学 | 基于力学试验结果的分子动力学模型验证方法 |
CN114334021A (zh) * | 2021-12-31 | 2022-04-12 | 华侨大学 | 基于透射电镜分析结果的分子动力学模型验证方法 |
CN115472248B (zh) * | 2022-09-23 | 2023-06-13 | 哈尔滨工业大学 | 一种CuZrAl非晶合金纳米压痕测试的分子动力学模拟计算方法 |
CN116822221B (zh) * | 2023-06-30 | 2024-02-23 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种基于相互侵入势的离散岩块间接触力计算方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009151215A (ja) * | 2007-12-21 | 2009-07-09 | Mitsui Chemicals Inc | 光学フィルムおよびその用途 |
CN101654222A (zh) * | 2009-09-04 | 2010-02-24 | 天津大学 | 基于分子动力学的三维纳米级切削加工模拟方法 |
CN101859336A (zh) * | 2010-05-20 | 2010-10-13 | 沈阳建筑大学 | 纳米加工中大规模分子动力学的并行优化方法 |
CN103047947A (zh) * | 2012-12-14 | 2013-04-17 | 北京工业大学 | 一种光刻技术和透射电子显微技术联合表征纳米薄膜微区形变的方法 |
CN105760598A (zh) * | 2016-02-15 | 2016-07-13 | 哈尔滨理工大学 | 一种基于分子动力学模拟的纳米材料塑性模量计算方法 |
CN106650021A (zh) * | 2016-11-24 | 2017-05-10 | 大连理工大学 | 一种脆性材料磨削过程建模仿真方法 |
JP2017105985A (ja) * | 2015-11-30 | 2017-06-15 | 東レ株式会社 | 光学フィルム製造用ポリエステルフィルム |
CN107356489A (zh) * | 2017-07-04 | 2017-11-17 | 湖北工业大学 | 一种基于纳米压痕试验测定花岗岩宏观力学性质的试验方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103207034B (zh) * | 2013-03-12 | 2015-02-18 | 清华大学 | 一种微机电系统桥膜结构残余应力的测试方法 |
CN104122152A (zh) * | 2014-07-22 | 2014-10-29 | 中国人民解放军装甲兵工程学院 | 基于仪器化Vickers压入O-P硬度确定材料维氏硬度的方法 |
-
2017
- 2017-12-20 CN CN201711386231.0A patent/CN108153956B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009151215A (ja) * | 2007-12-21 | 2009-07-09 | Mitsui Chemicals Inc | 光学フィルムおよびその用途 |
CN101654222A (zh) * | 2009-09-04 | 2010-02-24 | 天津大学 | 基于分子动力学的三维纳米级切削加工模拟方法 |
CN101859336A (zh) * | 2010-05-20 | 2010-10-13 | 沈阳建筑大学 | 纳米加工中大规模分子动力学的并行优化方法 |
CN103047947A (zh) * | 2012-12-14 | 2013-04-17 | 北京工业大学 | 一种光刻技术和透射电子显微技术联合表征纳米薄膜微区形变的方法 |
JP2017105985A (ja) * | 2015-11-30 | 2017-06-15 | 東レ株式会社 | 光学フィルム製造用ポリエステルフィルム |
CN105760598A (zh) * | 2016-02-15 | 2016-07-13 | 哈尔滨理工大学 | 一种基于分子动力学模拟的纳米材料塑性模量计算方法 |
CN106650021A (zh) * | 2016-11-24 | 2017-05-10 | 大连理工大学 | 一种脆性材料磨削过程建模仿真方法 |
CN107356489A (zh) * | 2017-07-04 | 2017-11-17 | 湖北工业大学 | 一种基于纳米压痕试验测定花岗岩宏观力学性质的试验方法 |
Non-Patent Citations (1)
Title |
---|
基于分子动力学的单层石墨烯纳米压痕模拟研究;贾瑞艳;《中国优秀硕士学位论文全文数据库工程科技I辑》;20160331(第03期);第B014-57页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108153956A (zh) | 2018-06-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108153956B (zh) | 基于分子动力学测定纳米材料压痕模量与硬度的计算方法 | |
Moser et al. | Strength and fracture of Si micropillars: A new scanning electron microscopy-based micro-compression test | |
Zhong et al. | Bending analysis and contact force modeling of soft pneumatic actuators with pleated structures | |
Gu et al. | Size-dependent deformation of nanocrystalline Pt nanopillars | |
Palacio et al. | Depth-sensing indentation of nanomaterials and nanostructures | |
Shen et al. | Reducing the rebound effect in micro-scale laser dynamic flexible forming through using plasticine as pressure-carrying medium | |
Dejeu et al. | Adhesion control for micro-and nanomanipulation | |
Hamdana et al. | Nanoindentation of crystalline silicon pillars fabricated by soft UV nanoimprint lithography and cryogenic deep reactive ion etching | |
Müller et al. | A hyperelastic model for simulating cells in flow | |
Tavakolian et al. | Pull-in instability of double clamped microbeams under dispersion forces in the presence of thermal and residual stress effects using nonlocal elasticity theory | |
Zhao et al. | Grain size effects on indentation-induced defect evolution and plastic deformation mechanism of ploycrystalline materials | |
Khatokar et al. | A study on improved methods in Micro-electromechanical systems technology | |
Sorgente et al. | Blow forming of AZ31 magnesium alloy at elevated temperatures | |
Lu et al. | Nanorobotic system iTRo for controllable 1D micro/nano material twisting test | |
Rinaldi et al. | Radial inner morphology effects on the mechanical properties of amorphous composite cobalt boride nanoparticles | |
Asher et al. | Bistable micro caps fabricated by sheet metal forming | |
Liang et al. | A variational size-dependent model for electrostatically actuated NEMS incorporating nonlinearities and Casimir force | |
Choi et al. | Hemispherical Arrays of Colloidal Crystals Fabricated by Transfer Printing | |
Leu et al. | Influence of punch shapes on the collar-drawing process of sheet steel | |
Wang et al. | Size effect in contact compression of nano-and microscale pyramid structures | |
CN110346103B (zh) | 一种mems悬臂梁的形变分析方法 | |
Ozaki et al. | Biaxial flexure testing of free-standing thin film membrane with nanoindentation system | |
Lin et al. | Novel microtensile method for monotonic and cyclic testing of freestanding copper thin films | |
Chen | An analysis of forming limit in the elliptic hole-flanging process of sheet metal | |
She et al. | A parametric study of compliant link design for safe physical human–robot interaction |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |