第35卷第7期 机 电 工 程 V01.35 No.7 2018年7月 Journal of Mechanical&Electrical Engineering Ju1.2018 DOI:10.3969/j.issn.1001—4551.2018.07.021 一种IGBT模块的新型建模与瞬态 仿真方法研究米 沈天浩 ,黄 瑞 ,杨 帆2,陈俊玄 ,冯 权 ,洪文华 ,俞小莉 (1.浙江大学能源工程学院,浙江杭州310000;2.上汽大众汽车有限公司,上海201805) 摘要:针对所采用的仿真方法无法全面评价IGBT模块可靠性的问题,提出了一种新型的IGBT模块仿真方法:对IGBT模块工作过 程中的瞬态损耗进行了分析和计算,并对IGBT模块进行三维建模与瞬态下流固耦合仿真。以FS800R07A2E3为例,采用一维仿真 分析计算其瞬态损耗,并将瞬态损耗等效为半正弦波,作为一维边界条件加载至三维仿真模型;同时基于Star—CCM+软件,采用瞬 态仿真与流固耦合仿结合的方法,仿真得到了模块内部温度场分布云图和芯片最高结温波动数据,并采用英飞凌实验数据验 证了仿真模型的准确性。研究结果表明:模块内部三相的不均匀性在2 ̄C左右,瞬态仿真下的最高结温比稳态仿真下高7℃左右; 该仿真方法更加接近IGBT工作的实际情况。 关键词:电机控制器;IGBT;结温;瞬态;流固耦合 中图分类号:TP391;TN32 文献标志码:A 文章编号:1001—4551(2018)07—0767—06 New modeling and transient simulation of IGBT SHEN Tian.hao ,HUANG Rui ,YANG Fan ,CHEN Jun.xuan , FENG Quan ,HONE Wen.hua ,YU Xiao.1i (1.College of Energy Engineering,Zhejiang University,Hangzhou 3 10000,China; 2.SAIC Volkswagen Co.Ltd.,Shanghai 201805,China) Abstract:Aiming at the lack of overall assessment on IGBT module simulation,a new type method was proposed:The transient loss of IGBT module in working process was analyzed and calculated and the fluid-solid coupling method was adopted.1-D simulation was used to calcu。 late the transient loss of Infineon FS800R07A2E3 Module.The transient loss was equivalent to half-sine wave and loaded into the 3-D model as boundary conditions.The fluctuation and distirbution of junction temperature were studied by the coupled method of trnasient simulation and 3-D simulation based 0n Star—CCM+software.Meanwhile.Infineon experimental data was adopted to verify the accuracy of the mode1. The resuhs indicate that a step change about 2。C in U/ ,W phases and the m ̄imum lunction temperature in transient simulation is about 7 oC higher than the steady state.The actual working situation of IGBT is better simulated by this new method. Key words:MCU;IGBT;junction temperature;transient;fluid—solid coupling 合适的IGBT模块对电动汽车电控模块的设计至关 0 引 言 重要。 在电动汽车开发的前期,需要对电池、电机、电控 IGBT在工作过程中会产生损耗导致其温度升高。 进行选型与匹配。IGBT模块是电控的关键部件,选择 IGBT模块由封装好的IGBT芯片和FWD(续流二极 收稿日期:2017一ll一24 基金项目:国家自然科学基金资助项目(51476143) 作者简介:沈天浩(1993一),男,江苏苏州人,硕士研究生,主要从事电动汽车及其热管理方面的研究。E—mail:shentianhao1993@126.corn 通信联系人:俞小莉,女,工学博士,教授。E—mail:yuxl@zju.edu.cn 机 电 管)芯片组成。结温是电子设备中半导体芯片PN结 的工作温度,当结温过高时会导致硅材料芯片性能降 低,乃至永久损坏。据统计,由温度过高引起的IGBT 失效占全部的55%以上H J。电动汽车设计前期只能 采用仿真的方法对IGBT模块进行结温预测,分析不 同工况下的发热与模块结温是否符合要求。 通过仿真对IGBT模块的结温进行预测是国内外 研究的重点。姚芳 、宋飞 和陈彦 搭建了IGBT 模块的热电联合仿真模型,得到IGBT与FWD的损耗 波形和结温波形,但其仅针对一维模型,未进行三维流 固耦合仿真,无法得到温度场分布;陈清 和王淑 旺 对某IGBT模块进行散热分析,得到其温度场和 流场分布,但其仅计算稳态工况,未计算瞬态结温。 IGBT模块工作时,其芯片发热会随时间作周期性 波动,芯片的结温也随之周期性变化。目前的IGBT 模块仿真方法包括一维和三维仿真,但无法同时得到 IGBT模块的结温分布云图和结温波动数据。 本文以型号为FS800RO7A2E3的Ir{ i;1 6 5 4 3 2 1 O 9 8 7 6 5 GBT模块为例,4 3 2 1 提出一种新型IGBT模块仿真方法。O 0 O O O 0 0 0 O O 0 O O O O 0 O O O O O 0 O O 0 O 0 0 O O O O O 1 损耗计算 IGBT模块在运行时会产生损耗,损耗主要包括 IGBT芯片和FWD芯片的损耗,是IGBT模块中的 热源。 1.1损耗计算理论 三相桥式整流电路如图1所示。 图1三相桥式整流电路 每相包含两个IGBT芯片和两个反向并联的续流 二极管(FWD),二者工作时经历周期性的开关状态并 产生损耗。 IGBT芯片的主要损耗是通态损耗 71和开关损 耗 8 J,IGBT的通态损耗是由于其正向导通压降产生 的,其开关损耗是因为在开通和关断的瞬间,电流和电 压有重叠期;FWD芯片的主要损耗是通态损耗和关断 损耗,FWD的通态损耗也是由于其正向导通电压而产 生的,其开关损耗则是因为二极管的反向恢复特性。 目前,电动汽车大多采用SVPWM(空间矢量)调 工 程 第35卷 制方式,其损耗可以通过计算得到 。 IGBT输出特性如图2所示。 I _,_一1‘ l | , I | 】 o℃l | | , |。 | ,// /, | ’f , ∥ O 0 0 2 0 4 0 6 0 8 1 0 1 2 1 4 1 6 1 8 2 0 2 2 2 4 2 6 V 图2 IGBT输出特性 为了简化损耗的计算,IGBT的导通压降 与 FWD的正向导通电压 与, 的关系可以近似用直线 表示: VcE(t)=VcE0+rcE:I:,c(t) (1) (t)= + ,c(t) (2) 式中: 。一IGBT门槛电压; 一FwD门槛电压; r 一IGBT通态等效电阻;r 一FwD通态等效电阻。 当采用传统7段SVPWM调制时,IGTB芯片和 FWD芯片的损耗通过计算获得。 (1)IGBT通态损耗。 当一詈≤ ≤詈时: 。 = +÷ +扣 ,Nc + f\ 48,rr ,…’ (3) 当 6≤ ≤詈或一号≤ ≤一詈时: 邶 =Vc瑚,N+ + c。s + f\ 481T 1,… 腑旺 (4) 式中: 一正弦波相电流幅值;M一调制比; 一电流 滞后电压角度。 (2)IGBT开关损耗。 Pswf sw( IN IGBT=)× × … + …(5) 式中 一开关频率;E……--IGBT在标称电压 和 第7期 沈天浩,等:一种IGBT模块的新型建模与瞬态仿真方法研究 ·769· 电况In 卜测得的升迪损稳; 0仃… --IGBT在标称电 1.2 Fs80o1 D7A2E3模块损耗计算 1.2.1 稳态损耗 压 。 和电流,n。 下测得的关断损耗; 一直流 电压。 FS800R07A2E3模块中IGBT芯片和FWD芯片的 特征参数分别如表1、表2所示(特征参数测试条件为 300 V,550 A,125 oC) 表1 FSS00R07A2E3模块IGBT芯片特征参数 参数 ∞(3)FWD通态损耗。 当一詈≤ ≤詈时: Pc0 』w。= ,N+ 1 r 一÷ 。s 一 1 , ” 数值 O.82 9.8×10一 /V f\ 48 二蝤a'r (6)rOE/12 当詈≤ ≤ ,f或一号≤ ≤一詈时: P。 。= ,N+ 1 rF 一÷ ,Nc。s 一 f\ 些 48叮『 1, 胁 一’ (7) (4)FWD开关损耗。 P.w VWD=f sw …om× × (8) 式中:E 。一 一FwD在标称电压 。 和电流,F… 下测 得的反向恢复损耗。 由式(3~8)可得IGBT芯片和FWD芯片在一个 相电流周期内的平均损耗,分别表示为: PIGBT=P。。 IGBT+P IGBT (9) PFwD=P 。 _FwD+P洲_FwD (1O) 上式计算的是IGBT模块工作时的平均损耗。三相 桥式逆变电路的基本工作方式是180。导电方式,即每 个桥臂的导电角度为180。,同一相(半桥)上下两个臂 交替导电,各相开始导电的角度依次相差120。。因此, 一个IGBT桥臂在一个周期内只有一半时间是导通的。 IGBT芯片与FWD芯片的损耗是不规则的曲线, 可以将两者的损耗函数简化为理想的半正弦波曲线。 两者瞬态损耗的函数为: f"trPIGBTsin(tot)0<£<÷ GBT i 。 … (11) i【f  ̄rPv0 wDs0 <in(art一 )一订)寺< < 寺 <T (12) 式中:P 。 (t)--IGBT芯片瞬态损耗;P w。(t)一FWD 芯片瞬态损耗;卜相电流周期。 E。…。 /mJ 10.5 0ff一/mJ 24.5 表2 FS800R07A2E3模块FWD芯片特征参数 参数 数值 I/V 1.04 rr/Q 5.5×10一 一 /mJ 12.5 IGBT模块的部分工作参数如表3所示。 表3 FSS00ROTA2E3模块工作参数 参数 数值 1 COS 0.8 Lw/kHz lO fo/Hz 50 注: 一调制比; 一电流滞后角 一开关频率 一相电 流频率 不同直流电压 和相电流幅值,N工况下(工况1 至工况3),IGBT芯片与FWD芯片的平均损耗如表4 所示。 表4 不同工况下FS800R07A2E3模块的损耗 由此可见:IGBT和FWD的平均损耗均随 和, 的增大而增大。 1.2.2 瞬态损耗 本研究取工况1至3中发热最大的工况3计算瞬 态损耗。根据一个周期内平均损耗功率相同的原理进 行等效,可知等效半正弦波的幅值为平均损耗的仃倍。 其他工况相同时,瞬态损耗的变化周期与相电流频率 有关。频率为50 Hz时,瞬态损耗拟合式为: P 。 ct =1 697·4si n 。。 ¨。.00 ’ 102 (13) 尸 。c :{47。.6 in 00。竹 一订 。.0 ‘ l02 (14) IGBT模块瞬态损耗IH1线如 3所示 3 IGB…1模块t瞬念捌牦曲线 2 维建模 在进行_t维建模前,埘 Brr模块的 维模型做 f… 下假设:(1)心it-片足一个均匀发热体,芯片rf1所 有层都是均匀的,无缺陷的;(2)芯片较薄,除r芯片 截面法向外,其他表面均视为绝热;(3)忽略辐射和 空气对流换热;(4)冷却液为不il『J匠缩流体且为湍流 流动 2.1 建模与网格划分 本研究以英 凌FS8OOR07A2E3模块为例,建立 其 维模型。通过测量及卡¨父资料得到模块的尺寸信 息;采 SolidWorks建 模块的一 维模型并进行一定 的简化,去除密封 槽,去除模块内部的连接线等刈‘散 热影响较小的结构;然后将I(;Bq 几何模型导入到 HyperMesh I}1进行网格划分,最小网格尺寸为0.1 111111 维模 和网格如 4所永. ..一 【:1)JL1- J懈’ Ib) 饼模型 4 蜒 凌FS8(X)R07A2 E3模块模 2.2 流固耦合计算模型 笔者采}{j Star—CCM+软件埘IGBq、模块进行瞬 态 维流同耦合仿真,汁箅 不 _l厂况、一定冷却条件 下的结温分布和波动. .计算模型巾的冷却液为 缩流体,满足质量 守 陋、动量守恒、能量守 程,即: 塑+ ( )=03 t 、 (15) 丛 =,) — ,J+ △ (16) 第35巷 : +A△7’+Pg (17) 『fI:p~流体的密度;£,一速度矢量;,,一流体f 力; 流体的动力 : 度 一流体的比热容;A一导热 系数;F~作川 流体I 的质量力;(『一流体所吸收的 热 ;卜流体或『IH体温度; 一能量耗散闲数 冷却液 f 机 制 中的流动为湍流,采川基1 雷诺时均法的数仳汁锋力‘法,引入湍流模型,计算巾采 J_Ij 一 ;流模 ,即: 苦( )+未( ,J)=毒[( + o' ̄/13丛:ri]+ + —p占 (1 8) p =未[( +(r ]13熊x]+C }(G +G! G )一c (19) 式小:IIi一£/ 力州㈨,J速度分量; 一湍动能; 一耗 散牢; ,一湍流 ;G 一m浮生力产生的湍流动能; G 一平均速度悌度f』J起的湍动能 的产生项;or 、 一 湍流普朗特数;C ,C、 一经验常数 Stm’一CCM+软件rf1,流体域采川 一 湍流模 ,边 层没 为2 J 2.3物性参数设置 小研究,将HyperMesh巾的面网 导人到Star— CCM+ ,没 卡f1 材料的物性参数。物性参数设 的对象包括1人】邴芯片、外壳材料、导热基板、冷却液等, 物性参数包括材料的密度、导热系数、比热容等 .『Ii于 1GBT芯片的多层结构,网格密度无法细小至此,笔者 采川复合导热系数的¨‘ 方法对其进行折算“ 芯片材料为 ,屣板材料为铜,外壳材料为r程塑 料,冷却液为50%水 Ij 50%乙二醇的混合物 2.4边界条件设置 FS8OOR07A2E3卡I!块为 相全桥整流,每个仑桥 南两个半轿绀成,·1\、 桥F}1 4个IGBT芯片和I 4个 FWD 片组成 桥内部的实际结构和简化 维模 型如 5所示. .( l】 II,Ji结 {b)简化模型 5 蜒 凌FS800R07A2E3模块 卜}乔 本研究将J 述IGBT和FWD芯片的发热模型作 第7期 沈天浩,等:一种IGBrr模块的新型建模与瞬态仿真方法研究 为边 条件,JJ【1载至 维模型,稳态加载为常数,瞬态 加载为时间的函数 电机控制器的入口温度一般要控制在65 以下, 本文计算最恶劣1_况,即冷却液入u温度为65 cc,流 量为l0 L/l11iII. 3 仿真结果 3.1稳态仿真结果 本研究将f:述稳态T况的发热功率设定为边界条 件,对其进行仿真 在迭代500步以后仿真收敛,得到 IGBT模块的温度分布云冈 :流同耦合计算得到的温 度分布如 6所示 , le|lIn。r“I1 ℃ 68‘)74 79 5 77 q0 l 8l l【】(J 78 ■■翻i l 。罔6稳态仿真温度分布云 【殳I 6 rfl,左侧为冷却液入口,右侧为 [] 定义冷 却液人Ej到 {【_1分别为U/WW相,U/I,VW三相的最 高温度存在温差,稳态下最高温度出现在 相的lG— BT芯片,最高结温约为122 oC 各相芯片的最高温度如表5所示 表5各相芯片最高温度 位 最高结温/ ̄C £卡f1 1 17.37 I 扎l l l9.39 卡日 l21.99 相温度的不均匀性在2 左存,分析其原凶可 能为£,相靠近冷却液入口,冷却液温度较低,I大1 带 走热量较多 3.2实验结果及分析 通过奁阅数据表r叮以获得不同冷却液流速下的热 阻,结合冷却液温度和散热功率,可以汁算模块的最高 结温一仿真数据与实验数据对比如表6所示 表6三维仿真模型验证 序号英飞凌试验最高结温/ : 维仿真最高结温/℃误差/% 表6中,本文三维仿真模型得到的稳态仿真结果 与厂家实验数据之间的误差在2%以内,验证了三维 模型的正确性。 3.3瞬态仿真结果 在稳态计算中, 相远离冷却液入口,温度最高, 本文仅计算 相最高结温波动。将上述瞬态损耗公 式加载至 相,£/相和 相的损耗仍采用稳态损耗加 载。瞬态时间步设置为0.001 S,每个时问步包含50 步迭代步数。由于从冷起动开始计算需要消耗大量计 算资源和时问,冈此将初始条件设置为稳态计算的 结果, r}1于瞬态热源的周期性波动, 相IGBT芯片和 FWD 片的温度也周期性波动。瞬态结温波动数据 如 7所示。 播 饱 啦 ,,S 7 相芯片结温波动 7中,相电流频率为50 Hz时,热源的周期为 0.02 S,结温波动的周期也为0.02 S。在前三分之一 周期内IGBT芯片温度上升,后 分之二周期内温度 下降,最高结温大约f【l现在周期的j分之一处。IGBT 芯片的最高结温达到129 oC左右。 由于同一相中IGBT芯片温度高于FWD芯片,本 研究考察不同频率下IGBT 片的结温波动 不同相 电流频率下IGBT芯片结温变化如图8所示。 I/S 8 不同相电流频率下结温波动 ·772· 机 电 工 程 electroniCS Journal,2000,31(9):781-785. 第35卷 图8中,相电流频率越高,最高结温越低。可以解 释为:在芯片发热的一个周期内,前部分周期芯片发 热,冷却液带走的热量比发热少,导致结温升高;后部 分周期内芯片不发热,冷却液持续带走热量,导致结温 降低。 [2] 姚芳,王少杰,李志刚.逆变器中IGBT功率模块的电热 联合仿真模型[J].半导体技术,2016,41(6):440445, 472. [3]宋[4]陈飞,梁哲兴,张彦,许伟.IGBT器件稳态及瞬态热模型仿 真分析[J].船电技术,2013,33(3):10—13. 鹏,李世乎.基于Matlab的IGBT能耗模型 及其在热仿真中的应用[J].大功率变流技术,2014(5): 33.35. 4 结束语 本文以英飞凌FS8OOR07A2E3模块为例,提出了 一种新型的IGBT模块的仿真方法,得到以下结论: [5] 陈清.电动汽车用IGBT模块液冷散热及封装可靠性 (1)IGBT模块内部三相芯片的最高结温存在2 ̄C 左右温差,一维瞬态计算无法模拟三相间温差与温度 分布云图; (2)IGBT模块进行选型时,仅采用稳态仿真是不 够的,其结温波动需要通过瞬态仿真得到; (3)IGBT的实际结温存在周期性波动,并且与相 电流的频率相关。在电机控制器设计阶段,可以考虑 研究[D].重庆:重庆大学汽车工程学院,2016. [6] 王淑旺,赵卫健,唐志国,等.纯电动汽车水冷电机控制器 的热仿真和热分析[J].微特电机,2013,41(9):11 13. [7] 白保东,陈德志,王鑫博.逆变器IGBT损耗计算及冷却装 置设计[J].电工技术学报,2013,28(8):97—106. [8]李志刚,梅[9]焦明亮,李霜,王少杰,等.IGBT模块开关损耗计算方 云,吴春冬,等.一种电动汽车逆变器IGBT 法综述[J].电子技术应用,2016,42(1):10—14,18. 功率损耗和结温的近似计算方法[J].大功率变流技术, 2015(5):58-63. 采用提高频率的方法以降低IGBT模块的最高结温; (4)采用瞬态仿真和流固耦合结合的方法,既能 仿真得到结温分布云图,也可以获得结温波动数据。 参考文献(References): [1] JANICKI M,NAPIERALSKI A.Modelling electronic circuit radiation cooling using analytical thermal model『J].Micro. [1O]朱元,李根生,吴志红,等.两种控制策略下电动汽车 逆变器功率损耗的近似计算[J].汽车安全与节能学报, 2013,4(1):87—92. [11] 张波.多层薄板的导热系数[J].大连轻]二业学院学 报,1994,13(1):121 124. [编辑:张豪] 本文引用格式: 沈天浩,黄瑞,杨帆,等.一种IGBT模块的新型建模与瞬态仿真方法研究[J].机电工程,2018,35(7):767—772. SHEN Tian。hao,HUANG Rui,YANG Fan,et a1.New modeling and transient simulation of IGBT[J].Joumal of Mechanical&Electrical Engineering.2018, 35(7):767—772t 《机电工程》杂志:http://www.mecm cI1m.(:n