试件几何形状
载荷与裂纹参数
材料参数
关键公式
$K_I = \sigma\,Y\sqrt{\pi a}$
塑性区：$r_p = \dfrac{1}{2\pi}\!\left(\dfrac{K_I}{\sigma_y}\right)^{\!2}$
临界裂纹：$a_c = \dfrac{1}{\pi}\!\left(\dfrac{K_{Ic}}{\sigma Y}\right)^{\!2}$
能量释放率：$G = K_I^2/E'$
什么是应力强度因子
🎓
简单来说，它就是一个用来“量化”裂纹有多危险的数字。想象一下玻璃上的裂缝，有的小裂缝没事，有的大裂缝一碰就碎。这个 $K_I$ 值就是工程师用来判断“这个裂缝在多大压力下会突然炸开”的关键指标。在实际工程中，比如飞机机翼的例行检查，就要算这个值来评估裂纹是否安全。
🧑🎓
诶，真的吗？那这个危险数字是怎么算出来的呢？
🎓
最核心的公式是 $K_I = \sigma Y \sqrt{\pi a}$。这里 $\sigma$ 是材料承受的应力，$a$ 是裂纹长度，$Y$ 是一个“形状修正系数”。你可以在这个模拟器里试试看：拖动“裂纹长度a”的滑块，把它从5毫米慢慢调到20毫米，你会发现计算出来的 $K_I$ 值并不是简单翻倍，而是增加得更快！这就是为什么长裂纹会特别危险的原因。
🧑🎓
那算出来之后，怎么知道安不安全？直接和材料的“断裂韧性” $K_{Ic}$ 比一比就行了吗？
🎓
基本思路是这样，但工程上要考虑更多。当 $K_I$ 接近 $K_{Ic}$ 时，裂纹尖端会进入塑性变形，形成一个“塑性区”。你可以在这个工具里把“屈服强度 $\sigma_y$”参数调低看看，塑性区半径 $r_p$ 会显著变大。如果这个塑性区太大，我们最初的线性弹性假设就不成立了，需要更复杂的分析。所以安全评估是综合看 $K_I$、$r_p$ 和临界裂纹长度 $a_c$ 的。
物理模型与关键公式
这是线性断裂力学的核心控制方程，描述了远场应力、裂纹尺寸和几何形状共同决定的裂纹尖端应力场强度。
$$K_I = \sigma \, Y(a/W) \, \sqrt{\pi a}$$
其中，$K_I$ 为 I 型（张开型）应力强度因子 (MPa√m)；$\sigma$ 为名义应力 (MPa)；$a$ 为裂纹长度 (m)；$Y$ 为几何形状因子，是无量纲数，通常是裂纹长度与试件宽度之比 $a/W$ 的函数。
为了评估断裂安全性和线性断裂力学的适用性，需要计算以下关键衍生参数：
$$r_p = \frac{1}{2\pi}\left( \frac{K_I}{\sigma_y}\right)^2, \quad a_c = \frac{1}{\pi}\left( \frac{K_{Ic}}{\sigma Y}\right)^2, \quad G = \frac{K_I^2}{E‘}$$
$r_p$：平面应力条件下的塑性区半径 (m)，表征裂纹尖端屈服区域大小。$a_c$：临界裂纹长度 (m)，表示在当前应力下导致失稳扩展的裂纹尺寸。$G$：能量释放率 (J/m²)，表示裂纹扩展单位面积所消耗的能量。$E‘$ 为等效弹性模量（平面应力为 $E$，平面应变为 $E/(1-\nu^2)$）。
现实世界中的应用
航空航天结构安全评估：飞机机身、机翼和起落架在循环载荷下会产生疲劳裂纹。工程师使用此类计算器定期评估检测到的裂纹的 $K_I$ 和 $a_c$，以制定合理的检修周期，确保在裂纹达到临界长度前进行修复或更换。
压力容器与管道完整性管理：在石油化工和能源行业，输送易燃易爆介质的压力容器和管道一旦破裂后果严重。通过计算不同缺陷尺寸下的应力强度因子，可以确定容器的最大安全工作压力，或判断已存在的焊接缺陷、腐蚀坑是否需要立即处理。
大型机械与桥梁的疲劳寿命预测：对于承受交变载荷的火车轮轴、风力发电机主轴和桥梁钢结构，初始缺陷会随着使用逐渐扩展。结合应力强度因子和疲劳裂纹扩展速率模型（如 Paris 定律），可以预测结构从初始裂纹扩展到临界尺寸的剩余寿命。
新材料研发与性能测试：在实验室中，通过制备标准裂纹试样（如紧凑拉伸CT试样），测量其断裂韧性 $K_{Ic}$。这个过程本身就依赖于精确的应力强度因子解（即几何因子 $Y$）。计算器帮助研究人员快速分析实验数据，评估陶瓷、复合材料或高强钢等新材料的抗断裂性能。
常见误解与注意事项
开始使用此工具时，有几个需要特别注意的要点。首先，“计算结果始终是基于理想化模型的参考值”。工具中使用的公式基本以无限大平板中的裂纹这一简化模型为基础。实际零件常存在棱角、孔洞或多个裂纹，此时修正系数Y会变得复杂。例如，直径100mm的轴存在表面裂纹的情况与宽大平板中心裂纹的情况相比，即使σ和a相同，K_I也完全不同。建议先通过工具掌握基本概念，在实际设计中务必通过更详细的标准（如JSME或ASTM规范）或有限元分析进行验证。
其次是“塑性区半径”的理解。若计算出的r_p大于裂纹长度a的1/10，则需引起警惕。这表明线弹性断裂力学的基本前提（裂纹尖端塑性区足够小）开始失效。例如，对于σ_y=300MPa的软钢，计算得到K_I=100 MPa√m时，r_p约为1.4mm。若裂纹长度a=5mm则问题不大，但若a=2mm则塑性区尺寸已不可忽略。此时需要采用弹塑性断裂力学的分析方法。
最后要注意“断裂韧性K_Ic会随条件变化”。材料表中给出的K_Ic值通常对应室温静态载荷条件。实际结构可能处于低温、循环载荷或腐蚀环境中。例如，同一钢材在-40℃时K_Ic可能显著低于室温值（脆化现象）。对于动态载荷，则需要另外考虑表征裂纹扩展速率的da/dN-ΔK曲线（疲劳裂纹扩展特性）。即使通过工具得出一个结果，也应始终追问：“实际使用环境如何？”
相关工程领域
应力强度因子K的概念不仅适用于此处学习的脆性断裂，更在众多工程领域发挥着核心作用。首当其冲的是“疲劳分析”。循环载荷下的裂纹扩展速率与应力强度因子范围ΔK (=K_max - K_min) 密切相关。通过这种关系（帕里斯定律），可以预测裂纹从亚临界长度扩展到临界长度a_c所需的循环次数，从而进行机械寿命设计。
另一个重要领域是“应力腐蚀开裂”。当应力与腐蚀环境共存时，裂纹可能在远低于材料固有K_Ic的应力强度因子下扩展。该临界值称为“应力腐蚀开裂门槛应力强度因子K_Iscc”。例如，铝合金和高强度钢在特定环境下可能面临K_Iscc问题。通过本工具计算K_I，并与该材料在对应环境下的K_Iscc进行比较，即可评估应力腐蚀开裂风险。
进一步拓展，该理论还可应用于“界面断裂”和“复合材料断裂”分析。要描述粘接界面或涂层剥离、碳纤维增强塑料中的裂纹扩展等问题，不仅需要考虑I型（张开型）应力强度因子，还需考虑II型（面内剪切）和III型（面外剪切）应力强度因子。本工具的基础正是学习这些更复杂断裂模式的第一步。
进阶学习指引
若想了解本工具计算公式的背景或进行深入学习，建议按以下顺序展开。首先是数学基础：建议追踪裂纹尖端应力场公式的推导过程。工具中使用的 $$K_I = \sigma Y \sqrt{\pi a}$$ 实际上源于基于复变函数理论的“韦斯特加德应力函数”弹性解。理解这推导过程，就能明白为何K_I与√a成正比。推荐参考《断裂力学基础》（如日本机械学会编写版本）的前几章。
接下来可学习“弹塑性断裂力学”。该领域以J积分和CTOD（裂纹尖端张开位移）为核心参数。这些概念对于评估裂纹尖端存在大范围塑性区的“韧性”材料（如管道用钢）的断裂行为至关重要。当本工具计算显示塑性区半径r_p较大时，正是EPFM的应用场景。
最后的实践阶段，建议通过“基于有限元法的断裂分析”进行体验。现代CAE软件（如Abaqus、ANSYS、FrontSTR等）都具备裂纹建模及直接计算J积分或应力强度因子的功能。若能掌握“简单形状用工具计算，复杂实际结构用FEA验证”的分工方法，便已成为合格的断裂力学实践者。