基于逐相旋转坐标变换的无刷直流电机转子磁场定向瞬时转矩控制技术

李珍国 王鹏磊 孙启航 贾益丞

(燕山大学电力电子节能与传动控制河北省重点实验室 秦皇岛 066004)

摘要 转子磁场定向无刷直流电机等效模型所具有的d轴转子磁链为恒值、q轴转子磁链始终为零的特征,是构建高性能无刷直流电机瞬时转矩控制系统的有利条件。但是传统的等功率3s/2r坐标变换无法使无刷直流电机得到转子磁场定向等效模型。为此,提出一种基于逐相旋转坐标变换的转子磁场定向方法,即对每相转子磁链逐一进行旋转坐标变换得到各相绕组的转子磁场定向等效模型,再以功率守恒为原则进行三相合成,得到无刷直流电机转子磁场定向等效模型。鉴于该转子磁场定向等效模型中电磁转矩仅与q轴电流成正比而与d轴电流无关,理论推导出直轴电枢反应磁动势、瞬时电磁转矩各自与d、q轴电流之间的关系式,构建了具有无增磁/去磁电枢反应的瞬时转矩控制系统。最后,通过Matlab仿真和DSP驱动实验,验证了所提控制技术的可行性和有效性。

关键词:无刷直流电机 转子磁场定向 逐相旋转坐标变换 无增磁和去磁 瞬时转矩控制

0 引言

无刷直流电机(Brushless DC Motor, BLDCM)因其结构简单、功率密度高、寿命长等优点,被广泛地应用于新能源汽车、医疗器械、航空航天等领域[1-3]。方波电流驱动方式凭借着实现简单、经济可靠等优点,在无刷直流电机控制系统中得到广泛应用,但是由于非理想反电动势的作用以及绕组电感的存在,使无刷直流电机在运行过程中普遍存在较大转矩脉动,其中换相转矩脉动最为严重,最大可达约50%的平均转矩[4],严重制约了其在高精度场合的推广和应用[5]

国内外相关学者针对如何有效抑制无刷直流电机的转矩脉动进行了大量研究,主要从优化电机本体设计[6-7]和优化控制策略两方面出发来抑制无刷直流电机的转矩脉动。其中基于优化控制策略的措施主要包括换相期间与非换相期间采用不同的调制策略[8-10]、在逆变器前端添加辅助电路[11-14]等,其主要思想是平衡换相期间开通/关断相电流的变化速率,维持非换相相电流稳定。此外将直接转矩控制方法应用到无刷直流电机的控制系统中,也能够取得良好的转矩脉动抑制效果[15-17]。考虑到换相期间三相绕组均参与能量转换,一些学者将坐标变换的思想引入无刷直流电机的换相过程控制中[18-19]。基于ab 坐标系对换相过程的分析,文献[18]提出一种电压矢量注入方法,通过注入与换相时刻反电动势矢量同向或反向的电压矢量,使换相期间的电流矢量切换轨迹正交于反电动势矢量,达到减小换相转矩脉动的目的;文献[19]也从Clarke变化入手,将换相过程的控制转化为ab 坐标系下两轴电流分量的控制,分别作用于转矩和换相时间,能够实现快速平稳换相。无刷直流电机控制系统的主电路一般采用三相电压型逆变器,除方波电流驱动方式外,在不改变逆变器拓扑结构的前提下,也可以采用三相导通方式驱动电机运行。具体选用何种驱动方式主要由无刷直流电机需要的性能决定。越来越多的研究发现,三相导通方式不仅不存在换相问题,而且在改善电机运行性能方面同样具备优势[20-24]。文献[20]分析了正弦波电流驱动方式对无刷直流电机机械特性和转矩性能的影响,结果表明,在高速区运行时正弦波电流驱动相比于方波电流驱动表现出更“硬”的机械特性;文献[21]采用高次谐波电流注入的方法来选择性地消除转矩中的谐波分量,从而减小转矩脉动,但该方法需要对反电动势进行傅里叶级数分解,较为复杂。为了避免复杂的复指数分解,一些学者提出从矢量的角度对无刷直流电机的运行过程进行分析。文献[22]以转矩脉动最小为控制目标,根据在ab 坐标系下得到的理想反电动势矢量轨迹,通过让参考电流矢量对齐反电动势矢量来得到最优参考电流;文献[23]通过推导dq旋转坐标系下转矩和绕组铜耗各自与d、q轴电流之间的表达式,利用拉格朗日定理直接求解出给定转矩下的最小电流矢量,抑制转矩脉动的同时也降低了绕组铜耗,但该控制方法未能做到d、q两轴之间的转矩解耦;文献[24]以绕组变换前后产生的电磁转矩和直轴电枢反应不变为依据,把三相绕组电流变化为d、q轴电流,分别用于控制直轴电枢反应和转矩,但并未给出直轴电枢反应模型。

本文提出一种基于逐相旋转坐标变换的无刷直流电机转子磁场定向瞬时转矩控制技术。该技术首先通过对每相绕组单独进行旋转坐标变换得到无刷直流电机的转子磁场定向等效模型。其次,在该转子磁场定向dq旋转坐标系下,根据电磁转矩仅与q轴电流成正比而与d轴电流无关的特点,理论推导出直轴电枢反应磁动势、瞬时电磁转矩各自与d、q轴电流之间的关系式,并以此构建了无刷直流电机转子磁场定向瞬时转矩控制系统。所提控制技术不存在换相过程和状态切换,具有较强的适应性和灵活性,并通过相应的Matlab仿真和DSP驱动实验,验证了其可行性和有效性。

1 转子磁场定向下无刷直流电机的等效模型

图1给出了无刷直流电机三相静止坐标系和转子磁场定向dq同步旋转坐标系下的等效模型示意图。图中,A、B、C轴分别处于三相绕组轴线方向;d轴始终置于转子永磁体N极所指方向,也是d轴绕组轴线方向;q轴始终置于空间上超前d轴90°电角度方向,也是q轴绕组轴线方向;d、q轴随转子永磁体同步旋转;q 为以A相绕组轴线为参考的当前时刻转子磁场位置角。由图1可知,在转子旋转过程中,无刷直流电机转子磁场定向dq同步旋转坐标系下的等效模型总是具有转子永磁体磁场与d轴绕组交链的转子磁链yd恒为最大值、q轴绕组交链的转子磁链yq始终为0的特征。

width=152.5,height=179.65

图1 三相静止坐标系与转子磁场定向dq同步旋转坐标系等效模型示意图

Fig.1 Schematic diagram of equivalent model in three-phase stationary coordinate system and rotor field oriented d-q synchronous rotating coordinate system

理想情况下,无刷直流电机转子永磁体产生的磁场沿定子内表面呈梯形波分布,导致在转子旋转过程中,其与A、B、C相绕组分别交链的转子磁链yAyByC为非正弦波。采用传统等功率3s/2r坐标变换矩阵,将该三相转子磁链坐标变换至d、q两轴转子磁链,则有

width=191,height=95 (1)

图2给出了理想梯形波磁场分布下的无刷直流电机三相转子磁链波形及通过式(1)得到的d、q轴转子磁链波形。图中,ymax表示转子永磁体磁场与每相绕组交链的转子磁链的最大值。由图2可知,采用式(1)中传统的等功率3s/2r坐标变换得到的d轴转子磁链yd大小非恒定、q轴转子磁链yq大小不等于0。显然,该两轴磁链特征与图1中转子磁场定向dq同步旋转坐标系下的两轴磁链特征不相符。由此可知,采用传统的等功率3s/2r坐标变换不能实现无刷直流电机转子磁场定向。

width=233.9,height=175.95

图2 理想梯形波磁场分布下三相转子磁链波形和d、q轴转子磁链波形

Fig.2 Three-phase rotor flux linkage waveform and d、q axis rotor flux linkage waveform under ideal trapezoidal wave magnetic field distribution

2 基于逐相旋转坐标变换的转子磁场定向瞬时转矩控制理论

2.1 逐相磁链旋转坐标变换

由前节分析可知,采用传统等功率3s/2r坐标变换不能实现无刷直流电机的转子磁场定向,其主要原因是无刷直流电机转子永磁体产生的磁场在空间上呈非正弦波分布,致使各相转子磁链为非正弦波的缘故。为此,本文拟采用一种逐相磁链旋转坐标变换来实现无刷直流电机转子磁场定向。

考虑到转子永磁体磁场与各相绕组交链的转子磁链为非正弦波,为了获取图1中dq同步旋转坐标系下yd恒为最大值、yq始终为0的特征,实现无刷直流电机转子磁场定向,本文对每相转子磁链单独进行旋转坐标变换。以A相为例,在yA大小不变的前提下,通过人为地调整A相绕组位置至width=13.95,height=12轴位置,使在qA角度下,yAymax之间满足余弦关系,其中qA表示从width=13.95,height=12轴方向到转子永磁体磁场方向的夹角。在qA角度下,width=13.95,height=12轴转子磁链width=17,height=15满足

width=96,height=15 (2)

建立dAqA旋转坐标系,取转子永磁体磁场方向为dA轴,空间上超前dA轴90°电角度的方向为qA轴。为将A相绕组等效旋转坐标变换至dA、qA轴绕组,虚拟构造出空间上超前width=13.95,height=12轴90°电角度的width=15,height=12轴及其绕组,两轴组成width=28,height=12坐标系。图3给出了某一转子位置下,width=28,height=12坐标系与dAqA旋转坐标系之间的位置关系示意图。为使虚拟构造的width=15,height=12轴绕组不产生额外的电磁转矩,故置width=15,height=12轴绕组电流width=13.95,height=15=0。由图3可知,在时间上width=15,height=12轴转子磁链width=18,height=15滞后width=13.95,height=12轴转子磁链90°电角度,因此,当前转子位置下width=18,height=15满足

width=72,height=15 (3)

width=157.3,height=139.7

图3 width=24.95,height=11坐标系与dA-qA旋转坐标系位置关系示意图

Fig.3 Diagram of relationship between width=24.95,height=11 coordinate system and dA-qA rotating coordinate system

width=28,height=12坐标系到dAqA旋转坐标系的变换矩阵CA表示为

width=100,height=33 (4)

采用式(4),对width=28,height=12坐标系下的两轴转子磁链进行坐标变换,可得到变换至dA、qA轴后的两轴转子磁链大小为

width=114.95,height=33 (5)

由式(5)可知,旋转变换至dA轴的转子磁链ydA恒为ymax、qA轴的转子磁链yqA始终为0。该两轴磁链特征与图1中转子磁场定向dq同步旋转坐标系下的两轴磁链特征相吻合。由此可知,建立的dAqA旋转坐标系与转子磁场定向dq同步旋转坐标系互相重合。

同理,对于B、C相也存在各自的width=26,height=12width=26,height=13坐标系和变换矩阵CBCC,其中从width=13,height=12width=13,height=13轴方向到转子永磁体磁场方向的夹角分别为qBqC。采用变换矩阵CBCC分别对width=26,height=12width=26,height=13坐标系下的两轴转子磁链进行坐标变换,也可得到旋转变换至dB、dC轴的转子磁链ydBydC恒为ymax;qB、qC轴的转子磁链yqByqC始终为0的结果。不难看出,dBqB、dCqC旋转坐标系也与转子磁场定向dq同步旋转坐标系互相重合。

2.2 转子磁场定向dq同步旋转坐标系下的瞬时转矩方程

以A相为例,width=28,height=12坐标系下的瞬时电磁转矩方程表示为

width=102,height=35 (6)

式中,np为电机极对数;width=13,height=15width=13.95,height=12轴绕组电流,与A相绕组电流iA在数值上相等;width=13.95,height=15width=15,height=12轴绕组电流,为使虚拟构造的width=15,height=12轴绕组不产生额外的电磁转矩,其大小始终为零。

采用式(4)给出的坐标变换矩阵CA,对式(6)进行坐标变换,可得到转子磁场定向dq同步旋转坐标系下的A相瞬时电磁转矩方程为

width=147,height=36 (7)

式(7)可进一步的表示为

width=180,height=35 (8)

式中,idAiqA分别为A相电流坐标变换至d、q轴后的电流;kA为dqA/dq

经推导,B、C相可得出类似结果,因此,A、B、C三相绕组在dq同步旋转坐标系下产生的合成瞬时电磁转矩T大小为

width=186.95,height=20 (9)

式中,iqBiqC分别为B、C相电流坐标变换至q轴后的电流;kB、kC分别为dqB/dqdqC/dqiq为q轴绕组电流。

2.3 dq两轴电流对应三相绕组电流计算

由式(9)可知,若能够做到iq为恒值则可以实现恒定瞬时转矩控制。根据式(9)中iqiqAiqBiqC之间的关系,可对应得到d轴绕组电流ididAidBidC之间的关系为

width=105,height=17 (10)

式中,idBidC分别为B、C相电流坐标变换至d轴后的电流。

根据2.1节中给出的坐标变换矩阵,在已知dAqA、dBqB、dCqC旋转坐标系下的两轴电流时,能够直接利用各相坐标变换矩阵的逆矩阵逐相求解出A、B、C相绕组电流。因此,求解三相绕组电流的关键在于得到各相旋转坐标系下的两轴电流。

根据式(9)、式(10)及各轴电流之间的约束关系,可得

width=175.95,height=78.95 (11)

由式(11)可知,为了能够计算获得各相旋转坐标系下的两轴电流,须同时给出d、q轴绕组电流。其中,iq可通过电磁转矩计算得到;id不参与产生电磁转矩,在恒定瞬时转矩控制中,作为一个额外的自由度可以给定为零,或根据其他控制需要得到。

在已知d、q轴电流的前提下,由式(11),可得到A、B、C相旋转坐标系下的两轴电流与d、q轴绕组电流之间的关系为

width=193.95,height=33(12)

其中

HAd=-kBsinqB+kCsinqC HAq=kBcosqB-kCcosqC

HBd=-kCsinqC+kAsinqA HBq=kCcosqC-kAcosqA

HCd=-kAsinqA+kBsinqB HCq=kAcosqA-kBcosqB

H=kAsinqAHAq+kBsinqBHBq+kCsinqCHCq

=kAcosqAHAd+kBcosqBHBd+kCcosqCHCd

采用CACBCC的逆矩阵分别对A、B、C相旋转坐标系下的两轴电流进行坐标反变换,即可得到静止坐标系下的各相绕组电流大小为

width=154,height=27 (13)

由式(13)可知,完全可以把每相绕组中的电流分为两部分,与id有关的第1项定义为可调电流分量,与iq有关的第2项定义为转矩电流分量,分别用ixdixq表示。

2.4 d轴电流给定值计算

电机负载运行时,电枢绕组电流产生的电枢反应磁动势,其直轴分量直接对由转子永磁体磁场主导的气隙磁场起到增磁或去磁作用。若能够做到该直轴分量等于零,则既可以避免因增磁引起的磁路饱和,又能够避免转子永磁体失磁风险。由2.3节分析可知,在恒定瞬时转矩控制中,id作为一个额外的自由度,可以用来实现直轴电枢反应磁动势等于零,即无增磁和去磁控制。

当给定id=0时,可以做到每相绕组电流中仅存在转矩电流分量,但不一定能够实现直轴电枢反应磁动势等于零。其理由如下:当三相绕组中分别流过大小为iAqiBqiCq的三相电流时,各相绕组电流产生的电枢反应磁动势在空间上按余弦规律分布,其大小为

width=114.95,height=101 (14)

式中,fAqfBqfCq分别为iAqiBqiCq产生的电枢反应磁动势;a 为沿气隙圆周的角度;N为每相绕组的有效串联匝数。

以各相绕组轴线方向作为fAqfBqfCq的空间矢量方向,对该三相磁动势进行空间矢量合成可得到三相合成磁动势,并将该合成磁动势投影到d轴即可得到直轴磁动势fd1大小为

width=87,height=31 (15)

其中

H11=HAqcosq+HBqcos(q-120°)+HCqcos(q+120°)

图4给出了id=0时,一个电周期内直轴磁动势fd1的变化波形。由图4可知,fd1大小非恒等于零,故id=0时,未能在任意转子位置下做到无增磁和去磁。在此基础之上,若通过灵活的给定id值使iAdiBdiCd产生的直轴电枢反应磁动势恰好能够抵消fd1的增磁或去磁效果,即可在任意转子位置下实现无增磁和去磁控制。

width=212.9,height=62.65

图4 id=0条件下的直轴磁动势

Fig.4 Direct axis magnetomotive force (id=0)

同理,当三相绕组中分别流过大小为iAdiBdiCd的三相电流时,该三相绕组电流产生的合成磁动势投影到d轴后的直轴磁动势fd2大小为

width=82,height=31 (16)

式中

H12=HAdcosq+HBdcos(q-120°)+HCdcos(q+120°)

若要在任意转子位置下实现无增磁和去磁控制,须满足fd1+fd2=0,定义kr为d、q轴电流的比例系数,则此时的kr满足

width=59,height=31 (17)

图5给出了恒转矩控制条件下,无增磁和去磁时所需的d轴电流波形及通过式(13)求得的三相电流波形。由图5可知,为了在任意转子位置下做到无增磁和去磁,逆变器需要工作在三相导通方式。

根据图5中的三相电流波形特点,可将一个周期的转子位置划分成Ⅰ~Ⅵ六个扇区,如图5所示,在每个扇区范围内,三相电流总是具有一相电流大小近似不变且方向与其他两相电流相反的特点。

width=209.15,height=72.25

图5 无增磁和去磁控制所需的d轴电流和三相电流

Fig.5 Three-phase current and d-axis current with without magnetization and demagnetization

3 转子磁场定向瞬时转矩控制系统的构建

基于以上理论分析,构建了如图6所示的无刷直流电机转子磁场定向瞬时转矩控制的整体系统框图。控制系统主要由三相电流参考值计算单元、扇区判断单元、组合电流/组合参考电流计算单元、电流PI控制器、转速计算单元、占空比补偿单元、开关管状态查询表和电压型功率变换电路等组成。

width=453.7,height=157.7

图6 整体系统框图

Fig.6 The block diagram of the control system

系统中,q轴参考电流width=10,height=19给定转矩width=13,height=13.95通过式(9)计算得到;qAqBqC由实际转子位置q 通过数据库查询得到;三相电流参考值计算单元以2.3节介绍的d、q两轴电流对应三相绕组电流计算为依据;在已知比例系数kr的前提下,为减轻微处理器的计算负担,可由width=10,height=19直接计算出三相参考电流width=11,height=17width=11,height=17width=11,height=17,其中width=10,height=19width=11,height=17width=11,height=17width=11,height=17之间的计算系数由kr和式(13)计算得到,可通过预先离线实验建立数据库。

除此之外,为减少三相导通过程中进行脉宽调制的开关管数量,根据图5中每个扇区范围内的三相电流波形特点,取电流大小近似不变相的相电流与剩余两相电流做差得到新的两相组合电流分别记为M相、T相,由于M相、T相电流与A、B、C相电流之间存在固定的对应关系,故可由两相组合电流闭环控制取代三相电流闭环控制。

图6中的组合参考电流计算单元和组合电流计算单元,分别用于获取当前转子位置所处扇区下M相、T相电流的参考值width=13,height=17width=11,height=17和实际值width=13,height=15width=11,height=15;电流PI控制器和占空比补偿单元,分别用于获取M相、T相电流的参考占空比width=20,height=17width=18,height=17和补偿占空比width=21,height=17width=19,height=17,两占空比对应相加,最终得到施加在进行脉宽调制的两相开关管上的占空比width=18,height=17width=16,height=17

图6中各扇区下的组合电流关系及开关管状态查询表,见表1。表中,开关管状态从左到右依次为A、B、C相上、下桥臂开关管的开关状态,0、1和width=18,height=17width=16,height=17分别表示开关管截止、导通和以占空比width=18,height=17width=16,height=17进行脉宽调制。

表1 各扇区电流组合和开关管状态查询表

Tab.1 Current combination and switching table of each sector

扇区S转子位置M相电流T相电流开关管状态 Ⅰ-180°~-120° Ⅱ-120°~-60° Ⅲ-60°~0° Ⅳ0°~60° Ⅴ60°~120° Ⅵ120°~180°

4 仿真与实验验证

为了验证本文所提控制方法的有效性,将一台5对极的三相无刷直流电机用作仿真和实验电机。仿真采用Matlab M-file文件,参数依据实验电机参数给定,搭建的以TMS320F28335为控制核心的无刷直流电机控制实验系统,如图7所示。仿真和实验验证中电机的主要参数见表2,主电路中直流母线电压为300V,开关管的开关频率为20kHz,电流采样和控制周期均为50ms,实验中的电机负载为转矩大小可调节的磁滞测功机。

width=226.8,height=126.5

图7 无刷直流电机控制实验系统

Fig.7 Control experiment system of BLDCM

考虑到实验用无刷直流电机转子永磁体产生的磁场沿定子内表面不可能为理想的梯形波,图8给出了实验用电机实际的三相转子磁链波形及额定转矩负载运行状态下无增磁和去磁时的三相参考电流波形。图9给出了不同转子位置下,实际三相转子磁链对应的各相旋转变换角qAqBqCkA、kB、kC波形。以A相为例,以下给出了由实际转子位置width=9,height=12获取旋转变换角qA的计算过程。根据式(2)可求得qAyA的关系为

表2 无刷直流电机参数

Tab.2 Parameters of BLDCM

参 数数 值 (型号) 电机型号92BL430 额定功率/W400 额定电压/V300 额定电流/A1.5 额定转矩/(N·m)1.27 极对数5 相电阻/W3.05 相电感/mH17 额定转速/(r/min)3 000 编码器分辨率2 048PPR/r

width=227.4,height=135.25

图8 实际的三相转子磁链和无增磁去磁控制所需三相电流波形

Fig.8 Three-phase current waveform of without magnetization and demagnetization under the actual three-phase rotor flux linkage

width=226.8,height=132

图9 旋转变换角qAqBqCkAkBkC

Fig.9 Rotation transform angle qAqBqC and kAkBkC

width=74,height=30 (18)

式中,yA大小仅与转子位置width=9,height=12有关,由式(19)计算得到。

width=131,height=15 (19)

式中,width=29,height=15为第k个转子位置处的A相转子磁链大小;width=31,height=15为随转子位置width=9,height=12变化的A相电动势常数,由事先进行的离线实验获取;width=17,height=12为从第k个转子位置到第k+1个转子位置电机旋转的电角度。考虑到width=9,height=12=-90°位置处的A相转子磁链为0,故以该位置作为A相转子磁链积分的初始位置。

图10、图11分别给出了额定转矩(1.27N×m)负载运行状态下,转速为600r/min时的仿真波形和实验波形。仿真和实验中所用到的控制框图如图6所示,q轴参考电流与三相参考电流之间的计算系数数据库通过预先离线实验得到,各扇区下的开关管状态查询表已在表1中给出。作为仿真结果给出了实际转矩和参考转矩、三相实际相电流和参考相电流及相电压;作为实验结果给出了实际转矩、A相参考电流和实际电流、A相电压和上桥臂开关管驱动信号。仿真和实验中给出的相电压大小均为主电路逆变器三相引出线相对于母线负极的电压。

width=230.15,height=240.1

图10 600r/min下的仿真结果(额定转矩)

Fig.10 Simulation results under 600r/min (rated torque)

由图10和图11可知,转矩脉动抑制效果明显,且不存在换相转矩脉动。仿真和实验中的参考电流波形均为本文提出的无增磁和去磁条件下的参考电流波形。图11中转矩脉动约为7.1%,由A相电压波形和A相上桥臂开关管驱动信号可看出,实验中的调制方式符合表1中给出的调制方式。以A相电流处于最小值时刻为例,结合图8和表1可知,此时无刷直流电机运行在第Ⅴ扇区,C相电流与A相电流相减构成M相电流,B相电流与A相电流相减构成T相电流,在该扇区下,A相的上桥臂开关管关断,B相的上桥臂开关管以占空比width=13.95,height=13.95进行脉宽调制,与图11中给出的开关管驱动信号相吻合。

width=205.55,height=234.9

图11 600r/min下的实验结果(额定转矩)

Fig.11 Experiment results under 600r/min (rated torque)

图12给出了额定转矩(1.27N×m)负载运行状态下,转速为600r/min时的A相参考电流和实际电流、A相电流跟踪误差和相电压波形。由图12可知,电流相对误差最大约为10%,实际电流的跟踪性能良好,表明实验电机实现了无直轴电枢反应控制。

width=226.8,height=111.25

图12 600r/min下的A相参考电流和实际电流、跟踪误差及A相电压(额定转矩)

Fig.12 Phase-A reference current and actual current, current tracking error and phase-A voltage under 600r/min (rated torque)

为了验证本文所提控制方法的动态特性,图13和图14分别给出了转速为1 200r/min运行条件下,转矩由轻载(0.2N×m)阶跃变化到额定转矩(1.27N×m)再阶跃变化到轻载时的仿真波形和实验波形。仿真结果给出了实际转矩和参考转矩、三相实际电流和参考电流及相电压;实验结果给出了实际转矩、A相参考电流和实际电流及A相电压。由图13、图14可知,由轻载阶跃变化到额定转矩时,实际转矩仍可以良好地跟踪参考转矩,转矩并没有出现较大脉动,实验中,稳定状态下的转矩脉动约为10%;在由额定转矩阶跃变化到轻载时,虽然实际转矩稍有滞后,但整体的跟踪控制效果良好。

width=229.8,height=240.7

图13 1 200r/min下的仿真结果(阶跃参考转矩)

Fig.13 Simulation results under 1 200r/min (step reference torque)

width=226.8,height=115.1

图14 1 200r/min下的实验结果(阶跃参考转矩)

Fig.14 Experiment results under 1 200r/min (step reference torque)

5 结论

考虑到传统的3s/2r坐标变换无法得到无刷直流电机的转子磁场定向等效模型,本文提出了一种基于逐相旋转坐标变换的无刷直流电机转子磁场定向瞬时转矩控制技术,并得出以下结论:

1)针对非正弦的三相转子磁链,通过对每相绕组单独进行旋转坐标变换,再以功率守恒为原则进行三相合成,可获得无刷直流电机转子磁场定向dq旋转坐标系下的等效模型。

2)转子磁场定向dq旋转坐标系下,能够做到d、q两轴之间的转矩解耦。其中,q轴电流用于实现恒定瞬时转矩控制;d轴电流可用于实现其他控制目标,如无增磁/去磁的电枢反应等。

3)所提的转子磁场定向瞬时转矩控制技术,不存在换相过程和状态切换;不需要添加额外的硬件电路,具有较强的实用性和灵活性。

参考文献

[1] 边春元, 邢海洋, 李晓霞, 等. 基于速度变化率的无位置传感器无刷直流电机风力发电系统换相误差补偿策略[J]. 电工技术学报, 2021, 36(11): 2374-2382.

Bian Chunyuan, Xing Haiyang, Li Xiaoxia, et al. Compensation strategy for commutation error of sensorless brushless DC motor wind power generation system based on speed change rate[J]. Transactions of China Electrotechnical Society, 2021, 36(11): 2374- 2382.

[2] Jahns T M, Soong W L. Pulsating torque minimi- zation techniques for permanent magnet AC motor drives-a review[J]. IEEE Transactions on Industrial Electronics, 1996, 43(2): 321-330.

[3] de Castro A G, Guazzelli P R U, dos Santos S T C A, et al. Zero sequence power contribution on BLDC motor drives. part I: a theoretical investigation[C]// 2018 13th IEEE International Conference on Industry Applications, Sao Paulo, Brazil, 2018: 1016-1023.

[4] 谭建成. 永磁无刷直流电机技术[M]. 北京: 机械工业出版社, 2011.

[5] Fang Jiancheng, Li Haitao, Han Bangcheng. Torque ripple reduction in BLDC torque motor with nonideal back EMF[J]. IEEE Transactions on Power Elec- tronics, 2012, 27(11): 4630-4637.

[6] 林黄达, 程思为, 王东, 等. 环形绕组无刷直流电机的换向过程分析及建模[J]. 电工技术学报, 2020, 35(8): 1603-1612.

Lin Huangda, Cheng Siwei, Wang Dong, et al. Analytical modeling of the circular winding brushless DC motor under load commutation[J]. Transactions of China Electrotechnical Society, 2020, 35(8): 1603- 1612.

[7] 贾喆武, 张庆湖, 王东. 环形绕组无刷直流电机的四象限运行分析[J]. 电工技术学报, 2020, 35(18): 3821-3829.

Jia Zhewu, Zhang Qinghu, Wang Dong. Analysis on 4-quadrant operation of circular winding brushless DC motor[J]. Transactions of China Electrotechnical Society, 2020, 35(18): 3821-3829.

[8] 林平, 韦鲲, 张仲超. 新型无刷直流电机换相转矩脉动的抑制控制方法[J]. 中国电机工程学报, 2006, 26(3): 153-158.

Lin Ping, Wei Kun, Zhang Zhongchao. A novel control scheme to suppress the commutation torque ripple in BLDCM[J]. Proceedings of the CSEE, 2006, 26(3): 153-158.

[9] 王大方, 卜德明, 朱成, 等. 一种减小无刷直流电机换相转矩脉动的调制方法[J]. 电工技术学报, 2014, 29(5): 160-166.

Wang Dafang, Bu Deming, Zhu Cheng, et al. A modulation method to decrease commutation torque ripple of brushless DC motors[J]. Transactions of China Electrotechnical Society, 2014, 29(5): 160- 166.

[10] 李珍国, 王江浩, 高雪飞, 等. 一种合成电流控制的无刷直流电机转矩脉动抑制系统[J]. 中国电机工程学报, 2015, 35(21): 5592-5599.

Li Zhenguo, Wang Jianghao, Gao Xuefei, et al. Torque ripple reduction systems for brushless DC motors with resultant current control[J]. Proceedings of the CSEE, 2015, 35(21): 5592-5599.

[11] 夏鲲, 朱琳玲, 曾彦能, 等. 基于准Z源网络的永磁无刷直流电机换相转矩脉动抑制方法[J]. 中国电机工程学报, 2015, 35(4): 971-978.

Xia Kun, Zhu Linling, Zeng Yanneng, et al. Researches on the method of suppressing com- mutation torque ripple for brushless DC motors based on a quasi-Z-source net[J]. Proceedings of the CSEE, 2015, 35(4): 971-978.

[12] Chen Wei, Liu Yapeng, Li Xinmin, et al. A novel method of reducing commutation torque ripple for brushless DC motor based on Cuk converter[J]. IEEE Transactions on Power Electronics, 2017, 32(7): 5497-5508.

[13] 姚绪梁, 赵继成, 王景芳, 等. 一种基于辅助升压前端的无刷直流电机换相转矩脉动抑制方法研究[J]. 中国电机工程学报, 2020, 40(9): 3021-3030.

Yao Xuliang, Zhao Jicheng, Wang Jingfang, et al. Research on suppressing commutation torque ripple of brushless DC motor based on an auxiliary step-up front end[J]. Proceedings of the CSEE, 2020, 40(9): 3021-3030.

[14] 曹彦飞, 陆海天, 李新旻, 等. 基于无电感升压拓扑的无刷直流电机电流控制策略[J]. 电工技术学报, 2021, 36(6): 1249-1258.

Cao Yanfei, Lu Haitian, Li Xinmin, et al. Current control strategy of brushless DC motor based on non-inductive Boost topology[J]. Transactions of China Electrotechnical Society, 2021, 36(6): 1249- 1258.

[15] Liu Yong, Zhu Z Q, Howe D. Commutation-torque- ripple minimization in direct-torque-controlled PM brushless DC drives[J]. IEEE Transactions on Industry Applications, 2007, 43(4): 1012-1021.

[16] Shi Tingna, Cao Yanfei, Jiang Guokai, et al. A torque control strategy for torque ripple reduction of brushless DC motor with nonideal back electromotive force[J]. IEEE Transactions on Industrial Electronics, 2017, 64(6): 4423-4433.

[17] 李珍国, 章松发, 周生海, 等. 考虑转矩脉动最小化的无刷直流电机直接转矩控制系统[J]. 电工技术学报, 2014, 29(1): 139-146.

Li Zhenguo, Zhang Songfa, Zhou Shenghai, et al. Direct torque control of brushless DC motor considering torque ripple minimization[J]. Transa- ctions of China Electrotechnical Society, 2014, 29(1): 139-146.

[18] 陆可, 蔡广瀚, 向南辉, 等. 基于电压矢量注入的无刷直流电机换相转矩脉动抑制方法[J]. 中国电机工程学报, 2021, 41(10): 3592-3601.

Lu Ke, Cai Guanghan, Xiang Nanhui, et al. Method for suppressing commutation torque ripple of brush- less DC motor based on voltage vector injection[J]. Proceedings of the CSEE, 2021, 41(10): 3592-3601.

[19] Jiang Weidong, Huang Hui, Wang Jinping, et al. Commutation analysis of brushless DC motor and reducing commutation torque ripple in the two-phase stationary frame[J]. IEEE Transactions on Power Electronics, 2017, 32(6): 4675-4682.

[20] Bertoluzzo M, Buja G, Keshri R K, et al. Sinusoidal versus square-wave current supply of PM brushless DC drives: a convenience analysis[J]. IEEE Transa- ctions on Industrial Electronics, 2015, 62(12): 7339- 7349.

[21] Le-Huy H, Perret R, Feuillet R. Minimization of torque ripple in brushless DC motor drives[J]. IEEE Transactions on Industry Applications, 1986, IA- 22(4): 748-755.

[22] Buja G, Bertoluzzo M, Keshri R K. Torque ripple-free operation of PM BLDC drives with petal-wave current supply[J]. IEEE Transactions on Industrial Electronics, 2015, 62(7): 4034-4043.

[23] Khazaee A, Zarchi H A, Markadeh G R A. Real-time maximum torque per ampere control of brushless DC motor drive with minimum torque ripple[J]. IEEE Transactions on Power Electronics, 2020, 35(2): 1194-1199.

[24] 李珍国, 孙启航, 王鹏磊, 等. 基于转子永磁体磁场定向的无刷直流电机转矩脉动抑制[J]. 电工技术学报, 2020, 35(14): 2987-2996.

Li Zhenguo, Sun Qihang, Wang Penglei, et al. Torque ripple reduction of brushless DC motor based on rotor permanent magnet field orientation[J]. Transactions of China Electrotechnical Society, 2020, 35(14): 2987-2996.

Rotor Field Oriented Instantaneous Torque Control Technology of Brushless DC Motor Based on Per-Phase Rotating Coordinate Transformation

Li Zhenguo Wang Penglei Sun Qihang Jia Yicheng

(Key Lab of Power Electronics for Energy Conservation and Motor Drive of Hebei Province Yanshan University Qinhuangdao 066004 China)

Abstract The equivalent model of a brushless DC motor (BLDCM) with the rotor field orientation has the characteristics that the d-axis rotor flux linkage is constant and the q-axis rotor flux linkage is zero. It isa favorable condition for building high-performance torque control systems. However, the traditional equal-power 3s/2r coordinate transformation cannot obtain the equivalent model of the rotor field orientation of BLDCM. Therefore, a rotor field orientation method based on per-phase coordinate transformation is proposed. That is, per-phase is transformed one by one to obtain the rotor field orientation equivalent model of per-phase winding, and then the three-phase is synthesized based on the principle of power conservation. Since the torque of the equivalent model is proportional to the q-axis current and independent of the d-axis current, the relationship between the d-axis armature reactive magnetomotiveforce, torque and the d- and q-axis currents is derived. Consequently, a BLDCM torque control system with non-magnetization/demagnetization is constructed. Finally, the feasibility and effectiveness of the proposed control technology are verified by Matlab simulation and DSP driving experiment.

Keywords:Brushless DC motor (BLDCM), rotor field orientation, per-phase rotating coordinate transformation, without magnetization and demagnetization, instantaneous torque control

中图分类号:TM351

DOI: 10.19595/j.cnki.1000-6753.tces.211495

国家自然科学基金资助项目(61873226)。

收稿日期 2021-09-23

改稿日期 2021-10-22

作者简介 李珍国 男,1973年生,博士,副教授,研究方向为电力电子与电力传动。

E-mail: lzg@ysu.edu.cn(通信作者)

王鹏磊 男,1995年生,硕士研究生,研究方向为电力电子与电力传动。

E-mail: 1971545348@qq.com

(编辑 郭丽军)