CN101549469B - 圆周铣削过程中刀具参数的标定方法 - Google Patents
圆周铣削过程中刀具参数的标定方法 Download PDFInfo
- Publication number
- CN101549469B CN101549469B CN2009100211266A CN200910021126A CN101549469B CN 101549469 B CN101549469 B CN 101549469B CN 2009100211266 A CN2009100211266 A CN 2009100211266A CN 200910021126 A CN200910021126 A CN 200910021126A CN 101549469 B CN101549469 B CN 101549469B
- Authority
- CN
- China
- Prior art keywords
- cutter
- milling force
- milling
- value
- parameter
- 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
Images
Landscapes
- Numerical Control (AREA)
Abstract
本发明公开了一种圆周铣削过程中切削参数的标定方法,特别是圆周铣削过程中刀具偏心四自由度参数和瞬时铣削力系数标定方法。其特点是将实测瞬时铣削力由笛卡尔坐标系(X-Y-Z)转换到局部坐标系(T-R-Z),然后通过数学推导建立局部坐标系下瞬时铣削力与瞬时未变形切屑厚度的关系,最后通过使用单纯形法实现对铣削力系数与刀具径向偏心标定。本发明通过两次铣削试验即可完成标定,并同时考虑了铣削力系数的尺寸效应以及刀具偏心沿轴线的变化,无需进行大量的铣削试验,将标定瞬时铣削力系数的工作量由18000个迭代步降低到18个迭代步;将标定刀具径向偏心参数φ和τ的工作量由21600个迭代步降低到38个迭代步,提高了工作效率,降低了试验成本。
Description
技术领域
本发明涉及一种参数标定方法,特别是圆周铣削过程中刀具偏心四自由度参数和瞬时铣削力系数标定方法。
背景技术
参照图1~3。文献1“S.Y.Liang,J.-J.J.Wang,Milling force convolution modeling foridentification of cutter axis offset,International Journal of Machine Tools and Manufacture 34(1994)1177-1190.”公开了一种利用实测平均铣削力和铣削力谐波分量标定常值铣削力系数和刀具偏心两自由度参数ρ1和λ1的方法,这种方法既没考虑铣削力系数的尺寸效应也没有考虑偏心沿刀具轴线的变化。
文献2“J.H.Ko,W.S.Yun,D.W.Cho,K.F.Ehmann,Development of a virtual machiningsystem-Part 1:approximation of the size effect for cutting force prediction.Journal of MachineTools and Manufacture 42(2002)1595-1605.”公开了一种利用实测瞬时铣削力标定瞬时铣削力系数和刀具偏心两自由度参数ρ1和λ1的方法,这种方法基于参数研究的思路进行,考虑了铣削力系数的尺寸效应但没有考虑偏心沿刀具轴线的变化。
但是,这种方法很耗费时间。例如,ρ1从0.001mm以0.001mm为步长变化到0.050mm,λ1从1度以1度为步长变化到360度,预测铣削力与实测铣削力之间偏差最小的结果作为标定结果,需要经过50×360=18000步才能完成标定。
文献3“L.Zheng,S.Y.Liang,Identification of cutter axis tilt in end milling,Transactions ofthe ASME Journal of Manufacturing Science and Engineering 119(1997)178-185.”公开了一种利用实测平均铣削力和铣削力谐波分量标定常值铣削力系数和刀具偏心四自由度参数ρ、λ、φ和τ的方法,这种方法考虑了偏心沿刀具轴线的变化但忽略了铣削力系数的尺寸效应。
发明内容
为了克服现有技术在进行铣削力系数与刀具偏心参数标定时不能同时考虑铣削力系数的尺寸效应以及刀具偏心沿轴线变化以及标定过程耗时大的不足,本发明提供一种圆周铣削过程中切削参数的标定方法,通过将测试瞬时铣削力转换到局部坐标系,并采用单纯形法,可以实现快速对瞬时铣削力系数和刀具偏心参数进行标定,同时考虑铣削力系数的尺寸效应以及刀具偏心沿轴线的变化,无需进行参数研究,可以提高参数标定效率。
本发明解决其技术问题所采用的技术方案:一种刀具瞬时铣削力系数标定方法,其特点是包括下述步骤:
(a)设定立铣刀的半径R、螺旋角β、刀齿数N、刀具安装到主轴后的悬臂长度L;设定标定试验的切削参数:单齿进给量f、轴向切削深度Rz、径向切削深度Rr、刀具主轴转速;Rz≤2mm,Rr=R;工件是长方体块;
(e)根据ρ1 (r)和λ1 (r)的值,采用下式标定铣削力系数:
[kt,mT,kr,mR,kz,mz]T=(BTB)-1BTd
其中,
(f)使用步骤(e)标定得到铣削力系数预测铣削力,并计算与所有采样点对应的预测铣削力与实测铣削力之间偏差的平方和Δ(r);
(g)设置r=r+1,并给ρ1 (r)和λ1 (r)重新赋值,重复步骤(e)和步骤(f),直到Δ(r)达到最小值,当Δ(r)取最小值时,标定得到的铣削力系数以及ρ1 (r)和λ1 (r)的值为最终的kT,mT,kR,mR,kz,mz和ρ1,λ1结果。
2、一种刀具偏心四自由度参数的标定方法,其特征在于包括下述步骤:
(h)设定立铣刀的半径R、螺旋角β、刀齿数N、刀具安装到主轴后的悬臂长度L;设定标定试验的切削参数:单齿进给量f、轴向切削深度Rz、径向切削深度Rr、刀具主轴转速;Rz>R,Rr=0.5~2mm;工件是长方体块;
(j)设置r=0,并给定偏心参数的初值,即设φ(r)=φ0,τ(r)=τ0;
(k)基于φ(r)和τ(r)的值,使用下式计算对应的ρ(r)和λ(r):
其中,
A=2(L-Rz/2)sinτ(r)cosφ(r)
式中,ρ1表示刀具旋转中心与刀具几何中心间的偏移量,λ1表示刀具偏心产生的方向与相邻最近的刀齿头部之间的夹角;
(1)使用权利要求1标定得到铣削力系数kT,mT,kR,mR,kZ,mZ以及步骤(1)中的ρ(r)、λ(r)和φ(r)、τ(r)的值预测铣削力,并计算与所有采样点对应的预测铣削力与实测铣削力之间偏差的平方和δ(r);
(m)设置r=r+1,并给φ(r)、τ(r)重新赋值,重复步骤(1)和步骤(m),直到δ(r)达到最小值,当δ(r)取最小值时,对应的ρ(r)、λ(r)以及φ(r)、τ(r)值为最终的刀具偏心ρ、λ、φ和τ结果。
本发明的有益效果是:由于通过将测试瞬时铣削力转换到局部坐标系,并采用单纯形法实现了快速对瞬时铣削力系数和刀具偏心参数进行标定,同时考虑了铣削力系数的尺寸效应以及刀具偏心沿轴线的变化,无需进行参数研究,相比现有技术提高了计算效率;利用两组试验的铣削力测试结果互相配合进行刀具偏心参数标定,不需要增加新的试验手段与测试工艺,降低了标定过程对试验数目或试验仪器的依赖性;本发明通过实测铣削力进行刀具偏心参数标定,实现了对铣削过程动态偏心参数的标定,克服了试验手段无法标定动态偏心的不足;本发明通过两次铣削试验即可完成标定,无需进行大量的铣削试验,提高了参数标定效率,降低了试验成本。
将标定刀具径向偏心参数ρ1和λ1的工作量由参数研究所需的18000个迭代步降低到18个迭代步;将标定刀具轴向偏心参数φ和τ的工作量由参数研究所需的21600个迭代步降低到38个迭代步,提高了工作效率,降低了试验成本。
下面结合附图和实施例对本发明作详细说明。
附图说明
图1是背景技术中刀具径向偏心参数τ的示意图。
图2是背景技术中刀具偏心参数ρ、λ、φ的示意图。
图3是背景技术中刀具径向偏心参数ρ1和λ1的示意图。
图中,■-刀具几何中心,●-刀具旋转中心,ρ-刀具旋转中心与刀具几何中心间的偏移量,τ-机床主轴轴线与刀具轴线间的夹角,φ-刀具主轴轴线倾斜方向与偏移量,ρ-产生的方向间的夹角,λ-刀具偏心产生的方向与相邻最近的刀齿头部之间的夹角。
具体实施方式
(1)选定半径R为8mm、螺旋角β为30度的三齿硬质合金立铣刀在三坐标立铣床上对铝合金Al7050进行顺铣切削,刀具装夹悬臂长度L为52.76mm。第一组切削参数:刀具主轴转速为1200RPM,单齿进给量0.1667mm/齿,轴向切削深度Rz等于2mm,径向切削深度Rr等于8mm;第二组切削参数:刀具主轴转速为1000RPM,单齿进给量0.1mm/齿,轴向切削深度Rz等于25mm,径向切削深度Rr等于0.8mm。
(2)根据步骤(1)设定的切削参数,进行铣削并测试铣削力,铣削时要求工件被加工面与刀具轴线垂直。用表示对应于第i个刀齿切削周期内的第j个采样点的相角,将对应于的瞬时铣削力记为当s=1时表示第一组参数所测的铣削力,当s=2时表不第二组参数所测的铣削力。
(3)根据步骤(2)测试得到的铣削力,进行铣削力系数和刀具偏心参数的标定。首先利用第一组测试结果标定瞬时铣削力系数kT,mT,kR,mR,kZ,mZ和对应的径向偏心参数ρ1和λ1。然后基于第一组测试的标定得到的kT,mT,kR,mR,kZ,mZ和ρ1,λ1,使用第二组测试结果,进行实际偏心参数ρ、λ、φ和τ的标定。标定过程如下:
(3.1)使用第一组测试结果标定铣削力系数kT,mT,kR,mR,kZ,mZ和径向偏心参数ρ1,λ1。
(3.1.2)设置r=0,并给定径向偏心参数的初值,即假设 r表示迭代步。λ0和ρ0的值可利用实测铣削力的平均值和谐波分量通过下式得到
式中
FYO M[1]表示Y向铣削力的一阶谐波分量,FX M表示X向铣削力的平均值,FY M表示Y向铣削力的平均值。
(3.1.3)根据ρ1 (r)和λ1 (r)的值,采用下式标定铣削力系数:
[kt,mT,kr,mR,kz,mZ]T=(BTB)-1BTd
其中:
(3.1.4)基于第一组测试的切削参数和ρ1 (r)和λ1 (r)的值,使用步骤(3.1.3)标定得到铣削力系数预测铣削力,并计算与所有采样点对应的预测铣削力与实测铣削力之间偏差的平方和Δ(r)。
(3.1.5)设置r=r+1,并给ρ1 (r)和λ1 (r)重新赋值,重复步骤(3.1.3)和(3.1.4),直到Δ(r)达到最小值。当Δ(r)取最小值时,标定得到的铣削力系数以及ρ1 (r)和λ1 (r)的值为最终的kT,mT,kR,mR,kZ,mZ和ρ1,λ1结果。在此步骤中,使用单纯形法给ρ1 (r)和λ1 (r)赋值,并判断Δ(r)是否达到最小的。r表示迭代步。
通过上面的步骤,标定得到的结果如下表所示:
ρ<sub>1</sub> (mm) | λ<sub>1</sub> (度) | k<sub>T</sub> (N/mm<sup>2</sup>) | m<sub>T</sub> | k<sub>R</sub> (N/mm<sup>2</sup>) | m<sub>R</sub> | k<sub>Z</sub> (N/mm<sup>2</sup>) | m<sub>Z</sub> |
0.03284 | 85.69 | 758.17 | -0.1723 | 86.10 | -0.6609 | 143.10 | -0.1555 |
采用以上步骤经过18个迭代步就达到了收敛;倘若采用参数研究的方法进行标定,即如果让ρ1从0.001mm以0.001mm为步长变化到0.050mm,让λ1从1度以1度为步长变化到360度,然后从中选择使得预测铣削力与实测铣削力之间偏差最小的结果作为标定结果,需要经过360×50=18000步才能完成标定;提高了工作效率,降低了试验成本。
(3.2)基于第二组测试结果,进行实际偏心参数ρ、λ、φ和τ的标定。
(3.2.1)设置r=0,并给定径向偏心参数的初值,即假设φ(r)=φ0,τ(r)=τ0。φ0和τ0一般设定为接近于0的数。
(3.2.2)并基于φ(r)和τ(r)的值,使用下式计算对应的ρ(r)和λ(r):
其中,
A=2(L-Rz/2)sinτ(r)cosφ(r)
ρ1=0.03284
λ1=65.69°
(3.2.3)基于第二组测试的切削参数,使用由第一组测试标定得到铣削力系数kT,mT,kR,mR,kZ,mZ以及步骤(3.2.2)中的ρ(r)、λ(r)和φ(r)、τ(r)的值预测铣削力,并计算与所有采样点对应的预测铣削力与实测铣削力之间偏差的平方和δ(r)。
(3.2.4)设置r=r+1,并给φ(r)、τ(r)重新赋值,重复步骤(3.2.2)和(3.2.3),直到δ(r)达到最小值。当δ(r)取最小值时,对应的ρ(r)、λ(r)以及φ(r)、τ(r)值为最终的刀具偏心ρ、λ、φ和τ结果。在此步骤中,使用单纯形法给φ(r)、τ(r)赋值,并判断δ(r)是否达到最小的。
通过上面的步骤,标定得到的最终偏心参数结果如下表所示:
ρ(mm) | λ(度) | φ(度) | τ(度) |
0.01450 | 131.95 | 70.85 | 0.02779 |
采用以上步骤经过38个迭代步就达到了收敛;倘若采用参数研究的方法进行标定,即如果让φ从从1度以1度为步长变化到360度,让τ从0.001度以0.001度为步长变化到0.060度,然后从中选择使得预测铣削力与实测铣削力之间偏差最小的结果作为标定结果,需要经过360×60=21600步才能完成标定;提高了工作效率,降低了试验成本。
Claims (2)
1.一种刀具瞬时铣削力系数标定方法,其特征在于包括下述步骤;
(a)设定立铣刀的半径R、螺旋角β、刀齿数N、刀具安装到主轴后的悬臂长度L;设定标定试验的切削参数:单齿进给量f、轴向切削深度Rz、径向切削深度Rr、刀具主轴转速;Rz≤2mm,Rr=R;工件是长方体块;
(e)根据ρ1 (r)和λ1 (r)的值,采用下式标定铣削力系数:
[kt,mT,kr,mR,kz,mZ]T=(BTB)-1BTd
其中,
(f)使用步骤(e)标定得到铣削力系数预测铣削力,并计算与所有采样点对应的预测铣削力与实测铣削力之间偏差的平方和Δ(r);
(g)设置r=r+1,并给ρ1 (r)和λ1 (r)重新赋值,重复步骤(e)和步骤(f),直到Δ(r)达到最小值,当Δ(r)取最小值时,标定得到的铣削力系数以及ρ1 (r)和λ1 (r)的值为最终的kT,mT,kR,mR,kZ,mZ和ρ1,λ1结果。
2.一种刀具偏心四自由度参数的标定方法,其特征在于包括下述步骤:
(h)设定立铣刀的半径R、螺旋角β、刀齿数N、刀具安装到主轴后的悬臂长度L;设定标定试验的切削参数:单齿进给量f、轴向切削深度Rz、径向切削深度Rr、刀具主轴转速;Rz>R,Rr=0.5~2mm;工件是长方体块;
(j)设置r=0,并给定偏心参数的初值,即设φ(r)=φ0,τ(r)=τ0;
(k)基于φ(r)和τ(r)的值,使用下式计算对应的ρ(r)和λ(r):
其中,
A=2(L-Rz/2)sinτ(r)cosφ(r)
式中,ρ1表示刀具旋转中心与刀具几何中心间的偏移量,λ1表示刀具偏心产生的方向与相邻最近的刀齿头部之间的夹角;
(l)使用权利要求1标定得到铣削力系数kT,mT,kR,mR,kZ,mZ以及步骤(1)中的ρ(r)、λ(r)和φ(r)、τ(r)的值预测铣削力,并计算与所有采样点对应的预测铣削力与实测铣削力之间偏差的平方和δ(r);
(m)设置r=r+1,并给φ(r)、τ(r)重新赋值,重复步骤(1)和步骤(m),直到δ(r)达到最小值,当δ(r)取最小值时,对应的ρ(r)、λ(r)以及φ(r)、τ(r)值为最终的刀具偏心ρ、λ、φ和τ结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100211266A CN101549469B (zh) | 2009-02-16 | 2009-02-16 | 圆周铣削过程中刀具参数的标定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2009100211266A CN101549469B (zh) | 2009-02-16 | 2009-02-16 | 圆周铣削过程中刀具参数的标定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101549469A CN101549469A (zh) | 2009-10-07 |
CN101549469B true CN101549469B (zh) | 2010-09-29 |
Family
ID=41154083
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2009100211266A Expired - Fee Related CN101549469B (zh) | 2009-02-16 | 2009-02-16 | 圆周铣削过程中刀具参数的标定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101549469B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102350655B (zh) * | 2011-07-07 | 2013-06-05 | 西北工业大学 | 平头立铣刀铣削过程径向偏心参数标定方法 |
CN108107840B (zh) * | 2017-12-20 | 2019-12-03 | 上海交通大学 | 一种铣削力系数及偏心参数的标定方法 |
-
2009
- 2009-02-16 CN CN2009100211266A patent/CN101549469B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN101549469A (zh) | 2009-10-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101491844B (zh) | 圆周铣削过程铣削力系数与刀具径向偏心标定方法 | |
Bissacco et al. | Modelling the cutting edge radius size effect for force prediction in micro milling | |
CN104778333B (zh) | 一种螺旋立铣刀正交车铣加工三维稳定性建模方法 | |
CN102004818B (zh) | 自由曲线外形零件圆周铣削过程中瞬时铣削力建模方法 | |
Chaves-Jacob et al. | Optimal strategy for finishing impeller blades using 5-axis machining | |
CN101590614B (zh) | 基于形状复制的数控铣削加工刀具磨损测量方法 | |
CN101853324A (zh) | 圆周铣过程中铣削力建模方法 | |
CN102129232A (zh) | 五轴侧铣加工工艺参数设计方法 | |
Perez et al. | An enhanced method for cutting force estimation in peripheral milling | |
CN105701323A (zh) | 一种插铣加工切削力预测建模方法 | |
CN109754332B (zh) | 基于切削力的机床铣削加工过程的能耗模型建模方法 | |
CN105873703A (zh) | 切削加工中的切削条件的设计方法 | |
CN101549469B (zh) | 圆周铣削过程中刀具参数的标定方法 | |
Rao et al. | Effect of workpiece curvature on cutting forces and surface error in peripheral milling | |
Liu et al. | Optimal fixture design in peripheral milling of thin-walled workpiece | |
CN108647413B (zh) | 一种微细表面位置误差与稳定性综合预测方法 | |
CN102350655B (zh) | 平头立铣刀铣削过程径向偏心参数标定方法 | |
CN102152071B (zh) | 罗茨鼓风机叶轮组件的加工工艺 | |
CN106077775B (zh) | 一种三轴机床加工倒扣类工件的加工方法 | |
CN103639494B (zh) | 一种加工航空发动机导向叶片小锁板的方法 | |
Li et al. | Research on surface morphology of the ruled surface in five-axis flank milling | |
Checchi et al. | Offline tool trajectory compensation for cutting forces induced errors in a portable machine tool | |
CN109304642B (zh) | 一种刀具工艺系统的动态误差补偿方法 | |
CN108145221B (zh) | 管线泵腹腔加工方法 | |
CN209849909U (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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20100929 Termination date: 20130216 |