In search of constitutive conditions in isotropic hyperelasticity: polyconvexity versus true-stress-true-strain monotonicity

本文表明，在各向同性超弹性理论中，仅靠多凸性或真应力-真应变单调性本身都无法保证物理上合理的行为，这暗示着尽管两者的结合是解决特鲁德尔尔主问题（Truesdell's Hauptproblem）的一个极具前景的方案，但目前尚未发现能够同时满足这两个条件的全局应变能函数。

要点

想象一下你是一位大师级建筑师，正在设计一座由特殊拉伸材料构成的建筑。你的目标是编写一套规则（即“本构方程”），能够精确预测这种材料在拉伸、挤压或扭转时的行为。你希望确保你的规则永远不会预测出一些不可能发生的情况，比如当你拉得更用力时，材料突然缩回，或者仅仅因为旋转了建筑，其行为就发生了剧烈的变化。

在物理学世界中，这就是“主问题”（Hauptproblem）：我们如何编写这些规则，使它们既在数学上严谨，又在物理上符合现实？

本文探讨了两套著名的规则，科学家们曾提出用这些规则来解决这个问题。作者 Wollner、Holzapfel 和 Neff 扮演着侦探的角色，通过相互测试这些规则来询问：“如果一种材料遵循规则 A，它是否会自动遵循规则 B？”

以下是他们使用简单类比进行的调查分解。

两位竞争者

1. 多凸性（Polyconvexity）（“数学安全网”）
将多凸性想象成一个严格的数学安全网。这个规则确保建筑不会坍塌进一个数学黑洞（即解不存在的情况）。它在计算机模拟中非常受欢迎，因为检查起来很容易。

• 承诺： 如果你使用这条规则，数学运算就能顺利进行，材料在方程中不会出现奇怪且不可能的行为。
• 代价： 作者发现，仅仅通过了这个“安全网”测试，并不意味着材料在每种情况下都能表现得像一种真实、合理的材料。

2. TSTS-M++（“常识单调性”）
将 TSTS-M++（真应力-真应变单调性）想象成一条“常识”规则。它规定：“如果你拉得更用力，所需的拉力应该持续增加。如果你扭得更厉害，阻力也应该持续增加。”这就像拉伸橡皮筋；随着拉伸距离的增加，它应该变得越来越难拉，而不是突然变得容易。

• 承诺： 这条规则保证了材料在特定的测试（如直线拉伸或扭转）中表现得具有可预测性。
• 代价： 这条规则也不是万能灵药。一种材料可以遵循这条规则，但在其他方式下仍然表现得异常。

调查过程：测试规则

作者设置了两个具体的挑战，以观察一种规则是否可以取代另一种。

挑战 1：拉伸测试（单轴延伸）

• 场景： 想象像拉扯太妃糖一样，将一块材料笔直地向外拉。
• 问题： 如果一种材料遵循“数学安全网”（多凸性），它是否总是会在拉伸时变得越来越难拉？
• 结果： 不是。 作者构建了一个特定的数学模型（一种“伪造材料”），它完美地通过了多凸性测试。然而，当他们模拟拉伸它时，所需的拉力先是上升，然后突然下降，接着又再次上升。
• 类比： 这就像一辆在数学上被保证是安全的汽车，但当你踩下油门时，它会自行加速，然后突然减速，然后再加速。这不符合真实的汽车（或真实的材料）的行为。
• 结论： 单靠多凸性不足以保证拉伸时的“常识”行为。

挑战 2：扭转测试（简单剪切）

• 场景： 想象在按住一副扑克牌底部不动的情况下，将顶部的牌向侧面滑动。这就是“剪切”。
• 问题： 如果一种材料遵循“常识”规则（TSTS-M++），它是否会在扭转程度增加时，变得越来越难扭转？
• 结果： 不是。 作者构建了另一种“伪造材料”，它完美地遵循了常识规则。但是，当他们模拟扭转它时，阻力先是上升，然后下降，接着又上升。
• 类比： 想象一个门铰链，它在推开时变得越来越难，然后突然变得松动且容易推，接着又变得难以推动。这违反了“数学安全网”（具体来说是一个称为 Legendre-Hadamard 椭圆性的条件，该条件确保了稳定性）。
• 结论： 单靠常识（TSTS-M++）不足以保证扭转时所需的数学稳定性。

大局观：缺失的环节

作者得出结论：这两条规则单独来看都不够强大。

• 你需要多凸性来确保数学上的稳定性（防止在扭转时出现剧烈波动）。
• 你需要 TSTS-M++ 来确保材料在拉伸时表现得合理（拉力随拉伸而持续增加）。

终极目标： 该领域的“圣杯”是找到一套能在所有可能的变形下同时满足这两个条件的单一规则。

• 现状： 作者非常努力地尝试寻找这种“完美材料”，但无法找到一个在全局范围内（对于所有的拉伸和扭转）都适用的模型。
• 部分成功： 他们确实找到了一些“链限制型”（chain-limited）解。可以将这些理解为表现完美的材料，但仅限于达到特定极限之前（比如一条橡皮筋在达到特定长度之前表现完美，一旦达到该长度，规则就会失效）。

面向大众的总结

这篇论文是对设计材料的科学家的现实提醒。它说：“不要仅仅依赖一种数学技巧来确保你的材料模型是优秀的。”

• 如果你只检查数学安全性（多凸性），你的材料在拉伸时可能会表现异常。
• 如果你只检查常识（TSTS-M++），你的材料在扭转时可能会表现不稳定。

要真正解决理想弹性材料建模的问题，我们可能需要结合这两种规则。然而，寻找一个能在每种情况下都完美满足这两个条件的单一公式仍然是一个未解之谜，尽管作者已经提供了新的工具和部分答案，以帮助未来的研究人员破解这一密码。

技术摘要

技术摘要：探寻各向同性超弹性中的构成条件：多凸性与真应力-真应变单调性的权衡

问题陈述
本文探讨了有限弹性理论中的“主要问题”（Hauptproblem）：即在理想各向同性超弹性中，确定合适的构成约束以确保能量函数 $W$ 具有符合物理常理的行为。虽然多凸性（polyconvexity）的数学特性已得到公认（确保了极小值的存在性和拟凸性），但其机械后果尚不明确。相反，真应力-真应变单调性（TSTS-M++）这一近期被重新发现的条件，虽然能确保特定变形模式下的应力-应变轨迹单调，但缺乏多凸性那样广泛的数学稳定性保证。核心问题在于：仅凭其中一种条件是否足以保证物理行为的合理性，或者是否需要两者的结合，以及这种结合是否可以在全局范围内实现。

研究方法
作者采用了基于各向同性超弹性理论的严谨分析方法。关键方法论要素包括：

• 不变性表示：能量函数使用一组特定的不变量 $K_i$（左柯西-格林张量 $\mathbf{B}$ 的主不变量的平方根）进行参数化，这与使用主伸长率相比，简化了构成不等式的表示。
• 推导充分条件：本文根据能量函数对不变量 $K_1, K_2, K_3$ 的导数，推导出了关于多凸性和 TSTS-M++ 的新的、显式的充分条件。
• 显式反例：为了测试这些条件的充分性，作者构建了特定的应变能函数族。这些函数旨在满足其中一个条件（例如多凸性），同时显式地违反另一个条件所对应的预期物理行为（例如在单轴拉伸中的单调性）。
• 解析证明：本研究依赖于符号微分、谱分解以及对常微分方程（特别是引理 5.15）的分析，从而在不诉诸大规模数值计算（除可视化外）的情况下，证明所构建能量函数的性质。

主要贡献与结果
本文对 Neff 等人 (2024) 提出的关于多凸性与 TSTS-M++ 相互作用的四个“挑战问题”给出了明确回答：

1. 多凸性在单轴拉伸中的不充分性（挑战 ii & iii）：

   • 可压缩情况：作者构建了一个多凸性的应变能函数（式 5.37），该函数在无约束单轴拉伸中产生了非单调的柯西应力响应。具体而言，尽管该函数是多凸的且满足 Legendre-Hadamard 椭圆性条件，但其柯西应力 $\sigma_{11}$ 在伸长率 $\lambda_1$ 增加时会达到最大值然后下降。
   • 不可压缩情况：虽然未发现反例，但作者证明了任何满足 Ball (1976) 所提出的多凸性充分条件的不可压缩应变能函数，都自动导致单调的真应力响应。这显著缩小了寻找潜在反例的搜索空间，但并未证明其不可能实现。

2. TSTS-M++ 在简单剪切中的不充分性（挑战 iv）：

   • 作者构建了一个满足全局 TSTS-M++ 的应变能函数（式 5.69），但该函数在简单剪切中表现出非单调的柯西剪切应力响应。
   • 这一结果意味着该函数不是秩一凸的（因此也不是多凸的），因为秩一凸性被证明是简单剪切中单调剪切应力的必要条件（命题 5.13）。

3. 链受限解（挑战 i）：

   • 作者尝试寻找一个在整个定义域 $GL^+(3)$ 上同时满足多凸性和 TSTS-M++ 的单一函数。他们未能识别出这样的函数。
   • 然而，他们成功构建了三个“链受限”（chain-limited）应变能函数族（命题 5.2, 5.4, 5.5），这些函数在受限的变形区域内（例如受限于线、面积或体积元平均值）同时满足这两个条件。

4. 新的充分条件：

   • 本文建立了关于不变量 $K_i$ 的多凸性和 TSTSM++ 的显式充分条件（定理 3.1 和 定理 4.4）。这些条件强调了两种约束在结构上的差异，特别是在能量函数及其导数的凸性要求方面。

意义与主张
本文得出结论：仅靠多凸性或 TSTS-M++ 均不足以保证理想各向同性弹性中具有物理合理的材料行为。

• 多凸性确保了稳定性（Legendre-Hadamard 椭圆性）和单调剪切应力，但在单轴拉伸中无法保证单调的应力响应。
• TSTS-M++ 确保了单轴拉伸中的单调应力，但在简单剪切中无法保证单调剪切应力（从而无法保证秩一凸性）。

作者指出，一个可行的 Truesdell “主要问题”解决方案可能需要结合这两者。然而，他们明确表示，迄今为止尚未发现任何在全局范围内同时满足这两个约束的各向同性应变能函数。构建此类函数仍是一个开放性问题。本文认为，这两个属性的推导出的充分条件在全局设定下可能是互斥的，但在受限（链受限）定义域内可以达成一致。这项工作为构成建模提供了一个警示，特别是对于数据驱动的方法（如神经网络），在这些方法中通常将多凸性作为唯一的稳定性约束，从而可能忽略了单轴拉伸中单调应力轨迹这一物理要求。