一种结构参数未知条件下的感应电机转速估计方法技术

本发明专利技术提供一种结构参数未知条件下的感应电机转速估计方法，所述方法包括：采集感应电动机的单相定子电流，将其傅里叶系数向量中最大值对应的频率识别为供电频率；将傅里叶系数向量在频域上展开，检测所有对应的局部频率峰值；基于检测得到的局部频率峰值向量，构造转速调制频率向量、齿槽调制频率向量；构造比率矩阵；初步筛选满足调制机理的转速调制频率与齿槽调制频率的组合；根据初步筛选的转速调制频率，计算对应的电机转频；根据电机转频的谐波特性，进一步筛选电机转频；将筛选得到的电机转频中的最小值，作为电机转速的最终估计值。本发明专利技术无需任何电机结构参数，且能够有效区分电机转频以及可能存在的机械故障频率，具有广泛的适用性。有广泛的适用性。有广泛的适用性。

【技术实现步骤摘要】
一种结构参数未知条件下的感应电机转速估计方法


[0001]本专利技术涉及感应电机转速估计
，尤其涉及一种结构参数未知条件下的感应电机转速估计方法。

技术介绍

[0002]感应电动机具有结构简单、可靠性高、成本较低等优势，作为动力源被广泛应用于诸多工业装备、运载工具、市政设施的传动链中，其安全稳定运行对于保障设备功能与质量至关重要。电机转速信息是辨识故障频率特征的基础，是有效开展传动链监测诊断的必备条件。然而，感应电动机的滑差率与转速通常是未知的，且额外加装转速传感器将增加总体成本并降低整体可靠性。无转速传感器的感应电机转速估计，是制约电传动链监测诊断技术发展的关键难题，一直以来备受工业界关注。
[0003]现有感应电机转速估计方法包含基于模型和观测器的方法、基于电流频域特征分析的方法两类。其中，基于模型和观测器的方法不仅需要同时测量供电电压与电流，还需要知道准确的电机内部结构参数，如定子绕组的电阻值。然而，这些参数不仅难以获取，而且可能随环境温度变化而变化，所产生的偏移对转速估计造成较大误差。因此，这类方法主要应用在电机控制系统的设计阶段中，对于内部物理结构参数未知的传动链而言，难以直接用于监测诊断。基于电流频域特征分析的转速估计方法以定子电流中包含的机械调制特征为基础，通过特征频率辨识反推转速信息。其优点是不依赖时变结构参数，且所需的单相定子电流容易从供电线路或控制系统中获取。因此，基于电流频域特征分析的转速估计方法在传动链监测诊断领域具有较大潜力。
[0004]中国专利技术专利CN 115586345 B公布了一种基于定子电流频域特征分析的转速估计方法。首先，根据定子电流信号频谱，识别幅值最大的频率成分为供电频率；其次，对定子电流信号进行低通滤波并计算幅值包络谱；再次，根据铭牌上的额定转速、基于已知极对数和供电频率计算得到的同步转速、人为给定的劣化因子转速，制定转速识别范围；最后，根据转速识别范围，从幅值包络谱中识别对应范围内的谱峰，作为转速频率。
[0005]上述方法的核心是在频域上识别电机转速对供电频率的调制作用。在工程实际中，转轴偏心/不对中或磁场不对称常常难以避免，因此转速调制作用通常容易被识别。然而，该方法对实际定子电流结构做了过于简单的假设。实际传动链中，尤其是在包含机械故障的状态下，存在周期性的扭矩波动，该波动频率同样对供电频率造成调制作用。因此，电流信号的幅值包络谱中通常包含多个特征频率谱峰。由于感应电机的滑移率和转速在不同负载条件下分布范围较广，因此无法将转速频率和故障特征频率区分开，容易得到错误的转速估计结果。此外，该方法需要已知电机极对数或特定工况下的额定转速，但在实际应用中常常无法获取。

技术实现思路

[0006]针对上述问题，本专利技术的目的在于提供一种结构参数未知条件下的感应电机转速
估计方法，能够在不提供任何电机结构参数的条件下，仅基于电流信号的频率结构，区分转速调制和其它特征频率调制，从而实现转速的准确辨识。
[0007]为解决上述技术问题，本专利技术提供如下技术方案：
[0008]一种结构参数未知条件下的感应电机转速估计方法，包括以下步骤：
[0009]S101、采集感应电动机的单相定子电流，将其傅里叶系数向量中最大值对应的频率识别为供电频率；
[0010]S102、将所得傅里叶系数向量在频域上展开，检测所有对应的局部频率峰值，得到局部频率峰值向量；
[0011]S103、基于检测得到的局部频率峰值向量，构造转速调制频率向量、齿槽调制频率向量；
[0012]S104、根据转速调制频率向量和齿槽调制频率向量，构造比率矩阵；
[0013]S105、根据比率矩阵的限制条件，初步筛选满足调制机理的转速调制频率与齿槽调制频率的组合；
[0014]S106、根据初步筛选的转速调制频率，计算对应的电机转频；
[0015]S107、根据电机转频的谐波特性，进一步筛选电机转频；
[0016]S108、将筛选得到的电机转频中的最小值，作为电机转速的最终估计值。
[0017]优选地，步骤S101包括：
[0018]采集感应电动机的单相定子电流信号x[t]；
[0019]对采集的单相定子电流信号x[t]加汉宁窗后进行离散傅里叶变换，得到傅里叶系数向量X[f]；其中f＝[F
s
/N:F
s
/N:F
s
]是与傅里叶系数对应的频率向量，F
s
代表电流信号采样频率，N代表电流信号数据点长度，t代表时间；
[0020]识别傅里叶系数向量X[f]中的最大值，将其对应的频率f识别为供电频率f
supply
，即
[0021]f
supply
＝f|
X[f]＝max{X[f]}
。
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)
[0022]优选地，步骤S102包括：
[0023]从傅里叶系数向量中提取所有局部频率峰值，得到局部频率峰值向量f
peak
；具体方法为：
[0024]a)对傅里叶系数向量X[f]求导，得到导数向量X
’
[f]＝X[f]‑
X[f
‑
1]，X[0]＝NaN；
[0025]b)对导数向量X
’
[f]进行符号标定，得到符号导数向量SX
’
[f][0026][0027]c)对符号导数向量SX
’
[f]求导，得到二阶导数向量X”[f]＝SX
’
[f]‑
SX
’
[f
‑
1]，SX
’
[0]＝0；
[0028]d)搜索二阶导数向量X”[f]中值为
‑
2的成分，将其对应的f作为谱峰频率，得到局部频率峰值向量f
peak
＝[f
peak1
,f
peak2
,...f
peakm
...f
peakM
]，其中m＝1,2,...M，M代表搜索得到的谱峰频率的个数。
[0029]优选地，步骤S103包括：
[0030]从局部频率峰值向量f
peak
中挑选转速调制频率向量f
rotate
和齿槽调制频率向量f
slot
，如下：
[0031][0032][0033]优选地，步骤S104包括：
[0034]令转速调制频率向量f
rotate
的长度为A，齿槽调制频率向量f
slot
的长度为B，计算比率矩阵
[0035]R(a,b)＝[f
slot
(b)
‑
f
supply
]/[f
rotate
(a)
‑
f
supply...

【技术保护点】

【技术特征摘要】
1.一种结构参数未知条件下的感应电机转速估计方法，其特征在于，包括以下步骤：S101、采集感应电动机的单相定子电流，将其傅里叶系数向量中最大值对应的频率识别为供电频率；S102、将所得傅里叶系数向量在频域上展开，检测所有对应的局部频率峰值，得到局部频率峰值向量；S103、基于检测得到的局部频率峰值向量，构造转速调制频率向量、齿槽调制频率向量；S104、根据转速调制频率向量和齿槽调制频率向量，构造比率矩阵；S105、根据比率矩阵的限制条件，初步筛选满足调制机理的转速调制频率与齿槽调制频率的组合；S106、根据初步筛选的转速调制频率，计算对应的电机转频；S107、根据电机转频的谐波特性，进一步筛选电机转频；S108、将筛选得到的电机转频中的最小值，作为电机转速的最终估计值。2.根据权利要求1所述的感应电机转速估计方法，其特征在于，步骤S101包括：采集感应电动机的单相定子电流信号x[t]；对采集的单相定子电流信号x[t]加汉宁窗后进行离散傅里叶变换，得到傅里叶系数向量X[f]；其中f＝[F
s
/N:F
s
/N:F
s
]是与傅里叶系数对应的频率向量，F
s
代表电流信号采样频率，N代表电流信号数据点长度，t代表时间；识别傅里叶系数向量X[f]中的最大值，将其对应的频率f识别为供电频率f
supply
，即f
supply
＝f|
X[f]＝max{X[f]}
。
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)3.根据权利要求2所述的感应电机转速估计方法，其特征在于，步骤S102包括：从傅里叶系数向量中提取所有局部频率峰值，得到局部频率峰值向量f
peak
；具体方法为：a)对傅里叶系数向量X[f]求导，得到导数向量X
’
[f]＝X[f]
‑
X[f
‑
1]，X[0]＝NaN；b)对导数向量X
’
[f]进行符号标定，得到符号导数向量SX
’
[f]c)对符号导数向量SX
’
[f]求导，得到二阶导数向量X”[f]＝SX
’
[f]
‑
SX
’
[f
‑
1]，SX
’
[0]＝0；d)搜索二阶导数向量X”[f]中值为
‑
2的成分，将其对应的f作为谱峰频率，得到局部频率峰值向量f
peak
＝[f
peak1
,f
peak2
,...f
peakm
...f
peakM
]，其...

【专利技术属性】
技术研发人员：陈小旺，冯志鹏，
申请(专利权)人：北京科技大学，
类型：发明
国别省市：