CN113283029B - 基于迭代滤波的实际齿廓上包络形状获取方法 - Google Patents
基于迭代滤波的实际齿廓上包络形状获取方法 Download PDFInfo
- Publication number
- CN113283029B CN113283029B CN202110556915.0A CN202110556915A CN113283029B CN 113283029 B CN113283029 B CN 113283029B CN 202110556915 A CN202110556915 A CN 202110556915A CN 113283029 B CN113283029 B CN 113283029B
- Authority
- CN
- China
- Prior art keywords
- filtering
- curve
- tooth profile
- filter
- filtered
- 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
Links
- 238000001914 filtration Methods 0.000 title claims abstract description 70
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000000877 morphologic effect Effects 0.000 claims description 3
- 238000005259 measurement Methods 0.000 abstract description 12
- 230000000694 effects Effects 0.000 abstract description 5
- 238000012545 processing Methods 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 3
- 238000005299 abrasion Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000004026 adhesive bonding Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000005461 lubrication Methods 0.000 description 1
- 230000002459 sustained effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
- G01B11/2416—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures of gears
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (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)
Abstract
本发明公开了基于迭代滤波的实际齿廓上包络形状获取方法,所采用的方法包括以下步骤:通过迭代滤波的方法获取实际齿廓上包络形状曲线;停止条件的设定和截止波长的选择;截止波长λC决定滤波的效果。截止波长选择较小会导致无法消除微小裂纹对滤波结果的影响,截止波长较大会使滤波结果与实测齿廓曲线脱离。所以需要根据测量曲线的实际情况选择合适的截止波长,以获得最佳滤波效果。现有的滤波手段无法获取齿廓测量曲线的上包络曲线,且针对于齿廓测量曲线无法通过传统信号处理技术中的方法获取上包络线。本方法针对于实际齿廓测量曲线提出通过迭代滤波的方法实现了实际齿廓上包络形状的获取。
Description
技术领域
本发明属于齿轮精密测量技术领域,具体为一种基于齿廓偏差的上包络形状获取方法。通过粗糙度轮廓仪对齿轮的齿廓进行测量,经过计算获取齿廓偏差,齿廓偏差通过迭代滤波的方式,获得实测齿廓偏差的上包络形状,可用于齿轮试验中磨损量的计算。
背景技术
齿轮作为传动机构的一种,有着悠久的历史,广泛应用于机械领域。齿轮具有传动比恒定、结构紧凑等多种优点。据统计,大量的机械零部件因磨损而失效。齿轮在服役过程中,由于设计不当、润滑失效、承受交变载荷及材料缺陷等原因会发生微点蚀、点蚀、胶合甚至断齿等失效。如何延缓齿轮失效一直是齿轮领域的研究热点问题之一。齿轮失效常表现为齿面磨损,即齿面材料的持续损伤和脱落,齿廓形状发生变化,这一过程贯穿齿轮的全服役周期。
齿轮磨损的精确测量是进行磨损机理研究的基础。目前常用的测量仪器有齿轮测量中心,三坐标测量机和粗糙度轮廓仪,但最适用于齿面磨损测量的是粗糙度轮廓仪。利用粗糙度轮廓仪测量齿轮的优点是可同时获得齿面的形状、波纹度和粗糙度,并可通过曲线拟合的方式确定被测齿廓在实际齿面的位置,使得利用粗糙度轮廓仪进行齿面磨损测量成为可能。
发明内容
本发明公开一种迭代滤波方法,对利用粗糙度轮廓仪获取的齿廓偏差曲线进行处理,逐渐逼近实测曲线的上表面包络轮廓。本发明所采用的方法包括以下步骤:
S1通过迭代滤波的方法获取实际齿廓上包络形状曲线;
滤波器选择由ISO 1328:2013齿轮精度国际标准推荐采用的高斯滤波器进行滤波演示,也可根据需要选择样条滤波器、形态学滤波器、巴特沃斯滤波器、切比雪夫滤波器或其他类型的数字低通滤波器。
设待滤波齿廓形状曲线为Fi(x,y),第一次滤波时将原始测量齿廓形状曲线作为待滤波齿廓形状曲线;经过滤波之后的测量齿廓形状曲线为Fi(x',y'),i=0,1,2...n,n为滤波迭代的次数。以每次滤波后的结果为比较基准,保留位于该次滤波结果曲线上方的原始曲线不变,将位于该次滤波结果曲线下方的点用滤波曲线代替形成新的待滤波齿廓形状曲线,作为下次滤波的待滤波曲线。其计算过程如式(1)所示。
Fi+1(xk,yk)=max{Fi(xk,yk),Fi(xk',yk')},k=1,2...K (1)
式中,K为齿廓曲线数据点个数。Fi(xk,yk)为第i次滤波时待滤波曲线上第k点的坐标,Fi(xk',yk')为第i次滤波时滤波结果曲线上第k点的坐标。Fi+1(xk,yk)为第i+1次滤波时待滤波结果曲线上第k点的坐标。
在每次滤波结束之后,判断是否满足停止条件。若满足停止条件,则迭代滤波结束。
S2停止条件的设定和截止波长的选择;
设迭代停止参数为mi,n为迭代次数,K为数据点总个数,且需满足Fi(xk',yk')>Fi(xk,yk)。设定mi<m1×T为停止条件,T根据实际情况选取。当满足停止条件时,结束运算,输出结果。
截止波长λC决定滤波的效果。截止波长选择较小会导致无法消除微小裂纹对滤波结果的影响,截止波长较大会使滤波结果与实测齿廓曲线脱离。所以需要根据测量曲线的实际情况选择合适的截止波长,以获得最佳滤波效果。
现有的滤波手段无法获取齿廓测量曲线的上包络曲线,且针对于齿廓测量曲线无法通过传统信号处理技术中的方法获取上包络线。本方法针对于实际齿廓测量曲线提出通过迭代滤波的方法实现了实际齿廓上包络形状的获取。
附图说明
图1为实验前齿廓偏差图。
图2为实验后齿廓偏差图。
图3为迭代滤波结果图。
图4为迭代滤波流程图。
具体实施方式
下面结合附图对本实施例作进一步说明:
通过迭代滤波的方法获取实际齿廓上包络形状
齿轮试验过程中,通过粗糙度轮廓仪测得的齿廓数据包含粗糙度、波纹度和齿廓形貌等信息。随着试验的进行,齿轮表面逐渐磨损,表现为被测齿廓偏差曲线的某些部位变化较为剧烈,如出现微裂纹等情况(如图1、图2所示),利用本发明的方法可获得齿廓偏差曲线的上包络轮廓,为进行齿面磨损分析奠定基础。
高斯滤波器是ISO 1328:2013齿轮精度国际标准推荐采用的滤波器,本发明以高斯滤波器为例进行迭代滤波演示,也可根据需要选择样条滤波器、形态学滤波器或其他类型的数字低通滤波器。高斯滤波器权函数定义如式所示;其频域相应函数如式所示。
式中:λ为波长;λc为高斯滤波器的截止波长,为了保证滤波器的通过率为50%,取
设待滤波齿廓形状曲线为Fi(x,y),第一次滤波时将原始测量齿廓形状曲线作为待滤波齿廓形状曲线;经过滤波之后的测量齿廓形状曲线为Fi(x',y'),i=0,1,2...n,n为滤波迭代的次数。n为一个不确定的数,随停止条件的变化而变化。在每次滤波结束之后,判断是否满足停止条件。若满足停止条件,则停止运算。当不满足停止条件时,以每次滤波后的结果为比较基准,保留位于该次滤波结果曲线上方的原始曲线不变,将位于该次滤波结果曲线下方的点用滤波曲线代替形成新的待滤波齿廓形状曲线,作为下次滤波的待滤波曲线。其计算过程如式所示。
Fi+1(xk,yk)=max{Fi(xk,yk),Fi(xk',yk')},k=1,2...K (3)
式中,K为齿廓曲线数据点个数。Fi(xk,yk)为第i次滤波时待滤波曲线上第k点的坐标,Fi(xk',yk')为第i次滤波时滤波结果曲线上第k点的坐标。Fi+1(xk,yk)为第i+1次滤波时待滤波结果曲线上第k点的坐标。
停止条件的设定和截止波长的选择
设迭代停止参数为mi,n为迭代次数,K为数据点总个数,且需满足Fi(xk',yk')>Fi(xk,yk)。设定mi<m1×T为停止条件,针对本发明案例,选择T=20%。当满足停止条件时,结束运算,输出结果。
截止波长λC决定了滤波的效果。通过分析齿轮偏差曲线的特征,本算例选择的截止波长位于0.05mm-0.15mm,对比后最终选择0.13mm作为本组曲线滤波迭代的截止波长。滤波迭代的结果如图3所示。
Claims (2)
1.一种基于迭代滤波的实际齿廓上包络形状获取方法,其特征在于:包括以下步骤,
S1通过迭代滤波的方法获取实际齿廓上包络形状曲线;
滤波器采用的高斯滤波器进行滤波演示,设待滤波齿廓形状曲线为Fi(x,y),第一次滤波时将原始测量齿廓形状曲线作为待滤波齿廓形状曲线;经过滤波之后的测量齿廓形状曲线为Fi(x',y'),i=0,1,2...n,n为滤波迭代的次数;以每次滤波后的结果为比较基准,保留位于该次滤波结果曲线上方的原始曲线不变,将位于该次滤波结果曲线下方的点用滤波结果曲线代替形成新的待滤波齿廓形状曲线,作为下次滤波的待滤波曲线;其计算过程如下式所示;
Fi+1(xk,yk)=max{Fi(xk,yk),Fi(xk',yk')},k=1,2...K
式中,K为齿廓曲线数据点个数;Fi(xk,yk)为第i次滤波时待滤波曲线上第k点的坐标,Fi(xk',yk')为第i次滤波时滤波结果曲线上第k点的坐标;Fi+1(xk,yk)为第i+1次滤波时待滤波结果曲线上第k点的坐标;
在每次滤波结束之后,判断是否满足停止条件;若满足停止条件,则迭代滤波结束;
S2停止条件的设定和截止波长的选择;
设迭代停止参数为mi,n为迭代次数,K为数据点总个数,且需满足Fi(xk',yk')>Fi(xk,yk);设定mi<m1×T为停止条件,T根据实际情况选取;当满足停止条件时,结束运算,输出结果;选择0.13mm作为曲线滤波迭代的截止波长。
2.根据权利要求1所述的基于迭代滤波的实际齿廓上包络形状获取方法,其特征在于:S1中的滤波器为样条滤波器、形态学滤波器、巴特沃斯滤波器、切比雪夫滤波器或数字低通滤波器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110556915.0A CN113283029B (zh) | 2021-05-21 | 2021-05-21 | 基于迭代滤波的实际齿廓上包络形状获取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110556915.0A CN113283029B (zh) | 2021-05-21 | 2021-05-21 | 基于迭代滤波的实际齿廓上包络形状获取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113283029A CN113283029A (zh) | 2021-08-20 |
CN113283029B true CN113283029B (zh) | 2023-11-21 |
Family
ID=77280785
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110556915.0A Active CN113283029B (zh) | 2021-05-21 | 2021-05-21 | 基于迭代滤波的实际齿廓上包络形状获取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113283029B (zh) |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102966721A (zh) * | 2012-11-29 | 2013-03-13 | 四川大学 | 高可靠精密滤波齿轮 |
CN106777407A (zh) * | 2015-11-19 | 2017-05-31 | 黑龙江恒能自控科技有限公司 | 一种滤波减速器的啮合冲击分析与修形方法 |
CN112257192A (zh) * | 2020-09-16 | 2021-01-22 | 北京工业大学 | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 |
-
2021
- 2021-05-21 CN CN202110556915.0A patent/CN113283029B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102966721A (zh) * | 2012-11-29 | 2013-03-13 | 四川大学 | 高可靠精密滤波齿轮 |
CN106777407A (zh) * | 2015-11-19 | 2017-05-31 | 黑龙江恒能自控科技有限公司 | 一种滤波减速器的啮合冲击分析与修形方法 |
CN112257192A (zh) * | 2020-09-16 | 2021-01-22 | 北京工业大学 | 一种基于实测齿廓偏差的一齿切向综合偏差计算方法 |
Non-Patent Citations (2)
Title |
---|
变速箱齿轮齿廓修形及其噪音振动特性研究;陈磊;知网;全文 * |
齿轮整体误差测量的基础理论及其应用研究;王笑一;知网;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN113283029A (zh) | 2021-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115388817B (zh) | 基于图像处理分析实现铸造件打磨质量检测的方法 | |
CN104296680B (zh) | 一种颗粒增强钛基复合材料磨削表面质量评价方法 | |
CN112903288A (zh) | 一种齿轮三维误差的特征线统一表征方法 | |
CN113283029B (zh) | 基于迭代滤波的实际齿廓上包络形状获取方法 | |
CN108168887B (zh) | 一种考察磨削工艺参数对齿轮啮合疲劳寿命影响的试验方法 | |
CN116713892B (zh) | 用于晶圆薄膜磨削的终点检测方法及设备 | |
Wang et al. | Optimization measurement for the ballscrew raceway profile based on optical measuring system | |
CN113749363A (zh) | 一种珠宝玉石及贵金属简化、系列标准化作业制备方法 | |
CN105043828A (zh) | 一种pcb金相切片精确磨抛到位的方法 | |
CN109099876B (zh) | 蜗轮的检查方法 | |
CN111633559B (zh) | 一种基于表面微观三维形貌的最小加工余量预测方法 | |
CN114119501B (zh) | 一种微铣削加工未变形切削厚度测量方法及*** | |
CN1453559A (zh) | 曲轴弯曲变形的测量方法 | |
Mutambi et al. | Application of digital image analysis method in metric screw thread metrology | |
CN106779322B (zh) | 服从二项分布的计件值工序能力指数评价方法 | |
CN110021027B (zh) | 一种基于双目视觉的切边点计算方法 | |
CN113740058A (zh) | 一种齿轮双面啮合检测方法 | |
CN114964042B (zh) | 一种曲线轮廓在线测量中数据异常点判别与辨识方法 | |
CN117741090B (zh) | 一种变速箱精密铸件质量缺陷检测方法 | |
CN113280708B (zh) | 一种面向复杂型面模具镶块加工精度的快速检测方法 | |
CN110068268A (zh) | 基于变焦测量的圆弧刃车刀几何参数三维检测***及方法 | |
CN117592313B (zh) | 一种形位误差测量不确定度仿真优化方法 | |
CN116956500A (zh) | 基于全坐标检测的齿轮误差评价方法 | |
CN110119567B (zh) | 一种齿轮齿面质量的热力图表示方法 | |
Kubátová et al. | Selecting software filter settings for evaluating periodic and aperiodic surface roughness |
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 |