Predicting the Brittle-to-Ductile Transition in Amorphous Polymers

想象你有一块塑料，比如一把透明的尺子或一个坚硬的塑料瓶。如果你在温暖的日子里慢慢拉扯它，它会拉伸、弯曲，最终变细然后断裂。这是一种延性行为——它会给你预警。但如果你在寒冷的冬日里快速拉扯同一块塑料，它会伴随着一声清脆的咔嚓声瞬间断裂，完全没有拉伸。这就是脆性行为。

塑料从“可拉伸”转变为“易断裂”的那个临界点，被称为脆 - 延转变（BDT）。

本文旨在建立一套简单的数学“规则手册”，以根据拉扯速度和温度高低，精确预测不同类型的塑料何时会发生这种转变。

以下是作者解决这一难题的过程，用通俗易懂的语言讲述：

1. 问题：为什么我们需要新的规则手册？

科学家们早就知道，塑料的行为会因温度和速度的不同而改变。然而，过去缺乏一种简单、通用的方法来精确预测某种特定塑料究竟何时会断裂、何时会拉伸。现有的模型要么过于复杂，要么与数据不太吻合。

作者希望找到一个“临界点”。他们问道：在什么拉扯速度下，塑料会停止拉伸并开始断裂？

2. 核心思想：一场“能量耗散”的赛跑

把拉扯塑料想象成一场与时间的赛跑。

输入： 你通过拉扯（应变速率）向塑料注入能量。
输出： 塑料试图通过流动和重新排列其分子来消耗这些能量（粘塑性流动）。

只要塑料能够以足够快的速度重新排列分子以消耗你注入的能量，它就会平滑地流动（表现为延性）。但如果你拉得太快，塑料就无法及时重新排列自身。能量开始积聚，材料无法承受应力，于是发生断裂（表现为脆性）。

作者提出，脆 - 延转变恰好发生在塑料没有足够“时间”重新排列自身的时刻。

3. 塑料内部的“双速时钟”

为了理解塑料重新排列的速度，作者观察了两个控制分子运动的内部“时钟”（弛豫时间）：

大时钟（α-弛豫）： 这是主聚合物链缓慢、沉重的运动。就像一头大象试图在狭小的房间里转身。这通常控制着塑料在“玻璃化转变”（从坚硬变为橡胶状）附近的行为。
小时钟（β-弛豫）： 这是一种更快、更微小的摆动。就像大象的尾巴摇摆或耳朵扇动。作者发现，即使那头“大象”（大时钟）被冻结了，它的“尾巴”（小时钟）仍然可以摆动。

关键发现： 作者意识到，塑料只有在其分子重新排列的速度快于你的拉扯速度时，才能流动（表现为延性）。然而，这里存在一个速度极限。即使你无限快地拉扯，分子的摆动速度也只能受限于它们的“小时钟”（β-弛豫）。如果你拉扯的速度超过了这个极限，塑料别无选择，只能断裂。

4. “玩具模型”：弹簧与阻尼器

为了验证这一想法，作者建立了一个简化的数学模型（即“玩具模型”）。将一块塑料想象成由两部分组成的组合：

弹簧： 代表弹性部分（它想要弹回）。
阻尼器（减震器）： 代表流体部分（它缓慢流动）。

他们增加了一个转折：让“弹簧”具有非线性特征。想象一个弹簧，在达到某个点之前拉伸起来越来越容易，但随后会撞上一个“天花板”，超过这个点如果不破裂就无法继续拉伸。

随后他们问道：如果我们以不同的速度拉扯这个弹簧 - 阻尼器系统，它何时停止流动并开始破裂？

通过求解数学方程，他们创建了一个包含三个区域的相图（图表）：

区域 1（脆性）： 你拉得太快。系统无法流动。它断裂了。
区域 2（带有“打嗝”的延性）： 你以中等速度拉扯。塑料拉伸，变得稍微变软（出现“应力过冲”），然后稳定流动。
区域 3（类液态）： 你拉得非常慢。塑料顺畅流动，没有任何“打嗝”。

5. 理论测试：聚苯乙烯、PMMA 和 PVC

作者用三种常见塑料的真实世界数据测试了他们的模型：

聚苯乙烯（PS）： 用于 CD 盒和一次性餐具的材料。
PMMA（有机玻璃）： 透明、抗破碎的玻璃替代品。
PVC： 用于管道和 plumbing 的材料。

他们发现，他们的模型效果出奇的好。

“效率”因素： 他们发现，不同的塑料在利用内部摆动（β-弛豫）来软化方面具有不同的“效率”。
- PMMA 和 PVC 效率很高。在受力时，它们几乎可以完全“熔化”其内部结构以进行流动。这使得它们不那么脆。
- 聚苯乙烯（PS） 效率较低。即使在受力时，其大部分结构仍然保持“冻结”和刚性。这就是为什么即使在相似的温度下，PS 比其他材料更脆、更容易断裂的原因。

6. 结论

该论文声称，预测塑料何时断裂并不需要复杂的断裂力学。相反，你只需要知道：

塑料分子能有多快摆动（β-弛豫时间）。
你拉扯它的速度有多快。

如果你拉扯的速度快于分子摆动的速度，塑料就会变脆。如果你拉扯的速度较慢，它就会流动。作者的模型成功预测了不同塑料的这种“临界点”，并与真实世界的实验结果相符。

简而言之： 该论文提供了一条简单、通用的规则：当你拉扯塑料的速度快于其分子摆动的速度时，塑料就会断裂。

技术摘要：预测非晶态聚合物的脆韧转变

问题陈述
脆韧转变（BDT）是非晶态和半结晶聚合物的关键特征，决定了材料是在低应变下发生脆性断裂，还是经历涉及应变软化、屈服和应变硬化的韧性变形。虽然实验数据常将 BDT 与 $\beta$ -转变（特别是 Johari-Goldstein 过程）联系起来，但目前尚无简单统一的模型能根据聚合物化学性质、样品历史、变形类型和应变率来预测 BDT。现有的理论方法通常依赖于比较屈服应力与脆性断裂应力，但这种方法因断裂固有的非均匀性以及估算断裂应力的困难而变得复杂。此外，近期对 Eyring 方法的批评指出，标量应力描述违反了对称性要求，因此需要基于弹性能量密度的模型。

方法论
作者提出了一种唯象的标量“玩具模型”（一维），用于描述粘弹塑性（VEP）剪切应力 - 应变曲线。该方法整合了以下几个关键组成部分：

广义麦克斯韦模型：材料被建模为麦克斯韦元件（弹簧和阻尼器串联）。弹性分量采用非线性 Frenkel 风格的余弦势来描述能景，允许在屈服前出现弹性应力最大值。粘性分量遵循麦克斯韦流动，其中应力松弛由依赖于应变率的松弛时间控制。
松弛时间降低：模型摒弃了传统的基于应力的 Eyring 降低方法，转而采用 Long 等人的方法，即松弛时间的降低与弹性应变能密度成正比。这确保了旋转不变性。
TS2 模型集成：模型利用“Sanchez-Lacombe 双态双（时）标”（SL-TS2）框架来计算 $\alpha$ -和 $\beta$ -松弛时间的温度依赖性。该框架将材料视为“固态”和“液态”区域的动态共存，其中固态分数（ $\psi$ ）取决于温度和老化历史（虚构温度）。
稳定性分析：BDT 的定义并非通过比较断裂应力和屈服应力，而是通过识别均匀、均质的粘塑性流动变得不稳定的点。这类似于在热力学中寻找旋节线。模型假设，当应变率超过材料最大能量耗散能力时，均匀流动就会失效，而该能力受限于 Johari-Goldstein $\beta$ -松弛时间。

主要贡献

相图构建：作者在 $(B, Y)$ $(B, Y)$ 空间中推导出了无量纲相图（其中 $B$ $B$ 与逆温度和相关体积有关， $Y$ $Y$ 与应变率有关）。该图划分了三个区域：
- I 型：脆性断裂（无法实现均匀流动）。
- II 型：具有应力过冲和屈服的韧性行为。
- III 型：具有屈服但无应力过冲的类液态行为。
统一的 BDT 判据：论文确立了均匀粘塑性流动的应变率上限与 $\beta$ -松弛时间成反比。由此导出了 BDT 活化能公式，该公式恢复了 Wu 的唯象方程，并将其解释为 $\alpha$ -和 $\beta$ -过程贡献之间的混合规则。
老化考量：模型通过计算虚构温度并调整固态分数（ $\psi$ ）和 $\alpha$ -松弛时间，明确考虑了物理老化，承认老化样品表现出不同的 BDT 行为。

结果
该模型应用于三种非晶态聚合物：聚苯乙烯（PS）、聚甲基丙烯酸甲酯（PMMA）和聚氯乙烯（PVC）。

与实验的一致性：模型预测的 BDT 温度随应变率的变化与 Wu 汇编的实验数据吻合良好。
材料特异性：一个关键发现是“效率”参数（ $x$ $x$ ）的差异，该参数代表变形过程中“熔化”的固态区域的比例。
- 对于 PMMA 和 PVC， $x \approx 1$ ，表明材料在屈服点附近表现为深度过冷液体，BDT 发生在更接近 $\beta$ -转变的位置。
- 对于 PS， $x \approx 0.5$ ，表明即使在屈服时仍有相当比例的固态区域存在。这导致 BDT 温度更接近玻璃化转变温度（ $T_g$ ），与实验观察一致，即 PS 比 PMMA 或 PVC 更脆。
活化能：计算得出的 BDT 活化能（ $E_b$ ）与实验值相当吻合（例如，PMMA 为 40 kcal/mol，实验值为 44 kcal/mol；PS 为 85 kcal/mol，实验值为 77 kcal/mol）。

意义与主张
该论文声称提供了一种简单的标量模型，能够在无需复杂断裂力学计算的情况下成功预测 BDT。通过将 BDT 框架化为均匀粘塑性流动的稳定性丧失（一种“旋节”方法），作者提供了一种比传统屈服应力与断裂应力比较法更高效的计算替代方案。该模型成功地将宏观机械失效与微观松弛过程（ $\alpha$ 和 $\beta$ ）及样品历史（老化）联系起来。

作者对模型的适用范围保持谦逊，指出它是一个局限于应变硬化前区域的“玩具模型”，并假设体积守恒。他们承认局限性，例如缺乏显式的体积变化效应，以及需要未来的工作来描述具有极低脆韧转变温度的聚合物（如聚碳酸酯），并模拟包含应变硬化的完整应力 - 应变曲线。其主要贡献在于证明，基于稳定性的方法结合 TS2 松弛时间框架，可以定量描述不同聚合物化学性质下的 BDT。