3. 热结构耦合有限元分析

在热结构协同拓扑优化的流程中，有限元分析承担着“把物理场翻译成设计灵敏度”的角色。热载荷不像机械载荷那样直接以力或压强给出，它先以温度场的形式存在，再通过材料的热膨胀系数转化为结构内部的等效热应变与热应力。因此，本节的核心问题是：温度场与位移场之间如何耦合求解、热载荷如何映射为结构载荷、热单元与结构单元如何匹配，以及如何用典型算例验证建模的正确性。

3.1 弱耦合与强耦合的求解策略

热结构耦合按双向作用强度分为弱耦合（顺序耦合）与强耦合（直接耦合）两类。判断依据是：温度场是否显著受结构变形影响。

弱耦合求解分两步：先独立求解热传导方程得到节点温度 T，再把 T 作为体载荷代入结构平衡方程求解位移 u。其控制方程为：

热传导：  K_t · T = Q
结构平衡：K_s · u = F_mech + F_th(T)
其中热载荷向量：F_th = ∫ Bᵀ · D · α · (T - T_ref) dV

弱耦合适用于热变形不改变传热路径的情形，例如散热板、热压电梁的小变形分析。它的优点是两个物理场可分别用成熟求解器处理，计算成本低、便于验证。

强耦合则把温度与位移放在同一方程组中迭代求解，耦合项来自热弹性本构与几何非线性。当结构变形引起接触状态变化、间隙改变或对流边界移动时，必须使用强耦合。工程上常用“交替迭代”实现：每步先算温度场，更新几何与边界，再算位移场，检查温度与位移的残差是否同时收敛。若残差在两次迭代间振荡，应减小松弛因子或改用牛顿-拉夫逊整体求解。

3.2 热载荷向结构载荷的映射

热载荷映射的本质是把节点温度转换为等效节点力。对每个单元，先由形函数插值得到单元内温度分布，再计算热应变 ε_th = α(T - T_ref)，最后通过虚功原理得到等效节点力：

F_e = ∫_V Bᵀ · D · α · (T - T_ref) dV

其中 B 为应变-位移矩阵，D 为弹性矩阵，α 为热膨胀系数，T_ref 为无应力参考温度。映射过程中有三个易错点：一是 T_ref 必须与材料数据来源一致，否则会引入虚假热应力；二是温度场与位移场的网格若不一致，需要做插值映射，插值误差会直接进入热载荷；三是各向异性材料中 α 是向量或张量，方向必须与材料坐标系对齐。

3.3 热单元与结构单元的匹配

热单元与结构单元的自由度不同：热单元每个节点只有温度自由度，结构单元每个节点有平动（必要时还有转动）自由度。匹配的关键是“网格一致、节点对应、材料一致”。

• 网格一致：同一区域的热网格与结构网格应尽量共用节点，避免插值误差。若必须不同，结构网格应比热网格更细，或采用同一套网格分别赋予热与结构属性。
• 单元类型对应：二维热单元对应二维结构单元，三维实体热单元对应三维实体结构单元；壳结构若需温度沿厚度变化，应使用分层壳或实体单元。
• 材料参数一致：热分析需要导热系数 k、密度 ρ、比热 c；结构分析需要弹性模量 E、泊松比 ν、热膨胀系数 α。同一材料在两套属性中必须使用同一组温度相关数据。

3.4 典型算例：热压电梁与散热板

热压电梁算例。取一层压电材料与一层金属粘结而成的悬臂梁，长度 100 mm、宽度 10 mm、总厚度 2 mm，压电层厚 0.5 mm。左端固定，初始温度 20 ℃，压电层施加 80 ℃ 温度载荷，金属层保持 20 ℃。建模时压电层用热-压电耦合单元，金属层用热-结构耦合单元，界面共用节点。弱耦合求解后，自由端挠度约 0.35 mm（示例量级，随材料参数变化）。验证方法是与解析双层梁热弯曲公式对比，若偏差超过 10%，应检查界面节点是否共用、α 方向是否正确。

散热板算例。取 200 mm × 200 mm × 3 mm 铝板，中心 40 mm × 40 mm 区域施加 20 W 热源，四边对流换热系数 15 W/(m²·K)，环境温度 25 ℃。先做稳态热分析得到温度场，最高温约 68 ℃（示例值），再映射为热载荷做结构分析，板中心面外位移约 0.12 mm。验证时检查热平衡：总发热量应等于对流散热量，误差应小于 2%。

3.5 验证方法与常见问题

验证热结构耦合模型，建议按三个层次进行。第一层是热平衡验证：稳态下总输入热量应等于总散热量，误差控制在 2% 以内。第二层是自由热膨胀验证：对无约束均质体施加均匀温升，理论上不产生热应力，若结果出现明显应力，说明 T_ref 或约束设置有问题。第三层是与解析解或文献算例对比：双层梁热弯曲、散热板中心挠度都有成熟解析公式，偏差应小于 10%。

完成上述建模与验证后，温度场与位移场之间的映射关系就建立起来了。后续拓扑优化中，设计变量改变材料分布，会同时影响导热路径与刚度路径，这正是热结构协同优化需要处理的核心矛盾。