4. 热-结构耦合分析：顺序耦合与直接耦合、热应变与热应力计算、耦合场单元选择

热-结构耦合分析，说白了就是算两件事：温度变化让零件胀了多少，胀不动的地方产生了多大的力。我刚开始做光机仿真时，总觉得把温度场算准了就行，结果结构变形算出来跟实测差一大截。后来才明白——耦合方式没选对，单元类型没用好，算出来的东西根本没法信。

这一章咱们就把耦合分析的门道捋清楚。我会结合自己踩过的坑，把顺序耦合、直接耦合、热应变计算和单元选择这几个关键点讲透。

4.1 顺序耦合 vs 直接耦合：什么时候该用哪个？

耦合分析有两种主流方法：顺序耦合和直接耦合。很多新手上来就问“哪个更准”，其实不存在绝对的优劣，关键看你的问题场景。

4.2 热应变与热应力计算：公式背后的物理直觉

热应变的核心公式很简单：

ε_thermal = α × (T - T_ref)

其中 α 是线膨胀系数，T 是当前温度，T_ref 是参考温度（通常是应力释放温度）。但实际工程中，坑全藏在细节里。

4.2.1 各向异性热膨胀：光机结构的噩梦

我遇到过最头疼的情况是碳纤维镜筒的热变形。碳纤维铺层方向不同，α 值可以差 10 倍以上。如果你用各向同性材料模型去算，结果完全不对。

正确的做法是定义正交各向异性热膨胀系数：

! ANSYS APDL 示例
MP, ALPX, 1, -0.5e-6   ! 纤维方向
MP, ALPY, 1, 12e-6     ! 垂直纤维方向
MP, ALPZ, 1, 30e-6     ! 厚度方向

4.2.2 热应力计算的本质

热应力不是直接由温度产生的，而是由“约束”产生的。你想想看：如果一根杆子自由膨胀，它不会有任何应力。只有当它被固定住，或者相邻材料膨胀量不一致时，才会产生应力。

所以热应力计算的本质是：

σ_thermal = E × (ε_total - ε_thermal)

其中 ε_total 是结构实际发生的应变（由位移场决定），ε_thermal 是自由热应变。两者的差值才是弹性应变，乘以弹性模量 E 就得到热应力。

4.3 耦合场单元选择：选错了，算出来就是废的

单元选择是热-结构耦合分析中最容易被忽视的环节。我见过太多人直接用结构单元加载温度体载荷，结果忽略了热-结构耦合单元的特殊要求。

4.3.1 常用耦合场单元一览

4.3.2 我的单元选择决策流程

我一般按这个逻辑来选：

1. 先判断耦合方式： 顺序耦合还是直接耦合？
2. 顺序耦合： 热分析用 SOLID90（20节点六面体热单元），结构分析用 SOLID186（20节点六面体结构单元）。网格可以共用，但单元类型不同。
3. 直接耦合： 首选 SOLID226（20节点六面体耦合单元）。如果涉及大变形接触，换 SOLID227。
4. 检查单元形状： 六面体优先，四面体慎用。耦合场四面体单元容易产生剪切闭锁，热应力结果偏硬。

4.4 本章知识体系：一张图看懂

下面这张图总结了热-结构耦合分析的核心逻辑，我建议你把它存下来，做项目时对照着看。

4.5 实战要点总结

最后，我把这一章的核心要点浓缩成几条，你写报告或做仿真时可以直接参考：

• 耦合方式： 光机结构分析优先选顺序耦合。除非你遇到摩擦生热或接触热阻随变形剧烈变化的情况，否则别碰直接耦合。
• 参考温度： 这是最容易出错的地方。T_ref 必须等于应力释放温度（如固化温度、装配温度），不是环境温度。
• 各向异性： 复合材料结构必须定义方向相关的热膨胀系数。用各向同性近似算碳纤维，结果会差一个数量级。
• 单元选择： 顺序耦合用 SOLID90 + SOLID186；直接耦合用 SOLID226。网格质量是命门，六面体 + 密实网格是唯一出路。
• 验证： 算完后用简单手算验证趋势。比如铝板均匀升温，自由膨胀量应该是 α × L × ΔT，如果仿真结果差太多，先查 T_ref 和材料参数。