玻尔百科计算机模拟已经彻底改变了工程学，使我们能够为从桥梁到飞机机翼的一切事物构建和测试数字孪生体。然而，在传统的连续介质力学理论中，潜藏着一个根本性的缺陷。当这些模型试图模拟材料失效——这恰恰是它们最常被需要解决的问题——它们可能会产生物理上荒谬的结果：模拟网格越精细，材料看起来就越脆弱。这场被称为“病态网格敏感性”的危机，源于一个错误的假设，即材料在某一点的行为与其周围环境无关。本文旨在通过介绍一种更复杂、物理基础更坚实的非局部损伤模型，来填补这一关键的知识空白。第一章“原理与机制”将诊断传统模型中“局部性的弊病”，并详述非局部公式所提供的“治愈之道”，探讨积分方法和梯度增强方法这两种不同但相关的途径。随后的“应用与跨学科联系”一章将展示这些理论如何转化为工程师和物理学家的强大工具，解释如尺寸效应等深奥现象，并揭示其跨越科学学科的深刻联系。
想象一下，你是一位正在设计桥梁的工程师。你想知道它在极端应力下的行为，直到可能形成裂纹的那一刻。你求助于功能强大的计算机，为你的桥梁构建了一个精美细致的数字孪生体——一个由数百万个微小互联点组成的模型，我们称之为有限元网格。你运行模拟，它向你展示了桥梁最薄弱的地方。“太好了，”你想，“现在让我们得到一个更精确的答案。”你细化网格，用数十亿个更小的点替换那数百万个点，然后再次运行模拟。
接着，一些令人深感不安的事情发生了。答案变了。不仅仅是微小的变化，而是一种根本性的改变。在你新的、更“精确”的模拟中，破坏桥梁所需的能量骤然下降。根据你的模型，网格划分得越精细，材料就越容易断裂。如果你能将网格划分得无限精细，那么将桥梁折成两段似乎根本不需要任何能量。这当然是彻头彻尾的胡说八道。它违背了我们所知的一切真实世界规律，包括热力学第一定律。你不可能免费得到一个断裂！
这不是软件中的一个错误。这是底层理论中一种深刻的弊病，一个我们称之为病态网格敏感性的问题。这场危机源于一个看似无害的假设，几个世纪以来，这一直是连续介质力学的基石：局部性假设。传统模型是局部的；它们假设材料在任意点的响应——其拉伸、屈服或断裂的决定——仅取决于该精确点的条件（如应力和应变）。这是一个每个点都是特立独行、对其邻居毫不知情的世界。
在材料开始失效之前，这套理论工作得非常完美。许多材料，从混凝土到金属再到骨骼，都表现出一种称为软化的行为。在达到峰值强度后，额外的拉伸实际上会导致其能承受的应力减小。想象一下拉扯一块太妃糖，直到它开始“颈缩”；颈缩区域就在软化。在局部模型中，一旦软化开始，控制方程的数学特性就会发生根本改变。问题变得不适定。这种数学上的不稳定性在模拟中会产生毁灭性的物理后果：所有的变形都会迅速集中到最薄弱的区域。这种现象称为应变局部化，它会形成一个失效带。并且由于局部模型没有内置的尺寸感，没有任何东西能阻止这个失效带塌缩成一条无限细的线。在计算机模拟中，“无限细”变成了“单个单元的宽度”。当你缩小单元时，你就缩小了失效区，而在这个趋于消失的体积内耗散的能量也虚假地降至零。
局部性假设，我们模型的基础，将我们引向了物理上的荒谬。这告诉我们，当材料开始断裂时，连续介质假设本身也正在崩溃。一个点不再是一个自成一体的世界。为了修复我们的模型，我们必须教会它们这个基本真理。我们必须教会它们进行非局部思考。
真实世界不是局部的。晶格中的原子感受到邻近原子的拉扯和推挤。混凝土块中的骨料相互挤压，将力分布在一个有限的区域上。岩石中的微裂纹形成复杂的网络，在远大于单个点的距离上传递应力。失效是一个集体事件。
因此，治愈局部性弊病的良方，就是将这种“邻里间的对话”直接构建到我们的数学模型中。我们必须创建非局部损伤模型。其核心思想在其简单性中体现出革命性：点 处的材料状态不应仅取决于 处的条件，而应取决于 周围一个小邻域内条件的加权平均值。
这在我们的物理学中引入了一个新的基本参数：内禀长度尺度，通常用 表示。这不仅仅是一个可供调整的数值旋钮；它是一个我们必须像测量密度或刚度一样去测量的新材料属性。长度 代表了材料微观结构组分相互“交谈”的特征距离。它可能是混凝土中沙粒的平均尺寸，合金中金属晶粒的直径，或者是可见裂纹前方微裂纹激烈形成区域（即“断裂过程区”）的大小。通过引入一个物理长度尺度，我们给了模型一把尺子。我们告诉它：“你不能将失效区域局部化到比这更小的范围。”
那么，我们如何用数学来编码这种“邻里间的对话”呢？两种优雅的思想流派应运而生。乍一看它们截然不同，但正如我们将看到的，它们密切相关。
第一种方法，积分型非局部模型，或许是非局部思想最直接的表达。想象一下，在每个点上，它不是单方面决定承受损伤，而是举行一次投票。它向由内禀长度 决定的半径范围内的邻居征求意见。它征求的量是某种驱动损伤的应变度量，我们称之为“局部等效应变”，记为 。
然后，点 处的非局部等效应变 的计算，不是看它自身的应变，而是通过计算来自其所有邻居 的投票的加权平均值：
在这里， 是一个权函数。它赋予近处的点强烈的发言权，而远处点的发言权则微弱且递减。该函数的设计使其影响在距离 之外实际上消失。正是这个平均应变 ，驱动了点 处损伤的演化。
这种平均化达到了什么效果？想象一个小裂纹尖端处有一个尖锐的应变峰值。局部模型会看到这个峰值并立即在那里引发灾难性的损伤。然而，积分模型会将这个高峰值与周围较低的应变进行平均。得到的非局部应变被平滑和钝化了。这具有强大的稳定作用。用信号处理的语言来说，非局部平均化就像一个低通滤波器，平滑了应变场中那些否则会导致模拟不稳定的“突兀”高频变化。这种简单的平均化行为引入了急需的内禀长度尺度，并确保模拟的失效区具有与 相关的、真实的有限宽度，从而恢复了问题的适定性。
第二种方法看起来完全不同。梯度增强模型不是进行平均，而是对损伤状态从一点到另一点的突变程度引入一个惩罚，或称一种“税”。这个思想被编码在材料的亥姆霍兹自由能 中，你可以将其视为储存的弹性势能。在梯度模型中，这个能量不仅依赖于应变 和损伤 ，还依赖于损伤的空间梯度 ：
看最后一项！它表明，无论何处损伤场具有陡峭的梯度——即损伤在空间中变化迅速的地方——能量就更高。大自然总是追求经济，倾向于寻找能量最低的状态。这个梯度项使得尖锐的裂纹在能量上变得“昂贵”，迫使系统更倾向于从无损到完全损坏的更平滑、更渐进的过渡。这个过渡区的宽度由我们的老朋友——内禀长度尺度 控制。更大的 意味着对尖锐度的“税”更高，迫使失效区变得更宽。
当我们从这个能量原理推导控制方程时，梯度项在损伤演化方程中产生了一个拉普拉斯算子（）。拉普拉斯算子在物理学中是著名的“平滑”或“扩散”算子。它的存在正是对抗病态局部化趋势的数学机制，同样确保了失效区具有有限宽度，并使问题得到正则化。
所以我们有了两种成功的策略：平均应变的积分模型和惩罚尖锐损伤过渡的梯度模型。一个基于积分，另一个基于微分。它们似乎是截然相反的。大自然真的会有两种不同的方式来做同一件事吗？
这里蕴含着一个真正科学之美的时刻。事实证明，这两个模型不仅仅是竞争对手；它们是近亲。如果你采用积分模型的定义，并假设应变场在长度尺度 上相对平滑，你可以在积分内部对该应变场进行泰勒级数展开。经过一番微积分运算，你会发现一个惊人的结果：
非局部积分应变约等于局部应变加上一个与其二阶导数成正比的项，其中比例常数 与权函数的几何形状和内禀长度 直接相关。这个方程看起来与梯度模型产生的亥姆霍兹型方程惊人地相似。事实上，可以证明，在二阶项的精度内，梯度模型不过是积分模型的一个微分近似。两条路径在此汇合。这种隐藏的统一性是稳健物理理论的标志，它让我们相信我们正走在正确的轨道上，并允许我们为特定问题选择更方便的公式。
让我们把这一点说得更具体些。这种非局部机制实际上对我们模拟一个断裂物体做了什么？让我们再次考虑那个带缺口的试件，那个给我们的局部模型带来如此多麻烦的试件。现在我们用一个包含内禀长度 的非局部模型来建模。我们进行一系列模拟，保持其他一切不变，只改变 的值。
极限 ： 当我们将内禀长度设为非常小时，我们的非局部模型退化回那个旧的、有问题的局部模型。模拟显示出灾难性的破坏：裂纹在缺口尖端出现并几乎瞬间扩展，承载能力急剧下降。材料表现出非常脆性的行为。
中等大小的 ： 现在，我们将 设置为一个物理上现实的值，比如说，我们混凝土试件中骨料的尺寸。一个迷人的变化发生了。试件在失效前能承受的峰值载荷增加了。为什么？因为非局部平均平滑了尖锐缺口尖端的强烈应变集中。只有当一个由 定义体积的整个邻域内的点都“同意”应变很高时，损伤才会开始。这种共同承担的责任延迟了失效的发生。
更大的 ： 当我们进一步增加 时，不仅峰值载荷继续增加（到一定程度），而且峰后的行为也变得温和得多。破坏更具延性。载荷不是突然断裂，而是逐渐减小。这是因为损伤被迫在一个更宽的过程区内扩散，该区域由更大的 决定。能量的释放是渐进和可控的，而不是灾难性的。
内禀长度尺度 不仅仅是一个数学上的修正。它是一个控制预测失效特征的参数：它支配着从脆性到延性行为的转变，并设定了结构的强度。它将我们的模型从一个简单的拉伸描述转变为一个丰富的断裂预测器。
在发现了一个更深刻的物理理论之后，一个最后的实际问题仍然存在：代价是什么？通常情况下，更高真实性的代价是更高的计算成本。
积分模型，虽然可能在物理上更直观，但在计算上可能要求很高。为了更新每个点的状态，计算机必须名副其实地向其在半径 内的所有邻居征求意见。对于一个精细的网格，这可能意味着每个点都要进行数百万次交互，导致算法的成本随着我们细化网格而急剧增加。
梯度模型，作为一个微分方程，更自然地融入了大多数有限元软件的结构中。它产生的代数系统是稀疏的——意味着每个点只与其直接的单元邻居相连。借助像多重网格法这样的现代数值求解器，这些系统可以以惊人的效率求解，其求解时间通常与模型中的点数成线性关系。
因此，我们面临着一个经典的工程权衡：积分方法概念上的直接性与梯度近似计算上的效率。但由于它们之间美妙的统一性，我们理解了它们之间的联系，可以做出明智的选择。通过摒弃局部性的虚构，拥抱物质的邻里特性，我们不仅修复了我们有问题的模拟，还发现了一种更深刻、更优美、更具预测性的方式来描述世界。
在上一章中，我们开始了一段奇特的旅程。我们看到，一个简单、直观的材料图像——将其视为一组独立的点，仅根据局部发生的情况做出决策——在材料开始失效时会导致一场数学和物理上的灾难。模型预测裂纹应在零厚度区域形成，耗散零能量，这个结果与我们所知的真实世界的一切都背道而驰。我们发现，解决方案是放弃这种极端的局部性。我们必须赋予材料点“与邻居交谈”的能力。通过允许一个点的状态受一个由特征“内禀长度” 定义的小邻域内的平均状态影响，灾难得以避免。
但这种非局部思想远不止是一个巧妙的数学补丁。它是一把钥匙，开启了广阔而美丽的物理学和工程学景观，揭示了初看之下毫无关联的现象之间的深刻联系。它解释了为什么一个微小的混凝土卵石比一座巨大的大坝按比例来说更坚固，指导工程师设计更安全的结构，甚至迫使我们重新思考我们关于“材料属性”这一最基本概念的理解。现在，让我们来探索这片非凡的景观。
想象一下，你是一名工程师，正在设计一个关键部件，比如说，用于飞机机翼。你需要确保它在应力下不会失效。今天，我们严重依赖计算机模拟来测试这些设计。但模拟的好坏取决于它所基于的模型。如果你的模型是那种病态的局部模型，你的预测不仅是错误的，而且是毫无意义的，完全取决于你模拟的“像素尺寸”或网格。
这就是非局部模型成为不可或缺工具的地方。通过引入内禀长度 ，模型的预测不再受计算网格的束缚。模型现在预测失效将发生在一个有限、可预测宽度的带内，从而抹平了困扰局部模型的不可能出现的尖锐应变峰值。模拟终于获得了客观性。
但是我们如何能确定我们新获得的客观模拟是正确的呢？科学家必须是怀疑论者。我们必须设计一个测试。一个关键的程序是“网格加密研究”。我们一遍又一遍地运行相同的模拟，每次都使计算网格更精细。然后我们观察哪些物理量会稳定下来并收敛到一个稳定的值。使用一个合适的非局部模型，不仅总的力-位移曲线会稳定下来，而且在形成断裂过程中耗散的总能量也会稳定下来。这个耗散的能量必须收敛到一个特定的值——材料的断裂能，记为 ，它是一个可测量的物理量，代表材料的韧性。当我们的模拟能量账单与真实材料的能量账单相符时，我们就可以开始建立信心。
这就引出了最实际的问题：这个神奇的数字 从何而来？它只是我们凭空捏造的吗？完全不是。内禀长度 是一个真正的材料属性，就像密度或刚度一样，我们可以测量它。最优雅的方法之一是通过“尺寸效应”试验。想象一下，在一块复合材料板上钻不同大小的孔，然后拉伸它直到断裂。一个简单的比例定律可能会表明，无论孔的大小如何，失效应力都应该相同，或者遵循一个简单的规则。但实验表明事实并非如此！带有小孔的板的名义强度按比例高于带有大孔的板。这种与简单比例律的偏离——尺寸效应——正是非局部模型最显著的标志。通过测量强度如何随尺寸变化，我们可以精确地校准该特定材料的 值。一旦校准，模型就不再仅仅是描述性的；它变得具有预测性。我们现在可以用它来自信地评估任何尺寸或几何形状的结构的强度。
尺寸效应不仅仅是一个校准工具；它本身就是一个深刻的物理现象。为什么小物体相对比大物体更坚固？两种不同失效模式之间的竞争是答案的核心，而内禀长度 则是决定胜负的裁判。
对于一个非常大的结构——例如一座混凝土桥——其特征尺寸，我们称之为 ，远大于材料的内禀长度 。在这个极限下，，生长裂纹尖端的微小微裂纹区域与整个结构相比可以忽略不计。失效由推进这个宏观裂纹所需的能量所控制。这是经典线性弹性断裂力学（LEFM）的世界，它预测几何相似结构的名义强度 随尺寸减小，其比例关系为 。
现在，考虑一个非常小的物体，比如一根单纤维或一个微小的实验室试件，其尺寸 远小于内禀长度，。在这里，“断裂过程区”比物体本身还大。整个试件都参与了失效过程。它的行为更像一根被拉断的均匀链条，而不是一个带有裂纹的结构。当平均应力达到材料的固有抗拉强度 时，就会发生失效。在这个区域，强度只是一个材料属性，名义失效应力 变得与尺寸无关。
非局部模型完美地捕捉了这整个过程。它描述了从强度主导的“小尺度”世界到断裂主导的“大尺度”世界的光滑过渡。内禀长度 充当了通用标尺。通过比较结构的尺寸 和 ，我们立即就能理解它将如何失效。
科学中一个真正基本思想的力量在于其普适性。非局部相互作用的概念并不仅限于混凝土和陶瓷的脆性断裂。它无处不在。
考虑一个被拉伸的延性金属。它的失效不是由单一的尖锐裂纹驱动，而是由数百万个微观空洞的形核和生长驱动。随着这些空洞连接起来，材料软化并最终撕裂。这个过程，通常由像 Gurson-Tvergaard-Needleman (GTN) 框架这样的模型描述，也遭受着局部化的病态问题。空洞倾向于在不可能的薄带中排列。解决方法是相同的：必须使用一个非局部的孔隙度度量，在一个特征长度尺度 上对其进行平均，以使问题正则化。无论损伤是裂纹还是生长的空洞，其底层的数学弊病及其非局部疗法都是相同的。
这个思想甚至可以应用于完全不同的断裂建模方式。工程师们有时不模拟连续体内部的损伤，而是预先定义一个裂纹可能扩展的路径，并描述该线或表面上的分离物理过程。这些被称为内聚区模型。然而，如果内聚律是纯粹局部的，它也可能受到人为因素的困扰，即裂纹路径会不自然地受到计算网格的偏置。解决方案是什么？一个非局部的内聚律，其中分离力取决于一个小邻域上的平均张开度。
也许这种统一性最美丽的例证是非局部积分模型与另一族流行的正则化模型——相场模型 之间的联系。相场模型将一个尖锐的裂纹表示为一个连续的、“模糊”的损伤带，通过将损伤场的梯度引入材料的能量中来描述。乍一看，这种“梯度损伤”方法似乎与我们一直在讨论的“积分损伤”方法大相径庭。一个涉及导数，另一个涉及积分。然而，它们是同一枚硬币的两面。深入的数学分析表明，对于缓慢变化的场，这两种形式是等价的。可以推导出积分模型的特征长度与相场模型的长度之间的精确关系。例如，通过匹配两种模型在一维情况下预测的断裂能，可以建立两者长度尺度参数之间的定量联系，其具体形式取决于所选的非局部核函数和相场势函数。发现这样的联系，就像发现两个看似不同的物种共享一个共同的祖先；它指向一个单一、强大的底层原理在起作用。
我们现在来到了非局部性最深远的含义，它挑战了我们对材料的定义本身。今天，计算科学的宏大目标之一是多尺度建模：通过同时模拟其内部微观结构的细节来预测大型结构的行为。这个被称为 FE² 的想法是，拥有一个大型结构的计算机模型，并在其内部的每一个点上放置一个虚拟显微镜，观察一个微小的“代表性体积单元”（RVE）的微观结构，以计算出局部应力。
这里存在一个惊人的悖论。这个复杂且计算昂贵的过程的结果是一个宏观模型，它再次是纯粹局部的！宏观点处的应力仅取决于该宏观点的应变。因此，如果微观结构可以软化，那么多尺度模型将以与我们最初开始的简单模型完全相同的方式失效。局部化的病态再次抬头，但这一次是在一个更宏大的舞台上。你可能已经猜到，解决方法是认识到涌现出的宏观模型本身也必须是一个非局部模型，拥有其自身的、源自微观结构物理的内禀长度尺度。
这把我们引向了最后一个，也是最深刻的问题。整个连续介质力学的大厦建立在代表性体积单元（RVE）的思想之上——即如果你取一块足够小的材料，它在统计上看起来与任何其他块都相同，并能代表整体。这个属性就是物理学家所说的遍历性。它使我们能够谈论钢的“属性”或骨骼的“属性”。
但是当材料开始失效时会发生什么？一个局部化带形成了。突然之间，材料不再是统计均匀的了。它在这里有一个高度损伤的区域，而在那里有一个近乎原始的区域。失效带的位置是任意的，是自发对称性破缺的结果。RVE 的概念本身就崩溃了。遍历性被打破了。我们再也不能谈论失效材料的“属性”，因为一个样本的平均属性现在将取决于其尺寸以及失效带恰好在其内部形成的位置。
这不仅仅是一个哲学难题；这是我们建模框架的根本性崩溃。而再一次，非局部原理前来救援。通过引入内禀长度 ，非局部模型驯服了局部化的剧烈不稳定性。它确保了失效过程具有一个特征性的、有限的宽度。这使我们能够复兴代表性体积的概念，我们现在可以称之为统计体积单元（SVE）。我们可以再次定义有意义的平均属性，但前提是我们的样本体积必须远大于微观结构非均匀性的尺度和失效过程本身的内禀长度 。
因此，我们的旅程回到了原点。我们从一个“简单”的数值问题开始——一个给出荒谬答案的模拟。寻求解决方案的探索引导我们从用于校准和验证的实用工程工具，走向对尺寸效应的深刻物理理解。它揭示了连接脆性裂纹、延性空洞和不同数学形式的统一原理。并最终，它迫使我们面对并解决材料理论中的一个基础性危机，将工程模拟的世界与遍历性这一深刻的统计物理概念联系起来。让一个点与它的邻居交谈这个简单而“不合理”的想法，结果证明是所有想法中最合理、最有效的想法之一。