📚 课程目录第3章损伤变量与有效应力概念

3. 损伤变量与有效应力概念

损伤力学要回答的第一个问题是：材料“坏到什么程度”如何用一个可计算、可代入本构方程的物理量来描述。这个物理量就是损伤变量。它不是一个凭经验打分的模糊指标，而是必须能与刚度、应力、应变建立明确数学关系的状态量。本节从损伤变量的标量、矢量、张量三种描述层次出发，引出有效应力与应变等效假设，并说明各向异性损伤的表示方法。

一、损伤变量的三种描述层次

损伤的本质是材料内部微缺陷（微裂纹、微孔洞）导致承载面积减少或承载能力退化。按缺陷的方向性是否显著，损伤变量可分为三个层次。

标量损伤用单个变量 D 描述，假设损伤在各方向相同（各向同性）。最经典的几何定义来自 Kachanov-Rabotnov：设试件初始横截面积为 A₀，受损后实际能传递应力的有效承载面积为 Ã，则

D = (A₀ − Ã) / A₀ = 1 − Ã/A₀

D = 0   完全无损
D = 1   完全破坏（承载面积为零）
0 < D < 1  中间损伤状态

标量损伤适用于微孔洞型损伤、蠕变损伤、韧性损伤等方向性不强的场合，计算量最小，工程上最常用。

矢量损伤用一组带方向的损伤分量描述，适合损伤主要集中在某一方向的情形，例如单向拉伸下垂直于拉应力方向的裂纹面。它比标量多保留了方向信息，但难以描述多方向耦合。

张量损伤用二阶张量 D_ij 描述，能同时表达不同方向的损伤程度及其耦合，是各向异性损伤的通用形式。三者的关系是：标量是张量的各向同性特例，矢量是张量在单一主方向占优时的简化。

二、有效应力与应变等效假设

有效应力的概念来自承载面积减少。名义应力 σ = F/A₀ 是按初始面积算的，但真正承担载荷的只有 Ã = A₀(1−D)，因此有效应力为

σ̃ = F/Ã = F/[A₀(1−D)] = σ/(1−D)

当 D 从 0 增大到接近 1 时，σ̃ 会急剧放大，这正是损伤后期材料迅速失稳的原因。

应变等效假设（Lemaitre 提出）是损伤力学的核心桥梁：受损材料的应变响应，等同于有效应力作用在完全无损材料上的应变响应。用一维线弹性为例，无损材料 σ = Eε；受损材料则把 σ 换成 σ̃：

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

即  σ = E(1−D)ε

定义损伤弹性模量  Ẽ = E(1−D)

这个式子给出了损伤变量最直接的实验测量方法：在加载—卸载循环中测量卸载刚度 Ẽ，则 D = 1 − Ẽ/E。例如某合金初始弹性模量 E = 210 GPa，加载到某损伤状态后测得卸载刚度 Ẽ = 168 GPa，则 D = 1 − 168/210 = 0.20，即已有 20% 的刚度退化。这一方法把抽象的损伤变量与可测的宏观量直接挂钩，是损伤力学能进入工程的关键。

应变等效假设的适用前提是损伤与弹性变形解耦、损伤在代表体积单元内均匀分布。当损伤局部化严重（如裂纹尖端）时，需改用局部化或梯度增强模型。

三、各向异性损伤的张量表示

当材料在不同方向退化程度不同时，标量 D 不再够用。以正交各向异性损伤为例，用二阶对称张量 D_ij 描述，在主方向上可写成对角形式：

D = diag(D₁, D₂, D₃)

D₁, D₂, D₃ 分别为三个材料主方向的损伤分量
D₁ = D₂ = D₃ = D 时退化为各向同性标量损伤

各向异性损伤下，有效应力不再是简单除以 (1−D)，而要通过损伤张量构造四阶损伤影响张量 M，使

σ̃ = M(D) : σ

一维简化：σ̃ = σ/(1−D)
三维各向异性：σ̃ᵢⱼ = Mᵢⱼₖₗ(D) σₖₗ

工程中常用的一种简化是“损伤算子”形式，把 M 写成基于 (I − D) 的逆或对称化形式，例如

M = (I − D)⁻¹  （直接求逆形式）
或 M = [(I − D)⁻¹ + (I − D)⁻ᵀ]/2  （对称化形式）

这样，损伤张量的主方向与材料刚度退化的主方向一致，D₁ 大表示该方向承载能力退化严重。例如单向纤维复合材料横向基体开裂时，D₂（横向）显著大于 D₁（纤维方向），宏观上表现为横向模量大幅下降而纵向模量基本不变，这正是各向异性损伤描述的价值所在。

四、损伤变量的选取与工程判断

实际建模时，损伤变量的选取应遵循三条原则：一是与可测宏观量对应，优先选能通过刚度、超声波速或电阻率间接测得的量；二是与损伤机制匹配，微孔洞型用标量，定向裂纹型用矢量或张量；三是与后续本构和失效判据兼容，损伤演化方程和失效阈值都建立在同一变量上。

需要特别注意的是损伤变量的取值范围与临界值。理论上 D ∈ [0, 1]，但工程中材料往往在 D 达到 0.2～0.5 时就已失去承载意义，因此常设临界损伤 D_c 作为失效判据。例如某韧性金属 D_c = 0.35，当计算或测得 D ≥ 0.35 时判定单元失效。D_c 的取值依赖材料与工况，应由实验标定，不能随意假定。

不要把损伤变量 D 与塑性应变、孔洞体积分数混为一谈。D 是刚度退化的度量，孔洞体积分数是几何量，二者在低孔洞率下近似相关，但高孔洞率时因孔洞合并而明显偏离，需分别建模。

小结：标量、矢量、张量损伤分别对应各向同性、单一方向占优和各向异性三种情形；有效应力 σ̃ = σ/(1−D) 把承载面积减少转化为应力放大；应变等效假设 ε = σ̃/E 建立了损伤与可测刚度的桥梁，使 D = 1 − Ẽ/E 成为最实用的损伤测量公式；各向异性损伤则通过损伤张量 D_ij 和影响张量 M 把方向性引入本构。理解这三层关系，是后续建立损伤演化方程和失效判据的基础。