杨洛鸿 刘 侃 胡 伟 周世超 黄 庆 陈泳丹
(1.湖南大学机械与运载工程学院 长沙 410082;2.株洲中车时代电气股份有限公司 株洲 412001;3.中国北方车辆研究所 北京 100072)
永磁同步电机(Permanent magnet synchronous motor,PMSM)因其具有体积小、功率密度高以及转矩密度高等优点[1],广泛应用于工业伺服控制、轨道交通以及风力发电等领域[2-4]。准确的转子位置信息是PMSM高效控制的基础[5],但电机拆装维修、位置传感器故障[6]、驱动控制器更换以及软件数据更新等问题都可能导致转子初始位置信息丢失,进而造成控制系统失效。
采用无位置传感器算法或位置传感器[7-9]是现有获取转子位置信息的两种主要方式。无位置传感器算法通常较复杂,且存在额外损耗、参数敏感等问题[10]。基于高频注入的无位置传感器控制算法通常适用于凸极式永磁同步电机的初始位置判定和中低速运行[11],而基于反电势的无位置传感器方法对电机参数随磁饱和的变化比较敏感[12]。现有的无位置传感器算法通常无法获取隐极式永磁同步电机的准确位置信息。因此,采用位置传感器仍是获取转子精确位置的主要方式。
其中,增量式编码器以其高精度、高可靠性、高性价比等优势[13],广泛应用于高精度控制等领域。增量式编码器通过计算脉冲数量来获取位置,无法直接获得转子的绝对位置信息,在使用前通常需要进行零位校正。目前,主要通过定子电流矢量法对增量式编码器进行零位校正[14-18]。
文献[19]中,通过定子电流矢量拖动转子到达零位,进而确定转子初始位置,对编码器信息清零后,通过首个z脉冲信息计算零位偏置角。文献[20]中,采用了六种校正模式,将转子分别拖动至六个确定位置,通过任一校正模式的编码器信息与实际转子位置信息间的偏差即可计算零位偏置角。然而,当系统软硬件故障导致零位偏置角出错或丢失时,利用现有方法通常无法实现快速零位校正,恢复系统正常运行。
针对该问题,本文提出了一种基于交轴磁链辨识的永磁同步电机位置传感器零位校正方法,该方法既可用于隐极式电机,也能用于凸极式电机。本文章节安排如下:第2节介绍了零位校正原理与位置误差的影响;第3节介绍了所提的零位校正方法原理与具体实现流程;第4节介绍了试验测试平台构成,并通过试验测试证明了所述方法的可行性与可靠性。试验结果表明,本文提出方法能够实现转子快速定位,进而能对增量式编码器进行零位校正,且启动过程不受转子初始位置误差的影响。
如图1所示,在静止坐标系下,当转子中心轴d轴与α轴重合时,此时转子的电气角度θf为0°,定义为零位1。
图1 零位校正示意图
增量式编码器输出信号通常包括两路正交计数脉冲与一路z脉冲。为了防止计数脉冲丢失产生累积误差,转子每旋转一周产生一个z脉冲,对正交脉冲计数清零。定义z脉冲出现时,转子中心轴位置在静止坐标系中的映射为d0轴。此时,编码器输出的电气角度θen为0°,定义为零位2。θen可以表示为
式中,θm为编码器输出的机械角度信息;Pn为极对数。
在增量式编码器安装过程中,通常无法将零位1与零位2对齐,这会导致零位偏置角θ0的存在。如图2所示,如果不进行零位校正,则系统位置信息θe就等于编码器输出角度θen,实际控制的同步旋转坐标系与转子中心轴不重合,影响控制效果。
图2 零位校正对控制效果的影响
只有在编码器输出角度θen的基础上校正零位偏置角θ0,才能使控制系统通过增量式编码器获得转子绝对位置信息,进而实现正确的坐标变换。当零位校正足够准确时,可以认为校正后的位置信息即为转子真实位置信息,即θe=θen–θ0=θf。
因此,零位校正的目的就是通过获取零位偏置角θ0对θen进行校正,进而获得转子绝对位置信息,实现PMSM的准确高效控制。
如图3所示,定义准确零位校正后的同步旋转坐标系为dq,此时θe=θen–θ0=θf,即位置信息准确工况。定义任意位置信息下的同步旋转坐标系为d'q′,此时θe=θb,即位置信息不准确工况。两种工况的同步旋转坐标系间存在位置误差θa。
图3 位置信息准确与不准确工况下的稳态相量图
dq坐标系下,永磁同步电机稳态电压方程为
式中,ud、uq为d、q轴电压;id、iq为d、q轴电流;Ld、Lq为d、q轴电感;ψf为转子磁链;R为定子电阻;ωe为电角速度;ψd、ψq为d、q轴磁链。
d ′q′坐标系下,永磁同步电机稳态电压方程为
式中,i'd、i'q为d'、q'轴电流。
id和iq可以分别表示为
将式(4)代入式(3)可以得到
式中,ψ'd、ψ'q为d'、q'轴磁链。
当id=id′= 0,iq=iq′时,由式(2)、式(5)可得到
通过式(6)可以看到,此时位置误差θa的存在导致两坐标系下交轴磁链大小不同,进而导致直轴电压大小不同,在转速相同的情况下,只有当θa=0时,直轴电压才能相等,此时,θb=θf则为u'd=ud的充分必要条件。
基于这一结论,当id=i'd=0时,如果已知dq坐标系下交轴磁链ψq的变化方程,通过式(6)则可以得到任意转速下的ud,在任意转子位置信息θb的控制下,将ud作为参考值,当u'd=ud时,则有θb=θf,即可获得准确的转子位置信息,实现转子定位的目的,对增量式编码器进行零位校正。因此,本文提出了基于交轴磁链辨识的位置传感器零位校正方法。
将式(2)采用欧拉方法离散化,可以得到
式中,k是系统的离散采样时刻。由于矢量控制系统通常采用脉宽调制技术进行电压输出,难以测量实际电压,所以通常采用参考电压代替。但实际电压通常会受到逆变器非线性因素的影响,因此,式(7)可以表示为
式中,Vdead Dd和Vdead Dq分别表示逆变器非线性因素引起的dq轴畸变电压。Vdead可以表示为[20]
结合式(9)可知,Vdead与以下因素有关:开关延时Ton和Toff、死区时间Td,以及IGBT和二极管导通压降Vce0和Vd0。
Dd和Dq与三相电流方向和系统电气角度有关[21],Dd和Dq可以表示为
其中
式(8)是秩数为 2 的方程组,但包含R、Vdead、ψd、ψq四个未知变量,无法直接求解得到ψd或ψq。因此,文献[22]提出了一种永磁同步电机dq轴磁链辨识的新方法,该方法的辨识模型消除了逆变器非线性因素以及电阻的影响,从而提高了辨识精度。
在文献[22]中,系统采用电流环控制模式,通过调整负载大小改变电机转速,实现连续变速工况,然后采集两组不同转速下的运行数据,并对基于智能算法的辨识模型进行最小化运算,实现dq轴磁链的准确辨识。然而,该方法的试验工况和模型求解过程较复杂,因此本文对二者进行了简化,采用的试验工况如图4所示,将原方法中的连续变速工况改为一次稳态变速工况,进而可以简化辨识模型和求解过程。下文仅对交轴磁链辨识过程进行介绍。
图4 稳态变速工况示意图
如图4所示,分别采集并存储两组采样时间同为t0,转速不同的稳态运行数据,称为data1和data2。因此,通过式(8)得到两个离散的d轴电压方程,可以表示为
式中,i、j分别为data1、data2的数据下标,且i=j,将式(12)中两方程相减,可以得到
data1、data2中每一组对应的数据点都满足式(13),因此将式(13)进行求和,可以得到
式中,N是data1和data2的数据长度。当采用电流环控制时,ψq、Vdead通常被认为是常数[23],并且通过式(10)可以看到,Dd仅与电流方向和电气角度有关。因此,可以得到
因为data1和data2均为稳态运行数据,所以式(14)可以简化为
如图5所示,通过定子电流矢量法对编码器进行零位校正,此时系统位置信息θe=θen–θ0,采用id=0控制,设定iq在0~b范围内,通过调整步长a改变iq大小。采用第3.1节所述方法,辨识不同iq下的磁链ψq。采用定子电流矢量法得到的零位偏置角θ0足够准确,可以认为校正后的系统位置信息即为转子真实位置信息,此时θe=θen–θ0=θf,进而可以获得准确定位工况下的q轴磁链曲线拟合方程ψq(iq)。通过该方程,可以得到dq坐标系下任意大小iq所对应的ψq。
图5 获取交轴磁链曲线拟合方程的方法流程图
如图6所示,为提出方法的系统控制框图,该方法采用电流环控制,且令id=0。如第2.2节所述,u′d=ud为θb=θf的充分必要条件,基于这一结论可采用两种控制模式,分别为模式1与模式2,实现基于交轴磁链辨识的增量式编码器零位校正。
图6 基于交轴磁链辨识的零位校正方法控制框图
在模式1中,将比例积分(PI)控制器的输出变量θpi1直接作为系统位置信息,即θe=θpi1;在模式2中,则将PI控制器的输出变量θpi2作为编码器位置信息θen的补偿,并将二者之差作为系统位置信息,即θe=θen–θpi2。在收敛过程中,θpi1与θpi2均大小任意,因此由图3可知,θe等价于θb,此时系统控制的同步旋转坐标系可以表示为d ′q′,则由式(6)可以得到
同时通过第3.2节所述方法,计算获得图3中dq坐标系下的q轴磁链曲线拟合方程ψq(iq),同样,由式(6)可以得到
式中,udfit为拟合电压,即准确定位工况下的d轴电压,iq为q轴实际电流,ωe通过增量式编码器获取。将udfit作为系统d轴电压ud*的参考值,并将ud*输入控制系统,与滤波后的udfit比较,将二者的误差Δud作为控制量输入到PI控制器,当Δud减小时,ud*逐渐收敛为udfit,输出量θpi1或θpi2滤波后参与坐标变换,同时θe逐渐收敛为θf,经变换得到的id、iq参与udfit和ud*的运算,从而实现闭环反馈控制。
如图7所示,当系统收敛稳定时,可以得到=udfit,此时在模式1中可以得到θpi1=θf,同时,通过增量式编码器不能获取准确的转子绝对位置信息,因此令 Δθ=θen–θpi1,将 Δθ称为偏置角,显然 Δθ即可作为新的零位偏置角,然后切换系统位置信息为θ1实现零位校正。在模式 2中可以得到θ2=θen–θpi2=θf,显然,此时θpi2即为新的零位偏置角。因此采用两种模式均可实现增量式编码器的零位校正。
图7 基于交轴磁链辨识的零位校正方法流程图
试验装置如图8所示,采用Speedgoat快速原型试验平台,主要包括一台三相PMSM,一台直流电机以及相应的控制系统,直流母线电压设定为40 V,采用2 500PPI分辨率的增量式编码器获取转子位置信息,PMSM参数如表1所示。
表1 永磁同步电机参数
图8 永磁同步电机快速原型试验平台
采用前文所述方法获取PMSM的q轴磁链曲线拟合方程,具体方法如下:PMSM采用电流环控制模式且id=0,iq从0递增至3 A,每增加0.5 A采集一组数据。采用直流电机控制转速,对应于每一个iq,先调整为低转速采集数据data1,然后提高转速采集数据data2,完成一组数据的采集,数据采集时间t0为 5 s。
需要说明的是,文献[22]采用连续变速工况试验,辨识结果不受转速波动影响,数据采集时间仅为1 s。本文采样频率与文献[22]相同,为防止速度波动对辨识结果造成影响,数据采集时间增加为5 s。
根据采集的稳态运行数据,通过式(16)计算对应iq下的q轴磁链ψq,并使用 Matlab中的 Curve Fitting Tool进行数据拟合。拟合方程ψq(iq)如式(19)所示,拟合曲线如图9所示,由图9可知,数据均匀分布在曲线附近,拟合效果良好。
图9 q轴磁链拟合曲线
零位校正试验采用电流环控制模式空载运行,设置iq=1 A,id=0 A。PMSM由静止状态起动,并加速至稳定运行。
4.3.1 模式1
图10为采用模式1进行零位校正时的试验波形。由图10可知,0~8 s为收敛阶段,同时为起动加速阶段,d轴参考电压ud*从0开始波动,逐渐收敛到拟合电压udfit。在8~12 s,ud*围绕udfit稳定波动,经计算得到平均值同时电流、转速稳定波动,说明8 s后系统已经收敛到稳定状态。
图10 模式1试验波形图
如图11所示,θf为采用定子电流矢量法得到的转子真实位置信息,θf=θen–θ0,在图中的表示为θen与θf间的稳定偏差。
如图11b所示,在收敛阶段,θpi1从0开始变化,且与θf间存在明显位置偏差,经过8 s逐渐收敛为θf。如图11c所示,在收敛稳定阶段,此时θpi1与θf十分接近,θpi1与θen间同样存在稳定偏差,即偏置角 Δθ。
图11 控制器输出位置信息θpi1、转子真实位置信息θf与编码器输出位置信息θen对比图,(b)、(c)为(a)的局部放大图
令θep=θf–θe,其中θep表示转子位置误差,为系统位置信息θe与转子真实位置信息θf间的误差,因此对于模式 1 有θep=θf–θpi1。
如图12所示,θep最终收敛在0附近,通过图中局部放大可以看到,θep在0附近波动,这一波动主要是由转速波动引起,且变化幅值较小,经计算可得,稳定阶段误差均值=1.1°,即机械角度0.275°,这一误差对控制效果的影响很小,可以认为此时θpi1完全收敛,通过PI控制器获得了SPMSM的转子真实位置信息,实现了隐极式电机转子定位。
图12 模式1转子位置误差θep
采用θpi1作为系统位置信息只能在特定工况下运行,无法适应全部工况,因此,仍需要校正后的编码器信息参与系统控制。在θf未知的情况下,可以通过模式 1获取θpi1,当收敛完成时,θpi1与θen间存在偏置角Δθ,显然Δθ即为新的零位偏置角,可以用于校正θen,从而获取准确的转子位置信息θ1,然后将系统位置信息θe从θpi1切换为θ1,完成模式1的增量式编码器零位校正。为消除转速波动的影响,求取Δθ的平均值作为新的零位偏置角,即θ1=θen–。
如图13所示,θ1为校正后位置信息,θen为校正前位置信息,θf作为位置信息参考,可以看到零位校正前,θen与θf间存在稳定位置偏差,零位校正后,θ1与θf几乎完全重合,校正误差即为=1.1°。综上所述,采用模式1可以实现编码器零位校正,且校正效果良好。
图13 模式1校正后位置信息θ1、转子真实位置信息θf与编码器输出位置信息θen对比图
4.3.2 模式2
图14为采用模式2零位校正时的试验波形。由图14可知,0~10 s为系统收敛阶段,相较于模式1,起动加速时间很短,同样d轴参考电压从0开始波动,逐渐收敛到拟合电压udfit。在10~14 s,围绕udfit稳定波动,经计算得到=udfit=–0.027 V,同时电流、转速稳定波动,说明10 s后系统已经收敛到稳定状态。
图14 模式2试验波形图
如图15所示,采用模式2时,θ2为校正后位置信息,θ2=θen–θpi2,θen为校正前位置信息,θf作为位置信息参考。
图15 控制器输出位置信息θpi2、转子真实位置信息θf与编码器输出位置信息θen对比图,(b)、(c)为(a)的局部放大图
如图15b所示,在收敛阶段,θpi2初始为0,θ2与θen十分接近,且与θf间存在明显偏差,随着θpi2收敛,θ2与θen间产生明显偏差,同时θ2逐渐收敛为θf。如图15c所示,在收敛稳定阶段,θ2与θf几乎完全重合。
如图16 所示,对于模式 2 有θep=θf–θ2,θep最终收敛在 0附近,经计算稳定阶段误差均值=1.6°,即机械角度 0.4°,且波动幅值仅为0.2°,可以认为此时θ2完全收敛。因此,采用模式2可以获取转子真实位置信息,同时也可以实现增量式编码器零位校正,θpi2即为新的零位偏置角。
图16 模式2转子位置误差θep
综上所述,当PMSM因为系统软硬件故障,导致原有的零位偏置角θ0出错或丢失时,即可采用本方法对增量式编码器进行快速零位校正。采用模式1时,空载运行8 s左右即可获取准确的转子位置信息θpi1,进而求取作为新的零位偏置角,从而校正编码器获取准确的转子位置信息θ1。采用模式 2时,空载运行10 s左右即可获取准确的转子位置信息θ2,θpi2即为新的零位偏置角,用于校正编码器信息θen。两种模式均可实现位置传感器的快速零位校正,试验证明了本方法的可行性。
由于 PI控制器没有设定初始值,在起动时刻θpi1=θpi2=0°,但转子初始位置是任意的,则可能存在由于转子初始位置误差太大,导致电机起动受到影响的情况。因此对于两种模式,分别调整转子位置误差θep为10°、50°、90°、130°和170°进行起动,测试初始位置误差对起动效果的影响。
4.4.1 模式1
如图17所示,10~14 s时θep均能收敛到0附近,电机起动及运行不受初始位置误差的影响。
图17 模式1转子位置误差θep,起动时刻θep分别为10°、50°、90°、130°和170°
由图17可知,在起动瞬间,转子位置误差突变到90°左右,这一位置突变主要由iq引起。在起动时刻,iq矢量在同步旋转坐标系中的位置为90°,因此在起动瞬间电流将转子牵引至90°附近,造成这一位置突变。这就导致了不同初始位置误差下,θep变化基本相同,同时也表明即使初始位置误差大于90°,电机仍然能够保持正转。
4.4.2 模式2
如图18所示,当转子初始位置误差分别为10°和50°时,采用模式2均能正常收敛,10~14 s时θep收敛到0附近,且电机正转。
图18 模式2转子位置误差θep,起动时刻θep分别为10°和50°
如图19所示,当转子初始位置误差分别为130°和170°时,采用模式2均能正常收敛,10~14 s时θep收敛到180°附近,说明此时转子真实位置信息θf与系统位置信息θe相差180°,此时电机反转。
图19 模式2转子位置误差θep,起动时刻θep分别为130°和170°
如图20所示,当转子初始位置误差为90°时,此时iq与转子同向,导致起动转矩很小,因此出现了以下三种结果:当电机正转时,θep收敛到0附近;当电机反转时,θep收敛到180°附近;当电机无法起动时,θep无法收敛。
图20 模式2转子位置误差θep,起动时刻θep为90°
经多次试验,当初始位置误差在 79°~103°范围内时,因起动转矩太小,均存在无法起动的可能,因此模式2存在起动死区,转子初始位置误差对起动过程有一定影响。
综上所述,采用本文提出的方法进行零位校正时,当使用模式1时,任意初始位置误差工况下,系统都能够正常起动并稳定收敛,但在起动时刻会存在转子位置突变的现象;当使用模式2时,在较大初始位置误差范围内,系统都能够正常起动并稳定收敛,但当电流矢量与转子间夹角过小,导致起动转矩太小时,存在无法起动的问题。试验验证了本方法的可靠性。
在id=0工况下,位置信息相同是直轴电压相等的充分必要条件。利用这一结论,本文通过两种控制模式实现了增量式编码器的快速零位校正。
(1) 对于模式1,无需位置传感器信息,利用控制器即可输出转子准确位置信息,进而可以对编码器进行零位校正。该模式起动不受转子初始位置误差影响,且对电机凸极性没有要求,既可用于凸极式电机,也可用于隐极式电机。
(2) 对于模式2,可以通过控制器直接获取零位偏置角,在线补偿转子位置信息,从而实现编码器零位校正。但该模式存在一定的起动死区,起动过程会受转子初始位置误差影响。
(3) 通过零位校正试验与任意初始位置误差起动试验分别说明了本文提出方法的可行性与可靠性。结果表明,该方法定位效果良好,零位校正速度快、精度高,可靠性高,具有一定的工程应用价值。