CN108649807B - 一种牵引变流器器件结温在线计算方法及系统 - Google Patents
一种牵引变流器器件结温在线计算方法及系统 Download PDFInfo
- Publication number
- CN108649807B CN108649807B CN201810490961.3A CN201810490961A CN108649807B CN 108649807 B CN108649807 B CN 108649807B CN 201810490961 A CN201810490961 A CN 201810490961A CN 108649807 B CN108649807 B CN 108649807B
- Authority
- CN
- China
- Prior art keywords
- indicate
- power consumption
- bridge arm
- diode
- junction temperature
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M5/00—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases
- H02M5/40—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc
- H02M5/42—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters
- H02M5/44—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters using discharge tubes or semiconductor devices to convert the intermediate dc into ac
- H02M5/453—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters using discharge tubes or semiconductor devices to convert the intermediate dc into ac using devices of a triode or transistor type requiring continuous application of a control signal
- H02M5/458—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters using discharge tubes or semiconductor devices to convert the intermediate dc into ac using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only
- H02M5/4585—Conversion of ac power input into ac power output, e.g. for change of voltage, for change of frequency, for change of number of phases with intermediate conversion into dc by static converters using discharge tubes or semiconductor devices to convert the intermediate dc into ac using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only having a rectifier with controlled elements
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02M—APPARATUS FOR CONVERSION BETWEEN AC AND AC, BETWEEN AC AND DC, OR BETWEEN DC AND DC, AND FOR USE WITH MAINS OR SIMILAR POWER SUPPLY SYSTEMS; CONVERSION OF DC OR AC INPUT POWER INTO SURGE OUTPUT POWER; CONTROL OR REGULATION THEREOF
- H02M7/00—Conversion of ac power input into dc power output; Conversion of dc power input into ac power output
- H02M7/42—Conversion of dc power input into ac power output without possibility of reversal
- H02M7/44—Conversion of dc power input into ac power output without possibility of reversal by static converters
- H02M7/48—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode
- H02M7/53—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal
- H02M7/537—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters
- H02M7/5387—Conversion of dc power input into ac power output without possibility of reversal by static converters using discharge tubes with control electrode or semiconductor devices with control electrode using devices of a triode or transistor type requiring continuous application of a control signal using semiconductor devices only, e.g. single switched pulse inverters in a bridge configuration
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Inverter Devices (AREA)
Abstract
本发明涉及电力电子技术领域,公开了一种牵引变流器器件结温在线计算方法及系统,以实现三电平牵引变流器变频调速系统中关键器件结温波动的在线计算,进一步检测器件的健康,保证系统的稳定性和鲁棒性;针对牵引变流器目标桥臂上的目标器件,建立目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型;建立目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型;在非第一周期的其他周期内,计算该周期所对应的目标器件到所在外壳之间的第一等效热阻抗、以及目标器件由外壳到外部环境网络之间的第二等效热阻抗;然后计算目标器件在该周期的结温值。
Description
技术领域
本发明涉及电力电子技术领域,尤其涉及一种牵引变流器器件结温在线计算方法及系统。
背景技术
牵引传动系统被誉为轨道交通装备的“心脏”,不仅是整个轨道交通装备的核心动力单元,更是其运行安全的关键系统之一,对牵引传动系统中关键器件的健康监测至关重要。牵引传动系统由牵引变压器、牵引变流器等主电路设备,以及牵引电机等装置构成,其中,牵引变流器包括脉冲整流器、中间直流环节以及牵引逆变器。因此,牵引变流器是这枚“心脏”的“起搏器”,由于长期在高温、高压等恶劣的工作环境下,牵引变流器是牵引传动系统的高发故障源。牵引变流器中关键器件,主要是开关管和二极管的性能退化带来的安全隐患若不能及时监测到并得到及时正确的处理,都有可能引发器件的损坏失效,进而导致系统发生连锁事故。因此,保证牵引变流器及其关键器件安全可靠地工作不仅是确保牵引传动系统安全运行的关键,也是确保轨道交通装备安全运行的关键。
由开关管和二极管的失效机理分析可知,55%以上的失效原因与器件的结温波动有关,即,器件结温的频繁波动将加速器件键合线脱落和焊料层出现空穴等老化问题,进而引起器件的性能退化,最后导致器件的损坏失效。因此,对关键器件进行结温波动的在线监测将为器件剩余寿命、老化程度的估计和失效机理的分析提供了重要的数据支撑。目前,直接在线测量开关管和二极管的结温的方案大多需要增加额外的测温传感器或模块,虽然精度较高,但成本也较高,且不符合当下及未来工业界对关键器件轻量化和集成化设计的需求。因此,建立系统中关键器件的功耗和结温等模型,通过系统中已有的传感器模块,间接评估器件结温波动成为一种有效的途径。但目前的研究成果主要为两电平的变流器拓扑结构,且集中于恒频工况和采用正弦脉宽调制(SinusoidalPulse Width Modulation,SPWM)策略的应用环境,如风力发电,电网能量传输等,难以适用于三电平,变频调速工况和采用空间矢量脉宽调制(Space Vector Pulse Width Modulation,SVPWM)策略的应用环境,如CRH2型高速列车牵引传动系统即为变频调速系统,采用的三电平牵引变流器拓扑,其中,脉冲整流器采用SPWM调制策略,牵引逆变器采用SVPWM调制策略。
因此,现需提供一种能实现三电平牵引变流器变频调速系统中关键器件结温波动的在线计算与监测的牵引变流器器件结温在线计算方法及系统。
发明内容
本发明目的在于提供一种牵引变流器器件结温在线计算方法及系统,以实现三电平牵引变流器变频调速系统中关键器件结温波动的在线计算,进一步检测器件的健康,保证系统的稳定性和鲁棒性。
为实现上述目的,本发明提供了一种牵引变流器器件结温在线计算方法,包括以下步骤:
针对牵引变流器目标桥臂上的目标器件,建立所述目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型;所述目标器件为二极管或开关管,所述功率损耗包括通断功耗和通态功耗;
考虑电热耦合,建立所述目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型;
在非第一周期的其他周期内,根据所述第一、第二及第三关系模型计算所述目标器件在该周期的平均功率,并计算该周期所对应的所述目标器件到所在外壳之间的第一等效热阻抗、以及所述目标器件由外壳到外部环境网络之间的第二等效热阻抗;然后根据所述平均功率、第一等效阻抗及第二等效阻抗计算所述目标器件在该周期的结温值。
与上述方法相对应地,本发明还提供一种牵引变流器器件结温在线计算系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。
本发明具有以下有益效果:
1、本发明提供的一种牵引变流器器件结温在线计算方法及系统,实现了计及电热耦合下牵引变流器关键器件结温波动的在线监测,能在线计算牵引变流器中目标器件的结温变化。
2、本发明提供的一种牵引变流器器件结温在线计算方法及系统,仅利用系统中已有的电流传感器数据和各开关管的驱动控制信号,不受系统调制策略和恒频工作条件的约束,能实现变频和不同调制策略条件下牵引变流器关键器件结温的在线计算,提高了器件结温在线计算的灵活性,为牵引变流器系统的健康监测与管理提供了有力的科学支撑。
下面将参照附图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例的牵引变流器器件结温在线计算方法的流程图;
图2是本发明优选实施例的三电平牵引变流器系统的主电路拓扑结构图;
图3是本发明优选实施例的牵引逆变器U相上开关管TU1和TU2每个计算周期内功耗的变化示意图;
图4是本发明优选实施例的牵引逆变器U相上开关管TU3和TU4每个计算周期内功耗的变化示意图;
图5是本发明优选实施例的牵引逆变器U相上二极管DU1和DU3每个计算周期内功耗的变化示意图;
图6是本发明优选实施例的牵引逆变器U相上二极管DU12和DU34每个计算周期内功耗的变化示意图;
图7是本发明优选实施例的牵引逆变器U相上开关管TU1和TU2每个计算周期平均功率变化示意图;
图8是本发明优选实施例的牵引逆变器U相上开关管TU3和TU4每个计算周期平均功率变化示意图;
图9是本发明优选实施例的牵引逆变器U相上开关管DU1和DU3每个计算周期平均功率变化示意图;
图10是本发明优选实施例的牵引逆变器U相上开关管DU12和DU34每个计算周期平均功率变化示意图;
图11是本发明优选实施例的牵引逆变器U相上开关管TU1和TU2的结温波动示意图;
图12是本发明优选实施例的牵引逆变器U相上开关管TU3和TU4的结温波动示意图;
图13是本发明优选实施例的牵引逆变器U相上开关管DU1和DU3的结温波动示意图;
图14是本发明优选实施例的牵引逆变器U相上开关管DU12和DU34的结温波动示意图。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以由权利要求限定和覆盖的多种不同方式实施。
实施例1
参见图1,本实施例提供一种牵引变流器器件结温在线计算方法,包括以下步骤:
针对牵引变流器目标桥臂上的目标器件,建立目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型;目标器件为二极管或开关管,功率损耗包括通断功耗和通态功耗;
考虑电热耦合,建立目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型;
在非第一周期的其他周期内,根据第一、第二及第三关系模型计算目标器件在该周期的平均功率,并计算该周期所对应的目标器件到所在外壳之间的第一等效热阻抗、以及目标器件由外壳到外部环境网络之间的第二等效热阻抗;然后根据平均功率、第一等效阻抗及第二等效阻抗计算目标器件在该周期的结温值。
需要说明的是,牵引变流器的目标桥臂,是指脉冲整流器A相、B相,以及牵引逆变器U相、V相、W相五相桥臂中任一相,桥臂上的目标器件,是指牵引变流器中的开关管和二极管,其中每一相桥臂由上往下由四个开关管组成,分别记为Tx1,Tx2,Tx3和Tx4,x=A,B,U,V,W,且每一相桥臂由上往下由四个续流二极管和两个钳位二极管组成,分别记为Dx1,Dx2,Dx3,Dx4和Dx12,Dx34。
优选地,本实施例参照CRH2型高速列车的三电平牵引变流器系统为仿真对象,其主电路拓扑结构如图2所示。并以牵引逆变器U相桥臂上的开关管和二极管的结温在线计算为例,仿真时长为1s,其中,三电平牵引变流器系统正常运行时。主要仿真参数如表1所示。
表1 三电平牵引变流器系统正常运行时主要仿真参数表
参数 | 数值 |
直流环节电压给定值 | 1300V |
给定转速 | 200km/h |
整流器开关周期 | 0.8ms |
逆变器开关周期 | 1ms |
控制器控制周期 | 40μs |
牵引电机极对数 | 2 |
输出功率 | 300kW |
转动惯量 | 0.4kg·m<sup>2</sup> |
首先,针对牵引变流器中牵引逆变器的U相桥臂,建立各开关管功耗与交流侧电流和各开关管驱动控制信号之间的关系模型,公式为:
式中,EIU1表示逆变器U相桥臂上开关管TU1的功耗,EIU2表示逆变器U相桥臂上开关管TU2的功耗,EIU3表示逆变器U相桥臂上开关管TU3的功耗,EIU4表示逆变器U相桥臂上开关管TU4的功耗,其中,表示桥臂交流侧电流极性的标志位,iU表示牵引逆变器U桥臂交流侧电流,在本实施例中该电流为牵引电机U相定子电流isa,sU1、sU2、sU3和sU4依次表示牵引逆变器U相桥臂上四个开关管的开关状态,受系统中控制器输出的开关管驱动控制信号控制,为1时代表导通状态,为0时代表关断状态;-、∧和∨分别表示取反、与和或运算;uce_U1表示逆变器U相桥臂上开关管TU1的通态压降,uce_U2表示逆变器U相桥臂上开关管TU2的通态压降,uce_U3表示逆变器U相桥臂上开关管TU3的通态压降,uce_U4表示逆变器U相桥臂上开关管TU4的通态压降,t和t0均表示时间,其中,t-t0到t表示一个控制周期的时间长度;γIU1表示逆变器U相桥臂上开关管TU1的通断功耗,γIU2表示逆变器U相桥臂上开关管TU2的通断功耗,γIU3表示逆变器U相桥臂上开关管TU3的通断功耗,γIU4表示逆变器U相桥臂上开关管TU4的通断功耗。
针对牵引变流器中牵引逆变器的U相桥臂,建立二极管功率损耗与交流侧电流和驱动控制信号之间的关系模型,公式为:
式中,EDU1表示逆变器U相桥臂上二极管DU1的功耗,EDU2表示逆变器U相桥臂上二极管DU2的功耗,EDU3表示逆变器U相桥臂上二极管DU3的功耗,EDU4表示逆变器U相桥臂上二极管DU4的功耗,EDU12表示逆变器U相桥臂上二极管DU12的功耗,EDU34表示逆变器U相桥臂上二极管DU34的功耗,ufd_U1表示逆变器U相桥臂上二极管DU1的通态压降,ufd_U3表示逆变器U相桥臂上二极管DU3的通态压降,ufd_U12表示逆变器U相桥臂上二极管DU12的通态压降,ufd_U34表示逆变器U相桥臂上二极管DU34的通态压降,γDU1表示逆变器U相桥臂上二极管DU1的通断功耗,γDU2表示逆变器U相桥臂上二极管DU1的通断功耗,γDU12表示逆变器U相桥臂上二极管DU12的通断功耗,γDU34表示逆变器U相桥臂上二极管DU34的通断功耗;
假定在一个控制周期Tctrl=40μs内(Tctrl远小于桥臂交流侧电流的响应时间常数),牵引变流器U相桥臂交流侧电流iU和开关管与二极管的通态压降近似为恒值,则公式(1)和(2)分别近似为公式(3)和(4):
进一步地,牵引变流器U相桥臂上开关管的通断功耗的计算公式为:
式中,Eon_IU1表示逆变器U相桥臂上开关管TU1的导通功耗,Eon_IU2表示逆变器U相桥臂上开关管TU2的导通功耗,Eon_IU3表示逆变器U相桥臂上开关管TU3的导通功耗,Eon_IU4表示逆变器U相桥臂上开关管TU4的导通功耗,Eoff_IU1表示逆变器U相桥臂上开关管TU1的关断功耗,Eoff_IU2表示逆变器U相桥臂上开关管TU2的关断功耗,Eoff_IU3表示逆变器U相桥臂上开关管TU3的关断功耗,Eoff_IU4表示逆变器U相桥臂上开关管TU4的关断功耗。
牵引变流器U相桥臂上二极管通断功耗的计算公式为:
式中,Erec_DU1表示逆变器U相桥臂上二极管DU1的导通功耗,Erec_DU3表示逆变器U相桥臂上二极管DU3的导通功耗,Erec_DU12表示逆变器U相桥臂上二极管DU12的导通功耗,Erec_DU34表示逆变器U相桥臂上二极管DU34的导通功耗。
需要说明的是,开关管的通断功耗包括导通功耗和关断功耗,且开关管的通态功耗受初始通态压降和通态等效电阻影响;二极管的通断功耗为反向恢复损耗,且二极管的通态功耗受初始通态压降和通态等效电阻影响;在实际的工作中,由于温度(热场)对电气模型中电气参数存在一定的影响,即,存在热敏效应,且根据温度的大小发生变化,而电气参数能对器件的功耗造成影响,所以,为了使计算结果更为准确,应当考虑考虑电热耦合,因此可以建立逆变器U相桥臂上的各开关管和二极管在当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型为:
式中,Tj_IUi表示逆变器U相桥臂上第i个开关管的芯片结温,uceo_Ui表示该开关管的初始通态压降,和rIon_Ui表示该开关管的通态等效电阻,i表示某相桥臂上开关管的编号,i=1,2,3,4;aceo_Ui和cceo_Ui表示开关管初始通态压降的拟合系数,aIon_Ui和cIon_Ui表示开关管通态等效电阻的拟合系数,表示开关管在结温为125℃时的导通功耗,表示开关管在结温为25℃时的导通功耗,表示开关管在结温为125℃时的关断功耗,表示开关管在结温为25℃时的关断功耗,表示开关管在结温为125℃时导通功耗的拟合系数,表示开关管在结温为25℃时导通功耗的拟合系数,表示开关管在结温为125℃时关断功耗的拟合系数, 表示开关管在结温为25℃时关断功耗的拟合系数。其中,开关管的各拟合系数由厂商用户数据手册中提供的数据拟合获得,具体数值参见如下表2所示。
表2 开关管的拟合系数的数值表
式中,Tj_DUk表示二极管的芯片结温,ufdo_Uk表示二极管的初始通态压降,ron_Uk表示二极管的通态等效电阻,k表示桥臂上开关管的编号,k=1,2,3,4,12,34,afdo_Uk和cfdo_Uk表示二极管初始通态压降的拟合系数,aon_Uk和con_Uk表示二极管通态等效电阻的拟合系数,表示二极管在结温为125℃时的反向恢复功耗,表示二极管在结温为25℃时的反向恢复功耗,表示二极管在结温为125℃时导通功耗的拟合系数,表示二极管在结温为25℃时导通功耗的拟合系数,式中各拟合系数由厂商用户数据手册中提供的数据拟合获得,具体数值参见如下表3所示。
表3 二极管的拟合系数的数值表
进一步地,获取牵引变流器系统中交流侧电流传感器的电流信号和各目标器件的驱动控制信号,以1/n个开关周期作为平均功率的计算周期计算一个计算周期内开关管的总功耗,得到开关管的平均功率,具体地,计算一个计算周期内开关管的总功耗,得到开关管的平均功率,公式为:
式中,Pavg_IUi表示一个计算周期Tcal内开关管的平均功率,Etot_IUi表示一个计算周期Tcal内开关管的总功耗,EIUi[q]表示一个计算周期Tcal中第q个控制周期Tctrl内开关管的功耗。
计算一个计算周期内二极管的总功耗,得到二极管的平均功率,公式为:
式中,Pavg_DUk表示一个计算周期Tcal内二极管的平均功率,Etot_DUk表示一个计算周期Tcal内二极管的总功耗,EDUk[q]表示一个计算周期Tcal中第q个控制周期Tctrl内二极管的功耗。
然后,计算逆变器U相桥臂上的开关管芯片到所在外壳之间的第一等效热阻抗为:
式中,Zjc_IUi表示开关管芯片到器件所在模块外壳之间的第一等效热阻抗,n表示芯-壳之间等效热阻抗的阶数,Rpjc_IUi表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻,Cpjc_IUi表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容,τpjc_DUk表示开关管芯-壳之间等效热网络的热时间常数。
本实施例中,以系统最小的开关周期0.8ms的1/5作为平均功率的计算周期Tcal,Tcal=0.8/5=0.16ms。
建立二极管芯片到所在外壳之间的第一等效热阻抗为:
式中,Zjc_DUk表示二极管芯片到器件所在模块外壳之间的第一等效热阻抗,n表示芯-壳之间等效热阻抗的阶数,Rpjc_DUk表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻,Cpjc_DUk表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容,τpjc_Dxk表示二极管芯-壳之间等效热网络的热时间常数;
建立开关管由外壳到外部环境网络之间的第二等效热阻抗,公式为:
式中,Zch_IUi表示开关管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗,m表示壳-散-环之间等效热阻抗的阶数,Rpch_IUi表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻,Cpch_IUi表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容,τpch_IUi表示开关管壳-散-环之间等效热网络的热时间常数;
建立二极管由外壳到外部环境网络之间的第二等效热阻抗,公式为:
式中,Zch_DUk表示二极管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗,m表示壳-散-环之间等效热阻抗的阶数,Rpch_DUk表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻,Cpch_DUk表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容,τpch_DUk表示二极管壳-散-环之间等效热网络的热时间常数。需要说明的是,上述第一等效热阻抗和第二等效热阻抗统称为开关管和二极管的等效热网络模型。
本实施例中,每相桥臂上的关键器件均同一厂商的产品,桥臂上不同位置的开关管或二极管芯-壳之间等效热网络模型的参数相同,具体赋值如表4所示。
表4 桥臂上各器件的等效热网络模型参数
进一步地,计算所述目标器件在该周期的结温值具体包括以下步骤:
在一个控制周期Tctrl内,计算开关管的结温,公式为:
式中,Tamb表示环境温度;
在一个控制周期Tctrl内,计算二极管的结温,公式为:
将所述开关管的结温和二极管的结温分别代入公式(13)-公式(17)中,以计算下一个平均功率的计算周期内牵引变流器关键器件的结温值。并根据各结温值得到器件的结温波动。
值得说明的是,第一周期的结温值根据统计经验设为默认值。基于该设定的默认值完成实时的牵引变流器器件结温在线计算。本实施例中,在第一个控制周期Tctrl内计算关键器件的结温中,关键器件平均功率的初始数值为Pavg0_Ixi和Pavg0_Dxk,结温的初始数值为Tj0_Ixi和Tj0_Dxk,环境温度的初始数值为Tamb0。作为可变换的实时方式,将结温的初始数值Tj0_IUi和Tj0_DUk都设为25℃,环境温度的初始数值为Tamb0为25℃牵引逆变器U相上开关管和二极管每个计算周期内功耗的变化,如图3至图6所示;牵引逆变器U相上开关管和二极管每个计算周期平均功率变化,如图7至图10所示;牵引逆变器U相上开关管和二极管的结温波动,如图11至图14所示。
作为本实施例优选的实施方式,等效热阻抗在用于计算器件结温前,针对连续控制系统需要进行拉普拉斯变换,针对离散控制系统需要进行Z变换。
作为本实施例优选的实施方式,在一个控制周期Tctrl内计算关键器件的结温中,关键器件的平均功率的数值通过计算所得,其中,零阶保持器的输出周期为一个计算周期Tcal。
实施例2
与上述方法实施例相对应地,本实施例提供一种牵引变流器器件结温在线计算系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述方法的步骤。
如上所述,本发明提供的一种牵引变流器器件结温在线计算方法及系统,实现了计及电热耦合下牵引变流器关键器件结温波动的在线监测,能在线计算牵引变流器中目标器件的结温变化。
进一步地,本发明仅利用系统中已有的电流传感器数据和各开关管的驱动控制信号,不受系统调制策略和恒频工作条件的约束,能实现变频和不同调制策略条件下牵引变流器关键器件结温的在线计算,提高了器件结温在线计算的灵活性,为牵引变流器系统的健康监测与管理提供了有力的科学支撑。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (9)
1.一种牵引变流器器件结温在线计算方法,其特征在于,包括以下步骤:
针对牵引变流器目标桥臂上的目标器件,建立所述目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型;所述目标器件为二极管或开关管,所述功率损耗包括通断功耗和通态功耗;
考虑电热耦合,建立所述目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型;
在非第一周期的其他周期内,根据所述第一、第二及第三关系模型计算所述目标器件在该周期的平均功率,并计算该周期所对应的所述目标器件到所在外壳之间的第一等效热阻抗、以及所述目标器件由外壳到外部环境网络之间的第二等效热阻抗;
根据所述平均功率、第一等效阻抗及第二等效阻抗计算所述目标器件在该周期的结温值。
2.根据权利要求1所述的牵引变流器器件结温在线计算方法,其特征在于,所述建立所述目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型具体包括以下步骤:
(1)建立目标桥臂上的各开关管的功率损耗与交流侧电流及驱动控制信号之间的关系模型,公式为:
式中,EIx1表示桥臂上开关管Tx1的功耗,EIx2表示桥臂上开关管Tx2的功耗,EIx3表示桥臂上开关管Tx3的功耗,和EIx4表示桥臂上开关管Tx4的功耗,其中,x表示所在桥臂,x=A,B,U,V,W;表示桥臂交流侧电流极性的标志位,ix表示桥臂交流侧电流,sx1表示桥臂上开关管Tx1的开关状态,sx2表示桥臂上开关管Tx2的开关状态,sx3表示桥臂上开关管Tx3的开关状态,sx4表示桥臂上开关管Tx4的开关状态;uce_x1表示桥臂上开关管Tx1的通态压降,uce_x2表示桥臂上开关管Tx2的通态压降,uce_x3表示桥臂上开关管Tx3的通态压降,和uce_x4表示桥臂上开关管Tx4的通态压降,t和t0表示时间;γIx1表示桥臂上开关管Tx1的通断功耗,γIx2表示桥臂上开关管Tx2的通断功耗,γIx3表示桥臂上开关管Tx3的通断功耗,γIx4表示桥臂上开关管Tx4的通断功耗;
(2)建立目标桥臂上的各二极管的功率损耗与交流侧电流及驱动控制信号之间的关系模型,公式为:
式中,EDx1表示桥臂上二极管Dx1的功耗,EDx2表示桥臂上二极管Dx2的功耗,EDx3表示桥臂上二极管Dx3的功耗,EDx4表示桥臂上二极管Dx4的功耗,EDx12表示桥臂上二极管Dx12的功耗,EDx34表示桥臂上二极管Dx34的功耗,ufd_x1表示桥臂上二极管Dx1的通态压降,ufd_x3表示桥臂上二极管Dx3的通态压降,ufd_x12表示桥臂上二极管Dx12的通态压降,ufd_x34表示桥臂上二极管Dx34的通态压降,γDx1表示桥臂上二极管Dx1的通断功耗,γDx3表示桥臂上二极管Dx3的通断功耗,γDx12表示桥臂上二极管Dx12的通断功耗,γDx34表示桥臂上二极管Dx34的通断功耗;
(3)假定在一个控制周期Tctrl内,桥臂交流侧电流ix和开关管与二极管的通态压降近似为恒值,则公式(1)和(2)分别近似为公式(3)和(4):
3.根据权利要求2所述的牵引变流器器件结温在线计算方法,其特征在于,所述桥臂上开关管的通断功耗的计算公式为:
式中,Eon_Ix1表示桥臂上开关管Tx1的导通功耗,Eon_Ix2表示桥臂上开关管Tx2的导通功耗,Eon_Ix3表示桥臂上开关管Tx3的导通功耗,Eon_Ix4表示桥臂上开关管Tx4的导通功耗,Eoff_Ix1表示桥臂上开关管Tx1的关断功耗,Eoff_Ix2表示桥臂上开关管Tx2的关断功耗,Eoff_Ix3表示桥臂上开关管Tx3的关断功耗,Eoff_Ix4表示桥臂上开关管Tx4的关断功耗;
所述桥臂上二极管通断功耗的计算公式为:
式中,Erec_Dx1表示桥臂上二极管Dx1的导通功耗,Erec_Dx3表示桥臂上二极管Dx3的导通功耗,Erec_Dx12表示桥臂上二极管Dx12的导通功耗,Erec_Dx34表示桥臂上二极管Dx34的导通功耗。
4.根据权利要求3所述的牵引变流器器件结温在线计算方法,其特征在于,开关管的通断功耗包括导通功耗和关断功耗,且开关管的通态功耗受初始通态压降和通态等效电阻影响;二极管的通断功耗为反向恢复损耗,且二极管的通态功耗受初始通态压降和通态等效电阻影响;所述建立所述目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型具体包括以下步骤:
式中,Tj_Ixi表示目标桥臂上第i个开关管的芯片结温,uceo_xi表示该开关管的初始通态压降,rIon_xi表示该开关管的通态等效电阻,i表示目标桥臂上开关管的编号,i=1,2,3,4,aceo_xi和cceo_xi表示开关管初始通态压降的拟合系数,aIon_xi和cIon_xi表示开关管通态等效电阻的拟合系数,表示开关管在结温为125℃时的导通功耗,表示开关管在结温为25℃时的导通功耗,表示开关管在结温为125℃时的关断功耗,表示开关管在结温为25℃时的关断功耗,表示开关管在结温为125℃时导通功耗的拟合系数,表示开关管在结温为25℃时导通功耗的拟合系数, 表示开关管在结温为125℃时关断功耗的拟合系数,表示开关管在结温为25℃时关断功耗的拟合系数;
式中,Tj_Dxk表示二极管的芯片结温,ufdo_xk表示二极管的初始通态压降,ron_xk表示二极管的初始通态等效电阻,k表示目标桥臂上开关管的编号,k=1,2,3,4,12,34,afdo_xk和cfdo_xk表示二极管初始通态压降的拟合系数,aon_xk和con_xk表示二极管通态等效电阻的拟合系数,表示二极管在结温为125℃时的反向恢复功耗,表示二极管在结温为25℃时的反向恢复功耗,表示二极管在结温为125℃时导通功耗的拟合系数,表示二极管在结温为25℃时导通功耗的拟合系数。
5.根据权利要求4所述的牵引变流器器件结温在线计算方法,其特征在于,所述根据所述第一、第二及第三关系模型计算所述目标器件在该周期的平均功率具体包括以下步骤:
获取牵引变流器系统中交流侧电流传感器的电流信号和各目标器件的驱动控制信号,以1/n个开关周期作为平均功率的计算周期计算一个计算周期内开关管的总功耗,得到开关管的平均功率,公式为:
式中,Pavg_Ixi表示一个计算周期Tcal内开关管的平均功率,Etot_Ixi表示一个计算周期Tcal内开关管的总功耗,EIxi[q]表示一个计算周期Tcal中第q个控制周期Tctrl内开关管的功耗;
计算一个计算周期内二极管的总功耗,得到二极管的平均功率,公式为:
式中,Pavg_Dxk表示一个计算周期Tcal内二极管的平均功率,Etot_Dxk表示一个计算周期Tcal内二极管的总功耗,EDxk[q]表示一个计算周期Tcal中第q个控制周期Tctrl内二极管的功耗。
6.根据权利要求5所述的牵引变流器器件结温在线计算方法,其特征在于,所述计算所述目标器件到所在外壳之间的第一等效热阻抗、以及所述目标器件由外壳到外部环境网络之间的第二等效热阻抗具体包括以下步骤:
建立开关管芯片到所在外壳之间的第一等效热阻抗,公式为:
式中,Zjc_Ixi表示开关管芯片到器件所在模块外壳之间的第一等效热阻抗,n表示芯-壳之间等效热阻抗的阶数,Rpjc_Ixi表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻,Cpjc_Ixi表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容,τpjc_Ixi表示开关管芯-壳之间等效热网络的热时间常数;
建立二极管芯片到所在外壳之间的第一等效热阻抗,公式为:
式中,Zjc_Dxk表示二极管芯片到器件所在模块外壳之间的第一等效热阻抗,n表示芯-壳之间等效热阻抗的阶数,Rpjc_Dxk表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻,Cpjc_Dxk表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容,τpjc_Dxk表示二极管芯-壳之间等效热网络的热时间常数;
建立开关管由外壳到外部环境网络之间的第二等效热阻抗,公式为:
式中,Zch_Ixi表示开关管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗,m表示壳-散-环之间等效热阻抗的阶数,Rpch_Ixi表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻,Cpch_Ixi表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容,τpch_Ixi表示开关管壳-散-环之间等效热网络的热时间常数;
建立二极管由外壳到外部环境网络之间的第二等效热阻抗,公式为:
式中,Zch_Dxk表示二极管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗,m表示壳-散-环之间等效热阻抗的阶数,Rpch_Dxk表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻,Cpch_Dxk表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容,τpch_Dxk表示二极管壳-散-环之间等效热网络的热时间常数。
7.根据权利要求6所述的牵引变流器器件结温在线计算方法,其特征在于,所述根据所述平均功率、第一等效阻抗及第二等效阻抗计算所述目标器件在该周期的结温值具体包括以下步骤:
在一个控制周期Tctrl内,计算开关管的结温,公式为:
Tj_Ixi=Tamb+(Zjc_Ixi+Zch_Ixi)*Pavg_Ixi (24)
式中,Tamb表示环境温度;
在一个控制周期Tctrl内,计算二极管的结温,公式为:
Tj_Dxk=Tamb+(Zjc_Dxk+Zch_Dxk)*Pavg_Dxk (25)
将所述开关管的结温和二极管的结温分别代入公式(13)-公式(17)中,以计算下一个平均功率的计算周期内牵引变流器关键器件的结温。
8.根据权利要求7所述的牵引变流器器件结温在线计算方法,其特征在于,所述第一周期的结温值根据统计经验设为默认值。
9.一种牵引变流器器件结温在线计算方法系统,包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现上述权利要求1至8任一所述方法的步骤。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810490961.3A CN108649807B (zh) | 2018-05-21 | 2018-05-21 | 一种牵引变流器器件结温在线计算方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810490961.3A CN108649807B (zh) | 2018-05-21 | 2018-05-21 | 一种牵引变流器器件结温在线计算方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108649807A CN108649807A (zh) | 2018-10-12 |
CN108649807B true CN108649807B (zh) | 2019-08-30 |
Family
ID=63757396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810490961.3A Active CN108649807B (zh) | 2018-05-21 | 2018-05-21 | 一种牵引变流器器件结温在线计算方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108649807B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110456264B (zh) * | 2019-08-16 | 2021-07-30 | 中国电子产品可靠性与环境试验研究所((工业和信息化部电子第五研究所)(中国赛宝实验室)) | 大功率器件功率循环试验结温监测方法、装置和系统 |
CN110795852B (zh) * | 2019-10-31 | 2023-05-09 | 兰州交通大学 | 考虑电热耦合影响的机车端电压预估方法 |
CN111371301B (zh) * | 2020-03-13 | 2021-11-02 | 中南大学 | 一种两电平牵引逆变器igbt结温控制方法及系统 |
CN113759226B (zh) * | 2020-06-01 | 2024-04-12 | 安徽威灵汽车部件有限公司 | 车辆、控制器及其功率器件结温估计方法和装置 |
CN113036725B (zh) * | 2021-04-14 | 2022-08-23 | 苏州汇川技术有限公司 | 伺服驱动器保护方法、装置、设备及存储介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103956887A (zh) * | 2014-05-15 | 2014-07-30 | 重庆大学 | 风电变流器igbt模块结温在线计算方法 |
CN106443400A (zh) * | 2016-09-14 | 2017-02-22 | 河北工业大学 | 一种igbt模块的电‑热‑老化结温计算模型建立方法 |
-
2018
- 2018-05-21 CN CN201810490961.3A patent/CN108649807B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103956887A (zh) * | 2014-05-15 | 2014-07-30 | 重庆大学 | 风电变流器igbt模块结温在线计算方法 |
CN106443400A (zh) * | 2016-09-14 | 2017-02-22 | 河北工业大学 | 一种igbt模块的电‑热‑老化结温计算模型建立方法 |
Non-Patent Citations (1)
Title |
---|
IGBT junction temperature measurement via combined TSEPs with collector current impact elimination;Xiang Wang等;《IEEE》;20170216;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN108649807A (zh) | 2018-10-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108649807B (zh) | 一种牵引变流器器件结温在线计算方法及系统 | |
US11539214B2 (en) | Modular power conversion system | |
CN103956887B (zh) | 风电变流器igbt模块结温在线计算方法 | |
CN105337523B (zh) | Npc型三电平逆变器可靠性在线监测系统及方法 | |
CN103472331B (zh) | 一种基于光伏物理模型的光伏发电故障诊断系统 | |
ES2859732T3 (es) | Inversores eólicos-solares híbridos | |
RU2534457C2 (ru) | Способ и устройство контроля показателей энергоэффективности оборудования | |
CN104020801A (zh) | 一种温度控制方法,温度计算方法及驱动器 | |
CN105974247A (zh) | 导电滑环的测试系统及方法 | |
CN107425533A (zh) | 一种电压不平衡时同步逆变器并网运行态势利导方法 | |
CN106019044B (zh) | 电力牵引交流传动两电平单相四象限脉冲整流器开关管开路故障诊断方法 | |
CN103399219B (zh) | 光伏发电站性能实时监测方法 | |
US20160369761A1 (en) | Direct Current Power Plant | |
BRPI0511697B1 (pt) | sistema de propulsão de baixa perda | |
CN110456201B (zh) | Mmc子模块的功率循环加速试验装置及其控制方法 | |
CN105093021A (zh) | 模块化多电平换流器用桥臂电抗器的温升试验方法 | |
CN112305350B (zh) | 一种新能源机组调频测试平台功率模块的故障检测装置 | |
US20240093663A1 (en) | Direct current power plant | |
CN206506299U (zh) | 基于mcu的svg与电容投切的集中控制电路 | |
KR101478658B1 (ko) | Pcm 진단장치 | |
CN206712669U (zh) | 变频器节能电路和变频节能系统 | |
CN205809616U (zh) | 分布式多路电缆热循环测控系统 | |
CN207516492U (zh) | 一种中压xlpe电缆及其附件多谐波源绝缘试验装置 | |
CN212675382U (zh) | 一种发电机组制冷中控系统 | |
CN218882575U (zh) | 一种通用的风扇测试系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |