玻尔百科理解材料在载荷下如何发生永久变形是现代工程学的基石。塑性理论为此提供了一个框架，但一个核心问题依然存在：一旦材料屈服，它会沿什么方向流动？一个优雅且在历史上具有重要意义的答案是关联流动法则，该法则假定塑性应变的方向由定义材料强度的同一函数唯一确定。这创造了一个优美简洁且统一的模型，对许多金属材料都非常有效。
然而，当面对其他关键材料（如土壤、岩石和混凝土）的行为时，这幅优雅的图景便会破碎。对于这些材料，关联流动法则对体积变化的预测与实验现实形成鲜明对比。这种差异揭示了一个根本性的知识空白，迫使我们在理论的简洁性与预测的准确性之间做出选择。本文通过引入非关联流动的概念来应对这一挑战，这是对塑性理论的一种实用而强大的修正。
在接下来的章节中，我们将探讨这一关键概念。首先，在“原理与机制”中，我们将解构理想的关联流动法则，识别其在剪胀性方面的局限性，并引入非关联模型作为解决方案，同时审视其中涉及的理论权衡。随后，在“应用与跨学科联系”中，我们将看到该理论如何成为从岩土力学到冶金学等领域不可或缺的工具，并揭示其对计算和结构分析的深远影响。
在我们探索材料如何屈服和流动的旅程中，我们现在到达了问题的核心。我们对屈服面有了一个大致的概念——这是一个位于应力抽象空间中的边界，它将弹性行为与永久的塑性变形区分开来。但一个关键问题依然存在：一旦材料决定屈服，它会沿什么方向流动？如果你挤压一块金属，它会变短变粗。但具体程度如何？是否存在一个普适的法则？这个问题的答案将我们引向一条引人入胜的道路，它始于一个异常优雅的想法，最终以一个反映真实世界优美复杂性的务实折衷而告终。
想象一下，你正站在一个连绵起伏的山丘上，把一个球放在山坡上。它会朝哪个方向滚动？它会沿着最陡峭的路径直滚而下。任何数学家都会告诉你，这条路径与该点的等高线垂直，即正交。在塑性理论的早期，一个类似的、优美的想法被提出：也许塑性流动的行为也是如此。“山丘”是应力的景观，“等高线”则是屈服面。这个被称为关联流动法则的想法指出，塑性应变率的方向在当前应力状态下总是垂直于屈服面。
这是一个非常简洁而强大的概念。这意味着定义屈服边界的单一数学函数，我们称之为屈服函数，同时也控制着流动的方向。这个方向就是屈服函数的梯度 。这就是为什么我们称之为“关联”——流动法则与屈服函数直接关联。
这个“正交法则”不仅仅是一个漂亮的猜测。它与稳定性的基本原理密切相关。Daniel C. Drucker 的一条公设表明，稳定材料应以某种可预测的方式行事。该公设的推论之一是，对于具有凸屈服面（意味着它没有任何向内弯曲的部分）的材料，正交法则必须成立 (2711718)。稳定性与关联流动之间的这种联系是深刻的。它意味着塑性耗散——在塑性流动过程中转化为热的能量——被最大化了。换句话说，对于给定的塑性应变增量，材料会“选择”屈服面上能做最多功的应力状态 (2654976)。在这幅理想的图景中，自然似乎遵循着最大努力原则。
这种优雅是引人注目的：一个函数 同时定义了强度极限和流动法则。对于许多材料，特别是当压力不显著影响其屈服的金属（如在 von Mises 屈服准则中），这幅美丽的图景表现得非常好 (2711713)。塑性流动是不可压缩的，意味着体积不发生变化，且理论与实验吻合。
然而，自然往往比我们最优雅的理论更为微妙。当我们将注意力从延性金属转向其他关键材料，如土壤、岩石、混凝土和颗粒状集合体时，这个美丽的体系出现了严重的裂痕。这些材料是压力相关的；它们的强度随着被挤压的程度而增加。想象一堆沙子：如果它被约束和压缩，你就更难将手指插进去。
像莫尔-库仑或德鲁克-普拉格这样的屈服准则被提出来描述这种摩擦行为。它们定义的屈服面在应力空间中不再是简单的圆柱体（如 von Mises 准则），而更像是圆锥体或棱锥体，其强度随静水压力增加而增加。那么，如果我们将关联流动法则应用于这些屈服面，会发生什么呢？
让我们在一个简化的二维应力空间中将其可视化，坐标为压力 和剪应力 。对于摩擦性材料，莫尔-库仑屈服面是一条斜线 (2559749)。这条线的斜率与材料的内摩擦角 有关。如果我们应用正交法则，塑性流动矢量必须与这条线垂直。因为这条线是倾斜的，所以法向矢量在剪切和压力方向上都有分量。压力方向上的流动对应于体积的变化。关联流动法则预测，当材料被剪切时，它的体积也必须膨胀。这种现象被称为剪胀性，而且它确实存在。如果你剪切一袋密实的沙子，它会膨胀。
问题就在这里：关联流动法则预测的剪胀量与摩擦角 直接耦合。摩擦力越大，材料被预测的膨胀就越多。而对于大多数土壤和岩石来说，这个预测是极其不准确的——其高估的幅度通常非常大 (2711713)。在面对这些常见且至关重要的材料的实验数据时，关联流动法则的优雅统一性就崩溃了。
当一个美丽的理论与事实不符时，我们面临一个选择。我们可以抛弃它，或者看看一个巧妙的修正能否挽救它。在这种情况下，物理学家和工程师选择了后者。问题在于，单一的函数 被要求同时完成两项工作：定义强度和定义流动。而它在第二项工作上失败了。
解决方案是引入第二个独立的函数，称为塑性势 (2671041)。新的法则是：
这被称为非关联流动法则，因为流动不再与屈服函数相关联。通过选择 ，我们可以将材料的强度与其流动行为解耦。
让我们回到 平面中莫尔-库仑材料的图像 (2559749)。
塑性流动矢量现在垂直于 这条线。如果我们选择一个小于摩擦角 的剪胀角 ，新的流动矢量将更偏向垂直方向。它在压力方向上的分量会更小，这意味着更少的体积膨胀。我们可以选择 来匹配实验观察到的剪胀性，同时保持摩擦角 来正确模拟材料的强度。我们用精度换取了优雅。非关联模型预测的体积应变与关联模型的体积应变之比，就是它们各自斜率参数的比值 (2893867)。这让工程师能够精确控制模型的预测。
这个巧妙的修正并非没有代价。通过引入第二个函数并打破屈服与流动之间的联系，我们失去了一些伴随关联模型而来的优美数学性质。
首先，我们明确违反了最大塑性耗散原理 (2655028)。材料不再遵循“最大努力”的路径。实际的应力状态不再是对于给定塑性应变做最多功的状态。这对宏观理论如极限分析产生了影响，其中坍塌载荷的上下限之间的清晰对应关系变得模糊 (2654976)。
其次，Drucker 公设所保证的材料稳定性也丧失了。在某些条件下，非关联模型可以预测不稳定的行为，例如应变局部化，即变形集中在非常窄的带中。
第三，对于计算工程师来说，一个非常实际的后果是，控制方程失去了一个关键性质：对称性。当使用有限元方法等方法模拟这些材料时，材料模型的内部刚度矩阵（它关联了应力率和应变率）变得非对称 (2631371)。这意味着用于求解方程的数值算法变得更复杂、效率更低，并需要更多内存。这是一个为获得更真实的材料描述而付出的巨大实际代价。
在付出了所有这些代价之后，一个深刻的问题出现了：非关联模型在物理上是否合法？它是否违反了某些基本的自然法则？这里的最终仲裁者是热力学第二定律，该定律要求在任何不可逆过程中，总耗散（损失为热的能量）不能为负。
对于塑性材料，这归结为塑性耗散率 必须为非负的条件。让我们为我们的非关联法则 检验这一点，其中 是一个表示塑性流动速率的正标量。耗散为 。为确保其为非负，我们需要 。
事实证明，这个条件可以在不要求 的情况下得到满足。一个充分条件是，塑性势 是一个凸函数，并且其在零应力状态下取最小值 (2711752)。这给了我们答案：非关联流动本质上不违反热力学第二定律。虽然它牺牲了稳定性公设的优雅和计算上的对称性，但它可以被构造成与基本能量定律完全一致的形式。
最终，非关联流动的故事是科学过程的一个完美范例。我们从一个简单、优美的理论开始。我们用现实检验它，并发现其局限性。然后我们修改它，使其更复杂，但也更强大、更准确，并在此过程中牺牲了一些优雅。与此同时，我们不断回顾物理学的基本定律，以确保我们的新创造，尽管务实而复杂，仍然建立在坚实的基础之上。
既然我们已经掌握了非关联流动的原理，你可能会想：“这仅仅是一个数学上的奇趣吗？”远非如此。这才是故事真正激动人心的地方。屈服函数 和塑性势 之间的区别不仅仅是学术上的细微差别；它是一个至关重要的工具，能让我们更深刻、更准确地理解我们周围的世界。它是将自然的卡通化简笔画与能够预测真实材料行为（从我们城市下方的土地到我们最先进机器的金属骨架）的模型区分开来的关键。让我们踏上一段旅程，看看这些想法将我们带向何方。
非关联流动最直接和最根本的应用是它能够描述材料的真实面貌。想象一堆密实的沙子。要使其“屈服”或开始流动，你需要用一定的剪应力去推它，而这种抵抗剪切的能力与颗粒间的摩擦有关。这正是屈服函数 所描述的。但是，一旦它开始流动，会发生什么呢？那些紧密堆积的颗粒必须相互推挤和滚动。为了做到这一点，它们必须相互推开，整个沙堆的体积就会膨胀。这种现象被称为剪胀性。
在一个简单的“关联”世界里，人们可能会猜测材料的膨胀量（剪胀性）与其内摩擦直接锁定。但实验告诉我们，对于包括土壤、岩石和混凝土在内的许多材料，情况并非如此。决定材料何时屈服的摩擦角，通常显著大于描述其变形时膨胀程度的剪胀角。关联流动法则，其流动方向由屈服函数本身决定，强制将这两个属性联系在一起。这常常导致对预测的体积膨胀量的大幅高估。
这时，非关联流动法则就来救场了。通过引入一个独立的塑性势 ，我们可以解耦这两种效应。屈服函数 继续定义“应力预算”——材料在屈服前能承受多少应力。但塑性势 拥有自己独立的参数来描述塑性流动的方向——特别是体积膨胀的量,。通过调整 中的参数（通常称为剪胀参数 或 ），我们可以精确匹配观察到的体积变化，而不会破坏我们对材料强度的描述。
但我们如何确定这不只是在玩一个数学游戏呢？我们可以测量它。想象一下，取一根对压力敏感的金属棒，在标准的单轴拉伸试验中拉伸它。当我们拉伸它时，它在轴向变长，在横向变细。对于塑性变形，我们可能天真地期望体积保持不变。然而，通过仔细测量轴向和横向应变，并将塑性部分从弹性部分中分离出来（这是材料测试中的标准程序），我们可以计算出塑性泊松比。非关联流动法则预测，这个比率将直接取决于塑性势中的参数 。如果我们的测量揭示了塑性体积变化，我们就在自己的实验室里当场捕获了非关联性！。
非关联塑性理论的经典应用领域是岩土力学。土壤、岩石和混凝土在载荷下的行为本质上是压力敏感和非关联的。当工程师设计地基、隧道或大坝时，他们的计算机模型必须考虑到土壤的摩擦角 不等于其剪胀角 的事实。像用于黏土的修正剑桥模型等先进模型明确地包含了非关联流动法则，以更好地预测土壤在不同加载条件下将如何压实或剪胀，这对于预测沉降或破坏至关重要。
但这一概念的影响力远不止于土木工程。考虑冶金学与损伤力学领域。当一个金属构件被拉伸到极限时，微小的孔洞会在其内部形成和长大。这些孔洞的生长导致材料膨胀，这种效应本质上是一种剪胀。相反，在高压下，这些孔洞可以被压碎，导致材料压实。为了模拟延性金属的失效，科学家和工程师使用如 Gurson-Tvergaard-Needleman (GTN) 模型等复杂模型。而在这些模型的核心是什么？你猜对了：一个非关联流动法则。屈服面描述了金属能承受的应力和孔隙含量的组合，而一个独立的塑性势则控制着孔洞生长或收缩的速率，最终导致断裂。这是一个统一原理的美妙例子：描述山体滑坡的同一个基本思想，也可以描述钢板的撕裂。
选择使用非关联流动法则并非一个简单的调整。它会在我们对物理学、数学和计算的理解中掀起深远的涟漪。
首先，考虑运动几何学。在理想塑性材料的经典塑性理论中，存在着“滑移线”，这是材料倾向于变形的自然路径。对于关联材料，会出现一个奇妙的简化：最大剪应力方向与塑性速度方向完全相同。应力特征线和速度特征线重合了。非关联性打破了这种优雅的对称性。材料仍然沿着最大剪切线（由屈服函数 控制）屈服，但它移动的方向是沿着一组不同的特征线（由塑性势 控制）。应力图谱和运动图谱被解耦，揭示出一种更复杂、更微妙的破坏编排。
其次，这一选择对计算力学构成了巨大的挑战。现代工程依赖有限元法（FEM）来模拟和设计几乎所有东西。当我们在 FEM 代码中实现一个材料模型时，过程的一个关键部分是验证代码是否确实解决了我们认为它在解决的方程。对于非关联模型，一个强有力的验证测试包括在一个简单、受控的载荷下对单个点进行模拟，并检查塑性体积变化与塑性剪切变形之比是否精确匹配塑性势 中的剪胀参数，而不是屈服函数 中的摩擦参数。
更深刻的是，采用更真实的非关联模型可能会导致我们的模拟变得不稳定。大多数 FEM 代码的数值引擎——牛顿-拉夫逊方法，在底层方程组具有对称正定结构时工作得最好。非关联流动法则破坏了这一点。由此产生的“切线矩阵”（它引导数值解）变得非对称。这不是一个程序错误；这是一个特性。它是一个深刻物理性质的数学标记。这种非对称性可能导致材料稳定性的丧失，即使材料仍在硬化。这正是导致剪切带等现象的不稳定性——这些是灾难性破坏前兆的强烈变形的窄带。数值算法在收敛上的挣扎，是物理学在告诉我们，一些戏剧性的事情即将发生,。工程师们已经开发出巧妙的数值策略，如线搜索或自适应正则化，来驾驭这些险恶的计算水域，并成功模拟这些复杂材料的行为。
最后，非关联性迫使我们重新思考结构工程中最优雅和最强大的工具之一：极限分析。经典的极限分析定理允许工程师用异常简单的方法计算结构坍塌载荷的严格上限和下限。特别是上限定理，依赖于一个称为最大塑性耗散原理的美妙思想，该原理指出，在所有可能的应力状态中，真实的应力状态是那个对变形的塑性材料做最多功的状态。然而，这个原理只保证对关联材料成立。
对于非关联材料，实际耗散的能量小于这个理论最大值。因此，使用这个最大值来估计坍塌载荷的经典上限定理，可能会给出一个过高的答案。“上限”不再是保证安全的估计，并可能危险地高估结构的强度。这是一个发人深省的深刻教训。即使是我们最美丽、最可信的定理也有其局限性，而理解非关联流动的微妙物理学，能让我们看到我们知识地图上的“警告：此处有龙”的警示牌。它提醒我们，我们理解自然的旅程是一个不断完善的过程，其中每一层新的复杂性都揭示出新的挑战和更深刻、更令人满意的真理。