try ai
科普
编辑
分享
反馈
  • 应变等效原理

应变等效原理

SciencePedia玻尔百科
核心要点
  • 应变等效原理指出,损伤材料的应变响应与一个假设的、承受更高“有效应力”的无损材料的应变响应相同。
  • 该原理通过假设材料的有效刚度由单一标量损伤变量D控制并均匀退化,从而简化了损伤建模。
  • 该框架具有热力学一致性,定义了作为材料退化和失效驱动力的损伤能量释放率。
  • 它提供了一种强大的方法来模拟损伤-塑性耦合现象,其核心假设是塑性流动由作用在完好材料上的有效应力决定。
  • 尽管局部版本的模型很简洁,但它存在局限性,例如模拟中的网格依赖性,这促使了更先进的非局部损伤模型的发展。

引言

预测材料何时会断裂是工程学和材料科学中最关键的挑战之一。虽然我们看到的材料是坚固且连续的,但在微观层面,它们是复杂的系统,其中损伤在最终失效前早已以空洞和裂纹的形式累积。核心问题在于,如何在不追踪每一个微小瑕疵的情况下,用数学方法描述这种渐进退化材料的行为。应变等效原理是连续介质损伤力学中的一个基本概念,它通过假设一个受损体与一个无损体之间的简单类比,为这一挑战提供了简洁而有力的解决方案。在本文中,我们将对这一原理进行全面的探索。第一章 ​​原理与机制​​ 将剖析该原理的理论基础,从有效应力的概念到其与热力学的深刻联系及其固有的局限性。随后,​​应用与跨学科联系​​ 章节将展示该原理的实际威力,展示其如何用于预测材料软化、诊断失效、模拟塑性等复杂行为,并构成先进计算模拟的基础。

原理与机制

想象一下你正握着一块金属。它感觉坚固、结实、连续。但如果你放大观察,深入到微观层面,你会看到一个由晶粒、晶体和微小缺陷构成的景象。现在,想象你开始拉伸这块金属。随着它的伸长,这些缺陷会增长。微观空洞出现,微小裂纹开始形成并连接起来。材料正在变得损伤。我们到底该如何描述这种复杂的、不断演变的“瑞士奶酪”材料的行为?我们是否必须追踪每一条裂纹和每一个空洞?这个任务似乎毫无希望。

但是,就像物理学中经常发生的那样,一个简单而有力的思想让我们能够见树又见林。这个思想就是我们所称的​​连续介质损伤力学​​的基础。

一个简单的图景:“瑞士奶酪”材料

让我们从一个简单的思想实验开始。考虑一根直杆在拉伸作用下。在宏观尺度上,它的横截面积为 A0A_0A0​。但在微观尺度上,它布满了空洞。实际仍然承载载荷的“有效”面积更小,我们称之为 AAA。

我们可以定义一个单一的数字来捕捉这种退化的本质:​​损伤变量​​,我们称之为 DDD。它就是损失面积的比例。

D=A0−AA0D = \frac{A_0 - A}{A_0}D=A0​A0​−A​

这个定义非常简洁。如果材料处于其原始、完好状态,没有面积损失,因此 A=A0A=A_0A=A0​ 且 D=0D=0D=0。如果材料即将断成两半,有效面积 AAA 趋近于零,而 DDD 趋近于 1。这个介于 0 和 1 之间的单一标量数字 DDD 告诉我们材料内部的健康状态。

两种应力的故事:我们所见的 vs. 材料所感的

这个简化的面积图景带来了一个关键的洞见。当我们对杆施加力 FFF 时,我们工程师和物理学家通常将应力计算为力除以初始面积。这是​​名义应力​​(或柯西应力),σ=F/A0\sigma = F/A_0σ=F/A0​。它是我们的仪器所测量的;它是对整个包含损伤和未损伤部分的横截面进行平均的应力。

但材料感受到的是什么呢?力 FFF 并非由空洞承载;它集中在材料的幸存部分上。这部分幸存区域上的应力,即​​有效应力​​ σ~\tilde{\sigma}σ~,必然更高。它是相同的力作用在更小的有效面积上,σ~=F/A\tilde{\sigma} = F/Aσ~=F/A。

现在,我们可以利用损伤变量将这两种应力联系起来。根据 DDD 的定义,有效面积为 A=A0(1−D)A = A_0(1-D)A=A0​(1−D)。将其代入可得:

σ~=FA0(1−D)=11−D(FA0)\tilde{\sigma} = \frac{F}{A_0(1-D)} = \frac{1}{1-D} \left(\frac{F}{A_0}\right)σ~=A0​(1−D)F​=1−D1​(A0​F​)

由于 σ=F/A0\sigma = F/A_0σ=F/A0​,我们得出了一个基石关系:

σ~=σ1−D\tilde{\sigma} = \frac{\sigma}{1-D}σ~=1−Dσ​

这个小小的方程式蕴含着丰富的物理意义。它告诉我们,材料完好部分感受到的应力总是大于我们测量的名义应力。随着损伤 DDD 的增加,即使我们保持施加的力恒定,有效应力 σ~\tilde{\sigma}σ~ 也会急剧上升。这就是为什么受损的东西会断裂!幸存的部分承受着巨大且不断增加的应力,直到它们也屈服。

等效原理:一个深刻的类比

所以我们有了“有效”应力的概念。它有什么用呢?毕竟,它只是一个数学构造。这时,一个天才的构想出现了,一个被称为​​应变等效原理​​的假设。

该原理指出:受损材料在真实应力 σ\boldsymbol{\sigma}σ 下的应变,与其未受损的对应物在有效应力 σ~\tilde{\boldsymbol{\sigma}}σ~ 作用下所经历的应变相同。

让我们来解读一下。这是一个类比。它说我们可以通过思考一个更简单的、想象中的物体来理解我们复杂的、受损的“瑞士奶酪”材料的行为:一个“有效的无损构型”。这个想象中的物体只是一块原始材料的完好块体。该原理宣称,如果我们用有效应力 σ~\tilde{\boldsymbol{\sigma}}σ~ 来拉伸这个想象中的完好块体,它的伸长量将与我们用名义应力 σ\boldsymbol{\sigma}σ 拉伸真实受损块体的伸长量完全相同。

这是一个深刻的简化。我们不需要为受损材料建立一套新的、复杂的理论,而可以重用旧的、简单的无损材料理论,但我们必须给它输入“正确”的应力——即材料实际感受到的有效应力。

损伤胡克定律:一个弱化世界的新法则

让我们看看这个优美思想的推论。对于一个健康的线性弹性材料,应力与应变的关系是胡克定律。在其通用形式中,使用初始刚度张量 C0\mathbf{C}_0C0​,它是 σvirgin=C0:εe\boldsymbol{\sigma}_{\text{virgin}} = \mathbf{C}_0 : \boldsymbol{\varepsilon}^eσvirgin​=C0​:εe,其中 εe\boldsymbol{\varepsilon}^eεe 是弹性应变。

应变等效原理告诉我们,只需将原始应力替换为有效应力:

σ~=C0:εe\tilde{\boldsymbol{\sigma}} = \mathbf{C}_0 : \boldsymbol{\varepsilon}^eσ~=C0​:εe

现在我们使用两个关键方程。我们将 σ~=σ/(1−D)\tilde{\boldsymbol{\sigma}} = \boldsymbol{\sigma}/(1-D)σ~=σ/(1−D) 代入上式:

σ1−D=C0:εe\frac{\boldsymbol{\sigma}}{1-D} = \mathbf{C}_0 : \boldsymbol{\varepsilon}^e1−Dσ​=C0​:εe

重新整理以求解真实的、可测量的应力 σ\boldsymbol{\sigma}σ,我们得到受损材料的本构律:

σ=(1−D)C0:εe\boldsymbol{\sigma} = (1-D) \mathbf{C}_0 : \boldsymbol{\varepsilon}^eσ=(1−D)C0​:εe

看!所有这些微观空洞和裂纹的复杂效应被归结为一个简单的乘法因子 (1−D)(1-D)(1−D)。受损材料的有效刚度 C(D)\mathbf{C}(D)C(D) 仅为 C(D)=(1−D)C0\mathbf{C}(D) = (1-D)\mathbf{C}_0C(D)=(1−D)C0​。对于一根具有杨氏模量 E0E_0E0​ 的简单杆,其损伤模量变为 Ed=(1−D)E0E_d = (1-D)E_0Ed​=(1−D)E0​。随着损伤 DDD 从 0 变为 1,刚度从 E0E_0E0​ 降至零。该模型完美地捕捉了材料的弱化过程。

有趣的是,这个简单的各向同性模型预测,虽然刚度下降,但泊松比——衡量材料拉伸时变薄程度的指标——保持不变。这是该模型的一个具体预测,可以在实验中进行检验。

我们也可以从柔度 S\mathbf{S}S 的角度来看待这个问题,它是刚度的倒数。损伤柔度变为 S(D)=S0/(1−D)\mathbf{S}(D) = \mathbf{S}_0 / (1-D)S(D)=S0​/(1−D)。这完全合乎逻辑:受损材料的刚度较低,这意味着在给定的施加应力下,它更具柔性,或更容易“拉伸”。

更深层的法则:为何热力学为王

这个等效原理可能看起来像是一种巧妙的唯象建模,一个方便的猜测。但它真正的美在于它与物理学最基本定律——热力学定律的深刻联系。

让我们思考能量。当我们使弹性材料变形时,我们在其中储存了势能,就像拉伸弹簧一样。这种储存的能量由一个名为​​亥姆霍兹自由能​​的函数 ψ\psiψ 来描述。在热力学的语言中,应力就是这个能量函数对畸变的导数:σ=∂ψ∂εe\boldsymbol{\sigma} = \frac{\partial\psi}{\partial\boldsymbol{\varepsilon}^e}σ=∂εe∂ψ​。

如果我们的损伤胡克定律 σ=(1−D)C0:εe\boldsymbol{\sigma} = (1-D) \mathbf{C}_0 : \boldsymbol{\varepsilon}^eσ=(1−D)C0​:εe 是正确的,那么能量函数必须是什么样的?我们可以通过积分来反向推导。结果同样优雅:

ψ(εe,D)=(1−D)(12εe:C0:εe)=(1−D)ψ0(εe)\psi(\boldsymbol{\varepsilon}^e, D) = (1-D) \left(\frac{1}{2} \boldsymbol{\varepsilon}^e : \mathbf{C}_0 : \boldsymbol{\varepsilon}^e \right) = (1-D) \psi_0(\boldsymbol{\varepsilon}^e)ψ(εe,D)=(1−D)(21​εe:C0​:εe)=(1−D)ψ0​(εe)

其中 ψ0\psi_0ψ0​ 是原始无损材料的自由能。这告诉我们,受损材料的储能能力就是原始材料的能力乘以一个因子 (1−D)(1-D)(1−D)。我们“瑞士奶酪”中的“孔洞”不能储存能量,所以只有完好的那部分材料贡献了能量。

值得注意的是,人们本可以从一个不同的起点出发,即​​能量等效假说​​,它直接假设了这种能量形式。对于这种简单的各向同性损伤情况,这两个原理——一个从应力和应变出发,另一个从能量出发——得出了完全相同的结果。不同物理论证的这种殊途同归,是稳健理论的标志。

失效的引擎:损伤能量释放率

是什么导致损伤增长?直观上,材料通过断裂来释放储存的能量。热力学将这一思想形式化。热力学第二定律(以克劳修斯-杜亥姆不等式的形式)要求任何不可逆过程,如损伤,都必须导致非负的能量耗散。

与损伤相关的耗散由乘积 D=YD˙\mathcal{D} = Y \dot{D}D=YD˙ 给出,其中 D˙\dot{D}D˙ 是损伤增长率,而 YYY 是其​​热力学共轭力​​,称为​​损伤能量释放率​​。第二定律要求 YD˙≥0Y \dot{D} \ge 0YD˙≥0。

这个驱动力 YYY 是什么?热力学框架给了我们一个精确的定义:它是自由能对损伤的负偏导数,Y=−∂ψ∂DY = -\frac{\partial\psi}{\partial D}Y=−∂D∂ψ​。让我们为我们的模型计算它:

Y=−∂∂D[(1−D)ψ0(εe)]=−[(−1)ψ0(εe)]=ψ0(εe)Y = - \frac{\partial}{\partial D} \Big[ (1-D) \psi_0(\boldsymbol{\varepsilon}^e) \Big] = - \Big[ (-1) \psi_0(\boldsymbol{\varepsilon}^e) \Big] = \psi_0(\boldsymbol{\varepsilon}^e)Y=−∂D∂​[(1−D)ψ0​(εe)]=−[(−1)ψ0​(εe)]=ψ0​(εe)

这个结果在其简洁性和物理丰富性上令人惊叹。在任何瞬间,驱动材料产生更多损伤的热力学力,等于其原始无损状态在相同应变下储存的弹性能。你拉伸材料越多,它“渴望”断裂以释放储存能量的意愿就越大。由于原始刚度 C0\mathbf{C}_0C0​ 是正定的,ψ0\psi_0ψ0​ 总是非负的。这意味着 Y≥0Y \ge 0Y≥0。结合损伤只能增长 (D˙≥0\dot{D} \ge 0D˙≥0) 的物理事实,热力学第二定律 YD˙≥0Y \dot{D} \ge 0YD˙≥0 总是得到满足。这个模型在热力学上是合理的。

了解局限:当简单模型不再足够

像科学中的任何模型一样,这种简单标量形式的应变等效原理完美地解释了很多现象,但不是所有。理解其局限性至关重要。

  • ​​拉伸 vs. 压缩:​​ 我们的模型使用单一数字 DDD 来描述材料的状态。它预测无论你是拉(拉伸)还是压(压缩),刚度都会同样降低。但对于许多材料,如混凝土,在拉伸时张开的裂纹在压缩时可以闭合并传递力。材料在这两种情况下的行为不同——这是一个我们的简单模型无法捕捉的​​单边效应​​。

  • ​​各向异性:​​ 如果我们的材料在所有方向上并非都相同怎么办?考虑一块木板或碳纤维复合材料。像纤维间开裂这样的损伤,在垂直于纤维的方向上对材料的削弱程度远大于沿纤维方向。这种损伤是​​各向异性​​的。单一标量 DDD 是不够的;我们需要一个更复杂的张量来描述损伤。

  • ​​其他耗散机制:​​ 一些材料,如土壤或颗粒材料,有其他耗散能量的方式。颗粒可以相互滑动,产生摩擦损失。这种摩擦耗散是一个与我们损伤变量 DDD 所模拟的新表面积产生不同的过程。

机器中的幽灵:局部化之谜

最后,这个优雅的局部模型在我们在计算机模拟中使用它时,隐藏了一个有趣且棘手的难题。方程预测,一旦材料开始“软化”(即当应力随着应变的增加而开始减小时),所有的变形都将迅速集中到一个无限薄的带中。

在使用有限元法的计算机模型中,这个“无限薄的带”变成了一个只有一个单元宽的带。当你为了获得更精确的答案而细化网格时,失效区域只会变得越来越窄。其病态的后果是,当网格尺寸趋于零时,预测破坏试件所需的总能量也趋于零!这显然是不符合物理的——产生新的裂纹表面需要有限的能量,即​​断裂能​​。

这种严重的​​网格敏感性​​表明,局部模型缺少了一个关键的物理要素:一个​​内禀长度尺度​​。该模型没有关于应该控制失效区宽度的微观结构特征尺寸的信息。为了解决这个问题,科学家们开发了“非局部”损伤模型,这是另一个话题了。这些先进的模型引入了一个特征长度,确保计算出的断裂能是真实的材料属性,与计算网格无关。这是一个绝佳的例子,说明一个简单模型的局限性如何为更深入、更完整的理论指明方向。

应用与跨学科联系

我们花了一些时间从其自身角度了解应变等效原理,探索了其优雅的热力学和力学基础。但物理学中的一个原理的价值在于它能描述的世界。一个真正伟大的原理不会孤立存在;它会延伸出去,连接不同的现象,并赋予我们理解、预测和创造的新能力。它的美不仅在于其公式,更在于其应用。

现在,我们将踏上一段旅程,看看这个简单而优雅的思想——即一个受损材料的行为就像一个健康的材料在虚构的“有效”应力下的行为——如何展现出一幅丰富的应用图景。我们将看到它如何成为设计桥梁的工程师、诊断失效的材料科学家以及模拟裂纹萌生的计算物理学家不可或缺的工具。这就是该原理焕发生机的地方。

工程师的工具箱:预测材料行为

普适的软化效应

应变等效原理最直接、最实际的推论是,损伤使材料变得“更软”或更具柔性。想象你有一块坚固的钢块。如果你用一定的力拉它,它会伸长一个小的、可预测的量,这个量由其杨氏模量 E0E_0E0​ 决定。现在,假设这块钢块布满了不可见的微裂纹网络。损伤变量 DDD 给了我们一个衡量材料受损程度的指标。当我们用相同的力拉伸这个受损的钢块时,它伸长得更多。为什么?

该原理告诉我们,观测到的应力 σ\sigmaσ 和应变 ε\varepsilonε 通过一个有效模量 E=(1−D)E0E = (1-D)E_0E=(1−D)E0​ 相关联。对于给定的应力,应变为 ε=σ/E=σ/((1−D)E0)\varepsilon = \sigma / E = \sigma / ((1-D)E_0)ε=σ/E=σ/((1−D)E0​)。由于 DDD 是一个正数,(1−D)(1-D)(1−D) 小于1,因此产生的应变比在未损伤材料中要大。材料实际上变得不那么刚硬了。

当我们考察更复杂的变形时,这个思想的真正力量就显现出来了。扭转(剪切)或静水压缩呢?应变等效原理,在其最简单的标量形式下,做出了一个大胆而有力的断言:所有的弹性性能都同步退化。剪切模量 G0G_0G0​ 变为 G=(1−D)G0G = (1-D)G_0G=(1−D)G0​,体积模量 K0K_0K0​ 变为 K=(1−D)K0K = (1-D)K_0K=(1−D)K0​。这意味着整个四阶刚度张量 C0\mathbf{C}_0C0​(包含了材料所有弹性响应的信息)被均匀地缩减:C(D)=(1−D)C0\mathbf{C}(D) = (1-D)\mathbf{C}_0C(D)=(1−D)C0​。这是一个巨大的简化!工程师不再需要追踪一般各向异性材料的21个独立弹性常数的演化,而是可以作为初步近似,追踪一个单一的标量损伤变量 DDD。

解读材料的“思想”:实验识别

这个原理不仅仅是预测已知损伤材料行为的单向通道。我们可以反向运用逻辑来做一些更深刻的事情:我们可以利用我们在外部测量到的行为来诊断内部累积的损伤。

想象一个材料测试实验室。我们将一个新复合材料的试样放入一台拉伸它的机器中,精确记录每一刻的应力 σ\sigmaσ 和应变 ε\varepsilonε。最初,在非常小的应变下,应力-应变图是一条直线,其斜率给出了原始杨氏模量 E0E_0E0​。当我们拉得更用力时,曲线开始弯曲。材料不再像以前那样刚硬了。这是因为损伤吗?

根据我们的原理,曲线上任意一点 (σ,ε)(\sigma, \varepsilon)(σ,ε) 的损伤可以直接计算出来。如果我们定义一个“割线”模量 Esec=σ/εE_{sec} = \sigma/\varepsilonEsec​=σ/ε,即从原点到该点的直线的斜率,我们的本构律 σ=(1−D)E0ε\sigma = (1-D)E_0 \varepsilonσ=(1−D)E0​ε 可以重新整理为:

D=1−σE0ε=1−EsecE0D = 1 - \frac{\sigma}{E_0 \varepsilon} = 1 - \frac{E_{sec}}{E_0}D=1−E0​εσ​=1−E0​Esec​​

突然之间,我们有了一个窥探材料内部状态的窗口。通过简单地测量应力-应变曲线,我们就可以绘制出材料加载过程中损伤变量 DDD 的演化图。曾经是一个抽象的内部变量,现在成了一个可测量的量,一个材料的“健康计”。

损伤与塑性:两种变形的故事

当材料的应力-应变曲线偏离其初始直线时,并不总是因为损伤。许多材料,尤其是金属,也可以发生塑性变形。想象一下弯曲一个回形针:它不会弹回原来的形状。这种永久变形就是塑性。那么,我们如何区分这两种现象呢?在初始加载曲线上,刚度下降和永久变形可能看起来相似。

在这里,一个巧妙的实验技术,结合我们的原理,让我们能够扮演侦探。我们可以将材料加载到一定程度,然后卸载。

  • 如果材料沿着与初始加载曲线斜率相同 (E0E_0E0​) 的直线卸载,但在零应力时留下永久应变,那么它经历了纯塑性变形。它没有失去刚度,只是被永久地重塑了。
  • 如果材料沿着更缓的斜率 (E<E0E < E_0E<E0​) 的直线卸载,并在零应变处回到零应力,那么它经历了纯损伤。它失去了刚度,但没有永久变形。
  • 如果它沿着更缓的斜率卸载并且在零应力时有永久应变,那么两种现象都在起作用:损伤和塑性同时发生。

应变等效原理为我们提供了解释的关键:卸载线的斜率告诉我们新的、损伤后的模量 E=(1−D)E0E = (1-D)E_0E=(1−D)E0​,而零应力处的应变偏移揭示了塑性应变 εp\varepsilon^pεp。这种解开两种复杂、重叠机制的能力是现代实验力学的基石。

搭建桥梁:与其他学科的联系

一个基本原理的真正考验是它与科学和工程其他领域建立联系,形成一个更完整、更强大的世界观的能力。

损伤与塑性的结合

我们已经看到损伤和塑性可以共存。应变等效原理提供了一个严谨且热力学一致的框架来共同模拟它们。由 Jean Lemaitre 等人形式化的关键洞见是,假设塑性定律不是在可测量的“柯西”应力 σ\boldsymbol{\sigma}σ 的世界中运作,而是在“有效”应力 σ~\tilde{\boldsymbol{\sigma}}σ~ 的世界中运作。

总应变 ε\boldsymbol{\varepsilon}ε 被分为弹性部分 εe\boldsymbol{\varepsilon}_eεe​ 和塑性部分 εp\boldsymbol{\varepsilon}_pεp​。弹性应变由损伤模量决定,但从有效应力的角度思考更清晰:

σ~=σ1−D和εe=S0:σ~\tilde{\boldsymbol{\sigma}} = \frac{\boldsymbol{\sigma}}{1-D} \quad \text{和} \quad \boldsymbol{\varepsilon}_e = \mathbf{S}_0 : \tilde{\boldsymbol{\sigma}}σ~=1−Dσ​和εe​=S0​:σ~

其中 S0\mathbf{S}_0S0​ 是无损材料的柔度。材料的屈服准则——决定其何时开始塑性流动的条件——被写成有效应力的函数,f(σ~,… )≤0f(\tilde{\boldsymbol{\sigma}}, \dots) \le 0f(σ~,…)≤0。

就好比材料的塑性流动机制对损伤是“盲目”的。它只“感受”到由其微观结构中仍然完好的部分所承载的应力。损伤变量 DDD 在这个有效应力的内部世界和我们测量的名义应力的外部世界之间,起到了一个面纱或过滤器的作用。这个优雅的框架使我们能够建立强大的预测模型,捕捉延性材料从屈服、硬化,到同时因损伤累积而退化,直至断裂点的完整应力-应变曲线。

结构完整性与微小缺陷的“陷阱”

让我们从一个小的材料样本转移到一个大的工程结构,比如飞机机翼或桥梁。这些结构不可避免地包含应力集中点:铆钉孔、尖角或焊缝。经典弹性理论告诉我们,在受拉平板中,圆孔边缘的应力可以达到远离孔洞处应力的三倍。现在,让我们看看当这个平板的材料存在一些预先存在的、均匀的损伤 DDD 时会发生什么。

乍一看,人们可能认为各处的应力都简单地减少了一个因子 (1−D)(1-D)(1−D)。但应变等效原理揭示了一个更危险的现实。受远场应力 σ0\sigma_0σ0​ 作用的损伤平板问题,在数学上等同于一个无损平板在更高的远场有效应力 σ~0=σ0/(1−D)\tilde{\sigma}_0 = \sigma_0 / (1-D)σ~0​=σ0​/(1−D) 下的问题。

因此,孔边缘的最大有效应力是 3σ~0=3σ0/(1−D)3 \tilde{\sigma}_0 = 3 \sigma_0 / (1-D)3σ~0​=3σ0​/(1−D)。而且因为应变与有效应力成正比,孔口处的最大应变也被这个因子放大了:

εθθmax⁡=3σ0E0(1−D)\varepsilon_{\theta\theta}^{\max} = \frac{3 \sigma_0}{E_0 (1-D)}εθθmax​=E0​(1−D)3σ0​​

这是一个关键结果。一个看似无害的背景损伤,比如说 D=0.5D=0.5D=0.5(意味着材料失去了一半的刚度),不仅仅是将强度减半。在应力集中点,它使给定施加载荷下的应变加倍!损伤极大地放大了几何缺陷的影响,这解释了为什么老化的结构会在以前被认为是安全的位置突然意外失效。

从微小裂纹到统一规律:微观力学一瞥

到目前为止,我们一直将损伤变量 DDD 视为一个唯象量。但它从何而来?我们能将它与材料中发生的实际物理变化联系起来吗?在这里,该原理与微观力学领域建立了一座美丽的桥梁。

考虑一个包含稀疏、随机分布的微小、不相互作用的微裂纹的固体。每个裂纹都使其附近的材料更具柔性。微观力学模型允许我们通过平均所有这些个体缺陷的影响来计算这个含裂纹固体的整体或“有效”刚度。令人惊奇的结果是,在随机裂纹取向的假设下,含裂纹固体的有效刚度张量可以一阶近似为:

Ceff≈C0−ηH\mathbf{C}_{\text{eff}} \approx \mathbf{C}_0 - \eta \mathbf{H}Ceff​≈C0​−ηH

其中 η\etaη 是裂纹密度参数,H\mathbf{H}H 是一个依赖于裂纹几何形状和固体性质的四阶张量。在裂纹对材料抗压和抗剪能力的影响恰好成比例的特殊(但不普遍)情况下,这个复杂的表达式简化为应变等效原理的精确形式,Ceff≈(1−D)C0\mathbf{C}_{\text{eff}} \approx (1-D)\mathbf{C}_0Ceff​≈(1−D)C0​。这为我们的唯象定律提供了物理基础,表明它可以被解释为大量微观缺陷的宏观表现。

驯服无穷:向非局部模型的飞跃

最后,我们进入计算力学的前沿。当我们尝试使用局部应变等效原理来模拟材料失效时,会遇到一个奇怪且灾难性的问题。随着损伤的累积,材料软化。应变会自然地集中在材料最薄弱的部分。这种应变集中导致更多的损伤,从而导致更多的软化,进而导致更多的应变集中。这种恶性循环导致损伤在我们的计算机模型中局部化到一个零厚度的带中,这是一个不符合物理的结果,会导致模拟崩溃。

解决方案是一个极其优雅的想法,它在不放弃其核心的情况下改进了我们的原理。我们保持应力定律的局部性:一个点的应力只取决于该点的应变和损伤。然而,我们改变了损伤演化的定律。点 xxx 处的损伤累积速率不再由点 xxx 处的应变决定,而是由 xxx 周围一个小邻域内应变的加权平均值决定:

εˉ(x)=∫Ωw(∣x−ξ∣)ε~(ξ) dξ∫Ωw(∣x−ξ∣) dξ\bar{\varepsilon}(x) = \frac{\int_{\Omega}w(|x-\xi|)\tilde{\varepsilon}(\xi)\,\mathrm{d}\xi}{\int_{\Omega}w(|x-\xi|)\,\mathrm{d}\xi}εˉ(x)=∫Ω​w(∣x−ξ∣)dξ∫Ω​w(∣x−ξ∣)ε~(ξ)dξ​

xxx 处的损伤现在取决于其邻居的经历。这种“非局部”相互作用阻止了损伤局部化到单一平面。它迫使失效区具有有限的宽度,该宽度与权重函数 www 的特征长度相关。这不仅解决了数值问题,而且在物理上也更现实,反映了内部失效过程(如微裂纹)涉及在微小但有限的距离上的相互作用。

一个力量惊人的原理

我们的旅程结束了。我们从一个简单的尺度假设开始,见证了它发展成为一个多功能且强大的工具。它使我们能够预测工程材料的行为,设计窥探其内部状态的实验,建立稳健的结构失效理论,并开发能够模拟复杂断裂过程的计算模型。

应变等效原理,以其清晰和实用性,完美地诠释了伟大的科学思想如何不仅提供答案,而且提供一种提出更好问题的新语言。它统一了微观与宏观、理论与实践,并在此过程中,揭示了一个更深刻、更关联的物质世界图景。