发明名称 电力牵引交流传动两电平单相四象限脉冲整流器故障建模仿真方法
摘要 本发明电力牵引交流传动两电平单相四象限脉冲整流器故障建模仿真方法,基于两电平单相四象限脉冲整流器的工作原理完成,具体方法包括:电力牵引交流传动两电平单相四象限脉冲整流器在每一种故障模式下所对应的开关函数计算以及基于开关函数计算结果的电力牵引交流传动两电平单相四象限脉冲整流器状态方程计算。本方法可实现离线仿真以及在线实时仿真系统中,两电平单相四象限脉冲整流器在正常工作状态下以及不同开关管故障状态下的仿真,并且基于状态方程的模型可以实现正常工作状态以及不同开关管故障状态下的切换,弥补了现有仿真模型中无故障仿真模型的技术问题。
申请公布号 CN103425826A 申请公布日期 2013.12.04
申请号 CN201310335516.7 申请日期 2013.08.05
申请人 西南交通大学 发明人 葛兴来;韩坤;冯晓云;熊成林;宋文胜;苟斌;崔恒斌
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 电力牵引交流传动两电平单相四象限脉冲整流器故障建模仿真方法,在Matlab/Simulink下进行两电平四象限脉冲整流器正常工作以及故障模式下的仿真,定义单相两电平脉冲整流器左边桥臂为A桥臂,右边桥臂为B桥臂,包括以下步骤:(1)、根据输入脉冲计算不同故障模式下开关状态函数:不同故障模式下开关状态函数计算通过电力牵引交流传动两电平单相四象限脉冲整流器的开关元器件控制脉冲P1、P2、P3、P4完成;不同故障模式下开关状态函数计算包括11种情况,定义SA、SB分别为A、B桥臂开关函数,在不同的脉冲整流器交流侧电流is下,每种情况下的开关状态函数计算方法如下:(a)正常工作情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(b)A桥臂上IGBT故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(c)A桥臂下IGBT故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(d)B桥臂上IGBT故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(e)B桥臂下IGBT故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(f)A桥臂上IGBT与A桥臂下IGBT同时故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(g)A桥臂上IGBT与B桥臂上IGBT同时故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(h)A桥臂上IGBT与B桥臂下IGBT同时故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(i)A桥臂下IGBT与B桥臂上IGBT同时故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mrow> <mn>1</mn> <mo>-</mo> <mi>P</mi> </mrow> <mn>4</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(j)A桥臂下IGBT与B桥臂下IGBT同时故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <msub> <mi>P</mi> <mn>3</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>(k)B桥臂上IGBT与B桥臂下IGBT同时故障情况的开关函数计算方法: <mrow> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>1</mn> <mo>-</mo> <msub> <mi>P</mi> <mn>2</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <msub> <mi>P</mi> <mn>1</mn> </msub> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow> <mrow> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>=</mo> <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>></mo> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>1</mn> </mtd> <mtd> <msub> <mi>i</mi> <mi>s</mi> </msub> <mo>&lt;</mo> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>其中is为两电平单相四象限脉冲整流器交流侧的电流值。(2)、将(1)所得开关状态函数,计算不同情况下的两电平单相四象限脉冲整流器状态方程:将(1)所得的开关函数值输入如下方程: <mfenced open='{' close=''> <mtable> <mtr> <mtd> <mi>R</mi> <mo>=</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>S</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>R</mi> <mi>p</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>ab</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>-</mo> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>U</mi> <mi>d</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>L</mi> <mi>s</mi> </msub> <mfrac> <msub> <mi>di</mi> <mi>s</mi> </msub> <mi>dt</mi> </mfrac> <mo>=</mo> <msub> <mi>u</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>u</mi> <mi>ab</mi> </msub> <mo>-</mo> <msub> <mi>Ri</mi> <mi>s</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>d</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <msub> <mi>S</mi> <mi>A</mi> </msub> <mo>-</mo> <msub> <mi>S</mi> <mi>B</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>i</mi> <mi>s</mi> </msub> </mtd> </mtr> </mtable> </mfenced>计算获得:is、id和uab;其中,Sp为预充电开关状态值、Rp为预充电电阻值、Rs为变压器二次侧漏电阻值、Ud为中间直流侧电压值、us为变压器二次侧输出交流电压值、Ls为变压器二次侧漏电感值,is为两电平单相四象限脉冲整流器交流侧的电流值,id为两电平单相四象限脉冲整流器直流侧的电流值;is为两电平单相四象限脉冲整流器交流侧的电流值;uab为脉冲整流器交流侧的电压值;(3)、将(2)所获得仿真值is、id和uab输出到牵引控制单元。
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处