CN112083349A - 一种永磁同步电机定子绕组匝间短路故障诊断方法 - Google Patents
一种永磁同步电机定子绕组匝间短路故障诊断方法 Download PDFInfo
- Publication number
- CN112083349A CN112083349A CN202010765090.9A CN202010765090A CN112083349A CN 112083349 A CN112083349 A CN 112083349A CN 202010765090 A CN202010765090 A CN 202010765090A CN 112083349 A CN112083349 A CN 112083349A
- Authority
- CN
- China
- Prior art keywords
- matrix
- turn
- variable matrix
- short circuit
- state variable
- 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.)
- Granted
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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/52—Testing for short-circuits, leakage current or ground faults
-
- 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/50—Testing of electric apparatus, lines, cables or components for short-circuits, continuity, leakage current or incorrect line connections
- G01R31/72—Testing of electric windings
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Electric Motors In General (AREA)
- Control Of Motors That Do Not Use Commutators (AREA)
Abstract
本发明提出的一种永磁同步电机定子绕组匝间短路故障诊断方法,首先,从传感器中获取k‑1时刻的三相定子电流,结合定子电阻模型构建当前状态变量矩阵X;再从编码器中得到永磁同步电机的电角度,计算电角速度和三相电压,构建输入变量矩阵u;然后,输入所设置的初始估计矩阵,根据k‑1时刻的状态量和误差协方差进行预估计得到更新后的k时刻的误差协方差和状态变量矩阵;最后,通过状态向量定义故障指标,当故障指标最大变化量超过阈值时,匝间短路故障诊定。本发明诊断精确度高,不局限于两相模型,并把纳入计算的误差独立存在,不受测量数据的影响,同时便于实时处理,还兼顾了断相检测。
Description
技术领域
本发明涉及一种永磁同步电机定子绕组匝间短路故障诊断方法,可应用于永磁同步电机定子绕组匝间短路故障诊断应用场合。
背景技术
电机是能够实现机电能量之间转换的设备,由固定的定子和转动的转子两个基本部分组成。在众多种类中,永磁同步电机以其自身性能好、体积小、结构简单、调速范围广等优点,已经被人们广泛的应用于电动汽车领域。在长期能量转换过程中,电机受供电情况、负载工况和运行环境等影响,定子绕组极易发生故障,其中匝间短路故障占有很大的比例。定子绕组匝间短路会导致电机中一相绕组阻抗变化,从而使三相电流失衡,电动机转矩降低,以至出现明显的振动现象。同时随着故障线圈温度升高,会影响更多线圈发生故障,从而导致永磁同步电机燃烧等严重后果。
目前,大多方法都是基于故障特征向量的判断方法,常采用小波分解、神经算法等,大量的计算不利于在线诊断,同时大多数诊断策略进行故障检测时都依赖于永磁同步电动机的两相模型(αβ或dq坐标系),但这些模型还不足以适用于定子绕组不平衡的所有情况。因此,提出一种基于三相模型可实现的诊断方法。首先,从传感器中获取k-1时刻的三相定子电流,结合定子电阻模型构建当前状态变量矩阵X;再通过从编码器中得到永磁同步电机的电角度,计算电角速度和三相电压,构建输入变量矩阵u;然后,输入所设置的初始估计矩阵,根据k-1时刻的状态量和误差协方差进行预估计得到更新后的k时刻的误差协方差和状态变量矩阵;最后,通过状态向量定义故障指标,当故障指标最大变化量超过阈值时,匝间短路故障诊定。该方法诊断精确度高,不局限于两相模型,并把纳入计算的误差独立存在,不受测量数据的影响,同时便于实时处理,还兼顾了断相检测。
发明内容
技术问题:针对上述说明,提出的一种永磁同步电机定子绕组匝间短路故障诊断方法,能够确保较高的精确度不局限于两相模型,并把纳入计算的误差独立存在,不受测量数据的影响,同时便于实时处理。
技术方案:一种永磁同步电机定子绕组匝间短路故障诊断方法,包括如下步骤:
步骤一、获取三相定子电流is(k-1)(s=a,b,c),构建状态变量矩阵;
从传感器中获取k-1时刻的三相定子电流is(k-1)(s=a,b,c),并结合定子电阻模型Rabc构建当前状态变量矩阵X;
步骤二、构建输入变量矩阵u;
联合从编码器中得到的永磁同步电机的电角度θ,计算电角速度ω和三相电压uabc,构建输入变量矩阵u;
步骤三、更新误差协方差和状态变量矩阵;
给定预估初始状态变量矩阵和初始误差协方差P0的初始估计矩阵,对k时刻的状态变量矩阵和误差协方差Pk -进行预估计并进行卡尔曼增益Kk计算,从而,根据测量残差和卡尔曼增益Kk分别更新误差协方差和状态变量矩阵;
步骤四、通过故障指标计算及比较,诊定出匝间短路故障;
通过状态变量矩阵中的电阻值计算故障指标,当故障指标最大变化量超过阈值时,匝间短路故障诊定。
进一步的,所述步骤一中的状态变量矩阵的搭建具体包括:将通过传感器测量所得的三相定子电流is(k-1)(s=a,b,c),与式(1)所示的定子电阻模型Rabc联合搭建出如式(2)所示的状态变量矩阵X;
X=[ia ib ic Ra Rb Rc]T (2)。
进一步的,所述步骤二中输入变量矩阵u的获取方法为:将编码器中得到的永磁同步电机的电角度θ通过式(3)计算电角速度ω,将三相电流is(k-1)(s=a,b,c)得到的测量变量矩阵z通过公式(5)计算相电压uabc,从而获得如式(6)所示的输入变量矩阵u;
z=[ia ib ic]T (4)
u=[ua ub uc ω θ]T (6)。
进一步的,所述步骤三中的误差协方差和状态变量矩阵更新方法为:给定初始估计状态变量矩阵和误差协方差P0=diag[0.5 0.5 0.5 0.5 0.50.5],通过式(7)、(8)对k时刻的状态变量矩阵和误差协方差Pk -进行预估计;再将得到的预估值Pk -带入式(9)进行卡尔曼增益Kk计算,从而,根据测量残差和卡尔曼增益Kk经过式(10)、式(11)分别更新得到最终的误差协方差和状态变量矩阵加了负号上标,表示忽略过程和测量噪声,根据当前先验状态得到的无噪声估计值;
式中,表示状态量的初始估计,uk表示k时刻的输入变量矩阵;令 A为上一状态到当前状态的转换矩阵,B为输入量到当前状态的转换矩阵,Pk-1表示误差协方差初始估计,F表示输出线性化,可定义为单位矩阵;Q表示过程噪声矩阵,可设置为Q=diag[10-4 10-4 10-4 10-10 10-10 10-10];
式中,I为单位矩阵,zk为k时刻的测量变量矩阵。
进一步的,所述步骤4故障指标计算与比较方法为:将经过卡尔曼观测器后得到的状态变量矩阵中的电阻阵通过式(12)、(13)计算得到故障指标δmm;设置阈值δTH,如果δmm超过阈值δTH,则匝间短路故障诊定;
δmm=max(δ)-min(δ) (13)。
有益效果:本发明是基于卡尔曼观测器的诊断方法,以k-1时刻的最优估计为准,预测k时刻的状态变量,同时又对该状态进行观测,得到观测变量z,再在预测与观测之间进行分析,或者说是以观测量对预测量进行修正,从而得到k时刻的最优状态估计因此该诊断方法能够确保较高的精确度,并把纳入计算的误差独立存在,不受测量数据的影响,同时便于实时处理。
附图说明
图1为永磁同步电机定子绕组匝间短路故障图;
图2为一种永磁同步电机定子绕组匝间短路故障诊断方法控制框图;
图3为永磁同步电机定子绕组a相匝间短路故障时三相电流波形;
图4为匝间短路故障时故障指标响应曲线。
具体实施方式
下面结合附图并通过实施例对本发明作进一步的详细说明,以下实施例是对本发明的解释而本发明并不局限于以下实施例。
永磁同步电机定子绕组匝间短路故障图如图1所示,假设a相发生匝间短路故障,则短路匝被接触电阻Rf短路,增加一个短路回路。
如图2所示,永磁同步电机定子绕组匝间短路故障诊断方法包括如下步骤:
步骤1:获取三相定子电流is(k-1)(s=a,b,c),构建状态变量矩阵:
将通过传感器测量所得的三相定子电流is(k-1)(s=a,b,c),与式(1)所示的定子电阻模型Rabc联合搭建出如式(2)所示的状态变量矩阵X。
X=[ia ib ic Ra Rb Rc]T (2)
步骤2:构建输入变量矩阵u;
将编码器中得到的永磁同步电机的电角度θ通过式(3)计算电角速度ω,将三相电流is(k-1)(s=a,b,c)得到的测量变量矩阵z通过公式(5)计算相电压uabc,从而获得如式(6)所示的输入变量矩阵u:
z=[ia ib ic]T (4)
u=[ua ub uc ω θ]T (6)
步骤3:更新误差协方差和状态变量矩阵;
所述误差协方差和状态变量矩阵的更新方法:具体如下:
给定初始估计状态变量矩阵和误差协方差P0=diag[0.5 0.5 0.5 0.5 0.5 0.5],通过式(7)、(8)对k时刻的状态变量矩阵和误差协方差Pk -进行预估计。再将得到的预估值Pk -带入式(9)进行卡尔曼增益Kk计算,从而,根据测量残差和卡尔曼增益Kk经过式(10)、式(11)分别更新得到最终的误差协方差和状态变量矩阵(加了负号上标,表示忽略过程和测量噪声,根据当前先验状态得到的无噪声估计值)
式中,表示状态量的初始估计,uk表示k时刻的输入变量矩阵;令 A为上一状态到当前状态的转换矩阵,B为输入量到当前状态的转换矩阵,Pk-1表示误差协方差初始估计,F表示输出线性化,可定义为单位矩阵;Q表示过程噪声矩阵,可设置为Q=diag[10-4 10-4 10-4 10-10 10-10 10-10]。
式中,I为单位矩阵,zk为k时刻的测量变量矩阵。
步骤4:通过故障指标计算及比较,诊定出匝间短路故障;
所述具体的故障指标计算与比较方法,如下:
将经过卡尔曼观测器后得到的状态变量矩阵中的电阻阵通过式(12)、(13)计算得到故障指标δmm。设置阈值δTH,如果δmm超过阈值δTH,则匝间短路故障诊定。
δmm=max(δ)-min(δ) (13)
本发明方法首先获取k-1时刻的三相定子电流is(k-1)(s=a,b,c),结合定子电阻模型Rabc构建当前状态变量矩阵X;再利用三相定子电流构为测量变量矩阵z,联合电角度θ,计算电角速度ω和三相电压uabc,构建输入变量矩阵u;然后,给定状态变量初始估计矩阵和初始误差协方差P0,对k时刻的状态变量矩阵和误差协方差Pk -进行预估计并进行卡尔曼增益Kk计算,从而,根据测量残差和卡尔曼增益Kk分别更新误差协方差和状态变量矩阵;最后,通过状态变量矩阵中的电阻值计算故障指标,当故障指标最大变化量超过阈值时,匝间短路故障诊定。
永磁同步电机定子绕组a相匝间短路故障时三相电流波形如图3所示,a相电流增加的幅度远大于b、c两相。图4为匝间短路故障时故障指标曲线,可以看出对于δTH在设置合理范围内,永磁同步电机的匝间短路故障总能被诊断,例如根据此时响应,δTH设置为0.005。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (5)
1.一种永磁同步电机定子绕组匝间短路故障诊断方法,其特征在于,包括如下步骤:
步骤一、获取三相定子电流is(k-1)(s=a,b,c),构建状态变量矩阵;
从传感器中获取k-1时刻的三相定子电流is(k-1)(s=a,b,c),并结合定子电阻模型Rabc构建当前状态变量矩阵X;
步骤二、构建输入变量矩阵u;
联合从编码器中得到的永磁同步电机的电角度θ,计算电角速度ω和三相电压uabc,构建输入变量矩阵u;
步骤三、更新误差协方差和状态变量矩阵;
给定预估初始状态变量矩阵和初始误差协方差P0的初始估计矩阵,对k时刻的状态变量矩阵和误差协方差Pk -进行预估计并进行卡尔曼增益Kk计算,从而,根据测量残差和卡尔曼增益Kk分别更新误差协方差和状态变量矩阵;
步骤四、通过故障指标计算及比较,诊定出匝间短路故障;
通过状态变量矩阵中的电阻值计算故障指标,当故障指标最大变化量超过阈值时,匝间短路故障诊定。
4.根据权利要求1所述的永磁同步电机定子绕组匝间短路故障诊断方法,其特征在于:所述步骤三中的误差协方差和状态变量矩阵更新方法为:给定初始估计状态变量矩阵和误差协方差P0=diag[0.5 0.5 0.5 0.5 0.5 0.5],通过式(7)、(8)对k时刻的状态变量矩阵和误差协方差Pk -进行预估计;再将得到的预估值Pk -带入式(9)进行卡尔曼增益Kk计算,从而,根据测量残差和卡尔曼增益Kk经过式(10)、式(11)分别更新得到最终的误差协方差和状态变量矩阵加了负号上标,表示忽略过程和测量噪声,根据当前先验状态得到的无噪声估计值;
式中,表示状态量的初始估计,uk表示k时刻的输入变量矩阵;令 A为上一状态到当前状态的转换矩阵,B为输入量到当前状态的转换矩阵,Pk-1表示误差协方差初始估计,F表示输出线性化,可定义为单位矩阵;Q表示过程噪声矩阵,可设置为Q=diag[10-4 10-4 10-410-10 10-10 10-10];
式中,I为单位矩阵,zk为k时刻的测量变量矩阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010765090.9A CN112083349B (zh) | 2020-08-01 | 2020-08-01 | 一种永磁同步电机定子绕组匝间短路故障诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010765090.9A CN112083349B (zh) | 2020-08-01 | 2020-08-01 | 一种永磁同步电机定子绕组匝间短路故障诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112083349A true CN112083349A (zh) | 2020-12-15 |
CN112083349B CN112083349B (zh) | 2023-04-07 |
Family
ID=73735941
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010765090.9A Active CN112083349B (zh) | 2020-08-01 | 2020-08-01 | 一种永磁同步电机定子绕组匝间短路故障诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112083349B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112688608A (zh) * | 2020-12-25 | 2021-04-20 | 北京航空航天大学 | 一种三相永磁同步电机控制系统的故障诊断方法 |
CN113064073A (zh) * | 2021-03-12 | 2021-07-02 | 合肥恒大江海泵业股份有限公司 | 基于剩余电流的永磁同步电机匝间短路故障诊断方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100194324A1 (en) * | 2009-02-03 | 2010-08-05 | General Electric Company | Stator turn fault detection apparatus and method for induction machine |
CN104777426A (zh) * | 2015-04-17 | 2015-07-15 | 河海大学 | 一种基于无迹变换强跟踪的发电机动态状态估计方法 |
CN106019045A (zh) * | 2016-05-16 | 2016-10-12 | 安徽大学 | 一种永磁同步电机匝间短路故障诊断的方法 |
CN107179442A (zh) * | 2017-04-06 | 2017-09-19 | 安徽大学 | 一种永磁同步电机定子绕组电阻不对称故障的诊断方法 |
CN107425548A (zh) * | 2017-09-11 | 2017-12-01 | 河海大学 | 一种插值h∞扩展卡尔曼滤波发电机动态状态估计方法 |
CN108134549A (zh) * | 2017-12-25 | 2018-06-08 | 西安理工大学 | 一种提高永磁同步电机转速估计稳定性的方法 |
US20180294751A1 (en) * | 2017-04-06 | 2018-10-11 | Korea Electronics Technology Institute | Apparatus, system and method of fault diagnosis for permanent magnet motor |
US10355634B1 (en) * | 2018-06-12 | 2019-07-16 | GM Global Technology Operations LLC | Isolation of particular fault conditions in an electric machine assembly |
CN110208642A (zh) * | 2019-05-22 | 2019-09-06 | 中南大学 | 永磁同步电机定子匝间短路故障退化过程模拟方法及系统 |
CN111123104A (zh) * | 2020-02-14 | 2020-05-08 | 东南大学 | 一种无需先验知识的永磁电机绕组故障诊断方法 |
-
2020
- 2020-08-01 CN CN202010765090.9A patent/CN112083349B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20100194324A1 (en) * | 2009-02-03 | 2010-08-05 | General Electric Company | Stator turn fault detection apparatus and method for induction machine |
CN104777426A (zh) * | 2015-04-17 | 2015-07-15 | 河海大学 | 一种基于无迹变换强跟踪的发电机动态状态估计方法 |
CN106019045A (zh) * | 2016-05-16 | 2016-10-12 | 安徽大学 | 一种永磁同步电机匝间短路故障诊断的方法 |
CN107179442A (zh) * | 2017-04-06 | 2017-09-19 | 安徽大学 | 一种永磁同步电机定子绕组电阻不对称故障的诊断方法 |
US20180294751A1 (en) * | 2017-04-06 | 2018-10-11 | Korea Electronics Technology Institute | Apparatus, system and method of fault diagnosis for permanent magnet motor |
CN107425548A (zh) * | 2017-09-11 | 2017-12-01 | 河海大学 | 一种插值h∞扩展卡尔曼滤波发电机动态状态估计方法 |
CN108134549A (zh) * | 2017-12-25 | 2018-06-08 | 西安理工大学 | 一种提高永磁同步电机转速估计稳定性的方法 |
US10355634B1 (en) * | 2018-06-12 | 2019-07-16 | GM Global Technology Operations LLC | Isolation of particular fault conditions in an electric machine assembly |
CN110208642A (zh) * | 2019-05-22 | 2019-09-06 | 中南大学 | 永磁同步电机定子匝间短路故障退化过程模拟方法及系统 |
CN111123104A (zh) * | 2020-02-14 | 2020-05-08 | 东南大学 | 一种无需先验知识的永磁电机绕组故障诊断方法 |
Non-Patent Citations (3)
Title |
---|
BRICE AUBERT 等: "Stator winding fault diagnosis in permanent magnet synchronous generators based on short‐circuited turns identification using extended Kalman filter", 《ACTA IMEKO》 * |
刘占伟: "一种参数模型的永磁同步电机故障诊断控制仿真", 《微电机》 * |
王延峰等: "永磁同步电机定子匝间短路故障阻抗参数分析", 《电机与控制应用》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112688608A (zh) * | 2020-12-25 | 2021-04-20 | 北京航空航天大学 | 一种三相永磁同步电机控制系统的故障诊断方法 |
CN112688608B (zh) * | 2020-12-25 | 2022-08-02 | 北京航空航天大学 | 一种三相永磁同步电机控制系统的故障诊断方法 |
CN113064073A (zh) * | 2021-03-12 | 2021-07-02 | 合肥恒大江海泵业股份有限公司 | 基于剩余电流的永磁同步电机匝间短路故障诊断方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112083349B (zh) | 2023-04-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110441643B (zh) | 永磁同步电机控制系统中逆变器功率管断路故障诊断方法 | |
CN111193448B (zh) | 基于扩展卡尔曼滤波器的永磁同步电机负载转矩观测方法 | |
CN112083349B (zh) | 一种永磁同步电机定子绕组匝间短路故障诊断方法 | |
CN109713971B (zh) | 一种永磁同步电机的扰动抑制方法 | |
CN113241985B (zh) | 无位置传感器磁悬浮飞轮电流自校正控制装置及方法 | |
Luukko et al. | Estimation of the flux linkage in a direct-torque-controlled drive | |
CN113676088B (zh) | 带有谐波抑制的永磁同步电机无速度传感器控制方法 | |
CN108521246A (zh) | 永磁同步电机单电流传感器电流预测控制的方法及装置 | |
El Sayed et al. | Comparative analysis between unscented and extended Kalman filters for PMSG inter-turn fault detection | |
CN114024382B (zh) | 一种基于负载转矩反馈的永磁同步电机单转速环调速装置 | |
Salvatore et al. | Improved rotor speed estimation using two Kalman filter-based algorithms | |
Uma et al. | Conventional controller design for switched reluctance motor drives | |
CN108649849A (zh) | 一种简单的无传感器永磁同步电机速度估测方法 | |
CN109962647B (zh) | 带有双模耦合结构观测器的电机控制系统及状态估计方法 | |
Meng et al. | A non-invasive dual-EKF-based rotor temperature estimation technique for permanent magnet machines | |
JP4281054B2 (ja) | 誘導電動機の速度センサレスベクトル制御方法及び装置 | |
Liu et al. | Sensorless Fault-Tolerant Control of A Nine-phase Permanent Magnet Synchronous Motor Under One-phase Open-Circuited Fault | |
CN115001342B (zh) | 感应电机的转速估算方法、装置及系统 | |
Kadrine et al. | High Gain Observer with Updated Gain for Sensorless Induction Motor | |
CN114465542B (zh) | 异步电机速度传感器故障容错控制方法 | |
CN112234888B (zh) | 一种用于电机的电流辨识方法、电机控制器及车辆 | |
CN113777485B (zh) | 基于价值函数误差的永磁同步电机开路故障诊断方法 | |
CN110460279B (zh) | 一种基于辅助变量的扩展转速自适应观测器低速控制方法 | |
CN117478015A (zh) | 一种固定型号电动执行机构转矩间接测量方法及系统 | |
CN117013918A (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 |