4、有限元热分析入门：网格划分原则、边界条件设定、稳态与瞬态分析
各位好，我是老张。今天咱们聊聊有限元热分析的基础操作。说实话，我刚入行那会儿，觉得只要把模型建好、网格一画、边界一加，软件跑完就完事了。结果呢？第一次做某空间相机的主镜热分析，算出来的温度场跟实测差了快8度。后来才发现，网格太粗、边界条件设错了。嗯，从那以后，我再也不敢小看这些“入门”步骤了。
4.1 网格划分：不是越密越好
网格划分是热分析的基石。你想想看，网格就像你测量温度时放的“温度计”，放得越多，测得越准，但计算量也越大。我个人的习惯是：关键区域加密，非关键区域放松。
核心原则：网格密度由温度梯度决定。梯度大的地方（如热源附近、接触界面），网格要密；温度均匀的区域（如远离热源的壳体），网格可以粗。
具体来说，我一般遵循这几条：
- 热源附近：网格尺寸取最小特征尺寸的 1/5 到 1/10。比如一个 10mm 的加热片，网格边长控制在 1-2mm。
- 接触界面：至少 2-3 层网格过渡。我曾经吃过亏，两个零件接触面只画了一层网格，结果热阻算出来完全不对。
- 薄壁结构：厚度方向至少 2 层网格。对于 1mm 厚的镜筒壁，我通常画 3 层。
- 整体控制：网格数量控制在 10 万到 50 万之间。太少了精度不够，太多了算不动。
| 结构类型 | 推荐网格尺寸 | 厚度方向层数 | 典型应用 |
|---|---|---|---|
| 主镜/透镜 | 2-5 mm | 3-5 层 | 光学元件热变形 |
| 镜筒/支撑结构 | 3-8 mm | 2-3 层 | 热传导路径 |
| 热源（加热片/芯片） | 0.5-2 mm | 2 层 | 局部热流加载 |
| 外壳/散热器 | 5-15 mm | 1-2 层 | 整体温度分布 |
小技巧：我习惯先用粗网格算一遍，看看温度梯度分布。然后在梯度大的区域局部加密，再算一遍。这样既省时间，又能保证关键位置精度。
4.2 边界条件设定：热分析的灵魂
边界条件设错了，网格再密也没用。说白了，边界条件就是告诉软件“热量从哪来、往哪去”。
常见的边界条件有四种：
- 温度边界：直接指定某个面的温度值。比如，安装基座固定在 20°C 的恒温台上。
- 热流边界：指定单位面积的热流量。比如，太阳照射到镜面的热流密度。
- 对流边界：指定对流换热系数和环境温度。比如，镜筒外表面与空气的自然对流。
- 辐射边界：指定辐射率和环境温度。比如，镜面与深空的热辐射交换。
注意：我曾经犯过一个低级错误——把辐射边界和对流边界同时加在同一个面上，结果热量被重复计算，温度高了 15°C。记住：一个面只能有一种热边界条件，除非你明确知道两者同时存在且互不干扰。
我个人做光机热分析时，边界条件设定遵循这个流程：
- 先确定热源：是内部发热（如电机、芯片）还是外部热源（如太阳、地球反照）？
- 再确定散热路径：热量通过传导、对流还是辐射散走？
- 最后确定环境条件：温度、压力、风速、背景辐射温度等。
4.3 稳态分析：看最终状态
稳态分析，说白了就是看系统“稳定下来”之后的温度分布。比如，卫星在轨道上运行一段时间后，各部件温度不再变化，这时候就是稳态。
稳态分析的应用场景：
- 评估光学系统在恒定热载荷下的热变形
- 确定热控系统的加热功率是否足够
- 验证热设计是否满足指标要求
我一般这样设置稳态分析：
1. 定义材料属性（导热系数、密度、比热容）
2. 划分网格
3. 施加边界条件（温度/热流/对流/辐射）
4. 设置求解器（直接法或迭代法）
5. 求解并查看结果经验之谈：稳态分析收敛很快，一般几分钟到半小时就能出结果。但要注意，如果模型中有辐射边界，收敛会慢一些。我习惯先不加辐射算一遍，再加辐射算第二遍，对比看看辐射的影响有多大。
4.4 瞬态分析：看变化过程
瞬态分析关注的是“温度随时间怎么变化”。比如，卫星从阴影区进入光照区，镜面温度在几分钟内从 -50°C 升到 +30°C，这个过程就需要瞬态分析。
瞬态分析的关键参数：
| 参数 | 说明 | 我的建议 |
|---|---|---|
| 时间步长 | 每次计算的时间间隔 | 取最小热时间常数的 1/10 到 1/5 |
| 总时间 | 分析的总时长 | 至少覆盖 3-5 个热时间常数 |
| 初始温度 | 起始时刻的温度分布 | 用稳态分析结果作为初始条件 |
| 输出步数 | 保存结果的时间点数量 | 至少 20-50 步，捕捉温度变化细节 |
避坑指南：我曾经做过一个瞬态分析，时间步长设得太大（0.1秒），结果温度曲线像锯齿一样，完全没法用。后来把步长缩小到 0.01 秒，曲线才平滑。记住：时间步长越小，精度越高，但计算时间也越长。一般先粗算找趋势，再细算定数值。
4.5 稳态 vs 瞬态：怎么选？
很多新手问我：“老张，我该用稳态还是瞬态？”我的回答是：看你的问题性质。
| 对比项 | 稳态分析 | 瞬态分析 |
|---|---|---|
| 适用场景 | 恒定热载荷、长时间运行 | 变化热载荷、短时间过程 |
| 计算时间 | 短（分钟级） | 长（小时级甚至天级） |
| 结果类型 | 最终温度分布 | 温度随时间变化曲线 |
| 是否需要初始条件 | 不需要 | 需要 |
| 典型应用 | 热设计验证、热平衡测试 | 热冲击分析、温度循环测试 |
我个人习惯是：先做稳态分析，摸清整体温度分布；如果发现温度变化剧烈或存在周期性热载荷，再补做瞬态分析。这样既高效又全面。
4.6 知识体系总览
下面这张图是我自己总结的有限元热分析知识体系，帮你理清思路：
好了，这一章的内容就这些。网格划分、边界条件、稳态与瞬态分析，这三块是有限元热分析的根基。你把这些搞明白了，后面做光机热耦合分析就顺了。记住：网格是骨架，边界是灵魂，分析类型是方向。三者缺一不可。