玻尔百科为什么有些材料会弯曲，而另一些则会折断？理解材料失效背后的机理是现代科学与工程的基石，对于设计从高弹性建筑到安全飞机的各种结构都至关重要。当延性金属被推向极限时，它们会表现出两种截然不同但又相互交织的行为：塑性，即形状的不可逆变化；以及损伤，即导致断裂的内部退化。孤立地处理这些现象的模型是不完整的；它们无法捕捉材料变形、弱化并最终断裂的全过程。为了真正预测失效，我们必须理解塑性与损伤在一个耦合系统中的协同作用。
本文将深入探讨耦合损伤-塑性的统一理论。第一部分原理与机理将奠定理论基础。我们将探讨塑性与损伤各自独特的物理特征，引入将它们在数学上联系起来的强大工具——应变等效原理，并考察变形驱动损伤、损伤反过来影响变形的反馈循环。第二部分应用与跨学科联系将展示该理论的实际应用价值。我们将看到这些模型如何通过真实世界的实验进行校准，并用于计算机模拟，以预测材料从循环疲劳到灾难性断裂的整个生命周期，揭示抽象方程与具体工程成果之间的深刻联系。
想象一下，你拿一个金属回形针来回弯折。它不会弹回原来的形状，而是保持弯曲。现在，再想象一下拿起一根粉笔并试图弯曲它。它不会弯曲，而是会折断。这些日常经验揭示了材料在受力过大时发生永久性变化的两种根本不同的方式。在材料科学领域，我们称这些变化为塑性和损伤。理解它们——无论是分开理解还是作为相互交织的伙伴来理解——是预测物体如何以及何时会断裂的关键。
让我们先分别看看这两个角色。
当我们的回形针弯曲并保持弯曲状态时，它经历了塑性。这是一种材料的不可逆流动。在微观层面，原子平面相互滑移，就像一副被剪切的扑克牌。塑性最显著的标志是在你松开力后留下的永久应变。如果你稍微拉伸一根金属棒，它的行为就像弹簧——它是弹性的。如果你拉伸得更远，它就开始塑性变形。当你释放它时，它会回弹一点，但会比原来永久性地变长。这部分多出来的伸长就是塑性应变，我们可以称之为。纯塑性的一个关键且可能令人惊讶的特点是，它不一定会使材料潜在的弹性性质“变弱”。如果你轻轻敲击弯曲的回形针，它的弹性刚度——其抵抗像弹簧一样被拉伸的内在能力——或多或少与之前相同。
损伤，另一方面，则是粉笔的故事。它关乎的不是流动，而是断裂。损伤代表了材料完整性的丧失。在微观层面，它是微小孔洞和裂纹的产生与增长。与塑性不同，损伤直接攻击材料承载载荷的能力。其最重要的后果是弹性刚度的降低。
想象一根由一百股绳索组成的粗绳。它的刚度是衡量其在给定重量下拉伸程度的指标。现在，如果你偷偷剪断了其中二十股，会发生什么？绳子被损伤了。下一次有人拉它时，它会更容易被拉伸，因为承载载荷的绳股变少了。它的刚度下降了。我们可以通过加载材料，然后在小的弹性范围内卸载并重新加载来测量这种刚度退化。那条应力-应变曲线的斜率会比原始材料的斜率更平缓。另一种“看到”损伤的方法是测量声波在材料中传播的速度。声波本质上是一种弹性压缩和膨胀的波。在受损、刚度较低的材料中，这种波传播得更慢，就像波在松弛的绳子上比在绷紧的绳子上跑得慢一样。
对于许多材料，比如粉笔，故事到此为止。但对于延性金属——构成飞机、汽车和建筑的材料——塑性与损伤是密不可分的舞伴。纯塑性模型无法预测断裂，而纯损伤模型则无法捕捉延性金属在失效前经历的巨大永久变形。要真正理解它们的行为，我们必须理解这两个过程是如何耦合的。
我们如何能构建一种单一、连贯的数学语言来描述一个既在塑性变形又在累积损伤的材料呢？答案在于一个极其简单而强大的思想：应变等效原理。
让我们回到那根100股的绳子。如果挂上100磅的重物，平均每股绳索承受1磅的力。但如果剪断了20股（即20%的面积损失），同样的100磅重物现在仅由80股支撑。作用在那些剩余、完好绳股上的应力现在是磅。作用在材料仍在工作的有效部分上的有效应力，高于你通过将重量除以绳子原始总面积计算出的平均应力或名义应力。
这正是该理论的精髓所在。我们定义一个标量变量，即损伤变量，对于原始材料，，对于完全失效的材料，。如果，意味着材料已损失了20%的有效承载面积。由名义应力与作用在材料未损伤骨架上的“真实”应力，即有效应力，通过一个简单的公式联系起来：
应变等效原理指出，只要我们用这个有效应力来表达，损伤材料的本构律（应力与应变之间的关系）与未损伤材料的本构律具有完全相同的形式。如果原始的、未损伤的材料遵循胡克定律，（其中是初始刚度张量，是弹性应变），我们现在写成：
代入我们对有效应力的定义，我们得到。重新整理得到损伤材料的应力-应变定律：
看看这个方程告诉了我们什么！材料的表观刚度不再是，而是一个退化了的刚度。这个从关于面积损失的简单物理概念中得出的优雅结果，完美地解释了我们之前讨论的卸载-再加载斜率的降低和声波速度的减慢。
现在我们已经为观察损伤与塑性如何相互作用搭好了舞台。这是一个反馈循环：塑性可以导致损伤，而损伤可以影响塑性。
首先，损伤如何影响塑性流动？当材料中的应力达到一个临界值，即屈服应力时，塑性开始。但是，是哪种应力呢？名义应力还是有效应力？物理上的屈服——原子平面的滑移——发生在材料的完好部分。因此，必须是有效应力达到屈服准则。
想象一下，屈服准则是在所有可能应力空间中的一个边界，一个“屈服面”。如果你的应力状态在这个边界内部，材料是弹性的。一旦应力状态达到边界，塑性流动就开始了。因为对于一个损伤的材料，大于，所以材料会在一个更低的名义应力水平上达到这个临界边界。换句话说，随着损伤的增加，在名义应力空间中所看到的屈服面缩小了！材料似乎变得“更软”，更容易屈服。损伤为更多的塑性流动准备了条件。
但损伤最初是从哪里来的呢？在延性金属中，正是塑性流动本身催生了损伤。这是反馈循环的另一面。对此的经典模型是Gurson模型，它将损伤看作是金属基体中微观球形孔洞的集合。
思考一下这个惊人的事实：如果你拿一个坚固、致密的金属块，并对其施加均匀的静水拉伸（在所有方向上均匀拉伸，就像深海压力的负值），它不会发生塑性屈服。但如果金属中有微观孔洞，它会屈服！孔洞周围的材料会开始流动，导致孔洞扩张。Gurson模型在宏观屈服行为和孔洞体积分数（损伤变量的一种物理表现）之间提供了一个优美的数学联系。它表明，随着孔隙率的增加，引起屈服所需的应力会急剧下降。这为抽象的损伤概念及其演化提供了具体、物理的来源。
损伤，就像打破一个鸡蛋，是一个不可逆的过程。损伤变量只能增加或保持不变；它永远不会减少（除非是奇异的自愈合材料）。我们写作。但一个完整的理论不仅必须说明损伤可以发生，还必须说明是什么驱动它。使微裂纹生长的“力”是什么？
一个诱人但最终有缺陷的想法是将损伤演化直接与应变联系起来。一个基于应变的模型可能会说，例如，一旦材料中的总应变超过某个临界值，损伤就开始增长。这种方法的问题在于，对于延性金属，总应变主要由塑性主导。该模型错误地暗示，仅仅是大量弯曲材料就会导致损伤，将“流动”与“断裂”混为一谈。
一种更为深刻且有物理基础的方法是基于能量的方法。在我们理论的热力学框架内，我们可以定义一个称为损伤能量释放率的量，用表示。你可以将看作是可用于“释放”以创建新裂纹表面的储存弹性能量的量。它是与损伤变量共轭的真正热力学力。基于能量的演化法则规定，只有当这个驱动力达到一个临界阈值时，损伤才会增长。
这不仅仅是理论上的优雅问题。它具有巨大的实际意义。创建单位面积新裂纹表面所需的能量是一个基本的材料属性，称为断裂能。基于能量的损伤模型是唯一可以直接且有意义地校准到这一物理属性的模型，从而能够对材料失效进行更可靠和更具预测性的模拟。
到目前为止，我们一直用一个单一的数字来处理损伤。这暗中假设损伤是各向同性的——即在所有方向上都相同。材料的刚度降低了，但在各个方向上降低的程度是相等的。但这现实吗？
想象一下，拿一块金属板，在x方向上施加的拉力是y方向的两倍。很自然地，我们预期形成的微观孔洞和裂纹会倾向于垂直于最大应力方向。损伤本身将是各向异性的。
如果我们随后探测该材料的刚度，我们应该会发现它在x方向上的降低远大于在y方向上。这个最初是各向同性的材料，由于其加载历史而变得各向异性。这正是精细实验所显示的。一个简单的标量损伤模型，其中，从根本上无法捕捉这种现象。如果你校准它以匹配x方向上测得的刚度，它在y方向上就会出错，反之亦然。
这不是理论的失败，而是对其进行改进的呼唤。它告诉我们，我们对损伤的描述过于简单。解决方案是将损伤变量从一个简单的标量提升为一个损伤张量。张量是一个数学对象（可以由矩阵表示），它具有针对不同方向的分量。通过损伤张量，我们可以描述这样一种状态，例如，x方向的损伤大于y方向的损伤。这样的模型可以成功预测和，完美地匹配了诱导各向异性的实验观察。
这段从简单的标量概念到更复杂张量必要性的旅程，是科学如何进步的完美例证。我们基于物理直觉建立简单、优美的模型，用观察来检验它们，当它们不足时，我们利用其失败的本质作为指导，来构建一个更丰富、更完整、更强大的对世界的描述。
既然我们已经熟悉了耦合损伤与塑性的原理和机理，我们可能会想坐下来欣赏该理论的数学优雅。但那样就完全错过了重点！这些思想的真正美妙之处不在于其抽象的表述，而在于它们描述、预测并最终控制我们周围世界的力量。这些方程不仅仅是黑板上的装饰；它们是我们用来提出并回答科学与工程中一些最关键问题的工具：“它会断裂吗？何时断裂？如何断裂？”
本理论使我们能够做一件了不起的事：为材料构件构建一个“虚拟孪生体”，将其置于计算机中，并在数小时或数分钟内对其进行一生的使用和滥用。现在，让我们踏上旅程，看看这个虚拟实验室能让我们发现什么。
在我们信任计算机模拟之前，我们必须首先让模型了解我们感兴趣的特定材料。理论是普适的，但材料是独特的。每种合金都有其自身的特性，其强度与脆性的故事。这个“教导”模型的过程，是理论家的方程与实验家的实验室工作台之间一场引人入胜的对话。
那么，我们究竟如何在一块真实的金属中“看到”塑性应变和损伤这些抽象的量呢？想象一下，我们拿起一根金属棒，在试验机中拉伸它。我们可以测量施加的力以及棒的伸长量。起初，棒的行为就像一个简单的弹簧。但如果我们拉得足够用力，情况就变了。如果我们随后释放载荷，会发现棒并没有恢复到其原始长度；它产生了一段永久的、不可逆的伸长。这个在零应力下测量的残余应变，正是塑性应变的物理表现。这是材料经历磨难后留下的永久伤疤。
但损伤呢？损伤更为微妙。它是微观孔洞和裂纹的累积，一种内部的衰退。探测它的最直接方法是测量材料的刚度。把它想象成给病人测脉搏。一个健康的材料有其特定的弹性刚度，我们称之为。当我们加载和卸载材料时，我们可以测量卸载过程中应力-应变曲线的斜率。这个卸载是一个弹性过程，其斜率揭示了材料的当前健康状况。如果我们发现这个卸载刚度下降了，就意味着材料变得“更软”了。这种刚度的降低正是损伤的标志性特征。损伤变量不过是对这种刚度损失的精确量化[@problemid:2924574]。
掌握了这些实验技术，我们就可以开始模型校准的复杂过程。我们从一个简单的拉伸试验开始。在试样开始出现明显“颈缩”之前的初始均匀拉伸阶段，材料通过应变硬化变得更强。通过将我们的测量值转换为“真实”应力和应变，我们可以分离出这种硬化行为，并拟合我们塑性模型的参数。然而，一旦棒材开始颈缩，情况就完全改变了。应力状态变得复杂且非均匀。我们在测量中看到的表观软化，现在是真实材料退化（损伤）与横截面收缩的简单几何效应的混乱混合。要解开这两者，需要一种更复杂的方法。通常，工程师必须采用一种“反演方法”，即一遍又一遍地运行整个拉伸试验的计算机模拟。计算机会智能地猜测损伤参数，运行模拟，将结果与实验数据进行比较，然后不断优化其猜测，直到模拟的力-位移曲线与真实曲线完美匹配。这个过程本身就是一门艺术，是实验数据与计算能力协同作用的证明。
一旦我们有了经过校准和验证的模型，真正的乐趣就开始了。我们现在可以用它来预测材料在远比简单拉伸复杂得多的情景下的响应。
任何模型的第一个测试，就是看它是否能重现用于校准它的实验。确实，一个正确配置的耦合损伤-塑性模型可以自行描绘出拉伸棒的整个生命故事：初始的弹性拉伸，屈服的开始，材料反抗的硬化阶段，最终的峰值载荷，以及当损伤占主导地位、材料软化走向完全失效时的最后悲剧性衰退。
更重要的是，我们可以探索材料在循环加载下的命运，这是疲劳分析的核心。为什么回形针来回弯折最终会断裂？每一次弯曲循环都会引起一点塑性变形，并产生微量的新损伤。在应力-应变图上，这个过程描绘出一个“滞回环”。这个环的面积代表了每个循环中损失的能量，这些能量转化为热量和产生新微裂纹所做的功。我们的模型可以精确计算这些滞回环，并通过分离热力学贡献，区分出由塑性“摩擦”耗散的能量和由产生新损伤耗散的能量。通过在数千或数百万次循环中累积这些耗散的能量，工程师可以预测构件的疲劳寿命。
同样的理论也解释了为什么结构在循环载荷下可以非常安全，只要载荷足够小。在许多模型中，损伤演化直接与塑性变形相关。如果应力保持在材料的弹性“安全区”内，就不会发生塑性流动。如果没有塑性，也就没有新的损伤。这是“无限寿命”或“安全寿命”设计背后的基本原则，它使我们能够建造能够承受数十亿次应力循环而不会失效的桥梁、飞机和发动机部件。
损伤与塑性的耦合还可以揭示更隐蔽的失效模式。考虑一个反复加压和减压但始终处于张力下的压力容器。人们可能认为，由于载荷从未反向，变形最终应会稳定下来。然而，如果每个压力循环都引入少量损伤，它实际上会降低材料的屈服强度。在下一个循环中，为了达到相同的峰值压力（应力），材料必须多屈服一点。这会产生一个微小的、额外的永久应变增量。循环往复，这种效应会导致结构缓慢但不可逆转地伸长，这种现象被称为“棘轮效应”。这种由损伤驱动的渐进式变形，与突然过载一样肯定会导致失效，而我们的耦合模型可以完美地捕捉和预测这种行为。
损伤力学最引人注目的应用是预测材料的最终灾难性失效。该理论的一个关键见解是，失效不仅仅是一个事件，而是一个失稳的过程。当损伤使材料软化时，它会达到一个临界点。一个稍弱的小区域变形更大，这反过来又导致其损伤增长更快，使其变得更弱。一个恶性反馈循环随之而来，曾经遍布整个材料的变形迅速“局部化”到一个窄带内的剧烈应变中。这就是裂纹的诞生。我们的理论使我们能够分析变形体的稳定性，并预测这种局部化发生的确切条件——临界应力和应变——甚至可以预测由此产生的剪切带形成的角度。
这导致了一个深刻的难题。如果软化是材料在数学点上的一个属性，那么一个没有体积的裂纹应该需要零能量来形成。这显然是无稽之谈；破坏东西是需要能量的！为了解决这个悖论，我们的模型必须与断裂的物理现实联系起来。工程师引入了断裂能的概念，这是一个基本的材料属性（如密度或熔点），它量化了创建单位面积新裂纹表面所需的能量。为了确保我们的计算机模拟尊重这一物理定律，我们必须采用一种正则化技术。其中最优雅的一种是“裂纹带”模型，它规定我们本构模型的参数（特别是软化速率）必须与我们模拟中的“像素”（有限元）大小相关联。这确保了无论我们如何划分虚拟构件的网格，模拟一次完整断裂所消耗的总能量始终等于物理测量的断裂能。这是连接连续介质方程的抽象世界与物理裂纹的离散现实的一个巧妙技巧。
对更强大、更优雅的断裂描述的追求，引导研究人员在其他看似不相关的物理领域之间架起了桥梁。其中最有成果的联系之一是与相变理论的联系。这催生了断裂的“相场”模型。
与其将损伤看作是材料的一种属性，不如想象它是一种弥漫在材料中的、独特的、连续的物质或“场”，就像雾气一样。在未损伤的区域，场的值为零。当材料断裂时，场的值增长，一个完全形成的裂纹是场值为一的区域。这个“损伤场”有其自身的能量和控制方程，描述了它如何扩散和局部化。这种方法借用自描述水如何变成冰或油水如何分离的理论，通过在理论中直接构建一个物理长度尺度，优雅地解决了网格依赖性问题。它能够实现对复杂、分支裂纹网络惊人逼真的模拟，而这用传统方法是极其困难的。
这种联系有力地提醒我们物理学固有的统一性。描述雪花形成的相同数学结构可以被调整来描述钢梁的断裂。这证明了抽象的力量，也为我们的旅程画上了一个圆满的句号，表明损伤和塑性的原理不是一个孤立的主题，而是宏伟科学画卷中至关重要的一部分。