各向同性损伤

玻尔百科

核心要点

有效应力概念，定义为 $\tilde{\sigma} = \frac{\sigma}{1-D}$ ，是其中心思想，解释了内部损伤如何放大了作用在材料剩余部分上的应力。
应变等效原理提供了一种方法，通过简单地使用有效应力，将无损材料的本构律应用于受损材料。
从热力学的角度来看，损伤增长的驱动力是材料在无损状态下所储存的应变能。
各向同性损伤理论将“裂纹带”中耗散的总能量与材料的断裂能 ( $G_f$ ) 联系起来，从而将连续介质力学与断裂力学联系在一起。

引言

理解材料如何以及为何失效是工程与科学领域的一项核心挑战。虽然失效始于微观裂纹和孔洞，但逐一追踪它们是不切实际的。连续介质损伤力学通过在宏观尺度上对材料劣化进行建模，提供了一种强大的替代方案。这种方法引入一个被称为损伤变量的状态变量来表示内部缺陷的集体效应，从而解决了微观缺陷与可观测的结构弱化之间的鸿沟。

本文全面概述了该领域的基础模型——各向同性损伤。您将首先探索赋予该理论预测能力的核心思想。随后，您将看到这个优雅的框架如何应用于解决不同学科的实际问题。第一章“原理与机理”将介绍有效应力、应变等效原理以及模型的热力学基础等基本概念。之后，“应用与跨学科联系”一章将展示如何利用这些原理来分析结构完整性、塑性、蠕变乃至地质构造。

原理与机理

为了理解材料如何失效，我们不一定需要追踪每一个微观裂纹和孔洞。那将是一项不可能完成的任务，就像试图通过追踪每个空气分子来描述天气一样。相反，我们可以退后一步，着眼于“大局”。连续介质损伤力学为我们提供了实现这一目标的工具，其核心思想非常优雅且强大：有效应力的概念。

一个简单的图景：有效应力

想象一下你正在拉一根金属杆。在宏观层面上，它似乎是完全实心的。但当你用力拉时，在材料深处，微小的缺陷——微孔洞和微裂纹——开始出现并生长。它们就像材料结构中的微小孔洞。最直接的后果是什么？能够实际承载载荷的材料部分变小了。

假设我们杆的原始横截面积是 $A_0$ 。经过一些拉伸后，这部分面积的一部分，我们称之为 $A_d$ ，现在被这些缺陷所占据。仍在工作的剩余“完好”面积为 $A_{\text{eff}} = A_0 - A_d$ 。你施加的力 $F$ 现在集中在这个更小的有效面积上。

虽然我们作为观察者可能会将名义应力计算为 $\sigma = F/A_0$ ，但材料本身感觉不到这种平均化的应力。它在仍然承载的部分上承受着更为强烈的应力。我们称之为有效应力 $\tilde{\sigma}$ 。其定义是简单的力学： $\tilde{\sigma} = F/A_{\text{eff}}$ 。

为了连接这两个世界——我们所见的和材料所感的——我们引入一个极其简单的数字：损伤变量 $D$ 。它被定义为因损伤而损失的横截面积的比例： $D = A_d / A_0$ 。一个原始、无损的材料 $D=0$ 。一个沿横截面完全失效的材料 $D=1$ 。

有了这个定义，有效面积可以写成 $A_{\text{eff}} = A_0 - D A_0 = A_0 (1-D)$ 。现在我们可以找到我们测量的应力（ $\sigma$ ）和材料感受到的应力（ $\tilde{\sigma}$ ）之间的关系。根据简单的力平衡，无论我们从外部还是内部看，总力都是相同的：

$F = \sigma A_0 = \tilde{\sigma} A_{\text{eff}}$

代入我们对 $A_{\text{eff}}$ 的表达式，我们得到：

$\sigma A_0 = \tilde{\sigma} A_0 (1-D)$

两边除以 $A_0 (1-D)$ ，我们就得到了有效应力的基本方程：

\tilde{\sigma} = \frac{\sigma}{1-D}

这个小小的方程是问题的核心。它告诉我们，随着损伤 $D$ 从 0 增长到 1，有效应力 $\tilde{\sigma}$ 会变得比名义应力 $\sigma$ 大得多。这解释了为什么即使在缓慢增加的载荷下，材料也会突然失效：在内部，剩余连接部分的应力正在急剧飙升。这个概念源于承载面积减小的简单图景，构成了我们整个讨论的基础。

一个统一的思想：应变等效原理

有效应力的概念可以从一个简单的图景提升为一个深刻的物理原理。这就是应变等效原理 (PSE)，它是损伤力学的一个基石。它做出了一个大胆而有力的论断：

在真实应力 $\boldsymbol{\sigma}$ 下，在受损材料中测得的应变，与原始无损材料在承受有效应力 $\tilde{\boldsymbol{\sigma}}$ 时所经历的应变相同。

想一想这意味着什么。它为我们提供了一座神奇的桥梁。我们不需要为受损材料发明全新的物理定律。我们可以简单地采用原始材料众所周知的本构律（如弹性力学的胡克定律），然后用我们的“有效”应力替换“真实”应力。

对于弹性材料，无损定律是 $\boldsymbol{\varepsilon} = \mathbb{S}_0 : \boldsymbol{\sigma}$ ，其中 $\mathbb{S}_0$ 是材料的初始柔度张量（其“弹性”）。PSE 告诉我们，对于受损材料，该定律仅仅是：

$\boldsymbol{\varepsilon} = \mathbb{S}_0 : \tilde{\boldsymbol{\sigma}}$

现在我们可以使用连接两个应力的公式 $\tilde{\boldsymbol{\sigma}} = \boldsymbol{\sigma}/(1-D)$ ，并将其代入：

$\boldsymbol{\varepsilon} = \mathbb{S}_0 : \left( \frac{\boldsymbol{\sigma}}{1-D} \right) = \frac{1}{1-D} (\mathbb{S}_0 : \boldsymbol{\sigma})$

这个结果太棒了。它告诉我们，受损材料的本构律是 $\boldsymbol{\varepsilon} = \mathbb{S}(D) : \boldsymbol{\sigma}$ ，其中新的、受损的柔度张量 $\mathbb{S}(D)$ 只是原始张量按比例放大：

\mathbb{S}(D) = \frac{1}{1-D} \mathbb{S}_0

随着损伤的增加，材料变得柔度更大——更软，或更具“延展性”——这在物理上是完全合理的。PSE 为预测这种软化提供了一个极其简单的方案。

“各向同性”的含义

我们一直在讨论各向同性损伤，但这个假设到底意味着什么？它意味着损伤没有优选方向。微裂纹和微孔洞被假定为在方向和空间上随机分布，形成一个弥散、无序的网络。材料的性质在所有方向上均等退化。

这会带来一些非常具体且可检验的推论。由于柔度张量 $\mathbb{S}_0$ 只是被一个数字缩放，其基本特性保持不变。对于各向同性材料，其弹性性质可以用两个常数来描述，例如杨氏模量 $E$ 和泊松比 $\nu$ 。

我们的结果 $\mathbb{S}(D) = \mathbb{S}_0 / (1-D)$ 意味着受损材料的有效杨氏模量 $E(D)$ 变为 $E(D) = (1-D) E_0$ 。刚度直接因损伤而降低。但泊松比呢？泊松比 $\nu$ 描述了当你在一个方向上拉伸材料时，它在横向方向上变薄的程度。由于刚度在所有方向上的降低都是相同的，这个比率保持不变！该模型一个惊人的预测是，表观泊松比不随损伤而改变： $\nu(D) = \nu_0$ 。

我们可以通过一个具体的例子看到这种方向无关性的作用。考虑一个薄板，在x方向上受到应力 $\sigma_1$ 的拉伸，在y方向上受到应力 $\sigma_2$ 的拉伸。对于无损材料，应变为： $\varepsilon_1^{(0)} = \frac{1}{E_0}(\sigma_1 - \nu \sigma_2)$ $\varepsilon_2^{(0)} = \frac{1}{E_0}(\sigma_2 - \nu \sigma_1)$ $\varepsilon_3^{(0)} = -\frac{\nu}{E_0}(\sigma_1 + \sigma_2)$

现在，如果我们应用应变等效原理，受损材料中的应变变为： $\varepsilon_1 = \frac{1}{1-D} \varepsilon_1^{(0)}$ $\varepsilon_2 = \frac{1}{1-D} \varepsilon_2^{(0)}$ $\varepsilon_3 = \frac{1}{1-D} \varepsilon_3^{(0)}$

每一个应变分量都被简单地放大了相同的因子 $1/(1-D)$ 。响应被均匀地、没有任何方向偏好地缩放，这正是各向同性的本质。

更深层的基础：热力学视角

到目前为止，我们的模型是基于直观的力学论证。但它是否尊重物理学的基本定律，特别是热力学第二定律？答案是肯定的，而探索这一点揭示了更深层次的美。

在热力学中，弹性材料的状态可以用其亥姆霍兹自由能 $\psi$ 来描述，它代表了储存的弹性势能。因此，应力不仅仅是一个力学量，而是一个热力学量，由该能量导出： $\boldsymbol{\sigma} = \frac{\partial \psi}{\partial \boldsymbol{\varepsilon}}$ 。

我们的受损材料的自由能是什么？一个自然而简单的假设，通常被称为能量等效假设，是损伤仅仅降低了材料储存能量的能力。所以，我们可以将受损自由能 $\psi(\boldsymbol{\varepsilon}, D)$ 写成原始自由能 $\psi_0(\boldsymbol{\varepsilon})$ 的一部分：

\psi(\boldsymbol{\varepsilon}, D) = (1-D) \psi_0(\boldsymbol{\varepsilon}) = (1-D) \left( \frac{1}{2} \boldsymbol{\varepsilon} : \mathbb{C}_0 : \boldsymbol{\varepsilon} \right)

其中 $\mathbb{C}_0$ 是无损刚度张量。让我们看看这对压力意味着什么：

$\boldsymbol{\sigma} = \frac{\partial \psi}{\partial \boldsymbol{\varepsilon}} = (1-D) \frac{\partial \psi_0}{\partial \boldsymbol{\varepsilon}} = (1-D) (\mathbb{C}_0 : \boldsymbol{\varepsilon})$

这正是我们之前从应变等效原理推导出的应力-应变定律！这两个原理，一个从有效应力出发，另一个从能量退化出发，对于各向同性损伤得出了相同的结果。它们是同一枚硬币的两面，揭示了该理论令人满意的统一性。

但热力学给了我们更多。损伤是一个不可逆过程——裂纹不会自发愈合。这意味着它必须产生熵，或耗散能量。驱动不可逆过程的热力学“力”是通过观察自由能如何随内变量变化来找到的。对于损伤，这个力被称为损伤能量释放率 $Y$ ：

$Y = - \frac{\partial \psi}{\partial D}$

使用我们对 $\psi$ 的表达式，我们发现：

$Y = - \frac{\partial}{\partial D} \left[ (1-D) \psi_0(\boldsymbol{\varepsilon}) \right] = (-1) \cdot (-1) \cdot \psi_0(\boldsymbol{\varepsilon})$

Y = \psi_0(\boldsymbol{\varepsilon}) = \frac{1}{2} \boldsymbol{\varepsilon} : \mathbb{C}_0 : \boldsymbol{\varepsilon}

这是一个优美而深刻的结果。它表明，产生更多损伤的热力学驱动力，恰恰是材料在无损状态下所储存的应变能。我们拉伸或变形材料越多，这种储存的能量就越高，材料产生损伤以释放这种能量的热力学“动机”就越大。这是大自然寻找更低能量状态的方式，即使这意味着破坏自身。

模型的局限性：当简单性不再足够时

这个各向同性损伤模型简单、优雅且强大。但像物理学中的所有模型一样，它是对现实的一种近似，建立在一组特定的假设之上。同样重要的是要理解模型在什么情况下不适用。

首先，连续损伤变量 $D$ 的整个概念依赖于代表性体积单元 (RVE) 的存在。我们必须能够放大到一个比单个微裂纹大，但仍比整体结构小的尺度。我们是在对微观的混沌进行平均。

其次，各向同性的假设至关重要。如果我们以一种会产生定向微裂纹的方式加载材料——例如，在一个方向上强力拉伸它——损伤将是各向异性的。材料在垂直于裂纹的方向上会比平行于裂纹的方向上更弱。单个标量 $D$ 无法捕捉到这一点，需要一个更复杂的张量损伤变量。

第三，该模型忽略了单边效应。它假设无论材料是受拉还是受压，刚度都降低 $(1-D)$ 。实际上，微裂纹在压缩下可以闭合并传递载荷，使得材料在压缩下的刚度比我们简单模型预测的要高。

那么，在实际实验中，我们如何知道我们的各向同性模型是否失效呢？我们必须设计能够探测量材料方向特性的测试。

一种方法是进行仔细的力学测试，以测量完整的损伤柔度张量 $\mathbb{S}^d$ 。如果各向同性模型成立，该张量的每个分量相对于无损张量 $\mathbb{S}_0$ 都应该按相同的因子 $1/(1-D)$ 进行缩放。如果我们发现需要不同值的 $D$ 来解释不同分量的变化——例如，x方向的刚度退化比剪切刚度更多——那么模型就失效了。损伤不是各向同性的。

另一种强大的技术是使用超声波。我们可以让声波穿过材料并测量其速度。对于各向同性材料，声速在所有方向上都是相同的。我们的模型预测损伤会使声速降低一个因子 $\sqrt{1-D}$ ，但它应该保持与方向无关。如果我们进行实验并发现波速现在沿一个轴比另一个轴快，我们就找到了诱导各向异性的明确证据，我们简单的标量模型就不再足够了。

这就是科学的循环。我们基于清晰的物理原理建立一个简单、优美的理论。我们探索它的预测，并惊叹于其内在的一致性。然后，我们用实验来检验它，将其推向极限，以发现它在何处失效。正是在这个空间——在理论的成功与失败的边界上——下一个更完整的理解诞生了。

应用与跨学科联系

现在我们已经熟悉了各向同性损伤的原理，我们准备踏上一段旅程，看看这个优雅的思想如何在世界上找到它的位置。像任何优秀的科学概念一样，其真正的价值不是体现在其抽象的公式中，而是在于其解释、预测和连接我们在工程和自然界中观察到的各种现象的能力。我们将看到，“有效应力”这个简单而核心的思想——即材料的完好部分必须更努力地工作以补偿受损部分——是一把用途广泛的钥匙，解开了许多学科中的奥秘。

强度与刚度工程

让我们从最直接的应用开始：结构工程领域。想象一根简单的钢筋支撑着一个重物。我们可以通过将力 $F$ 除以横截面积 $A$ 来计算其内部的应力。但如果材料不完美呢？如果它因制造过程或先前使用而布满微观孔洞和裂纹呢？总面积 $A$ 在几何上仍然存在，但并非所有面积都能承载载荷。损伤变量 $D$ 给了我们一个衡量这个不可用面积的尺度。实际起作用的有效面积仅为 $(1-D)A$ 。

那么，材料的完好骨架所感受到的真实应力是多少？它是力 $F$ 除以这个更小的有效面积。这就是有效应力， $\tilde{\sigma} = F / ((1-D)A)$ 。通过将其与名义应力 $\sigma = F/A$ 进行比较，我们得到了我们之前见过的基石关系： $\tilde{\sigma} = \sigma / (1-D)$ 。这不仅仅是一个数学技巧；这是一个物理陈述。对于任何程度的损伤，材料的健康部分所经历的应力都高于我们名义计算所显示的应力。正是这种被放大的应力决定了材料的命运。

这个简单的思想以优美的一致性向上扩展。它只适用于拉伸杆件吗？完全不是。考虑一个被扭转、承受剪切的材料块。材料抵抗剪切的能力由其剪切模量 $G_0$ 描述。在受损材料中，同样的逻辑也适用。我们观察到的刚度被有效地降低了，因为内部结构受到了损害。表观剪切模量变为 $G(D) = G_0(1-D)$ 。如果我们考虑最一般的三维应力状态，情况也是如此。整个四阶弹性张量 $\mathbb{C}_0$ ——包含了材料刚度的所有信息——被均匀地劣化了。受损的刚度张量就是 $\mathbb{C}_d = (1-D)\mathbb{C}_0$ 。各向同性损伤中的“各向同性”正是这个意思：材料在所有方向上被同等地削弱，这是将整个刚度张量乘以一个单一数字 $(1-D)$ 的直接结果。

这对现实世界的结构具有深远的影响。在工程设计中，我们痴迷于“应力集中”——即应力被自然放大的区域，比如孔洞或尖角的周围。考虑一个带有圆孔的无限大板，承受远场拉应力 $\sigma_0$ 。经典的弹性理论告诉我们，孔洞边缘的应力可以高达远场应力的三倍。现在，如果这个板已经包含一个均匀的背景损伤水平 $D$ ，会发生什么？分析揭示了一些非凡的现象。受损的板的行为与一个承受更高远场应力 $\sigma_0 / (1-D)$ 的无损板完全相同。因此，孔洞边缘的应变不仅被几何形状放大，也被损伤放大了。峰值应变增大了 $1/(1-D)$ 倍。这告诉我们，预先存在的损伤和几何特征共同作用，造成了比任何一个单独因素都危险得多的情况，为裂纹的萌生和扩展铺平了道路。

损伤与塑性和蠕变的共舞

到目前为止，我们只谈论了弹性材料，它们会恢复到原来的形状。但许多材料，如金属，也可以发生永久变形——这种行为被称为塑性。损伤如何与此相互作用？同样，有效应力的概念是我们的指南。

塑性屈服是由材料内部的应力状态决定的。假设一个无损金属在应力达到某个值 $\sigma_Y$ 时屈服。在受损材料中，名义应力 $\sigma$ 可能很低，但完好连接部分上的有效应力 $\tilde{\sigma}$ 要高得多。当这个有效应力达到 $\sigma_Y$ 时，屈服就会开始。这意味着引起塑性变形所需的名义应力降低到 $\sigma_Y^{damaged} = (1-D)\sigma_Y$ 。损伤有效地软化了材料，使其更容易永久变形。损伤与塑性的这种耦合对于精确模拟延性结构的失效至关重要，从碰撞中的汽车底盘到地震中的钢梁。在这样的过程中，总的能量耗散是两部分之和：用于塑性流动的能量和用于产生新微裂纹（损伤）的能量。

当我们考虑高温下的材料时，损伤的影响变得更加显著，这些材料会发生蠕变——在恒定载荷下缓慢、持续的变形。一个经典的蠕变曲线显示了三个阶段：变形减缓的第一阶段，速率稳定的第二阶段，以及变形加速直至断裂的第三阶段。几十年来，第三阶段有些神秘。为什么在恒定载荷下的材料会突然开始越来越快地拉伸？

连续介质损伤力学提供了一个优美的答案。在第二阶段，材料硬化和回复之间存在动态平衡。然而，在此期间，损伤 $D(t)$ 在悄悄累积。随着 $D(t)$ 的增加，有效应力 $\tilde{\sigma} = \sigma / (1-D(t))$ 稳步上升，即使名义应力 $\sigma$ 是恒定的。最终，有效应力的增加变得如此之快，以至于它压倒了任何硬化效应，导致应变率加速。这就是第三阶段。损伤为这种失控破坏提供了引擎，这对于设计必须在极端条件下安全运行多年的喷气发动机涡轮叶片或核电站组件来说，是一个至关重要的见解。

从无穷小到有限：连接损伤与断裂

连续介质损伤力学的哲学挑战之一是，它将损伤视为“弥散”在一个体积内。然而，我们知道失效通常最终会形成一个单一的、尖锐的裂纹。我们如何将连续场的世界与离散断裂的世界连接起来？

答案在于一个强大的思想，它将损伤与断裂力学的基石——断裂能 $G_f$ 联系起来。断裂能是一种材料属性，它量化了创建单位新裂纹面积所需的能量。在计算模拟中，当材料因损伤而软化时，变形倾向于集中在一个窄带中。“裂纹带”模型提出了一个简单而有力的假设：在特征宽度为 $l_c$ 的这个带内，由损伤过程耗散的总能量必须等于材料的断裂能 $G_f$ 。

这使我们能够校准我们的连续介质损伤模型。通过对整个损伤过程中单位体积耗散的能量 $\int Y dD$ 进行积分，并乘以带宽 $l_c$ ，我们可以强制这个量等于实验测得的 $G_f$ 。这个过程为我们提供了一种确定损伤演化律参数的方法，确保我们的模拟在断裂过程中耗散正确的能量，而不管我们计算机模型中单元的尺寸如何。这是一项杰出的理论工程，它将连续和离散的失效观点结合在一起。

自然界中的损伤：岩土力学

损伤力学的原理不仅限于金属和人造材料。它们对于理解我们脚下的地球同样至关重要。岩石、土壤和其他地质材料本质上是多孔的，并且常常是断裂的。考虑一个深埋地下的饱和土。其强度不是由上覆岩石的总压力（总应力）决定的，而是由*有效应力*——由固体颗粒骨架承载的那部分总应力——决定的。孔隙中水的压力 $p_f$ 抵消了总应力，有效地“松开”了颗粒，降低了它们的摩擦阻力。

这就是著名的多孔介质力学中的有效应力原理。但是当岩石骨架本身受损时会发生什么？受损的骨架柔度更大。根据多孔弹性理论，一个柔度更大的骨架会导致一个更大的 Biot 系数 $\alpha$ 。这个系数决定了孔隙压力屏蔽骨架免受总应力影响的有效程度。损伤 $D$ 的增加会增加 $\alpha$ ，这反过来又会在给定的总应力和孔隙压力下降低有效平均应力 $p'$ 。本质上，损伤给予了双重打击：它直接削弱了固体骨架，同时使其更容易受到孔隙压力削弱效应的影响。这种综合理解对于从石油和天然气开采、地热能到边坡和地基稳定性的各种应用都至关重要。

前沿：各向异性与热力学之美

与任何强大的模型一样，理解其局限性同样重要。各向同性损伤模型假设微裂纹和孔洞是随机分布的，在所有方向上同等地削弱材料。但这并不总是正确的。想象一下压缩一个多孔岩石。它很可能会产生优先垂直于压缩方向排列的微裂纹。材料变得各向异性——在一个方向上比另一个方向更弱。

一个简单的各向同性模型无法捕捉到这一点。例如，当对受损岩石施加围压时，微裂纹会闭合，材料会恢复部分刚度。这个恢复过程通常是各向异性的。一个使用单个标量 $D$ 的各向同性模型，与追踪不同方向平面上损伤的更复杂的各向异性模型（如微平面模型）相比，可能会高估或低估这种刚度恢复。这个局限性并没有否定各向同性模型；相反，它定义了其适用范围，并为处理各向异性占主导地位的问题指明了更先进理论的方向。

最后，这些模型的结构中隐藏着一种形式上的美。当一个损伤模型的构建方式与热力学定律一致时——具体来说，通过从自由能势推导出应力和其它内力——一个奇妙的数学性质常常会出现：将应变的无穷小变化与应力的无穷小变化联系起来的切线刚度算子，变得对称。这种对称性不仅仅是一个数学上的奇特现象；它是互易定理（如弹性力学中的 Betti 定理）的基础。这意味着A点的力对B点位移的影响与B点的相同力对A点位移的影响是相同的。对于热力学上一致且“关联的”（意味着损伤在与能量释放相关的方向上增长）损伤模型，这种对称性成立。对于不符合这些条件的模型，它则不成立。这在热力学的抽象定律、材料的实际行为以及我们计算算法的稳健性之间提供了深刻的联系。

从拉伸一根杆件到涡轮叶片的加速蠕变，从混凝土的断裂到山坡的稳定性，源于有效应力这个简单思想的各向同性损伤概念，提供了一个统一的框架。它证明了一个简单的物理思想能够为一个广阔而复杂的世界带来清晰度和预测能力的力量。