连续体、离散与混合断裂建模方法的比较研究

《Engineering Fracture Mechanics》：Continuum, discrete, and hybrid approaches to fracture modeling: A comparative study

　　

对连续体、离散与混合断裂建模方法进行了系统比较。计算断裂力学的发展沿着多条部分重叠的路径展开，其中一种广泛使用的分类视角是根据裂纹在数值模型中的表示方式来划分方法。在几何表示中，裂纹被视为显式实体，因此裂纹扩展需要对几何和/或离散化进行更新；而在非几何表示中，模型不显式包含裂纹，而是通过本构描述或强局部化等机制隐式捕捉断裂。本文针对准脆性材料中三种代表性断裂建模方法——变分相场法（Variational Phase-Field, VPF）、晶格单元法（Lattice Element Method, LEM）和有限-离散单元法（Finite-Discrete Element Method, FDEM）——开展了直接的基准比较。研究的创新之处在于采用相同的I型（mode-I）和复合型（mixed-mode）断裂问题集，对比了这三种不同建模选择所预测的裂纹轨迹和全局荷载-位移响应。据研究人员所知，本研究提供了这三个断裂建模框架之间独特的系统性比较，能够通过阐明它们各自的能力、局限性和计算影响，为针对特定断裂问题选择合适的方法提供支持。结果表明，比较研究可作为面向应用的方法选择指南。

论文解读

一、研究背景与问题提出

在工程材料的断裂预测研究中，计算断裂力学已成为处理复杂几何、材料非线性及混合边界条件问题的关键手段，因为这类问题通常无法通过解析解解决。然而，尽管连续体、离散体和混合断裂建模方法均得到了广泛发展，基于统一基准的系统性跨范式比较仍然十分匮乏。现有研究往往聚焦于单一方法的发展或其特定问题的应用，导致难以评估不同断裂描述方式如何从根本上影响裂纹演化和结构响应的预测结果。因此，研究人员亟需填补这一空白，通过直接的基准对比揭示各类方法的特性差异。

二、研究内容与结论

研究人员选取了变分相场法（Variational Phase-Field, VPF）、晶格单元法（Lattice Element Method, LEM）和有限-离散单元法（Finite-Discrete Element Method, FDEM）作为三种范式的代表方法，在统一的I型（mode-I）和复合型（mixed-mode）断裂基准问题框架下进行对比。研究发现，这三种方法在预测裂纹轨迹和全局荷载-位移响应方面表现出显著差异，从而揭示了不同断裂描述方式对模拟结果的直接影响。这项研究的重大意义在于：它首次提供了三个断裂建模框架之间的系统性比较，能够帮助研究者根据具体断裂问题的特点选择最合适的建模方法，同时明确了各方法的优势、局限性及其计算代价。

三、关键技术方法

研究人员采用了三种代表性的断裂建模方法：

1. 1.
   变分相场法（VPF）：采用Miehe等人提出的热力学一致框架，引入辅助标量相场变量d∈[0,1]区分完整材料（d=0）与完全断裂状态（d=1）。该方法将尖锐裂纹表面Γ替换为由内部长度尺度控制的弥散裂纹带，裂纹表面密度函数定义为γ(d,?d)=d2/(2l)+l/2|?d|2，其中l为长度尺度参数，G_c为材料临界能量释放率。
2. 2.
   晶格单元法（LEM）：将材料离散化为由一维结构单元构成的网络，断裂通过单个晶格单元的失效或渐进刚度退化来模拟。该方法特别适用于表征分布式微裂纹、断裂过程区（Fracture Process Zone, FPZ）的形成以及损伤向宏观裂纹的局部化演化。
3. 3.
   有限-离散单元法（FDEM）：由Munjiza等人提出的混合公式，将非线性断裂力学概念嵌入材料的连续体描述中，使材料能从完整状态过渡到完全非连续状态，显式模拟裂纹萌生、扩展、接触和碎裂全过程。

四、研究结果

2.1 VPF（变分相场法）

VPF将Griffith断裂理论重新表述为变分能量最小化问题，裂纹萌生与扩展源于体弹性能与断裂表面能的竞争。通过正则化裂纹泛函，尖锐裂纹表面被近似为弥散裂纹场参数，断裂演化遵循总能量泛函的最小化，从而在连续体框架内自然捕捉裂纹成核、扩展及复杂裂纹拓扑。

2.2 LEM（晶格单元法）

LEM将材料离散化为固定节点连接的一维结构单元网络。断裂在单个晶格单元层面建模，当预设准则达到时单元发生破坏。可通过单元移除模拟脆性断裂响应，也可通过渐进刚度退化模拟呈现峰后软化的渐进损伤过程。

2.3 FDEM（有限-离散单元法）

FDEM专为动态模拟可变形相互作用体而设计。它将非线性断裂力学概念嵌入材料的连续体描述中，使材料从完整状态向完全非连续状态的转变过程中，能够显式模拟断裂萌生、裂纹扩展及后续碎裂。

五、讨论与结论

讨论总结：尽管每种方法家族内部存在大量变体和扩展形式，但本研究选取的三种代表性方法——VPF、LEM和FDEM——分别代表了连续体、离散体和混合范式，各自满足广泛的断裂问题适用性、方法的成熟度和公认度、以及适合系统基准对比的既定公式和计算实现等选择标准。对于VPF，采用Miehe等人的公式因其清晰的变分结构、成熟的数值鲁棒性以及在基准和比较研究中的频繁使用；对于LEM和FDEM，则分别采用离散晶格单元法和Munjiza提出的原始FDEM框架。

研究结论翻译：本研究通过对代表性连续体、非连续体和混合方法（即VPF、LEM和FDEM）在统一基准框架下的直接比较，填补了现有研究中系统性基准对比的空白。研究的创新之处在于采用相同的I型和复合型断裂问题集考察这三种截然不同的建模选择，并比较其预测的裂纹轨迹和全局荷载-位移响应。据研究人员所知，本研究提供了这三个断裂建模框架之间独特的系统性比较，能够通过阐明它们各自的能力、局限性和计算含义，为针对特定断裂问题选择合适的方法提供支持。