CN104238455B - 一种利用球头刀铣削自由曲面的加工方法 - Google Patents

一种利用球头刀铣削自由曲面的加工方法 Download PDF

Info

Publication number
CN104238455B
CN104238455B CN201410526635.5A CN201410526635A CN104238455B CN 104238455 B CN104238455 B CN 104238455B CN 201410526635 A CN201410526635 A CN 201410526635A CN 104238455 B CN104238455 B CN 104238455B
Authority
CN
China
Prior art keywords
curved surface
point
points
degenerate
machining
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201410526635.5A
Other languages
English (en)
Other versions
CN104238455A (zh
Inventor
李迎光
李振瀚
刘旭
李强
马斯博
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing University of Aeronautics and Astronautics
Original Assignee
Nanjing University of Aeronautics and Astronautics
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Nanjing University of Aeronautics and Astronautics filed Critical Nanjing University of Aeronautics and Astronautics
Priority to CN201410526635.5A priority Critical patent/CN104238455B/zh
Publication of CN104238455A publication Critical patent/CN104238455A/zh
Application granted granted Critical
Publication of CN104238455B publication Critical patent/CN104238455B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Numerical Control (AREA)

Abstract

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

Description

一种利用球头刀铣削自由曲面的加工方法
技术领域
本发明涉及一种CAD/CAM技术,尤其是一种自由曲面的加工方法,具体地说是一种利用球头刀铣削自由曲面的加工方法。
背景技术
随着产品对气动性能要求的不断提高,复杂自由曲面大量出现航空航天、汽车及船舶等产品的零件中,由于几何性质的复杂性,复杂自由曲面大多采用球头刀铣削方式加工,该加工方式对曲面形状的适用性强,且适用于三轴及多轴机床加工,但球头刀铣削方式的加工效率较低,对于大型复杂曲面零件如飞机翼身融合曲面,其球头刀铣削精加工的时间通常在100小时以上。
切削宽度优化是优化曲面加工效率的重要方式之一,较大的切削宽度可以减少曲面加工轨迹的长度,从而缩短加工时间。目前面向切宽优化的自由曲面刀轨生成方法主要是美国北卡罗来纳州立大学Yuan-shin Lee教授提出的基于Machining Potential Field的方法,该方法首先在曲面上构造一系列曲线,选出平均切宽最大的曲线作为初始轨迹线,偏置初始轨迹线生成曲面的加工轨迹,当偏置得到的轨迹线的切削效率系数小于设定值时则停止轨迹计算,在未生成刀轨区域重新构造初始轨迹线进行轨迹计算。该方法能够一定程度上减少相邻刀轨间的重叠区域,实现切宽的优化,但该方法的本质是一种被动式的分区加工方法,是一种贪婪式算法,因此仅能得到局部优化的加工结果。
本发明公开了一种球头刀铣削自由曲面加工方法。该方法以张量的方式度量曲面上任意切触点处的切削宽度,构建切削宽度二阶张量场,提取张量场中的三分退化点,并以此为起点搜索曲面内部边界点,实现对曲面加工区域的划分,每个区域称为一个曲面加工特征。在每一个曲面加工特征中构建数条曲线,保证曲线上每点的切向与该点处切宽最大的进给方向重合,选出平均切宽最大的曲线为该区域内的初始加工轨迹线,采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。该方法能够保证所生成的刀轨尽可能遵循对应最大切宽的进给方向在曲面范围内的区域化分布规律,从而减少刀轨长度,实现加工效率的提高。
发明内容
本发明针对当前复杂曲面球头刀铣削加工时无法进行切宽全局优化的问题,提出一种利用球头刀铣削自由曲面的加工方法。
一种利用球头刀铣削自由曲面的加工方法,其特征是它主要包括以下实施步骤:
步骤1:对曲面进行网格划分,得到曲面包含的网格列表及每个网格包含的网格点;
步骤2:计算每个网格点处的切宽张量系数矩阵,构建切宽张量场;
步骤3:搜索切宽张量场中包含的退化点,提取三分退化点;
步骤4:以三分退化点为起点搜索曲面内部边界点,构造曲面内部边界,得到曲面加工特征。
步骤5:在每一个曲面加工特征中构建数条曲线,保证曲线上每点的切向与该点处切宽最大的进给方向重合,选出平均切宽最大的曲线为该区域内的初始加工轨迹线;
步骤6:采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。
上述步骤1中所述的曲面网格划分具体方法如下:将曲面投影至零件底面B上,如图2所示,提取包围该投影区域的最小矩形T,人工设置T在两个方向上的离散步距m和n,对T进行离散得到B上的一系列网格BG,将BG沿B的法向投影至曲面S上,得到曲面网格列表SG。
上述步骤2所述的切宽张量场构造方法如下:设u,v分别为曲面的两个参数,计算G中所有网格点处的法向n、并带入切宽张量系数矩阵H
其中r为球头刀半径,G为曲面第一标准型系数矩阵,且D为曲面第二基本型系数矩阵,且G-1为G的反向矩阵,在切宽张量场中的每一个曲面离散点处,切宽表示为Hαβduαduβ(α,β=0,or 1)其中u1和u2相当于u和v。H的最大特征方向即为切宽最大的进给方向,当图3所示的曲面选用直径为5mm的球头刀加工时,该曲面上各离散点的切宽最大进给方向可通过求解各点处H的最大特征方向获得,且在曲面范围内,切宽最大的进给方向呈现出区域化分布的趋势;
上述步骤3所述的搜索切宽张量场中包含的退化点具体实现方法如下:切宽张量场的退化点处满足各个进给方向上的切宽大小相同,因此在退化点处H具有两个相同大小的特征值,此时H满足如图4所示,对于网格SGi,首先判断其是否包含退化点,设其四个顶点分别记为Pi1、Pi2、Pi3和Pi4,各点处(H11-H22)分别记为ti1、ti2、ti3和ti4,H12分别记为si1、si2、si3和si4,若满足则SGi内包含退化点,其中tmin和tmax分别为ti1、ti2、ti3和ti4中的最小值与最大值,smin和smax分别为si1、si2、si3和si4中的最小值与最大值。当确定SGi内包含退化点之后,采用双线性插值方法求解去除落在[0,1]外的u、v解便得到SGi包含的退化点DPi,图5中所示的{P1,P2,P3,P4,P5,P6}为图3切宽张量场中包含的退化点;
上述步骤3所述的提取三分退化点的具体实现方法如下:计算各退化点处的δ=αλ-βγ,若δ<0,则该退化点为三分退化点,如图5中的点P1、P2、P5和P6
上述步骤4中所述的曲面内部边界构造方法如图6所示,其具体实现方法如下:
以三分退化点为起始点,搜索曲面内部边界点,构造曲面内部边界线,尽量确保曲面内部边界线上每一点的切向与该点处切宽最大的进给方向一致。对于三分退化点T_DPi,首先在其切平面内以该点为圆心做半径为R的圆Ci(R通常小于1mm,可由使用者自行设定),对Ci按角度步距Δω(Δω为一非常小的角度值)进行离散,得到离散点集合CPi,将CPi沿T_DPi在曲面上的法向Ni投影至曲面上得到PCPi,通过计算H的最大特征方向获取PCPi中所有点的切宽最大进给方向向量集合MVi,并将MVi沿Ni投影回切平面中得到PMVi,过T_DPi与CPi中所有点做向量得到可选方向集合PVi,计算所有PVi与其对应的PMVi间的夹角记作ANGi。对于任意三分退化点T_DPi,在CPi中可找到三个点使其对应的ANG近似为0,满足该条件的点称为种子点SP,将T_DPi与SPij连接成线段并沿Ni方向投影至曲面上便可以得到一条曲面内部边界曲线段CSij,SPij在曲面上的投影点PSPij为曲面上的内部边界点,因此对于一个三分退化点,其周围存在三个曲面内部边界点,即每一个三分退化点仅能被三条内部边界穿过。当从T_DPi进行曲面内部边界构造时,首先检查该三分退化点是否已被内部边界经过,若已由三条内部边界经过该点,则说明无需对该点进行边界构造,若小于三条,则按上述方法求出T_DPi临近的三个内部边界点,去掉落在穿过T_DPi内部边界上的边界点,利用余下的点进行搜索构造下一段曲线段。
对于曲面内部边界点PSPij,采取同样的方法构造半径为R的圆Cij并对其按Δω进行离散得到点集CPij,将CPij沿CPij处的曲面法向Nij投影至曲面上得到点集PCPij,通过计算H的最大特征方向获取PCPij中所有点的切宽最大进给方向向量集合MVij,并将MVij沿Nij投影回切平面中得到PMVij,过PSPij与CPij中所有点做向量得到可选方向集合PVij,计算所有PVij与其对应的PMVij间的夹角记作ANGij。此时CPij中可搜索得到两个种子点SPij1和SPij2使其对应的ANG近似为0,分别过SPij1和SPij2与PSPij做向量TVij1和TVij2,将上一段曲面内部边界曲线段沿Nij投影至PSPij的切平面中得到向量TSij,分别计算TVij1和TVij2与TSij的夹角θij1与θij2,提取两个夹角中较大的那个夹角对应的种子点记作SPij,将IPk与PSPij连接成线段并沿Nij方向投影至曲面上得到内部边界曲线段并添加到当前的内部边界曲线上,用SPij在曲面上的投影点替换PSPij继续进行边界的构造。
曲面内部边界点搜索停止的条件为:①当前曲面内部边界点PSPij所做的圆Cij有部分落在曲面外侧,此时过PSPij做上一曲面内部边界曲线段的切线向量,将该向量沿Nij投影至曲面上,提取该投影曲线段落在曲面内部的部分,将其作为这条曲面内部边界线的最后一段;②当所得到的曲面内部边界曲线段穿过某一三分退化点,或两者间最小距离小于某一给定的微小长度值,此时将该三分退化点看做是下一个曲面内部边界点,提取当前曲面内部边界点与该三分退化点间的最短曲线段作为这条曲面内部边界线的最后一段。
当所有三分退化点处都有三条曲面内部边界线经过则完成整个曲面区域的划分,每一个曲面区域称为一个曲面加工特征SMF。
上述步骤5中所述的各个曲面加工特征中的初始加工轨迹线构造具体实现方法如下:
如图7中流程所示,首先采用人工编辑的方法在各个曲面加工特征对应的曲面区域内创建一系列点,将每一个人工创建的点Pi以步骤4中的方法找到其对应的两个种子点SPi1和SPi2,过Pi分别与SPi1和SPi2做直线段,并沿Pi处的曲面法向Ni投影至曲面上得到曲线段CSi1和CSi2,将CSi1和CSi2保存为到曲线段集合CCi中,相同的方式得到两个种子点在曲面上对应的投影点PSPi1和PSPi2,选择PSPi1为起点按照步骤4中PSPij的搜索方法构造曲线段并添加到CCi中,构造终止条件为:①若当前搜索得到的曲线段穿过PSPi2或者两者间的最小距离小于某一给定的微小长度值时,则以提取PSPij与PSPi2间的曲线段作为CCi中的最后一段,从而形成一条封闭的曲线段;②若在当前PSPij的切平面内构造的半径为R的圆有部分落在该曲面加工特征外时,则过该点做上一曲线段的切向量,将该切向量沿该点处的曲面法向投影至曲面上,取出该投影线段落在该曲面加工特征内部的部分添加到CCi中,以SPi2为起点继续进行搜索,若以SPi2的起点的搜索已完成,则过Pi点的曲线构造完成。
对于曲面加工特征SMFs,按照上述方法在其对应的曲面区域内创建了数条曲线,对于其中的任意一条,通过以下方法计算其平均切宽CWa
其中α,β=0,or 1,u1和u2分别为曲面的两个参数,τ为曲线的弧长参数。
提取出CWa最大的那条曲线作为SMFs的初始刀具轨迹线。
本发明的有益效果是:
(1)本发明实现了曲面的分区加工,使刀轨能够遵循最大切宽进给方向在曲面上的区域化分布规律,尽可能使每一个刀触点处的进给方向与切宽最大的进给方向重合,从而缩短刀轨长度,实现加工效率的提高;
(2)本发明使用张量来度量球头刀加工曲面的切宽,从理论上突破了以切宽最大化为目标的曲面加工区域划分方法,能够自动生成准确的曲面内部边界,实现曲面加工特征的自动提取;
(3)本发明在选择每个曲面加工特征的初始轨迹线时,采用连续积分的方法,能够得到各条待选曲线的准确平均切宽,保证了初始轨迹线的优化性;
(4)本发明适用于曲率变化剧烈的复杂曲面加工,加工效率提升显著。
附图说明
图1为本发明的球头刀铣削自由曲面加工方法流程图。
图2为本发明的曲面网格划分方法示意图,其中m、n分别是Y和X方向上的离散步距。
图3为本发明的曲面S上各点利用切宽张量系数矩阵H求出的切削宽度最大进给方向分布图,其中V为切宽最大进给方向。
图4为本发明的曲面上的网格SGi,其中Pi1、Pi2、Pi3和Pi4为它的四个顶点,u和v为曲面的两个参数。
图5为图3中曲面S上的退化点分布,其中P1、P2、P5和P6为三分退化点。
图6为本发明的以三分退化点为起点的曲面内部边界构造算法流程图。
图7为本发明的各曲面加工特征对应曲面区域内的待选初始轨迹线构造算法流程图。
图8为本发明的球头刀铣削自由曲面加工方法的应用案例。图8中:图8a中T_DP1和T_DP2为两个三分退化点,IB11、IB12和IB13为从T_DP1出发搜索得到的三条内部边界,IB21、IB22和IB23为从T_DP2出发搜索得到的三条内部边界。该曲面被这6条内部边界线分成了5个加工特征,分别为SMF1、SMF2、SMF3、SMF4和SMF5;图8b中所示IPC1、IPC2、IPC3、IPC4和IPC5分别为这5个特征选出来的初始轨迹线,采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨,如图8b所示。
具体实施方式
下面是结合附图对本发明的技术方案进行详细说明。
如图1-8所示。
一种利用球头刀铣削自由曲面的加工方法,如图1所示,它主要包括以下实施步骤:
步骤1:对曲面进行网格划分,得到曲面包含的网格列表及每个网格包含的网格点;
步骤2:计算每个网格点处的切宽张量系数矩阵,构建切宽张量场;
步骤3:搜索切宽张量场中包含的退化点,提取三分退化点;
步骤4:以三分退化点为起点搜索曲面内部边界点,构造曲面内部边界,得到曲面加工特征。
步骤5:在每一个曲面加工特征中构建数条曲线,保证曲线上每点的切向与该点处切宽最大的进给方向重合,选出平均切宽最大的曲线为该区域内的初始加工轨迹线;
步骤6:采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。
上述步骤1中所述的曲面网格划分具体方法如下:将曲面投影至零件底面B上,如图2所示,提取包围该投影区域的最小矩形T,设置T在两个方向上的离散步距m和n,对T进行离散得到B上的一系列网格BG,将BG沿B的法向投影至曲面S上,得到曲面网格列表SG。
上述步骤2所述的切宽张量场构造方法如下:设u,v分别为曲面的两个参数,计算G中所有网格点处的法向n、并带入切宽张量系数矩阵H
其中r为球头刀半径,G为曲面第一标准型系数矩阵,且D为曲面第二基本型系数矩阵,且G-1为G的反向矩阵,在切宽张量场中的每一个曲面离散点处,切宽表示为Hαβduαduβ(α,β=0,or 1)其中u1和u2相当于u和v。H的最大特征方向即为切宽最大的进给方向,当图3所示的曲面选用直径为5mm的球头刀加工时,该曲面上各离散点的切宽最大进给方向可通过求解各点处H的最大特征方向获得,且在曲面范围内,切宽最大的进给方向呈现出区域化分布的趋势;
上述步骤3所述的搜索切宽张量场中包含的退化点具体实现方法如下:切宽张量场的退化点处满足各个进给方向上的切宽大小相同,因此在退化点处H具有两个相同大小的特征值,此时H满足如图4所示,对于网格SGi,首先判断其是否包含退化点,设其四个顶点分别记为Pi1、Pi2、Pi3和Pi4,各点处(H11-H22)分别记为ti1、ti2、ti3和ti4,H12分别记为si1、si2、si3和si4,若满足则SGi内包含退化点,其中tmin和tmax分别为ti1、ti2、ti3和ti4中的最小值与最大值,smin和smax分别为si1、si2、si3和si4中的最小值与最大值。当确定SGi内包含退化点之后,采用双线性插值方法求解去除落在[0,1]外的u、v解便得到SGi包含的退化点DPi,图5中所示的{P1,P2,P3,P4,P5,P6}为图3切宽张量场中包含的退化点;
上述步骤3所述的提取三分退化点的具体实现方法如下:计算各退化点处的δ=αλ-βγ,若δ<0,则该退化点为三分退化点,如图5中的点P1、P2、P5和P6
上述步骤4中所述的曲面内部边界构造方法如图6所示,其具体实现方法如下:
以三分退化点为起始点,搜索曲面内部边界点,构造曲面内部边界线,尽量确保曲面内部边界线上每一点的切向与该点处切宽最大的进给方向一致。对于三分退化点T_DPi,首先在其切平面内以该点为圆心做半径为R的圆Ci(R通常小于1mm,可由使用者自行设定),对Ci按角度步距Δω(Δω为一非常小的角度值)进行离散,得到离散点集合CPi,将CPi沿T_DPi在曲面上的法向Ni投影至曲面上得到PCPi,通过计算H的最大特征方向获取PCPi中所有点的切宽最大进给方向向量集合MVi,并将MVi沿Ni投影回切平面中得到PMVi,过T_DPi与CPi中所有点做向量得到可选方向集合PVi,计算所有PVi与其对应的PMVi间的夹角记作ANGi。对于任意三分退化点T_DPi,在CPi中可找到三个点使其对应的ANG近似为0,满足该条件的点称为种子点SP,将T_DPi与SPij连接成线段并沿Ni方向投影至曲面上便可以得到一条曲面内部边界曲线段CSij,SPij在曲面上的投影点PSPij为曲面上的内部边界点,因此对于一个三分退化点,其周围存在三个曲面内部边界点,即每一个三分退化点仅能被三条内部边界穿过。当从T_DPi进行曲面内部边界构造时,首先检查该三分退化点是否已被内部边界经过,若已由三条内部边界经过该点,则说明无需对该点进行边界构造,若小于三条,则按上述方法求出T_DPi临近的三个内部边界点,去掉落在穿过T_DPi内部边界上的边界点,利用余下的点进行搜索构造下一段曲线段。
对于曲面内部边界点PSPij,采取同样的方法构造半径为R的圆Cij并对其按Δω进行离散得到点集CPij,将CPij沿CPij处的曲面法向Nij投影至曲面上得到点集PCPij,通过计算H的最大特征方向获取PCPij中所有点的切宽最大进给方向向量集合MVij,并将MVij沿Nij投影回切平面中得到PMVij,过PSPij与CPij中所有点做向量得到可选方向集合PVij,计算所有PVij与其对应的PMVij间的夹角记作ANGij。此时CPij中可搜索得到两个种子点SPij1和SPij2使其对应的ANG近似为0,分别过SPij1和SPij2与PSPij做向量TVij1和TVij2,将上一段曲面内部边界曲线段沿Nij投影至PSPij的切平面中得到向量TSij,分别计算TVij1和TVij2与TSij的夹角θij1与θij2,提取两个夹角中较大的那个夹角对应的种子点记作SPij,将IPk与PSPij连接成线段并沿Nij方向投影至曲面上得到内部边界曲线段并添加到当前的内部边界曲线上,用SPij在曲面上的投影点替换PSPij继续进行边界的构造。
曲面内部边界点搜索停止的条件为:①当前曲面内部边界点PSPij所做的圆Cij有部分落在曲面外侧,此时过PSPij做上一曲面内部边界曲线段的切线向量,将该向量沿Nij投影至曲面上,提取该投影曲线段落在曲面内部的部分,将其作为这条曲面内部边界线的最后一段;②当所得到的曲面内部边界曲线段穿过某一三分退化点,或两者间最小距离小于某一给定的微小长度值,此时将该三分退化点看做是下一个曲面内部边界点,提取当前曲面内部边界点与该三分退化点间的最短曲线段作为这条曲面内部边界线的最后一段。
当所有三分退化点处都有三条曲面内部边界线经过则完成整个曲面区域的划分,每一个曲面区域称为一个曲面加工特征SMF。
上述步骤5中所述的各个曲面加工特征中的初始加工轨迹线构造具体实现方法如下:
如图7中流程所示,首先采用人工编辑的方法在各个曲面加工特征对应的曲面区域内创建一系列点,将每一个人工创建的点Pi以步骤4中的方法找到其对应的两个种子点SPi1和SPi2,过Pi分别与SPi1和SPi2做直线段,并沿Pi处的曲面法向Ni投影至曲面上得到曲线段CSi1和CSi2,将CSi1和CSi2保存为到曲线段集合CCi中,相同的方式得到两个种子点在曲面上对应的投影点PSPi1和PSPi2,选择PSPi1为起点按照步骤4中PSPij的搜索方法构造曲线段并添加到CCi中,构造终止条件为:①若当前搜索得到的曲线段穿过PSPi2或者两者间的最小距离小于某一给定的微小长度值时,则以提取PSPij与PSPi2间的曲线段作为CCi中的最后一段,从而形成一条封闭的曲线段;②若在当前PSPij的切平面内构造的半径为R的圆有部分落在该曲面加工特征外时,则过该点做上一曲线段的切向量,将该切向量沿该点处的曲面法向投影至曲面上,取出该投影线段落在该曲面加工特征内部的部分添加到CCi中,以SPi2为起点继续进行搜索,若以SPi2的起点的搜索已完成,则过Pi点的曲线构造完成。
对于曲面加工特征SMFs,按照上述方法在其对应的曲面区域内创建了数条曲线,对于其中的任意一条,通过以下方法计算其平均切宽CWa
其中α,β=0,or 1,u1和u2分别为曲面的两个参数,τ为曲线的弧长参数。
提取出CWa最大的那条曲线作为SMFs的初始刀具轨迹线。
根据以上所述方法对图8中的曲面进行加工轨迹生成,该曲面精加工所用球头刀的直径为6mm,表面最大残高为0.05mm。图8a中T_DP1和T_DP2为该曲面上的两个三分退化点,IB11、IB12和IB13为从T_DP1出发搜索得到的三条内部边界,IB21、IB22和IB23为从T_DP2出发搜索得到的三条内部边界。该曲面被这6条内部边界线分成了5个加工特征,分别为SMF1、SMF2、SMF3、SMF4和SMF5;图8b中所示IPC1、IPC2、IPC3、IPC4和IPC5分别为这5个特征选出来的初始轨迹线,采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨,如图b所示,经测算,利用本方法生成的轨迹线总长度为12921.6mm,而采用传统的方法,将整个曲面作为一个加工区域,采用Iso-parametric轨迹生成算法,以曲面最长的边界起始轨迹线,得到的轨迹总长度为14870.8mm,利用本方法生成的轨迹在总长度上缩短了13.12%。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。

Claims (7)

1.一种利用球头刀铣削自由曲面的加工方法,其特征是它主要包括以下实施步骤:
步骤1:对曲面进行网格划分,得到曲面包含的网格列表及每个网格包含的网格点;
步骤2:计算每个网格点处的切宽张量系数矩阵H,构建切宽张量场;
步骤3:搜索切宽张量场中包含的退化点,提取三分退化点;
步骤4:以三分退化点为起点搜索曲面内部边界点,构造曲面内部边界,得到曲面加工特征;
步骤5:在每一个曲面加工特征中构建数条曲线,保证曲线上每点的切向与该点处切宽最大的进给方向重合,选出平均切宽最大的曲线为该曲面加工特征的初始加工轨迹线;
步骤6:采用偏置初始加工轨迹线的方式计算各曲面加工特征的刀轨从而得到整个曲面的加工轨迹。
2.根据权利要求1所述的利用球头刀铣削自由曲面的加工方法,其特征是所述的对曲面进行网格划分是指:将曲面投影至零件底面B上,提取包围该投影区域的最小矩形T,人工设置T在两个方向上的离散步距m和n,对T进行离散得到B上的一系列网格BG,将BG沿B的法向投影至曲面S上,得到曲面网格列表SG。
3.根据权利要求1所述的利用球头刀铣削自由曲面的加工方法,其特征是所述的构建切宽张量场是指:设u,v分别为曲面的两个参数,计算曲面第一标准型系数矩阵G中所有网格点处的法向n、并带入切宽张量系数矩阵H,
H = 1 2 ( 2 r G - D ) G - 1
其中r为球头刀半径,S为曲面的参数表达式,D为曲面第二基本型系数矩阵,且G-1为G的反向矩阵,在切宽张量场中的每一个曲面离散点处,根据爱因斯坦求和约定,切宽可表示为Hαβduαduβ,其中α,β=0或者1,u1=u,u2=v;H的最大特征方向即为切宽最大的进给方向,且在曲面范围内,切宽最大的进给方向呈现出区域化分布的趋势。
4.根据权利要求1所述的利用球头刀铣削自由曲面的加工方法,其特征是所述的搜索切宽张量场中包含的退化点按以下方法:切宽张量场的退化点处满足各个进给方向上的切宽大小相同,因此在退化点处切宽张量系数矩阵H具有两个相同大小的特征值,此时H满足其中H11、H12、H21和H22分别为切宽张量系数矩阵H第一行第一列、第一行第二列、第二行第一列及第二行第二列;对于网格SGi,其中i表示第i个网格,首先判断其是否包含退化点,设其四个顶点分别记为Pi1、Pi2、Pi3和Pi4,各点处(H11-H22)分别记为ti1、ti2、ti3和ti4,H12分别记为si1、si2、si3和si4,若满足则SGi内包含退化点,其中tmin和tmax分别为ti1、ti2、ti3和ti4中的最小值与最大值,smin和smax分别为si1、si2、si3和si4中的最小值与最大值;当确定SGi内包含退化点之后,采用双线性插值方法求解其中u,v分别为曲面的两个参数,去除落在[0,1]外的u、v解便得到SGi包含的退化点DPi
5.根据权利要求1所述的利用球头刀铣削自由曲面的加工方法,其特征是所述的提取三分退化点按以下方法:计算各退化点处的δ=αλ-βγ,其中H11、H12、H21和H22分别为切宽张量系数矩阵H第一行第一、第二项以及第二行第一、第二项,u,v分别为曲面的两个参数,若δ<0,则该退化点为三分退化点。
6.根据权利要求1所述的利用球头刀铣削自由曲面的加工方法,其特征是所述的构造曲面内部边界是指:以三分退化点为起始点,搜索曲面内部边界点,构造曲面内部边界线,确保曲面内部边界线上每一点的切向与该点处切宽最大的进给方向一致;对于三分退化点T_DPi,其中i表示第i个三分退化点,i最小为1,最大为三分退化点的总数;首先在曲面的切平面内以该点为圆心做半径为R的圆Ci,R通常小于1mm,可由使用者自行设定,对Ci按角度步距Δω进行离散,Δω小于1°,得到离散点集合CPi,将CPi沿T_DPi在曲面上的法向Ni投影至曲面上得到PCPi,通过计算H的最大特征方向获取PCPi中所有点的切宽最大进给方向向量集合MVi,并将MVi沿Ni投影回切平面中得到PMVi,过T_DPi与CPi中所有点做向量得到可选方向集合PVi,计算所有PVi与其对应的PMVi间的夹角记作ANGi;对于任意三分退化点T_DPi,在CPi中可找到三个点使其对应的ANG为0,满足该条件的点称为种子点SP,将T_DPi与SP连接成线段并沿Ni方向投影至曲面上便可以得到一条曲面内部边界曲线段CS,SP在曲面上的投影点PSP为曲面上的内部边界点,因此对于一个三分退化点,其周围存在三个曲面内部边界点,即每一个三分退化点仅能被三条内部边界穿过;当从T_DPi进行曲面内部边界构造时,首先检查该三分退化点是否已被内部边界经过,若已由三条内部边界经过该点,则说明无需对该点进行边界构造,若小于三条,则按上述方法求出T_DPi临近的三个内部边界点,去掉落在穿过T_DPi内部边界上的边界点,利用余下的点进行搜索构造下一段曲线段;对于任一PSPij,其中PSPij表示T_DPi对应的第j个内部边界点,j为1到3的整数,采取上述同样的方法构造半径为R的圆Cij并对其按Δω进行离散得到点集CPij,将CPij沿CPij处的曲面法向Nij投影至曲面上得到点集PCPij,通过计算H的最大特征方向获取PCPij中所有点的切宽最大进给方向向量集合MVij,并将MVij沿Nij投影回切平面中得到PMVij,过PSPij与CPij中所有点做向量得到可选方向集合PVij,计算所有PVij与其对应的PMVij间的夹角记作ANGij;此时CPij中可搜索得到两个种子点SPij1和SPij2使其对应的ANG为0,分别过SPij1和SPij1与PSPij做向量TVij1和TVij2,将上一段曲面内部边界曲线段沿Nij投影至PSPij的切平面中得到向量TSij,分别计算TVij1和TVij2与TSij的夹角θij1与θij2,提取两个夹角中较大的那个夹角对应的种子点记作SPij,将SPij与PSPij连接成线段并沿Nij方向投影至曲面上得到内部边界曲线段并添加到当前的内部边界曲线上,用SPij在曲面上的投影点作为下一个内部边界点PSP继续进行边界的构造;曲面内部边界点搜索停止的条件为:①当前曲面内部边界点PSP所做的圆C有部分落在曲面外侧,此时过PSP做上一曲面内部边界曲线段的切线向量,将该向量沿法向N投影至曲面上,提取该投影曲线段落在曲面内部的部分,将其作为这条曲面内部边界线的最后一段;②当所得到的曲面内部边界曲线段穿过某一三分退化点,或两者间最小距离小于某一给定的微小长度值,此时将该三分退化点看做是下一个曲面内部边界点,提取当前曲面内部边界点与该三分退化点间的最短曲线段作为这条曲面内部边界线的最后一段;当所有三分退化点处都有三条曲面内部边界线经过则完成整个曲面区域的划分,每一个曲面区域称为一个曲面加工特征SMF。
7.根据权利要求1或6所述的利用球头刀铣削自由曲面的加工方法,其特征是各个曲面加工特征中的初始加工轨迹线构造方法为:首先采用人工编辑的方法在各个曲面加工特征对应的曲面区域内创建一系列点,找到任意点Pi对应的两个种子点SPi1和SPi2,使得SPi1和SPi2在曲面上的法向投影点与Pi连线形成的向量与种子点处的切宽最大进给方向在Pi切平面内的夹角为0,过Pi分别与SPi1和SPi2做直线段,并沿Pi处的曲面法向Ni投影至曲面上得到曲线段CSi1和CSi2,将CSi1和CSi2保存为到曲线段集合CCi中,相同的方式得到两个种子点在曲面上对应的投影点PSPi1和PSPi2,选择PSPi1为起点构造曲线段并添加到CCi中;构造终止条件为:①若当前搜索得到的曲线段穿过PSPi2或者两者间的最小距离小于某一给定的微小长度值时,则以提取PSPij与PSPi2间的曲线段作为CCi中的最后一段,从而形成一条封闭的曲线段;②若在当前PSPij的切平面内构造的半径为R的圆有部分落在该曲面加工特征外时,则过该点做上一曲线段的切向量,将该切向量沿该点处的曲面法向投影至曲面上,取出该投影线段落在该曲面加工特征内部的部分添加到CCi中,以SPi2为起点继续进行搜索,若以SPi2的起点的搜索已完成,则过Pi点的曲线构造完成;对于曲面加工特征SMFs,按照上述方法在其对应的曲面区域内创建了数条曲线,对于其中的任意一条,通过以下方法计算其平均切宽CWa
其中α,β=0或者1,Hαβ为切宽张量系数矩阵H对应的各项,u1和u2分别为曲面的两个参数,τ为曲线的弧长参数,提取出CWa最大的那条曲线作为SMFs的初始刀具轨迹线。
CN201410526635.5A 2014-10-09 2014-10-09 一种利用球头刀铣削自由曲面的加工方法 Active CN104238455B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410526635.5A CN104238455B (zh) 2014-10-09 2014-10-09 一种利用球头刀铣削自由曲面的加工方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410526635.5A CN104238455B (zh) 2014-10-09 2014-10-09 一种利用球头刀铣削自由曲面的加工方法

Publications (2)

Publication Number Publication Date
CN104238455A CN104238455A (zh) 2014-12-24
CN104238455B true CN104238455B (zh) 2017-03-15

Family

ID=52226767

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410526635.5A Active CN104238455B (zh) 2014-10-09 2014-10-09 一种利用球头刀铣削自由曲面的加工方法

Country Status (1)

Country Link
CN (1) CN104238455B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104656558B (zh) * 2015-01-15 2017-07-07 华中科技大学 一种获取刀具轨迹中加工特征点的方法及其应用
CN105785913B (zh) * 2016-04-06 2018-03-20 武汉工程大学 基于机床轴速度限制的刀具路径切削方向优化方法
CN106950916B (zh) * 2017-04-26 2019-07-16 山东理工大学 基于ab型五轴数控机床环形刀加工刀轴矢量光顺方法
CN108549325B (zh) * 2018-05-23 2020-10-16 合肥工业大学 一种自由曲面弧长参数曲线加工轨迹生成方法
CN110394487A (zh) * 2019-07-30 2019-11-01 哈尔滨理工大学 一种淬硬钢自由曲面试件及铣削工艺
CN112034785A (zh) * 2020-08-31 2020-12-04 江西江铃集团新能源汽车有限公司 一种加工路径生成方法、装置、存储介质及设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000123199A (ja) * 1998-10-14 2000-04-28 Ricoh Co Ltd 3次元形状の切断装置と切断方法及び切断処理プログラムを格納した記憶媒体
EP1241621A2 (en) * 2001-03-16 2002-09-18 Mitsubishi Denki Kabushiki Kaisha System and method for generating adaptively sampled distance fields with bounded distance trees
CN101908235A (zh) * 2010-06-23 2010-12-08 苏州科技学院 B样条曲面建模新方法
CN102637216A (zh) * 2011-12-14 2012-08-15 南京航空航天大学 一种复杂组合曲面的数控侧铣加工刀轨生成方法
CN102981454A (zh) * 2012-11-01 2013-03-20 武汉华中数控股份有限公司 一种基于势能场及能量泛函优化的刀具轨迹路径产生方法
CN103971411A (zh) * 2013-01-24 2014-08-06 岳天祥 利用三维物体的空间曲面采样点对空间曲面建模的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000123199A (ja) * 1998-10-14 2000-04-28 Ricoh Co Ltd 3次元形状の切断装置と切断方法及び切断処理プログラムを格納した記憶媒体
EP1241621A2 (en) * 2001-03-16 2002-09-18 Mitsubishi Denki Kabushiki Kaisha System and method for generating adaptively sampled distance fields with bounded distance trees
CN101908235A (zh) * 2010-06-23 2010-12-08 苏州科技学院 B样条曲面建模新方法
CN102637216A (zh) * 2011-12-14 2012-08-15 南京航空航天大学 一种复杂组合曲面的数控侧铣加工刀轨生成方法
CN102981454A (zh) * 2012-11-01 2013-03-20 武汉华中数控股份有限公司 一种基于势能场及能量泛函优化的刀具轨迹路径产生方法
CN103971411A (zh) * 2013-01-24 2014-08-06 岳天祥 利用三维物体的空间曲面采样点对空间曲面建模的方法

Also Published As

Publication number Publication date
CN104238455A (zh) 2014-12-24

Similar Documents

Publication Publication Date Title
CN104238455B (zh) 一种利用球头刀铣削自由曲面的加工方法
CN104238456B (zh) 一种非球头刀铣削加工自由曲面的方法
CN106969778B (zh) 一种多无人机协同施药的路径规划方法
CN108681787A (zh) 基于改进双向快速扩展随机树算法的无人机路径优化方法
CN107610223A (zh) 基于LiDAR点云的电力塔三维重建方法
CN105488562B (zh) 一种基于多因子粒子群算法的不规则件排样方法
CN107578472B (zh) 三维表面三角网络模型的角度优化方法及装置
CN107146285A (zh) 一种基于曲面拟合的任意自由曲面网格划分方法
CN102981454A (zh) 一种基于势能场及能量泛函优化的刀具轨迹路径产生方法
CN110773738B (zh) 基于多边形几何特征识别的激光扫描路径分区域规划方法
CN105654552B (zh) 一种面向任意分布大规模点云数据的快速Delaunay构网方法
CN106903889A (zh) 一种stl模型的快速切片算法
CN110442994A (zh) 一种基于图论的3d打印切片处理方法
CN108549327A (zh) 一种复杂结构件五轴增减材混合加工序列规划方法
CN106553760A (zh) 埋入式进气道唇口的设计方法
CN102915560A (zh) 一种阈值无关的机载激光雷达点云滤波方法及装置
CN102938165B (zh) 基于型面特征逼近的产品stl模型光顺方法
CN105488830A (zh) 一种二维轮廓区域单轨迹填充生成方法
CN102184538B (zh) 一种基于动态轮廓的合成孔径雷达sar图像自动分割方法
CN103037168B (zh) 基于复合型PCNN的平稳Surfacelet域多聚焦图像融合方法
CN103942366B (zh) 基于四段有理Bézier曲线表示的曲率连续的翼型及其生成方法
CN104615679A (zh) 一种基于人工免疫网络的多智能体数据挖掘方法
CN110001066B (zh) 一种三维打印最少分区填充方向确定方法
CN102930589A (zh) 一种非结构直角网格相贯线修正方法
CN105654498A (zh) 基于动态局部搜索和免疫克隆自动聚类的图像分割方法

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Li Yingguang

Inventor after: Li Zhenhan

Inventor after: Liu Xu

Inventor after: Li Qiang

Inventor after: Ma Sibo

Inventor before: Li Yingguang

Inventor before: Liu Xu

Inventor before: Li Qiang

Inventor before: Ma Sibo

COR Change of bibliographic data
C14 Grant of patent or utility model
GR01 Patent grant