CN102581364B - 计算球头铣刀铣削负载的方法 - Google Patents

计算球头铣刀铣削负载的方法 Download PDF

Info

Publication number
CN102581364B
CN102581364B CN201210055859.3A CN201210055859A CN102581364B CN 102581364 B CN102581364 B CN 102581364B CN 201210055859 A CN201210055859 A CN 201210055859A CN 102581364 B CN102581364 B CN 102581364B
Authority
CN
China
Prior art keywords
cutter
milling
infinitesimal
cutting
workpiece
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
CN201210055859.3A
Other languages
English (en)
Other versions
CN102581364A (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.)
HUST Wuxi Research Institute
Original Assignee
Huazhong University of Science and Technology
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 Huazhong University of Science and Technology filed Critical Huazhong University of Science and Technology
Priority to CN201210055859.3A priority Critical patent/CN102581364B/zh
Publication of CN102581364A publication Critical patent/CN102581364A/zh
Application granted granted Critical
Publication of CN102581364B publication Critical patent/CN102581364B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Numerical Control (AREA)

Abstract

本发明公开了一种计算球头铣刀铣削负载的方法,包括如下步骤:获得刀具与工件的相对位置参数模型,根据刀具与工件的相对位置参数模型确定参与切削的全部刀具微元,计算每个参与切削的侧倾铣削微元的切入角、切出角以及瞬时切削厚度,计算每个参与切削的前倾铣削微元的切入角、切出角以及瞬时切削厚度,将每个侧倾铣削微元的瞬时切削厚度和每个参与切削的前倾铣削微元的瞬时切削厚度进行叠加,以获得刀具微元的瞬时切削厚度,将参与切削的所有刀具微元的瞬时切削厚度进行求和,以获取球头铣刀的瞬时铣削负载。本发明能在刀具侧倾和前倾铣削时获得刀具的切削刃与工件的瞬时切削状态以及切削厚度,从而实现五轴铣削加工中切削力的预测。

Description

计算球头铣刀铣削负载的方法
技术领域
本发明涉及计算机数控加工技术领域,具体涉及一种五轴数控加工中计算球头铣刀铣削负载的方法。
背景技术
球头铣刀五轴铣削是航空发动机叶轮、叶片类空间曲面和模具类复杂曲面精加工的主要技术手段。由于零件的复杂性,通常需要改变刀具的轴线方向从而避免干涉。为预测五轴铣削过程的切削力,减小加工过程中的变形和振动,首先需要计算铣刀切削工件的瞬时切削厚度。
球头铣刀铣削力建模一般将刀具沿轴线方向分割为若干圆盘状的刀具微元,通过计算微元的切削状态和微元切削力,再将微元的切削力进行叠加而求出整体的铣削力。然而现有的模型一般都是在刀具与工件的相对位置比较特殊(只有侧倾或只有前倾)的情况下进行铣削力预测,并采用矢量计算方法获得切削负载。然而,这种方法无法反映微元在一个周期内切入工件和切出工件的真实过程,也无法计算球头铣刀既侧倾又前倾的铣削负载。
发明内容
本发明的目的在于针对现有铣削力建模中存在的缺陷,提供一种计算球头铣刀铣削负载的方法,其能在刀具侧倾和前倾铣削时获得刀具的切削刃与工件的瞬时切削状态以及切削厚度,从而实现五轴铣削加工中切削力的预测。
本发明是通过以下技术方案实现的:
一种计算球头铣刀铣削负载的方法,包括如下步骤:
(1)获得刀具与工件的相对位置参数模型;
(2)根据刀具与工件的相对位置参数模型确定参与切削的全部刀具微元,刀具微元包括侧倾铣削微元与前倾铣削微元;
(3)根据刀具与工件的相对位置参数模型计算每个参与切削的侧倾铣削微元的切入角、切出角以及瞬时切削厚度;
(4)根据刀具与工件的相对位置参数模型计算每个参与切削的前倾铣削微元的切入角、切出角以及瞬时切削厚度;
(5)将每个侧倾铣削微元的瞬时切削厚度和每个参与切削的前倾铣削微元的瞬时切削厚度进行叠加,以获得刀具微元的瞬时切削厚度;
(6)将参与切削的所有刀具微元的瞬时切削厚度进行求和,以获取球头铣刀的瞬时铣削负载。
本发明具有以下的优点和技术效果:通过将球头铣刀切削过程进行分解,揭示刀具不同轴向位置的刀具微元切入与切出工件的过程,以及微元瞬时切削的厚度,获得瞬时参与切削的切削点,从而计算刀具倾斜情况下的切削厚度,以预测五轴加工瞬时切削力。
附图说明
图1为本发明计算球头铣刀铣削负载的方法的流程图
图2为刀具与工件的相对空间位置示意图。
图3(a)至图3(c)为刀具球面参与切削的部分与进给方向的关系。
图4为刀具微元进给方向的分解示意图。
图5为侧倾铣削示意图。
图6(a)为前倾铣削示意图。
图6(b)至图6(e)为前倾铣削微元的不同铣削状态示意图。
图7(a)和(b)示出采用本发明方法计算得到的铣削负载结果。
具体实施方式
以下结合附图和实施例对本发明的接受方案作进一步的详细描述,但本实施例并不用于限制本发明。
本发明采用的球头铣刀的参数如下:直径10毫米,刃数为4,螺旋角30度,被加工表面为平面,切削方式采用铣槽方式,切削深度为1毫米,每齿进给为0.1毫米。
如图1所示,本发明计算球头铣刀铣削负载的方法包括以下步骤:
1、获得刀具与工件的相对位置参数模型;
如图2所示,T为刀具轴线方向,切削进给方向与X轴方向相同,刀具的倾斜方向和角度见右图,O为刀具的球心点,OA表示刀具轴线,α表示刀具轴线与加工表面法向所成的角度,φ表示刀具微元的转角,即刀具轴线在加工表面投影与进给方向X的夹角,A′表示A点在加工表面的投影(如图3(a)至(c)所示),当φ=0时,为前倾铣削;当φ=π/2时,为侧倾铣削。
2、根据刀具与工件的相对位置参数模型确定参与切削的全部刀具微元,刀具微元包括侧倾铣削微元和前倾铣削微元,并且刀具微元的进给运动可分解为侧倾铣削微元进给运动与前倾铣削微元进给运动;
刀具微元的进给运动如图4所示,刀具微元的进给可以分为侧倾铣削进给fc和前倾铣削进给fq,并由下式(1)计算:
f c = f sin φ f q = f cos φ - - - ( 1 )
式中f为球头铣刀的进给。
如图5所示,θ表示刀具微元的位置角,设参与切削的区域最小位置角为θs,最大的为θe,则
θ s = α - arccos | OH | R θ e = α + arccos | OH | R - - - ( 2 )
当θs<θ<θe时,则微元参与切削。其中R为球头铣刀的半径,|OH|表示球心O到垂足H的线段长度。
3、根据刀具与工件的相对位置参数模型计算每个参与切削的侧倾铣削微元的切入角、切出角以及瞬时切削厚度;
侧倾铣削截面如图5所示,进给方向垂直于纸面向外。对于参与切削的不同微元,其切入切出角是不相同的。位置角为θ处的刀具微元,其切削半径为|AD|,切出角Φex,切入角Φst表达式为:
Φ st = arccos ( - | AC | | AD | ) Φ ex = 2 π - Φ st - - - ( 3 )
式中,|AC|和|AD|的值可通过相似ΔBAC和ΔBHO:
| AC | = | OH | · | AB | | HB | = | OH | · | OB | - | OA | | OH · tan α | = | OH · sec α | - | R · cos θ | tan α | AD | = R · sin θ - - - ( 4 )
在侧倾铣削中,由于刀具轴线垂直于进给方向,当刀具微元的转角满足时,其切削厚度表达式为:
由于变化中可取正值和负值,所以hc可能为负值,此时表示该刀刃远离工件切削表面。
4、根据刀具与工件的相对位置参数模型计算每个参与切削的前倾铣削微元的切入角、切出角以及瞬时切削厚度;
如图6(a),在前倾铣削中,设当前参与切削的刀具微元的位置角为θ,而前一次切削该处的刀具微元位置角为θ′,垂直于刀轴方向的每齿进给量fθ可表示为:
±fθ=Rsinθ+fqcosα-Rsinθ′(6)
Rcosθ=Rcosθ′+fqsinα      (7)
即fθ与θ满足等式:
( R sin θ + f q cos α + ‾ f θ ) 2 + ( R cos θ - f q sin α ) 2 = R 2 - - - ( 8 )
式中fq为前倾铣削的每齿进给量,式(6)中,在深色区域取正号,在浅色区域取负号,在式(8)中,与此相反。
在前倾的铣削方式中,微元的瞬时未变形切削厚度如图6。图6(a)中虚线轮廓为一齿进给之前刀具的位置,实线为当前切削刀具的轮廓位置。前倾铣削时,刀具微元存在四个典型位置,如图6(a)中1、2、3、4所示的四个位置,图6(b)、(c)、(d)、(e)是A向的平面图。位置1,当前刀具的相应微元切削前,该处工件未被前一位置刀具切过,因此相应微元的切削状态如图6(b)所示,剖面线表示为为正切削;位置2,如图6(c),虚线表示被前一位置刀具所切,当前刀具的相应微元切削之前,该处工件已被前一位置刀具切过,但当前微元仍能切到工件表面,剖面线表示为正切削;位置3,相应微元的切削状态如图6(d)所示,当前相应微元切削之前,该处工件已被前一位置刀具切过,但当前微元不能切到工件表面,剖面线为正切削区域,刀齿切入和切出时刻为负切削区域;位置4时,如图6(e),当前刀具微元不能切到工件,刀齿处在负切削区域。
如图6(b)所示,在位置1处,刀具微元的位置角θ满足下式:
arccos R cos &theta; e + f q sin &alpha; R &le; &theta; < &theta; e
式中,θe由式(2)给出。
设点O2到水平线的距离为dθ,等同于式(4)中的|AC|:
d &theta; = | OH | &CenterDot; sec &alpha; - R cos &theta; tan &alpha;
切入切出角可表示为:
&Phi; st = &pi; - arccos d &theta; R sin &theta; &Phi; ex = &pi; + arccos d &theta; R sin &theta;
当刀具微元转角满足时,其切削厚度表达式为:
因为小于零,所以上式用正号。
同理可求出位置2、位置3和位置4刀具的切入切出角和瞬时切削厚度。
5、将每个侧倾铣削微元的瞬时切削厚度和每个参与切削的前倾铣削微元的瞬时切削厚度进行叠加,以获得刀具微元的瞬时切削厚度;
令α=45°,φ=45°,计算结果如图7(a),图中的曲线为刀具旋转一周一个刀刃的不同位置切削工件的状态,小于等于零表示没有切削工件,大于零表示切削工件的瞬时切削厚度。图中可见,微元的切削最大厚度发生在切入工件的瞬时,之后切削厚度随时间逐渐减小,直到切出工件。令α=45°,φ=10°,计算结果如图7(b),微元的切削最大厚度发生在切入工件后的某一时刻,同时可计算任意时刻刀刃上任意点的切削厚度。
6、将参与切削的所有刀具微元的瞬时切削厚度进行求和,以获取球头铣刀的瞬时铣削负载。

Claims (1)

1.一种计算球头铣刀五轴铣削负载的方法,其特征在于,包括如下步骤:
(1)获得刀具与工件的相对位置参数模型,在此步骤中,将球头铣刀的刀具球心点定义为O,刀具轴线定义为OA,刀具轴线OA在加工表面上的投影定义为OA′,刀具轴线与加工表面法向之间所成的角度表示为α,刀具轴线在加工表面上的所述投影OA′与进给方向之间的夹角表示为φ,以此方式来构建刀具与工件之间的相对位置参数模型;
(2)根据刀具与工件的上述相对位置参数模型来确定参与切削的全部刀具微元,所述刀具微元同时包括侧倾铣削微元与前倾铣削微元,由此将所有刀具微元的进给运动各自分解为由以下表达式(一)表示的侧倾铣削微元进给运动和前倾铣削微元进给运动:
    (一)
其中,fc表示侧倾铣削微元的进给运动量,fq表示前倾铣削微元的进给运动量,表示各个刀具微元的转角,也即其所对应的刀具轴线在加工表面上的投影与进给方向之间分别所形成的夹角,f表示球头铣刀的总进给;
(3)根据刀具与工件的上述相对位置参数模型计算所述每个参与切削的侧倾铣削微元的切入角和切出角,同时基于以下表达式(二)来计算得出各个侧倾铣削微元的瞬时切削厚度hc
    (二)
其中,fc表示各个侧倾铣削微元的进给运动量,其数值在步骤(2)中计算得出;
(4)根据所述刀具与工件的相对位置参数模型计算所述每个参与切削的前倾铣削微元的切入角和切出角,同时基于以下表达式(三)和(四)来计算得出各个前倾铣削微元的瞬时切削厚度hq
    (三)
d &theta; = | OH | &CenterDot; sec &alpha; - R cos &theta; tan &alpha; (四)
其中,R表示球头铣刀的半径,θ表示各个刀具微元的位置角,|OH|表示铣刀的刀具球心到其垂足之间的线段长度;
(5)将所述每个侧倾铣削微元的瞬时切削厚度和所述每个参与切削的前倾铣削微元的瞬时切削厚度进行叠加,以获得所述刀具微元的瞬时切削厚度;
(6)将参与切削的所有刀具微元的瞬时切削厚度进行求和,以获取球头铣刀的瞬时铣削负载。
CN201210055859.3A 2012-03-05 2012-03-05 计算球头铣刀铣削负载的方法 Active CN102581364B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210055859.3A CN102581364B (zh) 2012-03-05 2012-03-05 计算球头铣刀铣削负载的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210055859.3A CN102581364B (zh) 2012-03-05 2012-03-05 计算球头铣刀铣削负载的方法

Publications (2)

Publication Number Publication Date
CN102581364A CN102581364A (zh) 2012-07-18
CN102581364B true CN102581364B (zh) 2014-08-27

Family

ID=46470904

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210055859.3A Active CN102581364B (zh) 2012-03-05 2012-03-05 计算球头铣刀铣削负载的方法

Country Status (1)

Country Link
CN (1) CN102581364B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI562870B (en) * 2014-12-30 2016-12-21 Tongtai Machine & Tool Co Ltd Method for determining processing parameters of ball end mill in cutting process

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104462775B (zh) * 2014-11-12 2018-07-20 天津大学 一种五轴加工中心铣削力计算方法
CN104794337B (zh) * 2015-04-17 2017-10-20 北京航空航天大学 一种基于边界条件判断的正交车铣加工端面刃切削力建模方法
CN109656192B (zh) * 2018-12-18 2020-02-21 华中科技大学 一种基于球头铣刀铣削力的加工优化方法
CN111061218B (zh) * 2019-12-31 2021-07-27 华中科技大学 一种基于acis的复杂曲面加工切削力预测方法及设备

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5383131A (en) * 1992-07-14 1995-01-17 Mazda Motor Corporation Method of determining cutting trajectory in N/C machining system
US5827020A (en) * 1995-07-11 1998-10-27 Toshiba Kikai Kabushiki Kaisha Numerical control apparatus
US6824336B2 (en) * 2000-12-28 2004-11-30 Mazda Motor Corporation Method for controlling cutting machine
US6942436B2 (en) * 2002-03-11 2005-09-13 Yoshiaki Kakino NC program generating method, NC apparatus, computer memory product, and computer program product
CN101269423A (zh) * 2008-05-14 2008-09-24 北京航空航天大学 采用经线分划刀具的多坐标端铣加工刀位优化方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5383131A (en) * 1992-07-14 1995-01-17 Mazda Motor Corporation Method of determining cutting trajectory in N/C machining system
US5827020A (en) * 1995-07-11 1998-10-27 Toshiba Kikai Kabushiki Kaisha Numerical control apparatus
US6824336B2 (en) * 2000-12-28 2004-11-30 Mazda Motor Corporation Method for controlling cutting machine
US6942436B2 (en) * 2002-03-11 2005-09-13 Yoshiaki Kakino NC program generating method, NC apparatus, computer memory product, and computer program product
CN101269423A (zh) * 2008-05-14 2008-09-24 北京航空航天大学 采用经线分划刀具的多坐标端铣加工刀位优化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
富宏亚,张翔,韩振宇,孙雅洲.微径球头铣刀铣削力建模与仿真.《计算机集成制造***》.2011,第17卷(第7期),1448-1453. *
数控铣削加工物理建模与仿真;苑辉;《哈尔滨理工大学硕士学位论文》;20090331;8-16 *
苑辉.数控铣削加工物理建模与仿真.《哈尔滨理工大学硕士学位论文》.2009,8-16.

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI562870B (en) * 2014-12-30 2016-12-21 Tongtai Machine & Tool Co Ltd Method for determining processing parameters of ball end mill in cutting process

Also Published As

Publication number Publication date
CN102581364A (zh) 2012-07-18

Similar Documents

Publication Publication Date Title
CN102581364B (zh) 计算球头铣刀铣削负载的方法
CN102004818B (zh) 自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法
EP1868054B1 (en) A method for milling grooves in bored workpieces
CN102681488B (zh) 一种铣削加工工件表面形貌的建模方法
US9696707B2 (en) Method of controlling tool orientation and step-over distance in face milling of curvilinear surfaces
CN103042262B (zh) 槽特征内型转角一体化精加工方法
CN105234466A (zh) 一种槽特征侧铣加工刀轨生成方法
CN102637216A (zh) 一种复杂组合曲面的数控侧铣加工刀轨生成方法
CN113927378B (zh) 一种粗铣刀周齿分屑槽的数控磨削轨迹计算方法
CN105701323A (zh) 一种插铣加工切削力预测建模方法
CN103646141A (zh) 一种平底螺旋立铣刀正交车铣轴类零件的切削力建模方法
CN109909806A (zh) 一种提升机器人铣削稳定域的方法
CN102825308B (zh) 飞机结构件闭角清根方法
CN107335847A (zh) 一种切削效能约束刀具姿态的加工方法
CN104493305A (zh) 一种基于齿轮加工机床调整的车齿加工方法
CN107544433B (zh) 数控机床加工过程球头铣刀与工件接触区域半解析建模法
Zhaoyu et al. A variable-depth multi-layer five-axis trochoidal milling method for machining deep freeform 3D slots
CN105171077B (zh) 一种高精度控制表面粗糙度的复合斜面铣削加工方法
Gdula Adaptive method of 5-axis milling of sculptured surfaces elements with a curved line contour
CN105911957A (zh) 一种盘状铣刀展成法加工齿轮齿面残留高度的控制方法
CN111950100A (zh) 枞树型齿廓精拉刀的切削载荷空间重构方法
CN103042434A (zh) 超精密铣削加工表面形貌纹理控制方法
RU2429949C1 (ru) Способ обработки моноколес
CN105759717B (zh) 一种用于五轴数控加工的刀具轨迹防过切的方法
CN106077775B (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
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200615

Address after: No.329, Yanxin Road, Huishan Economic Development Zone, Huishan District, Wuxi City, Jiangsu Province

Patentee after: HUST-WUXI Research Institute

Address before: 430074 Hubei Province, Wuhan city Hongshan District Luoyu Road No. 1037

Patentee before: HUAZHONG University OF SCIENCE AND TECHNOLOGY