🔢 mathematics

A finite-strain logarithmic viscoelastic model for Antarctic ice shelves based on an additive split

本文提出了一种用于南极冰架的有限应变对数粘弹性模型，该模型通过采用亨尼（Hencky）应变率的加法拆分来将格伦（Glen）流变定律与各向同性弹性相结合，并通过有限元模拟展示了粘弹性与前缘形态如何共同支配冰架末端的张力。

原作者： Maxime Nutte, Sebastian Skatulla, Carlo Sansour

发布于 2026-09-09

📖 1 分钟阅读🧠 深度阅读

原作者： Maxime Nutte, Sebastian Skatulla, Carlo Sansour

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

南极洲巨大冰盖的漂浮边缘并非静态的冰墙；它们是动态的、会呼吸的结构，不断地发生位移、弯曲并最终破碎。这些被称为冰架的漂浮延伸部分起着至关重要的支柱作用，支撑着后方庞大的冰川。当它们破碎时，会发生被称为“崩解”的过程，将冰山释放到海洋中，这是这些冰盖流失质量的主要方式。几十年来，科学家们一直试图精确预测这种破碎发生的具体时间和地点。为此，他们依赖计算机模型来模拟冰在自身重量和海洋压力下的移动与变形。然而，这些模型传统上将冰视为两种形式之一：要么是受到挤压后会立即弹回的固体（弹性），要么是多年间缓慢流动的厚重流体（粘性）。实际上，冰架同时经历着两种时间尺度的作用力。潮汐和海浪在数小时内推动或拉扯冰层，而冰川的缓慢蠕动则发生在数十年间。一个忽略了冰的快速弹跳效应的模型可能会错过风暴期间形成的裂缝，而一个忽略了缓慢流动效应的模型则可能无法预测冰架的长期形状。

为了解决这个问题，研究团队开发了一种新的模拟冰架行为的方法，该方法能够同时捕捉快速的弹跳和缓慢的流动，即使在冰的形状发生显著变化时也是如此。他们创建了一个数学框架，将冰视为一种可以进行复杂拉伸和弯曲，但仍能记住其原始形状的材料。这种方法使他们能够追踪冰对涨潮等即时应力的反应，同时兼顾数个潮汐周期之间发生的多年缓慢拉伸。通过专注于冰架边缘（即冰与开阔海洋交界处）的具体几何形状，他们发现冰崖的形状以及冰层内部的排列方式，对于决定哪里最容易产生裂缝起着决定性作用。他们的工作表明，冰的弹回能力与缓慢流动倾向之间的相互作用，会在前端产生特定的张力模式，从而控制冰块如何以及何时脱落。

研究人员基于一个称为“麦克斯韦”（Maxwell）系统的概念构建了他们的模型，这是描述既具有弹簧特性又具有粘性流体特性的材料的标准方式。想象一下，一个弹簧与一个在蜂蜜中移动的活塞串联在一起；如果你拉动它们，弹簧会立即拉伸，而蜂蜜则会缓慢地拖动。在他们的模拟中，弹簧代表了冰的即时弹性响应，而蜂蜜则代表了遵循格伦流变定律（Glen's flow law）的缓慢粘性流动，该定律长期以来一直被用于描述冰在数个世纪中的运动。这项研究的创新之处在于他们如何处理冰发生大幅度拉伸时的数学计算。以前的模型在冰发生显著形状改变时往往难以处理，需要运行复杂的计算，这在计算机上非常困难。这种新方法使用了一种称为“对数应变”（logarithmic strain）的特定测量方式，它允许团队以一种直观的、加法性的方式将冰的总运动分解为弹性部分和粘性部分。这意味着他们可以将两种类型的运动相加，而无需追踪一个单独的、不可见的“中间”冰形，这使得他们的计算更加高效，并且更容易集成到现有的气候模型中。

为了测试这种新方法是否有效，团队首先对一个承受自身重量的垂直冰柱进行了模拟，并将其结果与一个成熟且更复杂的模型进行了对比。结果几乎完全吻合。在经过一年半的模拟时间后，新模型预测冰会下沉并向外凸出，其偏差与复杂模型相比仅为几厘米。冰内部的应力水平也几乎完全一致，在大多数区域的差异小于一千帕斯卡。这证实了他们简化的加法方法可以在不增加沉重计算成本的情况下，准确地重现冰的行为。然而，研究人员也谨慎地指出其方法可能失效的地方。他们发现，如果冰受到一种特定的扭转运动（即简单剪切，类似于扑克牌相互滑动）的影响，他们的模型会开始与更复杂的模型产生分歧。在这些极端的扭转场景中，新模型会预测冰的卸载或松弛过快。幸运的是，研究人员确定他们所研究的冰架并不会经历这种极端的扭转。相反，冰架前缘主要是受到重力和浮力的拉伸和弯曲，在这种类型的运动中，他们的新模型保持了高度的准确性。

在验证了模型之后，团队将其应用于一个理想化的冰舌（即一条漂浮在海洋上的长条状冰带），以观察不同的形状和内部属性如何影响前端的应力。他们加入了现实的细节，例如由于压缩作用，冰的底部比顶部密度更大，以及冰在温度较高时流动更容易。他们还测试了不同的冰崖形状，包括带有淹没“足部”的形状和带有凹陷基座的形状。模拟显示，冰崖的形状会剧烈改变张力集中的位置。一个淹没的足部实际上可以反转弯曲的方向，减轻表面的应力，从而可能使冰体更稳定。相反，一个凹陷的面会集中张力，并将最危险的应力点向内陆移动。研究还表明，通过冰层厚度的温度和密度变化会产生自身的内部弯曲力，这种力量独立于冰崖的形状。这些发现强调了崩解的风险不仅取决于冰的总重量，还在于由冰的几何形状和内部结构所产生的精确的力量平衡。

这项研究的最终目标是为预测南极洲冰架的未来提供一个更可靠的工具。通过捕捉对潮汐和波浪的快速弹性响应以及数年间的缓慢粘性蠕动，这种新模型为导致崩解的应力场提供了更清晰的图景。研究人员强调，虽然他们的方法并不是所有类型冰建模的通用替代品，但它特别适用于冰架边缘的具体情况。它允许科学家在模拟数十年间冰体形状的复杂有限变化时，不会失去观察短期事件即时影响的能力。这种能力对于理解冰架如何应对变化的气候至关重要，因为风暴频率的变化和海洋变暖的速度可能会改变维持这些庞大结构稳定的力量平衡。这项工作并不声称解决了崩解之谜，但它提供了一种稳健、高效且符合物理规律的方式来探索冰前缘的力学机制，为驱动冰块破碎的力量提供了更清晰的视角。

技术摘要：一种用于南极冰架的有限应变对数粘弹性模型

问题陈述
南极冰架的质量损失主要通过崩解（calving）过程实现，该过程受冰前缘附近的应力场控制。该应力场跨越了从弹性挠曲（秒至天）到粘性蠕变（年至十年）的时间尺度。虽然大规模冰盖模型通常将冰视为受格伦流变定律（Glen's flow law）支配的粘性流体，但这种方法忽略了准确模拟潮汐挠曲、海洋涌浪以及崩解事件后应力重分布所需的瞬时弹性响应。相反，纯弹性模型无法捕捉长期的蠕变过程。

现有的物理描述通常依赖于麦克斯韦（Maxwell）模型来结合弹性和粘性。然而，标准公式在应用于数年乃至数十年的冰架尺度时面临挑战。虽然冰中的弹性应变保持在较小水平（数量级为 $10^{-5}$ 至 $10^{-4}$ ），但累积的粘性拉伸和不断变化的冰前缘几何形状构成了有限变形。先前的有限应变麦克斯韦模型（例如 Schröder 等人，2026 年提出的乘法分解形式 $\mathbf{F} = \mathbf{F}^e \mathbf{F}^v$ ）需要涉及通过中间构型的张量传输的复杂算法切线。此外，线性应变假设在几年后就会失效，导致其预测的应力和应变极大值与有限应变模型相比出现显著差异（十年后差异可达 30%）。

方法论
作者提出了一种基于对数应变空间而非变形梯度乘法分解的加法分解应变率的有限应变粘弹性模型。

运动学与功共轭性：
该模型利用定义在固定参考构型上的 Hencky（对数）应变 $\boldsymbol{\alpha} = \ln \mathbf{U} = \frac{1}{2} \ln \mathbf{C}$ 。该应变率 $\dot{\boldsymbol{\alpha}}$ 被加法分解为弹性部分和粘性部分：
$\dot{\boldsymbol{\alpha}} = \dot{\boldsymbol{\alpha}}^e + \dot{\boldsymbol{\alpha}}^v$
这种分解避免了对中间构型的需求或对客观应力率的要求。由此导出的功共轭应力 $\boldsymbol{\Sigma}_{\ln}$ 满足内能功率 $\mathcal{W} = \boldsymbol{\Sigma}_{\ln} : \dot{\boldsymbol{\alpha}}$ 。
本构关系：
- 弹性分支： 建模为各向同性的 Hencky 弹性，通过 Lamé 模量（ $\lambda, \mu$ ）将对数应力率与对数应变率联系起来。
- 粘性分支： 建模为一种 Glen 型幂律。粘性应变率定义为 $\dot{\boldsymbol{\alpha}}^v = A q^{m-1} \text{dev} \boldsymbol{\Sigma}_{\ln}$ ，其中 $q$ 是偏应力的第二不变量。该公式确保了粘性流是等容的（ $\text{tr} \dot{\boldsymbol{\alpha}}^v = 0$ ），并在无穷小应变极限下精确还原为 Glen 流变定律。
时间积分：
模型采用中点评估法处理运动学量，并对 Maxwell 演化采用隐式后向欧拉方案。该积分将局部本构更新简化为关于等效偏应力 ( $q$ ) 的单个标量非线性方程。对于标准的 Glen 指数 ( $m=3$ )，这产生了一个可通过卡尔达诺公式（Cardano's formula）求解的降次三次方程；对于线性粘性 ( $m=1$ )，则得到闭式解。系统地推导了算法切线算子，以确保在有限元求解器中的二次收敛。
适用范围：
作者明确研究了加法对数分解的局限性。他们指出，在均匀简单剪切（非共轴路径）下，对数速率会随着剪切量的增加而减小，可能导致相对于乘法模型出现卸载伪影。然而，他们认为冰架变形主要由拉伸和挠曲（共轴或近共轴路径）主导，在这些路径下，加法和乘法公式保持高度一致。

关键结果

基准测试： 该模型通过线性麦克斯韦粘弹性 ( $m=1$ $m = 1$ ) 在自加权柱体基准测试（Schröder 等人，2026 年）中得到了验证。加法对数形式高保真地重现了乘法参考解：
- 顶部沉降和侧向隆起差异小于 3 cm。
- 控制点的冯·米塞斯应力差异仅为 0.3 kPa。
- 差异仅限于几何奇异点（固定角），这些位置的应力集中对边界条件极其敏感。
冰架应用： 该公式被应用于一个具有随深度变化的密度、随温度变化的流动性和各种悬崖几何形状（包括淹没足部和基底凹陷）的理想化冰舌。
- 结果展示了粘弹性与前缘形态如何控制末端的张力。
- 模型成功捕捉到了从弹性挠曲到粘性蠕变的转变，表明悬崖几何形状（例如“游离足部”机制）可以改变挠曲方向并降低表面张力极大值。

意义与主张
本文声称，所提出的加法对数分解为冰川学应用提供了一种计算高效且物理一致的替代方案，可替代乘法有限应变粘弹性。其主要意义在于：

速率兼容性： 由于模型完全基于速率进行构建，它能自然地与现有的基于 Glen 定律求解速度的冰盖求解器耦合，而无需引入单独的粘性变形梯度变量。
简洁性： 加法结构简化了向各向异性行为的扩展，并避免了乘法分解中繁琐的张量传输过程。
适用性： 作者断言，虽然该模型在存在持续大规模剪切的区域（如冰流剪切边缘）不能作为乘法粘弹性的通用替代品，但在冰架和崩解前缘特征性的拉伸和挠曲变形中非常有效。它允许在保留冰典型的小弹性应变的同时，模拟长达数十年的有限几何变化。

这项工作得出结论，该公式能够研究涵盖从潮汐周期到多年维持期在内的全频谱相关时间尺度的近前缘应力场和崩解力学，且不会产生复杂乘法方案带来的计算开销。

技术摘要：一种用于南极冰架的有限应变对数粘弹性模型

类似论文