CN112039389B - 矿用牵引永磁同步电机驱动控制方法 - Google Patents

矿用牵引永磁同步电机驱动控制方法 Download PDF

Info

Publication number
CN112039389B
CN112039389B CN202010918598.8A CN202010918598A CN112039389B CN 112039389 B CN112039389 B CN 112039389B CN 202010918598 A CN202010918598 A CN 202010918598A CN 112039389 B CN112039389 B CN 112039389B
Authority
CN
China
Prior art keywords
load torque
value
torque
feedback gain
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
CN202010918598.8A
Other languages
English (en)
Other versions
CN112039389A (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.)
Hunan University of Technology
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 CN202010918598.8A priority Critical patent/CN112039389B/zh
Priority to CN202210291754.1A priority patent/CN114844405A/zh
Priority to CN202210291698.1A priority patent/CN114865968A/zh
Publication of CN112039389A publication Critical patent/CN112039389A/zh
Application granted granted Critical
Publication of CN112039389B publication Critical patent/CN112039389B/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/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/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/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
    • 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
    • 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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02TCLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T10/00Road transport of goods or passengers
    • Y02T10/60Other road transportation technologies with climate change mitigation effect
    • Y02T10/72Electric energy management in electromobility

Abstract

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

Description

矿用牵引永磁同步电机驱动控制方法
技术领域
本发明涉及永磁同步电机技术领域,更具体地,尤其是涉及一种矿用牵引永磁同步电机驱动控制方法。
背景技术
永磁同步电机具有效率高、扭矩大、转速性能好等优点,被广泛应用在制造、电动汽车、工业生产等领域中。矿用牵引电机应用场合工作环境复杂多变,负载转矩变化大,调速范围广,且要求能够大力矩启动,基于PI控制器的永磁同步电机矢量控制方法的鲁棒性不如滑模控制方法,但滑模控制方法在发生负载扰动或内部参数摄动时,会使电机速度产生明显的抖振。
发明内容
本发明的目的是针对负载转矩变化大和调速范围广的情况,提供一种将观测的负载转矩前馈补偿至电流调节器中,提高负载转矩观测响应速度并减少转矩观测波动性的矿用牵引永磁同步电机驱动控制方法,包括:
永磁同步电机速度由滑模速度控制器控制,负载转矩观测器对负载转矩进行观测,负载转矩观测器的输出用于对滑模速度控制器的输出进行负载转矩补偿;其特征在于,负载转矩观测器根据负载转矩给定值和负载转矩观测值的变化对反馈增益进行调整,依据转子角速度ω和电流iq对负载转矩进行观测,得到新的负载转矩观测值;q轴转矩电流给定值为转矩电流给定分量与转矩电流补偿分量之和。
滑模速度控制器的状态变量为
Figure BDA0002665884180000011
其中,ω是转子角速度,ω*是给定转子角速度;滑模速度控制器的滑模面为s=cx1+x2,c为滑模面参数,且c>0。滑模速度控制器输出的负载转矩给定值
Figure BDA0002665884180000012
和转矩电流给定分量i′q
Figure BDA0002665884180000013
其中,J是转动惯量,p是电机极对数,ψf是永磁体磁链;系数k1、k2、k3、k4为速度滑模控制的指数趋近率系数,且k1>0,k2>0,1<k3<2,k4>0。
负载转矩观测器为
Figure BDA0002665884180000021
其中,
Figure BDA0002665884180000022
是转子角速度估计值,g是负载转矩观测器的反馈增益且g<0;
Figure BDA00026658841800000213
Figure BDA0002665884180000023
kg是负载转矩观测器的滑模增益且kg≤-|e2/J|,
Figure BDA0002665884180000024
为负载转矩观测误差,TL为负载转矩。
负载转矩观测器根据负载转矩给定值和负载转矩观测值的变化对反馈增益进行调整的方法是:
步骤1、计算
Figure BDA0002665884180000025
计算
Figure BDA0002665884180000026
步骤2、判断ΔT是否大于ε2;当ΔT大于ε2时,取反馈增益g等于gmin并进入步骤5;当ΔT小于等于ε2时,进入步骤3;
步骤3、判断ΔT是否小于ε1;当ΔT小于ε1时,取反馈增益g等于gmax并进入步骤5;当ΔT大于等于ε1时,进入步骤4;
步骤4、反馈增益g按照
Figure BDA0002665884180000027
进行计算;
步骤5、负载转矩观测器对负载转矩进行TL观测,得到负载转矩观测值
Figure BDA0002665884180000028
滑模速度控制器进行控制运算得到负载转矩给定值
Figure BDA0002665884180000029
其中,ε1为转矩变化低限比较阈值,ε2为转矩变化高限比较阈值,且0<ε1<ε2;gmax为反馈增益高值,gmin为反馈增益低值,且gmin<gmax<0。
转矩电流补偿分量i″q
Figure BDA00026658841800000210
q轴转矩电流给定值
Figure BDA00026658841800000211
Figure BDA00026658841800000212
选取gmin、gmax、ε1、ε2值的方法是:
步骤⑴,负载转矩观测器和滑模速度控制器均处于稳态且保持给定转子角速度不变和负载转矩不变;滑模速度控制器处于稳态指的是滑模速度控制器稳定地处于滑动模态;负载转矩观测器处于稳态指的是观测器稳定地处于滑动模态;
步骤⑵,令反馈增益g从一个较大值开始逐渐减小,例如,令反馈增益g从-0.01开始逐渐减小,当负载转矩观测的稳态误差达到负载转矩观测稳态误差限值时,确定此时的反馈增益g值为gmax
步骤⑶,保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmax,连续进行n次ΔT值的测量,并将此时n次测量中最大m个ΔT值之和的平均值作为转矩变化低限比较阈值ε1
步骤⑷,微调改变反馈增益g,负载转矩观测器和滑模速度控制器均处于稳态时,保持给定转子角速度不变和令负载转矩突变,在保证负载转矩观测器输出观测值的转矩观测跟踪超调量处于转矩观测跟踪超调限值之内的前提下,测量负载转矩观测器的跟踪调节时间;
步骤⑸,重复步骤⑷,选择跟踪调节时间最短的反馈增益g值为gmin值;
步骤⑹,再次保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmin,连续进行n次ΔT值的测量,并将此时n次测量中最大m个ΔT值之和的平均值作为转矩变化高限比较阈值ε2
所述n≥20,5≤m≤0.5n。
滑模增益kg按照
Figure BDA0002665884180000031
进行选择;其中,α≥1;优选地,1≤α≤5。
永磁同步电机速度控制的过程具体包括以下步骤:
步骤一、检测永磁同步电机的转子位置θ、转子角速度ω和三相电流ia、ib和ic
步骤二、依据三相电流ia、ib和ic对永磁同步电机进行Clark变换得到在α-β轴坐标系下的电流iα、电流iβ,依据电流iα、电流iβ和转子位置θ进行Park变换,得到在d-q轴坐标系下的电流id、电流iq
步骤三、负载转矩观测器的反馈增益g依据负载转矩给定值TL *和负载转矩观测值
Figure BDA0002665884180000032
的变化进行调整;
步骤四、负载转矩观测器依据转子角速度ω和电流iq对负载转矩进行观测,得到新的负载转矩观测值
Figure BDA0002665884180000033
和转矩电流补偿分量i″q
步骤五、滑模速度控制器依据输入的转子给定角速度ω*和转子角速度ω进行控制计算,得到负载转矩给定值
Figure BDA0002665884180000034
和转矩电流给定分量i′q
步骤六、依据转矩电流给定分量i′q和转矩电流补偿分量i″q计算得到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β;d轴转矩电流给定值id *等于0;
步骤七、将α-β轴坐标系下的控制电压Uα、Uβ作为SVPWM模块的输入,由SVPWM模块控制三相逆变器产生三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转。
上述周期性控制过程中,步骤三与步骤四、五的先后顺序可以互换,即可以先执行步骤四、五,后执行步骤三。
所述永磁同步电机速度控制由包括滑模速度控制器、负载转矩观测器、q轴电流控制器、d轴电流控制器、Clarke变换模块、位置与速度检测模块、Park变换模块、Park逆变换模块、SVPWM模块和三相逆变器的永磁同步电机速度控制系统实现。
本发明的有益效果是,将负载转矩观测值前馈补偿至电流调节器的给定值中,在不需要滑模速度控制器输出的给定电流部分产生较大调整的情况下,就能抵消负载发生扰动或者是系统参数发生变化所造成的相关影响,有效地削弱了系统的抖振。负载转矩观测采用反馈增益依据负载转矩给定值的变化量和负载转矩观测值的变化量进行自动调整的算法,避免了负载转矩观测器选择固定小反馈增益导致转矩观测波动大,选择固定大反馈增益导致收敛时间长的问题,能在系统的控制参数、模型参数等发生变化或者是负载发生扰动,导致负载转矩给定值发生变化或/和负载转矩观测值发生变化时,快速降低负载转矩的观测误差,改善了观测效果和电机速度控制的快速性与精确性。反馈增益在负载转矩给定值变化时即进行自动调整,能够在负载转矩观测值还没有发生较大变化,但因转子角速度给定值改变或/和转子角速度实际值改变使负载转矩给定值改变,或者是系统模型参数发生变化使负载转矩给定值发生改变,将引起负载转矩观测值有较大波动时,提前调整反馈增益,当负载转矩观测值真正产生观测误差时,加快观测器的响应速度,快速降低负载转矩观测值的观测误差,并进一步改善电机速度控制的快速性与精确性。
附图说明
图1为永磁同步电机速度控制系统实施例1框图;
图2为反馈增益自动调整方法实施例1流程图;
图3为永磁同步电机速度控制系统实施例2框图;
图4为反馈增益自动调整方法实施例2流程图。
具体实施方式
以下结合附图及实施例,对本发明进行进一步详细说明。
矿用牵引永磁同步电机驱动控制方法由永磁同步电机速度控制系统实现,图1为永磁同步电机速度控制系统实施例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轴转矩电流给定值i* q;q轴电流PI控制器输入q轴转矩电流给定值i* q和电流id,输出q轴坐标系下的控制电压Uq;d轴电流PI控制器输入q轴转矩电流给定值i* d和电流id,输出d轴坐标系下的控制电压Ud,d轴转矩电流给定值i* d等于0;Park逆变换模块输入d-q轴坐标系下的控制电压Ud、Uq,输出α-β轴坐标系下的控制电压Uα、Uβ;SVPWM模块(即空间矢量脉宽调制模块)输入控制电压Uα、Uβ,输出脉冲信号至三相逆变器,三相逆变器将直流电压Udc转换为三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转。
忽略铁心涡流与磁滞损耗等影响,采用id=0的PMSM转子磁场定向控制,建立PMSM在d-q轴旋转坐标系下的数学模型,电压方程为:
Figure BDA0002665884180000051
对于采用id=0控制方式的凸极式PMSM矢量控制系统,电磁转矩方程为:
Figure BDA0002665884180000052
PMSM运动方程为:
Figure BDA0002665884180000053
式(1)(2)(3)中,ud、uq分别是d-q轴的电压;id、iq分别是d-q轴的电流;Ld、Lq分别是d-q轴的电感;Te是电磁转矩;TL是负载转矩;R是定子的电阻;p是电机极对数;ωe是转子电角速度,即角频率;ω是转子角速度,即电机转子机械角速度;ψf是永磁体磁链;J是转动惯量;B是摩擦系数;t是时间。
令电机的转子角速度误差e=ω*-ω,ω*是电机的给定转子角速度。定义永磁同步电机速度控制系统实施例1的状态变量为:
Figure BDA0002665884180000054
由式(2)、(3)、(4)得:
Figure BDA0002665884180000061
对式(5)进行化简,令D=1.5pψf/J,
Figure BDA0002665884180000062
可得实施例1系统状态空间方程为:
Figure BDA0002665884180000063
选择滑模面函数为:
s=cx1+x2 (7)
式(7)中,s为滑模面,c为滑模面参数,且c>0。式(7)中c为转子角速度误差项的系数,其大小对控制作用的影响主要类似于PID控制中的比例系数,c的取值同时还兼顾平衡转子角速度误差和转子角速度误差变化率,通常情况下,在大于0且小于1000的范围内选取,例如,取c=60。对式(7)求导可得:
Figure BDA0002665884180000064
传统指数趋近律的表达式为:
Figure BDA0002665884180000065
式(9)中,sgn()为符号函数,-k1sgn(s)为等速趋近项,-k2s为指数趋近项,k1、k2两个系数分别决定滑模面的抖振和趋近过程的运动品质,且k1、k2均大于0。为提高系统响应速度,在传统指数趋近率的基础上改进,等速趋近项改变为变速趋近项,改进的趋近律为:
Figure BDA0002665884180000066
其中k1>0,k2>0,1<k3<2,k4>0。当|s|较大时,
Figure BDA0002665884180000067
变速趋近项趋近速度较大,可加快滑模趋近运动速度;当|s|较小时,
Figure BDA0002665884180000068
变速趋近项趋近速度较小,可减弱抖振。k4取值可以参照永磁同步电机额定负载启动时的转子角速度变化率,在该转子角速度变化率附近选取,进一步地,在该转子角速度变化率的80%至120%范围内取值;例如,设某永磁同步电机额定负载启动至额定转速1500r/min所用时间为0.21s,平均的转子角速度变化率是750rad/s2,建议此时k4在600~900范围内取值。k3越大,则变速越大,k3一般在1.05~1.3范围内取值。一般地,系数k1和系数k2的取值均小于2000;系数k2越大,系统状态能以较大速度趋近于滑动模态;系数k1决定到达切换面的速度,k1越小则穿越切换面的距离和抖动越小。k1和k2分别为变速趋近项系数和指数趋近项系数,由于
Figure BDA0002665884180000078
的值在1附近变化,因此,式(10)中变速趋近项系数k1和指数趋近项系数k2可以按照调整传统指数趋近率中等速趋近项系数和指数趋近项系数的方法进行整定。k3为变速系数,其大小改变变速快慢;k4是迁移系数,其大小改变变速临界点。
结合式(8)和(10),将计算得出的q轴给定电流作为转矩电流给定分量i′q,可得滑模速度控制器输出的负载转矩给定值TL *和转矩电流给定分量i′q为:
Figure BDA0002665884180000071
永磁同步电机速度控制系统实施例1滑模速度控制器输出中含有积分项,对控制量进行了滤波,可削弱系统抖振、减少系统的稳态误差。定义Lyapunov函数为:
Figure BDA0002665884180000072
由式(10)和式(12)可得:
Figure BDA0002665884180000073
式(13)中,k1>0,k2>0,s·sgn(s)≥0,
Figure BDA0002665884180000074
Figure BDA0002665884180000075
表明系统跟踪误差能在有限时间内收敛到零,系统能稳定运行。
设计滑模速度控制器时整定参数c、k1、k2、k3、k4的方法是,首先确定k3、k4的值;令q轴转矩电流给定值i* q只包括输入转矩电流给定分量i′q(即不进行负载转矩补偿控制),然后在系统的滑动模态从小到大调整滑模面参数c和变速趋近项系数k1的值,直到系统出现明显抖振,在此基础上兼顾抑制抖振和系统状态收敛速度,适当减小滑模面参数c和变速趋近项系数k1的值;最后,在兼顾抑制滑动模态抖振的情况下,主要依据系统到达段(例如,阶跃响应的电机启动阶段)的快速性来调整指数趋近项系数k2的值,并对滑模速度控制器的其他参数值进行适当微调。
根据PMSM电磁转矩与运动方程,对于恒定的阶跃性负载,在变化周期内可认为是一恒定值,即
Figure BDA0002665884180000076
将电机转子角速度与负载转矩作为状态变量,构成PMSM状态方程为:
Figure BDA0002665884180000077
在式(14)基础上,以负载转矩与电机转子角速度为观测对象,建立负载转矩观测器实施例1为:
Figure BDA0002665884180000081
式(15)中,
Figure BDA0002665884180000082
是负载转矩的观测值,
Figure BDA0002665884180000083
是转子角速度估计值,g是负载转矩观测器的反馈增益,
Figure BDA0002665884180000084
kg是负载转矩观测器实施例1的滑模增益,负载转矩观测器实施例1为滑模观测器。电机摩擦与负载转矩相比较,所占比重小,令B=0,忽略摩擦的影响,则式(15)的负载转矩观测器实施例1变为:
Figure BDA0002665884180000085
依据B=0时的(14)和式(16),得到负载转矩观测器实施例1的误差方程为:
Figure BDA0002665884180000086
式(17)中,
Figure BDA00026658841800000817
为转子角速度估算误差,
Figure BDA0002665884180000087
为负载转矩观测误差,并定义观测器滑模面为
Figure BDA0002665884180000088
根据滑模可达性条件,式(16)观测器系统稳定条件为kg≤-|e2/J|,且g<0。
在式(14)基础上,以负载转矩与电机转子角速度为观测对象,还可以建立负载转矩观测器实施例2为:
Figure BDA0002665884180000089
电机摩擦与负载转矩相比较,所占比重小,令B=0,忽略摩擦的影响,则式(18)的负载转矩观测器实施例2变为:
Figure BDA00026658841800000810
式(18)、(19)中,
Figure BDA00026658841800000811
是负载转矩的观测值,
Figure BDA00026658841800000812
是转子角速度估计值,g是负载转矩观测器的反馈增益,
Figure BDA00026658841800000813
kW是负载转矩观测器实施例2的比例增益,负载转矩观测器实施例2为状态观测器。依据B=0时的式(14)和式(19),得到负载转矩观测器实施例2误差方程为:
Figure BDA00026658841800000814
式(20)中,
Figure BDA00026658841800000815
为转子角速度估算误差,
Figure BDA00026658841800000816
为负载转矩观测误差。式(19)的状态观测器是一个自治的线性系统,在kW<0,且g<0时,该观测器是渐近稳定的。负载转矩观测器实施例1的式(15)和负载转矩观测器实施例2的式(18)均考虑了电机的摩擦因素,小摩擦阻尼的加入,会给系统响应的快速性带来不利影响,但可以使稳定性分别在式(16)和式(19)的基础上增加。
选择式(15)、(16)的观测器实施例1时,滑模增益kg的整定方法是,按照
Figure BDA0002665884180000091
进行选择。式(21)中,α≥1;一般地,α值在1~5的范围内选择,例如,选择α等于1.5。负载转矩观测器实施例1在观测负载转矩的过程中,kg的绝对值选择过小,当|e2|较大时观测器不能进入滑模状态;kg的绝对值选择足够大,可保证观测器进入滑模状态,但负载转矩的稳态观测波动变大;kg的值随着负载转矩观测误差的变化而变化,可以同时兼顾增加观测器稳定性和减小负载转矩的稳态观测波动。
选择式(18)、(19)的观测器实施例2时,比例增益kW的整定方法是,按照
Figure BDA0002665884180000092
进行选择。式(22)中,TN是电机的额定转矩,β>0;一般地,β值在1~20的范围内选择,例如,选择β=10。β选择增大时,负载转矩观测的稳态波动变大,但转矩观测跟踪超调量变小;β选择减小时,负载转矩观测的稳态波动变小,但转矩观测跟踪超调量变大。
在式(15)、(16)或者式(18)、(19)所表示的观测器中,反馈增益g的取值大小对负载转矩观测结果影响较大。反馈增益g越大,观测转矩波动性越小,但观测转矩辨识速度越慢;反馈增益g越小,观测转矩速度越快,但观测转矩波动性越大。出于对此问题的考虑,在传统负载转矩观测器中,综合考虑负载转矩的观测速度与波动性,将反馈增益g取一个折中值,但这会舍弃大反馈增益时的波动性小和小反馈增益时的观测速度快的优势。
电机滑模速度控制主要通过增大控制器中不连续项的幅值来抑制参数变化和外部负载扰动对系统的影响,但幅值增大会引起滑模固有抖振。为解决滑模控制系统抖振与抗扰性之间的矛盾,利用观测器实时观测负载扰动变化,将负载转矩观测值前馈补偿至电流调节器中,以降低滑模控制中不连续项幅值,削弱参数变化导致给定转矩变化,或者是负载扰动所引起系统抖振。为了充分利用反馈增益g在高、低值时的优势,根据两相邻时刻负载转矩观测值以及负载转矩给定值变化量的大小,在负载转矩给定值变化小且负载转矩观测值变化小时,给予反馈增益g较大值,使观测结果波动性小,稳定性更强;在负载转矩给定值变化大或者负载转矩观测值变化大时,给予反馈增益g较小值,使观测速度加快,最终通过对反馈增益g的调整,得到观测速度快和波动小、稳定性更强的综合结果。
负载转矩观测器实施例1,或者是负载转矩观测器实施例2用于图1永磁同步电机速度控制系统实施例1中时,负载转矩观测器根据负载转矩给定值
Figure BDA00026658841800001019
和负载转矩观测值
Figure BDA00026658841800001020
的变化对反馈增益g进行调整,依据转子角速度ω和电流iq对负载转矩进行观测,得到新的负载转矩观测值
Figure BDA00026658841800001018
图2为反馈增益自动调整方法实施例1流程图,在负载转矩观测器实施例1,或者是负载转矩观测器实施例2用于图1永磁同步电机速度控制系统实施例1中时,进行反馈增益自动调整。在一次电机速度的周期性控制过程中,图2中的(b)所示反馈增益g的调整晚于负载转矩观测和滑模速度控制器的输出计算,有:
步骤①、负载转矩观测器依据现有反馈增益g的值对负载转矩TL进行观测,得到负载转矩观测值
Figure BDA0002665884180000101
滑模速度控制器进行控制运算得到负载转矩给定值
Figure BDA0002665884180000102
此时的
Figure BDA0002665884180000103
Figure BDA0002665884180000104
Figure BDA0002665884180000105
Figure BDA0002665884180000106
至下一次进行反馈增益g调整时,该次
Figure BDA0002665884180000107
变为
Figure BDA0002665884180000108
变为
Figure BDA0002665884180000109
步骤②、计算
Figure BDA00026658841800001010
Figure BDA00026658841800001011
步骤③、判断ΔT是否大于ε2;当ΔT大于ε2时,取反馈增益g等于gmin并退出;当ΔT小于等于ε2时,进入步骤④;
步骤④、判断ΔT是否小于ε1;当ΔT小于ε1时,取反馈增益g等于gmax并退出;当ΔT大于等于ε1时,进入步骤⑤;
步骤⑤、反馈增益g按照
Figure BDA00026658841800001012
进行计算。
其中,ε1为转矩变化低限比较阈值,ε2为转矩变化高限比较阈值,且0<ε1<ε2;gmax为反馈增益高值,gmin为反馈增益低值,且gmin<gmax<0。
在一次电机速度的周期性控制过程中,图2中的(a)所示反馈增益g的调整先于负载转矩观测和滑模速度控制器的输出计算,其反馈增益g调整方法将上述步骤①变为步骤⑤,步骤②-⑤变为步骤①-④,各步骤中的退出改为进入步骤⑤,且
Figure BDA00026658841800001013
Figure BDA00026658841800001014
图2中,最近2次负载转矩给定值的变化量与负载转矩观测值的变化量之和
Figure BDA00026658841800001015
Figure BDA00026658841800001016
ΔTL *为最近2次负载转矩给定值之间的差值,
Figure BDA00026658841800001017
为最近2次负载转矩观测值之间的差值。当ΔT大于ε2时,表明负载转矩观测值波动大,或者是因系统模型参数变化、转子角速度给定值变化、转子角速度实际值变化,使负载转矩给定值变化大并将引起负载转矩观测值有较大波动,反馈增益g选择等于gmin进行负载转矩的快速辨识与观测;当ΔT小于ε1时,表明负载转矩给定值变化小且状态负载转矩观测值波动小,反馈增益g选择等于gmax进行以稳定为主的负载转矩辨识与观测;当ΔT大于等于ε1且小于等于ε2时,反馈增益g按照式(23)进行计算,使反馈增益g在此区间内随着ΔT的增加而减小,避免因ΔT的较小变化而使反馈增益g产生剧烈变化带来的对转矩观测器工作稳定性的不利影响。图2中,ε1、ε2的具体取值与滑模速度控制器的采样控制周期(周期时间)、永磁同步电机及其负载情况相关,ε2一般在小于额定转矩5%的范围内取值,例如,额定转矩为22N·m,取ε1=0.1N·m,ε2=0.6N·m。反馈增益g的取值满足gmin<gmax<0,一般情况下,gmin≥-5000。gmin取值应在负载转矩突变时,负载转矩观测器的转矩观测跟踪超调量处于转矩观测跟踪超调限值之内;gmax取值应在负载转矩不变,负载转矩观测器和滑模速度控制器均处于稳态时,最近2次负载转矩给定值的变化量与负载转矩观测值的变化量之和ΔT小于ε1;例如,选择反馈增益gmax=-0.5,gmin=-10。选取gmin、gmax、ε1、ε2值的具体方法是:
步骤⑴,负载转矩观测器和滑模速度控制器均处于稳态且保持给定转子角速度不变和负载转矩不变;
步骤⑵,令反馈增益g从一个较大值开始逐渐减小,例如,令反馈增益g从-0.01开始逐渐减小,当负载转矩观测的稳态误差达到负载转矩观测稳态误差限值时,确定此时的反馈增益g值为gmax
步骤⑶,保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmax,连续进行n次ΔT值的测量,并将此时n次测量中最大m个ΔT值之和的平均值作为转矩变化低限比较阈值ε1
步骤⑷,微调改变反馈增益g,负载转矩观测器和滑模速度控制器均处于稳态时,保持给定转子角速度不变和令负载转矩突变,在保证负载转矩观测器输出观测值的转矩观测跟踪超调量处于转矩观测跟踪超调限值之内的前提下,测量负载转矩观测器的跟踪调节时间;
步骤⑸,重复步骤⑷,选择跟踪调节时间最短的反馈增益g值为gmin值;通常情况下,转矩观测跟踪超调量接近转矩观测跟踪超调限值时,负载转矩观测器的跟踪调节时间较短;
步骤⑹,再次保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmin,连续进行n次ΔT值的测量,并将此时n次测量中最大m个ΔT值之和的平均值作为转矩变化高限比较阈值ε2
观测得到负载转矩观测值
Figure BDA0002665884180000111
后,将负载转矩的观测值
Figure BDA0002665884180000112
转换成转矩电流补偿分量i″q前馈补偿至q轴电流PI控制器的输入,对滑模速度控制器输出的转矩电流给定分量i′q进行补偿。q轴电流PI控制器的q轴转矩电流给定值i* q为:
Figure BDA0002665884180000121
式(24)中,kq=1/(1.5pψf)为转矩观测补偿系数。比较式(11)与式(24)可得,当负载发生扰动或者是系统参数发生变化时,式(11)中没有加入负载转矩补偿,需要选取较大的k1、k2值来提供足够大的给定电流变化量来抵消负载发生扰动或者是系统参数发生变化的相关影响,以保证电机转速能够快速恒定;而式(24)将负载转矩观测值前馈补偿至电流调节器中,在不需要较大k1、k2值的情况下就能在负载发生扰动或者是系统参数发生变化时,提供足够大的给定电流变化量来抵消负载发生扰动或者是系统参数发生变化的相关影响,减少滑模速度控制器的输出压力及不连续项的幅值,有效地削弱了系统的抖振。
固定反馈增益值时,反馈增益g越小,负载转矩观测的振荡幅值越大,波动性越强;反馈增益g越大,负载转矩观测的振荡幅值越小,观测精度越高。自动调整增益算法解决了负载转矩观测器中小反馈增益导致转矩观测波动大的问题,大反馈增益收敛时间长的问题,收敛时间和波动幅度指标都优于折中增益算法,能快速跟踪负载转矩变化值以及快速降低给定变化或者参数变化所带来的观测误差,且振荡幅值小,观测精度高,达到了较好的观测效果。
在额定负载转矩下改变给定转速时,尽管实际的负载转矩没有改变,但从式(15)、(16)或者式(18)、(19)所构造的负载转矩观测器可知,当转子角速度ω发生改变时,即使负载转矩没有改变,观测到的转矩观测值也会发生变化,造成观测误差。在额定负载转矩下改变给定转速时,永磁同步电机滑模控制系统的控制调节过程是,首先由滑模速度控制器依据给定速度变化,使输出的负载转矩给定值TL *改变,使转矩电流给定值i* q发生改变,进而使永磁同步电机的电磁转矩Te发生改变,带动电机使转子角速度ω发生改变;反馈增益g如果只依据负载转矩观测值的变化量
Figure BDA0002665884180000122
进行自动调整,则此时只有当转子角速度ω发生改变,使负载转矩观测值
Figure BDA0002665884180000123
发生了改变后,才对反馈增益g进行调整;反馈增益g依据负载转矩给定值的变化量ΔTL *和负载转矩观测值的变化量
Figure BDA0002665884180000124
的绝对值之和进行自动调整,当给定转速改变使负载转矩给定值TL *改变,负载转矩观测值
Figure BDA0002665884180000125
还没有发生变化时,就提前调整反馈增益g,当负载转矩观测值
Figure BDA0002665884180000126
真正产生观测误差时,能够加快观测器响应速度,尽快消除(减小)负载转矩观测值
Figure BDA0002665884180000127
的观测误差,进而改善电机速度控制的快速性与精确性。同样地,当系统模型参数发生变化使负载转矩给定值TL *先于负载转矩观测值
Figure BDA0002665884180000128
发生改变时,反馈增益g同时依据负载转矩给定值的变化量ΔTL *和负载转矩观测值的变化量
Figure BDA0002665884180000129
进行自动调整,能够提前调整反馈增益g,加快观测器响应速度,尽快消除(减小)负载转矩观测值
Figure BDA00026658841800001210
的观测误差,并进一步改善电机速度控制的快速性与精确性。当然,如果负载发生扰动导致观测值
Figure BDA00026658841800001211
发生改变时,
Figure BDA00026658841800001212
发生较大变化时,从图2可知,反馈增益g也能够进行自动调整,以尽快消除(减小)负载转矩观测值
Figure BDA0002665884180000131
的观测误差,使负载转矩观测值
Figure BDA0002665884180000132
尽快跟上负载转矩TL的变化。
在永磁同步电机速度控制系统实施例1的周期性控制过程中,将当次k时刻(或者是第k步)计算得到的负载转矩给定值TL *记为TL *(k),负载转矩观测值
Figure BDA0002665884180000133
记为
Figure BDA0002665884180000134
k-1时刻为k时刻的前一次周期性控制过程时刻,负载转矩给定值TL *记为TL *(k-1),负载转矩观测值
Figure BDA0002665884180000135
记为
Figure BDA0002665884180000136
k-2时刻为k-1时刻的前一次周期性控制过程时刻,负载转矩给定值TL *记为TL *(k-2),负载转矩观测值
Figure BDA0002665884180000137
记为
Figure BDA0002665884180000138
图2中的(b)先进行负载转矩观测和速度控制,后进行反馈增益自动调整,此时电机速度的周期性控制过程如下:
步骤一、检测永磁同步电机的转子位置θ、转子角速度ω和三相电流ia、ib和ic
步骤二、依据三相电流ia、ib和ic对永磁同步电机进行Clark变换得到在α-β轴坐标系下的电流iα、iβ,依据电流iα、iβ和转子位置θ进行Park变换,得到在d-q轴坐标系下的电流id、iq
步骤三、负载转矩观测器依据转子角速度ω和电流iq对负载转矩进行观测,得到负载转矩观测值
Figure BDA0002665884180000139
和转矩电流补偿分量i″q
步骤四、滑模速度控制器依据输入的转子给定角速度ω*和转子角速度ω进行控制计算,得到负载转矩给定值
Figure BDA00026658841800001310
和转矩电流给定分量i′q
步骤五、负载转矩观测器的反馈增益g依据负载转矩给定值TL *和负载转矩观测值
Figure BDA00026658841800001311
的变化进行调整;
步骤六、依据转矩电流给定分量i′q和转矩电流补偿分量i″q计算得到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β;d轴转矩电流给定值id *等于0;
步骤七、将α-β轴坐标系下的控制电压Uα、Uβ作为SVPWM模块的输入,由SVPWM模块控制三相逆变器产生三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转。
图2中的(a)先进行反馈增益自动调整,后进行负载转矩观测和速度控制,上述步骤中,应该将步骤五与步骤三、四的先后顺序互换,即先进行步骤五,后进行步骤三、四。
图3为实现基于负载转矩观测器的永磁同步电机滑模控制方法的永磁同步电机速度控制系统实施例2的框图。图3实施例2与图1实施例1的区别在于,滑模速度控制器采用积分滑模控制方式,负载转矩观测器的观测值
Figure BDA00026658841800001312
被送至滑模速度控制器,滑模速度控制器输出的q轴给定电流(转矩电流给定分量)中中已经包括有负载转矩观测值
Figure BDA00026658841800001313
的计算项(即补偿分量),因此,实施例2中滑模速度控制器输出的q轴给定电流(转矩电流给定分量)直接作为q轴转矩电流给定值,同样能够起到负载转矩补偿的作用;滑模速度控制器输出的负载转矩给定值TL Δ中同样已经包括有负载转矩观测值
Figure BDA0002665884180000141
的计算项,负载转矩观测器直接根据负载转矩给定值TL Δ的变化量进行反馈增益自动调整,其作用与前述反馈增益自动调整方法实施例1依据最近2次负载转矩给定值的变化量与负载转矩观测值的变化量之和ΔT进行反馈增益自动调整相同;
Figure BDA0002665884180000142
定义永磁同步电机速度控制系统实施例2的状态变量为:
Figure BDA0002665884180000143
选择滑模面函数为:
sy=cyy1+y2 (26)
式(26)中,cy为滑模面参数,且cy>0。式(26)中cy为转子角速度误差积分项的系数,其大小对控制作用的影响主要类似于PID控制中的比例系数,cy的取值同时还兼顾平衡转子角速度误差积分项和转子角速度误差项,通常情况下,cy在大于0且小于100的范围内选取。对式(26)求导可得:
Figure BDA0002665884180000144
在传统指数趋近律的基础上,采用新的趋近律为:
Figure BDA0002665884180000145
μ1、μ2、μ3、μ4为速度滑模控制的指数趋近率系数,其中,μ1>0,μ2>0,0<μ3<1,μ4>0。当|sy|较大时,
Figure BDA0002665884180000146
变速趋近项趋近速度较大,可加快滑模趋近运动速度;当|sy|较小时,
Figure BDA0002665884180000147
变速趋近项趋近速度较小,可减弱抖振。μ4取值可以参照永磁同步电机稳定运行时转子角速度的允许稳态误差,μ4取值不大于该允许稳态误差的倒数值,进一步地,在该倒数值的50%至100%的范围内取值;例如,设允许的永磁同步电机转子角速度稳态误差是5rad/s(弧度/秒),μ4建议取值不大于0.2,进一步地,μ4在0.1至0.2范围内取值。μ3一般在0.5附近取值,进一步地,μ3一般在0.4至0.6范围内取值。一般地,系数μ1和系数μ2的取值均小于5000。μ1和μ2分别为变速趋近项系数和指数趋近项系数,由于
Figure BDA0002665884180000148
的值在1附近变化,因此,式(28)中变速趋近项系数μ1和指数趋近项系数μ2可以按照调整传统指数趋近率中等速趋近项系数和指数趋近项系数的方法进行整定。μ3是迁移系数,其大小改变变速临界点;μ4为变速系数,其大小改变变速快慢。式(28)中的e是自然指数,即自然对数的底数。
结合式(2)、(3)、(27),得:
Figure BDA0002665884180000151
结合式(28)、(29),将计算得出的q轴给定电流直接作为q轴转矩电流给定值iΔ q,可得控制器输出的q轴转矩电流给定值iΔ q和负载转矩给定值TL Δ为:
Figure BDA0002665884180000152
式(30)中,负载转矩值TL使用负载转矩观测器的输出值
Figure BDA0002665884180000153
来代替。定义Lyapunov函数为:
Figure BDA0002665884180000154
由式(26)和式(28)可得:
Figure BDA0002665884180000155
式(32)中,μ1>0,μ2>0,
Figure BDA0002665884180000156
sy·sgn(sy)≥0,故
Figure BDA0002665884180000157
表明观测器的跟踪误差能在有限时间内收敛到零,系统能稳定运行。
设计滑模速度控制器时整定参数cy、μ1、μ2、μ3、μ4的方法是,首先确定μ3、μ4的值;令式(30)中的负载转矩观测器的输出值
Figure BDA0002665884180000158
(即不进行负载转矩补偿控制),然后在系统的滑动模态从小到大调整滑模面参数cy、变速趋近项系数μ1的值,直到系统出现明显抖振,在此基础上兼顾抑制抖振和系统状态收敛速度,适当减小滑模面参数cy和变速趋近项系数μ1的值;最后,在兼顾抑制滑动模态抖振的情况下,主要依据系统到达段(例如,阶跃响应的电机启动阶段)的快速性来调整指数趋近项系数μ2的值,并对滑模速度控制器的其他参数值进行适当微调。
图3永磁同步电机速度控制系统实施例2中的负载转矩观测器仍然采用前述的负载转矩观测器实施例1,或者是采用前述的负载转矩观测器实施例2;此时,负载转矩观测器根据滑模速度控制器输出的负载转矩给定值
Figure BDA0002665884180000159
的变化对反馈增益g进行调整,依据转子角速度ω和电流iq对负载转矩TL进行观测,得到负载转矩观测值
Figure BDA00026658841800001510
图4为反馈增益自动调整方法实施例2流程图,在负载转矩观测器实施例1,或者是负载转矩观测器实施例2用于图3永磁同步电机速度控制系统实施例2中时,进行反馈增益自动调整。图4中,ε1为转矩变化低限比较阈值,ε2为转矩变化高限比较阈值,且0<ε1<ε2;gmax为反馈增益高值,gmin为反馈增益低值,且gmin<gmax<0;ΔTL Δ为最近2次负载转矩给定值之间的差值。在一次电机速度的周期性控制过程中,图4中的(a)所示反馈增益g的调整先于负载转矩观测和滑模速度控制器的输出计算,具体方法是:
步骤㈠、计算
Figure BDA0002665884180000161
步骤㈡、判断
Figure BDA0002665884180000162
是否大于ε2;当
Figure BDA0002665884180000163
大于ε2时,取反馈增益g等于gmin并进入步骤㈤;当ΔT小于等于ε2时,进入步骤㈢;
步骤㈢、判断
Figure BDA0002665884180000164
是否小于ε1;当
Figure BDA0002665884180000165
小于ε1时,取反馈增益g等于gmax并进入步骤㈤;当
Figure BDA0002665884180000166
大于等于ε1时,进入步骤㈣;
步骤㈣、反馈增益g按照
Figure BDA0002665884180000167
进行计算后进入步骤㈤;
步骤㈤、负载转矩观测器对负载转矩TL进行观测,得到负载转矩观测值
Figure BDA0002665884180000168
滑模速度控制器进行控制运算得到
Figure BDA0002665884180000169
此时的
Figure BDA00026658841800001610
Figure BDA00026658841800001611
在下一次电机速度的周期性控制过程中变为
Figure BDA00026658841800001612
在一次电机速度的周期性控制过程中,图4中的(b)所示反馈增益g的调整晚于负载转矩观测和滑模速度控制器的输出计算,具体方法是:
步骤A、负载转矩观测器对负载转矩TL进行观测,得到负载转矩观测值
Figure BDA00026658841800001613
滑模速度控制器进行控制运算得到
Figure BDA00026658841800001614
此时的
Figure BDA00026658841800001615
Figure BDA00026658841800001616
在下一次电机速度的周期性控制过程中变为
Figure BDA00026658841800001617
步骤B、计算
Figure BDA00026658841800001618
步骤C、判断
Figure BDA00026658841800001619
是否大于ε2;当
Figure BDA00026658841800001620
大于ε2时,取反馈增益g等于gmin并退出;当ΔT小于等于ε2时,进入步骤D;
步骤D、判断
Figure BDA00026658841800001621
是否小于ε1;当
Figure BDA00026658841800001622
小于ε1时,取反馈增益g等于gmax并退出;当
Figure BDA00026658841800001623
大于等于ε1时,进入步骤E;
步骤E、反馈增益g按照式(33)进行计算后退出。
TL Δ的输出项中,既包括因系统参数变化,或者是转子角速度给定值变化,或者是转子角速度实际值变化而处于变化状态的给定分项
Figure BDA00026658841800001624
也包括补偿分项负载转矩观测值
Figure BDA00026658841800001625
当|ΔTL Δ|大于ε2时,表明负载转矩观测值波动大,或者是因系统模型参数变化、转子角速度给定值变化、转子角速度实际值变化,使TL Δ中的给定分项变化大并将引起负载转矩观测值有较大波动,反馈增益g选择等于gmin进行负载转矩的快速辨识与观测;当|ΔTL Δ|小于ε1时,表明负载转矩观测值波动小且将引起负载转矩观测值有较大波动的因素(即TL Δ中的给定分项)变化小,反馈增益g选择等于gmax进行以稳定为主的负载转矩辨识与观测;当|ΔTL Δ|大于等于ε1且小于等于ε2时,反馈增益g按照式(33)进行计算,使反馈增益g在此区间内随着|ΔTL Δ|的增加而减小,避免因|ΔTL Δ|的较小变化而使反馈增益g产生剧烈变化带来的对转矩观测器工作稳定性的不利影响。图4中,ε1、ε2的具体取值与滑模速度控制器的采样控制周期(周期时间)、永磁同步电机及其负载情况相关,ε2一般在小于额定转矩5%的范围内取值,例如,额定转矩为22N·m,取ε1=0.1N·m,ε2=0.6N·m。反馈增益g的取值满足gmin<gmax<0,一般情况下,gmin≥-5000。gmin取值应在负载转矩突变时,负载转矩观测器输出观测值的转矩观测跟踪超调量处于转矩观测跟踪超调限值之内;gmax取值应在负载转矩不变,负载转矩观测器和滑模速度控制器均处于稳态时,最近2次负载转矩给定值之间的差值|ΔTL Δ|小于ε1;例如,选择反馈增益gmax=-0.5,gmin=-10。
选取gmin、gmax、ε1、ε2值的具体方法是:
步骤1),负载转矩观测器和滑模速度控制器均处于稳态且保持给定转子角速度不变和负载转矩不变;
步骤2),令反馈增益g从一个较大值开始逐渐减小,例如,令反馈增益g从-0.01开始逐渐减小,当负载转矩观测的稳态误差达到负载转矩观测稳态误差限值时,确定此时的反馈增益g值为gmax
步骤3),保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmax,连续进行n次
Figure BDA0002665884180000171
值的测量,并将n次
Figure BDA0002665884180000172
测量中最大m个
Figure BDA0002665884180000173
测量值的平均值作为转矩变化比较阈值ε;
步骤4),微调改变反馈增益g,负载转矩观测器和滑模速度控制器均处于稳态时,保持给定转子角速度不变和令负载转矩突变,在保证负载转矩观测器输出观测值的转矩观测跟踪超调量处于转矩观测跟踪超调限值之内的前提下,测量负载转矩观测器的跟踪调节时间;
步骤5),重复步骤4),选择跟踪调节时间最短的反馈增益g值为gmin值;通常情况下,转矩观测跟踪超调量接近转矩观测跟踪超调限值时,负载转矩观测器的跟踪调节时间较短;
步骤6),再次保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmin,连续进行n次
Figure BDA0002665884180000174
值的测量,并将此时n次测量中最大m个
Figure BDA0002665884180000175
值之和的平均值作为转矩变化高限比较阈值ε2
在以上各选取gmin、gmax值和比较阈值ε1、ε2的具体方法中,均在滑模速度控制器中的参数已经完成整定,且在进行负载转矩补偿控制的情况下实现;建议n为大于等于20的整数,m为大于等于5且小于等于0.5n的整数。
式(30)滑模速度控制器的输出项中包括补偿分项负载转矩观测值
Figure BDA0002665884180000181
相当于式(24)中将负载转矩观测值前馈至电流调节器的给定值中,在不需要滑模速度控制器输出的给定电流部分产生较大调整的情况下,就能抵消负载发生扰动或者是系统参数发生变化所造成的相关影响,有效地削弱了系统的抖振。反馈增益g依据负载转矩给定值的变化量ΔTL Δ进行自动调整的算法,避免了负载转矩观测器选择固定小反馈增益导致转矩观测波动大,选择固定大反馈增益导致收敛时间长的问题,能在系统的控制参数、模型参数等发生变化或者是负载发生扰动,导致负载转矩给定值中的给定分项部分发生变化或/和负载转矩观测值部分发生变化时,快速降低负载转矩的观测误差,改善了观测效果和电机速度控制的快速性与精确性。反馈增益g在负载转矩给定值变化时即进行自动调整,能够在负载转矩观测值还没有发生较大变化,但因转子角速度给定值改变或/和转子角速度实际值改变使负载转矩给定值中的给定分项部分改变,或者是因系统模型参数发生变化使负载转矩给定值中的给定分项部分发生改变,将引起负载转矩观测值有较大波动时,提前调整反馈增益g,当负载转矩观测值真正产生观测误差时,加快观测器的响应速度,快速降低负载转矩观测值的观测误差,并进一步改善电机速度控制的快速性与精确性。
滑模速度控制器采用积分滑模控制方式的永磁同步电机速度控制系统实施例2,其反馈增益依据负载转矩给定值TL Δ的变化量进行自动调整的算法,作用与永磁同步电机速度控制系统实施例1中,反馈增益依据最近2次负载转矩给定值的变化量与负载转矩观测值的变化量之和
Figure BDA0002665884180000182
进行自动调整的算法相同,避免了负载转矩观测器选择固定小反馈增益导致转矩观测波动大,选择固定大反馈增益导致收敛时间长的问题,能在系统的控制参数、模型参数等发生变化或者是负载发生扰动,导致负载转矩给定值TL Δ发生变化(包括给定分项发生变化或/和补偿分项发生变化)时,快速降低负载转矩的观测误差,改善了观测效果和电机速度控制的快速性与精确性。反馈增益g依据负载转矩给定值变化TL Δ进行自动调整,能够在负载转矩观测值
Figure BDA0002665884180000183
还没有发生较大变化,但因转子角速度给定值变化或/和转子角速度实际值变化使负载转矩给定值的给定分项发生变化,或者是系统模型参数发生变化使负载转矩给定值的给定分项发生变化,将引起负载转矩观测值有较大波动时,提前调整反馈增益g,当负载转矩观测值
Figure BDA0002665884180000184
真正产生观测误差时,加快观测器的响应速度,快速降低负载转矩观测值
Figure BDA0002665884180000185
的观测误差,并进一步改善电机速度控制的快速性与精确性。
在永磁同步电机速度控制系统实施例2的周期性控制过程中,将当次k时刻(或者是第k步)计算得到的负载转矩给定值
Figure BDA0002665884180000186
记为
Figure BDA0002665884180000187
负载转矩观测值
Figure BDA0002665884180000188
记为
Figure BDA0002665884180000189
k-1时刻为k时刻的前一次周期性控制过程时刻,负载转矩给定值
Figure BDA00026658841800001810
记为
Figure BDA00026658841800001811
负载转矩观测值
Figure BDA0002665884180000191
记为
Figure BDA0002665884180000192
k-2时刻为k-1时刻的前一次周期性控制过程时刻,负载转矩给定值
Figure BDA0002665884180000193
记为
Figure BDA0002665884180000194
负载转矩观测值
Figure BDA0002665884180000195
记为
Figure BDA0002665884180000196
按照图4中的(a)进行反馈增益g调整时,控制过程包括以下步骤:
步骤一、检测永磁同步电机的转子位置θ、转子角速度ω和三相电流ia、ib和ic
步骤二、依据三相电流ia、ib和ic对永磁同步电机进行Clark变换得到在α-β轴坐标系下的电流iα、iβ,依据电流iα、iβ和转子位置θ进行Park变换,得到在d-q轴坐标系下的电流id、iq
步骤三、负载转矩观测器的反馈增益g依据负载转矩给定值
Figure BDA0002665884180000197
的变化进行调整;
步骤四、负载转矩观测器依据转子角速度ω和电流iq对负载转矩进行观测,得到负载转矩观测值
Figure BDA0002665884180000198
步骤五、滑模速度控制器依据输入的转子给定角速度ω*、转子角速度ω和负载转矩观测值
Figure BDA0002665884180000199
进行控制计算,得到负载转矩给定值
Figure BDA00026658841800001910
和q轴转矩电流给定值
Figure BDA00026658841800001911
步骤六、d轴电流控制器依据d轴转矩电流给定值id *与d轴坐标系下的电流id之间的差值进行PI控制运算,得到d轴坐标系下的控制电压Ud;q轴电流控制器依据q轴转矩电流给定值
Figure BDA00026658841800001912
与q轴坐标系下的电流iq之间的差值进行PI控制运算,得到q轴坐标系下的控制电压Uq;依据d-q轴坐标系下的控制电压Ud、Uq进行Park逆变换,得到α-β轴坐标系下的控制电压Uα、Uβ
步骤七、将α-β轴坐标系下的控制电压Uα、Uβ作为SVPWM模块的输入,由SVPWM模块控制三相逆变器产生三相交流电源Ua、Ub、Uc,从而驱动永磁同步电机运转。
按照图4中的(b)进行反馈增益g调整时,上述控制过程的步骤中,步骤四、五的内容先进行,步骤三的内容后进行。
在以上的各实施例中,转矩观测跟踪超调限值通常为电机额定转矩的1%~10%,具体来说,转矩观测跟踪超调限值是额定转矩的2%,或者是额定转矩的5%,或者是额定转矩的10%,等等。负载转矩从一个定值突变至另外一个定值,突变开始时刻至负载转矩观测器输出观测值稳定地进入负载转矩观测稳态误差限值范围时刻为转矩观测过渡过程,跟踪调节时间指的是该过渡过程时间;负载转矩观测稳态误差指的是负载转矩不变且负载转矩观测器处于稳态时观测转矩瞬时值与负载转矩之间的误差,该误差包括滑模观测器自身抖振引起的观测误差和负载波动之外干扰原因所造成的观测误差,或者是状态观测器因为转子角速度抖振引起的观测误差和负载波动之外干扰原因所造成的观测误差;负载转矩观测稳态误差限值是负载转矩观测器允许的负载转矩观测稳态误差最大绝对值;负载转矩观测稳态误差限值通常为电机额定转矩的1%~5%,具体来说,负载转矩观测稳态误差限值是额定转矩的1%,或者是额定转矩的2%,或者是额定转矩的5%,等等。转矩观测跟踪超调量指的是负载转矩从一个定值突变至另外一个定值,负载转矩观测器输出的观测值超过突变后负载转矩的最大偏差值。当负载转矩观测的稳态误差在负载转矩观测稳态误差限值的临近范围之内,例如,95%~105%范围内,或者是98%~102%范围内时,认为负载转矩观测的稳态误差增大到负载转矩观测稳态误差限值。滑模速度控制器处于稳态指的是滑模速度控制器稳定地处于滑动模态;转子角速度稳态误差指的是稳态时电机转子角速度瞬时值与稳态值之间的差值,转子角速度稳态误差限值为系统允许的转子角速度稳态误差最大绝对值。负载转矩观测器中,实施例1的滑模观测器处于稳态指的是滑模观测器稳定地处于滑动模态;实施例2的状态观测器处于稳态指的是状态观测器处于转矩观测过渡过程之后的工作状态。转子角速度稳态误差指的是稳态时电机转子角速度瞬时值与稳态值之间的差值,转子角速度稳态误差限值为系统允许的转子角速度稳态误差最大绝对值。
本发明的永磁同步电机速度控制系统及永磁同步电机驱动控制方法,除用于矿用牵引电机的控制之外,还可以用于其他永磁同步电机应用场合。
除说明书所述的技术特征外,本发明所涉及的其他技术均为本领域技术人员所掌握的常规技术。例如,q轴电流控制器、d轴电流控制器采用PI控制器进行控制及控制器参数的选择,滑模速度控制器控制参数的选择,位置与速度检测模块使用旋转变压器或光电编码器等进行永磁同步电机转子旋转角度与旋转速度检测,以及对于Clarke变换模块、Park变换模块、Park逆变换模块和SVPWM模块、三相逆变器的变换方法及应用方法,等等,均为本领域技术人员所掌握的常规技术。

Claims (6)

1.一种矿用牵引永磁同步电机驱动控制方法,永磁同步电机速度由滑模速度控制器控制,负载转矩观测器对负载转矩进行观测,负载转矩观测器的输出用于对滑模速度控制器的输出进行负载转矩补偿;其特征在于,负载转矩观测器根据负载转矩给定值和负载转矩观测值的变化对反馈增益进行调整,依据转子角速度ω和电流iq对负载转矩进行观测,得到负载转矩观测值;q轴转矩电流给定值为转矩电流给定分量与转矩电流补偿分量之和;
滑模速度控制器的状态变量为
Figure FDA0002665884170000011
其中,ω是转子角速度,ω*是给定转子角速度;滑模速度控制器的滑模面为s=cx1+x2,c为滑模面参数,且c>0;滑模速度控制器输出的负载转矩给定值
Figure FDA0002665884170000012
和转矩电流给定分量i′q
Figure FDA0002665884170000013
其中,J是转动惯量,p是电机极对数,ψf是永磁体磁链;系数k1、k2、k3、k4为速度滑模控制的指数趋近率系数,且k1>0,k2>0,1<k3<2,k4>0;
负载转矩观测器为
Figure FDA0002665884170000014
其中,
Figure FDA0002665884170000015
是转子角速度估计值,g是负载转矩观测器的反馈增益且g<0;
Figure FDA0002665884170000016
kg是负载转矩观测器的滑模增益且kg≤-|e2/J|,
Figure FDA0002665884170000017
为负载转矩观测误差,
Figure FDA0002665884170000018
为负载转矩观测值,TL为负载转矩;
负载转矩观测器根据负载转矩给定值和负载转矩观测值的变化对反馈增益进行调整的方法是:
步骤1、计算
Figure FDA0002665884170000019
计算
Figure FDA00026658841700000110
步骤2、判断ΔT是否大于ε2;当ΔT大于ε2时,取反馈增益g等于gmin并进入步骤5;当ΔT小于等于ε2时,进入步骤3;
步骤3、判断ΔT是否小于ε1;当ΔT小于ε1时,取反馈增益g等于gmax并进入步骤5;当ΔT大于等于ε1时,进入步骤4;
步骤4、反馈增益g按照
Figure FDA0002665884170000021
进行计算;
步骤5、负载转矩观测器对负载转矩进行TL观测,得到负载转矩观测值
Figure FDA0002665884170000022
滑模速度控制器进行控制运算得到负载转矩给定值
Figure FDA0002665884170000023
其中,ε1为转矩变化低限比较阈值,ε2为转矩变化高限比较阈值,且0<ε1<ε2;gmax为反馈增益高值,gmin为反馈增益低值,且gmin<gmax<0。
2.如权利要求1所述的矿用牵引永磁同步电机驱动控制方法,其特征在于,转矩电流补偿分量i″q
Figure FDA0002665884170000024
q轴转矩电流给定值
Figure FDA0002665884170000025
Figure FDA0002665884170000026
3.如权利要求2所述的矿用牵引永磁同步电机驱动控制方法,其特征在于,选取gmin、gmax、ε1、ε2值的方法是:
步骤⑴,负载转矩观测器和滑模速度控制器均处于稳态且保持给定转子角速度不变和负载转矩不变;
步骤⑵,令反馈增益g从一个较大值开始逐渐减小,当负载转矩观测的稳态误差达到负载转矩观测稳态误差限值时,确定此时的反馈增益g值为gmax
步骤⑶,保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmax,连续进行n次ΔT值的测量,并将此时n次测量中最大m个ΔT值之和的平均值作为转矩变化低限比较阈值ε1
步骤⑷,微调改变反馈增益g,负载转矩观测器和滑模速度控制器均处于稳态时,保持给定转子角速度不变和令负载转矩突变,在保证负载转矩观测器输出观测值的转矩观测跟踪超调量处于转矩观测跟踪超调限值之内的前提下,测量负载转矩观测器的跟踪调节时间;
步骤⑸,重复步骤⑷,选择跟踪调节时间最短的反馈增益g值为gmin值;
步骤⑹,再次保持给定转子角速度不变和负载转矩不变且令反馈增益g等于gmin,连续进行n次ΔT值的测量,并将此时n次测量中最大m个ΔT值之和的平均值作为转矩变化高限比较阈值ε2
4.如权利要求1所述的矿用牵引永磁同步电机驱动控制方法,其特征在于,滑模增益kg按照
Figure FDA0002665884170000031
进行选择;其中,α≥1。
5.如权利要求4所述的矿用牵引永磁同步电机驱动控制方法,其特征在于,1≤α≤5。
6.如权利要求1所述的矿用牵引永磁同步电机驱动控制方法,其特征在于,检测永磁同步电机的转子位置θ、和三相电流ia、ib和ic;依据三相电流ia、ib和ic对永磁同步电机进行Clark变换得到在α-β轴坐标系下的电流iα、电流iβ,依据电流iα、电流iβ和转子位置θ进行Park变换,得到在d-q轴坐标系下的电流id、电流iq
CN202010918598.8A 2020-09-04 2020-09-04 矿用牵引永磁同步电机驱动控制方法 Active CN112039389B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202010918598.8A CN112039389B (zh) 2020-09-04 2020-09-04 矿用牵引永磁同步电机驱动控制方法
CN202210291754.1A CN114844405A (zh) 2020-09-04 2020-09-04 一种永磁同步电机积分滑模控制系统
CN202210291698.1A CN114865968A (zh) 2020-09-04 2020-09-04 一种永磁同步电机负载转矩观测器

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010918598.8A CN112039389B (zh) 2020-09-04 2020-09-04 矿用牵引永磁同步电机驱动控制方法

Related Child Applications (2)

Application Number Title Priority Date Filing Date
CN202210291754.1A Division CN114844405A (zh) 2020-09-04 2020-09-04 一种永磁同步电机积分滑模控制系统
CN202210291698.1A Division CN114865968A (zh) 2020-09-04 2020-09-04 一种永磁同步电机负载转矩观测器

Publications (2)

Publication Number Publication Date
CN112039389A CN112039389A (zh) 2020-12-04
CN112039389B true CN112039389B (zh) 2022-06-07

Family

ID=73591377

Family Applications (3)

Application Number Title Priority Date Filing Date
CN202210291698.1A Pending CN114865968A (zh) 2020-09-04 2020-09-04 一种永磁同步电机负载转矩观测器
CN202210291754.1A Pending CN114844405A (zh) 2020-09-04 2020-09-04 一种永磁同步电机积分滑模控制系统
CN202010918598.8A Active CN112039389B (zh) 2020-09-04 2020-09-04 矿用牵引永磁同步电机驱动控制方法

Family Applications Before (2)

Application Number Title Priority Date Filing Date
CN202210291698.1A Pending CN114865968A (zh) 2020-09-04 2020-09-04 一种永磁同步电机负载转矩观测器
CN202210291754.1A Pending CN114844405A (zh) 2020-09-04 2020-09-04 一种永磁同步电机积分滑模控制系统

Country Status (1)

Country Link
CN (3) CN114865968A (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112953317A (zh) * 2021-03-13 2021-06-11 无锡信捷电气股份有限公司 一种基于观测器的负载扰动快速抑制方法

Citations (4)

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

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107070342A (zh) * 2017-02-20 2017-08-18 哈尔滨理工大学 一种带负载状态观测器的永磁同步电机控制系统
CN110492817B (zh) * 2019-08-05 2021-08-03 北方工业大学 一种永磁同步电机的直接速度预测控制方法及设备
CN111431460B (zh) * 2020-04-25 2023-10-31 南通大学 一种永磁同步电机无传感器模型预测磁链控制方法

Patent Citations (4)

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

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
Cascaded Observers for the Rotor Position, Speed and Load Torque of the PMSM using an Extended Dynamic Model;Mihai COMANESCU;《XVI-th International Conference on Electrical Machines, Drives and Power Systems ELMA 2019, 6-8 June 2019, Varna, Bulgaria》;20190630;全文 *
Low Speed Control of Hybrid Excitation Synchronous Machines Based on a Novel Load Torque Sliding Mode Observer;Bingchu Liu等;《Proceedings of the 39th Chinese Control Conference》;20200729;全文 *
基于负载转矩滑模观测的永磁同步电机滑模控制;张晓光 等;《中国电机工程学报》;20120125;第32卷(第3期);全文 *
永磁同步电机负载转矩观测器;郑泽东 等;《电工技术学报》;20100228;第25卷(第8期);全文 *

Also Published As

Publication number Publication date
CN114865968A (zh) 2022-08-05
CN112039389A (zh) 2020-12-04
CN114844405A (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
CN110429881B (zh) 一种永磁同步电机的自抗扰控制方法
CN112039390B (zh) 基于负载转矩观测的永磁同步电机滑模控制方法
CN110138297B (zh) 一种永磁同步直线电机速度和电流双闭环控制系统和控制方法
CN110572091B (zh) 一种优化永磁同步电机无传感器控制方法
CN112422004B (zh) 一种永磁同步电机弱磁控制模式下的扰动抑制方法
CN113206623B (zh) 基于快速积分终端滑模及干扰估计的永磁同步电机有限时间调速控制方法
CN110138298B (zh) 一种永磁同步电机滑模控制方法
CN110165953B (zh) 一种基于趋近律的pmsm调速控制方法
CN112290843B (zh) 一种变指数幂次趋近律及其pmsm控制应用
CN112054727B (zh) 电梯永磁同步电机驱动控制系统
Rong et al. A new PMSM speed modulation system with sliding mode based on active-disturbance-rejection control
CN112953335A (zh) 一种永磁同步电机有限时间自适应复合控制方法和系统
CN112039394A (zh) 一种基于模糊自抗扰的pmsm伺服控制系统
CN112039389B (zh) 矿用牵引永磁同步电机驱动控制方法
CN112039388B (zh) 工业机器人用永磁同步电机控制方法
CN112039391B (zh) 高压断路器电机驱动控制系统
CN112072974B (zh) 风机电动变桨距电机驱动控制系统
CN110943666A (zh) 一种永磁同步电机复合电流的约束控制系统及其构建方法
CN115378325A (zh) 基于动态权重因子的smpmsm驱动系统直接速度复合控制方法
CN113328672B (zh) 无位置传感器永磁电机无差拍电流预测的控制方法及系统
CN115133828A (zh) 一种永磁同步电机控制方法及系统
Han et al. Research on PMSM sensor-less system based on ADRC strategy
CN115700991A (zh) 一种永磁同步电机的高精度电流补偿系统及补偿方法
CN114244214A (zh) 基于改进型滑模控制的永磁同步电动机位置控制方法
CN114629401A (zh) 一种基于模型参考自适应的异步电动机的无速度传感器矢量控制系统建立方法

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