利用分子动力学模拟研究纳米压痕作用下二氧化硅的变形行为和能量分布

《Surfaces and Interfaces》:Deformation behavior and energy distribution of silica under nanoindentation using molecular dynamics simulation

【字体: 时间:2026年02月23日 来源:Surfaces and Interfaces 6.3

编辑推荐:

  硅石在压缩载荷下的微观变形与能量分布机制研究。通过分子动力学模拟的纳米压痕方法,分析含预制裂纹的硅石试体在压缩过程中的表面形貌演化、力学性能变化及能量分布规律。结果表明:预制裂纹在压痕区发生渐进性变形与愈合,裂纹扩展长度影响基体杨氏模量、硬度及总能量消耗;压痕区外基体呈现弹性变形,非晶态硅石塑性应变能占比高于晶态结构,但不同形态硅石的断裂能量占比稳定在31.2%左右。该研究揭示了矿物形态、预制裂纹长度及位置对压缩粉碎过程中能量转化的影响规律,为高能效破碎技术开发提供原子尺度理论支撑。

  
张传|潘永泰|曹星健|尹芳茜|张永强|李少卓
中国矿业大学(北京)化学与环境工程学院,北京100083,中华人民共和国;

摘要

压缩是粉碎过程中的关键加载方法。然而,大多数研究集中在压缩加载对产品的宏观影响上,缺乏对矿物在压缩过程中的微观变形行为和能量分布模式的关注。本研究采用基于分子动力学的纳米压痕模拟来研究具有预裂纹的二氧化硅在压缩下的行为,分析了表面积、机械性能和能量分布的变化。结果表明,在加载过程中,压头影响区域内发生了剪切变形,预裂纹逐渐变形直至完全愈合。压头影响区域内预裂纹的扩展降低了基底的杨氏模量、硬度以及总体能耗。在压头影响范围之外,基底和预裂纹经历了弹性变形。与晶体模型相比,非晶模型表现出更高比例的塑性能量和更低比例的弹性能量。所有模型之间的断裂能量比例没有显著差异,平均为31.20%。本研究从原子层面探讨了矿物形态、预裂纹长度和预裂纹位置如何影响二氧化硅粉碎过程中的能量演变,为提高粉碎过程的能量效率和减少能耗提供了理论见解。

引言

粉碎包括破碎和研磨[1],是矿物加工工程中的基础过程,用于有效释放矿物。这种尺寸减小过程通过能量密集型的颗粒破碎来促进有价值矿物与脉石的分离,产生大量的新鲜表面[2]。因此,内部矿石变形行为、能量分布以及最终产品新产生的表面积之间的相互关系成为粉碎研究的关键焦点[3]。研究这些方面有助于优化能量利用,从而提高整体粉碎效率[4]。
目前最广泛使用的粉碎方法是机械粉碎,主要的力施加方式包括破碎、劈开、断裂、研磨和冲击[5]。每种力施加方式都涉及在特定时间内施加的压缩载荷。尽管纯压缩载荷在粉碎过程中的能量利用效率低于拉伸载荷,但研究表明[6],矿物在压缩载荷下会发生变形,导致不同程度的内部损伤。在粉碎过程中,这些损伤由于外力而引发并扩展裂纹,影响产品的粒径分布和形状。这些发现表明,压缩机制在矿物释放过程中仍然起着重要作用。
因此,已经进行了大量研究来了解矿物在压缩载荷下的破碎机制,主要分为两个方向:一是研究不同类型设备(如高压研磨辊[7]、颚式破碎机[8]和圆锥破碎机[9])在压缩条件下的破碎性能;二是采用简单的实验室实验,对单个颗粒或颗粒集合体进行压缩,例如单轴压缩试验[10]、活塞压力试验[11]、落锤试验[12]和Split Hopkinson压力棒试验[13]。这些努力从宏观层面阐明了压缩力如何影响产品的粒径分布和形态,显著推动了粉碎理论的发展。
然而,这些研究缺乏对矿物在压缩过程中微观尺度上的变形行为和能量分布的详细分析,限制了它们揭示驱动宏观粉碎现象的基本因素的能力。在材料科学中,基于原位扫描电子显微镜的双悬臂梁试验可以观察样品在压缩载荷下的裂纹扩展[14]。然而,这类实验的高复杂性和成本限制了其在矿物粉碎研究中的应用。
随着计算机技术的进步,分子动力学(MD)模拟提供了一种方法,可以研究外部载荷与材料内部行为之间的机制关系,而不受外部实验条件的限制。MD模拟具有低成本和高精度的优势,使其成为研究材料微观尺度变形和失效机制的有效工具[15]。其中,基于MD的纳米压痕模拟通过检查材料的原子结构来分析压缩载荷下的变形行为和机械性能。这种方法已广泛应用于陶瓷材料[16]、合金材料[17]、单晶材料[18]、石墨烯[19]和涂层材料[20]的研究。这种方法有效地研究了矿物在原子尺度压缩载荷下的变形行为和能量演变机制。
粉碎本质上是通过能量消耗来减小材料粒径和削弱机械性能的过程。这一过程需要大量的能量消耗,但其能量效率仍然很低。粉碎过程消耗了全球约1.8%的电力[21],但其能量效率仅介于1%到5%之间[22]。压缩是机械粉碎中的基本力施加方式。研究矿物在压缩载荷下的能量分布对于提高效率和减少能耗至关重要。这有助于节能和环境的可持续性。
作为火成岩、变质岩和沉积岩的主要成分,二氧化硅(SiO2)占地球地壳质量的约12.6%[23]。这种地质上普遍存在的化合物主要以两种不同的结构形式存在:晶体(c-SiO2)和非晶体(a-SiO2)。这种普遍存在的矿物存在于多种矿床中,包括但不限于黄铜矿、磁铁矿、赤铁矿、白钨矿和含金矿床,其存在显著影响下游加工方法。作为一种典型的脆性材料,二氧化硅有效地反映了矿物在压缩载荷下的变形行为和机械响应[24]。此外,二氧化硅卓越的物理化学稳定性为其在多个工业领域的应用奠定了基础,从传统的二氧化硅玻璃和陶瓷制造到先进的技术实现,包括半导体级硅的生产、太阳能转换的光伏电池制造以及电信系统中的光纤开发[25]。这些多样的应用模式需要系统地研究其在压缩载荷下的机械行为。
我们之前的研究发现[26]表明,二氧化硅玻璃在压缩下的能量效率在中观和微观尺度上存在差异。这种差异与材料本身的内在裂纹有关。然而,裂纹如何影响各种能量形式的分布以及压缩过程中新表面的生成机制尚不清楚。此外,应该注意的是,晶体二氧化硅的性质更接近实际矿石,也应纳入粉碎机制的研究中。
因此,本研究选择了不同形式的二氧化硅作为基底。在基底中引入了不同位置和长度的预裂纹,以模拟天然矿物中存在的缺陷。采用MD模拟进行纳米压痕试验,分析了基底在压头加载过程中的变形行为。计算了机械性能,包括杨氏模量(E)和硬度(H),并阐明了主要能量成分的分布,如塑性应变能(Up)、弹性能(Ue)和断裂能(Uc)。结果揭示了矿物在压缩载荷下的变形行为和能量分布机制。
本研究阐明了创新粉碎技术(包括高压电脉冲和微波辐照)具有更高效率的机制。通过展示二氧化硅在压缩力作用下的变形和能量分布,我们的发现为优化能量分配策略提供了理论指导,以在粉碎过程中最大化新表面积的生成。这项研究为减少能耗和提高粉碎效率提供了重要见解。

部分片段

纳米压痕模型开发

本研究使用LAMMPS软件包[27]进行所有MD模拟,而OVITO软件[28]用于可视化结果。如图1所示,c-SiO2模型由基底和压头组成。α-石英单元格是基底的基本结构组件。为了清晰起见,根据晶体取向原理定义了一个笛卡尔坐标系,其中α-石英单元格的[210]方向对应于x轴,[010]方向

纳米压痕过程中的典型变形行为

如图4所示,以CIN模型为例,选择了xz平面上的截面快照来展示基底在压头加载-卸载过程中的典型变形行为。此外,使用不同的颜色可视化了基底的von Mises剪切应变[46]。如图4(a)-(c)所示,在加载过程中,压头正下方发生了显著的剪切变形,表明即使在纯压缩载荷下,

结论

为了研究压缩诱导粉碎过程中的微观变形机制和能量耗散,选择了二氧化硅作为目标材料。对具有不同二氧化硅形态、预裂纹长度和预裂纹位置的纳米压痕进行了MD模拟。主要发现总结如下:
  • 1.
    在压缩载荷下,基底内部也观察到了拉伸断裂。压头的主要影响区域位于其正下方,最大尺寸为
  • CRediT作者贡献声明

    张传:方法论、概念化、软件、写作 - 原始草稿、数据管理、验证。潘永泰:概念化、验证、写作 - 审稿与编辑、资金获取。曹星健:方法论、软件、研究。尹芳茜:方法论、软件。张永强:验证、正式分析。李少卓:监督、数据管理。

    数据可用性

    支持本研究发现的数据无法共享,因为将进一步进行研究。

    CRediT作者贡献声明

    张传:写作 – 原始草稿、验证、方法论、数据管理、概念化。潘永泰:写作 – 审稿与编辑、验证、资金获取、概念化。曹星健:软件、方法论、研究。尹芳茜:软件、方法论。张永强:验证、正式分析。李少卓:监督、数据管理。

    利益冲突声明

    作者声明他们没有已知的竞争性财务利益或个人关系可能影响本文报告的工作。

    致谢

    作者感谢国家自然科学基金(项目编号52074308)的支持。
    相关新闻
    生物通微信公众号
    微信
    新浪微博
    • 搜索
    • 国际
    • 国内
    • 人物
    • 产业
    • 热点
    • 科普

    知名企业招聘

    热点排行

      今日动态 | 人才市场 | 新技术专栏 | 中国科学人 | 云展台 | BioHot | 云讲堂直播 | 会展中心 | 特价专栏 | 技术快讯 | 免费试用

      版权所有 生物通

      Copyright© eBiotrade.com, All Rights Reserved

      联系信箱:

      粤ICP备09063491号