CN112464399A - 一种基于运动旋量的五轴加工刀具包络面计算方法 - Google Patents
一种基于运动旋量的五轴加工刀具包络面计算方法 Download PDFInfo
- Publication number
- CN112464399A CN112464399A CN202011260641.2A CN202011260641A CN112464399A CN 112464399 A CN112464399 A CN 112464399A CN 202011260641 A CN202011260641 A CN 202011260641A CN 112464399 A CN112464399 A CN 112464399A
- Authority
- CN
- China
- Prior art keywords
- coordinate system
- tool
- motion
- cutter
- feed
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- B23Q15/00—Automatic control or regulation of feed movement, cutting velocity or position of tool or work
- B23Q15/20—Automatic control or regulation of feed movement, cutting velocity or position of tool or work before or after the tool acts upon the workpiece
- B23Q15/22—Control or regulation of position of tool or workpiece
-
- 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
Abstract
本发明公开了一种基于运动旋量的五轴加工刀具包络面计算方法,属于铣削加工制造领域,包括:将刀具的空间运动分解为刀具以固定倾角沿刀路曲线所作的进给运动和刀具的倾角变化引起的姿态变动,利用Frenet标架对进给运动和姿态变动建模;基于进给运动和姿态变动的建模结果,对刀具的运动旋量进行建模,得到刀刃微元Q在刀具坐标系(TCS)下的运动旋量[VT];根据运动旋量[VT]计算刀刃微元Q在ti~ti+1时间段内做旋量运动的包络面,将其离散化为包络面点簇,并转换为工件坐标系(WCS)中刀具在ti~ti+1时间段内的包络面点簇;将WCS中刀具在各时间段内的包络面点簇集合起来,得到刀具的完整包络面。本发明在复杂曲面加工工况中也能精确生成刀具运动包络面。
Description
技术领域
本发明属于铣削加工制造领域,更具体地,涉及一种基于运动旋量的五轴加工刀具包络面计算方法。
背景技术
刀具运动包络面的精确求解是切削仿真理论的关键基础技术之一,可以将其应用于碰撞干涉检测、行距步距优化、刀具路径规划和工件表面更新等方面的工作。
为了对五轴加工刀具运动包络面进行求解,常用的方法主要是数值法和解析法。数值法主要包括扫掠微分方程法、隐式建模方法、包络理论法、雅可比降秩法等,这些方法需要求解复杂的微分方程,计算量很大。解析法可以快速求解刀具运动包络面,求解精度较高,但理论推导公式复杂,有学者采用双参数球族包络理论,可以快速解析得到侧铣加工刀具的运动包络面,但同样存在理论推导公式复杂的问题。
此外,在相关研究中刀具运动包络面的求解只针对特定的工况,对于刀路弯扭或刀具姿态变化程度较大工况下刀具运动包络面的求解还有待深入研究。
发明内容
针对现有技术的缺陷和改进需求,本发明提供了一种基于运动旋量的五轴加工刀具运动包络面计算方法,其目的在于,将运动旋量应用于描述五轴加工刀具的运动,并考虑刀路弯扭和刀具姿态变化,以解决现有的刀具运动包络面求解方法无法在刀路弯扭或刀具姿态变化程度较大的工况下精确计算刀具包络面的技术问题。
为实现上述目的,本发明提供了一种基于运动旋量的五轴加工刀具包络面计算方法,包括:
将刀具的空间运动分解为刀具以固定倾角沿刀路曲线所作的进给运动和刀具的倾角变化引起的姿态变动,利用Frenet标架对进给运动和姿态变动建模,得到进给坐标系中刀刃微元Q在ti时刻由于进给运动产生的第一瞬时速度和由于姿态变动产生的第二瞬时速度,将第一瞬时速度和第二瞬时速度之和作为进给坐标系中刀刃微元Q在运动旋量[VF]作用下的第三瞬时速度vF;
计算进给坐标系中刀刃微元Q在运动旋量[VF]作用下的角速度ωF,从而得到刀具在进给坐标系中的运动旋量[VF]为[VF]=(vF,ωF);将运动旋量[VF]变换至刀具坐标系下,得到刀刃微元Q在刀具坐标系下的运动旋量[VT];
在t∈[ti,ti+1]的情况下,计算刀具坐标系中刀刃微元Q在运动旋量[VT]的作用下经历时间t后所移动到的位置点坐标,以计算刀具在ti~ti+1时间段内做旋量运动的包络面,将其离散化,得到刀具坐标系中刀具在ti~ti+1时间段内的包络面点簇并将包络面点簇转换至工件坐标系中,得到工件坐标系统中刀具在ti~ti+1时间段内的包络面点簇Πi (WCS);将工件坐标系中刀具在各时间段内的包络面点簇集合起来,得到刀具在整条道路曲线上的完整包络面。
进一步地,进给坐标系中刀刃微元Q在ti时刻由于进给运动产生的第一瞬时速度为:
其中,vf为刀具的进给速率;Ki和τi分别表示刀路曲线在ti时刻的曲率和挠率;θi表示ZF(i)轴绕αi旋转至γi所历经的角度,αi和γi分别表示刀路曲线在OF(i)点的切矢和副法矢,OF(i)为ti时刻的刀触点,sθi和cθi分别表示sinθi和cosθi,表示的θi变化率;(xF,yF,zF,1)为刀刃微元Q在FCS(i)中的坐标,FCS(i)为OF(i)处的进给坐标系,ZF(i)表示进给坐标系FCS(i)的坐标轴。
进一步地,利用Frenet标架对进给运动建模,包括:
分别获得刀路曲线在刀触点OF(i)处的径矢ri、切矢αi、主法矢βi和副法矢γi,以及刀路曲线在刀触点OF(i+1)处的径矢ri+1、切矢αi+1、主法矢βi+1和副法矢γi+1,分别构成Frenet标架{ri;αi,βi,γi}和Frenet标架{ri+1;αi+1,βi+1,γi+1};OF(i)和OF(i+1)是刀路曲线上相邻的点,且分别是相邻的ti时刻和ti+1时刻的刀触点;构成Frenet标架的相关矢量满足如下Frenet公式:
其中,r、α、β、γ、K和τ分别表示同一刀触点处的径矢、切矢、主法矢、副法矢,以及刀路曲线的曲率和挠率;
利用Frenet标架{ri;αi,βi,γi}和Frenet标架{ri+1;αi+1,βi+1,γi+1}分别构建笛卡尔坐标系Fra(i)和Fra(i+1);
结合Frenet公式将刀路曲线在OF(i)点进行Taylor展开后表示为矩阵形式,并进行其次坐标变换,得到笛卡尔坐标系Fra(i+1)到Fra(i)的齐次坐标变换矩阵为:
其中,a11~a14、a21~a24、a31~a34均为含有Ki、τi和Δd的表达式,Δd表示刀刃微元Q从ti时刻到ti+1时刻的弧长;
在刀触点OF(i)和OF(i+1)处分别建立进给坐标系FCS(i)和FCS(i+1);XF(i)、YF(i)和ZF(i)分别表示进给坐标系FCS(i)的坐标轴,XF(i)轴与αi重合,ZF(i)轴表示设计表面在OF(i)点处的刀具侧单位法矢;XF(i+1)、YF(i+1)和ZF(i+1)分别表示进给坐标系FCS(i+1)的坐标轴,XF(i+1)轴与αi+1重合,ZF(i+1)轴表示设计表面在OF(i+1)点处的刀具侧单位法矢;
根据所建立的进给坐标系FCS(i)和FCS(i+1),计算进给坐标系从OF(i)点到OF(i+1)点的变化为:
进一步地,进给坐标系中刀刃微元Q在ti时刻由于姿态变动产生的第二瞬时速度为:
进一步地,利用Frenet标架对姿态变动建模,包括:
不考虑进给运动,将进给坐标系FCS(i)和FCS(i+1)统一记作FCS;
在进给坐标系中表示Qi点到Qi+1点的位移,并对时间求导,得到刀刃微元Q在Qi点的第二瞬时速度为:
进一步地,将运动旋量[VF]变换至刀具坐标系下,得到刀刃微元Q在刀具坐标系下的运动旋量[VT],包括:
其中,rT为刀具坐标系中的径矢,hT为刀具坐标系中的旋距,hF为进给坐标系中的旋距,rF为进给坐标系中的正交位置向量;表示进给坐标系到刀具坐标系的旋转变换矩阵,表示进给坐标系到刀具坐标系的齐次坐标变换矩阵。
进一步地,计算刀具坐标系中刀刃微元Q在运动旋量[VT]的作用下经历时间t后所移动到的位置点坐标,包括:
其中,rTx、rTy、rTz分别代表旋转运动轴线上一点在刀具坐标系下的x、y、z坐标;ωTx、ωTy、ωTz分别代表旋转运动轴线在刀具坐标系下的x、y、z方向矢量;Rx、Ry、Rz分别代表刀刃微元Q到旋转运动轴线上一点在刀具坐标系下x、y、z方向的距离位置矢量;c表示cos,s表示sin。
进一步地,计算刀具在ti~ti+1时间段内做旋量运动的包络面,包括:
进一步地,刀具坐标系到进给坐标系的旋转变换矩阵为:
进给坐标系到刀具坐标系的齐次坐标变换矩阵为:
进给坐标系到工件坐标系的其次坐标变换矩阵为:
其中,lead为前倾角,tilt为侧倾角,OF (TCS)表示刀触点OF在刀具坐标系中的坐标,XF、YF和ZF分别表示以刀触点为OF原点的进给坐标系的三个坐标轴,XF (WCS)、YF (WCS)和ZF (WCS)表示进给坐标系的对应坐标轴在工件坐标系中的坐标。
总体而言,通过本发明所构思的以上技术方案,能够取得以下有益效果:
(1)本发明将刀具空间运动分解为刀具以固定倾角沿刀路曲线所作的进给运动和刀具的倾角变化引起的姿态变动,并利用Frenet标架对进给运动和姿态变动建模,能够在计算刀具包络面的过程中考虑刀路弯扭和刀具姿态变化,从而在复杂曲面加工工况中也能精确生成刀具运动包络面;本发明在进给运动建模和姿态变动建模的基础上,解析地推导了能够精确描述刀具运动的运动旋量模型,由于旋量理论在描述刚体空间运动上具有简单清晰的表述方式且便于运动学求解,因此,本发明能够利用运动旋量在运动学描述方面的优势,精确计算刀具包络面。
(2)本发明所提供的基于运动旋量的五轴刀具包络面计算方法适用于包括球头刀、圆角刀等各种回转立铣刀,在五轴加工中碰撞干涉检测、行距步距优化、刀具路径规划和工件表面更新等方面具有较好的应用前景。
附图说明
图1为通用的立铣刀的几何模型;
图2为本发明实施例提供的刀具、工件相对位置示意图;
图3为本发明实施例提供的进给坐标系FCS到刀具坐标系TCS的旋转变换;
图4为本发明实施例提供的基于运动旋量的五轴加工刀具包络面计算方法流程图;
图5为本发明实施例提供的刀路曲线Γ上的坐标系;
图6为本发明实施例提供的刀具运动包络面生成示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
在本发明中,本发明及附图中的术语“第一”、“第二”等(如果存在)是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。
针对现有的刀具包络面求解方法无法在刀路弯扭或刀具变化程度较大的工况下精确计算刀具包络面的技术问题,本发明提供了一种基于运动旋量的五轴加工刀具包络面计算方法,其整体思路在于:将刀具空间运动分解为刀具的进给运动和姿态变动,根据曲线Frenet标架并通过齐次坐标变换推导分别对两种运动形式进行建模;结合旋量理论推导出描述刀具空间运动的运动旋量解析模型,实现对考虑刀路弯扭及刀具姿态变化的复杂曲面加工中刀具瞬时运动位姿的精确表示;最后根据加工刀路的运动旋量和刀具的几何参数,生成刀具运动包络面。
为了计算刀具包络面,需要使用到刀具的几何参数,作为一种可选的实施方式,本发明采用如图1所示的通用立铣刀几何模型对刀具进行几何建模。如图1所示,在刀具坐标系(TCS)中,通用立铣刀的几何形状可由七个独立参数完全确定:{H,Nz,Mz,α,β,Rr,r},其中刀具ZT轴的轴向高度H、Nz和Mz将刀具的包络面分为底刃区域、圆角刃区域和侧刃区域。在底刃区域,角度α为刀具下母线与XTOYT平面的夹角;在圆角刃区域,r为圆角半径,Rr为圆角圆心到刀具轴线的垂直距离;在侧刃区域,β为刀具上母线与ZT轴线的夹角。则刀具上任一刀刃微元点Q的坐标可表示为:
其中,z为Q点的轴向高度,r(z)为Q点到ZT轴的垂直距离,即为刀具的有效切削半径,它在刀具的三个区域可表示为:
其中,是指刀具的第jc条切削刃在轴向高度z处的径向接触角,即刀刃微元点Q到刀具ZT轴的垂线与刀具YT轴的夹角;κ是指刀具的m条切削刃在轴向高度z处的轴向接触角,即刀具扫掠面在刀刃微元点Q处的单位向外法向矢量n与刀具ZT轴的夹角,可分别表示为:
其中,jc表示Q点所在的刀刃序号,ψ为刀具旋转角,ψ(z)为径向滞后角。
刀具包络面在Q点的单位向外法向矢量n表示在刀具坐标系TCS中为:
为了在计算刀具包络面的过程中考虑刀路弯扭和刀具的姿态变化,本发明在实际计算过程中,涉及到不同坐标系之间的转换,坐标系之间的转换关系可按照如下方式进行计算:
图2所示为刀具与工件的相对位置示意图,其中D S表示设计表面,OW-XW-YW-ZW为工件坐标系(WCS)。OW-XW-YW-ZW为进给坐标系(RCS),其中OF点是刀具与DS面的相切点,即刀触点,ZF轴为DS在OF点的刀具侧幺法矢,XF轴为沿着瞬时进给速度方向的幺矢。
刀具轴线在FCS中的方向通常以前倾角lead和侧倾角tilt来描述,FCS到TCS的旋转变换如图3所示:FCS绕YF轴旋转,旋转角为lead,得到过渡坐标系O'-X'-Y'-Z',然后过渡坐标系绕XF轴旋转,旋转角为tilt,得到刀具坐标系TCS。则TCS到FCS的旋转变换矩阵为:
FCS到TCS的齐次坐标变换矩阵可表示为:
其中,OF (TCS)表示刀触点OF在TCS中的坐标,可通过OF点处的曲面法矢和刀具自旋面外法矢反向平行这一关系来求解;
FCS到WCS的齐次坐标变换矩阵可表示为:
其中,XF (WCS)、YF (WCS)和ZF (WCS)表示FCS的对应坐标轴在WCS中的坐标,OF (WCS)表示FCS的原点在WCS中的坐标。
基于以上刀具建模,以及坐标系之间的转换关系,本发明实施例提供的基于运动旋量的五轴加工刀具包络面计算方法,如图4所示,包括:
将刀具的空间运动分解为刀具以固定倾角沿刀路曲线所作的进给运动和刀具的倾角变化引起的姿态变动,利用Frenet标架对进给运动和姿态变动建模,得到进给坐标系中刀刃微元Q在ti时刻由于进给运动产生的第一瞬时速度和由于姿态变动产生的第二瞬时速度,将第一瞬时速度和第二瞬时速度之和作为进给坐标系中刀刃微元Q在运动旋量[VF]作用下的第三瞬时速度vF;
计算进给坐标系中刀刃微元Q在运动旋量[VF]作用下的角速度ωF,从而得到刀具在进给坐标系中的运动旋量[VF]为[VF]=(vF,ωF);将运动旋量[VF]变换至刀具坐标系下,得到刀刃微元Q在刀具坐标系下的运动旋量[VT];
在t∈[ti,ti+1]的情况下,计算刀具坐标系中刀刃微元Q在运动旋量[VT]的作用下经历时间t后所移动到的位置点坐标,以计算刀具在ti~ti+1时间段内做旋量运动的包络面,将其离散化,得到刀具坐标系中刀具在ti~ti+1时间段内的包络面点簇并将包络面点簇转换至工件坐标系中,得到工件坐标系统中刀具在ti~ti+1时间段内的包络面点簇Πi (WCS);将工件坐标系中刀具在各时间段内的包络面点簇集合起来,得到刀具在整条道路曲线上的完整包络面。
图5为刀路曲线Γ上各个坐标系的示意图,此曲线依附于DS上,可看作由不同时刻的刀触点OF组成。Γ的自然参数方程为r=r(d),OF(i)和OF(i)是Γ上邻近的点,且分别是相邻的ti时刻和ti+1时刻的刀触点;
本实施例中,按照如下方式完成对进给运动的建模,在此过程中考虑了刀路弯扭以及工件形状对进给运动的影响:
Γ在OF(i)点的Frenet标架表示为{ri;αi,βi,γi},其中的四个矢量分别表示Γ在OF(i)点的径矢、切矢、主法矢和副法矢;将{ri;αi,βi,γi}构成的笛卡尔坐标系记作Fra(i);同样地,Γ在OF(i+1)点的Frenet标架表示为{ri+1;αi+1,βi+1,γi+1},其中的四个矢量分别表示Γ在OF(i+1)点的径矢、切矢、主法矢和副法矢;将{ri+1;αi+1,βi+1,γi+1}构成的笛卡尔坐标系记作Fra(i+1);构成Frenet标架的相关矢量满足如下Frenet公式:
其中,r、α、β、γ、K和τ分别表示同一刀触点处的径矢、切矢、主法矢、副法矢,以及刀路曲线的曲率和挠率;
OF(i)-XF(i)-YF(i)-ZF(i)为在OF(i)点处的进给坐标系FCS(i),应当注意的是,XF(i)轴与αi重合,ZF(i)轴表示DS在OF(i)点处的刀具侧单位法矢,θi表示ZF(i)轴绕αi旋转至γi所历经的角度;在OF(i+1)点处的描述方法类似,OF(i+1)-XF(i+1)-YF(i+1)-ZF(i+1)为在OF(i+1)点处的进给坐标系FCS(i+1),应当注意的是,XF(i+1)轴与αi+1重合,ZF(i+1)轴表示DS在OF(i+1)点处的刀具侧单位法矢,θi+1表示ZF(i+1)轴绕αi+1旋转至γi+1所历经的角度。
因此,从OF(i)点到OF(i+1)点,FCS的变化可以描述为:
其中:
这里以sθ表示sinθ,cθ表示cosθ,下文中以此类推。
将刀路曲线Γ在OF(i)点处Taylor展开并结合Frenet公式,可得:
将上式改写成矩阵形式为:
其中,a11~a14、a21~a24、a31~a34均为含有Ki、τi和Δd的表达式,Δd表示刀刃微元Q从ti时刻到ti+1时刻的弧长;
利用齐次变换矩阵表示为:
由此便获得了为了将刀具进给运动以运动旋量的方式表示,还需要分析其对刀刃微元瞬时速度的影响。刀刃微元Q在ti时的位置记为Qi,Qi在FCS(i)中的坐标(xF,yF,zF,1),在FCS运动至FCS(i+1)后,Qi点移动到了Qi+1点,这两点满足如下表达式:
设FCS运动到FCS(i)和FCS(i+1)的时刻分别为ti和ti+1,且ti+1=ti+Δt,则Qi点在ti时刻的瞬时速度,即第一瞬时速度可表示为:
其中,vf为刀具进给速率,满足Δd=vfΔt。结合上述的推导,可化为如下形式:
本发明中,带有下标i的符号皆代表对应量在ti时刻的取值。
本实施例中,按照如下方式完成对姿态变动的建模:
刀具姿态变动建模不考虑进给运动,即将FCS固定,将FCS(i)和FCS(i)统一记作FCS,考虑TCS从TCS(i)运动到TCS(i+1)时刀刃微元的运动。同样以刀刃微元Qi作为研究对象,Qi在FCS中的坐标(xF,yF,zF,1),在刀具姿态改变后,Qi移动到了Qi+1,且它们在各自的刀具坐标系中的坐标保持不变,即:将Qi到Qi+1的位移[Δx2 Δy2 Δz2 0]T(FCS)表示在FCS中:
可进一步化为如下形式:
其中:
本实施例中,完后对刀具进给运动和姿态变动的建模后,按照如下方式对刀具的运动旋量进行建模:
将FCS中表征刀具空间运动的运动旋量记为[VF]=(vF,ωF),则有vF=rF×ωF+hFωF,其中:ωF为角速度,设为ωF=[ωFx ωFy ωFz]T,hF为旋距,rF为正交位置向量,设为rF=[rFx rFy rFz]T,满足ωF⊥rF,即:
ωFxrFx+ωFyrFy+ωFzrFz=0 (3-18)
则Qi(xF,yF,zF)点在[VF]作用下的瞬时速度表示为:
则[VF]可结合下式解出:
将TCS中表征刀具运动的运动旋量记为[VT]=(vT,ωT),其径矢为rT,旋距为hT,则有:
本实施例中,完成对刀具运动旋量的建模后,基于该建模结果,按照如下方式计算刀具包络面:
根据在TCS中的运动旋量[VT](如无特殊标记,以下所有矢量和坐标是均以TCS为参考),设任一刀刃微元Q在[VT]作用下历经时间t移动到Qs点,则Qs的坐标可表示为:
其中,表示矢量RP绕ω1旋转θ角度,ω1=[ω1xω1yω1z]T=ωT/|ωT|,RP=[Rx Ry Rz]T=Q-rT,θ=|ωT|/t。则可将Qs的坐标转换为与所述刀具的轴向高度z、径向接触角以及时间t相关的表达式如下:
其中,rTx、rTy、rTz分别代表旋转运动轴线上一点在刀具坐标系下的x、y、z坐标;ωTx、ωTy、ωTz分别代表旋转运动轴线在刀具坐标系下的x、y、z方向矢量;Rx、Ry、Rz分别代表刀刃微元Q到旋转运动轴线上一点在刀具坐标系下x、y、z方向的距离位置矢量;c表示cos,s表示sin;
其中,表示在ti时刻进给坐标系到刀具坐标系的齐次坐标变换矩阵,表示在ti时刻进给坐标系到工件坐标系的齐次坐标变换矩阵。将一条刀路中各个时间段的包络面的点簇集合起来:Π(WCS)=ΣiΠi (WCS),即可得到在整条刀路上刀具的包络面Π(WCS),刀具运动包络面生成如图6所示。
总体而言,本实施例针对考虑刀路弯扭及刀具姿态变化的通用立铣刀五轴加工,解析地推导了能够精确描述刀具运动的运动旋量模型,并在此基础上给出了刀具运动包络面的精确计算方法,在刀路弯扭和刀具姿态变动程度较大的工况下,本实施例同样可以精确地计算刀具包络面。并且,本实施例无需求解复杂的高阶常微分微分方程或超越方程,相比于传统的扫掠微分方程法、隐式建模方法、包络理论法、雅可比降秩法等常规的数值方法,大大减小了计算量。本实施例适用于包括球头刀、圆角刀等各种回转立铣刀,在五轴加工中碰撞干涉检测、行距步距优化、刀具路径规划和工件表面更新等方面具有较好的应用前景。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种基于运动旋量的五轴加工刀具包络面计算方法,其特征在于,包括:
将所述刀具的空间运动分解为刀具以固定倾角沿刀路曲线所作的进给运动和刀具的倾角变化引起的姿态变动,利用Frenet标架对所述进给运动和所述姿态变动建模,得到进给坐标系中刀刃微元Q在ti时刻由于进给运动产生的第一瞬时速度和由于姿态变动产生的第二瞬时速度,将所述第一瞬时速度和所述第二瞬时速度之和作为进给坐标系中刀刃微元Q在运动旋量[VF]作用下的第三瞬时速度vF;
计算进给坐标系中所述刀刃微元Q在所述运动旋量[VF]作用下的角速度ωF,从而得到所述刀具在进给坐标系中的运动旋量[VF]为[VF]=(vF,ωF);将所述运动旋量[VF]变换至刀具坐标系下,得到所述刀刃微元Q在刀具坐标系下的运动旋量[VT];
3.如权利要求2所述的基于运动旋量的五轴加工刀具包络面计算方法,其特征在于,利用Frenet标架对所述进给运动建模,包括:
分别获得所述刀路曲线在所述刀触点OF(i)处的径矢ri、切矢αi、主法矢βi和副法矢γi,以及所述刀路曲线在刀触点OF(i+1)处的径矢ri+1、切矢αi+1、主法矢βi+1和副法矢γi+1,分别构成Frenet标架{ri;αi,βi,γi}和Frenet标架{ri+1;αi+1,βi+1,γi+1};OF(i)和OF(i+1)是刀路曲线上相邻的点,且分别是相邻的ti时刻和ti+1时刻的刀触点;构成Frenet标架的相关矢量满足如下Frenet公式:
其中,r、α、β、γ、K和τ分别表示同一刀触点处的径矢、切矢、主法矢、副法矢,以及刀路曲线的曲率和挠率;
利用Frenet标架{ri;αi,βi,γi}和Frenet标架{ri+1;αi+1,βi+1,γi+1}分别构建笛卡尔坐标系Fra(i)和Fra(i+1);
结合Frenet公式将刀路曲线在OF(i)点进行Taylor展开后表示为矩阵形式,并进行其次坐标变换,得到笛卡尔坐标系Fra(i+1)到Fra(i)的齐次坐标变换矩阵为:
其中,a11~a14、a21~a24、a31~a34均为含有Ki、τi和Δd的表达式,Δd表示刀刃微元Q从ti时刻到ti+1时刻的弧长;
在刀触点OF(i)和OF(i+1)处分别建立进给坐标系FCS(i)和FCS(i+1);XF(i)、YF(i)和ZF(i)分别表示进给坐标系FCS(i)的坐标轴,XF(i)轴与αi重合,ZF(i)轴表示设计表面在OF(i)点处的刀具侧单位法矢;XF(i+1)、YF(i+1)和ZF(i+1)分别表示进给坐标系FCS(i+1)的坐标轴,XF(i+1)轴与αi+1重合,ZF(i+1)轴表示设计表面在OF(i+1)点处的刀具侧单位法矢;
根据所建立的进给坐标系FCS(i)和FCS(i+1),计算进给坐标系从OF(i)点到OF(i+1)点的变化为:
6.如权利要求1-5任一项所述的基于运动旋量的五轴加工刀具包络面的计算方法,其特征在于,将所述运动旋量[VF]变换至刀具坐标系下,得到所述刀刃微元Q在刀具坐标系下的运动旋量[VT],包括:
7.如权利要求6所述的基于运动旋量的五轴加工刀具包络面计算方法,其特征在于,计算刀具坐标系中刀刃微元Q在所述运动旋量[VT]的作用下经历时间t后所移动到的位置点坐标,包括:
其中,rTx、rTy、rTz分别代表旋转运动轴线上一点在刀具坐标系下的x、y、z坐标;ωTx、ωTy、ωTz分别代表旋转运动轴线在刀具坐标系下的x、y、z方向矢量;Rx、Ry、Rz分别代表刀刃微元Q到旋转运动轴线上一点在刀具坐标系下x、y、z方向的距离位置矢量;c表示cos,s表示sin。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011260641.2A CN112464399B (zh) | 2020-11-12 | 2020-11-12 | 一种基于运动旋量的五轴加工刀具包络面计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011260641.2A CN112464399B (zh) | 2020-11-12 | 2020-11-12 | 一种基于运动旋量的五轴加工刀具包络面计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464399A true CN112464399A (zh) | 2021-03-09 |
CN112464399B CN112464399B (zh) | 2022-04-22 |
Family
ID=74826642
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011260641.2A Active CN112464399B (zh) | 2020-11-12 | 2020-11-12 | 一种基于运动旋量的五轴加工刀具包络面计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464399B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113695645A (zh) * | 2021-09-03 | 2021-11-26 | 中国工程物理研究院激光聚变研究中心 | 一种微径金刚石球头铣刀竖直偏置校正方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001001020A1 (en) * | 1999-06-28 | 2001-01-04 | Visteon Global Technologies, Inc. | Helical and spur gear drive with double crowned pinion tooth surfaces and conjugated gear tooth surfaces |
JP2007286858A (ja) * | 2006-04-17 | 2007-11-01 | Toyota Motor Corp | 面モデルの作成装置と作成方法 |
CN104866655A (zh) * | 2015-05-07 | 2015-08-26 | 北京航空航天大学 | 一种基于包络理论和经线划分的回转型刀具包络特征线计算方法 |
CN106529045A (zh) * | 2016-11-15 | 2017-03-22 | 西北工业大学 | 基于旋量的多轴铣削加工刀具轴线建模方法 |
CN109062140A (zh) * | 2018-10-17 | 2018-12-21 | 清华大学 | 一种五轴侧铣加工过程中瞬时未变形切屑厚度计算方法 |
CN111061218A (zh) * | 2019-12-31 | 2020-04-24 | 华中科技大学 | 一种基于acis的复杂曲面加工切削力预测方法及设备 |
-
2020
- 2020-11-12 CN CN202011260641.2A patent/CN112464399B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2001001020A1 (en) * | 1999-06-28 | 2001-01-04 | Visteon Global Technologies, Inc. | Helical and spur gear drive with double crowned pinion tooth surfaces and conjugated gear tooth surfaces |
JP2007286858A (ja) * | 2006-04-17 | 2007-11-01 | Toyota Motor Corp | 面モデルの作成装置と作成方法 |
CN104866655A (zh) * | 2015-05-07 | 2015-08-26 | 北京航空航天大学 | 一种基于包络理论和经线划分的回转型刀具包络特征线计算方法 |
CN106529045A (zh) * | 2016-11-15 | 2017-03-22 | 西北工业大学 | 基于旋量的多轴铣削加工刀具轴线建模方法 |
CN109062140A (zh) * | 2018-10-17 | 2018-12-21 | 清华大学 | 一种五轴侧铣加工过程中瞬时未变形切屑厚度计算方法 |
CN111061218A (zh) * | 2019-12-31 | 2020-04-24 | 华中科技大学 | 一种基于acis的复杂曲面加工切削力预测方法及设备 |
Non-Patent Citations (1)
Title |
---|
王延忠等: "基于通用五坐标数控机床螺旋锥齿轮NC加工研究", 《中国机械工程》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113695645A (zh) * | 2021-09-03 | 2021-11-26 | 中国工程物理研究院激光聚变研究中心 | 一种微径金刚石球头铣刀竖直偏置校正方法 |
CN113695645B (zh) * | 2021-09-03 | 2022-07-01 | 中国工程物理研究院激光聚变研究中心 | 一种微径金刚石球头铣刀竖直偏置校正方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112464399B (zh) | 2022-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Li et al. | Mechanistic modeling of five-axis machining with a general end mill considering cutter runout | |
Zhu et al. | Global optimization of tool path for five-axis flank milling with a conical cutter | |
CN105425727B (zh) | 五轴侧铣加工刀具路径光顺方法 | |
Mali et al. | A comprehensive review of free-form surface milling–Advances over a decade | |
CN107577882B (zh) | 一种侧铣直纹曲面的表面形貌建模及成型的仿真方法 | |
Yan et al. | Multi-axis variable depth-of-cut machining of thin-walled workpieces based on the workpiece deflection constraint | |
CN103744349A (zh) | 一种平头立铣刀加工过渡曲面的无干涉刀具路径生成方法 | |
Xu et al. | Machine based energy-saving tool path generation for five-axis end milling of freeform surfaces | |
Warkentin et al. | Computer aided 5-axis machining | |
Rufeng et al. | Dual drive curve tool path planning method for 5-axis NC machining of sculptured surfaces | |
CN109343468B (zh) | 一种基于投影偏置的叶片多轴轨迹生成方法 | |
CN102866671A (zh) | 大圆弧直纹面数控加工刀具轨迹规划方法 | |
CN110032140B (zh) | 一种五轴加工中球形刀刀轴矢量规划方法 | |
CN101870073B (zh) | 基于工艺系统刚度特性的多轴数控加工刀具运动规划方法 | |
CN106202822A (zh) | 面向叶片自适应加工的b样条曲面模型重构方法 | |
CN109597357B (zh) | 一种面向叶片旋铣工艺的数控编程方法及装置 | |
CN109375579A (zh) | 一种基于运动学的五轴数控加工刀具姿态规划方法 | |
CN109570591A (zh) | 离心叶轮切削加工方法及装置和离心叶轮加工设备 | |
CN112464399B (zh) | 一种基于运动旋量的五轴加工刀具包络面计算方法 | |
Zhu et al. | High efficiency simulation of five-axis cutting force based on the symbolically solvable cutting contact boundary model | |
Chang et al. | Five-axis automated measurement by coordinate measuring machine | |
CN110989490B (zh) | 一种基于轮廓误差的工件最优安装位置的获取方法 | |
CN103324140A (zh) | 五轴加工中通用刀具扫描体的生成方法 | |
Sun et al. | A review on theories/methods to obtain surface topography and analysis of corresponding affecting factors in the milling process | |
Chen et al. | Modeling and simulation of the surface topography in ball-end milling based on biharmonic spline interpolation |
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 |