3. 损伤力学基础：从微观裂纹到宏观失效
各位工程师朋友，大家好。今天我们聊聊损伤力学。说实话，我刚接触这个领域时，也觉得它挺抽象的。什么损伤变量、有效应力，听着像玄学。但干了几十年叶片设计，我越来越觉得，不懂损伤力学，你根本搞不清复合材料叶片到底是怎么坏的。
这一节，我们不讲虚的。我会结合自己踩过的坑，把损伤力学的核心概念掰开揉碎了讲。你想想看，一个叶片在发动机里转，几万转的转速，上千度的高温，还有复杂的交变载荷。它不会突然就断了，而是内部慢慢积累损伤，最后才失效。损伤力学，就是描述这个过程的理论工具。
3.1 损伤变量与有效应力
先问个问题：怎么量化材料内部的损伤？
我个人习惯用损伤变量 D 来定义。它很简单：
D = 1 - (A_eff / A_0)其中，A_0 是初始截面积，A_eff 是扣除微裂纹、微孔洞后的有效承载面积。D=0 表示材料完好，D=1 表示完全失效。
嗯，这里要注意。这个定义在理论上很漂亮，但实际测量很难。你总不能把叶片切开数孔洞吧？
我在项目中遇到过类似问题。当时我们做某型叶片的疲劳试验，发现寿命预测总是不准。后来我改用有效应力的概念，才把问题理清楚。
有效应力 σ_eff 的定义是：
σ_eff = σ / (1 - D)说白了，就是名义应力除以剩余有效面积比。当 D 增大时，有效应力会急剧上升，加速材料失效。这个公式虽然简单，但非常实用。
核心要点：损伤变量 D 是连接微观缺陷与宏观力学行为的桥梁。有效应力才是材料真正承受的应力。
3.2 连续损伤力学（CDM）框架
连续损伤力学，英文叫 Continuum Damage Mechanics，简称 CDM。它把损伤当作一个内变量，嵌入到本构方程中。
为什么需要这个框架？
你想想看，传统的断裂力学只研究一条主裂纹。但复合材料叶片里，损伤是分散的——基体微裂纹、纤维断裂、界面脱粘，到处都是。CDM 正好处理这种分布式损伤。
CDM 的基本思路是这样的：
- 定义损伤变量：可以是标量（各向同性损伤），也可以是张量（各向异性损伤）。
- 建立损伤演化方程：描述 D 如何随载荷、时间变化。
- 耦合损伤的本构关系：把 D 代入应力-应变关系。
- 求解边值问题：用有限元等方法计算结构响应。
我建议你记住这个框架。它就像一张地图，告诉你每一步该做什么。
下面这张图是我自己整理的 CDM 分析流程，希望能帮你建立整体概念：
实用技巧：在 ABAQUS 中实现 CDM，通常用 UMAT 子程序。我建议你先从各向同性损伤开始，跑通后再加复杂的东西。别一上来就搞张量损伤，容易把自己绕晕。
3.3 损伤演化规律
有了损伤变量和框架，下一步就是回答：损伤怎么长大？
损伤演化规律，说白了就是 D 随循环次数 N 的变化关系。常见的模型有：
- 线性演化：dD/dN = 常数。简单，但精度有限。
- 幂律演化：dD/dN = C · (Δσ)^m。适用于高周疲劳。
- 指数演化：D = 1 - exp(-k · N)。适合描述加速损伤阶段。
我曾经吃过一次亏。当时用线性模型预测某复合材料叶片的寿命，结果试验件在 60% 寿命时就断了。后来发现，复合材料损伤在后期会加速，线性模型根本抓不住这个特征。
所以我的建议是：对于复合材料，至少用两阶段模型——初期线性，后期指数。这样既简单，又能反映真实情况。
注意：损伤演化参数强烈依赖于温度、应力比和加载频率。别指望一组参数打天下。每次换材料或工况，都要重新标定。
3.4 断裂力学基础
讲完分布式损伤，我们聊聊断裂力学。它处理的是已经形成的宏观裂纹。
两个核心概念：能量释放率 G 和 应力强度因子 K。
3.4.1 能量释放率 G
G 的定义很简单：裂纹扩展单位面积所释放的能量。
G = -dΠ / dA其中 Π 是系统的总势能，A 是裂纹面积。
判断准则：当 G ≥ G_c（临界能量释放率）时，裂纹失稳扩展。
我在做叶片冲击损伤评估时，经常用 G 来判断分层是否扩展。G_c 是材料常数，需要通过试验测定。不同铺层方向的 G_c 差别很大，这一点要特别注意。
3.4.2 应力强度因子 K
K 描述裂纹尖端的应力场强度。对于 I 型（张开型）裂纹：
K_I = σ · √(π · a) · Y其中 a 是裂纹长度，Y 是几何修正因子。
判断准则：当 K_I ≥ K_IC（断裂韧性）时，裂纹扩展。
你可能会问：G 和 K 有什么关系？
对于线弹性材料，它们可以互相转换：
G = K² / E'其中 E' 是有效弹性模量（平面应力时 E' = E，平面应变时 E' = E/(1-ν²)）。
| 参数 | 物理意义 | 单位 | 适用场景 |
|---|---|---|---|
| G | 能量释放率 | J/m² | 分层、脱粘 |
| K | 应力强度因子 | MPa·√m | 基体裂纹、纤维断裂 |
| J 积分 | 弹塑性断裂参数 | J/m² | 大范围屈服 |
我的经验：在复合材料叶片中，我更喜欢用 G 而不是 K。因为复合材料是各向异性的，K 的公式推导很复杂。而 G 可以直接从有限元计算得到，方便得多。
好了，这一节的内容就到这里。损伤力学和断裂力学是叶片损伤容限分析的两大支柱。前者管微观积累，后者管宏观扩展。两者结合，才能完整描述叶片从完好到失效的全过程。
记住：理论是死的，应用是活的。多动手算几个例子，比死记公式有用得多。
公众号:蓝海资料掘金营，微信deep3321