玻尔百科虽然我们的直觉认为材料在变形时会变强，即所谓的应变硬化现象，但许多材料表现出相反的行为：应变软化。这种反直觉的特性，即材料在超过某一变形点后会变弱，对经典力学和计算工程学构成了深远的挑战。它引入了一种传统理论难以预测的基本不稳定性，导致灾难性的失效和不可靠的模拟。本文旨在解决应变软化的复杂性，为工程师和科学家提供一个全面的概述。第一部分“原理与机制”将解构该行为背后的理论，解释材料不稳定性、应变局部化以及网格依赖性的计算危机等概念。随后，“应用与跨学科联系”部分将展示应变软化的现实世界影响，从在岩土力学中引发流动液化到在高分子材料中实现可控强化，并探讨那些已经成功驾驭这一复杂现象的先进建模技术。
在我们与材料的日常接触中，我们形成了一种特定的直觉。如果你拿一个金属回形针来回弯折，它会变得越来越难弯。铁匠锤打一块发光的钢，每一次捶打都使金属变得更强、更韧。这种被称为加工硬化或应变硬化的现象感觉很自然：使材料变形使其更能抵抗进一步的变形。但如果情况相反呢？如果超过某个点后，拉伸材料反而使其变得更弱呢？
这不是一个刁钻的问题。这是许多材料真实且极其重要的一种特性，称为应变软化。它描述了材料强度随变形增加而丧失的现象。想象一下应力-应变曲线，这是材料力学响应的基本“特征”。我们施加一个应力 ，并测量由此产生的应变 。对于许多材料，这条曲线最初是上升的——需要更大的应力才能产生更大的应变。在塑性区，曲线可能继续上升，这就是应变硬化。但对于表现出应变软化的材料，在达到峰值应力后，曲线开始向下倾斜。继续使材料变形所需的力量实际上在减小。
在某些高分子材料中可以清楚地看到这种现象的一个绝佳例子。取一块半结晶高分子，如聚乙烯，在略高于其玻璃化转变温度（）的温度下进行测试。在这个温度下，高分子的非晶部分呈橡胶态且易于移动。当你开始拉伸它时，它最初会弹性抵抗。然后它达到一个屈服点，接着发生了一些显著的变化。一个“颈缩”形成了，这是一个明显变细的局部区域。随着这个颈缩的扩展，保持样品拉伸所需的工程应力实际上会下降。材料正在软化。从数学上讲，我们说它的切线模量 变为负值。这个负斜率是应变软化的标志，正如我们即将看到的，它预示着一系列迷人而又具挑战性的问题的开始。
图上的负斜率似乎并不那么令人担忧。但在力学世界里，它是一个深层且根本的材料不稳定性形式的警示信号。几十年来，塑性理论——研究材料如何永久变形的科学——建立在一套确保材料以“稳定”和可预测方式表现的基础原则之上。其中最重要的原则之一是Drucker稳定性公设。
本质上，Drucker公设是关于功和能量的陈述。其一种形式是，为了发生塑性变形，应力的任何微小增加在所产生的塑性应变增量上所做的功必须为非负。用数学简写形式表示为 。这似乎很合理；它有助于保证一个结构问题的解是唯一的。遵守这一定律的材料是表现良好的。
然而，应变软化材料是“叛逆者”。它们公然违反了这条规则。在软化期间，应力在减小（），而塑性应变在增加（）。它们的乘积是负的。根据Drucker的定义，该材料是不稳定的。这不仅仅是一个学术上的奇闻。它告诉我们，我们传统的理论工具可能会失效。结构对载荷的响应存在单一、唯一方式的假设不再得到保证。材料进入了一种可以自己选择失效路径的状态。
那么，这种不稳定性在物理世界中是什么样子的呢？答案是一种叫做应变局部化的现象。
想象一条由一百个相同链环组成的链条。如果你拉动链条，每个链环都会分担载荷并略微伸长。现在，想象一条应变软化的链条。当你拉动时，一个链环——也许是由于微小的缺陷——首先开始软化。它的抵抗力下降了。由于每个链环必须承受相同的力，这个稍弱的链环被迫伸长得更多以作补偿。这种增加的伸长导致它进一步软化，使其变得更弱。一个恶性循环开始了。很快，整个链条中所有后续的变形都将“局部化”到这一个链环中，它会急剧伸长直到断裂，而其他99个链环几乎没有变形。
同样的原理也适用于连续材料。变形不再均匀分布，由软化引起的不稳定性迫使所有应变集中在一个极其狭窄的带中。在这个带之外，材料甚至可能松弛并弹性卸载。一个简单杆在拉伸下的控制方程是平衡条件，该条件要求应力 沿其长度是均匀的。当一个小区域开始软化时，其承载应力的能力下降。为了与邻近区域保持相同的应力，该区域内的应变必须急剧增加。
这种行为对经典工程理论造成了严重破坏。例如，极限分析是一种用于计算结构在坍塌前能承受的最大载荷的强大方法。它依赖于稳定塑性行为的假设。对于应变软化材料，极限分析可能会给出极其错误的答案。在一个引人注目的思想实验中，可以证明对于一个具有哪怕是无限小的、预软化区域的杆，其理论坍塌载荷会变为零。理论预测在任何载荷下都会发生灾难性失效，这仅仅是因为它无法处理局部化的病态性质。
当我们尝试使用有限元法（FEM）等工具在计算机上模拟这些材料时，问题变得更加离奇。有限元模拟将一个连续物体分解成一个由离散单元组成的网格。它如何处理应变局部化呢？
问题的关键在于：一个简单的、或称为“局部的”应变软化模型没有尺度感。本构律——即应力与应变关系的法则——在每一点都是相同的。它不包含任何关于特征长度的信息。因此，当材料变得不稳定并决定局部化时，应变会去哪里？它没有物理上的指引，所以它遵循它所拥有的唯一指引：数值网格。应变会集中到模拟所能表示的最小可能区域，这通常是一排有限元。
这意味着计算出的失效区宽度 不是材料的物理属性，而是由网格单元的大小 决定的。这被称为病态网格依赖性。如果你为了得到更“精确”的结果而细化网格（使用更小的单元），局部化带只会变得更窄。你计算的全局载荷-位移曲线会随着每个网格而改变。模拟永远不会收敛到一个唯一的、物理上正确的答案。
这场灾难背后深层的数学原因是，系统的控制偏微分方程（PDE）变得不适定（ill-posed）。一旦软化开始，系统就失去了一个称为椭圆性的关键性质。控制稳定弹性问题的椭圆方程倾向于使事物平滑化。失去椭圆性就像关闭了这种平滑属性；方程突然允许出现尖锐的、不连续的解，这恰恰是应变局部化的表现。
这种数字病态最惊人的后果与能量有关。完全断裂一种材料所需的能量，即其断裂能，应该是一种基本的材料属性，就像其密度或熔点一样。让我们尝试从我们的模拟中计算这个能量。
总耗散能量是单位体积耗散的比能（与应力-应变曲线下的面积有关）在局部化带的体积上积分。但我们刚刚确定，在局部软化模型中，这个带的宽度 与单元尺寸 成正比。因此，局部化区的体积也与 成正比。这导致了一个真正荒谬的结论：失效过程中总的耗散能量与网格尺寸 成正比。
这意味着当你细化网格，使 以期获得更精确的解时，计算出的断裂材料所需的总能量会虚假地消失！
模拟预测破坏该物体不耗费任何能量。模型创造了一个数字幽灵，一个看起来真实但没有能量实质的失效。这不仅仅是一个数值误差；它是一个迹象，表明我们的物理模型从根本上是不完整的。
问题出在哪里？我们简单的“局部”连续介质模型，即一个点的应力仅取决于该精确点的应变，是一个方便但有缺陷的理想化。真实的材料不是抽象的、点状的连续体。它们拥有微观结构：金属中的晶粒、混凝土中的骨料颗粒和孔隙、高分子中的缠结链。
失效过程，如微裂纹或剪切带的形成，发生在称为“断裂过程区”的有限物理区域内，其尺寸由这些微观结构特征决定。一个物理上忠实的连续介质模型必须以某种方式包含这些信息。它必须被赋予一个内禀长度尺度，一个通常用 表示的材料参数，它告诉模型底层物理过程的特征长度。
这是如何做到的呢？现代力学设计了几种优雅的策略，通常称为正则化技术：
梯度增强模型： 本构律被丰富，使得应力不仅依赖于应变，还依赖于其空间梯度（例如 ）。这引入了一个惩罚应变剧烈变化的项，有效地迫使任何局部化“弥散”在一个有限的带内，其宽度与 成正比。
非局部模型： 一个点的状态不再由局部决定。相反，它可能被定义为该点周围小邻域内状态的加权平均。这个平均邻域的半径由内禀长度尺度 决定。
通过引入这个物理长度尺度 ，我们对数学问题进行了正则化。即使在软化期间，控制方程也保持适定性。在模拟中，局部化带的宽度现在由材料的物理特性（）控制，而不是分析师的网格（）。随着网格的细化，计算出的载荷、位移以及——最重要的是——耗散能量的结果会收敛到一个唯一的、有限的、物理上有意义的答案。模拟变得具有网格客观性。
这个关于应变软化的故事——从一个简单的向下倾斜的曲线到一个深刻的计算危机及其优雅的物理解答——完美地展示了科学是如何进步的。我们模型中的一个悖论迫使我们更深入地探究材料的本质，最终导向一个更深刻的理论，它弥合了晶粒和空隙的微观世界与工程结构的宏观世界之间的鸿沟。
见一物，如见一世界。我们已经穿行于应变软化的抽象原理之中，观察到材料如何在其变形过程中悖论般地变得更弱。这个简单的想法，这种强度的丧失，似乎仅仅是失效的前奏。但当我们仔细观察时，我们发现它是一粒种子，从中生长出一片复杂、迷人且极其重要的现象森林。它的枝干从我们超级计算机的计算核心延伸到塑料的寂静缓慢老化，地球的剧烈震动，以及形成裂纹中原子的精妙舞蹈。
让我们从计算机世界开始我们的旅程，应变软化正是在这里首次显露其诡谲的本性。想象一下，你编写一个程序来模拟一根被拉断的简单金属棒。材料模型包含了软化。你运行模拟，它预测棒将会断裂。但随后你用更精细的网格，一个更详细的计算网格再次运行它。结果变了。你每次细化网格，结果都会改变。失效区域变得越来越窄，而预测的断裂棒所需的总能量——我们称之为断裂能，，本应是材料的基本属性——却虚假地趋近于零。这是一场被称为“病态网格依赖性”的灾难。这意味着我们的模拟告诉我们的是其自身的网格信息，而不是材料的物理特性。
这台机器中的幽灵是长度尺度的缺失。在简单的模型中，软化可以局部化到一个零厚度的区域，一条不符合物理的数学线。那么，解决方法就是告诉模型“一条裂纹的代价”。断裂能 是创造单位面积新的、破碎表面所需的真实物理能量。我们的模拟必须尊重这个预算，无论网格尺寸如何。“裂纹带模型”是一个非常实用的修正方法：我们测量计算单元的尺寸，比如说 ，然后巧妙地调整材料的软化率，以确保在宽度为 的局部化带内耗散的总能量始终等于 。这个过程被称为正则化，它确保了无论我们的模型使用理想化的线性软化定律还是更现实的指数定律，能量核算总是正确的。
这个想法可以用热力学的语言建立在更坚实的基础上。我们可以从亥姆霍兹自由能的角度来思考材料的状态，它是应变和追踪退化程度的“损伤变量” 的函数。损伤的驱动力是能量释放率 ，也就是可用于产生新损伤的能量。总耗散能量就是这个驱动力在整个损伤过程中的积分。裂纹带模型的核心策略，是确保这个热力学定义的耗散，在局部化带上积分后，与物理测量的断裂能 相匹配。 这是计算必要性与热力学原理的美妙结合。
软化的后果在任何地方都不如在我们脚下的大地中那般戏剧性。土壤、黏土和岩石不是简单的弹性固体；它们是复杂的颗粒材料，既可以硬化，关键是也可以软化。
考虑一个建在密实黏土层上的建筑基础。随着荷载增加，黏土发生变形。如果一片黏土区域开始软化，其承载能力就会下降。这部分荷载必须转移到相邻区域，而这些区域又可能被推向其自身的软化点。这可能导致全局的结构不稳定性。整个基础的承载能力可能会达到峰值，然后随着进一步的沉降而开始下降。在荷载对沉降的图上，曲线会向后弯曲，这种现象被称为“突弹”（snap-back）。一个缓慢增加荷载的简单模拟会在峰值处灾难性地失败，无法跟随结构走下不稳定的路径。为了追踪这种复杂的失效之舞，计算工程师必须采用复杂的“弧长法”，将荷载和位移都视为未知数，并沿着结构的平衡路径进行追踪，即使路径向后延伸。
在地震期间，风险变得更高。在地震波的快速、反复剪切作用下，一层饱和沙土可能表现出一种可怕的转变。松散的、收缩性的沙土倾向于压实，但由于其孔隙中的水无法迅速排出，孔隙水压力急剧上升。这个压力将沙粒推开，破坏了它们之间的有效应力，导致材料几乎丧失所有强度。此时，土壤自重产生的剪应力可能超过其残余强度，引发灾难性的、不受控制的变形——流滑。这就是流动液化，一种经典的应变软化不稳定性。 相比之下，更密实的、剪胀性的沙土表现则不同。虽然它也可能累积孔隙压力，但剪切时膨胀（剪胀）的趋势提供了一种救命的负反馈。当它变形时，它会反推水，降低压力并恢复部分强度。它不会灾难性地流动，而是在一个称为循环流动性的过程中经历大的循环变形。土壤的初始状态——它的密实度——决定了软化是导致完全坍塌，还是导致一场剧烈但能自我限制的舞蹈。
软化总是毁灭的预兆吗？大自然以其优雅，常常将明显的弱点转化为力量。高分子世界就是最好的例证。
拿一块简单的塑料，比如聚乙烯，然后慢慢拉伸它。在它屈服之后，会形成一个“颈缩”，材料在此处变薄。但它并不会立即断裂，这个颈缩通常会稳定下来，并以几乎恒定的力沿着样品传播。这种非凡的现象，“冷拉伸”，正是应变软化在起作用。未颈缩的材料软化并流入颈缩区，在那里的高分子链变得高度排列和取向。这种新的、拉伸过的材料实际上比原始材料更强、更硬。应变软化是实现从无序、较弱状态到有序、较强状态这种可控转变的机制。
当我们考虑像聚苯乙烯或聚碳酸酯这样的玻璃态高分子的隐藏生命时，故事变得更加深刻。一块塑料的性质不是固定不变的；它们会随时间演化。当高分子被冷却到其玻璃态时，它被困在一个高能量、无序的排列中。随着时间的推移，即使在室温下，它也会经历“物理老化”，这是一种向更密实、更低能量状态的缓慢结构弛豫。这个老化过程对材料的力学响应有深远的影响。一个高分子老化得越久，它对屈服的抵抗力就越强，因此其屈服应力增加。因为老化状态比变形状态有序得多，屈服后应力的下降——即应变软化的深度——变得更加显著。甚至对银纹（断裂的前兆）的抵抗力也随着老化而增加。我们所感知的静态材料实际上是一个动态的实体，它的历史被写在它应变软化行为的细节里。[@problem_g_id:2937928]
务实的“裂纹带”修正法虽然有用，但它承认了我们底层的局部模型是有缺陷的。更深层的真相是，材料行为不是局部的。材料中的一个点不会孤立地决定失效；它会与它的邻居交流。最先进的软化模型拥抱了这种非局部性。
积分型非局部模型假定，一个点的损伤趋势是其周围邻域平均应变状态的函数。这种平均化将一个物理长度尺度 直接引入到控制方程中。现在，局部化带的宽度由这个内禀材料参数决定，而不是人为的网格尺寸。近场动力学（Peridynamics）将这一思想更进一步，完全摒弃了应力的局部概念，并根据点与点之间的长程力重新构建了力学。 一种替代的、在数学上相关的途径是梯度增强模型。在这里，长度尺度是通过在方程中添加涉及应变或损伤的空间梯度的项来引入的。对于缓慢变化的场，可以证明积分方法和梯度方法是完美等价的，揭示了对同一底层物理的两种不同视角。
这些先进的模型不仅仅解决了一个数值问题；它们带来了更深刻的计算见解。损伤定律的精确数学表述——无论它是由能量类量还是应变类量驱动——决定了系统切线刚度矩阵的对称性。对称矩阵在计算上是“友好的”，允许使用像共轭梯度法这样高效的求解器。非对称矩阵则需要不同的、通常更昂贵的算法。当软化发生时，矩阵还会失去正定性，变为“不定”的，这又迫使求解策略发生改变。因此，软化的物理特性与模拟它所需的数值算法的选择紧密相连。
最后，软化的挑战在建模的各个尺度上泛起涟漪。在现代的“多尺度”模拟（FE²）中，我们试图通过在每个点上对其微观结构进行微小模拟来计算大型结构的行为。但如果微观结构本身表现出应变软化会发生什么？小尺度上的不稳定性可能会毒害整个事业。它打破了尺度分离的原则，并违反了基本的Hill-Mandel能量一致性条件，导致宏观模拟变得和我们最初的简单模型一样具有网格依赖性。一个局部的不稳定性可以产生真正全局性的后果，挑战我们最先进模拟技术的根本框架。
从一个计算上的小故障到一个创造更强材料的机制，从一个地基的坍塌到新数值方法的设计，在应变下变弱这个简单的行为，展现为一个极其丰富和复杂的主题。它提醒我们，在自然界中，如同在科学中一样，最深刻的洞见往往不是在力量中找到，而是在弱点那微妙而错综复杂的美中发现的。