玻尔百科预测材料如何以及何时断裂是现代工程学的基石，对于设计从安全车辆到耐用基础设施的各种产品至关重要。然而，用于完成此任务的传统工具——经典连续介质力学，却存在一个关键缺陷。当材料开始软化和失效时，标准的计算机模拟可能会产生物理上荒谬的结果，即计算出的断裂能取决于模拟的任意设置，而不是材料本身。这个被称为“病态网格敏感性”的悖论，使得这些模型在实际预测中并不可靠。
本文通过介绍非局部损伤模型的理论来解决这个根本性问题。我们将探讨这些先进的框架如何克服经典理论的局限性。在第一章“原理与机制”中，我们将深入研究局部模型的数学悖论，并揭示引入一个特征性的“内禀长度尺度”如何提供解决方案。我们将检验两种主要方法——积分型模型和梯度增强模型——它们教会材料点与其邻近点进行“交流”。随后，在“应用与跨学科联系”中，我们将看到这些理论概念如何成为强大的工程工具，实现客观的模拟，解释如尺寸效应等真实世界的现象，甚至在宏观结构与原子世界之间架起桥梁。读完本文，您将理解为什么非局部性不仅仅是一种数学修正，更是精确模拟材料失效的一个基本原则。
想象一下：你拿一根金属棒，拉它的两端。起初，它像一根非常硬的橡皮筋一样弹性伸长。但再用力拉，微小的裂纹开始在内部形成和增长。材料开始“软化”，失去了承载更多载荷的能力。最终，金属棒断裂。一个简单的问题出现了：断裂它需要多少能量？常识和热力学第一定律都要求这是一个确定且可测量的能量值。
现在，让我们尝试使用自19世纪以来一直沿用的经典连续介质力学定律，在计算机上模拟这个过程。我们将金属棒表示为点的集合，并在每个点上应用我们的应力和应变规则。然而，在这里，我们陷入了一场灾难。当模拟的材料开始软化时，模型发现，从能量上看最有利的做法是将所有失效集中到尽可能小的区域。在数学连续体中，最小的可能区域宽度为零。我们的计算机模拟试图忠实于这个数学原理，预测失效将局部化为仅一个单元宽度的裂纹。如果我们细化模拟网格，使单元更小，裂纹就会变得更窄。在无限细化网格的极限情况下，裂纹的宽度为零。
悖论就在这里：宽度为零的裂纹体积也为零。如果所有的断裂能都在这个不存在的体积内耗散，那么断裂整个金属棒所需的总能量就是……零！ 这是一场灾难。我们模拟的结果现在完全取决于我们选择的网格的任意尺寸。这种问题被称为病态网格敏感性。我们优美的数学模型给出了一个物理上毫无意义的答案。对此问题的数学诊断是椭圆性丧失；描述材料行为的控制方程在软化区域其性质发生了根本性改变，从而允许这些物理上不可能的、无限尖锐的解出现。
那么，问题出在哪里？缺陷深藏于经典连续介质理论的一个基本假设中：局部性。局部性假设空间中某一点的材料状态（如其应力）仅取决于该完全相同点上发生的情况（如其应变）。每个材料点都是一座孤岛，对其邻居的状态一无所知。
但真实的材料并非如此。在显微镜下观察，你会看到一个丰富而复杂的世界。金属由晶粒构成，混凝土是沙子、砾石和水泥浆的混合物，而骨骼则是纤维和矿物质的奇妙复合物。 在这个微观世界里，点与点之间在不断地相互作用。一个晶体中形成的微裂纹会向其邻居发送应力波。一个晶粒的滑移会推挤下一个晶粒。这里存在一个“通信范围”。
因此，要纠正我们这个充满悖论的模型，就必须放弃“孤点”假设。我们必须教会我们的材料点相互交谈。我们通过在理论中引入一个新的基本参数来实现这一点：一个内禀长度尺度，用希腊字母 表示。这个长度代表了材料微观结构的特征尺寸——平均晶粒大小、碎石的直径，或增强纤维之间的典型距离。正是在这个长度尺度上，简单的局部性图像失效，与邻近区域的相互作用变得至关重要。
一旦我们接受了材料点必须进行通信，问题就变成了：如何通信？物理学家和工程师们已经发展出两种主要且非常直观的方法来模拟这种非局部性。
第一种方法是积分型非局部模型。这个想法非常民主。我们不再让单一点的应变触发损伤，而是让其周围邻域内的应变进行“投票”来做决定。我们将非局部等效应变（我们称之为 ）定义为点 附近局部应变 的加权平均值。
在这里， 是一个权重函数或核函数。它像一个影响函数：附近的点 在确定 处的非局部应变时被赋予更大的“投票权”，而远处的点的影响则会衰减。该函数具有显著权重的特征距离就是我们的内禀长度 。
分母是这个谜题的关键部分；它是一个归一化因子。它确保了如果应变恰好在整个材料上完全均匀，非局部平均值就简单地等于那个均匀应变。我们这个更复杂的新模型正确地再现了简单的旧情况——这是一个至关重要的合理性检验。
这种平均化的效果是深远的。它像一个低通滤波器，平滑了应变场。如果在一个微小的点上出现了一个危险的高应变峰值，平均过程会削弱这个峰值，将其影响扩散开来。由于损伤现在是由这个平滑后的场驱动的，失效在物理上就不可能坍缩成一条零宽度的线。它被迫在一个有限的“过程区”内扩展，该区域的宽度与 相关。
第二种方法是梯度增强模型。这种方法采取了不同的哲学路径，但达到了相似的目的地。我们不是明确地进行平均，而是假设自然界厌恶尖锐的边缘。我们重新构建物理学，规定材料在损伤状态上产生急剧的空间变化时，必须“支付”一个能量代价。
我们通过修改材料储存的内能（其亥姆霍兹自由能 ）的方程来做到这一点。我们增加了一个取决于损伤梯度平方 的项。
在这里， 是一个材料模量， 是杨氏模量， 再次是我们的内禀长度尺度。 这个新项使得损伤从一点到另一点发生巨大而快速的变化在能量上变得昂贵。
当你有了这样一个能量项时会发生什么？当你通过数学（使用变分法）来寻找平衡状态时，一个神奇的算子出现在了损伤的控制方程中：拉普拉斯算子 。 这是物理学许多领域中的老朋友。它与控制热量在固体中扩散的数学算子相同。扩散的作用是什么？是平滑事物！一个热点不会一直保持热点状态；它通过加热周围环境来冷却。完全相同地，这个拉普拉斯项迫使任何高损伤区域“分散载荷”，防止损伤变得无限局部化。
乍一看，“邻里投票”（积分型）和“同伴压力”（梯度型）模型似乎截然不同。一个是有限区域上的积分；另一个是基于局部梯度的微分方程。但物理学常常有这样美妙而隐藏的统一性。
事实证明，作为许多梯度模型核心的亥姆霍兹型微分方程 ，可以被形式化地求解。在一个足够长的一维杆上，它的解正是一个积分平均，其中权重核函数 就是简单的指数函数 。 这揭示了梯度模型并非一个陌生的概念，而是积分模型的近亲——在某些情况下，数学上是等价的。它们是同一核心原理的两种表现形式：通过使材料点意识到其在特征长度 范围内的邻居，来强制实现平滑性。
所以，我们已经治愈了这个数学疾病。我们的模型现在给出了网格无关的结果：计算出的断裂一个物体的能量是一个有限的、可预测的材料属性，无论我们把模拟网格做得多细。 这是理论一致性的一个胜利。但任何物理理论的真正考验不仅在于解决悖论，还在于预测一些关于世界的新而真实的东西。
在这方面，非局部模型表现得非常出色。它们自然地解释了在混凝土、岩石、冰和先进陶瓷等材料中观察到的著名的尺寸效应。
想象你有一系列几何形状相同但尺寸不同的带缺口混凝土梁——一个足够小，可以放在实验室工作台上；另一个足够大，可以作为桥梁的一部分。一种朴素的材料强度方法会认为它们都在相同的名义应力下失效。但实验表明这是错误的。更大的梁按比例来说更弱，并且以一种更为脆性、灾难性的方式失效。
非局部理论完美地解释了这一点。结构的行为是其特征尺寸 （如梁的高度）与材料的内禀长度 之间竞争的结果。
对于非常大的结构 (): 决定断裂过程区大小的内禀长度尺度与梁的尺寸相比可以忽略不计。失效行为就像理想尖锐裂纹的扩展。失效时的名义强度 与尺寸的平方根倒数成比例，即 。这正是经典线性弹性断裂力学的预测。
对于非常小的结构 ( 与 的量级相当): 断裂过程区的大小现在与整个结构的尺寸相当。失效不再是一个尖锐的裂纹，而是一个更弥散的退化过程。失效由材料的体强度决定，名义强度 变得几乎恒定，与试件尺寸无关。
该理论提供了一条优美的、统一的曲线，将小尺寸物体的基于强度的行为与大尺寸结构的脆性、断裂力学行为联系起来。通过测试不同尺寸的试件，工程师可以将其数据拟合到这条曲线上，并测量出特定材料的 值。最初为解决悖论而发明的数学技巧，如今已成为一种强大的、可预测的工程工具，将“通信范围”的抽象概念转化为一个有形的数字，帮助我们设计更安全、更可靠的结构。
既然我们已经探讨了非局部模型的抽象原理，你可能会问一个很合理的问题：“我为什么要关心这些？”这是一个很好的问题。答案是，这些看似深奥的数学思想，正是让我们从仅仅描述材料失效转向准确预测它的工具。它们在数值计算和物理现实之间架起了一座桥梁。没有非局部性的概念，我们最强大的关于物体如何断裂的计算机模拟将产生——坦率地说——毫无意义的结果。
让我们踏上一段旅程，看看这个深刻的思想——某一点发生的事情取决于其邻域——如何在工程、材料科学领域开启更深的理解，甚至将工程结构的世界与原子的量子领域联系起来。
想象一下，你是一名工程师，任务是设计一辆新车。你希望确保它在碰撞中是安全的。你求助于你最强大的工具：有限元法（FEM），这是一种将复杂结构分解为数百万个微小、简单的部分，并在其上求解物理定律的方法。你构建了一个精美的汽车数字模型，并模拟了一次高速撞击。当金属部件开始变形和撕裂时，你的模拟中发生了奇怪的事情。损伤没有像现实中撕裂的边缘那样扩展，而是集中在一条无限细的线上，宽度恰好是一个单元。如果你用更精细、更小单元的网格再次运行模拟，这条线会变得更细，而断裂该部件所需的总力也改变了！结果不取决于材料的物理特性，而取决于你——工程师——如何绘制你的网格。这就是我们讨论过的“病态网格依赖性”，它使得模拟对于预测毫无用处。
这正是非局部模型被发明出来要解决的问题。回想一下我们那个中心有高应变的一维杆。一个局部模型会看到这个尖锐的峰值，并立即在该单一点上“断裂”。但非局部模型做的事情更为微妙和智能。为了决定是否在那个中心点引发损伤，它会观察其邻居。它执行一个加权平均——一种“邻里投票”——对应变场进行处理。一个点上的极高应变，如果被低应变的邻居包围，就会在平均过程中被“平滑掉”。那个尖锐的、不符合物理规律的应变峰值被展平成一个更宽、更平缓的凸起。
结果呢？模拟的损伤现在局部化在一个有限宽度的带内，而这个宽度不是由任意的网格尺寸决定的，而是由我们模型中内建的材料内禀长度尺度 决定的。现在，如果你细化网格，模拟会收敛到一个单一的、有物理意义的结果。数值噩梦被驯服了。这不仅仅是一个理论练习；对于应用复杂的材料模型，如Johnson-Cook法则，这是至关重要的。该法则是模拟车祸和弹道冲击中高应变率塑性和失效的行业标准。没有正则化，这些复杂的模型将建立在沙地之上。
起初，内禀长度 似乎只是一个解决数值问题的巧妙数学技巧。但故事远不止于此。这个长度尺度原来是一种真实的、物理的材料属性，通过将其纳入模型，我们的模型不仅能保持稳定，更能真正地进行预测。
考虑准脆性材料（如混凝土、岩石或飞机中使用的先进复合材料）中的“尺寸效应”。如果你取大小两块混凝土，一块小，一块巨大，并测试它们的强度，你会发现大块的按比例来说比小块的弱。经典的裂纹理论——线性弹性断裂力学无法解释这一点。但非局部模型可以。
想象一下测试一个带有圆孔的复合材料板，这是飞机机身中常见的情况。经典理论告诉我们，孔边缘的应力集中是一个固定因子，与孔的大小无关。但实验显示了不同的情况：孔越大，板材变得越弱。非局部模型完美地解决了这个悖论。模型理解失效并非由无穷小一点的应力决定，而是由一个尺寸为 的区域内的平均应力状态决定。对于小孔，高应力区域相对于 较小，因此非局部平均显著“钝化”了峰值应力，使板材看起来更强。对于非常大的孔，应力场相对于 变化缓慢，平均效应较小，其行为接近经典预测。
这意味着我们可以对不同孔径的试样进行几次实验，用数据来校准唯一的未知参数 ，然后用我们的非局部模型自信地预测任何尺寸孔或缺口的真实构件的强度！长度尺度 是连续介质用以表示“断裂过程区”（FPZ）的方式，这是扩展裂纹尖端微裂纹和能量耗散的区域。非局部模型甚至能够捕捉像R曲线这样的现象，即由于过程区的发展，扩展裂纹所需的能量随着裂纹开始增长而增加。即使在使用内聚区模型对复合材料分层进行建模时，引入非局部性对于获得客观结果也至关重要。长度尺度不再是一个“修正”；它是材料韧性的一个基本特征。
非局部概念的美妙之处在于其普适性。局部化问题不仅出现在我们模拟桥梁或飞机机翼时；它可以在任何使用连续介质描述的尺度上出现。这一点在计算多尺度建模领域变得惊人地清晰。
想象一下，试图预测一种由碳纤维嵌入聚合物基体中的新型复合材料的性能。整体行为取决于纤维和基体之间复杂的相互作用。一种称为FE²（有限元平方）的强大技术通过在模拟中运行模拟来解决这个问题。在你对大型构件的“宏观”模拟中的每一点，模型都会调用一个独立的“微观”模拟，该模拟针对材料内部结构的一个小的代表性体积单元（RVE）。这个微观模拟计算出局部的应力-应变响应，然后将其传递回宏观模拟。
但转折点在这里。如果你的RVE中的聚合物基体开始失效会怎样？如果你对基体使用一个简单的、局部的损伤模型，你的微观模拟将遭受病态网格依赖性并变得不适定！你可能已经猜到了，解决方案是使用一个正则化模型——如非局部、梯度或相场公式——在RVE内部以物理上有意义的方式捕捉微裂纹。
但故事并未就此结束！在运行了你适定的微观模拟之后，你可能会发现复合材料的整体、均匀化响应也表现出软化。如果你接着将这个均匀化定律用于对整个构件的标准、局部的“宏观”模拟中，那么宏观模拟现在将变得不适定。问题在上一级尺度再次出现！解决方案是相同的：宏观模型也必须是一个非局部或高阶理论。这种类似分形的级联效应揭示了，当涉及失效时，非局部性不是一个临时的修正，而是在尺度间一致地传递信息所必需的基本原则。
我们已经看到了非局部模型的几种形式。我们讨论了在邻域上对场进行平均的积分模型，也提到了对场的急剧变化进行惩罚的梯度或相场模型。这些似乎是非常不同的哲学。一种在精神上是全局的，另一种是局部的但对导数敏感。然而，该领域最优雅的发现之一是它们是深度相关的。对于在空间中缓慢变化的现象，泰勒展开表明，一个积分非局部模型，在主导近似下，等同于一个梯度模型[@problem_id:2667996, @problem_id:2700808]。它们是用不同的数学语言描述相同的基本物理：即一个材料点对其周围环境有“感知”。
这把我们带到了最后一个，也是最深刻的联系。这个内禀材料长度尺度 最终来自哪里？它来自原子。任何材料的行为最终都由其原子键的性质决定。有一种方法可以将宏观长度尺度 与这些基本量联系起来。一个从平衡失效时的弹性能与创建新表面所需能量得出的、极具洞察力的关系表明：
让我们欣赏一下这其中的美妙。这里， 是杨氏模量，衡量原子键的刚度。 是表面能，是打破那些键并创建一个新裂纹表面所需的能量。而 是理论内聚强度，是拉开一个完美、无瑕疵晶体所需的应力。这些都是根植于物质的原子和量子力学性质的量。
于是，我们回到了原点。一个始于使工程模拟正常工作的实际需求的理念，最终被证明是一个深刻的物理原理。非局部长度尺度 是原子键的量子世界、微裂纹的混乱介观尺度以及工程设计的宏观世界之间的握手。它是让我们从原子层面建立起对失效的理解，从而能够设计出塑造我们现代世界的复杂、可靠结构的连续性之线。而这确实是一件了不起的事情。