首页> 中国专利> 数控磨床凸轮磨削的速度优化控制方法

数控磨床凸轮磨削的速度优化控制方法

摘要

本发明涉及一种数控领域的凸轮磨削方法,具体涉及一种数控磨床凸轮磨削的速度优化控制方法。目的是利用本发明速度优化算法对凸轮进行数控磨削,实现了对凸轮升程曲线中斜率较大处和甚小处速度的明显优化,最终实现在保证凸轮磨削的精度等工艺指标前提下,提高凸轮磨削的效率。本发明利用反转法建立磨削过程运动学模型;在恒角速磨削的基础上利用经典公式对砂轮进给速度、加速度和凸轮旋转角速度、角加速度进行优化;利用优化后的值自动生成磨削G代码;利用G代码进行实际磨削得出凸轮轮廓误差对原升程进行补偿。该发明解决了常用的恒角速度磨削产生的非线性弹性变形和恒线速度磨削的过切或切削不充分现象,明显提高了凸轮磨削精度和加工效率。

著录项

  • 公开/公告号CN103802023A

    专利类型发明专利

  • 公开/公告日2014-05-21

    原文格式PDF

  • 申请/专利权人 吉林大学;

    申请/专利号CN201410081836.9

  • 发明设计人 王静;隋振;田彦涛;王勋龙;李琳;

    申请日2014-03-06

  • 分类号B24B51/00(20060101);

  • 代理机构22201 长春吉大专利代理有限责任公司;

  • 代理人朱世林;王寿珍

  • 地址 130012 吉林省长春市前进大街2699号

  • 入库时间 2024-02-19 23:06:30

法律信息

  • 法律状态公告日

    法律状态信息

    法律状态

  • 2016-04-27

    授权

    授权

  • 2014-06-25

    实质审查的生效 IPC(主分类):B24B51/00 申请日:20140306

    实质审查的生效

  • 2014-05-21

    公开

    公开

说明书

技术领域

本发明涉及一种数控领域的凸轮磨削方法,具体涉及一种数控磨床凸轮磨削的速度优化 控制方法。

背景技术

凸轮磨削在数控领域的发展日新月异,高精度、高速磨削依然是大家关注和追求的目标。 凸轮的磨削过程是靠旋转轴(C轴)和砂轮进给轴(X轴)的配合实现的,在两轴位置关系 一定的前提下,两轴的速度关系是随着凸轮转角的不同而时刻变化的。不同的速度关系,在 两轴上引起的动态误差是不同的,进而直接影响了凸轮的磨削精度。在整个加工系统中,砂 轮架质量最大,如果凸轮转速过快,就会导致磨头跟踪不到位而产生过切或切削不充分等现 象。尤其是在升程曲线中斜率较大的情况下,若速度和加速度依然大,就会造成两个轴的滞 后“量”不同,导致两轴不能同步,最后影响轮廓精度。若速度过慢,又会降低整个磨削过 程的效率,影响企业的产量。

对于此问题,国内的不少学者也做了一定的工作,并取得了相应的成果。他们对凸轮旋 转角速度的变化规律进行优化,使两轴的速度及加速度很好配合,以两轴的动态跟踪误差最 小为原则,对提高磨削精度的作用十分明显,这在人工速度优化时已得到验证。目前,国内 基于该理论的速度优化途径通常有两种:一是通过砂轮架允许的最大速度和加速度来反向调 节凸轮旋转的速度和加速度,王洪等人根据砂轮架X轴的运动速度和加速度来预测C轴的 旋转速度(中国专利:201010278922.0、“一种凸轮轴数控磨削加工方法”),该方法易于实 现程序自动化,但计算相对比较复杂;另一种是正向调节凸轮旋转轴的速度,计算砂轮架的 速度和加速度,并限制最大速度和加速度在允许范围内,上海大学的何永义等采用三角函数 逼近的方法来解决头架旋转加速度局部过大的问题(中国专利:201110058262.X)、湖南大 学的曹德芳和邓朝辉在文献“凸轮轴磨削加工速度优化调节与自动数控编程研究”中利用三 次样条函数拟合插值的方法来进行速度优化,实现了比较好的磨削效果,但这些优化算法都 比较依赖数控磨床的机械结构,算法的适应性不够强。与国内的数控磨削相比,国外的数控 磨削技术相对来说比较成熟,即使在同样的数控磨床上进行磨削加工,不论在磨削效率还是 磨削精度方面都远远超越了国内,但为了商机和盈利起见,国外公开的优化算法甚少。

本发明针对前述问题,提出一种新的凸轮数控磨削速度优化控制方法,即先计算出砂轮 架速度和加速度后,对砂轮架速度和加速度进行计算优化,然后反向调整凸轮转速。进一步 提出了基于恒角速度磨削,利用经典公式(本发明采用指数函数)优化砂轮速度和加速度, 进而调节凸轮旋转速度,以实现在凸轮升程曲线斜率较大处速度减缓、斜率较小处速度提高 的目的。

技术内容

本发明针对现有技术存在的上述不足,提供一种数控磨床凸轮磨削的速度优化控制方 法。其发明目的在于利用发明的速度优化算法对凸轮进行数控磨削,实现对凸轮升程曲线中 斜率较大处和甚小处速度的明显优化,最终实现在保证凸轮磨削的效率的前提下,提高凸轮 磨削精度。

为了达到这个目的,本发明采用反转法来描述两个轴的运动学关系,即假设凸轮并不转 动,而砂轮以数值相同而方向相反的角速度方向绕凸轮反转。在此基础上,本发明提出的速 度优化的磨削方法利用计算机的数控程序来控制整个磨削过程,具体步骤描述如下:

步骤一,向计算机输入凸轮升程表中的基本数据,包括:基圆半径、滚子半径、砂轮半 径、从动件的旋转角度以及对应的升程、初始的恒角速度;

步骤二,建立凸轮磨削过程运动学模型,计算砂轮架的进给量和凸轮的旋转角度;

步骤三,为了避免直接计算出来的凸轮旋转角度与砂轮进给量疏密化影响的严重性,对 其值进行三次样条插值计算,计算出砂轮圆心轨迹曲线上各点的值;

步骤四,在恒角速度基础上,基于经典公式设计速度调节系数,对砂轮进给速度、加速 度和凸轮旋转轴的速度和加速度进行优化;

步骤五,对步骤四中计算的修正后砂轮架进给速度、加速度和修正后的凸轮转速、加速 度进行三次样条曲线拟合,得到光滑曲线,由编程软件自动生成数控加工子程序G代码;

步骤六,利用以上得到的G代码对实际凸轮片进行实际磨削,并利用凸轮检测仪检测 加工的凸轮,计算轮廓误差,并对升程进行补偿,最后再由计算机自动得出数控加工G代 码,接着对凸轮进行磨削,直至磨削结果获得满意。

步骤二所述的磨削过程运动学模型建立步骤如下:

步骤一,按几何关系得出以下坐标表达式:

x=x1-r0·Vxy=y1-r0·Vy---(3)

式中,r为凸轮基圆半径,r0为滚轮半径,为凸轮的升程值,为从动件的旋转角度, (x1,y1)为滚轮圆心的轨迹坐标,即理论轮廓,(Vx,Vy)为滚轮圆心的单位法向量,(x,y)为实 际凸轮轮廓的表达式;

步骤二,利用反转法思想,假定凸轮并不转动,而砂轮一方面以大小相同但方向相反的 角速度绕凸轮反转,同时砂轮又以原有运动规律相对机架往复运动,此时砂轮中心点的运动 轨迹与凸轮轮廓轨迹相比,二者是间距为R的等距曲线,且二者的方向向量相同,这样砂轮 中心点的轨迹坐标可由以下公式获得:

x2=x+R·Vxy2=y+R·Vy---(4)

式中,R为砂轮半径,(x2,y2)为砂轮中心点的轨迹坐标;

步骤三,计算出砂轮进给量和凸轮旋转角度:

X=x22+y22---(5)

θ=arctany2x2---(6)

式中,X为砂轮进给量,θ为凸轮的旋转角度。

步骤四所述的基于经典公式设计速度调节系数,对砂轮进给速度、加速度和凸轮旋转轴 的速度和加速度进行优化的算法包括以下步骤:

步骤一,设定恒角速度磨削,其中n为常数;

步骤二,计算砂轮架进给速度和加速度:

vx=dXdt=dX·dt---(7)

ax=dvxdt=dvx·dt---(8)

步骤三,砂轮架相对于凸轮转角的速度:

vjx=dX---(9)

步骤四,给出整体速度调节系数:

m=p|vjx|×k---(10)

其中,p为减小速度变化率及减缓升程波动的速度调节系数(p>1),k为全局速度调节系 数;

步骤五,调节后的凸轮转速和砂轮架进给速度分别由以下式子获得:

wmc=wc/m                                    (11)

vmx=vx/m                                    (2)

步骤六,调整结束后,再次计算优化后的砂轮架以及凸轮旋转轴的加速度:

amx=dvmx·n---(13)

amc=dwmc·n---(14).

步骤六所述的对升程进行补偿是利用以下公式:

式中,α为补偿系数,取0.5—1.0,Δx为计算出来的轮廓误差。

步骤四中的所述基于经典公式设计速度调节系数,调节步骤有如下:

步骤一,给赋初值p=1(即为恒角速度磨削控制),在升程曲线斜率较大的情况下,会 导致速度和加速度很大,这样就容易引起砂轮进给轴的跟踪误差增大,使X-C两轴联动后的 凸轮轮廓误差变大,严重影响精度,这时,我们就需要调节系数p,使其增大,此时m值增 大,从而来降低速度和加速度,减小跟踪误差,从而提高凸轮轮廓精度;

步骤二,当p值固定,若发现整体的凸轮转速过快时,为了防止磨头跟踪不到位而产生 过切或切削不充分等现象,我们就需要调节k,使其减小;若整体速度过慢,为了提高效率, 便增大k值。

本发明的方法与现有技术相比,有以下几点优势:

1)利用反转法建立的凸轮磨削的运动学模型与目前通用的模型相比,易理解、计算方 便且简单;

2)利用指数函数所构造的速度调节系数,其中的调整参数物理意义很明确,能够有效 的对砂轮架及凸轮的速度、加速度进行优化;

3)本发明的方法不依赖于数控磨床的机械结构,适应性较强。针对不同的数控磨床, 进行适当的参数调整,即可取得很好的磨削效果;

4)本发明的速度优化方法既避免了非线性弹性变形的出现,又解决了升程曲线斜率较 大处速度过大造成的过切或切削不充分的问题,同时也可实现在升程曲斜率较小处速度提高 的目的。

总而言之,本发明在保证了加工效率的情况下,有效提高了磨削精度,降低了振纹现象 的出现。

附图说明

本发明将通过示例,参考下述附图以更进一步的阐述:

图1为速度优化后磨削的凸轮轮廓图。

图2为方案实施的流程图。

图3为凸轮磨削运动学模型示意图。

图4为砂轮架进给位移曲线。

图5为速度优化之前的凸轮旋转角度和砂轮进给速度曲线图。

图6为速度优化之后的凸轮旋转角度和砂轮进给速度曲线图。

图7为磨削后检测后得到的误差图。

图8为磨削后检测得到的相邻误差图。

具体实施方式

以下进一步说明本发明的具体内容及其实施方式:

本发明提出的速度优化的磨削方法利用计算机的数控程序来控制整个磨削过程,具体步 骤描述如下:

1)向计算机输入凸轮升程表中的基本数据,包括基圆半径、滚子半径、砂轮半径、从 动件的旋转角度以及对应的升程、初始的恒角速度;

2)建立凸轮磨削运动学模型,利用反转法及几何关系,推导并计算砂轮架的进给量和 凸轮的旋转角度:

x=x1-r0·Vxy=y1-r0·Vy---(3)

x2=x+R·Vxy2=y+R·Vy---(4)

X=x22+y22---(5)

θ=arctany2x2---(6)

式中,r为凸轮基圆半径,r0为滚轮半径,R为砂轮半径,为凸轮的升程值,为从动 件的旋转角度,(x1,y1)为滚轮圆心的轨迹坐标,即理论轮廓,(Vx,Vy)为滚轮圆心的单位法向 量,(x,y)为实际凸轮轮廓的表达式,(x2,y2)为砂轮中心点的轨迹坐标,X为砂轮进给量, θ为凸轮的旋转角度。

3)为了避免直接计算出来的旋转角度与砂轮进给量疏密化影响的严重性,对以上公式 (5)、(6)计算出来的C轴旋转角度及对应的X轴进给量进行三次样条插值计算,计算出 砂轮圆心轨迹曲线上各点的值。

4)在假定恒角速度情况下来计算速度、加速度,并利用经典公式设计速度调节系数, 修正速度和加速度。设定

wc=dt=·dt---(7)

vx=dXdt=dX·dt---(8)

ac=dwcdt=dwc·dt---(9)

ax=dvxdt=dvx·dt---(10)

vjx=dX---(11)

m=p|vjx|×k---(12)

wmc=wc/m                                  (13)

vmx=vx/m                                  (14)

amx=dvmx·n---(15)

amc=dwmc·n---(16)

式中,n为初始的恒角速度,wc为凸轮的转速,vx为砂轮架的进给速度,ac为凸轮的旋转 加速度,ax为砂轮架的进给加速度,vjx为砂轮架相对于凸轮转角的速度,vmx为修正后的砂 轮架进给速度,wmc为修正后的凸轮转速,amx为修正后的砂轮架进给加速度,amc为修正后 的凸轮旋转加速度,k为全局速度调节系数,p为减小变化率及减缓升程波动的速度调节系 数,m为整体的速度调节系数。

5)对上述步骤中计算的修正后砂轮架进给速度、加速度和修正后的凸轮转速、加速度 进行三次样条曲线拟合,得到光滑曲线,由编程软件自动生成数控加工子程序G代码。

6)利用以上得到的G代码对实际凸轮片进行实际磨削,并利用凸轮检测仪检测加工的 凸轮,计算轮廓误差,并对升程进行补偿,最后再由计算机自动得出数控加工G代码,接 着对凸轮进行磨削,直至磨削结果获得满意即可。

本发明一种数控磨床凸轮磨削的速度优化控制方法,以一个优选实施例并结合附图说明 及具体操作步骤如下:

利用计算机的数控程序来控制整个磨削过程,主要是控制作为X轴的砂轮架横向进给 量和作为C轴的头架的旋转角度来实现高速、高精度磨削。通常凸轮机构的从动件有三种不 同的形式:刀尖挺杆、滚子挺杆和平面挺杆。相对来说,滚子挺杆比较常见,且建立模型及 利用反转法计算较其他两种挺杆来说稍稍复杂,故本实施例选择滚子挺杆。在选择凸轮片方 面,本实施例选择了一个表面有凹凸、且不对称的凸轮,相对来说磨削加工比较复杂、困难, 图1为速度优化后的磨削轮廓图,其中外层和内层的曲线分别为国家标准的误差范围带 (±0.02mm)。方案实施的流程图如图2。

1)建立凸轮磨削运动学模型如图3,图3中标号1为挺杆滚子,2为砂轮,3为凸轮的 实际轮廓,4为理论轮廓,o为凸轮中心,o1为滚轮圆心,o2为砂轮中心,r为凸轮基圆 半径,r0为滚轮半径,R为砂轮半径,为从动件的旋转角度,θ为凸轮的旋转角度,P点 为凸轮轮廓与砂轮相切的切点。

2)根据某型号数控磨床提供的上述凸轮的升程表数据,向计算机输入初始参数:基圆 半径r=19(mm)、滚子半径r0=11(mm)、砂轮半径R=290(mm)、初始的角速度(100r/min) 及给定凸轮旋转的每一度所对应的升程(单位为:mm),总计360个数据(见表1),得出 砂轮架进给位移曲线(见图4)。

3)利用反转法及图3中的几何关系,推导并计算砂轮架的进给量和凸轮的旋转角度:

X=x22+y22---(1)

θ=arctany2x2---(2)

其中,(x2,y2)为砂轮中心点的轨迹坐标,X为砂轮进给量,θ为凸轮的旋转角度。

4)对以上公式(1)、(2)计算出来的C轴旋转角度及对应的X轴进给量进行三次样条 插值计算,计算出砂轮圆心轨迹曲线上各点的值。

5)在假定恒角速度情况下来计算速度、加速度,并利用经典公式设计速度调节系数, 修正速度和加速度。设定

vjx=dX---(3)

m=p|vjx|×k---(4)

wmc=wc/m                                       (5)

vmx=vx/m                                       (6)

amx=dvmx·n---(7)

amc=dwmc·n---(8)

式中,n为初始的恒角速度,vjx为砂轮架相对于凸轮转角的速度,vmx为修正后的砂轮架进 给速度,wmc为修正后的凸轮转速,amx为修正后的砂轮架进给加速度,amc为修正后的凸轮 旋转加速度,k为全局速度调节系数,p为减小变化率及减缓升程波动的速度调节系数,m 为整体的速度调节系数。

6)由给定的数据(见表1)及得出的砂轮架进给位移曲线(见图4),根据公式(3), 可计算出砂轮架相对于凸轮转角的速度。

7)根据|vjx|的值,来调节m值。

调节中,首先给赋初值p=1(即为恒角速度磨削控制),在升程突变的情况下即|vjx|较 大,我们就需要调节系数p,使其增大,此时m值增大,从而来降低速度和加速度,减小跟 踪误差,从而提高凸轮轮廓精度。根据给定的凸轮升程数据,最终调节p=10。

当p值固定,若发现整体的凸轮转速过快时,为了防止磨头跟踪不到位而产生过切或切 削不充分等现象,我们就需要调节k,使其减小;若整体速度过慢,为了提高效率,便增大 k值。针对本例,选择k=1.5。

8)根据公式(5)、(6)、(7)、(8),计算的修正后砂轮架进给速度vmx、加速度 amx和修正后的凸轮转速wmc、加速度amc进行三次样条曲线拟合,通过matlab软件对砂轮架 的速度和凸轮旋转角速度进行仿真,得到速度优化之前和之后的轮进给速度和凸轮旋转速度 对比曲线(见图5和图6)。

9)利用速度优化算法对速度进行优化,由编程软件自动生成数控加工子程序G代码。 利用以上得到的G代码对实际凸轮片进行实际磨削,同时还要利用凸轮检测仪监测加工的 凸轮,检测出误差Δx,并利用公式(17)对升程进行补偿,最后再由计算机自动得出数控 加工G代码,并利用所得的G代码对凸轮进行磨削,直至磨削结果获得满意即可。

式中,α为补偿系数,一般取为0.5—1.0。

10)实际磨削过程结束之后,可得出全程的轮廓误差图(见图7)和相邻误差图(见图 8)。结果表明经过速度优化及动态补偿后,实际磨削后得到的相邻误差小于0.004mm,小于 国家标准(≤0.005mm);轮廓误差均保持在0.012mm内,也已远远超越了国家标准 (≤0.025mm)。

按照上述步骤,即可利用经典公式对凸轮磨削的速度进行优化,从而实现高速、高精度 磨削的目的。

表1:某型号数控磨床提供的凸轮升程表数据(mm)

角度 升程 角度 升程 角度 升程 角度 升程 1 0.0000 63 17.0000 125 12.4640 187 3.6545 2 0.0069 64 17.0000 126 12.3230 188 3.5397 3 0.0276 65 17.0000 127 12.1820 189 3.4265 4 0.0622 66 17.0000 128 12.0400 190 3.3149 5 0.1108 67 17.0000 129 11.8960 191 3.2048 6 0.1735 68 17.0000 130 11.7520 192 3.0965 7 0.2501 69 16.9980 131 11.6070 193 2.9898 8 0.3422 70 16.9940 132 11.4610 194 2.8847 9 0.4487 71 16.9860 133 11.3140 195 2.7814 10 0.5705 72 16.9740 134 11.1660 196 2.6798 11 0.7080 73 16.9610 135 11.0180 197 2.5798 12 0.8616 74 16.9430 136 10.8690 198 2.4817 13 1.0320 75 16.9230 137 10.7190 199 2.3853 14 1.2196 76 16.8990 138 10.5700 200 2.2907 15 1.4253 77 16.8720 139 10.4190 201 2.1978 16 1.6499 78 16.8430 140 10.2680 202 2.1068 17 1.8942 79 16.8090 141 10.1170 203 2.0175 18 2.1592 80 16.7740 142 9.9659 204 1.9301 19 2.4462 81 16.7350 143 9.8143 205 1.8446 20 2.7564 82 16.6930 144 9.6626 206 1.7608 21 3.0914 83 16.6470 145 9.5107 207 1.6790 22 3.4527 84 16.5990 146 9.3588 208 1.5990 23 3.8424 85 16.5480 147 9.2069 209 1.5208 24 4.2626 86 16.4940 148 9.0550 210 1.4446 25 4.7052 87 16.4370 149 8.9032 211 1.3703 26 5.1508 88 16.3770 150 8.7516 212 1.2978 27 5.5986 89 16.3150 151 8.6000 213 1.2273 28 6.0484 90 16.2490 152 8.4490 214 1.1587 29 6.5000 91 16.1810 153 8.2981 215 1.0920 30 6.9532 92 16.1100 154 8.1476 216 1.0271

31 7.4077 93 16.0360 155 7.9975 217 0.9644 32 7.8632 94 15.9590 156 7.8478 218 0.9035 33 8.3194 95 15.8800 157 7.6986 219 0.8446 34 8.7761 96 15.7980 158 7.5499 220 0.7876 35 9.2331 97 15.7130 159 7.4019 221 0.7326 36 9.6900 98 15.6260 160 7.2545 222 0.6795 37 10.1470 99 15.5360 161 7.1078 223 0.6284 38 10.6030 100 15.4440 162 6.9617 224 0.5793 39 11.0580 101 15.3490 163 6.8165 225 0.5322 40 11.5120 102 15.2520 164 6.6720 226 0.4870 41 11.9650 103 15.1530 165 6.5284 227 0.4439 42 12.4160 104 15.0510 166 6.3856 228 0.4027 43 12.8650 105 14.9470 167 6.2438 229 0.3635 44 13.3120 106 14.8410 168 6.1030 230 0.3263 45 13.7470 107 14.7320 169 5.9631 231 0.2911 46 14.1530 108 14.6220 170 5.8243 232 0.2579 47 14.5290 109 14.5090 171 5.6865 233 0.2267 48 14.8790 110 14.3940 172 5.5499 234 0.1975 49 15.1990 111 14.2780 173 5.4144 235 0.1703 50 15.4930 112 14.1590 174 5.2801 236 0.1452 51 15.7600 113 14.0380 175 5.1470 237 0.1220 52 16.0000 114 13.9160 176 5.0151 238 0.1008 53 16.2140 115 13.7920 177 4.8845 239 0.0817 54 16.4020 116 13.6660 178 4.7552 240 0.0645 55 16.5640 117 13.5380 179 4.6272 241 0.0494 56 16.7000 118 13.4090 180 4.5006 242 0.0363 57 16.8110 119 13.2780 181 4.3753 243 0.0252 58 16.8960 120 13.1460 182 4.2515 244 0.0161 59 16.9560 121 13.0120 183 4.1292 245 0.0091 60 16.9900 122 12.8770 184 4.0082 246 0.0040 61 17.0000 123 12.7410 185 3.8888 247 0.0010 62 17.0000 124 12.6030 186 3.7709 248 0.0000

注:由于凸轮旋转角在249-360度之间,升程均为0mm,故不在表中列出。

去获取专利,查看全文>

相似文献

  • 专利
  • 中文文献
  • 外文文献
获取专利

客服邮箱:kefu@zhangqiaokeyan.com

京公网安备:11010802029741号 ICP备案号:京ICP备15016152号-6 六维联合信息科技 (北京) 有限公司©版权所有
  • 客服微信

  • 服务号