try ai
科普
编辑
分享
反馈
  • 非对称应力张量

非对称应力张量

SciencePedia玻尔百科
核心要点
  • 经典应力张量的对称性是简单、无结构连续体中角动量守恒的直接结果。
  • 当需要模拟具有内部微结构并能承受内力矩的材料(如泡沫、复合材料和颗粒介质)时,应力张量的非对称性就变得至关重要。
  • 微极(Cosserat)理论引入了新的物理量,如耦合应力和微旋转,以解释这些内力矩和旋转。
  • 非对称应力理论在解决有限元法中的数值问题以及解释手性超材料的独特性质方面具有实际应用。

引言

在经典连续介质力学中,应力张量的对称性是一项基本原则,它直接源自角动量守恒定律。这一优雅的概念完美地描述了简单、均匀材料内部的内力。然而,当面对内部结构对其行为起着至关重要作用的现代材料时,这个经典模型就达到了其极限。这就提出了一个关键问题:当一种材料的内部结构可以承受局部力矩时会发生什么?我们又该如何描述这类系统的物理学?本文将深入探讨非对称应力张量的世界,为从经典理论过渡到更全面的框架搭建一座桥梁。第一章“原理与机制”将剖析支持对称性的经典论证,然后打破它,引入微极理论的核心思想,如体力偶和耦合应力。接下来的“应用与跨学科联系”一章将展示这一看似抽象的理论如何在模拟复杂流体、设计先进复合材料,乃至解决计算工程中的长期挑战方面找到强大而实际的应用。

原理与机制

对称性的神圣性

在经典力学的世界里,有些真理几乎是神圣不可侵犯的。其中最为优雅的一条就是应力张量的对称性。当我们描述材料内部——一块钢、一柱水、这个房间里的空气——的内力状态时,我们使用一个称为​​Cauchy应力张量​​的数学对象,可以写作 σ\boldsymbol{\sigma}σ。你可以把它想象成一台机器,它能告诉你,在材料内部任意做一个假想切割时,该切面上的力矢量(面力)是多少。

那么,为什么这个张量必须是对称的呢?为什么应力分量 σ12\sigma_{12}σ12​(其法线方向为 x2x_2x2​ 的面上的 x1x_1x1​ 方向的力)必须等于 σ21\sigma_{21}σ21​(其法线方向为 x1x_1x1​ 的面上的 x2x_2x2​ 方向的力)呢?这并非为了数学上的便利,而是自然界基本定律——​​角动量守恒​​——的一个深刻推论。

想象一个漂浮在空间中的微小、无限小的材料立方体。周围材料的应力正在其六个面上拉伸和挤压。如果我们有 σ12≠σ21\sigma_{12} \ne \sigma_{21}σ12​=σ21​,那就意味着我们的小立方体上存在一个净扭转力,即​​力矩​​。σ12\sigma_{12}σ12​ 和 σ21\sigma_{21}σ21​ 这对力就像微小的手指试图让它旋转。如果立方体内部没有任何东西可以抵抗这个力矩,会发生什么?根据牛顿定律,这个立方体必须开始旋转。而且因为这个立方体是无限小的,它的转动惯量也是无限小的,所以它将以无限大的角加速度旋转!这当然是物理上荒谬的。自然界不允许在简单的“经典”材料中出现这种行为。

为了避免这场灾难,这些应力对立方体产生的净力矩必须为零。这直接迫使应力张量必须是对称的:σij=σji\sigma_{ij} = \sigma_{ji}σij​=σji​。任何非对称部分都会产生一个无抵抗的力矩,从而违反角动量守恒。

事实证明,应力张量中负责产生这个力矩的恰好是它的​​反对称部分​​,即 σskew=12(σ−σT)\boldsymbol{\sigma}^{\text{skew}} = \frac{1}{2}(\boldsymbol{\sigma} - \boldsymbol{\sigma}^{\mathsf{T}})σskew=21​(σ−σT)。实际上,由应力场产生的局部力矩密度(单位体积的力矩)可以用一个非常简洁的表达式表示,其第 kkk 个分量为 mk=ϵijkσijm_k = \epsilon_{ijk}\sigma_{ij}mk​=ϵijk​σij​,其中 ϵijk\epsilon_{ijk}ϵijk​ 是Levi-Civita符号,它优雅地处理了力矩的类叉乘性质。对于处于平衡状态的经典连续体,这个力矩密度必须为零。这再次告诉我们,应力张量必须是对称的。

微小扳手的世界:打破对称性

但是,如果我们处理的不是一个简单的、经典的连续体呢?物理学的进步常常始于“如果……会怎样?”的提问。如果一种材料能够承受一个净内力矩呢?如果我们那个无限小的立方体有某种方式可以抵抗被旋转呢?

想象一下,我们的材料不是均匀的果冻,而是充满了无数微观的磁性颗粒。现在,假设我们施加一个外部磁场。每个微小的磁铁都会试图与磁场对齐,在此过程中,它会感受到一个微小的力矩。现在,整个材料包含了一个分布式的力矩,一片在每一点上都在扭转的“微小扳手海洋”。我们称之为​​体力偶​​。

在这种情况下,一个非对称的应力张量就不再是自相矛盾的了。由应力张量非对称性产生的力矩可以与抵抗它的内部体力偶完美平衡。如果应力试图使一个单元顺时针旋转,而体力偶以同等大小的力矩试图使其逆时针旋转,那么该单元将保持旋转平衡。

对于静态情况,这种平衡可以用一个简洁的方程来描述:ϵijkσjk+ci=0\epsilon_{ijk} \sigma_{jk} + c_i = 0ϵijk​σjk​+ci​=0,其中 cic_ici​ 代表体力偶矢量的分量。一个非零的体力偶 cic_ici​ 允许,甚至要求,应力张量存在一个非零的反对称部分。例如,如果一种假想材料的应力状态为

σ=(120805095150756585100)Pa\boldsymbol{\sigma} = \begin{pmatrix} 120 & 80 & 50 \\ 95 & 150 & 75 \\ 65 & 85 & 100 \end{pmatrix} \text{Pa}σ=​1209565​8015085​5075100​​Pa

这种状态只有在存在一个体力偶密度 c=(10,−15,15) N/m2\mathbf{c} = (10, -15, 15) \text{ N/m}^2c=(10,−15,15) N/m2 来维持旋转平衡时,才是物理上可能的。应力张量的对称性并非所有物质的普适定律,而是一种特定、简单物质模型的结果。通过放宽这个模型,我们可以描述一个更丰富的世界。

更丰富的现实:微极世界

材料具有能够支持力偶的内部结构这一想法,并不仅仅是理论上的幻想。它是我们称之为​​微极​​(micropolar)或​​Cosserat理论​​的基础,该理论以一个多世纪前提出它的杰出的Cosserat兄弟命名。这个理论对于精确模拟那些材料“点”本身具有生命力的材料至关重要——例如颗粒状固体(如沙子)、泡沫、某些液晶,甚至骨骼。

在经典理论中,一个材料“点”就只是一个点:空间中一个没有特征的点,其状态由其位置描述。微极理论的基本假设是丰富这一图景。一个微极“点”被想象成不仅有位置,还有姿态。它就像一个微小的刚体。想想一个由滚珠轴承组成的流体与一种光滑、连续的液体。每个滚珠轴承不仅可以移动,还可以独立于其邻居旋转。

这赋予了我们的材料一个新的、独立的运动学自由度:​​微旋转矢量​​,通常用 φ\boldsymbol{\varphi}φ 或 χ\boldsymbol{\chi}χ 表示。这与你从速度场旋度中得到的熟悉的宏观旋转(涡量)不同;它是材料内部状态一个真正全新的、独立的属性。

正如力引起运动(位移和速度)的变化一样,力矩引起旋转的变化。这个新的旋转自由度带来了一整套新的物理量:

  • ​​耦合应力​​:正如应力(σ\boldsymbol{\sigma}σ)是单位面积传递的力,​​耦合应力​​(m\boldsymbol{m}m 或 μ\boldsymbol{\mu}μ)是单位面积传递的力矩。它描述了相邻微结构如何相互扭转。
  • ​​微惯量​​:正如质量(ρ\rhoρ)抵抗线性加速度,​​微惯量​​(jjj)抵抗这些内部结构的角加速度。它是材料微结构的转动惯量。

新的力与力矩平衡

有了这些新元素,我们可以写下一个新的、更完整的角动量平衡定律。它美轮美奂,证明了一个自洽的物理理论如何扩展以包含新现象。某一点的力偶局部平衡形式为:

ϵijkσjk+mji,j+ci=ρjχ¨i\epsilon_{ijk} \sigma_{jk} + m_{ji,j} + c_i = \rho j \ddot{\chi}_iϵijk​σjk​+mji,j​+ci​=ρjχ¨​i​

让我们逐项分析这个宏伟的方程。

  • ϵijkσjk\epsilon_{ijk} \sigma_{jk}ϵijk​σjk​:这是由力-应力张量的非对称性产生的力矩。它是试图旋转材料点的“引擎”。
  • mji,jm_{ji,j}mji,j​:这是耦合应力张量的散度。它代表了因周围耦合应力变化而作用在材料单元上的净力矩。它是将力矩从一点传递到另一点的内部“传动装置”。
  • cic_ici​:这是单位体积的外部体力偶,就像作用在我们磁性颗粒上的力矩。它是“外部驱动器”。
  • ρjχ¨i\rho j \ddot{\chi}_iρjχ¨​i​:这是微惯量项。它是对微结构角加速度的抵抗,是储存或释放自旋角动量的“飞轮”。

即使在没有外部体力偶和耦合应力的动态情况下,仅仅是微惯量的存在就可能导致应力张量变得非对称!如果微结构正在加速或减速旋转(χ¨i≠0\ddot{\chi}_i \neq 0χ¨​i​=0),自旋惯量 ρjχ¨i\rho j \ddot{\chi}_iρjχ¨​i​ 必须由力-应力产生的力矩 ϵijkσjk\epsilon_{ijk} \sigma_{jk}ϵijk​σjk​ 来平衡,这意味着 σ\boldsymbol{\sigma}σ 不可能是对称的。

这个丰富的框架从能量角度来看也是完全自洽的。​​虚功率原理​​表明,每一种作用都有其功共轭的伙伴。面力(ti=σjinjt_i = \sigma_{ji} n_jti​=σji​nj​)对位移(uiu_iui​)做功,而新引入的力偶面力(qi=mjinjq_i = m_{ji} n_jqi​=mji​nj​)对微旋转(ϕi\phi_iϕi​)做功。没有任何东西是孤立的;这一切都融入了一个关于能量和功的连贯结构中。

最后一个疑问:“主”应力何去何从?

任何学习经典力学的学生都会学到主应力:对于任何应力状态,都存在三个相互正交的平面,在这些平面上,力是纯法向的(没有剪切)。这产生了三个实数主应力值和三个正交的主方向。其数学基础是​​谱定理​​,该定理之所以适用,是因为经典应力张量是对称的。

当我们允许 σ\boldsymbol{\sigma}σ 非对称时,会发生什么?如果我们天真地应用相同的数学定义——寻找张量 σ\boldsymbol{\sigma}σ 的特征向量——我们就会遇到麻烦。对于一般的非对称张量,其特征值不保证是实数;它们可以是复数!而其特征向量,即使是实数,也不保证是正交的。

“复数应力”到底意味着什么?这种数学上的怪异现象是一个信号,一个标志,表明我们可能问错了问题。对于一个非对称张量,代数特征值问题可能不是物理上最有洞察力的问题。

让我们问一个更好的、更物理的问题。与其问“在哪个方向上面力矢量与法线平行?”,不如问“在哪个平面上,面力的法向分量达到最大或最小值?”这是一个良定的变分问题,其解极具启发性。 任何张量都可以分解为对称部分和反对称部分:σ=σsym+σskew\boldsymbol{\sigma} = \boldsymbol{\sigma}^{\text{sym}} + \boldsymbol{\sigma}^{\text{skew}}σ=σsym+σskew。当我们计算法向面力 n⋅σn\mathbf{n} \cdot \boldsymbol{\sigma} \mathbf{n}n⋅σn时,反对称部分完全消失了!这是因为反对称部分只产生有助于力矩的剪切力。拉伸和压缩完全由对称部分处理。

因此,使法向面力取极值的方向,正是应力张量对称部分 σsym\boldsymbol{\sigma}^{\text{sym}}σsym 的主方向。并且由于 σsym\boldsymbol{\sigma}^{\text{sym}}σsym 根据定义是对称的,它总是有实数特征值和正交特征向量。主应力和主方向的经典概念得以恢复,但它只适用于应力中负责法向力的那一部分。另一部分,即反对称部分,有着完全不同的工作:产生力矩,它在与材料内部结构的优美动态平衡中完成这项工作。我们再次看到了物理学的统一与优雅,即使在更复杂的理论中,各个部分的角色依然清晰分明。

应用与跨学科联系

现在我们已经了解了非对称应力和它们独特的伙伴——耦合应力这个奇特的新世界,你可能会问:“这仅仅是一场有趣的数学体操,还是自然界真的在玩这个游戏?” 这是一个合理的问题。源于角动量守恒的应力张量对称性原理,是经典力学最优雅、最稳固的支柱之一。提出它对于某些材料可能是不完整的,这听起来近乎异端。

然而,正如物理学中常有的情况,当我们认为我们已经拥有了最终、完美的定律时,大自然会向我们展示宇宙的一个角落,那里的规则有着微妙而奇妙的不同。关键在于知道去哪里寻找。经典理论对于那些在我们关心的尺度上“无结构”的材料——一块钢、一杯水——工作得非常出色。但是,当材料本身充满了可以翻滚、扭转和旋转的小结构时,会发生什么呢?如果材料更像一盒微小的、相互啮合的齿轮,而不是光滑、均匀的果冻呢?

正是在这里,经典假设开始动摇,非对称应力的物理学也活跃起来。那个引导我们得出对称应力张量的经典论证,依赖于一个将体积缩小到一点的极限过程。它含蓄地假设,在我们缩小的过程中,不会出现任何新奇而奇怪的东西。但对于某些材料,这并非事实。考虑一块开孔泡沫,它是由相互连接的支柱和巨大空隙组成的混乱体。如果我们的“无限小”体积只有一个孔洞那么大,那么在一点上平滑变化的应力这一概念本身就崩溃了。或者想一想现代的“超材料”,一种为具有奇异特性而设计的人造结构。如果这些结构是由能够将力矩传递给邻居的微小单元构成的——比如由手性(即具有“旋向性”)螺旋组成的晶格——那么一种新的力矩传递方式,即耦合应力,就变得重要起来。即使在纳米尺度下,在一束微小的硅梁中,由表面张力产生的力也可以像作用在我们假想控制体边缘的力线一样,而经典理论则忽略了这种可能性。正是在这些结构化材料中——泡沫、颗粒介质、复合材料、液晶、骨骼和工程超材料——非对称应力理论找到了它的归宿。它不是对经典理论的替代,而是对其必要而美丽的扩展。

扭曲的流体与旋转的固体

那么,放宽对称性条件 σij=σji\sigma_{ij} = \sigma_{ji}σij​=σji​ 最直接的物理后果是什么?这意味着在材料内部,存在着一个内在的力矩源。想象一个微小的材料立方体。如果顶面的剪应力 σyx\sigma_{yx}σyx​ 不等于侧面的剪应力 σxy\sigma_{xy}σxy​,立方体将经历一个使其想要旋转的净力矩。在经典流体中,这在平衡状态下是不可能的。但在微极流体中——一种其微观组分(如细长分子)可以旋转的流体——这种内力矩可以存在,前提是它被其他东西所平衡。这个“其他东西”就是我们所说的​​体力偶​​,一种可以对材料中每一点施加纯力矩的外部场(也许是电磁场)。

这不仅仅是一个抽象的平衡动作,它具有切实的效应。想象一个浸没在具有恒定但非对称应力状态的假想流体中的球体。在经典流体中,恒定的应力(如均匀压力)只能压碎球体,永远无法使其转动。但在我们这种奇怪的流体中,σxy\sigma_{xy}σxy​ 和 σyx\sigma_{yx}σyx​ 之间以及其他分量对之间的不匹配,在球体表面上累加起来,产生一个宏观的净力矩。这种流体甚至在没有流动的情况下,就能使球体旋转!这就是液晶、铁磁流体和其他内部结构至关重要的“复杂流体”的世界。

此外,故事甚至更加丰富。在这些材料中,不仅剪应力可以是非对称的,材料还可以在表面上传递纯力矩。这便产生了​​耦合应力张量​​ M\mathbf{M}M。可以这样想:标准的应力张量 σ\boldsymbol{\sigma}σ 告诉你表面上单位面积的力。耦合应力张量 M\mathbf{M}M 告诉你单位面积的力矩。如果你将一个物体浸没在具有旋转微结构的流体中,这种耦合应力可以对物体施加其自身的力矩,这与由非对称Cauchy应力引起的力矩是分开的。就好像流体颗粒不仅在推挤物体,还在主动地抓住并扭转它。

设计未来材料

这个框架不仅用于描述奇异的流体,它还是设计未来固体材料的强大工具。许多先进材料,从纤维增强复合材料到轻质金属泡沫,其强度和特殊性能都源于其复杂的内部微结构。对于一个大型物体来说,对每一根纤维或泡沫支柱进行建模在计算上是不可能的。这正是微极(或在固体中称为​​Cosserat​​)理论大放异彩的地方。它提供了一种“弥散化”的连续介质描述,捕捉了微结构的平均效应,包括其旋转能力。

考虑一种通过将微小的刚性球形颗粒嵌入柔软的弹性基体中制成的复合材料。由爱因斯坦本人开创的经典复合材料理论告诉我们,这些颗粒如何通过抵抗剪切来增加材料的刚度。但是,如果这些颗粒可以自由旋转,就像基体中的微小滚珠轴承一样呢?Cosserat模型可以解释这一点。它将材料视为一个点既可以平移也可以旋转的连续体。有趣的结果是,材料的整体行为取决于其组分的内部自旋。例如,复合材料的有效剪切模量最终是经典刚化效应和与嵌入颗粒微旋转相关的新效应的组合。该理论正确地直觉到,如果颗粒可以旋转,它会改变整个材料对扭曲的响应方式。

该理论还为审视金属强度和延展性的核心——位错——提供了一个更精细的视角。位错是晶格中的一种线缺陷,其运动是金属能够发生塑性变形的原因。经典弹性理论为位错周围的应力场提供了一个很好的初步近似,但它预测在位错核心处应力为无穷大,这在物理上是不现实的。Cosserat理论提供了一个漂亮的修正。通过引入一个与材料微结构(可能是晶粒尺寸或颗粒尺寸)相关的特征长度尺度,它使应力场正规化。材料“点”的旋转能力为适应位错核心处的极端变形提供了一种新方式,从而抹平了应力集中。由此产生的应力场不再像经典理论所要求的那样对称,反映了缺陷附近变形的复杂旋转性质。这弥合了晶格的离散世界与连续介质力学的平滑世界之间的鸿沟。

从物理理论到数字现实

这些思想的实际效用直接延伸到了计算工程领域。在有限元法(FEM)中,工程师通过将汽车、飞机和建筑物分解成小单元来建立它们的虚拟模型。对于像车身面板这样的薄结构,他们经常使用“壳”或“膜”单元。一个概念上的难题长期困扰着这些单元:“钻孔自由度”。想象一下模拟中的一个平面三角形单元。它具有上下、左右移动和倾斜的自由度。但是,关于在其自身平面内像风车一样旋转的自由度呢?这就是钻孔旋转。

在经典力学中,没有与这种纯旋转相关的物理能量。因此,在有限元模拟中,刚度矩阵可能变得奇异,导致灾难性的数值错误。多年来,计算力学家使用各种数学技巧或“罚函数法”来人为地增强这种旋转的刚度,使他们的模拟能够工作。但Cosserat理论为解决这个问题提供了一个严谨的、物理的基础。通过将壳体建模为二维微极连续体,钻孔旋转不再是数学上的虚构,它变成了材料点的物理微旋转。这种旋转受到材料耦合应力的抵抗。因此,该理论自然地为钻孔旋转提供了物理刚度,从而稳定了模拟,而无需任意的修正。这个曾被视为奇异和抽象的理论,在现代工程设计的核心找到了一个至关重要且实际的应用。

手性与对称性破缺

非对称应力理论所揭示的最深远的联系,或许与自然界最基本的属性之一有关:手性,或称“旋向性”。许多结构,从DNA的双螺旋到某些海贝壳的优雅螺旋,都是手性的——它们与它们的镜像不完全相同。这种结构不对称性也可以构建到人造超材料中。

在经典线性弹性理论中,有一条深刻而优美的原理,称为Betti互易定理。简单来说,它指出,对于一个给定的弹性体,一组力作用于由第二组力引起的位移上所做的功,等于第二组力作用于由第一组力引起的位移上所做的功。这是材料响应中深刻对称性的体现。

令人震惊的是,在手性微极材料中,这种神圣的对称性可以被打破。材料的旋向性在剪切和弯曲之间引入了一种不对称的耦合。例如,剪切材料可能会使其扭曲,但施加纯扭曲可能不会引起相应的剪应力。反映材料手性几何形状的非对称本构关系,直接导致了Betti定理的违背。当我们对这种材料的两种不同加载情况进行计算时,我们发现交叉功项不再相等。材料微观结构中镜像对称性的缺失,表现为其宏观响应中互易性的缺失。

这是一个完美的例子,说明了物理学为何如此引人入胜。微观尺度上的几何属性(手性)被本构关系中的一个数学特征(一个非对称矩阵)所捕捉,这反过来又导致了宏观尺度上一个深刻对称性原理(互易性)的破缺。非对称应力张量理论,最初只是一个奇特的智力练习,因此成为连接材料结构与物理世界基本对称性的桥梁。它提醒我们,每当我们推动材料的边界时,我们也必须准备好扩展我们物理定律的边界。