CN106208865A - 基于负载观测器的多永磁同步电机虚拟总轴控制方法 - Google Patents
基于负载观测器的多永磁同步电机虚拟总轴控制方法 Download PDFInfo
- Publication number
- CN106208865A CN106208865A CN201610662965.6A CN201610662965A CN106208865A CN 106208865 A CN106208865 A CN 106208865A CN 201610662965 A CN201610662965 A CN 201610662965A CN 106208865 A CN106208865 A CN 106208865A
- Authority
- CN
- China
- Prior art keywords
- motor
- load torque
- synchronous motors
- permanent magnet
- torque
- 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.)
- Granted
Links
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P21/00—Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
- H02P21/13—Observer control, e.g. using Luenberger observers or Kalman filters
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/04—Arrangements for controlling or regulating the speed or torque of more than one motor
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02P—CONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
- H02P6/00—Arrangements for controlling synchronous motors or other dynamo-electric motors using electronic commutation dependent on the rotor position; Electronic commutators therefor
- H02P6/04—Arrangements for controlling or regulating the speed or torque of more than one motor
- H02P2006/045—Control of current
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Control Of Multiple Motors (AREA)
Abstract
本发明属于多电机速度协同控制领域,为协调各从属轴之间的同步运行,也更真实地反映各轴之间动力学关系,同时为进一步提高***的鲁棒性,将负载转矩观测值前馈至从属轴电机的电流环。本发明采取以下技术方案:基于负载观测器的多永磁同步电机虚拟总轴控制方法,***中包括n台永磁同步电机和1台虚拟电机,使用微处理器内部AD转换接口采集整流桥输出的直流电压和相电流接着,通过转子位置传感器检测出第i台永磁同步电机的转子角位置θi,同时计算出第i台永磁同步电机的转速ωi;然后生成α、β轴参考电压;最后利用电压空间矢量脉宽调制SVPWM方法输出逆变器的6脉冲驱动控制信号。本发明主要应用于多电机速度协同控制场合。
Description
技术领域
本发明涉及基于负载观测器的多电机虚拟总轴控制方法,属于多电机速度协同控制领域。具体讲,涉及基于负载观测器的多永磁同步电机虚拟总轴控制方法。
背景技术
随着现代工业技术水平的不断发展,多电机同步驱动控制技术已经被应用到各个领域,例如造纸机、印刷机、纺织机、多自由度数控机床,甚至大功率的高速铁路牵引***等。近些年来在造纸和印刷行业,卷筒纸传输单元和印刷机印刷单元的机械轴驱动方式正大量地被无轴传动方式所取代,即由多个独立电机替代一个主电机同时驱动各印刷机组、接纸装置、折页装置或裁切装置。采用独立电机驱动具有诸多优点,如它避免了机械轴驱动带来的机械振动,其次独立电机驱动的印刷机配置非常灵活,不同于机械轴连接需要添加或拆除机械部件,它只需简单连线操作即可。目前,常用的多轴传动同步控制策略有主从控制、交叉耦合控制、偏差耦合控制和电子虚拟总轴(Electronic virtual line-shafting,ELS)控制策略等。而已有的资料已表明ELS控制更适合应用在卷筒纸传输和印刷机印刷单元。
ELS可以实现电机速度或位置调节,它为从属轴提供速度或位置的参考信号。稳态时,各从属轴跟随总轴,能够达到很好的同步效果;但传统ELS控制***在起动、负载发生扰动、停机的过程中,因耦合力矩的反馈存在时间滞后会引起各从属轴之间失同步的现象,这是ELS控制策略的技术难题。同时考虑到多电机控制***是一个多变量、参数时变、速度和张力强耦合***,传统比例积分微分(Proportional integral derivative,PID)控制并不能满足高性能控制的要求。
发明内容
传统ELS控制***在电机起动、负载发生扰动、停机的过程中,因耦合力矩反馈存在时间滞后而引起轴之间失同步的现象。针对上述问题,本发明方法设计了负载观测器,将观测的负载转矩作为反馈力矩反馈回虚拟总轴,从而协调了各从属轴之间的同步运行,也更真实地反映了各轴之间动力学关系,同时为进一步提高***的鲁棒性,将负载转矩观测值前馈至从属轴电机的电流环。
为了实现上述目的,本发明采取以下技术方案:基于负载观测器的多永磁同步电机虚拟总轴控制方法,***中包括n台永磁同步电机和1台虚拟电机,使用微处理器内部AD转换接口采集整流桥输出的直流电压Udc和相电流iai、ibi、ici,iai、ibi、ici分别为第i台永磁同步电机定子电流;由耦合力矩Ti得到电流的参考值取值为零,Ti是第i台永磁同步电机耦合力矩,分别为第i台永磁同步电机电流环定子电流直轴分量与交轴分量的参考值,此时永磁同步电机能够产生的最大转矩;接着,通过转子位置传感器检测出第i台永磁同步电机的转子角位置θi,同时计算出第i台永磁同步电机的转速ωi;然后通过电流传感器检测第i台永磁同步电机定子电流iai、ibi、ici,并经三相静止坐标变换成二相静止坐标Clarke变换和二相静止坐标变换成二相旋转坐标Park变换,得到第i台永磁同步电机电流环定子电流直轴分量与交轴分量的检测值idi、iqi;然后减去iqi得到q轴电流误差,减去idi得到d轴电流误差,电流误差量分别经PI调节器输出d、q轴上的电压值再经过反Park变换后生成电压值分别为第i台永磁同步电机经过PI调节器得到的d、q轴上的参考电压;分别为第i台永磁同步电机经过反Park变换后生成α、β轴参考电压;最后利用电压空间矢量脉宽调制SVPWM方法输出逆变器的6脉冲驱动控制信号;其中,采用改进的电子虚拟总轴ELS(Electronic virtual line-shafting)控制策略:利用负载观测器将观测的负载转矩反馈回虚拟总轴,虚拟总轴输出的参考信号即虚拟总轴输出的角速度ω、虚拟总轴输出的角位置θ发生变化,在ω、θ作用下,受扰电机和其它从属轴电机的转速差减小,使不同永磁同步电机转矩更加平衡;同时,将观测负载转矩前馈至从属轴电机的PI调节器,此时的力矩平衡关系式为
式中,T为虚拟总轴的驱动力矩,Js为虚拟总轴的转动惯量。
负载观测器的状态方程为:负载观测器和滑模控制器设计所用滑模趋近率如下
式中,s为自变量;–ks是指数趋近项,k为正常数;–εsgn(s)是等速趋近项,ε为正常数,sgn()是符号函数。
将第i台电机的角位置及转速作为***的状态变量有
式中,xi1、xi2分别代表第i台永磁同步电机的角位置和转速状态变量,分别为xi1、xi2的导数;是θi的导数;是ωi的导数,npi为第i台从属轴电机极对数,Ji为第i台电机转动惯量,ψfi为第i台电机转子磁链,TLi为第i台从属轴电机负载转矩。
定义变量如下
将的状态方程进一步整理得
在上式基础上,以第i台电机的角位置和转速为观测对象,建立滑模观测器为
式中,分别是第i台电机的角位置观测值的导数和转速观测值的导数;Wi1、Wi2是指数趋近律函数,即为
式中,ki1、ki2、εi1、εi2为正常数,用电机角位置和转速的观测值与实际值做差得观测误差为
式中,ei1和ei2分别是角位置和转速的观测误差,则观测误差ei1和ei2的导数为
定义滑模面为
式中,si1、si2是定义的滑模面,ci1>0为滑模面参数。
负载观测器进一步的具体工作方式是:
根据李雅普诺夫稳定性原理,取李雅普诺夫稳定性函数Vi1为
当滑模观测器进入滑动模态时,满足即电机的角位置误差又有ci1>0可得即得
ei2=Wi1=εi1sgn(ei1)+ki1ei1
同理滑模观测器进入滑动模态时,满足因si2=ci1ei1+ei2得又有得即得
LiTLi=Wi2=εi2sgn(ei2)+ki2ei2
定义为负载观测器值,则由上式可得的值为
为进一步提高***的抗扰性能和快速性,将观测负载转矩前馈至从属轴电机的电流环,即
式中,电流分量i′qi是负载转矩前馈至从属轴电机的电流环的分量,即将观测到的负载转矩前馈至从属轴电机的电流环,同时考虑更进一步地减小抖振幅值,用饱和函数sat()替代符号函数sgn(),即
式中,s为自变量;Δ为“边界层”,饱和函数sat()的本质:在边界层外,采用切换控制;在边界层内,采用线性化反馈控制。
在对第i台永磁同步电机建立数学模型时,作如下假设:1)假设相绕组中感应电动势波形为正弦,转子永磁磁场在气隙空间分布为标准的正弦波;2)忽略定子铁心饱和,不计涡流和磁滞损耗;3)永磁材料的电导率为零;4)转子上无阻尼绕组;采用最大转矩控制的永磁同步电机PMSM转子磁场定向控制方式,电压方程如下
式中,Ldi,Lqi分别为第i台电机定子直、交轴电感;Rsi为第i台电机定子绕组电阻;ψfi为第i台电机转子磁链。
PMSM的运动方程为
PMSM的转矩方程为
式中,Tei为第i台电机电磁转矩;TLi为第i台电机负载转矩;npi为第i台电机极对数;ωri为第i台电机的机械转速,ωri=ωi/npi;Ji为第i台电机转动惯量。
本发明的有益效果如下:
本发明对传统的ELS控制策略进行了改进,通过加入负载观测器解决了传统ELS控制***在起动、负载发生扰动、停机的过程中,因耦合力矩Ti(i=1,2,…,n;i代表电机台数)的反馈存在时间滞后而引起从属轴之间失同步的问题。同时考虑到多电机控制***是一个多变量、参数时变、速度和张力强耦合***,传统PID控制并不能满足高性能控制的要求,而SMC能够克服***的不确定性,对干扰和未建模***具有很强的鲁棒性等优点,所以本发明将指数趋近律算法应用在负载观测器设计部分,以进一步提高***的鲁棒性。
附图说明
图1为单台从属轴电机控制结构图。
图2为传统ELS控制结构图。
图3为耦合力矩结构图。
图4为改进型ELS控制结构图。
图5为负载转矩观测器原理图。
具体实施方式
本发明将观测的负载转矩作为反馈力矩反馈回虚拟总轴,从而协调了各从属轴之间的同步运行,也更真实地反映了各轴之间动力学关系,同时为进一步提高***的鲁棒性,将负载转矩观测值前馈至从属轴电机的电流环。
滑模变结构控制(sliding-mode variable structure control,SMC)能够克服***的不确定性,对干扰和未建模***具有很强的鲁棒性等优点,所以本发明将指数趋近律算法应用在负载观测器设计部分。
本发明采用的技术方案基于负载观测器的多永磁同步电机虚拟总轴控制,包括下列几个方面:
(1)***中包括n台永磁同步电机(permanent magnet synchronous motor,PMSM),在对第i(i=1,2,…,n)台电机建立数学模型时,作如下假设:1)假设相绕组中感应电动势波形为正弦,转子永磁磁场在气隙空间分布为标准的正弦波;2)忽略定子铁心饱和,不计涡流和磁滞损耗;3)永磁材料的电导率为零;4)转子上无阻尼绕组;采用最大转矩控制的PMSM转子磁场定向控制方式,电压方程如下
式中,udi、uqi分别为第i台电机定子端电压的直轴与交轴分量;idi,iqi分别为第i台电机定子电流的直轴分量与交轴分量;Ldi,Lqi分别为第i台电机定子直、交轴电感;Rsi为第i台电机定子绕组电阻;ωi为第i台电机转速;ψfi为第i台电机转子磁链。
PMSM的运动方程为
式中,Tei为第i台电机电磁转矩;TLi为第i台电机负载转矩;npi为第i台电机极对数;ωri为第i台电机的机械转速,ωri=ωi/npi;Ji为第i台电机转动惯量。
PMSM的转矩方程为
(2)改进的ELS控制策略
此时的力矩平衡关系式为
式中,是第i台电机的观测负载转矩;T为虚拟总轴的驱动力矩;Js为虚拟总轴的转动惯量;θ为虚拟总轴输出的角位置信号。
(3)设计负载观测器的状态方程
负载观测器和滑模控制器设计所用滑模趋近率如下
式中,s为自变量;–ks是指数趋近项,k为正常数;–εsgn(s)是等速趋近项,ε为正常数,sgn()是符号函数。
将第i台电机的角位置及转速作为***的状态变量有
式中,xi1、xi2分别代表第i台电机的角位置和转速状态变量,分别为xi1、xi2的导数;是第i台电机角位置的导数;是第i台电机转速的导数。
定义变量如下
将的状态方程进一步整理得
在上式基础上,以第i台电机的角位置和转速为观测对象,建立滑模观测器为
式中,分别是第i台电机的角位置观测值的导数和转速观测值的导数;Wi1、Wi2是指数趋近律函数,即为
式中,ki1、ki2、εi1、εi2为正常数。用电机角位置和转速的观测值与实际值做差得观测误差为
式中,ei1和ei2分别是角位置和转速的观测误差,则观测误差ei1和ei2的导数为
定义滑模面为
式中,si1、si2是定义的滑模面,ci1>0为滑模面参数。
(4)负载观测器的分析
根据李雅普诺夫稳定性原理,取李雅普诺夫稳定性函数Vi1为
当滑模观测器进入滑动模态时,满足即电机的角位置误差又有ci1>0可得即得
ei2=Wi1=εi1sgn(ei1)+ki1ei1
同理滑模观测器进入滑动模态时,满足因si2=ci1ei1+ei2得又有得即得
LiTLi=Wi2=εi2sgn(ei2)+ki2ei2
定义为负载观测器值,则由上式可得的值为
为进一步提高***的抗扰性能和快速性,将观测负载转矩前馈至从属轴电机的电流环,即
式中,电流分量i′qi是负载转矩前馈至从属轴电机的电流环的分量,即将观测到的负载转矩前馈至从属轴电机的电流环,同时考虑更进一步地减小抖振幅值,设计了用饱和函数sat()替代符号函数sgn(),即
式中,s为自变量;Δ为“边界层”,饱和函数sat()的本质:在边界层外,采用切换控制;在边界层内,采用线性化反馈控制。
下面根据附图和具体实例对本发明做进一步详述,所提出的控制方法具体包括如下步骤:
(1)***中包括3台永磁同步电机(permanent magnet synchronous motor,PMSM),图1为第i(i=1,2,3)台从属轴电机控制结构图,该结构图包含了从属轴电机电流环、电压生成算法及永磁同步电机等,图中ω为虚拟总轴输出的转速;图中ωi为第i台电机的转速;θ为虚拟总轴输出的角位置信号;θi是第i台电机的角位置;Ki为耦合机构的刚度系数;Bi为耦合机构的内部阻尼系数;Ti是耦合力矩;Kt是电机转矩系数;分别为第i台电机电流环定子电流直轴分量与交轴分量的参考值;idi,iqi分别为第i台电机定子电流的直轴分量与交轴分量;i′qi是第i台电机观测负载转矩前馈至从属轴电机的电流环的分量;K′t为前馈系数; 分别为第i台电机经过电流环调节器得到的d、q轴上的参考电压;分别为第i台电机经过反Park变换后生成α、β轴参考电压;Udc是整流桥输出的直流电压;iai、ibi、ici分别为第i台电机定子电流;PI为比例积分控制器;SVPWM是电压空间矢量脉宽调制技术的英文缩写。
在对第i(i=1,2,3)台电机建立数学模型时,作如下假设:1)假设相绕组中感应电动势波形为正弦,转子永磁磁场在气隙空间分布为标准的正弦波;2)忽略定子铁心饱和,不计涡流和磁滞损耗;3)永磁材料的电导率为零;4)转子上无阻尼绕组。采用最大转矩控制的PMSM转子磁场定向控制方式,电压方程如下
式中,udi、uqi分别为第i台电机定子端电压的直轴与交轴分量;Ldi,Lqi分别为第i台电机定子直、交轴电感;Rsi为第i台电机定子绕组电阻;ψfi为第i台电机转子磁链。
PMSM的运动方程为
式中,Tei为第i台电机电磁转矩;TLi为第i台电机负载转矩;npi为第i台电机极对数;ωri为第i台电机的机械转速(ωri=ωi/npi);Ji为第i台电机转动惯量。
PMSM的转矩方程为
在本实施例中,微处理器选用TI公司的TMS320F28335,使用微处理器内部AD转换接口采集直流电压Udc和相电流iai、ibi、ici。刚度系数Ki(i=1,2,3)取值3.00;阻尼系数Bi(i=1,2,3)取值0.03;电机转矩系数Kt取值0.98;第i台电机的转动惯量取值27.2Kg.cm2。首先,由耦合力矩Ti得到电流的参考值取值为零(此时电机能够产生的最大转矩);接着,通过转子位置传感器检测出转子角位置θi,同时计算出转速ωi;然后通过电流传感器检测定子电流iai、ibi、ici,并经Clarke变换(三相静止坐标变换成二相静止坐标)和Park变换(二相静止坐标变换成二相旋转坐标),得到检测值idi、iqi;然后减去iqi得到q轴电流误差,减去idi得到d轴电流误差,电流误差量分别经PI调节器输出d、q轴上的电压值再经过反Park变换后生成电压值最后利用SVPWM方法输出逆变器的6脉冲驱动控制信号。
三相电流iai、ibi、ici经过的Clarke变换(C3s/2s)矩阵为
两相电流iαi、iβi经过的Park变换(C2s/2r)矩阵为
经过的反Park变换(C2r/2s)矩阵为
PI控制器的数学描述为
u(t)=kpe(t)+ki∫e(t)dt (7)
式中,e(t)为输入量,u(t)为输出量,kp是比例系数,ki是积分常数。比例环节的作用是成比例地控制***的偏差e(t),偏差一旦产生,控制器立即产生控制作用,以减少偏差。积分环节的作用主要用于消除静差,提高***的跟踪能力。在本实施例中三台从属轴电机的电流环PI控制器参数为kp取值0.6,ki取值0.1。
针对具体情况举例说明,例如第i(i=1,2,3)台电机的负载转矩TLi突然增加,由式(2)可知电磁转矩不变,电机的转速ωi会减小,带动角位置θi会减小,转速、角位置通过负反馈引起耦合力矩Ti增加,Ti增加引起电机电流参考值增加;同时定子电流iai、ibi、ici瞬态情况下不变,而经Clarke变换和Park变换,得到电流反馈值idi、iqi也不变;从而减去iqi得到q轴电流误差会增大,经过电流PI控制器后增加,再经过反Park变换和空间矢量调制后,该电压会加在电机定子上;由式(1)可知此时电机定子电流iqi增大,由式(3)可知此时电机电磁转矩会增大,最后由式(2)可知,此时电机的转速ωi、角位置θi慢慢增大,最终电机输出的电磁转矩Tei与电机的负载转矩TLi平衡,电机的转速、角位置经历先减小后增加的过程,最终会跟随虚拟电机输出信号,保持稳定。
为进一步提高***的抗扰性能和快速性,将负载观测转矩前馈至从属轴电机的电流环,将iqi和θi输入到负载观测器模块,估算出负载转矩的值。即
式中,电流分量i′qi是前馈至从属轴电机的电流环的分量,K′t为前馈系数。在发生负载扰动时,将估算的i′qi和iqi运算之后作为第i(i=1,2,3)台电机的电流PI控制器输入量,从而缩短了电机的电流环对负载扰动的反应时间,进一步提高了***应对负载时变的能力。
(2)传统虚拟总轴(Electronic virtual line-shafting,ELS)控制策略的分析
图2是传统ELS控制结构图。图中ωref是***给定的参考转速;ω、θ分别为虚拟总轴输出的转速信号和角位置信号,分别作为各从属轴电机的转速参考信号和角位置参考信号;ωi(i=1,2,3)、θi(i=1,2,3)分别为各台电机的转速信号和角位置信号;Ki(i=1,2,3)为耦合机构的刚度系数;Bi(i=1,2,3)为耦合机构的内部阻尼系数;Ti(i=1,2,3)为各台电机反馈的耦合力矩;是积分符号。
ELS可以实现电机速度或位置的调节,它为从属轴电机提供速度ωi(i=1,2,3)或位置θi(i=1,2,3)的参考信号。稳态时,各从属轴电机跟随虚拟电机,能够达到很好的同步效果。当某一电机或多台电机受到干扰而偏离参考值时,通过耦合力矩Ti(i=1,2,3)的反馈,使虚拟电机感受这种变动,从而迫使从属轴电机跟随这种变动,实现了在瞬态时各从属轴电机之间的同步。
耦合力矩Ti的模型为
Ti=Ki∫(ω-ωi)dt+Bi(ω-ωi) (9)
耦合力矩Ti结构框图如图3所示。
虚拟总轴的力矩平衡关系式为
式中,T为虚拟总轴的驱动力矩;Js为虚拟总轴的转动惯量;θ为虚拟总轴输出的角位置信号。
在实际生产过程中,耦合机构的刚度系数Ki和内部阻尼系数Bi变化范围十分广泛,Ki和Bi取值具体可参考表1。而在ELS控制***中,Ki和Bi采用软件编程模拟机械轴之间的传动,省却了添加或拆除机械零部件,提高了***的灵活性。
表1参数选择
针对具体情况举例说明,例如第i(i=1,2,3)台电机的负载转矩TLi突然增加,则电机的转速ωi、角位置θi会减小,由式(9)知ωi、θi减小引起耦合力矩Ti变大。耦合力矩Ti反馈给虚拟电机,由式(10)可知在虚拟总轴驱动力矩T不变的情况下,Ti增大则虚拟电机输出的参考信号ω、θ减小。由步骤(1)可知第i(i=1,2,3)台电机受扰之后,其转速ωi、角位置θi信号先减小;而因参考信号ω、θ减小了,其他从属轴电机的转速ωi、角位置θi也跟着减小;所以最终受扰电机和其它从属轴电机的转速差减小,从而实现了在瞬态时各从属轴电机之间转速保持同步,转矩保持均衡。
传统ELS控制***在起动、负载发生扰动、停机的过程中,因耦合力矩Ti的反馈存在时间滞后而引起从属轴之间同步误差大,本发明针对这一问题对传统ELS控制策略进行了改进。
(3)改进的ELS控制策略
图4是改进的ELS控制结构图,iqi(i=1,2,3)为第i台电机定子电流的直轴分量;为第i台电机的观测负载转矩。当某台电机或多台电机的负载产生扰动时,负载观测器将观测的负载转矩反馈回虚拟总轴,虚拟总轴输出的参考信号ω、θ发生变化,在参考信号ω、θ作用下,受扰电机和其它从属轴电机的转速差减小,使不同电机转矩更加平衡;同时为更进一步提高***的鲁棒性,改进的ELS控制策略将观测的负载转矩前馈至从属轴电机的电流环;这两方面的作用能使***快速同步,此时的力矩平衡关系式为
在本实施例中,虚拟电机的转动惯量取值50.0Kg.cm2,虚拟电机的PI控制器参数为kp取值16.0,ki取值9.0,各从属轴电机参数取值在步骤(1)中已给出。
针对具体情况举例说明,例如第i(i=1,2,3)台电机的负载转矩TLi突然增加,一方面因观测的负载转矩前馈至从属轴电机的电流环,由步骤(1)可知电机的转速ωi、角位置θi开始减小,但前馈量使电压参考值更快地增大以平衡负载转矩;另一方面将观测的负载转矩取代耦合力矩Ti(i=1,2,3)反馈回虚拟总轴,由式(11)可知在虚拟总轴驱动力矩T不变的情况下,增大则虚拟电机输出的参考信号ω、θ会减小,而参考信号ω、θ减小,其他从属轴电机的转速ωi、角位置θi也跟着减小;所以最终受扰电机和其它从属轴电机的转速差减小,从而实现了在瞬态时各从属轴电机之间的快速同步。
(4)设计负载转矩观测器
将第i台电机的角位置及转速作为***的状态变量有
式中,xi1、xi2分别代表第i台电机的角位置和转速状态变量,分别为xi1、xi2的导数;是第i台电机角位置的导数;是第i台电机转速的导数。
定义变量如下
将的状态方程进一步整理得
负载观测器和滑模控制器设计所用滑模趋近率如下
式中,s为自变量;–ks是指数趋近项,k为正常数;–εsgn(s)是等速趋近项,ε为正常数,sgn()是符号函数。
在上式基础上,以第i台电机的角位置和转速为观测对象,建立滑模观测器为
式中,分别是第i台电机的角位置观测值的导数和转速观测值的导数;Wi1、Wi2是指数趋近律函数,即为
式中,ki1、ki2、εi1、εi2为正常数。用电机角位置和转速的观测值与实际值做差得观测误差为
式中,ei1和ei2分别是角位置和转速的观测误差,则观测误差ei1和ei2的导数为
定义滑模面为
式中,si1、si2是定义的滑模面,ci1>0为滑模面参数。
(5)负载观测器的分析
根据李雅普诺夫稳定性原理,取李雅普诺夫稳定性函数Vi1为
当滑模观测器进入滑动模态时,满足即电机的角位置误差又有ci1>0可得即得
ei2=Wi1=εi1sgn(ei1)+ki1ei1 (22)
同理当滑模观测器进入滑动模态时,满足因si2=ci1ei1+ei2得又有得即得
LiTLi=Wi2=εi2sgn(ei2)+ki2ei2 (23)
定义为负载观测器值,由式(23)得的值为
由上文对滑模观测器的分析可知,在εi1、εi2、ki1和ki2的选取范围内合理地选择参数,可保证滑模观测器稳定且可以准确地观测负载转矩。
本实施例中指数趋近律函数Wi1、Wi2的系数取值分别为εi1=2,ki1=5000,εi2=0.4,ki2=50。本实施例使用第i台电机的电流iqi和转子角θi估算出电机的负载转矩负载转矩观测原理如图5所示,式(16)是设计负载转矩观测器的基础,再由式(22)得到转速的观测误差ei2,最后由式(24)得到负载转矩的观测值。
考虑更进一步地减小抖振幅值,设计了用饱和函数sat()替代符号函数sgn(),即
式中,s为自变量;Δ为“边界层”。饱和函数sat(s)的本质:在边界层外,采用切换控制;在边界层内,采用线性化反馈控制。
综上所述,本发明对传统的ELS控制策略进行了改进,通过加入负载观测器解决了传统ELS控制***在起动、负载发生扰动、停机的过程中,因耦合力矩Ti的反馈存在时间滞后而引起从属轴之间失同步的问题。考虑到多电机控制***是一个多变量、参数时变、速度和张力强耦合***,传统PID控制并不能满足高性能控制的要求,而滑模变结构控制能够克服***的不确定性,对干扰和未建模***具有很强的鲁棒性等优点,所以本发明将指数趋近律算法应用在负载观测器设计部分,以进一步提高***的鲁棒性;同时为了克服负载时变对控制性能的影响,改进的ELS控制策略将观测的负载转矩前馈至从属轴电机的电流环。
Claims (4)
1.一种基于负载观测器的多永磁同步电机虚拟总轴控制方法,其特征是,***中包括n台永磁同步电机和1台虚拟电机,使用微处理器内部AD转换接口采集整流桥输出的直流电压Udc和相电流iai、ibi、ici,iai、ibi、ici分别为第i台永磁同步电机定子电流;由耦合力矩Ti得到电流的参考值取值为零,Ti是第i台永磁同步电机耦合力矩,分别为第i台永磁同步电机电流环定子电流直轴分量与交轴分量的参考值,此时永磁同步电机能够产生的最大转矩;接着,通过转子位置传感器检测出第i台永磁同步电机的转子角位置θi,同时计算出第i台永磁同步电机的转速ωi;然后通过电流传感器检测第i台永磁同步电机定子电流iai、ibi、ici,并经三相静止坐标变换成二相静止坐标Clarke变换和二相静止坐标变换成二相旋转坐标Park变换,得到第i台永磁同步电机电流环定子电流直轴分量与交轴分量的检测值idi、iqi;然后减去iqi得到q轴电流误差,减去idi得到d轴电流误差,电流误差量分别经PI调节器输出d、q轴上的电压值再经过反Park变换后生成电压值分别为第i台永磁同步电机经过PI调节器得到的d、q轴上的参考电压;分别为第i台永磁同步电机经过反Park变换后生成α、β轴参考电压;最后利用电压空间矢量脉宽调制SVPWM方法输出逆变器的6脉冲驱动控制信号;其中,采用改进的电子虚拟总轴ELS(Electronic virtual line-shafting)控制策略:利用负载观测器将观测的负载转矩反馈回虚拟总轴,虚拟总轴输出的参考信号即虚拟总轴输出的角速度ω、虚拟总轴输出的角位置θ发生变化,在ω、θ作用下,受扰电机和其它从属轴电机的转速差减小,使不同永磁同步电机转矩更加平衡;同时,将观测负载转矩前馈至从属轴电机的PI调节器,此时的力矩平衡关系式为:
式中,T为虚拟总轴的驱动力矩,Js为虚拟总轴的转动惯量。
2.如权利要求1所述的基于负载观测器的多永磁同步电机虚拟总轴控制方法,其特征是,负载观测器的状态方程为是,负载观测器和滑模控制器设计所用滑模趋近率如下:
式中,s为自变量;–ks是指数趋近项,k为正常数;–εsgn(s)是等速趋近项,ε为正常数,sgn()是符号函数。
将第i台电机的角位置及转速作为***的状态变量有
式中,xi1、xi2分别代表第i台永磁同步电机的角位置和转速状态变量,分别为xi1、xi2的导数;是θi的导数;是ωi的导数,npi为第i台从属轴电机极对数,Ji为第i台电机转动惯量,ψfi为第i台电机转子磁链,TLi为第i台从属轴电机负载转矩。
定义变量如下
将的状态方程进一步整理得
在上式基础上,以第i台电机的角位置和转速为观测对象,建立滑模观测器为
式中,分别是第i台电机的角位置观测值的导数和转速观测值的导数;Wi1、Wi2是指数趋近律函数,即为
式中,ki1、ki2、εi1、εi2为正常数,用电机角位置和转速的观测值与实际值做差得观测误差为
式中,ei1和ei2分别是角位置和转速的观测误差,则观测误差ei1和ei2的导数为
定义滑模面为
式中,si1、si2是定义的滑模面,ci1>0为滑模面参数。
3.如权利要求1所述的基于负载观测器的多永磁同步电机虚拟总轴控制方法,其特征是,负载观测器进一步的工作方式是,根据李雅普诺夫稳定性原理,取李雅普诺夫稳定性函数Vi1为:
当滑模观测器进入滑动模态时,满足即电机的角位置误差又有ci1>0可得即得
ei2=Wi1=εi1sgn(ei1)+ki1ei1
同理滑模观测器进入滑动模态时,满足因si2=ci1ei1+ei2得又有得即得
LiTLi=Wi2=εi2sgn(ei2)+ki2ei2
定义为负载观测器值,则由上式可得的值为
为进一步提高***的抗扰性能和快速性,将观测负载转矩前馈至从属轴电机的电流环,即
式中,电流分量iq′i是负载转矩前馈至从属轴电机的电流环的分量,即将观测到的负载转矩前馈至从属轴电机的电流环,同时考虑更进一步地减小抖振幅值,用饱和函数sat()替代符号函数sgn(),即
式中,s为自变量;Δ为“边界层”,饱和函数sat()的本质:在边界层外,采用切换控制;在边界层内,采用线性化反馈控制。
4.如权利要求1所述的基于负载观测器的多永磁同步电机虚拟总轴控制方法,其特征是,在对第i台永磁同步电机建立数学模型时,作如下假设:1)假设相绕组中感应电动势波形为正弦,转子永磁磁场在气隙空间分布为标准的正弦波;2)忽略定子铁心饱和,不计涡流和磁滞损耗;3)永磁材料的电导率为零;4)转子上无阻尼绕组;采用最大转矩控制的永磁同步电机PMSM转子磁场定向控制方式,电压方程如下:
式中,Ldi,Lqi分别为第i台电机定子直、交轴电感;Rsi为第i台电机定子绕组电阻;ψfi为第i台电机转子磁链。
PMSM的运动方程为
PMSM的转矩方程为
式中,Tei为第i台电机电磁转矩;TLi为第i台电机负载转矩;npi为第i台电机极对数;ωri为第i台电机的机械转速,ωri=ωi/npi;Ji为第i台电机转动惯量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610662965.6A CN106208865B (zh) | 2016-08-10 | 2016-08-10 | 基于负载观测器的多永磁同步电机虚拟总轴控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610662965.6A CN106208865B (zh) | 2016-08-10 | 2016-08-10 | 基于负载观测器的多永磁同步电机虚拟总轴控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106208865A true CN106208865A (zh) | 2016-12-07 |
CN106208865B CN106208865B (zh) | 2018-09-18 |
Family
ID=57515306
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610662965.6A Active CN106208865B (zh) | 2016-08-10 | 2016-08-10 | 基于负载观测器的多永磁同步电机虚拟总轴控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106208865B (zh) |
Cited By (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106992720A (zh) * | 2017-05-26 | 2017-07-28 | 西门子工厂自动化工程有限公司 | 基于位置同步的多轴耦合力矩平衡控制方法及装置 |
CN107070315A (zh) * | 2017-05-31 | 2017-08-18 | 湖南科技大学 | 一种多开关磁阻电机转速同步控制装置及控制方法 |
CN108173462A (zh) * | 2017-12-13 | 2018-06-15 | 天津津航计算技术研究所 | 一种两电机转矩平衡控制方法 |
CN108270379A (zh) * | 2018-02-02 | 2018-07-10 | 上海交通大学 | 一种多电机同步高精度滑模控制方法 |
CN108667374A (zh) * | 2018-03-30 | 2018-10-16 | 青岛大学 | 一种高刚度、免调试的变频调速***速度环的设计方法 |
CN109291803A (zh) * | 2018-08-21 | 2019-02-01 | 沈阳工业大学 | 基于四轮全驱电动汽车虚拟轮的稳定性控制方法 |
CN109474207A (zh) * | 2018-11-22 | 2019-03-15 | 天津大学 | 一种提高多轴运动***同步性能的控制方法 |
CN109510521A (zh) * | 2018-11-06 | 2019-03-22 | 湖南工业大学 | 一种基于虚拟总轴的多电机总量协同一致容错控制方法 |
CN109861374A (zh) * | 2019-01-31 | 2019-06-07 | 张欣 | 一种无需负载电流传感器的三相全桥不间断电源控制方法 |
CN110737193A (zh) * | 2019-09-19 | 2020-01-31 | 中建材创新科技研究院有限公司 | 一种虚拟轴控制多轴的同步算法控制*** |
CN111900906A (zh) * | 2020-08-28 | 2020-11-06 | 北京大豪工缝智控科技有限公司 | 主轴转速的控制方法、装置、缝纫机和可读存储介质 |
CN112015084A (zh) * | 2020-07-08 | 2020-12-01 | 湖南工业大学 | 一种基于虚拟总轴的总量协同一致抗饱和控制方法 |
CN112542965A (zh) * | 2020-11-26 | 2021-03-23 | 同济大学 | 一种多电机同步控制***及方法 |
CN108762088B (zh) * | 2018-06-20 | 2021-04-09 | 山东科技大学 | 一种迟滞非线性伺服电机***滑模控制方法 |
EP3809578A1 (en) * | 2019-10-15 | 2021-04-21 | ABB Schweiz AG | Method for starting a synchronous motor |
CN113271042A (zh) * | 2021-05-18 | 2021-08-17 | 湖南工业大学 | 一种多电机固定时间优化协同控制方法 |
CN114421818A (zh) * | 2022-01-26 | 2022-04-29 | 合肥倍豪海洋装备技术有限公司 | 一种船用全回转推进器回转***的负荷平衡控制方法 |
EP4296806A1 (en) | 2022-06-24 | 2023-12-27 | Arnold NextG GmbH | System to control and/or to synchronize the motion of two shafts |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080018296A1 (en) * | 2006-07-07 | 2008-01-24 | Sanyo Electric Co., Ltd. | Motor control device |
CN101976997A (zh) * | 2010-10-25 | 2011-02-16 | 四川省安普瑞自动化设备有限公司 | 带式运输机多电机同步控制*** |
CN102420553A (zh) * | 2011-12-09 | 2012-04-18 | 南京理工大学 | 基于改进型相邻交叉耦合的多电机比例同步控制算法 |
CN103647490A (zh) * | 2013-09-27 | 2014-03-19 | 天津大学 | 一种永磁电机的滑模控制策略 |
-
2016
- 2016-08-10 CN CN201610662965.6A patent/CN106208865B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20080018296A1 (en) * | 2006-07-07 | 2008-01-24 | Sanyo Electric Co., Ltd. | Motor control device |
CN101976997A (zh) * | 2010-10-25 | 2011-02-16 | 四川省安普瑞自动化设备有限公司 | 带式运输机多电机同步控制*** |
CN102420553A (zh) * | 2011-12-09 | 2012-04-18 | 南京理工大学 | 基于改进型相邻交叉耦合的多电机比例同步控制算法 |
CN103647490A (zh) * | 2013-09-27 | 2014-03-19 | 天津大学 | 一种永磁电机的滑模控制策略 |
Non-Patent Citations (1)
Title |
---|
肖园园: "基于虚拟总轴的同步控制技术及其应用研究", 《中国优秀硕士学位论文全文数据库》 * |
Cited By (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106992720A (zh) * | 2017-05-26 | 2017-07-28 | 西门子工厂自动化工程有限公司 | 基于位置同步的多轴耦合力矩平衡控制方法及装置 |
CN106992720B (zh) * | 2017-05-26 | 2019-12-17 | 西门子工厂自动化工程有限公司 | 基于位置同步的多轴耦合力矩平衡控制方法及装置 |
CN107070315A (zh) * | 2017-05-31 | 2017-08-18 | 湖南科技大学 | 一种多开关磁阻电机转速同步控制装置及控制方法 |
CN107070315B (zh) * | 2017-05-31 | 2023-05-12 | 湖南科技大学 | 一种多开关磁阻电机转速同步控制装置及控制方法 |
CN108173462B (zh) * | 2017-12-13 | 2020-06-26 | 天津津航计算技术研究所 | 一种两电机转矩平衡控制方法 |
CN108173462A (zh) * | 2017-12-13 | 2018-06-15 | 天津津航计算技术研究所 | 一种两电机转矩平衡控制方法 |
CN108270379A (zh) * | 2018-02-02 | 2018-07-10 | 上海交通大学 | 一种多电机同步高精度滑模控制方法 |
CN108270379B (zh) * | 2018-02-02 | 2020-06-19 | 上海交通大学 | 一种多电机同步高精度滑模控制方法 |
CN108667374A (zh) * | 2018-03-30 | 2018-10-16 | 青岛大学 | 一种高刚度、免调试的变频调速***速度环的设计方法 |
CN108667374B (zh) * | 2018-03-30 | 2021-02-19 | 青岛大学 | 一种高刚度、免调试的变频调速***速度环的设计方法 |
CN108762088B (zh) * | 2018-06-20 | 2021-04-09 | 山东科技大学 | 一种迟滞非线性伺服电机***滑模控制方法 |
CN109291803A (zh) * | 2018-08-21 | 2019-02-01 | 沈阳工业大学 | 基于四轮全驱电动汽车虚拟轮的稳定性控制方法 |
CN109510521B (zh) * | 2018-11-06 | 2022-01-14 | 湖南工业大学 | 一种基于虚拟总轴的多电机总量协同一致容错控制方法 |
CN109510521A (zh) * | 2018-11-06 | 2019-03-22 | 湖南工业大学 | 一种基于虚拟总轴的多电机总量协同一致容错控制方法 |
CN109474207A (zh) * | 2018-11-22 | 2019-03-15 | 天津大学 | 一种提高多轴运动***同步性能的控制方法 |
CN109861374A (zh) * | 2019-01-31 | 2019-06-07 | 张欣 | 一种无需负载电流传感器的三相全桥不间断电源控制方法 |
CN109861374B (zh) * | 2019-01-31 | 2022-09-02 | 张欣 | 一种无需负载电流传感器的三相全桥不间断电源控制方法 |
CN110737193B (zh) * | 2019-09-19 | 2021-11-23 | 中建材创新科技研究院有限公司 | 一种虚拟轴控制多轴的同步算法控制*** |
CN110737193A (zh) * | 2019-09-19 | 2020-01-31 | 中建材创新科技研究院有限公司 | 一种虚拟轴控制多轴的同步算法控制*** |
EP3809578A1 (en) * | 2019-10-15 | 2021-04-21 | ABB Schweiz AG | Method for starting a synchronous motor |
US11527973B2 (en) | 2019-10-15 | 2022-12-13 | Abb Schweiz Ag | Method for starting a synchronous motor |
CN112015084B (zh) * | 2020-07-08 | 2023-01-13 | 湖南工业大学 | 一种基于虚拟总轴的总量协同一致抗饱和控制方法 |
CN112015084A (zh) * | 2020-07-08 | 2020-12-01 | 湖南工业大学 | 一种基于虚拟总轴的总量协同一致抗饱和控制方法 |
CN111900906A (zh) * | 2020-08-28 | 2020-11-06 | 北京大豪工缝智控科技有限公司 | 主轴转速的控制方法、装置、缝纫机和可读存储介质 |
CN112542965A (zh) * | 2020-11-26 | 2021-03-23 | 同济大学 | 一种多电机同步控制***及方法 |
CN113271042A (zh) * | 2021-05-18 | 2021-08-17 | 湖南工业大学 | 一种多电机固定时间优化协同控制方法 |
CN113271042B (zh) * | 2021-05-18 | 2023-07-04 | 湖南工业大学 | 一种多电机固定时间优化协同控制方法 |
CN114421818A (zh) * | 2022-01-26 | 2022-04-29 | 合肥倍豪海洋装备技术有限公司 | 一种船用全回转推进器回转***的负荷平衡控制方法 |
EP4296806A1 (en) | 2022-06-24 | 2023-12-27 | Arnold NextG GmbH | System to control and/or to synchronize the motion of two shafts |
Also Published As
Publication number | Publication date |
---|---|
CN106208865B (zh) | 2018-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106208865B (zh) | 基于负载观测器的多永磁同步电机虚拟总轴控制方法 | |
CN106887976A (zh) | 考虑加速度的多永磁同步电机偏差耦合控制方法 | |
Kim et al. | Sensorless control of interior permanent-magnet machine drives with zero-phase lag position estimation | |
Guo et al. | Model reference adaptive control of five-phase IPM motors based on neural network | |
CN101247104B (zh) | 定子磁链定向的交流电机旋转速度和转子位置推测方法 | |
CN106330046B (zh) | 基于特定负载的五相容错永磁电机无位置传感器控制方法 | |
CN105048913B (zh) | 基于电流补偿的无轴承异步电机不平衡振动控制*** | |
CN107482982B (zh) | 一种基于铁损模型的异步电机矢量控制方法 | |
CN103051269A (zh) | 同步机控制装置 | |
CN105610369B (zh) | 一种基于滑模观测器的异步电机磁链观测方法 | |
CN101969289A (zh) | 改善电机在场削弱区域工作期间的扭矩线性的增益调节 | |
CN102197583A (zh) | 用于控制ac发电机的方法和控制设备 | |
CN107911057A (zh) | 一种用于飞轮储能***的转子位置鲁棒观测方法 | |
CN102386837A (zh) | 矢量控制装置以及电动机控制*** | |
CN105680754A (zh) | 一种永磁同步电机的直交轴电流矢量复合控制器 | |
CN103560735A (zh) | 一种电励磁同步电机的控制方法 | |
CN104333285A (zh) | 永磁同步电机准无传感器位置伺服控制装置及方法 | |
CN103701375B (zh) | 一种永磁同步电机启动***、方法及永磁同步电机 | |
CN103401493A (zh) | 永磁同步电主轴驱动控制***和方法 | |
CN103036496A (zh) | 自适应反推控制的永磁同步电机dtc***及其控制方法 | |
CN106788049A (zh) | 基于级联滑模观测器的无速度传感器转矩控制***及方法 | |
CN110350838A (zh) | 一种基于扩展卡尔曼滤波的无速度传感器bdfim直接转矩控制方法 | |
CN109728755A (zh) | 一种pmsm反演终端滑模控制方法 | |
Liu et al. | A MTPA control strategy for mono-inverter multi-PMSM system | |
Chi | Position-sensorless control of permanent magnet synchronous machines over wide speed range |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |