CN102357666B - 自由曲面叶片平底刀三坐标端铣加工方法 - Google Patents
自由曲面叶片平底刀三坐标端铣加工方法 Download PDFInfo
- Publication number
- CN102357666B CN102357666B CN 201110201226 CN201110201226A CN102357666B CN 102357666 B CN102357666 B CN 102357666B CN 201110201226 CN201110201226 CN 201110201226 CN 201110201226 A CN201110201226 A CN 201110201226A CN 102357666 B CN102357666 B CN 102357666B
- Authority
- CN
- China
- Prior art keywords
- partiald
- blade
- free
- sub
- form surface
- 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.)
- Expired - Fee Related
Links
- 238000003801 milling Methods 0.000 title claims abstract description 21
- 238000012545 processing Methods 0.000 claims abstract description 42
- 238000000034 method Methods 0.000 claims abstract description 29
- 239000013598 vector Substances 0.000 claims abstract description 11
- 239000011159 matrix material Substances 0.000 claims abstract description 4
- 238000003754 machining Methods 0.000 claims description 36
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000012805 post-processing Methods 0.000 claims description 3
- 241000237509 Patinopecten sp. Species 0.000 abstract 1
- 235000020637 scallop Nutrition 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 8
- 230000008901 benefit Effects 0.000 description 3
- 239000002131 composite material Substances 0.000 description 3
- 238000005299 abrasion Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000004880 explosion Methods 0.000 description 2
- 230000004075 alteration Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009194 climbing Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 239000008358 core component Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Numerical Control (AREA)
- Milling Processes (AREA)
Abstract
本发明提供了一种自由曲面叶片平底刀三坐标端铣加工方法,首先将叶片实体模型数据输入到计算机内;接着以双三次非均匀有理B样条矩阵形式定义叶片的自由曲面为S;计算出叶片自由曲面的法向量;定义三坐标刀轴矢量为[0,0,1]T,计算自由曲面的最陡曲率,得到最陡曲线群;结合最陡曲线群及叶轮叶片流线特征,计算出分割线将自由曲面叶片分割成为多个子曲面;确定工艺粗糙度(残留高度),结合各子曲面扭曲及流线特征,以中间子曲面作为基准面,遵循沿加工方向各曲面相邻处刀位轨迹数量相同及加工方向一致或者相反的原则确定各子曲面的刀位轨迹及加工方向;确定子曲面分块加工顺序并进行后置处理,在三轴数控机床上进行加工。
Description
技术领域
本发明属于机械加工领域,涉及一种自由曲面叶片平底刀三坐标端铣加工方法。
背景技术
自由曲面叶片作为航空叶盘、船舶螺旋桨、高性能透平叶轮等的核心部件,广泛应用于航空航天、汽车船舶、石油化工、模具加工等领域。鉴于其叶型的空间任意扭曲性,目前普遍采用多轴数控铣削技术。传统自由曲面叶轮叶片的数控加工,一般采用球头刀点铣加工,用点、线切削来拟合出所需要的复杂叶片曲面,这种加工方法能够比较好的拟合出所要加工的复杂叶片曲面,但缺点是加工表面的一致性差、表面硬化严重,同时由于加工路径密集造成刀具磨损较快、加工效率低、加工周期长;同时,由于自由曲面叶片具有的空间任意扭曲性,绝大多数数控加工采用五坐标数控机床来完成,加工的设备成本高,也不利于高效精密数控加工的实现。
文献“Chen,Z.Z.C.,Fu,Q.,2007.A practical approach to generating steepestascent tool-paths for three-axis finish milling of compound NURBS surfaces.Comput.Aided.Des.39(11),964-974.”尝试了一种基于空间曲面零件三坐标加工的最陡爬坡法加工轨迹规划方法,该方法首先完成自由曲面的造型工作并参数化,计算最陡曲线,然后进行三坐标数控加工。但是,该方法加工路径的计算完全依赖工件曲面的空间扭曲状态,加工轨迹比较繁杂混乱,不适合叶轮机械对叶片残留刻痕尽量与叶轮流线一致的要求;并且,该方法只适合球头刀对加工,加工路径密集、效率低、加工周期长。
发明内容
本发明所要解决的技术问题是提供一种自由曲面叶片平底刀三坐标端铣加工方法,该方法基于实体模型规划平底刀端铣自由曲面叶片三坐标数控加工刀位轨迹,可实现自由曲面叶片等零件的高效精密数控加工。
本发明自由曲面叶片平底刀三坐标端铣加工方法,包括以下步骤:
1)将叶片实体模型数据输入到计算机内,完成目标自由曲面的造型;
2)以双三次非均匀有理B样条矩阵形式定义叶片的自由曲面为S;
3)计算出叶片自由曲面的法向量;
4)定义三坐标刀轴矢量为[0,0,1]T,计算自由曲面的最陡曲率,得到最陡曲线群;
5)结合最陡曲线群及叶片流线特征,计算出分割线将自由曲面叶片分割成为若干个子曲面;
6)确定工艺粗糙度(残留高度),结合各子曲面扭曲及流线特征,以中间子曲面作为基准面,遵循沿加工方向各曲面相邻处刀位轨迹数量相同及加工方向一致或者相反的原则确定各子曲面的刀位轨迹及加工方向;
7)确定子曲面分块加工顺序并进行后置处理,在三轴数控机床上进行加工。
所述步骤2)中,所述自由曲面S的表达式为:
其中,U、V为空间变量参数,从进口边到出口边为u向、叶片顶部到底部方向为v向;W为控制顶点对应的权因子;X、Y和Z为空间加权控制顶点在X轴、Y轴、Z轴上的投影矩阵:
U=[u0,u1,...,um+4];
V=[v0,v1,…,vn+4];
其中(x(i,j),y(i,j),z(i,j))为控制顶点,i∈(0,m+4);j∈(0,n+4),m、n分别为在U、V方向上组成自由曲面的点的个数;
所述叶片自由曲面的法向量为:
其中,
所述叶片自由曲面的最陡曲率为
G=[0,0,1]T-([0,0,1]·n)·n=(j1 2+j2 2+j3 2)-1·[-j1·j3,-j2·j3,j1 2+j2 2]T;
所述在计算出分割线时,包括以下步骤:
1)在最陡曲线群中,保留与叶片流线基本一致或垂直的最陡曲线,其中,与叶片流线基本一致的最陡曲线为第一最陡曲线,与叶片流线垂直的最陡曲线为第四最陡曲线;
2)在V方向对原叶片自由曲面加密,形成密集流线群,然后寻找与第一最陡曲线最接近的曲线作为第一分割线,其中,U方向长度由第一最陡曲线约束;
3)固定第一分割线尾端在V方向的j值,在叶片自由曲面上确立第二分割线;
4)计算V方向的各个j值的平均值,固定此平均值,依此确定第三分割线;所述基准面刀位轨迹线数量的计算方法为:
根据公式计算平底刀加工带宽,其中,L为空间加工带宽,d为带圆角平底刀底部平面段半径,r为刀具边缘处圆弧半径,h为已知残留高度,然后计算基准自由曲面在投影面的总宽度,结合以上计算的平底刀加工带宽,得到基准面刀位轨迹线数量;
其余子曲面刀位轨迹线数量的计算方法为:根据刀位轨迹线数量,根据公式计算其余各子曲面的残留高度Δh,如果子曲面所有残留高度Δh≤0.025mm,则保持基准面刀位轨迹数量,如果子曲面有Δh≥0.025mm情况,添加刀位轨迹以达到Δh≤0.025mm。
本发明自由曲面叶片平底刀三坐标端铣加工方法至少具有以下优点:本发明采用平底刀代替球头刀精铣自由曲面叶片,在同等精度要求情况下提高数控加工效率40%以上,并且在端铣过程中,平底刀接近刀具底部切削处的容屑、排屑和切削能力都优于球头刀,很大程度上解决了球头刀端铣容屑空间小、刀具易磨损、工件表面质量难以提高等问题,加工轨迹也保持规则有序;本发明同时用三坐标数控机床代替大部分五坐标数控机床进行自由曲面的数控加工,便于类似自由曲面零件等数控加工的推广,可大量节约加工设备成本,具有很高的社会效益及推广价值。
附图说明
图1为某透平叶轮自由曲面叶片结构图;
图2为待加工叶片曲面投影图;
图3为最陡曲线示意图;
图4为分割线示意图;
图5为自由曲面叶片分割成4个子曲面的爆炸效果图;
图6为带圆角平底刀切削示意图;
图7为球头刀加工带宽计算示意图;
图8为自由曲面叶片分块加工刀位轨迹规划示意图。
具体实施方式
下面以某透平叶轮自由曲面叶片模具平底刀三坐标端铣加工为例,结合附图,对本发明自由曲面叶片平底刀三坐标端铣加工方法做详细描述:
1.建立模型
请参见图1所示,将叶片实体模型数据输入计算机,完成目标自由曲面的造型。
2.确定自由曲面空间位置、初步分析
3.计算最陡曲线
(1)用双三次非均匀有理B样条(NURBS)矩阵形式定义叶片自由曲面为S:则有
定义U、V为空间变量参数(从进口边到出口边为u向、叶片顶部到底部方向为v向);W为控制顶点对应的权因子;X、Y和Z为空间加权控制顶点在X轴、Y轴、Z轴上的投影矩阵:
U=[u0,u1,...,um+4];
V=[v0,v1,…,vn+4];
其中(x(i,j),y(i,j),z(i,j))为控制顶点,i∈(0,m+4);j∈(0,n+4),m、n分别为在U、V方向上组成自由曲面的点的个数。
(2)计算出叶片自由曲面的法向量
其中
(3)定义三坐标刀轴矢量为[0,0,1]T,则最陡曲率为
G=[0,0,1]T-([0,0,1]·n)·n=(j1 2+j2 2+j3 2)-1·[-j1·j3,-j2·j3,j1 2+j2 2]T
由此得到最陡曲线(4条),参见图3。
4.分析、选择最陡曲线。
从保证叶轮机械流道内部流场的可靠性及避免“涡”的出现,在叶片的加工过程要求叶片最终的成型残留刻痕尽可能的保证和流线形态一致(参见图1,叶片存在7条从进口边到出口边的流线)。
(1)对于最陡曲线1:由于其形态与流线基本一致,保留;
(2)对于最陡曲线2、3,与流线相差太多,去除;
(3)对于最陡曲线4,虽然与流线相差最多,但由于其形态与流线近似垂直,保留。
5.确立分割线
(1)在V方向对原叶片曲面加密,形成密集流线群,寻求与最陡曲线1最接近的一条曲线确立为分割线1,U方向长度也由最陡曲线1长度约束;
(2)固定分割线1尾端在V方向的j(j∈n+4)值,在叶片曲面依此确立分割线2;
(3)分析最陡曲线4,计算其V方向的各j(j∈n+4)值的平均值,固定此平均值,依此确立分割线3;(参见图4)
6.分割曲面
以分割曲线1、2、3为基准,将原叶片自由曲面分割成为4个子曲面,图5为叶片曲面分割后的爆炸图效果;
7.确定工艺参数
确定加工工艺要求残留误差h为0.025mm,加工刀具为φ20(mm)直柄带圆弧R10(mm)平底刀;
8.确定基准面、刀位轨迹及加工方向
(1)以中间子曲面为基准面;
(2)刀位轨迹线的排列遵循加工刀位轨迹线形态与气流通道流线“一致”原则;
(3)加工方向遵循上坡(拖刀)原则,与流线方向一样(从进口边到出口边)或者相反(从出口边到进口边);
(4)刀位轨迹线即可用以下公式表示:
P(u)=S(u,v(u))=[x(u,v(u)),y(u,v(u)),z(u,v(u))]T
9.计算基准面刀位轨迹线数量
(1)图6为加工所用带圆角平底刀切削示意图,残留刻痕的特点是形成“窄又高”和“宽又平”相间隔排布,与球头刀加工相比较,明显的优势是切削面积增大、效率提高,并且由于刀具平底部分的残留高度(“宽又平”)明显要低于刀具边缘部分(“窄又高”),所以,在计算过程中可以很简单的套用球头刀残留高度与加工带宽(两条刀位轨迹线之间的间隔距离,依此来确定刀位轨迹线数量);
(2)图7为球头刀加工带宽计算示意图,公式为:
其中,L为空间加工带宽,r为刀具边缘处圆弧半径,h为已知残留高度,θ为叶片自由曲面与投影面的夹角,δ为投影面加工带宽,本发明加工带宽以L为准,不涉及θ、δ。
(3)由于本发明使用带圆角的平底刀加工,加工带宽计算公式在球头刀基础上改进为:
其中d为带圆角平底刀底部平面段半径。
(4)计算基准自由曲面在投影面的总宽度,结合以上所求出的加工宽度,得到基准面刀位轨迹线数量;
10.计算其余子曲面刀位轨迹线数量
(1)在刀位轨迹线数量(加工带宽L′)已知的情况下,反算其余各子曲面的残留高度Δh,公式为:
(2)分析除基准面之外的三个子曲面残留高度Δh分布情况:
a)如果子曲面所有残留高度Δh≤0.025mm,则保持步骤9结果;
b)如果子曲面有Δh≥0.025mm情况,添加刀位轨迹以达到Δh≤0.025mm;
c)图8为步骤9、10的刀位轨迹计算示意图;
11.定义步骤9、10所得到的刀位轨迹线为加工切触点,计算得到平底刀刀具底部中心轨迹线并针对所使用数控机床对其进行后置处理;
12.统一规划,联接、修正各子曲面刀位轨迹并确定子曲面分块加工顺序等,在三坐标数控机床实践加工为加工成型自由曲面叶片。
经过实践加工并与球头刀端铣比较,在残留高度一致的情况下,提高加工效率40%以上。同时,用三坐标机床代替传统的五坐标数控机床,大量的节约了加工设备成本。
以上所述仅为本发明的一种实施方式,不是全部或唯一的实施方式,本领域普通技术人员通过阅读本发明说明书而对本发明技术方案采取的任何等效的变换,均为本发明的权利要求所涵盖。
Claims (7)
1.自由曲面叶片平底刀三坐标端铣加工方法,其特征在于:包括以下步骤:
1)将叶片实体模型数据输入到计算机内,完成目标自由曲面的造型;
2)以双三次非均匀有理B样条矩阵形式定义叶片的自由曲面为S;
3)计算出叶片自由曲面的法向量;
4)定义三坐标刀轴矢量为[0,0,1]T,计算自由曲面的最陡曲率,得到最陡曲线群;
5)结合最陡曲线群及叶片流线特征,计算出分割线将自由曲面叶片分割成为若干个子曲面;
6)确定工艺粗糙度,结合各子曲面扭曲及流线特征,以中间子曲面作为基准面,遵循沿加工方向各曲面相邻处刀位轨迹数量相同及加工方向一致或者相反的原则确定各子曲面的刀位轨迹及加工方向;
7)确定子曲面分块加工顺序并进行后置处理,在三轴数控机床上进行加工。
3.如权利要求1中所述的自由曲面叶片平底刀三坐标端铣加工方法,其特征在于:所述叶片自由曲面的法向量为:
其中,
4.根据权利要求3所述的自由曲面叶片平底刀三坐标端铣加工方法,其特征在于:所述叶片自由曲面的最陡曲率为
G=[0,0,1]T-([0,0,1]·n)·n=(j1 2+j2 2+j3 2)-1·[-j1·j3,-j2·j3,j1 2+j2 2]T。
5.根据权利要求1所述的自由曲面叶片平底刀三坐标端铣加工方法,其特征在于:在计算出分割线时,包括以下步骤:
1)在最陡曲线群中,保留与叶片流线基本一致或垂直的最陡曲线,其中,与叶片流线基本一致的最陡曲线为第一最陡曲线,与叶片流线垂直的最陡曲线为第四最陡曲线;
2)在V方向对原叶片自由曲面加密,形成密集流线群,然后寻找与第一最陡曲线最接近的曲线作为第一分割线,其中,U方向长度由第一最陡曲线约束;
3)固定第一分割线尾端在V方向的j值,在叶片自由曲面上确立第二分割线;
4)计算V方向的各个j值的平均值,固定此平均值,依此确定第三分割线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110201226 CN102357666B (zh) | 2011-07-18 | 2011-07-18 | 自由曲面叶片平底刀三坐标端铣加工方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110201226 CN102357666B (zh) | 2011-07-18 | 2011-07-18 | 自由曲面叶片平底刀三坐标端铣加工方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102357666A CN102357666A (zh) | 2012-02-22 |
CN102357666B true CN102357666B (zh) | 2013-04-17 |
Family
ID=45583078
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110201226 Expired - Fee Related CN102357666B (zh) | 2011-07-18 | 2011-07-18 | 自由曲面叶片平底刀三坐标端铣加工方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102357666B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102649178B (zh) * | 2012-05-17 | 2013-07-31 | 西安交通大学 | 自由曲面离心压气机叶轮五坐标整体铣高效粗加工方法 |
CN103586737B (zh) * | 2013-09-30 | 2016-01-20 | 沈阳黎明航空发动机(集团)有限责任公司 | 一种叶片型面高精度数控铣加工的补偿方法 |
CN103927426B (zh) * | 2014-05-06 | 2017-07-28 | 南京航空航天大学 | 航空发动机整体叶盘特征识别方法 |
CN103934726B (zh) * | 2014-05-13 | 2016-05-11 | 航天科工哈尔滨风华有限公司 | 一种封闭式非直纹曲面数控加工刀路轨迹规划方法 |
CN104907617B (zh) * | 2015-06-15 | 2017-04-26 | 西安交通大学 | 基于分区域切削的离心压缩机叶轮五坐标铣削方法 |
CN105022343B (zh) * | 2015-07-13 | 2017-12-15 | 上海交通大学 | 基于空间填充曲线的数控加工刀具轨迹生成方法 |
CN105955195B (zh) * | 2016-05-16 | 2018-05-22 | 哈尔滨理工大学 | 一种基于铣削力预测的曲面加工轨迹生成方法 |
CN113182565B (zh) * | 2021-03-30 | 2022-04-08 | 成都飞机工业(集团)有限责任公司 | 一种弱刚性型面分区域铣削方法 |
CN114799293B (zh) * | 2022-06-30 | 2023-03-21 | 中国空气动力研究与发展中心高速空气动力研究所 | 一种风洞复杂曲面收缩段的加工方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101062526A (zh) * | 2006-04-27 | 2007-10-31 | 中国兵器工业集团第七0研究所 | 涡轮增压器压气机叶轮的铣削加工方法 |
CN101249571A (zh) * | 2008-03-26 | 2008-08-27 | 西北工业大学 | 多曲面岛屿五轴螺旋加工方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR2828824B1 (fr) * | 2001-08-23 | 2003-12-05 | Snecma Moteurs | Procede de fabrication d'un disque aubage monobloc de rotor et disque correspondant |
JP4957102B2 (ja) * | 2006-07-14 | 2012-06-20 | 株式会社日立プラントテクノロジー | 加工方法及び動吸振器 |
JP2010184302A (ja) * | 2009-02-10 | 2010-08-26 | Toshiba Corp | 曲面形状部品の加工方法 |
-
2011
- 2011-07-18 CN CN 201110201226 patent/CN102357666B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101062526A (zh) * | 2006-04-27 | 2007-10-31 | 中国兵器工业集团第七0研究所 | 涡轮增压器压气机叶轮的铣削加工方法 |
CN101249571A (zh) * | 2008-03-26 | 2008-08-27 | 西北工业大学 | 多曲面岛屿五轴螺旋加工方法 |
Non-Patent Citations (2)
Title |
---|
复杂自由曲面数控规划及其数据管理技术研究;赵丽;《工程科技Ⅰ辑 》;20101215;全文 * |
赵丽.复杂自由曲面数控规划及其数据管理技术研究.《工程科技Ⅰ辑 》.2010, |
Also Published As
Publication number | Publication date |
---|---|
CN102357666A (zh) | 2012-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102357666B (zh) | 自由曲面叶片平底刀三坐标端铣加工方法 | |
CN102649178B (zh) | 自由曲面离心压气机叶轮五坐标整体铣高效粗加工方法 | |
CN109702567B (zh) | 一种圆弧头立铣刀圆弧刃前刀面的磨削轨迹求解方法 | |
CN105955195B (zh) | 一种基于铣削力预测的曲面加工轨迹生成方法 | |
Lasemi et al. | Recent development in CNC machining of freeform surfaces: A state-of-the-art review | |
CN100585523C (zh) | 直纹面叶轮刀具轨迹规划加工方法 | |
CN111413923A (zh) | 一种用于复杂曲面加工的高速精密加工系统及方法 | |
CN102637216B (zh) | 一种复杂组合曲面的数控侧铣加工刀轨生成方法 | |
CN104907617B (zh) | 基于分区域切削的离心压缩机叶轮五坐标铣削方法 | |
Chaves-Jacob et al. | Optimal strategy for finishing impeller blades using 5-axis machining | |
CN102794488B (zh) | 类直纹面整体叶轮曲面的侧铣加工方法 | |
CN106843140B (zh) | 一种闭式叶轮的精加工刀具轨迹规划方法 | |
CN102866671B (zh) | 大圆弧直纹面数控加工刀具轨迹规划方法 | |
TWI453078B (zh) | 一種五軸曲面側銑加工系統及其路徑規劃方法 | |
CN109343468B (zh) | 一种基于投影偏置的叶片多轴轨迹生成方法 | |
CN103777568B (zh) | 一种基于刃磨过程的整体式立铣刀容屑槽建模方法 | |
CN109304505A (zh) | 一种三元叶轮的粗铣加工方法 | |
CN104715110B (zh) | 一种精密复杂曲面零件型面磨料流等余量精准抛光模具设计方法 | |
CN115616983A (zh) | 一种五轴联动同步刀具路径插补方法和系统 | |
Wei et al. | Modeling and machining of integral impeller based on NURBS curve | |
CN109590559B (zh) | 一种五轴电火花加工闭式叶轮轨迹规划方法 | |
Chuang et al. | Integrated rough machining methodology for centrifugal impeller manufacturing | |
CN113065205B (zh) | 一种采用平行砂轮磨削圆弧头后刀面的轨迹求解方法 | |
CN110109415A (zh) | 一种基于密度聚类的多网格刀轴优化方法 | |
Dong et al. | An improved calculation method for cutting contact point and tool orientation analysis according to the CC points |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130417 Termination date: 20150718 |
|
EXPY | Termination of patent right or utility model |