3. 损伤力学基础与Abaqus材料模型
在连续介质损伤力学中，我们不再把材料视为“完好”或“完全破坏”的二值状态，而是引入一个连续的内变量——损伤变量，用它描述材料内部微裂纹、微孔洞的累积程度。Abaqus 中的损伤本构正是围绕这一变量展开：先定义损伤如何度量，再定义损伤如何影响刚度，最后定义损伤如何随载荷演化。
一、损伤变量与有效应力
最常用的损伤变量是标量 D，取值范围为 0 到 1。D=0 表示材料完好，D=1 表示材料完全失去承载能力。其物理含义可以用“有效承载面积”来理解：设材料初始截面积为 A，损伤后实际能传递应力的有效面积为 A_eff = A(1−D)。
由此引出有效应力的概念。名义应力 σ = F/A，而有效应力 σ_eff = F/A_eff = σ/(1−D)。也就是说，损伤使材料内部实际承受应力的面积减小，有效应力被放大。这一关系是损伤力学与 Abaqus 损伤本构之间的桥梁：Abaqus 在计算中通过刚度退化间接体现有效应力的放大效应。
关键关系：σ_eff = σ / (1 − D)，当 D→1 时 σ_eff→∞，意味着材料即将失去承载能力，数值计算中需要引入单元删除或黏性正则化来避免收敛困难。
二、刚度退化与损伤本构
损伤对材料最直接的影响是刚度退化。在弹性损伤框架下，损伤后的弹性模量可写为 E_d = (1−D)E₀，其中 E₀ 为初始弹性模量。相应地，应力—应变关系变为 σ = (1−D)E₀ε。这意味着随着损伤累积，同样的应变增量只能产生更小的应力增量，材料逐渐“变软”。
在 Abaqus 中，损伤通常与塑性或弹性本构耦合。以延性损伤为例，Abaqus 先按未损伤的弹塑性本构计算应力，再根据损伤演化变量对刚度进行折减。具体实现上，Abaqus 采用损伤因子 d 对刚度矩阵进行缩放，其效果与上述 (1−D) 一致。需要注意的是，Abaqus 中损伤变量与塑性应变、等效塑性应变或断裂能等内部变量绑定，而不是独立输入。
三、损伤演化与 Abaqus 输入框架
损伤演化描述 D 如何随加载历史增长。Abaqus 提供多种演化定义方式，常见的有：
- 基于等效塑性应变的演化：适用于延性金属，损伤起始后 D 随等效塑性应变线性或指数增长。
- 基于断裂能的演化：适用于韧性断裂，用能量释放率控制损伤增长，减少网格依赖性。
- 基于位移的演化：适用于界面或 cohesive 单元，损伤随特征位移增长。
在 Abaqus/CAE 中，损伤本构的输入框架通常分为三步：
- 定义损伤起始准则：在 Material → Damage for Ductile Metals 中设置等效塑性应变、应力三轴度等起始参数。
- 定义损伤演化：选择 Type 为 Energy 或 Displacement，输入断裂能或位移值，并设置软化方式（Linear、Exponential 等）。
- 指定单元删除：在 Section 中勾选 Element deletion，当 D 达到 1 时删除单元，模拟宏观裂纹。
下面这张图展示了损伤变量、有效应力、刚度退化与 Abaqus 输入框架之间的逻辑关系。图中左侧为损伤力学的物理量，右侧为 Abaqus 中对应的材料模型输入项，中间通过“刚度折减”这一核心机制连接。
四、教学示例：延性金属损伤参数设置
下面给出一个可用于教学的示例参数集，材料为某低碳钢，仅用于说明输入框架，不代表实测数据。假设弹性模量 E₀ = 210 GPa，泊松比 ν = 0.3，屈服应力 350 MPa。
*Material, name=Steel_Damage
*Elastic
210000., 0.3
*Plastic
350., 0.0
420., 0.05
480., 0.12
*Damage Initiation, criterion=DUCTILE
0.15, 0.0, 0.0
*Damage Evolution, type=ENERGY
50.0
*Section Controls, name=EC-1, hourglass=ENHANCED
1., 1., 1.
上述输入中，*Damage Initiation 的 0.15 表示等效塑性应变达到 0.15 时损伤起始；*Damage Evolution 的 50.0 表示断裂能密度为 50 mJ/mm³。在 Abaqus/CAE 中，这些参数分别对应 Ductile Damage 的 Fracture Strain 和 Damage Evolution 的 Fracture Energy。
操作提示：在 CAE 中设置损伤演化时，建议先选择 Type=Energy，再输入 Fracture Energy。若选择 Displacement，则需要输入损伤起始后的特征位移，适用于 cohesive 或界面单元。对于三维实体单元，Energy 类型通常更稳定，且对网格尺寸不敏感。
五、结果判读与常见问题
提交计算后，在 Visualization 模块中可查看 DAMAGEC（损伤变量）和 SDEG（刚度退化因子）云图。DAMAGEC 从 0 增长到 1 的过程对应损伤累积；当 DAMAGEC 达到 1 时，若启用了单元删除，该单元将从模型中移除，表现为裂纹扩展路径。
常见问题包括：损伤起始后计算不收敛、单元删除导致能量不守恒、网格依赖性过强。对应的处理方式为：引入黏性系数（Viscosity）进行正则化、使用能量型演化减少网格依赖、在 Section Controls 中设置合理的 Max Degradation 和 Element deletion 选项。若损伤变量在弹性阶段就异常增长，应检查损伤起始准则是否与塑性本构正确耦合。
注意：损伤变量 D 与 Abaqus 中的 SDEG 并不完全等同。SDEG 是刚度退化因子，其定义与具体损伤模型有关；在延性损伤中，SDEG 通常等于 D，但在某些模型中可能经过映射。查看结果时应以 Abaqus 文档中该材料模型的变量定义为准。
理解损伤变量、有效应力、刚度退化和损伤演化这四个概念之间的逻辑链条，是正确使用 Abaqus 损伤本构的前提。实际建模时，应先明确损伤起始判据，再选择演化方式，最后通过单元删除或刚度折减实现宏观失效。参数取值应结合材料试验或文献数据，并通过网格敏感性分析验证结果的稳定性。