[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN103904973A - 一种实现凸极永磁同步电机无速度传感器控制的方法 - Google Patents

一种实现凸极永磁同步电机无速度传感器控制的方法 Download PDF

Info

Publication number
CN103904973A
CN103904973A CN201410110944.4A CN201410110944A CN103904973A CN 103904973 A CN103904973 A CN 103904973A CN 201410110944 A CN201410110944 A CN 201410110944A CN 103904973 A CN103904973 A CN 103904973A
Authority
CN
China
Prior art keywords
msub
mrow
mtd
mtr
mover
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
Application number
CN201410110944.4A
Other languages
English (en)
Other versions
CN103904973B (zh
Inventor
张兴
郭磊磊
杨淑英
谢震
曹朋朋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Huangshan Development Investment Group Co ltd
Original Assignee
Hefei University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN201410110944.4A priority Critical patent/CN103904973B/zh
Publication of CN103904973A publication Critical patent/CN103904973A/zh
Application granted granted Critical
Publication of CN103904973B publication Critical patent/CN103904973B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Control Of Ac Motors In General (AREA)
  • Control Of Motors That Do Not Use Commutators (AREA)

Abstract

本发明公开了一种实现凸极永磁同步电机无速度传感器控制方法,属于电机调速领域。其特征在于,根据实时测得的定子电压和定子电流,通过建立的全阶状态滑模观测器估计得到有效反电动势,再通过锁相环获得转子位置和转速。本发明解决了凸极永磁同步电机无速度传感器控制的问题,使无速度传感器控制受电机参数影响更小,且在实现转子位置和转速同时估计时,提高了转子位置估计对转速估计偏差的鲁棒性,也避免了放大滑模噪声,保证了系统的稳定性。

Description

一种实现凸极永磁同步电机无速度传感器控制的方法
技术领域
本发明涉及一种电机无速度传感器控制的方法,尤其是用于风力发电场合的一种实现凸极永磁同步电机无速度传感器控制的方法。
背景技术
为了应对能源危机,新能源技术得到了蓬勃的发展。风能是一种取之不尽、用之不竭的可再生能源,风力发电技术已成为新能源技术的重中之重。凸极永磁同步电机因为高功率密度、高可靠性等优点而被广泛地应用于风力发电系统。为提高系统可靠性,风力发电系统常常要求实现无速度传感器控制。常用的无速度传感器控制方法有很多,如高频信号注入法、模型参考自适应法和滑模观测器法等。
发明专利申请《一种内置式永磁同步电机无位置传感器矢量控制装置》(CN102361430A)和《无位置传感器电机控制系统及控制方法》(CN102624322A)都提出了基于高频信号注入的永磁同步电机无速度传感器控制方法,但高频信号注入法只适用于低速区域,不适用于大功率风力发电系统。发明专利申请《一种永磁同步电机无位置传感器控制方法》(CN103051271A)设计了基于模型参考自适应原理的永磁同步电机无速度传感器控制方法,但受电机参数的影响较严重。发明专利申请《永磁同步电机的无传感器控制系统》(CN101964624A)提出了基于二阶滑模观测器的永磁同步电机无速度传感器控制算法,然而,由于估计的反电动势含有高频滑模噪声,必须进行滤波,可是滤波却带来了相位延迟,还必须进行相位补偿,导致系统实现复杂,角度估计受转速估计值影响严重。
2012年哈尔滨工业大学王高林等发表的会议文献《Rotor PositionEstimation with Full-Order Sliding-Mode Observer for SensorlesssIPMSM》提出了基于扩展反电动势的凸极永磁同步电机无速度传感器控制方法,解决了二阶滑模控制存在的问题,但是基于扩展反电动势的凸极永磁同步电机模型受电机定子电阻、d轴电感和q轴电感三个参数的影响。2012年Mihai Comanescu等学者发表的论文《A Family of Sensorless Observerswith Speed Estimate for Rotor Position Estimation of IM and PMSMDrives》设计了一个针对面装式永磁同步电机的全阶滑模观测器,研究了转子位置估计对转速估计误差的鲁棒性,但是为了提高其鲁棒性,需要增大滑模增益,可是增大滑模增益会放大滑模噪声,严重时会导致系统不稳定,且该方法不能直接应用于凸极永磁同步电机中。
发明内容
本发明要解决的技术问题为克服上述各种技术方案的局限性,提供一种将状态观测器和滑模观测器融为一体的一种实现凸极永磁同步电机无速度传感器控制的方法。
为解决本发明的技术问题,所采用的技术方案为:一种实现凸极永磁同步电机无速度传感器控制的方法包括定子电压uAB、uBC、uCA和定子电流iA、iB、iC的测量,特别是,
根据实时测得的定子电压uAB、uBC、uCA和定子电流iA、iB、iC,通过全阶状态滑模观测器得到估计的有效反电动势
Figure BDA0000480886130000021
估计的有效反电动势
Figure BDA0000480886130000022
再通过锁相环获得转子位置和转速;
所述全阶状态滑模观测器为:
p i ^ α = - R s L q i ^ α + u α L q - e ^ α L q - 1 L q s α p i ^ β = - R s L q i ^ β + u β L q - e ^ β L q - 1 L q s β p e ^ α = - ω ^ r e ^ α + Ns β p e ^ β = ω ^ r e ^ α + Ns β ,
其中, s α = Msgn ( i ^ α - i α ) + K ( i ^ α - i α ) , s β = Msgn ( i ^ β - i β ) + K ( i ^ β - i β ) , M为滑模增益,sgn()为符号函数,K为状态增益,N为有效反电动势eα、eβ滑模观测器的增益,
Figure BDA0000480886130000025
为在静止α-β坐标系上的估计电流值,
Figure BDA0000480886130000026
为估计的有效反电动势,
Figure BDA0000480886130000027
为估计的同步转速,Lq为电机q轴电感,Rs为电机定子电阻。
作为一种实现凸极永磁同步电机无速度传感器控制的方法的进一步改进:
优选地,得到估计的有效反电动势的步骤如下:
步骤1,如电机定子绕组为三角形接法,则执行步骤1.1,如电机定子绕组为星形接法,则执行步骤1.2,
步骤1.1,先采样电机定子电压uAB、uBC、uCA,再将uAB、uBC、uCA坐标变换到静止α-β坐标系上,如下式:
u α u β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u AB u BC u CA ,
其中,uα、uβ为定子电压在静止α-β坐标系上的分量,之后,先采样电机定子电流iA、iB,再计算相电流,如下式:
i Ax = 1 / 3 ( i A - i B ) i Bx = i Ax + i B i Cx = i Ax - i A
其中,iAx、iBx、iCx为三相相电流,最后,将三相相电流坐标变换到静止α-β坐标系上,如下式:
i α i β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i Ax i Bx i Cx ,
其中,iα、iβ为相电流在静止α-β坐标系上的分量;
步骤1.2,先采样电机定子电压uAB、uBC,再计算电机相电压,如下式:
u B = 1 / 3 ( u BC - u AB ) u A = u B + u AB u C = u B - u BC ,
其中,uA、uB、uC为三相相电压,之后,将三相相电压坐标变换到静止α-β坐标系上,如下式:
u α u β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u A u B u C ,
其中,uα、uβ为定子电压在静止α-β坐标系上的分量,最后,先采样电机定子电流iA、iB、iC,再经过坐标变换得到电流在静止α-β坐标系上的分量,如下式:
i α i β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i A i B i C ,
其中,iα、iβ为相电流在静止α-β坐标系上的分量;
步骤2,先用估计的电流
Figure BDA0000480886130000043
减去步骤1得到的实际电流iα、iβ,得到电流偏差
Figure BDA0000480886130000044
再分别经过符号函数sgn()运算和乘法运算,最后得到电流误差反馈量sα、sβ,如下式:
s α = Msgn ( i ‾ α ) + K ( i ‾ α ) , s β = Msgn ( i ‾ β ) + K ( i ‾ β ) ,
其中,M为滑模增益,K为状态增益,
步骤3,由步骤1得到的电压uα、uβ减去步骤2得到的电流误差反馈量sα、sβ,再减去第k次运算得到的估计的有效反电动势
Figure BDA0000480886130000047
得到第一组中间变量E、E,如下式:
E 1 α = u α - s α - e ^ α ( k ) E 1 β = u β - s β - e ^ β ( k ) ,
其中,E、E为第一组中间变量,k为正整数,(k=0,1,2········), e ^ β ( 0 ) = 0 ;
步骤4,用步骤3得到的第一组中间变量E、E除以电机q轴电感,得到第二组中间变量E、E,如下式:
E 2 α = E 1 α / L q E 2 β = E 1 β / L q ,
其中,Lq为电机q轴电感,E、E为第二组中间变量;
步骤5,根据估计的电机定子电流和电机定子电阻Rs以及q轴电感Lq,计算得到第三组中间变量E、E,如下式:
E 3 α = R s L q i ^ α E 3 β = R s L q i ^ β ,
其中,E、E为第三组中间变量,Rs为电机定子电阻;
步骤6,由步骤4得到的E、E减去步骤5得到的E、E,得到第四组中间变量E、E,如下式:
E 4 α = E 2 α - E 3 α E 4 β = E 2 β - E 3 β ,
其中,E、E为第四组中间变量;
步骤7,对步骤6得到的E、E进行积分,得到估计的电流
Figure BDA0000480886130000054
如下式:
i ^ α = i ^ α + T s E 4 α i ^ β = i ^ β + T s E 4 β ,
其中,Ts为采样周期;
步骤8,将步骤7得到的电流估计值
Figure BDA0000480886130000056
带入步骤2和步骤5,用于进一步估计定子电流;
步骤9,将步骤2得到的电流误差反馈量sα、sβ乘以增益N,得到第五组中间变量E、E,如下式:
E 5 α = Ns α E 5 β = Ns β ,
其中,E、E为第五组中间变量,N为有效反电动势eα、eβ滑模观测器的增益;
步骤10,估计的转速
Figure BDA0000480886130000058
分别乘以第k次运算得到的估计的有效反电动势
Figure BDA0000480886130000059
得到第六组中间变量E、E,如下式:
E 6 α = ω ^ r e ^ α ( k ) E 6 β = ω ^ r e ^ β ( k ) ,
其中,E、E为第六组中间变量;
步骤11,先由步骤9得到的E减去步骤10得到的E,得到第七个中间变量E,再由步骤9得到的E加上步骤10得到的E,得到第七个中间变量E,如下式:
E 7 α = E 5 α - E 6 β E 7 β = E 5 β + E 6 α ,
其中,E、E为第七组中间变量;
步骤12,分别对步骤11得到的E、E进行积分,得到第k+1次运算估计的有效反电动势
Figure BDA0000480886130000063
如下式:
e ^ α ( k + 1 ) = e ^ α ( k ) + T s E 7 α e ^ β ( k + 1 ) = e ^ β ( k ) + T s E 7 β ;
步骤13,先将步骤12得到的估计的有效反电动势
Figure BDA0000480886130000065
带入步骤3,进一步用于求解第一组中间变量E、E,再将步骤12得到的估计的有效反电动势
Figure BDA0000480886130000066
带入步骤10,用于计算第七组中间变量E、E
步骤14,将步骤12得到的估计的有效反电动势
Figure BDA0000480886130000067
带入锁相环,得到估计的转子位置
Figure BDA0000480886130000068
和估计的转速
Figure BDA0000480886130000069
步骤15,将估计的转速
Figure BDA00004808861300000610
带入步骤10,用于计算第七组中间变量E、E
步骤16,重复步骤1~步骤15,直至估计的电流值
Figure BDA00004808861300000611
等于实际的电流值iα、iβ,得到与实际的反电动势eα、eβ相符的估计的有效反电动势
Figure BDA00004808861300000612
和与实际的转子位置θ相符的估计的转子位置
Figure BDA00004808861300000613
优选地,锁相环根据估计的有效反电动势
Figure BDA00004808861300000614
获得转子位置
Figure BDA00004808861300000615
和估计的转速
Figure BDA00004808861300000616
的步骤如下:
步骤1,根据估计有效反电动势
Figure BDA00004808861300000617
计算有效反电动势q轴偏差,如下式:
e ^ q = - e ^ α ( k + 1 ) cos θ ^ - e ^ β ( k + 1 ) sin θ ^ ,
其中,
Figure BDA0000480886130000071
为有效反电动势q轴偏差;
步骤2,由步骤1得到的
Figure BDA0000480886130000072
经过比例积分调节器得到估计的转速
Figure BDA0000480886130000073
如下式:
ω ^ r = ( k i / s + k p ) e ^ q ,
其中,kp和ki分别为比例和积分系数,s为拉普拉斯算子;
步骤3,由步骤2得到的估计转速
Figure BDA0000480886130000075
通过积分得到估计的转子位置
Figure BDA0000480886130000076
如下式:
θ ^ = θ ^ + T s ω ^ r .
相对于现有技术的有益效果是:
本发明基于有效反电动势的概念设计了一个全阶状态滑模观测器。传统的基于有效磁链的凸极永磁同步电机全阶模型如下:
pi α = R s L q i α + u α L q - pψ α L q pi β = - R s L q i β + u β L q - pψ β L q pψ α = - ω r ψ β pψ β = ω r ψ α
ψ β ψ α = [ ψ f + ( L d - L q ) i d ] sin θ r cos θ r ;
其中,iα、iβ分别为定子电流在静止α-β坐标系上的两个分量,uα、uβ分别为定子电压在静止α-β坐标系上的两个分量,ψα、ψβ为有效磁链在静止α-β坐标系上的两个分量,ψf为永磁体磁链,Ld为d轴电感,Lq为q轴电感,Rs为定子电阻,ωr为同步转速,θr为转子位置。
本发明定义的有效反电动势eα、eβ如下:
e α = pψ α e β = pψ β .
认为磁链幅值不变可得:
pe α = - ω r e β pe β = ω r e α .
根据有效反电动势eα、eβ的定义,得到改进后的凸极永磁同步电机全阶模型如下:
pi α = - R s L q i α + u α L q - e α L q pi β = - R s L q i β + u β L q - e β L q pe α = - ω r e β pe β = ω r e α .
2012年Mihai Comanescu等学者发表的论文《A Family ofSensorless Observers with Speed Estimate for Rotor PositionEstimation of IM and PMSM Drives》设计了一个针对面装式永磁同步电机的全阶滑模观测器,其形式与上式一致。但是全阶滑模观测器还存在滑模抖振的问题。为此,基于本发明设计的改进型凸极永磁同步电机全阶模型,设计了全阶状态滑模观测器,引入状态误差项,可削弱滑模抖振,提高系统稳定性。
基于该模型建立的全阶状态滑模观测器如下:
p i ^ α = - R s L q i ^ α + u α L q - e ^ α L q - 1 L q s α p i ^ β = - R s L q i ^ β + u β L q - e ^ β L q - 1 L q s β p e ^ α = - ω ^ r e ^ β + Ns α p e ^ β = ω ^ r e ^ α + Ns β ;
其中, s α = Msgn ( i ^ α - i α ) + K ( i ^ α - i α ) , s β = Msgn ( i ^ β - i β ) + K ( i ^ β - i β ) , M为滑模增益,sgn()为符号函数,K为状态增益,N为有效反电动势eα、eβ滑模观测器的增益,
Figure BDA0000480886130000085
为在静止α-β坐标系上的估计电流值,
Figure BDA0000480886130000086
为估计的有效反电动势,
Figure BDA0000480886130000087
为估计的同步转速。
根据本发明建立的全阶状态滑模观测器,可实现有效反电动势
Figure BDA0000480886130000088
的估计,再通过锁相环可以得到转速和转子位置
Figure BDA0000480886130000092
本发明为了将全阶滑模观测器算法应用于凸极永磁同步电机中,首先改进了基于有效磁链的凸极永磁同步电机模型,提出了有效反电动势的概念;其次将状态观测器和滑模观测器融为一体,设计了状态滑模观测器,解决了提高转子位置估计对转速估计误差的鲁棒性时需要增大滑模增益,而抑制滑模噪声时却需要减小滑模增益这一对矛盾,实现了两者的统筹兼顾,从而提高了系统的稳定性和鲁棒性。经仿真实测,本发明可以自动实现转子位置的准确估计,实现稳定的无速度传感器控制。
本发明解决了凸极永磁同步电机无速度传感器控制的问题,使无速度传感器控制受电机参数影响更小,且在实现转子位置和转速同时估计时,提高了转子位置估计对转速估计偏差的鲁棒性,也避免了放大滑模噪声,保证了系统的稳定性。
附图说明
下面结合附图对本发明的优选方式作进一步详细的描述。
图1为电流状态滑模观测器的具体实施框图。
图2为有效反电动势状态滑模观测器的具体实施框图。
图3为锁相环具体实施图。
图4为M=100,N=50时反馈转速有偏差时的仿真波形图。
图5为M=100,N=500时反馈转速有偏差时的仿真波形图。
具体实施方式
参见图1、图2、图3、图4和图5,一种实现凸极永磁同步电机无速度传感器控制的方法如下:
首先进行定子电压uAB、uBC、uCA和定子电流iA、iB、iC的测量,再根据实时测得的定子电压uAB、uBC、uCA和定子电流iA、iB、iC,通过全阶状态滑模观测器得到估计的有效反电动势
Figure BDA0000480886130000093
估计的有效反电动势
Figure BDA0000480886130000094
再通过锁相环获得转子位置和转速;
所述全阶状态滑模观测器为:
p i ^ α = - R s L q i ^ α + u α L q - e ^ α L q - 1 L q s α p i ^ β = - R s L q i ^ β + u β L q - e ^ β L q - 1 L q s β p e ^ α = - ω ^ r e ^ β + Ns α p e ^ β = ω ^ r e ^ α + Ns β ,
其中, s α = Msgn ( i ^ α - i α ) + K ( i ^ α - i α ) , s β = Msgn ( i ^ β - i β ) + K ( i ^ β - i β ) , M为滑模增益,sgn()为符号函数,K为状态增益,N为有效反电动势eα、eβ滑模观测器的增益,
Figure BDA0000480886130000103
为在静止α-β坐标系上的估计电流值,
Figure BDA0000480886130000104
为估计的有效反电动势,
Figure BDA0000480886130000105
为估计的同步转速,Lq为电机q轴电感,Rs为电机定子电阻。
其中,本发明的具体步骤如下:
如图1所示,首先采样计算观测器所需的定子电压在静止α-β坐标系上的分量uα、uβ和相电流在静止α-β坐标系上的分量iα、iβ,再分两种情况:如果电机定子绕组是三角形接法,则执行步骤(a),如果电机定子绕组是星形接法,则执行步骤(b);
步骤(a),采样电机定子电压uAB、uBC、uCA,然后将uAB、uBC、uCA坐标变换到静止α-β坐标系上,如下:
u α u β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u AB u BC u CA ,
其中,uα、uβ为定子电压在静止α-β坐标系上的分量,
采样电机定子电流iA、iB,然后计算相电流,如下式:
i Ax = 1 / 3 ( i A - i B ) i Bx = i Ax + i B i Cx = i Ax - i A ,
其中,iAx、iBx、iCx为三相相电流,
然后将三相相电流坐标变换到静止α-β坐标系上,如下式:
i α i β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i Ax i Bx i Cx ,
其中,iα、iβ为相电流在静止α-β坐标系上的分量;
步骤(b),采样电机定子电压uAB、uBC,然后计算电机相电压,如下式:
u B = 1 / 3 ( u BC - u AB ) u A = u B + u AB u C = u B - u BC ,
其中,uA、uB、uC为三相相电压,
然后将三相相电压坐标变换到静止α-β坐标系上,如下式:
u α u β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u A u B u C ,
其中,uα、uβ为定子电压在静止α-β坐标系上的分量,
采样电机定子电流iA、iB、iC,然后经过坐标变换得到电流在静止α-β坐标系上的分量,如下式:
i α i β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i A i B i C ,
其中,iα、iβ为相电流在静止α-β坐标系上的分量;
在得到电压uα、uβ和电流iα、iβ后,如图1所示,再根据电压uα、uβ、第k次运算得到的估计的有效反电动势和电流误差反馈量sα、sβ,得到第一组中间变量E、E,如下式:
E 1 α = u α - s α - e ^ α ( k ) E 1 β = u β - s β - e ^ β ( k ) ,
其中,E、E为第一组中间变量,sα、sβ为电流误差反馈量,k为正整数,(k=0,1,2······),
Figure BDA0000480886130000121
sα、sβ的计算如下:
s α = Msgn ( i ^ α - i α ) + K ( i ^ α - i α ) , s β = Msgn ( i ^ β - i β ) + K ( i ^ β - i β ) ,
其中,M为滑模增益,K为状态增益,
Figure BDA0000480886130000123
估计的电流;
如图1所示,根据一组中间变量E、E计算得到第二组中间变量E、E,如下式:
E 2 α = E 1 α / L q E 2 β = E 1 β / L q ,
其中,Lq为电机q轴电感,E、E为第二组中间变量;
如图1所示,根据第二组中间变量E、E计算第四组中间变量E、E,如下式:
E 4 α = E 2 α - E 3 α E 4 β = E 2 β - E 3 β ,
其中,E、E为第三组中间变量,E、E为第四组中间变量;第三组中间变量E、E的计算如下式:
E 3 α = R s L q i ^ α E 3 β = R s L q i ^ β ,
其中,Rs为电机定子电阻;
如图1所示,对E、E进行积分,得到估计的电流
Figure BDA0000480886130000127
如下式:
i ^ α = i ^ α + T s E 4 α i ^ β = i ^ β + T s E 4 β ,
其中,Ts为采样周期;
如图2所示,电流误差反馈量sα、sβ乘以增益N得到第五组中间变量E、E,如下式:
E 5 α = Ns α E 5 β = Ns β ,
其中,E、E为第五组中间变量,N为有效反电动势eα、eβ滑模观测器的增益;
如图2所示,用估计的转速
Figure BDA0000480886130000131
分别乘以第k次运算得到的估计的有效反电动势得到第六组中间变量E、E,如下式:
E 6 α = ω ^ r e ^ α ( k ) E 6 β = ω ^ r e ^ β ( k ) ,
其中,E、E为第六组中间变量;
如图2所示,E减去E,得到第七个中间变量E,E加上E,得到第七个中间变量E,如下式:
E 7 α = E 5 α - E 6 β E 7 β = E 5 β + E 6 α ,
其中,E、E为第七组中间变量;
如图2所示,对E、E分别进行积分,得到第k+1次运算估计的有效反电动势 e ^ α ( k + 1 ) , e ^ β ( k + 1 ) , 如下式:
e ^ α ( k + 1 ) = e ^ α ( k ) + T s E 7 α e ^ β ( k + 1 ) = e ^ β ( k ) + T s E 7 β ;
如图3所示,根据第k+1次运算估计的有效反电动势
Figure BDA0000480886130000137
计算有效反电动势q轴偏差,如下式:
e ^ q = - e ^ α ( k + 1 ) cos θ ^ - e ^ β ( k + 1 ) sin θ ^ ,
其中,
Figure BDA0000480886130000139
为有效反电动势q轴偏差;
如图3所示,
Figure BDA00004808861300001310
经过比例积分调节器得到估计的转速
Figure BDA00004808861300001311
如下式:
ω ^ r = ( k i / s + k p ) e ^ q ,
其中,kp和ki分别为比例和积分系数,s为拉普拉斯算子;
如图3所示,估计的转速
Figure BDA00004808861300001313
通过积分得到估计的转子位置如下式:
θ ^ = θ ^ + T s ω ^ r .
为了验证本发明的有效性,进行了仿真验证。所用电机定子电阻为0.34欧,d轴电感为0.046H,q轴电感为0.135H,永磁体磁链为0.75Wb。采用本方法提出的全阶状态滑模观测器和锁相环实现转子位置估计。仿真时,给定励磁电流isd=0,isq=10A,给定电机转速为60Rad/s。首先从本方法建立的全阶状态滑模观测器来看,由于全阶状态滑模观测器不含d轴电感,故不受d轴电感的影响。为了进一步验证该方法转子位置估计对转速估计误差的鲁棒性,仿真时,在0.6s时人为在反馈速度上乘以一个误差系数0.8,即0.6s前采用估计的真实速度作为反馈,0.6s后采用错误的速度进行反馈,设计M=1,N=50,K=50时仿真结果如图4所示,可见当反馈速度出现偏差时,估计角度也出现了偏差,进一步增大N=500时,仿真结果如图5所示,可见,通过增大滑模增益N,可以提高转子位置估计对转速估计误差的鲁棒性。如果选择K=0,则本方法设计的全阶状态滑模观测器化简为全阶滑模观测器,此时为保证系统稳定,必须是当增大M,而为了保持滑模增益不变,即保持MN不变,则必须减小N,因此,转子位置估计对转速估计误差的鲁棒性也会降低。状态滑模观测器的优点正在于,通过引入状态误差反馈系数K,可允许减小滑模增益M,从而允许在一定的滑模噪声范围内,提高滑模增益N,以进一步提高转子位置估计对转速估计误差的鲁棒性。
显然,本领域的技术人员可以对本发明的一种实现凸极永磁同步电机无速度传感器控制的方法进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若对本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。

Claims (3)

1.一种实现凸极永磁同步电机无速度传感器控制的方法,包括定子电压uAB、uBC、uCA和定子电流iA、iB、iC的测量,其特征在于:
根据实时测得的定子电压uAB、uBC、uCA和定子电流iA、iB、iC,通过全阶状态滑模观测器得到估计的有效反电动势
Figure FDA0000480886120000011
估计的有效反电动势再通过锁相环获得转子位置和转速;
所述全阶状态滑模观测器为:
p i ^ α = - R s L q i ^ α + u α L q - e ^ α L q - 1 L q s α p i ^ β = - R s L q i ^ β + u β L q - e ^ β L q - 1 L q s β p e ^ α = - ω ^ r e ^ α + Ns β p e ^ β = ω ^ r e ^ α + Ns β ,
其中, s α = Msgn ( i ^ α - i α ) + K ( i ^ α - i α ) , s β = Msgn ( i ^ β - i β ) + K ( i ^ β - i β ) , M为滑模增益,sgn()为符号函数,K为状态增益,N为有效反电动势eα、eβ滑模观测器的增益,
Figure FDA0000480886120000015
为在静止α-β坐标系上的估计电流值,
Figure FDA0000480886120000016
为估计的有效反电动势,
Figure FDA0000480886120000017
为估计的同步转速,Lq为电机q轴电感,Rs为电机定子电阻。
2.根据权利要求1所述的一种实现凸极永磁同步电机无速度传感器控制的方法,其特征是得到估计的有效反电动势
Figure FDA0000480886120000018
的步骤如下:
步骤1,如电机定子绕组为三角形接法,则执行步骤1.1,如电机定子绕组为星形接法,则执行步骤1.2,
步骤1.1,先采样电机定子电压uAB、uBC、uCA,再将uAB、uBC、uCA坐标变换到静止α-β坐标系上,如下式:
u α u β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u AB u BC u CA ,
其中,uα、uβ为定子电压在静止α-β坐标系上的分量,之后,先采样电机定子电流iA、iB,再计算相电流,如下式:
i Ax = 1 / 3 ( i A - i B ) i Bx = i Ax + i B i Cx = i Ax - i A ,
其中,iAx、iBx、iCx为三相相电流,最后,将三相相电流坐标变换到静止α-β坐标系上,如下式:
i α i β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i Ax i Bx i Cx ,
其中,iα、iβ为相电流在静止α-β坐标系上的分量;
步骤1.2,先采样电机定子电压uAB、uBC,再计算电机相电压,如下式:
u B = 1 / 3 ( u BC - u AB ) u A = u B + u AB u C = u B - u BC ,
其中,uA、uB、uC为三相相电压,之后,将三相相电压坐标变换到静止α-β坐标系上,如下式:
u α u β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 u A u B u C ,
其中,uα、uβ为定子电压在静止α-β坐标系上的分量,最后,先采样电机定子电流iA、iB、iC,再经过坐标变换得到电流在静止α-β坐标系上的分量,如下式:
i α i β = 2 3 1 - 1 2 - 1 2 0 3 2 - 3 2 i A i B i C ,
其中,iα、iβ为相电流在静止α-β坐标系上的分量;
步骤2,先用估计的电流
Figure FDA0000480886120000031
减去步骤1得到的实际电流iα、iβ,得到电流偏差
Figure FDA0000480886120000032
再分别经过符号函数sgn()运算和乘法运算,最后得到电流误差反馈量sα、sβ,如下式:
s α = Msgn ( i ‾ α ) + K ( i ‾ α ) , s β = Msgn ( i ‾ β ) + K ( i ‾ β ) ,
其中,M为滑模增益,K为状态增益,
Figure FDA0000480886120000034
步骤3,由步骤1得到的电压uα、uβ减去步骤2得到的电流误差反馈量sα、sβ,再减去第k次运算得到的估计的有效反电动势
Figure FDA0000480886120000035
得到第一组中间变量E、E,如下式:
E 1 α = u α - s α - e ^ α ( k ) E 1 β = u β - s β - e ^ β ( k ) ,
其中,E、E为第一组中间变量,k为正整数,(k=0,1,2········), e ^ β ( 0 ) = 0 ;
步骤4,用步骤3得到的第一组中间变量E、E除以电机q轴电感,得到第二组中间变量E、E,如下式:
E 2 α = E 1 α / L q E 2 β = E 1 β / L q ,
其中,Lq为电机q轴电感,E、E为第二组中间变量;
步骤5,根据估计的电机定子电流和电机定子电阻Rs以及q轴电感Lq,计算得到第三组中间变量E、E,如下式:
E 3 α = R s L q i ^ α E 3 β = R s L q i ^ β ,
其中,E、E为第三组中间变量,Rs为电机定子电阻;
步骤6,由步骤4得到的E、E减去步骤5得到的E、E,得到第四组中间变量E、E,如下式:
E 4 α = E 2 α - E 3 α E 4 β = E 2 β - E 3 β ,
其中,E、E为第四组中间变量;
步骤7,对步骤6得到的E、E进行积分,得到估计的电流
Figure FDA0000480886120000042
如下式:
i ^ α = i ^ α + T s E 4 α i ^ β = i ^ β + T s E 4 β ,
其中,Ts为采样周期;
步骤8,将步骤7得到的电流估计值
Figure FDA0000480886120000044
带入步骤2和步骤5,用于进一步估计定子电流;
步骤9,将步骤2得到的电流误差反馈量sα、sβ乘以增益N,得到第五组中间变量E、E,如下式:
E 5 α = Ns α E 5 β = Ns β ,
其中,E、E为第五组中间变量,N为有效反电动势eα、eβ滑模观测器的增益;
步骤10,估计的转速
Figure FDA0000480886120000046
分别乘以第k次运算得到的估计的有效反电动势
Figure FDA0000480886120000047
得到第六组中间变量E、E,如下式:
E 6 α = ω ^ r e ^ α ( k ) E 6 β = ω ^ r e ^ β ( k ) ,
其中,E、E为第六组中间变量;
步骤11,先由步骤9得到的E减去步骤10得到的E,得到第七个中间变量E,再由步骤9得到的E加上步骤10得到的E,得到第七个中间变量E,如下式:
E 7 α = E 5 α - E 6 β E 7 β = E 5 β + E 6 α ,
其中,E、E为第七组中间变量;
步骤12,分别对步骤11得到的E、E进行积分,得到第k+1次运算估计的有效反电动势
Figure FDA00004808861200000410
如下式:
e ^ α ( k + 1 ) = e ^ α ( k ) + T s E 7 α e ^ β ( k + 1 ) = e ^ β ( k ) + T s E 7 β ,
步骤13,先将步骤12得到的估计的有效反电动势
Figure FDA0000480886120000052
带入步骤3,进一步用于求解第一组中间变量E、E,再将步骤12得到的估计的有效反电动势
Figure FDA0000480886120000053
带入步骤10,用于计算第七组中间变量E、E
步骤14,将步骤12得到的估计的有效反电动势
Figure FDA0000480886120000054
带入锁相环,得到估计的转子位置
Figure FDA0000480886120000055
和估计的转速
Figure FDA0000480886120000056
步骤15,将估计的转速
Figure FDA0000480886120000057
带入步骤10,用于计算第七组中间变量E、E
步骤16,重复步骤1~步骤15,直至估计的电流值等于实际的电流值iα、iβ,得到与实际的反电动势eα、eβ相符的估计的有效反电动势
Figure FDA00004808861200000522
和与实际的转子位置θ相符的估计的转子位置
Figure FDA0000480886120000059
3.根据权利要求2所述的一种实现凸极永磁同步电机无速度传感器控制的方法,其特征是锁相环根据估计的有效反电动势
Figure FDA00004808861200000510
获得转子位置和估计的转速
Figure FDA00004808861200000512
的步骤如下:
步骤1,根据估计有效反电动势计算有效反电动势q轴偏差,如下式:
e ^ q = - e ^ α ( k + 1 ) cos θ ^ - e ^ β ( k + 1 ) sin θ ^ ,
其中,
Figure FDA00004808861200000515
为有效反电动势q轴偏差;
步骤2,由步骤1得到的经过比例积分调节器得到估计的转速
Figure FDA00004808861200000517
如下式:
ω ^ r = ( k i / s + k p ) e ^ q ,
其中,kp和ki分别为比例和积分系数,s为拉普拉斯算子;
步骤3,由步骤2得到的估计转速
Figure FDA00004808861200000519
通过积分得到估计的转子位置
Figure FDA00004808861200000520
如下式:
θ ^ = θ ^ + T s ω ^ r .
CN201410110944.4A 2014-03-24 2014-03-24 一种实现凸极永磁同步电机无速度传感器控制的方法 Active CN103904973B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410110944.4A CN103904973B (zh) 2014-03-24 2014-03-24 一种实现凸极永磁同步电机无速度传感器控制的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410110944.4A CN103904973B (zh) 2014-03-24 2014-03-24 一种实现凸极永磁同步电机无速度传感器控制的方法

Publications (2)

Publication Number Publication Date
CN103904973A true CN103904973A (zh) 2014-07-02
CN103904973B CN103904973B (zh) 2016-04-20

Family

ID=50996155

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410110944.4A Active CN103904973B (zh) 2014-03-24 2014-03-24 一种实现凸极永磁同步电机无速度传感器控制的方法

Country Status (1)

Country Link
CN (1) CN103904973B (zh)

Cited By (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105119538A (zh) * 2015-09-17 2015-12-02 矽力杰半导体技术(杭州)有限公司 一种永磁同步电机的驱动电路及驱动方法
CN105186946A (zh) * 2015-09-30 2015-12-23 矽力杰半导体技术(杭州)有限公司 一种高效率的永磁同步电机的驱动电路及驱动方法
CN108069021A (zh) * 2018-01-04 2018-05-25 北京中航智科技有限公司 一种舵机及其控制系统
CN108390610A (zh) * 2018-01-24 2018-08-10 天津飞旋科技有限公司 一种永磁同步电机无位置传感器控制方法
CN109067285A (zh) * 2018-08-16 2018-12-21 哈尔滨工业大学 基于零矢量电流微分的永磁同步电机无位置传感器的控制方法
CN110112965A (zh) * 2019-05-07 2019-08-09 郑州轻工业学院 一种永磁同步电机反电动势观测方法
CN110149082A (zh) * 2019-05-09 2019-08-20 浙江大学 一种电机三相电感不对称情况下的无位置传感器控制方法
CN110581679A (zh) * 2019-09-18 2019-12-17 上海中科深江电动车辆有限公司 针对三角形接法永磁同步电机无位置传感器进行控制的方法及其验证系统
CN112737450A (zh) * 2020-12-24 2021-04-30 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106026834A (zh) * 2016-07-29 2016-10-12 扬州大学 一种永磁同步电机的无速度传感器控制方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1889358A (zh) * 2006-07-11 2007-01-03 合肥工业大学 基于自适应滑模变结构的永磁同步电机控制系统
CN103346726A (zh) * 2013-07-08 2013-10-09 合肥工业大学 基于扩展磁链观测器的永磁同步电机定子磁链观测方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1889358A (zh) * 2006-07-11 2007-01-03 合肥工业大学 基于自适应滑模变结构的永磁同步电机控制系统
CN103346726A (zh) * 2013-07-08 2013-10-09 合肥工业大学 基于扩展磁链观测器的永磁同步电机定子磁链观测方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
HASEGAWA M等: "Position sensorless control for interior permanent magnet synchronous motor using adaptive flux observer with inductance identification", 《IET ELECTR.POWER APPLICATIONS》 *
MIHAICOMANESCU: "A Family of Sensorless Observers with Speed Estimate for Rotor Position Estimation of IM and PMSM Drives", 《IEEE》 *
苏健勇等: "基于四阶混合滑模观测器的永磁同步电机无位置传感器控制", 《中国电机工程学报》 *

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105119538A (zh) * 2015-09-17 2015-12-02 矽力杰半导体技术(杭州)有限公司 一种永磁同步电机的驱动电路及驱动方法
CN105119538B (zh) * 2015-09-17 2018-06-26 矽力杰半导体技术(杭州)有限公司 一种永磁同步电机的驱动电路及驱动方法
CN105186946A (zh) * 2015-09-30 2015-12-23 矽力杰半导体技术(杭州)有限公司 一种高效率的永磁同步电机的驱动电路及驱动方法
CN105186946B (zh) * 2015-09-30 2017-12-05 矽力杰半导体技术(杭州)有限公司 一种高效率的永磁同步电机的驱动电路及驱动方法
CN108069021A (zh) * 2018-01-04 2018-05-25 北京中航智科技有限公司 一种舵机及其控制系统
CN108390610A (zh) * 2018-01-24 2018-08-10 天津飞旋科技有限公司 一种永磁同步电机无位置传感器控制方法
CN109067285A (zh) * 2018-08-16 2018-12-21 哈尔滨工业大学 基于零矢量电流微分的永磁同步电机无位置传感器的控制方法
CN110112965A (zh) * 2019-05-07 2019-08-09 郑州轻工业学院 一种永磁同步电机反电动势观测方法
CN110112965B (zh) * 2019-05-07 2020-09-04 郑州轻工业学院 一种永磁同步电机反电动势观测方法
CN110149082A (zh) * 2019-05-09 2019-08-20 浙江大学 一种电机三相电感不对称情况下的无位置传感器控制方法
CN110581679A (zh) * 2019-09-18 2019-12-17 上海中科深江电动车辆有限公司 针对三角形接法永磁同步电机无位置传感器进行控制的方法及其验证系统
CN110581679B (zh) * 2019-09-18 2021-09-24 上海中科深江电动车辆有限公司 针对三角形接法永磁同步电机无位置传感器进行控制的方法及其验证系统
CN112737450A (zh) * 2020-12-24 2021-04-30 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法
CN112737450B (zh) * 2020-12-24 2022-11-04 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法

Also Published As

Publication number Publication date
CN103904973B (zh) 2016-04-20

Similar Documents

Publication Publication Date Title
CN103904973B (zh) 一种实现凸极永磁同步电机无速度传感器控制的方法
EP0944164B1 (en) Sensorless control method and apparatus of permanent magnet synchronous motor
US9059653B2 (en) Rotating electrical machine control device
CN103117703B (zh) 一种永磁同步电机无传感器控制方法及其控制装置
CN102843091B (zh) 一种永磁同步电机转子初始位置的判断方法
CN103986393B (zh) 一种永磁同步电机转子初始位置检测方法
EP2784928A2 (en) Motor control apparatus and magnetic-pole position estimating method
CN102545740A (zh) 面贴式永磁同步电机的低速无位置传感器控制方法
US9379655B2 (en) Method of field weakening control of permanent magnet motor drivers
JP2007110837A (ja) 磁極位置推定方法及びモータ制御装置
CN103178769A (zh) 永磁同步电机静止状态下参数离线辨识方法
CN104901600A (zh) 宽转速范围内永磁同步电机的无位置传感器控制方法
JP2011176975A (ja) 交流電動機の駆動装置及び電動機車両
CN102647134B (zh) 一种永磁同步电机无角度传感器的效率优化控制方法
WO2008004417A1 (fr) Appareil de commande sans capteur de machine synchrone
CN106208872B (zh) 基于免疫算法的pmsm高频电压方波注入转速估计方法
CN106655942A (zh) 永磁同步电机无位置传感器控制方法
CN103944482B (zh) 定子磁链、电磁转矩观测方法及分别应用两种方法的装置
KR20110003659A (ko) 비상 운전용 고효율 영구자석 동기모터의 각도위치 센서리스 제어 방법
CN106374805B (zh) 永磁同步电机转子静止状态初始位置辨识方法及装置
US20060052972A1 (en) Shaft sensorless angular position and velocity estimation for a dynamoelectric machine based on extended rotor flux
CN106059435B (zh) 一种提高永磁同步电机转子位置估计精度的方法
CN105024615A (zh) 一种永磁同步电机低速无传感器控制方法及装置
CN112204869A (zh) 电力转换装置
CN104094517B (zh) 控制永磁电机

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20210210

Address after: 245000 No.50, Meilin Avenue, Huangshan Economic Development Zone, Huangshan City, Anhui Province

Patentee after: Huangshan Development Investment Group Co.,Ltd.

Address before: 230009 No. 193, Tunxi Road, Hefei, Anhui

Patentee before: Hefei University of Technology

TR01 Transfer of patent right