1. 引言
永磁直线同步电机(PMLSM)是一种通过动子直接驱动,将电能直接转换成直线形式的机械能而无需中间传动件的驱动装置。相对于旋转电机,PMLSM具有更大的驱动力、更高的定位精度以及更快的响应速度,在精加工领域备受关注 [1] [2] 。在单轴伺服系统中,直线电机与负载直接相连,没有中间传动部件,其易受到参数摄动、外界扰动等不利因素的影响 [3] ,需要对这些不利因素进行抑制,以提高系统的抗干扰能力。由于单轴直线电机受到上述因素的影响,并且传统直接力控制中的PID控制器存在对参数变化的适应性差以及无法兼顾响应速度与超调量等问题,故而难以实现系统高性能跟踪控制。因此,提高直线电机跟踪精度的控制方法受到越来越多的关注,其中,由韩京清教授提出的自抗扰控制技术 [4] ,由PID演变而来,通过对系统扰动进行观测,并补偿给控制输入,使系统具有强鲁棒性 [5] [6] 。
目前,自抗扰控制技术已被成功地应用在PMLSM控制、工业机器人以及无人机等领域。文献 [7] 针对PMLSM控制系统易受到负载扰动、推力波动等因素的影响,采用ADRC对控制系统中的总扰动进行观测,并补偿给系统控制输入,从而提高系统控制精度;文献 [8] 通过对ADRC的扩张状态观测器进行改进,以提高对高频转矩脉动的观测精度,从而有效抑制转矩脉动;文献 [9] 对ADRC的非线性函数进行改进,并在非线性函数的最外侧采用常数的方式,从而加快减小对系统的增益,提高了对PMLSM的控制精度。自抗扰控制技术作为一种针对控制系统的控制策略,是通过集成反馈回路和基于计算的估计器来有效地管理和抵消系统内外部扰动,与其他控制器相比,自抗扰控制器能够独立于控制对象,对提高控制系统的鲁棒性具有显著的作用。
为提高控制系统跟踪精度和抗干扰能力,本文提出一种基于改善型自抗扰控制的永磁直线同步电机单轴伺服系统直接力控制。首先搭建PMLSM数学模型,然后设计改进型自抗扰控制,针对改进后的自抗扰控制参数整定问题,阐述了其参数整定的一般规律。最后,将I-ADRC作用于单轴伺服系统的DFC中,通过MATLAB/Simulink建立仿真模型,仿真分析验证改进型自抗扰控制器的直接力控制,具有更好的跟踪精度和控制性能。
2. 永磁同步直线电机数学模型
在对直线电机分析之前,为了便于模型建立,给出以下假设 [10] :
(1) 忽略初级的铁芯饱和、涡流以及磁滞损耗;
(2) 气隙均匀分布,且反电动势为正弦;
(3) 忽略次级的阻尼绕组;
(4) 忽略端部效应的影响;
(5) 三相绕组均匀分布且参数相同。
搭建永磁同步电机在d-q坐标系下的数学模型。
电压方程可表示为:
(1.1)
式中,
、
为初级在d、q轴上的电压分量;
、
为初级在d、q轴上的磁链分量;
、
为初级在d、q轴上的电流分量;
为绕组的等效电阻;
为直线电机初级水平移动速度转换成的旋转电机动子角速度。
角速度
可表示为:
(1.2)
式中,
为极距;
为极对数;v为直线电机移动速度。
在d、q轴上,初级磁链可表示为:
(1.3)
式中,
、
为初级在d、q轴上的电感,并且对于表贴式永磁直线同步电机有
;
为永磁体磁链 [11] 。
永磁直线同步电机的总功率可表示为:
(1.4)
根据式(1.4)可知,永磁直线同步电机的总功率由三部分组成。第一部分
为三相绕组线圈上的电阻热损耗功率;第二部分
为三相绕组上电流的变化导致的磁场储能的变化;第三部分
为三相绕组输出能量中转化为机械能输出的功率。
因此由永磁直线同步电机机械能输出功率可得的推力方程为:
(1.5)
根据式(1.3)和式(1.5)可进一步分析得到表贴式永磁直线同步电机的推力方程为:
(1.6)
式中,
为推力系数。
由于初级磁链在d-q坐标系下的表达式为:
(1.7)
式中,
为
与d轴的夹角。
故根据式(1.3)以及式(1.7),永磁直线同步电机的推力公式还可表示为:
(1.8)
机械运动方程为:
(1.9)
式中,
为直线电机动子质量;
为永磁直线同步电机运动过程中所受总扰动;
为粘滞系数;
为电机所受负载阻力。
根据以上推导公式,最终可得d-q坐标系下永磁直线同步电机的数学模型框图如图1所示。
![](//html.hanspub.org/file/26-2571576x41_hanspub.png?20240520094819129)
Figure 1. Mathematical model diagram of permanent magnet linear synchronous motor in d-q coordinate system
图1. d-q坐标系下永磁直线同步电机的数学模型框图
3. 改进型自抗扰控制器设计
自抗扰控制是基于传统PID的新型主动抗扰技术,在保留PID控制优点的基础上对其进行了优化 [12] 。首先采用跟踪微分器(Tracking Differentiator, TD),为输入的跳变信号提供过渡过程。其次,为了减小扰动,扩张状态观测器将对内部和外部扰动进行观测,并将扰动观测值反馈给非线性状态误差反馈控制律形成闭环控制 [13] ,其结构如图2所示。
自抗扰控制中的TD采用的是最速控制综合函数fhan,且ESO和NLSEF是通过非线性函数fal实现的 [14] ,因此该控制器可以根据误差变化对比例环节的增益系数进行实时调节,从而即实现了较小的超调量,又保证了较好的跟踪性能。并且该控制方法可不依赖于控制对象的具体模型,而是通过扰动的观测及补偿来达到自抗扰的目的。
对于传统ADRC,其ESO以及NLSEF的非线性函数fal在分段点处存在抖动以及切换不平滑等问题,从而降低了ADRC对系统的控制性能。因此,针对该问题,以下分别对ESO以及NLSEF进行改进优化。改进型ADRC(I-ADRC)结构如图3所示。
3.1. 改进型ESO设计
ESO中最关键的部分为非线性函数,由此可见fal函数对系统控制性能的重要性。fal函数的曲线特性由非线性因子
和分段区间长度
决定,
的取值与函数曲线的非线性度有关,
的取值决定了函数曲线的拐点位置,由于传统fal函数在分段点
处连续但不可导,当取较小的
值时,fal函数的斜率在原点附近会发生突变,这就使得系统状态在原点处产生抖动,从而影响系统的抗扰性能,这对单轴PMLSM的位置跟踪性能有着不利的影响。因此本文在传统fal函数的基础上给出了一种连续且光滑的非线性函数Ifal,以改进ESO的观测性能。
1) 当
时,Ifal函数的表达式保持原公式:
(1.10)
2) 当
时,Ifal函数的表达式采用线性函数和正弦函数拟合的形式代替单一线性函数,从而更好的满足“小误差大增益”的原则,其表达式为:
(1.11)
为了使得函数在
处连续且可导,需使得函数在该点左右两侧的函数值以及导数相同,即:
(1.12)
3) 当
时,为了满足“大误差小增益”的原则,需要加快减小系统的增益,从而提高系统的抗干扰能力。Ifal函数的表达式采用反比例函数的形式:
(1.13)
选用反比例函数形式而非直接使用常数,是因为在
取较小值时,若使用常数形式,同样会因Ifal函数存在拐点,使得系统状态在原点附近产生抖动。而反比例函数则具有更好的平滑性。
为了使得Ifal函数在
处连续且可导,需使得函数在该点左右两侧的函数值以及导数相同,即:
(1.14)
最终得到Ifal函数表达式如下:
(1.15)
接下来对fal函数与Ifal函数的性能进行对比分析,分别取
、
、
,并在MATLAB中进行图形绘制,得到fal函数与Ifal函数的对比分析图如图4所示。
(a) fal与Ifal的特性曲线 (b) fal与Ifal的误差增益曲线
Figure 4. Comparative analysis chart of fal and Ifal
图4. fal与Ifal的对比分析图
由图4(a)可知,fal函数在分段点处存在明显拐点,而Ifal函数则表现出更好的连续性与光滑性。由图4(b)可知,当
时,Ifal函数对比传统fal函数具有更小的输出增益;当输出误差e在零点附近时,Ifal函数的误差增益明显大于fal函数。
经过上述分析表明,与传统fal函数相比,Ifal函数更能体现“大误差小增益,小误差大增益”的原则,当误差较小时采用大增益,以增加系统的快速性;当误差较大时采用小增益,从而避免超调现象。并且Ifal函数克服了传统fal函数因在原点附近存在斜率突变而导致系统高频振荡的问题,因此采用Ifal函数的ESO更能提高系统的抗干扰能力。
3.2. 改进型NLSEF设计
改进型NLSEF(INLSEF)在传统NLSEF的基础上引入了误差的积分
以及积分项的权重
,并利用非线性函数Ifal将TD的输出与IESO的输出之间的误差、误差的微分、误差的积分进行非线性组合,从而克服稳态下高频振荡的问题,增强系统的抗干扰能力,进一步提高系统的跟随精度。其表达式如下:
(1.16)
4. 改进型自抗扰控制参数整定规律
I-ADRC的性能由各个组成部分的参数决定,在设计了控制器结构之后,需要对其参数进行整定。I-ADRC的三个部分:TD、IESO和INLSEF之间相互独立,互不影响。以下分别分析此三个部分的调节参数对控制系统的影响,并给出参数整定的一般规律。
4.1. 微分跟踪器的参数整定
TD中需整定的参数有速度因子和滤波因子。
令滤波因子
以及给定值
,观察不同速度因子
对跟踪信号
和微分信号
的影响,分别如图5、图6所示。不难看出TD对输入信号的跟踪速度与速度因子
有关。
取值越小,越能够抑制超调现象,但过小的
会导致系统的响应速度过慢;
取值越大,系统的响应速度越快,但过大的
会导致微分信号产生较大波动,降低系统的稳定性。由于
的选取范围很大,在取值时一般只需要对数量级以及最高位数进行调整,可省略其它位上的数值,提高整定效率。
令速度因子
,以及给定值
,观察不同滤波因子
对跟踪信号
的影响,如图7所示。经过对比可知滤波因子
与TD对噪声的抑制性能有关。
取值越大,对噪声的抑制效果越好,但过大的
会严重损失跟踪信号的相位,跟踪速度也会急剧减慢。若输入信号未被污染,只需使滤波因子
小于或等于积分步长h即可,否则需要取更大的
对噪声进行过滤,可参考选取滤波因子
为积分步长h的1~6倍。
4.2. 改进型扩张状态观测器的参数整定
IESO含有的可调参数为分段区间的长度
,非线性因子
、
,观测器的增益系数
、
、
,补偿因子
以及常数
。
补偿因子
的取值与系统的总扰动估计有关,当系统所受扰动较大时,可取较大的
值,通常可在0.01~10区间内选取。
![](//html.hanspub.org/file/26-2571576x102_hanspub.png?20240520094819129)
Figure 5. The impact of different
on
图5. 不同的
对
的影响
![](//html.hanspub.org/file/26-2571576x107_hanspub.png?20240520094819129)
Figure 6. The impact of different
on
图6. 不同的
对
的影响
![](//html.hanspub.org/file/26-2571576x112_hanspub.png?20240520094819129)
Figure 7. The impact of different
on
图7. 不同的
对
的影响
分段区间的长度
的取值大小决定了fal函数线性区间的长度,若
取值过大,会使得线性区间远大于非线性区间,这样fal函数就失去了优越性。若
取值过小,则会使得系统在原点附近存在高频振荡。新的Ifal函数虽然解决了以上问题,但为了便于参数整定,此处
根据积分步长h的大小进行取值。
常数
的取值大小决定了Ifal函数加快减小误差增益的区间范围,若
取值过大,会导致系统在误差较大时,误差增益减小过慢,从而违背“大误差小增益”的原则;若
取值过小,会导致系统在误差较小时,误差增益过早减小,进而违背“小误差大增益”的原则。
可取分段区间长度
的2倍。
非线性因子
、
为Ifal函数的参数增益,不同取值范围的
对Ifal函数的影响效果如图8所示。为了更好的遵从“大误差小增益,小误差大增益”的原则,可在0~1范围内选取非线性因子
,一般可取
,
。
![](//html.hanspub.org/file/26-2571576x132_hanspub.png?20240520094819129)
Figure 8. The influence of
with different value ranges on the Ifal function
图8. 不同取值范围的
对Ifal函数的影响
观测器的增益系数
、
、
分别决定了IESO对系统状态变量x1、x2以及扰动x3的观测精度,因此,
、
、
的取值将直接影响系统的控制精度。若
取值过大,虽然能够抑制对扰动估计的滞后,但同时也会造成估计值的震颤。可通过调整
、
进行抑制,在调节
的同时,调整
、
以得到更好的输出效果。观测器的增益系数一般可通过如式(1.17)所示近似公式进行整定。
(1.17)
4.3. 改进型误差反馈控制律的参数整定
INLSEF含有的可调参数为非线性因子
、
、
,分段区间的长度
,常数
以及比例项、微分项、积分项的权重
、
、
。
非线性因子
、
、
、常数
以及分段区间的长度
分别与IESO中的非线性因子、常数以及分段区间的长度的作用一致,可参考IESO中的整定方法进行整定。
比例项、微分项、积分项的权重
、
、
类似于PID中的
、
、
。
取值越大,过渡过程越快,跟踪的精确度也越高,但会增加超调量。
取值越大,系统的动态性能越好,但会产生噪声。
取值越大,虽然能够提高系统的稳态性能,但会降低系统的动态性能。一般可采取
的调节顺序。
5. 实验验证
为了验证所给出的改进型自抗扰控制的性能,本节利用MATLAB/Simulink搭建如图9所示的基于改进型ADRC的单轴伺服系统直接力控制仿真模型,永磁直线同步电机的参数如表1所示,并给定初级的参考磁链
,正弦跟踪信号
,积分步长
。
针对单轴PMLSM对输入信号的跟踪问题,设计了空载状态下的正弦信号跟随仿真实验;为了验证单轴伺服系统的抗干扰能力,设计了单轴PMLSM变载环境下的仿真实验。为了使得仿真实验更加可靠,将选取传统PID控制器、传统ADRC控制器、非线性函数最外侧采用常数形式的ADRC控制器(C-ADRC) [15] 以及I-ADRC控制器进行仿真对比。
![](//html.hanspub.org/file/26-2571576x174_hanspub.png?20240520094819129)
Figure 9. Block diagram of direct force control for single axis servo system based on improved ADRC
图9.基于改进型ADRC的单轴伺服系统直接力控制框图
![](Images/Table_Tmp.jpg)
Table 1. Parameters of permanent magnet linear synchronous motor
表1. 永磁直线同步电机参数
由于I-ADRC的三个控制模块之间是相互独立的,因此可以根据I-ADRC的参数整定规律进行分开整定,整定后的参数如表2所示。ADRC和C-ADRC中与I-ADRC相同的参数取相同值。PID则采用粒子群算法进行整定,位置环PID:
,
,
;速度环PI:
,
。
5.1. 空载正弦信号跟踪仿真实验
空载状态下跟踪正弦信号的仿真结果如图10所示,整理系统性能指标数据后得到表3。
从图10(a)可以看出,进行对比的四个控制方式均能有效的跟踪正弦输入信号。由图10的(b)、(c)、(d)、(e)以及表3可知,相较于PID控制器,ADRC表现出更小的最大误差、稳态误差,以及更短的稳定时间。其中,I-ADRC在系统启动时的最大误差为210.84 μm,小于传统PID的301.15 μm,比次优的C-ADRC小了6.55 μm,并且I-ADRC在四者中具有最短的稳定时间,仅0.016 s就达到稳态。在系统稳定后,I-ADRC也具有更小的稳态误差,其值为19.55 μm,远小于PID控制的31.18 μm,比次优的C-ADRC小了1.47 μm。因此,在空载条件下,I-ADRC具有更高精度的跟踪性能。
5.2. 变载正弦信号跟踪仿真实验
在第二组变载实验中,系统输入信号与空载实验相同,将PMLSM在空载状态下启动,并在2 s时突
![](//html.hanspub.org/file/26-2571576x206_hanspub.png?20240520094819129)
(a) 位置跟踪响应曲线
![](//html.hanspub.org/file/26-2571576x207_hanspub.png?20240520094819129)
(b) PID位置跟踪误差 (c) ADRC位置跟踪误差![](//html.hanspub.org/file/26-2571576x209_hanspub.png?20240520094819129)
(d) C-ADRC位置跟踪误差 (e) I-ADRC位置跟踪误差
Figure 10. Simulation waveform under no-load conditions
图10. 空载条件下的仿真波形图
![](Images/Table_Tmp.jpg)
Table 3. Motor performance indicators under no-load sine signal input
表3. 空载正弦信号输入下的电机性能指标
![](//html.hanspub.org/file/26-2571576x211_hanspub.png?20240520094819129)
(a) 位置跟踪响应曲线
![](//html.hanspub.org/file/26-2571576x212_hanspub.png?20240520094819129)
(b) PID位置跟踪误差 (c) ADRC位置跟踪误差![](//html.hanspub.org/file/26-2571576x214_hanspub.png?20240520094819129)
(d) C-ADRC位置跟踪误差 (e) I-ADRC位置跟踪误差
Figure 11. Simulation waveform under variable load conditions
图11. 变载条件下的仿真波形图
![](Images/Table_Tmp.jpg)
Table 4. Motor performance indicators under variable load sine signal input
表4. 变载正弦信号输入下的电机性能指标
加50 N的负载,在4 s时突减50 N的负载。仿真结果如图11所示,整理系统性能指标数据后得到表4。
由图11(a)可知负载在2 s以及4 s突变时,四种控制方式未存在明显的位置突变峰值,这表明在正弦信号的变载跟踪过程中,四者均具有一定的抗干扰能力。通过图11的(b)、(c)、(d)、(e)以及表4可知在负载突增时,ADRC的误差突变峰值与调节时间远小于传统PID控制。其中,I-ADRC的误差突变峰值为0.94 μm,小于传统PID控制的8.74 μm,比次优的C-ADRC小了0.38 μm。由于PID控制下的跟踪误差突变较大,其调节到稳态所需的时间大于ADRC,并且I-ADRC的调节时间最短。在负载突减时,与负载突增时具有相类似的情况。因此,单轴伺服系统在对正弦信号进行跟踪时,I-ADRC表现出更好的抗干扰能力。
6. 结论
针对永磁同步直线电机单轴伺服系统的传统直接力控制中PID控制器对系统负载变化的适应性差,抗干扰能力弱,并且无法兼顾响应速度与超调量的问题,提出采用ADRC来代替PID控制器。同时针对传统ADRC中非线性函数在分段点处不平滑的问题,对ESO以及NLSEF进行改进。在空载与变载两种情况下,将I-ADRC与PID以及一些现有的自抗扰控制策略进行仿真对比实验。实验结果表明,相较于传统PID控制器、传统ADRC控制器、C-ADRC控制策略,I-ADRC控制策略有效提高了系统的跟随性能以及抗干扰能力。