try ai
科普
编辑
分享
反馈
  • 亚弹性模型

亚弹性模型

SciencePedia玻尔百科
关键要点
  • 亚弹性模型通过应力率和变形率之间的直接关系来定义材料的行为。
  • 该模型的一个朴素公式不具备客观性,需要使用共旋或客观应力率(如 Jaumann 率)来正确解释刚体旋转。
  • 与超弹性模型不同,亚弹性模型通常是路径依赖的,这意味着它可以在一个闭合的变形循环中虚假地产生或耗散能量,导致不符合物理规律的预测。
  • 尽管在纯弹性方面存在理论缺陷,但该框架对于解释某些非线性效应很有价值,并可作为通向塑性模型的计算桥梁。

引言

描述材料如何响应力的核心在于应力与应变之间的关系。亚弹性模型提供了最直观的起点之一:如果应力变化的速率与材料变形的速率成正比,会怎样?这个简单而优雅的想法构成了一种基于率的本构律,似乎是线性弹性在大变形情况下的自然延伸。然而,这种表面的简单性掩盖了深刻的物理和数学挑战,揭示了直观的率定律与热力学和客观性基本原理之间的关键知识鸿沟。

本文将引导您了解亚弹性模型的迷人故事。它首先考察其核心原理和机制,揭示其初始公式中的关键缺陷以及解决该问题的巧妙数学修正。您将了解到为什么这种修正虽然解决了一个问题,却暴露了更深层次、更麻烦的路径依赖性病态,这正是亚弹性与真弹性之间的区别所在。随后,我们将探讨该模型的应用和跨学科联系,展示它尽管存在局限性,却如何为复杂的材料行为提供宝贵的见解,提出独特的计算挑战,并作为通向塑性力学和损伤力学等更高级理论的重要桥梁。

原理与机制

想象一下,您正试图描述一块橡胶如何抵抗变形。最简单的想法,即弹簧胡克定律的美妙推广,是说材料的抵抗力——其应力——应与其拉伸量成正比。对于一个不仅会拉伸,还会流动、剪切和翻滚的复杂三维世界,我们或许可以完善这个想法:也许应力的变化率与拉伸率成正比。这便是​​亚弹性模型​​得以发展的那个极其直观的种子。这是一段始于优雅简洁,遭遇深刻悖论,并最终引导我们更深入理解材料本质的旅程。

一个看似简单的开端:应力率与应变率的相遇

让我们尝试将最初的直觉形式化。在连续介质力学中,我们有两个关键角色。首先是​​柯西应力​​张量 σ\boldsymbol{\sigma}σ,它是我们对材料内部任意点内力或“应力”的精密量度。其次是​​变形率​​张量 D\boldsymbol{D}D,它描述了该点材料如何拉伸或改变形状。它捕捉了应变的本质,并与纯刚体运动分离。

连接这两者最直接的方式是提出一个线性关系。最简单的猜测是,应力的物质时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙ 与变形率 D\boldsymbol{D}D 成正比。对于在所有方向上表现相同的材料(各向同性),数学中的表示定理告诉我们,这种关系必须采用我们熟悉的形式:

σ∘=2μD+λtr⁡(D)I\overset{\circ}{\boldsymbol{\sigma}} = 2\mu \boldsymbol{D} + \lambda \operatorname{tr}(\boldsymbol{D}) \boldsymbol{I}σ∘=2μD+λtr(D)I

此处,μ\muμ 和 λ\lambdaλ 是类似于线性弹性中拉梅参数的常数,代表材料的刚度,而 I\boldsymbol{I}I 是单位张量。目前,我们暂时将应力率 σ∘\overset{\circ}{\boldsymbol{\sigma}}σ∘ 等同于简单的时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙。这个方程 是亚弹性的核心。这是一个基于率的定律:告诉我你现在如何变形,我就会告诉你应力现在如何变化。这似乎是一个完美、符合常识的模型。但物理学中常有的情况是,当事物开始旋转时,常识可能是一个危险的向导。

翻滚带来的麻烦:物理学家的头痛问题

物理学中最基本的原理之一是,描述材料行为的定律不应依赖于观察者自身的运动。无论您是静止站立还是在旋转木马上旋转,材料的内在属性必须保持不变。这就是​​材料坐标系无关性原理​​,或称​​客观性​​。任何有效的本构律都必须遵守它。让我们来检验一下我们简单的亚弹性定律。

想象一下,我们取一块已经处于某种应力状态下的橡胶块——比如说,它已经被预压缩。现在,我们只是将整个块体进行刚性旋转,不施加任何额外的拉伸或挤压。由于没有变形,变形率张量 D\boldsymbol{D}D 为零。我们提出的定律 σ˙=C:D\dot{\boldsymbol{\sigma}} = \mathbb{C}:\boldsymbol{D}σ˙=C:D 将预测 σ˙=0\dot{\boldsymbol{\sigma}} = \boldsymbol{0}σ˙=0。这意味着,在我们固定的实验室坐标系中测量的应力张量分量不应改变。

但这显然是错误的!当块体旋转时,内力也随之旋转。压缩和拉伸的方向相对于我们实验室的坐标轴正在改变。应力张量 σ\boldsymbol{\sigma}σ 的分量必须改变以反映这种新的方向。我们简单的定律在明显需要变化时却预测没有变化,这一事实告诉我们,我们犯了一个根本性的错误。朴素的物质时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙ 不是一个客观量。它被材料的刚性翻滚所“污染”,未能区分材料应力状态的真实变化与该状态在空间中的纯粹重新定向。

共旋修正:与材料一同运动

为了挽救我们优美的基于率的想法,我们需要一种能够忽略纯旋转的应力率测量方法。解决方案是定义一种新的导数,这种导数是从一个与该点材料一同旋转的假想观察者的角度来计算的。这被称为​​共旋率​​,是​​客观应力率​​的一个例子。

为此,我们需要知道材料是如何旋转的。这个信息由​​自旋张量​​ W\boldsymbol{W}W 捕捉,它是速度梯度的反对称(或旋转)部分。然后,我们可以通过减去由这种自旋产生的部分来“修正”我们朴素的时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙。最常见的方法之一是定义​​Jaumann 率​​ σ˙J\dot{\boldsymbol{\sigma}}^{\mathrm{J}}σ˙J:

σ˙J=σ˙−Wσ+σW\dot{\boldsymbol{\sigma}}^{\mathrm{J}} = \dot{\boldsymbol{\sigma}} - \boldsymbol{W}\boldsymbol{\sigma} + \boldsymbol{\sigma}\boldsymbol{W}σ˙J=σ˙−Wσ+σW

项 (Wσ−σW)(\boldsymbol{W}\boldsymbol{\sigma} - \boldsymbol{\sigma}\boldsymbol{W})(Wσ−σW) 精确地解释了由于材料以自旋 W\boldsymbol{W}W 旋转而引起的应力分量的变化率。通过减去它,我们分离出了由实际变形引起的应力变化率。

现在,我们的亚弹性定律变为 σ˙J=2μD+λtr⁡(D)I\dot{\boldsymbol{\sigma}}^{\mathrm{J}} = 2\mu \boldsymbol{D} + \lambda \operatorname{tr}(\boldsymbol{D}) \boldsymbol{I}σ˙J=2μD+λtr(D)I。让我们重新审视我们的旋转块体实验。对于刚性旋转,D=0\boldsymbol{D}=\boldsymbol{0}D=0,因此我们的新定律正确地预测 Jaumann 率为零:σ˙J=0\dot{\boldsymbol{\sigma}}^{\mathrm{J}} = \boldsymbol{0}σ˙J=0。这意味着对于我们那个与材料一同运动的微小观察者来说,应力看起来是恒定的——这完全正确,因为没有产生新的应力。我们的定律现在是客观的了。危机似乎已经避免。

路径依赖的迷宫:亚弹性模型的黑暗秘密

共旋修正虽然在数学上是合理的,但却隐藏着一个深刻而令人不安的物理病态。真弹性的根本在于能量的储存和释放。如果你使一个完美的弹性物体变形,然后让它恢复到原始形状,它会归还你投入的全部能量。所做的功仅取决于最终状态,而与达到该状态所采取的路径无关。这意味着存在一个​​应变能势​​,拥有这种势的材料被称为​​超弹性​​材料。

我们的客观亚弹性模型描述的是这样的材料吗?它的率方程是否“可积”成一个能量势?令人震惊的答案是,对于大多数标准客观率(如 Jaumann 率)而言,答案是​​否​​。该模型在本质上是​​路径依赖​​的。

一个经典且极具说服力的证据来自​​单剪切​​的思想实验——想象一下将一本书的顶盖相对于底盖滑动。如果我们以恒定速率施加这种剪切运动,并要求我们的 Jaumann 率亚弹性模型预测剪切应力,它会给出一个奇怪的结果。它预测剪切应力不会稳定增长,而是会随着剪切的无限增加而像正弦波一样振荡、上升和下降。这仿佛是材料在你持续向同一方向剪切它时,魔术般地变弱,然后变强,再变弱。这在物理上是荒谬的,并且从未在简单的弹性固体中观察到。

这种不可积性意味着,如果我们能将一块这种假想的材料进行一次变形的“往返旅行”——比如,剪切它,旋转它,取消剪切,再取消旋转,使其精确地回到起始构型——应力状态可能不会恢复到初始状态。在这个闭合的弹性循环中所做的净功可能不为零。该材料要么创造了能量,要么销毁了能量,这对于一个纯弹性模型来说是致命的缺陷。

出路:超弹性与计算视角

亚弹性模型的病态使物理学家和工程师们认识到,虽然基于率的方法很直观,但它作为弹性的基础是存在缺陷的。模拟大弹性变形的正确方法是使用​​超弹性​​框架。超弹性不是为率假设一个定律,而是从一个标量​​应变能函数​​ WWW 开始,该函数定义了材料作为其变形状态(例如,作为变形梯度 F\boldsymbol{F}F 的函数)的函数所储存的能量。然后,应力就简单地从这个势函数中导出。

这种方法在几个方面具有内在的优越性:

  • ​​热力学一致性:​​ 它通过构造保证了弹性响应是路径无关和能量守恒的。它绝不会在一个闭合的变形回路上虚假地产生能量。
  • ​​内在客观性:​​ 因为能量取决于本身就是客观的变形度量(如右柯西-格林张量 C=FTF\boldsymbol{C}=\boldsymbol{F}^{\mathsf{T}}\boldsymbol{F}C=FTF),所以由此产生的应力-应变关系自动地与坐标系无关,无需明确选择一个客观率。
  • ​​物理真实性:​​ 超弹性模型,如 neo-Hookean 模型,预测在单剪切中剪切应力单调增加,这与实验观察一致,避免了 Jaumann 率模型不符合物理规律的振荡。

从计算的角度来看,差异是天壤之别。使用一个简单的数值方案来模拟一个有限旋转,一个亚弹性模型可能会产生完全错误的、非零的剪切应力,这纯粹是由于算法未能完美地复制模型所要求的旋转更新而产生的伪影。这种数值误差是模型底层物理不一致性的一个症状。相比之下,超弹性公式基于一个势函数,这导致了更稳健和可靠的数值算法,在复杂模拟中通常具有更好的收敛性。

那么,我们为什么还要研究亚弹性模型呢?因为通过其失败的历程,我们学到了弹性到底是什么。它迫使我们面对旋转和客观性的微妙之处,并揭示了为什么基于能量的超弹性描述不仅仅是一种替代方案,而是正确的物理基础。然而,将应力率与应变率关联起来的亚弹性思想并非完全没有价值;它在​​塑性力学​​的世界里找到了真正的用武之地,在那个世界里,材料会流动、耗散能量,而路径依赖性正是故事的全部意义所在。然而,对于弹性而言,亚弹性模型在通往更深层次物理真理的道路上,是一个美丽但最终具有警示意义的故事。

应用与跨学科联系

在经历了亚弹性模型抽象原理的旅程之后,我们现在到达了探索中最激动人心的部分:看到这些思想在现实世界中的应用。一个物理模型的优劣取决于它描述、预测和帮助我们改造周围世界的能力。您会发现,亚弹性框架尽管看似简单,却为我们理解一系列奇特的材料行为打开了一扇门,在计算工程中提出了深刻的挑战和见解,并成为我们迈向最先进物质理论的关键垫脚石。

应力与应变的舞蹈:捕捉无形之力

想象一下搅拌蜂蜜。你会感觉到对圆周运动的阻力——一种剪切应力。一个简单的牛顿流体模型说,这就是全部了。但如果流体在抵抗你搅拌的同时,还试图将勺子向上推出罐子呢?这种现象,即剪切材料会产生垂直于剪切方向的力,被称为“法向应力效应”。这是真实存在的。在实验室里,某些聚合物溶液会著名地沿着浸入其中的旋转杆向上爬升,似乎在挑战重力。

这正是那些更简单的模型无法捕捉,但亚弹性模型可以捕捉到的行为。因为亚弹性模型考虑了应力张量如何随着材料变形而旋转,所以它可以预测这些法向应力从纯剪切流中出现。在一个受控的思想实验中,通过让亚弹性材料经受单剪切,我们可以计算应力分量的演变,并看到这些法向应力差异从方程中自然地出现。这不仅仅是聚合物的一个怪癖;类似复杂的响应在地质力学中也至关重要,在构造运动或滑坡期间,岩石和土层的巨大剪切涉及的应力场远比简单直觉所暗示的要复杂得多。亚弹性模型为我们提供了一种数学语言来描述这些错综复杂的力学对话。

计算的熔炉:在计算机中构建材料

连续介质力学最深远的应用之一在于计算模拟。利用有限元法(FEM)等工具,工程师和科学家们构建并测试从飞机机翼到地质断层等一切事物的虚拟原型。为此,他们需要将材料行为的物理定律编码到计算机中。在这里,亚弹性不仅是一种理论,更是一种实用且具有挑战性的方法。

任何材料模型的一个基本要求是“客观性”。这是一个花哨的术语,用来表达一个简单而深刻的思想:支配材料的物理定律不能取决于你是从静止位置观察它,还是在坐旋转木马时观察它。材料本身不知道也不在乎你如何旋转。它只关心它是否被拉伸、压缩或剪切。将纯旋转与真实变形分离开来的数学工具就是*客观应力率*。

为了确保计算模型的客观性,工程师们会进行严格的测试。一个经典的测试是模拟纯刚体旋转——就像一块材料在空间中旋转而不改变形状——并验证模型预测应力变化为零。使用像 Jaumann 率这样的客观率的亚弹性模型完美地通过了这项测试;它正确地识别出没有变形(D=0D=0D=0),因此除了应力张量本身的预期旋转外,计算出的应力变化为零。缺乏客观率的模型会失败,会无中生有地从纯旋转中产生应力。

然而,这种客观性是有计算代价的。当这些模型在有限元软件中实现时,客观率中考虑自旋的项会产生一个微妙但关键的数学特性:它们使系统的“切线刚度矩阵”变得非对称。直观地说,这意味着材料对一个方向小推动的抵抗力不再保证与对另一个方向类似推动的抵抗力相同。这要求工程师使用更复杂、计算量更大的数值求解器。因此,亚弹性模型的影响从材料的物理特性延伸到了用于模拟它的软件的架构本身,这在模拟滑坡等极端变形问题的高级方法中(如物质点法 (MPM))是一个至关重要的考虑因素。

一个关于完整性的问题:路径依赖的世界

在这里,我们遇到了亚弹性模型一个优美微妙且在哲学上具有深度的方面。想象一下,你拿一根橡皮筋,拉伸它,扭转它,然后小心地反向执行每一个动作,使它精确地回到起始形状。你期望它再次变得没有应力,因为它已经归还了所有储存的能量。这种被称为“路径无关性”的特性,是真正弹性(或“超弹性”)材料的标志。这样的材料对其无应力状态有一个完美的“记忆”,由一个储存能量函数定义。

一般而言,亚弹性模型不具备此特性。它们是“失忆者”。它们的定义纯粹是增量式的——“如果你以这个速率拉伸我,我的应力就以那个速率变化”——而没有参照一个全局的能量图。如果你让一个亚弹性模型在变形空间中进行一次往返旅行——一个闭合回路——它通常会回到起点时带有残余应力,凭空创造或销毁了能量。这种不可积性,或称路径依赖性,是一个主要的理论局限。

这一发现引出了一个有趣的问题:在路径依赖性方面,所有客观率都是平等的吗?通过将不同类型的亚弹性模型(例如,使用 Jaumann 率与 Green-Naghdi 率)与一个“基准真相”的超弹性模型进行比较,我们发现它们并非如此。对于给定的变形路径,某些率产生的结果比其他率更接近能量守恒的理想状态。这表明选择亚弹性模型是一门艺术,是在简单性与物理保真度之间的权衡。

通往现代性之桥:塑性、损伤与稳定性

那么,在现代材料科学的世界里,亚弹性模型处于什么位置呢?它充当了一座至关重要的桥梁,将历史思想与我们最先进的理论联系起来。

例如,在金属塑性领域,工程师们长期以来一直使用直观的模型,其中总变形简单地是弹性部分和塑性(永久)部分的总和。事实证明,对于许多重要情况,比如金属构件的单剪切,在小弹性应变的假设下仔细分析,数学上严谨的亚弹性-塑性框架得到的结果与更简单的工程模型完全相同。这为那些更简单的模型为何如此有效提供了强有力的论证,并表明当其局限性被理解时,亚弹性可以是一个计算上方便的选择。

然而,当我们涉足更复杂的现象,如材料损伤——导致失效的微裂纹的形成和生长——亚弹性的临时特性就成了一个显著的缺点。现代物理学要求,一个关于材料的完整理论,特别是涉及像损伤这样的耗散过程的理论,必须以热力学为基础。这是通过从自由能势构建模型来实现的,这是超弹性的一个标志。这样的框架自动保证了热力学定律得到遵守,并且模型是客观的。在这种背景下,亚弹性被揭示为不是最终目的地,而是通往这些更强大、热力学上一致的理论的重要先驱。

最后,我们之前遇到的非对称切线矩阵还有最后一个、惊人的教训要告诉我们。在动力学中,对称的刚度矩阵会导致实值的振动频率。而非对称的刚度矩阵,正如亚弹性模型所产生的,则为复数频率打开了大门。复数频率对应于一种称为“颤振”的振荡不稳定性——与可能导致飞机机翼灾难性地撕裂自身的同一种不稳定性。因此,材料模型的选择直接影响到结构稳定性的预测。一个更简单的超弹性模型会产生一个对称系统,永远无法预测颤振,而一个亚弹性模型,无论好坏,都包含了使这种动态不稳定性成为可能的数学结构。

从捕捉聚合物奇特的爬升现象,到塑造超级计算机的算法,再到将材料行为与动态稳定性理论联系起来,亚弹性模型提供了一幅丰富而富有启发性的图景。它提醒我们,在科学中,即使是看起来最简单的想法也能引导我们踏上深刻的发现之旅,揭示世界的美丽复杂性以及描述它所需的微妙艺术。