CN107977528A - 一种用于丝杠旋风铣削加工区域的温度预测方法 - Google Patents

一种用于丝杠旋风铣削加工区域的温度预测方法 Download PDF

Info

Publication number
CN107977528A
CN107977528A CN201711366436.2A CN201711366436A CN107977528A CN 107977528 A CN107977528 A CN 107977528A CN 201711366436 A CN201711366436 A CN 201711366436A CN 107977528 A CN107977528 A CN 107977528A
Authority
CN
China
Prior art keywords
mrow
msub
msubsup
model
stage
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
Application number
CN201711366436.2A
Other languages
English (en)
Other versions
CN107977528B (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.)
Chongqing University
Original Assignee
Chongqing 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 Chongqing University filed Critical Chongqing University
Priority to CN201711366436.2A priority Critical patent/CN107977528B/zh
Publication of CN107977528A publication Critical patent/CN107977528A/zh
Application granted granted Critical
Publication of CN107977528B publication Critical patent/CN107977528B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
  • Automatic Control Of Machine Tools (AREA)

Abstract

本发明公开了一种用于丝杠旋风铣削切削加工区域的温度预测方法。首先，建立丝杠旋风铣削切削第一阶段与第二阶段瞬态未变形切屑厚度、宽度和面积模型；然后，对丝杠旋风铣削切削加工区域瞬态热源面积建模：主要包括第一变形区剪切瞬态带热源宽度模型和面积模型与第二变形区刀‑屑接触摩擦瞬态热源面积模型；最后，针对丝杠旋风铣削加工区域建立温度模型，并将所建立的关于未变形切屑方面的模型和关于热源面积方面的模型代入温度模型求解，最终获取切削加工区域中工件、切屑和刀具的瞬态温度分布。本发明能准确预测切削加工区域工件、切屑和刀具瞬态温度的分布情况，容易实现，操作方便，为丝杠旋风铣削参数优化提供了参考的依据。

Description

一种用于丝杠旋风铣削加工区域的温度预测方法
技术领域
本发明属于机械加工领域，尤其是用于预测旋风铣削加工区域温度的方法。
背景技术
传统的丝杠加工一般采用粗车、精车和磨削等加工工艺，其加工方式复杂且效率低下，并且在加工过程中采用的冷却液或切削液会对环境产生较大的污染。相比如传统的加工方式，丝杠旋风铣削具有一次成型的特点，并且具有更高的材料去除率和干式切削所带来的环保性特点。因此，采用旋风铣削的加工方式来加工丝杠能够更好的满足丝杠的生产要求和环境要求。丝杠旋风铣削的加工过程中主要有四种运动：如图1所示，刀盘1a的高速旋转运动；工件1c的进给旋转运动；刀盘1a相对工件1c的轴向进给运动；刀具1b相对工件1c的径向切削运动，丝杠旋风铣削加工过程比较复杂，导致加工过程中未变形切屑厚度、宽度和面积程瞬态变化的趋势。丝杠旋风铣削采用的是多刀断续切削的加工方式，切削加工可分为两个阶段，切削加工第一阶段与切削加工第二阶段；切削加工第一阶段为刀具开始切入工件到未变形切屑厚度达到最大厚度位置时刻，切削加工第二阶段为未变形切屑最大厚度位置时刻到刀具切出工件时刻。
不考虑刀具磨损情况下，切削加工过程中的热源主要分布在第一变形区与第二变形区，如图2所示，第一变形区主要为切削层金属发生弹性变形和塑性变形的区域，第二变形区为切屑与工件产生摩擦的区域。
丝杠旋风铣削采用干式切削，在加工过程中会产生大量的切削热，切削热的聚集会导致切削区域温度急剧上升，在促使工件表面软化利于切削的同时，切削机械应力所产生的剧烈塑性变形和摩擦，以及局部高热应力的综合作用又很大程度地影响着工件加工后的表面金相组织、残余应力分布和微观硬度，同时也加剧了刀具的磨损。虽然大量的切削热通过切屑从切削区域迅速被带走，但切削热在加工区域中的流动和传递仍将导致机床和工件局部热变形，并由此产生一定的加工误差。因此，揭示丝杠旋风铣削加工过程中切削加工区域切削热的产生及传递、掌握切削加工区域温度场的变化规律，对提升工艺性能至关重要。
目前，对于切削加工区域温度预测方法已有部分探索。例如中国专利CN105975766A公开的发明名称为《一种微铣削热力耦合解析建模方法》，考虑微铣削中存在的尺度效应，基于微铣削瞬时切削厚度模型，建立了切削温度模型。CN105930558A公开的发明名称为《一种高速切削第一和第二变形区温度的获取方法》，公开方案为通过人工热电偶方法获取刀具测量点的历史数据，根据所述温度历史建立刀具的三维传热模型、进行热分析，最后建立二维切削仿真模型，分析切削温度稳定状态下第一和第二变形区温度之间的对应关系，进而获得第一与第二变形区的温度分布情况。CN104268343A公开的发明名称为《一种用于端铣切削的切削力预测及温度预测的方法》，公开方案建立了空间任意位置有限长线热源温度场求解模型及空间任意位置有限长旋转运动线热源温度场求解模型，运用了有限元仿真方法，提出了嵌入式半人工热电偶方法，通过开展高速端铣削切削温度场分布测量对有限元仿真结果及理论计算结果进行验证及误差分析。
但是，丝杠旋风铣削区分于普通的车削与铣削加工，在加工过程中未变形切屑几何特征，包括未变形切屑厚度和宽度，呈瞬态变化趋势；这种变化情况将会导致热源热释放强度呈瞬态变化，最终引起切削温度的瞬态变化。因此，现有的热分析模型不能用来预测丝杠旋转铣削加工区域的温度分布。
发明内容
针对上述现有技术的不足，本发明提供一种用于丝杠旋风铣削加工区域的温度预测方法，解决现有技术无法针对丝杠旋风铣削加工区域进行温度预测的技术问题，能够针对旋风铣削加工的特性构建出温度预测模型，能够在不干涉切削过程的条件下进行加工区域的温度预测，摆脱对人工热电偶或半人工热电偶的依赖，只需输入相关参数即可实现温度预测，简便快速。
为解决上述技术问题，本发明的技术方案如下：一种用于丝杠旋风铣削加工区域的温度预测方法，包括以下步骤：
步骤1：分别建立切削加工的第一阶段瞬态未变形切屑厚度模型H1(θ)、第一阶段瞬态未变形切屑宽度模型w1以及第一阶段瞬态未变形切屑面积模型S1(θ)；
步骤2：分别建立切削加工的第二阶段瞬态未变形切屑厚度模型H2(θ)、第二阶段瞬态未变形切屑宽度模型w2、第二阶段瞬态未变形切屑面积模型S2(θ)以及第二阶段刀齿插入工件的最大深度模型H3(θ)；
步骤3：建立第一阶段热源模型，所述第一阶段热源模型包括第一阶段第一变形区剪切瞬态带热源宽度模型L1、第一阶段第一变形区剪切瞬态带热源面积As1、第一阶段刀-屑接触摩擦瞬态热源面积模型S1以及第一阶段刀-屑接触带热源长度L′1；
步骤4：建立第二阶段热源模型，所述第二阶段热源模型包括第二阶段第一变形区剪切瞬态带热源宽度模型L2、第二阶段第一变形区剪切瞬态带热源面积模型As2、第二阶段刀-屑接触摩擦瞬态热源面积模型S2以及第二阶段刀-屑接触带热源长度L′2；
步骤5：建立切削加工区域瞬态温度预测模型，包括工件瞬态温度模型θworkpiece、切屑瞬态温度预测模型θchip以及刀具瞬态温度预测模型θtool；
步骤6：预测切削加工区域温度分布情况：
切削加工区域中工件瞬态温度预测：
将步骤3中的第一阶段第一变形区剪切瞬态带热源宽度模型L1、第一阶段第一变形区剪切瞬态带热源面积As1以及步骤4中的第二阶段第一变形区剪切瞬态带热源宽度模型L2、第二阶段第一变形区剪切瞬态带热源面积模型As2输入步骤5中的工件瞬态温度模型θworkpiece，从而实现切削加工区域中工件的瞬态温度预测；
切削加工区域中切屑瞬态温度预测：
将步骤3中第一阶段刀-屑接触摩擦瞬态热源面积模型S1、第一阶段刀-屑接触带热源长度L′1以及步骤4中第二阶段刀-屑接触摩擦瞬态热源面积模型S2以及第二阶段刀-屑接触带热源长度L′2输入步骤5中切屑瞬态温度预测模型θchip，从而实现切削加工区域中切屑的瞬态温度的预测；
切削加工区域中刀具瞬态温度预测：
向步骤5中刀具瞬态温度预测模型θtool输入以下模型：
步骤1中第一阶段瞬态未变形切屑厚度模型H1(θ)、第一阶段瞬态未变形切屑宽度模型w1以及第一阶段瞬态未变形切屑面积模型S1(θ)；
步骤2中第二阶段瞬态未变形切屑厚度模型H2(θ)、第二阶段瞬态未变形切屑宽度模型w2、第二阶段瞬态未变形切屑面积模型S2(θ)；
步骤3中第一阶段刀-屑接触摩擦瞬态热源面积模型S1以及步骤4中第二阶段刀-屑接触摩擦瞬态热源面积模型S2；从而实现切削加工区域中刀具的瞬态温度的预测。
与现有技术相比，本发明具有以下有益效果：解决了目前关于丝杠旋风铣削加工过程中切削加工区域温度无法获取的问题。建立针对丝杠旋风铣削加工过程中切削加工第一阶段与第二阶段的瞬态未变形切屑厚度、未变形切屑宽度和未变形切屑面积的解析模型；基于此，分别建立切削加工第一阶段与第二阶段的瞬态热源模型，包括切削加工第一阶段与第二阶段的剪切瞬态带热源宽度与热源面积解析模型和切削加工第一阶段与第二阶段的刀-屑接触摩擦瞬态热源长度与热源面积解析模型；最终，将切削加工第一阶段与第二阶段的瞬态未变形切屑厚度、未变形切屑宽度、未变形切屑面积与瞬态热源模型引入切削加工过程中的瞬态温度模型，包括工件瞬态温度模型、切屑瞬态温度模型和刀具瞬态温度模型；通过瞬态温度模型可以求解丝杠旋风铣削加工过程中切削加工区域的温度分布情况。本发明能够预测丝杠旋风铣削加工过程中切削加工区域工件、切屑和刀具的温度分布，是工艺人员优化工件表面完整性、优化刀具寿命的重要技术。
附图说明
图1是丝杠旋风铣削切削加工示意图；
图2是切削加工过程中第一变形区与第二变形区分布图；
图3是刀-屑接触示意图；
图4是第n+1刀在切削第一阶段加工示意图；
图5切削第一阶段刀-屑接触区域建模；
图6第n+1刀在切削第二阶段加工示意图；
图7切削第二阶段刀-屑接触区域建模；
图8切削第一阶段摩擦热源区域边界建模；
图9切削第二阶段摩擦热源区域边界建模；
图10工件热模型；
图11切屑热模型；
图12刀具热模型
图13丝杠旋风铣削不同切削加工条件下切削加工区域温度变化情况。
具体实施方式
采用8米数控旋风铣床对丝杠进行切削加工，采用PCBN材料的刀具，刀具前角α为0°，无涂层，刀具材料PCBN导热系数λtool为0.044W/(m℃)；采用AISI 52100轴承钢的丝杠材料，密度为7.81g/cm3，热扩散系数achip为0.126J/(g℃)，导热系数λchip为0.466W/(m℃)。本文模型中所需的切削力F通过Kistler三向压电式切削测量仪获取；温度模型验证所需的实验仪器由FLIR公司提供的非接触式红外摄像仪。实验加工具体参数如下表所示：
一种用于丝杠旋风铣削加工区域的温度预测方法，包括以下步骤：
步骤1：分别建立切削加工的第一阶段瞬态未变形切屑厚度模型H1(θ)、第一阶段瞬态未变形切屑宽度模型w1以及第一阶段瞬态未变形切屑面积模型S1(θ)；
步骤2：分别建立切削加工的第二阶段瞬态未变形切屑厚度模型H2(θ)、第二阶段瞬态未变形切屑宽度模型w2、第二阶段瞬态未变形切屑面积模型S2(θ)以及第二阶段刀齿插入工件的最大深度模型H3(θ)；
步骤3：建立第一阶段热源模型，所述第一阶段热源模型包括第一阶段第一变形区剪切瞬态带热源宽度模型L1、第一阶段第一变形区剪切瞬态带热源面积As1、第一阶段刀-屑接触摩擦瞬态热源面积模型S1以及第一阶段刀-屑接触带热源长度L′1；
步骤4：建立第二阶段热源模型，所述第二阶段热源模型包括第二阶段第一变形区剪切瞬态带热源宽度模型L2、第二阶段第一变形区剪切瞬态带热源面积模型As2、第二阶段刀-屑接触摩擦瞬态热源面积模型S2以及第二阶段刀-屑接触带热源长度L′2；
步骤5：建立切削加工区域瞬态温度预测模型，包括工件瞬态温度模型θworkpiece、切屑瞬态温度预测模型θchip以及刀具瞬态温度预测模型θtool；
步骤6：预测切削加工区域温度分布情况：第一变形区的热源为剪切区热源，第二变形区的热源为摩擦热源，切削加工过程中第一变形区的剪切热源分布如图10至图11所示，第二变形区热源分布如图12所示；
切削加工区域中工件瞬态温度预测：
将步骤3中的第一阶段第一变形区剪切瞬态带热源宽度模型L1、第一阶段第一变形区剪切瞬态带热源面积As1以及步骤4中的第二阶段第一变形区剪切瞬态带热源宽度模型L2、第二阶段第一变形区剪切瞬态带热源面积模型As2输入步骤5中的工件瞬态温度模型θworkpiece，从而实现切削加工区域中工件的瞬态温度预测；
切削加工区域中切屑瞬态温度预测：
将步骤3中第一阶段刀-屑接触摩擦瞬态热源面积模型S1、第一阶段刀-屑接触带热源长度L′1以及步骤4中第二阶段刀-屑接触摩擦瞬态热源面积模型S2以及第二阶段刀-屑接触带热源长度L′2输入步骤5中切屑瞬态温度预测模型θchip，从而实现切削加工区域中切屑的瞬态温度的预测；
切削加工区域中刀具瞬态温度预测：
向步骤5中刀具瞬态温度预测模型θtool输入以下模型：
步骤1中第一阶段瞬态未变形切屑厚度模型H1(θ)、第一阶段瞬态未变形切屑宽度模型w1以及第一阶段瞬态未变形切屑面积模型S1(θ)；
步骤2中第二阶段瞬态未变形切屑厚度模型H2(θ)、第二阶段瞬态未变形切屑宽度模型w2、第二阶段瞬态未变形切屑面积模型S2(θ)；
步骤3中第一阶段刀-屑接触摩擦瞬态热源面积模型S1以及步骤4中第二阶段刀-屑接触摩擦瞬态热源面积模型S2；从而实现切削加工区域中刀具的瞬态温度的预测。
如图3至图5所示，本具体实施方式中，步骤1中各模型的表达式分别如下：
第一阶段瞬态未变形切屑厚度模型H1(θ)：
其中，θ为第n+1刀切削过程中切入点到任意切削位置时在刀尖回转中心上的夹角；Δ为刀具初始切入工件的位置角；坐标点(y1(Δ+θ),z1(Δ+θ))为辅助线ln+1与工件外圆圆周的交点P1；坐标点(y2(Δ+θ),z2(Δ+θ))为辅助线ln+1与第n+1刀刀盘回转圆周交点P2；辅助线ln+1是过点(mn+1,nn+1)与切削过程中任意切削位置所作的线段，点(mn+1,nn+1)为第n+1刀刀尖回转中心；
第一阶段瞬态未变形切屑宽度模型w1：
其中，rt为刀具圆弧半径，ρx(n+1)为第n+1刀刀尖圆弧插入工件部分所对应的圆心角，ρx(n+1)＝2arccos((rt-H1(θ))/rt)；
第一阶段瞬态未变形切屑面积模型S1(θ)：
第一阶段刀-屑接触带热源长度模型L′1：
φ为剪切角，α为前角。
如图6所示，本具体实施方式中，步骤2中各模型的表达式分别如下：
第二阶段瞬态未变形切屑厚度模型H2(θ)：
其中，坐标点(y3(Δ+θ),z3(Δ+θ))为辅助线ln+1与第n刀刀盘回转圆周的交点P3；辅助线ln是过点(mn,nn)与切削过程中任意切削位置所作的线段，点(mn,nn)为第n刀刀尖回转中心；
第二阶段刀齿插入工件的最大深度模型H3(θ)：
其中，坐标点(y4(Δ+θ),z4(Δ+θ))为辅助线ln与工件外圆圆周的交点P4；坐标点(y5(Δ+θ),z5(Δ+θ))为辅助线ln与第n刀刀盘回转圆周的交点P5；
如图7所示，第二阶段瞬态未变形切屑宽度模型w2：
其中，ρx(n)为第n刀刀尖圆弧插入工件部分所对应的圆心角，表达式为：
第二阶段瞬态未变形切屑面积模型S2(θ)：
第二阶段刀-屑接触带热源长度模型L′2：
本具体实施方式中，步骤3中各模型的表达式如下：
第一阶段第一变形区剪切瞬态带热源宽度模型L1：
其中，其中，φ为剪切角；
第一阶段第一变形区剪切瞬态带热源面积As1：
第一阶段刀-屑接触摩擦瞬态热源面积模型S1：
其中，w1为第一阶段瞬态未变形切屑宽度模型；δ11为第n+1刀在第一切削阶段过程中切削刃中心到左边界所对应的圆心角角度，δ12为第n+1刀在第一切削阶段过程中切削刃中心到右边界所对应的圆心角角度，δ11＝-ρx(n+1)/2，δ12＝ρx(n+1)/2；Lc(δ)为第n+1刀切削加工过程中刀-屑接触长度：
式中，φ为剪切角；α为前角；tc为任意位置未变形切屑厚度，tc＝rtcosδ-rt+H1(θ)，δ为第n+1刀切削刃插入工件所对应的半角。
本具体实施方式中，步骤4中各模型的表达式如下：
第二阶段第一变形区剪切瞬态带热源宽度模型L2：
第二阶段第一变形区剪切瞬态带热源面积模型As2：
第二阶段刀-屑接触摩擦瞬态热源面积模型S2：
其中，δ′11为第n+1刀在第二切削阶段过程中切削刃中心到左边界所对应的圆心角角度，δ′12为第n+1刀在第二切削阶段过程中切削刃中心到右边界所对应的圆心角角度，可分别表示为δ′21＝-ρx(n)/2，δ′22＝ρx(n)/2；L′c(δ)为第n刀切削加工过程中刀-屑接触长度，表达式如下：
δ′为第n刀切削刃插入工件所对应的半角。
本具体实施方式中，如图10所示，工件上任意一点的温升由移动带热源模型来评估，该模型基于第一变形区剪切热源与镜像热源而建立，所述工件瞬态温度模型θworkpiece：
其中，λworkpiece为工件导热系数；为倾斜角，aworkpiece为工件热扩散系数；Ko为零阶二类修正贝塞尔函数；X与Z为移动坐标系中温升点的坐标，从坐标系中取得；qs为剪平面热源热释放强度，qs＝FsVs/As，As为剪平面热源面积，Fs剪平面剪切力，Fs＝F cosφ-Ff sinφ，Vs为剪平面剪切速度，V为切削线速度，剪平面热源面积As的表达式如下：
其中，F为切削力，通过测力传感器获取；Ff为刀-屑接触面摩擦力，Ff＝0.4F。
本具体实施方式中，如图11所示，切屑温度主要由第一变形区剪切热源与第二变形区摩擦热源共同作用引起，如图8与图9所示，第二变形区的摩擦热源包括镜像热源与原热源，所述切屑瞬态温度预测模型θchip：θchip＝θchip-shear+θchip-frictional；
其中，θchip-shear为第一变形区剪平面引起的切削瞬态温度模型：
中，tch为已变形切屑厚度，tch＝tc cos(φ-α)/sinφ；
θchip-frictional为第二变形区剪平面热源引起的切屑瞬态温度模型：
其中，λchip为切屑导热系数，λchip＝λworkpiece；achip为切屑热扩散系数，achip＝aworkpiece；Bchip为刀-屑接触面上摩擦热源相对于切屑的热分配系数，取值为0.62；qr为摩擦面热源热释放强度，qr＝FfVch/Ar，Ff为刀-屑接触面摩擦力，Vch为切屑速度，Ar为刀屑接触摩擦面热源面积；其中,Ff＝0.4F，
本具体实施方式中，如图12所示，刀具温度主要由第二变形区刀屑接触摩擦热源引起，如图8与图9所示，摩擦热源包括镜像热源与原热源，刀具瞬态温度预测模型θtool包括第一阶段刀具瞬态温度预测模型θtool-1以及第二阶段刀具瞬态温度预测模型θtool-2：
第一阶段刀具瞬态温度预测模型θtool-1：
第二阶段刀具瞬态温度预测模型θtool-2：
其中，λtool为刀具导热系数；Btool为刀-屑接触面上摩擦热源相对于刀具的热分配系数；Btool为刀-屑接触面上摩擦热源相对于刀具的热分配系数，取值为0.38；f1(x,y)，…，f6(x,y)与f′1(x,y)，…，f′6(x,y)为摩擦热源边界方程，表达式分别如下：
本具体实施方式中，第n+1刀刀尖回转中心(mn+1,nn+1)，表示为：(mn+1,nn+1)＝(-ecos(Δ+η),-e sin(Δ+η))；其中，如图4所示，η为点(mn+1,nn+1)到第n+1刀切入点的连线与点(mn+1,nn+1)到工件圆心点(0,0)的连线的夹角，η通过运动过程中的三角函数关系确定：其中，e为刀尖回转中心与工件圆心的偏心距，R为刀盘回转半径，rw为工件外圆半径；
坐标点(y1(Δ+θ),z1(Δ+θ))通过以下方程组(1)获得：
坐标点(y2(Δ+θ),z2(Δ+θ))通过以下方程组(2)获得：
本具体实施方式中，第n刀刀尖回转中心(mn,nn)表示为：(mn,nn)＝(-e cos(θi+Δ+η),-e sin(θi+Δ+η))，其中，θi是从第n刀开始切削到第n+1刀开始切削这个时间段内工件转过的角度，即刀盘回转中心绕工件圆心旋转的角度，ωw为工件旋转速度，ωt为刀盘旋转速度，zt为刀盘上刀具的个数；
坐标点(y3(Δ+θ),z3(Δ+θ))通过以下方程组(3)获得：
坐标点(y4(Δ+θ),z4(Δ+θ))通过以下方程组(4)获得：
坐标点(y5(Δ+θ),z5(Δ+θ))通过以下方程组(5)获得：
其中，如图6所示，θn为当第n+1刀刀具切削工件后旋转了θ时，在相同位置处第n刀刀具切入工件后旋转的角度，表达式如下：
式中，为第n刀刀尖回转中心到坐标点(y1(Δ+θ),z1(Δ+θ))的距离：
采用本具体实施方式中用于丝杠旋风铣削加工区域的温度预测方法，所获取丝杠旋风铣削加工过程中切削加工区域切屑、工件和刀具的温度分布情况，如图13所示，图中(a)为组号1的切削加工区域温度分布情况，图中(b)为组号2的切削加工区域温度分布情况，图中(c)为组号3的切削加工区域温度分布情况，图中(d)为组号4的切削加工区域温度分布情况，切削加工区域最大温度值的实验值与预测值如下表所示：
组号 1 2 3 4
实验值(℃) 504 518 453 366
预测值(℃) 456 517 432 345
由上表可以得出通过本发明获取的加工区域最高温度随切削速度的变化规律与实验获取保持一致。切削线速度从160m/min-220m/min的变化区间内，加工区域的实验值与解析值最大温度均先增大后减小；最大温升值出现在切削线速度为180m/min的切削条件下。当切削线速度大于180m/min时，随着切削线速度的增加切削加工区域最大温度逐渐降低。由此可见，本发明所提出的用于预测丝杠旋风铣削加工过程中的温度模型具有可行性与可靠性。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (10)

1.一种用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：包括以下步骤：
步骤1：分别建立切削加工的第一阶段瞬态未变形切屑厚度模型H1(θ)、第一阶段瞬态未变形切屑宽度模型w1以及第一阶段瞬态未变形切屑面积模型S1(θ)；
步骤2：分别建立切削加工的第二阶段瞬态未变形切屑厚度模型H2(θ)、第二阶段瞬态未变形切屑宽度模型w2、第二阶段瞬态未变形切屑面积模型S2(θ)以及第二阶段刀齿插入工件的最大深度模型H3(θ)；
步骤3：建立第一阶段热源模型，所述第一阶段热源模型包括第一阶段第一变形区剪切瞬态带热源宽度模型L1、第一阶段第一变形区剪切瞬态带热源面积As1、第一阶段刀-屑接触摩擦瞬态热源面积模型S1以及第一阶段刀-屑接触带热源长度L1′；
步骤4：建立第二阶段热源模型，所述第二阶段热源模型包括第二阶段第一变形区剪切瞬态带热源宽度模型L2、第二阶段第一变形区剪切瞬态带热源面积模型As2、第二阶段刀-屑接触摩擦瞬态热源面积模型S2以及第二阶段刀-屑接触带热源长度L2′；
步骤5：建立切削加工区域瞬态温度预测模型，包括工件瞬态温度模型θworkpiece、切屑瞬态温度预测模型θchip以及刀具瞬态温度预测模型θtool；
步骤6：预测切削加工区域温度分布情况：
切削加工区域中工件瞬态温度预测：
将步骤3中的第一阶段第一变形区剪切瞬态带热源宽度模型L1、第一阶段第一变形区剪切瞬态带热源面积As1以及步骤4中的第二阶段第一变形区剪切瞬态带热源宽度模型L2、第二阶段第一变形区剪切瞬态带热源面积模型As2输入步骤5中的工件瞬态温度模型θworkpiece，从而实现切削加工区域中工件的瞬态温度预测；
切削加工区域中切屑瞬态温度预测：
将步骤3中第一阶段刀-屑接触摩擦瞬态热源面积模型S1、第一阶段刀-屑接触带热源长度L1′以及步骤4中第二阶段刀-屑接触摩擦瞬态热源面积模型S2以及第二阶段刀-屑接触带热源长度L2′输入步骤5中切屑瞬态温度预测模型θchip，从而实现切削加工区域中切屑的瞬态温度的预测；
切削加工区域中刀具瞬态温度预测θtool：
向步骤5中刀具瞬态温度预测模型θtool输入以下模型：
步骤1中第一阶段瞬态未变形切屑厚度模型H1(θ)、第一阶段瞬态未变形切屑宽度模型w1以及第一阶段瞬态未变形切屑面积模型S1(θ)；
步骤2中第二阶段瞬态未变形切屑厚度模型H2(θ)、第二阶段瞬态未变形切屑宽度模型w2、第二阶段瞬态未变形切屑面积模型S2(θ)；
步骤3中第一阶段刀-屑接触摩擦瞬态热源面积模型S1以及步骤4中第二阶段刀-屑接触摩擦瞬态热源面积模型S2；从而实现切削加工区域中刀具的瞬态温度的预测。
2.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：步骤1中各模型的表达式分别如下：
第一阶段瞬态未变形切屑厚度模型H1(θ)：
<mrow> <msub> <mi>H</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>1</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>1</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>;</mo> </mrow>
其中，θ为第n+1刀切削过程中切入点到任意切削位置时在刀尖回转中心上的夹角；Δ为刀具初始切入工件的位置角；坐标点(y1(Δ+θ),z1(Δ+θ))为辅助线ln+1与工件外圆圆周的交点；坐标点(y2(Δ+θ),z2(Δ+θ))为辅助线ln+1与第n+1刀刀盘回转圆周交点；辅助线ln+1是过点(mn+1,nn+1)与切削过程中任意切削位置所作的线段，点(mn+1,nn+1)为第n+1刀刀尖回转中心；
第一阶段瞬态未变形切屑宽度模型w1：
<mrow> <msub> <mi>w</mi> <mn>1</mn> </msub> <mo>=</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> <mo>:</mo> </mrow>
其中，rt为刀具圆弧半径，ρx(n+1)为第n+1刀刀尖圆弧插入工件部分所对应的圆心角，ρx(n+1)＝2arccos((rt-H1(θ))/rt)；
第一阶段瞬态未变形切屑面积模型S1(θ)：
<mrow> <msub> <mi>S</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
第一阶段刀-屑接触带热源长度模型L1′：
<mrow> <msubsup> <mi>L</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>H</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>;</mo> </mrow>
φ为剪切角，α为前角。
3.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：步骤2中各模型的表达式分别如下：
第二阶段瞬态未变形切屑厚度模型H2(θ)：
<mrow> <msub> <mi>H</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>;</mo> </mrow>
其中，坐标点(y3(Δ+θ),z3(Δ+θ))为辅助线ln+1与第n刀刀盘回转圆周的交点；辅助线ln是过点(mn,nn)与切削过程中任意切削位置所作的线段，点(mn,nn)为第n刀刀尖回转中心；
第二阶段刀齿插入工件的最大深度模型H3(θ)：
<mrow> <msub> <mi>H</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>5</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>y</mi> <mn>4</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>5</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>z</mi> <mn>4</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>;</mo> </mrow>
其中，坐标点(y4(Δ+θ),z4(Δ+θ))为辅助线ln与工件外圆圆周的交点；坐标点(y5(Δ+θ),z5(Δ+θ))为辅助线ln与第n刀刀盘回转圆周的交点；
第二阶段瞬态未变形切屑宽度模型w2：
<mrow> <msub> <mi>w</mi> <mn>2</mn> </msub> <mo>=</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> <mo>;</mo> </mrow>
其中，ρx(n)为第n刀刀尖圆弧插入工件部分所对应的圆心角，表达式为：
<mrow> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>=</mo> <mn>2</mn> <mi>a</mi> <mi>r</mi> <mi>c</mi> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mfrac> <mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>-</mo> <msub> <mi>H</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> </mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> </mfrac> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
第二阶段瞬态未变形切屑面积模型S2(θ)：
<mrow> <msub> <mi>S</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
第二阶段刀-屑接触带热源长度模型L2′：
<mrow> <msubsup> <mi>L</mi> <mn>2</mn> <mo>&amp;prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>H</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mi>c</mi> <mi>o</mi> <mi>s</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>.</mo> </mrow>
4.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：步骤3中各模型的表达式如下：
第一阶段第一变形区剪切瞬态带热源宽度模型L1：
其中，φ为剪切角；
第一阶段第一变形区剪切瞬态带热源面积As1：
<mrow> <msub> <mi>A</mi> <mrow> <mi>s</mi> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>)</mo> </mrow> <mo>/</mo> <mi>sin</mi> <mi>&amp;phi;</mi> <mo>;</mo> </mrow>
第一阶段刀-屑接触摩擦瞬态热源面积模型S1：
<mrow> <msub> <mi>S</mi> <mn>1</mn> </msub> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>&amp;delta;</mi> <mn>11</mn> </msub> <msub> <mi>&amp;delta;</mi> <mn>12</mn> </msub> </msubsup> <msub> <mi>w</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>L</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;delta;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;delta;</mi> <mo>;</mo> </mrow>
其中，w1为第一阶段瞬态未变形切屑宽度模型；δ11为第n+1刀在第一切削阶段过程中切削刃中心到左边界所对应的圆心角角度，δ12为第n+1刀在第一切削阶段过程中切削刃中心到右边界所对应的圆心角角度，δ11＝-ρx(n+1)/2，δ12＝ρx(n+1)/2；Lc(δ)为第n+1刀切削加工过程中刀-屑接触长度：
式中，φ为剪切角；α为前角；tc为任意位置未变形切屑厚度，tc＝rtcosδ-rt+H1(θ)，δ为第n+1刀切削刃插入工件所对应的半角。
5.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：步骤4中各模型的表达式如下：
第二阶段第一变形区剪切瞬态带热源宽度模型L2：
<mrow> <msub> <mi>L</mi> <mn>2</mn> </msub> <mo>=</mo> <msqrt> <mrow> <msup> <mrow> <mo>(</mo> <msub> <mi>y</mi> <mn>3</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>y</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>+</mo> <msup> <mrow> <mo>(</mo> <msub> <mi>z</mi> <mn>3</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>-</mo> <msub> <mi>z</mi> <mn>2</mn> </msub> <mo>(</mo> <mrow> <mi>&amp;Delta;</mi> <mo>+</mo> <mi>&amp;theta;</mi> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mrow> </msqrt> <mo>/</mo> <mi>sin</mi> <mi>&amp;phi;</mi> <mo>;</mo> </mrow>
第二阶段第一变形区剪切瞬态带热源面积模型As2：
<mrow> <msub> <mi>A</mi> <mrow> <mi>s</mi> <mn>2</mn> </mrow> </msub> <mo>=</mo> <mo>&amp;lsqb;</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <mfrac> <mn>1</mn> <mn>2</mn> </mfrac> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mrow> <mo>(</mo> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>-</mo> <msub> <mi>sin&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>)</mo> </mrow> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <mo>/</mo> <mi>s</mi> <mi>i</mi> <mi>n</mi> <mi>&amp;phi;</mi> <mo>;</mo> </mrow>
第二阶段刀-屑接触摩擦瞬态热源面积模型S2：
<mrow> <msub> <mi>S</mi> <mn>2</mn> </msub> <mo>=</mo> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>&amp;delta;</mi> <mn>11</mn> </msub> <msub> <mi>&amp;delta;</mi> <mn>12</mn> </msub> </msubsup> <msub> <mi>w</mi> <mn>1</mn> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>L</mi> <mi>c</mi> </msub> <mrow> <mo>(</mo> <mi>&amp;delta;</mi> <mo>)</mo> </mrow> <mi>d</mi> <mi>&amp;delta;</mi> <mo>-</mo> <msubsup> <mo>&amp;Integral;</mo> <msubsup> <mi>&amp;delta;</mi> <mn>21</mn> <mo>&amp;prime;</mo> </msubsup> <msubsup> <mi>&amp;delta;</mi> <mn>22</mn> <mo>&amp;prime;</mo> </msubsup> </msubsup> <msub> <mi>w</mi> <mn>2</mn> </msub> <mo>&amp;CenterDot;</mo> <msubsup> <mi>L</mi> <mi>c</mi> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>&amp;delta;</mi> <mo>)</mo> </mrow> <msup> <mi>d&amp;delta;</mi> <mo>&amp;prime;</mo> </msup> <mo>;</mo> </mrow>
其中，δ′11为第n+1刀在第二切削阶段过程中切削刃中心到左边界所对应的圆心角角度，δ′12为第n+1刀在第二切削阶段过程中切削刃中心到右边界所对应的圆心角角度，可分别表示为δ′21＝-ρx(n)/2，δ′22＝ρx(n)/2；Lc′(δ)为第n刀切削加工过程中刀-屑接触长度，表达式如下：
δ′为第n刀切削刃插入工件所对应的半角。
6.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：所述工件瞬态温度模型θworkpiece：
其中，λworkpiece为工件导热系数；为倾斜角，aworkpiece为工件热扩散系数；Ko为零阶二类修正贝塞尔函数；X与Z为移动坐标系中温升点的坐标；qs为剪平面热源热释放强度，qs＝FsVs/As，As为剪平面热源面积，Fs剪平面剪切力，Fs＝Fcosφ-Ffsinφ，Vs为剪平面剪切速度，V为切削线速度，剪平面热源面积As的表达式如下：
其中，F为切削力，通过测力传感器获取；Ff为刀-屑接触面摩擦力，Ff＝0.4F。
7.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：所述切屑瞬态温度预测模型θchip：θchip＝θchip-shear+θchip-frictional；
其中，θchip-shear为第一变形区剪平面引起的切削瞬态温度模型：
其中，tch为已变形切屑厚度，tch＝tc cos(φ-α)/sinφ；
θchip-frictional为第二变形区剪平面热源引起的切屑瞬态温度模型：
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> <mo>-</mo> <mi>f</mi> <mi>r</mi> <mi>i</mi> <mi>c</mi> <mi>t</mi> <mi>i</mi> <mi>o</mi> <mi>n</mi> <mi>a</mi> <mi>l</mi> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>q</mi> <mi>r</mi> </msub> <msub> <mi>B</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>&amp;pi;&amp;lambda;</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> </msub> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>-</mo> <msup> <mi>L</mi> <mo>&amp;prime;</mo> </msup> </mrow> </msubsup> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mrow> <mo>(</mo> <mi>X</mi> <mo>-</mo> <msub> <mi>l</mi> <mi>i</mi> </msub> <mo>-</mo> <mo>(</mo> <mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>-</mo> <msup> <mi>L</mi> <mo>&amp;prime;</mo> </msup> </mrow> <mo>)</mo> <mo>)</mo> </mrow> <msub> <mi>V</mi> <mi>c</mi> </msub> <mo>/</mo> <mn>2</mn> <msub> <mi>a</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> </msub> </mrow> </msup> <mo>&amp;lsqb;</mo> <msub> <mi>K</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msub> <mi>R</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> </msub> <msub> <mi>V</mi> <mrow> <mi>c</mi> <mi>h</mi> </mrow> </msub> <mo>/</mo> <mn>2</mn> <msub> <mi>a</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mn>0</mn> </msub> <mrow> <mo>(</mo> <msubsup> <mi>R</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> <msub> <mi>V</mi> <mrow> <mi>c</mi> <mi>h</mi> </mrow> </msub> <mo>/</mo> <mn>2</mn> <msub> <mi>a</mi> <mrow> <mi>c</mi> <mi>h</mi> <mi>i</mi> <mi>p</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>&amp;rsqb;</mo> <msub> <mi>dl</mi> <mi>i</mi> </msub> <mo>;</mo> </mrow>
其中，λchip为切屑导热系数，λchip＝λworkpiece；achip为切屑热扩散系数，achip＝aworkpiece；Bchip为刀-屑接触面上摩擦热源相对于切屑的热分配系数；qr为摩擦面热源热释放强度，qr＝FfVch/Ar，Ff为刀-屑接触面摩擦力，Vch为切屑速度，Ar为刀屑接触摩擦面热源面积；其中,Ff＝0.4F，
8.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：：刀具瞬态温度预测模型θtool包括第一阶段刀具瞬态温度预测模型θtool-1以及第二阶段刀具瞬态温度预测模型θtool-2：
第一阶段刀具瞬态温度预测模型θtool-1：
<mrow> <msub> <mi>&amp;theta;</mi> <mrow> <mi>t</mi> <mi>o</mi> <mi>o</mi> <mi>l</mi> <mo>-</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>q</mi> <mi>r</mi> </msub> <msub> <mi>B</mi> <mrow> <mi>t</mi> <mi>o</mi> <mi>o</mi> <mi>l</mi> </mrow> </msub> </mrow> <mrow> <mn>2</mn> <msub> <mi>&amp;pi;&amp;lambda;</mi> <mrow> <mi>t</mi> <mi>o</mi> <mi>o</mi> <mi>l</mi> </mrow> </msub> </mrow> </mfrac> <msubsup> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>sin</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> </mrow> <mrow> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>sin</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> </mrow> </msubsup> <mo>&amp;lsqb;</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <msub> <mi>f</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>f</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>/</mo> <msub> <mi>R</mi> <mrow> <mi>t</mi> <mi>o</mi> <mi>o</mi> <mi>l</mi> </mrow> </msub> </mrow> <mo>)</mo> </mrow> <msub> <mi>dx</mi> <mi>i</mi> </msub> <mo>+</mo> <msubsup> <mo>&amp;Integral;</mo> <mrow> <msubsup> <mi>f</mi> <mn>2</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msubsup> <mi>f</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> </mrow> </msubsup> <mrow> <mo>(</mo> <mrow> <mn>1</mn> <mo>/</mo> <msubsup> <mi>R</mi> <mrow> <mi>t</mi> <mi>o</mi> <mi>o</mi> <mi>l</mi> </mrow> <mo>&amp;prime;</mo> </msubsup> </mrow> <mo>)</mo> </mrow> <msub> <mi>dx</mi> <mi>i</mi> </msub> <mo>&amp;rsqb;</mo> <msub> <mi>dy</mi> <mi>i</mi> </msub> <mo>;</mo> </mrow>
第二阶段刀具瞬态温度预测模型θtool-2：
其中，λtool为刀具导热系数；Btool为刀-屑接触面上摩擦热源相对于刀具的热分配系数；f1(x,y)，…，f6(x,y)与f1′(x,y)，…，f6′(x,y)为摩擦热源边界方程，表达式分别如下：
<mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mn>1</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mn>2</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>cos</mi> <mi>&amp;delta;</mi> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>H</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mn>4</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mo>(</mo> <mrow> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>cos</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> </mrow> <mo>)</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>;</mo> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mn>5</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>H</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mn>6</mn> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mo>(</mo> <mrow> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>cos</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> </mrow> <mo>)</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> </mrow> </mtd> </mtr> </mtable> </mfenced>
<mrow> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msubsup> <mi>f</mi> <mn>1</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>f</mi> <mn>2</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <mfrac> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>cos</mi> <mi>&amp;delta;</mi> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mo>+</mo> <msub> <mi>H</mi> <mn>1</mn> </msub> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>f</mi> <mn>3</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>f</mi> <mn>4</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mfrac> <mrow> <mn>2</mn> <mo>(</mo> <mrow> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>cos</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> </mrow> <mo>)</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>f</mi> <mn>5</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mfrac> <mrow> <mn>2</mn> <msub> <mi>H</mi> <mn>3</mn> </msub> <mrow> <mo>(</mo> <mi>&amp;theta;</mi> <mo>)</mo> </mrow> <mi>cos</mi> <mrow> <mo>(</mo> <mrow> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> </mrow> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msubsup> <mi>f</mi> <mn>6</mn> <mo>&amp;prime;</mo> </msubsup> <mrow> <mo>(</mo> <mi>x</mi> <mo>,</mo> <mi>y</mi> <mo>)</mo> </mrow> <mo>=</mo> <mo>-</mo> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>+</mo> <mfrac> <mrow> <mn>2</mn> <mo>(</mo> <mrow> <msqrt> <mrow> <msubsup> <mi>r</mi> <mi>t</mi> <mn>2</mn> </msubsup> <mo>-</mo> <msup> <mi>y</mi> <mn>2</mn> </msup> </mrow> </msqrt> <mo>-</mo> <msub> <mi>r</mi> <mi>t</mi> </msub> <mi>cos</mi> <mfrac> <msub> <mi>&amp;rho;</mi> <mrow> <mi>x</mi> <mrow> <mo>(</mo> <mi>n</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> </mrow> </msub> <mn>2</mn> </mfrac> </mrow> <mo>)</mo> <mi>cos</mi> <mrow> <mo>(</mo> <mi>&amp;phi;</mi> <mo>-</mo> <mi>&amp;alpha;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <mi>sin</mi> <mi>&amp;phi;</mi> </mrow> </mfrac> <mo>+</mo> <mn>2</mn> <msub> <mi>r</mi> <mi>t</mi> </msub> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>.</mo> </mrow>
9.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：第n+1刀刀尖回转中心(mn+1,nn+1)，表示为：(mn+1,nn+1)＝(-ecos(Δ+η),-esin(Δ+η))；其中，η为点(mn+1,nn+1)到第n+1刀切入点的连线与点(mn+1,nn+1)到工件圆心点(0,0)的连线的夹角，η通过运动过程中的三角函数关系确定：其中，e为刀尖回转中心与工件圆心的偏心距，R为刀盘回转半径，rw为工件外圆半径；
坐标点(y1(Δ+θ),z1(Δ+θ))通过以下方程组(1)获得：
坐标点(y2(Δ+θ),z2(Δ+θ))通过以下方程组(2)获得：
10.根据权利要求1所述的用于丝杠旋风铣削加工区域的温度预测方法，其特征在于：第n刀刀尖回转中心(mn,nn)表示为：(mn,nn)＝(-e cos(θi+Δ+η),-e sin(θi+Δ+η))，其中，θi是从第n刀开始切削到第n+1刀开始切削这个时间段内工件转过的角度，即刀盘回转中心绕工件圆心旋转的角度，ωw为工件旋转速度，ωt为刀盘旋转速度，zt为刀盘上刀具的个数；
坐标点(y3(Δ+θ),z3(Δ+θ))通过以下方程组(3)获得：
坐标点(y4(Δ+θ),z4(Δ+θ))通过以下方程组(4)获得：
坐标点(y5(Δ+θ),z5(Δ+θ))通过以下方程组(5)获得：
其中，θn为当第n+1刀刀具切削工件后旋转了θ时，在相同位置处第n刀刀具切入工件后旋转的角度，表达式如下：
式中，为第n刀刀尖回转中心到坐标点(y1(Δ+θ),z1(Δ+θ))的距离：
CN201711366436.2A 2017-12-18 2017-12-18 一种用于丝杠旋风铣削加工区域的温度预测方法 Active CN107977528B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711366436.2A CN107977528B (zh) 2017-12-18 2017-12-18 一种用于丝杠旋风铣削加工区域的温度预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711366436.2A CN107977528B (zh) 2017-12-18 2017-12-18 一种用于丝杠旋风铣削加工区域的温度预测方法

Publications (2)

Publication Number Publication Date
CN107977528A true CN107977528A (zh) 2018-05-01
CN107977528B CN107977528B (zh) 2020-03-13

Family

ID=62006785

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711366436.2A Active CN107977528B (zh) 2017-12-18 2017-12-18 一种用于丝杠旋风铣削加工区域的温度预测方法

Country Status (1)

Country Link
CN (1) CN107977528B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109815563A (zh) * 2019-01-08 2019-05-28 武汉理工大学 一种基于镜像热源和非均匀热分配的三维热建模方法
CN111553032A (zh) * 2020-04-27 2020-08-18 湖北文理学院 整体叶盘铣削温度预测方法、装置、设备及存储介质
CN113177275A (zh) * 2021-04-25 2021-07-27 中国航发哈尔滨东安发动机有限公司 考虑磨削液润滑影响并联合fft技术的轴承跑道外圆磨削温度求解方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103324781A (zh) * 2013-05-27 2013-09-25 华中科技大学 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法
CN105930558A (zh) * 2016-04-13 2016-09-07 北京林业大学 一种高速切削第一和第二变形区温度的获取方法
CN105975766A (zh) * 2016-05-03 2016-09-28 大连理工大学 一种微铣削热力耦合解析建模方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103324781A (zh) * 2013-05-27 2013-09-25 华中科技大学 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法
CN105930558A (zh) * 2016-04-13 2016-09-07 北京林业大学 一种高速切削第一和第二变形区温度的获取方法
CN105975766A (zh) * 2016-05-03 2016-09-28 大连理工大学 一种微铣削热力耦合解析建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
RANGA KOMANDURI ET AL.: "Thermal modeling of the metal cutting process-Part III:temperature rise distribution due to the combined effects of shear plane heat source and the tool-chip interface frictional heat source", 《INTERNATIONAL JOURNAL OF MECHANICAL SCIENCES》 *
王立祥: "丝杠旋风铣削时变切削力的切削比能评估方法", 《中国优秀硕士学位论文全文数据库》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109815563A (zh) * 2019-01-08 2019-05-28 武汉理工大学 一种基于镜像热源和非均匀热分配的三维热建模方法
CN109815563B (zh) * 2019-01-08 2023-04-18 武汉理工大学 一种基于镜像热源和非均匀热分配的三维热建模方法
CN111553032A (zh) * 2020-04-27 2020-08-18 湖北文理学院 整体叶盘铣削温度预测方法、装置、设备及存储介质
CN111553032B (zh) * 2020-04-27 2023-09-15 湖北文理学院 整体叶盘铣削温度预测方法、装置、设备及存储介质
CN113177275A (zh) * 2021-04-25 2021-07-27 中国航发哈尔滨东安发动机有限公司 考虑磨削液润滑影响并联合fft技术的轴承跑道外圆磨削温度求解方法

Also Published As

Publication number Publication date
CN107977528B (zh) 2020-03-13

Similar Documents

Publication Publication Date Title
CN106156430B (zh) 一种基于刀具磨损效应的微铣削力建模方法
Lu et al. Model for the prediction of 3D surface topography and surface roughness in micro-milling Inconel 718
CN107977528A (zh) 一种用于丝杠旋风铣削加工区域的温度预测方法
CN103559550B (zh) 多模态耦合下的铣削稳定域预测方法
Xiao et al. A new and accurate mathematical model for computer numerically controlled programming of 4Y1 wheels in 2½-axis flute grinding of cylindrical end-mills
CN107168245B (zh) 一种考虑刃口效应的倒棱圆形车刀切削力精确预测方法
CN103324781A (zh) 一种平底螺旋立铣刀铣削工件端的三维温度场建模方法
Ferry et al. Cutter-workpiece engagement calculations by parallel slicing for five-axis flank milling of jet engine impellers
Xu et al. Time-efficient trochoidal tool path generation for milling arbitrary curved slots
Habibi et al. A generic and efficient approach to determining locations and orientations of complex standard and worn wheels for cutter flute grinding using characteristics of virtual grinding curves
Zhang et al. Finite element simulation analysis of bionic ball-end milling cutter
Fan et al. Rotary contact method for 5-axis tool positioning
Yue et al. 3D FEM simulation of milling force in corner machining process
Shen et al. Grinding wheel parametric design for machining arbitrary grooves on the helical rake face of the tool
Long∗ et al. Combined effects of flank and crater wear on cutting force modeling in orthogonal machining—part I: model development
Matsumura et al. Simulation of cutting process in peripheral milling by predictive cutting force model based on minimum cutting energy
Lasemi et al. Tool path re-planning in free-form surface machining for compensation of process-related errors
Matsumura et al. Analysis of cuttinng temperature in drilling process
CN105930654A (zh) 一种适用于难加工材料高速车削过程的切削力建模方法
Tao et al. Study on the design and cutting performance of stepped bi-directional milling cutters for hole making of CFRP
Zhang et al. Influence and prediction of tool wear on workpiece surface roughness based on milling topography analysis
Necpal et al. Finite element analysis of tool stresses, temperature and prediction of cutting forces in turning process
Yao et al. Computing cutter engagement values in milling tessellated free-form surfaces
Li et al. Modeling of instantaneous cutting force for large pitch screw with vibration consideration of the machine tool
Cheţan et al. Influence of radial finishing trajectories to the roughness obtained by milling of spherical surfaces

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