CN112257192B - 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 - Google Patents

一种基于实测齿廓偏差的一齿切向综合偏差计算方法 Download PDF

Info

Publication number
CN112257192B
CN112257192B CN202010973033.XA CN202010973033A CN112257192B CN 112257192 B CN112257192 B CN 112257192B CN 202010973033 A CN202010973033 A CN 202010973033A CN 112257192 B CN112257192 B CN 112257192B
Authority
CN
China
Prior art keywords
tooth profile
tooth
gear
deviation
tangential
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
CN202010973033.XA
Other languages
English (en)
Other versions
CN112257192A (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.)
Beijing University of Technology
Original Assignee
Beijing University of 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 Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202010973033.XA priority Critical patent/CN112257192B/zh
Publication of CN112257192A publication Critical patent/CN112257192A/zh
Application granted granted Critical
Publication of CN112257192B publication Critical patent/CN112257192B/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

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)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Length Measuring Devices With Unspecified Measuring Means (AREA)

Abstract

本发明公开了一种基于实际测量齿廓偏差曲线计算齿轮一齿切向综合偏差的方法,包括理想齿廓建模、实际齿廓偏差拟合、将齿廓偏差曲线与齿廓进行旋转与啮合、进行一齿切向综合偏差计算、如何划分旋转角度可以在保障计算精度的前提下,有效缩减计算时间等内容。该方法旨在仅通过齿轮测量仪器对实际齿廓进行测量,得到实际齿廓信息,进而计算该齿轮的一齿切向综合偏差。

Description

一种基于实测齿廓偏差的一齿切向综合偏差计算方法
技术领域
本发明属于齿轮精密测量技术领域,具体为一种基于实测齿廓偏差的一齿切向综合偏差计算方法。通过齿轮测量仪器对齿廓进行测量。并将测量出来的真实齿面信息输入到计算机中,对真实齿面进行建模。之后通过将真实齿廓形状信息经过坐标平移、坐标旋转等处理,计算出齿轮的一齿切向综合偏差。
背景技术
齿轮作为传动机构的一种,与轴、轴承等一样,在人类历史上,有着悠久的历史,广泛应用于机械领域。齿轮具有传动比恒定、结构紧凑等多种优点。其中,一齿切向综合偏差对于齿轮传动稳定性会产生较大影响,会直接影响传动***的振动噪声、工作精度等。一齿切向综合偏差作为振动噪声的重要影响因素,对其进行测量与计算,对于改进齿轮箱性能、降低汽车噪音,提升汽车驾驶舒适度有重要意义。
目前常用的齿轮测量手段有齿轮测量中心、三坐标测量机等,这些测量仪器不能直接实现对齿廓偏差、齿轮粗糙度等项目的测量。同时,这些仪器难以直接输出因齿廓偏差引起的一齿切向综合偏差结果。
目前,一齿切向综合偏差测量主要在单面啮合测量仪上进行,通过测量两个齿轮转动角度来获取一齿切向综合偏差。但是这种方法要求的齿轮装夹精度要求较高,并且由于实际齿轮传动的重合度大于1,无法针对某一个齿廓进行单独分析。齿轮整体误差测量技术能够获取一齿切向综合偏差,但是需要专门的仪器。
发明内容
本发明公开了一种方法,可以实现根据实测齿廓偏差曲线获取一齿切向综合偏差,其基本思想是让一个具有理想齿廓的齿轮与被测齿轮进行啮合,通过计算求得一齿切向综合偏差。在此之后,这些数据也可以通过其他的算法,用于计算该齿轮的滑动速度、赫兹接触应力、油膜厚度等参数,大大提高了数据的利用率。本发明所采用的方法包括以下步骤:
S1基于齿数、模数等信息对理论齿廓进行建模以坐标原点O(0,0)为齿轮A(主动轮)的中心点。以O2(a cosα,a sinα)为齿轮B(被动轮,a为两齿轮的中心距,α为压力角)的中心点,根据不同齿轮的齿数、模数、压力角、变位系数等参数建立理想齿廓模型,通过参数方程建立理论渐开线齿廓模型。于齿轮A基圆与x轴交点T1(rb,0)作x轴垂线,该垂线即为两个齿轮的啮合线。
S2获取真实齿廓信息,对其进行滤波处理,并将齿廓进行拟合通过测量获取被动齿轮齿廓的真实齿面信息。对其进行高斯低通滤波。在步骤1建模结果的基础之上,将被动轮的理论齿廓S2旋转到两个齿轮节点处,求解出S2上各个点与T2点的距离l2,将图2中各个点的偏差值Δl2与l2相加,得到L2。再将L2恢复成齿廓形式,即为拟合完成。将拟合后的齿廓记为S′2
S3基于算法实现主动轮理论齿廓与被动轮实测齿廓的啮合
设S1(主动轮的理论齿廓)与啮合线的交点为P1。在齿轮A转动的过程中,T1P1的长度是不断变化的。任意选取T1P1的值,令S1旋转到该处,在啮合计算开始之前,将步骤2中获得的拟合后的齿廓S′2逆时针旋转较大角度,使得S′2整条齿廓在S1下方。在后续的计算过程中,将S′2不断顺时针旋转微小角度。在每次旋转之后,在两个齿轮啮合区间(xb2′,xb1)上计算出Δyi=y′2i-y1i。当Δyi>0对于每一个xi都成立时,停止运算,此时的r2i即满足条件。即刻终止运算,此时认为两条曲线是相切的,输出S′2的旋转角度θn
S4求解一齿切向综合偏差
通过步骤3的方法,可以求解出S′2的旋转角度θn和两条理论齿廓啮合时,被动齿轮齿廓需要旋转的角度θt。计算出两个角度之间的差值Δθ=θnt,即可求解出一齿切向综合偏差。
本发明提供了一种根据实际测量齿廓计算出一齿切向综合偏差的方法。实现了仅通过测量一个齿的实际齿廓就可以输出一齿切向综合偏差。该方法可以实现对一个齿廓进行单独分析,降低了测量过程中齿轮装夹难度。不需要专门的仪器就可以获取一齿切向综合偏差。
附图说明
图1为两个齿廓啮合示意图。
图2为滤波之后的真实齿面偏差图。
图3为一齿切向综合偏差计算结果图。
图4为本发明的计算流程图。
图中:1、齿顶圆,2、分度圆,3、基圆,4、齿根圆,5、真实齿廓,6、理想齿廓,7、中心线,8、啮合线。
具体实施方式
下面结合附图对本实施例作进一步说明:
1.基于齿数、模数等信息对理论齿廓进行建模
以坐标原点O(0,0)为齿轮A(主动轮)的中心点。以O2(a cosα,a sinα)为齿轮B(被动轮,a为两齿轮的中心距,α为压力角)的中心点,根据不同齿轮的齿数、模数、压力角、变位系数等参数建立理想齿廓模型,通过参数方程建立理论渐开线齿廓模型。于齿轮A基圆与x轴交点T1(rb,0)作x轴垂线,该垂线即为两个齿轮的啮合线。通过此种方式进行建模,两齿轮的啮合线垂直于x轴,两个齿轮会沿着y轴方向进行啮合,为后续计算提供便利。两个齿轮模型的位置与形状、中心点的位置、中心线等如图1所示。
2.获取真实齿廓信息,对其进行滤波处理,并将齿廓进行拟合
i.通过粗糙度轮廓仪对齿轮进行测量,得到被动齿轮齿廓的真实齿面信息。将真实齿廓信息输入到计算机中。对其进行截至波长为0.5的高斯低通滤波,以去除粗糙度信息。
ii.图2代表的是真实齿廓相较于理论齿廓在齿廓各个点的偏差值。在进行真实齿廓与理论齿廓拟合的时候,在步骤1建模结果的基础之上,将被动轮的理论齿廓S2旋转到两个齿轮节点处,求解出S2上各个点与T2点的距离l2,然后只需要将图2中各个点的偏差值Δl2与l2相加,得到L2。再将L2恢复成齿廓形式,即为拟合完成。将拟合后的齿廓记为S′2
3.基于算法实现主动轮理论齿廓与被动轮实测齿廓的啮合设S1(主动轮的理论齿廓)与啮合线的交点为P1。在齿轮A转动的过程中,T1P1的长度是不断变化的。变化区间为(T1A,T1E),其中,T1点代表啮合线与齿轮A基圆的交点,A为两个齿轮啮合的起始点,E为两个齿轮啮合的结束点。任意选取T1P1的值,令S1旋转到该处,以进行后续的啮合计算;在啮合计算开始之前,将步骤2(ii)中获得的拟合后的齿廓S′2逆时针旋转较大角度,使得S′2整条齿廓在S1下方。在后续的计算过程中,将S′2不断顺时针旋转微小角度。在每次旋转之后,在两个齿轮啮合区间(xb2′,xb1)上计算出Δyi=y′2i-y1i。其中,xb2′为S′2齿顶处的横坐标,xb1为S1齿顶处的横坐标。设S1上点的坐标为(xi,y1i),设S′2上每个点的坐标为(xi,y′2i)。在S′2顺时针旋转的过程中,会依次出现三种情况,如式(1)所示。其中Δyi的详细计算方法,如式(2)~(6)所示。
Δx=xb1-xb2′ (4)
Δyi=y′2i-y1i,i=1,2,...,n (6)
其中,a是两个齿轮的中心距。θn为旋转角度,α为压力角。a×cos(α)、a×sin(α)分别代表S′2旋转中心的横、纵坐标。r2i代表将S′2旋转一定角度之后,S′2上各个点的坐标。Δx代表每次旋转之后,S1齿顶点横坐标与S′2上齿顶点的横坐标的差值。xi代表着,在啮合计算区间(xb2′,xb1)中依次取值。将取到的xi分别带入到S1和S′2中求解对应着的y1i和y′2i。之后再求取每一个xi所对应的Δyi。当Δyi>0对于每一个xi都成立时,停止运算,此时的r2i即满足条件。当啮合计算区间上xi取值足够密集,S′2的旋转角度的间隔足够小时,一旦不再出现y′2i-y1i<0的情况,即刻终止运算,此时认为两条曲线是相切的,输出S′2的旋转角度θn
在计算的过程中,如若一次性将旋转角度分得过于精细,无用计算部分的计算量就会大幅度增加,进而极大地增加计算成本和时间成本。针对此种情况,选取较大的旋转角度的间隔,先找到一个目标角度的大致区间。之后进一步将这个区间进行细分。如此便可保证在尽可能减少计算成本和时间成本的前提下,提高计算精度。
4.求解一齿切向综合偏差
在步骤3描述的方法下,可以求解出S′2的旋转角度θn。同时,利用相同的方法,也可以求解出两条理论齿廓啮合时,被动齿轮齿廓需要旋转的角度θt。此时,可以求解出两个角度之间的差值Δθ=θnt,在再行一定的单位转换运算,即可求解出一齿切向综合偏差。其计算结果如图3所示。

Claims (1)

1.一种基于实测齿廓偏差的一齿切向综合偏差计算方法,其特征在于:包括如下步骤,
S1基于齿数、模数信息对理论齿廓进行建模;
以坐标原点O(0,0)为齿轮A即主动轮的中心点;以O2(a cosα,a sinα)为齿轮B即被动轮的中心点,根据不同齿轮的齿数、模数、压力角、变位系数参数建立理论齿廓模型,通过参数方程建立理论渐开线齿廓模型;于齿轮A基圆与x轴交点T1(rb,0)作x轴垂线,即为两个齿轮的啮合线;a是两个齿轮的中心距,α为压力角;
S2获取真实齿廓信息,对其进行滤波处理,并将齿廓进行拟合;
通过测量获取被动齿轮齿廓的真实齿面信息;将真实齿廓信息输入到计算机并进行截至波长为0.5的高斯低通滤波,以去除粗糙度信息;在S1建模结果的基础之上,将被动轮的理论齿廓S2旋转到两个齿轮节点处,求解出S2上各个点与T2点的距离l2,将各个点的偏差值Δl2与l2相加,得到L2;再将L2恢复成齿廓形式,拟合完成;将拟合后的齿廓记为S′2
S3基于算法实现主动轮理论齿廓与被动轮实测齿廓的啮合;
设主动轮的理论齿廓S1与啮合线的交点为P1;在齿轮A转动的过程中,T1P1的长度不断变化;任意选取T1P1的值,令S1旋转到该处,在啮合计算开始之前,将S2中获得的拟合后的齿廓S′2逆时针旋转较大角度,使得S′2整条齿廓在S1下方;将S′2不断顺时针旋转微小角度;在每次旋转之后,在两个齿轮啮合区间(xb2′,xb1)上计算出Δyi=y′2i-y1i;当Δyi>0对于每一个xi都成立时,停止运算,此时的r2i即满足条件;即刻终止运算,此时认为两条曲线是相切的,输出S′2的旋转角度θn;xb2′为S′2齿顶处的横坐标,xb1为S1齿顶处的横坐标;y1i、y′2i分别为S1和S′2上每个点的纵坐标;r2i代表将S′2旋转一定角度之后,S′2上各个点的坐标;
S4求解一齿切向综合偏差;
通过S3求解出S′2的旋转角度θn;两条理论齿廓啮合时,通过与S3相同的方法求出被动齿轮齿廓需要旋转的角度θt;计算出两个角度之间的差值Δθ=θnt,即可求解出一齿切向综合偏差。
CN202010973033.XA 2020-09-16 2020-09-16 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 Active CN112257192B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010973033.XA CN112257192B (zh) 2020-09-16 2020-09-16 一种基于实测齿廓偏差的一齿切向综合偏差计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010973033.XA CN112257192B (zh) 2020-09-16 2020-09-16 一种基于实测齿廓偏差的一齿切向综合偏差计算方法

Publications (2)

Publication Number Publication Date
CN112257192A CN112257192A (zh) 2021-01-22
CN112257192B true CN112257192B (zh) 2024-03-08

Family

ID=74231890

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010973033.XA Active CN112257192B (zh) 2020-09-16 2020-09-16 一种基于实测齿廓偏差的一齿切向综合偏差计算方法

Country Status (1)

Country Link
CN (1) CN112257192B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113063589B (zh) * 2021-03-18 2023-07-04 重庆青山工业有限责任公司 一种基于神经网络的齿轮微观误差对振动的预测方法
CN113283029B (zh) * 2021-05-21 2023-11-21 北京工业大学 基于迭代滤波的实际齿廓上包络形状获取方法
CN113591238B (zh) * 2021-07-19 2024-03-22 北京工业大学 基于齿轮点云数据的齿轮切向综合偏差计算方法
CN113819874B (zh) * 2021-09-01 2024-04-02 宜昌长机科技有限责任公司 一种检验变速比及变速比鼓形齿齿扇加工精度的方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107588742A (zh) * 2017-10-25 2018-01-16 北京工业大学 一种基于线结构光的圆柱齿轮齿廓偏差测量方法
CN109655020A (zh) * 2019-01-14 2019-04-19 北方民族大学 齿轮一齿径向综合偏差及切向综合偏差检测方法
CN110455244A (zh) * 2019-08-28 2019-11-15 北京工业大学 一种利用粗糙度轮廓仪进行圆柱齿轮齿距偏差测量的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107588742A (zh) * 2017-10-25 2018-01-16 北京工业大学 一种基于线结构光的圆柱齿轮齿廓偏差测量方法
CN109655020A (zh) * 2019-01-14 2019-04-19 北方民族大学 齿轮一齿径向综合偏差及切向综合偏差检测方法
CN110455244A (zh) * 2019-08-28 2019-11-15 北京工业大学 一种利用粗糙度轮廓仪进行圆柱齿轮齿距偏差测量的方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
渐开线圆柱齿轮精度设计的研究;李海燕;科技视界;全文 *
齿轮单面啮合测量技术的发展及其应用;石照耀;工具技术;全文 *

Also Published As

Publication number Publication date
CN112257192A (zh) 2021-01-22

Similar Documents

Publication Publication Date Title
CN112257192B (zh) 一种基于实测齿廓偏差的一齿切向综合偏差计算方法
Suh et al. Geometric error measurement of spiral bevel gears using a virtual gear model for STEP-NC
WO2006109838A1 (ja) コルヌ螺旋歯形歯車
CN110008610B (zh) 摆线轮齿廓分段优化设计方法
CN106989157B (zh) 基于运动学原理的渐开线圆柱齿轮双面啮合测量仿真方法
CN111666643B (zh) 一种用于确定复杂齿面接触性能的方法
Huang et al. Influence of system-inherent phase in spur gear multi-state mesh based on nonlinear dynamics
CN113221289A (zh) 一种基于区间模型的谐波减速器动态传动误差近似方法
CN101339009A (zh) 齿轮圆柱或球测量方法和***
CN110940283B (zh) 基于白光干涉仪的微小齿轮齿距偏差和齿廓偏差高精度测量方法
CN113446960B (zh) 一种齿面点云理论分布建模方法及测量方法
CN109341629B (zh) 滚刀安装轴交角误差对加工齿轮表面误差影响的分析方法
CN115143922B (zh) 一种面向传动性能分析的齿轮齿面粗糙度测量与评定方法
CN115221655B (zh) 一种斜齿非圆齿轮的传动精度评价方法
CN113566772B (zh) 一种基于坐标测量的局部齿面定位方法
CN113408082B (zh) 一种外啮合直齿圆柱齿轮副动态啮合力的计算方法
CN109190324A (zh) 一种弧齿锥齿轮含误差的齿面载荷接触变形数值计算方法
Niu et al. Geometrical design of variable ratio tooth profile based on Boolean subtraction operation and a novel modification method
KR20030050469A (ko) 스파이럴 베벨기어의 형상오차 측정방법
Tsuji et al. Machining and running test of high-performance face gear set
Shih et al. Precision evaluation for cycloidal gears
CN113591238B (zh) 基于齿轮点云数据的齿轮切向综合偏差计算方法
Suzuki et al. Design of high-reduction hypoid gears meshing in plane of action
CN115098988B (zh) 一种基于传递误差的弧齿锥齿轮内源激励识别方法
Ghionea et al. Analysis on the measurement results of the precision grades and their influence on the performance of the hydraulic pumps with spur gears

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