Active disturbance rejection control of permanent magnet synchronous motor based on improved faln function
-
摘要:
针对干扰对永磁同步电机控制系统的影响,提出了一种基于改进faln函数的永磁同步电机自抗扰控制方法。设计了永磁同步电机自抗扰控制器,通过分析fal函数对自抗扰控制器的影响,提出了一种限幅的改进faln函数,进而设计了扩张状态观测器,证明了改进后的观测器的收敛性。通过对比试验证明了基于改进faln函数的永磁同步电机自抗扰控制方法相对于传统自抗扰控制具有更好的抗干扰性能。
Abstract:In view of the influence of interference on permanent magnet synchronous motor control system, an improved faln function based active disturbance rejection control for permanent magnet synchronous motor was proposed. An active disturbance rejection controller of permanent magnet synchronous motor is designed. By analyzing the influence of fal function on the active disturbance rejection controller, an improved faln function with limiting amplitude is proposed, and then an extended state observer is designed. The experiment results show that the improved faln function based on permanent magnet synchronous motor auto-rejection control method has better anti-interference performance than the traditional auto-rejection control.
-
永磁同步电机(PMSM)的优点非常明显,其拥有的体积小、功率密度高以及可靠性强等突出优点被广泛应用于汽车行业[1]。电机的控制技术是永磁同步电机能高性能运转的核心,实现好的控制技术不仅与电机本身的结构和参数有关,更与好的控制算法紧密相联[2]。
为了满足现代化工业对控制精度的要求,国内外研究学者提出了很多的控制算法并将其应用于电机的控制系统中,取得了一定的进展。目前针对永磁同步电机的控制算法有PID 控制[3]、自适应控制[4]、滑模控制[5]和智能控制[6]等。永磁同步电机作为一个非线性、不稳定的系统,传统的PID控制器针对其还是具有难以解决的问题,针对不足,韩京清教授首先提出了自抗扰(ADRC)控制器[7]。
王斯墨等[8]将自抗扰技术和负载转矩观测器合并应用于永磁同步电机转速环,设计了改进的二阶自抗扰控制器,有效提高了调速系统的动静态特征。刘鉴等[9]设计了降维转矩观测器与自抗扰技术结合形成基于转矩前馈的自抗扰控制策略,改善了系统动态响应速度。张兴华等[10]将永磁同步电机ia=0的矢量控制系统与自适应反步控制方法组合,设计了一种高性能的电机控制方案,当电机参数改变时,系统具有良好的鲁棒性。李德同等[11]提出了一种变论域双模糊自抗扰控制方法,对于系统内外部不确定性扰动具有良好的鲁棒性。上述方法都是将传统自抗扰技术运用到电机控制中,本文将进一步对自抗扰控制技术内在机理进行分析研究, 提出一种限幅的改进faln函数用在自抗扰控制器中。
本文针对永磁同步电机工作中存在外部不确定干扰的问题,将传统的非线性函数进行改进,证明改进faln 函数的扩张状态观测器的收敛性,通过实验证明永磁同步电机使用本文设计控制器的抗扰性好。
1. 永磁同步电机自抗扰控制器设计
1.1 永磁同步电机数学模型
永磁同步电机凭借着其结构简单、响应快和精度高的优点被广泛使用,根据其结构不同分为表贴式和内嵌式,本文选择表贴式电动机作为研究对象,在不考虑永磁体的阻尼作用和磁滞损耗前提下,PMSM在dq旋转坐标系的数学模型可以表示为
$$ U_d=R_si_d+L_s\frac{\mathrm{d}}{\mathrm{d}t}i_d-\omega_eLi_q $$ (1) $$ U_q=R_si_q+L_s\frac{\mathrm{d}}{\mathrm{d}t}i_q+\omega_e(L_si_d+\psi_f) $$ (2) $$ {T}_{e}=1.5p{\psi }_{f}{i}_{q} $$ (3) $$ T_e-T_l=\eta\frac{\mathrm{d}}{\mathrm{d}t}\omega_m+B\omega_m $$ (4) 式中:$ {U}_{d}\mathrm{、}{U}_{q}\mathrm{、}{i}_{d}\mathrm{、}{i}_{q} $分别是d轴和q轴上的定子电压和定子电流;$ L_s\mathrm{、}\psi_f $和$ R $s分别是定子电感、转子磁链和定子电阻;$ p\mathrm{、}{T}_{l}\mathrm{、}{\omega }_{m}\mathrm{、}\eta \mathrm{、}{\omega }_{e}\mathrm{、}{T}_{e} $和$ B $分别是极对数、负载转矩、机械速度、机器惯性、电角速度、电磁转矩和黏性摩擦因数。
1.2 自抗扰控制器设计
自抗扰控制器(ADRC)由非线性微分跟踪器(tracking differentiator, TD)、扩张观测器(extend-ed state observer, ESO)、 非线性误差反馈控制律(nonlinear law state error feedback, NLSEF)3部分组成,其结构如图1所示。
TD作用是给定信号设置过渡过程并提取信号的微分值;ESO作为自抗扰技术中关键核心,可以观测出状态变量和对系统扰动进行估计并加以补偿;NLSEF将TD和ESO输出进行非线性组合,结合扰动补偿一起构成系统的控制量。本文设计一种简化的一阶自抗扰控制器,一阶ADRC系统具体可表示为
TD:
$$ \left\{\begin{array}{l}{e}_{0}={v}^{*}-{v}_{1}\\ {v}_{1}=-{r}_{0}{f}_{al}({e}_{0},\;{\alpha }_{0},\;{\delta }_{0})\end{array}\right. $$ (5) ESO:
$$ \left\{\begin{array}{l}e={{\textit{z}}}_{1}-y\\ {\dot{{\textit{z}}}}_{1}={{\textit{z}}}_{2}-{\beta }_{1}{f}_{al}\left({e}_{1},\;{\alpha }_{1},\;{\delta }_{1}\right)+{b}_{0}u\left(t\right)\\ {\dot{{\textit{z}}}}_{2}=-{\beta }_{2}{f}_{al}\left({e}_{1},\;{\alpha }_{1},\;{\delta }_{1}\right)\end{array}\right. $$ (6) NLSEF:
$$ \left\{\begin{array}{l}{e}_{1}={v}_{1}-{{\textit{z}}}_{1}\\ {u}_{0}\left(t\right)=k{f}_{al}\left({e}_{2},\;{\alpha }_{2},\;{\delta }_{2}\right)\\ u={u}_{0}\left(t\right)-\dfrac{{{\textit{z}}}_{2}}{{b}_{0}}\end{array}\right. $$ (7) 式中:e0、e、e1为误差信号;v*为TD模块的输入信号;v1为v跟踪信号;r0为速度因子;y为被控对象的输出信号;z1为y的跟踪信号;z2为扰动观测值;α为跟踪因子;δ为滤波因子;β1、β2为ESO输出误差的矫正增益;b0为补偿因子;k为调节器增益;fal为最优控制函数。
$ {f}_{al} $表达式为
$$ f_{al}\left(e,\; \alpha,\; \delta\right)=\left\{\begin{array}{l}\left|e\right|^{\alpha}\mathrm{sin}\ e,\; \left|e\right| > \delta \\ \dfrac{e}{\delta^{1-\alpha}},\; \left|e\right| < \delta\end{array}\right. $$ (8) 在现代电机自抗扰控制领域,选择合适的函数形式以提高电机性能一直是一个重要的研究课题,本文将fal函数应用于永磁同步电机的控制中[12]。这种函数的特点:当误差较小的时候,它可以避免扩张状态观测器产生高频振荡,提高系统的控制精度,从而保证系统的稳定性。
2. fal函数改进及分析
2.1 fal函数改进
永磁同步电机在工作过程中是一个不稳定、非线性的系统,在稳态行驶时误差基本小于1,当行驶在特殊工况时会产生较大的扰动,进而使误差变大,导致系统不稳定。因此本文提出将限幅的改进faln函数应用于永磁同步电机的自抗扰控制器中,改进faln函数表达式见式(9)。
$$ f_{aln}\left(e,\; \alpha,\; \delta\right)=\left\{\begin{array}{l}\left|e\right|^{\alpha}\mathrm{sin}\ e,\; \delta < \left|e\right| < 1 \\ \dfrac{e}{\delta^{1-\alpha}},\; \left|e\right|\leqslant\delta \\ 1,\; \left|e\right|\geqslant1\end{array}\right. $$ (9) 图2所示为$ \delta $取0.01,$ \alpha $分别取0.01、0.15、0.3时fal和faln函数图像。
由图2可得,当误差大于1时,改进faln函数的值被限制在1,明显比fal函数值小,降低了大误差产生的增益量。
2.2 faln函数稳定性分析
分析faln函数的稳定性,首先要知道扩张观测器的误差方程${\dot{e}} $,其表达式为
$$ {\dot{e}}=-A\left(e_1\right)\boldsymbol{e} $$ (10) 式中:e为观测器对状态量的观测误差;e1为对状态量x1的观测误差信号。
要证明使用faln函数的观测器可以完成搭建的永磁同步电机系统的状态观测,引用文献[13]中提出的系统稳定性理论来证明,其引理具体为
引理1 如果存在主对角元素为正数的斜对称矩阵D
$$ {{\boldsymbol{D}}}=\left[\begin{array}{ccc}{{d}}_{{11}}& {{d}}_{{12}}& {{d}}_{{13}}\\ {{-d}}_{{12}}& {{d}}_{{22}}& {{d}}_{{23}}\\ {{-d}}_{{13}}& -{{d}}_{{23}}& {{d}}_{{33}}\end{array}\right] $$ (11) 使DA(e)对称正定,则系统是 Lyapunov渐近稳定的。首先设定了一些参数,例d11=d22=d33=ε>0,来满足主对角线元素为正数的要求。通过一系列计算,得到矩阵DA(e)的表达式,即
$$\boldsymbol{DA} (e)=\left[\begin{array}{lll} D_{11} & -d_{11} & -d_{12} \\ D_{21} & \;\;d_{12} & -d_{22} \\ D_{31} & \;\;d_{31} & \;\;d_{23} \end{array}\right] $$ (12) 其中,
$$ {D}_{11}={d}_{11}{\beta }_{01}+{d}_{12}{\beta }_{02}F+{d}_{13}{\beta }_{03}F $$ (13) $$ {D}_{21}={-d}_{12}{\beta }_{01}+{d}_{22}{\beta }_{02}F+{d}_{23}{\beta }_{03}F $$ (14) $$ {D}_{31}={-d}_{13}{\beta }_{01}-{d}_{23}{\beta }_{02}F+{d}_{33}{\beta }_{03}F $$ (15) $$ F=\frac{{f}_{aln}(e,\alpha ,\delta )}{e} $$ (16) 若DA(e)为对称正定矩阵,需要满足以下条件:
$$ {D}_{11}={\beta }_{01}+\varepsilon {\beta }_{02}F+\varepsilon {\beta }_{03}F > 0 $$ (17) $$ \frac{{\beta }_{01}{\beta }_{02}B+{{\beta }_{02}^{3}}F}{{B}^{2}} > 1 $$ (18) $$ {D}_{11}\left({d}_{12}{d}_{23}-{\varepsilon }^{2}\right)-2{d}_{12}\varepsilon -{{d}_{12}^{3}}-{d}_{23} > 0 $$ (19) 其中,
$$ B=({\beta }_{01}{\beta }_{02}-{\beta }_{03})$$ (20) 若要使上述的不定式成立,需要证明D矩阵存在,在引用引理1时便要考虑矩阵的存在性,而式(17)存在的必要条件为
$$ \mathop {\lim }\nolimits_{\varepsilon \to {0^ + }} {D}_{11}={\beta }_{01}+\frac{{\beta }_{02}^{2}F}{B} > 0 $$ (21) 若式(18)成立,需满足B>0或者是$ {\beta }_{01}+ \dfrac{{\beta }_{02}^{2}F}{B} < 0 $,但由于矩阵DA(e)为对称正定矩阵,与条件$ {\beta }_{01}+\dfrac{{\beta }_{02}^{2}F}{B} < 0 $相互矛盾,因此只需满足B>0这一条件即可。
式(18)和式(19)可以根据矩阵的对称性要求可整理为
$$ \frac{{\beta }_{01}{\beta }_{02}B+{{\beta }_{02}^{3}}F}{{B}^{2}}=1+\frac{{\beta }_{02}^{3}F+{\beta }_{03}B}{{B}^{2}} > 1 $$ (22) $$ {\beta }_{01}{\beta }_{02}-B > 0 $$ (23) 式中:β01、β02、β03为非线性函数的系数。
因为$ B=({\beta }_{01}{\beta }_{02}-{\beta }_{03}) > 0 $,所以以上公式成立。系统状态改变不会影响系统的稳定性,这可以证明使用faln函数的扩张状态观测器搭建的永磁同步电机系统是稳定的。
3. 结果分析
3.1 试验平台
为了验证所提出算法的有效性,本文搭建了图3 所示的试验平台,主要由上位机、RTU-BOX201实时控制器、双向直流电源、驱动模块、PMSM、负载电机、转接板卡等组成。上位机与控制器之间通过以太网通信连接,双向直流电源为系统供电,驱动模块通过控制器控制电机的驱动电压。
将PMSM一端各线与编码器一一对应并与转接卡板连接,另一端接扭矩传感器与负载电机,电机部分参数见表1 。本试验依照参数整定分离原则,将自抗扰控制器3个部分进行分析,通过固定其他参数值,调节单个参数来分析该参数对试验的影响,最终确定自抗扰控制器和改进自抗扰控制器参数为α=0.05,δ=0.01,$ {\mathit{r}}_{0} $=
8000 ,$ {\mathit{b}}_{0} $=720,k=3,$ {\mathit{\beta }}_{0} $=200。表 1 永磁同步电机参数电机参数 数值 磁链/Wb 0.018 极对数 4 电感/mH 0.65 电阻/Ω 0.19 阻尼系数/((N·s)/m) 0.000 1 转动惯量/(kg·m2) 0.000 342 3.2 空载启动
为了验证本文提出的控制方法的有效性,给定永磁同步电机1 000 r/min的额定转速,并使电机处于空载启动状态。当电机运转到达稳态后,记录电机转速和三相电流并绘图,图4所示为使用传统自抗扰控制器和改进自抗扰控制器的电机转速波形图。由图4可以看出,与传统自抗扰控制器相比,使用改进自抗扰控制器的电机转速波形图波动较小。因此在无负载情况下,电机稳态运行时使用改进自抗扰控制器精度更高、稳定性更好。
图5所示为使用传统自抗扰控制器和改进自抗扰控制器的电机三相电流波形图。图5a和图5b分别为使用原fal函数和改进faln函数的三相电流图,为了进一步分析试验结果,对其进行误差分析计算,在转速波动误差的分析结果中,使用传统自抗扰控制器的电机三相电流ia、ib、ic误差分别为4.149、2.186、3.291,而使用改进自抗扰控制器的电机三相电流ia、ib、ic误差则分别为4.121、2.169、3.286。在无外部干扰的情况下,这两种控制策略的运行情况基本一致。但是总的来看,改进后的自抗扰控制器在系统稳定性方面实现了小程度的提升。
3.3 突加负载
为了进一步验证改进后的ADRC具有更出色的抗扰性能,设定1 000 r/min的额定转速,并使电机处于空载启动状态。在达到稳定后,在3 s时施加0.5 N∙m的负载。图6所示为突加0.5 N∙m负载后,转速变化的局部放大图。通过该试验,可以更直观地了解改进ADRC在实际应用中的性能表现,进一步证实其抗扰效果好。
对于外部干扰,改进faln函数和原fal函数设计的自抗扰控制器都能使永磁同步电机恢复稳定运行。由图6可知,采用改进型自抗扰控制器的永磁同步电机转速下降了36 r/min,而使用传统自抗扰控制器的永磁同步电机转速则下降了47 r/min。经过改进的自抗扰控制器波动幅度较小,因此在抗干扰性能方面更好。
以上结果证明,本文设计的改进faln函数自抗扰控制器在永磁同步电机上使用的可行性。相较于传统的自抗扰控制器,使用改进faln函数的自抗扰控制器可以使永磁同步电机具有更好的稳定性,并且抵抗外部干扰的能力更强。
4. 结语
本文的研究重点在于提高使用永磁同步电机在工作时抗干扰能力。基于此目标,本文将使用改进faln函数的自抗扰控制器与永磁同步电机结合。对原fal函数和改进faln函数形式进行了对比分析,并通过数学理论证明了改进faln函数形式的稳定性。进一步通过试验验证基于改进faln函数的永磁同步电机自抗扰控制器在稳定性和抵抗外部干扰能力方面的优势。试验结果表明,改进后的控制器有更好的抗干扰能力。
-
表 1 永磁同步电机参数
电机参数 数值 磁链/Wb 0.018 极对数 4 电感/mH 0.65 电阻/Ω 0.19 阻尼系数/((N·s)/m) 0.000 1 转动惯量/(kg·m2) 0.000 342 -
[1] 王鑫,李伟力,程树康. 永磁同步电动机发展展望[J]. 微电机,2007(5):69-72. DOI: 10.3969/j.issn.1001-6848.2007.05.020 [2] 刘计龙,肖飞,沈洋,等. 永磁同步电机无位置传感器控制技术研究综述[J]. 电工技术学报,2017,32(16):76-88. [3] Xia C L,Ji B N,Yan Y. Smooth speed control for low-speed high-torque permanent-magnet synchronous motor using proportional–integral–resonant controller[J]. IEEE Transactions on Industrial Electronics,2014,62(4):2123-2134.
[4] Liang D ,Li J,Qu R,et al. Adaptive second-order sliding-mode observer for PMSM sensorless control considering VSI nonlinearity[J]. IEEE Transactions on Power Electronics,2017,33(10):8994-9004.
[5] 谭草,鲁应涛,葛文庆,等. 直驱式永磁直线电机深度模糊滑模-自抗扰控制[J]. 西安交通大学学报,2023,57(1):185-194. DOI: 10.7652/xjtuxb202301018 [6] Liu L,Han Y,Fang Y M,et al. Neural network dynamic surface backstepping control for the speedand tension system of reversible cold strip rolling mill[J]. Asian Journal of Control,2018,20(4):1452-1463. DOI: 10.1002/asjc.1628
[7] 韩京清. 自抗扰控制技术[J]. 前沿科学,2007(1):24-31. DOI: 10.3969/j.issn.1673-8128.2007.01.004 [8] 王斯墨,范文超. 改进二阶自抗扰控制器在永磁同步电机上的仿真应用[J]. 火力与指挥控制,2020,45(3):162-166. DOI: 10.3969/j.issn.1002-0640.2020.03.029 [9] 刘鉴,诸德宏,费城. 基于转矩补偿的永磁同步电机自抗扰控制研究[J]. 微电机,2020,53(10):38-41. DOI: 10.3969/j.issn.1001-6848.2020.10.008 [10] 张兴华,唐其太. 考虑参数和负载不确定性的内置式永磁同步电机自适应反步控制[J]. 控制与决策,2016,31(8):1509-1512. [11] 李德同,葛文庆,谭草,等. 电磁直线执行器变论域双模糊自抗扰控制[J]. 制造技术与机床,2023(3):101-106. [12] Zhang Y W,Yang L P,Zhu Y W. Estimate error analysis of the nonlinear third order extended state observer[C]. Intelligent Control and Automation,2012:1621-1627.
[13] Lozgachev G I. On a method of construction of lyapunovfunctions[J]. Automation and Remote Control,1998,59(10):1365-1368.
-
期刊类型引用(1)
1. 姜家锋,李建贵,凌楚冲,徐世康. 基于新型高次谐波滤波器的SPMSM龙伯格观测器. 制造技术与机床. 2024(12): 91-97 . 本站查看
其他类型引用(0)