try ai
科普
编辑
分享
反馈
  • 同旋率

同旋率

SciencePedia玻尔百科
关键要点
  • 材料坐标系无关性原理要求使用客观应力率,因为标准导数无法区分材料变形与刚体转动。
  • 同旋率是一种客观度量,通过在一个随局部材料区域一同旋转的参考系中观察应力变化来构建。
  • 同旋率的选择(例如,Jaumann率与对数率)对有限应变问题中应力演化、能量守恒和稳定性的预测有至关重要的影响。
  • 这些率是计算力学(FEM)和材料科学的基础,用于构建大变形下材料的物理精确且数值稳健的模拟。

引言

描述变形材料内部的力与应力如何变化,是物理学和工程学的基石。然而,当一个物体经历既包含拉伸又包含旋转的大变形时,一个根本性问题便随之出现:我们对于变化的直观度量可能会产生误导。物理学家常用的计时工具——标准时间导数,无法区分由真实材料变形引起的应力变化和由材料自身旋转引起的表观变化。这会导致模型违反物理学的一条神圣准则:材料坐标系无关性原理,该原理指出物理规律必须独立于观察者。

本文旨在通过引入同旋率这一概念来弥补这一关键知识空白。同旋率是一套精密的工具,旨在提供一种客观的、或称观察者无关的应力变化度量。通过两个章节,您将对这一重要主题获得深刻理解。“原理与机制”一章将剖析标准导数为何不适用,并揭示构建同旋率背后的精妙逻辑。随后的“应用与跨学科联系”一章将阐明这些概念的深远而实际的影响,展示率的选择如何影响从有限元法中的计算机模拟精度到我们对材料的热力学理解等方方面面。这段旅程始于审视那些要求我们以更复杂的视角看待变化率的核心原理。

原理与机制

想象一下,你正试图描述一块正在被拉伸和扭转的太妃糖的状态。它在某些地方伸长,在另一些地方压缩,并且整个过程都在旋转。你将如何描述其内部的力?更重要的是,你将如何描述这些力是如何变化的?事实证明,你的答案完全取决于你的视角。如果你在太妃糖旋转时静止站立,你会看到力的变化,这仅仅是因为太妃糖的不同部分在你面前旋转而过。但如果你是一只骑在太妃糖上的小瓢虫,随之一起旋转,你对力变化的描述将会截然不同。你只会注意到由于太妃糖实际拉伸或压缩而引起的变化。

这个简单的思想实验触及了物理学中一个深刻而优美的原理:​​材料坐标系无关性原理​​(Principle of Material Frame Indifference),或称​​客观性​​(objectivity)。这个听起来很复杂的名称背后是一个简单的思想:支配一种材料的物理规律不能依赖于观察者是谁,或者观察者如何运动。无论是地球上的工程师,还是在旋转空间站中的宇航员,描述钢材在负载下行为的方程必须是相同的。这是关于物理规律对称性的一个基本论断。

作为物理学家和工程师,我们的任务是写下本构律——即关联材料变形与其所受应力的规则。挑战在于,书写这些定律的方式必须尊重客观性原理。

观察者的困境:视角问题

让我们说得更精确一些。我们使用一个称为​​Cauchy应力张量​​的数学对象来度量材料的内力,我们用符号 σ\boldsymbol{\sigma}σ 来表示它。你可以把它想象成一台机器,它能告诉你穿过材料的任何假想平面上的面力(单位面积上的力)。现在,如果我们想建立一个描述材料行为的定律——比如一种粘弹性流体或一种正在锻造的金属——我们需要将这种应力的变化率与材料的形变率联系起来。

测量变化率最显而易见的方法就是用秒表——我们直接取标准的物质时间导数,σ˙\dot{\boldsymbol{\sigma}}σ˙。所以,一个直观的本构律猜想可能是这样的:σ˙=形变率的某个函数\dot{\boldsymbol{\sigma}} = \text{形变率的某个函数}σ˙=形变率的某个函数。让我们看看为什么这个简单的想法会彻底失败。

考虑一块材料,比如一块已经承受某种应力但现在正在进行纯刚体转动的果冻,其形状和大小没有变化。从物理角度看,既然材料没有变形,其内在的应力状态就不应该改变,不应产生新的应力。然而,一个在实验室里静止的观察者看到这块果冻在旋转。从他们固定的角度来看,应力张量的分量是在变化的。例如,一个原本沿x轴的拉伸现在指向了不同的方向。所以,对于这个静止的观察者,σ˙\dot{\boldsymbol{\sigma}}σ˙ 不为零!我们这个幼稚的本构律会错误地预测纯旋转会产生应力,这个结论与物理现实背道而驰。

这背后的数学原理非常清晰。如果一个“带星号”的观察者相对于一个“不带星号”的观察者以自旋张量 Ω\boldsymbol{\Omega}Ω 旋转,他们测量的应力 σ∗\boldsymbol{\sigma}^{\ast}σ∗ 通过一个简单的旋转与原始应力相关联:σ∗=QσQT\boldsymbol{\sigma}^{\ast} = \boldsymbol{Q}\boldsymbol{\sigma}\boldsymbol{Q}^{\mathsf{T}}σ∗=QσQT,其中 Q\boldsymbol{Q}Q 是旋转张量。但是当我们看时间导数时,情况就变得更复杂了。一点微积分知识表明其关系是:

σ˙∗=Qσ˙QT+Ωσ∗−σ∗Ω\dot{\boldsymbol{\sigma}}^{\ast} = \boldsymbol{Q}\dot{\boldsymbol{\sigma}}\boldsymbol{Q}^{\mathsf{T}} + \boldsymbol{\Omega}\boldsymbol{\sigma}^{\ast} - \boldsymbol{\sigma}^{\ast}\boldsymbol{\Omega}σ˙∗=Qσ˙QT+Ωσ∗−σ∗Ω

一个客观量应该能简单地变换为 Rate∗=Q(Rate)QT\text{Rate}^{\ast} = \boldsymbol{Q}(\text{Rate})\boldsymbol{Q}^{\mathsf{T}}Rate∗=Q(Rate)QT。物质导数 σ˙\dot{\boldsymbol{\sigma}}σ˙ 未能通过这个测试,就是因为多了 Ωσ∗−σ∗Ω\boldsymbol{\Omega}\boldsymbol{\sigma}^{\ast} - \boldsymbol{\sigma}^{\ast}\boldsymbol{\Omega}Ωσ∗−σ∗Ω 这两项。这个被称为对易子的项,代表了我们的测量被观察者自身自旋所带来的“污染”。我们的秒表坏了,因为它无法区分由于实际材料变形引起的应力变化和由于旋转视角引起的表观变化。

与材料同行:同旋思想

我们如何修理这个“坏掉的秒表”?解决方法既优雅又直观:如果从静止坐标系观察得到错误的答案,那么就让我们从一个随材料一起运动的坐标系来观察。具体来说,我们将从一个与材料局部区域一同旋转的微小观察者的角度来测量应力的变化率。这就是​​同旋率​​(corotational rate)背后的核心思想。

其几何意义非常强大。对于我们那块进行纯刚体转动的果冻,一个与它一同旋转的观察者会看到应力张量是完全恒定的。他们计算出的变化率为零,这正是物理上正确的答案,因为没有发生变形。

在数学上,我们通过取“坏的”物质导数 σ˙\dot{\boldsymbol{\sigma}}σ˙ 并减去那些导致问题的项来构造这样一个率,我们称之为 σ∘\overset{\circ}{\boldsymbol{\sigma}}σ∘。一个通用的同旋率定义为:

σ∘=σ˙+σω−ωσ\overset{\circ}{\boldsymbol{\sigma}} = \dot{\boldsymbol{\sigma}} + \boldsymbol{\sigma}\boldsymbol{\omega} - \boldsymbol{\omega}\boldsymbol{\sigma}σ∘=σ˙+σω−ωσ

这里,ω\boldsymbol{\omega}ω 是一个反对称张量,代表了我们决定“乘坐”的局部材料区域的自旋或角速度。这个新构造的率 σ∘\overset{\circ}{\boldsymbol{\sigma}}σ∘,根据其构造方式,就是客观的。其非客观部分奇迹般地抵消了。

有了这个工具,我们现在可以写出正确的本构律了。例如,一个描述类弹性材料的简单定律,称为次弹性定律(hypoelastic law),将客观应力率与形变率 D\boldsymbol{D}D 联系起来:

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

这个方程现在是站得住脚的了。它将一个客观量(同旋应力率)与另一个客观量(形变率)关联起来。它正确地预测了纯刚体转动(此时 D=0\boldsymbol{D}=\boldsymbol{0}D=0)的应力变化为零。我们成功地建立了一种独立于观察者的材料行为描述。

自旋的多种选择:该乘坐哪座旋转木马?

我们解决了一个问题,但一个新的、更微妙的问题出现了。我们所说的“材料的自旋”究竟指的是什么?对于像旋转陀螺这样的刚体,答案很简单——每个部分都有相同的角速度。但对于像我们那块太妃糖一样的变形体,情况就没那么简单了。不同部分可以相互拉伸和旋转。对于一个变形介质的“那个”自旋,没有一个上帝赋予的、唯一的定义。

这种模糊性意味着我们可以自由选择用于定义同旋率的自旋 ω\boldsymbol{\omega}ω。这导致了一整套不同的客观率的“动物园”的发展,每一种都对应于一种不同的物理或数学上的选择,即与材料一同旋转意味着什么。

让我们来认识一下这个动物园里几位最著名的成员:

  1. ​​Jaumann率:​​ 这也许是最直接的选择。它将自旋 ω\boldsymbol{\omega}ω 定义为空间速度梯度 L=∇v\boldsymbol{L}=\nabla\boldsymbol{v}L=∇v 的反对称部分。我们称之为​​自旋张量,W\boldsymbol{W}W​​。在物理上,它代表了一个材料点无限小邻域的平均瞬时角速度。它是一个纯粹的局部自旋度量。

  2. ​​Green-Naghdi率:​​ 这个率对旋转采取了一种更综合、更“整体”的观点。得益于力学中一个优美的定理——​​极分解​​,任何变形,由变形梯度张量 F\boldsymbol{F}F 描述,都可以被唯一地分解为一个纯拉伸,然后是一个纯刚体旋转。这个旋转由一个张量 R\boldsymbol{R}R 捕捉。Green-Naghdi率的定义是通过选择自旋为该物质旋转张量的变化率,即 ω=R˙RT\boldsymbol{\omega} = \dot{\boldsymbol{R}}\boldsymbol{R}^{\mathsf{T}}ω=R˙RT。它在一个与材料主变形轴一同旋转的坐标系中测量变化率。

  3. ​​Truesdell率:​​ 这是另一个重要的率,它源于一个与不同应力度量(第二Piola-Kirchhoff应力)相关的不同视角。它的形式更复杂,事实上,它不是我们一直在讨论的那种简单的同旋率形式,但它仍然是客观的并被广泛使用。

为了保证客观性,我们选择的任何自旋张量 ω\boldsymbol{\omega}ω 都必须满足一个特定的变换规则:在观察者变换下,它必须变换为 ω∗=QωQT+Ω\boldsymbol{\omega}^{\ast} = \boldsymbol{Q}\boldsymbol{\omega}\boldsymbol{Q}^{\mathsf{T}} + \boldsymbol{\Omega}ω∗=QωQT+Ω,其中 Ω\boldsymbol{\Omega}Ω 是观察者的自旋。Jaumann自旋 W\boldsymbol{W}W 和Green-Naghdi自旋 R˙RT\dot{\boldsymbol{R}}\boldsymbol{R}^{\mathsf{T}}R˙RT 都满足这个关键要求。

选择的后果:当不同视角至关重要时

这种对自旋的选择——选择乘坐哪座旋转木马——真的重要吗?答案是一个引人入胜的“视情况而定”。

对于绝大多数日常工程问题——设计变形量极小的建筑物、桥梁或机器零件——这完全不重要。在​​无穷小应变​​的世界里,旋转非常小,以至于客观率中的修正自旋项是高阶小量,可以安全地忽略。所有客观率基本上都等同于简单的物质时间导数,经典的线性弹性理论工作得很好。

但在​​有限应变​​的世界里——在涉及橡胶、塑料、金属成形或生物组织的问题中——这个选择至关重要,并具有深远的物理意义。

  • ​​轴与值的舞蹈:​​ Jaumann率和Green-Naghdi率之间有什么区别?数学揭示了一个美妙的秘密。这两个率之间的差是一个对易子,[σ,R˙RT−W][\boldsymbol{\sigma}, \dot{\boldsymbol{R}}\boldsymbol{R}^{\mathsf{T}}-\boldsymbol{W}][σ,R˙RT−W]。当我们在与应力主方向对齐的坐标系中写出它时,这个差值张量是纯非对角的。用通俗的话说,这意味着对于给定的变形,使用Jaumann率的模型和使用Green-Naghdi率的模型将预测​​主应力值(拉伸和压缩的大小)的变化率完全相同​​。它们的区别在于对​​主应力轴旋转速率​​的预测。客观率的选择决定了应力方向如何演变,而不是大小。

  • ​​虚假振荡:​​ 这种预测旋转的差异具有实际后果。考虑一块经历大的、连续的“简单剪切”变形的材料。一个使用Jaumann率的简单次弹性模型著名地预测,随着变形的增加,剪切应力会以一种奇怪的、非物理的方式振荡。而一个使用Green-Naghdi率的模型,它在运动学上与整体材料旋转联系更紧密,通常会避免这种虚假行为,并给出更合理的结果。

  • ​​通往真正弹性的道路:​​ 我们选择的最深远后果在于它与能量和热力学的联系。一种“真正的”弹性材料,我们称之为​​超弹性​​材料,应该像一个完美的弹簧一样储存变形能。你为使其变形所做的功,在卸载时应该能完全收回;在变形的闭合循环上所做的功必须为零。一个次弹性(率型)模型并不能自动保证这一点。事实上,基于Jaumann率的模型通常是不可积的——它可能预测在一个闭合的变形循环中能量被创造或毁灭,这是非物理的。这揭示了Jaumann率模型充其量只是一个耗散过程的模型,而不是纯弹性过程的模型。于是问题出现了:是否存在任何一种率的选择是可积的?答案是肯定的。研究表明,存在一种非常特殊的同旋率,称为​​对数率​​,它是根据对数应变定义的。一个使用这种特定率的次弹性定律是可积的,并且精确地对应于一个定义明确的超弹性模型,即Hencky模型。

这段从一个关于观察者的简单问题到热力学深层结构的旅程,展示了力学的非凡统一性和美感。对同旋率的需求不仅仅是一个数学上的细微差别;它是使我们的物理定律保持一致和有意义的基本要求,揭示了运动学、动力学和材料定义本身之间错综复杂的舞蹈。

应用与跨学科联系

在我们之前的讨论中,我们进入了运动学的微妙世界,发现在物体既变形又旋转时,我们日常简单的“变化率”概念变得不足。我们学会了构建新的数学工具——客观的,或称同旋的率——它们可以巧妙地将纯旋转从真实变形中分离出来。这无疑是一场愉悦的数学练习,但你可能想知道,“为什么要费这么大劲?” 答案是,这个概念不仅仅是一个数学上的奇物;它是一个基础支柱,支撑着我们在众多学科中对物理世界的理解。它是确保我们物理定律诚实可靠的守门人。

物理学中至高无上的法则是,一个完善的原理必须对每一个观察者都成立。无论你是静止站立还是在旋转木馬上旋转,自然法则都必须适用。这就是​​材料坐标系无关性原理​​,它是不容置疑的。如果你试图用简单的物质时间导数 σ˙\dot{\boldsymbol{\sigma}}σ˙ 来书写材料如何响应力的定律,你会发现你测量到的变化中有一部分只是由旋转引起的幻觉。一个纯粹的刚体转动,没有任何变形,竟然也会表现出应力率! 这在物理上是荒谬的。客观率正是我们的补救措施。它们的构造就是为了对刚性旋转“视而不见”,确保我们写下的定律是普适的,并且只依赖于材料真实的、物理的变形。

剪切的考验:正应力的奇特现象

让我们来探索一个看似简单的实验。想象一块厚厚的橡胶块,或者像玉米糖浆一样非常粘稠的流体,夹在两块板之间。现在,以恒定的速度侧向滑动顶板。这被称为简单剪切。我们的直觉告诉我们,为了保持顶板移动,我们必须施加一个平行于板的剪切应力。但垂直于板的力呢?材料是应该推开板还是将它们拉拢?

这正是我们选择的客观率变得极为重要的地方。如果我们用经典的​​Zaremba–Jaumann率​​为我们的材料建立一个简单的“次弹性”模型,我们会得到一个怪异的预测:随着剪切变形的增加,垂直于板的正应力会振荡,交替地推和拉。对于许多真实材料,这根本不是发生的情况。这是一个数学上的假象,是我们微积分工具选择所产生的非物理幽灵。

现在,让我们尝试一个不同的工具:​​对数同旋率​​。当在同一模型中使用这个更复杂的率时,它预测正应力将稳定增长并将板推开。这个预测与我们在许多非牛顿流体和聚合物中观察到的现象——即所谓的Weissenberg效应——更为一致。Weissenberg效应是指被搅拌的聚合物会沿着旋转的搅拌棒向上爬。数学“导数”的选择完全改变了物理预测,将一个振荡的幻影变成了一种真实、可测量的力。 显然,这不仅仅是一场符号游戏;这是一场寻求能讲出自然真相的数学语言的探索。

热力学约束:避免能量无中生有

宇宙对能量有严格的规定,其中最著名的一条是:你不能无中生有。一种理想的弹性材料,就像一根完美的弹簧,在变形时应该储存能量,并在释放时全部归还。如果你使其变形并通过一个闭合循环回到原始形状,所做的净功必须为零。任何其他结果都意味着你要么创造了能量,要么毁灭了能量,从而造出了一台永动机。

在这里,我们再次发现大多数简单的客观率都有一个不可告人的秘密。如果你使用Jaumann率或Green-Naghdi率来构建一个弹性定律,并让材料经历一个包含旋转的闭合弹性变形过程,模型可能会预测出能量的净增加或减少。 这是物理一致性的一个深刻失败。该模型缺乏“热力学良知”。

再一次,一个特定的选择前来拯救我们。​​对数率​​,当与“Hencky”对数应变配对时,是特殊的。这种组合在数学上保证了对于各向同性材料是​​可积的​​,意味着它可以从一个真实的弹性势能函数推导出来。 以这种方式建立的模型是一个​​超弹性​​模型,它在闭合的弹性循环上总能得到零净功。它正确地守恒了能量。这揭示了运动的几何学与热力学定律之间的深刻联系,表明了对运动学的谨慎选择如何能满足物理学最基本的原理之一。

数字宇宙:在计算机内部构建材料

在现代世界,从汽车到飞机再到生物医学植入物,许多工程设计都在计算机内完成。工程师们使用基于有限元方法(FEM)的强大软件,在切割任何一块金属之前,构建并测试虚拟原型。这些模拟依赖于本构模型来告诉计算机材料如何表现。在这里,客观率的抽象概念具有直接、实际且代价高昂的后果。

想象一个涉及大变形的金属成形过程的模拟。软件通过小的时间步长推进来求解运动方程。在每一步中,它都必须更新材料中的应力。如果程序员使用一个幼稚的算法来积分一个简单的次弹性定律,一个奇怪的假象可能会出现:​​人为硬化​​。 即使模型是针对“理想塑性”材料(这种材料不应随更多变形而变得更强),模拟也会显示应力随着剪切持续上升。这是由于在许多时间步长上不当累积旋转误差而引起的数值幻觉。模拟在撒谎,可能导致有缺陷的设计。

解决方法是使用更复杂的积分算法,这些算法正是建立在我们讨论的原理之上的。这些算法的工作方式是在每个时间步中将问题转换到​​同旋坐标系​​,执行仅涉及变形的计算部分,然后将结果旋转回来。这通常通过一个数学上精确的操作,称为​​指数映射​​来实现。 这正确地将变形与旋转分离开来,并驱除了人为硬化的幽灵。

此外,客观率的选择影响了模拟的数值效率。计算上稳健的率,如对数率,通常会得到一个“对称算法切线模量”,这是一个数学特性,使得复杂的方程组能够更快、更稳定地求解。一个糟糕的选择会导致非对称系统,这对数值求解器来说是个头痛的问题。 因此,正确的同旋率选择不仅能给出正确的物理答案,还能帮助计算机高效地找到那个答案。

深入晶格:当自旋变为现实

到目前为止,在不同的自旋——Jaumann自旋W\mathbf{W}W、Green-Naghdi自旋R˙RT\dot{\mathbf{R}}\mathbf{R}^TR˙RT等——之间的选择,似乎只是数学品味或便利性的问题。但当材料本身具有内在的方向性时,会发生什么呢?

考虑一个单晶金属。与均匀的橡胶块不同,它的性质是高度​​各向异性​​的;它在某些方向上比其他方向更强。这种各向异性由其晶格中原子的排列决定。晶体的弹性性质不是在固定的实验室坐标系中保持不变,而是在一个随该晶格旋转的坐标系中保持不变。

突然之间,自旋的选择不再是任意的。物理情境要求我们使用一个与​​晶格本身​​同旋的客观率。我们率方程中的自旋ω\boldsymbol{\omega}ω找到了一个直接的物理意义:它必须是​​弹性自旋​​We\mathbf{W}^eWe,它描述了材料变形时晶格的旋转。 这是一个美妙的统一时刻。我们为了满足一个高层原理而引入的抽象数学概念——自旋,现在被识别为一个具体的、微观的物理过程。

这个思想延伸到材料如何“记忆”其历史。在塑性力学中,我们使用像​​背应力张量​​α\boldsymbol{\alpha}α这样的内变量来模拟屈服面在应力空间中的移动,这种现象称为随动硬化。这个背应力代表了材料结构中的微观残余应力。为了具有物理意义,它的演化也必须是客观的,并且必须随着它所代表的底层材料结构被正确地传递和旋转。 那个适用于应力的同旋框架,也为我们提供了处理这些内部记忆的工具。

统一的视角

我们的探索之旅从对观察者无关定律的抽象需求,走到了预测流体流动、确保能量守恒、编写稳健的计算机代码以及理解晶体微观行为的具体挑战。在整个过程中,同旋率的概念一直是我们忠实的向导。这证明了物理学非凡的力量和统一性,一个单一、优雅的思想能够提供一种连贯的语言,来描述跨越如此广阔的尺度和学科的现象,确保我们的模型不仅仅是数学上一致的,更是物理上真实的。