CN101491844A - 圆周铣削过程铣削力系数与刀具径向偏心标定方法 - Google Patents
圆周铣削过程铣削力系数与刀具径向偏心标定方法 Download PDFInfo
- Publication number
- CN101491844A CN101491844A CNA2009100209069A CN200910020906A CN101491844A CN 101491844 A CN101491844 A CN 101491844A CN A2009100209069 A CNA2009100209069 A CN A2009100209069A CN 200910020906 A CN200910020906 A CN 200910020906A CN 101491844 A CN101491844 A CN 101491844A
- Authority
- CN
- China
- Prior art keywords
- tan
- lambda
- rho
- cos
- milling force
- 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.)
- Granted
Links
Images
Abstract
本发明公开了一种圆周铣削过程铣削力系数与刀具径向偏心标定方法,其特点是将实测瞬时铣削力由笛卡尔坐标系(X-Y-Z)转换到局部坐标系(T-R-Z),然后通过数学推导建立局部坐标系下瞬时铣削力与瞬时未变形切屑厚度的关系,最后通过使用非线性最小二乘法实现对铣削力系数与刀具径向偏心标定。由于本发明方法无需进行参数研究,因此将铣削力系数计算工作量由10800~18000个迭代步降低到13~26个迭代步;将刀具径向偏心计算工作量由10800~18000个迭代步降低到18~33个迭代步,提高了工作效率,降低了试验成本。
Description
技术领域
本发明涉及一种参数标定方法,特别是圆周铣削过程铣削力系数与刀具径向偏心标定方法,适用于金属材料圆周铣削过程。
背景技术
圆周铣削主要是通过立铣刀在材料表面切除多余的材料层来获得理想的工件形状、尺寸以及表面光洁度的机械加工方法。铣削过程涉及到刀刃与切削层材料之间的剪切作用,前刀面与切屑之间、后刀面与已加工表面之间的摩擦作用,刀具磨损等物理现象。导致这些物理现象发生的功率消耗在力学上即表现为铣削力的作用,根据Merchant正交切削和斜角切削模型,铣削过程的铣削力常被表示为:“铣削力=铣削力系数×切屑载荷”的形式。使用该公式时,需要知道铣削力系数以及切屑载荷的大小。切屑载荷一般根据刀具偏心和预设定的进给率计算得到,但由于铣削过程中刀具和工件的动态耦合作用,刀具偏心很难通过试验直接测试得到。因而,为了实现对切削力的有效建模,对铣削力系数和刀具偏心的有效标定是关键。
文献1“E.Budak,Y.Altintas,E.J.A.Armarego,Prediction of milling force coefficients fromorthogonal cutting data,Transactions of the ASME Journal of Manufacturing Science andEngineering 118(1996)216-224.”公开了一种利用实测平均铣削力标定常值铣削力系数的方法,这种方法忽略了刀具偏心且需要做大量的铣削试验来保证标定精度,试验成本高。
参照附图,文献2“J.H.Ko,W.S.Yun,D.W.Cho,K.F.Ehmann,Development of a virtualmachining system-Part 1:approximation of the size effect for cutting force prediction.Journal ofMachine Tools and Manufacture 42(2002)1595-1605.”公开了一种利用实测瞬时铣削力标定瞬时铣削力系数和刀具径向偏心参数ρ和λ的方法,该方法的基本步骤如下:
(1)设置k=0,l=0。k、l表示迭代步。
(2)假设刀具径向偏心参数值ρk=kΔρ、λl=lΔλ,依据给定的切削参数和ρk、λk计算切屑载荷。Δρ和Δλ表示迭代步长。
(3)使用步骤(1)得到的切屑载荷和试验实测瞬时铣削力,计算与ρk、λk对应的瞬时铣削力系数,并计算所得到的瞬时铣削力系数的标准方差。
(4)设置k=k+1。
(5)如果ρk≤ρcri,重复步骤(2)到步骤(4);如果ρk>ρcri,设置k=0,l=l+1。ρcri表示刀具偏心量ρ可能达到的最大上限值。
(6)如果λl≤2π,重复步骤(2)到步骤(5);如果λl>2π,迭代结束。
(7)选择标准方差最小的瞬时铣削力系数作为标定得到的铣削力系数值,对应的ρk及λk作为标定得到的刀具径向偏心参数值。
从以上步骤可以看出,该方法是通过参数研究进行优选实现标定,这种方法需要扫描每一种可能的参数组合情况,计算工作量大,比如,若ρcri=50μm、Δρ=1μm、 按照上面的步骤,需要50×360=18000个迭代步才能完成标定。
发明内容
为了克服现有技术在进行铣削过程铣削力标定时计算工作量大的不足,本发明提供一种圆周铣削过程铣削力系数与刀具径向偏心标定方法,首先将实测瞬时铣削力由笛卡尔坐标系(X-Y-Z)转换到局部坐标系(T-R-Z),然后通过数学推导建立局部坐标系下瞬时铣削力与瞬时未变形切屑厚度的关系,最后通过使用非线性最小二乘法实现对铣削力系数与刀具径向偏心标定,无需进行参数研究,可以提高工作效率,降低试验成本。本发明主要针对表达形式为“铣削力=铣削力系数1×瞬时未变形切屑厚度-铣削力系数2×瞬时未变形切屑厚度×切削宽度”的铣削力模型展开,其中铣削力系数1用切向分量kT、径向分量kR以及轴向分量kZ表示;铣削力系数2用切向分量mT、径向分量mR以及轴向分量mZ表示。
本发明解决其技术问题所采用的技术方案:一种圆周铣削过程铣削力系数与刀具径向偏心标定方法,其特点是包括下述步骤:
(a)确定立铣刀的半径R、螺旋角β、刀齿数N,设定标定试验的切削参数单齿进给量f、轴向切削深度Rz、径向切削深度Rr;要求:Rz≤2mm,Rr=R,工件是长方体块;
(d)根据步骤(c)的结果,将与各个采样点对应的局部分量表示为如下形式:
(e)采用非线性最小二乘法求解由步骤(d)得到铣削力系数kT,mT,kR,mR,kZ,mZ以及ξ1,ξ2…,ξN;
(f)使用由步骤(e)得到ξ1,ξ2…,ξN,建立如下关系式
(g)使用非线性最小二乘法求解步骤(f),得到刀具偏心参数ρ和λ的值。
本发明的有益效果是:由于将实测瞬时铣削力由笛卡尔坐标系转换到局部坐标系,然后通过数学推导建立局部坐标系下瞬时铣削力与瞬时未变形切屑厚度的关系,最后通过使用非线性最小二乘法实现对铣削力系数与刀具径向偏心标定,无需进行参数研究,将铣削力系数计算工作量由10800~18000个迭代步降低到13~26个迭代步;将刀具径向偏心计算工作量由10800~18000个迭代步降低到18~33个迭代步,提高了工作效率,降低了试验成本。
下面结合附图和实施例对本发明作详细说明。
附图说明
附图是背景技术中刀具径向偏心参数ρ和λ的示意图。
图中,■-刀具旋转中心,●-刀具几何中心。
具体实施方式
实施例1:(1)选定半径R为8mm、螺旋角β为30度的三齿硬质合金立铣刀,在四坐标立铣床上对长方体块铝合金Al7050进行顺铣切削,刀具主轴转速为1200RPM,单齿进给量0.15mm/齿,轴向切削深度Rz等于2mm,径向切削深度Rr等于8mm。
(4)根据步骤(3)的结果,将与各个采样点对应的局部分量表示为如下形式:
ξi=ρcos[λ-tan(π/6)/8-2(i-1)π/3]-ρcos[λ-tan(π/6)/8-2(i-2)π/3],i=1,2,3。
(5)采用非线性最小二乘法求解由步骤(4)得到的等式中的kT,mT,kR,mR,kZ,mZ,ξ1,ξ2,ξ3。
(6)使用由步骤(5)得到ξ1,ξ2,ξ3建立如下关系式
(7)采用非线性最小二乘法求解由步骤(6)得到的等式中的ρ和λ。
由步骤(5)和步骤(7)即标定得到了铣削力系数和刀具偏心参数的值,其标定结果如下表所示。
ρ(μm) | λ(度) | kT(N/mm2) | mT | kR(N/mm2) | mR | kZ(N/mm2) | mZ |
32.30 | 91.95 | 894.867 | 0.0745 | 65.043 | 0.618 | 125.592 | 0.187 |
值得说明的是,步骤(5)通过13步迭代就收敛了,步骤(7)通过5步迭代也收敛了。如果采用参数研究,取ρcri=50μm、Δρ=1μm、 按照文献1的步骤,需要50×360=18000个迭代步才能完成标定。以上数据说明采用本发明的方法,大幅度降低了计算工作量,降低了试验成本。
实施例2:(1)选定半径为5mm、螺旋角β为30度的三齿高速钢立铣刀,在三坐标立铣床上对长方体块铝合金Al7050进行顺铣切削,刀具主轴转速为4000RPM,单齿进给量0.075mm/齿,轴向切削深度Rz等于2mm,径向切削深度Rr等于5mm。
(4)根据步骤(3)的结果,将与各个采样点对应的局部分量表示为如下形式:
ξi=ρcos[λ-tan(π/6)/8-2(i-1)π/3]-ρcos[λ-tan(π/6)/8-2(i-2)π/3],i=1,2,3。
(5)采用非线性最小二乘法求解由步骤(4)得到的等式中的kT,mT,kR,mR,kZ,mZ,ξ1,ξ2,ξ3。
(6)使用由步骤(5)得到ξ1,ξ2,ξ3建立如下关系式
(7)采用非线性最小二乘法求解由步骤(6)得到的等式中的ρ和λ。
由步骤(5)和步骤(7)即标定得到了铣削力系数和刀具偏心参数的值,其标定结果如下表所示。
ρ(μm) | λ(度) | kT(N/mm2) | mT | kR(N/mm2) | mR | kZ(N/mm2) | mZ |
26.98 | 69.86 | 113.512 | 0.686 | 13.873 | 0.945 | 39.034 | 0.607 |
值得说明的是,步骤(5)通过26步迭代就收敛了,步骤(7)通过7步迭代也收敛了。如果采用参数研究,取ρcri=30μm、Δρ=1μm、 按照文献1的步骤,需要30×360=10800个迭代步才能完成标定。以上数据说明采用本发明的方法,大幅度降低了计算工作量,降低了试验成本。
Claims (1)
1、一种圆周铣削过程铣削力系数与刀具径向偏心标定方法,其特征在于包括下述步骤:
(a)确定立铣刀的半径R、螺旋角β、刀齿数N,设定标定试验的切削参数单齿进给量f、轴向切削深度Rz、径向切削深度Rr;要求:Rz≤2mm,Rr=R,工件是长方体块;
(d)根据步骤(c)的结果,将与各个采样点对应的局部分量表示为如下形式:
(e)采用非线性最小二乘法求解由步骤(d)得到铣削力系数kT,mT,kR,mR,kZ,mZ以及ξ1,ξ2…,ξN;
(f)使用由步骤(e)得到ξ1,ξ2…,ξN,建立如下关系式
(g)使用非线性最小二乘法求解步骤(f),得到刀具偏心参数ρ和λ的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100209069A CN101491844B (zh) | 2009-01-15 | 2009-01-15 | 圆周铣削过程铣削力系数与刀具径向偏心标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100209069A CN101491844B (zh) | 2009-01-15 | 2009-01-15 | 圆周铣削过程铣削力系数与刀具径向偏心标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101491844A true CN101491844A (zh) | 2009-07-29 |
CN101491844B CN101491844B (zh) | 2010-08-18 |
Family
ID=40922733
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100209069A Expired - Fee Related CN101491844B (zh) | 2009-01-15 | 2009-01-15 | 圆周铣削过程铣削力系数与刀具径向偏心标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101491844B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103894884A (zh) * | 2014-03-21 | 2014-07-02 | 西北工业大学 | 铣刀刀尖频响函数的预测方法 |
CN104392090A (zh) * | 2014-09-26 | 2015-03-04 | 北京理工大学 | 铝合金材料端铣切削力及切削加工变形模型的构建方法 |
CN104476327A (zh) * | 2014-12-31 | 2015-04-01 | 无锡透平叶片有限公司 | 一种基于稳定性分析的立铣刀结构参数优化方法 |
CN105436821A (zh) * | 2015-11-11 | 2016-03-30 | 无锡清杨机械制造有限公司 | 一种铝合金加工方法 |
CN107065765A (zh) * | 2016-12-26 | 2017-08-18 | 西北工业大学 | 基于刀具轴线运动模型的切屑厚度计算方法 |
CN109834512A (zh) * | 2017-11-27 | 2019-06-04 | 上海电力学院 | 两齿螺旋立铣刀刀齿半径误差计算方法 |
CN110625437A (zh) * | 2019-09-04 | 2019-12-31 | 中国工程物理研究院材料研究所 | 一种铣削过程中刀具偏心参数标定方法 |
CN111428366A (zh) * | 2020-03-24 | 2020-07-17 | 湖北文理学院 | 三齿交错盘铣刀铣削力建模方法、设备、存储介质及装置 |
CN111618657A (zh) * | 2020-05-21 | 2020-09-04 | 西安交通大学 | 一种铣刀偏心参数的通用在机非接触标定方法 |
CN113553727A (zh) * | 2021-08-12 | 2021-10-26 | 东北大学 | 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 |
-
2009
- 2009-01-15 CN CN2009100209069A patent/CN101491844B/zh not_active Expired - Fee Related
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103894884A (zh) * | 2014-03-21 | 2014-07-02 | 西北工业大学 | 铣刀刀尖频响函数的预测方法 |
CN103894884B (zh) * | 2014-03-21 | 2016-04-27 | 西北工业大学 | 铣刀刀尖频响函数的预测方法 |
CN104392090A (zh) * | 2014-09-26 | 2015-03-04 | 北京理工大学 | 铝合金材料端铣切削力及切削加工变形模型的构建方法 |
CN104476327A (zh) * | 2014-12-31 | 2015-04-01 | 无锡透平叶片有限公司 | 一种基于稳定性分析的立铣刀结构参数优化方法 |
CN104476327B (zh) * | 2014-12-31 | 2017-02-01 | 无锡透平叶片有限公司 | 一种基于稳定性分析的立铣刀结构参数优化方法 |
CN105436821A (zh) * | 2015-11-11 | 2016-03-30 | 无锡清杨机械制造有限公司 | 一种铝合金加工方法 |
CN107065765A (zh) * | 2016-12-26 | 2017-08-18 | 西北工业大学 | 基于刀具轴线运动模型的切屑厚度计算方法 |
CN107065765B (zh) * | 2016-12-26 | 2019-04-19 | 西北工业大学 | 基于刀具轴线运动模型的切削厚度计算方法 |
CN109834512A (zh) * | 2017-11-27 | 2019-06-04 | 上海电力学院 | 两齿螺旋立铣刀刀齿半径误差计算方法 |
CN110625437A (zh) * | 2019-09-04 | 2019-12-31 | 中国工程物理研究院材料研究所 | 一种铣削过程中刀具偏心参数标定方法 |
CN111428366A (zh) * | 2020-03-24 | 2020-07-17 | 湖北文理学院 | 三齿交错盘铣刀铣削力建模方法、设备、存储介质及装置 |
CN111428366B (zh) * | 2020-03-24 | 2023-08-25 | 湖北文理学院 | 三齿交错盘铣刀铣削力建模方法、设备、存储介质及装置 |
CN111618657A (zh) * | 2020-05-21 | 2020-09-04 | 西安交通大学 | 一种铣刀偏心参数的通用在机非接触标定方法 |
CN113553727A (zh) * | 2021-08-12 | 2021-10-26 | 东北大学 | 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101491844B (zh) | 2010-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101491844B (zh) | 圆周铣削过程铣削力系数与刀具径向偏心标定方法 | |
CN103258095B (zh) | 平底立铣刀通用铣削力建模方法 | |
CN102136021B (zh) | 钛合金tc18铣削过程铣削力建模方法 | |
CN103567815A (zh) | 基于铣削小孔的数控机床切削热误差测试和评价方法 | |
CN101412122A (zh) | 液体火箭发动机喷管冷却通道立式加工方法 | |
CN101853324A (zh) | 圆周铣过程中铣削力建模方法 | |
CN102004818A (zh) | 自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法 | |
CN101840457A (zh) | 大尺寸枞树型叶根轮槽型面残余应力场重构的方法 | |
CN103433807A (zh) | 一种铣削力模型工艺参数的优化方法 | |
Daniyan et al. | MODELLING AND OPTIMIZATION OF THE CUTTING FORCES DURING TI6AL4V MILLING PROCESS USING THE RESPONSE SURFACE METHODOLOGY AND DYNAMOMETER. | |
CN109754332A (zh) | 基于切削力的机床铣削加工过程的能耗模型建模方法 | |
CN106863009B (zh) | 基于刀杆两点变形的切削力测量方法 | |
CN102059381A (zh) | 周铣过程中加工误差预测方法 | |
CN103358183A (zh) | 在线测量铣削刀具变形的装置及方法 | |
CN103105332B (zh) | 用正交切削法测定材料断裂韧性的方法 | |
CN205342042U (zh) | 一种可调式刮面刀具 | |
CN108647413B (zh) | 一种微细表面位置误差与稳定性综合预测方法 | |
Klocke et al. | New concepts of force measurement systems for specific machining processes in aeronautic industry | |
CN202726640U (zh) | 一种用于检测机床切削能力的测试试件 | |
CN105930654A (zh) | 一种适用于难加工材料高速车削过程的切削力建模方法 | |
CN104889821A (zh) | 模拟工况下的数控刀架可靠性切削试验方法 | |
CN101908084B (zh) | 1000mw超临界汽轮机转子轮槽面温度场重构的方法 | |
CN101549469B (zh) | 圆周铣削过程中刀具参数的标定方法 | |
Sheng | A modeling method for turning parameters coupling based on minimum cutting tool wear | |
Liu et al. | Geometry features of breakage section and variation of cutting force for end mills after brittle breakage |
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: 20100818 Termination date: 20130115 |
|
CF01 | Termination of patent right due to non-payment of annual fee |