IMDEA 材料研究所提出双梯度校正方法，解决相场断裂模拟中的裂纹面积高估问题
来自 IMDEA 材料研究所（IMDEA Materials）和马德里技术大学（Universidad Politécnica de Madrid）的研究团队开发了一种新的数学模型，旨在修正材料破裂数值模拟中长期存在的裂纹面积高估现象。这项研究成果发表于科学期刊《计算力学》（Computational Mechanics），为在不增加计算成本的前提下提高结构故障预测的准确性提供了新途径。
背景：应变局部化导致的数值偏差
基于相场（Phase-field）的计算模型是现代工程中预测裂纹传播及材料在外部载荷下失效行为的重要工具。然而，传统的数值实现方案存在一个广为人知的缺陷，即“应变局部化”（strain localization）。这种数值伪影会导致损伤变量在模拟区域内出现人工扩散，扭曲计算结果，使得软件预测出的破裂区域远大于实际物理尺寸。
DGCM 方法原理：利用断裂能量等分原理
为了纠正这一偏差，研究人员提出了“双梯度校正方法”（Double Gradient Correction Method, DGCM）。该方法的核心在于利用一个基本的物理性质——断裂能量的等分原理。研究作者米格尔·卡斯蒂隆博士（Dr. Miguel Castillón）、哈维尔·塞古拉多教授（Prof. Javier Segurado）和伊格纳西奥·罗梅罗（Ignacio Romero）证明，尽管数值错误会严重影响直接的损伤变量，但它们对与损伤变量空间梯度相关的能量影响微乎其微。
通过重复计算与损伤变量空间梯度相关的贡献项，DGCM 能够获得更准确的断裂表面估计。验证测试表明，与传统模拟方法相比，这种方法显著减少了偏差。
应用前景与开源代码发布
卡斯蒂隆指出：“这种新的数学方法使我们能够更可靠地分析接近其故障极限的材料的行为，在减少数值错误的同时不增加计算复杂性，并且可以应用于三维模型。”
目前，实施该方法所需的代码已在开源许可证下发布，以促进其集成到标准的工业模拟工具中。该存储库托管于 GitHub，并在 Zenodo 上进行了永久存档。相关数字模拟使用了名为 "PhaseFieldX" 的开源库进行，该库基于 FEniCSx 有限元素框架构建。所有研究结果均可通过提供的材料完全重现。
论文信息：
- 标题：A correction method for crack area overestimation in phase-field fracture
- 作者：Miguel Castillón et al.
- 期刊：Computational Mechanics (2026)
- DOI：10.1007/s00466-026-02834-2
