圆柱齿轮插齿加工切削力预测的方法技术

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

【技术实现步骤摘要】
圆柱齿轮插齿加工切削力预测的方法
本专利技术涉及圆柱齿轮插齿加工领域，更具体地说，涉及一种圆柱齿轮插齿加工切削力预测的方法。
技术介绍
齿轮是动力和运动传递的机械关键基础零件。圆柱齿轮是最为普遍的一种齿轮样式，广泛应用于航空齿轮箱、高铁齿轮箱、风电齿轮箱、高档汽车变速器等高端装备。插齿加工是圆柱齿轮重要的加工制造方法，插齿加工过程中切削力对插齿刀具寿命、加工效率、加工质量等具有重要影响，如何精确预测切削力是圆柱齿轮插齿加工需要解决的一个关键技术难题。目前，切削力的预测方法主要有切削力经验公式、机械学说模型和理论切削力模型三种，以上三种方法难以精确预测圆柱齿轮插齿加工中的切削力。
技术实现思路
本专利技术要解决的技术问题在于，提供一种可以精确预测圆柱齿轮插齿加工切削力的方法。本专利技术解决其技术问题所采用的技术方案是：构造一种圆柱齿轮插齿加工切削力预测的方法，包括以下步骤：步骤一、计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点，根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵，将圆环内的点经过坐标变换矩阵变换为插齿刀坐标系中的点；步骤二、插齿刀每进行一步切削时，将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角大小为判据，筛选出工件在插齿刀内的点，点的数量即为切屑面积的大小a；步骤三、在插齿刀的齿顶圆和左、右渐开线上等距离分布数据点；步骤四、插齿刀每进行一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点，点的数量即为接触长度的大小b；步骤五，基于单位切削力公式：切削方向Ft＝Ktca+Kteb，径向方向Ff＝Kfca+Kfeb，轴向方向Fr＝Krca+Kreb，计算圆柱齿轮插齿加工切削力的大小，其中Ktc、Kfc、Krc为切削系数，Kte、Kfe、Kre为边缘系数。上述方案中，步骤一中所述圆柱齿轮与插齿刀加工时的接触区域圆环的表达式为：其中，x1、y1分别表示齿顶圆的横、纵坐标，x2、y2分别表示齿根圆的横、纵坐标；x表示变化角，r1、r2分别表示工件齿顶圆、齿根圆的半径。上述方案中，在公式(1)、公式(2)接触区域内分布数据点的坐标方程式为：其中，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)点在插齿刀坐标系下的极角与极径上述方案中，所述步骤二中的判据为：pij≤r(6)其中，sita表示插齿刀每个齿对应的圆心角，ct2表示齿顶圆上渐开线上展角，r表示为对应点与插齿刀圆心连线与插齿刀齿廓交点的距离。上述方案中，步骤二中实际插齿加工过程中，上一步切除的点将不会在下一步切削加工中出现，所以还需要将原有的圆柱齿轮等距线上的点集Nij(xij,yij)中剔除上一步切削的点，被切除的点集V经过矩阵Mn变换获得，为了使原点集Nij(xij,yij)的点能够剔除掉，因此还需要将矩阵V经过逆矩阵Mn′变换得到V1，再从原点集Nij(xij,yij)中剔除V1，具体的剔除程序如下：其中，vpa(V1,5)表示在矩阵V1中每个元素精确到小数点5位，setdiff函数用来求矩阵V1与矩阵Nij中不同的元素；通过以上判据公式(5)、(6)和剔除程序公式(7)筛选的点数为n，圆环内每一点代表的面积为s，计算该步圆柱齿轮插齿加工筛选的点的面积为a＝n*s；如果插齿刀沿工件的旋转的角度x≥2π，再进行径向进给，每次径向进给的距离k＝k+0.1，初始值k＝0；重复上述步骤，计算圆柱齿轮插齿加工每一步筛选的点的面积，直到插齿刀加工完圆柱齿轮为止。上述方案中，步骤三中所述插齿刀的齿顶圆等距离分布点的坐标为：其中fail表示插齿刀齿顶圆的范围角，ra表示插齿刀齿顶圆的半径；由于插齿刀齿顶圆上的点需要进行矩阵M1S变换，公式(8)中的点变换后得到插齿刀齿顶圆上的横纵坐标(x3,y3)为：步骤三中所述插齿刀的左渐开线上等距离分布点的坐标为：由于插齿刀左渐开线上的点需要进行矩阵M1S变换，公式(10)中的点变换后得到插齿刀左渐开线上点的横纵坐标(x4,y4)为：其中，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)。上述方案中，步骤三中所述插齿刀的右渐开线上等距离分布点的坐标为：公式(12)中的点经过矩阵M1S变换后得到插齿刀右渐开线上点的横纵坐标(x5,y5)为：其中，p1表示插齿刀右渐开线上点的极角；将公式(9)、(11)、(13)中的点集也即插齿刀齿顶圆、左右渐开线上分布的均匀点储存在矩阵V2中，并用(x6,y6)表示插齿刀分布的点集的横纵坐标；矩阵V2中的点到工件圆心的距离用p1表示为p1＝(x62+y62)(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.一种圆柱齿轮插齿加工切削力预测的方法，其特征在于，包括以下步骤：步骤一、计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点，根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵，将圆环内的点经过坐标变换矩阵变换为插齿刀坐标系中的点；步骤二、插齿刀每进行一步切削时，将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角大小为判据，筛选出工件在插齿刀内的点，点的数量即为切屑面积的大小a；步骤三、在插齿刀的齿顶圆和左、右渐开线上等距离分布数据点；步骤四、插齿刀每进行一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点，点的数量即为接触长度的大小b；步骤五，基于单位切削力公式：切削方向Ft＝Ktca+Kteb，径向方向Ff＝Kfca+Kfeb，轴向方向Fr＝Krca+Kreb，计算圆柱齿轮插齿加工切削力的大小，其中Ktc、Kfc、Krc为切削系数，Kte、Kfe、Kre为边缘系数。

【技术特征摘要】
1.一种圆柱齿轮插齿加工切削力预测的方法，其特征在于，包括以下步骤：步骤一、计算圆柱齿轮与插齿刀加工时的接触区域圆环的面积，并在圆环中均匀分布足够多的数据点，根据圆柱齿轮与插齿刀加工时的位置关系求解出坐标变换矩阵，将圆环内的点经过坐标变换矩阵变换为插齿刀坐标系中的点；步骤二、插齿刀每进行一步切削时，将坐标变换后的点进行判别，以插齿刀轮齿实体位置的极径与极角大小为判据，筛选出工件在插齿刀内的点，点的数量即为切屑面积的大小a；步骤三、在插齿刀的齿顶圆和左、右渐开线上等距离分布数据点；步骤四、插齿刀每进行一步切削时，将插齿刀与工件接触的点进行判别，以插齿刀齿廓上分布的等距离点到工件中心的距离为判据，筛选出插齿刀与工件接触的点，点的数量即为接触长度的大小b；步骤五，基于单位切削力公式：切削方向Ft＝Ktca+Kteb，径向方向Ff＝Kfca+Kfeb，轴向方向Fr＝Krca+Kreb，计算圆柱齿轮插齿加工切削力的大小，其中Ktc、Kfc、Krc为切削系数，Kte、Kfe、Kre为边缘系数。2.根据权利要求1所述的圆柱齿轮插齿加工切削力预测的方法，其特征在于，步骤一中所述圆柱齿轮与插齿刀加工时的接触区域圆环的表达式为：其中，x1、y1分别表示齿顶圆的横、纵坐标，x2、y2分别表示齿根圆的横、纵坐标；x表示变化角，r1、r2分别表示工件齿顶圆、齿根圆的半径。3.根据权利要求2所述的圆柱齿轮插齿加工切削力预测的方法，其特征在于，在公式(1)、公式(2)接触区域内分布数据点的坐标方程式为：其中，r表示极径，i表示第i条等距圆环，通过公式(3)在每条等距曲线上均匀取j个点得到点集Nij(xij,yij)。4.根据权利要求1所述的圆柱齿轮插齿加工切削力预测的方法，其特征在于，步骤一中圆柱齿轮固定在坐标系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任一一点的坐标。5.根据权利要求4所述的圆柱齿轮插齿加工切削力预测的方法，其特征在于，步骤二中将圆柱齿轮上的等距曲线上的点Nij(xij,yij)通过所述坐标变换矩阵Mn转化成插齿刀坐标系下的点N'ij(xsij,ysij)，并求出N'ij(xsij,ysij)点在插齿刀坐标系下的极角与极径6.根据权利要求5所述的圆柱齿轮插齿加工切削力预测的方法，其特征在于，所述步骤二中的判据为：pij≤r(6)其中，sita表示插齿刀每个齿对应的圆心角，ct2表示齿顶圆上渐开线上展角，r表示为对应点与插齿刀圆心连线与插齿刀齿廓交点的距离。7.根据权利要...

【专利技术属性】
技术研发人员：韩星会，王利波，华林，
申请(专利权)人：武汉理工大学，
类型：发明
国别省市：湖北,42