五轴侧铣加工刀具路径光顺方法
【技术领域】
[0001] 本发明涉及计算机图形学和五轴数控侧铣加工技术领域，具体地，涉及一种五轴 侧铣加工刀具路径光顺方法。
【背景技术】
[0002] 复杂曲面零件的加工是在机床上进行的。因此，在规划刀具加工路径的时候，需要 考虑机床的动态性能，避免机床的旋转轴角度发生很大的变化。以往的侧铣加工刀具路径 规划是在工件坐标系下进行的，通过使刀轴方向在工件坐标系下光顺从而希望机床各轴的 运动也光顺。但由于五轴机床的后置处理是非线性模型，在工件坐标系下光顺的五轴刀具 加工轨迹不代表机床的各轴运动也光顺，尤其是当刀具经过机床的奇异区域时，即使刀具 的轴线发生很小的变化，机床的旋转轴会产生很大的摆动角度。因此在规划刀具加工路径 时需要考虑机床的特性。对于侧铣加工来说，刀轴方向的微小变化可能会造成比较大的几 何误差，所以在规划五轴侧铣加工刀具路径时还需要考虑几何误差。
[0003]目前没有发现同本发明类似技术的说明或报道，也尚未收集到国内外类似的资 料。
【发明内容】
[0004] 针对现有技术中的缺陷，本发明的目的在于利用双参数球族包络理论，提供一种 五轴侧铣加工刀具路径光顺方法。
[0005] 本发明是通过以下技术方案实现的。
[0006] -种五轴侧铣加工刀具路径光顺方法，包括如下步骤：
[0007] 步骤S1 :计算初始刀位，根据选用的机床，通过后置处理算法计算初始刀位刀轴 方向对应的机床旋转轴角度；
[0008] 步骤S2 :插值刀具刀尖点位置和机床旋转轴角度，得到刀具加工路径的刀尖点位 置轨迹线和刀轴矢量曲线；
[0009] 步骤S3:计算刀具路径对应的机床旋转轴旋转角度曲线的刚度矩阵；
[0010] 步骤S4 :根据微分演化算法和点到曲面的距离函数计算设计面上任一点到刀具 包络面的最近距离；
[0011] 步骤S5 :根据加权最小二乘法建立五轴侧铣加工刀具路径光顺模型；
[0012] 步骤S6 :根据高斯-牛顿法迭代求解光顺模型；
[0013] 步骤S7 :得到优化后的刀具路径的刀尖点位置轨迹线和光顺的机床旋转轴旋转 角度曲线，进而输出NC(数控）文件。
[0014] 优选地，所述步骤S2包括如下步骤：
[0015] 步骤S2. 1 :刀具的刀轴方向采用机床的两个旋转角度表示。
[0016] 优选地，针对五轴双摆头机床，刀轴方向0(t)表示为：
[0017]
(2)
[0018] 其中CA(t)，Cc(t)分别表示刀路对应的机床旋转角度曲线，Ojt)表示刀轴矢量 〇(t)的X坐标，0j(t)表示0(t)的y坐标，0k(t)表示0(t)的Z坐标；侧铣加工刀具路径 采用工件坐标系下的刀尖点轨迹P(t)与刀轴矢量曲线0(t)表示；因此刀具路径S(w;a，t) 表示为：
[0019] S(w；a,t) =P(t)+a·Η· 0 (t) (2)
[0020] 其中κ= eHw表示刀具路径的控制点集合，包括刀具刀尖点轨迹线控制 点和机床旋转角度曲线的控制点，a、t分别表示刀具路径的参数，Η表示刀具刃长；如果刀 尖点位置曲线的控制点个数为1，则m= 31+1+1 ;利用公式（2)表示侧铣加工刀具路径，直 接优化机床的两个旋转角度。
[0021] 优选地，所述步骤S3包括如下步骤：
[0022] 步骤S3. 1 :令机床旋转角度曲线的控制点分别为其中1为刀尖点 位置曲线的棹制点个教.W弄示1維这:教空間，刹用公忒⑶庶量刀县路径的光顺件F_ - th:
[0023]
[0024] 其矩阵形式为：
[0025]
[0026] 其中K表示刚度矩阵，刚度矩阵K中的元素计算公式为：
[0027]
[0028] 优选地，所述步骤S4包括如下步骤：
[0029] 步骤S4. 1:根据双参数球族包络法，无需构造包络面即可计算设计面上任一点p 到刀具包络面X(a，t)的距离；设计面上任一点到刀具包络面X(a，t)的距离dp,x(w)计算 公式为：
[0030]
[0031] 其中r(w;a，t)表示刀具包络球族的半径；
[0032] 步骤S4.2 :利用微分演化算法，根据公式（4)计算设计面上任一点p到刀具包络 面的最近距离，并得到最近点在刀具路径S(a，t)上的参数（a，t)。
[0033] 优选地，所述步骤S5包括如下步骤：
[0034] 步骤S5. 1 :对刀具路径进行光顺优化，需要调整刀具路径的控制点，这样会增大 刀具包络面与设计面之间的几何偏差，因此，光顺刀具路径时需要同时考虑几何偏差，通过 加权最小二乘法可以建立光顺机床旋转角度的侧铣加工刀具路径规划模型：
[0035]
[0036] 其中λ是光顺权重。
[0037] 优选地，所述步骤S6包括如下步骤：
[0038] 步骤6. 1:公式（5)中的模型是非线性最小二乘问题，采用高斯-牛顿法迭代求 解，求解过程如下：
[0039]气
?公式（5)表示为：
[0040]
[0041] 其中Y表示由刚度矩阵K组成的矩阵；!Τ表示刀具路径的控制点集合。
[0042] 设wk为当前解，将目标函数在w作一阶线性泰勒展开，得到相应的线性最小二 乘问题：
[0043]
[0044] 其中Δwk表示最优解，则最优解Δwk为：
[0045] Awk= [(Ak)T(Ak) +AY] 1 ((Ak)V+(wk)TY)
[0046] 求得Δwk后，令wk+1 =wk+Δwk，并进行下一次的迭代，当迭代次数超过最大迭代次 数或者Awk的变化小于设定值，则算法结束。
[0047] 优选地，根据输出刀位个数或者刀尖点位置弓高误差要求，在优化后的刀尖点位 置轨迹线和机床旋转角度曲线上采样，然后输出NC文件。
[0048] 与现有技术相比，本发明具有如下的有益效果：
[0049] 1、本发明解决了五轴机床侧铣加工曲面时出现的机床旋转轴剧烈变化的问题，通 过优化算法使机床的旋转轴运动平滑，同时考虑了加工后零件的几何精度要求。
[0050] 2、本发明适用于自由曲面、直纹面或类直纹面曲面五轴侧铣加工。
【附图说明】
[0051] 通过阅读参照以下附图对非限制性实施例所作的详细描述，本发明的其它特征、 目的和优点将会变得更明显：
[0052] 图1为本发明的流程图；
[0053] 图2为本发明中刀具加工路径刀尖点位置轨迹线示意图；
[0054] 图3为本发明中刀具加工路径刀轴方向示意图；
[0055] 图4为本发明中五轴双摆头机床示意图；
[0056] 图5为本发明中S形设计曲面示意图；
[0057]图6为本发明中优化前后刀具路径对应的机床C轴变化曲线；
[0058]图7为本发明中优化前后刀具路径对应的机床A轴变化曲线。
【具体实施方式】
[0059] 下面对本发明的实施例作详细说明：本实施例在以本发明技术方案为前提下进行 实施，给出了详细的实施方式和具体的操作过程。应当指出的是，对本领域的普通技术人员 来说，在不脱离本发明构思的前提下，还可以做出若干变形和改进，这些都属于本发明的保 护范围。
[0060] 实施例
[0061] 本实施例提供了一种五轴侧铣加工刀具路径光顺方法，下面结合附图对本实施例 进行详细描述。
[0062] 本实施例中，如图2、图3、图4以及图5所示，本实施例根据双参数球族包络理论 计算几何偏差，利用加权最小二乘法建立五轴侧铣加工刀具路径光顺模型。
[0063] 五轴侧铣刀具路径可以采用工件坐标系下的刀尖点轨迹P(t)与刀轴矢量曲线 〇(t)表示
[0064]S(w；a,t) =P(t)+a·Η·0(t)
[0065] 其中Wt_f1?,.....f€ur表示刀具路径的控制点集合，包括刀具刀尖点轨迹线控制 点和机床旋转角度曲线的控制点，Η表示刀具刃长。设计面上一点到刀具包络面X(a，t)的 距离的计算公式为
[0066]
[0067]