热有限元分析

我做过一个平板稳态导热分析项目，客户想验证一块受热平板内的温度分布是否符合傅里叶导热定律。这个案例非常适合作为热仿真的入门教学。

左侧图是平板温度场云图。平板左侧是热壁，温度 T₁=500 K，颜色偏红；右侧是冷壁，温度 T₂=300 K，颜色偏蓝；中间从左到右形成均匀的温度梯度。这是一个典型的一维稳态导热问题，没有内热源，材料导热系数 k=20 W/(m·K) 为常数。

右侧图是沿平板厚度方向的温度分布 T(x)。横轴是位置 x，从 0 到 8 m；纵轴是温度，从 300 K 到 500 K。红色实线是一条严格的直线，说明温度随距离线性下降，斜率为 −(T₁−T₂)/L = −25 K/m。图中标注了热流密度 q=12.5 W/cm²，这个值可以用傅里叶定律直接算出来：q = −k·dT/dx = 20 × 25 = 500 W/m² = 12.5 W/cm²（注意单位换算）。

这个案例我用 ANSYS Mechanical 做的。建模时取长 8 m、宽 1 m、厚 0.1 m 的平板，左右两端分别设 500 K 和 300 K 温度边界，其余面绝热。材料属性里只要输入导热系数 k=20 W/(m·K)。网格可以粗一些，一维稳态导热用少量单元就能得到精确解。求解后输出温度云图和热通量矢量。

踩坑经验：第一个是边界条件。如果侧面不是绝热而是对流，温度分布就不是直线，需要加对流换热系数。第二个是单位换算：W/m² 和 W/cm² 差 10000 倍，客户一开始把 q 写成 12.5 W/m²，差了四个数量级，我发现后纠正。第三个是材料各向异性：如果 k 随方向变化，温度场会出现弯曲，但这个案例假设各向同性。

工程经验：稳态导热看似简单，但它是判断网格、边界条件和材料参数是否正确的第一道关卡。理论解 q = k·ΔT/L 必须和仿真结果一致，否则后续复杂热分析都不可信。我建议每个热仿真项目都先用一个简单解析解做验证，再推广到复杂几何。

更多 ANSYS 稳态/瞬态热分析和传热仿真的实战案例，我汇总在 [热有限元分析](https://www.keyanxueshu.com/category/fea/) 栏目里。[热有限元分析](https://www.keyanxueshu.com/)

这个平板案例虽然简单，但可以延伸出很多工程判断。比如多层材料串联导热时，总热阻等于各层热阻之和 R_total = Σ(L_i/k_iA)，热流密度 q 处处相等，温差按热阻比例分配。如果某一层 k 很小，它两端的温差就很大。工程中常把隔热层放在高温侧还是低温侧有争议，实际上对稳态总热阻没有影响，但会影响温度梯度分布和材料耐温要求。并联导热时则按热导加权分配热流。我建议客户把这套“热阻网络”思维和有限元结果对照，复杂结构的等效导热系数可以用 FE 结果反推，这对封装、散热设计很有用。