接触熱抵抗
理论与物理
什么是接触热阻
老师，两块金属贴在一起，界面处温度真的会跳跃吗？
是真的。实际的固体表面在微观下凹凸不平，真实接触面积只有表观面积的1~5%左右。剩余的间隙中充满了空气或导热膏。这种不完全接触导致了温度不连续。这就是接触热阻。
只有百分之几的面积接触吗？
是的。因此接触热阻 $R_c$ 定义为界面温度跳跃 $\Delta T$ 与热流密度 $q''$ 的比值。
单位是 [m$^2$ K/W]。其倒数是接触热导 $h_c = 1/R_c$ [W/(m$^2$ K)]。
主导因素
接触热阻的大小由什么决定？
主要有四个因素。
| 因素 | 影响 | 典型值 |
|---|---|---|
| 表面粗糙度 $\sigma_s$ | 越粗糙 $R_c$ 越大 | Ra 0.1〜10 μm |
| 接触压力 $P$ | 越高 $R_c$ 越小 | 0.1〜100 MPa |
| 界面介质 | 使用导热膏可将 $R_c$ 降低至1/10 | — |
| 材料硬度 $H$ | 越软真实接触面积越大 | — |
Cooper-Mikic-Yovanovich模型被广泛使用。
其中 $k_s$ 是调和平均热导率 $k_s = 2k_1 k_2/(k_1+k_2)$，$m$ 是表面斜率，$\sigma$ 是合成粗糙度 $\sigma = \sqrt{\sigma_1^2 + \sigma_2^2}$，$H$ 是显微硬度。
也就是说几乎与接触压力的0.95次方成正比。螺栓紧固力很关键啊。
正是如此。在电子设备散热片安装中，螺栓的扭矩管理直接关系到温度。有时将紧固扭矩加倍，接触热导也会增加约一倍。
微观接触模型的三代演进
接触热阻的理论模型经历了三代演进：(1)GreenWood-Williamson（1966年：弹性接触）→(2)Majumdar-Bhushan（1991年：分形表面）→(3)Persson（2001年：全尺度弹性）。Persson模型源于汽车轮胎摩擦领域，但在热接触方面的应用受到关注，并在2010年代被精密设备设计广泛采用。
各项的物理意义
- 蓄热项 $\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中如何处理接触热阻？
在界面处设置一个薄的虚拟层，或者直接在接触单元上指定热导值。基本的离散化通过界面节点间的热导矩阵表示。
在单元级别组装后，得到界面的对流传热矩阵。
和对流边界条件的Robin条件在数学上是相同的形式呢。
很好的观察。在实现上也与Robin条件处理相同。但需要额外进行接触面的配对（主/从面的对应关系）。
间隙热导的设置
在实际工作中，根据接触状态使用可变热导的情况很多。
| 接触状态 | 热导 $h_c$ [W/(m$^2$ K)] | 应用场景 |
|---|---|---|
| 完全接触 | $10^5$〜$10^6$ | 焊接部、过盈配合 |
| 导热膏填充 | $10^3$〜$10^4$ | 散热片安装 |
| 金属间直接接触 | $10^2$〜$10^4$ | 螺栓紧固面 |
| 存在空气间隙 | $10^0$〜$10^2$ | 松配合 |
有4个数量级以上的差异呢。这里的估算如果出错，结果会完全不同吧。
正是如此。常规做法是通过灵敏度分析，将 $h_c$ 设为一半和两倍，确认温度差异。如果结果强烈依赖于 $h_c$，则应考虑获取实测值。
非线性接触热阻
压力依赖或温度依赖的接触热阻会成为非线性问题。通过结构分析耦合求解接触压力，根据该压力用CMY模型计算 $h_c$，再反馈给热分析。重复此迭代直至收敛。
需要结构和热的耦合分析呢。
Ansys Mechanical或Abaqus中已将此耦合自动化。在Ansys中，只需在接触单元CONTA174/TARGE170上设置TCC（热接触导率）即可。
用激光闪光法实测
热接触阻力的精密测量采用ASME D5470标准的激光闪光法。该方法通过层叠样品并改变压力来反算接触面的TCR，Netzsch LFA 467（2020年发售）能以±2%精度测量0.5〜50 mm²·K/W的范围。广泛应用于功率半导体封装的评估。
线性单元 vs 二次单元
在热传导分析中，线性单元通常足以获得足够的精度。在温度梯度陡峭的区域（如热冲击），推荐使用二次单元。
热流密度的评估
根据单元内的温度梯度计算。与节点应力类似，有时需要进行平滑处理。
对流-扩散问题
当佩克莱数较高（对流主导）时，需要迎风稳定化（如SUPG）。纯热传导问题则不需要。
瞬态分析的时间步长
时间步长应远小于热扩散的特征时间 $\tau = L^2 / \alpha$（$\alpha$: 热扩散率）。对于急剧的温度变化，自动时间步长控制有效。
非线性收敛
由温度依赖物性值引起的非线性通常较温和，Picard迭代（直接替换法）通常足够。对于辐射的强非线性，推荐牛顿法。
稳态分析的判定
当所有节点的温度变化低于阈值（如 $|\Delta T| / T_{max} < 10^{-5}$）时判定为收敛。
显式法与隐式法的比喻
显式法是“仅凭当前信息预测未来的天气预报”——计算快，但时间步长大时不稳定（会漏掉风暴）。隐式法是“考虑未来状态的预测”——即使时间步长大也能稳定，但每个步长都需要解方程，计算量大。对于没有急剧温度变化的问题，使用隐式法配合较大的时间步长更高效。
实践指南
分析流程
请告诉我包含接触热阻的分析步骤。
标准流程如下。
1. 识别接触面: 在CAD上明确接触的面组
2. 确定热导值: 根据实测值、文献值或CMY模型计算
3. 设置接触对: 定义主/从面并分配 $h_c$
4. 网格协调: 调整网格使接触面上的节点相对
5. 求解·验证: 确认界面温度跳跃在物理上是否合理
网格协调是必须的吗？
非协调网格也能求解（如Mortar法或GGI连接），但对于接触面，协调网格通常精度更高。特别是当温度跳跃较小时，网格不协调产生的伪影可能会显现出来。
典型的接触热阻值
| 界面 | $R_c$ [m$^2$ K/W] | 条件 |
|---|---|---|
| Al-Al（抛光面，无导热膏） | $2 \times 10^{-4}$ | P=1 MPa |
| Al-Al（涂有导热膏） | $5 \times 10^{-6}$ | k=5 W/(m K) 导热膏 |
| Cu-Cu（抛光面） | $1 \times 10^{-4}$ | P=1 MPa |
| Si-散热片（带TIM） | $1 \times 10^{-5}$ | TIM厚度50μm |
| 螺栓紧固法兰 | $10^{-4}$〜$10^{-3}$ | 螺栓附近与远处有波动 |
使用导热膏能改善两个数量级啊。
因为导热膏填充了原本充满空气（k=0.026 W/(m K)）的间隙。但需要考虑导热膏的长期劣化（泵出、干燥），在长期可靠性评估中应使用劣化后的值。
结果的验证
接触热阻分析的验证要点如下。
- 确认温度跳跃: 绘制接触面两侧的节点温度，检查是否与 $\Delta T = q'' \cdot R_c$ 一致
- 能量守恒: 通过接触面的热量在上游和下游是否一致
- 灵敏度分析: 将 $h_c$ 变动±50%，量化对结果的影响
なった
詳しく
報告