第4章仿真理论基础

第一章仿真理论基础：CFD、传热、凝固与应力

各位同行，欢迎来到铸造仿真学习的第一站。说实话，很多人一上来就急着学软件操作，结果算出来的结果自己都不敢信。为什么？因为底子没打好。今天我们就来聊聊铸造仿真背后的四大理论支柱：流体力学、传热学、凝固理论和应力分析。

我个人习惯，学任何仿真工具之前，先花30%的时间理解原理。这就像开车前先看仪表盘——你知道每个数字代表什么，才不会开翻车。

核心观点：铸造仿真不是黑箱游戏。你输入的每个参数，软件计算的每个方程，背后都有物理意义。理解这些，你才能判断结果对不对。

说白了，CFD就是解流体运动的方程。铸造里最常见的是充型过程——金属液怎么流进型腔，有没有卷气，有没有冷隔。这些都得靠CFD来算。

核心方程有三个：

你可能会问：「这些方程我能手算吗？」当然不能。实际铸造过程是三维瞬态的，方程极其复杂。所以我们需要数值方法——把连续的空间和时间离散成网格，用计算机迭代求解。

我的经验：网格质量直接影响CFD结果。我在项目中遇到过，网格太粗，充型顺序都算反了。网格太细，算三天三夜不出结果。建议：浇口附近、薄壁区域加密，远离浇口处可以粗一些。

铸造CFD里有个特殊问题——自由表面。金属液和空气之间有个界面，这个界面怎么追踪？常用的有VOF法（流体体积法）。简单说，每个网格里标记一下：金属液占多少比例，空气占多少。比例从0到1，0.5就是界面位置。

我记得第一次用VOF法算充型，看到界面在屏幕上一点点推进，那种感觉——嗯，挺有成就感的。

铸造本质上是个散热过程。高温金属液浇进冷型腔，热量怎么散掉？这决定了凝固顺序、缩孔位置、甚至铸件性能。

传热有三种方式：

传热控制方程是傅里叶导热定律的推广形式。你想想看，一个三维物体，温度随时间变化，还要考虑内部热源（比如凝固潜热）——这就是传热方程要解的。

这里有个坑，我踩过：界面换热系数。铸件和铸型之间不是完美接触的，有气隙、有涂料层。这个换热系数怎么设？

避坑指南：我曾经直接用了软件默认的换热系数，结果算出来的冷却曲线和实测差了50℃。后来老老实实做了热物性测试，才把模型校准好。记住：默认参数仅供参考，关键数据一定要自己验证。

凝固是铸造的核心。金属液从液态变成固态，这个过程决定了铸件的微观组织和宏观缺陷。

凝固理论要掌握几个关键概念：

凝固仿真的核心是判断缩孔缩松。怎么判断？有个经典判据：Niyama判据。它计算的是温度梯度G和冷却速率R的比值G/√R。这个值越小，缩松倾向越大。

我刚开始用Niyama判据时，总觉得阈值设多少合适。后来发现，不同合金、不同工艺，阈值都不一样。铸铁和铸钢的阈值能差一个数量级。所以——

我的建议：别迷信某个固定阈值。先算一个案例，和实际剖切结果对比，反推出适合你工艺的阈值。这叫「标定」，是仿真工程师的基本功。

铸件冷却过程中，不同部位收缩不一致，就会产生应力。应力大了，要么热裂，要么冷裂，要么变形超差。

应力分析的理论基础是弹性力学和塑性力学。核心方程包括：

铸造应力分析有个特殊点：热-力耦合。温度场影响应力场（热膨胀），应力场反过来也会影响温度场（变形导致接触变化）。所以严格来说，这是个双向耦合问题。

不过实际工程中，我们常用「顺序耦合」——先算温度场，再把温度结果加载到应力模型里。这样计算量小很多，精度也够用。

注意：应力分析对网格质量要求很高。尤其是应力集中区域（比如圆角、截面突变处），网格必须足够密。我见过有人用粗网格算应力，结果最大应力位置都算错了——这会导致你误判哪里会开裂。

下面这张图，是我自己总结的铸造仿真理论框架。四个模块相互关联，缺一不可。

看到这张图，你应该明白了：铸造仿真不是孤立地算某个物理场，而是四个模块的耦合。CFD算充型，传热算温度，凝固算相变，应力算变形——它们环环相扣。

好了，第一章的内容就到这里。这些理论听起来可能有点枯燥，但它们是后续所有章节的基础。下一章我们会进入实际操作——怎么搭建一个靠谱的仿真模型。到时候你会发现，今天学的这些理论，每一个都会用到。

公众号:蓝海资料掘金营，微信deep3321