接触熱抵抗
理论与物理
什么是接触热阻
为什么两个固体压在一起，界面处温度会发生跳跃？
从微观层面看，由于表面粗糙度，接触面实际上只是点接触（真实接触点）。实际接触面积仅为表观面积的1%左右，其余部分是空气间隙。这种不完全接触导致了温度跳跃。
基本公式
接触热阻的定义式是这样的。
接触热导是 $h_c = 1/R_c$ [W/(m2K)]。
典型值大概是多少？
根据条件不同，数值会相差几个数量级。
| 接触条件 | $h_c$ [W/(m2K)] |
|---|---|
| 铝与铝・抛光面・高压力 | 10000〜25000 |
| 钢与钢・机加工面・中压力 | 2000〜5000 |
| 使用导热硅脂 | 5000〜50000 |
| 空气间隙（0.1mm） | 250 |
| 真空中・低压力 | 100〜500 |
Cooper-Mikic-Yovanovich模型
理论模型的代表是Cooper-Mikic-Yovanovich (CMY) 相关式。
其中 $k_s = 2k_1k_2/(k_1+k_2)$ 是调和平均热导率，$m$ 是表面斜率，$\sigma$ 是复合粗糙度，$P$ 是接触压力，$H_c$ 是显微硬度。
提高压力会降低接触热阻呢。
没错。在螺栓连接部位的设计中，紧固扭矩决定了热路径的性能。扭矩不足会成为热设计的致命风险。
接触热阻的发现源自NASA
接触热阻（TCR）在工程上受到关注始于1950年代的太空开发。由于真空中没有对流，螺栓连接部位的TCR成为主导性的热阻，NASA的Glenn研究中心在1959年构建了第一个系统性的实测数据库。至今其数据仍作为设计初期的参考值使用。
各项的物理意义
- 蓄热项 $\rho c_p \partial T/\partial t$：单位体积的热能储存率。【日常示例】铁锅不易热也不易冷，而铝锅易热易冷——这是密度 $\rho$ 与比热 $c_p$ 的乘积（热容）的差异。热容大的物体温度变化缓慢。水的比热非常大（4,186 J/(kg·K)），因此沿海地区气温比内陆稳定。在非稳态分析中，此项决定了温度随时间的变化速率。
- 热传导项 $\nabla \cdot (k \nabla T)$：基于傅里叶定律的热传导。与温度梯度成比例的热流。【日常示例】将金属勺子放入热锅，手柄也会变热——因为金属的热导率 $k$ 高，热量能快速从高温侧传到低温侧。木勺不烫是因为 $k$ 小。隔热材料（如玻璃棉）的 $k$ 极小，即使有温度梯度也难以传热。这是将“有温差的地方就有热流”这一自然趋势公式化的结果。
- 对流项 $\rho c_p \mathbf{u} \cdot \nabla T$：伴随流体运动的热输送。【日常示例】吹风扇感到凉爽，是因为风（流体流动）带走了体表附近的暖空气，并供应了新鲜的冷空气——这就是强制对流。暖气使房间天花板附近变暖，是因为受热空气因浮力上升的自然对流。PC的CPU散热器风扇也是通过强制对流散热。对流是比热传导高效几个数量级的热输送手段。
- 热源项 $Q$：内部发热（焦耳热、化学反应热、辐射吸收等）。单位: W/m³。【日常示例】微波炉通过食品内部的微波吸收（体积发热）加热。电热毯的加热线通过焦耳发热（$Q = I^2 R / V$）变暖。锂离子电池充放电时的发热、刹车片的摩擦热也需要作为热源在分析中考虑。与从外部向“表面”施加热量的边界条件不同，热源项表示“内部”的能量生成。
假设条件与适用范围
数值解法与实现
FEM中的接触热阻建模
在FEM中如何表现接触热阻？
在Ansys Mechanical中通过接触单元的TCC (Thermal Contact Conductance) 设置。在Abaqus中可通过*GAP CONDUCTANCE 关键字定义压力相关表格。
如果要考虑压力依赖性，就需要与结构分析耦合吧。
是的。在Ansys Workbench中，可以通过「Static Structural → Steady-State Thermal」的链接传递接触压力，并引用CMY模型的$h_c(P)$表格。这是螺栓连接部位热路径设计必备的工作流程。
TIM的建模
所以硅脂的涂抹方式和压力很重要呢。
没错。即使是相同的硅脂，根据装配条件，实际的 $R_{TIM}$ 也会有2〜5倍的变化。
加压与TCR的实验相关性
Cooper-Mikic-Yovanovich（CMY）模型（1969年）根据接触压力p和表面粗糙度σ预测接触传热系数hc。对于金属间接触，压力加倍通常会使hc增加约1.5倍。该模型被ISO/TS 22007引用，并成为CPU散热器安装压力设计的理论基础。
线性单元 vs 二次单元
热传导分析中，线性单元通常也能获得足够的精度。温度梯度陡峭的区域（如热冲击）推荐使用二次单元。
热流评估
根据单元内的温度梯度计算。有时需要像节点应力那样进行平滑处理。
对流-扩散问题
当佩克莱数高（对流主导）时，需要迎风稳定化（如SUPG）。纯热传导问题则不需要。
非稳态分析的时间步长
时间步长应远小于热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$: 热扩散率）。对于急剧的温度变化，自动时间步长控制有效。
非线性收敛
由温度相关物性值引起的非线性通常较温和，Picard迭代（直接替换法）通常足够。对于辐射的强非线性，推荐使用牛顿法。
稳态分析判定
当所有节点的温度变化低于阈值（如 $|\Delta T| / T_{max} < 10^{-5}$）时，判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算快，但时间步长大了就不稳定（会错过风暴）。隐式法是“也考虑未来状态的预测”——即使时间步长大也能稳定，但每一步都需要解方程，计算量大。对于没有急剧温度变化的问题，使用隐式法配合较大的时间步长更高效。
实践指南
测量方法
如何获取接触热阻的值？
最可靠的是实测。按照ASTM D5470标准的稳态法测量。
1. 在两个铜块之间夹入TIM样品
2. 一侧加热，另一侧冷却
3. 从各铜块内的温度梯度外推界面温差
4. $R_{TIM} = \Delta T_{interface} / q''$
市售TIM数据表上的值可以直接用吗？
需要注意。数据表上的值通常是在理想条件（高压力、完全润湿）下测量的。实际装配条件（压力、表面粗糙度、涂布量）下的值可能会恶化1.5〜3倍。初期设计时，按数据表值的50%估算比较安全。
螺栓连接部位的热设计
以散热板的螺栓连接为例。M4螺栓4根，紧固扭矩1.5 Nm 时：
- 螺栓轴向力: 约2600N/根
- 接触面压: 座面（φ8mm）处约52 MPa
- 根据CMY模型，$h_c \approx 8000$ W/(m2K)（铝・抛光面）
可以从接触面压推算$h_c$呢。
是的。但螺栓之间的区域面压会降低，因此更精确的做法是先通过结构分析求得面压分布，再将$h_c$映射上去。
TIM的热导率与实际TCR
Intel的Core i9-13900K（2022年）在CPU芯片与集成散热器（IHS）之间采用了InFusion液态金属TIM，将接触热阻相比固体硅脂降低了约50%。液态金属（GaInSn系）的热导率约为40 W/m·K，远高于顶级硅脂（约12 W/m·K）。但因其会腐蚀铝制IHS，故已改用铜制。
分析流程的比喻
热分析的流程可以想象成“浴缸的循环加热设计”。确定浴缸形状（分析对象），设定初始水温（初始条件）和室外气温（边界条件），调整循环加热的功率（热源）。计算预测“2小时后会不会变凉？”——这就是非稳态热分析的本质。
初学者容易掉入的陷阱
“可以忽略辐射吗？”——室温附近通常可以。但超过几百度就另当别论了。辐射传热与温度的四次方成正比，高温时会压倒对流。您有过晴天时向阳处和背阴处体感温度完全不同的经验吧？那就是辐射的威力。在工业炉或发动机周围的分析中忽略辐射，就像在酷暑天坚持说“阳光没关系”一样。
相关主题
なった
詳しく
報告