CN113467376B - 一种面向多加工场景的多轴轨迹压缩方法 - Google Patents

Abstract

本发明公开了一种面向多加工场景的多轴轨迹压缩方法，属于数控加工领域。本发明采用C2连续且保证单调性的三次有理样条来进行刀尖轨迹与刀轴轨迹之间的参数同步，C2连续保证刀轴矢量的变化连续，并且加速度连续，能使旋转轴运动达到加速度连续；单调保证刀轴点样条参数随刀尖点样条参数增大而增大或不变，不会出现刀轴点回退的情况。本发明实现了误差计算与轨迹压缩方法的解耦，当新的加工场景出现时，只需要设计新的允差空间和对应OTE计算方法即可，而无需修改轨迹压缩流程中的任何一步。本发明利用目标压缩比结合特征曲率积分，不需要迭代，能达到快速选取初始拟合节点的效果，提高压缩效率。

Description

一种面向多加工场景的多轴轨迹压缩方法

技术领域

本发明属于数控加工领域，更具体地，涉及一种面向多加工场景的多轴轨迹压缩方法。

背景技术

在数控加工领域，要得到多轴加工程序，首先需要由CAM(Computer AidedManufacturing)软件根据零件模型并结合工艺参数生成加工轨迹，然后对其离散，所以，多轴加工程序通常是由G01小线段组成，这会导致以下缺陷：相邻小线段之间只达到了G0连续性，这会导致加工过程中进给速度、加速度波动，使机床振动，加工质量变差。同时，为了达到较高的加工精度，有时会让G01小线段非常短，这会使得加工效率较低，并且程序信息过大，难以满足我们对高速加工的诉求。因此，将原始小线段程序压缩为连续性更好，单位程序段信息量更大的参数曲线后再进行加工，可以有效地改善上述问题。但是，这对于多轴***而言绝非易事。

多轴机床相较三轴机床，由于旋转轴的引入，不仅可以控制刀具位置，还可以控制刀具姿态，可以实现更复杂的加工，因此有了更多的加工场景。目前除了广泛应用于汽车、航空工业中的叶轮、叶片和模具的自由曲面铣削，还应用于航天结构件等外壳的复合材料蒙皮铺放、缠绕加工。加工场景的不同，意味着加工特点的不同，即误差计算方法可以不同，以最大程度地利用具体的加工特征。目前，在中国发明专利CN105425725A中，使用Hausdorff距离描述拟合轨迹与程序折线段的偏差，在非专利文献《五轴刀具轨迹的六维B样条拟合算法研究》(闵康,2018)中，进一步定义了五轴轨迹的Hausdorff距离，并且从刀具的角度建立了误差计算方法。但是这些方法均是以金属切削加工为前提的，并且与轨迹压缩算法耦合太紧，导致轨迹压缩方法无法适用多种加工场景。因此，需要提出一套误差计算标准接口，与轨迹压缩方法分离，不同的加工场景设计不同的误差计算方法，然后配合轨迹压缩方法使用。

对于多轴加工轨迹，除了需要对刀尖点进行压缩，还需要对刀轴矢量进行压缩，以平滑旋转轴运动，因此，若这两条轨迹并非在同一参数域下进行的压缩，则还需要构建它们的参数同步关系。在非专利文献《Dual NURBS Path Smoothing for 5-Axis Linear Pathof Flank Milling》(Li D,Zhang W,Zhou W.International Journal of PrecisionEngineering and Manufacturing,2018)中，使用了线性插值的方式来进行两条轨迹之间的同步，这种同步方式只有G0的连续性，无法让旋转轴运动高阶连续。因此，需要构建一个双样条参数同步模型，满足高阶连续和单调非递减的约束。

为了实现多轴轨迹的快速压缩，初始节点的快速选取非常重要，在非专利文献《Fast Automatic Knot Placement Method for Accurate B-spline Curve Fitting》(Yeh R,Nashed Y S G,Peterka T.Computer-Aided Design,2020)中，利用离散点的高阶导数积分，快速选取拟合节点，但是该方法只适用于三维轨迹，无法用于五轴加工轨迹数据。在非专利文献《An Adaptive Curvature-Guided Approach for the Knot-PlacementProblem in Fitted Splines》(Aguilar E,Elizalde H,Diego Cárdenas.Journal ofComputing and Information Science in Engineering,2018)中，利用曲率指导节点选取，但是其存在迭代，效率不高，并且同样也只适用于三维轨迹，无法用于五轴加工轨迹数据。因此，需要提出一个基于曲率的，快速的，可用于五轴轨迹压缩的初始节点选取方法。

综上所述，目前的轨迹压缩方法现状如下：①误差计算模型普遍更加适用于金属切削加工，并且与轨迹压缩方法耦合严重。②多轴轨迹压缩中，若刀尖轨迹和刀轴轨迹不在同一参数域下，则目前方法难以保证旋转轴运动的高阶连续性。③对于压缩的初始节点选取，目前的方法不能高效且很好地用在多轴轨迹压缩中。

发明内容

针对现有技术的以上缺陷或改进需求，本发明提供了一种面向多加工场景的多轴轨迹压缩方法，其目的在于使得误差计算、初始节点选取、刀尖与刀轴矢量的参数同步能够适用多种压缩场景。

为实现上述目的，本发明提供了一种面向多加工场景的多轴轨迹压缩方法，包括：

S1.对原始刀位点进行区间划分；

S2.选取每个刀位点子区间的初始压缩节点序列；

S3.对每个初始压缩节点序列中的刀位点分别计算其节点矢量，并采用C2连续且单调的三次有理样条构建刀尖曲线与刀轴曲线的参数同步关系；

S4.根据每个子区间初始压缩节点序列中的节点、节点矢量和首末端点切矢，反算控制点，进而得到刀尖B样条轨迹，再结合刀尖曲线与刀轴曲线的参数同步关系得到刀轴B样条轨迹；

S5.计算每个刀尖点在刀尖点B样条轨迹上的垂足到刀尖点的误差，若所有误差均小于0，则进入步骤S6，否则选择误差最大的刀位点，将其加入到初始压缩节点序列，返回执行步骤S3；

S6.计算压缩后所有刀尖点对应的刀轴角度偏差；若所有刀轴角度偏差均小于给定允许值，则进入步骤S7，否则选择刀轴角度偏差最大的刀位点，将其加入到初始压缩节点序列，返回执行步骤S3；

S7.遍历所有压缩节点，若相邻压缩节点之间没有其它刀位点，则计算相邻压缩节点之间的弦高误差，若弦高误差均满足设定要求，则输出步骤S4压缩得到的B样条轨迹；否则，针对弦高误差最大的一段，重新进行区间划分，各个区间再依次进入步骤S2中进行压缩。

进一步地，步骤S5具体为，计算每个刀尖点在刀尖点B样条轨迹上的垂足到刀尖点所构成的允差空间的超允差误差，若所有超允差误差均小于0，则进入步骤S6，否则选择超允差误差最大的刀位点，将其加入到初始压缩节点序列，返回执行步骤S3；其中，允差空间为刀位点在满足误差要求下可以调整的位置所构成的空间；超允差误差为任意一空间点到允差空间外表面的距离。

进一步地，对于减材制造，其允差空间是以刀位点为球心，以允许轮廓误差大小为半径的球型空间。

进一步地，超允差误差为任意一空间点到球形空间外表面的距离。

进一步地，对于增材加工，其允差空间为圆柱体，所述圆柱体底面圆半径为切向允差绝对值，高为正允差与负允差绝对值之和；其中，切向允差为加工轮廓允差，方向为允差向量的正交方向；允差向量为刀位点在工件表面上的正交投影点指向刀位点的单位矢量；正允差为允差向量方向的允差；负允差为允差向量反方向的允差。

进一步地，超允差误差的计算过程具体为，

若空间任意一点P与第i个刀尖点C_i的连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的底面，且α小于π/2，超允差误差为：

||C_iP||表示C_i与P之间的距离，D_n表示正允差，α表示连线C_iP与允差向量N_i之间夹角，

表示刀轴矢量在轨迹切矢方向上的投影模长，

表示刀轴矢量，

表示第i个刀尖点切矢；

若连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的圆柱面，且α小于π/2，超允差误差为：

D_t表示切向允差；

若连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的圆柱面，且α大于π/2，超允差误差为：

若连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的底面，且α大于π/2，超允差误差为：

D_-n表示负允差。

进一步地，步骤S2具体为根据目标压缩比和特征曲率选取每个刀位点子区间的初始压缩节点序列，包括以下步骤：

将子区间中第一个刀位点加入压缩节点序列Q_comp中；

计算每个刀尖点的离散曲率

每个刀轴点的离散曲率

以及每个刀位点的特征曲率

Q_i表示第i个刀位点，C_i表示第i个刀尖点，a≤i≤b，a表示初始压缩节点序列第一个刀位点在原始刀位点集合中的下标，b表示初始压缩节点序列最后一个刀位点在原始刀位点集合中的下标；

计算每个刀位点的特征曲率积分

计算相邻刀位点特征曲率积分的变化量

R为目标压缩比；

遍历压缩区间中除第一个刀位点和最后一个刀位点的中间刀位点，若G(j)-ΔG大于等于压缩节点序列Q_comp中最后一个刀位点的特征曲率积分，则将第j个刀位点加入压缩节点序列；

将最后一个刀位点加入到压缩节点序列，得到初始压缩节点序列。

进一步地，步骤S1包括：

计算所有相邻刀尖点之间的距离，与给定长直线长度l_max比较；

若距离大于l_max，则判定相邻刀尖点之间连线为长直线，按照长直线左侧、长直线自身、长直线右侧的规则划分为三个区间。

进一步地，步骤S7所述针对弦高误差最大的一段，重新进行区间划分，具体为，若左半边弦高误差更大，则选择左边的刀位点为新的区间划分点，若右半边弦高误差更大，则选择右边的刀位点为新的区间划分点，新的区间划分点将当前压缩区间一分为二。

总体而言，通过本发明所构思的以上技术方案与现有技术相比，能够取得下列有益效果。

1、本发明采用C2连续且保证单调性的三次有理样条来进行刀尖轨迹与刀轴轨迹之间的参数同步，C2连续保证刀轴矢量的变化连续，即刀轴点速度连续，并且加速度连续，能使旋转轴运动达到加速度连续；单调保证刀轴点样条参数随刀尖点样条参数增大而增大或不变，不会出现刀轴点回退的情况，即刀具前后摆动的现象。

2、本发明实现了误差计算与轨迹压缩方法的解耦，当新的加工场景出现时，只需要设计新的允差空间和对应OTE计算方法即可，而无需修改轨迹压缩流程中的任何一步。其中，增材加工误差计算过程，能够更有效地利用增材加工的特点，即刀位点在圆柱型允差空间内产生偏移，也有相同的加工效果，扩大了刀位点可以调整的范围，从而提高压缩比。

3、现有方法多采用迭代选取初始压缩节点，而本发明利用目标压缩比结合特征曲率积分，不需要迭代，能达到快速选取初始拟合节点的效果，提高压缩效率。

附图说明

图1为金属切削加工的允差空间及超允差误差示意图；

图2为增材加工的允差空间示意图；

图3为允差向量定义示意图；

图4为刀位轨迹在刀位点处的法平面示意图；

图5为刀轴矢量与轨迹切矢所在的平面示意图；

图6为允差向量所在位置示意图；

图7为空间点在允差空间上底面时的示意图；

图8为空间点在允差空间侧面的上半部的示意图；

图9为空间点在允差空间侧面的下半部的示意图；

图10为空间点在允差空间下底面时的示意图；

图11为轨迹压缩总体流程图；

图12为根据长直线段划分区间示意图；

图13为对单个区间进行压缩的流程图；

图14为相邻压缩节点之间没有其它数据点时的示意图；

图15为相邻压缩节点之间左半段弦高误差过大时的示意图，此时选择左边的刀位点为新的区间划分点；

图16为相邻压缩节点之间右半段弦高误差过大时的示意图，此时选择右边的刀位点为新的区间划分点。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。此外，下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。

本发明的主要改进点体现在以下几方面：

(1)首先定义允差空间(Tolerance Space,TS)为刀位点在满足误差要求下可以调整的位置所构成的空间。定义超允差误差(Over Tolerance Error,OTE)为任意一空间点到允差空间外表面的距离，若空间点在允差空间内，则OTE为0。不同的允差空间需要提供统一的计算OTE的接口，然后在轨迹压缩方法中，所有需要用到点到点的误差计算，均使用该接口，而无需了解具体的允差空间设计。这里给出常用的金属切削减材加工的允差空间和增材加工的允差空间。

金属切削为减材制造，针对一个刀位点而言，无论往哪个方向偏离，对加工轨迹的影响都是相同的，因此其允差空间为圆心是刀位点，半径是允许轮廓误差大小的球型空间，其形状以及OTE的计算，均可见附图1。

对于铺放这种增材加工而言，刀位点在材料增厚方向的偏离不会造成加工误差，在材料增厚方向的两侧偏移则会造成加工误差，于是有如附图2所示的圆柱型允差空间，其中允差向量定义为刀位点在工件表面上的正交投影点指向刀位点的单位矢量，如附图3，即增材加工材料增厚的方向；正允差D_n则为允差向量方向的允差；负允差D_-n为允差向量反方向的允差；切向允差D_t为加工轮廓允差，方向为允差向量的正交方向。由于刀位文件中没有原始工件表面信息，所以无法直接获得允差向量，接下来给出根据刀位点计算允差向量的方法：

①使用Boehm法计算各刀尖点切矢

②计算刀位轨迹在刀位点处的法平面。一个点和平面在该点处的法矢即可描述一个平面，此点即为刀位点，刀位点法平面的法矢即为刀路轨迹在此处的切矢，所以有如图4平面：

其中，C_i为第i个刀位点，

为第i个刀尖点轨迹切矢。

③计算允差向量。由于铺放加工中刀轴矢量与轨迹切矢组成的平面正交于工件表面，又刀位点在该平面上，所以允差向量也必然在刀轴矢量与轨迹切矢组成的平面上，所以允差向量所在直线为该平面与轨迹法平面的交线，如图5。显然刀轴矢量不会是指向工件内部的，所以允差向量方向为与刀轴矢量夹角较小的方向，如图6所示。所以刀轴矢量在轨迹切矢方向上的投影模长为

允差向量为

接下来给出圆柱型允差空间的OTE计算方法，首先有第i个刀尖点C_i，空间任意一点P，于是有C_iP与允差向量N_i之间夹角为

然后分为4种情况：

①如图7所示，空间点与刀尖点的连线C_iP与允差空间的交点在允差空间圆柱体的底面，且α小于π/2。此时OTE为：

②如图8所示，空间点与刀尖点的连线C_iP与允差空间的交点在允差空间圆柱体的圆柱面，且α小于π/2。此时OTE为：

③如图9所示，空间点与刀尖点的连线C_iP与允差空间的交点在允差空间圆柱体的圆柱面，且α大于π/2。此时OTE为：

④如图10所示，空间点与刀尖点的连线C_iP与允差空间的交点在允差空间圆柱体的底面，且α大于π/2。此时OTE为：

(2)本发明率先采用C2连续且保证单调性的三次有理样条来进行刀尖轨迹与刀轴轨迹之间的参数同步。

首先有刀尖轨迹的节点矢量U＝{0,0,0,0,...,u_i,...,1,1,1,1}，刀轴轨迹的节点矢量W＝{0,0,0,0,...,w_i,...,1,1,1,1}，为简化过程表达式，特令h_i＝u_i+1-u_i，

有参数同步关系为

p_i(u)＝w_i(1-θ)₃+(w_i(3+γ_i)+h_id_i)θ(1-θ)²+(w_i+1(3+γ_i)-h_id_i+1)θ²(1-θ)+w_i+1θ³

q_i(u)＝(1-θ)²+θ(1-θ)(2+γ_i)+θ²

其中

d_i则通过求解以下线性方程组得到：

h_id_i-1+(h_i(2+γ_i-1)+h_i-1(2+γ_i))d_i+h_i-1d_i+1＝h_i-1(3+γ_i)Δ_i+h_i-1(3+γ_i-1)Δ_i-1

(3)本发明提供一种根据目标压缩比和曲率快速得到初始压缩节点序列的方法。首先计算刀尖点和刀轴点的离散曲率。

刀尖点离散曲率为

刀轴点离散曲率为

其中C_i为第i个刀尖点，O_i为第i个刀轴点。定义特征曲率(Feature Curvature，FC)为

定义特征曲率积分(Cumulative Feature Curvature，CFC)为特征曲率在刀尖点弧长上的积分：

其中，n为刀位点数目。

目标压缩比即为期望达到的压缩比，显然对于原始刀位点不能无限压缩，因此目标压缩比不能无限大；给定的目标压缩比为R，计算所有点的CFC，每增加

就选取当前节点为压缩的初始节点，最后将第一个点也加入压缩的初始节点，即完成了初始节点的选取。

本发明选取初始节点的方法不同于现有一般方法，现有方法多采用迭代选取初始拟合节点，而本发明利用目标压缩比结合特征曲率积分，不需要迭代，能达到快速选取初始拟合节点的效果。

本发明结合前三点，提供了一种面向多加工场景的多轴轨迹压缩方法，总体流程如图11，包括以下步骤：

(1)进行区间划分

原始刀位点为{Q₀,...,Q_n}，n表示刀位点个数，其中第i个刀位点Q_i包含第i个刀尖点C_i和第i个刀轴点O_i，即

给定长直线长度为l_max，计算所有相邻刀尖点之间的距离l_i＝||C_i-1C_i||,1≤i≤n，若l_i≥l_max，则判定C_i-1C_i为长直线，然后将长直线左侧、长直线自身、长直线右侧划分为三个区间。使用该方法可以将原始刀位点{Q₀,...,Q_n}划分为若干个子区间{Q₀,...,Q_a},{Q_a,...,Q_b},...,{Q_z,...,Q_n}。不失一般性地，下述步骤均以子区间{Q_a,...,Q_b}为例。

上述区间划分方式仅作为实施例，不同的区间划分方式可以得到一系列子区间，而区间划分方式并不影响对单个子区间实施本发明所述的多轴轨迹压缩方法。

(2)选取每个子区间的初始压缩节点；

Step1：计算每个刀尖点的离散曲率

Step2：计算每个刀轴点的离散曲率

Step3：计算每个刀位点的特征曲率

Step4：计算

R为目标压缩比，n为原始刀位点个数；

Step5：将第一个刀位点加入压缩节点序列Q_comp中，则Q_comp＝{Q_a}；

Step6：计算每个刀位点的特征曲率积分

Step7：遍历压缩区间的刀位点，若G(j)-ΔG大于等于压缩节点序列Q_comp中最后一个刀位点的特征曲率积分，则将第j个刀位点加入压缩节点序列，则Q_comp＝{Q_a,...,Q_j}。

Step8：将最后一个刀位点加入到压缩节点序列，则得到了初始压缩节点序列Q_comp＝{Q_a,...,Q_b}。

(3)对压缩节点序列Q_comp中的刀位点使用累计弦长参数化方法分别计算每个刀尖曲线和刀轴曲线节点矢量。

(4)构建刀尖曲线与刀轴曲线的参数同步关系

Step1：有刀尖曲线的节点矢量U＝{0,0,0,0,...,u_i,...,1,1,1,1}，刀轴曲线的节点矢量W＝{0,0,0,0,...,w_i,...,1,1,1,1}，为简化过程表达式，特令h_i＝u_i+1-u_i，

γ_i＝1，h_i、θ、Δ_i、γ_i为过程变量，无实际物理意义；

Step2：求解线性方程组：

h_id_i-1+(h_i(2+γ_i-1)+h_i-1(2+γ_i))d_i+h_i-1d_i+1＝h_i-1(3+γ_i)Δ_i+h_i-1(3+γ_i-1)Δ_i-1

d₀＝d_b-a＝1

可以得到的d_i的值，d_i为过程变量，无实际物理意义。判断

是否满足，若满足则进入Step3；否则令

然后重复Step2；

Step3：得到参数同步关系为

p_i(u)＝w_i(1-θ)₃+(w_i(3+γ_i)+h_id_i)θ(1-θ)²+(w_i+1(3+γ_i)-h_id_i+1)θ²(1-θ)+w_i+1θ³

q_i(u)＝(1-θ)²+θ(1-θ)(2+γ_i)+θ²

目前普遍使用的参数同步方法有两种，一是直接将刀尖点样条的参数作为刀轴点样条的参数，该方法会在每个节点处存在误差e_i＝w_i-u_i；另外一种方法是基于线性插值的参数同步模型，即让u_i和w_i之间为线性映射，这会导致映射关系在节点处一阶导数不连续，最终导致刀轴矢量的变化不连续。

本发明使用C2连续且单调的有理三次样条来进行刀尖轨迹与刀轴轨迹之间的参数同步，C2连续保证刀轴矢量的变化连续，即刀轴点速度连续，并且加速度连续，单调保证刀轴点样条参数随刀尖点样条参数增大而增大或不变，不会出现刀轴点回退的情况，即刀具前后摆动的现象。

(5)使用Boehm法计算压缩子区间的首末点端矢。

(6)将压缩节点序列Q_comp作为插值点，根据节点、节点矢量、首末端矢，反算得到B样条控制点，于是能得到刀尖B样条轨迹Q_c(u)、刀轴B样条轨迹Q_o(w)。

(7)找到第j个刀尖点C_j在刀尖点曲线上的垂足F_j，然后在刀尖点处构建允差空间，然后计算垂足到刀尖点所构成的允差空间的OTE＝||C_iF_i||-d，其中d为允差。若加工类型为增材加工，对应的OTE计算方法可见发明内容中的(1)。若所有OTE均小于0，即没有刀尖点超差，则进入(8)；否则选择OTE最大的刀位点Q_i，加入到压缩节点序列Q_comp中，得到下一次迭代的压缩节点序列Q_comp＝{Q_a,...Q_i...,Q_b}，然后进入步骤(3)。

(8)根据第j个刀尖点在刀尖点曲线上垂足的参数u_j，通过参数同步关系w＝f_u-w(u)计算得到刀轴点的参数w_j，于是可以得到压缩后第j个刀尖点对应的刀轴矢量T_j＝Q_o(w_j)-Q_c(u_j)，而原始第j个刀尖点的原始刀轴矢量为O_j-C_j，计算它们之间的夹角，即可得到刀轴角度偏差。若所有刀轴角度偏差均小于给定允许值

则进入步骤(9)；否则选择刀轴角度偏差最大的刀位点Q_i，加入到压缩节点序列Q_comp中，得到下一次迭代的压缩节点序列Q_comp＝{Q_a,...Q_i...,Q_b}，然后进入步骤(3)。

(9)遍历所有压缩节点，若存在相邻压缩节点之间如图14没有其它刀位点，则计算相邻压缩节点之间的弦高误差，若弦高误差均满足要求，则输出当前压缩得到的样条曲线；否则，找到弦高误差最大的一段，判断其是如图15左半边弦高误差更大，还是如图16右半边弦高误差更大，若是左半边弦高误差更大，则选择左边的刀位点为新的区间划分点，若是右半边弦高误差更大，则选择右边的刀位点为新的区间划分点。新的区间划分点将当前压缩区间一分为二，然后再依次进入步骤(2)中进行压缩。

通过以上具体实施方法可以看出，本算法计算简单，效率较高，可以适用于多种加工场景的多轴轨迹压缩，实用性高。

本领域的技术人员容易理解，以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (9)

1.一种面向多加工场景的多轴轨迹压缩方法，其特征在于,包括：

S1.对原始刀位点进行区间划分；

S2.选取每个刀位点子区间的初始压缩节点序列；

S3.对每个初始压缩节点序列中的刀位点分别计算其节点矢量，并采用C2连续且单调的三次有理样条构建刀尖曲线与刀轴曲线的参数同步关系；

S4.根据每个子区间初始压缩节点序列中的节点、节点矢量和首末端点切矢，反算控制点，进而得到刀尖B样条轨迹，再结合刀尖曲线与刀轴曲线的参数同步关系得到刀轴B样条轨迹；

S5.计算每个刀尖点在刀尖点B样条轨迹上的垂足到刀尖点的误差，若所有误差均小于0，则进入步骤S6，否则选择误差最大的刀位点，将其加入到初始压缩节点序列，返回执行步骤S3；

S6.计算压缩后所有刀尖点对应的刀轴角度偏差；若所有刀轴角度偏差均小于给定允许值，则进入步骤S7，否则选择刀轴角度偏差最大的刀位点，将其加入到初始压缩节点序列，返回执行步骤S3；

S7.遍历所有压缩节点，若相邻压缩节点之间没有其它刀位点，则计算相邻压缩节点之间的弦高误差，若弦高误差均满足设定要求，则输出步骤S4压缩得到的B样条轨迹；否则，针对弦高误差最大的一段，重新进行区间划分，各个区间再依次进入步骤S2中进行压缩。

2.根据权利要求1所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，步骤S5具体为，计算每个刀尖点在刀尖点B样条轨迹上的垂足到刀尖点所构成的允差空间的超允差误差，若所有超允差误差均小于0，则进入步骤S6，否则选择超允差误差最大的刀位点，将其加入到初始压缩节点序列，返回执行步骤S3；其中，允差空间为刀位点在满足误差要求下可以调整的位置所构成的空间；超允差误差为任意一空间点到允差空间外表面的距离。

3.根据权利要求2所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，对于减材制造，其允差空间是以刀位点为球心，以允许轮廓误差大小为半径的球型空间。

4.根据权利要求3所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，超允差误差为任意一空间点到球形空间外表面的距离。

5.根据权利要求2所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，对于增材加工，其允差空间为圆柱体，所述圆柱体底面圆半径为切向允差绝对值，高为正允差与负允差绝对值之和；其中，切向允差为加工轮廓允差，方向为允差向量的正交方向；允差向量为刀位点在工件表面上的正交投影点指向刀位点的单位矢量；正允差为允差向量方向的允差；负允差为允差向量反方向的允差。

6.根据权利要求5所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，超允差误差的计算过程具体为，

若空间任意一点P与第i个刀尖点C_i的连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的底面，且α小于π/2，超允差误差为：

||C_iP||表示C_i与P之间的距离，D_n表示正允差，α表示连线C_iP与允差向量N_i之间夹角，

表示刀轴矢量在轨迹切矢方向上的投影模长，

表示刀轴矢量，

表示第i个刀尖点切矢；

若连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的圆柱面，超允差误差为：

D_t表示切向允差；

若连线C_iP与允差空间圆柱体的交点在允差空间圆柱体的底面，且α大于π/2，超允差误差为：

D_-n表示负允差。

7.根据权利要求1-6任一项所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，步骤S2具体为根据目标压缩比和特征曲率选取每个刀位点子区间的初始压缩节点序列，包括以下步骤：

将子区间中第一个刀位点加入压缩节点序列Q_comp中；

计算每个刀尖点的离散曲率

每个刀轴点的离散曲率

以及每个刀位点的特征曲率

O_i为第i个刀轴点，C_i表示第i个刀尖点，a≤i≤b，a表示初始压缩节点序列第一个刀位点在原始刀位点集合中的下标，b表示初始压缩节点序列最后一个刀位点在原始刀位点集合中的下标；

计算每个刀位点的特征曲率积分

计算相邻刀位点特征曲率积分的变化量

R为目标压缩比；

遍历压缩区间中除第一个刀位点和最后一个刀位点的中间刀位点，若G(j)-△G大于等于压缩节点序列Q_comp中最后一个刀位点的特征曲率积分，则将第j个刀位点加入压缩节点序列；

将最后一个刀位点加入到压缩节点序列，得到初始压缩节点序列。

8.根据权利要求1所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，步骤S1包括：

计算所有相邻刀尖点之间的距离，与给定长直线长度l_max比较；

若距离大于l_max，则判定相邻刀尖点之间连线为长直线，按照长直线左侧、长直线自身、长直线右侧的规则划分为三个区间。

9.根据权利要求1所述的一种面向多加工场景的多轴轨迹压缩方法，其特征在于，步骤S7所述针对弦高误差最大的一段，重新进行区间划分，具体为，若左半边弦高误差更大，则选择左边的刀位点为新的区间划分点，若右半边弦高误差更大，则选择右边的刀位点为新的区间划分点，新的区间划分点将当前压缩区间一分为二。

Images