1.一種五軸數(shù)控制齒機(jī)床幾何誤差實(shí)際逆向運(yùn)動(dòng)學(xué)補(bǔ)償方法,其特征在于,包括以下步驟:
(1)設(shè)定機(jī)床局部坐標(biāo)系,確定各運(yùn)動(dòng)軸之間的相互耦合關(guān)系和運(yùn)動(dòng)傳遞鏈,建立機(jī)床運(yùn)動(dòng)學(xué)模型;
(2)將加工軌跡刀位數(shù)據(jù)表示為(x,y,z,i,j,k),其中x、y、z指加工軌跡的位置坐標(biāo),i、j、k指加工軌跡的方向坐標(biāo),由逆向運(yùn)動(dòng)學(xué)算法求得對(duì)應(yīng)X軸、Y軸、Z軸、A軸和C軸的理想加工代碼X、Y、Z、A、C,并換算到機(jī)床坐標(biāo)系下;
(3)將辨識(shí)得到的機(jī)床位置相關(guān)誤差與機(jī)床坐標(biāo)系下運(yùn)動(dòng)軸位置X、Y、Z、A、C擬合得到對(duì)應(yīng)的機(jī)床位置相關(guān)誤差與運(yùn)動(dòng)軸位置之間的函數(shù)關(guān)系δ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預(yù)測(cè)待補(bǔ)償?shù)臋C(jī)床位置相關(guān)誤差;
(4)建立齒輪/刀具24項(xiàng)位姿偏差與齒輪誤差的耦合映射關(guān)系并進(jìn)行解耦,以工件坐標(biāo)系為基準(zhǔn),將齒輪/刀具的24項(xiàng)位姿偏差轉(zhuǎn)化為刀具相對(duì)于工件的12項(xiàng)位姿偏差(x/y/z/α/β/γ方向6項(xiàng)靜態(tài)位姿偏差和6項(xiàng)運(yùn)動(dòng)位姿偏差):Δx,Δx(p),Δy,Δy(p),Δz,Δz(p),Δi,Δi(p),Δj,Δj(p),Δk,Δk(p),p指機(jī)床運(yùn)動(dòng)軸位置X、Y、Z、A、C;根據(jù)理想刀位數(shù)據(jù)與刀具相對(duì)于齒輪的位姿偏差,求得修正后的刀位數(shù)據(jù)(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)建立機(jī)床空間幾何誤差模型,確定含誤差的實(shí)際刀具坐標(biāo)系與實(shí)際工件坐標(biāo)系之間的變換矩陣(RBX·XBC·CBW)-1RBY·YBZ·ZBA·ABT (1);
(6)實(shí)際逆向運(yùn)動(dòng)學(xué)算法具體是應(yīng)用齊次坐標(biāo)變換矩陣的可逆特性、乘法結(jié)合律、矩陣分塊計(jì)算原理、運(yùn)動(dòng)變換不變特征等將誤差模型中機(jī)床運(yùn)動(dòng)代碼X、Y、Z、A、C與機(jī)床幾何誤差進(jìn)行解耦,得到修正后的機(jī)床運(yùn)動(dòng)代碼解析表達(dá)式;
(7)將修正后的刀位數(shù)據(jù)(x’、y’、z’、i’、j’、k’)、機(jī)床幾何誤差δ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代入上述推導(dǎo)得到的解析表達(dá)式中,可求得補(bǔ)償后的運(yùn)動(dòng)代碼,實(shí)現(xiàn)機(jī)床幾何誤差的補(bǔ)償,對(duì)刀具/齒輪位姿偏差進(jìn)行修正。
2.根據(jù)權(quán)利要求1所述的一種五軸數(shù)控制齒機(jī)床幾何誤差實(shí)際逆向運(yùn)動(dòng)學(xué)補(bǔ)償方法,其特征在于所述步驟(6)中求得補(bǔ)償后的運(yùn)動(dòng)代碼與修正的刀位數(shù)據(jù)、機(jī)床幾何誤差項(xiàng)的解析函數(shù)關(guān)系式為:
3.根據(jù)權(quán)利要求2所述的一種五軸數(shù)控制齒機(jī)床幾何誤差實(shí)際逆向運(yùn)動(dòng)學(xué)補(bǔ)償方法,其特征在于所述步驟(6)中考慮旋轉(zhuǎn)軸運(yùn)動(dòng)與直線軸運(yùn)動(dòng)之間的耦合關(guān)系,首先對(duì)旋轉(zhuǎn)軸運(yùn)動(dòng)代碼解析表達(dá)式進(jìn)行推導(dǎo)得到:
式中:
再對(duì)直線軸運(yùn)動(dòng)代碼解析表達(dá)式進(jìn)行推導(dǎo)得到:
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.根據(jù)權(quán)利要求1所述的一種五軸數(shù)控制齒機(jī)床幾何誤差實(shí)際逆向運(yùn)動(dòng)學(xué)補(bǔ)償方法,其特征在于所述步驟(5)中指C軸實(shí)際坐標(biāo)系到工件實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換;
指含C軸位置無(wú)關(guān)誤差δCx、δCy、αCY、βCX及位置相關(guān)誤差δCx(C)、δCy(C)、δCz(C)、αC(C)、βC(C)、γC(C)的X軸實(shí)際坐標(biāo)系到實(shí)際C軸實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換;
指含X軸位置相關(guān)誤差δXx(X)、δXy(X)、δXz(X)、αX(X)、βX(X)、γX(X)的參考坐標(biāo)系到X軸實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換;
指含Y軸、X軸垂直度誤差γYX及Y軸位置相關(guān)誤差δYx(Y)、δYy(Y)、δYz(Y)、αY(Y)、βY(Y)、γY(Y)的參考坐標(biāo)系到Y(jié)軸實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換;
指含Z軸、Y軸,Z軸、X軸垂直度誤差αZY、βZX以及Z軸位置相關(guān)誤差δZx(Z)、δZy(Z)、δZz(Z)、αZ(Z)、βZ(Z)、γZ(Z)的Y軸實(shí)際坐標(biāo)系到Z軸實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換;
指的是含A軸位置無(wú)關(guān)誤差δAy、δAz、γAY、βAZ及位置相關(guān)誤差δAx(A)、δAy(A)、δAz(A)、αA(A)、βA(A)、γA(A)的Z軸實(shí)際坐標(biāo)系到A軸實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換;
指A軸實(shí)際坐標(biāo)系到刀具實(shí)際坐標(biāo)系的運(yùn)動(dòng)變換。
5.根據(jù)權(quán)利要求1所述的一種五軸數(shù)控制齒機(jī)床幾何誤差實(shí)際逆向運(yùn)動(dòng)學(xué)補(bǔ)償方法,其特征在于所述步驟(6)中最終推導(dǎo)得到誤差補(bǔ)償?shù)倪\(yùn)動(dòng)解析表達(dá)式,具體步驟如下:
假設(shè)修正后的刀位數(shù)據(jù)為(x,、y,、z,、i,、j,、k’),幾何誤差補(bǔ)償后的機(jī)床運(yùn)動(dòng)代碼為(X’,Y‘,Z’,A’,C’),則根據(jù)式(1)可建立如下關(guān)系式:
[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]指實(shí)際刀具坐標(biāo)系下的刀軸矢量齊次坐標(biāo),[0;0;0;1]指實(shí)際刀具坐標(biāo)系下的刀具中心點(diǎn)位置齊次坐標(biāo);
考慮到所有平動(dòng)變換不會(huì)對(duì)矢量方向產(chǎn)生影響,在求解旋轉(zhuǎn)軸運(yùn)動(dòng)代碼解析表達(dá)式時(shí),將所有平動(dòng)值設(shè)為零,同時(shí)將式(2)推導(dǎo)轉(zhuǎn)換到如式(4)、式(5)所示的方程形式,我們發(fā)現(xiàn)方程左右兩側(cè)計(jì)算結(jié)果均為一個(gè)4×1的實(shí)數(shù)矩陣,同時(shí)觀察到式(4)左側(cè)第3行元素是不含C軸運(yùn)動(dòng)的,式(5)右側(cè)第一行元素是不含A軸運(yùn)動(dòng)的,這是由于僅有C軸旋轉(zhuǎn)運(yùn)動(dòng)不會(huì)對(duì)變換前后的Z軸坐標(biāo)有任何影響,同樣僅有A軸旋轉(zhuǎn)運(yùn)動(dòng),不會(huì)對(duì)變換前后的X軸坐標(biāo)有任何影響(這一原理同樣適用于B軸的旋轉(zhuǎn)運(yùn)動(dòng)),可根據(jù)這一特征,將特征元素A’,C‘分離開(kāi)來(lái),得到僅關(guān)于A’或者僅關(guān)于C‘的超越方程,從而可求得其解析表達(dá)式:
直線軸代碼解析表達(dá)式的求解主要運(yùn)用齊次變換矩陣分塊計(jì)算原理推導(dǎo)得到:
式中R3×1,3×1均指3×1實(shí)數(shù)矩陣,指平動(dòng)齊次變換矩陣,I3×3指3×3單位矩陣,利用上式計(jì)算原理、矩陣乘法結(jié)合律以及齊次變換矩陣均可逆的性質(zhì),可將式(3)推導(dǎo)為如下形式:
由式(6)將直線軸運(yùn)動(dòng)代碼分離開(kāi)來(lái),進(jìn)而得到關(guān)于直線軸運(yùn)動(dòng)代碼的簡(jiǎn)單線性方程組;忽略高階項(xiàng),則由式(4)、(5)、(6)推導(dǎo)得到的運(yùn)動(dòng)軸代碼解析表達(dá)式為:
式中: