4. 热结构直接耦合与间接耦合方法
在ANSYS中做热结构耦合，核心分歧只有一个：温度场和位移场是同时求解，还是分两步先后求解。前者叫直接耦合（Direct Coupling），后者叫间接耦合（Sequential Coupling，也称载荷传递法）。选错方法，轻则算得慢，重则结果失真甚至不收敛。
一、两种方法的本质区别
直接耦合使用一个包含温度自由度和位移自由度的耦合单元（如SOLID226、SOLID227、PLANE223），在同一套方程中组装热传导矩阵与结构刚度矩阵，并引入热膨胀耦合项。其矩阵形式可简写为：
[K_tt 0 ] {T} {Q}
[K_ut K_uu] {u} = {F}
其中K_ut为热-结构耦合项，代表温度变化引起的等效热载荷。由于T和u在同一迭代步内更新，直接耦合天然处理强耦合问题：温度改变立即影响变形，变形（如接触传热、间隙变化）又反过来影响温度。
间接耦合则把问题拆开：先做纯热分析得到节点温度，再把温度作为体载荷读入结构分析。温度场一旦算出就固定不变，结构分析不再反馈给热分析。因此它只能处理单向耦合——热影响结构，结构不影响热。
二、适用条件：什么时候必须用直接耦合
判断依据是耦合是否双向。以下三类问题必须用直接耦合：
- 接触传热问题：两个零件接触面之间的热传导系数随接触压力、间隙变化。温度改变间隙，间隙又改变传热，形成闭环。
- 瞬态热冲击：如淬火、激光加热，温度梯度在毫秒级剧烈变化，热应变与温度场强烈耦合，分步求解会引入时间滞后误差。
- 大变形热问题：结构变形量级足以改变热边界条件（如辐射角系数、对流面积）。
反之，以下情况用间接耦合更划算：稳态热应力分析、电子封装热循环（焊点蠕变但温度场由外部决定）、管道热膨胀、以及结构变形远小于特征尺寸的常规工况。
三、操作流程对比
直接耦合流程（以SOLID226为例）：
ET,1,SOLID226,11 ! KEYOPT(1)=11 耦合场单元
MP,EX,1,2e11 ! 弹性模量
MP,ALPX,1,1.2e-5 ! 热膨胀系数
MP,KXX,1,45 ! 导热系数
MP,C,1,460 ! 比热
MP,DENS,1,7850 ! 密度
! 施加温度与位移边界
D,ALL,TEMP,25 ! 初始温度
BF,ALL,HGEN,1e6 ! 热生成
SOLVE关键点：耦合单元需要同时定义热物性和结构物性，边界条件也需同时施加温度约束和位移约束。
间接耦合流程：
! 第1步：热分析
ET,1,SOLID278
MP,KXX,1,45
SOLVE
FINISH
! 第2步：结构分析
ET,1,SOLID185
MP,EX,1,2e11
MP,ALPX,1,1.2e-5
LDREAD,TEMP,,,t,,,rth ! 读入温度场
SOLVE注意LDREAD命令中的时间点参数必须与热分析结果对应，瞬态问题需逐载荷步循环读取。
四、精度与效率对比
| 对比维度 | 直接耦合 | 间接耦合 |
|---|---|---|
| 耦合方向 | 双向 | 单向 |
| 求解规模 | 大（自由度翻倍） | 小（分步求解） |
| 计算时间 | 长，约为间接的2~5倍 | 短 |
| 收敛难度 | 高，需同时满足热与力收敛 | 低，各步独立收敛 |
| 网格要求 | 热与结构网格必须一致 | 可分别划分，热网格可粗 |
| 精度 | 强耦合问题精度高 | 弱耦合问题精度足够 |
| 结果文件 | 单一结果文件 | 需传递.rth文件 |
LDREAD读入温度，或时间点参数写错，导致结构分析中温度载荷为零，结果只剩机械载荷效应。务必检查结构分析输出中的温度分布是否与热分析一致。
五、选择决策
实际工程中可按以下顺序判断：先问“结构变形是否显著改变温度场”，若否则选间接耦合；若是，再问“是否涉及接触传热或瞬态热冲击”，若是则必须直接耦合。对于大多数常规热应力问题，间接耦合在精度和效率之间提供了最佳平衡。
需要特别说明的是，间接耦合并非“低精度”的代名词。在弱耦合条件下，它的结果与直接耦合几乎一致，而计算成本可能只有后者的三分之一。真正需要警惕的是误用——把强耦合问题交给间接耦合，得到的温度场和应力场都会偏离真实解。