try ai
科普
编辑
分享
反馈
  • 有限应变公式:一种大变形框架

有限应变公式:一种大变形框架

SciencePedia玻尔百科
核心要点
  • 变形梯度 (F\mathbf{F}F) 是有限应变理论中的基本量,它完整描述了材料的局部拉伸、剪切和旋转。
  • 像格林-拉格朗日张量 (E\mathbf{E}E) 这样的应变度量至关重要,因为它们将真实变形与刚体运动分离开来,为本构律提供了基础。
  • 全拉格朗日 (TL) 公式(参考初始状态)和更新拉格朗日 (UL) 公式(参考当前状态)之间的选择,是一个取决于具体问题的关键计算决策。
  • 对于塑性,变形梯度的乘法分解 (F=FeFp\mathbf{F} = \mathbf{F}_e \mathbf{F}_pF=Fe​Fp​) 是一个核心概念,它将总变形分解为可恢复的弹性和永久的塑性两部分。

引言

当橡胶等材料拉伸至其原始长度的数倍,或金属被锻造成新形状时,线性力学的简单规则便不再适用。描述这些大变形需要一个更复杂、更强大的框架,即有限应变公式。本文旨在解决超越无穷小假设的挑战,以精确模拟变形世界中复杂的现实情况。它为支配材料经历显著形状和方向变化的原理提供了一份概念性指南。

首先,在“原理与机制”部分,我们将深入探讨基本运动学,定义变形梯度、客观应变张量及其共轭应力伙伴等关键概念。我们将探究这些元素如何构成全拉格朗日和更新拉格朗日公式等计算策略的基础。在这一理论基础之后,“应用与跨学科联系”部分将展示该框架如何应用于解决现实世界的问题,从模拟金属塑性和软组织生物力学,到理解土壤和先进材料的行为。

原理与机制

想象一下观看铁匠锻造一块钢材。炽热的金属在锤下流动,从一个简单的方块变成一个复杂的涡卷形。或者想象一根橡皮筋,拉伸到其原始长度的许多倍然后弹回。描述这种看似混沌的物质之舞是连续介质力学的一个核心挑战。我们无法追踪每一个原子。取而代之,我们进行了一次强大而优雅的想象飞跃:我们假装材料是一个​​连续介质​​,一种光滑无缝的物质织物,我们可以在其每一点上定义密度和温度等属性。这是我们进入有限变形世界的起点。

描述运动:世界并非刚体

我们的首要任务是创建一张运动的地图。假设我们的物体——橡皮筋、钢块——最初占据空间中一个我们称之为​​参考构型​​ B0\mathcal{B}_0B0​ 的区域。我们可以用一个位置矢量 X\mathbf{X}X 来标记这个原始状态下的每一个材料点。随着物体的变形,它移动到一个新的形状,即​​当前构型​​ Bt\mathcal{B}_tBt​。原来位于 X\mathbf{X}X 的材料点现在处于一个新的空间位置 x\mathbf{x}x。整个运动被一个优美的数学函数——​​变形映射​​ φ\boldsymbol{\varphi}φ 所捕捉,它告诉我们每个点去向何方:x=φ(X,t)\mathbf{x} = \boldsymbol{\varphi}(\mathbf{X}, t)x=φ(X,t)。

思考这种变化最简单的方式是​​位移矢量​​ u=x−X\mathbf{u} = \mathbf{x} - \mathbf{X}u=x−X,它就是从一个粒子原始位置指向其新位置的箭头。这告诉我们物体移动到了哪里,但并未揭示全部情况。它没有告诉我们材料本身是如何被拉伸、剪切或旋转的。为此,我们需要一个更强大的工具。

变形梯度:一个局部放大镜

想象一下,在我们未变形的橡皮筋表面上取一个无穷小的矢量,比如一个微小的绘制箭头 dXd\mathbf{X}dX。拉伸后,这个微小的箭头在当前构型中变成一个新的矢量 dxd\mathbf{x}dx。它的长度和指向很可能都发生了变化。那么,原始箭头和新箭头之间是否存在一种简单的关系呢?答案是肯定的,而且非常显著。对于平滑变形,这种关系是线性的。存在一个张量,一种数学机器,称为​​变形梯度​​ F\mathbf{F}F,它将旧矢量转换为新矢量:

dx=FdXd\mathbf{x} = \mathbf{F} d\mathbf{X}dx=FdX

这个 F\mathbf{F}F 定义为运动相对于原始材料坐标的梯度,即 F=∇Xx\mathbf{F} = \nabla_X \mathbf{x}F=∇X​x。它就像一个局部放大镜,包含了关于一个点周围材料邻域如何变形的所有信息。它告诉我们物体被拉伸了多少、剪切了多少,以及旋转了多大的角度。

变形梯度与位移之间存在一个非常简洁而精确的关系。通过对位移定义 x=X+u\mathbf{x} = \mathbf{X} + \mathbf{u}x=X+u 求梯度,我们发现:

F=I+∇Xu\mathbf{F} = \mathbf{I} + \nabla_X \mathbf{u}F=I+∇X​u

其中 I\mathbf{I}I 是单位张量。这并非针对小位移的近似;它是一个纯粹的运动学恒等式,对任何变形都成立,无论多大。它优雅地将总变形描述 F\mathbf{F}F 与更直观的位移场 u\mathbf{u}u 联系起来。

变形梯度具有深刻的物理意义。它的行列式 J=det⁡(F)J = \det(\mathbf{F})J=det(F) 告诉我们局部体积的变化情况。如果 J=1J=1J=1,则运动是保体积的。如果 J>1J > 1J>1,则材料发生了膨胀。而且至关重要的是,我们必须始终有 J>0J > 0J>0。为什么?因为 J=0J=0J=0 意味着一个有限体积被压缩为零,而 J0J 0J0 则意味着材料被由内向外翻转——这些都是我们必须在理论中排除的物理不可能性。

测量应变:它被拉伸了多少?

材料并不会仅仅因为在空间中进行刚性平移或旋转就产生应力。放在旋转转盘上的咖啡杯什么也感觉不到。应力只源于拉伸和剪切——即真实的变形。因此,我们需要一种能够将这种拉伸与刚性旋转分离开来的度量,一种纯粹的​​应变​​度量。

让我们回到那个微小的材料箭头 dXd\mathbf{X}dX。它的原始长度平方是 (dS)2=dX⋅dX(dS)^2 = d\mathbf{X} \cdot d\mathbf{X}(dS)2=dX⋅dX。变形后,它的新长度平方是 (ds)2=dx⋅dx(ds)^2 = d\mathbf{x} \cdot d\mathbf{x}(ds)2=dx⋅dx。利用我们的新工具,我们可以将其写为 (ds)2=(FdX)⋅(FdX)=dX⋅(FTF)dX(ds)^2 = (\mathbf{F} d\mathbf{X}) \cdot (\mathbf{F} d\mathbf{X}) = d\mathbf{X} \cdot (\mathbf{F}^T \mathbf{F}) d\mathbf{X}(ds)2=(FdX)⋅(FdX)=dX⋅(FTF)dX。于是,长度平方的变化为:

(ds)2−(dS)2=dX⋅(FTF−I)dX(ds)^2 - (dS)^2 = d\mathbf{X} \cdot (\mathbf{F}^T \mathbf{F} - \mathbf{I}) d\mathbf{X}(ds)2−(dS)2=dX⋅(FTF−I)dX

这个表达式启发了固体力学中最重要的量之一的定义:​​格林-拉格朗日应变张量​​ E\mathbf{E}E:

E=12(FTF−I)\mathbf{E} = \frac{1}{2}(\mathbf{F}^T \mathbf{F} - \mathbf{I})E=21​(FTF−I)

这个张量堪称精美。如果变形是纯刚性旋转,那么 F\mathbf{F}F 就是一个正交张量 Q\mathbf{Q}Q,满足 QTQ=I\mathbf{Q}^T\mathbf{Q}=\mathbf{I}QTQ=I。在这种情况下,E=12(I−I)=0\mathbf{E} = \frac{1}{2}(\mathbf{I} - \mathbf{I}) = \mathbf{0}E=21​(I−I)=0。格林-拉格朗日应变完全“看不见”刚性旋转;它只“看到”拉伸。此外,它是一个材料量,是相对于原始未变形状态定义的,这使其成为​​全拉格朗日 (TL) 公式​​的完美应变度量——这是一种我们总是返回到原始参考构型进行计算的策略。

当然,物理学不应依赖于我们的观察角度。我们也可以从当前变形状态的角度来定义一个应变度量。这就引出了​​欧拉-阿尔曼西应变张量​​ e=12(I−b−1)\mathbf{e} = \frac{1}{2}(\mathbf{I} - \mathbf{b}^{-1})e=21​(I−b−1),其中 b=FFT\mathbf{b} = \mathbf{F}\mathbf{F}^Tb=FFT 是左柯西-格林张量。这是​​更新拉格朗日 (UL) 公式​​的自然度量,在该公式中,当前状态成为下一步的参考。这两个应变度量 E\mathbf{E}E 和 e\mathbf{e}e 并非竞争对手,而是亲戚,从不同角度描述同一个物理现实。它们通过一个简单的变换联系在一起:e=F−TEF−1\mathbf{e} = \mathbf{F}^{-T} \mathbf{E} \mathbf{F}^{-1}e=F−TEF−1,这优美地展示了底层运动学的统一性。

应力与应变的舞蹈:功与能

应变并非凭空产生;它是材料对应力的响应。应力与应变是一场能量之舞的伙伴。当你拉伸一根橡皮筋时,你对它做功,这个功以势能的形式储存起来。这种联系通过​​功共轭​​的概念被形式化。内功密度由一个应力度量与其共轭应变度量的变分之积给出。

在全拉格朗日公式的材料世界中,优雅的格林-拉格朗日应变 E\mathbf{E}E 与一个同样优雅的应力度量配对:​​第二皮奥拉-基尔霍夫 (PK2) 应力​​ S\mathbf{S}S。功密度为 S:δE\mathbf{S} : \delta\mathbf{E}S:δE。这对组合在描述像橡胶这类被称为​​超弹性​​材料时异常强大。对于这些材料,储存的能量仅仅是应变的函数,W(E)W(\mathbf{E})W(E)。应力则直接从能量函数中得出: S=∂W∂E\mathbf{S} = \frac{\partial W}{\partial \mathbf{E}}S=∂E∂W​ 这是一种简单的、路径无关的关系,这也是TL公式在这类材料中如此普遍的原因之一。

在更新拉格朗日公式的空间世界中,舞台是当前构型。在这里,我们熟悉的​​柯西应力​​ σ\boldsymbol{\sigma}σ(单位当前面积上的力)扮演主角,与变形率 d\mathbf{d}d 共舞。为了完善演员阵容,我们还有​​第一皮奥拉-基尔霍夫 (PK1) 应力​​ P\mathbf{P}P,这是一个连接两个世界的“混合”度量,它将当前构型中的力与参考构型中的面积联系起来。

当事情变得复杂:塑性与率

超弹性理论很美,但当你弯曲一个回形针时会发生什么?它会保持弯曲状态。变形是永久性的。这就是​​塑性​​,为了描述它,我们需要一个更深层次的思想。由 Lee 和 Liu 首次提出的绝妙见解是,想象总变形 F\mathbf{F}F 可以被分解为两个连续的步骤:

F=FeFp\mathbf{F} = \mathbf{F}_e \mathbf{F}_pF=Fe​Fp​

在这里,Fp\mathbf{F}_pFp​ 代表塑性变形——材料微观结构的永久性重排,比如金属中的位错滑移。这个映射将你从原始构型带到一个假想的、局部无应力的​​中间构型​​。然后,Fe\mathbf{F}_eFe​ 代表随后的弹性变形——原子晶格从这个松弛状态到最终受力构型的拉伸。这种乘法分解是现代固体力学中最深刻的概念之一。它使我们能够将运动中的可恢复(弹性)部分与永久(塑性)部分分离开来。对于大多数金属,塑性流动是一个保持体积的剪切过程,这一事实可以通过简单的约束 det⁡(Fp)=1\det(\mathbf{F}_p) = 1det(Fp​)=1 来体现。

然而,当我们考虑时间时,这个框架带来了新的挑战。对于许多材料,应力不仅取决于应变,还取决于应变率。这需要一个率相关的本构律。但什么是“应力率”呢?你可能认为是材料时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙。但这个简单的选择隐藏着一个微妙的缺陷。想象一个已经承受应力的物体受到纯刚性旋转。物理上,应力状态应该只是随物体旋转,而不应产生新的应力。材料时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙ 未能通过这个测试;它错误地预测了虚假的应力。这违反了一个称为​​材料框架无关性原理​​或客观性的基本原则。

解决方案是发明一种新的导数,即​​客观应力率​​,例如 Jaumann 率或 Green-Naghdi 率。它们被巧妙地构造出来,减去由纯旋转引起的应力变化部分,只留下由真实变形引起的应力变化。它们确保了我们的物理定律不依赖于一个旋转的观察者。这种微妙之处是有限应变理论所独有的;在无穷小应变的世界里,旋转效应是高阶项,可以被安全地忽略。

从理论到计算:拉格朗日公式

我们如何将这个优美的理论结构转化为一种实用的工具,用计算机(例如使用有限元法)来解决现实世界的问题?运动学框架的选择导致了两种主要的计算策略。

​​全拉格朗日 (TL) 公式​​将其目光固定在初始的、未变形的状态上。所有关于运动、应变和应力的方程都被“拉回”到原始、不变的计算网格上求解。对于那些初始状态是自然且恒定参考的问题,比如模拟一个橡胶支座(其行为总是相对于其制造形状来定义),这种方法异常优雅和高效。

相比之下,​​更新拉格朗日 (UL) 公式​​则像一个游牧者。它将当前构型作为计算下一个运动增量的临时参考。对于那些当前状态最为重要的问题,这种方法是不可或缺的。想象一下汽车碰撞模拟或金属零件的锻造过程:新的接触面不断形成,压力载荷施加在不断变化的几何体上。UL 公式以其自然的优雅处理这些复杂性。

即使在这里,理论的统一性也熠熠生辉。一个材料的本构律可能在材料坐标系中定义最为自然(例如,关联 S\mathbf{S}S 和 E\mathbf{E}E)。要在一个运行于空间坐标系的 UL 代码中使用这个本构律,我们必须将材料刚度转换为空间刚度。这是通过一个​​前推​​操作完成的,这是一个数学桥梁,确保计算模型与底层物理一致,正确地将材料的内在响应转换到我们观察到的变形世界中。从最简单的运动描述到最先进的计算算法,有限应变理论的原理为理解我们这个可变形世界的力学提供了一个连贯、强大且极具美感的框架。

应用与跨学科联系

在遍历了有限应变运动学的复杂景观之后,你可能会有一种数学上的成就感,但可能还有一个挥之不去的问题:“这一切究竟是为了什么?” 这是一个合理的问题。既然简单的线性理论在很多时候似乎行之有效,我们为什么还要投入如此多的精力来定义变形梯度 F\mathbf{F}F 和各种应变张量等量呢?

答案,以及这个框架的真正美妙之处在于,世界并非线性的。物体以宏大而复杂的方式弯曲、扭转、拉伸和压缩。当变形巨大时,简单的规则就会失效,我们需要一种更强大的语言来描述现实。有限应变原理不仅仅是一个数学上的奇趣;它们是我们理解真实变形世界的基础。它是材料科学家锻造新合金、地球物理学家预测滑坡、生物工程师设计人造组织,乃至纺织工程师创造智能织物时所共同使用的统一语法。现在,让我们来探讨其中一些对话。

金属的世界:成形与断裂

考虑一个普通的回形针。你可以轻微弯曲它,它会弹回——这是我们熟悉的弹性范畴。但如果弯曲得更厉害,它就会保持弯曲状态。它发生了塑性变形。发生了什么?材料经历了永久、不可逆的变化。为了描述这一点,我们必须将总变形分为两部分:可恢复的弹性部分和永久的塑性部分。这正是乘法分解所捕捉到的洞见: F=FeFp\mathbf{F} = \mathbf{F}_e \mathbf{F}_pF=Fe​Fp​ 想象材料首先发生塑性变形,形成一个新的、无应力的“中间”形状,然后从这个新形状弹性拉伸到其最终的加载形态。这个优雅的思想是现代金属塑性理论的基石。它使工程师能够模拟复杂的工艺,如将一块平板钢材冲压成车门,或分析结构在车祸中的响应,在这些情况下,材料被推向远超其弹性极限的状态。

但故事并非止于弯曲。如果你反复来回弯曲回形针,它最终会折断。为什么?如果你在显微镜下观察金属,你会发现它并非一个完美、均匀的连续体。它充满了微观的空洞和缺陷。当材料发生塑性变形时,这些微小的空洞开始生长、拉伸并连接在一起。这个“损伤”过程逐渐从内部削弱材料。像 Gurson-Tvergaard-Needleman (GTN) 理论这样的模型为这个过程提供了数学描述。材料的屈服行为不再仅仅是应力的函数,还与“空隙体积分数”——一种衡量这种内部损伤的度量——有关。通过将此纳入有限应变框架,我们不仅可以预测材料将如何弯曲,还可以预测它何时以及如何断裂。这对于确保从飞机部件到桥梁等一切事物的安全性和可靠性具有极其重要的意义。

柔软与柔韧:橡胶、凝胶和组织

现在让我们把注意力从坚硬、强韧的金属转向软物质世界。想象一个橡胶气球、一块明胶或生物组织。你可以毫不费力地将它们拉伸和扭曲成非凡的形状,但试着将它们压缩成更小的体积——这几乎是不可能的。这些材料被归类为近乎不可压缩。

模拟这种行为提出了一个有趣的挑战。抵抗形状变化的刚度非常低,但抵抗体积变化的刚度(体积模量)却巨大。仅使用位移的标准计算方法可能导致数值病态,这个问题被称为“体积锁定”。有限应变框架通过所谓的混合公式提供了一个漂亮的解决方案。我们不再仅仅求解位移场,而是引入第二个独立的场:压力 ppp。系统的势能随后被写成位移和压力的函数。位移部分描述剪切和扭曲,而压力部分则强制执行近乎不可压缩的约束。这种方法优雅地避开了数值问题,并能够对从橡胶气象气球的充气到心脏跳动的力学等各种现象进行准确而稳定的模拟。有限应变理论与生物力学之间的这种联系是现代研究中一个充满活力且能拯救生命的领域。

时间的流动:聚合物与粘弹性

到目前为止,我们讨论的材料都是瞬时响应力的。但许多材料有记忆。它们根据你变形的速度表现出不同的行为。想想“傻瓜橡皮泥”:慢慢拉它,它像太妃糖一样伸展;猛地一拉,它像脆性固体一样断裂。这种时间依赖性行为被称为*粘弹性*。

为了模拟这类材料,我们必须描述应力如何随时间松弛,或者应变在恒定载荷下如何蠕变。这通常通过将材料想象成“麦克斯韦单元”的集合——一个弹簧(弹性部分)和一个阻尼器(粘性部分)串联而成。但在有限应变背景下,我们遇到了一个微妙而深刻的问题。如果一个物体在变形的同时也在旋转,我们如何定义“应力的变化率”?简单的时间导数是不够的,因为它被旋转所污染。我们需要一种方法来测量一个与材料元素一同旋转的观察者所看到的应力变化率。我们需要一个*客观应力率*。

连续介质力学提供了几种这样的量,如 Jaumann 率或上随钻率。通过使用这些客观应力率来构建我们的粘弹性演化定律,我们可以建立与物理学基本原理一致的模型。这些模型通常基于多重松弛机制的乘法分解,对于设计和加工聚合物、理解生物流体的流变学以及预测塑料部件在负载下的长期行为至关重要。

我们脚下的土地:地质力学

同样的基本有限应变原理也适用于更宏大的尺度,支配着我们脚下土地的行为。土壤、沙子和粘土是引人入胜的材料,其行为完全取决于围压。一把干沙像液体一样流动,但深层沙堆底部的沙子却坚硬如固体。

临界状态土力学利用在压力 ppp 和剪应力 qqq 空间中定义的塑性模型来捕捉这种行为。例如,著名的修正剑桥模型在这个空间中使用了一个优雅的椭圆形屈服面。随着土壤被压缩,它会硬化,屈服椭圆会变大。当被剪切时,它会发生塑性变形,直到达到一个“临界状态”,此时它可以在恒定应力和恒定体积下继续变形。

此外,这些材料的颗粒性质引入了另一层复杂性:各向异性。沙粒不是完美的球体,它们的堆积方式创造了一种“组构”,赋予材料方向性属性。随着沙子变形,这种组构会演变——颗粒旋转和重新排列。先进的有限应变模型结合了*组构张量*的演化,这是一个追踪材料内部结构的客观量。这使得工程师能够建立更真实的模型来预测斜坡的稳定性、地震中地基的行为以及复杂的地质过程。

我们生活的织物:多物理场在行动

或许,有限应变理论力量最优雅的展示之一在于其连接不同物理领域的能力。考虑一块潮湿天气下的机织物。你可能会注意到它感觉更重或悬挂方式不同。这是一个真实的物理效应:纺织纤维吸收了空气中的水分子,导致它们膨胀。

这种现象是化学-力学耦合的完美例子。化学环境(湿度)的变化驱动了物理变形(膨胀)。我们可以使用乘法分解概念的一个优美扩展来对此建模。一根纱线的总拉伸 λ\lambdaλ 可以分解为承载应力的力学部分 λm\lambda_mλm​ 和由局部水分含量决定的膨胀部分 λsw\lambda_{\mathrm{sw}}λsw​:λ=λswλm\lambda = \lambda_{\mathrm{sw}} \lambda_mλ=λsw​λm​。如果一块织物的边缘被固定,然后暴露在湿度梯度中,这种差异性膨胀将引起复杂的内部应力状态。一些纱线将被拉伸,另一些将被压缩,所有这一切都没有施加任何外力。有限应变分析使我们能够计算这些内力并预测织物最终的翘曲形状。同样的原理也适用于热膨胀、凝胶的溶胀,甚至生物组织的生长。

从地壳中的巨大压力到衬衫中纤维的精细相互作用,有限应变的语言提供了一个稳健而统一的框架。它证明了对几何和运动的深入细致的思考如何能够开启对物质世界所有丰富多样复杂性的深刻理解。