パラメータ設定
既定値は軟鋼相当（σ_0=100 MPa、k_H=600 MPa·√μm）に d=25 μm の中粒組織、σ_target=500 MPa の高強度鋼目標。粒径 d は平均直径として扱い、ホール・ペッチが成立する目安は d ≥ 100 nm（ナノ粒では逆ホール・ペッチに注意）。
多結晶構造の模式図
不規則な多角形セルで多結晶を表現。粒径 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 は √10≈3.16 倍になるから、強化項 k_H/√d は 120 → 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)²=1.5²=2.25 μm。これは TMCP（熱間制御圧延・加速冷却）で得られる典型的な粒径レンジで、SM490 や HT780 級の高強度鋼を粒界強化だけで作る場合の目標値だ。実務では析出強化や転位密度との合算で達成するけど、本ツールは「粒径だけでどこまで行けるか」を切り分けて見せてくれる。σ_target を σ_0 以下にすれば d_req は表示されない（粒界強化なしで達成可能）から試してみて。
🙋
グラフが直線なのが面白いです。横軸を 1/√d にすると σ_y が直線になるって、何か便利なんですか？
🎓
めちゃくちゃ便利だよ。実験で粒径と降伏応力を測ったとき、横軸 1/√d、縦軸 σ_y にプロットして直線になれば「ホール・ペッチが成立している」、y 切片から σ_0、傾きから k_H が求められる。もし直線が湾曲したり下に折れたら、粒径 100 nm 以下で起こる「逆ホール・ペッチ」（粒界滑りで軟化）に入った可能性が高い。研究者はこのプロットを Hall-Petch plot と呼んで、合金の強化機構を切り分ける標準ツールにしている。本ツールでも k_H を 50→1500 にスライドすると直線の傾きが急峻になり、強い粒界障壁を持つ材料（HCP の Mg、bcc の 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$ は平均結晶粒径で、線分法（intercept method）または 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 から d=5 μm に変えると、σ_y が 220 MPa → 368 MPa に上昇することが確認できます。これは LNG タンカー、海洋構造物、X70/X80 ラインパイプの高強度化を支える基本原理で、合金元素を増やさずに強度を上げるためコストと溶接性の両方で有利です。
ECAP/HPT による超微細粒材料：等チャンネル角押出（ECAP）や高圧捻り（HPT）といった激しい塑性加工（SPD）は、純銅・純チタン・アルミ合金の結晶粒を 100〜500 nm まで微細化します。本ツールの粒径下限 0.1 μm（=100 nm）に近い領域で σ_y を確認すると、粒界強化だけで降伏応力が 5〜10 倍に跳ね上がることがわかります。航空宇宙用 Ti-6Al-4V では従来 20 μm 粒で σ_y=900 MPa 程度ですが、SPD で 200 nm 化すると 1500 MPa 級が達成可能で、軽量化と強度を両立できます。
HSLA 鋼の合金設計：High-Strength Low-Alloy 鋼（API 5L X70 など）は、Nb、V、Ti などの微量元素を析出物（NbC、VC、TiN）として粒成長を抑制し、最終粒径 5〜10 μm を実現します。本ツールで k_H を変えずに d を 25→5 μm に縮めると σ_y が約 1.7 倍になり、これがパイプライン用鋼の経済的な高強度化の核心であることがわかります。Nb 量 0.05% の添加で粒径制御コストは大きく下がり、大量生産にも適合します。
結晶粒径測定と非破壊評価：粒径は EBSD（電子線後方散乱回折）や光学顕微鏡による線分法で測定しますが、超音波の減衰係数からも推定できます。発電所配管の経年劣化評価では、減衰係数 → 粒径 → ホール・ペッチで残存強度を逆算する手法が使われ、本ツールのような σ_y vs 1/√d 直線が現場の余寿命診断ツールとして実用化されています。
よくある誤解と注意点
最も多い誤解は、「粒径を小さくすればするほど無限に強くなる」というものです。粒径が 10〜30 nm を下回ると、転位機構が粒界滑り（grain boundary sliding）や Coble クリープに置き換わり、σ_y はむしろ低下する「逆ホール・ペッチ」現象が観察されます。本ツールの粒径下限 0.1 μm は通常の塑性加工で得られる最小レンジに対応していますが、ナノ結晶（10 nm 以下）の評価には別途分子動力学や複合モデルが必要です。100 nm 以上のサブミクロン粒径ならホール・ペッチがよく成り立ちます。
次に多いのが、「k_H は材料固有の定数」という誤解です。実際は同じ材料でも温度、変形速度、合金組成、固溶元素濃度で k_H が変動します。例えば bcc の鉄は低温で k_H ≈ 700 MPa·√μm ですが、室温の純銅では 110 程度。窒素や炭素の固溶量が 0.01% 増えるだけで k_H が 1.5 倍になる例もあり、文献値を使うときは測定温度・組成・前処理を必ず確認してください。本ツールでは k_H を広範囲（50〜1500）にスライドできるため、各種材料の比較が可能です。
最後に、「σ_y の予測精度がホール・ペッチだけで十分」という思い込みです。実際の鋼や合金では、粒界強化に加えて固溶強化（σ_ss = K·c^n）、析出強化（オロワン機構）、転位強化（σ_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（熱間制御圧延・加速冷却）で典型的に得られる粒径レンジで、High-Strength Low-Alloy 鋼の合金設計指針となります。σ_target が σ_0 以下なら粒界強化なしで達成できるため d_req は表示されません。