不同接触零点对纳米压痕测试影响的有限元解析
马尧, 董晓龙, 李泽君, 沈林, 赵宏伟
吉林大学 机械科学与工程学院, 长春 130022
通信作者:赵宏伟(1976-),男,教授,博士生导师.研究方向:材料微观力学性能原位测试技术与应用. E-mail:hwzhao@jlu.edu.cn

作者简介:马尧(1992-),男.研究方向:材料微观性能与3D打印制造.E-mail:myao1992@sina.com

摘要

为了有效估计由于粗糙度因素造成纳米压痕测试结果的偏差,对纯铝的纳米压痕进行了有限元仿真研究。利用二维轮廓模型,采用模型轮廓最大高度Ry作为评定样品表面粗糙度的参数,并引入轮廓的波峰与波谷之间的水平距离D和波谷的中垂线与相邻凸面切线之间的夹角θ仿真发现:相比光滑表面,在不同的接触零点,测试结果存在差异:接触零点为粗糙表面波峰时,最大载荷测试结果偏小;接触零点为粗糙表面波谷时,最大载荷测试结果偏大;压头压向波谷,压头棱边与凸起接触为接触零点时,测试结果引起的偏差相对较小。

关键词: 纳米压痕; 表面粗糙度; 接触零点; 有限元仿真
中图分类号:TG156 文献标志码:A 文章编号:1671-5497(2014)05-1366-05
Finite element simulation of the effect of different initial contact on nano-indentation
MA Yao, DONG Xiao-long, LI Ze-jun, SHEN Lin, ZHAO Hong-wei
College of Materials Science and Engineering, Jilin University, Changchun 130022,China
Abstract

To effectively estimate the test deviation of nano-indentation due to surface roughness, a finite element simulation of nano-indentation is carried on pure aluminium. A two dimension model is adopted and the maximum heightRy of the model contour is regarded as the assessment of the sample surface roughness parameters. ParameterD represents the horizontal distance between the peaks and valleys of the contour. The angleθ between the valley's midperpendicular and the tangent of adjacent convex is introduced. Compared with the smooth surface, for different initial contact, simulation results show obvious diversity: for initial contact at the peak of the rough surface, the maximum load is smaller; for initial contact at the valley of the rough surface, the maximum load is bigger; for initial contact at the edge of indenter and the bulge of the contour, the result is relatively close to the smooth one.

Keyword: nano-indentation; surface roughness; initial contact; finite element method(FEM)
0 引言

纳米压痕技术(Nano-indentation),又称为深度敏感压痕技术,是在特定形状的刚性压头上施加一定的载荷,压入被测样品的表面,同时通过位移传感器采集压入样品的深度信号。纳米压痕测试可以在纳米尺度下测量出材料的力学性能,如硬度、弹性模量、断裂韧性、应变硬化效应以及蠕变行为等[ 1, 2, 3, 4, 5, 6]。在纳米压痕测试中压头压入的深度一般控制在微纳米的尺度下,在实际加载过程中会带来一些误差,特别是在压入深度测量当中存在着偏差。测量误差主要来自测试环境(温度、噪声等)、测试仪器(压头的形状及尖端缺陷、位移和载荷的分辨力、机架柔度等)、被测材料的性质(凹陷、凸起)和材料表面特征(粗糙度等)[ 7, 8]。对于商用的纳米压痕仪来说,载荷和位移分辨力、设备的柔度、噪声等可以看作是固定的,因此,被测材料的性质和表面特征对纳米压痕实验的精度影响尤为重要。Donnelly等[ 9]研究了多孔骨质材料的纳米压痕测试,发现材料表面粗糙度的增加会造成测量数据的分散性的增加,并降低测试精确度。Jiang等[ 10]发现由于材料实际表面粗糙不平,压头的尺寸、压入的深度以及样品的表面形貌都会对纳米压痕结果造成偏差。

有限元仿真方法是分析接触力学的一种有效手段[ 11]。在纳米压痕过程的仿真中,有限元法可以不受实验条件的影响,研究纳米压痕过程中被测试样的应力应变分布,并能够构造单参量变化的材料特性,得出仿真结果进行指导性的试验分析,更好地分析被测材料在压痕实验中的影响因素。Walter等[ 12, 13]分别通过二维和三维有限元模型探究了纳米压痕测试,发现过大的粗糙度对测试硬度和杨氏模量会引入偏差。Chen等[ 14]通过有限元软件ABAQUS模拟了真实表面的粗糙模型,分析了在纳米压痕测试中考虑表面粗糙度的参数识别。本文在前人的基础上,利用有限元软件ABAQUS对塑性材料纯铝的纳米压痕过程进行仿真,研究了不同接触零点对测试结果的影响规律。

1 理想表面的有限元仿真
1.1 模型建立

选用纯铝作为仿真被测材料。铝为典型的塑性材料,在室温下,杨氏模量为71 700 MPa,泊松比为0.33,密度为2.7 g/cm3。假设材料为均匀的、各向同性的。根据Bouzakis等[ 15]的研究结果,在其他因素相同的情况下,纳米压痕测试结果与被选压头(维氏压头和玻氏压头)尖端形状无关。本文采用理想形状的玻氏压头和深度控制模式,压深分别为200、300、400 nm。被测试样模型采用棱边长18 μm的正立方体,把压头与试样的接触区域进行了网格细化。将被测试样底面固定作为边界条件,同时给波氏压头施加垂直于试件表面的直线位移。图1(a)为仿真所用的玻氏压头模型,图1(b)为被测试样模型。

图1 有限元模型Fig.1 Finite element model

1.2 仿真结果分析

图2(a)(b)分别是最大压深为200、400 nm纯铝的应力分布的剖面图,可以看出应力影响区域随着最大压深的增大而增大。图3为压深200 nm的残余表面应力分布图。最大压深分别为200、300、400 nm的深度-载荷曲线图如图4所示,从图中卸载曲线可以看出,对于塑性材料纯铝的压痕测试,卸载后变形恢复不明显。

图2 应力分布图Fig.2 Stress distribution

图3 残余表面应力分布图Fig.3 Residual stress distribution

图4 深度-载荷曲线图Fig.4 Depth-load graph

2 粗糙表面的有限元仿真

为检测样品表面粗糙度对纳米压痕测试的影响,Jiang等[ 10]采用预制的连续直凹槽来模拟粗糙度的效果,Berke等[ 16]引入正弦函数来模拟样品二维粗糙度形貌,Chen等[ 14]分别采用了单水平的正弦函数和多水平的正弦函数来模拟样品的二维粗糙度形貌,仿真结果证实了以上模型均能在纳米压痕测试中较好地反映真实样品的表面形貌。本文采用轮廓最大高度 Ry作为评定样品表面粗糙度的参数,将被测表面看作为由若干个连续的波峰和波谷组成的光滑曲线的集合。轮廓波峰与波谷之间的幅值为Ry,定义轮廓相邻的波峰与波谷之间的水平距离为D,图5所示。当 D hmax,或 Ry hmax时,材料表面近似为一个平面,对仿真的结果并不会有明显的影响,因此仿真中所讨论的 D Ry hmax三个参数的变化量是在一定的尺度范围内的。

图5 粗糙表面示意图Fig.5 Schematic diagram of the rough surface

接触零点的确定在纳米压痕测试中是关键的影响因素之一,压头与样品的初始接触位置会对压痕测试结果有着不同的影响[ 14, 15, 16]。在ISO14577[ 17]中推荐两种方法来确定接触零点。方法1:通过拟合函数来确定接触零点(例如二次多项式外推);方法2:通过接触刚度和测试载荷的首次增加量来确定接触零点。本文选用第二种方法。考虑玻氏压头的形状尺寸,玻氏压头棱边与中心线的夹角选为77.05°,引入波谷的垂线与相邻凸面切线之间的夹角 θ,分三种情况分析压头与样品初始接触位置,如图6所示。

图6 压入位置示意图Fig.6 Schematic diagram of indent position

2.1 压头尖端与粗糙表面高峰接触

图6(a)所示,设定压头尖端与粗糙表面波峰接触时即为接触零点。当 D=100 nm, Ry=50 nm, hmax=200 nm时,接触零点在不同表面波峰的深度-载荷曲线图如图7所示。从图中可以看出,接触零点为粗糙表面波峰时,粗糙表面在相同位移下的载荷测试结果比光滑表面的测试结果小。由纳米压痕测试的硬度公式 H=Pmax/A,说明压头初始接触在粗糙表面波峰时,测试得到的材料硬度值偏小。此效应可以通过在载荷-位移曲线上设定新的初始点来削弱[ 14]

图7 与粗糙表面波峰接触Fig.7 Contact at the peak

图8为接触零点在不同粗糙度表面波峰仿真得出的深度-载荷曲线图( hmax=300 nm, D=300 nm)。从图中可以看出,相同压入深度下,光滑表面的载荷比粗糙表面的大,粗糙度 Ry的值越减小,其载荷值越接近光滑表面的载荷。

图8 不同粗糙度表面的深度-载荷曲线图Fig.8 Depth-load graph of different Ry

2.2 压头压向波谷,其棱边与突起接触处为接触零点

当玻氏压头尖端沿粗糙表面的波谷法线垂直压入,波谷的垂线与相邻凸面切线之间的夹角 θ<77.05°时,如图6(b),压头的尖端还未与粗糙表面底谷接触,其一条棱边就已经与波谷相邻的凸起处发生了接触,载荷开始增加,此时设定压头的棱边与突起接触处为接触零点。当 D=100 nm, Ry=50 nm, hmax=200 nm时,深度-载荷曲线图如图9所示。从图中可以看出,相比光滑表面,此种接触状态下的深度-载荷曲线的偏差较小,两条曲线的吻合性好。

图9 棱边与粗糙表面凸起接触Fig.9 Contact at the edge of indenter and the bulge

2.3 压头尖端与粗糙表面的波谷接触

当玻氏压头尖端沿粗糙表面波谷的法线压入,且 θ>77.05°时,设定接触零点在材料表面轮廓的波谷处,如图6(c)所示。由于 θ的大小与 DRy的取值有关,Ry值不变,通过D的值来改变θ的大小,测得当 D=300 nm, Ry=50 nm时, θ=79°;当 D=500 nm, Ry=50 nm时, θ=83°。 θ值均大于77.05°,对这两个模型进行了仿真分析。

hmax=200 nm时,接触零点在不同表面波谷时仿真得出的深度-载荷曲线图如图10所示。从图中可以看出,相比光滑表面的深度-载荷曲线,接触零点在粗糙表面波谷时的最大载荷值比光滑表面的大,说明在此种情况下测得的硬度值偏大。值得注意的是,与光滑表面的曲线相比, θ=79°模型的曲线在深度为20 nm左右时,载荷显现出明显的增加;而 θ=83°模型的曲线在深度为50 nm时载荷才出现明显的增加。推测其原因为:当压头向下移动时,对于 θ=79°的模型,在压深为20 nm时,压头的一条棱边已与波谷相邻的凸起处发生了接触;对于 θ=83°的模型,压深50 nm左右时,压头的一条棱边才与波谷相邻的凸起处发生了接触。

图10 与粗糙表面波谷接触Fig.10 Contact at the valley

3 结论

通过对纯铝试样粗糙表面的纳米压痕测试过程进行仿真,发现当模型轮廓最大高度 Ry和最大压入深度 hmax相同时,在不同的初始压入位置,测试结果存在显著的不同:接触零点为粗糙表面波峰时所得的硬度值比光滑样品的测试值要小;接触零点为粗糙表面低谷时所得的硬度值偏大;压头压向低谷,压头棱边与突起接触处设为接触零点时,测试结果引起的偏差相对较小。

材料表面粗糙度对纳米压痕的测试结果有显著的影响,越光滑的测试表面越能够精确地反映出材料的真实纳米力学性能。本文通过有限元仿真分析,发现在不同的接触零点,纳米压痕测试结果显示出了不同的效果,仿真结果有助于辨识由于粗糙度因素造成测量参量的偏差,也为合理设计压痕试验提供了指导。

The authors have declared that no competing interests exist.

参考文献
[1] 文玉华, 周富信, 刘曰武. 纳米材料的研究进展[J]. 力学进展, 2001, 31(1): 47-61.
Wen Yu-hua, Zhou Fu-xin, Liu Yue-wu. Research progress of nano materials[J]. Advances In Mechanics, 2001, 31(1): 47-61. [本文引用:1] [CJCR: 1.226]
[2] 张泰华, 杨业敏, 赵亚博, . MEMS材料力学性能的测试技术[J]. 力学进展, 2002, 32(4): 545-562.
Zhang Tai-hua, Yang Ye-min, Zhao Ya-bo, et al. Mechanics properties testing technology of MEMS materials[J]. Advances in Mechanics, 2002, 32(4): 545-562. [本文引用:1] [CJCR: 1.226]
[3] 黎明, 温诗铸. 纳米压痕技术理论基础[J]. 机械工程学报, 2003, 39(3): 142-145.
Li Ming, Wen Shi-zhu. Nano indentation theory basis[J]. Chinese Journal of Mechanical Engineering, 2003, 39(3): 142-145. [本文引用:1] [CJCR: 1.505]
[4] Lucas B N, Oliver W C. Indentation power-law creep of high-purity indium[J]. Metallurgical and Materials Transactions A, 1999, 30(3): 601-610. [本文引用:1]
[5] Poisl W H, Oliver W C, Fabes B D. The relationship between indentation and uniaxial creep in amorphous selenium[J]. Materials Research Soeiety, 1995, 10(8): 2024-2032. [本文引用:1]
[6] Lawn B R, Evans A G, Marshall D B. Elastic/plastic indentation damage in ceramics: the median/radial crack system[J]. Journal of the American Ceramic Society, 1980, 63(9-10): 574-581. [本文引用:1] [JCR: 2.107]
[7] 张泰华. 影响纳米压入测试结果的因素[J]. 实验力学, 2004, 19(4): 437-442.
Zhang Tai-hua. The influencing factors of nanoindentation test results[J]. Journal of Experimental Mechanics, 2004, 19(4): 437-442. [本文引用:1] [CJCR: 0.526]
[8] 张泰华, 杨业敏. 纳米硬度计及其在微机电系统中的应用[S]. 现代科学仪器, 2002(1): 32-37.
Zhang Tai-hua, Yang Ye-min. Nano-hardness tester and its application in MEMS[J]. Modern Scientific Instruments, 2002(1): 32-37. [本文引用:1] [CJCR: 0.3394]
[9] Eve Donnelly, Shefford P Baker, Adele L Boskey, et al. Effects of surface roughness and maximum load on the mechanical properties of cancellous bone measured by nanoindentation[J]. Journal of Biomedical Materials Research Part A, 2006, 77A(2): 426-435. [本文引用:1]
[10] Jiang Wu-gui, Su Jian-jun, Feng Xi-qiao. Effect of surface roughness on nanoindentation test of thin films[J]. Engineering Fracture Mechanics, 2008, 75(17): 4965-4972. [本文引用:2] [JCR: 1.413]
[11] 温卫东, 高德平. 接触问题数值分析方法的研究现状与发展[J]. 南 京 航 空 航 天 大 学 学 报, 1994, 26(5): 664-675.
Wen Wei-dong, Gao De-ping. Research situation and development of contact problems numerical methods[J]. Journal of Nanjing University of Aeronautics & Astronautics, 1994, 26(5): 664-675. [本文引用:1] [CJCR: 0.3373]
[12] Walter C, Antretter T, Daniel R, et al. Finite element simulation of the effect of surface roughness on nanoindentation of thin films with spherical indenters[J]. Surface & Coatings Technology, 2007, 202(4-7): 1103-1107. [本文引用:1]
[13] Walter C, Mitterer C. 3D versus 2D finite element simulation of the effect of surface roughness on nanoindentation of hard coatings[J]. Surface & Coatings Technology, 2009, 203(20-21): 3286-3290. [本文引用:1]
[14] Chen Zhao-yu, Stefan Diebels. Modelling and parameter re-identification of nanoindentation of soft polymers taking into account effects of surface roughness[J]. Computers & Mathematics with Applications, 2012, 64(9): 2775-2786. [本文引用:4]
[15] Bouzakis K D, Michailidis N, Hadjiyiannis S, et al. The effect of specimen roughness and indenter tip geometry on the determination accuracy of thin hard coatings stress-strain laws by nanoindentation[J]. Materials Characterization, 2002, 49(2): 149-156. [本文引用:2] [JCR: 1.88]
[16] Berke P, El Houdaigui F, Massart T J. Coupled friction and roughness surface effects in shallow spherical nanoindentation[J]. Wear, 2010, 268(1/2): 223-232. [本文引用:2] [JCR: 1.262]
[17] ISO/14577-2002. Metallic Materials-Instrumented Indentation Test for Hardness and Materials Parameters[S]. [本文引用:1]