CN113341293B - 一种基于参数辨识的地铁牵引整流器二极管结温提取方法 - Google Patents
一种基于参数辨识的地铁牵引整流器二极管结温提取方法 Download PDFInfo
- Publication number
- CN113341293B CN113341293B CN202110639218.1A CN202110639218A CN113341293B CN 113341293 B CN113341293 B CN 113341293B CN 202110639218 A CN202110639218 A CN 202110639218A CN 113341293 B CN113341293 B CN 113341293B
- Authority
- CN
- China
- Prior art keywords
- diode
- thermal
- temperature
- rectifier
- radiator
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
- G01R31/26—Testing of individual semiconductor devices
- G01R31/2607—Circuits therefor
- G01R31/2632—Circuits therefor for testing diodes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Algebra (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Testing Of Individual Semiconductor Devices (AREA)
Abstract
本发明公开了一种基于参数辨识的地铁牵引整流器二极管结温提取方法,具体为:采集地铁牵引变电所整流器二极管正向导通电流,散热器温度和整流柜外环境温度;根据二极管产品数据手册提供的正向伏安特性曲线,求得二极管各时刻正向导通电压;根据瞬态热阻抗曲线,拟合初始状态下二极管热阻,二极管热容,散热器热阻和散热器热容,并计算初始状态下的辨识量;计算单个采样周期内的二极管导通损耗;使用最小二乘法辨识得到t时刻的辨识量,修正t时刻的热网络模型参数,并以此计算二极管结温。本发明利用整流器温度保护中散热器温度的监控数据,对热网络模型参数进行辨识,得到二极管当前状态下的热网络模型参数,提高了结温计算的精度。
Description
技术领域
本发明属于二极管结温计算技术领域,尤其涉及一种基于参数辨识的地铁牵引整流器二极管结温提取方法。
背景技术
根据电力电子系统可靠性调研,功率器件是变流系统中失效率最高的部件,约占34%,在各类失效因素中,约55%的电力电子系统失效主要由温度因素诱发。由于地铁运行中频繁地启动与制动,导致整流器输出电流波动快、摆幅大,从而造成其二极管结温波动剧烈;同时,地铁整流柜柜内无风机,装置的冷却方式为空气自然对流换热,散热条件较差。这些因素都会引起整流器二极管的疲劳损伤,甚至导致二极管的开、短路故障。因此有必要结合地铁牵引整流器实际运行工况,研究准确有效的二极管结温提取方法,对二极管进行状态监测。
由于在实际中二极管会随着服役时间增加而发生老化,进而导致其热阻增大。若一直使用二极管产品数据手册所提供的热阻进行结温计算,将会导致结果偏小,从而无法正确反应二极管当前的状态。因此需要对二极管热阻进行实时修正,提高结温计算的精度。
发明内容
为能辨识出二极管老化造成的热阻变化,修正二极管当前状态下的热网络模型参数,提高结温计算精度。本发明提供一种基于参数辨识的地铁牵引整流器二极管结温提取方法。
本发明的一种基于参数辨识的地铁牵引整流器二极管结温提取方法,包括以下步骤:
步骤1:以Δt为采样间隔,采集地铁牵引变电所整流器二极管正向导通电流iF(t)、散热器温度Th(t)和整流柜外环境温度Ta(t)。
步骤2:根据二极管产品数据手册提供的正向伏安特性曲线,求得二极管正向导通电压iF(t)所对应的正向导通电压uF(t);根据二极管产品数据手册提供的瞬态热阻抗曲线,拟合得到初始状态下二极管热阻R1(0),二极管热容C1(0),散热器热阻Rh(0),散热器热容Ch(0);根据式(1)计算得到初始状态下的辨识量θ1(0)、θ2(0)、θ3(0)和θ4(0)。
步骤3:根据式(2)计算t时刻的二极管导通损耗P(t)。
P(t)=iF(t)uF(t) (2)
步骤4:使用0到t时刻二极管导通损耗、散热器温度和整流柜外环境温度数据,根据式(3),利用最小二乘法辨识得到t时刻的辨识量θ1(t)、θ2(t)、θ3(t)和θ4(t)。
Th(t)=θ1(t)Th(t-1)+θ2(t)Th(t-2)+θ3(t)P(t-1)+θ4(t)Ta(t-1) (3)
步骤5:根据式(4)计算t时刻二极管热阻改变量系数g1(t)和散热器热阻改变量系数gh(t)。
步骤6:根据式(5)计算t时刻的热网络模型参数:二极管热阻R1(t),二极管热容C1(t),散热器热阻Rh(t),散热器热容Ch(t)。
步骤7:根据t时刻的热网络模型参数计算二极管结温,具体如式(6):
进一步的,采样间隔Δt为0.1s。
本发明的有益技术效果为:
本发明利用整流器温度保护中散热器温度的监控数据,对热网络模型参数进行辨识,得到二极管当前状态下的热网络模型参数,提高了结温计算的精度。
附图说明
图1为本发明适用的地铁牵引整流器二极管二阶连续热网络模型。
图2为仿真实验中二极管正向导通电流信号曲线图。
图3为仿真实验中散热器温度信号曲线图。
图4为仿真实验中二极管导通损耗曲线图。
图5为仿真实验中二极管结温曲线图。
具体实施方式
下面结合附图和仿真实验对本发明作进一步的详细说明。
本发明的一种基于参数辨识的地铁牵引整流器二极管结温提取方法,包括以下步骤:
步骤1:以Δt为采样间隔,采集地铁牵引变电所整流器二极管正向导通电流iF(t)、散热器温度Th(t)和整流柜外环境温度Ta(t)。
步骤2:根据二极管产品数据手册提供的正向伏安特性曲线,求得二极管正向导通电压iF(t)所对应的正向导通电压uF(t);根据二极管产品数据手册提供的瞬态热阻抗曲线,拟合得到初始状态下二极管热阻R1(0),二极管热容C1(0),散热器热阻Rh(0),散热器热容Ch(0);根据式(1)计算得到初始状态下的辨识量θ1(0)、θ2(0)、θ3(0)和θ4(0)。
步骤3:根据式(2)计算t时刻的二极管导通损耗P(t)。
P(t)=iF(t)uF(t) (2)
步骤4:使用0到t时刻二极管导通损耗、散热器温度和整流柜外环境温度数据,根据式(3),利用最小二乘法辨识得到t时刻的辨识量θ1(t)、θ2(t)、θ3(t)和θ4(t)。
Th(t)=θ1(t)Th(t-1)+θ2(t)Th(t-2)+θ3(t)P(t-1)+θ4(t)Ta(t-1) (3)
步骤5:根据式(4)计算t时刻二极管热阻改变量系数g1(t)和散热器热阻改变量系数gh(t)。
步骤6:根据式(5)计算t时刻的热网络模型参数:二极管热阻R1(t),二极管热容C1(t),散热器热阻Rh(t),散热器热容Ch(t)。
步骤7:根据t时刻的热网络模型参数计算二极管结温,具体如式(6):
进一步的,采样间隔Δt为0.1s。
仿真实验:
为验证本发明方法的有效性,进行以下仿真实验。
利用PLECS软件建立如图1所示的地铁牵引整流器二极管二阶连续热网络模型。
以0.1s为采样间隔,采集地铁牵引变电所整流器二极管正向导通电流iF(t),散热器温度Th(t),分别如图2和图3所示。假设整流柜外环境温度Ta(t)为恒定值25℃。
根据产品数据手册提供的瞬态热阻抗曲线,拟合得到获取初始状态下二极管热阻R1(0)=0.1K/W,二极管热容C1(0)=10J/K,散热器热阻Rh(0)=0.2K/W,散热器热容Ch(0)=20J/K。计算得到初始状态下的辨识量θ1(0)=1.8489,θ2(0)=-0.8511,θ3(0)=4.2553×10-4,θ4(0)=2.1277×10-3。
计算单个采样周期内的二极管导通损耗P(t)如图4所示。使用0到t时刻二极管导通损耗、散热器温度和整流柜外环境温度数据,利用最小二乘法辨识得到t时刻的辨识量θ1(t)=1.8532,θ2(t)=-0.8551、θ3(t)=4.0997×10-4,θ4(t)=1.9300×10-3。
计算t时刻二极管热阻改变量系数g1(t)=1.1068,散热器热阻改变量系数gh(t)=1.1075.
计算t时刻的热网络模型参数R1(t)=0.11068K/W,二极管热容C1(t)=10J/K,散热器热阻Rh(t)=0.2215K/W,散热器热容Ch(t)=20J/K。
根据t时刻的热网络模型参数计算二极管结温如图5所示。
Claims (2)
1.一种基于参数辨识的地铁牵引整流器二极管结温提取方法,其特征在于,包括以下步骤:
步骤1:以Δt为采样间隔,采集地铁牵引变电所整流器二极管正向导通电流iF(t)、散热器温度Th(t)和整流柜外环境温度Ta(t);
步骤2:根据二极管产品数据手册提供的正向伏安特性曲线,求得二极管正向导通电流iF(t)所对应的正向导通电压uF(t);根据二极管产品数据手册提供的瞬态热阻抗曲线,拟合得到初始状态下二极管热阻R1(0),二极管热容C1(0),散热器热阻Rh(0),散热器热容Ch(0);根据式(1)计算得到初始状态下的辨识量θ1(0)、θ2(0)、θ3(0)和θ4(0);
步骤3:根据式(2)计算t时刻的二极管导通损耗P(t):
P(t)=iF(t)uF(t) (2)
步骤4:使用0到t时刻二极管导通损耗、散热器温度和整流柜外环境温度数据,根据式(3),利用最小二乘法辨识得到t时刻的辨识量θ1(t)、θ2(t)、θ3(t)和θ4(t);
Th(t)=θ1(t)Th(t-1)+θ2(t)Th(t-2)+θ3(t)P(t-1)+θ4(t)Ta(t-1) (3)
步骤5:根据式(4)计算t时刻二极管热阻改变量系数g1(t)和散热器热阻改变量系数gh(t);
步骤6:根据式(5)计算t时刻的热网络模型参数:二极管热阻R1(t),二极管热容C1(t),散热器热阻Rh(t),散热器热容Ch(t);
步骤7:根据t时刻的热网络模型参数计算二极管结温,具体如式(6):
2.根据权利要求1所述的一种基于参数辨识的地铁牵引整流器二极管结温提取方法,其特征在于,所述采样间隔Δt为0.1s。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110639218.1A CN113341293B (zh) | 2021-06-08 | 2021-06-08 | 一种基于参数辨识的地铁牵引整流器二极管结温提取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110639218.1A CN113341293B (zh) | 2021-06-08 | 2021-06-08 | 一种基于参数辨识的地铁牵引整流器二极管结温提取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113341293A CN113341293A (zh) | 2021-09-03 |
CN113341293B true CN113341293B (zh) | 2022-09-09 |
Family
ID=77475396
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110639218.1A Active CN113341293B (zh) | 2021-06-08 | 2021-06-08 | 一种基于参数辨识的地铁牵引整流器二极管结温提取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113341293B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101490921A (zh) * | 2006-05-22 | 2009-07-22 | Rt专利有限公司 | 用于使车辆交流发电机的输出最大化的系统和方法 |
CN103234656A (zh) * | 2013-05-17 | 2013-08-07 | 厦门大学 | 一种发光二极管结温的测量方法 |
CN203395974U (zh) * | 2013-07-17 | 2014-01-15 | 张冠夫 | 一种水暖型电暖气 |
CN104091203A (zh) * | 2014-07-24 | 2014-10-08 | 重庆大学 | 一种风力发电用变流器的实时可靠性评估方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110295558A1 (en) * | 2010-05-27 | 2011-12-01 | General Electric Company | Method and system for detecting faults in a brushless exciter for a generator |
-
2021
- 2021-06-08 CN CN202110639218.1A patent/CN113341293B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101490921A (zh) * | 2006-05-22 | 2009-07-22 | Rt专利有限公司 | 用于使车辆交流发电机的输出最大化的系统和方法 |
CN103234656A (zh) * | 2013-05-17 | 2013-08-07 | 厦门大学 | 一种发光二极管结温的测量方法 |
CN203395974U (zh) * | 2013-07-17 | 2014-01-15 | 张冠夫 | 一种水暖型电暖气 |
CN104091203A (zh) * | 2014-07-24 | 2014-10-08 | 重庆大学 | 一种风力发电用变流器的实时可靠性评估方法 |
Non-Patent Citations (4)
Title |
---|
Study on Stray Current Distribution of the Metro Depot;Liping Zhao 等;《Proceedings of the 4th International Conference on Electrical and Information Technologies for Rail Transportation (EITRT) 2019 》;20200408;第377–385页 * |
地铁整流器二极管配置及均流均压方法;孙才勤 等;《电气化铁道》;20080831(第4期);第42-44+47页 * |
大功率二极管(晶闸管)瞬态电流持续时的瞬态结温评定;F.Profumo 等;《变流技术与电力牵引》;20000430(第2期);第20-26页 * |
大型发电机整流器励磁系统;H.J.Humphries等;《大电机技术》;19731231(第03期);第94-104+115页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113341293A (zh) | 2021-09-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102608511B (zh) | 一种金属氧化物半导体管的结温和热阻测量方法 | |
CN109738773A (zh) | 一种非平稳工况下igbt模块寿命预测方法 | |
CN103956887B (zh) | 风电变流器igbt模块结温在线计算方法 | |
CN105242189A (zh) | 基于集射极饱和压降与焊料层空洞率的igbt健康状态监测方法 | |
CN112084651A (zh) | 计及疲劳损伤的多尺度风电igbt可靠性评估方法及系统 | |
CN107621599A (zh) | 一种测量igbt在高温反偏试验中结温变化的方法 | |
CN111060797A (zh) | 一种基于热网络自然频率的igbt模块健康状态监测方法 | |
CN113219315B (zh) | 基于结温归一化的逆变器igbt老化在线监测方法及系统 | |
Yang et al. | A temperature-dependent Cauer model simulation of IGBT module with analytical thermal impedance characterization | |
CN113341293B (zh) | 一种基于参数辨识的地铁牵引整流器二极管结温提取方法 | |
CN110514978A (zh) | 补偿变流器可靠性在线监测方法 | |
CN111090940B (zh) | 一种基于ansys的mmc子模块压接式igbt短期失效分析方法 | |
US20230089918A1 (en) | Method and apparatus for controlling charging, based on monitored communication signals associated with a charging session | |
CN115753880B (zh) | 一种基于综合温升因子的油浸式车载牵引变压器散热性能的评估方法 | |
CN107703432A (zh) | 一种电力电子开关器件结温在线检测方法及检测电路 | |
CN110008488A (zh) | 一种提取晶体管的通用型热阻和热容特性的方法 | |
CN113343463B (zh) | 考虑老化进程的地铁牵引整流器二极管剩余寿命预测方法 | |
CN105717382A (zh) | 一种变压器使用寿命预估系统 | |
CN113420407B (zh) | 一种igct水冷散热器建模及其结温计算方法 | |
CN110197012B (zh) | 考虑牵引传动系统故障影响的支撑电容寿命评估方法 | |
CN113987747A (zh) | 一种考虑键合线断裂的igbt模块老化失效分析方法 | |
CN104218898B (zh) | 基于多参数物理模型的三相程控功率源过热保护方法 | |
CN112858378A (zh) | 一种集成电路散热性能测试方法 | |
CN113987998B (zh) | 基于电力工况的全控型功率半导体模块功率循环方法 | |
CN108647453A (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 |