AM残留応力解析
理論と物理
概要
先生！ 今日はAM残留応力解析の話なんですよね？ どんなものなんですか？
積層造形(AM)プロセスにおける急速加熱・冷却サイクルに起因する残留応力の予測。固有ひずみ法(Inherent Strain Method)やlumped layer法による効率的な解析手法が実用的。
先生の説明分かりやすい！ 積層造形のモヤモヤが晴れました。
支配方程式
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
固有ひずみ法：
待って待って、固有ひずみ法ってことは、つまりこういうケースでも使えますか？
理論的基盤
「理論的基盤」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
AM残留応力解析のシミュレーションは、熱力学・材料力学・流体力学の連成問題として定式化される。製造プロセスの物理現象は複数の時間・空間スケールにまたがるため、マクロスケールの連続体モデルとメゾ/ミクロスケールの材料モデルの適切な組合せが求められる。プロセスパラメータ（温度、速度、荷重等）と製品品質（寸法精度、欠陥、機械特性）の因果関係を定量的に予測することが目標なんだ。
製造プロセスの支配方程式
数式は苦手なんですけど…AM残留応力解析の式の「意味」を教えてもらえますか？
製造プロセスシミュレーションは、熱力学・流体力学・固体力学の連成問題として定式化される。
熱伝導方程式（エネルギー保存）
熱伝導方程式って、具体的にはどういうことですか？
ここで $T$ は温度、$\mathbf{v}$ は材料の速度場、$k$ は熱伝導率、$Q$ は内部発熱（ジュール熱、潜熱、摩擦熱等）なんだ。
先輩が「製造プロセスシミュレだけはちゃんとやれ」って言ってた意味が分かりました。
凝固・相変化
「凝固・相変化」について教えてください！
凝固過程では潜熱の放出/吸収が温度場に大きく影響する。エンタルピー法による定式化：
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
ここで $L$ は潜熱、$f_l(T)$ は液相率（固液共存域で0から1の間の値をとる）。
塑性変形の構成則
塑性変形の構成則って、具体的にはどういうことですか？
金属の塑性変形はJohnson-Cook構成則等で記述される：
$A$: 初期降伏応力、$B$: 硬化係数、$n$: 硬化指数、$C$: 歪み速度感度、$m$: 温度軟化指数。
ここまで聞いて、製造プロセスシミュレがなぜ重要か、やっと腹落ちしました！
流動解析（充填・鋳造）
次は流動解析の話ですね。どんな内容ですか？
溶融金属や樹脂の流動はナビエ-ストークス方程式に従うが、高粘性・非ニュートン流体特性を考慮する必要がある。射出成形ではCross-WLFモデルが標準的：
なるほど…製造プロセスシミュレって一見シンプルだけど、実はすごく奥が深いんですね。
仮定と適用限界
この式って万能じゃないんですか？ 使えない場面ってどんなとき？
無次元パラメータと支配的スケール
先生、「無次元パラメータと支配的スケール」について教えてください！
解析対象の物理現象を支配する無次元パラメータの理解は、適切なモデル選択とパラメータ設定の基盤となる。
あっ、そういうことか！ 解析対象の物理現象をってそういう仕組みだったんですね。
次元解析による検証
「次元解析による検証」について教えてください！
解析結果のオーダー推定には、バッキンガムのΠ定理に基づく次元解析が効果的なんだ。代表長さ $L$、代表速度 $U$、代表時間 $T = L/U$ を用いて、各物理量のオーダーを事前に推定し、解析結果の妥当性を確認する。
なるほど。じゃあ解析対象の物理現象をができていれば、まずは大丈夫ってことですか？
境界条件の分類と数学的特徴
境界条件って、ここを間違えると全部ダメになるって聞いたんですけど…
| 種類 | 数学的表現 | 物理的意味 | 例 |
|---|---|---|---|
| ディリクレ条件 | $u = u_0$ on $\Gamma_D$ | 変数値の指定 | 固定壁、温度指定 |
| ノイマン条件 | $\partial u/\partial n = g$ on $\Gamma_N$ | 勾配（フラックス）の指定 | 熱流束、力 |
| ロビン条件 | $\alpha u + \beta \partial u/\partial n = h$ | 変数と勾配の線形結合 | 対流熱伝達 |
| 周期境界条件 | $u(x) = u(x+L)$ | 空間的周期性 | 単位セル解析 |
適切な境界条件の選択は解の一意性と物理的妥当性に直結するんだよ。不足した境界条件は不適切な問題となり、過剰な境界条件は矛盾を生じさせる。
AM残留応力解析の全体像がつかめました！ 明日から実務で意識してみます。
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
AM残留応力の「引張と圧縮のサンドイッチ構造」
AM残留応力の理論を学ぶと、まず驚くのが「表面が引張、内部が圧縮」という直感に反した分布パターンです。これはLPBFの冷却過程を想像するとわかりやすくて、上層のビードが急冷で収縮しようとすると、すでに固化した下層がそれを拘束する。結果として上層には引張応力、下層（内部）には圧縮応力がたまる——ちょうどカステラを乾燥させたときに表面だけひび割れるイメージです。この応力分布がわかると、なぜ「表面研磨で疲労強度が改善するか」が理論的に説明できます。引張応力の強い表面を削り落とすことで、圧縮が表面に露出して亀裂進展を抑制するんですね。
各項の物理的意味
- 保存量の時間変化項：対象とする物理量の時間的変化率を表す。定常問題では零となる。【イメージ】浴槽にお湯を張るとき、水位が時間と共に上がる——この「時間あたりの変化速度」が時間変化項。バルブを閉じて水位が一定になった状態が「定常」であり、時間変化項はゼロ。
- フラックス項（流束項）：物理量の空間的な輸送・拡散を記述する。対流と拡散の2種類に大別される。【イメージ】対流は「川の流れがボートを運ぶ」ように流れに乗って物が運ばれること。拡散は「インクが静止した水中で自然に広がる」ように濃度差で物が移動すること。この2つの輸送メカニズムの競合が多くの物理現象を支配する。
- ソース項（生成・消滅項）：物理量の局所的な生成または消滅を表す外力・反応項。【イメージ】部屋の中でヒーターをつけると、その場所に熱エネルギーが「生成」される。化学反応で燃料が消費されると質量が「消滅」する。外部から系に注入される物理量を表す項。
仮定条件と適用限界
- 連続体仮定が成立する空間スケールであること
- 材料・流体の構成則（応力-歪み関係、ニュートン流体則等）が適用範囲内であること
- 境界条件が物理的に妥当かつ数学的に適切に定義されていること
次元解析と単位系
| 変数 | SI単位 | 注意点・換算メモ |
|---|---|---|
| 代表長さ $L$ | m | CADモデルの単位系と一致させること |
| 代表時間 $t$ | s | 過渡解析の時間刻みはCFL条件・物理的時定数を考慮 |
数値解法と実装
数値手法の詳細
具体的にはどんなアルゴリズムでAM残留応力解析を解くんですか？
AM残留応力解析のシミュレーションに用いる数値手法を解説する。
離散化手法
大変形を伴う製造プロセスでは、Updated Lagrangian法またはALE(Arbitrary Lagrangian-Eulerian)法が一般的に用いられる。接触問題にはペナルティ法またはラグランジュ乗数法を使うんだ。Euler法を用いた定常流れ場定式化は鍛造・押出しなどの定常プロセスに有効なんだ。
時間積分
先生、「時間積分」について教えてください！
準静的問題では陰解法(Newton-Raphson)、高速変形・衝撃問題では陽解法(中心差分法)を選択する。質量スケーリングにより陽解法の時間ステップ制限を緩和できるが、運動エネルギーが内部エネルギーの5-10%以下であることを監視する必要がある。
メッシュ管理
メッシュって細かければ細かいほどいいんですよね？ …あれ、違いますか？
大変形によるメッシュ歪みに対してリメッシング(r-adaptivity)やALEメッシュスムージングを使うんだ。SPH法やMPM(Material Point Method)等のメッシュフリー手法も選択肢となる。
待って待って、大変形によるメッシュってことは、つまりこういうケースでも使えますか？
接触・摩擦のモデリング
「接触・摩擦のモデリング」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
製造プロセスでは工具と被加工材の接触が不可避であり、接触アルゴリズムの選択が解の精度と安定性を左右する。クーロン摩擦、せん断摩擦、温度依存摩擦モデルを工程に応じて使い分ける。接触検出のペナルティパラメータやセグメント対セグメント法の設定が計算安定性に大きく影響する。
先輩が「製造プロセスでは工具だけはちゃんとやれ」って言ってた意味が分かりました。
数値解法の実装詳細
先生、「数値解法の実装詳細」について教えてください！
メッシュ要件
メッシュ要件って、具体的にはどういうことですか？
製造プロセスシミュレーションでは、移動する界面（固液界面、自由表面）を追跡する必要があるため、メッシュ戦略がすごく大事なんだ。
「接触・摩擦のモデリング」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
製造プロセスでは工具と被加工材の接触が不可避であり、接触アルゴリズムの選択が解の精度と安定性を左右する。クーロン摩擦、せん断摩擦、温度依存摩擦モデルを工程に応じて使い分ける。接触検出のペナルティパラメータやセグメント対セグメント法の設定が計算安定性に大きく影響する。
先輩が「製造プロセスでは工具だけはちゃんとやれ」って言ってた意味が分かりました。
先生、「数値解法の実装詳細」について教えてください！
メッシュ要件って、具体的にはどういうことですか？
製造プロセスシミュレーションでは、移動する界面（固液界面、自由表面）を追跡する必要があるため、メッシュ戦略がすごく大事なんだ。
| 手法 | 概要 | 適用 |
|---|---|---|
| ALE法 | メッシュが材料と共に移動 | 鍛造、圧延 |
| オイラー法 | 固定メッシュ上で材料が流動 | 鋳造充填 |
| VOF法 | 体積分率で自由表面を追跡 | 鋳造、射出成形 |
| CEL法 | 結合オイラー-ラグランジュ | 衝撃加工 |
| SPH法 | 粒子法、メッシュフリー | AM溶融池 |
熱源モデル（溶接・AM）
熱源モデルって、具体的にはどういうことですか？
Goldak二重楕円体モデル：
ここで $P$ はレーザー/アーク出力、$\eta$ は吸収効率、$a,b,c$ は楕円体の半軸長。
つまりメッシュ要件のところで手を抜くと、後で痛い目を見るってことですね。肝に銘じます！
時間積分
時間積分って、具体的にはどういうことですか？
つまりメッシュ要件のところで手を抜くと、後で痛い目を見るってことですね。肝に銘じます！
連成ソルバー戦略
次は連成ソルバー戦略の話ですね。どんな内容ですか？
あっ、そういうことか！ メッシュ要件ってそういう仕組みだったんですね。
誤差評価と精度検証
「誤差評価と精度検証」って聞いたことはあるんですけど、ちゃんと理解できてないかもしれません…
離散化誤差の評価
離散化誤差の評価って、具体的にはどういうことですか？
リチャードソン外挿法による離散化誤差の推定：
ここで $f_h$ はメッシュ幅 $h$ での解、$r$ はメッシュ比、$p$ は離散化の次数。
GCI（Grid Convergence Index）
「GCI」について教えてください！
ASME V&V 20-2009に基づくメッシュ収束性の定量評価：
ここまで聞いて、離散化誤差の評価がなぜ重要か、やっと腹落ちしました！
これを数式で表すとこうなるよ。
うーん、式だけだとピンとこないです… 何を表してるんですか？
安全係数 $F_s = 1.25$（3水準以上のメッシュ比較時）。GCI < 5% を収束の目安とする。
先輩が「離散化誤差の評価だけはちゃんとやれ」って言ってた意味が分かりました。
検証ベンチマーク問題
「検証ベンチマーク問題」について教えてください！
解析結果の信頼性を担保するため、以下のベンチマーク問題との比較を推奨：
| 分野 | ベンチマーク | 参照解 |
|---|---|---|
| 構造 | パッチテスト | 一様応力場の再現 |
| 構造 | Scordelis-Loの屋根 | 参照変位 |
| 流体 | 蓋駆動キャビティ | Ghia et al. (1982) |
| 熱 | 1D解析解 | $T(x) = T_0 + (T_1-T_0)x/L$ |
高速化手法
先生、「高速化手法」について教えてください！
AM残留応力解析の全体像がつかめました！ 明日から実務で意識してみます。
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。
「固有歪み法」——なぜ航空宇宙メーカーはこれを使うのか
AM残留応力の数値解法として、フルFEM熱機械連成解析は精度が高い反面、大型部品への適用では計算コストが膨大になります。そこで現場でよく使われるのが「固有歪み法（Inherent Strain Method）」。メルトプール近傍のミクロスケール解析で得た固有歪みをビード単位・層単位のマクロモデルに入力として与えることで、計算時間をフル解析の1/100以下に短縮できます。Boeing、Airbus、GEアビエーションといった大手が採用しているのも、この理由です。精度はフル連成に劣りますが「反り予測の精度90%以上、計算時間2〜3時間以内」という実用的なバランスを取っています。数値解法の選択は、正確さと速さのトレードオフそのものです。
実践ガイド
実践ガイド
先生、「実践ガイド」について教えてください！
AM残留応力解析の実務的な解析手順とベストプラクティスを解説する。
待って待って、残留応力解析の実務的ってことは、つまりこういうケースでも使えますか？
解析フロー
最初の一歩から教えてください！ 何から始めればいいですか？
1. 工程条件の定義: プロセスパラメータ（温度、速度、荷重、時間）の整理と範囲設定
2. 材料データの準備: 温度・ひずみ速度依存の構成則パラメータを試験データから同定
4. プロセスシミュレーション実行: 段階的な複雑度の増加と収束性の確認
5. 結果の検証: 実験データとの比較（寸法精度、荷重履歴、温度分布、欠陥位置）
あっ、そういうことか！ 工程条件の定義ってそういう仕組みだったんですね。
ベストプラクティス
先生、「ベストプラクティス」について教えてください！
品質管理と文書化
教科書には載ってない「現場の知恵」みたいなものってありますか？
解析条件（材料データ出典、境界条件の根拠、メッシュ設定の妥当性）を体系的に文書化する。解析結果のレビュープロセスを確立し、実験との比較による精度評価を定量的に記録する。NAFEMS等のベンチマーク問題を用いたソルバー検証も定期的に実施する。
実務的な解析手順
実務でAM残留応力解析を使うときに、いちばん気をつけるべきことは何ですか？
鋳造シミュレーションのワークフロー
鋳造シミュレーションのワークフって、具体的にはどういうことですか？
1. CADモデル準備: 製品形状 + 湯口系 + 押湯 + 冷し金の3Dモデル
2. メッシュ生成: ヘキサドミナント要素推奨。薄肉部は最低3層以上
4. 境界条件: 鋳型-金属間の熱伝達係数（IHTC）。型温の初期設定
5. 充填解析: 注湯速度・温度を設定。空気巻き込みの監視
6. 凝固解析: 充填完了後の温度場解析。引け巣予測
7. 応力解析: 凝固後の残留応力、型抜き後の変形
射出成形シミュレーションのパラメータ設定
次は射出成形シミュレーションのパラの話ですね。どんな内容ですか？
| パラメータ | 典型値 | 影響 |
|---|---|---|
| 樹脂温度 | 200-300°C | 流動性、表面品質 |
| 金型温度 | 40-100°C | 冷却時間、結晶化度 |
| 射出速度 | 50-200 mm/s | ゲート圧、せん断応力 |
| 保圧 | 50-100 MPa | 収縮補償、寸法精度 |
| 冷却時間 | 10-60 s | 生産性、そり変形 |
先生の説明分かりやすい！ 鋳造シミュレーションのモヤモヤが晴れました。
AM（積層造形）シミュレーションの注意点
次は積層造形の話ですね。どんな内容ですか？
品質保証チェックリスト
品質保証チェックリストって、具体的にはどういうことですか？
あっ、そういうことか！ 鋳造シミュレーションってそういう仕組みだったんですね。
プロジェクト管理とワークフロー自動化
全体の流れをざっくり把握したいんですけど、ステップごとに教えてもらえますか？
ディレクトリ構成の推奨
次はディレクトリ構成の推奨の話ですね。どんな内容ですか？
```
project/
├── cad/ # CADモデル
├── mesh/ # メッシュファイル
├── setup/ # 解析設定ファイル
├── results/ # 計算結果
│ ├── case01/
│ ├── case02/
│ └── ...
├── postprocess/ # 後処理スクリプト・画像
├── report/ # レポート
└── validation/ # 検証データ
```
自動化スクリプトの活用
次は自動化スクリプトの活用の話ですね。どんな内容ですか？
パラメトリックスタディやメッシュ収束性確認は、Pythonスクリプトで自動化することで再現性と効率を大幅に向上できるんだよ。
なるほど。じゃあディレクトリ構成の推ができていれば、まずは大丈夫ってことですか？
レビューチェックリスト
「レビューチェックリスト」について教えてください！
1. 入力データ: 材料定数の単位系、CADの寸法精度、メッシュ品質指標
2. 境界条件: 物理的妥当性、過拘束/拘束不足のチェック
3. ソルバー設定: 収束判定基準、時間刻み、出力頻度
4. 結果検証: 力の釣り合い、エネルギーバランス、理論解との比較
5. 感度分析: メッシュ依存性、境界条件の影響、材料パラメータの不確かさ
つまりディレクトリ構成の推のところで手を抜くと、後で痛い目を見るってことですね。肝に銘じます！
報告書作成のポイント
先生、「報告書作成のポイント」について教えてください！
AM残留応力解析の全体像がつかめました！ 明日から実務で意識してみます。
うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。