3. 热-力耦合分析：温度场与应力场的耦合机制、顺序耦合与直接耦合方法、耦合分析软件操作流程
说到热-力耦合，我得先跟你聊聊一个真实案例。几年前我处理过一套大型铸钢件，材质是ZG35CrMo。浇注完冷却到室温，一探伤，裂纹！当时我第一反应是「应力太大了」。但仔细一查，问题出在哪儿？温度场和应力场没算准。说白了，你只算温度，或者只算应力，都抓不住问题的根。
这一节，我就把热-力耦合的核心逻辑掰开揉碎了讲。你想想看，铸件从液态到固态，温度从一千多度降到室温，这中间材料性能在变，约束在变，应力也在变。它们不是孤立的，是互相影响的。
3.1 耦合机制：温度场与应力场如何互相影响
先搞清楚一个基本问题：温度场怎么影响应力场？反过来，应力场又怎么影响温度场？
温度场对应力场的影响，这个好理解。温度变化引起热膨胀或收缩，如果收缩被约束住，就产生热应力。我习惯把这个过程分成三步：
- 热膨胀/收缩：温度升高，材料膨胀；温度降低，材料收缩。这是最直接的物理效应。
- 材料性能随温度变化：弹性模量、屈服强度、热导率、比热容……这些参数在高温下和室温下完全不一样。我记得有一次算一个高锰钢件，室温下屈服强度500多MPa，到了800℃直接掉到几十MPa。你不考虑这个，算出来的应力就是错的。
- 相变潜热与体积变化：凝固过程中，液态变固态，体积会收缩（有的合金还会膨胀）。这个体积变化会产生额外的应力。我曾经遇到一个球墨铸铁件，就是因为石墨化膨胀没算对，导致缩松和裂纹同时出现。
应力场对温度场的影响，这个相对弱一些，但也不能忽略：
- 塑性功生热：材料发生塑性变形时，一部分机械能会转化成热能。虽然这个量通常不大，但在高速变形或大变形情况下，不能忽略。
- 接触热阻变化：铸件和铸型之间，如果应力导致接触状态改变（比如脱开或压紧），热传导系数会变化。这个我在实际项目中吃过亏——没考虑接触热阻变化，算出来的冷却曲线和实测差了30%。
核心要点：热-力耦合的本质是「温度→热应变→应力→变形→接触状态→热传导→温度」这个闭环。你把这个闭环想清楚了，后面的方法选择就简单了。
3.2 顺序耦合与直接耦合：两种方法怎么选
搞清楚了耦合机制，接下来就是怎么算的问题。目前主流的方法有两种：顺序耦合和直接耦合。我分别说说它们的原理和适用场景。
3.2.1 顺序耦合（Sequential Coupling）
顺序耦合，说白了就是「先算温度，再算应力」。流程是这样的：
- 先做一次瞬态热分析，得到整个冷却过程中每个时间步的温度场。
- 把温度场结果作为载荷，映射到结构分析中。
- 在结构分析中，计算热应力和变形。
这个方法的好处是：计算量小，容易收敛。我个人的习惯是，对于大多数常规铸件，优先用顺序耦合。为什么呢？因为温度场和应力场的时间尺度差异很大——温度变化慢，应力响应快。顺序耦合正好能抓住这个特点。
我的经验：顺序耦合适合以下场景——
- 铸件结构不太复杂，没有剧烈的局部变形
- 材料性能随温度变化不是特别剧烈
- 你主要关心最终的残余应力分布，而不是过程中的瞬态应力
3.2.2 直接耦合（Direct Coupling）
直接耦合就厉害了。它是在同一个求解器中，同时求解温度场和应力场。每个时间步，温度场和应力场互相迭代，直到收敛。
这个方法精度高，但代价也大——计算量是顺序耦合的3到5倍，而且收敛性差。我记得有一次算一个大型汽轮机缸体，用直接耦合，算了整整一周才出结果。中间还崩了两次，气得我差点砸电脑。
注意：直接耦合不是万能的。以下情况才值得用——
- 铸件有剧烈的局部塑性变形（比如热裂区域）
- 接触状态变化很大（比如铸件和铸型之间反复脱开和压紧）
- 你需要精确捕捉某个时间点的瞬态应力（比如开箱时刻）
3.3 耦合分析软件操作流程
理论讲完了，咱们来点实际的。我以ProCAST为例，说说热-力耦合分析的具体操作流程。其他软件（如Abaqus、MAGMA）逻辑类似，只是菜单名称不同。
3.3.1 前处理阶段
- 网格划分：热分析和力分析对网格的要求不同。热分析可以粗一些，力分析需要细一些。我建议在应力集中区域（比如圆角、壁厚变化处）加密网格。
- 材料参数输入：这是最关键的步骤。你需要输入随温度变化的：弹性模量、泊松比、热膨胀系数、屈服强度、硬化模量、热导率、比热容、密度、潜热……
- 边界条件设置：热边界（换热系数、环境温度）和力学边界（约束、接触定义）都要设。这里有个坑——力学约束不能过约束，否则算出来的应力会偏大。
避坑指南：我曾经犯过一个错误——把铸件底部完全固定，结果算出来的应力比实测大了两倍。后来改成「只约束法向位移」，结果就对了。记住：铸件在冷却过程中是可以自由收缩的，你约束得越多，应力越大。
3.3.2 求解设置
如果你用顺序耦合，操作流程是这样的：
// ProCAST 顺序耦合设置示例
// 第一步：热分析
1. 选择分析类型：Thermal
2. 设置时间步长：初始0.01s，最大10s
3. 设置总时间：到室温为止
4. 求解并保存温度场结果
// 第二步：应力分析
1. 选择分析类型：Stress
2. 导入温度场结果文件
3. 设置力学边界条件
4. 求解如果你用直接耦合，操作更简单，但计算更慢：
// ProCAST 直接耦合设置示例
1. 选择分析类型：Thermal + Stress (Coupled)
2. 设置耦合迭代次数：每个时间步3-5次
3. 设置收敛容差：温度0.1℃，应力1MPa
4. 求解3.3.3 后处理分析
算完之后，重点看这几个指标：
| 指标 | 含义 | 判断标准 |
|---|---|---|
| 等效米塞斯应力 | 综合应力水平 | 超过屈服强度则可能产生塑性变形 |
| 最大主应力 | 拉应力大小 | 超过抗拉强度则可能开裂 |
| 残余变形量 | 铸件最终变形 | 根据公差要求判断 |
| 热裂判据 | 基于应变或应力的热裂预测 | 超过临界值则可能热裂 |
我的习惯：后处理时，我会把温度场和应力场放在同一个时间轴上对比。比如，看某个时刻的温度梯度大不大，对应的应力集中区域在哪里。这样能快速定位问题。
3.4 知识体系总览
为了让你更直观地理解这一章的内容，我画了一张流程图。它把热-力耦合的核心逻辑、方法选择、操作流程串在了一起。
这张图把这一章的核心逻辑串起来了。你从耦合机制出发，根据实际需求选择方法，然后按照操作流程一步步执行。记住：没有最好的方法，只有最合适的方法。
好了，这一节就到这里。热-力耦合这块内容，说白了就是「温度算准了，应力才能算准」。你回去可以拿一个简单的铸件试试顺序耦合，先跑通流程，再慢慢优化。有什么问题，咱们下一节接着聊。