4. 断裂力学基本概念

在有限元损伤断裂分析中，断裂力学基本概念决定了“用什么量描述裂纹、用什么判据判断裂纹是否扩展、以及网格与后处理该提取哪些结果”。本节围绕裂纹类型、应力强度因子、能量释放率、J积分、裂纹尖端场与断裂准则展开，并给出可在有限元中直接落地的提取与判断方法。

一、裂纹类型：按受力与位移模式分类

断裂力学通常按裂纹面相对位移将裂纹分为三种基本类型：

  • I型（张开型）：裂纹面沿法向张开，受垂直于裂纹面的拉应力驱动，是最常见也最危险的类型。
  • II型（滑开型）：裂纹面在自身平面内沿切向相对滑移，由面内剪应力驱动。
  • III型（撕开型）：裂纹面沿裂纹前沿方向相对错动，由面外剪应力驱动。

工程中多数裂纹为复合型（I+II或I+II+III）。在有限元中，裂纹类型决定了裂纹前沿局部坐标系的建立方式：通常令裂纹前沿切向为局部x3，裂纹面法向为局部x2，裂纹扩展方向为局部x1，再据此分离KI、KII、KIII。

二、应力强度因子K：裂纹尖端奇异场的强度

线弹性断裂力学给出裂纹尖端应力场具有r-1/2奇异性，其强度由应力强度因子K表征。以I型为例，裂纹尖端附近应力可写为：

σ_ij(r,θ) = K_I / sqrt(2πr) · f_ij(θ) + 高阶项

其中r为到裂纹尖端的距离，θ为相对裂纹面的极角，fij(θ)为角分布函数。K的量纲为应力×长度1/2，常用MPa·mm1/2或MPa·m1/2。KI越大，裂纹尖端应力场越强，越接近起裂。

在有限元中提取K的常用方法有三类：

  1. 位移外推法：读取裂纹面附近节点位移，利用裂纹面张开位移与K的关系外推到裂纹尖端。
  2. 应力外推法：读取裂纹前沿前方应力，按r-1/2关系外推。
  3. 相互作用积分法：在裂纹尖端周围建立环形积分域，通过辅助场分离KI、KII、KIII，对网格不敏感，适合复合型裂纹。

三、能量释放率G：从能量角度判断扩展驱动力

能量释放率G定义为单位裂纹面积扩展所释放的弹性能，量纲为能量/面积，常用N/mm或J/m2。在线弹性条件下，G与K存在明确换算关系：

G = (K_I² + K_II²)/E' + K_III²/(2μ)
平面应力：E' = E
平面应变：E' = E/(1-ν²)

其中E为弹性模量，ν为泊松比，μ为剪切模量。该式说明：K描述局部应力场强度，G描述全局能量驱动力，两者在线弹性范围内等价。有限元中可通过虚拟裂纹闭合法（VCCT）直接计算G，也可先算K再换算G。

四、J积分：弹塑性断裂的路径无关积分

当裂纹尖端发生塑性变形时，K主导的奇异场不再适用，此时采用J积分。J积分定义为围绕裂纹尖端的回路积分：

J = ∫_Γ (W dy - T_i ∂u_i/∂x ds)

其中W为应变能密度，Ti为回路上的面力，ui为位移，Γ为围绕裂纹尖端的逆时针回路。J积分具有路径无关性，且在弹塑性条件下可作为裂纹扩展驱动力参数。有限元中J积分通常通过域积分法计算，需在裂纹尖端设置足够细的网格，并保证积分域完全包含塑性区。

五、裂纹尖端场与断裂准则

裂纹尖端场分为三类区域：

  • K主导区：靠近裂纹尖端，应力由K控制，适用于线弹性或小范围屈服。
  • J主导区：弹塑性条件下，由J积分控制，需满足小应变假设。
  • 大范围屈服区：塑性区尺寸与裂纹长度相当，K、J均失效，需采用CTOD或完整弹塑性分析。

常用断裂准则包括：

  • K准则：KI ≥ KIC时起裂，KIC为材料平面应变断裂韧度。
  • G准则：G ≥ GC时起裂，GC为临界能量释放率。
  • J准则：J ≥ JIC时起裂，JIC为弹塑性断裂韧度。
  • 复合型准则：如最大周向应力准则、最小应变能密度准则，用于判断复合型裂纹的扩展方向与起裂条件。
在有限元后处理中，建议同时输出K、G、J三个量并交叉验证：线弹性阶段三者应满足换算关系；若偏差过大，通常说明网格过粗、积分域设置不当或材料模型已进入塑性。
裂纹尖端场分区与断裂参数关系 裂纹尖端局部场 裂纹面 裂纹尖端 K主导区 J主导区 大范围屈服区 x₁ x₂ 断裂参数换算与判据 K 应力强度因子 G 能量释放率 J J积分 线弹性 弹塑性 起裂判据 K ≥ K_IC | G ≥ G_C | J ≥ J_IC 复合型：最大周向应力 / 最小应变能密度准则 有限元中K / G / J的提取与验证流程 ① 建立局部坐标系 裂纹前沿切向/法向 ② 设置积分域 环形域/VCCT ③ 提取K/G/J 相互作用积分/位移外推 ④ 交叉验证与判据 K↔G换算 / 与K_IC比较 验证要点：线弹性阶段 K、G、J 应满足换算关系；若偏差大，检查网格密度、积分域半径与材料模型是否进入塑性。 建议积分域半径取裂纹尖端塑性区尺寸的2~3倍，且不超过裂纹长度的1/10。

六、有限元实操中的关键参数与判断

以I型裂纹为例，假设材料E=210 GPa，ν=0.3，KIC=50 MPa·m1/2，平面应变。若有限元提取得到KI=42 MPa·m1/2，则KI<KIC，裂纹不扩展；若KI=55 MPa·m1/2，则超过KIC，判定起裂。此时可换算G：

E' = E/(1-ν²) = 210000/(1-0.09) ≈ 230769 MPa
G = K_I²/E' = 55²/230769 ≈ 0.0131 N/mm

若采用J积分，弹塑性条件下JIC需由材料试验给出，有限元中J值随载荷增加而增大，当J达到JIC时判定起裂。注意：J积分要求裂纹尖端网格足够细，且积分域应避开塑性区边界，通常取塑性区尺寸的2~3倍。

常见错误：在弹塑性分析中仍用K准则判断起裂，或J积分域设置过小导致路径无关性失效。建议先做网格收敛性检查，再对比不同积分域半径下的J值，若变化小于5%方可认为结果稳定。

综上，裂纹类型决定局部坐标系与K分量分离方式；K、G、J分别从应力场强度、能量释放和弹塑性驱动力三个角度描述裂纹尖端状态；裂纹尖端场分区决定适用哪种参数；断裂准则则将提取结果与材料韧度比较，给出起裂判断。有限元分析中应同时输出K、G、J并交叉验证，确保网格、积分域与材料模型匹配，才能得到可靠的断裂判断。

← 上一章 📖 返回目录 下一章 →