玻尔百科您是否曾想过，一个回形针在反复弯折几次后为什么会断裂？它并非在第一次弯折时就失效。恰恰相反，每一次循环都会造成无形的内部损伤，这些损伤不断累积，直到材料的完整性丧失。这种在最终断裂前很久就发生的渐进退化过程，被称为损伤。预测这一过程的科学被称为损伤建模，这是一个确保从桥梁到飞机等一切事物安全性和可靠性的关键领域。然而，捕捉这一复杂现象是一项重大挑战；简单的模型往往无法预测真实世界的行为，导致不准确甚至危险的结论。本文对损伤建模进行了全面的介绍。首先，“原理与机制”部分将深入探讨其核心理论，探索损伤在数学上是如何定义的、从简单的标量概念到先进的各向异性模型和非局部模型的演变，以及控制其增长的热力学定律。接着，“应用与跨学科联系”部分将展示这些模型的巨大实用价值，不仅在疲劳和蠕变等传统工程问题中，还在生物学和环境科学等出人意料的多样化领域中，揭示了损伤累积的普适性。
想象一下，你手里拿着一个回形针。你来回弯折它。起初，它只是改变了形状；这是一种我们称之为塑性的永久变形。此时回形针依然坚固。但继续弯折它，你会感觉它在变弱，某种意义上变得更硬，但也更脆了。最终，在最后一次弯折时，它“啪”地一声断了。从最初的弯折到最后的断裂，中间发生了什么？材料不仅仅是在变形，它内部正在死亡。这种内部退化，即微观孔洞和裂纹的累积，在材料断裂前早已损害了其完整性，我们称之为损伤。对这种“死亡”过程进行建模，就是损伤建模的艺术与科学。
要理解损伤，关键在于将其与它的“近亲”——塑性——区分开来。当材料发生塑性变形时，其原子会永久性地移动位置，就像一群人挤过一扇门——他们最终到达了新的位置，但这群人仍然存在。材料有了新的形状，但其固有的弹性特性，如刚度，基本保持不变。如果你敲击它，它仍然会发出同样清脆的声音。
损伤则不同。它是材料完整性的真实丧失。这就像人群中出现了空洞。材料承载载荷的能力下降了。它的弹性刚度降低了。如果你敲击一个受损的物体，声音会很沉闷，回响含糊不清。在我们的模型中，关键的区别在于：塑性引入了永久应变，但保持弹性刚度不变，而损伤直接降低了材料的弹性刚度。
我们如何用数学来捕捉这个想法呢？最简单、最直观的方法是创造一个数字。我们称之为，一个标量损伤变量，其取值范围从（代表原始无损材料）到（代表完全断裂的材料）。这个数字如何影响材料的行为？科学家J. Lemaitre提出了一个绝妙的想法——应变等效原理。该原理认为，在特定应力下，受损材料中观察到的应变，与未受损材料在更高的“有效”应力下所产生的应变相同。
可以这样理解：如果一个杆件的横截面有10%的面积被微孔洞占据（），那么剩余的90%的材料必须承担全部载荷。这部分完好材料所“感受”到的应力实际上比我们施加的名义应力要高。最简单的表达方式是 。遵循这个逻辑，受损材料的刚度，比如它的杨氏模量，就变成了，其中是初始的无损模量。这个简单的法则抓住了刚度损失的本质。
这个标量损伤模型非常简单。但在科学中，只有当简单能反映现实时，它才是一种美德。让我们来检验一下我们的模型。想象一下拉伸一张材料薄板。它在拉伸方向上伸长，但在垂直方向上也会变薄一些。这种侧向收缩与前向伸长之比是一个基本属性，称为泊松比，。
我们基于应变等效原理的简单模型，对泊松比随损伤增加会发生什么变化作何预测？当我们进行数学推导时，会发现一个惊人的结果：受损材料的表观泊松比与初始无损值完全相同。该模型预测，无论损伤累积多少，。
这是一个清晰、明确的预测。现在我们进入实验室。我们取一块混凝土或金属合金，对其施加拉力，并仔细测量其在微裂纹形成和扩展过程中的尺寸变化。我们经常发现，表观泊松比确实发生了变化！这个模型，在其优雅的简洁性中，忽略了损伤本质的某些基本方面。这不是失败，而是一个发现！它告诉我们，我们最初的假设——即损伤是一个简单的、无方向性的量，可以用一个数字来表示——是不完整的。
为什么简单的模型会失败？想一想损伤本身。如果你拉伸一种材料，微小的裂纹很可能垂直于拉伸方向形成。材料在那个方向上受拉时会变得弱得多，但其在横向的刚度受到的影响可能小得多。损伤通常不是各向同性的（即在所有方向上都相同）；它是各向异性的（即具有方向性）。
为了捕捉这一点，一个单一的数字是不够的。我们需要一个更复杂的数学对象，比如一个损伤张量，它可以描述在不同方向上存在多少损伤。例如，在一个双轴试验中，我们用不相等的力在两个方向上拉伸一块薄板，我们观察到刚度在应力较大的方向上退化得更厉害。标量损伤模型无法解释这一点，因为它会预测两个方向上刚度损失相等。然而，一个各向异性的张量损伤模型，可以自然地捕捉到如果应力大于，则损伤分量大于的情况，这与实验观察完全吻合。
构建这样一个模型的一种特别优美的方法是使用谱分解这个数学工具。我们可以将任何应变状态分解为其主拉伸及其方向的集合。由于裂纹和孔洞主要是由拉伸而非压缩引起的，我们可以设计一个模型，其中损伤只在拉伸主应变下才降低材料的性能，而保持其在压缩下的响应不变。这是一个绝佳的例子，说明了物理直觉（拉伸导致损伤）如何能被转化为优雅而强大的数学。
损伤是如何增长的？它并非随机发生。这里有规则，而这些规则受物理学最深刻的定律之一——热力学第二定律的支配。第二定律告诉我们，自然界中的过程倾向于是不可逆的；你不能让炒熟的鸡蛋变回生鸡蛋。损伤是一个耗散过程，它将有用的机械能转化为无用的热量。因此，损伤只能累积或保持不变；它无法自行愈合。在数学上，这意味着损伤变化率必须为非负：。这就是不可逆性条件。
损伤的增长必须有足够的“推力”。这个推力是一种热力学力，通常称为能量释放率，记为。它代表了如果损伤稍有增加，将释放到材料中的能量。损伤演化就像一个棘轮：
这个由所谓的Karush-Kuhn-Tucker（KKT）条件支配的框架，确保了损伤是不可逆的，并且只在材料受到足够载荷时发生。它提供了区分简单瞬时失效（如Tsai-Wu准则中的情况）和导致循环载荷下疲劳失效的损伤逐渐累积的机制。
现在让我们尝试在计算机模拟中使用这些优雅的模型。我们取一个虚拟杆件，拉伸它，并观察损伤的演化。一件奇怪而灾难性的事情发生了。随着材料软化，所有的应变和损伤都集中在计算机网格所能允许的最窄区域：一排单元。如果我们为了得到更精确的答案而细化网格，这个局部化区域会变得更薄。在无限精细网格的极限下，损伤区域的宽度为零，而破坏杆件所需的总能量悖论性地降至零！这被称为病态网格依赖性。我们的模拟给出的答案取决于我们选择的网格，这在物理上是荒谬的。
哪里出错了？我们的模型假设一个点的材料行为与其邻近点无关。这就是连续介质假设，在面对软化时，它失效了。真实的材料具有微观结构——晶粒、纤维、骨料。一个点的失效过程受到其周围一个小体积内发生情况的影响。
解决方法是通过教会我们的模型关于邻近点的信息来修正它。我们可以使用非局部模型，其中一个点的损伤驱动变量是其周围小区域内应变的平均值。或者，我们可以使用梯度增强模型，它会对损伤场的急剧变化施加惩罚。这两种方法都引入了一个新的基本参数：内禀长度尺度，。这个长度不是一个数值技巧；它是一种材料属性，代表了微观结构特征或断裂过程区的大小。它为失效带强制设定了最小宽度，从而对问题进行正则化，并使模拟结果独立于网格。材料本身告诉了模拟它必须如何失效。这个长度尺度是局部性失效的标志，是关于材料如何失效的深刻真理。
微观世界与我们宏观模型之间的这种联系是最后也是最关键的一环。我们的唯象模型，及其诸如的退化函数，可能看起来很抽象。但我们可以通过比较它们在小损伤极限下的预测与微观力学（研究包含显式孔洞和裂纹材料的学科）的结论，将它们与现实联系起来。通过将我们的损伤定律的初始斜率与从底层物理学推导出的斜率相匹配，我们可以为我们的连续介质模型奠定基础，确保它们不仅有效，而且忠实地反映了材料深处发生的复杂过程。从一个简单的数字到跨越尺度的深刻联系，损伤建模揭示了事物如何土崩瓦解的那个美丽、复杂而最终可预测的故事。
在我们完成了对损伤基本原理和机制的探索之后，你可能会有一种……所以呢？的感觉。我们有这些优雅的方程，这些演化的变量，但它们有什么用处？它们解决了什么问题？这是一个合理的问题，而且我认为，你会发现答案相当美妙。损伤建模的故事不仅仅是防止桥梁坍塌或飞机失事——尽管它确实如此！它也是一个关于科学思想非凡统一性的故事，其中源自研究金属断裂的思想可以阐明我们自身DNA的奥秘，甚至指导我们如何思考地球的健康。
让我们从一个熟悉的场景开始。你拿一个金属回形针，弯折它。你再把它弯回去。你又弯了一次。你凭直觉知道，你不能永远这样做下去。每一次弯折都会造成一点“损伤”，一种无形的伤害，不断累积，直到突然之间，回形针断了。最简单的思考方式就是计数。如果你知道一个全新的回形针能承受（比如说）次弯折才会失效，那么每次弯折必定消耗了其“寿命”的。如果你在一个角度下弯折次，然后在另一个角度下弯折次，著名的Palmgren-Miner法则表明，当损伤分数之和达到1时，就会发生失效。
这种线性的损伤计数是一个非常好的起点，但自然界一如既往地要更调皮一些。假设你让一个构件经受几次非常大的载荷循环，然后是许多次较轻的载荷循环。然后，你拿一个完全相同的构件，反转顺序：先是轻载荷，然后是重载荷。根据我们简单的计数规则，总寿命应该是相同的。顺序不应该有影响。但在现实世界中，它往往有影响！初始的重载荷可以产生压缩残余应力或改变微观结构，从而实际上使材料对后续的轻载荷更具抵抗力。这种“序列效应”揭示了损伤不仅仅是一个抽象的会计数字；它是一个物理过程，与材料的历史紧密相连。为了捕捉这一点，我们需要超越简单的计数，转向那些损伤本身会影响未来损伤如何累积的模型，例如通过一个非线性关系，如，或者基于每个加载循环中耗散的能量的模型。
这才是真正乐趣的开始。我们开始将损伤不视为整个物体的全局属性，而是一个场，一个随空间和时间变化的量，就像温度或压力一样。这就是连续介质损伤力学（CDM）的核心。CDM中最优雅的思想是“有效应力”的概念。想象一下，随着材料受损，它内部充满了微观的孔洞和裂纹。施加在横截面上的总力现在由一个更小的、未受损的面积来承载。为了保持平衡，这个剩余“有效”区域上的应力必须更高。我们可以用优美的简洁性来写出这一点：如果是整个区域的平均应力，那么未受损部分的有效应力是。随着损伤从0增长到接近1，有效应力急剧上升，这反过来又加速了进一步的损伤。这是一个反馈循环，一个失控过程，完美地描述了失效的最终灾难性阶段。
有了这个强大的思想，我们可以建立起惊人准确的预测模型。在高性能计算领域，工程师们将这些损伤定律嵌入到有限元模拟中。在桥梁或发动机涡轮机模型的每一个微小的计算单元（一个“高斯点”）上，计算机运行一个小算法：它计算应变，检查应变是否高到足以引起新的损伤，如果是，它就更新损伤变量，降低材料的刚度，并计算新的应力。这个“返回映射算法”是将我们的连续理论转化为具体工程预测的主力军。
世界充满了复杂的材料，我们的模型也必须同样丰富。想想现代复合材料，比如飞机上使用的碳纤维。它们不是均匀的材料；它们有强壮的纤维嵌入在较弱的基体中。失效不是一个单一事件。基体可能首先在一个方向上开裂，而纤维保持完好。一个真正具有预测能力的损伤模型必须足够智能，能够区分这两者。先进的复合材料模型正是这样做的，它们使用像Hashin失效模式这样的准则来识别是纤维还是基体失效，然后有选择地仅在适当的方向上降低材料的刚度。损伤变量不再是一个简单的标量，而变成了一个张量，捕捉了材料伤口的定向性。
当我们考虑恶劣环境时，比如喷气发动机或发电厂的内部，故事变得更加错综复杂。在高温下，材料不仅仅是疲劳；它们还会“蠕变”——在恒定载荷下缓慢而永久地变形，就像冰川从山上流下。蠕变和疲劳会以毁灭性的方式相互作用。例如，焊接会产生一个具有细晶微观结构的“热影响区”，该区域特别容易受到蠕变损伤。先进的模型可以通过将粘塑性（蠕变）变形与一个不仅对应力水平敏感，而且对应力状态（例如高静水拉应力，它会拉开孔洞）甚至材料晶粒尺寸敏感的损伤演化定律耦合起来，来捕捉这一点。为确保这些复杂模型在物理上是合理的，它们通常建立在热力学的严格基础上，从一个自由能势出发，并确保所有过程要么储存能量，要么将其作为热量耗散，但绝不会无中生有地创造能量。
但损伤从何而来？它始于微观尺度——一个断裂的原子键，一个位错堆积，一个微观孔洞。力学领域的一大挑战是弥合这个微观世界与我们观察到的宏观行为之间的鸿沟。这就是多尺度建模的领域。在这里，我们可以模拟材料的一个微小的、“代表性体积单元”（RVE），其微观结构一应俱全。我们可以将无数微裂纹的产生建模为RVE内的一个“弥散”的连续损伤场，或者我们可以明确地插入描述表面分离的“内聚界面”。通过求解这个微小RVE的力学问题，我们可以计算出宏观材料点的有效应力-应变-损伤行为。这种方法优美地展示了材料复杂的宏观软化和失效是如何从其微小组成部分的集体行为中涌现出来的。这里出现了一个迷人的微妙之处：如果处理不当，软化模型可能导致非物理的结果，即所有损伤都局部化到一条无限细的线上。为了解决这个问题，建模者必须引入一个特征“长度尺度”，这证明了损伤本质上是一个非局部现象。
到目前为止，我们谈论的都是工程结构。但“损伤累积”的概念远比这更为普适。让我们跨越到生物学的世界。想象一下，我们想存储大量信息——比如世界上所有的书籍——达数千年之久。一个不可思议的想法是将其编码在合成DNA中。但它能持续多久？像任何结构一样，DNA也会受到损伤，在这种情况下是来自宇宙辐射的损伤。我们可以用与材料疲劳惊人相似的方式来对此问题建模！我们可以将宇宙射线的撞击视为随机、独立的事件——一个泊松过程。每种类型的核苷酸，嘌呤或嘧啶，都有不同的“损伤截面”，就像对疲劳有不同的敏感性一样。通过对随时间变化的辐射通量进行建模，我们可以计算出每个核苷酸预期接收到的“撞击”平均次数，并由此得出它完好无损地存活下来的概率。通过对整个序列的这些概率求和，我们可以预测我们无价档案的半衰期。用于预测钢梁寿命的数学工具，同样可以用来预测编码生命本身的分子寿命。
让我们再做最后一次、甚至更大的跨越。我们可以谈论整个生态系统或人类社会的“损伤”吗？在生命周期评估（LCA）领域，这正是其目标所在。LCA分析产品从摇篮到坟墓的环境影响。一个LCA首先要清点所有的排放物和资源使用（例如，公斤级的，立方米的水）。这些是“中点”指标。为了使它们更有意义，像ReCiPe这样的方法试图将它们转化为“终点”指标，即对受保护领域实际危害或损伤的度量：人类健康、生态系统完整性和资源可用性。对人类健康的损伤可以用“伤残调整生命年”（DALYs）来衡量，这是从公共卫生中借来的单位。对生态系统的损伤可以用“物种年”来衡量，代表物种随时间的丧失。当然，将一公斤转换成一小部分DALY是一个复杂的过程，充满了建模、不确定性，以及至关重要的、关于不同影响相对重要性的价值判断。但概念上的联系是清晰的：我们再次在模拟一个导致系统功能退化的事件链。
从一个回形针到一个星球，损伤的概念提供了一个强大的镜头来观察世界。它教会我们，失效很少是突发事件，而是一个随时间书写的故事。通过学习阅读和理解这个故事，我们不仅能够建造更安全、更有弹性的东西，而且还能更深刻地欣赏到那些支配着系统完整性的复杂而相互关联的过程，无论是生命系统还是非生命系统。