try ai
科普
编辑
分享
反馈
  • 有限应变力学

有限应变力学

SciencePedia玻尔百科
核心要点
  • 变形梯度张量(F\boldsymbol{F}F)是有限应变力学中的基本数学工具,它精确描述了局部的拉伸、剪切和旋转。
  • 物质标架无关性原理(客观性)要求本构律依赖于客观张量,例如右柯西-格林张量(C\boldsymbol{C}C),以确保物理定律与观察者无关。
  • 变形梯度的乘法分解(F=FeFp\boldsymbol{F} = \boldsymbol{F}_e \boldsymbol{F}_pF=Fe​Fp​)为模拟经历弹性(可恢复)和塑性(永久)变形的材料提供了一个强大的框架。
  • 有限应变力学对于精确模拟各种现象至关重要,从橡胶和生物组织的非线性行为到先进智能材料的设计。

引言

当橡皮筋被拉伸到极限或金属被永久弯曲时,入门力学中的简单规则便不再适用。这些通常是不可逆的大变形属于有限应变力学的范畴,这是一种对现代工程和科学至关重要的、更通用、更强大的理论。在分析软材料、制造过程和生物系统时,线性小应变近似的不足造成了巨大的知识鸿沟。本文通过对该主题进行全面介绍来弥合这一鸿沟。文章首先建立基本的“原理与机制”,探讨大变形的数学语言,包括变形梯度、客观应变和应力张量,以及能量和功的统一概念。随后,文章在“应用与跨学科联系”中展示了该框架的强大功能,说明了它如何解释惊人的物理效应,并如何能够模拟从生物组织到先进智能材料的各种事物。

原理与机制

想象一下,拉伸一根即将断裂的橡皮筋,揉一个面团,或者在慢动作视频中观看汽车保险杠的褶皱变形。这些都不是你在入门课程中可能学到的、桥梁在交通负载下那种温和、几乎看不见的变形。这些是巨大的、剧烈的,有时甚至是永久的形状变化。为了建立关于这些现象的科学,我们不能只做微小近似。我们需要一套新的工具,一种新的思维方式,它要像变形本身一样灵活而强大。这就是有限应变力学的世界。

运动的映射:变形梯度

我们的首要任务是描述物体如何移动和变形。让我们考虑一个处于其原始、未变形状态的物体——我们称之为​​参考构型​​。我们可以用其位置矢量 X\boldsymbol{X}X 来标记该物体中的每一个粒子。现在,我们施加一些力,物体移动并扭曲成某个稍后时刻的新形状。这个新形状就是​​当前构型​​。原来位于 X\boldsymbol{X}X 的粒子现在处于一个新的位置 x\boldsymbol{x}x。

物理学是关于寻找支配变化的定律,因此我们需要一个函数,一个“映射”,将“之前”与“之后”联系起来:x=χ(X)\boldsymbol{x} = \boldsymbol{\chi}(\boldsymbol{X})x=χ(X)。这个映射告诉我们每个粒子去了哪里。

但我们通常更关心局部发生的情况,即一个点附近的微小邻域。如果我们取一个从点 X\boldsymbol{X}X 开始的无穷小矢量 dXd\boldsymbol{X}dX,它在变形后的物体中会变成什么?它会变成一个新的矢量 dxd\boldsymbol{x}dx。由于我们观察的是一个无穷小的区域,这个映射看起来是线性的,就像一个曲面在足够放大后看起来是平的一样。这个局部线性映射是我们这场戏中无可争议的主角:​​变形梯度​​ F\boldsymbol{F}F。

它被定义为位置映射的梯度:F=∂x∂X\boldsymbol{F} = \frac{\partial \boldsymbol{x}}{\partial \boldsymbol{X}}F=∂X∂x​。这个偏导数矩阵告诉我们关于局部变形的一切——拉伸、剪切和旋转。它将无穷小矢量从参考物体转换到当前物体:dx=FdXd\boldsymbol{x} = \boldsymbol{F} d\boldsymbol{X}dx=FdX。

你可能会想考虑位移 u=x−X\boldsymbol{u} = \boldsymbol{x} - \boldsymbol{X}u=x−X,它是指一个粒子从原来的位置指向现在位置的矢量。它的梯度 ∇u\nabla \boldsymbol{u}∇u 当然与变形有关。实际上,稍作代数运算可知 F=I+∇u\boldsymbol{F} = \boldsymbol{I} + \nabla \boldsymbol{u}F=I+∇u,其中 I\boldsymbol{I}I 是单位矩阵。对于非常小的变形,其中 ∇u\nabla \boldsymbol{u}∇u 的分量非常小,F\boldsymbol{F}F 几乎与 I\boldsymbol{I}I 相同,∇u\nabla \boldsymbol{u}∇u 捕捉了大部分的物理现象。但在我们这个大变形的世界里,情况并非如此。变形梯度 F\boldsymbol{F}F 是更基本的量,因为它描述了映射本身的最终状态,而不仅仅是达到该状态的微小变化。

五花八门的应变:我们如何测量变形?

知道微小矢量如何变换是一回事,但我们如何得到一个数字或一个张量来表示“它应变了多少”?让我们看看我们的小矢量 dXd\boldsymbol{X}dX 的长度平方的变化。最初,它是 ∣dX∣2=dX⋅dX|d\boldsymbol{X}|^2 = d\boldsymbol{X} \cdot d\boldsymbol{X}∣dX∣2=dX⋅dX。变形后,它的新长度平方是 ∣dx∣2=dx⋅dx|d\boldsymbol{x}|^2 = d\boldsymbol{x} \cdot d\boldsymbol{x}∣dx∣2=dx⋅dx。使用我们的新工具,我们可以写出:

∣dx∣2=(FdX)⋅(FdX)=dXTFTFdX|d\boldsymbol{x}|^2 = (\boldsymbol{F} d\boldsymbol{X}) \cdot (\boldsymbol{F} d\boldsymbol{X}) = d\boldsymbol{X}^T \boldsymbol{F}^T \boldsymbol{F} d\boldsymbol{X}∣dx∣2=(FdX)⋅(FdX)=dXTFTFdX

看看中间的那个组合:FTF\boldsymbol{F}^T \boldsymbol{F}FTF。它非常重要,以至于有自己的名字:​​右柯西-格林变形张量​​,C=FTF\boldsymbol{C} = \boldsymbol{F}^T \boldsymbol{F}C=FTF。它是一个对称张量,巧妙地包含了关于一个点邻域内长度和角度如何变化的所有信息。

长度平方的变化量则为 ∣dx∣2−∣dX∣2=dXT(C−I)dX|d\boldsymbol{x}|^2 - |d\boldsymbol{X}|^2 = d\boldsymbol{X}^T (\boldsymbol{C} - \boldsymbol{I}) d\boldsymbol{X}∣dx∣2−∣dX∣2=dXT(C−I)dX。这自然引出了​​格林-拉格朗日应变张量​​的定义:

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

这个张量非常优美,因为它当且仅当局部没有任何变化时才为零。而且因为它是在参考构型(即原始、未变形的物体)上定义的,所以对于许多理论来说,它是一个自然的选择。

那么体积变化呢?这也隐藏在 F\boldsymbol{F}F 之中。当前构型中的一个小体积元 dVdVdV 与其原始体积 dV0dV_0dV0​ 之比由 F\boldsymbol{F}F 的行列式给出,即​​雅可比行列式​​ J=det⁡(F)J = \det(\boldsymbol{F})J=det(F)。对于任何真实材料,我们必须有 J>0J > 0J>0,因为你不能让物质消失或占据负体积。

现在,一个好奇的物理学家可能会问:为什么要基于初始构型来定义应变?我们不能基于最终构型来定义吗?是的,我们可以!这引出了另一种度量,即​​欧拉-阿尔曼西应变张量​​ e=12(I−B−1)\boldsymbol{e} = \frac{1}{2}(\boldsymbol{I} - \boldsymbol{B}^{-1})e=21​(I−B−1),其中 B=FFT\boldsymbol{B} = \boldsymbol{F}\boldsymbol{F}^TB=FFT 是​​左柯西-格林张量​​。对于小应变,E\boldsymbol{E}E 和 e\boldsymbol{e}e 几乎相同。但对于大应变,它们讲述的是不同的故事。

这个选择重要吗?哦,它至关重要!让我们做一个思想实验。如果我们试图将一种材料压碎至虚无,会发生什么?这对应于拉伸比 λ\lambdaλ 趋近于零。让我们看看我们的应变度量如何反应:

  • 格林-拉格朗日应变 EEE 趋近于一个有限值 −1/2-1/2−1/2。
  • 欧拉-阿尔曼西应变 eee 骤降至 −∞-\infty−∞,其行为类似于 −1/(2λ2)-1/(2\lambda^2)−1/(2λ2)。
  • 另一种度量,​​亨基(或对数)应变​​ HHH 也趋向于 −∞-\infty−∞,但更温和,像 ln⁡(λ)\ln(\lambda)ln(λ)。

这不仅仅是一个数学游戏;它具有深刻的物理后果。如果我们基于格林-拉格朗日应变 E\boldsymbol{E}E 建立材料能量理论,那么当体积坍缩到零时,能量将趋于一个有限值。在我们的理论中,材料将“感觉”不到即将到来的灾难!然而,一个基于 e\boldsymbol{e}e 或 H\boldsymbol{H}H 的理论会看到能量飙升至无穷大,形成一个强大的屏障,反映了你不能将物质压碎至消失的物理现实。数学定义的选择,就是对你希望描述的物理现象的选择。

客观性:物理学对每个人都是相同的

物理学中最深刻的原则之一是,自然定律不应依赖于观察者。如果一种材料只是静静地放在那里,它的内部应力状态不应因为你开始旋转着观察它而突然改变。这就是​​物质标架无关性原理​​,或称​​客观性​​。

这对我们的理论有巨大影响。如果我们用一个旋转矩阵 Q\boldsymbol{Q}Q 来旋转变形后的物体,新的变形梯度是 F∗=QF\boldsymbol{F}^* = \boldsymbol{Q}\boldsymbol{F}F∗=QF。它改变了!所以,F\boldsymbol{F}F 是非客观的。这意味着我们不能将我们的材料定律,比如能量函数,仅仅建立在 F\boldsymbol{F}F 之上。

但请看这个魔法。右柯西-格林张量 C\boldsymbol{C}C 会发生什么变化? C∗=(F∗)TF∗=(QF)T(QF)=FTQTQF=FTIF=C\boldsymbol{C}^* = (\boldsymbol{F}^*)^T \boldsymbol{F}^* = (\boldsymbol{Q}\boldsymbol{F})^T (\boldsymbol{Q}\boldsymbol{F}) = \boldsymbol{F}^T \boldsymbol{Q}^T \boldsymbol{Q} \boldsymbol{F} = \boldsymbol{F}^T \boldsymbol{I} \boldsymbol{F} = \boldsymbol{C}C∗=(F∗)TF∗=(QF)T(QF)=FTQTQF=FTIF=C 它保持不变!C\boldsymbol{C}C 是一个客观张量。这就是为什么弹性材料的材料定律几乎总是用 C\boldsymbol{C}C 或其​​主不变量​​(I1=tr CI_1 = \mathrm{tr}\,\boldsymbol{C}I1​=trC、I2I_2I2​、I3=det⁡CI_3 = \det \boldsymbol{C}I3​=detC)来表示的原因,它们也是客观标量。数学尊重物理。

应力、功率和功共轭

现在我们能描述变形了,让我们来谈谈材料内部的力——应力。应力是单位面积上的力。但是哪个面积?是原始、未变形物体中的面积,还是当前、变形后的面积?这种模糊性导致了一系列的应力张量。

  • ​​柯西应力​​ σ\boldsymbol{\sigma}σ 是“真实”应力。它是作用在当前物体中某个面积上的、位于当前物体中的力。它是对称的,也是工程师会尝试测量的量。
  • ​​第一皮奥拉-基尔霍夫应力​​ P\boldsymbol{P}P 是一个混合量。它将当前物体中的力与参考物体中的面积联系起来。它很有用,因为我们是在固定的参考物体上积分,但它通常是不对称的。
  • ​​第二皮奥拉-基尔霍夫应力​​ S\boldsymbol{S}S 是一个更抽象的度量,但它非常优美。它是一个完全拉回到参考构型的对称应力张量。

我们为什么需要这么多应力?答案在于能量和功率。单位体积的应力做功速率是一个所有观察者都必须认同的物理量。单位当前体积的功率由 σ:d\boldsymbol{\sigma}:\boldsymbol{d}σ:d 给出,其中 d\boldsymbol{d}d 是速度梯度的对称部分。如果我们将这同一个功率用参考体积来表示,我们会发现它呈现出不同的形式:

单位参考体积的功率=P:F˙=12S:C˙\text{单位参考体积的功率} = \boldsymbol{P}:\dot{\boldsymbol{F}} = \frac{1}{2}\boldsymbol{S}:\dot{\boldsymbol{C}}单位参考体积的功率=P:F˙=21​S:C˙

这揭示了一种深刻的对称性。应力-应变对 (P,F)(\boldsymbol{P}, \boldsymbol{F})(P,F) 和 (S,E)(\boldsymbol{S}, \boldsymbol{E})(S,E)(或等效地,(S,C)(\boldsymbol{S}, \boldsymbol{C})(S,C))是​​功共轭​​的。它们是描述能量的“正确”应力和应变对。功共轭原理是允许我们在不同运动学和动力学描述之间进行转换的“罗塞塔石碑”,最重要的是,它使我们能够从能量推导出应力。

材料行为:游戏规则

有了这个框架,我们现在可以描述材料的实际行为了。

超弹性:完美的弹簧

​​超弹性​​材料,就像理想的橡皮筋,是一种将其上所做的所有功都储存为可恢复势能的材料。存在一个​​储存能函数​​ Ψ\PsiΨ,它只依赖于变形。由于客观性原理,我们必须将其写成一个客观张量(如 C\boldsymbol{C}C)的函数:Ψ=Ψ(C)\Psi = \Psi(\boldsymbol{C})Ψ=Ψ(C)。对于各向同性材料(在所有方向上看起来都一样),这可以简化为不变量的函数:Ψ(I1,I2,I3)\Psi(I_1, I_2, I_3)Ψ(I1​,I2​,I3​)。

现在,功共轭的魔力显现了:因为 S\boldsymbol{S}S 与 E\boldsymbol{E}E 是功共轭的,所以应力就是能量的导数! S=∂Ψ∂E=2∂Ψ∂C\boldsymbol{S} = \frac{\partial \Psi}{\partial \boldsymbol{E}} = 2\frac{\partial \Psi}{\partial \boldsymbol{C}}S=∂E∂Ψ​=2∂C∂Ψ​ 这是力学的一个基石。给定一个描述材料的能量函数,我们可以推导出任何变形下的应力响应。思想实验,比如考虑一个应力功率恒为零的材料,甚至可以揭示能量函数本身的基本属性,例如其齐次性。

通常,将变形分解为改变体积的部分(体积变形)和只改变形状的部分(等容变形)是很有用的。对于像橡胶这样易于扭曲但难以压缩的材料,这一点至关重要。数学提供了优雅的方法来做到这一点,例如通过定义一个对体积变化不敏感的修正或“等容”变形张量。

塑性:不归点

当你弯曲一个回形针时会发生什么?它会保持弯曲。这种变形是永久性的。这就是​​塑性​​。这是一个复杂得多的世界,因为材料的内部状态正在以不可逆的方式改变。

处理这个问题的关键思想是由 E. H. Lee 提出的变形梯度的​​乘法分解​​。该理论假设总变形 F\boldsymbol{F}F 可以被看作一个两步过程:

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

首先,材料经历一个不可逆的塑性变形 Fp\boldsymbol{F}_pFp​。这将其从参考构型带到一个概念性的​​中间构型​​。这个构型是无应力的,但其内部一团糟——它包含了所有永久性的重排,比如晶体中的位错。它通常是“不协调的”,意味着你无法切出一块让它平放。然后,从这个新的重排状态,材料通过 Fe\boldsymbol{F}_eFe​ 进行弹性变形,达到最终的当前构型。所有储存的弹性能和应力都来自 Fe\boldsymbol{F}_eFe​。

当我们考察变化率时,我们发现塑性流动由​​塑性速度梯度​​ Lp=F˙pFp−1\boldsymbol{L}_p = \dot{\boldsymbol{F}}_p \boldsymbol{F}_p^{-1}Lp​=F˙p​Fp−1​ 控制。它的对称部分 Dp\boldsymbol{D}_pDp​ 代表塑性应变率,能量正是在这里损失或​​耗散​​掉的,通常以热的形式。热力学基本定律要求这种耗散必须为非负。这就是​​Drucker稳定性公设​​的核心:塑性变形总是消耗能量。这个用功共轭量以客观方式陈述的原理,无论人们对响应的弹性部分做出何种具体的建模选择,都同样成立。

从点之间映射的简单想法出发,我们建立了一个丰富而强大的结构,能够描述物质在大变形下的复杂舞蹈,从完美的弹性回弹到不可逆的塑性流动,所有这些都被能量和客观性这两个优美而统一的原则联系在一起。

应用与跨学科联系

想象一下拉伸一根橡皮筋。它当然会变长。但它还会发生别的变化吗?它会变热吗?它的宽度会以一种简单的方式变化吗?在微小、温和的推拉世界——即线性弹性的领域——中,答案是异常简单的。但是当你真正地推动物体时,当你拉伸、扭转和剪切它们直到它们“叫苦不迭”时,一种全新的、远为丰富的物理学便浮现出来。这就是有限应变力学的世界。它不仅仅是一种数学上的奇观;它是描述我们周围世界真实行为的基本语言,从冰川的流动到心脏的跳动。在上一章建立了理论框架之后,现在让我们走出去,看看这个框架能让我们理解和构建哪些宏伟的结构。

大变形的几何学:惊人的后果

让我们从一个简单的思想实验开始。想象一下剪切一副扑克牌。顶部相对于底部滑动。在小应变的简单世界里,故事到此结束。但在现实世界中,还会发生别的事情。如果你能以这种方式剪切一块橡胶,你会发现它还需要在垂直于剪切的方向上拉伸!这种被称为 Poynting 效应的显著现象,是大变形几何学的直接结果,。忽略“小量平方”项的小应变理论对此效应是“视而不见”的,因为它产生的正应力与剪切量的平方 γ2\gamma^2γ2 成正比。这是一个绝佳的例子,说明了自然界如何将秘密隐藏在我们常常倾向于忽略的高阶项中。

为什么要在我们之前遇到的所有复杂张量上费心呢?为什么不坚持使用简单的线性图像?答案是,简单的图像常常是错误的,而且错得离谱。对于简单的剪切变形,使用线性化应变而非正确的格林-拉格朗日应变所产生的误差会随着剪切量的平方而增长。对于小剪切,误差可能微不足道,但对于大剪切——这在橡胶或生物组织等软材料中很常见——线性理论是完全不够用的。有限应变理论是进入大变形世界的入场券。

此外,当一个物体被严重变形时,即使是“它应变了多少?”这个问题也变得微妙起来。我们应该相对于物体原始、未变形的形状来测量应变吗?还是应该相对于其新的、变形后的形状来测量?这两种观点分别产生了不同的数学度量,如拉格朗日应变和欧拉应变。两者都并非比对方更“正确”;它们只是用于不同工作的不同工具,一个更适合跟踪材料历史,另一个更适合描述当前的应力状态。

物质世界:从原子到构造

材料抵抗变形的方式——即其刚度——并非某种神奇的属性。它源于其构成原子之间无数次的相互作用。想象两个由弹簧连接的原子。这太简单了。真实的原子键更像一个在压缩时变得更硬、在拉伸过远时会断裂的弹簧。力-距离曲线不是一条直线。有限应变理论使我们能够将这种原子间势的精确非线性形状与晶体的宏观弹性常数联系起来,包括那些只在大变形期间才起作用的高阶常数。这在固态物理和连续介质力学之间提供了一个深刻的联系,使我们能够从第一性原理预测材料的行为。

但是当变形是永久性的时候会发生什么?当你弯曲一个回形针时,它不会弹回。它经历了塑性变形。为了描述这一点,我们需要将变形分解为两部分:一个可恢复的“弹性”部分和一个永久的“塑性”部分。变形梯度的乘法分解 F=FeFp\boldsymbol{F} = \boldsymbol{F}_e \boldsymbol{F}_pF=Fe​Fp​ 是实现这种分离的优美数学工具。它构成了现代塑性理论的基石,使工程师能够以惊人的准确性模拟复杂的制造过程,如锻造、轧制和挤压。对于被推到极限的材料,比如在灼热温度下发光的喷气发动机涡轮中的合金,该框架被扩展以包含时间和温度的影响,这被称为粘塑性模型。当材料开始失效时,微小的空洞和微裂纹会累积,使材料软化。这种“损伤”也可以被编织到我们的有限应变框架中,使我们能够预测从混凝土桥梁到橡胶轮胎等各种事物的断裂起始。

连接世界的桥梁:多尺度建模与计算

如果没有能力求解这些复杂理论的方程,它们将几乎没有实际用途。这就是计算能力发挥作用的地方。有限元法(FEM)是现代工程的“主力军”,能够预测复杂结构对载荷的响应。在处理有限应变的非线性方程时,计算机实质上是通过一系列小的修正步骤“悄悄地”逼近解,这个过程被称为 Newton-Raphson 法。在每一步中,它不仅需要知道当前状态离平衡态有多远(“残差”),还需要知道如何调整变形以更接近平衡(“切线刚度”)。这两个量都直接源于我们讨论过的有限应变本构模型。

但如果材料本身是复杂的,是像碳纤维复合材料或骨骼一样由不同组分混合而成,那该怎么办?对每一根纤维或晶体进行建模是不可能的。取而代之,我们使用一个强大的思想,称为均匀化。我们非常详细地分析材料的一个小的“代表性体积元”(RVE)。挑战在于确保这个小单元的平均响应能准确反映大块材料的行为。Hill-Mandel 条件是保证这种联系的关键原则。它是一个能量一致性的陈述:作用于宏观材料上的功必须等于其微观组分内部所做功的平均值。这使我们能够建立对材料的“自下而上”的理解,通过设计其微观结构来实现所需的宏观性能。

生命世界及其他:生物力学、地质力学和智能材料

也许有限应变力学在生命研究中比在任何其他地方都更为重要。软生物组织——皮肤、肌肉、软骨——通常会经历对传统工程材料而言是灾难性的变形。要理解它们的功能,我们必须拥抱大应变的语言。例如,通过在软体蠕虫表面撒上斑点,并在其移动时用相机跟踪这些图案,生物学家可以使用数字图像相关(DIC)技术来获得其体壁上的全场应变图。用正确的有限应变运动学工具分析这些数据,揭示了该动物如何通过其静水骨架产生运动。有限应变度量是“客观的”——即对蠕虫扭动时的平凡刚体旋转不敏感——这一事实对于提取有意义的信息至关重要。

许多这些组织,如软骨,也是充满流体的多孔固体。土壤和岩石也是如此。这些材料的力学行为由变形的固体基质与流经其孔隙的流体压力之间微妙的相互作用所支配。有限应变多孔介质力学为模拟这种耦合提供了统一的框架,解释了从我们关节的润滑到山体滑坡的力学等一切现象。

有限应变力学的旅程并不仅仅止于描述世界本来的样子;它还延伸到创造一个可能的世界。想象一下,将一个微观的层状结构——一个一维光子晶体——嵌入到一个柔软、可拉伸的弹性体中。该结构反射的颜色取决于其层间距,这由布拉格定律决定。当你拉伸弹性体时,有限变形会以精确、可预测的方式改变这个间距,从而改变其颜色。这就是力致变色材料背后的原理,这种材料会响应机械力而改变颜色。有限应变力学与光学的这种美妙结合,为创造智能传感器、自适应伪装以及能够通过视觉传达其力学状态的软体机器人系统打开了大门。从纯粹的变形几何学出发,我们最终实现了功能、响应,甚至一种新的美。