CN112818477B - 一种整体式平底立铣刀的刀具失效界限图建立方法及系统 - Google Patents

一种整体式平底立铣刀的刀具失效界限图建立方法及系统 Download PDF

Info

Publication number
CN112818477B
CN112818477B CN202110002359.2A CN202110002359A CN112818477B CN 112818477 B CN112818477 B CN 112818477B CN 202110002359 A CN202110002359 A CN 202110002359A CN 112818477 B CN112818477 B CN 112818477B
Authority
CN
China
Prior art keywords
wear
cutting
model
face
end mill
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
CN202110002359.2A
Other languages
English (en)
Other versions
CN112818477A (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.)
Harbin University of Science and Technology
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN202110002359.2A priority Critical patent/CN112818477B/zh
Publication of CN112818477A publication Critical patent/CN112818477A/zh
Application granted granted Critical
Publication of CN112818477B publication Critical patent/CN112818477B/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
    • 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/08Thermal analysis or thermal optimisation
    • 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)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Numerical Control (AREA)
  • Milling Processes (AREA)

Abstract

本发明涉及一种整体式平底立铣刀的刀具失效界限图建立方法及系统,根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型、前刀面月牙洼磨损深度变化率模型和临界断裂时许用应力模型,获得满足条件的第一切削参数、第二切削参数和第三切削参数,在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。本发明通过结合三种失效形式不同的失效标准,建立刀具失效界限图,为刀具在铣削钛合金时,提供一个安全区域,避免刀具过快失效。

Description

一种整体式平底立铣刀的刀具失效界限图建立方法及系统
技术领域
本发明涉及钛合金切削加工领域,特别是涉及一种整体式平底立铣刀的刀具失效界限图建立方法及系统。
背景技术
钛合金Ti6Al4V因其具有密度小、强度-重量比高、耐高温以及耐腐蚀等一系列优良特性,而广泛应用在航空、航天、能源以及生物医疗等制造业领域。为满足上述领域对于加工效率和工件表面加工质量的更高应用要求,采用高速铣削加工方式显得尤为重要。但是由于自身导热率较低以及高温切削加工情况下化学活性较高导致其成为典型的难加工材料之一。硬质合金立铣刀因其硬度高、耐磨性好和抗氧化性等优良性能而广泛用于高速铣削加工钛合金材料。然而,在加工钛合金过程中,硬质合金立铣刀极易出现前后刀面磨损及切削刃破损等失效现象,将导致刀具寿命缩短、加工效率降低以及加工表面质量变差等问题。
因此,钛合金加工用硬质合金立铣刀失效的机理与解析建模预测研究显得尤为重要。
发明内容
本发明的目的是提供一种整体式平底立铣刀的刀具失效界限图建立方法及系统,为切削刀具提供了一个安全的操作区域,以保证安全区域内的切削用量组合不会在预定更换刀具时间之前发生故障。
为实现上述目的,本发明提供了如下方案:
一种整体式平底立铣刀的刀具失效界限图建立方法,所述方法包括:
获取预设参数范围内整体式平底立铣刀的多个初始切削参数;每个初始切削参数包括切削厚度和切削速度,每个切削参数设置不同的切削厚度和切削速度;
根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型;
确定每个后刀面磨损宽度变化率模型的磨损宽度的最大值,并筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型,获得筛选出的后刀面磨损宽度变化率模型对应的初始切削参数,确定为第一切削参数;
根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型;
确定每个前刀面月牙洼磨损深度模型的磨损深度的最大值,并筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型,获得筛选出的前刀面月牙洼磨损深度模型对应的初始切削参数,确定为第二切削参数;
根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型;
确定每个临界断裂时许用应力模型的临界断裂时许用应力,并筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型,获得筛选出的临界断裂时许用应力模型对应的初始切削参数,确定为第三切削参数;
在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线;所述同一坐标系的横坐标为切削速度,纵坐标为切削宽度;
将所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。
可选的,所述根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型,具体包括:
根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000031
Figure BDA0002882022360000032
建立基于几何的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000033
将基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型的体积等价于基于几何的整体式平底立铣刀的后刀面磨损体积模型的体积,获得每个初始切削参数的整体式平底立铣刀的后刀面磨损宽度变化率模型为
Figure BDA0002882022360000034
其中,
Figure BDA0002882022360000035
为第z个基于磨损机理的整体式平底立铣刀的后刀面磨损体积,σn-f(x")为后刀面磨损带位置x″处的正应力,w为切屑宽度,vz为第z个初始切削参数的切削速度,Δtz为第z个切削有效时间,
Figure BDA0002882022360000036
fz为每齿进给量,R为刀具半径,
Figure BDA0002882022360000037
为第z个初始切削参数的切削宽度,H0为硬质合金硬度,a为与温度有关的系数,K为切削刃上剪应力和工件上剪切流动应力的比值,
Figure BDA0002882022360000038
为第z个刀具后刀面和工件接触面的温度,Ht为铣刀表面硬度值,n为三体磨粒磨损常数,VB为铣刀后刀面磨损宽度,d为粘结磨损系数,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数,
Figure BDA0002882022360000039
为第z个基于几何的整体式平底立铣刀的后刀面磨损体积,λs为刀具刃倾角,
Figure BDA00028820223600000310
为第z个后刀面磨损宽度变化率,ΔVBz为第z个铣刀后刀面磨损宽度变化量,α为后角,γn为前角。
可选的,所述根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型,具体包括:
建立整体式平底立铣刀的后刀面磨损带正应力分布模型为
Figure BDA0002882022360000041
根据每个初始切削参数的切削速度,建立整体式平底立铣刀的后刀面磨损带温度分布函数模型为
Figure BDA0002882022360000042
根据每个初始切削参数的切削厚度,计算一个刀齿切入到切出过程的切削有效时间为
Figure BDA0002882022360000043
根据后刀面磨损带正应力分布模型、后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于磨粒磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面磨粒磨损体积模型为
Figure BDA0002882022360000044
根据后刀面磨损带正应力分布模型、后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于粘结磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面粘结磨损体积模型为
Figure BDA0002882022360000045
根据后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于扩散磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面扩散磨损体积模型为
Figure BDA0002882022360000046
根据每个切削参数建立的后刀面磨粒磨损体积模型、后刀面粘结磨损体积模型和后刀面扩散磨损体积模型,分别建立每个初始切削参数的基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000051
Figure BDA0002882022360000052
其中,σtip-flank为微元刀尖上靠近后刀面磨损带的正应力,
Figure BDA0002882022360000053
ρ是材料密度,φn为剪切角,ηw为后刀面磨损带上摩擦力的滑移线场角,VBCR为弹性接触区域的宽度,μc为后刀面摩擦系数,k是材料热导率,c是材料比热容,i为第一累加运算符号,j为第二累加运算符号,
Figure BDA0002882022360000054
Figure BDA0002882022360000055
均为排列组合运算,T0为室温,
Figure BDA0002882022360000056
为泰勒展开式系数,t为切削时间,
Figure BDA0002882022360000057
为第z个后刀面磨粒磨损体积,Cflank-abrasion为无量纲磨粒磨损系数,
Figure BDA0002882022360000058
为第z个粘结磨损导致的铣刀后刀面磨损体积,Cflank-adhension为无量纲粘结磨损系数,
Figure BDA0002882022360000059
为第z个扩散磨损机理导致的铣刀后刀面磨损体积,Cflank-diffusion为无量纲扩散磨损系数,
Figure BDA00028820223600000510
为第z个刀具后刀面总磨损体积。
可选的,所述根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型,具体包括:
根据每个初始切削参数的切削厚度,分别建立整体式平底立铣刀的前刀面法向接触应力分布模型为
Figure BDA00028820223600000511
根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面温度分布函数模型为
Figure BDA00028820223600000512
根据每个初始切削参数的切削速度,基于磨粒磨损机理,分别建立整体式平底立铣刀的前刀面磨粒磨损率模型为
Figure BDA00028820223600000513
根据每个前刀面法向接触应力分布模型和每个前刀面温度分布函数模型,基于粘结磨损机理,分别建立整体式平底立铣刀的前刀面粘结磨损率模型为
Figure BDA00028820223600000514
根据每个前刀面温度分布函数模型和每个前刀面温度分布函数模型对应的切削速度,基于扩散磨损机理,分别建立整体式平底立铣刀的前刀面扩散磨损率模型为
Figure BDA0002882022360000061
根据每个切削参数建立的前刀面磨粒磨损率模型、前刀面粘结磨损率模型和前刀面扩散磨损率模型,分别建立每个初始切削参数的整体式平底立铣刀的前刀面月牙洼深度变化率模型为
Figure BDA0002882022360000062
其中,
Figure BDA0002882022360000063
为第z个初始切削参数的切削厚度对应的前刀面法向接触应力,σtip-rake是刀尖处最大法向接触应力,
Figure BDA0002882022360000064
是第z个初始切削参数的切削厚度对应的刀屑接触长度,
Figure BDA0002882022360000065
Figure BDA0002882022360000066
为第z个初始切削参数的切削厚度对应的最大切屑厚度,
Figure BDA0002882022360000067
φn为剪切角,γn为前角,βn为摩擦角,ξ1为相关幂指数,ηc为切屑流出角,R为刀具半径,f为每齿进给量,ξ为应力指数参数,x为刀屑接触距离,
Figure BDA0002882022360000068
为第z个初始切削参数对应的前刀面温度,
Figure BDA0002882022360000069
是刀屑平均摩擦系数,
Figure BDA00028820223600000610
为第z个初始切削参数的切削速度对应的切屑速度,
Figure BDA00028820223600000611
λs为刀具刃倾角,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容,T0为室温,
Figure BDA00028820223600000612
为第z个磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率,Qtool-chip为刀屑接触面之间的接触载荷,Hrake-face为刀具前刀面的硬度,θabrasion为磨粒平均粗糙角,
Figure BDA00028820223600000613
为第z个粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率,N1为第一常数,N2为第二常数,
Figure BDA00028820223600000614
为第z个扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率,Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数,KT(x,t)为第z个初始切削参数的前刀面月牙洼深度变化率,t为切削时间。
可选的,所述根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型,具体包括:
根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型为
Figure BDA0002882022360000071
其中,
Figure BDA0002882022360000072
为第z个临界断裂时许用应力,Fz为第z个初始切削参数的切削速度对应的导致刀齿破损的冲击力,
Figure BDA0002882022360000073
fz为每齿进给量,ρ为材料密度,vz为切削速度,ap是切削深度,mA为断裂面积相关比系数,VB为铣刀后刀面磨损宽度,τ为许用应力,
Figure BDA0002882022360000074
为第z个初始切削参数的切削厚度对应的损伤弹性模量值,γ为硬质合金的切应变,bv为泊松比。
可选的,所述将所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,具体包括:
确定所述切削刃破损失效曲线上最小切削速度所在的第一坐标点,建立经过所述第一坐标点且平行于纵坐标轴的第一直线;
确定所述后刀面失效曲线上最大切削速度所在的第二坐标点,建立经过所述第二坐标点且平行于横坐标轴的第二直线;
将所述第一直线、所述第二直线、所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线形成的封闭区域确定为安全切削区域。
一种整体式平底立铣刀的刀具失效界限图建立系统,所述系统包括:
初始切削参数获取模块,用于获取预设参数范围内整体式平底立铣刀的多个初始切削参数;每个初始切削参数包括切削厚度和切削速度,每个切削参数设置不同的切削厚度和切削速度;
后刀面磨损宽度变化率模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型;
第一切削参数确定模块,用于确定每个后刀面磨损宽度变化率模型的磨损宽度的最大值,并筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型,获得筛选出的后刀面磨损宽度变化率模型对应的初始切削参数,确定为第一切削参数;
前刀面月牙洼磨损深度变化率模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型;
第二切削参数确定模块,用于确定每个前刀面月牙洼磨损深度模型的磨损深度的最大值,并筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型,获得筛选出的前刀面月牙洼磨损深度模型对应的初始切削参数,确定为第二切削参数;
临界断裂时许用应力模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型;
第三切削参数确定模块,用于确定每个临界断裂时许用应力模型的临界断裂时许用应力,并筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型,获得筛选出的临界断裂时许用应力模型对应的初始切削参数,确定为第三切削参数;
失效曲线建立模块,用于在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线;所述同一坐标系的横坐标为切削速度,纵坐标为切削宽度;
刀具失效界面图获得模块,用于将所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。
可选的,所述后刀面磨损宽度变化率模型建立模块,具体包括:
基于磨损机理的后刀面磨损体积模型建立子模块,用于根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000081
Figure BDA0002882022360000082
基于几何的后刀面磨损体积模型建立子模块,用于建立基于几何的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000083
后刀面磨损宽度变化率模型获得子模块,用于将基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型的体积等价于基于几何的整体式平底立铣刀的后刀面磨损体积模型的体积,获得每个初始切削参数的整体式平底立铣刀的后刀面磨损宽度变化率模型为
Figure BDA0002882022360000091
其中,
Figure BDA0002882022360000092
为第z个基于磨损机理的整体式平底立铣刀的后刀面磨损体积,σn-f(x")为后刀面磨损带位置x″处的正应力,w为切屑宽度,vz为第z个初始切削参数的切削速度,Δtz为第z个切削有效时间,
Figure BDA0002882022360000093
fz为每齿进给量,R为刀具半径,
Figure BDA0002882022360000094
为第z个初始切削参数的切削宽度,H0为硬质合金硬度,a为与温度有关的系数,K为切削刃上剪应力和工件上剪切流动应力的比值,
Figure BDA0002882022360000095
为第z个刀具后刀面和工件接触面的温度,Ht为铣刀表面硬度值,n为三体磨粒磨损常数,VB为铣刀后刀面磨损宽度,d为粘结磨损系数,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数,
Figure BDA0002882022360000096
为第z个基于几何的整体式平底立铣刀的后刀面磨损体积,λs为刀具刃倾角,
Figure BDA0002882022360000097
为第z个后刀面磨损宽度变化率,ΔVBz为第z个铣刀后刀面磨损宽度变化量,α为后角,γn为前角。
可选的,所述前刀面月牙洼磨损深度变化率模型建立模块,具体包括:
前刀面法向接触应力分布模型建立子模块,用于根据每个初始切削参数的切削厚度,分别建立整体式平底立铣刀的前刀面法向接触应力分布模型为
Figure BDA0002882022360000098
前刀面温度分布函数模型建立子模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面温度分布函数模型为
Figure BDA0002882022360000099
前刀面磨粒磨损率模型建立子模块,用于根据每个初始切削参数的切削速度,基于磨粒磨损机理,分别建立整体式平底立铣刀的前刀面磨粒磨损率模型为
Figure BDA0002882022360000101
前刀面粘结磨损率模型建立子模块,用于根据每个前刀面法向接触应力分布模型和每个前刀面温度分布函数模型,基于粘结磨损机理,分别建立整体式平底立铣刀的前刀面粘结磨损率模型为
Figure BDA0002882022360000102
前刀面扩散磨损率模型建立子模块,用于根据每个前刀面温度分布函数模型和每个前刀面温度分布函数模型对应的切削速度,基于扩散磨损机理,分别建立整体式平底立铣刀的前刀面扩散磨损率模型为
Figure BDA0002882022360000103
前刀面月牙洼深度变化率模型建立子模块,用于根据每个切削参数建立的前刀面磨粒磨损率模型、前刀面粘结磨损率模型和前刀面扩散磨损率模型,分别建立每个初始切削参数的整体式平底立铣刀的前刀面月牙洼深度变化率模型为
Figure BDA0002882022360000104
其中,
Figure BDA0002882022360000105
为第z个初始切削参数的切削厚度对应的前刀面法向接触应力,σtip-rake是刀尖处最大法向接触应力,
Figure BDA0002882022360000106
是第z个初始切削参数的切削厚度对应的刀屑接触长度,
Figure BDA0002882022360000107
Figure BDA0002882022360000108
为第z个初始切削参数的切削厚度对应的最大切屑厚度,
Figure BDA0002882022360000109
φn为剪切角,γn为前角,βn为摩擦角,ξ1为相关幂指数,ηc为切屑流出角,R为刀具半径,f为每齿进给量,ξ为应力指数参数,x为刀屑接触距离,
Figure BDA00028820223600001010
为第z个初始切削参数对应的前刀面温度,
Figure BDA00028820223600001011
是刀屑平均摩擦系数,
Figure BDA00028820223600001012
为第z个初始切削参数的切削速度对应的切屑速度,
Figure BDA00028820223600001013
λs为刀具刃倾角,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容,T0为室温,
Figure BDA0002882022360000111
为第z个磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率,Qtool-chip为刀屑接触面之间的接触载荷,Hrake-face为刀具前刀面的硬度,θabrasion为磨粒平均粗糙角,
Figure BDA0002882022360000112
为第z个粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率,N1为第一常数,N2为第二常数,
Figure BDA0002882022360000113
为第z个扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率,Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数,KT(x,t)为第z个初始切削参数的前刀面月牙洼深度变化率,t为切削时间。
可选的,所述临界断裂时许用应力模型建立模块,具体包括:
临界断裂时许用应力模型建立子模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型为
Figure BDA0002882022360000114
其中,
Figure BDA0002882022360000115
为第z个临界断裂时许用应力,Fz为第z个初始切削参数的切削速度对应的导致刀齿破损的冲击力,
Figure BDA0002882022360000116
fz为每齿进给量,ρ为材料密度,vz为切削速度,ap是切削深度,mA为断裂面积相关比系数,VB为铣刀后刀面磨损宽度,τ为许用应力,
Figure BDA0002882022360000117
为第z个初始切削参数的切削厚度对应的损伤弹性模量值,γ为硬质合金的切应变,bv为泊松比。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种整体式平底立铣刀的刀具失效界限图建立方法,根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型,筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型及对应的第一切削参数,根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型,筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型及对应的第二切削参数,根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型,筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型及对应的第三切削参数,在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线,将后刀面失效曲线、前刀面失效曲线和切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。本发明通过结合三种失效形式不同的失效标准,建立刀具失效界限图,为刀具在铣削钛合金时,提供一个安全区域,避免刀具过快失效,保证安全区域内的切削用量组合不会在预定更换刀具时间之前发生故障。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种整体式平底立铣刀的刀具失效界限图建立方法的流程图;
图2为本发明提供的整体式平底立铣刀的后刀面失效界限原理图;
图3为本发明提供的整体式平底立铣刀的后刀面磨损带正应力示意图;
图4为本发明提供的基于几何的刀具后刀面磨损体积示意图;
图5为本发明提供的整体式平底立铣刀的前刀面月牙洼失效界限原理图;
图6为本发明提供的前刀面月牙洼磨损建模示意图;
图7为本发明提供的切削刃破损失效原理图;
图8为本发明提供的具有安全切削区域的刀具失效界面图的示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种整体式平底立铣刀的刀具失效界限图建立方法及系统,为切削刀具提供了一个安全的操作区域,以保证安全区域内的切削用量组合不会在预定更换刀具时间之前发生故障。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明提供了一种整体式平底立铣刀的刀具失效界限图建立方法,如图1所示,方法包括:
S101,获取预设参数范围内整体式平底立铣刀的多个初始切削参数;每个初始切削参数包括切削厚度和切削速度,每个切削参数设置不同的切削厚度和切削速度;
S102,根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型;
S103,确定每个后刀面磨损宽度变化率模型的磨损宽度的最大值,并筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型,获得筛选出的后刀面磨损宽度变化率模型对应的初始切削参数,确定为第一切削参数;
S104,根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型;
S105,确定每个前刀面月牙洼磨损深度模型的磨损深度的最大值,并筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型,获得筛选出的前刀面月牙洼磨损深度模型对应的初始切削参数,确定为第二切削参数;
S106,根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型;
S107,确定每个临界断裂时许用应力模型的临界断裂时许用应力,并筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型,获得筛选出的临界断裂时许用应力模型对应的初始切削参数,确定为第三切削参数;
S108,在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线;同一坐标系的横坐标为切削速度,纵坐标为切削宽度;
S109,将后刀面失效曲线、前刀面失效曲线和切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。
具体过程如下:
步骤S102,根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型,如图2所示,具体包括:
根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000141
Figure BDA0002882022360000142
建立基于几何的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000143
将基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型的体积等价于基于几何的整体式平底立铣刀的后刀面磨损体积模型的体积,获得每个初始切削参数的整体式平底立铣刀的后刀面磨损宽度变化率模型为
Figure BDA0002882022360000144
其中,
Figure BDA0002882022360000145
为第z个基于磨损机理的整体式平底立铣刀的后刀面磨损体积,σn-f(x")为后刀面磨损带位置x″处的正应力,w为切屑宽度,vz为第z个初始切削参数的切削速度,Δtz为第z个切削有效时间,
Figure BDA0002882022360000146
fz为每齿进给量,R为刀具半径,
Figure BDA0002882022360000147
为第z个初始切削参数的切削宽度,H0为硬质合金硬度,a为与温度有关的系数,K为切削刃上剪应力和工件上剪切流动应力的比值,
Figure BDA0002882022360000151
为第z个刀具后刀面和工件接触面的温度,Ht为铣刀表面硬度值,n为三体磨粒磨损常数,VB为铣刀后刀面磨损宽度,d为粘结磨损系数,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数,
Figure BDA0002882022360000152
为第z个基于几何的整体式平底立铣刀的后刀面磨损体积,λs为刀具刃倾角,
Figure BDA0002882022360000153
为第z个后刀面磨损宽度变化率,ΔVBz为第z个铣刀后刀面磨损宽度变化量,α为后角,γn为前角。
步骤S102中,根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型,具体包括:
建立整体式平底立铣刀的后刀面磨损带正应力分布模型为
Figure BDA0002882022360000154
根据每个初始切削参数的切削速度,建立整体式平底立铣刀的后刀面磨损带温度分布函数模型为
Figure BDA0002882022360000155
根据每个初始切削参数的切削厚度,计算一个刀齿切入到切出过程的切削有效时间为
Figure BDA0002882022360000156
根据后刀面磨损带正应力分布模型、后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于磨粒磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面磨粒磨损体积模型为
Figure BDA0002882022360000157
根据后刀面磨损带正应力分布模型、后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于粘结磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面粘结磨损体积模型为
Figure BDA0002882022360000158
根据后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于扩散磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面扩散磨损体积模型为
Figure BDA0002882022360000161
根据每个切削参数建立的后刀面磨粒磨损体积模型、后刀面粘结磨损体积模型和后刀面扩散磨损体积模型,分别建立每个初始切削参数的基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000162
Figure BDA0002882022360000163
其中,σtip-flank为微元刀尖上靠近后刀面磨损带的正应力,
Figure BDA0002882022360000164
ρ是材料密度,φn为剪切角,ηw为后刀面磨损带上摩擦力的滑移线场角,VBCR为弹性接触区域的宽度,μc为后刀面摩擦系数,k是材料热导率,c是材料比热容,i为第一累加运算符号,j为第二累加运算符号,
Figure BDA0002882022360000165
Figure BDA0002882022360000166
均为排列组合运算,T0为室温,
Figure BDA0002882022360000167
为泰勒展开式系数,t为切削时间,
Figure BDA0002882022360000168
为第z个后刀面磨粒磨损体积,Cflank-abrasion为无量纲磨粒磨损系数,
Figure BDA0002882022360000169
为第z个粘结磨损导致的铣刀后刀面磨损体积,Cflank-adhension为无量纲粘结磨损系数,
Figure BDA00028820223600001610
为第z个扩散磨损机理导致的铣刀后刀面磨损体积,Cflank-diffusion为无量纲扩散磨损系数,
Figure BDA00028820223600001611
为第z个刀具后刀面总磨损体积。
步骤S103,磨损宽度阈值为0.3mm。
步骤S104,根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型,如图5所示,具体包括:
根据每个初始切削参数的切削厚度,分别建立整体式平底立铣刀的前刀面法向接触应力分布模型为
Figure BDA00028820223600001612
根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面温度分布函数模型为
Figure BDA0002882022360000171
根据每个初始切削参数的切削速度,基于磨粒磨损机理,分别建立整体式平底立铣刀的前刀面磨粒磨损率模型为
Figure BDA0002882022360000172
根据每个前刀面法向接触应力分布模型和每个前刀面温度分布函数模型,基于粘结磨损机理,分别建立整体式平底立铣刀的前刀面粘结磨损率模型为
Figure BDA0002882022360000173
根据每个前刀面温度分布函数模型和每个前刀面温度分布函数模型对应的切削速度,基于扩散磨损机理,分别建立整体式平底立铣刀的前刀面扩散磨损率模型为
Figure BDA0002882022360000174
根据每个切削参数建立的前刀面磨粒磨损率模型、前刀面粘结磨损率模型和前刀面扩散磨损率模型,分别建立每个初始切削参数的整体式平底立铣刀的前刀面月牙洼深度变化率模型为
Figure BDA0002882022360000175
其中,
Figure BDA0002882022360000176
为第z个初始切削参数的切削厚度对应的前刀面法向接触应力,σtip-rake是刀尖处最大法向接触应力,
Figure BDA0002882022360000177
是第z个初始切削参数的切削厚度对应的刀屑接触长度,
Figure BDA0002882022360000178
Figure BDA0002882022360000179
为第z个初始切削参数的切削厚度对应的最大切屑厚度,
Figure BDA00028820223600001710
φn为剪切角,γn为前角,βn为摩擦角,ξ1为相关幂指数,ηc为切屑流出角,R为刀具半径,f为每齿进给量,ξ为应力指数参数,x为刀屑接触距离,
Figure BDA00028820223600001711
为第z个初始切削参数对应的前刀面温度,
Figure BDA00028820223600001712
是刀屑平均摩擦系数,
Figure BDA00028820223600001713
为第z个初始切削参数的切削速度对应的切屑速度,
Figure BDA00028820223600001714
λs为刀具刃倾角,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容,T0为室温,
Figure BDA00028820223600001715
为第z个磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率,Qtool-chip为刀屑接触面之间的接触载荷,Hrake-face为刀具前刀面的硬度,θabrasion为磨粒平均粗糙角,
Figure BDA0002882022360000181
为第z个粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率,N1为第一常数,N2为第二常数,
Figure BDA0002882022360000182
为第z个扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率,Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数,KT(x,t)为第z个初始切削参数的前刀面月牙洼深度变化率,t为切削时间。
步骤S105,磨损深度阈值为0.1mm。
步骤S106,根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型,具体包括:
根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型为
Figure BDA0002882022360000183
其中,
Figure BDA0002882022360000184
为第z个临界断裂时许用应力,Fz为第z个初始切削参数的切削速度对应的导致刀齿破损的冲击力,
Figure BDA0002882022360000185
fz为每齿进给量,ρ为材料密度,vz为切削速度,ap是切削深度,mA为断裂面积相关比系数,VB为铣刀后刀面磨损宽度,τ为许用应力,
Figure BDA0002882022360000186
为第z个初始切削参数的切削厚度对应的损伤弹性模量值,γ为硬质合金的切应变,bv为泊松比。
步骤S109,将后刀面失效曲线、前刀面失效曲线和切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,具体包括:
确定切削刃破损失效曲线上最小切削速度所在的第一坐标点,建立经过第一坐标点且平行于纵坐标轴的第一直线;
确定后刀面失效曲线上最大切削速度所在的第二坐标点,建立经过第二坐标点且平行于横坐标轴的第二直线;
将第一直线、第二直线、后刀面失效曲线、前刀面失效曲线和切削刃破损失效曲线形成的封闭区域确定为安全切削区域。
本发明通过分析热-力耦合作用下刀具损伤机理及失效行为,建立了刀具损伤解析模型,利用Matlab软件预测高速干式切削刀具失效界限,建立刀具失效界限图,为刀具在铣削钛合金时,提供一个安全区域,避免刀具过快失效。刀具失效界限图是通过叠加不同类型的失效边界为切削刀具提供了一个安全的操作区域,以保证安全区域内的切削用量组合不会在预定更换刀具时间之前发生故障,为硬质合金刀具铣削钛合金加工领域提供借鉴意义。
本发明针对一种整体式平底立铣刀的刀具失效界限图建立方法,还提供了一个具体实施例。
建立整体式平底立铣刀后刀面磨损宽度VB变化率模型,具体包括:
建立由于磨粒磨损机理导致的整体式平底立铣刀后刀面磨粒磨损体积模型具体包括:
建立三体磨粒磨损体积模型:
Figure BDA0002882022360000191
其中,Wthree-body为三体磨粒磨损体积,Nflank为刀具后刀面磨损带正压力,L是磨粒滑动距离,θ是磨粒平均粗糙度角,Ht为铣刀表面硬度值,Ha为磨粒硬度值,
Figure BDA0002882022360000192
为磨粒硬度值的n-1次方,
Figure BDA0002882022360000193
为刀具硬度值的n次方,n、K为三体磨粒磨损常数,n与K则有:
Figure BDA0002882022360000194
建立刀具后刀面磨损带正压力Nflank模型:
Nflank=σflankAslide (3)
其中,σflank为刀具后刀面磨损带正应力,Aslide为磨粒在刀具后刀面磨损带上滑动过的面积。
刀具后刀面磨损带上的正压力Nflank与作用在所有磨粒上的正压力存在关系。
Nflankpflank-abrasion-particles%=Fflank-particlenflank-particlesAslide (4)
其中,pflank-abration-particles%为磨粒上的正压力占平底立铣刀后刀面磨损带正压力的比例,Fflank-particle为作用在所有磨粒上的正压力,nflank-particles为单位接触面上的磨粒数。
联立公式(3)和公式(4),可得到作用在平底立铣刀后刀面磨损带上每个磨粒上的正压力Fflank-particle,如式(4)所示。
Figure BDA0002882022360000201
计算一段加工时间Δt之后,磨粒在刀具后刀面磨损带上滑动过的面积Aslide
Figure BDA0002882022360000202
其中,w为切屑宽度,w'为斜角切削时的切屑厚度,λs为刀具刃倾角。
进一步得到所形成的滑动面积上的磨粒总数Nflank-particles
Nflank-particles=nflank-particlesAslide=nflank-particleswvΔt (7)
联立公式(1)至公式(7),可得由于磨粒磨损机理导致的铣刀后刀面磨损体积Vflank-abrasion-wear
Figure BDA0002882022360000203
进一步对Vflank-abrasion-wear简化变形。
Figure BDA0002882022360000204
其中,Cflank-abrasion=0.5Pflank-abrasion-particles%tanθ是一个无量纲磨粒磨损系数。
建立由于粘结磨损机理导致的整体式平底立铣刀后刀面粘结磨损体积模型具体包括:
计算加工工件的硬度Hw
Figure BDA0002882022360000211
其中,VB为铣刀后刀面磨损宽度,nflank为微焊点个数,
Figure BDA0002882022360000212
为微焊点横截面的平均面积。
进一步变形:
Figure BDA0002882022360000213
计算微焊点的体积Vflank-m
Figure BDA0002882022360000214
其中,hflank-m为微焊点的高度。
刀具材料硬度Ht与微焊点的高度hflank-m存在关系:
Figure BDA0002882022360000215
其中,K0为系数。
联立公式(10)-(13),可得刀由于粘结磨损导致铣刀后刀面平均磨损量
Figure BDA0002882022360000216
Figure BDA0002882022360000217
计算经过一段铣削加工时间Δt之后,由于粘结磨损导致铣刀后刀面磨损量Vflank-adhension-wear
Figure BDA0002882022360000218
其中,Anew为粘结磨损面积,Anew=wvΔt,p0为常数。
刀具后刀面和工件接触面的硬度取决于接触面的温度、应变以及应变率等因素。为简化计算,本文只考虑接触面温度因素。
因此,建立关于温度的指数函数模型为:
Figure BDA0002882022360000219
其中,c为粘结磨损系数,d为粘结磨损系数。
进一步对Vflank-adhension-wear变形:
Figure BDA0002882022360000221
其中,Cflank-adhension=c-1p0K0(1-pflank-adration-particles%)是一个无量纲粘结磨损系数。
建立由于扩散磨损机理导致的整体式平底立铣刀后刀面扩散磨损体积模型具体包括:
计算沿后刀面和工件接触面上任意位置z″处的浓度扩散梯度:
Figure BDA0002882022360000222
图3中,X″和Z″建立后刀面和工件的接触面,X″为横坐标轴,Z″为纵坐标轴,x″为横坐标轴X″上一点,表示铣刀后刀面磨损宽度。
其中,元素扩散系数D与温度Tflank(VB)关系为:
Figure BDA0002882022360000223
其中,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数。
联立公式(18)-(19),计算沿后刀面和工件接触面上任意位置处元素扩散通量Jy=0
Figure BDA0002882022360000224
计算后刀面和工件接触面上,新形成的接触点以速度v在其上滑动,所需要的滑动时间ts
Figure BDA0002882022360000225
计算经过一段铣削加工时间Δt之后,由于扩散磨损机理导致的平底立铣刀后刀面磨损体积Vflank-diffusion-wear
Figure BDA0002882022360000231
其中,mCo为Co元素的分子质量,PCo为Co元素的硬度,ρt为时间t内的扩散密度。
进一步变形:
Figure BDA0002882022360000232
其中,Cflank-diffusion=-2C0·mCo/PCoρt是一个无量纲扩散磨损系数。
建立整体式平底立铣刀后刀面磨损影响的正应力分布模型:
Figure BDA0002882022360000233
其中,VBCR为塑性流动区域和弹性接触区域临界点,此处亦为弹性接触区域的宽度,σtip-flank为微元刀尖上靠近后刀面磨损带的正应力,x″为图3中横坐标轴X″上的一点。
计算微元刀尖上靠近后刀面磨损带的正应力σtip-flank
Figure BDA0002882022360000234
其中,ηw为后刀面磨损带上摩擦力的滑移线场角,ηw=0.5cos-1(mw),φn为剪切角,K为切削刃上剪应力和工件上剪切流动应力的比值,它在数值上等于靠近后刀面的切削刃摩擦系数(μc),由于后刀面磨损带与已加工表面接触的粘结本质,认为μc是均匀的,mw为后刀面摩擦系数,它的数值大小与靠近后刀面的切削刃摩擦系数μc相等。
建立整体式平底立铣刀后刀面温度分布函数模型:
Figure BDA0002882022360000241
其中:ρ是材料密度,k是材料热导率,c是材料比热容,μc为后刀面摩擦系数,
Figure BDA0002882022360000242
为泰勒展开式系数,等于2;i和j均为累加运算符号,
Figure BDA0002882022360000243
Figure BDA0002882022360000244
均为排列组合运算,T0为室温。
计算一个刀齿切入到切出过程的有效时间:
铣削过程中,刀齿切入角、切出角与切削有效时间关系为:
Figure BDA0002882022360000245
其中,θin为切入角,θout为切出角,Δt为切削有效时间,fz为每齿进给量,ae为切削宽度,t为切削时间。
建立基于磨损机理的整体式平底立铣刀后刀面磨损体积模型:
假定上述建立的三种磨损模型互不影响,因此刀具后刀面总磨损体积Vflank-wear-mechanism为三种磨损模型的叠加。
Figure BDA0002882022360000246
其中,
Figure BDA0002882022360000247
H0为硬质合金硬度,a为与温度有关的系数。
建立基于几何的整体式平底立铣刀后刀面磨损体积模型:
Figure BDA0002882022360000251
如图4所示,γn为前角,α为后角,w为切屑厚度,w′为斜角切削时的切屑厚度,A、B、C、D、E、F、G、H分别表示构成后刀面磨损侧面的各个顶点,FED所在平面为后刀面磨损平面,ABC所在平面为t时间后的后刀面磨损平面。
为简化式(29),此处忽略高次项|AE|2,进一步变形:
Figure BDA0002882022360000252
建立整体式平底立铣刀后刀面磨损宽度变化率模型:
将建立的基于机理的整体式平底立铣刀后刀面磨损体积Vflank-wear-mechanism等价于基于几何的整体式平底立铣刀后刀面磨损体积Vflank-wear-geometry:
Figure BDA0002882022360000253
进一步,变形可得整体式平底立铣刀后刀面磨损宽度VB随时间的变化率:
Figure BDA0002882022360000254
建立整体式平底立铣刀前刀面月牙洼磨损模型,具体包括:
建立整体式平底立铣刀前刀面法向接触应力分布模型:
σrake(x)=σtip-rake(1-x/ltool-chip)ξ (33)
其中:σtip-rake是刀尖处最大法向接触应力,ltool-chip是刀屑接触长度,ξ是应力指数参数,x为图6中横坐标轴X上一点,表示刀屑接触距离。
计算刀具前刀面刀尖处最大正应力σtip-rake
Figure BDA0002882022360000261
其中:φn是剪切角,γn是前角,βn是摩擦角,τh为剪切区出口处的剪应力,b是相关指数,v为切削速度。
考虑在最大切屑厚度hmax位置处,刀具承受最大的应力和温度,此位置刀具磨损最严重,计算刀屑接触长度ltool-chip
Figure BDA0002882022360000262
其中:ξ1是相关幂指数,ηc是切屑流出角。
计算最大切屑厚度hmax
Figure BDA0002882022360000263
其中:R是刀具半径,f是每齿进给量,ae是切削厚度。
建立与切屑速度、刀屑接触长度有关的前刀面温度分布函数模型:
Figure BDA0002882022360000264
其中:
Figure BDA0002882022360000265
是刀屑平均摩擦系数,vc为切屑速度,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容。
计算流过前刀面的切屑速度:
Figure BDA0002882022360000266
建立由于磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率模型:
计算刀屑接触面之间的接触载荷为Qtool-chip
Figure BDA0002882022360000271
其中,Hrake-face为刀具前刀面的硬度,Sabrasion为磨粒与前刀面接触表面面积,dabrasion为磨粒与前刀面接触划过凹坑表面直径。
计算磨粒在前刀面磨损带滑动L距离时,划过的单位磨粒磨损体积:
Figure BDA0002882022360000272
其中:θabrasion为磨粒平均粗糙角,vc为相对滑移速度。
进一步,可得于磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率模型:
Figure BDA0002882022360000273
建立由于粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率模型:
铣刀前刀面粘结磨损体积dVrake-adhesion-wear在滑移距离dl为:
Figure BDA0002882022360000274
其中:c、b分别为切屑上硬质颗粒的高度和两相邻硬质颗粒的平均间距,由于尺寸效应的存在,c/b可以近似看作是常数,Z为接触界面上硬质颗粒产生粘结磨损的Holm概率,Hrake-particle为粗糙颗粒的硬度。
由于粗糙颗粒硬度h与刀-屑接触表面对的整体性能有关,即它主要取决于接触表面上的应变、应变率以及温度,又因为应变和应变率也受温度的影响,忽略硬质合金刀具实际切削条件范围内的应变和应变速率的变化所以h可以用刀-屑接触表面上的温度来表示为:
Figure BDA0002882022360000275
其中:A1、A2为常数,T为接触表面温度。
由于粘结磨损的发生需要依赖温度激活,所以Holm概率Z表示为:
Figure BDA0002882022360000276
其中:B1为常数,E为激活能,λ为玻尔兹曼常数。
联立以上三式可得:
Figure BDA0002882022360000277
其中:N1为常数,(E+λA2)/λ主要取决于接触表面的扩散层结构和元素浓度,如果切削条件的变化有限,我们可以认为它是近似恒定的,记为N2
进一步,将N2带入式(50),可得由于粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率模型:
Figure BDA0002882022360000281
建立由于扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率模型:
沿刀屑接触面上任意位置处的C0元素浓度扩散梯度为:
Figure BDA0002882022360000282
计算扩散系数Drake是与前刀面温度Trake(x)有关的关系式:
Figure BDA0002882022360000283
根据FICK扩散定律,计算前刀面与刀屑接触界面任意一点沿y方向的扩散通量:
Figure BDA0002882022360000284
进一步,可得由于扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率模型:
Figure BDA0002882022360000285
其中:Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数。
建立整体式立铣刀前刀面月牙洼深度模型:
刀具前刀面磨损是多种磨损机理共同作用的结果,只考虑一种磨损机理并不能反映刀具磨损的真实情况。因此综合考虑三种磨损机理,假定他们彼此无关联,且不考虑切削过程中的非正常磨损,则前刀面复合磨损率模型为:
Figure BDA0002882022360000291
在某一切削时间t内的单位切削厚度的月牙洼磨损量Vrake为该点磨损率与单位距离dx的乘积,如图6所示,则有:
Figure BDA0002882022360000292
进一步,得到与刀屑接触界面温度,切削时间相关的前刀面不同位置处月牙洼磨损深度仿真模型:
Figure BDA0002882022360000293
建立所述整体式平底立铣刀临界断裂时许用应力模型:
如图7所示,建立作用在切削刃上的等效单轴的压应力模型:
采用损伤当量应力这一概念,可将复杂的三轴应力状态简化为单轴应力状态,并保证在同一损伤值下,这两种应力状态可以产生同等的弹性形变能应变能密度。这里假定原三轴应力状态以及转化后的单轴应力状态对刀具材料的损伤演变作用是相同的。将裂纹闭合的准单边条件考虑进去,则单轴应力模型为:
Figure BDA0002882022360000294
其中,σc是损伤当量应力,σ是三轴状态下的应力,bv是刀具材料的泊松比,D是刀具材料的损伤值,h是相关系数可以取为0.2。
假设刀具材料单元预先存在N个初始的微裂纹。那么,如果不考虑微裂纹中间相互的作用,总的损伤应变就可以表达,选定的刀具材料单元的的总的应变可以表达为:
Figure BDA0002882022360000295
进一步,令ε1=ε,对上式化简就可以得到:
Figure BDA0002882022360000296
其中,S11、S12、S21和S22分别为无量纲矩阵的元素,S11、S12、S21和S22均为常量,且S12=S21;ε1和ε2为刀具材料单元的两个方向的应变,
Figure BDA0002882022360000301
Figure BDA0002882022360000302
为无损伤刀具材料的弹性应变,Δε1和Δε2为由初始微裂纹滑动以及滑动裂纹扩展所导致的总的损伤应变,k由泊松比决定参数,对于平面应变条件k=3-4bv,而对于平面应力条件k=[(3-bv)(1+bv)],bv为泊松比,E为弹性模量。
计算损伤弹性模量值ED
在单轴的压应力下,损伤材料的应变是由弹性模量ED控制的,并且弹性模量这作为刀具的固有的参数,随着且切削不断地进行,刀具材料不断地进行损伤累积,因此随着用刀具材料退化的损伤模量来定义材料单元的损伤:
Figure BDA0002882022360000303
其中,E是弹性模量,N是初始的微裂纹数。
对刀具材料的损伤值与退化的弹性模量进行求解,可以分别求解张开裂纹扩展所导致的弹性应变能We,初始微裂纹滑动所释放的摩擦能Wf,施加的单轴压应力所做的功W1,通过建立处于压应力作用下还有单个滑动裂纹体的能量平衡方程进行求解:
W1作为单轴压应力所做的功,可以表达为
Figure BDA0002882022360000304
式中,σ1和σ2为平面内两个方向的正应力,a=b是选定单元的长度与裂纹长度为同一数量级,单位是μm。
而进一步的,根据Li和Tada等人的研究,裂纹阵列的应力强度因子可以通过以下的式子进行计算
Figure BDA0002882022360000305
Figure BDA0002882022360000306
其中,l1为第1次冲击裂纹长度,l1=0.27c,l1的作用是对较小裂纹长度进行修正,2c是初始裂纹长度,2c=d,其中d是刀具材料的粒径。ln是第n次冲击裂纹长度,初始裂纹间距2w是4倍的裂纹长度就是在裂纹扩展初期也能是应力强度因子的计算结果保持准确。
目前在工程中应用最为广泛的方法依然是1963年由Paris和Erdogan在实验基础上提出的疲劳裂纹扩展公式,这就是著名的Paris公式,它建立了应力强度因子和裂纹扩展速率之间的关系,是当今工程应用中预测疲劳裂纹扩展寿命理论的基础,其形式为:
Figure BDA0002882022360000311
其中,n是应力循环次数也就是刀具的切削次数,C、m都是材料常数,可以根据实验数据拟合得到。
根据裂纹长度随着载荷的循环次数的变化规律,裂纹扩展速率dl/dn与应力强度因子ΔK之间的关系可以根据Paris公式进行描述,通过三点弯曲试验可以得到硬质合金材料的裂纹扩展速率可以表示为:
Figure BDA0002882022360000312
Figure BDA0002882022360000313
Wf=2cτfχ (64)
其中,We是裂纹扩展的弹性应变能,Wf是初始滑动裂纹的摩擦能,τf为单轴压应力引起的应力法向分量的裂纹表而切向分量,χ为初始微裂纹表面的滑动距离,k为相关系数。
Figure BDA0002882022360000314
Figure BDA0002882022360000315
式中,dln/dn表示刀具材料承受第n次冲击时的裂纹扩展速度。
Figure BDA0002882022360000316
和l0分别为刀具材料的初始损伤值和初始裂纹扩展速率。在刀具材料承受n次冲击后,刀具材料损伤值为Dn,张开裂纹扩展长度为ln。确定损伤值Dn和张开裂纹扩展长度ln是刀具材料承受(n+1)次冲击时的初始条件。在(n+1)次冲击前,张开裂纹扩展速度为0。
2We+Wf=W1 (67)
其中,We为张开裂纹扩展所导致的弹性应变能,Wf是有初始微裂纹滑动所释放的摩擦能,W1是施加的单轴压应力所做的功。
将式(63)、(64)带入到式(67)中即可得到式(57)中的S11,并进一步获得刀具材料的损伤值和损伤后的弹性模量。
在进行刀具冲击破损计算时,铣刀作为断续切削导致刀齿破损的冲击力F:
F=fzapv2ρ (68)
其中,L为冲击行程,fz为每齿进给量,ρ为材料密度,v为切削速度。
建立刀具损伤以后的冲击断裂失效标准:
假设冲击过程中没有能量损失,即冲击过程中动能全部转化成用于刀齿断裂时能量,那么硬质合金刀具冲击断裂临界条件为:
F≥τA=GγA (69)
其中,G为硬质合金刀具的切变模量,γ为硬质合金的切应变,A为断裂时横截面积。
进一步,根据剪切胡克定律可知:
Figure BDA0002882022360000321
其中,ED为弹性模量,bv为泊松比。
失效标准为刀具临界断裂时应力为ξB
Figure BDA0002882022360000322
其中,ap是切削深度,mA是由实验获得断裂面积相关比系数。
建立所述整体式立铣刀失效界限图具体包括:
刀具失效界限图是通过叠加不同类型的失效边界为切削刀具提供了一个安全的操作区域,以保证安全区域内的切削用量组合不会在预定更换刀具时间之前发生故障。本阶段通过分析热-力耦合作用下刀具损伤机理及失效行为,建立了刀具损伤解析模型,利用Matlab软件预测高速干式切削刀具失效界限。
基于钛合金的加工特性,失效界限图主要考虑了三种刀具失效的原因:后刀面磨损、月牙洼磨损和切削刃破损。
(1)后刀面磨损失效标准为后刀面磨损量VBmax为0.3mm,建立后刀面失效界限。
(2)前刀面月牙洼磨损失效标准为月牙洼深度KTmax为0.1mm,建立前刀面月牙洼失效界限。
(3)切削刃破损,失效标准为刀具临界断裂时许用应力为ξB时,刀具因破损失效,建立切削刃破损失效界限。
目前铣削Ti6Al4V的刀具寿命至少应为30分钟,因此当刀具在30分钟以内达到某一失效形式的失效标准时,认为该刀具提前失效。如图8所示,横坐标为热载荷用切削速度表征,纵坐标为机械载荷用切削宽度表征的二维,根据三种失效形式不同的失效标准,分别获得三条失效界限曲线,以三种失效界限交点为中心,进行合理区间取值,建立了在不同切削参数下的刀具失效界限图。
本发明还提供了一种整体式平底立铣刀的刀具失效界限图建立系统,系统包括:
初始切削参数获取模块,用于获取预设参数范围内整体式平底立铣刀的多个初始切削参数;每个初始切削参数包括切削厚度和切削速度,每个切削参数设置不同的切削厚度和切削速度;
后刀面磨损宽度变化率模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型;
第一切削参数确定模块,用于确定每个后刀面磨损宽度变化率模型的磨损宽度的最大值,并筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型,获得筛选出的后刀面磨损宽度变化率模型对应的初始切削参数,确定为第一切削参数;
前刀面月牙洼磨损深度变化率模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型;
第二切削参数确定模块,用于确定每个前刀面月牙洼磨损深度模型的磨损深度的最大值,并筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型,获得筛选出的前刀面月牙洼磨损深度模型对应的初始切削参数,确定为第二切削参数;
临界断裂时许用应力模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型;
第三切削参数确定模块,用于确定每个临界断裂时许用应力模型的临界断裂时许用应力,并筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型,获得筛选出的临界断裂时许用应力模型对应的初始切削参数,确定为第三切削参数;
失效曲线建立模块,用于在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线;同一坐标系的横坐标为切削速度,纵坐标为切削宽度;
刀具失效界面图获得模块,用于将后刀面失效曲线、前刀面失效曲线和切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。
后刀面磨损宽度变化率模型建立模块,具体包括:
基于磨损机理的后刀面磨损体积模型建立子模块,用于根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000341
Figure BDA0002882022360000342
基于几何的后刀面磨损体积模型建立子模块,用于建立基于几何的整体式平底立铣刀的后刀面磨损体积模型为
Figure BDA0002882022360000343
后刀面磨损宽度变化率模型获得子模块,用于将基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型的体积等价于基于几何的整体式平底立铣刀的后刀面磨损体积模型的体积,获得每个初始切削参数的整体式平底立铣刀的后刀面磨损宽度变化率模型为
Figure BDA0002882022360000344
其中,
Figure BDA0002882022360000345
为第z个基于磨损机理的整体式平底立铣刀的后刀面磨损体积,σn-f(x")为后刀面磨损带位置x″处的正应力,w为切屑宽度,vz为第z个初始切削参数的切削速度,Δtz为第z个切削有效时间,
Figure BDA0002882022360000351
fz为每齿进给量,R为刀具半径,
Figure BDA0002882022360000352
为第z个初始切削参数的切削宽度,H0为硬质合金硬度,a为与温度有关的系数,K为切削刃上剪应力和工件上剪切流动应力的比值,
Figure BDA0002882022360000353
为第z个刀具后刀面和工件接触面的温度,Ht为铣刀表面硬度值,n为三体磨粒磨损常数,VB为铣刀后刀面磨损宽度,d为粘结磨损系数,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数,
Figure BDA0002882022360000354
为第z个基于几何的整体式平底立铣刀的后刀面磨损体积,λs为刀具刃倾角,
Figure BDA0002882022360000355
为第z个后刀面磨损宽度变化率,ΔVBz为第z个铣刀后刀面磨损宽度变化量,α为后角,γn为前角。
前刀面月牙洼磨损深度变化率模型建立模块,具体包括:
前刀面法向接触应力分布模型建立子模块,用于根据每个初始切削参数的切削厚度,分别建立整体式平底立铣刀的前刀面法向接触应力分布模型为
Figure BDA0002882022360000356
前刀面温度分布函数模型建立子模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面温度分布函数模型为
Figure BDA0002882022360000357
前刀面磨粒磨损率模型建立子模块,用于根据每个初始切削参数的切削速度,基于磨粒磨损机理,分别建立整体式平底立铣刀的前刀面磨粒磨损率模型为
Figure BDA0002882022360000358
前刀面粘结磨损率模型建立子模块,用于根据每个前刀面法向接触应力分布模型和每个前刀面温度分布函数模型,基于粘结磨损机理,分别建立整体式平底立铣刀的前刀面粘结磨损率模型为
Figure BDA0002882022360000359
前刀面扩散磨损率模型建立子模块,用于根据每个前刀面温度分布函数模型和每个前刀面温度分布函数模型对应的切削速度,基于扩散磨损机理,分别建立整体式平底立铣刀的前刀面扩散磨损率模型为
Figure BDA0002882022360000361
前刀面月牙洼深度变化率模型建立子模块,用于根据每个切削参数建立的前刀面磨粒磨损率模型、前刀面粘结磨损率模型和前刀面扩散磨损率模型,分别建立每个初始切削参数的整体式平底立铣刀的前刀面月牙洼深度变化率模型为
Figure BDA0002882022360000362
其中,
Figure BDA0002882022360000363
为第z个初始切削参数的切削厚度对应的前刀面法向接触应力,σtip-rake是刀尖处最大法向接触应力,
Figure BDA0002882022360000364
是第z个初始切削参数的切削厚度对应的刀屑接触长度,
Figure BDA0002882022360000365
Figure BDA0002882022360000366
为第z个初始切削参数的切削厚度对应的最大切屑厚度,
Figure BDA0002882022360000367
φn为剪切角,γn为前角,βn为摩擦角,ξ1为相关幂指数,ηc为切屑流出角,R为刀具半径,f为每齿进给量,ξ为应力指数参数,x为刀屑接触距离,
Figure BDA0002882022360000368
为第z个初始切削参数对应的前刀面温度,
Figure BDA0002882022360000369
是刀屑平均摩擦系数,
Figure BDA00028820223600003610
为第z个初始切削参数的切削速度对应的切屑速度,
Figure BDA00028820223600003611
λs为刀具刃倾角,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容,T0为室温,
Figure BDA00028820223600003612
为第z个磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率,Qtool-chip为刀屑接触面之间的接触载荷,Hrake-face为刀具前刀面的硬度,θabrasion为磨粒平均粗糙角,
Figure BDA00028820223600003613
为第z个粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率,N1为第一常数,N2为第二常数,
Figure BDA00028820223600003614
为第z个扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率,Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数,KT(x,t)为第z个初始切削参数的前刀面月牙洼深度变化率,t为切削时间。
临界断裂时许用应力模型建立模块,具体包括:
临界断裂时许用应力模型建立子模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型为
Figure BDA0002882022360000371
其中,
Figure BDA0002882022360000372
为第z个临界断裂时许用应力,Fz为第z个初始切削参数的切削速度对应的导致刀齿破损的冲击力,
Figure BDA0002882022360000373
fz为每齿进给量,ρ为材料密度,vz为切削速度,ap是切削深度,mA为断裂面积相关比系数,VB为铣刀后刀面磨损宽度,τ为许用应力,
Figure BDA0002882022360000374
为第z个初始切削参数的切削厚度对应的损伤弹性模量值,γ为硬质合金的切应变,bv为泊松比。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (4)

1.一种整体式平底立铣刀的刀具失效界限图建立方法,其特征在于,所述方法包括:
获取预设参数范围内整体式平底立铣刀的多个初始切削参数;每个初始切削参数包括切削厚度和切削速度,每个切削参数设置不同的切削厚度和切削速度;
根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型,具体包括:根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure FDA0003810554180000011
建立基于几何的整体式平底立铣刀的后刀面磨损体积模型为
Figure FDA0003810554180000012
将基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型的体积等价于基于几何的整体式平底立铣刀的后刀面磨损体积模型的体积,获得每个初始切削参数的整体式平底立铣刀的后刀面磨损宽度变化率模型为
Figure FDA0003810554180000013
其中,
Figure FDA0003810554180000014
为第z个基于磨损机理的整体式平底立铣刀的后刀面磨损体积,σn-f(x")为后刀面磨损带位置x″处的正应力,w为切屑宽度,vz为第z个初始切削参数的切削速度,Δtz为第z个切削有效时间,
Figure FDA0003810554180000015
fz为每齿进给量,R为刀具半径,
Figure FDA0003810554180000016
为第z个初始切削参数的切削厚度,H0为硬质合金硬度,a为与温度有关的系数,K为切削刃上剪应力和工件上剪切流动应力的比值,
Figure FDA0003810554180000017
为第z个刀具后刀面和工件接触面的温度,Ht为铣刀表面硬度值,n为三体磨粒磨损常数,VB为铣刀后刀面磨损宽度,d为粘结磨损系数,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数,
Figure FDA0003810554180000021
为第z个基于几何的整体式平底立铣刀的后刀面磨损体积,λs为刀具刃倾角,
Figure FDA0003810554180000022
为第z个后刀面磨损宽度变化率,ΔVBz为第z个铣刀后刀面磨损宽度变化量,α为后角,γn为前角;
确定每个后刀面磨损宽度变化率模型的磨损宽度的最大值,并筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型,获得筛选出的后刀面磨损宽度变化率模型对应的初始切削参数,确定为第一切削参数;
根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型,具体包括:根据每个初始切削参数的切削厚度,分别建立整体式平底立铣刀的前刀面法向接触应力分布模型为
Figure FDA0003810554180000023
根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面温度分布函数模型为
Figure FDA0003810554180000024
根据每个初始切削参数的切削速度,基于磨粒磨损机理,分别建立整体式平底立铣刀的前刀面磨粒磨损率模型为
Figure FDA0003810554180000025
根据每个前刀面法向接触应力分布模型和每个前刀面温度分布函数模型,基于粘结磨损机理,分别建立整体式平底立铣刀的前刀面粘结磨损率模型为
Figure FDA0003810554180000026
根据每个前刀面温度分布函数模型和每个前刀面温度分布函数模型对应的切削速度,基于扩散磨损机理,分别建立整体式平底立铣刀的前刀面扩散磨损率模型为
Figure FDA0003810554180000027
根据每个切削参数建立的前刀面磨粒磨损率模型、前刀面粘结磨损率模型和前刀面扩散磨损率模型,分别建立每个初始切削参数的整体式平底立铣刀的前刀面月牙洼深度变化率模型为
Figure FDA0003810554180000031
其中,
Figure FDA0003810554180000032
为第z个初始切削参数的切削厚度对应的前刀面法向接触应力,σtip-rake是刀尖处最大法向接触应力,
Figure FDA0003810554180000033
是第z个初始切削参数的切削厚度对应的刀屑接触长度,
Figure FDA0003810554180000034
Figure FDA0003810554180000035
为第z个初始切削参数的切削厚度对应的最大切屑厚度,
Figure FDA0003810554180000036
φn为剪切角,γn为前角,βn为摩擦角,ξ1为相关幂指数,ηc为切屑流出角,R为刀具半径,f为每齿进给量,ξ为应力指数参数,x为刀屑接触距离,
Figure FDA0003810554180000037
为第z个初始切削参数对应的前刀面温度,
Figure FDA0003810554180000038
是刀屑平均摩擦系数,
Figure FDA0003810554180000039
为第z个初始切削参数的切削速度对应的切屑速度,
Figure FDA00038105541800000310
λs为刀具刃倾角,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容,T0为室温,
Figure FDA00038105541800000311
为第z个磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率,Qtool-chip为刀屑接触面之间的接触载荷,Hrake-face为刀具前刀面的硬度,θabrasion为磨粒平均粗糙角,
Figure FDA00038105541800000312
为第z个粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率,N1为第一常数,N2为第二常数,
Figure FDA00038105541800000313
为第z个扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率,Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数,KT(x,t)为第z个初始切削参数的前刀面月牙洼深度变化率,t为切削时间;
确定每个前刀面月牙洼磨损深度模型的磨损深度的最大值,并筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型,获得筛选出的前刀面月牙洼磨损深度模型对应的初始切削参数,确定为第二切削参数;
根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型,具体包括:根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型为
Figure FDA00038105541800000314
其中,
Figure FDA00038105541800000315
为第z个临界断裂时许用应力,Fz为第z个初始切削参数的切削速度对应的导致刀齿破损的冲击力,
Figure FDA00038105541800000316
fz为每齿进给量,ρ为材料密度,vz为切削速度,ap是切削深度,mA为断裂面积相关比系数,VB为铣刀后刀面磨损宽度,τ为许用应力,
Figure FDA0003810554180000041
为第z个初始切削参数的切削厚度对应的损伤弹性模量值,γ为硬质合金的切应变,bv为泊松比;
确定每个临界断裂时许用应力模型的临界断裂时许用应力,并筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型,获得筛选出的临界断裂时许用应力模型对应的初始切削参数,确定为第三切削参数;
在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线;所述同一坐标系的横坐标为切削速度,纵坐标为切削厚度;
将所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。
2.根据权利要求1所述的整体式平底立铣刀的刀具失效界限图建立方法,其特征在于,所述根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型,具体包括:
建立整体式平底立铣刀的后刀面磨损带正应力分布模型为
Figure FDA0003810554180000042
根据每个初始切削参数的切削速度,建立整体式平底立铣刀的后刀面磨损带温度分布函数模型为
Figure FDA0003810554180000043
根据每个初始切削参数的切削厚度,计算一个刀齿切入到切出过程的切削有效时间为
Figure FDA0003810554180000044
根据后刀面磨损带正应力分布模型、后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于磨粒磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面磨粒磨损体积模型为
Figure FDA0003810554180000051
根据后刀面磨损带正应力分布模型、后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于粘结磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面粘结磨损体积模型为
Figure FDA0003810554180000052
根据后刀面磨损带温度分布函数模型、每个初始切削参数的切削速度以及每个切削有效时间,基于扩散磨损机理,分别建立每个初始切削参数的整体式平底立铣刀的后刀面扩散磨损体积模型为
Figure FDA0003810554180000053
根据每个切削参数建立的后刀面磨粒磨损体积模型、后刀面粘结磨损体积模型和后刀面扩散磨损体积模型,分别建立每个初始切削参数的基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure FDA0003810554180000054
其中,σtip-flank为微元刀尖上靠近后刀面磨损带的正应力,
Figure FDA0003810554180000055
ρ是材料密度,φn为剪切角,ηw为后刀面磨损带上摩擦力的滑移线场角,VBCR为弹性接触区域的宽度,μc为后刀面摩擦系数,k是材料热导率,c是材料比热容,i为第一累加运算符号,j为第二累加运算符号,
Figure FDA0003810554180000056
Figure FDA0003810554180000057
均为排列组合运算,T0为室温,
Figure FDA0003810554180000058
为泰勒展开式系数,t为切削时间,
Figure FDA0003810554180000059
为第z个后刀面磨粒磨损体积,Cflank-abrasion为无量纲磨粒磨损系数,
Figure FDA00038105541800000510
为第z个粘结磨损导致的铣刀后刀面磨损体积,Cflank-adhension为无量纲粘结磨损系数,
Figure FDA00038105541800000511
为第z个扩散磨损机理导致的铣刀后刀面磨损体积,Cflank-diffusion为无量纲扩散磨损系数,
Figure FDA00038105541800000512
为第z个刀具后刀面总磨损体积。
3.根据权利要求1所述的整体式平底立铣刀的刀具失效界限图建立方法,其特征在于,所述将所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,具体包括:
确定所述切削刃破损失效曲线上最小切削速度所在的第一坐标点,建立经过所述第一坐标点且平行于纵坐标轴的第一直线;
确定所述后刀面失效曲线上最大切削速度所在的第二坐标点,建立经过所述第二坐标点且平行于横坐标轴的第二直线;
将所述第一直线、所述第二直线、所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线形成的封闭区域确定为安全切削区域。
4.一种整体式平底立铣刀的刀具失效界限图建立系统,其特征在于,所述系统包括:
初始切削参数获取模块,用于获取预设参数范围内整体式平底立铣刀的多个初始切削参数;每个初始切削参数包括切削厚度和切削速度,每个切削参数设置不同的切削厚度和切削速度;
后刀面磨损宽度变化率模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的后刀面磨损宽度变化率模型,具体包括:基于磨损机理的后刀面磨损体积模型建立子模块,用于根据每个初始切削参数,分别建立基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型为
Figure FDA0003810554180000061
基于几何的后刀面磨损体积模型建立子模块,用于建立基于几何的整体式平底立铣刀的后刀面磨损体积模型为
Figure FDA0003810554180000062
后刀面磨损宽度变化率模型获得子模块,用于将基于磨损机理的整体式平底立铣刀的后刀面磨损体积模型的体积等价于基于几何的整体式平底立铣刀的后刀面磨损体积模型的体积,获得每个初始切削参数的整体式平底立铣刀的后刀面磨损宽度变化率模型为
Figure FDA0003810554180000071
其中,
Figure FDA0003810554180000072
为第z个基于磨损机理的整体式平底立铣刀的后刀面磨损体积,σn-f(x")为后刀面磨损带位置x″处的正应力,w为切屑宽度,vz为第z个初始切削参数的切削速度,Δtz为第z个切削有效时间,
Figure FDA0003810554180000073
fz为每齿进给量,R为刀具半径,
Figure FDA0003810554180000074
为第z个初始切削参数的切削厚度,H0为硬质合金硬度,a为与温度有关的系数,K为切削刃上剪应力和工件上剪切流动应力的比值,
Figure FDA0003810554180000075
为第z个刀具后刀面和工件接触面的温度,Ht为铣刀表面硬度值,n为三体磨粒磨损常数,VB为铣刀后刀面磨损宽度,d为粘结磨损系数,D0为方程式系数,Q为扩散磨损过程活化能,RQ为气体常数,
Figure FDA0003810554180000076
为第z个基于几何的整体式平底立铣刀的后刀面磨损体积,λs为刀具刃倾角,
Figure FDA0003810554180000077
为第z个后刀面磨损宽度变化率,ΔVBz为第z个铣刀后刀面磨损宽度变化量,α为后角,γn为前角;
第一切削参数确定模块,用于确定每个后刀面磨损宽度变化率模型的磨损宽度的最大值,并筛选出磨损宽度的最大值大于或等于磨损宽度阈值的后刀面磨损宽度变化率模型,获得筛选出的后刀面磨损宽度变化率模型对应的初始切削参数,确定为第一切削参数;
前刀面月牙洼磨损深度变化率模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面月牙洼磨损深度变化率模型,具体包括:前刀面法向接触应力分布模型建立子模块,用于根据每个初始切削参数的切削厚度,分别建立整体式平底立铣刀的前刀面法向接触应力分布模型为
Figure FDA0003810554180000078
前刀面温度分布函数模型建立子模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的前刀面温度分布函数模型为
Figure FDA0003810554180000081
前刀面磨粒磨损率模型建立子模块,用于根据每个初始切削参数的切削速度,基于磨粒磨损机理,分别建立整体式平底立铣刀的前刀面磨粒磨损率模型为
Figure FDA0003810554180000082
前刀面粘结磨损率模型建立子模块,用于根据每个前刀面法向接触应力分布模型和每个前刀面温度分布函数模型,基于粘结磨损机理,分别建立整体式平底立铣刀的前刀面粘结磨损率模型为
Figure FDA0003810554180000083
前刀面扩散磨损率模型建立子模块,用于根据每个前刀面温度分布函数模型和每个前刀面温度分布函数模型对应的切削速度,基于扩散磨损机理,分别建立整体式平底立铣刀的前刀面扩散磨损率模型为
Figure FDA0003810554180000084
前刀面月牙洼深度变化率模型建立子模块,用于根据每个切削参数建立的前刀面磨粒磨损率模型、前刀面粘结磨损率模型和前刀面扩散磨损率模型,分别建立每个初始切削参数的整体式平底立铣刀的前刀面月牙洼深度变化率模型为
Figure FDA0003810554180000085
其中,
Figure FDA0003810554180000086
为第z个初始切削参数的切削厚度对应的前刀面法向接触应力,σtip-rake是刀尖处最大法向接触应力,
Figure FDA0003810554180000087
是第z个初始切削参数的切削厚度对应的刀屑接触长度,
Figure FDA0003810554180000088
Figure FDA0003810554180000089
为第z个初始切削参数的切削厚度对应的最大切屑厚度,
Figure FDA00038105541800000810
φn为剪切角,γn为前角,βn为摩擦角,ξ1为相关幂指数,ηc为切屑流出角,R为刀具半径,f为每齿进给量,ξ为应力指数参数,x为刀屑接触距离,
Figure FDA00038105541800000811
为第z个初始切削参数对应的前刀面温度,
Figure FDA00038105541800000812
是刀屑平均摩擦系数,
Figure FDA00038105541800000813
为第z个初始切削参数的切削速度对应的切屑速度,
Figure FDA00038105541800000814
λs为刀具刃倾角,ρ是钛合金密度,k是钛合金热导率,c是钛合金比热容,T0为室温,
Figure FDA00038105541800000815
为第z个磨粒磨损机理导致的整体式立铣刀前刀面磨粒磨损率,Qtool-chip为刀屑接触面之间的接触载荷,Hrake-face为刀具前刀面的硬度,θabrasion为磨粒平均粗糙角,
Figure FDA0003810554180000091
为第z个粘结磨损机理导致的整体式立铣刀前刀面粘结磨损率,N1为第一常数,N2为第二常数,
Figure FDA0003810554180000092
为第z个扩散磨损机理导致的整体式立铣刀前刀面扩散磨损率,Co为扩散物质的浓度,ρ为刀具材料密度,D0为方程式系数,Q为激活能,R1为气体常数,KT(x,t)为第z个初始切削参数的前刀面月牙洼深度变化率,t为切削时间;
第二切削参数确定模块,用于确定每个前刀面月牙洼磨损深度模型的磨损深度的最大值,并筛选出磨损深度的最大值大于或等于磨损深度阈值的前刀面月牙洼磨损深度模型,获得筛选出的前刀面月牙洼磨损深度模型对应的初始切削参数,确定为第二切削参数;
临界断裂时许用应力模型建立模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型,具体包括:临界断裂时许用应力模型建立子模块,用于根据每个初始切削参数,分别建立整体式平底立铣刀的临界断裂时许用应力模型为
Figure FDA0003810554180000093
其中,
Figure FDA0003810554180000094
为第z个临界断裂时许用应力,Fz为第z个初始切削参数的切削速度对应的导致刀齿破损的冲击力,
Figure FDA0003810554180000095
fz为每齿进给量,ρ为材料密度,vz为切削速度,ap是切削深度,mA为断裂面积相关比系数,VB为铣刀后刀面磨损宽度,τ为许用应力,
Figure FDA0003810554180000096
为第z个初始切削参数的切削厚度对应的损伤弹性模量值,γ为硬质合金的切应变,bv为泊松比;
第三切削参数确定模块,用于确定每个临界断裂时许用应力模型的临界断裂时许用应力,并筛选出临界断裂时许用应力大于或等于临界断裂时许用应力阈值的临界断裂时许用应力模型,获得筛选出的临界断裂时许用应力模型对应的初始切削参数,确定为第三切削参数;
失效曲线建立模块,用于在同一坐标系中分别建立第一切削参数的后刀面失效曲线、第二切削参数的前刀面失效曲线和第三切削参数的切削刃破损失效曲线;所述同一坐标系的横坐标为切削速度,纵坐标为切削厚度;
刀具失效界面图获得模块,用于将所述后刀面失效曲线、所述前刀面失效曲线和所述切削刃破损失效曲线作为三条失效界限曲线,根据三条失效界限曲线确定安全切削区域,并获得具有安全切削区域的刀具失效界面图。
CN202110002359.2A 2021-01-04 2021-01-04 一种整体式平底立铣刀的刀具失效界限图建立方法及系统 Active CN112818477B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110002359.2A CN112818477B (zh) 2021-01-04 2021-01-04 一种整体式平底立铣刀的刀具失效界限图建立方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110002359.2A CN112818477B (zh) 2021-01-04 2021-01-04 一种整体式平底立铣刀的刀具失效界限图建立方法及系统

Publications (2)

Publication Number Publication Date
CN112818477A CN112818477A (zh) 2021-05-18
CN112818477B true CN112818477B (zh) 2022-10-04

Family

ID=75856884

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110002359.2A Active CN112818477B (zh) 2021-01-04 2021-01-04 一种整体式平底立铣刀的刀具失效界限图建立方法及系统

Country Status (1)

Country Link
CN (1) CN112818477B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105893760A (zh) * 2016-04-01 2016-08-24 北京航空航天大学 基于广义磨损量的刀具磨损评价和预测方法及系统
CN109732404A (zh) * 2018-12-28 2019-05-10 哈尔滨理工大学 一种基于后刀面磨损的球头铣刀铣削力确定方法及系统

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103927450A (zh) * 2014-01-22 2014-07-16 黑龙江科技大学 一种切削刀具与高温合金组合加工时的切削参数确定方法
CN106650001B (zh) * 2016-11-15 2019-08-20 大连理工大学 一种微铣刀早期破损的预测方法
CN108549320B (zh) * 2018-03-28 2020-11-06 西北工业大学 一种基于粗糙度的钛合金铣削参数和刀具磨损控制方法
US20220176472A1 (en) * 2018-08-02 2022-06-09 Us Synthetic Corporation Cutting tool with pcd inserts, systems incorporating same and related methods
CN109571141A (zh) * 2018-11-01 2019-04-05 北京理工大学 一种基于机器学习的刀具磨损状态监测方法
CN109571142B (zh) * 2018-12-28 2020-01-17 哈尔滨理工大学 一种整体式平头立铣刀前刀面的瞬时温度确定方法及系统
CN109746766B (zh) * 2019-01-21 2022-01-25 哈尔滨理工大学 一种整体式立铣刀后刀面磨损带温度场确定方法及系统
CN110293451A (zh) * 2019-07-03 2019-10-01 哈尔滨理工大学 一种整体式平底立铣刀后刀面磨损率确定方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105893760A (zh) * 2016-04-01 2016-08-24 北京航空航天大学 基于广义磨损量的刀具磨损评价和预测方法及系统
CN109732404A (zh) * 2018-12-28 2019-05-10 哈尔滨理工大学 一种基于后刀面磨损的球头铣刀铣削力确定方法及系统

Also Published As

Publication number Publication date
CN112818477A (zh) 2021-05-18

Similar Documents

Publication Publication Date Title
Sharman et al. Workpiece surface integrity considerations when finish turning gamma titanium aluminide
Arokiadass et al. Prediction and optimization of end milling process parameters of cast aluminium based MMC
Li et al. Performance and wear analysis of polycrystalline diamond (PCD) tools manufactured with different methods in turning titanium alloy Ti-6Al-4V
Chen et al. Effects of edge preparation and feed when hard turning a hot work die steel with polycrystalline cubic boron nitride tools
Chen et al. Design, optimization and manufacturing of polycrystalline diamond micro-end-mill for micro-milling of GH4169
Ping et al. Surface integrity and tool wear mechanism of 7050-T7451 aluminum alloy under dry cutting
Li et al. Investigation and modeling of flank wear process of different PCD tools in cutting titanium alloy Ti6Al4V
CN110293451A (zh) 一种整体式平底立铣刀后刀面磨损率确定方法
Kumar Investigations into the surface quality and micro-hardness in the ultrasonic machining of titanium (ASTM GRADE-1)
CN107168245A (zh) 一种考虑刃口效应的倒棱圆形车刀切削力精确预测方法
Dou et al. Prediction model, simulation, and experimental validation on thrust force and torque in drilling SiCp/Al6063
Umer et al. On modeling tool performance while machining aluminum-based metal matrix composites
Xiao et al. Wear behavior of alumina abrasive belt and its effect on surface integrity of titanium alloy during conventional and creep-feed grinding
Xuan-Truong et al. Effect of cutting condition on tool wear and surface roughness during machining of Inconel 718
CN112818477B (zh) 一种整体式平底立铣刀的刀具失效界限图建立方法及系统
Gaitonde et al. Some studies in metal matrix composites machining using response surface methodology
Das et al. Machinability analysis of lead free brass in high speed micro turning using minimum quantity lubrication
Bhopale et al. Cutting forces during orthogonal machining process of AISI 1018 steel: Numerical and experimental modeling
Setia et al. Experimental investigation and comparison of output cutting regimes with finite element machining models for SiCp based Al7075 composite
Zahedi et al. An analytical force and surface roughness model for cylindrical grinding of brittle materials
Ge et al. Tool-chip contact characteristics during micro-cutting of compacted graphite iron (CGI) with textured tools
Necpal et al. Finite element analysis of tool stresses, temperature and prediction of cutting forces in turning process
Bharat et al. A study on conventional and non-conventional machining behaviour of metal matrix composites: a review
Zhygalov et al. Investigation of the carbide-tipped tool wear hardened by method of aerodynamic impact
Nieslony et al. 3D FEM simulation of titanium machining

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