一种椭圆超声振动辅助铣削切削力预测方法与流程

1.本发明涉及切削力建模领域，更具体地，涉及一种椭圆超声振动辅助铣削切削力预测方法。


背景技术：

2.切削过程中振动会导致切削行为不断发生改变，引起切削层参数的不断变化，导致切削力频繁变动，从而引起铣削过程中的切削载荷发生变化，最终引起刀具磨损和加工表面质量降低。同时，振动与切削力产生交互作用，并随着切削的进行，这种交互作用不断加强，最终导致切削失稳。
3.公开号为cn109408866a的中国发明专利公开了一种振动辅助铣削切削力建模方法，包括如下步骤：步骤一：通过振动辅助铣削运动学方程获得刀尖轨迹；步骤二：根据刀尖轨迹判断刀具和工件是否发生分离，并获得瞬时切削厚度；步骤三：根据步骤二获取的瞬时切削厚度，基于切削厚度的切削力计算模型得出获得振动辅助条件下的铣削力模型。上述技术方案在判断刀具和工件的分离情况时，需要计算刀具切削刃的位置，与刀具在t时刻之前所产生的轨迹的交点，过程复杂。
4.椭圆超声振动辅助铣削在切削刀具的刀柄上增加椭圆超声振动装置，使刀尖相对于工件作椭圆轨迹超声振动，从而实现超声频断续振动铣削。与常规铣削相比，椭圆超声振动辅助铣削使刀尖与工件频繁分离，可降低切削力，有利于切削液顺利进入切削区内部，降低切削区温度，减少刀具磨损，进而延长刀具使用寿命。目前，椭圆超声振动辅助铣削被广泛应用于航空难加工材料的切削加工。
5.与常规铣削类似，切削力也是影响椭圆超声振动铣削加工质量与效率的关键因素之一。因此进行准确的切削力预测，并在此基础上选择合适的超声振动铣削参数，对于实现高质量超声振动铣削至关重要。与常规铣削相比，椭圆超声振动辅助铣削刀尖运动轨迹更为复杂，切削形成机理也更为复杂，目前未见有关椭圆超声振动辅助铣削切削解析预测方法的报道。


技术实现要素：

6.本发明的目的在于针对现有技术中的不足，提供一种椭圆超声振动辅助铣削切削力预测方法，运用z-map模型来表示工件，进而根据刀尖运动轨迹方程计算瞬时切削厚度及动态切削力，简化了计算过程。
7.本发明的目的通过以下技术方案实现：
8.一种椭圆超声振动辅助铣削切削力预测方法，包括以下步骤：
9.s1.采用z-map模型对工件进行描述；
10.s2.建立刀具中心、刀尖运动轨迹方程以及速度方程；
11.s3.更新工件z-map模型并计算瞬时切削厚度；
12.s4.运用瞬时刚性力模型，计算椭圆超声振动辅助铣削在x、y和z三个方向动态切
削力。
13.进一步地，步骤s1中具体步骤如下：在工件的x和y方向划分坐标网格，网格的z坐标值用于表示工件的加工深度。
14.进一步地，步骤s2包括以下步骤：
15.s21.建立椭圆超声振动辅助铣削刀尖运动轨迹方程；
16.s22.建立椭圆超声振动辅助铣削刀尖运动速度方程。
17.进一步地，步骤s21中所述椭圆超声振动辅助铣削刀尖运动轨迹的建模过程如下：
18.椭圆超声振动辅助铣削的运动由铣刀绕主轴的回转运动、绕铣刀轴心椭圆振动及铣刀的直线进给运动复合而成，其中铣刀的进给运动方程为：
[0019][0020]
其中，vf为进给速度，f
t
为铣刀每齿进给量，n为主轴转速，n为刀齿数，t为时间；
[0021]
铣刀绕主轴回转运动方程为：
[0022][0023]
其中，r为铣刀半径，ω为铣刀回转角频率，刀尖在oxy平面内作超声椭圆振动，
[0024]
椭圆振动的轨迹方程为：
[0025][0026]
其中，a、b分别为超声椭圆振动在椭圆长轴和短轴两个方向上的振幅，f为超声振动频率；当铣刀旋转ωt时，椭圆长轴方向与ouv坐标系中的u轴方向一致；
[0027]
将ouv坐标系中椭圆方程转换为oxy坐标系中的轨迹方程，其坐标变换可表示为：
[0028][0029]
刀尖点运动轨迹方程为：
[0030][0031]
进一步地，
[0032]
步骤s22中所述椭圆超声振动辅助铣削刀尖运动速度方程的建模过程如下：
[0033]
对刀尖运动轨迹方程(5)求导可得刀尖运动速度方程为：
[0034][0035]
合成速度及其方向为：
[0036][0037]
切向和径向的切削速度分量分别为：
[0038][0039]
进一步地，步骤s3包括以下步骤：
[0040]
s31.利用公式(4)和公式(5)计算当前时刻刀具中心及刀尖点位置；
[0041]
s32.利用s31得到的先前时刻和当前时刻刀具中心及刀尖点坐标，构造出一个四边形用于近似表示该时间段内铣刀所切削的区域，并以dx和dy间距对该区域进行网格化处理；
[0042]
s33.在当前刀具中心点与刀尖点确定的线段上，利用二分法，快速确定该直线与工件的交点坐标，该点到刀尖点的距离即为瞬时切削厚度；
[0043]
s34.利用s32中四边形x、y方向上的最小最大坐标值，构造一个矩形。对矩形内的每个网格点进行判断，若该网格点位于四边形内部且该点的z坐标值未被更新，则更新该点z坐标值，利用matlab软件提供的库函数inploygon来确定某点是否在一个指定的四边形内部。
[0044]
进一步地，步骤s3还包括以下步骤：s35.重复s31～s34直至仿真结束。
[0045]
进一步地，步骤s4中具体过程如下：沿轴向将铣刀切削刃分割成若干微元，并对每个微元施用瞬时刚性力模型，进而获得瞬时切削力。
[0046]
进一步地，步骤s4包括以下步骤：
[0047]
s41.沿轴向将铣刀分隔成n个切削刃微元，切削刃上z高度处的径向滞后角为：
[0048]
ψ＝ztanβ/r(9)
[0049]
对于刀齿均布的铣刀，其齿间角为φ
p
＝2π/n，假设第一个刀齿刀尖处的角位移为φ
10
，那么刀齿j上的切削刃微元l处的角位移可表示为：
[0050]
φ
jl
＝φ
10
(j-1)φ
p
l dztanβ/r
ꢀꢀꢀ
(10)
[0051]
s42.施用瞬时刚性力模型，作用在刀齿j、微元l上的瞬时切削力可表示为：
[0052][0053]
式中，h为切削刃微元的瞬时切削厚度，可由s33.获得；k
tc
、k
rc
、k
ac
分别为切向、径向和轴向切削力系数；k
te
、k
re
、k
ae
分别为切向、径向和轴向刃口力系数。
[0054]
s43.通过坐标变换，可得到作用在机床坐标系中的切削力分量如下：
[0055][0056]
沿轴向对切削刃微元上的切削力积分，对每个刀齿的切削力求和，可得到三向瞬时切削力为：
[0057]
[0058]
进一步地，所述预测方法还包括以下步骤：s5.重复s3～s4，直至仿真结束。
[0059]
与现有技术相比，本发明的有益效果如下：
[0060]
本发明采用z-map模型来表示工件的切削过程，有效避免了切削刃与复杂空间曲线之间的求交点问题，大大简化了瞬时切削厚度以及瞬时切削力的计算过程；采用网格化与二分法相结合方法，快速有效实现了瞬时切削厚度的计算。
[0061]
本发明对切削力具有较高的预测精度，仿真结果也反映了相对于常规铣削，椭圆超声振动辅助铣削可显著提高切削速度，降低切削力。
附图说明
[0062]
图1为椭圆超声振动辅助铣削工件z-map模型表示；
[0063]
图2为椭圆超声振动辅助铣削刀尖运动合成示意图；
[0064]
图3椭圆振动轨迹随铣刀旋转变化示意图
[0065]
图4为椭圆超声振动辅助铣削刀尖运动轨迹图；
[0066]
图5为不同主轴转速下椭圆超声振动辅助铣削与常规铣削切削速度对比图；
[0067]
图6为椭圆超声振动辅助铣削瞬时切削厚度计算流程图；
[0068]
图7为椭圆超声振动辅助铣削与常规铣削瞬时切削厚度对比图；
[0069]
图8为椭圆超声振动辅助铣削切削力验证试验装置示意图；
[0070]
图9为椭圆超声振动辅助铣削切削力仿真值与实测值对比图；
[0071]
图10为椭圆超声振动辅助铣削与常规铣削瞬时切削力对比图。
具体实施方式
[0072]
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0073]
实施例1
[0074]
本实施例提供一种椭圆超声振动辅助铣削切削力预测方法，具体包括以下步骤：
[0075]
s1.采用z-map模型对工件进行描述；
[0076]
s2.建立刀具中心、刀尖运动轨迹方程以及速度方程；
[0077]
s3.更新工件z-map模型并计算瞬时切削厚度；
[0078]
s4.运用瞬时刚性力模型，计算椭圆超声振动辅助铣削在x、y和z三个方向动态切削力。
[0079]
其中，图1所示为椭圆超声振动辅助铣削工件z-map模型，工件在x和y方向被划分为一系列网格，网格的z坐标值用于表示工件的加工深度，通过计算椭圆超声振动辅助铣削加工过程中当前时刻和dt时刻后铣刀中心及刀尖点位置坐标，并用一个封闭四边形来表示该时间段内刀具切除的面积，进而通过计算封闭区域内网格数在刀具半径方向的投影得到切削厚度值。
[0080]
本实施例中步骤s2包括以下步骤：
[0081]
s21.建立椭圆超声振动辅助铣削刀尖运动轨迹方程；
[0082]
s22.建立椭圆超声振动辅助铣削刀尖运动速度方程。
[0083]
其中，s21为椭圆超声振动辅助铣削刀尖运动轨迹建模过程，具体如下：
[0084]
如图2所示，椭圆超声振动辅助铣削的运动由铣刀绕主轴的回转运动、绕铣刀轴心椭圆振动及铣刀的直线进给运动复合而成。其中铣刀的进给运动方程为：
[0085][0086]
式中：vf为进给速度，f
t
为铣刀每齿进给量，n为主轴转速，n为刀齿数，t为时间。铣刀绕主轴回转运动方程为：
[0087][0088]
式中：r为铣刀半径，ω为铣刀回转角频率。如图3所示，刀尖在oxy平面内作超声椭圆振动，椭圆振动的轨迹方程可表示为：
[0089][0090]
式中：a、b分别为超声椭圆振动在椭圆长轴和短轴两个方向上的振幅；f为超声振动频率。
[0091]
当铣刀旋转ωt时，椭圆长轴方向与ouv坐标系中的u轴方向一致。将ouv坐标系中椭圆方程转换为oxy坐标系中的轨迹方程，其坐标变换可表示为：
[0092][0093]
刀尖点运动轨迹方程为：
[0094][0095]
图4所示为仿真得到的椭圆超声振动辅助铣削刀尖轨迹示意图。
[0096]
s22为刀尖运动速度建模过程，具体如下：
[0097]
对刀尖运动轨迹方程求导可得刀尖运动速度方程为：
[0098][0099]
合成速度及其方向为：
[0100][0101]
切向和径向的切削速度分量分别为：
[0102][0103]
如图5所示，相同主轴转速下，椭圆超声振动辅助铣削的切削速度远高于常规铣
削，切削试验过程中的工艺参数如下：刀尖做椭圆振动时的长袖a＝0.012mm、刀尖做椭圆振动时的短轴b＝0.008mm、振动频率f＝20000hz、刀具半径r＝4mm、每齿的进给量f
t
＝0.1mm/tooth、刀齿数n＝2。
[0104]
在步骤s3中，可以根据刀尖运动位置，实现工件z-map模型的更新并计算瞬时切削厚度，具体步骤如下：
[0105]
s31.利用公式(4)和公式(5)计算当前时刻刀具中心及刀尖点为止；
[0106]
s32.利用s31得到的先前时刻和当前时刻刀具中心及刀尖点坐标，构造出一个四边形用于近似表示该时间段内铣刀所切削的区域，并以dx和dy间距对该区域进行网格化处理；
[0107]
s33.在当前刀具中心点与刀尖点确定的线段上，利用二分法，快速确定该直线与工件的交点坐标，该点到刀尖点的距离即为瞬时切削厚度；
[0108]
s34.利用s32中四边形x、y方向上的最小最大坐标值，构造一个矩形。对矩形内的每个网格点进行判断，若该网格点位于四边形内部且该点的z坐标值未被更新，则更新该点z坐标值；利用matlab软件提供的库函数inploygon来确定某点是否在一个指定的四边形内部。
[0109]
椭圆超声振动辅助铣削瞬时切削力仿真流程图如图6所示，不同主轴转速下椭圆超声振动辅助铣削与常规铣削切削厚度对比如图7所示，图7切削试验过程中的工艺参数如下：刀尖做椭圆振动时的长袖a＝0.012mm、刀尖做椭圆振动时的短轴b＝0.008mm、振动频率f＝20000hz、刀具半径r＝4mm、每齿的进给量f
t
＝0.2mm、刀齿数n＝2。
[0110]
本实施例步骤s4中，未考虑刀具螺旋角对切削力的影响，沿轴向将铣刀切削刃分割成若干微元，并对每个微元施用瞬时刚性力模型，进而获得瞬时切削力，其具体步骤如下：
[0111]
s41.沿轴向将铣刀分隔成n个切削刃微元，切削刃上z高度处的径向滞后角为：
[0112]
ψ＝ztanβ/r(9)
[0113]
对于刀齿均布的铣刀，其齿间角为φ
p
＝2π/n，假设第一个刀齿刀尖处的角位移为φ
10
，那么刀齿j上的切削刃微元l处的角位移可表示为：
[0114]
φ
jl
＝φ
10
(j-1)φ
p
l dztanβ/r
ꢀꢀꢀ
(10)
[0115]
s42.施用瞬时刚性力模型，作用在刀齿j、微元l上的瞬时切削力可表示为：
[0116][0117]
式中，h为切削刃微元的瞬时切削厚度，可由s33.获得；k
tc
、k
rc
、k
ac
分别为切向、径向和轴向切削力系数；k
te
、k
re
、k
ae
分别为切向、径向和轴向刃口力系数。
[0118]
s43.通过坐标变换，可得到作用在机床坐标系中的切削力分量如下：
[0119][0120]
沿轴向对切削刃微元上的切削力积分，对每个刀齿的切削力求和，可得到三向瞬
时切削力为：
[0121][0122]
椭圆超声振动辅助铣削切削力验证实验装置如图8所示。椭圆超声振动辅助铣削切削力仿真结果与实测值对比如图9所示，两者吻合良好，说明本发明提出的切削力预测方法有效；图9中切削试验过程中的工艺参数如下：转速n＝600r/min、每齿的进给量f
t
＝0.20mm、径向切深ae＝0.2mm、轴向切深a
p
＝5.0mm；刀尖做椭圆振动时的长袖a＝0.012mm、刀尖做椭圆振动时的短轴b＝0.008mm，振动频率f＝17880hz。椭圆超声振动辅助铣削与常规铣削切削力对比如图10所示，该图显示，椭圆超声振动辅助铣削可显著降低切削力；图9中切削试验过程中的工艺参数如下：每齿的进给量f
t
ft＝0.10mm、径向切深ae＝12.0mm、轴向切深a
p
＝2.0mm；刀尖做椭圆振动时的长袖a＝0.012mm、刀尖做椭圆振动时的短轴b＝0.008mm，振动频率f＝17880hz。
[0123]
本实施例的有益效果如下：
[0124]
本技术基于工件的z-map模型表示来实现椭圆超声振动辅助铣削过程中瞬时切削厚度的预测，并在此基础上通过借鉴常规铣削的切削力预测方案实现椭圆超声振动辅助铣削动态切削力的预测，避免了复杂的解析计算，实验验证结果证实，该方法具有较高的预测精度。仿真结果也客观反映了相对于常规铣削，椭圆超声振动辅助铣削可显著提高切削速度，降低切削力。
[0125]
实施例2
[0126]
本实施例参照实施例1，提供一种椭圆超声振动辅助铣削切削力预测方法，与实施例1不同之处在于，在步骤s3中还包括以下步骤，s35.重复s31～s34直至仿真结束。
[0127]
实施例3
[0128]
本实施例参照实施例1，提供一种椭圆超声振动辅助铣削切削力预测方法，与实施例1不同之处在于，本实施例中的预测方法还包括以下步骤，s5.重复s3～s4,直至仿真结束。
[0129]
需要说明，若本发明实施例中有涉及方向性指示(诸如上、下、左、右、前、后
……
)，则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之间的相对位置关系、运动情况等，如果该特定姿态发生改变时，则该方向性指示也相应地随之改变。
[0130]
另外，若本发明实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。