损伤演化律：原理、机理及应用

玻尔百科

核心要点

热力学第二定律为损伤演化提供了根本依据，它规定了变形过程中耗散的能量是材料退化的驱动力。
有效应力的概念解释了微观空洞如何减少材料的承载面积，巧妙地将微观尺度的损伤与宏观尺度的力学行为联系起来。
损伤演化律经过实验数据标定后，对于预测复杂工程结构的寿命和失效模式的计算机仿真至关重要。
损伤力学原理具有高度的通用性，构成了一种描述失效的通用语言，适用于生物力学、复合材料科学乃至原子物理学等不同领域。

引言

材料很少在瞬间失效。从混凝土桥梁产生微裂纹，到金属部件在高温下缓慢拉伸，失效是一个渐进的内部退化过程，我们称之为损伤。理解并预测这一过程是工程学和材料科学中最关键的挑战之一。其难点在于，我们需要超越简单的观察，构建一种预测性理论——一种损伤演化律——它不仅是描述性的，而且深深植根于物理学的基本定律。本文通过对连续介质损伤力学理论的全面介绍来应对这一挑战。文章首先探讨其核心原理和机理，展示热力学第二定律如何支配材料退化的不可逆过程。随后，通过综述其在从预测关键部件寿命到模拟生物力学和复合材料等不同领域的失效过程中的多样化应用，展示了该理论巨大的实际应用价值。我们的探索之旅将从揭示材料为何以及如何断裂的热力学核心开始。

原理与机理

如果你拉一根橡皮筋，它会伸长。再用力拉，它就会断裂。但在它断裂前的瞬间发生了什么？或者想想一座混凝土桥梁，在发生任何灾难性破坏之前，它会产生一张几乎看不见的微小裂纹网。这些日常观察指向一个深刻的真理：材料很少一次性完全失效。它们经历一个渐进的内部退化过程，即微裂纹、空洞和断裂化学键的累积。我们将这个过程称为损伤。我们的任务是为这个过程建立一种理论——一种损伤演化律——不仅仅是通过描述我们所看到的现象，而是要将其根植于物理学最基本的原理。对于任何涉及变化、能量和不可逆性的过程，我们的基石是热力学第二定律。

问题的核心：能量与耗散

想象一下拉伸一块金属。你对它做功，大部分功以弹性势能的形式储存起来，就像压缩弹簧一样。当你松手时，金属会弹回，释放这些储存的能量。但如果你拉伸过度，它会开始产生内部损伤。你所做的功不再完全可恢复。其中一部分会不可挽回地损失掉，转化为热量以及形成新表面——即微观裂纹的表面——所需的能量。这种损失的能量称为耗散。热力学第二定律坚称，对于任何自发过程，这种耗散必须是正的。物体会损坏，但不会自发地“修复”。这是刻在物质核心的时间之箭。

为了精确地表达这一思想，物理学家使用了亥姆霍兹自由能的概念，用希腊字母 $\psi$ (psi) 表示。这是系统内能中可以“自由”转化为功的部分。单位体积内对材料所做的功的速率由应力张量 $\boldsymbol{\sigma}$ 与应变率张量 $\dot{\boldsymbol{\varepsilon}}$ 的缩并给出。储存能量的变化率是 $\dot{\psi}$ 。克劳修斯-杜亥姆不等式，即热力学第二定律的局部表述，宣称耗散率 $\mathcal{D}$ 必须为非负：

\mathcal{D} = \boldsymbol{\sigma}:\dot{\boldsymbol{\varepsilon}} - \dot{\psi} \ge 0

这个简单而强大的不等式是我们的指南针。为了使用它，我们需要描述材料的状态。我们引入一个单一而优美的思想：一个标量损伤变量 $D$ 。我们定义 $D$ 为一个介于 0 和 1 之间的数，其中 $D=0$ 代表原始、未损伤的材料，而 $D=1$ 代表完全失效、不具备任何承载能力的材料。

损伤如何影响自由能 $\psi$ ？最直接的假设是，损伤降低了材料储存弹性能的能力。如果未损伤材料储存的能量密度为 $\psi_0$ ，我们可以提出，损伤材料储存的能量为：

\psi = (1-D) \psi_0

对于一个简单的一维杆，弹性能为 $\psi_0 = \frac{1}{2} E \varepsilon^2$ ，其中 $E$ 是杨氏模量， $\varepsilon$ 是应变。因此，自由能变为 $\psi(\varepsilon, D) = \frac{1}{2}(1-D)E\varepsilon^2$ 。这种优美的形式抓住了退化过程的本质。

毁灭的驱动力

现在，让我们看看我们的指南针——耗散不等式——将我们引向何方。通过对 $\dot{\psi}$ 应用链式法则并将其代入不等式，经过一番数学整理后，会发生一个显著的分离，正如在和等问题的详细分析中所示。耗散被揭示为两个项的简单乘积：

\mathcal{D} = Y \dot{D} \ge 0

在这里， $\dot{D}$ 是损伤累积的速率。因为我们知道损伤是不可逆的（物体不会自我修复），所以必须有 $\dot{D} \ge 0$ 。这迫使另一项 $Y$ 也必须是正的。但是 $Y$ 是什么？当我们进行数学推导时，我们发现 $Y$ 是自由能对损伤的负偏导数，即 $Y = -\frac{\partial\psi}{\partial D}$ 。对于我们的简单模型，这得到：

Y = \frac{1}{2} E \varepsilon^2 = \psi_0

这是一个惊人的结果。我们称之为损伤能量释放率的量 $Y$ ，恰好是材料在未损伤状态下所储存的弹性能密度。它是驱动损伤演化的热力学力。你试图通过应变在材料中储存的能量越多，通过产生新的微裂纹来释放该能量的热力学“压力”就越大。储存的能量本身成为了材料自我毁灭的引擎。

裂纹的剖析：有效应力

让我们从一个更形象、更力学的角度来探讨同一个问题。损伤变量 $D$ 在物理上代表什么？想象一下在显微镜下观察一根钢筋的横截面。当它受损时，并非整个材料都变弱了，而是微观空洞和裂纹出现并扩展。名义上的横截面积 $A$ 保持不变，但实际承载载荷的固体材料面积——有效面积 $\tilde{A}$ ——在缩小。一个简单而有力的对 $D$ 的物理解释是，它是因这些缺陷而损失的面积分数： $\tilde{A} = (1-D)A$ 。

我们从外部测量的应力，即柯西应力，是作用力 $F$ 分布在整个面积上的结果： $\sigma = F/A$ 。但剩余固体部分中的原子并不知道空洞的存在；它们只感觉到集中在它们身上的力。因此，它们所经受的应力，即有效应力 $\tilde{\sigma}$ ，要高得多： $\tilde{\sigma} = F/\tilde{A}$ 。

片刻的代数运算将这两个观点以一个优美而简单的公式联系起来：

\tilde{\sigma} = \frac{F}{(1-D)A} = \frac{\sigma}{1-D}

这种关系是应变等效原理的基础。它假设，在柯西应力 $\sigma$ 作用下，受损材料的力学响应与一个假想的、未受损的材料在更高的有效应力 $\tilde{\sigma}$ 作用下的响应相同。从某种意义上说，材料被“欺骗”了。其剩余的、完整的部分的行为与它们一贯的表现完全相同，但它们承受的应力比外部世界所感知的要强烈得多。这一由 Kachanov 和 Lemaitre 等先驱倡导的概念，巧妙地弥合了微观空洞的现实与我们观察到的宏观行为之间的鸿沟。

失效定律：何时开始？速度多快？

我们现在有了一个驱动力 ( $Y$ ) 和一个清晰的物理图像（有效应力），但我们仍然需要一个预测性定律。我们需要回答两个问题：损伤何时开始增长，以及一旦开始增长，其速度有多快？这就是损伤演化律的作用。

损伤萌生

材料可以承受一定量的应变而没有任何永久性损伤。就像一个正在被注水的玻璃杯，直到驱动力 $Y$ 的“水平”达到一个临界阈值，损伤才开始。我们可以通过在热力学力空间中定义一个损伤面来将其形式化，通常采用简单的形式，如 $f(Y, \kappa) = Y - \kappa \le 0$ ，其中 $\kappa$ 代表材料当前的损伤抵抗能力。只要 $Y$ 小于 $\kappa$ ，材料就表现为弹性行为，不会发生新的损伤。只有当状态处于损伤面上 ( $f=0$ ) 并且被进一步推动时，损伤才会增长。这种被称为库恩-塔克条件的逻辑结构确保了不可逆性，并且它也是用来描述塑性屈服开始的那个强大的框架。

损伤演化

一旦满足该准则， $D$ 的增长速度有多快？答案取决于材料的性质。

对于像金属这样的延性材料，损伤与塑性流动密不可分。构成损伤的微观空洞随着固态金属基体的流动而被扩大和拉伸。因此，损伤率 $\dot{D}$ 必须与塑性变形速率耦合，通常用累积塑性应变率 $\dot{p}$ 来衡量。Jean Lemaitre 提出的著名模型提出了一个幂律关系，它巧妙地将热力学驱动力与塑性机理结合起来：

\dot{D} = \left(\frac{Y}{S}\right)^{s} \dot{p}

在这里， $S$ 是一个代表材料固有损伤抗性的材料参数， $s$ 是一个控制响应非线性程度的指数。该定律优美地指出，只有在存在塑性流动（ $\dot{p} > 0$ ）时损伤才会增长，并且其增长速率会被能量释放率 $Y$ 放大。

对于其他现象，例如蠕变——材料在高温下缓慢且随时间变化的变形——则需要另一种定律。Kachanov 和 Rabotnov 最初的唯象模型提出，损伤率直接取决于时间和应力水平。使用我们的有效应力概念表示的一种典型形式是：

\dot{D} = B \tilde{\sigma}^m = B \left(\frac{\sigma}{1-D}\right)^m

这类定律对于预测高温环境下部件（如喷气发动机涡轮或发电厂锅炉）的断裂时间非常实用。

细节中的魔鬼：超越最简模型

我们简单的标量模型 $\psi = (1-D)\psi_0$ 功能强大，但现实总是更为丰富。它的局限性迫使我们进行更深入的思考，并揭示新的物理现象。

拉伸与压缩： 如果你推一根混凝土柱，它会非常坚固。但如果你拉它，它会轻易断裂。原因在于，在压缩作用下，其内部的微裂纹被压闭，裂纹表面相互摩擦，使得柱子能够传递载荷。然而，我们简单的 $(1-D)$ 因子对受拉还是受压是“视而不见”的。为了捕捉这种单边效应，我们必须采用更复杂的方法。例如，Mazars' model 巧妙地定义了一个“等效拉伸应变” $\varepsilon^*$ ，它仅由应变张量的正（拉伸）部分构成。如果一种材料处于纯静水压缩状态，其所有主应变为负值，因此 $\varepsilon^*$ 为零，根据该模型，不会有损伤演化。这是一个根据具体观察到的物理行为来调整通用热力学框架的绝佳范例。

方向性损伤： 想象一种由强纤维（如碳纤维复合材料）增强的材料。如果产生裂纹，它们很可能会沿着与纤维垂直的方向扩展，从而在该方向上对材料的削弱程度远大于沿纤维方向。这种损伤是各向异性的。一个单一的数字 $D$ 无法描述这种方向偏好。为此，我们必须将损伤变量从标量提升为张量，这是一种更复杂的、携带方向信息的数学对象。

连续介质的失效： 当材料开始软化，即其所能承受的应力随着应变的增加而开始减小时，可能会出现一个最微妙和深刻的挑战。在这个阶段，变形倾向于局部化到一个无限薄的带中。当我们试图用计算机通过有限元法模拟这一过程时，会出现一种病态现象：失效带的宽度缩小到单个计算单元的尺寸。当我们使用更精细的网格以获得更准确的答案时，预测的失效变得越来越脆，总耗散能趋于零。模拟结果变得具有病态的网格依赖性。

这个数值计算上的噩梦是一个深层次数学问题的症状：问题的控制方程失去了一种称为椭圆性的性质。纯粹的局部模型——即某一点的应力仅取决于该点的应变——缺少了一个关键的物理要素：一个内在的长度尺度。为了解决这个问题，我们必须在模型中引入一个长度尺度。我们可以通过非局部理论来实现这一点，即某一点的损伤演化取决于其周围一个小邻域内应变的平均值。这种“涂抹”效应防止了局部化现象收缩成零厚度的线，恢复了问题的适定性，并产生了随着网格细化而收敛的、符合物理现实的结果。这展示了物理学、数学和计算之间美妙的相互作用，说明一个领域的崩溃如何能揭示另一个领域中缺失的原理。

归根结底，对物体如何损坏的研究绝非一件病态之事。它是一段旅程，带领我们从裂纹和断裂的现实世界，深入到能量和熵的基本原理，再回到现代计算工程的复杂挑战。损伤演化理论是物理学统一性的证明，展示了如何建立、调整和完善一个单一、连贯的框架，来描述自然界中最复杂和最重要的过程之一。

应用与跨学科联系

既然我们已经熟悉了损伤演化的内部机制——即材料具有内部生命，或者说，一个缓慢走向死亡的进程——现在是时候提出最重要的问题了：这又如何？损伤变量 $D$ 这个抽象概念究竟能带我们走向何方？你会欣喜地发现，答案是：几乎无处不在。它是一把万能钥匙，能解开一系列惊人多样的问题，从预测喷气发动机的寿命到理解为何髋关节置换会失败。在本章中，我们将把我们的理论机器投入运行，看看一个材料从内部逐渐变弱这个简单的想法，如何提供一种统一的语言来描述各种形式的失效。

预测的艺术：寿命与可靠性

损伤力学最直接、最有价值的应用之一是回答工程师最迫切的问题：“它还能用多久才会坏？”我们的基础设施和技术的安全取决于我们能否自信地回答这个问题。

想象一下喷气发动机涡轮内部的一片叶片。它以每分钟数千转的速度在足以熔化钢铁的温度下旋转，同时被巨大的离心力向外拉扯。即使是最先进的合金也无法永远承受这样的折磨。在这种持续的应力和高温下，材料开始“蠕变”——它会随着时间的推移缓慢地永久拉伸。我们的损伤框架为我们提供了一个观察这个过程的窗口。在金属内部，微观空洞开始形核和生长，主要沿着晶界分布。我们可以通过建立损伤累积速率 $\frac{d\omega}{dt}$ 与*有效应力*——即尚未布满孔洞的材料所感受到的应力——的关系来模拟这个过程。

随着损伤 $\omega$ 的增加，可用于承载的横截面积减小。这意味着剩余完整材料上的应力会增加，这反过来又使损伤增长得更快。这就形成了一个恶性反馈循环。它解释了在实验中观察到的蠕变的一个关键特征：第三阶段，即变形速率开始急剧加速，迅速导致断裂。通过将损伤演化律从其原始状态（ $\omega=0$ ）积分到最终失效状态（ $\omega=1$ ），我们可以计算出总断裂时间 $t_r$ 。这为预测关键部件的工作寿命提供了一个强大的工具，从而能够在灾难发生前及时更换。

但材料不仅仅在极端的恒定载荷下失效。它们也会因疲劳而失效，这是导致从桥梁到飞机机身等各种结构失效的无声杀手。一根反复弯折的回形针最终会断裂，尽管每一次单独的弯曲都是无害的。在这里，损伤的“时钟”不是时间，而是加载循环次数 $N$ 。在每个加载和卸载循环中，微量的能量——即滞回能量，也就是应力-应变滞回环内的面积——在材料内部耗散。这种能量为微裂纹的生长提供了动力。我们可以建立一个基于循环次数的损伤演化律 $\frac{dD}{dN}$ ，由这种循环能量输入驱动。这使我们能够估算结构的疲劳寿命，这是现代机械设计的基石。

连接理论与现实：从实验室到计算机

这些数学模型虽然优美，但如果不能与真实材料联系起来，它们将是无用的奇珍异品。我们如何找到我们方程中各种材料参数的值——那些指数和系数？我们通过倾听材料本身来找到它们。

在实验室里，我们可以取一个材料样本并将其拉开，同时仔细测量应力随应变的增加。由此产生的应力-应变曲线就像是材料的“签名”。它告诉我们在载荷下的行为，包括达到其峰值强度并开始软化或失去承载能力的关键时刻。这种峰后软化行为是损伤累积的直接体现。通过将我们基于热力学的损伤模型与这条实验曲线进行拟合，我们可以精确地确定我们方程所需的材料参数。这个标定过程确保了我们抽象的模型不仅仅是一个数学游戏，而是对真实材料特性的忠实再现。

一旦经过标定，这些损伤定律就成为强大的计算机模拟（最著名的是有限元法，FEM）的核心。在数字世界中，我们可以建立汽车、建筑或喷气发动机的虚拟原型，并对其施加模拟的运行载荷。通过引入损伤演化，我们可以在屏幕上观察微观裂纹的萌生和扩展，从而准确地看到结构最薄弱的地方。现代模型可以将损伤与其他复杂行为（如塑性（永久变形）和粘塑性（变形依赖于加载速率））耦合起来。这使得工程师能够模拟高度动态和灾难性的事件，如车祸或鸟撞飞机机翼，并设计出以更安全、更可预测的方式失效的结构。我们之前探讨的严格的热力学框架确保了这些复杂的模拟不会违反基本的物理定律，从而赋予它们真正的预测能力。

一种通用的失效语言：拓宽视野

也许损伤演化律最美妙之处在于其纯粹的通用性。其基本概念——一个系统随着其组成部分的失效和载荷重新分配而变弱——并非金属所独有。它是一个普遍的原理，适用于惊人广泛的材料和科学学科。

考虑一个由纤维增强复合材料而非金属制成的现代飞机机翼。在这里，“损伤”不是空洞的生长，而是单根高强度纤维的断裂。当一根纤维断裂时，它所承受的载荷会转移到相邻的纤维上，使它们更容易断裂。我们可以为这个过程写一个损伤演化律，但有所不同。纤维的强度不是均匀的；它遵循一个统计分布。通过将我们的损伤框架与强度统计理论（如 Weibull 模型）相结合，我们可以建立极其复杂的模型，来预测复合材料的疲劳寿命和可靠性。

同样的想法甚至适用于我们自身。在生物力学中，连续介质损伤力学被用来理解骨组织的退化。在骨科植入物（如髋关节置换物）周围，骨骼可能会经历异常的应力集中。随着时间的推移，这可能导致假体周围骨中微裂纹的累积——即损伤——从而可能导致植入物松动。通过基于局部应变能对骨的损伤演化进行建模，研究人员可以预测高风险区域，并设计出能与人体更安全、更持久地结合的植入物。这是利用机械工程原理改善人类健康的一个深刻例子。

我们旅程的最后一站将我们从工程结构的宏大规模带到单个原子的领域，揭示了失效物理学中深刻的统一性。人们可能会想，我们的连续介质损伤变量 $D$ 是否只是一个方便的虚构。惊人的答案是否定的。利用强大的超级计算机，我们可以运行分子动力学（MD）模拟，跟踪一小块材料中每个原子的运动。我们可以亲眼看到，在应力和热振动的影响下，原子如何重新排列，形成最初的萌生缺陷。这个过程发生的速度可以用过渡态理论的基本物理学来描述。神奇之处在于，当我们对数百万个原子的这种原子行为进行平均时。其集体结果是一个宏观的损伤演化律，其数学形式与我们一直使用的唯象定律完全相同。事实证明，我们的工程模型是原子统计力学的直接结果。

从涡轮机的缓慢蠕变到桥梁的循环疲劳，从实验室中模型的标定到计算机中失效的模拟，从先进复合材料的断裂到我们自身骨骼的健康，最后，将这一切与原子的基本舞蹈联系起来——损伤演化概念提供了一个单一、连贯的叙述。它证明了一个简单的物理思想在为复杂且看似混乱的材料失效世界带来清晰度和预测能力方面的强大力量。