参数设置
默认值对应于软钢（σ_0=100 MPa、k_H=600 MPa·√μm），中等晶粒组织 d=25 μm，高强度钢目标 σ_target=500 MPa。晶粒尺寸 d 作为平均直径处理，霍尔-佩奇方程式成立的范围约为 d ≥ 100 nm（纳米晶粒时需注意逆霍尔-佩奇软化）。
多晶结构示意图
用不规则多边形单元表示多晶体。晶粒尺寸 d 较小时显示许多小晶粒，d 较大时显示少数大晶粒。晶界线（红色）是转位的障碍。晶粒数根据 d 自动调整。
霍尔-佩奇直线（σ_y vs 1/√d）
横轴为 1/√d (×10⁻¹ /√μm)，纵轴为 σ_y (MPa)。直线的 y 截距为 σ_0，斜率为 k_H。黄点表示当前工作点 (1/√d, σ_y)，红虚线表示 σ_target，橙色点对应 d_req。
理论与主要公式
多晶金属的屈服应力可表示为摩擦应力和晶界强化项的和。
霍尔-佩奇方程式：
$$\sigma_{y} = \sigma_{0} + \frac{k_{H}}{\sqrt{d}}$$
达成目标屈服应力所需的晶粒尺寸：
$$d_{\mathrm{req}} = \left(\frac{k_{H}}{\sigma_{\mathrm{target}} - \sigma_{0}}\right)^{2}$$
晶粒细化至 1/n 时的强化倍率：
$$\frac{\sigma_{y}(d/n)}{\sigma_{y}(d)} = \frac{\sigma_{0} + \sqrt{n}\,k_{H}/\sqrt{d}}{\sigma_{0} + k_{H}/\sqrt{d}}$$
其中 $\sigma_{0}$ 是摩擦应力（Peierls-Nabarro 应力）[MPa]，$k_{H}$ 是霍尔-佩奇系数 [MPa·√μm]，$d$ 是平均晶粒尺寸 [μm]。$\sigma_{0}$ 代表单晶的基础强度，$k_{H}/\sqrt{d}$ 代表晶界对转位堆积的阻碍作用。该方程在晶粒尺寸 d ≥ 100 nm 时成立较好，更小尺寸时出现逆霍尔-佩奇软化现象。
霍尔-佩奇方程式模拟器介绍
🙋
使用默认值时 σ_y=220 MPa，相当于软钢的强度。为什么仅通过输入晶粒尺寸就能决定屈服应力呢？
🎓
问得好。霍尔-佩奇方程式 σ_y = σ_0 + k_H/√d 是从试验观察得出的经验规律，说明多晶金属的屈服应力可分解为「单晶基础强度 σ_0」和「晶界强化贡献 k_H/√d」两部分。在晶界处，晶体取向发生变化，转位无法穿过，必须堆积或被阻挡。晶粒越小，每个晶粒内的转位堆积长度越短，为了推动转位就需要施加更高的应力。默认值 σ_0=100、k_H=600、d=25 μm 代入，得到 σ_y=100+600/√25=100+120=220 MPa，正好对应软钢的屈服点。
🙋
我看到 d/10 时 σ_y 从 220 升到 479 MPa。仅仅是晶粒细化 10 倍就能强化这么多？
🎓
这就是晶粒细化强化的威力所在。将 d 缩小到 1/10，1/√d 就增加到约 3.16 倍，强化项 k_H/√d 从 120 MPa 提升到 379 MPa，接近 3 倍。基础强度 σ_0=100 保持不变，所以总应力达到 479 MPa，强化倍率为 2.18×。实际上，通过 ECAP（等通道转角挤压）或 HPT（高压扭转）使纯铜晶粒从 25 μm 细化到 250 nm，屈服应力可从 70 MPa 跳跃到 400 MPa 以上，这正是晶界强化的典型应用。
🙋
计算出 d_req=2.25 μm，这是说为了达成 σ_target=500 MPa 必须把粒径细化到这个值吗？
🎓
完全正确。d_req=(k_H/(σ_target−σ_0))²=(600/400)²=2.25 μm。这个粒径范围恰好是 TMCP（热轧控制-加速冷却）工艺可以达到的，是 SM490 和 HT780 级高强度钢的典型设计目标。实际生产中会结合析出强化和位错强化来共同达成，但本工具可以将「粒界强化单独的贡献」清晰地分离出来。试试把 σ_target 改到 σ_0 以下，看 d_req 会怎样。
🙋
图中曲线是直线很有意思。为什么把横轴改成 1/√d 就能得到直线？
🎓
这是个非常实用的技巧。在实验中测量粒径和屈服应力后，如果用 1/√d 和 σ_y 作坐标轴画图，得到直线说明「霍尔-佩奇规律成立」；y 截距给出 σ_0，斜率给出 k_H。如果直线弯曲或向下折，往往表示进入了纳米晶范围（<100 nm），引发逆霍尔-佩奇软化。研究人员把这种图叫做 Hall-Petch 图表，是分离合金强化机制的标准工具。在本工具中改变 k_H（从 50 到 1500），你会看到直线斜率随之改变，可以对比不同晶体结构材料（如 Mg、Mo）的特性。
物理模型与主要公式
多晶金属的屈服应力遵循霍尔（E. O. Hall, 1951）和佩奇（N. J. Petch, 1953）独立发现的经验规律。
$$\sigma_{y} = \sigma_{0} + \frac{k_{H}}{\sqrt{d}}$$
$\sigma_{0}$ 是摩擦应力（Peierls-Nabarro 应力），代表单晶的基础强度；$k_{H}$ 是霍尔-佩奇系数，反映晶界对转位的障碍强度。典型值包括：纯铁 σ_0 ≈ 50 MPa、k_H ≈ 700 MPa·√μm；低碳钢 σ_0 ≈ 100 MPa、k_H ≈ 600 MPa·√μm；纯铜 σ_0 ≈ 25 MPa、k_H ≈ 110 MPa·√μm。系数取决于晶体结构（bcc 最大，fcc 中等，hcp 变化范围大）和转位堆积的几何特征。
$d$ 是平均晶粒尺寸，通常用截距法或 Heyn 法测得。粒径分布宽时可采用对数平均或面积加权平均。本工具假设均匀粒径，可直接输入 d 值。
目标屈服应力 σ_target 所需的粒径由反向求解 $d_{\mathrm{req}} = (k_{H}/(\sigma_{\mathrm{target}}-\sigma_{0}))^2$ 得出，是合金设计和加工工艺规划的基本指导式。当 σ_target ≤ σ_0 时，无需晶界强化即可达成，因此 d_req 无定义。
实际应用
TMCP 钢（造船、桥梁、管线）：热轧控制与加速冷却相结合的 TMCP 工艺可将钢的晶粒从 30～50 μm 细化到 5～10 μm，屈服应力从 250 MPa 提升到 450～500 MPa。在本工具中输入 σ_0=100 MPa、k_H=600 MPa·√μm，将 d 从 25 μm 改为 5 μm，可看到 σ_y 从 220 MPa 升至 368 MPa。这一机制支撑了 LNG 运输船、海洋结构物和 X70/X80 级管线的高强度化，无需增加合金元素就能提高强度，降低成本和焊接困难。
ECAP/HPT 超细晶材料：等通道转角挤压（ECAP）和高压扭转（HPT）等严重塑性变形（SPD）手段可将纯铜、纯钛、铝合金的晶粒细化到 100～500 nm 范围。在本工具粒径下限 0.1 μm（即 100 nm）附近，仅晶界强化的贡献就能使屈服应力增加 5～10 倍。航空航天用 Ti-6Al-4V 合金在常规 20 μm 粒径时 σ_y ≈ 900 MPa，经 SPD 细化至 200 nm 后可达 1500 MPa 级，实现轻量化和高强度的双赢。
HSLA 钢合金设计：高强度低合金（HSLA）钢如 API 5L X70 通过添加 Nb、V、Ti 等微量元素，利用析出物（NbC、VC、TiN）抑制晶粒长大，最终获得 5～10 μm 的细晶组织。本工具中保持 k_H 不变，仅将 d 从 25 μm 缩减到 5 μm，可看到 σ_y 上升约 1.7 倍，这正是经济高效强化管线钢的核心原理。仅添加 0.05% 的 Nb 就能显著降低粒径控制成本，适合大规模生产。
晶粒尺寸测量与无损评估：晶粒大小通常用 EBSD（电子背散射衍射）或光学显微镜线截法测得，也可从超声衰减系数反推。在发电站管道老化评估中，利用衰减系数 → 粒径 → 霍尔-佩奇方程推断残余强度，这样的 σ_y vs 1/√d 直线已成为现场寿命诊断工具，被广泛工业化应用。
常见误解与注意事项
最常见的误解是「晶粒越细越强，无限细化可无限增强」。当晶粒尺寸降至 10～30 nm 以下时，会观察到称为逆霍尔-佩奇（inverse Hall-Petch）的软化现象，σ_y 反而下降。这是因为晶界滑移（grain boundary sliding）或 Coble 蠕变接管了主要变形机制，转位强化机理失效。实用范围是 100 nm 以上的亚微米晶粒，本工具的粒径范围（0.1～500 μm）正好覆盖这个有效区域。处理更细的纳米晶时需借助分子动力学模拟或特殊蠕变模型。
第二个常见误解是「k_H 是材料的绝对常数」。实际上 k_H 随温度、应变速率、合金成分、固溶元素浓度而变化。同一钢材在室温下 k_H ≈ 700 MPa·√μm，在高温时可能大幅下降；氮气或碳固溶增加 0.01%，k_H 可增加 1.5 倍。使用文献数据时必须核实测试条件。本工具允许 k_H 在 50～1500 范围内调节，可对比不同材料的相对强化能力。
第三个误解是「仅用霍尔-佩奇就能精确预测实际屈服应力」。实际合金的强度来自多个机制叠加：固溶强化（σ_ss = K·c^n）、析出强化（Orowan 机制）、位错强化（σ_d = αGb√ρ）、相变强化（马氏体）等。实用合金采用重叠模型将各项相加，霍尔-佩奇单独成立的情况仅限于纯金属或晶界外机制微小的简单合金。本工具是用来「分离晶粒效应」的分析工具，需要理解其局限性。
常见问题
σ_y = σ_0 + k_H/√d 表明多晶金属的屈服应力 σ_y 可以写成摩擦应力 σ_0 与晶界强化项 k_H/√d 的和。σ_0 是对应于单晶 Peierls-Nabarro 应力的基础分量，k_H 是霍尔-佩奇系数，表示晶界阻碍转位运动的强度。当晶粒尺寸 d 减小时，1/√d 增大，屈服应力单调增加。本工具使用默认值（σ_0=100 MPa、k_H=600 MPa√μm、d=25 μm）时，得到 σ_y=220 MPa，这对应于纯铁或软钢的典型屈服应力范围。
金属的塑性变形由转位运动引起，但在晶界处晶体取向改变，转位无法直接穿过，必须激活相邻晶粒的滑移系或在晶界处堆积（堆积）。晶粒越小，每个晶粒内转位堆积的长度越短，为了释放转位所需的应力就越高。这就是 k_H/√d 项的物理起源。本工具将 d 从 25 μm 减小到 2.5 μm 时，σ_y 从 220 增加到 479 MPa，约 2.18 倍，可以量化地验证通过加工热处理或 ECAP（等通道转角挤压）实现的超细晶粒化强度倍增机制。
当晶粒尺寸降低到约 10～30 nm 以下时，观察到称为逆霍尔-佩奇（inverse Hall-Petch）的软化现象，σ_y 反而下降。这是因为晶界滑移（grain boundary sliding）或三重点处的 Coble 蠕变成为主导，转位机制失效。实际上，对于 100 nm 以上的亚微米晶粒尺寸，霍尔-佩奇方程式成立得很好。本工具的晶粒尺寸范围（0.1～500 μm）对应其有效区域。处理 100 nm 以下的情况时，需要结合原子尺度模拟或特殊蠕变模型。
d_req = (k_H/(σ_target − σ_0))² 是达成 σ_y = σ_target 所需的平均晶粒尺寸，用于设计新合金或加工工艺。例如，对于 σ_0=100 MPa、k_H=600 MPa√μm 的钢，若要达成 σ_target=500 MPa，本工具计算得 d_req=2.25 μm。这是 TMCP（热轧控制-加速冷却）工艺通常达到的晶粒尺寸范围，也是高强度低合金钢的合金设计指标。当 σ_target ≤ σ_0 时无需晶界强化即可达成，因此不显示 d_req。