
我们如何超越对物体形状变化的简单描述,并用数学的精度来量化它?答案在于应变的概念,即材料变形的度量。虽然对于微小变化来说,这似乎很简单,但大变形的世界——比如在柔性电子设备、软生物组织和金属成形中看到的变形——提出了一个重大的挑战。在这里,简单的公式不再适用,事实证明,单一、普适的应变“标尺”是难以捉摸的。这种模糊性导致了各种各样的应变度量,每种度量都有其特定的目的和有效范围,形成了一个可能令人困惑的领域。
本文旨在作为一份全面的指南,揭开应变度量世界的神秘面纱。我们将弥合抽象理论与实际应用之间的鸿沟,为针对任何给定问题选择和应用正确的度量提供必要的清晰度。本文的结构旨在帮助您从头开始建立理解,从基本概念入手,逐步走向真实世界的场景。
首先,在“原理与机制”一节中,我们将剖析变形的数学核心。我们将探讨变形梯度,理解为分离拉伸与旋转而对客观性的迫切需求,并介绍主要的应变度量族,包括 Green-Lagrange 应变、Euler-Almansi 应变和“真实”对数应变。在这一理论基础之后,“应用与跨学科联系”一节将展示为什么这些区别至关重要。我们将看到,应变度量的选择如何影响从飞机机翼的安全分析和实验数据的解读,到工程、生物学和地球物理学等不同领域中复杂计算机模拟的稳定性等方方面面。
想象一下,你拿一块粘土,把它滚成一根又长又细的金属丝。你使它变形了。你改变了它的形状。在物理学和工程学的世界里,我们希望用数学的精度而不是模糊的词语来描述这种变化。我们该怎么做呢?
第一步是将物体看作点的集合。我们为物体在其初始、舒适的状态下拍一张“快照”——我们称之为参考构型。这个参考构型中的每一点都有一个“名字”或坐标,我们可以标记为 。现在,我们施加一些力,我们推拉它,物体移动并改变形状,进入一个新的状态,即当前构型。每个点 都移动到了一个新位置,我们称之为 。
整个变形过程由一个规则,一个函数来捕捉,它告诉我们每一个点去了哪里。我们称之为变形映射,,因此有 。这个映射完整地讲述了物体中每个粒子的运动历程。
但通常我们不关心整个历程,而是关心一个点附近的局部情况。材料在那里被拉伸或剪切了多少?要回答这个问题,我们需要观察当我们从 移开一小步时,映射 是如何变化的。这就引出了我们故事的主角:变形梯度,用符号 表示。它被定义为变形映射的梯度:
不要被这个符号吓到。你可以把 看作一个小的、局部的机器(数学上,它是一个矩阵或二阶张量)。如果你给它输入一个代表参考构型中无穷小纤维的微小向量 ,它就会输出当前构型中新的、变形后的向量 。这个小机器 包含了关于材料局部拉伸、剪切和旋转的所有信息。
现在,为了使这个图像完美运作,使得 在每一点都有良好定义,我们通常在经典物理学中假设我们的变形映射 是非常光滑且表现良好的——即所谓的 -diffeomorphism( 微分同胚)。然而,大自然并不总是那么井然有序。当我们用计算机模拟现实世界的现象,如金属撕裂或岩石断裂时,我们会遇到完全不光滑的情况。现代力学,特别是计算力学,已经发展出一种更强大、更现实的数学框架,使用 Sobolev 空间。这个框架允许可能存在扭结或尖角的映射,在“弱”意义上定义变形梯度 ,使其“几乎处处”成立,而不是处处成立。正是这种推广赋予了有限元法等方法模拟复杂的真实世界变形的强大能力。
这里我们遇到了一个至关重要的问题。想象一下建筑物中一根长的钢制工字梁。如果一台起重机将这根梁吊起,只是将它旋转到位而没有弯曲或拉伸它,那么这根梁“应变”了吗?你的直觉告诉你没有。它的长度没有改变,形状也一样。然而,梁中的每一点都移动了,变形梯度 肯定不是单位矩阵。在这种情况下, 纯粹是一个旋转矩阵,我们称之为 。
这个简单的思想实验揭示了任何合理的变形或应变度量的一个关键要求:它必须完全忽略刚体旋转。如果一个运动只包含旋转(和平移),我们的应变度量必须给出零的结果。这是一个基本原则,称为材料标架无关性,或更简单地称为客观性。
那么,我们如何从数学上从变形梯度 中“过滤掉”旋转,以得到“纯粹”的应变呢?大自然通过数学的语言,提供了一个惊人优雅的工具,称为极分解。它告诉我们,任何可逆的变形梯度 都可以唯一地分解为两种方式:
在这个分解中, 是一个正常正交张量,代表变形的纯旋转部分。另外两个张量, 和 ,是对称正定的,它们代表变形的纯拉伸部分。它们是应变的本质,剥离了任何旋转污染。
那么 和 是什么呢?
因为 存在于参考坐标系中,所以由它构建的度量被称为 Lagrangian 度量。因为 存在于当前的、空间的坐标系中,所以由它构建的度量被称为 Eulerian 度量。这两种“风格”的拉伸的存在,是我们即将探讨的不同应变度量族的根源。
现在我们已经分离出纯拉伸 和 ,我们如何量化它呢?事实证明,没有唯一的、天赐的答案。有很多方法可以定义一个“标尺”来测量应变,每种方法都有其优点。让我们通过一个简单的例子来探讨最著名的几种方法:一根杆从初始长度 拉伸到最终长度 。比率 是一维的“拉伸”。
工程应变: 这可能是你首先想到的。它是长度变化量除以原始长度:。它简单、直观,对于微小变化非常有用。
Green-Lagrange 应变: 这个度量定义为 。为什么是平方项?它自然地产生于右 Cauchy-Green 张量 。如果将 代入,你会发现 。由于 是我们客观的拉伸度量, 也是客观的,而由它构建的应变度量 也巧妙地继承了这种客观性。它是典型的 Lagrangian 度量,将变形后的几何状态与参考状态进行比较。
Euler-Almansi 应变: 这是 Green-Lagrange 应变在空间构型上的对应物,定义为 。它由左 Cauchy-Green 张量 构建。其张量定义为 。这是一个 Eulerian 度量,从最终的、当前构型的角度描述应变。
真实(或对数)应变: 这个度量有一种特殊的吸引力。想象一下不是一次性拉伸杆,而是一系列微小的增量。在每一步,增量应变是微小的长度变化 除以当前长度 。真实应变就是你将所有这些微小的分数变化加起来得到的结果:。它的张量形式,Hencky 应变 ,通常被认为是弹性应变最“自然”的度量。
对于小变形,其中 非常接近 1,所有这些不同的标尺给出的答案几乎相同。如果你用泰勒级数展开它们,你会发现它们都以 开始,仅在高阶项上有所不同。这就是为什么对于许多日常工程问题,简单的工程应变就足够了。
但对于“有限”或大应变呢?让我们以将材料拉伸 20% 为例,所以 。Green-Lagrange 应变为 ,而对数应变为 。差异约为 21%!这不仅仅是学术上的好奇心。如果你正在用一种可以大幅拉伸的现代聚合物或生物组织设计一个零件,选择错误的应变度量可能会导致危险的错误预测。
那么,为什么我们会有这么多的应变度量呢?第一个原因,正如我们所见,是客观性。任何由拉伸张量 或 构建的度量(如 Green-Lagrange 应变、Euler-Almansi 应变或 Hencky 应变)都能正确地忽略刚体旋转,并在纯旋转情况下给出零值,这是必须的。通常在入门课程中教授的简单无穷小应变张量 ,是出了名的不客观。它混淆了旋转和应变,对于纯粹的刚体旋转会给出非零的“应变”。它仅在应变和旋转都非常小的情况下才可靠。
第二个更深层次的原因在于应力与应变之间错综复杂的共舞。应变度量的选择并非独立于你如何度量应力。为了使本构律(描述材料应力与应变关系的定律)具有物理意义,应力和应变度量必须是功共轭的。这意味着应力张量与应变张量率的乘积给出了每单位体积耗散或储存的功率。不同的应变度量在这场共舞中有不同的应力伙伴:
这就是为什么不同领域使用不同的配对。计算固体力学代码可能会使用 对在其参考坐标系中进行计算。而进行拉伸试验的材料科学家可能会绘制 Cauchy 应力与对数应变的关系图,因为它通常能更直接地揭示材料的行为。
这就引出了一个常见且具有挑战性的问题:分析经历大旋转但只有小应变的结构,例如弯曲的钓鱼竿或旋转的直升机桨叶。在这里,无穷小应变张量 是无用的。必须使用像 这样的客观度量。一个非常聪明的替代方法是使用协同旋转列式(corotational formulation)。这是一种计算技巧,在变形的每一步,你都在数学上将单元“反向旋转”回一个中性方向,并在那个旋转的参考系中测量小应变。它巧妙地将大的、几何上复杂的旋转与小的、物理上简单的应变分离开来。
有限应变的世界具有优美的数学结构,但它要求严谨和细心。即使是像 Mohr 圆这样分析无穷小应变的常用工具,也不能被轻率地应用。我们不能简单地为总的 Green-Lagrange 应变张量画一个 Mohr 圆来寻找当前构型中的主方向。然而,我们可以在给定的瞬间正确地将其应用于变形率 ,因为你是在一个单一、固定的空间坐标系内分析一个无穷小的变化。这表明,来自线性世界的概念可以被小心翼翼地、正确地应用于大变形的非线性世界,展示了连续介质力学的深刻一致性和优雅。
在遍历了应变的基本原理之后,我们可能会觉得我们的探索已经完成。我们已经定义了术语,并考察了数学机制。但这就像只学习了国际象棋的规则却从未下过一盘棋。一个科学概念的真正美和力量,不是体现在其定义中,而是在其应用中。为什么我们需要这么多不同的方法来测量变形?什么时候简单的近似就足够了,什么时候它又会把我们引向灾难性的歧途?
这些问题的答案将我们从宏大的工程项目中心带到活细胞的精妙之舞,从地壳的挤压深处带到计算机模拟的无声世界。应变不仅仅是教科书中的一个主题;它是一种通用的语言,用于探究、预测和塑造物理世界。现在让我们来探索这个世界,看看“应变度量”的选择如何成为一个具有深远影响的重大决定。
最直接地说,应变是一个安全问题。想象一下飞机机翼在湍流中弯曲。工程师首要关心的是材料的应力没有超过其极限。但我们怎么知道呢?我们看不到应力。然而,我们可以看到它的效果:变形。通过将一个小型、灵敏的电阻应变片粘合到金属梁的表面,我们可以以惊人的精度测量局部的拉伸。一个读数,比如说 (或 350“微应变”),可能看起来微不足道,但通过知道材料的杨氏模量——其固有的刚度——我们可以使用胡克定律将这个微小的拉伸直接转化为维持机翼整体的内部拉应力的度量。这个简单的想法,即在测量的应变和计算的应力之间建立直接的桥梁,是所有结构完整性分析的基石,从摩天大楼到航天器皆是如此。
只要变形保持微小且材料能弹回原状,这种方法就非常有效。但是,如果我们想知道材料的最终极限呢?当我们拉伸一根金属棒直到它永久变形并开始“颈缩”然后断裂时,会发生什么?在这里,我们简单的近似开始失效。随着杆的拉伸,其横截面积会缩小。为了准确地描述材料真正承受的应力,我们必须将力除以瞬时面积,而不是原始面积,这个瞬时面积是不断变化的。
同样,应变的定义本身也变得模棱两可。如果一根10米长的杆伸长了1米,“工程应变”是 。但如果我们把它看作是两根5米长的杆,每根伸长了半米,那么每根的应变是 。这似乎没问题。但一个真正有洞察力的度量,即对数或真实应变,认识到杆的每个无穷小部分都是相对于其当前长度拉伸的。通过积分这些微小的变化,我们得到了一个度量,,它正确地捕捉了变形的累积过程。对于经历大塑性变形的材料,我们可以假设其体积保持不变,此时真实应变和瞬时面积之间有着优美而密不可分的联系。进行拉伸试验以表征新合金的材料科学家必须使用这些更复杂的度量——真实应力和真实应变——来建立远超弹性极限的精确材料行为模型。
测量应变是一门精妙的艺术。我们选择的仪器反映了不同的测量哲学,每种都有其优缺点。夹式引伸计是一种机械装置,它在两个点上物理夹持试样,提供了这两点之间应变的可靠平均度量。它非常直接,并且因为它只关心其接触点之间的距离,所以它完全不受任何刚体运动的影响——如果整个试样摆动或旋转,引伸计的读数不受影响。应变片,正如我们所见,提供了其粘合的小块表面上平均应变的电气读数。它也基本上对刚体运动不敏感,因为它只响应其所附着的材料的变形。
一项革命性的现代技术是数字图像相关法(DIC)。通过追踪“之前”和“之后”图像之间表面上随机散斑图案的运动,DIC可以生成一个完整的、连续的位移场图。从这个位移图中,通过进行空间求导来计算应变场。这种方法非常强大,提供了单点测量仪无法想象的丰富数据。
然而,这种强大功能伴随着一种深刻的微妙之处,这正触及了我们需要不同应变度量的核心原因。想象一个平面物体仅仅旋转了,比如说, 度,而形状没有任何改变。应变是多少?物理上,它必须是零。像 Green-Lagrange 应变这样的标架无关度量在这种情况下正确地报告为零。但一个线性化的应变度量——那种只是位移场的简单导数——会变得非常混乱。它看到物体上的点在移动,并错误地将这种运动解释为均匀的压缩!。这不是一个小错误;这是测量哲学的根本性失败。如果 DIC 使用线性化应变计算进行天真的处理,它会为一个简单的旋转报告巨大的、虚假的应变。这迫使我们面对“几何非线性”,不是作为一个数学上的好奇心,而是作为一个避免被我们自己的工具所愚弄的物理必要性。
此外,从 DIC 中的位移数据计算应变是一种权衡。为了得到某一点的应变,我们必须在某个小区域上对位移场进行微分。如果我们选择一个大区域,我们就会平均掉细节,失去空间分辨率。如果我们选择一个小区域,我们的计算就会对位移测量中的任何噪声变得极其敏感,因为微分天生会放高频噪声。这给实验者带来了经典的偏差-方差困境:减少随机误差的代价是模糊了人们希望看到的特征。
这些细致实验的数据本身并非终点。它们是构建能够预测复杂系统行为的计算模型的原始成分。为了模拟延性金属,工程师使用像 von Mises 屈服准则这样的框架;对于像土壤这样的颗粒材料,他们可能会使用 Mohr-Coulomb 准则。但这些模型有参数——杨氏模量、泊松比、屈服应力、内聚力、摩擦角——这些都不是先验已知的。它们必须被测量。一套最基本的实验室测试,比如对金属的单轴拉伸试验和对土壤的一对三轴压缩试验,提供了校准这些模型并在计算机内使其生效所需的应力-应变数据。
在计算机内部,在有限元法(FEM)的世界里,应变度量的选择具有深远的影响。最简单的结构单元,桁架杆,通常使用无穷小应变 来建模。然而,这仅仅是更完整的、几何非线性的 Green-Lagrange 应变 的一阶近似。这为什么重要?如果一个结构经历大旋转(想象一根柔性钓鱼竿弯曲),简单的线性近似会失效,原因与它在我们的 DIC 例子中失效的原因相同:它无法区分旋转和应变。许多软件包中使用的一个巧妙解决方案是“协同旋转”列式。模拟会单独跟踪单元的大刚体旋转,然后在单元自己的、局部的、旋转的参考系中使用简单的小应变理论。这是一个绝妙的技巧,它兼具了计算的简便性和物理的准确性。
对于涉及真正巨大变形的问题,比如金属零件的锻造或软橡胶的建模,应变度量的选择直接影响模拟的稳定性和效率。基于 Green-Lagrange 应变的模型可能导致计算出的刚度随拉伸呈二次方增长。在模拟中,这意味着模型中高度拉伸的部分可能在数值上变得“刚硬”,导致迭代的 Newton-Raphson 求解器难以收敛或失败。相比之下,使用对数应变(就是我们在材料测试中遇到的那个“真实应变”)制定的模型,会得到一个表现好得多的刚度,即使在极端变形下也保持有界。这可能决定了一个模拟是平滑收敛还是完全失败。在理论的最高层次,当模拟像高温下塑性这样的复杂行为时,其框架是由热力学定律决定的。有限应变塑性的标准模型需要变形梯度的乘法分解,,这个概念比小应变理论中使用的简单加法分解要复杂得多。
也许应变最令人惊奇的方面是它的原理如何将看似毫不相干的科学领域统一起来。
想想卑微的蚯蚓。它通过收缩肌肉来改变形状,从而移动,经历巨大的轴向缩短和径向膨胀——40%或更多的变形很常见。希望模拟这种运动的生物学家必须同时应对大拉伸和大刚体旋转,因为蠕虫会转弯。如果他们使用线性化的应变模型,他们将陷入我们之前指出的陷阱:模型会预测蠕虫每次简单地旋转其身体节段时都在承受应力。为了正确模拟这种静水骨骼的物理特性,必须采用一种标架无关的、有限应变度量,如 Green-Lagrange 应变。几何非线性不是一个学术注脚;它是蠕虫如何爬行的基本物理学。
现在,从蠕虫的尺度转换到行星的尺度。一位模拟脆弱上地壳的地球物理学家想要了解断层和地震。在地质时间尺度上,刚体运动是巨大的——大陆漂移数千公里,岩体发生显著旋转。然而,岩石在断裂前能够承受的实际弹性应变非常小,通常不到百分之一。这里我们有一个悖论:巨大的旋转,但微小的应变。一个完整的有限应变公式是正确的,但计算成本高昂。一个简单的小应变公式计算成本低,但由于大旋转而物理上是错误的。解决方案与我们在计算工程中找到的相同:一个协同旋转框架。通过使用一个考虑了旋转的客观应力率,地球物理学家可以准确地模拟在一个随缓慢搅动的构造板块一起旋转的参考系内,那些引起断裂的小应变。
从确保飞机安全的工程师,到解码运动的生物学家,再到模拟行星的地球物理学家,应变的概念提供了一种通用的、强大的语言。从简单的工程应变到连续介质力学的完整张量体系的旅程,是一个直面复杂性的故事。我们被迫采用更复杂的度量,不是为了数学本身,而是因为物理世界——在其转动、拉伸和流动中——要求我们这样做。