基于最小热阻力法则的油浸式变压器绕组等效导热系数计算方法及应用.pdf

第 50 卷第 4 期2023 年 7 月华 北 电 力 大 学 学 报Journal of North China Electric Power UniversityVol.50,No.4Jul.,2023doi:10.3969/j.ISSN.1007-2691.2023.04.07基于最小热阻力法则的油浸式变压器绕组等效导热系数计算方法及应用沈 楠1,刘 刚1,王晓晗2,武卫革3(1.河北省输变电设备安全防御重点实验室(华北电力大学),河北 保定 071003;2.国网保定电力公司,河北 保定 071003;3.保定天威保变电气股份有限公司 河北省输变电装备电磁与结构性能重点实验室,河北 保定 071056)摘要:准确计算变压器绕组温升及热点需要考虑绕组匝间绝缘的影响,而考虑匝间绝缘的精细建模的建模工作量和剖分网格量都非常大。针对此问题,提出了一种基于最小热阻力法则的等效导热系数计算方法,分别计算油浸式电力变压器绕组的轴向等效导热系数和径向等效导热系数,从而可以获得等效导热系数下的饼式绕组来模拟考虑匝间绝缘的实际绕组结构。为了验证所提方法的有效性,建立了油浸式变压器绕组温升仿真模型,将等效模型的仿真结果与考虑匝间绝缘的仿真结果进行温升和热点的对比分析。结果表明二者的温升和热点位置基本一致,且等效模型减少了 72%的网格剖分量,验证了所提方法的有效性。将所提方法应用于换流变压器网侧和阀侧绕组的温升和热点分析,得到了换流变压器的绕组温升和热点分布。关键词:油浸式变压器;等效导热系数;最小热阻力法则;热点中图分类号:TM401 文献标识码:A 文章编号:1007-2691(2023)04-0060-09收稿日期:2021-12-27.基金项目:国家重点研发计划项目(2021YFB2401703);河北省重点研发计划项目(20312101D).Calculation and Application of Equivalent Thermal Conductivity of Oil-immersed Transformer Windings Based on the Law of Least Thermal ResistanceSHEN Nan1,LIU Gang1,WANG Xiaohan2,WU Weige3(1.Hebei Provincial Key Laboratory of Power Transmission Equipment Security Defense(North China Electric Power University),Baoding 071003,China;2.State Grid Baoding Electric Power Supply Company,Baoding 071003,China;3.Hebei Key Laboratory of Electromagnetic&Structural Performance of Power Transmission and Transformation Equipment,Baoding Tianwei Baobian Electric Co.,Ltd,Baoding 071056,China)Abstract:Accurate calculation of transformer winding temperature rise and hot spot need to consider the influence of winding interturn insulation,but the workload of fine modeling and the amount of meshing considering winding interturn insulation are very large.To solve this problem,this paper presented a method based on the minimum thermal resist-ance rule to calculate the equivalent thermal conductivity of oil-immersed power transformer windings.The axial equiva-lent thermal conductivity and radial equivalent thermal conductivity of oil-immersed power transformer windings were calculated respectively.In this way,the cake winding with equivalent thermal conductivity can be obtained to simulate the actual winding structure considering the interturn insulation.In order to verify the effectiveness of the proposed method,a simulation model of the temperature rise of the oil-immersed transformer windings was established,and the simulation results on temperature rise and hot spot of equivalent model were compared with the simulation results of ac-tual model considering interturn insulation.The results show that the temperature rise and hot spot location are basically the same,and the equivalent model reduces the mesh profile component by 72%,which verifies the effectiveness of the proposed method.The proposed method is applied to the temperature rise and hot spot analysis of the grid side and 第 4 期沈 楠,等:基于最小热阻力法则的油浸式变压器绕组等效导热系数计算方法及应用valve side windings of the converter transformer,and the winding temperature rise and hot spot distribution of the con-verter transformer are obtained.Key words:oil-immersed transformer;equivalent thermal conductivity;minimum thermal resistance rule;hot spot0 引 言 电力变压器是电力系统中的关键设备,电力变压器的运行可靠性直接影响电力系统运行的稳定性和可靠性1,2。电力变压器绕组温度较铁心和其它部件高,并且温度分布不均匀,局部过热导致绝缘老化,绝缘性能下降,降低电力变压器的使用寿命,遵从“六度法则”。因此变压器绕组热点温度是影响变压器使用寿命的决定性因素,也是影响变压器负载能力的最主要限制因素。准确计算电力变压器绕组热点对于预防电力变压器故障,预测电力变压器使用寿命,电力变压器的设计和研制等都有重要意义3-5。电力变压器绕组热点温度可以通过测量法6得到,测量法通过将光纤探头或者热电偶放置于电力变压器绕组的热点区域,测量热点温度,这种方法需要在变压器制造过程中将温度传感器放入预估的热点位置中,对于已经投运的变压器温升无法采用此方法。其余获得热点温度的常见方法的有热路模型法7-11、经验公式法12,13、人工智能算法14,15和多物理场建模计算方法16-18。热路模型法大多基于传热学原理,通过底层油温、顶层油温等通过热电类比法将变压器绕组等效为热路模型7-11,根据公式求得绕组热点温度和绕组热点位置,该方法受模型参数影响较大,而有些输入参数需要通过大量实验测量;经验公式法基于经验公式可以粗略推算特定变压器在特定运行情况下的绕组热点温度12,13,准确性不高,适用性差,受变压器运行状况影响较大;人工智能算法通过神经网络算法等进行变压器热点温度反演和预测14,15,可以考虑多种因素影响,有较好的适用性,但是该方法对初始模型参数和数据有较大的依赖性;多物理场是采用电磁、流、热耦合方法进行耦合场分析16-18,其中温升及热点分析时,往往是采用有限元或者迎风有限元进行分析,但是在温度场分析时,由于绕组间存在匝间绝缘,而匝间绝缘的厚度只有 mm 级,因此建模时需要大量的时间,同时采用有限元分析时,由于绕组和匝间绝缘尺寸相差很大,导致网格的剖分量很大,从而导致计算量大大增加。等效导热系数表征复合材料导热能力的特征参数,目前可根据经验公式法、最小热阻法、热阻网络法、数值模拟法、渐进均匀化法等方法对复合材料的等效导热系数进行预测19,20。陈则韶21通过最小热阻法以及比等效导热系数法则计算与复合材料拥有相同比等效热阻的单元体等效导热系数,该单元体等效导热系数与复合材料的等效导热系数相同。李垣明等22通过均匀化理论以及ABAQUS 仿真软件,建立了四种不同相体积下的燃料等效导热系数,通过该方法初步建立了燃料等效导热系数计算模型,但该模型考虑因素不充分,需要进一步研究。王亮亮23根据逾渗理论,建立了复合材料的等效逾渗导热方程,提出了导热聚合物两相体系的“海岛-网络”模型,该方法得到的高、低填充型复合材料的等效导热系数符合实际。目前在进行绕组温升及热点仿真时,为了减少建模工作量,常常将分匝的绕组简化为线饼24,25,但是这种方法忽略了绕组间的匝间绝缘,绕组温升和热点计算结果与实际情况会有很大偏差。因此,为了更加准确的计算变压器绕组温升及热点,同时兼顾绕组建模的工作量和仿真时的计算量,本文尝试基于最小热阻力法则21,26,27,将绕组等效为热路模型,分别计算横向热阻无穷大、横向热阻无穷小时的等效导热系数,然后取两者平均值,求得绕组轴向等效导热系数和径向等效导热系数。为验证所提方法的有效性,本文以大型油浸式变压器绕组二维模型为例,分别采用等效方法与精细分匝模型进行绕组温升和热点仿真,并对绕组温升和热点的仿真结果进行了对比分析。在有效性验证的基础上,将所提方法应用于换流变压器网侧和阀侧绕组的温升和热点分析。1 基于最小热阻力法则的等效导热系数计算方法 油浸式电力变压器绕组模型如图 128所示,变压器单个绕组线饼模型切面图如图 2 所示。像电流沿着电阻力最小的电路流动,同样,在16华 北 电 力 大 学 学 报2023 年图 1 变压器绕组模型Fig.1 Model of transformer winding图 2 变压器绕组线饼Fig.2 The disc of transformer winding传递热量的过程中,热量同样会沿着热阻力最小的热路流动,即最小热阻力法则,相应的热路为最小热阻。Tr=Rq(1)对于并联的热路通道,热阻力是相同的,此时热流最大,热阻最小,为等效热阻值。其中 r为热阻力,R 为热阻,q 是单位时间内传导的热量,也称导热速率。根据傅里叶公式,对于均匀介质:R=L/(A)(2)其中,L 为热阻通道长度,A 为导热面积;为导热系数。对于油浸式电力变压器绕组这样的复合材料,则有公式Re=L/(Ae)(3)图 3 为单匝导线等效热阻图,图中 R1代表绝缘层 1 的热阻,R5代表铜线匝的热阻,R代表绝缘层微元块和铜线匝微元块接点之间的热阻,即横向热阻。将单匝导线分割成为 9 个微元块,根据公式(1)计算微元块的平均热阻,借用电路原理的原则和方法,将热阻比作电阻,应用电路串并联等效原理,求得等效热阻。同时,为了简化计算,将同一热流方向上的具有相同导热系数的材料等效为一个通道。根据数学逼近法,可以做出两种极限假设即横向热阻无穷大和横向热阻无穷小,轴向对应两个热路图。图 4 为横向热阻无穷大时的等效热阻图。可知横向热阻无穷大时,根据戴维南等效定则,轴向等效热阻是用先计算各通道串联等效热图 3 单匝导线等效热阻图Fig.3 Single zone equivalent thermal resistance图 4 横向热阻无穷大Fig.4 Infinite transverse thermal resistance阻,然后计算各通道并联等效热阻的计算方法得到的。横向热阻无穷大时,轴向等效热阻 Rvt-b可利用公式(1)求得:Rvt-b=1ni=11Ri=1ni=1nj=1Aiij/Lij(4)其中,i 代表热流通道序号,j 代表微元所处分层的行号;Lij为第 i 通道第 j 层微元的长度;Ai为第 i通道的横截面积。图 5 为横向热阻无穷小时的等效热阻图。可知当横向热阻无穷小时,根据戴维南等效定则,轴向等效热阻是由先计算各通道并联等效热阻,然后计算各通道串联等效热阻的计算方法得到的。横向热阻无穷小时,轴向等效热阻 Rvb-t可利用公式(2)求得:Rvb-t=ni=1Ri=ni=1nj=1Lij/Aij(5)则轴向等效热阻 Rv可表示为26第 4 期沈 楠,等:基于最小热阻力法则的油浸式变压器绕组等效导热系数计算方法及应用图 5 横向热阻无穷小Fig.5 Infinitesimal transverse thermal resistanceRv=Rvt-b+Rvb-t2(6)轴向等效导热系数 v可表示为v=2vt-bvb-tvt-b+vb-t=vt-b+vb-t2(7)式中:vt-b为横向热阻无穷大时求得的轴向导热系数;vb-t为横向热阻无穷小时求得的轴向导热系数。变压器径向同样对应两个等效热路图,即横向热阻无穷大和横向热阻无穷小,如图 6、图 7 所示。图 6 横向热阻无穷大Fig.6 Infinite transverse thermal resistance图 7 横向热阻无穷小Fig.7 Infinitesimal transverse thermal resistance由轴向公式,推广到径向,可得当横向热阻无穷小时,径向等效热阻 Rht-b可表示为Rht-b=1ni=11Ri=1ni=1nj=1Lij/Aij(8)当横向热阻无穷大时,径向等效热阻 Rhb-t可表示为Rhb-t=ni=1Ri=ni=1nj=1Lij/Aij(9)则径向等效热阻 Rh可表示为Rh=Rht-b+Rhb-t2(10)径向等效导热系数 h可表示为h=2ht-bhb-tht-b+hb-t=ht-b+hb-t2(11)式中:ht-b为横向热阻无穷大时求得的径向导热系数;hb-t为横向热阻无穷小时求得的径向导热系数。2 油浸式电力变压器绕组模型 采用文献29中的油浸式电力变压器绕组模型,绕组模型如图 8 所示。对导向分区从上到下依次进行编号,共有 8 个导向分区;对线饼从上到下依次进行编号,共有 66 个线饼,其中前 3 个分区包含 7 个线饼,后 5 个分区包含 9 个线饼;从内轴到外轴依次进行编号,每个线饼共包含 30 根导线。图 8 变压器绕组模型Fig.8 Model of transformer winding为验证等效模型的有效性,本文采用考虑变压器绕组匝间绝缘的分匝模型作为对照组,其剖分网格单元数为 883 989;采用基于最小热阻力法则的等效导热系数方法作为实验组,剖分网格单元数为 211 613。模型的物性参数如表 1 所示。使用基于最小热阻力法则的等效导热系数计算方法在仿真时不考虑匝间绝缘,将分匝的绕组等效为线饼,线饼的材料属性设置为铜,同时将线饼的导热系数设置36华 北 电 力 大 学 学 报2023 年为通过最小热阻力法则计算所得的轴向等效导热系数和径向等效导热系数,通过计算得到的轴向等效导热系数为1.642 W(mK)-1、径向等效导热系数为 0.781 3 W(mK)-1。该模型的入口设置为速度和温度边界条件,速度方向为竖直方向,大小为 0.02 m/s,对两种模型施加相同的热源,温度设置为 330 K。出口设置为压力边界条件,热流密度为 0。绕组、绝缘筒设置为无滑移的壁面边界条件和绝热边界条件;绕组的边界为流-固耦合的无滑移壁面边界条件。表 1 材料物性参数表Tab.1 Physical parameters of model materials材料参数类型数值变压器油密度/kgm-31 098.72-0.712T定压比热容/J(kgK)-1807.163+3.58 T导热系数/W(mK)-10.1 509-7.10110-5 T动力粘度/Pas0.0 846-410-4 T+510-7 T2绝缘纸密度/kgm-3980定压比热容/J(kgK)-12 000导热系数/W(mK)-10.25铜密度/kgm-38 960定压比热容/J(kgK)-1390导热系数/W(mK)-13383 有效性验证 考虑匝间绝缘和基于最小热阻力法则的等效模型的温度场仿真结果如图 9 所示,从左往右依次是分匝模型、等效模型的温度场仿真结果图。图 9 温度场云图Fig.9 Temperature distribution contour由温度场云图可以看出,两种方法得到的仿真结果基本一致,差异并不明显。提取模型的轴向中心线的仿真结果如图 10 所示。在轴向垂直中心线方向,热点位置一致,分匝模型与等效模型最大温度相对误差百分比的绝对值为 1.54%,对应温度误差绝对值为 5.67 K。本文提出的等效导热系数方法能很好近似等效前的仿真结果,仿真得到的等效模型的几个线饼误差如下:线饼 12 的最大相对误差百分比为0.57%,对应的相对误差绝对值为 2 K;线饼 20 的最大相对误差百分比为 0.53%,对应的相对误差绝对值为 1.85 K;线饼 30 的最大相对误差百分比为 0.51%,对应的相对误差绝对值为 1.74 K;线饼38 的最大相对误差百分比为 0.5%,对应的相对误差绝对值为 1.72 K。等效模型的温度场仿真结果在径向上与分匝模型的最大误差发生在线饼 12,线饼 12 的径向温度对比如图 11 所示。图 10 轴向温度结果对比图Fig.10 Comparison chart of axial temperature results图 11 线饼 12 径向温度对比图Fig.11 Comparison chart of line pie 12 radial temperature等效方法的热点位置与分匝模型的热点位置46第 4 期沈 楠,等:基于最小热阻力法则的油浸式变压器绕组等效导热系数计算方法及应用基本一致;表 2 给出了分匝模型与等效模型下几个线饼的热点温度,充分验证了本文所提等效导热系数方法的有效性和准确性。表 2 各线饼的热点温度Tab.2 Hot spot temperature of disc线饼 12线饼 20线饼 30线饼 38分匝模型364.24 K364.95 K371.87 K362.61 K等效模型364.23 K365.07 K371.84 K362.64 K4 应 用 在有效性验证的基础上,将本文方法应用于实际换流变压器的绕组温升和热点分析。建立换流变压器仿真模型,如图 12 所示,由于换流变压器基本符合轴对称结构,因此在图中给出部分简化后的变压器几何结构图。图 12 简化后部分变压器几何模型图Fig.12 Geometric model diagram of simplified parts of transformer换流变压器铁芯、绝缘纸板的物性参数如表 3所示,通过计算得到的网侧绕组和阀侧绕组的轴向等效导热系数为 2.588 W(mK)-1,径向导热系数为 0.965 5 W(mK)-1,将入口速度设置为 0.02 m/s,出口设为压力出口边界条件。网侧绕组和阀侧绕组施加的热源参考文献30,网侧为 67 000 Wm-3,阀侧为 37 500 Wm-3。表 3 材料的物性参数Tab.3 Physical parameters of model materials材料参数类型数值铁芯密度/kgm-37 650定压比热容/J(kgK)-1460导热系数/W(mK)-145绝缘纸板密度/kgm-31 200定压比热容/J(kgK)-12 000导热系数/W(mK)-10.03采用等效模型的仿真结果如图 13 所示。对换流变压器模型的网侧绕组和阀侧绕组从上到下依次进行编号,提取各线饼中心的温度场数据进行对比,得到图 14 的结果。从图 14 中可以看出,等效模型下的换流变压器网侧绕组的温度明显高于阀侧绕组的温度,这主要是由于网侧绕组的热源密度高,导致网侧的整体温度高。网侧绕组和阀侧绕组各分区的最高温度都出现在每分区的最后一个线饼中,这是因为挡板的存在阻碍了变压器油的流动,最后一个线饼和挡板之间的油流速度较小,阻碍了线饼中热量的散出,温度较高。图 13 等效模型温度场云图Fig.13 Temperature distribution contour图 14 绕组轴向温度对比图Fig.14 Comparison chart of axial temperature of winding文献7中基于绕组热分布的油浸式变压器绕组热点温度改进计算模型,绕组的热点位置大约出现在绕组高度的 88%92%,由图 14 可以看出,本文换流变压器的热点温度出现在第一分区的最后一个线饼中,位于整个绕组高度的 92%,符合文献中的结论。5 结 论 论文针对准确计算变压器绕组温升及热点时建模工作量大、计算代价高的问题,基于传热学基本原理给出了基于最小热阻力法则的等效导热系56华 北 电 力 大 学 学 报2023 年数计算方法,并基于大型变压器绕组模型验证了温升及热点计算的有效性,计算结果表明:(1)采用考虑匝间绝缘的精细建模时剖分网格数量是 883 989,而采用等效导热系数方法建模时,剖分网格数量为 211 613,网格数量减少了76%。(2)综合对比等效模型和分匝模型的温度场仿真结果,热点位置基本一致,各个线饼间的温度分布和热点温度基本一致,验证了等效导热系数方法的准确性。(3)将所提方法应用于换流变压器绕组温升及热点分析中,得到了换流变压器绕组的温度场分布。仿真结果表明网侧绕组和阀侧绕组中的热点位置均位于绕组高度的 92%,与基于绕组热分布的改进热点温度计算模型计算结果一致。参考文献:1 陈伟根,奚红娟,苏小平,等.广义回归神经网络在变压器绕组热点温度预测中的应用 J.高电压技术,2012,38(1):16-21.CHEN Weigen,XI Hongjuan,SU Xiaoping,et al.Ap-plication of generalized regression neural network to transformer winding hot spot temperature forecasting J.High Voltage Engineering,2012,38(1):16-21(in Chinese).2 张明泽,刘骥,齐朋帅,等.变压器油纸绝缘老化含水率频域法定量评估 J.高电压技术,2019,45(10):3221-3230.ZHANG Mingze,LIU Ji,QI Pengshuai,et al.Quanti-tative assessment of moisture content in frequency do-main response for transformers oil-paper insulation aging J.High Voltage Engineering,2019,45(10):3221-3230(in Chinese).3 周利军,王健,王路伽,等.强油导向结构变压器绕组区域温度建模及热点定位 J.高电压技术,2020,46(11):3896-3904.ZHOU Lijun,WANG Jian,WANG Lujia,et al.Ther-mal modeling and hot spot locating for transformer wind-ing in oil forced and directed cooling mode J.High Voltage Engineering,2020,46(11):3896-3904(in Chinese).4 刘云鹏,李昕烨,李欢,等.内置分布式光纤传感的 35 kV 油浸式变压器研制 J.高电压技术,2020,46(6):1886-1894.LIU Yunpeng,LI Xinye,LI Huan,et al.Development of 35kV oil-immersed transformer with built-in distribu-ted optical fiber J.High Voltage Engineering,2020,46(6):1886-1894(in Chinese).5 苏小平,陈伟根,奚红娟,等.采用 Kalman 滤波算法预测变压器绕组热点温度 J.高电压技术,2012,38(8):1909-1916.SU Xiaoping,CHEN Weigen,XI Hongjuan,et al.Winding hot-spot temperature prediction for oil-im-mersed power transformer based on Kalman filter algo-rithm J.High Voltage Engineering,2012,38(8):1909-1916(in Chinese).6 徐肖伟,付骁,董涛,等.油浸式配电变压器不平衡运行条件下的温升实验研究 J.云南电力技术,2020,48(6):81-85.XU Xiaowei,FU Xiao,DONG Tao,et al.Study on temperature-rising experiment of oil-immersed distribu-tion transformer in unbalanced operation J.Yunnan Electric Power,2020,48(6):81-85(in Chi-nese).7 王丰华,周翔,高沛,等.基于绕组热分布的改进油浸式变压器绕组热点温度计算模型 J.高电压技术,2015,41(3):895-901.WANG Fenghua,ZHOU Xiang,GAO Pei,et al.Im-proved thermal circuit model of hot spot temperature in oil-immersed transformers based on heat distribution of winding J.High Voltage Engineering,2015,41(3):895-901(in Chinese).8 罗汉武,陈连凯,姜国义,等.计及环境条件的电力变压器热路模型及其应用 J.高电压技术,2018,4(11):3561-3568.LUO Hanwu,CHEN Liankai,JIANG Guoyi,et al.Thermal circuit model of the transformer considering en-vironment factors and its application J.High Voltage Engineering,2018,4(11):3561-3568(in Chi-nese).9 李永建,闫鑫笑,张长庚,等.基于磁-热-流耦合模型的变压器损耗计算和热点预测 J.电工技术学报,2020,35(21):4483-4491.LI Yongjian,YAN Xinxiao,ZHANG Changgeng,et al.Numerical prediction of losses and local overheating in transformer windings based on magnetic-thermal-fluid model J.Transactions of China Electrotechnical So-ciety,2020,35(21):4483-4491(in Chinese).10 魏本刚,徐延顺,黄华,等.分体式变压器顶层油温升计算模型 J.变压器,2019,56(1):57-62.WEI Bengang,XU Yanshun,HUANG Hua,et al.Top oil temperature calculation model of split-type trans-former based on thermal circuit and krill herd algorithm 66第 4 期沈 楠,等:基于最小热阻力法则的油浸式变压器绕组等效导热系数计算方法及应用J.Transformer,2019,56(1):57-62(in Chi-nese).11 丁玉琴,张乔根,高萌,等.油浸式配电变压器分布式热路模型 J.高电压技术,2019,45(3):968-974.DING Yuqin,ZHANG Qiaogen,GAO Meng,et al.Distributed thermal circuit model of oil-immersed distri-bution transformers J.High Voltage Engineering,2019,45(3):968-974(in Chinese).12 NARBUTOVSKIH P.Simplified graphical method of computing thermal transients J.Transactions of the American Institute of Electrical Engineers,1947,66(1):78-83.13 BLAKE J,KELLY E.Oilimmersed power transformer overload calculations by computer J.IEEE Transac-tions on Power Apparatus and Systems,1969,PAS-88(8):1205-1215.14 骆小满,阮江军,邓永清,等.基于多物理场计算和模糊神经网络算法的变压器热点温度反演 J.高电压技术,2020,46(3):860-866.LUO Xiaoman,RUAN Jiangjun,DENG Yongqing,et al.Transformer hot-spot temperature inversion based on multi-physics calculation and fuzzy neural network algo-rithm J.High Voltage Engineering,2020,46(3):860-866(in Chinese).15 HELL M,COSTA J,GOMIDE F.Recurrent neuro fuzzy net work in thermal modeling of power transform-ers J.IEEE Transactions on Power Delivery,2007,22(2):904-910.16 王永强,马伦,律方成,等.基于有限差分和有限体积法相结合的油浸式变压器三维温度场计算 J.高电压技术,2014,40(10):3179-3185.WANG Yongqiang,MA Lun,L Fangcheng,et al.Calculation of 3D temperature field of oil immersed transformer by the combination of the finite element and finite volume method J.High Voltage Engineering,2014,40(10):3179-3185(in Chinese).17 苏小平,陈伟根,胡启元,等.基于解析-数值技术的变压器绕组温度分布计算 J.高电压技术,2014,40(10):3164-3170.SU Xiaoping,CHEN Weigen,HU Qiyuan,et al.Cal-culation for transformer winding temperature distribution by numerical analytical technology J.High Voltage Engineering,2014,40(10):3164-3170(in Chi-nese).18 邓永清,牛曙,晋涛,等.66kV 油浸式变压器瞬态热点温度计算 J.变压器,2020,57(7):7-11.DENG Yongqing,NIU Shu,JIN Tao,et al.Calcula-tion of transient hot spot temperature of 66kV oil-im-mersed transformer J.Transformer,2020,57(7):7-11(in Chinese).19 钱立波,余红星,孙玉发,等.固-固二元复合材料等效导热系数模型研究综述及评价 J.原子能科学技术,2020,54(3):409-420.QIAN Libo,YU Hongxing,SUN Yufa,et al.Review and evaluation on equivalent thermal conductivity model of solid-solid binary composites J.Atomic Energy Science and Technology,2020,54(3):409-420(in Chinese).20 董其伍,刘琳琳,刘敏珊.预测聚合物基复合材料导热系数方法综述 J.化学推进剂与高分子材料,2007,(6):36-39.DONG Qiwu,LIU Linlin,LIU Minshan.Review of methods for predicting thermal conductivity of polymer matrix composites J.Chemical Propellants&Poly-meric Materials,2007,(6):36-39(in Chinese).21 陈则韶,钱军,叶一火.复合材料等效导热系数的理论推算 J.中国科学技术大学学报,1992,(4):416-424.CHEN Zeshao,QIAN Jun,YE Yihuo.Theoretical cal-culation of equivalent thermal conductivity of composite materials J.Journal of University of Science and Technology of China,1992,(4):416-424(in Chi-nese).22 李垣明,唐昌兵,余红星,等.锆基弥散微封装燃料等效传热系数数值模拟研究 J.核动力工程,2018,39(2):76-79.LI Yuanming,TANG Ch