CN109571142A - 一种整体式平头立铣刀前刀面的瞬时温度确定方法及*** - Google Patents
一种整体式平头立铣刀前刀面的瞬时温度确定方法及*** Download PDFInfo
- Publication number
- CN109571142A CN109571142A CN201811627476.2A CN201811627476A CN109571142A CN 109571142 A CN109571142 A CN 109571142A CN 201811627476 A CN201811627476 A CN 201811627476A CN 109571142 A CN109571142 A CN 109571142A
- Authority
- CN
- China
- Prior art keywords
- milling cutter
- temperature
- rake face
- deformation zone
- temperature field
- 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.)
- Granted
Links
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B23—MACHINE TOOLS; METAL-WORKING NOT OTHERWISE PROVIDED FOR
- B23Q—DETAILS, COMPONENTS, OR ACCESSORIES FOR MACHINE TOOLS, e.g. ARRANGEMENTS FOR COPYING OR CONTROLLING; MACHINE TOOLS IN GENERAL CHARACTERISED BY THE CONSTRUCTION OF PARTICULAR DETAILS OR COMPONENTS; COMBINATIONS OR ASSOCIATIONS OF METAL-WORKING MACHINES, NOT DIRECTED TO A PARTICULAR RESULT
- B23Q17/00—Arrangements for observing, indicating or measuring on machine tools
- B23Q17/09—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool
- B23Q17/0952—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining
- B23Q17/0985—Arrangements for observing, indicating or measuring on machine tools for indicating or measuring cutting pressure or for determining cutting-tool condition, e.g. cutting ability, load on tool during machining by measuring temperature
Landscapes
- Engineering & Computer Science (AREA)
- Mechanical Engineering (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种整体式平头立铣刀前刀面的瞬时温度确定方法及***。建立所述铣刀与工件之间的瞬时接触角的模型;在所述铣刀和所述工件均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度模型;根据所述未变形切屑厚度模型计算未变形切屑厚度;建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型;建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型;根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。采用分别对第一变形区热源和第二变形区热源的温度场模型的建立,提高了切削过程中铣刀温度监测的准确性。
Description
技术领域
本发明涉及整体式平头立铣刀领域,特别是涉及一种整体式平头立铣刀前刀面的瞬时温度确定方法及***。
背景技术
钛合金具有良好的比强度、比刚度、耐热性、耐腐蚀性,所有广泛应用在航空、航天、能源、建筑和军事等领域。钛合金的导热率低,因而铣削加工钛合金时,铣刀和钛合金工件之间会产生大量的切削热,瞬间导致铣刀温度急剧上升,从而加剧铣刀磨损。切削温度的变化会直接影响工件表面的质量和加工精度,所以对铣刀温度的研究尤为重要。
现有技术中针对平面周铣加工提出了一个刀-屑接触面平均温度解析计算的预报模型,并通过热电偶测量得到硬质合金铣刀前刀面的平均瞬态温度来验证模型的有效性。采用有限差分法建立直齿铣削加工过程中确定刀具温度分布的数值模型。应用有限元法对高速铣削铝合金薄壁件过程中的工件与刀具接触面温度及工件内部的温度分布进行了仿真。采用移动热源法,考虑真实的铣削加热和冷却时间,建立了符合高速铣削特点的有限元温度仿真模型。难切削材料铣削过程中刀-屑接触面的温度和热分配,利用单根热电偶被用来测量工刀-屑接触面的温度,采用响应面法对刀-屑接触面的热量分配比进行优化。还提出了一种考虑后刀面摩擦因素来描述端铣加工中工件温度循环变化的热模型,从理论和实验上表明了不同加工条件、后刀面磨损宽度和刀具位置对工件温度变化的影响。并且实现了工件温度变化的精确仿真,算法的计算时间明显短于有限元法。采用半解析的方法首次考虑了铣刀底刃对于切削温度的影响的情况下建立了一种新的端面铣削热模型,模型仿真结果和实验结果表明,温度场对不同钴粘结浓度的WC刀具的扩散磨损有非常重要的影响。采用一系列有限元数值仿真和实验并且考虑尺寸效应,开展了切削热对于微端铣加工铝合金过程影响的研究,发现随着刀尖半径的增加,刀刃上的平均切削温度会降低。考虑第一变形区的剪切热源和第二变形区的刀屑摩擦热源,提出了基于有限差分技术建立的数值分析方法的刀具温度场,发现碳化钨刀具在铣削加工时,刀具中钴含量对于其后刀面磨损影响较明显。
根据目前机械加工中刀具温度建模研究现状,研究工作主要集中在具有断续加工特点的可转位铣刀温度建模,没有涉及对于螺旋刃立铣刀前刀面瞬态温度场建模的相关研究,无法准确地确定铣刀前刀面的温度场。
发明内容
本发明的目的是提供一种能够准确确定铣刀前刀面的温度场的整体式平头立铣刀前刀面的瞬时温度确定方法及***。
为实现上述目的,本发明提供了如下方案:
一种整体式平头立铣刀前刀面的瞬时温度确定方法,所述确定方法包括:
将坐标系XYZ转换到铣刀前刀面对应的坐标系X'Y'Z';
其中,γn为法向前角,lc为刀屑接触长度;
建立所述铣刀与工件之间的瞬时接触角θ的模型;
在所述铣刀加工所述工件过程中,两者均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度h模型;
根据所述未变形切屑厚度h模型计算未变形切屑厚度h;
建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型,所述第一变形区热源为所述切屑在剪切面内的剪切滑移中产生的热量;
建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型,所述第二变形区热源为所述切屑与所述铣刀的前刀面接触产生的热量;
建立铣刀前刀面对应的坐标系X'Y'Z'下的温度场模型,获得铣刀前刀面温度场模型;
根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。
可选的,所述建立所述铣刀与工件之间的瞬时接触角θ的模型具体包括:
将所述铣刀沿着轴向切削方向离散为m个厚度为w=dz=Ap/m的圆片,Ap为铣销深度,当每个所述圆片的厚度足够小时,所述螺纹为直线BD,以所述直线BD为斜边,建立直角三角形ABD,水平方向的直角边为AB,竖直方向的直角边为AD;
根据公式(2)、(3)、(4)计算所述铣刀与工件之间的瞬时接触角θ;
其中,Ap为铣削深度,β为所述铣刀的螺旋角,R为所述铣刀的半径。
可选的,所述建立未变形切屑厚度h模型具体包括:
切削刃每次铣削均从距离已加工表面处切入
其中,
其中,vf为进给速度,Ae为铣削宽度,fz为每齿进给量,z为铣刀刃数,n为机床转速,θen为切入角,θex为切出角。
可选的,所述建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型具体包括:
建立第一变形区原生热源点dl与所述第一变形区原生热源点的镜像热源点dl'对沿着X方向上任意点PI(0,0,z)的温升模型,分别计算第一变形区原生热源点的温度升高值ΔTI(x,0,z)和所述镜像热源点的温度升高值ΔTI'(x,0,z);
计算沿着X方向上所述第一变形区的任意点PI(0,0,z)的温度升高值ΔTrake-I(0,0,z);
ΔTrake-I(x,0,z)=ΔTI(x,0,z)+ΔTI'(x,0,z) (9)
其中,λt为所述铣刀的热导率,h(θ)为未变形切屑厚度,φn法向剪切角,ηc为切屑流出角,αw为所述工件的热扩散率;K0为第二类修正贝塞尔函数的零阶,Rl为所述第一变形区原生热源的热源点dl到点PI(x,0,0)的距离,Rl'为所述第一变形区镜像热源的热源点dl到点PI(x,0,0)的距离;
建立的第一变形区热强度qI模型
Vsh=Vcosλs (12)
其中,A为参考应变率和参考温度下的初始屈服应力,B为所述工件的应***化模量,εABCD-P为剪切面ABCD的有效塑性应变,n为所述工件的应***化指数,C为所述工件的应变率强化参数,为剪切面ABCD的有效塑性应变率,为参考应变率,T为当前温度,Tr为参考温度,Tm为所述工件的熔化温度,m为所述工件的热软化指数。
可选的,所述建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型具体包括:
建立第二变形区原生热源点与第二变形区镜像热源点对沿着X'方向上任意点PII(x',0,0)的温度升高模型,分别计算所述第二变形区原生热源点的温升ΔTII(x',0,0)和所述第二变形区镜像热源点ΔTII'(x',0,0);
计算沿着X'方向上所述第二变形区热源内的任意点PII(x',0,0)的温升ΔTrake-II(x',0,0);
ΔTrake-II(x',0,0)=ΔTII(x',0,0)+ΔTII'(x',0,0) (16)
其中,Ri-2为所述第二变形区热源点到点PII(x',0,0)的距离,Ri-2'为所述第二变形区热源的镜像热源上的热源点到点PII(x',0,0)的距离,w为斜角车削宽度,λs为刃倾角;
建立所述第一变形区热强度qI模型;
Achip(θ)=lc(θ)wc (20)
其中,σtip为所述铣刀的刀尖处的正应力, 为所述铣刀和所述切屑接触面平均摩擦角,Achip(θ)为所述铣刀和所述切屑接触面积;
计算第二变形区热源流向前刀面的热量分配比
可选的,所述建立铣刀前刀面对应的坐标系X'Y'Z'下的温度场模型,获得铣刀前刀面温度场模型具体包括:
计算所述建立铣刀前刀面对应的坐标系X'Y'Z'下沿着X'方向上任意点P(x',0,0)的温度Trake(x',0,0)
Trake(x',0,0)=ΔTrake-I(lc+Z,0,0)+ΔTrake-II(x',0,0)+T0 (23)
一种整体式平头立铣刀前刀面的瞬时温度确定***,所述确定***包括:
坐标系建立模块,用于将坐标系XYZ转换到铣刀前刀面对应的坐标系X'Y'Z';
瞬时接触角模型建立模块,用于建立所述铣刀与工件之间的瞬时接触角θ的模型;
未变形切屑厚度模型建立模块,在所述铣刀加工所述工件过程中,两者均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度h模型;
未变形切屑厚度计算模块,用于根据所述未变形切屑厚度h模型计算未变形切屑厚度h;
第一温度场模型建立模块,用于建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型,所述第一变形区热源为所述切屑在剪切面内的剪切滑移中产生的热量;
第二温度场模型建立模块,用于建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型,所述第二变形区热源为所述切屑与所述铣刀的前刀面接触产生的热量;
铣刀前刀面温度场模型建立模块,用于建立铣刀前刀面对应的坐标系下的温度场模型,获得铣刀前刀面温度场模型;
瞬时温度计算模块,用于根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明公开了一种整体式平头立铣刀前刀面的瞬时温度确定方法及***,采用分别对第一变形区热源和第二变形区热源的温度场模型的建立,提高了切削过程中的铣刀温度监测的准确性,延长了所述铣刀的使用寿命,同时提高了加工工件表面的质量和加工的精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的整体式平头立铣刀前刀面的瞬时温度确定方法的流程图;
图2为本发明提供的第一变形区热源和第二变形区热源的热量分布示意图;
图3为本发明提供的所述铣刀微元刀刃示意图;
图4为本发明提供的未变形切屑厚度示意图;
图5为本发明提供的第一变形区热源示意图;
图6为本发明提供的第二变形区热源示意图;
图7为本发明提供的所述铣刀前刀面瞬态温度场的示意图;
图8为本发明提供的所述铣刀前刀面瞬态温度场随铣销时间变化示意图;
图9为本发明提供的整体式平头立铣刀前刀面的瞬时温度确定***的组成框图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种能够准确确定铣刀前刀面的温度场的整体式平头立铣刀前刀面的瞬时温度确定方法及***。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,一种整体式平头立铣刀前刀面的瞬时温度确定方法,所述确定方法包括:
步骤100:将坐标系XYZ转换到铣刀前刀面对应的坐标系X'Y'Z';
其中,γn为法向前角,lc为刀屑接触长度;
步骤200:建立所述铣刀与工件之间的瞬时接触角θ的模型;
步骤300:在所述铣刀加工所述工件过程中,两者均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度h模型;
步骤400:根据所述未变形切屑厚度h模型计算未变形切屑厚度h;
步骤500:建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型,所述第一变形区热源为所述切屑在剪切面内的剪切滑移中产生的热量;
步骤600:建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型,所述第二变形区热源为所述切屑与所述铣刀的前刀面接触产生的热量;
步骤700:建立铣刀前刀面对应的坐标系X'Y'Z'下的温度场模型,获得铣刀前刀面温度场模型;
步骤800:根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。
所述建立所述铣刀与工件之间的瞬时接触角θ的模型具体包括:
如图3所示,将所述铣刀沿着轴向切削方向离散为m个厚度为w=dz=Ap/m的圆片,Ap为铣销深度,当每个所述圆片的厚度足够小时,所述螺纹为直线BD,以所述直线BD为斜边,建立直角三角形ABD,水平方向的直角边为AB,竖直方向的直角边为AD;
根据公式(2)、(3)和(4)计算所述铣刀与工件之间的瞬时接触角θ;
其中,Ap为铣削深度,β为所述铣刀的螺旋角,R为所述铣刀的半径。
如图4所示,所述建立未变形切屑厚度h模型具体包括:
切削刃每次铣削均从距离已加工表面处切入
其中,
其中,vf为进给速度,Ae为铣削宽度,fz为每齿进给量,z为铣刀刃数,n为机床转速,θen为切入角,θex为切出角。
如图2、图5和图6所示,所述建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型具体包括:
建立第一变形区原生热源点dl与所述第一变形区原生热源点的镜像热源点dl'对沿着X方向上任意点PI(0,0,z)的温升模型,分别计算第一变形区原生热源点的温度升高值ΔTI(x,0,z)和所述镜像热源点的温度升高值ΔTI'(x,0,z);
计算沿着X方向上所述第一变形区的任意点PI(0,0,z)的温度升高值ΔTrake-I(0,0,z);
ΔTrake-I(x,0,z)=ΔTI(x,0,z)+ΔTI'(x,0,z) (9)
其中,λt为所述铣刀的热导率,h(θ)为未变形切屑厚度,φn法向剪切角,ηc为切屑流出角,αw为所述工件的热扩散率;K0为第二类修正贝塞尔函数的零阶,Rl为所述第一变形区原生热源的热源点dl到点PI(x,0,0)的距离,Rl'为所述第一变形区镜像热源的热源点dl到点PI(x,0,0)的距离;
建立的第一变形区热强度qI模型
Vsh=Vcosλs (12)
其中,A为参考应变率和参考温度下的初始屈服应力,B为所述工件的应***化模量,εABCD-P为剪切面ABCD的有效塑性应变,n为所述工件的应***化指数,C为所述工件的应变率强化参数,为剪切面ABCD的有效塑性应变率,为参考应变率,T为当前温度,Tr为参考温度,Tm为所述工件的熔化温度,m为所述工件的热软化指数。
所述建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型具体包括:
建立第二变形区原生热源点与第二变形区镜像热源点对沿着X'方向上任意点PII(x',0,0)的温度升高模型,分别计算所述第二变形区原生热源点的温升ΔTII(x',0,0)和所述第二变形区镜像热源点ΔTII'(x',0,0);
计算沿着X'方向上所述第二变形区热源内的任意点PII(x',0,0)的温升ΔTrake-II(x',0,0);
ΔTrake-II(x',0,0)=ΔTII(x',0,0)+ΔTII'(x',0,0) (16)
其中,Ri-2为所述第二变形区热源点到点PII(x',0,0)的距离,Ri-2'为所述第二变形区热源的镜像热源上的热源点到点PII(x',0,0)的距离,w为斜角车削宽度,λs为刃倾角;
建立所述第一变形区热强度qI模型;
Achip(θ)=lc(θ)wc (20)
其中,σtip为所述铣刀的刀尖处的正应力, 为所述铣刀和所述切屑接触面平均摩擦角,Achip(θ)为所述铣刀和所述切屑接触面积;
如图7和图8所示,计算第二变形区热源流向前刀面的热量分配比
所述建立铣刀前刀面对应的坐标系X'Y'Z'下的温度场模型,获得铣刀前刀面温度场模型具体包括:
计算所述建立铣刀前刀面对应的坐标系X'Y'Z'下沿着X'方向上任意点P(x',0,0)的温度Trake(x',0,0)
Trake(x',0,0)=ΔTrake-I(lc+Z,0,0)+ΔTrake-II(x',0,0)+T0 (23)
如图9所示,一种整体式平头立铣刀前刀面的瞬时温度确定***,所述确定***包括:
坐标系建立模块1,用于将坐标系XYZ转换到铣刀前刀面对应的坐标系X'Y'Z';
瞬时接触角模型建立模块2,用于建立所述铣刀与工件之间的瞬时接触角θ的模型;
未变形切屑厚度模型建立模块3,在所述铣刀加工所述工件过程中,两者均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度h模型;
未变形切屑厚度计算模块4,用于根据所述未变形切屑厚度h模型计算未变形切屑厚度h;
第一温度场模型建立模块5,用于建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型,所述第一变形区热源为所述切屑在剪切面内的剪切滑移中产生的热量;
第二温度场模型建立模块6,用于建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型,所述第二变形区热源为所述切屑与所述铣刀的前刀面接触产生的热量;
铣刀前刀面温度场模型建立模块7,用于建立铣刀前刀面对应的坐标系下的温度场模型,获得铣刀前刀面温度场模型;
瞬时温度计算模块8,用于根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (7)
1.一种整体式平头立铣刀前刀面的瞬时温度确定方法,其特征在于,所述确定方法包括:
将坐标系XYZ转换到铣刀前刀面对应的坐标系X'Y'Z';
其中,γn为法向前角,lc为刀屑接触长度;
建立所述铣刀与工件之间的瞬时接触角θ的模型;
在所述铣刀加工所述工件过程中,两者均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度h模型;
根据所述未变形切屑厚度h模型计算未变形切屑厚度h;
建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型,所述第一变形区热源为所述切屑在剪切面内的剪切滑移中产生的热量;
建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型,所述第二变形区热源为所述切屑与所述铣刀的前刀面接触产生的热量;
建立铣刀前刀面对应的坐标系X'Y'Z'下的温度场模型,获得铣刀前刀面温度场模型;
根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。
2.根据权利要求1所述的一种整体式平头立铣刀前刀面的瞬时温度确定方法,其特征在于,所述建立所述铣刀与工件之间的瞬时接触角θ的模型具体包括:
将所述铣刀沿着轴向切削方向离散为m个厚度为w=dz=Ap/m的圆片,Ap为铣销深度,当每个所述圆片的厚度足够小时,所述螺纹为直线BD,以所述直线BD为斜边,建立直角三角形ABD,水平方向的直角边为AB,竖直方向的直角边为AD;
根据公式(2)、(3)和(4)计算所述铣刀与工件之间的瞬时接触角θ;
其中,Ap为铣削深度,β为所述铣刀的螺旋角,R为所述铣刀的半径。
3.根据权利要求1所述的一种整体式平头立铣刀前刀面的瞬时温度确定方法,其特征在于,所述建立未变形切屑厚度h模型具体包括:
切削刃每次铣削均从距离已加工表面处切入
其中,
其中,vf为进给速度,Ae为铣削宽度,fz为每齿进给量,z为铣刀刃数,n为机床转速,θen为切入角,θex为切出角。
4.根据权利要求1所述的一种整体式平头立铣刀前刀面的瞬时温度确定方法,其特征在于,所述建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型具体包括:
建立第一变形区原生热源点dl与所述第一变形区原生热源点的镜像热源点dl'对沿着X方向上任意点PI(0,0,z)的温升模型,分别计算第一变形区原生热源点的温度升高值ΔTI(x,0,z)和所述镜像热源点的温度升高值ΔTI'(x,0,z);
计算沿着X方向上所述第一变形区的任意点PI(0,0,z)的温度升高值ΔTrake-I(0,0,z);
ΔTrake-I(x,0,z)=ΔTI(x,0,z)+ΔTI'(x,0,z) (9)
其中,λt为所述铣刀的热导率,h(θ)为未变形切屑厚度,φn法向剪切角,ηc为切屑流出角,αw为所述工件的热扩散率;K0为第二类修正贝塞尔函数的零阶,Rl为所述第一变形区原生热源的热源点dl到点PI(x,0,0)的距离,R'l为所述第一变形区镜像热源的热源点dl到点PI(x,0,0)的距离;
建立的第一变形区热强度qI模型
Vsh=Vcosλs (12)
其中,A为参考应变率和参考温度下的初始屈服应力,B为所述工件的应***化模量,εABCD-P为剪切面ABCD的有效塑性应变,n为所述工件的应***化指数,C为所述工件的应变率强化参数,为剪切面ABCD的有效塑性应变率,为参考应变率,T为当前温度,Tr为参考温度,Tm为所述工件的熔化温度,m为所述工件的热软化指数。
5.根据权利要求1所述的一种整体式平头立铣刀前刀面的瞬时温度确定方法,其特征在于,所述建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型具体包括:
建立第二变形区原生热源点与第二变形区镜像热源点对沿着X'方向上任意点PII(x',0,0)的温度升高模型,分别计算所述第二变形区原生热源点的温升ΔTII(x',0,0)和所述第二变形区镜像热源点ΔTII'(x',0,0);
计算沿着X'方向上所述第二变形区热源内的任意点PII(x',0,0)的温升ΔTrake-II(x',0,0);
ΔTrake-II(x',0,0)=ΔTII(x',0,0)+ΔTII'(x',0,0) (16)
其中,Ri-2为所述第二变形区热源点到点PII(x',0,0)的距离,Ri-2'为所述第二变形区热源的镜像热源上的热源点到点PII(x',0,0)的距离,w为斜角车削宽度,λs为刃倾角;
建立所述第二变形区热强度qII(θ、x')模型;
Achip(θ)=lc(θ)wc (20)
其中,σtip为所述铣刀的刀尖处的正应力, 为所述铣刀和所述切屑接触面平均摩擦角,Achip(θ)为所述铣刀和所述切屑接触面积;
计算第二变形区热源热量流向前刀面的热量分配比BII-rake(θ)
6.根据权利要求1所述的一种整体式平头立铣刀前刀面的瞬时温度确定方法,其特征在于,所述建立铣刀前刀面对应的坐标系X'Y'Z'下的温度场模型,获得铣刀前刀面温度场模型具体包括:
计算所述建立铣刀前刀面对应的坐标系X'Y'Z'下沿着X'方向上任意点P(x',0,0)的温度Trake(x',0,0)
Trake(x',0,0)=ΔTrake-I(lc+Z,0,0)+ΔTrake-II(x',0,0)+T0 (23)
7.一种整体式平头立铣刀前刀面的瞬时温度确定***,其特征在于,所述确定***包括:
坐标系建立模块,用于将坐标系XYZ转换到铣刀前刀面对应的坐标系X'Y'Z';
瞬时接触角模型建立模块,用于建立所述铣刀与工件之间的瞬时接触角θ的模型;
未变形切屑厚度模型建立模块,在所述铣刀加工所述工件过程中,两者均未变形和振动,切削刃每次铣削均从距离已加工表面处切入的条件下,建立未变形切屑厚度h模型;
未变形切屑厚度计算模块,用于根据所述未变形切屑厚度h模型计算未变形切屑厚度h;
第一温度场模型建立模块,用于建立第一变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第一温度场模型,所述第一变形区热源为所述切屑在剪切面内的剪切滑移中产生的热量;
第二温度场模型建立模块,用于建立第二变形区热源对所述铣刀前刀面产生温升影响的温度场模型,获得第二温度场模型,所述第二变形区热源为所述切屑与所述铣刀的前刀面挤压和摩擦而产生的热量;
铣刀前刀面温度场模型建立模块,用于建立铣刀前刀面对应的坐标系下的温度场模型,获得铣刀前刀面温度场模型;
瞬时温度计算模块,用于根据所述第一温度场模型、所述第二温度场模型和所述铣刀前刀面温度场模型计算所述铣刀前刀面的瞬时温度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811627476.2A CN109571142B (zh) | 2018-12-28 | 2018-12-28 | 一种整体式平头立铣刀前刀面的瞬时温度确定方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811627476.2A CN109571142B (zh) | 2018-12-28 | 2018-12-28 | 一种整体式平头立铣刀前刀面的瞬时温度确定方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109571142A true CN109571142A (zh) | 2019-04-05 |
CN109571142B CN109571142B (zh) | 2020-01-17 |
Family
ID=65933433
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811627476.2A Active CN109571142B (zh) | 2018-12-28 | 2018-12-28 | 一种整体式平头立铣刀前刀面的瞬时温度确定方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109571142B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111553032A (zh) * | 2020-04-27 | 2020-08-18 | 湖北文理学院 | 整体叶盘铣削温度预测方法、装置、设备及存储介质 |
CN112818477A (zh) * | 2021-01-04 | 2021-05-18 | 哈尔滨理工大学 | 一种整体式平底立铣刀的刀具失效界限图建立方法及*** |
CN113255120A (zh) * | 2021-05-13 | 2021-08-13 | 上海市高级技工学校 | 一种断续磨削工件边缘的温度场模型的建立方法 |
CN113369991A (zh) * | 2021-05-26 | 2021-09-10 | 武汉理工大学 | 一种后刀面磨损刀具切削温度预测方法、***及装置 |
CN114091194A (zh) * | 2021-11-16 | 2022-02-25 | 哈尔滨工程大学 | 摆线铣削中的切削温度预测方法 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0905492A2 (en) * | 1997-09-30 | 1999-03-31 | Sumitomo Electric Industries, Ltd. | Noncontact temperature distribution measuring apparatus |
CN102436527A (zh) * | 2011-11-10 | 2012-05-02 | 西北工业大学 | 平头立铣刀铣削过程铣削力建模方法 |
CN103324781A (zh) * | 2013-05-27 | 2013-09-25 | 华中科技大学 | 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法 |
CN104268343A (zh) * | 2014-09-28 | 2015-01-07 | 北京理工大学 | 一种用于端铣切削的切削力预测及温度预测的方法 |
CN105328511A (zh) * | 2015-11-18 | 2016-02-17 | 江苏科技大学 | 一种用于多孔材料切削加工在线温度采集方法 |
CN106484961A (zh) * | 2016-09-20 | 2017-03-08 | 天津大学 | 一种确定立铣刀铣削过程切削热分配比例的方法 |
DE202017105356U1 (de) * | 2017-09-05 | 2017-09-18 | Zuleadshin Enterprise Co., Ltd. | Aufbau eines intelligenten Induktionswerkzeugs |
CN107330138A (zh) * | 2017-05-17 | 2017-11-07 | 武汉工程大学 | 一种平底螺旋立铣刀的铣削力解析建模方法 |
-
2018
- 2018-12-28 CN CN201811627476.2A patent/CN109571142B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0905492A2 (en) * | 1997-09-30 | 1999-03-31 | Sumitomo Electric Industries, Ltd. | Noncontact temperature distribution measuring apparatus |
CN102436527A (zh) * | 2011-11-10 | 2012-05-02 | 西北工业大学 | 平头立铣刀铣削过程铣削力建模方法 |
CN103324781A (zh) * | 2013-05-27 | 2013-09-25 | 华中科技大学 | 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法 |
CN104268343A (zh) * | 2014-09-28 | 2015-01-07 | 北京理工大学 | 一种用于端铣切削的切削力预测及温度预测的方法 |
CN105328511A (zh) * | 2015-11-18 | 2016-02-17 | 江苏科技大学 | 一种用于多孔材料切削加工在线温度采集方法 |
CN106484961A (zh) * | 2016-09-20 | 2017-03-08 | 天津大学 | 一种确定立铣刀铣削过程切削热分配比例的方法 |
CN107330138A (zh) * | 2017-05-17 | 2017-11-07 | 武汉工程大学 | 一种平底螺旋立铣刀的铣削力解析建模方法 |
DE202017105356U1 (de) * | 2017-09-05 | 2017-09-18 | Zuleadshin Enterprise Co., Ltd. | Aufbau eines intelligenten Induktionswerkzeugs |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111553032A (zh) * | 2020-04-27 | 2020-08-18 | 湖北文理学院 | 整体叶盘铣削温度预测方法、装置、设备及存储介质 |
CN111553032B (zh) * | 2020-04-27 | 2023-09-15 | 湖北文理学院 | 整体叶盘铣削温度预测方法、装置、设备及存储介质 |
CN112818477A (zh) * | 2021-01-04 | 2021-05-18 | 哈尔滨理工大学 | 一种整体式平底立铣刀的刀具失效界限图建立方法及*** |
CN113255120A (zh) * | 2021-05-13 | 2021-08-13 | 上海市高级技工学校 | 一种断续磨削工件边缘的温度场模型的建立方法 |
CN113255120B (zh) * | 2021-05-13 | 2022-10-28 | 上海市高级技工学校 | 一种断续磨削工件边缘的温度场模型的建立方法 |
CN113369991A (zh) * | 2021-05-26 | 2021-09-10 | 武汉理工大学 | 一种后刀面磨损刀具切削温度预测方法、***及装置 |
CN114091194A (zh) * | 2021-11-16 | 2022-02-25 | 哈尔滨工程大学 | 摆线铣削中的切削温度预测方法 |
CN114091194B (zh) * | 2021-11-16 | 2024-04-30 | 哈尔滨工程大学 | 摆线铣削中的切削温度预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109571142B (zh) | 2020-01-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109571142A (zh) | 一种整体式平头立铣刀前刀面的瞬时温度确定方法及*** | |
CN110270705B (zh) | 考虑刀具跳动影响的微细铣削加工切削力仿真预测方法 | |
Ding et al. | Thermal and mechanical modeling analysis of laser-assisted micro-milling of difficult-to-machine alloys | |
Klocke et al. | Interdisciplinary modelling of the electrochemical machining process for engine blades | |
CN109746766A (zh) | 一种整体式立铣刀后刀面磨损带温度场确定方法及*** | |
Yan et al. | Modeling and analysis of coated tool temperature variation in dry milling of Inconel 718 turbine blade considering flank wear effect | |
Agmell et al. | The influence of tool micro-geometry on stress distribution in turning operations of AISI 4140 by FE analysis | |
CN107168245B (zh) | 一种考虑刃口效应的倒棱圆形车刀切削力精确预测方法 | |
Matsumura et al. | Cutting force model in milling with cutter runout | |
Niesłony et al. | Meshing strategies in FEM simulation of the machining process | |
CN110293451A (zh) | 一种整体式平底立铣刀后刀面磨损率确定方法 | |
Fu et al. | An analytical force model for ball-end milling based on a predictive machining theory considering cutter runout | |
Li et al. | Development of a hybrid cutting force model for micromilling of brass | |
CN108549320A (zh) | 一种基于粗糙度的钛合金铣削参数和刀具磨损控制方法 | |
CN111859753A (zh) | 一种负前角切削力和温度场分布预测方法及*** | |
US20050267623A1 (en) | Method for modeling material constitutive behavior | |
CN102566492A (zh) | 一种用于预测金属难加工材料插铣最大铣削力的方法 | |
Wan et al. | Investigations on the effects of different tool edge geometries in the finite element simulation of machining | |
Cheng et al. | Optimizing the geometric parameters of cutting edge for rough machining Fe-Cr-Ni stainless steel | |
Sharma et al. | In metal turning, effect of various parameters on cutting tool: A Review | |
Akyildiz | Evaluating of cutting forces in thread machining | |
Deng et al. | Finite element simulation for burr formation near the exit of orthogonal cutting | |
Duan et al. | A thermo-mechanical coupled model of derivative cutting of microtextured tools | |
Luan et al. | Comprehensive effects of tool paths on energy consumption, machining efficiency, and surface integrity in the milling of alloy cast Iron | |
Kundrák et al. | Numerical investigation of the influence of tool rake angle on residual stresses in precision hard turning |
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 |