基于XFEM研究纳米/微米划痕测试中划痕尖端尺寸效应与摩擦系数变化的学术报告
本报告旨在介绍Kwangmin Lee等研究者在期刊Tribology International第120卷（2018年）发表的题为“Scratch-tip size effect and change of friction coefficient in nano/micro scratch tests using XFEM”的原创性研究。该研究针对划痕测试领域长期存在的数值模拟难题，通过创新性地结合实验与先进的扩展有限元方法，系统探究了划痕尖端尺寸效应（STSE）及其对摩擦行为的影响。
本研究主要由韩国西江大学（Sogang University）和延世大学（Yonsei University）的研究人员合作完成。主要作者包括Kwangmin Lee, Karuppasamy Pandian Marimuthu, Chang-Lae Kim以及通讯作者Hyungyil Lee。论文于2018年1月在线发表于Tribology International期刊。
划痕测试是评估材料摩擦学性能的常用方法。然而，相比于发展成熟的压痕测试，划痕测试的数值分析因其复杂性而进展缓慢。理解划痕过程中的摩擦机制——尤其是摩擦系数的变化——是极具挑战性的任务，因为摩擦力受多种因素影响，包括失效模式、塑性变形、表面粗糙度以及材料特性等。其中，划痕尖端尺寸效应（Scratch-Tip-Size Effect, STSE）是决定脆性材料在划痕过程中发生屈服还是开裂的关键因素，是压痕尺寸效应（Indentation Size Effect, ISE）在划痕领域的延伸。
尽管已有一些关于划痕测试的实验研究，但能够综合考虑STSE、摩擦系数变化、多种失效模式、材料屈服与剪切强度、弹性回复（Elastic recovery）及堆积（Pile-up）等现象的数值模型十分有限。特别是，对于非晶态材料（如玻璃），描述其弹塑性行为的本构模型（如经典的von Mises或Tresca模型）往往与实验存在偏差，而考虑压力的Drucker-Prager（DP）模型虽被使用，却未能同时兼顾ISE现象。
因此，本研究的主要目标是： 1. 开发一个基于扩展有限元法（Extended Finite Element Method, XFEM）的划痕测试仿真模型，该模型需集成STSE、摩擦系数动态变化、多种失效模式、基于ISE校准的材料屈服强度与剪切强度、弹性回复及堆积效应。 2. 通过在钠钙玻璃上进行纳米和微米尺度的划痕实验，验证所建立的XFEM模型的有效性。 3. 利用模型解析划痕过程中摩擦系数的组成分量，并研究拉伸裂纹以及DP模型参数对摩擦系数的影响，从而深化对划痕损伤机理的理解。
本研究是一项集实验、材料参数校准、有限元建模与仿真分析于一体的综合性工作。其核心流程可分为以下三个主要部分：
1. 实验研究：材料性能评估与划痕测试 * 研究对象与样本量：研究对象为商业钠钙玻璃（Soda-lime glass）。所有测试均在室温条件下进行，每个划痕实验重复十次以确保结果的可靠性。 * 实验步骤与方法： * 纳米压痕与维氏硬度测试：首先，采用纳米压痕仪（Berkovich压头）获取材料在纳米尺度下的弹性模量、硬度等力学性能。同时，通过维氏硬度计获取材料在微米尺度下的硬度与断裂韧性。实验数据（如表1所示）明确揭示了钠钙玻璃的ISE现象：纳米硬度（7.8 GPa）显著高于微米硬度（5.4 GPa）。这些数据为后续DP模型参数的校准提供了关键基准。 * 纳米/微米划痕测试：为了研究STSE，研究团队选择了两种不同尺寸的划痕尖端。在纳米划痕测试中，使用半径R=2 μm (< 过渡半径Rt=0.218 mm)的球锥形金刚石尖端；在微米划痕测试中，使用半径R=0.5 mm (> Rt)的半球形不锈钢尖端。两种测试均采用渐进加载模式，具体参数见表2。实验过程中，同步记录法向载荷、切向载荷，并计算得到总体摩擦系数（μt）。 * 形貌与失效模式观察：划痕结束后，利用光学显微镜和扫描电子显微镜（SEM）详细观察划痕路径及其损伤形貌。结果显示，由于STSE，两种测试引发了截然不同的失效模式：纳米划痕（小尖端）主要产生屈服和随后的切屑（Chipping）失效；而微米划痕（大尖端）则在划痕尖端后方产生拉伸裂纹（Tensile cracks），属于典型的脆性断裂。
2. 有限元模型建立与材料模型校准 * 模型开发与创新方法： * 划痕测试XFEM模型：使用商业软件Abaqus 6.13建立了三维有限元模型（纳米划痕模型约含11万节点，微米划痕模型约含38万节点）。模型创新性地集成了多种复杂物理过程：1) 采用Drucker-Prager（DP）屈服准则来描述玻璃的压力依赖性塑性行为；2) 利用XFEM技术模拟微米划痕中观察到的拉伸裂纹萌生与扩展；3) 采用内聚力模型（CZM）定义损伤起始与演化；4) 实现了包含临界剪切应力限制的摩擦模型，以考虑材料发生整体滑移（Gross sliding）的条件；5) 通过用户子程序（fric_coef）实现了摩擦系数（对应于库仑摩擦）的非单调变化，以精确匹配实验观测的摩擦行为。 * DP模型参数校准——一种新方法：本研究的一个关键创新在于校准DP模型参数时，同时考虑了载荷-深度曲线和硬度-深度曲线，以反映ISE。传统方法仅拟合载荷-深度曲线。研究者建立了纳米压痕和微米压痕的FE仿真模型，通过反复调整DP模型中的零压力屈服强度（σc）和材料摩擦角（β），使仿真结果同时匹配从文献中提取的实验载荷-深度与硬度-深度数据（如图10，图11所示）。最终确定了适用于钠钙玻璃的线性DP模型：y = 1.757 + p * tan(20°) [GPa]。此外，材料的剪切强度基于Mohr-Coulomb准则和高压实验数据确定：τ = 0.495 + 0.652p [GPa]。
3. 摩擦系数分量解析方法与仿真分析 * 摩擦系数分量提取方法：研究者发展了一种结合理论与有限元分析的摩擦系数分量提取方法。总体摩擦系数μt被分解为界面摩擦系数μs（包括微凸体变形μd和粘着μa）和犁削摩擦系数μp。核心步骤如下：首先，建立包含库仑摩擦的FE模型（ITF模型），使其μt变化曲线与实验一致。然后，建立另一个库仑摩擦为零的FE模型（NITF模型）。在假设弹性回复和堆积效应不受库仑摩擦显著影响的前提下，NITF模型计算出的摩擦系数即近似为纯犁削贡献的μp。最后，根据公式μs = μt(实验) - μp(FE)，从实验总摩擦系数中减去FE获得的犁削分量，即可分离出界面摩擦分量μs。 * 验证与参数化研究：将校准后的材料属性和DP模型代入划痕XFEM模型进行仿真。结果显示，仿真与实验的法向载荷曲线以及摩擦系数变化曲线高度吻合（如图15，图16所示），验证了模型的准确性。利用验证后的模型，研究者：1) 成功解析了纳米划痕过程中μs和μp随划痕距离的演变（如图12所示），并将过程分为不同阶段进行机理阐释；2) 通过对比ITF和NITF模型中的弹性回复角ω和堆积高度，论证了所提摩擦分量提取方法的适用条件；3) 分析了微米划痕中XFEM模拟的拉伸裂纹，发现其对摩擦系数变化影响甚微；4) 进行了DP模型参数（σc和β）的敏感性分析，揭示了材料屈服特性对摩擦系数演变趋势的影响规律（如图20，图21所示）。
本研究的核心结论是，通过结合实验与先进的数值模拟技术，成功构建了一个能够精确模拟纳米/微米划痕测试复杂物理过程的集成化XFEM框架。该模型不仅成功复现了划痕尖端尺寸效应导致的截然不同的失效模式，还能精确预测摩擦系数的演变。
其科学价值在于： * 方法论贡献：提供了一套校准非晶材料DP模型参数的新思路（同时拟合载荷-深度与硬度-深度曲线），以及一种结合FE仿真与理论解析来分离摩擦系数分量的有效方法。 * 机理洞察：深化了对划痕过程中摩擦起源的理解，定量揭示了界面摩擦和犁削摩擦在不同损伤阶段的贡献变化，并明确了拉伸裂纹对宏观摩擦行为影响较小的新认识。 * 模型普适性：所开发的模型框架具有广泛的适用性，不仅适用于块体非晶材料，其考虑摩擦和剪切应力的特点使其尤其适用于模拟涂层材料的划痕行为，因为涂层的剥离往往由界面处的II型（剪切）断裂所引发。
本研究在附录中详细推导了球形划痕尖端下摩擦系数分量的理论解，并阐述了在有限元分析中应用摩擦分量提取方法所需满足的条件（即弹性回复和堆积效应在有无界面摩擦的模型间应相似）。这些内容为其他研究者理解和应用该方法提供了坚实的理论基础和实用指南。此外，论文也坦诚指出了当前模型的局限性，如未考虑切屑的形成机制、加载速率、划痕速度、环境湿度与温度等因素，并为未来的研究方向提出了建议，包括建立更完善的脆性材料延性断裂准则以及能够完全反映ISE的屈服模型。