4. 导热与接触热阻建模
电子设备散热仿真中，固体域的热量传递路径往往比流体域更“隐蔽”：芯片到封装外壳、外壳到散热器、散热器到机箱，每一段都依赖材料导热和界面接触。若只给固体赋一个导热系数、把装配面默认成理想贴合，仿真常会明显高估散热能力。本节围绕四个关键点展开：固体导热系数的取值与方向性、界面材料（TIM）的等效处理、接触热阻的设置方法，以及薄壁结构用壳导热替代实体网格的等效原则。
一、固体导热系数：先分清各向同性还是各向异性
Fluent 的固体域材料需要给定导热系数。对铝、铜、普通塑料等各向同性材料，只需一个 Thermal Conductivity 值，例如纯铝约 202 W/(m·K)、纯铜约 387 W/(m·K)、ABS 塑料约 0.2 W/(m·K)。但电子设备中大量存在各向异性材料：PCB 的 FR-4 基材在平面内（x、y）导热约 0.3～0.5 W/(m·K)，而厚度方向（z）仅约 0.2～0.3 W/(m·K)；含铜箔层后，平面内等效导热可升到 5～20 W/(m·K)，厚度方向仍很低。石墨片、导热垫、热管等效层也常呈强各向异性。
在 Fluent 中，各向异性导热通过材料属性面板选择 orthotropic，再分别填入三个方向的导热系数。方向必须与建模坐标系一致：若 PCB 在建模时厚度方向是 z，则 z 方向填厚度导热，x、y 填平面导热。方向填错，热扩散路径会完全走偏。对于层叠结构，也可用“等效导热系数”把多层压成一个固体域，但等效值应按热阻串联/并联公式计算，而不是简单取平均。
二、界面材料与接触热阻：把“贴合”变成可计算的薄层
两个固体表面即使压紧，微观上也只有少量凸点接触，真实接触面积可能不到名义面积的 10%。热量在缝隙中只能通过空气或界面材料传导，形成接触热阻 Rc，单位常用 K·m²/W 或 K/W。Fluent 提供两种处理方式：一是用 Thermal Contact Resistance 在固体-固体交界面直接给定单位面积热阻；二是把 TIM 建成一层很薄的实体或壳，赋予其导热系数和厚度。
若用薄层实体，热阻为 R = t / k。例如导热硅脂 k = 3 W/(m·K)、厚度 0.05 mm，则 R = 0.05×10⁻³ / 3 ≈ 1.67×10⁻⁵ K·m²/W；导热垫 k = 1.5 W/(m·K)、厚度 0.5 mm，则 R ≈ 3.33×10⁻⁴ K·m²/W，比硅脂高约 20 倍。这说明同样面积下，厚导热垫会显著抬高界面温降。若已知热阻而不知厚度，可直接在交界面填 Thermal Contact Resistance，避免为 0.05 mm 薄层划分极小网格。
三、薄壁热阻等效：壳导热与实体网格的取舍
机箱壁、散热器鳍片、屏蔽罩、电池外壳等薄壁结构，厚度常在 0.2～2 mm，而整体尺寸可达数百毫米。若按实体划分网格，厚度方向需要至少 2～3 层单元，全局网格量会急剧上升。Fluent 的壳导热（Shell Conduction）允许用无厚度或给定厚度的壳单元参与导热计算，只需输入壁厚和材料导热系数，即可在面网格上求解厚度方向热阻。
壳导热的等效条件是：壳的面内导热由材料导热系数和厚度共同决定，厚度方向热阻为 t/k。因此壳的 Wall Thickness 必须填真实厚度，材料 k 填真实导热系数，二者不能互相“抵消”成任意组合。壳两侧若分别与流体和固体耦合，需确认壳的耦合面设置正确，否则会出现一侧绝热、另一侧温度异常的假象。
四、在 Fluent 中的设置要点与检查
实际操作时，建议按以下顺序推进。第一，在 Materials 中建立各固体材料，各向异性材料选择 orthotropic 并核对方向。第二，在 Cell Zone Conditions 中为每个固体域指定材料，确认没有误把流体域设为固体。第三，对固体-固体交界面，在 Thermal 选项卡中设置 Coupled，若需要接触热阻则勾选 Thermal Contact Resistance 并填入单位面积热阻；若 TIM 已建成薄层实体，则交界面保持 Coupled，由薄层材料承担热阻。第四，薄壁区域在 Wall 边界中启用 Shell Conduction，填写 Wall Thickness 和材料，并确认壳的两侧耦合面都正确。
判断设置是否合理，可以观察三个量：芯片结温与散热器基座温度的差值是否与 Q×R 估算一致；界面两侧温差是否集中在 TIM 层；壳导热区域的面内温度分布是否连续、无异常跳变。若界面温降远小于手算值，多半是接触热阻漏设或 TIM 厚度填错；若壳两侧温度完全对称且无梯度，则要检查壳是否真正参与了导热耦合。
最后给出一个便于复用的参数记录表，建议在每次建模时填写，避免材料、厚度、热阻三者混用。
| 位置 | 材料/界面 | k (W/m·K) | 厚度 (mm) | R (K·m²/W) | Fluent 设置 |
|---|---|---|---|---|---|
| 芯片-外壳 | 导热硅脂 | 3 | 0.05 | 1.67e-5 | 接触热阻或薄层实体 |
| 外壳-散热器 | 导热垫 | 1.5 | 0.5 | 3.33e-4 | 接触热阻或薄层实体 |
| PCB 平面内 | FR-4+铜 | 10（x,y） | — | — | orthotropic |
| PCB 厚度向 | FR-4 | 0.3（z） | — | — | orthotropic |
| 机箱壁 | 铝合金 | 202 | 1.0 | 4.95e-6 | Shell Conduction |
把上述热阻按串联关系累加，即可在仿真前快速估算总温降，作为结果合理性的第一道校验。固体导热与接触热阻建模的核心，不是把每个参数填得“看起来对”，而是让每一段热阻都有明确来源、明确单位和明确方向，这样仿真结果才具备工程判断价值。