技术领域
[0001] 本
发明涉及齿轮检测方法，具体为一种基于坐标法的
非圆齿轮检测方法。
背景技术
[0002] 非圆齿轮机构以变
传动比形式传递空间任意两轴间的运动和动
力，兼有
凸轮机构和圆柱齿轮机构的优点，能够实现变速回转传动，随着非圆齿轮
啮合理论研究的深入和制造方法的现代化，非圆齿轮在工业中的应用日益广泛。
[0003] 同圆柱齿轮一样，非圆齿轮的
质量控制中也要求对齿廓偏差、齿距偏差、齿距累积偏差、齿厚偏差、径向跳动等误差项目进行检测，但由于非圆齿轮的特殊性，其测量方法与普通圆柱齿轮和圆
锥齿轮的测量方法也有很大的不同。
[0004] 传统齿轮的测量是在万能齿轮测量机、渐开线检查仪等一些专用的齿轮测量仪器上进行的，这些仪器的工作原理主要基于展开法，要求在整个齿轮测量过程中
基圆半径保持固定不变，但由于非圆齿轮的各齿向径不同，对应的基圆半径也不同，因此基于此原理设计的齿轮检测仪器均不能检测非圆齿轮。
[0005] 如目前世界知名的齿轮检测仪器供应商德国克林贝格(Klingelnberg)、美国格里森(Gleason)、德国温泽(Wenzel)、德国
马尔(Mahr)等公司的齿轮检测仪器以及国内较知名的如哈量、成量、爱德华等专业齿轮检测仪器公司生产的齿轮检测仪器均只能检测圆柱齿轮、圆锥齿轮等，而不能检测非圆齿轮。因此，非圆齿轮的精密检测目前在世界上尚处于空白状态。
发明内容
[0006] 本发明为了解决现有齿轮测量方法存在无法对非圆齿轮进行检测的问题，提供了一种基于坐标法的非圆齿轮检测方法。
[0007] 本发明是采用如下技术方案实现的：一种基于坐标法的非圆齿轮检测方法，采用如下步骤：a、在非圆齿轮的齿廓上进行坐标点采集；b、根据非圆齿轮节曲线方程R＝f(φ)，并用 公式化为直
角坐标，用当量齿轮法建立理论齿廓的曲线模型；c、作齿廓实测点到对应当量齿轮基圆的切线，求出实测点到切点的长度，即实测展长Lai，再求出切线或其沿长线与对应理论齿廓曲线的交点到切点的长，即理论展长Lni，实测展长与理论展长之差，即为该齿廓实测点偏离理论齿廓的偏差值eFαi，每个
齿面上所有实测点的eFαi中最大和最小值的代数差即为该齿面的实测齿廓总偏差Fα；d、将计算得到的齿面所有实测点的理论展长Lni作为xi坐标，对应的偏差值eFαi作为yi坐标，按最小二乘法拟合为直线再求出所有偏差值eFαi与对应实测点理论展长Lni在拟合直线上对应的 值之差，则所有差值中最大和最小值的代数差即为齿面的实测齿廓形状偏差ffα；e、将齿面有效齿廓起始点对应的理论展长Lni与有效齿廓终止点对应的理论展长Lni分别作为x坐标值代入d步骤中拟合得到的直线中，分别求出对应的起始点y值和终止点y值，起始点y值和终止点y值的代数差即为该齿面的实测齿廓倾斜偏差fHα；f、求出所有实测齿廓曲线与节曲线的交点，并计算相邻同侧齿面上两交点间的节曲线弧长即实测齿距，实测齿距与公称齿距的代数差即为该齿的该齿面齿距偏差fpti，所有齿的该齿面单齿齿距偏差fpti中绝对值最大的偏差值即为被测齿轮该齿面的单个齿距偏差fpt；g、求出所有实测齿廓曲线与节曲线的交点，并计算某齿与其后第k(k为自然数)齿的同侧齿面上两交点间的节曲线弧长即为该齿面k齿实测齿距，k齿实测齿距与k齿公称齿距的代数差即为该齿的该齿面k齿齿距偏差，所有齿的该齿面k齿齿距偏差中绝对值最大的偏差值即为被测齿轮该齿面的k齿齿距累积偏差Fpk；h、求出所有实测齿廓曲线与节曲线的交点，并计算每个齿与第1齿的同侧齿面上两交点间的节曲线弧长即该齿相对第1齿的实测齿距，再计算该实测齿距与公称齿距的代数差即每个齿相对第1齿的齿距偏差，所有同侧齿面相对第1齿的齿距偏差中最大和最小值的代数差即为即为被测齿轮该齿面的齿距累积总偏差Fp；i、求出所有实测齿廓曲线与节曲线的交点，并计算每个齿左、右齿面上该交点间的节曲线弧长即实测齿厚，则实测齿厚与公称齿厚的代数差即为该齿的单齿齿厚偏差Esn，所有偏差中的最大值为被测齿轮的实测齿厚上偏差Esns，最小值为被测齿轮的实测齿厚下偏差Esni；j、求出所有实测齿廓曲线与节曲线的交点，并计算每个
齿槽左、右齿面上该交点间的节曲线弧长即实测齿槽宽，再计算出实测齿槽宽与公称齿槽宽的差值即齿槽宽偏差 即为该齿槽相对理论正确
位置的单齿径向跳动偏差Fri，所有齿的单齿径向跳动偏差Fri中最大和最小值的代数差即为被测齿轮的径向跳动Fr。
[0008] 当量齿轮法建立理论齿廓的曲线模型采用如下步骤：a、计算出非圆齿轮节曲线的周长L；b、将非圆齿轮节曲线平均分为4×z份，得到4×z个节曲线均分点，每4个均分点为一组，分别记为右齿面点PR、齿中点Pm、左齿面点PL和齿槽中点Pa；c、用下式计算齿中点Pm处节曲线的
曲率中心坐标o(a,b)和
曲率半径ρ：d、将曲率中心坐标o(a,b)作为该齿的当量齿轮基圆圆心，曲率半径ρ作为该齿的当量齿轮
节圆半径，并用公式rb＝ρcosα计算出基圆半径，然后应用下式分别计算出该齿的左右齿面理论齿廓曲线坐标点： 式中 为滚动角，当 由0逐
渐增大时，得到右齿面理论齿廓；当 由0逐渐减小时，得到左齿面理论齿廓；e、将计算得到的左右齿面理论齿廓分别绕当量齿轮基圆圆心o(a,b)旋转，使左齿面理论齿廓通过非圆齿轮节曲线上该齿的左齿面点PL；右齿面理论齿廓通过非圆齿轮节曲线上该齿的右齿面点PR；
f、将旋转后的左右齿面理论齿廓线超出齿顶曲线和齿根曲线的部分舍去，剩余部分即为该齿的左右齿面理论齿廓曲线；j、重复a-f步骤，建立所有齿的理论齿廓曲线。
[0009] 本发明所述的非圆齿轮检测方法填补了国内外专业齿轮测量仪器均无法检测非圆齿轮的空白，而且大幅提高了检测
精度，同时只需采集一次数据，即可得出齿廓总偏差、齿廓形状偏差、齿廓倾斜偏差、齿距偏差、齿距累积偏差、齿距累积总偏差、齿厚偏差、径向跳动等多项偏差值，检测效率较高。
附图说明
[0010] 图1为本发明中非圆齿轮的结构示意图；
[0011] 图2为本发明中非圆齿轮节圆均分点；
[0012] 图3为本发明中建立理论齿廓曲线模型的示意图；
[0013] 图4为本发明中齿廓偏差的示意图。
具体实施方式
[0014] 一种基于坐标法的非圆齿轮检测方法，采用如下步骤：
[0015] a、在非圆齿轮的齿廓上进行坐标点采集；
[0016] b、根据非圆齿轮节曲线方程R＝f(φ)，并用 公式化为直角坐标，用当量齿轮法建立理论齿廓的曲线模型；
[0017] c、作齿廓实测点到对应当量齿轮基圆的切线，求出实测点到切点的长度，即实测展长Lai，再求出切线或其沿长线与对应理论齿廓曲线的交点到切点的长，即理论展长Lni，实测展长与理论展长之差，即为该齿廓实测点偏离理论齿廓的偏差值eFαi，每个齿面上所有实测点的eFαi中最大和最小值的代数差即为该齿面的实测齿廓总偏差Fα；
[0018] d、将计算得到的齿面所有实测点的理论展长Lni作为xi坐标，对应的偏差值eFαi作为yi坐标，按最小二乘法拟合为直线 再求出所有偏差值eFαi与对应实测点理论展长Lni在拟合直线上对应的 值之差，则所有差值中最大和最小值的代数差即为齿面的实测齿廓形状偏差ffα；
[0019] e、将齿面有效齿廓起始点对应的理论展长Lni与有效齿廓终止点对应的理论展长Lni分别作为x坐标值代入d步骤中拟合得到的直线中，分别求出对应的起始点y值和终止点y值，起始点y值和终止点y值的代数差即为该齿面的实测齿廓倾斜偏差fHα；
[0020] f、求出所有实测齿廓曲线与节曲线的交点，并计算相邻同侧齿面上两交点间的节曲线弧长即实测齿距，实测齿距与公称齿距的代数差即为该齿的该齿面单齿齿距偏差fpti，所有齿的该齿面单齿齿距偏差fpti中绝对值最大的偏差值即为被测齿轮该齿面的单个齿距偏差fpt；
[0021] g、求出所有实测齿廓曲线与节曲线的交点，并计算某齿与其后第k(k为自然数)齿的同侧齿面上两交点间的节曲线弧长即为该齿面k齿实测齿距，k齿实测齿距与k齿公称齿距的代数差即为该齿的该齿面k齿齿距偏差，所有齿的该齿面k齿齿距偏差中绝对值最大的偏差值即为被测齿轮该齿面的k齿齿距累积偏差Fpk；
[0022] h、求出所有实测齿廓曲线与节曲线的交点，并计算每个齿与第1齿的同侧齿面上两交点间的节曲线弧长即该齿相对第1齿的实测齿距，再计算该实测齿距与公称齿距的代数差即每个齿相对第1齿的齿距偏差，所有同侧齿面相对第1齿的齿距偏差中最大和最小值的代数差即为即为被测齿轮该齿面的齿距累积总偏差Fp；
[0023] i、求出所有实测齿廓曲线与节曲线的交点，并计算每个齿左、右齿面上该交点间的节曲线弧长即实测齿厚，则实测齿厚与公称齿厚的代数差即为该齿的单齿齿厚偏差Esn，所有偏差中的最大值为被测齿轮的实测齿厚上偏差Esns，最小值为被测齿轮的实测齿厚下偏差Esni；
[0024] j、求出所有实测齿廓曲线与节曲线的交点，并计算每个齿槽左、右齿面上该交点间的节曲线弧长即实测齿槽宽，再计算出实测齿槽宽与公称齿槽宽的差值即齿槽宽偏差即为该齿槽相对理论正确位置的单齿径向跳动偏差Fri，所有齿的单齿径向跳动偏差Fri中最大和最小值的代数差即为被测齿轮的径向跳动Fr。
[0025] 当量齿轮法建立理论齿廓的曲线模型采用如下步骤：
[0026] a、计算出非圆齿轮节曲线的周长L；
[0027] b、将非圆齿轮节曲线平均分为4×z份，得到4×z个节曲线均分点，每4个均分点为一组，分别记为右齿面点PR、齿中点Pm、左齿面点PL和齿槽中点Pa；
[0028] c、用下式计算齿中点Pm处节曲线的曲率中心坐标o(a,b)和曲率半径ρ：
[0029]
[0030]
[0031] d、将曲率中心坐标o(a,b)作为该齿的当量齿轮基圆圆心，曲率半径ρ作为该齿的当量齿轮节圆半径，并用公式rb＝ρcosα计算出基圆半径，然后应用下式分别计算出该齿的左右齿面理论齿廓曲线坐标点：
[0032]
[0033] 式中 为滚动角，当 由0逐渐增大时，得到右齿面理论齿廓；当 由0逐渐减小时，得到左齿面理论齿廓；
[0034] e、将计算得到的左右齿面理论齿廓分别绕当量齿轮基圆圆心o(a,b)旋转，使左齿面理论齿廓通过非圆齿轮节曲线上该齿的左齿面点PL；右齿面理论齿廓通过非圆齿轮节曲线上该齿的右齿面点PR；
[0035] f、将旋转后的左右齿面理论齿廓线超出齿顶曲线和齿根曲线的部分舍去，剩余部分即为该齿的左右齿面理论齿廓曲线；
[0036] j、重复a-f步骤，建立所有齿的理论齿廓曲线。
[0037] 具体实施过程中，非圆齿轮模数m＝1、齿数z＝62、齿轮
压力角α＝20°，其节曲线方程为 其中齿轮偏心率k＝0.2，节曲线长轴半径a＝30.40926，直线中k是指斜率。
[0038] 一种基于坐标法的非圆齿轮检测方法，采用如下步骤：
[0039] a、用坐标测量机以连续扫描采点方式在被测齿轮的齿廓上进行坐标点采集，得到所有齿的实测齿廓坐标点，以下列出61#齿左齿面的实测齿廓坐标点：
[0040]
[0041]
[0042] b、用 公式将非圆齿轮方程化为直角坐标，用当量齿轮法建立理论齿廓的曲线模型。
[0043] c、求出第1个实测点(33.9645,6.7188)到61#齿对应的当量齿轮基圆的切线(位于渐开线内侧的切线)的切点坐标，求出实测点到切点的长度(即实测展长)Lai＝0.5756；求出切线(或其沿长线)与61#齿左齿面理论齿廓曲线的交点到切点的长(即理论展长)Lni＝0.5796；求出实测展长与理论展长之差(即该实测点偏离理论齿廓的偏差值)eFαi＝Lai-Lni＝-0.0040；求出该齿面所有点的eFα列表如下：
[0044]
[0045]
[0046] 则该齿面的实测齿廓总偏差Fα＝max(eFαi)-min(eFαi)＝0.0340；
[0047] d、将c步骤中计算得到的61#齿左齿面所有实测点的理论展长Lni作为xi坐标，对应的偏差值eFαi作为yi坐标，按最小二乘法拟合为直线 再求出所有偏差值eFαi与对应实测点理论展长在拟合直线上对应的 值之差 列表如下：
[0048]
[0049] 则61#齿左齿面齿廓形状偏差ffα＝max(eFαi-yi)-min(eFαi-yi)＝0.0290；
[0050] e、将61#齿左齿面第1点对应的理论展长Ln1＝0.5796与最后一点(第20点)对应的理论展长Ln20＝8.9663，分别作为x坐标值代入d步骤中拟合得到的直线中，求出对应的y1＝-0.0092，y20＝0.0101，则这两个y值之代数差
即为该齿面的实测齿廓倾斜偏差fHα＝y20-y1＝0.0193；
[0051] f、求出所有实测齿廓曲线与节曲线的交点(如第61#齿左齿面实测齿廓曲线与节曲线的交点为35.1575，6.8961)，并计算相邻同侧齿面上两交点间的节曲线弧长即实测齿距(如61#齿与60#齿左齿面实测齿距为3.1369)，则其与公称齿距3.1416的代数差即为该齿的该齿面单齿齿距偏差fpti(如61#齿与60#齿左齿面的单齿齿距偏差fpti＝-0.0047)，所有齿的该齿面单齿齿距偏差中绝对值最大的偏差值即为齿轮该齿面的单个齿距偏差fpt＝max(|fpti|)＝0.0251；
[0052] g、求出所有实测齿廓曲线与节曲线的交点(f步骤已求出，可直接使用)，并计算每个齿与其后第4个齿的同侧齿面上两交点间的节曲线弧长即该齿面4齿实测齿距(如61#齿与2#齿左齿面实测齿距为12.5589)，则其与4齿公称齿距12.5664的代数差即为即为该齿的该齿面4齿齿距偏差Fp4i(如61#齿的2#齿齿距偏差Fp4i＝-0.0075)，所有齿的该齿面4齿齿距偏差中绝对值最大的偏差值即为该齿面4齿齿距累积偏差Fp4＝max(|Fp4i|)＝0.0188；
[0053] h、求出所有实测齿廓曲线与节曲线的交点(f步骤已求出，可直接使用)，并计算每个齿与1#齿的同侧齿面上两交点间的节曲线弧长即该齿相对1#齿的实测齿距(如61#齿与1#齿左齿面实测齿距为6.2953)，再计算该实测齿距与公称齿距(如61#齿与1#齿左齿面的公称齿距为6.2832)的代数差即每个齿相对1#齿的齿距偏差Fpi(如61#齿相对1#齿的齿距偏差为Fpi＝0.0121)，所有同侧齿面相对1#齿的齿距偏差中最大和最小值的代数差即为该齿面齿距累积总偏差Fp＝max(Fpi)-min(Fpi)＝0.0492；
[0054] i、求出所有实测齿廓曲线与节曲线的交点(f步骤已求出，可直接使用)，并计算每个齿左、右齿面上该交点间的节曲线弧长即实测齿厚(如61#齿实测齿厚为3.1432)，则其与公称齿厚3.1416的代数差即为该齿的单齿齿厚偏差Esn(如61#齿的单齿齿厚偏差为Esn＝0.0016)，所有偏差中的最大值为实测齿厚上偏差Esns＝0.0195，最小值为实测齿厚下偏差Esni＝-0.0394；
[0055] j、求出所有实测齿廓曲线与节曲线的交点(f步骤已求出，可直接使用)，并计算每个齿槽左、右齿面上该交点间的节曲线弧长即实测齿槽宽(如61#齿实测齿槽宽为3.1310)，再计算出实测齿槽宽与公称齿槽宽3.1416的差值即齿槽宽偏差e(如61#齿实测齿槽宽偏差e＝-0.0106)，则 即为该齿槽相对理论正确位置的单齿径向跳动偏差Fri(如61#齿单齿径向跳动偏差Fr61＝0.0146)，所有齿的单齿径向跳动偏差中最大和最小值的代数差即为该齿轮的径向跳动Fr＝max(Fri)-min(Fri)＝0.0493。
[0056] 当量齿轮法建立理论齿廓的曲线模型采用如下步骤：
[0057] a、计算出被测齿轮的节曲线周长L＝πmz＝194.7787；
[0058] b、将非圆齿轮节曲线平均分为4×62＝248份，得到248个节曲线均分点，部分均分点列表如下：
[0059]
[0060] c、以61#齿为例，计算齿中点Pm处的节曲线的曲率中心坐标o(a,b)和曲率半径ρ：
[0061]
[0062]
[0063] d、将c步骤中计算得到的曲率中心o(17.0884,-0.2657)作为该齿的当量齿轮基圆圆心，曲率半径ρ＝19.4366作为该齿的当量齿轮节圆半径，并计算出基圆半径rb＝ρcosα＝18.2644，然后应用 公式分别计算出61#齿的左齿面理
论齿廓(渐开线)曲线坐标点列表如下：
[0064]点位号 X坐标 Y坐标 点位号 X坐标 Y坐标
1 35.3528 -0.2657 11 36.4432 -0.5257
2 35.4550 -0.2729 12 36.5526 -0.5668
3 35.5650 -0.2874 13 36.6619 -0.6102
4 35.6749 -0.3064 14 36.7711 -0.6555
5 35.7849 -0.3292 15 36.8801 -0.7030
6 35.8948 -0.3552 16 36.9890 -0.7524
7 36.0046 -0.3842 17 37.0977 -0.8037
8 36.1144 -0.4158 18 37.2063 -0.8569
9 36.2241 -0.4501 19 37.3146 -0.9119
10 36.3337 -0.4867 20 37.4228 -0.9688
[0065] e、将计算得到的左齿面理论齿廓分别绕当量齿轮基圆圆心o旋转，使左齿面理论齿廓通过非圆齿轮节曲线上61#齿的左齿面点PL(35.1565,6.8986)；
[0066] f、将旋转后的左齿面理论齿廓线超出齿顶曲线和齿根曲线的部分舍去，剩余部分即为61#齿的左齿面理论齿廓曲线，坐标点列表如下：
[0067]点位号 X坐标 Y坐标 点位号 X坐标 Y坐标
1 33.9645 6.7188 11 35.0715 6.8956
2 34.0617 6.7512 12 35.1883 6.8994
3 34.1689 6.7799 13 35.3059 6.9012
4 34.2778 6.8044 14 35.4241 6.9010
5 34.3880 6.8254 15 35.5430 6.8989
6 34.4995 6.8434 16 35.6625 6.8949
7 34.6121 6.8586 17 35.7826 6.8891
8 34.7256 6.8714 18 35.9032 6.8814
9 34.8401 6.8817 19 36.0244 6.8720
10 34.9554 6.8898 20 36.1461 6.8608
[0068] j、重复a-f步骤，建立所有齿的理论齿廓曲线。
