💻 computer science

Phase-Field Peridynamics

本文提出了一种新型相场非局部动力学方法，通过引入键相场参数连续退化能量贡献并采用独立的运动学退化函数，解决了传统键删除法导致的数值不稳定问题，同时保持了非局部变形梯度计算的精度与热力学一致性。

原作者： Kai Partmann, Christian Wieners, Michael Ortiz, Kerstin Weinberg

发布于 2026-04-13

📖 1 分钟阅读☕ 轻松阅读

原作者： Kai Partmann, Christian Wieners, Michael Ortiz, Kerstin Weinberg

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

这篇论文介绍了一种名为**“相场近场动力学”（Phase-Field Peridynamics, PFPD）**的新方法，旨在更稳定、更准确地模拟材料是如何断裂的。

为了让你更容易理解，我们可以把这篇论文的核心思想想象成**“修补一张破网”**的故事。

1. 背景：材料为什么会断？

想象你手里拿着一张巨大的渔网（这就好比一块金属或混凝土材料）。

传统方法（有限元法）： 就像用剪刀剪断网线。一旦剪断，那个地方就彻底空了。如果裂缝蔓延，你需要不断重新画网格，非常麻烦，而且如果剪得太快，计算容易出错。
近场动力学（Peridynamics）： 这是一种更聪明的方法。它不把材料看作连续的固体，而是看作无数个小点，每个点都通过“弹簧”（称为键，Bond）和周围的点相连。如果某个弹簧拉得太长，它就断了。这种方法天生就能处理断裂，不需要重新画网格。

2. 问题：现有的“剪断”方法有个大毛病

在现有的近场动力学中，当弹簧（键）拉断时，计算机的做法是直接把它从计算中删掉（就像把断掉的绳子直接扔掉）。

比喻： 想象你在玩一个由无数根绳子连接的气球网络。如果为了模拟气球破裂，你直接剪断并扔掉绳子，网络的结构就会变得支离破碎。
后果： 当裂缝变大，周围被“扔掉”的绳子太多时，剩下的绳子无法再准确支撑起整个结构的形状计算。这就像搭积木，如果你把底部的几块突然拿走，上面的积木就会乱套，导致计算崩溃（数值不稳定）。特别是在材料受剪切力（像剪刀一样错动）时，这种“直接删除”的方法经常会让电脑算不出来。

3. 解决方案：不要“剪断”，要“软化”

这篇论文提出的新方法（PFPD）做了一个巧妙的改变：不再直接剪断绳子，而是让绳子慢慢变软。

核心创新 1：能量降解（相场参数）
想象那根绳子不是突然消失，而是慢慢变成了“果冻”。它还在连接着两个点，但它传递力量的能力越来越弱，直到最后几乎为零。
- 好处： 绳子还在，网络的结构完整性就还在，计算机不会“晕头转向”。
核心创新 2：运动学保护（Kinematic Degradation）
这是论文最精彩的地方。虽然绳子变软了（能量没了），但为了保持网络形状的准确性，作者设计了一个特殊的“保护机制”。
- 比喻： 想象这根绳子虽然变成了果冻（不再传递拉力），但它依然是一根**“隐形的骨架”**，帮助维持周围点的相对位置。只有当绳子彻底烂掉（完全断裂）时，这个骨架才慢慢消失。
- 结果： 即使材料已经裂开了，计算机依然能准确地算出材料是怎么变形的，不会因为裂缝太大而算错。

4. 数学上的“魔法”：归一化常数

为了让这种“变软”的过程符合物理定律（即断裂消耗的能量要符合格里菲斯断裂理论），作者推导出了一个数学公式。

比喻： 就像你在调一杯特制鸡尾酒，你需要知道加多少“软度剂”才能让味道（能量）刚刚好。作者通过复杂的数学推导，找到了这个完美的配方（归一化常数），确保无论用什么样的“绳子”（核函数），计算出的断裂能量都是准确的。

5. 实验证明：它真的好用

作者用这个新方法做了几个著名的测试：

拉断测试（Mode I）： 像拉橡皮筋一样把材料拉断。结果：裂缝直直地走，非常准。
剪切测试（Mode II）： 像用剪刀剪东西。
- 旧方法： 算到一半就崩溃了，因为绳子被删太多，结构乱了。
- 新方法： 稳稳地算出了裂缝弯曲的路径，非常稳定。
冲击实验（Kalthoff-Winkler）： 模拟子弹打钢板。
- 结果： 新方法算出的裂缝角度（约 70 度）和真实实验观察到的完全一致，而且裂缝扩展的速度曲线非常平滑，没有那种忽高忽低的“抖动”。

总结

这篇论文就像给材料断裂模拟装上了一个**“智能减震器”**。

以前： 遇到裂缝，直接切断连接，导致结构不稳，计算容易崩。
现在： 遇到裂缝，先让连接“软化”并保留“隐形骨架”，既模拟了断裂，又保证了计算过程的稳定。

这使得科学家和工程师能够更放心地模拟复杂的动态断裂场景（如汽车碰撞、飞机零件疲劳、地震中的建筑倒塌），而不用担心电脑程序会因为算不出裂缝而崩溃。

这是一份关于论文《Phase-Field Peridynamics》（相场近场动力学）的详细技术总结。该论文提出了一种将相场断裂理论与近场动力学（Peridynamics, PD）相结合的新方法，旨在解决传统近场动力学在断裂模拟中遇到的数值稳定性问题。

1. 研究背景与问题 (Problem)

断裂模拟的挑战：预测裂纹的萌生和扩展是计算力学的核心难题。传统的有限元法（FEM）需要显式追踪裂纹，导致网格重划分成本高昂；扩展有限元法（XFEM）和相场法（Phase-Field）虽然有所改进，但相场法通常需要极细的网格来解析扩散损伤带，且基于经典连续介质力学（偏微分方程），在裂纹处导数未定义。
近场动力学（PD）的优势与局限：PD 通过积分微分方程描述动量平衡，无需空间导数，天然适合处理不连续性（如裂纹）。
- 对应关系公式（Correspondence Formulation）：允许直接使用经典本构模型，消除了表面效应，是精度最高的 PD 公式之一。
- 键关联方法（Bond-Associated Approach）：通过计算键相关的变形梯度，消除了零能模式（Zero-energy modes），进一步提高了精度。
核心痛点：传统的 PD 断裂模拟采用**不可逆的键删除（Bond Deletion）机制。当键的变形超过临界值时，键被永久移除。在基于对应关系的公式中，随着裂纹扩展，大量键被删除会导致计算变形梯度所需的形状张量（Shape Tensor）**病态甚至奇异，进而破坏运动学近似，引发数值不稳定。这在剪切断裂（Mode II）等复杂工况下尤为严重。

2. 方法论 (Methodology)

论文提出了一种**相场近场动力学（PFPD）**方法，将相场断裂理论嵌入到键关联对应关系框架中，核心思想是用“连续退化”替代“瞬时删除”。

2.1 核心机制

键相场参数 ( $s$ )：引入一个键相关的相场变量 $s \in [0, 1]$ 来描述损伤状态（0 为无损，1 为完全损伤）。
双重退化函数：
1. 能量退化函数 $g(s)$ ：连续降低键的能量贡献（应力响应）。通常采用二次函数 $g(s) = (1-s)^2$ ，确保应力随损伤平滑下降至零。
2. 运动学退化函数 $h(s)$ ：这是该方法的创新点。为了保持形状张量和变形梯度近似的数值稳定性，引入一个独立的运动学退化函数。
  - 当损伤 $s$ 低于阈值 $s_c$ 时， $h(s) = 1$ （保持运动学贡献完整）。
  - 当 $s > s_c$ 时， $h(s)$ 逐渐减小。
  - 作用：即使键完全断裂（ $s=1$ ），其运动学贡献也不会立即消失，从而防止形状张量奇异，保证了大变形和严重损伤下的数值稳定性。

2.2 能量泛函与演化

总能量泛函：包含动能、退化后的应变能密度以及基于 Griffith 理论的裂纹表面能。
损伤演化：基于变分原理最小化总能量。
- 定义裂纹驱动力 $Y$ （基于拉伸应变能密度或最大主应力）。
- 引入历史变量 $Y$ 确保损伤不可逆。
- 相场演化方程： $s = \min(1, Y / (Y + Y_c))$ ，其中 $Y_c$ 是临界裂纹驱动力。

2.3 归一化常数的解析推导

为了确保模型符合 Griffith 断裂理论（即单位裂纹面积的能量释放率等于 $G_c$ ），必须确定归一化常数 $c_0$ 。
作者针对一般球形核函数，在无限域平面裂纹假设下，通过解析推导得出了 $c_0$ 的精确闭式解。
结果： $c_0$ 被表示为两个一维积分的比值：
$c_0 = \frac{\int_0^1 \omega_s(\rho\delta) \rho^3 d\rho}{2 \int_0^1 \omega_s(\rho\delta) \rho^2 d\rho}$
其中 $\omega_s$ 是标量核函数， $\delta$ 是视界半径。这使得该常数可以针对任意球形核函数精确计算，无需数值拟合。

3. 主要贡献 (Key Contributions)

提出 PFPD 框架：成功将相场法与键关联对应关系 PD 结合，解决了传统键删除法在对应关系公式中导致的数值不稳定问题。
运动学退化策略：创新性地引入了独立的运动学退化函数 $h(s)$ ，在保持能量耗散物理意义的同时，维持了形状张量的良态性，显著提高了剪切断裂等复杂工况下的稳定性。
归一化常数的通用解析解：推导出了适用于任意球形核函数的归一化常数 $c_0$ 的精确解析表达式，消除了以往依赖特定核函数或数值校准的局限性。
算法实现：详细描述了离散化方案，并对比了传统 BAQP 模型（键删除）与 PFPD 模型（相场退化）的算法流程，证明了 PFPD 在保持热力学一致性的同时具有更好的鲁棒性。

4. 数值结果 (Results)

论文通过多个基准算例验证了 PFPD 模型的有效性：

I 型和 II 型裂纹扩展：
- 在 I 型（拉伸）和 II 型（剪切）裂纹扩展中，PFPD 均能准确预测裂纹路径。
- 关键对比：在 II 型剪切加载下，传统的 BAQP 模型（键删除）因形状张量退化而失稳，导致计算失败；而 PFPD 模型保持稳定，成功模拟了裂纹偏转。
边界拉伸测试 (BTT) 与归一化验证：
- 使用不同的视界比（ $m_\delta$ ）和核函数（三次 B 样条、线性核）进行模拟。
- 结果显示，无论核函数和视界比如何变化，裂纹分叉的角度和位置高度一致，验证了推导出的归一化常数 $c_0$ 的正确性和通用性。
- 裂纹尖端速度平滑，未超过理论极限（半瑞利波速），且比传统键基模型更平滑。
Kalthoff-Winkler 冲击实验：
- 模拟了钢板的冲击断裂实验。
- PFPD 预测的裂纹扩展角度约为 $70^\circ$ ，与实验观测值高度吻合。
- 相比普通状态基（OSB）模型，PFPD 的裂纹尖端速度曲线更平滑，且裂纹萌生时间更准确，体现了连续损伤演化的优势。

5. 意义与结论 (Significance & Conclusion)

理论突破：该工作证明了相场法与高精度近场动力学（键关联对应关系）可以无缝融合，既保留了 PD 处理不连续性的天然优势，又利用了相场法处理复杂裂纹拓扑（如分叉、合并）的能力。
解决稳定性难题：通过“连续退化”替代“键删除”，特别是引入运动学退化函数，彻底解决了对应关系公式在断裂模拟中因形状张量奇异而导致的数值不稳定问题。
工程应用价值：PFPD 模型在动态断裂、剪切断裂及冲击载荷下表现出卓越的鲁棒性和准确性，为复杂工程结构（如航空航天、汽车安全）的断裂预测提供了一个可靠、热力学一致且无需显式追踪裂纹的计算工具。
未来展望：该方法为开发更复杂的材料断裂模型（如各向异性、多物理场耦合）奠定了坚实基础，且其解析归一化常数可直接应用于不同的核函数选择中。

总结：这篇论文通过引入相场概念和创新的运动学退化机制，成功克服了传统近场动力学在断裂模拟中的稳定性瓶颈，提供了一种高精度、高稳定性且物理意义明确的动态断裂模拟新范式。