try ai
科普
编辑
分享
反馈
  • 非局域模型

非局域模型

SciencePedia玻尔百科
核心要点
  • 当材料的微观结构长度与物理场变化的尺度相当时,经典的局域理论会失效,导致断裂力学中的悖论,并且无法解释尺寸效应。
  • 非局域模型通过引入一个内禀长度尺度来解决这些问题,将材料点的状态定义为其邻域内状态的加权平均。
  • 虽然基于积分的非局域模型在物理上是基础性的,但其在数学上更简单的微分(梯度)对应形式可能会引入新的悖论,并需要谨慎处理边界条件。
  • 非局域框架提供了一种统一的方法,能够精确模拟多种现象,包括材料失效、纳米尺度器件的力学行为以及化学中的溶剂屏蔽效应。

引言

在经典物理学和工程学中,我们对物质世界的理解建立在局域性原理之上:某一点的性质仅取决于该精确点上发生的情况。这个被称为连续介质假设的强大思想,使我们能够用优美的微分方程描述从桥梁到行星的一切事物。然而,这种观点是一种近似,当我们研究极小尺度或极端条件下的材料时,它会显著失效。当材料失效集中于尖锐裂纹,或当构件缩小至纳米尺度时,局域理论会预测出物理上的荒谬结果,揭示了我们经典模型中的一个根本性缺陷。

本文通过探索非局域模型的世界来应对这一挑战,这是一个理论框架,其中一个点的行为受其整个邻域的影响。首先,在“原理与机制”一章中,我们将剖析非局域性的基本思想,对比其积分和微分形式,并探讨其中出现的微妙之处和悖论。随后,“应用与跨学科联系”一章将展示这些模型的预测能力,说明它们如何解决断裂力学中长期存在的难题,解释纳米材料的独特性质,甚至在材料科学与理论化学之间架起桥梁。

原理与机制

想象一下,你正试图描述一个房间里的温度。在大多数情况下,你可以在房间里任选一点,在那里放一个温度计,就能得到那一点的温度。那个确切位置的空气分子决定了读数。这就是​​局域​​理论的精髓。一个点上的性质——应力、温度、电场——仅取决于那个确切点上发生的情况。这个思想,被称为​​连续介质假设​​,是经典物理学和工程学的基石。这是一个极其强大的近似,使我们能够使用优美的微积分工具来描述世界。但从根本上说,它是一个近似。和所有近似一样,它也有其局限性。

当世界不再局域

这种舒适的局域图景何时开始崩塌?当我们“放大”得太近时,无论是观察非常微小的事物,还是观察在空间中变化非常迅速的现象,这种情况就会发生。连续介质假设依赖于一个称为​​尺度分离​​的原则。它假设材料微观结构的特征长度,我们称之为 lml_mlm​(可以想象成金属中的晶粒尺寸,或复合材料中纤维的间距),远小于物理场(如应变)变化的长度尺度,我们称之为 LgL_gLg​。用数学符号简写,局域世界在 lm≪Lgl_m \ll L_glm​≪Lg​ 时有效。

然而,自然界频繁地违反这个条件。

考虑一块被拉伸直至开始失效的金属。曾经均匀分布的变形会突然集中到一个称为​​剪切带​​的非常狭窄的区域。在这个带内,应变在与金属晶粒尺寸相当的距离上急剧变化。在这里,LgL_gLg​(剪切带的宽度)与 lml_mlm​ 处于同一数量级!。试图描述这种情况的经典局域理论会遇到一个灾难性的问题。它预测剪切带应该是无限薄的,这会导致物理上荒谬的结果,比如破坏材料所需的能量为零。这不仅仅是一个数学上的奇特现象;在使用局域模型的计算机模拟中,预测的失效区域会随着计算网格的尺寸而缩小,给出的答案取决于模拟的设置,而不是材料的物理性质。

或者想一想微机电系统(MEMS)中的一个微小构件。如果该构件本身只有几十个材料晶粒的宽度,那么其整体尺寸 HHH 就成为特征长度 LgL_gLg​。同样,lm≪Lgl_m \ll L_glm​≪Lg​ 的条件不再成立。实验证实,在这种情况下,越小可能越强;材料的测量强度开始取决于被测试样的尺寸,这是一种局域理论根本无法解释的​​尺寸效应​​。

即使在大型物体中,快速振动也能揭示局域性的失效。当在材料中传播的波的波长 λ\lambdaλ 与微观结构长度 lml_mlm​ 相当时,波开始“感觉”到材料的离散性质。这导致了​​色散​​现象,即波速取决于其频率,这是纯局域模型所忽略的另一个现象。在所有这些情况中——剧烈梯度、小尺寸和短波长——一个点只关心自身的舒适假设不再成立。我们需要一个新的思想。

材料点的民主

这个新思想非常直观:一个点上发生的事情不仅由该点的状态决定,还由其整个邻域内状态的加权平均值决定。这就是​​非局域性​​的原理。一个材料点在决定感受多大应力之前,实际上会“征求其邻居的意见”。

表达这一思想最直接的方式是使用积分。在经典局域弹性理论中,点 x\mathbf{x}x 处的应力 σ\boldsymbol{\sigma}σ 与同一点的应变 ε\boldsymbol{\varepsilon}ε 成正比:σ(x)=C:ε(x)\boldsymbol{\sigma}(\mathbf{x}) = \mathbf{C}:\boldsymbol{\varepsilon}(\mathbf{x})σ(x)=C:ε(x),其中 C\mathbf{C}C 是刚度张量。这个思想最简单的非局域推广是,点 x\mathbf{x}x 处的应力是物体中所有其他点 x′\mathbf{x}'x′ 处局域弹性响应的空间平均值:

σ(x)=∫Ωα(∣x−x′∣;ℓ) (C:ε(x′)) dVx′\boldsymbol{\sigma}(\mathbf{x}) = \int_{\Omega} \alpha(|\mathbf{x}-\mathbf{x}'|; \ell) \, (\mathbf{C}:\boldsymbol{\varepsilon}(\mathbf{x}')) \, \mathrm{d}V_{\mathbf{x}'}σ(x)=∫Ω​α(∣x−x′∣;ℓ)(C:ε(x′))dVx′​

其中的奥妙在于函数 α(r;ℓ)\alpha(r; \ell)α(r;ℓ),它被称为​​影响函数​​或​​核函数​​。它充当一个加权因子。离 x\mathbf{x}x 很近的点 x′\mathbf{x}'x′ 影响很大,而远处的点影响则减弱。该函数引入了一个新的基本材料参数 ℓ\ellℓ,即​​内禀长度尺度​​。这个长度尺度决定了材料点之间“对话的范围”,并且在物理上与我们之前讨论的微观结构长度 lml_mlm​ 相关联。

为了使这个思想保持一致,核函数必须满足一个简单的常识性规则:如果应变处处均匀,即 ε(x′)=ε0\boldsymbol{\varepsilon}(\mathbf{x}') = \boldsymbol{\varepsilon}_0ε(x′)=ε0​,非局域模型必须能退化为经典结果。这意味着非局域应力也应该是均匀的,且等于 C:ε0\mathbf{C}:\boldsymbol{\varepsilon}_0C:ε0​。这只有在影响函数在整个空间上的积分等于1时才能成立:

∫R3α(r;ℓ) dr=1\int_{\mathbb{R}^3} \alpha(r; \ell) \, \mathrm{d}\mathbf{r} = 1∫R3​α(r;ℓ)dr=1

这种归一化确保了“民主投票”过程不会人为地夸大或削弱共识。如果邻域中的每个人都同意,那么平均意见就是那个意见。在内禀长度 ℓ\ellℓ 趋于零的极限情况下,一个恰当定义的核函数会演变成狄拉克δ函数 δ(x−x′)\delta(\mathbf{x}-\mathbf{x}')δ(x−x′),该函数除了在 x′=x\mathbf{x}' = \mathbf{x}x′=x 处外处处为零。在此极限下,积分坍缩,我们恢复了纯粹的局域定律。因此,经典理论被巧妙地作为一种特殊情况包含在非局域框架之内。

从积分到微分,再返回

虽然积分形式是非局域性最基本、物理上最直观的描述,但处理积分可能非常困难。如果我们能用一个更简单的微分方程来捕捉同样的物理现象,那就太好了。

令人惊讶的是,对于一种非常特殊的核函数选择,这是可能的!如果我们在一个无限空间中工作(以避免复杂的边界效应),并选择一个特定的指数衰减核函数,即所谓的​​Yukawa核​​或​​Helmholtz核​​,复杂的积分定律就变得与一个看起来简单得多的微分方程完全等价。在一维情况下,这个神奇的核函数是 α(x)=12ℓexp⁡(−∣x∣/ℓ)\alpha(x) = \frac{1}{2\ell} \exp(-|x|/\ell)α(x)=2ℓ1​exp(−∣x∣/ℓ)。使用这个核函数,积分定律

σ(x)=∫−∞∞12ℓexp⁡(−∣x−ξ∣/ℓ) Eε(ξ) dξ\sigma(x) = \int_{-\infty}^{\infty} \frac{1}{2\ell} \exp(-|x-\xi|/\ell) \, E \varepsilon(\xi) \, d\xiσ(x)=∫−∞∞​2ℓ1​exp(−∣x−ξ∣/ℓ)Eε(ξ)dξ

在数学上与以下微分方程完全相同:

σ(x)−ℓ2d2σdx2=Eε(x)\sigma(x) - \ell^2 \frac{d^2\sigma}{dx^2} = E \varepsilon(x)σ(x)−ℓ2dx2d2σ​=Eε(x)

这是一个优美的结果,是连接积分卷积世界和微分算子世界的一座桥梁。这种微分形式通常被称为​​梯度模型​​,因为如果你将其重新整理以求解应变,会发现应变取决于应力及其二阶导数(梯度的梯度)。其一般形式常写作 (I−ℓ2∇2)σ=C:ε(\mathbb{I} - \ell^2 \nabla^2) \boldsymbol{\sigma} = \mathbf{C} : \boldsymbol{\varepsilon}(I−ℓ2∇2)σ=C:ε。

这种等价性提供了一个强大的捷径。然而,这个魔法附带了细则。这种等价性仅对这种特定的核函数和仅在无限域中是精确的。使用不同的核函数(如高斯核),或者更重要的是,处理一个有边界的现实物体时,会破坏这种精确的等价性。微分形式于是变成了一种近似,必须非常小心地处理。

边缘与尖角的麻烦

微分形式的便利性带来了高昂的代价,揭示了从局域物理到非局域物理过渡中的深层微妙之处。问题,正如物理学中常见的那样,出在边界上。

微分形式 (I−ℓ2∇2)σ=C:ε(\mathbb{I} - \ell^2 \nabla^2) \boldsymbol{\sigma} = \mathbf{C} : \boldsymbol{\varepsilon}(I−ℓ2∇2)σ=C:ε 涉及到比经典力学中任何方程都更高阶的应力场导数。在数学上,一个更高阶的微分方程需要更多的边界条件才能获得唯一解。但这些额外的边界条件从何而来?经典力学只告诉我们如何在边界上指定力或位移,而不是“超曳力”或其他高阶量。相比之下,积分模型不需要这样的额外条件;经典的边界条件就足够了。这种差异表明,微分形式不仅仅是积分的简单替代品;它是一种不同的数学实体,提出了新的物理和数学问题。

其后果可能是戏剧性的,导致了所谓的​​悬臂梁悖论​​。考虑一个教科书中的问题:一个简单的悬臂梁,一端固定,另一端受到一个点力 PPP 的作用。如果你应用微分形式的非局域模型,会发生一件奇怪的事。梁的平衡方程要求,在没有分布载荷的区域,弯矩的二阶导数必须为零,即 M′′(x)=0M''(x)=0M′′(x)=0。当你把这个代入非局域本构律 M(x)−ℓ2M′′(x)=EIκ(x)M(x) - \ell^2 M''(x) = EI \kappa(x)M(x)−ℓ2M′′(x)=EIκ(x) 时,非局域项 ℓ2M′′(x)\ell^2 M''(x)ℓ2M′′(x) 完全消失了!该模型矛盾地预测出一个与经典局域理论完全相同的响应,完全没有显示出任何尺寸效应。你试图引入的非局域性就这样凭空消失了。这个警示故事教给我们一个深刻的教训:物理的平衡定律和材料的本构定律是紧密交织在一起的,如果这种相互作用没有得到充分的尊重,一个看似复杂的模型也可能以意想不到的方式失败。

非局域思想家族

简单微分模型带来的挑战激励了物理学家和工程师们发展出一个丰富的非局域理论家族,每种理论都有其自身的优点和微妙之处。

  • ​​应变驱动与应力驱动模型:​​ 我们讨论的原始公式,其中应变决定了非局域应力,被称为​​应变驱动​​模型。有些悖论可以通过字面上颠倒本构逻辑来解决。在​​应力驱动​​模型中,是通过对应力场进行平均来确定局域应变:ε(x)=∫α(∣x−x′∣) (C−1:σ(x′)) dVx′\varepsilon(\mathbf{x}) = \int \alpha(|\mathbf{x}-\mathbf{x}'|) \, (\mathbf{C}^{-1}:\boldsymbol{\sigma}(\mathbf{x}')) \, \mathrm{d}V_{\mathbf{x}'}ε(x)=∫α(∣x−x′∣)(C−1:σ(x′))dVx′​。这看似只是一个微小的代数重排,但它导致了一个具有不同且有时更理想的数学性质的模型。

  • ​​强非局域与弱非局域模型:​​ 这将我们引向一个有用的分类。基于显式积分卷积的模型通常被称为​​强非局域​​模型。它们以最纯粹的形式体现了“超距作用”原理。相比之下,通过将应变(或其他内变量)的梯度直接纳入材料能量函数来构建非局域性的模型被称为​​弱非局域​​或​​应变梯度​​理论。例如,能量密度可能同时依赖于 ε\boldsymbol{\varepsilon}ε 和 ∇ε\nabla\boldsymbol{\varepsilon}∇ε。这些模型也引入了内禀长度尺度 ℓ\ellℓ 并且能够捕捉尺寸效应,但它们没有显式的积分。它们付出的数学代价是要求位移场“更光滑”,因为能量现在依赖于位移的二阶导数,而积分模型对光滑性的要求则不那么严格。

从简单的局域世界到相互关联的非局域世界的这一旅程展示了科学的实践过程。一个优美而强大的理论(经典连续介质理论)在面对新实验时遇到了它的极限。一个新的、直观的思想(非局域性)被提出,引发了一系列不同的数学表述,每一种都有其自身的美感、力量和意想不到的陷阱。探索一个材料点如何与其邻居交流的征途仍在继续,推动着物理学和工程学的前沿。

应用与跨学科联系

空间中的一个点“知道”什么?如果你信奉经典物理学那优美而简单的世界观,一个点只知道其自身无限小位置上的条件——场、力、应变。这就是局域性原理。这是一个极其强大且成功的思想,是大部分经典力学和场论赖以建立的基石。这就像通过一个针孔观察世界;你可以构建一幅完整的画面,但只能通过拼合无数独立的、局域化的快照来实现。

但当这种近似失效时会发生什么?如果一个点对其周围环境并非一无所知呢?如果它能感受到邻居的影响,不仅是无限近的邻居,而是跨越一个微小但有限的距离的邻居呢?这就是非局域模型的世界。进入这个世界不仅仅是一次数学练习;这是我们必须踏上的旅程,以解决悖论、解释新实验,并统一我们对横跨众多学科的现象的理解。我们发现,这个简单的问题,“一个点知道什么?”,将我们引向对自然更深刻、更强大的描述。

当物体断裂时:失效物理学

让我们从一些戏剧性的事情开始:物体断裂的方式。想象一下拉伸一根金属棒直到它开始失效。一个经典的、局域的材料模型假设棒中的每个点仅根据该精确点的应变来决定是否软化和屈服。如果你在计算机模拟中使用这个模型,你会遇到一场灾难。当材料开始失效时,模拟会预测所有的变形都将集中到一个厚度为零的裂纹中。这会导致诸如无限应变以及预测的断裂能不取决于材料本身而取决于你选择的计算网格细度的荒谬结果!该模型已经失去了与物理现实的联系。

这种“病态”行为的根源在于局域假设。一个真实的裂纹不是一条宽度为零的数学线;它是一个过程区,在这个微小但有限的体积内,化学键被拉伸和断裂。非局域模型通过将这一物理事实直接构建到数学中,将我们从悖论中解救出来。它们引入了一个*内禀材料长度尺度*,我们称之为 ℓ\ellℓ。在这些模型中,一个点的应力或损伤不是由局域应变决定的,而是由一个大小为 ℓ\ellℓ 的邻域内应变的加权平均值决定的。

这一个改变带来了深远的影响。数学问题再次变得适定,裂纹现在有了一个与 ℓ\ellℓ 相关的自然的、有限的宽度。计算出的破坏材料所需的能量现在收敛到一个真实的物理值。这个长度尺度 ℓ\ellℓ 不仅仅是一个凑数因子;它是一个真实的材料属性,可以通过诸如 ℓ∝EΓ/σth2\ell \propto E\Gamma/\sigma_{\mathrm{th}}^{2}ℓ∝EΓ/σth2​ 的关系与材料刚度 EEE、表面能 Γ\GammaΓ 和理论强度 σth\sigma_{\mathrm{th}}σth​ 等基本量联系起来。

在研究薄的、有延展性的金属板断裂时,我们可以看到这一原理的实际应用。实验表明,测得的板材韧性可能取决于其厚度,如果韧性被认为是一种内禀材料属性,这就令人费解了。一个非局域模型完美地解决了这个难题。它认识到存在两个相互竞争的长度尺度:板材的几何厚度 ttt,以及断裂过程的内禀材料长度尺度 ℓ\ellℓ。观察到的行为取决于它们的比值 t/ℓt/\ellt/ℓ。通过使用非局域框架,我们可以理清这些效应,并提取出材料真实的、与尺寸无关的断裂属性,这是纯粹的局域理论无法完成的任务。

纳米尺度的世界:尺寸决定一切

当我们将整个世界缩小到纳米尺度时,对内禀长度尺度的需求变得更加明显。考虑一根微小的悬臂梁,比人类头发细数千倍,在真空中像微型跳水板一样振动。这类器件是纳米机电系统(NEMS)的核心。如果我们使用经典的局域弹性理论来预测其共振频率,我们会得到一个特定的数值。但是当我们进行实验时,我们通常会测量到一个更低的频率。这根梁似乎比我们的经典理论预测的更“软”或更柔韧。

再一次,非局域性是关键。在真实的晶体中,原子由化学键连接。一个点的应力不仅仅是该处晶格变形的函数,还与几个原子之外的晶格变形有关。一个非局域弹性模型,比如由 Eringen 提出的模型,通过将应力定义为应变场在一个小邻域上的积分来捕捉这一点。原子间的这种“超距作用”提供了一个局域理论中所没有的额外柔度机制。结果如何?非局域模型预测了梁的软化和更低的共振频率,与实验观察结果精确匹配。在局域世界观中看似异常的现象,在非局域世界观中变成了一个自然的预测。

连接世界:从分子到材料

非局域思想是连接不同物理尺度之间鸿沟的强大工具。我们如何描述一个体材料的性质,明知它是由复杂、异质的微观结构组成的?想象一种由硬质陶瓷纤维嵌入软质聚合物基体中的复合材料。远离任何边缘,我们可以使用局域的“均匀化”理论来找到复合材料的有效平均刚度。但是在表面处会发生什么呢?纤维和基体的重复模式被突然切断。

经典的局域模型对这种截断是视而不见的。而非局域模型则内在地捕捉到了这一点。在像近场动力学(Peridynamics)这样的模型中,每个点都与其有限“作用范围” δ\deltaδ 内的邻居相互作用,靠近表面的点与之相互作用的邻居较少。其“相互作用邻域”是不完整的。这自然而然地自动改变了表面附近薄边界层内的有效刚度,导致了局域均匀化完全忽略的尺寸依赖效应。

同样的空间依赖响应思想出现在一个完全不同的领域:理论化学。考虑一个溶解在水中的离子。最简单的“隐式溶剂”模型将水视为一种均匀的介电胶体,由一个单一的数字——相对介电常数 ϵs≈80\epsilon_s \approx 80ϵs​≈80 来表征。这是一个局域模型。但水是由分子组成的,其屏蔽电场的能力取决于与离子的距离。在离子附近,其电场变化迅速,水分子无法完美地取向以屏蔽电荷。在远处,场弱且变化缓慢,它们则可以。

一个非局域介电模型通过使介电常数成为波矢的函数 ϵ(k)\epsilon(\mathbf{k})ϵ(k) 来完美地捕捉这一点。波矢 k\mathbf{k}k 与波长成反比;大的 ∣k∣|\mathbf{k}|∣k∣ 对应于快速变化的场(短距离),而小的 ∣k∣|\mathbf{k}|∣k∣ 对应于缓慢变化的场(长距离)。非局域模型正确地指出,当 ∣k∣→∞|\mathbf{k}| \to \infty∣k∣→∞(非常靠近离子)时,屏蔽很弱,ϵ(k)→1\epsilon(\mathbf{k}) \to 1ϵ(k)→1。当 ∣k∣→0|\mathbf{k}| \to 0∣k∣→0(远离离子)时,屏蔽很强,ϵ(k)→ϵs\epsilon(\mathbf{k}) \to \epsilon_sϵ(k)→ϵs​。这种被称为“介电平滑”的现象,导致了离子周围极化电荷的更真实、有限的分布,并为计算溶剂化自由能提供了更精确的计算,这是计算生物化学的基石。

一个统一的思想:运动和计算中的非局域性

非局域框架的力量甚至延伸到流体力学的混沌世界。湍流,即流体中旋转、不可预测的运动,是非局域现象的缩影。一个位置的涡流受到其周围流动的整个历史和结构的影响。Prandtl 著名的“混合长度”假设是一个杰出的局域近似,用于模拟驱动湍流的雷诺剪应力。但它仍然是一个近似。

我们可以通过定义一个非局域混合长度来构建一个更复杂的模型,其中一个点的有效混合是周围区域的加权平均值。这种积分形式提供了一种系统性的方法,可以将更多的物理学内容纳入模型。例如,我们可以通过仔细选择平均核函数的形式,来强制模型与著名的“壁面律”——在固体边界附近观察到的普适速度剖面——保持一致。这展示了非局域性的数学语言如何成为构建更好物理理论的强大工具箱。

最后,在看到了这些模型的预测能力之后,我们必须问一个实际问题:我们如何用它们进行计算?一个将每个点连接到其邻域内其他点的积分听起来像是一场计算噩梦,尤其是在拥有数千个处理器的现代超级计算机上。然而,在这里,非局域性的物理原理也引导我们找到了一个优雅的解决方案。

在大多数物理应用中,非局域相互作用虽然是有限的,但其作用范围并非无限。平均核函数具有“紧支集”或有限的作用范围 δ\deltaδ。一个点只关心其作用范围内的邻居。为了并行化计算,我们可以使用一种称为“区域分解”与“晕圈交换”的策略。计算机将问题划分为子域,每个处理器负责一个。为了计算自己域内的值,一个处理器只需要知道其邻居域中一个薄的“晕圈”区域内的系统状态。在每个计算步骤之前,处理器只需与其直接邻居交换这个晕圈数据。通信保持局域性,问题变得易于处理。短程非局域性的物理性质决定了最优的并行计算策略。

从钢中裂纹扩展的灾难性优雅,到离子周围水分子的微妙舞蹈,再到湍急河流的旋转混沌,非局域性原理提供了一条统一的线索。它提醒我们,有时,要真正理解一个点,我们必须欣赏它所处的邻域。这种视角的转变不仅仅是对我们旧理论的修正;它是一扇通往对物理世界更丰富、更准确、更相互关联的理解的大门。