第二章 铸造仿真基础理论：传热学基础、流体力学基础、凝固理论、应力场理论
各位同行，大家好。我是老张，干铸造仿真这行有十几年了。今天咱们聊聊铸造仿真最底层的那些理论。说实话，很多人觉得理论枯燥，但我要说——不懂这些，你做的仿真就是空中楼阁。
我个人习惯，每次接手一个新铸件，第一件事不是急着建模，而是先想想：这个件主要靠什么原理成型？热量怎么走？金属液怎么流？凝固顺序对不对？应力会不会裂？想通了这些，仿真才有方向。
核心观点：铸造仿真不是“点一下按钮就出结果”的黑盒子。传热、流动、凝固、应力，这四大理论是仿真的四根柱子。缺一根，你的结果就可能翻车。
2.1 传热学基础——热量怎么跑，铸件就怎么长
传热学，说白了就是研究热量怎么从高温区跑到低温区。铸造过程里，金属液从浇注温度降到室温，热量必须及时散走。散快了，铸件容易裂；散慢了，缩松缩孔找上门。
传热有三种方式：
- 热传导：热量在固体内部传递。比如砂型内部的热量传递。傅里叶定律是基础：q = -λ·grad(T)。λ是导热系数，不同材料差很多。钢的λ大概30-50 W/(m·K)，而砂型只有0.5-1.0。嗯，这里要注意，砂型的导热能力差，所以铸件冷却慢。
- 热对流：流体（金属液、空气）带着热量跑。牛顿冷却定律：q = h·(T_w - T_f)。h是对流换热系数，这个值很难精确，我一般靠经验公式估算。
- 热辐射：高温物体向外辐射能量。斯特藩-玻尔兹曼定律：q = ε·σ·T⁴。温度越高，辐射越强。浇注初期，金属液温度1400℃以上，辐射占比很大。
我在项目中遇到过一个问题：一个大型铸钢件，仿真出来的冷却曲线和实测差很多。后来发现，是我把辐射换热忽略了。高温阶段，辐射占主导，忽略它，冷却时间能差30%。
避坑指南：我曾经在模拟薄壁件时，把网格画得太粗，结果热传导计算失真。薄壁件厚度方向至少要有3层网格，否则温度梯度算不准。
2.2 流体力学基础——金属液怎么流，缺陷就怎么来
金属液在型腔里流动，不是简单的“倒进去就行”。它受重力、惯性力、粘性力共同作用。流体力学的基础方程是纳维-斯托克斯方程（N-S方程），但说实话，手算这个不现实，仿真软件里都是数值求解。
铸造流动有几个关键概念：
- 雷诺数（Re）：判断流动是层流还是湍流。Re = ρ·v·L/μ。铸造充型过程，Re通常很大，属于湍流。湍流模拟需要额外的湍流模型，比如k-ε模型。
- 自由表面：金属液和空气的界面。VOF方法（流体体积法）是主流。每个网格里，金属液占多少体积，空气占多少，算出来就知道自由表面在哪。
- 充型时间：浇注速度要合适。太快了，容易卷气；太慢了，冷隔、浇不足。我一般用伯努利方程估算浇口速度，再用仿真验证。
你想想看，为什么有些铸件容易产生气孔？说白了，就是金属液流动时把空气裹进去了。流动路径设计不好，空气排不出去，就留在铸件里了。
警告：湍流模型不是万能的。k-ε模型对强旋流、大曲率流动预测不准。遇到复杂流道，建议用k-ω SST模型或LES（大涡模拟）。
2.3 凝固理论——从液态到固态，决定铸件命运
凝固，是铸造最核心的环节。金属液从液态变成固态，体积会收缩（大约3-7%）。收缩得不到补缩，就会产生缩松缩孔。
凝固理论有几个关键点：
- 温度梯度：决定了凝固前沿的形态。正温度梯度（从型壁到中心温度升高）有利于柱状晶生长。负温度梯度会导致等轴晶。
- 固相率：当固相率达到0.7左右时，金属液基本失去流动性，补缩通道被切断。这个点叫“补缩临界固相率”。
- 凝固时间：用Chvorinov法则估算：t = C·(V/A)²。V是体积，A是表面积，C是常数。模数（V/A）越大，凝固越慢。
我记得有一次做铝合金缸盖的仿真，缩松预测总是不准。后来发现，是我用的凝固模型太简单了。铝合金凝固区间宽（液相线和固相线温差大），必须用“糊状区”模型，考虑枝晶间流动。
经验之谈：判断缩松位置，最简单的方法是看“最后凝固区域”。用仿真软件算一下凝固时间分布，最晚凝固的地方，大概率是缩松区。然后在这个位置放冒口，或者加冷铁。
2.4 应力场理论——铸件裂不裂，应力说了算
铸件冷却过程中，不同部位收缩不一致，就会产生内应力。应力超过材料强度，就开裂。应力场理论，就是算这个的。
应力来源主要有三种：
- 热应力：温度不均匀导致的热胀冷缩不一致。这是最主要的应力来源。
- 相变应力：固态相变（比如奥氏体→马氏体）伴随体积变化。铸铁的石墨化膨胀就是典型。
- 机械约束应力：砂型、砂芯对铸件收缩的阻碍。
应力计算的本构关系：σ = E·(ε_total - ε_thermal - ε_phase)。E是弹性模量，ε_total是总应变，ε_thermal是热应变，ε_phase是相变应变。
我建议，做应力仿真时，材料参数一定要准。尤其是高温下的弹性模量和屈服强度。很多材料库里的数据是室温的，高温下差很多。我曾经吃过这个亏，算出来的应力偏小，结果铸件实际裂了。
避坑指南：我曾经在模拟大型铸钢件时，忽略了砂型的退让性。砂型在高温下会软化，对铸件的约束没那么强。如果不考虑这个，算出来的应力会偏大，导致误判。
2.5 四大理论的耦合关系
这四大理论不是孤立的。传热影响凝固，凝固影响流动，流动又影响传热。应力场更是和温度场、相变场耦合在一起。
我画了一张图，帮你理清它们的关系：
从图上你能看到，传热和流体是双向耦合的。温度影响金属液粘度，粘度影响流动；流动又影响热量分布。凝固和传热也是双向的。凝固释放潜热，影响温度场；温度场又决定凝固速度。应力场则受温度场和凝固过程共同影响。
做仿真时，我建议先做流动+传热的充型分析，再做凝固分析，最后做应力分析。但要注意，有些软件支持全耦合，有些是顺序耦合。全耦合更准，但计算量大。顺序耦合效率高，但精度略差。
总结：传热学告诉你热量怎么走，流体力学告诉你金属液怎么流，凝固理论告诉你铸件怎么长，应力场理论告诉你铸件会不会裂。把这四个吃透了，你的仿真水平就上了一个大台阶。