CN104238456B - 一种非球头刀铣削加工自由曲面的方法 - Google Patents

Abstract

一种非球头刀铣削加工自由曲面的方法，其特征是采用摆定前倾角和零度旋转角的加工方式，构造刀具有效切削曲面，以张量的方式度量曲面上任意切触点处的切削宽度，构建非球头刀切削宽度二阶张量场，提取张量场中的三分退化点，并以此为起点搜索曲面内部边界点，实现对曲面加工区域的划分，每个区域称为一个曲面加工特征。在每一个曲面加工特征中构建数条曲线，使曲线上每点的切向与该点处切宽最大的进给方向重合，选出平均切宽最大的曲线为该区域内的初始加工轨迹线，采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。本发明能够减少刀轨长度，实现加工效率的提高。

Description

一种非球头刀铣削加工自由曲面的方法

技术领域

本发明涉及一种CAD/CAM技术，尤其是一种自由曲面的CAD/CAM技术，具体地说是一种非球头刀铣削加工自由曲面的方法

背景技术

随着产品整体化、轻量化设计以及对气动性能要求的不断提升，越来越多的复杂自由曲面出现在航空航天、汽车及船舶等产品的零件中，由于这些复杂曲面几何性质复杂，曲率变化大，传统上采用球头刀铣削的加工方式，但在实际生产中球头刀的切削宽度较小，加工效率较低，例如在大型模具的复杂自由曲面加工中，其精加工时间往往需要100小时以上。相比之下五轴联动非球头刀的切削通过改变刀具姿态可以实现较大的切宽，加工效率明显提高，因而成为复杂自由曲面加工的高效加工方式。但非球头刀的切削宽度受刀具姿态以及进给方向的影响，而五轴加工刀具姿态的确定依赖人的经验且缺乏灵活性，在实际加工中往往切削宽度达不到预期效果，仍然需要进行切削宽度优化。

切削宽度优化是优化曲面加工效率的重要方式之一，较大的切削宽度可以减少曲面加工轨迹的长度，在不考虑机床运动学限制的情况下可以有效缩短加工时间，提高加工效率。一些学者通过以曲面最长边界作为初始刀轨进行偏置生产刀轨，或者以曲面主曲率最大的对应方向作为进给方向从而实现切削宽度的优化，还有目前较为成熟的美国北卡罗来纳州立大学Yuan-shin Lee教授提出的基于Machining Potential Field的方法，即首先在曲面上构造一系列曲线，选出平均切宽最大的曲线作为初始轨迹线，偏置初始轨迹线生成曲面的加工轨迹，当偏置得到的轨迹线的切削效率系数小于设定值时则停止轨迹偏置，在未生成刀轨区域重新构造初始轨迹线进行轨迹计算。该方法能够一定程度上减少相邻刀轨间的重叠区域，实现切宽的优化，但该方法的本质是一种被动式的分区加工方法，仍然受人经验的影响，是一种贪婪式算法，因此仅能得到局部优化的加工结果。

本发明公开了一种非球头刀铣削加工自由曲面方法。针对自由曲面五轴联动铣削加工中采用摆定前倾角和零度旋转角的加工方式，构造刀具有效切削曲面，以张量的方式度量曲面上任意切触点处的切削宽度，构建非球头刀切削宽度二阶张量场，提取张量场中的三分退化点，并以此为起点搜索曲面内部边界点，实现对曲面加工区域的划分，每个区域称为一个曲面加工特征。在每一个曲面加工特征中构建数条曲线，使曲线上每点的切向与该点处切宽最大的进给方向重合，选出平均切宽最大的曲线为该区域内的初始加工轨迹线，采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。该方法能够保证所生成的平底刀加工轨迹尽可能遵循对应最大切宽的进给方向在曲面范围内的区域化分布规律，从而减少刀轨长度，实现加工效率的提高。

发明内容

本发明的目的是针对当前复杂自由曲面五轴联动铣削加工中，采用摆定前倾角和零度旋转角的加工方式下，非球头刀铣削无法进行切宽全局优化的问题，发明一种非球头刀铣削加工自由曲面的方法。

本发明的技术方案是：

一种非球头刀铣削加工自由曲面的方法，其特征是它包括以下步骤：

步骤1：依据最大弦高比值，对加工曲面进行网格划分，得到加工曲面包含的网格列表及每个网格包含的网格点；

步骤2：构造刀具有效切削曲面表达方式；

步骤3：基于加工曲面几何特性及刀具有效切削曲面表达方式，计算每个网格点处的非球头刀切宽张量系数矩阵，构建切宽张量场；

步骤4：搜索切宽张量场中包含的退化点，提取三分退化点；

步骤5：以三分退化点为起点搜索曲面内部边界点，构造曲面内部边界，得到曲面加工特征。

步骤6：在每一个曲面加工特征中构建数条曲线，保证曲线上每点的切向与该点处切宽最大的进给方向重合，选出平均切宽最大的曲线为该区域内的初始加工轨迹线；

步骤7：采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。

上述步骤1中所述的加工曲面网格划分具体方法如下：将加工曲面投影至零件底面B上，如图2所示，提取包围该投影区域的最小矩形T，设置T在两垂直边方向上的离散步距m和n，其中m、n的初始默认值为各自所对应边长的百分之一，以初始m和n对T进行离散，得到底面上的一系列网格BG，将网格BG沿底面的法向投影至加工曲面S上，得到曲面网格列表SG和网格点。并计算出所有两点间的最大弦高比，并将其与弦高比预设值进行比较，其中弦高比预设值可默认为N_XG，并可根据曲面加工精度进行调整，如果最大弦高比小于N_XG则调节对应m、n值并继续迭代计算，直到满足条件，即最大弦高比大于等于N_XG，m、n确定后，记录下对应的曲面网格列表SG。

上述步骤2中所述的构建抽象刀具曲面表达方式，具体是指，由于在实际加工中，切削宽度的度量与刀具和零件的几何信息相关，因此首先需要对刀具的曲面进行描述，对于APT刀具(通用自动编程刀具)而言，如图3所示，其刀具几何表达为：

其中a₁∈[0,1],a₂∈[0,1],λ∈[-π,+π]，并且在本方法中，采用成熟的摆定前倾角λ₀、旋转角ω为零的切削方式，如图4所示，对于任意刀触点，在其局部坐标系下，刀具的几何轮廓可以表达为：

在实际加工中，刀具可沿刀触点所在的加工曲面法平面上的任意方向切削，因此刀具的有效切削曲面ECS应该为刀具有效切削轮廓ECP绕Z方向旋转180度，扫掠出的回转轮廓面，其中刀具有效切削轮廓ECP为刀具前进方向法平面与刀具曲面的交线，在摆定前倾角和零度旋转角的加工方式中，根据广义刀具几何轮廓，在局部坐标系下，刀具有效切削轮廓ECP可表示为

从而绕Z方向旋转得到的刀具有效切削曲面ECS可表示为

上述步骤3所述的基于加工曲面几何特性及刀具有效切削曲面表达方式，构建切宽张量场方法如下：

设u¹和u²分别为待加工曲面的通过参数化得到的两个参数，切削宽度w的张量表示为：

w＝H_αβdu^αdu^β,α,β＝0,or1

其中切宽张量系数矩阵H_αβ可表示为：

其中P＝(A^TB)^-1G，A和B分别为加工曲面S和刀具有效切削曲面ECS的雅阁比(Jacobian)矩阵，D和W分别为S和ECS的曲面第二标准型系数矩阵，G为S的第一标准型系数矩阵，G^-1为G的反向矩阵。在切宽张量场中的每一个曲面离散点处，H_αβ的最大特征方向即为切宽最大的进给方向，当图5所示的曲面选用直径为6mm的平底刀，前倾角设为5度进行加工时，该曲面上各离散点的切宽最大进给方向可通过求解各点处H_αβ的最大特征方向获得，且在曲面范围内，切宽最大的进给方向呈现出区域化分布的趋势；

上述步骤4所述的搜索切宽张量场中包含的退化点具体实现方法如下：切宽张量场的退化点处满足各个进给方向上的切宽大小相同，因此在退化点处H_αβ具有两个相同大小的特征值，此时H_αβ满足如图6所示，对于网格SG_i，首先判断其是否包含退化点，设其四个顶点分别记为P_i1、P_i2、P_i3和P_i4，各点处(H₁₁-H₂₂)分别记为t_i1、t_i2、t_i3和t_i4,H₁₂分别记为s_i1、s_i2、s_i3和s_i4，若满足则SG_i内包含退化点，其中t_min和t_max分别为t_i1、t_i2、t_i3和t_i4中的最小值与最大值，s_min和s_max分别为s_i1、s_i2、s_i3和s_i4中的最小值与最大值。当确定SG_i内包含退化点之后，采用双线性插值方法求解去除落在[0,1]外的u、v解便得到SG_i包含的退化点DP_i，图6中所示的{P₁,P₂,P₃,P₄}为图5切宽张量场中包含的退化点；

上述步骤4所述的提取三分退化点的具体实现方法如下：计算各退化点处的和δ＝αλ-βγ，若δ<0，则该退化点为三分退化点，即有三条内部边界通过的退化点，如图7中的点P₁、P₂；

上述步骤5中所述的曲面内部边界构造方法如图8所示，其具体实现方法如下：

以三分退化点为起始点，搜索曲面内部边界点，构造曲面内部边界线，尽量确保曲面内部边界线上每一点的切向与该点处切宽最大的进给方向一致。对于三分退化点T_DP_i，首先在其切平面内以该点为圆心做半径为R的圆C_i(R通常小于1mm，可由使用者自行设定)，对C_i按角度步距Δω(Δω为一非常小的角度值)进行离散，得到离散点集合CP_i，将CP_i沿T_DP_i在曲面上的法向N_i投影至曲面上得到PCP_i，通过计算H_αβ的最大特征方向获取PCP_i中所有点的切宽最大进给方向向量集合MV_i，并将MV_i沿N_i投影回切平面中得到PMV_i，过T_DP_i与CP_i中所有点做向量得到可选方向集合PV_i，计算所有PV_i与其对应的PMV_i间的夹角记作ANG_i。对于任意三分退化点T_DP_i，在CP_i中可找到三个点使其对应的ANG近似为0，满足该条件的点称为种子点SP，将T_DP_i与SP_ij连接成线段并沿N_i方向投影至曲面上便可以得到一条曲面内部边界曲线段CS_ij，SP_ij在曲面上的投影点PSP_ij为曲面上的内部边界点，因此对于一个三分退化点，其周围存在三个曲面内部边界点，即每一个三分退化点仅能被三条内部边界穿过。当从T_DP_i进行曲面内部边界构造时，首先检查该三分退化点是否已被内部边界经过，若已由三条内部边界经过该点，则说明无需对该点进行边界构造，若小于三条，则按上述方法求出T_DP_i临近的三个内部边界点，去掉落在穿过T_DP_i内部边界上的边界点，利用余下的点进行搜索构造下一段曲线段。

对于曲面内部边界点PSP_ij，采取同样的方法构造半径为R的圆C_ij并对其按Δω进行离散得到点集CP_ij，将CP_ij沿CP_ij处的曲面法向N_ij投影至曲面上得到点集PCP_ij，通过计算H_αβ的最大特征方向获取PCP_ij中所有点的切宽最大进给方向向量集合MC_ij，并将MV_ij沿N_ij投影回切平面中得到PMV_ij，过PSP_ij与CP_ij中所有点做向量得到可选方向集合PV_ij，计算所有PV_ij与其对应的PMV_ij间的夹角记作ANG_ij。此时CP_ij中可搜索得到两个种子点SP_ij1和SP_ij2使其对应的ANG近似为0，分别过SP_ij1和SP_ij2与PSP_ij做向量TV_ij1和TV_ij2，将上一段曲面内部边界曲线段沿N_ij投影至PSP_ij的切平面中得到向量TS_ij，分别计算TV_ij1和TV_ij2与TS_ij的夹角θ_ij1与θ_ij2，提取两个夹角中较大的那个夹角对应的种子点记作SP_ij，将IP_k与PSP_ij连接成线段并沿N_ij方向投影至曲面上得到内部边界曲线段并添加到当前的内部边界曲线上，用SP_ij在曲面上的投影点替换PSP_ij继续进行边界的构造。

曲面内部边界点搜索停止的条件为：当前曲面内部边界点PSP_ij所做的圆C_ij有部分落在曲面外侧，此时过PSP_ij做上一曲面内部边界曲线段的切线向量，将该向量沿N_ij投影至曲面上，提取该投影曲线段落在曲面内部的部分，将其作为这条曲面内部边界线的最后一段；②当所得到的曲面内部边界曲线段穿过某一三分退化点，或两者间最小距离小于某一给定的微小长度值，此时将该三分退化点看做是下一个曲面内部边界点，提取当前曲面内部边界点与该三分退化点间的最短曲线段作为这条曲面内部边界线的最后一段。

当所有三分退化点处都有三条曲面内部边界线经过则完成整个曲面区域的划分，每一个曲面区域称为一个曲面加工特征SMF。

上述步骤6中所述的各个曲面加工特征中的初始加工轨迹线构造具体实现方法如下：

如图9中流程所示，首先采用人工编辑的方法在各个曲面加工特征对应的曲面区域内创建一系列点，将每一个人工创建的点P_i以步骤4中的方法找到其对应的两个种子点SP_i1和SP_i2，过P_i分别与SP_i1和SP_i2做直线段，并沿P_i处的曲面法向N_i投影至曲面上得到曲线段CS_i1和CS_i2，将CS_i1和CS_i2保存为到曲线段集合CC_i中，相同的方式得到两个种子点在曲面上对应的投影点PSP_i1和PSP_i2，选择PSP_i1为起点按照步骤4中PSP_ij的搜索方法构造曲线段并添加到CCi中，构造终止条件为：①若当前搜索得到的曲线段穿过PSP_i2或者两者间的最小距离小于某一给定的微小长度值时，则以提取PSP_ij与PSP_i2间的曲线段作为CC_i中的最后一段，从而形成一条封闭的曲线段；②若在当前PSP_ij的切平面内构造的半径为R的圆有部分落在该曲面加工特征外时，则过该点做上一曲线段的切向量，将该切向量沿该点处的曲面法向投影至曲面上，取出该投影线段落在该曲面加工特征内部的部分添加到CC_i中，以SP_i2为起点继续进行搜索，若以SP_i2的起点的搜索已完成，则过P_i点的曲线构造完成。

对于曲面加工特征SMF_s，按照上述方法在其对应的曲面区域内创建了数条曲线，对于其中的任意一条，通过以下方法计算其平均切宽CW_a：

其中α,β＝0,or 1,u¹和u²分别为待加工曲面的通过参数化得到的两个参数，τ为曲线的弧长参数。

提取出CW_a最大的那条曲线作为SMF_s的初始刀具轨迹线。

本发明的有益效果是：

(1)基于特征的加工方法实现了曲面的分区加工，使刀轨能够遵循最大切宽进给方向在曲面上的区域化分布规律，尽可能使每一个刀触点处的进给方向与切宽最大的进给方向重合，从而缩短刀轨长度，实现加工效率的提高；

(2)用张量来度量非球头刀加工曲面的切宽，从理论上突破了以切宽最大化为目标的曲面加工区域划分方法，能够自动生成准确的曲面内部边界，消除人为因素对分区的影响，实现曲面加工特征的自动提取；

(3)该方法在选择每个曲面加工特征的初始轨迹线时，采用连续积分的方法，能够得到各条待选曲线的准确平均切宽，保证了初始轨迹线的优化性；

(4)该方法尤其适用于曲率变化剧烈的复杂曲面加工，加工效率提升显著。

附图说明

图1为本发明的非球头刀铣削加工自由曲面方法流程图。

图2为本发明加工曲面网格划分方法示意图，其中m、n分别是Y和X方向上的离散步距，L₁为步距选取对应的弦高值，L₂为步距选取对应的弦长值，B为底面，T为最小包围矩形，S为曲面，BG为投影底面上的网格，SG为曲面网格列表；

图3为标准APT刀具示意图。

图4为本发明的平底刀铣削前倾角示意图，其中X、Y、Z为刀触点CCpoint的局部坐标系，X轴为进给方向，Z为刀触点处曲面外法向，T为刀轴方向，λ为前倾角示意，ω为旋转角示意，λ₀为实际前倾角，ECP为当前刀具切削部分与刀具前进方向的法平面的交线，箭头所示为ECP绕Z轴旋转180度，ECS为旋转后得到的刀具有效切削曲面。

图5为本发明加工曲面S上各点利用切宽张量系数矩阵H求出的切削宽度最大进给方向分布图，其中每条线段为所在点处切宽最大进给方向。

图6为本发明加工曲面上的网格SG_i，其中P_i1、P_i2、P_i3和P_i4为它的四个顶点，u和v为加工曲面的两个参数。

图7为图4中加工曲面S上的退化点分布，其中P₁、P₂为三分退化点。

图8为以三分退化点为起点的曲面内部边界构造算法流程图。

图9为各曲面加工特征对应曲面区域内的待选初始轨迹线构造算法流程图。

图10为平底刀铣削自由曲面加工方法的应用案例，图a中T_DP₁和T_DP₂为两个三分退化点，IB₁₁、IB₁₂和IB₁₃为从T_DP₁出发搜索得到的三条内部边界，IB₂₁、IB₂₂和IB₂₃为从T_DP₂出发搜索得到的三条内部边界，其中IB₁₃和IB₂₁重合，算作一条内部边界。该曲面被以上5条内部边界线分成了4个加工特征，分别为SMF₁、SMF₂、SMF₃和SMF₄；图b中所示IPC₁、IPC₂、IPC₃和IPC₄分别为这4个特征选出来的初始轨迹线，采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨，如图c所示。

具体实施方式

下面是结合附图和具体实施例对本发明的技术方案进行详细说明。

如图1-10所示。

一种非球头刀铣削加工自由曲面的方法，它主要包括以下步骤：

步骤1：依据最大弦高比值，对加工曲面进行网格划分，得到加工曲面包含的网格列表及每个网格包含的网格点；

步骤2：构造刀具有效切削曲面表达方式；

步骤3：基于加工曲面几何特性及刀具有效切削曲面表达方式，计算每个网格点处的非球头刀切宽张量系数矩阵，构建切宽张量场；

步骤4：搜索切宽张量场中包含的退化点，提取三分退化点；

步骤5：以三分退化点为起点搜索曲面内部边界点，构造曲面内部边界，得到曲面加工特征。

步骤6：在每一个曲面加工特征中构建数条曲线，保证曲线上每点的切向与该点处切宽最大的进给方向重合，选出平均切宽最大的曲线为该区域内的初始加工轨迹线；

步骤7：采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。

详述如下：

上述步骤1中所述的加工曲面网格划分具体方法如下：将加工曲面投影至零件底面B上，如图2所示，提取包围该投影区域的最小矩形T，设置T在两垂直边方向上的离散步距m和n，其中m、n的初始默认值为各自所对应边长的百分之一，以初始m和n对T进行离散，得到底面上的一系列网格BG，将网格BG沿底面的法向投影至加工曲面S上，得到曲面网格列表SG和网格点。并计算出所有两点间的最大弦高比，并将其与弦高比预设值进行比较，其中弦高比预设值可默认为N_XG，并可根据曲面加工精度进行调整，如果最大弦高比小于N_XG则调节对应m、n值并继续迭代计算，直到满足条件，即最大弦高比大于等于N_XG，将m、n确定后，记录下对应的曲面网格列表SG。

步骤2中所述的构建抽象刀具曲面表达方式，具体是指，由于在实际加工中，切削宽度的度量与刀具和零件的几何信息相关，因此首先需要对刀具的曲面进行描述，对于APT刀具(通用自动编程刀具)而言，如图3所示，其刀具几何表达为：

其中a₁∈[0,1],a₂∈[0,1],λ∈[-π,+π]，并且在本方法中，采用成熟的摆定前倾角λ₀、旋转角ω为零的切削方式，如图4所示，对于任意刀触点，在其局部坐标系下，刀具的几何轮廓可以表达为：

在实际加工中，刀具可沿刀触点所在的加工曲面法平面上的任意方向切削，因此刀具的有效切削曲面ECS应该为刀具有效切削轮廓ECP绕Z方向旋转180度，扫掠出的回转轮廓面，其中刀具有效切削轮廓ECP为刀具前进方向法平面与刀具曲面的交线，在摆定前倾角和零度旋转角的加工方式中，根据广义刀具几何轮廓，在局部坐标系下，刀具有效切削轮廓ECP可表示为：

从而绕Z方向旋转得到的刀具有效切削曲面ECS可表示为：

步骤3所述的基于加工曲面几何特性及刀具有效切削曲面表达方式，构建切宽张量场方法如下：

设u¹和u²分别为待加工曲面的通过参数化得到的两个参数，切削宽度w的张量表示为：

w＝H_αβdu^αdu^β,α,β＝0,or1

其中切宽张量系数矩阵H_αβ可表示为：

其中P＝(A^TB)^-1G，A和B分别为加工曲面S和刀具有效切削曲面ECS的雅阁比(Jacobian)矩阵，D和W分别为S和ECS的曲面第二标准型系数矩阵，G为S的第一标准型系数矩阵，G^-1为G的反向矩阵。对于摆定前倾角λ₀，旋转角为零的平底刀而言，β₁＝β₂＝0，R₂＝0，R₁＝R，α＝90°其刀具轮廓曲面为：

其中a₁,a₂为1。其刀具有效切削轮廓为：

其中a与相对应。刀具有效切削曲面为：

此时

在切宽张量场中的每一个曲面离散点处，H_αβ的最大特征方向即为切宽最大的进给方向，当图5所示的曲面选用直径为6mm的平底刀，前倾角设为5度进行加工时，该曲面上各离散点的切宽最大进给方向可通过求解各点处H_αβ的最大特征方向获得，且在曲面范围内，切宽最大的进给方向呈现出区域化分布的趋势；

步骤4所述的搜索切宽张量场中包含的退化点具体实现方法如下：切宽张量场的退化点处满足各个进给方向上的切宽大小相同，因此在退化点处H_αβ具有两个相同大小的特征值，此时H_αβ满足如图6所示，对于网格SG_i，首先判断其是否包含退化点，设其四个顶点分别记为P_i1、P_i2、P_i3和P_i4，各点处(H₁₁-H₂₂)分别记为t_i1、t_i2、t_i3和t_i4,H₁₂分别记为s_i1、s_i2、s_i3和s_i4，若满足则SG_i内包含退化点，其中t_min和t_max分别为t_i1、t_i2、t_i3和t_i4中的最小值与最大值，s_min和s_max分别为s_i1、s_i2、s_i3和s_i4中的最小值与最大值。当确定SG_i内包含退化点之后，采用双线性插值方法求解去除落在[0,1]外的u、v解便得到SG_i包含的退化点DP_i，图7中所示的{P₁,P₂,P₃,P₄}为图5切宽张量场中包含的退化点。

步骤4所述的提取三分退化点的具体实现方法如下：计算各退化点处的和δ＝αλ-βγ，若δ<0，则该退化点为三分退化点，即有三条内部边界通过的退化点，如图7中的点P₁、P₂；

步骤5中所述的曲面内部边界构造方法如图8所示，其具体实现方法如下：

以三分退化点为起始点，搜索曲面内部边界点，构造曲面内部边界线，尽量确保曲面内部边界线上每一点的切向与该点处切宽最大的进给方向一致。对于三分退化点T_DP_i，首先在其切平面内以该点为圆心做半径为R的圆C_i(R通常小于1mm，可由使用者自行设定)，对C_i按角度步距Δω(Δω为一非常小的角度值)进行离散，得到离散点集合CP_i，将CP_i沿T_DP_i在曲面上的法向N_i投影至曲面上得到PCP_i，通过计算H_αβ的最大特征方向获取PCP_i中所有点的切宽最大进给方向向量集合MV_i，并将MV_i沿N_i投影回切平面中得到PMV_i，过T_DP_i与CP_i中所有点做向量得到可选方向集合PV_i，计算所有PV_i与其对应的PMV_i间的夹角记作ANG_i。对于任意三分退化点T_DP_i，在CP_i中可找到三个点使其对应的ANG近似为0，满足该条件的点称为种子点SP，将T_DP_i与SP_ij连接成线段并沿N_i方向投影至曲面上便可以得到一条曲面内部边界曲线段CS_ij，SP_ij在曲面上的投影点PSP_ij为曲面上的内部边界点，因此对于一个三分退化点，其周围存在三个曲面内部边界点，即每一个三分退化点仅能被三条内部边界穿过。当从T_DP_i进行曲面内部边界构造时，首先检查该三分退化点是否已被内部边界经过，若已由三条内部边界经过该点，则说明无需对该点进行边界构造，若小于三条，则按上述方法求出T_DP_i临近的三个内部边界点，去掉落在穿过T_DP_i内部边界上的边界点，利用余下的点进行搜索构造下一段曲线段。

对于曲面内部边界点PSP_ij，采取同样的方法构造半径为R的圆C_ij并对其按Δω进行离散得到点集CP_ij，将CP_ij沿CP_ij处的曲面法向N_ij投影至曲面上得到点集PCP_ij，通过计算H_αβ的最大特征方向获取PCP_ij中所有点的切宽最大进给方向向量集合MV_ij，并将MV_ij沿N_ij投影回切平面中得到PMV_ij，过PSP_ij与CP_ij中所有点做向量得到可选方向集合PV_ij，计算所有PV_ij与其对应的PMV_ij间的夹角记作ANG_ij。此时CP_ij中可搜索得到两个种子点SP_ij1和SP_ij2使其对应的ANG近似为0，分别过SP_ij1和SP_ij2与PSP_ij做向量TV_ij1和TV_ij2，将上一段曲面内部边界曲线段沿N_ij投影至PSP_ij的切平面中得到向量TS_ij，分别计算TV_ij1和TV_ij2与TS_ij的夹角θ_ij1与θ_ij2，提取两个夹角中较大的那个夹角对应的种子点记作SP_ij，将IP_k与PSP_ij连接成线段并沿N_ij方向投影至曲面上得到内部边界曲线段并添加到当前的内部边界曲线上，用SP_ij在曲面上的投影点替换PSP_ij继续进行边界的构造。

曲面内部边界点搜索停止的条件为：当前曲面内部边界点PSP_ij所做的圆C_ij有部分落在曲面外侧，此时过PSP_ij做上一曲面内部边界曲线段的切线向量，将该向量沿N_ij投影至曲面上，提取该投影曲线段落在曲面内部的部分，将其作为这条曲面内部边界线的最后一段；当所得到的曲面内部边界曲线段穿过某一三分退化点，或两者间最小距离小于某一给定的微小长度值，此时将该三分退化点看做是下一个曲面内部边界点，提取当前曲面内部边界点与该三分退化点间的最短曲线段作为这条曲面内部边界线的最后一段。

当所有三分退化点处都有三条曲面内部边界线经过则完成整个曲面区域的划分，每一个曲面区域称为一个曲面加工特征SMF。

步骤6中所述的各个曲面加工特征中的初始加工轨迹线构造具体实现方法如下：

如图9中流程所示，首先采用人工编辑的方法在各个曲面加工特征对应的曲面区域内创建一系列点，将每一个人工创建的点P_i以步骤4中的方法找到其对应的两个种子点SP_i1和SP_i2，过P_i分别与SP_i1和SP_i2做直线段，并沿P_i处的曲面法向N_i投影至曲面上得到曲线段CS_i1和CS_i2，将CS_i1和CS_i2保存为到曲线段集合CC_i中，相同的方式得到两个种子点在曲面上对应的投影点PSP_i1和PSP_i2，选择PSP_i1为起点按照步骤4中PSP_ij的搜索方法构造曲线段并添加到CC_i中，构造终止条件为：①若当前搜索得到的曲线段穿过PSP_i2或者两者间的最小距离小于某一给定的微小长度值时，则以提取PSP_ij与PSP_i2间的曲线段作为CC_i中的最后一段，从而形成一条封闭的曲线段；②若在当前PSP_ij的切平面内构造的半径为R的圆有部分落在该曲面加工特征外时，则过该点做上一曲线段的切向量，将该切向量沿该点处的曲面法向投影至曲面上，取出该投影线段落在该曲面加工特征内部的部分添加到CC_i中，以SP_i2为起点继续进行搜索，若以SP_i2的起点的搜索已完成，则过P_i点的曲线构造完成。

对于曲面加工特征SMF_s，按照上述方法在其对应的曲面区域内创建了数条曲线，对于其中的任意一条，通过以下方法计算其平均切宽CW_a：

其中α,β＝0,or 1,u¹和u²分别为将曲面参数化后得到的两参数，τ为曲线的弧长参数。

提取出CW_a最大的那条曲线作为SMF_s的初始刀具轨迹线。

根据以上所述方法对图10中的曲面进行加工轨迹生成，该曲面精加工所用平底刀的直径为6mm，前倾角设为8度，表面最大残高为0.05mm。图10 a中T_DP₁和T_DP₂为该曲面上的两个三分退化点，IB₁₁、IB₁₂和IB₁₃为从T_DP₁出发搜索得到的三条内部边界，IB₂₁、IB₂₂和IB₂₃为从T_DP₂出发搜索得到的三条内部边界，其中IB₁₃与IB₂₁重合，算为一条边界。该曲面被这5条内部边界线分成了4个加工特征，分别为SMF₁、SMF₂、SMF₃和SMF₄；图10 b中所示IPC₁、IPC₂、IPC₃和IPC₄分别为这4个特征选出来的初始轨迹线，采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨，如图10c所示，经测算，利用本方法生成的轨迹线总长度为1366.5mm，而采用传统的方法，将整个曲面作为一个加工区域，采用Iso-parametric轨迹生成算法，以曲面最长的边界起始轨迹线，得到的轨迹总长度为1607.6mm，利用本方法生成的轨迹在总长度上缩短了14.99％。

本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。

Claims (7)

1.一种非球头刀铣削加工自由曲面的方法，其特征是它包括以下步骤：

步骤1：依据最大弦高比值，对加工曲面S进行网格划分，得到加工曲面包含的网格列表及每个网格包含的网格点；

步骤2：构造刀具有效切削曲面表达方式；由于在实际加工中，切削宽度的度量与刀具和零件的几何信息相关，因此首先需要对刀具的曲面进行描述，对于APT刀具而言，其刀具几何表达为：

其中a₁和a₂分别为刀具径向和轴向的长度系数，且a₁∈[0,1],a₂∈[0,1]，R₁为刀底半径，L为刀具柱体的轴向高度，R₂为底角半径，β₁为刀底面与刀具径向平面夹角，β₂为刀具柱面与刀具轴向夹角，α为底角圆弧位置角，且 为刀具径向平面内的旋转位置角，且采用摆定前倾角λ₀、旋转角ω为零的切削方式，对于任意刀触点，在其局部坐标系下，刀具的几何轮廓可以表达为：

由于刀具沿刀触点所在的加工曲面法平面上的任意方向切削，因此刀具的有效切削曲面ECS应该为刀具有效切削轮廓ECP绕Z轴方向旋转180度，扫掠出的回转轮廓面，其中刀具有效切削轮廓ECP为刀具前进方向法平面与刀具曲面的交线，在摆定前倾角和零度旋转角的加工方式中，根据广义刀具几何轮廓，在局部坐标系下，刀具有效切削轮廓ECP可表示为：

从而绕Z轴方向旋转得到的刀具有效切削曲面ECS可表示为：

步骤3：基于加工曲面几何特性及刀具有效切削曲面表达方式，计算每个网格点处的非球头刀切宽张量系数矩阵，构建切宽张量场；

步骤4：搜索切宽张量场中包含的退化点，提取三分退化点；

步骤5：以三分退化点为起点搜索曲面内部边界点，构造曲面内部边界，得到曲面加工特征；

步骤6：在每一个曲面加工特征中构建数条曲线，保证曲线上每点的切向与该点处切宽最大的进给方向重合，选出平均切宽最大的曲线为该区域内的初始加工轨迹线；

步骤7：采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。

2.根据权利要求1所述的方法，其特征是所述的加工曲面网格划分方法为：将加工曲面投影至零件底面上，提取包围该投影区域的最小矩形T，设置T在两垂直边方向上的离散步距m和n，其中m、n的初始默认值为各自所对应边长的百分之一，以初始m和n对T进行离散，得到底面上的一系列网格BG，将网格BG沿底面的法向投影至加工曲面S上，得到曲面网格列表SG和网格点，并计算出所有两点间的最大弦高比，将其与弦高比预设值进行比较，其中弦高比预设值可默认为N_XG，并可根据曲面加工精度进行调整，如果最大弦高比小于N_XG则调节对应m、n值并继续迭代计算，直到满足条件最大弦高比大于等于N_XG，将m、n确定后，记录下对应的曲面网格列表SG。

3.根据权利要求1所述的方法，其特征是所述的构建切宽张量场的方法为：

设u¹和u²分别为待加工曲面的通过参数化得到的两个参数，切削宽度w的张量表示为：

w＝H_αβdu^αdu^β,α,β＝0,or 1

其中切宽张量系数矩阵H_αβ可表示为：

$H_{\mathrm{\&alpha;}\mathrm{\&beta;}}=\frac{1}{2}(P^{T}WP-D)G^{-1}$

其中P＝(A^TB)^-1G，A和B分别为加工曲面S和刀具有效切削曲面ECS的雅阁比(Jacobian)矩阵，D和W分别为加工曲面S和刀具有效切削曲面ECS的曲面第二标准型系数矩阵，G为S的第一标准型系数矩阵，G^-1为G的反向矩阵；在切宽张量场中的每一个曲面离散点处，H_αβ的最大特征方向即为切宽最大的进给方向。

4.根据权利要求1所述的方法，其特征是所述的搜索切宽张量场中包含的退化点的方法为：切宽张量场的退化点处满足各个进给方向上的切宽大小相同，因此在退化点处H_αβ具有两个相同大小的特征值，此时H_αβ满足式中H₁₁、H₂₂、H₁₂、H₂₁为H_αβ二阶矩阵表达的对应项

5.根据权利要求1所述的方法，其特征是所述的提取三分退化点的方法为：计算各退化点处的和δ＝αλ-βγ，若δ<0，则该退化点为三分退化点，即有三条内部边界通过的退化点，式中u和v分别为待加工曲面的通过参数化得到的两个参数。

6.根据权利要求1所述的方法，其特征是所述的曲面内部边界构造方法如下：

以三分退化点为起始点，搜索曲面内部边界点，构造曲面内部边界线，确保曲面内部边界线上每一点的切向与该点处切宽最大的进给方向一致；对于三分退化点T_DP_i，首先在其切平面内以该点为圆心做半径为R的圆C_i，对C_i按角度步距Δω进行离散，得到离散点集合CP_i，将CP_i沿T_DP_i在曲面上的法向N_i投影至曲面上得到PCP_i，通过计算H_αβ的最大特征方向获取PCP_i中所有点的切宽最大进给方向向量集合MV_i，并将MV_i沿N_i投影回切平面中得到PMV_i，过T_DP_i与CP_i中所有点做向量得到可选方向集合PV_i，计算所有PV_i与其对应的PMV_i间的夹角记作ANG_i；对于任意三分退化点T_DP_i，在CP_i中可找到三个点使其对应的ANG近似为0，满足该条件的点称为种子点SP，将T_DP_i与SP_ij连接成线段并沿N_i方向投影至曲面上便可以得到一条曲面内部边界曲线段CS_ij，SP_ij在曲面上的投影点PSP_ij为曲面上的内部边界点，因此对于一个三分退化点，其周围存在三个曲面内部边界点，即每一个三分退化点仅能被三条内部边界穿过；当从T_DP_i进行曲面内部边界构造时，首先检查该三分退化点是否已被内部边界经过，若已由三条内部边界经过该点，则说明无需对该点进行边界构造，若小于三条，则按上述方法求出T_DP_i临近的三个内部边界点，去掉落在穿过T_DP_i内部边界上的边界点，利用余下的点进行搜索构造下一段曲线段；

对于曲面内部边界点PSP_ij，采取同样的方法构造半径为R的圆C_ij并对其按Δω进行离散得到点集CP_ij，将CP_ij沿CP_ij处的曲面法向N_ij投影至曲面上得到点集PCP_ij，通过计算H_αβ的最大特征方向获取PCP_ij中所有点的切宽最大进给方向向量集合MV_ij，并将MV_ij沿N_ij投影回切平面中得到PMV_ij，过PSP_ij与CP_ij中所有点做向量得到可选方向集合PV_ij，计算所有PV_ij与其对应的PMV_ij间的夹角记作ANG_ij；此时CP_ij中可搜索得到两个种子点SP_ij1和SP_ij2使其对应的ANG近似为0，分别过SP_ij1和SP_ij2与PSP_ij做向量TV_ij1和TV_ij2，将上一段曲面内部边界曲线段沿N_ij投影至PSP_ij的切平面中得到向量TS_ij，分别计算TV_ij1和TV_ij2与TS_ij的夹角θ_ij1与θ_ij2，提取两个夹角中较大的那个夹角对应的种子点记作SP_ij，将IP_k与PSP_ij连接成线段并沿N_ij方向投影至曲面上得到内部边界曲线段并添加到当前的内部边界曲线上，用SP_ij在曲面上的投影点替换PSP_ij继续进行边界的构造；

曲面内部边界点搜索停止的条件为：当前曲面内部边界点PSP_ij所做的圆C_ij有部分落在曲面外侧，此时过PSP_ij做上一曲面内部边界曲线段的切线向量，将该向量沿N_ij投影至曲面上，提取该投影曲线段落在曲面内部的部分，将其作为这条曲面内部边界线的最后一段；当所得到的曲面内部边界曲线段穿过某一三分退化点，或两者间最小距离小于某一给定的微小长度值，此时将该三分退化点看做是下一个曲面内部边界点，提取当前曲面内部边界点与该三分退化点间的最短曲线段作为这条曲面内部边界线的最后一段；

当所有三分退化点处都有三条曲面内部边界线经过则完成整个曲面区域的划分，每一个曲面区域称为一个曲面加工特征SMF。

7.根据权利要求1所述的方法，其特征是所述的各个曲面加工特征中的初始加工轨迹线构造方法如下：

首先采用人工编辑的方法在各个曲面加工特征对应的曲面区域内创建一系列点，将每一个人工创建的点P_i以步骤4中的方法找到其对应的两个种子点SP_i1和SP_i2，过P_i分别与SP_i1和SP_i2做直线段，并沿P_i处的曲面法向N_i投影至曲面上得到曲线段CS_i1和CS_i2，将CS_i1和CS_i2保存为到曲线段集合CC_i中，相同的方式得到两个种子点在曲面上对应的投影点PSP_i1和PSP_i2，选择PSP_i1为起点按照步骤4中PSP_ij的搜索方法构造曲线段并添加到CC_i中，构造终止条件为：①若当前搜索得到的曲线段穿过PGP_i2或者两者间的最小距离小于某一给定的微小长度值时，则以提取PSP_ij与PSP_i2间的曲线段作为CC_i中的最后一段，从而形成一条封闭的曲线段；②若在当前PSP_ij的切平面内构造的半径为R的圆有部分落在该曲面加工特征外时，则过该点做上一曲线段的切向量，将该切向量沿该点处的曲面法向投影至曲面上，取出该投影线段落在该曲面加工特征内部的部分添加到CC_i中，以SP_i2为起点继续进行搜索，若以SP_i2的起点的搜索已完成，则过P_i点的曲线构造完成；

对于曲面加工特征SMF_s，按照上述方法在其对应的曲面区域内创建了数条曲线，对于其中的任意一条，通过以下方法计算其平均切宽CW_a：

其中α,β＝0,or 1,u¹和u²分别为待加工曲面的通过参数化得到的两个参数，τ为曲线的弧长参数；

提取出CW_a最大的那条曲线作为SMF_s的初始刀具轨迹线。