一种五轴加工中球形刀刀轴矢量规划方法

摘要

本发明提出一种五轴加工中球形刀刀轴矢量规划方法。该方法采用等距螺旋走刀运动，首先将复杂曲面网格化后提取刀触点，并确定该点所在切平面。然后通过分类计算所在切平面法向矢量，在空间中经过罗德里格旋转变化后转化为可用于实际加工的球形刀刀轴矢量。最后分析加工中相邻刀轴矢量旋转运动带来的非线性误差，利用四元数球面线性插值法控制旋转误差角的范围，保证复杂曲面零件加工质量。本发明刀轴矢量规划方法操作简单、易懂，规划的刀轴矢量变化均匀，加工轨迹平滑、连续，整体非线性误差角满足精度要求，应用范围较广，适用于各种类型的五轴数控机床。

权利要求

1.一种五轴加工中球形刀刀轴矢量规划方法，其特征在于，包括如下步骤：(1)将复杂曲面整体分割成三角网格，三角网格化的弦高误差δ和角度误差β由该复杂曲面的加工精度决定，通常设定0.02mm≤δ≤0.04mm，0.5°≤β≤2.5°，确定每一个三角网格Tj的三个顶点Aj、Bj、Cj的位置坐标，及三条边AjBj、BjCj、AjCj所代表的线性方程；

(2)采用等距螺旋线的走刀方式确定所有刀触点的空间坐标Pi(xi,yi,zi)，其中xi、yi、zi为坐标Pi在X、Y、Z轴上的坐标，相邻刀触点Pi之间距离由球形刀刀具半径和加工精度确定；

(3)通过计算刀触点Pi与三角网格三个顶点Aj、Bj、Cj间向量的叉积，来确定刀触点Pi所属的三角形网格坐标及其在该三角形网格中的位置，并根据Pi在三角形中位置分为三类：第一类：刀触点Pi在三角网格内；第二类：刀触点Pi在三角网格顶点上；第三类：刀触点Pi在三角网格边上；

(4)根据刀触点Pi所属分类计算该点法向矢量，第第二类、第三类类情况需考虑到刀触点Pi周围三角网格法向矢量对其影响，引入形状权重因子进行补偿，精确得出每个刀触点Pi在复杂曲面中的法向矢量ni(nxi,nyi,nzi)，其中nxi、nyi、nzi为法向矢量ni在X、Y、Z轴上的坐标；

(5)将刀触点Pi的法向矢量ni绕空间中某一单位向量Fi(ai,bi,ci)，其中ai、bi、ci为单位向量Fi在X、Y、Z轴上的坐标，进行罗德里格旋转角度 后得到适用于球形刀加工的刀轴矢量Nki(Nxi,Nyi,Nzi),其中Nxi、Nyi、Nzi为刀轴矢量Nki在X、Y、Z轴上的坐标；

(6)确定五轴机床各运动轴之间的耦合关系和运动传递链，建立机床运动学模型，并根据工件坐标系和机床坐标系之间的转换，确定相邻刀轴矢量Nki、Nki+1之间由加工运动所产生的最大非线性误差角αi，及αi与相邻刀轴矢量Nki、Nki+1间夹角θi关系；

(7)利用四元数球面线性插补法缩小相邻刀轴矢量间夹角θi，从而控制非线性误差角αi，并确保该复杂曲面的加工精度；

(8)将所有满足加工精度要求的刀触点P和刀轴矢量Nk转换为五轴机床加工的各轴代码{X,Y,Z,A,C}。

2.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于所述步骤(3)中刀触点Pi与任意三角网格Tj的三个顶点Aj、Bj、Cj所构成向量计算三个向量之间的叉积如下：

3.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于所述步骤(4)中第二类中刀触点Pi的法向矢量计算如下：其中Wj、Aj和Nj分别为周围三角网格Tj的形状权重因子、面积和法向矢量。

4.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于所述步骤(4)中第三类中Pi点的法向矢量ni计算如下：ni＝(1-u)Nm+uNn+u(1-u)Ncomp，u∈[0,1]其中，Nm、Nn为法向矢量ni所在边PmPn的两个顶点Pm、Pn的法向矢量，Na、Nb为相邻两个网格Ta、Tb的法向矢量，Ncomp为Ta、Tb的补偿矢量，u为Pi点到Pn点的长度。

5.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于所述步骤(5)中法向矢量ni经过罗德里格旋转得到Nki的式如下：Νki＝Ri·ni

6.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于步骤(5)中单位向量Fi(ai,bi,ci)由刀触点Pi(xi,yi,zi)确定：

7.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于所述步骤(6)中Nki至Nki+1过程中非线性误差角αi表示如下：α(t)＝arccos(U(t))

其中，Nkp为Nki至Nki+1过程中某一位置理论刀轴矢量，Nkp’为该刀矢在实际加工中机床坐标系的表达，θi为相邻刀轴矢量Nki和Nki+1之间夹角，tθi为刀轴矢量Nkp与Nki的夹角，t∈[0,

1]。

8.根据权利要求1所述的一种五轴加工中球形刀刀轴矢量规划方法，其特征在于所述步骤(7)中四元数球面线性插补法确定Nki和Nki+1间新刀位点的刀轴矢量Nki+t：其中ti为相邻刀轴矢量Nki和Nki+1间非线性误差角αi最大位置处的t值，即t＝ti。

说明书

一种五轴加工中球形刀刀轴矢量规划方法

技术领域

[0001] 本发明属于五轴机床刀具轨迹规划领域，具体是一种五轴加工中球形刀轴矢量规划方法，是一 种将复杂曲面网格化后提取刀触点法向矢量，并通过空间旋转为球形刀刀轴矢量的规划方法。

背景技术

[0002] 五轴机床的加工精度和加工质量较于普通三轴铣床有很大提升，但两个旋转轴的增加使机床刀 具位置和运动姿态变得较复杂，容易在实际加工过程中产生奇异问题和干涉碰撞，甚至影响机床整 体性能。目前常用的刀具轨迹规划方法生成的刀轨与曲面模型本身有一定差距，由于曲面变化复杂 还会出现加工轨迹分段不断退刀的问题，从而降低加工效率。五轴加工中还要考虑机床本身的动力 学特点，避免运动轨迹突变和不连续造成不断退刀，对加工表面质量和机床本身性能造成影响。

[0003] 五轴机床刀轴矢量规划不合理还容易在实际加工中产生非线性误差，而造成加工表面误差过 大。现有的控制方法主要通过分析五轴机床加工时各轴的运动姿态，通过综合二维投影下的弦高误 差和旋转带来的误差，在相邻刀位点间直接插入新点来控制误差；或建立新的坐标系，将相邻间刀 轴矢量构成一个扇面进行线性插补来控制误差。以上方法非线性误差模型整体不够准确，对旋转轴 走刀方向垂直平面内非线性误差考虑较少，对于曲率变化较大的模型作用不明显，控制方法效果差。 因此亟需一种在空间上直接控制旋转误差角的球形刀刀轴矢量规划方法，并满足复杂曲面加工的精 度要求。

发明内容

[0004] 为了能够满足五轴加工中刀具轨迹规划在精度方面的要求，特别是旋转轴运动，本发明提出了 一种五轴加工中球形刀刀轴矢量规划方法。该方法利用三角网格的几何特点，将复杂曲面网格化后 提取刀触点和该点所属切平面的法向矢量，通过空间旋转后转化为球形刀刀轴矢量，并通过控制加 工中刀轴旋转的非线性误差，该方法简单易操作，满足复杂曲面加工的精度要求。

[0005] 一种五轴加工中球形刀刀轴矢量规划方法，包括如下步骤：

[0006] (1)将复杂曲面整体分割成三角网格。三角网格化的弦高误差δ和角度误差β由该复杂曲面的加工 精度决定，通常设定0.02mm≤δ≤0.04mm，0.5°≤β≤2.5°。确定每一个三角网格Tj的三个顶 点Aj、Bj、Cj的位置坐标，及三条边AjBj、BjCj、AjCj所代表的线性方程。

[0007] (2)采用等距螺旋线的走刀方式确定所有刀触点的空间坐标Pi(xi,yi,zi)，其中xi、yi、zi为坐标Pi在X、 Y、Z轴上的坐标，相邻刀触点Pi之间距离由球形刀刀具半径和加工精度确定。

[0008] (3)通过计算刀触点Pi与三角网格三个顶点Aj、Bj、Cj间向量的叉积，来确定刀触点Pi所属的三角 形网格坐标及其在该三角形网格中的位置，并根据Pi在三角形中位置分为三类：①刀触点Pi在三 角网格内；②刀触点Pi在三角网格顶点上；③刀触点Pi在三角网格边上。

[0009] (4)根据刀触点Pi所属分类计算该点法向矢量，第②、③类情况需考虑到刀触点Pi周围三角网格法 向矢量对其影响，引入形状权重因子进行补偿，精确得出每个刀触点Pi在复杂曲面中的法向矢量 ni(nxi,nyi,nzi),其中nxi、nyi、nzi为法向矢量ni在X、Y、Z轴上的坐标。

[0010] (5)将刀触点Pi的法向矢量ni绕空间中某一单位向量Fi(ai,bi,ci),其中ai、bi、ci为单位向量Fi在X、 Y、Z轴上的坐标，进行罗德里格旋转角度 后得到适用于球形刀加工的刀轴矢量Nki(Nxi,Nyi,Nzi),其 中Nxi、Nyi、Nzi为刀轴矢量Nki在X、Y、Z轴上的坐标。

[0011] (6)确定五轴机床各运动轴之间的耦合关系和运动传递链，建立机床运动学模型。并根据工件坐标 系和机床坐标系之间的转换，确定相邻刀轴矢量Nki、Nki+1之间由加工运动所产生的最大非线性误 差角αi，及αi与相邻刀轴矢量Nki、Nki+1间夹角θi关系。

[0012] (7)利用四元数球面线性插补法缩小相邻刀轴矢量间夹角θi，从而控制非线性误差角αi，并确保该 复杂曲面的加工精度。

[0013] (8)将所有满足加工精度要求的刀触点P和刀轴矢量Nk转换为五轴机床加工的各轴代码 {X,Y,Z,A,C}。

[0014] 所述步骤(3)中刀触点Pi与任意三角网格Tj的三个顶点Aj、Bj、Cj所构成向量计算三个向量之间的叉积如下：

[0015]

[0016]

[0017]

[0018] 所述步骤(4)中第②类中刀触点Pi的法向矢量计算如下：

[0019]

[0020] 其中Wj、Aj和Nj分别为周围三角网格Tj的形状权重因子、面积和法向矢量。

[0021] 所述步骤(4)中第③类中Pi点的法向矢量ni计算如下：

[0022] ni＝(1-u)Nm+uNn+u(1-u)Ncomp，u∈[0,1]

[0023]

[0024] 其中，Nm、Nn为法向矢量ni所在边PmPn的两个顶点Pm、Pn的法向矢量，Na、Nb为相邻两个 网格Ta、Tb的法向矢量。Ncomp为Ta、Tb的补偿矢量，u为Pi点到Pn点的长度。

[0025] 所述步骤(5)中法向矢量ni经过罗德里格旋转得到Nki的式如下：

[0026] Νki＝Ri·ni

[0027]

[0028] 步骤(5)中单位向量Fi(ai,bi,ci)由刀触点Pi(xi,yi,zi)确定：

[0029]

[0030]

[0031] 所述步骤(6)中Nki至Nki+1过程中非线性误差角αi表示如下：

[0032] α(t)＝arccos(U(t))

[0033]

[0034] 其中，Nkp为Nki至Nki+1过程中某一位置处理论刀轴矢量，Nkp’为该刀矢在实际加工中机床坐 标系的表达。θi为相邻刀轴矢量Nki和Nki+1之间夹角。tθi为刀轴矢量Nkp与Nki的夹角，t∈[0,1]。

[0035] 所述步骤(7)中四元数球面线性插补法确定Nki和Nki+1间新刀位点的刀轴矢量Nki+t：

[0036]

[0037] 其中ti为相邻刀轴矢量Nki和Nki+1间非线性误差角αi最大位置处的t值，即t＝ti。

[0038] 本发明的有益效果是：

[0039] 本发明通过将复杂曲面三角网格化的方法来计算曲面的法向矢量，简单易操作，且提取曲面的 法向矢量精度较高。利用罗德里格空间旋转矩阵，直接将法向矢量转化为球形刀具刀轴矢量，实施 性强，可直接用于实际加工。为了控制在实际加工中五轴机床运动产生的非线性误差，通过计算相 邻刀轴矢量间的非线性误差角αi及其位置，并在该位置插入新的刀轴矢量来缩小误差角。有效控 制旋转引起非线性加工误差，避免误差过大令刀轴矢量接近奇异点的矢量，使一个旋转轴运动速率 发生突变而产生奇异问题，保证曲面加工质量和加工精度。同时可运用于各构型五轴机床加工，适 用范围较广。

附图说明

[0040] 图1为本发明的五轴加工中球形刀刀轴矢量规划方法流程图。

[0041] 图2为本实施例的AC摆动/回转型五轴机床结构及运动链示意图。

[0042] 图3a为本实施例刀触点在三角网格位置分类情况之一示意图。

[0043] 图3b为本实施例刀触点在三角网格位置分类情况之二示意图。

[0044] 图3c为本实施例刀触点在三角网格位置分类情况之三示意图。

[0045] 图4为本发明法向矢量ni在空间中围绕单位向量Fi旋转角度 的示意图。

[0046] 图5为本实施例相邻刀轴矢量Nki和Nki+1之间运动示意图。

具体实施方式

[0047] 为使本发明的技术方案更加清楚，以下结合附图对本发明技术内容作详细说明。此处描述的实 例仅为解释本发明，并非对本发明的构思及应用范围作限定。

[0048] 如图1至图5，本发明通过将复杂曲面三角网格化，采用等距螺旋走刀方式排列并提取刀触点， 并根据刀触点所属三角网格中位置分类计算法向矢量，利用罗德里格旋转矩阵将法向矢量转化为刀 轴矢量，最后通过四元数球面线性插补法缩小相邻刀矢间夹角，控制加工中的非线性误差，保证加 工精度。

[0049] 如图2所示为AC式刀具与工作台分开回转型机床，本发明以该机床结构为例进行具体实施方 式叙述。

[0050] 为了计算出复杂曲面中刀触点Pi所在切平面的法向矢量ni，将复杂曲面整体碎片化，分割成 许多三角网格。三角网格的弦高误差δ和角度误差β由该复杂曲面的加工精度要求决定。弦高公差 δ指复杂曲面模型设计表面与三角网格表面的最大距离，弦高越小，三角表面越小，表面曲率也越 精确。角度误差β表示三角网格平面与曲面的切平面之间最大夹角，影响相邻三角网格法线之间角 度，减小角度将增加曲面网格化的分辨率，但在后期计算方面需要大量的时间。当三角网格足够小 时，复杂曲面中刀触点Pi的位置可近似于某一三角网格Tj中，则该点的切平面为所属三角网格所 在平面。本案例设定弦高误差δ为0.04mm，角度误差β为2.5°。根据三角网格Tj的顶点Aj、Bj、Cj的坐标位置，计算出垂直于该三角网格的法向矢量ni。

[0051] 结合零件加工表面质量和加工效率等方面要求，本发明采用等距螺旋线的走刀方式。在复杂曲 面表面选取点进行曲线拟合后，按等弦长进行分段取刀触点，使点之间步长和步距小于刀具半径， 保证加工精度。最后生成等距螺旋刀触点轨迹，得到刀触点集P。

[0052] 确定复杂曲面三角网格中所有刀触点位置后，计算每一个刀触点的法向矢量，用于后续球形刀 刀轴矢量的计算。依据刀触点Pi在三角网格Tj中位置，如图3主要分为三类：①刀触点Pi在三角 网格Tj内；②刀触点Pi在三角网格Tj顶点上；③刀触点Pi在三角网格Tj边上。

[0053] 其中，判断刀触点Pi所属的三角形网格Tj的坐标及其在该三角形网格中的位置方法如下：

[0054] 以某一刀触点Pi为例，该刀触点与任意三角网格Tj的三个顶点Aj、Bj、Cj构成向量分别计算三个向量之间的叉积。

[0055]

[0056]

[0057]

[0058] 若t1、t2、t3同号，则刀触点Pi在三角形网格Tj内部；若t1、t2、t3其中一个为0，另两个同号， 则刀触点Pi在边AjBj或BjCj或AjCj上；若t1、t2、t3其中两个为0，则刀触点Pi在顶点Aj或Bj或 Cj上；否则不在三角形网格Tj上。

[0059] 分类后得到点集P中每一个刀触点所属的三角形网格集T，及相应的顶点坐标集A、B、C和 t1、t2、t3值。依据三角形顶点坐标及分类情况计算刀触点Pi的法向矢量ni，计算方法如下：

[0060] 如图3a所示，对于第①类中刀触点Pi的法向矢量ni即为所属三角网格Tj所确定切平面的法向 矢量Nj。

[0061] 如图3b所示，第②类中Pi点的法向矢量ni受周围三角网格Tj～Tj+5的法向矢量Nj～Nj+5影响。 引入形状权重因子Wj，Wj依据三角形形状分布进行取值，求出法向矢量ni：

[0062]

[0063] 其中Aj和Nj分别为周围三角网格Tj的面积和法向矢量。

[0064] 如图3c所示，第③类中Pi点的法向矢量ni受所属三角网格边PmPn的两个顶点Pm、Pn位置的 法向矢量Nm、Nn影响，受相邻两个网格Ta、Tb的法向矢量Na、Nb影响。所以法向矢量ni计算如 下：

[0065] ni＝(1-u)Nm+uNn+u(1-u)Ncomp，u∈[0,1]

[0066]

[0067] 其中Ncomp为Ta、Tb的补偿矢量，u为Pi点到Pn点的长度。

[0068] 由于球形刀底部切削速度低、容屑空间小、切削条件较差，所以在实际加工中，通常刀尖点不 能直接用作刀触点，刀轴矢量Nki不能与零件表面法向矢量ni重合。在已经确定复杂曲面法向矢量 ni后，将法矢ni在空间中围绕某单位向量Fi(ai,bi,ci)旋转角度 即进行Ri转换(Ri为罗德里格旋 转矩阵)，得到球形刀的刀轴矢量Nki(Nxi,Nyi,Nzi)。

[0069] Νki＝Ri·ni

[0070]

[0071] 其中，单位向量Fi(ai,bi,ci)确定方法如下：

[0072] 将刀触点Pi(xi,yi,zi)投影到XOY面上，即(zi＝0)，Pi与坐标系O点相连得向量Ei(-xi,-yi,0)，将 向量Ei单位化后围绕Z轴逆时针旋转90°后得到单位向量Fi：

[0073]

[0074]

[0075] 通过空间旋转得到的刀轴矢量集Nk可能存在相邻刀轴矢量夹角过大，导致旋转非线性误差过 大，影响加工质量，造成过切或表面奇异等问题。因此需要对刀轴旋转运动而带来的非线性误差进 行控制，保证复杂曲面表面加工精度。

[0076] 确定相邻刀位点之间加工运动所产生的最大非线性误差角αi方法如下：

[0077] 根据本发明实施例中AC式五轴机床结构及运动链传递如图2所示，该AC式五轴机床由三个 直线运动轴(X、Y、Z轴)和两个旋转运动轴(A、C轴)构成，A轴绕X轴旋转，C轴绕Z轴 旋转，待加工工件安装在C轴工作台上。其中OWXWYWZW为工件固联坐标系，OTXTYTZT为刀具 固联坐标系，OAXAYAZA为与旋转轴A固联的坐标系，OCXCYCZC为与旋转轴C固联的坐标系，rA为OAXAYAZA在刀具坐标系中位置，rC为OCXCYCZC在工件坐标系中位置，L为刀具长度。在生成 工件坐标系中刀具加工数据，即刀触点集P和刀轴矢量集Nk后，根据五轴机床结构，经过平移旋 转逆运动变换后，反求刀具坐标系中运动坐标。其中，该五轴机床正向运动链为刀具、A轴旋转、 机床床身、Y轴、Z轴、X轴、C轴旋转、工件。假设刀具坐标系中初始刀触点和刀轴矢量分别为 (0,0,0)和(0,-1,0)，变换后的工件坐标系刀触点和刀轴矢量分别为(x,y,z)和(Nx,Ny,Nz)，机床各直线运 动轴相对于初始状态位置为rs，旋转轴A、C相对于初始角度为a和c。正向变换如下：

[0078]

[0079] 其中得到的刀轴矢量(Nx,Ny,Nz)为：

[0080]

[0081] 本发明分析选取加工中某一刀触点Pi运动至下一刀触点Pi+1，计算其刀轴矢量Nki至Nki+1过 程中Nkp处产生的最大非线性误差角αi，如图5所示。Nkp’为刀矢Nkp在实际加工中机床坐标系的 表达，即Nkp与Nkp’的最大夹角。

[0082] Nkp利用空间四元数球面线性法表示如下：

[0083]

[0084] 如图5所示，θi为相邻刀轴矢量Nki和Nki+1之间夹角。tθi为刀轴矢量Nkp与Nki的夹角，t∈[0,1]。

[0085] 实际加工中刀具坐标系中刀轴矢量Nkp’为：

[0086] Nkp'(t)＝(-cos(a(t))sin(c(t)) -cos(a(t))cos(c(t)) -sin(a(t)))

[0087] 设U(t)＝cos(α(t))，则非线性误差角αi表示如下：

[0088]

[0089] Nk1·Nkp'(t)＝cos(ai)cos(a(t))cos(t(ci+1-ci))+sin(ai)sin(a(t))

[0090] Nk2·Nkp'(t)＝cos(ai+1)cos(a(t))cos((1-t)(ci+1-ci))+sin(ai+1)sin(a(t))[0091] 当U(t)求得最小值时，非线性误差角α(t)＝arccos(U(t))求得最大值，并通过t值确定非线性误差 角最大处的位置。对U(t)关于t求导，为了简便计算，假设Pi位于初始位置，ai＝-π/2，ci＝0，则 cos(θi)＝cos(Δa)。利用泰勒公式展开求导式至三次方，当t＝1/2时，  所以当刀具运动至Pi和Pi+1中间位置时， 刀轴矢量Nki至Nki+1间的误差角αi最大。

[0092] 且当t＝1/2时，U(t)和α(t)结果如下：

[0093]

[0094]

[0095] 所以当相邻刀轴矢量之间夹角θi过大时，实际加工中非线性误差角αi也会过大，从而影响复 杂曲面表面加工质量。

[0096] 通过四元数球面线性插补法缩小θi而控制误差角αi的方法如下：

[0097] 在相邻刀轴矢量Nki和Nki+1间非线性误差角αi最大位置处，即t＝ti＝1/2，利用四元数球面线性 插补法插入新的刀轴矢量Nki+t，控制刀轴矢量Nki、Nki+t、Nki+1之间的夹角小于θmax。θmax为相邻 刀轴矢量间允许最大夹角，由允许最大非线性误差角αmax计算得到。

[0098]

[0099] 将新得到的刀触点集P和刀轴矢量集Nk转化为可用于球形刀加工的刀位点集O和刀轴矢量集 Nk，通过工件坐标系和机床坐标系转换，得到五轴机床加工的各轴代码{X,Y,Z,A,C}。

[0100] 以上所述的实施例仅是对本发明的较佳实例进行阐述，并非限制本发明，在不脱离本发明刀轴 矢量规划方法及思路的前提下，本领域的工程技术人员对本发明技术方案作出的各种变形和改进， 均应落入本发明的保护范围内。