细分曲面切削加工仿真技术

王彩年

通用技术集团大连机床有限责任公司 辽宁大连 116620

1 序言

在工业领域中，为实现机械特殊的功能，具有复杂形状零件的应用越来越广泛，这对零件设计和加工都提出了更高的要求。对比传统NURBS建模方法，细分曲面建模方法在复杂形状零件的设计和加工方面更具优势。

数控加工仿真时，需要首先将加工模型进行离散化处理，从最初的离散化模型[1]，逐渐发展为基于实体模型表示的数控加工仿真技术[2]。HUNT等[3]最早将CSG与数控加工仿真相结合，SU等[4]将CSG用于实时碰撞干涉检测。B-rep表示法最早由BRAID[5]提出。在B-rep表示法的基础上，为提高仿真效率，FLEISIG等[6]创建了加速几何建模方法。但是由于细分曲面加工模型往往具有复杂的形状和大量的数据，直接采用现有方法会导致数控加工仿真效率降低[7-10]。综上所述，本文通过研究毛坯模型构建、切削刀位点判断和误差分析等数控加工技术的关键环节，提出高性能数控加工仿真技术方法。

2 毛坯模型网格细分

毛坯模型如图1所示。一般数控加工都是由上至下分层切削加工，因此毛坯上表面常常被设置为加工部分。毛坯模型上表面的网格离散非常重要。实际上，粗加工和精加工对毛坯模型表示精度有着不同的需求。粗加工的加工模型的表示精度要远低于精加工。因此，可以先构建一个较粗的网格作为粗加工的毛坯模型。具体方法如下。

图1 毛坯模型

细分曲面初始控制网格，通过细分n次获得粗加工模型的极限网格，以ds表示同一细分曲面面片中网格边的长度。进一步进行平均值计算。计算所求的平均值dm即为粗加工毛坯上表面离散网格间距，具体计算公式为

式中，m表示精加工模型的细分次数，m>n。

粗加工仿真后，再将毛坯模型上表面细分（m-n）次，即可获得精加工毛坯模型。

3 切削刀位点判断

3.1 遍历毛坯节点

如图2所示，毛坯上表面区域范围由大矩形表示，刀具扫描面区域范围由剖面线矩形表示。毛坯节点的几何坐标位于扫描面内部，表示方法为

图2 毛坯节点与扫描面的关系

式中，Xmax（Xmin）为毛坯顶点在X坐标轴向的最大（最小）坐标分量（mm）；Ymax（Ymin）为毛坯顶点在Y坐标轴向的最大（最小）坐标分量（mm）；xmax（xmin）为刀具扫描面在X坐标轴向的最大（最小）坐标分量（mm）；ymax（ymin）为刀具扫描面在Y坐标轴向的最大（最小）坐标分量（mm）。

同理可得扫描面与毛坯不存在交集满足

3.2 判断刀具扫描域

遍历毛坯节点后，分析毛坯节点与刀具扫描面（见图3）的位置关系。如果毛坯节点位于刀具扫描面内部，判断毛坯节点能否进行下一步的计算分析，主要依据是节点能否先后符合式（2）～式（5）的要求。

图3 毛坯节点与刀具扫描面

设两刀位点分别是M（XM、YM、ZM）和N（XN、YN、ZN），则以这两个刀位点为基础构成的X轴、Y轴以及Z轴的方向向量，分别可以用如下形式表示：a=XN－XM，b=YN－YM，c=ZN－ZM。由于以上3个方向向量具有对应的比例关系，依照此关系进一步分析研究可得

式中：x、y、z分别表示毛坯节点笛卡尔坐标系下的对应数值；Xc、Yc、Zc分别表示毛坯节点向刀具轨迹做垂线垂心的三维坐标值。

当t＝0时，可判定刀位点M（XM，YM，ZM）与毛坯节点的几何坐标发生重叠；当t＝1时，代表刀位点N（XN，YN，ZN）与毛坯节点的几何坐标发生重叠。此时表明辅助坐标轴已经完成构建。

当t＝0时，表示辅助坐标轴的原点可以用M点表示，根据式（7）能够实现对原始坐标轴中的坐标进行转化，转化后的数值结果即为辅助坐标轴的三维坐标值。

3.3 切削毛坯节点位置计算

以环形铣刀为例获得刀具的扫描域形状，将扫描域分为若干部分来进行毛坯节点位置的数学表达。其他刀具同理。

如图4所示，环形铣刀的刀具扫描域可以分为三部分：以M点和N点为分界，M点以左的区域表示为第一部分；M点和N点之间的区域表示为第二部分；N点以右的区域表示为第三部分。其中大径半圆的尺寸用R表示；小径半圆的尺寸用r表示，两半圆间的距离用Rc表示。以此可确定判断环形铣刀的刀具扫描域材料切削的数学表达式为

用P（x、y、z）对毛坯中的任意节点进行表示。当其符合式（8）时，可对二者位置关系进行判定，确定刀具扫描体此时位于刀具扫描域的第一部分。依据环形铣刀的自身特点，毛坯节点的位置可进一步判定为两种情况：一是毛坯节点的几何坐标处在圆环之间；二是毛坯节点的几何坐标处在环形铣刀平面圆内的区域。具体计算公式为

倘若点P符合式（9），可判定毛坯节点的几何坐标位于底面圆的内部区域，以此可确定变更P节点处在扫描面内的Z坐标，同时X坐标和Y坐标保持不变，计算公式为

倘若P点符合式（11），即可判定毛坯节点处在圆环区域中。

当毛坯节点位于刀具扫描体内M点和N点之间的第二部分时，毛坯节点P（x、y、z）满足式（12）。中间部位可进一步分为两种情况：一是毛坯节点位于刀具扫描域的圆环区域内；二是毛坯节点位于刀具扫描域的矩形部位，即毛坯节点在刀具底部平面区域内。

倘若P点符合式（13），可判定P点几何坐标处于加工刀具底面扫描域的区域内，则保持扫描面的X坐标和Y坐标不变的同时，进一步确定该毛坯节点P的Z坐标为

倘若P点不能同时符合式（11）及式（13），可判定P点几何坐标处于刀具扫描域的圆环区域内，对刀具扫描体中四分之一球的圆心计算是此时的研究难点。依据式（15）可以实现对球心Ot（xt、yt、zt）的计算

进而根据公式（16）得到z的值

倘若P（x、y、z）点符合式（17），可判定刀具扫描体位于刀具扫描域中N点以右的第三部分。依据环形铣刀的自身特点，毛坯节点的位置可进一步判定为两种情况：一是毛坯节点的位置处在圆环间的区域；二是毛坯节点位置处在环形铣刀平面圆内的区域。

倘若P点符合式（18），可判定毛坯节点位置处在底部的平面圆区域内，保持扫描面毛坯节点P的X、Y坐标原数值，Z坐标修改为式（19）

倘若P点符合式（20），可判定毛坯的节点几何坐标处在圆环的几何区域范围内，保持扫描面的X坐标和Y坐标原数值的同时，可依据环形铣刀圆角的公式模型进行计算，以此确定P点的Z坐标。

分析比较毛坯节点的z值以及式（19）计算得出的Z值。倘若Z

4 毛坯加工误差分析

在加工模型上搜索毛坯模型顶点的最近距离点对是毛坯加工误差分析的基础。基于细分曲面多分辨力采样方法，可以实现最近距离点对的快速搜索。图5表示多分辨率采样搜索方法，其中正方形表示择优点；黑色圆点表示采样点。最后在黑色区域择优选取最近点，其所在的空间范围为下次采样范围。

图5 多分辨率采样搜索方法

因此择优点可表示为

式中，表示第n次选取的择优点；表示择优点的行索引；代表择优点的列索引。

对于第n次搜索的4个采样点、、、，第n-1次择优点索引值为、，则可以得到从第2次到第n－1次获得采样点的公式

5 加工实例验证

采用图6所示的细分曲面模型作为加工模型，粗加工采用环形铣刀；半精加工和精加工采用球头铣刀对加工模型进行数控加工试验。粗加工、半精加工和精加工阶段的刀具轨迹如图7所示。

图6 加工模型

图7 3个加工阶段的刀具轨迹

各加工阶段的加工仿真和加工试验结果如图8所示。

图8 3个加工阶段的仿真和试验结果

从加工误差分析色图（见图9）可以看出，颜色由蓝到黄表示加工误差由小变大。图中误差分析结果表明，该加工试验的加工误差较小。

图9 加工误差分析色图

6 结束语

根据环形铣刀数学模型，建立刀具的扫描域，搜索毛坯模型中位于扫描域内部的数据，最终改变被切削点的数值，利用该方式能够完成数控加工切削的高效仿真。

利用细分曲面多分辨率采样方法，实现了最近距离点对的快速搜索。该方法为加工仿真的快速实现和误差分析的高效计算奠定了基础。