铸造缺陷预测
铸造缺陷预测的理论基础
概要
老师！ 今天是关于铸造缺陷预测的内容吧？ 具体是什么？
缩孔(shrinkage porosity)、气孔、热裂纹(高温裂纹)等铸造缺陷的仿真预测。Niyama准则、Hot Spot分析、HTI(Hot Tearing Indicator)的定量评估。
等等，缩孔是指，也就是说这样的情况也能用到吗？
支配方程
用式子来表示就是这样的。
嗯，只看式子还是很难理解… 能说说是代表什么吗？
由Darcy流造成的补充流动：
理论的基础
我听说过「理论的基础」，但可能没有真正理解…
铸造缺陷预测的仿真是作为热力学、材料力学、流体力学的耦合问题来定义的。制造工艺的物理现象涉及多个时间和空间尺度，因此需要宏观尺度的连续体模型和中观/微观尺度材料模型的恰当组合。定量预测工艺参数（温度、速度、荷重等）与产品质量（尺寸精度、缺陷、机械特性）的因果关系是目标。
制造工艺的支配方程
我不太擅长公式…能教我铸造缺陷预测的式子的「含义」吗？
制造工艺仿真被定义为热力学、流体力学、固体力学的耦合问题。
热传导方程（能量守恒）
热传导方程是什么？
这里 $T$ 是温度，$\mathbf{v}$ 是材料的速度场，$k$ 是热导率，$Q$ 是内热（焦耳热、潜热、摩擦热等）。
我理解了师兄说「制造工艺仿真一定要好好做」的含义。
凝固·相变
「凝固·相变」请告诉我！
在凝固过程中，潜热的释放/吸收对温度场产生很大影响。基于焓方法的定式化：
用式子来表示就是这样的。
嗯，只看式子还是很难理解… 能说说是代表什么吗？
这里 $L$ 是潜热，$f_l(T)$ 是液相率（在固液共存域中取0到1之间的值）。
塑性变形的本构关系
塑性变形的本构关系是什么？
金属的塑性变形用Johnson-Cook本构关系等来表述：
$A$: 初始屈服应力，$B$: 硬化系数，$n$: 硬化指数，$C$: 应变速率敏感系数，$m$: 温度软化指数。
现在我理解了为什么制造工艺仿真很重要！
流动解析（充填·铸造）
接下来是流动解析的话题。内容是什么？
熔融金属或树脂的流动遵循Navier-Stokes方程，但需要考虑高粘度和非牛顿流体特性。注射成形中Cross-WLF模型是标准的：
确实，制造工艺仿真看似简单，实际上非常深邃。
假设和适用极限
这个式子不是万能的吧？ 不能用的时候有哪些？
无量纲参数和支配性尺度
老师，「无量纲参数和支配性尺度」请告诉我！
理解支配物理现象的无量纲参数是选择恰当模型和参数设置的基础。
哎呀，原来是这样！ 分析对象的物理现象就是这个原理。
量纲分析的验证
「量纲分析的验证」请告诉我！
基于Buckingham Π定理的量纲分析是分析结果数量级估计的有效方法。使用代表长度 $L$、代表速度 $U$、代表时间 $T = L/U$，提前估计各物理量的数量级，确认分析结果的合理性。
好吧，那么只要分析对象的物理现象能做好，基本上就没问题了？
边界条件的分类和数学特征
边界条件，听说这个搞错了就全毁了…
| 类型 | 数学表达 | 物理意义 | 示例 |
|---|---|---|---|
| Dirichlet条件 | $u = u_0$ on $\Gamma_D$ | 变量值的指定 | 固定墙、指定温度 |
| Neumann条件 | $\partial u/\partial n = g$ on $\Gamma_N$ | 梯度（流量）的指定 | 热流束、力 |
| Robin条件 | $\alpha u + \beta \partial u/\partial n = h$ | 变量和梯度的线性组合 | 对流热传达 |
| 周期边界条件 | $u(x) = u(x+L)$ | 空间周期性 | 单位胞分析 |
适当的边界条件选择与解的唯一性和物理合理性直接相关。不足的边界条件会导致不适定问题，过度的边界条件会产生矛盾。
铸造缺陷预测的整体情况我都理解了！ 从明天开始在实务中意识到这些。
好的，进展很顺利！ 实际动手是最好的学习方法。有问题的时候随时问我。
树枝晶创造的「堵塞」——凝固收缩的真面目
金属凝固时，树枝状结晶（树晶）首先成长形成骨架。这个骨架的空隙被熔融金属试图补给，但树晶之间相互缠绕形成「堵塞」。补给不足的空间变成缩孔。铝合金的凝固收缩率约为6%，这个体积分的熔湯无法补给就会产生缺陷。这个物理用Darcy准则来描述是凝固收缩仿真的本质，树晶的渗透率（Permeability）精度大大影响预测的好坏。
铸造缺陷预测的数值计算方法
数值方法的详细说明
具体用什么算法来解铸造缺陷预测呢？
解释铸造缺陷预测仿真采用的数值方法。
离散化方法
大变形的制造工艺一般用Updated Lagrangian法或ALE(Arbitrary Lagrangian-Eulerian)法。接触问题用罚函数法或Lagrange乘数法。Euler法用于鍛造、挤出等定常工艺流场定式。
时间积分
老师，「时间积分」请告诉我！
准静态问题用隐式法(Newton-Raphson)，高速变形、冲击问题用显式法(中心差分法)。通过质量缩放可以缓解显式法的时间步长限制，但需要监视运动能量不超过内部能的5-10%。
网格管理
网格越细越好对吧？ …不对吗？
针对大变形引起的网格失真，使用重新网格化(r-adaptivity)或ALE网格光滑处理。SPH法或MPM(Material Point Method)等无网格方法也是选项。
等等，大变形引起的网格是指，也就是说这样的情况也能用到吗？
接触与摩擦建模
「接触与摩擦建模」我听说过，但可能没有真正理解…
制造工艺中工具和被加工材的接触不可避免，接触算法的选择影响解的精度和稳定性。根据工序采用Coulomb摩擦、剪切摩擦、温度相关摩擦模型。接触检测的罚函数参数或分段对分段法的设置对计算稳定性有很大影响。
我理解了师兄说「制造工艺的工具一定要好好做」的意思。
数值求解的实现细节
老师，「数值求解的实现细节」请告诉我！
网格需求
网格需求是什么？
制造工艺仿真需要追踪移动界面（固液界面、自由表面），网格策略非常重要。
| 方法 | 概要 | 适用 |
|---|---|---|
| ALE法 | 网格随材料一起运动 | 鍛造、压延 |
| 欧拉法 | 材料在固定网格上流动 | 铸造充填 |
| VOF法 | 用体积分数追踪自由表面 | 铸造、注塑成形 |
| CEL法 | 耦合欧拉-Lagrange | 冲击加工 |
| SPH法 | 粒子法，无网格 | AM熔融池 |
热源模型（焊接·AM）
热源模型是什么？
Goldak二重椭球体模型：
这里 $P$ 是激光/电弧功率，$\eta$ 是吸收效率，$a,b,c$ 是椭球体的半轴长。
也就是说，网格需求的地方要是不仔细的话，后面会很痛啊。铭记在心！
时间积分
时间积分是什么？
也就是说，网格需求的地方要是不仔细的话，后面会很痛啊。铭记在心！
耦合求解器策略
接下来是耦合求解器策略。什么内容？
哎呀，原来是这样！ 网格需求就是这个原理。
误差评估和精度验证
「误差评估和精度验证」我听说过，但可能没有真正理解…
离散化误差的评估
离散化误差的评估是什么？
Richardson外推法估计离散化误差：
这里 $f_h$ 是网格宽度 $h$ 的解，$r$ 是网格比，$p$ 是离散化的阶。
GCI（Grid Convergence Index）
「GCI」请告诉我！
基于ASME V&V 20-2009的网格收敛性定量评价：
现在我理解了为什么离散化误差的评估很重要！
用式子来表示就是这样的。
嗯，只看式子还是很难理解… 能说说是代表什么吗？
安全系数 $F_s = 1.25$（3水准以上网格比较时）。GCI < 5% 作为收敛的参考。
我理解了师兄说「离散化误差的评估一定要好好做」的意思。
验证基准问题
「验证基准问题」请告诉我！
为了保证分析结果的可信度，推荐与以下基准问题进行比较：
| 领域 | 基准 | 参考解 |
|---|---|---|
| 结构 | 单元补丁测试 | 一致应力场的重现 |
| 结构 | Scordelis-Lo屋顶 | 参考位移 |
| 流体 | 盖驱动型空腔 | Ghia et al. (1982) |
| 热 | 1D解析解 | $T(x) = T_0 + (T_1-T_0)x/L$ |
加速方法
老师，「加速方法」请告诉我！
铸造缺陷预测的整体情况我都理解了！ 从明天开始在实务中意识到这些。
好的，进展很顺利！ 实际动手是最好的学习方法。有问题的时候随时问我。
Niyama准则值——蕴含在简洁式子里的现场智慧
Niyama准则值用一个简单的式子 G/√Ṫ（温度梯度÷冷却速度的平方根）评估缩孔风险。1976年日本的二山氏提出这个指标，把复杂的凝固物理全部跳过，简化为「梯度不足且冷却快的地方危险」的现场经验法则进行数式化。现代仿真器要解Navier-Stokes和Darcy定律，但最终的缺陷判定仍在使用这个近50年前的经典经验式。简洁的智慧能够活到现在，这很有意思。
铸造缺陷预测的实务应用
实践指南
老师，「实践指南」请告诉我！
解释铸造缺陷预测的实务分析步骤和最佳实践。
我理解了师兄说「铸造缺陷预测的实务就是…」的意思。
分析流程
请从最开始的一步教我！ 应该从什么开始？
1. 工艺条件的定义: 整理和范围设定工艺参数（温度、速度、荷重、时间）
2. 材料数据的准备: 从试验数据同定温度、应变速率相关的本构关系参数
4. 工艺仿真实行: 逐步增加复杂度并确认收敛性
5. 结果验证: 与试验数据比较（尺寸精度、荷重历程、温度分布、缺陷位置）
哎呀，原来是这样！ 工艺条件的定义就是这个原理。
最佳实践
老师，「最佳实践」请告诉我！
品质管理和文档化
有没有教科书上没有的「现场智慧」这样的东西呢？
把分析条件（材料数据出处、边界条件的根据、网格设定的妥当性）系统地文档化。建立分析结果的审查流程，进行试验比较的精度评价并定量记录。还要定期进行NAFEMS等基准问题的求解器验证。
实务的分析步骤
在实务中使用铸造缺陷预测时，最应该注意的是什么？
铸造仿真的工作流
铸造仿真的工作流是什么？
1. CAD模型准备: 产品形状 + 湯道系 + 冒口 + 冷铁的3D模型
2. 网格生成: 推荐以六面体占主导的网格。薄壁部分最少3层以上
4. 边界条件: 铸型-金属间的热传达系数（IHTC）。型温的初期设定
5. 充填解析: 设定注湯速度·温度。監視空气卷入
6. 凝固解析: 充填完成后的温度场解析。缩孔预测
7. 应力解析: 凝固后的残余应力、脱模后的变形
注塑成形仿真的参数设定
接下来是注塑成形仿真的参数的话题。什么内容？
| 参数 | 典型值 | 影响 |
|---|---|---|
| 树脂温度 | 200-300°C | 流动性、表面品质 |
| 模具温度 | 40-100°C | 冷却时间、结晶度 |
| 注射速度 | 50-200 mm/s | 浇口压、剪切应力 |
| 保压 | 50-100 MPa | 收缩补偿、尺寸精度 |
| 冷却时间 | 10-60 s | 生产性、翘曲变形 |
老师的说明很清楚！ 铸造仿真的疑惑消散了。
AM（积层造形）仿真的注意点
接下来是积层造形的话题。什么内容？
品质保证检查表
品质保证检查表是什么？
哎呀，原来是这样！ 铸造仿真就是这个原理。
项目管理和工作流自动化
想要粗略把握整个流程，能按步骤教我吗？
目录结构的推荐
接下来是目录结构的推荐的话题。什么内容？
```
project/
├── cad/ # CAD模型
├── mesh/ # 网格文件
├── setup/ # 分析设定文件
├── results/ # 计算结果
│ ├── case01/
│ ├── case02/
│ └── ...
├── postprocess/ # 后处理脚本·图像
├── report/ # 报告
└── validation/ # 验证数据
```
自动化脚本的活用
接下来是自动化脚本的活用的话题。什么内容？
可以用Python脚本把参数化研究和网格收敛性确认自动化，大幅提高可再现性和效率。
好吧，那么目录结构的推荐能做好的话，基本上就没问题了？
审查检查表
「审查检查表」请告诉我！
1. 输入数据: 材料常数的单位系，CAD的尺寸精度，网格品质指标
2. 边界条件: 物理合理性，过约束/约束不足的检查
3. 求解器设定: 收敛判定基准，时间步，输出频率
4. 结果验证: 力的平衡，能量守恒，与理论解的比较
5. 敏感性分析: 网格依赖性，边界条件的影响，材料参数的不确定性
也就是说，目录结构的推荐的地方要是不仔细的话，后面会很痛啊。铭记在心！
报告书制作的要点
老师，「报告书制作的要点」请告诉我！
铸造缺陷预测的整体情况我都理解了！ 从明天开始在实务中意识到这些。
好的，进展很顺利！ 实际动手是最好的学习方法。有问题的时候随时问我。
「不要砸模具」——仿真机曾被铸造现场讨厌的年代
2000年代初期，仿真软件开始进入铸造现场，资深职人中反对「电脑会毁模具」的声音很大。这些失败案例是通过推崇仿真结果改变浇道位置后缺陷反而增加的事例口口相传传开的。大多数原因是材料物性值输入错误，铸铁的固相线温度只要搞错10℃就会让缩孔预测偏离很大。输入支持工具充实后「仿真比现场」的文化才得以改变。
铸造缺陷预测的软件比较
商用工具比较
有各种各样的软件吧？ 分别有什么特点呢？
比较铸造缺陷预测对应的主要商用仿真工具。
主要工具
有各种各样的软件吧？ 分别有什么特点呢？
| 工具 | 供应商 | 优点 |
|---|---|---|
| MAGMASOFT | MAGMA | 铸造工艺全般的统合解析 |
| Moldflow | Autodesk | 注塑成形的业界标准工具 |
| Simufact | Hexagon | 焊接、AM、塑性加工的统合 |
| DEFORM | SFTC | 锻造、压延实績丰富 |
| AutoForm | AutoForm | 板金成形的高速解析特化 |
| PAM-STAMP | ESI | 冲压成形的详细解析 |
| Amphyon/Netfabb | Oqton/Autodesk | AM向的工艺最优化 |
| ProCAST | ESI | 铸造的高精度耦合解析 |
选择标准
结果到底选哪个，告诉我判断基准？
评价对象工艺的专业性，材料数据库的充实度，与现有CAD/PLM的统合性，技术支持的质量。推荐用试用许可事前验证。
老师的说明很清楚！ 对象工艺的专业性疑惑消散了。
商用工具比较矩阵
那么，做铸造缺陷预测的话，有什么软件能用？
铸造仿真
铸造仿真是什么？
| 工具 | 开发方 | 主要功能 | 特点 |
|---|---|---|---|
| MAGMASOFT | MAGMA | 充填·凝固·应力·组织 | 铸造专用，世界份额No.1 |
| ProCAST | ESI Group | 充填·凝固·电磁搅拌 | 多物理场对应 |
| FLOW-3D CAST | Flow Science | 自由表面流动 | 用VOF法高精度充填解析 |
注塑成形仿真
接下来是注塑成形仿真的话题。什么内容？
| 工具 | 开发方 | 主要功能 |
|---|---|---|
| Moldflow | Autodesk | 充填·保压·冷却·翘曲·纤维取向 |
| Moldex3D | CoreTech | 真正的3D解析，IC封装对应 |
| Sigmasoft | SIGMA | 虚拟DOE、多周期解析 |
焊接·AM 仿真
焊接是什么？
| 工具 | 开发方 | 主要功能 |
|---|---|---|
| Simufact Welding | Hexagon | 焊接变形·残余应力 |
| Ansys Additive | Ansys | L-PBF/DED热力学解析 |
| Amphyon | Additive Works | AM翘曲补偿 |
| Virfac | Geonx | 焊接·AM热力学 |
塑性加工
接下来是塑性加工的话题。什么内容？
| 工具 | 开发方 | 主要功能 |
|---|---|---|
| AutoForm | AutoForm | 冲压成形、板成形 |
| DEFORM | Scientific Forming | 锻造、挤出、压延 |
| LS-DYNA | Ansys/LST | 冲击、板成形、通用显式法 |
| FORGE | Transvalor | 锻造、压延 |
许可形态和总拥有成本（TCO）
接下来是「许可形态和总拥有成本（TCO）」！ 什么内容？
商用工具的成本结构
商用工具的成本结构是什么？
| 项目 | 年额目安 | 备考 |
|---|---|---|
| 节点锁定许可 | 100-500万日元 | 1台PC固定 |
| 浮动许可 | 150-800万日元 | 网络内共享 |
| HPC令牌 | 50-300万日元 | 按并行核数的从量制 |
| 支持·保养 | 许可的15-25% | 包括版本升级 |
| 培训 | 30-80万日元/课程 | 初期导入必须 |
TCO比较的要点
比较的要点是什么？
供应商的技术支持比较
「供应商的技术支持比较」请告诉我！
导入流程和迁移策略
接下来是「导入流程和迁移策略」！ 什么内容？
供应商选定的步骤
「供应商选定的步骤」请告诉我！
1. 需求定义: 必要的解析功能、规模、精度需求明确化
2. 候选表清单作成: 3-5社进行筛选
3. 基准评价: 用各工具解析自社的典型问题
4. TCO算出: 5年间总拥有成本（许可+HPC+教育+支持）
5. PoC（概念验证）: 实业务的试用期间（3-6个月）
6. 最终选定: 技术评价+成本+支持+将来性的综合评价
工具迁移时的注意点
「工具迁移时的注意点」请告诉我！
铸造缺陷预测的整体情况我都理解了！ 从明天开始在实务中意识到这些。
好的，进展很顺利！ 实际动手是最好的学习方法。有问题的时候随时问我。
ProCAST vs MAGMASOFT——缺陷预测的「作法」的区别
铸造缺陷分析的两大工具ProCAST（ESI）和MAGMASOFT（MAGMA）虽然都解析同一部件，但缺陷图的细部会有所不同。原因是Niyama准则值的计算时机和补间方法的差异。MAGMASOFT倾向于在凝固末期的特定温度区间评估，ProCAST则倾向于按时间步累积。没有什么「哪个对」的说法，更多是「习惯该工具的作法」在现场准确度更高，因为相同工具积累的合否判定阈值可以再利用。中途更换工具的话过去数据的整合性会崩。
铸造缺陷预测的先进研究
先进话题
铸造缺陷预测领域，今后怎么发展呢？
讲述铸造缺陷预测领域的最新研究动向和今后展望。
等等，铸造缺陷预测领域是指，也就是说这样的情况也能用到吗？
最新研究动向
铸造缺陷预测领域，今后怎么发展呢？
数字孪生通过实时工艺监视·控制在快速进展。传感器内部数据（热成像、AE传感、力传感等）和仿真的融合，实现了制造中的品质预测和自适应制御。
国际标准化和标准符合
接下来是「国际标准化和标准符合」！ 什么内容？
制造工艺仿真的结果用于品质保证时，需要符合ISO、ASTM、JIS等相关规范。面向仿真结果的认证(Certification by Analysis)，国际的框架整备在进展，V&V指南(