玻尔百科物体是如何断裂的？这个简单的问题是工程学和物理学中最深刻的挑战之一。几百年来，经典连续介质力学以其精妙的准确性描述了材料的行为，但它对微分方程的依赖使其在面对撕裂或裂纹等突变性不连续问题时显得力不从心。一种名为近场动力学（peridynamics）的革命性替代方案重新定义了这个问题，它提出某一点的力并不取决于其无限小的直接周围环境，而是取决于一个有限邻域内的相互作用。虽然这种非局部方法天然地适应了断裂问题，但其最初的表述，即基于键的近场动力学，带有一个显著的局限：它被锁定在预测一个单一、固定的材料属性，即泊松比，这一限制无法捕捉大量真实世界材料的行为。
本文旨在探索从这一初始模型演变而来的更强大、更通用的基于态的近场动力学框架。我们将回顾解决“泊松比问题”并开启材料建模新纪元的关键概念性飞跃。第一章“原理与机制”将剖析其核心理论，将基于键模型的简单的“弹簧宇宙”与基于态的公式化的集体感知进行对比，其中包括构建通往经典世界桥梁的优雅的对应模型。随后的“应用与跨学科联系”章节将展示这一稳健的理论如何应用于解决实际问题，从模拟金属的复杂失效到从零开始设计新材料，再到实现大规模、高性能的计算机模拟。
想象一下，你试图理解一大群人是如何移动的。你可以尝试为整个人群建立方程，将其视为一种连续的流体。或者，你可以观察每一个个体，根据他们紧邻的人来模拟他们的决定。近场动力学提供了第三种引人入胜的视角：它表明，要理解任何一个点的行为，你必须同时考虑它与整个邻域中其他所有点的关系。这就是非局部观点，它带来了深远的影响。
让我们从固体材料最简单、最直观的图像开始。想象它是一个巨大的三维点阵，代表原子或物质微粒。现在，用微小的、无形的弹簧将每个点与其邻居连接起来。这就是基于键的近场动力学的精髓。任何一个弹簧中的力仅取决于该弹簧被拉伸或压缩的程度。它不关心其他任何弹簧发生了什么；它的世界只有它连接的两个点。这就是我们所说的中心力模型——力沿着连接两点的直线中心作用。
这个“弹簧宇宙”是一个极其简单的想法。它很容易在计算机上模拟，并且有一个显著的内在优势：如果你将材料拉伸到足以使一个弹簧断裂的程度，你就创造了一条裂纹。你不需要任何特殊的断裂方程；断裂是从模型中自然出现的。
但这种简单性是有代价的。让我们想一想当你拉伸一根橡皮筋时会发生什么。它当然会变长，但也会变薄。它变薄的程度与变长的程度之比是一个基本的材料属性，称为泊松比，用 表示。对于像橡胶这样非常“柔软”的材料，泊松比接近 。刚性材料的泊松比可能较小。著名的软木，其泊松比接近于零——当你把它塞进酒瓶时，它不会向侧面鼓出太多。因此，一个好的材料理论应该允许我们选择任何物理上合理的泊松比。
这就是简单弹簧模型的致命缺陷所在。当你进行数学计算并对所有这些独立弹簧的行为进行平均时，你会发现材料整体表现得好像它有一个固定的泊松比。对于任何以此方式建模的三维各向同性材料，该理论预测 必须恰好是 。永远如此。无论你使用什么样的弹簧都无关紧要。基于键的模型既不能描述橡胶，也不能描述软木。这是一个美好的想法，却撞上了物理现实的壁垒。
我们如何摆脱 的束缚？正如我们所见，问题在于每个键都是孤立作用的。一个键中的力对邻近键的变形一无所知。真实的材料要复杂得多。例如，晶体中的原子相互作用的方式取决于键与键之间的角度，而不仅仅是它们的长度。材料是集体响应的。
这正是引出基于态的近场动力学的绝妙概念飞跃。我们需要给我们的物质点更多的信息。我们不再仅仅考虑单个键的伸长，而是定义了两个新的、更强大的概念：
变形态：对于任何给定的点，我们观察其延伸到一定距离（作用范围）内的所有键族。变形态，记作 ，是该族中所有新的、变形后的键向量的集合。它是整个邻域如何被扭曲的完整快照。
力态：现在是关键的一步。与单个键相关联的力向量 不再是该键伸长的简单函数。相反，它成为整个变形态 的函数。一个键中的力现在可以取决于其邻域中所有其他键发生的情况。它获得了“集体意识”。
这一改变带来了戏剧性的后果：一个键中的力不再必须与键本身平行。它可以有垂直于键的分量，用以抵抗邻域的剪切或扭转。中心力的魔咒被打破了。
一旦我们从变形态中获得了这种集体信息，我们就可以开始做一些巧妙的事情。我们可以问：“总体而言，这个邻域是在被挤压（经历体积变化）还是仅仅在被剪切（改变其形状）？”
这就是常态的基于态 (OSB) 模型背后的原理。我们可以为每个点定义一个单一的数值，称为体胀 ，它是其族中所有键伸长量的加权平均值。这个数值告诉我们该点的总体积变化。然后，我们可以观察每个单独的键，计算出其伸长中有多少贡献给了这个体积变化，又有多少是剩余的——即它的偏量伸长，它对应于形状变化。
现在，我们可以建立一个更复杂的力定律。一个键中的力可以有两个部分：
一个键 的完整力态可以表示为这样的形式： 因为我们现在有两个独立的旋钮可以调节——体积模量和剪切模量——我们不再被锁定在它们之间的固定关系中。我们可以选择它们的值来表示任何物理上可能的泊松比！“软木问题”解决了。这个模型被称为“常态”，因为一个键中的力向量仍然被约束为指向原始键的方向，但它的大小现在取决于邻域的集体状态。
OSB 模型是向前迈出的一大步，但还有一种更优雅、更强大的方法。这就是非常态的基于态模型，也被称为对应模型。它建立在一个非常实用的想法之上：我们有超过两个世纪的经典连续介质力学经验。我们有非常成功的模型，可以描述从钢铁到橡胶再到土壤的一切。如果我们能建立一个可以利用所有这些现有知识的近场动力学理论呢？
对应模型正是通过在近场动力学的非局部世界和经典力学的局部世界之间建立一座桥梁来实现这一点的。其过程如下：
寻找局部图像：在点 处，我们观察其变形的邻域（即变形态 ）。我们问：“能够近似这整个扭曲键集合的、由一个变形梯度张量 描述的最佳‘经典’变形是什么？”这有点像用一个平面对一个凹凸不平的表面进行拟合。答案是通过加权最小二乘法最小化得到的，这为我们提供了这个非局部变形梯度的绝佳公式： 这里， 是一个影响函数，它对不同键的重要性进行加权，而 是至关重要的形状张量。形状张量是材料参考构型的一个几何属性，它说明了邻域中物质点的分布情况。
使用经典定律：一旦我们有了张量 ，我们就进入了经典世界！我们可以使用任何经过时间考验的本构定律来计算经典应力张量，例如，第一皮奥拉-基尔霍夫应力 。你想模拟线性弹性材料吗？使用胡克定律：。想模拟橡胶吗？使用超弹性定律。金属塑性？使用塑性模型。大门是敞开的。
返回非局部世界：现在，我们需要将这个经典应力 转换回近场动力学的力态 。对应原理为我们提供了最后一块拼图： 这非同寻常。力态是使用经典应力定义的。现在的力明确是非中心的；向量 通常不沿着键 的方向。这个非中心部分正是表达全部材料行为所需要的。而且这个特定的公式不仅仅是一个猜测；它经过精心构造，以确保像线动量和角动量守恒这样的基本物理定律得到完美满足。
对应模型让我们两全其美：一个能够自然处理像裂纹这样的不连续性的非局部框架，由经典力学中稳健且被充分理解的本构定律驱动。
至此，你可能会想：“如果对应模型最终只是使用经典定律，那么影响函数 的选择还重要吗？”毕竟，如果我们用两个不同的影响函数 和 来校准两个不同的模型，以匹配相同的杨氏模量和泊松比，它们难道不应该给出相同的答案吗？
对于简单的均匀拉伸，它们会给出相同答案。但如果我们施加一个非均匀的变形，比如一个具有特定波长的缓和波浪呢？在这里，近场动力学更深层的魔力就显现出来了。事实上，这两个模型将预测出略有不同的应力响应。
原因在于非局部模型不仅对位移的一阶导数（它给我们应变）敏感，而且对高阶导数（它关系到变形的曲率）也敏感。结果表明，对于波浪状变形，近场动力学的应力略低于经典理论的预测。这种现象被称为频散软化。这种软化的程度取决于影响函数的高阶矩——即其形状的更精细细节。
这告诉我们一些深刻的东西。影响函数的选择不仅仅是方便与否的问题；它是一个关于材料如何响应快速变化的变形的物理陈述。它控制着不同长度的波如何在材料中传播，这是经典理论完全忽略的真实物理效应。
而且，其通用性还不止于此。通过设计一个不仅依赖于键的长度，还依赖于其方向的影响函数，我们可以将各向异性直接构建到模型的结构中。这使我们能够捕捉像木材这样的材料的行为，它沿着纹理方向比横跨纹理方向要硬得多，或者现代的纤维增强复合材料。
从一个简单但有缺陷的弹簧模型出发，我们已经抵达了一个复杂而强大的理论。通过拥抱集体态的思想并建立与经典世界的对应关系，基于态的近场动力学提供了一个不仅在数学上优雅，而且与构成我们世界的材料丰富而复杂的现实紧密相连的框架。
在我们迄今的旅程中，我们已经探索了基于态的近场动力学的优雅架构，这一理论并非用点的微分语言，而是用相互作用的积分语言重写了固体力学定律。我们已经看到这种非局部视角如何自然地容纳了断裂的突变性不连续。但是，一个物理理论，无论多么优美，只有当它与世界联系起来时——当它能解释我们所见，预测我们所未见，并赋予我们能力去建造我们所能想象之物时——才真正焕发生机。现在，让我们踏上旅程的下一段，探索基于态的近场动力学如何作为一座强大的桥梁：一座通往成熟的经典力学世界的桥梁，一座通往其无法触及的现象的桥梁，一座从微观到宏观的桥梁，以及最后，一座从抽象方程到具体、大规模计算机模拟的桥梁。
物理学中的一个新理论并不会简单地抛弃旧理论；它必须恭敬地包含它。在我们能够自信地使用近场动力学去探索断裂这个未知前沿之前，我们必须首先验证它在经典力学占主导地位的那些平静、已充分勘测的领域中表现正常。想象一块无限大、均匀的橡胶块被轻柔地剪切。每个点都移动了一点，但变形是平滑且处处相同的。在经典力学中，应力在整个块体中是恒定的，因此其散度——任何一点的净内力——为零。块体处于平衡状态。近场动力学对此有何看法？如果我们煞费苦心地对一个点的邻域内（其作用范围内）所有微小的拉力和推力进行积分，我们会发现一个极其简单的结果：总力恰好为零。来自所有方向的相互竞争的力完美地抵消了，正如对称性所要求的那样。这不是一个微不足道的结果；它是一致性的关键标志，向我们保证，非局部公式并非某种数学幻想，而是牢固地植根于力学的基础原理之上。
这种一致性检验令人安心，但在连接这两个世界方面，真正的神来之笔是为非常态的基于态模型发展的对应原理。想一想过去两个世纪里发展起来的庞大材料模型库——用于钢铁、混凝土、木材以及飞机上使用的复杂各向异性复合材料的模型。我们必须抛弃所有这些知识从头开始吗？对应原理说：不。它扮演了一个出色的“万能翻译器”的角色。
这种翻译背后的机制是一个纯粹的几何对象，称为形状张量，通常用 表示。该张量表征了一个点邻域的几何形状。通过使用它，对应框架可以获取一个点周围的非局部变形模式，并计算出一个等效的局部变形梯度。这就是关键：它将近场动力学的非局部语言翻译成传统本构模型所理解的经典应变的局部语言。然后，我们可以将此应变输入到任何经典模型中——例如，一个用于横观各向同性材料（如纤维增强聚合物）的模型——以获得经典应力。框架接着使用完全相同的形状张量将此应力翻译回近场动力学的键力语言。这是一个惊人强大的想法。它允许工程师将他们信任的、经过良好校准的经典材料定律直接“插入”到近场动力学引擎中。这意味着我们可以使用近场动力学模拟一个复杂的、真实的各向异性材料的行为，既能利用其断裂模拟能力，又无需放弃数十年来建立的材料科学。
建立这些通往经典世界的桥梁的真正目的，是拥有一个坚实的基础，以便 venturing into the unknown——即经典理论失效的材料失效领域。虽然基于键的近场动力学通过简单地断开键来处理脆性断裂，但基于态的模型，特别是通过对应框架，可以描述更为复杂的失效过程。
考虑金属的行为。在断裂之前，它会发生永久变形——它会屈服并流动。这就是塑性现象。经典力学对塑性有丰富而成功的理论，涉及屈服面和流动法则等概念。使用我们刚才讨论的相同对应原理，我们可以将这整个经典塑性框架嵌入到近场动力学中。我们计算等效的局部应变，将其输入经典塑性模型以确定应力和任何塑性流动，然后将产生的力转换回近场动力学模拟中。这使我们能够模拟延性构件的整个生命周期：它弹性变形，开始屈服，塑性流动，最终，累积的损伤可能导致裂纹的形核和生长，而近场动力学对此能自然处理。这种统一的方法对于碰撞安全分析、金属成形等制造过程以及预测极端载荷下的失效至关重要。
到目前为止，我们已经讨论了近场动力学如何再现已知材料的行为。但也许其最激动人心的应用在于其预测新材料行为的能力。在这里，近场动力学成为一种创造的工具，将材料结构的微观世界与工程性能的宏观世界联系起来。
想象你是一位材料设计师。你可以构想一种具有特定微观结构的材料。例如，如果你能创造一种材料，其中粒子间的微小键在水平方向上的强度是垂直方向的两倍，会怎么样？这是一种各向异性的“微观模量”。由这种材料构成的块体的宏观属性会是什么？近场动力学提供了答案。通过在键的层面上定义这些简单的局部规则，近场动力学积分会自然地将它们均匀化。我们可以在我们的虚拟材料上进行模拟的“单轴拉伸试验”，在不同方向上拉伸它，模拟将揭示出每个方向上涌现出的宏观杨氏模量。例如，我们可以发现，当沿着其强微观键的方向拉伸时，材料要硬得多。这种能力将计算机变成了一个用于材料发现的实验室。它为超材料、复合材料和仿生结构的理性设计打开了大门，使我们能够在实验室合成之前就设计出具有定制属性的材料——前所未有的刚度重量比、新颖的振动阻尼或特定的断裂行为。
一个物理理论，无论其威力如何，除非我们能用它来进行定量预测，否则它仍然只是一个奇观。对于像近场动力学这样计算要求极高的理论，从黑板走向高性能计算机是一个巨大的挑战，本身就是一个引人入胜的故事。
“非局域性的诅咒”在于每个粒子都与大量的邻居相互作用。一种天真的“暴力”方法，即检查系统中每一对粒子之间的距离以确定它们是否是邻居，会导致计算成本按粒子数的平方 增长。对于一个拥有一百万个粒子 () 的模拟，这在每个时间步中将需要一万亿次计算，这是一项计算上无法承受的任务。这种规模的增长会将近场动力学局限于玩具问题。
计算科学的胜利在于开发出驯服这种复杂性的算法。通过巧妙地将粒子组织成一个“单元”网格（一种称为单元链接列表的方法），我们可以将对给定粒子邻居的搜索限制在其自身的单元和紧邻的单元中。这个简单而深刻的想法改变了问题。计算成本不再由系统中的粒子总数决定，而仅由局部密度决定。规模变为线性，即 。正是这一算法上的飞跃，使得模拟数百万甚至数十亿个粒子成为可能，从而能够分析现实的工程系统。
但是，对于运行如此复杂代码的大规模模拟，一个关键问题出现了：“计算机的答案正确吗？”我们如何信任输出？在这里，科学方法通过验证与确认进入了计算领域。一种强大的技术是精确解方法。我们为位移场发明一个平滑的、已知的数学解，并将其代入近场动力学方程中，从而解析地计算出相应的力场。然后我们使用这个位移场运行我们的计算机模拟，并将数值计算出的力与精确的解析答案进行比较。当我们细化我们的模拟网格，使其越来越精细时，计算结果与真实解之间的误差应该以一种可预测的方式减小。对于一个表现良好的方法，这个误差通常随着网格间距的平方减小，这被称为二阶精度。在数值实验中看到这种预期的收敛速度，让我们相信我们的代码是底层数学理论的忠实实现，是科学发现的可靠工具。
通过这一系列桥梁，基于态的近场动力学展现了其真实特质：不是一个孤立的断裂理论，而是一个统一的框架。它无缝地连接到经典世界，同时又超越了其局限。它将微观结构的设计与宏观系统的功能联系起来。通过物理学和计算机科学的结合，它已成为现代科学家和工程师的强大而实用的工具，准备好应对当今的挑战，并帮助设计未来的材料。