CN112039390B - 基于负载转矩观测的永磁同步电机滑模控制方法 - Google Patents

基于负载转矩观测的永磁同步电机滑模控制方法 Download PDF

Info

Publication number
CN112039390B
CN112039390B CN202010918601.6A CN202010918601A CN112039390B CN 112039390 B CN112039390 B CN 112039390B CN 202010918601 A CN202010918601 A CN 202010918601A CN 112039390 B CN112039390 B CN 112039390B
Authority
CN
China
Prior art keywords
load torque
value
current
torque
sliding mode
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
CN202010918601.6A
Other languages
English (en)
Other versions
CN112039390A (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.)
Ningbo Fenghua Wanyun Motor Co ltd
Shenzhen Dragon Totem Technology Achievement Transformation Co ltd
Original Assignee
Hunan 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 Hunan University of Technology filed Critical Hunan University of Technology
Priority to CN202010918601.6A priority Critical patent/CN112039390B/zh
Publication of CN112039390A publication Critical patent/CN112039390A/zh
Application granted granted Critical
Publication of CN112039390B publication Critical patent/CN112039390B/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/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • H02P21/0007Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control using sliding mode 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
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/05Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation specially adapted for damping motor oscillations, e.g. for reducing hunting
    • 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/13Observer control, e.g. using Luenberger observers or Kalman filters
    • 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/20Estimation of torque
    • 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
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022Synchronous motors
    • 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
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

本发明公开了一种基于负载转矩观测的永磁同步电机滑模控制方法,速度环采用滑模控制方式,负载转矩观测采用反馈增益依据负载转矩给定值的变化量和负载转矩观测值的变化量进行自动调整的算法,能在系统参数发生变化或者是负载发生扰动,导致负载转矩给定值发生变化或/和负载转矩观测值发生变化时,快速降低负载转矩的观测误差并将负载转矩观测值前馈补偿至电流调节器中,有效的削弱了系统的抖振,且动态响应速度快,鲁棒性高,提高了永磁同步电机的控制精度。

Description

基于负载转矩观测的永磁同步电机滑模控制方法
技术领域
本发明涉及永磁同步电机技术领域,更具体地,尤其是涉及一种基于负载转矩观测的永磁同步电机滑模控制方法。
背景技术
永磁同步电机具有效率高、扭矩大、转速性能好等优点,被广泛应用在制造、电动汽车、工业生产等领域中。永磁同步电机控制系统是参数时变的复杂系统,在发生负载扰动或内部参数摄动导致转矩偏差增大(包括控制器输出的给定转矩变化导致,或者是负载转矩发生变化导致)时,会影响永磁同步电机速度控制的稳态误差、动态性能和调速范围等性能参数。
发明内容
本发明的目的是将观测的负载转矩前馈补偿至电流调节器中时,为提高对负载转矩观测响应速度并减少转矩观测的波动性,提供一种基于负载转矩观测的永磁同步电机滑模控制方法,具体包括以下步骤:
步骤一、检测永磁同步电机的转子位置θ、转子角速度ω和三相电流ia、ib和ic
步骤二、依据三相电流ia、ib和ic对永磁同步电机进行Clark变换得到在α-β轴坐标系下的电流iα、电流iβ,依据电流iα、电流iβ和转子位置θ进行Park变换,得到在d-q轴坐标系下的电流id、电流iq
步骤三、负载转矩观测器依据负载转矩给定值TL *、负载转矩观测值
Figure BDA0002665882440000011
转子角速度ω和电流iq对负载转矩进行观测,得到新的负载转矩观测值
Figure BDA0002665882440000012
步骤四、滑模速度控制器依据输入的转子给定角速度ω*和转子角速度ω进行控制计算,得到负载转矩给定值TL *
步骤五、依据负载转矩给定值TL *和新的负载转矩观测值
Figure BDA0002665882440000013
计算得到q轴转矩电流给定值iq *
步骤六、d轴电流控制器依据d轴转矩电流给定值id *与d轴坐标系下的电流id之间的差值进行PI控制运算,得到d轴坐标系下的控制电压Ud;q轴电流控制器依据q轴转矩电流给定值iq *与q轴坐标系下的电流iq之间的差值进行PI控制运算,得到q轴坐标系下的控制电压Uq;依据d-q轴坐标系下的控制电压Ud、Uq进行Park逆变换,得到α-β轴坐标系下的控制电压Uα、Uβ
步骤七、将α-β轴坐标系下的控制电压Uα、Uβ作为SVPWM模块的输入,由SVPWM模块控制三相逆变器产生三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转。
所述步骤三中,得到新的负载转矩观测值
Figure BDA0002665882440000021
的具体方法是:
依据负载转矩观测器
Figure BDA0002665882440000022
计算得到新的负载转矩观测值
Figure BDA0002665882440000023
其中,U=kg·sgn(e1),
Figure BDA0002665882440000024
为转子角速度估算误差,
Figure BDA0002665882440000025
是转子角速度估计值,g是观测器的反馈增益,kg是观测器的滑模增益,ψf是永磁体磁链,J是转动惯量,p是电机极对数。
定义负载转矩观测器滑模面为
Figure BDA0002665882440000026
定义
Figure BDA0002665882440000027
为负载转矩观测误差;滑模增益kg的取值条件为:
Figure BDA0002665882440000028
反馈增益g的取值方法为:
步骤1、计算ΔTL *=TL *(k-1)-TL *(k-2),
Figure BDA0002665882440000029
步骤2、判断|ΔTL *|是否大于等于ε;当|ΔTL *|大于等于ε时,取反馈增益g等于gmin并退出;当|ΔTL *|小于ε时,进入步骤3;
步骤3、判断
Figure BDA00026658824400000210
是否大于等于ε;当
Figure BDA00026658824400000211
大于等于ε时,取反馈增益g等于gmin并退出;当
Figure BDA00026658824400000212
小于ε时,取反馈增益g等于gmax并退出。
其中,gmin<gmax<0,ε为转矩变化比较阈值且ε>0;TL *(k-1)为k-1时刻滑模速度控制器计算得到的负载转矩给定值TL *,TL *(k-2)为k-2时刻滑模速度控制器计算得到的负载转矩给定值TL *
Figure BDA00026658824400000213
为k-1时刻负载转矩观测器计算得到(观测得到)的负载转矩观测值
Figure BDA00026658824400000214
Figure BDA00026658824400000215
为k-2时刻负载转矩观测器计算得到的负载转矩观测值
Figure BDA00026658824400000216
步骤四中,滑模速度控制器的具体控制方法为:
定义状态变量为
Figure BDA00026658824400000217
定义滑模速度控制器滑模面为s=cx1+x2;其中c为滑模面参数,且c>0;滑模速度控制器的输出负载转矩给定值TL *
TL *=J·∫(cx2+k1|x1|sgn(s)+k2s)dt
其中,系数k1、k2均为常数且k1>0,k2>0。
步骤五中,计算得到q轴转矩电流给定值iq *的具体方法是:
Figure BDA0002665882440000031
其中,kq为转矩观测补偿系数,用于将转矩补偿转换为电流补偿,且按照式
Figure BDA0002665882440000032
进行计算。
所述基于负载转矩观测的永磁同步电机滑模控制方法由永磁同步电机速度控制系统实现,永磁同步电机速度控制系统包括滑模速度控制器、负载转矩观测器、q轴电流控制器、d轴电流控制器、Clarke变换模块、位置与速度检测模块、Park变换模块、Park逆变换模块、SVPWM模块和三相逆变器。
在步骤一至步骤七的永磁同步电机速度控制系统的周期性控制过程中,将当次k时刻(或者是第k步)滑模速度控制器计算得到的负载转矩给定值TL *记为TL *(k),负载转矩观测器计算得到的负载转矩观测值
Figure BDA0002665882440000033
记为
Figure BDA0002665882440000034
k-1时刻为k时刻的前一次周期性控制过程时刻,k-2时刻为k-1时刻的前一次周期性控制过程时刻。
本发明的有益效果是,将负载转矩观测值前馈补偿至电流调节器中,在不需要滑模速度控制器输出较大幅值的负载转矩给定值变化量的情况下,就能提供足够大的给定电流变化量来抵消负载发生扰动或者是系统参数发生变化的相关影响,减少滑模速度控制器的输出压力及不连续项的幅值,有效的削弱了系统的抖振。负载转矩观测采用反馈增益依据负载转矩给定值的变化量和负载转矩观测值的变化量进行自动调整的算法,避免了负载转矩观测器选择固定小反馈增益导致转矩观测波动大,选择固定大反馈增益导致收敛时间长的问题,能在系统参数发生变化或者是负载发生扰动,导致负载转矩给定值发生变化或/和负载转矩观测值发生变化时,快速消除负载转矩的观测误差,改善了观测效果和电机速度控制的快速性与精确性。反馈增益在负载转矩给定值变化时即进行自动调整,能够在给定转速改变使负载转矩给定值改变,或者是系统参数发生变化使负载转矩给定值先于负载转矩观测值发生改变时,提前调整反馈增益,当负载转矩观测值真正产生观测误差时,加快观测器响应速度,快速减小负载转矩观测值的观测误差,并进一步改善电机速度控制的快速性与精确性。
附图说明
图1为永磁同步电机速度控制系统框图;
图2为基于反馈增益自动调整的负载转矩观测流程图;
图3为反馈增益g取固定值gmin的负载转矩观测器仿真结果;
图4为反馈增益g取固定值gmed的负载转矩观测器仿真结果;
图5为反馈增益g取固定值gmax的负载转矩观测器仿真结果;
图6为反馈增益g自动调整时的负载转矩观测器仿真结果;
图7为固定反馈增益和自动调整增益两种控制策略下的转速响应仿真结果;
图8为转速变化时自动调整增益负载转矩观测器的转矩观测波形。
具体实施方式
以下结合附图及实施例,对本发明进行进一步详细说明。
图1为实现基于负载转矩观测的永磁同步电机滑模控制方法的永磁同步电机速度控制系统框图。图1中,Clarke变换模块输入永磁同步电机(即PMSM)的三相电流ia、ib和ic,输出两相静止α-β轴坐标系下的电流iα、iβ;位置与速度检测模块中的位置传感器检测永磁同步电机的转子位置θ后转换为转子角速度ω输出;Park变换模块输入电流iα、iβ和转子位置θ,输出旋转d-q轴坐标系下的电流id、iq;滑模速度控制器SMC输入转子给定角速度ω*和转子角速度ω,输出负载转矩给定值TL *和转矩电流给定分量i′q;负载转矩观测器输入负载转矩给定值TL *、转子角速度ω和电流iq,输出转矩电流补偿分量i″q;转矩电流给定分量i′q和转矩电流补偿分量i″q相加后,作为q轴转矩电流给定值iq *;q轴电流PI控制器输入q轴转矩电流给定值iq *和电流id,输出q轴坐标系下的控制电压Uq;d轴电流PI控制器输入q轴转矩电流给定值id *和电流id,输出d轴坐标系下的控制电压Ud,d轴转矩电流给定值id *等于0;Park逆变换模块输入d-q轴坐标系下的控制电压Ud、Uq,输出α-β轴坐标系下的控制电压Uα、Uβ;SVPWM模块(即空间矢量脉宽调制模块)输入控制电压Uα、Uβ,输出脉冲信号至三相逆变器,三相逆变器将直流电压Udc转换为三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转。
忽略铁心涡流与磁滞损耗等影响,采用id=0的PMSM转子磁场定向控制,建立PMSM在d-q轴旋转坐标系下的数学模型,电压方程为:
Figure BDA0002665882440000041
对于采用id=0控制方式的凸极式PMSM矢量控制系统,电磁转矩方程为:
Figure BDA0002665882440000042
PMSM运动方程为:
Figure BDA0002665882440000043
式(1)(2)(3)中,ud、uq分别是d-q轴的电压;id、iq分别是d-q轴的电流;Ld、Lq分别是d-q轴的电感;Te是电磁转矩;TL是负载转矩;R是定子的电阻;p是电机极对数;ωe是转子电角速度,即角频率;ω是转子角速度,即转子机械角速度;ψf是永磁体磁链;J是转动惯量;B是摩擦系数。
定义PMSM系统的状态变量为:
Figure BDA0002665882440000051
其中,ω*是电机的给定转子角速度。由式(2)、(3)、(4)得:
Figure BDA0002665882440000052
对式(5)进行化简,令D=1.5pψf/J,
Figure BDA0002665882440000053
可得系统状态空间方程为:
Figure BDA0002665882440000054
选择滑模面函数为:
s=cx1+x2 (7)
式(7)中,c为滑模面参数,且c>0。对式(7)求导可得:
Figure BDA0002665882440000055
传统指数趋近律的表达式为:
Figure BDA0002665882440000056
式(9)中,sgn()为符号函数,-k1sgn(s)为等速趋近项,-k2s为指数趋近项,k1、k2两个系数分别决定滑模面的抖振和趋近过程的运动品质,且k1、k2均大于0。为提高系统响应速度,在传统指数趋近率的基础上改进,改进的趋近律为:
Figure BDA0002665882440000057
在趋近律等速趋近项中加入转速误差绝对值|x1|,使系统控制量在趋近滑模面过程中的速度与系统转速误差的绝对值|x1|相关联。结合式(8)和(10),得滑模速度控制器输出的负载转矩给定值TL *和转矩电流给定分量i′q为:
Figure BDA0002665882440000058
滑模速度控制器输出中含有积分项,对输出的控制量进行滤波,可削弱系统抖振、减少系统的稳态误差。定义Lyapunov函数为:
Figure BDA0002665882440000059
由式(10)和式(12)可得:
Figure BDA0002665882440000061
式(13)中,k1>0,k2>0,s·sgn(s)≥0,故
Figure BDA0002665882440000062
表明系统跟踪误差能在有限时间内收敛到零,系统能稳定运行。此时,k1和k2分别为变速趋近项系数和指数趋近项系数,系数k2越大,系统状态能以较大速度趋近于滑动模态;系数k1决定到达切换面的速度,k1越小则穿越切换面的距离和抖动越小。根据经验,进行小功率永磁同步电机控制时,c取值通常小于100,系数k1取值小于1000,系数k2取值小于5000。
根据PMSM电磁转矩与运动方程,对于恒定的阶跃性负载,在变化周期内可认为是一恒定值,即
Figure BDA0002665882440000063
将电机角速度与负载转矩作为状态变量,构成PMSM状态方程为:
Figure BDA0002665882440000064
电机摩擦与负载转矩相比较,所占比重小,在式(14)基础上,令B=0,忽略摩擦的影响,且以负载转矩与电机角速度为观测对象,建立负载转矩观测器为:
Figure BDA0002665882440000065
式(15)中,
Figure BDA0002665882440000066
是负载转矩的观测值,
Figure BDA0002665882440000067
是转子角速度估计值,g是负载转矩观测器的反馈增益,
Figure BDA0002665882440000068
kg是负载转矩观测器的滑模增益。
依据式(14)和式(15)得到负载转矩观测器误差方程为:
Figure BDA0002665882440000069
式(16)中,
Figure BDA00026658824400000610
为转子角速度估算误差,
Figure BDA00026658824400000611
为负载转矩观测误差,并定义观测器滑模面为
Figure BDA00026658824400000612
根据滑模可达性条件,观测器系统稳定条件为kg≤-|e2/J|,且g<0,负载转矩观测误差为:
Figure BDA00026658824400000613
式(17)中,c1为常数,观测误差e2随时间t变化以指数方式逐渐趋近到零。反馈增益g的取值大小对负载转矩观测结果影响较大。反馈增益g越大,观测转矩波动性越小,但观测转矩辨识速度越慢;反馈增益g越小,观测转矩速度越快,但观测转矩波动性越大。出于对此问题的考虑,在传统负载转矩观测器中,综合考虑负载转矩的观测速度与波动性,将反馈增益g取一个折中值,但这会舍弃大反馈增益时的波动性小和小反馈增益时的观测速度快的优势。
滑模控制主要通过增大控制器中不连续项的幅值来抑制参数变化和外部负载扰动对系统的影响,但幅值增大会引起滑模固有抖振。为解决滑模控制系统抖振与抗扰性之间的矛盾,利用观测器实时观测负载扰动变化,将负载转矩观测值前馈补偿至电流调节器中,以降低滑模控制中不连续项幅值,削弱参数变化导致给定转矩变化,或者是负载扰动所引起系统抖振。为了充分利用反馈增益g在高、低值时的优势,根据两相邻时刻负载转矩观测值变化量以及负载转矩给定值变化量的大小,在负载转矩给定值变化小且负载转矩观测值变化小时,给予反馈增益g较大值,使观测结果波动性小,稳定性更强;在负载转矩给定值变化大或者负载转矩观测值变化大时,给予反馈增益g较小值,使观测速度加快,最终通过对反馈增益g参数的自动调整,得到观测速度快和波动小、稳定性更强的综合结果。
图2为基于反馈增益自动调整的负载转矩观测流程图。ε为转矩变化比较阈值,将负载转矩给定值和负载转矩观测值前后时刻差值与ε进行比较来判断当前给定值和观测值的变化波动情况。图2中,ΔTL *为最近2次负载转矩给定值之间的差值,
Figure BDA0002665882440000071
为最近2次负载转矩观测值之间的差值;当|ΔTL *|大于等于ε时,表明负载转矩给定值因参数变化,或者是因转子角速度给定值变化等因素而处于较大变化状态,反馈增益g选择等于gmin进行转矩辨识与观测;当|ΔTL *|小于ε,且
Figure BDA0002665882440000072
大于等于ε时,表明负载转矩观测值波动大,反馈增益g选择等于gmin进行转矩辨识与观测;当|ΔTL *|小于ε且
Figure BDA0002665882440000073
小于ε时,表明负载转矩给定值和负载转矩观测值的变化波动均小,反馈增益g选择等于gmax进行转矩辨识与观测。图2中,gmin<gmax<0,比较阈值ε>0。
基于负载转矩观测的永磁同步电机滑模控制方法是离散化后,通过计算机,或者是具有CPU等控制核心中运行的程序来实现。在周期性的控制过程中,将当次k时刻(或者是第k步)计算得到的负载转矩给定值TL *记为TL *(k),计算得到的负载转矩观测值
Figure BDA0002665882440000074
记为
Figure BDA0002665882440000075
k-1时刻为k时刻的前一次周期性控制过程时刻,计算得到的负载转矩给定值TL *记为TL *(k-1),计算得到的负载转矩观测值
Figure BDA0002665882440000076
记为
Figure BDA0002665882440000077
k-1时刻为k-1时刻的前一次周期性控制过程时刻,计算得到的负载转矩给定值TL *记为TL *(k-2),计算得到的负载转矩观测值
Figure BDA0002665882440000078
记为
Figure BDA0002665882440000079
图2中,ΔTL *=TL *(k-1)-TL *(k-2),
Figure BDA00026658824400000710
观测得到负载转矩观测值
Figure BDA00026658824400000711
后,将负载转矩的观测值
Figure BDA00026658824400000712
转换成转矩电流补偿分量i″q前馈补偿至q轴电流PI控制器的输入,对滑模速度控制器输出的转矩电流给定分量i′q进行补偿。q轴电流PI控制器的q轴转矩电流给定值iq *为:
Figure BDA00026658824400000713
式(18)中,kq=2/(3pψf)为转矩观测补偿系数。比较式(11)与式(18)可得,当负载发生扰动或者是系统参数发生变化时,式(11)中没有加入负载转矩补偿,需要选取较大的k1、k2值来提供足够大的给定电流变化量来抵消负载发生扰动或者是系统参数发生变化的相关影响,以保证电机转速能够快速恒定;而式(18)将负载转矩观测值前馈补偿至电流调节器中,在不需要较大k1、k2值的情况下就能在负载发生扰动或者是系统参数发生变化时,提供足够大的给定电流变化量来抵消负载发生扰动或者是系统参数发生变化的相关影响,减少滑模速度控制器的输出压力及不连续项的幅值,有效地削弱了系统的抖振。
滑模增益kg按照
Figure BDA0002665882440000081
进行选择。式(19)中,α≥1;一般地,α值在1~5的范围内选择,例如,选择α等于1.5。kg的绝对值选择过小,当|e2|较大时观测器不能进入滑模状态;kg的绝对值选择足够大,可保证观测器进入滑模状态,但负载转矩的稳态观测波动变大;kg的值随着负载转矩观测误差的变化而变化,可以同时兼顾增加观测器稳定性和减小负载转矩的稳态观测波动。
ε在大于0,且一般小于额定转矩5%的范围内取值,例如,额定转矩为22N·m,取ε=0.2N·m。同样地,反馈增益g的取值满足gmin<gmax<0,一般情况下,gmin≥-5000;gmin取值应使负载转矩观测器输出的观测值的波动处在合理范围之内,例如,观测值的波动在额定转矩的10%之内,或者是在额定转矩的5%之内,等等;gmax取值应使负载转矩不变时,负载转矩观测器输出的观测值波动明显小于ε,例如,选择反馈增益gmax=-0.5,gmin=-10。选取gmin、gmax、ε值的具体方法是:
步骤①,负载转矩观测器和滑模速度控制器均处于稳态且保持给定转子角速度不变和负载转矩不变;
步骤②,令反馈增益g从一个较大值开始逐渐减小,例如,令反馈增益g从-0.01开始逐渐减小,当负载转矩观测的稳态误差达到负载转矩观测稳态误差限值时,确定此时的反馈增益g值为gmax
步骤③,保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmax,连续进行u次
Figure BDA0002665882440000082
值的测量,并将u次
Figure BDA0002665882440000083
测量中最大m个
Figure BDA0002665882440000084
测量值的平均值作为转矩变化比较阈值ε;u≥10;2≤m≤0.5u;
步骤④,微调改变反馈增益g,负载转矩观测器和滑模速度控制器均处于稳态时,保持给定转子角速度不变且令负载转矩突变,在保证转矩观测跟踪超调量处于转矩观测跟踪超调限值之内的前提下,测量负载转矩观测器的跟踪调节时间;
步骤⑤,重复步骤④,选择跟踪调节时间最短的反馈增益g值为gmin值。通常情况下,转矩观测跟踪超调量接近转矩观测跟踪超调限值时,负载转矩观测器的跟踪调节时间较短。
转矩观测跟踪超调限值通常为电机额定转矩的1%~10%,具体来说,转矩观测跟踪超调限值是额定转矩的2%,或者是额定转矩的5%,或者是额定转矩的10%,等等。负载转矩从一个定值突变至另外一个定值,突变开始时刻至负载转矩观测器输出观测值稳定地进入负载转矩观测稳态误差限值范围时刻为转矩观测过渡过程,跟踪调节时间指的是该过渡过程时间;负载转矩观测稳态误差指的是负载转矩不变且负载转矩观测器处于稳态时观测转矩瞬时值与负载转矩之间的误差,该误差包括滑模观测器自身抖振引起的观测误差和负载波动之外干扰原因所造成的观测误差,或者是状态观测器因为转子角速度抖振引起的观测误差和负载波动之外干扰原因所造成的观测误差;负载转矩观测稳态误差限值是负载转矩观测器允许的负载转矩观测稳态误差最大绝对值;负载转矩观测稳态误差限值通常为电机额定转矩的1%~5%,具体来说,负载转矩观测稳态误差限值是额定转矩的1%,或者是额定转矩的2%,或者是额定转矩的5%,等等。转矩观测跟踪超调量指的是负载转矩从一个定值突变至另外一个定值,负载转矩观测器输出的观测值超过突变后负载转矩的最大偏差值。当负载转矩观测的稳态误差在负载转矩观测稳态误差限值的临近范围之内,例如,95%~105%范围内,或者是98%~102%范围内时,认为负载转矩观测的稳态误差增大到负载转矩观测稳态误差限值。滑模速度控制器处于稳态指的是滑模速度控制器稳定地处于滑动模态;滑模观测器处于稳态指的是滑模观测器稳定地处于滑动模态。
依据图1的系统框图在Simulink中搭建仿真模型,选取的永磁同步电机参数见表1。忽略摩擦,取摩擦系数B为0。
表1永磁同步电机参数
Figure BDA0002665882440000091
选取反馈增益高值为gmax=-0.5,反馈增益低值为gmin=-20,反馈增益折中值gmed=-3,负载从0N·m突变至22N·m时对永磁同步电机的负载转矩进行观测,反馈增益g分别取固定值gmin、gmed、gmax的负载转矩观测器仿真结果分别如图3、图4、图5所示;反馈增益g按照前述方法自动调整时的负载转矩观测器仿真结果如图6所示,设置转矩变化比较阈值ε为0.2。图3至图6中,虚线为实际转矩值,实线为观测转矩值,不同反馈增益取法下观测效果对比见表2。设观测转矩值的允许误差范围是0.25N·m,收敛时间为负载突变至观测转矩值稳定进入误差范围的时间,滑模态观测波动为负载转矩观测器进入滑模态后的观测波动误差。
表2不同反馈增益取法下观测效果对比
Figure BDA0002665882440000101
从表2可以看出,固定反馈增益值时,反馈增益g越小,负载转矩观测的滑模态观测波动振荡幅值越大,波动性越强;反馈增益g越大,负载转矩观测的滑模态观测波动振荡幅值越小,观测精度越高。自动调整增益算法解决了负载转矩观测器中小反馈增益导致转矩观测波动大的问题,大反馈增益收敛时间长的问题,收敛时间和波动幅度指标都优于折中增益算法。可变增益算法能快速跟踪负载转矩变化值,将传统负载转矩观测器的固定反馈增益优点保留,对缺点进行改进,达到了较好的观测效果。
在额定负载转矩下,设系统初始转速为1500r/min,0.2s时转速突减至1200r/min,之后在0.4s时突增转速至1400r/min,图7为固定反馈增益和自动调整增益两种控制策略下转速响应仿真结果,其中的虚线为固定反馈增益前馈补偿控制(即传统SMC前馈补偿控制)时的转速变化曲线,实线为自动调整增益前馈补偿控制(即改进SMC前馈补偿控制)时的转速变化曲线;可以看出在转速给定突变时,固定反馈增益前馈补偿控制超调小,但调节时间长,需经过16ms左右恢复至给定转速;而自动调整增益前馈补偿控制超调较小,调节时间短,响应速度快,只需经过6ms左右恢复至给定转速。图8为转速变化时自动调整增益负载转矩观测器对转矩的观测波形,虚线为实际转矩值,实线为观测转矩值,可以看出当转速稳定时能很好观测负载转矩,在转速变化时,观测转矩出现波动,但观测值误差不超过0.004N·m,响应时间在20ms以内,观测效果较好。
图7中的纵轴是电机转速n,单位是r/min;转子角速度ω=2πn/60rad/s。
在额定负载转矩下改变给定转速时,尽管实际的负载转矩没有改变,但从式(15)-(17)所构造的负载转矩观测器可知,当转子角速度ω发生改变时,即使负载转矩没有改变,观测到的转矩观测值也会发生变化,造成观测误差。在额定负载转矩下改变给定转速时,永磁同步电机滑模控制系统的控制调节过程是,首先由滑模速度控制器依据给定速度变化,使输出的负载转矩给定值TL *改变,以及使转矩电流给定值iq *发生改变,进而使永磁同步电机的电磁转矩Te发生改变,带动电机使转子角速度ω发生改变;反馈增益g如果只依据负载转矩观测值的变化量
Figure BDA0002665882440000111
进行自动调整,则此时只有当转子角速度ω发生改变,使负载转矩观测值
Figure BDA0002665882440000112
发生了改变后,才对反馈增益g进行调整;反馈增益g同时依据负载转矩给定值的变化量ΔTL *和负载转矩观测值的变化量
Figure BDA0002665882440000113
进行自动调整,当给定转速改变使负载转矩给定值TL *改变,负载转矩观测值
Figure BDA0002665882440000114
还没有发生变化时,就提前调整反馈增益g,当负载转矩观测值
Figure BDA0002665882440000115
真正产生观测误差时,能够加快观测器响应速度,尽快消除负载转矩观测值
Figure BDA0002665882440000116
的观测误差,进而改善电机速度控制的快速性与精确性。同样地,当系统参数发生变化使负载转矩给定值TL *先于负载转矩观测值
Figure BDA0002665882440000117
发生改变时,反馈增益g同时依据负载转矩给定值的变化量ΔTL *和负载转矩观测值的变化量
Figure BDA0002665882440000118
进行自动调整,能够提前调整反馈增益g,加快观测器响应速度,尽快消除负载转矩观测值
Figure BDA0002665882440000119
的观测误差,并进一步改善电机速度控制的快速性与精确性。
除说明书所述的技术特征外,本发明所涉及的其他技术均为本领域技术人员所掌握的常规技术。例如,q轴电流控制器、d轴电流控制器采用PI控制器进行控制及控制器参数的选择,滑模速度控制器控制参数的选择,位置与速度检测模块使用旋转变压器或光电编码器等进行永磁同步电机磁极位置、转子旋转角度与旋转速度检测,以及对于Clarke变换模块、Park变换模块、Park逆变换模块和SVPWM模块、三相逆变器的变换方法及应用方法,等等,均为本领域技术人员所掌握的常规技术。

Claims (3)

1.一种基于负载转矩观测的永磁同步电机滑模控制方法,其特征在于,包括以下步骤:
步骤一、检测永磁同步电机的转子位置θ、转子角速度ω和三相电流ia、ib和ic
步骤二、依据三相电流ia、ib和ic对永磁同步电机进行Clark变换得到在α-β轴坐标系下的电流iα、电流iβ,依据电流iα、电流iβ和转子位置θ进行Park变换,得到在d-q轴坐标系下的电流id、电流iq
步骤三、负载转矩观测器依据负载转矩给定值TL *、负载转矩观测值
Figure FDA0003622769670000011
转子角速度ω和电流iq对负载转矩进行观测,得到新的负载转矩观测值
Figure FDA0003622769670000012
步骤四、滑模速度控制器依据输入的转子给定角速度ω*和转子角速度ω进行控制计算,得到负载转矩给定值TL *
步骤五、依据负载转矩给定值TL *和新的负载转矩观测值
Figure FDA0003622769670000013
计算得到q轴转矩电流给定值iq *
步骤六、d轴电流控制器依据d轴转矩电流给定值id *与d轴坐标系下的电流id之间的差值进行PI控制运算,得到d轴坐标系下的控制电压Ud;q轴电流控制器依据q轴转矩电流给定值iq *与q轴坐标系下的电流iq之间的差值进行PI控制运算,得到q轴坐标系下的控制电压Uq;依据d-q轴坐标系下的控制电压Ud、Uq进行Park逆变换,得到α-β轴坐标系下的控制电压Uα、Uβ
步骤七、将α-β轴坐标系下的控制电压Uα、Uβ作为SVPWM模块的输入,由SVPWM模块控制三相逆变器产生三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转;
所述步骤三中,得到新的负载转矩观测值
Figure FDA0003622769670000014
的具体方法是:
依据负载转矩观测器
Figure FDA0003622769670000015
计算得到新的负载转矩观测值
Figure FDA0003622769670000016
其中,U=kg·sgn(e1),
Figure FDA0003622769670000017
为转子角速度估算误差,
Figure FDA0003622769670000018
是转子角速度估计值,g是观测器的反馈增益,kg是观测器的滑模增益,ψf是永磁体磁链,J是转动惯量,p是电机极对数;
定义观测器滑模面为
Figure FDA0003622769670000019
定义
Figure FDA00036227696700000110
为负载转矩观测误差;滑模增益kg的取值条件为:
Figure FDA00036227696700000111
反馈增益g的取值方法为:
步骤1、计算ΔTL *=TL *(k-1)-TL *(k-2),
Figure FDA0003622769670000021
步骤2、判断|ΔTL *|是否大于等于ε;当|ΔTL *|大于等于ε时,取反馈增益g等于gmin并退出;当|ΔTL *|小于ε时,进入步骤3;
步骤3、判断
Figure FDA0003622769670000022
是否大于等于ε;当
Figure FDA0003622769670000023
大于等于ε时,取反馈增益g等于gmin并退出;当
Figure FDA0003622769670000024
小于ε时,取反馈增益g等于gmax并退出;
其中,gmin<gmax<0,ε为转矩变化比较阈值且ε>0。
2.如权利要求1所述的基于负载转矩观测的永磁同步电机滑模控制方法,其特征在于,步骤四中,滑模速度控制器的具体控制方法为:
定义状态变量为
Figure FDA0003622769670000025
定义滑模速度控制器滑模面为s=cx1+x2;其中c为滑模面参数,且c>0;滑模速度控制器的输出负载转矩给定值TL *
TL *=J·∫(cx2+k1|x1|sgn(s)+k2s)dt;
其中,系数k1、k2均为常数且k1>0,k2>0。
3.如权利要求2所述的基于负载转矩观测的永磁同步电机滑模控制方法,其特征在于,步骤五中,计算得到q轴转矩电流给定值iq *的具体方法是:
Figure FDA0003622769670000026
其中,kq为转矩观测补偿系数,且按照式
Figure FDA0003622769670000027
进行计算。
CN202010918601.6A 2020-09-04 2020-09-04 基于负载转矩观测的永磁同步电机滑模控制方法 Active CN112039390B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010918601.6A CN112039390B (zh) 2020-09-04 2020-09-04 基于负载转矩观测的永磁同步电机滑模控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010918601.6A CN112039390B (zh) 2020-09-04 2020-09-04 基于负载转矩观测的永磁同步电机滑模控制方法

Publications (2)

Publication Number Publication Date
CN112039390A CN112039390A (zh) 2020-12-04
CN112039390B true CN112039390B (zh) 2022-07-26

Family

ID=73592381

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010918601.6A Active CN112039390B (zh) 2020-09-04 2020-09-04 基于负载转矩观测的永磁同步电机滑模控制方法

Country Status (1)

Country Link
CN (1) CN112039390B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113078861B (zh) * 2021-04-07 2023-04-25 长春工业大学 一种永磁同步电机滑模控制方法、系统、介质及应用
CN113131814B (zh) * 2021-04-25 2022-11-18 中南大学 用于低精度编码器场合的异步电机负载转矩观测方法
CN114172426B (zh) * 2021-10-27 2024-06-11 北京自动化控制设备研究所 一种永磁同步电机系数自整定速度补偿控制方法
CN114421836A (zh) * 2022-01-25 2022-04-29 中国船舶重工集团公司第七二四研究所 一种基于转矩观测的永磁同步电机控制方法
CN114726273A (zh) * 2022-03-30 2022-07-08 电子科技大学 一种基于改进滑模观测器的pcb盘式永磁同步电机无感控制系统及方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101714844B (zh) * 2009-11-10 2012-05-02 哈尔滨工业大学 一种内置式永磁同步电机转子磁极初始位置检测方法
CN107070335B (zh) * 2016-12-30 2019-05-24 湖南大学 双pwm永磁电力驱动系统转矩前馈控制方法及其控制装置
CN107359837B (zh) * 2017-09-06 2024-01-09 深圳市立三机电有限公司 基于滑模观测器和自抗扰的永磁同步电机转矩控制方法
CN110557070A (zh) * 2019-09-30 2019-12-10 山东深川变频科技股份有限公司 基于二阶滑模观测器的永磁同步电机参数辨识方法

Also Published As

Publication number Publication date
CN112039390A (zh) 2020-12-04

Similar Documents

Publication Publication Date Title
CN112039390B (zh) 基于负载转矩观测的永磁同步电机滑模控制方法
CN107046387B (zh) 一种永磁同步电机的变pid参数电流环启动方法
CN112422004B (zh) 一种永磁同步电机弱磁控制模式下的扰动抑制方法
CN110752806B (zh) 改进趋近律的内置式永磁同步电机的滑模转速控制方法
CN104601071B (zh) 基于扰动观测器的永磁同步电机电流环滑模控制系统
CN110165953B (zh) 一种基于趋近律的pmsm调速控制方法
CN110138298B (zh) 一种永磁同步电机滑模控制方法
CN112003525B (zh) 一种矢量控制永磁同步风力发电系统转速环自抗扰调节器
CN112290843B (zh) 一种变指数幂次趋近律及其pmsm控制应用
Rong et al. A new PMSM speed modulation system with sliding mode based on active-disturbance-rejection control
CN113206623A (zh) 基于快速积分终端滑模及干扰估计的永磁同步电机有限时间调速控制方法
CN109167538B (zh) 基于双层结构抑制永磁无刷直流电机转矩脉动的控制方法
CN115276491A (zh) 一种永磁同步电机速度优化控制方法
CN105071723A (zh) 三步法设计的有刷直流电机复合控制方法
CN116526919A (zh) 一种永磁同步电机伺服系统及其电流预测控制方法、装置
CN104977850A (zh) 一种基于分数阶预测器的无时滞鲁棒伺服电机控制方法
CN108448983B (zh) 一种抑制伺服系统极低速时非线性干扰的方法
CN112039388B (zh) 工业机器人用永磁同步电机控制方法
CN114900086A (zh) 一种永磁同步电机负载转矩滑模观测方法
CN112039389B (zh) 矿用牵引永磁同步电机驱动控制方法
CN115459667B (zh) 基于改进型趋近律的永磁同步电机无传感器滑模控制方法
CN110943666A (zh) 一种永磁同步电机复合电流的约束控制系统及其构建方法
Li et al. Composite fractional order sliding mode control of permanent magnet synchronous motor based on disturbance observer
Bai et al. Integrated Sliding Mode Velocity Control of Linear Permanent Magnet Synchronous Motor with Thrust Ripple Compensation
CN115378325A (zh) 基于动态权重因子的smpmsm驱动系统直接速度复合控制方法

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
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20221221

Address after: No. 268 Binhai Avenue, Binhai New Area, Fenghua Economic Development Zone, Ningbo, Zhejiang 315500

Patentee after: Ningbo Fenghua Wanyun Motor Co.,Ltd.

Address before: Room 2202, 22 / F, Wantong building, No. 3002, Sungang East Road, Sungang street, Luohu District, Shenzhen City, Guangdong Province

Patentee before: Shenzhen dragon totem technology achievement transformation Co.,Ltd.

Effective date of registration: 20221221

Address after: Room 2202, 22 / F, Wantong building, No. 3002, Sungang East Road, Sungang street, Luohu District, Shenzhen City, Guangdong Province

Patentee after: Shenzhen dragon totem technology achievement transformation Co.,Ltd.

Address before: Department of science and technology, Hunan University of technology, No. 88, Taishan West Road, Zhuzhou City, Hunan Province

Patentee before: HUNAN University OF TECHNOLOGY