ガラス成形シミュレーション

ガラス成形の物理と理論基礎

粘度が10桁動く材料

ガラス成形シミュレーションの主役は粘度の温度依存性です。ソーダ石灰ガラスで、溶解域では水飴程度（10² Pa·s）、成形の作業域で10³〜10⁷ Pa·s、徐冷点で10¹² Pa·s、歪点で10¹³·⁵ Pa·s——プロセス全体で10桁以上変化します。温度依存はアレニウス型では表せず、VFT（Vogel-Fulcher-Tammann）式が標準です。

$$ \log_{10} \eta = A + \frac{B}{T - T_0} $$

この3定数（組成ごとに実測）がプロセス解析の物性の心臓部で、粘度曲線の誤差はそのまま肉厚分布・成形可否の誤差になります。基準点（軟化点=10⁶·⁶ Pa·s、徐冷点=10¹² Pa·s等）での粘度実測からVFT定数をフィットするのが実務の出発点です。

半透明媒体としての熱——内部放射

高温のガラスは近赤外を内部透過する半透明（参与）媒体で、冷却は表面からの伝導・対流だけでなく体積からの放射で進みます。光学的に厚い（肉厚が吸収長より十分大きい）場合はロスランド近似で「放射伝導率」として熱伝導に繰り込め、実効熱伝導率が高温で数倍に増える形になります。薄肉・低吸収の場合はこの近似が破れ、帯域放射輸送（DOM等）が必要——肉厚と波長帯で放射の扱いを切り替える判断が、ガラス熱解析の最初の分岐です。

ガラス転移は「速度過程」——構造緩和

数値解析の手法

成形段階——超高粘度・非等温の自由表面流

プレス・ブロー成形の充填過程は、粘度10³〜10⁷ Pa·sの非等温ストークス流（慣性はほぼ無視できる：Re≪1）として解きます。数値上の要点は、①自由表面と金型接触の追跡（ALE/界面追跡）、②温度-粘度の強連成（薄くなった部分が冷えて硬くなり、さらに伸びなくなる正帰還）、③金型との接触熱伝達——接触コンダクタンスは接触圧・離型剤・表面粗さ依存で、数百〜数千 W/m²Kの幅を持つ校正パラメータです。ブロー成形では圧力駆動の膜伸展として殻近似を使う実装もあり、肉厚分布の予測が主目的になります。

冷却・応力段階——粘弾性とTNMの二本柱

冷却過程の残留応力は2つの機構の和です。①熱粘弾性応力——表面が先に固化し内部が後から収縮する「焼き入れ」機構（強化ガラスはこれを意図的に使う）。応力緩和関数（Prony級数）と温度シフト（熱レオロジー単純性）で計算。②構造緩和応力——TNMモデルの仮想温度分布の差による密度差起因。光学素子の複屈折・寸法安定性はここが支配します。徐冷（アニール）スケジュール設計は「徐冷点付近をゆっくり通す」の定量化で、この2機構を解けば徐冷時間とロットの残留応力・複屈折が予測できます。

数値的な急所

• 粘度の桁変化 — 要素間で粘度が数桁違う系は行列が悪条件化。粘度の上限クリップ（固化域は弾性体へ切替）や領域の段階的な固体化処理が定石
• 接触熱伝達の不確かさ — 金型hは第一級の校正パラメータ。実測温度（金型埋込み熱電対）でのフィットを省略しない
• 放射モデルの選択 — 光学厚さで近似の適用可否を先に判定。薄物にロスランドを使うと冷却を過大評価
• 時間スケールの分離 — 成形は秒、徐冷は時間〜日。段階ごとに解析を分割し、界面でデータを渡す

実務ワークフロー

標準ワークフロー

1. 物性整備 — VFT粘度曲線・熱物性・吸収スペクトル（放射用）・緩和関数/TNMパラメータ。組成が変わったら全て取り直し
2. 成形解析 — ゴブ（素地）温度・金型温度・プレス速度/ブロー圧→充填・肉厚分布・表面欠陥リスク
3. 冷却解析 — 搬送・冷却条件→温度履歴（内部放射込み）
4. 応力・品質解析 — 残留応力・複屈折・仮想温度分布→徐冷条件の設計
5. 校正と検証 — 肉厚実測・応力（偏光計測）・温度実測でモデルを校正。条件変更の予測に展開

肉厚分布が最初の勝負

容器・管球のブロー成形では、肉厚分布の予測精度がシミュ導入価値の大半を占めます。肉厚を支配するのはゴブの初期温度分布と金型接触のタイミング——「先に触れた場所は冷えて伸びなくなる」の積み重ねです。実測肉厚マップとの比較でゴブ温度分布・接触hを校正すれば、金型形状変更・ゴブ形状変更の効果を試作前に順位付けできます。絶対精度±10%でも、条件間の優劣比較には十分という使い方は、他の成形シミュ（溶接等）と同じ実務哲学です。

品質不良と解析の対応表

ツールの選択肢

ツール比較

選定の軸

①主目的が成形（肉厚・充填）か品質（応力・複屈折・徐冷）か——前者はPolyflow系の流動、後者はAbaqus系の粘弾性+TNMが本流。②内部放射の精度要求（光学ガラスの精密冷却なら帯域放射対応）。③物性データの整備体制——どのツールでもVFT・緩和関数・吸収係数は自社材料での実測が前提で、ツール内蔵の「一般ソーダガラス」で精密議論はできません。

先端動向

精密ガラスモールド成形（GMP）

非球面レンズのモールド成形では、サブμmの形状転写精度が要求され、シミュレーションの焦点も「肉厚」から「金型形状補正量の予測」へ移ります。冷却収縮・構造緩和による屈折率分布・金型弾性変形まで含めた高精度連成で、金型補正の試行回数を減らす取り組みが実用段階です。TNMパラメータの高精度同定（示差走査熱量計＋ビーム曲げ実験）が精度の律速になっています。

溶解槽〜成形のプロセスチェーン連成

溶解槽（燃焼・電気ブースト・泡・均質化）のCFDと成形解析をつなぎ、ゴブの温度・組成の履歴を上流から供給するプロセスチェーンシミュレーションが進んでいます。省エネ（溶解は極めてエネルギー集約的）と品質（脈理・泡）の同時改善が動機で、燃焼の参与媒体放射・電気加熱のジュール熱・原料溶け込みまで含む大規模連成が研究・産業の前線です。

機械学習の活用

成形条件→肉厚分布のサロゲートによる金型設計探索、炉の運転データからの品質予測（ソフトセンサー）、徐冷スケジュールの最適化など、校正済み物理モデル＋データ駆動補正のグレーボックス構成が現実的な適用形です。物性（VFT・TNM）の逆同定に実測とベイズ推定を組み合わせる研究も、パラメータ不確かさの定量化（ベイズ較正）として整備されつつあります。

トラブル対応

症状別の原因と対策

導入立ち上げの順序

関連記事：産業熱プロセスの記事一覧、放射熱交換の基礎、リフロー炉シミュレーション。