try ai
科普
编辑
分享
反馈
  • 非局部力学

非局部力学

SciencePedia玻尔百科
核心要点
  • 非局部力学通过假设某一点的应力取决于其整个邻域内的应变场,而不仅仅是该点本身,从而重新定义了材料行为。
  • 它通过引入一个特征性的内禀长度尺度,解决了在裂纹尖端和位错核心处出现无限应力的非物理预测。
  • 该理论解释了纳米技术中的尺寸依赖现象,例如纳米梁振动频率的降低以及石墨烯等材料中的波色散。
  • 作为一种广义连续介质理论,非局部力学在离散原子模型和经典连续介质力学之间架起了一座桥梁。

引言

几个世纪以来,经典连续介质力学一直是工程学的基石,使我们能够建造桥梁、设计飞机,并以惊人的精度理解宏观世界。其核心假设是局部性:材料在任意一点的行为仅取决于该精确点的条件。然而,当我们的科学雄心进入纳米尺度,并探究材料完整性的极限时,这个优雅的简化开始显现出其裂痕。经典理论预测在裂纹尖端存在非物理的无穷大值,也无法解释为什么10纳米的金属丝与其宏观对应物的行为不同。我们理解上的这一差距要求我们对物质进行更深层次的描述——一种承认材料内部结构固有相互关联性的描述。

本文介绍​​非局部力学​​,这是一个革命性的框架,它通过纳入这些至关重要的长程相互作用来扩展经典理论。通过超越局部性的限制,它为我们观察力学世界提供了一个更真实、更强大的视角。在接下来的章节中,我们将对这个引人入胜的课题进行全面探索。第一章​​“原理与机制”​​将奠定理论基础,解构非局部本构律、内禀长度尺度的作用,及其与其他广义连续介质理论的关系。随后,关于​​“应用与跨学科联系”​​的章节将展示该理论的实际应用威力,演示它如何抑制断裂力学中的无穷大值,并为纳米材料的动态行为解锁新的见解。

原理与机制

点的终结:一个充满邻居的宇宙

几个世纪以来,我们对物理世界的描述一直是深刻地局部的。想象一下,你想知道一根钢梁内部某一点的应力——即内部的推与拉。经典物理学,沿袭 Isaac Newton 及其后继者的传统,给出了一个极其简单的答案:该精确点的应力仅取决于该完全相同点的应变(局部变形)。在这个宇宙里,每个点都是一个自恋者,只关心自己的状态。这里的应力由这里的应变决定。仅此而已。这就是​​连续介质假设​​,它是一个极其强大的思想,曾用于建造桥梁和飞行飞机。

但如果我们看得更近一些呢?如果我们放大到能看到构成钢的基石——原子呢?单个原子并非存在于真空中。它感受到邻居的推和拉。它的存在状态是与整个邻域集体协商的结果。那么,为什么我们的连续介质描述就应该有所不同呢?

这就是​​非局部力学​​的出发点。它提出了一个激进但又非常直观的偏离:一个点的应力不仅仅是局部应变的函数,而是所有周围点应变的加权平均。这是一个社群的世界,而非孤立个体的世界。一个点能感受到其整个材料邻域的影响,尽管这种影响会随着距离的增加而减弱。

在数学上,这个优雅的思想表现为一个积分形式。我们不再使用像 σ=Eε\sigma = E \varepsilonσ=Eε 这样简单的代数关系,而是拥有了更丰富的内容。点 xxx 处的非局部应力 σij\sigma_{ij}σij​ 是通过对物体中所有其他点 ξ\xiξ 的贡献进行积分得到的:

σij(x)=∫Ωα(∣x−ξ∣,ℓ) Cijkl εkl(ξ) dξ\sigma_{ij}(x) = \int_{\Omega} \alpha(|x-\xi|, \ell) \, C_{ijkl} \, \varepsilon_{kl}(\xi) \, d\xiσij​(x)=∫Ω​α(∣x−ξ∣,ℓ)Cijkl​εkl​(ξ)dξ

这是 Eringen 非局部理论的核心。不要被这些符号吓到。它所表达的只是:要找到我们关注的点 xxx 的应力,我们需要遍历整个物体。在每一站 ξ\xiξ,我们找到那里存在的局部应力(Cijklεkl(ξ)C_{ijkl} \varepsilon_{kl}(\xi)Cijkl​εkl​(ξ)),将其乘以一个权重因子 α\alphaα,然后加到我们的总和中。

这个权重因子 α\alphaα 就是​​衰减核函数​​,它是秘密武器。它告诉我们一个点 ξ\xiξ 对我们关注的点 xxx 有多大的影响。对于一个行为良好的材料,这种影响应该只取决于它们之间的距离 ∣x−ξ∣|x-\xi|∣x−ξ∣,以及一个新的、至关重要的材料属性:​​内禀长度尺度​​ ℓ\ellℓ。

建立在坚实基础上:保留什么,舍弃什么?

当一个新理论出现时,人们很自然会问:我们是否要抛弃所有学过的一切?对于非局部力学,答案是令人安心的“不”。该理论是一个谨慎而深刻的修正,而非全盘推翻。

力学的结构建立在三大支柱之上:

  1. ​​运动学​​:运动和变形的几何学(我们如何根据位移定义应变)。
  2. ​​平衡定律​​:游戏的基本、普适规则(质量、动量和能量守恒)。
  3. ​​本构律​​:特定材料的“个性”(特定材料如何响应力)。

非局部力学保留了前两大支柱。我们将应变作为位移梯度来测量的方式保持不变。线性动量平衡方程 ∇⋅σ+b=ρu¨\nabla \cdot \boldsymbol{\sigma} + \mathbf{b} = \rho \ddot{\mathbf{u}}∇⋅σ+b=ρu¨,也就是连续介质的牛顿第二定律,依然成立。质量守恒定律也未动摇。

所改变的——也正是其核心所在——是第三根支柱:​​本构律​​。我们只是用新的非局部规则“应力取决于各处的应变场”取代了旧的局部规则“应力取决于此处的应变”。物理学的基本定律的宏伟性得以保留;我们只是赋予了材料一种更复杂、也更真实的“个性”。

机器的灵魂:核函数及其长度

衰减核函数 α\alphaα 及其特征长度 ℓ\ellℓ 是非局部这台机器的灵魂。它们决定了非局部性的特征。核函数就像一个影响范围。离 xxx 很远的点,其 α\alphaα 值会很小,对 xxx 处的应力贡献甚微。而邻近点的 α\alphaα 值会很大,贡献也很大。

长度尺度 ℓ\ellℓ 设定了该影响范围的大小。具有较大 ℓ\ellℓ 值的材料“非局部性更强”——它对其远邻正在做什么有更长的记忆。而 ℓ\ellℓ 值很小的材料则几乎是局部的,只关心其直接邻域。事实上,为使理论保持一致,核函数的定义必须使得当 ℓ→0\ell \to 0ℓ→0 时,它会演变成著名的​​狄拉克δ函数​​——一个除了单一点外处处为零的数学尖峰。在此极限下,积分公式优美地退化回经典的局部理论,这正是它应有的表现。

为了感受这一点,让我们考虑一个具体的核函数,一种与物理学中的屏蔽势相关的形式:

α(r)=14πℓ2r exp⁡(−rℓ)\alpha(r) = \frac{1}{4\pi \ell^{2} r}\,\exp\left(-\frac{r}{\ell}\right)α(r)=4πℓ2r1​exp(−ℓr​)

指数项 exp⁡(−r/ℓ)\exp(-r/\ell)exp(−r/ℓ) 清楚地表明,当距离 rrr 大于 ℓ\ellℓ 时,影响迅速减小。如果你进行计算,你会发现使用该核函数的材料的平均相互作用距离是 2ℓ2\ell2ℓ,均方根距离是 ℓ6\ell\sqrt{6}ℓ6​。长度尺度 ℓ\ellℓ 不仅是一个抽象参数;它直接衡量了​​原子间作用力的范围​​。

核函数的具体数学形式也具有深远的影响。考虑两种流行的选择:形如钟形曲线的高斯核函数,和更为尖锐的双指数核函数。当我们研究这些核函数如何与不同长度的波相互作用时(通过傅里叶变换),我们发现它们的行为不同。它们充当低通滤波器,对极短波长变形的抑制作用大于长波长变形。但这些滤波器的“截止频率”因核函数形状而异。这告诉我们,原子间作用力随距离衰减的精细细节对材料的宏观动态行为有可测量的影响。

我们为什么需要它:驯服无穷大与聆听波

你可能会说,这一切都非常优雅,但它有什么实际用途吗?答案是响亮的“有”。非局部力学不仅仅是一个数学上的奇趣之物;它是一个强大的工具,解决了长期存在的难题,并描述了经典理论完全忽略的现象。

驯服无穷大

经典力学最引人注目的失败之一发生在断裂问题中。如果你使用线性弹性理论计算一个尖锐裂纹尖端的应力,你会得到一个令人不寒而栗的答案:无穷大。这当然在物理上是不可能的。没有材料是无限强的。原子在键断裂前只能被拉开一定的距离。这个无穷大表明我们的局部模型在裂纹尖端附近的小尺度上失效了。

非局部力学应运而生。本构律的积分特性“抹平”了应力。裂纹尖端的应力不再由尖端处的无限应变决定,而是由一个大小为 ℓ\ellℓ 的邻域内应变的加权平均值决定。这个平均过程巧妙地驯服了无穷大,得到了一个有限的、物理上合理的应力。现在的最大应力大约为 σmax⁡∼KI/2πℓ\sigma_{\max} \sim K_I / \sqrt{2\pi\ell}σmax​∼KI​/2πℓ​,其中 KIK_IKI​ 是经典理论的应力强度因子。非物理的奇异性消失了,取而代之的是由内禀长度尺度支配的新物理学。

尺寸很重要与波色散

在人类尺度上,一根1厘米宽的钢丝和一根1米宽的钢梁,其行为本质上是相同的。它们的材料属性是恒定的。但在纳米尺度上,情况不再如此。一根10纳米厚的金属丝可能比同种材料制成的100纳米厚的金属丝刚度更高。这种​​尺寸效应​​在经典力学中无法解释,因为它没有内在的长度尺度。而非局部力学,凭借其内置的长度 ℓ\ellℓ,自然地捕捉了这些效应。当一个物体的尺寸变得与 ℓ\ellℓ 相当时,其行为就开始改变。

这同样体现在波的传播方式上。在经典连续介质中,声速是一个常数,与波的频率或波长无关。它是非色散的。但在真实的原子晶格或非局部连续介质中,事实并非如此。高频(短波长)的波更敏锐地感受到材料离散、颗粒状的本质。非局部理论的一个关键预测是​​色散​​:波速随波长而变化。

对于一种常见的非局部模型,我们发现波的相速度 cpc_pcp​ 与经典声速 c0c_0c0​ 的关系如下:

cp(k)c0=11+(kℓ)2\frac{c_p(k)}{c_0} = \frac{1}{\sqrt{1+(k\ell)^2}}c0​cp​(k)​=1+(kℓ)2​1​

其中 kkk 是波数(与波长成反比)。对于长波(kkk 很小),速度就是 c0c_0c0​。但随着波长变短(kkk 变大),波速下降。材料对于这些短促而剧烈的振动显得更“软”。这种被称为​​软化​​的现象,正是在许多原子模拟和实验中观察到的情况。

宇宙中不孤单:一个非局部思想的家族

Eringen 的理论是一个基石,但它是一个更大的“广义”连续介质理论家族的一部分,所有这些理论都试图将长度尺度纳入力学中。其中两个最重要的亲缘理论是​​应变梯度弹性理论​​和​​近场动力学​​。

应变梯度理论走的是一条不同的哲学路径。它不采用积分,而是认为某一点的能量(以及应力)不仅取决于应变,还取决于应变的梯度——即应变在附近的變化速率。这也引入了一个长度尺度,但其数学结构有根本的不同。在波的世界(傅里叶空间)中,Eringen 的模型将经典刚度乘以一个类似 (1+k2ℓ2)−1(1+k^2\ell^2)^{-1}(1+k2ℓ2)−1 的项,而一个简单的应变梯度模型则乘以 (1+k2ℓg2)(1+k^2\ell_g^2)(1+k2ℓg2​)。一个是分数,另一个是多项式。除了在 ℓ\ellℓ 和 ℓg\ell_gℓg​ 都为零的平凡局部极限情况下,它们无法等同。它们是描述具有内部结构的世界的两种独特而强大的方式。

​​近场动力学​​甚至更具革命性。它抛弃了应力和应变梯度的概念。取而代之的是,它假定物体中的任意两点由一个“键”连接,如果它们的相对距离发生变化,这个键就会产生力。一个点上的总力是其视界 δ\deltaδ 内所有邻居的键力之和。如果说 Eringen 的模型是一个修正经典应力-应变关系的“改良派”,那么近场动力学就是一个从积分运动方程的新基础上重建力学的“革命派”。一个有趣的结果是,最简单的“键基”近场动力学模型实际上比经典弹性理论的通用性要差,在三维情况下其泊松比被固定为 1/41/41/4。这凸显了在构建物理模型时所涉及的微妙选择和后果。

美丽的陷阱:简化的诱惑与危险

非局部力学的积分公式虽然强大,但求解起来可能很困难。这催生了一种诱人的简化:微分形式。对于一种特殊的核函数(亥姆霍兹核函数),该积分定律等价于一个看似简单的微分方程:

(1−ℓ2∇2)σij=σijlocal(1 - \ell^2 \nabla^2) \sigma_{ij} = \sigma^{\text{local}}_{ij}(1−ℓ2∇2)σij​=σijlocal​

这在数值上处理起来容易得多。然而,这种简化是一个美丽的陷阱。虽然在无限介质中完全有效,但当轻率地应用于具有边界的有限体时,它可能导致悖论。

考虑一个在自由端受载的简单悬臂梁。如果你应用微分模型,奇怪的事情发生了。在梁的内部,经典弯矩是线性的,其二阶导数为零。项 ℓ2∇2σ\ell^2 \nabla^2 \sigmaℓ2∇2σ 消失了,非局部方程退化回局部方程!这个为非局部性而设计的模型,却预测出纯粹的局部行为,没有表现出任何尺寸效应。这是因为微分形式在边界处是对积分的糟糕近似,而所有有趣的非局部物理现象都发生在边界处。

此外,这种微分形式以一种会产生新的、非经典边界层效应的方式修改了控制方程。自平衡载荷的影响,在经典理论中其衰减距离与载荷宽度相当(圣维南原理),现在伴随着与 ℓ\ellℓ 直接相关的第二种更快的衰减模式。衰减指数不再仅仅是 kkk,而是 k2+1/ℓ2\sqrt{k^2 + 1/\ell^2}k2+1/ℓ2​。这不是一个错误,而是该模型的一个特征,告诉我们非局部材料比经典材料能更快地屏蔽掉边界扰动。

这是一个对任何科学家都深刻的教训。我们的模型很强大,但它们有其有效范围。从完整的积分描述到简化的微分形式的旅程,铺满了隐藏的假设。认识到这些微妙之处不是理论的失败,而是真正理解的标志,揭示了我们试图描述的物理学那深刻而复杂的美。

应用与跨学科联系

既然我们已经探索了非局部力学的基本原理,你可能会问一个完全合理的问题:“这一切都是为了什么?”我希望你会发现,答案既深刻又优美。从经典力学到非局部力学的旅程不仅仅是一次数理练习;它是解决我们旧理论中深层悖论,并开启对物理世界——尤其是最小尺度世界——全新、更精确理解的必要一步。让我们踏上一段旅程,看看这个新视角将我们带到哪些非凡之地,从材料的灾难性失效到纳米技术的精微振动。

驯服无穷大

任何新物理理论的伟大胜利之一就是它能够攻克其前辈理论的顽疾。对于经典连续介质力学来说,最可怕的顽疾之一是奇异性——即在单一点上预测出无限的力或应力。这不仅仅是美学上的瑕疵;它表明理论正在崩溃。毕竟,自然界不会产生无穷大。

考虑断裂问题。如果你观察一个含有尖锐裂纹的材料,经典线性弹性理论预测裂纹尖端的应力是无限的,其随距尖端距离 rrr 按 1/r1/\sqrt{r}1/r​ 的规律变化。但物理应力怎能是无限的呢?不可能。材料必然会屈服,或者键会断裂,但理论本身并没有提供一个有限而合理的答案。非局部力学应运而生。通过承认裂纹尖端的应力取决于其周围一个小邻域内材料的状态,该理论施展了一点奇妙的魔法。尖锐的、无限的奇异性被“抹平”或正则化到一个由内禀长度 ℓ\ellℓ 表征的区域上。应力不再飙升至无穷大;相反,它在尖端处上升到一个很高但有限的峰值。这个峰值应力由载荷(经典应力强度因子 KIK_IKI​)和材料自身的内禀长度决定,其量级与 KI/ℓK_I/\sqrt{\ell}KI​/ℓ​ 成正比。突然之间,断裂力学变得更加物理,非物理的无穷大被一个可测量的、有限的强度所取代。

这种对无穷大的驯服不仅限于裂纹。当我们在材料科学中研究位错时,也出现了同样的问题。位错是晶体原子点阵中的线状缺陷,是金属弯曲和变形方式的基础。经典理论预测螺位错的应力场在其核心处也有一个奇异性,按 1/r1/r1/r 变化。这同样是一个非物理的人为产物。而非局部力学再次漂亮地解决了它。通过将经典应力场与非局部核函数进行卷积,该理论用一个光滑、有限的应力分布取代了位错核心周围的尖锐奇异性。这使得对位错能量及其与晶体中其他特征相互作用的计算更加真实,这是现代材料物理学的基石。

纳米世界中波与振动的舞蹈

现在让我们从静态的悖论转向动态的波与振动世界。想象一根吉他弦。它的音高或频率取决于其长度、张力和质量。经典理论为任何给定的弦提供了一组谐波频率。但如果这根“弦”是一个纳米梁,一块只有几百个原子厚度的微小硅片,被设计成纳米传感器或计算机的组件呢?在这个尺度上,经典图像还成立吗?

非局部力学断然告诉我们,它不成立。它预测了一种被称为“软化”的现象。非局部相互作用就像材料内部额外的长程连接,使其在弯曲时更具柔顺性或更“软”。因此,一个非局部纳米梁的振动频率比其经典对应物更低。这种效应被一个极其简洁优美的公式所捕捉,该公式给出了非局部频率 ωn\omega_nωn​ 与经典频率 ωn(0)\omega_{n}^{(0)}ωn(0)​ 之比:

ωnωn(0)=11+(nπ)2(ℓL)2\frac{\omega_{n}}{\omega_{n}^{(0)}} = \frac{1}{\sqrt{1 + (n\pi)^{2} \left(\frac{\ell}{L}\right)^{2}}}ωn(0)​ωn​​=1+(nπ)2(Lℓ​)2​1​

其中 LLL 是梁的长度, nnn 是模态数。你可以立刻看到,因为分母总是大于1,所以非局部频率总是更低。此外,当内禀长度 ℓ\ellℓ 是梁长 LLL 的一个显著部分时,这种效应最为明显,而这恰恰是纳米技术中的情况。

这不仅仅是一个抽象的公式;它对纳米机电系统 (NEMS) 有着深远的影响。这些微型谐振器是超灵敏质量探测器和高频信号处理器的构建模块。正确预测它们的共振频率对其设计至关重要,而非局部理论提供了必要的修正。

这种现象不仅限于一维梁。考虑石墨烯,一种非凡的单原子厚度的碳片。当我们研究声波(声学声子)通过这个二维薄膜的传播时,非局部理论预测了一种有趣的行为。在经典世界里,声速是恒定的,与其频率无关。在非局部世界里,这不再成立——材料表现出色散性。对于短波长的波,当波长与内禀长度 ℓ\ellℓ 相当时,相速度和群速度会减小。频率不再随波数 kkk 线性增加;相反,它会接近一个最大的饱和值。这正是在真实原子晶格中发生的情况,其中振动的最大频率由原子间距和键强度决定。非局部力学,一个连续介质理论,却习得了离散原子世界的智慧。

弯曲、拉伸与边界的微妙艺术

非局部性的后果并不总是像驯服无穷大或改变波速那样引人注目。有时,它们在结构的静态响应中以更微妙但同样重要的方式表现出来,尤其是在边界附近。一个关键主题是应力的重新分布。

想象一下拉伸一根简单的纳米棒。经典的均匀弹性理论会告诉你应力处处恒定。非局部理论提供了一幅更细致的图景。虽然在材料内部深处(“体相”)应力确实几乎是恒定的,但它会形成独特的“边界层”,其中应力分布迅速变化以满足自由端的条件。这是材料解释表面附近原子比中间原子拥有更少相互作用邻居的方式。控制这种现象的方程通常是亥姆霍兹型方程,它自然会产生这些指数衰减的边界层解。

然而,我们必须小心。尺寸依赖力学的世界丰富而复杂,不同的模型有时可能导致不同——甚至看似矛盾——的预测。例如,Eringen 非局部模型最简单的*微分*形式虽然有用,但却表现出某些“悖论”。对于一个承受均布体力的杆 或一个承受端部载荷的悬臂梁,这个特定模型预测的响应与经典响应完全相同!非局部效应似乎消失了。这是一个数学上的人为产物,因为模型的算子 (1−ℓ2∇2)(1 - \ell^2 \nabla^2)(1−ℓ2∇2) 作用于一个在这些情况下二阶导数为零的应力场,从而抵消了非局部修正项。这是一个重要的教训:我们的数学模型是近似的,我们必须理解它们的局限性。

这种复杂性为一场引人入胜的科学辩论打开了大门。虽然非局部弹性的积分形式几乎总是预测软化行为(结构变得更柔顺),但也有其他理论被提出来捕捉不同的尺寸效应。例如,惩罚应变急剧变化的*应变梯度弹性理论,预测了硬化*效应——纳米梁变得比其经典对应物更刚硬。哪个是正确的?答案取决于材料和在纳米尺度上占主导地位的物理机制。自然界比任何单一模型都更为复杂,而这种各种相互竞争的理论之间活跃的相互作用推动着该领域向前发展。

连接理论、实验与计算

此时,你可能想知道这一切是否只是一个美丽的理论建构,一座空中楼阁。我们如何将它与现实世界联系起来?这正是科学方法大放异彩之处,它将理论、实验和计算结合在一起。

一个理论只有在其可被检验的能力范围内才强大。非局部理论引入了一个新的材料参数,即内禀长度 ℓ\ellℓ。我们能测量它吗?答案是肯定的。考虑一个实验,我们对一系列不同长度 LLL 的微梁进行弯曲测试。非局部理论预测了中点挠度、载荷、梁长和内禀长度 ℓ\ellℓ 之间的特定关系。通过巧妙地绘制实验数据——例如,绘制挠度除以长度与长度平方的关系图——理论公式就变成了一条直线。从这条直线的斜率和截距,我们可以通过实验确定 ℓ\ellℓ 的值。曾经纯粹的理论参数变成了一个可测量的材料属性,就像其密度或杨氏模量一样。

最后,这些思想融入了工程师用来设计从飞机到微芯片等各种产品的强大计算工具中。有限元法 (FEM) 是现代工程的基石,它将复杂的结构分解为更简单的“单元”网格。为了创建一个“非局部有限元”,我们可以修改每个微小单元的刚度。通过做一个简单的近似,即每个单元内的变形由单一特征波长主导,我们可以推导出一个包含非局部软化效应的“有效刚度”。这使我们能够将非局部性的深层物理学直接构建到预测下一代纳米尺度设备行为的软件中。

归根结底,进入非局部力学的旅程是物理学故事中引人入胜的一章。它向我们展示了如何修正珍贵旧理论的不足,揭示了一个更丰富、更准确的现实图景。它是一个统一的概念,触及断裂力学、材料科学、动力学和纳米技术,为描述那个尺寸真正重要的奇妙世界提供了一种通用语言。