关键词:
IGBT
动态模型
双极输运方程
热模型
网络拓扑法
降阶方法
摘要:
电力电子装置与系统是实现电能有效变换和有效利用的关键设备;而大容量电力电子装置与系统更是我国重大装备的关键基础性技术之一,是面向经济社会发展的重大需求。绝缘栅双极型晶体管(IGBT)则是现代电力电子装置与系统的核心电力半导体器件。(动态)性能的精确分析和计算是IGBT设计的前提和理论与技术基础。因此,为提高我国IGBT等场控器件的设计与制造水平,以及相关电力电子产品在国际市场上的竞争力,亟需开展IGBT等场控器件多物理场计算模型和计算方法研究。IGBT器件及其电力电子系统的分析和设计所涉及的多物理场十分复杂。首先,IGBT导通、关断过程是极小时间尺度的电流场,其控制方程为载流子输运方程和电荷守恒定律。其次,极小时间尺度下电磁过程的位移电流与集肤效应的分析与计算等必须借助于三维涡流场的计算理论和计算方法。第三,电流场、涡流场和温度场相互耦合与相互影响。第四,从场的物理本质来看,这些多物理场(电流场、涡流场和温度场)都是复杂的三维场。第五,不仅不同物理场的时间尺度不同,即使是同一物理场,如载流子输运方程对应的电流场,在开通、关断瞬间以及持续开通时段内的时间尺度也不一样。而上述这些关键理论和问题是目前计算电磁学的研究热点,也是迄今尚未完全解决的难点。需要说明的是,目前电力电子装置和系统主要聚焦于长时间尺度下的热-电磁、热-力和力-电磁场之间的相互作用关系分析,极少涉及中、小时间尺度下的多物理量的互动关系研究,导致多时间尺度下的多物理场耦合机制不清楚。有鉴于此,为解决我国高性能IGBT等场控器件开发设计中的关键理论和技术,结合国家自然科学基金重大课题“多时间尺度下大容量电力电子混杂系统运行匹配规律研究”,本文重点解决IGBT模块及其系统多时间尺度多物理场的计算模型和计算方法。主要研究内容和成果如下:1、研究了 IGBT模块动态多时间尺度电流场数值计算的计算模型和计算方法。基于半导体物理的载流子输运方程(ADE),在国内外率先考虑了载流子二维分布以及非线性等复杂因素的影响;计算模型能够处理载流子全注入和真动态等复杂运行工况,实现了“真动态、全工况,以及二维载流子非线性运动”等复杂过程中IGBT模块内部载流子分布的分析与计算。为考虑小时间尺度下外部电路杂散电感和杂散电容的影响,建立了ADE-外电路耦合模型,并提出了相应的求解策略。此外,为满足不同工程问题的分析和计算需要,在上述模型和算法的基础上进一步提出了不同复杂度的简化模型和计算方法。2、研究了基于网络拓扑法的IGBT模块三维电场和温度场数值分析的计算模型和计算方法。绝缘栅双极晶体管(IGBT)具有高功率密度和低损耗等优点,在高频、高功率等场合获得了广泛应用。结温是预测IGBT可靠性的关键参数,也是导致IGBT失效的关键应力源。因此,建立IGBT的动态温度场模型并准确预测器件结温具有现实的应用价值。为此,首先提出了适合IGBT温度场分析计算的网络拓扑法,包括不同媒质及边界条件的处理方法。随后,构建了基于网络拓扑法的中时间尺度IGBT温度场数值分析的热模型,同时考虑了媒质特性参数的非线性和三维热分布的影响。最后,应用前述的理论和方法实现了典型IGBT模块的动态温度场-电流场的耦合分析和计算,并与商用软件ANSYS的计算结果进行了对比分析,验证了本文热模型和电热模型的准确性与快速性。3、研究了宽频ENOR降阶理论和算法。无论是有限元法,还是有限差分法,亦或是本文的网络拓扑法,其计算多物理场的核心是将连续域的偏微分方程的定解问题转化为大型代数方程组,然后求解。对于复杂的三维多物理场问题,其离散代数方程组的规模巨大,计算机求解需要庞大的计算资源和计算时间。因此,现有数值方法从某种程度上来说不便于工程问题的分析与计算。故,近年来大型代数方程组(模型)降阶理论和算法成为计算电磁学的研究热点。为此,基于IGBT瞬态多物理场动态响应的“宽频”特征,本文提出了模型降阶的宽频ENOR降阶理论和算法,并应用于前述的IGBT热模型。本文的降阶理论和算法(宽频ENOR降阶算法)基于经典的ENOR降阶算法,通过引入工作频域拓宽和自适应调节算法,实现了以简单(降阶)模型复现宽频动态响应特性的计算目标。4、研究了多时间尺度IGBT三维涡流场计算的有限元-分布参数耦合模型和计算方法。随着IGBT开关频率的增加,基于IGBT的电力电子能量变换装置的电磁过程覆盖的频域越来越宽,其动态过程的时间尺度从毫秒级上升到微秒级直至纳秒级。相应地,IGBT模块(包括本体)极小时间尺度下电磁过程的位移电流与集肤效应,以及引线的邻近效应等影响将不可忽略,即现有的集中参数电路模型将不完全适用于多时间尺度,尤其是极小时间尺度电磁过程的分析与计算。为此,本文研究了IGBT多时间尺度,尤其是极小时间尺度动态电磁行为的数值计算模型和计算方法。为考虑