1. 断裂力学导论与课程地图

工程结构在服役中突然断裂，往往并非因为材料强度不足，而是因为结构中存在裂纹或类裂纹缺陷。经典强度设计以“应力小于许用应力”为判据，但当构件含有裂纹时，裂纹尖端的应力会急剧集中，局部应力可远超材料屈服强度，此时强度判据失效。断裂力学正是为回答“含裂纹结构还能不能安全使用、还能用多久”而发展起来的学科分支。

先看两个典型失效场景。其一是焊接压力容器或管道的脆性断裂：结构整体应力水平远低于材料屈服强度，却在低温或冲击载荷下瞬间开裂，断口平齐、无明显塑性变形。其二是飞机蒙皮、桥梁拉索或风电叶片的疲劳裂纹扩展：裂纹在循环载荷下从微小缺陷缓慢长大，达到临界尺寸后失稳断裂。前者关注“裂纹是否立即失稳”，后者关注“裂纹在循环载荷下如何逐周长大”。这两类问题分别对应断裂力学的两大核心：断裂判据与裂纹扩展规律。

为回答第一问，需要引入应力强度因子 K。对于张开型（I型）裂纹，裂纹尖端附近的应力场可写为 σ_ij = K_I / √(2πr) · f_ij(θ)，其中 r 为距裂纹尖端的距离，θ 为极角。K_I 的单位是 MPa·√m，它把外载、裂纹尺寸和构件几何统一到一个参数中。当 K_I 达到材料的断裂韧度 K_IC 时，裂纹失稳扩展。K_IC 是材料抵抗裂纹扩展能力的度量，通常由标准试样（如三点弯曲 SE(B)、紧凑拉伸 CT）测得。

第二问的判据可写成 K_I ≥ K_IC，或等价地用能量释放率 G_I ≥ G_IC。对于线弹性材料，G_I 与 K_I 的关系为 G_I = K_I² / E′，其中平面应力下 E′ = E，平面应变下 E′ = E / (1 − ν²)。这一关系把能量观点与应力强度观点联系起来，也是后续数值仿真中提取断裂参数的基础。

第三问涉及疲劳裂纹扩展。Paris 公式给出 da/dN = C (ΔK)^m，其中 a 为裂纹长度，N 为循环次数，ΔK 为应力强度因子幅值，C 与 m 为材料常数。对金属材料，m 常在 2～4 之间。若已知初始裂纹 a₀、临界裂纹 a_c 和载荷谱，可对上式积分估算剩余寿命。例如某钢制构件 C = 1.5×10⁻¹¹（da/dN 单位 m/cycle，ΔK 单位 MPa·√m），m = 3，a₀ = 2 mm，a_c = 20 mm，ΔK 近似取 15 MPa·√m，则寿命量级可由积分估算，用于判断检修周期是否足够。这里的数据为教学示例，实际工程必须采用材料实测或规范数据。

当裂纹尖端塑性区不可忽略时，线弹性断裂力学不再适用，需要引入弹塑性断裂力学参数，如 J 积分和裂纹尖端张开位移 CTOD。J 积分定义为围绕裂纹尖端的回路积分，具有能量释放率的物理意义；CTOD 则直接描述裂纹尖端的张开程度。它们与 K 的关系在特定条件下可相互换算，但在大范围屈服时需谨慎使用。

本课程的技术对象是“含裂纹的工程结构与材料”，主线是“裂纹—断裂参数—扩展规律—数值仿真—工程判断”。学习路径分为四个层次：第一层建立裂纹类型、应力强度因子、断裂韧度等基本概念；第二层掌握 K、G、J、CTOD 的计算与测试方法；第三层学习疲劳裂纹扩展、环境辅助开裂等扩展规律；第四层进入裂纹扩展仿真，包括有限元建模、裂纹尖端网格、参数提取、扩展路径模拟与结果判读。仿真不是替代理论，而是把理论参数放到真实几何与载荷下求解，因此每一层都需要与前一层对照验证。