接触热阻
接触热阻的理论基础
什么是接触热阻
老师，两块金属接触在一起也会在界面处出现温度跳变吗？
是的。实际固体表面从微观上看凹凸不平，真实接触面积仅为表观面积的1～5%左右。剩余间隙中充满空气或润滑脂。这种不完全接触产生温度不连续。这就是接触热阻。
只有几个百分点在接触？
对。所以接触热阻 $R_c$ 定义为界面温度跳变 $\Delta T$ 与热流密度 $q''$ 的比值。
单位为 [m$^2$ K/W]。其倒数为接触导热系数 $h_c = 1/R_c$ [W/(m$^2$ K)]。
支配因子
接触热阻的大小由什么决定？
主要有4个因子。
| 因子 | 影响 | 典型值 |
|---|---|---|
| 表面粗糙度 $\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年代被精密设备设计广泛采用。
接触热阻的数值计算方法
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。被广泛用于半导体器件封装评估。
接触热阻的实务应用
解析流程
请告诉我包含接触热阻的解析步骤。
标准流程是这样的。
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}$ | 螺栓附近 vs 远处变动 |
导热润滑脂能改善两个数量级。
因为润滑脂用高导热系数的材料替代空气（k=0.026 W/(m K)）来填充间隙。但也要考虑润滑脂的长期劣化（泵出、干燥），长期可靠性评估应使用劣化后的值。
结果验证
接触热阻解析的验证重点是这样的。
- 温度跳变的确认: 绘制接触面两侧的节点温度，检查是否满足 $\Delta T = q'' \cdot R_c$
- 能量守恒: 通过接触面的热量在上下游是否一致
- 灵敏度分析: 改变 $h_c$ 的±50%来量化结果敏感度
灵敏度分析要做到什么程度？
最少要检查 $h_c$ 的上下限如何影响最高温度。变化在规范范围内就没问题。超出规范时建议获取实测值。
EV逆变器的结点热阻管理
特斯拉Model 3（2017年）的SiC逆变器中，IGBT模块底面与散热器的接触热阻被严格控制在0.05 K/W以下，采用相变TIM薄片。生产线上进行了扭矩管理的螺栓紧固，表面粗糙度精加工到Ra<0.8 µm以下，配合使用保证批量品质。
接触热阻的软件比较
商用工具的实现
处理接触热阻时，各个软件怎样设置？
让我比较一下主要工具的设置方法。
| 工具 | 设置方法 | 关键词/GUI操作 |
|---|---|---|
| ANSYS Mechanical | Contact单元指定TCC值 | KEYOPT(1)=1为热接触，TCC值由realConstant定义 |
| Abaqus | *GAP CONDUCTANCE | 压力相关表由*SURFACEINTERACTION定义 |
| COMSOL | Thermal Contact节点 | Heat Transfer > Thermal Contact，$h_c$ 可用表达式输入 |
| ANSYS Icepak | Interface Resistance | 通过GUI直接输入部件间接触阻 |
ANSYS的APDL怎么写？
接触对设置的例子是这样的。
```
! 接触对定义（热解析）
ET,2,TARGE170
ET,3,CONTA174
KEYOPT,3,1,1 ! 热接触
R,2,,,,,TCC,5000 ! h_c = 5000 W/(m2K)
```
Abaqus可以设置压力相关的间隙导热系数表。
```
*GAP CONDUCTANCE
5000., 0.0, 1.0E6 ! h_c=5000 at clearance=0, P=1MPa
100., 0.001, 0.0 ! h_c=100 at clearance=1mm, P=0
```
能用压力相关表很方便。
COMSOL可在GUI上直接用数式输入 $h_c(P,T)$，能直接实现CMY模型，灵活性很高。
工具选择
| 用途 | 推荐工具 |
|---|---|
| 结构-热耦合的压力相关接触阻 | ANSYS Mechanical, Abaqus |
| 电子设备界面阻 | ANSYS Icepak, FloTHERM |
| 自定义模型实现 | COMSOL |
电子设备用Icepak或FloTHERM更方便。
Icepak的TIM（热界面材料）库很充实，所以很多情况下不需要逐个细致设置接触阻。
ANSYS Icepak接触阻的自动计算
ANSYS Icepak在2016年的18.0版本中，增加了基于JEDEC JESD51-14标准自动计算部件-PCB间接触热阻的功能。在此之前用户必须从数据表手工输入，对于100个以上部件的电路板来说输入工作量曾是个问题。
接触热阻的先进研究
多尺度接触模型
最近的研究采用什么方法？
关注表面粗糙度分形性的多尺度模型备受关注。Majumdar-Bhushan（MB）模型用分形维数 $D$ 和表面参数 $G$ 来描述接触热阻。
其中 $a_L$ 是最大接触斑点面积，$A_n$ 是表观接触面积，$a_s$ 是最小斑点面积。
用分形描述表面很有意思。
相比传统高斯粗糙度模型，在更广条件范围内精度更高。不过需要测量 $D$ 和 $G$，白色干涉仪或AFM进行表面形状测量是前提。
方向性接触热阻
非对称表面对（如研磨面和喷砂面）会因热流方向而产生接触热阻的整流效应（热整流）。温度差较大时，材料热膨胀差导致真实接触面积变化也是方向依存性的原因。
像二极管一样只允许热在一个方向流动？
正是这样。热二极管在航天器热控研究中进展快。放热时阻低，入热时阻高，只让热向一个方向逃逸。
纳米尺度接触
接触斑点进入纳米尺度时，声子弹道输运占主导，傅立叶定律失效。半导体器件的芯片-封装界面等处用声学失配模型（AMM）或漫散射失配模型（DMM）。
其中 $c_j$ 是比热容，$v_j$ 是声子群速度，$\alpha_{1\to 2}$ 是透射率。
到声子级别的考虑工作很大。
实务上先用DMM粗估，用NEMD（非平衡分子动力学）仿真验证的方法被采用。
真空中辐射引起的并联传热
航天器的接触热阻不仅需考虑固体传导，接触面间的微小间隙中的辐射也要作为并联路径。国际空间站太阳电池板臂中，NASA喷气推进实验室在2005年报告接触部辐射分量占全热传递最多15%，导致低温低压环境设计标准修订。
接触热阻的故障排除
常见问题及应对
接触热阻解析遇到问题时，应该怀疑什么？
让我整理常见问题。
1. 界面温度连续（无跳变）
原因: 接触单元导热系数过大，或节点已合并。$h_c = 10^{10}$ 这样的值实际上等于完全接触。
对策: 改正为物理上合理的 $h_c$ 值。检查节点是否有合并。
2. 界面温度差与实测值偏差大
解析显示 $\Delta T = 2$ ℃但实测有10℃这样的情况。
检查清单如下。
| 检查项目 | 常见原因 |
|---|---|
| $h_c$ 的值 | 文献值与实际表面状态偏差 |
| TIM涂布不均 | 局部残留空气层 |
| 螺栓紧固力 | 扭矩不足导致接触压力低 |
| 辐射贡献 | 高温时存在辐射热传递 |
| 经年劣化 | 润滑脂泵出或干燥 |
3. 不收敛
原因: 压力相关的 $h_c$ 表中零压力时导热系数为零。接触离开时变绝热，产生尖峰温度。
对策: 在零压力处也设最小导热系数（如 $h_c = 10$ W/(m$^2$ K)）。物理上通过辐射或空气层也存在微小热传递。
4. 接触对设置错误
主面/从面的指定有影响吗？
影响很大。原则上粗网格侧为主，细网格侧为从。反过来会产生穿透，温度精度下降。ANSYS中刚体侧用TARGE170，变形体侧用CONTA174是基本原则。
接触热阻是看不见的现象，更需要仔细验证。
与实测比较是最可靠的验证。用热像仪拍摄界面附近的温度分布，与解析结果叠加检查是工业上常用做法。
螺栓重新紧固改善TCR
在半导体制造设备的晶圆舞台中，因出现温度不均而排查发现舞台裏面螺栓因热循环而松动。重新紧固后TCR降低约40%，面内温度差从±5℃改善到±1.5℃。将紧固扭矩确认加入定期维护检查清单后防止了再发生。
价值
详细
错误