技术特征:1.一种五轴数控制齿机床几何误差实际逆向运动学补偿方法,其特征在于,包括以下步骤:
(1)设定机床局部坐标系,确定各运动轴之间的相互耦合关系和运动传递链,建立机床运动学模型;
(2)将加工轨迹刀位数据表示为(x,y,z,i,j,k),其中x、y、z指加工轨迹的位置坐标,i、j、k指加工轨迹的方向坐标,由逆向运动学算法求得对应X轴、Y轴、Z轴、A轴和C轴的理想加工代码X、Y、Z、A、C,并换算到机床坐标系下;
(3)将辨识得到的机床位置相关误差与机床坐标系下运动轴位置X、Y、Z、A、C拟合得到对应的机床位置相关误差与运动轴位置之间的函数关系δXx(X)、δXy(X)、δXz(X)、αX(X)、βX(X)、γX(X),δYx(Y)、δYy(Y)、δYz(Y)、αY(Y)、βY(Y)、γY(Y),δZx(Z)、δZy(Z)、δZz(Z)、αZ(Z)、βZ(Z)、γZ(Z),δAx(A)、δAy(A)、δAz(A)、αA(A)、βA(A)、γA(A),δCx(C)、δCy(C)、δCz(C)、αC(C)、βC(C)、γC(C),由X、Y、Z、A、C预测待补偿的机床位置相关误差;
(4)建立齿轮/刀具24项位姿偏差与齿轮误差的耦合映射关系并进行解耦,以工件坐标系为基准,将齿轮/刀具的24项位姿偏差转化为刀具相对于工件的12项位姿偏差(x/y/z/α/β/γ方向6项静态位姿偏差和6项运动位姿偏差):Δx,Δx(p),Δy,Δy(p),Δz,Δz(p),Δi,Δi(p),Δj,Δj(p),Δk,Δk(p),p指机床运动轴位置X、Y、Z、A、C;根据理想刀位数据与刀具相对于齿轮的位姿偏差,求得修正后的刀位数据(x’,y’,z’,i’,j’,k’)=(x-Δx-Δx(p),y-Δy-Δy(p),z-Δz-Δz(p),i-Δi-Δi(p),j-Δj-Δj(p),k-Δk-Δk(p));
(5)建立机床空间几何误差模型,确定含误差的实际刀具坐标系与实际工件坐标系之间的变换矩阵(RBX·XBC·CBW)-1RBY·YBZ·ZBA·ABT (1);
(6)实际逆向运动学算法具体是应用齐次坐标变换矩阵的可逆特性、乘法结合律、矩阵分块计算原理、运动变换不变特征等将误差模型中机床运动代码X、Y、Z、A、C与机床几何误差进行解耦,得到修正后的机床运动代码解析表达式;
(7)将修正后的刀位数据(x’、y’、z’、i’、j’、k’)、机床几何误差δXx(X)、δXy(X)、δXz(X)、αX(X)、βX(X)、γX(X),δYx(Y)、δYy(Y)、δYz(Y)、αY(Y)、βY(Y)、γY(Y),δZx(Z)、δZy(Z)、δZz(Z)、αZ(Z)、βZ(Z)、γZ(Z),δAx(A)、δAy(A)、δAz(A)、αA(A)、βA(A)、γA(A),δCx(C)、δCy(C)、δCz(C)、αC(C)、βC(C)、γC(C),δCx、δCy、αCY、βCX,δAy、δAz、γAY、βAZ,γYX,αZY、βZX代入上述推导得到的解析表达式中,可求得补偿后的运动代码,实现机床几何误差的补偿,对刀具/齿轮位姿偏差进行修正。
2.根据权利要求1所述的一种五轴数控制齿机床几何误差实际逆向运动学补偿方法,其特征在于所述步骤(6)中求得补偿后的运动代码与修正的刀位数据、机床几何误差项的解析函数关系式为:
3.根据权利要求2所述的一种五轴数控制齿机床几何误差实际逆向运动学补偿方法,其特征在于所述步骤(6)中考虑旋转轴运动与直线轴运动之间的耦合关系,首先对旋转轴运动代码解析表达式进行推导得到:
式中:
再对直线轴运动代码解析表达式进行推导得到:
Z'=z'+δWz-δAz-L+δCz(C)+δXz(X)-δYz(Y)-δZz(Z)+(αW+αC(C))y'-(βW+βC(C))x'+
(αCY+αX(X))y'cos(C')-(βCX+βX(X))x'cos(C')+(αCY+αX(X))x'sin(C')+
(βCX+βX(X))y'sin(C')-(δTy-δAy+δAy(A)+(αA(A)+αY(Y)+αZY+αZ(Z))L)sin(A')+
(L+δAz-δTz-δAz(A))cos(A')
Y'=(δTz-δAz-L+δAz(A))sin(A')+δCy-δAy+δXy(X)-δYy(Y)-δZy(Z)+(αY(Y)+αZY+αZ(Z))L-(δTy-δAy+δAy(A)+(αA(A)+αY(Y)+αZY+αZ(Z))L)cos(A')-(αCY+αX(X))z'+
(δWy+δCy(C)+y'-(αW+αC(C))z'+(γW+γC(C)+γX(X))x')cos(C')+
(δWx+δCx(C)+x'+(βW+βC(C))z'-(γW+γC(C)+γX(X))y')sin(C')+(αY(Y)+αZY)Z'
X'=-{(βAZ+βY(Y)+βZX+βZ(Z))Lcos(A')+δCx-δTx-δAx(A)+δXx(X)-δYx(Y)-δZx(Z)+
(βA(A)-βAZ-βY(Y)-βZX-βZ(Z))L+(βCX+βX(X))z'+
(δWx+δCx(C)+x'+(βW+βC(C))z'-(γW+γC(C)+γX(X))y')cos(C')-
(δWy+δCy(C)+y'-(αW+αC(C))z'+(γW+γC(C)+γX(X))x')sin(C')+
(γAY+γYX+γY(Y)+γZ(Z))Lsin(A')}+(βY(Y)+βZX)Z'-γYXY'。
4.根据权利要求1所述的一种五轴数控制齿机床几何误差实际逆向运动学补偿方法,其特征在于所述步骤(5)中指C轴实际坐标系到工件实际坐标系的运动变换;
指含C轴位置无关误差δCx、δCy、αCY、βCX及位置相关误差δCx(C)、δCy(C)、δCz(C)、αC(C)、βC(C)、γC(C)的X轴实际坐标系到实际C轴实际坐标系的运动变换;
指含X轴位置相关误差δXx(X)、δXy(X)、δXz(X)、αX(X)、βX(X)、γX(X)的参考坐标系到X轴实际坐标系的运动变换;
指含Y轴、X轴垂直度误差γYX及Y轴位置相关误差δYx(Y)、δYy(Y)、δYz(Y)、αY(Y)、βY(Y)、γY(Y)的参考坐标系到Y轴实际坐标系的运动变换;
指含Z轴、Y轴,Z轴、X轴垂直度误差αZY、βZX以及Z轴位置相关误差δZx(Z)、δZy(Z)、δZz(Z)、αZ(Z)、βZ(Z)、γZ(Z)的Y轴实际坐标系到Z轴实际坐标系的运动变换;
指的是含A轴位置无关误差δAy、δAz、γAY、βAZ及位置相关误差δAx(A)、δAy(A)、δAz(A)、αA(A)、βA(A)、γA(A)的Z轴实际坐标系到A轴实际坐标系的运动变换;
指A轴实际坐标系到刀具实际坐标系的运动变换。
5.根据权利要求1所述的一种五轴数控制齿机床几何误差实际逆向运动学补偿方法,其特征在于所述步骤(6)中最终推导得到误差补偿的运动解析表达式,具体步骤如下:
假设修正后的刀位数据为(x,、y,、z,、i,、j,、k’),几何误差补偿后的机床运动代码为(X’,Y‘,Z’,A’,C’),则根据式(1)可建立如下关系式:
[i';j';k';0]=(RBX(X')·XBC(C')·CBW)-1RBY(Y')·YBZ(Z')·ZBA(A')·ABT·[0;0;1;0] (2)
[x';y';z';1]=(RBX(X')·XBC(C')·CBW)-1RBY(Y')·YBZ(Z')·ZBA(A')·ABT·[0;0;0;1] (3)
其中[0;0;1;0]指实际刀具坐标系下的刀轴矢量齐次坐标,[0;0;0;1]指实际刀具坐标系下的刀具中心点位置齐次坐标;
考虑到所有平动变换不会对矢量方向产生影响,在求解旋转轴运动代码解析表达式时,将所有平动值设为零,同时将式(2)推导转换到如式(4)、式(5)所示的方程形式,我们发现方程左右两侧计算结果均为一个4×1的实数矩阵,同时观察到式(4)左侧第3行元素是不含C轴运动的,式(5)右侧第一行元素是不含A轴运动的,这是由于仅有C轴旋转运动不会对变换前后的Z轴坐标有任何影响,同样仅有A轴旋转运动,不会对变换前后的X轴坐标有任何影响(这一原理同样适用于B轴的旋转运动),可根据这一特征,将特征元素A’,C‘分离开来,得到仅关于A’或者仅关于C‘的超越方程,从而可求得其解析表达式:
直线轴代码解析表达式的求解主要运用齐次变换矩阵分块计算原理推导得到:
式中R3×1,3×1均指3×1实数矩阵,指平动齐次变换矩阵,I3×3指3×3单位矩阵,利用上式计算原理、矩阵乘法结合律以及齐次变换矩阵均可逆的性质,可将式(3)推导为如下形式:
由式(6)将直线轴运动代码分离开来,进而得到关于直线轴运动代码的简单线性方程组;忽略高阶项,则由式(4)、(5)、(6)推导得到的运动轴代码解析表达式为:
式中: