On modeling fracture of soft polymers

该论文提出了一种基于临界应力功的损伤模型，通过结合能量驱动的损伤起始准则与广义梯度损伤框架，成功统一描述了多种软聚合物在不同加载条件和几何构型下的断裂行为，并实现了从宏观参数到微观链解离能的可预测性。

要点

这篇论文主要解决了一个困扰材料科学界的难题：如何准确预测那些“软绵绵”的材料（如橡胶、凝胶、生物组织）在什么情况下会断裂？

想象一下，你手里拿着一块口香糖（软聚合物）。如果你慢慢拉它，它会被拉得很长很长，甚至像融化的糖一样流动；但如果你猛地一扯，它可能“啪”地一下就断了。以前的模型很难统一解释这种“慢拉不断，快拉就断”的现象，也很难预测它到底会在哪里断、断成什么形状。

这篇论文提出了一套新的“通用法则”，就像给这些软材料装上了一个智能的“断裂预警系统”。

以下是用通俗语言和比喻对这篇论文核心内容的解读：

1. 核心难题：以前的尺子不好使了

• 过去的困境：以前科学家测量材料“抗断能力”（断裂韧性）时，用的尺子（公式）是固定的。但这就像用一把尺子去量水和量石头一样，对于软材料来说，这把尺子会“失灵”。
  • 如果你拉得快，尺子读数是一个值；拉得慢，读数又变了。
  • 如果你把材料切成三角形，读数是一个值；切成圆形，读数又变了。
  • 这导致科学家无法用一个统一的公式来预测所有情况下的断裂。
• 比喻：这就好比你想用“重量”来衡量一个人的“力气”。举重运动员举得重，但让他去跑马拉松，力气再大也跑不动。以前的模型就像只看重“举重”，忽略了“跑步”（加载速度）和“姿势”（材料形状）的影响。

2. 新发现：一把神奇的“能量尺” ($W_{cr}$)

• 新方案：作者提出，不管你是快拉还是慢拉，不管材料是什么形状，只要计算**“把材料彻底弄坏所需的总能量密度”（即单位体积内消耗的能量），这个数值对于某一类软材料来说，几乎是恒定不变**的。
• 比喻：想象你要把一根橡皮筋彻底扯断。
  • 慢慢拉，橡皮筋会发热、变形，消耗很多能量。
  • 快速拉，橡皮筋来不及变形，直接崩断，但也消耗了能量。
  • 作者发现，虽然过程不同，但**“彻底扯断它所需的总能量”**是一个固定的“门槛值”。只要输入的能量超过了这个门槛，材料就会坏。这个门槛值就是论文中的 $W_{cr}$。

3. 新模型：给材料装上“智能损伤计数器”

为了预测断裂，作者开发了一个新的数学模型，它像是一个多机制的“损伤计数器”：

• 双保险机制：软材料内部有两种“骨架”在起作用：
  1. 分子间的摩擦（像拥挤的人群）：分子之间互相推挤、摩擦，产生阻力。
  2. 化学键的拉伸（像弹簧网）：分子之间通过化学键连接，像弹簧一样被拉长。
  • 模型认为，只有当这两种“骨架”都彻底损坏（弹簧断了、人群散架了），材料才会真正断裂。
• 损伤的“蔓延”：
  • 当局部能量超过门槛，损伤（Damage）就开始产生。
  • 这个损伤不会瞬间传遍全身，而是像**“病毒”或“裂纹”**一样，从一点开始慢慢扩散。
  • 模型还考虑了“梯度”，即损伤边缘是模糊的，而不是像刀切一样锋利，这更符合真实物理现象。

4. 惊人的预测能力：从“喇叭口”到“抛物线”

这个模型最厉害的地方在于，它能预测出材料断裂时的形状，而且这种形状会随着拉伸速度变化：

• 慢速拉伸（像拉口香糖）：
  • 材料有足够时间“思考”和重组。
  • 现象：裂纹尖端会变钝（Blunting），然后像吹喇叭一样，在钝头前面又形成一个新的小缺口。
  • 比喻：就像你慢慢撕开一张湿纸巾，它会先变宽，然后像**“喇叭花”**一样裂开。
• 快速拉伸（像扯断橡皮筋）：
  • 材料来不及重组，直接崩断。
  • 现象：裂纹尖端保持尖锐，直接裂开。
  • 比喻：就像你快速撕开干纸，裂口是**“抛物线”**形状的，干脆利落。
• 模型成果：以前的模型很难同时预测这两种截然不同的形状，而这个新模型通过调整参数，完美复现了从“喇叭口”到“抛物线”的所有变化。

5. 实际应用：从微观到宏观的“翻译官”

• 微观洞察：这个模型不仅能算宏观的断裂，还能反推微观的分子行为。
  • 作者通过宏观的断裂数据，成功估算出了微观层面上**“化学键断裂所需的能量”**。
  • 比喻：就像通过观察大楼倒塌的方式，反推出大楼里每一根钢筋的强度。
• 验证广泛：他们用这个模型测试了多种材料：
  • PBS（一种像果冻一样的粘性聚合物）：能自愈合，能流动。
  • EPDM、SBR、NR（各种橡胶）：常见的轮胎、密封圈材料。
  • EPS25（一种特殊的玻璃态聚合物）。
  • 结果发现，无论材料是像橡胶一样有弹性，还是像蜂蜜一样粘稠，这个模型都能准确预测。

总结

这篇论文就像给软材料领域发明了一套**“通用语言”**。

以前，科学家面对不同的软材料、不同的拉伸速度，需要换不同的公式，而且经常算不准。现在，他们只需要知道**“破坏这个材料需要多少总能量”，就能用一个统一的模型，精准预测它什么时候断、在哪里断、以及断成什么形状**。

这对于设计更安全的人造血管、柔性机器人、防弹衣、甚至自愈合的电子设备都至关重要。它告诉我们，无论材料是快是慢、是软是硬，只要掌握了那个核心的“能量门槛”，我们就能掌控它们的命运。

技术摘要

软聚合物断裂建模研究技术总结

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

软聚合物（包括弹性体、粘弹性聚合物、聚合物凝胶及生物组织）在自然界和工程应用中无处不在。然而，针对这类材料的断裂建模面临以下核心挑战：

• 速率依赖性难题：软聚合物的力学行为从速率无关（如弹性体）到强速率相关（如动态交联聚合物）变化巨大。现有的断裂韧性指标（如能量释放率或裂纹尖端应力强度因子）对于速率相关材料缺乏唯一性，其数值随加载路径和试件几何形状的不同而显著变化。
• 模型局限性：现有的断裂模型大多局限于特定的预裂纹几何形状，或仅适用于弹性体（基于梯度损伤框架）。缺乏一个通用的、统一的模型来预测从粘性聚合物到弹性体的损伤萌生、扩展及完全断裂，且该模型需独立于加载历史和试件几何。
• 物理现象捕捉不足：现有变分相场方法难以捕捉实验中观察到的特定失效面（failure surfaces），且难以同时预测脆性与延性响应以及不同的断裂轮廓（如粘性聚合物的“喇叭形”轮廓与弹性体的“抛物线形”轮廓）。

2. 方法论 (Methodology)

作者提出并验证了一种基于**广义梯度损伤框架（Generalized Gradient-Damage Framework）和虚功率原理（Method of Virtual Power）**的断裂模型。

2.1 核心理论创新

• 临界应力功 ($W_{cr}$)：提出将“临界应力功”作为衡量特定类别软聚合物断裂阻力的唯一物理参数。$W_{cr}$ 定义为在足够快的加载速率下，导致任意材料点完全失效（无法承受任何机械应力）所需的单位参考体积总能量密度。
  • 假设：对于给定的加载模式（如拉伸），$W_{cr}$ 是加载速率和试件几何形状的独立常数。
  • 优势：$W_{cr}$ 的单位是 $J/m^3$（能量密度），不同于传统断裂韧性的 $J/m^2$（单位面积能量），使其具有几何无关性。
• 损伤萌生准则：将 $W_{cr}$ 中的能量贡献部分（$\psi^+_{cr}$）作为损伤萌生准则。该准则显式地考虑了材料强度，能够捕捉失效面。
• 多机制梯度损伤框架：
  • 采用 Kröner-Lee 乘法分解和并联多机制方法，将变形分为分子间阻力（机制 A，包含弹塑性）和网络阻力（机制 B，包含动态交联链的拉伸）。
  • 引入损伤变量 $d$ 和梯度项，通过虚功率原理推导宏观和微观力平衡方程。
  • 损伤演化方程结合了历史函数 $H(t)$，确保损伤仅在驱动能量超过临界值时发生，并考虑了粘性效应导致的网络重构。

2.2 数值实现

• 模型在有限元软件 ABAQUS 中实现。
• 开发了用户自定义单元（VUEL）用于处理损伤变量 $d$ 的演化，以及用户自定义材料子程序（VUMAT）用于处理力学响应。
• 采用显式积分格式以适应短时间尺度的模拟，并通过全局变量在 VUEL 和 VUMAT 之间传递数据。

3. 关键贡献 (Key Contributions)

1. 提出并验证了 $W_{cr}$ 作为断裂阻力度量：通过实验证明，对于特定类别的软聚合物，$W_{cr}$ 在宽范围的加载速率（跨越一个数量级）和不同加载模式（单轴拉伸、纯剪切、双轴拉伸）下近似为常数。
2. 开发了统一的损伤与断裂模型：建立了一个基于虚功率原理的多机制梯度损伤模型。该模型能够：
   • 预测损伤萌生、扩展及完全失效。
   • 适用于从粘性软聚合物（如 PBS、EPS25 玻璃体）到弹性体（如 EPDM、SBR、NR）的广泛材料。
   • 独立于加载历史和试件几何形状。
3. 微观 - 宏观跨尺度关联：展示了如何从宏观模型参数（$\psi^+_{cr}$）推导出微观结构物理量（如亚链解离能），并通过独立的流变学实验进行了验证。
4. 物理现象的预测能力：模型成功预测了软聚合物断裂中的关键物理现象，包括：
   • 脆性与延性响应的转变。
   • 不同加载速率下的断裂轮廓差异（粘性聚合物的“喇叭形”轮廓 vs. 弹性体的“抛物线形”轮廓）。
   • 裂纹的萌生、扩展及自止（Arrest）。

4. 实验结果与验证 (Results)

模型在多种材料和实验条件下进行了严格验证：

• Polyborosiloxane (PBS)：
  • 参数标定：利用单轴拉伸实验标定 $W_{cr}$ 和损伤参数。
  • 独立验证：在单边缺口拉伸、中心缺口拉伸及带孔拉伸（使用插入件）等多种几何形状和加载速率（2 mm/s 至 60 mm/s）下，模型准确预测了力 - 位移曲线。
  • 断裂轮廓：模型成功复现了快速加载下的“喇叭形”断裂轮廓（由粘性效应导致的二次缺口形成）和慢速加载下的高延展性（无断裂）。
  • 微观验证：基于 $\psi^+_{cr}$ 估算的亚链解离能（15 kJ/mol）与通过流变学实验测得的活化能（24 kJ/mol）处于同一数量级，验证了模型的物理一致性。
• EPDM 弹性体：
  • 利用 V 型缺口拉伸实验数据，模型独立预测了不同缺口深度和半径下的断裂位移，与实验结果吻合良好。
• EPS25 玻璃体 (Vitrimer)：
  • 在纯剪切加载下，模型预测了加载速率依赖的“流动 - 断裂”转变。特别是在中间加载速率下，模型成功预测了裂纹萌生、扩展随后停止（自止）的现象，这与实验观察一致。
• SBR 和 NR (天然橡胶)：
  • 模型预测的失效面（Failure Surfaces）与实验数据高度吻合，克服了传统变分相场模型在多轴应力状态下预测失效应力过大或不准确的缺陷。

5. 意义与展望 (Significance)

• 理论突破：该研究解决了速率相关材料断裂韧性度量非唯一性的长期难题，提供了一个基于能量密度的通用断裂阻力度量标准。
• 工程应用：所提出的模型具有广泛的适用性，可用于设计软聚合物在生物医学（药物递送、健康监测）、冲击防护、致动器及粘合剂等领域的应用，特别是在涉及复杂几何形状和动态加载的场景中。
• 物理洞察：通过将宏观断裂参数与微观化学键能联系起来，该模型为理解软材料断裂的微观机制提供了新的视角，实现了从连续介质力学到微观物理过程的跨尺度桥梁。
• 未来方向：未来的工作将集中在更复杂的试件几何形状、加载速率和加载模式下的进一步实验，以全面验证临界应力功的独立性，并继续完善软材料损伤与断裂的数值模拟工具。

综上所述，这篇论文通过提出临界应力功概念和开发统一的梯度损伤模型，显著推进了对软聚合物（特别是速率相关材料）断裂行为的理解和预测能力。