复杂曲面侧铣加工的刀具优化和加工轨迹生成方法技术

本发明专利技术公开了一种复杂曲面侧铣加工的刀具优化和加工轨迹生成方法，所示方法包括如下步骤：步骤一、交互设计初始曲线A(t)和B(t),生成初始的刀具加工路径；步骤二、生成初始的刀具形状；步骤三、优化加工路径和刀具形状；步骤四、误差计算；步骤五、刀具曲面的计算。本发明专利技术的主要贡献是提出一种计算框架，能够直接计算复杂曲面侧铣加工的刀具形状的连续参数化曲面表示，具有较高的精度和较高的计算效率。跟已有方法不同，本发明专利技术的方法直接优化一个B样条表示的刀具曲面，得到连续表示的刀具形状，且该方法能够处理过切控制和刀具的抖动控制。

【技术实现步骤摘要】
复杂曲面侧铣加工的刀具优化和加工轨迹生成方法
本专利技术涉及一种面向侧铣加工的刀具形状优化方法，具体涉及一种复杂曲面高精度侧铣加工的刀具优化和与所生成刀具对应的加工轨迹生成方法。
技术介绍
在航空发动机的叶轮叶片等复杂曲面的加工中，常用的加工技术有端铣加工和侧铣加工。端铣加工技术存在容易产生干涉以及加工效率低等问题。侧铣加工属于线接触加工，具有较宽的加工行距，因此加工效率远大于端铣加工。侧铣加工具有较高的加工精度，一般应用于精加工阶段。侧铣加工要求在每一个时刻刀具和工件曲面相切于一条曲线，这样严格的几何约束使得侧铣加工的路径规划比较困难，因此目前侧铣加工主要用于直纹面的加工。对非可展直纹面或者自由曲面，采用标准刀具难以得到高精度的加工结果，需要针对设计曲面计算出合适的刀具曲面形状以及相应的加工轨迹，然后定制刀具用于实际加工。已有方法可以得到离散的刀具曲面的表示，表示精度不够高，且已有侧铣加工刀具形状优化方法没有考虑运动抖动和过切的控制。
技术实现思路
本专利技术的目的是提供一种复杂曲面高精度侧铣加工的刀具优化和与所生成刀具对应的加工轨迹生成方法。跟已有方法不同，本专利技术的方法直接优化一个B样条表示的刀具曲面，得到连续参数曲面表示的刀具形状，同时得到刀具运动路径的连续表示，且该方法能够控制过切和刀具抖动。本专利技术的目的是通过以下技术方案实现的：一种复杂曲面侧铣加工的刀具优化和加工轨迹生成方法，包括如下步骤：步骤一、交互设计初始曲线A(t)和B(t)，生成初始加工路径：(1)用户交互指定K+1个线段，每一个线段的两个端点位于曲面S上，假设得到线段序列L0＝(A0，B0)，...，LK＝(AK，BK)，其中Ai,Bi,i＝0,...,K表示线段的两个端点，每一个线段的长度等于L，L为刀具的长度；(2)将(1)中得到的线段的端点Ai或者Bi沿着曲面法向的方向分别移动距离RA和RB，即更新为：Ai＝Ai+N(Ai)·RA；Bi＝Bi+N(Bi)·RB；其中，N(P)表示P点的曲面法向；(3)基于端点Ai、Bi，定义刀具轴线两端的轨迹曲线为：其中，Ni,d(t)是d次的B样条基函数；连接A(t)和B(t)上相同参数点的线段表示刀具的轴线，轴线上的点定义的轴迹面R(t,s)如下：R(t,s)＝A(t)(1-s)+B(t)st∈[0,1],s∈[0,L]；步骤二、生成初始的刀具形状：(1)在轴迹面R(t,s)上取采样点Ri,j：Ri,j＝R(ti,sj)＝A(ti)(1-sj)+B(ti)sj；其中，sj∈[0,1],j＝0,...,W，是参数域s内的采样值，W+1是采样点的个数；ti∈[0,1],j＝0,...,Y，是参数域t内的采样值，Y+1是采样点的个数；(2)计算每一个采样点Ri,j到设计曲面S的最近距离di,j，在平面直角坐标系X-Y内定义数据点(sj,di,j)；(3)将数据点(sj,di,j)拟合成如下形式的函数：其中，fi,i＝0,...,Nf是1维控制点，控制点个数是Nf+1；Ni,d(s)是d次B样条基函数；假设刀具曲面的轴线是L(s)，r(s)定义了一个平面曲线C，曲线C绕X轴旋转360度得到的旋转曲面作为初始刀具曲面；步骤三、优化加工路径和刀具形状：(1)在轴迹面R(t,s)上采样得到数据点R(ti,sj)：R(ti,sj),i＝0,...,N1,j＝0,...,N2；其中，N1和N2是整数，分别代表了在参数t和参数s上的采样点数目；计算采样点R(ti,sj)在设计曲面S上的垂足点S(ui,j,vi,j)，计算设计曲面上点S(ui,j,vi,j)处的单位法向量N(ui,j,vi,j)；(2)求解下式定义的优化问题：其中，(A′(ti)-B′(ti))·(B(ti)-A(ti))＝0,i＝0,...,N是刚体运动约束，意味着刀具轴线在运动中保持长度固定；FApprox是考虑过切控制的距离误差函数，FEnergy是刀头运动的稳定性控制函数；步骤四、误差计算：计算R(ti,sj)在设计曲面S上的垂足点S(ui,j,vi,j)，计算加工误差其中G是采样点的数目；如果误差满足要求，则输出确定加工路径的函数A(t)、B(t)和确定刀具形状的函数r(s)，转到步骤五，否则转到步骤三；步骤五、刀具曲面的计算；以坐标系的X轴作为z轴，以任意垂直于z的向量作为x轴，以z×x(×表示向量的叉乘)作为y轴建立直角坐标系，在该坐标系下刀具曲面的参数化如下式定义：其中，s和θ是刀具曲面的参数。本专利技术中，优化问题中的变量包括曲线A(t)的控制点Ai,i＝0,…,K,B(t)的控制点Bi,i＝0,…,K(K代表控制点的数目)，以及函数r(s)的控制点fi,i＝0,...,Nf(Nf代表控制点的数目)。相比于现有技术，本专利技术具有如下优点：(1)本专利技术的主要贡献是提出一种计算框架，能够直接计算复杂曲面侧铣加工的刀具形状的连续参数化曲面表示，具有较高的精度和较高的计算效率。已有的方法只能计算刀具曲面的离散表示，若想得到较高的精度需要较多的变量，而较多的变量也会增加计算量，降低计算效率。(2)本专利技术的计算框架直接优化刀具的加工路径，通过引入刀具抖动相关的函数，实现了对刀具抖动的控制。(3)本专利技术的计算框架可以对加工中的过切进行控制，能生成无过切的加工结果。附图说明图1为设计曲面上的线段序列；图2为将线段序列沿曲面法向抬高；图3为生成初始轴迹面；图4为一般刀具的形状定义；图5为初始刀具形状和运动；图6为本专利技术方法的流程图；图7为本专利技术方法生成的加工轨迹(左)和对应的近似圆锥曲面的刀具形状(右)；图8为本专利技术方法生成的加工轨迹(左)和对应的凸曲面的刀具形状(右)；图9为本专利技术方法生成的加工轨迹(左)和对应的凹曲面的刀具形状(右)。具体实施方式下面结合附图对本专利技术的技术方案作进一步的说明，但并不局限于此，凡是对本专利技术技术方案进行修改或者等同替换，而不脱离本专利技术技术方案的精神和范围，均应涵盖在本专利技术的保护范围中。在五轴数控加工技术中，侧铣加工采用刀具的侧面对材料进行去除。高速旋转的刀头可以看作一个旋转曲面，因此高精度的侧铣加工对应刀具旋转曲面运动生成的包络曲面高精度逼近加工曲面。本专利技术适用的刀具可以是铣削刀具或者是研磨刀具。刀具曲面的形状和加工轨迹是紧密联系的两个问题，需要同时求解。已有的刀头形状方法不能直接得到连续的刀具曲面。本专利技术提出了一种针对给定设计曲面S的刀具曲面和加工轨迹的计算方法，如图6所示，具体包括如下步骤：1、求解刀具曲面和加工轨迹假设设计曲面为S，刀具曲面φ沿着加工运动轨迹生本文档来自技高网...

【技术保护点】
1.一种复杂曲面侧铣加工的刀具优化和加工轨迹生成方法，其特征在于所述方法包括如下步骤：/n步骤一、交互设计曲线A(t)和B(t)，生成初始加工路径：/n(1)用户交互指定K+1个线段，每一个线段的两个端点位于曲面S上，假设得到线段序列L

【技术特征摘要】
1.一种复杂曲面侧铣加工的刀具优化和加工轨迹生成方法，其特征在于所述方法包括如下步骤：
步骤一、交互设计曲线A(t)和B(t)，生成初始加工路径：
(1)用户交互指定K+1个线段，每一个线段的两个端点位于曲面S上，假设得到线段序列L0＝(A0，B0)，...，LK＝(AK，BK)，其中Ai,Bi,i＝0,...,K表示线段的两个端点，每一个线段的长度等于L，L为刀具的长度；
(2)将(1)中得到的线段的端点Ai或者Bi沿着曲面法向的方向分别移动距离RA和RB，即更新为：
Ai＝Ai+N(Ai)·RA；
Bi＝Bi+N(Bi)·RB；
其中，N(P)表示P点的曲面法向；
(3)基于端点Ai、Bi，定义刀具轴线两端的轨迹曲线为：






其中，Ni,d(t)是d次的B样条基函数；
连接A(t)和B(t)上相同参数点的线段表示刀具的轴线，轴线上的点定义的轴迹面R(t,s)如下：
R(t,s)＝A(t)(1-s)+B(t)st∈[0,1],s∈[0,L]；
步骤二、生成初始的刀具形状：
(1)在轴迹面R(t,s)上取采样点Ri,j：
Ri,j＝R(ti,sj)＝A(ti)(1-sj)+B(ti)sj；
其中，sj∈[0,1],j＝0,...,W，是参数域s内的采样值，W+1是采样点的个数；ti∈[0,1],j＝0,...,Y，是参数域t内的采样值，Y+1是采样点的个数；
(2)计算每一个采样点Ri,j到设计曲面S的最近距离di,j，在平面直角坐标系X-Y内定义数据点(sj,di,j)；
(3)将数据点(sj,di,j)拟合成如下形式的函数：



其中，fi,i＝0,...,Nf是1维控制点，控制点个数是Nf+1；Ni,d(s)是d次B样条基函数；
假设刀具曲面的轴线是L(s)，r(s)定义了一个平面曲线C，曲线C绕X轴旋转360度得到的旋转曲面作为初始刀具曲面；
步骤三、优化加工路径和刀具形状：
(1)在轴迹面R(t,s)上采样得到数据点R(ti,sj)：
R(ti,sj),i＝0,...,N1,j＝0,...,N2；
其中，N1和N2是整数，分别代表了在参数t和参数s上的采样点数目；
计算采样点R(ti,sj)在设...

【专利技术属性】
技术研发人员：伯彭波，张晓苗，
申请(专利权)人：哈尔滨工业大学威海，
类型：发明
国别省市：山东;37