玻璃成形模拟
玻璃成形的理论基础
概要
老师！今天我们要讨论玻璃成形模拟，这是什么呢？
玻璃粘性流动。吹塑成形、压制成形。VFT粘度模型。
原来如此。那么掌握了玻璃粘性流动就基本没问题了吧？
支配方程
离散化手法
这些方程在计算机上实际如何求解呢？
使用有限元法（FEM）进行空间离散化。组装单元刚度矩阵，构建全局刚度方程。
矩阵求解算法
矩阵求解算法具体是什么意思呢？
通过直接法（LU分解、Cholesky分解）或迭代法（CG法、GMRES法）求解联立方程组。对于大规模问题，预处理迭代法效果显著。
| 求解法 | 分类 | 内存使用 | 适用规模 |
|---|---|---|---|
| LU分解 | 直接法 | O(n²) | 小~中规模 |
| Cholesky分解 | 直接法（对称正定） | O(n²) | 小~中规模 |
| PCG法 | 迭代法 | O(n) | 大规模 |
| GMRES法 | 迭代法 | O(n·m) | 大规模·非对称 |
| AMG预处理 | 预处理 | O(n) | 超大规模 |
也就是说，有限元法那个地方要是偷懒了，后面会吃大亏啊。我一定会牢记在心！
商用工具中的实现
做玻璃成形模拟的话，可以用什么样的软件呢？
| 工具名称 | 开发商/现所有者 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
供应商系谱及产品整合经历
各个软件的历史背景好像挺有戏剧性的？
Ansys Mechanical (旧ANSYS Structural)
请给我讲讲「Ansys Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys Parametric Design Language）。
现所有权：Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是什么呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，并整合到SIMULIA品牌。
现所有权：Dassault Systèmes SIMULIA
听到这里，我终于明白为什么开发历史那么重要了！
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。以MATLAB联合的FEMLAB开始，后改名为COMSOL。在多物理场方面有优势。
现所有权：COMSOL AB
哇~，开发的历史故事超级有意思！请再多讲讲。
文件格式与互操作性
不同软件之间互相传输数据的时候有什么需要注意的吗？
| 格式 | 扩展名 | 类别 | 概要 |
|---|---|---|---|
| STEP | .stp/.step | 中立CAD | 符合ISO 10303的3D CAD数据交换格式。支持形状+PMI。 |
| IGES | .igs/.iges | 中立CAD | 早期的CAD数据交换规范。曲面数据兼容性有问题。逐步向STEP迁移。 |
在不同求解器间转换模型时，需要注意单元类型的对应关系、材料模型的兼容性、荷载和边界条件的表达差异。特别是对于高阶单元和特殊单元（内聚力单元、用户定义单元等），多数情况下无法在求解器间直接转换。
原来…格式看起来很简单，其实深度有这么深啊。
实务注意事项
教科书里没有的「现场经验」之类的有吗？
网格收敛性确认、边界条件的合理性验证、材料参数的敏感性分析非常重要。
玻璃成形模拟的全貌我搞清楚了！明天开始在实务中试试。
嗯，很好！实际动手做才是最好的学习。有不明白的地方随时来问我。
玻璃粘度与温度的关系
Vogel-Fulcher-Tammann（VFT）方程在1921~26年间由3位科学家独立发现，用来描述玻璃粘度与温度的关系。其形式为log(η) = A + B/(T-T₀)，对于钠钙玻璃，常数为A=-2.6、B=4,160、T₀=261。该公式已成为玻璃成形CAE中粘度输入的实际标准，ANSYS Polyflow的Glass Flow模块已默认集成。
玻璃成形的数值计算方法
数值方法详述
具体用什么算法来求解玻璃成形模拟呢？
离散化的定式化
使用形状函数 $N_i$ 来近似未知量：
这样就用数学式表示了。
基本方程的离散形式
这样就用数学式表示了。
嗯…只看公式的话，我有点不理解…这表示什么呢？
连续体的支配方程离散化后，可得如下代数方程组：
这里$[K]$是全局刚度矩阵（或等效系统矩阵），$\{u\}$是未知节点变量向量，$\{F\}$是外力向量。
啊，这样啊！连续体的支配方程离散化就是这么一回事啊。
单元技术
我听过「单元技术」这个说法，但好像没真正理解…
| 单元类型 | 阶数 | 节点数(3D) | 精度 | 计算成本 |
|---|---|---|---|---|
| 四面体1阶 | 线性 | 4 | 低（剪切锁定） | 低 |
| 四面体2阶 | 二次 | 10 | 高 | 中 |
| 六面体1阶 | 线性 | 8 | 中 | 中 |
| 六面体2阶 | 二次 | 20 | 非常高 | 高 |
| 棱柱 | 线性/二次 | 6/15 | 中~高 | 中 |
积分方案
积分方案具体是什么意思呢？
听到这里，我终于明白为什么单元类型那么重要了！
收敛性和稳定性
如果收敛不了的话，应该先检查什么？
收敛速度：对于二阶单元，误差以 $O(h^2)$ 的阶减少（光滑解的情况）
原来…网格细分看似简单，其实深度有这么深啊。
求解器设置的建议
具体用什么算法来求解玻璃成形模拟呢？
| 参数 | 推荐值 | 备注 |
|---|---|---|
| 迭代法收敛判定 | $10^{-6}$ | 残差范数标准 |
| 预处理方法 | ILU(0) 或 AMG | 按问题规模 |
| 最大迭代次数 | 1000 | 不收敛时需重新调整设置 |
| 内存模式 | In-core | 尽可能使用 |
线性单元 vs 2阶单元
在热传导分析中，线性单元往往能获得足够精度。温度梯度急剧的区域（热冲击等）推荐使用2阶单元。
热流量的评估
从单元内温度梯度算出。和节点应力一样，有时需要光滑化处理。
对流-扩散问题
当Peclet数高时（对流支配），需要迎风稳定化（SUPG等）。纯热传导问题不需要。
非定常分析的时间步长
热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$：热扩散率）相比应设置足够小的步长。快速温度变化可使用自适应时间步长控制。
非线性收敛
由温度相关物性导致的非线性通常较温和，采用Picard迭代（直接代入法）往往足够。放射的强非线性推荐使用Newton法。
定常分析的判定
当全部节点温度变化低于阈值（$|\Delta T| / T_{max} < 10^{-5}$等）时判定为收敛。
玻璃成形的实务应用
实务指南
老师，请给我讲讲「实务指南」！
讲解玻璃成形模拟的实务分析流程和注意事项。
分析流程
从最开始的第一步教我吧！应该从什么开始呢？
1. 预处理 (Pre-processing)
- 导入CAD数据并简化形状
- 定义材料特性
- 网格生成（确定单元类型和尺寸）
- 设置边界条件和荷载条件
2. 求解 (Solving)
- 设置求解器参数（求解方法、收敛基准、输出控制）
- 提交作业并执行计算
- 监控收敛过程
3. 后处理 (Post-processing)
- 结果可视化（位移、应力、其他物理量）
- 结果验证和合理性确认
- 报告编制
网格生成的最佳实践
怎样判断网格的好坏呢？
单元品质指标
请给我讲讲「单元品质指标」！
| 指标 | 理想值 | 允许范围 | 影响 |
|---|---|---|---|
| 纵横比 | 1.0 | < 5.0 | 精度降低 |
| Jacobian比 | 1.0 | > 0.3 | 单元退化 |
| 翘曲 | 0° | < 15° | 精度降低 |
| 歪斜度 | 0° | < 45° | 收敛性恶化 |
| 锥形比 | 0 | < 0.5 | 精度降低 |
网格密度的确定
网格密度的确定具体是什么意思呢？
边界条件的设置指南
我听说边界条件要是弄错了，全部都完蛋了…
啊，这样啊！过度约束的注意就是这么一回事啊。
按商用工具的实现步骤
有那么多软件吧？各自有什么特点？请告诉我！
| 工具名称 | 开发商/现所有者 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
Ansys Mechanical (旧ANSYS Structural)
请给我讲讲「Ansys Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys Parametric Design Language）。
现所有权：Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是什么呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，并整合到SIMULIA品牌。
现所有权：Dassault Systèmes SIMULIA
听到这里，我终于明白开发历史为什么那么重要了！
常见失误和对策
初学者容易犯什么样的错误？想提前知道！
| 症状 | 原因 | 对策 |
|---|---|---|
| 计算不收敛 | 网格品质不良、不适当的边界条件 | 网格改善、约束条件检查 |
| 应力异常大 | 应力奇点、网格依存 | 避免奇点、局部网格细分 |
| 位移不现实 | 材料常数错误、单位系不一致 | 确认输入数据 |
| 计算时间过长 | 不必要的细分、非高效求解法 | 网格优化、并行计算 |
质量保证检查清单
教科书里没有的「现场经验」之类的有吗？
玻璃成形模拟的全貌我搞清楚了！明天开始在实务中试试。
嗯，很好！实际动手做才是最好的学习。有不明白的地方随时来问我。
旭硝子（AGC）的玻璃成形优化
AGC（原旭硝子）在2009年利用Polyflow优化了智能手机盖板玻璃「Dragonfly（后改名为Dragonfly Eye）」的成形工艺。在650~720℃压制成形温度范围内计算了50个条件，实现了厚度偏差±0.05mm以内，这项成就成为获得iPhone 3GS(2009年)供应合同的技术背景，后来技术报告中有明确提及。
玻璃成形的软件对比
商用工具对比
有那么多软件吧？各自有什么特点？请告诉我！
讲解支持玻璃成形模拟的主要商用CAE工具的功能对比和各产品的历史背景。
支持工具列表
做玻璃成形模拟的话，可以用什么样的软件呢？
| 工具名称 | 开发商/现所有者 | 主要文件格式 |
|---|---|---|
| Ansys Mechanical (旧ANSYS Structural) | Ansys Inc. | .cdb, .rst, .db, .ans, .mac |
| Abaqus FEA (SIMULIA) | Dassault Systèmes SIMULIA | .inp, .odb, .cae, .sta, .msg |
| COMSOL Multiphysics | COMSOL AB | .mph |
| Ansys Fluent | Ansys Inc. | .cas, .dat, .msh, .jou |
Ansys Mechanical (旧ANSYS Structural)
请给我讲讲「Ansys Mechanical」！
1970年由Swanson Analysis Systems Inc. (SASI) 开发。基于APDL（Ansys Parametric Design Language）。
现所有权：Ansys Inc.
Abaqus FEA (SIMULIA)
Abaqus FEA具体是什么呢？
1978年由HKS (Hibbitt, Karlsson & Sorensen) 开发。2005年被Dassault Systèmes收购，并整合到SIMULIA品牌。
现所有权：Dassault Systèmes SIMULIA
听到这里，我终于明白开发历史为什么那么重要了！
COMSOL Multiphysics
请给我讲讲「COMSOL Multiphysics」！
1986年在瑞典成立。以MATLAB联合的FEMLAB开始，后改名为COMSOL。在多物理场方面有优势。
现所有权：COMSOL AB
Ansys Fluent
现在轮到Ansys Fluent了。有什么内容吗？
Fluent Inc.开发。2006年被Ansys收购。非结构网格的通用CFD求解器。
现所有权：Ansys Inc.
啊，这样啊！开发历史就是这么一回事啊。
功能对比矩阵
预算和时间都有限，最具性价比的是哪个？
| 功能 | Ansys Mechanical | Abaqus | COMSOL | Fluent |
|---|---|---|---|---|
| 基本功能 | ○ | ○ | ○ | ○ |
| 高级功能 | ○ | ○ | ○ | △ |
| 自动化/脚本 | ○ | ○ | ○ | ○ |
| 并行计算 | ○ | ○ | ○ | ○ |
| GPU支持 | △ | △ | △ | ○ |
转换时的风险
转换时的风险具体是什么意思呢？
啊，这样啊！不同工具间的转换就是这么一回事啊。
许可证形式
我听过「许可证形式」，但好像没真正理解…
| 工具 | 许可证 | 特点 |
|---|---|---|
| 商用FEA | 节点锁定/浮动 | 价格高但有官方支持 |
| OpenFOAM | GPL | 免费但支持有偿 |
| COMSOL | 节点锁定/浮动 | 按模块购买 |
| Code_Aster | GPL | EDF开发的开源求解器 |
选择指南
最后到底选哪个呢？给我一个判断标准吧？
在玻璃成形模拟工具的选择中，应考虑以下几点：
玻璃成形模拟的全貌我搞清楚了！明天开始在实务中试试。
嗯，很好！实际动手做才是最好的学习。有不明白的地方随时来问我。
CPFD Software公司的玻璃炉分析
美国CPFD Software（2004年成立）的Barracuda VR是针对流化床、玻璃熔融炉专化的CFD求解器，采用MP-PIC法处理固液气三相。在AGC的板玻璃熔融炉(1,500t/天热容量)应用中，炉内温度均一性的计算精度较既有CFD提高了25%，直接实现了燃料消耗量优化。从2020年起以SaaS形式通过云提供。
玻璃成形的先进研究
前沿话题和研究动向
玻璃成形模拟这个领域，今后会怎样发展呢？
了解玻璃成形模拟最新研究动向和先进方法。
最新数值方法
下面是最新数值方法的话题吧。有什么内容？
嗯…只看公式的话，我有点不理解…这表示什么呢？
高性能计算(HPC)的支持
| 并行化方法 | 概要 | 适用求解器 |
|---|---|---|
| MPI (区域分割) | 分布式内存型。大规模问题的标准 | 全主要求解器 |
| OpenMP | 共享内存型。节点内并行 | 多数求解器 |
| GPU (CUDA/OpenCL) | GPGPU活用。对显式法特别有效 | LS-DYNA、Fluent等 |
| 混合 MPI+OpenMP | 节点间+节点内并行 | 大规模HPC环境 |
玻璃成形的故障处理
故障排查
常见错误和对策
老师也因为玻璃成形模拟彻夜调试过吗？（笑）
1. 收敛失败
收敛失败具体是什么意思呢？
症状：求解器在指定迭代次数内无法收敛异常终止
可能原因：
- 网格品质不足（过度歪斜单元）
- 材料参数设置不当
- 初始条件不当
- 非线性过强（荷载步数不足）
对策：
- 进行网格品质检查（纵横比、Jacobian）
- 确认材料参数单位制
- 将荷载分为多个步骤（增加子步数）
- 放松收敛判定基准（注意精度）
也就是说，收敛失败那个地方要是偷懒了，后面会吃大亏啊。我一定会牢记在心！
2. 非物理的结果
下面是非物理的结果的话题吧。有什么内容？
症状：应力/位移/温度等出现物理上不现实的数值
可能原因：
- 边界条件设置错误
- 单位制混乱（SI单位与工程单位混淆）
- 不当的单元类型选择
- 应力奇点的存在
对策：
- 确认反力总和（力的平衡）
- 确认单位制一致性
- 重新检查单元类型的适当性
- 消除奇点或进行子建模
前辈说「收敛失败这块一定要搞好」，我现在总算理解意思了。
3. 计算时间超过
计算时间超过具体是什么意思呢？
症状：计算耗时远超预期
对策：
- 优化网格粗密分布
- 活用对称性（1/2、1/4模型）
- 优化求解器设置（迭代法、预处理选择）
- 利用并行计算
4. 内存不足
请给我讲讲「内存不足」！
症状：Out of Memory 错误
前辈说「收敛失败这块一定要搞好」，我现在总算理解意思了。
对策：
- 使用核外求解法
- 减小网格规模
- 确认64bit版求解器的使用
- 增加内存分配
哇~，收敛失败的话题超级有意思！请再多讲讲。
Nastran代表错误
代表错误具体是什么意思呢？