4. 线弹性断裂力学与能量释放率G
前三章我们从应力强度因子K的角度描述了裂纹尖端的奇异性。但断裂力学的另一条主线来自能量视角：裂纹为什么会扩展？Griffith在1920年代研究玻璃强度时提出，裂纹扩展的本质是能量竞争——裂纹扩展释放弹性应变能，同时要消耗表面能。这一思想后来被Irwin和Rice发展为现代断裂力学的核心概念：能量释放率G。
一、Griffith能量平衡
考虑一块含中心穿透裂纹的无限大薄板，远场受拉应力σ，裂纹长度为2a。Griffith假设：裂纹扩展dA面积时，系统总势能Π的变化由两部分组成——弹性应变能的释放和新生表面的表面能消耗。
含裂纹体的总势能可写为：
Π = U_elastic + U_surface - W_external对于线弹性体，Irwin指出，裂纹扩展释放的弹性应变能等于裂纹存在导致的应变能增量。Griffith用Inglis的椭圆孔解得到：
U_elastic = -πσ²a²B / E' （释放的应变能，负号表示能量降低）
U_surface = 4aBγ_s （两个裂纹面，每个面面积2aB）其中B为板厚，γ_s为单位面积表面能，E'为有效弹性模量：平面应力下E'=E，平面应变下E'=E/(1-ν²)。
总能量对裂纹半长a求导，令dΠ/da=0，得到Griffith临界条件：
σ_c = sqrt(2E'γ_s / (πa))这就是著名的Griffith强度公式。它表明：裂纹越长，断裂强度越低；材料越硬（E'越大）或表面能越高，强度越高。对于玻璃等脆性材料，Griffith公式与实验吻合良好。但对于金属，实测断裂强度远高于Griffith预测，因为金属在裂纹尖端存在塑性耗散，远大于表面能。Orowan和Irwin将γ_s替换为γ_s+γ_p（塑性耗散），使能量准则适用于金属。
二、能量释放率G的定义
Griffith的全局能量平衡可以局部化。定义能量释放率G为单位裂纹面积扩展所释放的弹性应变能：
G = -dΠ/dA = - (1/B) dΠ/da （对于穿透裂纹，dA = B da）G的量纲是能量/面积，即N/m或J/m²，与应力强度因子K的量纲MPa·√m不同，但两者有确定关系。G的物理意义是“裂纹扩展的驱动力”：当G达到材料的临界值G_c（即断裂韧度）时，裂纹开始扩展。
对于Griffith裂纹（中心穿透裂纹，远场拉伸），代入应变能表达式可得：
G = πσ²a / E'令G=G_c，即得临界应力σ_c，与Griffith公式一致。这说明能量释放率G与Griffith能量平衡是同一物理规律的两种表述。
三、G与K的关系
Irwin在1957年证明了能量释放率G与应力强度因子K之间的等价关系。对于线弹性材料，裂纹尖端附近的应力场和位移场由K唯一确定，而能量释放率可由裂纹闭合积分计算。结果如下：
平面应力：G = K² / E
平面应变：G = K²(1-ν²) / E对于I型、II型、III型复合加载，总能量释放率为：
G = G_I + G_II + G_III
G_I = K_I² / E'，G_II = K_II² / E'，G_III = K_III² / (2μ)其中μ为剪切模量。这个关系极为重要：它把应力强度因子K（局部场参数）与能量释放率G（全局能量参数）统一起来。在有限元仿真中，我们既可以提取K，也可以计算G，两者应一致。如果仿真结果中G与K²/E'不吻合，通常说明网格不够细、裂纹尖端奇异性未捕捉好，或材料模型有误。
四、柔度法测G
能量释放率还可以通过柔度法实验测定。考虑含裂纹试样，裂纹长度a变化时，试样柔度C（位移/载荷）随a变化。Irwin和Kies证明：
G = (P² / 2B) · (dC/da)其中P为外载荷，B为厚度，dC/da为柔度对裂纹长度的导数。这个公式的优点是：不需要知道裂纹尖端的局部应力场，只需测量宏观的载荷-位移曲线和柔度变化，就能得到G。
柔度法的典型实施步骤：
- 制备一组相同材料、不同初始裂纹长度a的试样（如DCB双悬臂梁、CT紧凑拉伸试样）。
- 对每个试样施加单调载荷，记录载荷-位移曲线，得到柔度C=δ/P。
- 拟合C与a的关系曲线，通常C∝a^n（n为试样几何相关指数）。
- 对C(a)求导得到dC/da，代入G=(P²/2B)(dC/da)计算各裂纹长度下的G。
- 当G达到临界值G_c时，裂纹失稳扩展，对应的G_c即为材料断裂韧度。
以DCB试样为例，其柔度理论解为C=8a³/(E'Bh³)，其中h为悬臂梁半高。则dC/da=24a²/(E'Bh³)，代入得G=12P²a²/(E'B²h³)。这个表达式在复合材料层间断裂韧性测试中广泛使用。
五、能量视角下的裂纹扩展驱动力
综合来看，线弹性断裂力学提供了两条等价路径：
- 局部场路径：计算裂纹尖端应力强度因子K，与材料K_c比较，判断是否扩展。
- 能量路径：计算能量释放率G，与材料G_c比较，判断是否扩展。
两者通过G=K²/E'（平面应力）或G=K²(1-ν²)/E（平面应变）联系。在仿真中，K可通过位移外推法、J积分或相互作用积分提取；G可通过虚拟裂纹闭合技术（VCCT）或J积分计算。对于线弹性问题，J积分等于G。
能量视角的独特价值在于：它不依赖于裂纹尖端应力场的精确解，而是从全局能量平衡出发，适用于复杂几何、多材料界面和复合加载。柔度法则提供了实验测定G的实用手段，是连接理论与实验的桥梁。
在仿真实践中，建议同时输出K和G进行交叉验证。例如，在ABAQUS中通过VCCT计算G，同时用相互作用积分提取K，检查G是否等于K²/E'。若偏差超过5%，应检查裂纹尖端网格、材料参数和平面假设设置。这种交叉验证是保证断裂力学仿真可靠性的重要手段。