2. 上海航天设备制造总厂有限公司, 上海 200245
2. Shanghai Aerospace Equipment Manufacturing General Factory, Shanghai 200245, China
相较于传统系统,直驱进给系统具有传动误差小、能耗低、噪声小、定位精度高以及动态响应性能好等优点,近年来在高速高精密加工领域得到广泛应用[1-5]。
传统进给系统存在机械传动结构,其系统刚度主要取决于传动部分机械刚度;直驱进给系统省去机械传动环节,其系统刚度主要取决于伺服刚度[6],由于负载扰动无缓冲地作用在直驱进给系统上,同时系统运行时受非线性因素影响,因此伺服控制具有较高难度[7]。对于直驱进给系统,不合理的伺服参数设置会导致系统抗干扰能力差、精度低以及系统动态特性差等问题,影响工件加工精度和表面质量[8]。伺服刚度优化研究对提高直驱式伺服进给系统动态特性具有重要意义。
目前,关于直驱式伺服进给系统的伺服刚度优化研究受到了国内外学者的广泛关注。相关研究可大致分为3类:1) 对于直驱式伺服进给系统结构建模优化及部件匹配优化[9-10];2) 基于具体优化目标,如跟随性能和伺服动刚度等,给出针对性的电机控制算法,通过算法实时自动调整伺服控制参数对系统进行优化[11-12];3) 通过各种方式开展不同的频率响应实验,以获得实际机械系统动刚度以及各控制参数对系统动刚度的影响[13-14]。
然而,目前对直驱式伺服进给系统伺服刚度优化的研究大多基于单个参数对系统伺服刚度的影响,给出相应的优化方法。本文基于多个参数对伺服动刚度的影响,对系统Routh判稳条件进行化简得到多个参数之间的约束关系,提出了一种同时调整多个伺服控制参数以提高系统伺服动刚度的优化方法。该优化方法对于实际工程应用中提高直驱式伺服进给系统动刚度具有指导意义,并为后续寻求基于伺服动刚度的全局最优控制参数整定方案的理论研究奠定了基础。
1 力矩电机直驱伺服进给系统机电耦合建模 1.1 力矩电机建模目前数控机床直驱进给系统多采用环形永磁力矩电机,是一种永磁同步电机(permanent magnet synchronous motor, PMSM)。该类电机建模复杂,控制参数整定困难。PMSM电流环在以控制结构为重点研究对象时可等效为直流电机电流环,其具体等效原理本文中不再赘述,其物理模型如图 1所示[15]。
图 1中:ei(t)为电机输入电压,La为电枢电感,Ra为电枢电阻,em(t)为电机工作时产生的反电动势,ia(t)为回路电流,T(t)为电机输出转矩, θ0(t)为电机输出转角, Je为系统及负载在电机轴上的等效转动惯量, MC(t)为负载转矩。
由Kirchhoff定律及电磁感应定律有:
$ \begin{gathered} e_{\mathrm{i}}(t)=R_{\mathrm{a}} i_{\mathrm{a}}(t)+L_{\mathrm{a}} \frac{\mathrm{d} i_{\mathrm{a}}(t)}{\mathrm{d} t}+e_{\mathrm{m}}(t), \\ e_{\mathrm{m}}(t)=K_{\mathrm{e}} \frac{\mathrm{d} \theta_{0}(t)}{\mathrm{d} t}. \end{gathered} $ | (1) |
其中Ke为电机反电动势系数。
根据力矩平衡方程以及Lorentz电磁定律有:
$ T(t)-M_{\mathrm{C}}(t)-D_{\mathrm{m}} \frac{\mathrm{d} \theta_{0}(t)}{\mathrm{d} t}=J_{\mathrm{e}} \frac{\mathrm{d}^{2} \theta_{0}(t)}{\mathrm{d} t^{2}}, $ | (2) |
$ T(t)=K_{\mathrm{t}} i_{\mathrm{a}}(t). $ | (3) |
其中:Kt为电动机力矩系数;Dm为系统在电机轴上的等效黏滞阻尼系数。联立式(2)和(3)并进行Laplace变换可得
$ \begin{gathered} \varTheta_{0}(s)=\frac{K_{\mathrm{t}} E_{\mathrm{i}}(s)}{L_{\mathrm{a}} J_{\mathrm{e}} s^{3}+R_{\mathrm{a}} J_{\mathrm{e}} s^{2}+K_{\mathrm{t}} K_{\mathrm{e}} s}- \\ \frac{\left(L_{\mathrm{a}} s+R_{\mathrm{a}}\right) M_{\mathrm{C}}(s)}{L_{\mathrm{a}} J_{\mathrm{e}} s^{3}+R_{\mathrm{a}} J_{\mathrm{e}} s^{2}+K_{\mathrm{t}} K_{\mathrm{e}}s}. \end{gathered} $ | (4) |
由式(4)可得电机系统模型框图如图 2所示。由式(4)不难看出,电机输出角度取决于Ei(s)与MC(s)两部分系统输入,2个输入变量与输出分别构成单输入-单输出系统,且2个系统均为典型的二阶振荡系统。
1.2 伺服进给控制系统建模
如图 3所示,根据图 2并加入电流环、速度环和位置环后即可建立电机的控制模型,电机控制模型有2个输入,即输入电压Ei(s)和负载力矩MC(s),输出为电机转角Θ0(s)。采用比例-积分-微分(proportion-integration-differentiation, PID)控制器,采用PI控制的电流环和速度环,以及P控制的位置环,就可以保证控制系统具有较好的稳定性、较高的精度以及较快的响应。
从电流环开始,以输出角速度作为输出,以输入电流作为输入,从内向外解算该控制系统的传递函数。解算各环的推导过程在文[15]中已给出,在此仅给出最后得到的系统整体传递函数为
$ \varTheta_{0}(s)=\varTheta_{\mathrm{i}}(s) G_{\mathrm{p}}(s)-M_{\mathrm{C}}(s) G_{\mathrm{M}}(s). $ | (5) |
其中:
$ \begin{gathered} G_{\mathrm{p}}(s)=\left[K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} K_{\mathrm{pp}}\left(T_{\mathrm{ii}} T_{\mathrm{iv}} s^{2}+\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right) s+1\right)\right] / \\ \quad\left\{T_{\mathrm{iv}} J_{\mathrm{e}} T_{\mathrm{ii}} L_{\mathrm{a}} s^{5}+T_{\mathrm{iv}}\left[J_{\mathrm{e}} T_{\mathrm{ii}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)+D_{\mathrm{m}} L_{\mathrm{a}} T_{\mathrm{ii}}\right] s^{4}+\right. \\ T_{\mathrm{iv}}\left(K_{\mathrm{t}} K_{\mathrm{e}} T_{\mathrm{ii}}+J_{\mathrm{e}} K_{\mathrm{pi}}+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{ii}}+D_{\mathrm{m}} R_{\mathrm{a}} T_{\mathrm{ii}}+K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}}\right) s^{3}+ \\ \left(K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{iv}}+K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}} T_{\mathrm{iv}}\right) s^{2}+ \\ \left.K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(K_{\mathrm{pp}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+1\right) s+K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\right\}. \end{gathered} $ |
$ \begin{gathered} G_{\mathrm{M}}(s)=\left[T_{\mathrm{ii}} L_{\mathrm{a}} T_{\mathrm{iv}} s^{3}+T_{\mathrm{ii}} T_{\mathrm{iv}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right) s^{2}+K_{\mathrm{pi}} T_{\mathrm{iv}} s\right] / \\ \left\{T_{\mathrm{iv}} J_{\mathrm{e}} T_{\mathrm{ii}} L_{\mathrm{a}} s^{5}+T_{\mathrm{iv}}\left[J_{\mathrm{i}} T_{\mathrm{ii}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)+D_{\mathrm{m}} L_{\mathrm{a}} T_{\mathrm{ii}}\right] s^{4}+\right. \\ T_{\mathrm{iv}}\left(K_{\mathrm{t}} K_{\mathrm{e}} T_{\mathrm{ii}}+J_{\mathrm{e}} K_{\mathrm{pi}}+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{ii}}+D_{\mathrm{m}} R_{\mathrm{a}} T_{\mathrm{ii}}+K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}}\right) s^{3}+ \\ \left(K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{iv}}+K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}} T_{\mathrm{i}}\right) s^{2}+ \\ \left.K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(K_{\mathrm{pp}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+1\right) s+K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\right\}. \end{gathered} $ |
式(5)给出了以完整的位置指令与干扰力矩为双输入、以电机实际转角为输出的双输入单输出PID伺服进给控制系统表达式。
2 伺服控制参数对伺服刚度的影响 2.1 伺服刚度的定义对于任一系统其刚度具有各向异性,本文所研究的刚度指A、C双摆头系统圆周进给方向的刚度。在后文中若未加特别说明,均指该圆周进给方向刚度。
直驱式力矩电机机床伺服系统抵抗外界干扰引起输出偏差的能力取决于伺服控制系统以及A、C双摆头组成的系统整体所体现出的在外载荷作用下保持位置不产生偏差的能力,即伺服刚度。伺服刚度单位为N·m/rad,又分为静刚度与动刚度。
伺服静刚度的定义为伺服电机为消除位置偏差而产生的输出力矩与位置偏差之比,即在恒定外载荷作用下进给驱动系统抵抗位置偏差的能力[12],
$ K_{\mathrm{sj}}=\frac{T_{\mathrm{dis}}}{\Delta \theta}. $ | (6) |
其中:Tdis为系统所受干扰力矩,Δθ为干扰力矩引起的伺服系统输出转角增量。
伺服动刚度的定义为在动态载荷下动态干扰力矩与伺服进给系统为消除干扰所输出的转角增量之比,即在动态外载荷作用下伺服进给系统抵抗位置偏差的能力[12],
$ K_{\mathrm{sd}}(\mathrm{j} \omega)=\frac{T_{\mathrm{dis}}(\mathrm{j} \omega)}{\Delta \theta(\mathrm{j} \omega)}. $ | (7) |
其中:Tdis(jω)为动态干扰载荷,Δθ(jω)为随动态干扰载荷变化的伺服系统输出转角增量。
因Ksd(jω)为复数,故用其幅值衡量动刚度大小,
$ K_{\mathrm{sd}}=\left|K_{\mathrm{sd}}(\mathrm{j} \omega)\right| \text {. } $ | (8) |
首先验证节1.2中所建伺服进给系统控制模型在描述系统伺服刚度方面的有效性及准确性。依据文[16]中用于仿真分析及实验验证的伺服系统控制参数及电机结构参数,通过节1.2中图 3所示模型得到各参数对系统伺服刚度的影响,如图 4—8所示。
![]() |
图 4 位置环增益Kpp对伺服动刚度的影响 |
![]() |
图 5 速度环增益Kpv对伺服动刚度的影响 |
![]() |
图 6 电流环增益Kpi对伺服动刚度的影响 |
![]() |
图 7 速度环时间常数Tⅳ对伺服动刚度的影响 |
![]() |
图 8 电流环时间常数Tii对伺服动刚度的影响 |
将得到的各参数对伺服动刚度影响的仿真曲线与文[16]中的实验结果进行对比,可验证本文节1.2所建立的直驱式伺服进给系统模型在描述伺服刚度方面的合理性、有效性以及较高的准确性。
3 各伺服控制参数约束关系分析本节从伺服控制系统稳定性的整体角度,给出了各控制参数之间的约束关系,并提出了以伺服动刚度为优化目标的伺服系统参数优化策略。
3.1 基于Routh判据的各参数约束关系分析完整的三闭环PID伺服控制系统十分复杂。考虑到力矩电机实际结构、初步三环PID控制器控制参数整定及节2.2中各系统参数对伺服刚度的影响,为使基于Routh判据的各参数约束关系分析能对于本文后续基于伺服刚度的伺服系统参数优化具有更实际的指导意义,作出如下假设:
1) 由于实际电机结构中电枢电感La及电枢电阻Ra较小,与除速度环积分时间常数Tⅳ及电流环积分时间常数Tii外其他的系统参数相比要小多个量级,因此在化简时将电枢电感La及电枢电阻Ra作为小量处理;
2) 由节2.2中速度环积分时间常数Tⅳ与电流环积分时间常数Tii分别对伺服动刚度的影响曲线(图 7和8)可以看出,两者越大,系统的伺服动刚度越小;且在系统三环PID控制器参数初步整定时Tⅳ和Tii一般均设定较小,与除电枢电感La及电枢电阻Ra外其他的系统参数相比要小多个量级,因此在化简时将速度环积分时间常数Tⅳ与电流环积分时间常数Tii作为小量处理。
设H为系统传递函数的特征多项式,由式(5)可得
$ \begin{gathered} H=T_{\mathrm{iv}} J_{\mathrm{e}} T_{\mathrm{ii}} L_{\mathrm{a}} s^{5}+ \\ T_{\mathrm{iv}}\left[J_{\mathrm{e}} T_{\mathrm{ii}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)+D_{\mathrm{m}} L_{\mathrm{a}} T_{\mathrm{ii}}\right] s^{4}+ \\ T_{\mathrm{iv}}\left(K_{\mathrm{t}} K_{\mathrm{e}} T_{\mathrm{ii}}+J_{\mathrm{e}} K_{\mathrm{pi}}+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{ii}}+D_{\mathrm{m}} R_{\mathrm{a}} T_{\mathrm{ii}}+\right. \\ \left.K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}}\right) s^{3}+\left(K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+\right. \\ \left.D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{iv}}+K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}} T_{\mathrm{iv}}\right) s^{2}+ \\ K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(K_{\mathrm{pp}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+1\right) s+K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} . \end{gathered} $ |
H可化为
$ H=a_{0} s^{5}+a_{1} s^{4}+a_{2} s^{3}+a_{3} s^{2}+a_{4} s+a_{5} \text {. } $ | (9) |
其中:
$ a_{0}=T_{\mathrm{iv}} J_{\mathrm{e}} T_{\mathrm{ii}} L_{\mathrm{a}}, $ |
$ a_{1}=T_{\mathrm{iv}}\left[J_{\mathrm{e}} T_{\mathrm{ii}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)+D_{\mathrm{m}} L_{\mathrm{a}} T_{\mathrm{ii}}\right], $ |
$ \begin{gathered} a_{2}=T_{\mathrm{iv}}\left(K_{\mathrm{t}} K_{\mathrm{e}} T_{\mathrm{ii}}+J_{\mathrm{e}} K_{\mathrm{pi}}+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{ii}}+\right. \\ \left.D_{\mathrm{m}} R_{\mathrm{a}} T_{\mathrm{ii}}+K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}}\right), \end{gathered} $ |
$ \begin{gathered} a_{3}=K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+D_{\mathrm{m}} K_{\mathrm{pi}} T_{\mathrm{iv}}+ \\ K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}} T_{\mathrm{iv}}, \end{gathered} $ |
$ \begin{gathered} a_{4}=K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}}\left(K_{\mathrm{pp}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)+1\right), \\ a_{5}=K_{\mathrm{pp}} K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} . \end{gathered} $ |
由Routh判据可得,该闭环系统稳定的充分必要条件为:
$ a_{0}>0, $ | (10) |
$ a_{1}>0, $ | (11) |
$ \frac{a_{1} a_{2}-a_{0} a_{3}}{a_{1}}>0, $ | (12) |
$ \frac{a_{1} a_{2} a_{3}-a_{2} a_{3}^{2}-a_{1}^{2} a_{4}+a_{0} a_{1} a_{5}}{a_{1} a_{2}-a_{0} a_{3}}>0, $ | (13) |
$ \frac{\left(a_{1} a_{2} a_{3}-a_{0} a_{3}^{2}-a_{1}^{2} a_{4}+a_{0} a_{1} a_{5}\right)\left(a_{1} a_{4}-a_{0} a_{5}\right)-\left(a_{1} a_{2}-a_{0} a_{3}\right)^{2} a_{5}}{\left(a_{1} a_{2}-a_{0} a_{3}\right) a_{1} a_{3}-a_{1}^{3} a_{4}+a_{0} a_{1}^{2} a_{5}}>0 , $ | (14) |
$ a_{5}>0 . $ | (15) |
由伺服进给系统性质可知,所涉及的各参数均为正数,因此可知式(10)、(11)及(15)均成立。
在式(11)成立的条件下,式(12)两边同乘a1后可化为
$ a_{1} a_{2}>a_{0} a_{3}. $ |
上式两边同时除以a0并代入式(9)中各阶次系数,可得
$ \begin{gathered} T_{\mathrm{ii}} T_{\mathrm{iv}}+\frac{D_{\mathrm{m}} T_{\mathrm{iv}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pp}}}+\frac{T_{\mathrm{ii}}+T_{\mathrm{iv}}}{K_{\mathrm{pp}}}> \\ \left(\frac{K_{\mathrm{pi}}+R_{\mathrm{a}}}{L_{\mathrm{a}}}+\frac{D_{\mathrm{m}}}{J_{\mathrm{e}}}\right)\left(\frac{K_{\mathrm{e}} T_{\mathrm{ii}} T_{\mathrm{iv}}}{K_{\mathrm{pi}} K_{\mathrm{pv}} K_{\mathrm{pp}}}+\frac{J_{\mathrm{e}} T_{\mathrm{iv}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pp}}}+\right. \\ \left.\frac{D_{\mathrm{m}} T_{\mathrm{ii}} T_{\mathrm{iv}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pp}}}+\frac{D_{\mathrm{m}} R_{\mathrm{a}} T_{\mathrm{ii}} T_{\mathrm{iv}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pi}} K_{\mathrm{pp}}}+\frac{T_{\mathrm{ii}} T_{\mathrm{iv}}}{K_{\mathrm{pp}}}\right). \end{gathered} $ |
在上面不等式的右式中,忽略高阶小量并进行合理放缩,不等式可化为
$ \begin{aligned} &\left(\frac{K_{\mathrm{pi}}+R_{\mathrm{a}}}{L_{\mathrm{a}}}+\frac{D_{\mathrm{m}}}{J_{\mathrm{e}}}\right) \frac{J_{\mathrm{e}} T_{\mathrm{iv}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pp}}}> \\ &T_{\mathrm{ii}} T_{\mathrm{iv}}+\frac{D_{\mathrm{m}} T_{\mathrm{iv}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pp}}}+\frac{T_{\mathrm{ii}}+T_{\mathrm{iv}}}{K_{\mathrm{pp}}}. \end{aligned} $ | (16) |
对不等式(16)进行移项整理后可得
$ \frac{J_{\mathrm{e}} K_{\mathrm{pi}}+K_{\mathrm{t}} K_{\mathrm{pi}} K_{\mathrm{pv}} T_{\mathrm{ii}}}{L_{\mathrm{a}} K_{\mathrm{t}} T_{\mathrm{ii}} K_{\mathrm{pv}}}-\frac{T_{\mathrm{ii}}+T_{\mathrm{iv}}}{T_{\mathrm{ii}} T_{\mathrm{iv}}}>K_{\mathrm{pp}}. $ | (17) |
由式(17)不难看出,存在系统各参数的取值使得式(17)成立。在式(12)成立的条件下,不等式(13)两边同乘a1a2-a0a3并移项可得
$ a_{1} a_{2} a_{3}+a_{0} a_{1} a_{5}>a_{2} a_{3}^{2}+a_{1}^{2} a_{4}. $ | (18) |
考虑式(18)中各项在代入各阶次系数后高阶小量的阶次,式(18)右式中消去a2a32并将不等式两边同时除以a1可得
$ a_{2} a_{3}+a_{0} a_{5}>a_{1} a_{4} $ | (19) |
代入各阶次系数,不等式(19)两边同时除以JeKtKpiKpvKpp可得
$ \begin{gathered} \frac{K_{\mathrm{pi}} D_{\mathrm{m}} T_{\mathrm{ii}}}{K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pp}}}+\frac{K_{\mathrm{pi}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)}{K_{\mathrm{pp}}}+K_{\mathrm{pi}} T_{\mathrm{ii}} T_{\mathrm{iv}}+L_{\mathrm{a}} T_{\mathrm{ii}}> \\ T_{\mathrm{ii}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)\left(\frac{1}{K_{\mathrm{pp}}}+T_{\mathrm{ii}}+T_{\mathrm{iv}}\right) . \end{gathered} $ |
忽略高阶小量,式(19)移项整理后可简化为
$ \frac{K_{\mathrm{pi}} T_{\mathrm{iv}}}{T_{\mathrm{ii}}\left(K_{\mathrm{pi}} T_{\mathrm{ii}}-L_{\mathrm{a}}\right)}>K_{\mathrm{pp}}. $ | (20) |
分析式(17)及(20)不难得出,存在系统各参数的取值使得两式同时成立。在式(11)、(17)及(20)同时成立的条件下,式(14)两边同时乘以(a1a2-a0a3)a1a3-a13a4+a0a12a5,移项整理后可得
$ \begin{gathered} \left(a_{1} a_{2} a_{3}-a_{0} a_{3}^{2}-a_{1}^{2} a_{4}+a_{0} a_{1} a_{5}\right)\left(a_{1} a_{4}-a_{0} a_{5}\right)> \\ \left(a_{1} a_{2}-a_{0} a_{3}\right)^{2} a_{5}. \end{gathered} $ |
去括号后整理不等式有
$ \begin{gathered} a_{1}^{2} a_{2} a_{3} a_{4}+2 a_{0} a_{1}^{2} a_{4} a_{5}+a_{0} a_{1} a_{2} a_{3} a_{5}-a_{0} a_{1} a_{3}^{2} a_{4}- \\ a_{1}^{3} a_{4}^{2}-a_{0}^{2} a_{1} a_{5}^{2}>a_{1}^{2} a_{2}^{2} a_{5}. \end{gathered} $ |
考虑上式中各项在代入各阶次系数后高阶小量的阶次,式左消去a0a1a32a4,不等式两边同时除以a0a1a5可得
$ \frac{a_{1} a_{2} a_{3} a_{4}}{a_{0} a_{5}}+a_{2} a_{3}>\frac{a_{1} a_{2}^{2}}{a_{0}}+\left(\frac{a_{1} a_{4}}{a_{0} a_{5}}-2\right) a_{1} a_{4}+a_{0} a_{5} . $ | (21) |
代入各阶次系数,忽略高阶小量,不等式(21)左边可化为
$ \frac{J_{\mathrm{e}} K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pi}}^{2} T_{\mathrm{iv}}}{L_{\mathrm{a}}}\left[\frac{K_{\mathrm{pi}}}{K_{\mathrm{pp}}}+K_{\mathrm{pi}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)\right]\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right), $ |
不等式(21)右边可化为
$ \begin{gathered} \frac{J_{\mathrm{e}} K_{\mathrm{t}} K_{\mathrm{pv}} K_{\mathrm{pi}}^{2} T_{\mathrm{iv}}}{L_{\mathrm{a}}}\left[\frac{J_{\mathrm{e}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)+D_{\mathrm{m}} L_{\mathrm{a}}}{K_{\mathrm{t}} K_{\mathrm{pv}}}+\right. \\ \left.\left[\frac{K_{\mathrm{pi}}}{K_{\mathrm{pp}}}+K_{\mathrm{pi}}\left(T_{\mathrm{ii}}+T_{\mathrm{iv}}\right)\right] T_{\mathrm{i}}\right]. \end{gathered} $ |
综上,有
$ \frac{K_{\mathrm{pi}}}{K_{\mathrm{pp}}}+K_{\mathrm{pi}} T_{\mathrm{iv}}>\frac{J_{\mathrm{e}}\left(K_{\mathrm{pi}}+R_{\mathrm{a}}\right)+D_{\mathrm{m}} L_{\mathrm{a}}}{K_{\mathrm{t}} K_{\mathrm{pv}}}. $ | (22) |
分析式(17)、(20)及(22)不难得出,存在系统各参数的取值使得这3个公式同时成立,即存在系统各参数取值使得该系统稳定。式(17)、(20)及(22)给出了在系统稳定条件下各参数之间的约束关系。
3.2 系统稳定性判据验证以某一装有A、C双摆头的五轴联动机床中A轴伺服进给系统为例。该系统初步整定的各伺服控制参数及电机结构参数如表 1所示。
参数 | 物理意义 | 单位 | 值 |
Kpp | 位置环比例增益 | s-1 | 20.851 |
Kpv | 速度环比例增益 | A·s/rad | 30.257 |
Kpi | 电流环比例增益 | V/A | 10.521 |
Tⅳ | 速度环积分时间常数 | s | 0.006 |
Tii | 电流环积分时间常数 | s | 0.002 |
Kt | 电机力矩常数 | N·m/A | 30 |
Ke | 电机反电动势常数 | V·s/rad | 18.52 |
La | 电枢绕组电感 | H | 0.003 5 |
Ra | 电枢绕组电阻 | Ω | 0.052 |
Je | 系统等效转动惯量 | kg·m2 | 20 |
Dm | 系统等效黏滞阻尼系数 | kg·m2/s | 0.3 |
基于表 1中实际系统的各参数值对节3.1中所得Routh判据表达式的准确性进行验证。基于节3.1中在进行Routh判据表达式简化前所作假设,在验证时不妨令Tⅳ、Tii及其他电机结构参数保持不变;为验证理论模型的准确性,暂不考虑系统实际工作时电机力矩、转速及电流等限制。
对于绝大多数采用PID控制的数控机床进给系统,其各环比例增益的整定范围均在0~200之间(各比例增益数值单位同表 1),而分析式(17)与(20)不难看出,当各环比例增益在该范围内变化时,式(17)及(20)始终成立。
因此,本节仅给出式(22)的准确性验证过程。保持Kpi不变,令Kpv=50 A·s/rad以使得在验证式(22)时式(17)及(20)始终成立。由式(22)可得,当Kpp>136.364 s-1时系统失稳即式(22)不再成立。
逐渐增大Kpp并通过Simulink得到系统阶跃响应仿真曲线,得到当Kpp=137.5 s-1时,系统响应出现振荡失稳,仿真曲线如图 9所示。综上,理论模型中式(22)的准确性得到仿真验证。
![]() |
图 9 Kpp=137.5 s-1、Kpv=50 A·s/rad时系统阶跃响应曲线 |
4 基于Routh判据的伺服刚度优化 4.1 伺服系统多参数优化方法
由节2.2中伺服系统各参数对伺服刚度的影响曲线可知,减小Tⅳ及Tii可提高伺服动刚度;由PID控制原理可知,增大Tⅳ及Tii会极大地影响系统的动态响应性能。综合考虑以上两方面,结合实际工程应用中对于伺服系统的性能要求,在对Tⅳ及Tii进行优化时供优化搜索的范围很小,导致其数值变化对于伺服刚度及系统动态响应性能影响较小,因此在优化中暂不考虑对Tⅳ及Tii的优化整定,而是将它们设成较小的常数值。
在对伺服系统中各环比例增益进行优化整定时,将Kpp、Kpv、Kpi的优化搜索范围设置为Kpp∈ (0, kpp lim], Kpv∈(0, kpv lim], Kpi∈(0, kpi lim]。各增益数值单位与表 1中相同;其中,kpp lim、kpv lim以及kpi lim的值取决于系统给定输入以及实际电机运行时电机的最大转速、最大电流和最大输出转矩。
首先,取出Kpv搜索范围内所有整数,
$ K_{\mathrm{pv}_{i}}=i \quad\left(i \in\left(0, k_{\mathrm{pv} \lim }\right], i \in {\bf{Z}}^{+}\right). $ | (23) |
将Kpvi代入式(20)和(22)中得到在系统稳定条件下各Kpvi对应的位置环比例增益Kpp取值范围,将该范围与Kpp的优化搜索范围取交集后取整,与Kpv同理,即
$ K_{\mathrm{pp}_{i j}}=j \quad\left(j \in\left(0, k_{\mathrm{pp} \lim }^{\prime}\right], j \in {\bf{Z}}^{+}\right). $ | (24) |
其中kpp lim′为取交集后新的取值上限。为缩短优化方法计算时间,考虑到节2.2中所得位置环比例增益Kpp增大到一定程度后对伺服刚度的影响变化不大,将式(24)改写为
$ \left\{\begin{array}{l} K_{\mathrm{pp}_{i j}}=j \\ \ \ \ \ \ \ \ \ \left(j \in\left(0, \frac{k_{\mathrm{pp\lim }}^{\prime} }{2}\right], j \in {\bf{Z}}^{+}\right), \\ K_{\mathrm{pp}_{i j}}=2 j-\left\lceil {\frac{k_{\mathrm{pp\lim}}^{\prime}}{2}} \right\rceil \\ \ \ \ \ \ \ \ \ \left(j \in\left(\frac{k_{\mathrm{pp \lim}}^{\prime} }{2}, \frac{3 k_{\mathrm{pp} \lim}^{\prime}}{4}\right], j \in {\bf{Z}}^{+}\right). \end{array}\right. $ | (25) |
其中
将Kpvi及其对应的j个Kppij分别代入式(17)中得到系统稳定条件下各组位置环比例增益Kpp与速度环比例增益Kpv对应的电流环比例增益取值范围。由节2.2可知Kpi对于伺服刚度的影响与Kpp相似,因此为缩短优化方法计算时间,对Kpi的数据处理方式与Kpp相同,可得
$ \left\{\begin{array}{l} K_{\mathrm{pi}_{i j h}}=\left\lceil h+k_{\mathrm{pi\min}}\right\rceil \\ \left(h \in\left(0, \frac{k_{\mathrm{pi} \lim}-k_{\mathrm{pimin}}}{2}\right], h \in {\bf{Z}}^{+}\right), \\ K_{\mathrm{pi}_{i j h}}=2\left\lceil h+k_{\mathrm{pi\min}}\right\rceil-\left\lceil\frac{k_{\mathrm{pi\min}}+k_{\mathrm{pi\lim}}}{2}\right\rceil \\ \left(h \in\left(\frac{k_{\mathrm{pi} \min }+k_{\mathrm{pi\lim}}}{2}, \frac{3 k_{\mathrm{pi\min}}+3 k_{\mathrm{pi\lim}}}{4}\right], h \in {\bf{Z}}^{+}\right). \end{array}\right. $ |
综上,即可在精度为1的条件下遍历在系统稳定前提下各环比例增益在优化搜索范围内的所有组合。
选取伺服动刚度的倒数即伺服动柔度作为伺服动刚度高低的评价指标。伺服动柔度的表达式即为干扰力矩输入(即电机实际转角输出)的系统传递函数,也就是式(5)中的GM(s)。伺服动柔度的最大值对应于伺服动刚度的最小值。除此之外,为避免优化结果一味地增大各环比例增益而导致系统出现超调过大、振荡频繁等诸多动态性能问题,选取系统阶跃响应的调节时间Tstable为另一评价指标。综上,采用式(26)作为目标优化函数,
$ f=\frac{1}{\frac{\lambda}{K_{\mathrm{sd}}}+T_{\text {stable }}}. $ | (26) |
其中λ为根据系统实际性能需求设定的权重系数,为常数。
通过以上步骤即可得到精度为1的全局最优解,此时再根据所需精度,通过搜索区间长度为2的局部寻优,即可得到精度更高的解。可见,本节采用的优化方法在优化整定各环比例增益以提高了伺服刚度的同时,对系统的动态性能影响较小。
4.2 伺服系统动刚度优化方法验证基于表 1中的系统初始参数作Bode图,如图 10所示,得到优化前伺服动柔度的频域特性。从图 10中可知,优化前系统的固有频率约88 rad/s,动柔度振幅峰值约为-90 dB。
![]() |
图 10 优化前伺服动柔度Bode图 |
该系统采用的力矩电机额定转矩为6 000 N·m,额定转速为200 r/min,额定电流为200 A。将表 1中除伺服控制参数Kpp、Kpv、Kpi外的各参数代入式(17)、(20)及(22)中可得该系统的稳定性条件。
由于在实际加工过程对该机床A轴伺服控制系统的性能要求是动态响应快、跟随误差小且抵抗切削力干扰的能力强,因此在输入为0.1 rad的阶跃响应下,设定权重系数为1010。
设定整定精度为0.1。通过节4.1中方法优化后伺服系统动柔度的频域特性如图 11所示,系统伺服动刚度得到明显提高,固有频率约127 rad/s,动柔度振幅峰值约为-110 dB,且全频域内整体伺服动刚度均有所提高。此时,Kpp=25.3 s-1,Kpv=65.2 A·s/rad,Kpi=40.5 V/A,其他参数与初始值相同。可见,本文优化方法对于系统伺服刚度的优化效果较好。
![]() |
图 11 优化后伺服动柔度Bode图 |
5 结论
本文提出了一种对采用力矩电机直驱的数控机床进给系统的伺服刚度进行优化的方法。首先通过电流环等效简化了力矩电机模型,并建立了完整的伺服进给系统机电耦合模型;根据伺服系统各环PID控制参数对伺服刚度的影响,对系统稳定性Routh判据不等式进行了简化求解,得到各控制参数间的约束关系式并通过Simulink仿真验证了关系式的准确性;基于上述分析结果,提出了以伺服刚度为优化目标且能同时考虑系统调节时间的伺服控制参数优化方法。最后,以某五轴联动数控机床A轴力矩电机直驱进给系统为例,进行了优化应用验证,Bode图结果表明本文优化方法能有效提高系统伺服刚度,且操作难度低,易于实现,对实际工程实践具有指导意义。
[1] |
王国彪, 赖一楠, 范大鹏, 等. 新型精密传动机构设计与制造综述[J]. 中国机械工程, 2010, 21(16): 1891-1897. WANG G B, LAI Y N, FAN D P, et al. Summary of new type precision transmission design and manufacture[J]. China Mechanical Engineering, 2010, 21(16): 1891-1897. (in Chinese) |
[2] |
邓朝晖, 刘战强, 张晓红. 高速高效加工领域科学技术发展研究[J]. 机械工程学报, 2010, 46(23): 106-120. DENG Z H, LIU Z Q, ZHANG X H. Research of the science and technology development in high-speed and efficient processing field[J]. Journal of Mechanical Engineering, 2010, 46(23): 106-120. (in Chinese) |
[3] |
ALTINTAS Y, VERL A, BRECHER C, et al. Machine tool feed drives[J]. CIRP Annals, 2011, 60(2): 779-796. |
[4] |
HOSSEINI M S, KHAMENEH M M, VAEZ-ZADEH S. Design of thrust ripple minimization in wound secondary linear synchronous motors by response surface methodology (RSM)[C]//2011 2nd Power Electronics, Drive Systems and Technologies Conference. Tehran, Iran, 201l: 56-61.
|
[5] |
TAVANA N R, SHOULAIE A. Pole-shape optimization of permanent-magnet linear synchronous motor for reduction of thrust ripple[J]. Energy Conversion and Management, 2011, 52(1): 349-354. DOI:10.1016/j.enconman.2010.07.007 |
[6] |
徐进宝. 直接驱动数控转台伺服刚度研究[D]. 沈阳: 沈阳工业大学, 2006. XU J B. The stiffness analysis of direct drive NC machine tool rotary table servo system[D]. Shenyang: Shenyang University of Technology, 2006. (in Chinese) |
[7] |
汪木兰, 张崇巍, 林健, 等. 数控机床用直线电机伺服试验平台开发及应用[J]. 中国机械工程, 2012, 23(3): 274-278. WANG M L, ZHANG C W, LIN J, et al. Development of an experimental platform for linear motors in CNC machine tools[J]. China Mechanical Engineering, 2012, 23(3): 274-278. DOI:10.3969/j.issn.1004-132X.2012.03.005 (in Chinese) |
[8] |
WANG H Y, WANG F J, ZHAO X Y, et al. Dynamic characteristics of a high speed precision positioning platform driven by linear motors[J]. Applied Mechanics and Materials, 2012, 197: 46-49. DOI:10.4028/www.scientific.net/AMM.197.46 |
[9] |
PORTMAN V T, CHAPSKY V S, SHNEOR Y. Evaluation and optimization of dynamic stiffness values of the PKMs: Collinear stiffness value approach[J]. Mechanism and Machine Theory, 2014, 74: 216-244. DOI:10.1016/j.mechmachtheory.2013.12.009 |
[10] |
MI L, YIN G F, SUN M N, et al. Effects of preloads on joints on dynamic stiffness of a whole machine tool structure[J]. Journal of Mechanical Science and Technology, 2012, 26(2): 495-508. DOI:10.1007/s12206-011-1033-4 |
[11] |
SHEN B H, TSAI M C. Robust dynamic stiffness design of linear servomotor drives[J]. Control Engineering Practice, 2006, 14(11): 1325-1336. DOI:10.1016/j.conengprac.2005.09.001 |
[12] |
彭波. 数控机床双摆铣头力矩电机高刚度控制策略研究[D]. 沈阳: 沈阳工业大学, 2011. PENG B. High stiffness control strategies of NC double swivel milling head driven by torque motors[D]. Shenyang: Shenyang University of Technology, 2011. (in Chinese) |
[13] |
尹明, 范丽鹏. 基于Bode图的数控转台伺服动刚度分析方法[J]. 制造技术与机床, 2013(1): 97-99. YIN M, FAN L P. Dynamic servo stiffness analysis method of NC rotary table based on Bode diagram[J]. Manufacturing Technology and Machine Tool, 2013(1): 97-99. DOI:10.3969/j.issn.1005-2402.2013.01.032 (in Chinese) |
[14] |
BAI Y H, QUAN L. New method to improve dynamic stiffness of electro-hydraulic servo systems[J]. Chinese Journal of Mechanical Engineering, 2013, 26(5): 997-1005. |
[15] |
莫艽. 基于S试件的加工原理误差和机床动态误差研究[D]. 北京: 清华大学, 2018. MO J. Research on the machining theoretical error and dynamic error of machine tool base on the S-shaped test piece[D]. Beijing: Tsinghua University, 2018. (in Chinese) |
[16] |
马平, 陈振环, 李劼科, 等. 零传动机床的高速直线进给单元的伺服动刚度研究[J]. 中国机械工程, 2004, 15(7): 575-577, 581. MA P, CHEN Z H, LI J K, et al. Dynamic stiffness of high velocity feed unit of the zero transmission machine tools[J]. China Mechanical Engineering, 2004, 15(7): 575-577, 581. DOI:10.3321/j.issn:1004-132X.2004.07.004 (in Chinese) |