鋳造凝固解析 — 数値解法と実装

鋳造凝固解析のシミュレーションに用いる数値手法を解説する。

大変形を伴う製造プロセスでは、Updated Lagrangian法またはALE(Arbitrary Lagrangian-Eulerian)法が一般的に用いられる。接触問題にはペナルティ法またはラグランジュ乗数法を使うんだ。Euler法を用いた定常流れ場定式化は鍛造・押出しなどの定常プロセスに有効なんだ。

準静的問題では陰解法(Newton-Raphson)、高速変形・衝撃問題では陽解法(中心差分法)を選択する。質量スケーリングにより陽解法の時間ステップ制限を緩和できるが、運動エネルギーが内部エネルギーの5-10%以下であることを監視する必要がある。

大変形によるメッシュ歪みに対してリメッシング(r-adaptivity)やALEメッシュスムージングを使うんだ。SPH法やMPM(Material Point Method)等のメッシュフリー手法も選択肢となる。

製造プロセスでは工具と被加工材の接触が不可避であり、接触アルゴリズムの選択が解の精度と安定性を左右する。クーロン摩擦、せん断摩擦、温度依存摩擦モデルを工程に応じて使い分ける。接触検出のペナルティパラメータやセグメント対セグメント法の設定が計算安定性に大きく影響する。

製造プロセスシミュレーションでは、移動する界面（固液界面、自由表面）を追跡する必要があるため、メッシュ戦略がすごく大事なんだ。

Goldak二重楕円体モデル：

式にするとこう。一つずつ見ていこう。

ここで $P$ はレーザー/アーク出力、$\eta$ は吸収効率、$a,b,c$ は楕円体の半軸長。

• 陽解法: CFL条件による時間刻み制限。衝撃問題に向いているよ。
• 陰解法: 無条件安定。大きな時間刻みが可能だが各ステップで連立方程式を解く。
• 半陰的手法: 対流項を陽的、拡散項を陰的に扱う。

熱-力学連成：各時間ステップで温度場→応力場を逐次的に解く（弱連成）か、同時に解く（強連成）。射出成形では流動-冷却-構造の3場連成が必要になるんだ。

リチャードソン外挿法による離散化誤差の推定：

数学的に書くと、こんな形になるんだ。

ここで $f_h$ はメッシュ幅 $h$ での解、$r$ はメッシュ比、$p$ は離散化の次数。

ASME V&V 20-2009に基づくメッシュ収束性の定量評価：

これを数式で表すとこうなるよ。

安全係数 $F_s = 1.25$（3水準以上のメッシュ比較時）。GCI < 5% を収束の目安とする。

解析結果の信頼性を担保するため、以下のベンチマーク問題との比較を推奨：

• マルチグリッド（AMG）前処理: 大規模問題のスケーラビリティ向上
• GPU並列化: 行列-ベクトル積のGPUオフロード
• ドメイン分割法: MPI並列による分散メモリ計算
• 縮約基底法（ROM）: パラメータスタディの高速化

うん、いい調子だよ！ 実際に手を動かしてみることが一番の勉強だからね。分からないことがあったらいつでも聞いてくれ。