CN106292629B - 一种高铁牵引电机复合故障模拟方法 - Google Patents
一种高铁牵引电机复合故障模拟方法 Download PDFInfo
- Publication number
- CN106292629B CN106292629B CN201610691990.7A CN201610691990A CN106292629B CN 106292629 B CN106292629 B CN 106292629B CN 201610691990 A CN201610691990 A CN 201610691990A CN 106292629 B CN106292629 B CN 106292629B
- Authority
- CN
- China
- Prior art keywords
- traction motor
- speed rail
- rail traction
- phase
- stator
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 86
- 238000004804 winding Methods 0.000 claims abstract description 48
- 238000004088 simulation Methods 0.000 claims abstract description 9
- 230000003068 static effect Effects 0.000 claims description 94
- 230000004907 flux Effects 0.000 claims description 55
- 239000002131 composite material Substances 0.000 claims description 49
- 150000001875 compounds Chemical class 0.000 claims description 23
- 230000009466 transformation Effects 0.000 claims description 16
- 238000013178 mathematical model Methods 0.000 claims description 8
- 230000001131 transforming effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 9
- 238000011160 research Methods 0.000 description 7
- 238000005457 optimization Methods 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 229910052742 iron Inorganic materials 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000008646 thermal stress Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B23/00—Testing or monitoring of control systems or parts thereof
- G05B23/02—Electric testing or monitoring
- G05B23/0205—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults
- G05B23/0218—Electric testing or monitoring by means of a monitoring system capable of detecting and responding to faults characterised by the fault detection method dealing with either existing or incipient faults
-
- 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
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
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)
- Pure & Applied Mathematics (AREA)
- Operations Research (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明公开了一种高铁牵引电机复合故障模拟方法,通过建立高铁牵引电机定子匝间短路与转子断条复合故障模型进行高铁牵引电机定子匝间短路与转子断条复合故障的模拟,其中,建立高铁牵引电机复合故障模型的具体步骤为:根据高铁牵引电机转子断条与定子匝间短路复合故障发生原理,引入复合故障系数矩阵,建立高铁牵引复合故障三相坐标系下的模型;将三相坐标系下的高铁牵引复合故障模型通过坐标变换到两相静止坐标系下,进而得到简化的高铁牵引复合故障模型。本发明解决了高铁牵引多回路复合故障模型计算的复杂、维数大、参数获取困难、效率低的问题,只需要设置复合故障系数矩阵,就可以模拟高铁牵引定子匝间短路与转子断条复合故障。
Description
技术领域
本发明涉及一种高铁牵引电机复合故障模拟方法,属于复合故障技术领域。
背景技术
异步电机作为能量转换装置,在很多领域得到了广泛的应用。随着中国高速铁路向高速和重载方向发展,作为高速列车牵引系统的心脏的异步电机,受到了高度的重视。异步电机是高速列车牵引交流传动系统的核心部件之一,在高速列车前进时,异步电机将电能转换为机械能,在高速列车制动时,牵引电机作为发电机使用,将机械能转换为电能,回馈电网。但是,因为牵引电机的负载会随着需求的不同而频繁变化,工作的环境恶劣等,这些因素会增大其故障发生的概率。当列车发生故障时,若没有及时发现,会造成严重的伤亡事故,巨大的经济损失和社会影响。
由于高铁牵引电机转子导条与端环之间的连接处很薄弱,运行环境恶劣,以及转子负载波动,牵引电机转子需要承受很大的热应力与机械应力,所以电机转子导条很容易断裂。导条断裂产生的铜屑会损坏定子绕组的绝缘层,造成定子绕组匝间短路。所以在实际的工程中,牵引电机存在转子断条与定子匝间短路复合故障。因此研究转子断条与定子匝间短路复合故障建模方法是有必要的,以便于异步电机该复合故障诊断的研究。
国内外学者围绕异步电机单故障建模已经开展了很多工作,而对于复合故障建模的研究还很少。很多研究都是在三相坐标系子下建立异步电机的数学模型,通过增加故障相电阻阻值来模拟转子断条故障,减少故障相定子电阻阻值来模拟定子匝间短路故障。而对于异步电机定子匝间短路与转子断条的复合故障模型,相关研究很少,主要是多回路复合故障模型。但是多回路复合模型维数大,计算复杂,很多参数难以获得。
发明内容
针对现有多回路复合故障模型繁琐复杂的不足,本发明提供一种高铁牵引电机定子匝间短路与转子断条复合故障模拟方法,通过坐标变换矩阵,将三相坐标系下的异步电机复合故障模型转换到两相静止坐标系下,通过引入复合故障系数矩阵,注入复合故障,从而获得异步电机复合故障模型。该方法简单明了,物理意义明确,易于编程实现,计算量少,适合于实际应用,本发明解决了高铁牵引多回路复合故障模型计算的复杂、维数大、参数获取困难、效率低的问题,只需要设置复合故障系数矩阵,就可以模拟高铁牵引定子匝间短路与转子断条复合故障,以便于后续高铁牵引转子断条与定子匝间短复合故障的研究。
本发明为解决上述技术问题采用以下技术方案:
本发明提供一种高铁牵引电机复合故障模拟方法,通过建立高铁牵引电机定子匝间短路与转子断条复合故障模型进行高铁牵引电机定子匝间短路与转子断条复合故障的模拟。
其中,建立高铁牵引电机定子匝间短路与转子断条复合故障模型的具体步骤如下:
步骤1,根据高铁牵引电机定子匝间短路与转子断条复合故障的发生原理,计算复合故障系数矩阵;
步骤2,将复合故障系数矩阵引入三相坐标系下的高铁牵引电机数学模型,建立三相坐标系下的高铁牵引电机复合故障模型;
步骤3,采用坐标变换,将三相坐标系下的高铁牵引电机复合故障模型转化到两相静止坐标系下,从而得到两相静止坐标系下的高铁牵引电机复合故障模型。
作为本发明的进一步优化方案,步骤2中高铁牵引电机数学模型由电压方程、磁链方程、转矩方程和运动方程组成。
作为本发明的进一步优化方案,步骤1中根据高铁牵引电机定子匝间短路与转子断条复合故障的发生原理,计算复合故障系数矩阵,具体为:
和
其中,Fs为三相坐标系下定子匝间短路故障系数矩阵,Fr为三相坐标系下转子断条故障系数矩阵;Δs为定子匝间短路系数,表征定子匝间短路故障的严重程度,Δr为转子断条系数,表征转子断条的严重程度,n1为定子绕组匝间短路匝数,N1为定子绕组每相串联匝数,n2为转子导条断裂根数,N2为转子导条数。
作为本发明的进一步优化方案,步骤2中建立三相坐标系下的高铁牵引电机复合故障模型,具体为:
电压方程为:
其中,usABC=[uA uB uC]T为高铁牵引电机定子三相电压矩阵,uA为高铁牵引电机定子A相电压,uB为高铁牵引电机定子B相电压,uC为高铁牵引电机定子C相电压;urabc=[ua ubuc]T为高铁牵引电机转子三相电压矩阵,ua为高铁牵引电转子a相电压,ub为高铁牵引电机转子b相电压,uc为高铁牵引电机转子c相电压;Rs为高铁牵引电机定子三相电阻矩阵,R1表示高铁牵引电机定子每相电阻;Rr为高铁牵引电机转子三相电阻矩阵,R2表示高铁牵引电机转子每相电阻;isABC为高铁牵引电机定子三相电流矩阵,isABC=[iA iB iC]T,iA为高铁牵引电机定子A相电流,iB为高铁牵引电机定子B相电流,iC为高铁牵引电机定子C相电流;irabc为高铁牵引电机转子三相电流矩阵,irabc=[ia ib ic]T,ia为高铁牵引电机转子a相电流,ib为高铁牵引电机转子b相电流,ic为高铁牵引电机转子c相电流;p为微分算子;为高铁牵引电机定子三相磁链矩阵, 为高铁牵引电机定子A相磁链,为高铁牵引电机定子B相磁链,为高铁牵引电机定子C相磁链;为高铁牵引电机转子三相磁链矩阵, 为高铁牵引电机转子a相磁链,为高铁牵引电机转子b相磁链,为高铁牵引电机转子c相磁链;
磁链方程为:
其中,Lss为高铁牵引电机定子三相绕组的自感矩阵,Lm为高铁牵引电机定子绕组与转子绕组互感,Ll1为高铁牵引电机定子绕组漏感;Lrr为高铁牵引电机转子三相绕组的自感矩阵,Ll2为高铁牵引电机转子绕组漏感;Msr为高铁牵引电机定子绕组与转子绕组的互感矩阵,
θ为三相坐标系中高铁牵引电机定子A轴与转子a轴的夹角;
转矩方程为:
Te=-npLm[(iAia+iBib+iCic)sinθ+(iAib+iBic+iCia)sin(θ+2π/3)+(iAic+iBia+iCib)sin(θ-2π/3)]
运动方程为:
ω=pθ
其中,Te为电磁转矩,Tl表示负载转,J表示转动惯量,np表示极对数,ω为电机角速度。
作为本发明的进一步优化方案,步骤3中三相坐标变换到两相静止坐标的变换矩阵为:
其中,C3s/2s表示高铁牵引电机定子三相坐标系变换到两相静止αβ坐标系中的转换矩阵,C3r/2s表示高铁牵引电机转子三相坐标转换到两相静止αβ坐标系中的转换矩阵。
作为本发明的进一步优化方案,步骤3中采用坐标变换,将三相坐标系下的高铁牵引电机复合故障模型转化到两相静止坐标系下,具体为:
Te=-npLm[(iAia+iBib+iCic)sinθ+(iAib+iBic+iCia)sin(θ+2π/3)+(iAic+iBia+iCib)sin(θ-2π/3)]
其中,usαβ=[usα usβ]T为高铁牵引电机定子两相静止坐标系下电压矩阵,usα为高铁牵引电机定子两相静止坐标系下α轴电压,usβ为高铁牵引电机定子两相静止坐标系下β轴电压;urαβ=[urα urβ]T为高铁牵引电机转子两相静止坐标系下电压矩阵,urα为高铁牵引电机转子两相静止坐标系下α轴电压,urβ为高铁牵引电机转子两相静止坐标系下β轴电压;为高铁牵引电机定子两相静止坐标系下磁链矩阵,为高铁牵引电机定子两相静止坐标系下α轴磁链,为高铁牵引电机定子两相静止坐标系下β轴磁链;为高铁牵引电机转子两相静止坐标系下磁链矩阵,为高铁牵引电机转子两相静止坐标系下α轴磁链,为高铁牵引电机转子两相静止坐标系下β轴磁链;isαβ=[isαisβ]T为高铁牵引电机定子两相静止坐标系下电流矩阵,isα为高铁牵引电机定子两相静止坐标系下α轴电流,isβ为高铁牵引电机定子两相静止坐标系下β轴电流;irαβ=[irα irβ]T为高铁牵引电机转子两相静止坐标系下电流矩阵,irα为高铁牵引电机转子两相静止坐标系下α轴电流,irβ为高铁牵引电机转子两相静止坐标系下β轴电流;
两相静止坐标系下的高铁牵引电机复合故障模型具体为:
电压方程为:
其中,为高铁牵引电机两相静止坐标系下定子匝间短路故障系数矩阵,fsα为高铁牵引电机两相静止坐标系下定子α轴故障系数、fsβ为高铁牵引电机两相静止坐标系下定子β轴故障系数;为高铁牵引电机两相静止坐标系下转子断条故障系数矩阵,frα为高铁牵引电机两相静止坐标系下转子α轴故障系数、frβ为高铁牵引电机两相静止坐标系下转子β轴故障系数;
磁链方程为:
其中,为高铁牵引电机两相静止坐标系下定子绕组电感;为高铁牵引电机两相静止坐标系下转子绕组电感;为高铁牵引电机两相静止坐标系下定子绕组与转子绕组互感;
转矩方程为:
Te=npM(irαisβ-isαirβ)
运动方程为:
ω=pθ。
本发明采用以上技术方案与现有技术相比,具有以下技术效果:
(1)复合故障模型所需的参数容易获取;(2)降低了模型的维数,减少了计算量。(3)表达式清晰简洁,具有明确的物理意义,易于编程实现。
附图说明
图1是C3s/2s坐标转换矩阵框图。
图2是高铁牵引电机复合故障定子电压方程框图。
图3是高铁牵引电机复合故障转子电压方程框图。
图4是高铁牵引电机复合故障磁链方程框图。
图5是高铁牵引电机复合故障转矩方程框图。
图6是高铁牵引电机复合故障运动方程框图。
图7是高铁牵引电机复合故障模型整体框图。
图8是无故障高铁牵引电机转速图。
图9是高铁牵引电机复合故障转速图。
具体实施方式
下面结合附图对本发明的技术方案做进一步的详细说明:
本发明提供了一种高铁牵引电机定子匝间短路与转子断条复合故障模拟方法,模拟高铁牵引电机定子匝间短路与转子断条复合故障,得到复合故障数据,为复合诊断研究提供数据;解决高铁牵引电机定子匝间短路与转子断条复合故障如何模拟的技术问题;通过设置复合故障系数矩阵,可以模拟高铁牵引电机复合故障。所述定子匝间短路是指因定子绕组绝缘层损坏而造成短路的现象,所述转子断条是指转子导条断裂。
高铁牵引电机的数学模型主要由电压方程、磁链方程、转矩方程和运动方程部分组成。
(1)根据高铁牵引电机复合故障发生原理,计算复合故障系数矩阵。
其中,Δr为转子断条系数,表示转子断条的严重程度,Δs为定子匝间短路系数,表示定子匝间短路故障的严重程度,Fs为定子匝间短路故障系数矩阵,Fr为转子断条故障系数矩阵,N2为转子导条数,N1为定子绕组每相串联匝数,n2为转子导条断裂根数,n1为定子绕组匝间短路匝数。
(2)建立高铁牵引电机复合故障三相坐标系下的数学模型;
电压方程为:
其中,usABC=[uA uB uC]T为高铁牵引电机定子三相电压矩阵,uA为高铁牵引电机定子A相电压,uB为高铁牵引电机定子B相电压,uC为高铁牵引电机定子C相电压;urabc=[ua ubuc]T为高铁牵引电机转子三相电压矩阵,ua为高铁牵引电转子a相电压,ub为高铁牵引电机转子b相电压,uc为高铁牵引电机转子c相电压;Rs为高铁牵引电机定子三相电阻矩阵,R1表示高铁牵引电机定子每相电阻;Rr为高铁牵引电机转子三相电阻矩阵,R2表示高铁牵引电机转子每相电阻;isABC为高铁牵引电机定子三相电流矩阵,isABC=[iA iB iC]T,iA为高铁牵引电机定子A相电流,iB为高铁牵引电机定子B相电流,iC为高铁牵引电机定子C相电流;irabc为高铁牵引电机转子三相电流矩阵,irabc=[ia ib ic]T,ia为高铁牵引电机转子a相电流,ib为高铁牵引电机转子b相电流,ic为高铁牵引电机转子c相电流;p为微分算子;为高铁牵引电机定子三相磁链矩阵, 为高铁牵引电机定子A相磁链,为高铁牵引电机定子B相磁链,为高铁牵引电机定子C相磁链;为高铁牵引电机转子三相磁链矩阵, 为高铁牵引电机转子a相磁链,为高铁牵引电机转子b相磁链,为高铁牵引电机转子c相磁链;ur=0。
磁链方程为:
其中,Lss为高铁牵引电机定子三相绕组的自感矩阵,Lm为高铁牵引电机定子绕组与转子绕组互感,Ll1为高铁牵引电机定子绕组漏感;Lrr为高铁牵引电机转子三相绕组的自感矩阵,Ll2为高铁牵引电机转子绕组漏感;Msr为高铁牵引电机定子绕组与转子绕组的互感矩阵,
θ为三相坐标系中高铁牵引电机定子A轴与转子a轴的夹角。
转矩方程为:
Te=-npLm[(iAia+iBib+iCic)sinθ+(iAib+iBic+iCia)sin(θ+2π/3)+(iAic+iBia+iCib)sin(θ-2π/3)] (0.7)
运动方程为:
ω=pθ (0.9)
其中,Te为电磁转矩,Tl表示负载转,J表示转动惯量,np表示极对数,ω为电机角速度。
(3)如图1所示的坐标转换矩阵框图,采用坐标变换矩阵将三相坐标系下的复合故障数学模型转换到两相静止坐标系下。
a)坐标变换矩阵,C3s/2s如图1所示:
其中,C3s/2s表示高铁牵引电机定子三相坐标系变换到两相静止αβ坐标系中的转换矩阵,C3r/2s表示高铁牵引电机转子坐标转换到两相静止αβ坐标系中的转换矩阵。
b)两相静止坐标系下复合故障模型推导过程如下:
Te=-npLm[(iAia+iBib+iCic)sinθ+(iAib+iBic+iCia)sin(θ+2π/3)+(iAic+iBia+iCib)sin(θ-2π/3)] (0.15)
其中,usαβ=[usα usβ]T为高铁牵引电机定子两相静止坐标系下电压矩阵,usα为高铁牵引电机定子两相静止坐标系下α轴电压,usβ为高铁牵引电机定子两相静止坐标系下β轴电压;urαβ=[urα urβ]T为高铁牵引电机转子两相静止坐标系下电压矩阵,urα为高铁牵引电机转子两相静止坐标系下α轴电压,urβ为高铁牵引电机转子两相静止坐标系下β轴电压;为高铁牵引电机定子两相静止坐标系下磁链矩阵,为高铁牵引电机定子两相静止坐标系下α轴磁链,为高铁牵引电机定子两相静止坐标系下β轴磁链;为高铁牵引电机转子两相静止坐标系下磁链矩阵,为高铁牵引电机转子两相静止坐标系下α轴磁链,为高铁牵引电机转子两相静止坐标系下β轴磁链;isαβ=[isαisβ]T为高铁牵引电机定子两相静止坐标系下电流矩阵,isα为高铁牵引电机定子两相静止坐标系下α轴电流,isβ为高铁牵引电机定子两相静止坐标系下β轴电流;irαβ=[irα irβ]T为高铁牵引电机转子两相静止坐标系下电流矩阵,irα为高铁牵引电机转子两相静止坐标系下α轴电流,irβ为高铁牵引电机转子两相静止坐标系下β轴电流。
(4)最后得到两相静止坐标系下的高铁牵引电机定子匝间短路与转子断条复合故障模型。
如图2和图3所示,电压方程为:
其中,为高铁牵引电机两相静止坐标系下定子匝间短路故障系数矩阵,fsα为高铁牵引电机两相静止坐标系下定子α轴故障系数、fsβ为高铁牵引电机两相静止坐标系下定子β轴故障系数;为高铁牵引电机两相静止坐标系下转子断条故障系数矩阵,frα为高铁牵引电机两相静止坐标系下转子α轴故障系数、frβ为高铁牵引电机两相静止坐标系下转子β轴故障系数;urα=0,urβ=0。
如图4所示,磁链方程为:
其中,为高铁牵引电机两相静止坐标系下定子绕组电感;为高铁牵引电机两相静止坐标系下转子绕组电感;为高铁牵引电机两相静止坐标系下定子绕组与转子绕组互感。
如图5所示,转矩方程为:
Te=npM(irαisβ-isαirβ) (0.18)
如图6所示,运动方程为:
ω=pθ。 (1.29)
高铁牵引电机定子匝间短路与转子断条复合故障模型框图如图7所示,仿真结果如图8和图9所示。
以上所述,仅为本发明中的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可理解想到的变换或替换,都应涵盖在本发明的包含范围之内,因此,本发明的保护范围应该以权利要求书的保护范围为准。
Claims (6)
1.一种高铁牵引电机复合故障模拟方法,其特征在于,通过建立高铁牵引电机定子匝间短路与转子断条复合故障模型进行高铁牵引电机定子匝间短路与转子断条复合故障的模拟,其中,建立高铁牵引电机定子匝间短路与转子断条复合故障模型的具体步骤如下:
步骤1,根据高铁牵引电机定子匝间短路与转子断条复合故障的发生原理,计算复合故障系数矩阵;
步骤2,将复合故障系数矩阵引入三相坐标系下的高铁牵引电机数学模型,建立三相坐标系下的高铁牵引电机复合故障模型;
步骤3,采用坐标变换,将三相坐标系下的高铁牵引电机复合故障模型转化到两相静止坐标系下,从而得到两相静止坐标系下的高铁牵引电机复合故障模型。
2.根据权利要求1所述的一种高铁牵引电机复合故障模拟方法,其特征在于,步骤2中高铁牵引电机数学模型由电压方程、磁链方程、转矩方程和运动方程组成。
3.根据权利要求1所述的一种高铁牵引电机复合故障模拟方法,其特征在于,步骤1中根据高铁牵引电机定子匝间短路与转子断条复合故障的发生原理,计算复合故障系数矩阵,具体为:
和
其中,Fs为三相坐标系下定子匝间短路故障系数矩阵,Fr为三相坐标系下转子断条故障系数矩阵;Δs为定子匝间短路系数,表征定子匝间短路故障的严重程度,Δr为转子断条系数,表征转子断条的严重程度,n1为定子绕组匝间短路匝数,N1为定子绕组每相串联匝数,n2为转子导条断裂根数,N2为转子导条数。
4.根据权利要求3所述的一种高铁牵引电机复合故障模拟方法,其特征在于,步骤2中建立三相坐标系下的高铁牵引电机复合故障模型,具体为:
电压方程为:
其中,usABC=[uA uB uC]T为高铁牵引电机定子三相电压矩阵,uA为高铁牵引电机定子A相电压,uB为高铁牵引电机定子B相电压,uC为高铁牵引电机定子C相电压;urabc=[ua ub uc]T为高铁牵引电机转子三相电压矩阵,ua为高铁牵引电转子a相电压,ub为高铁牵引电机转子b相电压,uc为高铁牵引电机转子c相电压;Rs为高铁牵引电机定子三相电阻矩阵,R1表示高铁牵引电机定子每相电阻;Rr为高铁牵引电机转子三相电阻矩阵,R2表示高铁牵引电机转子每相电阻;isABC为高铁牵引电机定子三相电流矩阵,isABC=[iA iB iC]T,iA为高铁牵引电机定子A相电流,iB为高铁牵引电机定子B相电流,iC为高铁牵引电机定子C相电流;irabc为高铁牵引电机转子三相电流矩阵,irabc=[ia ib ic]T,ia为高铁牵引电机转子a相电流,ib为高铁牵引电机转子b相电流,ic为高铁牵引电机转子c相电流;p为微分算子;为高铁牵引电机定子三相磁链矩阵, 为高铁牵引电机定子A相磁链,为高铁牵引电机定子B相磁链,为高铁牵引电机定子C相磁链;为高铁牵引电机转子三相磁链矩阵, 为高铁牵引电机转子a相磁链,为高铁牵引电机转子b相磁链,为高铁牵引电机转子c相磁链;
磁链方程为:
其中,Lss为高铁牵引电机定子三相绕组的自感矩阵,Lm为高铁牵引电机定子绕组与转子绕组互感,Ll1为高铁牵引电机定子绕组漏感;Lrr为高铁牵引电机转子三相绕组的自感矩阵,Ll2为高铁牵引电机转子绕组漏感;Msr为高铁牵引电机定子绕组与转子绕组的互感矩阵,
θ为三相坐标系中高铁牵引电机定子A轴与转子a轴的夹角;
转矩方程为:
Te=-npLm[(iAia+iBib+iCic)sinθ+(iAib+iBic+iCia)sin(θ+2π/3)+(iAic+iBia+iCib)sin(θ-2π/3)]
运动方程为:
ω=pθ
其中,Te为电磁转矩,Tl表示负载转,J表示转动惯量,np表示极对数,ω为电机角速度。
5.根据权利要求4所述的一种高铁牵引电机复合故障模拟方法,其特征在于,步骤3中三相坐标变换到两相静止坐标的变换矩阵为:
其中,C3s/2s表示高铁牵引电机定子三相坐标系变换到两相静止αβ坐标系中的转换矩阵,C3r/2s表示高铁牵引电机转子三相坐标转换到两相静止αβ坐标系中的转换矩阵。
6.根据权利要求5所述的一种高铁牵引电机复合故障模拟方法,其特征在于,步骤3中采用坐标变换,将三相坐标系下的高铁牵引电机复合故障模型转化到两相静止坐标系下,具体为:
Te=-npLm[(iAia+iBib+iCic)sinθ+(iAib+iBic+iCia)sin(θ+2π/3)+(iAic+iBia+iCib)sin(θ-2π/3)]
其中,usαβ=[usα usβ]T为高铁牵引电机定子两相静止坐标系下电压矩阵,usα为高铁牵引电机定子两相静止坐标系下α轴电压,usβ为高铁牵引电机定子两相静止坐标系下β轴电压;urαβ=[urα urβ]T为高铁牵引电机转子两相静止坐标系下电压矩阵,urα为高铁牵引电机转子两相静止坐标系下α轴电压,urβ为高铁牵引电机转子两相静止坐标系下β轴电压;为高铁牵引电机定子两相静止坐标系下磁链矩阵,为高铁牵引电机定子两相静止坐标系下α轴磁链,为高铁牵引电机定子两相静止坐标系下β轴磁链;为高铁牵引电机转子两相静止坐标系下磁链矩阵,为高铁牵引电机转子两相静止坐标系下α轴磁链,为高铁牵引电机转子两相静止坐标系下β轴磁链;isαβ=[isαisβ]T为高铁牵引电机定子两相静止坐标系下电流矩阵,isα为高铁牵引电机定子两相静止坐标系下α轴电流,isβ为高铁牵引电机定子两相静止坐标系下β轴电流;irαβ=[irα irβ]T为高铁牵引电机转子两相静止坐标系下电流矩阵,irα为高铁牵引电机转子两相静止坐标系下α轴电流,irβ为高铁牵引电机转子两相静止坐标系下β轴电流;
两相静止坐标系下的高铁牵引电机复合故障模型具体为:
电压方程为:
其中,为高铁牵引电机两相静止坐标系下定子匝间短路故障系数矩阵,fsα为高铁牵引电机两相静止坐标系下定子α轴故障系数、fsβ为高铁牵引电机两相静止坐标系下定子β轴故障系数;为高铁牵引电机两相静止坐标系下转子断条故障系数矩阵,frα为高铁牵引电机两相静止坐标系下转子α轴故障系数、frβ为高铁牵引电机两相静止坐标系下转子β轴故障系数;
磁链方程为:
其中,为高铁牵引电机两相静止坐标系下定子绕组电感;为高铁牵引电机两相静止坐标系下转子绕组电感;为高铁牵引电机两相静止坐标系下定子绕组与转子绕组互感;
转矩方程为:
Te=npM(irαisβ-isαirβ)
运动方程为:
ω=pθ。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610691990.7A CN106292629B (zh) | 2016-08-18 | 2016-08-18 | 一种高铁牵引电机复合故障模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610691990.7A CN106292629B (zh) | 2016-08-18 | 2016-08-18 | 一种高铁牵引电机复合故障模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106292629A CN106292629A (zh) | 2017-01-04 |
CN106292629B true CN106292629B (zh) | 2018-08-24 |
Family
ID=57661584
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610691990.7A Expired - Fee Related CN106292629B (zh) | 2016-08-18 | 2016-08-18 | 一种高铁牵引电机复合故障模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106292629B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108279362B (zh) * | 2017-12-15 | 2020-04-10 | 上海电力学院 | 海上双馈电机转子绕组匝间短路故障辨识方法 |
CN109061475B (zh) * | 2018-11-01 | 2021-09-10 | 上海科梁信息科技股份有限公司 | 永磁同步电机匝间短路故障模拟方法、装置、系统、介质 |
CN113378352A (zh) * | 2021-05-10 | 2021-09-10 | 电子科技大学 | 一种异步电机的多故障建模方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201000476Y (zh) * | 2007-01-19 | 2008-01-02 | 中广核工程有限公司 | 一种发电机转子绕组匝间短路在线诊断仪 |
EP2541217A1 (en) * | 2011-06-29 | 2013-01-02 | ABB Research Ltd. | A method for identifying a fault in an electrical machine |
CN103530442A (zh) * | 2013-09-16 | 2014-01-22 | 西南交通大学 | 交流异步牵引电机定子绕组匝间短路的故障建模仿真方法 |
CN103955556A (zh) * | 2014-03-27 | 2014-07-30 | 北京交通大学 | 高速铁路列车运行控制车载系统故障逻辑建模方法 |
CN104615121A (zh) * | 2014-12-04 | 2015-05-13 | 深圳市永达电子股份有限公司 | 一种列车故障诊断方法及系统 |
CN105628406A (zh) * | 2015-12-28 | 2016-06-01 | 中南大学 | 高速列车牵引传动控制系统故障注入方法及仿真系统 |
-
2016
- 2016-08-18 CN CN201610691990.7A patent/CN106292629B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN201000476Y (zh) * | 2007-01-19 | 2008-01-02 | 中广核工程有限公司 | 一种发电机转子绕组匝间短路在线诊断仪 |
EP2541217A1 (en) * | 2011-06-29 | 2013-01-02 | ABB Research Ltd. | A method for identifying a fault in an electrical machine |
CN103530442A (zh) * | 2013-09-16 | 2014-01-22 | 西南交通大学 | 交流异步牵引电机定子绕组匝间短路的故障建模仿真方法 |
CN103955556A (zh) * | 2014-03-27 | 2014-07-30 | 北京交通大学 | 高速铁路列车运行控制车载系统故障逻辑建模方法 |
CN104615121A (zh) * | 2014-12-04 | 2015-05-13 | 深圳市永达电子股份有限公司 | 一种列车故障诊断方法及系统 |
CN105628406A (zh) * | 2015-12-28 | 2016-06-01 | 中南大学 | 高速列车牵引传动控制系统故障注入方法及仿真系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106292629A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106292629B (zh) | 一种高铁牵引电机复合故障模拟方法 | |
WO2019007354A1 (zh) | 一种考虑不同转差的双馈风力发电系统短路电流计算方法 | |
Hao et al. | Theoretical analysis and calculation model of the electromagnetic torque of nonsalient-pole synchronous machines with interturn short circuit in field windings | |
CN103346728B (zh) | 双馈风力发电机定转子电流传感器故障判断方法 | |
CN103427751A (zh) | 永磁同步电机静态参数在线辨识的装置与方法 | |
CN109061475B (zh) | 永磁同步电机匝间短路故障模拟方法、装置、系统、介质 | |
CN103530442A (zh) | 交流异步牵引电机定子绕组匝间短路的故障建模仿真方法 | |
Jiang et al. | Design and verification of high speed permanent magnet synchronous motor for electric car | |
CN105938515A (zh) | 一种双馈电机三相短路电流实用计算方法 | |
Kallio | Modeling and parameter estimation of double-star permanent magnet synchronous machines | |
Aliyev et al. | Synthesis of the Exact Parameters of the Electromagnetic Brake of a Wind Electric Installation. | |
CN108599660A (zh) | 永磁同步电机定子绕组不对称故障的矢量控制方法 | |
Zhang et al. | Equivalent circuit of DFIG with crowbar protection for relay protection analysis | |
CN102570952B (zh) | 双馈式风力发电系统用于低电压穿越的一种去磁控制方法 | |
CN102081691A (zh) | 考虑频率变化的电动机实用化建模方法 | |
CN105720874B (zh) | 基于分布参数的电机气隙磁场建模方法及其应用 | |
Jung et al. | Stator winding fault influence on the field coil of a 10 MW superconducting synchronous generator | |
Wei et al. | Short circuit current analysis of DFIG with crowbar under unsymmetrical grid fault | |
CN105550496B (zh) | 一种双馈风力发电机定子绕组匝间短路数学建模方法 | |
Xiao et al. | Stator winding inter-turn short circuit and rotor eccentricity diagnosis of permanent magnet synchronous generator | |
CN110346721B (zh) | 双支路交流永磁电机损耗的测试方法 | |
Dorji et al. | DQ mathematical modelling and simulation of three-phase induction motor for electrical fault analysis | |
CN101340121A (zh) | 一种发电机内部故障暂态仿真简化模型及其建模方法和应用 | |
El Menzhi et al. | Induction motor fault diagnosis using voltage spectrum of auxiliary winding and lissajous curve of its park components | |
Ligang et al. | Optimal speed tracking for double fed wind generator via switching control |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180824 |