CN104462775B - 一种五轴加工中心铣削力计算方法 - Google Patents
一种五轴加工中心铣削力计算方法 Download PDFInfo
- Publication number
- CN104462775B CN104462775B CN201410636591.1A CN201410636591A CN104462775B CN 104462775 B CN104462775 B CN 104462775B CN 201410636591 A CN201410636591 A CN 201410636591A CN 104462775 B CN104462775 B CN 104462775B
- Authority
- CN
- China
- Prior art keywords
- cutting
- indicate
- infinitesimal
- axis
- tool
- 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
Links
Abstract
本发明公开了一种五轴加工中心铣削力计算方法,所述方法包括以下步骤:计算未变形切削厚度;通过未变形切削厚度建立铣削力模型,通过铣削力模型计算五轴加工中心铣削力。本发明基于刀具切削刃微分化的思想,揭示了五轴加工中心刀具真实运动轨迹,以及微元瞬时未变形切削厚度,考虑到耕犁力对铣削力的影响,从而提高了计算五轴加工瞬时切削力的精度。
Description
技术领域
本发明涉及数控加工技术领域,尤其涉及一种五轴加工中心铣削力计算方法。
背景技术
在航空航天、汽车、船舶和各种高技术装备中,具有复杂空间型面几何特征的零件应用日益广泛,在实现系统力学性能、光学性能、流体性能等物理性能要求方面扮演重要角色。五轴联动数控技术在复杂型面的加工中具有独特的优势,即可以通过调整刀具的位姿使刀具保持最优的切削状态并避免刀具干涉,从而提高了零件的加工制造精度。铣削力为加工过程中的参数优化,加工精度的提升提供了重要的理论依据,因此针对复杂型面零件的铣削加工进行铣削力的精确计算具有十分重要的意义。
(1)曹清园对复杂型面加工中的铣削力和加工误差补偿进行了研究。建立了五轴机床球头铣刀曲面加工中的铣削力模型,基于此模型又提出了加工误差补偿的方法(参见基于铣削力建模的复杂曲面加工误差补偿研究[D].山东大学,2011)。
(2)Huang T,Zhang X和Ding H对铣削厚度进行解耦分析,建立了五轴加工中心铣削力模型(参见Huang T,Zhang X,Ding H(2013)Decoupled chip thickness calculationmodel for cutting force prediction in five-axis ball-end milling.TheInternational Journal of Advanced Manufacturing Technology 69(5-8):1203-1217)。
(3)Azeem A和Feng H-Y通过对刀具真实切削轨迹的分析,建立了铣削力的预测模型,并通过实验验证该模型的准确性(参见Azeem A,Feng H-Y(2013)Cutting forceprediction for ball-end mills with non-horizontal and rotational cuttingmotions.The International Journal of Advanced Manufacturing Technology 67(5-8):1833-1845)。
文献(1)(2)中所建立的铣削力模型,其中瞬时未变形切削厚度是以圆弧近似计算获得,然而此方法仅适用于刀具每齿进给量与刀具直径的比值极小时,对于实际的加工而言,所建立的力学模型并不精确。
文献(3)中提出了刀具类摆线运动的铣削力预测模型,分析了刀具的真实运动轨迹,但此模型仅适用于三轴机床。
发明内容
本发明旨在通过对五轴加工中心刀具真实运动轨迹的分析,考虑到耕犁力对铣削力的影响,提出一种能够快速准确计算五轴铣削力的计算方法,为高精密加工复杂曲面零件提供理论依据,详见下文描述:
一种五轴加工中心铣削力计算方法,所述方法包括以下步骤:
计算未变形切削厚度;
通过未变形切削厚度建立铣削力模型,通过铣削力模型计算五轴加工中心铣削力;
其中,所述未变形切削厚度为:
式中,和分别表示M,N点在工件坐标系下的三维坐标;M,N分别为第一、二条刀刃上某一微元点。
所述计算未变形切削厚度的步骤为:
建立刀具切削刃线的参数模型;根据五轴机床类型,建立刀具在切削工件时的轨迹方程;
采用分段三次NURBS(非均匀有理B样条曲线)曲线表示刀具切削刃,通过NURBS曲线和被切除材料实体之间的求交运算,获得刀刃微元在切削过程中的有效切削区间,计算刀具微元瞬时未变形切削厚度。
所述铣削力模型为:
在某一时刻,任意位置角作用在刀具上的切削合力通过对所有参与切削的切削微元的受力进行累加得到,其各向分力如下:
式中,分别表示刀具整体切削力在工件坐标系Xw,Yw,Zw三个方向上的分力,dz表示切削微元在Z方向的切削厚度;Nf表示刀具切削刃个数;表示刀具切削刃微元的有效切削范围,分别表示微元在工件坐标系Xw,Yw,Zw三个方向上的分力。
本发明提供的技术方案的有益效果是:本发明基于刀具切削刃微分化的思想,建立刀具切削刃线的参数模型。根据五轴机床类型,建立刀具在切削工件时的运动学模型。基于商用软件UG NX8.0的实体造型技术,用分段三次NURBS曲线表示刀具切削刃,通过NURBS曲线和被切除材料实体之间的求交运算,从而获得刀刃微元在切削过程中的有效切削区间,计算刀具微元瞬时未变形切削厚度。最后通过试切实验对切削力参数进行识别,由各微元局部切削力建立铣刀的整体铣削力预测模型。本发明基于刀具切削刃微分化的思想,揭示了五轴加工中心刀具真实运动轨迹,以及微元瞬时未变形切削厚度,考虑到耕犁力对铣削力的影响,从而提高了计算五轴加工瞬时切削力的精度。
附图说明
图1为球头铣刀切削刃几何模型及刀具受力分析示意图;
图2为未变形切削厚度示意图;
图3为采用本发明方法计算得到的铣削力与测量数据对比图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面对本发明实施方式作进一步地详细描述。
以下结合附图,以B摆头C转台五轴的联动机床球头铣削为例,详细说明本发明的具体实施。
101:计算未变形切削厚度;
通过分析切削刃线几何模型,发现在切削刃位置角一定的情况下,切削刃微元与螺旋滞后角、轴向位置角、水平位置角、轴向位置等参数是一一对应的,即根据其中任一参数就能确定切削刃微元位置。因此,建立刀具坐标系[Oc,Xc,Yc,Zc],Oc为坐标系原点位于刀具底端,以切削刃微元轴向位置z为参数,切削刃曲线描述为:
切削刃微元的切削半径可由式(2)计算获得。
公式(1)、(2)中,分别切削刃曲线上任意一点的坐标,R(z)表示切削刃微元在轴向位置z处的切削半径,R表示球头铣刀半径,表示第i号切削刃上微元水平位置角,θi表示第i号切削刃位置角,Nf表示刀具切削刃个数,ψ(z)表示微元相对于刀尖位置的螺旋滞后角,如图1所示。
在五轴机床加工过程中,刀具进行着绕刀轴的高速旋转和沿进给方向平动的复合运动与工件相对于机床的旋转运动,刀刃不断切入-切出工件材料形成切屑。针对B摆头C转台五轴联动机床进行运动学分析,建立切削刃在工件坐标系[Ow,Xw,Yw,Zw]下的轨迹方程,其中Ow为工件坐标系原点位于工件上表面中心。
第一条刀刃的轨迹方程:
第二条刀刃的轨迹方程:
其中T为坐标转换矩阵,将刀具坐标系转换到工件坐标系。
公式(3)、(4)中,和分别表示M,N点在工件坐标系下的三维坐标,zM,zN∈(0,R)分别表示第一、二条刀刃上某一微元点M,N在刀具坐标系下的轴向位置,fx,fy,fz分别表示刀具沿Xc,Yc,Zc方向的每齿进给量。θB表示刀具绕Yw轴的旋转角度,θC表示工件绕Zw轴的回转角度,L表示刀尖点与机床旋转轴B之间的距离,D表示刀尖初始状态下在工件坐标系下的轴向距离。
当R=12mm,f=0x,fy=0.1,fz=0,θB=π/6,θC=0时,得到如图2所示的刀具刃线轨迹。
基于微元刀刃切削过程为斜交切削的假设,微元刀刃对应的法向未变形切屑厚度可表示为:过球心点用与刀刃微元的线段与刀刃轨迹交点间的距离来度量。在实际加工过程中,刀具的中心也随之沿切削方向运动,假设切削至N点时球头铣刀球心位置为Ci,则M点的瞬时切削厚度为tn。空间直线MCi可由公式(5)获得,其中Ci点坐标由公式(6)计算获得,分别表示球心点Ci在工件坐标系中的三维坐标,k为比例系数。如图2所示。
基于商用软件UG NX8.0的实体造型技术表示工件、刀具、刀具扫描体以及被切除材料实体,用分段三次NURBS曲线表示刀具切削刃,通过NURBS曲线和被切除材料实体之间的求交运算,抽取参与切削的切削刃片段,从而获得刀刃微元在切削过程中的有效切削区间:切入和切出角轴向极限位置
由于刀具的真实运动轨迹为类摆线运动,将公式(5)、(6)带入公式(7),从而刀刃微元的瞬时未变形切削厚度tn可由M和N之间的距离表示。如图2所示。
102:通过未变形切削厚度建立铣削力模型,通过铣削力模型计算五轴加工中心铣削力。
球头铣刀加工复杂型面切削力模型,是将刀具切削刃沿轴向分割为许多切削微元,刀具受到的切削力为所有参加切削的微元受力之和。如图2所示。作用在切削微元上的空间铣削力可以分解为微切向力dFt、微径向力dFr和微轴向力dFa,其表达式分别为:
公式(8)中,tn表示瞬时切削厚度,Kte,Kre,和Kae表示切向、径向和轴向切削力系数,Ktc,Krc,和Kac表示切向、径向和轴向耕犁力系数,通过试切测量实验,按照转速-每齿进给量-切深-切宽设计切削力实验参数,采用非线性最小二乘的方法求解相应系数。dS,db分别为切削刃微元的切削宽度和长度。为了便于铣削力的分析和测量,微切向力dFt、微径向力dFr和微轴向力dFa分解到工件坐标系下,得到微元在工件坐标系Xw,Yw,Zw三个方向上的分力
其中,κ表示微元点在刀具坐标系下的位置角。
在某一时刻,即在任意位置角作用在刀具上的切削合力可通过对所有参与切削的切削微元的受力进行累加得到,其各向分力如下:
公式(10)中,分别表示刀具整体切削力在在工件坐标系Xw,Yw,Zw三个方向上的分力,dz表示切削微元在Z方向的切削厚度。
为了验证本方法的铣削力计算能力,对实际应用中的加工过程进行仿真。铣削条件如下:机床为DMU-60五轴加工中心,刀具为直径12mm的球头铣刀,切削刃数Nf=2,主轴转速n=1500r/min,每齿进给量Fz=0.2mm/min,B=-15°,C=5°,铣削类型为槽铣加工,切深ap=2mm。
图3为仿真实验在走刀过程中的铣削力F随刀具旋转角的变化图。由图3可见,仿真与实测的各铣削力分量比较吻合。通过该实验验证了本方法的可行性。
本领域技术人员可以理解附图只是一个优选实施例的示意图,上述本发明实施例序号仅仅为了描述,不代表实施例的优劣。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种五轴加工中心铣削力计算方法,其特征在于,所述方法包括以下步骤:
计算未变形切削厚度;
通过未变形切削厚度建立铣削力模型,通过铣削力模型计算五轴加工中心铣削力;
其中,未变形切削厚度定义为:过球心点与刀刃微元构成一线段,该线段与刀刃轨迹交于一点,交点与微元点之间的距离即为未变形切削厚度;
所述未变形切削厚度为:
式中,和分别表示M,N点在工件坐标系下的三维坐标;M,N分别为第一、二条刀刃上某一微元点;
所述铣削力模型为:
在某一时刻,任意位置角作用在刀具上的切削合力通过对所有参与切削的切削微元的受力进行累加得到,其各向分力如下:
式中,分别表示刀具整体切削力在工件坐标系Xw,Yw,Zw三个方向上的分力,dz表示切削微元在Z方向的切削厚度;Nf表示刀具切削刃个数;表示刀具切削刃微元的有效切削范围,分别表示微元在工件坐标系Xw,Yw,Zw三个方向上的分力;
第一条刀刃的轨迹方程:
第二条刀刃的轨迹方程:
其中T为坐标转换矩阵,将刀具坐标系转换到工件坐标系;
和分别表示M,N点在工件坐标系下的三维坐标,zM,zN∈(0,R)分别表示第一、二条刀刃上某一微元点M,N在刀具坐标系下的轴向位置,fx,fy,fz分别表示刀具沿Xc,Yc,Zc方向的每齿进给量,θB表示刀具绕Yw轴的旋转角度,θC表示工件绕Zw轴的回转角度,L表示刀尖点与机床旋转轴B之间的距离,D表示刀尖初始状态下在工件坐标系下的轴向距离;表示第i号切削刃上微元点M的水平位置角,表示第i+1号切削刃上微元点N的水平位置角,R表示球头铣刀半径;
R=12mm,fx=0,fy=0.1,fz=0,θB=π/6,θC=0。
2.根据权利要求1所述的一种五轴加工中心铣削力计算方法,其特征在于,所述计算未变形切削厚度的步骤为:
建立刀具切削刃线的参数模型;根据五轴机床类型,建立刀具在切削工件时的轨迹方程;
采用分段三次非均匀有理B样条曲线表示刀具切削刃,通过非均匀有理B样条曲线和被切除材料实体之间的求交运算,获得刀刃微元在切削过程中的有效切削区间,计算刀具微元瞬时未变形切削厚度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410636591.1A CN104462775B (zh) | 2014-11-12 | 2014-11-12 | 一种五轴加工中心铣削力计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410636591.1A CN104462775B (zh) | 2014-11-12 | 2014-11-12 | 一种五轴加工中心铣削力计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104462775A CN104462775A (zh) | 2015-03-25 |
CN104462775B true CN104462775B (zh) | 2018-07-20 |
Family
ID=52908804
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410636591.1A Active CN104462775B (zh) | 2014-11-12 | 2014-11-12 | 一种五轴加工中心铣削力计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104462775B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104699919A (zh) * | 2015-03-30 | 2015-06-10 | 济南大学 | 一种利用断续车削仿真多周期铣削的建模方法 |
CN104794337B (zh) * | 2015-04-17 | 2017-10-20 | 北京航空航天大学 | 一种基于边界条件判断的正交车铣加工端面刃切削力建模方法 |
CN105426697B (zh) * | 2015-12-30 | 2018-12-07 | 西安交通大学 | 一种螺旋铣刀五轴加工铣削力精确预测方法 |
CN106326647B (zh) * | 2016-08-18 | 2018-08-21 | 大连理工大学 | 碳纤维增强复合材料铣削力预测模型的建立方法 |
CN107728577B (zh) * | 2017-09-28 | 2019-09-27 | 大连理工大学 | 基于薄壁曲面加工变形的瞬时切削量规划方法 |
CN110989503A (zh) * | 2019-11-29 | 2020-04-10 | 华中科技大学 | 一种让刀误差约束复杂曲面铣削加工进给速度的控制方法 |
CN111428366B (zh) * | 2020-03-24 | 2023-08-25 | 湖北文理学院 | 三齿交错盘铣刀铣削力建模方法、设备、存储介质及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
HU4590U (en) * | 2011-05-09 | 2016-01-28 | Pécsi Tudományegyetem | Computer controlled milling machine for optimizing cutting forces in a milling process |
CN102581364B (zh) * | 2012-03-05 | 2014-08-27 | 华中科技大学 | 计算球头铣刀铣削负载的方法 |
-
2014
- 2014-11-12 CN CN201410636591.1A patent/CN104462775B/zh active Active
Non-Patent Citations (3)
Title |
---|
基于真实刀刃轨迹的立铣刀切削厚度模型;闫雪等;《机械工程学报》;20110105;第47卷(第1期);第182-186页 * |
复杂曲面高性能侧铣加工技术与方法研究;郭强;《中国博士学位论文全文数据库》;20130815(第08期);第27-98页 * |
多轴联动数控加工球头铣刀铣削力建模与仿真;丁云鹏;《中国优秀硕士学位论文全文数据库》;20140615(第06期);第18-32页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104462775A (zh) | 2015-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104462775B (zh) | 一种五轴加工中心铣削力计算方法 | |
Ming et al. | Four-axis trochoidal toolpath planning for rough milling of aero-engine blisks | |
CN105159228B (zh) | 五轴联动数控机床实现rtcp功能的五轴标定方法 | |
CN109158617A (zh) | 控制点驱动投影生成自由曲面车削刀具路径的方法 | |
CN110262397A (zh) | 车铣加工空间螺旋次摆线运动轨迹及瞬时切削力预测模型 | |
Huang et al. | Decoupled chip thickness calculation model for cutting force prediction in five-axis ball-end milling | |
CN101870073B (zh) | 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法 | |
CN109597357B (zh) | 一种面向叶片旋铣工艺的数控编程方法及装置 | |
CN104460516A (zh) | 一种基于后置处理五轴刀具半径补偿方法 | |
CN110032140B (zh) | 一种五轴加工中球形刀刀轴矢量规划方法 | |
CN109696884A (zh) | 一种非正交双转台型五轴机床后置处理方法 | |
Yu et al. | Post-processing algorithm of a five-axis machine tool with dual rotary tables based on the TCS method | |
CN112883505B (zh) | 考虑刀具工件相对振动的超精密端面车削表面建模方法 | |
Zhang et al. | A tool path generation method based on smooth machine rotary angle and tilt angle in five-axis surface machining with torus cutters | |
Tang et al. | Post-processor development for a turning and milling composite machine tool | |
Yang et al. | Research on multi-axis CNC programming in machining large hydraulic turbine's blades based on UG | |
Lu et al. | Smooth flank milling tool path generation for blade surfaces considering geometric constraints | |
Song et al. | Postprocessor algorithm and feedrate optimization for nine-axis milling machine tool with twin cutters | |
Qin et al. | Estimation and compensation for deviation of contour in slow tool servo precision turning for complicated curved surface | |
Shan et al. | Three half-axis tool orientation optimization for spiral machining of blades | |
Lin et al. | Trajectory analyses for five-axis machine tools | |
Deng et al. | The theoretical and experimental research of machining equal base circle bevel gears with a pot-shaped milling cutter | |
Cho et al. | Computation of optimal tool length for 5-axis ball-ended milling of molding die | |
CN105373662A (zh) | 一种三轴加工切屑几何与切削力集成建模方法 | |
Ma et al. | A non-singular five-axis trochoidal milling process method for 3D curved slots |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CP02 | Change in the address of a patent holder | ||
CP02 | Change in the address of a patent holder |
Address after: 300350 Haijing garden, Haihe Education Park, Jinnan, Tianjin, 135, Tianjin University. Patentee after: Tianjin University Address before: 300072 Tianjin City, Nankai District Wei Jin Road No. 92 Patentee before: Tianjin University |