边缘凸起式超显微硬度压痕的共聚焦测量

边缘凸起式超显微硬度压痕的共聚焦测量 Confocal microscopic measurement of ultramicro hardness considering indentation pile-up 1 first-author 石　伟（1972-），女，北京人，高级工程师，2006年于北京航空航天大学分别获得硕士学位，主要从事力学硬度精密计量测试仪器研制、校准及溯源技术研究及标准制修订等的研究。E-mail： swlltm@126.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848301&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848296&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848301&type= 18.89099121 28.04547119 石　伟 （1972-），女，北京人，高级工程师，2006年于北京航空航天大学分别获得硕士学位，主要从事力学硬度精密计量测试仪器研制、校准及溯源技术研究及标准制修订等的研究。E-mail： swlltm@126.com swlltm@126.com 2 corresp E-mail： 13811357882@163.com E-mail： 13811357882@163.com 张凯林（1994-），男，河北衡水人，硕士，2021年于北京航空航天大学获得硕士学位，主要从事光学计量技术及硬度测量技术的研究。E-mail：13811357882@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848313&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848309&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848313&type= 20.53198242 28.85673523 张凯林 （1994-），男，河北衡水人，硕士，2021年于北京航空航天大学获得硕士学位，主要从事光学计量技术及硬度测量技术的研究。E-mail： 13811357882@163.com 13811357882@163.com 3 中国航空工业集团公司 北京长城计量测试技术研究所，北京 100095 中国航空工业集团公司 北京长城计量测试技术研究所，北京 100095 Changcheng Institute of Metrology and Measurement， Aviation Industry Corporation of China， Ltd.， Beijing 100095， China Changcheng Institute of Metrology and Measurement， Aviation Industry Corporation of China， Ltd.， Beijing 100095， China 北京航空航天大学 仪器科学与光电工程学院， 北京 100191 北京航空航天大学 仪器科学与光电工程学院， 北京 100191 School of Instrumentation and Optoelectronic Engineering， Beihang University， Beijing 100191， China School of Instrumentation and Optoelectronic Engineering， Beihang University， Beijing 100191， China 东莞材料基因高等理工研究院 应力测量技术研究所， 广东 东莞 523808 东莞材料基因高等理工研究院 应力测量技术研究所， 广东 东莞 523808 Department of Stress Analysis， Centre of Excellence for Advanced Materials， Dongguan 523808， China Department of Stress Analysis， Centre of Excellence for Advanced Materials， Dongguan 523808， China 为了精确测量边缘产生凸起的超显微硬度压痕，采用高分辨率且能获得压痕三维形貌的激光扫描共聚焦显微镜，获得压痕边缘复杂的三维凸起结构。对压痕的三维形貌数据进行分析，利用压痕对角线方法提取压痕的4个角点，计算超显微硬度。然而，压痕对角线法并不能反映残余压痕的复杂结构，尤其是边缘凸起结构。在测量压痕三维形貌的基础上，为寻找更符合接触力学和压痕形成规律的压头与试样的接触面积，提出了基准平面的概念。通过对压痕面积函数进行微分，寻找最优的基准平面位置，既考虑压痕的边缘凸起，也能很好地反映接触面积。实验结果表明，通过接触力学规律寻找到的最优基准平面测量硬度，其测量误差小于±1.5%、稳定性小于1%。基准平面法可以稳定且准确地得到更符合硬度定义的硬度值。 To accurately measure the ultra-micro-hardness indentation by considering the pile-up， a high-resolution laser scanning confocal microscope was used to obtain the 3D morphology of the residual indent in this study. The four corners of indentation were then extracted to calculate the ultra-micro-hardness value using the indentation diagonal method. However， this traditional method cannot accurately reflect the complex structure of the residual indenter， in particular， the pile-up. Taking advantage of the knowledge on the 3D shape of the indent， a datum plane concept is proposed to determine the true contact area in accordance with the contact mechanics. It was found that the contact area determined by the datum plane at the inflexion point of the area function differentiated by the depth is a good indication of the contact area， even considering pile-up. The experimental results show that the hardness values measured by the optimal datum plane obtained by the contact mechanics law are <±1.5% error and <1% stable. The reference plane method can obtain a more stable and accurate hardness value that is more consistent with the definition of hardness. 共聚焦显微镜 凸起式 超显微硬度 接触力学 基准平面 confocal microscope pile-up ultra-micro hardness contact mechanics datum plane 1 引　言 1 随着科技的发展及新材料的不断涌现，产品的多样化发展和质量控制水平对硬度测量及校准提出了越来越高的要求，其中对新材料微小力及微压痕超显微硬度领域的测试问题尤为突出。微观硬度由于特别适用于微小、薄型试件及脆硬试件的测试，对金属组织、金属表面加工层、电镀层、硬化层、热处理试件、碳化试件、光纤维、碳纤维和陶瓷等材料的校准测试有着特殊的需求和研究意义。 超显微硬度压痕是介于显微压痕和纳米压痕之间的硬度压痕，超显微硬度压痕的特点为：其对角线长度在5~20 μm之间；载荷、加载速度及保持时间可控，载荷在10 gf~300 gf。 随着超显微硬度测量需求的增多，要求也越来越高。目前，常用的显微硬度计采用传统的光学成像来测量压痕大小，由于成像分辨率的限制，往往不能满足显微硬度压痕尤其是超显微硬度压痕的测量。另外，传统的光学显微镜仅能够得到压痕的二维投影数据，然后通过测量压痕对角线计算试样的硬度；但是在一些情况下，试样材料的不均匀或杂质的存在使得压痕边缘产生不规则的变化，呈现不规则的凸起结构，对测量结果的影响较大 ［ 1 1 ］ 。 在科研和工程实践中，对小试件尤其是薄膜、涂层等材料的显微硬度测量的常用方法是仪器化压入试验（通常被叫做纳米压痕仪）。和常规的硬度测量相比，仪器化压入技术可以记录整个压入过程中的试验力-位移曲线，并根据Oliver-Pharr方法利用卸载部分的载荷-压深关系曲线计算出等效接触面积，从而求出材料的压痕硬度 ［ 2 2 ］ 。仪器化压入试验的最大特点是不需要直接测定压痕大小，仅通过载荷-压深关系曲线和压头面积函数即可得到硬度和其它材料的力学参数，尤其适用于受限于光学成像分辨率的纳米硬度测量 ［ 3 3 ］ 。但是仪器化压入试验机相对于传统硬度仪的设备成本和维护成本都更高，而且对于相关技术人员的要求也高，其工业应用受到了一定程度的限制。除此之外，仪器化压入测量中的接触面积是通过压头面积函数和压入深度计算得到的，一般不会考虑压入过程中压痕周围所产生的凸起或凹陷，对于一些特定材料的硬度测量会存在误差 ［ 4 4 ］ 。因此，有必要针对超显微压痕的三维形貌进行精确测量以获得微观硬度，并在此基础上和仪器化压入测量得到的硬度值进行对比分析。由此可见，常用的微小压痕测量装置包括显微维氏测量装置和纳米压痕仪都不能满足边缘凸起式的超显微压痕的测量。为实现对角线小于20 μm的边缘凸起式超显微压痕的测量，本文采用分辨率更高的激光扫描共聚焦显微镜（Laser Scanning Confocal Microscopy，LSCM）对压痕进行测量，得到了压痕的全面三维形貌数据，更准确地表征材料的硬度。 在压痕等微观结构的形貌测量中，人们已经应用LSCM辅助观察压痕的形貌。Cheng ［ 5 5 ］ 等通过LSCM观察硅样品和碳化硅样品，比较压痕在微观结构形态上的不同。Li ［ 6 6 ］ 等研究球压头尺寸对布氏硬度的影响，利用LSCM观察球形压痕的三维形貌。Jong ［ 7 7 ］ 通过LSCM对压痕进行测量，并提出了通过体积参数计算洛氏硬度的模型。Hou ［ 8 8 ］ 等比较了光学显微镜、原子力显微镜和LSCM的测量结果的优缺点和不确定度，并提出了一种压痕处理方法来区分球形压痕的接触内点区域和外点区域。综上所述，目前共聚焦显微技术一是辅助观察不同材料压痕的结构形态，二是用于测量球形压痕；而对于边缘凸起不规则的维氏四棱锥压痕的测量研究却很少。 使用高分辨率的LSCM测量边缘凸起不规则四棱锥压痕的一个关键问题是基准平面的确定。本文对不同基准平面的判定标准开展研究，计算相对应的压痕有效区域和硬度，并与对角线测量方法得到的超显微硬度值进行对比，得到了符合硬度测量规律的基准平面判定方法。通过比较不同硬度值的硬度块的测量结果，验证了该基准平面的有效性。 2 测量原理 1 LSCM是在传统光学显微镜上加装激光扫描装置，对试样焦平面上的每一点进行扫描，逐点逐层测量试样上的点。物镜的衍射分辨率为： <math display="block" id="M1"><mi>σ</mi><mo>=</mo><mfrac><mrow><mi>k</mi><mi>λ</mi></mrow><mrow><mi>N</mi><mi>A</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848324&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848321&type= 13.71599960 7.61999989 ， 式中： <math id="M2"><mi>λ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848333&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848328&type= 1.77800000 2.79399991 为波长， NA 为物镜的数值孔径， k 为系数。在这个光路中，只有处于焦平面位置的光能穿过探测针孔，焦平面以外的光则被小孔过滤掉，使得非观察点的背景呈黑色，反差增大，成像清晰，可减小 k 的值，有效提高LSCM的分辨率 ［ 9 9 - 11 11 9-11 ］ 。此外，LSCM采用逐层扫描的设计，不仅具有较高的横向分辨率（可达到0.3 μm），还具有较高的纵向分辨率（可达到20 nm），所以各层扫描的图像也较为清晰。 LSCM测量的压痕为三维立体结构，且试样并不是理想的均匀材料，在压痕边缘会产生不规则的凸起，因此，确定三维压痕的边缘非常重要。 2.1　超显微硬度测量原理 2 超显微硬度属于压入硬度法，将顶部两相对面为136°的正四棱锥体金刚石压头用试验力压入试样表面，保持规定时间后，卸除试验力，测量试样表面压痕。超显微硬度通过测量压痕对角线长度 d 进而计算压痕表面积，最后求出压痕表面积上的平均压力，即有： <math display="block" id="M3"><mi>H</mi><mi>V</mi><mo>=</mo><mn mathvariant="normal">0.102</mn><mo>×</mo><mfrac><mrow><mn mathvariant="normal">2</mn><mi>F</mi><mo>⋅</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfrac><mrow><mn mathvariant="normal">136</mn><mo>°</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac></mrow><mrow><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0.189</mn><mtext> </mtext><mn mathvariant="normal">1</mn><mo>×</mo><mfrac><mrow><mi>F</mi></mrow><mrow><msup><mrow><mi>d</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848341&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848337&type= 69.17266846 11.59933376 ， 式中： HV 为超显微硬度； F 为压入载荷； d 为压痕对角线长度 ［ 12 12 ］ 。 2.2　压痕对角线法 2 为提高超显微硬度的测量准确性，充分发挥LSCM测量三维压痕的优势，测量三维压痕的角点位置，通过试验力与三微压痕的对角线来计算硬度。 以压痕最低点所在轴为平面的旋转轴，取矩形截面，该截面与压痕的交线即为压痕的截面曲线。三维压痕的截面如 图1 图1 所示，横轴为压痕点与轴线的距离 <math id="M4"><mi>X</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848351&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848346&type= 2.62466669 2.79399991 ，纵轴为压痕点的高度 <math id="M5"><mi>Z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848505&type= 2.37066650 2.79399991 ，提取截面中的极值点，极值点位置可从一阶微分 图2 图2 中得到，即微分快速变化并接近零值的点。极值点左边的部分为压痕部分（即压头部分和被测材料的接触部分），极值点右边的部分没有和压头直接接触，去除极值点右边与压痕无关的部分，实现压痕的提取，选择压痕截面的极值点作为压痕的边缘。充分利用LSCM在 <math id="M6"><mi>Z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848505&type= 2.37066650 2.79399991 轴高分辨率的特点，可以精确取到压痕的边缘最高点，作为压痕的边缘，之后进一步从压痕边缘中筛选出4个角点 <math id="M7"><mo stretchy="false">(</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msub><mo>,</mo><msub><mrow><mi>p</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848379&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848374&type= 20.99733353 4.65666676 ，即可根据超显微硬度的定义计算硬度值。 2.3　基准平面法 2 压痕产生过程如 图3 图3 所示。原始试样平面为 h 0 ，压头加载到最大载荷时的压痕深度为 h m ，此时的试样与压头的接触位置的高度为 h c 。当载荷卸载后压痕产生弹性回复，残余压痕各面更接近球形曲面，其中残余压痕的凸起高度为 h b 。 超显微硬度的传统计算中，假定压痕为理想的标准四棱锥，但是实际测试得到的压痕会由于塑性变形的不均匀性而偏离此假定（如 图3 图3 所示），而且压痕边缘往往会产生凸起或凹陷。此方法还采用最高边缘位置来确定压痕大小，但是实际的压头和被测材料的最后接触点（即 h c 所对应的接触点）往往低于最高边缘，这也意味着传统的基于压痕对角线的硬度测量方法和硬度定义本身有着不小的偏离：对于压痕边缘产生凸起的材料，测量得到的硬度值低于实际硬度值，而对于边缘产生凹陷的材料则会相反。综上所述，通过对角线计算硬度值的方法过于简单，忽略了压痕本身的形貌和边缘信息，所确定的硬度值与硬度的基本定义有所差别 ［ 13 13 ］ 。LSCM测量硬度压痕的方法可以获得压痕的三维数据信息，为进一步分析压痕接触面积提供数据支持，因此，尝试寻找更符合接触力学规律的压痕接触面积计算方法。 本文针对具有边缘凸起特征的压痕测量，基准平面定义为确定接触面积所对应的平面，即压头和被测材料最高接触点所对应的平面，见 图4 图4 。假设最优的基准平面为 P r 平面， P r 下方的区域为有效压痕区域，若 P r 上方仍有压痕数据，则 P r 上方的区域为打压试样过程中产生的无效凸起，即无效压痕区域。本文的目标就是求得符合接触力学规律的基准平面 P r 。 根据接触力学规律，均匀介质的接触面积函数应均匀变化，试样与压头接触部分随着压入深度 h 的增加，接触部分的面积应呈现一定的规律，即接触面积函数。当 h 到达基准平面以上的部分时，面积函数呈现不同的规律，所以可以通过微分方程确定接触面积增长变化的拐点，而该拐点即为最优基准平面的位置 ［ 8 8 ， 14 14 ］ 。 测量的压痕为三维立体结构，在压痕边缘会产生不规则的凸起。使用压痕点云分割方法将压痕区域分割出来。设定基准平面，逐层切削三维压痕模型，可得到不同高度对应的面积函数。基准平面与试样平面平行，表示如下： <math display="block" id="M8"><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>x</mi><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mi>y</mi><mo>+</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mi>z</mi><mo>-</mo><mi>h</mi><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848442&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848438&type= 39.28533173 4.65666676 ， 式中： <math id="M9"><mo stretchy="false">(</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>,</mo><msub><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848447&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848443&type= 15.74800014 4.48733330 为试样平面的法向量，以压痕尖端的点为零点， h 为基准平面与尖端的距离。寻找最优基准平面的问题可以转化为寻找最优基准平面高度的问题。 基准平面 h 的迭代公式为： <math display="block" id="M10"><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>h</mi></mrow><mrow><mi>n</mi></mrow></msub><mo>=</mo><msub><mrow><mi>h</mi></mrow><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><mi>α</mi></mtd></mtr><mtr><mtd><mi>α</mi><mo>=</mo><msub><mrow><mi>α</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mi>σ</mi><mo>×</mo><mi>i</mi></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848455&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848451&type= 25.56933212 12.36133385 ， 式中： <math id="M11"><mi>α</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848462&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848459&type= 1.94733346 2.79399991 为学习率， <math id="M12"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848469&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848466&type= 1.69333339 2.79399991 为衰减因子。随着迭代次数的增加，学习率不断下降，所以开始时应使用较大的学习率，以加快收敛速度，防止发生震荡；而在接近压痕基准平面位置的时候采用较小的学习率，以更准确更快地得到不同高度对应的压痕面积。 3 实验结果及分析 1 3.1　实验及压痕提取 2 测试试样选择硬度检定中常用的钢试样，硬度块选择硬度值在200~800 HV的不同硬度值钢硬度块，硬度块为硬度测量中常用的GCr15钢材料标准块。利用共聚焦显微镜测量不同硬度值的压痕形貌。显微压痕实验载荷采用500，300，200和100 gf，获得压痕对角线大于等于20 μm的压痕，从而获得试样的基准硬度。力加载装置为北京长城计量所的MV-1000显微硬度计，如 图5 图5 所示。超显微压痕实验载荷根据硬度的不同采用200，100，50，25和10 gf打压试样，可以得到5~20 μm的压痕。 数据采集设备采用奥林帕斯OLS-5000型LSCM（如 图6 图6 所示）。测量使用的物镜的数值孔径NA为0.95，测量视场为128 μm×128 μm，采集点数为1 024×1 024，测量的层析高度设置为0.12 μm。测量过程中首先要确定包含整个压痕形貌的扫描高度范围，然后通过激光逐点逐层扫描试样采集数据，获得各层的数据之后对所有数据进行包络拟合，得到压痕的三维数据集合，即测量区域范围内（包含压痕）的每一点的三维坐标。 激光扫描共聚焦显微镜测量得到的压痕为试样表面上每一点的三维数据，包括每一点的坐标，如 图7 图7 所示。压痕提取的实现过程首先进行滤波，去除压痕数据中明显偏离试样的离散点云数据。其次， 图7 图7 所示以压痕最底端的四棱锥尖端点 <math id="M13"><mi>P</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848501&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848497&type= 2.45533323 2.79399991 所在的 <math id="M14"><mi>Z</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848507&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848505&type= 2.37066650 2.79399991 轴为旋转轴，以垂直于压痕表面的平面作为旋转平面，旋转360°。每一个旋转平面与压痕点云数据的交线，提取该交线的极大值点即拐点作为分割点，分割点左边的数据即为压痕区域数据，分割点右边的数据即为试样平面区域数据。对压痕数据与旋转平面的每一个交线做分割操作 ［ 15 15 ］ ，即可提取完整压痕区域数据。 该方法提取的压痕区域数据为试样中整个压痕数据，如 图8 图8 所示。压痕区域包括压痕的内点和压痕的边缘点，通过提取压痕边缘确定压痕边缘上的4个角点，得到压痕对角线的长度，从而与传统的显微硬度测量方式一致，实现了小于20 μm的超显微硬度的测量 ［ 15 15 ］ 。 3.2　压痕对角线法测量结果 2 通过试样与垂直平面的截面提取压痕区域，通过提取边缘获得压痕4个角点的位置，进而得到压痕的对角线长度。通过 式（1） 式（1） 的施加载荷和对角线长度得到试样的超显微硬度，分别进行不同硬度以及不同载荷的超显微测量实验，可以看到实验结果与显微硬度结果基本一致。但在个别点有差别，当压痕较小时硬度变化较大，其原因主要包括压痕尺寸效应、材料不均匀和测量误差等 ［ 16 16 ］ 。 如 表1 表1 所示，显微硬度值与超显微硬度值的差别较小，可以大体表征试样的硬度值。造成该偏差的根本原因是测量方式不同， 图9 图9 （a）所示为显微硬度测量压痕的2D投影，利用像素之间的对比度来提取边缘进而测量对角线长度；共聚焦测量压痕的3D形状，截面高度对比如 图9 图9 （b）所示，提取极值点位置作为边缘。虽然都是通过对角线计算硬度值，但成像原理和分辨率的不同从根本上决定了显微硬度和超显微硬度的差别。 编号 载荷/g 显微硬度/HV 超显微硬度/HV 相对偏差/% 1 300 393 395.59 0.66 2 200 393 392.21 -0.20 3 100 393 408.23 3.87 4 50 393 398.62 1.43 5 25 393 391.85 -0.29 6 10 393 399.41 1.63 7 300 556 571.91 2.86 8 200 556 561.98 1.08 9 100 556 560.58 0.82 10 50 556 559.27 0.59 11 25 556 549.21 -1.22 12 500 829 832.34 0.40 13 300 829 838.93 1.20 14 200 829 835.95 0.84 15 100 829 852.28 2.81 16 50 829 861.28 3.89 17 25 829 865.36 4.39 18 10 829 873.68 5.39 3.3　基准平面法测量结果 2 测量的压痕为三维立体结构，由于试样并不是理想的均匀材料，压痕边缘会产生不规则的凸起，如 图7 图7 所示。受边缘的凸起和下沉以及打压过程中的弹性变形和塑性变形的影响，上述超显微测量中对角线计算压痕表面积的方法必然具有误差。因此，为准确获得有效的接触面积，设计基准平面的概念如 图4 图4 所示，并寻找最优的基准平面来准确获得压痕的有效接触面积。 实验采取的压头是夹角为136°的四棱锥，与试样接触的平面均为平面。根据接触力学的规律，一个受载的钢塑性材料包含发生塑形流动的区域和基于刚性假设的无变形区域，被四棱锥挤出的塑性流动区域的材料凸起在压头的两边，所以形成如 图6 图6 所示的边缘凸起，但凸起部分并不都与压头接触。假定压头为标准的四棱锥，那么与压头接触的面积应呈现一定的函数关系，称为接触面积函数，即有： <math display="block" id="M15"><mi>S</mi><mo>=</mo><mi>f</mi><mo stretchy="false">(</mo><msub><mrow><mi>h</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848557&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848552&type= 14.64733315 4.65666676 ， 式中： S 为压痕面积， h i 为基准平面的高度。 通过 公式（4） 公式（4） 逐步迭代截取不同高度的压痕，获得不同高度所对应的接触面积，从而拟合出接触面积函数。当截取高度为压痕的非接触部分时，接触面积会发生突变，这一高度为接触部分的临界节点，所对应的截取平面为最优的基准平面，即更符合接触力学规律的基准平面。 设定不同高度的平面截取压痕，得到高度与面积的函数曲线如 图10 图10 所示。可以看到，随着高度的上升，面积呈二次多项式增长，但存在一个拐点使得高度上升，面积的增长趋势发生变化，该拐点即为最优基准平面的位置。对该函数进行一阶微分，得到： <math display="block" id="M16"><mi>S</mi><mi mathvariant="normal">'</mi><mo>=</mo><mi>f</mi><mi mathvariant="normal">'</mi><mo stretchy="false">(</mo><msub><mrow><mi>h</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848566&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848562&type= 15.23999977 5.07999992 . 一阶微分函数如 图11 图11 所示。 图11 图11 中为面积函数散点图的一阶导数，可以很明显看到在高度小于拐点（2.38，196.3）的位置，面积函数的一阶导数基本呈现线性变化，该变化符合面积函数为二次多项式的变化规律，在大于高度 h 1 点时面积函数的一阶导数变化增大。 对面积函数求二阶导数，得到： <math display="block" id="M17"><mi>S</mi><mi mathvariant="normal">'</mi><mi mathvariant="normal">'</mi><mo>=</mo><mi>f</mi><mtext> </mtext><mi mathvariant="normal">'</mi><mi mathvariant="normal">'</mi><mo stretchy="false">(</mo><msub><mrow><mi>h</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848605&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848601&type= 18.20333290 5.07999992 . 二阶导数如 图12 图12 所示，可以看到在同样高度二阶导数的变化也较剧烈，说明该点的高度为符合接触力学变化规律的接触面积的基准平面高度。 分别选取不同硬度的硬度块在相同条件下进行重复测量。在同一硬度块上施加相同的载荷打压试样，在相同条件下使用共聚焦显微镜测量压痕的最优基准平面的位置和相对应的压痕的表面积，试验结果如 表2 表2 所示。对于不同硬度块进行多次测量，测量值与平均值的差值在±1.5%以内，不同硬度块的测量实验的相对标准差在1%以内。实验表明，最优基准平面测量得到的硬度值的误差小于±1.5%，且测量稳定性小于1%。实验结果表明，基准平面法可以较稳定地测量得到压痕的硬度。目前，关于小于20 μm的压痕，并没有硬度测量真值，而基准平面法从定义出发，且实验稳定性较好，该实验结果可以更好地反映微小压痕条件下的试样硬度。 编号 最优基准平面高度 h d /μm 压痕表面积 S/ μm 2 硬度值/HV 平均值/HV 相对偏差/% 相对标准差/% 200-1 2.74 218.20 229.15 226.76 1.06 0.60 200-2 2.77 221.00 226.24 226.76 -0.21 0.60 200-3 2.76 220.80 226.45 226.76 -0.12 0.60 200-4 2.76 221.90 225.33 226.76 -0.62 0.60 200-5 2.78 221.30 225.94 226.76 -0.35 0.60 200-6 2.78 219.80 227.48 226.76 0.33 0.60 600-1 1.52 73.17 683.34 676.13 1.07 1.03 600-2 1.51 73.61 679.26 676.13 0.46 1.03 600-3 1.52 72.89 685.97 676.13 1.45 1.03 600-4 1.54 74.96 667.02 676.13 -1.35 1.03 600-5 1.53 74.60 670.24 676.13 -0.87 1.03 600-6 1.52 74.16 674.22 676.13 -0.28 1.03 600-7 1.51 74.31 672.86 676.13 -0.48 1.03 1000-1 1.15 47.01 1 063.60 1 056.58 0.66 0.87 1000-2 1.18 47.94 1 042.97 1 056.58 -1.29 0.87 1000-3 1.15 47.65 1 049.32 1 056.58 -0.69 0.87 1000-4 1.14 47.55 1 051.52 1 056.58 -0.48 0.87 1000-5 1.15 46.73 1 069.98 1 056.58 1.27 0.87 1000-6 1.16 47.26 1 057.98 1 056.58 0.13 0.87 1000-7 1.17 47.14 1 060.67 1 056.58 0.39 0.87 通过计算压痕的最优基准平面得到接触表面积的方法虽然可以准确得到硬度值，但该方法操作复杂，且和传统的硬度测量结果偏差较大。因此，考虑到操作复杂度以及硬度结果的重复性，在实际工程中对于超显微压痕的测量建议仍采用四角点测量对角线，进而计算硬度值的方法。最优基准平面测量硬度值的方法作为更接近真实硬度值的理论方法，可在后续工作中进一步深入研究。 4 结　论 1 本文在LSCM测量得到的三维数据的基础上，提出了一种有别于传统测量方法的显微硬度测量方法。该方法充分发挥了LSCM高精度、高分辨率、较高的纵向分辨率和轴向分辨率的特点，通过对压痕三维形貌数据进行滤波和边缘提取，提取压痕角点，并计算超显微硬度值。基于共聚焦的超显微三维压痕测量过程如下：（1）选择载荷 F 打压试样，共聚焦测量压痕的三维数据；（2）提取压痕区域；（3）提取压痕的4个角点，计算压痕的对角线长度 d ；（4）根据载荷 F 和对角线 d 计算硬度值，即为试样的超显微硬度值。实验结果表明，超显微硬度与标准显微硬度的偏差大部分在±3%以内，可以有效表征试样的超显微硬度。 本文从硬度的定义出发，寻找更符合接触力学规律以及压痕产生过程的接触面积的最优基准平面。通过微分方法准确获得了接触面积所对应基准平面的位置，寻找到更符合硬度定义的硬度值。基准平面法得到的硬度值的测量误差小于±1.5%，测量稳定性小于1%，该值在理论上比传统对角线方法更接近硬度值的真值，对硬度真值的测量和硬度与其他力学性能的关系的研究具有重大意义。本文以GCr15钢材料作为试验样品，基准平面法可以准确地测量该样品的微小压痕硬度值。不过，该方法是以压痕边缘凸起为前提，所以理论上仅适用于边缘产生凸起的维氏四棱锥压痕，而对于其他材料的适用结果以及测试效果有待于进一步的研究。 参考文献 ESTEBAN B . Indentation Hardness Measurements at Macro- ， Micro-， and Nanoscale： A Critical Overview ［EB/OL］. Springerlink .com， 2017 . Indentation Hardness Measurements at Macro- OLIVER W ， PHARR G . An improved technique for determining hardness and elastic modulus using load and displacement sensing indentation experiments ［J］. Journal of Materials Research . 1992 ， 7 （ 6 ）： 1564 - 1583 . An improved technique for determining hardness and elastic modulus using load and displacement sensing indentation experiments OLIVER W ， PHARR G . Measurement of hardness and elastic modulus by instrumented indentation： advances in understanding and refinements to methodology ［J］. Journal of Materials Research ， 2004 ， 19 （ 1 ）： 3 - 20 . Measurement of hardness and elastic modulus by instrumented indentation： advances in understanding and refinements to methodology CHENG Y T ， CHENG C M . What is indentation hardness ［J］ . Surface and Coatings Technology ， 2000 ： 417 - 424 . What is indentation hardness FAN C ， ZHANG S ， ZHAO Q ， et al . Study on Indentation Experiments and Mechanical Behavior of RB-SiC ［C］. Matec Web of Conferences ， 2017 ： 128 . Study on Indentation Experiments and Mechanical Behavior of RB-SiC MA L ， SAM L ， SONG J . Investigation of brinell indentation diameter from confocal microscope measurement and fea modeling ［J］. Hardmeko ， 2007 ： 65 - 70 . Investigation of brinell indentation diameter from confocal microscope measurement and fea modeling JONG D K ， MOON C Y ， JEEYOUL R . Rockwell hardness modeling using indented volume ［J］. Applied Mechanics and Materials ， 2012 ： 152 - 154 . Rockwell hardness modeling using indented volume HOU X D ， BUSHBY A J ， JENNETT N M . Direct measurement of surface shape for validation of indentation deformation and plasticity length-scale effects： a comparison of methods ［J］. Measurement Science and Technology ， 2010 ， 21 ： 1 - 7 . Direct measurement of surface shape for validation of indentation deformation and plasticity length-scale effects： a comparison of methods 张一 ， 余卿 ， 张昆 ， 等 . 基于数字微镜器件的并行彩色共聚焦测量系统 ［J］. 光学 精密工程 ， 2020 ， 28 （ 4 ）： 859 - 866 . 基于数字微镜器件的并行彩色共聚焦测量系统 ZHANG Y ， YU Q ， ZHANG K ， et al . Parallel color confocal measurement system based on digital Micromirror devices ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 4 ）： 859 - 866 . （in Chinese） . Parallel color confocal measurement system based on digital Micromirror devices 李叶 ， 等 . 激光扫描共聚焦显微镜 ［J］． 实验室研究与探索 ， 2015 ， 34 （ 7 ）： 262 - 269 . 激光扫描共聚焦显微镜 LI Y ， et al .. Laser scanning confocal microscopy ［J］. Laboratory Research and Exploration ， 2015 ， 34 （ 7 ）： 262 - 269 . （in Chinese） . Laser scanning confocal microscopy 徐宝腾 ， 杨西斌 ， 刘家林 ， 等 . 高速扫描激光共聚焦显微内窥镜图像校正 ［J］. 光学 精密工程 ， 2020 ， 28 （ 1 ）： 60 - 68 . 高速扫描激光共聚焦显微内窥镜图像校正 XU B T ， YANG X B ， LIU J L ， et al . Image correction for high-speed scanning laser confocal microscopy ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 1 ）： 60 - 68 . （in Chinese） Image correction for high-speed scanning laser confocal microscopy ISO 6507-1：2018 Metallic materials-Vickers hardness test-Part 1： Test method [S]. Switzerland : International Standards Organization , 2018 . BECK C E ， HOFMANN F ， ELIASON J K ， et al . Correcting for contact area changes in nanoindentation using surface acoustic waves ［J］. Scripta Materialia ， 2017 ， 128 ： 83 - 86 . Correcting for contact area changes in nanoindentation using surface acoustic waves GAYLORD G ， GUILLAUME K ， JEAN-MICHEL B . A new method to determine the true projected contact area using nanoindentation testing ［J］. Comptes Rendus Mecanique ， 2015 ， 343 ： 410 - 418 . A new method to determine the true projected contact area using nanoindentation testing 陈华伟 ， 袁小翠 ， 吴禄慎 . 基于曲率突变分析的点云特征线自动提取 ［J］. 光学 精密工程 . 2019 ， 27 （ 5 ）： 1218 - 1228 . 基于曲率突变分析的点云特征线自动提取 CHEN H W ， YUAN X C ， WU L SH . Automatic extraction of point cloud feature lines based on curvature mutation analysis ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 5 ）： 1218 - 1228 . （in Chinese） Automatic extraction of point cloud feature lines based on curvature mutation analysis 金宏平 ， 陈建国 . 基于能量法计算材料的弹性模量和压痕硬度 ［J］. 机械工程材料 ， 2013 ， 37 （ 9 ）： 84 - 95 . 基于能量法计算材料的弹性模量和压痕硬度 JIN H P ， CHEN J J . Calculation of elastic modulus and indentation hardness of materials based on energy method ［J］. Materials for Mechanical Engineering ， 2013 ， 37 （ 9 ） ： 84 - 95 . （in Chinese） Calculation of elastic modulus and indentation hardness of materials based on energy method 10.37188/OPE.20212905.1014 TB932 A 1004-924X（2021）05-1014-10 石伟,张凯林,侯晓东.边缘凸起式超显微硬度压痕的共聚焦测量[J].光学精密工程,2021,29(05):1014-1023. 石伟,张凯林,侯晓东.边缘凸起式超显微硬度压痕的共聚焦测量[J].光学精密工程,2021,29(05):1014-1023. SHI Wei,ZHANG Kai-lin,HOU Xiao-dong.Confocal microscopic measurement of ultramicro hardness considering indentation pile-up[J].Optics and Precision Engineering,2021,29(05):1014-1023. SHI Wei,ZHANG Kai-lin,HOU Xiao-dong.Confocal microscopic measurement of ultramicro hardness considering indentation pile-up[J].Optics and Precision Engineering,2021,29(05):1014-1023. 国家“十三五”技术基础项目（No.JSJL2016205A014）；广东省引进创新团队项目（No.2016ZT06G025） 2020-10-14 2021-01-14 2021-05-15 2021-06-24 光学精密工程 5 29