试件几何形状
载荷与裂纹参数
材料参数
理论与主要公式
$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曲线（疲劳裂纹扩展特性）。即使通过工具得出一个结果，也应始终追问：“实际使用环境如何？”