高功率密度SiC静止无功补偿器强迫风冷散热综合建模及优化设计方法

林弘毅 伍 梁 郭 潇 陈国柱

(浙江大学电气工程学院 杭州 310027)

摘要 为了提高大功率高功率密度电力电子装置热设计的准确性和设计效率,该文综合热传导、对流换热与流体力学理论,针对电力电子装置的典型强迫风冷散热系统,提出基于截面积二次方根为无量纲特征长度的综合热模型,并提出一种强迫风冷散热系统体积最优的优化设计方法。以380V/50kvar高功率密度静止无功补偿器(SVG)为例,比较所提出的综合热模型与温升实验的热阻,综合热模型平均热阻误差在7%以内,即在SVG满载运行时,散热器平均表面温升的绝对误差DTe在1.3℃以内。该文提出的综合热模型相对于传统热模型精度提升62%,证明了综合热模型的准确性。基于所提出的优化方法设计的散热系统体积4.03L,对比传统方法设计的体积5.7L,体积缩小30%,从而证明了该文提出优化设计方法的可行性。

关键词:热设计 热模型 强迫风冷 静止无功补偿器(SVG)

0 引言

电力电子器件在工作中会产生损耗,这些损耗会转化为热量。若热设计不合理,电力电子器件的结温过高,将导致电力电子器件的失效率增大,较高的过温还会造成器件烧毁,直接影响电力电子装置的寿命和可靠性。随着大功率电力电子装置向高功率密度发展,电力电子器件的散热问题越来越突出,进而影响了电力电子装置的可靠性和稳定性,成为电力电子装置功率密度进一步提高的瓶颈[1]

随着宽禁带器件的发展,电力电子装置的开关频率得以提升,无源器件的体积显著减小[2-3]。对基于宽禁带器件的强迫风冷电力电子装置而言,散热系统(包括散热器和风扇)占装置总体积的25%以上。因此,散热系统体积优化对提高电力电子装置的功率密度起着关键作用。

目前,电力电子装置的热设计主要依赖工程实践经验。有些学者基于实践经验提出了一些经验公式,但经验公式通常误差较大,且不具有普遍适用性[4]。文献[5]采用HyperMesh和Fluent软件联合仿真进行散热设计,但是热设计过程缺乏理论支撑。文献[6]基于传统热模型,提出了散热系统的优化设计方法,对15V/2kA的高频、大功率开关电源样机进行散热设计,温度误差在3.8%以内。但是其在优化过程中仅利用热仿真软件对散热器几何尺寸进行优化,由于散热器变量较多且互相耦合,设计的效率较低。文献[7]基于传统热模型对20kW电动汽车SiC逆变器的散热器进行优化设计,热阻误差在15%以内,但是仅对散热器几何尺寸进行优化,忽略了风扇的影响。文献[8]基于传统热模型,对强迫风冷散热系统进行建模分析,并提出了一种衡量散热效果的新指标,但是其模型误差较大。

在电力电子装置设计之初,散热设计应该和电路设计、结构设计同步规划开展。散热系统设计流程如图1所示。值得注意的是,热设计包括理论计算、热仿真验证、模拟实验验证等流程,耗时较长。因此,在热设计完成之前,可以先利用传统经验公式对散热进行概要设计,初步验证设计方案的可行性,并指导结构设计和功率PCB布局。在热设计的过程中,电路测试可以同步进行,以提升设计效率。

本文首先介绍了380V/50kvar静止无功补偿器(Static Var Generator, SVG)系统的电路拓扑和主要参数,并分析了SVG系统的损耗,介绍了典型强迫风冷散热系统的结构及其热阻网络;然后利用传统经验公式对散热系统进行概要设计;接着综合热传导、对流换热与流体力学理论,针对电力电子装置的典型强迫风冷散热系统,提出了一种基于截面积二次方根为无量纲特征长度的综合热模型。同时,提出一种优化方法对风扇和散热器总体积进行优化;最后给出了散热系统有限元仿真和温升实验的结果。本文提出的综合热模型平均热阻误差在7%以内,经过折算,在SVG满载运行时,散热器平均表面温升的绝对误差DTe在1.3℃以内。而传统热模型[7-8]平均热阻误差为69%,证明了本文提出散热系统综合热模型的准确性。同时,对比了传统设计方法和本文提出优化设计方法的散热系统体积,基于优化方法设计的散热系统体积为4.03L,较传统经验公式设计的体积5.7L,缩小30%,证明了本文提出优化设计方法的可行性。

width=137.5,height=239.5

图1 系统设计流程

Fig.1 The design process of power system

1 SiC-SVG损耗分析与散热器热阻网络

1.1 SiC-SVG电路设计与损耗分析

SVG系统的电路拓扑如图2所示。考虑成本,系统由两个模块并联而成,为了增加通流能力,每个开关由两个TO-247封装的SiC-MOSFET(型号SCT50N120)并联。系统的主要电路参数见表1。

MOSFET损耗由通态损耗PMOS,cond和开关损耗width=30,height=15组成。体二极管的损耗主要由通态损耗Pdiode,cond和反向恢复损耗Pdiode,rr组成。根据文献[9-10]和表1提供的电路参数可得到每个MOSFET的损耗分布,SVG满载时,平均每个MOSFET损耗Ps=36W,单个MOSFET损耗分布如图3所示。本文中共有24个MOSFET,总损耗功率P=864W。

width=220.8,height=160.7

图2 SVG电路拓扑

Fig.2 The structure of SVG

表1 SVG主要电路参数

Tab.1 Circuit parameters of SVG

参 数数 值 电网线电压Vd/V380 SVG直流侧电压Vdc/V750 SVG额定功率Sc/kvar50 最大单相输出电流有效值/A100 每个MOSFET输出电流有效值Ix/A25 开关频率fsw/kHz50 电网频率f0/Hz50 开关损耗常数项系数Esw,const/J17×10-6 开关损耗一次项系数Esw,k/J140×10-6 体二极管正向压降Vf/V4 二极管反向充电电荷Qrr/nC230 导通电阻Ron/mW60 死区时间Td/ns200

width=177.6,height=103.7

图3 单个MOSFET损耗分布

Fig.3 The losses distribution of single MOSFET

1.2 典型强迫风冷散热系统结构及其热阻网络

典型强迫风冷散热系统结构示意图如图4所示,典型强迫风冷散热系统结构由轴流风扇、系统风道和翅片型散热器组成。

width=204.5,height=164.05

图4 典型强迫风冷散热系统结构示意图

Fig.4 Typical structural of forced air cooling system

翅片型散热器的宽为b,长为L,基板厚为d,散热器基板的面积为Ahs=Lb,散热系统风道长度为Lduct。散热器有n个风道,共n+1个翅片,每个翅片的长为c、厚为t,相邻翅片间表面间距即散热器风道宽为s。散热系统需要k个风扇,其中,风扇厚为bfan,高为Hfan。考虑风扇后,散热系统的总体积可定义

width=121.95,height=17 (1)

散热系统热阻模型如图5所示,Tj为MOSFET的结温,Tc为MOSFET的壳温,Th为散热器的表面温度,Ta为环境温度,Rth,j-c为MOSFET结壳之间的热阻,Rth,TIM为MOSFET与散热器之间热界面材料(Thermal Interface Material, TIM)的热阻,Rth,h-a为散热器与空气之间的热阻[11]

width=190.3,height=112.9

图5 散热系统热阻模型

Fig.5 The thermal resistance of the cooling system

TIM的热阻Rth,TIM与材料的热导率、材料硬度、表面光洁度、安装压力等因素相关,很难精确计算。由于固体之间接触并非理想接触,接触面之间存在微小的空气空隙。因此,需要在固体导热垫片两面均匀涂上导热硅脂,填充空气间隙,如图5所示,以提高TIM与散热器的接触面积,提高散热效率。TO-247封装TIM的热阻值见表2。

表2 常用TIM及其热阻

Tab.2 Common TIM and its thermal resistance

TIM热导率/[W/(m·℃)]厚度(约)/mm热阻/(℃/W) 陶瓷(绝缘)301.50.11 优质硅胶(绝缘)1.60.30.39 普通硅胶(绝缘)0.90.30.70 云母(绝缘)0.710.30.88 导热硅脂1~50.50.21~1.05

装置的拓扑如图2所示,其对应的等效热阻网络如图6所示。热阻为在物体两端温度差与热源的功率之间的比值。基于等效热阻网络和热阻定义可得

width=96.95,height=65 (2)

width=127.8,height=131.15

图6 SVG散热系统等效热阻网络

Fig.6 The thermal resistance network of cooling system

根据热阻定义与式(2)可推导得到电力电子装置允许的最大热阻Rth,h-a,max

width=181,height=31 (3)

本文中设定器件的最大结温Tj,max=120℃,环境温度Ta,max=40℃。Rth,j-c=0.55℃/W,陶瓷垫片热阻Rth,TIM≈0.6℃/W。考虑散热器表面温度不均,温差为10℃,可求出满足电力电子装置允许的最大热阻Rth,h-a,max=0.030 6℃/W,即允许SVG满载时,散热器最大温升DT=26.4℃,允许MOSFET壳最大温升DTc=48℃。

2 强迫风冷散热系统概要设计

本节利用传统经验设计方法对强迫风冷散热系统进行概要设计,经验设计方法的公式简单,但是准确性较差,容易导致裕度过大。概要设计的目的有:①论证项目的合理性,如装置体积、功率密度是否可以基本达到指标要求;②热设计过程较为繁琐,需要理论设计、热仿真、模拟实验等步骤,在完成完整的热设计之前,概要设计结果可以指导装置结构设计和PCB布局,提高设计效率。

2.1 风机概要选型

根据能量守恒原则,系统稳定时产生的热量与其散发的热量是相同的,由此可以得到流量V0、系统产生的热量Q与空气温升DTair的关系为

width=74,height=30 (4)

式中,V0为实际通过散热器风道空气的流量,即散热系统工作点;rair为空气的密度,rair=1.23kg/m3Cair为空气的比热容,Cair=1 005J/(kg·℃);DTair为散热器出风口处的流体相对入风口的温度差。一般认为DTair<5℃。

由于风道并非理想封闭,散热器对流体存在阻力,风机提供的最大流量将会大于实际通过散热器风道空气流量V0。经验上认为风机提供的最大流量Vfan,max≈(1.5~2)V0

根据式(4)可知,本文SVG装置散热系统最小流量Vmin=0.139m3/s。本文初选风扇型号为AGB08038_ 24H,数量5个,高度Hfan=80mm。可提供最大风量为5×0.04m3/s =0.2m3/s,风扇详细参数见表3。

表3 备选风扇型号参数

Tab.3 The fan list and parameters

编号风扇型号Hfan/mmbfan/mm最大压强/Pa最大流量/ (m3/s) FAN1AGB08038_24H8038221.50.039 8 FAN2AGE08025_12U8025101.90.030 6 FAN3AGE08020_24U802070.560.024 1 FAN4AGE07025_24U7025125.40.025 7 FAN5AGB06038_24M6038515.50.031 6 FAN6AGE06025_12U6025127.40.016 4

2.2 散热器概要设计

散热器的散热面积Ahs须略大于功率半导体器件的总散热面积。翅片厚t和散热器风道宽s可根据经验在2~5mm之间选择。基板厚d对系统总热阻影响不大,可以在满足安装的情况下取一个较小的值。为了结构设计方便和提高散热系统体积利用率,可令kHfanbHfan=c+d。风道长度为Lduct可取30~50mm。散热器的热阻主要取决于散热器风道流体换热的热阻Rth,A估算,即

width=60,height=28 (5)

式中,h为流体传热系数,在强迫风冷的条件下可取60~80。功率PCB上的功率器件分布示意图如图7所示。散热面积长80mm,宽为350mm。因此,本文取散热器宽b=400mm,翅片长c=60mm,基板厚度d=20mm,翅片厚和散热器风道宽可令其分别为t =3mm、s=2mm,可得翅片数为n=b/(s+t)=80。风道长度可取Lduct=40mm。令Rth,A=0.030 6℃/W,取h=70,根据式(5)估算得到散热器的长L≈100mm。根据式(1)可知,概要设计散热系统总体积为5.70L。基于经验公式概要设计的散热系统参数见表4。

width=217.7,height=94.95

图7 功率PCB布局

Fig.7 Power PCB layout

表4 概要设计与优化设计散热器参数

Tab.4 The cooling system parameters

散热器参数数值(型号) 概要设计优化设计 长L/mm10090 宽b/mm400400 基板厚度d/mm2020 翅片长c/mm6040 翅片宽t/mm32 翅片个数n8080 风扇型号AGB08038_24HAGB06038_24M 体积/L5.704.03 综合热模型热阻Rth,h-a/(℃/W)0.020 50.023 8 传统热模型[7-8]热阻Rth,h-a/(℃/W)0.036 60.039 4

3 强迫风冷散热系统建模及优化设计

强迫风冷散热系统建模中,关键要得到表面传热系数h。表面传热系数h是一个过程量,与流体状态、流速、表面摩擦因数等物理量紧密相关。工程热物理学中,很难直接从理论推导出其表达式,通常会采用实验的方法总结归纳得出表达式,这种表达式又称之为实验关联式。

传统的实验关联式通常以水力直径作为无量纲特征长度,其适用的范围较广,但所造成的误差/不确定度较大,常常达到±25%[7-8]。文献[7-8]提出基于传统充分发展层流和湍流的实验关联式建立的传统热模型,所得到的平均热阻误差为69%。对电力电子装置的强迫风冷散热器而言,传统热模型误差偏大,准确性较低。

Yovanovich的物理实验结果表明,典型强迫风冷散热器风道窄而短,流体常常不能够充分发展,因此流体的状态通常为入口段层流[12-14]。由此,Yovanovich提出了适用于入口段层流的基于截面积二次方根为无量纲特征长度的实验关联式,误差/不确定度减小至±10%[13-14]。针对传统热模型误差偏大的情况,本文基于截面积二次方根为无量纲特征长度的实验关联式,综合热传导、对流换热与流体力学等理论进行建模,得到散热器热阻的平均热阻误差为7%。因此,本文提出的热模型更为准确,更适用于描述电力电子强迫风冷散热器。

3.1 散热器热阻模型

热设计最为关键的就是求出散热器配合风扇的热阻Rth,h-a。散热器单个风道热阻模型如图8所示。散热器热阻Rth,h-a由基板热传导热阻Rth,d、流体传热热阻Rth,conv和流体略过散热器风道引起流体温度变化的热阻Rth,fluid组成,即

width=130,height=16 (6)

width=92.3,height=120.5

图8 散热器风道热阻模型

Fig.8 Thermal resistance model of a heat sink air duct

基板热传导热阻Rth,d

width=60.95,height=30 (7)

式中,lhs为铝材料散热器的热导率,lhs=210W/(m∙℃)。根据风道热阻模型可以列出散热器风道的流体传热热阻Rth,conv表达式为

width=155,height=29 (8)

式中,Rth,FIN为散热器翅片传导热阻;Rth,aRth,A分别为散热器风道与流体换流的热阻。文献[15]提出了一种简化热阻模型,同时考虑了流体传热热阻Rth,conv、流量热阻Rth,fluid表面传热的有效面积Aeff、散热器效率h 以及表面传热系数h,分别表示为

width=186.95,height=63 (9)

width=77,height=15 (10)

width=99,height=64 (11)

width=67,height=87 (12)

式中,lair为空气的热导率,lair=0.03W/(m∙℃);V为通过散热器风道的空气流量;dh为散热器单个风道的水力直径;width=26,height=17为基于截面积二次方根为无量纲特征长度的努赛尔(Nusselt, Nu)数,表示为对流传热和传导散热的比值[12-14],可表示为

width=184.85,height=148.85(13)

式中,width=24.95,height=17为基于截面积二次方根为无量纲特征长度的雷诺(Reynolds, Re)数,用于表示流体的流动状态(湍流或层流);普朗特(Prandtl, Pr)为联系气体流体黏度和热导率的无量纲数[13],空气的Pr≈0.7。

width=26,height=17与流体的Pr具有强函数关联。因此,可表示为基于Pr的特征系数f(Pr)实验关联式[13],即

width=142,height=44 (14)

在假设散热器表面温度均匀的情况下[13],可以得出各项系数的值为

width=93,height=132.95 (15)

无量纲热充分发展长度z*和散热器风道结构特征参数width=8,height=9分别表示为

width=55,height=55 (16)

式中,nair为空气的运动黏度,nair=2.1×10-5m2/s。传热学中通常将摩擦因子fRe的相乘作为一项实验关联式[14]表示为

width=159,height=42.95 (17)

width=155,height=35 (18)

式中,width=33,height=18为基于截面积二次方根为无量纲特征长度的气体充分流动情况下的雷诺数。综上可知,在散热器结构参数和工作点V0已知的情况下,根据式(6)~式(18)可计算出散热器的热阻Rth,h-a

3.2 散热器与风扇流体力学模型

散热系统工作点V0通过风扇的静压降特征式Dpfan,k(V)和静压阻抗特征式Dptot(V)方程联立求解得

width=102,height=16 (19)

风扇并联的特点为:在相同的压强下,并联风扇的流量是单个风扇流量的k倍。据此,可以得到多个风扇并联系统的静压降特征曲线。利用多阶多项式拟合该曲线,得到多个风扇并联系统的静压降特征式Dpfan,k(V)为

width=207,height=18(20)

式中,aii=0, 1, 2, 3, …)为拟合多项式的系数。为了得到散热系统工作点V0,还需建立散热器的流体力学的模型,以求出散热器的静压阻抗特征式Dptot(V)。Dptot(V)主要由两部分构成:翅片散热器的静压阻抗Dphs(V)和系统风道的静压阻抗Dpduct(V)[14],分别表示为

width=128,height=15 (21)

翅片散热器的静压阻抗Dphs(V)、散热器出口处的截面积突增系数Kse、截面积散热器入口处的截面积突减系数Ksc、空气与散热器风道表面的摩擦因数fapp(V)以及系统风道的等效水力半径dh,duct分别表

[9, 14]

width=186.95,height=34 (22)

width=173,height=39 (23)

width=204,height=37(24)

width=107,height=30 (25)

width=49.95,height=30 (26)

系统风道的静压阻抗Dpduct(V)、空气与系统风道表面的摩擦因数fapp,duct(V)分别表示[9, 14]

width=216,height=35(27)

width=204.95,height=40(28)

width=178,height=44 (29)

式中,Kventuri为效应摩擦因数,Kventuri≈0.2。根据式(6)~式(18)可知散热器热阻Rth,h-a(V),根据式(19)~式(29)可求出散热系统的工作点V0,散热器工作点热阻的求解过程如图9所示。根据 图9,可知概要设计的散热系统热阻Rth,h-a= 0.020 5℃/W。采用文献[7-8]所提出的传统热模型,计算得到热阻Rth,h-a=0.036 6℃/W,如图9所示。

width=219.5,height=277.7

图9 散热器工作点热阻的求解过程

Fig.9 Trocess of solving the thermal resistance of heatsinlss at operating point

3.3 基于综合热模型的优化设计方法

散热器优化设计流程如图10所示。

(1)计算SVG满载情况下系统损耗Q,计算满足系统要求的最大热阻Rth,h-a,max。然后,根据整机结构要求、发热器件的分布、PCB布局,确定散热器结构参数取值范围。本文中,基板厚度d应尽可能设置为更小的值,本文取d=20mm。根据图7中发热器件的分布,取散热器宽400mm≤b≤450mm,长90mm≤L≤120mm。散热器翅片风道宽s,翅片厚度t的取值范围取决于散热器厂家的工艺,本文取2mm≤s≤20mm,2mm≤t≤20mm。

(2)对备选的风扇进行扫描,本文中备选风扇型号见表3。为了提高系统空间的利用率,令散热器高c+d等于风扇的高度Hfan,风扇的个数k=b/Hfan向下取整。接着,根据式(20)拟合得到风扇的静压降特征曲线Dpfan,k(V)。

(3)对散热器长L、宽b、风道宽s、翅片厚度t参数进行扫描,翅片参数的上、下限可分别表 示为

width=227.5,height=523.9

图10 散热器优化设计流程

Fig.10 The optimized design process of the cooling system

width=90,height=60.95 (30)

(4)根据式(19)~式(29)得到散热系统的工作点V0.

(5)根据式(6)~式(18)求出散热器配合风扇的热阻Rth,h-a

(6)当热阻Rth,h-a小于系统要求的最大热阻Rth,h-a,max时,输出散热系统参数、风扇型号以及散热系统的总体积Vh

(7)输出该风扇型号对应的符合要求的散热器的体积Vh,以及该体积Vh下的最小热阻Rth,h-a(即翅片厚s,风道宽t最优)。对所有参数和风扇型号扫描后,选择最优的散热器-风扇组合。

图11为风扇FAN1~FAN6对应的散热系统体积Vh与最小热阻Rth,h-a曲线,从图11可知,Rth,h-aRth,h-a,max的点均满足散热系统设计要求,其中优化后的散热器为A,其对应的参数见表4。与概要设计的散热器B相比,A的体积缩小了30%。

width=207.85,height=130.7

图11 散热系统体积Vh与最小热阻Rth,h-a曲线

Fig.11 Thermal resistance characteristics curves Rth,h-a and volume Vh with different cooling system

4 散热系统仿真与实验分析

4.1 散热系统仿真验证

为验证本文提出优化设计方法的实用性与本文提出综合热模型的准确性,使用热仿真软件Flothem建立散热系统仿真模型,利用仿真软件中的热源模拟每个SiC-MOSFET的损耗。设置环境温度为25℃,发热功率为864W,TIM热阻为0.6℃/W。图12为散热系统仿真温度场稳态分布。概要设计和优化设计的散热系统热仿真结果见表5。表中,本文模型中假设散热器表面温度均匀,实际中散热器表面温度并不均匀,为验证模型的准确性,记录仿真中散热器表面温度的最大值和最小值,并用式(2)计算散热器的最大热阻和最小热阻。定义散热器热阻平均误差为

width=160,height=35 (31)

式中,Rth,sim,max为仿真中得到的散热器最大热阻;Rth,sim,min为仿真中得到的散热器最小热阻。但是由于散热器表面温度分布不均,高温度附近的功率器件将会承受更高的热应力,因此将散热器最大热阻误差定义为

width=111,height=35 (32)

width=225.2,height=237.2

图12 散热系统仿真温度场稳态分布(Ta=25℃)

Fig.12 The steady-state simulated temperature field of cooling system (Ta=25℃)

表5 散热系统仿真结果

Tab.5 Simulation results of cooling system

设计方法温升/℃散热器热阻/(℃/W) MOSFET散热器 传统概要设计42.214.10~23.80.016 4~0.027 5 优化设计44.915.46~25.30.017 9~0.029 2

由图12仿真结果可以看出,两个散热系统的MOSFET壳最大温升均小于48℃,满足设计要求。本文提出综合热模型的理论计算热阻在仿真热阻值的区间内,平均热阻误差save=6.6%,即SVG满载运行时,散热器平均表面温升的绝对误差DTe在1.2℃以内。最大热阻误差smax=25%,即SVG满载运行时,散热器最大表面温升的绝对误差为DTe为5.6℃。仿真结果证明了本文提出的综合热模型的准 确性。

文献[7-8]的传统热模型得到的散热器热阻均高于仿真热阻的最大值,平均热阻误差save=67%,即SVG满载运行时,散热器平均表面温升的绝对误差DTe高达13.7℃。最大热阻误差smax=35%,即SVG满载运行时,散热器最大表面温升的绝对误差DTe为8.8℃。模型准确度较低。

4.2 散热系统实验测试

为验证本文提出优化设计方法的实用性与综合热模型的准确性,搭建380V/50kvar SVG测试平台,如图13所示。采用热电偶和温度巡检仪检测散热器表面位置的温度,温度稳定后,取散热器表面温度的最大值和最小值计算散热器的热阻,并记录温度最高MOSFET的壳温。

width=219,height=100.8

图13 散热系统测试平台

Fig.13 The test platform of cooling system

令SVG运行在开环有源逆变模式,向电网注入10~60kvar(即单相无功电流15~90A)的无功能量,通过功率分析仪(YOKOGAWA WT2030)测量SVG的输出无功功率和效率,即可计算出功率器件的总损耗。

不同功率下本文提出的综合热模型、文献[7-8]传统热模型、有限元仿真、温升实验测试的结果,如图14~图17所示。同时将本文提出的综合热模型和文献[7-8]传统热模型的理论温升DT,一并绘在图14、图16中。在不同发热功率下,两个散热系统的MOSFET壳最大温升均小于DTc,max=48℃,满足设计要求。由于散热器表面温度不均,高温度附近的功率器件将会承受更高的热应力,因此,设计时需为此预留裕度。

width=197.4,height=115.2

图14 概要设计散热器理论计算、热仿真和实验温升

Fig.14 Temperature rising of rough design model, simulation, and test

width=203.4,height=117.1

图15 概要设计散热器理论计算、热仿真和实验热阻

Fig.15 Thermal resistance of the model, simulation, and test

width=197.9,height=115.7

图16 优化设计散热器理论计算、热仿真和实验温升

Fig.16 Temperature rising of the model, simulation, and test

width=203.4,height=117.1

图17 优化设计散热器理论计算、热仿真和实验热阻

Fig.17 Thermal resistance of the model, simulation, and test

图15、图17中,本文提出的综合热模型理论计算热阻在实验测试热阻值的区间内,平均热阻误差save在7.0%以内,即SVG满载运行时,散热器平均表面温升的绝对误差DTe在1.3℃以内,较传统散热模型提高62%。最大热阻误差smax=26%,即SVG满载运行时,散热器最大表面温升的绝对误差DTe在5.8℃以内。

文献[7-8]的传统热模型得到的散热器热阻均高于实验测试热阻的最大值,平均热阻误差save=69%,即SVG满载运行时,散热器平均表面温升的绝对误差DTe高达14℃,最大热阻误差为smax=32%,即SVG满载运行时,散热器最大表面温升的绝对误差DTe为8.2℃。

实验结果验证了本文提出综合热模型的准确性,证明了该模型更适用于电力电子装置强迫风冷散热设计。

5 结论

电力电子装置热设计对高功率密度大功率电力电子装置的可靠性起着重要作用。为了提高热设计的准确性和设计效率,本文综合传热学与流体力学理论,提出了一种基于截面积二次方根为无量纲特征长度的综合热模型。同时,提出了一种针对典型强迫风冷散热系统的体积最优的优化方法,该优化方法亦可推广到质量最优、损耗最优等优化设计当中。以380V/50kvar高功率密度SVG为例,利用本文提出的方法设计的散热系统,较传统概要方法体积可缩小30%。对所设计的散热系统进行热仿真和温升测试,本文提出的综合热模型、热仿真、实验测试平均热阻误差save在7%以内,即SVG满载运行时,散热器平均表面温升的绝对误差DTe=1.3℃。相较于传统热模型[7-8]平均热阻误差save=69%,本文提出的综合热模型准确性有较大的提升。

参考文献

[1] Laloya E, Lucía Ó, Sarnago H, et al. Heat manage- ment in power converters: from state of the art to future ultrahigh efficiency systems[J]. IEEE Transa- ctions on Power Electronics, 2016, 31(11): 7896- 7908.

[2] 周林, 李寒江, 解宝, 等. SiC MOSFET的Saber建模及其在光伏并网逆变器中的应用和分析[J]. 电工技术学报, 2019, 34(20): 4251-4263.

Zhou Lin, Li Hanjiang, Xie Bao, et al. Saber modeling of SiC MOSFET and its application and analysis in photovoltaic grid-connected inverter[J]. Transactions of China Electrotechnical Society, 2019, 34(20): 4251-4263.

[3] 肖龙, 伍梁, 李新, 等. 高频LLC变换器平面磁集成矩阵变压器的优化设计[J]. 电工技术学报, 2020, 35(4): 758-766.

Xiao Long, Wu Liang, Li Xin, et al. Optimal design of planar magnetic integrated matrix transformer for high frequency LLC converter[J]. Transactions of China Electrotechnical Society, 2020, 35(4): 758-766.

[4] 胡建辉, 李锦庚, 邹继斌, 等. 变频器中的IGBT模块损耗计算及散热系统设计[J]. 电工技术学报, 2009, 24(3): 159-163.

Hu Jianhui, Li Jingeng, Zou Jibin, et al. Losses calculation of IGBT module and heat dissipation system design of inverters[J]. Transactions of China Electrotechnical Society, 2009, 24(3): 159-163.

[5] 丁杰, 张平. 电机控制器用IGBT风冷散热器的热仿真与实验[J]. 电源学报, 2015, 13(2): 38-44.

Ding Jie, Zhang Ping. Thermal analysis and experi- mental of IGBT air-cooled radiator for motor controller[J]. Journal of Power Supply, 2015, 13(2): 38-44.

[6] 何文志, 丘东元, 肖文勋, 等. 高频大功率开关电源结构的热设计[J]. 电工技术学报, 2013, 28(2): 185-191, 218.

He Wenzhi, Qiu Dongyuan, Xiao Wenxun, et al. Thermal design of high frequency high power switched- mode power supply[J]. Transactions of China Electro- technical Society, 2013, 28(2): 185-191, 218.

[7] 刘超, 贾晓宇, 胡长生, 等. 电动汽车SiC MOSFET风冷逆变器的散热器设计[J]. 电源学报, 2018, 16(3): 151-157.

Liu Chao, Jia Xiaoyu, Hu Changsheng, et al. Design of heat sink for electric vehicle SiC MOSFET air- cooled inverter[J]. Journal of Power Supply, 2018, 16(3): 151-157.

[8] Drofenik U, Stupar A, Kolar J W. Analysis of theoretical limits offorced-air cooling using advanced composite materials with high thermal conducti- vities[J]. IEEE Transactions on Components, Pack- aging and Manufacturing Technology, 2011, 4(11): 528-535.

[9] Kolar J W, Zach F C, Casanellas F. Losses in PWM inverters using IGBTs[J]. IEE Proceedings-Electric Power Applications, 1995, 142(4): 285-288.

[10] 白保东, 陈德志, 王鑫博. 逆变器IGBT损耗计算及冷却装置设计[J]. 电工技术学报, 2013, 28(8): 97- 106.

Bai Baodong, Chen Dezhi, Wang Xinbo. Loss calculation of inverter IGBT and design of cooling device[J]. Transactions of China Electrotechnical Society, 2013, 28(8): 97-106.

[11] 传热学. 电力电子器件热管理[M]. 北京: 机械工业出版社, 2013.

[12] Gammeter C, Krismer F, Kolar J W. Weight optimi- zation of a cooling system composed of fan and extruded-fin heat sink[J]. IEEE Transactions on Industry Applications, 2015, 51(1): 2193-2200.

[13] Muzychka Y S, Yovanovich M M. Laminar forced convection heat transfer in the combined entry region of non-circular ducts[J]. Journal of Heat Transfer, 2004, 126(1): 54-61.

[14] Muzychka Y S, Yovanovich M M. Pressure drop in laminar developing flow in noncircular ducts: a scaling and modeling approach[J]. Journal of Fluids Engineering, 2009, 131(11): 105-111.

[15] Simons R E. A simple thermal resistance model isoflux versus isothermal[J]. Electronics Cooling, 2006, 13(6): 10-14.

A Comprehensive Model of Forced Air Cooling and Optimal Design Method of High Power Density SiC-Static Var Generator

Lin Hongyi Wu Liang Guo Xiao Chen Guozhu

(College of Electrical Engineering Zhejiang University Hangzhou 310027 China)

Abstract Base on a novel characteristic length scale (the square root of cross-sectional), a comprehensive thermal model for the typical forced air cooling system of power electronic devices was proposed integrating heat conduction, convection heat transfer and fluid mechanics theories. In order to improve the accuracy and efficiency of the thermal design, an optimized design method for the optimal volume was proposed for the high-power and high-density power electronic device. Taking a 380V/50kvar high power density static var generator (SVG) as an example, compared with the experimental results, the average thermal resistance error of the comprehensive thermal model is within 7%, that is, the corresponding surface temperature rise is within 1.3℃ when SVG is running at full load. The accuracy of the comprehensive thermal model is improved by 62% compared with the traditional thermal models. Consenquently, the demention volume of the cooling part designed based on the proposed method is 4.03L, which is 30% smaller than the volume of 5.7L designed by the traditional methods.

keywords:Thermal design, thermal model, forced air cooling, static var generator (SVG)

DOI: 10.19595/j.cnki.1000-6753.tces.200821

中图分类号:TM762; O551

作者简介

林弘毅 男,1996年生,博士研究生,主要研究方向为电能质量、宽禁带器件应用。E-mail: lhy2007.11@qq.com

陈国柱 男,1967年生,教授,博士生导师,主要研究方向为大功率电力电子装置及其数字控制、有源电能质量控制技术及可再生能源发电并网技术。E-mail: gzchen@zju.edu.cn(通信作者)

国家自然科学基金资助项目(51777186)。

收稿日期 2020-07-09

改稿日期 2020-09-12

(编辑 陈 诚)