CN113553727A - 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 - Google Patents

一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 Download PDF

Info

Publication number
CN113553727A
CN113553727A CN202110922315.1A CN202110922315A CN113553727A CN 113553727 A CN113553727 A CN 113553727A CN 202110922315 A CN202110922315 A CN 202110922315A CN 113553727 A CN113553727 A CN 113553727A
Authority
CN
China
Prior art keywords
milling
force
amplitude
milling force
axis
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.)
Pending
Application number
CN202110922315.1A
Other languages
English (en)
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.)
Northeastern University China
Original Assignee
Northeastern University China
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 Northeastern University China filed Critical Northeastern University China
Priority to CN202110922315.1A priority Critical patent/CN113553727A/zh
Publication of CN113553727A publication Critical patent/CN113553727A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Numerical Control (AREA)

Abstract

本发明提供一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,首先通过设计槽铣加工实验获取一个周期内力的均值向量和幅值向量,然后设计槽铣加工实验中铣削力均值、幅值的理论计算模型,最后根据铣削力均值、幅值的理论计算模型以及实验测得的均值向量和幅值向量,建立包含待求解铣削力系数的优化目标,通过优化算法求解出最优解得到铣削力系数;本发明方法具有较高精度,且相比于其他利用平均力求解系数的方法,可以减少所需要的槽铣实验组数,在保证精度的同时,提高工作效率,降低试验成本。

Description

一种同时考虑铣削力幅值与均值的铣削力系数辨识方法
技术领域
本发明属于铣削力预测技术领域,具体涉及一种同时考虑铣削力幅值与均值的铣削力系数辨识方法。
背景技术
球头铣刀铣削作为一种重要的加工方法,广泛应用于模具制造、航空航天等工业。在铣削加工中,铣削力对刀具寿命、工件加工精度和机床振动等有显著影响,所以精准和高效的预测铣削力对铣削工艺的优化起着至关重要的作用。
铣削力的预测是一切关于铣削加工机理研究的基础,铣削力系数的精准辨识又是铣削力预测的基础。近年来,国内外学者对铣削力建模进行了大量研究。参照附图2-3,文献1“Lee,P.,and Y.Altintas.1996.Prediction of ball-end milling forces fromorthogonal cutting data.International Journal Machine Tools and Manufacturing1996:1059-1072.”提出了一种铣削力计算模型应用最为广泛,其将切削刃离散为微元,然后对微元铣削力进行积分得到最终总铣削力。微元铣削力模型如下:
Figure BDA0003207863960000011
dFt,dFr,dFa为微元切向力、微元径向力、微元轴向力;Ktc(k),Krc(k),Kac(k)为切向、径向和轴向剪切力切削力系数;Kte(k),Kre(k),Kae(k)为切向、径向和轴向耕犁力切削力系数;db为切削宽度,dS为切削刃长度,k为轴向位置角。
剪切力系数表示为切削刃微元轴向位置角的多项式函数,耕犁力系数表示为常数。
Ktc(k)=Ktc0+Ktc1k+Ktc2k2+…+Ktcmkm
Krc(k)=Krc0+Krc1k+Krc2k2+…+Krcmkm
Kac(k)=Kac0+Kac1k+Kac2k2+…+Kacmkm
其中Kic0,Kic1,Kic2,…,Kicm(i=t,r,a)为常数。
将微元铣削力转移到刀具坐标系中:
Figure BDA0003207863960000012
其中θ为微元周向位置角。将微元铣削力积分得到刀具坐标系中总铣削力:
Figure BDA0003207863960000021
文献2“Lamikiz,A.,L.N.Lopez de Lacalle,J.A.Sanchez,andU.Bravo.2005.Calculation of the specific cutting coefficients and geometricalaspects in sculptured surface machining.Machining Science and Technology 9:411-436.”公开了一种利用平均力标定铣削力系数的方法,该方法的基本步骤如下:
(1)设置槽切实验。
(2)取实验中铣削力数据的平均值。
(3)建立平均力理论值的计算模型。
(4)利用最小二乘法将实验铣削力平均值与理论平均值拟合,得出铣削力系数值。
通过上述步骤可以看出,该方法只考虑了实验数据的平均值,槽切组数n与待求切削力系数的个数p必须满足3n>p,需要的实验组数多。
发明内容
基于上述问题,为提高铣削力预测精度和效率,本发明提出一种同时考虑铣削力均值与幅值的铣削力系数辨识方法。首先进行槽铣加工实验,测试并记录实验过程中的瞬时铣削力;然后计算铣削力信号的均值与幅值,推导铣削力均值的理论模型,计算每组实验中Y、Z轴方向上铣削力最大值与最小值可能出现的周向位置角范围;最后以实验平均力向量与理论平均力向量差的最小二范数为目标函数,以实验幅值向量与理论幅值向量差的二范数小于一定值为约束函数,通过最优化算法得出铣削力系数。
本发明提出的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,包括:
步骤1:设计n组槽铣加工实验获取铣削力的瞬时值,并计算出各实验中一个周期内力的均值向量和幅值向量;
步骤2:设计槽铣加工实验中铣削力均值、幅值的理论计算模型;
步骤3:根据铣削力均值、幅值的理论计算模型以及实验测得的均值向量和幅值向量,建立包含待求解铣削力系数的优化目标,通过优化算法求解出最优解得到铣削力系数。
所述步骤1包括:
步骤1.1:进行槽铣加工前确定加工参数,包括刀具半径D,公称螺旋角β,刀具齿数N,每齿进给量fz,轴向切削深度dr;
步骤1.2:进行槽铣加工实验时,保持铣刀轴线与工件表面垂直;
步骤1.3:进行n次槽铣加工实验,每次实验时测量铣削力在X轴、Y轴、Z轴的瞬时值;
步骤1.4:根据测得的瞬时值计算出铣削力在一个周期内的均值,并将铣削力在X轴、Y轴、Z轴的均值组成均值向量,均值向量记为
Figure BDA0003207863960000031
其中,
Figure BDA0003207863960000032
分别表示第i次槽铣加工中一个周期内铣削力在X轴、Y轴、Z轴的均值,T表示矩阵的转置,i=1,2,…,n;
步骤1.5:根据测得的瞬时值计算出铣削力在一个周期内的幅值,即最大值和最小值之差,并将Y轴、Z轴方向的幅值构成幅值向量,幅值向量记为
Figure BDA0003207863960000033
其中,
Figure BDA0003207863960000034
分别表示第i次槽铣加工中一个周期内铣削力幅值在Y轴、Z轴的幅值。
所述步骤2中铣削力均值理论计算模型的构建,包括:
步骤2.1.1:根据公式(1)确定微元铣削力模型:
Figure BDA0003207863960000035
式中,dFt、dFr、dFa分别为微元切向力、微元径向力、微元轴向力;Ktc(k)、Krc(k)、Kac(k)分别为切向、径向和轴向的剪切力切削系数;Kte(k)、Kre(k)、Kae(k)分别为切向、径向和轴向的耕犁力切削系数;db为切削宽度,db=Rdk,k为切削刃微元轴向位置角,dS为切削刃长度,
Figure BDA0003207863960000036
θ为切削刃微元位置角,R为刀具半径,fn为瞬时切屑厚度;
步骤2.1.2:根据公式(2)计算切削时的瞬时切屑厚度fn(θ):
Figure BDA0003207863960000037
式中,fz为每齿进给量;
步骤2.1.3:刀具旋转一周相当于单个切削刃上的每一个微元都经过一次切削区域,结合公式(1)~(2)推导出铣削力均值向量的理论计算模型为:
Figure BDA0003207863960000041
其中,
Figure BDA0003207863960000042
式中,m为剪切力切削系数的最高次,N为刀具切削刃个数,kup为参与切削的切削刃轴向位置角上限,K为由待求解铣削力系数构成的向量;
Figure BDA0003207863960000043
则铣削力均值的理论计算模型简化为:
Figure BDA0003207863960000044
所述步骤2中铣削力幅值理论计算模型的构建,包括:
步骤2.2.1:确定Y轴和Z轴方向上铣削力最大值的理论切削刃位置角范围,其中最大值的位置角下限ψia定义为:
Figure BDA0003207863960000051
dr为槽切的轴向切削深度,β为公称螺旋角;
最大值的位置角ψib上限定义为:
Figure BDA0003207863960000052
步骤2.2.2:确定Y轴和Z轴方向上铣削力最小值的理论切削刃位置角范围,其中最小值的位置角下限ψic定义为:
Figure BDA0003207863960000053
最小值的位置角上限ψid定义为:
Figure BDA0003207863960000054
步骤2.2.3:从位置角[ψia,ψib]中等间距取m个角度值,并求每个角度值对应的铣削力
Figure BDA0003207863960000055
j=1,2,…,m,从位置角[ψic,ψid]中等间距取m个角度值求其对应的
Figure BDA0003207863960000056
计算
Figure BDA0003207863960000057
的最大值
Figure BDA0003207863960000058
Figure BDA0003207863960000059
的最小值
Figure BDA00032078639600000510
则第i组的理论幅值为
Figure BDA00032078639600000511
所述步骤3中所述包含待求解铣削力系数的优化目标表述为:
Figure BDA00032078639600000512
式中,Xi为第i组实验参数下的矩阵X,ai=[aiy,aiz]T为幅值的理论值构成的向量,aiy、aiz分别为第i组实验下Y、Z轴方向上幅值的理论值,t为铣削力幅值容许的误差值。
所述步骤1设计槽铣加工实验时,为了尽量减少所需要的槽铣加工实验组数,针对相同的轴向切深,至多取两组不同的进给速度;槽铣加工实验组数n与待求切削力系数的个数p满足如下条件:5n>p。
所述铣削力幅值容许的误差值t具体表示为:当机床运转但铣刀未接触工件时测得铣削力的幅值。
本发明的有益效果是:
本发明提出了一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,推导出了槽铣加工中铣削力均值、幅值的理论计算模型,建立了包含待求解铣削力系数的优化目标,通过优化算法求解最优解得到铣削力系数,本发明方法具有较高精度,且相比于其他利用平均力求解系数的方法,可以减少所需要的槽铣实验组数,在保证精度的同时,提高工作效率,降低试验成本。
附图说明
图1为本发明实施例中同时考虑铣削力幅值与均值的铣削力系数辨识方法流程图;
图2为本发明实施例中球头铣刀切削刃轨迹示意图,图中,刀具坐标系x轴方向为刀具进给方向,z方向为刀轴方向,y方向与x方向和z方向符合右手定则,ω为主轴旋转速度,k为轴向位置角,θ为微元周向位置角,
Figure BDA0003207863960000061
为微元螺旋滞后角,ψ为切削刃位置角;
图3为本发明实施例中微元瞬时切屑示意图,图中fn为瞬时切削厚度,db为微元宽度;
图4为本发明实施例中实验值与模拟值的对比图,其中,(a)为x轴方向实验值与模拟值的对比图;(b)为y轴方向实验值与模拟值的对比图;(c)为z轴方向实验值与模拟值的对比图。
具体实施方式
下面结合附图和具体实施实例对发明做进一步说明。为了克服现有技术在进行铣削过程铣削力系数辨识中的实验组数需求量大的问题,本发明提供一种铣削力系数辨识方法,首先设置参数,计算实测铣削力一个周期内的平均值与幅值;然后通过数学推导建立理论平均铣削力的计算模型,并求得y,z方向铣削力最大值与最小值的周向角范围,最后通过优化算法得到铣削力系数。同时考虑平均值与幅值的铣削力系数辨识方法,可以提高工作效率,降低试验成本。
本发明主要针对文献1(Lee,P.,and Y.Altintas.1996.Prediction of ball-endmilling forces from orthogonal cutting data.International Journal MachineTools and Manufacturing1996:1059-1072)中两齿球头铣刀的铣削力计算形式,待求的各铣削力系数构成的向量表示为K=[Ktc0,Krc0,Kac0,Ktc1,Krc1,Kac1,…,Ktcm,Krcm,Kacm,Kte,Kre,Kae]T,m为大于1的正整数。
为了达到上述目的,本发明提出一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,如图1所示,包括:
步骤1:设计n组槽铣加工实验获取铣削力的瞬时值,并计算出各实验中一个周期内力的均值向量和幅值向量;包括:
步骤1.1:槽铣加工前确定加工参数,包括刀具半径D,公称螺旋角β,刀具齿数N,每齿进给量fz,轴向切削深度dr;为了尽量减少所需要的槽铣加工实验组数,针对相同的轴向切深,至多取两组不同的进给速度;槽铣加工实验组数n与待求切削力系数的个数p满足如下条件:5n>p。
步骤1.2:进行槽铣加工实验时,保持铣刀轴线与工件表面垂直;
步骤1.3:进行n次槽铣加工实验,每次实验时测量铣削力在X轴、Y轴、Z轴的瞬时值;
步骤1.4:根据测得的瞬时值计算出铣削力在一个周期内的均值,并将铣削力在X轴、Y轴、Z轴的均值组成均值向量,均值向量记为
Figure BDA0003207863960000071
其中,
Figure BDA0003207863960000072
分别表示第i次槽铣加工实验中一个周期内铣削力在X轴、Y轴、Z轴的均值,T表示矩阵的转置,i=1,2,…,n;
步骤1.5:根据测得的瞬时值计算出铣削力在一个周期内的幅值,即最大值和最小值之差,并将Y轴、Z轴方向的幅值构成幅值向量,幅值向量记为
Figure BDA0003207863960000073
其中,
Figure BDA0003207863960000074
分别表示第i次槽铣加工实验中一个周期内铣削力幅值在Y轴、Z轴的幅值。
设计槽铣加工实验时,采用某立式升降台铣床对一块宽8cm,长16cm,厚3cm的45号钢板进行铣削加工,选用两切削刃硬质合金球头铣刀为切削刀具,刀具直径D=10mm,公称螺旋角β=30°;设计四组槽铣加工实验,测试并记录实验过程中的瞬时铣削力;计算铣削力实验值的均值与幅值,具体切削参数和实验结果表1所示:
表1切削参数和实验结果表
Figure BDA0003207863960000075
步骤2:设计槽铣加工实验中铣削力均值、幅值的理论计算模型;
其中,铣削力均值理论计算模型的构建,包括:
步骤2.1.1:根据公式(1)确定微元铣削力模型:
Figure BDA0003207863960000081
式中,dFt、dFr、dFa分别为微元切向力、微元径向力、微元轴向力;Ktc(k)、Krc(k)、Kac(k)分别为切向、径向和轴向的剪切力切削系数;Kte(k)、Kre(k)、Kae(k)分别为切向、径向和轴向的耕犁力切削系数;db为切削宽度,db=Rdk,k为切削刃微元轴向位置角,dS为切削刃长度,
Figure BDA0003207863960000082
θ为切削刃微元位置角,R为刀具半径,fn为瞬时切屑厚度;
步骤2.1.2:根据公式(2)计算切削时的瞬时切屑厚度fn(θ):
Figure BDA0003207863960000083
式中,fz为每齿进给量;
步骤2.1.3:刀具旋转一周相当于单个切削刃上的每一个微元都经过一次切削区域,结合公式(1)~(2)推导出铣削力均值向量的理论计算模型为:
Figure BDA0003207863960000084
其中,
Figure BDA0003207863960000085
式中,m为剪切力切削系数的最高次,m=3,N为刀具切削刃个数,kup为参与切削的切削刃轴向位置角上限,K为由待求解铣削力系数构成的向量,设K=[Ktc0,Krc0,Kac0,Ktc1,Krc1,Kac1,…,Ktcm,Krcm,Kacm,Kte,Kre,Kae]T
Figure BDA0003207863960000091
则铣削力均值的理论计算模型简化为:
Figure BDA0003207863960000092
β=30°,将四组槽铣加工实验得到的每齿进给量分别代入上式,计算得到四组X的具体值。
铣削力幅值理论计算模型的构建,包括:
步骤2.2.1:确定Y轴和Z轴方向上铣削力最大值的理论切削刃位置角范围,其中最大值的位置角下限ψia定义为:
Figure BDA0003207863960000093
dr为槽切的轴向切削深度,β为公称螺旋角;
最大值的位置角ψib上限定义为:
Figure BDA0003207863960000094
步骤2.2.2:确定Y轴和Z轴方向上铣削力最小值的理论切削刃位置角范围,其中最小值的位置角下限ψic定义为:
Figure BDA0003207863960000095
最小值的位置角上限ψid定义为:
Figure BDA0003207863960000096
步骤2.2.3:从位置角[ψia,ψib]中等间距取m个角度值,并求每个角度值对应的铣削力
Figure BDA0003207863960000097
j=1,2,…,m,从位置角[ψic,ψid]中等间距取m个角度值求其对应的
Figure BDA0003207863960000098
计算
Figure BDA0003207863960000099
的最大值
Figure BDA00032078639600000910
Figure BDA00032078639600000911
的最小值
Figure BDA00032078639600000912
则第i组的理论幅值为
Figure BDA00032078639600000913
β=30°,将四组槽铣加工实验得到的轴向切削深度分别代入上式得到四组实验中Y,Z方向最大值的理论周向位置角范围,计算结果如表2所示:
表2理论周向位置角取值范围表
Figure BDA0003207863960000101
步骤3:根据铣削力均值、幅值的理论计算模型以及实验测得的均值向量和幅值向量,建立包含待求解铣削力系数向量K的优化目标,通过优化算法求解出K的最优解得到铣削力系数,其中包含待求解铣削力系数向量K的优化目标表述为:
Figure BDA0003207863960000102
式中,Xi为第i组实验参数下的矩阵X,ai=[aiy,aiz]T为幅值的理论值构成的向量,aiy、aiz分别为第i组实验下Y、Z轴方向上幅值的理论值,t为铣削力幅值容许的误差值,是指当机床运转但铣刀未接触工件时测得铣削力的幅值,通常取值为机床未启动时传感器的跳动幅值。
当t=7时,利用matlab中的fmincon优化函数求解得到的各参数值如下:
Figure BDA0003207863960000103
为了验证切削力模型的准确性,取切削深度dr=0.3,主轴转速r=600r/min,进给速度v=30mm/min,再对45号钢进行槽切加工,并将加工参数和上述铣削力系数代入本发明中的铣削力预测模型,得到模拟切削力。将模型计算的切削力值与实验获得的切削力信号进行对比,如图4所示;通过模拟值与实验值对比图可以看出,本方法利用四组铣削力实验得出来了高精度的铣削力系数。如果通过像文献2中的传统平均力法辨识铣削力系数,则至少需要五组实验。以上结果说明本发明方法具有高精度且节省实验成本。

Claims (7)

1.一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,包括:
步骤1:设计n组槽铣加工实验获取铣削力的瞬时值,并计算出各实验中一个周期内力的均值向量和幅值向量;
步骤2:设计槽铣加工实验中铣削力均值、幅值的理论计算模型;
步骤3:根据铣削力均值、幅值的理论计算模型以及实验测得的均值向量和幅值向量,建立包含待求解铣削力系数的优化目标,通过优化算法求解出最优解得到铣削力系数。
2.根据权利要求1所述的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,所述步骤1包括:
步骤1.1:槽铣加工前确定加工参数,包括刀具半径D,公称螺旋角β,刀具齿数N,每齿进给量fz,轴向切削深度dr;
步骤1.2:进行槽铣加工实验时,保持铣刀轴线与工件表面垂直;
步骤1.3:进行n次槽铣加工实验,每次实验时测量铣削力在X轴、Y轴、Z轴的瞬时值;
步骤1.4:根据测得的瞬时值计算出铣削力在一个周期内的均值,并将铣削力在X轴、Y轴、Z轴的均值组成均值向量,均值向量记为
Figure FDA0003207863950000011
其中,
Figure FDA0003207863950000012
分别表示第i次槽铣加工中一个周期内铣削力在X轴、Y轴、Z轴的均值,T表示矩阵的转置,i=1,2,…,n;
步骤1.5:根据测得的瞬时值计算出铣削力在一个周期内的幅值,并将Y轴、Z轴方向的幅值构成幅值向量,幅值向量记为
Figure FDA0003207863950000013
其中,
Figure FDA0003207863950000014
分别表示第i次槽铣加工中一个周期内铣削力幅值在Y轴、Z轴的幅值。
3.根据权利要求1所述的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,所述步骤2中铣削力均值理论计算模型的构建,包括:
步骤2.1.1:根据公式(1)确定微元铣削力模型:
Figure FDA0003207863950000015
式中,dFt、dFr、dFa分别为微元切向力、微元径向力、微元轴向力;Ktc(k)、Krc(k)、Kac(k)分别为切向、径向和轴向的剪切力切削系数;Kte(k)、Kre(k)、Kae(k)分别为切向、径向和轴向的耕犁力切削系数;db为切削宽度,db=Rdk,k为切削刃微元轴向位置角,dS为切削刃长度,
Figure FDA0003207863950000016
θ为切削刃微元位置角,R为刀具半径,fn为瞬时切屑厚度;
步骤2.1.2:根据公式(2)计算切削时的瞬时切屑厚度fn(θ):
Figure FDA0003207863950000021
式中,fz为每齿进给量;
步骤2.1.3:刀具旋转一周相当于单个切削刃上的每一个微元都经过一次切削区域,结合公式(1)~(2)推导出铣削力均值向量的理论计算模型为:
Figure FDA0003207863950000022
其中,
Figure FDA0003207863950000023
式中,m为剪切力切削系数的最高次,N为刀具切削刃个数,kup为参与切削的切削刃轴向位置角上限,K为由待求解铣削力系数构成的向量;
Figure FDA0003207863950000024
则铣削力均值的理论计算模型简化为:
Figure FDA0003207863950000031
4.根据权利要求1所述的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,所述步骤2中铣削力幅值理论计算模型的构建,包括:
步骤2.2.1:确定Y轴和Z轴方向上铣削力最大值的理论切削刃位置角范围,其中最大值的位置角下限ψia定义为:
Figure FDA0003207863950000032
dr为槽切的轴向切削深度,β为公称螺旋角;
最大值的位置角ψib上限定义为:
Figure FDA0003207863950000033
步骤2.2.2:确定Y轴和Z轴方向上铣削力最小值的理论切削刃位置角范围,其中最小值的位置角下限ψic定义为:
Figure FDA0003207863950000034
最小值的位置角上限ψid定义为:
Figure FDA0003207863950000035
步骤2.2.3:从位置角[ψia,ψib]中等间距取m个角度值,并求每个角度值对应的铣削力
Figure FDA0003207863950000036
从位置角[ψic,ψid]中等间距取m个角度值求其对应的
Figure FDA0003207863950000037
计算
Figure FDA0003207863950000038
的最大值
Figure FDA0003207863950000039
Figure FDA00032078639500000310
的最小值
Figure FDA00032078639500000311
则第i组的理论幅值为
Figure FDA00032078639500000312
5.根据权利要求1所述的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,所述步骤3中所述包含待求解铣削力系数的优化目标表述为:
Figure FDA00032078639500000313
式中,Xi为第i组实验参数下的矩阵X,ai=[aiy,aiz]T为幅值的理论值构成的向量,aiy、aiz分别为第i组实验下Y、Z轴方向上幅值的理论值,t为铣削力幅值容许的误差值。
6.根据权利要求1所述的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,所述步骤1设计槽铣加工实验时,为了尽量减少所需要的槽铣加工实验组数,针对相同的轴向切深,至多取两组不同的进给速度;槽铣加工实验组数n与待求切削力系数的个数p满足如下条件:5n>p。
7.根据权利要求5所述的一种同时考虑铣削力幅值与均值的铣削力系数辨识方法,其特征在于,所述铣削力幅值容许的误差值t具体表示为:当机床运转但铣刀未接触工件时测得铣削力的幅值。
CN202110922315.1A 2021-08-12 2021-08-12 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法 Pending CN113553727A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110922315.1A CN113553727A (zh) 2021-08-12 2021-08-12 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110922315.1A CN113553727A (zh) 2021-08-12 2021-08-12 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法

Publications (1)

Publication Number Publication Date
CN113553727A true CN113553727A (zh) 2021-10-26

Family

ID=78105544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110922315.1A Pending CN113553727A (zh) 2021-08-12 2021-08-12 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法

Country Status (1)

Country Link
CN (1) CN113553727A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114633152A (zh) * 2022-05-05 2022-06-17 陕西科技大学 一种机床铣削加工工艺下的虚拟刀装置及加载方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101491844A (zh) * 2009-01-15 2009-07-29 西北工业大学 圆周铣削过程铣削力系数与刀具径向偏心标定方法
CN107480352A (zh) * 2017-07-28 2017-12-15 东北大学 一种铣削加工工艺参数的可靠性优化方法
CN107944149A (zh) * 2017-11-28 2018-04-20 清华大学 一种普通立铣刀圆弧段切削力系数辨识方法
CN107944176A (zh) * 2017-12-06 2018-04-20 上海复合材料科技有限公司 一种球头铣刀多轴铣削钛合金铣削力预测方法
CN109396955A (zh) * 2017-08-16 2019-03-01 山东大学 一种面向旋风包络铣削工艺的切削力预测方法及系统
CN110281087A (zh) * 2019-06-11 2019-09-27 北京航空航天大学 一种旋转超声振动磨削切削力预测方法
US20200164447A1 (en) * 2018-11-22 2020-05-28 Qingdao university of technology Milling system and method under different lubrication conditions
CN112949148A (zh) * 2021-04-21 2021-06-11 江苏电子信息职业学院 一种基于有限元模拟的刀具偏心跳动辨识方法
CN113158371A (zh) * 2021-04-21 2021-07-23 江苏电子信息职业学院 一种高速铣削加工动态切削力预测系统及参数优化方法
CN113177280A (zh) * 2021-05-20 2021-07-27 重庆大学 球刀切削力系数标定方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101491844A (zh) * 2009-01-15 2009-07-29 西北工业大学 圆周铣削过程铣削力系数与刀具径向偏心标定方法
CN107480352A (zh) * 2017-07-28 2017-12-15 东北大学 一种铣削加工工艺参数的可靠性优化方法
CN109396955A (zh) * 2017-08-16 2019-03-01 山东大学 一种面向旋风包络铣削工艺的切削力预测方法及系统
CN107944149A (zh) * 2017-11-28 2018-04-20 清华大学 一种普通立铣刀圆弧段切削力系数辨识方法
CN107944176A (zh) * 2017-12-06 2018-04-20 上海复合材料科技有限公司 一种球头铣刀多轴铣削钛合金铣削力预测方法
US20200164447A1 (en) * 2018-11-22 2020-05-28 Qingdao university of technology Milling system and method under different lubrication conditions
CN110281087A (zh) * 2019-06-11 2019-09-27 北京航空航天大学 一种旋转超声振动磨削切削力预测方法
CN112949148A (zh) * 2021-04-21 2021-06-11 江苏电子信息职业学院 一种基于有限元模拟的刀具偏心跳动辨识方法
CN113158371A (zh) * 2021-04-21 2021-07-23 江苏电子信息职业学院 一种高速铣削加工动态切削力预测系统及参数优化方法
CN113177280A (zh) * 2021-05-20 2021-07-27 重庆大学 球刀切削力系数标定方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
CHANGLI WANG ET AL.: "A method for predicting ball-end cutter milling force and its probabilistic characteristics", 《MECHANICS BASED DESIGN OF STRUCTURES AND MACHINES》, vol. 51, no. 6, 23 May 2021 (2021-05-23), pages 3416 - 3433 *
XIAOJUN LIN ET AL.: "The identification of the cutting force coefficients for ball-end finish milling", 《THE INTERNATIONAL JOURNAL OF ADVANCED MANUFACTURING TECHNOLOGY》, vol. 102, 7 April 2019 (2019-04-07), pages 4121 - 4135, XP036790959, DOI: 10.1007/s00170-019-03481-z *
姚琦 等: "考虑动力学影响的异构工艺系统切削力系数辨识方法", 《中国优秀硕士学位论文全文数据库 工程科技Ⅰ辑》, no. 7, 15 July 2021 (2021-07-15), pages 022 - 469 *
宋相强 等: "铣削加工表面位置误差的可靠性灵敏度分析", 《机械设计与制造》, no. 9, 8 September 2019 (2019-09-08), pages 86 - 88 *
黎柏春 等: "基于瞬时铣削力的球头铣刀铣削力系数辨识", 《东北大学学报(自然科学版)》, vol. 37, no. 5, 15 May 2016 (2016-05-15), pages 678 - 682 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114633152A (zh) * 2022-05-05 2022-06-17 陕西科技大学 一种机床铣削加工工艺下的虚拟刀装置及加载方法

Similar Documents

Publication Publication Date Title
CN106808320B (zh) 一种考虑刀具齿间角、螺旋角及偏心的铣削力预测方法
Zhenyu et al. Influence of dynamic effects on surface roughness for face milling process
Daniyan et al. Design and optimization of machining parameters for effective AISI P20 removal rate during milling operation
CN110032794B (zh) 一种振动作用下的铣刀动态切削力模型构建与验证方法
Uriarte et al. Mechanistic modelling of the micro end milling operation
Subrahmanyam et al. Cutting force prediction for ball nose milling of inclined surface
Daniyan et al. Modelling and optimization of the cutting forces during Ti6Al4V milling process using the response surface methodology and dynamometer
CN110262397A (zh) 车铣加工空间螺旋次摆线运动轨迹及瞬时切削力预测模型
CN107451382B (zh) 高速切削加工工件表面形貌的控制方法
WO2023138327A1 (zh) 一种旋风铣削中螺纹工件表面形貌的预测方法
CN113553727A (zh) 一种同时考虑铣削力幅值与均值的铣削力系数辨识方法
CN114186175A (zh) 振动作用下高能效铣刀主切削力能耗动态特性的解算方法
Akyildiz Evaluating of cutting forces in thread machining
CN102744450A (zh) 一种检测机床切削能力的测试试件及其应用
Gao et al. Effects of machining inclination angles on microgroove quality in micro ball end milling of Ti-6Al-4V
Fomin Microgeometry of surfaces after profile milling with the use of automatic cutting control system
Chen et al. Multiple regression prediction model for cutting forces and surface roughness in micro-milling of TA2
Li et al. A numerical study of the effects of cutter runout on milling process geometry based on true tooth trajectory
Matras et al. Precision milling of hardened steel with CBN tools
Bolar et al. Experimental investigation on surface quality and dimensional accuracy during curvilinear thin-wall machining
Lin et al. The identification of the cutting force coefficients for ball-end finish milling
CN202726640U (zh) 一种用于检测机床切削能力的测试试件
CN102059597B (zh) 一种硬齿面插齿刀负倒棱的加工方法
Liu et al. Improved dynamic cutting force model in ball-end milling. Part I: theoretical modelling and experimental calibration
CN113626953B (zh) 高能效铣削加工误差动态分布特性识别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination