玻璃成形模拟

玻璃成形的物理与理论基础

粘度跨十个数量级的材料

玻璃成形仿真的主角是粘度的温度依赖性。以钠钙玻璃为例，熔化区大致相当于糖稀（10² Pa·s），成形的作业区为10³～10⁷ Pa·s，退火点为10¹² Pa·s，应变点为10¹³·⁵ Pa·s——整个工艺跨越十个数量级以上。这种温度依赖性无法用阿伦尼乌斯型表达，标准做法是VFT（Vogel-Fulcher-Tammann）式。

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

这三个常数（按组成实测）是工艺分析中物性的心脏，粘度曲线的误差会原样变成壁厚分布与能否成形的误差。在基准点（软化点=10⁶·⁶ Pa·s、退火点=10¹² Pa·s等）实测粘度并拟合出VFT常数，是实务的出发点。

作为半透明介质的热——内部辐射

高温玻璃是内部可透过近红外的半透明（参与性）介质，冷却不仅靠表面的导热与对流，还靠来自体积内部的辐射推进。当光学厚度大（壁厚远大于吸收长度）时，可用罗斯兰德近似把它折算成"辐射导热系数"并入热传导，表现为有效导热系数在高温下增大数倍。薄壁或低吸收时该近似失效，需要分波段辐射输运（DOM等）——按壁厚与波段切换辐射的处理方式，正是玻璃热分析的第一个分岔点。

玻璃化转变是"速率过程"——结构弛豫

数值分析的方法

成形阶段——超高粘度·非等温的自由表面流

压制与吹制成形的充填过程，按粘度10³～10⁷ Pa·s的非等温斯托克斯流求解（惯性基本可忽略：Re≪1）。数值上的要点是：①自由表面与模具接触的追踪（ALE／界面追踪）；②温度-粘度的强耦合（变薄的部位冷却变硬、于是更拉不动的正反馈）；③与模具之间的接触传热——接触热导取决于接触压力、脱模剂与表面粗糙度，是一个跨越数百至数千 W/m²K 的标定参数。吹制成形中也有把它当作压力驱动的膜片伸展、采用壳近似的实现，此时主要目的就是预测壁厚分布。

冷却·应力阶段——粘弹性与TNM两根支柱

冷却过程的残余应力是两种机构之和。①热粘弹性应力——表面先固化、内部随后收缩的"淬火"机构（钢化玻璃就是有意利用它）。用应力松弛函数（Prony级数）与温度移位（热流变简单性）计算。②结构弛豫应力——源自TNM模型中虚拟温度分布之差带来的密度差。光学元件的双折射与尺寸稳定性由这一项支配。退火工艺的设计就是把"在退火点附近慢慢通过"定量化，解出这两种机构，就能预测退火时间以及批次的残余应力与双折射。

数值上的要害

• 粘度的数量级变化——单元之间粘度相差数个数量级的体系会使矩阵病态。粘度上限截断（固化区切换为弹性体）与区域的分阶段固化处理是定式
• 接触传热的不确定度——模具h是第一等的标定参数。不要省掉用实测温度（模具内埋热电偶）做拟合这一步
• 辐射模型的选择——先用光学厚度判定近似是否适用。对薄件套用罗斯兰德近似会高估冷却
• 时间尺度的分离——成形是秒级，退火是小时到天级。要按阶段拆分分析，在接口处传递数据

实务工作流程

标准工作流程

1. 物性整备——VFT粘度曲线、热物性、吸收光谱（供辐射用）、松弛函数/TNM参数。组成一变就要全部重新测
2. 成形分析——料滴（gob）温度·模具温度·压制速度/吹制压力→充填、壁厚分布、表面缺陷风险
3. 冷却分析——输送与冷却条件→温度历程（含内部辐射）
4. 应力·质量分析——残余应力·双折射·虚拟温度分布→退火条件的设计
5. 标定与验证——用实测壁厚、应力（偏光测量）、实测温度标定模型。再推广到条件变更的预测

壁厚分布是第一战场

在容器与管泡的吹制成形中，壁厚分布的预测精度占了引入仿真价值的大半。支配壁厚的是料滴的初始温度分布与模具接触的时机——"先碰到的地方就冷下来、拉不动了"的层层累积。用实测壁厚图去标定料滴温度分布与接触h，就能在开模之前把模具形状变更、料滴形状变更的效果排出优劣顺序。即便绝对精度只有±10%，对条件之间的优劣比较也够用——这种用法与其他成形仿真（如焊接）的实务哲学是一样的。

质量缺陷与分析量的对照表

工具的选择

工具比较

选型的着眼点

①主要目的是成形（壁厚·充填）还是质量（应力·双折射·退火）——前者以Polyflow系的流动为主流，后者以Abaqus系的粘弹性+TNM为主流。②对内部辐射的精度要求（光学玻璃的精密冷却就需要分波段辐射）。③物性数据的整备体制——无论哪种工具，VFT·松弛函数·吸收系数都以自家材料的实测为前提，靠工具内置的"通用钠钙玻璃"是无法做精密讨论的。

前沿动向

精密玻璃模压成形（GMP）

非球面透镜的模压成形要求亚微米级的形状转写精度，仿真的焦点也从"壁厚"转向"预测模具形状的补偿量"。把冷却收缩、结构弛豫带来的折射率分布乃至模具弹性变形都纳入的高精度耦合分析，已进入实用阶段，用于减少模具修正的试验次数。而TNM参数的高精度辨识（差示扫描量热＋梁弯曲实验）成了精度的瓶颈。

熔窑～成形的工艺链耦合

把熔窑（燃烧·电助熔·气泡·均化）的CFD与成形分析连起来、从上游供给料滴的温度与组成历程的工艺链仿真正在推进。动机是节能（熔化是极其耗能的环节）与质量（条纹·气泡）的同时改善；包含燃烧的参与性介质辐射、电加热的焦耳热乃至配合料溶入在内的大规模耦合，是研究与产业的最前线。

机器学习的应用

用"成形条件→壁厚分布"的代理模型做模具设计探索、从炉子的运行数据做质量预测（软测量）、优化退火工艺等，标定过的物理模型＋数据驱动修正的灰箱结构是现实的应用形态。把实测与贝叶斯推断结合起来做物性（VFT·TNM）逆辨识的研究，也正作为参数不确定度的定量化（贝叶斯标定）逐步成形。

故障排查

按症状的原因与对策

导入起步的顺序

相关文章：工业热过程文章一览、辐射换热基础、回流焊炉仿真。