CN112380646B - 不同磨损刀具车削温度变化与车削振动耦合特性研究方法 - Google Patents

不同磨损刀具车削温度变化与车削振动耦合特性研究方法 Download PDF

Info

Publication number
CN112380646B
CN112380646B CN202011238945.9A CN202011238945A CN112380646B CN 112380646 B CN112380646 B CN 112380646B CN 202011238945 A CN202011238945 A CN 202011238945A CN 112380646 B CN112380646 B CN 112380646B
Authority
CN
China
Prior art keywords
turning
vibration
temperature change
parameters
temperature
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
Application number
CN202011238945.9A
Other languages
English (en)
Other versions
CN112380646A (zh
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.)
Shandong Yijite Cnc Machine Tool Co ltd
Original Assignee
Jiangsu Normal University
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 Jiangsu Normal University filed Critical Jiangsu Normal University
Priority to CN202011238945.9A priority Critical patent/CN112380646B/zh
Publication of CN112380646A publication Critical patent/CN112380646A/zh
Application granted granted Critical
Publication of CN112380646B publication Critical patent/CN112380646B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Abstract

不同磨损刀具车削温度变化与车削振动耦合特性研究方法,包括:制定车削试验方案;搭建试验系统;对不同磨损状态的刀具进行车削试验,得到相应车削参数下刀尖附近的车削温度和三向车削振动数据;从车削温度和三向车削振动数据中提取出温度特征值车削温度变化均值和振动特征值三向加速度均方根。根据温度特征值、振动特征值和车削参数,通过灰色系统理论分别计算三种磨损状态刀具的车削温度变化与车削参数的关联度、车削振动与车削参数的关联度、车削振动与车削温度变化的关联度。基于最小二乘法建立车削温度变化关于车削振动及车削参数预测模型和车削振动关于车削温度变化及车削参数预测模型,分析不同磨损状态刀具车削温度和车削振动的耦合特性。

Description

不同磨损刀具车削温度变化与车削振动耦合特性研究方法
技术领域
本发明属于车削加工领域,具备涉及一种不同刀具磨损状态下车削温度变化与车削振动耦合特性的研究方法。
背景技术
车削加工过程中,刀具与工件之间会产生相对摩擦,不仅导致刀具和工件的温度不断上升,而且还会产生严重的车削振动,加快刀具的磨损,而刀具的磨损会使车削温度上升速度加快,缩短刀具的寿命,降低加工工件的强度以及其表面的质量。在柔性生产系统中,如果不了解刀具磨损情况或者不能及时更换刀具都会导致生产效率的降低。现有的文献中有很多研究刀具的材料、几何参数、加工条件和车削参数等方面对车削温度的影响,或以车削振动为研究对象,通过试验和建立回归模型研究车削参数、车削力、机床结构和加工条件等因素对振动的影响,或通过试验研究不同的加工条件对刀具磨损的影响。但现有的研究中缺乏不同刀具磨损状态下车削参数、车削温度和车削振动三者的耦合关系的研究,尤其缺乏关于车削温度变化对车削振动影响的研究成果。
发明内容
本发明的目的在于提供一种不同刀具磨损状态下车削温度变化与车削振动耦合特性的研究方法,通过不同刀具磨损状态下的车削试验,得到相应车削参数下刀尖附近的车削温度和三向车削振动试验数据,基于灰色系统理论分别计算三种磨损状态刀具的(1)车削温度变化与车削参数的关联度;(2)车削振动与车削参数的关联度;(3)车削振动与车削温度变化的关联度。基于试验数据及最小二乘法建立(1)车削温度变化关于车削振动及车削参数预测模型,(2)车削振动关于车削温度变化及车削参数预测模型,分析对不同刀具磨损状态下车削温度变化与车削振动之间的耦合特性。
为实现上述发明目的,本发明的技术方案具体如下:
不同刀具磨损状态下车削温度变化与车削振动耦合特性的研究方法,包括以下步骤:
S1:制定车削试验方案,搭建试验系统;
S2:对不同磨损状态刀具进行车削试验,得到相应车削参数下刀尖附近的车削温度和三向车削振动数据;从车削温度和三向车削振动数据中提取出升温特征值和振动特征值;
S3:根据车削温度变化特征值、车削振动特征值和车削参数建立关于预测升温均值
Figure BDA0002767729800000021
拟合公式,计算不同磨损状态刀具的车削温度变化均值拟合值
Figure BDA0002767729800000022
与实际车削温度变化测量值
Figure BDA0002767729800000023
的相关性。根据车削温度变化特征值、车削振动特征值和车削参数建立各向振动特征值关于车削温度变化的拟合公式,计算不同磨损状态刀具的三向振动加速度均值拟合值
Figure BDA0002767729800000024
与实际测量值
Figure BDA0002767729800000025
的相关性。分析出不同磨损状态刀具对车削温度和车削振动及其相关性的影响。
作为优选,所述步骤S2具体包括:将各时刻实际测量温度与各次试验初始温度之差的平均值
Figure BDA0002767729800000026
作为温度信号分析的特征值,选择三向振动加速度轴向
Figure BDA0002767729800000027
径向
Figure BDA0002767729800000028
和切向
Figure BDA0002767729800000029
均方根值作为振动特征值。
作为优选,所述步骤S3包括:
S3-1:采用单一变量分析方法,对比分析出不同刀具磨损状态在相同水平主轴转速、进给速度和切削深度下对车削温度变化均值影响;对比分析出不同刀具磨损状态相同车削参数下对轴向、径向、切向的加速度均方根值影响;对比分析出不同刀具磨损状态车削参数对车削温度变化均值影响。
S3-2:建立车削温度变化关于车削振动及车削参数的拟合模型,分析车削振动及车削参数对车削温度变化影响;建立车削振动关于车削温度变化及车削参数的拟合模型,分析车削温度变化及车削参数对车削振动的影响。
作为优选,所述步骤S3-1包括:
根据试验设定的车削参数及在各次试验中的车削温度变化均值,计算出不同刀具磨损状态下车削温度变化与车削参数之间的灰色关联度,分析出何种车削参数对刀具温升的影响最显著;根据试验设定的车削参数及不同刀具磨损状态下各次试验中加速度均方根值,计算不同刀具磨损状态对应的加速度均方根值与车削参数之间的灰色关联度,分析出何种车削参数对加速度均方根值的影响最显著;基于灰色关联度分析出给定车削参数下,不同刀具磨损状态车削温度变化与车削三向振动之间的相关性。
作为优选,所述步骤3-2包括:
S5-1:以三向振动中某个方向的加速度均方根值
Figure BDA00027677298000000210
以及三个车削参数为自变量,建立预测车削温度变化均值
Figure BDA00027677298000000211
的拟合模型,分析车削振动及车削参数对车削温度变化影响:
Figure BDA0002767729800000031
其中,
Figure BDA0002767729800000032
为车削温度变化拟合值,C,x,y,z,w为拟合公式中的待定系数。
首先基于最小二乘法建立车削温度变化均值的回归模型,然后对车削温度变化关于车削参数和车削振动对的相关性进行分析。通过使各次试验中平均温度升高的拟合值与实测值的残差平方和∏为最小值,由下面的方程联立求解待定系数C,x,y,z,w:
Figure BDA0002767729800000033
Figure BDA0002767729800000034
将车削参数及某向加速度特征值、实测温升均值
Figure BDA0002767729800000035
代入到公式(2)~(4)中,通过在 Matlab里编程并即可求出回归系数C,x,y,z,w的值。
将回归系数带回到公式(1)分别得到与轴向振动、径向振动、切向振动对应的拟合温升均值,分别为
Figure BDA0002767729800000036
(轴向),
Figure BDA0002767729800000037
(径向),
Figure BDA0002767729800000038
(切向).
Figure BDA0002767729800000039
其中Ci,xi,yi,zi,wi,i=1,2,3,为与三向振动的回归系数。通过Matlab的 [R,P]=corrcoef(x,y)程序计算相关系数,其中R代表实测温升均值与拟合温升均值相关系数,P代表相关性为零的概率值;
S5-2:以车削温度变化均值
Figure BDA00027677298000000310
以及三个车削参数为自变量,建立预测车削振动加速度均方根值
Figure BDA00027677298000000311
的拟合模型,以分析车削温度变化及车削参数对车削振动的影响:
Figure BDA00027677298000000312
然后三向加速度均方根预测公式为:
Figure BDA0002767729800000041
其中,
Figure BDA0002767729800000042
为单向加速度的均方根值的拟合值,且分别对应于轴向振动,径向振动和切向振动。根据S5-1方法计算公式(7)中的回归系数 x′i,y′i,z′i,w′i,C′i(i=1,2,3)。
作为优选,同样的,计算预测值和实测三向加速度均方根值
Figure BDA0002767729800000043
的相关系数,并进行相关性分析。
与现有技术相比,本发明的有益效果:
本发明的一种不同刀具磨损状态下车削温度变化与车削振动耦合特性的研究方法,通过搭建车削温度与车削振动同步测试系统,采集三种磨损状态刀具刀尖附近的车削温度和车削振动,基于试验数据及灰色关联度理论,分析了不同磨损状态刀具车削温度变化与车削振动的耦合特性;通过灰色系统理论分别计算三种磨损状态刀具的车削温度变化与车削参数的关联度、车削振动与车削参数的关联度和车削温度变化与车削振动的关联度,并基于试验数据及最小二乘法分别建立车削温度变化关于车削振动及车削参数回归模型和车削振动关于车削温度变化及车削参数回归模型,分析车削温度变化与车削振动的互相关特性。
附图说明
图1是本发明的研究方法流程图;
图2是实验系统示意图;
图3是测试方案流程图;
图4是三种刀具车削温度的时域曲线;
图5是刀具D1,D2和D3在n=800r/min,vf=80mm/min,ap=0.3mm时三向振动加速度随时间的变化曲线;
图6不同水平主轴转速、进给速度和切削深度下刀具D1、D2和D3的车削温度变化对比图;
图7相同车削参数下各次试验中的轴向、径向、切向的加速度均方根值对比图。
具体实施方式:
下面将结合附图和具体实施例,对本发明的技术方案进行清楚、完整地描述。显然,所描述的实施例仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
实施例1:
通过不同刀具磨损状态下的车削试验,得到相应车削参数下刀尖附近的车削温度和三向车削振动,基于试验数据通过灰色系统理论计算三种刀具的车削参数、车削温度和车削振动之间的相对关联度。然后基于试验建立车削温度与车削振动的互相关模型,来分析不同刀具磨损状态下车削温度变化与车削振动、车削参数之间的耦合性。
研究流程如图1所示,包括:制定车削试验方案;搭建试验系统;对三种不同磨损状态刀具进行车削试验,得到相应车削参数下刀尖附近的车削温度和三向车削振动数据;从车削温度和三向车削振动数据中提取出温度特征值车削温度变化均值和振动特征值三向加速度均方根;根据温度特征值、振动特征值和车削参数建立关于预测车削温度变化均值
Figure BDA0002767729800000051
拟合公式,计算不同磨损状态刀具的车削温度变化均值拟合值
Figure BDA0002767729800000052
与实际车削温度变化测量值
Figure BDA0002767729800000053
的相关性;以车削温度变化均值
Figure BDA0002767729800000054
以及三个车削参数为自变量,建立预测车削振动加速度均方根值
Figure BDA0002767729800000055
的拟合模型,以分析车削温度变化及车削参数对车削振动的影响。进而分析不同磨损状态刀具对车削温度和车削振动及其相关性的影响。
1.试验设备
试验采用CJK6136-MATE-360×570数控车床,试验工件为45mm铝合金棒料。温度采集仪器采用OS523E-2手持红外测温仪、通过配套的OMEGASOFT OS5xx Access采集软件和计算机组成车削温度测量系统,可以实时测量和显示车削温度。采用YD-21 压电三向加速度传感器及WS-2402振动信号采集仪、DAQ信号采集处理软件和计算机组成振动测量与分析系统。试验采用D1(新刀片)、D2(中度磨损刀片)和D3(严重磨损刀片)三种磨损状态的硬质合金涂层刀片。试验系统如图2所示。
2.试验方案
测试方案的流程如图3所示,在此次的试验设计中,所用刀具为D1,D2和D3三种硬质合金刀具,分别对直径为45mm的铝棒进行干车削。
主轴转速设计800r/min,1200r/min,1600r/min,2000r/min这4个水平,进给速度设计3个水平,分别为40mm/min,80mm/min,160mm/min,切削深度0.3mm,0.5mm,0.8mm这3个水平。为了防止失误,每组车削参数做两次车削试验,每次车削试验时间为1分钟,具体参数如表1所示。
表1车削测试方案
Figure BDA0002767729800000061
3.试验步骤
车削时同时采集车削温度和刀具三向加速度,试验步骤如下:
(1)准备工件:选择粗加工后直径45mm的铝棒作为试验工件;
(2)连接和检测测温设备:将红外测温仪用专用数据线连接到个人计算机,打开红外测温仪开关,使其保持采集锁定状态。红外测温仪发射率为0.95。通过测温仪及计算机中测温软件采集记录刀尖附近的温度;
(3)连接和检测三向测振设备:将三向加速度传感器通过磁力座固定在刀柄下表面刀尖处,三向传感器的三个输出端分别接入WS-2402测振仪的9、10、11通道,通过专用数据线将WS-2402测振仪与计算机相连接.连通测振仪电源,打开计算机上相应的振动信号采集软件采集记录三向加速度信号;
(4)完成车削试验:按表1设定的车削参数进行车削试验,同步采集车削温度和车削三向振动加速度信号,并保存试验数据。试验结束后关闭电源。
4.试验结果与分析
通过时域分析得到各次试验中车削温度升高的平均值。图4代表刀具D1,D2和 D3在n=1200r/min,vf=80mm/min,ap=0.3mm时车削温度随时间的变化曲线。
利用Vib’SYS振动信号的采集和与分析系统得到三向振动加速度的时域曲线,通过时域分析,得到三向加速度、最大值、均方根等。图5代表刀具D1,D2和D3在n=800r/min,vf=80mm/min,ap=0.3mm时三向振动加速度随时间的变化曲线。
试验中发现:严重磨损刀具车削时的温度急剧升高,多次发生火星溅射现象并烧坏了传感器传输线,导致试验中断,故本论文主要给出D1和D2的实验数据,严重磨损刀具的试验数据只有一部分,如表2所示。每次试验中工件由于切削加工导致工件的直径不同,车削速度也不同,可以通过公式得到车削速度。
表2不同刀具的车削试验数据
Figure BDA0002767729800000071
Figure BDA0002767729800000081
由于车削试验是连续进行的,每次车削试验结束后刀具的温度都高于室内环境温度,因而每次车削试验中刀具的初始温度均不同。为此,提取各时刻实际测量温度与各次试验初始温度之差的平均值
Figure BDA0002767729800000082
作为温度信号分析的特征值,选择三向振动加速度轴向
Figure BDA0002767729800000083
径向
Figure BDA0002767729800000084
和切向
Figure BDA0002767729800000085
均方根值作为振动特征值,结合车削参数建立关于预测升温均值
Figure BDA0002767729800000086
拟合公式,通过Matlab计算D1、D2和D3三种刀具的温升均值拟合值
Figure BDA0002767729800000087
与实际温升测量值
Figure BDA0002767729800000088
的相关性。再以车削温度变化均值
Figure BDA0002767729800000089
以及三个车削参数为自变量,建立预测车削振动加速度均方根值
Figure BDA00027677298000000810
的拟合模型,以分析车削温度变化及车削参数对车削振动的影响进一步分析刀具磨损状态下对车削温度和车削振动及其相关性的影响。
5.不同刀具磨损状态下车削参数对车削温度变化影响
根据表2可以对比分析3种刀具D1,D2和D3在不同水平主轴转速、进给速度和切削深度下的车削温度变化均值变化规律,采用单一变量分析方法。图6(a)-(c)分别为不同水平主轴转速、进给速度和切削深度下D1、D2和D3的车削温度变化对比图。
根据图6(a)可知,随着主轴转速的增加,D3刀具温度急剧升高,远远大于D1和D2的温升值。D1和D2的温度升高对比不太显著。
根据图6(b)可知,随着进给速度的增加,D1和D2的车削温度变化均值曲线呈缓慢上升趋势,D3的温升均值远高于D1和D2。
根据图6(c)可知,3种刀具的温升均值随着切削深度的增加而急剧增加,且磨损越严重的刀具温度上升最快,最大的温升均值达到了153℃。
综上所述,D3刀具车削温度变化比D1和D2上升幅度更大,且在切削深度比较大时,磨损严重的刀具温升更加明显。
6.刀具磨损状态对车削振动的影响
对3种刀具在相同车削参数下各次试验中的轴向、径向、切向的加速度均方根值进行比较,如图7(a)-(c)所示。根据图7可知:3种刀具在轴向、径向和切向的加速度均方值曲线图的形状基本相似。大多数试验中D3和D1的加速度均方根大于D2的加速度均方根值。 D1刀具磨损量小,刀具过渡面与工件接触面积小,加快刀具磨损从而导致相对激烈的振动,D3因为磨损严重而导致严重的车削振动。
7.车削温度变化、车削振动、车削参数三者之间的灰色关联度分析
根据序列相对于起始点的变化速率之间的联系,灰色关联度分析可以判断序列间关系是否密切。变化速率越接近,相对关联度越大。结合试验数据,可以通过灰色关联度来寻找车削参数中对车削温度变化和车削振动影响最大的因素以及车削温度变化与车削振动二者间的相关度。首先需要通过灰色相对关联度理论进行计算得到灰色关联度值,与3种刀具相关的灰色关联度分别用RD1、RD2、RD3表示。先分析车削参数分别与车削温度变化及车削振动之间的灰色关联度,再分析车削温度变化与车削振动之间的灰色关联度。
7.1车削温度变化与车削参数的灰色关联度分析
根据试验设定的车削参数及3种刀具在各次试验中的车削温度变化均值,计算3种刀具的车削温度变化与车削参数之间的灰色关联度,分析哪种车削参数对刀具温升的影响最显著,计算结果如表3所示。
表3车削温度与车削参数的灰色关联度
Figure BDA0002767729800000091
根据表3:(1)3种刀具的车削温度变化均值
Figure BDA0002767729800000092
与车削参数的关联度系数均大于0.75,说明车削参数对车削温度的变化有很大影响;(2)相同车削参数下,可以比较3种刀具的与车削参数的相关性,由表3可知:RD3>RD2>RD1。即严重磨损刀具的平均车削温度变化和车削参数的关联度最高,因而车削参数的变化会显著影响严重磨损刀具的车削温度变化,但而对新刀具的车削温度变化的影响相对较最小;(3)在三个车削参数中,对D1和D2刀具,温升均值
Figure BDA0002767729800000101
与进给速度vf的关联度最高,与车削速度v的关联度最低,故进给速度对D1 和D2刀具车削温度的影响最大,车削速度对温度变化的影响最小。对D3刀具,温升均值
Figure BDA0002767729800000102
与车削深度ap与温升均值的关联度最高,与进给速度vf车削速度v的关联度最低。故车削深度对D3刀具车削温度影响最大,而进给速度对车削温度的影响最小的是进给速度。
7.2车削振动与车削参数的灰色关联度分析
根据试验设定的车削参数及3种刀具在各次试验中加速度均方根值,计算3种刀具的加速度均方根值与车削参数之间的灰色关联度,分析哪种车削参数对加速度均方根值的影响最显著,计算结果如表4所示。
表4加速度均方根值与车削参数之间的灰色关联度
Figure BDA0002767729800000103
根据表4:(1)相同车削参数下,可以比较3种刀具的车削振动与车削参数的相关性,可知:RD3>RD1>RD2,严重磨损刀具的加速度均方根值
Figure BDA0002767729800000104
和车削参数的关联度最高,而中度磨损的刀具的关联度最低;(2)对同一把刀具,可以比较哪个方向的振动加速度与同一车削参数的关联度最高,可知:三把刀具几乎轴向的振动加速度与同一车削参数的关联度最低,切向的关联度最高;(3)对同一把刀具某个方向的振动,还可以比较该向振动与哪个车削参数的关联度最高。可知:对于D1、D2刀具,三向振动与车削速度v的关联度最高,与进给速度vf的关联度最低。而D3刀具,其轴向的振动、径向的振动都与进给速度vf的关联度最高,与车削深度ap的关联度最低,切向的振动都与车削速度v的关联度最高,与进给速度vf的关联度最低。
7.3车削温度变化与车削振动之间的灰色关联度分析
基于灰色关联度分析给定车削参数下,3种刀具车削温度变化与车削三向振动之间的相关性,计算结果如表5所示。
表5车削温度与车削振动的灰色关联度
Figure BDA0002767729800000111
根据表5可知:(1)相同车削参数下,比较3种刀具的车削温度变化均值
Figure BDA0002767729800000112
与加速度均方根值
Figure BDA0002767729800000113
的相关性,可知:RD3>RD1>RD2。严重磨损刀具的车削温度变化均值
Figure BDA0002767729800000114
和加速度均方根值
Figure BDA0002767729800000115
的关联度最高,因而车削振动的变化会显著影响严重磨损刀具的车削温度,但是对中度磨损的刀具的影响很小。(2)对同一把刀具,可以比较哪个方向的振动加速度与车削温度变化的关联度最高,可知:对于D2、D3刀具,车削温度变化
Figure BDA0002767729800000116
与径向振动
Figure BDA0002767729800000117
的关联度最高,与轴向振动
Figure BDA0002767729800000118
的关联度最低。而D1刀具,车削温度变化
Figure BDA0002767729800000119
与切向振动
Figure BDA00027677298000001110
的关联度最高,与轴向振动
Figure BDA00027677298000001111
的关联度最低
8.不同磨损状态的刀具的车削温度变化与车削振动的耦合特性
8.1车削温度变化关于车削振动及车削参数的拟合模型
为了分析车削振动及车削参数对车削温度变化影响,下面以三向振动中某个方向的加速度均方根值
Figure BDA00027677298000001112
以及三个车削参数为自变量,建立预测车削温度变化均值
Figure BDA00027677298000001113
的拟合模型:
Figure BDA00027677298000001114
其中,
Figure BDA00027677298000001115
为车削温度变化拟合值,C,x,y,z,w为拟合公式中的待定系数。
首先基于最小二乘法建立车削温度变化均值的回归模型,然后对车削温度变化关于车削参数和车削振动对的相关性进行分析。通过使各次试验中平均温度升高的拟合值与实测值的残差平方和∏为最小值,由下面的方程联立求解待定系数C,x,y,z,w:
Figure BDA00027677298000001116
Figure BDA0002767729800000121
Figure BDA0002767729800000122
Figure BDA0002767729800000123
Figure BDA0002767729800000124
Figure BDA0002767729800000125
Figure BDA0002767729800000126
可以将等式(8)-(12)转换为有关于x,y,z,w,η的线性方程,然后把三个车削参数,实测车削温度变化均值
Figure BDA0002767729800000127
和单向加速度均方根值
Figure BDA0002767729800000128
带入到公式(8)-(12),因此回归系数x, y,z,w,C(C=eη)可以通过MATLAB软件编程求解。
然后将获得的回归系数C,x,y,z,w代入公式(2)中,可以得出拟合温度升高
Figure BDA0002767729800000129
(轴向),
Figure BDA00027677298000001210
(径向),
Figure BDA00027677298000001211
(切向),分别对应于轴向振动,径向振动和切向振动,我们有
Figure BDA0002767729800000131
其中Ci,xi,yi,zi,wi,i=1,2,3,为与三向振动的回归系数。通过Matlab的[R,P]=corrcoef(x,y) 程序计算相关系数,其中R代表实测温升均值与拟合温升均值相关系数,P代表相关性为零的概率值。3种刀具的计算结果如表6所示。
表6三种刀具的回归模型参数和相关系数
Figure BDA0002767729800000132
根据表6可知:D1的相关系数要比D3和D2的相关系数要高,而且D1相关性为0的概率P值小于D2和D3,这代表D1的车削振动与车削温度变化有着非常好的相关性,可以非常好的预测刀具的实际车削温度。三把刀具径向振动拟合的相关性较高,说明刀具的车削温度变化会对车削径向振动产生一定的影响。
8.2车削振动关于车削温度变化及车削参数的拟合模型
为了分析车削温度变化及车削参数对车削振动的影响,下文以车削温度变化均值
Figure BDA0002767729800000133
以及三个车削参数为自变量,建立预测车削振动加速度均方根值
Figure BDA0002767729800000134
的拟合模型:
Figure BDA0002767729800000135
然后三向加速度均方根预测公式为:
Figure BDA0002767729800000136
其中,
Figure BDA0002767729800000137
为单向加速度的均方根值的拟合值,且分别对应于轴向振动,径向振动和切向振动。根据8.2节介绍的方法计算公式(7)中的回归系数 x′i,y′i,z′i,w′i,C′i(i=1,2,3)。
预测值和实测三向加速度均方根值
Figure BDA0002767729800000141
的相关系数,计算结果如表7所示。
表7三种刀具的回归模型参数和相关系数
Figure BDA0002767729800000142
根据表7可知:D1和D2刀具的相关系数比较低,说明车削温度的变化对车削振动的影响不显著;而对于D3刀具,相关性系数均大于0.6,且以轴向振动拟合的相关性较高,说明严重磨损刀具的车削温度变化会对车削轴向振动产生一定的影响。
总结
对三种不同磨损状态的刀具在完成了相同车削参数、相同工件尺寸下的车削实验,同步采集刀尖附近的车削温度及车削振动,研究不同磨损状态下刀具车削温度变化与车削振动的相互耦合特性,研究表明:
(1)严重磨损刀具的车削温度与车削参数的关联度最高,新刀具的关联度最低,说明磨损越严重的刀具受到车削温度影响越显著。
(2)车削参数对严重磨损刀具的车削振动的影响最显著,而中度磨损刀具的受到影响最小。
(3)基于车削振动及车削参数建立了车削温度升高均值的拟合模型,该模型可以较好地预测刀具在给定车削参数下的车削温度变化均值;同时也基于车削温度变化及车削参数建立了车削振动的预测模型,通过相关性计算表明D1和D2刀具的车削温度变化对其车削振动的影响很弱。
(4)实验数据中车削温度变化
Figure BDA0002767729800000143
存在急剧上升和急剧下降的现象,且特别是磨损严重的刀具采用较大车削参数时,车削温度过高,出现冒烟现象,试验被迫停止,可能原因为刀具车削工件时切屑没有及时脱落导致积屑从而使得车削温度变化较大。
(5)新刀与严重磨损的刀具其车削温度变化与车削振动的相关性高于中度磨损刀具。

Claims (2)

1.不同磨损刀具车削温度变化与车削振动耦合特性研究方法,其特征在于,包括以下步骤:
S1:制定车削试验方案,搭建试验系统;
S2:对不同磨损状态刀具进行车削试验,得到相应车削参数下刀尖附近的车削温度和三向车削振动数据;从车削温度和三向车削振动数据中提取出升温特征值和振动特征值;
S3:分析不同磨损状态刀具对车削温度和车削振动及其相关性的影响;
根据车削温度变化特征值、车削振动特征值和车削参数建立关于预测车削温度变化均值
Figure FDA0003526712370000011
的拟合公式,计算不同磨损状态刀具的车削温度变化均值拟合值
Figure FDA0003526712370000012
与实际车削温度变化均值
Figure FDA0003526712370000013
的相关性;
根据车削温度变化特征值、车削振动特征值和车削参数建立各向振动特征值关于车削温度变化的拟合公式,计算不同磨损状态刀具的三向振动加速度均值拟合值
Figure FDA0003526712370000014
与实际测量值
Figure FDA0003526712370000015
的相关性;
所述步骤S2具体包括:将车削温度变化均值
Figure FDA0003526712370000016
作为温度信号分析的特征值,
Figure FDA0003526712370000017
是各时刻实际测量温度与各次试验初始温度之差的平均值,选择三向振动加速度轴向
Figure FDA0003526712370000018
径向
Figure FDA0003526712370000019
和切向
Figure FDA00035267123700000110
均方根值作为振动特征值;
所述步骤S3包括:
S3-1:采用单一变量分析方法,对比分析出不同刀具磨损状态在相同水平主轴转速、进给速度和切削深度下对车削温度变化均值影响;对比分析出不同刀具磨损状态在相同车削参数下对轴向、径向、切向的加速度均方根值影响;对比分析出不同刀具磨损状态车削振动对车削温度变化均值影响;
S3-2:建立车削温度变化关于车削振动及车削参数的拟合模型,分析车削振动及车削参数对车削温度变化影响;建立车削振动关于车削温度变化及车削参数的拟合模型,分析车削温度变化及车削参数对车削振动的影响;
所述步骤S3-1包括:
根据试验设定的车削参数及在各次试验中的车削温度变化均值,计算出不同刀具磨损状态下车削温度变化与车削参数之间的灰色关联度,分析出何种车削参数对刀具温升的影响最显著;根据试验设定的车削参数及不同刀具磨损状态下各次试验中加速度均方根值,计算不同刀具磨损状态对应的加速度均方根值与车削参数之间的灰色关联度,分析出何种车削参数对加速度均方根值的影响最显著;基于灰色关联度分析出给定车削参数下,不同刀具磨损状态车削温度变化与车削三向振动之间的相关性;
所述步骤3-2包括:
S5-1:以三向振动中某个方向的加速度均方根值
Figure FDA0003526712370000021
以及三个车削参数为自变量,建立预测车削温度变化均值
Figure FDA0003526712370000022
的拟合模型,分析车削振动及车削参数对车削温度变化影响:
Figure FDA0003526712370000023
其中,
Figure FDA0003526712370000024
为车削温度变化均值拟合值,C,x,y,z,w为拟合公式中的待定系数,v是车削速度,vf是进给速度,ap是车削深度;首先基于最小二乘法建立车削温度变化均值的回归模型,然后对车削温度变化关于车削参数和车削振动对的相关性进行分析;通过使各次试验中车削温度变化均值拟合值与实测车削温度变化均值的残差平方和∏为最小值,由下面的方程联立求解待定系数C,x,y,z,w:
Figure FDA0003526712370000025
Figure FDA0003526712370000026
将车削参数及某向加速度特征值、实测车削温度变化均值
Figure FDA0003526712370000027
代入到公式(2)~(4)中,通过在Matlab里编程并即可求出回归系数C,x,y,z,w的值;将回归系数带回到公式(2)分别得到与轴向振动、径向振动、切向振动对应的车削温度变化均值拟合值
Figure FDA0003526712370000028
Figure FDA0003526712370000029
其中Ci,xi,yi,zi,wi是三向振动的回归系数,i=1,2,3;通过Matlab的[R,P]=corrcoef(x,y)程序计算相关系数,其中R代表实测车削温度变化均值与车削温度变化均值拟合值相关系数,P代表相关性为零的概率值;
S5-2:以车削温度变化均值
Figure FDA0003526712370000031
以及三个车削参数为自变量,建立预测车削振动加速度均方根值
Figure FDA0003526712370000032
的拟合模型,以分析车削温度变化及车削参数对车削振动的影响:
Figure FDA0003526712370000033
三向加速度均方根预测公式为:
Figure FDA0003526712370000034
其中,
Figure FDA0003526712370000035
分别为轴向、径向、切向的加速度均方根值的拟合值,且分别对应于轴向振动,径向振动和切向振动;根据S5-1方法计算公式(7)中的回归系数x′i,y′i,z′i,w′i,C′i,i=1,2,3。
2.根据权利要求1所述的研究方法,其特征在于,所述步骤S5-2还包括,计算预测值和实测三向加速度均方根值
Figure FDA0003526712370000036
的相关系数并进行相关性分析。
CN202011238945.9A 2020-11-09 2020-11-09 不同磨损刀具车削温度变化与车削振动耦合特性研究方法 Active CN112380646B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011238945.9A CN112380646B (zh) 2020-11-09 2020-11-09 不同磨损刀具车削温度变化与车削振动耦合特性研究方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011238945.9A CN112380646B (zh) 2020-11-09 2020-11-09 不同磨损刀具车削温度变化与车削振动耦合特性研究方法

Publications (2)

Publication Number Publication Date
CN112380646A CN112380646A (zh) 2021-02-19
CN112380646B true CN112380646B (zh) 2022-05-03

Family

ID=74579014

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011238945.9A Active CN112380646B (zh) 2020-11-09 2020-11-09 不同磨损刀具车削温度变化与车削振动耦合特性研究方法

Country Status (1)

Country Link
CN (1) CN112380646B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115082433B (zh) * 2022-07-21 2022-12-09 深圳市信润富联数字科技有限公司 微铣削刀工作参数确定方法、装置、电子设备及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107066775A (zh) * 2016-05-30 2017-08-18 江苏师范大学 一种刀具车削温升均值的预测方法
CN109746465A (zh) * 2018-09-26 2019-05-14 江苏师范大学 一种车削振动-车削变形-车削温度实时监测与分析系统
CN110728049A (zh) * 2019-10-09 2020-01-24 江苏师范大学 一种刀具车削温度变化均值的组合预测模型建立方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2311265C1 (ru) * 2006-03-31 2007-11-27 Рыбинская государственная авиационная технологическая академия им. П.А. Соловьева Способ стабилизации температуры резания при точении
WO2018176337A1 (zh) * 2017-03-30 2018-10-04 深圳市红昌机电设备有限公司 车削床的控制方法及系统
CN209550635U (zh) * 2018-09-26 2019-10-29 江苏师范大学 一种车削振动-车削变形-车削温度实时监测与分析系统

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107066775A (zh) * 2016-05-30 2017-08-18 江苏师范大学 一种刀具车削温升均值的预测方法
CN109746465A (zh) * 2018-09-26 2019-05-14 江苏师范大学 一种车削振动-车削变形-车削温度实时监测与分析系统
CN110728049A (zh) * 2019-10-09 2020-01-24 江苏师范大学 一种刀具车削温度变化均值的组合预测模型建立方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
刀具磨损状态对车削温度影响的试验研究;袁冠雷等;《煤矿机械》;20151130;第36卷(第11期);107-109 *
基于灰色关联分析的微小孔振动钻削参数优化;鄢春艳等;《机械设计与制造》;20151231(第12期);19-22 *
基于试验设计与灰色关联分析的数控车削参数优化;王灵玲等;《工艺与检测》;20150731(第7期);正文第2节 *

Also Published As

Publication number Publication date
CN112380646A (zh) 2021-02-19

Similar Documents

Publication Publication Date Title
Yang et al. Research on the milling tool wear and life prediction by establishing an integrated predictive model
Orhan et al. Tool wear evaluation by vibration analysis during end milling of AISI D3 cold work tool steel with 35 HRC hardness
Şahinoğlu et al. Optimization of cutting parameters with respect to roughness for machining of hardened AISI 1040 steel
Rao et al. Influence of cutting parameters on cutting force and surface finish in turning operation
Shankar et al. Tool condition monitoring while using vegetable based cutting fluids during milling of inconel 625
Arslan et al. Tool condition monitoring in turning using statistical parameters of vibration signal
CN112757052B (zh) 不同磨损刀具的车削热与多元影响因子的相关性分析方法
Chibane et al. Optimal milling conditions for carbon/epoxy composite material using damage and vibration analysis
Klaic et al. Tool wear classification using decision treesin stone drilling applications: A preliminary study
Bhuiyan et al. Tool condition monitoring using acoustic emission and vibration signature in turning
CN112380646B (zh) 不同磨损刀具车削温度变化与车削振动耦合特性研究方法
Cui et al. Identification of the optimum cutting parameters in intermittent hard turning with specific cutting energy, damage equivalent stress, and surface roughness considered
Pereira et al. Topographic wear monitoring of the interface tool/workpiece in milling AISI H13 steel
Zhong et al. Correlation analysis of cutting force and acoustic emission signals for tool condition monitoring
CN111761410A (zh) 一种基于振动信号的刀具磨损状态监测装置及其方法
CN109357826B (zh) 一种车削刀具-工件非线性振动系统耦合特征的研究方法
Mutalib et al. Evaluation of tool wear and machining performance by analyzing vibration signal in friction drilling
Xiuli et al. Research on predictive model surface roughness in high speed milling for aluminum alloy 7050-T7451
CN112364454B (zh) 一种基于模态质量分布矩阵的机床薄弱构件识别方法
Li et al. Tool condition monitoring based on vibration signal from an on-rotor sensor in CNC turning process
Roszkowski et al. Study on the impact of cutting tool wear on machine tool energy consumption
Prasad et al. Analysis of real-time vibration assisted tool condition monitoring in drilling
Khaisar et al. Prediction of tool wear during turning of EN9 work material against coated carbide insert using vibration signal processing
Shah et al. Tool control monitoring applied to drilling
Zhang et al. Face milling tool wear condition monitoring based on wavelet transform and Shannon entropy

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
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20240131

Address after: 250000 Tianchen Road, High tech Zone, Jinan City, Shandong Province

Patentee after: Tianhong (Jinan) Intelligent Equipment Industry Research Co.,Ltd.

Country or region after: China

Address before: No. 101, Shanghai Road, Tongshan District, Xuzhou City, Jiangsu Province

Patentee before: Jiangsu Normal University

Country or region before: China

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20240401

Address after: Room 434, Comprehensive Office Building, Rendong Pharmaceutical Industrial Park, No. 1728 Tian'an 1st Road, Xingren Street, High tech Zone, Zaozhuang City, Shandong Province, 277000

Patentee after: Shandong Yijite CNC Machine Tool Co.,Ltd.

Country or region after: China

Address before: 250000 Tianchen Road, High tech Zone, Jinan City, Shandong Province

Patentee before: Tianhong (Jinan) Intelligent Equipment Industry Research Co.,Ltd.

Country or region before: China