摘要:在国内外电力电子领域对功率二极管模型研究成果的基础上,利用计算机仿真技术,从数学物理模型和电路模型两方面,对快速功率二极管的反向恢复特性进行了较深入的研究,获得了　　　　可正确描述正反向恢复过程的功率二极管仿真模型.该模型克服了标准二极管模型完全忽视正向恢复效应,对二极管反向恢复现象的模拟也会产生错误振荡的缺陷,因此具有较好的现实意义.

1 引言

　　在弧焊逆变器中,大功率快恢复二极管的瞬态过渡过程和特性对于高频整流和主开关器件的正常工作具有至关重要的影响.在igbt开通瞬间,功率二极管处于续流状态,会引起变压器次级短路,对 igbt形成电流冲击,不利于igbt的可靠运行.

　　功率二极管存在显著的电导调制和电荷存储效应,其开、关状态的改变需要一定的时间[1].在正向恢复过程中,开通初期出现高出正常通态压降2~10倍的电压尖峰;而在反向恢复过程中,关断初期不能承受反向电压而产生较高的反向电流.为准确表征弧焊逆变器的动态过程,功率二极管的仿真模型必须能正确反映正向和反向的恢复特性.

　　spice标准二极管模型完全忽视了正向恢复效应,对二极管反向恢复现象的模拟也会产生错误的振荡[1].国际上许多文献对功率二极管的实用仿真模型进行了大量研究”[1~4],其中基于集中电荷的概念,根据功率二极管内部载流子的简化输运方程,成功地给出了可以同时正确描述二极管正、反向恢复瞬态过程的仿真模型.国内针对spice标准二极管模型的不足而展开深入讨论的文献不多[5].以下将以 c l ma和p o lauritizen的系列研究为基础,对功率二极管的电路仿真模型进行深入探讨.

2 功率二极管的数学物理方程

　功率二极管的数学物理方程由反向恢复、正向恢复、发射区复合方程以及包括接触电阻和结电容的总变量方程等4部分组成.

2.1 反向恢复的方程

　　扩散电流可由下式决定:

　　(1)

　　式中 qo——pi界面附近的电荷变量

　　q2——i/2区的电荷变量

　　t12——扩散穿越时间常数

　　对于q2,电荷控制的连续性方程为:

　　(2)



　　式中 ——电荷变化率

　　qo——pn结的注入电荷水平变量

　　——电荷在寿命,期间的复合率

　　——p+i结到q2的扩散电荷

　　由pn结定律可得:

　　　(3)

　　式中iso——二极管的饱和电流常数

　　t——电荷寿命

　　vt——热电势常量

　　ve——结电压

　　由式(1),式(2)和式(3)组成的方程组就是完整描述功率二极管反向恢复过程的数学物理模型.

2.2正向恢复的方程

　　在二极管正向过渡期间,由于i区的初始导电性低,在二极管两端建立起较高的电压.随着注入载流子浓度的增加,i区上的电压很快降低到二极管正常稳态正向电压.按照载流子漂移导电方程,

　　可得:

　　(4)

　　式中vm——i/2区的电压

　　rmo——i区的初始电阻常量

　　ivd——二极管总电流

　　式(4)是功率二极管的正向恢复方程,说明i区上的电压直接取决于该区的总电流ivd和电荷q2.

2.3发射区的载流子复合

　　当功率二极管的导通电流很高时,由于载流子注入到重搀杂的p+和n+发射区.此时,发射区的载流子复合电流ie不可忽略.根据发射区的注入电子浓度变化方程,可得:

　　 (5)

　　式中ise——电流常数

　　通常,ise比iso小几个数量级,也就是说,ie仅在非常大的电流下才成为总电流的主要部分.

2.4 包括接触电阻和结电容的总变量方程

　　功率二极管完整的数学物理模型应该将接触电阻rs和结电容cj也考虑进去.由此可得,二极管端口电压v和总电流i的输出变量方程分别为:

　　(6)

　　综上所述,由式(2)~(6)构成的方程组,就是完整描述功率二极管动态过程的数学物理模型.

3二极管正反向恢复spice仿真模型

　　运用非线性受控源b元素,spice3.0以上版本具有直接将数学方程转化为电路仿真模型的功能.

　　例如,假设设计的仿真模型中节点1的电压表示数学方程(3)中的变量q0,节点2的电压表示数学方程中变量ve,则式(3)对应的仿真模型为:b10v={iso}*|τ*(exp(v(2)vt)-1).

　　由此可得图1所示的二极管正、反向恢复特性的spice仿真子电路模型.

　　其模型清单为:

　　subckt drec 1 9

　　* connection ac

　　* following components include for space charge capacitance

　　dmodel 1 2dcap

　　model dcap d(is=1e-21rs=0 tt=o cjo={cap})

　　bdl21=(v(5,6)/{tm})+({ise} x(exp((v(1)-v(2))/{vta})-1))

　　* following components model reverse recovery

　　be50v={isl}*{tau}*(exp((v(1)-v(2))/(2*{vta}))-1)

　　|;set v(5)=qe

　　re 50 1e6

　　bm 6 0 v