CN115096239A - 一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法 - Google Patents
一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法 Download PDFInfo
- Publication number
- CN115096239A CN115096239A CN202210660207.6A CN202210660207A CN115096239A CN 115096239 A CN115096239 A CN 115096239A CN 202210660207 A CN202210660207 A CN 202210660207A CN 115096239 A CN115096239 A CN 115096239A
- Authority
- CN
- China
- Prior art keywords
- tooth
- gear
- waviness
- tooth profile
- information
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B21/00—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant
- G01B21/20—Measuring arrangements or details thereof, where the measuring technique is not covered by the other groups of this subclass, unspecified or not relevant for measuring contours or curvatures, e.g. determining profile
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Length Measuring Devices With Unspecified Measuring Means (AREA)
Abstract
本发明公开了一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法,该方法能够基于齿轮测量中心的测量结果,高效、可靠地获取单个轮齿齿面齿廓波纹度信息和连续齿面齿廓波纹度信息,并得到其阶次信息。具体包括以下方法步骤:测量齿轮的所有轮齿齿面齿廓信息,在各个轮齿齿面相同位置选择三条齿廓线进行测量,测量设备为齿轮测量中心。本发明的技术优势是除获取单齿齿面齿廓波纹度信息之外,还可以将齿轮中轮齿波度信息连续起来,从而获得整个齿轮的连续波度信息。
Description
技术领域
本发明涉及一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法,属于齿轮传动技术领域。
背景技术
齿轮作为传动的关键零部件,因其具有传动平稳、效率高、使用寿命长等特点,被广泛应用于高铁、汽车、精密机床等领域。人们对装备中齿轮传动***使用要求的提高,使得齿轮向着高精度、低振动、低噪声和高可靠性的方向发展。
齿轮传动的振动和噪声是反映齿轮传动性能的重要指标,齿轮的振动和噪声水平影响着用户的使用感受,而齿轮表面形貌影响着齿轮齿轮传动的振动和噪声水平。目前,大量的研究聚焦于齿轮齿面粗糙度和齿形误差对齿轮传动噪声的影响,并通过提高齿轮表面质量和利用齿面修形的手段改善齿轮的传动噪声。专门针对介于宏观齿形偏差和微观齿面粗糙度之间的齿面波纹度,研究其与齿轮传动振动和噪声的关联性的工作较少。因此针对齿轮齿面波纹度开展研究工作具有重要的理论和实际价值,里面涉及的关键问题之一就是齿面波纹度的有效提取与分析。
因此,本发明基于上述关键技术问题:提出一种用于获取渐开线直齿圆柱齿轮齿廓波纹度的方法,并对齿廓波纹度信息进行分析,该方法能够基于齿轮测量中心的测量结果,高效、可靠地获取齿轮齿廓波纹度信息,并得到其阶次信息。
发明内容
本发明采用的技术方案为一种用于获取及分析渐开线直齿圆柱齿轮齿廓波纹度的方法,该方法能够基于齿轮测量中心的测量结果,高效、可靠地获取单个轮齿齿面齿廓波纹度信息和连续齿面齿廓波纹度信息,并得到其阶次信息。具体包括以下方法步骤:
测量齿轮的所有轮齿齿面齿廓信息,在各个轮齿齿面相同位置选择三条齿廓线进行测量,测量设备为齿轮测量中心。
(1)单个轮齿齿面齿廓波纹度的获取
采用测得的三条齿廓线的信息开展分析,测量结果中包含表面粗糙度、波纹度和形状误差三部分,通过对齿廓信息进行滤波处理,即可获取齿廓波纹度信息。
齿轮测量中心的直接测量结果是以测量时齿轮的转动角度为横坐标,即展角值,对测量结果进行坐标变换,转换公式为式1。
l:渐开线长度值Rb:基圆半径θ:齿轮转角
通过坐标变换,将齿廓的测量结果转变为由渐开线长度值为横坐标。
接下来提取单齿齿面波纹度,使用滤波器法需要满足截止波长大于采样长度,情况一:单齿齿面实际采集齿廓渐开线展开长度大于采样长度时,直接使用滤波器法滤除齿廓信息中的形状误差。情况二:单齿齿面实际采集齿廓渐开线展开长度小于采样长度时,可采用多项式拟合法滤除形状误差。
情况一下直接使用高斯滤波器滤除形状误差。
情况二下由于齿廓形状误差通常为简单的曲线,采用最小二乘法非线性拟合去除齿廓形状误差。
高斯滤波是目前最常用的过滤粗糙度的方法,且发展较为成熟,因此接下来使用高斯滤波器对表面粗糙度进行滤波处理。表面粗糙度在采集的信号中为高频成分,采用高斯滤波器进行低通滤波,即可去除采集信号中的表面粗糙度成分。
(2)齿轮齿廓波纹度连续
齿轮传动是一个连续的过程,产生的齿轮传动噪声是连续信号,而单齿齿面信息为独立的数据,因此接下来需要对齿轮各个轮齿齿面信息进行连续。
依据展成原理,将整个齿轮的单齿齿面对应位置的齿廓信息依次连续,对连续过程中,齿面信息横坐标重合的部分,选取幅值高的值作为提取的信息,即可以得到齿轮360度连续的完整齿面信息。
对连续的齿面信息,可以直接进行滤波处理,利用高斯滤波器去除形状误差和表面粗糙度,即可得到整个齿轮连续齿面齿廓的波纹度。
(3)阶次分析
对滤波后的数据进行阶次分析,阶次与频率转化关系如式2所示。
其中,频率单位为Hz,转速单位为r/min。
对得到的单齿齿面齿廓波纹度和连续齿面齿廓波纹度信息依次进行傅里叶变换即可得到对应的齿廓波纹度的阶次信息,阶次信息高幅值前三阶进行统计分析。
与现有技术相比,本发明的技术优势是除获取单齿齿面齿廓波纹度信息之外,还可以将齿轮中轮齿波度信息连续起来,从而获得整个齿轮的连续波度信息。
附图说明
下面基于附图更详细的描述本发明,附图示出了示例性实施例。在示例附图中:
图1显示了单轮齿齿面齿廓测量位置;
图2显示齿面齿廓信息连续原理;a)齿廓信息连续之前;b)齿廓信息连续之后;
图3显示齿轮齿廓波纹度提取过程;a)滤除形状误差,b)波纹度与表面粗糙度,c)滤除表面粗糙度,d)波纹度。
图4显示齿轮齿廓波纹度阶次分析结果;a)单齿齿廓波纹度阶次信息;b)连续齿廓波纹度阶次信息。
具体实施方式
测量齿轮的单齿齿面齿廓信息,选用克林贝格P26齿轮测量中心,选用合适的红宝石测头,按如图1所示测量路径进行测量(3条)。
测量后得到三组测量数据,根据式(1)进行坐标变换,将横坐标值由转动角度变换为渐开线长度值。
情况一下直接使用高斯滤波器滤除形状误差。
情况二下通过最小二乘法进行多项式拟合,采用拟合的多项式进行滤波,去除采集信号中的形状误差。
具体为:原始数据有n个数据点,对应的坐标为Pi(xi,yi),i=1,2,…,n。设定目标函数型为低次多项式,如式(3):
f(x)=ax2+bx+c#(3)
式中a、b、c为常数。
采用最小二乘法构造评价函数,如式(4):
求得使得评价函数值最小的a,b,c,即可得到拟合的多项式。采用拟合的多项式进行滤波,即可去除采集信号中的形状误差。
对去除形状误差的测量数据,进行高斯滤波,去除高频的粗糙度成分。
具体为:采集到的齿面轮廓信息函数为P(x),其中高频信号为h(x),低频信号为l(x),则:
P(x)=h(x)+l(x)#(5)
采用高斯滤波来去除粗糙度,高斯分布的表达式为:
其中μ是均值,σ是标准差。
高斯滤波的权函数为:
将(5)进行傅里叶变换得:
式中λ是波长,λc是高斯滤波器的截止波长,α是常数。
通过齿面轮廓信息函数和高斯权函数做卷积可以得到粗糙度为:
通过低通滤波过滤掉表面粗糙度信息,即可得到齿面的波纹度信息。
依次测量各个齿面信息,采用展成原理,将单齿齿面对应位置齿廓信息依次连续,得到连续的完整齿面信息,连续原理示意如图2所示,将图2a)横坐标重叠部分保留纵坐标较大值,即可得到图2b)。对连续后的信息直接进行滤波处理,利用高斯滤波器去除形状误差和表面粗糙度,即可得到整个齿轮齿面齿廓的波纹度。对实际测量获得信息进行齿廓波纹度的提取与分析,过程如图3所示。
对齿面齿廓的波纹度信息进行傅里叶变换,即可得到波纹度阶次信息,以齿面一个位置的信息为例,分析结果如图4所示。
Claims (1)
1.一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法,其特征在于:具体包括以下方法步骤,
(1)单个轮齿齿面齿廓波纹度的获取;
采用测得的三条齿廓线的信息开展分析,测量结果中包含表面粗糙度、波纹度和形状误差三部分,通过对齿廓信息进行滤波处理,即可获取齿廓波纹度信息;
齿轮测量中心的直接测量结果是以测量时齿轮的转动角度为横坐标,即展角值,对测量结果进行坐标变换,转换公式为式1;
l:渐开线长度值Rb:基圆半径θ:齿轮转角
通过坐标变换,将齿廓的测量结果转变为由渐开线长度值为横坐标;
提取单齿齿面波纹度,使用滤波器法需要满足截止波长大于采样长度,情况一:单齿齿面实际采集齿廓渐开线展开长度大于采样长度时,直接使用滤波器法滤除齿廓信息中的形状误差;情况二:单齿齿面实际采集齿廓渐开线展开长度小于采样长度时,采用多项式拟合法滤除形状误差;
情况一下直接使用高斯滤波器滤除形状误差;
情况二下由于齿廓形状误差通常为简单的曲线,采用最小二乘法非线性拟合去除齿廓形状误差;
使用高斯滤波器对表面粗糙度进行滤波处理;表面粗糙度在采集的信号中为高频成分,采用高斯滤波器进行低通滤波,即可去除采集信号中的表面粗糙度成分;
(2)齿轮齿廓波纹度连续;
依据展成原理,将整个齿轮的单齿齿面对应位置的齿廓信息依次连续,对连续过程中,齿面信息横坐标重合的部分,选取幅值高的值作为提取的信息,即得到齿轮360度连续的完整齿面信息;
对连续的齿面信息,直接进行滤波处理,利用高斯滤波器去除形状误差和表面粗糙度,即可得到整个齿轮连续齿面齿廓的波纹度;
(3)阶次分析;
对滤波后的数据进行阶次分析,阶次与频率转化关系如式2所示;
其中,频率单位为Hz,转速单位为r/min;
对得到的单齿齿面齿廓波纹度和连续齿面齿廓波纹度信息依次进行傅里叶变换即可得到对应的齿廓波纹度的阶次信息,阶次信息高幅值前三阶进行统计分析。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210660207.6A CN115096239B (zh) | 2022-06-13 | 2022-06-13 | 一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210660207.6A CN115096239B (zh) | 2022-06-13 | 2022-06-13 | 一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115096239A true CN115096239A (zh) | 2022-09-23 |
CN115096239B CN115096239B (zh) | 2023-08-18 |
Family
ID=83291199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210660207.6A Active CN115096239B (zh) | 2022-06-13 | 2022-06-13 | 一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115096239B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0663819A (ja) * | 1992-08-20 | 1994-03-08 | Daihatsu Motor Co Ltd | 歯車検査装置 |
JPH0911085A (ja) * | 1995-06-27 | 1997-01-14 | Kobe Steel Ltd | 歯車誤差計測方法及び補正方法並びに計測/補正装置 |
JP2004098191A (ja) * | 2002-09-06 | 2004-04-02 | Ntn Corp | ボールねじ及びそのねじ溝研削方法 |
CN102706306A (zh) * | 2012-03-01 | 2012-10-03 | 北京工业大学 | 渐开线齿轮波度样板 |
CN203964850U (zh) * | 2014-07-02 | 2014-11-26 | 成都飞机工业(集团)有限责任公司 | 检测平面阶差与波纹度的测量工具 |
JP2018185166A (ja) * | 2017-04-24 | 2018-11-22 | 本田技研工業株式会社 | 歯面うねり評価方法 |
US20190368863A1 (en) * | 2018-05-29 | 2019-12-05 | Klingelnberg Gmbh | Method for analyzing surface waviness |
CN113434817A (zh) * | 2021-05-06 | 2021-09-24 | 北京工业大学 | 一种齿轮单项拓扑误差图谱的分析方法 |
CN114216677A (zh) * | 2021-12-15 | 2022-03-22 | 重庆青山工业有限责任公司 | 齿面波纹阶次检测方法 |
-
2022
- 2022-06-13 CN CN202210660207.6A patent/CN115096239B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0663819A (ja) * | 1992-08-20 | 1994-03-08 | Daihatsu Motor Co Ltd | 歯車検査装置 |
JPH0911085A (ja) * | 1995-06-27 | 1997-01-14 | Kobe Steel Ltd | 歯車誤差計測方法及び補正方法並びに計測/補正装置 |
JP2004098191A (ja) * | 2002-09-06 | 2004-04-02 | Ntn Corp | ボールねじ及びそのねじ溝研削方法 |
CN102706306A (zh) * | 2012-03-01 | 2012-10-03 | 北京工业大学 | 渐开线齿轮波度样板 |
CN203964850U (zh) * | 2014-07-02 | 2014-11-26 | 成都飞机工业(集团)有限责任公司 | 检测平面阶差与波纹度的测量工具 |
JP2018185166A (ja) * | 2017-04-24 | 2018-11-22 | 本田技研工業株式会社 | 歯面うねり評価方法 |
US20190368863A1 (en) * | 2018-05-29 | 2019-12-05 | Klingelnberg Gmbh | Method for analyzing surface waviness |
CN110542389A (zh) * | 2018-05-29 | 2019-12-06 | 科令志因伯格有限公司 | 用于分析表面波纹度的方法 |
CN113434817A (zh) * | 2021-05-06 | 2021-09-24 | 北京工业大学 | 一种齿轮单项拓扑误差图谱的分析方法 |
CN114216677A (zh) * | 2021-12-15 | 2022-03-22 | 重庆青山工业有限责任公司 | 齿面波纹阶次检测方法 |
Non-Patent Citations (2)
Title |
---|
T REAVIE等: "A method to quantify waviness measurement capability of a gear measurement machine" * |
王泽贵: "齿轮表面波纹度阶次噪声机理分析及改进", vol. 34, no. 34 * |
Also Published As
Publication number | Publication date |
---|---|
CN115096239B (zh) | 2023-08-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Meltzer et al. | Fault detection in gear drives with non-stationary rotational speed-part I: the time-frequency approach | |
US4831872A (en) | Method of measuring gear errors through meshing and separating and a machine thereof | |
CN105300691B (zh) | 基于最佳安装距的锥齿轮传动误差测量方法 | |
CN108844733B (zh) | 一种基于kl散度与均方根值的齿轮状态监测指标提取方法 | |
US4704799A (en) | Master gear error compensation | |
CN108007326A (zh) | 一种采用齿轮测量中心测量评定齿面粗糙度的方法 | |
CN106153339B (zh) | 一种基于变分模式分解滤波的包络分析方法 | |
CN112257192B (zh) | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 | |
CN111238813A (zh) | 一种强干扰下提取滚动轴承故障特征的方法 | |
Lin et al. | Determination of proper frequency bandwidth for 3D topography measurement using spectral analysis. Part I: isotropic surfaces | |
Li et al. | Research on a signal separation method based on Vold-Kalman filter of improved adaptive instantaneous frequency estimation | |
Przyklenk et al. | Holistic evaluation of involute gears | |
DE69712749D1 (de) | Verfahren zur berechnung eines werkstücks zur bearbeitung | |
CN115096239A (zh) | 一种渐开线直齿圆柱齿轮齿廓波纹度的获取及分析方法 | |
CN114714157A (zh) | 一种基于时变滤波经验模态分解与瞬时能量比的磨削颤振监测方法 | |
CN111323233B (zh) | 一种用于低速旋转机械故障诊断的局部均值分解方法 | |
CN106153333B (zh) | 一种基于小波分解滤波的包络分析方法 | |
CN109099876B (zh) | 蜗轮的检查方法 | |
CN113705421B (zh) | 一种磨削工件表面振纹在线监测方法及*** | |
Srinivasan et al. | A form tolerancing theory using fractals and wavelets | |
CN116203270A (zh) | 一种利用非均匀反光斑马带的转子轴系转速测量方法 | |
CN113446960B (zh) | 一种齿面点云理论分布建模方法及测量方法 | |
Stanciu et al. | MASTER WHEEL ERRORS ELIMINATION IN FORCED GEAR ENGAGEMENT TESTING MACHINE USING HARMONIC ANALYSIS THROUGH FAST FOURIER TRANSFORM | |
CN106053061B (zh) | 一种基于非线性模式分解和谱峭度的包络分析方法 | |
CN110595751B (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 |