Endowing variational phase-field fracture models with custom strength criteria

本文提出了一种新型变分相场断裂框架，该框架通过引入状态相关耗散势，在保持变分结构的同时，实现了对任意弹性域和自定义强度准则的整合，从而允许在多轴应力状态下对弹性退化和裂纹成核进行独立控制。

要点

想象一下，你正试图预测一块玻璃或混凝土在受到挤压、拉伸或扭转时会如何产生裂纹。长期以来，科学家们一直使用一种被称为“相场”（phase-field）法的巧妙数学工具来模拟这一过程。你可以将这种方法想象成一张高科技的裂纹天气图：它不是画出一条尖锐、锯齿状的线来表示裂纹出现的位置，而是涂抹出一个模糊、柔和的区域，这个区域从“固体”逐渐过渡到“破碎”。

然而，旧的地图存在一个问题。它们就像是一件“一码通”的成衣：假设无论你是拉伸还是挤压材料，材料的破坏方式都是一样的。但实际上，材料是非常挑剔的。例如，混凝土讨厌被拉开，但在受压时却非常坚韧。旧模型无法在不破坏使其运作的数学规则的前提下，轻松区分这些不同的应力“个性”。

新思路：一件定制西装

本文作者提出了一种构建这些模型的新方法。他们称之为“赋予变分相场断裂模型自定义强度准则”。用通俗的话说，他们想出了如何为这些裂纹预测模型打造一件量身定制的西装，以适应任何特定材料的规则。

以下是他们如何实现的，使用了简单的类比：

两部分系统：夹克与护盾

想象一下，某种材料穿着两层衣服：

1. 夹克（自由能）： 这一层代表材料的刚度。随着材料受损（就像夹克上出现了洞），它会变得越来越弱、越来越不硬。在旧模型中，夹克和破裂的规则是粘在一起的。如果你改变了夹克，就会不小心改变了破裂的规则。
2. 护盾（耗散势）： 这一层代表材料的强度或其“破坏点”。它决定了需要多少力量才能开始产生裂纹。

创新之处：
作者意识到，可以让护盾根据你推或拉材料的方式来改变形状，而不会干扰到夹克。

• 旧方法： 如果你想让材料在受压时比受拉时更强，你必须重写整个夹克的数学公式。这既混乱又容易破坏维持模拟稳定的“变分结构”（内部逻辑）。
• 新方法： 他们让护盾变成了“状态相关型”。这意味着，根据受力的方向，护盾可以呈现为圆形、椭圆形或奇形怪状的团块。
  • 如果你拉伸材料，护盾可能很小（容易破碎）。
  • 如果你挤压材料，护盾可能很大（难以破碎）。
  • 至关重要的是，夹克（刚度）保持完全不变。两者现在是独立的。

“弹性域”地图

论文多次提到“弹性域”。想象一张安全区地图。只要作用在材料上的力保持在这个区域内，材料就是安全的，不会产生裂纹。

• 在旧模型中，这个安全区总是完美的对称圆（或半圆）。
• 在新模型中，作者可以绘制出任何形状的安全区。
  • 他们可以制作一个双椭圆（像花生形状一样），以处理拉伸与挤压之间的不同极限。
  • 他们可以制作一个 Drucker-Prager 圆锥体（像冰淇淋甜筒一样），来模拟在压力下表现不同的岩石或土壤。
  • 他们可以制作一个 Huber 形状，允许材料在受压时不产生裂纹（非互穿性），但在被拉开时仍容易破碎。

这为什么重要（根据论文所述）

作者通过几种不同的“配方”（模型 M1 到 M5）测试了他们的新方法。他们模拟了一个在不同角度下受到拉力和压力作用的圆盘材料。

1. 灵活性： 他们展示了如何创建一个模型，使材料在被拉开时容易破碎，而在被挤压时非常坚固，反之亦然，同时保持数学逻辑的简洁与稳定。
2. 独立性： 他们证明了你可以分别调节“刚度”（弯曲程度）和“强度”（何时破碎）。在此之前，改变其中一个往往会迫使你改变另一个。
3. 准确性： 模拟结果显示，裂纹恰好在自定义“安全区”地图所指示的位置开始，能够匹配各种复杂的加载条件（如同时进行扭转和挤压）。

核心结论

本文并不声称能立即治愈疾病或建造新桥梁。相反，它提供了一个更灵活的数学工具箱。它允许科学家构建能够尊重不同材料（如混凝土、岩石或生物组织）特定且古怪规则的计算机模拟，而不会破坏使模拟可靠的物理基本定律。这就像是从使用通用的预制地图升级到了一个可以为任何地形绘制自定义路线的 GPS。

技术摘要

技术摘要：赋予变分相场断裂模型自定义强度准则

问题陈述
尽管相场断裂建模已成为捕捉复杂裂纹模式的成熟理论，但仍存在一个显著局限：无法任意规定强度准则，且难以在多轴应力状态下一致地描述相应的弹性域。现有方法大致分为三类：(a) 通过分裂自由能，这种方法保留了变分结构，但在定义任意弹性域方面缺乏灵活性；(b) 直接修改相场驱动力或阻力，这种方法具有通用性，但牺牲了变分结构；(c) 与塑性或非线性应变耦合，这种方法保持了变分一致性，但引入了额外的内变量和复杂的公式化过程。至关重要的是，在大多数现有方法中，弹性退化行为与强度准则是强耦合的，这导致无法独立调节刚度退化与材料强度。

方法论
作者提出了一种在广义标准材料（GSM）能量框架内构建的根本不同的策略。其核心创新在于引入了一个状态相关耗散势。不同于耗散势为常数或仅取决于相场变量的标准模型，该方法允许耗散势取决于材料的当前状态（具体为应变或应力不变量）。

关键方法组成部分包括：

• 机制分离： 自由能密度 ($\psi$) 仅负责编码弹性退化（刚度损失），而状态相关耗散势 ($\phi$) 则负责编码强度准则和弹性域的形状。这种解耦使得这两个截然不同的材料属性可以被独立控制。
• 断裂能依赖性： 断裂能 $G_c$ 被视为加载方向（通过体积-偏应变或应力平面中的角度 $\vartheta$ 表示）的函数。这种依赖性通过耗电势中的断裂函数 $G_f(\varepsilon, \alpha)$ 被纳入其中。
• 变分一致性： 通过从全局能量平衡和稳定性条件（能量公式）中推导演化方程，所提出的模型保留了问题的变分结构，从而保留了相关的理论和计算优势。
• 数值实现： 采用交替求解方案，交替求解平衡问题和损伤问题。在迭代过程中，状态相关的断裂函数在上一收敛载荷步被冻结，以近似状态相关的耗散能。

主要贡献与材料模型
论文通过构建并分析五个特定的材料模型（M1–M5）展示了该框架的灵活性，这些模型结合了不同的退化律和强度准则：

1. M1 (标准 AT1)： 全弹性退化 (FED) 结合标准的各向同性强度准则。
2. M2 (双椭圆)： FED 结合可调的双椭圆强度准则，允许任意的体积和偏应力极限。
3. M3 (Drucker–Prager)： FED 结合 Drucker–Prager 准则，适用于压力敏感型材料。
4. M4 (带非互穿性的 Drucker–Prager)： 部分弹性退化 (PED)，其中压缩体模系数不发生退化，但仍能恢复 Drucker–Prager 强度准则。
5. M5 (带非互穿性的 Huber)： PED 结合 Huber 型准则，允许独立调节拉伸-体积极限与偏应力极限，同时确保压缩下的非互穿性。

结果
论文展示了在各种加载方向 ($\vartheta$) 下，双轴圆盘问题的解析和数值结果：

• 均匀响应： 解析解证实，这些模型成功重现了不同的弹性域演化。例如，M3 和 M4 具有相同的弹性域形状，但由于是否存在体模系数退化，表现出不同的力学响应。M5 展示了非对称的退化行为，区分了拉伸体积应变和压缩体积应变。
• 裂纹成核： 多轴加载下的裂纹成核数值模拟表明，这些模型能够准确预测与预设强度面相对应的成核点。结果说明，所提出的框架可以捕捉到与每个模型定义的特定强度准则相一致的混合模式成核行为（例如，各向同性膨胀、单轴拉伸、纯剪切、单轴压缩）。
• 特征独立性： 结果验证了弹性退化（由自由能控制）与强度准则（由耗散势控制）可以被独立调节，这是以往方法难以实现的特性。

意义与主张
作者声称，这项工作通过提供一个灵活且具有变分一致性的框架，为将任意强度准则引入相场断裂模型铺平了道路。其主要意义在于：

1. 在引入复杂的、状态相关的强度准则的同时，保留了问题的变分结构。
2. 将刚度退化与强度解耦，从而实现更真实的、可独立控制的材料建模。
3. 系统地构建复杂材料（如准脆性、摩擦性或压力敏感型材料）的模型，而无需借助非变分源项或过多的内变量。

论文结论较为谦逊，指出虽然该框架能很好地处理成核和稳定演化，但捕捉不稳定的裂纹扩展事件将需要未来集成时间正则化策略（如惯性或粘性项）。当前的工作为这些扩展奠定了理论基础。