第二章 仿真理论基础：传热、流动与凝固
各位工程师朋友，大家好。我是老张，干铸造仿真这行有十几年了。今天咱们聊聊仿真背后的那些理论。说实话，刚入行那会儿，我也觉得理论枯燥，但后来吃过亏才明白——不懂理论，仿真就是瞎蒙。
这一章，我打算从三个维度展开：传热学、流体力学、凝固理论。这三者就像三根柱子，撑起了铸造仿真的整个框架。咱们一个一个来。
2.1 传热学基础：温度场怎么算出来的？
传热学，说白了就是研究热量怎么跑。在铸造里，金属液从浇注到凝固，温度一直在变。这个温度场算不准，后面所有分析都是白搭。
传热有三种方式：热传导、热对流、热辐射。我一个个说。
2.1.1 热传导
热传导是热量在物体内部传递。你想想看，把一根铁棒一端烧红，另一端也会烫手——这就是热传导。在铸造仿真里，我们用傅里叶定律来描述：
q = -k · ∇T其中 q 是热流密度，k 是导热系数，∇T 是温度梯度。负号表示热量从高温往低温跑。
2.1.2 热对流
热对流发生在流体和固体表面之间。金属液在型腔里流动，热量被带走——这就是对流换热。牛顿冷却定律是基础：
q = h · (T_s - T_f)h 是对流换热系数，T_s 是固体表面温度，T_f 是流体温度。
2.1.3 热辐射
热辐射不需要介质，真空中也能传热。在铸造中，高温金属液表面会向外辐射热量。斯蒂芬-玻尔兹曼定律：
q = ε · σ · (T^4 - T_env^4)ε 是发射率，σ 是斯蒂芬-玻尔兹曼常数。注意，这里温度是四次方关系——温度越高，辐射越强。
2.2 流体力学基础：金属液怎么流动？
金属液在型腔里流动，不是随便乱跑的。它遵循流体力学的基本规律。我重点讲两个概念：层流与湍流，以及纳维-斯托克斯方程。
2.2.1 层流与湍流
你倒一杯水，水流是平稳的——这是层流。你开水龙头开到最大，水流变得混乱——这是湍流。判断标准是雷诺数 Re：
Re = ρ · v · L / μρ 是密度，v 是流速，L 是特征长度，μ 是动力粘度。
- Re < 2300：层流，流动平稳，热量传递慢
- 2300 < Re < 4000：过渡区，不稳定
- Re > 4000：湍流，流动混乱，传热快
2.2.2 纳维-斯托克斯方程
这个方程，搞流体的都绕不开。它描述了流体运动的基本规律。形式如下：
ρ · (∂v/∂t + v · ∇v) = -∇p + μ · ∇²v + ρ · g左边是惯性力，右边依次是压力项、粘性力项、重力项。说白了，就是牛顿第二定律在流体上的应用。
嗯，这里要注意——N-S 方程是非线性的，解析解很难求。仿真软件用的是数值方法，比如有限体积法。我刚开始做仿真时，总想手算验证，后来发现不现实。软件算出来的结果，只要网格够密、边界条件设对，基本可信。
2.3 凝固理论：金属怎么从液态变固态？
凝固过程，是铸造的核心。金属液从液态变成固态，不是一瞬间的事。它经历了形核、长大、枝晶生长，最后可能留下缩松缩孔。
2.3.1 枝晶生长
金属凝固时，晶体不是均匀长大的。它会长成树枝状——这就是枝晶。为什么会长成树枝？因为热量和溶质扩散不均匀。
枝晶生长有几个关键参数：
| 参数 | 含义 | 对枝晶的影响 |
|---|---|---|
| 温度梯度 G | 固液界面前沿的温度变化率 | G 越大，枝晶越细 |
| 生长速度 R | 固液界面推进速度 | R 越大，枝晶间距越小 |
| 过冷度 ΔT | 实际温度与平衡温度之差 | ΔT 越大，形核率越高 |
2.3.2 缩松缩孔形成机理
缩松和缩孔，是铸造的常见缺陷。它们是怎么来的？说白了，就是金属凝固时体积收缩，液体补缩不及时。
金属从液态到固态，体积会缩小。比如铝合金，体积收缩率约 6-7%。如果最后凝固的区域没有液体补充，就会形成空洞。
- 缩孔：集中在一个大区域，形状不规则，通常出现在热节处
- 缩松：分散在枝晶之间，呈海绵状，肉眼可能看不见
判断缩松缩孔，常用的是 Niyama 判据：
Ny = G / √(R)G 是温度梯度，R 是冷却速率。Ny 值越小，缩松风险越大。一般 Ny < 1 时，就要警惕了。
公众号:蓝海资料掘金营，微信deep3321