辅助资料2C数控机床闭环进给伺服系统运动误差的研究

1补充内容：数控机床闭环进给伺服系统轮廓误差分析1．闭环位置控制的数学模型典型X轴位置闭环控制简化数学模型如图1所示。图1简化位置闭环控制数学模型其中，Kp为控制器内部的软件位置增益，可用于调整系统位置环的开环增益；Kda（V/P）为数模转换系数；Ke（P/R）为位置编码器的脉冲数；Kn（V/RPM）为速度反馈系数；Ks为速度调节器增益；Kb（V/RPM）为反电势常数；Ra（Ω）为电枢回路电阻；Kq（Nm/A）为电机力矩常数；J（Kg•m2）为转动惯量。速度闭环控制模型可简化为：1)(sKsWm其中Km为每伏电压对应电机转速nsbsmKKKKKτ为时间常数)(nsbqaKKKKJR位置闭环控制的开环增益emdapKKKKK。当不考虑速度闭环控制模型时间常数τ时，位置控制开环传递函数为sKGk/，其闭环传递函数为如下一阶模型：11)(TssG(1)其中KT/1当考虑速度闭环控制模型时间常数τ时，位置控制开环传递函数为)]1(/[ssKGk，其闭环传递函数为二阶模型：2222)(nnnsssG(2)2其中K121Kn在数控系统中X、Y轴通常具有同样的上述位置闭环控制数学模型，因此由X、Y两轴独立位置闭环控制组成的系统如图2所示。图2两轴闭环位置系统其中Xi、X分别是X轴的输入输出，Yi、Y分别是Y轴的输入输出，下面模型参数以下标x或y区分其为X或Y轴参数。2．直线插补运动的轮廓误差图1所示闭环系统的误差传递函数：)(11)()(1)()()(sGsXsXsXsEsGkioie(3)直线插补时，根据拉氏变换终值定理)(11lim)(lim)(00sGsssEekss(4)当进行XY轴直线联动插补时，对应X、Y轴的指令为斜坡输入，即tVtXxi)(tVtYyi)(由式4，对一阶模型和二阶模型，其稳态误差均为xxxKVE/yyyKVE/设运动直线如图3所示，与X轴的夹角为θ，合成进给速度为V，则有运动轮廓误差为VKVVVKVVVKVVVKVVEEExyxyxyxyxycossinsincossincosyxyxKKKKV2)2sin((5)图3直线运动的轮廓误差3由式5可知（1）当yxKK，即两轴增益相等（既nyynxx//）时，两轴跟随误差的滞后效应抵消，稳态轮廓误差为零，Jx、Jy的不匹配对稳态轮廓误差没有影响；（2）当sin(2θ)=0，即θ=0或90º时，则E=0，这具有明显的物理意义；（3）实际系统中很难保证xK、yK完全相等，可以看出，只要xK、yK足够大，轮廓误差就会很小，因此尽可能选择较高的增益或使两轴增益尽可能一致均可减小轮廓误差；（4）轮廓误差与运动速度成正比；3．圆弧插补运动的轮廓误差圆弧插补运动的轮廓误差主要包括伺服系统有限带宽引起的圆弧半径误差和运动轴性能不匹配引起的椭圆误差。考虑沿半径为R的圆弧进行角速度为ω的运动，圆心坐标为（0，0），X、Y轴的指令输入及初始条件为)cos()(tRtxi且Rxi)0(，0)0(ix；)sin()(tRtyi且0)0(iy，Ryi)0(。稳态输出为]sinsincos[cos)cos()(xxxxxttRtRtx(6)]sincoscos[sin)sin()(yyyyyttRtRty(7)这里xR、yR为输出幅值，x、y为滞后相位。由式6、式7得)](sin)(cos[)cos(1costyRtxRRRtxxyyyxyx(8))](cos)(sin[)cos(1sintyRtxRRRtxxyyyxyx(9)且输出轨迹为下列方程)(cos)sin(22222222yxyxyxyxxyRRxyRRyRxR(10)对一般二次曲线方程有022233231322212211ayaxayaxyaxa(11)将坐标系oxy绕原点o旋转角度θ到新坐标系o'x'y'，则有cos'sin'sin'cos'yxyyxx(12)将式2.12代入2.11得0'''2''2'''''2''33231322212211ayaxayayxaxa(13)其中：43333231323231313222122112222121122122221221111'cossin'sincos'coscossin2sin')sin(coscossin)('sinsincos2cos'aaaaaaaaaaaaaaaaaaaa由式13可以看出如下规律：（1）式13方程二次项系数仅与原方程二次项系数及旋转角θ有关，与一次项系数和常数无关。（2）式13方程一次项系数仅与原方程一次项系数及旋转角θ有关，与二次项系数和常数无关。（3）当原方程一次项系数不为零时，通过坐标轴旋转不能完全消除一次项系数；同样当原方程一次项系数为零时，坐标轴旋转后一次项仍为零。（4）坐标轴旋转后方程的常数项不变。（5）当012a时，通过坐标轴旋转可消除交叉乘积项。既当122211221aaaarcctg时)sin(coscossin)('2212112212aaaa02cos2sin)(21121122aaa(14)对照式10、式11，当坐标系旋转角度)sin(22122yxyxxyRRRRarcctg(15)式10变换为标准的椭圆方程。由此可见，一般地说圆弧实际运动的轨迹为以原点为中心，长轴或短轴在角度θ处的椭园。3．1一阶模型圆弧运动误差分析当图2中Gx(s)、Gy(s)为式1所示一阶模型时，沿半径为R的圆弧进行角速度为ω的运动，其稳态输出为)cos()(1)(2xxtTRtxxxTtan(16))cos()(1)(2yytTRtyyyTtan(17)当TTTyx时yx，且2)(1'TRRRRyx由式10可知此时轨迹方程为222'Ryx圆的半径误差为5RTFTRTRRRR2)(21])(111['2222(18)可见，当两轴匹配时，实际运动轨迹为圆，其半径总是比指令圆弧半径值小，半径误差随时间常数和进给速度的增大而增大。显然，这是由伺服系统有限带宽引起的圆弧半径误差。3．2二阶模型圆弧运动误差分析当图2中Gx(s)、Gy(s)为式2所示二阶模型时，沿半径为R的圆弧进行角速度为ω的运动，其稳态输出为)cos()cos()2()()(22222xxxnxxnxnxtRtRtx其中222tannxnxxx)cos()cos()2()()(22222yyynyynynytRtRty其中222tannynyyy当yx，nnynx时，则yx由式10可知其运动轨迹为圆，其轨迹方程为222'Ryx42222222)()(2)2(1)2()('nnnnnnyxRRRRR(19)可见，当两轴匹配时，实际稳态运动轨迹为圆，但显然圆的半径存在误差。图4为两轴匹配且145n时，对不同的ξ值，RR/'随运动角速度ω的变化曲线6图4两轴匹配时RR/'随ω的变化曲线当707.0时，二阶系统的谐振频率ωr和谐振峰值Mr分别为221nr2121rM因此，圆的轨迹误差是由于幅频特性不等于1引起的。因此当阻尼比较小时，通常情况下圆的运动半径大于指令半径。由于当707.0时，系统不产生谐振。显然，此时圆的半径小于指令半径。通过调整系统增益使得707.0，从而在不产生谐振的条件下，得到更宽的频带。从式19可以看出，当707.0时4)(1'nyxRRRR(20)当n时，由于4次幂的作用，实际产生的误差较小。因此使两轴匹配且707.0可取得较小的圆弧运动误差。但在大多数的应用中，很难保证两轴完全匹配。例如在立式数控铣床和加工中心中，X轴的惯性负载通常远大于Y轴，即yxJJ。在此情况下，通过调整位置增益满足两轴匹配条件需要同时满足nynx，yx即yxyxJJKK(21)xyyxJJKK(22)显然式21与22不可能同时成立，而使两轴匹配条件完全满足。由上可见，虽然通过精心调整，常规位置控制算法可以满足通常数控机床轮廓加工的需要。但随着现代机械加工对加工精度和加工效率要求的不断提高，常规位置闭环控制算法已不相适应，因此对新的控制算法进行研究就显得非常必要。4．仿真研究4．1仿真数学模型图5为数控机床所采用的典型双环(电流环、速度环)进给伺服单元的传递函数数学模型。其中In1为负载干扰力矩，In2为数控装置位置闭环控制电压，Out1为转速输出。7图5典型双闭环伺服单元数学模型取Kx、Ky为软件增益，由16位DAC所产生的转换系数为0.0001526，Kxe、Kye为位置传感器（8mm/r）的检测系数，ScopeError为轮廓误差值。4．2伺服系统及传动机构对圆弧运动精度的影响分析加工整圆（标准圆）并检查其误差分布是数控机床精度综合检验最常用的方法。设顺时针加工半径为100mm的整圆，运动角频率为，)sin(100tXi，)cos(100tYi。下面讨论位置开环增益对圆弧运动精度的影响。(1)当62yxKK时，随角频率ω的变化轮廓误差的变化如表1，所对应的关系曲线如图6b。图6a为05.0时轮廓误差ScopeError(t)的变化曲线。表1：角频率ω与轮廓误差的关系ω(1/s)0.050.10.150.20.250.30.350.5ScopeError(μm)-0.7-2.84-6.37-11-17.5-25.5-34.7-70.5a.轮廓误差的变化曲线b.随角频率的变化轮廓误差的变化图6随角频率ω的变化轮廓误差变化曲线(2)当yxKK，且1.0时，轮廓运动误差ScopeError(t)如图7；各轴增益与轮廓误差的关系如表2；图8为当62xK时对于不同的yK，轮廓误差最大值与最小值的变化曲线。8表2：yxKK时各轴增益与轮廓误差的关系Kx626262626262626262Ky585960616263646566ScopeErrorMax2316103-2.8391520ScopeErrorMin-29-22-15-9-2.8-9-14-20-25a.yxKK时轮廓误差的变化b.yxKK时轮廓误差的变化图7yxKK时轮廓误差的变化曲线图8随yK的变化轮廓误差最大值与最小值的变化曲线根据上述曲线分析，可得出如下结论：——当yxKK时，圆的半径误差为一恒定值，且随着角频率的增加而增加；——当yxKK时，圆的半径误差明显增加且呈正弦规律变化，即实际运动轨迹为一椭圆；9——当yxKK时，正最大误差处在椭圆的长轴，即大约在45°、225°位置上；——当yxKK时，正最大误差处在椭圆的长轴，即大约在135°、315°位置上。可见，上述仿真结果与前述圆弧插补运动轮廓误差的分析是吻合的。5．结论目前数控机床进给运动大多采用各轴独立的跟随控制，本文对此类常规跟随控制算法由于有限带宽和各轴特性不匹配所引起的轮廓误差特性进行了理论分析，并讨论了其内在的规律。事实上，还有许多因素均会对进给运动精度产生重要影响，如导轨的非线性摩擦特性等。因此有必要进一步研究各种先进的控制与补偿算法，以提高伺服运动轴的动态性能，从而达到改善轮廓运动精度的目的。此外，从本质上看，跟随控制仍然是轮廓开环控制系统；通过对轮廓误差的计算或估计，向各轴提供附加轮廓误差补偿信息，从而实现系统的轮廓闭环控制也是重要的研究方向。参考文献[1]GeneF.Frankl