CN102873628B - 一种用于数控小工具抛光的螺旋线式加工路径的生成方法 - Google Patents

一种用于数控小工具抛光的螺旋线式加工路径的生成方法 Download PDF

Info

Publication number
CN102873628B
CN102873628B CN201210363814.2A CN201210363814A CN102873628B CN 102873628 B CN102873628 B CN 102873628B CN 201210363814 A CN201210363814 A CN 201210363814A CN 102873628 B CN102873628 B CN 102873628B
Authority
CN
China
Prior art keywords
machining path
helix
polar
polishing
represent
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.)
Expired - Fee Related
Application number
CN201210363814.2A
Other languages
English (en)
Other versions
CN102873628A (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.)
Tsinghua University
Original Assignee
Tsinghua 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 Tsinghua University filed Critical Tsinghua University
Priority to CN201210363814.2A priority Critical patent/CN102873628B/zh
Publication of CN102873628A publication Critical patent/CN102873628A/zh
Application granted granted Critical
Publication of CN102873628B publication Critical patent/CN102873628B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

一种用于数控小工具抛光的螺旋线式加工路径的生成方法,属于光学加工技术领域。为了解决现有螺旋线抛光路径容易导致中频误差较大的问题,通过引入权值因子构造了一种生成过程中不再保持恒定的角度增量和半径增量的新的离散螺旋线加工路径,该螺旋线加工路径在以θ为极角,r为半径坐标值的极坐标系θ-r下,由n个离散点构成,并且由中心向外依次生成,到中心的距离也依次变大,有利于减小工件面形的同心圆状的波纹误差和辐射状波纹误差这两种中频误差,进一步提高光学工件精度。

Description

一种用于数控小工具抛光的螺旋线式加工路径的生成方法
技术领域
本发明属于光学加工技术领域,特别涉及一种用于数控小工具抛光的螺旋线式加工路径的生成方法。 
背景技术
随着人们对光学系统精度的要求越来越高,光学工件明显呈现出向大尺度、高精度和非球面方向发展的态势。抛光作为光学工件成形的最后工序,对现代光学工件的制造至关重要。 
然而传统光学工件的抛光技术使用比工件更大或近似同等大小的抛光工具,在对大尺度工件进行抛光时尤为不便;传统抛光方法加工过程中,工件表面形状(面形)误差收敛较慢,甚至很难收敛至高精度光学系统的允许范围;使用传统抛光方法对非球面加工过程中,存在加工精度低、产品质量不稳定、对工人技术要求高等问题。 
近年来出现的数控小研抛盘抛光、气囊抛光、磁流变抛光、射流抛光和离子束抛光等技术都属于数控小工具抛光技术,解决了传统光学抛光方法所遇到的问题,大大提高了工件成形精度。数控小工具抛光技术使用较小的抛光工具对工件表面材料进行局部去除,然后在数控技术的辅助下按照一定的运动路径和变化的速率遍历整个表面,快速实现对误差的针对性修正,从而达到提高工件面形精度的目的。 
光学系统中大量使用圆形光学工件,目前在数控小工具抛光技术中主要对这类工件使用阿基米德螺旋线进行路径规划。如图1所示。 
阿基米德螺旋线通常在极坐标(θ-r)下表示为: 
r=c·θ        (1) 
其中,c为常数。阿基米德螺旋线的相邻两圈的间隔Δr为定值,即Δr=2πc。 
为了计算小工具在工件表面各点的加工时间(驻留时间)分布和数控编程,必须将连续的阿基米德螺旋线离散化。目前通常大多采用角度增量Δθ恒定的方式进行离散。 
ri=c·Δθ·i        (2) 
其中,i为整数。离散后的阿基米德螺旋线如图2所示。 
在路径规划时采用此种螺旋线容易导致加工后的工件面形出现两种形状的中频误差:同心圆状的波纹误差和辐射状波纹误差。由于螺旋线的相邻两圈的间隔Δr为定值,也就是说小工具在工件表面运动时与相邻的内圈或外圈的距离总是完全相等,所以容易导致同心圆状的波纹误差的出现;按照角度增量Δθ恒定方式对螺旋线进行离散,势必导致离散后的点阵列呈现放射状,如图2所示,明显可以看出,最外圈上相邻点之间的间距比最内圈上的相邻点之间的间距大得多,过大的间距意味着对外圈驻留时间的计算精度和控制精度不足,所以容易导致辐射状波纹误差的出现。 
为了解决对中心区域加工时工件转动速度过快的问题,中国发明专利文献公开了一种螺 旋线抛光路径(申请号200810030958.X),如图3所示。该螺旋线在极坐标(θ-r)下的方程为r=bθ1/2,其中b为常数。该螺旋线的特点在于其所包围面积的增大速率保持恒定,即dA/dθ=πb2。 
为了满足面积的增大速率保持恒定,相邻两圈的间隔Δr不再恒定,而是由内向外逐渐变小,有利于避免同心圆状的波纹误差;但该螺旋线依然是在角度增量Δθ恒定的条件下进行离散的(呈现放射状),依然容易导致辐射状波纹误差的出现。 
发明内容
本发明的目的是提供一种用于数控小工具抛光的螺旋线式加工路径的生成方法,以解决现有路径容易导致同心圆状的波纹误差和辐射状波纹误差两种中频误差出现的问题,为减小小工具抛光技术的中频误差提供一种更优的加工路径。 
本发明的发明内容如下:一种用于数控小工具抛光的螺旋线式加工路径的生成方法,该螺旋线加工路径在以θ为极角,r为半径坐标值的极坐标系θ-r下,由n个离散点构成,即{P1,P2,P3…Pk,Pk+1…Pn},并且由中心向外依次生成,到中心的距离也依次变大;螺旋线加工路径的起点P1位于极坐标系θ-r的极点,点P1极坐标为(θ1,r1)=(0,0);其特征在于,螺旋线加工路径遵循如下几何关系: 
k=1,2,…n-1 n,k均为正整数, 
θk+1=θk+Δθk k=1,2,…n-1 n,k均为正整数, 
Δθ k = arccos ( ( r k ) 2 + ( r k + 1 ) 2 - Δl 2 2 r k · r k + 1 ) k ≥ 2 Δθ k = Δθ 1 k = 1 k=1,2,…n-1 n,k均为正整数, 
其中:k,rk)和(θk+1,rk+1)分别为螺旋线加工路径上的第k个点Pk以及与Pk相邻的第k+1个点Pk+1的极坐标;rk表示点Pk与极点的距离值,θk表示点Pk的极角,即极轴按逆时针方向转动到射线P1Pk所在位置经过的角度值;rk+1表示点Pk+1与极点的距离值,θk+1表示点Pk+1的极角,即极轴按逆时针方向转动到射线P1Pk+1所在位置经过的角度值;rk<rk+1;Δl表示点Pk与点Pk+1之间的距离,当k≥2时,0<Δl<(rk+rk+1);Δθk表示点Pk+1与点Pk极角的差值;由于k=1时,点P2的极角θ2=θ1+Δθ1=Δθ1,故Δθ1与点P2的极角相等;Δθ1、Δl为给定的已知数,ak、bk分别为给定的已知序列a、b的第k个值;ak>0,bk>0,ak+bk=1。 
本发明具有以下优点及突出性效果:本发明所述的一种用于数控小工具抛光的螺旋线式加工路径的生成方法,为了解决现有螺旋线抛光路径容易导致中频误差较大的问题,引入了权值因子,具有可变性强,适应性广的优点。该螺旋线加工路径在以θ为极角,r为半径坐标值的极坐标系θ-r下,由n个离散点构成,并且由中心向外依次生成,到中心的距离也依次变大。通过给定不同的已知参数,可实现螺旋线在生成过程中不再保持恒定的角度增量和半 径增量,从而有利于减小工件面形的同心圆状的波纹误差和辐射状波纹误差这两种中频误差,进一步提高光学工件精度。 
附图说明
图1为阿基米德螺旋线。 
图2为离散化的阿基米德螺旋线。 
图3为包围面积的增大速率保持恒定的螺旋线。 
图4为本发明所述螺旋线加工路径示意图。 
图5为本发明所述螺旋线加工路径相邻两离散点的几何关系。 
图6为本发明所述螺旋线加工路径离散点分布(局部)。 
具体实施方式
图4为本发明所述螺旋线加工路径,该螺旋线加工路径在以θ为极角,r为半径坐标值的极坐标系θ-r下,由n个离散点构成,即{P1,P2,P3…Pk,Pk+1…Pn},并且由中心向外依次生成,到中心的距离也依次变大;螺旋线加工路径的起点P1位于极坐标系θ-r的极点,点P1极坐标为(θ1,r1)=(0,0);其特征在于,螺旋线加工路径遵循如下几何关系: 
k=1,2,…n-1 n,k均为正整数, 
θk+1=θk+Δθk k=1,2,…n-1 n,k均为正整数, 
&Delta;&theta; k = arccos ( ( r k ) 2 + ( r k + 1 ) 2 - &Delta;l 2 2 r k &CenterDot; r k + 1 ) k &GreaterEqual; 2 &Delta;&theta; k = &Delta;&theta; 1 k = 1 k=1,2,…n-1 n,k均为正整数, 
其中:k,rk)和(θk+1,rk+1)分别为螺旋线加工路径上的第k个点Pk以及与Pk相邻的第k+1个点Pk+1的极坐标;rk表示点Pk与极点的距离值,θk表示点Pk的极角,即极轴按逆时针方向转动到射线P1Pk所在位置经过的角度值;rk+1表示点Pk+1与极点的距离值,θk+1表示点Pk+1的极角,即极轴按逆时针方向转动到射线P1Pk+1所在位置经过的角度值;rk<rk+1;Δl表示点Pk与点Pk+1之间的距离,当k≥2时,0<Δl<(rk+rk+1);Δθk表示点Pk+1与点Pk极角的差值;由于k=1时,点P2的极角θ2=θ1+Δθ1=Δθ1,故Δθ1与点P2的极角相等;Δθ1、Δl为给定的已知数,ak、bk分别为给定的已知序列a、b的第k个值;ak>0,bk>0,ak+bk=1。 
图5为本发明所述螺旋线加工路径相邻两离散点的几何关系,线段P1Pk、P1Pk'和P1Pk+1的长度分别为点Pk、Pk'和Pk+1在极坐标系中的半径值:rk、rk'和rk+1。 
以工件表面80mm×80mm范围为研究对象,令Δθ1=0.05,Δl=1mm,n=10204,ak=0.6+(k·(0.8-0.6)/(n-1)),bk=0.4+(k·(0.2-0.4)/(n-1))。 
使用计算机依次生成(θ1,r1),(θ2,r2)…(θk,rk),(θk+1,rk+1)…(θn,rn),形成离散抛光路径,由于点数过多过密,在图6中仅显示了一部分。以磁流变抛光为例,利用测量仪器测得的单 位去除函数和面形初始误差,求解抛光路径各点所对应的驻留时间。为了排除实验中不确定性因素的影响,采用仿真算法求解按照本发明所述螺旋线抛光路径和所得驻留时间进行加工后的面形误差,没有发现明显的同心圆状的波纹误差和辐射状波纹误差。 
令本发明所述的螺旋线方程中的ak≡0,bk≡1,其他参数不变,则本发明所述的螺旋线退化为以勾股定理逐点生成的螺旋线抛光路径。按照上述方法计算加工后的面形误差,从误差分布图中可以发现同心圆状的波纹误差。 
使用角度增量Δθ恒定的方式对阿基米德螺旋线进行离散,形成抛光路径。按照上述方法计算加工后的面形误差,从误差分布图中可以发现辐射状波纹误差。 
可见,按照本发明所述方法生成的螺旋线加工路径,对同心圆状的波纹误差和辐射状波纹误差这两种中频误差由一定的抑制作用。由本发明所述方法引入了由ak,bk构成的权值因子,使螺旋线在生成过程中具有角度增量变化和半径增量变化的特征,从而使中频误差得以减弱。 

Claims (1)

1.一种用于数控小工具抛光的螺旋线式加工路径的生成方法,该螺旋线加工路径在以θ为极角,r为半径坐标值的极坐标系θ-r下,由n个离散点构成,即{P1,P2,P3…Pk,Pk+1…Pn},并且由中心向外依次生成,到中心的距离也依次变大;螺旋线加工路径的起点P1位于极坐标系θ-r的极点,点P1极坐标为(θ1,r1)=(0,0);其特征在于,螺旋线加工路径遵循如下几何关系:
k=1,2,…n-1 n,k均为正整数,
θk+1=θk+Δθk k=1,2,…n-1 n,k均为正整数,
&Delta;&theta; k = arccos ( ( r k ) 2 + ( r k + 1 ) 2 - &Delta;l 2 2 r k &CenterDot; r k + 1 ) k &GreaterEqual; 2 &Delta;&theta; k = &Delta;&theta; 1 k = 1 k=1,2,…n-1 n,k均为正整数,
其中:k,rk)和(θk+1,rk+1)分别为螺旋线加工路径上的第k个点Pk以及与Pk相邻的第k+1个点Pk+1的极坐标;rk表示点Pk与极点的距离值,θk表示点Pk的极角,即极轴按逆时针方向转动到射线P1Pk所在位置经过的角度值;rk+1表示点Pk+1与极点的距离值,θk+1表示点Pk+1的极角,即极轴按逆时针方向转动到射线P1Pk+1所在位置经过的角度值;rk<rk+1;Δl表示点Pk与点Pk+1之间的距离,当k≥2时,0<Δl<(rk+rk+1);Δθk表示点Pk+1与点Pk极角的差值;由于k=1时,点P2的极角θ2=θ1+Δθ1=Δθ1,故Δθ1与点P2的极角相等;Δθ1、Δl为给定的已知数,ak、bk分别为给定的已知序列a、b的第k个值;ak>0,bk>0,ak+bk=1。
CN201210363814.2A 2012-09-26 2012-09-26 一种用于数控小工具抛光的螺旋线式加工路径的生成方法 Expired - Fee Related CN102873628B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210363814.2A CN102873628B (zh) 2012-09-26 2012-09-26 一种用于数控小工具抛光的螺旋线式加工路径的生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210363814.2A CN102873628B (zh) 2012-09-26 2012-09-26 一种用于数控小工具抛光的螺旋线式加工路径的生成方法

Publications (2)

Publication Number Publication Date
CN102873628A CN102873628A (zh) 2013-01-16
CN102873628B true CN102873628B (zh) 2015-02-18

Family

ID=47475260

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210363814.2A Expired - Fee Related CN102873628B (zh) 2012-09-26 2012-09-26 一种用于数控小工具抛光的螺旋线式加工路径的生成方法

Country Status (1)

Country Link
CN (1) CN102873628B (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106625032A (zh) * 2016-11-03 2017-05-10 天津津航技术物理研究所 一种螺旋正弦式小工具抛光去除金刚石刀痕的方法
CN109922923B (zh) * 2016-11-08 2020-10-23 Abb瑞士股份有限公司 抛光工件的方法以及使用该方法的系统
CN107520683B (zh) * 2017-08-31 2019-03-19 中国工程物理研究院激光聚变研究中心 子口径抛光的等高线路径规划方法
DE102018101293A1 (de) * 2018-01-22 2019-07-25 Rud. Starcke Gmbh & Co. Kg Verfahren zum Schleifen und/oder Polieren einer Fehlstelle sowie Vorrichtung zur Durchführung des Verfahrens
CN110456734A (zh) * 2019-08-24 2019-11-15 南通苏尼康自动化技术有限公司 一种电熔管件加工设备用控制器及其控制方法
CN111070028B (zh) * 2019-11-26 2022-06-14 天津津航技术物理研究所 一种非旋转对称表面光学加工轨迹设计方法
CN111079889B (zh) * 2019-12-13 2022-07-05 吉林大学 改进的基于分解的多目标粒子群规划螺旋线抛光轨迹方法
CN111015671B (zh) * 2019-12-30 2022-03-08 南京埃斯顿机器人工程有限公司 一种机器人平面螺旋线打磨轨迹的规划方法
CN113370071B (zh) * 2021-06-18 2022-08-09 深圳大学 一种轴对称元器件的研磨方法及装置
CN114474075B (zh) * 2022-03-28 2023-03-28 法奥意威(苏州)机器人系统有限公司 机器人螺旋轨迹控制方法、装置、存储介质及电子设备
CN115302357B (zh) * 2022-08-05 2023-05-16 中国人民解放军空军工程大学航空机务士官学校 一种基于评价函数的螺旋抛光路径规划方法

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4128968A (en) * 1976-09-22 1978-12-12 The Perkin-Elmer Corporation Optical surface polisher
US5067085A (en) * 1989-05-15 1991-11-19 Southwest Research Institute Optical robotic canopy polishing system
JP2002535151A (ja) * 1998-12-01 2002-10-22 ユニヴァーシティ カレッジ ロンドン 研磨装置および方法
GB2452091B (en) * 2007-08-24 2013-01-02 Zeeko Ltd Computer controlled work tool apparatus and method
CN100589054C (zh) * 2008-03-31 2010-02-10 中国人民解放军国防科学技术大学 离子束抛光工艺中面形收敛精度的控制方法
CN101274822B (zh) * 2008-03-31 2010-08-11 中国人民解放军国防科学技术大学 一种离子束抛光路径的规划方法
CN101456681B (zh) * 2009-01-06 2011-05-18 中国人民解放军国防科学技术大学 能修除光学元件局部误差的离子束极轴加工方法
CN102059638B (zh) * 2010-11-19 2012-09-05 苏州大学 一种计算机控制小工具抛光方法

Also Published As

Publication number Publication date
CN102873628A (zh) 2013-01-16

Similar Documents

Publication Publication Date Title
CN102873628B (zh) 一种用于数控小工具抛光的螺旋线式加工路径的生成方法
CN111061213B (zh) 一种基于Bezier曲线转角平滑过渡算法的加工方法
CN108062073B (zh) 一种用于高质量加工的圆弧平滑压缩插补方法
CN107562013B (zh) 等残余高度的曲面分区变刀具加工轨迹规划方法
CN110426992B (zh) 一种曲面子区域环切加工刀具轨迹规划方法
CN101539769A (zh) 基于二次b样条曲线对g01代码的拟合及插补方法
CN104238456B (zh) 一种非球头刀铣削加工自由曲面的方法
CN106875439B (zh) 基于三维点云模型的单晶硅棒外形尺寸测量方法
CN102019572A (zh) 采用组合螺旋式抛光路径的抛光工艺
CN110076680B (zh) 一种近轴端远轴端等厚离轴非球面加工方法
CN105069271A (zh) 一种锥齿轮齿面加工误差修正方法
CN102222138B (zh) 一种基于曲面截型线分划的最短距离线对获取方法
CN109270892B (zh) 一种数控加工中非圆曲线的最小二乘螺旋线逼近算法
CN106774153A (zh) 一种s形球头铣刀的轮廓度补偿方法
CN109597354A (zh) 一种三角网格模型的多约束数控加工轨迹生成方法
CN108000731A (zh) 一种圆盘锯切割石材的圆弧和椭圆轮廓的方法
CN103809510B (zh) 一种面向高精加工的自由曲面往复式刀具轨迹规划方法
CN102778859B (zh) 基于双螺旋空间填充曲线的数控加工刀具轨迹生成方法
CN107103136B (zh) 一种基于距离场的轨道车辆轮轨接触点计算方法
CN105867298B (zh) 螺旋曲面展成加工中确定刀具廓形的点矢量二次包络法
CN101113573B (zh) 基于几何关系渐进迭代的三角股钢丝绳的配丝方法
CN102902849A (zh) 一种公路缓和曲线的简化设计方法
CN105785919A (zh) 点云五轴无全局干涉刀轴矢量可行域计算方法
CN110497261A (zh) 一种基于像素法获取整体立铣刀容削槽端面截形的方法
CN114036668A (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
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20150218

Termination date: 20200926