CN108681631B - 圆柱齿轮插齿加工切削力预测的方法 - Google Patents

圆柱齿轮插齿加工切削力预测的方法 Download PDF

Info

Publication number
CN108681631B
CN108681631B CN201810433979.XA CN201810433979A CN108681631B CN 108681631 B CN108681631 B CN 108681631B CN 201810433979 A CN201810433979 A CN 201810433979A CN 108681631 B CN108681631 B CN 108681631B
Authority
CN
China
Prior art keywords
cutter
gear
point
points
pinion
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
CN201810433979.XA
Other languages
English (en)
Other versions
CN108681631A (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.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
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 Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN201810433979.XA priority Critical patent/CN108681631B/zh
Publication of CN108681631A publication Critical patent/CN108681631A/zh
Application granted granted Critical
Publication of CN108681631B publication Critical patent/CN108681631B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • B23F—MAKING GEARS OR TOOTHED RACKS
    • B23F1/00—Making gear teeth by tools of which the profile matches the profile of the required surface
    • B23F1/04—Making gear teeth by tools of which the profile matches the profile of the required surface by planing or slotting

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computational Mathematics (AREA)
  • Mechanical Engineering (AREA)
  • Gear Processing (AREA)

Abstract

本发明涉及一种圆柱齿轮插齿加工切削力预测的方法，包括以下步骤：计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点；根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵；将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角的大小为判据，筛选出工件在插齿刀内的点；在插齿刀的齿顶圆和渐开线上等距离分布足够多的数据点；插齿刀每一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点；基于单位切削力的公式可求出切削力的大小。本发明能够准确预测圆柱齿轮插齿加工切削力，具有计算精度高、计算效率高等优点。

Description

圆柱齿轮插齿加工切削力预测的方法
技术领域
本发明涉及圆柱齿轮插齿加工领域，更具体地说，涉及一种圆柱齿轮插齿加工切削力预测的方法。
背景技术
齿轮是动力和运动传递的机械关键基础零件。圆柱齿轮是最为普遍的一种齿轮样式，广泛应用于航空齿轮箱、高铁齿轮箱、风电齿轮箱、高档汽车变速器等高端装备。插齿加工是圆柱齿轮重要的加工制造方法，插齿加工过程中切削力对插齿刀具寿命、加工效率、加工质量等具有重要影响，如何精确预测切削力是圆柱齿轮插齿加工需要解决的一个关键技术难题。
目前，切削力的预测方法主要有切削力经验公式、机械学说模型和理论切削力模型三种，以上三种方法难以精确预测圆柱齿轮插齿加工中的切削力。
发明内容
本发明要解决的技术问题在于，提供一种可以精确预测圆柱齿轮插齿加工切削力的方法。
本发明解决其技术问题所采用的技术方案是：构造一种圆柱齿轮插齿加工切削力预测的方法，包括以下步骤：
步骤一、计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点，根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵，将圆环内的点经过坐标变换矩阵变换为插齿刀坐标系中的点；
步骤二、插齿刀每进行一步切削时，将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角大小为判据，筛选出工件在插齿刀内的点，点的数量即为切屑面积的大小a；
步骤三、在插齿刀的齿顶圆和左、右渐开线上等距离分布数据点；
步骤四、插齿刀每进行一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点，点的数量即为接触长度的大小b；
步骤五，基于单位切削力公式：切削方向Ft＝Ktca+Kteb，径向方向Ff＝Kfca+Kfeb，轴向方向Fr＝Krca+Kreb，计算圆柱齿轮插齿加工切削力的大小，其中Ktc、Kfc、Krc为切削系数，Kte、Kfe、Kre为边缘系数。
上述方案中，步骤一中所述圆柱齿轮与插齿刀加工时的接触区域圆环的表达式为：
Figure BDA0001654153380000021
Figure BDA0001654153380000022
其中，x1、y1分别表示齿顶圆的横、纵坐标，x2、y2分别表示齿根圆的横、纵坐标；x表示变化角，r1、r2分别表示工件齿顶圆、齿根圆的半径。
上述方案中，在公式(1)、公式(2)接触区域内分布数据点的坐标方程式为：
Figure BDA0001654153380000023
其中，r表示极径，i表示第i条等距圆环，通过公式(3)在每条等距曲线上均匀取j个点得到点集Nij(xij,yij)。
上述方案中，步骤一中圆柱齿轮固定在坐标系O0-x0y0z0上，坐标系O1-x1y1z1与插齿刀中心固定，其x轴通过插齿刀分度圆与圆柱齿轮节曲线的切点，指向节曲线在该点的法线反方向，坐标系Os-xsyszs固定在插齿刀上，相对于坐标系O1-x1y1z1转过的角度为θ，插齿刀中间坐标系O1-x1y1z1到圆柱齿轮坐标系O0-x0y0z0下的坐标变换矩阵为M01，插齿刀固连坐标系Os-xsyszs到中间坐标系O1-x1y1z1的坐标变换矩阵为M1s，圆柱齿轮坐标系O0-x0y0z0到插齿刀坐标系Os-xsyszs的坐标变换矩阵为：
Mn＝(M01M1s)′＝[(2ty-nx-ny)*cosθ+(2tx+nx+ny)*sinθ+xn1+yn1+1]′ (4)
其中nx,ny为圆柱齿轮节曲线单位法向量的分量，tx,ty为圆柱齿轮节曲线单位切向量的分量，xn1,yn1为坐标系O1-x1y1z1任一一点的坐标。
上述方案中，步骤二中将圆柱齿轮上的等距曲线上的点Nij(xij,yij)通过所述坐标变换矩阵Mn转化成插齿刀坐标系下的点N'ij(xsij,ysij)，并求出N'ij(xsij,ysij)点在插齿刀坐标系下的极角
Figure BDA0001654153380000031
与极径
Figure BDA0001654153380000032
上述方案中，所述步骤二中的判据为：
Figure BDA0001654153380000033
pij≤r (6)
其中，sita表示插齿刀每个齿对应的圆心角，ct2表示齿顶圆上渐开线上展角，r表示为对应点与插齿刀圆心连线与插齿刀齿廓交点的距离。
上述方案中，步骤二中实际插齿加工过程中，上一步切除的点将不会在下一步切削加工中出现，所以还需要将原有的圆柱齿轮等距线上的点集Nij(xij,yij)中剔除上一步切削的点，被切除的点集V经过矩阵Mn变换获得，为了使原点集Nij(xij,yij)的点能够剔除掉，因此还需要将矩阵V经过逆矩阵Mn′变换得到V1，再从原点集Nij(xij,yij)中剔除V1，具体的剔除程序如下：
Figure BDA0001654153380000041
其中，vpa(V1,5)表示在矩阵V1中每个元素精确到小数点5位，setdiff函数用来求矩阵V1与矩阵Nij中不同的元素；
通过以上判据公式(5)、(6)和剔除程序公式(7)筛选的点数为n，圆环内每一点代表的面积为s，计算该步圆柱齿轮插齿加工筛选的点的面积为a＝n*s；
如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0；重复上述步骤，计算圆柱齿轮插齿加工每一步筛选的点的面积，直到插齿刀加工完圆柱齿轮为止。
上述方案中，步骤三中所述插齿刀的齿顶圆等距离分布点的坐标为：
Figure BDA0001654153380000042
其中fail表示插齿刀齿顶圆的范围角，ra表示插齿刀齿顶圆的半径；由于插齿刀齿顶圆上的点需要进行矩阵M1S变换，公式(8)中的点变换后得到插齿刀齿顶圆上的横纵坐标(x3,y3)为：
Figure BDA0001654153380000043
步骤三中所述插齿刀的左渐开线上等距离分布点的坐标为：
Figure BDA0001654153380000044
由于插齿刀左渐开线上的点需要进行矩阵M1S变换，公式(10)中的点变换后得到插齿刀左渐开线上点的横纵坐标(x4,y4)为：
Figure BDA0001654153380000045
其中，i1表示插齿刀径向方向上与左渐开线的交点范围，p＝V1(i2,1)表示插齿刀左渐开线上点的极角，i2＝length(V1)表示矩阵V1的长度；用矩阵V1把左渐开线上分布点的极角储存起来，表示为V1＝[V1；bi]；其中，插齿刀左渐开线上分布点在坐标系上的极角为bi＝(i3-0.25)*sita+ct-ct1；插齿刀左渐开线上分布点在左渐开线上的展角为ct＝tan(a1)-a1；插齿刀等距圆与左渐开线上的交点的压力角为a1＝acos(rb/i1)。
上述方案中，步骤三中所述插齿刀的右渐开线上等距离分布点的坐标为：
Figure BDA0001654153380000051
公式(12)中的点经过矩阵M1S变换后得到插齿刀右渐开线上点的横纵坐标(x5,y5)为：
Figure BDA0001654153380000052
其中，p1表示插齿刀右渐开线上点的极角；
将公式(9)、(11)、(13)中的点集也即插齿刀齿顶圆、左右渐开线上分布的均匀点储存在矩阵V2中，并用(x6,y6)表示插齿刀分布的点集的横纵坐标；矩阵V2中的点到工件圆心的距离用p1表示为
p1＝(x6 2+y6 2)(1/2) (14)
根据插齿刀齿顶圆、左右渐开线的长度均匀分布点，总长度为s1，均匀分布的点数为n1，每个点代表的长度为d＝s1/n1。
上述方案中，步骤四中所述判据为：
p1≤ra (15)
通过公式(14)、(15)中的判据筛选的点数为n2，计算该步圆柱齿轮插齿加工接触长度的大小为b＝n2*d；
如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0；重复上述步骤，计算圆柱齿轮插齿加工每一步接触长度，直到插齿刀加工完圆柱齿轮为止。
实施本发明的圆柱齿轮插齿加工切削力预测的方法，具有以下有益效果：
本发明基于插齿加工包络原理，利用判据判断得出插齿加工圆柱齿轮每一步的切屑面积和接触长度，进而预测求出切削力的大小，具有计算精度高、效率高等优点。
附图说明
下面将结合附图及实施例对本发明作进一步说明，附图中：
图1为插齿刀加工圆柱齿轮等距曲线的某一瞬时；
图2为圆柱齿轮与插齿刀坐标位置关系示意图；
图3为插齿刀切削工件计算切屑面积时第一步筛选点的示意图；
图4为插齿刀切削工件计算切屑面积时第二步筛选点的示意图；
图5为插齿刀切削工件第二次径向进给筛选点的示意图；
图6为在插齿刀齿廓上分布均匀点集的模型图；
图7为插齿刀切削工件计算接触长度时第一步筛选点的示意图；
图8为插齿刀切削工件计算接触长度时第二次径向进给筛选点的示意图；
图9a-9c为插齿刀切削工件预测的切削方向、径向方向、轴向方向上的切削力图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解，现对照附图详细说明本发明的具体实施方式。
在此以已知的圆柱齿轮传动为例，按照上述方法预测出插齿加工切削力的方法。齿轮的基本参数如表1所示：
表1圆柱齿轮基本参数
Figure BDA0001654153380000071
圆柱齿轮节曲线弧长
Figure BDA0001654153380000072
插齿刀相对转角
Figure BDA0001654153380000073
R为插齿刀分度圆半径。
插齿刀选用标准渐开线插刀，渐开线插齿刀基本参数为：模数m与计算出来的圆柱齿轮模数一致，齿数Z0＝18，压力角α，节圆齿厚S，节圆半径R和基圆半径Rb。
请参见图1-图9c，本发明方法具体步骤如下：
步骤一，计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点，根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵，将圆环内的点经过坐标变换矩阵变换为插齿刀坐标系中的点。
步骤二，插齿刀每一步切削时，将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角大小为判据，筛选出工件在插齿刀内的点，点的数量即为切屑面积的大小a。
步骤三，在插齿刀的齿顶圆和左、右渐开线上等距离分布足够多的数据点。
步骤四，插齿刀每一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点，点的数量即为接触长度的大小b。
步骤五，基于单位切削力公式Ft＝Ktca+Kteb，Ff＝Kfca+Kfeb，Fr＝Krca+Kreb，计算圆柱齿轮插齿加工切削力的大小。
步骤一中圆柱齿轮与插齿刀加工时的接触区域圆环的表达式为：
Figure BDA0001654153380000081
Figure BDA0001654153380000082
其中，x1、y1表示齿顶圆的横、纵坐标，x2、y2表示齿根圆的横、纵坐标。x表示变化角。
步骤一中在公式(1)、(2)接触区域内分布数据点的坐标方程式为：
Figure BDA0001654153380000083
其中，r表示极径，i表示第i条等距圆环，通过公式(3)在每条等距曲线上均匀取j个点得到点集Nij(xij,yij)。
步骤一中圆柱齿轮固定在坐标系O0-x0y0z0上。坐标系O1-x1y1z1与插齿刀中心固定，其x轴通过插齿刀分度圆与圆柱齿轮节曲线的切点，指向节曲线在该点的法线反方向。坐标系Os-xsyszs固定在插齿刀上，相对于坐标系O1-x1y1z1转过的角度为θ。插齿刀中间坐标系O1-x1y1z1到圆柱齿轮坐标系O0-x0y0z0下的坐标变换矩阵为M01，插齿刀固连坐标系Os-xsyszs到中间坐标系O1-x1y1z1的坐标变换矩阵为M1s。圆柱齿轮坐标系O0-x0y0z0到插齿刀坐标系Os-xsyszs的坐标变换矩阵为：
Mn＝(M01M1s)′＝[(2ty-nx-ny)*cosθ+(2tx+nx+ny)*sinθ+xn1+yn1+1]′ (4)
其中nx,ny为圆柱齿轮节曲线单位法向量的分量，tx,ty为圆柱齿轮节曲线单位切向量的分量，xn1,yn1为坐标系O1-x1y1z1任一一点的坐标。
步骤二中将圆柱齿轮上的等距曲线上的点Nij(xij,yij)通过上述坐标变换矩阵Mn转化成插齿刀坐标系下的点N'ij(xsij,ysij)，并求出N'ij(xsij,ysij)点在插齿刀坐标系下的极角
Figure BDA0001654153380000091
与极径
Figure BDA0001654153380000092
步骤二中的判据具体为：
Figure BDA0001654153380000093
pij≤r (6)
其中，sita表示插齿刀每个齿对应的圆心角，ct2表示齿顶圆上渐开线上展角，r表示为对应点与插齿刀圆心连线与插齿刀齿廓交点的距离。
步骤二中实际插齿加工过程中，上一步切除的点将不会在下一步切削加工中出现，所以还需要将原有的圆柱齿轮等距线上的点集Nij(xij,yij)中剔除上一步切削的点，被切除的点集V经过矩阵Mn变换获得，为了使原点集Nij(xij,yij)的点能够剔除掉，因此还需要将矩阵V经过逆矩阵Mn′变换得到V1，再从原点集Nij(xij,yij)中剔除V1。具体的剔除程序如下：
Figure BDA0001654153380000094
其中，vpa(V1,5)表示在矩阵V1中每个元素精确到小数点5位，setdiff函数用来求矩阵V1与矩阵Nij中不同的元素。
通过以上判据公式(5)、(6)和剔除程序公式(7)筛选的点数为18291，圆环内每一点代表的面积为0.153mm2，计算该步圆柱齿轮插齿加工筛选的点的面积为a＝0.153n。
如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0。重复上述步骤，计算圆柱齿轮插齿加工每一步筛选的点的面积，直到插齿刀加工完圆柱齿轮为止。
步骤三中插齿刀的齿顶圆等距离分布点的坐标为：
Figure BDA0001654153380000101
其中fail表示插齿刀齿顶圆的范围角，由于插齿刀齿顶圆上的点需要进行矩阵M1S变换，公式(8)中的点变换后得到插齿刀齿顶圆上的横纵坐标(x3,y3)为：
Figure BDA0001654153380000102
步骤三中插齿刀的左渐开线上等距离分布点的坐标为：
Figure BDA0001654153380000103
由于插齿刀左渐开线上的点需要进行矩阵M1S变换，公式(10)中的点变换后得到插齿刀左渐开线上点的横纵坐标(x4,y4)为：
Figure BDA0001654153380000104
其中，i1表示插齿刀径向方向上与左渐开线的交点范围为i1＝32:0.1:40，p＝V1(i2,1)表示插齿刀左渐开线上点的极角，i2＝length(V1)表示矩阵V1的长度。用矩阵V1把左渐开线上分布点的极角储存起来，表示为V1＝[V1；bi]。其中，插齿刀左渐开线上分布点在坐标系上的极角为bi＝(i3-0.25)*sita+ct-ct1；插齿刀左渐开线上分布点在左渐开线上的展角为ct＝tan(a1)-a1；插齿刀等距圆与左渐开线上的交点的压力角为a1＝acos(rb/i1)。
步骤三中所述插齿刀的右渐开线上等距离分布点的坐标为：
Figure BDA0001654153380000105
公式(12)中的点经过矩阵M1S变换后得到插齿刀右渐开线上点的横纵坐标(x5,y5)为：
Figure BDA0001654153380000111
其中，p1表示插齿刀右渐开线上点的极角，计算极角p1的方法与上述左渐开线上点的极角计算方法类似。
将公式(9)、(11)、(13)中的点集也即插齿刀齿顶圆、左右渐开线上分布的均匀点储存在矩阵V2中，并用(x6,y6)表示插齿刀分布的点集的横纵坐标。矩阵V2中的点到工件圆心的距离用p1表示为：
p1＝(x6 2+y6 2)(1/2) (14)
根据插齿刀齿顶圆、左右渐开线的长度均匀分布点，总长度为17.4mm，均匀分布的点数为147，每个点代表的长度为d＝0.1186mm。
步骤四中所述判据为：
p1＜＝54 (15)
通过公式(15)中的判据筛选的点数为n2，计算该步圆柱齿轮插齿加工接触长度的大小为b＝0.1186*n2。
如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0。重复上述步骤，计算圆柱齿轮插齿加工每一步接触长度，直到插齿刀加工完圆柱齿轮为止。
步骤五，基于单位切削力的公式Ft＝Ktca+Kte b，Ff＝Kfc a+Kfe b，Fr＝Krca+Kreb，即可求出切削力的大小。其中切削系数Ktc、Kfc、Krc和边缘系数Kte、Kfe、Kre可以从切削实验中直接通过刀具-工件进行校核得到，都为一常数。经试验得到Ktc＝1636.88、Kfc＝455.48、Krc＝0、Kte＝28、Kfe＝2.3、Kre＝0。
将得到的切屑面积a、接触长度b带入以上公式中，得到预测的切削方向Ft、径向方向Ff、轴向方向Fr上的切削力如下图9a-9c所示。
上面结合附图对本发明的实施例进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，而不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可做出很多形式，这些均属于本发明的保护之内。

Claims (2)

1.一种圆柱齿轮插齿加工切削力预测的方法，其特征在于，包括以下步骤：
步骤一、计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点，根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵，将圆环内的点经过坐标变换矩阵变换为插齿刀坐标系中的点；
步骤二、插齿刀每进行一步切削时，将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角大小为判据，筛选出工件在插齿刀内的点，点的数量即为切屑面积的大小a；
步骤三、在插齿刀的齿顶圆和左、右渐开线上等距离分布数据点；
步骤四、插齿刀每进行一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点，点的数量即为接触长度的大小b；
步骤五，基于单位切削力公式：切削方向Ft＝Ktca+Kteb，径向方向Ff＝Kfca+Kfeb，轴向方向Fr＝Krca+Kreb，计算圆柱齿轮插齿加工切削力的大小，其中Ktc、Kfc、Krc为切削系数，Kte、Kfe、Kre为边缘系数；
步骤一中所述圆柱齿轮与插齿刀加工时的接触区域圆环的表达式为：
Figure FDA0002249795200000011
Figure FDA0002249795200000012
其中，x1、y1分别表示齿顶圆的横、纵坐标，x2、y2分别表示齿根圆的横、纵坐标；x表示变化角，r1、r2分别表示工件齿顶圆、齿根圆的半径；
在公式(1)、公式(2)接触区域内分布数据点的坐标方程式为：
Figure FDA0002249795200000021
其中，r表示极径，i表示第i条等距圆环，通过公式(3)在每条等距曲线上均匀取j个点得到点集Nij(xij,yij)；
步骤一中圆柱齿轮固定在坐标系O0-x0y0z0上，坐标系O1-x1y1z1与插齿刀中心固定，其x轴通过插齿刀分度圆与圆柱齿轮节曲线的切点，指向节曲线在该点的法线反方向，坐标系Os-xsyszs固定在插齿刀上，相对于坐标系O1-x1y1z1转过的角度为θ，插齿刀中间坐标系O1-x1y1z1到圆柱齿轮坐标系O0-x0y0z0下的坐标变换矩阵为M01，插齿刀固连坐标系Os-xsyszs到中间坐标系O1-x1y1z1的坐标变换矩阵为M1s，圆柱齿轮坐标系O0-x0y0z0到插齿刀坐标系Os-xsyszs的坐标变换矩阵为：
Mn＝(M01M1s)′＝[(2ty-nx-ny)*cosθ+(2tx+nx+ny)*sinθ+xn1+yn1+1]′ (4)
其中nx,ny为圆柱齿轮节曲线单位法向量的分量，tx,ty为圆柱齿轮节曲线单位切向量的分量，xn1,yn1为坐标系O1-x1y1z1任意一点的坐标；
步骤二中将圆柱齿轮上的等距曲线上的点Nij(xij,yij)通过所述坐标变换矩阵Mn转化成插齿刀坐标系下的点N'ij(xsij,ysij)，并求出N'ij(xsij,ysij)点在插齿刀坐标系下的极角
Figure FDA0002249795200000022
与极径
Figure FDA0002249795200000023
所述步骤二中的判据为：
Figure FDA0002249795200000024
pij≤r (6)
其中，sita表示插齿刀每个齿对应的圆心角，ct2表示齿顶圆上渐开线上展角，r表示为对应点与插齿刀圆心连线与插齿刀齿廓交点的距离；
步骤二中实际插齿加工过程中，上一步切除的点将不会在下一步切削加工中出现，所以还需要将原有的圆柱齿轮等距线上的点集Nij(xij,yij)中剔除上一步切削的点，被切除的点集V经过矩阵Mn变换获得，为了使原点集Nij(xij,yij)的点能够剔除掉，因此还需要将矩阵V经过逆矩阵Mn′变换得到V1，再从原点集Nij(xij,yij)中剔除V1，具体的剔除程序如下：
Figure FDA0002249795200000031
其中，vpa(V1,5)表示在矩阵V1中每个元素精确到小数点5位，setdiff函数用来求矩阵V1与矩阵Nij中不同的元素；
通过以上判据公式(5)、(6)和剔除程序公式(7)筛选的点数为n，圆环内每一点代表的面积为s，计算该步圆柱齿轮插齿加工筛选的点的面积为a＝n*s；
如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0；重复上述步骤，计算圆柱齿轮插齿加工每一步筛选的点的面积，直到插齿刀加工完圆柱齿轮为止；
步骤三中所述插齿刀的齿顶圆等距离分布点的坐标为：
Figure FDA0002249795200000032
其中fail表示插齿刀齿顶圆的范围角，ra表示插齿刀齿顶圆的半径；由于插齿刀齿顶圆上的点需要进行矩阵M1S变换，公式(8)中的点变换后得到插齿刀齿顶圆上的横纵坐标(x3,y3)为：
Figure FDA0002249795200000041
步骤三中所述插齿刀的左渐开线上等距离分布点的坐标为：
Figure FDA0002249795200000042
由于插齿刀左渐开线上的点需要进行矩阵M1S变换，公式(10)中的点变换后得到插齿刀左渐开线上点的横纵坐标(x4,y4)为：
Figure FDA0002249795200000043
其中，i1表示插齿刀径向方向上与左渐开线的交点范围，p＝V1(i2,1)表示插齿刀左渐开线上点的极角，i2＝length(V1)表示矩阵V1的长度；用矩阵V1把左渐开线上分布点的极角储存起来，表示为V1＝[V1；bi]；其中，插齿刀左渐开线上分布点在坐标系上的极角为bi＝(i3-0.25)*sita+ct-ct1；插齿刀左渐开线上分布点在左渐开线上的展角为ct＝tan(a1)-a1；插齿刀等距圆与左渐开线上的交点的压力角为a1＝acos(rb/i1)；
步骤三中所述插齿刀的右渐开线上等距离分布点的坐标为：
Figure FDA0002249795200000044
公式(12)中的点经过矩阵M1S变换后得到插齿刀右渐开线上点的横纵坐标(x5,y5)为：
Figure FDA0002249795200000045
其中，p1表示插齿刀右渐开线上点的极角；
将公式(9)、(11)、(13)中的点集也即插齿刀齿顶圆、左右渐开线上分布的均匀点储存在矩阵V2中，并用(x6,y6)表示插齿刀分布的点集的横纵坐标；矩阵V2中的点到工件圆心的距离用p1表示为
p1＝(x6 2+y6 2)(1/2) (14)
根据插齿刀齿顶圆、左右渐开线的长度均匀分布点，总长度为s1，均匀分布的点数为n1，每个点代表的长度为d＝s1/n1。
2.根据权利要求1所述的圆柱齿轮插齿加工切削力预测的方法，其特征在于，步骤四中所述判据为：
p1≤ra (15)
通过公式(14)、(15)中的判据筛选的点数为n2，计算该步圆柱齿轮插齿加工接触长度的大小为b＝n2*d；
如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0；重复上述步骤，计算圆柱齿轮插齿加工每一步接触长度，直到插齿刀加工完圆柱齿轮为止。
CN201810433979.XA 2018-05-08 2018-05-08 圆柱齿轮插齿加工切削力预测的方法 Active CN108681631B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810433979.XA CN108681631B (zh) 2018-05-08 2018-05-08 圆柱齿轮插齿加工切削力预测的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810433979.XA CN108681631B (zh) 2018-05-08 2018-05-08 圆柱齿轮插齿加工切削力预测的方法

Publications (2)

Publication Number Publication Date
CN108681631A CN108681631A (zh) 2018-10-19
CN108681631B true CN108681631B (zh) 2020-05-05

Family

ID=63805204

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810433979.XA Active CN108681631B (zh) 2018-05-08 2018-05-08 圆柱齿轮插齿加工切削力预测的方法

Country Status (1)

Country Link
CN (1) CN108681631B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110653403B (zh) * 2019-10-16 2020-11-03 盐城工业职业技术学院 一种圆角铣削每齿切削面积快速预测方法
CN111112755B (zh) * 2019-12-11 2021-01-26 武汉理工大学 一种变传动比齿条插齿力的预测方法
CN112518424B (zh) * 2020-12-01 2021-11-09 湘潭大学 螺纹车削的切削力预测方法和装置
CN114101807B (zh) * 2021-10-28 2022-12-09 西安交通大学 一种齿轮工件端面连续等距倒角切削刀具的刃型设计方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106151452A (zh) * 2016-08-27 2016-11-23 武汉理工大学 一种非圆柱齿轮齿廓的设计方法
WO2017189256A1 (en) * 2016-04-29 2017-11-02 Tri Tool Inc. Portable lathe machine
CN107992939A (zh) * 2017-12-06 2018-05-04 湖北工业大学 基于深度增强学习的等切削力齿轮加工方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017189256A1 (en) * 2016-04-29 2017-11-02 Tri Tool Inc. Portable lathe machine
CN106151452A (zh) * 2016-08-27 2016-11-23 武汉理工大学 一种非圆柱齿轮齿廓的设计方法
CN107992939A (zh) * 2017-12-06 2018-05-04 湖北工业大学 基于深度增强学习的等切削力齿轮加工方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Cutting force prediction in hypoid gear machining;Shinichi SHIRAISHI等;《Journal of Advanced Mechanical Design,Systems,and Manufacturing》;20161231;第10卷(第5期);第1-9页 *
VC++编程计算剃前插齿刀倒棱渐开线的分度圆压力角;廖映华 等;《工具技术》;20050430;第39卷(第4期);第46-48页 *

Also Published As

Publication number Publication date
CN108681631A (zh) 2018-10-19

Similar Documents

Publication Publication Date Title
CN108681631B (zh) 圆柱齿轮插齿加工切削力预测的方法
CN108629120B (zh) 非圆柱齿轮插齿加工切削力的预测方法
CN101774048A (zh) 一种锥齿轮加工方法
US20130171912A1 (en) Method for producing periodic tooth flank modifications, machine tool, and computer-readable medium
CN107942931B (zh) 一种正弦柱面超精密车削螺旋刀具轨迹生成方法
Álvarez et al. Large spiral bevel gears on universal 5-axis milling machines: a complete process
CN109190313A (zh) 一种基于曲面离散化的螺旋锥齿轮齿顶倒棱计算方法
CN113798599A (zh) 一种基于近似蜗杆砂轮的斜齿面齿轮磨削方法
Guo et al. A correction method for power skiving of cylindrical gears lead modification
CN110193708A (zh) 超大直径多段式齿圈加工方法
CN111666640B (zh) 一种考虑纤维切削角与切深变化的碳纤维增强复合材料瞬时铣削力计算方法
CN110587038A (zh) 一种刮齿加工齿廓误差补偿方法
CN114563481B (zh) 一种金属材料结构的损伤监测及在线维修系统
Janßen et al. Validation of the plane-based penetration calculation for gear skiving
Kimme et al. Simulation of error-prone continuous generating production processes of helical gears and the influence on the vibration excitation in gear mesh
CN116956601A (zh) 一种基于端面刃形的小模数斜齿车齿刀具设计方法
CN111112755B (zh) 一种变传动比齿条插齿力的预测方法
CN104708119A (zh) 一种蜗轮副的加工方法
RU2165837C1 (ru) Способ размерного микрошлифования изделий, устройство для его осуществления и приспособление для крепления обрабатываемых изделий
CN115026354B (zh) 一种复杂齿形的车齿刀具逆向包络设计方法
CN108153243A (zh) 用于螺纹铣削的数控机床及螺纹铣加工方法
CN106041224A (zh) 一种螺旋锥齿轮加工方法
CN105798396B (zh) 一种新型的圆柱内齿轮车削方法
CN104722851A (zh) 直齿锥齿轮包络刨削成形方法
CN110102829B (zh) 一种锥齿轮加工工艺的对比方法

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