CN109358568B - 基于矢量场的曲面分区加工轨迹拓扑形状设计方法 - Google Patents

Abstract

本发明基于矢量场的曲面分区加工轨迹拓扑形状设计方法属于复杂曲面零件高精高效铣削加工技术领域，涉及一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法。该方法建立弦高误差与残高误差协同约束的走刀矢量场，构建归一化两类误差双目标优化模型，通过评价函数求解各点最优走刀方向，建立空间矢量场；计算投影后平面矢量场散度与旋度，根据其是否为零实现加工区域初分割；分析不同矢量场特征，拟合空间矢量场流线，计算机床旋转轴沿流线进给运动的运动学参数，根据其是否突变完成加工区域细分；结合走刀方式实现曲面不同加工区域加工轨迹拓扑形状设计。该方法适用于复杂曲面分区加工，可减小加工误差，提高进给运动平稳性。

Description

基于矢量场的曲面分区加工轨迹拓扑形状设计方法

技术领域

本发明属于复杂曲面零件高精高效铣削加工技术领域，涉及一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法。

背景技术

当前，航空航天、能源动力、汽车船舶等领域的发展对复杂曲面零件的加工精度提出了极高的要求，传统采用统一工艺全域加工的方法适用度降低，易导致误差分布不均及局部超差，难以适用于复杂曲面零件高精度加工。为满足此类零件的加工要求，常采用分区加工方法，将复杂曲面划分成多个加工区域，对不同加工区域分别进行适宜的加工轨迹及加工工艺规划。传统分区方法多围绕曲面曲率、法向量等制定分区准则，其分区结果对于避免加工干涉与提高加工效率等方面具有一定的指导意义。然而，现有分区准则仅针对曲面几何特征，缺乏对加工走刀方式与加工过程中进给运动平稳性的综合考虑，且加工区域分割与轨迹规划相互独立，不同加工区域的加工轨迹拓扑形状设计缺乏理论依据，易诱发局部误差超差及显著切削振动，严重制约曲面加工质量的全面提高。由此，迫切需要研究一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法，以满足复杂曲面零件分区加工的高精度加工要求。

郝小忠等人的专利“考虑误差区域性分布的数控加工刀轨分区映射方法”，公开号CN107145126A。该专利根据工件理论面生成理论刀轨，将工件实际面与理论面对比得到误差分布规律，根据误差分布规律对实际面进行区域划分，针对不同误差类型的区域选择相应的刀轨映射方法将理论刀轨映射到工件实际面上完成刀轨映射过程。然而，该方法的分区结果易受测量误差的影响，且分区过程中实际面与理论面的对比方法缺乏理论支撑，因此该分区方法具有较大局限性。Tuong等人的文献“A practical approach forpartitioning free-form surfaces”，International Journal of Computer IntegratedManufacturing，2010，23(11)，992-1001，提出了一种基于曲面曲率的复杂曲面分区方法。该方法通过计算曲面各点的高斯曲率和平均曲率，将加工区域划分为凹区、凸区及鞍区，并对不同区域规划不同加工工艺，实现了加工效率的大幅度提高。然而，该方法未考虑加工过程中进给运动平稳性，加工轨迹拓扑形状设计缺乏理论依据，具有一定的局限性。

发明内容

本发明针对现有技术缺陷，发明了一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法。该方法针对复杂曲面分区加工过程中，加工区域分割与轨迹规划缺乏对走刀方式与进给运动平稳性综合考虑，造成的局部误差超差及显著切削振动等问题，通过建立弦高误差与残高误差协同约束的走刀矢量场，实现基于走刀矢量分布的加工区域初分割及确保进给运动平稳的加工区域细分，根据分区结果，最终完成复杂曲面零件分区加工轨迹拓扑形状设计。由此有效提高了曲面分区加工的加工精度及进给运动的平稳性，为复杂曲面高精高效加工提供技术支撑。

本发明的技术方案是一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法，其特征在于，该方法建立弦高误差与残高误差协同约束的走刀矢量场，经归一化构建两类误差的双目标优化模型，通过构造评价函数求解各点最优走刀方向，建立空间矢量场；计算投影后平面矢量场散度与旋度，根据其是否为零对矢量场分类，实现加工区域初分割；分析不同矢量场特征，拟合空间矢量场流线，计算机床旋转轴沿流线进给运动的运动学参数，根据其是否突变完成加工区域细分；结合走刀方式，实现曲面不同加工区域加工轨迹拓扑形状设计；方法的具体步骤如下：

步骤1，建立弦高误差与残高误差协同约束的走刀矢量场

取空间直角坐标系O-xyz中的曲面加工区域为S(u,v)＝(x(u,v),y(u,v),z(u,v))，P为曲面加工轨迹上任意一点，在行距d与步长f确定时，弦高误差e计算为：

其中，k_f为走刀方向法曲率。

取变量e^*计算为：

e^*＝k_f ²(2)

分别求解弦高误差e与变量e^*的单调性知，两者呈正相关关系，为降低后续计算难度，使用变量e^*表示弦高误差。

残高误差h计算为：

其中，k_d为行距方向法曲率，r为刀具有效半径。

由微分几何知，行距方向与走刀方向互相垂直，满足式(4)：

k_d+k_f＝2H(4)

其中，H为平均曲率，对于曲面任意一点，H为定值。

联立式(3)-(4)，残高误差表示为：

由式(2)与式(5)知，步长及行距确定时，弦高误差与残高误差仅由走刀方向法曲率决定。为求解各点最优走刀方向，以两类误差为目标函数构造双目标优化模型。

对两类误差归一化处理为：

其中，e_n和h_n分别表示归一化的弦高误差和残高误差，

和

表示弦高误差的最小值和最大值，h_min和h_max表示残高误差的最小值和最大值，其求解如下：

以k₁和k₂表示曲面主曲率(k₁>k₂)，其计算为：

其中，K为高斯曲率，对于曲面任意一点，K为定值。

由式(2)与式(5)得，残高误差的最值为：

弦高误差的最大值为：

弦高误差的最小值为：

根据式(6)-(10)得，归一化的残高误差h_n为：

归一化的弦高误差e_n为：

为方便描述，后文出现的弦高误差与残高误差均为归一化的弦高误差与残高误差，以两类误差为目标函数建立双目标优化模型为：

运用线性加权法构造评价函数Q为：

其中，α₁和α₂分别为弦高误差和残高误差在评价函数中所占权重，其计算如下：

其中，h_n1为弦高误差取最小值时对应的残高误差，e_n1为残高误差取最小值时对应的弦高误差，h_nmin为残高误差最小值，e_nmin为弦高误差最小值。

根据主曲率取值范围，计算双目标优化模型最优解如下：

(1)当k₂≥0时，取k_f＝k₂，e_n和h_n均取得最小值，因而，k_f＝k₂为双目标优化模型最优解，由欧拉定理知：

k_f＝k₁cos²θ+k₂sin²θ(16)

其中，θ为走刀方向与k₁对应主方向的夹角，此时，

(2)当-k₁≤k₂<0时，评价函数Q为：

取

Q取得最小值，因而，

为双目标优化模型最优解，此时，

(3)当0≤k₁<-k₂时，评价函数Q为：

取

Q取得最小值，因而，

为双目标优化模型最优解，此时，

(4)当k₁<0时，评价函数Q为：

取

Q取得最小值，因而，

为双目标优化模型最优解，此时，

以(du:dv)表示曲面任意切方向，沿该方向的曲面第一基本形式I和第二基本形式II分别为：

其中，E、F、G为曲面第一类基本量，L、M、N为曲面第二类基本量，表达式如下：

其中，S_u、S_v为曲面S(u,v)的一阶偏导数，S_uu、S_uv、S_vv为曲面S(u,v)的二阶偏导数，n为曲面S(u,v)的法向量。

以(du¹:dv¹)表示主曲率k₁对应主方向，其计算为：

以(du^*:dv^*)表示最优走刀方向，其计算为：

其中，a、b、c分别表示为：

联立式(22)-(24)，求解获得曲面各点最优走刀方向，沿此方向建立单位走刀矢量，从而形成弦高误差与残高误差协同约束的走刀矢量场。

步骤2，基于走刀矢量分布的加工区域初分割

经步骤1获得的空间走刀矢量场在空间直角坐标系O-xyz中可表示为F(γ)＝(g₁(x,y,z),g₂(x,y,z),g₃(x,y,z))，其中γ为曲面各点沿最优走刀方向的单位走刀矢量。将该矢量场向xOy平面投影，沿各点投影后的走刀方向建立单位走刀矢量γ^*，获得平面矢量场A(γ^*)＝(f₁(x,y),f₂(x,y))。在复杂曲面生成的走刀矢量场中，走刀矢量分布呈现区域性。针对具有区域性走刀矢量分布的矢量场，开展基于走刀矢量分布的加工区域初分割。

根据矢量场的散度和旋度，可以确定该矢量场内的走刀矢量分布。在平面矢量场A(γ^*)＝(f₁(x,y),f₂(x,y))中任意一点的散度divA(γ^*)计算为：

式中，

为哈密顿算子。

任意一点的旋度rotA(γ^*)计算为：

式中，k为沿空间直角坐标系z轴正向的单位矢量。

计算曲面各点散度和旋度，根据散度和旋度是否为零，将平面矢量场分为以下四类：

(1)无散无旋场

此类平面矢量场散度与旋度均为零，各点走刀矢量分布呈现平行状态。

(2)无散有旋场

此类平面矢量场散度为零，旋度不为零，各点走刀矢量分布呈现涡旋状态。

(3)有散无旋场

此类平面矢量场散度不为零，旋度为零，各点走刀矢量分布呈现辐射状态。

(4)有散有旋场

此类平面矢量场散度与旋度均不为零，各点走刀矢量分布呈现无规则状态。

通过对平面矢量场各点散度与旋度求解，确定不同加工区域对应矢量场类型，将平面矢量场中散度或旋度突变点映射到曲面上，拟合曲线作为不同加工区域的边界，完成基于走刀矢量分布的加工区域初分割。

步骤3，确保进给运动平稳的加工区域细分

结合步骤2的加工区域初分割结果，对不同加工区域的空间矢量场进行流线拟合，根据机床旋转轴沿矢量场流线进给运动的运动学参数，开展确保进给运动平稳的加工区域细分。

以AC双转台型五轴数控机床为例，其旋转工作台A的转轴与机床坐标系的X轴平行，旋转工作台C的转轴与机床坐标系的Z轴平行。取加工区域空间矢量场的流线为r(ξ)，其中，ξ为流线方程的弧长参数。以曲面上流线处的法向量n(ξ)＝(n_x(ξ),n_y(ξ),n_z(ξ))作为机床旋转轴沿流线进给运动的刀轴矢量，计算为：

其中，S_u(u(ξ),v(ξ))、S_v(u(ξ),v(ξ))为曲面S(u,v)在流线处的一阶偏导数。

根据流线上各点刀轴矢量，求解机床A、C轴转角β_A(ξ)、β_C(ξ)分别为：

结合式(27)-(28)，计算机床A、C轴沿矢量场流线进给运动的角速度ω_A(ξ)和ω_C(ξ)分别为：

其中，t表示加工时间，v表示设定的进给速度。

计算机床A、C轴沿矢量场流线进给运动的角加速度α_A(ξ)和α_C(ξ)分别为：

由此可得到机床旋转轴沿矢量场流线进给运动的角速度ω(ξ)与角加速度α(ξ)的值为：

为避免机床旋转轴沿流线进给运动时出现明显切削振动，各加工区域内角速度和角加速度需满足：

其中，σ表示设定的平稳系数，σ越小则进给运动越平稳，其取值范围一般为1-5，s表示流线总弧长。

根据式(32)对进给运动是否平稳进行判断，以不满足该式的点作为旋转轴角速度或角加速度突变点，在加工区域初分割的基础上，求解得到各加工区域突变点并拟合曲线作为新增分区边界，完成确保进给运动平稳的加工区域细分。

步骤4，曲面不同加工区域加工轨迹拓扑形状设计

通过步骤2与步骤3将复杂曲面分割为走刀矢量分布规则且沿矢量场流线进给运动平稳的不同加工区域，结合走刀方式，对各加工区域进行适宜的加工轨迹拓扑形状设计。

常用走刀方式主要有行切走刀与环切走刀两种，针对不同加工区域对应的平面矢量场类型，分别在不同加工区域对应的平面上设计二维加工轨迹拓扑形状如下：

(1)无散无旋场

此类平面矢量场各点走刀矢量分布呈现平行状态，因此采用行切走刀加工。以各点走刀矢量为加工轨迹切向进行流线拟合，生成平行状的二维加工轨迹。

(2)无散有旋场

此类平面矢量场各点走刀矢量分布呈现涡旋状态，当区域边界曲线与矢量场流线平行时，则流线为完整涡旋线，采用环切走刀加工；当区域边界曲线与矢量场流线不平行时，则流线为局部涡旋线，采用行切走刀加工。以各点走刀矢量为加工轨迹切向进行流线拟合，生成涡旋状的二维加工轨迹。

(3)有散无旋场

此类平面矢量场各点走刀矢量分布呈现辐射状态，因此采用行切走刀加工。以各点走刀矢量为加工轨迹切向进行流线拟合，生成辐射状的二维加工轨迹。

(4)有散有旋场

此类平面矢量场各点走刀矢量分布呈现无规则状态，不宜采用行切或环切走刀加工，因此将有散有旋场转化为无散无旋场，鉴于无散无旋场走刀矢量分布呈现平行状态，故对有散有旋场走刀矢量调整如下：

取有散有旋场B(γ^*)＝(f₁(x,y),f₂(x,y))内部m个走刀矢量，分别表示为

其中i＝1,2,…,m。以走刀矢量平均值

作为调整后无散无旋场对应的走刀矢量，计算为：

以各点调整后的走刀矢量为加工轨迹切向进行流线拟合，生成平行状的二维加工轨迹。

在不同加工区域对应的平面上完成二维加工轨迹拓扑形状设计后，将二维加工轨迹映射到曲面上获得三维加工轨迹，从而实现基于矢量场的曲面分区加工轨迹拓扑形状设计。

本发明的显著效果和益处是提出了一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法，该方法针对复杂曲面分区加工过程中，加工区域分割与轨迹规划缺乏对走刀方式与进给运动平稳性综合考虑造成的局部误差超差及显著切削振动等问题，通过建立弦高误差与残高误差协同约束的走刀矢量场，完成基于走刀矢量分布的加工区域初分割及确保进给运动平稳的加工区域细分，并最终实现复杂曲面零件分区加工轨迹拓扑形状设计，有效减小了曲面分区加工的加工误差，提高了进给运动的平稳性，为复杂曲面高精高效加工提供技术支撑。

附图说明

图1—方法整体流程图。

图2—弦高误差与残高误差协同约束的走刀矢量场。

图3—经基于走刀矢量分布的加工区域初分割与确保进给运动平稳的加工区域细分后的曲面分区结果图。

图4—单一加工区域设计的加工轨迹拓扑形状。

图5—a)表示采用传统整体加工方法的加工表面粗糙度，b)表示采用本发明所述分区加工方法的加工表面粗糙度；Ra表示加工表面粗糙度(μm)。

图6—三坐标测量机测量的分别采用传统整体加工方法与本发明所述分区加工方法的加工表面轮廓对比图。

具体实施方式

结合技术方案与附图详细说明本发明的具体实施方式。

复杂曲面分区加工过程中，由于加工区域分割与轨迹规划缺乏对走刀方式与进给运动平稳性的综合考虑，易诱发局部误差超差及显著切削振动，严重影响曲面加工质量。针对这一情况，为了提高曲面分区加工的加工精度及进给运动的平稳性，发明了一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法，整体流程如附图1所示。

实施例采用AC双转台型五轴数控机床，以具有高陡度、频凹凸、曲率突变特征的复杂曲面为例，借助UG软件和MATLAB软件，详细说明本发明实施过程。

首先，利用UG软件对该复杂曲面进行建模，以局部陡度60°的对称回转面为主体，在该回转面上设计两个不规则凹坑，进而建立具有高陡度、频凹凸、曲率突变特征的复杂曲面，如图3所示。根据式(1)-(5)建立曲面各点弦高/残高误差与走刀方向法曲率/行距方向法曲率的函数关系，经式(6)-(12)完成两类误差的归一化处理，利用MATLAB软件，通过式(13)-(24)构建并求解两类误差的双目标优化模型，获得各点最优走刀方向，进而建立弦高误差与残高误差协同约束的走刀矢量场参见附图2。

其次，对该空间向量场投影获得平面向量场，通过MATLAB软件，由式(25)-(26)计算平面矢量场各点散度与旋度，其中散度的最大值与最小值分别为0.2299和0，旋度的最大值与最小值分别为0.1191和0.0045。根据各点散度和旋度是否为零对平面矢量场进行分类，以散度或旋度突变点拟合的曲线作为分区边界，完成该复杂曲面的加工区域初分割。其中，曲面凹坑部分对应有散有旋场，其余部分对应无散有旋场，

然后，利用MATLAB软件中streamslice函数对不同加工区域的空间矢量场进行流线拟合，通过式(27)-(31)计算机床旋转轴沿流线进给运动的角速度和角加速度，其中曲面凹坑处最大角速度为16.56°/s，最大角加速度为864.29°/s²。根据式(32)设定平稳系数为2，求解流线上旋转轴角速度或角加速度突变点并拟合曲线作为新增分区边界，完成该复杂曲面的加工区域细分，最终分区结果参见附图3。

最后，结合常见的行切与环切走刀方式，对不同加工区域加工轨迹拓扑形状进行设计，该复杂曲面凹坑处为有散有旋场，根据式(33)调整各点走刀矢量，生成相应的行切加工轨迹参见附图4，其余部分为无散有旋场，综合考虑区域边界与矢量场流线，分别生成相应的行切与环切加工轨迹，完成基于矢量场的曲面分区加工轨迹拓扑形状设计。

为验证此方法的有效性，进行传统整体加工与本发明所述分区加工对比实验。对比实验采用相同的加工工艺参数，分别为主轴转速5000r/min，切深0.25mm，进给速度300mm/min，行距0.2mm，步长1.8mm。实验结果表明，经本发明所述分区加工得到的曲面与经传统整体加工得到的曲面相比，分区加工的加工质量有显著提升。以工件高度40mm处的加工区域为代表，对其表面粗糙度及轮廓度进行测量。经传统整体加工得到的曲面表面粗糙度为3.5668μm，经本发明所述分区加工得到的曲面表面粗糙度为2.7193μm，较整体加工降低了23.76％，所示结果参见附图5，其中，图5a)表示采用传统整体加工方法的加工表面粗糙度，图5b)表示采用本发明所述分区加工方法的加工表面粗糙度。经传统整体加工得到的曲面轮廓度为153.6μm，经本发明所述分区加工得到的曲面轮廓度为91.84μm，较整体加工降低了40.21％，所示结果参见附图6。判定结果与实验结果一致，说明利用本发明的基于矢量场的曲面分区加工轨迹拓扑形状设计方法可有效减小曲面加工的加工误差，提高复杂曲面零件的加工质量，对工程实际中的加工轨迹规划具有重要的指导作用。

Claims (1)

1.一种基于矢量场的曲面分区加工轨迹拓扑形状设计方法，其特征是，该方法建立弦高误差与残高误差协同约束的走刀矢量场，经归一化构建两类误差的双目标优化模型，通过构造评价函数求解各点最优走刀方向，建立空间矢量场；计算投影后平面矢量场散度与旋度，根据其是否为零对矢量场分类，实现加工区域初分割；分析不同矢量场特征，拟合空间矢量场流线，计算机床旋转轴沿流线进给运动的运动学参数，根据其是否突变完成加工区域细分；结合走刀方式，实现曲面不同加工区域加工轨迹拓扑形状设计；方法的具体步骤如下：

步骤1，弦高误差与残高误差协同约束的走刀矢量场建立

取空间直角坐标系O-xyz中的曲面加工区域为S(u,v)＝(x(u,v),y(u,v),z(u,v))，P为曲面加工轨迹上任意一点，在行距d与步长f确定时，弦高误差e计算为：

其中，k_f为走刀方向法曲率；

取变量e^*计算为：

e^*＝k_f ² (2)

分别求解弦高误差e与变量e^*的单调性知，两者呈正相关关系，为降低后续计算难度，使用变量e^*表示弦高误差；

残高误差h计算为：

其中，k_d为行距方向法曲率，r为刀具有效半径；

由微分几何知，行距方向与走刀方向互相垂直，满足式(4)：

k_d+k_f＝2H (4)

其中，H为平均曲率，对于曲面任意一点，H为定值；

联立式(3)-(4)，残高误差表示为：

由式(2)与式(5)知，步长及行距确定时，弦高误差与残高误差仅由走刀方向法曲率决定；为求解各点最优走刀方向，以两类误差为目标函数构造双目标优化模型；

对两类误差归一化处理为：

其中，e_n和h_n分别表示归一化的弦高误差和残高误差，

和

表示弦高误差的最小值和最大值，h_min和h_max表示残高误差的最小值和最大值，其求解如下：

以k₁和k₂表示曲面主曲率(k₁＞k₂)，其计算为：

其中，K为高斯曲率，对于曲面任意一点，K为定值；

由式(2)与式(5)得，残高误差的最值为：

弦高误差的最大值为：

弦高误差的最小值为：

根据式(6)-(10)得，归一化的残高误差h_n为：

归一化的弦高误差e_n为：

弦高误差与残高误差均为归一化的弦高误差与残高误差，以两类误差为目标函数建立双目标优化模型为：

运用线性加权法构造评价函数Q为：

其中，α₁和α₂分别为弦高误差和残高误差在评价函数中所占权重，其计算如下：

其中，h_n1为弦高误差取最小值时对应的残高误差，e_n1为残高误差取最小值时对应的弦高误差，h_{n min}为残高误差最小值，e_{n min}为弦高误差最小值；

根据主曲率取值范围，计算双目标优化模型最优解如下：

(1)当k₂≥0时，取k_f＝k₂，e_n和h_n均取得最小值，因而，k_f＝k₂为双目标优化模型最优解，由欧拉定理知：

k_f＝k₁cos²θ+k₂sin²θ (16)

其中，θ为走刀方向与k₁对应主方向的夹角，此时，

(2)当-k₁≤k₂＜0时，评价函数Q为：

取

Q取得最小值，因而，

为双目标优化模型最优解，此时，

(3)当0≤k₁＜-k₂时，评价函数Q为：

取

Q取得最小值，因而，

为双目标优化模型最优解，此时，

(4)当k₁＜0时，评价函数Q为：

取

Q取得最小值，因而，

为双目标优化模型最优解，此时，

以(du:dv)表示曲面任意切方向，沿该方向的曲面第一基本形式I和第二基本形式II分别为：

其中，E、F、G为曲面第一类基本量，L、M、N为曲面第二类基本量，表达式如下：

其中，S_u、S_v为曲面S(u,v)的一阶偏导数，S_uu、S_uv、S_vv为曲面S(u,v)的二阶偏导数，n为曲面S(u,v)的法向量；

以(du¹:dv¹)表示主曲率k₁对应主方向，其计算为：

以(du^*:dv^*)表示最优走刀方向，其计算为：

其中，a、b、c分别表示为：

联立式(22)-(24)，求解获得曲面各点最优走刀方向，沿此方向建立单位走刀矢量，从而形成弦高误差与残高误差协同约束的走刀矢量场；

步骤2，基于走刀矢量分布的加工区域初分割

经步骤1获得的空间走刀矢量场在空间直角坐标系O-xyz中可表示为F(γ)＝(g₁(x,y,z),g₂(x,y,z),g₃(x,y,z))，其中γ为曲面各点沿最优走刀方向的单位走刀矢量；将该矢量场向xOy平面投影，沿各点投影后的走刀方向建立单位走刀矢量γ^*，获得平面矢量场A(γ^*)＝(f₁(x,y),f₂(x,y))；在复杂曲面生成的走刀矢量场中，走刀矢量分布呈现区域性；针对具有区域性走刀矢量分布的矢量场，开展基于走刀矢量分布的加工区域初分割；

根据矢量场的散度和旋度，可以确定该矢量场内的走刀矢量分布；在平面矢量场A(γ^*)＝(f₁(x,y),f₂(x,y))中任意一点的散度divA(γ^*)计算为：

式中，▽为哈密顿算子；

任意一点的旋度rotA(γ^*)计算为：

式中，k为沿空间直角坐标系z轴正向的单位矢量；

计算曲面各点散度和旋度，根据散度和旋度是否为零，将平面矢量场分为以下四类：

(1)无散无旋场

此类平面矢量场散度与旋度均为零，各点走刀矢量分布呈现平行状态；

(2)无散有旋场

此类平面矢量场散度为零，旋度不为零，各点走刀矢量分布呈现涡旋状态；

(3)有散无旋场

此类平面矢量场散度不为零，旋度为零，各点走刀矢量分布呈现辐射状态；

(4)有散有旋场

此类平面矢量场散度与旋度均不为零，各点走刀矢量分布呈现无规则状态；

通过对平面矢量场各点散度与旋度求解，确定不同加工区域对应矢量场类型，将平面矢量场中散度或旋度突变点映射到曲面上，拟合曲线作为不同加工区域的边界，完成基于走刀矢量分布的加工区域初分割；

步骤3，确保进给运动平稳的加工区域细分

结合步骤2的加工区域初分割结果，对不同加工区域的空间矢量场进行流线拟合，根据机床旋转轴沿矢量场流线进给运动的运动学参数，开展确保进给运动平稳的加工区域细分；

采用AC双转台型五轴数控机床，其旋转工作台A的转轴与机床坐标系的X轴平行，旋转工作台C的转轴与机床坐标系的Z轴平行；取加工区域空间矢量场的流线为r(ξ)，其中，ξ为流线方程的弧长参数；以曲面上流线处的法向量n(ξ)＝(n_x(ξ),n_y(ξ),n_z(ξ))作为机床旋转轴沿流线进给运动的刀轴矢量，计算为：

其中，S_u(u(ξ),v(ξ))、S_v(u(ξ),v(ξ))为曲面S(u,v)在流线处的一阶偏导数；

根据流线上各点刀轴矢量，求解机床A、C轴转角β_A(ξ)、β_C(ξ)分别为：

结合式(27)-(28)，计算机床A、C轴沿矢量场流线进给运动的角速度ω_A(ξ)和ω_C(ξ)分别为：

其中，t表示加工时间，v表示设定的进给速度；

计算机床A、C轴沿矢量场流线进给运动的角加速度α_A(ξ)和α_C(ξ)分别为：

由此可得到机床旋转轴沿矢量场流线进给运动的角速度ω(ξ)与角加速度α(ξ)的值为：

为避免机床旋转轴沿流线进给运动时出现明显切削振动，各加工区域内角速度和角加速度需满足：

其中，σ表示设定的平稳系数，σ越小则进给运动越平稳，其取值范围为1-5，s表示流线总弧长；

根据式(32)对进给运动是否平稳进行判断，以不满足该式的点作为旋转轴角速度或角加速度突变点，在加工区域初分割的基础上，求解得到各加工区域突变点并拟合曲线作为新增分区边界，完成确保进给运动平稳的加工区域细分；

步骤4，曲面不同加工区域加工轨迹拓扑形状设计

通过步骤2与步骤3将复杂曲面分割为走刀矢量分布规则且沿矢量场流线进给运动平稳的不同加工区域，结合走刀方式，对各加工区域进行适宜的加工轨迹拓扑形状设计；

常用走刀方式主要有行切走刀与环切走刀两种，针对不同加工区域对应的平面矢量场类型，分别在不同加工区域对应的平面上设计二维加工轨迹拓扑形状如下：

(1)无散无旋场

此类平面矢量场各点走刀矢量分布呈现平行状态，因此采用行切走刀加工；以各点走刀矢量为加工轨迹切向进行流线拟合，生成平行状的二维加工轨迹；

(2)无散有旋场

此类平面矢量场各点走刀矢量分布呈现涡旋状态，当区域边界曲线与矢量场流线平行时，则流线为完整涡旋线，采用环切走刀加工；当区域边界曲线与矢量场流线不平行时，则流线为局部涡旋线，采用行切走刀加工；以各点走刀矢量为加工轨迹切向进行流线拟合，生成涡旋状的二维加工轨迹；

(3)有散无旋场

此类平面矢量场各点走刀矢量分布呈现辐射状态，因此采用行切走刀加工；以各点走刀矢量为加工轨迹切向进行流线拟合，生成辐射状的二维加工轨迹；

(4)有散有旋场

此类平面矢量场各点走刀矢量分布呈现无规则状态，不宜采用行切或环切走刀加工，因此将有散有旋场转化为无散无旋场，鉴于无散无旋场走刀矢量分布呈现平行状态，故对有散有旋场走刀矢量调整如下：

取有散有旋场B(γ^*)＝(f₁(x,y),f₂(x,y))内部m个走刀矢量，分别表示为

其中i＝1,2,···,m；以走刀矢量平均值

作为调整后无散无旋场对应的走刀矢量，计算为：

以各点调整后的走刀矢量为加工轨迹切向进行流线拟合，生成平行状的二维加工轨迹；

在不同加工区域对应的平面上完成二维加工轨迹拓扑形状设计后，将二维加工轨迹映射到曲面上获得三维加工轨迹，从而实现基于矢量场的曲面分区加工轨迹拓扑形状设计。

Images