CN113224992B - 一种抑制永磁同步电机单电阻采样电流畸变的方法 - Google Patents

一种抑制永磁同步电机单电阻采样电流畸变的方法 Download PDF

Info

Publication number
CN113224992B
CN113224992B CN202110460647.2A CN202110460647A CN113224992B CN 113224992 B CN113224992 B CN 113224992B CN 202110460647 A CN202110460647 A CN 202110460647A CN 113224992 B CN113224992 B CN 113224992B
Authority
CN
China
Prior art keywords
current
axis
lpf
sampling
motor
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
Application number
CN202110460647.2A
Other languages
English (en)
Other versions
CN113224992A (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.)
Guangdong University of Technology
Original Assignee
Guangdong 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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN202110460647.2A priority Critical patent/CN113224992B/zh
Publication of CN113224992A publication Critical patent/CN113224992A/zh
Application granted granted Critical
Publication of CN113224992B publication Critical patent/CN113224992B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/22Current control, e.g. using a current control loop
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P27/00Arrangements or methods for the control of AC motors characterised by the kind of supply voltage
    • H02P27/04Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage
    • H02P27/06Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters
    • H02P27/08Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation
    • H02P27/12Arrangements or methods for the control of AC motors characterised by the kind of supply voltage using variable-frequency supply voltage, e.g. inverter or converter supply voltage using dc to ac converters or inverters with pulse width modulation pulsing by guiding the flux vector, current vector or voltage vector on a circle or a closed curve, e.g. for direct torque control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2209/00Indexing scheme relating to controlling arrangements characterised by the waveform of the supplied voltage or current
    • H02P2209/11Sinusoidal waveform

Abstract

本发明公开了一种抑制永磁同步电机单电阻采样电流畸变的方法,一方面通过对矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线进行平滑处理,减少实际运行过程中T1、T2不规则变化引起的相电流采样误差,另一方面通过解决由于单电阻电流采样误差导致转子位置角估算误差被放大的问题,进一步抑制相电流波形畸变,使电机相电流尽可能接近正弦波形。

Description

一种抑制永磁同步电机单电阻采样电流畸变的方法
技术领域
本发明涉及电机控制的技术领域,尤其涉及到一种抑制永磁同步电机单电阻采样电流畸变的方法。
背景技术
由于永磁同步电机具有结构紧凑、功率密度和电机效率高且调速范围宽等优点,目前已被广泛应用于家电、风机泵类产品等领域。小功率永磁同步电机及其驱动系统产品要求较高的性价比,在位置检测方面通常所采用的HALL位置传感器或光电编码位置传感器逐渐被无位置传感器控制算法所替代;在电流采样方面,常用的电流采样方法是采用两个电流传感器测量相电流,这种方法成本较高,在整机硬件成本中占有不可忽视的比例,相电流的电阻采样方案可以简化系统结构,降低硬件电路成本,目前较常用的方案是双电阻采样方案。电机相电流单电阻采样方案可以进一步简化硬件电路布局、降低硬件电路成本,近几年来无位置传感器、单电阻采样永磁同步电机矢量控制方案越来越受到生产厂家的重视。
由于逆变器的非线性以及永磁同步电机气隙磁场畸变,在运行过程中永磁同步电机定子电流中含有大量的谐波分量,导致电流波形发生畸变,进一步引起电磁噪音,同时也影响电机的效率。李毅拓,陆海峰等在《电机工程学报》(2014,NO.34)提出了基于谐振调节器永磁同步电机电流谐波抑制方法;李生民等在《电机与控制应用》(2019,NO.46)提出了基于闭环电流平均值的谐波电压注入法来降低电流波形的畸变率及5、7次谐波含量。这些方法属于传统电流畸变抑制方法,主要思路是采用谐振调节器或低通滤波器提取并消除相电流中的5、7次谐波分量。
在减小单电阻采样永磁同步电机电流波形畸变的研究方面,黄科元等人在《电力系统及其自动化学报》(VOL.30,NO.9)上提出了一种单电阻采样永磁同步电机相电流重构策略,在低调制不可观测区和中调制不可观测区插入测量脉冲,在高调制不可观测区采用电压矢量近似的方法,但这种方法存在算法执行时间长、实施过程复杂的缺点;Saritha B等人在《IEEE Trans on Industrial Electronics》(VOL.54,NO.5)针对不可观测区采用正弦曲线拟合观测器,使估计的电流趋近参考正弦三相电流,但这种方法依赖电机参数,且不能解决低速误差问题。
除逆变器和气隙磁场的非线性原因外,无位置传感器、单电阻采样永磁同步电机控制系统电流畸变的原因还包括:①通过无位置传感器算法计算出的转子位置角与转子实际位置角之间存在误差;②单电阻采样需要对驱动系统母线电流进行三相电流重构,电流重构过程中存在盲区,即不可观测区,从而产生电流采样误差。现有技术同时针对这两种原因进行电流畸变抑制的策略目前还没有相关的文献记载。
发明内容
本发明的目的在于克服现有技术的不足,提供一种抑制永磁同步电机单电阻采样电流畸变的方法,通过同时减小转子位置角估算和电流采样两方面误差,来抑制永磁同步电机控制系统相电流波形的畸变,使相电流尽可能接近正弦波。
为实现上述目的,本发明所提供的技术方案为:
一种抑制永磁同步电机单电阻采样电流畸变的方法,应用在单电阻采样永磁同步电机无位置传感器控制系统中,该控制系统包括电流采样、转子位置估算、Clarke和PARK变换、最大转矩电流比控制MTPA、速度环、dq轴电流环、PARK逆变换、SVPWM计算、三相PWM逆变器在内的单元;
针对该系统,一方面通过对矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线进行平滑处理,减少实际运行过程中T1、T2不规则变化引起的相电流采样误差,另一方面通过解决由于单电阻电流采样误差导致转子位置角估算误差被放大的问题,进一步抑制相电流波形畸变,使电机相电流尽可能接近正弦波形,具体包括如下步骤:
S1、当永磁同步电机在设定转速nset稳定运行时,计算矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线;
S2、根据T1、T2变化曲线,确定每一次矢量变频控制迭代运算过程中T1、T2变化最大幅值,在运行过程中对T1、T2的变化进行限幅处理,减小T1、T2波形的毛刺;
S3、在稳定运行时,对d轴电流id进行低通滤波得到id_LPF,对q轴电流id进行低通滤波得到iq_LPF
S4、根据id_LPF、iq_LPF对d、q轴电流进行限幅处理得到
Figure BDA0003042111410000031
Figure BDA0003042111410000032
S5、根据
Figure BDA0003042111410000033
Figure BDA0003042111410000034
进行转子位置角和转速的估算,减小单电阻电流采样误差导致的转子位置角估算误差。
进一步地,所述步骤S1中,设永磁同步电机矢量控制七段SVPWM对应的空间电压矢量,有八种开关工作状态,包括六个非零电压矢量V1~V6和两个零电压矢量V0、V7,其将电压空间平面分成六个扇区,每个扇区中任意电压矢量Vref均可由对应扇区的两个非零电压矢量和零矢量共同合成产生,则第一扇区中对应的两个非零电压矢量作用时间T1、T2的计算公式如下:
Figure BDA0003042111410000035
Ts=T0+T1+T2; (2)
Figure BDA0003042111410000036
Figure BDA0003042111410000037
其中,Ts为一个矢量控制运算周期,第一扇区对应的非零电压矢量为V4、V6,T0为零电压矢量V0、V7的作用时间总和,Vdc为母线电压,θref为电压矢量Vref与零电压矢量V4之间的夹角。
进一步地,所述步骤S2的具体过程如下:
S2-1、根据设定转速nset计算出每个相电流的电周期
Figure BDA0003042111410000041
其中pn为电机的极对数;每个相电流周期包含六个扇区,进一步计算出每一个扇区控制迭代运算的次数
Figure BDA0003042111410000042
其中Ts为一个矢量控制运算周期;
S2-2、计算每一次矢量变频控制迭代运算T1、T2变化的最大幅值
Figure BDA0003042111410000043
Figure BDA0003042111410000044
式(5)和式(6)中,T1_max、T2_max分别为T1、T2的最大值,T1_min、T2_min分别为T1、T2的最小值,ξ为放大系数;
S2-3、根据T1、T2处于上升期或下降期的情况,对T1、T2进行限幅处理:
若T1处于上升期:
T1(n)<T1(n-1)+ΔT1_max; (7)
若T1处于下降期:
T1(n)>T1(n-1)-ΔT1_max; (8)
若T2处于上升期:
T2(n)<T2(n-1)+ΔT2_max; (9)
若T2处于下降期:
T2(n)>T2(n-1)-ΔT2_max。 (10)
进一步地,所述步骤S3中,采用的低通滤波的离散方程为:
id_LPF(n+1)=id_LPF(n)+KLPF{id(n+1)-id_LPF(n)}; (11)
iq_LPF(n+1)=iq_LPF(n)+KLPF{iq(n+1)-iq_LPF(n)}; (12)
式(11)和式(12)中,id(n+1)、iq(n+1)分别为d轴、q轴电流在采样点n+1处的瞬时值,id_LPF(n)、iq_LPF(n)分别为d轴、q轴电流在采样点n处的低通滤波值,id_LPF(n+1)、iq_LPF(n+1)分别为d轴、q轴电流在采样点n+1处的低通滤波值,KLPF为低通滤波系数,d轴和q轴采用相同的低通滤波系数。
进一步地,所述步骤S4根据id_LPF、iq_LPF对d、q轴电流进行限幅处理得到
Figure BDA0003042111410000051
Figure BDA0003042111410000052
具体为:
Figure BDA0003042111410000053
Figure BDA0003042111410000054
Figure BDA0003042111410000055
Figure BDA0003042111410000056
Figure BDA0003042111410000057
Figure BDA0003042111410000058
Figure BDA0003042111410000059
Figure BDA00030421114100000510
式(13)至式(16)中,
Figure BDA00030421114100000511
分别为d、q轴电流在采样点n+1处的限幅值,Δid_max、Δiq_max为每一次迭代运算d、q轴电流变化的最大幅值。
进一步地,所述步骤S5中,进行转子位置角估算的过程如下:
S5-1、建立dq坐标系下的定子电压方程:
Figure BDA0003042111410000061
式(17)中,ud、uq分别为定子绕组d、q轴电压;id、iq分别为d、q轴电流;Rs为定子电阻;Ld、Lq分别为d、q轴电感;e为电机反电动势;ω代表dq坐标系的旋转角速度;p为微分算子,p=d/dt;
S5-2、建立γδ估算坐标系下的定子电压方程:
Figure BDA0003042111410000062
式(18)中,uγ、uδ分别为γ、δ轴定子电压分量,iγ、iδ分别为γ、δ轴电流分量,ωM代表γδ估算坐标系的旋转角速度;Δθ为γδ估算坐标系与dq坐标系的夹角,即位置角估算误差;
S5-3、求出采样点n+1处电机的实际电流iγ(n+1)、iδ(n+1):
Figure BDA0003042111410000063
式(19)中,T为离散点的采样时间,iγ(n)、iδ(n)为采样点n处电机的实际电流,uγ(n)、uδ(n)为采样点n处电机的实际电压;
S5-4、求出采样点n+1处估算电流i(n+1)、i(n+1):
Figure BDA0003042111410000064
S5-5、从式(19)和式(20)求出采样点n+1处估算电流误差Δiγ(n+1)、Δiδ(n+1):
Figure BDA0003042111410000065
S5-6、根据式(21)求出的电流误差Δiδ推算出电机反电动势:
eM(n+1)=eM(n)-KδΔiδ(n+1) (22)
式(22)中,eM(n+1)、eM(n)分别为采样点n+1、采样点n处的电机反电动势,Δiδ(n+1)为采样点n+1处的δ轴电流误差,Kδ为反电势估算系数;
S5-7、根据式(22)推算的电机反电动势推算出电机转子位置角:
Figure BDA0003042111410000071
式(23)中,θM(n+1)、θM(n)分别为采样点n+1、采样点n处的电机转子位置角,Δiγ(n+1)为采样点n+1处的γ轴电流误差,KE为电机反电势系数,Kθ为转子位置角补偿系数。
与现有技术相比,本方案原理及优点如下:
本方案一方面通过对矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线进行平滑处理,减少实际运行过程中T1、T2不规则变化引起的相电流采样误差,另一方面通过解决由于单电阻电流采样误差导致转子位置角估算误差被放大的问题,进一步抑制相电流波形畸变,使电机相电流尽可能接近正弦波形。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的服务作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为单电阻采样永磁同步电机三相电流重构示意图;
图2为永磁同步电机空间电压矢量图;
图3为γδ假定坐标系和dq坐标系;
图4为永磁同步电机系统控制框图;
图5为限幅处理前后T1、T2曲线波形对比图;
图6为优化前对d、q轴电流进行低通滤波的波形图;
图7为d、q轴电流限幅处理前后位置角估算误差对比图;
图8为电流畸变抑制前后的波形对比图。
具体实施方式
下面结合具体实施例对本发明作进一步说明:
本实施例所采用的永磁同步电机为一台应用于家用空气净化器的外转子风扇电机,其中永磁同步电机的参数为:额定电压DC 24V,最小运行转速nset_min=500转/分(rpm),最高运行转速nset_max=3000转/分,极对数pn=2,定子电阻Rs=1.6Ω,定子直轴电感Ld=1.0mH,交轴电感Lq=1.2mH,反电势系数ke=4.2V/krpm。根据系统硬件参数及功率管死区时间设定,单电阻采样的最小采样时间Tmin设为4μs。
而对应的单电阻采样永磁同步电机系统采用无位置传感器矢量控制,如图4为系统控制框图,包括电流采样、转子位置估算、Clarke和PARK变换、最大转矩电流比控制(MTPA)、速度环、dq轴电流环、PARK逆变换、SVPWM计算、三相PWM逆变器等单元。图4中转子位置估算模块的输入电流由原先的id、iq,改成了通过限幅处理之后的
Figure BDA0003042111410000081
针对该系统,本实施例所述的一种抑制永磁同步电机单电阻采样电流畸变的方法,具体包括如下步骤:
S1、当永磁同步电机在设定转速nset稳定运行时,计算矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线;
单电阻采样永磁同步电机三相电流重构的基本原理就是利用一个PWM周期内在不同的时刻采样母线电流,通过相电流重构得到各个相电流,如图1为单电阻采样永磁同步电机三相电流重构示意图。电机控制器采用SVPWM调制方式控制,如图2所示为永磁同步电机矢量控制七段SVPWM对应的空间电压矢量图,有八种开关工作状态,包括六个非零电压矢量V1~V6和两个零电压矢量V0、V7,其将电压空间平面分成六个扇区,每个扇区中任意电压矢量Vref均可由对应扇区的两个非零电压矢量和零矢量共同合成产生,分析电流重构不可观测区首先需要计算每个扇区对应的非零电压矢量的作用时间T1、T2,以第一扇区为例,可以按如下的公式计算T1、T2
Figure BDA0003042111410000091
Ts=T0+T1+T2; (2)
Figure BDA0003042111410000092
Figure BDA0003042111410000093
其中,Ts为一个矢量控制运算周期,第一扇区对应的非零电压矢量为V4、V6,T0为零电压矢量V0、V7的作用时间总和,Vdc为母线电压,θref为电压矢量Vref与零电压矢量V4之间的夹角。其他五个扇区T1、T2的计算可以按第一扇区类似的方法获得。
S2、根据T1、T2变化曲线,确定每一次矢量变频控制迭代运算过程中T1、T2变化最大幅值,在运行过程中对T1、T2的变化进行限幅处理,减小T1、T2波形的毛刺,消除采样误差的影响,具体过程如下:
S2-1、根据设定转速nset计算出每个相电流的电周期
Figure BDA0003042111410000094
其中pn为电机的极对数;每个相电流周期包含六个扇区,进一步计算出每一个扇区控制迭代运算的次数
Figure BDA0003042111410000095
其中Ts为一个矢量控制运算周期;
S2-2、计算每一次矢量变频控制迭代运算T1、T2变化的最大幅值
Figure BDA0003042111410000101
Figure BDA0003042111410000102
式(5)和式(6)中,T1_max、T2_max分别为T1、T2的最大值,T1_min、T2_min分别为T1、T2的最小值,ξ为放大系数;
S2-3、根据T1、T2处于上升期或下降期的情况,对T1、T2进行限幅处理:
若T1处于上升期:
T1(n)<T1(n-1)+ΔT1_max; (7)
若T1处于下降期:
T1(n)>T1(n-1)-ΔT1_max; (8)
若T2处于上升期:
T2(n)<T2(n-1)+ΔT2_max; (9)
若T2处于下降期:
T2(n)>T2(n-1)-ΔT2_max。 (10)
S3、在稳定运行时,对d轴电流id进行低通滤波得到id_LPF,对q轴电流id进行低通滤波得到iq_LPF,采用的低通滤波的离散方程为:
id_LPF(n+1)=id_LPF(n)+KLPF{id(n+1)-id_LPF(n)}; (11)
iq_LPF(n+1)=iq_LPF(n)+KLPF{iq(n+1)-iq_LPF(n)}; (12)
式(11)和式(12)中,id(n+1)、iq(n+1)分别为d轴、q轴电流在采样点n+1处的瞬时值,id_LPF(n)、iq_LPF(n)分别为d轴、q轴电流在采样点n处的低通滤波值,id_LPF(n+1)、iq_LPF(n+1)分别为d轴、q轴电流在采样点n+1处的低通滤波值,KLPF为低通滤波系数,d轴和q轴采用相同的低通滤波系数。
S4、根据id_LPF、iq_LPF对d、q轴电流进行限幅处理得到
Figure BDA0003042111410000111
Figure BDA0003042111410000112
具体为:
Figure BDA0003042111410000113
Figure BDA0003042111410000114
Figure BDA0003042111410000115
Figure BDA0003042111410000116
Figure BDA0003042111410000117
Figure BDA0003042111410000118
Figure BDA0003042111410000119
Figure BDA00030421114100001110
式(13)至式(16)中,
Figure BDA00030421114100001111
分别为d、q轴电流在采样点n+1处的限幅值,Δid_max、Δiq_max为每一次迭代运算d、q轴电流变化的最大幅值。当系统受到负载扰动或输入电压突变时,为了满足系统动态响应能力要求,可适当增大Δid_max、Δiq_max的取值。
S5、根据
Figure BDA00030421114100001112
Figure BDA00030421114100001113
进行转子位置角和转速的估算,减小单电阻电流采样误差导致的转子位置角估算误差。通过下面的推导可以解释对d、q轴电流进行限幅处理减小转子位置角估算误差的有效性,如图3为永磁同步电机无位置传感器矢量控制所采用的dq坐标系和为了进行转子位置估算而建立的γδ假定坐标系,本实施例进行转子位置角估算的算法如下:
S5-1、建立dq坐标系下的定子电压方程:
Figure BDA00030421114100001114
式(17)中,ud、uq分别为定子绕组d、q轴电压;id、iq分别为d、q轴电流;Rs为定子电阻;Ld、Lq分别为d、q轴电感;e为电机反电动势;ω代表dq坐标系的旋转角速度;p为微分算子,p=d/dt;
S5-2、建立γδ估算坐标系下的定子电压方程:
Figure BDA0003042111410000121
式(18)中,uγ、uδ分别为γ、δ轴定子电压分量,iγ、iδ分别为γ、δ轴电流分量,ωM代表γδ估算坐标系的旋转角速度;Δθ为γδ估算坐标系与dq坐标系的夹角,即位置角估算误差;
S5-3、求出采样点n+1处电机的实际电流iγ(n+1)、iδ(n+1):
Figure BDA0003042111410000122
式(19)中,T为离散点的采样时间,iγ(n)、iδ(n)为采样点n处电机的实际电流,uγ(n)、uδ(n)为采样点n处电机的实际电压;
S5-4、求出采样点n+1处估算电流i(n+1)、i(n+1):
Figure BDA0003042111410000123
S5-5、从式(19)和式(20)求出采样点n+1处估算电流误差Δiγ(n+1)、Δiδ(n+1):
Figure BDA0003042111410000124
S5-6、根据式(21)求出的电流误差Δiδ推算出电机反电动势:
eM(n+1)=eM(n)-KδΔiδ(n+1) (22)
式(22)中,eM(n+1)、eM(n)分别为采样点n+1、采样点n处的电机反电动势,Δiδ(n+1)为采样点n+1处的δ轴电流误差,Kδ为反电势估算系数;
S5-7、根据式(22)推算的电机反电动势推算出电机转子位置角:
Figure BDA0003042111410000131
式(23)中,θM(n+1)、θM(n)分别为采样点n+1、采样点n处的电机转子位置角,Δiγ(n+1)为采样点n+1处的γ轴电流误差,KE为电机反电势系数,Kθ为转子位置角补偿系数。
从公式(21)可见,位置角估算误差Δθ与γ轴电流误差Δiγ(n+1)是成正比的,而从图3可看出,γ轴电流分量非常接近d轴电流分量,因此在本实施例的无位置传感器矢量控制系统中,由于单电阻采样的相电流产生的误差,通过CLARK和PARK变换会产生d、q轴电流误差,而其中的d轴电流分量误差引起γ轴电流误差Δiγ(n+1),γ轴电流误差Δiγ(n+1)通过公式(21)将直接产生位置角估算误差Δθ。而位置角估算误差Δθ在下一次迭代运算的CLARK和PARK变换中,和单电阻采样的相电流误差一起,将使d、q轴电流误差成分变得复杂。本实施例对d、q轴电流进行限幅处理以减小转子位置角估算误差的方法具有合理性。
如图5为实施方案中限幅处理前后T1、T2曲线波形对比,图中T1、T2的变化范围都是(0~20)ns,图5(a)进行限幅处理前,T1、T2曲线波形存在大量毛刺,这时单电阻电流重构既可能进入扇区过渡不可观测区域,也可能进入低调制不可观测区域,为了获得足够的采样时间,控制算法需要进行频繁的脉冲移相,这导致PWM输出波形极不对称,增大了采样误差,也增加了相电流谐波成分,引起电流畸变;图5(b)进行限幅处理前,T1、T2曲线波形毛刺大大减少,这里取ΔT1_min=ΔT2_min=0.04ns、ξ=1.3。通过限幅处理,单电阻电流重构稳定运行情况下不进入低调制不可观测区域,而只进入扇区过渡不可观测区域,这样可以大大减小脉冲移相的次数,减小电流采样误差,抑制电流畸变。
如图6为优化前对d、q轴电流进行低通滤波的波形,d轴和q轴电流都存在丰富的谐波,对d轴电流id进行低通滤波得到id_LPF在x轴附近波动,对q轴电流id进行低通滤波得到iq_LPF在0.42A附近波动,这里取Δid_max=Δiq_max=0.02A。如图7为对d、q轴电流限幅处理前后位置角估算误差对比,d、q轴电流限幅处理前位置估算误差约为5°,限幅处理后位置估算误差约为1.5°,对d、q轴电流进行限幅处理可以大大减小转子位置角估算误差。
图8为本实施例方案应用前后,电机在nset=500rpm处相电流波形的对比,从图8(a)可见优化前相电流波形存在畸变,畸变的电流波形关于X轴不对称,并且可以看出明显不是因为5、7次谐波分量过高而引起的畸变,实验表明采用通常的消除5、7次谐波的方法并不能抑制本实施例中的电流畸变。进一步分析发现是由于脉冲移相次数频繁,图8(a)中的相电流在波峰处电流不光滑,存在大量毛刺。本实施例的电流畸变抑制策略应用后的相电流波形如图8(b)所示,从图中可见电流波形的畸变得到明显抑制,同时由于脉冲移相次数大大减少,相电流波形中的毛刺也大大减少,电流波峰处比较光滑。
图8表明本实施例所提出的电流畸变抑制策略是有效的。
以上所述之实施例子只为本发明之较佳实施例,并非以此限制本发明的实施范围,故凡依本发明之形状、原理所作的变化,均应涵盖在本发明的保护范围内。

Claims (5)

1.一种抑制永磁同步电机单电阻采样电流畸变的方法,其特征在于,应用在单电阻采样永磁同步电机无位置传感器控制系统中,该控制系统包括电流采样、转子位置估算、Clarke和PARK变换、最大转矩电流比控制MTPA、速度环、dq轴电流环、PARK逆变换、SVPWM计算、三相PWM逆变器在内的单元;
针对该系统,一方面通过对矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线进行平滑处理,减少实际运行过程中T1、T2不规则变化引起的相电流采样误差,另一方面通过解决由于单电阻电流采样误差导致转子位置角估算误差被放大的问题,进一步抑制相电流波形畸变,使电机相电流尽可能接近正弦波形,具体包括如下步骤:
S1、当永磁同步电机在设定转速nset稳定运行时,计算矢量变频控制每一个扇区对应的两个非零电压矢量作用时间T1、T2的变化曲线;
S2、根据T1、T2变化曲线,确定每一次矢量变频控制迭代运算过程中T1、T2变化最大幅值,在运行过程中对T1、T2的变化进行限幅处理,减小T1、T2波形的毛刺;
S3、在稳定运行时,对d轴电流id进行低通滤波得到id_LPF,对q轴电流id进行低通滤波得到iq_LPF
S4、根据id_LPF、iq_LPF对d、q轴电流进行限幅处理得到
Figure FDA0003496518930000011
Figure FDA0003496518930000012
S5、根据
Figure FDA0003496518930000013
Figure FDA0003496518930000014
进行转子位置角和转速的估算,减小单电阻电流采样误差导致的转子位置角估算误差;
所述步骤S4根据id_LPF、iq_LPF对d、q轴电流进行限幅处理得到
Figure FDA0003496518930000015
Figure FDA0003496518930000016
具体为:
Figure FDA0003496518930000017
Figure FDA0003496518930000021
Figure FDA0003496518930000022
Figure FDA0003496518930000023
Figure FDA0003496518930000024
Figure FDA0003496518930000025
Figure FDA0003496518930000026
Figure FDA0003496518930000027
式(13)至式(16)中,
Figure FDA0003496518930000028
分别为d、q轴电流在采样点n+1处的限幅值,Δid_max、Δiq_max为每一次迭代运算d、q轴电流变化的最大幅值,id_LPF(n+1)、iq_LPF(n+1)分别为d轴、q轴电流在采样点n+1处的低通滤波值。
2.根据权利要求1所述的一种抑制永磁同步电机单电阻采样电流畸变的方法,其特征在于,所述步骤S1中,设永磁同步电机矢量控制七段SVPWM对应的空间电压矢量,有八种开关工作状态,包括六个非零电压矢量V1~V6和两个零电压矢量V0、V7,其将电压空间平面分成六个扇区,每个扇区中任意电压矢量Vref均可由对应扇区的两个非零电压矢量和零矢量共同合成产生,则第一扇区中对应的两个非零电压矢量作用时间T1、T2的计算公式如下:
Figure FDA0003496518930000029
Ts=T0+T1+T2 (2)
Figure FDA00034965189300000210
Figure FDA00034965189300000211
其中,Ts为一个矢量控制运算周期,第一扇区对应的非零电压矢量为V4、V6,T0为零电压矢量V0、V7的作用时间总和,Vdc为母线电压,θref为电压矢量Vref与非零电压矢量V4之间的夹角。
3.根据权利要求1所述的一种抑制永磁同步电机单电阻采样电流畸变的方法,其特征在于,所述步骤S2的具体过程如下:
S2-1、根据设定转速nset计算出每个相电流的电周期
Figure FDA0003496518930000031
其中pn为电机的极对数;每个相电流周期包含六个扇区,进一步计算出每一个扇区控制迭代运算的次数
Figure FDA0003496518930000032
其中Ts为一个矢量控制运算周期;
S2-2、计算每一次矢量变频控制迭代运算T1、T2变化的最大幅值
Figure FDA0003496518930000033
Figure FDA0003496518930000034
式(5)和式(6)中,T1_max、T2_max分别为T1、T2的最大值,T1_min、T2_min分别为T1、T2的最小值,ξ为放大系数;
S2-3、根据T1、T2处于上升期或下降期的情况,对T1、T2进行限幅处理:
若T1处于上升期:
T1(n)<T1(n-1)+ΔT1_max (7)
若T1处于下降期:
T1(n)>T1(n-1)-ΔT1_max (8)
若T2处于上升期:
T2(n)<T2(n-1)+ΔT2_max (9)
若T2处于下降期:
T2(n)>T2(n-1)-ΔT2_max (10)。
4.根据权利要求1所述的一种抑制永磁同步电机单电阻采样电流畸变的方法,其特征在于,所述步骤S3中,采用的低通滤波的离散方程为:
id_LPF(n+1)=id_LPF(n)+KLPF{id(n+1)-id_LPF(n)} (11)
iq_LPF(n+1)=iq_LPF(n)+KLPF{iq(n+1)-iq_LPF(n)} (12)
式(11)和式(12)中,id(n+1)、iq(n+1)分别为d轴、q轴电流在采样点n+1处的瞬时值,id_LPF(n)、iq_LPF(n)分别为d轴、q轴电流在采样点n处的低通滤波值,id_LPF(n+1)、iq_LPF(n+1)分别为d轴、q轴电流在采样点n+1处的低通滤波值,KLPF为低通滤波系数,d轴和q轴采用相同的低通滤波系数。
5.根据权利要求1所述的一种抑制永磁同步电机单电阻采样电流畸变的方法,其特征在于,所述步骤S5中,进行转子位置角估算的过程如下:
S5-1、建立dq坐标系下的定子电压方程:
Figure FDA0003496518930000041
式(17)中,ud、uq分别为定子绕组d、q轴电压;id、iq分别为d、q轴电流;Rs为定子电阻;Ld、Lq分别为d、q轴电感;e为电机反电动势;ω代表dq坐标系的旋转角速度;p为微分算子,p=d/dt;
S5-2、建立γδ估算坐标系下的定子电压方程:
Figure FDA0003496518930000042
式(18)中,uγ、uδ分别为γ、δ轴定子电压分量,iγ、iδ分别为γ、δ轴电流分量,ωM代表γδ估算坐标系的旋转角速度;Δθ为γδ估算坐标系与dq坐标系的夹角,即位置角估算误差;
S5-3、求出采样点n+1处电机的实际电流iγ(n+1)、iδ(n+1):
Figure FDA0003496518930000051
式(19)中,T为离散点的采样时间,iγ(n)、iδ(n)为采样点n处电机的实际电流,uγ(n)、uδ(n)为采样点n处电机的实际电压;
S5-4、求出采样点n+1处估算电流i(n+1)、i(n+1):
Figure FDA0003496518930000052
S5-5、从式(19)和式(20)求出采样点n+1处估算电流误差Δiγ(n+1)、Δiδ(n+1):
Figure FDA0003496518930000053
S5-6、根据式(21)求出的电流误差Δiδ推算出电机反电动势:
eM(n+1)=eM(n)-KδΔiδ(n+1) (22)
式(22)中,eM(n+1)、eM(n)分别为采样点n+1、采样点n处的电机反电动势,Δiδ(n+1)为采样点n+1处的δ轴电流误差,Kδ为反电势估算系数;
S5-7、根据式(22)推算的电机反电动势推算出电机转子位置角:
Figure FDA0003496518930000054
式(23)中,θM(n+1)、θM(n)分别为采样点n+1、采样点n处的电机转子位置角,Δiγ(n+1)为采样点n+1处的γ轴电流误差,KE为电机反电势系数,Kθ为转子位置角补偿系数。
CN202110460647.2A 2021-04-27 2021-04-27 一种抑制永磁同步电机单电阻采样电流畸变的方法 Active CN113224992B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110460647.2A CN113224992B (zh) 2021-04-27 2021-04-27 一种抑制永磁同步电机单电阻采样电流畸变的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110460647.2A CN113224992B (zh) 2021-04-27 2021-04-27 一种抑制永磁同步电机单电阻采样电流畸变的方法

Publications (2)

Publication Number Publication Date
CN113224992A CN113224992A (zh) 2021-08-06
CN113224992B true CN113224992B (zh) 2022-04-22

Family

ID=77089261

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110460647.2A Active CN113224992B (zh) 2021-04-27 2021-04-27 一种抑制永磁同步电机单电阻采样电流畸变的方法

Country Status (1)

Country Link
CN (1) CN113224992B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104682811A (zh) * 2013-05-12 2015-06-03 英飞凌科技股份有限公司 同步电动机优化控制
CN110323986A (zh) * 2019-06-04 2019-10-11 广东工业大学 一种永磁同步电机转子位置角估算方法
CN111835258A (zh) * 2020-07-27 2020-10-27 广东希塔变频技术有限公司 电机驱动控制方法、控制装置和驻车空调控制器
CN112688611A (zh) * 2020-12-14 2021-04-20 泛仕达机电股份有限公司 一种抑制单电阻采样永磁同步电机低速噪音的方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4988374B2 (ja) * 2007-02-15 2012-08-01 三洋電機株式会社 モータ制御装置
JP2011125107A (ja) * 2009-12-09 2011-06-23 Sanyo Electric Co Ltd モータ制御装置、モータ駆動システム及びインバータ制御装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104682811A (zh) * 2013-05-12 2015-06-03 英飞凌科技股份有限公司 同步电动机优化控制
CN110323986A (zh) * 2019-06-04 2019-10-11 广东工业大学 一种永磁同步电机转子位置角估算方法
CN111835258A (zh) * 2020-07-27 2020-10-27 广东希塔变频技术有限公司 电机驱动控制方法、控制装置和驻车空调控制器
CN112688611A (zh) * 2020-12-14 2021-04-20 泛仕达机电股份有限公司 一种抑制单电阻采样永磁同步电机低速噪音的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
永磁同步电机无位置传感器矢量控制的实验研究;闫茂德 等;《微特电机》;20160831;第44卷(第8期);第99-104页 *

Also Published As

Publication number Publication date
CN113224992A (zh) 2021-08-06

Similar Documents

Publication Publication Date Title
JP3582505B2 (ja) モーター制御装置
Jevremovic et al. Speed-sensorless control of induction motor based on reactive power with rotor time constant identification
JP2003061386A (ja) 同期電動機駆動システム
CN112688611B (zh) 一种抑制单电阻采样永磁同步电机低速噪音的方法
CN107508521B (zh) 一种永磁同步电机的无速度传感器控制方法和系统
Chen et al. Design and implementation of a novel high-performance sensorless control system for interior permanent magnet synchronous motors
CN110022106A (zh) 一种基于高频信号注入的永磁同步电机无位置传感器控制方法
CN112671298B (zh) 改进的用于永磁同步电机控制的pll无感控制算法
CN111181458A (zh) 基于扩展卡尔曼滤波器的表贴式永磁同步电机转子磁链观测方法
CN111769779A (zh) 基于改进型Luenberger观测器的PMSM直接转矩控制方法
CN111800056A (zh) 一种基于新型开关表的永磁同步电机三矢量模型预测转矩控制方法
JP3586078B2 (ja) 電力変換装置
CN112688614A (zh) 一种新型同步磁阻电机转速估计方法
CN112039386A (zh) 一种基于模糊准比例谐振的永磁同步电机转矩脉动抑制方法
CN113224992B (zh) 一种抑制永磁同步电机单电阻采样电流畸变的方法
Chen et al. Implementation of a predictive controller for a sensorless interior permanent-magnet synchronous motor drive system
Wang et al. Comparative study of low-pass filter and phase-locked loop type speed filters for sensorless control of AC drives
Kawamura et al. Position Sensorless Adaptive Positioning Servo System with Simplified Differential Calculation and High-Frequency Voltage Injection Strategy Considering Acoustic Noise Suppression
CN115149866A (zh) 一种永磁同步电机全速域无位置传感器矢量控制方法
CN114499327A (zh) 永磁同步电机磁链补偿无位置传感器控制方法及控制系统
CN113992087A (zh) 一种电机全速域无传感位置估计与控制方法及系统
Zhang et al. Fundamental PWM Excitation Based Low-Speed Sensorless Control Method for PMSM Drives
CN111654223A (zh) 一种混合励磁磁通切换电机无位置传感器复合控制方法
Mitsui et al. Motor Current Reconstruction Method Using Single Shunt Resistance by High-Frequency Voltage Injection
Zhai et al. An Improved Position Observer for Sensorless FOC 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