第二章：仿真技术基础——有限元、流体与热力耦合

各位工程师朋友，大家好。今天咱们聊聊仿真技术的三大基石：有限元分析（FEA）、计算流体动力学（CFD）和热-力耦合分析。说实话，这三块内容单独拎出来都能讲三天三夜，但咱们这门课的核心是多材质铸造，所以我只挑最实用的部分讲。

我个人习惯把仿真技术比作「三把刀」——FEA是切应力场的刀，CFD是切流场的刀，热-力耦合是把两把刀绑在一起用。你想想看，铸造过程里金属液流动、凝固收缩、热应力开裂，哪个环节离得开这三样？

2.1 有限元分析（FEA）基础——别把网格当儿戏

有限元分析，说白了就是把一个连续体切成小块块，然后近似求解。我在项目中遇到过不少工程师，一上来就加密网格，觉得网格越细越准。其实这是个误区。

FEA的核心思想就三个字：离散化。 把无限自由度的问题变成有限自由度。你想想看，一个铸件形状再复杂，切成几万个小四面体，每个小四面体上的位移用形函数插值，最后组装成一个大方程组——这就是FEA的底牌。

FEA的数学基础其实不复杂：

// 线性静力分析的核心方程
[K]{u} = {F}

其中：
[K] —— 整体刚度矩阵（由单元刚度矩阵组装）
{u} —— 节点位移向量（待求解）
{F} —— 节点载荷向量（外力+边界条件）

嗯，这里要注意：刚度矩阵[K]是奇异的，必须施加足够的位移约束才能求解。我见过新手忘了加约束，结果求解器报错「矩阵奇异」，急得满头大汗。

对于铸造仿真，FEA主要用在：

• 应力分析： 凝固收缩产生的热应力、残余应力
• 变形预测： 铸件冷却后的翘曲、尺寸偏差
• 模态分析： 砂型/金属型的振动特性（高速造型机用得上）

2.2 计算流体动力学（CFD）基础——流动的奥秘

CFD，说白了就是解纳维-斯托克斯方程。但你别被这个高大上的名字吓到，在铸造仿真里，我们关心的是：金属液能不能顺利充满型腔？会不会卷气？冷隔出现在哪里？

我记得刚入行时，师傅跟我说：「CFD算流动，三分模型七分边界。」当时不理解，后来自己算了一个复杂薄壁件的充型过程，才发现边界条件设错了，结果全白算。

CFD的核心控制方程：

连续性方程：∂ρ/∂t + ∇·(ρv) = 0
动量方程：∂(ρv)/∂t + ∇·(ρvv) = -∇p + ∇·τ + ρg
能量方程：∂(ρh)/∂t + ∇·(ρvh) = ∇·(k∇T) + S

在铸造充型仿真中，我们通常做以下简化：

• 假设金属液为不可压缩牛顿流体（大部分铝合金、铸铁适用）
• 采用VOF方法追踪自由液面（就是金属液和空气的界面）
• 考虑湍流效应（高流速浇注系统必须开湍流模型）

CFD在铸造中的典型应用：

应用场景	关键参数	常见问题
充型过程模拟	充型时间、速度场、压力场	卷气、冷隔、浇不足
凝固过程模拟	温度梯度、固相分数、缩松缩孔	热节、偏析、裂纹
浇注系统设计	流量分配、雷诺数、压头损失	冲砂、紊流、氧化夹杂

2.3 热-力耦合分析基础——温度与应力的纠缠

热-力耦合，说白了就是温度场和应力场互相影响。铸造过程中，金属液从高温冷却到室温，温度变化导致体积收缩，收缩受到约束就产生应力——这就是热-力耦合的底层逻辑。

我做过一个双金属复合铸造的案例，外层是高铬铸铁，内层是球墨铸铁。两种材料的热膨胀系数不一样，冷却过程中界面处产生了巨大的热应力。如果不做热-力耦合分析，根本算不准界面结合强度。

热-力耦合的两种方式：

• 顺序耦合： 先算温度场，再把温度结果作为载荷施加到结构分析中。简单、稳定，大部分铸造仿真用这个。
• 完全耦合： 同时求解温度和位移。理论上更精确，但计算量巨大，一般只在特殊场合用（比如热冲击、快速加热）。

热-力耦合分析在铸造中的难点：

1. 材料非线性： 高温下材料的力学性能（弹性模量、屈服强度）随温度剧烈变化
2. 接触非线性： 铸件与铸型之间的接触状态（接触/分离、摩擦）不断变化
3. 相变潜热： 凝固过程中释放的潜热会影响温度场，进而影响应力场

最后说一句，这三块技术不是孤立的。FEA算应力需要温度场，CFD算流动需要考虑凝固潜热，热-力耦合又需要FEA和CFD的结果——它们是一个整体。多材质铸造仿真，玩的就是这三者的协同。

好了，这一章的基础就讲到这里。记住：理论是死的，工程是活的。多动手算几个案例，比看十遍书都管用。