金属積層造形の熱解析 — 理論と支配方程式
概要
先生！ 今日は金属積層造形の熱解析の話なんですよね？ どんなものなんですか？
レーザーPBF/DEDの温度場予測。溶融池、残留応力。
本記事では金属積層造形の熱解析の理論的基礎、支配方程式、離散化手法、および主要商用ツールでの実装について詳しく見ていこう。
へぇ〜！ レーザーについてだいぶ理解が深まりました。メモメモ…📝
支配方程式
いよいよ数式ですね…！ 金属積層造形の熱解析ではどんな方程式が出てくるんですか？
金属積層造形の熱解析の基本となる方程式をこんな感じだよ。
数学的に書くと、こんな形になるんだ。
えっと…各項はどんな物理現象を表してるんですか？
ここで各変数は問題に応じた物理量を表す。上記の支配方程式は、適切な境界条件（Dirichlet条件、Neumann条件、混合条件）と初期条件のもとで一意解を持つ。
ここまで聞いて、金属積層造形の熱解析がなぜ重要か、やっと腹落ちしました！
離散化手法
この方程式を、コンピュータで実際にどうやって解くんですか？
有限要素法（FEM）による空間離散化を使うんだ。要素剛性マトリクスを組み立て、全体剛性方程式を構築する。
へぇ〜！ 有限要素法についてだいぶ理解が深まりました。メモメモ…📝
行列解法アルゴリズム
行列解法アルゴリズムって、具体的にはどういうことですか？
直接法（LU分解、Cholesky分解）または反復法（CG法、GMRES法）により連立方程式を解く。大規模問題では前処理付き反復法が効果的なんだ。
へぇ〜！ 有限要素法についてだいぶ理解が深まりました。メモメモ…📝
| 解法 | 分類 | メモリ使用量 | 適用規模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小〜中規模 |
| Cholesky分解 | 直接法（対称正定値） | O(n²) | 小〜中規模 |
| PCG法 | 反復法 | O(n) | 大規模 |
| GMRES法 | 反復法 | O(n·m) | 大規模・非対称 |
| AMG前処理 | 前処理 | O(n) | 超大規模 |
つまり有限要素法のところで手を抜くと、後で痛い目を見るってことですね。肝に銘じます！
商用ツールにおける実装
で、金属積層造形の熱解析をやるにはどんなソフトが使えるんですか？
| ツール名 | 開発元/現在 | 主要ファイル形式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
ベンダーの系譜と製品統合の経緯
各ソフトの成り立ちって、結構ドラマチックだったりしますか？
Ansys Mechanical (旧ANSYS Structural)
「Ansys Mechanical」について教えてください！
1970年にSwanson Analysis Systems Inc. (SASI) が開発。APDL（Ansys Parametric Design Language）ベース。
現在の所属: Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEAって、具体的にはどういうことですか？
1978年にHKS (Hibbitt, Karlsson & Sorensen) が開発。2005年にDassault Systèmesが買収し、SIMULIAブランドに統合。
現在の所属: Dassault Systèmes SIMULIA
ここまで聞いて、が開発がなぜ重要か、やっと腹落ちしました！
COMSOL Multiphysics
「COMSOL Multiphysics」について教えてください！
1986年スウェーデンで設立。MATLAB連携のFEMLABとして開始、後にCOMSOLに改名。マルチフィジックスに強み。
現在の所属: COMSOL AB
おお〜、が開発の話、めちゃくちゃ面白いです！ もっと聞かせてください。
ファイル形式と相互運用性
異なるソルバー間でモデルを変換する際は、要素タイプの対応関係、材料モデルの互換性、荷重・境界条件の表現差異に注意が必要になるんだ。特に高次要素や特殊要素（コヒーシブ要素、ユーザー定義要素等）はソルバー間で直接変換できない場合が多い。
なるほど…フォーマットって一見シンプルだけど、実はすごく奥が深いんですね。
実務上の注意点
教科書には載ってない「現場の知恵」みたいなものってありますか？
メッシュ収束性の確認、境界条件の妥当性検証、材料パラメータの感度分析がすごく大事なんだ。
- メッシュ依存性の検証: 少なくとも3水準のメッシュ密度で収束性を確認
- 境界条件の妥当性: 物理的に意味のある拘束条件の設定
- 結果の検証: 理論解、実験データ、既知ベンチマーク問題との比較
いやぁ、金属積層造形の熱解析って奥が深いですね… でも先生の説明のおかげでだいぶ整理できました！
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
ムーアの法則と冷却の戦い
CPUの集積度は2年で2倍になる（ムーアの法則）。しかし発熱密度もほぼ同じペースで増加。最新のCPUは数百ワットを数cm²の面積で発熱しており、単位面積あたりの発熱密度はホットプレートを超えています。電子機器の熱設計CAEは、まさに「ムーアの法則との終わりなき競争」なのです。
各項の物理的意味
- 蓄熱項 $\rho c_p \partial T/\partial t$：単位体積あたりの熱エネルギー蓄積率。【日常の例】鉄のフライパンは熱しにくく冷めにくいが、アルミ鍋は熱しやすく冷めやすい——これは密度 $\rho$ と比熱 $c_p$ の積（熱容量）の違い。熱容量が大きい物体は温度変化が緩やかになる。水は比熱が非常に大きい（4,186 J/(kg·K)）ため、海沿いの気温は内陸より安定する。非定常解析ではこの項が温度の時間変化速度を決める。
- 熱伝導項 $\nabla \cdot (k \nabla T)$：フーリエの法則に基づく熱伝導。温度勾配に比例した熱流束。【日常の例】金属スプーンを熱い鍋に入れると持ち手まで熱くなる——金属は熱伝導率 $k$ が高いため、高温側から低温側へ素早く熱が伝わる。木製スプーンが熱くならないのは $k$ が小さいから。断熱材（グラスウール等）は $k$ が極めて小さく、温度勾配があっても熱が伝わりにくい。「温度差のあるところに熱が流れる」という自然の傾向を数式化したもの。
- 対流項 $\rho c_p \mathbf{u} \cdot \nabla T$：流体の運動に伴う熱輸送。【日常の例】扇風機に当たると涼しく感じるのは、風（流体の流れ）が体表面近くの暖かい空気を運び去り、新鮮な冷たい空気を供給するから——これが強制対流。暖房で部屋の天井付近が暖かくなるのは、暖められた空気が浮力で上昇する自然対流。PCのCPUクーラーのファンも強制対流で放熱している。対流は熱伝導よりも桁違いに効率的な熱輸送手段。
- 熱源項 $Q$：内部発熱（ジュール熱、化学反応熱、放射線吸収等）。単位: W/m³。【日常の例】電子レンジは食品内部のマイクロ波吸収（体積発熱）で加熱する。電気毛布のヒーター線はジュール発熱（$Q = I^2 R / V$）で暖かくなる。リチウムイオン電池の充放電時の発熱、ブレーキパッドの摩擦熱も熱源として解析で考慮される。外部から「表面」に熱を与える境界条件とは異なり、熱源項は「内部」でのエネルギー生成を表す。
仮定条件と適用限界
数値例：平板の定常熱伝導（厚み10mm, 鋼k=50W/(m·K), 表面100°C/裏面20°C）
熱流束 q = k×ΔT/L = 50×80/0.01 = 400,000 W/m² 各位置の温度は線形分布
材料別の熱伝導率の比較（数値が大きいほど熱を伝えやすい）：
銅は空気の約15,000倍も熱を伝えやすい！ ヒートシンクに銅やアルミが使われる理由がこのグラフで一目瞭然です。
熱解析の境界条件設定は経験と試行錯誤の繰り返し。 — Project NovaSolverは、実務者の知見を活かしやすい解析環境の実現を研究しています。
金属積層造形の熱解析の実務で感じる課題を教えてください
Project NovaSolverは、CAEエンジニアが日々直面する課題——セットアップの煩雑さ、計算コスト、結果の解釈——の解決を目指しています。あなたの実務経験が、より良いツール開発の原動力になります。
実務課題アンケートに回答する →