《一种笼型异步电动机转子断条故障检测方法.pdf》由会员分享,可在线阅读,更多相关《一种笼型异步电动机转子断条故障检测方法.pdf(20页完整版)》请在专利查询网上搜索。
1、(10)申请公布号 CN 102944842 A (43)申请公布日 2013.02.27 CN 102944842 A *CN102944842A* (21)申请号 201210502826.9 (22)申请日 2012.11.30 G01R 31/34(2006.01) G01R 23/16(2006.01) (71)申请人 华北电力大学 (保定) 地址 071003 河北省保定市永华北大街 619 号 (72)发明人 许伯强 孙丽玲 (74)专利代理机构 石家庄冀科专利商标事务所 有限公司 13108 代理人 李羡民 高锡明 (54) 发明名称 一种笼型异步电动机转子断条故障检测方法 (。
2、57) 摘要 一种笼型异步电动机转子断条故障检测方 法, 它首先对按一定频率采集的定子电流瞬时信 号应用希尔伯特变换, 得到其希尔伯特模量 ; 然 后通过减去其平均值而滤除该希尔伯特模量中的 直流分量 ; 再应用 ESPRIT 进行频谱分析, 获得其 ESPRIT频谱图 ; 最后根据ESPRIT频谱图中是否存 在 2sf1分量、 4sf1分量谱峰判断是否存在转子断 条故障 : 若同时存在 2sf1分量和 4sf1分量谱峰, 则转子断条, 否则转子正常。 本发明有效克服了负 荷波动的影响, 仅需采样时间很短的定子电流信 号即可高灵敏度、 高可靠性地检测异步电动机转 子断条故障, 特别适用于工程中。
3、广泛采用的大型 异步电动机。 (51)Int.Cl. 权利要求书 3 页 说明书 11 页 附图 5 页 (19)中华人民共和国国家知识产权局 (12)发明专利申请 权利要求书 3 页 说明书 11 页 附图 5 页 1/3 页 2 1. 一种笼型异步电动机转子断条故障检测方法, 其特征是, 它首先对按一定频率采集 的定子电流瞬时信号应用希尔伯特变换, 得到其希尔伯特模量 ; 然后通过减去其平均值而 滤除该希尔伯特模量中的直流分量 ; 再应用 ESPRIT 进行频谱分析, 获得其 ESPRIT 频谱图 ; 最后根据 ESPRIT 频谱图中是否存在 2sf1分量、 4sf1分量谱峰判断是否存在转。
4、子断条故障 : 若同时存在 2sf1分量和 4sf1分量谱峰, 则转子断条, 否则转子正常。 2. 根据权利要求 1 所述笼型异步电动机转子断条故障检测方法, 其特征是, 具体按以 下步骤进行 : a. 测取一相定子电流瞬时信号: 对于高压电机, 采用一只电流钳在电流互感器 CT 二次侧测取一相定子电流瞬时信号 ; 对于低压电机, 采用一只电流钳直接在电机接线端子处测取一相定子电流瞬时信号 ; b. 采用单工频周期滑动窗方法计算定子电流瞬时信号的有效值, 通过分析有效 值的变化趋势, 提取其最平稳亦即波动最小的一段数据, 记为; c. 对提取的定子电流瞬时信号应用希尔伯特变换, 得到其希尔伯特。
5、模量 Mod ; d. 滤除该希尔伯特模量中的直流分量, 获得待分析信号 M, M=Mod-mean(Mod), mean(Mod) 表示 Mod 的平均值 ; e. 应用 ESPRIT 对待分析信号 M 进行频谱分析, 获得 ESPRIT 频谱图 ; f. 根据 ESPRIT 频谱图中是否存在 2sf1分量、 4sf1分量谱峰判断是否存在转子断条故 障 : 若同时存在 2sf1分量、 4sf1分量谱峰, 则转子断条, 否则转子正常。 3. 根据权利要求 2 所述笼型异步电动机转子断条故障检测方法, 其特征是, 采用单工 频周期滑动窗方法计算定子电流瞬时信号的有效值的方法是 : 选取定子电流瞬。
6、时信号中的连续 20 点, 计算其有效值 ; 对于所选 取的中的连续20点, 保留后19点, 顺序递补中的后面一点, 即第21点, 从而再次获得 中的连续 20 点, 再次计算其有效值, 以此类推, 确定的有效值变化趋势。 4. 根据权利要求 3 所述笼型异步电动机转子断条故障检测方法, 其特征是, 所述定子 电流瞬时信号的采样频率设定为 1000Hz、 采样时长设定为 10 s, 从中提取的最平稳的一 段数据的时长为 4s 。 5. 根据权利要求 4 所述笼型异步电动机转子断条故障检测方法, 其特征是, 应用 ESPRIT 对待分析信号 M 进行频谱分析的方法包括以下步骤 : a. 构造相关。
7、矩阵、, 具体如下 ; 不失一般性, 待分析信号 M 可以表示为一系列余弦谐波分量之组合, 如下式所示。 , 其中, 表示采样周期 ; 表示采样点数 ; 表示谐波个数 ; 、 、 分别表示第 个 权 利 要 求 书 CN 102944842 A 2 2/3 页 3 谐波的幅值、 频率、 初相角 ; 所述信号 M 主要包括 2 个频率分量, 即p 2 ; 定义, 引入以下阶矩阵, 其中,m应远大于p, 此处取m=20 : , , 则的自相关矩阵为 , 而和的互相关矩阵为 , 其中, 表示数学期望, 表示共轭转置 ; b. 对进行特征值分解, 确定其最小特征值; c. 计算, 表示阶单位阵 ; d。
8、. 计算, 为阶阵, 此处, 表示阶单位阵 ; e. 对进行奇异值分解, 此处,由个主 奇异值组成,; f. 计算矩阵; g. 对进行广义特征值分解, 确定个广义特征值, 其余 个广义特征值恒等于 0 ; h. 根据广义特征值, 确定采样信号各个分量的频率, 、分别表示特征值的虚部、 实部 ; i. 计算矩阵; j. 计算矩阵, 此处为一列矢量, 而为列矢量 ; 权 利 要 求 书 CN 102944842 A 3 3/3 页 4 k. 确定采样信号各个分量的幅值, 。 权 利 要 求 书 CN 102944842 A 4 1/11 页 5 一种笼型异步电动机转子断条故障检测方法 技术领域 0。
9、001 本发明涉及一种能够检测笼型异步电动机转子断条故障的方法, 属电机技术领 域。 背景技术 0002 笼型异步电动机在运行过程中, 转子导条受到径向电磁力、 旋转电磁力、 离心力、 热弯曲挠度力等交变应力的作用, 加之转子制造缺陷, 可能导致断条故障, 此种故障发生概 率约为 15% 。 0003 转子断条是典型的渐进性故障, 初期通常 1、 2 根导条断裂, 而后逐渐发展以至电 机出力下降甚至停机, 因此, 必须实施转子断条故障检测。 0004 笼型异步电动机发生转子断条故障之后, 在其定子电流中将出现 (12ks)f1频率 的附加电流分量 (s为转差率,f1为供电频率,k=1, 2, 。
10、) , 可以将其作为转子断条故障特 征。而定子电流信号易于采集, 因此基于快速傅里叶变换 (FFT) 的定子电流信号频谱分析 方法被广泛应用于转子断条故障检测。 0005 最初的转子断条故障检测方法是对稳态定子电流信号直接进行 FFT 频谱分 析, 根据频谱图中是否存在 (12ks)f1频率分量判断转子有无断条。由于转子轻微 断条时, (12ks)f1分量的幅值相对于f1分量非常小, 而异步电动机运行时转差率 很小, (12ks)f1与f1这两个频率数值接近, 如果直接做 FFT 频谱分析, 则 (12ks)f1分 量可能被f1分量的泄漏所淹没。这是此方法的不足之处。 0006 为了弥补此方法。
11、之不足, 发展形成了基于希尔伯特变换的笼型异步电动机转子断 条故障检测方法, 其核心在于 : 将转子断条故障时的定子电流信号视为 2ksf1分量 (调制信 号) 对f1分量 (载波信号) 调制的结果, 而调制信号显然包含转子断条故障特征 ; 因此, 通过 适当方法进行 “解调” 而获取调制信号, 提取并分析其中的 2ksf1分量即可实现转子断条故 障检测。 0007 上述希尔伯特变换检测方法, 它包括以下步骤 : a. 测取一相定子电流瞬时信号; b. 计算定子电流瞬时信号的希尔伯特变换,t表示时间、 表示时 间延迟 ; c. 计算希尔伯特模量 Mod,; 该希尔伯特模量的频率成分非常复杂, 。
12、包括 2sf1、 4sf1、 6sf1、 8sf1频率分量以及直流分 量, 说明如下。 0008 异步电动机发生转子断条故障后, 在其定子电流中将出现 (12ks)f1频率的附加 电流分量。 因此, 异步电动机在转子断条故障情况下的定子电流信号可采用式(1)模拟。 此 处, 取k=1, 2(k2 时的附加电流分量幅值过小, 忽略不计) 。 说 明 书 CN 102944842 A 5 2/11 页 6 0009 (1) 其中,Im1、ImL、ImR、ImLL、ImRR分别代表f1、 (1-2s)f1、 (1+2s)f1、 (1-4s)f1、 (1+4s)f1分量 的幅值 ;f 1、f L、f 。
13、R、f LL、f RR分别代表各分量的初相角。 0010 推导其希尔伯特变换, 如式 (2) 所示 : (2) 进一步推导希尔伯特模量 Mod, 可得 : (3) 根据式(3)可知, 该希尔伯特模量包括2sf1、 4sf1、 6sf1、 8sf1频率分量以及直流分量, 非 常复杂。 但是, 对于实际电机而言, 存在如下数值关系 :Im1 ImL、Im1 ImR、ImL ImLL、ImRImRR。 因此, 结合式 (3) 可知, 该希尔伯特模量的最主要频率成分是直流分量、 2sf1分量、 4sf1分 量。 0011 d. 滤除希尔伯特模量 Mod 中的直流分量, 获得待分析信号 M, M=Mod。
14、-mean(Mod), mean(Mod) 表示 Mod 的平均值 ; 显然, 待分析信号 M 的主要频率成分是 2sf1分量、 4sf1分量。 0012 e. 对待分析信号 M 做 FFT 频谱分析, 根据频谱图中是否存在 2sf1分量、 4sf1分量 谱峰而进行转子断条故障检测 : 若同时存在2sf1分量、 4sf1分量谱峰, 则转子断条 ; 否则, 转 说 明 书 CN 102944842 A 6 3/11 页 7 子正常。 0013 至此可知, 从 FFT 频谱分析的角度考虑, 上述希尔伯特变换方法本质上是将定子 电流信号中的f1分量变换为直流分量, 同时将 (12ks)f1分量变换为。
15、 2ksf1分量。但是, 上述方法是以 “直流分量可以通过减去平均值而理想滤除” 为前提的。这就要求希尔 伯特模量在采集期间保持平稳, 进而要求定子电流信号、 负荷在采集期间保持平稳。 具体而 言, 式 (13) 中的Im1、ImL、ImR、ImLL、ImRR在采集期间均应恒定。如此, 希尔伯特模量中的直流 分量在采集期间才是恒定的, 方可通过减去平均值而理想滤除。 否则, 该直流分量将是波动的, 无法通过减去平均值而完全滤除。从 FFT 频谱分析的角度考 虑, 这将导致直流分量 “泄漏” 而使 FFT 频谱趋于复杂甚至混淆转子断条故障特征2sf1 分量、 4sf1分量, 影响转子断条故障检测。
16、的可靠性。 但是, 在工程实际中,“希尔伯特模量、 定 子电流信号、 负荷在采集期间保持平稳” 这一要求是无法满足的, 原因在于对于实际电 机而言, 一定程度的负荷波动是不可避免的。因此, 在工程实际中, 上述希尔伯特变换方法 欠缺可靠性。 0014 另外, 上述希尔伯特方法基于 FFT 频谱分析而受限于频率分辨力 (采样时长之倒 数) , 解释如下。 0015 在工程实际中广泛采用大型异步电动机, 并且一般在 40%60% 的负荷率下运行, 转差率s数值很小 (p) (5) (6) 则的自相关矩阵为 (7) 而和的互相关矩阵为 (8) 式 (7)、 (8) 中, 表示数学期望, 表示共轭转置。
17、。 0041 ESPRIT 算法如下 : a 构造相关矩阵、; b 对进行特征值分解, 确定其最小特征值; c 计算, 表示阶单位阵 ; d 计算, 为阶阵,( 此处, 表示阶单位阵 ) ; e 对进行奇异值分解, 此处,(由个主 奇异值组成 ),; f 计算矩阵; g 对进行广义特征值分解, 确定个广义特征值( 其余 个广义特征值恒等于 0) ; h 根据广义特征值, 确定采样信号各个分量的频率, 、分别表示特征值的虚部、 实部 ; 说 明 书 CN 102944842 A 12 9/11 页 13 i 计算矩阵; j 计算矩阵, 此处为一列矢量, 而为列矢量 ; k 确定采样信号各个分量的。
18、幅值, 。 0042 根据上述 ESPRIT 的基本概念、 步骤, 可以推断与 FFT 逐频点进行频谱分析不 同, ESPRIT 着眼于全频段、 通过信号相关矩阵的特征值分解而进行频谱分析, 这可在一定程 度上抑制频谱泄漏。并且, 与 FFT 对比而言, ESPRIT 具备原信号外推能力, 因而其频率分辨 力可以摆脱采样时长的限制, 即使针对短时信号, 亦可达到高频率分辨力。因此, 将 ESPRIT 应用于转子断条故障检测具备可行性。首先, 即使引入负荷波动, ESPRIT 亦可在一定程度 上抑制直流分量泄漏, 保证转子断条故障检测的可靠性。其次, 由于仅需采集短时信号, ESPRIT 较 F。
19、FT 具备更大可能而回避负荷波动, 这亦是有利的。 0043 应用该方法对一台实验电机进行转子断条故障检测, 效果令人满意。 0044 接线示于图 2。电机采用一台 Y100L-2 型三相异步电动机 (3kW、 380V、 50Hz) , 除 正常转子外, 另行配备一故障转子以模拟断条故障, 该故障转子存在一根断裂导条 (距端环 10mm 处钻孔, 直径 6mm、 深度 10mm) 。 0045 在实验中, 电机负荷接近空载且不做调节以期保持恒定、 转差率s 约为 0.33% 。 为了尽可能回避负荷波动, 尽量缩短信号采集时长, 选择为4s 。 这样, 信号采集时长大于待 提取转子断条故障特征。
20、2s f1分量的一个持续周期。此时,2s f10.33Hz ( f150Hz), 对应 的一个持续周期约为 3s 。 0046 图 3、 图 4 表示电机在转子正常与故障时的定子电流有效值时变曲线。 0047 图 5、 图 6 表示电机在转子正常与故障时的定子电流希尔伯特模量 FFT 频谱, 具体 数据示于表 1, 这是采用现有的希尔伯特变换方法所获得的。 0048 图7、 图8表示电机在转子正常与故障时的定子电流希尔伯特模量ESPRIT频谱, 具 体数据示于表 2, 这是采用本发明所获得的。 0049 图9、 图10表示信号采集时长10s情况下, 电机在转子正常与故障时的定子电流希 尔伯特模。
21、量 FFT 频谱, 具体数据示于表 3, 这是采用现有的希尔伯特变换方法所获得的。 0050 根据图 3、 图 4 可知 : 在转子正常与故障情况下, 尽管在实验过程中, 电机负荷接近 空载且不做调节以力图保持恒定, 但定子电流有效值时变曲线均包含一定程度的波动, 这 说明对于实际电机而言, 一定程度的负荷波动是不可避免的, 当前异步电动机转子断 条故障检测的希尔伯特变换方法的缺陷即来源于此。 0051 对比图 5、 图 7 并结合表 1、 表 2, 可以发现 : 在正常情况下, FFT 频谱中即包含故障 特征谱峰2sf1、 4sf1分量, 这将导致误判 “转子断条” ; 而 ESPRIT 频。
22、谱中并无故障特征谱 峰, 故可排除误判 “转子断条” 的可能, 这是一显著进步。 0052 对比图 6、 图 8 并结合表 1、 表 2, 可以发现 : 在断条情况下, FFT 频谱中出现明显的 说 明 书 CN 102944842 A 13 10/11 页 14 故障特征谱峰2sf1分量, 但另一故障特征谱峰 4sf1分量并不存在, 这将导致 “转子正 常” 的错误判断, 这归因于采用 4s 信号, FFT 频谱频率分辨力 (0.25Hz) 不足以分辨转 子断条故障特征 ; 但是, 在ESPRIT频谱中, 故障特征谱峰2sf1、 4sf1分量均非常明显, 据 此即可实现转子断条故障的可靠检测。
23、。 0053 另外, 根据图 5、 图 6 并结合表 1, 可知 : 对于仅持续 4s 的电机信号, 现有的希尔伯 特变换方法无法正确检测转子断条故障。延长信号采集时间至 10s, 该问题得以解决, 参阅 图 9、 图 10 与表 3。此时, FFT 频谱频率分辨力提高至 0.1Hz, 图 10 中包含明显的故障特征 谱峰2sf1、 4sf1分量。但是, 对于正常转子, FFT 频谱中亦包含故障特征谱峰2sf1、 4sf1分量, 参阅图 9, 同样可能导致 “转子断条” 的错误检测结果。 0054 综上可知, 当前异步电动机转子断条故障检测的希尔伯特变换方法存在局限 负荷波动导致解调信号中直流。
24、分量波动、 泄漏, 致使 FFT 频谱趋于复杂而影响转子断条故 障检测的可靠性。由于 FFT 频谱分析受限于频率分辨力而需要足够时长的电机信号, 意味 着以更高概率引入负荷波动, 致使上述问题愈发突出。 特别是, 对于工程实际中广泛采用的 大型异步电动机, 上述希尔伯特变换方法面临严峻挑战, 甚至失效。本发明引入 ESPRIT 以 克服上述局限, 并且仅需短时信号即可高可靠性地检测转子断条故障, 适用于工程实际中 广泛采用的大型异步电动机, 这是本发明的显著优势。 0055 相对于当前的希尔伯特变换方法, 本发明做出了重大改进, 说明如下。 0056 首先, 当前的希尔伯特变换方法是以 FFT。
25、 频谱分析为基础的, 本发明则颠覆了这 一基础而采用 ESPRIT 频谱分析。 0057 其次, 当前的希尔伯特变换方法基于 FFT 频谱分析而受限于负荷波动 ( 这在工程 实际中是不可避免的 )、 直流分量泄漏以及频率分辨力, 在异步电动机低转差率运行情况下 因需要持续采集足够时长的电机信号而失效 ; 本发明则基于 ESPRIT 频谱分析而摆脱了上 述限制, 在异步电动机低转差率运行情况下亦可高可靠性地检测转子断条故障, 因而适用 于工程实际中广泛采用的大型异步电动机 ( 低转差率运行 ), 具备重大工程价值、 广阔应用 前景。 0058 表 1 定子电流希尔伯特模量 FFT 频谱数据 00。
26、59 表 2 定子电流希尔伯特模量 ESPRIT 频谱数据 说 明 书 CN 102944842 A 14 11/11 页 15 0060 表 3 定子电流希尔伯特模量 FFT 频谱数据 (信号采集时间 10s) 。 说 明 书 CN 102944842 A 15 1/5 页 16 图 1 图 2 说 明 书 附 图 CN 102944842 A 16 2/5 页 17 图 3 图 4 说 明 书 附 图 CN 102944842 A 17 3/5 页 18 图 5 图 6 说 明 书 附 图 CN 102944842 A 18 4/5 页 19 图 7 图 8 说 明 书 附 图 CN 102944842 A 19 5/5 页 20 图 9 图 10 说 明 书 附 图 CN 102944842 A 20 。