接触热阻
接触热阻的理论基础
接触热阻是什么
为什么将两个固体紧压在一起，界面处温度仍会发生跳跃？
从显微镜级别看，接触面由于表面粗糙度只在点接触（真实接触点）处接触。实际接触面积仅为表观面积的1%左右，其余为空气间隙。这种接触不完全性产生温度跳跃。
基本公式
接触热阻的定义公式如下。
接触导热系数为 $h_c = 1/R_c$ [W/(m²K)]。
典型值是多少？
取决于条件，差异很大。
| 接触条件 | $h_c$ [W/(m²K)] |
|---|---|
| 铝同对、研磨面、高压力 | 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年建立了第一个系统性的实测数据库。今天仍将该数据用作设计初期的参考值。
接触热阻的数值计算方法
FEM中的接触热阻建模
在FEM中如何表达接触热阻？
主要有三种方法。
| 方法 | 概述 | 精度 |
|---|---|---|
| 间隙导热系数 | 在面间设置 $h_c$ | 实用 |
| 薄层单元 | 用等效k、t将TIM建模 | 需要厚度 |
| 压力相关导热系数 | $h_c(P)$ 通过结构耦合算出 | 高精度 |
在Ansys Mechanical中通过contact element的TCC (Thermal Contact Conductance)设置。在Abaqus中可通过*GAP CONDUCTANCE关键词定义压力相关表格。
如果做成压力相关的话需要结构分析耦合吧。
是的。在Ansys Workbench中通过"Static Structural → Steady-State Thermal"的链接转移接触压力，参考CMY模型的$h_c(P)$表。这是螺栓紧固部位热路径设计的必要工作流程。
TIM的建模
是体导热阻与两个表面接触阻之和。对于热界面润滑脂，体导热系数k=3 W/(mK)、厚度50um时体导热阻较小，通常两个表面的接触阻起支配作用。
所以润滑脂的涂抹方式和压力很重要啊。
没错。同一种润滑脂，根据实装条件不同，实效的 $R_{TIM}$ 会有2～5倍的差异。
加载力和TCR的实验相关性
Cooper-Mikic-Yovanovich（CMY）模型（1969年）能根据接触压力p和表面粗糙度σ预测接触热传导率hc。在金属间接触中，压力翻倍时hc约增加1.5倍。该模型被参考在ISO/TS 22007中，成为CPU散热器安装压力设计的理论依据。
接触热阻的实务应用
测量方法
接触热阻的值如何获取？
最可靠的是实测。按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/(m²K)（铝·研磨面）
能从接触面压推算$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。
接触热阻的软件比较
按工具的设置方法
请告诉我各软件中接触热阻的设置方法。
主要工具的设置方法如下。
| 工具 | 设置位置 | 关键词 |
|---|---|---|
| Ansys Mechanical | Contact Region > Thermal Conductance | 直接输入TCC值 |
| Abaqus | *GAP CONDUCTANCE | 可设压力相关表 |
| COMSOL | Thermal Contact (Pair) | 也可用薄层近似 |
| Ansys Icepak | Resistance设置 | 单位面积热阻 |
| FloTHERM | Interface Resistance | SmartPart接触设置 |
Icepak和FloTHERM输入形式不同吗？
Icepak是$R$ [m²K/W]，FloTHERM也是。Ansys Mechanical输入$h_c$ [W/(m²K)]，所以要注意转换 $h_c = 1/R$。倒数弄反是常见错误。
压力传导耦合的自动化
Ansys Workbench的自动化流程：
1. Static Structural建模螺栓紧固
2. 获取接触面压分布
3. 通过ACT Extension或APDL片段应用 $h_c = f(P)$
4. 用Steady-State Thermal求解温度场
用Python脚本全自动化的话，可以高效进行螺栓数量或紧固扭矩的参数化研究。
能做结构热的双向耦合吗？
一般单向就够了。接触压力几乎不受温度影响（除了线膨胀导致的微小变化），所以结构→热的单向传递就没问题。
Mentor Flootherm的接触热阻模型
Mentor Graphics（现为Siemens EDA）的FloTHERM XT（2014年）在PCB实装部件间的接触热阻中采用"JEDEC JESD51"标准模型的自动生成功能。在此之前，用户需要单独输入TCR值，遗漏输入导致分析错误频繁发生。
接触热阻的先端研究
纳米级接触热阻
在纳米级接触热阻的物理会改变吗？
当声子平均自由程（Si约40nm）与接触点尺寸接近时，基于傅里叶定律的连续体模型失效。需要考虑弹道声子输运的Landauer形式的接触导纳模型。
$M$ 是通过界面的声子模式数。
时间相关接触阻
接触面的长期变化（氧化、蠕变、磨损）会使$h_c$随时间变化。电源模块的焊接点因热循环导致空隙成长，$R_c$逐渐增大。
可靠性设计需要预测长期劣化啊。
是的。用Coffin-Manson规律或Norris-Landzberg式预测寿命，建模空隙率增加导致的$R_c$增大。与Ansys Sherlock等可靠性工具联动很有效。
相变TIM（PC-TIM）
在工作温度下熔融并变薄的PC-TIM存在泵出（熔融TIM流出）问题。通过CFD分析预测熔融TIM的流动行为，评估长期可靠性的研究在进行中。
TIM的世界也很深啊。
下一代TIM如液体金属（镓合金，k=20～30 W/(mK)）、碳纳米管阵列、石墨烯薄片等在研究中。界面阻的低减都是关键。
声子与电子的接触热传导
在金属间接触中，热主要由电子和声子承载，但在穿越表面间微小空隙时输运机制会切换。MIT的Evan等人2009年发表的AMD（声学失配）和DMD（扩散失配）模型的统合理论表明，在Cu-Al接触中声子成分的贡献约为电子的3倍。该知识被用于功率器件的长期可靠性设计。
接触热阻的故障排除
常见问题
接触热阻设置容易出错的地方是什么？
1. R与hc与G的混淆
问题：文献和数据表中R [m²K/W]、hc [W/(m²K)]、G [W/K]混淆，单位转换出错。
- $R$：单位面积热阻 [m²K/W]
- $h_c = 1/R$：接触导热系数 [W/(m²K)]
- $G = h_c \cdot A$：总导热系数 [W/K]
对策：统一为SI单位，输入前进行量纲检查。
2. 接触面积评估过高
直接用CAD面积不行吗？
问题：在实际螺栓紧固面中面压不均匀，有效接触面积小于CAD面积。特别是螺栓间距大时，中央部分实际上是非接触状态。
对策：用结构分析求面压分布，只把面压>0.1MPa的区域作为有效接触面积。
3. 分析值与实测值的偏差
接触热阻是不确定性最大的参数。分析值与实测值偏差2倍是常见的。
对策：必须做$h_c$±50%的变动敏感性分析，评估对温度的影响。如影响大则考虑实测校准。
最后接触热阻是最不确定的啊。
正是。可以说热分析的精度由接触条件的设置决定。积累实测数据是改善长期分析精度的关键。
表面氧化膜急剧增加TCR
铝合金接触面形成的氧化膜（Al₂O₃，λ≈1 W/m·K）会使TCR相比清洁面增加5～10倍。NASA詹森航天中心2003年发表的调查显示，航天器长期运行中螺栓紧固部位的TCR发生经年变化，导致预期外的温度上升。
更详细
错误