GIS设备用表带触指电连接结构接触电阻的分形理论计算与分析.pdf

面向双碳目标的环保型电力设备DOI:10.15961/j.jsuese.202200890GIS设备用表带触指电连接结构接触电阻的分形理论计算与分析杨为1，朱太云1，任汀2，王青于2*，靳守锋2(1.国网安徽省电力有限公司电力科学研究院，安徽 合肥 230601；2.西安交通大学 电气绝缘与电力设备国家重点实验室，陕西 西安 710049)摘要：电连接结构的接触电阻（electrical contact resistance，ECR）是评价电接触系统可靠性的重要指标。表带触指是环保气体GIS设备中重要的电连接部件，其ECR不仅影响回路总电阻，通流过大导体振动、安装不当等因素也会导致电接触不良，引发触指局部过热，引起绝缘放电等故障，还会关系到环保气体GIS设备的运行稳定性。为研究触指ECR的产生机理和变化规律，采用MATLAB和激光共聚焦显微镜捕捉接触区域灰度图像，获得触指接触区域的名义接触面积Aa和分形维数D；运用分形理论模型，数值计算了表带触指电连接结构的ECR，并通过ECR试验测量，验证了应用分形理论模型解析计算表带触指电连接结构ECR的合理性和准确性；研究了表带触指电接触过程中的弹性变形、第一弹塑性变形、第二弹塑性变形和塑性变形4个阶段，对比分析了不同分形维数对触指电接触实际接触面积的影响规律。研究发现：分形理论模型可用于GIS设备用表带触指的接触电阻理论计算；触指电接触表面形貌越复杂，分形维数越大，触点变形越先进入弹塑性阶段，也越易进入第二弹塑性变形阶段，达到相同承载能力所需实际接触面积也越小。研究成果有助于环保气体GIS电连接结构的优化设计，可为环保气体GIS设备触指电连接结构的多物理场计算和电接触性能的状态评估提供理论基础。关键词：分形模型；表带触指；接触电阻；弹塑性变形；实际接触面积中图分类号：TM85文献标志码：A文章编号：2096-3246（2023）04-0011-10Calculation and Analysis of the Electrical Contact Resistance of the Strap Contacts Used inGIS Equipment Based on the Fractal ModelYANG Wei1，ZHU Taiyun1，REN Ting2，WANG Qingyu2*，JIN Shoufeng2(1.State Grid Anhui Electric Power Research Inst.,Hefei 230601,China;2.State Key Lab.of Electrical Insulation and Power Equipment,Xian Jiaotong Univ.,Xian 710049,China)Abstract:Contact resistance of an electrical connection structure(ECR)is an important indicator to evaluate the reliability of the electrical con-tact system.The strap contact is an important electrical connection component in eco-friendly insulating gas GIS equipment.Its ECR will affectthe total resistance of the equipment.Factors such as the excessive current flow,the conductor vibration,and the improper installation can alsoleads to poor electrical contact,local overheating of the strap contacts,insulation discharge and other failures,which is directly related to the oper-ation stability of environmental protection gas GIS equipment.To study the mechanism and regularity of ECR,MATLAB and laser confocal mi-croscope were used to capture the grayscale image of the contact area in this paper.The nominal contact area Aa and fractal dimension D of thecontact area were obtained.Using the fractal theoretical model,the contact resistance of the strap contact was numerically calculated.The ration-ality and accuracy of applying the fractal theoretical model to analytically calculate the contact resistance of the electrical connection structure was收稿日期:2022 08 24基金项目:国网安徽省电力有限公司科技项目（B3120521000Z）作者简介:杨为（1984），男，高级工程师.研究方向：变电设备运维检测技术与标准修制订.E-mail：hfut_*通信作者:王青于，E-mail：网络出版时间:2023 07 06 10:53:59 网络出版地址:https:/ http:/http:/ 第 55 卷 第 4 期工 程 科 学 与 技 术Vol.55 No.42023 年 7 月ADVANCED ENGINEERING SCIENCESJuly 2023verified through ECR test measurements.Since ECR is closely related to the actual contact area,the elastic deformation,the first elastic-plasticdeformation,the second elastic-plastic deformation and the plastic deformation during the electrical contact process of the strap contact were stud-ied.Whats more,the effect of different fractal dimensions on the actual contact area of the electrical contact was studied.The research found thatthe fractal theoretical model can be used for the theoretical calculation of ECR for the strap contact of the GIS equipment.It was obtained that thelarger the fractal dimension of the electrical contact surface,the earlier the contact deformation enters the elastic-plastic stage,and the easier it isto enter the second elastic stage.In the plastic deformation stage,the more complex the surface morphology,the smaller the actual contact area re-quired to achieve the same bearing capacity.The research results are helpful for the optimal design of the eco-friendly insulating gas GIS electric-al connection structure,and can provide a theoretical basis for the multi-physics calculation of the eco-friendly insulating gas GIS equipment elec-trical connection structure and the state evaluation of the electrical contact performance.Key words:fractal model;strap contact;contact resistance;elastic-plastic deformation;actual contact area“双碳”目标和背景下，环保气体绝缘开关设备GIS发展前景广阔1。GIS设备载流导体长径比大，金属外壳与载流导体膨胀收缩产生相对位移差，通常采用插接式滑动触头进行补偿，故而电接触部件成为电连接结构的关键部件。然而，GIS设备运行条件复杂苛刻，且随着中国特高压输电技术的不断发展，一次设备的运行电流大幅增加，触指电连接结构成为影响电力设备安全可靠运行的关键环节和薄弱部位24。GIS中心载流导体分段处常采用弹簧、表带和梅花等类型的触指电连接结构。表带触指多圈同向安装在公头外侧或母头内壁上，触指叶片之间靠不锈钢龙骨连接并保持一定间距。表带触指能够在一定范围接触压力条件下保持良好的机械和电接触性能。GIS设备的公头与母头会因膨胀收缩产生周期往复运动2；同时，GIS导体会因电磁力作用发生振动，引起各个触指叶片通流不均匀，严重时将引发局部过热和电接触故障。近年来，GIS设备中大量使用的触指电连接结构因电接触区的性能劣化导致的事故时有发生，严重威胁了电网的安全稳定运行（图1）。因此，电连接结构的可靠性至关重要，其中，接触电阻是评价此可靠性的一项重要指标。接触电阻源于不规则触指表面与导电杆的接触面之间通电后产生的导电斑点5。当电连接结构通流，电流流过触指表面众多尺寸较小的导电斑点时，发生大幅收缩，产生收缩电阻5。除此之外，触指与导体的显微表面凹凸不平，接触面是众多导电斑点构成的集合，并非理想平面接触，故触指电连接结构的实际接触面积与名义接触面积有所不同，其实际接触面积与触指和导体表面的粗糙度、接触压力载荷等因素相关5。由于导电斑点微观结构不规则，经典数学难以表征，随着非线性科学的发展，分形理论被应用于微观层面，可用于建立更为精准的粗糙表面接触模型68。Mandelbrot于1975年描述了分形几何，后续建立了分形理论用以研究分形特征和应用9。Majum-dar等1011开发了一种新的各向同性粗糙表面之间的电接触理论模型，发现接触面积小于临界接触面积的所有接触点都是塑性接触；当接触载荷增加时，这些塑性变形点会连接到一起而形成弹性接触，但所用分形粗糙参数与尺度无关。杨红平等12基于分形理论，表征了粗糙微凸体表面，计算获得了微凸体在不同塑性指数下的接触刚度和接触载荷间的关系。Liou等13建立了球体和圆柱体表面微凸体的分形参数模型，推导了弹性、弹塑性和塑性变形3个阶段的尺寸分布函数，得到接触微凸体弹性、弹塑性和完全塑性的变形机制与经典（统计）理论描述一致，在各种情况下预测的接触载荷和实际接触面积结果与使用统计方法获得的结果非常吻合。综上所述，国内外研究学者针对接触电阻计算的分形理论模型开展了大量工作，为分形理论模型在电接触领域的应用奠定了坚实的基础14。然而，有关实际GIS设备用表带触指接触电阻的理论计算报道较少。本文基于分形理论模型，运用MATLAB识别触指接触区域灰度图像，获得了表带触指接触面的分形维数D，试验对比分析了分形理论应用于接触电阻理论计算的合理性，并研究了不同分形维数D对表带触指电接触面实际接触面积和变形状态的影响。研究成果有助于环保气体GIS电连接结构的优化设计，可为环保气体GIS设备触指电连接结构的多物理 机械损伤异常发热磨损放电烧蚀刮痕2刮痕1疲劳变形A相触指有电灼烧痕迹图 1GIS设备用触指电连接结构故障现场Fig.1Fault scene of the electrical connection structure ofthe contact finger for GIS equipment 12工程科学与技术第 55 卷场计算和电接触性能的状态评估提供理论基础。1 接触电阻分形计算 1.1 分形理论模型粗糙面表面轮廓曲线测试表明粗糙面具有复杂而精细的结构，经典的几何语言很难描述其局部和整体。分形理论通过定义分形维数D提供了一种描述表面粗糙度的方法，相比较统计方法，其对粗糙表面进行表征和建模的准确性较好。分形维数D表示图形的维数，其计算公式为：D=lg Nlg（1）式中，N表示各个区域的像素数量，表示计算区域和原图的相对大小。针对粗糙表面的研究，分形维数的计算方法众多，其中以均方根和结构函数法较为常见，后者计算精度高9。然而，表带触指的接触区域十分不规则，形状类似圆柱面，通过计算式（1）无法表征触指触点的真实形貌，因此，通过结合图像识别技术，获取更接近真实形貌的分形维数D。但结构完全不同的表面也有可能具有相近甚至相同的分形维数D，说明仅用分形维数D无法唯一确定表面。Majumdar等15运用结构参数G表示粗糙表面的特征尺度参数，以进一步完善对表面微观形貌的描述。在获取分形维数D和结构参数G的基础上，Kogut等16提出一种接触电阻计算的综合理论，获得了已知材料特性和各向同性，均匀导电表面的接触电阻与接触载荷和接触面积之间的关系，为计算触指静态接触电阻提供了重要的理论依据。1.2 面积分布函数实际触指电连接结构是触指粗糙表面和公头/母头导体粗糙表面之间的电接触，为简化计算，可将两个粗糙表面的接触等效为一个刚性的光滑平面接触与一个粗糙平面接触17。MB分形接触模型18给出了在分形粗糙表面接触理想的刚性平面时，微凸体接触面积的分布密度函数为：n(a)=D2(2D)/2a1D/2a(D+2)/2（2）实际接触面积Ar为：Ar=wa10n(a)ada=D2D(2D)/2a1（3）式（2）（3）中，a为接触面积，al为微凸体接触面积最大值，为分形区域扩展系数。在已知分形维数的情况下，分形区域扩展系数可通过方程（4）计算求得10：(2D)/2(1+D/2)(2D)/D(2D)/D=1（4）1.3 微凸体变形微凸体结构的变形方式可分为3个阶段：弹性阶段、弹塑性阶段和完全塑性变形阶段。图2为单个微凸体在承载接触压力前后的变形轮廓示意图。图2中，l为微凸体底部的尺寸（基底长度），l为微凸体变形量为时变形后微凸体与初始未变形微时凸体交点间的长度，lr为微凸体变形量为时与刚性平面的实际接触长度。电连接结构的表面均进行过机加工，表面形貌精细到原子尺度仍具有精细的自相似结构。分形几何学的WeierstrassMandelbrot分形函数（简称WM函数）满足粗糙表面自仿射性及其他相关数学特性，因此可以利用2维WM函数表达某一具体的粗糙表面轮廓，数学模型如式（5）所示9：z(x)=GD1l2Dcosxl,l2 x l2（5）根据式（5）可推导出微凸体顶端曲率半径为：=?(1+z2)3/2z?x=0=lD2GD1（6）微凸体的高度为：=z(0)=GD1l2D（7）微凸体的实际变形量为：=z(l2)=GD1l2D1cos(l2l)（8）1.3.1 弹性变形阶段当微凸体承载压力较小，根据Herz弹性接触理论19，微凸体的最大接触压力P0为：P0=2E()12（9）式中，E为等效弹性模量，与电连接两个接触体材料 zxlrOll图 2承载接触压力前后微凸体轮廓示意图Fig.2Schematic diagram of asperity profile before andafter bearing contact pressure 第 4 期杨为，等：GIS设备用表带触指电连接结构接触电阻的分形理论计算与分析13的弹性模量和泊松比有关。由Herz理论还可得到微凸体发生弹性变形时，实际接触面积a为：a=（10）接触载荷F为：F=43E1232（11）根据文献15，半径为的微凸体与光滑刚性平面接触，临界弹性接触变形量为：ec=(K2)2（12）=H/E式中：为衡量微凸体变形能力的参数，H为较软材料的硬度；K为材料的磨损系数，K=2.8。将式（6）和（12）代入式（10）中可以得到微凸体临界弹性接触面积为：ec=1(K2)2(lDGD1)2（13）将式（6）、（11）和（13）结合，可以得到当单个微凸体发生弹性变形时（临界弹性变形也适用），接触载荷F与接触面积a之间的关系为：F=4E1/2GD13lDa3/2（14）1.3.2 弹塑性变形阶段ec 6ec6ec 110ecKogut等16仿真计算了不同分形维数D和结构参数G时粗糙表面的接触电阻，得到微凸体第一弹塑性变形区为，第二弹塑性变形区为；在这两个区域形成的关键临界接触面积分别为：aepc=0.93ec(6ecec)1.136=7.119 7aec（15）apc=0.94ec(110ecec)1.146=205.382 7aec（16）式中，aepc 和apc分别为第一弹塑性、第二弹塑性的临界接触面积。1.3.3 完全塑性变形阶段ec当微凸体的实际变形量大于110时，微凸体发生完全塑性变形，接触载荷F与面积a之间的关系为：F=Ha（17）此时，接触面积为：a=2（18）1.4 分形条件的接触电阻计算分析第1.3节可知，触指粗糙表面和电极之间的接触载荷与实际接触面积分别是所有微凸体接触载荷之和与接触面积之和，而根据微凸体尺寸和承受载荷的不同，微凸体的变形分为弹性变形、第一弹塑性变形、第二弹塑性变形和完全塑性变形。因此，实际接触面积是以上变形接触面积的总和。当实际接触面积处于相应区间范围内时，即可唯一求得此时的实际接触面积和接触载荷。实际接触面积的计算方法为：Ar=Are+Arep1+Arep2+Arp（19）Are=waec0n(a)ada=D2D(2D)/2aec(2D)/2aD/21（20）Arep1=waepcaecn(a)ada=D2D(2D)/2a(2D)/2epca(2D)/2ecaD/21（21）Arep2=wapcaepcn(a)ada=D2D(2D)/2a(2D)/2pca(2D)/2epcaD/21（22）Arp=wa1apcn(a)ada=D2D(2D)/2a(2D)/21a(2D)/2pcaD/2（23）式中，Are为弹性变形接触面积，Arep1为第一弹塑性变形接触面积，Arep2为第二弹塑性变形接触面积，Arp为完全塑性变形接触面积。实际接触载荷Fr的计算方法为：Fr=Fre+Frep1+Frep2+Frp（24）式中，Fre为弹性变形接触载荷，Frep1为第一弹塑性变形接触载荷，Frep2为第二弹塑性变形接触载荷，Frp为完全塑性变形接触载荷。需要说明的是，由于微凸体基底长度l随着接触面积变化而变化，即对于不同的微凸体，基底长度l不同，且基底长度l与接触面积之间的关系并非已知线性关系，接触微凸体总基底长度的积分无法计算。G2/(K/2)2/(D1)为了解决上述难题，令临界弹性接触面积等于MB分形接触模型17给出的临界接触面积，可求得临界基底长度lec。这样一来，在计算所有弹性接触微凸体的接触载荷时，将基底长度视为常数lec，此时的临界基底长度相当于一个基准值描述所有微凸体，只是不同尺寸的微凸体的基底长度与临界基底长度相差某一个倍数，在计算基底长度积分时，可以等效为lec的积分。通过计算发现，的取值只影响载荷数量级，并不影响载荷变化趋势，此时，只需要辅以试验确定的取值，就能求解出实际接触面积和载荷之间的关系：Fre=waec0Fn(a)ada=D3D(2D)/24E1/2GD13lDa(3D)/2ecaD/21（25）14工程科学与技术第 55 卷Frep1=waepcaecFn(a)da=D3KH(2D)/21.281 6aec0.254 41.4250.568D(a1.254 4epc0.5Da1.254 4ec0.5D)aD/21（26）Frep1=wapcaepcFn(a)da=D3KH(2D)/21.717a0.102 1ec1.2630.573Da1.102 10.5Dpca1.102 1epc0.5DaD/21（27）Frp=wa1apcFn(a)da=HD2D(2D)/2a(2D)/21a(2D)/2pcaD/21（28）l=lmaxAalmaxAr=ArAaaec=aecAa为消除量纲影响，将式（25）（28）进行归一化处理，两边同时除以EAa。则有：，为最大微凸体尺寸；。因此，当材料属性、分形维数D、结构参数G保持不变时，统一量纲后的接触载荷是最大接触面积和实际接触面积的函数。对以上不同变形阶段接触载荷加和，即可得到电接触表面总接触载荷与实际接触面积的函数表达式：Fr=41/2G(D1)3lDg1(D)a(3D)/2ecArD/2+D3Kg2(D)1.281 6a0.254 4ec1.4250.568D(7.119 7aec)(1.254 40.5D)a(1.254 40.5D)ecArD/2+D3Kg2(D)1.717a0.102 1ec1.2630.573D(205.382 7aec)(1.102 10.5D)(7.119 7aec)(1.102 10.5D)ArD/2+g3(D)g3(D)Ar(2D)/2(205.382 7aec)(2D)/2ArD/2,1 D 2（29）g1(D)=D3D(2DD)D/2(D2)42g2(D)=(2DD)D/2(D2)24g3(D)=(2DD)D/2(D2)24式中，。综上，结合分形理论及经验公式推导得到了接触载荷Fr、实际接触面积Ar的内在关系。由于施加载荷P和接触载荷Fr相等，因此，通过与式（29）形成无量纲等式，可以求解触指电接触面各个变形阶段的接触载荷，最终实现接触电阻的理论计算。结合Holm接触电阻计算经验公式6，推导出分形模型计算接触电阻的公式如下：Rce=2Are/=(2D)2DaD/2la(1D/2)ec(1D/2)（30）Rcep1=(2D)2DaD/2l(a(1D/2)epca(1D/2)ec)(1D/2)（31）Rcep2=(2D)2DaD/2l(a(1D/2)pca(1D/2)epc)(1D/2)（32）Rcp=(2D)2DaD/2l(a(1D/2)la(1D/2)pc)(1D/2)（33）Rc=11/Rce+1/Rcep1+1/Rcep2+1/Rcp（34）式中，Rce为弹性变形产生的接触电阻，Rcep1为第一弹塑性变形产生的接触电阻，Rcep2为第二弹塑性变形产生的接触电阻，Rcp为完全塑性变形产生的接触电阻，Rc为总接触电阻。2 表带触指电接触微观形貌参数和接触电阻试验测量 2.1 表带触指电接触微观形貌参数测量由上述计算过程可知，接触面参数不变时，接触载荷直接决定了电接触面的变形区域及接触面积的大小，且二者存在一一对应关系。由此设计并开展了GIS设备用表带触指触点区域的微观形貌参数测量，为表带触指电连接结构接触电阻的理论计算提供基础参数。试验对象为GIS设备用LA-CUD型表带触指，材料为铍铜表面镀银，镀层厚20 m20。图3为激光共聚焦显微镜（型号为奥林巴斯LEXT OLS4100）拍摄得到的接触区域3维实拍图与高度图。其中，图3（a）为 500 m500 m表带触指顶部触点触点实拍图触点亮度图Y/mZ/m10886430318637X/m4468921 3391 785(a)触指实物和触点放大图(b)触指接触区域高度和轴测图图 3表带触指实物和接触区域3维实拍图Fig.33D real images of the physical object and contact area 第 4 期杨为，等：GIS设备用表带触指电连接结构接触电阻的分形理论计算与分析15触指实物照片和触点放大图，图3（b）为激光共聚焦显微镜自带软件对扫描结果分析得到的高度和轴测图。由图3（a）可知：接触过程中，触点高度越高，表面越平滑，在灰度图像中更接近白色；当表面出现凹坑或斜坡时，在灰度图像中就更接近黑色，因此可以根据灰度图像获取接触区域的分形维数D。由图3（b）可知，触点接触前后均为不规则表面，整个接触区域存在22 m的高度差，并且测得单片触指接触区域的投影面积约为1.3107 m2。由于现有手段难以获得准确的名义接触面积，且投影面积与名义接触面积数量级相近，故将投影面积近似为名义接触面积Aa。分形维数D的获取步骤如下：1）运用激光共聚焦显微镜自带软件对扫描结果进行分析，获得实拍、彩色、灰度图像；2）截取256256像素的灰度图像，运用MATLAB软件对图像进行识别，绘制对应的3维立体图；3）运用差分盒子法计算分形维数D。差分盒子法的原理是将一幅大小为MN的图像分为大小相等的若干kk子块8，图像(x,y)处的灰度值为f(x,y)，总的灰度级为Q。此时，将图像视为3维物体的表面灰度集(x,y，f(x,y)；XY平面上是kk的网格，每个网格上有若干个盒子叠加，盒子高度h=(Q1)k/min(M,N)。若在第(i,j)个网格中，第m个盒子中所包含网格内的灰度为最小值，第l个盒子中所包含网格内的灰度为最大值，则覆盖第(i,j)个网格的盒子数nr(i,j)=lm+1，覆盖全部网格的盒子数Nr=nr(i,j)。可求得分形维数为：D=limlg Nrlg(1/r)（35）式中，r=k/min(M,N)。通过改变网格k的大小计算1组Nr，对lg(1/r)和lg Nr进行线性回归，其斜率即分形维数D。根据接触区域多个位置的灰度图像获取分形维数D，取平均值后得到本文研究对象表带触指接触区域分形维数D=1.5。2.2 表带触指电连接结构接触电阻测量 2.2.1 测量系统搭建表带触指接触电阻测量系统示意图和实物图分别如图4、5所示。测量系统主要由表带触指、上下极板、压力加载装置和直流低电阻测量装置组成。在下极板开槽，将表带触指放入槽内，分别与上下极板电接触，触指与触指之间并联连接。压力加载系统是指在表带触指和电极之间施加正压力的系统，其加载器采用气缸结构，可通过上支撑杆对上电极施加最大值为800 N的正压力。上支撑杆与上电极采用滚动轴承活动连接，可保证上电极在往复运动的过程中受到均匀向下的正压力。压力传感器串联在上支撑杆和上电极之间，可对正压力的变化进行实时测量。由于接触电阻为微欧级别，选取QJ36S-2直流低电阻测试仪进行测量。整个系统放置在一个可封闭的不锈钢腔体中。具体试验步骤为：1）将下电极安装在基座中，拧紧螺丝完全固定，将210片表带触指安装在卡槽中，上电极安装在水平作动杆上，电子水平仪保证水平。2）将热电偶布置在上下电极槽内，测量位置靠近触指触点。3）打开气缸，调节上支撑杆，调整压力传感器显示数值为20 N。4）关闭不锈钢腔门，抽真空后向罐体内注入表压为0的氮气作为保护气，避免因金属电极发热而加速接触区域氧化。5）打开直流低电阻测试仪电源，将正负极短接，预热15 min。6）预热完成后，将正负极连接在接线端子上，待示数稳定后记录电阻值。7）改变载荷从20 N逐级升高至220 N，间隔为10 N，表带触指压力传感器固定支撑直流低电阻测试仪载荷滚动轴承图 4接触电阻测量系统示意图Fig.4Schematic diagram of contact resistance measure-ment system 压力加载系统上电极下电极图 5接触电阻测量系统实物Fig.5Physical of contact resistance measurement system 16工程科学与技术第 55 卷即平均到每片触指载荷的变化范围为1至11 N，间隔为0.5 N。每次升高正压力后间隔30 min测量1次数据。试验测得的电阻值为上电极体电阻、下电极体电阻、表带触指体电阻和总接触电阻之和。由于表带触指体电阻可忽略不计21，因此，将测得的总电阻减去上、下电极的体电阻，即可得到总的接触电阻。2.2.2 测量误差消除触指接触电阻测量对精度要求高，测量过程中的误差来源有3个方面：一是螺栓和电极连接接触面的接触电阻，这是无法避免但可以被抵消的；二是直流低电阻测试仪自身的系统误差；三是通流后，接触电阻随接触面温度升高而变化。为消除以上误差，首先，利用相同的接线装置测量上、下电极的体电阻，两者相减后测量误差和毫伏表的系统误差也被减去，从而避免结果的分散性。其次，每次试验前保证直流低电阻测试仪充分预热，测量过程间隔保证一致。但是，在电接触理论中，触点的温升难以通过试验测量，为消除温度误差，通常采用VT关系计算触点的温升22。Konchits等23建立了电接触中的VT关系，在电热平衡且与周围环境电绝缘和热隔离的情况下，该关系对于同类材料的任意形状电接触皆有效：T3=U2c4L+T21（36）式中：Uc为接触电势；T1为基体温度；L为Wiedemann-Franz常数，其值为2.44108 V2/K2；T3为触点温度。通过计算得到触点相比表带触指基体温度升高了0.22，结合热电偶测量基体温度校正之前相减得到的接触电阻，可得到最终的接触电阻测量值。3 表带触指电连接结构接触电阻理论分析 3.1 试验与理论对比表带触指材料为铜基表面镀银，电极材料为铝合金，理论计算时基于接触材料特性设置两个接触表面的相关参数，见表1。由于单片触指额定载荷约为5 N20，因此设定载荷为110 N。采用本文分形理论模型数值计算得到了一组不同载荷下单片触指的接触电阻Rc，并将第2节试验测量的表带触指接触电阻与之对比，结果如图6所示。工程实际中，表带触指单片合适压力为2.55.0 N之间，当单片压力为5 N时，触指额定ECR为320 20。由于测量误差等原因，试验得到的触指ECR为355，相对误差11%，表明表带触指电连接结构的接触电阻试验测量结果合理。由图6（a）可知，试验测量和分形模型计算得到的表带触指ECR变化趋势相近，分形理论模型计算结果整体高于试验测量值。当所受载荷为110 N时，分形理论模型计算得到表带触指电接触区域处于第一弹塑性变形阶段，接触电阻由弹性和弹塑性变形两个阶段组成。当单片触指载荷为13 N时，整体接触电阻偏高，之后随着载荷增大，接触面积增加，接触电阻缓慢减小，变化趋势与触指接触电阻试验测量结果基本一致。由图6（b）可知，分形理论计算结果相对误差约为20%。当载荷为13 N时，相对误差较小，约为10%；表 1分形数值计算用粗糙表面参数Tab.1 Rough surface parameters required by fractal cal-culation 接触面材料弹性模量E/GPa泊松比电阻率/(m)屈服强度/MPa分形维数D特征尺度G/m名义接触面积Aa/m2铍铜1280.34 1.711081401.51107211.3107铝合金720.35 2.73108310银730.38 1.5910835 024681002004006008001 000单片触指接触电阻Rc/单片触指所受载荷/N分形理论计算结果试验测量结果024681010681216202414182226误差/%单片触指所受载荷/N(b)分形理论和试验测量结果误差(a)计算值和试验值对比图 6单片触指ECR分形计算和试验测量结果Fig.6Single-piece ECR fractal calculation and test meas-urement results 第 4 期杨为，等：GIS设备用表带触指电连接结构接触电阻的分形理论计算与分析17随着载荷从4 N增大至10 N，由于试验过程中表带触指发生了偏转，导致接触面积增大，ECR减小，因此相对误差略微增大，在10%20%之间波动。综上可知，接触电阻试验结果验证了分形理论模型数值计算单片表带触指不同压力下接触电阻的合理性，可为触指电连接结构的电磁热流体场的仿真计算提供理论基础。3.2 分形维数对实际接触面积的影响从式（30）（33）可以看出，接触电阻与实际接触面积Ar*有关，Ar*=Ar/Aa。为了探究不同表带触指表面微观形貌、接触载荷对接触电阻和实际接触面积Ar*的影响规律，通过改变触指表面分形维数D和载荷范围，获得了弹性变形、第一弹塑性变形两个阶段电接触量纲为1的实际接触面积Ar*与载荷Fr*的关系，如图7所示。由图7（a）可以看出：当载荷较小时，触指表面微凸体只发生弹性变形，Fr*与Ar*成正比。由于弹性变形属于可恢复的初始阶段，微小的载荷变化将导致接触面积出现较大波动。不同分形维数下实际接触面积的变化趋势相近，因此随着D增大，实际接触面积指数级增大，数值差异较大。由图7（b）可以看出：分形维数D相同时，Fr*与Ar*之间近似呈线性关系，斜率有小幅度改变；随着分形维数D的增大，Fr*与 Ar*的直线坡度也随之增大，且增大幅度也随之增大。这是由于：D越大，粗糙表面微观形态越复杂，单位面积上的微凸体数量越多，在实际接触面积同样的情况下，粗糙面微凸体承载的接触压力减小，变形量减小，单个微凸体的承载力增加。图8为当分形维数D不同，仅改变载荷P的大小时，弹性变形接触面积在实际接触面积中的占比Are/Ar与量纲为1的实际接触面积Ar/Aa间的关系。由图（8）可知：在载荷逐渐增加，接触面积逐渐增大的过程中，触指电接触粗糙表面的接触从最初的弹性阶段过渡到弹塑性变形阶段。当实际接触面积Ar/Aa相同时，分形维数D越大，Are/Ar越大。这是因为分形维数D表示表面复杂程度，D越大则相同面积下微凸体越多，其抵抗变形的能力就越强，弹性面积在实际接触面积中的占比就越大。当分形维数D相同时，随着实际接触面积Ar/Aa增大，弹塑性由临界点开始进入，此时Are/Ar从第一个临界点开始减小。表2给出了不同分形维数下，电接触从弹性变形过渡到第一弹塑性变形和从第一弹塑性变形过渡到第二弹塑性变形阶段弹性面积占比变化时的实际接触面积。由于分形维数衡量的是粗糙表面的微观结构复杂程度24，由每一次弹性面积占比变化的临界点可以看出，触点的分形维数越大，变形越先进入弹塑性阶段，同理也容易进入第二弹塑性变形阶段。即：分形维数越大，在实际接触面积较小时，就从弹性变形过渡到第一弹塑性变形阶段，或是从第一弹塑性变形过渡到第二弹塑性变形阶段。024681016101410101061041021011081012101510111071051031091013Ar*Fr*/1017(a)弹性变形阶段D=1.3D=1.5D=1.7D=1.9051020150468210Ar*/103 Fr*/105(b)第一弹塑性变形阶段D=1.3D=1.5D=1.7D=1.9图 7量纲为1的载荷与实际接触面积的关系Fig.7Relationship between the load and the actual con-tact area with dimension one 10141012101010810610410200.40.60.80.21.0AreAr1ArAa1D=1.9D=1.7D=1.5D=1.3图 8弹性面积占比与实际接触面积的关系Fig.8Relationship between the percentage of el