CN109116803B - 一种剖分式等基圆锥齿轮切齿运动轨迹计算方法 - Google Patents

一种剖分式等基圆锥齿轮切齿运动轨迹计算方法 Download PDF

Info

Publication number
CN109116803B
CN109116803B CN201810911599.2A CN201810911599A CN109116803B CN 109116803 B CN109116803 B CN 109116803B CN 201810911599 A CN201810911599 A CN 201810911599A CN 109116803 B CN109116803 B CN 109116803B
Authority
CN
China
Prior art keywords
coordinate system
coordinate
tool
machine tool
vector
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
Application number
CN201810911599.2A
Other languages
English (en)
Other versions
CN109116803A (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.)
Henan University of Science and Technology
Original Assignee
Henan 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 Henan University of Science and Technology filed Critical Henan University of Science and Technology
Priority to CN201810911599.2A priority Critical patent/CN109116803B/zh
Publication of CN109116803A publication Critical patent/CN109116803A/zh
Application granted granted Critical
Publication of CN109116803B publication Critical patent/CN109116803B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/19Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by positioning or contouring control systems, e.g. to control position from one programmed point to another or to control movement along a programmed continuous path
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/35Nc in input of data, input till input file format
    • G05B2219/35349Display part, programmed locus and tool path, traject, dynamic locus

Landscapes

  • Engineering & Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)
  • Milling Processes (AREA)

Abstract

一种剖分式等基圆锥齿轮切齿运动轨迹计算方法,本方案是基于等基圆曲线齿锥齿轮理论,根据其齿面成形特点,首先建立剖分体切齿加工坐标系,分析齿面每个瞬时的刀具位置姿态,通过刀轴矢量、刀心坐标求解及其坐标变换,得到了轮坯坐标系下的刀位。其次,结合通用机床的结构布局、以及剖分轮坯与机床的安装位置关系,计算切齿加工中每一个瞬时的刀位,再通过后置处理,最终得到了数控切齿时机床运动坐标的表达函数。本发明改变了了在通用数控机床上,当等基圆锥齿轮剖分轮坯的回转中心落在机床之外,无法与机床的某一回转轴线重合时,无法进行加工的问题。

Description

一种剖分式等基圆锥齿轮切齿运动轨迹计算方法
技术领域
本发明涉及一种锥齿轮切齿方法,尤其是涉及一种剖分式等基圆锥齿轮切齿运动轨迹计算方法。
背景技术
特大型锥齿轮是大型重点装备的关键零部件,在国民经济和国防建设中具有不可替代性。目前,大型重点装备中所用的特大型锥齿轮均为直齿锥齿轮。齿轮副的大轮直径大于3000mm以上时,由于其整体轮坯结构特点呈现为环形薄壁件,所以其整体式大轮结构刚性差,加工、运输、装配等环节容易变形,且整体加工常受限于现有制齿设备规格,加工后难以进行公路、铁路运输等,所以,现有特大型直齿锥齿轮的大轮常采用剖分式结构。与曲线齿锥齿轮相比较,直齿锥齿轮承载能力弱、传动平稳性差、冲击大、噪音高,采用剖分式曲线齿锥齿轮替代现有的剖分式直齿锥齿轮,已是行业发展、产业升级的重要方向。
等基圆锥齿轮是我国特有的一种曲线齿锥齿轮,大、小轮均可使用指形铣刀成形加工,所用机床结构简单且多为通用机床,刀具价格低廉,两轴联动即可实现加工,其剖分后单齿加工较易于实现。此外,等基圆锥齿轮螺旋角较小,适于沿齿槽方向对其剖分。因此,剖分式等基圆锥齿轮是替代现有特大型剖分式直齿锥齿轮的理想曲线齿锥齿轮。
现有的等基圆曲线齿锥齿轮切齿理论基于其整体加工而言,加工中齿轮的回转中心必定与机床的某一回转轴重合,这样才能完成切齿运动及切齿时的分度运动。针对特大型等基圆曲线齿锥齿轮,现有机床设备无法加工时,将轮坯进行剖分,然后在较小规格的机床上进行切齿加工时,轮坯剖分体的回转中心不会与机床的某一回转轴重合,原有的切齿运动控制理论显然已经无法适用,并且切齿中的分齿运动已经无法轻松实现。
发明内容
本发明的目的是为解决在通用数控机床上,当等基圆锥齿轮剖分轮坯的回转中心落在机床之外,无法与机床的某一回转轴线重合时,机床无法进行加工的问题,提供一种剖分式等基圆锥齿轮切齿加工方法。
本发明为解决上述技术问题的不足,所采用的技术方案是:
一种剖分式等基圆锥齿轮切齿运动轨迹计算方法,包括如下步骤:
步骤1:将相互啮合的一对特大型等基圆锥齿轮的大轮划分为若干锥齿轮剖分体,建立轮坯坐标系和空间固定坐标系,定义轮坯坐标系
Figure GDA0002671659150000021
轴和空间固定坐标系
Figure GDA0002671659150000022
轴的夹角为e
定义被加工的第k个左旋凹齿面的夹角e=emkl
Figure GDA0002671659150000023
定义被加工的第k个左旋凸齿面的夹角e=enkl
Figure GDA0002671659150000024
定义被加工的第k个右旋凹齿面的夹角e=emkr
Figure GDA0002671659150000025
定义被加工的第k个右旋凸齿面的夹角e=enkr
Figure GDA0002671659150000026
式(01)~(04)中,“±”在沿逆时针方向依次加工齿面时取“+”,“在沿顺时针方向依次切齿时取“-”;
Figure GDA0002671659150000027
Figure GDA0002671659150000028
k=1,2,3...,
z为锥齿轮齿数;
Re为大端锥距;
βe是Re处对应的齿线螺旋角
R是小端锥距到大端锥距之间任意值,
β——对应R处的齿线螺旋角;
Figure GDA0002671659150000031
δi—锥齿轮分锥角,其中i=1或2,当i=1表示相互啮合的特大型锥齿轮中的小轮,当i=2表示相互啮合的特大型锥齿轮中的大轮;
r0—理论齿线与其等距线之距离,由刀具尺寸确定;
s—沿理论齿线法向修形量;
步骤2:推导刀心坐标函数
在固定坐标系下根据刀具加工位置建立对应的刀具坐标下,定义固定坐标系中由原点O指向刀具坐标系的原点Oc的向量为
Figure GDA0002671659150000032
向量
Figure GDA0002671659150000033
的坐标值即为固定坐标系下的刀心坐标,将
Figure GDA0002671659150000034
的坐标通过坐标换算公式转换到轮坯坐标系下,得到在轮坯坐标系下的刀心向量
Figure GDA0002671659150000035
刀心向量
Figure GDA0002671659150000036
的坐标即为轮坯坐标系下刀心坐标,定义刀心向量
Figure GDA0002671659150000037
所述的轮坯坐标系下刀心坐标的换算公式为:
对于第k个左旋凹齿面:
Figure GDA0002671659150000038
对于第k个左旋凸齿面:
Figure GDA0002671659150000039
对于第k个右旋凹齿面:
Figure GDA0002671659150000041
对于第k个右旋凸齿面:
Figure GDA0002671659150000042
式中的emkl、enkl、emkr、enkr,通过步骤1计算获得,
R是小端锥距到大端锥距之间任意值,
β为对应R处的齿线螺旋角;
Figure GDA0002671659150000043
Re为大端锥距;
βe是Re处对应的齿线螺旋角
δi—锥齿轮分锥角,i=1或2,当i=1表示相互啮合的特大型锥齿轮中的小轮,当i=2表示相互啮合的特大型锥齿轮中的大轮;
r0—理论齿线与其等距线之距离;
s—沿理论齿线法向修形量;
步骤3:求解刀轴矢量表达式
定义步骤2中刀具坐标系为Sc:[Oc-ic,jc,kc],因为在刀具坐标系下铣刀始终绕着ic轴回转,因此初始的刀具坐标系刀轴矢量定义为
Figure GDA0002671659150000044
将刀轴矢量
Figure GDA0002671659150000045
从刀具坐标系Sc转换到空间固定坐标系S下,定义空间固定坐标系S刀轴矢量为
Figure GDA0002671659150000051
通过坐标变换可得:
Figure GDA0002671659150000052
其中Moc为刀具坐标系Sc到空间固定坐标系S的坐标变换矩阵,将其代入并运算化简后可得:
Figure GDA0002671659150000053
固定坐标系S下的刀轴矢量
Figure GDA0002671659150000054
转换到轮坯坐标系下刀轴矢量,定义轮坯坐标系下的刀轴矢量为
Figure GDA0002671659150000055
则:
Figure GDA0002671659150000056
为空间固定坐标系S到轮坯坐标系Si的转换矩阵,,计算得出轮坯坐标系的刀轴矢量分别为:
对于第k个左旋凹齿面:
Figure GDA0002671659150000057
对于第k个左旋凸齿面:
Figure GDA0002671659150000058
对于第k个右旋凹齿面:
Figure GDA0002671659150000059
对于第k个右旋凸齿面:
Figure GDA00026716591500000510
步骤4:求解机床直线移动坐标
在步骤三中计算得到的刀轴矢量在轮坯坐标系的三个坐标轴上的投影都不恒为零,根据后置处理理论,采用五坐标联动方式进行加工,根据加工所使用的机床建立机床坐标系取机床坐标系下刀位信息对应的X、Y、Z三个移动坐标与三个转动坐标A、B、C中的任意二个组成的五坐标联动。根据步骤2中确定数控机床的X、Y、Z三个移动坐标在轮坯坐标系的数据信息,在通过后置处理将轮坯坐标系下的刀心坐标转换为机床坐标系下的刀心坐标,根据机床的具体结构确定对应的后置处理算法,定义五坐标联动为X、Y、Z、A、B五坐标,其中转动坐标A、B的运动由刀具的摆动来实现,剖分轮坯在机床工作台上只有平移运动,
建立轮坯与机床的坐标位置关系图,相对于空间固定坐标系,轮坯坐标系是转动的,因此固定坐标系到轮坯坐标系的转换矩阵会根据转动角度e的取值不同而发生改变,定义e的初值为固定值e0
e0=θc/sinδi (14)
对应的转换矩阵设为Moi,Moi是矩阵Mio中的参数取e0时的逆矩阵,
此处的θc=θ±θd,“+”用于凹齿面,“-”用于凸齿面;θ是刀具处在R时对应的刀具中心极角
θd是冠伦平面上理论齿线极径与刀具中心轨迹极径的夹角
轮坯坐标系Si到过渡坐标系So的齐次坐标变换矩阵为:
Figure GDA0002671659150000061
过渡坐标系So变换到工件安装坐标系Sw的坐标转换矩阵为:
Figure GDA0002671659150000062
工件安装坐标系Sw到机床坐标系Sm的坐标变换矩阵为:
Figure GDA0002671659150000071
已知轮坯坐标系下的刀心向量
Figure GDA0002671659150000072
将其从轮坯坐标系经过渡坐标系、工件安装坐标系,最终转化到机床坐标系,设机床坐标系下,刀心向量为
Figure GDA0002671659150000073
则:
Figure GDA0002671659150000074
计算可得:
Figure GDA0002671659150000075
其中,(xi,yi,zi)的取值由步骤2中的式(05)至(08)确定,根据不同情况将(xi,yi,zi)代入式(18)中,可分别求得对应于左旋凹齿面、左旋凸齿面、右旋凹齿面、右旋凸齿面的机床直线移动坐标X、Y、Z;
步骤5:机床回转运动求解
由步骤3中的式(10)到式(13),在轮坯坐标系下求得刀轴矢量
Figure GDA0002671659150000076
根据后置处理原理,将其当作自由矢量处理,只回转不平移,即可求得机床的回转角度,据此,先将其从轮坯坐标系下转换至过渡坐标系下,然后将其转换至工件的安装坐标系下,但此时只回转不平移,所以此时式中的Δz取零。刀轴矢量转换至此,由于工件安装坐标系各轴和机床坐标系各对应轴平行,所以不必继续变换到机床坐标系下,据此就可以求解机床的回转角度,定义工件安装坐标系下的刀轴矢量为
Figure GDA0002671659150000077
则:
Figure GDA0002671659150000078
将各式代入上式,化简可得:
Figure GDA0002671659150000079
根据刀轴矢量的坐标通过后置处理原理计算A、B角,定义A、B角为剖分式数控加工中机床旋转轴的实际回转角度,
将刀轴矢量的起点移动到工件安装坐标系的原点Ow,然后将刀轴矢量绕工件安装坐标系的Yw轴顺时针旋转到YwOwZw平面上,旋转角度为B;再将刀轴矢量绕工件安装坐标系Xw轴顺时针旋转到与Zw轴方向一致,旋转角度为A,其中A、B角的正负号根据右手螺旋定则来确定,
计算获得A角如下:
Figure GDA0002671659150000081
B角如下:
Figure GDA0002671659150000082
当azw=0时,令
Figure GDA0002671659150000083
至此,加工剖分式等基圆锥齿轮时机床刀具的回转运动A、B求解完成,结合步骤4确定机床直线移动坐标X、Y、Z,依据五坐标联动方式加工各个锥齿轮剖分体上的齿面,根据X、Y、Z和A、B获得刀具的运动轨迹。
本方案是基于等基圆曲线齿锥齿轮理论,根据其齿面成形特点,首先建立剖分体切齿加工坐标系,分析齿面每个瞬时的刀具位置姿态,通过刀轴矢量、刀心坐标求解及其坐标变换,得到了轮坯坐标系下的刀位。其次,结合通用机床的结构布局、以及剖分轮坯与机床的安装位置关系,计算切齿加工中每一个瞬时的刀位,再通过后置处理,最终得到了数控切齿时机床运动坐标的表达函数。
本发明的有益效果是:本发明改变了了在在通用数控机床上,当等基圆锥齿轮剖分轮坯的回转中心落在机床之外,无法与机床的某一回转轴线重合时,通过本方法可以利用X、Y、Z、A、B五坐标运动计算数学模型,可精确计算机床运动轨迹,控制刀具姿态,实现剖分式等基圆曲线齿锥齿轮的自由切齿加工。
附图说明
图1为加工坐标系。
图2为轮坯与机床的坐标关系图。
图3为刀轴矢量回转示意图。
图4为算例齿轮的刀具与齿面的位置关系三维绘图。
具体实施方式
图中所示,具体实施方式如下:
下面以表1所示的齿轮副中大轮的加工为例,结合附图对本发明一种剖分式等基圆锥齿轮切齿运动轨迹计算方法作进一步详细说明。本实施例中大轮剖分体的加工齿数为3个。
表1齿轮副基本几何参数
Figure GDA0002671659150000091
步骤1:求解不同加工顺序、不同齿面时的夹角e
在等基圆锥齿轮切齿坐标系中(附图1),轮坯坐标系的
Figure GDA0002671659150000092
轴与空间固定坐标系
Figure GDA0002671659150000093
轴的夹角e的取值决定了每个凸、凹齿面在空间固定坐标下的相对位置关系。针对表1中的齿轮副设计参数,以加工大轮剖分体的不同凹、凸齿面为例,根据等基圆锥齿轮理论推得如下算式。
对于第1个被加工的凹齿面:
Figure GDA0002671659150000101
对于第2个被加工的凹齿面:
Figure GDA0002671659150000102
对于第3个被加工的凹齿面:
Figure GDA0002671659150000103
对于第1个被加工的凸齿面:
Figure GDA0002671659150000104
对于第2个被加工的凸齿面:
Figure GDA0002671659150000105
对于第3个被加工的凸齿面:
Figure GDA0002671659150000106
步骤2:推导刀心坐标函数
由附图1所示的加工坐标系图可知,固定坐标系中的向量
Figure GDA0002671659150000107
由原点O指向刀具坐标系的原点Oc,向量
Figure GDA0002671659150000108
的模长Rc可求解确定。该矢量随着加工的进行,也就是随着坐标原点Oc的不同,会得到不同取值,矢量
Figure GDA0002671659150000109
的一系列坐标值就是整个加工过程中的刀心坐标。将
Figure GDA00026716591500001010
都转换到轮坯坐标系下,就得到缠绕着轮坯的一系列刀心向量
Figure GDA0002671659150000111
设轮坯坐标系下的刀心向量为
Figure GDA0002671659150000112
根据齿轮的旋向及齿面的凹凸性不同,经计算、化简后,轮坯坐标系下的刀心坐标计算公式分别为:
对于大轮剖分体上被加工的第1个凹齿面:
Figure GDA0002671659150000113
对于大轮剖分体上被加工的第2个凹齿面:
Figure GDA0002671659150000114
对于大轮剖分体上被加工的第3个凹齿面:
Figure GDA0002671659150000115
对于大轮剖分体上被加工的第1个凸齿面:
Figure GDA0002671659150000116
对于大轮剖分体上被加工的第2个凸齿面:
Figure GDA0002671659150000117
对于大轮剖分体上被加工的第3个凸齿面:
Figure GDA0002671659150000121
步骤3:求解刀轴矢量表达式
在等基圆锥齿轮加工坐标系中(附图1),刀具坐标系Sc:[Oc-ic,jc,kc]下铣刀始终绕着ic轴回转,因此初始的刀轴矢量可定义为
Figure GDA0002671659150000122
然后,将刀轴矢量
Figure GDA0002671659150000123
从刀具坐标系Sc转换到空间固定坐标系S下,此时的向量定义为
Figure GDA0002671659150000124
通过坐标变换可得:
Figure GDA0002671659150000125
其中Moc为刀具坐标系Sc到空间固定坐标系S的坐标变换矩阵,将其代入并运算化简后可得:
Figure GDA0002671659150000126
然后将向量
Figure GDA0002671659150000127
转换到轮坯坐标系下,设轮坯坐标系下的刀轴矢量为
Figure GDA0002671659150000128
则:
Figure GDA0002671659150000129
Mio为空间固定坐标系S到轮坯坐标系Si的转换矩阵,
Figure GDA00026716591500001210
,将其代入运算、化简后,轮坯坐标系的刀轴矢量分别为:
对于大轮剖分体上被加工的第1个凹齿面:
Figure GDA00026716591500001211
对于大轮剖分体上被加工的第2个凹齿面:
Figure GDA00026716591500001212
对于大轮剖分体上被加工的第3个凹齿面:
Figure GDA0002671659150000131
对于大轮剖分体上被加工的第1个凸齿面:
Figure GDA0002671659150000132
对于大轮剖分体上被加工的第2个凸齿面:
Figure GDA0002671659150000133
对于大轮剖分体上被加工的第3个凸齿面:
Figure GDA0002671659150000134
分析式(35)到式(40)求得的刀轴矢量,发现其在轮坯坐标系的三个坐标轴上的投影都不恒为零,根据后置处理理论,要实现其数控加工,需要五坐标联动。
步骤4:求解机床直线移动坐标
一般来说,五坐标联动是指数控机床的X、Y、Z三个移坐标和绕X、Y、Z轴旋转的三个转动坐标A、B、C中的任意五个坐标的线性插补运动。通常是X、Y、Z与三个转动坐标A、B、C中的任意二个组成的五坐标联动。步骤1、步骤2求得指形铣刀的刀位之后,在后置处理过程中,机床的具体结构不同,后置处理算法不同。不失一般性,以X、Y、Z、A、B五坐标数控机床为例,进行其后置处理计算,并且转动坐标A、B的运动由刀具的摆动来实现,剖分轮坯在机床工作台上只有平移运动。
根据以上分析,建立轮坯与机床的坐标位置关系图(附图2),相对于空间固定坐标系,轮坯坐标系是转动的,因此固定坐标系到轮坯坐标系的转换矩阵会根据转动角度e的取值不同而发生改变。为了确定转动的轮坯坐标系与固定坐标系的关系,也就是需要确定e的大小,当夹角e的取值固定时,轮坯和固定坐标系的关系此时就固定不变了。根据e的初值,求得此时的夹角为e0为:
e0=θc/sinδi (41)
对应的转换矩阵设为Moi,Moi是矩阵Mio中的参数取e0时的逆矩阵。
上述各坐标系之间存在相应的齐次变换关系,其中,轮坯坐标系Si到过渡坐标系So的齐次坐标变换矩阵为:
Figure GDA0002671659150000141
过渡坐标系So变换到工件安装坐标系Sw的坐标转换矩阵为:
Figure GDA0002671659150000142
工件安装坐标系Sw到机床坐标系Sm的坐标变换矩阵为:
Figure GDA0002671659150000143
已知轮坯坐标系下的刀心向量
Figure GDA0002671659150000144
将其从轮坯坐标系经过渡坐标系、工件安装坐标系,最终转化到机床坐标系。设机床坐标系下,刀心向量为
Figure GDA0002671659150000145
则:
Figure GDA0002671659150000146
将各式代入上式,化简可得:
Figure GDA0002671659150000147
其中,(xi,yi,zi)的取值由步骤2中的式(28)至(33)确定,根据不同情况将(xi,yi,zi)代入式(45)中,可分别求得对应于左旋凹齿面、左旋凸齿面、右旋凹齿面、右旋凸齿面的机床直线移动坐标。
至此,指形铣刀加工剖分式等基圆锥齿轮时机床的直线移动坐标求解完成。
步骤5:机床回转运动求解
由步骤3中的式(35)到式(40),在轮坯坐标系下求得刀轴矢量
Figure GDA0002671659150000151
根据后置处理原理,将其当作自由矢量处理,只回转不平移,即可求得机床的回转角度。据此,先将其从图中的轮坯坐标系下转换至过渡坐标系下,然后将其转换至工件的安装坐标系下,但此时只回转不平移,所以此时式中的Δz取零。刀轴矢量转换至此,由于工件安装坐标系各轴和机床坐标系各对应轴平行,所以不必继续变换到机床坐标系下,据此就可以求解机床的回转角度。设工件安装坐标系下的刀轴矢量为
Figure GDA0002671659150000152
则:
Figure GDA0002671659150000153
将各式代入上式,化简可得:
Figure GDA0002671659150000154
在轮坯与机床的坐标位置关系图(附图2)中,A、B角可根据刀轴矢量的坐标,通过后置处理原理计算得到。此时的A、B角就是剖分式数控加工中机床旋转轴的实际回转角度。
如附图3所示,将刀轴矢量的起点移动到工件安装坐标系的原点Ow,然后将刀轴矢量绕Yw轴顺时针旋转到YwOwZw平面上,旋转角度为B;再将刀轴矢量绕Xw轴顺时针旋转到与Zw轴方向一致,旋转角度为A,其中A、B角的正负号根据右手螺旋定则来确定。
结合附图3,由三角函数关系可得:
Figure GDA0002671659150000155
B角的计算公式汇总如下:
Figure GDA0002671659150000161
特殊情况:当azw=0时,令
Figure GDA0002671659150000162
至此,加工剖分式等基圆锥齿轮时机床刀具的回转运动A、B求解完成。
针对本实施例,将表1中的齿轮参数,代入公式(28)-(33),可求得刀心运动轨迹,再将刀心轨迹坐标代入式(45),得到机床直线移动坐标。将表1中的齿轮参数代入公式(35)-(40),求得刀轴矢量,然后将刀轴矢量代入式(47)、式(48),求得机床的回转坐标。根据求得的数据,在matlab环境下,进行实施例齿轮的刀具与齿面的空间位置关系三维绘图,输出结果如附图4所示。由附图4可见,剖分体上每个齿的加工轨迹各不相同,每个瞬时刀具与齿面的相互位置姿态各异。
本发明所列举的技术方案和实施方式并非是限制,与本发明所列举的技术方案和实施方式等同或者效果相同方案都在本发明所保护的范围内。

Claims (1)

1.一种剖分式等基圆锥齿轮切齿运动轨迹计算方法,其特征在于:包括如下步骤:
步骤1:将相互啮合的一对特大型等基圆锥齿轮的大轮划分为若干锥齿轮剖分体,建立轮坯坐标系和空间固定坐标系,定义轮坯坐标系
Figure FDA0002671659140000018
轴和空间固定坐标系
Figure FDA0002671659140000017
轴的夹角为e
定义被加工的第k个左旋凹齿面的夹角e=emkl
Figure FDA0002671659140000011
定义被加工的第k个左旋凸齿面的夹角e=enkl
Figure FDA0002671659140000012
定义被加工的第k个右旋凹齿面的夹角e=emkr
Figure FDA0002671659140000013
定义被加工的第k个右旋凸齿面的夹角e=enkr
Figure FDA0002671659140000014
式(01)~(04)中,“±”在沿逆时针方向依次加工齿面时取“+”,“在沿顺时针方向依次切齿时取“-”;
Figure FDA0002671659140000015
Figure FDA0002671659140000016
k=1,2,3...,
z为锥齿轮齿数;
Re为大端锥距;
βe是Re处对应的齿线螺旋角
R是小端锥距到大端锥距之间任意值,
β——对应R处的齿线螺旋角;
Figure FDA0002671659140000021
δi—锥齿轮分锥角,其中i=1或2,当i=1表示相互啮合的特大型锥齿轮中的小轮,当i=2表示相互啮合的特大型锥齿轮中的大轮;
r0—理论齿线与其等距线之距离,由刀具尺寸确定;
s—沿理论齿线法向修形量;
步骤2:推导刀心坐标函数
在固定坐标系下根据刀具加工位置建立对应的刀具坐标下,定义固定坐标系中由原点O指向刀具坐标系的原点Oc的向量为
Figure FDA0002671659140000022
向量
Figure FDA0002671659140000023
的坐标值即为固定坐标系下的刀心坐标,将
Figure FDA0002671659140000024
的坐标通过坐标换算公式转换到轮坯坐标系下,得到在轮坯坐标系下的刀心向量
Figure FDA0002671659140000025
刀心向量
Figure FDA0002671659140000026
的坐标即为轮坯坐标系下刀心坐标,定义刀心向量
Figure FDA0002671659140000027
所述的轮坯坐标系下刀心坐标的换算公式为:
对于第k个左旋凹齿面:
Figure FDA0002671659140000028
对于第k个左旋凸齿面:
Figure FDA0002671659140000031
对于第k个右旋凹齿面:
Figure FDA0002671659140000032
对于第k个右旋凸齿面:
Figure FDA0002671659140000033
式中的emkl、enkl、emkr、enkr,通过步骤1计算获得,
R是小端锥距到大端锥距之间任意值,
β为对应R处的齿线螺旋角;
Figure FDA0002671659140000034
Re为大端锥距;
βe是Re处对应的齿线螺旋角
δi—锥齿轮分锥角,i=1或2,当i=1表示相互啮合的特大型锥齿轮中的小轮,当i=2表示相互啮合的特大型锥齿轮中的大轮;
r0—理论齿线与其等距线之距离;
s—沿理论齿线法向修形量;
步骤3:求解刀轴矢量表达式
定义步骤2中刀具坐标系为Sc:[Oc-ic,jc,kc],因为在刀具坐标系下铣刀始终绕着ic轴回转,因此初始的刀具坐标系刀轴矢量定义为
Figure FDA0002671659140000041
将刀轴矢量
Figure FDA0002671659140000042
从刀具坐标系Sc转换到空间固定坐标系S下,定义空间固定坐标系S刀轴矢量为
Figure FDA0002671659140000043
通过坐标变换可得:
Figure FDA0002671659140000044
其中Moc为刀具坐标系Sc到空间固定坐标系S的坐标变换矩阵,将其代入并运算化简后可得:
Figure FDA0002671659140000045
固定坐标系S下的刀轴矢量
Figure FDA0002671659140000046
转换到轮坯坐标系下刀轴矢量,定义轮坯坐标系下的刀轴矢量为
Figure FDA0002671659140000047
则:
Figure FDA0002671659140000048
Mio为空间固定坐标系S到轮坯坐标系Si的转换矩阵,,计算得出轮坯坐标系的刀轴矢量分别为:
对于第k个左旋凹齿面:
Figure FDA0002671659140000049
对于第k个左旋凸齿面:
Figure FDA00026716591400000410
对于第k个右旋凹齿面:
Figure FDA0002671659140000051
对于第k个右旋凸齿面:
Figure FDA0002671659140000052
步骤4:求解机床直线移动坐标
在步骤三中计算得到的刀轴矢量在轮坯坐标系的三个坐标轴上的投影都不恒为零,根据后置处理理论,采用五坐标联动方式进行加工,根据加工所使用的机床建立机床坐标系取机床坐标系下刀位信息对应的X、Y、Z三个移动坐标与三个转动坐标A、B、C中的任意二个组成的五坐标联动;根据步骤2中确定数控机床的X、Y、Z三个移动坐标在轮坯坐标系的数据信息,在通过后置处理将轮坯坐标系下的刀心坐标转换为机床坐标系下的刀心坐标,根据机床的具体结构确定对应的后置处理算法,定义五坐标联动为X、Y、Z、A、B五坐标,其中转动坐标A、B的运动由刀具的摆动来实现,剖分轮坯在机床工作台上只有平移运动,
建立轮坯与机床的坐标位置关系图,相对于空间固定坐标系,轮坯坐标系是转动的,因此固定坐标系到轮坯坐标系的转换矩阵会根据转动角度e的取值不同而发生改变,定义e的初值为固定值e0
e0=θc/sinδi (14)
对应的转换矩阵设为Moi,Moi是矩阵Mio中的参数取e0时的逆矩阵,
此处的θc=θ±θd,“+”用于凹齿面,“-”用于凸齿面;θ是刀具处在R时对应的刀具中心极角
θd是冠伦平面上理论齿线极径与刀具中心轨迹极径的夹角
轮坯坐标系Si到过渡坐标系So的齐次坐标变换矩阵为:
Figure FDA0002671659140000061
过渡坐标系So变换到工件安装坐标系Sw的坐标转换矩阵为:
Figure FDA0002671659140000062
工件安装坐标系Sw到机床坐标系Sm的坐标变换矩阵为:
Figure FDA0002671659140000063
已知轮坯坐标系下的刀心向量
Figure FDA0002671659140000064
将其从轮坯坐标系经过渡坐标系、工件安装坐标系,最终转化到机床坐标系,设机床坐标系下,刀心向量为
Figure FDA0002671659140000065
则:
Figure FDA0002671659140000066
计算可得:
Figure FDA0002671659140000067
其中,(xi,yi,zi)的取值由步骤2中的式(05)至(08)确定,根据不同情况将(xi,yi,zi)代入式(18)中,可分别求得对应于左旋凹齿面、左旋凸齿面、右旋凹齿面、右旋凸齿面的机床直线移动坐标X、Y、Z;
步骤5:机床回转运动求解
由步骤3中的式(10)到式(13),在轮坯坐标系下求得刀轴矢量
Figure FDA0002671659140000071
根据后置处理原理,将其当作自由矢量处理,只回转不平移,即可求得机床的回转角度,据此,先将其从轮坯坐标系下转换至过渡坐标系下,然后将其转换至工件的安装坐标系下,但此时只回转不平移,所以此时式中的Δz取零;刀轴矢量转换至此,由于工件安装坐标系各轴和机床坐标系各对应轴平行,定义工件安装坐标系下的刀轴矢量为
Figure FDA0002671659140000072
则:
Figure FDA0002671659140000073
将各式代入上式,化简可得:
Figure FDA0002671659140000074
根据刀轴矢量的坐标通过后置处理原理计算A、B角,定义A、B角为剖分式数控加工中机床旋转轴的实际回转角度,
将刀轴矢量的起点移动到工件安装坐标系的原点Ow,然后将刀轴矢量绕工件安装坐标系的Yw轴顺时针旋转到YwOwZw平面上,旋转角度为B;再将刀轴矢量绕工件安装坐标系Xw轴顺时针旋转到与Zw轴方向一致,旋转角度为A,其中A、B角的正负号根据右手螺旋定则来确定,
计算获得A角如下:
Figure FDA0002671659140000081
B角如下:
Figure FDA0002671659140000082
当azw=0时,令
Figure FDA0002671659140000083
至此,加工剖分式等基圆锥齿轮时机床刀具的回转运动A、B求解完成,结合步骤4确定机床直线移动坐标X、Y、Z,依据五坐标联动方式加工各个锥齿轮剖分体上的齿面,根据X、Y、Z和A、B获得刀具的运动轨迹。
CN201810911599.2A 2018-08-10 2018-08-10 一种剖分式等基圆锥齿轮切齿运动轨迹计算方法 Active CN109116803B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810911599.2A CN109116803B (zh) 2018-08-10 2018-08-10 一种剖分式等基圆锥齿轮切齿运动轨迹计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810911599.2A CN109116803B (zh) 2018-08-10 2018-08-10 一种剖分式等基圆锥齿轮切齿运动轨迹计算方法

Publications (2)

Publication Number Publication Date
CN109116803A CN109116803A (zh) 2019-01-01
CN109116803B true CN109116803B (zh) 2020-11-17

Family

ID=64853022

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810911599.2A Active CN109116803B (zh) 2018-08-10 2018-08-10 一种剖分式等基圆锥齿轮切齿运动轨迹计算方法

Country Status (1)

Country Link
CN (1) CN109116803B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111299716B (zh) * 2019-12-05 2021-03-16 河南科技大学 一种等基圆锥齿轮专用的指形铣刀定位工装以及定位方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101391322A (zh) * 2008-10-30 2009-03-25 吉林大学 球面渐开线齿形弧齿锥齿轮切齿法与机床
CN101391324A (zh) * 2008-10-30 2009-03-25 吉林大学 球面渐开线齿形阿基米德螺线齿锥齿轮切齿法与机床
DE102010026989A1 (de) * 2010-07-13 2012-01-19 Marina Deitert Verfahren und Vorrichtung zur spannenden Fertigung von Bogenverzahnungen an Kegelrädern
CN103418851A (zh) * 2012-04-02 2013-12-04 克林格伦贝格股份公司 采用插齿工艺制造锥形齿轮或准双曲面齿轮的方法
CN103934513A (zh) * 2014-04-04 2014-07-23 西安交通大学 一种摆线齿锥齿轮的加工方法
CN105127519A (zh) * 2015-09-18 2015-12-09 长安大学 一种螺旋锥齿轮的小轮粗切方法
CN108319221A (zh) * 2018-01-24 2018-07-24 东莞长盈精密技术有限公司 刀具轨迹计算方法及中框加工方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2336839A4 (en) * 2008-09-16 2013-11-06 Shin Nippon Koki Co Ltd NUMERICAL CONTROL

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101391322A (zh) * 2008-10-30 2009-03-25 吉林大学 球面渐开线齿形弧齿锥齿轮切齿法与机床
CN101391324A (zh) * 2008-10-30 2009-03-25 吉林大学 球面渐开线齿形阿基米德螺线齿锥齿轮切齿法与机床
DE102010026989A1 (de) * 2010-07-13 2012-01-19 Marina Deitert Verfahren und Vorrichtung zur spannenden Fertigung von Bogenverzahnungen an Kegelrädern
CN103418851A (zh) * 2012-04-02 2013-12-04 克林格伦贝格股份公司 采用插齿工艺制造锥形齿轮或准双曲面齿轮的方法
CN103934513A (zh) * 2014-04-04 2014-07-23 西安交通大学 一种摆线齿锥齿轮的加工方法
CN105127519A (zh) * 2015-09-18 2015-12-09 长安大学 一种螺旋锥齿轮的小轮粗切方法
CN108319221A (zh) * 2018-01-24 2018-07-24 东莞长盈精密技术有限公司 刀具轨迹计算方法及中框加工方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Digitized Modeling and Machining Simulation of Equal Base Bevel Gear;Bin Wang 等;《ITM Web Conf.》;20180202;第17卷;1-8 *
基于数控技术的长锥距弧齿锥齿轮的加工方法研究;姚荣达;《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》;20160215(第2期);B022-719 *
螺旋锥齿轮数控加工及精度控制关键技术研究;第6期;《中国博士学位论文全文数据库工程科技Ⅰ辑》;20160615;B022-40 *

Also Published As

Publication number Publication date
CN109116803A (zh) 2019-01-01

Similar Documents

Publication Publication Date Title
CN111008441B (zh) 整体平头立铣刀端齿直纹型后刀面的磨削轨迹求解方法
CN104500654B (zh) 一种减变速一体化面齿轮副及其加工方法
CN112541235B (zh) 一种准双曲面齿轮副通用设计方法
CN109604738B (zh) 一种基于尼曼蜗轮数学模型的高效侧刃精铣加工方法
CN106312850A (zh) 螺杆磨削刀具的设计方法
Mohan et al. Simulation of whirling process and tool profiling for machining of worms
CN107273612A (zh) 螺旋锥齿轮凹模齿面设计方法及螺旋锥齿轮加工方法
Lu et al. Digitization modeling and CNC machining for enveloping surface parts
CN111715947B (zh) 一种线接触渐缩齿弧齿锥齿轮副成形方法
CN110193708A (zh) 超大直径多段式齿圈加工方法
CN112705794A (zh) 一种用于摆线轮加工的剐齿刀具及其设计方法
Shen et al. Research on shaving processing of spiroid face gear
CN109116803B (zh) 一种剖分式等基圆锥齿轮切齿运动轨迹计算方法
CN111967096A (zh) 金刚滚轮及蜗杆砂轮的设计方法
CN106392199A (zh) 一种汽车转向器摇臂轴齿扇齿廓修形方法
CN109396567B (zh) 修形摆线齿轮展成磨削中确定蜗杆砂轮廓形的数字包络法
Wang et al. The mathematical model of spiral bevel gears-A review
Zhang et al. Method of processing and an analysis of meshing and contact of circular arc tooth trace cylindrical gears
CN110802280B (zh) 一种渐开线螺旋锥齿轮齿面设计方法
Lin et al. A five-axis CNC machining method of orthogonal variable transmission ratio face gear
Zheng et al. Avoidance of cutter retracting interference in noncircular gear shaping through 4-linkage model
CN114769737A (zh) 一种少齿数齿轮成型磨削加工方法
CN103934727B (zh) 一种等角螺旋线分区变参数控精密加工方法
Qiu et al. A processing method for rough cutting hourglass worm on general NC lathe
Zhengyang et al. Research on machining nonorthogonal face gears by power skiving with tooth flank modification based on a six-axis machine tool

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