第2章模拟基础理论

第二章：模拟基础理论——充型凝固的物理本质

各位工程师朋友，大家好。我是老张，干铸造CAE这行快二十年了。今天咱们聊聊模拟的基础理论。说实话，很多人觉得理论枯燥，但我要说——不懂这些，你算出来的结果就是一堆彩色云图，根本没法指导生产。

这一章，我重点讲三个核心：充型过程的流体力学、凝固过程的传热学，还有那个让人头疼的相变潜热。嗯，咱们一个一个来。

充型，说白了就是金属液怎么流进型腔。你想想看，这跟水管里流水有啥区别？区别大了！金属液是高温、高粘度、带自由表面的流体，而且流动过程中还在散热、凝固。

核心控制方程：Navier-Stokes方程组

我个人习惯把N-S方程拆成三部分来理解：

我在项目中遇到过一个问题：一个铝合金壳体件，模拟充型总是卷气。后来发现是浇注系统设计不合理，金属液在横浇道里形成了紊流。紊流这东西，说白了就是流体乱窜，裹挟着空气进型腔。怎么判断？看雷诺数Re：

Re = ρ * v * D / μ

其中：
ρ — 金属液密度 (kg/m³)
v — 流速 (m/s)
D — 特征尺寸 (m)
μ — 动力粘度 (Pa·s)

当Re > 2300，紊流就来了。我一般控制在层流范围，实在不行就用紊流模型，比如k-ε模型。但记住，紊流模型只是近似，别太迷信。

关键点：充型模拟中，自由表面的追踪是难点。VOF方法（Volume of Fluid）是目前的主流，它把每个网格单元里金属液的体积分数算出来，0是空气，1是金属，中间值就是自由表面。

凝固，就是热量往外跑的过程。传热有三种方式：传导、对流、辐射。在精密铸造里，传导是主角，对流和辐射也不能忽视。

傅里叶导热定律

这个公式简单，但很实用：

q = -k * (dT/dx)

q — 热流密度 (W/m²)
k — 导热系数 (W/(m·K))
dT/dx — 温度梯度 (K/m)

负号表示热量从高温向低温传递。嗯，这里要注意：金属液的导热系数随温度变化很大，尤其是固液两相区。我习惯用温度相关的物性参数，别偷懒用常数。

凝固过程中的热传导方程

三维非稳态导热方程长这样：

ρ * Cp * (∂T/∂t) = k * (∂²T/∂x² + ∂²T/∂y² + ∂²T/∂z²) + Q

Q — 内热源项，这里就是相变潜热释放的热量

这个方程，说白了就是：金属液降温放出的热量 = 传导出去的热量 + 相变潜热。我曾经遇到一个案例，模拟出来的凝固时间比实际短了30%，后来发现是忽略了铸件和铸型之间的界面热阻。这个界面热阻，说白了就是接触不紧密，有气隙，热量传得慢。

我的经验：界面换热系数不是常数。凝固初期，金属液和铸型紧密接触，换热系数大；凝固后期，铸件收缩产生气隙，换热系数急剧下降。我一般用温度相关的换热系数曲线，或者用热接触模型。

相变潜热，这是凝固模拟里最让人头疼的。金属从液态变成固态，要释放大量的热量。这个热量如果不处理好，温度场算出来就是错的。

等效比热法

我比较喜欢用这个方法。思路很简单：把潜热折算到比热容里，让比热容在凝固温度区间内变大。

Cp_eff = Cp + L / (T_liquidus - T_solidus)

Cp_eff — 等效比热容
L — 潜热 (J/kg)
T_liquidus — 液相线温度
T_solidus — 固相线温度

举个例子：铝合金A356，潜热约389 kJ/kg，固液相线温差约60°C。算下来，等效比热容比实际比热容大了将近6倍。你想想看，这个热量有多大！

警告：等效比热法有个坑——如果固液相线温差很小（比如共晶合金），等效比热容会变得非常大，导致计算不稳定。这时候我建议用热焓法。

热焓法

热焓法更直接，把潜热作为内热源处理。每个时间步，根据温度计算固相分数，然后更新热焓：

H = Cp * T + L * (1 - f_s)

H — 热焓 (J/kg)
f_s — 固相分数 (0~1)

然后能量方程变成：

ρ * (∂H/∂t) = k * ∇²T

这个方法的好处是，不需要显式处理潜热，数值稳定性好。我一般在做大型铸件模拟时用热焓法，小件用等效比热法，省时间。

说了这么多，咱们用一张图把整个知识体系串起来。这张图是我自己画的，把充型、凝固、潜热处理的关系理清楚了。

这张图把三大模块的关系说清楚了。充型流体力学提供初始温度场和流场，凝固传热学计算温度变化，相变潜热处理保证能量守恒。三者耦合，才能得到准确的模拟结果。

最后，我分享几个实战中踩过的坑：

总结一下：充型看流场，凝固看温度，潜热是关键。三者搞明白了，你的模拟结果就能指导生产了。嗯，今天就聊到这儿，下次咱们讲网格划分的实战技巧。