CN108649807B - 一种牵引变流器器件结温在线计算方法及系统 - Google Patents

Abstract

本发明涉及电力电子技术领域，公开了一种牵引变流器器件结温在线计算方法及系统，以实现三电平牵引变流器变频调速系统中关键器件结温波动的在线计算，进一步检测器件的健康，保证系统的稳定性和鲁棒性；针对牵引变流器目标桥臂上的目标器件，建立目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型；建立目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型；在非第一周期的其他周期内，计算该周期所对应的目标器件到所在外壳之间的第一等效热阻抗、以及目标器件由外壳到外部环境网络之间的第二等效热阻抗；然后计算目标器件在该周期的结温值。

Description

一种牵引变流器器件结温在线计算方法及系统

技术领域

本发明涉及电力电子技术领域，尤其涉及一种牵引变流器器件结温在线计算方法及系统。

背景技术

牵引传动系统被誉为轨道交通装备的“心脏”，不仅是整个轨道交通装备的核心动力单元，更是其运行安全的关键系统之一，对牵引传动系统中关键器件的健康监测至关重要。牵引传动系统由牵引变压器、牵引变流器等主电路设备，以及牵引电机等装置构成，其中，牵引变流器包括脉冲整流器、中间直流环节以及牵引逆变器。因此，牵引变流器是这枚“心脏”的“起搏器”，由于长期在高温、高压等恶劣的工作环境下，牵引变流器是牵引传动系统的高发故障源。牵引变流器中关键器件，主要是开关管和二极管的性能退化带来的安全隐患若不能及时监测到并得到及时正确的处理，都有可能引发器件的损坏失效，进而导致系统发生连锁事故。因此，保证牵引变流器及其关键器件安全可靠地工作不仅是确保牵引传动系统安全运行的关键，也是确保轨道交通装备安全运行的关键。

由开关管和二极管的失效机理分析可知，55％以上的失效原因与器件的结温波动有关，即，器件结温的频繁波动将加速器件键合线脱落和焊料层出现空穴等老化问题，进而引起器件的性能退化，最后导致器件的损坏失效。因此，对关键器件进行结温波动的在线监测将为器件剩余寿命、老化程度的估计和失效机理的分析提供了重要的数据支撑。目前，直接在线测量开关管和二极管的结温的方案大多需要增加额外的测温传感器或模块，虽然精度较高，但成本也较高，且不符合当下及未来工业界对关键器件轻量化和集成化设计的需求。因此，建立系统中关键器件的功耗和结温等模型，通过系统中已有的传感器模块，间接评估器件结温波动成为一种有效的途径。但目前的研究成果主要为两电平的变流器拓扑结构，且集中于恒频工况和采用正弦脉宽调制(SinusoidalPulse Width Modulation，SPWM)策略的应用环境，如风力发电，电网能量传输等，难以适用于三电平，变频调速工况和采用空间矢量脉宽调制(Space Vector Pulse Width Modulation，SVPWM)策略的应用环境，如CRH2型高速列车牵引传动系统即为变频调速系统，采用的三电平牵引变流器拓扑，其中，脉冲整流器采用SPWM调制策略，牵引逆变器采用SVPWM调制策略。

因此，现需提供一种能实现三电平牵引变流器变频调速系统中关键器件结温波动的在线计算与监测的牵引变流器器件结温在线计算方法及系统。

发明内容

本发明目的在于提供一种牵引变流器器件结温在线计算方法及系统，以实现三电平牵引变流器变频调速系统中关键器件结温波动的在线计算，进一步检测器件的健康，保证系统的稳定性和鲁棒性。

为实现上述目的，本发明提供了一种牵引变流器器件结温在线计算方法，包括以下步骤：

针对牵引变流器目标桥臂上的目标器件，建立所述目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型；所述目标器件为二极管或开关管，所述功率损耗包括通断功耗和通态功耗；

考虑电热耦合，建立所述目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型；

在非第一周期的其他周期内，根据所述第一、第二及第三关系模型计算所述目标器件在该周期的平均功率，并计算该周期所对应的所述目标器件到所在外壳之间的第一等效热阻抗、以及所述目标器件由外壳到外部环境网络之间的第二等效热阻抗；然后根据所述平均功率、第一等效阻抗及第二等效阻抗计算所述目标器件在该周期的结温值。

与上述方法相对应地，本发明还提供一种牵引变流器器件结温在线计算系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法的步骤。

本发明具有以下有益效果：

1、本发明提供的一种牵引变流器器件结温在线计算方法及系统，实现了计及电热耦合下牵引变流器关键器件结温波动的在线监测，能在线计算牵引变流器中目标器件的结温变化。

2、本发明提供的一种牵引变流器器件结温在线计算方法及系统，仅利用系统中已有的电流传感器数据和各开关管的驱动控制信号，不受系统调制策略和恒频工作条件的约束，能实现变频和不同调制策略条件下牵引变流器关键器件结温的在线计算，提高了器件结温在线计算的灵活性，为牵引变流器系统的健康监测与管理提供了有力的科学支撑。

下面将参照附图，对本发明作进一步详细的说明。

附图说明

构成本申请的一部分的附图用来提供对本发明的进一步理解，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：

图1是本发明优选实施例的牵引变流器器件结温在线计算方法的流程图；

图2是本发明优选实施例的三电平牵引变流器系统的主电路拓扑结构图；

图3是本发明优选实施例的牵引逆变器U相上开关管T_U1和T_U2每个计算周期内功耗的变化示意图；

图4是本发明优选实施例的牵引逆变器U相上开关管T_U3和T_U4每个计算周期内功耗的变化示意图；

图5是本发明优选实施例的牵引逆变器U相上二极管D_U1和D_U3每个计算周期内功耗的变化示意图；

图6是本发明优选实施例的牵引逆变器U相上二极管D_U12和D_U34每个计算周期内功耗的变化示意图；

图7是本发明优选实施例的牵引逆变器U相上开关管T_U1和T_U2每个计算周期平均功率变化示意图；

图8是本发明优选实施例的牵引逆变器U相上开关管T_U3和T_U4每个计算周期平均功率变化示意图；

图9是本发明优选实施例的牵引逆变器U相上开关管D_U1和D_U3每个计算周期平均功率变化示意图；

图10是本发明优选实施例的牵引逆变器U相上开关管D_U12和D_U34每个计算周期平均功率变化示意图；

图11是本发明优选实施例的牵引逆变器U相上开关管T_U1和T_U2的结温波动示意图；

图12是本发明优选实施例的牵引逆变器U相上开关管T_U3和T_U4的结温波动示意图；

图13是本发明优选实施例的牵引逆变器U相上开关管D_U1和D_U3的结温波动示意图；

图14是本发明优选实施例的牵引逆变器U相上开关管D_U12和D_U34的结温波动示意图。

具体实施方式

以下结合附图对本发明的实施例进行详细说明，但是本发明可以由权利要求限定和覆盖的多种不同方式实施。

实施例1

参见图1，本实施例提供一种牵引变流器器件结温在线计算方法，包括以下步骤：

针对牵引变流器目标桥臂上的目标器件，建立目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型；目标器件为二极管或开关管，功率损耗包括通断功耗和通态功耗；

考虑电热耦合，建立目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型；

在非第一周期的其他周期内，根据第一、第二及第三关系模型计算目标器件在该周期的平均功率，并计算该周期所对应的目标器件到所在外壳之间的第一等效热阻抗、以及目标器件由外壳到外部环境网络之间的第二等效热阻抗；然后根据平均功率、第一等效阻抗及第二等效阻抗计算目标器件在该周期的结温值。

需要说明的是，牵引变流器的目标桥臂，是指脉冲整流器A相、B相，以及牵引逆变器U相、V相、W相五相桥臂中任一相，桥臂上的目标器件，是指牵引变流器中的开关管和二极管，其中每一相桥臂由上往下由四个开关管组成，分别记为T_x1，T_x2，T_x3和T_x4，x＝A,B,U,V,W，且每一相桥臂由上往下由四个续流二极管和两个钳位二极管组成，分别记为D_x1，D_x2，D_x3，D_x4和D_x12，D_x34。

优选地，本实施例参照CRH2型高速列车的三电平牵引变流器系统为仿真对象，其主电路拓扑结构如图2所示。并以牵引逆变器U相桥臂上的开关管和二极管的结温在线计算为例，仿真时长为1s，其中，三电平牵引变流器系统正常运行时。主要仿真参数如表1所示。

表1 三电平牵引变流器系统正常运行时主要仿真参数表

参数	数值
		直流环节电压给定值	1300V
给定转速	200km/h
		整流器开关周期	0.8ms
逆变器开关周期	1ms
		控制器控制周期	40μs
牵引电机极对数	2
		输出功率	300kW
转动惯量	0.4kg·m<sup>2</sup>

首先，针对牵引变流器中牵引逆变器的U相桥臂，建立各开关管功耗与交流侧电流和各开关管驱动控制信号之间的关系模型，公式为：

式中，E_IU1表示逆变器U相桥臂上开关管T_U1的功耗，E_IU2表示逆变器U相桥臂上开关管T_U2的功耗，E_IU3表示逆变器U相桥臂上开关管T_U3的功耗，E_IU4表示逆变器U相桥臂上开关管T_U4的功耗，其中，表示桥臂交流侧电流极性的标志位，i_U表示牵引逆变器U桥臂交流侧电流，在本实施例中该电流为牵引电机U相定子电流i_sa，s_U1、s_U2、s_U3和s_U4依次表示牵引逆变器U相桥臂上四个开关管的开关状态，受系统中控制器输出的开关管驱动控制信号控制，为1时代表导通状态，为0时代表关断状态；-、∧和∨分别表示取反、与和或运算；u_{ce_U1}表示逆变器U相桥臂上开关管T_U1的通态压降，u_{ce_U2}表示逆变器U相桥臂上开关管T_U2的通态压降，u_{ce_U3}表示逆变器U相桥臂上开关管T_U3的通态压降，u_{ce_U4}表示逆变器U相桥臂上开关管T_U4的通态压降，t和t₀均表示时间，其中，t-t₀到t表示一个控制周期的时间长度；γ_IU1表示逆变器U相桥臂上开关管T_U1的通断功耗，γ_IU2表示逆变器U相桥臂上开关管T_U2的通断功耗，γ_IU3表示逆变器U相桥臂上开关管T_U3的通断功耗，γ_IU4表示逆变器U相桥臂上开关管T_U4的通断功耗。

针对牵引变流器中牵引逆变器的U相桥臂，建立二极管功率损耗与交流侧电流和驱动控制信号之间的关系模型，公式为：

式中，E_DU1表示逆变器U相桥臂上二极管D_U1的功耗，E_DU2表示逆变器U相桥臂上二极管D_U2的功耗，E_DU3表示逆变器U相桥臂上二极管D_U3的功耗，E_DU4表示逆变器U相桥臂上二极管D_U4的功耗，E_DU12表示逆变器U相桥臂上二极管D_U12的功耗，E_DU34表示逆变器U相桥臂上二极管D_U34的功耗，u_{fd_U1}表示逆变器U相桥臂上二极管D_U1的通态压降，u_{fd_U3}表示逆变器U相桥臂上二极管D_U3的通态压降，u_{fd_U12}表示逆变器U相桥臂上二极管D_U12的通态压降，u_{fd_U34}表示逆变器U相桥臂上二极管D_U34的通态压降，γ_DU1表示逆变器U相桥臂上二极管D_U1的通断功耗，γ_DU2表示逆变器U相桥臂上二极管D_U1的通断功耗，γ_DU12表示逆变器U相桥臂上二极管D_U12的通断功耗，γ_DU34表示逆变器U相桥臂上二极管D_U34的通断功耗；

假定在一个控制周期T_ctrl＝40μs内(T_ctrl远小于桥臂交流侧电流的响应时间常数)，牵引变流器U相桥臂交流侧电流i_U和开关管与二极管的通态压降近似为恒值，则公式(1)和(2)分别近似为公式(3)和(4)：

进一步地，牵引变流器U相桥臂上开关管的通断功耗的计算公式为：

式中，E_{on_IU1}表示逆变器U相桥臂上开关管T_U1的导通功耗，E_{on_IU2}表示逆变器U相桥臂上开关管T_U2的导通功耗，E_{on_IU3}表示逆变器U相桥臂上开关管T_U3的导通功耗，E_{on_IU4}表示逆变器U相桥臂上开关管T_U4的导通功耗，E_{off_IU1}表示逆变器U相桥臂上开关管T_U1的关断功耗，E_{off_IU2}表示逆变器U相桥臂上开关管T_U2的关断功耗，E_{off_IU3}表示逆变器U相桥臂上开关管T_U3的关断功耗，E_{off_IU4}表示逆变器U相桥臂上开关管T_U4的关断功耗。

牵引变流器U相桥臂上二极管通断功耗的计算公式为：

式中，E_{rec_DU1}表示逆变器U相桥臂上二极管D_U1的导通功耗，E_{rec_DU3}表示逆变器U相桥臂上二极管D_U3的导通功耗，E_{rec_DU12}表示逆变器U相桥臂上二极管D_U12的导通功耗，E_{rec_DU34}表示逆变器U相桥臂上二极管D_U34的导通功耗。

需要说明的是，开关管的通断功耗包括导通功耗和关断功耗，且开关管的通态功耗受初始通态压降和通态等效电阻影响；二极管的通断功耗为反向恢复损耗，且二极管的通态功耗受初始通态压降和通态等效电阻影响；在实际的工作中，由于温度(热场)对电气模型中电气参数存在一定的影响，即，存在热敏效应，且根据温度的大小发生变化，而电气参数能对器件的功耗造成影响，所以，为了使计算结果更为准确，应当考虑考虑电热耦合，因此可以建立逆变器U相桥臂上的各开关管和二极管在当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型为：

式中，T_{j_IUi}表示逆变器U相桥臂上第i个开关管的芯片结温，u_{ceo_Ui}表示该开关管的初始通态压降，和r_{Ion_Ui}表示该开关管的通态等效电阻，i表示某相桥臂上开关管的编号，i＝1,2,3,4；a_{ceo_Ui}和c_{ceo_Ui}表示开关管初始通态压降的拟合系数，a_{Ion_Ui}和c_{Ion_Ui}表示开关管通态等效电阻的拟合系数，表示开关管在结温为125℃时的导通功耗，表示开关管在结温为25℃时的导通功耗，表示开关管在结温为125℃时的关断功耗，表示开关管在结温为25℃时的关断功耗，表示开关管在结温为125℃时导通功耗的拟合系数，表示开关管在结温为25℃时导通功耗的拟合系数，表示开关管在结温为125℃时关断功耗的拟合系数， 表示开关管在结温为25℃时关断功耗的拟合系数。其中，开关管的各拟合系数由厂商用户数据手册中提供的数据拟合获得，具体数值参见如下表2所示。

表2 开关管的拟合系数的数值表

式中，T_{j_DUk}表示二极管的芯片结温，u_{fdo_Uk}表示二极管的初始通态压降，r_{on_Uk}表示二极管的通态等效电阻，k表示桥臂上开关管的编号，k＝1,2,3,4,12,34，a_{fdo_Uk}和c_{fdo_Uk}表示二极管初始通态压降的拟合系数，a_{on_Uk}和c_{on_Uk}表示二极管通态等效电阻的拟合系数，表示二极管在结温为125℃时的反向恢复功耗，表示二极管在结温为25℃时的反向恢复功耗，表示二极管在结温为125℃时导通功耗的拟合系数，表示二极管在结温为25℃时导通功耗的拟合系数，式中各拟合系数由厂商用户数据手册中提供的数据拟合获得，具体数值参见如下表3所示。

表3 二极管的拟合系数的数值表

进一步地，获取牵引变流器系统中交流侧电流传感器的电流信号和各目标器件的驱动控制信号，以1/n个开关周期作为平均功率的计算周期计算一个计算周期内开关管的总功耗，得到开关管的平均功率，具体地，计算一个计算周期内开关管的总功耗，得到开关管的平均功率，公式为：

式中，P_{avg_IUi}表示一个计算周期T_cal内开关管的平均功率，E_{tot_IUi}表示一个计算周期T_cal内开关管的总功耗，E_IUi[q]表示一个计算周期T_cal中第q个控制周期T_ctrl内开关管的功耗。

计算一个计算周期内二极管的总功耗，得到二极管的平均功率，公式为：

式中，P_{avg_DUk}表示一个计算周期T_cal内二极管的平均功率，E_{tot_DUk}表示一个计算周期T_cal内二极管的总功耗，E_DUk[q]表示一个计算周期T_cal中第q个控制周期T_ctrl内二极管的功耗。

然后，计算逆变器U相桥臂上的开关管芯片到所在外壳之间的第一等效热阻抗为：

式中，Z_{jc_IUi}表示开关管芯片到器件所在模块外壳之间的第一等效热阻抗，n表示芯-壳之间等效热阻抗的阶数，R_{pjc_IUi}表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻，C_{pjc_IUi}表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容，τ_{pjc_DUk}表示开关管芯-壳之间等效热网络的热时间常数。

本实施例中，以系统最小的开关周期0.8ms的1/5作为平均功率的计算周期T_cal，T_cal＝0.8/5＝0.16ms。

建立二极管芯片到所在外壳之间的第一等效热阻抗为：

式中，Z_{jc_DUk}表示二极管芯片到器件所在模块外壳之间的第一等效热阻抗，n表示芯-壳之间等效热阻抗的阶数，R_{pjc_DUk}表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻，C_{pjc_DUk}表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容，τ_{pjc_Dxk}表示二极管芯-壳之间等效热网络的热时间常数；

建立开关管由外壳到外部环境网络之间的第二等效热阻抗，公式为：

式中，Z_{ch_IUi}表示开关管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗，m表示壳-散-环之间等效热阻抗的阶数，R_{pch_IUi}表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻，C_{pch_IUi}表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容，τ_{pch_IUi}表示开关管壳-散-环之间等效热网络的热时间常数；

建立二极管由外壳到外部环境网络之间的第二等效热阻抗，公式为：

式中，Z_{ch_DUk}表示二极管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗，m表示壳-散-环之间等效热阻抗的阶数，R_{pch_DUk}表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻，C_{pch_DUk}表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容，τ_{pch_DUk}表示二极管壳-散-环之间等效热网络的热时间常数。需要说明的是，上述第一等效热阻抗和第二等效热阻抗统称为开关管和二极管的等效热网络模型。

本实施例中，每相桥臂上的关键器件均同一厂商的产品，桥臂上不同位置的开关管或二极管芯-壳之间等效热网络模型的参数相同，具体赋值如表4所示。

表4 桥臂上各器件的等效热网络模型参数

进一步地，计算所述目标器件在该周期的结温值具体包括以下步骤：

在一个控制周期T_ctrl内，计算开关管的结温，公式为：

式中，T_amb表示环境温度；

在一个控制周期T_ctrl内，计算二极管的结温，公式为：

将所述开关管的结温和二极管的结温分别代入公式(13)-公式(17)中，以计算下一个平均功率的计算周期内牵引变流器关键器件的结温值。并根据各结温值得到器件的结温波动。

值得说明的是，第一周期的结温值根据统计经验设为默认值。基于该设定的默认值完成实时的牵引变流器器件结温在线计算。本实施例中，在第一个控制周期T_ctrl内计算关键器件的结温中，关键器件平均功率的初始数值为P_{avg0_Ixi}和P_{avg0_Dxk}，结温的初始数值为T_{j0_Ixi}和T_{j0_Dxk}，环境温度的初始数值为T_amb0。作为可变换的实时方式，将结温的初始数值T_{j0_IUi}和T_{j0_DUk}都设为25℃，环境温度的初始数值为T_amb0为25℃牵引逆变器U相上开关管和二极管每个计算周期内功耗的变化，如图3至图6所示；牵引逆变器U相上开关管和二极管每个计算周期平均功率变化，如图7至图10所示；牵引逆变器U相上开关管和二极管的结温波动，如图11至图14所示。

作为本实施例优选的实施方式，等效热阻抗在用于计算器件结温前，针对连续控制系统需要进行拉普拉斯变换，针对离散控制系统需要进行Z变换。

作为本实施例优选的实施方式，在一个控制周期T_ctrl内计算关键器件的结温中，关键器件的平均功率的数值通过计算所得，其中，零阶保持器的输出周期为一个计算周期T_cal。

实施例2

与上述方法实施例相对应地，本实施例提供一种牵引变流器器件结温在线计算系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述方法的步骤。

如上所述，本发明提供的一种牵引变流器器件结温在线计算方法及系统，实现了计及电热耦合下牵引变流器关键器件结温波动的在线监测，能在线计算牵引变流器中目标器件的结温变化。

进一步地，本发明仅利用系统中已有的电流传感器数据和各开关管的驱动控制信号，不受系统调制策略和恒频工作条件的约束，能实现变频和不同调制策略条件下牵引变流器关键器件结温的在线计算，提高了器件结温在线计算的灵活性，为牵引变流器系统的健康监测与管理提供了有力的科学支撑。

以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种牵引变流器器件结温在线计算方法，其特征在于，包括以下步骤：

针对牵引变流器目标桥臂上的目标器件，建立所述目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型；所述目标器件为二极管或开关管，所述功率损耗包括通断功耗和通态功耗；

考虑电热耦合，建立所述目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型；

在非第一周期的其他周期内，根据所述第一、第二及第三关系模型计算所述目标器件在该周期的平均功率，并计算该周期所对应的所述目标器件到所在外壳之间的第一等效热阻抗、以及所述目标器件由外壳到外部环境网络之间的第二等效热阻抗；

根据所述平均功率、第一等效阻抗及第二等效阻抗计算所述目标器件在该周期的结温值。

2.根据权利要求1所述的牵引变流器器件结温在线计算方法，其特征在于，所述建立所述目标器件功率损耗与交流侧电流及驱动控制信号之间的第一关系模型具体包括以下步骤：

(1)建立目标桥臂上的各开关管的功率损耗与交流侧电流及驱动控制信号之间的关系模型，公式为：

式中，E_Ix1表示桥臂上开关管T_x1的功耗，E_Ix2表示桥臂上开关管T_x2的功耗，E_Ix3表示桥臂上开关管T_x3的功耗，和E_Ix4表示桥臂上开关管T_x4的功耗，其中，x表示所在桥臂，x＝A,B,U,V,W；表示桥臂交流侧电流极性的标志位，i_x表示桥臂交流侧电流，s_x1表示桥臂上开关管T_x1的开关状态，s_x2表示桥臂上开关管T_x2的开关状态，s_x3表示桥臂上开关管T_x3的开关状态，s_x4表示桥臂上开关管T_x4的开关状态；u_{ce_x1}表示桥臂上开关管T_x1的通态压降，u_{ce_x2}表示桥臂上开关管T_x2的通态压降，u_{ce_x3}表示桥臂上开关管T_x3的通态压降，和u_{ce_x4}表示桥臂上开关管T_x4的通态压降，t和t₀表示时间；γ_Ix1表示桥臂上开关管T_x1的通断功耗，γ_Ix2表示桥臂上开关管T_x2的通断功耗，γ_Ix3表示桥臂上开关管T_x3的通断功耗，γ_Ix4表示桥臂上开关管T_x4的通断功耗；

(2)建立目标桥臂上的各二极管的功率损耗与交流侧电流及驱动控制信号之间的关系模型，公式为：

式中，E_Dx1表示桥臂上二极管D_x1的功耗，E_Dx2表示桥臂上二极管D_x2的功耗，E_Dx3表示桥臂上二极管D_x3的功耗，E_Dx4表示桥臂上二极管D_x4的功耗，E_Dx12表示桥臂上二极管D_x12的功耗，E_Dx34表示桥臂上二极管D_x34的功耗，u_{fd_x1}表示桥臂上二极管D_x1的通态压降，u_{fd_x3}表示桥臂上二极管D_x3的通态压降，u_{fd_x12}表示桥臂上二极管D_x12的通态压降，u_{fd_x34}表示桥臂上二极管D_x34的通态压降，γ_Dx1表示桥臂上二极管D_x1的通断功耗，γ_Dx3表示桥臂上二极管D_x3的通断功耗，γ_Dx12表示桥臂上二极管D_x12的通断功耗，γ_Dx34表示桥臂上二极管D_x34的通断功耗；

(3)假定在一个控制周期T_ctrl内，桥臂交流侧电流i_x和开关管与二极管的通态压降近似为恒值，则公式(1)和(2)分别近似为公式(3)和(4)：

3.根据权利要求2所述的牵引变流器器件结温在线计算方法，其特征在于，所述桥臂上开关管的通断功耗的计算公式为：

式中，E_{on_Ix1}表示桥臂上开关管T_x1的导通功耗，E_{on_Ix2}表示桥臂上开关管T_x2的导通功耗，E_{on_Ix3}表示桥臂上开关管T_x3的导通功耗，E_{on_Ix4}表示桥臂上开关管T_x4的导通功耗，E_{off_Ix1}表示桥臂上开关管T_x1的关断功耗，E_{off_Ix2}表示桥臂上开关管T_x2的关断功耗，E_{off_Ix3}表示桥臂上开关管T_x3的关断功耗，E_{off_Ix4}表示桥臂上开关管T_x4的关断功耗；

所述桥臂上二极管通断功耗的计算公式为：

式中，E_{rec_Dx1}表示桥臂上二极管D_x1的导通功耗，E_{rec_Dx3}表示桥臂上二极管D_x3的导通功耗，E_{rec_Dx12}表示桥臂上二极管D_x12的导通功耗，E_{rec_Dx34}表示桥臂上二极管D_x34的导通功耗。

4.根据权利要求3所述的牵引变流器器件结温在线计算方法，其特征在于，开关管的通断功耗包括导通功耗和关断功耗，且开关管的通态功耗受初始通态压降和通态等效电阻影响；二极管的通断功耗为反向恢复损耗，且二极管的通态功耗受初始通态压降和通态等效电阻影响；所述建立所述目标器件的当前周期的通断功耗与上一周期结温之间的第二关系模型、以及所述目标器件当前周期的通态功耗与上一周期结温之间的第三关系模型具体包括以下步骤：

式中，T_{j_Ixi}表示目标桥臂上第i个开关管的芯片结温，u_{ceo_xi}表示该开关管的初始通态压降，r_{Ion_xi}表示该开关管的通态等效电阻，i表示目标桥臂上开关管的编号，i＝1,2,3,4，a_{ceo_xi}和c_{ceo_xi}表示开关管初始通态压降的拟合系数，a_{Ion_xi}和c_{Ion_xi}表示开关管通态等效电阻的拟合系数，表示开关管在结温为125℃时的导通功耗，表示开关管在结温为25℃时的导通功耗，表示开关管在结温为125℃时的关断功耗，表示开关管在结温为25℃时的关断功耗，表示开关管在结温为125℃时导通功耗的拟合系数，表示开关管在结温为25℃时导通功耗的拟合系数， 表示开关管在结温为125℃时关断功耗的拟合系数，表示开关管在结温为25℃时关断功耗的拟合系数；

式中，T_{j_Dxk}表示二极管的芯片结温，u_{fdo_xk}表示二极管的初始通态压降，r_{on_xk}表示二极管的初始通态等效电阻，k表示目标桥臂上开关管的编号，k＝1,2,3,4,12,34，a_{fdo_xk}和c_{fdo_xk}表示二极管初始通态压降的拟合系数，a_{on_xk}和c_{on_xk}表示二极管通态等效电阻的拟合系数，表示二极管在结温为125℃时的反向恢复功耗，表示二极管在结温为25℃时的反向恢复功耗，表示二极管在结温为125℃时导通功耗的拟合系数，表示二极管在结温为25℃时导通功耗的拟合系数。

5.根据权利要求4所述的牵引变流器器件结温在线计算方法，其特征在于，所述根据所述第一、第二及第三关系模型计算所述目标器件在该周期的平均功率具体包括以下步骤：

获取牵引变流器系统中交流侧电流传感器的电流信号和各目标器件的驱动控制信号，以1/n个开关周期作为平均功率的计算周期计算一个计算周期内开关管的总功耗，得到开关管的平均功率，公式为：

式中，P_{avg_Ixi}表示一个计算周期T_cal内开关管的平均功率，E_{tot_Ixi}表示一个计算周期T_cal内开关管的总功耗，E_Ixi[q]表示一个计算周期T_cal中第q个控制周期T_ctrl内开关管的功耗；

计算一个计算周期内二极管的总功耗，得到二极管的平均功率，公式为：

式中，P_{avg_Dxk}表示一个计算周期T_cal内二极管的平均功率，E_{tot_Dxk}表示一个计算周期T_cal内二极管的总功耗，E_Dxk[q]表示一个计算周期T_cal中第q个控制周期T_ctrl内二极管的功耗。

6.根据权利要求5所述的牵引变流器器件结温在线计算方法，其特征在于，所述计算所述目标器件到所在外壳之间的第一等效热阻抗、以及所述目标器件由外壳到外部环境网络之间的第二等效热阻抗具体包括以下步骤:

建立开关管芯片到所在外壳之间的第一等效热阻抗，公式为：

式中，Z_{jc_Ixi}表示开关管芯片到器件所在模块外壳之间的第一等效热阻抗，n表示芯-壳之间等效热阻抗的阶数，R_{pjc_Ixi}表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻，C_{pjc_Ixi}表示开关管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容，τ_{pjc_Ixi}表示开关管芯-壳之间等效热网络的热时间常数；

建立二极管芯片到所在外壳之间的第一等效热阻抗，公式为：

式中，Z_{jc_Dxk}表示二极管芯片到器件所在模块外壳之间的第一等效热阻抗，n表示芯-壳之间等效热阻抗的阶数，R_{pjc_Dxk}表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热阻，C_{pjc_Dxk}表示二极管芯-壳之间等效网络的第pjc阶等效热阻抗的等效热容，τ_{pjc_Dxk}表示二极管芯-壳之间等效热网络的热时间常数；

建立开关管由外壳到外部环境网络之间的第二等效热阻抗，公式为：

式中，Z_{ch_Ixi}表示开关管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗，m表示壳-散-环之间等效热阻抗的阶数，R_{pch_Ixi}表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻，C_{pch_Ixi}表示开关管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容，τ_{pch_Ixi}表示开关管壳-散-环之间等效热网络的热时间常数；

建立二极管由外壳到外部环境网络之间的第二等效热阻抗，公式为：

式中，Z_{ch_Dxk}表示二极管所在模块外壳到散热器再到外部环境之间的第二等效热阻抗，m表示壳-散-环之间等效热阻抗的阶数，R_{pch_Dxk}表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热阻，C_{pch_Dxk}表示二极管壳-散-环之间等效网络的第pch阶等效热阻抗的等效热容，τ_{pch_Dxk}表示二极管壳-散-环之间等效热网络的热时间常数。

7.根据权利要求6所述的牵引变流器器件结温在线计算方法，其特征在于，所述根据所述平均功率、第一等效阻抗及第二等效阻抗计算所述目标器件在该周期的结温值具体包括以下步骤：

在一个控制周期T_ctrl内，计算开关管的结温，公式为：

T_{j_Ixi}＝T_amb+(Z_{jc_Ixi}+Z_{ch_Ixi})*P_{avg_Ixi} (24)

式中，T_amb表示环境温度；

在一个控制周期T_ctrl内，计算二极管的结温，公式为：

T_{j_Dxk}＝T_amb+(Z_{jc_Dxk}+Z_{ch_Dxk})*P_{avg_Dxk} (25)

将所述开关管的结温和二极管的结温分别代入公式(13)-公式(17)中，以计算下一个平均功率的计算周期内牵引变流器关键器件的结温。

8.根据权利要求7所述的牵引变流器器件结温在线计算方法，其特征在于，所述第一周期的结温值根据统计经验设为默认值。

9.一种牵引变流器器件结温在线计算方法系统，包括存储器、处理器以及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现上述权利要求1至8任一所述方法的步骤。