# 裂纹学：萌生、扩展与断裂力学

版本：2026-10-28

一条很小的裂纹为什么能改变结构强度？先画清裂纹和受力，再计算尖端驱动力，最后讨论材料阻力、扩展速率和模型是否适用。

这是原创课程阅读稿。自编案例和数值为教学设定；公式的假设写在正文中。网页交互不包含在 Markdown 文件内。部分原文需机构权限，标准以项目采用版次为准。

## 目录

- 01 · 裂纹的语言：长度、深度、前缘与三种加载模式
- 02 · 裂纹从哪里来：滑移、颗粒、晶界与残余应力
- 03 · Griffith 能量平衡：裂纹为什么能够扩展
- 04 · 应力强度因子 K：把载荷、裂纹与几何连起来
- 05 · 塑性区、厚度与小范围屈服
- 06 · J 积分、CTOD 与阻力曲线
- 07 · 韧度试验与稳定性：材料参数怎样得到
- 08 · 疲劳裂纹扩展：ΔK、R 比与 Paris 区间
- 09 · 从扩展速率积分到循环数
- 10 · 短裂纹、闭合与过载：为何真实曲线偏离幂律
- 11 · 环境、时间与混合模式裂纹
- 12 · 裂纹数据怎么分析：从 a–N 到 da/dN–ΔK
- 13 · 数值模拟入门：网格、J 轮廓与裂纹模型
- 14 · 把检测与力学连接：从“发现裂纹”到可审查判断

## 01 · 裂纹的语言：长度、深度、前缘与三种加载模式

先画清几何与坐标，再区分缺口、裂纹、损伤和断裂，避免把不同长度直接代入同一个公式。

### 裂纹不是所有暗线的统称

裂纹是材料内部或表面形成的分离面，尖端或前缘附近具有特殊的应力与变形场。划痕、晶界腐蚀沟、孔洞和夹杂界面都可能在图上呈暗线，需要通过焦深、截面、连续切片或其他检测复核。裂纹闭合时，开口可能很小，影像上不可见不等于完全没有裂纹。

缺口具有有限根部半径，裂纹理想化为尖锐缺陷。缺口应力集中系数 Kt 是无量纲的局部应力与名义应力之比；应力强度因子 K 有 MPa√m 等量纲。二者物理意义不同，不能把 Kt 当成断裂韧度。

### 先约定 a 表示什么

无限板中央贯穿裂纹通常把总长度写成 2a；边裂纹常把从边缘到尖端的长度写为 a；半椭圆表面裂纹常用 a 表示深度、c 表示表面半长。相同的 a 在不同几何中含义不同，代公式前应画图。三维裂纹前缘不同点的 K 也可能不同。

裂纹开口位移、裂纹深度和裂纹投影面积描述不同性质。XCT 得到的张开宽度可能受分辨率限制，超声得到的指示长度也可能包含声束宽度效应。把检测量转为力学模型的 a 前，需要说明几何等效方式和不确定度。

### Mode I、II、III

I 型为张开，裂面沿法向相对分离；II 型为面内滑移，位移方向在裂纹面内且垂直于前缘；III 型为反平面撕裂，位移方向平行于前缘。应以局部裂纹坐标定义，不能仅凭零件整体受拉、受弯或受扭给整个弯曲裂纹贴上一个模式。

混合模式常见。裂纹偏转后，原来的 II 型加载可能在新方向产生张开分量。把 KI、KII、KIII 直接相加没有通用依据，扩展判据与方向预测需要相应模型。

### 学习路径

先学材料力学中的应力、应变、弹性模量、泊松比与屈服，再读 Griffith 能量平衡、Irwin 裂尖场和 Rice 的 J 积分。疲劳、环境扩展和检测结果的工程使用，都建立在清楚的几何与载荷定义上。

练习：中央裂纹总长 4 mm 与边裂纹长 4 mm 的 a 分别为 2 mm 和 4 mm。即使名义应力相同，几何因子也不同，不能仅因图上长度相同就判断裂尖驱动力相等。

### 本章来源

- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [Irwin（1957）：Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate](https://doi.org/10.1115/1.4011547)：应力强度因子与裂尖场。

## 02 · 裂纹从哪里来：滑移、颗粒、晶界与残余应力

理解裂纹萌生与后续扩展的区别，认识组织、表面和制造缺陷如何影响早期损伤。

### 局部变形先于可见裂纹

金属循环加载中，滑移可能集中在部分晶粒和滑移带，表面挤出、侵入及局部应变累积可促成小裂纹。夹杂或第二相破裂、界面脱粘、孔隙边缘及加工缺口也可能成为起点。裂纹不是必然从最大颗粒或最粗晶粒处起裂，还取决于位置、取向和局部载荷。

宏观平均应力较低时，晶粒间各向异性和约束仍会产生局部高应力。拉伸试验的均匀延伸率也不能直接代表缺口根部或界面的局部变形能力。观察到滑移带只是候选证据，需确认它是否与裂纹空间对应。

### 表面为何常见起裂

表面受加工、磨损、腐蚀和残余应力影响，也可能处在弯曲最大拉应力位置。相对内部，表面约束不同、环境接触更直接，因此常见表面起裂。但内部缺陷、夹杂群或某些超高周疲劳条件同样可以导致内部起裂。

表面压缩残余应力可降低某些载荷下的张开驱动力，但可能随循环或温度松弛。不能仅凭一次残余应力测量承诺寿命提高多少，也不能把压应力区域解释为永远不会起裂。

### 萌生寿命的定义依赖检测尺度

“萌生循环数”可能指第一个微观裂纹出现，也可能指裂纹达到 0.1 mm 或仪器可检出的尺寸。不同定义下，萌生与扩展寿命的分配不同。必须报告判定尺寸和检测频率，不能直接比较定义不同的实验。

例如每 10,000 周次检查一次，首次发现裂纹在第 50,000 周次，只能把达到该检出尺寸的时刻约束在前次检查到本次之间，还需考虑前次漏检可能。它不是精确的原子尺度萌生时刻。

### 形成可验证假设

教学假设：表面加工沟槽附近先出现裂纹。可比较沟槽深度、方向、局部应力及起点位置，并查看同批无该沟槽件。若只有一张照片，没有载荷和对照，应保留“沟槽相关应力集中可能参与起裂”的表述。

学习时把起裂原因、扩展机制和最终失稳分成三个问题。它们可能由不同因素控制，不能从最终韧窝形貌推断最早的微裂纹一定由同一机制产生。

### 本章来源

- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASM Handbook, Volume 12: Fractography（2024）](https://doi.org/10.31399/asm.hb.v12.9781627083874)：金属及其他材料断口的宏观、微观形貌与检查方法。
- [Paris 与 Erdogan（1963）：A Critical Analysis of Crack Propagation Laws](https://doi.org/10.1115/1.3656900)：疲劳裂纹扩展速率的经典研究。

## 03 · Griffith 能量平衡：裂纹为什么能够扩展

从总势能变化推导能量释放率 G，理解表面能、有效断裂能和稳定性的区别。

### 把系统的能量账列出来

在准静态弹性条件下，可把结构势能写成 Π=U−W，其中 U 为储存的弹性应变能，W 为外力势的相应项。裂纹面积 A 增大时，结构释放的能量率为 G=−∂Π/∂A，单位 J/m²，也可写 N/m。A 的定义必须与文献约定一致，常以裂纹在原平面的增量面积计。

形成两片新表面需要能量。对理想脆性材料，在采用上述面积约定时临界条件可写 G=2γ，γ 为一片表面的表面能。真实金属还耗散塑性能，不能用纯表面能代表其全部断裂阻力。通常用实验表征的阻力 R 或 Gc。

### 无限板中央裂纹的简化解

对总长 2a 的中央贯穿裂纹，在均匀远场应力 σ、各向同性线弹性假设下，G=πaσ²/E′。平面应力 E′=E，平面应变 E′=E/(1−ν²)。理想脆性起裂条件给出 σc=√(2E′γ/(πa))。这说明缺陷变大时达到扩展所需的远场应力下降。

公式没有包含有限宽度、复杂边界、三维表面裂纹或大范围塑性，不能直接用于任意零件。真实材料的 R 还可能随扩展量变化；“达到 G=R”只说明可能开始扩展，后续是否稳定要比较驱动力和阻力随 a 的变化。

### 手算与单位

教学设定 E′=70 GPa，γ=1 J/m²，a=0.1 mm。代入 SI 单位得到 σc≈21.1 MPa。若 a 增为 0.4 mm，σc 减半约为 10.6 MPa。这些数值不对应某种已认证玻璃，是用来检查平方根关系与单位的算例。

常见错误是把 0.1 mm 直接当作 0.1 m 代入，导致应力差 √1000 倍；或把总裂纹长度 2a 误当成 a。每次计算先把几何图、单位和变量表写在公式旁。

### 从历史论文到实际材料

Griffith 1921 年工作的价值在于用能量描述含缺陷固体的断裂，而不是提供所有材料的通用强度常数。对于有塑性耗散的金属，需要进一步学习 Irwin 的 K 与弹塑性断裂参量。

练习：若弹性结构在固定力和固定位移下扩展，外部系统做功不同，G 随裂纹长度的变化可能不同。由此可见，稳定性不能只看材料，也与加载装置和边界条件有关。

### 本章来源

- [Griffith（1921）：The phenomena of rupture and flow in solids](https://doi.org/10.1098/rsta.1921.0006)：裂纹能量平衡的经典起点。
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。

## 04 · 应力强度因子 K：把载荷、裂纹与几何连起来

解释裂尖奇异场、几何因子和 K 的单位，并用可核算例区分载荷效应与材料韧度。

### 近尖端为什么出现平方根

Irwin 的线弹性裂尖场可写为 σij≈K/√(2πr)·fij(θ)，其中 r、θ 是以尖端为原点的局部极坐标，fij 随加载模式而异。r 趋近零时公式发散，意味着理想尖裂纹的线弹性解失效，不代表真实材料能承受无限应力。实际尖端附近会发生塑性、钝化或其他损伤。

K 表示这一近场的强度。它由载荷、几何及裂纹尺寸决定；KIc 则是在规定条件和有效试验下得到的材料断裂韧度。不能将某次零件计算得到的 KI 称为材料 KIc。

### 常见形式与几何因子

I 型常写 KI=Yσ√(πa)。无限板中央贯穿裂纹 Y=1；半无限板边缘小裂纹的经典近似约为 1.12，但有限宽度、裂纹较长、弯曲、孔边和表面裂纹都需要相应解。Y 有时随 a/W 和前缘角变化，不是一个随意选择的安全系数。

若 σ 用 MPa、a 用 m，K 的单位是 MPa√m。若 a 用 mm，则结果是 MPa√mm；1 MPa√m=√1000 MPa√mm。混用会造成约 31.62 倍的误差。

### 手算：裂纹增大后的变化

教学设定无限板中央裂纹半长 a=2 mm，远场 σ=100 MPa，Y=1，则 KI≈7.93 MPa√m。a 增为 8 mm，K 变为约 15.85 MPa√m；应力加倍也使 K 加倍。平方根关系和线性关系应分清。

若假定同条件的临界韧度为 30 MPa√m，理想无限板给出 ac=(30/100)²/π≈28.65 mm。这个 ac 是半长，且只有当无限板近似、线弹性及对应韧度适用时才有意义；它不是任何真实零件的允许裂纹尺寸。

### 计算前的四项检查

检查 a 的定义、Y 的适用几何、应力来自何处、韧度是否适用于当前厚度/温度/取向/速率。残余应力可能沿厚度不均匀，不能总用一个常数叠加到名义应力。线弹性叠加虽可用，但要先求每项载荷对应的 K。

本页互动只采用无限板中央裂纹的 Y=1 模型，并公开参数与假设。改变 a 和 σ 可观察 K 增大，但页面不输出工程合格判定。

### 本章来源

- [Irwin（1957）：Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate](https://doi.org/10.1115/1.4011547)：应力强度因子与裂尖场。
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。

## 05 · 塑性区、厚度与小范围屈服

用塑性区数量级理解线弹性断裂力学的适用范围，区分平面应力和高约束状态。

### 尖端塑性区如何估算

用线弹性裂尖正应力达到屈服强度的近似，可得到 I 型平面应力下 rp≈(K/σy)²/(2π)。常用的平面应变估算为 rp≈(K/σy)²/(6π)。这些是数量级公式，实际塑性区形状受角度、屈服准则、硬化和约束影响，不是完美圆盘。

例如教学参数 K=20 MPa√m、σy=400 MPa，平面应力 rp≈0.398 mm，平面应变估算约 0.133 mm。K 与 σy 使用相同应力单位，平方后自然得到长度 m。

### 什么是小范围屈服

小范围屈服要求塑性区相对裂纹尺度、剩余韧带、厚度和其他控制几何足够小，使外部仍存在由 K 控制的区域。不能只把 rp 与一个尺寸比较，也不能认为出现任何塑性就完全不能使用 K。

如果裂纹长度仅 0.2 mm，而估算塑性区约 0.4 mm，简单 K 主导假设就很可疑。此时需要弹塑性分析或针对短裂纹的处理。一个很小的 r 点落在塑性区内，不能继续使用发散的弹性应力解释真实应力值。

### 厚度改变约束

薄板表面附近更接近平面应力，厚板内部在一定条件下可接近平面应变。约束增加通常限制塑性侧向变形，并影响测得的断裂阻力。实际三维裂纹前缘上约束可以变化，厚度本身不是唯一判据。

有效的平面应变韧度 KIc 需要满足相应试样、预裂、载荷及尺寸有效性条件。试验所得 KQ 是候选值；未经核查不能改名为 KIc。具体数值判据应按采用版 ASTM E399 查证，而非只用本章的 rp 估算替代。

### 把近似用在正确位置

塑性区估算适合做前期适用性检查，也帮助理解为什么薄片可能比厚件表现出更大的表观断裂阻力。但它不是完整试验规程，更不能直接用来宣布一个构件安全。

练习：屈服强度加倍、K 不变时，rp 变为原来的四分之一；K 加倍则 rp 增四倍。理解这个平方关系后，再判断某些高强材料虽然塑性区小，却可能对尖锐缺陷更敏感。

### 本章来源

- [Irwin（1957）：Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate](https://doi.org/10.1115/1.4011547)：应力强度因子与裂尖场。
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASTM E399](https://store.astm.org/standards/e399)：金属线弹性平面应变断裂韧度试验

## 06 · J 积分、CTOD 与阻力曲线

从能量和裂尖张开认识弹塑性断裂，分清参量定义、起裂值和稳定扩展阻力。

### J 积分的定义与条件

二维形式可写 J=∫Γ(W nₓ−Ti ∂ui/∂x)ds。W 是应变能密度，Ti 是路径上的面力，ui 为位移，x 为裂纹扩展方向，Γ 环绕裂尖。对于相应弹性或形变理论条件，J 可与能量释放率联系并具有路径无关性。

Rice 1968 年论文的条件很重要。存在卸载、复杂循环塑性、材料非均匀、体力或其他源项时，简单形式可能不再直接路径无关，需要扩展定义和数值处理。有限元中几条轮廓的 J 相近只是检查之一，不自动证明整个模型正确。

### 与 K 的联系

在线弹性 I 型情形，J=G=KI²/E′。若 K=20 MPa√m、E′=200 GPa，则 J=2000 J/m²=2 kJ/m²。这里把 MPa 换为 Pa 后平方计算，避免出现百万倍单位错误。混合模式各分量的贡献要采用对应公式。

进入显著塑性后，J 仍常用于适当条件下的裂尖表征，但不能用弹性 K 公式取代全部塑性功。实验中 J 的弹性与塑性部分如何估算，取决于试样、载荷位移曲线和标准程序。

### CTOD 与 J–R 曲线

CTOD 是裂纹尖端张开位移，单位为长度。其具体测量或构造定义应注明。常见量级关系 δ≈J/(mσflow) 中，m 与约束及模型有关，σflow 的定义也需约定，因此它不是无需校准的通用精确换算。

J–R 曲线记录裂纹增量 Δa 与断裂阻力的关系。起裂韧度、某一扩展量的阻力和最终失稳点不是同一个量。裂纹可以先稳定扩展，再在驱动力增长超过阻力增长时失稳。

### 实验记录的必要信息

保存试样方向、厚度、预裂长度、载荷—位移数据、裂纹扩展测量和卸载顺应性等信息。试样边沟、约束、温度和加载速率可能影响结果。标准中的资格判定需要逐项执行，不能从一条漂亮的曲线直接报告有效 JIc。

练习：把 J=2 kJ/m² 写为 2 N/mm；两者数值相等。再说明为何 CTOD 不能写成 kJ/m²：它是位移，不是单位面积能量。

### 本章来源

- [Rice（1968）：A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks](https://doi.org/10.1115/1.3601206)：J 积分及其路径无关条件。
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASTM E1820](https://store.astm.org/standards/e1820)：金属断裂韧度测定，涉及 J、CTOD 等

## 07 · 韧度试验与稳定性：材料参数怎样得到

理解疲劳预裂、载荷位移曲线、裂纹测量与有效性条件，避免把冲击功和断裂韧度混同。

### 为何需要尖锐预裂纹

机加工缺口具有有限根部半径，不能完全代表尖裂纹。断裂力学试样常通过疲劳预裂形成更接近实际尖端的裂纹，再进行后续加载。预裂载荷过高可能产生较大塑性区并影响测量，所以其过程也受规程约束。

紧凑拉伸 C(T) 和单边缺口弯曲 SE(B) 试样的 K、J 计算包含特定几何函数。不要把无限板的 Y=1 代入这些试样。试验前后裂纹前缘可能弯曲，需要按标准指定的测点与方式测量。

### 载荷曲线不是全部证据

应记录力、位移和必要的裂纹长度信息，检查机器柔度、传感器标定、对中和温度。偏载、夹具间隙、塑性变形和裂纹跳跃都可能改变曲线。只截一段最大载荷附近曲线会丢失试验有效性证据。

ASTM E399 主要针对金属线弹性平面应变韧度，E1820 涉及更广的韧度测定。具体采用哪一个，由材料行为、研究问题与项目要求决定；不能先得到某个数字，再挑一个能让它“合格”的标准。

### 冲击功为什么不能直接变成 KIc

Charpy 冲击功是特定缺口试样和动态条件下吸收的能量，包含多种贡献。KIc 是有效条件下的应力强度临界值，两者的试样、缺口、加载速率和量纲不同。某些材料存在经验相关式，但适用范围有限，不能作为普适换算。

同样，抗拉强度高不代表断裂韧度高。选材需分别考虑强度、韧度、环境和缺陷容限。材料手册中的数值还应核对取向、温度、厚度及统计性质。

### 稳定性要结合结构

在阻力随裂纹扩展上升的材料中，G=R 可以对应稳定扩展。若继续扩展时驱动力增加更快，可能失稳。载荷控制与位移控制、结构柔度和机器柔度都会影响这场竞争。

教学练习：两块材料拥有相同的起裂 J，却有不同的 J–R 曲线斜率。它们开始扩展的条件可能相近，随后抵抗稳定撕裂的能力却不同；只列一个起裂值无法完整描述全过程。

### 本章来源

- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASTM E399](https://store.astm.org/standards/e399)：金属线弹性平面应变断裂韧度试验
- [ASTM E1820](https://store.astm.org/standards/e1820)：金属断裂韧度测定，涉及 J、CTOD 等

## 08 · 疲劳裂纹扩展：ΔK、R 比与 Paris 区间

建立从载荷范围到扩展速率的计算过程，说明门槛区、中间区和快速扩展区不能用同一条直线代替。

### 先定义循环载荷

Kmax 与 Kmin 是同一循环的最大、最小应力强度因子，ΔK=Kmax−Kmin，R=Kmin/Kmax。在几何与线弹性关系相同且没有额外复杂项时，R 可对应载荷比。若裂面发生闭合或压缩阶段，名义 ΔK 与实际有效张开范围要分开。

相同 ΔK、不同 R 的试验可能具有不同扩展速率，因为 Kmax、闭合和环境过程都可能变化。报告扩展数据时不能只写一个 ΔK，至少还需载荷比、波形、频率、环境、温度及材料方向。

### Paris 幂律的用途

常用关系 da/dN=C(ΔK)^m 描述某些材料长裂纹的中间扩展区。da/dN 是每循环裂纹增长量，C 的数值与采用单位有关，m 为无量纲指数。门槛附近和接近失稳时通常偏离该幂律。

经典论文 Paris 与 Erdogan（1963）强调检验裂纹扩展规律的证据。拟合时若混入门槛区、过载后的瞬态或最终失稳数据，单一直线可能掩盖机制变化。不能把某材料手册的 C、m 用到另一热处理、R 比或环境。

### 一个可核算的速率例子

教学设定 C=1×10⁻¹¹，配套单位为 da/dN 使用 m/cycle、ΔK 使用 MPa√m，m=3。在 ΔK=10 MPa√m 时，速率为 1×10⁻⁸ m/cycle=0.01 μm/cycle；ΔK 加倍，速率增为八倍。C 只是演算参数，不代表某一合金。

若把速率单位改成 mm/cycle，C 的数值应随之改变。若 ΔK 单位也改变，C 还要按幂次换算。抄录一个没有单位的 C 数值，是寿命估算中很危险的错误来源。

### 怎样理解门槛

ΔKth 与试验定义、R 比、环境和加载历史有关，不能理解为一切尺寸裂纹的永不扩展界限。短裂纹可能在长裂纹门槛以下扩展，腐蚀环境也会改变行为。检测不到裂纹与裂纹不扩展属于两个不同判断。

对数据拟合应保存原始 a–N、载荷历史和取点方法。只提供拟合系数而没有有效区间和数据散布，会使后续使用者无法判断该模型是否适合。

### 本章来源

- [Paris 与 Erdogan（1963）：A Critical Analysis of Crack Propagation Laws](https://doi.org/10.1115/1.3656900)：疲劳裂纹扩展速率的经典研究。
- [Elber（1970）：Fatigue crack closure under cyclic tension](https://doi.org/10.1016/0013-7944(70)90028-7)：拉伸循环中的裂纹闭合及断口证据。
- [ASTM E647](https://store.astm.org/standards/e647)：疲劳裂纹扩展速率试验

## 09 · 从扩展速率积分到循环数

逐步计算理想恒幅条件下的扩展循环数，并说明初始裂纹、终止条件和变幅载荷的影响。

### 积分前先设定边界

若 ΔK=YΔσ√(πa)，且 Y、Δσ、C、m 恒定，则 N=∫a₀ᵃᶠ da/[C(YΔσ√π)^m a^(m/2)]。a₀ 是开始计算的裂纹尺寸，af 是计算终点，二者都必须落在模型的有效区间。N 只包含这段扩展，不含此前萌生。

对 m≠2，N=[af^(1−m/2)−a₀^(1−m/2)]/[(1−m/2)C(YΔσ√π)^m]；m=2 时为 ln(af/a₀)/[C(YΔσ√π)²]。这两个表达在边界一致时为正，若算出负寿命，应检查积分上下限和符号。

### m=3 的手算

教学参数：Y=1、Δσ=100 MPa、C=10⁻¹¹、m=3，a 从 1 mm 增至 4 mm。令 A=C(Δσ√π)³≈5.5683×10⁻⁵，则 N=2/A·(a₀⁻¹ᐟ²−af⁻¹ᐟ²)≈5.679×10⁵ cycles。a 使用 m，C 与上一章的配套单位相同。

这只是理想幂律积分。要将它用于材料，必须确认对应 ΔK 区间、R 比、频率和环境有实测依据，并确认 Kmax 没有先达到失稳条件。演算不能证明任何零件还有约 57 万次寿命。

### 为什么初始尺寸很重要

对 m>2 的理想模型，小裂纹阶段对积分有较大贡献，但短裂纹又可能偏离长裂纹 Paris 规律。初始裂纹取自检测时，还受到检出概率和尺寸误差影响。把检测下限当作“所有未检出裂纹的最大尺寸”并无统计依据。

实际有限宽板的 Y 常随 a 改变，应数值积分并更新几何。若存在多个裂纹合并、分叉或非半椭圆前缘，单变量 a 模型可能不再合适。

### 变幅不能只取平均应力

过载可通过塑性与闭合影响后续扩展，欠载和环境暴露也可能产生顺序效应。两段含同样循环数但排序不同的载荷谱，不保证有相同扩展。把全部载荷压成一个平均 Δσ，因幂律非线性也会丢失信息。

教学检查：若其他条件固定且始终在同一 m=3 区间，Δσ 加倍使积分循环数变为八分之一。这个比例关系比背诵单一寿命数更值得掌握。

### 本章来源

- [Paris 与 Erdogan（1963）：A Critical Analysis of Crack Propagation Laws](https://doi.org/10.1115/1.3656900)：疲劳裂纹扩展速率的经典研究。
- [ASTM E647](https://store.astm.org/standards/e647)：疲劳裂纹扩展速率试验
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。

## 10 · 短裂纹、闭合与过载：为何真实曲线偏离幂律

理解微观尺度和加载历史对裂纹扩展的影响，避免把长裂纹模型外推到萌生阶段。

### “短”有多种含义

微观短裂纹的长度可与晶粒或组织尺度相当；力学短裂纹相对塑性区或局部应力梯度较小；物理短裂纹还常用某种约定尺寸描述。这些定义不是同一个固定毫米数。短裂纹扩展会受到晶界、滑移面、表面及局部应力影响。

一个短裂纹可能先快后慢，在晶界停滞后再次扩展。将这种曲线与均匀材料的长裂纹 Paris 曲线直接比较，可能误以为仪器出了问题。需要同时记录组织位置与裂纹路径。

### 裂纹闭合的意义

裂面接触可能来自塑性尾迹、粗糙度或氧化产物等。常用 ΔKeff=Kmax−Kop 表示完全打开后的一部分驱动力范围，其中 Kop 为开口对应的应力强度；具体定义和测量方法需要说明，也须处理 Kop 与 Kmin 的关系。

Elber 的经典结果说明拉伸循环中也可以闭合。闭合修正可能改善某些数据归一化，但不消除所有 R 比、环境或短裂纹影响。不同开口测量方法可给出不同结果，不能把 Kop 当作没有误差的常数。

### 过载为什么可能产生迟滞

一次较大的拉伸载荷可能形成较大的塑性区，改变残余应力和闭合状态，使后续裂纹暂时减速；结果依赖材料、厚度、过载比、基线载荷和后续历史。不能把“过载必然延寿”当作普遍规则，过载本身也可能直接造成危险扩展。

建模需要相应迟滞或闭合模型及验证数据。只调小 Paris 系数来适配一个过载试验，通常不能解释其他载荷序列。

### 如何读取实验曲线

将 a–N 曲线、载荷谱、局部组织与测量分辨率放在同一时间线上。区分真实停滞与仪器在一段循环内尚不能分辨微小变化。差分会放大噪声，少量抖动可能使瞬时 da/dN 看似为负。

报告可同时给原始曲线、采用的平滑或增量多项式方法、窗口大小及其敏感性。短裂纹研究尤其需要说明裂纹长度是表面投影、真实弧长还是等效尺寸。

### 本章来源

- [Elber（1970）：Fatigue crack closure under cyclic tension](https://doi.org/10.1016/0013-7944(70)90028-7)：拉伸循环中的裂纹闭合及断口证据。
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASTM E647](https://store.astm.org/standards/e647)：疲劳裂纹扩展速率试验

## 11 · 环境、时间与混合模式裂纹

区分按循环和按时间描述的扩展，理解应力腐蚀、蠕变与混合模式所需的附加信息。

### 循环数并非唯一时间变量

疲劳常用 da/dN，应力腐蚀或持续载荷下的扩展常用 da/dt；两者单位不同。若环境过程在保载时仍继续，增加每循环保载时间可能显著改变每循环增长量。把不同频率的数据只按循环数归一，会隐藏时间效应。

某些工况下用疲劳项与时间项相加是工程近似，但相互作用可能使简单加法失效。需要对目标载荷波形、温度和介质验证，不能任意拼接两份不同来源的速率公式。

### 环境阈值具有条件性

应力腐蚀扩展需要特定材料/环境/拉应力组合，相关阈值或速率曲线依赖试验条件。腐蚀疲劳、氢相关开裂和高温蠕变还涉及不同过程。断口形貌可帮助判断，但不足以单独确定速率模型。

记录介质、溶解组分、温度、电化学条件、加载波形、频率与保载时长。材料状态、残余应力及取样储存也重要。一个没有环境信息的“裂纹扩展速率”数值通常无法复用。

### 混合模式与扩展方向

当 KI、KII、KIII 同时存在时，可使用最大周向应力、能量释放率或其他准则预测起裂方向和条件，但它们有各自假设。材料各向异性、界面裂纹、塑性和动态效应可能要求更专门的模型。

三维前缘局部模式可能不同；裂纹扭转后，面上摩擦接触与约束改变。不能根据载荷方向在二维草图中简单把所有模式化成一个“等效 K”，除非说明采用的准则和验证范围。

### 工程练习

设焊接件在含氯环境中承受循环载荷，并在每周期峰值处保载。应首先建立循环与时间双重记录，核对焊接残余应力、材料敏感性及真实裂纹路径。用空气、无保载的 Paris 数据直接预测，缺少关键条件。

本课给出分析变量与推理顺序，不给未验证的环境修正倍数。若没有适用扩展数据，应报告模型缺口并设计对照试验，而不是用一个看似保守的系数替代证据。

### 本章来源

- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASM Handbook, Volume 12: Fractography（2024）](https://doi.org/10.31399/asm.hb.v12.9781627083874)：金属及其他材料断口的宏观、微观形貌与检查方法。
- [ASTM E647](https://store.astm.org/standards/e647)：疲劳裂纹扩展速率试验

## 12 · 裂纹数据怎么分析：从 a–N 到 da/dN–ΔK

给出原始数据校验、差分、拟合与不确定度检查的顺序，并解释与断口测量的区别。

### 保存原始测量与校准

裂纹长度可通过光学、柔度、电位法或其他技术估计，各有校准关系和检测尺度。记录原始信号、换算函数、标尺、载荷同步及时间/循环索引。试后断面检查可用于复核实际前缘形状，但不能无说明地改写全部原始数据。

测量表至少含 N、a、Pmax、Pmin、温度、频率和异常标记。N 应递增，单位一致；仪器量程改变、暂停和重新对中应作为事件保存。异常值先查原因，不能因为它破坏直线拟合就删除。

### 手算割线速率

教学数据：N=0、10,000、20,000 cycles 时，a=1.00、1.12、1.30 mm。两个区间割线速率分别为 1.2×10⁻⁸ 和 1.8×10⁻⁸ m/cycle。可把速率与区间中点长度配对作教学近似，但正式分析应遵循选定标准的方法。

如果两次长度读数各有约 0.01 mm 的独立标准不确定度，则差值的不确定度约为 √2×0.01 mm。较短测量区间虽提供更多点，却可能使长度增量接近噪声，从而放大速率误差。

### 计算 ΔK 与拟合

对每个有效区间按实际试样几何求 ΔK，核对 a/W 范围。对中间稳定区用 log(da/dN)=log C+m log ΔK 可作线性拟合，但需关注坐标变换后的误差结构；并非所有数据都适合等权最小二乘。

应展示散点、拟合区间、样本数、残差和适当的不确定度。不同试样、载荷比或环境的数据先分别标识。重复测量同一个试样的许多点，不能当作许多独立材料批次。

### 与断口条带和检测量对照

条带间距是特定位置的局部断口量，a–N 常是某种宏观裂纹长度的演化。它们可能在适用条件下相互支持，但存在前缘弯曲、投影和不可见条带区等差异。不能要求两种量在所有位置完全相等。

XCT 或超声测到的裂纹尺寸也应保留误差范围。若用它们更新增长模型，需说明测量几何如何转成模型 a，以及未检出结果如何进入概率分析。

### 本章来源

- [ASTM E647](https://store.astm.org/standards/e647)：疲劳裂纹扩展速率试验
- [Paris 与 Erdogan（1963）：A Critical Analysis of Crack Propagation Laws](https://doi.org/10.1115/1.3656900)：疲劳裂纹扩展速率的经典研究。
- [ASM Handbook, Volume 12: Fractography（2024）](https://doi.org/10.31399/asm.hb.v12.9781627083874)：金属及其他材料断口的宏观、微观形貌与检查方法。

## 13 · 数值模拟入门：网格、J 轮廓与裂纹模型

建立有限元断裂分析的基本检查顺序，区分裂尖参量计算、损伤演化模型和图像展示。

### 模型必须先回答明确问题

若只需求一个已知裂纹的 KI，可用适用解析解、权函数或有限元。若需模拟稳定撕裂、界面脱粘或路径选择，可能采用黏聚区、XFEM、相场或其他模型。模型复杂度增加并不会自动减少参数不确定度。

几何、边界、载荷、材料本构和初始裂纹缺一不可。输入实际测量轮廓前要决定哪些细节影响研究问题；过度平滑会改变局部曲率，逐体素阶梯边界也可能制造虚假应力集中。

### 线弹性裂尖计算的检查

先用有解析解的简单几何验证单位、载荷和 K 的提取。检查网格收敛、裂尖单元处理、积分轮廓位置及边界距离。多个轮廓结果一致可增强信心，但错误边界条件同样可能产生一致结果。

不要把裂尖最大节点应力当成网格无关的材料参数。理想尖裂纹应力随网格细化继续增大；应比较 K、J 或合适的平均量。若塑性显著，还要确认所求参量在当前模型中的意义。

### 演化模型需要额外参数

黏聚模型需要牵引—分离关系及断裂能等参数；相场模型通常涉及能量、长度尺度和分解选择。参数可能需要独立标定。仅让一条载荷位移曲线贴合，不能证明裂纹路径与尺寸效应都正确。

模型输出的损伤变量不一定就是可直接测量的裂纹开口。将模拟结果和 XCT、DIC 或超声图比较时，需说明测量分辨率与模型量之间的映射。不能把彩色损伤图当作实测断口。

### 报告最小可复现集

保存模型版本、单位、材料来源、载荷历程、几何、网格和收敛检查。分别列出已知输入、标定参数与独立验证数据。用另一几何或载荷条件验证，比继续增加同一拟合曲线的小数位更有说服力。

本网站动画采用解析表达与示意运动，没有求解真实弹塑性裂纹扩展。它用于理解变量关系，不能代替有限元验证或构件试验。

### 本章来源

- [Rice（1968）：A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks](https://doi.org/10.1115/1.3601206)：J 积分及其路径无关条件。
- [Irwin（1957）：Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate](https://doi.org/10.1115/1.4011547)：应力强度因子与裂尖场。
- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。

## 14 · 把检测与力学连接：从“发现裂纹”到可审查判断

学习如何整理几何、载荷、材料阻力和检测不确定度，形成完整的裂纹评估输入表。

### 四组输入缺一不可

第一组是裂纹几何：位置、深度、表面长度、取向、前缘形状及其误差。第二组是载荷：膜应力、弯曲、残余应力、温度、载荷谱。第三组是材料阻力：适用韧度、屈服强度、扩展规律及环境条件。第四组是检测能力：覆盖、盲区、定量误差和检出概率。

一个 XCT 孔隙团并不总能直接当成平面裂纹。一个超声高幅指示也不能直接赋值为真实裂纹面积。应先确认缺陷类别与几何模型的合理性，再讨论等效化方式。

### 教学计算与实际评估的差别

课堂可用 Y=1、恒定应力和给定韧度说明 K 随 a 增长。实际评估可能涉及有限几何、塑性坍塌、裂纹相互作用、残余应力和断裂评定图。选择评价方法应遵循项目采用的规范及其适用范围。

不能把计算 K 小于某个手册 KIc 就当成全工况合格。温度、加载速率、材料方向、韧度散布和极端载荷都可能改变结论。安全裕度的定义也应由相应工程流程给出。

### 检出能力与复检间隔

检出概率 POD 随缺陷大小、位置、方向和操作条件变化。a90/95 常用于描述在给定统计约定下，POD 的 95% 下置信界达到 90% 时对应的缺陷尺寸；这并不是绝对检出下限。其数值不能从一张漂亮的扫描图得出。

复检安排通常需要裂纹增长分布、检测能力和可接受风险共同支持。教学积分的单一寿命数字没有覆盖这些不确定性，因此不应自动转成检查周期。

### 完成一份学习记录

最后写出“已知什么、计算了什么、假设什么、还缺什么”。例如：检测给出表面裂纹投影长度，但深度未可靠测得，那么下一步应优先补充深度信息或使用有依据的边界分析，而不是继续精修一个随意假设的 a。

四个专题在此汇合：断口分析提供失效后的机制证据，裂纹力学解释驱动力，XCT 和超声提供特定分辨率与覆盖条件下的几何及指示。它们相互补充，各自的误差和适用条件仍需保留。

### 本章来源

- [Anderson：Fracture Mechanics，第 4 版（2017）](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)：线弹性、弹塑性断裂力学与试验应用。
- [ASTM E399](https://store.astm.org/standards/e399)：金属线弹性平面应变断裂韧度试验
- [ASTM E1820](https://store.astm.org/standards/e1820)：金属断裂韧度测定，涉及 J、CTOD 等
- [NIST：XCT 缺陷检出概率的实验与模拟（2021）](https://doi.org/10.1016/j.ndteint.2021.102416)：将参考尺寸不确定度、阈值、位置与取向纳入检出分析。

## 分析工作流

### 裂纹力学工作流：几何、驱动力与适用性

结构驱动力与材料阻力在对应条件下比较；模型选择先于计算。

1. 画几何：标出 a、2a、c、厚度、宽度、前缘与载荷方向；注明检测量如何转换。

2. 整理载荷：区分膜、弯曲、残余应力、温度和循环谱，保留来源与不确定度。

3. 选模型：检查几何函数和 a/W 范围；用塑性区、剩余韧带和约束判断 K 主导假设。

4. 计算并手验：用 SI 或明确配套单位，先以一个可解析极限验证代码与量纲。

5. 核对材料：韧度及扩展参数对应材料状态、方向、温度、速率、环境与载荷比。

6. 敏感性分析：对尺寸、载荷和材料散布分别扰动，避免隐藏主导不确定因素。

7. 形成计算书：列出输入、公式、来源、适用范围与未覆盖工况，不把教学值用作许可值。

解释限制：线弹性近似不适合任意大塑性或所有短裂纹。；构件合格性需相应工程评价流程，单个 K 比值不足以涵盖全部失效方式。

### 裂纹扩展数据工作流：a–N、速率与拟合

扩展速率由长度增量和循环增量获得；ΔK 按实际试样几何与载荷计算。

1. 准备原始表：N、a、Pmax/Pmin、频率、温度、环境和异常标记；保留校准信号。

2. 数据核查：处理单位、重复索引和仪器事件；不因偏离预期直线而直接删除点。

3. 选择速率估计：按标准或有依据方法进行割线/局部拟合，记录窗口和长度分辨率。

4. 逐区间计算 ΔK 与 R，检查试样几何函数、预裂及塑性有效性。

5. 分区拟合：仅在支持幂律的区间求 C、m，报告单位、散布和残差。

6. 验证模型：独立试样或不同条件作预测检查；积分时更新几何并设置适用终点。

解释限制：差分会放大长度噪声。；长裂纹拟合不能无验证外推至萌生、短裂纹、变幅或新环境。

## 术语

- **a / c**：裂纹尺度，具体是半长、深度或其他定义由几何解规定。
- **K**：裂尖弹性场强度；常用 MPa√m。
- **G / J**：单位裂纹面积的能量参量；常用 J/m² 或 N/mm。
- **CTOD**：裂纹尖端张开位移，需注明定义。
- **R 比**：Kmin/Kmax 或相应载荷比，按条件区分。
- **da/dN**：单位循环的裂纹增量。
- **小范围屈服**：塑性区相对控制几何小，外部存在适用的 K 主导区域。

## 自测与解析

### 中央裂纹总长 4 mm，公式中 a 是多少？

在本课无限板中央裂纹约定中 a=2 mm；必须先核对所用几何解。

### Y 和应力不变，a 增四倍，K 怎么变？

K 与 √a 成正比，因此增两倍。

### K=20 MPa√m，E′=200 GPa，线弹性 J 是多少？

J=K²/E′=2000 J/m²=2 kJ/m²，先把应力单位统一。

### 测得一个 KQ，可以直接写成 KIc 吗？

不能。应核查试样、预裂、载荷和尺寸等有效性条件。

### Paris 积分给出的是全部寿命吗？

只给所选 a₀ 到 af 的模型扩展循环数，不包含此前萌生，且须满足模型区间。

### 长裂纹门槛以下，小裂纹一定停止吗？

不一定，微观短裂纹、闭合和环境等因素可能使长裂纹规律不适用。

## 教材、论文与官方文档

### ASM Handbook, Volume 12: Fractography（2024）

Schroeder, C. J., Parrington, R. J., Maciejewski, J. O. & Lane, J. F. (eds.). ASM International, 2024.

[原始来源](https://doi.org/10.31399/asm.hb.v12.9781627083874)

先读 Types and Appearances of Fracture 与 Fracture Examination，再按材料查阅。1987 版与 2024 版编排不同；本课程不虚构对应页码。出版社目录可公开访问，章节可能需机构权限。

### Anderson：Fracture Mechanics，第 4 版（2017）

Anderson, T. L. Fracture Mechanics: Fundamentals and Applications. CRC Press, 2017.

[原始来源](https://www.routledge.com/Fracture-Mechanics-Fundamentals-and-Applications/Anderson/p/book/9781315370293)

先读能量平衡与裂尖场，再读材料行为、韧度试验和疲劳扩展；做题时先画清裂纹几何与载荷。公开书目用于核对版本，正文讲解是独立组织的课程。

### Griffith（1921）：The phenomena of rupture and flow in solids

Griffith, A. A. Philosophical Transactions of the Royal Society A 221, 163–198 (1921).

[原始来源](https://doi.org/10.1098/rsta.1921.0006)

抓住弹性能释放与新生表面能的竞争。玻璃纤维研究说明缺陷尺寸与强度有关；不要把理想脆性表面能模型原样当作金属韧度模型。

### Irwin（1957）：Analysis of Stresses and Strains Near the End of a Crack Traversing a Plate

Irwin, G. R. Journal of Applied Mechanics 24, 361–364 (1957).

[原始来源](https://doi.org/10.1115/1.4011547)

注意近场表达的适用区域、平面应力与平面应变的区别。K 的量纲是应力乘长度平方根。

### Paris 与 Erdogan（1963）：A Critical Analysis of Crack Propagation Laws

Paris, P. & Erdogan, F. Journal of Basic Engineering 85, 528–533 (1963).

[原始来源](https://doi.org/10.1115/1.3656900)

重点是如何检验扩展规律，而不只是记住幂律。C、m 必须连同单位、R 比、环境与拟合区间报告。

### Rice（1968）：A Path Independent Integral and the Approximate Analysis of Strain Concentration by Notches and Cracks

Rice, J. R. Journal of Applied Mechanics 35, 379–386 (1968).

[原始来源](https://doi.org/10.1115/1.3601206)

原论文讨论弹性及形变理论弹塑性材料。卸载、循环塑性、非比例加载或非均匀材料中，不能直接照搬最简单的路径无关结论。作者存档：esag.harvard.edu/rice/015_Rice_PathIndepInt_JAM68.pdf。

### Elber（1970）：Fatigue crack closure under cyclic tension

Elber, W. Engineering Fracture Mechanics 2, 37–45 (1970).

[原始来源](https://doi.org/10.1016/0013-7944(70)90028-7)

裂纹可能在外载仍为拉力时接触闭合。读懂开口载荷，再理解名义 ΔK 与有效驱动力为何不必相同。

### NIST：XCT 缺陷检出概率的实验与模拟（2021）

Merging experiments and computer simulations in X-ray Computed Tomography probability of detection analysis of additive manufacturing flaws. NDT & E International (2021).

[原始来源](https://doi.org/10.1016/j.ndteint.2021.102416)

先理解检出率随缺陷尺寸变化，再读统计模型。一次扫描没有发现指示，不能据此给出 100% 无缺陷的结论。NIST 官方书目已核对。

### ASTM E399

ASTM E399，以项目采用的有效版次为准。

[原始来源](https://store.astm.org/standards/e399)

这里引用标准用途与官方入口，不转录付费条款，也不提供未经核实的验收限值。正式试验须核对适用材料、试样、校准、有效性条件和采用版次。

### ASTM E1820

ASTM E1820，以项目采用的有效版次为准。

[原始来源](https://store.astm.org/standards/e1820)

这里引用标准用途与官方入口，不转录付费条款，也不提供未经核实的验收限值。正式试验须核对适用材料、试样、校准、有效性条件和采用版次。

### ASTM E647

ASTM E647，以项目采用的有效版次为准。

[原始来源](https://store.astm.org/standards/e647)

这里引用标准用途与官方入口，不转录付费条款，也不提供未经核实的验收限值。正式试验须核对适用材料、试样、校准、有效性条件和采用版次。