Caveats on formulating finite elasto-plasticity in curvilinear coordinates

本文针对有限弹塑性在曲线坐标系下的数值实现，通过显式基底变换而非微分几何框架，阐明了轴对称大变形中变形梯度、雅可比行列式及移位算子等关键量的处理难点，并提出了一套适用于有限元一致线性化的实用步骤，以解决非笛卡尔基底引入的额外项及非弹性效应带来的复杂性。

要点

这篇论文就像是在教工程师们如何**“在弯曲的表面上画直线”**，同时确保他们画的图在数学上是完全准确的。

为了让你更容易理解，我们可以把这篇论文的核心内容想象成**“在一个巨大的、会变形的气球上玩橡皮泥”**的故事。

1. 背景：为什么要玩这个游戏？

在工程界（比如造大坝、挖隧道、做心脏支架），我们经常需要模拟物体在受力后发生的巨大变形。

• 通常的做法（笛卡尔坐标系）： 就像在一张平整的方格纸上画图。所有的线都是直的，格子都是正方形。在这种纸上，计算非常简单，就像做小学数学题一样。
• 现实的问题（曲线坐标系）： 但很多物体是圆柱形或球形的（比如隧道、储油罐、血管）。在圆柱体上画图，就像在卷起来的纸或者气球表面上画线。这里的“格子”不再是正方形，而是随着位置变化，有的变长，有的变宽。

论文的核心痛点：
很多工程师想偷懒，直接把“方格纸”上的公式套用到“卷纸”上。但这就像试图用直尺去量地球的周长，虽然看起来差不多，但在大变形（比如把气球吹得很大）时，误差会像滚雪球一样越来越大，导致计算结果完全错误，甚至让模拟程序崩溃。

2. 核心发现：三个被忽视的“隐形助手”

论文通过一个“吹气球”（空腔膨胀）的简单实验，指出了在弯曲坐标下，有三个关键概念必须小心处理，否则就会出错：

A. 变形梯度（Deformation Gradient）：不仅仅是“拉伸”

• 通俗解释： 想象你在气球上画了一个小正方形。当你吹大气球时，这个正方形被拉长了。
• 陷阱： 在平纸上，你只需要看它变长了多少。但在气球上，除了变长，方向也变了，而且格子的密度也变了。
• 比喻： 就像你在地图上画路。在平地上，路变长就是距离增加。但在地球仪上，如果你沿着经线走，路变长的同时，你离北极的距离也在变，这会影响你对“方向”的判断。论文指出，如果不修正这种“方向”和“格子密度”的变化，算出来的变形就是错的。

B. 雅可比行列式（Jacobian）：体积的“隐形税”

• 通俗解释： 这是用来计算物体变形后体积变了多少的。
• 陷阱： 在平纸上，体积变化就是长×宽×高。但在圆柱体上，靠近中心的地方和靠近边缘的地方，同样的“长度变化”代表的实际体积是不一样的（就像切蛋糕，靠近圆心的那一小块和靠近边缘的一大块，虽然角度一样，但面积不同）。
• 比喻： 这就像**“汇率”**。在平地上，1 美元就是 1 美元。但在弯曲的坐标里，1 个“单位长度”在中心可能只值 0.5 个“真实体积”，在边缘可能值 2 个。论文强调，必须时刻记得加上这个“汇率转换”，否则算出来的体积（比如气球里的气量）就是错的。

C. 移位器（Shifter）：连接过去与现在的“翻译官”

• 通俗解释： 这是一个数学工具，用来把物体变形前的状态（参考系）和变形后的状态（当前系）联系起来。
• 陷阱： 在平纸上，过去和现在的方向是一样的，不需要翻译。但在圆柱体上，变形前指向“东”的线，变形后可能指向了“东北”。
• 比喻： 就像翻译官。如果你直接拿中文（变形前）去对英文（变形后），意思就乱了。移位器就是那个确保“东”还是“东”，“北”还是“北”的翻译官。论文指出，如果忽略了这个翻译官，所有的应力计算都会南辕北辙。

3. 弹塑性：当橡皮泥变成“记忆金属”

论文还讨论了一个更复杂的情况：弹塑性。

• 弹性（Elastic）： 像橡皮筋，拉长了松手就缩回去。
• 塑性（Plastic）： 像捏橡皮泥，捏扁了就回不去了。
• 难点： 当物体既发生弹性变形又发生塑性变形时，体积的变化变得很复杂。论文提出了一种聪明的方法，把体积变化拆分成“可恢复的”和“不可恢复的”两部分，并给出了具体的计算公式。这就像在计算你钱包里的钱时，要分清哪些是暂时借来的（弹性），哪些是真正花掉的（塑性），不能混为一谈。

4. 最终成果：更聪明、更省钱的“模拟软件”

这篇论文不仅仅是理论探讨，它的最终目的是教工程师怎么写代码。

• 以前的做法： 为了模拟一个圆柱体，工程师不得不把圆柱体切成无数个小方块（3D 模拟），计算量巨大，电脑跑起来像老牛拉破车。
• 论文的方法： 利用圆柱体的对称性，只算一个“切片”（2D 模拟），但通过上述的“修正公式”（移位器、雅可比等），让这 2D 切片能完美代表整个 3D 圆柱体。
• 结果： 就像**“用一张纸的算力，算出了整个地球的数据”**。论文通过对比实验证明，这种新方法算出来的结果和笨重的 3D 模拟几乎一模一样，但速度快得多，省资源得多。

总结

这篇论文就像是一本**“弯曲世界生存指南”**。它告诉工程师们：

“别以为在圆柱体上算东西和在平地上一样简单！如果你忽略了格子的变形、体积的汇率和方向的翻译，你的计算就会像在没有指南针的迷宫里乱撞。只要按照我们给出的步骤，你就能用更少的电脑资源，算出更精准的工程结果。”

这对于设计更安全的隧道、更耐用的轮胎、甚至更精准的人体血管模型，都有着非常重要的实际意义。

技术摘要

这是一份关于论文《Caveats on formulating finite elasto-plasticity in curvilinear coordinates》（有限弹塑性在曲线坐标系中 formulated 的注意事项）的详细技术总结。

1. 研究背景与问题 (Problem)

• 核心挑战：连续介质力学通常基于张量分析，具有坐标不变性。然而，有限元（FE）数值实现依赖于用户选择的基底下的矩阵表示。当使用非笛卡尔、非正交归一的曲线坐标系（如圆柱坐标系）时，如果直接套用笛卡尔坐标系的公式，会导致严重的定义错误。
• 具体痛点：
  • 大变形（有限应变）：在涉及大变形的问题中（如空腔扩张、金属成型），轴对称模型能显著降低计算成本，但曲线坐标系下的运动学量（如变形梯度、雅可比行列式、移位张量）定义变得复杂。
  • 弹塑性耦合：当引入非弹性效应（如塑性变形）时，必须区分弹性部分和非弹性部分。在曲线坐标系下，柯西应力（Cauchy stress）不仅取决于当前的弹性状态，还取决于构型变化，这使得一致线性化（consistent linearisation）变得极其复杂。
  • 常见误区：许多文献或代码在将笛卡尔公式“ naive（天真）”地扩展到曲线坐标时，忽略了度量系数（metric coefficients）和基矢量的变化，导致体积变换、应力计算和平衡方程出现根本性错误。
• 研究目标：澄清在轴对称有限应变弹塑性问题中，变形梯度、雅可比行列式和移位张量的正确定义与角色，并提供一套稳健的有限元实现步骤。

2. 方法论 (Methodology)

本文并未采用复杂的黎曼流形微分几何框架，而是坚持在标准笛卡尔表示的基础上，通过显式的**基底变换（change of basis）**来推导曲线坐标形式，以确保工程实现的实用性。

• 运动学基础：

  • 变形梯度 ($\mathbf{F}$)：明确区分了笛卡尔分量与曲线分量。指出在圆柱坐标系下，变形梯度的面外分量（out-of-plane component）并不总是 1，且其行列式（体积比）不能直接由曲线分量计算，必须引入度量系数的行列式修正。
  • 移位张量 (Shifter, $\mathbf{S}$)：引入了移位张量来处理参考构型与当前构型之间基矢量的映射。在笛卡尔坐标系中它通常被忽略（为单位矩阵），但在曲线坐标系中，它是连接不同构型张量分量的关键。
  • 雅可比行列式 ($J$)：推导了曲线坐标系下的体积变换公式，指出 $J$ 的计算必须包含当前和参考构型的度量系数行列式。
  • 弹塑性分解：基于 Bilby-Kröner-Lee 的多重分解理论 ($\mathbf{F} = \mathbf{F}^e \mathbf{F}^p$)，定义了中间无应力构型。针对各向同性弹塑性材料，提出了弹性雅可比 ($J_e$) 和塑性雅可比 ($J_p$) 的分离计算方法，特别是利用对数应变的迹来线性化 $J_e$，以解决中间构型度量未知的难题。

• 本构与平衡方程：

  • 采用了基于 Eshelby-zeta 应力 ($\zeta = J_e \sigma$) 的本构关系，该应力测度能更好地保持中间无应力构型，并满足热力学第二定律（Clausius-Duhem 不等式）。
  • 推导了曲线坐标系下的强形式（强平衡方程）和弱形式（虚功原理），并特别处理了散度算子中的克里斯托费尔符号（Christoffel symbols）。

• 轴对称简化：

  • 将三维弱形式特化为轴对称二维问题。详细说明了如何将体积积分和表面积分从三维域映射到 R-Z 平面，并处理了面外分量（$\Theta, \theta$）对变形梯度和位移的缩放影响。

• 数值实现：

  • 提供了更新拉格朗日（UL）和总拉格朗日（TL）两种格式的算法流程。
  • 详细列出了在牛顿 - 拉夫逊（Newton-Raphson）迭代中需要更新的量，包括形状函数的协变导数、克里斯托费尔符号以及变形梯度的逆和转置。

3. 关键贡献 (Key Contributions)

1. 澄清了曲线坐标系下的运动学陷阱：

   • 明确指出在圆柱坐标系下，变形梯度的面外分量（$F^{\theta}_{\Theta}$）在轴对称变形下等于 1，但这并不意味着体积比可以直接由分量乘积得出。
   • 揭示了直接套用笛卡尔公式计算雅可比行列式会导致体积变换不守恒的错误。

2. 提出了稳健的线性化方案：

   • 针对弹塑性大变形，提出了一种基于对数弹性应变迹（$\text{tr}(\epsilon^e)$）的雅可比线性化方法（$J_e = \exp(\text{tr}(\epsilon^e))$），避免了直接对复杂的中间构型度量求导，显著简化了切线刚度矩阵的推导。

3. 构建了完整的轴对称有限元框架：

   • 从运动学定义、应力测度、平衡方程到离散化，提供了一套完整的、自洽的轴对称有限应变弹塑性有限元实现指南。
   • 明确区分了移位张量在平面内和面外分量的不同作用，解决了位移场在曲线坐标下的定义歧义。

4. 算法流程的透明化：

   • 提供了详细的算法伪代码（Algorithm 1 & 2），对比了 UL 和 TL 格式在曲线坐标下的具体实现差异（如 Christoffel 符号的更新时机），为开发者提供了可直接参考的代码逻辑。

4. 数值算例与结果 (Results)

• 算例设置：

  • 对象：厚壁圆筒的空腔扩张问题（Cavity expansion）。
  • 材料：解耦的 Neo-Hookean 超弹性材料，结合修正剑桥模型（Modified Cam-Clay）的塑性硬化。
  • 对比：将本文提出的轴对称二维公式与标准的三维有限元分析进行对比。

• 主要发现：

  • 精度验证：在负柯西压力、偏应力不变量、弹性/塑性雅可比行列式等关键指标上，轴对称公式的结果与三维分析结果高度吻合（差异极小，主要源于三维网格离散化带来的几何近似误差）。
  • 收敛性：牛顿 - 拉夫逊迭代在两种格式下表现出相似的收敛行为（平均迭代次数约为 3.5 次），证明了所提出的一致线性化方案是正确且稳健的。
  • 计算效率：轴对称方法仅需二维离散化，显著降低了计算成本，同时保持了三维分析的精度。

5. 意义与结论 (Significance & Conclusion)

• 工程应用价值：该论文为岩土工程（如桩基、隧道）、金属成型（如挤压、拉拔）和生物力学（如动脉生长）中的轴对称大变形问题提供了可靠的数值模拟基础。
• 理论修正：纠正了以往在曲线坐标系下处理有限应变时可能存在的“直觉性”错误，强调了移位张量和度量系数在张量分量转换中的核心作用。
• 指导意义：论文强调，任何对底层运动学量定义的错误都会导致后续导数计算和应力测度的连锁错误，进而导致更新拉格朗日和总拉格朗日格式之间的不一致。因此，在开发有限元代码时，必须严格遵循曲线坐标下的张量变换规则。

总结：这是一篇兼具理论深度与工程实用性的论文，它填补了有限应变弹塑性在轴对称曲线坐标系下数值实现细节的空白，确保了此类高非线性问题在计算力学中的准确性和鲁棒性。