工程科学学报，第41卷，第10期：1240-1248,2019年10月 Chinese Journal of Engineering,Vol.41,No.10:1240-1248,October 2019 D0I:10.13374/j.issn2095-9389.2018.10.05.001;http:/journals.ustb.ed.cn 宏观接触热阻研究综述 汪献伟)四，王兆亮2)，何庆”，李秀莲) 1)江苏理工学院机械工程学院，常州2130012)中国科学院等离子体物理研究所，合肥230031 ☒通信作者，E-mail:wangxw(@jsut.cdu.cm 摘要基于国内外学者在宏观接触热阻领域的研究成果，阐述了宏观接触热阻的理论研究和工程应用现状，介绍了宏观接 触热阻的理论计算方法、实验测量方法和数值模拟方法，明确了各种方法的优缺点.简述了宏观接触热阻的影响因素，基于中 国聚变工程实验堆(C℉ETR)低温超导线圈的降温实验，创建界面接触热阻分析模型，详细分析了热流方向、温度、压力等因素 对超导磁体不锈钢铠甲、介电绝缘材料、失超保护材料接触热阻的影响，并进一步探究了温度、压力等因素导致接触热阻发生 变化的热力学原因.最后，结合工程计算中对准确性和便捷性的要求，指出了接触热阻的未来研究方向. 关键词接触热阻：接触热导：微凸体：粗糙度：热流通量 分类号TL62·2 Research overview of macroscopic thermal contact resistance WANG Xian-wei,WANG Zhao-liang?),HE Qing),LI Xiu-lian) 1)School of Mechanical Engineering,Jiangsu University of Technology,Changzhou 213001,China 2)Institute of Plasma Physics,Chinese Academy of Sciences,Heifei 230031,China Corresponding author,E-mail:wangxw@jsut.edu.cn ABSTRACT Given the limitations of manufacturing precision,manufacturing cost,and other factors,contact clearance is always in- evitable at the contact interface of components.The contact clearance leads to the reduction of heat flux during the heat transfer process.The effect of thermal contact resistance is significant,particularly in the fields of aerospace,microelectrical technology,and cryogenic superconductor that are closely related to the operating temperature.The thermal contact resistance is affected by many fac- tors,such as size,shape,space of asperity,mechanical properties of the material,external pressure,and temperature.Moreover, these factors usually interact with each other and need to be coupled.Thus,how to describe the thermal contact resistance accurately and build the appropriate prediction model are the key problems that should be resolved during engineering calculations.On the basis of the research results of domestic and international scholars,the current research state of thermal contact resistance during theoretical cal- culations and engineering applications were presented.The theoretical calculation,experimental measurement,and digital simulation methods used to analyze macroscopic thermal contact resistance were summarized,and the advantages and disadvantages of these meth- ods were indicated.The effects of different factors on macroscopic thermal contact resistance were briefly discussed.On the basis of the combined cooling experiments of low-temperature superconducting magnet coils of the China Fusion Engineering Test Reactor,the effects of heat flow direction,temperature,and pressure on the thermal contact resistance of superconducting magnet components,such as stainless-steel jacket,dielectric insulation material,and quench protection material,were analyzed.Moreover,the reason why the effects of temperature and external pressure result in the change of thermal contact resistance was investigated from the perspective of thermal mechanics.Finally,given the accuracy and convenience requirements for the calculation of thermal contact resistance during engineering practice,the future research direction was indicated. KEY WORDS thermal contact resistance;thermal contact conductance;asperity;roughness;heat flux 收稿日期：2018-10-05 基金项目：国家磁约束聚变能发展研究专项资助项目(2014GB105002)
工程科学学报,第 41 卷,第 10 期:1240鄄鄄1248,2019 年 10 月 Chinese Journal of Engineering, Vol. 41, No. 10: 1240鄄鄄1248, October 2019 DOI: 10. 13374 / j. issn2095鄄鄄9389. 2018. 10. 05. 001; http: / / journals. ustb. edu. cn 宏观接触热阻研究综述 汪献伟1) 苣 , 王兆亮2) , 何 庆1) , 李秀莲1) 1) 江苏理工学院机械工程学院, 常州 213001 2) 中国科学院等离子体物理研究所, 合肥 230031 苣通信作者, E鄄mail: wangxw@ jsut. edu. cn 摘 要 基于国内外学者在宏观接触热阻领域的研究成果,阐述了宏观接触热阻的理论研究和工程应用现状,介绍了宏观接 触热阻的理论计算方法、实验测量方法和数值模拟方法,明确了各种方法的优缺点. 简述了宏观接触热阻的影响因素,基于中 国聚变工程实验堆(CFETR)低温超导线圈的降温实验,创建界面接触热阻分析模型,详细分析了热流方向、温度、压力等因素 对超导磁体不锈钢铠甲、介电绝缘材料、失超保护材料接触热阻的影响,并进一步探究了温度、压力等因素导致接触热阻发生 变化的热力学原因. 最后,结合工程计算中对准确性和便捷性的要求,指出了接触热阻的未来研究方向. 关键词 接触热阻; 接触热导; 微凸体; 粗糙度; 热流通量 分类号 TL62 + 2 收稿日期: 2018鄄鄄10鄄鄄05 基金项目: 国家磁约束聚变能发展研究专项资助项目(2014GB105002) Research overview of macroscopic thermal contact resistance WANG Xian鄄wei 1) 苣 , WANG Zhao鄄liang 2) , HE Qing 1) , LI Xiu鄄lian 1) 1) School of Mechanical Engineering, Jiangsu University of Technology, Changzhou 213001, China 2) Institute of Plasma Physics, Chinese Academy of Sciences, Heifei 230031, China 苣Corresponding author, E鄄mail: wangxw@ jsut. edu. cn ABSTRACT Given the limitations of manufacturing precision, manufacturing cost, and other factors, contact clearance is always in鄄 evitable at the contact interface of components. The contact clearance leads to the reduction of heat flux during the heat transfer process. The effect of thermal contact resistance is significant, particularly in the fields of aerospace, microelectrical technology, and cryogenic superconductor that are closely related to the operating temperature. The thermal contact resistance is affected by many fac鄄 tors, such as size, shape, space of asperity, mechanical properties of the material, external pressure, and temperature. Moreover, these factors usually interact with each other and need to be coupled. Thus, how to describe the thermal contact resistance accurately and build the appropriate prediction model are the key problems that should be resolved during engineering calculations. On the basis of the research results of domestic and international scholars, the current research state of thermal contact resistance during theoretical cal鄄 culations and engineering applications were presented. The theoretical calculation, experimental measurement, and digital simulation methods used to analyze macroscopic thermal contact resistance were summarized, and the advantages and disadvantages of these meth鄄 ods were indicated. The effects of different factors on macroscopic thermal contact resistance were briefly discussed. On the basis of the combined cooling experiments of low鄄temperature superconducting magnet coils of the China Fusion Engineering Test Reactor, the effects of heat flow direction, temperature, and pressure on the thermal contact resistance of superconducting magnet components, such as stainless鄄steel jacket, dielectric insulation material, and quench protection material, were analyzed. Moreover, the reason why the effects of temperature and external pressure result in the change of thermal contact resistance was investigated from the perspective of thermal mechanics. Finally, given the accuracy and convenience requirements for the calculation of thermal contact resistance during engineering practice, the future research direction was indicated. KEY WORDS thermal contact resistance; thermal contact conductance; asperity; roughness; heat flux
汪献伟等：宏观接触热阻研究综述 ·1241· 由于机械加工精度的限制，两固体表面的接触 断涌现，近年来不少学者基于已有的研究成果，在机 往往只发生在一些离散的局部区域，实际接触面积 械制造、航空航天、微电子等诸多领域开展了接触热 远小于名义接触面积，从而影响了界面之间的热传 阻的应用研究，比如：基于接触热阻对笔记本电脑的 递.在微电子领域，电子产品不断朝小型化、轻量化 散热性能进行分析优化]：接触热阻对双H型鳍 方向发展，如何在产品尺寸、重量减小的同时保证热 片管传热系数的影响[2]：接触热阻对高速精密电主 量的及时耗散：在低温超导设备中，超导材料对温度 轴热特性的影响9]：以及接触热阻对燃料电池温度 具有高度灵敏性，如何避免接触界面因温度跳变引 场分布的影响[30]等 发的失超破坏：在航空航天领域如何降低一些接收 部件因表面粗糙度导致的局部热梯度.为了将温度 2接触热阻的理论及试验研究 控制在合理的范围之内，一方面可以通过设计特定 2.1单点接触热阻的预测模型 的降温系统将热量及时耗散掉，另一方面需要采取 如果把离散的局部接触简化为一个个的点，从 必要的措施降低零部件接触表面的温度梯度.当前 而开展接触热阻研究称之为点理论，接触热阻的点 普遍采用的方案是改变材料接触界面的接触状态， 理论分为单点接触理论和多点接触理论.无论接触 减小接触界面的接触热阻.所谓接触热阻是指由接 面的几何形貌如何复杂，总的接触热阻总可以简化 触表面之间的接触间隙引起的热流线收缩而产生的 为单点接触热阻的并联，单点接触热阻研究的核心 附加阻)].开展接触热阻研究对于建立接触热 内容是确定接触点的形状并获取接触点的尺寸信 阻的经验公式，揭示接触传热的物理机理，解决实际 息，其理想的传热模型为接触点与半无限大空间或 的工程技术问题具有重要的现实意义.接触热阻的 半无限长流管之间的接触换热 研究主要包括基础研究和工程应用，基础研究主要 对单点接触热阻的计算通常将接触，点简化为圆 从宏观与微观方面开展理论、试验、数值模拟等方面 柱（圆盘）、半球体、截锥体等模型.Gibson3] 的预测计算，探索可用于工程实践的接触热阻半经 Fenech与Rohsenow】、Cooper等]建立了圆盘计 验关联式以及接触热阻的影响因素.接触热阻的工 算模型；Bahrami等3-]建立了半球体计算模型： 程应用主要是在接触间隙填充导热性能各异的填充 Willimams与Major[6]建立了截锥体计算模型； 材料来实现减小或增大接触热阻的目的.本文重点 Yovanovich!详细比较了不同形状接触点的接触换 综述宏观接触热阻的影响因素和国内外学者在基础 热系数.对于圆盘计算模型如图1(b)所示，假定接 研究领域取得的成果，并对各种预测模型和计算方 触点的实际半径为a,名义接触半径为b,名义接触 法进行评价，指出进一步的研究重点.对于微观接 面积为A,两接触材料的导热率分别为入1和入2，则 触热阻的研究以及新型的表征方法和技术可参阅文 接触点的接触热阻R[38]为 献[3]. R=AFx) (1) 1接触热阻国内外研究进展 Ta入. 11.1 (2) 接触热阻研究最早起源于Kapitza热阻4]，随后 元元+元 众多学者对接触热阻进行了广泛的理论和试验研 f(x）=1-1.40918x+0.33801x3+0.06792x+… 究s-io).Lambert与Fletcher【-iay,Yovanovich与 (3) Marotta2]对接触热阻的预测模型和基本机理进行了 式中，x为实际接触半径与名义接触半径的比值. 研究：Fletcher),Kraus与Bar-Cohen[对接触热 对于当量球模型如图1(a)所示，Greenwood与 阻进行了试验研究；Madhusudana与Fletchertis、 Williamson!39]、Yovanovich[ao]、Burde与Yovanov- Fletchert6],Swartz与Pohl[o)和Gmelin等[对接触 ich]首先把接触热阻细分为宏观接触热阻、微观 热阻的研究进展做了综述.国内学者对接触热阻的 接触热阻和薄膜接触热阻.其中薄膜接触热阻可以 研究起步较晚，90年代南京航空航天大学顾慰兰开 通过实验来进行评估，宏观接触热阻通过当量尺寸 发了一套接触热阻测试装置，研究了温度和接触压 的单个球型接触点进行计算，假定接触面的实际接 力对接触热阻的影响[18-)]：徐烈课题组[20-21]和王惠 触半径和名义接触半径分别为α，b,两接触体的热 龄课题组[2-2]对低温真空环境下的接触热阻进行 导率分别为k,k2,则交界面的宏观接触热阻为 了研究.此外，应济等2-2)和韩玉阁等]也对接触 R=(x) (4) 热阻进行了探索研究.随着接触热阻研究成果的不 2ak
汪献伟等: 宏观接触热阻研究综述 由于机械加工精度的限制,两固体表面的接触 往往只发生在一些离散的局部区域,实际接触面积 远小于名义接触面积,从而影响了界面之间的热传 递. 在微电子领域,电子产品不断朝小型化、轻量化 方向发展,如何在产品尺寸、重量减小的同时保证热 量的及时耗散;在低温超导设备中,超导材料对温度 具有高度灵敏性,如何避免接触界面因温度跳变引 发的失超破坏;在航空航天领域如何降低一些接收 部件因表面粗糙度导致的局部热梯度. 为了将温度 控制在合理的范围之内,一方面可以通过设计特定 的降温系统将热量及时耗散掉,另一方面需要采取 必要的措施降低零部件接触表面的温度梯度. 当前 普遍采用的方案是改变材料接触界面的接触状态, 减小接触界面的接触热阻. 所谓接触热阻是指由接 触表面之间的接触间隙引起的热流线收缩而产生的 附加阻力[1鄄鄄2] . 开展接触热阻研究对于建立接触热 阻的经验公式,揭示接触传热的物理机理,解决实际 的工程技术问题具有重要的现实意义. 接触热阻的 研究主要包括基础研究和工程应用,基础研究主要 从宏观与微观方面开展理论、试验、数值模拟等方面 的预测计算,探索可用于工程实践的接触热阻半经 验关联式以及接触热阻的影响因素. 接触热阻的工 程应用主要是在接触间隙填充导热性能各异的填充 材料来实现减小或增大接触热阻的目的. 本文重点 综述宏观接触热阻的影响因素和国内外学者在基础 研究领域取得的成果,并对各种预测模型和计算方 法进行评价,指出进一步的研究重点. 对于微观接 触热阻的研究以及新型的表征方法和技术可参阅文 献[3]. 1 接触热阻国内外研究进展 接触热阻研究最早起源于 Kapitza 热阻[4] ,随后 众多学者对接触热阻进行了广泛的理论和试验研 究[5鄄鄄10] . Lambert 与 Fletcher [11鄄鄄12] , Yovanovich 与 Marotta [2]对接触热阻的预测模型和基本机理进行了 研究;Fletcher [13] , Kraus 与 Bar鄄Cohen [14]对接触热 阻进 行 了 试 验 研 究; Madhusudana 与 Fletcher [15] 、 Fletcher [16] 、Swartz 与 Pohl [6]和 Gmelin 等[17]对接触 热阻的研究进展做了综述. 国内学者对接触热阻的 研究起步较晚,90 年代南京航空航天大学顾慰兰开 发了一套接触热阻测试装置,研究了温度和接触压 力对接触热阻的影响[18鄄鄄19] ;徐烈课题组[20鄄鄄21]和王惠 龄课题组[22鄄鄄23]对低温真空环境下的接触热阻进行 了研究. 此外,应济等[24鄄鄄25]和韩玉阁等[26]也对接触 热阻进行了探索研究. 随着接触热阻研究成果的不 断涌现,近年来不少学者基于已有的研究成果,在机 械制造、航空航天、微电子等诸多领域开展了接触热 阻的应用研究,比如:基于接触热阻对笔记本电脑的 散热性能进行分析优化[27] ;接触热阻对双 H 型鳍 片管传热系数的影响[28] ;接触热阻对高速精密电主 轴热特性的影响[29] ;以及接触热阻对燃料电池温度 场分布的影响[30]等. 2 接触热阻的理论及试验研究 2郾 1 单点接触热阻的预测模型 如果把离散的局部接触简化为一个个的点,从 而开展接触热阻研究称之为点理论,接触热阻的点 理论分为单点接触理论和多点接触理论. 无论接触 面的几何形貌如何复杂,总的接触热阻总可以简化 为单点接触热阻的并联,单点接触热阻研究的核心 内容是确定接触点的形状并获取接触点的尺寸信 息,其理想的传热模型为接触点与半无限大空间或 半无限长流管之间的接触换热. 对单点接触热阻的计算通常将接触点简化为圆 柱( 圆 盘)、 半 球 体、 截 锥 体 等 模 型. Gibson [31] 、 Fenech 与Rohsenow [32] 、Cooper 等[33] 建立了圆盘计 算模型;Bahrami 等[34鄄鄄35] 建立了半球体计算模型; Willimams 与 Major [36]建 立 了 截 锥 体 计 算 模 型; Yovanovich [37]详细比较了不同形状接触点的接触换 热系数. 对于圆盘计算模型如图 1(b)所示,假定接 触点的实际半径为 a,名义接触半径为 b,名义接触 面积为 A,两接触材料的导热率分别为 姿1 和 姿2 ,则 接触点的接触热阻 R [38]为 R = A 仔a姿s f(x) (1) 1 姿s = 1 姿1 + 1 姿2 (2) f(x) = 1 - 1郾 40918x + 0郾 33801x 3 + 0郾 06792x 5 + … (3) 式中,x 为实际接触半径与名义接触半径的比值. 对于当量球模型如图 1( a)所示,Greenwood 与 Williamson [39] 、 Yovanovich [40] 、 Burde 与 Yovanov鄄 ich [41]首先把接触热阻细分为宏观接触热阻、微观 接触热阻和薄膜接触热阻. 其中薄膜接触热阻可以 通过实验来进行评估,宏观接触热阻通过当量尺寸 的单个球型接触点进行计算,假定接触面的实际接 触半径和名义接触半径分别为 a,b,两接触体的热 导率分别为 k1 ,k2 ,则交界面的宏观接触热阻为 R = g(x) 2akm (4) ·1241·
.1242· 工程科学学报.第41卷，第10期 式中，x=a/b,km为接触面的平均热导率 1=1+1 (8) 2k k2 mm m2 kmka (5) H=min(H,H) (9) g(x）=1-1.40925x+0.2959x3+0.05254x5+ 单个粗糙峰发生接触时，其接触半径α：，等效硬 度H,接触载荷F如下： 0.02105x7+… (6) 对于截锥体如图1（©）所示，假定截锥体的斜率 F:THa (10) 分别为m1,m2,两接触界面的硬度分别为H1,H,两 式中，a:=8,/m,由此可计算截锥体单热流通道的接 接触面的粗糙度为各自粗糙峰高度的均方差σ， 触热阻2-43]： σ2，同时将两相接触的粗糙表面等效为一个随机峰 1 T0.685m+(1.945+0.431m)E1. 表面和一个刚性表面相接触，则等效表面的尺寸参 Ra=4.8mk:ai (1-e2) (11) 数，诸如等效粗糙度σ、等效斜率m和等效硬度H 计算如下： 式中，￡=1-d,在上述几种计算模型中，截锥体 ma: =√+ (7) 考虑了锥角对热阻的影响，因而更加接近实际情况 (a) c 实体1 实体2 图1接触热阻计算模型.(a)半球型模型：(b)圆盘模型：(c)截椎体模型 Fig.1 Calculation model of thermal contact resistance:(a)hemisphere model;(b)disk model;(c)frustum of the cone model 2.2多点接触热阻的预测模型 触热阻.Greenwood和Williansont39)通过测量发现， 多点接触热阻计算的核心问题是确定接触点的 界面粗糙峰的高度近似服从高斯分布，建立了接触 形状尺寸和空间分布信息，从而进一步计算接触点 凸点的尺寸与密度数学模型. 的总接触面积.Yovanovich!9]与Bahrami]提出了 为了便于计算，可将两粗糙表面相接触如图3 接触表面微凸体（如图2所示）基于表观接触面积 所示，等效为一随机峰表面和一个刚性表面相接触， 随机分布的CMY(Cooper-Mikic--Yovanovich)塑性 假定等效表面的等效粗糙度为σ，粗糙峰的高度为 接触热阻模型，Mandelbrot创立分形几何以后， Z,粗糙峰的平均高度为d,名义接触面积为Am,实 Zahouani等[44]和Xu等[4把接触点密度和接触面 际接触面积为A,表面粗糙峰的密度为η，粗糙峰的 积与分形维数相结合进行了接触热阻的理论和试 斜率为m,则粗糙峰高度Z的概率密度函数p(z): 验研究.分形理论假设表面粗糙度的高度为非平 稳随机分布，以粗糙度高度的自相似和自仿射性， 推演接触点的尺度和空间分布特征[46]，但是在利 用分形维数来描述接触点密度时注意与机械加工 表面的凹坑、凸点的尺度分布相区分[).Singhal 与Garimella4s]将粗糙度曲线和粗糙蜂的变形相关 图2表面微凸体示意图 联，由单个峰的变形类型和接触热阻来计算总的接 Fig.2 Schematic of surface asperity
工程科学学报,第 41 卷,第 10 期 式中,x = a / b,km 为接触面的平均热导率. km = 2k1 k2 k1 + k2 (5) g(x) = 1 - 1郾 40925x + 0郾 2959x 3 + 0郾 05254x 5 + 0郾 02105x 7 + … (6) 对于截锥体如图 1(c)所示,假定截锥体的斜率 分别为 m1 ,m2 ,两接触界面的硬度分别为 H1 ,H2 ,两 接触面的粗糙度为各自粗糙峰高度的均方差 滓1 , 滓2 ,同时将两相接触的粗糙表面等效为一个随机峰 表面和一个刚性表面相接触,则等效表面的尺寸参 数,诸如等效粗糙度 滓、等效斜率 m 和等效硬度 H 计算如下: 滓 = 滓 2 1 + 滓 2 2 (7) 1 m = 1 m1 + 1 m2 (8) H = min(H1 ,H2 ) (9) 单个粗糙峰发生接触时,其接触半径 ai,等效硬 度 H,接触载荷 Fi 如下: Fi = 仔Ha 2 i (10) 式中,ai = 啄i / m,由此可计算截锥体单热流通道的接 触热阻[42鄄鄄43] : Rci = 1 4郾 8仔kia [ i 0郾 685m + (1郾 945 + 0郾 431m)着 ] 着 · (1 - 着 2 ) (11) 式中,着 = 1 - d mai ,在上述几种计算模型中,截锥体 考虑了锥角对热阻的影响,因而更加接近实际情况. 图 1 接触热阻计算模型. (a) 半球型模型;(b) 圆盘模型;(c) 截椎体模型 Fig. 1 Calculation model of thermal contact resistance: (a) hemisphere model; (b) disk model; (c) frustum of the cone model 2郾 2 多点接触热阻的预测模型 多点接触热阻计算的核心问题是确定接触点的 形状尺寸和空间分布信息,从而进一步计算接触点 的总接触面积. Yovanovich [9] 与 Bahrami [34] 提出了 接触表面微凸体(如图 2 所示)基于表观接触面积 随机分布的 CMY (Cooper鄄鄄 Mikic鄄鄄 Yovanovich)塑性 接触 热 阻 模 型, Mandelbrot 创 立 分 形 几 何 以 后, Zahouani 等[44]和 Xu 等[45]把接触点密度和接触面 积与分形维数相结合进行了接触热阻的理论和试 验研究. 分形理论假设表面粗糙度的高度为非平 稳随机分布,以粗糙度高度的自相似和自仿射性, 推演接触点的尺度和空间分布特征[46] ,但是在利 用分形维数来描述接触点密度时注意与机械加工 表面的凹坑、凸点的尺度分布相区分[47] . Singhal 与Garimella [48]将粗糙度曲线和粗糙峰的变形相关 联,由单个峰的变形类型和接触热阻来计算总的接 触热阻. Greenwood 和 Willianson [39] 通过测量发现, 界面粗糙峰的高度近似服从高斯分布,建立了接触 凸点的尺寸与密度数学模型. 图 2 表面微凸体示意图 Fig. 2 Schematic of surface asperity 为了便于计算,可将两粗糙表面相接触如图 3 所示,等效为一随机峰表面和一个刚性表面相接触, 假定等效表面的等效粗糙度为 滓,粗糙峰的高度为 Z,粗糙峰的平均高度为 d,名义接触面积为 Am ,实 际接触面积为 Ac,表面粗糙峰的密度为 浊,粗糙峰的 斜率为 m,则粗糙峰高度 Z 的概率密度函数 渍(z): ·1242·
汪献伟等：宏观接触热阻研究综述 ·1243· p(a)=1e总 (12) 触针式轮廓测量仪操作简单，测量范围大，但是测量 √2T 精度低并且会对接触面产生一定程度上的损伤.采 任意粗糙峰高度Z与刚性表面发生接触的概 用光学和激光扫描测量法不破坏表面形貌，测量精 率P: 度高，但是受精度限制，测量范围较小[3] P.={z>d=9(a)d (13) 2.4接触热阻的数值模拟 尽管众多学者对接触热阻进行了理论研究并建 理论接触面积上发生接触的粗糙峰的总数N。: 立了一些理论模型，但是一些参数是根据实验结果 N.=n19(e)d (14) 确定的，在实际的工程应用中具有一定的难度.此 实际总接触面积A。: 外，接触热阻的试验研究需要基于不同的试验条件 4.T0A 进行大量的数据测量，测量后的数据需要进一步的 2 丁d (z-d)2p(z）dz (15) 2m2 分析处理，往往费时费力.而采用数值模拟只需要 在此基础上基于单点接触热阻理论，即可确定 建立接触模型并施加相应的热边界条件即可完成接 接触界面的总接触热阻R,: 触热阻的计算，整个模拟过程简单、高效，模拟结果 可用于同理论和试验数据对比，具有一定的理论意 1 RT= (16) 义和工程应用价值 使用仿真软件开展接触热阻的模拟计算主要包 人界面1的微凸体 含两个方面的内容：在接触热阻未知的情况下，通过 在接触界面创建微凸体等效模型来模拟界面的微凸 界面2的微凸体 体，并基于一定的边界条件通过数值计算获得边界 的温差△1和热流通量q,据此，由公式(17)计算界 777777777777777777777777777777777777777 图3界面微凸体接触形式 面的接触热阻，但该计算结果的准确性受有限元模 Fig.3 Contact form of interface asperity 型中微凸体的尺寸参数和密度分布的影响. 2.3接触热阻的试验研究 R=△ (17) 9 对于一维传热，接触热阻的测试可基于一维稳 基于上述理念，陆敏恂等[5]使用长方体微凸体 态传热模型如图4所示，在测试样品的周围布置绝 来模拟MOSFET功率管和散热片之间的接触峰，并 缘材料，假定热量从样品轴向流过，测量两接触体的 计算了界面的接触热阻：刘冬欢等]利用呈周期性 温度，由傅里叶定律计算两接触界面的接触热阻. 分布的矩形微凸体模拟航空航天设备热防护结构的 浙江大学应济等据此设计了一套接触热阻测量装 接触热阻.余远峰等6]等通过反演得到界面接触 置，开展了试验测量并与理论计算结果进行了对 热阻并与实际的接触热阻进行了对比.在接触热阻 比25】，上海交通大学的徐烈等9将实验设备置于 已知的情况下，可以通过有限元分析软件中的接触 液氨中测试了低温下的固体接触热阻.该测量方法 单元来模拟接触热阻对温度场分布的影响，该方法 虽然简单便于操作，但是测量需要的时间较长，并且 相对于接触区域过渡薄层法更加简便，解决了接触 当热电偶嵌入测量试件时会影响温度场分布和测量 区域过渡薄层网格划分困难、计算量巨大的难题. 的准确性[o] 杨华威等s7-s]使用ANSYS的接触单元模拟了接触 为了避免接触式测量引起的温度误差积累，可 热阻对空间行波管收集极温度场部分的影响：曹骏 以采用非接触式的超声1)、激光[s2]等测量方法.激 等[9]模拟了接触热阻对机床主轴系统温度场分布 光光热测量法在测量材料的热物性时，利用了材料 的影响 的光学反射率随温度而变化的特性.该方法具有非 3接触热阻的影响因素 接触、响应快、测量周期短等优点，在许多领域得到 了广泛的应用.除了上述通过测量温度来间接计算 接触热阻的变化受接触面的几何形貌、接触压 接触热阻，也可以通过测量接触面的三维几何形貌， 力、接触变形特性、温度梯度、材料的热物性、接触界 获得接触点的尺度、密度和接触面积信息从而计算 面的填充情况、滑动状态、热流方向等因素的影响， 接触热阻.对于样品几何形貌的测量可以分为接触 属于几何学，材料学、力学、传热学等多学科交叉问 式的触针测量和非接触式的光学、激光扫描测量. 题.湛利华等[60]，顾慰兰等18-19通过试验的方法定
汪献伟等: 宏观接触热阻研究综述 渍(z) = 1 2仔 e - z2 2滓2 (12) 任意粗糙峰高度 Z 与刚性表面发生接触的概 率 Pr: Pr = {z > d} = 乙 肄 d 渍(z)dz (13) 理论接触面积上发生接触的粗糙峰的总数 Nc: Nc = 浊Am 乙 肄 d 渍(z)dz (14) 实际总接触面积 Ac: Ac = 仔浊Am 2m 2 乙 肄 d (z - d) 2 渍(z)dz (15) 在此基础上基于单点接触热阻理论,即可确定 接触界面的总接触热阻 RT : RT = 移 n i = 1 1 Ri (16) 图 3 界面微凸体接触形式 Fig. 3 Contact form of interface asperity 2郾 3 接触热阻的试验研究 对于一维传热,接触热阻的测试可基于一维稳 态传热模型如图 4 所示,在测试样品的周围布置绝 缘材料,假定热量从样品轴向流过,测量两接触体的 温度,由傅里叶定律计算两接触界面的接触热阻. 浙江大学应济等据此设计了一套接触热阻测量装 置,开展了试验测量并与理论计算结果进行了对 比[25] ,上海交通大学的徐烈等[49] 将实验设备置于 液氮中测试了低温下的固体接触热阻. 该测量方法 虽然简单便于操作,但是测量需要的时间较长,并且 当热电偶嵌入测量试件时会影响温度场分布和测量 的准确性[50] . 为了避免接触式测量引起的温度误差积累,可 以采用非接触式的超声[51] 、激光[52]等测量方法. 激 光光热测量法在测量材料的热物性时,利用了材料 的光学反射率随温度而变化的特性. 该方法具有非 接触、响应快、测量周期短等优点,在许多领域得到 了广泛的应用. 除了上述通过测量温度来间接计算 接触热阻,也可以通过测量接触面的三维几何形貌, 获得接触点的尺度、密度和接触面积信息从而计算 接触热阻. 对于样品几何形貌的测量可以分为接触 式的触针测量和非接触式的光学、激光扫描测量. 触针式轮廓测量仪操作简单,测量范围大,但是测量 精度低并且会对接触面产生一定程度上的损伤. 采 用光学和激光扫描测量法不破坏表面形貌,测量精 度高,但是受精度限制,测量范围较小[53] . 2郾 4 接触热阻的数值模拟 尽管众多学者对接触热阻进行了理论研究并建 立了一些理论模型,但是一些参数是根据实验结果 确定的,在实际的工程应用中具有一定的难度. 此 外,接触热阻的试验研究需要基于不同的试验条件 进行大量的数据测量,测量后的数据需要进一步的 分析处理,往往费时费力. 而采用数值模拟只需要 建立接触模型并施加相应的热边界条件即可完成接 触热阻的计算,整个模拟过程简单、高效,模拟结果 可用于同理论和试验数据对比,具有一定的理论意 义和工程应用价值. 使用仿真软件开展接触热阻的模拟计算主要包 含两个方面的内容:在接触热阻未知的情况下,通过 在接触界面创建微凸体等效模型来模拟界面的微凸 体,并基于一定的边界条件通过数值计算获得边界 的温差 驻t 和热流通量 q,据此,由公式(17)计算界 面的接触热阻,但该计算结果的准确性受有限元模 型中微凸体的尺寸参数和密度分布的影响. R = 驻t q (17) 基于上述理念,陆敏恂等[54]使用长方体微凸体 来模拟 MOSFET 功率管和散热片之间的接触峰,并 计算了界面的接触热阻;刘冬欢等[55] 利用呈周期性 分布的矩形微凸体模拟航空航天设备热防护结构的 接触热阻. 余远峰等[56] 等通过反演得到界面接触 热阻并与实际的接触热阻进行了对比. 在接触热阻 已知的情况下,可以通过有限元分析软件中的接触 单元来模拟接触热阻对温度场分布的影响,该方法 相对于接触区域过渡薄层法更加简便,解决了接触 区域过渡薄层网格划分困难、计算量巨大的难题. 杨华威等[57鄄鄄58]使用 ANSYS 的接触单元模拟了接触 热阻对空间行波管收集极温度场部分的影响;曹骏 等[59]模拟了接触热阻对机床主轴系统温度场分布 的影响. 3 接触热阻的影响因素 接触热阻的变化受接触面的几何形貌、接触压 力、接触变形特性、温度梯度、材料的热物性、接触界 面的填充情况、滑动状态、热流方向等因素的影响, 属于几何学,材料学、力学、传热学等多学科交叉问 题. 湛利华等[60] ,顾慰兰等[18鄄鄄19]通过试验的方法定 ·1243·
.1244. 工程科学学报.第41卷，第10期 ∠E之∠∠之之∠之之∠∠2之∠∠之∠∠之∠∠之之∠之之∠∠之之∠之2∠∠ 接触件温度K 热流从压缩 实体1 实体2 AAAAAAA 接触件 盘流向缓种 10.5556 777777777777777777777777777777777777777 16.1111 21.6667 27.2222 微凸体温度K 32.7778 38.3333 43.8889 子 494444 44.6145 55 45.6637 46.7129 微凸体 47.7621 48.8113 位置 49.8606 50.9098 图4交界面温度变化示意图 51.9590 53.0082 Fig.4 Schematic of temperature variation on the interface 54.0574 性研究了各因素对接触热阻的影响： 图5热流从316LN流向G10时的温度分布 1)表征表面几何形状的粗糙度对接触热阻的 Fig.5 Temperature distribution during heat transfer from 316LN to 影响.粗糙度越大接触热阻越大，这是由于表面粗 G10 糙度的大小直接影响两固体表面的接触状态，当粗 糙度增大时，热流线收缩加剧，影响热流的传递 接触件 接触件温度K 热流从压缩 效率 盘流向缓冲 10.5556 2)为了研究热流方向对接触热阻的影响，中南 16.1111 216667 大学湛利华等人以钢-铝接触为研究对象，分析了 27.2222 32.7778 热流方向对接触热阻的影响，当热流从导热率较高 38.3333 微凸体 43.8889 温度K 的钢流向导热率较低的铝时接触热阻减小，反之，接 49.4444 12.5435 55 14.4821 触热阻增大，但是该分析结果与Rogers等的研究结 16.4208 微凸体 18.3594 果正好相反 2029R0 22.2367 为了进一步验证热流方向对接触热阻的影响， 24.1753 26.1140 本项目组针对中国聚变工程实验堆(Chinese fusion 28.0526 29.9912 engineering test reactor,CFETR)中心螺线管模型线 圈降温以及磁体预紧支撑部件漏热过程，采用多点 图6热流从G10流向316LN时的温度分布 接触热阻的圆盘等效模型计算热流方向改变时压缩 Fig.6 Temperature distribution during heat transfer from G10 to 316LN 盘不锈钢316LN和缓冲区高性能绝缘材料G10之 间的接触热阻.热分析的边界条件采用降温过程所 一定减小，接触热阻的变化需要和具体的材料属性 允许的极限温差50K,微凸体的接触热阻采用等效 以及相应的工作温度相结合.当材料的工作温度区 接触热阻，微凸体的粗糙度取两接触面粗糙度均值 间较大且热导率随温度变化明显时，接触热阻的计 38m,接触间隙比0.1.当热流从高热导率的压缩 算需要考虑热导率随温度的变化状态.对于超导磁 盘流向缓冲区时接触热阻为1.85×10-3m2.K· 体线圈常用的铠装式馆内线缆导体(CCC),励磁前 W-1,图5为相应的温度场分布.当热流从较低热 需要从室温降至4.5K的超低温，图7给出了接触 导率的缓冲区流向较高热导率的压缩盘时接触热阻 体之间最大温差△：=2K时，铠装式馆内线缆导体 为3.14×10-3m2.K·W-1,图6为相应的温度场分 各组件的接触热阻随温度的变化情况.显然温度升 布，该分析结果与湛利华等的试验结果一致. 高时316LN和G10的接触热阻减小，然而对于超导 3)温度对接触热阻的影响表现为温度升高以 电缆中的铜芯，由于其热导率随温度的升高而降低， 后接触热阻减小〔61-62】，研究表明温度并不是导致接 当温度升高以后铜导体C10200的接触热阻不仅没 触热阻减小的直接因素，温度的变化使接触表面发 有减小，反而出现了些许增大，这与铜导体在超低温 生蠕变，从而间接导致接触热阻发生变化.该结论 下对超导线圈的失超保护功能相吻合.为了进一步 基于特定的实验条件，并不具有普遍的适用性 验证理论计算的可靠性，对铜-石墨在100~400℃ 本文研究表明，温度升高材料的接触热阻并不 的界面接触热阻进行试验研究（此温度区间内，铜-
工程科学学报,第 41 卷,第 10 期 图 4 交界面温度变化示意图 Fig. 4 Schematic of temperature variation on the interface 性研究了各因素对接触热阻的影响: 1)表征表面几何形状的粗糙度对接触热阻的 影响. 粗糙度越大接触热阻越大,这是由于表面粗 糙度的大小直接影响两固体表面的接触状态,当粗 糙度增大时,热流线收缩加剧,影响热流的传递 效率. 2)为了研究热流方向对接触热阻的影响,中南 大学湛利华等人以钢鄄鄄 铝接触为研究对象,分析了 热流方向对接触热阻的影响,当热流从导热率较高 的钢流向导热率较低的铝时接触热阻减小,反之,接 触热阻增大,但是该分析结果与 Rogers 等的研究结 果正好相反. 为了进一步验证热流方向对接触热阻的影响, 本项目组针对中国聚变工程实验堆(Chinese fusion engineering test reactor, CFETR)中心螺线管模型线 圈降温以及磁体预紧支撑部件漏热过程,采用多点 接触热阻的圆盘等效模型计算热流方向改变时压缩 盘不锈钢 316LN 和缓冲区高性能绝缘材料 G10 之 间的接触热阻. 热分析的边界条件采用降温过程所 允许的极限温差 50 K,微凸体的接触热阻采用等效 接触热阻,微凸体的粗糙度取两接触面粗糙度均值 38 滋m,接触间隙比 0郾 1. 当热流从高热导率的压缩 盘流向缓冲区时接触热阻为 1郾 85 伊 10 - 3 m 2·K· W - 1 ,图 5 为相应的温度场分布. 当热流从较低热 导率的缓冲区流向较高热导率的压缩盘时接触热阻 为 3郾 14 伊 10 - 3 m 2·K·W - 1 ,图 6 为相应的温度场分 布,该分析结果与湛利华等的试验结果一致. 3)温度对接触热阻的影响表现为温度升高以 后接触热阻减小[61鄄鄄62] ,研究表明温度并不是导致接 触热阻减小的直接因素,温度的变化使接触表面发 生蠕变,从而间接导致接触热阻发生变化. 该结论 基于特定的实验条件,并不具有普遍的适用性. 本文研究表明,温度升高材料的接触热阻并不 图 5 热流从 316LN 流向 G10 时的温度分布 Fig. 5 Temperature distribution during heat transfer from 316LN to G10 图 6 热流从 G10 流向 316LN 时的温度分布 Fig. 6 Temperature distribution during heat transfer from G10 to 316LN 一定减小,接触热阻的变化需要和具体的材料属性 以及相应的工作温度相结合. 当材料的工作温度区 间较大且热导率随温度变化明显时,接触热阻的计 算需要考虑热导率随温度的变化状态. 对于超导磁 体线圈常用的铠装式馆内线缆导体(CICC),励磁前 需要从室温降至 4郾 5 K 的超低温,图 7 给出了接触 体之间最大温差 驻t = 2 K 时,铠装式馆内线缆导体 各组件的接触热阻随温度的变化情况. 显然温度升 高时 316LN 和 G10 的接触热阻减小,然而对于超导 电缆中的铜芯,由于其热导率随温度的升高而降低, 当温度升高以后铜导体 C 10200 的接触热阻不仅没 有减小,反而出现了些许增大,这与铜导体在超低温 下对超导线圈的失超保护功能相吻合. 为了进一步 验证理论计算的可靠性,对铜鄄鄄 石墨在 100 ~ 400 益 的界面接触热阻进行试验研究(此温度区间内,铜鄄鄄 ·1244·
汪献伟等：宏观接触热阻研究综述 ·1245· 石墨的热导率均随温度的升高而降低)，图8的实 20 -0.6 MPa 验数据表明温度升高，接触热导减小 -0*1.4MPa 18 ---1.8MPa 1.72 716 --…2.1Pa ---2.4MPa 1.70 2.8 MPa 16 7.5 --◆-·6.1Pa 316L.N接触热阴1.66 14 --C10200接角触热阻1.64 3 14 6.5 。G10接触热阻 1.62 12 1.60 00 6.0 1.58 12 1.56 5.5 10 1.54 1.52 9 80100120140160180200220240260 1953333934343453473 493 界面平均温度K 界面平均温度K 图7不同材料随温度变化的接触热阻 图8在真空度为1.5×102Pa和特定接触压力下接触热导随 Fig.7 Changes of thermal contact resistance with the temperature of 温度变化 different materials Fig.8 Thermal contact conductance variation versus temperature un- der the vacuum of 1.5x10-2 Pa and specific contact pressure 4)接触界面在外部接触载荷的作用下，接触点 的可能变形形式为弹性、塑性和弹塑性，不同的变形 700 ·T=418K,T=361K 形式对应不同的接触热阻.对于均匀加载，在加载 600oT=513K.T=401K 过程中随着加载压力的增大，接触热阻呈减小趋 T=693K,T=473K 500 试验测量 —T-418K,T-361K 势[3-].接触压力的增大，实际上增加了接触界面 400 --T,=513K,T=401K 的有效接触面积，有助于热量的传递.对于弹性变 --T=693K,T-473K 理论计算 形，卸载后接触热阻无明显变化，对于塑性变形卸载 此 里 后接触热阻减小.Cooper等[)、Johson]、Sridgar -Y 与Yovanovich[]基于半球型接触峰模型给出了在 不同变形下，接触压力和接触面积之间的关系 4 6 接触压力/MPa H. 塑性变形 图9真空度为6Pa时接触热导与接触压力关系 A P Fig.9 Relationship between thermal contact conductance and contact 弹塑性变形 A ee (18) pressure at a vacuum of 6 Pa T 弹性变形 流：真空接触间隙之间的热辐射，对接触表面进行表 K 面处理可以改变导热状态，例如在界面间隙填充导 式中，A,为实际接触面积，A。为名义接触面积；P.为 热膏、导热垫、导热凝胶以及具有良好适应性和导热 外部压力载荷，H。为塑性变形接触面的平均屈服压 性的液态介质：对接触表面进行焊接处理：在金属基 力；H,为弹塑性变形的屈服压力；σ，为表面微凸体 体表面涂覆铂、铟、镍、石墨等材料都可以有效的减 高度的标准差，K为表面的均方根斜率.本项目组 小接触热阻，川.此外，真空度也是影响接触热阻 对铜-石墨开展了真空环境下接触热阻测试分析， 的因素之一，真空度的变化不仅影响接触间隙的热 图9的实验测量和理论计算均表明随着接触压力的 对流，而且导致空气和接触材料热导率的变化，关于 增大接触热导近似呈指数增长. 不同真空度下材料界面接触热阻的变化状态，则需 5)螺栓连接的接触热阻属于在非均匀压力作 要进一步的深入研究 用下的传导问题，Roca与Mikic[68-9]给出了基于高 4 结论 斯分布的粗糙平面螺栓连接的比较完善的接触热导 计算理论，但是该理论不便于工程设计人员分析使 迄今为止，国内外学者从理论到试验对接触热 用.Bevans等[0对螺栓连接模型进行了简化，提出 阻开展了大量研究，部分学者将已有的研究成果应 了一种更简便的计算方法 用于工程实践有效解决了一些热力学问题.为了更 6)接触界面的热量传递方式包括接触面凸峰 加准确的预测接触热阻，下述工作需要深入开展： 之间的热传导，接触间隙介质之间的热传导和热对 (1)当前关于接触热阻的半经验公式均是建立
汪献伟等: 宏观接触热阻研究综述 石墨的热导率均随温度的升高而降低),图 8 的实 验数据表明温度升高,接触热导减小. 图 7 不同材料随温度变化的接触热阻 Fig. 7 Changes of thermal contact resistance with the temperature of different materials 4)接触界面在外部接触载荷的作用下,接触点 的可能变形形式为弹性、塑性和弹塑性,不同的变形 形式对应不同的接触热阻. 对于均匀加载,在加载 过程中随着加载压力的增大,接触热阻呈减小趋 势[63鄄鄄65] . 接触压力的增大,实际上增加了接触界面 的有效接触面积,有助于热量的传递. 对于弹性变 形,卸载后接触热阻无明显变化,对于塑性变形卸载 后接触热阻减小. Cooper 等[33] 、 Johson [66] 、 Sridgar 与 Yovanovich [67]基于半球型接触峰模型给出了在 不同变形下,接触压力和接触面积之间的关系. Ar Aa = Pa Hp 塑性变形 Pa Hep 弹塑性变形 仔 滓sKs ì î í ï ï ï ï ï ï ï ï 弹性变形 (18) 式中,Ar 为实际接触面积,Aa 为名义接触面积;Pa 为 外部压力载荷,Hp 为塑性变形接触面的平均屈服压 力;Hep为弹塑性变形的屈服压力;滓s 为表面微凸体 高度的标准差,Ks 为表面的均方根斜率. 本项目组 对铜鄄鄄石墨开展了真空环境下接触热阻测试分析, 图 9 的实验测量和理论计算均表明随着接触压力的 增大接触热导近似呈指数增长. 5)螺栓连接的接触热阻属于在非均匀压力作 用下的传导问题,Roca 与 Mikic [68鄄鄄69] 给出了基于高 斯分布的粗糙平面螺栓连接的比较完善的接触热导 计算理论,但是该理论不便于工程设计人员分析使 用. Bevans 等[70]对螺栓连接模型进行了简化,提出 了一种更简便的计算方法. 6)接触界面的热量传递方式包括接触面凸峰 之间的热传导,接触间隙介质之间的热传导和热对 图 8 在真空度为 1郾 5 伊 10 - 2 Pa 和特定接触压力下接触热导随 温度变化 Fig. 8 Thermal contact conductance variation versus temperature un鄄 der the vacuum of 1郾 5 伊 10 - 2 Pa and specific contact pressure 图 9 真空度为 6 Pa 时接触热导与接触压力关系 Fig. 9 Relationship between thermal contact conductance and contact pressure at a vacuum of 6 Pa 流;真空接触间隙之间的热辐射,对接触表面进行表 面处理可以改变导热状态,例如在界面间隙填充导 热膏、导热垫、导热凝胶以及具有良好适应性和导热 性的液态介质;对接触表面进行焊接处理;在金属基 体表面涂覆铂、铟、镍、石墨等材料都可以有效的减 小接触热阻[9,71] . 此外,真空度也是影响接触热阻 的因素之一,真空度的变化不仅影响接触间隙的热 对流,而且导致空气和接触材料热导率的变化,关于 不同真空度下材料界面接触热阻的变化状态,则需 要进一步的深入研究. 4 结论 迄今为止,国内外学者从理论到试验对接触热 阻开展了大量研究,部分学者将已有的研究成果应 用于工程实践有效解决了一些热力学问题. 为了更 加准确的预测接触热阻,下述工作需要深入开展: (1)当前关于接触热阻的半经验公式均是建立 ·1245·
·1246· 工程科学学报.第41卷，第10期 在一定的理论假设和试验测量之上，在工程应用上 [11]Lambert MA,Fletcher L S.Review of models for thermal contact 具有一定的局限性.同时，试验测量系统的精确性 conductance of metals.J Thermophys Heat Transfer,1997,11 和试验方法的科学性直接影响接触热阻理论计算的 (2):129 12]Lambert M A,Fletcher L S.Thermal contact conductance of non- 准确性，如何更为高效、准确的测量接触表面的形貌 flat,rough,metallic coated metals.J Heat Transfer,2002,124 特征，能否建立一套具有普适性的接触热阻理论是 (3):405 未来的努力方向. [13] Fletcher L.S.A review of thermal enhancement techniques for (2)数值模拟法计算接触热阻具有快捷、高效 electronic systems.IEEE Trans Comp Hybrids Mfg,1990,13 等优点，但是仿真模型往往人为假定接触峰形状和 (4):1012 [14]Kraus A D,Bar-Cohen A.Thermal Analysis Control of Elec- 接触间隙比，实际上由于接触表面形貌的复杂性和 tronic Equipment.Florida:MeGraw-Hill,1983 接触点分布的随机性，为了更加准确的获知接触点 [15]Madhusudana C V,Fletcher L S.Contact heat transfer-the last 大小和空间分布信息，需要开展更为精确的三维表 decade.AlAA J,1986,24(3):510 面测量和表征. [16]Fletcher L S.Recent developments in contact conductance heat (3)接触热阻的各影响因素往往是相互耦合 transfer.J Heat Transfer,1988,110(4B):1059 的，基于多场耦合理论借助大型仿真软件采用直接 [17]Gmelin E,Asen-Palmer M,Reuther M,et al.Thermal boundary resistance of mechanical contacts between solids at sub-ambient 耦合法开展多耦合分析有助于更加的准确预测接触 temperatures.J Phys D Appl Phys,1999,32(6):R19 热阻. [18]Gu W L.Experiment on thermal contact resistance.J Nanjing (4)接触热阻是一个多学科交叉的跨尺度问 4 eron Inst,1992,24(1):46 题，在微尺度领域傅里叶定律不再适用，为了准确预 (顾慰兰.接触热阻的试验研究.南京航空学院学报，1992， 测接触热阻需要宏观和微观相结合，进一步深入开 24(1):46) 展跨尺度研究 [19]Gu WL.Yang Y S.The influence of temperature on contact re- sistance of metallic surfaces temperature.J Nanjing Unir Aeron 参考文献 4sron,1994,26(3):342 (顾慰兰，杨燕生.温度对接触热阻的影响.南京航空航天 [1]Lambert M A,Fletcher L S.Thermal contact conductance of 大学学报，1994,26(3)：342) spherical rough metals.J Heat Transfer,1997,119(4):684 [20]Xu L,Yang J,Xu J M,et al.The research of contact thermal [2]Adrian Bejan,Allan D K.Heat Transfer Handbook.New Jersey: resistance of solid surface at low temperature.Cryog Supercond, Wiley.2003 1996,24(1):53 [3]Song C Y,Luan T,Wu J B,et al.Research progress on the inter- (徐烈，杨军，徐佳梅，等.低温下固体表面接触热阻的研 face thermal resistance of polymer based thermal interface materi- 究.低温与超导，1996,24(1)：53) als.J Integr Technol,2019.8(1):56 [21]Xu L,Zhang T,Xiong W,et al.The effect of contact interfaces (宋成轶，栾添，邬剑波，等.聚合物基热界面材料界面接触 on thermal contact resistance at low temperature and vacuum.Vac 热阻的研究进展.集成技术，2019,8(1)：56) Cog,1998,4(1)1) [4]Dter Haar.Collected Papers of P.L.Kapitza.London:Pergamon (徐烈，张涛，熊炜，等.真空低温下接触表面对接触热阻的 Press,1965 影响.真空与低温，1998,4(1)：1) [5]Thomas T R,Probert S D.Correlations for thermal contact con- [22]Xu R P,Feng H D,Zhao L P,et al.Experimental investigation ductance in vacuo.J Heat Transfer,1972.94(3):276 of thermal contact conductance at low temperature based on fractal [6] Swartz E T,Pohl R 0.Thermal boundary resistance.Rer Mod description.Int Commun Heat Mass Transfer,2006,33(7):811 Phs,1989,61(3):605 [23]Chen J.Inrestigation on the Computer Simulation of Boundary [7]Hmina N,Scudeller Y.Thermal interface resistance and subsur- Heat Transfer Process betwceen Solids at Low Temperature [Disser- face effusivity of submicron metallic films on dielectric substrates: tation ]Wuhan:Huazhong University of Science and Technolo- an experimental method for simultaneous determination.Int Heat ,2004 Mass Transfer,1998,41(18):2781 (陈进.低温固体界面传热过程计算机仿真研究[学位论 [8]Rostami AA,Hassan A Y,Lim PC.Parametric study of thermal 文].武汉：华中科技大学，2004) constriction resistance.Heat Mass Transfer,2001,37(1):5 [24]Ying J.Chen WQ.Chen ZC.Effect of heat transfer on the esti- [9]Yovanovich MM.Four decades of research on thermal contact, mation of contact resistance.JZhejiang Unie Nat Sci,1996,30 gap,and joint resistance in microelectronies.IEEE Trans Compon (4):415) Packag Technol,2005,28(2):182 (应济，陈伟球，陈子辰.换热影响下的接触热阻估计.浙江 [10]Prasher R.Thermal boundary resistance and thermal conductivity 大学学报：自然科学版，1996,30(4)：415) of multiwalled carbon nanotubes.Phys Rev B,2008.77(7): [25]Ying J.Jia Y.Chen Z C.et al.Theoretical and experimental re- 075424 search on the contact thermal resistance between real surface
工程科学学报,第 41 卷,第 10 期 在一定的理论假设和试验测量之上,在工程应用上 具有一定的局限性. 同时,试验测量系统的精确性 和试验方法的科学性直接影响接触热阻理论计算的 准确性,如何更为高效、准确的测量接触表面的形貌 特征,能否建立一套具有普适性的接触热阻理论是 未来的努力方向. (2)数值模拟法计算接触热阻具有快捷、高效 等优点,但是仿真模型往往人为假定接触峰形状和 接触间隙比,实际上由于接触表面形貌的复杂性和 接触点分布的随机性,为了更加准确的获知接触点 大小和空间分布信息,需要开展更为精确的三维表 面测量和表征. (3)接触热阻的各影响因素往往是相互耦合 的,基于多场耦合理论借助大型仿真软件采用直接 耦合法开展多耦合分析有助于更加的准确预测接触 热阻. (4)接触热阻是一个多学科交叉的跨尺度问 题,在微尺度领域傅里叶定律不再适用,为了准确预 测接触热阻需要宏观和微观相结合,进一步深入开 展跨尺度研究. 参 考 文 献 [1] Lambert M A, Fletcher L S. Thermal contact conductance of spherical rough metals. J Heat Transfer, 1997, 119(4): 684 [2] Adrian Bejan, Allan D K. Heat Transfer Handbook. New Jersey: Wiley, 2003 [3] Song C Y, Luan T, Wu J B, et al. Research progress on the inter鄄 face thermal resistance of polymer based thermal interface materi鄄 als. J Integr Technol, 2019, 8(1): 56 (宋成轶, 栾添, 邬剑波, 等. 聚合物基热界面材料界面接触 热阻的研究进展. 集成技术, 2019, 8(1): 56) [4] Dter Haar. Collected Papers of P. L. Kapitza. London: Pergamon Press, 1965 [5] Thomas T R, Probert S D. Correlations for thermal contact con鄄 ductance in vacuo. J Heat Transfer, 1972, 94(3): 276 [6] Swartz E T, Pohl R O. Thermal boundary resistance. Rev Mod Phys, 1989, 61(3): 605 [7] Hmina N, Scudeller Y. Thermal interface resistance and subsur鄄 face effusivity of submicron metallic films on dielectric substrates: an experimental method for simultaneous determination. Int J Heat Mass Transfer, 1998, 41(18): 2781 [8] Rostami A A, Hassan A Y, Lim P C. Parametric study of thermal constriction resistance. Heat Mass Transfer, 2001, 37(1): 5 [9] Yovanovich M M. Four decades of research on thermal contact, gap, and joint resistance in microelectronics. IEEE Trans Compon Packag Technol, 2005, 28(2): 182 [10] Prasher R. Thermal boundary resistance and thermal conductivity of multiwalled carbon nanotubes. Phys Rev B, 2008, 77 ( 7 ): 075424 [11] Lambert M A, Fletcher L S. Review of models for thermal contact conductance of metals. J Thermophys Heat Transfer, 1997, 11 (2): 129 [12] Lambert M A, Fletcher L S. Thermal contact conductance of non鄄 flat, rough, metallic coated metals. J Heat Transfer, 2002, 124 (3): 405 [13] Fletcher L S. A review of thermal enhancement techniques for electronic systems. IEEE Trans Comp Hybrids Mfg, 1990, 13 (4): 1012 [14] Kraus A D, Bar鄄Cohen A. Thermal Analysis & Control of Elec鄄 tronic Equipment. Florida: McGraw鄄Hill, 1983 [15] Madhusudana C V, Fletcher L S. Contact heat transfer鄄the last decade. AIAA J, 1986, 24(3): 510 [16] Fletcher L S. Recent developments in contact conductance heat transfer. J Heat Transfer, 1988, 110(4B): 1059 [17] Gmelin E, Asen鄄Palmer M, Reuther M, et al. Thermal boundary resistance of mechanical contacts between solids at sub鄄ambient temperatures. J Phys D Appl Phys, 1999, 32(6): R19 [18] Gu W L. Experiment on thermal contact resistance. J Nanjing Aeron Inst, 1992, 24(1): 46 (顾慰兰. 接触热阻的试验研究. 南京航空学院学报, 1992, 24(1): 46) [19] Gu W L, Yang Y S. The influence of temperature on contact re鄄 sistance of metallic surfaces temperature. J Nanjing Univ Aeron Astron, 1994, 26(3): 342 (顾慰兰, 杨燕生. 温度对接触热阻的影响. 南京航空航天 大学学报, 1994, 26(3): 342) [20] Xu L, Yang J, Xu J M, et al. The research of contact thermal resistance of solid surface at low temperature. Cryog Supercond, 1996, 24(1): 53 (徐烈, 杨军, 徐佳梅, 等. 低温下固体表面接触热阻的研 究. 低温与超导, 1996, 24(1): 53) [21] Xu L, Zhang T, Xiong W, et al. The effect of contact interfaces on thermal contact resistance at low temperature and vacuum. Vac Cryog, 1998, 4(1): 1) (徐烈, 张涛, 熊炜, 等. 真空低温下接触表面对接触热阻的 影响. 真空与低温, 1998, 4(1): 1) [22] Xu R P, Feng H D, Zhao L P, et al. Experimental investigation of thermal contact conductance at low temperature based on fractal description. Int Commun Heat Mass Transfer, 2006, 33(7): 811 [23] Chen J. Investigation on the Computer Simulation of Boundary Heat Transfer Process between Solids at Low Temperature [Disser鄄 tation]. Wuhan: Huazhong University of Science and Technolo鄄 gy, 2004 (陈进. 低温固体界面传热过程计算机仿真研究[ 学位论 文]. 武汉: 华中科技大学, 2004) [24] Ying J, Chen W Q, Chen Z C. Effect of heat transfer on the esti鄄 mation of contact resistance. J Zhejiang Univ Nat Sci, 1996, 30 (4): 415) (应济, 陈伟球, 陈子辰. 换热影响下的接触热阻估计. 浙江 大学学报:自然科学版, 1996, 30(4): 415) [25] Ying J, Jia Y, Chen Z C, et al. Theoretical and experimental re鄄 search on the contact thermal resistance between real surface. J ·1246·
汪献伟等：宏观接触热阻研究综述 ·1247. Zhejiang Unie Nat Sci,1997,31(1):104 [41]Burde SS,Yovanovich MM.Thermal resistance at smooth- (应济，贾显，陈子辰，等.粗糙表面接触热阻的理论和实验 sphere/rough-flat contacts:Theoretical analysis /2nd AlAA/ 研究.浙江大学学报：自然科学版，1997,31(1)：104) ASME Thermophysics and Heat Transfer Conference.Palo Alto. [26]Han Y G,Xuan Y M,Tang T F.Heat transfer on the contact 1978:83 surfaces with frictional heating.J Nanjing Unig Sci Technol, [42]Shai I,Santo M.Heat transfer with contact resistance.Int Heat 1998,22(3):260 Mass Transfer,1982,25(4):465 (韩玉阁，宜益民，汤瑞蜂.摩擦接触界面传热规律研究.南 [43]Williams A.Heat flow through single point of metallic contacts of 京理T大学学报，1998,22(3)：260) simple shapes.AlAA/ASME Thermo-physics and Heat Transfer [27]Lin PZ,Sun B B,Zhao T.Thermal analysis and optimization of Conference.Boston,1974:1 notebook computer based on thermal contact resistance.South- [44]Zahouani H,Vargiolu R,Loubet J L.Fractal models of surface east Unin Nat Sci Ed,2014,44(1):93 topography and contact mechanics.Math Comput Modell,1998, (林潘忠，孙蓓蓓，赵天.基于接触热阻模型的笔记本电脑 28(4-8):517 散热性能分析与优化.东南大学学报：自然科学版，2014， [45]Xu R P,Feng H D,Zhao L P,et al.Experimental investigation 44(1):93) of thermal contact conductance at low temperature based on fractal [28]Ge M,Zhao L J,Shu S X,et al.Impact of thermal contact re- description.Int Commun Heat Mass Transfer,2006,33(7):811 sistance on heat transfer coefficient of double H-type finned [46]Wang A L,Zhao J F.Review of prediction for thermal contact re- tubes.Electric Power,2017,50(2):57 sistance.Sci China Technol Sci,2011,41(5):545 (葛铭，赵利杰，舒少辛，等.接触热阻对双H型鳍片管传 (王安良，赵剑锋.接触热阻预测的研究综述.中国科学：技 热系数的影响.中国电力，2017,50(2)：57) 术科学，2011,41(5)：545) [29]Liu Z F,Ma C Y,Zhao Y S,et al.Analysis of high speed motor- [47]Wang A L,Wu Y T,Yang C X.Size distribution characteristics ized spindle thermal characteristics based on thermal contact re- of cavities on a pool boiling surface./Eng Therm Energy Power, sistance.J Beijing Univ Technol,2016,42(1):17 2003,18(3):291 (刘志峰，马澄宇，赵永胜，等。基于接触热阻的高速精密电 (王安良，吴玉庭，杨春信。沸腾表面凹坑的尺度分布特征 主轴热特性分析.北京工业大学学报，2016,42(1)：17) 热能动力工程，2003,18(3)：291) [30] Cao T F,Mu Y T,Ding J,et al.Numerical investigation the [48]Singhal V,Garimella S V.Prediction of Thermal Contact Con- temperature distribution of PEM fuel cell with the effect of ther- ductance by Surface Deformation Analysis /International Me- mal contact resistance.J Eng Thermophys,2015,36(12):2661 chanical Engineering Congress and Exposition.New York,2000: (曹涛锋，母玉同，丁靖，等.接触热阻对燃料电池温度分布 43 影响的数值模拟.工程热物理学报，2016,36(12)：2661) [49]Xu L,Zhang T,Zhao L P,et al.Using double heat flux meter [31]Gibson R D.The contact resistance for a semi-infinite cylinder in method to measure the thermal contact resistance of solid material a vacuum.Appl Energy,1976.2(1):57 at low temperature and vacuum.Cryogenics,1999(4):185 [32]Fenech H,Rohsenow W M.Prediction of thermal conductance of (徐烈，张涛，赵兰萍，等。双热流法测定低温真空下固体界 metallic surfaces in contact.J Heat Transfer,1963,85(1):15 面的接触热阻.低温工程，1999(4)：185) [33]Cooper M G,Mikic BB,Yovanovich MM.Thermal contact con- [50]Fieberg C,Kneer R.Determination of thermal contact resistance ductance.Int Heat Mass Transfer,1969,12(3):279 from transient temperature measurements.Int Heat Mass Trans- [34]Bahrami M.Modeling of Thermal Joint Resistance for Sphere-Flat fr,2008,51(5-6):1017 Contacts in a Vacuum [Dissertation].Ontario:University of Wa- [51]Cong PZ,Zhang X,Fujii M.Estimation of thermal contact re- terloo,2004 sistance using ultrasonic waves.Int J Thermophys,2006,27(1): [35] Bahrami M,Yovanovich MM,Culham J R.A compact model for 171 spherical rough contacts.J Tribol,2005,127(4):884 [52] Shi L Investigation of Heat Transport on the Solid-Solid Contact [36]Willimams A.Major S.The solution of a steady state conduction Interface at Low Temperature Dissertation].Wuhan:Huazhong heat transfer problem using and electrolytic tank analogue.Me- University of Science and Technology,2006 chanical Eng Trans,1977,ME2:7 (石零.低温固-固接触界面热传输研究[学位论文]，武汉： [37]Yovanovich MM.Thermal constriction resistance of contacts on a 华中科技大学，2006) half-space:Integral formulation /AlAA 10th Thermophysics Con- [53]Zhang P,Xuan Y M,Li Q.Development on thermal contact re- ference.Denver,1975:397 sistance.CIESC J,2012,63(2):335 [38]Gibson R D.The contact resistance for a semi-infinite cylinder in (张平，宣益民，李强.界面接触热阻的研究进展.化工学 a vacuum.Appl Energy,1976.2(1):57 报，2012,63(2)：335) [39]Greenwood J A,Williamson J B P.Contact of nominally flat sur- [54]Lu MX.Nan G T.Zhou A G.Finite element analysis of thermal faces.Proc R Soc London Ser A,1966,295(1442):300 contact resistance based on ANSYS.Fluid Power Trans Control, [40]Yovanovich MM.Overall constriction resistance between contac- 2013(6):34 ting rough,wavy surfaces.Int J Heat Mass Transfer,1969,12 (陆敏恂，南国腾，周爱国.基于ANSYS的接触热阻的有限 (11):1517 元分析.流体传动与控制，2013(6)：34)
汪献伟等: 宏观接触热阻研究综述 Zhejiang Univ Nat Sci, 1997, 31(1): 104 (应济, 贾昱, 陈子辰, 等. 粗糙表面接触热阻的理论和实验 研究. 浙江大学学报:自然科学版, 1997, 31(1): 104) [26] Han Y G, Xuan Y M, Tang T F. Heat transfer on the contact surfaces with frictional heating. J Nanjing Univ Sci Technol, 1998, 22(3): 260 (韩玉阁, 宣益民, 汤瑞峰. 摩擦接触界面传热规律研究. 南 京理工大学学报, 1998, 22(3): 260) [27] Lin P Z, Sun B B, Zhao T. Thermal analysis and optimization of notebook computer based on thermal contact resistance. J South鄄 east Univ Nat Sci Ed, 2014, 44(1): 93 (林潘忠, 孙蓓蓓, 赵天. 基于接触热阻模型的笔记本电脑 散热性能分析与优化. 东南大学学报: 自然科学版, 2014, 44(1): 93) [28] Ge M, Zhao L J, Shu S X, et al. Impact of thermal contact re鄄 sistance on heat transfer coefficient of double H鄄type finned tubes. Electric Power, 2017, 50(2): 57 (葛铭, 赵利杰, 舒少辛, 等. 接触热阻对双 H 型鳍片管传 热系数的影响. 中国电力, 2017, 50(2): 57) [29] Liu Z F, Ma C Y, Zhao Y S, et al. Analysis of high speed motor鄄 ized spindle thermal characteristics based on thermal contact re鄄 sistance. J Beijing Univ Technol, 2016, 42(1): 17 (刘志峰, 马澄宇, 赵永胜, 等. 基于接触热阻的高速精密电 主轴热特性分析. 北京工业大学学报, 2016, 42(1): 17) [30] Cao T F, Mu Y T, Ding J, et al. Numerical investigation the temperature distribution of PEM fuel cell with the effect of ther鄄 mal contact resistance. J Eng Thermophys, 2015, 36(12): 2661 (曹涛锋, 母玉同, 丁靖, 等. 接触热阻对燃料电池温度分布 影响的数值模拟. 工程热物理学报, 2016, 36(12): 2661) [31] Gibson R D. The contact resistance for a semi鄄infinite cylinder in a vacuum. Appl Energy, 1976, 2(1): 57 [32] Fenech H, Rohsenow W M. Prediction of thermal conductance of metallic surfaces in contact. J Heat Transfer, 1963, 85(1): 15 [33] Cooper M G, Mikic B B, Yovanovich M M. Thermal contact con鄄 ductance. Int J Heat Mass Transfer, 1969, 12(3): 279 [34] Bahrami M. Modeling of Thermal Joint Resistance for Sphere鄄Flat Contacts in a Vacuum [Dissertation]. Ontario: University of Wa鄄 terloo, 2004 [35] Bahrami M, Yovanovich M M, Culham J R. A compact model for spherical rough contacts. J Tribol, 2005, 127(4): 884 [36] Willimams A. Major S. The solution of a steady state conduction heat transfer problem using and electrolytic tank analogue. Me鄄 chanical Eng Trans, 1977, ME2: 7 [37] Yovanovich M M. Thermal constriction resistance of contacts on a half鄄space: Integral formulation / / AIAA 10th Thermophysics Con鄄 ference. Denver, 1975: 397 [38] Gibson R D. The contact resistance for a semi鄄infinite cylinder in a vacuum. Appl Energy, 1976, 2(1): 57 [39] Greenwood J A, Williamson J B P. Contact of nominally flat sur鄄 faces. Proc R Soc London Ser A, 1966, 295(1442): 300 [40] Yovanovich M M. Overall constriction resistance between contac鄄 ting rough, wavy surfaces. Int J Heat Mass Transfer, 1969, 12 (11): 1517 [41] B俟rde S S, Yovanovich M M. Thermal resistance at smooth鄄 sphere / rough鄄flat contacts: Theoretical analysis / / 2nd AIAA / ASME Thermophysics and Heat Transfer Conference. Palo Alto, 1978: 83 [42] Shai I, Santo M. Heat transfer with contact resistance. Int J Heat Mass Transfer, 1982, 25(4): 465 [43] Williams A. Heat flow through single point of metallic contacts of simple shapes. AIAA / ASME Thermo鄄physics and Heat Transfer Conference. Boston, 1974:1 [44] Zahouani H, Vargiolu R, Loubet J L. Fractal models of surface topography and contact mechanics. Math Comput Modell, 1998, 28(4鄄8): 517 [45] Xu R P, Feng H D, Zhao L P, et al. Experimental investigation of thermal contact conductance at low temperature based on fractal description. Int Commun Heat Mass Transfer, 2006, 33(7): 811 [46] Wang A L, Zhao J F. Review of prediction for thermal contact re鄄 sistance. Sci China Technol Sci, 2011, 41(5): 545 (王安良, 赵剑锋. 接触热阻预测的研究综述. 中国科学: 技 术科学, 2011, 41(5): 545) [47] Wang A L, Wu Y T, Yang C X. Size distribution characteristics of cavities on a pool boiling surface. J Eng Therm Energy Power, 2003, 18(3): 291 (王安良, 吴玉庭, 杨春信. 沸腾表面凹坑的尺度分布特征. 热能动力工程, 2003, 18(3): 291) [48] Singhal V, Garimella S V. Prediction of Thermal Contact Con鄄 ductance by Surface Deformation Analysis / / International Me鄄 chanical Engineering Congress and Exposition. New York, 2000: 43 [49] Xu L, Zhang T, Zhao L P, et al. Using double heat flux meter method to measure the thermal contact resistance of solid material at low temperature and vacuum. Cryogenics, 1999(4): 185 (徐烈, 张涛, 赵兰萍, 等. 双热流法测定低温真空下固体界 面的接触热阻. 低温工程, 1999(4): 185) [50] Fieberg C, Kneer R. Determination of thermal contact resistance from transient temperature measurements. Int J Heat Mass Trans鄄 fer, 2008, 51(5鄄6): 1017 [51] Cong P Z, Zhang X, Fujii M. Estimation of thermal contact re鄄 sistance using ultrasonic waves. Int J Thermophys, 2006, 27(1): 171 [52] Shi L. Investigation of Heat Transport on the Solid鄄鄄 Solid Contact Interface at Low Temperature [Dissertation]. Wuhan: Huazhong University of Science and Technology, 2006 (石零. 低温固鄄鄄固接触界面热传输研究[学位论文]. 武汉: 华中科技大学, 2006) [53] Zhang P, Xuan Y M, Li Q. Development on thermal contact re鄄 sistance. CIESC J, 2012, 63(2): 335 (张平, 宣益民, 李强. 界面接触热阻的研究进展. 化工学 报, 2012, 63(2): 335) [54] Lu M X, Nan G T, Zhou A G. Finite element analysis of thermal contact resistance based on ANSYS. Fluid Power Trans Control, 2013(6): 34 (陆敏恂, 南国腾, 周爱国. 基于 ANSYS 的接触热阻的有限 元分析. 流体传动与控制, 2013(6): 34) ·1247·
·1248· 工程科学学报.第41卷，第10期 [55]Liu D H,Wang F,Zeng F W,et al.Finite element simulation fastening rotor.Turbine Technol,2018,60(4):267 method of high temperature thermal contact resistance.Eng (李俊南，张锁怀，吕超颖，等。分布式拉杆转子轮盘结合面 Mech.2012.29(9):375 接触热阻建模与分析.汽轮机技术，2018,60(4)：267) (刘冬欢，王飞，曾凡文，等.高温接触热阻的有限元模拟方 [63] Jiang L,Jiang B Y,Wu W Q,et al.Research on thermal con- 法.工程力学，2012,29(9)：375) tact resistance between polymer and mold during thin-wall injec- [56]Yu Y F,Li Z W,Zheng X Y.The inverse problem of thermal tion molding.J Cent South Univ Sci Technol,2018.49(6): contact resistance between rough surfaces.Chin Theor Appl 1381 Mech,2018,50(3):479 (蒋礼，蒋炳炎，吴旺青，等.薄壁注射成型中聚合物与模具 (余远锋，李泽卫，郑晓亚.粗糙表面之间接触热阻反问题 间接触热阻研究.中南大学学报（自然科学版），2018,49 研究.力学学报，2018,50(3)：479) (6):1381) [57]Yang H W,Yuan GJ,Xiao L.Application of ANSYS contact el- [64]Hu Y,Wang C,Xu F X,et al.Research on thermal contact re- ement in simulation of thermal contact resistance.Microuares, sistance for boron steel sheet in hot stamping process.Forg Stamp 2012(Suppl2):241 Technol,.2018,43(7):147 (杨华威，袁广江，肖刘.ANSYS接触单元在接触热阻仿真 (胡翼，王超，徐峰祥，等.硼钢板热冲压成形过程接触热阻 中的应用.微波学报，2012（增刊2）：241) 研究.锻压技术，2018,43(7)：147) [58]Yang H W,Yuan G J,Xiao L,et al.Simulation of thermal con- [65]Wang C H,Zhou J,Zhang X Y.The research on thermal contact tact resistance in collector of traveling wave tube.Chin Vac Sci resistance and contact resistance of solid at low temperature. Technol,.2013,33(6):517 Cryog Supercond,2018,46(3):1 (杨华威，袁广江，肖刘，等.使用ANSYS接触单元模拟研 (王存洪，周军，张兴义.低温下固体接触热阻与接触电阻 究行波管收集极的接触热阻.真空科学与技术学报，2013， 的实验研究.低温与超导，2018,46(3)：1) 33(6):517) [66]Johnson K L.Contact Mechanics.Cambridge:Cambridge Univer- [59]Cao J,Hu P J,Ying J.Finite element analysis of spindle thermal sity Press,1987 character based on thermal contact resistance.Mech Electr Eng [67]Sridgar M R,Yovanovich MM.Elastoplastic Constriction Resist- Mag,2008,25(2):20 ance of Sphere/flat Contacts.J Heat Trans.,1996,118(1): (曹骏，胡佩俊，应济.基于接触热阻的主轴热特性有限元 202 分析.机电工程，2008,25(2)：20) [68]Roca R T,Mikic R E.Thermal Contact Resistance in a Non-ideal [60]Zhan L H,Li X Q,Hu S C.Experiment research for the interfa- Joint.Massachusetts:Massachusetts Institute of Technology, cial thermal contact resistance.Light Alloy Fabric Technol,2002, 1971 30(9):40 [69]Albin K J,Hasselstrom U,Nileeon E.Thermal Contact Conduct- (湛利华，李晓谦，胡仕成.界面接触热阻影响因素的实验 ance in a Bolted Joint [Dissertation ]Gothenburg:Chalmers U- 研究.轻合金加工技术，2002.30(9)：40) niversity of Technology,2012 [61]Wang X W,Hu P.Liu Y.Study on experimental device and [70]Bevans JT,Ishimoto T,Loya B R,et al.Prediction of Space Ve- measurements for interface thermal resistance of adhesive struc- hiclet Thermal Characteristics.Ohio:Air Force Flight Dynamics ture of CFRP.J Unin Sci Technol China,2018,48(7):600 Laboratory,1965 (王向往，胡芃，刘阳.界面热阻实验装置与碳纤维增强板 [71]Abadi P PS S,Leong C K,Chung DD L.Factors that govern 胶接热阻研究.中国科学技术大学学报，2018,48(7)：600) the performance of thermal interface materials.Electron Mater, [62]Li J N,Zhang S H,Lii C Y,et al.Modeling and analysis of 2009,38(1):175 thermal contact resistance on joint interfaces of distributed rod
工程科学学报,第 41 卷,第 10 期 [55] Liu D H, Wang F, Zeng F W, et al. Finite element simulation method of high temperature thermal contact resistance. Eng Mech, 2012, 29(9): 375 (刘冬欢, 王飞, 曾凡文, 等. 高温接触热阻的有限元模拟方 法. 工程力学, 2012, 29(9): 375) [56] Yu Y F, Li Z W, Zheng X Y. The inverse problem of thermal contact resistance between rough surfaces. Chin J Theor Appl Mech, 2018, 50(3): 479 (余远锋, 李泽卫, 郑晓亚. 粗糙表面之间接触热阻反问题 研究. 力学学报, 2018, 50(3): 479) [57] Yang H W, Yuan G J, Xiao L. Application of ANSYS contact el鄄 ement in simulation of thermal contact resistance. J Microwaves, 2012(Suppl 2): 241 (杨华威, 袁广江, 肖刘. ANSYS 接触单元在接触热阻仿真 中的应用. 微波学报, 2012(增刊 2): 241) [58] Yang H W, Yuan G J, Xiao L, et al. Simulation of thermal con鄄 tact resistance in collector of traveling wave tube. Chin J Vac Sci Technol, 2013, 33(6): 517 (杨华威, 袁广江, 肖刘, 等. 使用 ANSYS 接触单元模拟研 究行波管收集极的接触热阻. 真空科学与技术学报, 2013, 33(6): 517) [59] Cao J, Hu P J, Ying J. Finite element analysis of spindle thermal character based on thermal contact resistance. Mech Electr Eng Mag, 2008, 25(2): 20 (曹骏, 胡佩俊, 应济. 基于接触热阻的主轴热特性有限元 分析. 机电工程, 2008, 25(2): 20) [60] Zhan L H, Li X Q, Hu S C. Experiment research for the interfa鄄 cial thermal contact resistance. Light Alloy Fabric Technol, 2002, 30(9): 40 (湛利华, 李晓谦, 胡仕成. 界面接触热阻影响因素的实验 研究. 轻合金加工技术, 2002, 30(9): 40) [61] Wang X W, Hu P, Liu Y. Study on experimental device and measurements for interface thermal resistance of adhesive struc鄄 ture of CFRP. J Univ Sci Technol China, 2018, 48(7): 600 (王向往, 胡芃, 刘阳. 界面热阻实验装置与碳纤维增强板 胶接热阻研究. 中国科学技术大学学报, 2018, 48(7): 600) [62] Li J N, Zhang S H, L俟 C Y, et al. Modeling and analysis of thermal contact resistance on joint interfaces of distributed rod fastening rotor. Turbine Technol, 2018, 60(4): 267 (李俊南, 张锁怀, 吕超颖, 等. 分布式拉杆转子轮盘结合面 接触热阻建模与分析. 汽轮机技术, 2018, 60(4): 267) [63] Jiang L, Jiang B Y, Wu W Q, et al. Research on thermal con鄄 tact resistance between polymer and mold during thin鄄wall injec鄄 tion molding. J Cent South Univ Sci Technol, 2018, 49 ( 6 ): 1381 (蒋礼, 蒋炳炎, 吴旺青, 等. 薄壁注射成型中聚合物与模具 间接触热阻研究. 中南大学学报(自然科学版), 2018, 49 (6): 1381) [64] Hu Y, Wang C, Xu F X, et al. Research on thermal contact re鄄 sistance for boron steel sheet in hot stamping process. Forg Stamp Technol, 2018, 43(7): 147 (胡翼, 王超, 徐峰祥, 等. 硼钢板热冲压成形过程接触热阻 研究. 锻压技术, 2018, 43(7): 147) [65] Wang C H, Zhou J, Zhang X Y. The research on thermal contact resistance and contact resistance of solid at low temperature. Cryog Supercond, 2018, 46(3):1 (王存洪, 周军, 张兴义. 低温下固体接触热阻与接触电阻 的实验研究. 低温与超导, 2018, 46(3):1) [66] Johnson K L. Contact Mechanics. Cambridge: Cambridge Univer鄄 sity Press, 1987 [67] Sridgar M R, Yovanovich M M. Elastoplastic Constriction Resist鄄 ance of Sphere / flat Contacts. J Heat Trans. , 1996, 118 (1): 202 [68] Roca R T, Mikic R E. Thermal Contact Resistance in a Non鄄ideal Joint. Massachusetts: Massachusetts Institute of Technology, 1971 [69] Albin K J, Hasselstrom U, Nileeon E. Thermal Contact Conduct鄄 ance in a Bolted Joint [Dissertation]. Gothenburg:Chalmers U鄄 niversity of Technology, 2012 [70] Bevans J T, Ishimoto T, Loya B R, et al. Prediction of Space Ve鄄 hiclet Thermal Characteristics. Ohio: Air Force Flight Dynamics Laboratory, 1965 [71] Abadi P P S S, Leong C K, Chung D D L. Factors that govern the performance of thermal interface materials. J Electron Mater, 2009, 38(1): 175 ·1248·