发明名称 交流异步牵引电机定子绕组匝间短路的故障建模仿真方法
摘要 本发明涉及交流异步牵引电机定子绕组匝间短路的故障建模仿真方法,基于两相静止坐标系下交流异步牵引电机定子绕组匝间短路的数学模型完成,具体方法包括:计算交流异步牵引电机定子绕组匝间短路时参数,并基于此通过建立的交流异步牵引电机定子绕组匝间短路故障状态方程计算交流异步牵引电机状态变量,进而计算当前时刻异步牵引电机输出电流、电磁转矩以及转速值,从而完成交流异步牵引电机定子绕组匝间短路的故障建模仿真。本方法可实现离线仿真以及在线实时仿真系统中,交流异步牵引电机匝间短路故障的仿真。该方法中基于两相静止坐标系下交流异步牵引电机定子绕组匝间短路状态方程的建立,使得该方法易于计算机仿真实现。
申请公布号 CN103530442A 申请公布日期 2014.01.22
申请号 CN201310421982.7 申请日期 2013.09.16
申请人 西南交通大学 发明人 葛兴来;韩坤;冯晓云;熊成林;宋文胜;刘卉圻;苟斌
分类号 G06F17/50(2006.01)I 主分类号 G06F17/50(2006.01)I
代理机构 成都信博专利代理有限责任公司 51200 代理人 张澎
主权项 交流异步牵引电机定子绕组匝间短路的故障建模仿真方法,通过建立异步牵引电机定子绕组匝间短路故障模型进行异步牵引电机定子绕组匝间短路故障仿真,包括以下步骤:(1)确定并计算异步牵引电机定子短路绕组系数μ、参数σ、参数ζ、以及参数矩阵N: <mrow> <mi>&mu;</mi> <mo>=</mo> <msub> <mi>sa</mi> <mn>2</mn> </msub> <mo>/</mo> <mrow> <mo>(</mo> <msub> <mi>sa</mi> <mn>1</mn> </msub> <mo>+</mo> <msub> <mi>sa</mi> <mn>2</mn> </msub> <mo>)</mo> </mrow> </mrow> <mrow> <mi>&sigma;</mi> <mo>=</mo> <msub> <mi>L</mi> <mi>s</mi> </msub> <msub> <mi>L</mi> <mi>r</mi> </msub> <mo>-</mo> <msub> <mi>L</mi> <mi>m</mi> </msub> <msub> <mi>L</mi> <mi>m</mi> </msub> </mrow> <mrow> <mi>&zeta;</mi> <mo>=</mo> <msub> <mi>L</mi> <mi>ls</mi> </msub> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&mu;</mi> <msub> <mi>L</mi> <mi>m</mi> </msub> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&mu;</mi> <msub> <mi>L</mi> <mi>s</mi> </msub> </mrow> <mrow> <mi>N</mi> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mfrac> <msub> <mi>L</mi> <mi>r</mi> </msub> <mi>&sigma;</mi> </mfrac> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&mu;</mi> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msub> <mi>L</mi> <mi>r</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msub> <mi>L</mi> <mi>s</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msub> <mi>L</mi> <mi>s</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&mu;</mi> </mfrac> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> </mtr> </mtable> </mfenced> </mrow>其中,sa2为交流异步牵引电机单相定子绕组短路绕组匝数,sa1为交流异步牵引电机单相定子绕组正常绕组匝数,Ls为交流异步牵引电机定子电感参数,Lr为交流异步牵引电机转子电感参数,Lm为交流异步牵引电机定转子互感参数,Lls为交流异步牵引电机定子漏感参数;(2)建立基于两相静止坐标系的交流异步牵引电机定子绕组匝间短路的状态方程,并进行求解,计算异步牵引电机状态变量ψsα、ψsβ、ψrα、ψrβ、ψsα2: <mrow> <mi>p</mi> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>s&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>s&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>r&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>r&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mrow> <mi>sa</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mo>-</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <msub> <mi>N</mi> <mn>2</mn> </msub> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&mu;</mi> <msub> <mi>R</mi> <mi>s</mi> </msub> <msub> <mi>N</mi> <mn>5</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>R</mi> <mi>s</mi> </msub> <msub> <mi>N</mi> <mn>1</mn> </msub> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>R</mi> <mi>r</mi> </msub> <msub> <mi>N</mi> <mn>3</mn> </msub> <mo>-</mo> <msub> <mi>&omega;</mi> <mi>r</mi> </msub> <mo>[</mo> <mn>00010</mn> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <msub> <mi>R</mi> <mi>r</mi> </msub> <msub> <mi>N</mi> <mn>4</mn> </msub> <mo>+</mo> <msub> <mi>&omega;</mi> <mi>r</mi> </msub> <mo>[</mo> <mn>00100</mn> <mo>]</mo> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mi>&mu;</mi> <msub> <mi>R</mi> <mi>s</mi> </msub> <msub> <mi>N</mi> <mn>2</mn> </msub> <mo>+</mo> <mrow> <mo>(</mo> <mi>&mu;</mi> <msub> <mi>R</mi> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>R</mi> <mi>f</mi> </msub> <mo>)</mo> </mrow> <msub> <mi>N</mi> <mn>5</mn> </msub> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>s&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>s&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>r&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>r&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mrow> <mi>sa</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>+</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>u</mi> <mi>s&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>u</mi> <mi>s&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> </mtr> </mtable> </mfenced> </mrow>其中,p为微分算子,ψsα为交流异步牵引电机定子磁链在两相静止坐标系α轴上的分量,ψsβ为交流异步牵引电机定子磁链在两相静止坐标系β轴上的分量,ψrα为交流异步牵引电机转子磁链在两相静止坐标系α轴上的分量,ψrβ为交流异步牵引电机转子磁链在两相静止坐标系β轴上的分量,ψsα2为交流异步牵引电机定子短路绕组磁链在两相静止坐标系α轴上的分量,Nn为参数矩阵N中的第n行,Rs为交流异步牵引电机定子电阻参数,Rr为交流异步牵引电机转子电阻参数,Rf为交流异步牵引电机定子绕组短路电阻,usα为交流异步牵引电机输入电压在两相静止坐标系α轴上的分量,usβ为交流异步牵引电机输入电压在两相静止坐标系β轴上的分量;(3)根据异步牵引电机状态变量值计算当前时刻交流异步牵引电机输出电流isα、isβ、irα、irβ,交流异步牵引电机定子绕组短路电流if,交流异步牵引电机输出电磁转矩Te以及转速值ωr; <mrow> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>i</mi> <mi>s&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>s&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>r&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>r&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>i</mi> <mi>f</mi> </msub> </mtd> </mtr> </mtable> </mfenced> <mo>=</mo> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <mfrac> <msub> <mi>L</mi> <mi>r</mi> </msub> <mi>&sigma;</mi> </mfrac> <mo>+</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mi>&mu;</mi> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <mn>2</mn> <mn>3</mn> </mfrac> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msub> <mi>L</mi> <mi>r</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msub> <mi>L</mi> <mi>s</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <msub> <mi>L</mi> <mi>m</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mfrac> <msub> <mi>L</mi> <mi>s</mi> </msub> <mi>&sigma;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> </mtr> <mtr> <mtd> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mn>0</mn> </mtd> <mtd> <mo>-</mo> <mfrac> <mn>1</mn> <mi>&mu;</mi> </mfrac> <mfrac> <mn>1</mn> <mi>&zeta;</mi> </mfrac> </mtd> </mtr> </mtable> </mfenced> <mfenced open='[' close=']'> <mtable> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>s&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>s&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>r&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mi>r&beta;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <msub> <mi>&psi;</mi> <mrow> <mi>sa</mi> <mn>2</mn> </mrow> </msub> </mtd> </mtr> </mtable> </mfenced> </mrow> <mfenced open='' close=''> <mtable> <mtr> <mtd> <msub> <mi>T</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mn>3</mn> <mn>2</mn> </mfrac> <msub> <mi>n</mi> <mi>p</mi> </msub> <msub> <mi>L</mi> <mi>m</mi> </msub> <mrow> <mo>(</mo> <msub> <mi>i</mi> <mi>s&beta;</mi> </msub> <msub> <mi>i</mi> <mi>r&alpha;</mi> </msub> <mo>-</mo> <msub> <mi>i</mi> <mi>s&alpha;</mi> </msub> <msub> <mi>i</mi> <mi>r&beta;</mi> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>n</mi> <mi>p</mi> </msub> <mi>&mu;</mi> <msub> <mi>L</mi> <mi>m</mi> </msub> <msub> <mi>i</mi> <mi>f</mi> </msub> <msub> <mi>i</mi> <mi>r&alpha;</mi> </msub> </mtd> </mtr> <mtr> <mtd> <mfrac> <mrow> <mi>d</mi> <msub> <mi>&omega;</mi> <mi>r</mi> </msub> </mrow> <mi>dt</mi> </mfrac> <mo>=</mo> <mfrac> <msub> <mi>n</mi> <mi>p</mi> </msub> <mi>J</mi> </mfrac> <mrow> <mo>(</mo> <msub> <mi>T</mi> <mi>e</mi> </msub> <mo>-</mo> <msub> <mi>T</mi> <mi>L</mi> </msub> <mo>)</mo> </mrow> </mtd> </mtr> </mtable> </mfenced>其中,isα为交流异步牵引电机定子电流在两相静止坐标系α轴上的分量,isβ为交流异步牵引电机定子电流在两相静止坐标系β轴上的分量,irα为交流异步牵引电机转子电流在两相静止坐标系α轴上的分量,irβ为交流异步牵引电机转子电流在两相静止坐标系β轴上的分量,np为交流异步牵引电机极对数参数,J为交流异步牵引电机转动惯量参数,TL为交流异步牵引电机负载转矩参数;(4)将(3)所得计算数值输出至后续处理设备。
地址 610031 四川省成都市二环路北一段111号西南交通大学科技处