为什么经典断裂力学算不准？——从"无限尖裂纹"到"真实物理过程"的范式转变

经典断裂力学假设裂纹尖端是数学上的"无限尖点"，导致应力/应变出现非物理的奇异性，且完全忽略缺陷尺度对承载能力的影响。新理论通过"均匀化能量密度"框架，证明裂纹尖端变形实际上是非奇异的，并能客观预测缺陷尺寸效应，为准脆性材料的极限承载能力评估提供了物理上一致的方法。

一、经典断裂力学的"阿喀琉斯之踵"

1.1 数学尖点 vs 物理现实

1913年，Inglis分析了含裂纹无限大板的应力集中问题，奠定了线弹性断裂力学（LEFM）的基础。其核心假设是：裂纹尖端是数学上的尖点（半径为零）。

这个假设带来了严重的物理矛盾：

• 应力奇异性：根据LEFM，裂纹尖端的应力随  发散，理论上趋于无穷大，裂纹扩展将无任何阻力应变奇异性：位移梯度在尖端处同样发散

1.2 无法解释的"尺寸效应"

经典理论的另一个致命缺陷是无法考虑缺陷尺度对承载能力的影响：

这些现象表明，当结构特征尺寸（孔径、裂纹长度、梁厚度）与材料的特征微观尺度（晶粒尺寸、骨料粒径、高分子链尺度）相当时，经典理论彻底失效。

二、均匀化能量密度理论：从"点"到"体"的范式革命

2.1 核心洞察：RVE不是数学点

均匀化能量密度模型（Homogenized Energy Density Model）从根本上改变了这一局面。

关键认识：经典理论将"代表性体积单元"（RVE）视为无限小的数学点，但实际上，真实的RVE具有有限的物理尺寸 h （如金属晶粒尺寸、混凝土骨料粒径、聚合物回转半径）。

当变形场在RVE内呈非线性分布时，体积平均值 ≠ 几何中心值：

其中 Δp 是拉普拉斯算子，描述场的"弯曲程度"。经典理论只保留了第一项，忽略了  和  项——这就是"均匀化误差"的来源。

2.2 退化应变能密度：损伤与变形的统一描述

对于准脆性断裂问题，理论引入了退化均匀化应变能密度（Degraded Homogenized SED）：

其中  是局部退化应变能密度，d 为损伤变量。

这一公式的意义在于：

1. 单参数描述：仅需RVE尺寸 h 一个尺度参数，具有明确物理意义自然正则化：高阶项自动正则化应变/损伤局部化，无需人为引入宽度参数客观预测尺寸效应：h 与结构特征尺寸的比值直接决定尺寸效应的强弱

三、裂纹尖端变形是非奇异的

在2025年发表于Computer Methods in Applied Mechanics and Engineering的论文中通过渐近分析验证了裂纹尖端的变形是非奇异的：

I型裂纹尖端曲率半径对应变的影响

当裂纹尖端半径 R 从  mm减小到  mm 时：

• 经典局部理论：应变  随  发散，趋于无穷（黑色方块）新理论预测：应变趋于有限值，稳定在约0.025（红色圆点）

3.2 物理机制：高阶变形主导

关键发现：随着裂纹变尖，宏观应变能密度 U 几乎不变，而总均匀化应变能密度  显著增加。

这意味着：

裂纹尖端的内部功主要由高阶微观变形（应变梯度项）贡献，而非传统宏观变形。因此，宏观应变被约束在有限值，奇异性被自然消除。

这一结果与经典内聚力模型（CZM）有本质区别——CZM通过人为引入"过程区"长度参数来消除奇异性，而新理论从能量均匀化的物理本质出发，无需额外假设即可得到非奇异解。

四、拉压不对称性的能量描述

2026年发表于International Journal of Engineering Science的论文进一步解决了准脆性材料的核心特征——拉压不对称性。

4.1 原始模型的局限

第一篇论文采用修正von Mises等效应变作为损伤准则：

其中 k 是拉压强度比，需要通过实验标定。这带来了问题：

• 当高阶变形显著时，局部拉压状态复杂，k 的标定变得模糊不同变形模式下，k 可能需要取不同值

4.2 谱分解：从"强度准则"到"能量准则"

改进后的模型采用谱分解（Spectral Decomposition）策略，将应变能密度分解为拉伸和压缩部分：

其中：

•  —— 仅当主应变为正时激活 —— 仅当主应变为负时激活

关键改进：损伤仅退化拉伸部分的能量，压缩部分保持完好：

这一改进的物理意义：

1. 消除经验参数：不再需要标定拉压强度比 k物理一致性：裂纹扩展由拉伸变形驱动，压缩变形提供约束——这与混凝土、岩石等准脆性材料的实际破坏机制完全一致高阶项的拉压不对称：高阶均匀化误差项同样进行谱分解，确保微观尺度上的拉压不对称性被正确传递至宏观

4.3 验证：复杂裂纹路径预测

在非对称缺口梁三点弯曲试验中，改进模型展现出显著优势：

这是因为裂纹尖端经历了复杂的拉-压-剪混合状态，谱分解模型能准确捕捉拉伸变形主导的损伤演化，而强度参数模型无法适应这种复杂的应力状态变化。

五、理论框架的完整图景

5.1 两尺度耦合机制

微观尺度（RVE内部）          宏观尺度（结构响应）
     ↓                           ↓
 应变/损伤非均匀分布  ←——h——→  均匀化能量密度
     ↓                           ↓
 高阶梯度项（∇ε, ∇²ε）      退化刚度矩阵
     ↓                           ↓
 拉伸/压缩谱分解            损伤演化准则
     ↓                           ↓
 非局部相互作用              客观尺寸效应预测

5.2 与现有方法的对比

六、工程应用价值

6.1 带孔板极限承载力预测

实验表明，PMMA带孔板的名义断裂强度随孔径减小而显著增加，均匀化能量理论能定量捕捉这一尺寸效应。

方形孔尺寸对名义拉伸强度的影响

关键发现：当孔径与RVE尺寸相当时，高阶变形项显著贡献于总能量，导致表观"强化"效应。

6.2 混凝土结构裂缝扩展

在L型板、双边缺口试件等经典benchmark问题中，改进模型准确预测了：

• 裂纹从缺口萌生的位置拉-剪混合模式下的裂纹偏转峰值荷载后的软化行为

特别地，模型成功捕捉了DEN试件中法向荷载由拉转压的复杂过程，这是检验拉压不对称处理能力的严苛测试。

(a) 实验 （b)基于能量分解准则 (c)基于等效应变准则

力-位移曲线

(a)实验 （b)基于能量分解准则 (c)基于等效应变准则

(a) 拉伸力 (b)剪切力

七、结语：从"数学技巧"到"物理真实"

均匀化能量密度理论代表了断裂力学研究范式的转变：

从"如何数学上处理奇异性"转向"如何从物理上消除奇异性"。

通过认识到RVE的有限尺寸和内部场的不均匀性，理论自然导出了：

1. 非奇异的裂纹尖端变形——符合物理真实客观的缺陷尺寸效应——无需经验拟合拉压不对称的能量描述——基于变形本质而非强度假设

这一框架不仅为准脆性材料的极限承载能力预测提供了可靠工具，更重要的是，它建立了微观结构特征（RVE尺寸）与宏观结构响应之间的客观联系，为材料-结构一体化设计奠定了理论基础。

参考文献：

Cao YH, Zhang CY. A Two-scale High-order Damaged Elasticity Theory and Solution Procedure for Quasi-brittle Fracture. Computer Methods in Applied Mechanics and Engineering, 445 (2025) 118206.Cao YH, Zhang CY. Damaged Homogenized Energy Model based on a Decomposition of Higher-order Strain Energy Density. International Journal of Engineering Science, 223 (2026) 104498.