CN106326575A - 一种基于坐标法的非圆齿轮检测方法 - Google Patents

一种基于坐标法的非圆齿轮检测方法 Download PDF

Info

Publication number
CN106326575A
CN106326575A CN201610741223.2A CN201610741223A CN106326575A CN 106326575 A CN106326575 A CN 106326575A CN 201610741223 A CN201610741223 A CN 201610741223A CN 106326575 A CN106326575 A CN 106326575A
Authority
CN
China
Prior art keywords
tooth
flank
pitch
prime
deviation
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.)
Granted
Application number
CN201610741223.2A
Other languages
English (en)
Other versions
CN106326575B (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.)
Shanxi Pingyang Industry Machinery Co Ltd
Original Assignee
Shanxi Pingyang Industry Machinery Co Ltd
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 Shanxi Pingyang Industry Machinery Co Ltd filed Critical Shanxi Pingyang Industry Machinery Co Ltd
Priority to CN201610741223.2A priority Critical patent/CN106326575B/zh
Publication of CN106326575A publication Critical patent/CN106326575A/zh
Application granted granted Critical
Publication of CN106326575B publication Critical patent/CN106326575B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/18Manufacturability analysis or optimisation for manufacturability

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)
  • Gears, Cams (AREA)

Abstract

本发明具体为一种基于坐标法的非圆齿轮检测方法,解决了现有齿轮测量方法存在无法对非圆齿轮进行检测的问题。一种基于坐标法的非圆齿轮检测方法,采用如下步骤:a、在非圆齿轮的齿廓上进行坐标点采集;b、建立理论齿廓的曲线模型;c、计算处理齿廓总偏差Fα;d、计算处理齿廓形状偏差f;e、计算处理齿廓倾斜偏差f;f、计算处理单个齿距偏差fpt;g、计算处理齿距累积偏差Fpk;h、计算处理齿距累积总偏差Fp;i、计算处理齿厚偏差Esn;j、计算处理径向跳动Fr。本发明大幅提高了检测精度,同时检测效率较高。

Description

一种基于坐标法的非圆齿轮检测方法
技术领域
本发明涉及齿轮检测方法,具体为一种基于坐标法的非圆齿轮检测方法。
背景技术
非圆齿轮机构以变传动比形式传递空间任意两轴间的运动和动力,兼有凸轮机构和圆柱齿轮机构的优点,能够实现变速回转传动,随着非圆齿轮啮合理论研究的深入和制造方法的现代化,非圆齿轮在工业中的应用日益广泛。
同圆柱齿轮一样,非圆齿轮的质量控制中也要求对齿廓偏差、齿距偏差、齿距累积偏差、齿厚偏差、径向跳动等误差项目进行检测,但由于非圆齿轮的特殊性,其测量方法与普通圆柱齿轮和圆锥齿轮的测量方法也有很大的不同。
传统齿轮的测量是在万能齿轮测量机、渐开线检查仪等一些专用的齿轮测量仪器上进行的,这些仪器的工作原理主要基于展开法,要求在整个齿轮测量过程中基圆半径保持固定不变,但由于非圆齿轮的各齿向径不同,对应的基圆半径也不同,因此基于此原理设计的齿轮检测仪器均不能检测非圆齿轮。
如目前世界知名的齿轮检测仪器供应商德国克林贝格(Klingelnberg)、美国格里森(Gleason)、德国温泽(Wenzel)、德国马尔(Mahr)等公司的齿轮检测仪器以及国内较知名的如哈量、成量、爱德华等专业齿轮检测仪器公司生产的齿轮检测仪器均只能检测圆柱齿轮、圆锥齿轮等,而不能检测非圆齿轮。因此,非圆齿轮的精密检测目前在世界上尚处于空白状态。
发明内容
本发明为了解决现有齿轮测量方法存在无法对非圆齿轮进行检测的问题,提供了一种基于坐标法的非圆齿轮检测方法。
本发明是采用如下技术方案实现的:一种基于坐标法的非圆齿轮检测方法,采用如下步骤:a、在非圆齿轮的齿廓上进行坐标点采集;b、根据非圆齿轮节曲线方程R=f(φ),并用公式化为直角坐标,用当量齿轮法建立理论齿廓的曲线模型;c、作齿廓实测点到对应当量齿轮基圆的切线,求出实测点到切点的长度,即实测展长Lai,再求出切线或其沿长线与对应理论齿廓曲线的交点到切点的长,即理论展长Lni,实测展长与理论展长之差,即为该齿廓实测点偏离理论齿廓的偏差值eFαi,每个齿面上所有实测点的eFαi中最大和最小值的代数差即为该齿面的实测齿廓总偏差Fα;d、将计算得到的齿面所有实测点的理论展长Lni作为xi坐标,对应的偏差值eFαi作为yi坐标,按最小二乘法拟合为直线再求出所有偏差值eFαi与对应实测点理论展长Lni在拟合直线上对应的值之差,则所有差值中最大和最小值的代数差即为齿面的实测齿廓形状偏差f;e、将齿面有效齿廓起始点对应的理论展长Lni与有效齿廓终止点对应的理论展长Lni分别作为x坐标值代入d步骤中拟合得到的直线中,分别求出对应的起始点y值和终止点y值,起始点y值和终止点y值的代数差即为该齿面的实测齿廓倾斜偏差f;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
当量齿轮法建立理论齿廓的曲线模型采用如下步骤: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步骤,建立所有齿的理论齿廓曲线。
本发明所述的非圆齿轮检测方法填补了国内外专业齿轮测量仪器均无法检测非圆齿轮的空白,而且大幅提高了检测精度,同时只需采集一次数据,即可得出齿廓总偏差、齿廓形状偏差、齿廓倾斜偏差、齿距偏差、齿距累积偏差、齿距累积总偏差、齿厚偏差、径向跳动等多项偏差值,检测效率较高。
附图说明
图1为本发明中非圆齿轮的结构示意图;
图2为本发明中非圆齿轮节圆均分点;
图3为本发明中建立理论齿廓曲线模型的示意图;
图4为本发明中齿廓偏差的示意图。
具体实施方式
一种基于坐标法的非圆齿轮检测方法,采用如下步骤:
a、在非圆齿轮的齿廓上进行坐标点采集;
b、根据非圆齿轮节曲线方程R=f(φ),并用公式化为直角坐标,用当量齿轮法建立理论齿廓的曲线模型;
c、作齿廓实测点到对应当量齿轮基圆的切线,求出实测点到切点的长度,即实测展长Lai,再求出切线或其沿长线与对应理论齿廓曲线的交点到切点的长,即理论展长Lni,实测展长与理论展长之差,即为该齿廓实测点偏离理论齿廓的偏差值eFαi,每个齿面上所有实测点的eFαi中最大和最小值的代数差即为该齿面的实测齿廓总偏差Fα
d、将计算得到的齿面所有实测点的理论展长Lni作为xi坐标,对应的偏差值eFαi作为yi坐标,按最小二乘法拟合为直线再求出所有偏差值eFαi与对应实测点理论展长Lni在拟合直线上对应的值之差,则所有差值中最大和最小值的代数差即为齿面的实测齿廓形状偏差f
e、将齿面有效齿廓起始点对应的理论展长Lni与有效齿廓终止点对应的理论展长Lni分别作为x坐标值代入d步骤中拟合得到的直线中,分别求出对应的起始点y值和终止点y值,起始点y值和终止点y值的代数差即为该齿面的实测齿廓倾斜偏差f
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
当量齿轮法建立理论齿廓的曲线模型采用如下步骤:
a、计算出非圆齿轮节曲线的周长L;
b、将非圆齿轮节曲线平均分为4×z份,得到4×z个节曲线均分点,每4个均分点为一组,分别记为右齿面点PR、齿中点Pm、左齿面点PL和齿槽中点Pa
c、用下式计算齿中点Pm处节曲线的曲率中心坐标o(a,b)和曲率半径ρ:
a = x - y ′ ( x ′ 2 + y ′ 2 ) x ′ y ′ ′ - x ′ ′ y ′ b = y + x ′ ( x ′ 2 + y ′ 2 ) x ′ y ′ ′ - x ′ ′ y ′
ρ = ( x ′ 2 + y ′ 2 ) 3 2 | x ′ y ′ ′ - x ′ ′ y ′ |
d、将曲率中心坐标o(a,b)作为该齿的当量齿轮基圆圆心,曲率半径ρ作为该齿的当量齿轮节圆半径,并用公式rb=ρcosα计算出基圆半径,然后应用下式分别计算出该齿的左右齿面理论齿廓曲线坐标点:
式中为滚动角,当由0逐渐增大时,得到右齿面理论齿廓;当由0逐渐减小时,得到左齿面理论齿廓;
e、将计算得到的左右齿面理论齿廓分别绕当量齿轮基圆圆心o(a,b)旋转,使左齿面理论齿廓通过非圆齿轮节曲线上该齿的左齿面点PL;右齿面理论齿廓通过非圆齿轮节曲线上该齿的右齿面点PR
f、将旋转后的左右齿面理论齿廓线超出齿顶曲线和齿根曲线的部分舍去,剩余部分即为该齿的左右齿面理论齿廓曲线;
j、重复a-f步骤,建立所有齿的理论齿廓曲线。
具体实施过程中,非圆齿轮模数m=1、齿数z=62、齿轮压力角α=20°,其节曲线方程为其中齿轮偏心率k=0.2,节曲线长轴半径a=30.40926,直线中k是指斜率。
一种基于坐标法的非圆齿轮检测方法,采用如下步骤:
a、用坐标测量机以连续扫描采点方式在被测齿轮的齿廓上进行坐标点采集,得到所有齿的实测齿廓坐标点,以下列出61#齿左齿面的实测齿廓坐标点:
b、用公式将非圆齿轮方程化为直角坐标,用当量齿轮法建立理论齿廓的曲线模型。
c、求出第1个实测点(33.9645,6.7188)到61#齿对应的当量齿轮基圆的切线(位于渐开线内侧的切线)的切点坐标,求出实测点到切点的长度(即实测展长)Lai=0.5756;求出切线(或其沿长线)与61#齿左齿面理论齿廓曲线的交点到切点的长(即理论展长)Lni=0.5796;求出实测展长与理论展长之差(即该实测点偏离理论齿廓的偏差值)eFαi=Lai-Lni=-0.0040;求出该齿面所有点的e列表如下:
则该齿面的实测齿廓总偏差Fα=max(eFαi)-min(eFαi)=0.0340;
d、将c步骤中计算得到的61#齿左齿面所有实测点的理论展长Lni作为xi坐标,对应的偏差值eFαi作为yi坐标,按最小二乘法拟合为直线再求出所有偏差值eFαi与对应实测点理论展长在拟合直线上对应的值之差列表如下:
则61#齿左齿面齿廓形状偏差f=max(eFαi-yi)-min(eFαi-yi)=0.0290;
e、将61#齿左齿面第1点对应的理论展长Ln1=0.5796与最后一点(第20点)对应的理论展长Ln20=8.9663,分别作为x坐标值代入d步骤中拟合得到的直线中,求出对应的y1=-0.0092,y20=0.0101,则这两个y值之代数差即为该齿面的实测齿廓倾斜偏差f=y20-y1=0.0193;
f、求出所有实测齿廓曲线与节曲线的交点(如第61#齿左齿面实测齿廓曲线与节曲线的交点为35.1575,6.8961),并计算相邻同侧齿面上两交点间的节曲线弧长即实测齿距(如61#齿与60#齿左齿面实测齿距为3.1369),则其与公称齿距3.1416的代数差即为该齿的该齿面单齿齿距偏差fpti(如61#齿与60#齿左齿面的单齿齿距偏差fpti=-0.0047),所有齿的该齿面单齿齿距偏差中绝对值最大的偏差值即为齿轮该齿面的单个齿距偏差fpt=max(|fpti|)=0.0251;
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;
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;
i、求出所有实测齿廓曲线与节曲线的交点(f步骤已求出,可直接使用),并计算每个齿左、右齿面上该交点间的节曲线弧长即实测齿厚(如61#齿实测齿厚为3.1432),则其与公称齿厚3.1416的代数差即为该齿的单齿齿厚偏差Esn(如61#齿的单齿齿厚偏差为Esn=0.0016),所有偏差中的最大值为实测齿厚上偏差Esns=0.0195,最小值为实测齿厚下偏差Esni=-0.0394;
j、求出所有实测齿廓曲线与节曲线的交点(f步骤已求出,可直接使用),并计算每个齿槽左、右齿面上该交点间的节曲线弧长即实测齿槽宽(如61#齿实测齿槽宽为3.1310),再计算出实测齿槽宽与公称齿槽宽3.1416的差值即齿槽宽偏差e(如61#齿实测齿槽宽偏差e=-0.0106),则即为该齿槽相对理论正确位置的单齿径向跳动偏差Fri(如61#齿单齿径向跳动偏差Fr61=0.0146),所有齿的单齿径向跳动偏差中最大和最小值的代数差即为该齿轮的径向跳动Fr=max(Fri)-min(Fri)=0.0493。
当量齿轮法建立理论齿廓的曲线模型采用如下步骤:
a、计算出被测齿轮的节曲线周长L=πmz=194.7787;
b、将非圆齿轮节曲线平均分为4×62=248份,得到248个节曲线均分点,部分均分点列表如下:
c、以61#齿为例,计算齿中点Pm处的节曲线的曲率中心坐标o(a,b)和曲率半径ρ:
a = x - y ′ ( x ′ 2 + y ′ 2 ) x ′ y ′ ′ - x ′ ′ y ′ = 17.0884 b = y + x ′ ( x ′ 2 + y ′ 2 ) x ′ y ′ ′ - x ′ ′ y ′ = - 0.2657
ρ = ( x ′ 2 + y ′ 2 ) 3 2 | x ′ y ′ ′ - x ′ ′ y ′ | = 19.4366
d、将c步骤中计算得到的曲率中心o(17.0884,-0.2657)作为该齿的当量齿轮基圆圆心,曲率半径ρ=19.4366作为该齿的当量齿轮节圆半径,并计算出基圆半径rb=ρcosα=18.2644,然后应用公式分别计算出61#齿的左齿面理论齿廓(渐开线)曲线坐标点列表如下:
点位号 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
e、将计算得到的左齿面理论齿廓分别绕当量齿轮基圆圆心o旋转,使左齿面理论齿廓通过非圆齿轮节曲线上61#齿的左齿面点PL(35.1565,6.8986);
f、将旋转后的左齿面理论齿廓线超出齿顶曲线和齿根曲线的部分舍去,剩余部分即为61#齿的左齿面理论齿廓曲线,坐标点列表如下:
点位号 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
j、重复a-f步骤,建立所有齿的理论齿廓曲线。

Claims (2)

1.一种基于坐标法的非圆齿轮检测方法,其特征在于:采用如下步骤:
a、在非圆齿轮的齿廓上进行坐标点采集;
b、根据非圆齿轮节曲线方程R=f(φ),并用公式化为直角坐标,用当量齿轮法建立理论齿廓的曲线模型;
c、作齿廓实测点到对应当量齿轮基圆的切线,求出实测点到切点的长度,即实测展长Lai,再求出切线或其沿长线与对应理论齿廓曲线的交点到切点的长,即理论展长Lni,实测展长与理论展长之差,即为该齿廓实测点偏离理论齿廓的偏差值eFαi,每个齿面上所有实测点的eFαi中最大和最小值的代数差即为该齿面的实测齿廓总偏差Fα
d、将计算得到的齿面所有实测点的理论展长Lni作为xi坐标,对应的偏差值eFαi作为yi坐标,按最小二乘法拟合为直线再求出所有偏差值eFαi与对应实测点理论展长Lni在拟合直线上对应的值之差,则所有差值中最大和最小值的代数差即为齿面的实测齿廓形状偏差f
e、将齿面有效齿廓起始点对应的理论展长Lni与有效齿廓终止点对应的理论展长Lni分别作为x坐标值代入d步骤中拟合得到的直线中,分别求出对应的起始点y值和终止点y值,起始点y值和终止点y值的代数差即为该齿面的实测齿廓倾斜偏差f
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、求出所有实测齿廓曲线与节曲线的交点,并计算每个齿槽左、右齿面上该交点间的节曲线弧长即实测齿槽宽,再计算出实测齿槽宽与公称齿槽宽的差值即齿槽宽偏差e,即为该齿槽相对理论正确位置的单齿径向跳动偏差Fri,所有齿的单齿径向跳动偏差Fri中最大和最小值的代数差即为被测齿轮的径向跳动Fr
2.根据权利要求1所述的一种基于坐标法的非圆齿轮检测方法,其特征在于:当量齿轮法建立理论齿廓的曲线模型采用如下步骤:
a、计算出非圆齿轮节曲线的周长L;
b、将非圆齿轮节曲线平均分为4×z份,得到4×z个节曲线均分点,每4个均分点为一组,分别记为右齿面点PR、齿中点Pm、左齿面点PL和齿槽中点Pa
c、用下式计算齿中点Pm处节曲线的曲率中心坐标o(a,b)和曲率半径ρ:
a = x - y ′ ( x ′ 2 + y ′ 2 ) x ′ y ′ ′ - x ′ ′ y ′ b = y + x ′ ( x ′ 2 + y ′ 2 ) x ′ y ′ ′ - x ′ ′ y ′
ρ = ( x ′ 2 + y ′ 2 ) 3 2 | x ′ y ′ ′ - x ′ ′ y ′ |
d、将曲率中心坐标o(a,b)作为该齿的当量齿轮基圆圆心,曲率半径ρ作为该齿的当量齿轮节圆半径,并用公式rb=ρcosα计算出基圆半径,然后应用下式分别计算出该齿的左右齿面理论齿廓曲线坐标点:
式中为滚动角,当由0逐渐增大时,得到右齿面理论齿廓;当由0逐渐减小时,得到左齿面理论齿廓;
e、将计算得到的左右齿面理论齿廓分别绕当量齿轮基圆圆心o(a,b)旋转,使左齿面理论齿廓通过非圆齿轮节曲线上该齿的左齿面点PL;右齿面理论齿廓通过非圆齿轮节曲线上该齿的右齿面点PR
f、将旋转后的左右齿面理论齿廓线超出齿顶曲线和齿根曲线的部分舍去,剩余部分即为该齿的左右齿面理论齿廓曲线;
j、重复a-f步骤,建立所有齿的理论齿廓曲线。
CN201610741223.2A 2016-08-26 2016-08-26 一种基于坐标法的非圆齿轮检测方法 Active CN106326575B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610741223.2A CN106326575B (zh) 2016-08-26 2016-08-26 一种基于坐标法的非圆齿轮检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610741223.2A CN106326575B (zh) 2016-08-26 2016-08-26 一种基于坐标法的非圆齿轮检测方法

Publications (2)

Publication Number Publication Date
CN106326575A true CN106326575A (zh) 2017-01-11
CN106326575B CN106326575B (zh) 2019-08-02

Family

ID=57791083

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610741223.2A Active CN106326575B (zh) 2016-08-26 2016-08-26 一种基于坐标法的非圆齿轮检测方法

Country Status (1)

Country Link
CN (1) CN106326575B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108627133A (zh) * 2018-05-07 2018-10-09 中国航发哈尔滨轴承有限公司 一种外圈滚道三瓣波形滚子轴承的径向游隙测量方法
CN114800048A (zh) * 2022-05-07 2022-07-29 西安理工大学 一种齿轮在机测量过程中基于齿距测量的径跳检测方法
CN115221655A (zh) * 2022-09-13 2022-10-21 南通市中吕齿轮有限公司 一种斜齿非圆齿轮的传动精度评价方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101419066A (zh) * 2008-06-24 2009-04-29 中国计量学院 一种两轴随动控制检测齿轮齿廓和螺旋线的方法
CN102567596A (zh) * 2012-03-05 2012-07-11 西北工业大学 一种面齿轮轮齿参数化建模方法
CN103712567A (zh) * 2013-05-27 2014-04-09 北京理工大学 小模数圆弧齿轮非接触式光学精密测量方法
CN105889456A (zh) * 2016-05-10 2016-08-24 武汉理工大学 曲线齿非圆齿轮的设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101419066A (zh) * 2008-06-24 2009-04-29 中国计量学院 一种两轴随动控制检测齿轮齿廓和螺旋线的方法
CN102567596A (zh) * 2012-03-05 2012-07-11 西北工业大学 一种面齿轮轮齿参数化建模方法
CN103712567A (zh) * 2013-05-27 2014-04-09 北京理工大学 小模数圆弧齿轮非接触式光学精密测量方法
CN105889456A (zh) * 2016-05-10 2016-08-24 武汉理工大学 曲线齿非圆齿轮的设计方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHAO LIN等: "Pitch deviation measurement and analysis of curve-face gear pair", 《ELSEVIER》 *
MING-FENG TSAY等: "Study on the Generalized Mathematical Model of Noncircular Gears", 《MATHEMATICAL AND COMPUTER MODELLING》 *
曹喜军等: "端曲面齿轮的误差检测及分析研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 *
李静: "非圆齿轮渐开线齿廓数学模型的建立", 《辽宁石油化工大学学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108627133A (zh) * 2018-05-07 2018-10-09 中国航发哈尔滨轴承有限公司 一种外圈滚道三瓣波形滚子轴承的径向游隙测量方法
CN114800048A (zh) * 2022-05-07 2022-07-29 西安理工大学 一种齿轮在机测量过程中基于齿距测量的径跳检测方法
CN114800048B (zh) * 2022-05-07 2023-07-18 西安理工大学 一种齿轮在机测量过程中基于齿距测量的径跳检测方法
CN115221655A (zh) * 2022-09-13 2022-10-21 南通市中吕齿轮有限公司 一种斜齿非圆齿轮的传动精度评价方法

Also Published As

Publication number Publication date
CN106326575B (zh) 2019-08-02

Similar Documents

Publication Publication Date Title
CN102636097B (zh) 基于双面啮合的齿轮齿廓偏差测量方法
RU2015119233A (ru) Способ для определения местоположения эвольвент в зубчатых колесах
CN110297462B (zh) 一种考虑机床几何误差影响的磨齿精度预测建模方法
CN106326575A (zh) 一种基于坐标法的非圆齿轮检测方法
CN103344210A (zh) 一种齿轮误差多自由度评定方法
CN108645323A (zh) 安装误差影响下的渐开线圆柱直齿轮齿廓偏差评定方法
CN108645322B (zh) 安装误差影响下的渐开线圆柱直齿轮齿距偏差评定方法
CN102645200B (zh) 基于双面啮合的齿轮齿距偏差测量方法
CN106556366A (zh) 一种解决微小齿轮测量的柱形测头及测量方法
CN107860313A (zh) 一种基于线结构光的斜齿轮齿向偏差测量方法
CN103791870A (zh) 平面二次包络环面蜗杆测量方法
CN105134910B (zh) 同向渐开线齿轮副啮合传动装置
CN107617794B (zh) 一种人字齿内齿圈加工方法
CN106980724A (zh) 一种基于复杂齿面网格规划原则的齿面精确模型截面放样构建方法
CN103791869A (zh) 一种平面二次包络环面滚刀的测量方法
CN105127521B (zh) 一种用于加工双渐开线齿轮的滚刀及加工方法
CN111666643A (zh) 一种用于确定复杂齿面接触性能的方法
CN105241415B (zh) 一种接触式渐开线蜗杆齿形的测量方法
CN105277165A (zh) 一种齿面等间隔采样提高齿轮测量精度的方法
CN115221655B (zh) 一种斜齿非圆齿轮的传动精度评价方法
CN102797829B (zh) 准双导程锥蜗杆蜗轮设计方法
CN104330065A (zh) 一种螺旋锥齿轮齿面曲率干涉检验方法
CN104864800A (zh) 滚刀前角检测仪及使用方法
CN108562256B (zh) 一种齿形中部不完整圆弧端齿齿顶高测量方法
CN105179600B (zh) 双大负变位渐开线齿轮传动装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant