CN114429064A - 高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法 - Google Patents

高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法 Download PDF

Info

Publication number
CN114429064A
CN114429064A CN202210005504.7A CN202210005504A CN114429064A CN 114429064 A CN114429064 A CN 114429064A CN 202210005504 A CN202210005504 A CN 202210005504A CN 114429064 A CN114429064 A CN 114429064A
Authority
CN
China
Prior art keywords
cutter
face
tooth
boundary
friction
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.)
Pending
Application number
CN202210005504.7A
Other languages
English (en)
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 CN202210005504.7A priority Critical patent/CN114429064A/zh
Publication of CN114429064A publication Critical patent/CN114429064A/zh
Pending legal-status Critical Current

Links

Images

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
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • 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
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Theoretical Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Numerical Control (AREA)

Abstract

高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，属于铣刀刀齿后刀面边界检测识别技术领域。本发明包括步骤S1、构建高能效铣刀坐标系及其刀齿后刀面方程；步骤S2、构建高能效铣刀刀齿后刀面摩擦边界测量坐标系；步骤S3、对刀齿后刀面摩擦边界分形特征参数进行解算；步骤S4、表征刀齿后刀面与加工过渡表面瞬时几何接触关系；步骤S5、构建刀齿后刀面瞬时法向应力、切向应力、温度场判据；步骤S6、构建刀齿后刀面瞬时摩擦边界特征点判据；步骤S7、解算刀齿后刀面瞬时摩擦边界分形特征参数动态分布特性。本发明建立了刀齿与加工过渡表面瞬时接触关系模型，解决了已有方法忽略瞬时刀工接触关系对刀齿后刀面瞬时摩擦边界形成的影响。

Description

高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法
技术领域
本发明涉及一种刀齿后刀面摩擦边界分形特征的识别方法，属于铣刀刀齿后刀面边界检测识别技术领域。
背景技术
高能效铣刀刀齿后刀面摩擦边界的分形特征的识别是揭示铣刀刀齿后刀面摩擦边界形成过程和评价刀齿后刀面摩擦磨损的重要指标。
建立高能效铣刀刀齿后刀面摩擦边界分形特征识别方法，揭示其动态分布特性，可为高能效铣刀刀齿后刀面摩擦边界形成过程的演变机理具有指导意义。
已有关于铣刀刀齿后刀面摩擦边界解算方法，通过实验利用刀齿后刀面累积磨损最大宽度来表征。忽略了铣刀刀齿摩擦边界的不同位置处的信息，并且无法定量的揭示出整个铣削行程中铣刀刀齿后刀面摩擦边界的动态变化特征，铣刀刀齿后刀面摩擦瞬时边界与累积边界形成的过程。
因此需要一种高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，以实现铣刀刀齿后刀面摩擦边界分形特征识别，达到准确揭示出整个铣削行程中铣刀刀齿后刀面摩擦边界的动态变化特征，铣刀刀齿后刀面摩擦瞬时边界与累积边界形成的过程。
发明内容
本发明为了解决上述技术问题，在下文中给出了关于本发明的简要概述，以便提供关于本发明的某些方面的基本理解。应当理解，这个概述并不是关于本发明的穷举性概述。它并不是意图确定本发明的关键或重要部分，也不是意图限定本发明的范围。
本发明提出一种铣刀刀齿后刀面摩擦边界分形特征的解算方法，该方法通过构建铣刀刀齿及其后刀面摩擦边界测量方法；利用铣刀振动与刀齿误差对铣刀刀齿与加工过渡表面瞬时接触关系的影响，构建刀齿后刀面瞬时摩擦边界几何判据。采用有限元仿真识别刀齿后刀面瞬时摩擦边界特征点的法向应力、切向应力、温度场热力耦合判据，构建刀齿后刀面瞬时摩擦边界，并提出刀齿后刀面累积摩擦边界的解算方法；通过对铣刀刀齿瞬时摩擦边界几何特征的表征，利用分形维数与尺度系数定量的描述了铣刀刀齿后刀面瞬时摩擦边界的动态分布特性。
本发明的技术方案：
高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于，包括以下步骤：
步骤S1、构建高能效铣刀坐标系及其刀齿后刀面方程
对铣刀及其刀齿结构与后刀面方程进行求解，获得刀齿后刀面方程和刀齿切削刃方程；
步骤S2、构建高能效铣刀刀齿后刀面摩擦边界测量坐标系；
步骤S3、对刀齿后刀面摩擦边界分形特征参数进行解算；
在步骤S2构建的高能效铣刀刀齿后刀面摩擦边界测量坐标系下，截取刀齿后刀面投影面，对瞬时摩擦边界特征点进行提取，对提取变量特征参数进行解算；
步骤S4、表征刀齿后刀面与加工过渡表面瞬时几何接触关系；
揭示振动作用下铣刀刀齿后刀面瞬时几何摩擦边界的形成过程，表征不同位置处铣刀与工件的瞬时接触关系；
步骤S5、构建刀齿后刀面瞬时法向应力、切向应力、温度场判据；
步骤S6、构建刀齿后刀面瞬时摩擦边界特征点判据；
步骤S7、解算刀齿后刀面瞬时摩擦边界分形特征参数动态分布特性。
本发明具有以下有益效果：
本发明利用铣刀刀齿坐标系与其后刀面测量坐标系关系，解决了已有方法无法定量表征刀齿后刀面摩擦边界的问题。建立了刀齿与加工过渡表面瞬时接触关系模型，解决了已有方法忽略瞬时刀工接触关系对刀齿后刀面瞬时摩擦边界形成的影响。
本发明构建了热力耦合场仿真模型，识别出刀齿后刀面瞬时摩擦边界的判据，解决了刀齿后刀面瞬时摩擦边界无法获取的问题，该方法提出了刀齿后刀面瞬时摩擦边界分形特征的动态分布特性，解决了刀齿后刀面瞬时摩擦边界描述的不确定性问题，该方法提出了刀齿后刀面累积摩擦边界的形成问题，解决了已有方法无法定量描述累积摩擦边界形成问题，该方法可用于揭示铣刀切削过程中刀齿后刀面摩擦动态形成过程，并为高能效铣刀刀齿后刀面摩擦边界形成机理提供依据。
本发明与与已经公开的技术不同之处：
已有关于铣刀刀齿后刀面摩擦边界的解算方法，利用铣刀刀齿后刀面摩擦边界的最大磨损宽度表征摩擦边界，并基于铣刀切削整个铣削行程后进行摩擦边界的度量。
上述方法以刀齿后刀面最大磨损宽度表征摩擦结果，忽略了铣刀刀齿摩擦边界的不同位置处的信息，并且无法定量的揭示出整个铣削行程中铣刀刀齿后刀面摩擦边界的动态变化特征，铣刀刀齿后刀面摩擦瞬时边界与累积边界形成的过程。
本方法考虑了刀齿与铣刀的结构与铣刀动态切削过程中刀工瞬时接触状态，提出高能效铣刀刀齿后刀面摩擦边界分形特征参数的解算方法，构建铣刀刀齿后刀面摩擦边界表征与检测模型，并利用高能效铣刀铣削有限元热力耦合仿真，解算铣刀从切入到切出过程中刀齿后刀面不同位置处几何特征参数、法向应力、切向应力、温度场动态分布特性，利用铣刀刀齿后刀面瞬时摩擦上下边界的识别判据，并通过刀齿后刀面摩擦边界分形特征参数表征方法定量的揭示了铣刀刀齿后刀面摩擦边界动态分布特性。
附图说明
图1是高能效铣刀及刀齿结构示图；
图2是刀齿后刀面的结构示图；
图3是图2中A-A向剖视图；
图4是图3中B-B向剖视图；
图5是图2中C向示意图；
图6是高能效铣刀刀齿后刀面摩擦边界测量方法注释图；
图7是刀齿后刀面瞬时摩擦边界曲线图；
图8是刀齿后刀面瞬时摩擦边界分形特征参数的双对数函数曲线图；
图9是铣削加工现场图；
图10是铣削行程5m的高能效铣刀刀齿底刃后刀面磨损状态(刀齿一后刀面磨损状态)；
图11是铣削行程5m的高能效铣刀刀齿底刃后刀面磨损状态(刀齿二后刀面磨损状态)；
图12是铣削行程5m的高能效铣刀刀齿底刃后刀面磨损状态(刀齿三后刀面磨损状态)；
图13是铣削行程5m的高能效铣刀刀齿后刀面瞬时摩擦上边界曲线图；
图14是铣削行程5m的高能效铣刀刀齿后刀面瞬时摩擦下边界曲线图；
图15是振动作用下铣刀及其刀齿的瞬时位姿图；
图16是刀齿后刀面与加工过渡表面瞬时接触关系图；
图17是铣刀刀齿后刀面瞬时几何接触边界分布；
图18是铣刀后刀面上特征点的有限元仿真结果与等效应力分解模型图；
图19是铣刀刀齿后刀面瞬时法向应力沿U方向场分布；
图20是铣刀刀齿后刀面瞬时切向应力沿U方向场分布；
图21是铣刀刀齿后刀面瞬时温度沿U方向场分布；
图22是铣刀刀齿后刀面瞬时上边界节点判据图；
图23是铣刀刀齿后刀面瞬时法向应力变化率；
图24是铣刀刀齿后刀面瞬时切向应力变化率；
图25是铣刀刀齿后刀面瞬时切向应力变化率；
图26是铣刀刀齿后刀面瞬时摩擦上边界曲线；
图27是铣刀刀齿后刀面瞬时下边界判据突变节点边界曲线；
图28是t1时刻刀齿一后刀面边界形成过程；
图29是t2时刻刀齿二后刀面边界形成过程；
图30是t3时刻刀齿三后刀面边界形成过程；
图31是刀齿后刀面上边界分形特征参数D(t)变化过程；
图32是刀齿后刀面上边界分形特征参数G(t)变化过程；
图33是刀齿后刀面下边界分形特征参数D(t)变化过程；
图34是刀齿后刀面下边界分形特征参数G(t)变化过程；
图35是刀齿后刀面累积摩擦上边界解算方法；
图36是刀齿后刀面累积摩擦下边界解算方法；
图37是刀齿后刀面仿真累积摩擦边界；
图38是t1时刻刀齿一后刀面边界形成过程；
图39是t2时刻刀齿二后刀面边界形成过程；
图40是t3时刻刀齿三后刀面边界形成过程；
图41是刀齿后刀面上边界分形特征参数D(t)形成过程；
图42是刀齿后刀面上边界分形特征参数G(t)形成过程；
图43是刀齿后刀面下边界分形特征参数D(t)形成过程；
图44是刀齿后刀面下边界分形特征参数G(t)形成过程；
图45是方案2刀齿后刀面仿真累积摩擦边界解算结果；
图46是本发明的方法流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面通过附图中示出的具体实施例来描述本发明。但是应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
具体实施方式一：高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，具体包括以下步骤：
步骤S1、构建高能效铣刀坐标系及其刀齿后刀面方程；
对铣刀及其刀齿结构与后刀面方程进行求解，获得刀齿后刀面方程和刀齿切削刃方程；
步骤S2、构建高能效铣刀刀齿后刀面摩擦边界测量坐标系；
步骤S3、对刀齿后刀面摩擦边界分形特征参数进行解算；
在步骤S2构建的高能效铣刀刀齿后刀面摩擦边界测量坐标系下，截取刀齿后刀面投影面，对瞬时摩擦边界特征点进行提取，对提取变量特征参数进行解算；
步骤S4、表征刀齿后刀面与加工过渡表面瞬时几何接触关系；
揭示振动作用下铣刀刀齿后刀面瞬时几何摩擦边界的形成过程，表征不同位置处铣刀与工件的瞬时接触关系；
步骤S5、构建刀齿后刀面瞬时法向应力、切向应力、温度场判据；
步骤S6、构建刀齿后刀面瞬时摩擦边界特征点判据；
步骤S7、解算刀齿后刀面瞬时摩擦边界分形特征参数动态分布特性。
在本实施例的步骤S1中，高能效铣刀坐标系及刀齿后刀面方程构建方法，具体是：
(1)为揭示振动作用下刀齿后刀面与加工过渡表面的瞬时动态特性，首先对铣刀及其刀齿结构与后刀面方程进行求解，如图1-图5所示；
图1-图5中，D为铣刀刀柄直径，l为铣刀悬伸量，L为铣刀总长度；ezmin为铣刀刀齿轴向最低刀尖点，e1为铣刀回转半径最大的刀尖点，rmax为铣刀刀齿最大回转半径；o0-x0y0z0为铣刀结构坐标系，其中，o0为铣刀轴向最低刀尖点的回转中心，与铣刀轴向最低刀尖点共面；x0轴与回转半径最大刀尖点切削速度方向平行；y0轴与回转半径最大刀尖点的径向平行；zd轴为铣刀回转轴，并指向刀柄方向。
oi-xiyizi为刀齿i坐标系，其中oi为刀齿刀尖点外接圆r1的中心，yi轴为过原点oi垂直与刀齿底面的方向；zi轴通过oi并指向刀尖点；xi经过oi并同时垂直于yi与zi，η为zi与z0轴之间的安装角；e为铣刀第i个刀齿后刀面上的任一点，ri为刀齿i的回转半径，Δri为刀齿i的径向误差，Δzi为刀齿i的轴向误差；
Figure BDA0003456527400000051
为y0轴与o0e1在x0o0y0平面的夹角。
em为刀齿切削刃中点，R0为em点在xioizi平面上的外切圆半径，oz为em点外切圆的中心点；在xioiyi平面上，s为刀齿后刀面上任一点e到xi轴的距离，s0为刀齿切削刃中点em到xi轴的距离，s1为刀尖点e1到xi轴的距离；z0为切削刃中点外接圆中心到刀齿坐标系原点oi在xioizi平面上的距离；γ为oie1与oiem在xioizi平面上夹角，R为B-B剖面下刀尖点在xioizi平面上的外切圆半径；δ为B-B剖面下直线oze与ozoi在xioizi平面上的夹角，δ0为B-B剖面下刀尖到oz与ozoi在xioizi平面上的夹角，α在yioizi平面上的后角；R1为切削刃在xioiyi平面上所在弧的外切圆，β为C视图下oi'e与oi'em在xioiyi平面上的夹角，β0为oi'e1与oi'em在xioiyi平面上的夹角。
(2)刀齿i坐标系oi-xiyizi与铣刀结构坐标系o0-x0y0z0的旋转矩阵Q1、Q2、平移矩阵M1分别为：
Figure BDA0003456527400000061
Figure BDA0003456527400000062
其中ri＝rmax-Δri
(3)刀齿坐标系下后刀面方程求解方法为：
在刀齿坐标系中，刀齿i的后刀面方程为：
Figure BDA0003456527400000063
其中，
0≤s≤s0,-δ0≤δ≤δ0 (4)
R＝R0-s0·tanα+s·tanα (5)
Figure BDA0003456527400000064
由式(3)满足如下条件可获得刀齿i的切削刃方程：
R·sinδ＝R1·sinβ (7)
s＝R1 cosβ-(R1-s0) (8)
s1≤s≤s0，-β0≤β≤β0 (9)
Figure BDA0003456527400000065
。
在本实施例的步骤S2中，高能效铣刀刀齿后刀面摩擦边界测量坐标系构建方法，具体是：
(1)刀齿坐标系下刀齿后刀面摩擦边界测量方法，如图6所示：图中，om-UVW为刀齿后刀面测量坐标系，过切削刃中点em做垂直于刀齿底部所在平面Ⅱ与yi轴交于点y1，其中om的xi方向坐标与刀尖点en的xi方向坐标一致，yi方向坐标与的y1方向坐标一致，zi方向坐标与的oi方向坐标一致；U轴过y1并平行于xi轴，V轴过om并平行于zi轴，W轴过om并平行于yi轴；f(U，V)＝0为切削刃原始边界，f(U，Vs)＝0为摩擦上边界，f(U，Vp)＝0为摩擦上边界；e0为切削刃中点em在投影面Ⅰ上的投影。
(2)刀齿测量坐标系om-UVW与刀齿i坐标系oi-xiyizi的旋转矩阵Q7、Q8、平移矩阵M4分别为：
Figure BDA0003456527400000071
Figure BDA0003456527400000072
。
在本实施例的步骤S3中，对刀齿后刀面摩擦边界分形特征参数进行解算方法，具体是：
分形特征参数是用来描述摩擦边界变化稳定状态的参数，表征摩擦边界变化的波动程度。
(1)刀齿后刀面投影面摩擦边界为：
在图7中后刀面摩擦边界测量坐标系下，做若干平行于V轴且间距为ΔU的平行线截取刀齿后刀面投影面，α为ΔU的整数倍；分别与原始刃、上边界、下边界交于一点，其中e1、e2、e3、e4分别为切削刃原始边界、摩擦上边界、摩擦下边界、摩擦区域上的点；在投影坐标系下，Ui0为刀齿后刀面摩擦边界U向最小值为Uie与Ve分别为切削刃中点e0的U，V轴的坐标，Uia与Uib分别为沿U轴方向的任意位置取样线的坐标，Uin与Vn分别为切削刃原始边界最大U向坐标值与刀尖点V向坐标值。
(2)刀齿后刀面摩擦边界分形特征参数解算方法为：
对瞬时摩擦边界特征点进行提取，将铣刀刀齿后刀面摩擦边界曲线U-V视作一个空间序列，对提取变量特征参数进行解算：
Figure BDA0003456527400000073
式中，S(γ)是结构函数的测度，即表示差方算数平均根[V(U0+γ)-V(U0)]2/γ，γ是采样间隔的尺度，是数据间隔的任意取值，S(ω)是功率谱密度函数；D是分形维数，G为尺度系数。
针对若干γ对曲线的离散信号，计算出相应的S(γ)，在双对数坐标logS(γ)-logγ中，通过matlab运用最小二乘法对其进行拟合，直线斜率k与截距A，换算得到表面轮廓曲线的分形维数D和尺度系数G，其中D与斜率k的转换关系为式(14)。
D＝2-k/2 (14)
如图8刀齿后刀面瞬时摩擦边界分形特征参数的双对数函数曲线，其中分形维数D反映轮廓变化的相似性与复杂性，尺度系数G反映曲线轮廓的不平整度。
(3)实验中采用的Walter公司生产的三齿等齿距直径为32mm的可转位高进给铣刀，铣削方案一如表1所示，采用型号为TC4的钛合金材料，钛合金工件的长宽高为250×100×20mm，铣削方案如表1所示。
表1高能效铣刀铣削钛合金实验方案
Figure BDA0003456527400000081
表中，Δzi为刀齿轴向误差，Δri为刀齿径向误差，Δθi为齿间夹角误差。
采用实验方案一，实验结束后，使用超景深显微镜，依次对铣刀各刀齿后刀面摩擦边界进行检测，如图10、图11、图12、图13和图14所示；
分形维数D反映刀齿后刀面摩擦边界的波动程度，尺度系数G反映出摩擦边界曲线幅值变化的不平整度；由表2可知，通过分形特征参数的表征方法可以看出根据刀齿的误差不同，铣刀刀齿后刀面摩擦边界呈现出差异性。
表2刀齿后刀面累积摩擦边界分形特征参数实验结果
Figure BDA0003456527400000082
。
在本实施例的步骤S4中，表征刀齿后刀面与加工过渡表面瞬时几何接触关系，高能效铣刀刀齿后刀面摩擦边界判据识别方法，具体是：
1.刀齿后刀面与加工过渡表面瞬时接触关系几何判据解算方法
(1)铣刀及其刀齿与加工过渡表面的瞬时接触关系直接影响铣刀刀齿后刀面摩擦边界的演变特性。为揭示振动作用下铣刀刀齿后刀面瞬时几何摩擦边界的形成过程，表征不同位置处铣刀与工件的瞬时接触关系，如图15振动作用下铣刀及其刀齿的瞬时位姿所示；
图中，L、S、H分别为工件的长度、宽度、高度，o-xyz为工件坐标系；n为铣刀的转速，vf为铣刀名义进给速度，ae为铣刀的切削宽度，ap为铣刀的切削深度，od-xdydzd为无振动作用下的铣刀切削坐标系；Ax(t)、Ay(t)、Az(t)分别为铣刀振动沿x、y、z三个方向上的位移，od'-xd'yd'zd'为铣削振动引起的坐标系od-xdydzd偏置后的铣刀切削坐标系；od'与o0重合，xd'轴为铣刀中心o0运动轨迹的切矢量方向，zd'轴与z0轴重合；od'(x,y,z)为坐标原点od'的运动轨迹，vx、vy、vz分别为od'(x,y,z)的切矢量沿x轴、y轴、z轴的分量；θ(t)为t时刻z0轴与zd轴的夹角，θs(t)为t时刻zi轴与刀齿中心轴的夹角。
(2)铣刀结构坐标系o0-x0y0z0与振动作用下的铣刀切削坐标系od'-xd'yd'zd'的瞬时旋转矩阵Q3为：
Figure BDA0003456527400000091
其中，
Figure BDA0003456527400000092
为y0轴与yd'轴在xd'od'yd'平面内的瞬时夹角：
Figure BDA0003456527400000093
式中，
Figure BDA0003456527400000094
为铣刀初始切入t＝0时刻y0轴与yd'轴在xd'od'yd'平面内的夹角。
(3)振动作用下的铣刀切削坐标系od'-xd'yd'zd'与无振动作用下的铣刀切削坐标系od-xdydzd的瞬时旋转矩阵Q4、Q5，瞬时平移矩阵M2分别为：
Figure BDA0003456527400000095
Figure BDA0003456527400000096
其中，θ1(t)为铣刀瞬时姿态角在xdodzd面的投影，θ2(t)为铣刀瞬时姿态角在ydodzd面的投影与xo轴的瞬时夹角，θ(t)为zd'轴与zd轴的瞬时夹角：
Figure BDA0003456527400000097
zd'轴与zd轴的瞬时夹角θ(t)为：
Figure BDA0003456527400000098
(4)无振动作用下的铣刀切削坐标系od-xdydzd与工件坐标系o-xyz的瞬时平移矩阵M3为：
Figure BDA0003456527400000101
(5)刀齿i的坐标系oi-xiyizi与工件坐标系o-xyz变换矩阵Φi为：
Φi＝M3·M2·Q5·Q4·Q3·M1·Q2·Q1 (22)
(6)刀齿i的中心点oi在工件坐标系中的运动轨迹oi(x,y,z)为：
oi(x,y,z)＝[x y z 1]T＝Φi·[0 0 0 1]T (23)
根据刀齿切削刃在工件中瞬时运动关系，解算刀齿从切入到切出加工过渡表面，其刀齿与加工过渡表面的瞬时接触关系如图16所示；
(7)在工件坐标系中，刀齿i的切削刃方程I(x,y,z)为：
Figure BDA0003456527400000102
(8)解算刀齿i瞬时切入至切出加工过渡表面方程Ki为：
Figure BDA0003456527400000103
其中，I(x,y,z)为刀齿i切削过程中，由t1 i至t2 i时段内，在工件坐标系下的切削刃方程，t1 i为刀齿t2 i切入待加工表面的初始时刻，为刀齿i切出待加工表面时刻；ei为刀齿后刀面摩擦区域中任一点，。
(9)图16中e(x1，y1，z1)点为铣削加工过程中刀齿后刀面上任一点，做其后刀面所在的切平面的法矢量
Figure BDA0003456527400000104
与加工过渡表面相交与p(x2，y2，z2)点为：
Figure BDA0003456527400000105
如图16所示，在刀齿后刀面上选取点e(x1，y1，z1)，过e所在的切削平面向加工过渡表面方向做法矢量
Figure BDA0003456527400000106
与加工过渡表面交于点p(x2，y2，z2)，设
Figure BDA0003456527400000107
以刀齿后刀面的特征点指向刀齿外的方向为正方向，其大小为
Figure BDA0003456527400000108
那么其瞬时接触关系判定：
1)当
Figure BDA0003456527400000109
时，
Figure BDA00034565274000001010
为正方向，e与p处相切，属于接触临界点；
2)当
Figure BDA00034565274000001011
时，
Figure BDA00034565274000001012
为负方向，e与p未产生接触；
3)当
Figure BDA00034565274000001013
时，
Figure BDA00034565274000001014
为正方向，e与p处产生摩擦，磨损。
对刀齿后刀面上法矢量几何判定边界进行计算，结果如图17所示。
在本实施例的步骤S5中，构建刀齿后刀面瞬时法向应力、切向应力、温度场判据，具体方法是：
(1)通过Deform仿真后，对后刀面上特征点的有限元仿真结果进行提取，其仿真结果与等效应力分解模型如图18所示；
图18中，e为刀齿后刀面上任一点，σq 1，σq 2，σq 3分别为等效应力沿正四面体的三个方向力；
Figure BDA0003456527400000111
为过e所在的切削平面向加工过渡表面方向做法矢量，σ与τ分别为刀齿后刀面点e处的所受的法向应力与切向应力，
Figure BDA0003456527400000112
为切向应力τ方向为相对运动速度vvn方向在刀齿后刀面与加工过渡表面的公切面投影的方向，vn为摩擦速度；
Figure BDA0003456527400000113
分别为刀齿后刀面上特征点e处等效应力三个分量与法向应力σ的矢量方向的矢量夹角，
Figure BDA0003456527400000114
Figure BDA0003456527400000115
分别为分别为刀齿后刀面上特征点e处等效应力三个分量与法向应力τ的矢量方向的矢量夹角。
(2)在工件坐标系下，对仿真结果提取特征节点的等效应力，将其等效应力沿正四面体的三个方向力σq 1，σq 2，σq 3进行分解。并将等效应力三个分量在特征点ei(x1，y1，z1)，沿着其切削平面的法矢量
Figure BDA00034565274000001112
的投影的方向，在刀齿坐标系下，刀齿后刀面上特征点e(x1，y1，z1)处等效应力三个分量与法向应力σ的矢量方向的矢量夹角：
Figure BDA0003456527400000116
其中，
Figure BDA0003456527400000117
为等效应力沿正四面体三个方向分力，
Figure BDA0003456527400000118
为等效应力三个方向分量与法向应力间的夹角。法向应力σ为：
Figure BDA0003456527400000119
式中
Figure BDA00034565274000001110
σj为等效应力三个方向分量在法向应力方向上的分量；
(3)铣刀沿od(x,y,z)运动轨迹的切矢量x，y，z三个坐标方向上的运动速度的向量表达
Figure BDA00034565274000001111
为：
Figure BDA00034565274000001113
其中，vx(t)、vy(t)、vz(t)分别为后刀面任意点的运动速度大小vvn(t)在x、y、z三个坐标方向上的分量。
(4)将摩擦速度通过矩阵变换转换到铣刀刀齿坐标系内为：
Figure BDA0003456527400000121
(5)铣刀刀齿后刀面的任一点的切向应力τ方向为相对运动速度vvn方向在刀齿后刀面与加工过渡表面的公切面投影的方向
Figure BDA0003456527400000122
为：
Figure BDA0003456527400000123
在工件坐标系下，刀齿后刀面特征点e(x1，y1，z1)处等效应力三个分量与切向应力τ的矢量方向的矢量夹角如式(32)为：
Figure BDA0003456527400000124
其中，
Figure BDA0003456527400000125
为等效应力三个方向分量与切向应力τ间的夹角。切向应力τ为：
Figure BDA0003456527400000126
式中
Figure BDA0003456527400000127
τj为等效应力三个方向分量在切向应力方向上的分量。
(6)铣刀刀齿切削刃的瞬时法向应力、切向应力、温度场分布，如图19、图20和图21所示。
在本实施例的步骤S6中，构建刀齿后刀面瞬时摩擦边界特征点判据，即刀齿后刀面瞬时摩擦边界特征点判据的识别方法是：
(1)如图22，刀齿后刀面瞬时摩擦上边界识别以刀齿材料等效应力与屈服强度临界值作为判据为：
σq≥σs
其中，σq为刀齿后刀面等效应力特征点数值，σs为屈服强度的大小。
(2)刀齿后刀面瞬时摩擦下边界识别出后刀面各截面位置几何特征、温度场、法向应力、切向应力的变化率突变的特征节点为：
Figure BDA0003456527400000128
其中，T、σ、τ分别温度、法向应力、切向应力，T(U,V)、σ(U,V)、τ(U,V)分别为刀齿后刀面同一时刻瞬时温度、法向应力、切向应力特征点数值；T'、σ'、τ'分别为齿后刀面特征点瞬时温度、法向应力、切向应力变化率，如图23、图24和图25所示。
在本实施例的步骤S7中，解算刀齿后刀面瞬时摩擦边界分形特征参数动态分布特性，具体方法是：
1.刀齿后刀面瞬时摩擦边界的解算方法
(1)在后刀面摩擦磨损测量坐标系中，截平面上特征节点的屈服强度所对应位置即为刀齿后刀面瞬时摩擦磨损上边界，识别出t时刻刀齿后刀面上所有截平面上等效应力大于屈服强度的临界特征节点作为刀齿后刀面瞬时摩擦上边界如图26所示。
Figure BDA0003456527400000131
其中，σq为刀齿后刀面等效应力特征点数值，σs为屈服强度的大小；Vs(U，t)为刀齿后刀面瞬时摩擦上边界，Vσs(U，t)为刀齿后刀面瞬时摩擦边界屈服强度曲线Vs(U，t-Δt)为上一时刻铣刀刀齿后刀面瞬时摩擦边界。
(2)根据判据结果识别出铣刀刀齿t时刻后刀面不同位置处几何接触关系、法向应力、切向应力、温度变化率发生突变的节点位置，并在测量坐标系下分别表征出瞬时特征点的变化曲线。在测量坐标系下，识别出沿切削刃特征点切矢量方向的不同判据的最大值Vmax，作为刀齿后刀面瞬时摩擦下边界判据曲线如图27所示。
Vξ(U，t)＝max{VN(U，t),Vσ(U，t),VT(U，t),Vτ(U，t)} (36)
其中，V0(U)为刀齿后刀面原始刃边界，VN(U，t)，Vσ(U，t)，Vτ(U，t)，VT(U，t)分别为t时刻刀齿后刀面上几何接触、法向应力、切向应力、温度判据识别曲线，Vξ(U，t)为t时刻铣刀刀齿后刀面瞬时摩擦边界。
(3)通过分形特征参数表征方法将几何接触、法向应力、切向应力、温度场判据与Vξ(U，ti)的分形特征参数进行解算，结果如表3所示。
表3刀齿后刀面累积摩擦边界分形特征参数表征
Figure BDA0003456527400000132
(4)根据上述瞬时摩擦边界求解方法，取切削行程5m最后一层铣刀刀齿切入时段、切削中段、切出时段三个连续周期t1＝18.6s，t2＝12.4s，t3＝24.8s瞬时摩擦边界结果如图28、图29和图30所示。
2.刀齿后刀面瞬时摩擦边界分形特征参数动态分布特性的解算方法
在表1实验方案下识别出切削5m行程最后一次切削行程下，铣刀三个刀齿切入到切出不同周期相同位置角的瞬时边界变化曲线，根据分形特征参数表征方法分别对其进行表征，识别出铣刀三个刀齿后刀面瞬时摩擦边界的分形特征参数的变化过程，如图31-图34所示。
由图31-图34可知，由于振动与误差的影响，各个刀齿后刀面瞬时摩擦边界分形维数与尺度系数随刀齿位置角变化曲线总体呈现出相近的变化趋势，但其变化特性有所不同；刀齿后刀面瞬时摩擦上边界分形维数D(t)逐渐增大，说明切削刃的损伤程度与其力热耦合场复杂，靠近切削变形区导致其多变性；刀齿后刀面瞬时摩擦下边界分形维数D(t)逐渐将减小，说明摩擦下边界由于热力耦合场的扩展不断趋于平稳；刀齿后刀面瞬时摩擦边界尺度系数G(t)均逐渐增大，说明其边界的变动幅度逐渐增大，导致摩擦边界具有不确定性。
具体实施方式二：高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，还包括铣刀刀齿后刀面累积摩擦边界解算及其分形特征验证方法，具体实施方法是：
(1)将刀齿i分别在铣刀切入、切中、切出分别取不同刻每个周期内瞬时摩擦最大摩擦边界形成的上边界曲线进行识别与表征如图35所示。
根据上述分析，铣刀刀齿后刀面累积摩擦上边界为式(37)所示。
Vs(U)＝Vs(U,t) (37)
式(37)中，Vs(U,t)为t时刻铣刀刀齿后刀面瞬时摩擦上边界曲线在U位置处V坐标值，Vs(U)为铣刀刀齿后刀面累积摩擦上边界曲线。
(2)提取不同时刻下刀齿后刀面瞬时摩擦下边界曲线，铣刀刀齿后刀面瞬时摩擦边界取相交特征点的最小值Vmin作为铣刀刀齿后刀面累积摩擦下边界，如图36所示。
Vp(U)＝minVξ(U,t) (38)
式(38)中，Vξ(U,t)为t时刻铣刀刀齿后刀面瞬时摩擦下边界曲线在U位置处V坐标值，Vp(U)为铣刀刀齿后刀面累积摩擦下边界曲线。
根据累积摩擦边界的解算方法求得在实验方案一下切削行程5m三个刀齿的解算结果如图37所示。
(3)为验证累积摩擦边界分形特征表征方法的正确性，通过式(13)、(14)、(35)、(36)、(37)、(38)，利用分形特征参数表征方法对实验与解算的摩擦边界分性特征参数，如式(39)、(40)所示。
D(t)＝{Ds(t),Dp(t)},G(t)＝{Gs(t),Gp(t)} (39)
D＝{Ds,Dp},G＝{Gs,Gp} (40)
其中，D(t)与G(t)分别为刀齿后刀面瞬时摩擦边界分形特征参数；Ds(t),Dp(t)分别为刀齿后刀面瞬时摩擦上下边界分形维数；Gs(t),Gp(t)分别为刀齿后刀面瞬时摩擦上下边界尺度系数。D、G分别为刀齿后刀面累积磨损边界分形特征参数解算结果，Ds、Gs分别为刀齿后刀面累积磨损上边界分形特征参数解算结果；Dp、Gp分别为刀齿后刀面累积磨损下边界分形特征参数解算结果。
(4)对解算的累积摩擦边界与实验累积摩擦边界的分形特征参数进行相对误差计算，如式(41)所示。
Figure BDA0003456527400000151
其中，D0、G0分别为刀齿后刀面累积摩擦边界分形特征参数实验结果。
表4刀齿后刀面累积摩擦边界分形特征参数解算结果
Figure BDA0003456527400000152
表5刀齿后刀面累积摩擦边界分形特征参数相对误差
Figure BDA0003456527400000153
由表5可知，铣刀各刀齿后刀面实验与解算累积摩擦上下边界曲线分形特征参数平均相对误差值均小于15％，其结果表明铣刀刀齿后刀面摩擦边界解算结果与实验的吻合度较高。综上所述，采用铣刀刀齿后刀面摩擦边界分形特征参数识别方法可以定量的揭示出热力耦合场下铣削振动与误差对刀齿后刀面摩擦边界的变化特性。
具体实施方式三：高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法；
为了验证高能效铣刀刀齿后刀面摩擦边界分形表征识别方法的有效性与适用性。设计并进行高能效铣刀采用与方案一相同的铣刀、工件、安装方式、切削方式和检测方法，保持切削深度和切削宽度不变，通过提高铣刀转速，改变刀齿误差分布，以改变铣削振动特性，相应减少每齿进给量以保持切削效率不变，进行刀齿后刀面摩擦边界分形特征参数的解算，其方案二如表6所示。
表6高能效铣刀铣削钛合金实验方案
Figure BDA0003456527400000154
Figure BDA0003456527400000161
依据表6和式(13)、式(14)、式(34)、式(35)、式(36)、式(37)，获得铣刀刀齿后刀面不同刀齿瞬时摩擦边界曲线、分形特征参数曲线与累积摩擦边界曲线，如图28-图45所示。
根据式(13)和(36)采用实验二的工艺方案，对解算的累积摩擦边界与实验累积摩擦边界的分形特征参数进行相对误差计算，其解算与实验的分形参数与其相对误差结果如表7与8所示。
表7刀齿后刀面累积摩擦边界分形特征参数解算结果
Figure BDA0003456527400000162
表8刀齿后刀面累积摩擦边界分形特征参数相对误差
Figure BDA0003456527400000163
由表8可知，铣刀各刀齿后刀面实验与解算累积摩擦上下边界曲线分形特征参数平均相对误差值均小于15％，其结果表明采用上述模型和方法，可有效识别出铣刀刀齿后刀面摩擦边界分形特征，为铣刀刀齿后刀面摩擦磨损状态评判提供依据。
需要说明的是，在以上实施例中，只要不矛盾的技术方案都能够进行排列组合，本领域技术人员能够根据排列组合的数学知识穷尽所有可能，因此本发明不再对排列组合后的技术方案进行一一说明，但应该理解为排列组合后的技术方案已经被本发明所公开。
本实施方式只是对本专利的示例性说明，并不限定它的保护范围，本领域技术人员还可以对其局部进行改变，只要没有超出本专利的精神实质，都在本专利的保护范围内。

Claims (10)

1.高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于，包括以下步骤：
步骤S1、构建高能效铣刀坐标系及其刀齿后刀面方程
对铣刀及其刀齿结构与后刀面方程进行求解，获得刀齿后刀面方程和刀齿切削刃方程；
步骤S2、构建高能效铣刀刀齿后刀面摩擦边界测量坐标系；
步骤S3、对刀齿后刀面摩擦边界分形特征参数进行解算；
在步骤S2构建的高能效铣刀刀齿后刀面摩擦边界测量坐标系下，截取刀齿后刀面投影面，对瞬时摩擦边界特征点进行提取，对提取变量特征参数进行解算；
步骤S4、表征刀齿后刀面与加工过渡表面瞬时几何接触关系；
揭示振动作用下铣刀刀齿后刀面瞬时几何摩擦边界的形成过程，表征不同位置处铣刀与工件的瞬时接触关系；
步骤S5、构建刀齿后刀面瞬时法向应力、切向应力、温度场判据；
步骤S6、构建刀齿后刀面瞬时摩擦边界特征点判据；
步骤S7、解算刀齿后刀面瞬时摩擦边界分形特征参数动态分布特性。
2.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：所述步骤S1中，构建高能效铣刀坐标系及其刀齿后刀面方程的具体方法是：
步骤S11、确定刀齿i坐标系oi-xiyizi与铣刀结构坐标系o0-x0y0z0的旋转矩阵Q1、Q2、平移矩阵M1分别为：
Figure FDA0003456527390000011
Figure FDA0003456527390000012
其中ri＝rmax-Δri，η为铣刀结构坐标系zi轴与z0轴之间的安装角，ri为刀齿i的回转半径，o0e1为o0与刀尖点e1的连线，
Figure FDA0003456527390000013
为铣刀结构坐标系y0轴与o0e1在x0o0y0平面的夹角；
步骤S12、确定刀齿坐标系下后刀面方程，在刀齿坐标系中，刀齿i的后刀面方程为：
Figure FDA0003456527390000014
其中，
0≤s≤s0,δ-0≤δ≤δ (4)
R＝R0-s0·tanα+s·tanα (5)
Figure FDA0003456527390000021
式(3)、(4)、(5)和(6)中，δ为B-B剖面下直线oze与ozoi在xioizi平面上的夹角，s为刀齿后刀面上任一点e到xi轴的距离，z0为切削刃中点外接圆中心到刀齿坐标系原点oi在xioizi平面上的距离，δ0为B-B剖面下刀尖到oz与ozoi在xioizi平面上的夹角，s0为刀齿切削刃中点em到xi轴的距离，s为刀齿后刀面上任一点e到xi轴的距离，γ为直线oie1与直线oiem在xioizi平面上夹角，R为B-B剖面下刀尖点在xioizi平面上的外切圆半径，R0为em点在xioizi平面上的外切圆半径，α为在yioizi平面上的后角；
由式(3)满足如下条件可获得刀齿i的切削刃方程：
R·sinδ＝R1·sinβ (7)
s＝R1cosβ-(RI-s0) (8)
s1≤s≤sc，-β0≤β≤β0 (9)
Figure FDA0003456527390000022
式(7)、(8)、(9)和(10)中，β为C视图下直线oi'e与直线oi'em在xioiyi平面上的夹角，β0为直线oi'e1与直线oi'em在xioiyi平面上的夹角，s1为刀尖点e1到xi轴的距离；z0为切削刃中点外接圆中心到刀齿坐标系原点oi在xioizi平面上的距离。
3.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S2中，构建后的高能效铣刀刀齿后刀面摩擦边界测量坐标系的旋转矩阵Q7、Q8、平移矩阵M4分别为：
Figure FDA0003456527390000023
Figure FDA0003456527390000024
。
4.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S3中，对刀齿后刀面摩擦边界分形特征参数进行解算的具体方法是：
步骤S31、在高能效铣刀刀齿后刀面摩擦边界测量坐标系下，做若干平行于V轴且间距为ΔU的平行线截取刀齿后刀面投影面；分别与原始刃、上边界、下边界交于一点，其中e1、e2、e3、e4分别为切削刃原始边界、摩擦上边界、摩擦下边界、摩擦区域上的点；在投影坐标系下，Ui0为刀齿后刀面摩擦边界U向最小值为Uie与Ve分别为切削刃中点e0的U，V轴的坐标，Uia与Uib分别为沿U轴方向的任意位置取样线的坐标，Uin与Vn分别为切削刃原始边界最大U向坐标值与刀尖点V向坐标值；
步骤S32、对瞬时摩擦边界特征点进行提取，将铣刀刀齿后刀面摩擦边界曲线U-V视作一个空间序列，对提取变量特征参数进行解算，具体解算公式为：
Figure FDA0003456527390000031
式中，S(γ)是结构函数的测度，即表示差方算数平均根[V(U0+γ)-V(U0)]2/γ，γ是采样间隔的尺度，是数据间隔的任意取值，S(ω)是功率谱密度函数；D是分形维数，G为尺度系数；
针对若干γ对铣刀刀齿后刀面摩擦边界曲线的离散信号，计算出相应的S(α)，在双对数坐标logS(γ)-logγ中，通过matlab运用最小二乘法对其进行拟合得到直线斜率k与截距A，换算得到表面轮廓曲线的分形维数D和尺度系数G，其中D与斜率k的转换关系为：
D＝2-k/2 (14)。
5.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S4中，表征刀齿后刀面与加工过渡表面瞬时几何接触关系的具体方法是：
步骤S41、确定铣刀结构坐标系o0-x0y0z0与振动作用下的铣刀切削坐标系od'-xd'yd'zd'的瞬时旋转矩阵Q3为：
Figure FDA0003456527390000032
其中，
Figure FDA0003456527390000033
为y0轴与yd'轴在xd'od'yd'平面内的瞬时夹角：
Figure FDA0003456527390000034
式中，
Figure FDA0003456527390000035
为铣刀初始切入t＝0时刻y0轴与yd'轴在xd'od'yd'平面内的夹角；
步骤S42、振动作用下的铣刀切削坐标系od'-xd'yd'zd'与无振动作用下的铣刀切削坐标系od-xdydzd的瞬时旋转矩阵Q4、Q5，瞬时平移矩阵M2分别为：
Figure FDA0003456527390000036
Figure FDA0003456527390000037
其中，θ1(t)为铣刀瞬时姿态角在xdodzd面的投影，θ2(t)为铣刀瞬时姿态角在ydodzd面的投影与xo轴的瞬时夹角，θ(t)为zd'轴与zd轴的瞬时夹角：
Figure FDA0003456527390000041
zd'轴与zd轴的瞬时夹角θ(t)为：
Figure FDA0003456527390000042
步骤S43、无振动作用下的铣刀切削坐标系od-xdydzd与工件坐标系o-xyz的瞬时平移矩阵M3为：
Figure FDA0003456527390000043
刀齿i的坐标系oi-xiyizi与工件坐标系o-xyz变换矩阵Φi为：
Φi＝M3·M2·Q5·Q4·Q3·M1·Q2·Q1 (22)
刀齿i的中心点oi在工件坐标系中的运动轨迹oi(x,y,z)为：
oi(x，y，z)＝[x y z 1]T＝Φi·[0 0 0 1]T (23)
步骤S44、根据刀齿切削刃在工件中瞬时运动关系，解算刀齿从切入到切出加工过渡表面，其刀齿与加工过渡表面的瞬时接触关系；
在工件坐标系中，刀齿i的切削刃方程I(x,y,z)为：
Figure FDA0003456527390000044
解算刀齿i瞬时切入至切出加工过渡表面方程Ki为：
Figure FDA0003456527390000045
其中，I(x,y,z)为刀齿i切削过程中，由t1 i至t2 i时段内，在工件坐标系下的切削刃方程，t1 i为刀齿t2 i切入待加工表面的初始时刻，为刀齿i切出待加工表面时刻；e为刀齿后刀面中任一点；
e(x1，y1，z1)点为铣削加工过程中刀齿后刀面上的任一点，做其后刀面所在的切平面的法矢量
Figure FDA0003456527390000046
与加工过渡表面相交与p(x2，y2，z2)点为：
Figure FDA0003456527390000051
在刀齿后刀面上选取点e(x1，y1，z1)，过e所在的切削平面向加工过渡表面方向做法矢量
Figure FDA0003456527390000052
与加工过渡表面交于点p(x2，y2，z2)，设
Figure FDA0003456527390000053
以刀齿后刀面的特征点指向刀齿外的方向为正方向，其大小为
Figure FDA0003456527390000054
其瞬时接触关系判定：
1)当
Figure FDA0003456527390000055
时，
Figure FDA0003456527390000056
为正方向，e与p相切，属于接触临界点；
2)当
Figure FDA0003456527390000057
时，
Figure FDA0003456527390000058
为负方向，e与p未产生接触；
3)当
Figure FDA0003456527390000059
时，
Figure FDA00034565273900000510
为正方向，e与p处产生摩擦，磨损。
6.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S5中构建刀齿后刀面瞬时法向应力、切向应力、温度场判据，具体方法是：
步骤S51、通过Deform仿真，对后刀面上特征点的有限元仿真结果进行提取；
步骤S52、在工件坐标系下，对仿真结果提取特征节点的等效应力，将其等效应力沿正四面体的三个方向力σq 1，σq 2，σq 3进行分解，并将等效应力三个分量在特征点e(x1，y1，z1)，沿着其切削平面的法矢量
Figure FDA00034565273900000511
的投影的方向，在刀齿坐标系下，刀齿后刀面上特征点e(x1，y1，z1)处等效应力三个分量与法向应力σ的矢量方向的矢量夹角：
Figure FDA00034565273900000512
其中，
Figure FDA00034565273900000513
为等效应力沿正四面体三个方向分力，
Figure FDA00034565273900000514
为等效应力三个方向分量与法向应力间的夹角，法向应力σ为：
Figure FDA00034565273900000515
式中
Figure FDA00034565273900000516
σj为等效应力三个方向分量在法向应力方向上的分量；
铣刀沿od(x,y,z)运动轨迹的切矢量x，y，z三个坐标方向上的运动速度的向量表达
Figure FDA00034565273900000517
为：
Figure FDA00034565273900000518
其中，vx(t)、vy(t)、vz(t)分别为后刀面任意点的运动速度大小vvn(t)在x、y、z三个坐标方向上的分量；
步骤S53、将摩擦速度通过矩阵变换转换到铣刀刀齿坐标系内为：
Figure FDA0003456527390000061
其中
Figure FDA0003456527390000062
为Φi的逆矩阵。
步骤S54、铣刀刀齿后刀面的任一点的切向应力τ方向为相对运动速度vvn方向在刀齿后刀面与加工过渡表面的公切面投影的方向
Figure FDA0003456527390000063
为：
Figure FDA0003456527390000064
在工件坐标系下，刀齿后刀面特征点e(x1，y1，z1)处等效应力三个分量与切向应力τ的矢量方向的矢量夹角如式(32)为：
Figure FDA0003456527390000065
其中，
Figure FDA0003456527390000066
为等效应力三个方向分量与切向应力τ间的夹角，切向应力τ为：
Figure FDA0003456527390000067
式中
Figure FDA0003456527390000068
τj为等效应力三个方向分量在切向应力方向上的分量。
7.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S6中、刀齿后刀面瞬时摩擦边界特征点判据为：
刀齿后刀面瞬时摩擦下边界识别出后刀面各截面位置几何特征、温度场、法向应力、切向应力的变化率突变的特征节点为：
Figure FDA0003456527390000069
其中，T、σ、τ分别温度、法向应力、切向应力，T(U,V)、σ(U,V)、τ(U,V)分别为刀齿后刀面同一时刻瞬时温度、法向应力、切向应力特征点数值；T'、σ'、τ'分别为齿后刀面瞬时温度、法向应力、切向应力变化率突变的特征点。
8.根据权利要求1所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S7中，刀齿后刀面瞬时摩擦边界的解算方法为：
步骤S71、在后刀面摩擦磨损测量坐标系中，截平面上特征节点的屈服强度所对应位置，为刀齿后刀面瞬时摩擦磨损上边界，识别出t时刻刀齿后刀面上所有截平面上等效应力大于屈服强度的临界特征节点作为刀齿后刀面瞬时摩擦上边界曲线：
Figure FDA0003456527390000071
其中，σq为刀齿后刀面等效应力特征点数值，σs为屈服强度的大小；Vs(U，t)为刀齿后刀面瞬时摩擦上边界，Vσs(U，t)为刀齿后刀面瞬时摩擦边界屈服强度曲线Vs(U，t-Δt)为上一时刻铣刀刀齿后刀面瞬时摩擦边界；
步骤S72、根据判据结果识别出铣刀刀齿t时刻后刀面不同位置处几何接触关系、法向应力、切向应力、温度变化率发生突变的节点位置，并在测量坐标系下分别表征出瞬时特征点的变化曲线，在测量坐标系下，识别出沿切削刃特征点切矢量方向的不同判据的最大值Vmax，作为刀齿后刀面瞬时摩擦下边界判据曲线
Vξ(U，t)＝max{VN(U，t)，Vσ(U，t)，VT(U，t)，Vτ(U，t)} (36)
其中，V0(U)为刀齿后刀面原始刃边界，VN(U，t)，Vσ(U，t)，Vτ(U，t)，VT(U，t)分别为t时刻刀齿后刀面上几何接触、法向应力、切向应力、温度判据识别曲线，Vξ(U，t)为t时刻铣刀刀齿后刀面瞬时摩擦边界；
步骤S73、在铣削行程中识别出铣刀三个刀齿切入到切出不同周期相同位置角的瞬时摩擦边界变化曲线，根据分形特征参数表征方法分别对其进行表征，识别出铣刀不同刀齿后刀面瞬时摩擦边界的分形特征参数的变化过程。
9.根据权利要求1-8任意一项所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：还包括步骤S8、刀齿后刀面累积摩擦边界解算及其分形特征参数的验证方法。
10.根据权利要求9所述的高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法，其特征在于：步骤S8中，刀齿后刀面累积摩擦边界解算及其分形特征参数的验证方法包括：
步骤S81、将刀齿i分别在铣刀切入、切中、切出分别取不同刻每个周期内瞬时摩擦最大摩擦边界形成的上边界曲线进行识别与表征，铣刀刀齿后刀面累积摩擦上边界为式(37)所示：
Vs(U)＝Vs(U，t) (37)
式中，Vs(U,t)为t时刻铣刀刀齿后刀面瞬时摩擦上边界曲线在U位置处V坐标值，Vs(U)为铣刀刀齿后刀面累积摩擦上边界曲线；
步骤S82、提取不同时刻下刀齿后刀面瞬时摩擦下边界曲线，铣刀刀齿后刀面瞬时摩擦边界取相交特征点的最小值Vmin作为铣刀刀齿后刀面累积摩擦下边界，具体表达式为：
Vp(U)＝minVξ(U，t) (38)
式中，Vξ(U，t)为t时刻铣刀刀齿后刀面瞬时摩擦下边界曲线在U位置处V坐标值，Vp(U)为铣刀刀齿后刀面累积摩擦下边界曲线；
步骤S83、为验证累积摩擦边界分形特征表征方法的正确性，通过式(13)、(14)、(35)、(36)、(37)、(38)，利用分形特征参数表征方法对实验与解算的摩擦边界分性特征参数，如式(39)、(40)所示：
D(t)＝{Ds(t)，Dp(t)}，G(t)＝{Gs(t)，Gp(t)} (39)
D＝{Ds，Dp}，G＝{Gs，Gp} (40)
其中，D(t)与G(t)分别为刀齿后刀面瞬时摩擦边界分形特征参数；Ds(t)，Dp(t)分别为刀齿后刀面瞬时摩擦上下边界分形维数；Gs(t)，Gp(t)分别为刀齿后刀面瞬时摩擦上下边界尺度系数。D、G分别为刀齿后刀面累积磨损边界分形特征参数解算结果，Ds、Gs分别为刀齿后刀面累积磨损上边界分形特征参数解算结果；Dp、Gp分别为刀齿后刀面累积磨损下边界分形特征参数解算结果；
步骤S84、对解算的累积摩擦边界与实验累积摩擦边界的分形特征参数进行相对误差计算，如式(41)所示：
Figure FDA0003456527390000081
其中，D0、G0分别为刀齿后刀面累积磨损边界分形特征参数实验结果。
CN202210005504.7A 2022-01-05 2022-01-05 高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法 Pending CN114429064A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210005504.7A CN114429064A (zh) 2022-01-05 2022-01-05 高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210005504.7A CN114429064A (zh) 2022-01-05 2022-01-05 高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法

Publications (1)

Publication Number Publication Date
CN114429064A true CN114429064A (zh) 2022-05-03

Family

ID=81311780

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210005504.7A Pending CN114429064A (zh) 2022-01-05 2022-01-05 高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法

Country Status (1)

Country Link
CN (1) CN114429064A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115647440A (zh) * 2022-11-10 2023-01-31 哈尔滨理工大学 方肩铣刀主副切削刃铣削微元能耗特征参数解算方法
CN116306169A (zh) * 2023-05-08 2023-06-23 哈尔滨理工大学 高效铣刀刀齿后刀面非稳态摩擦特性识别方法
CN117034725A (zh) * 2023-08-07 2023-11-10 哈尔滨理工大学 振动作用下铣刀后刀面摩擦力和热传导的热力学熵值解算方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115647440A (zh) * 2022-11-10 2023-01-31 哈尔滨理工大学 方肩铣刀主副切削刃铣削微元能耗特征参数解算方法
CN116306169A (zh) * 2023-05-08 2023-06-23 哈尔滨理工大学 高效铣刀刀齿后刀面非稳态摩擦特性识别方法
CN116306169B (zh) * 2023-05-08 2023-09-19 哈尔滨理工大学 高效铣刀刀齿后刀面非稳态摩擦特性识别方法
CN117034725A (zh) * 2023-08-07 2023-11-10 哈尔滨理工大学 振动作用下铣刀后刀面摩擦力和热传导的热力学熵值解算方法
CN117034725B (zh) * 2023-08-07 2024-04-12 哈尔滨理工大学 振动作用下铣刀后刀面摩擦力和热传导的热力学熵值解算方法

Similar Documents

Publication Publication Date Title
CN114429064A (zh) 高能效铣刀刀齿后刀面摩擦边界分形特征的识别方法
Altıntas et al. Mechanics and dynamics of ball end milling
Engin et al. Mechanics and dynamics of general milling cutters.: Part I: helical end mills
CN112372372B (zh) 一种高效铣刀累积摩擦磨损边界识别与验证方法
Abrari et al. On the dynamics of ball end milling: modeling of cutting forces and stability analysis
CN111002104B (zh) 一种高进给铣刀刀齿后刀面摩擦磨损边界检测与解算方法
Tsay et al. Accurate 5-axis machining of twisted ruled surfaces
Wei et al. Prediction of cutting force in ball-end milling of sculptured surface using improved Z-map
CN110032794B (zh) 一种振动作用下的铣刀动态切削力模型构建与验证方法
CN113962105B (zh) 一种无颤振精加工铣削过程的高效参数优化方法
CN109857061A (zh) 一种基于热力耦合模型的工件表面残余应力调控方法
CN116306169B (zh) 高效铣刀刀齿后刀面非稳态摩擦特性识别方法
Ma et al. Tool posture dependent chatter suppression in five-axis milling of thin-walled workpiece with ball-end cutter
Ghorbani et al. Extraction of surface curvatures from tool path data and prediction of cutting forces in the finish milling of sculptured surfaces
Scippa et al. Milled surface generation model for chip thickness detection in peripheral milling
Zhu et al. High efficiency simulation of five-axis cutting force based on the symbolically solvable cutting contact boundary model
CN108255134A (zh) 一种考虑倒棱几何的难加工材料高速车削切削力预测方法
Engin et al. Generalized modeling of milling mechanics and dynamics: Part I—helical end mills
CN109299581B (zh) 一种结合曲面插值的端面铣刀铣削力预测方法
Li et al. A numerical study of the effects of cutter runout on milling process geometry based on true tooth trajectory
Hilligardt et al. A new approach for local cutting force modeling enabling the transfer between different milling conditions and tool geometries
Quan et al. Cutting path planning for ruled surface impellers
CN106599389B (zh) 一种平头立铣刀宽行加工的刀刃接触判别方法
Reddy et al. A mechanistic force model for contour turning
CN116604401A (zh) 振动作用下方肩铣刀副后刀面非线性摩擦动力学识别方法

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