CN111510041A - 一种永磁同步电机的运行状态评估方法及*** - Google Patents
一种永磁同步电机的运行状态评估方法及*** Download PDFInfo
- Publication number
- CN111510041A CN111510041A CN202010528730.4A CN202010528730A CN111510041A CN 111510041 A CN111510041 A CN 111510041A CN 202010528730 A CN202010528730 A CN 202010528730A CN 111510041 A CN111510041 A CN 111510041A
- Authority
- CN
- China
- Prior art keywords
- permanent magnet
- equation
- stator
- determining
- pmsm
- 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
- 230000001360 synchronised effect Effects 0.000 title claims abstract description 118
- 238000011156 evaluation Methods 0.000 title claims abstract description 19
- 230000004907 flux Effects 0.000 claims abstract description 57
- 238000000034 method Methods 0.000 claims abstract description 34
- 230000005284 excitation Effects 0.000 claims abstract description 15
- 238000004804 winding Methods 0.000 claims description 39
- 230000009466 transformation Effects 0.000 claims description 26
- 238000006243 chemical reaction Methods 0.000 claims description 22
- 230000009467 reduction Effects 0.000 claims description 19
- 238000012545 processing Methods 0.000 claims description 17
- 238000013016 damping Methods 0.000 claims description 12
- 238000004422 calculation algorithm Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 239000011159 matrix material Substances 0.000 description 9
- 230000008569 process Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 230000036541 health Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000011426 transformation method Methods 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/14—Estimation or adaptation of machine parameters, e.g. flux, current or voltage
- H02P21/141—Flux estimation
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Ac Motors In General (AREA)
Abstract
本发明涉及一种永磁同步电机的运行状态评估方法及***。本发明利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链,基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程,对不可测参数永磁体磁链进行精确评估的同时,将不可测参数永磁体磁链融入到永磁同步电机的状态观测方程中,从而提高永磁同步电机的运行状态的评估精度。
Description
技术领域
本发明涉及工程技术领域,特别是涉及一种永磁同步电机的运行状态评估方法及***。
背景技术
为进一步对工业设备以及军事装备的可靠性和其他关键指标进行全面监控,大量的工程技术人员侧重于研究设备运行过程中变量的变化规律和***的运行状态。随着传感器技术的飞速发展,对设备运行过程中的大量参数的测量成为了可能。然而,针对一些特殊环境条件下或特定设备时,现有的传感器仍然不能满足对这些参数的测量。在***实际运行过程中,这些参数的变化会直接或间接影响到***正常工作和产品的最终质量,因此需要设计一种方法对这些参数进行辨识估计,作为评估***运行状态的一个重要依据。
由于不同设备的模型和工作原理不同,对其不可测参数的观测和辨识估计方法也不相同。实际上,各种对设备不可测参数观测和辨识估计方法的基本思想都大致相同,都是通过对可测的输入和输出之间的关系建立模型模拟未知参数在***中的变化规律,其本质上大多是黑箱建模。
以三相永磁同步电机(Permanent Magnet Synchronous Motor,PMSM)为例,电机运行在不同的工况时,磁饱和程度和温度等参数的变化会影响到定子电阻、电感和永磁体磁链变化,从而影响电机运行状态,为保证电机在运行中的平稳和可靠,需要对其各个参数的变化进行监测。此外,通过对永磁同步电机各个参数的观测还能够将这些观测数据与基于数据驱动的故障检测方法结合起来,对永磁同步电机进行健康状态监测和故障诊断等等。研究表明,永磁同步电机在运行过程中随着电机运行时间的增加,定子电阻大小会随着温度的升高变化较大,从而影响到电机的转子位置。针对该类具有不可测参数的电机设备,必须对其搭建状态观测器对其进行观测,首先利用辨识估计算法对其他不可测参数进行辨识估计,之后再将辨识估计得到的数据引入到设备模型的状态观测器中,利用观测得到的数据对电机进行性能分析,从而保证电机运行的安全性和可靠性。但现有的PMSM模型无法对不可测参数永磁体磁链大小进行精确辨识评估,且辨识估计的数据无法融入到永磁同步电机的状态观测方程中,导致永磁同步电机的运行状态评估精度低。
发明内容
本发明的目的是提供一种永磁同步电机的运行状态评估方法及***,以解决现有的PMSM模型评估永磁同步电机的运行状态时,评估精度低的问题。
为实现上述目的,本发明提供了如下方案:
一种永磁同步电机的运行状态评估方法,包括:
获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程;
获取三相永磁同步电机的电机参数;所述电机参数包括电机的机械角速度、电磁转矩、阻尼系数以及转动惯量;
基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程;
对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型;所述不可测参数包括转子永磁体产生的励磁磁场以及三相绕组的定子电阻;所述不可测参数影响下的PMSM模型包括同步旋转坐标系下的PMSM模型以及同步旋转坐标系下的机械运动方程;
利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链;
基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程;
获取所述三相永磁同步电机的实际电压以及实际电流;
以所述实际电压为所述离散化后的定子电阻状态观测方程的输入,以定子电阻为输出,根据所述定子电阻状态观测方程确定定子电阻观测值;
根据所述实际电压以及所述定子电阻确定观测电流;
将所述实际电流以及所述观测电流的电流误差作为补偿,输入到所述定子电阻状态观测方程中,确定下一个定子电阻观测值,直至所述电流误差为0;
根据所述定子电阻观测值评估所述永磁同步电机的运行状态;所述运行状态包括健康状态以及故障状态。
可选的,所述获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程,具体包括:
可选的,所述基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程,具体包括:
可选的,所述对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型,具体包括:
根据公式对所述三相电压方程进行降阶和解耦变换处理,确定同步旋转坐标系下的PMSM模型;其中,ud为定子电压d轴分量;uq为定子电压q轴分量;Rs是PMSM的定子电阻;id为定子电流d轴分量;iq为定子电流q轴分量;ψd为定子磁链d轴分量;ψq是定子磁链q轴分量;wr是电角速度,p为微分算子;
可选的,所述基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程,具体包括:
为实现上述目的,本发明还提供了如下方案:
一种永磁同步电机的运行状态评估***,包括:
三相电压方程获取模块,用于获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程;
电机参数获取模块,用于获取三相永磁同步电机的电机参数;所述电机参数包括电机的机械角速度、电磁转矩、阻尼系数以及转动惯量;
机械运动方程确定模块,用于基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程;
不可测参数影响下的PMSM模型确定模块,用于对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型;所述不可测参数包括转子永磁体产生的励磁磁场以及三相绕组的定子电阻;所述不可测参数影响下的PMSM模型包括同步旋转坐标系下的PMSM模型以及同步旋转坐标系下的机械运动方程;
永磁体磁链辨识估计模块,用于利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链;
定子电阻状态观测方程确定模块,用于基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程;
实际电压以及实际电流获取模块,用于获取所述三相永磁同步电机的实际电压以及实际电流;
定子电阻观测值确定模块,用于以所述实际电压为所述离散化后的定子电阻状态观测方程的输入,以定子电阻为输出,根据所述定子电阻状态观测方程确定定子电阻观测值;
观测电流确定模块,用于根据所述实际电压以及所述定子电阻确定观测电流;
补偿模块,用于将所述实际电流以及所述观测电流的电流误差作为补偿,输入到所述定子电阻状态观测方程中,确定下一个定子电阻观测值,直至所述电流误差为0;
运行状态评估模块,用于根据所述定子电阻观测值评估所述永磁同步电机的运行状态;所述运行状态包括健康状态以及故障状态。
可选的,所述三相电压方程获取模块具体包括:
三相电压方程获取单元,用于根据公式对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程;其中,us为三相绕组的定子电压,Rs为三相绕组的定子电阻,is为三相绕组的定子电流,ψs为电机的定子磁链。
可选的,所述机械运动方程确定模块具体包括:
可选的,所述不可测参数影响下的PMSM模型确定模块具体包括:
同步旋转坐标系下的PMSM模型确定单元,用于根据公式对所述三相电压方程进行降阶和解耦变换处理,确定同步旋转坐标系下的PMSM模型;其中,ud为定子电压d轴分量;uq为定子电压q轴分量;Rs是PMSM的定子电阻;id为定子电流d轴分量;iq为定子电流q轴分量;ψd为定子磁链d轴分量;ψq是定子磁链q轴分量;wr是电角速度,p为微分算子;
可选的,所述定子电阻状态观测方程确定模块具体包括:
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供了一种永磁同步电机的运行状态评估方法及***,本发明利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链,基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程,对不可测参数永磁体磁链进行精确评估的同时,将不可测参数永磁体磁链融入到永磁同步电机的状态观测方程中,从而提高永磁同步电机的运行状态的评估精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明所提供的永磁同步电机的运行状态评估方法流程图;
图2为本发明所提供的三相两极式永磁同步电机物理模型示意图;
图3为本发明所提供的各个坐标系之间的关系图;
图4为本发明所提供的Park变换矩阵仿真模型示意图;
图5为本发明所提供的基于滑模观测器的PMSM定子电阻观测器示意图;
图6为本发明设计的PMSM定子电阻状态观测器结构框图;
图7为本发明所提供的永磁同步电机的运行状态评估***结构图;
图8为本发明所提供的PMSM的三相电压曲线图;
图9为本发明所提供的PMSM的三相电流曲线图;
图10为本发明所提供的PMSM的d-q轴电压曲线图;
图11为本发明所提供的PMSM的转速曲线图;
图12为本发明所提供的PMSM的磁链估计曲线图;
图13为本发明所提供的PMSM的定子电阻观测值曲线图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种永磁同步电机的运行状态评估方法及***,能够提高永磁同步电机的运行状态的评估精度。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的永磁同步电机的运行状态评估方法流程图,如图1所示,一种永磁同步电机的运行状态评估方法,包括:
步骤101:获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程。
三相永磁同步电机是一个强耦合、复杂的非线性***,图2为本发明所提供的三相两极式永磁同步电机物理模型示意图,为了简化电机模型,在下列假设的前提下对PMSM模型简化后进行建模:
(1)假设PMSM忽略铁芯饱和。
(2)不计涡流与磁滞损耗。
(3)永磁体的电导率为零和相绕组中感应电动势的波形为正弦波形。
简化后自然坐标系下PMSM的三相电压方程为:
其中,us为三相绕组的定子电压,Rs为三相绕组的定子电阻,is为三相绕组的定子电流,ψs为电机的定子磁链,并具有如下形式:
其中,Ls为PMSM的定子电感,ψf为转子永磁体产生的励磁磁场,其在静止坐标系下的相位取决于电角度θr,上述各参数满足:
其中,iA、iB、iC分别为PMSM的三相电流,R为PMSM的三相绕组的电阻,ψA、ψB、ψC分别为三相绕组的磁链。
其中,uA、uB、uC分别为PMSM的三相电压。
其中,Lm3为定子互感;Ll3为定子漏感。根据公式(1)和公式(2)可得如下关系:
步骤102:获取三相永磁同步电机的电机参数;所述电机参数包括电机的机械角速度、电磁转矩、阻尼系数以及转动惯量。
步骤103:基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程。
根据机电能量转换关系,可得电机的机械运动方程为
其中,wm为电机的机械角速度,Te为电磁转矩,TL为负载转矩,B为阻尼系数,J为转动惯量。公式(7)即为PMSM在自然坐标系下的基本数学模型。可以看出三相PMSM的模型是一个较为复杂的多变量***,为了便于设计电机控制器需要对该模型进行降阶和解耦变换。常用的坐标变换方法有静止坐标变换(Clark变换)和同步旋转坐标变换(Park变换),本发明选用Park变换对电机模型进行处理。
步骤104:对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型;所述不可测参数包括转子永磁体产生的励磁磁场以及三相绕组的定子电阻;所述不可测参数影响下的PMSM模型包括同步旋转坐标系下的PMSM模型以及同步旋转坐标系下的机械运动方程。
将自然坐标系ABC变换到静止坐标系α-β的坐标变换成为Clark变换,将静止坐标系α-β变换到同步旋转坐标系d-q的坐标变换成为Park变换,图3为本发明所提供的各个坐标系之间的关系图。首先将自然坐标系下的变量转换为静止坐标系下的变量,可以使用以下坐标变换公式:
[fαfβf0]Τ=T3s/2s[fAfBfC]Τ (8)
其中,公式(8)的左侧为静止坐标系下的各变量,公式(8)的右侧则为自然坐标系下的各变量。f代表电机的电压、电流或磁链等变量,T3s/2s为坐标变换矩阵,具体表示方式如下:
将静止坐标系α-β变换到同步旋转坐标系d-q可表示为:
[fdfq]Τ=T2s/2r[fαfβ]Τ (10)
其中,公式(10)的左侧为同步旋转坐标系下的各变量,公式(10)的右侧则为静止坐标系下的各变量。T2s/2r为坐标变换矩阵,可表示为
由公式(8)-(11)可得自然坐标系ABC变换到同步旋转坐标系d-q的坐标变换关系如下,图4为本发明所提供的Park变换矩阵仿真模型示意图:
[fd fq f0]Τ=T3s/2r[fA fB fC]Τ (12)
其中,T3s/2r为坐标变换矩阵,可表示为
需要特别说明的是对于三相对称***而言,在计算时零序分量f0可以忽略不计。由此,得到了同步旋转坐标系下的PMSM模型:
其中,ud、uq分别是定子电压d-q轴分量,Rs是PMSM的定子电阻,id、iq分别是定子电流d-q轴分量,ψd、ψq是定子磁链的d-q轴分量,wr是电角速度,p为微分算子,且有
Te=pn(ψdiq-ψqid) (16)
其中,pn表示电机极对数。至此,本发明得到了在不可测参数Rs和ψf影响下的PMSM模型描述。
步骤105:利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链。
步骤106:基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程。
通过对上述考虑不可测参数下的PMSM模型描述,结合式(14)和式(15)推导关系可以得到如下的PMSM电压方程式:
方程(18)中的各个变量均为瞬时值,ud、uq、id、iq分别为定子绕组的d、q轴分量的电压和电流,而Rs是定子电阻,Rs为不可测且为本发明需要对其观测的参数,Ld和Lq分别是定子绕组的d、q轴的电感,wr为电角速度,p为微分算子。由公式(18)可以得到如下关系:
将式(19)进一步化简可以得到:
为了能够计算,还应该对定子相电阻的状态观测方程进行离散化处理,取足够小的采样时间Ts,那么会有如下公式:
结合式(18)和式(22)、(23)可得永磁同步电机离散的动态方程如下式:
iq(k)=aiq(k-1)+b[uq(k)+uq(k-1)]+c[wr(k)id(k)+wr(k-1)id(k-1)]+d[wr(k)+wr(k-1)]
(24)
为了能够对定子电阻进行观测,掌握电机运行过程中的健康状态,需要首先对其他的不可测参数利用最小二乘及其扩展算法进行辨识估计。本发明采用多新息最小二乘算法对PMSM的永磁体磁链ψf进行辨识估计,多新息最小二乘算法的基本思想如下式所示:
其中p为多新息,对算法进行初值设定,一般初值设定为如下关系:
P(0)=αI (27)
α为充分大的正实数(104~106);ε为零向量。其中:
Y(p,k)=[y(k)y(k-1)…y(k-p+1)](32)
其中,λ为遗忘因子,P(k)为协方差矩阵。在本发明中PMSM的输入矩阵和输出矩阵分别具有如下形式:
Y(p,k)=[iq(k-1)iq(k-1)…iq(k-p+1)]Τ (34)
采用多新息最小二乘算法对永磁体磁链辨识估计的对应算法为:
最后,得出多新息最小二乘算法辨识估计后的永磁体磁链如下:
步骤107:获取所述三相永磁同步电机的实际电压以及实际电流。
步骤108:以所述实际电压为所述离散化后的定子电阻状态观测方程的输入,以定子电阻为输出,根据所述定子电阻状态观测方程确定定子电阻观测值。
步骤109:根据所述实际电压以及所述定子电阻确定观测电流。
步骤110:将所述实际电流以及所述观测电流的电流误差作为补偿,输入到所述定子电阻状态观测方程中,确定下一个定子电阻观测值,直至所述电流误差为0。
步骤111:根据所述定子电阻观测值评估所述永磁同步电机的运行状态;所述运行状态包括健康状态以及故障状态。
本发明通过使用RLS算法对永磁体磁链进行辨识估计后,将得到的PMSM永磁体磁链估计数据加入到状态观测器中,用来获得更准确的观测数据。为实现对PMSM不可测参数的观测,减小不可测参数对***运行过程中产生的影响,本发明设计了基于滑模观测器的PMSM定子电阻观测器,其输入输出关系如图5所示。
本发明的Vs作为输入驱动电机运行,同时将电压Vs作为观测器的输入,根据状态观测方程得出定子电阻观测值然后将定子电阻作为输入,根据PMSM定子电压、电流和电阻三者之间的关系方程计算出观测电流大小;再通过将实际的电流值与观测电流值做差后得到电流误差将输入到滑模控制器内;其中,k为常数,sgn(x)为符号函数:
图7为本发明所提供的永磁同步电机的运行状态评估***结构图,如图7所示,一种永磁同步电机的运行状态评估***,包括:
三相电压方程获取模块701,用于获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程。
所述三相电压方程获取模块701具体包括:根据公式对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程;其中,us为三相绕组的定子电压,Rs为三相绕组的定子电阻,is为三相绕组的定子电流,ψs为电机的定子磁链。
电机参数获取模块702,用于获取三相永磁同步电机的电机参数;所述电机参数包括电机的机械角速度、电磁转矩、阻尼系数以及转动惯量。
机械运动方程确定模块703,用于基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程。
所述机械运动方程确定模块703具体包括:基于机电能量转换关系,根据公式确定所述三相永磁同步电机的机械运动方程;其中,J为转动惯量;Te为电磁转矩;TL为负载转矩;B为阻尼系数;wm为电机的机械角速度。
不可测参数影响下的PMSM模型确定模块704,用于对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型;所述不可测参数包括转子永磁体产生的励磁磁场以及三相绕组的定子电阻;所述不可测参数影响下的PMSM模型包括同步旋转坐标系下的PMSM模型以及同步旋转坐标系下的机械运动方程。
所述不可测参数影响下的PMSM模型确定模块704具体包括:同步旋转坐标系下的PMSM模型确定单元,用于根据公式对所述三相电压方程进行降阶和解耦变换处理,确定同步旋转坐标系下的PMSM模型;其中,ud为定子电压d轴分量;uq为定子电压q轴分量;Rs是PMSM的定子电阻;id为定子电流d轴分量;iq为定子电流q轴分量;ψd为定子磁链d轴分量;ψq是定子磁链q轴分量;wr是电角速度,p为微分算子;同步旋转坐标系下的机械运动方程确定单元,用于基于机电能量转换关系,根据公式对所述机械运动方程进行降阶和解耦变换处理,得到或确定同步旋转坐标系下的机械运动方程。
永磁体磁链辨识估计模块705,用于利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链;
定子电阻状态观测方程确定模块706,用于基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程。
实际电压以及实际电流获取模块707,用于获取所述三相永磁同步电机的实际电压以及实际电流。
定子电阻观测值确定模块708,用于以所述实际电压为所述离散化后的定子电阻状态观测方程的输入,以定子电阻为输出,根据所述定子电阻状态观测方程确定定子电阻观测值。
所述定子电阻状态观测方程确定模块708具体包括:定子电阻状态观测方程确定单元,用于根据公式确定定子电阻状态观测方程;其中,为定子电阻观测值;Ld为定子绕组的d轴电感;Lq为定子绕组的q轴电感;ψf为辨识估计后的永磁体磁链。
观测电流确定模块709,用于根据所述实际电压以及所述定子电阻确定观测电流。
补偿模块710,用于将所述实际电流以及所述观测电流的电流误差作为补偿,输入到所述定子电阻状态观测方程中,确定下一个定子电阻观测值,直至所述电流误差为0。
运行状态评估模块711,用于根据所述定子电阻观测值评估所述永磁同步电机的运行状态;所述运行状态包括健康状态以及故障状态。
下面以一个具体的实际应用来阐述本发明所要保护的技术方案:
三相永磁同步电机被广泛应用于工业生产的各个领域,发挥了巨大的作用。永磁同步电机作为电驱动伺服机构中的关键组成部件,其性能会直接影响到整个***的运作,因此,保证永磁同步电机在运行中的安全性和稳定性变得至关重要。为了对电机的运行状态进行全方位的监控,需要对电机的各个参数的变化进行采样,然而,由于永磁同步电机的结构和安装以及环境限制等原因,无法对电机的所有变量都使用传感器进行数据采集。所以,对于一些无法采集的参数需要进行观测、辨识和估计,将这些不可测参数进行观测、辨识估计后得到的数据与传感器采集得到的数据融合可以建立更具代表性的模型,大大提高了永磁同步电机的监控性能。
本发明基于三相永磁同步电机数据进行实验,验证本发明所提供的不可测参数观测、辨识估计方法。本发明的每组变量都包含1000个采样数据,采样时间间隔为Ts=1×10-5,电机选用380V的三相电压,图8和图9分别为PMSM的三相电压曲线图、电流曲线图。
图10为本发明所提供的PMSM的d-q轴电压曲线图,PMSM的d-q轴电压由三相电压经过坐标变换所得;图11为本发明所提供的PMSM的转速曲线图,本发明实验将电机的转速控制在600r/min,负载转矩初始值设置为0。
图12为本发明所提供的PMSM的磁链估计曲线图,本发明使用的永磁同步电机磁链为ψf=0.3Wb,通过递推最小二乘算法估计后的值收敛于0.3Wb,可以看出利用本发明的RLS辨识估计算法得出的辨识结果有效。
图13为本发明所提供的PMSM的定子电阻观测值曲线图,本发明使用的永磁同步电机定子电阻值为R=2.875Ω,可以看出利用状态观测器所观测出的电阻值有所增大,但是比较接近电阻的实际值,表明本发明所搭建的状态观测器有效。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种永磁同步电机的运行状态评估方法,其特征在于,包括:
获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程;
获取三相永磁同步电机的电机参数;所述电机参数包括电机的机械角速度、电磁转矩、阻尼系数以及转动惯量;
基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程;
对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型;所述不可测参数包括转子永磁体产生的励磁磁场以及三相绕组的定子电阻;所述不可测参数影响下的PMSM模型包括同步旋转坐标系下的PMSM模型以及同步旋转坐标系下的机械运动方程;
利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链;
基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程;
获取所述三相永磁同步电机的实际电压以及实际电流;
以所述实际电压为所述离散化后的定子电阻状态观测方程的输入,以定子电阻为输出,根据所述定子电阻状态观测方程确定定子电阻观测值;
根据所述实际电压以及所述定子电阻确定观测电流;
将所述实际电流以及所述观测电流的电流误差作为补偿,输入到所述定子电阻状态观测方程中,确定下一个定子电阻观测值,直至所述电流误差为0;
根据所述定子电阻观测值评估所述永磁同步电机的运行状态;所述运行状态包括健康状态以及故障状态。
6.一种永磁同步电机的运行状态评估***,其特征在于,包括:
三相电压方程获取模块,用于获取三相永磁同步电机PMSM模型,并对所述PMSM模型进行简化,确定自然坐标系下PMSM的三相电压方程;
电机参数获取模块,用于获取三相永磁同步电机的电机参数;所述电机参数包括电机的机械角速度、电磁转矩、阻尼系数以及转动惯量;
机械运动方程确定模块,用于基于机电能量转换关系,根据所述电机参数确定所述三相永磁同步电机的机械运动方程;
不可测参数影响下的PMSM模型确定模块,用于对所述三相电压方程以及所述机械运动方程进行降阶和解耦变换处理,确定不可测参数影响下的PMSM模型;所述不可测参数包括转子永磁体产生的励磁磁场以及三相绕组的定子电阻;所述不可测参数影响下的PMSM模型包括同步旋转坐标系下的PMSM模型以及同步旋转坐标系下的机械运动方程;
永磁体磁链辨识估计模块,用于利用多新息最小二乘法对所述转子永磁体产生的励磁磁场进行辨识估计,确定辨识估计后的永磁体磁链;
定子电阻状态观测方程确定模块,用于基于所述辨识估计后的永磁体磁链,根据所述不可测参数影响下的PMSM模型确定定子电阻状态观测方程;
实际电压以及实际电流获取模块,用于获取所述三相永磁同步电机的实际电压以及实际电流;
定子电阻观测值确定模块,用于以所述实际电压为所述离散化后的定子电阻状态观测方程的输入,以定子电阻为输出,根据所述定子电阻状态观测方程确定定子电阻观测值;
观测电流确定模块,用于根据所述实际电压以及所述定子电阻确定观测电流;
补偿模块,用于将所述实际电流以及所述观测电流的电流误差作为补偿,输入到所述定子电阻状态观测方程中,确定下一个定子电阻观测值,直至所述电流误差为0;
运行状态评估模块,用于根据所述定子电阻观测值评估所述永磁同步电机的运行状态;所述运行状态包括健康状态以及故障状态。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010528730.4A CN111510041B (zh) | 2020-06-11 | 2020-06-11 | 一种永磁同步电机的运行状态评估方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010528730.4A CN111510041B (zh) | 2020-06-11 | 2020-06-11 | 一种永磁同步电机的运行状态评估方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111510041A true CN111510041A (zh) | 2020-08-07 |
CN111510041B CN111510041B (zh) | 2021-07-20 |
Family
ID=71877120
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010528730.4A Active CN111510041B (zh) | 2020-06-11 | 2020-06-11 | 一种永磁同步电机的运行状态评估方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111510041B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110286332A (zh) * | 2019-06-17 | 2019-09-27 | 南通大学 | 一种基于多新息理论的电动汽车动力电池soc估计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150194914A1 (en) * | 2013-03-15 | 2015-07-09 | Texas Instruments Incorporated | Automated motor control |
CN107425775A (zh) * | 2017-07-04 | 2017-12-01 | 武汉理工大学 | 一种基于改进最小二乘法的永磁同步电机参数辨识*** |
US20180076749A1 (en) * | 2013-10-03 | 2018-03-15 | Texas Instruments Incorporated | Stator Resistance Estimation for Electric Motors |
CN109728758A (zh) * | 2018-08-08 | 2019-05-07 | 济南大学 | 一种永磁同步电机三相定子电阻的辨识方法 |
CN110176895A (zh) * | 2019-05-28 | 2019-08-27 | 上海联影医疗科技有限公司 | 一种电机参数辨识方法及装置、电机转速调节装置 |
CN111064406A (zh) * | 2019-12-31 | 2020-04-24 | 上海应用技术大学 | 基于限定记忆最小二乘法的改进模型参考自适应控制*** |
-
2020
- 2020-06-11 CN CN202010528730.4A patent/CN111510041B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150194914A1 (en) * | 2013-03-15 | 2015-07-09 | Texas Instruments Incorporated | Automated motor control |
US20180076749A1 (en) * | 2013-10-03 | 2018-03-15 | Texas Instruments Incorporated | Stator Resistance Estimation for Electric Motors |
CN107425775A (zh) * | 2017-07-04 | 2017-12-01 | 武汉理工大学 | 一种基于改进最小二乘法的永磁同步电机参数辨识*** |
CN109728758A (zh) * | 2018-08-08 | 2019-05-07 | 济南大学 | 一种永磁同步电机三相定子电阻的辨识方法 |
CN110176895A (zh) * | 2019-05-28 | 2019-08-27 | 上海联影医疗科技有限公司 | 一种电机参数辨识方法及装置、电机转速调节装置 |
CN111064406A (zh) * | 2019-12-31 | 2020-04-24 | 上海应用技术大学 | 基于限定记忆最小二乘法的改进模型参考自适应控制*** |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110286332A (zh) * | 2019-06-17 | 2019-09-27 | 南通大学 | 一种基于多新息理论的电动汽车动力电池soc估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111510041B (zh) | 2021-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rafaq et al. | A comprehensive review of state-of-the-art parameter estimation techniques for permanent magnet synchronous motors in wide speed range | |
Liu et al. | Global identification of electrical and mechanical parameters in PMSM drive based on dynamic self-learning PSO | |
Jones et al. | A state observer for the permanent-magnet synchronous motor | |
Akin et al. | Simple derivative-free nonlinear state observer for sensorless AC drives | |
CN114006557B (zh) | 基于扩展滑模观测器的永磁同步电机机械参数辨识方法 | |
CN113131817A (zh) | 一种永磁同步电机在线参数辨识***及方法 | |
Ezzat et al. | Sensorless speed control of permanent magnet synchronous motor by using sliding mode observer | |
CN103259479A (zh) | 一种永磁同步电机神经网络左逆状态观测方法 | |
CN104779873B (zh) | 一种用于pmsm伺服***的预测函数控制方法 | |
CN112149274B (zh) | 一种带有死区输入非线性的多轴雕刻机***在线建模方法 | |
Liu et al. | Extended-Kalman-filter-based magnet flux linkage and inductance estimation for PMSM considering magnetic saturation | |
CN111510041B (zh) | 一种永磁同步电机的运行状态评估方法及*** | |
Jannati et al. | Speed sensorless fault-tolerant drive system of 3-phase induction motor using switching extended kalman filter | |
CN109687792B (zh) | 面向矢量控制***的牵引电机转子参数在线辨识优化方法 | |
Khov et al. | Monitoring of turn short-circuit faults in stator of PMSM in closed loop by on-line parameter estimation | |
CN115664283A (zh) | 一种基于广义参数估计观测器的滑模控制方法及*** | |
Chen et al. | [Retracted] Nonlinear Identification of PMSM Rotor Magnetic Linkages Based on an Improved Extended Kalman Filter | |
Marwali et al. | Implementation of indirect vector control on an integrated digital signal processor-based system | |
Hezzi et al. | Sensorless control strategy for five-phase pmsm based on uio for linear parameter varying system | |
Vesely et al. | Stator resistance identification of PMSM | |
Ramana et al. | Parameter estimation of permanent magnet synchronous motor-a review | |
CN116317791B (zh) | 一种同步磁阻电机无位置传感器电感辨识方法及装置 | |
Ghanes et al. | Novel controller for induction motor without mechanical sensor and experimental validation | |
Huazhong | PMSM parameter identification motor test system framework based on virtual instrument technology | |
Cai et al. | Offline Identification of Induction Motor Rotor Time Constant Based on Least Squares Algorithm |
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 |