射出成形繊維配向解析 — 理論と支配方程式

繊維強化樹脂(GF/CF)の射出成形における繊維配向テンソルの予測。Folgar-Tucker方程式による配向分布の計算と、異方性機械特性のマッピング。

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

射出成形繊維配向解析のシミュレーションは、熱力学・材料力学・流体力学の連成問題として定式化される。製造プロセスの物理現象は複数の時間・空間スケールにまたがるため、マクロスケールの連続体モデルとメゾ/ミクロスケールの材料モデルの適切な組合せが求められる。プロセスパラメータ（温度、速度、荷重等）と製品品質（寸法精度、欠陥、機械特性）の因果関係を定量的に予測することが目標なんだ。

製造プロセスシミュレーションの精度は材料モデルの忠実度に大きく依存する。弾塑性構成則、クリープ則、相変態モデルなどを温度・ひずみ速度の関数として適切に定義する必要がある。材料試験(引張、圧縮、ねじり)から得られたデータをフィッティングし、外挿範囲での妥当性を検証する。JMatProやThermo-Calcなどの熱力学データベースも活用する。

製造プロセスシミュレーションは、熱力学・流体力学・固体力学の連成問題として定式化される。

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

ここで $T$ は温度、$\mathbf{v}$ は材料の速度場、$k$ は熱伝導率、$Q$ は内部発熱（ジュール熱、潜熱、摩擦熱等）なんだ。

凝固過程では潜熱の放出/吸収が温度場に大きく影響する。エンタルピー法による定式化：

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

ここで $L$ は潜熱、$f_l(T)$ は液相率（固液共存域で0から1の間の値をとる）。

金属の塑性変形はJohnson-Cook構成則等で記述される：

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

$A$: 初期降伏応力、$B$: 硬化係数、$n$: 硬化指数、$C$: 歪み速度感度、$m$: 温度軟化指数。

溶融金属や樹脂の流動はナビエ-ストークス方程式に従うが、高粘性・非ニュートン流体特性を考慮する必要がある。射出成形ではCross-WLFモデルが標準的：

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

• 連続体力学の仮定が成立するスケール（粒子径 >> 分子間距離）
• 相変化の温度幅が十分大きい場合、mushy zone のモデル化が精度に影響
• 高速変形（衝撃鍛造等）では慣性効果の考慮が必要
• 微細組織予測にはフェーズフィールド法やセルラーオートマトンの追加が必要

解析対象の物理現象を支配する無次元パラメータの理解は、適切なモデル選択とパラメータ設定の基盤となる。

• ペクレ数 Pe: 対流と拡散の相対的重要性。Pe >> 1 で対流支配（安定化手法が必要）
• レイノルズ数 Re: 慣性力と粘性力の比。流体問題の基本パラメータ
• ビオ数 Bi: 内部伝導と表面対流の比。Bi < 0.1 で集中熱容量法が適用可能
• クーラン数 CFL: 数値安定性の指標。陽解法では CFL ≤ 1 が必要

解析結果のオーダー推定には、バッキンガムのΠ定理に基づく次元解析が効果的なんだ。代表長さ $L$、代表速度 $U$、代表時間 $T = L/U$ を用いて、各物理量のオーダーを事前に推定し、解析結果の妥当性を確認する。

適切な境界条件の選択は解の一意性と物理的妥当性に直結するんだよ。不足した境界条件は不適切な問題となり、過剰な境界条件は矛盾を生じさせる。

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