3. 损伤力学基本概念

损伤力学研究的是材料在载荷、温度或环境作用下，内部微缺陷（微裂纹、微孔洞）的萌生、扩展与聚合过程，以及这些缺陷如何导致材料宏观刚度和承载能力的逐步退化。与断裂力学关注“一条主裂纹是否失稳扩展”不同，损伤力学关注的是“材料整体还有多少有效承载能力”。在有限元工程中，损伤变量就是把这种退化引入本构方程、从而让仿真结果能反映刚度衰减和破坏路径的关键桥梁。

损伤变量的定义

损伤变量 D 是一个无量纲的内变量，用来量化材料内部缺陷造成的承载面积损失。最直观的定义来自 Kachanov-Rabotnov 的连续损伤力学框架：设材料某个截面上的名义（表观）面积为 A，其中真正能传递应力的有效承载面积为 A_eff，则

D = (A − A_eff) / A，   0 ≤ D ≤ 1

当 D = 0 时材料完好无损；D = 1 时有效承载面积为零，材料完全失去承载能力，对应宏观裂纹形成或单元失效。工程中通常取一个临界损伤值 D_c（如 0.2～0.5，视材料和判据而定）作为失效阈值，而不是等到 D = 1。损伤变量可以是标量（各向同性损伤）、矢量或二阶张量（各向异性损伤），其选择取决于缺陷的取向性和载荷路径。

有效应力与应变等效假设

损伤的力学后果通过“有效应力”体现。名义应力 σ = F/A，而真正作用在有效承载面积上的应力为

σ̃ = F / A_eff = σ / (1 − D)

这就是有效应力公式。它说明：随着损伤累积，同样的外载会在剩余材料上产生越来越高的局部应力，形成“损伤—应力升高—损伤加速”的正反馈，这正是损伤演化后期急剧失稳的根源。

应变等效假设（Lemaitre 提出）进一步给出本构关系的处理方式：受损材料的应变响应，等同于用有效应力 σ̃ 作用在无损材料上所产生的应变。以一维线弹性为例，无损时 ε = σ/E；受损后

ε = σ̃ / E = σ / [E(1 − D)]

等价于把弹性模量退化为 Ẽ = E(1 − D)。这一假设的价值在于：无需重新推导整套本构，只需把无损本构中的应力替换为有效应力，即可得到损伤耦合本构。在有限元中，这直接体现为刚度矩阵随 D 增大而折减，例如各向同性损伤下 K_eff = (1 − D)K₀。

各向同性损伤与各向异性损伤

当微缺陷在空间上近似均匀分布、无明显取向时，可用标量 D 描述，称为各向同性损伤。它的优点是本构简单、参数少、数值实现稳定，适合蠕变损伤、韧性损伤初期、多晶金属的均匀孔洞化等场景。此时刚度退化对各个方向相同：C_eff = (1 − D)C₀。

当缺陷具有明显取向（如层合复合材料的分层、单向纤维的基体开裂、轧制板材的方向性损伤），标量 D 无法区分不同方向的退化程度，必须采用张量型损伤变量 D_ij 或 D_ijkl。工程上常用二阶张量 D，其主值对应三个主损伤方向，本构中通过损伤张量对刚度张量进行定向折减。各向异性损伤精度更高，但参数标定复杂、计算量大，且需要更多实验（如不同方向的模量衰减测试）来识别。

损伤演化描述

损伤变量不是常数，它随载荷历史单调增长。演化方程通常写成损伤驱动力（如等效塑性应变、能量释放率、蠕变应变）的函数。以 Lemaitre 韧性损伤模型为例，其演化可写为

D_dot = (Y / S) * (ε_p_dot)  当 ε_p ≥ ε_D
D = 0                          当 ε_p < ε_D

其中 Y 为损伤应变能释放率，S 为材料损伤参数，ε_D 为损伤起始阈值。工程实现时，损伤演化常以分段形式给出：损伤起始判据（如等效塑性应变达到阈值）、演化律（线性、指数或幂律）、失效判据（D 达到 D_c）。在有限元中，每个积分点独立更新 D，再通过刚度折减影响整体响应，因此损伤场往往呈现局部化特征，需要配合网格正则化（如特征长度、黏性正则化）避免网格依赖。

理解这四个概念之间的逻辑链——损伤变量定义退化程度，有效应力把退化转化为应力放大，应变等效假设把有效应力接入本构，各向同性/各向异性决定描述精度，演化方程决定退化如何随时间发展——是后续在有限元中正确设置损伤模型、解读刚度衰减和破坏路径的基础。