CN102222138B - 一种基于曲面截型线分划的最短距离线对获取方法 - Google Patents
一种基于曲面截型线分划的最短距离线对获取方法 Download PDFInfo
- Publication number
- CN102222138B CN102222138B CN 201110144357 CN201110144357A CN102222138B CN 102222138 B CN102222138 B CN 102222138B CN 201110144357 CN201110144357 CN 201110144357 CN 201110144357 A CN201110144357 A CN 201110144357A CN 102222138 B CN102222138 B CN 102222138B
- Authority
- CN
- China
- Prior art keywords
- cutter
- line
- curved surface
- point
- molded line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Landscapes
- Numerical Control (AREA)
Abstract
一种基于曲面截型线分划的最短距离线对获取方法,该方法有三大步骤:步骤一:曲面截型线离散;步骤二:截型线到刀具曲面最短距离点的获取;步骤三:最短距离线对的获取。本发明是将加工曲面沿进给方向离散成若干条截型线,根据包络原理在加工过程中的任意刀位得到每条截型线与刀具间的最短距离来获得该刀位的最短距离线对的方法,本方法构思科学,计算精确,实施容易,效率高,在曲面数控加工技术领域里具有较好的实用价值和广阔的应用前景。
Description
技术领域
本发明涉及一种数控加工中最短距离线对的获取方法,尤其涉及一种基于曲面截型线分划的最短距离线对获取方法,具体地说,是指将加工曲面沿进给方向离散成若干条截型线,根据包络原理在加工过程中的任意刀位求出每条截型线与刀具间的最短距离来获得该刀位的最短距离线对的方法,属于曲面数控加工技术领域。
背景技术
在五坐标数控加工中,通过优化刀位以增大行宽已经成为目前的研究热点问题,已有的刀位优化算法中最典型的有刀轴倾斜法和曲率匹配算法,但它们都是利用微分几何原理研究刀具和工件在一点邻域的接触状况,不足以描述二者在邻域外的相对位置关系。因此,利用离散曲面或刀具的方法研究二者的接触关系成为必然的发展趋势。
20世纪末,北京航空航天大学的张洪等提出并利用“最短距离线对”原理以直接建立任意刀位下的理论加工误差分布函数,即沿行宽方向加工误差的分布规律。随后又发展出一系列最短距离线对的获取方法,主要有离散刀具的环心圆法、经线法、纬线法和离散工件的点阵法。前三种方法的特点均是将刀具曲面按照不同的方法进行离散,然后通过求取刀具离散点与工件曲面间的最短距离来获得刀位误差,刀具离散点与工件曲面间的最短距离求取采用西交大提出的数值迭代算法或UG软件内置算法。其中数值迭代算法计算效率较高,但是迭代条件在工件曲面边界区域时不充分,计算不稳定。UG软件内置算法计算稳定,但是计算速度很慢、效率较低。点阵法的特点是将工件曲面离散成若干截型线,再将截型线离散成若干点,然后利用刀具曲面的回转特性直接用解析法求取工件曲面离散点与刀具曲面间的最短距离,进而获得最短距离线对。这样,曲线与自由曲面间的最短距离求解问题转换为点与回转曲面间的最短距离求解问题,求解稳定可靠。但是该方法在获得截型线与刀具曲面间的最短距离时,对截型线如何离散、离散成多少点等问题不好把握,当离散点数较少时,求得的截型线与刀具曲面间的最短距离误差较大,当离散点数较多时,计算时间也会随之增加,计算效率大大降低。
发明内容
1、目的:本发明的目的是提供一种基于曲面截型线分划的最短距离线对获取方法,它克服了现有技术的不足,从而解决数控加工中刀位优化计算的问题。
2、技术方案:
1)本发明一种基于曲面截型线分划的最短距离线对获取方法,它通过获取加工过程中任意刀位的最短距离线对来得到该刀位的误差分布,从而进行刀位优化。该方法适用的加工刀具为圆弧母线的回转刀具,如圆环面刀,鼓形刀等。
在曲面加工过程中存在三张曲面:工件设计曲面、刀具工作面和被刀具包络出来的已加工表面。而已加工表面在工件设计曲面的一定公差范围之内,所以数控加工中可以根据这个约束条件来研究工件设计曲面和刀具曲面的关系。在每个刀位,刀具曲面和已加工表面存在一条切线,称为刀具包络特征线,刀具包络特征线也是刀具曲面上的最短距离线。从刀具包络特征线上任意点往工件设计曲面求最短距离,在设计曲面上可以得到一个垂足点,将所有垂足点连成线便得到设计曲面上的最短距离线。刀具曲面上的最短距离线和设计曲面上的最短距离线构成最短距离线对,如图1所示。可以看出,最短距离线对是用来描述刀具曲面与工件曲面在相对运动过程中相互接近、接触或干涉的概念,它是两者沿进给方向贴合程度的反映。
2)本发明一种基于曲面截型线分划的最短距离线对获取方法,该方法的具体步骤如下:
步骤一:曲面截型线离散
将工件曲面沿进给方向离散成n条截型线,记工件曲面截型线的集合可表示为{Mi|i=1,2,...,n}。任意一条截型线可用参数方程表示为
Mi=ρ(ui,v) (1)
其中,ui可表示为第i条截型线的u值,为某一定值。
步骤二:截型线到刀具曲面最短距离点的获取
当求截型线i上到刀具曲面的最短距离时,
如图2所示,在刀具加工过程中,工件曲面截型线与刀具曲面间的最短距离线段对应着截型线和刀具曲面上各一个相应点,分别称为截型线上的最短距离点和刀具曲面上的最短距离点;刀具曲面上的最短距离点必然在刀具曲面最短距离线上,而截型线上的最短距离点也一定在工件曲面最短距离线上;记刀具曲面最短距离线上的最短距离点的集合为{Pi T|i=1,2,...},工件曲面最短距离线上的最短距离点的集合为{Pi M|i=1,2,…},点集内的点Pi T与Pi M对应着工件曲面第i条截型线Mi到刀具曲面的最短距离。取nPi T为刀具在点Pi T处的单位内法矢,τPi M为截型线在点Pi M处沿进给方向的单位切线矢量,则有
nPi T·τPi M≈0 (2)
因此,可知工件曲面截型线到刀具曲面的最短距离点应满足式(2);反之,可由式(2)求出工件曲面截型线上到刀具曲面距离最短的点,即获得各截型线到刀具曲面的最短距离点;
步骤三:最短距离线对的获取
根据求出的截型线Mi上到刀具的最短距离点,记为Pi M。过点Pi M与刀具回转轴的截平面截刀具与工件曲面,可得截平面内的两条刀具母线,在截平面内连接点Pi M与靠近工件曲面的刀具母线圆弧的圆心Oi T,直线Pi MOi T与截平面内靠近工件曲面的刀具母线交于点Pi T,点Pi T即为工件曲面截型线Mi上的最短距离点Pi M到刀具曲面的最短距离点,|Pi MPi T|的长度值即为截型线Mi与刀具的最短距离值。依此,可求出所有截型线到刀具曲面的最短距离值及对应刀具曲面上的最短距离点。分别将所有刀具曲面上的最短距离点和工件曲面截型线上的最短距离点依序连接成线,便获得该刀位下刀具与工件曲面间的最短距离线对,根据最短距离线对间的有向长度便可获得该刀位的误差分布。
其中,步骤二中所述的nPi T·τPi M≈0(2),是根据求解精度来控制,即给定控制精度为ξ时,当|nPi T·τPi M|≤ξ时可获得所求解。
其中,步骤二中所述的截型线到刀具曲面最短距离点的求取,是指当求截型线i上到刀具曲面的最短距离时,截型线i上的初始点跟据当前刀位的切触点来定,若当前刀位切触点为Ps(us,vs),则截型线i上选取的初始点为Pi0 M(u,v),其中初始点的u值与截型线i的u值一致,初始点的v值与切触点的v值一致,即u=ui,v=vs。
3、优点及功效:本发明一种基于曲面截型线分划的最短距离线对获取方法:
1)可根据包络原理求出每条截型线到刀具的最短距离点,获得最短距离线对。
2)根据最短距离线对可获得刀位误差分布,进而进行刀位优化调整。
3)通过实例计算,该方法的计算效率是原点阵法的15倍以上。
附图说明
图1为加工过程中最短距离线对示意图
图2为工件曲面截型线到刀具曲面距离最短点的获取示意图
图3为截型线与刀具曲面间距离最短点获取的具体实施示意图
图4为截型线到刀具曲面最短距离点的计算流程图
图5为实施例中计算曲面示意图
图6为实施例中某刀位最短距离线对计算结果的二维分布图
图中符号说明如下:
图2:S(u,v)表示工件曲面,u,v分别表示曲面参数;Pi T为刀具特征线上的最短距离点,i表示序号;Pi M为工件曲面特征线上的最短距离点,i表示序号;nPi T为刀具曲面在点Pi T处的单位法矢;vfPi T为刀具在点Pi T处的单位进给速度矢量;τPi M为截型线在点Pi M处的单位切线矢量。
图3:S(u,v)表示工件曲面,u,v分别表示曲面参数;E为截平面;Oi T为截平面内靠近工件曲面的刀具母线圆弧的圆心;R为刀具母线圆弧半径;Ps为该刀位在工件曲面的切触点;Pi0 M为第i条截型线上选取的初始点,该点的u值与截型线i的u值一致,v值与切触点的v值一致;Pi0 T为工件曲面上点Pi0 M到刀具曲面的最短距离点;nPi0 T为刀具曲面在点Pi0 T处的单位法矢;vfPi0 T为刀具在点Pi0 T处的单位进给速度;τPi0 M为工件曲面在点Pi0 M处沿参数v方向的单位切矢。
图4:us,vs为切触点Ps的u,v值;ui为截型线i的u值;ξ为求解控制精度;v0为截型线i的初始值;Δv为迭代搜索步长;v为截型线的当前v值;vmin为截型线的最小v值,vmax为截型线的最大v值;vi为截型线i上所求的v值。
图5:x,y,z为工件曲面的坐标系;u,v为工件曲面的参数方向。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步说明:
见图3、图4,本发明一种基于曲面截型线分划的最短距离线对获取方法,该方法的具体步骤如下:
步骤一:曲面截型线离散
将工件曲面NU RBS参数化表示为S(u,v),为便于讨论又不失一般性,假设刀具沿曲面v方向进给。因此,可以将工件曲面沿v方向离散成n条截型线,记工件曲面截型线的集合表示为{Mi|i=1,2,…,n}。任意一条截型线可用参数方程表示为
Mi=ρ(ui,v)
其中,ui可表示为第i条截型线的u值,为某一定值。
步骤二:截型线到刀具曲面最短距离点的获取
在刀具加工过程中,任意刀位的切触点可用点Ps(us,vs)表示。在任意刀位处,刀具附近任意截型线i到刀具曲面的最短距离点可按如下流程获取:
1)刀具附近任意截型线可表示为第i条截型线,在第i条截型线上选取初始点Pi0 M(u,v),其中u值与截型线i的u值一致,v值与切触点的v值一致,即u=ui,v=vs。
2)用过点Pi0 M与刀具回转轴的截平面E截刀具与工件曲面,可得截平面内的两条刀具母线,在截平面内连接点Pi0 M与靠近工件曲面的刀具母线圆弧的圆心Oi T,直线Pi0 MOi T与截平面内靠近工件曲面的刀具母线交于点Pi0 T,点Pi0 T即为工件曲面上点Pi0 M到刀具曲面的最短距离点。
4)判断|nPi0 T·τPi0 M|≤ξ是否成立,ξ是一个给定的足够小的控制精度。如果成立,则点Pi0 M即为截型线i到刀具曲面距离最短的点Pi M,最短距离值为
di=|Pi0 MPio T| (7)
由上可知,点Pi0 T与Oi T决定于点Pi0 M,点Pi0 M在截型线i上,截型线上任意点的坐标又取决于参数ui、v。因此点Pi0 T与Oi T均可由参数ui、v表示,截型线上任意点处沿参数v方向的切矢也可由参数ui、v表示。即向量nPi0 T与τPi0 M可表示为参数ui、v的方程:
nPi0 T=G(ui,v)
τPi0 M=H(ui,v)
对于固定的截型线,ui为一定值。因此截型线到刀具曲面最短距离点的获取可转化为方程
A(v)=G(ui,v)·H(ui,v)≈0 (8)
的根的求解。
方程(8)的近似根可按图4中的计算流程求解,求解选取初始值v0=vs,求解范围v∈[vmin,vmax],求解控制精度为ξ,即|A(v)|≤ξ时为所求解。
步骤三:最短距离线对的获取
记工件曲面截型线上到刀具曲面距离最短的点为Pi M,记工件曲面截型线上的最短距离点Pi M对应刀具曲面上的最短距离点为Pi T,其中i表示第i条截型线,i∈[1,n]。将工件曲面上所有的点Pi M顺序连接成线构成工件曲面上的最短距离线,将刀具曲面上的所有点Pi T顺序连接成线构成刀具曲面上的最短距离线,工件曲面上的最短距离线和刀具曲面上的最短距离线便构成最短距离线对。根据最短距离线对间的有向长度便可获得该刀位的误差分布,进而优化刀位。图1为加工过程中最短距离线对示意图,图2为工件曲面截型线到刀具曲面距离最短点的获取示意图。
实施例1:
下面结合附图5所示曲面对曲面内任意一个刀位进行具体数值运算,以鼓形刀侧铣加工为例,鼓形刀刀具母线半径Re=700mm,最大回转半径R=25mm,刀具长度L=50mm,具体步骤如下:
步骤一:曲面截型线离散
将图5所示曲面NURBS参数化表示为S(u,v),曲面参数u、v的取值范围为[0,1]。选择刀具进给方向沿曲面v方向,因此将曲面沿v方向等参数离散成201条截型线,记工件曲面截型线的集合表示为{Mi|i=1,2,…,201}。任意一条截型线可用参数方程表示为
Mi=ρ(ui,v)
其中
步骤二:截型线到刀具曲面最短距离点的获取
任意选择工件曲面切触点为us=0.4,vs=0.6,设刀具在此刀位的侧偏角为0度,刀具切触点在刀具母线的正中间处,因此在该刀位下,刀具投影范围内的曲面截型线可表示为{Mi|i=31,32,…,131}。
首先,对第31条截型线进行计算,选取第31条截型线上的初始点为u31=30/200=0.15,v0=vs=0.6。结合附图4,按照具体实施方式中步骤二的方法可求出第31条截型线上到刀具距离最短的点为u31=0.145,v31=0.614。相同的方法依此可求出第32条到第131条截型线到刀具曲面的最短距离点。
步骤三:最短距离线对的获取
结合图3,按照实施方案中步骤三的方法,依此可求出各曲面截型线上最短距离点对应刀具曲面上的最短距离点及最短距离。将曲面截型线上的最短距离点和刀具曲面上的最短距离点分别依此连成线便得到该刀位下的最短距离线对。以曲面各截型线的u值为横坐标,各截型线到刀具的有向最短距离值为纵坐标,可绘制出该刀位下的最短距离线对的二维分布图。对于该例中最短距离线对的二维分布如图6所示。
Claims (3)
1.一种基于曲面截型线分划的最短距离线对获取方法,其特征在于:该方法的具体步骤如下:
步骤一:曲面截型线离散
将工件曲面沿进给方向离散成n条截型线,记工件曲面截型线的集合表示为{Mi|i=1,2,…,n},任意一条截型线用参数方程表示为
Mi=ρ(ui,v) (1)
其中,ui表示为第i条截型线的u值,为某一定值;
步骤二:截型线到刀具曲面最短距离点的获取
在刀具加工过程中,工件曲面截型线与刀具曲面间的最短距离线段对应着截型线和刀具曲面上各一个相应点,分别称为截型线上的最短距离点和刀具曲面上的最短距离点;刀具曲面上的最短距离点必然在刀具曲面最短距离线上,而截型线上的最短距离点也一定在工件曲面最短距离线上;记刀具曲面最短距离线上的最短距离点的集合为{Pi T|i=1,2,…,n},工件曲面最短距离线上的最短距离点的集合为{Pi M|i=1,2,…,n},点集内的点Pi T与Pi M对应着工件曲面第i条截型线Mi到刀具曲面的最短距离;取nPi T为刀具在点Pi T处的单位法矢,方向指向刀具内部,τPi M为截型线在点Pi M处沿进给方向的单位切线矢量,则有
nPi T·τPi M≈0 (2)
因此,工件曲面截型线到刀具曲面的最短距离点应满足式(2);反之,由式(2)求出工件曲面截型线上到刀具曲面距离最短的点,即获得各截型线到刀具曲面的最短距离点;
步骤三:最短距离线对的获取
根据求出的截型线Mi上到刀具的最短距离点,记为Pi M;过点Pi M与刀具回转轴的截平面截刀具与工件曲面,得截平面内的两条刀具母线,在截平面内连接点Pi M与靠近工件曲面的刀具母线圆弧的圆心Oi T,直线Pi MOi T与截平面内靠近工件曲面的刀具母线交于点Pi T,点Pi T即为工件曲面截型线Mi上的最短距离点Pi M到刀具曲面的最短距离点,|Pi MPi T|的长度值即为截型线Mi与刀具的最短距离值;依此,求出所有截型线到刀具曲面的最短距离值及对应刀具曲面上的最短距离点;分别将所有刀具曲面上的最短距离点和工件曲面截型线上的最短距离点依序连接成线,便获得当前刀位下刀具与工件曲面间的最短距离线对,根据最短距离线对间的有向长度便获得该刀位的误差分布。
2.根据权利要求1所述的一种基于曲面截型线分划的最短距离线对获取方法,其特征在于,步骤二中所述的nPi T·τPi M≈0,是根据求解精度来控制,即给定控制精度为ξ时,当|nPi T·τPi M|≤ξ时获得所求解。
3.根据权利要求1所述的一种基于曲面截型线分划的最短距离线对获取方法,其特征在于,步骤二中所述的截型线到刀具曲面最短距离点的获取,是指当求截型线i上到刀具曲面的最短距离时,截型线i上的初始点根据当前刀位的切触点来定,若当前刀位切触点为Ps(us,vs),则截型线i上选取的初始点为Pi0 M(u,v),其中初始点的u值与截型线i的u值一致,初始点的v值与切触点的v值一致,即u=ui,v=vs。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110144357 CN102222138B (zh) | 2011-05-31 | 2011-05-31 | 一种基于曲面截型线分划的最短距离线对获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110144357 CN102222138B (zh) | 2011-05-31 | 2011-05-31 | 一种基于曲面截型线分划的最短距离线对获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102222138A CN102222138A (zh) | 2011-10-19 |
CN102222138B true CN102222138B (zh) | 2013-10-23 |
Family
ID=44778688
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110144357 Active CN102222138B (zh) | 2011-05-31 | 2011-05-31 | 一种基于曲面截型线分划的最短距离线对获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102222138B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102322407B (zh) * | 2011-08-30 | 2013-06-12 | 重庆大学 | 风力机叶片气动外形协同设计方法 |
CN103034743A (zh) * | 2011-09-30 | 2013-04-10 | 鸿富锦精密工业(深圳)有限公司 | 曲面距离计算***及方法 |
CN103207589B (zh) * | 2013-04-19 | 2015-11-18 | 北京航空航天大学 | 一种环形刀宽行滚切清根加工方法 |
CN103577636B (zh) * | 2013-10-22 | 2016-10-19 | 北京航空航天大学 | 一种基于曲率特征的叶片曲面统一离散方法 |
CN104866655B (zh) * | 2015-05-07 | 2018-04-27 | 北京航空航天大学 | 一种基于包络理论和经线划分的回转型刀具包络特征线计算方法 |
CN108145865B (zh) * | 2017-11-16 | 2019-11-15 | 华中科技大学 | 倒角加工方法及相应装置 |
CN109977999A (zh) * | 2019-02-19 | 2019-07-05 | 山东理工大学 | 一种通过离散工件曲面求解非回转工具有效特征线的方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS57169814A (en) * | 1981-04-10 | 1982-10-19 | Fanuc Ltd | Forming method of curved surface |
CN100582975C (zh) * | 2008-06-19 | 2010-01-20 | 上海交通大学 | 五轴数控加工安全最短刀具长度的规划方法 |
CN101412123B (zh) * | 2008-11-27 | 2011-06-29 | 北京航空航天大学 | 单参数变曲率鼓形刀具及其位姿确定方法和应用 |
-
2011
- 2011-05-31 CN CN 201110144357 patent/CN102222138B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN102222138A (zh) | 2011-10-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102222138B (zh) | 一种基于曲面截型线分划的最短距离线对获取方法 | |
CN102794488B (zh) | 类直纹面整体叶轮曲面的侧铣加工方法 | |
CN101615022B (zh) | 基于曲面重构的三轴数控机床刀具的半径补偿方法 | |
CN102419570B (zh) | 数控机床高速加工的加减速前瞻控制方法 | |
CN102527554B (zh) | 一种自由曲面喷涂机器人的喷枪轨迹规划方法 | |
CN102637216B (zh) | 一种复杂组合曲面的数控侧铣加工刀轨生成方法 | |
CN104615084B (zh) | 加工进给速度优化的刀轨曲线轮廓误差补偿方法 | |
CN105955195B (zh) | 一种基于铣削力预测的曲面加工轨迹生成方法 | |
CN102799145A (zh) | 拐角多周期恒加加速度过渡的s曲线加减速直线插补方法 | |
CN103949705A (zh) | 槽特征腹板摆线螺旋复合铣加工方法 | |
CN102722140A (zh) | 基于s曲线加减速控制的多周期拐角小直线段插补方法 | |
CN102147600A (zh) | 实时生成曲率连续路径的数控插补*** | |
CN109597357B (zh) | 一种面向叶片旋铣工艺的数控编程方法及装置 | |
CN104155916B (zh) | 一种高精度快速空间圆弧插补方法 | |
CN102357666B (zh) | 自由曲面叶片平底刀三坐标端铣加工方法 | |
CN104090528A (zh) | 一种适用于摆线高速铣削的加工路径的方法 | |
CN104007696B (zh) | 飞机结构件型腔复杂转角残留区域刀轨生成方法 | |
CN106647623A (zh) | 一种几何精度及衔接速度最优化的五轴联动平滑插补方法 | |
CN108563848B (zh) | 一种平底螺旋立铣刀的铣削力建模方法 | |
CN108549326B (zh) | 一种类摆线铣削开槽刀轨生成方法 | |
CN103268081A (zh) | 一种用于数控机床的基于精度误差控制的轨迹段转接处理算法 | |
CN110032140B (zh) | 一种五轴加工中球形刀刀轴矢量规划方法 | |
CN103454973B (zh) | 参数化数控立铣刀螺旋槽数控加工编程方法 | |
CN105414616A (zh) | 螺旋铣孔过程中切削力预报和稳定性判别方法 | |
CN102528554A (zh) | 一种五轴加工奇异区域的轨迹优化方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |