針金を何度も折り曲げていると、だんだん固くなって最後にポキッと折れる — この経験は誰にでも覚えがあるでしょう。最初は軽い力で曲がるのに、途中から急に抵抗感が増す。これが 加工硬化（ひずみ硬化、work hardening / strain hardening）です。
「変形するたびに硬くなる」という現象は一見不思議に思えます。変形しても原子は無くならないし、何かが加わるわけでもない。では、なぜ硬くなるのでしょうか？答えは金属の内部に隠れています。結晶格子の中を動き回る「転位（dislocation）」という線状の欠陥が、変形とともに急増し、互いに絡まり合って動けなくなる — これが加工硬化の本質です。
加工硬化の理解は、実用上も非常に重要です。主な応用先を2つ挙げましょう。
プレス加工・深絞り成形: 自動車の車体パネルや飲料缶のアルミは、プレスによる塑性変形で最終形状に仕上げられます。このとき「加工硬化指数 $n$」が大きい材料ほど、成形中に局所的な薄肉化（くびれ）が起きにくく、複雑な形状に成形しやすいのです。
コールドワーキング（冷間加工）による強化: 鋼線やピアノ線は、室温での伸線加工によって引張強さを飛躍的に高めます。熱処理を使わずに転位を増殖させるだけで、純鉄の約200 MPaを1,000 MPa超に強化できます。これはトン単位の材料を扱う産業では、エネルギーと時間の大幅な節約につながります。
本記事では、加工硬化の現象を応力-ひずみ曲線で確認するところから始め、転位論に基づく物理的メカニズムを掘り下げ、Taylor式とHollomon式という2つの重要な定量関係を導出します。さらに、加工硬化を「逆転」させる焼鈍の仕組みと、他の強化機構（固溶強化・析出強化・結晶粒微細化）との比較まで解説します。
前提知識
この記事を読む前に、以下の記事を読んでおくと理解が深まります。
応力-ひずみ曲線で見る加工硬化
まず、金属の「応力-ひずみ曲線」を改めて見てみましょう。この曲線には、材料が変形する過程で経験するすべての段階が凝縮されています。
引張試験を行うと、変形の初期は 弾性域 です。応力 $\sigma$ はひずみ $\varepsilon$ に比例し（フックの法則 $\sigma = E \varepsilon$）、荷重を取り除けば元に戻ります。この段階では転位はほとんど動きません。
ひずみが増すと 降伏点 に達します。軟鋼では上降伏点と下降伏点の間に「降伏棚」と呼ばれる水平部分が現れ、応力がほぼ一定のままひずみが進みます。これは刃状転位の先端に固定されていた炭素原子が外れて、転位が一斉に動き始める現象（Cottrellの雰囲気）です。
降伏棚を過ぎると、曲線は再び上昇します。これが 加工硬化域（ひずみ硬化域） です。変形が進むほど応力は増し続けます。この「さらに変形させるのに、より大きな力が必要になる」という現象が加工硬化そのものです。
そして曲線はやがてピーク（引張強さ、UTS）に達し、その後ネッキング（くびれ）が始まって応力が低下し、最終的に破断します。
グラフを見ると、降伏後のひずみ範囲（赤いシェーディング部）で応力が継続的に上昇していることがわかります。降伏点（約400 MPa）から引張強さ（約780 MPa）まで、ひずみ 0.004 から 0.30 の広い範囲で硬化が続いています。このエネルギー蓄積域が、冷間加工による強化のエンジンです。
それでは、なぜ変形すると硬くなるのか？その答えを求めて、結晶の内部へと降りていきましょう。
加工硬化の物理的メカニズム：転位の増殖と絡み合い
転位とは何か
金属の塑性変形は、「すべり面」上での「転位」の移動によって起こります。転位とは結晶格子の線状の乱れで、その端から端へと移動することで、巨視的には結晶の一部がすべり方向に1つのバーガースベクトル $\bm{b}$ 分だけずれる効果をもたらします。
未変形の金属にも転位は存在します。その密度は $\rho \approx 10^{10}$ m$^{-2}$ 程度です。$\rho$ は「単位体積あたりの転位線の総長さ」で、単位はm/m³ = m$^{-2}$です。
転位の増殖：Frank-Readソース
変形が進むとき、転位はどこから来るのでしょうか？一見すると「あるだけ動かせば尽きてしまう」ように思えますが、実際には変形中に転位が 爆発的に増殖 します。
主要な増殖機構が Frank-Readソース です。両端がピン止めされた転位セグメント（長さ $L$）に応力をかけると、中央部がふくらんで円弧を描き、最終的にループを形成して新たな転位を生み出します。1サイクルごとに転位が1本ずつ生成され、繰り返すことで多数の転位ループが連続生成されます。生成に必要な臨界せん断応力は
$$ \tau_c = \frac{G b}{L} $$
です。ここで $G$ はせん断剛性、$b$ はバーガースベクトルの大きさです。この式から、ピン間距離 $L$ が短いほど（＝すでに密に転位が存在するほど）、新たな転位を動かすのにより大きな応力が必要になることがわかります。
転位の絡み合い（転位の森）
増殖した転位はすべて同じ方向に動くわけではありません。金属の結晶には複数のすべり系（すべり面とすべり方向の組み合わせ）があり、異なるすべり系の転位が互いに 障害物 となります。
2本の転位が交差すると、交差点で ジョグ（jog） やキンク（kink）という段差が生まれ、転位の移動が妨げられます。特にジョグは転位がすべり面外に飛び出した点なので、そこから先へ転位が動けなくなります（ジョグが固着点となる）。
このようにして変形が進むにつれて転位は増殖し、互いに絡まり合い、固着されて動けない転位が増えていきます。この状態を 転位の森（forest of dislocations） と呼びます。動ける転位が動くためには、この「森」を突き抜けなければなりません。森が濃くなるほど、突き抜けるのに必要な応力も増す — これが加工硬化の本質です。
3つのパネルは変形の進行に伴う転位の状態変化を模式的に示しています。左（初期）では転位の数が少なく互いの距離も遠い状態です。中央（変形中）では転位が増殖して密度が増し、右（加工硬化後）では無数の転位が複雑に絡まり合った「転位の森」が形成されています。この段階では、どの転位も自由には動けません。
転位密度の増加
変形度と転位密度の関係を定量的に見てみましょう。未変形の金属では $\rho \sim 10^{10}$ m$^{-2}$ 程度ですが、冷間加工を繰り返すと
$$ \rho \sim 10^{12}\text{（軽加工）} \sim 10^{15}\text{（重加工）}\quad\text{m}^{-2} $$
へと指数的に増加します。対数スケールで4〜5桁の増加です。
グラフでは、冷間加工度が増えるにつれて転位密度が急激に増加する様子が示されています。特に加工度0〜30%の初期段階で増加が著しく、その後は上昇が緩やかになります（転位同士が消滅・合体することで純増が鈍化する）。橙色のライン（$10^{12}$）と赤色のライン（$10^{15}$）は、軽加工と重加工の典型的な転位密度の参照値です。
ここまでで、変形によって転位密度が増加することがわかりました。次は、この転位密度の増加が「なぜ応力の増加につながるのか」を定量的に説明する Taylor 式を導きます。
Taylor式：転位密度と流動応力の関係
物理的背景
1つの動く転位が転位の森を突き抜けるには、「森の転位」を切断しなければなりません。転位同士の弾性的な相互作用（反発・引力）がこの障壁を生みます。
転位の弾性エネルギーは、その周囲の応力場から来ます。転位1本の単位長さあたりの弾性エネルギーは
$$ E_{\text{line}} \approx \frac{G b^2}{2} $$
です。この「線張力」が転位の移動に対する抵抗の源です。
平均障害物間距離
転位密度が $\rho$ のとき、動く転位が次の「森の転位」に出会うまでの平均距離（平均自由行程） $l$ は
$$ l \approx \frac{1}{\sqrt{\rho}} $$
と見積もれます。これは転位線が平面上にランダムに分布していると仮定したときの、最近接距離の期待値に相当します。
臨界分解せん断応力の導出
動く転位が障害転位を乗り越えるには、力のバランスからせん断応力が必要です。線張力 $T \approx G b^2 / 2$ を持つ転位が曲率半径 $R = l/2 = 1/(2\sqrt{\rho})$ で曲がるときに釣り合う応力は
$$ \tau \cdot b = \frac{T}{R} = \frac{G b^2 / 2}{1/(2\sqrt{\rho})} = G b^2 \sqrt{\rho} $$
両辺を $b$ で割り、幾何学的因子 $\alpha$（通常 0.2〜0.5）を加えると
$$ \boxed{\tau = \tau_0 + \alpha G b \sqrt{\rho}} $$
これが Taylor式 です。$\tau_0$ は格子摩擦応力（ペイエルスナバロ応力）で転位のない完全結晶でも存在する固有の抵抗です。この式は「流動応力は転位密度の平方根に比例して増加する」という、実験でよく確認される関係を明快に予測します。
Taylor式の各項の意味を整理しておきましょう。$\tau_0$ は材料固有の「ベースの強さ」、$\alpha G b$ は強化の感度（材料定数）、$\sqrt{\rho}$ が加工硬化の駆動変数です。 $G$ と $b$ は材料で決まりますが、 $\rho$ は加工によって変化します。
Taylor式に数値を代入してみましょう。銅（$G = 48$ GPa, $b = 2.56 \times 10^{-10}$ m, $\alpha = 0.3$, $\tau_0 = 5$ MPa）の場合：
- 未変形（$\rho = 10^{10}$ m$^{-2}$）: $\tau = 5 + 0.3 \times 48000 \times 2.56 \times 10^{-10} \times 10^5 \approx 5 + 3.7 = 8.7$ MPa
- 重加工（$\rho = 10^{15}$ m$^{-2}$）: $\tau = 5 + 0.3 \times 48000 \times 2.56 \times 10^{-10} \times \sqrt{10^{15}} \approx 5 + 117 = 122$ MPa
転位密度が5桁増えると、強度は約14倍になる計算です。実際の実験値とよく対応しています。
グラフでは横軸を $\sqrt{\rho}$（対数スケール）にとって Taylor式のプロファイルを示しています。転位密度 $10^{10}$（未変形）での臨界せん断応力は小さく、密度が増すにつれて $\sqrt{\rho}$ に比例して線形に増加する様子が読み取れます。このグラフにより、転位密度を制御することが強化の直接的な手段であることが視覚的に理解できます。
Taylor式が定性的・定量的な骨格を与えてくれましたが、エンジニアが実際に設計計算で使うのは真応力-真ひずみの関係式です。次はその Hollomon 式を見ていきましょう。
Hollomon式：加工硬化指数 $n$
真応力と公称応力の違い
まず、「真応力」と「公称応力」の区別が重要です。
- 公称応力 $\sigma_n = F / A_0$: 初期断面積 $A_0$ に対する荷重
- 真応力 $\sigma = F / A$: 変形後の実際の断面積 $A$ に対する荷重
塑性変形では体積が保存されます（$V = A_0 L_0 = A L$）。したがって変形が進んで断面積が減ると、公称応力は過小評価になります。真応力と公称ひずみ $e$ の間には
$$ \sigma = \sigma_n (1 + e) $$
の関係があります（小ひずみでは差が小さいですが、大変形では重要）。
同様に、真ひずみ $\varepsilon = \ln(1 + e)$ は公称ひずみとは異なります。
Hollomon式の形と意味
真応力-真ひずみ曲線の塑性域（加工硬化域）は、多くの金属材料で次の Hollomon式（べき乗則）でよく近似されます:
$$ \boxed{\sigma = K \varepsilon^n} $$
ここで: – $\sigma$: 真応力 [MPa] – $\varepsilon$: 真ひずみ（塑性成分） – $K$: 強度係数 [MPa]（$\varepsilon = 1$ での強度に相当） – $n$: 加工硬化指数（0 〜 1 の無次元数）
$n = 0$ は硬化なし（完全塑性体）、$n = 1$ は線形硬化（弾性に相当）です。実用金属では $n = 0.1 \sim 0.5$ の範囲に収まります。
両対数プロットで $n$ を決める
Hollomon式の両辺の対数をとると
$$ \ln \sigma = \ln K + n \ln \varepsilon $$
となり、$\ln \sigma$ vs $\ln \varepsilon$ の直線になります。傾きが $n$、切片が $\ln K$ です。引張試験データを両対数プロットして直線フィットすることで $n$ と $K$ を実験的に求められます。
左パネルは通常スケールで、Hollomon式のべき乗的な上昇がわかります。右パネルは両対数スケールで、データ点が直線に乗っており傾きが $n = 0.22$ であることが読み取れます。両対数プロットが直線になることは Hollomon式の成立を示す重要な実験的確認です。また、切片（$\varepsilon = 1$ での外挿値）から $K$ が求まります。
$n$ の物理的意味：ネッキング安定性
$n$ には深い物理的意味があります。引張変形中にある断面で細くなり始める（くびれが生じる）と、その箇所では応力が増します。もし材料が十分に加工硬化するなら、細くなった箇所がより硬化して変形が他の部分に分散します。この「自己安定化」がどれだけ効くかを表すのが $n$ です。
Considèreの条件から、ネッキングが始まる真ひずみは $\varepsilon = n$ です。したがって:
- $n$ が大きい材料 → ネッキング開始が遅い → 大きな均一伸びが得られる → 深絞り・プレス成形に有利
- $n$ が小さい材料 → ネッキングが早期に始まる → 成形性は低いが、加工の初期段階での精密加工に向く
このグラフでは5種類の材料を比較しています。純銅（$n = 0.54$）は初期から急速に硬化しながら広いひずみ範囲で変形が続くのに対し、高強度鋼（$n = 0.08$）はほとんど硬化せず早期に飽和することが読み取れます。ステンレス鋼（SUS304, $n = 0.45$）は純銅に次いで高い硬化能を示し、加工による強化が大きいことがわかります。アルミ合金（Al6061, $n = 0.20$）と低炭素鋼（$n = 0.22$）はほぼ同等の硬化挙動を示します。
冷間加工と硬さの変化
冷間加工とは
冷間加工（cold working） とは、材料の再結晶温度以下（金属の種類にもよりますが、鉄では約450°C以下、アルミでは約150°C以下）での塑性加工を指します。室温加工が典型例です。
代表的な冷間加工には以下があります： – 伸線（wire drawing）: 金属線を小さなダイスを通して引き伸ばし、細くする – 圧延（cold rolling）: 2本のロールで金属板を薄く伸ばす – 鍛造（cold forging）: 型に圧入して形状を変える – 引抜き（extrusion）: ダイスから押し出す
これらの操作により、材料内の転位密度が増大し、Hollomon式に従って強度（降伏応力・引張強さ）が上昇します。硬さ（ビッカース硬さHVやロックウェル硬さ）も同様に上昇します。
加工度と硬さの実験的関係
加工度（reduction in area, RA）は次で定義されます:
$$ \text{加工度} [\%] = \frac{A_0 – A}{A_0} \times 100 $$
加工度と硬さの関係を純銅と純アルミで見てみましょう。
グラフから読み取れる重要な点は2つです。第一に、どちらの材料も加工初期（0〜30%）に急速に硬化します。これは Frank-Read ソースによる転位の急増期に対応しています。第二に、加工度が70%を超えると硬化が飽和し、グラフが水平に近づきます。これは転位の増殖と消滅（転位間の反応や回復）がバランスする状態に近づくためです。純銅は純アルミより絶対的な硬さが高く、またより高い加工度まで硬化が続きます。
冷間加工の弊点：加工限界
冷間加工を続けると加工硬化が進む一方で、延性と靭性が低下 します。転位が過密になると局所的な応力集中が生じ、微小き裂が発生しやすくなります。針金を繰り返し折り曲げると最終的に折れる（加工誘起破壊）のはこのためです。
また、過度の冷間加工は内部残留応力を高め、部品の寸法精度や後加工性にも悪影響を与えます。
このジレンマ（強化したいが延性も欲しい）を解決するのが「焼鈍」です。次のセクションで見ていきましょう。
焼鈍による回復・再結晶・粒成長
焼鈍とは
焼鈍（アニーリング, annealing） とは、冷間加工した金属を適切な温度に加熱・保持して、加工によって蓄積したひずみエネルギーを解放し、軟化させる熱処理です。加工硬化の「逆転操作」と言えます。
焼鈍には温度に応じて3つの段階があります：回復（recovery）→ 再結晶（recrystallization）→ 粒成長（grain growth）。
回復（recovery）
比較的低温（例: 銅で150〜350°C）での加熱段階です。転位は動ける程度の熱エネルギーを得て、互いに消滅（正負の転位が出会って消える）したり、規則正しい「転位セル構造（subgrain）」を形成したりします。
回復では： – 硬さの低下は緩やか（20〜30%程度） – 結晶粒の形状は変わらない（引き伸ばされたまま） – 内部残留応力が大幅に低下する
残留応力を取り除きながら強度を維持したい場合に有効な処理です。
再結晶（recrystallization）
温度が上がる（銅で350〜550°C程度）と、加工されて引き伸ばされた結晶粒の中に 新たな無欠陥の結晶粒（核）が発生し、それが成長 して元の組織を置き換えます。これが再結晶です。
再結晶が完了すると： – 転位密度が未加工レベル（$10^{10}$ m$^{-2}$）に激減 – 硬さが急激に低下（もとの軟らかい状態に近づく） – 延性・靭性が回復
再結晶温度 $T_r$ は経験的に融点 $T_m$（絶対温度）の約0.3〜0.5倍です：$T_r \approx (0.3 \sim 0.5) T_m$
粒成長（grain growth）
再結晶後もさらに加熱を続けると、小さな粒が大きな粒に飲み込まれ、粒が粗大化します（粒成長）。粒成長は表面エネルギーを減らす方向に進みます。
粒成長すると： – 粒界面積が減り、粒界による強化効果（Hall-Petch効果）が弱まる – 硬さ・強度がさらに低下 – 粒が粗大すぎると表面粗さが悪化（橙皮状組織）
上パネルは硬さ、下パネルは平均粒径の温度変化を示しています。橙色の「回復」域では硬さが緩やかに低下し粒径はほぼ不変です。赤色の「再結晶」域で硬さが急落し、緑色の「粒成長」域では粒径が急増することが読み取れます。設計では、再結晶温度のやや上で止めると（軟化させつつ過度な粒成長を抑える）、バランスの取れた材料特性が得られます。
金属の強化機構：4つのアプローチ
加工硬化は強化機構の一つにすぎません。金属の強化には大きく分けて以下の4つの戦略があり、実用合金ではこれらが重複して使われます。
1. 固溶強化（solid solution strengthening）
基材金属の格子に異種原子を溶かし込む（固溶させる）方法です。置換型固溶（サイズが近い原子）と侵入型固溶（格子間に入る小さな原子）があります。
異種原子は格子を歪め、その歪み場が転位の移動を妨げます。有名な例は炭素を格子間に含む鉄（鋼）で、わずか0.8%Cで降伏応力が大幅に上昇します。
$$ \Delta \tau_{\text{solid}} \propto c^{1/2} \text{ または } c^{2/3} $$
（$c$: 溶質濃度）
2. 析出強化（precipitation strengthening / age hardening）
合金を高温から急冷して過飽和固溶体を作り、その後低温で保持すると微細な析出物（例: アルミニウム合金中のMgSi、Al-Cu合金中のθ’相）が析出します。
転位がこの析出物（障害物）を切断する（剪断機構）か回りを抜ける（Orowan機構）かによって強化されます。Orowan機構の増加応力は
$$ \Delta \tau_{\text{Orowan}} = \frac{G b}{L_p} $$
（$L_p$: 析出物間の平均距離）
析出強化は最も効果的な強化手段の一つで、7000系アルミ合金（航空機材）や時効硬化ステンレス鋼などで活用されています。
3. 結晶粒微細化（grain refinement） ― Hall-Petch則
粒界は転位の通過に対して障壁になります。転位が粒界に達すると、隣の粒でスリップを開始するために応力集中が必要です。粒が細かいほど粒界が多く、転位の平均自由行程が短くなります。
この関係を定式化したのが Hall-Petch則:
$$ \boxed{\sigma_y = \sigma_0 + \frac{k}{\sqrt{d}}} $$
ここで: – $\sigma_y$: 降伏応力 [MPa] – $\sigma_0$: 格子摩擦応力（単結晶の降伏応力）[MPa] – $k$: Hall-Petch係数（材料定数）[MPa·m$^{1/2}$] – $d$: 平均結晶粒径 [m]
$d$ を4分の1にすれば（粒径を半分にすれば2倍に）$\sigma_y$ が $k/\sqrt{d}$ の項で2倍増加します。ナノ結晶材料では粒径を数十nmまで小さくすることで著しい強化が得られます（ただし極端に微細化するとHall-Petch則が崩れる「逆Hall-Petch」領域があることも知られています）。
4. 加工硬化（work hardening）
本記事の主題です。Taylor式 $\tau = \tau_0 + \alpha G b \sqrt{\rho}$ に従い、転位密度 $\rho$ を増やすことで強化します。
これら4つの機構は独立ではなく、重畳して作用します。Taylor式の $\tau_0$ の中に固溶強化や析出強化の寄与を組み込んで議論することもあります。
棒グラフは各強化機構の降伏応力への寄与を積み上げて示しています。純金属のベース（50 MPa）から出発し、固溶強化（+60 MPa）、析出強化（+120 MPa）、結晶粒微細化（+80 MPa）、加工硬化（+150 MPa）が積み上がり、合計460 MPaになっています。実際の高強度アルミ合金（7075-T6など）はこのような複合強化で500 MPa超の降伏応力を達成しています。加工硬化の寄与が最大であることに注目してください。
Hall-Petch則の詳細：粒径と強度
Hall-Petch則は「粒径を小さくすれば強くなる」という、シンプルかつ強力な関係式です。材料設計では非常に重要な指針となります。
導出の背景：ピールアップモデル
Hall-Petch則の物理的背景を簡単に説明します。あるすべり系で動いてきた転位が粒界に積み重なる（pile-up）と考えます。$n$ 本の転位が積み重なったとき、最先端の転位が感じる応力は印加応力の $n$ 倍になります。粒径 $d$ と転位のpile-up長さは比例するため
$$ n \propto \frac{\tau \cdot d}{G b} $$
この応力集中が隣の粒で降伏を引き起こすための臨界値を $\tau_c$ とすれば
$$ \tau_c = \frac{\text{const}}{d^{1/2}} $$
これがHall-Petch関係の定性的な導出です。正確な係数は実験的に決まります。
Pythonで Hall-Petch則を可視化・フィットする
実際の実験データから $\sigma_0$ と $k$ を決める手順をコードで示します。
import numpy as np
import matplotlib.pyplot as plt
from scipy.optimize import curve_fit
# 模擬実験データ: 低炭素鋼の結晶粒径 vs 降伏応力
# d [μm], σy [MPa]
d_data = np.array([200, 100, 50, 25, 10, 5, 2, 1]) # μm
sigma_data = np.array([168, 183, 200, 220, 257, 295, 358, 428]) # MPa
# Hall-Petch 式フィット
def hall_petch(d, sigma0, k):
return sigma0 + k / np.sqrt(d)
popt, _ = curve_fit(hall_petch, d_data, sigma_data, p0=[150, 100])
sigma0_fit, k_fit = popt
print(f"フィット結果: sigma0 = {sigma0_fit:.1f} MPa, k = {k_fit:.1f} MPa*sqrt(μm)")
# 結果: sigma0 = 160.8 MPa, k = 275.9 MPa*sqrt(μm)
# Hall-Petch プロット (1/√d vs σy)
inv_sqrt_d = 1 / np.sqrt(d_data)
inv_sqrt_fit = np.linspace(0, 1.1, 200)
sigma_fit_line = sigma0_fit + k_fit * inv_sqrt_fit
# R^2 計算
residuals = sigma_data - hall_petch(d_data, *popt)
ss_res = np.sum(residuals**2)
ss_tot = np.sum((sigma_data - sigma_data.mean())**2)
r2 = 1 - ss_res / ss_tot
print(f"R^2 = {r2:.4f}")
# R^2 = 0.9919
フィット結果は $\sigma_0 = 160.8$ MPa、$k = 275.9$ MPa·$\mu$m$^{1/2}$ であり、$R^2 = 0.992$ と非常によく直線に乗っています。このことは、Hall-Petch則が広い粒径範囲で実験をよく記述することを示しています。$1/\sqrt{d}$ vs $\sigma_y$ の線形プロットが直線になれば、データが Hall-Petch 則に従うことの実験的証明となります。
左パネルでは粒径が細かくなるほど降伏応力が急上昇する様子（特に $d < 10$ μm での強化が顕著）がわかります。右パネルでは $1/\sqrt{d}$ vs $\sigma_y$ が完全な直線になり、Hall-Petch 則の妥当性が視覚的に確認できます。切片が $\sigma_0 = 50$ MPa（格子摩擦応力）、傾きが $k = 20$ MPa·$\mu$m$^{1/2}$（Hall-Petch係数）に対応しています。
数値計算でまとめる：加工硬化の総合シミュレーション
ここまで学んだ Taylor式・Hollomon式・Hall-Petch則を Python で実際に計算し、各パラメータの影響を確認しましょう。
import numpy as np
import matplotlib
import matplotlib.pyplot as plt
for cand in ["Hiragino Sans", "Yu Gothic", "Noto Sans CJK JP", "IPAexGothic", "Meiryo"]:
if any(cand == f.name for f in matplotlib.font_manager.fontManager.ttflist):
plt.rcParams["font.family"] = cand
break
plt.rcParams["axes.unicode_minus"] = False
# ==============================
# 1. Taylor式: 転位密度 → 強度
# ==============================
tau0 = 10 # MPa (ペイエルス応力、格子摩擦)
alpha = 0.3 # 幾何学的因子 (典型値)
G = 80000 # MPa (せん断剛性、軟鋼的)
b = 2.5e-10 # m (バーガースベクトル)
rho_range = np.logspace(10, 15, 200) # m^-2
tau_taylor = tau0 + alpha * G * b * np.sqrt(rho_range)
print("=== Taylor式 ===")
for rho_val in [1e10, 1e12, 1e14]:
tv = tau0 + alpha * G * b * np.sqrt(rho_val)
print(f" rho = {rho_val:.0e}: tau = {tv:.1f} MPa")
# rho = 1e10: tau = 10.6 MPa
# rho = 1e12: tau = 16.0 MPa
# rho = 1e14: tau = 70.0 MPa
# ==============================
# 2. Hollomon式: ひずみ → 応力
# ==============================
K_steel, n_steel = 800, 0.22 # 低炭素鋼
K_copper, n_copper = 500, 0.54 # 純銅
eps_range = np.linspace(0.005, 0.45, 300)
sigma_steel = K_steel * eps_range**n_steel
sigma_copper = K_copper * eps_range**n_copper
# Considere条件: ネッキング開始ひずみ = n
eps_neck_steel = n_steel
eps_neck_copper = n_copper
print("\n=== Hollomon式 (Considere条件) ===")
print(f" 低炭素鋼: ネッキング開始ひずみ = {eps_neck_steel:.2f}")
print(f" 純銅: ネッキング開始ひずみ = {eps_neck_copper:.2f}")
# 低炭素鋼: ネッキング開始ひずみ = 0.22
# 純銅: ネッキング開始ひずみ = 0.54
# ==============================
# 3. Hall-Petch則: 粒径 → 降伏応力
# ==============================
sigma0_hp = 50 # MPa
k_hp = 20 # MPa*μm^0.5
d_range = np.array([1, 5, 10, 25, 50, 100, 200]) # μm
sigma_y_hp = sigma0_hp + k_hp / np.sqrt(d_range)
print("\n=== Hall-Petch則 ===")
for d_val, sy in zip(d_range, sigma_y_hp):
print(f" d = {d_val:3d} μm: σy = {sy:.1f} MPa")
# d = 1 μm: σy = 70.0 MPa
# d = 10 μm: σy = 56.3 MPa
# d = 100 μm: σy = 52.0 MPa
# d = 200 μm: σy = 51.4 MPa
# ==============================
# 4. 複合強化の計算
# ==============================
sigma_solid = 60 # MPa (固溶強化)
sigma_precip = 120 # MPa (析出強化)
sigma_grain = 80 # MPa (結晶粒微細化)
sigma_work = 150 # MPa (加工硬化)
sigma_base = 50 # MPa (純金属)
sigma_total = sigma_base + sigma_solid + sigma_precip + sigma_grain + sigma_work
print(f"\n=== 複合強化の合計降伏応力 ===")
print(f" 基本強度: {sigma_base} MPa")
print(f" 固溶強化: +{sigma_solid} MPa")
print(f" 析出強化: +{sigma_precip} MPa")
print(f" 粒微細化: +{sigma_grain} MPa")
print(f" 加工硬化: +{sigma_work} MPa")
print(f" 合計: {sigma_total} MPa")
# 合計: 460 MPa
計算結果を確認します。Taylor式では転位密度が $10^{10}$ から $10^{14}$ に増えると、せん断応力が10.6 MPaから70.0 MPaへと約6.6倍上昇しています。Hollomon式からは、純銅は低炭素鋼の2.5倍近い真ひずみまでネッキングなく変形できることがわかります（均一伸びが大きい）。複合強化の計算では、各機構が独立に加算されると仮定すると460 MPaという高い降伏応力が得られ、高強度合金の現実的な強度レベルと一致しています。
実用材料の加工硬化指数と強度係数
主要な工業材料の Hollomon パラメータと特徴をまとめます。
| 材料 | $K$ [MPa] | $n$ | 特徴 |
|---|---|---|---|
| 純銅（アニール） | 530 | 0.54 | 高 $n$、深絞り・成形に優秀 |
| 低炭素鋼 | 800 | 0.22 | バランス型、汎用プレス材 |
| ステンレス鋼 SUS304 | 1300 | 0.45 | 高 $n$、高強度×高成形性 |
| Al合金 AA6061-O | 400 | 0.20 | 軽量、成形後T6処理で析出強化 |
| チタン合金 Ti-6Al-4V | 1400 | 0.10 | 低 $n$、高強度・難加工 |
| 高強度低合金鋼 HSLA | 1500 | 0.08 | 極低 $n$、加工硬化より固溶・析出強化主体 |
SUS304のような低合金ステンレス鋼が高い $n$ を持つのは、変形誘起マルテンサイト変態（TRIP効果）が関与するためです。変形中にオーステナイトがマルテンサイトに変態し、追加の転位増殖を引き起こします。これを利用したのがTRIP鋼（Transformation Induced Plasticity steel）で、自動車構造材に広く使われています。
まとめ
本記事では、加工硬化（ひずみ硬化）を転位論から定量的に解説しました。
- 応力-ひずみ曲線 の降伏後領域（加工硬化域）では、変形とともに流動応力が上昇する。これは転位の増殖と絡み合いによって転位密度 $\rho$ が増大するためである
- Taylor式 $\tau = \tau_0 + \alpha G b \sqrt{\rho}$ は、流動応力が $\sqrt{\rho}$ に比例して増加することを示す。転位の森を突き抜けるコストが硬化の物理的本質
- Hollomon式 $\sigma = K \varepsilon^n$ は工学的な定量ツール。両対数プロットで直線になり、傾きの $n$（加工硬化指数）はネッキング安定性の指標でもある（Considère条件: ネッキング開始ひずみ $\varepsilon_{\text{neck}} = n$）
- 冷間加工 により硬さと強度は上昇するが、延性は低下する。加工限界を超えると加工誘起破壊が生じる
- 焼鈍 により回復（残留応力低減）→ 再結晶（軟化）→ 粒成長（粗大化）という3段階で加工硬化を解除できる
- 強化機構の整理: 固溶強化・析出強化・結晶粒微細化（Hall-Petch則）・加工硬化は独立かつ重畳して働く。実用高強度合金はこれらを組み合わせて設計されている
加工硬化の次のステップとして、疲労（繰り返し加工硬化と軟化のサイクル）や破壊力学（加工硬化域でのき裂進展）に進むと、材料の変形と破損を統一的に理解できます。