玻尔百科当材料受到外力作用时，它们的响应方式十分复杂，并最终决定了材料的性能和寿命。两种基本但截然不同的现象主导着材料在载荷下的行为：塑性，即形状的永久改变；以及损伤，即削弱材料的渐进式内部退化。虽然我们可以直观地理解这些概念——比如弯曲回形针与它折断的区别——但它们之间错综复杂的相互作用，为试图预测材料失效的科学家和工程师带来了巨大挑战。本文旨在通过提出一个损伤与塑性耦合的统一理论框架来应对这一挑战。我们将首先在“原理与机制”一章中探究其核心物理机理，利用有效应力和热力学一致性等概念，厘清如何在实验上区分这些过程并对其进行数学建模。随后，“应用与跨学科联系”一章将展示该理论卓越的预测能力，演示其在预测机械零件的疲劳寿命、模拟灾难性材料失效等众多领域的应用。
想象一下，你拿一个金属回形针来回弯折。起初，它只是弯曲，并保持新的形状。如果松手，它不会完全弹回；它产生了一些永久变形。这就是塑性，一种形状的改变。但如果你继续弯折，你会感觉到它越来越软，最终折断。在最终断裂之前，金属内部在微观层面已经开始撕裂，布满了微小的孔洞和裂纹。这种内部退化，这种完整性的丧失，就是损伤。塑性是关于形状的改变；损伤是关于结构的破坏。这两种现象常常交织在一起，是我们故事的主角。我们的任务是理解它们各自独特的特性、它们之间秘密的联系，以及支配它们行为的普适规律。
作为科学家，我们如何区分这两个过程？它们都发生在载荷作用下，并且都对材料的失效有贡献。关键在于设计一种巧妙的探究方法。让我们把材料样本放入一台机器中，缓慢地拉伸它，同时精确地测量力（应力）和伸长量（应变）。所得的应力-应变曲线包含了所有的线索。
对于一个原始无损的材料，这条曲线的初始部分是一条直线。这是弹性区域，材料在此表现得像一个完美的弹簧。如果从这个区域卸载，它会恢复到原始形状，你输入的能量也会被释放。但如果你拉伸得足够远，它就开始屈服。此时，如果你卸载至零应力，你会发现应变并未归零，而是保留了一段永久的伸长。这个残余应变就是塑性无可置疑的印记。它是微观晶体平面相互滑移，从而永久改变材料形状的宏观证据。
那么，损伤又该如何识别呢？损伤更为隐蔽。它是一个内在的敌人，从内部削弱材料。它的标志不是永久应变，而是刚度的损失。假设我们加载样本直至产生了一定的损伤。现在，如果我们对其进行小幅卸载，然后在一个小范围内重新加载（一个“微循环”），我们就在探测它当前的弹性特性。我们会发现这条卸载-再加载线的斜率比原始的、无损状态下的弹性斜率要小。材料变得更“松软”了。初始刚度，即杨氏模量 ，已经降低到一个新的、受损的数值 。我们可以用一个简单的标量变量 来量化损伤（其中 代表原始无损状态， 代表完全失效），其关系式优雅地表示为：
这为我们提供了一个区分这两种“元凶”的绝佳策略。在试验的任何时刻，我们可以：
通过使用这些简单的卸载循环，我们可以在一次实验中独立地追踪塑性和损伤的演化过程，从而厘清它们耦合在一起的效应。
要从观察走向预测性理论，我们必须使用物理学的语言。在材料的理论描述中，核心角色是亥姆霍兹自由能，记作 。你可以把它想象成弹性储存在材料中的能量，就像一个被压缩弹簧中的势能。宇宙的一条基本法则，即热力学第二定律，规定任何真实过程要么守恒此能量，要么耗散它（通常以热的形式）。能量绝不能凭空产生。
为了描述材料的状态，我们需要数学变量。总变形由应变张量 捕捉。但我们现在知道还存在隐藏的变量，即内变量。我们引入塑性应变张量 来追踪永久的形状改变，并引入标量损伤变量 来追踪完整性的丧失。
总应变可以分解为两部分：可恢复的弹性部分 和永久的塑性部分 。
一个绝妙的洞见在于，储存的自由能 仅依赖于应变的弹性部分，因为塑性部分代表了一种永久的、已耗散的构型。但损伤是如何影响的呢？最简单也最强大的思想，即所谓的Lemaitre 模型的核心，是损伤仅仅降低了材料储存能量的能力。如果原始无损材料的弹性应变能为 ，那么对于一个截面只有 部分仍然完好的受损材料，它能储存的能量为：
这个简单的假设带来了一个深远的推论。我们实际测量的应力，即柯西应力 ，也是一个“稀释”版本，它对应于材料中仍保持完好的骨架所真实承受的应力。我们称后一个应力为有效应力 。它们之间的关系非常直观：
这就是著名的应变等效原理（一个密切相关的公式被称为应力等效，但就我们的目的而言，它们导向同一个主导定律）。这是一个概念上的突破。它意味着我们可以用与无损材料相同的规则来构建所有关于屈服和塑性流动的复杂理论，只要我们使用有效应力 而非真实应力 。我们已经将损伤的效应（即 因子）与材料固有的、现在在“有效”世界中运作的塑性响应分离开来。
我们知道了什么是塑性和损伤，并且有了一个描述它们对应力和能量影响的框架。但它们在何时以及如何发生？是什么决定了材料何时永久弯曲，或何时形成新的微裂纹？
在物理学中，变化是由力驱动的。对于我们的内变量来说，这些不是你能用手感觉到的机械力，而是从能量公式中自然产生的抽象的热力学力。塑性的驱动力与有效应力有关。损伤的驱动力是一个被称为损伤能量释放率的量，记为 。而这个力 是什么呢？它正是材料自身的弹性应变能密度。这非常直观：材料中储存的弹性应变能越多，通过产生新的开裂表面来释放能量的“压力”就越大，而这正是损伤增长的定义。
然而，驱动力并非总是足够的。桌上的一本书受到向下的引力，但它没有移动，因为桌子给它向上的支持力。只有当力超过某个阈值时，它才会移动。类似地，材料对其状态的改变也有一种内在的抵抗力。这由一个加载函数（或屈服函数） 来描述。只要热力学力的组合足够小以至于 ，材料就处于一个“安全”的弹性状态。不可逆的事情不会发生。但是，当载荷增加，使得这些力达到一个临界组合，即 时，材料便处在了临界点。它达到了其弹性极限。
接下来发生的事情由力学中最优雅的数学结构之一所支配，这套规则被称为Karush-Kuhn-Tucker (KKT) 条件。这些规则源于将热力学第二定律与最优化原理相结合，支配着塑性和损伤的演化。简单来说，它们规定：
同样的数学框架，同样的 KKT 条件，既能描述金属的延性流动，又能描述陶瓷的脆性开裂，这一事实证明了支配材料行为的物理定律深刻的统一性与美感。
到目前为止，我们一直使用一个单一的数字，一个标量 ，来描述损伤。这假设材料在所有方向上都以相同的方式弱化——这个过程被称为各向同性损伤。对于某些情况，比如受拉伸的延性金属中大致球形的孔洞的生长，这是一个合理的近似。
但是，一块木头呢？或者一片来自现代飞机的碳纤维复合材料板呢？这些材料天生具有方向性，即各向异性。它们沿纹理或纤维方向的强度和刚度与横跨该方向的性质大相径庭。一个垂直于纤维方向的裂纹远比一个平行于纤维方向的裂纹更具破坏性。一个单一的数字 完全不足以描述这一现实。
为了捕捉这种方向性的弱化，我们需要一个更复杂的数学对象。我们引入一个损伤张量 ，而不是标量，它可以被看作一个 的矩阵。这个张量有其自身的主方向和主值，分别代表最大损伤的方向及其大小。例如，在一个纤维间形成裂纹的复合材料中，损伤张量会显示在垂直于纤维的方向上损伤值很大，而沿着纤维方向的损伤值则很小。这使我们能够模拟真实世界高性能材料中发生的复杂、各向异性的刚度退化。该理论的优雅之处得以保留：驱动这种张量损伤的热力学力本身也是一个张量，从而在内变量与其共轭力之间保持了完美的对称性。这种在保持基本原理不变的情况下，丰富数学描述以捕捉日益复杂的物理现实的能力，是一个强大而优美的科学理论的标志。
现在我们已经熟悉了损伤与塑性之间错综复杂的舞蹈，你可能会好奇：“这套复杂的机制有什么用？”一个物理理论，无论多么优雅，其价值在于它描述、预测并最终帮助我们改造世界的能力。在本章中，我们将踏上一段旅程，看看我们刚刚学到的损伤与塑性耦合理论，如何不仅仅是一种抽象的数学游戏，而是一个在广泛的科学和工程学科中使用的强大而通用的工具。我们将看到，我们刚刚学到的同样的基本原理，可以解释材料在反复应力作用下的缓慢死亡，可以解释一块金属中剪切带的剧烈形成，甚至可以解释陨石撞击造成的陨石坑的灾难性诞生。这正是物理学固有的美：它的统一力量。
我们的理论最直接、或许也是最重要的应用在于结构工程领域。当我们建造一座桥梁、一架飞机或一栋摩天大楼时，最重要的问题是：“它能撑住吗？”要回答这个问题，我们需要超越入门物理课程中的简单弹性模型。我们不仅需要知道材料何时开始弯曲，还需要知道它何时开始永久变形，以及最重要的是，它何时会断裂。
想象一下拉伸一根简单的金属棒。我们的理论可以预测它从最初的轻微拉伸到最终灾难性断裂的完整生命历程。起初，材料像弹簧一样弹性拉伸。随着拉力的增加，塑性流动开始，材料开始硬化。但与此同时，另一个更险恶的过程也在进行：损伤正在累积。微观孔洞正在形核和生长。在某个时刻，这种损伤的弱化效应将开始超过塑性硬化的强化效应。继续拉伸杆件所需的应力将达到一个峰值，然后开始下降。这种“软化”是材料的绝唱，是一个预示着完全失效即将来临的明确信号。通过结合有效应力、塑性硬化和损伤演化定律的原理，我们可以构建出材料完整的应力-应变曲线，从而精确预测这个峰值载荷以及随后的失效路径。这不仅仅是一个学术练习；它是现代安全分析和抗失效设计的基础。
当然，如果没有正确的数字输入，模型是无用的。我们如何确定一种真实材料（比如一种新型钢材）中控制损伤和塑性的参数呢？我们不能只是查阅资料；我们必须通过测量来获得。这是实验力学的领域，它既是一门科学，也是一门艺术。在这里，我们化身侦探，设计巧妙的实验，让材料揭示其秘密。
一个关键的挑战是损伤和塑性是深度交织在一起的。如果我们只是拉伸一个样本，这两个过程会同时发生。我们如何分离它们的影响？在这里，我们可以利用我们理论中的一个深刻见解：损伤的演化对应力状态的性质高度敏感。特别是一个称为*应力三轴度*的量——粗略地说，就是材料从各个方向被拉伸的程度，就像带缺口杆件的中心那样——起着主导作用。高的正三轴度，如在简单的拉伸试验中，会促进微孔洞的生长，从而加速损伤。
相反，纯剪切状态（你可以通过扭转一个薄壁管来想象）的三轴度为零。在这种状态下，塑性流动可以发生，而孔洞生长却非常少。一个聪明的实验者可以利用这一点。他们可以对材料进行剪切试验，以在近乎无损伤的条件下分离和表征其塑性硬化行为。然后，有了这些知识，他们可以进行一个损伤显著的单轴拉伸试验。通过在拉伸试验中测量材料刚度的退化，他们可以精确地量化损伤的演化过程，，并校准模型的损伤部分。这种系统性的解耦是理论与实验协同作用的一个绝佳例子，使我们能够建立稳健的、可预测的材料行为模型。
为什么有些材料，如铜，在断裂前会大幅拉伸，而其他材料却几乎没有预兆地突然断裂？答案在于塑性硬化和损伤引起的软化之间的微妙竞争。材料在断裂前能够大量变形的能力——其延性或韧性——是其最宝贵的属性之一，尤其对于飞机机身或压力容器等安全关键应用。
我们的理论揭示了塑性硬化是这个故事中的英雄。当材料变形时，应力倾向于集中在微观缺陷或孔洞的尖端。如果材料没有硬化能力，所有后续的塑性流动都将集中在这些点上，迅速导致裂纹形成。然而，作功硬化使得应力集中处的材料变得更强，迫使塑性变形扩散到周围的材料中。这个过程“钝化”了初生的裂纹，并将应变分散到更大的体积中，从而延迟了孔洞连接形成宏观断裂的时刻。因此，具有高硬化能力的材料具有更好的损伤容限；它们可以优雅地适应变形而不会屈服于灾难性失效。
但这个故事也有其黑暗的一面。当平衡被打破时会发生什么？当由损伤快速增长驱动的材料软化，压倒了作功硬化的稳定效应时，材料的切线刚度可能变为负值。在这个临界点，平衡的控制方程失去了一个称为椭圆性的数学性质。其物理后果是一场灾难：材料再也无法以平滑、均匀的方式变形。取而代之的是，所有变形突然集中到一个极窄的区域，称为剪切带。这种现象被称为应变局部化，标志着从延性撕裂到突然的、类似脆性的剪切失效的转变。理解这种不稳定性对于预测金属成形操作的极限和分析某些类型的高速断裂至关重要。
并非所有的失效都源于单一的、剧烈的事件。我们世界中的许多部件，从汽车的轮轴到喷气发动机的涡轮叶片，都承受着数百万次微小、重复的加载循环。在这种无情的侵袭下，材料可能会在远低于其单次拉伸所能承受的应力水平下失效。这个隐蔽的过程被称为疲劳。
损伤和塑性的框架为理解这一现象提供了一个绝佳的视角。我们可以区分两种主要的疲劳机制。在高周疲劳 (HCF)中，施加的载荷足够小，以至于大部分材料都呈弹性响应。损伤累积极其缓慢，可能源于晶界或夹杂物处的微观塑性滑移。在这种机制下，失效循环次数与施加的应力幅值关联性最好。
然而，在低周疲劳 (LCF)中，载荷足够大，以至于在每个循环中都会引起显著的、广泛的塑性变形。材料的应力-应变路径形成一个闭合的环路，称为滞回环。该环路所包围的面积代表了在每个循环中转化为热量的塑性功——它是材料耗散的能量。我们的理论告诉我们，这种耗散能是正在发生的塑性损伤的直接度量。在 LCF 机制下，决定材料寿命的不是应力幅值，而是塑性应变幅值或每循环的耗散能。对于任何必须在循环载荷下承受长久使用寿命的部件的设计和分析而言，这种区分是根本性的。
在现代世界，我们不能总是通过建造和破坏实物来检验其安全性。我们必须对它们进行模拟。计算力学使我们能够在计算机内部建立一个“虚拟实验室”，在那里我们可以在各种条件下测试设计。我们讨论过的损伤和塑性模型正是这些模拟的核心与灵魂。
将这些理论在一个计算机程序中实现，通常使用有限元法 (FEM)，是一次从物理到算法的非凡旅程。我们理论的结构，特别是应变等效原理，指导了数值更新过程的设计，即所谓的“回映算法”。在一个小的时间步内，首先假设一切都是弹性的，计算出一个“试探”状态。如果这个试探状态违反了屈服条件，一个“校正”步骤会将状态拉回到屈服面上，并计算出塑性流动。然后，后续步骤根据新状态更新损伤。
然而，一个深刻而微妙的陷阱等待着粗心的程序员。对软化材料模型进行直接的“局部”实现——即某一点的状态仅取决于该点的历史——会导致对有限元网格中单元尺寸的病态依赖。随着网格的细化，模拟的断裂区变得无限薄，而产生失效所耗散的能量会虚假地趋近于零！这意味着模拟结果不具有物理意义。
这个危机的解决是现代计算力学的伟大胜利之一。问题在于局部模型缺少一个基本的材料属性：内禀长度尺度。为了使模型正则化并获得与网格无关的结果，我们必须引入这个长度尺度。这可以通过高级方法实现，如非局部或梯度增强损伤模型，或通过实用方法，如“裂纹带模型”，该模型根据单元尺寸来调整材料的软化行为。这些技术确保了模拟中每单位断裂面积耗散一个恒定的、与网格无关的能量，这个量被称为断裂能 ()。这些正则化方法将一个简单的本构模型转变为一个用于模拟失效的稳健、可预测的工具。
我们框架的多功能性确实令人惊叹。同样的核心思想可以被应用于描述截然不同尺度上的现象。
考虑两种材料之间的界面，比如粘合接头或复合材料的层间。我们可以将这个界面建模为一个拥有自身本构律的无限薄层。在这里，“损伤”可以代表界面的渐进脱粘和分离，而“塑性”可以代表摩擦滑移。这类*内聚区模型*是现代断裂力学的基石，用于模拟复杂材料和结构中裂纹的萌生和扩展。
最后，让我们将目光从微小和缓慢转向宏大和灾难性的事物。当一颗以每秒数公里速度飞行的陨石撞击行星时，会发生什么？这是一个几乎无法想象的剧烈事件，涉及巨大的变形、极端的压力和极高的应变率。然而，同样的基本原理也适用。我们可以建立一个计算模型，将大变形运动学与一个包含 Johnson-Cook 塑性模型（用于应变率相关性）、渐进损伤以及——至关重要的——热力耦合的材料定律结合起来。撞击过程中巨大的塑性功转化为热量，导致岩石温度飙升。这种热软化削弱了材料，如果温度达到熔点，就会发生相变。通过将所有这些物理要素整合到一个大规模的显式动力学模拟中，我们可以模拟整个成坑过程：初始的冲击波、物质的挖掘以及陨石坑的最终形状。
从不起眼的拉伸试验到宇宙级的碰撞，损伤与塑性的耦合理论为理解材料如何变形和失效提供了一个统一而强大的框架。它证明了在世界纷繁复杂的表象之下，蕴含着一种非凡而优美的简单性。