CN106909076B - 一种pmsm永磁体局部退磁故障程度评估与数值化描述方法 - Google Patents
一种pmsm永磁体局部退磁故障程度评估与数值化描述方法 Download PDFInfo
- Publication number
- CN106909076B CN106909076B CN201710110932.5A CN201710110932A CN106909076B CN 106909076 B CN106909076 B CN 106909076B CN 201710110932 A CN201710110932 A CN 201710110932A CN 106909076 B CN106909076 B CN 106909076B
- Authority
- CN
- China
- Prior art keywords
- pmsm
- permanent magnet
- fault
- degree
- local demagnetization
- 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.)
- Active
Links
- 230000005347 demagnetization Effects 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000000605 extraction Methods 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 230000003044 adaptive effect Effects 0.000 claims abstract description 12
- 238000011156 evaluation Methods 0.000 claims abstract description 12
- 238000005070 sampling Methods 0.000 claims abstract description 10
- 238000013507 mapping Methods 0.000 claims abstract description 5
- 238000012512 characterization method Methods 0.000 claims description 23
- 230000008569 process Effects 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 5
- 238000011478 gradient descent method Methods 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000001914 filtration Methods 0.000 claims description 2
- 230000008030 elimination Effects 0.000 claims 1
- 238000003379 elimination reaction Methods 0.000 claims 1
- 230000008859 change Effects 0.000 description 9
- 230000009467 reduction Effects 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 4
- 230000007423 decrease Effects 0.000 description 3
- 230000004907 flux Effects 0.000 description 3
- 230000001360 synchronised effect Effects 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000005389 magnetism Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 230000001133 acceleration Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 238000004141 dimensional analysis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 210000003739 neck Anatomy 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000010349 pulsation Effects 0.000 description 1
- 238000010183 spectrum analysis Methods 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
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Abstract
本发明公开了一种PMSM永磁体局部退磁故障程度评估与数值化描述方法,其步骤为:利用电流采样电路采集PMSM定子电流并滤除逆变器高频谐波;采用自适应基波提取算法提取采集的PMSM定子电流的基波;计算表征PMSM永磁体局部退磁程度的微弱故障特征谐波的分形盒维数;分析PMSM驱动系统运行工况对故障特性谐波分形盒维数计算结果的影响,建立故障程度与其故障特征谐波分形盒维数之间的唯一映射关系。本发明以故障特征谐波为分析对象,消除了逆变器高频谐波及幅值较大的电流基波对微弱故障特征谐波的影响,实现故障特征谐波分形盒维数的准确计算,实现了独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。
Description
技术领域
本发明涉及永磁同步电机(Permanent Magnet Synchronous Motor,PMSM)技术领域,具体涉及一种PMSM永磁体局部退磁故障程度评估与数值化描述方法。
背景技术
永磁同步电机(PMSM)具有结构简单、故障率低和运行效率高等技术优势,在工业伺服驱动、新能源汽车以及新能源发电等领域获得了广泛应用。然而,对于上述应用领域而言,PMSM功率密度要求高,散热条件差,运行工况复杂,多处于加减速、变负荷运行以及最大转矩电流比与弱磁联合控制下,存在较强的电枢反应,加之自然老化等因素,容易出现永磁体均匀退磁或局部退磁故障,导致电机输出转矩减少及转矩脉动,直接影响电驱动系统的转矩控制精度与运行可靠性。
针对PMSM永磁体退磁故障的研究目前多集中于基于数据驱动、模型驱动及高频信号注入方法的故障定性诊断领域,而对其故障程度的评估及有效描述方面的研究较为少见。鉴于PMSM永磁体均匀退磁故障并没有打破其转子等效物理结构的对称性,仅将导致永磁体磁链幅值的降低。因此,基于PMSM动态数学模型并采用动态数据处理技术即可实现永磁体磁链幅值的在线观测,进而实现永磁体均匀退磁故障程度的准确评估。而PMSM永磁体局部退磁故障不仅会降低永磁体磁链幅值,亦会因转子等效物理结构对称性的破坏而在PMSM定子电流中出现特定次故障特征谐波ffault=fs(1±k/p)(式中,fs为定子电流基波频率,p为PMSM极对数,k取正整数),从而导致PMSM电磁转矩的降低及其驱动系统控制性能的下降。因此,应从上述两个维度对PMSM永磁体局部退磁故障程度进行评估。
鉴于永磁体局部退磁故障所导致的PMSM定子电流中的故障特征谐波能够对PMSM电磁转矩的降低及其驱动系统性能的下降进行综合表征,故可将上述故障特征谐波作为分析对象,通过计算表征其复杂程度的分形盒维数实现PMSM永磁体局部退磁故障程度的评估与数值化表征,以便采用积极有效的容错控制方法,确保PMSM驱动系统的安全可靠运行。对于此思路,目前尚未发现有效的文献报道。
鉴于此,PMSM永磁体局部退磁故障程度的准确评估与数值化描述已成为确保PMSM驱动系统安全可靠运行亟待解决的关键技术。
发明内容
针对PMSM永磁体局部退磁故障导致PMSM电磁转矩的降低及其驱动系统控制性能下降的技术问题,本发明公开一种PMSM永磁体局部退磁故障程度评估与数值化描述方法,以故障特征谐波为分析对象,基于自适应基波提取算法与分形盒维数相结合的方法,消除逆变器谐波及幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。
为了解决上述技术问题,本发明的技术方案是:一种PMSM永磁体局部退磁故障程度评估与数值化描述方法,其特征在于,以故障特征谐波为分析对象,基于自适应基波提取算法与分形盒维数相结合的方法,消除逆变器谐波及幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响,实现了独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述,其步骤如下:
步骤一:利用电流采样电路采集PMSM定子电流,并消除逆变器高频谐波对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响;
步骤二:采用自适应基波提取算法提取采集的PMSM定子电流的基波,消除幅值较大的基波电流对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响;
步骤三:计算表征PMSM永磁体局部退磁程度的微弱故障特征谐波的分形盒维数;
步骤四:分析PMSM驱动系统运行工况对故障特征信号分形盒维数计算结果的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的评估与数值化描述。
所述采集PMSM定子电流并滤除逆变器高频谐波的方法是:在PMSM电流内环、转速外环双闭环驱动系统中采用电流采样电路采集PMSM定子电流,并通过低通滤波电路滤除逆变器高频谐波,消除其对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响。
所述的自适应基波电流提取算法是:设从PMSM定子电流i(t)中提取到的基波电流为iext(t),定义代价函数为式中,θ为表征所提取基波电流的幅值、频率和相位瞬时值的参数向量;采用梯度下降法使代价函数J(t,θ)收敛于最小值的未知参数向量θ的调整方法为:该收敛过程可以生成一组表征基波电流的幅值、频率及相位瞬时值提取过程的非线性微分方程,表达式为:
式中,I(t)、ω(t)、φ(t)分别表示提取基波电流iext(t)的幅值、频率及相位的瞬时值,e(t)为提取误差,μ1、μ2、μ3为正常数;对上述微分方程进行离散处理并求解,即可实现PMSM定子基波电流的提取,消除幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响。
所述分形盒维数的计算方法为:设离散信号Y是n维欧式空间Rn上的闭集,用尽可能细的ε网格划分Rn,Nε是集合Y的网格计数;以ε尺寸网格作为基准,逐步放大到kε尺寸网格,其中k取正整数;令Nkε为离散空间上集合Y的网格计数,则有:
式中,i=1,2,...,N/k,N为采样点数,k=1,2,...,M,M<N;yk(i-1)+1,yk(i-1)+2,...yk(i-1)+k+1表示离散空间上的集合Y的样本,则网格计数Nkε=P(kε)/ε+1,P(kε)为上述公式的计算结果;在lg(kε)-lgNkε双对数图中确定线性度较好的一段直线,令其起点和终点分别为k1和k2,则存在lgNkε=αlg(kε)+b,k1≤k≤k2,采用最小二乘法确定该直线的斜率,即获得离散信号的分形盒维数,表达式为:
所述分析PMSM驱动系统运行工况对故障特征信号分形盒维数计算结果的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的评估与数值化描述的方法是:分析PMSM驱动系统负载及运行速度变化对故障特征谐波分形盒维数计算结果的影响,确定不同永磁体局部退磁故障程度与其故障特征谐波分形盒维数的唯一映射关系,并存储为数据表格,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。
针对PMSM永磁体局部退磁故障程度的评估与数值化描述问题,本发明以PMSM定子电流中的故障特征谐波为分析对象,采用滤波电路及基于自适应基波提取算法与故障特征谐波分形盒维数相结合的方法消除逆变器高频谐波及幅值较大的基波电流对微弱故障特征谐波的影响,实现表征PMSM永磁体局部退磁故障程度的故障特征谐波分形盒维数的准确计算;同时,通过分析PMSM驱动系统不同运行工况对表征永磁体局部退磁程度的故障特征谐波分形盒维数计算结果的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明的PMSM驱动系统结构框图。
图2为本发明的自适应基波提取算法结构框图。
图3为本发明设定参考转速750rpm、负载转矩50N·m的PMSM驱动系统运行工况下、单个永磁体50%局部退磁时,PMSM定子电流波形为图3(a)、基波提取电流波形图为图3(b)及故障特征谐波分形盒维数的线性区双对数图为图3(c)。
图4为本发明PMSM驱动系统动态运行工况下,设定负载转矩50N·m,转速动态如图4(a)所示;单个永磁体50%局部退磁程度的PMSM定子电流波形如图4(b)所示、基波提取电流波形如图4(c)所示及故障特征谐波分形盒维数的线性区双对数图如图4(d)所示。
图5为本发明单个永磁体局部退磁50%、参考转速750转/分、负载转矩从20N·m到100N·m变化时的故障特征谐波分形盒维数变化曲线。
图6为本发明单个永磁体退磁50%、负载转矩50N·m、转速从300转/分到900转/分变化时的故障特征谐波分形盒维数变化曲线。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
一种PMSM永磁体局部退磁故障程度评估与数值化描述方法,以故障特征谐波为分析对象,采用滤波电路及基于自适应基波提取算法与分形盒维数相结合的方法,消除逆变器高频谐波及幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述;步骤如下:
步骤一、利用电流采样电路采集PMSM定子电流,并采用滤波电路消除逆变器高频谐波对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响。
具体地,根据图1所示的PMSM电流内环、转速外环双闭环驱动系统的结构框图,采用电流采样电路采集PMSM定子电流,并通过低通滤波电路消除逆变器高频谐波,消除其对表征PMSM永磁体局部退磁程度的微弱故障特征谐波分形盒维数计算值的影响;PMSM双闭环驱动系统采用转速控制模式,转速偏差经PI调节后用于产生指令电流,该电流指令值与实际采样电流比较,产生误差信号并经PI调节后生成SVPWM信号,控制逆变器工作。
步骤二、采用自适应基波提取算法提取采集的PMSM定子电流的基波,消除幅值较大的基波电流对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响。
具体地,定义i(t)为PMSM定子电流,t为时间,其包括基波电流io(t)和各种谐波电流i1(t),且i(t)=io(t)+i1(t)。定义从PMSM定子电流i(t)中提取到的基波成份为iext(t),提取过程中采用梯度下降法减小PMSM定子电流i(t)与所提取基波电流iext(t)之间的最小平方误差。PMSM实际的定子基波电流与所提取的基波电流不能采用相同的符号,二者之间存在误差,采用梯度下降法可以减小该误差但不能完全消除。定义代价函数为式中,θ为表征所提取基波电流的幅值I(t)、频率ω(t)和相位φ(t)瞬时值的参数向量。梯度下降法提供了一种使代价函数J(t,θ)收敛于最小值的未知参数向量θ的调整方法,调整过程由式描述,该收敛过程可以生成一组表征基波电流的幅值、频率及相位瞬时值提取过程的非线性微分方程,表达式为:
式中,I(t)、ω(t)、φ(t)分别代表提取基波电流iext(t)的幅值、频率及相位的瞬时值,e(t)为提取误差,μ1、μ2、μ3为正常数。
对上述微分方程进行离散处理并求解,即可实现PMSM定子电流基波的提取,消除幅值较大的基波电流对PMSM永磁体局部退磁微弱故障特征谐波的影响,实现表征永磁体局部退磁故障程度的微弱故障特征信号分形盒维数的准确计算,基波提取算法框图如图2所示。
步骤三、计算表征PMSM永磁体局部退磁程度的微弱故障特征谐波的分形盒维数。
具体地,设X为n维欧式空间Rn的非空有界子集,N(X,ε)表示最大直径为ε且能够覆盖X集合的最少个数,定义X的分形盒维数dimBX为:dimBX为一定义量,为X的分形盒维数)。由于定义式中的极限无法按照定义求出,故采用近似方法计算离散信号盒维数。设离散信号Y是n维欧式空间Rn上的闭集,用尽可能细的ε网格划分Rn,Nε是集合Y的网格计数。以ε尺寸网格作为基准,逐步放大到kε尺寸网格,其中k取正整数,令Nkε为离散空间上的集合Y的网格计数,则有:
式中,i=1,2,···,N/k,N为采样点数,k=1,2,···,M,M<N;yk(i-1)+1,yk(i-1)+2,……,yk(i-1)+k+1表示离散空间上的集合Y的样本。则网格计数Nkε为:Nkε=P(kε)/ε+1,P(kε)为上述公式的计算结果。在lg(kε)-lgNkε双对数图中确定线性度较好的一段直线,令其起点和终点分别为k1和k2,则存在lgNkε=αlg(kε)+b,k1≤k≤k2,采用最小二乘法确定该直线的斜率,即获得离散信号的分形盒维数,表达式为:
用于计算表征永磁体局部退磁程度的故障特征信号的分形盒维数。据此计算不同永磁体局部退磁程度的故障特征谐波的分形盒维数,即可建立不同永磁体局部退磁程度与其故障特征谐波分形盒维数之间的映射关系。
步骤四:分析PMSM驱动系统运行工况对故障特征信号分形盒维数计算结果的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的评估与数值化描述。
具体地,分析PMSM驱动系统负载及运行速度变化对故障特征谐波分形盒维数计算结果的影响,确定永磁体局部退磁故障程度与其故障特征谐波分形盒维数的唯一映射关系,并存储为数据表格,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。所建立的映射关系即为描述永磁体局部退磁程度与故障特征信号分形盒维数计算值的之间关系的数据表格,实际应用中通过查询此表格即可反向获取永磁体局部退磁故障程度。
对本发明进行仿真验证,验证思路为:首先建立具有单个永磁体局部退磁故障的PMSM有限元模型,再对PMSM空载时永磁体径向气隙磁通密度进行频谱分析和磁链计算,根据计算结果建立永磁体局部退磁故障的PMSM数学模型及系统仿真模型,基于此系统仿真模型进行仿真验证,电机设计参数如表1所示。
表1 电机设计参数
图3为本发明设定参考转速750rpm、负载转矩50N·m的PMSM驱动系统运行工况下,单个永磁体50%局部退磁时,PMSM定子电流波形如图3(a)所示,基波提取电流波形如图3(b)所示,故障特征谐波分形盒维数的线性区双对数图如图3(c)所示,该系统运行工况下表征不同永磁体退磁程度的故障特征谐波的分形盒维数计算结果如表2所示。
图4为本发明PMSM驱动系统动态运行工况下,设定负载转矩50N·m,转速动态如图4a所示,单个永磁体50%局部退磁程度的PMSM定子电流波形如图4(b)所示,基波提取电流波形如图4(c)所示,故障特征谐波分形盒维数的线性区双对数图如图4(d)所示,该系统运行工况下表征不同永磁体退磁程度的故障特征谐波的分形盒维数计算结果如表3所示。
表2 不同永磁体退磁程度的故障特征谐波的分形盒维数
表3 不同永磁体退磁程度的故障特征谐波的分形盒维数
基于上述仿真结果,本发明从负载转矩及电机转速变化两个维度分析基于故障特征谐波分形盒维数,实现PMSM永磁体局部退磁故障程度评估的可行性。图5为单个永磁体局部退磁50%,电机转速设定为750转/分,负载转矩从20N·m到100N·m变化时的故障特征谐波分形盒维数变化曲线,图6为单个永磁体局部退磁50%,负载转矩设定为50N·m,转速从300转/分到900转/分变化时的故障特征谐波分形盒维数变化曲线。由图5和图6知,相同故障程度下的故障特征谐波分形盒维数受负载变化的影响较小,而受转速变化的影响较为明显,其原因在于系统运行速度的变化会导致PMSM定子电流基波及各故障特征谐波频率发生改变,导致用来衡量故障特征谐波波形复杂程度的分形盒维数发生变化;而负载变化时,系统稳定运行条件下的定子电流基波及故障特征谐波只会等比缩放,波形结构并不发生改变,故不会导致表征永磁体局部退磁程度的故障特征谐波分形盒维数计算结果的变化。
因此,以故障特征谐波为分析对象,采用本发明的方法,可以消除逆变器谐波及幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响,准确计算其分形盒维数,在PMSM驱动系统定速度点下实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (5)
1.一种PMSM永磁体局部退磁故障程度评估与数值化描述方法,其特征在于,以故障特征谐波为分析对象,基于自适应基波提取算法与分形盒维数相结合的方法,消除逆变器谐波及幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响,实现了独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述,其步骤如下:
步骤一:利用电流采样电路采集PMSM定子电流,并消除逆变器高频谐波对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响;
步骤二:采用自适应基波提取算法提取采集的PMSM定子电流的基波,消除幅值较大的基波电流对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响;
步骤三:计算表征PMSM永磁体局部退磁程度的微弱故障特征谐波的分形盒维数;
步骤四:分析PMSM驱动系统运行工况对故障特征信号分形盒维数计算结果的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的评估与数值化描述。
2.根据权利要求1所述的PMSM永磁体局部退磁故障程度评估与数值化描述方法,其特征在于,所述采集PMSM定子电流并滤除逆变器高频谐波的方法是:在PMSM电流内环、转速外环双闭环驱动系统中采用电流采样电路采集PMSM定子电流,并通过低通滤波电路滤除逆变器高频谐波,消除其对表征PMSM永磁体局部退磁程度的微弱故障特征谐波的影响。
3.根据权利要求1所述的PMSM永磁体局部退磁故障程度评估与数值化描述方法,其特征在于,所述的自适应基波电流提取算法是:设从PMSM定子电流i(t)中提取到的基波电流为iext(t),定义代价函数为式中,θ为表征所提取基波电流的幅值、频率和相位瞬时值的参数向量;采用梯度下降法使代价函数J(t,θ)收敛于最小值的未知参数向量θ的调整方法为:该收敛过程可以生成一组表征基波电流的幅值、频率及相位瞬时值提取过程的非线性微分方程,表达式为:
式中,I(t)、ω(t)、φ(t)分别表示提取基波电流iext(t)的幅值、频率及相位的瞬时值,e(t)为提取误差,μ1、μ2、μ3为正常数;对上述微分方程进行离散处理并求解,即可实现PMSM定子基波电流的提取,消除幅值较大的基波电流对表征永磁体局部退磁程度的微弱故障特征谐波的影响。
4.根据权利要求1所述的PMSM永磁体局部退磁故障程度评估与数值化描述方法,其特征在于,所述分形盒维数的计算方法为:设离散信号Y是n维欧式空间Rn上的闭集,用尽可能细的ε网格划分Rn,Nε是集合Y的网格计数;以ε尺寸网格作为基准,逐步放大到kε尺寸网格,其中k取正整数;令Nkε为离散空间上集合Y的网格计数,则有:
式中,i=1,2,···,N/k,N为采样点数,k=1,2,···,M,M<N;yk(i-1)+1,yk(i-1)+2,···yk(i-1)+k+1表示离散空间上的集合Y的样本,则网格计数Nkε=P(kε)/ε+1,P(kε)为上述公式的计算结果;在lg(kε)-lgNkε双对数图中确定线性度较好的一段直线,令其起点和终点分别为k1和k2,则存在lgNkε=αlg(kε)+b,k1≤k≤k2,采用最小二乘法确定该直线的斜率,即获得离散信号的分形盒维数,表达式为:
5.根据权利要求1所述的PMSM永磁体局部退磁故障程度评估与数值化描述方法,其特征在于,所述分析PMSM驱动系统运行工况对故障特征信号分形盒维数计算结果的影响,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的评估与数值化描述的方法是:分析PMSM驱动系统负载及运行速度变化对故障特征谐波分形盒维数计算结果的影响,确定不同永磁体局部退磁故障程度与其故障特征谐波分形盒维数的唯一映射关系,并存储为数据表格,实现独立于PMSM驱动系统负载变化的永磁体局部退磁故障程度的准确评估与数值化描述。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710110932.5A CN106909076B (zh) | 2017-02-28 | 2017-02-28 | 一种pmsm永磁体局部退磁故障程度评估与数值化描述方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710110932.5A CN106909076B (zh) | 2017-02-28 | 2017-02-28 | 一种pmsm永磁体局部退磁故障程度评估与数值化描述方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106909076A CN106909076A (zh) | 2017-06-30 |
CN106909076B true CN106909076B (zh) | 2019-08-06 |
Family
ID=59208128
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710110932.5A Active CN106909076B (zh) | 2017-02-28 | 2017-02-28 | 一种pmsm永磁体局部退磁故障程度评估与数值化描述方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106909076B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115079575B (zh) * | 2022-07-19 | 2023-03-03 | 北京理工大学 | 一种基于有限控制集模型预测控制的永磁同步电机模拟器 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003235286A (ja) * | 2002-02-13 | 2003-08-22 | Nissan Motor Co Ltd | 同期機の制御装置 |
CN103454585A (zh) * | 2013-08-27 | 2013-12-18 | 河海大学 | 一种基于失电残压的永磁同步电机失磁故障诊断方法 |
CN104462762A (zh) * | 2014-11-04 | 2015-03-25 | 西南交通大学 | 一种输电线路的模糊故障分类方法 |
CN105974312A (zh) * | 2016-06-06 | 2016-09-28 | 安徽大学 | 永磁同步电机退磁故障诊断方法 |
-
2017
- 2017-02-28 CN CN201710110932.5A patent/CN106909076B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003235286A (ja) * | 2002-02-13 | 2003-08-22 | Nissan Motor Co Ltd | 同期機の制御装置 |
CN103454585A (zh) * | 2013-08-27 | 2013-12-18 | 河海大学 | 一种基于失电残压的永磁同步电机失磁故障诊断方法 |
CN104462762A (zh) * | 2014-11-04 | 2015-03-25 | 西南交通大学 | 一种输电线路的模糊故障分类方法 |
CN105974312A (zh) * | 2016-06-06 | 2016-09-28 | 安徽大学 | 永磁同步电机退磁故障诊断方法 |
Non-Patent Citations (1)
Title |
---|
电动汽车PMSM退磁故障诊断及容错控制;陈涛;《万方学位论文》;20170523;全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106909076A (zh) | 2017-06-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ullah et al. | A comprehensive review of winding short circuit fault and irreversible demagnetization fault detection in PM type machines | |
Romeral et al. | Modeling of surface-mounted permanent magnet synchronous motors with stator winding interturn faults | |
Zhang et al. | A novel detection method of motor broken rotor bars based on wavelet ridge | |
CN106597276B (zh) | 一种pmsm永磁体退磁故障诊断与故障模式识别方法 | |
CN102510263B (zh) | 基于抛载试验和数值差分的同步发电机实用参数辨识方法 | |
Wang et al. | A review of Permanent Magnet Synchronous Motor fault diagnosis | |
CN102998591B (zh) | 一种发电机转子绕组动态匝间短路故障的定位方法 | |
El Sayed et al. | Fault diagnosis of PMSG stator inter-turn fault using extended Kalman filter and unscented Kalman filter | |
Zhang et al. | Shifting window average method for phasor measurement at offnominal frequencies | |
Zhao et al. | Probabilistic principal component analysis assisted new optimal scale morphological top-hat filter for the fault diagnosis of rolling bearing | |
CN107656203B (zh) | 一种永磁同步发电机的失磁故障检测方法 | |
CN106909076B (zh) | 一种pmsm永磁体局部退磁故障程度评估与数值化描述方法 | |
Liu et al. | A review of modeling and diagnostic techniques for eccentricity fault in electric machines | |
Zafarani et al. | Modeling and dynamic behavior analysis of magnet defect signatures in permanent magnet synchronous motors | |
Gyftakis et al. | Demagnetization monitoring and identification in PM generators with concentrated windings during transient conditions | |
Chang et al. | A synchrosqueezing transform-based hybrid method for voltage fluctuations assessment | |
CN106324505A (zh) | 一种电动机谐波分析方法 | |
CN113691183A (zh) | 永磁同步电机的电感辨识方法、系统、介质及终端 | |
Gyftakis et al. | Detection and identification of demagnetization in PM generators during transient conditions | |
Guan et al. | Adaptive surge detection of magnetic suspension centrifugal blower based on rotor radial displacement signal and SOGI-FLL with prefilter | |
Yang et al. | Fast nonlinear chirplet dictionary-based sparse decomposition for rotating machinery fault diagnosis under nonstationary conditions | |
Liu et al. | Diagnosis of air‐gap eccentricity and partial demagnetisation of an interior permanent magnet synchronous motor based on inverse transient complex inductance vector theory | |
DOĞAN et al. | A diagnosis of stator winding fault based on empirical mode decomposition in PMSMs | |
CN106610474B (zh) | 基于同步发电机异步旋转频率响应试验的参数测试方法 | |
Navarro-Navarro et al. | Application of transient analysis techniques to fault diagnosis in low-and medium-power synchronous machines |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20231216 Address after: 250000 Hengda Yueting, No. 58, Gongye North Road, Licheng District, Jinan City, Shandong Province Patentee after: Dongdai (Jinan) Intelligent Technology Co.,Ltd. Address before: 451191 Xianghe Road, Longhu Town, Zhengzhou City, Henan Province Patentee before: HENAN INSTITUTE OF ENGINEERING |