AM微細組織シミュレーション
理論：GとRが組織を決める
概要
AMの微細組織シミュレーションでは、結晶の一つひとつまで計算するのですか？
目的によって3段階ある。いちばん簡単なのは、部品スケールの熱解析で凝固界面の温度勾配Gと凝固速度Rを求め、材料ごとの凝固マップに当てはめて「どこが柱状晶で、どこが等軸晶か」「組織の細かさはどの程度か」を予測する方法で、計算は軽く、設計の比較に向く。次がセル・オートマトン法（CA）で、結晶粒の核生成と成長をセルの状態の変化として計算し、数mm角の領域で粒の形や向きを表す。最も細かいのがフェーズフィールド法で、デンドライトの枝や元素の偏析まで解くが、扱えるのは数十〜数百µm角にとどまる。どの方法でも出発点は温度場なので、まず熱解析でGとRを正しく求めることが基本になる。
冷却速度・凝固速度・セル間隔
Ṫ は冷却速度、v は走査速度、θ は凝固界面の法線と走査方向のなす角、λ はセル（一次デンドライト）間隔。G×R が大きいほど組織は細かく、G/R が大きいほど平滑・セル状に、小さいほど等軸晶になりやすい。
溶接の研究から積層造形へ
溶融池が移動しながら固まるという点で、レーザー積層造形はレーザー溶接とよく似ている。冷却速度とセル間隔の関係は、1980年代にステンレス鋼のレーザー溶接の研究で整理されたものが、そのまま積層造形の組織の予測にも使われている。ただし積層造形では、同じ場所が後の層で何度も再加熱されるため、一度できた組織が焼き戻されたり変態したりする点が溶接と大きく違う。
計算例
セル間隔は示した経験式への代入値で、結晶粒径や強度の直接予測ではない。砂型の1 K/sは式の外挿例であり、適用範囲外の測定値ではない。G/Rだけでは柱状晶・等軸晶は決まらず、核生成や合金組成も必要。
例1：プロセスによる冷却速度とセル間隔（316L）
| プロセス | G（目安） | R（目安） | 冷却速度 | セル間隔 |
|---|---|---|---|---|
| LPBF | 10⁶ K/m | 1 m/s | 10⁶ K/s | 0.84 µm |
| DED | 10⁵ K/m | 0.01 m/s | 10³ K/s | 8.2 µm |
| 砂型鋳造 | — | — | 約1 K/s | 80 µm |
GとRは桁の目安。同じ316Lでも、LPBFの組織は鋳造より約100倍細かく、これがLPBF材の高い降伏強さの一因になっている。
例2：溶融池の中でのG/Rの変化（走査速度1 m/s）
| 位置 | 界面の向き θ | 凝固速度 R | G（目安） | G/R |
|---|---|---|---|---|
| 溶融池の底 | 89° | 0.017 m/s | 10⁷ K/m | 5.7×10⁸ K·s/m² |
| 中ほど | 45° | 0.71 m/s | 3×10⁶ K/m | 4.2×10⁶ K·s/m² |
| 上部（後端） | 0° | 1.0 m/s | 10⁶ K/m | 1.0×10⁶ K·s/m² |
凝固速度は界面の向きで決まり、溶融池の底では走査速度のごく一部しかない。G/Rは底から上部へ約570倍小さくなり、底ではセル状・柱状、上部ほど等軸晶が出やすくなる。
上部で等軸晶が出やすいなら、なぜ完成品の断面は造形方向に長い柱状晶ばかりなのですか？
次の層が上部を溶かし直すからだ。LPBFでは溶融池の深さが層の厚さ（数十µm）の2〜3倍あることが多く、上部にできた等軸晶の層は次の層の走査でほとんど溶けてなくなる。残るのは各層の溶融池の下の部分、つまりG/Rが大きく柱状に育った部分で、それが次の層の結晶の土台になって同じ向きのまま上へ上へとつながる（エピタキシャル成長）。だから組織を予測するには、1層だけでなく複数層の再溶融を含めて計算する必要がある。等軸晶を増やしたい場合は、溶融池を浅く広くする条件や、結晶の核になる粒子を加えた合金が使われ、シミュレーションはその条件探しに使われる。
実務での進め方
- 溶融池スケールの熱解析で、溶融池の形を断面写真に合わせて熱源を調整する。
- 凝固界面でのGとRを取り出し、冷却速度とG/Rの分布を求める。
- 凝固マップや経験式で組織の形と細かさを予測し、組織写真と比べる。
- 結晶粒の形や向きが必要なら、複数層の再溶融を含めてセル・オートマトン法で解く。
- 組織から強さの異方性を見積もり、部品の向きや造形条件の選択に使う。
「溶融池の形がすべての出発点」
組織予測がうまくいかないとき、原因の多くは組織のモデルではなく、その前の熱解析にある。レーザーの吸収率は粉末か固体か、溶けているかで大きく変わり、溶融池の深さが2倍違えば底でのGとRも、再溶融される量も変わる。単一トラックを造形して切断し、溶融池の幅と深さを測って熱源を合わせる作業を最初に行うと、その後の組織予測がずっと安定する。
よくある間違い
間違いと対策
| 間違い | 起きること | 対策 |
|---|---|---|
| 熱源を調整しない | GとRが誤る | 単一トラックの断面に合わせる |
| 1層だけ計算 | 柱状晶を見逃す | 複数層の再溶融を含める |
| Rを走査速度とする | 底の凝固速度を過大評価 | 界面の向きで補正 |
| 吸収率を一定 | 溶融池の深さがずれる | 状態ごとの吸収率 |
| 再加熱による変態を無視 | 硬さ・相の予測がずれる | 温度履歴で相変態を評価 |
関連する内容も知りたいです。
関連するページはLPBF(レーザー粉末床融合)シミュレーション、DED(指向性エネルギー堆積)シミュレーション、AM残留応力解析、鋳造凝固解析が参考になる。
よくある質問
AMの微細組織は何で決まりますか？
凝固界面での温度勾配Gと凝固速度（界面の進む速さ）Rの2つでほぼ決まります。積G×Rは冷却速度で組織の細かさを、比G/Rは組織の形（平滑界面・セル・柱状デンドライト・等軸晶）を決めます。熱解析でGとRを求め、凝固マップに当てはめるのが基本の考え方です。
LPBFの組織はどれくらい細かいですか？
LPBF（レーザー粉末床溶融結合）では冷却速度が約10⁶ K/sに達し、316Lステンレスのセル間隔は経験式で約0.84 µmです。冷却速度10³ K/s程度のDED（指向性エネルギー堆積）では約8 µm、砂型鋳造（約1 K/s）では約80 µmで、同じ材料でもプロセスで100倍違います。
なぜAMでは柱状晶が多くなるのですか？
溶融池の底では凝固速度Rが小さくG/Rが大きいため、下の層の結晶を引き継いで上向きに伸びる柱状晶になりやすいからです。走査速度1 m/sの例では、G/Rは溶融池の底で約5.7×10⁸、上部で約1×10⁶ K·s/m²と約570倍違います。造形方向に伸びた柱状晶は強さの異方性の原因になります。
参考資料・検算の位置づけ
掲載値は本文の仮定に基づく教育用の参考計算。資料は理論・設定の参照先であり、本記事の数値を各ソフトで再現した証明ではない。設計・製造の判断には対象条件での検証が必要。
関連トピック
なった
詳しく
報告